Chromatin accessibility analysis reveals that TFAP2A promotes angiogenesis in acquired resistance to anlotinib in lung cancer cells

TFAP2A downregulation mediates tumor-suppressive effect of miR-8072 in triple-negative breast cancer via inhibiting SNAI1 transcription

BACKGROUND

Triple-negative breast cancer (TNBC) represents a highly aggressive subset of breast malignancies characterized by its challenging clinical management and unfavorable prognosis. While TFAP2A, a member of the AP-2 transcription factor family, has been implicated in maintaining the basal phenotype of breast cancer, its precise regulatory role in TNBC remains undefined.

METHODS

In vitro assessments of TNBC cell growth and migratory potential were conducted using MTS, colony formation, and EdU assays. Quantitative PCR was employed to analyze mRNA expression levels, while Western blot was utilized to evaluate protein expression and phosphorylation status of AKT and ERK. The post-transcriptional regulation of TFAP2A by miR-8072 and the transcriptional activation of SNAI1 by TFAP2A were investigated through luciferase reporter assays. A xenograft mouse model was employed to assess the in vivo growth capacity of TNBC cells.

RESULTS

Selective silencing of TFAP2A significantly impeded the proliferation and migration of TNBC cells, with elevated TFAP2A expression observed in breast cancer tissues. Notably, TNBC patients exhibiting heightened TFAP2A levels experienced abbreviated overall survival. Mechanistically, TFAP2A was identified as a transcriptional activator of SNAI1, a crucial regulator of epithelial-mesenchymal transition (EMT) and cellular proliferation, thereby augmenting the oncogenic properties of TFAP2A in TNBC. Moreover, miR-8072 was unveiled as a negative regulator of TFAP2A, exerting potent inhibitory effects on TNBC cell growth and migration. Importantly, the tumor-suppressive actions mediated by the miR-8072/TFAP2A axis were intricately associated with the attenuation of AKT/ERK signaling cascades and the blockade of EMT processes.

CONCLUSIONS

Our findings unravel the role and underlying molecular mechanism of TFAP2A in driving tumorigenesis of TNBC. Targeting the TFAP2A/SNAI1 pathway and utilizing miR-8072 as a suppressor represent promising therapeutic strategies for treating TNBC.

Potential prognosis and immunotherapy predictor TFAP2A in pan-cancer

BACKGROUND

TFAP2A is critical in regulating the expression of various genes, affecting various biological processes and driving tumorigenesis and tumor development. However, the significance of TFAP2A in carcinogenesis processes remains obscure.

METHODS

In our study, we explored multiple databases including TCGA, GTEx, HPA, cBioPortal, TCIA, and other well-established databases for further analysis to expound TFAP2A expression, genetic alternations, and their relationship with the prognosis and cellular signaling network alternations. GO term and KEGG pathway enrichment analysis as well as GSEA were conducted to examine the common functions of TFAP2A. RT-qPCR, Western Blot and Dual Luciferase Reporter assay were employed to perform experimental validation.

RESULTS

TFAP2A mRNA expression level was upregulated and its genetic alternations were frequently present in most cancer types. The enrichment analysis results prompted us to investigate the changes in the tumor immune microenvironment further. We discovered that the expression of TFAP2A was significantly associated with the expression of immune checkpoint genes, immune subtypes, ESTIMATE scores, tumor-infiltrating immune cells, and the possible role of TFAP2A in predicting immunotherapy efficacy. In addition, high TFAP2A expression significantly correlated with several ICP genes, and promoted the expression of PD-L1 on mRNA and protein levels through regulating its expression at the transcriptional level. TFAP2A protein level was upregulated in fresh colon tumor tissue samples compared to that in the adjacent normal tissues, which essentially positively correlated with the expression of PD-L1.

CONCLUSIONS

Our study suggests that targeting TFAP2A may provide a novel and effective strategy for cancer treatment.

A rapid and stable spontaneous reprogramming system of Spermatogonial stem cells to Pluripotent State

BACKGROUND

The scarcity of pluripotent stem cells poses a major challenge to the clinical application, given ethical and biosafety considerations. While germline stem cells commit to gamete differentiation throughout life, studies demonstrated the spontaneous acquisition of pluripotency by spermatogonial stem cells (SSCs) from neonatal testes at a low frequency (1 in 1.5 × 107). Notably, this process occurs without exogenous oncogenes or chemical supplementation. However, while knockout of the p53 gene accelerates the transformation of SSCs, it also increases risk and hampers their clinical use.

RESULTS

We report a transformation system that efficiently and stably convert SSCs into pluripotent stem cells around 10 passages with the morphology similar to that of epiblast stem cells, which convert to embryonic stem (ES) cell-like colonies after change with ES medium. Epidermal growth factor (EGF), leukemia inhibitory factor (LIF) and fresh mouse embryonic fibroblast feeder (MEF) are essential for transformation, and addition of 2i (CHIR99021 and PD0325901) further enhanced the pluripotency. Transcriptome analysis revealed that EGF activated the RAS signaling pathway and inhibited p38 to initiate transformation, and synergically cooperated with LIF to promote the transformation.

CONCLUSION

This system established an efficient and safe resource of pluripotent cells from autologous germline, and provide new avenues for regenerative medicine and animal cloning.

Single-cell data analysis of malignant epithelial cell heterogeneity in lung adenocarcinoma for patient classification and prognosis prediction

Lung cancer is one of the leading causes of cancer-related death. Most advanced lung adenocarcinoma (LUAD) patients have poor survival because of drug resistance and relapse. Neglecting intratumoral heterogeneity might be one of the reasons for treatment insensitivity, while single-cell RNA sequencing (scRNA-seq) technologies can provide transcriptome information at the single-cell level. Herein, we combined scRNA-seq and bulk RNA-seq data of LUAD and identified a novel cluster of malignant epithelial cells - KRT81+ malignant epithelial cells - associated with worse prognoses. Further analysis revealed that the hypoxia and EMT pathways of these cells were activated to predispose them to differentiate into metastatic lung adenocarcinoma cells. Finally, we also studied the role of these tumor cells in the immune microenvironment and their role in the classification and prognosis prediction of lung adenocarcinoma patients.

Crucial role of the transcription factors family activator protein 2 in cancer: current clue and views

The transcription factor family activator protein 2 (TFAP2) is vital for regulating both embryonic and oncogenic development. The TFAP2 family consists of five DNA-binding proteins, including TFAP2A, TFAP2B, TFAP2C, TFAP2D and TFAP2E. The importance of TFAP2 in tumor biology is becoming more widely recognized. While TFAP2D is not well studied, here, we mainly focus on the other four TFAP2 members. As a transcription factor, TFAP2 regulates the downstream targets directly by binding to their regulatory region. In addition, the regulation of downstream targets by epigenetic modification, posttranslational regulation, and interaction with noncoding RNA have also been identified. According to the pathways in which the downstream targets are involved in, the regulatory effects of TFAP2 on tumorigenesis are generally summarized as follows: stemness and EMT, interaction between TFAP2 and tumor microenvironment, cell cycle and DNA damage repair, ER- and ERBB2-related signaling pathway, ferroptosis and therapeutic response. Moreover, the factors that affect TFAP2 expression in oncogenesis are also summarized. Here, we review and discuss the most recent studies on TFAP2 and its effects on carcinogenesis and regulatory mechanisms.

TGF-β Signaling Activation Confers Anlotinib Resistance in Gastric Cancer

BACKGROUND

Gastric cancer (GC) has always been a great threat to human health due to its aggressiveness and lethality. Anlotinib, a novel multi-target tyrosine kinase inhibitor (TKI), has been certified its anti-tumor effects on various tumors. Nonetheless, there are few studies on applying anlotinib as a treatment for GC. The underlying mechanism of acquired resistance during anlotinib administration remains unclear.

METHODS

We investigated the toxicologic effects of anlotinib on GC cells through CCK8, colony-forming, and flow cytometry assays in vitro and xenograft models in vivo. Anlotinib-resistant GC cells, AGS-R and MGC803-R, were generated and characterized by cell proliferation and apoptosis assays. The signaling pathways involved in anlotinib resistance were probed using Cignal™ Finder 10-Pathway Reporter Array. Western blot and dual-luciferase reporter assays were performed to confirm the relationships. The TGF-β inhibitor LY364947 was introduced to demonstrate the importance of TGF-β signaling in anlotinib resistance via a series of functional assays.

RESULTS

Anlotinib suppressed cell growth and induced apoptosis in vitro and inhibited tumorigenesis and metastasis in vivo, while its anti-tumor effects were impaired in anlotinib-resistant cells. The results of dual-luciferase reporter assays and western blot indicated TGF-β signaling was activated in anlotinib-resistant GC cells. LY364947 combined with Anlotinib exerted a better antineoplastic effect than monotherapy and considerably reversed the anlotinib resistance in GC.

CONCLUSIONS

Our findings suggested that TGF-β signaling may take a significant part in anlotinib resistance in GC. The suppression of TGF-β signaling may be a possible and promising approach for the GC oncotherapy when combined with anlotinib.

The oncogenic role of TFAP2A in bladder urothelial carcinoma via a novel long noncoding RNA TPRG1-AS1/DNMT3A/CRTAC1 axis

BACKGROUND

Overexpression of TFAP2A has been linked to increased lymph node metastasis in basal-squamous bladder cancer. However, its downstream targets in bladder urothelial carcinoma (BLCA), the most malignant cancer of the urinary tract, remain unclear. In the current study, we aim to explore the function and mechanism of TFAP2A in BLCA.

METHODS

TFAP2A expression and the prognostic significance in BLCA was analyzed using TCGA and GTEX projects. TFAP2A was knocked-down in BLCA cells to study its impact on glucose uptake, lactate and ATP production, expression of HK2, and the number of vascular meshes formed by HUVEC. The target long noncoding RNAs (lncRNAs) of TFAP2A were predicted by bioinformatics tools, followed by ChIP-qPCR and luciferase assays. The downstream targets of TPRG1-AS1 were analyzed by microarray analysis. Rescue experiments were conducted for validation.

RESULTS

TFAP2A upregulation in BLCA predicted dismal survival of patients. Loss of TFAP2A inhibited glycolysis (as evidenced by reduced glucose uptake, lactate, ATP production, and the expression of HK2) and angiogenesis (decreased number of vascular meshes formed by HUVEC). TFAP2A promoted the transcription of TPRG1-AS1. TPRG1-AS1 reversed the inhibitory effect of TFAP2A knockdown on glycolysis and angiogenesis in BLCA cells. TPRG1-AS1 inhibited the transcription of CRTAC1 by recruiting a DNA methyltransferase to the promoter of CRTAC1 and increasing the DNA methylation of its promoter. CRTAC1 inhibited glycolysis and angiogenesis in BLCA cells. TFAP2A silencing curbed tumor growth in vivo via the TPRG1-AS1/CRTAC1 axis.

CONCLUSION

TFAP2A reduces CRTAC1 expression by promoting TPRG1-AS1 transcription, thereby expediting BLCA glycolysis and angiogenesis.

FGFR2-ERC1: A Subtype of FGFR2 Oncogenic Fusion Variant in Lung Adenocarcinoma and the Response to Anlotinib

Background

Fibroblast growth factor receptor (FGFR) fusions in non-small cell lung cancer (NSCLC) are small genomic events. At present, there is no standard treatment strategy for patients with NSCLC carrying an FGFR fusion.

Case Presentation

We report the case of a 45-year-old female patient who was diagnosed with lung adenocarcinoma and underwent right upper lobectomy and postoperative adjuvant chemotherapy. After 13 months, the patient’s lung lesions progressed. Next-generation sequencing of venous blood and lung tissues confirmed an FGFR2-ERC1 fusion, and she received chemotherapy and immunotherapy. Two months later, the patient’s lung lesions progressed again. Based on the target effect of anlotinib on FGFR, the patient was subsequently treated with anlotinib, and the progression-free survival interval exceeded 8.0 months.

Conclusion

These findings showed that patients with lung adenocarcinoma carrying an FGFR2-ERC1 fusion gene may benefit from anlotinib. This case provided evidence to support the use of anlotinib in the treatment of NSCLC patients with FGFR fusion gene subtypes.

Analysis of chromatin accessibility in p53 deficient spermatogonial stem cells for high frequency transformation into pluripotent state

OBJECTIVES

Spermatogonial stem cells (SSCs), the germline stem cells (GSCs) committed to spermatogenesis in niche, can transform into pluripotent state in long-term culture without introduction of exogenous factors, typically in p53 deficiency condition. As the guardian for genomic stability, p53 is associated with epigenetic alterations during SSCs transformation. However, the mechanism is still unknown, since complicated roles of p53 baffle our understanding of the regulating process.

MATERIALS AND METHODS

The chromatin accessibility and differentially expressed genes (DEGs) were analysed in p53+/+ and p53-/- SSCs using the Assay for Transposase-Accessible Chromatin with high-throughput Sequencing (ATAC-seq) and RNA-sequencing (RNA-seq), to explore the connection of p53 and cell fate at chromosomal level.

RESULTS

Several transcription factors (TFs), such as CTCF, SMAD3 and SOX2, were predicted as important factors mediating the transformation. Molecular evidence suggested that SMAD3 efficiently promoted pluripotency-associated gene expression both in fresh and long-term cultured SSCs. However, p53 knockout (KO) is insufficient to induce SMAD3 expression in SSCs.

CONCLUSIONS

These observations indicate that SMAD3 is a key factor for SSCs transformation, and an unknown event is required to activate SMAD3 as the prerequisite for SSCs reprogramming, which may occur in the long-term culture of SSCs. This study demonstrates the connection of p53 and pluripotency-associated factors, providing new insight for understanding the mechanisms of SSCs reprogramming and germline tumorigenesis.

Integrated Chromatin Accessibility and Transcriptome Landscapes of 5-Fluorouracil-Resistant Colon Cancer Cells

Background: 5-Fluorouracil (5-FU) is one of the most effective and widely used chemotherapeutic drugs in the treatment of colon cancer, yet chemoresistance is a common feature of colon cancer treatment, resulting in poor prognosis and short survival. Dynamic reprogramming of chromatin accessibility is crucial for proper regulation of gene transcription associated with cancer drug resistance by providing the gene regulatory machinery with rapid access to the open genomic DNA.

Methods: Here, we explored the global chromatin accessibility and transcription changes by the assay for transposase-accessible chromatin using sequencing (ATAC-seq) in combination with transcriptome sequencing of both parental and 5-FU-resistant HCT15 cells, followed by integrative analysis to better understand the regulatory network underlying 5-FU resistance in colon cancer cells.

Results: A total of 3,175 differentially expressed mRNAs (DEGs), lncRNAs (DELs), and miRNAs (DEMs) related to 5-FU resistance were identified, including significantly upregulated IL33, H19, and miR-17-5p; the downregulated AKR1B10, LINC01012, and miR-125b-5p; and chromatin modifiers such as INO80C, HDAC6, and KDM5A. The construction of the ceRNA regulatory network revealed that H19, HOXA11-AS, and NEAT1 might function as ceRNAs associated with 5-FU resistance in HCT15 cells. Moreover, 9,868 differentially accessible regions (DARs) were obtained, which were positively (r = 0.58) correlated with their nearest DEGs and DELs. The upregulated genes related to 4,937 hyper-accessible regions were significantly enriched in signaling pathways of MAPK, FOX, and WNT, while the 4,931 hypo-accessible regions were considered to be involved in declined biosynthesis of amino acids and nucleotide sugars, signaling pathways of Notch, and HIF-1. Analyses of the DAR sequences revealed that besides the AP-1 family, the TF motifs of FOX and KLF family members were highly enriched in hyper- and hypo-accessible regions, respectively. Finally, we obtained several critical TFs and their potential targets associated with DARs and 5-FU resistance, including FOXA1 and KLF3.

Conclusion: These data provided clear insights and valuable resources for an improved understanding of the non-genetic landscape of 5-FU-resistant colon cancer cells based on chromatin accessibility and transcript levels, which allowed for genome-wide detection of TF binding sites, potential cis-regulatory elements and therapeutic targets.

TP53 Mutation Status and Biopsy Lesion Type Determine the Immunotherapeutic Stratification in Non-Small-Cell Lung Cancer

Immunotherapy, a chemotherapy-free process, has emerged as a promising therapeutic strategy to prolong the overall survival (OS) of patients with non-small-cell lung cancer (NSCLC). However, effective stratification factors for immunotherapy remain unclear. The purpose of this study was to discuss the potential stratification factors of NSCLC immunotherapy using immune checkpoint inhibitors (ICIs) by integrating genomic profiling and tumor lesion–type information. In this study, 344 patients with NSCLC, whose clinical and tissue (including metastatic and primary lesions) mutation information was available, were included. The potential gene mutation status for predicting the outcomes of immunotherapy was screened by comparing the difference in mutation frequency between responders and non-responders. Our results indicated that the potential predictors of immunotherapy were significantly different, especially between patients with TP53(+) (including metastatic and primary lesions) and TP53(−) (including metastatic and primary lesions). According to this classification, patients with NSCLC who suggested immunotherapy had a higher OS than those who did not (25 months vs. 7 months, P < 0.0001, hazard ratio = 0.39). Collectively, this study provides a new perspective for screening immunotherapy predictors in NSCLC, suggesting that the TP53 mutation status and source of biopsy tissue should be considered during the development of immunotherapy biomarkers.

ctDNA-Profiling-Based UBL Biological Process Mutation Status as a Predictor of Atezolizumab Response Among TP53-Negative NSCLC Patients

Atezolizumab, an immune checkpoint inhibitor, has been approved for use in clinical practice in non-small cell lung cancer (NSCLC) patients, but potential biomarkers for response stratification still need further screening. In the present study, a total of 399 patients with high-quality ctDNA profiling results were included. The mutation status of ubiquitin-like conjugation (UBL) biological process genes (including ABL1, APC, LRP6, FUBP1, KEAP1, and TOP2A) and clinical information were further integrated. The results suggested that the patients with the clinical characteristics of male or history of smoking had a higher frequency of UBL mutation positivity [UBL (+)]; the patients who were UBL (+) had shorter progression-free survival (PFS) (1.69 vs. 3.22 months, p = 0.0007) and overall survival (8.61 vs. 16.10 months, p < 0.0001) than those patients with UBL mutation negativity [UBL (–)]; and more promising predictive values were shown in the smoker subgroup and ≤ 3 metastasis subgroup. More interestingly, we found the predictor has more performance in TP53-negative cohorts [training in an independent POPLAR and OAK cohorts (n = 200), and validation in an independent MSKCC cohort (n = 127)]. Overall, this study provides a predictor, UBL biological process gene mutation status, not only for identifying NSCLC patients who may respond to atezolizumab therapy but also for screening out the potential NSCLC responders who received other immune checkpoint inhibitors.

Integrated Chromatin Accessibility and Transcriptome Landscapes of Doxorubicin-Resistant Breast Cancer Cells

BACKGROUND

Doxorubicin is one of the most effective chemotherapeutic drugs for breast cancer while its common drug resistance leads to poor patient prognosis and survival. Growing evidence indicate dynamically reorganized chromatin allows rapid access of the gene regulatory machinery to open genomic regions facilitating subsequent gene expression through direct transcription factor (TF) activation and regulatory element binding.

METHODS

To better understand the regulatory network underlying doxorubicin resistance in breast cancer cells, we explored the systematic alterations of chromatin accessibility and gene expression by the assay for transposase-accessible chromatin using sequencing (ATAC-seq) in combination with RNA sequencing, followed by integrative analysis to identify potential regulators and their targets associated with differentially accessible regions (DARs) in doxorubicin-resistant MCF7 (MCF7-DR) cells.

RESULTS

A total of 3,963 differentially expressed genes (DEGs) related to doxorubicin resistance were identified, including dramatically up-regulated MT1E, GSTP1, LDHB, significantly down-regulated TFF1, UBB, DSCAM-AS1, and histone-modifying enzyme coding genes HDAC2, EZH2, PRMT5, etc. By integrating with transcriptomic datasets, we identified 18,228 DARs in MCF7-DR cells compared to control, which were positively correlated with their nearest DEGs (r = 0.6). There were 11,686 increased chromatin-accessible regions, which were enriched in up-regulated genes related to diverse KEGG pathways, such as the cell cycle, regulation of actin cytoskeleton, signaling pathways of MAPK, PI3K/Akt and Hippo, which play essential roles in regulating cell apoptosis, proliferation, metabolism, and inflammatory responses. The 6,542 decreased chromatin-accessible regions were identified for the declined doxorubicin-associated biological processes, for instance, endocrine and insulin resistance, central carbon metabolism, signaling pathways of TGF-beta and P53. Combining data from TCGA, analyses of the DAR sequences associated with the DNA-binding motifs of significantly enriched TF families including AP-1, TEAD and FOX, indicated that the loss-function of FOXA1 might play a critical role in doxorubicin-resistant breast cancer cells (DOX-R BCCs).

CONCLUSION

These data exhibit the non-genetic landscape of chromatin accessibility and transcript levels in the DOX-R BCCs, and provide clear insights and resources for the detection of critical TFs and potential cis-regulatory elements-based putative therapeutic targets.

Systematic Chromatin Accessibility Analysis Based on Different Immunological Subtypes of Clear Cell Renal Cell Carcinoma

Background

Recent research of clear cell renal cell carcinoma (ccRCC) is focused on the tumor immune microenvironment (TIME). Chromatin accessibility is critical for regulation of gene expression. However, its role in different immunological subtypes of ccRCC based on immune cell infiltration has not been systematically studied.

Methods

Five hundred thirty patient data from The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) were adopted to estimate immune cell infiltration. Twenty-four types of immune cells were evaluated with single-sample Gene Set Enrichment Analysis (ssGSEA). Patients were divided into two clusters based on immune cell infiltration. Systematic chromatin accessibility analysis was conducted based on the two clusters.

Results

We compared the relative expression of the immune gene signatures among 530 patients of TCGA-KIRC using ssGSEA. Overall survival (OS) analysis revealed 10 types of immune cells were significantly associated with prognosis. Patients were divided into two clusters based on 24 types of immune cell infiltration. Immune cell signals as well as PD-1/PD-L1 signal were higher in cluster 1. Among the two clusters, 2,400 differential peaks were found in TCGA-KIRC Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) data. The distribution of differential peaks and prognosis-related immune cells in 23 chromosomes are essentially the same. There is no peak distribution downstream. The proportion of peaks upstream of the 5’ transcription start site decreases, and both sides of binding regions of the TSS 0.1-1 kb becomes smaller. Enrichment analysis of GO and KEGG of these differential peaks showed that they are remarkably related to the immune regulation in tumor microenvironment. Known motifs and de novo motifs were found by linking motif annotations to different peaks. Survival analysis of related motif transcription factors were prognostic. The GSEA enrichment analysis showed that high SP1 expression positively correlates with TGF-beta signaling and inflammatory response, while negatively correlates with TNF-alpha signaling via NFKB. High KLF12 expression negatively correlates with interferon gamma response, IL2-STAT5 signaling, TNF-alpha signaling via NFKB, IL6-JAK-STAT3 signaling.

Conclusion

The abnormality of chromatin accessibility may play an important regulatory role in ccRCC immunity.

PlGF knockdown attenuates hypoxia-induced stimulation of cell proliferation and glycolysis of lung adenocarcinoma through inhibiting Wnt/β-catenin pathway

Background

Angiogenic placental growth factor (PlGF) plays a role in hypoxia-induced angiogenesis. Here, we aimed to investigate the biological roles of PlGF in cell proliferation and glycolysis of lung adenocarcinoma (LUAD) and the underlying molecular mechanisms.

Methods

PlGF was knocked down in H358 and H1975 cells by lentiviruses, which were then cultured under hypoxia (90% N₂, 5%CO₂ and 5%O₂) for 24 h. PlGF was overexpressed in PC9 cells treated with XAV939, inhibitor of Wnt/β-catenin signaling pathway. PlGF-silencing H1975 cells were implanted into mice, and tumor xenografts were harvested and analyzed.

Results

Hypoxia treatment led to up-regulation of PlGF, C-myc, lactate dehydrogenase A (LDHA), and β-catenin, promotion of cell proliferation and glycolysis in H358 and H1975 cells, which were obviously reversed by knocking down PlGF. In tumors, PlGF knockdown significantly prohibited cell proliferation and glycolysis, and decreased expression of C-myc, LDHA, and β-catenin. PlGF overexpression markedly strengthened cell proliferation, which was inhibited by β-catenin knockdown. Consistently, XAV939, inhibitor of Wnt/β-catenin pathway, also inhibited PlGF-induced cell proliferation, glycolysis, and β-catenin expression in PC9 cells.

Conclusion

PlGF knockdown inhibited the stimulatory effect of hypoxia on cell proliferation and glycolysis of LUAD through deactivating Wnt/β-catenin pathway.

A Prognostic Model Based on the Immune-related Genes in Colon Adenocarcinoma

Background: Immune-related genes (IRGs) are critically involved in the tumor microenvironment (TME) of colon adenocarcinoma (COAD). Here, the study was mainly designed to establish a prognostic model of IRGs to predict the survival of COAD patients. Methods: The Cancer Genome Atlas (TCGA), Immunology Database and Analysis Portal (ImmPort) database, and Cistrome database were utilized for extracting data regarding the expression of immune gene- and tumor-related transcription factors (TFs), aimed at the identification of differentially expressed genes (DEGs), differentially expressed IRGs (DEIRGs), and differentially expressed TFs (DETFs). Univariate Cox regression analysis was subsequently performed for the acquisition of prognosis-related IRGs, followed by establishment of TF regulatory network for uncovering the possible molecular regulatory association in COAD. Subsequently, multivariate Cox regression analysis was conducted to further determine the role of prognosis-related IRGs for prognostic prediction in COAD. Finally, the feasibility of a prognostic model with immunocytes was explored by immunocyte infiltration analysis. Results: A total of 2450 DEGs, 8 DETFs, and 79 DEIRGs were extracted from the corresponding databases. Univariate Cox regression analysis revealed 11 prognosis-related IRGs, followed by establishment of a regulatory network on prognosis-related IRGs at transcriptional levels. Functionally, IRG GLP2R was negatively modulated by TF MYH11, whereas IRG TDGF1 was positively modulated by TF TFAP2A. Multivariate Cox regression analysis was subsequently performed to establish a prognostic model on the basis of seven prognosis-related IRGs (GLP2R, ESM1, TDGF1, SLC10A2, INHBA, STC2, and CXCL1). Moreover, correlation analysis of immunocyte infiltration also revealed that the seven-IRG prognostic model was positively associated with five types of immunocytes (dendritic cell, macrophage, CD4 T cell, CD8 T cell, and neutrophil), which may directly reflect tumor immune state in COAD. Conclusions: Our present findings indicate that the prognostic model based on prognosis-related IRGs plays a crucial role in the clinical supervision and prognostic prediction of COAD patients at both molecular and cellular levels.