Trastuzumab and first-line taxane chemotherapy in metastatic breast cancer patients with a HER2-negative tumor and HER2-positive circulating tumor cells: a phase II trial

PURPOSE

HER2 overexpressing circulating tumor cells (CTCs) are observed in up to 25% of HER2-negative metastatic breast cancer patients. Since targeted anti-HER2 therapy has drastically improved clinical outcomes of patients with HER2-positive breast cancer, we hypothesized that patients with HER2 overexpressing CTCs might benefit from the addition of trastuzumab to chemotherapy.

METHODS

In this single-arm, phase II trial, patients with HER2-positive CTCs received trastuzumab as addition to first-line treatment with taxane chemotherapy. Patients with detectable CTCs but without HER2 overexpression that received taxane chemotherapy only, were used as control group. The primary outcome measure was progression-free rate at 6 months (PFR6), with a target of 80%. In November 2022, the study was terminated early due to slow patient accrual.

RESULTS

63 patients were screened, of which eight patients had HER2-positive CTCs and were treated with trastuzumab. The median number of CTCs was 15 per 7.5 ml of blood (range 1-131) in patients with HER2-positive CTCs, compared to median 5 (range 1-1047) in the control group. PFR6 was 50% in the trastuzumab group and 54% in the taxane monotherapy group, with no significant difference in median PFS (8 versus 9 months, p = 0.51).

CONCLUSION

No clinical benefit of trastuzumab was observed, although this study was performed in a limited number of patients. Additionally, we observed a strong correlation between the number of evaluable CTCs and the presence of HER2-positive CTCs. We argue that randomized studies investigating agents that are proven to be solely effective in the HER2-positive patient group in patients with HER2-positive CTCs and HER2-negative tissue are currently infeasible. Several factors contribute to this impracticality, including the need for more stringent thresholds, and the rapidly evolving landscape of cancer treatments.

Whole-genome mapping of APOBEC mutagenesis in metastatic urothelial carcinoma identifies driver hotspot mutations and a novel mutational signature

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) enzymes mutate specific DNA sequences and hairpin-loop structures, challenging the distinction between passenger and driver hotspot mutations. Here, we characterized 115 whole genomes of metastatic urothelial carcinoma (mUC) to identify APOBEC mutagenic hotspot drivers. APOBEC-associated mutations were detected in 92% of mUCs and were equally distributed across the genome, while APOBEC hotspot mutations (ApoHMs) were enriched in open chromatin. Hairpin loops were frequent targets of didymi (twins in Greek), two hotspot mutations characterized by the APOBEC SBS2 signature, in conjunction with an uncharacterized mutational context (Ap[C>T]). Next, we developed a statistical framework that identified ApoHMs as drivers in coding and non-coding genomic regions of mUCs. Our results and statistical framework were validated in independent cohorts of 23 non-metastatic UCs and 3,744 samples of 17 metastatic cancers, identifying cancer-type-specific drivers. Our study highlights the role of APOBEC in cancer development and may contribute to developing novel targeted therapy options for APOBEC-driven cancers.

Tumor-Infiltrating Lymphocytes in Triple-Negative Breast Cancer

Importance

The association of tumor-infiltrating lymphocyte (TIL) abundance in breast cancer tissue with cancer recurrence and death in patients with early-stage triple-negative breast cancer (TNBC) who are not treated with adjuvant or neoadjuvant chemotherapy is unclear.

Objective

To study the association of TIL abundance in breast cancer tissue with survival among patients with early-stage TNBC who were treated with locoregional therapy but no chemotherapy.

Design, Setting, and Participants

Retrospective pooled analysis of individual patient-level data from 13 participating centers in North America (Rochester, Minnesota; Vancouver, British Columbia, Canada), Europe (Paris, Lyon, and Villejuif, France; Amsterdam and Rotterdam, the Netherlands; Milan, Padova, and Genova, Italy; Gothenburg, Sweden), and Asia (Tokyo, Japan; Seoul, Korea), including 1966 participants diagnosed with TNBC between 1979 and 2017 (with follow-up until September 27, 2021) who received treatment with surgery with or without radiotherapy but no adjuvant or neoadjuvant chemotherapy.

Exposure

TIL abundance in breast tissue from resected primary tumors.

Main Outcomes and Measures

The primary outcome was invasive disease-free survival [iDFS]. Secondary outcomes were recurrence-free survival [RFS], survival free of distant recurrence [distant RFS, DRFS], and overall survival. Associations were assessed using a multivariable Cox model stratified by participating center.

Results

This study included 1966 patients with TNBC (median age, 56 years [IQR, 39-71]; 55% had stage I TNBC). The median TIL level was 15% (IQR, 5%-40%). Four-hundred seventeen (21%) had a TIL level of 50% or more (median age, 41 years [IQR, 36-63]), and 1300 (66%) had a TIL level of less than 30% (median age, 59 years [IQR, 41-72]). Five-year DRFS for stage I TNBC was 94% (95% CI, 91%-96%) for patients with a TIL level of 50% or more, compared with 78% (95% CI, 75%-80%) for those with a TIL level of less than 30%; 5-year overall survival was 95% (95% CI, 92%-97%) for patients with a TIL level of 50% or more, compared with 82% (95% CI, 79%-84%) for those with a TIL level of less than 30%. At a median follow-up of 18 years, and after adjusting for age, tumor size, nodal status, histological grade, and receipt of radiotherapy, each 10% higher TIL increment was associated independently with improved iDFS (hazard ratio [HR], 0.92 [0.89-0.94]), RFS (HR, 0.90 [0.87-0.92]), DRFS (HR, 0.87 [0.84-0.90]), and overall survival (0.88 [0.85-0.91]) (likelihood ratio test, P < 10e-6).

Conclusions and Relevance

In patients with early-stage TNBC who did not undergo adjuvant or neoadjuvant chemotherapy, breast cancer tissue with a higher abundance of TIL levels was associated with significantly better survival. These results suggest that breast tissue TIL abundance is a prognostic factor for patients with early-stage TNBC.

The effect of (neo)adjuvant chemotherapy on long-term survival outcomes in patients with invasive lobular breast cancer treated with endocrine therapy: A retrospective cohort study

BACKGROUND

Despite histological and molecular differences between invasive lobular carcinoma (ILC) and invasive carcinoma of no special type, according to national treatment guidelines no distinction is made regarding the use of (neo)adjuvant chemotherapy. Studies on the long-term outcome of chemotherapy in patients with ILC are scarce and show inconclusive results.

METHODS

All patients with estrogen receptor (ER)-positive, human epidermal growth factor receptor 2 (HER2)-negative ILC with an indication for chemotherapy treated with adjuvant endocrine therapy were selected from the Erasmus Medical Center Breast Cancer database. Cox proportional hazards models were used to estimate the effect of chemotherapy on recurrence-free survival (RFS), breast cancer-specific survival (BCSS), and overall survival (OS).

RESULTS

A total of 520 patients were selected, of whom 379 were treated with chemotherapy and 141 were not. Patients in the chemotherapy group were younger (51 vs. 61 years old; p < .001), had a higher T status (T3+, 33% vs. 14%; p < .001), and more often had lymph node involvement (80% vs. 49%; p < .001) in comparison to the no-chemotherapy group. After adjusting for confounders, chemotherapy treatment was not associated with better RFS (hazard ratio [HR], 1.20; 95% confidence interval [CI], 0.63-2.31), BCSS (HR, 1.24; 95% CI, 0.60-2.58), or OS (HR, 0.97; 95% CI, 0.56-1.66). This was also reflected by adjusted Cox survival curves in the chemotherapy versus no-chemotherapy group for RFS (75% vs. 79%), BCSS (80% vs. 84%), and OS (72% vs. 71%).

CONCLUSIONS

Chemotherapy is not associated with improved RFS, BCSS, or OS for patients with ER+/HER2- ILC treated with adjuvant endocrine therapy and with an indication for chemotherapy.

Functional Homologous Recombination (HR) Screening Shows the Majority of BRCA1/2-Mutant Breast and Ovarian Cancer Cell Lines Are HR-Proficient

Tumors with a pathogenic BRCA1/2 mutation are homologous recombination (HR)-deficient (HRD) and consequently sensitive to platinum-based chemotherapy and Poly-[ADP-Ribose]-Polymerase inhibitors (PARPi). We hypothesized that functional HR status better reflects real-time HR status than BRCA1/2 mutation status. Therefore, we determined the functional HR status of 53 breast cancer (BC) and 38 ovarian cancer (OC) cell lines by measuring the formation of RAD51 foci after irradiation. Discrepancies between functional HR and BRCA1/2 mutation status were investigated using exome sequencing, methylation and gene expression data from 50 HR-related genes. A pathogenic BRCA1/2 mutation was found in 10/53 (18.9%) of BC and 7/38 (18.4%) of OC cell lines. Among BRCA1/2-mutant cell lines, 14/17 (82.4%) were HR-proficient (HRP), while 1/74 (1.4%) wild-type cell lines was HRD. For most (80%) cell lines, we explained the discrepancy between functional HR and BRCA1/2 mutation status. Importantly, 12/14 (85.7%) BRCA1/2-mutant HRP cell lines were explained by mechanisms directly acting on BRCA1/2. Finally, functional HR status was strongly associated with COSMIC single base substitution signature 3, but not BRCA1/2 mutation status. Thus, the majority of BRCA1/2-mutant cell lines do not represent a suitable model for HRD. Moreover, exclusively determining BRCA1/2 mutation status may not suffice for platinum-based chemotherapy or PARPi patient selection.

Expression and Localization of Ferritin-Heavy Chain Predicts Recurrence for Breast Cancer Patients with a BRCA1/2 Mutation

The ferritin-heavy chain (FTH1) is the catalytic subunit of the ferroxidase ferritin, which prevents oxidative DNA damage via intracellular iron storage. FTH1 was shown to be a prognostic marker for triple-negative breast cancer (BC) patients and associated with an enrichment of CD8+ effector T cells. However, whether the expression and localization of FTH1 are also associated with clinical outcome in other BC subtypes is unknown. Here, we investigated the association of FTH1 with time to survival in BCs from 222 BRCA1/2 mutation carriers by immunohistochemistry on tissue microarrays. In addition, for 51 of these patients, the association between FTH1 and specific subsets of T cells was evaluated on whole slides using automatic scoring algorithms. We revealed that nuclear FTH1 (nFTH1) expression, in multivariable analyses, was associated with a shorter disease-free (HR = 2.71, 95% CI = 1.49-4.92, p = 0.001) and metastasis-free survival (HR = 3.54, 95% CI = 1.45-8.66, p = 0.006) in patients carrying a BRCA1/2 mutation. However, we found no relation between cytoplasmic FTH1 expression and survival of BRCA1/2 mutation carriers. Moreover, we did not detect an association between FTH1 expression and the amount of CD45+ (p = 0.13), CD8+ (p = 0.18), CD4+ (p = 0.20) or FOXP3+ cells (p = 0.17). Consequently, the mechanism underlying the worse recurrence-free survival of nFTH1 expression in BRCA1/2 mutation carriers needs further investigation.

Integrative whole-genome and transcriptome analysis of HER2-amplified metastatic breast cancer

BACKGROUND

In breast cancer, the advent of anti-HER2 therapies has made HER2+ tumors a highly relevant subgroup. However, the exact characteristics which prohibit clinical response to anti-HER2 therapies and drive disease progression are not yet fully known. Integrative whole-genome and transcriptomic sequencing data from both primary and metastatic HER2-positive breast cancer will enhance our understanding of underlying biological processes.

METHODS

Here, we used WGS and RNA sequencing data of 700 metastatic breast tumors, of which 68 being HER2+, to search for specific genomic features of HER2+ disease and therapy resistance. Furthermore, we integrated results with transcriptomic data to associate tumors exhibiting a HER2+-specific gene expression profile with ERBB2 mutation status, prior therapy and relevant gene expression signatures.

RESULTS

Overall genomic profiles of primary and metastatic HER2+ breast cancers were similar, and no specific acquired genomics traits connected to prior anti-HER2 treatment were observed. However, specific genomic features were predictive of progression-free survival on post-biopsy anti-HER2 treatment. Furthermore, a HER2-driven expression profile grouped HER2-amplified tumors with ERBB2-mutated cases and cases without HER2 alterations. The latter were reported as ER positive in primary disease, but the metastatic biopsy showed low ESR1 expression and upregulation of the MAPK pathway, suggesting transformation to ER independence.

CONCLUSIONS

In summary, although the quantity of variants increased throughout HER2-positive breast cancer progression, the genomic composition remained largely consistent, thus yielding no new major processes beside those already operational in primary disease. Our results suggest that integrated genomic and transcriptomic analyses may be key in establishing therapeutic options.

Improved enrichment of circulating tumor cells from diagnostic leukapheresis product

The median number of circulating tumor cells (CTCs) detected in 7.5 mL of peripheral blood by CellSearch (PB-CS) in patients with metastatic prostate cancer is in the order of 1-10, which means many samples have insufficient tumor cells for comprehensive characterization. A significant increase is obtained through diagnostic leukapheresis (DLA), however, only 2%-3% of the DLA product can be processed per CellSearch test, limiting the gain. We processed aliquots from 30 DLA products of metastatic prostate cancer patients consisting of 0.2 × 109 leukocytes using CellSearch (DLA-CS) as well as the newly introduced reduced enrichment reagent protocol (RER), which uses 10-fold less enrichment reagents than DLA-CS. The number of tumor cells and the total number of captured cells were determined using the CellTracks Analyzer. Additionally, for six DLA samples, a 1.0 × 109 leukocyte aliquot was processed (RER+), using twofold less enrichment reagents than DLA-CS. A median 2.7-fold reduction in leukocyte co-enrichment was found between DLA-CS and RER methods without any loss in tumor cell recovery (Wilcoxon Signed Ranks Test, p = 0.953). Using 1.0 × 109 leukocyte aliquots a fourfold increase in tumor cells was found compared to DLA-CS and a 19-fold increase compared to PB-CS was obtained. The here-introduced RER protocol results in a higher final sample purity without any loss in tumor cell recovery while using 10-fold less CellSearch capture reagent. With this improved method, 26% of the leukapheresis sample can now be processed using reagents from a single CellSearch test, enabling the obtainment of a sufficient number of CTCs for comprehensive characterization in most metastatic prostate cancer patients.

Circulating Tumour DNA as Biomarker for Colorectal Liver Metastases: A Systematic Review and Meta-Analysis

Circulating tumour DNA (ctDNA) is a potential biomarker that could contribute to more judicious patient selection for personalised treatment. This review and meta-analysis gives an overview of the current knowledge in the literature investigating the value of ctDNA in patients with colorectal liver metastases (CRLM). A systematic search was conducted in electronic databases for studies published prior to the 26th of May 2023. Studies investigating the association between ctDNA and oncological outcomes in patients undergoing curative-intent local therapy for CRLM were included. Meta-analyses were performed to pool hazard ratios (HR) for the recurrence-free survival (RFS) and overall survival (OS). A total of eleven studies were included and nine were eligible for meta-analyses. Patients with detectable ctDNA after surgery experienced a significantly higher chance of recurrence (HR 3.12, 95% CI 2.27-4.28, p < 0.000010) and shorter OS (HR 5.04, 95% CI 2.53-10.04, p < 0.00001) compared to patients without detectable ctDNA. A similar association for recurrence was found in patients with detectable ctDNA after the completion of adjuvant therapy (HR 6.39, 95% CI 2.13-19.17, p < 0.0009). The meta-analyses revealed no association between detectable ctDNA before surgery and the RFS and OS. These meta-analyses demonstrate the strong association between detectable ctDNA after treatment and oncological outcomes in CRLM patients.

Liquid Biopsies for Colorectal Cancer and Advanced Adenoma Screening and Surveillance: What to Measure?

Colorectal cancer (CRC) colonoscopic surveillance is effective but burdensome. Circulating tumor DNA (ctDNA) analysis has emerged as a promising, minimally invasive tool for disease detection and management. Here, we assessed which ctDNA assay might be most suitable for a ctDNA-based CRC screening/surveillance blood test. In this prospective, proof-of-concept study, patients with colonoscopies for Lynch surveillance or the National Colorectal Cancer screening program were included between 7 July 2019 and 3 June 2022. Blood was drawn, and if advanced neoplasia (adenoma with villous component, high-grade dysplasia, ≥10 mm, or CRC) was detected, it was analyzed for chromosomal copy number variations, single nucleotide variants, and genome-wide methylation (MeD-seq). Outcomes were compared with corresponding patients' tissues and the MeD-seq results of healthy blood donors. Two Lynch carriers and eight screening program patients were included: five with CRC and five with advanced adenomas. cfDNA showed copy number variations and single nucleotide variants in one patient with CRC and liver metastases. Eight patients analyzed with MeD-seq showed clustering of Lynch-associated and sporadic microsatellite instable lesions separate from microsatellite stable lesions, as did healthy blood donors. In conclusion, whereas copy number changes and single nucleotide variants were only detected in one patient, cfDNA methylation profiles could discriminate all microsatellite instable advanced neoplasia, rendering this tool particularly promising for LS surveillance. Larger studies are warranted to validate these findings.

Genomic Alterations Associated with Estrogen Receptor Pathway Activity in Metastatic Breast Cancer Have a Differential Impact on Downstream ER Signaling

Mutations in the estrogen receptor gene (ESR1), its transcriptional regulators, and the mitogen-activated protein kinase (MAPK) pathway are enriched in patients with endocrine-resistant metastatic breast cancer (MBC). Here, we integrated whole genome sequencing with RNA sequencing data from the same samples of 101 ER-positive/HER2-negative MBC patients who underwent a tumor biopsy prior to the start of a new line of treatment for MBC (CPCT-02 study, NCT01855477) to analyze the downstream effects of DNA alterations previously linked to endocrine resistance, thereby gaining a better understanding of the associated mechanisms. Hierarchical clustering was performed using expression of ESR1 target genes. Genomic alterations at the DNA level, gene expression levels, and last administered therapy were compared between the identified clusters. Hierarchical clustering revealed two distinct clusters, one of which was characterized by increased expression of ESR1 and its target genes. Samples in this cluster were significantly enriched for mutations in ESR1 and amplifications in FGFR1 and TSPYL. Patients in the other cluster showed relatively lower expression levels of ESR1 and its target genes, comparable to ER-negative samples, and more often received endocrine therapy as their last treatment before biopsy. Genes in the MAPK-pathway, including NF1, and ESR1 transcriptional regulators were evenly distributed. In conclusion, RNA sequencing identified a subgroup of patients with clear expression of ESR1 and its downstream targets, probably still benefiting from ER-targeting agents. The lower ER expression in the other subgroup might be partially explained by ER activity still being blocked by recently administered endocrine treatment, indicating that biopsy timing relative to endocrine treatment needs to be considered when interpreting transcriptomic data.

mFast-SeqS-based aneuploidy score in circulating cell-free DNA is a prognostic biomarker in prostate cancer

Multiple prognostic biomarkers, including circulating tumour cell (CTC) counts, exist in metastatic castration-resistant prostate cancer (mCRPC) patients, but none of them have been implemented into daily clinical care. The modified fast aneuploidy screening test-sequencing system (mFast-SeqS), which yields a genome-wide aneuploidy score, is able to reflect the fraction of cell-free tumour DNA (ctDNA) within cell-free DNA (cfDNA) and may be a promising biomarker in mCRPC. In this study, we investigated the prognostic value of dichotomized aneuploidy scores (< 5 vs. ≥ 5) as well as CTC counts (< 5 vs. ≥ 5) in 131 mCRPC patients prior to treatment with cabazitaxel. We validated our findings in an independent cohort of 50 similarly treated mCRPC patients. We observed that, similar to the dichotomized CTC count [HR: 2.92; 95% confidence interval (CI);1.84-4.62], dichotomized aneuploidy scores (HR: 3.24; CI: 2.12-4.94) significantly correlated with overall survival in mCRPC patients. We conclude that a dichotomized aneuploidy score from cfDNA is a prognostic marker for survival in mCRPC patients within our discovery cohort and in an independent mCRPC validation cohort. Therefore, this easy and robust minimally-invasive assay can be readily implemented as a prognostic marker in mCRPC. A dichotomized aneuploidy score might also be used as a stratification factor in clinical studies to account for tumour load.

Liquid Biopsies for Early Response Evaluation of Radium-223 in Metastatic Prostate Cancer

PURPOSE

Reliable biomarkers for response monitoring during radium-223 treatment in patients with metastatic castration-resistant prostate cancer (mCRPC) are lacking. Circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA), obtained from liquid biopsies, are shown to have prognostic value in mCRPC. The aim of this study was to determine the value of CTCs and ctDNA for response evaluation of radium-223.

METHODS

In this prospective multicenter study, longitudinal blood draws and imaging were performed in 28 patients with mCRPC and predominantly bone disease, who were treated with radium-223. CTCs were counted (CELLSEARCH CTC test), while fraction of ctDNA was estimated by measuring aneuploidy of cell-free DNA (cfDNA; modified Fast Aneuploidy Screening Test-Sequencing System). CTC counts and aneuploidy score (AS) were categorized as low (<5) and high (≥5). Primary and secondary clinical end points were failure-free survival (FFS), and overall survival (OS) and development of extraosseous metastases, respectively. Additionally, CTC count and AS were related to alkaline phosphatase (ALP) and total tumor volume in bone (TTVbone) on positron emission tomography-computed tomography with 68gallium prostate-specific membrane antigen.

RESULTS

FFS was longer in patients with a low CTC count or AS either at baseline or after 12 weeks, whereas for OS, only a significant association with CTC count was observed. Liquid biopsy results correlated well with ALP and TTVbone at baseline, but not with change in both parameters after three cycles of radium-223. AS and CTC count were significantly correlated.

CONCLUSION

CTC count and AS of cfDNA at baseline and during treatment predict clinical response to radium-223 in patients with mCRPC, warranting future evaluation of their value in treatment guidance.

Obesity-associated changes in molecular biology of primary breast cancer

Obesity is associated with an increased risk of developing breast cancer (BC) and worse prognosis in BC patients, yet its impact on BC biology remains understudied in humans. This study investigates how the biology of untreated primary BC differs according to patients' body mass index (BMI) using data from >2,000 patients. We identify several genomic alterations that are differentially prevalent in overweight or obese patients compared to lean patients. We report evidence supporting an ageing accelerating effect of obesity at the genetic level. We show that BMI-associated differences in bulk transcriptomic profile are subtle, while single cell profiling allows detection of more pronounced changes in different cell compartments. These analyses further reveal an elevated and unresolved inflammation of the BC tumor microenvironment associated with obesity, with distinct characteristics contingent on the estrogen receptor status. Collectively, our analyses imply that obesity is associated with an inflammaging-like phenotype. We conclude that patient adiposity may play a significant role in the heterogeneity of BC and should be considered for BC treatment tailoring.

Tumor-agnostic ctDNA levels by mFAST-SeqS in first-line HR-positive, HER2 negative metastatic breast cancer patients as a biomarker for survival

This prospective cohort study reports aneuploidy score by mFast-SeqS as a strong prognostic marker in MBC patients. mFAST-SeqS is an affordable and easily implementable method for the assessment of total ctDNA levels and, as such, provides an alternative prognostic tool. One mixed cohort (cohort A, n = 45) starting any type of treatment in any line of therapy and one larger cohort (cohort B, n = 129) consisting of patients starting aromatase inhibitors (AI) as first-line therapy were used. mFAST-SeqS was performed using plasma of blood in which CTCs (CellSearch) were enumerated. The resulting aneuploidy score was correlated with categorized CTC count and associated with outcome. The aneuploidy score was significantly correlated with CTC count, but discordance was observed in 31.6% when applying cut-offs of 5. In both cohorts, aneuploidy score was a significant prognostic marker for both PFS and OS. In the Cox regression models, the HR for aneuploidy score for PFS was 2.52 (95% CI: 1.56-4.07), and the HR for OS was 2.37 (95% CI: 1.36-4.14). Results presented here warrant further investigations into the clinical utility of this marker in MBC patients.

Machine learning-based somatic variant calling in cell-free DNA of metastatic breast cancer patients using large NGS panels

Next generation sequencing of cell-free DNA (cfDNA) is a promising method for treatment monitoring and therapy selection in metastatic breast cancer (MBC). However, distinguishing tumor-specific variants from sequencing artefacts and germline variation with low false discovery rate is challenging when using large targeted sequencing panels covering many tumor suppressor genes. To address this, we built a machine learning model to remove false positive variant calls and augmented it with additional filters to ensure selection of tumor-derived variants. We used cfDNA of 70 MBC patients profiled with both the small targeted Oncomine breast panel (Thermofisher) and the much larger Qiaseq Human Breast Cancer Panel (Qiagen). The model was trained on the panels' common regions using Oncomine hotspot mutations as ground truth. Applied to Qiaseq data, it achieved 35% sensitivity and 36% precision, outperforming basic filtering. For 20 patients we used germline DNA to filter for somatic variants and obtained 245 variants in total, while our model found seven variants, of which six were also detected using the germline strategy. In ten tumor-free individuals, our method detected in total one (potentially germline) variant, in contrast to 521 variants detected without our model. These results indicate that our model largely detects somatic variants.

Plasma ESR1 mutations and outcome to first-line paclitaxel and bevacizumab in patients with advanced ER-positive/HER2-negative breast cancer

BACKGROUND

ESR1 mutations have been identified as mechanism for endocrine resistance and are also associated with a decreased overall survival. We assessed ESR1 mutations in circulating tumor DNA (ctDNA) for impact on outcome to taxane-based chemotherapy in advanced breast cancer patients.

METHODS

ESR1 mutations were determined in archived plasma samples from patients treated with paclitaxel and bevacizumab (AT arm, N = 91) in the randomized phase II ATX study. Samples collected at baseline (n = 51) and at cycle 2 (n = 13, C2) were analyzed using a breast cancer next-generation sequencing panel. This study was powered to detect a benefit in progression-free survival (PFS) at six months for patients treated with paclitaxel/bevacizumab compared to historical trials with fulvestrant. PFS, overall survival (OS), and ctDNA dynamics were exploratory analyses.

RESULTS

PFS at six months was 86% (18/21) in patients with an ESR1 mutation detected and 85% (23/27) in wildtype ESR1 patients. In our exploratory analysis, median progression-free survival (PFS) was 8.2 months [95% CI, 7.6-8.8] for ESR1 mutant patients versus 8.7 months [95% confidence interval (CI), 8.3-9.2] for ESR1 wildtype patients [p = 0.47]. The median overall survival (OS) was 20.7 months [95% CI, 6.6-33.7] for ESR1 mutant patients versus 28.1 months [95% confidence interval (CI), 19.3-36.9] for ESR1 wildtype patients [p = 0.27]. Patients with ≥ two ESR1 mutations had a significantly worse OS, but not PFS, compared to those who did not [p = 0.003]. Change in ctDNA level at C2 was not different between ESR1 and other mutations.

CONCLUSIONS

Presence of ESR1 mutations in baseline ctDNA might not be associated with inferior PFS and OS in advanced breast cancer patients treated with paclitaxel/bevacizumab.

Breast cancer genomes from CHEK2 c.1100delC mutation carriers lack somatic TP53 mutations and display a unique structural variant size distribution profile

BACKGROUND

CHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC.

METHODS

We sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases.

RESULTS

BC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile.

CONCLUSIONS

CHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance.

Aberrant APOBEC3B Expression in Breast Cancer Is Linked to Proliferation and Cell Cycle Phase

APOBEC3B (A3B) is aberrantly overexpressed in a subset of breast cancers, where it associates with advanced disease, poor prognosis, and treatment resistance, yet the causes of A3B dysregulation in breast cancer remain unclear. Here, A3B mRNA and protein expression levels were quantified in different cell lines and breast tumors and related to cell cycle markers using RT-qPCR and multiplex immunofluorescence imaging. The inducibility of A3B expression during the cell cycle was additionally addressed after cell cycle synchronization with multiple methods. First, we found that A3B protein levels within cell lines and tumors are heterogeneous and associate strongly with the proliferation marker Cyclin B1 characteristic of the G₂/M phase of the cell cycle. Second, in multiple breast cancer cell lines with high A3B, expression levels were observed to oscillate throughout the cell cycle and again associate with Cyclin B1. Third, induction of A3B expression is potently repressed throughout G₀/early G₁, likely by RB/E2F pathway effector proteins. Fourth, in cells with low A3B, induction of A3B through the PKC/ncNF-κB pathway occurs predominantly in actively proliferating cells and is largely absent in cells arrested in G₀. Altogether, these results support a model in which dysregulated A3B overexpression in breast cancer is the cumulative result of proliferation-associated relief from repression with concomitant pathway activation during the G₂/M phase of the cell cycle.

A Compendium of AR Splice Variants in Metastatic Castration-Resistant Prostate Cancer

Treatment-induced AR alterations, including AR alternative splice variants (AR-Vs), have been extensively linked to harboring roles in primary and acquired resistance to conventional and next-generation hormonal therapies in prostate cancer and therefore have gained momentum. Our aim was to uniformly determine recurrent AR-Vs in metastatic castration-resistant prostate cancer (mCRPC) using whole transcriptome sequencing in order to assess which AR-Vs might hold potential diagnostic or prognostic relevance in future research. This study reports that in addition to the promising AR-V7 as a biomarker, AR45 and AR-V3 were also seen as recurrent AR-Vs and that the presence of any AR-V could be associated with higher AR expression. With future research, these AR-Vs may therefore harbor similar or complementary roles to AR-V7 as predictive and prognostic biomarkers in mCRPC or as proxies for abundant AR expression.

Methylated Cell-Free DNA Sequencing (MeD-seq) of LpnPI Digested Fragments to Identify Early Progression in Metastatic Renal Cell Carcinoma Patients on Watchful Waiting

According to the current guidelines, watchful waiting (WW) is a feasible option for patients with good or intermediate prognosis renal-cell carcinoma (RCC). However, some patients rapidly progress during WW, requiring the initiation of treatment. Here, we explore whether we can identify those patients using circulating cell-free DNA (cfDNA) methylation. We first defined a panel of RCC-specific circulating methylation markers by intersecting differentially methylated regions from a publicly available dataset with known RCC methylation markers from the literature. The resulting RCC-specific methylation marker panel of 22 markers was subsequently evaluated for an association with rapid progression by methylated DNA sequencing (MeD-seq) in serum from 10 HBDs and 34 RCC patients with a good or intermediate prognosis starting WW in the IMPACT-RCC study. Patients with an elevated RCC-specific methylation score compared to HBDs had a shorter progression-free survival (PFS, p = 0.018), but not a shorter WW-time (p = 0.15). Cox proportional hazards regression showed that only the International Metastatic Renal Cell Carcinoma Database Consortium (IMDC) criteria were significantly associated with WW time (HR 2.01, p = 0.01), whereas only our RCC-specific methylation score (HR 4.45, p = 0.02) was significantly associated with PFS. The results of this study suggest that cfDNA methylation is predictive of PFS but not WW.

Circulating tumour DNA as biomarker for rectal cancer: A systematic review and meta-analyses

Background

Circulating tumour DNA (ctDNA) has been established as a promising (prognostic) biomarker with the potential to personalise treatment in cancer patients. The objective of this systematic review is to provide an overview of the current literature and the future perspectives of ctDNA in non-metastatic rectal cancer.

Methods

A comprehensive search for studies published prior to the 4^th of October 2022 was conducted in Embase, Medline, Cochrane, Google scholar, and Web of Science. Only peer-reviewed original articles and ongoing clinical trials investigating the association between ctDNA and oncological outcomes in non-metastatic rectal cancer patients were included. Meta-analyses were performed to pool hazard ratios (HR) for recurrence-free survival (RFS).

Results

A total of 291 unique records were screened, of which 261 were original publications and 30 ongoing trials. Nineteen original publications were reviewed and discussed, of which seven provided sufficient data for meta-analyses on the association between the presence of post-treatment ctDNA and RFS. Results of the meta-analyses demonstrated that ctDNA analysis can be used to stratify patients into very high and low risk groups for recurrence, especially when detected after neoadjuvant treatment (HR for RFS: 9.3 [4.6 - 18.8]) and after surgery (HR for RFS: 15.5 [8.2 - 29.3]). Studies investigated different types of assays and used various techniques for the detection and quantification of ctDNA.

Conclusions

This literature overview and meta-analyses provide evidence for the strong association between ctDNA and recurrent disease. Future research should focus on the feasibility of ctDNA-guided treatment and follow-up strategies in rectal cancer. A blueprint for agreed-upon timing, preprocessing, and assay techniques is needed to empower adaptation of ctDNA into daily practice.

Cost-Effectiveness of Abemaciclib in Early Breast Cancer Patients: One Size Fits All or Tailoring to Patients' Needs?

Purpose

The addition of two years of abemaciclib treatment to standard adjuvant endocrine therapy in all patients with high risk ER+, HER2- early breast cancer (EBC) has been approved by the US Food and Drug Administration (FDA). Pre-selection of patients with an immediate risk of recurrence within the group of clinically high risk patients using detection of minimal residual disease (MRD) using patient-informed circulating tumor DNA assays during follow-up could enhance efficacy. Here, we investigate the cost-effectiveness of the addition of two years abemaciclib in all high risk HR+, HER2- patients and in MRD-guided high risk patients only.

Methods

Two semi-Markov models were developed to evaluate the cost-effectiveness of adding two years of abemaciclib compared to "standard treatment": 1) "abemaciclib all" and 2) "MRD-guided abemaciclib" using MRD-guidance. Data of the MonarchE trial were used to model the invasive disease-free survival (iDFS). Since iDFS and overall survival (OS) data of abemaciclib were currently limited, abemaciclib effects were extrapolated using a favorable, intermediate and unfavorable effect scenario.

Results

The addition of abemaciclib in all high-risk EBC patients prolonged iDFS slightly (0.04 additional quality adjusted life years (QALYs)) and led to higher costs compared to standard ET, leading to a high incremental cost effectiveness ratio (ICER) of €1,551,876/QALY. Neither the favorable effect scenario (additional 1.09 QALYs) was cost-effective (ICER €62,935/QALY), using a willingness-to-pay threshold of €50,000/QALY. The "MRD-guided abemaciclib" strategy resulted in lower costs and an increase in QALYs (1.27) compared to "standard treatment" in the unfavorable effect scenario.

Conclusion

The addition of abemaciclib to adjuvant endocrine therapy in all high-risk ER+, HER2- EBC patients is not cost-effective. However, using MRD detection to justify the addition of abemaciclib treatment dominates standard treatment in this cost-effectiveness analysis. Further evaluation of MRD detection in EBC by means of prospective clinical trials assessing clinical utility is recommended and promising in terms of cost-effectiveness.

GRHL2-controlled gene expression networks in luminal breast cancer

Grainyhead like 2 (GRHL2) is an essential transcription factor for development and function of epithelial tissues. It has dual roles in cancer by supporting tumor growth while suppressing epithelial to mesenchymal transitions (EMT). GRHL2 cooperates with androgen and estrogen receptors (ER) to regulate gene expression. We explore genome wide GRHL2 binding sites conserved in three ER⍺/GRHL2 positive luminal breast cancer cell lines by ChIP-Seq. Interaction with the ER⍺/FOXA1/GATA3 complex is observed, however, only for a minor fraction of conserved GRHL2 peaks. We determine genome wide transcriptional dynamics in response to loss of GRHL2 by nascent RNA Bru-seq using an MCF7 conditional knockout model. Integration of ChIP- and Bru-seq pinpoints candidate direct GRHL2 target genes in luminal breast cancer. Multiple connections between GRHL2 and proliferation are uncovered, including transcriptional activation of ETS and E2F transcription factors. Among EMT-related genes, direct regulation of CLDN4 is corroborated but several targets identified in other cells (including CDH1 and ZEB1) are ruled out by both ChIP- and Bru-seq as being directly controlled by GRHL2 in luminal breast cancer cells. Gene clusters correlating positively (including known GRHL2 targets such as ErbB3, CLDN4/7) or negatively (including TGFB1 and TGFBR2) with GRHL2 in the MCF7 knockout model, display similar correlation with GRHL2 in ER positive as well as ER negative breast cancer patients. Altogether, this study uncovers gene sets regulated directly or indirectly by GRHL2 in luminal breast cancer, identifies novel GRHL2-regulated genes, and points to distinct GRHL2 regulation of EMT in luminal breast cancer cells. Video Abstract.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

CABA-V7: a prospective biomarker selected trial of cabazitaxel treatment in AR-V7 positive prostate cancer patients

BACKGROUND

Metastatic castration-resistant prostate cancer (mCRPC) patients with positive AR-V7 expression in their circulating tumour cells (CTCs) rarely derive benefit from abiraterone and enzalutamide.

DESIGN

We performed a prospective, multicenter, single arm phase II clinical trial (CABA-V7) in mCRPC patients previously treated with docetaxel and androgen deprivation therapy.

OBJECTIVE

In this trial, we investigated whether cabazitaxel treatment resulted in clinically meaningful PSA response rates in patients with positive CTC-based AR-V7 expression and collected liquid biopsies for genomic profiling.

RESULTS

Cabazitaxel was found to be modestly effective, with only 12% of these patients obtaining a PSA response. Genomic profiling revealed that CTC-based AR-V7 expression was not associated with other known mCRPC-associated alterations. CTC-based AR-V7 status and dichotomised CTC counts were observed as independent prognostic markers at baseline.

CONCLUSIONS

AR-V7 positivity predicted poor overall survival (OS). However, cabazitaxel-treated AR-V7 positive patients and those lacking AR-V7 positivity, who received cabazitaxel as standard of care, appeared to have similar OS. Therefore, despite the low response rate, cabazitaxel may still be an effective treatment in this poor prognosis, AR-V7 positive patient population.

Copy number variants as modifiers of breast cancer risk for BRCA1/BRCA2 pathogenic variant carriers

The contribution of germline copy number variants (CNVs) to risk of developing cancer in individuals with pathogenic BRCA1 or BRCA2 variants remains relatively unknown. We conducted the largest genome-wide analysis of CNVs in 15,342 BRCA1 and 10,740 BRCA2 pathogenic variant carriers. We used these results to prioritise a candidate breast cancer risk-modifier gene for laboratory analysis and biological validation. Notably, the HR for deletions in BRCA1 suggested an elevated breast cancer risk estimate (hazard ratio (HR) = 1.21), 95% confidence interval (95% CI = 1.09-1.35) compared with non-CNV pathogenic variants. In contrast, deletions overlapping SULT1A1 suggested a decreased breast cancer risk (HR = 0.73, 95% CI 0.59-0.91) in BRCA1 pathogenic variant carriers. Functional analyses of SULT1A1 showed that reduced mRNA expression in pathogenic BRCA1 variant cells was associated with reduced cellular proliferation and reduced DNA damage after treatment with DNA damaging agents. These data provide evidence that deleterious variants in BRCA1 plus SULT1A1 deletions contribute to variable breast cancer risk in BRCA1 carriers.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

Functional RECAP (REpair CAPacity) assay identifies homologous recombination deficiency undetected by DNA-based BRCAness tests

Germline BRCA1/2 mutation status is predictive for response to Poly-[ADP-Ribose]-Polymerase (PARP) inhibitors in breast cancer (BC) patients. However, non-germline BRCA1/2 mutated and homologous recombination repair deficient (HRD) tumors are likely also PARP-inhibitor sensitive. Clinical validity and utility of various HRD biomarkers are under investigation. The REpair CAPacity (RECAP) test is a functional method to select HRD tumors based on their inability to form RAD51 foci. We investigated whether this functional test defines a similar group of HRD tumors as DNA-based tests. An HRD enriched cohort (n = 71; 52 primary and 19 metastatic BCs) selected based on the RECAP test (26 RECAP-HRD; 37%), was subjected to DNA-based HRD tests (i.e., Classifier of HOmologous Recombination Deficiency (CHORD) and BRCA1/2-like classifier). Whole genome sequencing (WGS) was carried out for 38 primary and 19 metastatic BCs. The RECAP test identified all bi-allelic BRCA deficient samples (n = 15) in this cohort. RECAP status partially correlated with DNA-based HRD test outcomes (70% concordance for both RECAP-CHORD and RECAP-BRCA1/2-like classifier). RECAP selected additional samples unable to form RAD51 foci, suggesting that this functional assay identified deficiencies in other DNA repair genes, which could also result in PARP-inhibitor sensitivity. Direct comparison of these HRD tests in clinical trials will be required to evaluate the optimal predictive test for clinical decision making.

Liquid Biopsy Based Circulating Biomarkers in Metastatic Prostate Cancer

Prostate cancer is the most dominant male malignancy worldwide. The clinical presentation of prostate cancer ranges from localized indolent to rapidly progressing lethal metastatic disease. Despite a decline in death rate over the past years, with the advent of early diagnosis and new treatment options, challenges remain towards the management of metastatic prostate cancer, particularly metastatic castration sensitive prostate cancer (mCSPC) and castration resistant prostate cancer (mCRPC). Current treatments involve a combination of chemotherapy with androgen deprivation therapy and/or androgen receptor signalling inhibitors. However, treatment outcomes are heterogeneous due to significant tumor heterogeneity indicating a need for better prognostic biomarkers to identify patients with poor outcomes. Liquid biopsy has opened a plethora of opportunities from early diagnosis to (personalized) therapeutic disease interventions. In this review, we first provide recent insights about (metastatic) prostate cancer and its current treatment landscape. We highlight recent studies involving various circulating biomarkers such as circulating tumor cells, genetic markers, circulating nucleic acids, extracellular vesicles, tumor-educated platelets, and the secretome from (circulating) tumor cells and tumor microenvironment in metastatic prostate cancer. The comprehensive array of biomarkers can provide a powerful approach to understanding the spectrum of prostate cancer disease and guide in developing improved and personalized treatments for patients.

ESR1 Methylation Measured in Cell-Free DNA to Evaluate Endocrine Resistance in Metastatic Breast Cancer Patients

ESR1 methylation was proposed as mechanism for endocrine resistance in metastatic breast cancer patients. To evaluate its potential as a minimally invasive biomarker, we investigated the feasibility of measuring ESR1 methylation in cell-free DNA (cfDNA) and its association with endocrine resistance. First, we provided evidence that demethylation in vitro restores ER expression. Subsequently, we found that ESR1 methylation in cfDNA was not enriched in endocrine-resistant versus endocrine-sensitive patients. Interestingly, we found a correlation between ESR1 methylation and age. Publicly available data confirm an age-related increase in ESR1 methylation in leukocytes, confounding the determination of the ESR1 methylation status of tumors using cfDNA.

Genome-wide aneuploidy detected by mFast-SeqS in circulating cell-free DNA is associated with poor response to pembrolizumab in patients with advanced urothelial cancer

Second-line treatment with immune checkpoint inhibition in patients with metastatic urothelial cancer (mUC) has a low success rate (~ 20%). Circulating tumour-derived DNA (ctDNA) levels may guide patient stratification, provided that an affordable and robust assay is available. Here, we investigate whether the modified fast aneuploidy screening test-sequencing system (mFast-SeqS) may provide such an assay. To this end, mFast-SeqS was performed on cell-free DNA (cfDNA) from 74 patients with mUC prior to treatment with pembrolizumab. Results were associated with corresponding tissue-based profiles, plasma-based variant allele frequencies (VAFs) and clinical response. We found that plasma-derived mFast-SeqS-based aneuploidy scores significantly correlated with those observed in the corresponding tumour tissue as well as with the ctDNA level in the plasma. In multivariate logistic regression analysis, a high aneuploidy score was independently associated with lack of clinical benefit from treatment with pembrolizumab. In conclusion, mFast-SeqS provides a patient-friendly, high-throughput and affordable method to estimate ctDNA level. Following independent validation, this test could be used to stratify mUC patients for response prior to the initiation of treatment with pembrolizumab.

The prognostic and predictive value of ESR1 fusion gene transcripts in primary breast cancer

Background

In breast cancer (BC), recurrent fusion genes of estrogen receptor alpha (ESR1) and AKAP12, ARMT1 and CCDC170 have been reported. In these gene fusions the ligand binding domain of ESR1 has been replaced by the transactivation domain of the fusion partner constitutively activating the receptor. As a result, these gene fusions can drive tumor growth hormone independently as been shown in preclinical models, but the clinical value of these fusions have not been reported. Here, we studied the prognostic and predictive value of different frequently reported ESR1 fusion transcripts in primary BC.

Methods

We evaluated 732 patients with primary BC (131 ESR1-negative and 601 ESR1-positive cases), including two ER-positive BC patient cohorts: one cohort of 322 patients with advanced disease who received first-line endocrine therapy (ET) (predictive cohort), and a second cohort of 279 patients with lymph node negative disease (LNN) who received no adjuvant systemic treatment (prognostic cohort). Fusion gene transcript levels were measured by reverse transcriptase quantitative PCR. The presence of the different fusion transcripts was associated, in uni- and multivariable Cox regression analysis taking along current clinico-pathological characteristics, to progression free survival (PFS) during first-line endocrine therapy in the predictive cohort, and disease- free survival (DFS) and overall survival (OS) in the prognostic cohort.

Results

The ESR1-CCDC170 fusion transcript was present in 27.6% of the ESR1-positive BC subjects and in 2.3% of the ESR1-negative cases. In the predictive cohort, none of the fusion transcripts were associated with response to first-line ET. In the prognostic cohort, the median DFS and OS were respectively 37 and 93 months for patients with an ESR1-CCDC170 exon 8 gene fusion transcript and respectively 91 and 212 months for patients without this fusion transcript. In a multivariable analysis, this ESR1-CCDC170 fusion transcript was an independent prognostic factor for DFS (HR) (95% confidence interval (CI): 1.8 (1.2–2.8), P = 0.005) and OS (HR (95% CI: 1.7 (1.1–2.7), P = 0.023).

Conclusions

Our study shows that in primary BC only ESR1-CCDC170 exon 8 gene fusion transcript carries prognostic value. None of the ESR1 fusion transcripts, which are considered to have constitutive ER activity, was predictive for outcome in BC with advanced disease treated with endocrine treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09265-1.

Apolipoprotein B mRNA-Editing Catalytic Polypeptide-Like-Induced Protein Changes in Estrogen Receptor-Positive, Human Epidermal Growth Factor Receptor 2-Negative Breast Cancer Throughout Disease Progression

PURPOSE

Apolipoprotein B mRNA-Editing Catalytic Polypeptide-like (APOBEC) enzymes are mutagenic factors contributing to tumor progression and therapy resistance. However, the effects of APOBEC-induced protein changes have not been systematically assessed. Here, we describe the effects of APOBEC on the coding sequence in primary and metastatic estrogen receptor–positive (ER+)/human epidermal growth factor receptor 2–negative (HER2–) breast cancer (BC).

METHODS

We determined the enrichment of amino acid (AA) changes resulting from APOBEC mutagenesis in 323 primary BC tumors and 424 metastatic breast cancer (mBC) lesions via comparison with a simulated mutational genomic landscape not under selection pressure. We subsequently explored genes with recurrent APOBEC-associated AA changes and investigated the clonality of individual APOBEC-associated mutations. Using public sequencing data from an independent primary BC and mBC cohort, we further confirm our findings by reporting genes having these enriched AA changes in an APOBEC context.

RESULTS

Our analysis demonstrated that several APOBEC-derived AA changes are significantly enriched compared with a simulated AA change distribution drawn at random. Among the enriched AA changes, Glutamate(E)>Lysine(K) and Glutamate(E)>Glutamine(Q) were mostly found at hotspots in oncogenes, whereas termination codons (Glutamine[Q]>STOP[X] and Serine[S]>STOP[X]) occurred in tumor suppressor genes and, mostly, not at hotspot locations. These mutations are found in genes contributing to BC initiation, eg, introduction of termination codons in TP53, MAP3K1, and CDH1 (in lobular BC) and two oncogenic hotspots in PIK3CA (p.E542K and p.E545K). In endocrine-resistant BC, we observed APOBEC-induced termination codons at ER-modulating genes, KMT2C, ARID1A, NF1, and ZFHX3. Finally, in mBC, compared with single PIK3CA mutations, dual PIK3CA mutations occurred more frequently in an APOBEC context (Fisher's exact P < .001).

CONCLUSION

Our results show that APOBEC mutagenesis recurrently targets various known drivers of BC initiation, progression, and endocrine resistance.

APOBEC drives early, metastatic and endocrine resistant disease in ER+/HER2- breast cancer

A pipeline for copy number profiling of single circulating tumor cells to assess intra-patient tumor heterogeneity

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike‐in experiments using MCF‐7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI‐g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient‐derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient‐derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

Fusion transcripts and their genomic breakpoints in polyadenylated and ribosomal RNA-minus RNA sequencing data

BACKGROUND

Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non-poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints.

RESULTS

We have developed an algorithm, Dr. Disco, that searches for fusion transcripts by taking an entire reference genome into account as search space. This includes exons but also introns, intergenic regions, and sequences that do not meet splice junction motifs. Using 1,275 RNA-seq samples, we investigated to what extent genomic breakpoints can be extracted from RNA-seq data and their implications regarding poly(A)-enriched and ribosomal RNA-minus RNA-seq data. Comparison with whole-genome sequencing data revealed that most genomic breakpoints are not, or minimally, transcribed while, in contrast, the genomic breakpoints of all 32 TMPRSS2-ERG-positive tumours were present at RNA level. We also revealed tumours in which the ERG breakpoint was located before ERG, which co-existed with additional deletions and messenger RNA that incorporated intergenic cryptic exons. In breast cancer we identified rearrangement hot spots near CCND1 and in glioma near CDK4 and MDM2 and could directly associate this with increased expression. Furthermore, in all datasets we find fusions to intergenic regions, often spanning multiple cryptic exons that potentially encode neo-antigens. Thus, fusion transcripts other than classical gene-to-gene fusions are prominently present and can be identified using RNA-seq.

CONCLUSION

By using the full potential of non-poly(A)-enriched RNA-seq data, sophisticated analysis can reliably identify expressed genomic breakpoints and their transcriptional effects.

SPEN is required for Xist upregulation during initiation of X chromosome inactivation

At initiation of X chromosome inactivation (XCI), Xist is monoallelically upregulated from the future inactive X (Xi) chromosome, overcoming repression by its antisense transcript Tsix. Xist recruits various chromatin remodelers, amongst them SPEN, which are involved in silencing of X-linked genes in cis and establishment of the Xi. Here, we show that SPEN plays an important role in initiation of XCI. Spen null female mouse embryonic stem cells (ESCs) are defective in Xist upregulation upon differentiation. We find that Xist-mediated SPEN recruitment to the Xi chromosome happens very early in XCI, and that SPEN-mediated silencing of the Tsix promoter is required for Xist upregulation. Accordingly, failed Xist upregulation in Spen-/- ESCs can be rescued by concomitant removal of Tsix. These findings indicate that SPEN is not only required for the establishment of the Xi, but is also crucial in initiation of the XCI process.

High-throughput and affordable genome-wide methylation profiling of circulating cell-free DNA by methylated DNA sequencing (MeD-seq) of LpnPI digested fragments

Background

DNA methylation detection in liquid biopsies provides a highly promising and much needed means for real-time monitoring of disease load in advanced cancer patient care. Compared to the often-used somatic mutations, tissue- and cancer-type specific epigenetic marks affect a larger part of the cancer genome and generally have a high penetrance throughout the tumour. Here, we describe the successful application of the recently described MeD-seq assay for genome-wide DNA methylation profiling on cell-free DNA (cfDNA). The compatibility of the MeD-seq assay with different types of blood collection tubes, cfDNA input amounts, cfDNA isolation methods, and vacuum concentration of samples was evaluated using plasma from both metastatic cancer patients and healthy blood donors (HBDs). To investigate the potential value of cfDNA methylation profiling for tumour load monitoring, we profiled paired samples from 8 patients with resectable colorectal liver metastases (CRLM) before and after surgery.

Results

The MeD-seq assay worked on plasma-derived cfDNA from both EDTA and CellSave blood collection tubes when at least 10 ng of cfDNA was used. From the 3 evaluated cfDNA isolation methods, both the manual QIAamp Circulating Nucleic Acid Kit (Qiagen) and the semi-automated Maxwell® RSC ccfDNA Plasma Kit (Promega) were compatible with MeD-seq analysis, whereas the QiaSymphony DSP Circulating DNA Kit (Qiagen) yielded significantly fewer reads when compared to the QIAamp kit (p < 0.001). Vacuum concentration of samples before MeD-seq analysis was possible with samples in AVE buffer (QIAamp) or water, but yielded inconsistent results for samples in EDTA-containing Maxwell buffer. Principal component analysis showed that pre-surgical samples from CRLM patients were very distinct from HBDs, whereas post-surgical samples were more similar. Several described methylation markers for colorectal cancer monitoring in liquid biopsies showed differential methylation between pre-surgical CRLM samples and HBDs in our data, supporting the validity of our approach. Results for MSC, ITGA4, GRIA4, and EYA4 were validated by quantitative methylation specific PCR.

Conclusions

The MeD-seq assay provides a promising new method for cfDNA methylation profiling. Potential future applications of the assay include marker discovery specifically for liquid biopsy analysis as well as direct use as a disease load monitoring tool in advanced cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01177-4.

Spatial immunophenotypes predict response to anti-PD1 treatment and capture distinct paths of T cell evasion in triple negative breast cancer

Only a subgroup of triple-negative breast cancer (TNBC) responds to immune checkpoint inhibitors (ICI). To better understand lack of response to ICI, we analyze 681 TNBCs for spatial immune cell contextures in relation to clinical outcomes and pathways of T cell evasion. Excluded, ignored and inflamed phenotypes can be captured by a gene classifier that predicts prognosis of various cancers as well as anti-PD1 response of metastatic TNBC patients in a phase II trial. The excluded phenotype, which is associated with resistance to anti-PD1, demonstrates deposits of collagen-10, enhanced glycolysis, and activation of TGFβ/VEGF pathways; the ignored phenotype, also associated with resistance to anti-PD1, shows either high density of CD163+ myeloid cells or activation of WNT/PPARγ pathways; whereas the inflamed phenotype, which is associated with response to anti-PD1, revealed necrosis, high density of CLEC9A+ dendritic cells, high TCR clonality independent of neo-antigens, and enhanced expression of T cell co-inhibitory receptors.

Liquid biopsy in esophageal cancer: a case report of false-positive circulating tumor DNA detection due to clonal hematopoiesis

Circulating tumor DNA (ctDNA) analysis is a promising non-invasive technique for active surveillance after chemoradiotherapy for locally advanced resectable esophageal carcinoma. In other malignancies false-positive results in ctDNA analysis have been reported due to clonal hematopoiesis. In this case, we present a 66-year-old male who had adenocarcinoma of the gastroesophageal junction for which he received neoadjuvant chemoradiotherapy and underwent a transhiatal esophagectomy. Postoperatively our patient received follow-up with ctDNA analysis using next generation sequencing (NGS) and droplet digital PCR (ddPCR). This case report illustrates a number of the current challenges in ctDNA diagnostics in esophageal carcinoma. Firstly, the TP53 c.524G>A; p.R175H mutation that was found in preoperative tumor biopsies became detectable in ctDNA only after distant metastases had already been confirmed by clinical symptoms and standard imaging- and biopsy techniques. Secondly our patient repeatedly had false-positive outcomes of ctDNA analysis. Genomic analysis of white blood cells revealed that the origin of these discordant mutations lies in clonal hematopoiesis. Failure to detect TP53 c.524G>A; p.R175H in cell-free DNA (cfDNA) is most likely due to the amount of ctDNA in the cfDNA fraction being below the limit of detection for NGS and ddPCR analyses. Clinicians should be aware of the possibility of finding mutations originating from clonal hematopoiesis when using ctDNA analysis during active surveillance for esophageal carcinoma. We recommend correlation of mutations in cfDNA with mutations in tumor biopsies.

A panel of DNA methylation markers for the classification of consensus molecular subtypes 2 and 3 in patients with colorectal cancer

Consensus molecular subtypes (CMSs) can guide precision treatment of colorectal cancer (CRC). We aim to identify methylation markers to distinguish between CMS2 and CMS3 in patients with CRC, for which an easy test is currently lacking. To this aim, fresh‐frozen tumor tissue of 239 patients with stage I‐III CRC was analyzed. Methylation profiles were obtained using the Infinium HumanMethylation450 BeadChip. We performed adaptive group‐regularized logistic ridge regression with post hoc group‐weighted elastic net marker selection to build prediction models for classification of CMS2 and CMS3. The Cancer Genome Atlas (TCGA) data were used for validation. Group regularization of the probes was done based on their location either relative to a CpG island or relative to a gene present in the CMS classifier, resulting in two different prediction models and subsequently different marker panels. For both panels, even when using only five markers, accuracies were > 90% in our cohort and in the TCGA validation set. Our methylation marker panel accurately distinguishes between CMS2 and CMS3. This enables development of a targeted assay to provide a robust and clinically relevant classification tool for CRC patients.

Association of germline genetic variants with breast cancer-specific survival in patient subgroups defined by clinic-pathological variables related to tumor biology and type of systemic treatment

BACKGROUND

Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients.

METHODS

We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15).

RESULTS

Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy.

CONCLUSIONS

We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

Transcriptomic Properties of HER2+ Ductal Carcinoma In Situ of the Breast Associate with Absence of Immune Cells

The identification of transcriptomic alterations of HER2+ ductal carcinoma in situ (DCIS) that are associated with the density of tumor-infiltrating lymphocytes (TILs) could contribute to optimizing choices regarding the potential benefit of immune therapy. We compared the gene expression profile of TIL-poor HER2+ DCIS to that of TIL-rich HER2+ DCIS. Tumor cells from 11 TIL-rich and 12 TIL-poor DCIS cases were micro-dissected for RNA isolation. The Ion AmpliSeq Transcriptome Human Gene Expression Kit was used for RNA sequencing. After normalization, a Mann-Whitney rank sum test was used to analyze differentially expressed genes between TIL-poor and TIL-rich HER2+ DCIS. Whole tissue sections were immunostained for validation of protein expression. We identified a 29-gene expression profile that differentiated TIL-rich from TIL-poor HER2+ DCIS. These genes included CCND3, DUSP10 and RAP1GAP, which were previously described in breast cancer and cancer immunity and were more highly expressed in TIL-rich DCIS. Using immunohistochemistry, we found lower protein expression in TIL-rich DCIS. This suggests regulation of protein expression at the posttranslational level. We identified a gene expression profile of HER2+ DCIS cells that was associated with the density of TILs. This classifier may guide towards more rationalized choices regarding immune-mediated therapy in HER2+ DCIS, such as targeted vaccine therapy.

Prospective Evaluation of a Circulating Tumor Cell Sensitivity Profile to Predict Response to Cisplatin Chemotherapy in Metastatic Breast Cancer Patients

Background

Cisplatin (cDDP) has regained interest for metastatic breast cancer (MBC) patients, given the platinum sensitivity in subtypes and better manageable toxicity. Here, the primary aim was to determine whether molecular characteristics of circulating tumor cells (CTCs) could identify patients responding to cDDP and to describe the outcomes to cDDP monotherapy in a large group of MBC patients pretreated with anthracycline- and taxane-based treatments.

Methods

Based on cell line data, a CTC-cDDP-sensitivity profile was generated. Applying an A’Herns single-stage phase II design, further investigation was considered worthwhile if 5/10 patients with a favorable profile responded to cDDP. Patients received 70mg/m² cDDP every three weeks, CTCs were enumerated and the CTC-cDDP-sensitivity profile was determined. In total, 65 heavily pretreated MBC patients (77% received ≥2 lines of previous chemotherapy for MBC) were eligible for the per-protocol analysis. Primary endpoint was response rate, secondary endpoints included best observed response, progression-free survival (PFS) and overall survival (OS).

Results

The best observed response during cDDP therapy was a partial response in 7% and stable disease in 56% of the patients. None of the patients with a favorable CTC-cDDP-sensitivity profile had a response. The median baseline CTC count was 8 (range 0-3254). Patients with <5 CTCs had a better PFS and OS than patients with ≥5 CTCs (median PFS 4.5 months (95%CI 2.38-6.62) vs. 2.1 months [(95%CI 1.34-2.80)(p=0.009)] and median OS 13.1 months (95%CI 9.89-16.33) vs. 5.6 months [(95%CI 3.60-7.64)(p=0.003)]. No other factors than CTC count were associated with outcome to cDDP therapy, including triple-negative breast cancer versus ER-positive tumors.

Conclusions

The CTC-cDDP-sensitivity profile was unable to select patients responding to cDDP monotherapy. In an unselected group of heavily pretreated MBC patients, cDDP yields outcomes comparable to other chemotherapeutic regimens for heavily pretreated MBC patients. CTC count was the only factor associated with outcome in these patients.

Clinical Trial Registration

(https://www.trialregister.nl/trial/3885, identifier NTR4046)

Estrogens and Progestogens in Triple Negative Breast Cancer: Do They Harm?

Triple-negative breast cancers (TNBC) occur more frequently in younger women and do not express estrogen receptor (ER) nor progesterone receptor (PR), and are therefore often considered hormone-insensitive. Treatment of premenopausal TNBC patients almost always includes chemotherapy, which may lead to premature ovarian insufficiency (POI) and can severely impact quality of life. Hormone replacement therapy (HRT) is contraindicated for patients with a history of hormone-sensitive breast cancer, but the data on safety for TNBC patients is inconclusive, with a few randomized trials showing increased risk-ratios with wide confidence intervals for recurrence after HRT. Here, we review the literature on alternative pathways from the classical ER/PR. We find that for both estrogens and progestogens, potential alternatives exist for exerting their effects on TNBC, ranging from receptor conversion, to alternative receptors capable of binding estrogens, as well as paracrine pathways, such as RANK/RANKL, which can cause progestogens to indirectly stimulate growth and metastasis of TNBC. Finally, HRT may also influence other hormones, such as androgens, and their effects on TNBCs expressing androgen receptors (AR). Concluding, the assumption that TNBC is completely hormone-insensitive is incorrect. However, the direction of the effects of the alternative pathways is not always clear, and will need to be investigated further.

Optimizing Nanopore sequencing-based detection of structural variants enables individualized circulating tumor DNA-based disease monitoring in cancer patients

Here, we describe a novel approach for rapid discovery of a set of tumor-specific genomic structural variants (SVs), based on a combination of low coverage cancer genome sequencing using Oxford Nanopore with an SV calling and filtering pipeline. We applied the method to tumor samples of high-grade ovarian and prostate cancer patients and validated on average ten somatic SVs per patient with breakpoint-spanning PCR mini-amplicons. These SVs could be quantified in ctDNA samples of patients with metastatic prostate cancer using a digital PCR assay. The results suggest that SV dynamics correlate with and may improve existing treatment-response biomarkers such as PSA. https://github.com/UMCUGenetics/SHARC .