Background
Elevated tumor-infiltrating lymphocytes (TILs) within the tumor microenvironment is a known positive prognostic factor in colorectal cancer (CRC). We hypothesized that since cytotoxic T ...cells release cytolytic proteins such as perforin (PRF1) and pro-apoptotic granzymes (GZMA) to attack cancer cells, a cytolytic activity score (CYT) would be a useful tool to assess anticancer immunity.
Methods
Genomic expression data were obtained from 456 patients from The Cancer Genome Atlas (TCGA). CYT was defined by GZMA and PRF1 expression, and CIBERSORT was used to evaluate intratumoral immune cell composition.
Results
High CYT was associated with high microsatellite instability (MSI-H), as well as high levels of activated memory CD4+T cells, gamma-delta T cells, and M1 macrophages. CYT-high CRC patients had improved overall survival (
p
= 0.019) and disease-free survival (
p
= 0.016) compared with CYT-low CRC patients, especially in TIL-positive tumors. Multivariate analysis demonstrated that CYT- high associates with improved survival independently after controlling for age, lymphovascular invasion, colonic location, microsatellite instability, and TIL positivity. The levels of immune checkpoint molecules (ICMs)—programmed death-1 (PD-1), programmed death-ligand 1 (PD-L1), cytotoxic T-lymphocyte-associated protein 4 (CTLA4), lymphocyte-activation gene 3 (LAG3), T cell immunoglobulin and mucin domain 3 (TIM3), and indoleamine 2,3-dioxygenase 1 (IDO1)—correlated significantly with CYT (
p
< 0.0001); with improved survival in CYT-high and ICM-low patients, and poorer survival in ICM-high patients.
Conclusions
High CYT within CRC is associated with improved survival, likely due to increased immunity and cytolytic activity of T cells and M1 macrophages. High CYT is also associated with high expression of ICMs; thus, further studies to elucidate the role of CYT as a predictive biomarker of the efficacy of immune checkpoint blockade are warranted.
Background
Intratumor heterogeneity implies that subpopulations of cancer cells that differ in genetic, phenotypic, or behavioral characteristics coexist in a single tumor (Ma in Breast Cancer Res ...Treat 162(1):39–48,
2017
; Martelotto in Breast Cancer Res 16(3):210,
2014
). Tumor heterogeneity drives progression, metastasis and treatment resistance, but its relationship with tumor infiltrating immune cells is a matter of debate, where some argue that tumors with high heterogeneity may generate neoantigens that attract immune cells, and others claim that immune cells provide selection pressure that shapes tumor heterogeneity (McGranahan et al. in Science 351(6280):1463–1469,
2016
; McGranahan and Swanton in Cell 168(4):613–628,
2017
). We sought to study the association between tumor heterogeneity and immune cells in a real-world cohort utilizing The Cancer Genome Atlas.
Methods
Mutant allele tumor heterogeneity (MATH) was calculated to estimate intratumoral heterogeneity, and immune cell compositions were estimated using CIBERSORT. Survival analyses were demonstrated using Kaplan–Meir curves.
Results
Tumors with high heterogeneity (high MATH) were associated with worse overall survival (
p
= 0.049), as well as estrogen receptor-positive (
p
= 0.011) and non-triple-negative tumors (
p
= 0.01). High MATH tumors were also associated with less infiltration of anti-tumor CD8 (
p
< 0.013) and CD4 T cells (
p
< 0.00024), more tumor-promoting regulatory T cells (
p
< 4e−04), lower expression of T-cell exhaustion markers, specifically PDL-1 (
p
= 0.0031), IDO2 (
p
= 0.34), ADORA2A (
p
= 0.018), VISTA (
p
= 0.00013), and CCR4 (
p
< 0.00001), lower expression of cytolytic enzymes granzyme A (
p
= 0.0056) and perforin 1 (
p
= 0.053), and low cytolytic activity score (
p
= 0.0028).
Conclusions
High heterogeneity tumors are associated with less immune cell infiltration, less activation of the immune response, and worse survival in breast cancer. Our results support the notion that tumor heterogeneity is shaped by selection pressure of tumor-infiltrating immune cells.
Pancreatic ductal adenocarcinoma (PDAC) is known for its hypovascularity. Bevacizumab, an anti-angiogenic drug, added to standard chemotherapy demonstrated no improvement in outcome for PDAC. ...Therefore, we hypothesized that increased vascularity may be associated with improved outcomes in PDAC possibly due to better delivery of tumor specific immune cells. To test this hypothesis, PDAC patients were classified into either high or low CD31 expression groups utilizing mRNA expression from RNA-sequence data in The Cancer Genome Atlas (TCGA) pancreatic cancer cohort. High expression of CD31, which indicates presence of more vascular endothelial cells, was associated with significantly better OS (p = 0.002). Multivariate analysis demonstrated that residual tumor (R1, 2; p = 0.026) and CD31 low expression (p = 0.007) were the only independent predictors that negatively impacted OS. Vascular stability as well as immune response related pathways were significantly upregulated in the CD31 high expressing tumors. Furthermore, there were higher proportions of anti-cancer immune cells infiltration, including activated memory CD4+ T cells (p = 0.038), CD8+ T cells (p = 0.027), gamma-delta T cells (p < 0.001) as well as naïve B cells (p = 0.006), whereas lower proportions of regulatory T cell fractions (p = 0.009), which induce an immune tolerant microenvironment, in the CD31 high expressing tumors. These findings imply that stable vessels supply anti-cancer immune cells, which are at least partially responsible for better OS in the CD31 high expressing tumors. In conclusion, CD31 high expressing PDACs have better OS, which may be due to stable vessels that supply anti-cancer immune cells.
Partially paired data, either with incompleteness in one or both arms, are common in practice. For testing equality of means of two arms, practitioners often use only the portion of data with ...complete pairs and perform paired tests. Although such tests (referred as 'naive paired tests') are legitimate, their powers might be low as only partial data are utilized. The recently proposed 'P-value pooling methods', based on combining P-values from two tests, use all data, have reasonable type-I error control and good power property. While it is generally believed that 'P-value pooling methods' are superior to 'naive paired tests' in terms of power as the former use more data than the latter, no detailed power comparison has been done. This paper aims to compare powers of 'naive paired tests' and 'P-value pooling methods' analytically and our findings are counterintuitive, i.e. the 'P-value pooling methods' do not always outperform the naive paired tests in terms of power. Based on these results, we present guidance on how to select the best test for testing equality of means with partially paired data.
Due to the loss of DNA repair mechanisms in colorectal cancer (CRC) with microsatellite instability (MSI), somatic mutations accumulate within DNA; making them more prone to attack by tumor ...infiltrating lymphocytes (TIL) and macrophages. We hypothesize that MSI-High (MSI-H) patients have favorable survival due to increased tumor immunogenicity. The Cancer Genome Atlas (TCGA) was used to evaluate gene expression from 283 patients with CRC, comparing MSI-H and microsatellite stable (MSS) patients. CIBERSORT algorithm estimated the fraction of immune cell types. We found that low expression of DNA repair genes (MLH1, MLH3, PMS1, PMS2, ATR, PRKDC, ATM, BRCA2) associated with MSI-H. MSI-H was directly associated with Helper T-cells (p = 0.034) and M1 macrophages (p < 0.0001). MSI-H tumors associated with diminished intra-tumoral heterogeneity as well as higher expression of checkpoint molecules PD-1, PD-L1, CTLA4, LAG3 and TIM3 (p < 0.0001). Improved OS was seen in patients with low ATM, PMS2 and MLH3. In the TCGA CRC cohort, decreased expression of DNA repair genes associated with MSI-H. MSI-H patients had improved survival, likely due to higher TIL and M1 macrophage infiltration as well as lower intra-tumoral heterogeneity. MSI-H also associates with expression of immune checkpoint molecules with potential for development of therapeutic targets.
Some microRNAs (miRNAs) are known to suppress breast cancer. However, whether the expressions of these tumor suppressive miRNAs translate to patient survival were not investigated in large cohort. ...Nine miRNAs (miR-30a, miR-30c, miR-31, miR-126, miR-140, miR-146b, miR-200c, miR-206, and miR-335) known to be tumor suppressive miRNAs in breast cancer were investigated in Genomic Data Common data portal miRNA-Seq dataset and The Cancer Genome Atlas (TCGA) (n = 1052). Of the 9 miRNAs, miR-30a, miR-30c, miR-126, miR-140, miR-206, and miR-335 were found to have significantly lower expression in breast cancer tissues compared to paired normal breast tissue. High expression of miR-30a or miR-200c was associated with significantly better overall survival (OS). Gene Set Enrichment Analysis (GSEA) demonstrated that low expression levels of miR-30a had the tendency to associate with gene enrichment of EMT, while miR-200c did not, in TCGA cohort, and our findings support the need of validation using large cohort to use miRNA as prognostic biomarker for patients with breast cancer.
Abstract
Background
Blacks tend to have a stronger inflammatory immune response than Whites. We hypothesized that racial differences in host immunity also manifest in the tumor microenvironment, ...constituting part of a distinct aggressive tumor biology underlying higher mortality in Black women.
Methods
Pathological and gene expression profiling approaches were used for characterizing infiltrating immune cells in breast tumor microenvironment from 1315 patients from the Women’s Circle of Health Study. Racial differences in tumor immune phenotypes were compared, with results validated in a publicly accessible dataset. Prognostic associations of immune phenotypes were assessed in 3 independent cohorts.
Results
We found marked and consistent differences in tumor immune responses between Black and White patients. Not only did tumors from Blacks display a stronger overall immune presence but also the composition and quality of immune infiltrates differed, regardless of tumor subtypes. Black patients had a stronger CD4+ and B-cell response, and further, a more exhausted CD8+ T-cell profile. A signature indicating a higher ratio of exhausted CD8+ T cells to total CD8+ T cells (ExCD8-r) was consistently associated with poorer survival, particularly among hormone receptor–positive patients. Among hormone receptor–negative patients, combinations of the absolute fraction of CD8+ T cells and ExCD8-r signature identified the CD8lowExCD8-rhigh subgroup, the most prevalent among Blacks, with the worst survival.
Conclusions
Our findings of a distinct exhausted CD8+ T-cell signature in Black breast cancer patients indicate an immunobiological basis for their more aggressive disease and a rationale for the use of immune checkpoint inhibitors targeting the exhaustion phenotype.
Background
The use of biomarkers that allow early therapeutic intervention or intensive follow-up evaluation is expected to be a powerful means for reducing breast cancer mortality. MicroRNAs ...(miRNAs) are known to play major roles in cancer biology including metastasis. This study aimed to develop a novel miRNA risk score to predict patient survival and metastasis in breast cancer.
Methods
An integrated unbiased approach was applied to derive a composite risk score for prognosis based on miRNA expression in primary breast tumors in 1051 breast cancer patients from The Cancer Genome Atlas (TCGA). Further analysis of the risk score with metastasis/recurrence was performed using the TCGA data set and validated in a separate patient population using small RNA sequencing.
Results
The three-miRNAs risk score (miR-19a, miR-93, and miR-106a) was developed using the TCGA cohort, which predicted poor prognosis (
p
= 0.0005) independently of known clinical risk factors. The prognostic value was validated in another three following independent cohorts: GSE19536 (
p
= 0.0009), GSE22220 (
p
= 0.0003), and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) (
p
= 0.0023). The three-miRNAs risk score predicted bone recurrence in TCGA (
p
= 0.0052), and the findings were validated in another independent population of patients who experienced bone recurrence and age/stage-matched patients without any recurrence. The three-miRNAs risk score enriched multiple metastasis-related gene sets such as angiogenesis and epithelial mesenchymal transition in a gene-set-enrichment analysis.
Conclusions
The authors developed the novel miRNA-based risk score, which is a promising biomarker for prediction of worse survival and bone recurrence potential in breast cancer.
Estrogen receptor (ER) positive breast cancer (BC), the most abundant BC subtype, is notorious for poor response to neoadjuvant chemotherapy (NAC). The androgen receptor (AR) was reported to support ...estradiol-mediated ER activity in an in vitro system. Recently, ER-positive BC with fewer tumor infiltrating lymphocytes (TILs) was shown to have a better prognosis, opposite to the trend seen with ER-negative BC. We hypothesized that ER-positive BC with high expression of AR will have fewer TILs and an inferior response to NAC, but with a better prognosis. In both TCGA and METABRIC cohorts, AR expression was significantly higher in ER-positive BCs compared to ER-negatives (
< 0.001,
< 0.001, respectively) and it correlated with ER expression (
= 0.630,
= 0.509, respectively). In ER-positive tumors, AR high tumors enriched UV response down (NES = 2.01,
< 0.001), and AR low tumors enriched DNA repair (NES = -2.02,
< 0.001). AR high tumors were significantly associated with procancer regulatory T-cells, and AR low tumors were associated with anticancer immune cells, such as CD4, CD8, and Gamma-Delta T-cells and memory B-cells in ER-positive BC (
< 0.01). Further, cytolytic activity was significantly lower in AR high BC in both cohorts. Finally, AR high tumors had a significantly lower rate of attaining pathological complete response to NAC (GSE22358), but better survival. In conclusion, our results demonstrated that high AR has fewer tumor infiltrating lymphocytes as well as cytolytic activity and an inferior response to NAC, but better survival in ER-positive BC.
Sphingosine kinase 2 (SphK2) is known to phosphorylate the nuclear sphingolipid metabolite to generate sphingosine‐1‐phosphate (S1P). Nuclear S1P is involved in epigenetic regulation of gene ...expression; however, the underlying mechanisms are not well understood. In this work, we have identified the role of nuclear S1P and SphK2 in regulating hypoxia‐responsive master transcription factors hypoxia‐inducible factor (HIF)‐1α/2α, and their functions in breast cancer, with a focus on triple‐negative breast cancer (TNBC). We have shown SphK2 is associated with HIF‐1α in protein complexes, and is enriched at the promoters of HIF target genes, including vascular endothelial growth factor (VEGF), where it enhances local histone H3 acetylation and transcription. S1P specifically binds to the PAS domains of HIF‐1α. SphK2, and HIF‐1α expression levels are elevated in metastatic estrogen receptor‐positive (ER+) and TNBC clinical tissue specimens compared to healthy breast tissue samples. To determine if S1P formation in the nucleus by SphK2 is a key regulator of HIF functions, we found using a preclinical TNBC xenograft mouse model, and an existing selective SphK2 inhibitor K‐145, that nuclear S1P, histone acetylation, HIF‐1α expression, and TNBC tumor growth were all reduced in vivo. Our results suggest that S1P and SphK2 in the nucleus are linked to the regulation of HIF‐1α/2α functions associated with breast cancer progression, and may provide potential therapeutic targets.