Enhanced DNA Repair Pathway is Associated with Cell Proliferation and Worse Survival in Hepatocellular Carcinoma (HCC)

Angiogenesis Is Associated With Aggressive Biology That Counterbalances With Tumor Immunogenicity in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is an arterialized tumor; thus, anti-angiogenesis targeted therapy is in clinical practice. Herein, we hypothesized that HCC with high angiogenesis is biologically aggressive with worse survival.

Methods

Angiogenesis score (AS) was derived from the Molecular Signatures Database (MSigDB) Hallmark Angiogenesis Gene Set, and median was used to divide high versus low groups. Transcriptome of HCC patients of The Cancer Genome Atlas (TCGA, n = 386) and GSE76427 (n = 115) cohorts were analyzed.

Results

High AS correlated with angiogenesis-related gene expressions. Both microvascular and lymphatic endothelial cell infiltrations were higher in high angiogenesis HCC. Surprisingly, no survival difference was seen with varying levels of angiogenesis. High angiogenesis significantly enriched tumor aggravating signaling pathways: glycolysis, Notch, Hedgehog, KRAS, epithelial mesenchymal transition, and transforming growth factor-beta (TGF-β) in Gene Set Enrichment Analysis (GSEA), but also infiltrated less CD8+ T cells and T-helper 1 cells, and higher M1 macrophages and conventional dendritic cells (cDCs) with elevated cytolytic activity score in both cohorts. In agreement, immune response-related gene sets: inflammatory response, tumor necrosis factor-alpha (TNF-α) signaling, allograft rejection, interferon-alpha, and interferon-gamma were all enriched to high angiogenesis HCC. Programmed cell death protein 1 (PD1), programmed death ligand 1 (PD-L1), programmed death ligand 2 (PD-L2), and cytotoxic T lymphocyte-associated protein 4 (CTLA-4) were higher in high angiogenesis HCC in TCGA, but not in GSE76427 cohort.

Conclusions

Angiogenesis quantified using transcriptome of HCC patients demonstrated that it is associated with aggressive biology but also with tumor immunogenicity and immune response that counterbalance and did not reflect in survival. Given high expression of immune checkpoint molecules, we cannot help but speculate that immunotherapy may be useful for high angiogenesis HCC patients.

Heterogeneity in polyamine metabolism dictates prognosis and immune checkpoint blockade response in hepatocellular carcinoma

Immune checkpoint blockade holds promise in hepatocellular carcinoma (HCC) treatment, but its efficacy remains limited. Dysregulated polyamine metabolism and its interaction with oncogenic pathways promote tumor progression. However, the heterogeneity of polyamine metabolism and its effects on the immune microenvironment and response to immunotherapy in HCC remain unclear. Here, we aimed to investigate the prognostic and immunotherapeutic implications of polyamine metabolism in HCC. Based on polyamine-related genes, HCC patients were categorized into two clusters with distinct survival outcomes. We developed a polyamine-related signature, termed PAscore, which was found to be a strong predictor of both poor prognosis and reduced immunocyte infiltration. Notably, a high PAscore was also associated with decreased sensitivity to immunotherapy. Within the HCC microenvironment, malignant cells exhibited polyamine metabolic heterogeneity, those with high polyamine metabolic activity showed altered hallmark pathway signatures and increased communication with myeloid cells. In vitro experiments suggested that FIRRE, the gene with the greatest impact on the PAscore, significantly contributed to HCC proliferation and metastasis. This study underscores the potential of our polyamine-related signature in predicting the prognosis and immunotherapy response in HCC patients, and also reveals the polyamine metabolic heterogeneity among HCC cells that influences their crosstalk with infiltrating myeloid cells.

An Integrated Framework to Identify Prognostic Biomarkers and Novel Therapeutic Targets in Hepatocellular Carcinoma-Based Disabilities

Hepatocellular carcinoma (HCC) is one of the most prevalent malignant tumors globally, significantly affecting liver functions, thus necessitating the identification of biomarkers and effective therapeutics to improve HCC-based disabilities. This study aimed to identify prognostic biomarkers, signaling cascades, and candidate drugs for the treatment of HCC through integrated bioinformatics approaches such as functional enrichment analysis, survival analysis, molecular docking, and simulation. Differential expression and functional enrichment analyses revealed 176 common differentially expressed genes from two microarray datasets, GSE29721 and GSE49515, significantly involved in HCC development and progression. Topological analyses revealed 12 hub genes exhibiting elevated expression in patients with higher tumor stages and grades. Survival analyses indicated that 11 hub genes (CCNB1, AURKA, RACGAP1, CEP55, SMC4, RRM2, PRC1, CKAP2, SMC2, UHRF1, and FANCI) and three transcription factors (E2F1, CREB1, and NFYA) are strongly linked to poor patient survival. Finally, molecular docking and simulation identified seven candidate drugs with stable complexes to their target proteins: tozasertib (-9.8 kcal/mol), tamatinib (-9.6 kcal/mol), ilorasertib (-9.5 kcal/mol), hesperidin (-9.5 kcal/mol), PF-562271 (-9.3 kcal/mol), coumestrol (-8.4 kcal/mol), and clofarabine (-7.7 kcal/mol). These findings suggest that the identified hub genes and TFs could serve as valuable prognostic biomarkers and therapeutic targets for HCC-based disabilities.

Enhanced cancer cell proliferation and aggressive phenotype counterbalance in breast cancer with high BRCA1 gene expression

PURPOSE

While comprehensive research exists on the mutation of the DNA repair gene BRCA1, limited information is available regarding the clinical significance of BRCA1 gene expression. Given that cancer cell proliferation is aggrevated by DNA repair, we hypothesized that high BRCA1 gene expression breast cancer (BC) might be linked with aggressive tumor biology and poor clinical outcomes.

METHODS

The cohorts: The Cancer Genome Atlas (TCGA, n = 1069), METABRIC (n = 1903), and SCAN-B (n = 3273) were utilzed to obtain data of 6245 BC patients.

RESULTS

BC patients without BRCA1 mutation exhibited higher BRCA1 expression, which was associated with DNA repair functionality. However, no such correlation was observed with BRCA2 expression. The association of high BRCA1 expression with cancer cell proliferation was evidenced by significant enrichment of cell proliferation-related gene sets, higher histological grade, and proliferation score. Furthermore, increased levels of homologous recombination deficiency, intratumoral heterogeneity, and altered fractions were associated with high BRCA1 expression. Moreover, BC with high BRCA1 expression exhibited reduced infiltration of dendritic cells and CD8 T-cells, while showing increased infiltration of Th1 cells. Surprisingly, BRCA1 expression was not associated with the survival of BC irrespective of the subtypes. Conversely, BC with low BRCA1 expression enriched cancer aggravating pathway gene sets, such as Cancer Stem Cell-related signaling (NOTCH and HEDGEHOG), Angiogenesis, Epithelial-Mesenchymal Transition, Inflammatory Response, and TGF-beta signaling.

CONCLUSION

Despite being linked to heightened proliferation of cancer cells and unassertive phenotype, BRCA1 expression did not show any association with survival in BC.

Breast cancer in adolescents and young adults has a specific biology and poor patient outcome compared with older patients

BACKGROUND

We aimed to clarify the features of adolescents and young adults (AYA: younger than 40 years old) breast cancer (BC) compared with other age groups in estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative BC, given the effects of age-related hormonal status.

METHODS

The cohorts analyzed were divided into AYA (15-39 years old), perimenopausal (40-54 years old), menopausal (55-64 years old), and old (65+ years old). Clinicopathological and biological features were analyzed using gene set variation analysis and xCell algorithm using transcriptome profiles from large public databases of ER-positive/HER2-negative BC (METABRIC; n = 1353, SCAN-B; n = 2381).

RESULTS

In the ER-positive/HER2-negative subtype, pathological lymph node positivity, and Nottingham grade 3 were higher among AYA (all P < 0.001). AYA patients had a trend toward worse disease-specific and overall survival, particularly compared with the perimenopausal group. Estrogen response late signaling decreased with age (all P ≤ 0.001 in both METABRIC and SCAN-B cohorts). AYA was associated with significantly higher BRCAness and DNA repair than the other groups (all P < 0.05 in both cohorts). AYA significantly enriched cell proliferation-related and procancerous gene sets [mTORC1, unfolded protein response, and phosphoinositide 3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/AKT/mTOR) signaling] when compared with the others (all P < 0.03 in both cohorts). Interestingly, these features have also been observed in tumors <2 cm. Infiltration of CD8+, regulatory, T helper type 2 cells, and M1 macrophages was higher, while M2 macrophages were lower in AYA (all P < 0.03 in both cohorts). Finally, ER-positive/HER2-negative BC in AYA patients has different features of gene mutations, including AHNAK2, GATA3, HERC2, and TG, which were observed at a higher rate in AYA, and KMT2C, which was observed at a lower rate in AYA, compared with other age groups.

CONCLUSIONS

ER-positive/HER2-negative BC in AYA was highly proliferative with high immune cell infiltration compared with the other age groups.

RAD51 High-Expressed Hepatocellular Carcinomas Are Associated With High Cell Proliferation

INTRODUCTION

RAD51 is a pivotal DNA repair gene managing double-stranded DNA break recognition and repair. RAD51 high expression was associated with adverse outcomes in other cancer types. This study aims to investigate the tumor microenvironment and immune landscape in the RAD51 high-expressed Hepatocellular Carcinoma (HCCs).

METHODS

A total of 467 patients from two large independent cohorts with clinical and transcriptomic data were obtained. The cohort was dichotomized based on the median RAD51 gene expression. xCell and Gene Set Enrichment Analysis (GSEA) were used.

RESULTS

RAD51 high-expressed HCCs were associated with worse recurrence-free, progression-free, disease-specific, and overall survival (all P < 0.05). While RAD51 high-expressed HCCs were associated with intratumoral heterogeneity, homologous recombination deficiency, and fraction altered scores, mutation or neoantigens were not increased in this group. xCell analysis demonstrated inconsistent immune cell infiltration between two cohorts. Cytolytic activity as well as GSEA with immune-related gene sets also demonstrated inconsistent results between two cohorts as well. On the other hand, RAD51 expression was significantly increased in higher-grade tumors, larger tumors, and higher clinical stages. RAD51 high-expressed HCCs were found to have elevated proliferation score. Furthermore, GSEA exhibited significant enrichment of all the cell proliferation-related gene sets in the Hallmark collection, including E2F targets, G2M checkpoint, Mitotic spindle, MYC targets, and MTORC1 signaling consistently in both cohorts (all false discovery rate < 0.25).

CONCLUSIONS

RAD51 high-expressed HCCs were associated with worse survival and with increased cell proliferation and were not necessarily associated with immune infiltration or inflammation.

KIN17 functions in DNA damage repair and chemosensitivity by modulating RAD51 in hepatocellular carcinoma

The limited response of hepatocellular carcinoma (HCC) to chemotherapy drugs has always been a bottleneck in therapy. DNA damage repair is a major reason for chemoresistance. Previous studies have confirmed that KIN17 affects chemosensitivity. In this study, we examined the impact of KIN17 on chemotherapy response and DNA repair in HCC cells treated with oxaliplatin (L-OHP). We evaluated the expression and biological roles of KIN17 in HCC using bioinformatic analysis. The correlation between KIN17 and RAD51, particularly their nuclear expression levels, was evaluated using immunofluorescence, immunoblotting after nucleocytoplasmic separation in HCC cells, and immunohistochemistry of mouse xenograft tumors and human HCC tissues. The results indicated a significant increase in KIN17 expression in HCC tissues compared to normal tissues. The GSEA analysis revealed that upregulation of KIN17 was significantly associated with DNA damage repair. Knockdown of KIN17 led to increased DNA damage and reduced cellular survival after exposure to L-OHP. On the other hand, overexpression of KIN17 was linked to decreased DNA damage and improved cell survival following L-OHP treatment. Further experiments indicated that KIN17 affects the expression of RAD51, particularly in the nucleus. KIN17 plays a crucial role in influencing the sensitivity of HCC to chemotherapy by triggering the DNA repair response. Increased expression of KIN17 is associated with a poor prognosis for HCC patients, indicating that KIN17 could serve as a prognostic marker and therapeutic target for HCC.

Higher inflammatory response in hepatocellular carcinoma is associated with immune cell infiltration and a better outcome

BACKGROUND AND AIMS

Hepatocellular carcinoma (HCC) often develops from chronic liver inflammation. Inflammation within a tumor can either promote cancer progression or activate an immune response against it. This study aims to determine the clinical significance of enhanced inflammation in HCC.

METHODS

Data from 655 HCC patients across four cohorts (TCGA, GSE6764, GSE76427, GSE89377) were examined. Inflammatory response was quantified using a scoring system derived from the gene set variation analysis of the "INFLAMMATORY_RESPONSE" gene set.

RESULTS

A stepwise increase in inflammatory response was noted from normal liver to cirrhosis, with consistently lower levels in HCC across both GSE6764 and GSE89377 cohorts (both p < 0.001). Similar trends were observed in interferon response, pathways such as IL6/JAK/STAT3 and complement signaling, coagulation cascade, and allograft rejection (all p < 0.02). HCCs with high inflammatory response were associated with increased immune cell infiltrations (p < 0.01) and cytolytic activity (p < 0.001). Interestingly, these HCCs had reduced mutation rates, no relationship with cell proliferation, and displayed both immune responses and pro-cancerous signals including epithelial-mesenchymal transition, KRAS, and hypoxia. Further, a high inflammatory score correlated with improved disease-free survival in TCGA (p = 0.034) and overall survival in GSE76427 (p = 0.008).

CONCLUSION

HCC with higher levels of inflammatory response demonstrated increased immune cell infiltration, enhanced immune-related and other pro-cancerous-related signaling, and showed a trend toward a better patient prognosis.

Eribulin induces micronuclei and enhances the nuclear localization of cGAS in triple-negative breast cancer cells

Eribulin (ERI), clinically utilized for locally advanced or metastatic breast tumors, has shown potential links to the immune system. Notably, the cGAS-STING pathway, a key component of innate immunity, has gained prominence. Yet, limited reports explore ERI's effects on the cGAS-STING pathway. Additionally, the nuclear presence of cGAS remains poorly understood. This study uniquely delves into ERI's impact on both the cytosolic cGAS-STING pathway and nuclear cGAS. ERI enhances nuclear localization of cGAS, resulting in hyper-activation of the cGAS-STING pathway in triple-negative breast cancer cells. Reduction of cGAS heightened both cell proliferation and ERI sensitivity. In clinical data using ERI in a neo-adjuvant setting, patients with low cGAS cases exhibited reduced likelihood of achieving pathological complete response after ERI treatment. These findings illuminate the potential of cGAS and IFNβ as predictive biomarkers for ERI sensitivity, providing valuable insights for personalized breast cancer treatment strategies.

A novel subtype based on driver methylation-transcription in lung adenocarcinoma

AIMS

To identify driver methylation genes and a novel subtype of lung adenocarcinoma (LUAD) by multi-omics and elucidate its molecular features and clinical significance.

METHODS

We collected LUAD patients from public databases, and identified driver methylation genes (DMGs) by MethSig and MethylMix algrothms. And novel driver methylation multi-omics subtypes were identified by similarity network fusion (SNF). Furthermore, the prognosis, tumor microenvironment (TME), molecular features and therapy efficiency among subtypes were comprehensively evaluated.

RESULTS

147 overlapped driver methylation were identified and validated. By integrating the mRNA expression and methylation of DMGs using SNF, four distinct patterns, termed as S1-S4, were characterized by differences in prognosis, biological features, and TME. The S2 subtype showed unfavorable prognosis. By comparing the characteristics of the DMGs subtypes with the traditional subtypes, S3 was concentrated in proximal-inflammatory (PI) subtype, and S4 was consisted of terminal respiratory unit (TRU) subtype and PI subtype. By analyzing TME and epithelial mesenchymal transition (EMT) features, increased immune infiltration and higher expression of immune checkpoint genes were found in S3 and S4. While S4 showed higher EMT score and expression of EMT associated genes, indicating S4 may not be as immunosensitive as the S3. Additionally, S3 had lower TIDE and higher IPS score, indicating its increased sensitivity to immunotherapy.

CONCLUSION

The driver methylation-related subtypes of LUAD demonstrate prognostic predictive ability that could help inform treatment response and provide complementary information to the existing subtypes.

Enhanced epithelial-mesenchymal transition signatures are linked with adverse tumor microenvironment, angiogenesis and worse survival in gastric cancer

Epithelial-mesenchymal transition (EMT) is a crucial mechanism that facilitates cancer cell metastasis. Despite its importance, the clinical significance of EMT in gastric cancer (GC) patients has yet to be clearly demonstrated. For gauging the extent of EMT in GC, we employed gene set variation analysis to score 807 patient samples from two large cohorts: TCGA and GSE84437. In both cohorts, EMT high GC showed a significant association with worse overall survival (hazard ratio (HR) = 1.74, p = 0.011 and HR = 2.01, p < 0.001, respectively). This association was stronger when considering the EMT signature score compared to the individual expressions of EMT-related genes (CDH1, CDH2, VIM, and FN1). While the EMT signature level did not differ among various cancers, high EMT signature specifically correlated with survival in GC alone. Mucinous and diffuse histological types exhibited higher EMT levels compared to others (p < 0.001), and the EMT signature level was correlated with tumor depth and AJCC stage (all p < 0.001). Interestingly, the EMT score was an independent factor for overall and disease-specific survival (multivariate; p = 0.006 and 0.032, respectively). EMT high GC displayed a lower fraction of Th1 cells and a higher fraction of dendritic cells, M1 macrophages and several stromal cells. EMT high GC exhibited an inverse correlation with cell proliferation-related gene sets. While they significantly enriched multiple pro-cancerous gene sets, such as TGF-β signaling, hypoxia, and angiogenesis. The presence of EMT signature in a bulk tumor was linked to TGF-β signaling, hypoxia, and angiogenesis, and was also associated with poorer survival outcomes in GC patients.

Spinster Homologue 2 Expression Correlates With Improved Patient Survival in Hepatocellular Carcinoma Despite Association With Lymph-Angiogenesis

BACKGROUND

Spinster homologue 2 (SPNS2) is a transporter of sphingosine-1-phosphate (S1P), a bioactive lipid linked to cancer progression. We studied the link between SPNS2 gene expression, tumor aggressiveness, and outcomes in patients with hepatocellular carcinoma (HCC).

METHODS

Gene expression in patients with HCC was analyzed from the Cancer Genome Atlas (TCGA) (n = 350) and GSE76427 (n = 115) as a validation cohort, as well as liver tissue cohort GSE6764 (n = 75).

RESULTS

High-SPNS2 HCC was significantly associated with high level of lymph-angiogenesis-related factors. SPNS2 expression was significantly higher in normal liver and early HCC versus advanced HCC (P < 0.02). High SPNS2 levels enriched immune response-related gene sets; inflammatory, interferon (IFN)-α, IFN-γ responses, and tumor necrosis factor (TNF)-α, interleukin (IL)-6/Janus kinase/signal transducer and activator of transcription (JAK/STAT3) signaling, complement and allograft rejection, but did not significantly infiltrate specific immune cells nor cytolytic activity score. High-SPNS2 HCC enriched tumor aggravating pathway gene sets such as KRAS (Kirsten rat sarcoma virus) signaling, but inversely correlated with Nottingham histological grade, MKI67 (marker of proliferation Ki-67) expression, and cell proliferation-related gene sets. Further, high-SPNS2 HCC had significantly high infiltration of stromal cells, showing that low-SPNS2 HCC is highly proliferative. Finally, high-SPNS2 HCC was associated with better disease-free, disease-specific, and overall survival (P = 0.031, 0.046, and 0.040, respectively).

CONCLUSIONS

Although SPNS2 expression correlated with lymph-angiogenesis and other cancer-promoting pathways, it also enriched immune response. SPNS2 levels were higher in normal liver compared to HCC, and inversely correlated with cancer cell proliferation and better survival. SPNS2 expression may be beneficial in HCC patients despite detrimental in-vitro effects.

High Ki67 Gene Expression Is Associated With Aggressive Phenotype in Hepatocellular Carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) with high Ki67 protein expression, the most commonly used cell proliferation marker, is associated with an aggressive biologic phenotype; however, conventional immunostaining is hampered by variability in institutional protocol, specific antibody probe, and by assessor subjectivity. To this end, we hypothesized that Ki67 gene (MKi67) expression would identify highly proliferative HCC, and clarify its association with oncologic outcome, tumor progression, and immune cell population in the tumor microenvironment (TME). Furthermore, we sought to identify the cell-cycle gene expression profile that confers this aggressive phenotype.

METHODS

A total of 473 HCC patients with clinicopathological data associated with transcriptome were selected for this study: 358 patients from The Cancer Genome Atlas (TCGA) as the testing cohort, and 115 from GSE76427 as the validation cohort. Each cohort was divided into a highly proliferative group (MKi67-high) and the low MKi67 group (MKi67-low) by the median of Ki67 gene (MKi67) expression levels.

RESULTS

MKi67-high HCC patients had worse disease-free survival (DFS), disease-specific survival (DSS), and overall survival (OS) independent of histological grade in the TCGA cohort. MKi67 expression correlated with histological grade and tumor size. MKi67 expression increased throughout the HCC carcinomatous sequence from normal liver, cirrhotic liver, early HCC, and advanced HCC. MKi67-high HCC was associated with higher intratumor heterogeneity, homologous recombination deficiency, and altered fraction as well as intratumoral infiltration of T helper type 1 (Th1) and Th2 cells, but lower interferon-gamma response and M2 macrophage infiltration. Cell proliferation-related gene sets in the Hallmark collection (E2F targets, G2M checkpoint, Myc target v1 and mitotic spindle), MTORC1 signaling, DNA repair, PI3K MTOR signaling, and unfolded protein response were all enriched in the MKi67-high HCC (false discovery rate (FDR) < 0.25).

CONCLUSIONS

High MKi67 gene expression identified highly proliferative HCC with aggressive biology involving classical pathways in cell cycle regulation and DNA repair, as well as poor overall oncologic outcomes. This suggests potential for personalized treatment strategies, but validation and refinement of these observations require further research to elucidate the underlying mechanisms and validate therapeutic targeting of these pathways in MKi67-high HCC tumors.

Higher Inflammatory Response in Hepatocellular Carcinoma is Associated with Immune Cell Infiltration and a Better Outcome

Background & Aims

Hepatocellular carcinoma (HCC) often develops from chronic liver inflammation. Inflammation within a tumor can either promote cancer progression or activate an immune response against it. This study aims to determine the clinical significance of enhanced inflammation in HCC.

Methods

Data from 655 HCC patients across four cohorts (TCGA, GSE6764, GSE76427, GSE89377) were examined. Inflammatory response was quantified using a scoring system derived from the gene set variation analysis of the "INFLAMMATORY_RESPONSE" gene set.

Results

A stepwise increase in inflammatory response was noted from normal liver to cirrhosis, with consistently lower levels in HCC across both GSE6764 and GSE89377 cohorts (both p<0.001). Similar trends were observed in interferon response, pathways such as IL6/JAK/STAT3 and complement signaling, coagulation cascade, and allograft rejection (all p<0.02). HCCs with high inflammatory response were associated with increased immune cell infiltrations (p<0.01) and cytolytic activity (p<0.001). Interestingly, these HCCs had reduced mutation rates, no relationship with cell proliferation, and displayed both immune responses and pro-cancerous signals including epithelial-mesenchymal transition, KRAS, and hypoxia. Further, a high inflammatory score correlated with improved disease-free survival in TCGA (p=0.034) and overall survival in GSE76427 (p=0.008).

Conclusion

HCC with higher levels of inflammatory response demonstrated increased immune cell infiltration, enhanced immune-related and other pro-cancerous-related signaling, and better patient prognosis.

E2F target score is associated with cell proliferation and survival of patients with hepatocellular carcinoma

BACKGROUND

E2F target genes are essential for the cell cycle. A score that quantifies its activity is expected to reflect the aggressiveness and prognosis of hepatocellular carcinoma.

METHODS

Cohorts of hepatocellular carcinoma patients (total n = 655) from The Cancer Genome Atlas, GSE89377, GSE76427, and GSE6764 were analyzed. The cohorts were divided into high versus low by the median.

RESULTS

All the Hallmark cell proliferation-related gene sets were consistently enriched in hepatocellular carcinoma with high E2F targets score, and E2F score was associated with grade, tumor size, American Joint Committee on Cancer staging, proliferation score, and MKI67 expression, as well as with less abundance of hepatocytes and stromal cells. E2F targets enriched DNA repair, mTORC1 signaling, glycolysis, and unfolded protein response gene sets and were significantly associated with the higher intratumoral genomic heterogeneity, homologous recombination deficiency, and progression of hepatocellular carcinoma. On the other hand, there was no relationship between E2F targets and mutation rates or neoantigens. High E2F hepatocellular carcinoma did not enrich any of the immune-response-related gene sets but was associated with high infiltration of Th1, Th2 cells, and M2 macrophage; however, there was no difference in cytolytic activity. In both early (I and II) and late (III and IV) stages of hepatocellular carcinoma, a high E2F score was associated with worse survival and was an independent prognostic factor for overall and disease-specific survival in patients with hepatocellular carcinoma.

CONCLUSION

The E2F target score, associated with cancer aggressiveness and worse survival, could be used as a prognostic biomarker in patients with hepatocellular carcinoma.

Accelerated glycolysis in tumor microenvironment is associated with worse survival in triple-negative but not consistently with ER+/HER2- breast cancer

Metabolic reprogramming to sustain immortality is a hallmark of cancer and glycolysis is an important way to attain this. Thus, we investigate the association of glycolysis and associated pathways in the survival of breast cancer. A total of 5,176 breast cancer patients from multiple independent cohorts were analyzed. We determined the glycolytic signaling score by the degree of enrichment by Gene Set Variant Analysis and the median was used to divide each cohort into high vs low score groups. Glycolysis high breast cancer significantly enriched the hallmark cell proliferation-related gene sets (E2F targets, G2M checkpoint, and MYC targets v1 and v2) and was associated with high MKI67 expression. In all cohorts, triple-negative breast cancer (TNBC) was associated with the highest glycolysis score. It was found that in TNBC, glycolysis high breast cancer was associated with worse survival but in ER-positive/HER2-negative breast cancer this was not observed consistently. The glycolysis high TNBC enriched multiple pro-cancerous gene sets and was infiltrated with a low level of B-cells and anti-cancerous immune cells, and significantly associated with a decreased level of cytolytic activity. It was also observed that the glycolysis was higher in the metastatic sites than in the primary breast cancer and the survival was not affected by the metastatic sites. In conclusion, accelerated glycolysis is associated with cancer cell proliferation and worse survival in TNBC.

Autophagy orchestrates resistance in hepatocellular carcinoma cells

Treatment resistance is one of the major barriers for therapeutic strategies in hepatocellular carcinoma (HCC). Many studies have indicated that chemotherapy and radiotherapy induce autophagy machinery (cell protective autophagy) in HCC cells. In addition, many experiments report a remarkable crosstalk between treatment resistance and autophagy pathways. Thus, autophagy could be one of the key factors enabling tumor cells to hinder induced cell death after medical interventions. Therefore, extensive research on the molecular pathways involved in resistance induction and autophagy have been conducted to achieve the desired therapeutic response. The key molecular pathways related to the therapy resistance are TGF-β, MAPK, NRF2, NF-κB, and non-coding RNAs. In addition, EMT, drug transports, apoptosis evasion, DNA repair, cancer stem cells, and hypoxia could have considerable impact on the hepatoma cell's response to therapies. These mechanisms protect tumor cells against various treatments and many studies have shown that each of them is connected to the molecular pathways of autophagy induction in HCC. Hence, autophagy inhibition may be an effective strategy to improve therapeutic outcome in HCC patients. In this review, we further highlight how autophagy leads to poor response during treatment through a complex molecular network and how it enhances resistance in primary liver cancer. We propose that combinational regimens of approved HCC therapeutic protocols plus autophagy inhibitors may overcome drug resistance in HCC therapy.

Cancer Stem Cell and Hepatic Stellate Cells in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the most common liver cancer. It is highly lethal and has high recurrence. Death among HCC patients occur mainly due to tumor progression, recurrence, metastasis, and chemoresistance. Cancer stem cells (CSCs) are cell subpopulations within the tumor that promote invasion, recurrence, metastasis, and drug resistance. Hepatic stellate cells (HSCs) are important components of the tumor microenvironment (TME) responsible for primary secretory ECM proteins during liver injury and inflammation. These cells promote fibrogenesis, infiltrate the tumor stroma, and contribute to HCC development. Interactions between HSC and CSC and their microenvironment help promote carcinogenesis through different mechanisms. This review summarizes the roles of CSCs and HSCs in establishing the TME in primary liver tumors and describes their involvement in HCC chemoresistance.

Development of a novel BRCAness score that predicts response to PARP inhibitors

BACKGROUND

BRCAness is a characteristic feature of homologous recombination deficiency (HRD) mimicking BRCA gene mutation in breast cancer. We hypothesized that a measure to quantify BRCAness that causes synthetic lethality in BRCA mutated tumors will identify responders to PARP inhibitors.

METHODS

A total of 6753 breast cancer patients from 3 large independent cohorts were analyzed. A score was generated by transcriptomic profiling using gene set variation analysis algorithm on 34 BRCA1-mutation related genes selected by high AUC levels in ROC curve between BRCA1 mutation and wildtype breast cancer.

RESULTS

The score was significantly associated with BRCA1 mutation, high mutation load and intratumoral heterogeneity as expected, as well as with high HRD, DNA repair and MKi67 expression regardless of BRCA mutations. High BRCAness tumors enriched not only DNA repair, but also all five Hallmark cell proliferation-related gene sets. High BRCAness tumors were significantly associated with higher cytolytic activity and with higher anti-cancerous immune cell infiltration. Not only did the breast cancer cell lines with BRCA-mutation show high score, but even the other cells in human breast cancer tumor microenvironment were contributing to the score. The BRCAness score was the highest in triple-negative breast cancer consistently in all 3 cohorts. BRCAness was associated with response to chemotherapy and correlated strongly with response to PARP inhibitor in both triple-negative and ER-positive/HER2-negative breast cancer.

CONCLUSIONS

We established a novel BRCAness score using BRCA-mutation-related gene expressions and found that it associates with DNA repair and predicts response to PARP inhibitors regardless of BRCA mutation.

Analyses of Transcriptomics Cell Signalling for Pre-Screening Applications in the Integrated Approach for Testing and Assessment of Non-Genotoxic Carcinogens

With recent rapid advancement of methodological tools, mechanistic understanding of biological processes leading to carcinogenesis is expanding. New approach methodologies such as transcriptomics can inform on non-genotoxic mechanisms of chemical carcinogens and can be developed for regulatory applications. The Organisation for the Economic Cooperation and Development (OECD) expert group developing an Integrated Approach to the Testing and Assessment (IATA) of Non-Genotoxic Carcinogens (NGTxC) is reviewing the possible assays to be integrated therein. In this context, we review the application of transcriptomics approaches suitable for pre-screening gene expression changes associated with phenotypic alterations that underlie the carcinogenic processes for subsequent prioritisation of downstream test methods appropriate to specific key events of non-genotoxic carcinogenesis. Using case studies, we evaluate the potential of gene expression analyses especially in relation to breast cancer, to identify the most relevant approaches that could be utilised as (pre-) screening tools, for example Gene Set Enrichment Analysis (GSEA). We also consider how to address the challenges to integrate gene panels and transcriptomic assays into the IATA, highlighting the pivotal omics markers identified for assay measurement in the IATA key events of inflammation, immune response, mitogenic signalling and cell injury.

Identification of DAPK1 as an autophagy-related biomarker for myotonic dystrophy type 1

Myotonic dystrophy type I (DM1), a CTG repeat expansion hereditary disorder, is primarily characterized by myotonia. Several studies have reported that abnormal autophagy pathway has a close relationship with DM1. However, the underlying key regulatory molecules dictating autophagy disturbance still remains elusive. Previous studies mainly focused on finding targeted therapies for DM1, but the clinical heterogeneity of the DM1 is rarely addressed. Herein, to identify potential regulator genes related to autophagy and cross-correlation among clinical symptoms, we performed weighted gene co-expression network analysis (WGCNA) to construct the co-expression network and screened out 7 core autophagy-related genes (DAPK1, KLHL4, ERBB3, SESN3, ATF4, MEG3, and COL1A1) by overlapping within differentially expressed genes (DEG), cytoHubba, gene significance (GS) and module membership (MM) score. Meanwhile, we here analyzed autophagy-related molecular subtypes of DM1 in relation to the clinical phenotype. Our results show that three genes (DAPK1, SESN3, and MEG3) contribute to distinguish these two molecular subtypes of DM1. We then develop an analysis of RNA-seq data from six human skin fibroblasts (3 DM1, 3 healthy donors). Intriguingly, of the 7 hallmark genes obtained, DAPK1 is the only confirmed gene, and finally identified in vitro by RT-PCR. Furthermore, we assessed the DAPK1 accuracy diagnosis of DM1 by plotting a receiver operating characteristic curve (ROC) (AUC = 0.965). In this study, we first validated autophagy status of DM1 individuals exhibits a clearly heterogeneity. Our study identified and validated DAPK1 serve as a novel autophagy-related biomarker that correlate with the progression of DM1.

Comprehensive analysis of the clinical significance, immune infiltration, and biological role of MARCH ligases in HCC

The membrane‐associated RING‐CH (MARCH) family, a member of the E3 ubiquitin ligases, has been confirmed by a growing number of studies to be associated with immune function and has been highlighted as a potential immunotherapy target. In our research, hepatocellular carcinoma (HCC) patients were divided into C1 and C2 MARCH ligase-related patterns by the non-negative matrix factorization (NMF) algorithm. Multiple analyses revealed that the MARCH ligase-related cluster was related to prognosis, clinicopathological characteristics, and the tumor immune microenvironment (TIME). Next, the signature (risk score) of the MARCH prognosis was constructed, including eight genes associated with the MARCH ligase (CYP2C9, G6PD, SLC1A5, SPP1, ANXA10, CDC20, PON1, and FTCD). The risk score showed accuracy and stability. We found that the correlations between risk score and TIME, tumor mutation burden (TMB), prognosis, and clinicopathological characteristics were significant. Additionally, the risk score also had important guiding significance for HCC treatment, including chemotherapy, immunotherapy, and transarterial chemoembolization (TACE).

DNA Damage Repair Classifier Defines Distinct Groups in Hepatocellular Carcinoma

Simple Summary

DNA repair pathways have been implicated in hepatocellular carcinoma outcomes. We found that hepatocellular carcinomas (HCC) could be separated into two groups (high and low) based on the overall expression of genes involved in DNA repair. Among the low repair group, there were three subgroups, one of which shared features of the high repair group. Given the important role of liver in metabolism and detoxification and its regenerative capacity, proliferation and DNA damage responses are critical in subdividing major biological categories of liver tumors. High repair samples showed more proliferative and regenerative signatures and had poorer outcomes versus the low repair that were more associated with the genes involved in normal liver biology. These biological groups suggest that dysregulation in endogenous liver processes promotes a pro-tumorigenic microenvironment that may facilitate tumor progression or identify tumors that require more substantial clinical intervention.

Abstract

DNA repair pathways have been associated with variability in hepatocellular carcinoma (HCC) clinical outcomes, but the mechanism through which DNA repair varies as a function of liver regeneration and other HCC characteristics is poorly understood. We curated a panel of 199 genes representing 15 DNA repair pathways to identify DNA repair expression classes and evaluate their associations with liver features and clinicopathologic variables in The Cancer Genome Atlas (TCGA) HCC study. We identified two groups in HCC, defined by low or high expression across all DNA repair pathways. The low-repair group had lower grade and retained the expression of classical liver markers, whereas the high-repair group had more clinically aggressive features, increased p53 mutant-like gene expression, and high liver regenerative gene expression. These pronounced features overshadowed the variation in the low-repair subset, but when considered separately, the low-repair samples included three subgroups: L1, L2, and L3. L3 had high DNA repair expression with worse progression-free (HR 1.24, 95% CI 0.81–1.91) and overall (HR 1.63, 95% CI 0.98–2.71) survival. High-repair outcomes were also significantly worse compared with the L1 and L2 groups. HCCs vary in DNA repair expression, and a subset of tumors with high regeneration profoundly disrupts liver biology and poor prognosis.

A score of DNA damage repair pathway with the predictive ability for chemotherapy and immunotherapy is strongly associated with immune signaling pathway in pan-cancer

DNA damage repair (DDR) is critical in maintaining normal cellular function and genome integrity and is associated with cancer risk, progression, and therapeutic response. However, there is still a lack of a thorough understanding of the effects of DDR genes’ expression level in cancer progression and therapeutic resistance. Therefore, we defined a tumor-related DDR score (TR-DDR score), utilizing the expression levels of 20 genes, to quantify the tumor signature of DNA damage repair pathways in tumors and explore the possible function and mechanism for the score among different cancers. The TR-DDR score has remarkably predictive power for tumor tissues. It is a more accurate indicator for the response of chemotherapy or immunotherapy combined with the tumor-infiltrating lymphocyte (TIL) and G2M checkpoint score than the pre-existing predictors (CD8 or PD-L1). This study points out that the TR-DDR score generally has positive correlations with patients of advanced-stage, genome-instability, and cell proliferation signature, while negative correlations with inflammatory response, apoptosis, and p53 pathway signature. In the context of tumor immune response, the TR-DDR score strongly positively correlates with the number of T cells (CD4+ activated memory cells, CD8+ cells, T regs, Tfh) and macrophages M1 polarization. In addition, by difference analysis and correlation analysis, COL2A1, MAGEA4, FCRL4, and ZIC1 are screened out as the potential modulating factors for the TR-DDR score. In summary, we light on a new biomarker for DNA damage repair pathways and explore its possible mechanism to guide therapeutic strategies and drug response prediction.

Higher intra-tumoral expression of pro-coagulation genes is a predictor of angiogenesis, epithelial mesenchymal transition and worse patient survival in gastric cancer

Coagulation regulates angiogenesis in cancer, and is associated with tumor development and metastasis. To date, there have been no studies quantifying the state of intra-tumoral coagulation. We measured intra-tumoral coagulation gene expression using the "Hallmark-COAGULATION" gene set in the MSigDB, performing gene set variation analysis and then assigning a "coagulation score" to quantify gene expression. Clinical, histologic, and genetic data were analyzed in 807 gastric cancer patients from the TCGA_STAD and GSE84437 databases. Tumors with increased expression of pro-coagulation genes were consistently associated with higher AJCC T-categories (p = 0.018), lymph node metastasis (p = 0.036), and stage (p = 0.006) in both cohorts. Patients with high coagulation scores were found to have worse disease-specific survival and overall survival (OS) (p = 0.019 and 0.011, respectively) in TCGA, and worse OS in GSE84437 cohort (p = 0.012). Higher expression of pro-coagulation genes correlated with increased intra-tumoral angiogenesis, as well as increased proportions of lymphatic and microvascular endothelial cells, endothelial cells, and pericytes, calculated by xCell algorithm. High coagulation scores were significantly associated with low tumor mutation burden, but not with intratumor heterogeneity and homologous recombination deficiency. Gastric cancers with high coagulation scores contained higher amounts of M1 macrophages and dendritic cells, and low numbers of Th1 cells (all P<0.001). Genes for epithelial mesenchymal transition (EMT), myogenesis, apical junction, transforming growth factor (TGF)-β signaling, and angiogenesis were enriched in high coagulation score-gastric cancers (all false discovery rate <0.25). In conclusion, gastric cancers expressing higher levels of pro-coagulation genes demonstrate increased angiogenesis, EMT, TGF-β signaling and worse patient prognosis.

Enhanced immune response outperform aggressive cancer biology and is associated with better survival in triple-negative breast cancer

Although the value of tumor-infiltrating lymphocytes is well known, the clinical relevance of an increased immune response, specifically in breast cancer, has not been investigated across large cohorts of patients using computational algorithms. Our hypothesis stated that an enhanced immune response is associated with an improvement in outcomes. To quantify the immune response, we utilized the allograft rejection score correlated with cytolytic activity and with all the other Hallmark immune-related gene sets. The score reflected the amount of infiltrating immune cells that correlated with the immune checkpoint molecule expressions, including CD4⁺ and CD8⁺ T cells, T helper type 1 (Th1) and type 2 (Th2) cells, M1 macrophages, B cells, and plasmacytoid dendritic cells (pDC). A high score was associated with high levels of intratumor heterogeneity, homologous recombination defects, mutation rate, histological grade, advanced stage, and lymph node metastasis. Breast malignancy with a high score enriched immune-related gene sets and pro-cancer-related gene sets, including epithelial–mesenchymal transition and KRAS pathway, in ER-positive/HER2-negative and triple-negative breast cancer (TNBC) groups. TNBC had the highest score compared to other subtypes, and was associated with better survival. In conclusion, we found that breast cancer with a high immune response is associated with aggressive cancer biology, but with better survival in TNBC.

The upregulation of stromal antigen 3 expression suppresses the phenotypic hallmarks of hepatocellular carcinoma through the Smad3-CDK4/CDK6-cyclin D1 and CXCR4/RhoA pathways

BACKGROUND

The stromal antigen 3 (STAG3) gene encodes an adhesion complex subunit that can regulate sister chromatid cohesion during cell division. Chromosome instability caused by STAG3 gene mutation may potentially promote tumor progression, but the effect of STAG3 on hepatocellular carcinoma (HCC) and the related molecular mechanism are not reported in the literature. The mechanism of the occurrence and development of HCC is not adequately understood. Therefore, the biological role of STAG3 in HCC remains to be studied, and whether STAG3 might be a sensitive therapeutic target in HCC remains to be determined.

METHODS

The expression and clinical significance of STAG3 in HCC tissues and cell lines were determined by RT-qPCR and immunohistochemistry analyses. The biological functions of STAG3 in HCC were determined through in vitro and in vivo cell function tests. The molecular mechanism of STAG3 in HCC cells was then investigated by western blot assay.

RESULTS

The mRNA expression of STAG3 was lower in most HCC cells than in normal cells. Subsequently, an immunohistochemical analysis of STAG3 was performed with 126 samples, and lower STAG3 expression was associated with worse overall survival in HCC patients. Moreover, cytofunctional tests revealed that the lentivirus-mediated overexpression of STAG3 in HCC cells inhibited cell proliferation, migration, and invasion; promoted apoptosis; induced G1/S phase arrest in vitro; and inhibited tumor growth in vivo. Furthermore, studies of the molecular mechanism suggested that the overexpression of STAG3 increased Smad3 expression and decreased CDK4, CDK6, cyclin D1, CXCR4 and RhoA expression.

CONCLUSION

STAG3 exhibits anticancer effects against HCC, and these effects involve the Smad3-CDK4/CDK6-cyclin D1 and CXCR4/RhoA pathways. STAG3 is a tumor-suppressor gene that may serve as a potential target for molecular therapy, which provides a new idea for the treatment of HCC.

Abundance of reactive oxygen species (ROS) is associated with tumor aggressiveness, immune response, and worse survival in breast cancer

PURPOSE

Reactive oxygen species (ROS) are oxygen-containing molecules that have high reactivity and play roles in protection or harm the cancer cells. We aimed to clarify the clinical relevance of ROS in breast cancer (BC) tumor microenvironment (TME). We hypothesized that it is associated with worse BC patient outcomes.

METHODS

ROS score was generated by Gene Set Variation Analysis of Hallmark ROS pathway gene set and a total of 6245 BC patients were analyzed.

RESULTS

High ROS BC significantly enriched cell proliferation-related gene sets (MYC targets v1 and v2, G2M checkpoint, E2F targets), pro-cancer-related gene sets (DNA repair, unfolded protein response, MTORC1 signaling, PI3K/AKT/MTOR signaling, glycolysis, and oxidative phosphorylation), immune-related gene sets (inflammatory response, allograft rejection, interferon-α and γ responses, complement, and IL6/JAK/STAT3 signaling), and infiltrated immune cells (CD4⁺ memory and CD8⁺ T cells, Th1 and Th2, dendritic cells, Tregs, M1 and M2 macrophages) and B cells, as well as elevated cytolytic activity consistently in both METABRIC and GSE96058 cohorts. Cancer cells were the major source of ROS in BC TME of single-cell sequence (GSE75688) cohort. High ROS was associated with intratumor heterogeneity, homologous recombination defects, mutation rates, and neoantigens, and with clinical aggressiveness in AJCC stage, Nottingham grade and Ki67 expression, as well as worse overall survival in both GSE96058 and METABRIC, and with worse disease-specific survival in METABRIC.

CONCLUSION

Abundant ROS in BC patients is associated with abundant mutations, aggressive cancer biology, immune response, and worse survival.

The clinical relevance of unfolded protein response signaling in breast cancer

Protein homeostasis regulated by the Endoplasmic Reticulum (ER) is a recognized process involved in cancer progression. ER stress activates the Unfolded Protein Response (UPR) and has been implicated in a variety of cancers. Given the role of the UPR activation in carcinogenesis, we hypothesized that UPR activation could be associated with pathological progression, higher clinical stage, and worse survival in breast cancer. A total of 4,416 breast cancer patients from multiple independent cohorts were analyzed. We defined the UPR pathway score by the degree of enrichment by Gene Set Variant Analysis and median was used to divide high vs. low score groups in each cohort. High UPR breast cancer significantly enriched not only cell proliferation-related but also other pro-cancerous gene sets consistently in both METABIC and GSE96058 cohort. Majority of UPR pathway score high cells in the bulk tumor were tumor cells compared to other cells, including stromal, T-, B-, and myeloid-cells (P<0.001). UPR score was significantly associated with advanced stage, high grade, and triple negative breast cancer (TNBC) (all P<0.001). High UPR breast cancer was associated with worse patient survival in both cohorts (all P<0.001). Among breast cancer subtype, ER-positive/HER2-negative breast cancer with high UPR was significantly associated with worse survival, but neither HER-positive nor TNBC. High UPR ER-positive/HER2-negative breast cancer was infiltrated with high level of Th1 and Th2 cells, M1 macrophage, and plasma cells. On the other hand, they were significantly infiltrated with high level of several types of stromal cells in tumor microenvironment (all P<0.001). Finally, high UPR metastatic breast cancer was also associated with worse patient survival (P=0.041). UPR signaling is associated with cancer aggressiveness, and worse survival, especially ER-positive/HER2-negative breast cancer subtype.

Diagnostic accuracy of contrast-enhanced ultrasound for microvascular invasion of hepatocellular carcinoma: A protocol for systematic review and meta-analysis

BACKGROUND

Microvascular invasion is an independent risk factor for the recurrence of hepatocellular carcinoma (HCC). Early detection and timely treatment can reduce the recurrence and prolong the overall survival of HCC. Contrast-enhanced ultrasound (CEUS) has been validated for their predictive potential of microvascular invasion in HCC patients, although the conclusion remains controversial. Therefore, a meta-analysis is necessary to be conducted to validate the diagnostic value of CEUS for microvascular invasion in HCC, thus supporting guideline development and clinical practice.

METHODS

Relevant studies reporting the potential diagnostic value of CEUS for microvascular invasion in HCC patients published before February 2022 will be searched in the PubMed, EMBASE, Cochrane Library, and Web of Science. Data will be extracted by 2 researchers independently. The risk of bias will be assessed by the Quality Assessment of Diagnostic Accuracy Studies-2. All of the above statistical analysis will be carried out with Stata 14.0.

RESULTS

The results of this study will be published in a peer-reviewed journal.

CONCLUSION

This study will comprehensively evaluate the diagnostic accuracy of CEUS for microvascular invasion in HCC patients, thus providing high-quality evidence to support clinical practice and guideline development.

Drug Resistance and Endoplasmic Reticulum Stress in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most common and deadly cancers worldwide. It is usually diagnosed in an advanced stage and is characterized by a high intrinsic drug resistance, leading to limited chemotherapeutic efficacy and relapse after treatment. There is therefore a vast need for understanding underlying mechanisms that contribute to drug resistance and for developing therapeutic strategies that would overcome this. The rapid proliferation of tumor cells, in combination with a highly inflammatory microenvironment, causes a chronic increase of protein synthesis in different hepatic cell populations. This leads to an intensified demand of protein folding, which inevitably causes an accumulation of misfolded or unfolded proteins in the lumen of the endoplasmic reticulum (ER). This process is called ER stress and triggers the unfolded protein response (UPR) in order to restore protein synthesis or—in the case of severe or prolonged ER stress—to induce cell death. Interestingly, the three different arms of the ER stress signaling pathways have been shown to drive chemoresistance in several tumors and could therefore form a promising therapeutic target. This review provides an overview of how ER stress and activation of the UPR contributes to drug resistance in HCC.

Low RUFY3 expression level is associated with lymph node metastasis in older women with invasive breast cancer

PURPOSE

Sentinel lymph node biopsy is omitted in older women (≥ 70 years old) with clinical lymph node (LN)-negative hormone receptor-positive breast cancer as it does not influence adjuvant treatment decision-making. However, older women are heterogeneous in frailty while the chance of recurrence increase with improving longevity. Therefore, a biomarker that identifies LN metastasis may facilitate treatment decision-making. RUFY3 is associated with cancer progression. We evaluated RUFY3 expression level as a biomarker for LN-positive breast cancer in older women.

METHODS

Clinical and transcriptomic data of breast cancer patients were obtained from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC, n = 1903) and The Cancer Genome Atlas (TCGA, n = 1046) Pan-cancer study cohorts.

RESULTS

A total of 510 (METABRIC) and 211 (TCGA) older women were identified. LN-positive breast cancer, which represented 51.4% (METABRIC) and 48.4% (TCGA), demonstrated worse disease-free, disease-specific, and overall survival. RUFY3 levels were significantly lower in LN-positive tumors regardless of age. The area under the curve for the receiver operator characteristic (AUC-ROC) curves showed RUFY3-predicted LN metastasis. Low RUFY3 enriched oxidative phosphorylation, DNA repair, MYC targets, unfolded protein response, and mtorc1 signaling gene sets, was associated with T helper type 1 cell infiltration, and with intratumor heterogeneity and fraction altered. Low RUFY3 expression was associated with LN-positive breast cancer and with worse disease-specific survival among older women.

CONCLUSION

Older women with breast cancers who had low expression level of RUFY3 were more frequently diagnosed with LN-positive tumors, which translated into worse prognosis.

Low intratumoral genetic neutrophil-to-lymphocyte ratio (NLR) is associated with favorable tumor immune microenvironment and with survival in triple negative breast cancer (TNBC)

Patients with triple negative breast cancer (TNBC) have a poor prognosis. A novel prognostic biomarker may guide management by appropriately selecting patients for particular treatments. Peripheral blood neutrophil-to-lymphocyte ratio (NLR) was reported to associate with cancer progression, thus we hypothesized that intratumor genetic NLR will reflect tumor immune microenvironment (TIME) and breast cancer biology. The intratumoral genetic NLR previously defined as the ratio of CD66b (CEACAM8) and CD8 (CD8A) gene expressions was utilized to analyze total of 2,994 patients from METABRIC, TCGA, GSE21094, GSE22358, GSE25088, GSE32646, and GSE2603 cohorts. Intratumoral genetic NLR did not correlate with cancer stage nor clinical parameters of cancer cell proliferation such as Nottingham histological grade or MKI67 expression levels in neither the METABRIC or TCGA cohorts. Intratumoral genetic NLR-high breast cancer was not associated with pathologic complete response (pCR) after neoadjuvant chemotherapy in 5 independent cohorts with different regimens. Despite these results, intratumoral genetic NLR-high TNBC demonstrated worse disease-free, disease-specific, and overall survival. Intratumoral genetic NLR-low TNBC enriched multiple immune-related gene sets, was associated with higher favorable immune-related scores and with a favorable TIME, whereas no gene sets enriched to NLR-high TNBC. In conclusion, intratumoral genetic NLR-low TNBC was associated with favorable TIME and with better survival.

Development of KAM score to predict metastasis and worse survival in breast cancer

Some may think that prediction of metastasis is meaningless since metastatic breast cancer is currently incurable. We argue that effective identification of developing metastasis will enable us to design and conduct clinical trials specifically targeting those patients at high risk. The current study sought to generate the KAM score by 4 genes (BRSK2, EYA1, SIGLEC15, and AGTR1) overexpressed in primary breast cancer that developed metastasis to bone compared with matched controls without metastasis longer than 10 years. A high KAM score was prognostic of poor overall (OS), disease free survival (DFS), and disease specific survival (DSS) in the METABRIC, and OS in the GSE96058 cohorts. A high KAM score was significantly associated with clinical aggressiveness, such as high American Joint Committee Cancer (AJCC) stage, lymph node metastasis, Nottingham pathological grade, and triple negative breast cancer (TNBC). Subgroup analysis revealed that a high KAM score was associated with worse OS in ER-positive/HER2-negative breast cancer in both cohorts. A high KAM breast cancer enriched all 5 cell proliferation-related gene sets of the Hallmark collection and interferon (IFN)-γ response gene sets. Furthermore, a high KAM breast cancer was significantly infiltrated with a high fraction of not only anti-cancer but also pro-cancer immune cells and associated with high level of cytolytic activity. Finally, a high KAM breast cancer was significantly associated with lung metastasis. In conclusion, we developed KAM score using 4 gene expressions that predict lung metastasis and patient survival in breast cancer.

CCNB1 promotes the development of hepatocellular carcinoma by mediating DNA replication in the cell cycle

In our studies, cyclin B1 (CCNB1) mRNA and protein were overexpressed in hepatocellular carcinoma (HCC) tissues compared with non-HCC tissues. Moreover, CCNB1 was overexpressed in the serum of HCC patients. The expression of CCNB1 was associated with several crucial clinicopathologic characteristics, and the HCC patients with overexpressed CCNB1 had worse overall survival outcomes. In the screening of interactional genes, a total of 266 upregulated co-expression genes, which were positively associated with CCNB1, were selected from the datasets, and 67 downregulated co-expression genes, which were negatively associated with CCNB1, were identified. The key genes might be functionally enriched in DNA replication and the cell cycle pathways. CDC20, CCNA2, PLK1, and FTCD were selected for further research because they were highly connected in the protein-protein interaction networks. Upregulated CDC20, CCNA2, and PLK1 and downregulated FTCD might result in undesirable overall survival outcomes for HCC patients. The univariate Cox analysis results showed that CDC20 and PLK1 might be two independent risk factors, while FTCD might be protective in HCC. Therefore, CCNB1 may participate in the cell cycle of HCC by regulating DNA replication, and CCNB1 may provide a direction for the diagnosis of early-stage HCC and targeted HCC therapy.

Conflicting roles of EGFR expression by subtypes in breast cancer

Epidermal growth factor receptor (EGFR) is one of the receptors that belong to the epidermal growth factor family of receptor tyrosine kinases (ErbBs). Several malignancies including breast cancer that express EGFR have poor prognosis. Our study examined the EGFR expression among 5176 breast cancer patients from GSE96058 and METABRIC cohorts and the contribution of tumor immune microenvironment in different subtypes. We found that among different breast cancer subtypes, EGFR expression in TNBC was the highest compared to other subtypes. EGFR high ER-positive/HER2-negative breast cancer had significantly higher survival compared to EGFR low ER-positive/HER2-negative breast cancer. It was also associated with high level of intratumor heterogeneity and homologous recombination defects (HRD). This group was also enriched in immune-related gene sets. On the other hand, low EGFR tumor was enriched in cell proliferation-related gene sets. However, these findings were not observed in TNBC. Interestingly, there was a greater infiltration of anti-cancer immune cells in high EGFR ER-positive/HER2-negative breast cancers, while, TNBC with higher EGFR expression had lower fraction of immune cells along with low level of cytolytic activity. Tumor cells have significantly higher EGFR expression compared to immune cells in single cell sequencing data. There was higher expression of immune checkpoint molecules in high EGFR ER-positive/HER2-negative breast cancer but lower expression in TNBC. High EGFR metastatic tumor was significantly associated with worse survival, but no association with infiltrating immune cells was observed. Our study shows that higher EGFR expression in ER-positive/HER2-negative breast cancer is associated with improved outcomes and an anti-cancer immune microenvironment.

Sphingosine 1-phosphate (S1P) produced by sphingosine kinase 1 (SphK1) and exported via ABCC1 is related to hepatocellular carcinoma (HCC) progression

Sphingosine-1-Phosphate (S1P) is produced by Sphingosine Kinase 1 (SphK1) in the cell and is transported out of the cells by ABCC1 transporter. S1P induces inflammation, angiogenesis and modulates tumor immune microenvironment (TIME) in autocrine and paracrine manner. We hypothesized that high S1P export is associated with hepatocellular carcinoma (HCC) progression and worse survival. Transcriptome linked with clinical data were obtained from a total of 533 patients from TCGA (The Cancer Genome Atlas)-HCC (n = 350), GSE6764 (n = 75), and GSE89377 (n = 108) cohorts. Both SphK1 and ABCC1 were expressed higher in aggressive HCC than normal liver or cirrhosis and correlated with MKi67 expression. High S1P export by high expression of both SphK1 and ABCC1 enriched gene sets related with cell proliferation (E2F targets, G2M checkpoint, MYC targets), inflammation (Inflammatory response, TNFα, IL6), angiogenesis, metastasis (TGF-β, epithelial-mesenchymal transition), and immune response (allograft rejection, complement, interferon-gamma) in gene set enrichment analysis. High S1P export was associated with elevation of HGF, HSP90AA1, TRAF2, and AKR1B10. It was also associated with high intratumor heterogeneity, leucocyte fraction, macrophage regulation and lymphocyte infiltration, as well as T helper type2 cells, macrophages, dendritic cells, CD4⁺ T memory activated cells, B-cells and cytolytic activity score in TIME. High S1P export was associated with significantly worse disease specific survival (P = 0.034) and overall survival (P = 0.004) compared to low S1P export group. In conclusion, simultaneous high expression of SphK1 and ABCC1 that reflect S1P export is associated with enhancement of both HCC progression and immune response. Given that S1P export was also associated with worse survival, we cannot help but speculate that pro-cancer pathways activated by S1P may overwhelm the anti-cancer immune response mediated by S1P.

Increased intratumor heterogeneity, angiogenesis and epithelial to mesenchymal transition pathways in metaplastic breast cancer

Metaplastic breast cancer (MBC) constitutes a rare but unique histologic entity with poor prognosis. We hypothesized that MBC possesses unique genetic profile and tumor immune microenvironment. MBC cases were identified from a total of 10827 breast cancer entries in the Cancer Genome Atlas Data Set (TCGA) and the AACR-GENIE (Genomics Evidence Neoplasia Information Exchange) cohorts. Tumor infiltrated immune cells were estimated by xCell. Baseline clinical characteristics were compared, and gene set enrichment analysis (GSEA) was performed. MBC comprised 0.66% of the cohorts (1.2% of TCGA and 0.6% of GENIE). MBC cases were predominantly triple-negative (TNBC) (8 (61.5%) vs 151 (14.4%), P<0.001), and high Nottingham histological grade (8 (61.5%) vs 222 (21.1%), P=0.02) compared to non-MBC in the TCGA cohort. Increased infiltration of M1 macrophages (P=0.012), dendritic cells (P<0.001) and eosinophils (P=0.036) was noted in the MBC cohort however there was no difference in cytolytic activity (P=0.806), CD4 memory (P=0.297) or CD8 T-cells (P=0.864). Tumor mutation burden was lower in the MBC compared to the non-MBC, median: 0.4 vs 1.6/Mb in the TCGA-TNBC cohort (P=0.67) and 3.0 vs 4.0/Mb (P=0.1) in the GENIE-cohort. MBC had increased intratumor heterogeneity (P<0.001), macrophage regulation (P=0.008) and TGF-beta response (P<0.001). Disease-specific survival was decreased in MBC (P=0.018). Angiogenesis and epithelial-to-mesenchymal transition pathways were enriched in triple-negative MBC by GSEA (P=0.004 and P<0.001, respectively). Our results suggest that high intratumor heterogeneity, enriched angiogenesis and EMT pathway expression represent possible mechanisms leading to worse disease-specific survival found in metaplastic breast cancer.

The Unfolded Protein Response Is Associated with Cancer Proliferation and Worse Survival in Hepatocellular Carcinoma

Simple Summary

We studied the association between the unfolded protein response (UPR) and carcinogenesis, cancer progression, and survival in hepatocellular carcinoma (HCC). We studied 655 HCC patients from 4 independent cohorts using an UPR score. The UPR was enhanced as normal liver became cancerous and as HCC advanced in stage. The UPR was correlated with cancer cell proliferation that was confirmed by multiple parameters. Significantly, a high UPR score was associated with worse patient survival. Interestingly, though UPR was associated with a high mutational load, it was not associated with immune response, immune cell infiltration, or angiogenesis. To our knowledge, this is the first study to investigate the clinical relevance of the unfolded protein response in HCC.

Abstract

Hepatocellular carcinoma is a leading cause of cancer death worldwide. The unfolded protein response (UPR) has been revealed to confer tumorigenic capacity in cancer cells. We hypothesized that a quantifiable score representative of the UPR could be used as a biomarker for cancer progression in HCC. In this study, a total of 655 HCC patients from 4 independent HCC cohorts were studied to examine the relationships between enhancement of the UPR and cancer biology and patient survival in HCC utilizing an UPR score. The UPR correlated with carcinogenic sequence and progression of HCC consistently in two cohorts. Enhanced UPR was associated with the clinical parameters of HCC progression, such as cancer stage and multiple parameters of cell proliferation, including histological grade, mKI67 gene expression, and enrichment of cell proliferation-related gene sets. The UPR was significantly associated with increased mutational load, but not with immune cell infiltration or angiogeneis across independent cohorts. The UPR was consistently associated with worse survival across independent cohorts of HCC. In conclusion, the UPR score may be useful as a biomarker to predict prognosis and to better understand HCC.

Construction of a prognostic signature of autophagy-related lncRNAs in non-small-cell lung cancer

Background

Autophagy inhibits tumorigenesis by limiting inflammation. LncRNAs regulate gene expression at various levels as RNAs; thus, both autophagy and lncRNAs are closely related to the occurrence and development of tumours.

Methods

A total of 232 autophagy-related genes were used to construct a coexpression network to extract autophagy-related lncRNAs. A prognostic signature was constructed by multivariate regression analysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was applied to analyse enrichment in cancer-related pathways. Immune infiltration analysis was used to analyse the relationship between the prognostic signature and the tumour microenvironment.

Results

Nine autophagy-related lncRNAs were used to construct a prognostic model for non-small-cell lung cancer. The median risk score was used to discriminate the high- and low-risk groups, and the low-risk group was found to have better survival. Because KEGG pathway analysis showed that the prognostic signature was enriched in some immune pathways, further analysis of immune infiltration was conducted, and it was found that the prognostic signature did play a unique role in the immune microenvironment. Additionally, the prognostic signature was associated with clinical factors.

Conclusion

We constructed a prognostic model of autophagy-related lncRNAs that can predict the prognosis of non-small-cell lung cancer.

A novel five-gene score to predict complete pathological response to neoadjuvant chemotherapy in ER-positive/HER2-negative breast cancer

Neoadjuvant Chemotherapy (NAC) is not frequently used in ER-positive/HER2-negative breast cancer (BC) because around 10% patients achieve pathological complete response (pCR). Since NAC can result in cancer downstaging both in the breast and axilla and prevent a morbid surgery, thus a score to predict pCR in this population will be crucial to identify patients who can benefit from this approach. A total of 4038 patients from cohorts; GSE25066, GSE20194, Hess, GSE20181, TCGA-BRCA and METBRIC were analyzed. The score was generated by the 5 most highly expressed genes in the Hallmark E2F targets gene set amongst patients in the GSE25066 cohort with ER-positive/HER2-negative BC who achieved pCR. The area under the curve was significantly higher in the score than that for the E2F targets score. High score ER-positive/HER2-negative BCs were significantly associated with higher Nottingham pathological grade, AJCC cancer stage, MKI67 expression levels, intratumor heterogeneity, homologous recombination defects, mutation burden, neoantigen load, and infiltration of anti-cancer immune cells (CD4+, T helper type1, plasmacytoid dendritic cells, M1 macrophages). They also expressed lower abundance of stromal cells including fibroblasts, lymphatic endothelial cells, pericytes and adipocytes consistently in GSE25066, TCGA and METABRIC cohorts. All cell proliferation-related gene sets, G2M checkpoint, E2F targets, MYC targets v1 and v2, Mitotic Spindle, were strongly enriched in high score BCs consistently in 3 cohorts. The gene score was significantly associated with high pCR rate consistently in the GSE25066 (38%, P < 0.001), GSE20194 (16%, P = 0.006), and Hess cohort (23%, P = 0.037). In conclusion, the 5-gene score reflects cancer cell proliferation and immune cell infiltration, and predicts pCR after NAC in ER-positive/HER2-negative breast cancer.

Increased apoptosis is associated with robust immune cell infiltration and cytolytic activity in breast cancer

Tumor infiltrating immune cells plays a critical role in cancer progression. Apoptosis is an autonomous cell death that counteracts tumor growth. To this end, we hypothesized that increased apoptosis in breast cancer is associated with immune cell killing. Apoptosis score of MSigDB Hallmark collection was used to analyze METABRIC cohort (n=1904) and TCGA (n=1069) as validation cohort. High apoptosis tumors enriched cancer promoting signaling pathways; hypoxia, KRAS, TGF-β, PI3K signaling, and was associated with low MKI67 expression and less cell proliferation gene sets, less homologous recombination defects, and less altered fraction. High apoptosis tumors also enriched angiogenesis and high infiltration of vascular endothelial cells, pericytes and stromal cells and significantly enriched inflammation and immune response-related gene sets and high infiltration of CD8, CD4 memory, dendritic cells, M1 and M2 macrophages and significant elevation of cytolytic activity and immune checkpoint molecules, consistently in both cohorts. In conclusion, breast cancer patients with high apoptosis are associated with angiogenesis, immune response, high immune cell infiltration and cytolytic activity. To the best of our knowledge, this is the first study to utilize in silico translational approach to demonstrate the clinical relevance of apoptosis in breast cancer patients in large cohorts.

Low expression of miR-29a is associated with aggressive biology and worse survival in gastric cancer

Advanced gastric cancer (GC) is one of the most lethal cancer types, thus a better understanding of its biology in patients is urgently needed. MicroRNA (miR)-29a is a known tumor suppressive miR that is related to metastasis, but its clinical relevance in GC remains ambiguous. Here, using a large GC patient cohort we hypothesized that low expression of miR-29a in GC is associated with aggressive cancer biology and worse survival. We demonstrated that low miR-29a GC enriched cell proliferation, apoptosis, metastasis, and angiogenesis related gene sets, as well as the higher expression of related genes. Low miR-29a GC was associated with less anti-cancer immune cell infiltration as well as immune related scoring. Low miR-29a GC demonstrated a worse overall survival (OS) as well as disease specific survival (DSS) compared with high expressing miR-29a GC. Notably, low miR-29a expression was the only factor, other than residual tumor status, to be an independent prognostic biomarker of worse OS and DSS. In conclusion, low miR-29a GC was associated with aggressive cancer biology and worse OS as well as DSS. Additionally, low expression of miR-29a was an independent prognostic biomarker of OS and DSS in gastric cancer patients.

G2M checkpoint pathway alone is associated with drug response and survival among cell proliferation-related pathways in pancreatic cancer

Given the severe side effects of the treatments and poor survival, prognostic and predictive biomarkers to guide management of pancreatic cancer are in critical need. We hypothesized that cell proliferation-related pathways are associated with drug response and survival in pancreatic cancer. Six Hallmark cell proliferation-related gene sets (G2M Checkpoint, E2F Targets, MYC Targets V1 and V2, Mitotic Spindle, p53 pathway) defined by MSigDB in gene set variant analysis were evaluated in 3 independent cohorts- TCGA-PAAD (n = 176), GSE57495 (n = 63), and GSE62452 (n = 69). G2M and E2F, as well as Mitotic and p53 pathway correlated highly with other gene sets. All pathways were significantly correlated with MKI67 expression and its proliferation score, but none with cytolytic activity and the rate of pathologically complete resection (R0). All pathways were significantly associated with high alteration and expression of KRAS gene except for MYC v1. G2M, E2F, and p53 pathway were significantly associated with high alteration of TP53 gene. Interestingly, different pathways correlated with the AUC of different cancer therapeutics, such as Gemcitabine (Mitotic: r = 0.706 [P = 0.01]), Paclitaxel (MYC v2: r = -0.636 [P < 0.05]), Apatinib (Mitotic: r = -0.556 [P = 0.03]), Palbociclib (E2F: r = 0.675 [P < 0.01]), and Sorafenib (G2M: r = -0.593 [P = 0.03]). Among all six pathways, only G2M was consistently associated with worse patient survival in all three cohorts. In conclusion, each cell proliferation-related pathway was predictive of a unique agent, and the G2M score alone predicts survival in pancreatic cancer.

Octogenarians' Breast Cancer Is Associated with an Unfavorable Tumor Immune Microenvironment and Worse Disease-Free Survival

Elderly patients are known to have a worse prognosis for breast cancer. This is commonly blamed on their medical comorbidities and access to care. However, in addition to these social issues, we hypothesized that the extreme elderly (octogenarians-patients over 80 years old) have biologically worse cancer with unfavorable tumor immune microenvironment. The Cancer Genomic Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were analyzed. The control (aged 40-65) and octogenarians numbered 668 and 53 in TCGA and 979 and 118 in METABRIC, respectively. Octogenarians had significantly worse breast cancer-specific survival in both cohorts (p < 0.01). Octogenarians had a higher ER-positive subtype rate than controls in both cohorts. Regarding PAM50 classification, luminal-A and -B subtypes were significantly higher in octogenarians, whereas basal and claudin-low subtypes were significantly lower (p < 0.05) in octogenarians. There was no difference in tumor mutation load, intratumor heterogeneity, or cytolytic activity by age. However, the octogenarian cohort was significantly associated with high infiltration of pro-cancer immune cells, M2 macrophage, and regulatory T cells in both cohorts (p < 0.05). Our results demonstrate that octogenarians' breast cancer is associated with worse survival and with an unfavorable tumor immune microenvironment.

Enhanced Thermogenesis in Triple-Negative Breast Cancer Is Associated with Pro-Tumor Immune Microenvironment

Mild cold stress induced by housing mice with a 4T1 triple-negative breast cancer (TNBC) cell implantation model at 22 °C increases tumor growth rate with a pro-tumorigenic immune microenvironment (lower CD8 +T cells, higher myeloid-derived suppressor cells (MDSCs) and regulatory T-cells (Tregs)). Since cold stress also activates thermogenesis, we hypothesized that enhanced thermogenesis is associated with more aggressive cancer biology and unfavorable tumor microenvironment (TME) in TNBC patients. A total of 6479 breast cancer patients from METABRIC, TCGA, GSE96058, GSE20194, and GSE25066 cohorts were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) thermogenesis score. High-thermogenesis TNBC was associated with a trend towards worse survival and with angiogenesis, adipogenesis, and fatty acid metabolism pathways. On the other hand, low-thermogenesis TNBC enriched most of the hallmark cell-proliferation-related gene sets (i.e., mitotic spindle, E2F targets, G2M checkpoint, MYC targets), as well as immune-related gene sets (i.e., IFN-α and IFN-γ response). Favorable cytotoxic T-cell-attracting chemokines CCL5, CXCL9, CXCL10, and CXCL11 were lower; while the MDSC- and Treg-attracting chemokine CXCL12 was higher. There were higher M2 but lower M1 macrophages and Tregs. In conclusion, high-thermogenesis TNBC is associated with pro-tumor immune microenvironment and may serve as biomarker for testing strategies to overcome this immunosuppression.

High BRCA2 Gene Expression is Associated with Aggressive and Highly Proliferative Breast Cancer

BACKGROUND

Mutations of BRCA genes are the most studied in breast cancer, but the clinical relevance of BRCA2 gene expression has been less well studied. Given that BRCA2 is a DNA repair gene, we hypothesized that high BRCA2 expression is associated with highly proliferative and aggressive biology in breast cancer.

MATERIALS AND METHODS

A total of 4342 breast cancer patients were analyzed from The Cancer Genome Atlas (TCGA, n = 1069) as the testing cohort and Gene Expression Omnibus (GEO) dataset GSE96058 (n = 3273) as a validation cohort.

RESULTS

There was no relationship between BRCA2 mutation and BRCA2 gene expression. BRCA2 high expression breast cancer was associated with higher Nottingham grade (p < 0.001), with high proliferation (MKI-67, p < 0.001), and with higher intratumor heterogeneity, homologous recombination deficiency, mutation rate, fraction altered, and neoantigens (all p < 0.001). BRCA2 high expression was associated with E2F1, RB1, PALB2, and PARP, as well as cell proliferation-related gene sets, E2F targets, G2M checkpoints, and mitotic spindle, and with less ESR1 and ER response early and late gene sets. They were associated with significantly reduced number of stroma cells and with high infiltration of both favorable and unfavorable immune cells. BRCA2 high expression significantly correlated with response to olaparib, a PARP inhibitor, and inversely with cyclophosphamide in ER-positive/HER2-negative tumors, but not in TNBC.

CONCLUSIONS

BRCA2 high gene expression was associated with highly proliferative and aggressive breast cancer that was highly immunogenic with better response to PARP inhibitors in ER-positive patients. BRCA2 gene expression may become a candidate marker for aggressive biology and to tailor specific treatment strategies in the future.

Annexin A1 Expression Is Associated with Epithelial-Mesenchymal Transition (EMT), Cell Proliferation, Prognosis, and Drug Response in Pancreatic Cancer

Annexin A1 (ANXA1) is a calcium-dependent phospholipid-binding protein overexpressed in pancreatic cancer (PC). ANXA1 expression has been shown to take part in a wide variety of cancer biology, including carcinogenesis, cell proliferation, invasion, apoptosis, and metastasis, in addition to the initially identified anti-inflammatory effect in experimental settings. We hypothesized that ANXA1 expression is associated with cell proliferation and survival in PC patients. To test this hypothesis, we analyzed 239 PC patients in The Cancer Genome Atlas (TCGA) and GSE57495 cohorts. ANXA1 expression correlated with epithelial-mesenchymal transition (EMT) but weakly with angiogenesis in PC patients. ANXA1-high PC was significantly associated with a high fraction of fibroblasts and keratinocytes in the tumor microenvironment. ANXA1 high PC enriched multiple malignant gene sets, including hypoxia, tumor necrosis factor (TNF)-α signaling via nuclear factor-kappa B (NF-kB), and MTORC1, as well as apoptosis, protein secretion, glycolysis, and the androgen response gene sets consistently in both cohorts. ANXA1 expression was associated with TP53 mutation alone but associated with all KRAS, p53, E2F, and transforming growth factor (TGF)-β signaling pathways and also associated with homologous recombination deficiency in the TCGA cohort. ANXA1 high PC was associated with a high infiltration of T-helper type 2 cells in the TME, with advanced histological grade and MKI67 expression, as well as with a worse prognosis regardless of the grade. ANXA1 expression correlated with a sensitivity to gemcitabine, doxorubicin, and 5-fluorouracil in PC cell lines. In conclusion, ANXA1 expression is associated with EMT, cell proliferation, survival, and the drug response in PC.