Whole-genome sequencing of triple-negative breast cancers in a population-based clinical study

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Luminal androgen receptor subtype and tumor-infiltrating lymphocytes groups based on triple-negative breast cancer molecular subclassification

In our previous study, we developed a triple-negative breast cancer (TNBC) subtype classification that correlated with the TNBC molecular subclassification. In this study, we aimed to evaluate the predictor variables of this subtype classification on the whole slide and to validate the model's performance by using an external test set. We explored the characteristics of this subtype classification and investigated genomic alterations, including genomic scar signature scores. First, TNBC was classified into the luminal androgen receptor (LAR) and non-luminal androgen receptor (non-LAR) subtypes based on the AR Allred score (≥ 6 and < 6, respectively). Then, the non-LAR subtype was further classified into the lymphocyte-predominant (LP), lymphocyte-intermediate (LI), and lymphocyte-depleted (LD) groups based on stromal tumor-infiltrating lymphocytes (TILs) (< 20%, > 20% but < 60%, and ≥ 60%, respectively). This classification showed fair agreement with the molecular classification in the test set. The LAR subtype was characterized by a high rate of PIK3CA mutation, CD274 (encodes PD-L1) and PDCD1LG2 (encodes PD-L2) deletion, and a low homologous recombination deficiency (HRD) score. The non-LAR LD TIL group was characterized by a high frequency of NOTCH2 and MYC amplification and a high HRD score.

Kataegis in clinical and molecular subgroups of primary breast cancer

Kataegis is a hypermutation phenomenon characterized by localized clusters of single base pair substitution (SBS) reported in multiple cancer types. Despite a high frequency in breast cancer, large-scale analyses of kataegis patterns and associations with clinicopathological and molecular variables in established breast cancer subgroups are lacking. Therefore, WGS profiled primary breast cancers (n = 791) with associated clinical and molecular data layers, like RNA-sequencing data, were analyzed for kataegis frequency, recurrence, and associations with genomic contexts and functional elements, transcriptional patterns, driver alterations, homologous recombination deficiency (HRD), and prognosis in tumor subgroups defined by ER, PR, and HER2/ERBB2 status. Kataegis frequency was highest in the HER2-positive(p) subgroups, including both ER-negative(n)/positive(p) tumors (ERnHER2p/ERpHER2p). In TNBC, kataegis was neither associated with PAM50 nor TNBC mRNA subtypes nor with distant relapse in chemotherapy-treated patients. In ERpHER2n tumors, kataegis was associated with aggressive characteristics, including PR-negativity, molecular Luminal B subtype, higher mutational burden, higher grade, and expression of proliferation-associated genes. Recurrent kataegis loci frequently targeted regions commonly amplified in ER-positive tumors, while few recurrent loci were observed in TNBC. SBSs in kataegis loci appeared enriched in regions of open chromatin. Kataegis status was not associated with HRD in any subgroup or with distinct transcriptional patterns in unsupervised or supervised analysis. In summary, kataegis is a common hypermutation phenomenon in established breast cancer subgroups, particularly in HER2p subgroups, coinciding with an aggressive tumor phenotype in ERpHER2n disease. In TNBC, the molecular implications and associations of kataegis are less clear, including its prognostic value.

PALB2-mutated human mammary cells display a broad spectrum of morphological and functional abnormalities induced by increased TGFβ signaling

Heterozygous mutations in any of three major genes, BRCA1, BRCA2 and PALB2, are associated with high-risk hereditary breast cancer susceptibility frequently seen as familial disease clustering. PALB2 is a key interaction partner and regulator of several vital cellular activities of BRCA1 and BRCA2, and is thus required for DNA damage repair and alleviation of replicative and oxidative stress. Little is however known about how PALB2-deficiency affects cell function beyond that, especially in the three-dimensional setting, and also about its role during early steps of malignancy development. To answer these questions, we have generated biologically relevant MCF10A mammary epithelial cell lines with mutations that are comparable to certain clinically important PALB2 defects. We show in a non-cancerous background how both mono- and biallelically PALB2-mutated cells exhibit gross spontaneous DNA damage and mitotic aberrations. Furthermore, PALB2-deficiency disturbs three-dimensional spheroid morphology, increases the migrational capacity and invasiveness of the cells, and broadly alters their transcriptome profiles. TGFβ signaling and KRT14 expression are enhanced in PALB2-mutated cells and their inhibition and knock down, respectively, lead to partial restoration of cell functions. KRT14-positive cells are also more abundant with DNA damage than KRT14-negative cells. The obtained results indicate comprehensive cellular changes upon PALB2 mutations, even in the presence of half dosage of wild type PALB2 and demonstrate how PALB2 mutations may predispose their carriers to malignancy.

Longitudinal profiling identifies co-occurring BRCA1/2 reversions, TP53BP1, RIF1 and PAXIP1 mutations in PARP inhibitor-resistant advanced breast cancer

BACKGROUND

Resistance to therapies that target homologous recombination deficiency (HRD) in breast cancer limits their overall effectiveness. Multiple, preclinically validated, mechanisms of resistance have been proposed, but their existence and relative frequency in clinical disease are unclear, as is how to target resistance.

PATIENTS AND METHODS

Longitudinal mutation and methylation profiling of circulating tumour (ct)DNA was carried out in 47 patients with metastatic BRCA1-, BRCA2- or PALB2-mutant breast cancer treated with HRD-targeted therapy who developed progressive disease-18 patients had primary resistance and 29 exhibited response followed by resistance. ctDNA isolated at multiple time points in the patient treatment course (before, on-treatment and at progression) was sequenced using a novel >750-gene intron/exon targeted sequencing panel. Where available, matched tumour biopsies were whole exome and RNA sequenced and also used to assess nuclear RAD51.

RESULTS

BRCA1/2 reversion mutations were present in 60% of patients and were the most prevalent form of resistance. In 10 cases, reversions were detected in ctDNA before clinical progression. Two new reversion-based mechanisms were identified: (i) intragenic BRCA1/2 deletions with intronic breakpoints; and (ii) intragenic BRCA1/2 secondary mutations that formed novel splice acceptor sites, the latter being confirmed by in vitro minigene reporter assays. When seen before commencing subsequent treatment, reversions were associated with significantly shorter time to progression. Tumours with reversions retained HRD mutational signatures but had functional homologous recombination based on RAD51 status. Although less frequent than reversions, nonreversion mechanisms [loss-of-function (LoF) mutations in TP53BP1, RIF1 or PAXIP1] were evident in patients with acquired resistance and occasionally coexisted with reversions, challenging the notion that singular resistance mechanisms emerge in each patient.

CONCLUSIONS

These observations map the prevalence of candidate drivers of resistance across time in a clinical setting, information with implications for clinical management and trial design in HRD breast cancers.

Free drug and ROS-responsive nanoparticle delivery of synergistic doxorubicin and olaparib combinations to triple negative breast cancer models

Combinations of the topoisomerase II inhibitor doxorubicin and the poly (ADP-ribose) polymerase inhibitor olaparib offer potential drug-drug synergy for the treatment of triple negative breast cancers (TNBC). In this study we performed in vitro screening of combinations of these drugs, administered directly or encapsulated within polymer nanoparticles, in both 2D and in 3D spheroid models of breast cancer. A variety of assays were used to evaluate drug potency, and calculations of combination index (CI) values indicated that synergistic effects of drug combinations occurred in a molar-ratio dependent manner. It is suggested that the mechanisms of synergy were related to enhancement of DNA damage as shown by the level of double-strand DNA breaks, and mechanisms of antagonism associated with mitochondrial mediated cell survival, as indicated by reactive oxygen species (ROS) generation. Enhanced drug delivery and potency was observed with nanoparticle formulations, with a greater extent of doxorubicin localised to cell nuclei as evidenced by microscopy, and higher cytotoxicity at the same time points compared to free drugs. Together, the work presented identifies specific combinations of doxorubicin and olaparib which were most effective in a panel of TNBC cell lines, explores the mechanisms by which these combined agents might act, and shows that formulation of these drug combinations into polymeric nanoparticles at specific ratios conserves synergistic action and enhanced potency in vitro compared to the free drugs.

Intragenic Rearrangement Burden Associates with Immune Cell Infiltration and Response to Immune Checkpoint Blockade in Cancer

Immune checkpoint blockade (ICB) can induce durable cancer remission. However, only a small subset of patients gains benefits. While tumor mutation burden (TMB) differentiates responders from nonresponders in some cases, it is a weak predictor in tumor types with low mutation rates. Thus, there is an unmet need to discover a new class of genetic aberrations that predict ICB responses in these tumor types. Here, we report analyses of pan-cancer whole genomes which revealed that intragenic rearrangement (IGR) burden is significantly associated with immune infiltration in breast, ovarian, esophageal, and endometrial cancers, particularly with increased M1 macrophage and CD8+ T-cell signatures. Multivariate regression against spatially counted tumor-infiltrating lymphocytes in breast, endometrial, and ovarian cancers suggested that IGR burden is a more influential covariate than other genetic aberrations in these cancers. In the MEDI4736 trial evaluating durvalumab in esophageal adenocarcinoma, IGR burden correlated with patient benefits. In the IMVigor210 trial evaluating atezolizumab in urothelial carcinoma, IGR burden increased with platinum exposure and predicted patient benefit among TMB-low, platinum-exposed tumors. Altogether, we have demonstrated that IGR burden correlates with T-cell inflammation and predicts ICB benefit in TMB-low, IGR-dominant tumors, and in platinum-exposed tumors.

Blocking LBH expression causes replication stress and sensitizes triple-negative breast cancer cells to ATR inhibitor treatment

Triple-negative (ER-PR-HER2-) breast cancers (TNBC) are highly aggressive and difficult to treat. TNBC exhibit high genomic instability, which enables them to adapt and become resistant to chemo/radiation therapy, leading to rapid disease relapse and mortality. The pro-survival factors that safeguard genome integrity in TNBC cells are poorly understood. LBH is an essential mammary stem cell-specific transcription regulator in the WNT pathway that is aberrantly overexpressed in TNBC, correlating with poor prognosis. Herein, we demonstrate a novel role for LBH in promoting TNBC cell survival. Depletion of LBH in multiple TNBC cell models triggered apoptotic cell death both in vitro and in vivo and led to S-G2M cell cycle delays. Mechanistically, LBH loss causes replication stress due to DNA replication fork stalling, leading to ssDNA breaks, ɣH2AX and 53BP1 nuclear foci formation, and activation of the ATR/CHK1 DNA damage response. Notably, ATR inhibition in combination with LBH downmodulation had a synergistic effect, boosting TNBC cell killing and blocking in vivo tumor growth. Our findings demonstrate, for the first time, that LBH protects the genome integrity of cancer cells by preventing replicative stress. Importantly, they uncover new synthetic lethal vulnerabilities in TNBC that could be exploited for future multi-modal precision medicine.

Accurate and sensitive mutational signature analysis with MuSiCal

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

Somatic mutations in four novel genes contribute to homologous recombination deficiency in breast cancer: a real-world clinical tumor sequencing study

Breast cancers involving mutations in homologous recombination (HR) genes, most commonly BRCA1 and BRCA2 (BRCA1/2), respond well to PARP inhibitors and platinum-based chemotherapy. However, except for these specific HR genes, it is not clear which other mutations contribute to homologous recombination defects (HRD). Here, we performed next-generation sequencing of tumor tissues and matched blood samples from 119 breast cancer patients using the OncoScreen Plus panel. Genomic mutation characteristics and HRD scores were analyzed. In the HR genes, we found that BRCA1/2 and PLAB2 mutations were related to HRD. HRD was also detected in a subset of patients without germline or somatic mutations in BRCA1/2, PLAB2, or other HR-related genes. Notably, LRP1B, NOTCH3, GATA2, and CARD11 (abbreviated as LNGC) mutations were associated with high HRD scores in breast cancer patients. Furthermore, functional experiments demonstrated that silencing CARD11 and GATA2 impairs HR repair efficiency and enhances the sensitivity of tumor cells to olaparib treatment. In summary, in the absence of mutations in the HR genes, the sensitivity of tumor cells to PARP inhibitors and platinum-based chemotherapy may be enhanced in a subset of breast cancer patients with LNGC somatic mutations.

Digital PCR quantification of ultrahigh ERBB2 copy number identifies poor breast cancer survival after trastuzumab

HER2/ERBB2 evaluation is necessary for treatment decision-making in breast cancer (BC), however current methods have limitations and considerable variability exists. DNA copy number (CN) evaluation by droplet digital PCR (ddPCR) has complementary advantages for HER2/ERBB2 diagnostics. In this study, we developed a single-reaction multiplex ddPCR assay for determination of ERBB2 CN in reference to two control regions, CEP17 and a copy-number-stable region of chr. 2p13.1, validated CN estimations to clinical in situ hybridization (ISH) HER2 status, and investigated the association of ERBB2 CN with clinical outcomes. 909 primary BC tissues were evaluated and the area under the curve for concordance to HER2 status was 0.93 and 0.96 for ERBB2 CN using either CEP17 or 2p13.1 as reference, respectively. The accuracy of ddPCR ERBB2 CN was 93.7% and 94.1% in the training and validation groups, respectively. Positive and negative predictive value for the classic HER2 amplification and non-amplification groups was 97.2% and 94.8%, respectively. An identified biological "ultrahigh" ERBB2 ddPCR CN group had significantly worse survival within patients treated with adjuvant trastuzumab for both recurrence-free survival (hazard ratio, HR: 3.3; 95% CI 1.1-9.6; p = 0.031, multivariable Cox regression) and overall survival (HR: 3.6; 95% CI 1.1-12.6; p = 0.041). For validation using RNA-seq data as a surrogate, in a population-based SCAN-B cohort (NCT02306096) of 682 consecutive patients receiving adjuvant trastuzumab, the ultrahigh-ERBB2 mRNA group had significantly worse survival. Multiplex ddPCR is useful for ERBB2 CN estimation and ultrahigh ERBB2 may be a predictive factor for decreased long-term survival after trastuzumab treatment.

EXO1 protects BRCA1-deficient cells against toxic DNA lesions

Inactivating mutations in the BRCA1 and BRCA2 genes impair DNA double-strand break (DSB) repair by homologous recombination (HR), leading to chromosomal instability and cancer. Importantly, BRCA1/2 deficiency also causes therapeutically targetable vulnerabilities. Here, we identify the dependency on the end resection factor EXO1 as a key vulnerability of BRCA1-deficient cells. EXO1 deficiency generates poly(ADP-ribose)-decorated DNA lesions during S phase that associate with unresolved DSBs and genomic instability in BRCA1-deficient but not in wild-type or BRCA2-deficient cells. Our data indicate that BRCA1/EXO1 double-deficient cells accumulate DSBs due to impaired repair by single-strand annealing (SSA) on top of their HR defect. In contrast, BRCA2-deficient cells retain SSA activity in the absence of EXO1 and hence tolerate EXO1 loss. Consistent with a dependency on EXO1-mediated SSA, we find that BRCA1-mutated tumors show elevated EXO1 expression and increased SSA-associated genomic scars compared with BRCA1-proficient tumors. Overall, our findings uncover EXO1 as a promising therapeutic target for BRCA1-deficient tumors.

Atezolizumab in Combination With Carboplatin and Survival Outcomes in Patients With Metastatic Triple-Negative Breast Cancer: The TBCRC 043 Phase 2 Randomized Clinical Trial

Importance

Agents targeting programmed death ligand 1 (PD-L1) have demonstrated efficacy in triple-negative breast cancer (TNBC) when combined with chemotherapy and are now the standard of care in patients with PD-L1-positive metastatic disease. In contrast to microtubule-targeting agents, the effect of combining platinum compounds with programmed cell death 1 (PD-1)/PD-L1 immunotherapy has not been extensively determined.

Objective

To evaluate the efficacy of atezolizumab with carboplatin in patients with metastatic TNBC.

Design, Setting, and Participants

This phase 2 randomized clinical trial was conducted in 6 centers from August 2017 to June 2021.

Interventions

Patients with metastatic TNBC were randomized to receive carboplatin area under the curve (AUC) 6 alone or with atezolizumab, 1200 mg, every 3 weeks until disease progression or unacceptable toxic effects with a 3-year duration of follow-up.

Main Outcome and Measures

The primary end point was investigator-assessed progression-free survival (PFS). Secondary end points included overall response rate (ORR), clinical benefit rate (CBR), and overall survival (OS). Other objectives included correlation of response with tumor PD-L1 levels, tumor-infiltrating lymphocytes (TILs), tumor DNA- and RNA-sequenced biomarkers, TNBC subtyping, and multiplex analyses of immune markers.

Results

All 106 patients with metastatic TNBC who were enrolled were female with a mean (range) age of 55 (27-79) years, of which 12 (19%) identified as African American/Black, 1 (1%) as Asian, 73 (69%) as White, and 11 (10%) as unknown. Patients were randomized and received either carboplatin (n = 50) or carboplatin and atezolizumab (n = 56). The combination improved PFS (hazard ratio [HR], 0.66; 95% CI, 0.44-1.01; P = .05) from a median of 2.2 to 4.1 months, increased ORR from 8.0% (95% CI, 3.2%-18.8%) to 30.4% (95% CI, 19.9%-43.3%), increased CBR at 6 months from 18.0% (95% CI, 9.8%-30.1%) to 37.5% (95% CI, 26.0%-50.6%), and improved OS (HR, 0.60; 95% CI, 0.37-0.96; P = .03) from a median of 8.6 to 12.6 months. Subgroup analysis showed PD-L1-positive tumors did not benefit more from adding atezolizumab (HR, 0.62; 95% CI, 0.23-1.65; P = .35). Patients with high TILs (HR, 0.12; 95% CI, 0.30-0.50), high mutation burden (HR, 0.50; 95% CI, 0.23-1.06), and prior chemotherapy (HR, 0.59; 95% CI, 0.36-0.95) received greater benefit on the combination. Patients with obesity and patients with more than 125 mg/dL on-treatment blood glucose levels were associated with better PFS (HR, 0.35; 95% CI, 0.10-1.80) on the combination. TNBC subtypes benefited from adding atezolizumab, except the luminal androgen receptor subtype.

Conclusions and Relevance

In this randomized clinical trial, the addition of atezolizumab to carboplatin significantly improved survival of patients with metastatic TNBC regardless of PD-L1 status. Further, lower risk of disease progression was associated with increased TILs, higher mutation burden, obesity, and uncontrolled blood glucose levels.

Trial Registration

ClinicalTrials.gov Identifier: NCT03206203.

Improved detection of clinically relevant fusion transcripts in cancer by machine learning classification

BACKGROUND

Genomic rearrangements in cancer cells can create fusion genes that encode chimeric proteins or alter the expression of coding and non-coding RNAs. In some cancer types, fusions involving specific kinases are used as targets for therapy. Fusion genes can be detected by whole genome sequencing (WGS) and targeted fusion panels, but RNA sequencing (RNA-Seq) has the advantageous capability of broadly detecting expressed fusion transcripts.

RESULTS

We developed a pipeline for validation of fusion transcripts identified in RNA-Seq data using matched WGS data from The Cancer Genome Atlas (TCGA) and applied it to 910 tumors from 11 different cancer types. This resulted in 4237 validated gene fusions, 3049 of them with at least one identified genomic breakpoint. Utilizing validated fusions as true positive events, we trained a machine learning classifier to predict true and false positive fusion transcripts from RNA-Seq data. The final precision and recall metrics of the classifier were 0.74 and 0.71, respectively, in an independent dataset of 249 breast tumors. Application of this classifier to all samples with RNA-Seq data from these cancer types vastly extended the number of likely true positive fusion transcripts and identified many potentially targetable kinase fusions. Further analysis of the validated gene fusions suggested that many are created by intrachromosomal amplification events with microhomology-mediated non-homologous end-joining.

CONCLUSIONS

A classifier trained on validated fusion events increased the accuracy of fusion transcript identification in samples without WGS data. This allowed the analysis to be extended to all samples with RNA-Seq data, facilitating studies of tumor biology and increasing the number of detected kinase fusions. Machine learning could thus be used in identification of clinically relevant fusion events for targeted therapy. The large dataset of validated gene fusions generated here presents a useful resource for development and evaluation of fusion transcript detection algorithms.

Olaparib and Ceralasertib (AZD6738) in Patients with Triple-Negative Advanced Breast Cancer: Results from Cohort E of the plasmaMATCH Trial (CRUK/15/010)

PURPOSE

Approximately 10% to 15% of triple-negative breast cancers (TNBC) have deleterious mutations in BRCA1 and BRCA2 and may benefit from PARP inhibitor treatment. PARP inhibitors may also increase exogenous replication stress and thereby increase sensitivity to inhibitors of ataxia telangiectasia and Rad3-related (ATR) protein. This phase II study examined the activity of the combination of PARP inhibitor, olaparib, and ATR inhibitor, ceralasertib (AZD6738), in patients with advanced TNBC.

PATIENTS AND METHODS

Patients with TNBC on most recent biopsy who had received 1 or 2 lines of chemotherapy for advanced disease or had relapsed within 12 months of (neo)adjuvant chemotherapy were eligible. Treatment was olaparib 300 mg twice a day continuously and celarasertib 160 mg on days 1-7 on a 28-day cycle until disease progression. The primary endpoint was confirmed objective response rate (ORR). Tissue and plasma biomarker analyses were preplanned to identify predictors of response.

RESULTS

70 evaluable patients were enrolled. Germline BRCA1/2 mutations were present in 10 (14%) patients and 3 (4%) patients had somatic BRCA mutations. The confirmed ORR was 12/70; 17.1% (95% confidence interval, 10.4-25.5). Responses were observed in patients without germline or somatic BRCA1/2 mutations, including patients with mutations in other homologous recombination repair genes and tumors with functional homologous recombination deficiency by RAD51 foci.

CONCLUSIONS

The response rate to olaparib and ceralasertib did not meet prespecified criteria for activity in the overall evaluable population, but responses were observed in patients who would not be expected to respond to olaparib monotherapy.

Chromatin insulation orchestrates matrix metalloproteinase gene cluster expression reprogramming in aggressive breast cancer tumors

BACKGROUND

Triple-negative breast cancer (TNBC) is an aggressive subtype that exhibits a high incidence of distant metastases and lacks targeted therapeutic options. Here we explored how the epigenome contributes to matrix metalloprotease (MMP) dysregulation impacting tumor invasion, which is the first step of the metastatic process.

METHODS

We combined RNA expression and chromatin interaction data to identify insulator elements potentially associated with MMP gene expression and invasion. We employed CRISPR/Cas9 to disrupt the CCCTC-Binding Factor (CTCF) binding site on an insulator element downstream of the MMP8 gene (IE8) in two TNBC cellular models. We characterized these models by combining Hi-C, ATAC-seq, and RNA-seq with functional experiments to determine invasive ability. The potential of our findings to predict the progression of ductal carcinoma in situ (DCIS), was tested in data from clinical specimens.

RESULTS

We explored the clinical relevance of an insulator element located within the Chr11q22.2 locus, downstream of the MMP8 gene (IE8). This regulatory element resulted in a topologically associating domain (TAD) boundary that isolated nine MMP genes into two anti-correlated expression clusters. This expression pattern was associated with worse relapse-free (HR = 1.57 [1.06 - 2.33]; p = 0.023) and overall (HR = 2.65 [1.31 - 5.37], p = 0.005) survival of TNBC patients. After CRISPR/Cas9-mediated disruption of IE8, cancer cells showed a switch in the MMP expression signature, specifically downregulating the pro-invasive MMP1 gene and upregulating the antitumorigenic MMP8 gene, resulting in reduced invasive ability and collagen degradation. We observed that the MMP expression pattern predicts DCIS that eventually progresses into invasive ductal carcinomas (AUC = 0.77, p < 0.01).

CONCLUSION

Our study demonstrates how the activation of an IE near the MMP8 gene determines the regional transcriptional regulation of MMP genes with opposing functional activity, ultimately influencing the invasive properties of aggressive forms of breast cancer.

The landscape of human SVA retrotransposons

SINE-VNTR-Alu (SVA) retrotransposons are evolutionarily young and still-active transposable elements (TEs) in the human genome. Several pathogenic SVA insertions have been identified that directly mutate host genes to cause neurodegenerative and other types of diseases. However, due to their sequence heterogeneity and complex structures as well as limitations in sequencing techniques and analysis, SVA insertions have been less well studied compared to other mobile element insertions. Here, we identified polymorphic SVA insertions from 3646 whole-genome sequencing (WGS) samples of >150 diverse populations and constructed a polymorphic SVA insertion reference catalog. Using 20 long-read samples, we also assembled reference and polymorphic SVA sequences and characterized the internal hexamer/variable-number-tandem-repeat (VNTR) expansions as well as differing SVA activity for SVA subfamilies and human populations. In addition, we developed a module to annotate both reference and polymorphic SVA copies. By characterizing the landscape of both reference and polymorphic SVA retrotransposons, our study enables more accurate genotyping of these elements and facilitate the discovery of pathogenic SVA insertions.

Low dose PFDA induces DNA damage and DNA repair inhibition by promoting nuclear cGAS accumulation in ovarian epithelial cells

Per- and polyfluoroalkyl substances (PFAS), the versatile anthropogenic chemicals, are popular with the markets and manufactured in large quantities yearly. Accumulation of PFAS has various adverse health effects on human. Albeit certain members of PFAS were found to have genotoxicity in previous studies, the mechanisms underlying their effects on DNA damage repair remain unclear. Here, we investigated the effects of Perfluorodecanoic acid (PFDA) on DNA damage and DNA damage repair in ovarian epithelial cells through a series of in vivo and in vitro experiments. At environmentally relevant concentration, we firstly found that PFDA can cause DNA damage in primary mouse ovarian epithelial cells and IOSE-80 cells. Moreover, nuclear cGAS increased in PFDA-treated cells, which leaded to the efficiency of DNA homologous recombination (HR) decreased and DNA double-strand breaks perpetuated. In vivo experiments also verified that PFDA can induce more DNA double-strand breaks lesions and nuclear cGAS in ovarian tissue. Taken together, our results unveiled that low dose PFDA can cause deleterious effects on DNA and DNA damage repair (DDR) in ovarian epithelial cells and induce genomic instability.

Comparison of SP142 and 22C3 PD-L1 assays in a population-based cohort of triple-negative breast cancer patients in the context of their clinically established scoring algorithms

BACKGROUND

Immunohistochemical (IHC) PD-L1 expression is commonly employed as predictive biomarker for checkpoint inhibitors in triple-negative breast cancer (TNBC). However, IHC evaluation methods are non-uniform and further studies are needed to optimize clinical utility.

METHODS

We compared the concordance, prognostic value and gene expression between PD-L1 IHC expression by SP142 immune cell (IC) score and 22C3 combined positive score (CPS; companion IHC diagnostic assays for atezolizumab and pembrolizumab, respectively) in a population-based cohort of 232 early-stage TNBC patients.

RESULTS

The expression rates of PD-L1 for SP142 IC ≥ 1%, 22C3 CPS ≥ 10, 22C3 CPS ≥ 1 and 22C3 IC ≥ 1% were 50.9%, 27.2%, 53.9% and 41.8%, respectively. The analytical concordance (kappa values) between SP142 IC+ and these three different 22C3 scorings were 73.7% (0.48, weak agreement), 81.5% (0.63) and 86.6% (0.73), respectively. The SP142 assay was better at identifying 22C3 positive tumors than the 22C3 assay was at detecting SP142 positive tumors. PD-L1 (CD274) gene expression (mRNA) showed a strong positive association with all two-categorical IHC scorings of the PD-L1 expression, irrespective of antibody and cut-off (Spearman Rho ranged from 0.59 to 0.62; all p-values < 0.001). PD-L1 IHC positivity and abundance of tumor infiltrating lymphocytes were of positive prognostic value in univariable regression analyses in patients treated with (neo)adjuvant chemotherapy, where it was strongest for 22C3 CPS ≥ 10 and distant relapse-free interval (HR = 0.18, p = 0.019). However, PD-L1 status was not independently prognostic when adjusting for abundance of tumor infiltrating lymphocytes in multivariable analyses.

CONCLUSION

Our findings support that the SP142 and 22C3 IHC assays, with their respective clinically applied scoring algorithms, are not analytically equivalent where they identify partially non-overlapping subpopulations of TNBC patients and cannot be substituted with one another regarding PD-L1 detection. Trial registration The Swedish Cancerome Analysis Network - Breast (SCAN-B) study, retrospectively registered 2nd Dec 2014 at ClinicalTrials.gov; ID NCT02306096.

Precision cancer medicine: Concepts, current practice, and future developments

Precision cancer medicine is a multidisciplinary team effort that requires involvement and commitment of many stakeholders including the society at large. Building on the success of significant advances in precision therapy for oncological patients over the last two decades, future developments will be significantly shaped by improvements in scalable molecular diagnostics in which increasingly complex multilayered datasets require transformation into clinically useful information guiding patient management at fast turnaround times. Adaptive profiling strategies involving tissue- and liquid-based testing that account for the immense plasticity of cancer during the patient's journey and also include early detection approaches are already finding their way into clinical routine and will become paramount. A second major driver is the development of smart clinical trials and trial concepts which, complemented by real-world evidence, rapidly broaden the spectrum of therapeutic options. Tight coordination with regulatory agencies and health technology assessment bodies is crucial in this context. Multicentric networks operating nationally and internationally are key in implementing precision oncology in clinical practice and support developing and improving the ecosystem and framework needed to turn invocation into benefits for patients. The review provides an overview of the diagnostic tools, innovative clinical studies, and collaborative efforts needed to realize precision cancer medicine.

Integrated Multimodal Analyses of DNA Damage Response and Immune Markers as Predictors of Response in Metastatic Triple-Negative Breast Cancer in the TNT Trial (NCT00532727)

PURPOSE

The TNT trial (NCT00532727) showed no evidence of carboplatin superiority over docetaxel in metastatic triple-negative breast cancer (mTNBC), but carboplatin benefit was observed in the germline BRCA1/2 mutation subgroup. Broader response-predictive biomarkers are needed. We explored the predictive ability of DNA damage response (DDR) and immune markers.

EXPERIMENTAL DESIGN

Tumor-infiltrating lymphocytes were evaluated for 222 of 376 patients. Primary tumors (PT) from 186 TNT participants (13 matched recurrences) were profiled using total RNA sequencing. Four transcriptional DDR-related and 25 immune-related signatures were evaluated. We assessed their association with objective response rate (ORR) and progression-free survival (PFS). Conditional inference forest clustering was applied to integrate multimodal data. The biology of subgroups was characterized by 693 gene expression modules and other markers.

RESULTS

Transcriptional DDR-related biomarkers were not predictive of ORR to either treatment overall. Changes from PT to recurrence were demonstrated; in chemotherapy-naïve patients, transcriptional DDR markers separated carboplatin responders from nonresponders (P values = 0.017; 0.046). High immune infiltration was associated with docetaxel ORR (interaction P values < 0.05). Six subgroups were identified; the immune-enriched cluster had preferential docetaxel response [62.5% (D) vs. 29.4% (C); P = 0.016]. The immune-depleted cluster had preferential carboplatin response [8.0% (D) vs. 40.0% (C); P = 0.011]. DDR-related subgroups were too small to assess ORR.

CONCLUSIONS

High immune features predict docetaxel response, and high DDR signature scores predict carboplatin response in treatment-naïve mTNBC. Integrating multimodal DDR and immune-related markers identifies subgroups with differential treatment sensitivity. Treatment options for patients with immune-low and DDR-proficient tumors remains an outstanding need. Caution is needed using PT-derived transcriptional signatures to direct treatment in mTNBC, particularly DDR-related markers following prior chemotherapy.

A novel marker integrating multiple genetic alterations better predicts platinum sensitivity in ovarian cancer than HRD score

Introduction: Platinum-based chemotherapy is the first-line treatment strategy for ovarian cancer patients. The dismal prognosis of ovarian cancer was shown to be stringently associated with the heterogeneity of tumor cells in response to this therapy, therefore understanding platinum sensitivity in ovarian cancer would be helpful for improving patients' quality of life and clinical outcomes. HRDetect, utilized to characterize patients' homologous recombination repair deficiency, was used to predict patients' response to platinum-based chemotherapy. However, whether each of the single features contributing to HRD score is associated with platinum sensitivity remains elusive. Methods: We analyzed the whole-exome sequencing data of 196 patients who received platinum-based chemotherapy from the TCGA database. Genetic features were determined individually to see if they could indicate patients' response to platinum-based chemotherapy and prognosis, then integrated into a Pt-score employing LASSO regression model to assess its predictive performance. Results and discussion: Multiple genetic features, including bi-allelic inactivation of BRCA1/2 genes and genes involved in HR pathway, multiple somatic mutations in genes involved in DNA damage repair (DDR), and previously reported HRD-related features, were found to be stringently associated with platinum sensitivity and improved prognosis. Higher contributions of mutational signature SBS39 or ID6 predicted improved overall survival. Besides, arm-level loss of heterozygosity (LOH) of either chr4p or chr5q predicted significantly better disease-free survival. Notably, some of these features were found independent of HRD. And SBS3, an HRD-related feature, was found irrelevant to platinum sensitivity. Integrated all candidate markers using the LASSO model to yield a Pt-score, which showed better predictive ability compared to HRDetect in determining platinum sensitivity and predicting patients' prognosis, and this performance was validated in an independent cohort. The outcomes of our study will be instrumental in devising effective strategies for treating ovarian cancer with platinum-based chemotherapy.

TEQUILA-seq: a versatile and low-cost method for targeted long-read RNA sequencing

Long-read RNA sequencing (RNA-seq) is a powerful technology for transcriptome analysis, but the relatively low throughput of current long-read sequencing platforms limits transcript coverage. One strategy for overcoming this bottleneck is targeted long-read RNA-seq for preselected gene panels. We present TEQUILA-seq, a versatile, easy-to-implement, and low-cost method for targeted long-read RNA-seq utilizing isothermally linear-amplified capture probes. When performed on the Oxford nanopore platform with multiple gene panels of varying sizes, TEQUILA-seq consistently and substantially enriches transcript coverage while preserving transcript quantification. We profile full-length transcript isoforms of 468 actionable cancer genes across 40 representative breast cancer cell lines. We identify transcript isoforms enriched in specific subtypes and discover novel transcript isoforms in extensively studied cancer genes such as TP53. Among cancer genes, tumor suppressor genes (TSGs) are significantly enriched for aberrant transcript isoforms targeted for degradation via mRNA nonsense-mediated decay, revealing a common RNA-associated mechanism for TSG inactivation. TEQUILA-seq reduces the per-reaction cost of targeted capture by 2-3 orders of magnitude, as compared to a standard commercial solution. TEQUILA-seq can be broadly used for targeted sequencing of full-length transcripts in diverse biomedical research settings.

Copy number alteration features in pan-cancer homologous recombination deficiency prediction and biology

Homologous recombination deficiency (HRD) renders cancer cells vulnerable to unrepaired double-strand breaks and is an important therapeutic target as exemplified by the clinical efficacy of poly ADP-ribose polymerase (PARP) inhibitors as well as the platinum chemotherapy drugs applied to HRD patients. However, it remains a challenge to predict HRD status precisely and economically. Copy number alteration (CNA), as a pervasive trait of human cancers, can be extracted from a variety of data sources, including whole genome sequencing (WGS), SNP array, and panel sequencing, and thus can be easily applied clinically. Here we systematically evaluate the predictive performance of various CNA features and signatures in HRD prediction and build a gradient boosting machine model (HRD_CNA) for pan-cancer HRD prediction based on these CNA features. CNA features BP10MB[1] (The number of breakpoints per 10MB of DNA is 1) and SS[ > 7 & <=8] (The log10-based size of segments is greater than 7 and less than or equal to 8) are identified as the most important features in HRD prediction. HRD_CNA suggests the biallelic inactivation of BRCA1, BRCA2, PALB2, RAD51C, RAD51D, and BARD1 as the major genetic basis for human HRD, and may also be applied to effectively validate the pathogenicity of BRCA1/2 variants of uncertain significance (VUS). Together, this study provides a robust tool for cost-effective HRD prediction and also demonstrates the applicability of CNA features and signatures in cancer precision medicine.

Molecular characteristics of breast tumors in patients screened for germline predisposition from a population-based observational study

BACKGROUND

Pathogenic germline variants (PGVs) in certain genes are linked to higher lifetime risk of developing breast cancer and can influence preventive surgery decisions and therapy choices. Public health programs offer genetic screening based on criteria designed to assess personal risk and identify individuals more likely to carry PGVs, dividing patients into screened and non-screened groups. How tumor biology and clinicopathological characteristics differ between these groups is understudied and could guide refinement of screening criteria.

METHODS

Six thousand six hundred sixty breast cancer patients diagnosed in South Sweden during 2010-2018 were included with available clinicopathological and RNA sequencing data, 900 (13.5%) of which had genes screened for PGVs through routine clinical screening programs. We compared characteristics of screened patients and tumors to non-screened patients, as well as between screened patients with (n = 124) and without (n = 776) PGVs.

RESULTS

Broadly, breast tumors in screened patients showed features of a more aggressive disease. However, few differences related to tumor biology or patient outcome remained significant after stratification by clinical subgroups or PAM50 subtypes. Triple-negative breast cancer (TNBC), the subgroup most enriched for PGVs, showed the most differences between screening subpopulations (e.g., higher tumor proliferation in screened cases). Significant differences in PGV prevalence were found between clinical subgroups/molecular subtypes, e.g., TNBC cases were enriched for BRCA1 PGVs. In general, clinicopathological differences between screened and non-screened patients mimicked those between patients with and without PGVs, e.g., younger age at diagnosis for positive cases. However, differences in tumor biology/microenvironment such as immune cell composition were additionally seen within PGV carriers/non-carriers in ER + /HER2 - cases, but not between screening subpopulations in this subgroup.

CONCLUSIONS

Characterization of molecular tumor features in patients clinically screened and not screened for PGVs represents a relevant read-out of guideline criteria. The general lack of molecular differences between screened/non-screened patients after stratification by relevant breast cancer subsets questions the ability to improve the identification of screening candidates based on currently used patient and tumor characteristics, pointing us towards universal screening. Nevertheless, while that is not attained, molecular differences identified between PGV carriers/non-carriers suggest the possibility of further refining patient selection within certain patient subsets using RNA-seq through, e.g., gene signatures.

TRIAL REGISTRATION

The Sweden Cancerome Analysis Network - Breast (SCAN-B) was prospectively registered at ClinicalTrials.gov under the identifier NCT02306096.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

BRCA1-methylated triple negative breast cancers previously exposed to neoadjuvant chemotherapy form RAD51 foci and respond poorly to olaparib

Background

About 15% of Triple-Negative-Breast-Cancer (TNBC) present silencing of the BRCA1 promoter methylation and are assumed to be Homologous Recombination Deficient (HRD). BRCA1-methylated (BRCA1-Me) TNBC could, thus, be eligible to treatment based on PARP-inhibitors or Platinum salts. However, their actual HRD status is discussed, as these tumors are suspected to develop resistance after chemotherapy exposure.

Methods

We interrogated the sensitivity to olaparib vs. carboplatin of 8 TNBC Patient-Derived Xenografts (PDX) models. Four PDX corresponded to BRCA1-Me, of which 3 were previously exposed to NeoAdjuvant-Chemotherapy (NACT). The remaining PDX models corresponded to two BRCA1-mutated (BRCA1-Mut) and two BRCA1-wild type PDX that were respectively included as positive and negative controls. The HRD status of our PDX models was assessed using both genomic signatures and the functional BRCA1 and RAD51 nuclear foci formation assay. To assess HR restoration associated with olaparib resistance, we studied pairs of BRCA1 deficient cell lines and their resistant subclones.

Results

The 3 BRCA1-Me PDX that had been exposed to NACT responded poorly to olaparib, likewise BRCA1-WT PDX. Contrastingly, 3 treatment-naïve BRCA1-deficient PDX (1 BRCA1-Me and 2 BRCA1-mutated) responded to olaparib. Noticeably, the three olaparib-responsive PDX scored negative for BRCA1- and RAD51-foci, whereas all non-responsive PDX models, including the 3 NACT-exposed BRCA1-Me PDX, scored positive for RAD51-foci. This suggested HRD in olaparib responsive PDX, while non-responsive models were HR proficient. These results were consistent with observations in cell lines showing a significant increase of RAD51-foci in olaparib-resistant subclones compared with sensitive parental cells, suggesting HR restoration in these models.

Conclusion

Our results thus support the notion that the actual HRD status of BRCA1-Me TNBC, especially if previously exposed to chemotherapy, may be questioned and should be verified using the BRCA1- and RAD51-foci assay.

Artificial Diets with Selective Restriction of Amino Acids and Very Low Levels of Lipids Induce Anticancer Activity in Mice with Metastatic Triple-Negative Breast Cancer

Patients with metastatic triple negative breast cancer (TNBC) need new therapies to improve the low survival rates achieved with standard treatments. In this work, we show for the first time that the survival of mice with metastatic TNBC can be markedly increased by replacing their normal diet with artificial diets in which the levels of amino acids (AAs) and lipids are strongly manipulated. After observing selective anticancer activity in vitro, we prepared five artificial diets and evaluated their anticancer activity in a challenging model of metastatic TNBC. The model was established by injecting 4T1 murine TNBC cells into the tail vein of immunocompetent BALB/cAnNRj mice. First-line drugs doxorubicin and capecitabine were also tested in this model. AA manipulation led to modest improvements in mice survival when the levels of lipids were normal. Reducing lipid levels to 1% markedly improved the activity of several diets with different AA content. Some mice fed the artificial diets as monotherapy lived much longer than mice treated with doxorubicin and capecitabine. An artificial diet without 10 non-essential AAs, with reduced levels of essential AAs, and with 1% lipids improved the survival not only of mice with TNBC but also of mice with other types of metastatic cancers.

EXO1-mediated DNA repair by single-strand annealing is essential for BRCA1-deficient cells

Deficiency for the repair of DNA double-strand breaks (DSBs) via homologous recombination (HR) leads to chromosomal instability and diseases such as cancer. Yet, defective HR also results in vulnerabilities that can be exploited for targeted therapy. Here, we identify such a vulnerability and show that BRCA1-deficient cells are dependent on the long-range end-resection factor EXO1 for survival. EXO1 loss results in DNA replication-induced lesions decorated by poly(ADP-ribose)-chains. In cells that lack both BRCA1 and EXO1, this is accompanied by unresolved DSBs due to impaired single-strand annealing (SSA), a DSB repair process that requires the activity of both proteins. In contrast, BRCA2-deficient cells have increased SSA, also in the absence of EXO1, and hence are not dependent on EXO1 for survival. In agreement with our mechanistic data, BRCA1-mutated tumours have elevated EXO1 expression and contain more genomic signatures of SSA compared to BRCA1-proficient tumours. Collectively, our data indicate that EXO1 is a promising novel target for treatment of BRCA1-deficient tumours.

Systemic Therapy for Hereditary Breast Cancers

Approximately 5% to 10% of all breast cancers are hereditary; many of which are caused by pathogenic variants in genes required for homologous recombination, including BRCA1 and BRCA2. Here we discuss systemic treatment for such breast cancers, including approved chemotherapeutic approaches and also targeted treatment approaches using poly-(ADP ribose) polymerase inhibitors. We also discuss experimental approaches to treating hereditary breast cancer, including new small molecule DNA repair inhibitors and also immunomodulatory agents. Finally, we discuss how drug resistance emerges in patients with hereditary breast cancer, how this might be delayed or prevented, and how biomarker-adapted treatment is molding the future management of hereditary breast cancer.

Retrotransposon insertion as a novel mutational cause of spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder resulting from biallelic alterations of the SMN1 gene: deletion, gene conversion or, in rare cases, intragenic variants. The disease severity is mainly influenced by the copy number of SMN2, a nearly identical gene, which produces only low amounts of full-length (FL) mRNA. Here we describe the first example of retrotransposon insertion as a pathogenic SMN1 mutational event. The 50-year-old patient is clinically affected by SMA type III with a diagnostic odyssey spanning nearly 30 years. Despite a mild disease course, he carries a single SMN2 copy. Using Exome Sequencing and Sanger sequencing, we characterized a SINE-VNTR-Alu (SVA) type F retrotransposon inserted in SMN1 intron 7. Using RT-PCR and RNASeq experiments on lymphoblastoid cell lines, we documented the dramatic decrease of FL transcript production in the patient compared to subjects with the same SMN1 and SMN2 copy number, thus validating the pathogenicity of this SVA insertion. We described the mutant FL-SMN1-SVA transcript characterized by exon extension and showed that it is subject to degradation by nonsense-mediated mRNA decay. The stability of the SMN-SVA protein may explain the mild course of the disease. This observation exemplifies the role of retrotransposons in human genetic disorders.

Development of Homologous Recombination Functional Assays for Targeting the DDR

Identification of tumours that have homologous recombination deficiency (HRD) has become of increasing interest following the licensing of PARP inhibitors. Potential methods to assess HRD status include; clinical selection for platinum sensitive disease, mutational/methylation status, genomic scars/signature and functional RAD51 assays. Homologous recombination (HR) is a dynamic process with the potential to evolve over a disease course, particularly in relation to previous treatment. This is one of the major drawbacks of genomic scars/signatures, as they only demonstrate historic HR status. Functional HR assays have the benefit of giving a real time HR status readout and therefore have the potential for clearer identification of patients who may benefit from PARP inhibitors at that specific time point. However, the development of RAD51 foci assays ready for clinical practice has been challenging. Pre-clinical considerations have included; controlling for variation in tumour proliferation, tissue type and whether DNA damage induction is required. Furthermore, the assays require correlation with clinical outcomes, an understanding of how they complement current testing modalities and validation of test performance in large cohorts. Despite these challenges, given the profound benefit from PARP inhibitors seen in those with an HRD phenotype to date, the ongoing development and validation of these functional HR assays remains of high clinical importance.

GenoREC: A Recommendation System for Interactive Genomics Data Visualization

Interpretation of genomics data is critically reliant on the application of a wide range of visualization tools. A large number of visualization techniques for genomics data and different analysis tasks pose a significant challenge for analysts: which visualization technique is most likely to help them generate insights into their data? Since genomics analysts typically have limited training in data visualization, their choices are often based on trial and error or guided by technical details, such as data formats that a specific tool can load. This approach prevents them from making effective visualization choices for the many combinations of data types and analysis questions they encounter in their work. Visualization recommendation systems assist non-experts in creating data visualization by recommending appropriate visualizations based on the data and task characteristics. However, existing visualization recommendation systems are not designed to handle domain-specific problems. To address these challenges, we designed GenoREC, a novel visualization recommendation system for genomics. GenoREC enables genomics analysts to select effective visualizations based on a description of their data and analysis tasks. Here, we present the recommendation model that uses a knowledge-based method for choosing appropriate visualizations and a web application that enables analysts to input their requirements, explore recommended visualizations, and export them for their usage. Furthermore, we present the results of two user studies demonstrating that GenoREC recommends visualizations that are both accepted by domain experts and suited to address the given genomics analysis problem. All supplemental materials are available at https://osf.io/y73pt/.

Update on triple-negative breast cancers - highlighting subtyping update and treatment implication

Triple-negative breast cancer (TNBC) remains a major challenge in breast cancer management. Continuing research in the past years aimed at understanding the biology of this tumour and developing more effective therapeutic options. It is now clear that TNBC is vastly heterogeneous with diverse histological, molecular, immunological profiles and clinical differences. Current evidence suggested the existence of at least four predominant subtypes based on expression profiling across studies. These subtypes exhibited specific genomic alterations and tumour microenvironment. Subtype-specific therapeutic strategies were identified. Recognising these subtypes allows not only an improved prognostication but also a better treatment decision. Herein, we provide an overview of the recent findings on TNBC heterogeneity at different levels and corresponding subtyping. The characteristic of subtypes and the implication of these subtypings in therapeutic approaches are also discussed.

Preclinical evaluation of Insulin-like growth factor receptor 1 (IGF1R) and Insulin Receptor (IR) as a therapeutic targets in triple negative breast cancer

Triple Negative Breast Cancer (TNBC), a subtype of breast cancer, has fewer successful therapeutic therapies than other types of breast cancer. Insulin-like growth factor receptor 1 (IGF1R) and the Insulin receptor (IR) are associated with poor outcomes in TNBC. Targeting IGF1R has failed clinically. We aimed to test if inhibiting both IR/IGF1R was a rationale therapeutic approach to treat TNBC. We showed that despite IGF1R and IR being expressed in TNBC, their expression is not associated with a negative survival outcome. Furthermore, targeting both IR/IGF1R with inhibitors in multiple TNBC cell lines did not inhibit cell growth. Linsitinib, a small molecule inhibitor of both IGF1R and IR, did not block tumour formation and had no effect on tumour growth in vivo. Cumulatively these data suggest that while IGF1R and IR are expressed in TNBC, they are not good therapeutic targets. A potential reason for the limited anti-cancer impact when IR/IGF1R was targeted may be because multiple signalling pathways are altered in TNBC. Therefore, targeting individual signalling pathways may not be sufficient to inhibit cancer growth.

Optimized detection of homologous recombination deficiency improves the prediction of clinical outcomes in cancer

Homologous recombination DNA-repair deficiency (HRD) is a common driver of genomic instability and confers a therapeutic vulnerability in cancer. The accurate detection of somatic allelic imbalances (AIs) has been limited by methods focused on BRCA1/2 mutations and using mixtures of cancer types. Using pan-cancer data, we revealed distinct patterns of AIs in high-grade serous ovarian cancer (HGSC). We used machine learning and statistics to generate improved criteria to identify HRD in HGSC (ovaHRDscar). ovaHRDscar significantly predicted clinical outcomes in three independent patient cohorts with higher precision than previous methods. Characterization of 98 spatiotemporally distinct metastatic samples revealed low intra-patient variation and indicated the primary tumor as the preferred site for clinical sampling in HGSC. Further, our approach improved the prediction of clinical outcomes in triple-negative breast cancer (tnbcHRDscar), validated in two independent patient cohorts. In conclusion, our tumor-specific, systematic approach has the potential to improve patient selection for HR-targeted therapies.

Single-cell genomic variation induced by mutational processes in cancer

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Genome Protection by DNA Polymerase θ

DNA polymerase θ (Pol θ) is a DNA repair enzyme widely conserved in animals and plants. Pol θ uses short DNA sequence homologies to initiate repair of double-strand breaks by theta-mediated end joining. The DNA polymerase domain of Pol θ is at the C terminus and is connected to an N-terminal DNA helicase-like domain by a central linker. Pol θ is crucial for maintenance of damaged genomes during development, protects DNA against extensive deletions, and limits loss of heterozygosity. The cost of using Pol θ for genome protection is that a few nucleotides are usually deleted or added at the repair site. Inactivation of Pol θ often enhances the sensitivity of cells to DNA strand-breaking chemicals and radiation. Since some homologous recombination-defective cancers depend on Pol θ for growth, inhibitors of Pol θ may be useful in treating such tumors.

Aberrant hyperexpression of the RNA binding protein FMRP in tumors mediates immune evasion

Many human cancers manifest the capability to circumvent attack by the adaptive immune system. In this work, we identified a component of immune evasion that involves frequent up-regulation of fragile X mental retardation protein (FMRP) in solid tumors. FMRP represses immune attack, as revealed by cancer cells engineered to lack its expression. FMRP-deficient tumors were infiltrated by activated T cells that impaired tumor growth and enhanced survival in mice. Mechanistically, FMRP’s immunosuppression was multifactorial, involving repression of the chemoattractant C-C motif chemokine ligand 7 (CCL7) concomitant with up-regulation of three immunomodulators—interleukin-33 (IL-33), tumor-secreted protein S (PROS1), and extracellular vesicles. Gene signatures associate FMRP’s cancer network with poor prognosis and response to therapy in cancer patients. Collectively, FMRP is implicated as a regulator that orchestrates a multifaceted barrier to antitumor immune responses.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications.

SIGNIFICANCE

Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

The triple negative breast cancer drugs graveyard: a review of failed clinical trials 2017-2022

INTRODUCTION

Triple-negative breast cancer (TNBC) accounts for 15-20% of breast cancers (BC) and has the worst prognosis. It is characterized by the absence of both hormone receptor (HR) and human epidermal growth factor receptor 2 (HER2). TNBC has more limited therapeutic options compared to other subtypes, meaning that there is still a long way to go to discover target treatments.

AREAS COVERED

Our review aims to summarize phase II/III clinical trials enrolling patients with TNBC that have been published between 2017 and 2022 but failed to reach their primary endpoint. We here try to emphasize the limitations and weaknesses noted in negative studies and to point out unexpected results which might be useful to enhance the therapeutic approach to TNBC disease.

EXPERT OPINION

A deeper understanding of the mechanisms behind TNBC heterogeneity allowed to enhance the knowledge of new prognostic and predictive biomarkers of response. However, it is also through several failed clinical trials that we were able to define new therapeutic approaches which improved TNBC patients' clinical outcomes. Nowadays, we still need to overcome several difficulties to fully recognize different intracellular and extracellular pathways that crosstalk in TNBC and the mechanisms of resistance to identify novel tailored-patients' therapies.

Two subtypes of cutaneous melanoma with distinct mutational signatures and clinico-genomic characteristics

Background: To decipher mutational signatures and their associations with biological implications in cutaneous melanomas (CMs), including those with a low ultraviolet (UV) signature.

Materials and Methods: We applied non-negative matrix factorization (NMF) and unsupervised clustering to the 96-class mutational context of The Cancer Genome Atlas (TCGA) cohort (N = 466) as well as other publicly available datasets (N = 527). To explore the feasibility of mutational signature-based classification using panel sequencing data, independent panel sequencing data were analyzed.

Results: NMF decomposition of the TCGA cohort and other publicly available datasets consistently found two mutational signatures: UV (SBS7a/7b dominant) and non-UV (SBS1/5 dominant) signatures. Based on mutational signatures, TCGA CMs were classified into two clusters: UV-high and UV-low. CMs belonging to the UV-low cluster showed significantly worse overall survival and landmark survival at 1-year than those in the UV-high cluster; low or high UV signature remained the most significant prognostic factor in multivariate analysis. The UV-low cluster showed distinct genomic and functional characteristic patterns: low mutation counts, increased proportion of triple wild-type and KIT mutations, high burden of copy number alteration, expression of genes related to keratinocyte differentiation, and low activation of tumor immunity. We verified that UV-high and UV-low clusters can be distinguished by panel sequencing.

Conclusion: Our study revealed two mutational signatures of CMs that divide CMs into two clusters with distinct clinico-genomic characteristics. Our results will be helpful for the clinical application of mutational signature-based classification of CMs.

Tumor purity adjusted beta values improve biological interpretability of high-dimensional DNA methylation data

A common issue affecting DNA methylation analysis in tumor tissue is the presence of a substantial amount of non-tumor methylation signal derived from the surrounding microenvironment. Although approaches for quantifying and correcting for the infiltration component have been proposed previously, we believe these have not fully addressed the issue in a comprehensive and universally applicable way. We present a multi-population framework for adjusting DNA methylation beta values on the Illumina 450/850K platform using generic purity estimates to account for non-tumor signal. Our approach also provides an indirect estimate of the aggregate methylation state of the surrounding normal tissue. Using whole exome sequencing derived purity estimates and Illumina 450K methylation array data generated by The Cancer Genome Atlas project (TCGA), we provide a demonstration of this framework in breast cancer illustrating the effect of beta correction on the aggregate methylation beta value distribution, clustering accuracy, and global methylation profiles.

Efficacy of different neoadjuvant treatment regimens in BRCA-mutated triple negative breast cancer: a systematic review and meta-analysis

Purpose

Triple negative breast cancer (TNBC) is an aggressive breast cancer strongly associated with BRCA mutation. Standard neoadjuvant chemotherapy remains the standard of care for early stage TNBC, the optimal chemotherapy regimen is still a matter of discussion. Other agents, such as poly-ADP-ribosyl polymerase inhibitors (PARPi) and anti-vascular endothelial growth factor (VEGF) antibodies were evaluated in the neoadjuvant setting. This systematic review and meta-analysis intend to evaluate the impact of neoadjuvant treatments in pCR rates in TNBC gBRCA mutation, beyond traditional standard chemotherapy.

Methods

PubMed, Clinicaltrials.gov, Cochrane CENTRAL, Embase and key oncological meetings for trials were searched for studies reporting neoadjuvant chemo-immunotherapy in BRCA positive TNBC.

Results

Out of 1238 records reviewed, thirty-one trials were included, resulting in a total 619 BRCA-mutated TNBC patients. In BRCA mutated TNBC patients who received cisplatin in monotherapy the proportion of patients who achieved pCR was 0.53 (95%CI [0.30, 0.76]), and when treatment combined standard chemotherapy and platin derivatives the proportion of pCR increased to 0.62 (95% CI [0.48, 0.76]). The group of patients treated with platin derivatives, anthracyclines ± taxanes achieved the highest proportion of pCR, 0.66. Patients treated with PARPi alone show a pCR proportion of 0.55 (95% CI [0.30, 0.81]); and when standard chemotherapy and platin derivatives were combined with PARPi the proportion of pCR did not vary.

Conclusions

Patients with BRCA mutated TNBC treated with cisplatin in monotherapy demonstrate inferior proportion in the pCR achievement when compared with standard chemotherapy plus platin derivates. The best pCR was achieved with platin derivates in association with anthracyclines ± taxanes. No difference in pCR was found between PARPi alone vs PARPi with standard chemotherapy.

Discordance in receptor status between primary and metastatic breast cancer and overall survival: A single-center analysis

BACKGROUND

The tumor phenotype may change between primary and metastatic breast cancer. We compared the expression of estrogen receptor (ER), progesterone receptor (PR), and HER2 in a series of primary breast carcinomas (PBC) with their metastatic relapses and analyzed the impact of any changes on survival.

MATERIALS AND METHODS

It was a single-center retrospective study, collecting consecutive cases of metastatic breast carcinoma diagnosed in the pathology and medical oncology departments at Habib Bourguiba University Hospital in Sfax, Tunisia. An immunohistochemical study was used to assess ER, PR, and HER2 expression. Overall survival (OS) and post-metastasis survival (PMS) were evaluated using multivariable Cox regression analysis.

RESULTS

Our study included 68 patients. ER and PR status changed in 29.4 % and 39.7 % of cases, respectively. Conversions were mainly from positive to negative status (22 % and 23.5 % for ER and PR, respectively). Differences in HER2 status were observed in 19.6 % of cases, with loss of overexpression in 6 patients (10.7 %). Adjuvant trastuzumab therapy and PBC molecular subtype (HR-, HER2+) were associated with HER2 status discordance (p = 0.02 and 0.03, respectively). On multivariable analysis, HR-negative conversion tumors were significantly associated with a worse OS (p = 0.042) and PMS (p < 0.001), compared to HR-concordant positive tumors.

CONCLUSION

This study establishes that HR and HER2 status discordance between primary and metastatic breast carcinoma has a prognostic impact on patient outcome. Analyzing these receptors' status in all newly diagnosed cases of metastatic breast carcinoma is strongly recommended and would provide information for changing treatment strategies.

Targeting the DNA Damage Response Pathways and Replication Stress in Colorectal Cancer

PURPOSE

Genomic instability is a hallmark of cancer and targeting DNA damage response (DDR) is emerging as a promising therapeutic strategy in different solid tumors. The effectiveness of targeting DDR in colorectal cancer has not been extensively explored.

EXPERIMENTAL DESIGN

We challenged 112 cell models recapitulating the genomic landscape of metastatic colorectal cancer with ATM, ATR, CHK1, WEE1, and DNA-PK inhibitors, in parallel with chemotherapeutic agents. We focused then on ATR inhibitors (ATRi) and, to identify putative biomarkers of response and resistance, we analyzed at multiple levels colorectal cancer models highly sensitive or resistant to these drugs.

RESULTS

We found that around 30% of colorectal cancers, including those carrying KRAS and BRAF mutations and unresponsive to targeted agents, are sensitive to at least one DDR inhibitor. By investigating potential biomarkers of response to ATRi, we found that ATRi-sensitive cells displayed reduced phospho-RPA32 foci at basal level, while ATRi-resistant cells showed increased RAD51 foci formation in response to replication stress. Lack of ATM and RAD51C expression was associated with ATRi sensitivity. Analysis of mutational signatures and HRDetect score identified a subgroup of ATRi-sensitive models. Organoids derived from patients with metastatic colorectal cancer recapitulated findings obtained in cell lines.

CONCLUSIONS

In conclusion, a subset of colorectal cancers refractory to current therapies could benefit from inhibitors of DDR pathways and replication stress. A composite biomarker involving phospho-RPA32 and RAD51 foci, lack of ATM and RAD51C expression, as well as analysis of mutational signatures could be used to identify colorectal cancers likely to respond to ATRi.