A Novel Four-Gene Score to Predict Pathologically Complete (R0) Resection and Survival in Pancreatic Cancer

Predictive biomarkers of pancreatic cancer metastasis: A comprehensive review

This review provides a comprehensive overview of predictive biomarkers associated with metastasis in pancreatic cancer (PC), one of the most aggressive malignancies characterized by late-stage diagnosis and poor prognosis. Metastasis, particularly to the liver, lungs, and lymph nodes, significantly worsens patient outcomes by compromising organ function and promoting disease progression. Reliable biomarkers for predicting and detecting metastasis at early stages are critical for improving survival rates and guiding personalized therapies. This paper highlights both general and specific biomarkers, including genetic mutations, protein expression changes, and carbohydrate tumor markers such as CA19-9. Immunological factors, including PD-L1, inflammatory cytokines, and chemokines, further influence the metastatic process within the tumor microenvironment (TME). Specific biomarkers play pivotal roles in promoting metastasis through mechanisms such as epithelial-to-mesenchymal transition (EMT), tumor microenvironment remodeling, and immune evasion. Emerging markers such as circulating tumor cells (CTCs) and volatile organic compounds (VOCs) offer promising non-invasive tools for metastasis detection and monitoring. This review not only consolidates existing knowledge but also highlights the mechanisms through which specific biomarkers facilitate metastasis. Despite recent progress, challenges such as biomarker standardization, technical variability, and clinical validation remain, and addressing these hurdles is essential for integrating predictive biomarkers into clinical practice. Ultimately, this review contributes to advancing early detection strategies, personalized treatment options, and improved prognosis for PC patients.

Enhanced Error Suppression for Accurate Detection of Low-Frequency Variants

The identification of low-frequency variants remains challenging due to the inevitable high error rates of next-generation sequencing (NGS). Numerous promising strategies employ unique molecular identifiers (UMIs) for error suppression. However, their efficiency depends highly on redundant sequencing and quality control, leading to tremendous read waste and cost inefficiency. Here, we describe a novel approach, enhanced error suppression strategy (EES), that addresses these challenges by (1) optimizing data utilization and reducing read waste by utilizing single-read correction that reserves abundant single reads that complement other single reads or single-strand consensus sequences (SSCSs), and (2) effectively enhancing the accuracy of NGS by employing Bayes' theorem. EES significantly improves variant detection accuracy, achieving a background error rate of less than 4.4 × 10-5 per base pair. Additionally, the data utilization rate is dramatically increased, with a 22.9-fold enhancement in duplex consensus sequence (DCS) recovery compared to traditional methodologies. Furthermore, EES demonstrates superior error suppression performance across various base substitutions. In conclusion, EES represents a significant advancement in detecting low-frequency variants by improving data utilization and reducing sequencing errors. It potentially enhances the sensitivity and accuracy of NGS applications, proving highly valuable in clinical and research contexts where precise variant detection is critical.

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.

Integrated analysis of single-cell sequencing and weighted co-expression network identifies a novel signature based on cellular senescence-related genes to predict prognosis in glioblastoma

BACKGROUND

Glioblastoma (GBM) is a highly aggressive cancer with heavy mortality rates and poor prognosis. Cellular senescence exerts a pivotal influence on the development and progression of various cancers. However, the underlying effect of cellular senescence on the outcomes of patients with GBM remains to be elucidated.

METHODS

Transcriptome RNA sequencing data with clinical information and single-cell sequencing data of GBM cases were obtained from CGGA, TCGA, and GEO (GSE84465) databases respectively. Single-sample gene set enrichment analysis (ssGSEA) analysis was utilized to calculate the cellular senescence score. WGCNA analysis was employed to ascertain the key gene modules and identify differentially expressed genes (DEGs) associated with the cellular senescence score in GBM. The prognostic senescence-related risk model was developed by least absolute shrinkage and selection operator (LASSO) regression analyses. The immune infiltration level was calculated by microenvironment cell populations counter (MCPcounter), ssGSEA, and xCell algorithms. Potential anti-cancer small molecular compounds of GBM were estimated by "oncoPredict" R package.

RESULTS

A total of 150 DEGs were selected from the pink module through WGCNA analysis. The risk-scoring model was constructed based on 5 cell senescence-associated genes (CCDC151, DRC1, C2orf73, CCDC13, and WDR63). Patients in low-risk group had a better prognostic value compared to those in high-risk group. The nomogram exhibited excellent predictive performance in assessing the survival outcomes of patients with GBM. Top 30 potential anti-cancer small molecular compounds with higher drug sensitivity scores were predicted.

CONCLUSION

Cellular senescence-related genes and clusters in GBM have the potential to provide valuable insights in prognosis and guide clinical decisions.

Gastric cancer with enhanced myogenesis is associated with less cell proliferation, enriched epithelial-to-mesenchymal transition and angiogenesis, and poor clinical outcomes

Gastric cancer (GC) remains a lethal disease, with over 26,000 new cases and more than 11,000 deaths annually in the US. Thus, a deeper understanding of GC biology is critical to improve survival. Myogenesis is the formation of muscle fibers, which is a mesodermal tissue. In cancer, epithelial-to-mesenchymal transition (EMT) is a known phenomenon that promotes metastasis and poor survival. Given that myogenesis produces mesenchymal cells, we hypothesized that GC with increased myogenesis is linked to aggressive tumor behaviors and less favorable outcomes. In this study, three GC patient cohorts: TCGA (n=375), GSE26253 (n=432), and GSE84437 (n=482), were analyzed. The "MYOGENESIS" set in the Hallmark collection which comprises 200 myogenesis-related genes was analyzed to perform gene set variation analysis to create a score to quantify the myogenesis activity. Our results showed that T category of AJCC cancer staging that reflects the tumor invasion to stomach wall consistently correlated with myogenesis activity in two GC cohorts. High myogenesis GC was associated with lower cell proliferation, evidenced by reduced proliferation scores, decreased Ki67 gene expression, and less enrichment of E2F Targets, G2M checkpoint, MYC Targets V1, and V2 gene sets. High myogenesis tumors showed increased stromal cells (fibroblasts and adipocytes) infiltration within the tumor microenvironment, as well as less silent and non-silent mutation rates and copy number alterations. Higher lymphocyte infiltration, leukocyte fraction, T-cell receptor richness, and B-cell receptor richness were associated with high myogenesis GC. However, infiltration of CD4 cells, T helper type 1 and 2 cells, Natural Killer cells, regulatory T cells, and plasma cells was lower, with increased infiltration of dendritic cells in high myogenesis GC. High myogenesis GC enriched EMT, Hedgehog, TGF-β, and KRAS gene sets. Furthermore, it was associated with enhanced angiogenesis, evidenced by enrichment of Angiogenesis, Coagulation, and Hypoxia gene sets, and increased infiltration of microvascular and lymphatic endothelial cells and pericytes. High myogenesis GC consistently correlated with worse overall survival in all three cohorts, and worse disease-specific and progression-free survival in the TCGA cohort. Hence, our findings suggest that GC with enhanced myogenesis is associated with decreased cell proliferation, increased EMT and angiogenesis, and worse prognosis.

Robinin inhibits pancreatic cancer cell proliferation, EMT and inflammation via regulating TLR2-PI3k-AKT signaling pathway

PURPOSE

To investigate the anti-tumor effect of Robinin (Toll-like receptor 2 inhibitor) in pancreatic cancer cells via regulating tumor microenvironment.

METHODS

The effects of Robinin on cell proliferation or migration in Mia-PACA2 and PANC-1 were determined, using CCK8 or wound healing assay, respectively. The typical markers of EMT (αSMA and snail) and the inflammation markers (IL-6 and TNF-α) were all detected by western blot. CU-T12-9 (TLR2 agonist) was used to rescue Robinin's effect. PI3k-p85α and Phosphorylated-AKT (p-AKT) were evaluated, compared to the β-actin and AKT, using western blot.

RESULTS

Robinin significantly inhibited cell proliferation and migration in Mia-PACA2 and PANC-1, compared to HPNE (**P < 0.01). Robinin also attenuated the expression of α-SMA and snail in Mia-PACA2, and PANC-1 (**P < 0.01). Besides, it was found that expression of IL-6 and TNF-α were diminished in presence of Robinin in Mia-PACA2, and PANC-1 (**P < 0.01). Western blot confirmed that Robinin could target on TLR2, and further downregulated PI3k-AKT signaling pathway to exert biological function.

CONCLUSIONS

Robinin exerts anti-tumor effect perhaps via downregulating inflammation and EMT in pancreatic cancer cell through inhibiting TLR2-PI3k-AKT signaling pathway. Robinin may be a novel agent in adjuvant therapy of pancreatic cancer.

Pancreatic ductal adenocarcinoma with a high expression of alcohol dehydrogenase 1B is associated with less aggressive features and a favorable prognosis

Alcohol dehydrogenase (ADH) oxidizes alcohol into acetaldehyde (AA), which is a known carcinogen. Aldehyde dehydrogenase (ALDH) oxidizes AA into acetate. Therefore, pancreatic cancer that expresses a high level of ADH1B that generates more AA is expected to be associated with aggressive cancer. On the other hand, given that the differentiated cells that retain their cellular functions typically exhibit lower proliferation rates, it remains unclear whether pancreatic adenocarcinoma (PDAC) with high ADH1B gene expression is linked to aggressive features in patients. The Cancer Genome Atlas (n = 145) was used to obtain data of PDAC patients and GSE62452 cohort (n = 69) was used as a validation cohort. PDAC with high ADH1B expression was associated with less cancer cell proliferation as evidenced by lower MKI67 expression and lower histological grade; with a higher fraction of stromal cells consistent with less proliferative cancer. PDAC with high ADH1B expression also had lower homologous recombination deficiency and mutation rates, lower KRAS and TP53 mutation rates. ADH1B expression correlated with ALDH2 expression in PDAC, but not with DNA repair genes. High ADH1B expression PDAC was associated with high infiltration of anti-cancerous CD8+ T cells and pro-cancerous M2 macrophages but with lower levels of Th1 T cells, with a higher cytolytic activity. PDAC patients with a high ADH1B expression had better disease-specific survival (DSS) and overall survival (OS) and ADH1B was an independent prognostic biomarker for both DSS (HR = 0.89, 95% CI = 0.80-0.99, P = 0.045) and OS (HR = 0.90, 95% CI = 0.82-0.99, P = 0.044) in multivariate analysis. In conclusion, PDAC with high ADH1B expression had less cell proliferation and malignant features, along with higher immune cell infiltration, and had a better 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.

GALNT1 Expression Is Associated with Angiogenesis and Is a Prognostic Biomarker for Breast Cancer in Adolescents and Young Adults (AYA)

It is well established that genetic information differs amongst the adolescent and young adult population (AYA) and older patients. Although several studies on genetic information have been conducted, no current prognostic biomarker exists to help differentiate survival outcomes amongst AYA patients. The GALNT family of genes have been associated with several cancer etiologies, such as the Tn antigen and epithelial-mesenchymal transition (EMT); however, the clinical significance of GALNT1 expression in breast cancer (BC) remains unclear. We investigated the clinical relevance of GALNT1 expression in BC using two large independent cohorts. We found that, although triple-negative BC (TNBC) had the highest GALNT1 expression compared to ER-positive/HER2-negative BC, GALNT1 levels in BC were not associated with clinical aggressiveness, including histological grade, AJCC stage and N-category, and patient survival, consistently in both the METABRIC and GSE96058 cohorts. There was also no biological difference between low- and high-GALNT1 expression BC, as analyzed by hallmark gene sets via gene set enrichment analysis (GSEA). Further, no significant difference was found in GALNT1 expression levels among AYAs and older patients. However, high GALNT1 expression was associated with significantly worse survival in AYA patients, in both cohorts. Furthermore, high GALNT1 expression was found to be an independent factor among several clinical features, including subtype, histological grade, AJCC T and N-category, in AYA patients. In both cohorts, BC with high GALNT1 expression demonstrated low levels of CD8+ T-cell infiltration, but not other anti-cancerous or pro-cancerous immune cells. Finally, high levels of GALNT1 BC demonstrated increased EMT, angiogenesis, and protein secretion in the AYA population, but not in older patients. In conclusion, our findings demonstrate that GALNT1 expression was found to be associated with angiogenesis and EMT, and may have potential as prognostic biomarker, specifically in AYA patients.

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.

High miR-99b expression is associated with cell proliferation and worse patient outcomes in breast cancer

Although miR-99b is a known suppressive microRNA (miRNA) in several cancers, its role in breast cancer has not been elucidated. In this study, we examined the clinical relevance of miR-99b expression in breast cancer. We analyzed miRNA and mRNA expression and their relationships with clinical parameters in 1,961 breast cancer samples from two independent large cohorts, the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA). Several algorithms, including gene set enrichment analysis (GSEA) and xCell, have been used to investigate biological functions and the tumor microenvironment. High miR-99b expression significantly enriched the mTORC1 signaling gene set in breast cancer (NES = 1.63, FDR = 0.03, and NES = 1.58, FDR = 0.10, in METABRIC and TCGA, respectively). No other mechanisms, including the epithelial mesenchymal transition, NFκB, and TGF-β signaling, were consistently enriched in both cohorts. MiR-99b-high breast cancer was associated with high homologous recombination deficiencies, intratumor heterogeneity, and high rates of mutation and neoantigens. In agreement, miR-99b-high breast cancer was associated with increased cell proliferation, correlating with Nottingham histological grade, and significant enrichment of E2F targets, G2/M checkpoint, and mitotic spindle gene sets consistently in both cohorts (P = 0.01, P < 0.001). High miR-99b levels were also associated with low stromal cell fractions in the tumor microenvironment, including adipocytes, keratinocytes, and lymphatic endothelial cells (P < 0.001). However, in both cohorts, miR-99b expression was not associated with significant infiltration of immune cells, except dendritic cells (P = 0.006, 0.020). Finally, in both cohorts, breast cancer with high miR-99b expression was significantly associated with worse disease-free survival (DSS) and overall survival (OS), particularly in estrogen receptor (ER)-positive/human epidermal growth factor (HER)2-negative breast cancer (DSS hazard ratio (HR) 1.29, 95% confidence interval (CI) 1.10-1.51, P < 0.001 in the METABRIC cohort and HR 1.82, 95% CI 1.12-2.98, P = 0.017 in the TCGA cohort). In conclusion, breast cancer with high miR-99b expression was significantly associated with mTORC1 signaling, cell proliferation, and decreased patient survival, particularly in the ER-positive/HER2-negative subtype.

Low adipocyte hepatocellular carcinoma is associated with aggressive cancer biology and with worse survival

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide, and non-alcoholic fatty liver disease is strongly associated with its development. To explore the role of adipocytes in HCC, we investigated intratumoral adipocytes, also known as cancer-associated adipocytes (CAA). Based on our prior breast cancer findings, we hypothesized that low intratumoral adipocytes would be associated with aggressive cancer biology, worse tumor microenvironment (TME), and clinical outcomes. The Cancer Genome Atlas (TCGA) was used and validated by the Gene Expression Omnibus (GEO) cohort. xCell algorithm was used to quantify intratumoral adipocytes and top 90% were defined as adipocyte high (AH) and bottom 10% as adipocyte low (AL). We found that AL-HCC was significantly associated with worse disease-free survival (DFS), disease-specific survival (DSS), and overall survival (OS). AL-HCC were higher-grade, had high MKI67 expression, enriched cell proliferation-related gene sets, and had increased altered fraction, aneuploidy, and homologous recombination defects. Also, anti-cancer immune cells, CD8, Th1, and M1 cells, as well as pro-cancer Th2 cells were increased in AL-HCC. Micro-RNAs miR-122 (associated with cholesterol metabolism) and miR-885 (associated with liver pathologies) were significantly increased in the AL TME. In conclusion, we found that AL-HCC has worse patient outcomes and is biologically more aggressive with enhanced cell proliferation. Our findings take initial steps to clarify the role of adipocytes in HCC.

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.

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.

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.

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.

Identification of Novel Metabolism-Associated Subtypes for Pancreatic Cancer to Establish an Eighteen-Gene Risk Prediction Model

Pancreatic cancer (PanC) is an intractable malignancy with a high mortality. Metabolic processes contribute to cancer progression and therapeutic responses, and histopathological subtypes are insufficient for determining prognosis and treatment strategies. In this study, PanC subtypes based on metabolism-related genes were identified and further utilized to construct a prognostic model. Using a cohort of 171 patients from The Cancer Genome Atlas (TCGA) database, transcriptome data, simple nucleotide variants (SNV), and clinical information were analyzed. We divided patients with PanC into metabolic gene-enriched and metabolic gene-desert subtypes. The metabolic gene-enriched subgroup is a high-risk subtype with worse outcomes and a higher frequency of SNVs, especially in KRAS. After further characterizing the subtypes, we constructed a risk score algorithm involving multiple genes (i.e., NEU2, GMPS, PRIM2, PNPT1, LDHA, INPP4B, DPYD, PYGL, CA12, DHRS9, SULT1E1, ENPP2, PDE1C, TPH1, CHST12, POLR3GL, DNMT3A, and PGS1). We verified the reproducibility and reliability of the risk score using three validation cohorts (i.e., independent datasets from TCGA, Gene Expression Omnibus, and Ensemble databases). Finally, drug prediction was completed using a ridge regression model, yielding nine candidate drugs for high-risk patients. These findings support the classification of PanC into two metabolic subtypes and further suggest that the metabolic gene-enriched subgroup is associated with worse outcomes. The newly established risk model for prognosis and therapeutic responses may improve outcomes in patients with PanC.

Liquid Biopsy, ctDNA Diagnosis through NGS

Liquid biopsy with circulating tumor DNA (ctDNA) profiling by next-generation sequencing holds great promise to revolutionize clinical oncology. It relies on the basis that ctDNA represents the real-time status of the tumor genome which contains information of genetic alterations. Compared to tissue biopsy, liquid biopsy possesses great advantages such as a less demanding procedure, minimal invasion, ease of frequent sampling, and less sampling bias. Next-generation sequencing (NGS) methods have come to a point that both the cost and performance are suitable for clinical diagnosis. Thus, profiling ctDNA by NGS technologies is becoming more and more popular since it can be applied in the whole process of cancer diagnosis and management. Further developments of liquid biopsy ctDNA testing will be beneficial for cancer patients, paving the way for precision medicine. In conclusion, profiling ctDNA with NGS for cancer diagnosis is both biologically sound and technically convenient.

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.

Identification of Prognostic miRNAs Associated With Immune Cell Tumor Infiltration Predictive of Clinical Outcomes in Patients With Non-Small Cell Lung Cancer

Background

A detailed means of prognostic stratification in patients with non-small cell lung cancer (NSCLC) is urgently needed to support individualized treatment plans. Recently, microRNAs (miRNAs) have been used as biomarkers due to their previously reported prognostic roles in cancer. This study aimed to construct an immune-related miRNA signature that effectively predicts NSCLC patient prognosis.

Methods

The miRNAs and mRNA expression and mutation data of NSCLC was obtained from The Cancer Genome Atlas (TCGA). Immune-associated miRNAs were identified using immune scores calculated by the ESTIMATE algorithm. LASSO-penalized multivariate survival models were using for development of a tumor immune-related miRNA signature (TIM-Sig), which was evaluated in several public cohorts from the Gene Expression Omnibus (GEO) and the CellMiner database. The miRTarBase was used for constructing the miRNA-target interactions.

Results

The TIM-Sig, including 10 immune-related miRNAs, was constructed and successfully predicted overall survival (OS) in the validation cohorts. TIM-Sig score negatively correlated with CD8+ T cell infiltration, IFN-γ expression, CYT activity, and tumor mutation burden. The correlation between TIM-Sig score and genomic mutation and cancer chemotherapeutics was also evaluated. A miRNA-target network of 10 miRNAs in TIM-Sig was constructed. Further analysis revealed that these target genes showed prognostic value in both lung squamous cell carcinoma and adenocarcinoma.

Conclusions

We concluded that the immune-related miRNAs demonstrated a potential value in clinical prognosis.

Adipogenesis in triple-negative breast cancer is associated with unfavorable tumor immune microenvironment and with worse survival

Cancer-associated adipocytes are known to cause inflammation; however, the role of adipogenesis, the formation of adipocytes, in breast cancer is unclear. We hypothesized that intra-tumoral adipogenesis reflects a different cancer biology than abundance of intra-tumoral adipocytes. The Molecular Signatures Database Hallmark adipogenesis gene set of gene set variant analysis was used to quantify adipogenesis. Total of 5,098 breast cancer patients in multiple cohorts (training; GSE96058 (n = 3273), validation; TCGA (n = 1069), treatment response; GSE25066 (n = 508) and GSE20194 (n = 248)) were analyzed. Adipogenesis did not correlate with abundance of adipocytes. Adipogenesis was significantly lower in triple negative breast cancer (TNBC). Elevated adipogenesis was significantly associated with worse survival in TNBC, but not in the other subtypes. High adipogenesis TNBC was significantly associated with low homologous recombination deficiency, but not with mutation load. High adipogenesis TNBC enriched metabolism-related gene sets, but neither of cell proliferation- nor inflammation-related gene sets, which were enriched to adipocytes. High adipogenesis TNBC was infiltrated with low CD8⁺ T cells and high M2 macrophages. Although adipogenesis was not associated with neoadjuvant chemotherapy response, high adipogenesis TNBC was significantly associated with low expression of PD-L1 and PD-L2 genes, and immune checkpoint molecules index. In conclusion, adipogenesis in TNBC was associated with cancer metabolism and unfavorable tumor immune microenvironment, which is different from abundance of adipocytes.

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.

Development and Validation of a 7-Gene Prognostic Signature to Improve Survival Prediction in Pancreatic Ductal Adenocarcinoma

Background: Previous prognostic signatures of pancreatic ductal adenocarcinoma (PDAC) are mainly constructed to predict the overall survival (OS), and their predictive accuracy needs to be improved. Gene signatures that efficaciously predict both OS and disease-free survival (DFS) are of great clinical significance but are rarely reported. Methods: Univariate Cox regression analysis was adopted to screen common genes that were significantly associated with both OS and DFS in three independent cohorts. Multivariate Cox regression analysis was subsequently performed on the identified genes to determine an optimal gene signature in the MTAB-6134 training cohort. The Kaplan-Meier (K-M), calibration, and receiver operating characteristic (ROC) curves were employed to assess the predictive accuracy. Biological process and pathway enrichment analyses were conducted to elucidate the biological role of this signature. Results: Multivariate Cox regression analysis determined a 7-gene signature that contained ASPH, DDX10, NR0B2, BLOC1S3, FAM83A, SLAMF6, and PPM1H. The signature had the ability to stratify PDAC patients with different OS and DFS, both in the training and validation cohorts. ROC curves confirmed the moderate predictive accuracy of this signature. Mechanically, the signature was related to multiple cancer-related pathways. Conclusion: A novel OS and DFS prediction model was constructed in PDAC with multi-cohort and cross-platform compatibility. This signature might foster individualized therapy and appropriate management of PDAC patients.

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.