Lymphovascular invasion in breast cancer is associated with gene expression signatures of cell proliferation but not lymphangiogenesis or immune response

TCGA-Reports: A machine-readable pathology report resource for benchmarking text-based AI models

In cancer research, pathology report text is a largely untapped data source. Pathology reports are routinely generated, more nuanced than structured data, and contain added insight from pathologists. However, there are no publicly available datasets for benchmarking report-based models. Two recent advances suggest the urgent need for a benchmark dataset. First, improved optical character recognition (OCR) techniques will make it possible to access older pathology reports in an automated way, increasing the data available for analysis. Second, recent improvements in natural language processing (NLP) techniques using artificial intelligence (AI) allow more accurate prediction of clinical targets from text. We apply state-of-the-art OCR and customized post-processing to report PDFs from The Cancer Genome Atlas, generating a machine-readable corpus of 9,523 reports. Finally, we perform a proof-of-principle cancer-type classification across 32 tissues, achieving 0.992 average AU-ROC. This dataset will be useful to researchers across specialties, including research clinicians, clinical trial investigators, and clinical NLP researchers.

Profiling of Lymphovascular Space Invasion in Cervical Cancer Revealed PI3K/Akt Signaling Pathway Overactivation and Heterogenic Tumor-Immune Microenvironments

Lymphovascular space invasion (LVSI) is the presence of tumor emboli in the endothelial-lined space at the tumor body's invasive edge. LVSI is one of three Sedlis criteria components-a prognostic tool for early cervical cancer (CC)-essential for indicating poor prognosis, such as lymph node metastasis, distant metastasis, or shorter survival rate. Despite its clinical significance, an in-depth comprehension of the molecular mechanisms or immune dynamics underlying LVSI in CC remains elusive. Therefore, this study investigated tumor-immune microenvironment (TIME) dynamics of the LVSI-positive group in CC. RNA sequencing included formalin-fixed paraffin-embedded (FFPE) slides from 21 CC patients, and differentially expressed genes (DEGs) were analyzed. Functional analysis and immune deconvolution revealed aberrantly enriched PI3K/Akt pathway activation and a heterogenic immune composition with a low abundance of regulatory T cells (Treg) between LVSI-positive and LVSI-absent groups. These findings improve the comprehension of LSVI TIME and immune mechanisms, benefiting targeted LVSI therapy for CC.

Development of a Machine Learning-Based Prognostic Model for Hormone Receptor-Positive Breast Cancer Using Nine-Gene Expression Signature

Background

Determining the prognosis of hormone receptor positive (HR+) breast cancer (BC), which accounts for 80% of all BCs, is critical in improving survival outcomes. Stratifying individuals at high risk of BC-related mortality and improving prognosis has been the focus of research for over a decade. However, these tools are not universal as they are limited to clinical factors. We hypothesized that a new framework for predicting prognosis in HR+ BC patients can develop using artificial intelligence.

Methods

A total of 2,338 HR+ human epidermal growth factor receptor 2 negative (HER2-) BC cases were analyzed from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), The Cancer Genome Atlas (TCGA), and Gene Expression Omnibus (GEO) cohorts. Groups were then divided into high- and low-risk categories utilizing a recurrence prediction model (RPM). An RPM was created by extracting nine prognosis-related genes from over 18,000 genes using a logistic progression model.

Results

Risk classification by RPM was significantly stratified in both the discovery cohort and validation cohort. In the time-dependent area under the curve analysis, there was some variation depending on the cohort, but accuracy was found to decline significantly after about 10 years. Cell cycle related gene sets, MYC, and PI3K-AKT-mTOR signaling were enriched in high-risk tumors by the Gene Set Enrichment Analysis. High-risk tumors were associated with high levels of immune cells from the lymphoid and myeloid lineage and immune cytolytic activity, as well as low levels of stem cells and stromal cells. High-risk tumors were also associated with poor therapeutic effects of chemotherapy and endocrine therapy.

Conclusions

This model was able to stratify prognosis in multiple cohorts. This is because the model reflects major BC therapeutic target pathways and tumor immune microenvironment and, further is supported by the therapeutic effect of chemotherapy and endocrine therapy.

Benchmark Pathology Report Text Corpus with Cancer Type Classification

In cancer research, pathology report text is a largely un-tapped data source. Pathology reports are routinely generated, more nuanced than structured data, and contain added insight from pathologists. However, there are no publicly-available datasets for benchmarking report-based models. Two recent advances suggest the urgent need for a benchmark dataset. First, improved optical character recognition (OCR) techniques will make it possible to access older pathology reports in an automated way, increasing data available for analysis. Second, recent improvements in natural language processing (NLP) techniques using AI allow more accurate prediction of clinical targets from text. We apply state-of-the-art OCR and customized post-processing to publicly available report PDFs from The Cancer Genome Atlas, generating a machine-readable corpus of 9,523 reports. We perform a proof-of-principle cancer-type classification across 32 tissues, achieving 0.992 average AU-ROC. This dataset will be useful to researchers across specialties, including research clinicians, clinical trial investigators, and clinical NLP researchers.

Correlation Analysis of Pathological Features and Axillary Lymph Node Metastasis in Patients with Invasive Breast Cancer

Objective

To investigate the risk factors of axillary lymph node metastasis in patients with invasive breast cancer.

Methods

This study retrospectively included 122 cases of invasive breast cancer patients admitted to the First Medical Center of PLA General Hospital from January 2019 to September 2020. According to postoperative pathological results, axillary lymph node metastasis was divided into axillary lymph node metastasis (ALNM) group (n =40) and non-axillary lymph node metastasis (NALNM) group (n =82). General demographic information was collected and compared between the two groups. Collected pathological results included lymphovascular invasion (LVI) and the expression of estrogen receptor (ER), progestogen receptor (PR), human epidermal growth factor receptor 2 (HER-2), and Ki-67 detected by immunohistochemistry. Imaging parameters of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) including apparent diffusion coefficient (ADC), early enhanced rate, and time-intensity curve (TIC) were also included into univariate analysis. The variables with differences between the two groups were compared by univariate analysis, and the related factors of axillary lymph node metastasis were analyzed by logistic regression model.

Results

There was no significant difference in general demographic information between the two groups. No significant differences were found in the positive rates of HER-2, ER, PR, Ki-67, pathological types, and clavicular lymph node metastasis and skin chest wall invasion between the two groups (P > 0.05). The proportion of LVI in ALNM group was significantly higher than that in NALNM group (37.50% vs. 6.10%, P < 0.001). The proportion of breast cancer on the left side in the ALNM group was higher than that in the NALNM group, and the difference was statistically significant (70.00% vs. 47.56%, P = 0.019). There were no significant differences in the imaging parameters obtained by DCE-MRI between the two groups. Binary logistics regression analysis showed that LVI (OR =12.258, 95% CI =3.681-40.812, P < 0.001) and left breast cancer (OR =3.598, 95% CI =1.404-9.219, P = 0.008) were risk factors for axillary lymph node metastasis in patients with invasive breast cancer.

Conclusion

The formation of vascular tumor thrombi in breast cancer tissue and left breast cancer are risk factors for axillary lymph node metastasis in invasive breast cancer and might be helpful for preoperative detailed assessment of the patient's condition.

The ambivalent nature of the relationship between lymphatics and cancer

Do lymphatic vessels support cancer cells? Or are they vessels that help suppress cancer development? It is known that the lymphatic system is a vehicle for tumor metastasis and that the lymphangiogenic regulator VEGF-C supports the tumor. One such role of VEGF-C is the suppression of the immune response to cancer. The lymphatic system has also been correlated with an increase in interstitial fluid pressure of the tumor microenvironment. On the other hand, lymphatic vessels facilitate immune surveillance to mount an immune response against tumors with the support of VEGF-C. Furthermore, the activation of lymphatic fluid drainage may prove to filter and decrease tumor interstitial fluid pressure. In this review, we provide an overview of the dynamic between lymphatics, cancer, and tumor fluid pressure to suggest that lymphatic vessels may be used as an antitumor therapy due to their capabilities of immune surveillance and fluid pressure drainage. The application of this potential may help to prevent tumor proliferation or increase the efficacy of drugs that target cancer.

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.

NR2F1, a Tumor Dormancy Marker, Is Expressed Predominantly in Cancer-Associated Fibroblasts and Is Associated with Suppressed Breast Cancer Cell Proliferation

BACKGROUND

Tumor dormancy is a crucial mechanism responsible for the late recurrence of breast cancer. Thus, we investigated the clinical relevance of the expression of NR2F1, a known dormancy biomarker.

METHODS

A total of 6758 transcriptomes of bulk tumors from multiple breast cancer patient cohorts and two single-cell sequence cohorts were analyzed.

RESULTS

Breast cancer (BC) with high NR2F1 expression enriched TGFβ signaling, multiple metastases, and stem cell-related pathways. Cell proliferation-related gene sets were suppressed, and MKi67 expression was lower in high NR2F1 BC. In tumors with high Nottingham grade, NR2F1 expression was found to be lower. There was no consistent relationship between NR2F1 expression and metastasis or survival. Cancer mutation rates, immune responses, and immune cell infiltrations were lower in high NR2F1 tumors, whereas the infiltration of stromal cells including cancer-associated fibroblasts (CAFs) was higher. NR2F1 was predominantly expressed in CAFs, particularly inflammatory CAFs, rather than in cancer cells, consistently in the two single-cell sequence cohorts.

CONCLUSIONS

NR2F1 expression in breast cancer is associated with tumor dormancy traits, and it is predominantly expressed in CAFs in the tumor microenvironment.

Clinical Relevance of Estrogen Reactivity in the Breast Cancer Microenvironment

Purpose

Estrogen signals play an important role in the phenotype of estrogen receptor-positive breast cancer. However, comprehensive analyses of the effect of responsiveness to estrogen signals on the tumor microenvironment and survival in large cohorts of primary breast cancer patients have been lacking. We aimed to test the hypothesis that estrogen reactivity affects gene expression and immune cell infiltration profiles in the tumor microenvironment and survival.

Methods

A total of 3,098 breast cancer cases were analyzed: 1,904 from the Molecular Taxonomy of Breast Cancer (METABRIC) cohort, 1,082 from The Cancer Genome Atlas (TCGA) cohort, and 112 from the Hokkaido University Hospital cohort. We divided the group into estrogen reactivity-high and estrogen reactivity-low groups utilizing the scores of ESTROGEN_RESPONSE_EARLY and ESTROGEN_RESPONSE_LATE in Gene Set Variation Analysis.

Results

Breast cancer with high estrogen reactivity was related to Myc targets, metabolism-related signaling, cell stress response, TGF-beta signaling, androgen response, and MTORC1 signaling gene sets in the tumor microenvironment. Low estrogen reactivity was related to immune-related proteins, IL2-STAT5 signaling, IL6-JAK-STAT3 signaling, KRAS signaling, cell cycle-related gene sets, and EMT. In addition, breast cancer with high levels of estrogen reactivity had low immune cytolytic activity and low levels of immunostimulatory cells. It also had low levels of stimulatory and inhibitory factors of the cancer immunity cycle. Patients with high estrogen reactivity were also associated with a better prognosis.

Conclusion

We demonstrated the relationship between estrogen reactivity and the profiles of immune cells and gene expression, as well as survival.

An evaluation of lymphovascular invasion in relation to biology and prognosis according to subtypes in invasive breast cancer

Lymphovascular invasion (LVI) is associated with a poor outcome in breast cancer. The purpose of the present study was to evaluate the clinical significance of LVI in primary breast cancer and to investigate disease-free survival as a prognostic marker according to the breast cancer subtypes. This study examined 4,652 consecutive cases of invasive breast cancer excluding the patients with non-invasive cancer, stage IV and those who underwent neo-adjuvant therapy from February 2002 to February 2021. The clinicopathological characteristics and prognosis of LVI-positive and -negative tumors were compared. LVI was evaluated in H&E staining specimens from surgically resected samples. The LVI expression rates were 29.2% (low, 19.7%; high, 9.5%) in all primary cases. The LVI-positive rate was significantly associated with specimens with the following characteristics: ER/PgR-negative, HER2-positive, p53 overexpression, higher Ki-67 index values, higher nuclear grade, positive nodes and larger tumors. Moreover, the subtypes were significantly associated with LVI positivity; 20% in Luminal A, 34.6% in Luminal B, 40.9% in Lumina/HER2, 38.1% in HER2-enriched and 29.8% in triple negative (TN). There were significant differences in disease-free survival between LVI status in Luminal A, Luminal B and TN subtypes, but there was no difference in the Luminal/HER2 and HER2-enriched subtypes. A multivariate analysis revealed that LVI was a significant factor in Luminal B and TN subtypes. Overall, LVI was significantly associated with the advanced and aggressive characteristics in breast cancer. Luminal A type had a lower LVI rate, and HER2 type had a higher LVI rate. Moreover, LVI was a significant prognostic factor in Luminal B and TN subtypes. These data suggested that the LVI status was useful in predicting the prognosis in HER2 negative breast cancer cases.

Gastric cancer with enhanced apical junction pathway has increased metastatic potential and worse clinical outcomes

Excessive intercellular connection at confluency may be limiting further cell growth or a sign of aggressive biology in the cell culture. As apical junction complex is a main component of cell-to-cell connection, we aimed to investigate gastric cancer biology using Apical Junction Pathway score that we generated using Gene set variant analysis (GSVA) of the "Hallmark Apical Junction" gene set. 1,239 gastric cancer patients from the Cancer Genome Atlas (TCGA) and two GSE cohorts were included in this study. The cohorts were dichotomized using the median of the score. Apical Junction Pathway score high gastric cancer was not consistently associated with increased cell proliferation or immune cell infiltration. On the other hand, Apical Junction Pathway score high gastric cancer was associated with significantly higher infiltration of stromal cells, such as endothelial cells; hence, increased neovascularization and angiogenesis in the tumor microenvironment (TME) were speculated. Gene set enrichment analysis (GSEA) confirmed increased expression of epithelial mesenchymal transition (EMT) and angiogenesis in the high Apical Junction Pathway score group (false discovery rate (FDR) <0.25). Lastly, the high Apical Junction Pathway score group was associated with more aggressive clinicopathological characteristics, such as significantly higher American Joint Committee on Cancer (AJCC) T-category and higher pathological stage, leading to worse disease-specific survival and overall survival (P<0.05, respectively). In conclusion, enhanced Apical Junction Pathway score gastric cancer was associated with aggressive clinical characteristics leading to shorter survival likely due to increased metastatic potential from EMT and angiogenesis.

Feasibility of Shear Wave Elastography Imaging for Evaluating the Biological Behavior of Breast Cancer

Objective

To explore the feasibility of shear wave elastography (SWE) parameters for assessing the biological behavior of breast cancer.

Materials and Methods

In this prospective study, 224 breast cancer lesions in 216 female patients were examined by B-mode ultrasound and shear wave elastography in sequence. The maximum size (S_max) of the lesion was measured by B-mode ultrasound, and then shear wave elastography was performed on this section to obtain relevant parameters, including maximum elasticity (E_max), mean elasticity (E_mean), standard deviation of elasticity (SD), and the area ratio of shear wave elastography to B-mode ultrasound (AR). The relationship between SWE parameters and pathological type, histopathological classification, histological grade, lymphovascular invasion status (LVI), axillary lymph node status (ALN), and immunohistochemistry of breast cancer lesions was performed according to postoperative pathology.

Results

In the univariate analysis, the pathological type and histopathological classification of breast cancer were not significantly associated with SWE parameters; with an increase in the histological grade of invasive ductal carcinoma (IDC), SD (p = 0.016) and S_max (p = 0.000) values increased. In the ALN-positive group, S_max (p = 0.004) was significantly greater than in the ALN-negative group; S_max (p = 0.003), E_max (p = 0.034), and SD (p = 0.045) were significantly higher in the LVI-positive group than in the LVI-negative group; SD (p = 0.043, p = 0.047) and S_max (p = 0.000, p = 0.000) were significantly lower in the ER⁺ and PR⁺ groups than in the ER^- and PR^- groups, respectively; AR (p = 0.032) was significantly higher in the ER⁺ groups than in the ER^- groups, and S_max (p = 0.002) of the HER2⁺ group showed higher values than that of the HER2^- group; S_max (p = 0.000), SD (p = 0.006), and E_max (p = 0.004) of the Ki-67 high-expression group showed significantly higher values than those of the Ki-67 low-expression group. In the multivariate analysis, Ki-67 was an independent factor of S_max (p = 0.005), E_max (p = 0.004), and SD (p = 0.006); ER was an independent influencing factor of S_max (p = 0.000) and AR (p = 0.032). LVI independently influences S_max (p = 0.006).

Conclusions

The SWE parameters E_max, SD, and AR can be used to evaluate the biological behavior of breast cancer.

Different members of the APOBEC3 family of DNA mutators have opposing associations with the landscape of breast cancer

APOBEC enzymes are strong mutagenic factors. In breast cancer, expression of APOBEC3B is increased and associated with mutation load and poor outcome. Other APOBEC3s can also mutate DNA but their clinical significance in breast cancer and its underpinnings have not been comprehensively studied. In our examination of 1,091 breast carcinoma cases, high expression of APOBEC3A or APOBEC3B genes was associated with greater tumor burden of mutations and other genomic aberrations. Expression of none of the five APOBEC3C-H genes had any correlation with these features, including T[C-T/G]W mutations, but their high expression levels indicated a robust anti-cancer immune response within tumors, with elevated CD8⁺ T cell abundance, T cell receptor diversity, and immune cytolytic activity. Concordantly, survival analyses of this and two other cohorts with > 3,000 patients each showed favorable prognostic benefit of high APOBEC3C-H expression for both cancer progression and mortality. A detrimental prognostic value was observed for APOBEC3A and APOBEC3B. Single-cell data revealed cancer epithelial and stromal immune cells as major sources of APOBEC3B and APOBEC3C-H expression in tumors, respectively. These observations on opposing associations with breast cancer of different APOBEC3s highlight the contrasting roles of these enzymes, promoting cancer through mutagenesis while antagonizing it through immune response.

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.

Preoperative Serum Hypersensitive-c-Reactive-Protein (Hs-CRP) to Albumin Ratio Predicts Survival in Patients with Luminal B Subtype Breast Cancer

Objective

To evaluate the clinical prognostic significance of preoperative serum hypersensitive-c-reactive-protein (Hs-CRP) to albumin ratio (CAR) in patients with luminal B subtype breast cancer.

Methods

A total of 199 patients with luminal B subtype breast cancer enrolled in this study were analyzed retrospectively. The optimal cutoff value of CAR was performed by the receiver operating characteristic curve (ROC). The associations between luminal B subtype breast cancer and clinicopathological variables by CAR were performed by chi-square test. Kaplan–Meier and log rank method were used for survival analysis. The independent prognostic factors were determined by univariate and multivariate Cox’s proportional hazards regression model.

Results

The patients were divided into low CAR group (CAR<0.044) and high CAR group (CAR≥0.044) by ROC. CAR was the independent factor by univariate and multivariate analysis, and the mean DFS and OS in the low CAR group survived longer than those in the high CAR group (p<0.05). According to the endocrine therapy with aromatase inhibitors, the mean survival time of DFS and OS in the low CAR group was significantly higher than that in the high CAR group (p<0.05). Moreover, patients with pathological I+II stage survived longer than those with pathological III stage, and the mean survival time of DFS and OS in the low CAR group was significantly higher than that in the high CAR group (p<0.05). Patients without lymph vessel invasion survived longer than those with lymph vessel invasion (p<0.05), and the mean survival time of DFS and OS in low the CAR group was significantly higher than that in the high CAR group (p<0.05).

Conclusion

Preoperative CAR was significantly associated with survival and prognosis of breast cancer, and it can be used as a routine prognostic indicator to predict the prognosis of luminal B subtype breast 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.

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

Simple Summary

We studied the relationship between enhancement of DNA repair and cancer aggressiveness, tumor immune microenvironment, and patient survival in 749 hepatocellular carcinoma (HCC) patients from 5 cohorts using a DNA repair pathway score. We show that the DNA repair pathway was enhanced by the stepwise carcinogenic process of HCC, notably in grade 3 compared to grade 1 or 2 HCC. DNA repair high HCC was associated with worse survival, elevated intratumor heterogeneity, and mutation load, but not with the fraction of immune cell infiltration nor cytolytic activity. The expression of proliferation- and other cancer aggressiveness-related gene sets was also increased. Interestingly, these features were more pronounced in low-grade compared to high-grade HCC. In conclusion, the DNA repair score may be used to understand the role of DNA repair pathways in patient prognosis and treatment sensitivity and be used to improve patient outcome. To our knowledge, this is the first study using DNA repair pathway-related gene set expression data to examine and validate the clinical relevance of DNA repair pathway activity in HCC.

Abstract

Hepatocellular carcinoma (HCC) is one of the most common malignancies and a leading cause of cancer-related deaths worldwide. In this study, a total of 749 HCC patients from 5 cohorts were studied to examine the relationships between enhancement of DNA repair and cancer aggressiveness, tumor immune microenvironment, and patient survival in HCC, utilizing a DNA repair pathway score. Our findings suggest that the DNA repair pathway was not only enhanced by the stepwise carcinogenic process of HCC, but also significantly enhanced in grade 3 HCC compared with grade 1 and 2 tumors. DNA repair high HCC was associated with worse survival, elevated intratumor heterogeneity, and mutation load, but not with the fraction of immune cell infiltration nor immune response. HCC tumors with a DNA repair high score enriched the cell proliferation- and other cancer aggressiveness-related gene sets. Interestingly, these features were more pronounced in grade 1 and 2 HCC compared to grade 3 HCC. To our knowledge, this is the first study to use DNA repair pathway-related gene set expression data to examine and validate the clinical relevance of DNA repair pathway activity in HCC. The DNA repair score may be used to better understand and predict prognosis in HCC.

Th2 cell infiltrations predict neoadjuvant chemotherapy response of estrogen receptor-positive breast cancer

Background

High infiltration of Th2 is linked to breast cancer progression and metastasis through the induction of cytokine release and T-cell anergy. The estrogen receptor (ER)-positive subtype, which accounts for 70% of breast cancer, is known to respond less to neoadjuvant chemotherapy (NAC) due to its low potential for proliferation. We hypothesized that Th2 high tumors are highly proliferative, and thus more likely to respond to NAC in ER-positive breast cancer.

Methods

We obtained clinicopathological data and overall survival information on 1,069 breast cancer patients from The Cancer Genome Atlas (TCGA). Computational algorithms and CIBERSORT were used to estimate immune cell infiltration. Additionally, xCell was used for validation.

Results

Th2 high tumors did not consistently associate with an unfavorable immune cell composition and tumor immune microenvironment but were found to be significantly elevated in the cancer stage. Th2 high tumors also correlated with high Nottingham pathological grade, as well as with Ki-67 and proliferation score in ER-positive subtypes. High Th2 tumors achieved a pathological complete response (pCR) significantly higher in ER-positive breast cancer.

Conclusions

In conclusion, high levels of Th2 are associated with aggressive features of breast cancer. Th2 levels may be a biomarker in patient selection for NAC in ER-positive breast cancer.

The Impact of Immunofunctional Phenotyping on the Malfunction of the Cancer Immunity Cycle in Breast Cancer

Simple Summary

The cancer-immunity cycle (CIC) is a series of self-sustaining stepwise events to fight cancer growth by the immune system. We hypothesized that immunofunctional phenotyping that represent the malfunction of the CIC is clinically relevant in breast cancer (BC) utilizing total of 2979 BC cases; 1075 from TCGA cohort, 1904 from METABRIC cohort were analyzed. The immunofunctional phenotype was classified as follows: hot T-cell infiltrated, high immune cytolytic activity (CYT), cold T-cell infiltrated, high frequency of CD8+ T cells and low CYT, and non-inflamed, low frequency of CD8+ T cells and low CYT. We demonstrated that immunofunctional phenotyping not only indicated the degree of anti-cancer immune dysfunction, but also served as a prognostic biomarker and HTI was inversely related to estrogen response.

Abstract

The cancer-immunity cycle (CIC) is a series of self-sustaining stepwise events to fight cancer growth by the immune system. We hypothesized that immunofunctional phenotyping that represent the malfunction of the CIC is clinically relevant in breast cancer (BC). Total of 2979 BC cases; 1075 from TCGA cohort, 1904 from METABRIC cohort were analyzed. The immunofunctional phenotype was classified as follows: hot T-cell infiltrated (HTI), high immune cytolytic activity (CYT), Cold T-cell infiltrated (CTI), high frequency of CD8+ T cells and low CYT, and non-inflamed, low frequency of CD8+ T cells and low CYT. The analysis of tumor immune microenvironment in the immunofunctional phenotype revealed that not only immunostimulatory factors, but also immunosuppressive factors were significantly elevated and immunosuppressive cells were significantly decreased in HTI. Patients in HTI were significantly associated with better survival in whole cohort and patients in CTI were significantly associated with worse survival in triple negative. Furthers, HTI was inversely related to estrogen responsive signaling. We demonstrated that immunofunctional phenotype not only indicated the degree of anti-cancer immune dysfunction, but also served as a prognostic biomarker and HTI was inversely related to estrogen response.

High Expression of NRF2 Is Associated with Increased Tumor-Infiltrating Lymphocytes and Cancer Immunity in ER-Positive/HER2-Negative Breast Cancer

Nuclear factor erythroid 2-related factor 2 (NRF2) is a key modifier in breast cancer. It is unclear whether NRF2 suppresses or promotes breast cancer progression. We studied the clinical relevance of NRF2 expression by conducting in silico analyses in 5443 breast cancer patients from several large patient cohorts (METABRIC, GSE96058, GSE25066, GSE20194, and GSE75688). NRF2 expression was significantly associated with better survival, low Nottingham pathological grade, and ER-positive/HER2-negative and triple negative breast cancer (TNBC). High NRF2 ER-positive/HER2-negative breast cancer enriched inflammation- and immune-related gene sets by GSEA. NRF2 expression was elevated in immune, stromal, and cancer cells. High NRF2 tumors were associated with high infiltration of immune cells (CD8+, CD4+, and dendritic cells (DC)) and stromal cells (adipocyte, fibroblasts, and keratinocytes), and with low fraction of Th1 cells. NRF2 expression significantly correlated with area under the curve (AUC) of several drug response in multiple ER-positive breast cancer cell lines, however, there was no significant association between NRF2 and pathologic complete response (pCR) rate after neoadjuvant chemotherapy in human samples. Finally, high NRF2 breast cancer was associated with high expression of immune checkpoint molecules. In conclusion, NRF2 expression was associated with enhanced tumor-infiltrating lymphocytes in ER-positive/HER2-negative breast cancer.

Inflammation Is Associated with Worse Outcome in the Whole Cohort but with Better Outcome in Triple-Negative Subtype of Breast Cancer Patients

Inflammation has been linked with cancer, but whether it is part of the problem or part of the solution remains to be a matter of debate in breast cancer. Our group and others have demonstrated that inflammation aggravates cancer progression; however, some claim that inflammation may support immune cell infiltration and suppress cancer. We defined the gene set variation analysis of the Molecular Signatures Database Hallmark inflammatory response gene set as the inflammatory pathway score and analyzed 3632 tumors in total from 4 breast cancer cohorts (METABRIC, TCGA, GSE25066, and GSE21094). In the whole breast cancer cohort, high-score tumors were associated with aggressive clinical characteristics, such as worse disease specific survival, higher Nottingham histological grade, and younger age. Inflammatory score was significantly higher in triple-negative (TNBC) as well as basal and normal subtypes compared with the other subtypes, which suggest that the detrimental effect of high level of inflammation may be because it includes a more aggressive subtype. On the contrary, high score within TNBC was significantly associated with better survival. TNBC with high score enriched not only IFN-α, IFN-γ response, IL-2/STAT5 signaling, Allograft rejection, Complement, p53 pathway, Reactive Oxygen, and Apoptosis but also TNF-α signaling, IL6-JAK-STAT signaling, TGF-β signaling, Coagulation, Angiogenesis, EMT, KRAS signaling, and PI3K-AKT-MTOR signaling gene sets. High score was associated with mainly favorable anticancerous immune cell infiltration as well as Leukocyte fraction, TIL regional fraction, Lymphocyte infiltration, IFN-γ response, TGF-β response, and cytolytic activity scores. Although the inflammatory pathway score was not associated with neoadjuvant treatment response, it associated with expressions of immune checkpoint molecules. In conclusion, inflammation was associated with worse outcome in the whole breast cancer cohort, but with better outcome in TNBC, which was associated with favorable anticancerous immune response and immune cell infiltrations.

Plasmacytoid Dendritic Cell (pDC) Infiltration Correlate with Tumor Infiltrating Lymphocytes, Cancer Immunity, and Better Survival in Triple Negative Breast Cancer (TNBC) More Strongly than Conventional Dendritic Cell (cDC)

Dendritic cells (DC) represent a major antigen-presenting cell type in the tumor immune microenvironment (TIME) and play an essential role in cancer immunity. Conventional DC (cDC) and plasmacytoid DC (pDC) were defined by the xCell algorithm and a total of 2968 breast cancer patients (TCGA and METABRIC) were analyzed. We found that triple-negative breast cancer (TNBC) had a high fraction of cDC and pDC compared to the other subtypes. In contrast to cDC, high pDC in TNBC was significantly associated with better disease-specific and disease-free survival consistently in both cohorts. High cDC TNBC tumors enriched not only inflammation and immune-related, but also metastasis-related gene sets in Gene Set Enrichment Analysis, whereas high pDC TNBC enriched inflammation and immune -related gene sets including IFN-γ signaling more strongly than cDC. pDC TNBC correlated with CD8+, CD4+ memory, IFN-γ score, and cytolytic activity stronger than cDC TNBC. High pDC TNBC were associated with a high fraction of anti-cancer immune cells and high expression of all the immune check point molecules examined. In conclusion, pDC levels correlated with the infiltration of immune cells and patient survival in TNBC more strongly than cDC; this is the first study suggesting the clinical relevance of pDC infiltration in TNBC.

Abundance of Regulatory T Cell (Treg) as a Predictive Biomarker for Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer

Regulatory CD4+ T cell (Treg), a subset of tumor-infiltrating lymphocytes (TILs), are known to suppress anticancer immunity but its clinical relevance in human breast cancer remains unclear. In this study, we estimated the relative abundance of Tregs in breast cancer of multiple patient cohorts by using the xCell algorithm on bulk tumor gene expression data. In total, 5177 breast cancer patients from five independent cohorts (TCGA-BRCA, GSE96058, GSE25066, GSE20194, and GSE110590) were analyzed. Treg abundance was not associated with cancer aggressiveness, patient survival, or immune activity markers, but it was lower in metastatic tumors when compared to matched primary tumors. Treg was associated with a high mutation rate of TP53 genes and copy number mutations as well as with increased tumor infiltration of M2 macrophages and decreased infiltration of T helper type 1 (Th1) cells. Pathological complete response (pCR) after neoadjuvant chemotherapy (NAC) was significantly associated with low Treg abundance in triple negative breast cancer (TNBC) but not in ER-positive/Her2-negative subtype. High Treg abundance was significantly associated with high tumor expression of multiple immune checkpoint inhibitor genes. In conclusion, Treg abundance may have potential as a predictive biomarker of pCR after NAC in TNBC.

M1 Macrophage and M1/M2 ratio defined by transcriptomic signatures resemble only part of their conventional clinical characteristics in breast cancer

Tumor associated macrophages (TAMs) play a critical role in biology of various cancers, including breast cancer. In the current study, we defined "M1" macrophage and "M1"/"M2" ratio by transcriptomic signatures using xCell. We investigated the association between high level of "M1" macrophage or "M1"/"M2" ratio and the tumor immune microenvironment by analyzing the transcriptome of publicly available cohorts, TCGA and METABRIC. We found that "M1" high tumors were not associated with prolonged survival compared with "M1" low tumors, or with the response to neoadjuvant chemotherapy. "M1" high tumors were associated with clinically aggressive features and "M1" high tumors enriched the cell proliferation and cell cycle related gene sets in GSEA. At the same time, "M1" high tumors were associated with high immune activity and favorable tumor immune microenvironment, as well as high expression of immune check point molecules. Strikingly, all these results were mirrored in "M1"/"M2" ratio high tumors. In conclusion, transcriptomically defined "M1" or "M1"/"M2" high tumors were associated with aggressive cancer biology and favorable tumor immune microenvironment but not with survival benefit, which resembled only part of their conventional clinical characteristics.

High G2M Pathway Score Pancreatic Cancer is Associated with Worse Survival, Particularly after Margin-Positive (R1 or R2) Resection

Pancreatic cancer is highly mortal due to uncontrolled cell proliferation. The G2M checkpoint pathway is an essential part of the cell cycle. We hypothesized that a high G2M pathway score is associated with cell proliferation and worse survival in pancreatic cancer patients. Gene set variation analysis using the Hallmark G2M checkpoint gene set was used as a score to analyze a total of 390 human pancreatic cancer patients from 3 cohorts (TCGA, GSE62452, GSE57495). High G2M score tumors enriched other cell proliferation genes sets as well as MKI67 expression, pathological grade, and proliferation score. Independent of other prognostic factors, G2M score was predictive of disease-specific survival in pancreatic cancer. High G2M tumor was associated with high mutation rate of KRAS and TP53 and significantly enriched these pathway gene sets, as well as high infiltration of Th2 cells. High G2M score consistently associated with worse overall survival in 3 cohorts, particularly in R1/2 resection, but not in R0. High G2M tumor in R1/2 highly enriched metabolic and cellular components' gene sets compared to R0. To our knowledge, this is the first study to use gene set variation analysis as a score to examine the clinical relevancy of the G2M pathway in pancreatic cancer.

ITPKC as a Prognostic and Predictive Biomarker of Neoadjuvant Chemotherapy for Triple Negative Breast Cancer

Triple negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with higher mortality than the others. Pathological complete response (pCR) to neoadjuvant chemotherapy (NAC) is considered as a surrogate to predict survival. Inositol 1,4,5-trisphosphate 3-kinase C (ITPKC) is a negative regulator of T cell activation, and reduction in ITPKC function is known to promote Kawasaki disease. Given the role of tumor infiltrating lymphocytes in NAC and since TNBC has the most abundant immune cell infiltration in breast cancer, we hypothesized that the ITPKC expression level is associated with NAC response and prognosis in TNBC. The ITPKC gene was expressed in the mammary gland, but its expression was highest in breast cancer cells among other stromal cells in a bulk tumor. ITPKC expression was highest in TNBC, associated with its survival, and was its independent prognostic factor. Although high ITPKC was not associated with immune function nor with any immune cell fraction, low ITPKC significantly enriched cell proliferation-related gene sets in TNBC. TNBC with low ITPKC achieved a significantly higher pCR rate after NAC. To the best of our knowledge, this is the first report to demonstrate that ITPKC gene expression may be useful as a prognostic and predictive biomarker in TNBC.

CD8 T Cell Score as a Prognostic Biomarker for Triple Negative Breast Cancer

CD8 T cell is an essential component of tumor-infiltrating lymphocytes (TIL) and tumor immune microenvironment (TIME). Using the xCell CD8 T cell score of whole tumor gene expression data, we estimated these cells in total of 3837 breast cancer patients from TCGA, METABRIC and various GEO cohorts. The CD8 score correlated strongly with expression of CD8 genes. The score was highest for triple-negative breast cancer (TNBC), and a high score was associated with high tumor immune cytolytic activity and better survival in TNBC but not other breast cancer subtypes. In TNBC, tumors with a high CD8 score had enriched expression of interferon (IFN)-α and IFN-γ response and allograft rejection gene sets, and greater infiltration of anti-cancerous immune cells. The score strongly correlated with CD4 memory T cells in TNBC, and tumors with both a high CD8 score and high CD4 memory T cell abundance had significantly better survival. Finally, a high CD8 score was significantly associated with high expression of multiple immune checkpoint molecules. In conclusion, a high CD8 T cell score is associated with better survival in TNBC, particularly when tumor CD4 memory T cells were elevated. Our findings also suggest a possible use of the score as a predictive biomarker for response to immune checkpoint therapy.

Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation

Cancer-associated adipocytes are known to cause inflammation, leading to cancer progression and metastasis. The clinicopathological and transcriptomic data from 2256 patients with breast cancer were obtained based on three cohorts: The Cancer Genome Atlas (TCGA), GSE25066, and a study by Yau et al. For the current study, we defined the adipocyte, which is calculated by utilizing a computational algorithm, xCell, as "intratumoral adipocyte". These intratumoral adipocytes appropriately reflected mature adipocytes in a bulk tumor. The amount of intratumoral adipocytes demonstrated no relationship with survival. Intratumoral adipocyte-high tumors significantly enriched for metastasis and inflammation-related gene sets and are associated with a favorable tumor immune microenvironment, especially in the ER+/HER2- subtype. On the other hand, intratumoral adipocyte-low tumors significantly enriched for cell cycle and cell proliferation-related gene sets. Correspondingly, intratumoral adipocyte-low tumors are associated with advanced pathological grades and inversely correlated with MKI67 expression. In conclusion, a high amount of intratumoral adipocytes in breast cancer was associated with inflammation, metastatic pathways, cancer stemness, and favorable tumor immune microenvironment. However, a low amount of adipocytes was associated with a highly proliferative tumor in ER-positive breast cancer. This cancer biology may explain the reason why patient survival did not differ by the amount of adipocytes.

Transcriptomic Profile of Lymphovascular Invasion, a Known Risk Factor of Pancreatic Ductal Adenocarcinoma Metastasis

Lymphovascular invasion (LVI) is an aggressive pathologic feature and considered a risk factor for distant metastasis. We hypothesized that pancreatic ductal adenocarcinomas (PDACs) with LVI are associated with shorter survival, as well as aggressive cancer biology and lymphangiogenesis in transcriptomic analysis. Utilizing the cancer genome atlas (TCGA)-PDAC cohort, we found that positive LVI was significantly associated with positive perineural invasion (PNI) (p = 0.023), and higher American Joint Committee on Cancer (AJCC) T (p = 0.017) and N (p < 0.001) categories. Furthermore, positive LVI was associated with shorter overall survival (OS) (p = 0.014) and was an independent risk factor of poor OS. Although there was no association between LVI status and lymphangiogenesis, epithelial-mesenchymal transition (EMT), or metastasis-related genes, Gene Set Enrichment Analysis revealed a strong association with cell-proliferation-related gene sets such as mitotic spindles (Normalized enrichment score (NES) = 1.76, p = 0.016) and G2/M checkpoints (NES = 1.75, p = 0.036), as well as with transforming growth factor beta (TGF-beta) signaling (NES = 1.61, p = 0.043), which is a known mechanism of metastasis in PDACs. In conclusion, positive LVI was an independent risk factor of poor OS in PDACs. We found that PDACs with LVI were possibly associated with accelerated cell proliferation and enhanced TGF-beta signaling independent of lymphangiogenesis. Transcriptomic profiling elucidates more precise tumor biology of LVI-positive PDACs.

The E2F Pathway Score as a Predictive Biomarker of Response to Neoadjuvant Therapy in ER+/HER2- Breast Cancer

E2F transcription factors play critical roles in the cell cycle. Therefore, their activity is expected to reflect tumor aggressiveness and responsiveness to therapy. We scored 3905 tumors of nine breast cancer cohorts for this activity based on their gene expression for the Hallmark E2F targets gene set. As expected, tumors with a high score had an increased expression of cell proliferation-related genes. A high score was significantly associated with shorter patient survival, greater MKI67 expression, histological grade, stage, and genomic aberrations. Furthermore, metastatic tumors had higher E2F scores than the primary tumors from which they arose. Although tumors with a high score had greater infiltration by both pro- and anti-cancerous immune cells, they had an increased expression of immune checkpoint genes. Estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative cancer with a high E2F score achieved a significantly higher pathological complete response (pCR) rate to neoadjuvant chemotherapy. The E2F score was significantly associated with the expression of cyclin-dependent kinase (CDK)-related genes and strongly correlated with sensitivity to CDK inhibition in cell lines. In conclusion, the E2F score is a marker of breast cancer aggressiveness and predicts the responsiveness of ER-positive/HER2-negative patients to neoadjuvant chemotherapy and possibly to CDK and immune checkpoint inhibitors.

Contribution of Immune Cells to Glucocorticoid Receptor Expression in Breast Cancer

Breast cancer (BC) patients experience increased stress with elevated cortisol levels, increasing risk of cancer recurrence. Cortisol binds to a cytoplasmic receptor, glucocorticoid receptor (GR) encoded by GR gene (NR3C1). We hypothesized that not only cancer cells, but even immune cells in the tumor microenvironment (TME) may contribute to GR expression in bulk tumor and influence prognosis. To test this, mRNA expression data was accessed from METABRIC and TCGA. "High" and "low" expression was based on highest and lowest quartiles of NR3C1 gene expression, respectively. Single-cell sequencing data were obtained from GSE75688 and GSE114725 cohorts. Computer algorithms CIBERSORT, Gene Set Enrichment Analysis and TIMER were used. GR-high BC has better median disease-free and disease-specific survival. Single cell sequencing data showed higher GR expression on immune cells compared to cancer and stromal cells. Positive correlation between GR-high BC and CD8+ T-cells was noted. In GR-high tumors, higher cytolytic activity (CYT) with decreased T-regulatory and T-follicular helper cells was observed. High GR expression was associated with lower proliferation index Ki67, enriched in IL-2_STAT5, apoptosis, KRAS, TGF-β signaling, and epithelial-to-mesenchymal transition. Immune cells significantly contribute to GR expression of bulk BC. GR-high BC has a favorable TME with higher CYT with favorable outcomes.

Fibroblasts as a Biological Marker for Curative Resection in Pancreatic Ductal Adenocarcinoma

Achievement of microscopic tumor clearance (R0) after pancreatic ductal adenocarcinoma (PDAC) surgery is determined by cancer biology rather than operative technique. Fibroblasts are known to play pro-cancer roles; however, a small subset was recently found to play anti-cancer roles. Therefore, we hypothesized that intratumor fibroblasts contribute to curative resection and a better survival of PDAC. Utilizing a large, publicly available PDAC cohort, we found that fibroblast composition was associated with R0 curative resection. A high amount of fibroblasts in PDACs was significantly associated with a higher amount of mature vessels, but not with blood angiogenesis. A high amount of fibroblasts was also associated with a higher infiltration of anti-cancer immune cells, such as CD8+ T-cells and dendritic cells, together with higher inflammatory signaling, including IL2/STAT5 and IL6/JAK/STAT3 signaling. Further, the fibroblast composition was inversely associated with cancer cell composition in the bulk tumor, along with an inverse association with proliferative characteristics, such as MYC signaling and glycolysis. The patients with high-fibroblast PDACs showed an improved prognosis. In conclusion, we found that PDACs with high fibroblasts were associated with a higher R0 resection rate, resulting in a better prognosis. These findings may be due to less aggressive biology with a higher vascularity and anti-cancer immunity, and a low cancer cell component.