High Expression of microRNA-143 is Associated with Favorable Tumor Immune Microenvironment and Better Survival in Estrogen Receptor Positive Breast Cancer

Stable Dual miR-143 and miR-506 Upregulation Inhibits Proliferation and Cell Cycle Progression

The mainstays of lung cancer pathogenesis are cell cycle progression dysregulation, impaired apoptosis, and unregulated cell proliferation. While individual microRNA (miR) targeting or delivering is a promising approach that has been extensively studied, combination of miR targeting can enhance therapeutic efficacy and overcome limitations present in individual miR regulations. We previously reported on the use of a miR-143 and miR-506 combination via transient transfections against lung cancer. In this study, we evaluated the effect of miR-143 and miR-506 under stable deregulations in A549 lung cancer cells. We used lentiviral transductions to either up- or downregulate the two miRs individually or in combination. The cells were sorted and analyzed for miR deregulation via qPCR. We determined the miR deregulations' effects on the cell cycle, cell proliferation, cancer cell morphology, and cell motility. Compared to the individual miR deregulations, the combined miR upregulation demonstrated a miR-expression-dependent G2 cell cycle arrest and a significant increase in the cell doubling time, whereas the miR-143/506 dual downregulation demonstrated increased cellular motility. Furthermore, the individual miR-143 and miR-506 up- and downregulations exhibited cellular responses lacking an apparent miR-expression-dependent response in the respective analyses. Our work here indicates that, unlike the individual miR upregulations, the combinatorial miR treatment remained advantageous, even under prolonged miR upregulation. Finally, our findings demonstrate potential advantages of miR combinations vs. individual miR treatments.

Identification of stably expressed microRNAs in plasma from high-grade serous ovarian carcinoma and benign tumor patients

BACKGROUND

Ovarian cancer is a lethal gynecological cancer and no reliable minimally invasive early diagnosis tools exist. High grade serous ovarian carcinoma (HGSOC) is often diagnosed at advanced stages, resulting in poorer outcome than those diagnosed in early stage. Circulating microRNAs have been investigated for their biomarker potential. However, due to lack of standardization methods for microRNA detection, there is no consensus, which microRNAs should be used as stable endogenous controls. We aimed to identify microRNAs that are stably expressed in plasma of HGSOC and benign ovarian tumor patients.

METHODS AND RESULTS

We isolated RNA from plasma samples of 60 HGSOC and 48 benign patients. RT-qPCR was accomplished with a custom panel covering 40 microRNAs and 8 controls. Stability analysis was performed using five algorithms: Normfinder, geNorm, Delta-Ct, BestKeeper and RefFinder using an R-package; RefSeeker developed by our study group [1]. Among 41 analyzed RNAs, 13 were present in all samples and eligible for stability analysis. Differences between stability rankings were observed across algorithms. In HGSOC samples, hsa-miR-126-3p and hsa-miR-23a-3p were identified as the two most stable miRNAs. In benign samples, hsa-miR-191-5p and hsa-miR-27a-3p were most stable. In the combined HGSOC and benign group, hsa-miR-23a-3p and hsa-miR-27a-3p were identified by both the RefFinder and Normfinder analysis as the most stable miRNAs.

CONCLUSIONS

Consensus regarding normalization approaches in microRNA studies is needed. The choice of endogenous microRNAs used for normalization depends on the histological content of the cohort. Furthermore, normalization also depends on the algorithms used for stability analysis.

Metastasis Inhibition

Cancer metastasis is a common biological phenomenon observed in malignant tumors that can lead to death in affected individuals [...].

Strategies for data normalization and missing data imputation and consequences for potential diagnostic microRNA biomarkers in epithelial ovarian cancer

MicroRNAs (miRNAs) are small non-coding RNA molecules regulating gene expression with diagnostic potential in different diseases, including epithelial ovarian carcinomas (EOC). As only a few studies have been published on the identification of stable endogenous miRNA in EOC, there is no consensus which miRNAs should be used aiming standardization. Currently, U6-snRNA is widely adopted as a normalization control in RT-qPCR when investigating miRNAs in EOC; despite its variable expression across cancers being reported. Therefore, our goal was to compare different missing data and normalization approaches to investigate their impact on the choice of stable endogenous controls and subsequent survival analysis while performing expression analysis of miRNAs by RT-qPCR in most frequent subtype of EOC: high-grade serous carcinoma (HGSC). 40 miRNAs were included based on their potential as stable endogenous controls or as biomarkers in EOC. Following RNA extraction from formalin-fixed paraffin embedded tissues from 63 HGSC patients, RT-qPCR was performed with a custom panel covering 40 target miRNAs and 8 controls. The raw data was analyzed by applying various strategies regarding choosing stable endogenous controls (geNorm, BestKeeper, NormFinder, the comparative ΔCt method and RefFinder), missing data (single/multiple imputation), and normalization (endogenous miRNA controls, U6-snRNA or global mean). Based on our study, we propose hsa-miR-23a-3p and hsa-miR-193a-5p, but not U6-snRNA as endogenous controls in HGSC patients. Our findings are validated in two external cohorts retrieved from the NCBI Gene Expression Omnibus database. We present that the outcome of stability analysis depends on the histological composition of the cohort, and it might suggest unique pattern of miRNA stability profiles for each subtype of EOC. Moreover, our data demonstrates the challenge of miRNA data analysis by presenting various outcomes from normalization and missing data imputation strategies on survival analysis.

Molecular Mechanisms of Anti-Estrogen Therapy Resistance and Novel Targeted Therapies

Breast cancer (BC) is the most commonly diagnosed cancer in women, constituting one-third of all cancers in women, and it is the second leading cause of cancer-related deaths in the United States. Anti-estrogen therapies, such as selective estrogen receptor modulators, significantly improve survival in estrogen receptor-positive (ER+) BC patients, which represents about 70% of cases. However, about 60% of patients inevitably experience intrinsic or acquired resistance to anti-estrogen therapies, representing a major clinical problem that leads to relapse, metastasis, and patient deaths. The resistance mechanisms involve mutations of the direct targets of anti-estrogen therapies, compensatory survival pathways, as well as alterations in the expression of non-coding RNAs (e.g., microRNA) that regulate the activity of survival and signaling pathways. Although cyclin-dependent kinase 4/6 and phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) inhibitors have significantly improved survival, the efficacy of these therapies alone and in combination with anti-estrogen therapy for advanced ER+ BC, are not curative in advanced and metastatic disease. Therefore, understanding the molecular mechanisms causing treatment resistance is critical for developing highly effective therapies and improving patient survival. This review focuses on the key mechanisms that contribute to anti-estrogen therapy resistance and potential new treatment strategies alone and in combination with anti-estrogen drugs to improve the survival of BC patients.

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.

Development and validation of an extended Cox prognostic model for patients with ER/PR+ and HER2- breast cancer: a retrospective cohort study

Background

The purpose of this study was to explore a new estrogen receptor (ER) and/or progesterone receptor (PR)+ and human epidermal growth factor receptor 2 (HER2)− breast cancer prognostic model, called the extended Cox prognostic model, for determining the cutoff values for multiple continuous prognostic factors and their interaction via the new model concept and variable selection method.

Methods

A total of 335 patients with ER/PR+ and HER2− breast cancer were enrolled for the final analysis. The primary endpoint was breast cancer-specific mortality (BCSM). Prognostic factors (histological grade, histological type, stage, T, N, lymphovascular invasion (LVI), P53, Ki67, ER, PR, and age) were included in this study. The four continuous variables (Ki67, ER, PR, and age) were partitioned into a series of binary variables that were fitted in the multivariate Cox analysis. A smoothly clipped absolute deviation (SCAD) variable selection method was used. Model performance was expressed in discrimination and calibration.

Results

We developed an extended Cox model with a time threshold of 164-week (more than 3 years) postoperation and developed a user-friendly nomogram based on our extended Cox model to facilitate clinical application. We found that the cutoff values for PR, Ki67, and age were 20%, 60%, and 41–55 years, respectively. There was an interaction between age and PR for patients aged ≥ 41 years and PR ≥ 20% at 164-week postoperation: the older the patients with ER/PR+, HER2−, and PR ≥ 20% were, the lower the survival and more likely to recur and metastasize exceeding 164 weeks (more than 3 years) after surgery.

Conclusions

Our study offers guidance on the prognosis of patients with ER/PR+ and HER2− breast cancer in China. The new concept can inform modeling and the determination of cutoff values of prognostic factors in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12957-022-02790-0.

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.

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

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

An immune-inflamed tumor microenvironment as defined by CD8 score is associated with favorable oncologic outcomes in hepatocellular carcinoma independent of measures of tumor mutational burden

Despite low mutational burden, immune checkpoint inhibitors have demonstrated promising results in a significant minority of hepatocellular carcinoma (HCC) patients with advanced disease. We hypothesized that HCC patients with higher levels of CD8+ T cell infiltration reflect an immune-inflamed cohort which has improved oncologic outcomes. 355 HCC patients with clinical and transcriptome data in the Cancer Genome Atlas (TCGA) and 151 HCC patients from cohort GSE7624 were analyzed. xCell computational algorithm was used to analyze immune cell infiltration in these patients. Each cohort was divided into high and low expression by the highest 2 terciles value. Gene Set Enrichment Analysis was performed to identify enriched gene sets. High CD8 score associated with improved overall survival in both cohorts (both P < 0.05). High score correlates with early BCLC stage (P = 0.035) but not AJCC stage. High CD8 also correlated with increased IFN-γ response (p = 0.038), lymphocyte infiltration (P < 0.001), and leukocyte fraction (P < 0.001). It was associated with increased polyclonality of T cell (P < 0.001) and B cell response (P = 0.017). High CD8 score correlated with increased cytolytic activity score (P < 0.001) and expression of multiple immune checkpoints including PD-1, PD-L1, CTLA-4 and Lag3 (all P < 0.001). There was no correlation to tumor mutational burden and neoantigens. GSEA demonstrated upregulation of several gene sets involved in inflammatory response and IFN-γ response. In conclusion, HCC patients with high CD8 score demonstrated favorable oncologic outcomes, which may be due to immune-mediated tumor cell attack. Furthermore, CD8 score may be a potentially useful biomarker to select patients for immune checkpoint inhibition.

A novel immune-related microRNA signature for prognosis of thymoma

INTRODUCTION

Immune microenvironment and microRNAs serve as common predictors for diagnosis and prognosis of tumors.

METHODS

Expression of 122 genes and 126 microRNAs in thymoma was obtained from TCGA database. The proportion of tumor-infiltrating cells was calculated, and IMRS was constructed. TREM2hi score was calculated before functional enrichment analysis on gene sets.

RESULTS

IMRS3, TREM2hi score, and CD8+ T lymphocyte abundance were significantly different among WHO classifications. WHO classification, Masaoka staging, and miR-130b-5p, miR-1307-3p, miR-425-5p, CD8, CD68, and CCL18 expression were prognostic factors for relapse-free survival and overall survival. IMRS3 upregulation polarized macrophages into M2, which rejected CD8+ T and other effector lymphocytes to promote thymoma malignant progression.

CONCLUSIONS

BRRS may present a novel immune-related microRNA signature for TET prognosis.

Epigenetic Regulation of Estrogen Receptor Genes' Expressions in Adipose Tissue in the Course of Obesity

Estrogen affects adipose tissue function. Therefore, this study aimed at assessing changes in the transcriptional activity of estrogen receptor (ER) α and β genes (ESR1 and ESR2, respectively) in the adipose tissues of obese individuals before and after weight loss and verifying whether epigenetic mechanisms were involved in this phenomenon. ESR1 and ESR2 mRNA and miRNA levels were evaluated using real-time PCR in visceral (VAT) and subcutaneous adipose tissue (SAT) of 78 obese (BMI > 40 kg/m2) and 31 normal-weight (BMI = 20−24.9 kg/m2) individuals and in 19 SAT samples from post-bariatric patients. ESR1 and ESR2 methylation status was studied using the methylation-sensitive digestion/real-time PCR method. Obesity was associated with a decrease in mRNA levels of both ERs in SAT (p < 0.0001) and ESR2 in VAT (p = 0.0001), while weight loss increased ESR transcription (p < 0.0001). Methylation levels of ESR1 and ESR2 promoters were unaffected. However, ESR1 mRNA in the AT of obese subjects correlated negatively with the expression of hsa-miR-18a-5p (rs = −0.444), hsa-miR-18b-5p (rs = −0.329), hsa-miR-22-3p (rs = −0.413), hsa-miR-100-5p (rs = −0.371), and hsa-miR-143-5p (rs = −0.289), while the expression of ESR2 in VAT correlated negatively with hsa-miR-576-5p (rs = −0.353) and in SAT with hsa-miR-495-3p (rs = −0.308). In conclusion, obesity-associated downregulation of ER mRNA levels in adipose tissue may result from miRNA interference.

Identifying and Validating of an Autophagy-Related Gene Signature for the Prediction of Early Relapse in Breast Cancer

BACKGROUND

Compelling evidence has demonstrated the pivotal role of autophagy in the prognosis of breast cancer. Breast cancer (BC) patients with early relapse consistently exhibited worse survival.

METHODS

The autophagy-related genes were derived from the Human Autophagy Database (HADb) and high-sequencing data were obtained from The Cancer Genome Atlas (TCGA). Discrepantly expressed autophagy genes (DEAGs) between early relapse and long-term survival groups were performed using the Linear Models for Microarray data (LIMMA) method. Lasso Cox regression analysis was conducted for the selection of the 4-gene autophagy-related gene signature. GSE42568 and GSE21653 databases were enrolled in this study for the external validation of the signature. Then patients were divided into high and low-risk groups based on the specific score formula. GSEA was used to discover the related signaling pathway. The Kaplan-Meier curves and the receiver operating characteristic (ROC) curves were used to evaluate the discrimination and accuracy of the 4-gene signature.

RESULTS

A signature composed of four autophagy-related mRNA including APOL1, HSPA8, SIRT1, and TP73, was identified as significantly associated with the early relapse in BC patients. Time-dependent receiver-operating characteristic at 1 year suggested remarkable accuracy of the signature [area under the curve (AUC = 0.748)]. The risk score model based on the autophagy-related signature showed favorable predicting value in 1-, 2-, and 3-year relapse-free survival (RFS) in training and two validating cohorts. The GSEA displayed gene sets were remarkably enriched in carcinogenic activation pathways and autophagy-related pathways. The nomogram involving three variables (progesterone receptor status, T stage, and 4-gene signature) exhibited relatively good discrimination with a C-index of 0.766.

CONCLUSIONS

Our study establishes an autophagy-related 4-gene signature that can effectively stratify the high-risk and low-risk BC patients for early relapse. Combined with the clinicopathological variables, the signature could significantly help oncologists tailor more efficient treatment strategies for BC patients.

Intratumoral lymphatic endothelial cell infiltration reflecting lymphangiogenesis is counterbalanced by immune responses and better cancer biology in the breast cancer tumor microenvironment

Lymphangiogenesis, the generation of new lymphatic vessels from existing ones, results from the dynamic interactions of lymphatic endothelial cells and the tumor microenvironment (TME). It is well known that lymphangiogenesis occurs during the initial stage of metastasis in various types of malignant tumors. However, it is currently not used as a biomarker partially because gold standard method to quantify it is labor and cost intensive. We hypothesized that the quantity of intratumoral lymphatic endothelial cells (iLECs) in the TME is an indicator of lymphangiogenesis and a predictor of metastatic potential and overall survival in breast cancer. We analyzed a total of 4145 breast cancer patients from the Cancer Genome Atlas (TCGA) and GSE96058 by quantifying their iLECs using the xCell algorithm, and correlated these scores with patient survival, tumor grade, and cancer stage. We also assessed various pro- and anti-cancer gene sets for each tumor to characterize tumor behavior and aggressiveness. As we expected, high-iLEC breast cancer demonstrated enriched lymphoangiogenesis and angiogenesis gene sets and was associated with increased expressions of related genes. Also enriched were inflammatory response and immune response-related gene sets; IL2/STAT5 pathway, IL6/JAK/STAT3 pathway, TNFα pathway, allograft rejection, and complement as well as cancer stemness related gene sets like Notch signaling, Hedgehog signaling, epithelial mesenchymal transition, and Wnt beta-catenin signaling. Tumors with high-iLEC showed higher proportions of stromal cells and fewer anti-cancer immune cells. On the other hand, iLEC score did not correlate with patient survival or lymph node metastasis. Surprisingly, breast cancers with fewer iLECs demonstrated enriched E2F Targets, G2M Checkpoint, MYC Targets v1, and MTORC1 signaling which are cancer cell proliferation-related gene sets and exhibited an abundance of pro-cancer immune cells. The amount of iLEC correlated inversely with Ki67 expression and histological grade, which is in agreement that low-iLEC breast cancer was associated with enhanced cancer cell proliferation. In conclusion, while iLECs can be used as a surrogate for lymphangiogenesis in breast cancer, low-iLEC tumors also exhibit features which correspond to aggressive tumor biology, which may explain why the amount of iLECs was not associated with patient survival in our cohorts.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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

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

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

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

The Role of MicroRNA as Clinical Biomarkers for Breast Cancer Surgery and Treatment

Breast cancer is the most common cancer diagnosed in women. In recent times, survival outcomes have improved dramatically in accordance with our enhanced understanding of the molecular processes driving breast cancer proliferation and development. Refined surgical approaches, combined with novel and targeted treatment options, have aided the personalisation of breast cancer patient care. Despite this, some patients will unfortunately succumb to the disease. In recent times, translational research efforts have been focused on identifying novel biomarkers capable of informing patient outcome; microRNAs (miRNAs) are small non-coding molecules, which regulate gene expression at a post-transcriptional level. Aberrant miRNA expression profiles have been observed in cancer proliferation and development. The measurement and correlation of miRNA expression levels with oncological outcomes such as response to current conventional therapies, and disease recurrence are being investigated. Herein, we outline the clinical utility of miRNA expression profiles in informing breast cancer prognosis, predicting response to treatment strategies as well as their potential as therapeutic targets to enhance treatment modalities in the era of precision oncology.

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

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

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

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

Low expression of miR-195 is associated with cell proliferation, glycolysis and poor survival in estrogen receptor (ER)-positive but not in triple negative breast cancer

MiR-195 is a tumor suppressive microRNA in breast cancer. Its clinical relevance remains debatable as it has only been studied via in vitro experiments or small cohort studies. We analyzed a total of 2,038 patients in the TCGA and METABRIC cohorts to assess whether low miR-195 expressing tumors are associated with aggressive cancer characteristics and poor prognostic outcomes. The median cutoff of miR-195 expression was used to split the groups into miR-195 high and low groups. Low miR-19 expressing tumors demonstrated high cell proliferating features by enriching the gene sets associated with cell proliferation, MKI67 expression and pathological grade. One-third of the top target miR-195 genes were related to cell proliferation. Low miR-195 expressing tumors were associated with both pro-cancerous and anti-cancerous immune cells. Low miR-195 expressing tumors were associated with enhanced glycolysis and poor survival in ER-positive tumors, but not other subtypes of breast cancer. In conclusion, low expression of miR-195 in ER-positive breast cancer was associated with enhanced cancer cell proliferation, glycolysis, and worse overall survival.

Urine as a Source of Liquid Biopsy for Cancer

Simple Summary

Tissue biopsy is essential for diagnosis and characterization of a tumor. Recently circulating tumor cells and other tumor-derived nucleic acid can be detected from blood, which is called liquid biopsy. Now this concept has been expanded to many other body fluids including urine. Urine is the least invasive method to obtain a liquid biopsy and can be done anywhere, which allows longitudinal repeated sampling. Here, we review the latest update on urine liquid biopsy in urological and non-urological cancers.

Abstract

Tissue biopsy is the gold standard for diagnosis and morphological and immunohistochemical analyses to characterize cancer. However, tissue biopsy usually requires an invasive procedure, and it can be challenging depending on the condition of the patient and the location of the tumor. Even liquid biopsy analysis of body fluids such as blood, saliva, gastric juice, sweat, tears and cerebrospinal fluid may require invasive procedures to obtain samples. Liquid biopsy can be applied to circulating tumor cells (CTCs) or nucleic acids (NAs) in blood. Recently, urine has gained popularity due to its less invasive sampling, ability to easily repeat samples, and ability to follow tumor evolution in real-time, making it a powerful tool for diagnosis and treatment monitoring in cancer patients. With the development and advancements in extraction methods of urinary substances, urinary NAs have been found to be closely related to carcinogenesis, metastasis, and therapeutic response, not only in urological cancers but also in non-urological cancers. This review mainly highlights the components of urine liquid biopsy and their utility and limitations in oncology, especially in non-urological cancers.

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.

Angiogenesis is associated with an attenuated tumor microenvironment, aggressive biology, and worse survival in gastric cancer patients

Angiogenesis is a cornerstone of cancer as it allows tumors to receive oxygen and nutrients. A high level of angiogenesis within a tumor may therefore be indicative of its aggressiveness. In this study, we examined this hypothesis in gastric cancer. Gene set variation analysis was used to measure the level of angiogenesis in tumors in 1,348 gastric cancer patients using the Hallmark_angiogenesis gene set to score tumor transcriptomes. As we predicted, there was a significant correlation between angiogenesis score and expression of angiogenesis-related genes. The score moderately correlated with abundance of vessel-related stromal cells, fibroblasts and chondrocytes in the tumor microenvironment (TME). Tumors with high score had low infiltration of T helper type 1 and 2 cells but a greater infiltration of M1 macrophages and dendritic cells. They also had enriched expression of gene sets for coagulation, hypoxia, epithelial mesenchymal transition (EMT), and TGF-β signaling. High angiogenesis score was significantly associated with advanced AJCC stage and higher T- but not N-parameters in the TNM staging system. Patients with a high score also had shorter survival. In conclusion, bulk tumor transcriptome-based quantification of tumor angiogenesis using a computational algorithm may serve to identify patients with worse survival in gastric cancer.

Prognostic and Functional Significant of Heat Shock Proteins (HSPs) in Breast Cancer Unveiled by Multi-Omics Approaches

Simple Summary

In this study, we investigated the expression pattern and prognostic significance of the heat shock proteins (HSPs) family members in breast cancer (BC) by using several bioinformatics tools and proteomics investigations. Our results demonstrated that, collectively, HSPs were deregulated in BC, acting as both oncogene and onco-suppressor genes. In particular, two different HSP-clusters were significantly associated with a poor or good prognosis. Interestingly, the HSPs deregulation impacted gene expression and miRNAs regulation that, in turn, affected important biological pathways involved in cell cycle, DNA replication, and receptors-mediated signaling. Finally, the proteomic identification of several HSPs members and isoforms revealed much more complexity of HSPs roles in BC and showed that their expression is quite variable among patients. In conclusion, we elaborated two panels of HSPs that could be further explored as potential biomarkers for BC progression and prognosis.

Abstract

Heat shock proteins (HSPs) are a well-characterized molecular chaperones protein family, classified into six major families, according to their molecular size. A wide range of tumors have been shown to express atypical levels of one or more HSPs, suggesting that they could be used as biomarkers. However, the collective role and the possible coordination of HSP members, as well as the prognostic significance and the functional implications of their deregulated expression in breast cancer (BC) are poorly investigated. Here, we used a systematic multi-omics approach to assess the HSPs expression, the prognostic value, and the underlying mechanisms of tumorigenesis in BC. By using data mining, we showed that several HSPs were deregulated in BC and significantly correlated with a poor or good prognosis. Functional network analysis of HSPs co-expressed genes and miRNAs highlighted their regulatory effects on several biological pathways involved in cancer progression. In particular, these pathways concerned cell cycle and DNA replication for the HSPs co-expressed genes, and miRNAs up-regulated in poor prognosis and Epithelial to Mesenchymal Transition (ETM), as well as receptors-mediated signaling for the HSPs co-expressed genes up-regulated in good prognosis. Furthermore, the proteomic expression of HSPs in a large sample-set of breast cancer tissues revealed much more complexity in their roles in BC and showed that their expression is quite variable among patients and confined into different cellular compartments. In conclusion, integrative analysis of multi-omics data revealed the distinct impact of several HSPs members in BC progression and indicate that collectively they could be useful as biomarkers and therapeutic targets for BC management.

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.

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.

Degree of Early Estrogen Response Predict Survival after Endocrine Therapy in Primary and Metastatic ER-Positive Breast Cancer

Endocrine therapy is the gold-standard treatment for ER-positive/HER2-negative breast cancer. Although its clear benefit, patient compliance is poor (50-80%) due to its long administration period and adverse effects. Therefore, a predictive biomarker that can predict whether endocrine therapy is truly beneficial may improve patient compliance. In this study, we use estrogen response early gene sets of gene set enrichment assay algorithm as the score. We hypothesize that the score could predict the response to endocrine therapy and survival of breast cancer patients. A total of 6549 breast cancer from multiple patient cohorts were analyzed. The score was highest in ER-positive/HER2-negative compared to the other subtypes. Earlier AJCC stage, as well as lower Nottingham pathological grade, were associated with a high score. Low score tumors enriched only allograft rejection gene set, and was significantly infiltrated with immune cells, and high cytolytic activity score. A low score was significantly associated with a worse response to endocrine therapy and worse survival in both primary and metastatic breast cancer patients. The hazard ratio was double that of ESR1 expression. In conclusion, the estrogen response early score predicts response to endocrine therapy and is associated with survival in primary and metastatic breast cancer.

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.

Intra-Tumoral Angiogenesis Is Associated with Inflammation, Immune Reaction and Metastatic Recurrence in Breast Cancer

Angiogenesis is one of the hallmarks of cancer. We hypothesized that intra-tumoral angiogenesis correlates with inflammation and metastasis in breast cancer patients. To test this hypothesis, we generated an angiogenesis pathway score using gene set variation analysis and analyzed the tumor transcriptome of 3999 breast cancer patients from The Cancer Genome Atlas Breast Cancer (TCGA-BRCA), Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), GSE20194, GSE25066, GSE32646, and GSE2034 cohorts. We found that the score correlated with expression of various angiogenesis-, vascular stability-, and sphingosine-1-phosphate (S1P)-related genes. Surprisingly, the angiogenesis score was not associated with breast cancer subtype, Nottingham pathological grade, clinical stage, response to neoadjuvant chemotherapy, or patient survival. However, a high score was associated with a low fraction of both favorable and unfavorable immune cell infiltrations except for dendritic cell and M2 macrophage, and with Leukocyte Fraction, Tumor Infiltrating Lymphocyte Regional Fraction and Lymphocyte Infiltration Signature scores. High-score tumors had significant enrichment for unfavorable inflammation-related gene sets (interleukin (IL)6, and tumor necrosis factor (TNF)α- and TGFβ-signaling), as well as metastasis-related gene sets (epithelial mesenchymal transition, and Hedgehog-, Notch-, and WNT-signaling). High score was significantly associated with metastatic recurrence particularly to brain and bone. In conclusion, using the angiogenesis pathway score, we found that intra-tumoral angiogenesis is associated with immune reaction, inflammation and metastasis-related pathways, and metastatic recurrence in breast cancer.

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.

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.