Screening and Identification of Differentially Expressed Genes Expressed among Left and Right Colon Adenocarcinoma

A Comprehensive Pan-Cancer Analysis Reveals Cyclin-Dependent Kinase Inhibitor 2A Gene as a Potential Diagnostic and Prognostic Biomarker in Colon Adenocarcinoma

Introduction Cyclin-dependent kinase inhibitor 2A (CDKN2A) is a suppressor carcinogenic gene that is upregulated across various types of cancer including breast, liver, thyroid, and bile duct cancer due to its crucial role in cell cycle regulation and cell division. Nevertheless, it is mostly investigated at the genetic level, but it is still poorly studied on pan-cancer analysis as a biomarker and this study shows its significant potential diagnostic and prognostic characteristics. However, this study aims to investigate the role of CDKN2A as a diagnostic and prognostic biomarker across various types of cancer focusing primarily on colon adenocarcinoma (COAD). Methods We investigated CDKN2A gene expression in a pan-cancer analysis across different types of cancer to show its diagnostic potential characteristics by using various bioinformatic tools, including Tumor Immune Estimation Resource (TIMER) 2.0, Gene Expression Profiling Interactive Analysis (GEPIA), and University of Alabama at Birmingham Cancer Data Analysis Portal (UALCAN) database. TIMER was used to profile gene expression across 32 types of cancer composed of 10,000 RNA-seq samples obtained from the Cancer Genome Atlas (TCGA) and to analyze the tumor-infiltrating immune cells. In addition, GEPIA and UALCAN were further used to analyze gene expression, in terms of gene regulation, pathological stages, and clinical parameters, including gender, age, and race. Therefore, we used GEPIA, UALCAN, and Kaplan-Meier plotter particularly across adenocarcinoma to investigate CDKN2A prognosis by studying its high expression association with the patient's overall survival rate to show the tumor progression. Then, we looked into the genetic alteration of CDKN2A by using the cBio Cancer Genomics Portal (cBioPortal), including 10 pan-cancer studies. We concluded the analysis with gene validation by using a public cohort in Gene Expression Omnibus (GEO). Results CDKN2A showed a trend of upregulation in most cancers and it was significantly upregulated in five cancers, which were commonly identifiable in three databases, including breast invasive carcinoma (p < 0.001), kidney chromophobe (p < 0.001), kidney renal clear cell carcinoma (p < 0.001), kidney renal papillary cell carcinoma (p < 0.001), and COAD (p < 0.001). The upregulation was significantly different in association with pathogenic stages II and III (pr(>F) = 0.00234) which was identifiable significantly in COAD more than in other cancers. The gene showed a high upregulation in association with poor prognosis of patient survival in three cancers, including COAD (log-rank p = 0.011), mesothelioma (log-rank p = 5.9e-07), and liver hepatocellular carcinoma (log-rank p = 0.0045). Therefore, COAD was the only comprehensively analyzed tumor to show a diagnostic and prognostic potential characteristic during high upregulation of CDKN2A. Furthermore, CDKN2A displayed a rare mutation in the form of deep deletion (9%) and revealed an upregulation associated with CD4+ T cells (p = 0.0108), macrophage (p = 0.0073), and neutrophils (p = 0.0272) as immune cells infiltrating COAD.  Conclusion Our study demonstrates the pan-cancer relevance of CDKN2A and revealed a novelty in showing CDKN2A underscores its potential as a diagnostic prognostic biomarker in COAD since CDKN2A is mostly studied at a genetic level across COAD.

Identification and validation of a novel stress granules-related prognostic model in colorectal cancer

Aims: A growing body of evidence demonstrates that Stress granules (SGs), a non-membrane cytoplasmic compartments, are important to colorectal development and chemoresistance. However, the clinical and pathological significance of SGs in colorectal cancer (CRC) patients is unclear. The aim of this study is to propose a new prognostic model related to SGs for CRC on the basis of transcriptional expression. Main methods: Differentially expressed SGs-related genes (DESGGs) were identified in CRC patients from TCGA dataset by limma R package. The univariate and Multivariate Cox regression model was used to construct a SGs-related prognostic prediction gene signature (SGPPGS). The CIBERSORT algorithm was used to assess cellular immune components between the two different risk groups. The mRNA expression levels of the predictive signature from 3 partial response (PR) and 6 stable disease (SD) or progress disease (PD) after neoadjuvant therapy CRC patients' specimen were examined. Key findings: By screening and identification, SGPPGS comprised of four genes (CPT2, NRG1, GAP43, and CDKN2A) from DESGGs is established. Furthermore, we find that the risk score of SGPPGS is an independent prognostic factor to overall survival. Notably, the abundance of immune response inhibitory components in tumor tissues is upregulated in the group with a high-risk score of SGPPGS. Importantly, the risk score of SGPPGS is associated with the chemotherapy response in metastatic colorectal cancer. Significance: This study reveals the association between SGs related genes and CRC prognosis and provides a novel SGs related gene signature for CRC prognosis prediction.

TRIM29 is differentially expressed in colorectal cancers of different primary locations and affects survival by regulating tumor immunity based on retrospective study and bioinformatics analysis

Background

In colorectal cancer (CRC) patients, different primary tumor locations caused distinct prognosis and clinicopathological features. It is necessary to identify specific tumor markers according tumor site. Our previous work has identified differentially expressed genes between CRC and adjacent normal tissues, in which only TRIM29 was differently expressed between right colon cancer (RCC) and left colon cancer (LCC) patients. Rectal cancer (RECC) was not included in this latter study and the effects of TRIM29 on the survival with RCC and LCC patients were not investigated. This study further verified TRIM29 expression through Gene Expression Omnibus (GEO) database and our retrospective study. The role of TRIM29 on survival according tumor sites was also explored. Furthermore, the molecular mechanisms of TRIM29 were explored.

Methods

The GEO dataset was used to confirm the differential expression of TRIM29 in proximal and distal cancers. Moreover, TRIM29 were assess using immunohistochemistry (IHC) in 227 cases to observe the correlation between TRIM29 and tumor site. The relationship between TRIM29 and the clinicopathologic features was investigated according tumor sites. Furthermore, the disease-free survival (DFS) and overall survival (OS) was analyzed using the Kaplan-Meier method to assess the prognostic value of TRIM29. Finally, bioinformatics analysis was used to explore the molecular mechanisms. The Tumor-Immune System Interactions and Drug Bank database (TISIDB) was used to analyze the correlations between TRIM29 expression and tumor immune functions. The correlation of TRIM29 with tumor infiltrating lymphocytes or mismatch-repair-proficient/mismatch-repair-deficient (pMMR/dMMR) status was also investigated.

Results

TRIM29 expression was significantly higher in patients with RCC (P<0.001). RCC patients with high TRIM29 tended to be older, male, in stage III-IV, with N+ staging, and intestinal obstruction (P<0.001, P<0.001, P<0.001, P<0.001, and P=0.010, respectively). High TRIM29 expression was associated with an increased risk of recurrence/metastasis and death, only in RCC patients (P=0.020 and P<0.001). Functional annotations and immune activity analysis showed that TRIM29 is related to tumor infiltrating lymphocytes and immune dysfunction.

Conclusions

TRIM29 plays varying roles in patients with different tumor sites. TRIM29 is correlated with the clinicopathological features and prognosis in RCC patients. Indeed, TRIM29 may serve as a new biomarker for RCC patients.

The p53 Pathway Related Genes Predict the Prognosis of Colon Cancer

Background

Colon cancer is a common gastrointestinal malignancy. This study aimed to explore the relationship between p53 pathway-related genes and prognosis of colon cancer.

Methods

The mRNA datasets of colon cancer and adjacent tissues were downloaded from The Cancer Genome Atlas (TCGA) database, and the differential expression of genes in two groups was analyzed. Then, P53 pathway-related genes were intersected with differentially expressed genes (DEGs) to obtain P53 pathway-related differentially expressed genes. Then, overall survival (OS), disease-specific survival (DSS), and progression-free survival (PFS) in clusters were compared by consistent cluster analysis. Univariate and multivariate Cox regression analysis of DEGs was performed to obtain survival-related DEGs. Risk scores were calculated for each sample based on survival-related DEGs, and patients were divided into high/low risk scores. Prognostic differences, tumor immune cell infiltration levels, and immune pathway activation status were compared between the two groups.

Results

We identified 28 DEGs and two clusters. There are significant differences in PFS between the two clusters (P=0.011), and no significant difference between OS and DSS. We obtained 3 DEGs (CDKN2A, BAK1, BTG1) that were significantly related to PFS, and CDKN2A was considered an independent prognostic factor. PFS showed statistically significant difference between high/low risk score groups (P=0.015). There were significant differences in immune cell infiltration level and immune pathway activity between two groups.

Conclusion

The p53 pathway-related genes are significantly related to PFS in colon cancer patients and play an important role in regulating the tumor immune microenvironment.

Identification and validation of hub genes for diabetic retinopathy

Background

Diabetic retinopathy (DR) is characterized by a gradually progressive alteration in the retinal microvasculature that leads to middle-aged adult acquired persistent blindness. Limited research has been conducted on DR pathogenesis at the gene level. Thus, we aimed to reveal novel key genes that might be associated with DR formation via a bioinformatics analysis.

Methods

The GSE53257 dataset from the Gene Expression Omnibus was downloaded for gene co-expression analysis. We identified significant gene modules via the Weighted Gene Co-expression Network Analysis, which was conducted by the Protein-Protein Interaction (PPI) Network via Cytoscape and from this we screened for key genes and gene sets for particular functional and pathway-specific enrichments. The hub gene expression was verified by real-time PCR in DR rats modeling and an external database.

Results

Two significant gene modules were identified. Significant key genes were predominantly associated with mitochondrial function, fatty acid oxidation and oxidative stress. Among all key genes analyzed, six up-regulated genes (i.e., SLC25A33, NDUFS1, MRPS23, CYB5R1, MECR, and MRPL15) were highly and significantly relevant in the context of DR formation. The PCR results showed that SLC25A33 and NDUFS1 expression were increased in DR rats modeling group.

Conclusion

Gene co-expression network analysis highlights the importance of mitochondria and oxidative stress in the pathophysiology of DR. DR co-expressing gene module was constructed and key genes were identified, and both SLC25A33 and NDUFS1 may serve as potential biomarker and therapeutic target for DR.