Background
Circulating tumor DNA (ctDNA) has predictive and prognostic value in localized and metastatic cancer. This study analyzed the prognostic value of baseline and on-treatment ctDNA in ...metastatic gastroesophageal cancer (mGEC) using a region-specific next generation sequencing (NGS) panel.
Methods
Cell free DNA was isolated from plasma of patients before start of first-line palliative systemic treatment and after 9 and 18 weeks. Two NGS panels were designed comprising the most frequently mutated genes and targetable mutations in GEC. Tumor-derived mutations in matched metastatic biopsies were used to validate that the sequencing panels assessed true tumor-derived variants. Tumor volumes were calculated from baseline CT scans and correlated to variant allele frequency (VAF). Survival analyses were performed using univariable and multivariable Cox-regression analyses.
Results
ctDNA was detected in pretreatment plasma in 75% of 72 patients and correlated well with mutations in metastatic biopsies (86% accordance). The VAF correlated with baseline tumor volume (Pearson’s
R
0.53,
p
< 0.0001). Detection of multiple gene mutations at baseline in plasma was associated with worse overall survival (OS, HR 2.16, 95% CI 1.10–4.28;
p
= 0.027) and progression free survival (PFS, HR 2.71, 95% CI 1.28–5.73;
p
= 0.009). OS and PFS were inferior in patients with residual detectable ctDNA after 9 weeks of treatment (OS: HR 4.95, 95% CI 1.53–16.04;
p
= 0.008; PFS: HR 4.08, 95% CI 1.31–12.75;
p
= 0.016).
Conclusion
Based on our NGS panel, the number of ctDNA mutations before start of first-line chemotherapy has prognostic value. Moreover, residual ctDNA after three cycles of systemic treatment is associated with inferior survival.
Deconvolution of bulk gene expression profiles into the cellular components is pivotal to portraying tissue's complex cellular make-up, such as the tumor microenvironment. However, the inherently ...variable nature of gene expression requires a comprehensive statistical model and reliable prior knowledge of individual cell types that can be obtained from single-cell RNA sequencing. We introduce BLADE (Bayesian Log-normAl Deconvolution), a unified Bayesian framework to estimate both cellular composition and gene expression profiles for each cell type. Unlike previous comprehensive statistical approaches, BLADE can handle > 20 types of cells due to the efficient variational inference. Throughout an intensive evaluation with > 700 simulated and real datasets, BLADE demonstrated enhanced robustness against gene expression variability and better completeness than conventional methods, in particular, to reconstruct gene expression profiles of each cell type. In summary, BLADE is a powerful tool to unravel heterogeneous cellular activity in complex biological systems from standard bulk gene expression data.
Circulating tumor DNA (ctDNA) is assumed to reflect tumor burden and has been suggested as a tool for prognostication and follow‐up in patients with metastatic pancreatic ductal adenocarcinoma ...(mPDAC). However, the prognostic value of ctDNA and its relation with tumor burden has yet to be substantiated, especially in mPDAC. In this retrospective analysis of prospectively collected samples, cell‐free DNA from plasma samples of 58 treatment‐naive mPDAC patients was isolated and sequenced using a custom‐made pancreatobiliary NGS panel. Pathogenic mutations were detected in 26/58 (44.8%) samples. Cross‐check with droplet digital PCR showed good agreement in Bland–Altman analysis (p = 0.217, nonsignificance indicating good agreement). In patients with liver metastases, ctDNA was more frequently detected (24/37, p < 0.001). Tumor volume (3D reconstructions from imaging) and ctDNA variant allele frequency (VAF) were correlated (Spearman's ρ = 0.544, p < 0.001). Median overall survival (OS) was 3.2 (95% confidence interval CI 1.6–4.9) versus 8.4 (95% CI 1.6–15.1) months in patients with detectable versus undetectable ctDNA (p = 0.005). Both ctDNA VAF and tumor volume independently predicted OS after adjustment for carbohydrate antigen 19.9 and treatment regimen (hazard ratio HR 1.05, 95% CI 1.01–1.09, p = 0.005; HR 1.00, 95% CI 1.01–1.05, p = 0.003). In conclusion, our study showed that ctDNA detection rates are higher in patients with larger tumor volume and liver metastases. Nevertheless, measurements may diverge and, thus, can provide complementary information. Both ctDNA VAF and tumor volume were strong predictors of OS.
What's new?
Circulating tumor DNA (ctDNA) attracts much interest as a possible prognostic tool for cancer. Here, the authors showed that the quantity of ctDNA correlated strongly with tumor volume in metastatic pancreatic ductal adenocarcinoma (mPDAC). They conducted a retrospective analysis using samples collected from 58 untreated mPDAC patients. For this study, the authors designed a pancreatobiliary NGS panel, which they used to test the patients’ cell‐free DNA, along with droplet digital PCR. Both ctDNA variant allele frequency and tumor volume predicted overall survival, they found.
Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis of all common cancers. However, divergent outcomes exist between patients, suggesting distinct underlying tumor biology. Here, we ...delineated this heterogeneity, compared interconnectivity between classification systems, and experimentally addressed the tumor biology that drives poor outcome. RNA-sequencing of 90 resected specimens and unsupervised classification revealed four subgroups associated with distinct outcomes. The worst-prognosis subtype was characterized by mesenchymal gene signatures. Comparative (network) analysis showed high interconnectivity with previously identified classification schemes and high robustness of the mesenchymal subtype. From species-specific transcript analysis of matching patient-derived xenografts we constructed dedicated classifiers for experimental models. Detailed assessments of tumor growth in subtyped experimental models revealed that a highly invasive growth pattern of mesenchymal subtype tumor cells is responsible for its poor outcome. Concluding, by developing a classification system tailored to experimental models, we have uncovered subtype-specific biology that should be further explored to improve treatment of a group of PDAC patients that currently has little therapeutic benefit from surgical treatment.
Pancreatic ductal adenocarcinoma (PDAC) is a dismal disease with a poor clinical prognosis and unsatisfactory treatment options. We previously found that the transcription factor ...CCAAT/Enhancer-Binding Protein Delta (C/EBPδ) is lowly expressed in PDAC compared to healthy pancreas duct cells, and that patient survival and lymph node involvement in PDAC is correlated with the expression of C/EBPδ in primary tumor cells. C/EBPδ shares a homologous DNA-binding sequence with other C/EBP-proteins, leading to the presumption that other C/EBP-family members might act redundantly and compensate for the loss of C/EBPδ. This implies that patient stratification could be improved when expression levels of multiple C/EBP-family members are considered simultaneously. In this study, we assessed whether the quantification of C/EBPβ or C/EBPγ in addition to that of C/EBPδ might improve the prediction of patient survival and lymph node involvement using a cohort of 68 resectable PDAC patients. Using Kaplan-Meier analyses of patient groups with different C/EBP-expression levels, we found that both C/EBPβ and C/EBPγ can partially compensate for low C/EBPδ and improve patient survival. Further, we uncovered C/EBPβ as a novel predictor of a decreased likelihood of lymph node involvement in PDAC, and found that C/EBPβ and C/EBPδ can compensate for the lack of each other in order to reduce the risk of lymph node involvement. C/EBPγ, on the other hand, appears to promote lymph node involvement in the absence of C/EBPδ. Altogether, our results show that the redundancy of C/EBP-family members might have a profound influence on clinical prognoses and that the expression of both C/EPBβ and C/EBPγ should be taken into account when dichotomizing patients according to C/EBPδ expression.
Low grade serous ovarian carcinoma (LGSOC) demonstrates unique clinical and molecular features compared to other ovarian cancer types. The relationship between common histological features of LGSOC ...and molecular events, such as hormone receptor expression patterns and MAPK gene mutation status, remains poorly understood. Recent data suggest some of these molecular features may be biomarkers of response to recently introduced biologically-targeted therapies, namely endocrine therapy and MEK inhibitors. We utilize a cohort of 63 pathologically-confirmed LGSOC cases with whole exome sequencing and hormone receptor expression data to investigate these relationships. LGSOC cases demonstrated uniformly high oestrogen receptor (ER) expression, but variable progesterone receptor (PR) expression intensity. 60% and 37% of cases demonstrated micropapillary and macropapillary patterns of stromal invasion, respectively. 63% of cases demonstrated desmoplasia, which was significantly associated with advanced disease stage and visible residual disease after cytoreductive surgery. MAPK-mutant cases (KRAS, BRAF, NRAS) more frequently demonstrated macropapillary stromal invasion, while Chr1p loss was associated with desmoplasia and low PR expression. Presence of micropapillary stromal invasion and low PR expression were associated with significantly poorer survival after accounting for stage and residual disease status. Together, these data identify novel relationships between histopathological features and molecularly-defined subgroups in LGSOC.
Pancreatic ductal adenocarcinoma (PDAC) is characterized by an abundance of stroma. Multiple molecular classification efforts have identified a mesenchymal tumor subtype that is consistently ...characterized by high‐grade growth and poor clinical outcome. The relation between PDAC stroma and tumor subtypes is still unclear. Here, we aimed to identify how PDAC cells instruct the main cellular component of stroma, the pancreatic stellate cells (PSCs). We found in primary tissue that high‐grade PDAC had reduced collagen deposition compared to low‐grade PDAC. Xenografts and organotypic co‐cultures established from mesenchymal‐like PDAC cells featured reduced collagen and activated PSC content. Medium transfer experiments using a large set of PDAC cell lines revealed that mesenchymal‐like PDAC cells consistently downregulated ACTA2 and COL1A1 expression in PSCs and reduced proliferation. We identified colony‐stimulating factor 1 as the mesenchymal PDAC‐derived ligand that deactivates PSCs, and inhibition of its receptor CSF1R was able to counteract this effect. In conclusion, high‐grade PDAC features stroma that is low in collagen and activated PSC content, and targeting CSF1R offers direct options to maintain a tumor‐restricting microenvironment.
Synopsis
This study reveals that PDAC cells counteract mechanical restrictions of tumor stroma by secreting CSF‐1, which could be a potential therapeutic target.
Poorly differentiated, high‐grade pancreatic cancers, corresponding with the mesenchymal subtype, feature stroma that is low in collagen and activated pancreatic stellate cell content.
These mesenchymal subtype cancer cells secrete high amounts of CSF‐1.
CSF‐1 downregulates stromal activation markers and reduces proliferation via STAT3 in PSCs.
Depletion of CSF‐1 in high‐grade PDAC abolishes PSC deactivation.
High‐grade CSF‐1‐positive PDAC patient samples exhibit deactivated stroma.
This study reveals that PDAC cells counteract mechanical restrictions of tumor stroma by secreting CSF‐1, which could be a potential therapeutic target.
Pancreatic ductal adenocarcinoma (PDAC) is considered to be a poorly immunogenic cancer type that combines a low mutation burden with a strong immunosuppressive tumor microenvironment. Regulatory T ...cells (Tregs) are major drivers of immune suppression but their prognostic role, particularly in gastrointestinal malignancies, remains controversial. Lymphocytic infiltration in 122 PDAC samples was assessed by multispectral immunofluorescence with anti-Keratin, -CD3, -CD8, -FOXP3 and -CD163 antibodies. Differential infiltration by Tregs was analyzed in the context of transcriptomic profiles that were available for 65 tumors. High infiltration of CD3+CD8− (mainly CD4+) T cells and, especially, of the subset expressing FOXP3 (Tregs) was associated with improved patient survival, whilst cytotoxic CD3+CD8+ T cell infiltration did not have an impact on overall survival. Transcriptomic analysis revealed three signatures in PDAC tumors comprising of epithelial-mesenchymal transition (EMT)/stromal, metabolic, and secretory/pancreatic signature. However, none of these signatures explained differences in Treg infiltration. We show that Tregs associate with improved overall survival in PDAC patients. This effect was independent of cytotoxic T cell infiltration and the transcriptomic profiles of their respective tumors. These findings provide a new layer of complexity in the study of PDAC tumor microenvironment that must be considered when developing immunotherapeutic interventions for this disease.
Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with a limited number of known driver mutations but considerable cancer cell heterogeneity. Phosphoproteomics provides a direct ...read‐out of aberrant signaling and the resultant clinically relevant phenotype. Mass spectrometry (MS)‐based proteomics and phosphoproteomics were applied to 42 PDAC tumors. Data encompassed over 19 936 phosphoserine or phosphothreonine (pS/T; in 5412 phosphoproteins) and 1208 phosphotyrosine (pY; in 501 phosphoproteins) sites and a total of 3756 proteins. Proteome data identified three distinct subtypes with tumor intrinsic and stromal features. Subsequently, three phospho‐subtypes were apparent: two tumor intrinsic (Phos1/2) and one stromal (Phos3), resembling known PDAC molecular subtypes. Kinase activity was analyzed by the Integrative iNferred Kinase Activity (INKA) scoring. Phospho‐subtypes displayed differential phosphorylation signals and kinase activity, such as FGR and GSK3 activation in Phos1, SRC kinase family and EPHA2 in Phos2, and EGFR, INSR, MET, ABL1, HIPK1, JAK, and PRKCD in Phos3. Kinase activity analysis of an external PDAC cohort supported our findings and underscored the importance of PI3K/AKT and ERK pathways, among others. Interestingly, unfavorable patient prognosis correlated with higher RTK, PAK2, STK10, and CDK7 activity and high proliferation, whereas long survival was associated with MYLK and PTK6 activity, which was previously unknown. Subtype‐associated activity profiles can guide therapeutic combination approaches in tumor and stroma‐enriched tissues, and emphasize the critical role of parallel signaling pathways. In addition, kinase activity profiling identifies potential disease markers with prognostic significance.
Despite the treatment advancements in pancreatic ductal adenocarcinoma, 5‐year survival rates are still below 10%. Analyzing the phosphoproteome of 42 tumors, three subtypes with intrinsic and stromal features were revealed, each displaying kinase activity patterns with potential treatment options. Kinase activities associated with specific mutations, and correlated with clinical prognosis. This study provides insights for tailored and combination strategies.
Irreversible electroporation causes cell death through low frequency, high voltage electrical pulses and is increasingly used to treat non-resectable cancers. A recent systematic review revealed that ...tissue damage through irreversible electroporation is time-dependent, but the impact of time on the ablation zone size remains unknown. Irreversible electroporation ablations were performed hourly during 24 consecutive hours in the peripheral liver of 2 anaesthetized domestic pigs using clinical treatment settings. Immediately after the 24th ablation, the livers were harvested and examined for tissue response in time based on macroscopic and microscopic pathology. The impact of time on these outcomes was assessed with Spearman rank correlation test. Ablation zones were sharply demarcated as early as 1 hour after treatment. During 24 hours, the ablation zones showed a significant increase in diameter (rs = 0.493, P = .014) and total surface (rs = 0.499, P = .013), whereas the impact of time on the homogeneous ablated area was not significant (rs = 0.172, P = .421). Therefore, the increase in size could mainly be attributed to an increase in the transition zone. Microscopically, the ablation zones showed progression in cell death and inflammation. This study assessed the dynamics of irreversible electroporation on the porcine liver during 24 consecutive hours and found that the pathological response (ie, cell death/inflammation), and ablation size continue to develop for at least 24 hours. Consequently, future studies on irreversible electroporation should prolong their observation period.