Mutational Signatures in Breast Cancer: The Problem at the DNA Level

Plasma THBS1 as a predictive biomarker for poor prognosis and brain metastasis in patients with HER2-enriched breast cancer

BACKGROUND

Thrombospondin-1 (THBS1) is a secretory adhesive glycoprotein involved in the progression of multiple malignancies, including breast cancer. However, the clinical significance and prognostic role of plasma THBS1 in breast cancer have yet to be clarified.

METHODS

Plasma THBS1 levels in 627 breast cancer patients were analyzed by enzyme-linked immunosorbent assay. Bone marrow blood was drawn from the anterior/posterior superior iliac spine to detect the presence of disseminated tumor cells (DTCs). The effects of plasma THBS1 on the clinicopathological characteristics and survival prediction of breast cancer patients were explored.

RESULTS

Plasma THBS1 did not correlate with overall survival, breast cancer-specific survival (BCSS), and distant disease-free survival (DDFS) in the entire breast cancer cohort. Notably, HER2-enriched patients with high-plasma THBS1 levels had significantly shorter BCSS (P = 0.027) and DDFS (P = 0.011) than those with low levels. Multivariate analyses revealed that plasma THBS1 was an independent prognostic marker of BCSS (P = 0.026) and DDFS (P = 0.007) in HER2-enriched patients. THBS1 levels were 24% higher in positive DTC patients than in negative DTC patients (P = 0.031), and high levels were significantly associated with poor BCSS in positive DTC patients (HR 2.08, 95% CI 1.17-3.71; P = 0.019). Moreover, high-plasma THBS1 levels were specifically associated with an increased occurrence of brain metastasis in HER2-enriched patients (P = 0.041).

CONCLUSION

These findings suggest that plasma THBS1 may be serving as an unfavorable prognosis predictor for HER2-enriched breast cancer and justifies the need for further research.

Nonlobular Invasive Breast Carcinomas with Biallelic Pathogenic CDH1 Somatic Alterations: A Histologic, Immunophenotypic, and Genomic Characterization

CDH1 encodes for E-cadherin, and its loss of function is the hallmark of invasive lobular carcinoma (ILC). Albeit vanishingly rare, biallelic CDH1 alterations may be found in nonlobular breast carcinomas (NL-BCs). We sought to determine the clinicopathologic characteristics and repertoire of genetic alterations of NL-BCs harboring CDH1 biallelic genetic alterations. Analysis of 5842 breast cancers (BCs) subjected to clinical tumor-normal sequencing with an FDA-cleared multigene panel was conducted to identify BCs with biallelic CDH1 pathogenic/likely pathogenic somatic mutations lacking lobular features. The genomic profiles of NL-BCs with CDH1 biallelic genetic alterations were compared with those of ILCs and invasive ductal carcinomas (IDCs), matched by clinicopathologic characteristics. Of the 896 CDH1-altered BCs, 889 samples were excluded based on the diagnosis of invasive mixed ductal/lobular carcinoma or ILC or the detection of monoallelic CDH1 alterations. Only 7 of the 5842 (0.11%) BCs harbored biallelic CDH1 alterations and lacked lobular features. Of these, 4/7 (57%) cases were ER-positive/HER2-negative, 1/7 (14%) was ER-positive/HER2-positive, and 2/7 (29%) were ER-negative/HER2-negative. In total, 5/7 (71%) were of Nottingham grade 2, and 2/7 (29%) were of grade 3. The NL-BCs with CDH1 biallelic genetic alterations included a mucinous carcinoma (n = 1), IDCs with focal nested growth (n = 2), IDC with solid papillary (n = 1) or apocrine (n = 2) features, and an IDC of no special type (NST; n = 1). E-cadherin expression, as detected by immunohistochemistry, was absent (3/5) or aberrant (discontinuous membranous/cytoplasmic/granular; 2/5). However, NL-BCs with CDH1 biallelic genetic alterations displayed recurrent genetic alterations, including TP53, PIK3CA (57%, 4/7; each), FGFR1, and NCOR1 (28%, 2/7, each) alterations. Compared with CDH1 wild-type IDC-NSTs, NL-BCs less frequently harbored GATA3 mutations (0% vs 47%, P = .03), but no significant differences were detected when compared with matched ILCs. Therefore, NL-BCs with CDH1 biallelic genetic alterations are vanishingly rare, predominantly comprise IDCs with special histologic features, and have genomic features akin to luminal B ER-positive BCs.

Non-homogeneous Poisson and renewal processes as spatial models for cancer mutation

Advances in sequencing technology assisted biologists in revealing signatures of DNA cancer mutation process and in demonstrating the mutagenesis behind. However, most of these signatures proposed by majority of work focus only on the type and frequency of mutations, without considering spatial information which is non-negligible in exploring mechanisms of mutation occurrence, e.g., Kataegis. Statistical characterization of the distance between consecutive mutations can give us relative spatial information; however, it ignores location information which is as important as distance information. In this work, we assume that DNA cancer mutations are location-dependent and that integrating the two variables, location and inter-distance, is beneficial to study DNA cancer mutation processes more accurately. Particularly, instead of following a specific distribution over the whole DNA sequence, we found out that the distribution of distance between successive mutations alternates between exponential and power-law distributions. Apart from this, the parameters of either of the two distributions vary with DNA locations. The cancers with kataegis phenomenon, a specific mutation pattern caused by abnormal activity of APOBEC protein family, are more likely to be accompanied by higher parameter values of distance distribution, implying higher occurrence rate of mutation. Therefore, we propose non-homogeneous Poisson and Renewal processes to spatially model DNA cancer mutations and to describe mutation patterns quantitatively and more accurately through a statistical perspective.

Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

APOBEC3 activity promotes the survival and evolution of drug-tolerant persister cells during acquired resistance to EGFR inhibitors in lung cancer

APOBEC mutagenesis is one of the most common endogenous sources of mutations in human cancer and is a major source of genetic intratumor heterogeneity. High levels of APOBEC mutagenesis are associated with poor prognosis and aggressive disease across diverse cancers, but the mechanistic and functional impacts of APOBEC mutagenesis on tumor evolution and therapy resistance remain relatively unexplored. To address this, we investigated the contribution of APOBEC mutagenesis to acquired therapy resistance in a model of EGFR-mutant non-small cell lung cancer. We find that inhibition of EGFR in lung cancer cells leads to a rapid and pronounced induction of APOBEC3 expression and activity. Functionally, APOBEC expression promotes the survival of drug-tolerant persister cells (DTPs) following EGFR inhibition. Constitutive expression of APOBEC3B alters the evolutionary trajectory of acquired resistance to the EGFR inhibitor gefitinib, making it more likely that resistance arises through de novo acquisition of the T790M gatekeeper mutation and squamous transdifferentiation during the DTP state. APOBEC3B expression is associated with increased expression of the squamous cell transcription factor ΔNp63 and squamous cell transdifferentiation in gefitinib-resistant cells. Knockout of ΔNp63 in gefitinibresistant cells reduces the expression of the p63 target genes IL1a/b and sensitizes these cells to the thirdgeneration EGFR inhibitor osimertinib. These results suggest that APOBEC activity promotes acquired resistance by facilitating evolution and transdifferentiation in DTPs, and suggest that approaches to target ΔNp63 in gefitinib-resistant lung cancers may have therapeutic benefit.

Clinical Implications of APOBEC3-Mediated Mutagenesis in Breast Cancer

Over recent years, members of the APOBEC3 family of cytosine deaminases have been implicated in increased cancer genome mutagenesis, thereby contributing to intratumor and intertumor genomic heterogeneity and therapy resistance in, among others, breast cancer. Understanding the available methods for clinical detection of these enzymes, the conditions required for their (dysregulated) expression, the clinical impact they have, and the clinical implications they may offer is crucial in understanding the current impact of APOBEC3-mediated mutagenesis in breast cancer. Here, we provide a comprehensive review of recent developments in the detection of APOBEC3-mediated mutagenesis and responsible APOBEC3 enzymes, summarize the pathways that control their expression, and explore the clinical ramifications and opportunities they pose. We propose that APOBEC3-mediated mutagenesis can function as a helpful predictive biomarker in several standard-of-care breast cancer treatment plans and may be a novel target for treatment.

Unexpected moves: a conformational change in MutSα enables high-affinity DNA mismatch binding

The DNA mismatch repair protein MutSα recognizes wrongly incorporated DNA bases and initiates their correction during DNA replication. Dysfunctions in mismatch repair lead to a predisposition to cancer. Here, we study the homozygous mutation V63E in MSH2 that was found in the germline of a patient with suspected constitutional mismatch repair deficiency syndrome who developed colorectal cancer before the age of 30. Characterization of the mutant in mouse models, as well as slippage and repair assays, shows a mildly pathogenic phenotype. Using cryogenic electron microscopy and surface plasmon resonance, we explored the mechanistic effect of this mutation on MutSα function. We discovered that V63E disrupts a previously unappreciated interface between the mismatch binding domains (MBDs) of MSH2 and MSH6 and leads to reduced DNA binding. Our research identifies this interface as a 'safety lock' that ensures high-affinity DNA binding to increase replication fidelity. Our mechanistic model explains the hypomorphic phenotype of the V63E patient mutation and other variants in the MBD interface.

An update on the pathological classification of breast cancer

Breast cancer (BC) is a heterogeneous disease, encompassing a diverse spectrum of tumours with varying morphological, biological, and clinical phenotypes. Although tumours may show phenotypic overlap, they often display different biological behaviour and response to therapy. Advances in high-throughput molecular techniques and bioinformatics have contributed to improved understanding of BC biology and refinement of molecular taxonomy with the identification of specific molecular subclasses. Although the traditional pathological morphological classification of BC is of paramount importance and provides diagnostic and prognostic information, current interest focusses on the use of a single gene and multigene assays to stratify BC into distinct groups to guide decisions on systemic therapy. This review considers approaches to the classification of BC, including their limitations, and with particular emphasis on the fundamental role of morphology in establishing an accurate diagnosis of primary invasive carcinoma of breast origin. This forms the basis for further morphological characterization and for all other approaches to BC classification that are used to provide prognostic and therapeutic predictive information.

RNA-Based Classification of Homologous Recombination Deficiency in Racially Diverse Patients with Breast Cancer

BACKGROUND

Aberrant expression of DNA repair pathways such as homologous recombination (HR) can lead to DNA repair imbalance, genomic instability, and altered chemotherapy response. DNA repair imbalance may predict prognosis, but variation in DNA repair in diverse cohorts of breast cancer patients is understudied.

METHODS

To identify RNA-based patterns of DNA repair expression, we performed unsupervised clustering on 51 DNA repair-related genes in the Cancer Genome Atlas Breast Cancer [TCGA BRCA (n = 1,094)] and Carolina Breast Cancer Study [CBCS (n = 1,461)]. Using published DNA-based HR deficiency (HRD) scores (high-HRD ≥ 42) from TCGA, we trained an RNA-based supervised classifier. Unsupervised and supervised HRD classifiers were evaluated in association with demographics, tumor characteristics, and clinical outcomes.

RESULTS

: Unsupervised clustering on DNA repair genes identified four clusters of breast tumors, with one group having high expression of HR genes. Approximately 39.7% of CBCS and 29.3% of TCGA breast tumors had this unsupervised high-HRD (U-HRD) profile. A supervised HRD classifier (S-HRD) trained on TCGA had 84% sensitivity and 73% specificity to detect HRD-high samples. Both U-HRD and S-HRD tumors in CBCS had higher frequency of TP53 mutant-like status (45% and 41% enrichment) and basal-like subtype (63% and 58% enrichment). S-HRD high was more common among black patients. Among chemotherapy-treated participants, recurrence was associated with S-HRD high (HR: 2.38, 95% confidence interval = 1.50-3.78).

CONCLUSIONS

HRD is associated with poor prognosis and enriched in the tumors of black women.

IMPACT

RNA-level indicators of HRD are predictive of breast cancer outcomes in diverse populations.

Homologous Recombination Deficiency Scar: Mutations and Beyond-Implications for Precision Oncology

Homologous recombination deficiency (HRD) is a prevalent in approximately 17% of tumors and is associated with enhanced sensitivity to anticancer therapies inducing double-strand DNA breaks. Accurate detection of HRD would therefore allow improved patient selection and outcome of conventional and targeted anticancer therapies. However, current clinical assessment of HRD mainly relies on determining germline BRCA1/2 mutational status and is insufficient for adequate patient stratification as mechanisms of HRD occurrence extend beyond functional BRCA1/2 loss. HRD, regardless of BRCA1/2 status, is associated with specific forms of genomic and mutational signatures termed HRD scar. Detection of this HRD scar might therefore be a more reliable biomarker for HRD. This review discusses and compares different methods of assessing HRD and HRD scar, their advances into the clinic, and their potential implications for precision oncology.

Loss of genes in chromosome arms 5q and 16q in breast cancer

INTRODUCTION

Losses of genetic material from chromosomes 5q and 16q commonly occur in sub-sets of breast cancer. Their significance from a pathophysiologic point of view is not well-defined.

METHODS

This study uses publicly available genomic data from extensive breast cancer datasets to define the landscape of losses in chromosomal arms 5q and 16q in the two sub-types of breast cancer they most commonly occur, basal-like and luminal A cancers, respectively.

RESULTS

It is shown that dozens of genes from these chromosomal arms are putatively hemi-deleted in few samples each. No individual gene from either 5q or 16q shows an incidence of deep deletion above 10% in the cohorts with basal-like and luminal A cancers or in the whole cohorts. A few tumor suppressor genes are deleted in a small number of samples, less than 5% in each cohort. Losses of 5q or 16q confer no survival advantage in either the basal-like or the luminal A cohorts from TCGA that harbor them.

CONCLUSION

Results suggest that there are no individual genes in chromosomes 5q and 16q whose loss can be implicated in a dominant pathophysiologic sequence of events in breast cancer or its sub-sets.

Somatic DNA Damage Response and Homologous Repair Gene Alterations and Its Association With Tumor Variant Burden in Breast Cancer Patients With Occupational Exposure to Pesticides

Homologous recombination is a crucial pathway that is specialized in repairing double-strand breaks; thus, alterations in genes of this pathway may lead to loss of genomic stability and cell growth suppression. Pesticide exposure potentially increases cancer risk through several mechanisms, such as the genotoxicity caused by chronic exposure, leading to gene alteration. To analyze this hypothesis, we investigated if breast cancer patients exposed to pesticides present a different mutational pattern in genes related to homologous recombination (BRCA1, BRCA2, PALB2, and RAD51D) and damage-response (TP53) concerning unexposed patients. We performed multiplex PCR-based assays and next-generation sequencing (NGS) of all coding regions and flanking splicing sites of BRCA1, BRCA2, PALB2, TP53, and RAD51D in 158 unpaired tumor samples from breast cancer patients on MiSeq (Illumina) platform. We found that exposed patients had tumors with more pathogenic and likely pathogenic variants than unexposed patients (p = 0.017). In general, tumors that harbored a pathogenic or likely pathogenic variant had a higher mutational burden (p < 0.001). We also observed that breast cancer patients exposed to pesticides had a higher mutational burden when diagnosed before 50 years old (p = 0.00978) and/or when carrying BRCA1 (p = 0.0138), BRCA2 (p = 0.0366), and/or PALB2 (p = 0.00058) variants, a result not found in the unexposed group. Our results show that pesticide exposure impacts the tumor mutational landscape and could be associated with the carcinogenesis process, therapy response, and disease progression. Further studies should increase the observation period in exposed patients to better evaluate the impact of these findings.

Kindlin-2 promotes Src-mediated tyrosine phosphorylation of androgen receptor and contributes to breast cancer progression

Androgen receptor (AR) signaling plays important roles in breast cancer progression. We show here that Kindlin-2, a focal adhesion protein, is critically involved in the promotion of AR signaling and breast cancer progression. Kindlin-2 physically associates with AR and Src through its two neighboring domains, namely F1 and F0 domains, resulting in formation of a Kindlin-2-AR-Src supramolecular complex and consequently facilitating Src-mediated AR Tyr-534 phosphorylation and signaling. Depletion of Kindlin-2 was sufficient to suppress Src-mediated AR Tyr-534 phosphorylation and signaling, resulting in diminished breast cancer cell proliferation and migration. Re-expression of wild-type Kindlin-2, but not AR-binding-defective or Src-binding-defective mutant forms of Kindlin-2, in Kindlin-2-deficient cells restored AR Tyr-534 phosphorylation, signaling, breast cancer cell proliferation and migration. Furthermore, re-introduction of phosphor-mimic mutant AR-Y534D, but not wild-type AR reversed Kindlin-2 deficiency-induced inhibition of AR signaling and breast cancer progression. Finally, using a genetic knockout strategy, we show that ablation of Kindlin-2 from mammary tumors in mouse significantly reduced AR Tyr-534 phosphorylation, breast tumor progression and metastasis in vivo. Our results suggest a critical role of Kindlin-2 in promoting breast cancer progression and shed light on the molecular mechanism through which it functions in this process.

Alterations in homologous recombination repair genes in prostate cancer brain metastases

Improved survival rates for prostate cancer through more effective therapies have also led to an increase in the diagnosis of metastases to infrequent locations such as the brain. Here we investigate the repertoire of somatic genetic alterations present in brain metastases from 51 patients with prostate cancer brain metastases (PCBM). We highlight the clonal evolution occurring in PCBM and demonstrate an increased mutational burden, concomitant with an enrichment of the homologous recombination deficiency mutational signature in PCBM compared to non-brain metastases. Focusing on known pathogenic alterations within homologous recombination repair genes, we find 10 patients (19.6%) fulfilling the inclusion criteria used in the PROfound clinical trial, which assessed the efficacy of PARP inhibitors (PARPi) in homologous recombination deficient prostate cancer. Eight (15.7%) patients show biallelic loss of one of the 15 genes included in the trial, while 5 patients (9.8%) harbor pathogenic alterations in BRCA1/2 specifically. Uncovering these molecular features of PCBM may have therapeutic implications, suggesting the need of clinical trial enrollment of PCBM patients when evaluating potential benefit from PARPi.

The Significance of Tumor Microenvironment Score for Breast Cancer Patients

Purpose

This study was designed to clarify the prognostic value of tumor microenvironment score and abnormal genomic alterations in TME for breast cancer patients.

Method

The TCGA-BRCA data were downloaded from TCGA and analyzed with R software. The results from analyses were further validated using the dataset from GSE96058, GSE124647, and GSE25066.

Results

After analyzing the TCGA data and verifying it with the GEO data, we developed a TMEscore model based on the TME infiltration pattern and validated it in 3273 breast cancer patients. The results suggested that our TMEscore model has high prognostic value. TME features with the TMEscore model can help to predict breast cancer patients' response to immunotherapy and provide new strategies for breast cancer treatment. Signature 24 was first found in breast cancer. In focal SCNAs, a total of 95 amplified genes and 169 deletion genes in the TMEscore high group were found to be significantly related to the prognosis of breast cancer patients, while 61 amplified genes and 174 deletion genes in the TMEscore low group were identified. LRRC48, CFAP69, and cg25726128 were first discovered and reported to be related to the survival of breast cancer patients. We identified specific mutation signatures that correlate with TMEscore and prognosis.

Conclusion

TMEscore model has high predictive value regarding prognosis and patients' response to immunotherapy. Signature 24 was first found in breast cancer. Specific mutation signatures that correlate with TMEscore and prognosis might be used for providing additional indicators for disease evaluation.

Genomic analysis of focal nodular hyperplasia with associated hepatocellular carcinoma unveils its malignant potential: a case report

Background

Focal nodular hyperplasia (FNH) is typically considered a benign tumor of the liver without malignant potential. The co-occurrence of FNH and hepatocellular carcinoma (HCC) has been reported in rare cases. In this study we sought to investigate the clonal relationship between these lesions in a patient with FNH-HCC co-occurrence.

Methods

A 74-year-old female patient underwent liver tumor resection. The resected nodule was subjected to histologic analyses using hematoxylin and eosin stain and immunohistochemistry. DNA extracted from microdissected FNH and HCC regions was subjected to whole exome sequencing. Clonality analysis were performed using PyClone.

Results

Histologic analysis reveals that the nodule consists of an FNH and two adjoining HCC components with distinct histopathological features. Immunophenotypic characterization and genomic analyses suggest that the FNH is clonally related to the HCC components, and is composed of multiple clones at diagnosis, that are likely to have progressed to HCC through clonal selection and/or the acquisition of additional genetic events.

Conclusion

To the best of our knowledge, our work is the first study showing a clonal relationship between FNH and HCC. We show that FNH may possess the capability to undergo malignant transformation and to progress to HCC in very rare cases.

Defining genomic events involved in the evolutionary trajectories of myeloma and its precursor conditions

All patients with a diagnosis of multiple myeloma (MM) have a preceding, asymptomatic expansion of clonal plasma cells, clinically recognized as monoclonal gammopathy of undetermined significance or smoldering multiple myeloma (SMM). While most patients with monoclonal gammopathy of undetermined significance have a very small rate of progression, SMM is a widely heterogeneous condition where a fraction of patients will progress to symptomatic MM rather quickly, while others will experience an indolent clinical course. The differentiation between progressive and stable precursor condition thus represents one of the most important unmet clinical needs in the MM community. The ability to identify patients at high-risk of progression before major clonal expansion and onset of end-organ damage would enable strategies for early prevention and perhaps more effective intervention. All proposed criteria to predict the progression of myeloma precursor conditions are built around indirect markers of disease burden and, therefore, are generally able to accurately identify only a small fraction of patients in whom progression to MM is already occurring. Leveraging whole genome and exome sequencing, it has been shown that patients with stable myeloma precursor conditions are characterized by either absence or lower prevalence of distinct genomic events that are detectable in progressive precursor condition years before the progression. In this review, we discuss evolving genomic concepts and tools; and their ability to differentiate myeloma precursor conditions into two distinct entities: one benign (monoclonal gammopathy of benign significance) and another malignant (asymptomatic multiple myeloma).

Identifying Transcripts with Tandem Duplications from RNA-Sequencing Data to Predict BRCA1-Type Primary Breast Cancer

Simple Summary

Homologous recombination repair deficiency (HRD) is a biomarker for the response to PARP inhibitor anti-cancer treatment. Therefore, methods that detect the HRD phenotype in cancers in a (cost-)effective manner are pivotal. In this respect, the HRDetect and CHORD algorithms were developed to classify (the type of) HRD cancers from whole genome sequencing data. In addition, functional assays have also been established, but these require fresh cancer tissue. Here we present a novel method to specifically classify BRCA1-type HRD from RNA-sequencing data with high sensitivity. BRCA1-type cancers typically display small (<10 kb) tandem duplications, in contrast to BRCA2-type cancers. By detecting these small TDs among transcripts, we increase the toolbox for detecting HRD with a method that does not require whole genome sequencing of both tumor and normal tissue.

Abstract

Patients with cancers that are deficient for homologous recombination repair (HRD) may benefit from PARP inhibitor treatment. Therefore, methods that identify such cancers are crucial. Using whole genome sequencing data, specific genomic scars derived from somatic mutations and genomic rearrangements can identify HRD tumors, with only BRCA1-like HRD cancers profoundly displaying small (<10 kb) tandem duplications (TDs). In this manuscript we describe a method of detecting BRCA1-type HRD in breast cancer (BC) solely from RNA sequencing data by identifying TDs surfacing in transcribed genes. We find that the number of identified TDs (TD-score) is significantly higher in BRCA1-type vs. BRCA2-type BCs, or vs. HR-proficient BCs (p = 2.4 × 10⁻⁶ and p = 2.7 × 10⁻¹², respectively). A TD-score ≥2 shows an 88.2% sensitivity (30 out of 34) to detect a BRCA1-type BC, with a specificity of 64.7% (143 out of 221). Pathway enrichment analyses showed genes implicated in cancer to be affected by TDs of which PTEN was found significantly more frequently affected by a TD in BRCA1-type BC. In conclusion, we here describe a novel method to identify TDs in transcripts and classify BRCA1-type BCs with high sensitivity.

Pathogenesis of Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) encompasses a heterogeneous group of fundamentally different diseases with different histologic, genomic, and immunologic profiles, which are aggregated under this term because of their lack of estrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2 expression. Massively parallel sequencing and other omics technologies have demonstrated the level of heterogeneity in TNBCs and shed light into the pathogenesis of this therapeutically challenging entity in breast cancer. In this review, we discuss the histologic and molecular classifications of TNBC, the genomic alterations these different tumor types harbor, and the potential impact of these alterations on the pathogenesis of these tumors. We also explore the role of the tumor microenvironment in the biology of TNBCs and its potential impact on therapeutic response. Dissecting the biology and understanding the therapeutic dependencies of each TNBC subtype will be essential to delivering on the promise of precision medicine for patients with triple-negative disease.

Multi-omics analysis identifies therapeutic vulnerabilities in triple-negative breast cancer subtypes

Triple-negative breast cancer (TNBC) is a collection of biologically diverse cancers characterized by distinct transcriptional patterns, biology, and immune composition. TNBCs subtypes include two basal-like (BL1, BL2), a mesenchymal (M) and a luminal androgen receptor (LAR) subtype. Through a comprehensive analysis of mutation, copy number, transcriptomic, epigenetic, proteomic, and phospho-proteomic patterns we describe the genomic landscape of TNBC subtypes. Mesenchymal subtype tumors display high mutation loads, genomic instability, absence of immune cells, low PD-L1 expression, decreased global DNA methylation, and transcriptional repression of antigen presentation genes. We demonstrate that major histocompatibility complex I (MHC-I) is transcriptionally suppressed by H3K27me3 modifications by the polycomb repressor complex 2 (PRC2). Pharmacological inhibition of PRC2 subunits EZH2 or EED restores MHC-I expression and enhances chemotherapy efficacy in murine tumor models, providing a rationale for using PRC2 inhibitors in PD-L1 negative mesenchymal tumors. Subtype-specific differences in immune cell composition and differential genetic/pharmacological vulnerabilities suggest additional treatment strategies for TNBC.

Triple negative breast cancer can be divided into additional subtypes. Here, using omics analyses, the authors show that in the mesenchymal subtype expression of MHC-1 is repressed and that this can be restored by using drugs that target subunits of the epigenetic modifier PRC2.

Genetic heterogeneity during breast cancer progression in young patients

Background

Because a number of years may be required for normal cells to develop into carcinoma, genes involved in tumorigenesis and progression might differ among breast cancers in young women and those in older women. The present study sought to analyze subclonality during breast cancer evolution as well as diversity within each individual in our young patients’ cohort.

Methods

A total of 13 women aged <35 years at diagnosis with early breast cancer were recruited. Serial sections of breast samples consisting of synchronous invasive carcinoma, adjacent ductal carcinoma in situ (DCIS), normal breast tissue, and metastatic lymph nodes were collected and prepared for immunohistochemical analysis of estrogen receptor, progesterone receptor, HER2, and Ki67, and for extraction of genomic DNA. Germline and somatic gene alterations of genomic DNA were examined by targeted sequencing.

Results

Genomic DNA from 13 blood samples and 36 breast tissues consisting of 14 invasive carcinomas, nine adjacent DCIS, 11 normal breast tissues, and two metastatic lymph nodes were successfully sequenced. Germline gene alterations including pathogenic variants and gene alterations that were not yet evaluated for their clinical significance were detected in all patients but one. Somatic gene alterations were identified in eight invasive carcinomas, five DCIS, and one metastatic lymph node. Different somatic gene alterations between invasive carcinoma and DCIS were detected in two patients. Somatic gene mutations were present in non-neoplastic tissues in three patients. No two patients had the same gene alterations.

Conclusion

Our results reveal diversity within each individual during breast cancer progression.

•

Subclonality during breast cancer progression in young patients were analyzed.

•

No two patients had the same germline and/or somatic gene alterations.

•

Somatic gene mutations were present in normal breast tissue.

The Evolution of Clinically Aggressive Triple-Negative Breast Cancer Shows a Large Mutational Diversity and Early Metastasis to Lymph Nodes

In triple-negative breast cancer (TNBC), only 30% of patients treated with neoadjuvant chemotherapy achieve a pathological complete response after treatment and more than 90% die due to metastasis formation. The diverse clinical responses and metastatic developments are attributed to extensive intrapatient genetic heterogeneity and tumor evolution acting on this neoplasm. In this work, we aimed to evaluate genomic alterations and tumor evolution in TNBC patients with aggressive disease. We sequenced the whole exome of 16 lesions from four patients who did not respond to therapy, and took several follow-up samples, including samples from tumors before and after treatment, as well as from the lymph nodes and skin metastases. We found substantial intrapatient genetic heterogeneity, with a variable tumor mutational composition. Early truncal events were MCL1 amplifications. Metastatic lesions had deletions in RB1 and PTEN, along with TERT, AKT2, and CCNE1 amplifications. Mutational signatures 06 and 12 were mainly detected in skin metastases and lymph nodes. According to phylogenetic analysis, the lymph node metastases occurred at an early stage of TNBC development. Finally, each patient had three to eight candidate driving mutations for targeted treatments. This study delves into the genomic complexity and the phylogenetic and evolutionary development of aggressive TNBC, supporting early metastatic development, and identifies specific genetic alterations associated with a response to targeted therapies.

Integrated analysis reveals prognostic value of HLA-I LOH in triple-negative breast cancer

Background

Triple-negative breast cancers (TNBCs), especially those non-immune-inflamed tumors, have a poor prognosis and limited therapies. Human leukocyte antigen (HLA)-I not only contributes to antitumor immune response and the phenotype of the tumor microenvironment, but also is a negative predictor of outcomes after immunotherapy. However, the importance of HLA functional status in TNBCs remains poorly understood.

Methods

Using the largest original multiomics datasets on TNBCs, we systematically characterized the HLA-Ⅰ status of TNBCs from the perspective of HLA-Ⅰ homogeneity and loss of heterozygosity (LOH). The prognostic significance of HLA-I status was measured. To explain the potential mechanism of prognostic value in HLA-Ⅰ status, the mutational signature, copy number alteration, neoantigen and intratumoral heterogeneity were measured. Furthermore, the correlation between HLA-Ⅰ functional status and the tumor immune microenvironment was analyzed.

Results

LOH and homogeneity in HLA-I accounted for 18% and 21% of TNBCs, respectively. HLA-I LOH instead of HLA-I homogeneity was an independent prognostic biomarker in TNBCs. In particular, for patients with non-immune-inflamed tumors, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH. Furthermore, integrated genomic and transcriptomic analysis showed that HLA-I LOH was accompanied by upregulated scores of mutational signature 3 and homologous recombination deficiency scores, which implied the failure of DNA double-strand break repair. Moreover, HLA-I LOH had higher mutation and neoantigen loads and more subclones than HLA-I non-LOH. These results indicated that although HLA-I LOH tumors with failure of DNA double-strand break repair were prone to produce neoantigens, their limited capacity for antigen presentation finally contributed to poor immune selection pressure.

Conclusion

Our study illustrates the genomic landscape of HLA-I functional status and stresses the prognostic significance of HLA-I LOH in TNBCs. For “cold” tumors in TNBCs, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH.

Genomic profiling using the UltraSEEK panel identifies discordancy between paired primary and breast cancer brain metastases and an association with brain metastasis-free survival

PURPOSE

Brain metastases (BM) are an increasing clinical problem. This study aimed to assess paired primary breast cancers (BC) and BM for aberrations within TP53, PIK3CA, ESR1, ERBB2 and AKT utilising the MassARRAY® UltraSEEK® technology (Agena Bioscience, San Diego, USA).

METHODS

DNA isolated from 32 paired primary BCs and BMs was screened using the custom UltraSEEK® Breast Cancer Panel. Data acquisition and analysis was performed by the Agena Bioscience Typer software v4.0.26.74.

RESULTS

Mutations were identified in 91% primary BCs and 88% BM cases. TP53, AKT1, ESR1, PIK3CA and ERBB2 genes were mutated in 68.8%, 37.5%, 31.3%, 28.1% and 3.1% respectively of primary BCs and in 59.4%, 37.5%, 28.1%, 28.1% and 3.1% respectively of BMs. Differences in the mutations within the 5 genes between BC and paired BM were identified in 62.5% of paired cases. In primary BCs, ER-positive/HER2-negative cases harboured the most mutations (70%), followed by ER-positive/HER2-positive (15%) and triple-negatives (13.4%), whereas in BMs, the highest number of mutations was observed in triple-negative (52.5%), followed by ER-positive/HER2-negative (35.6%) and ER-negative/HER2-positive (12%). There was a significant association between the number of mutations in the primary BC and breast-to-brain metastasis-free survival (p = 0.0001) but not with overall survival (p = 0.056).

CONCLUSION

These data demonstrate the discordancy between primary BC and BM, as well as the presence of clinically important, actionable mutations in BCBM. The UltraSEEK® Breast Cancer Panel provides a tool for BCBM that can be utilised to direct more tailored treatment decisions and for clinical studies investigating targeted agents.

Durable Complete Response With Immune Checkpoint Inhibitor in Breast Cancer With High Tumor Mutational Burden and APOBEC Signature

Increasing data support the importance of preexisting host immune response and neoantigen burden for determining response to immune checkpoint inhibitors (ICIs). In lung cancer and melanoma, tumor mutational burden (TMB) has emerged as an independent biomarker for ICI response. However, the significance of TMB in breast cancer, particularly in the context of PD-L1 negativity, remains unclear. This report describes a patient with HER2-negative breast cancer with high TMB and an apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) trinucleotide signature; her disease was refractory to multiple lines of treatments but achieved durable complete response using ICIs and capecitabine. Additional analysis of the tumor revealed a low amount of stromal tumor-infiltrating lymphocytes (sTILs) and PD-L1 negativity, reflecting a poor preexisting host immune response. In collaboration with Foundation Medicine, comprehensive genomic profiling from 14,867 patients with breast cancer with the FoundationOne test was evaluated. Using the cutoff of ≥10 mutations/megabase (mut/Mb) for high TMB, PD-L1 positivity and TMB-high populations were not significantly overlapping (odds ratio, 1.02; P=.87). Up to 79% of TMB-high tumors with >20 mut/Mb were PD-L1-negative. Our study highlights that despite having low TILs and PD-L1 negativity, some patients may still experience response to ICIs.

APOBEC3B expression in breast cancer cell lines and tumors depends on the estrogen receptor status

Increased exposure to estrogen is associated with an elevated risk of breast cancer. Considering estrogen as a possible mutagen, we hypothesized that exposure to estrogen alone or in combination with the DNA-damaging chemotherapy drug, cisplatin, could induce expression of genes encoding enzymes involved in APOBEC-mediated mutagenesis. To test this hypothesis, we measured the expression of APOBEC3A (A3A) and APOBEC3B (A3B) genes in two breast cancer cell lines treated with estradiol, cisplatin or their combination. These cell lines, T-47D (ER+) and MDA-MB-231 (ER-), differed by the status of the estrogen receptor (ER). Expression of A3A was not detectable in any conditions tested, while A3B expression was induced by treatment with cisplatin and estradiol in ER+ cells but was not affected by estradiol in ER- cells. In The Cancer Genome Atlas, expression of A3B was significantly associated with genotypes of a regulatory germline variant rs17000526 upstream of the APOBEC3 cluster in 116 ER- breast tumors (P = 0.006) but not in 387 ER+ tumors (P = 0.48). In conclusion, we show that in breast cancer cell lines, A3B expression was induced by estradiol in ER+ cells and by cisplatin regardless of ER status. In ER+ breast tumors, the effect of estrogen may be masking the association of rs17000526 with A3B expression, which was apparent in ER- tumors. Our results provide new insights into the differential etiology of ER+ and ER- breast cancer and the possible role of A3B in this process through a mitogenic rather than the mutagenic activity of estrogen.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Impact of clinical targeted sequencing on endocrine responsiveness in estrogen receptor-positive, HER2-negative metastatic breast cancer

Clinical targeted sequencing allows for the selection of patients expected to have a better treatment response, and reveals mechanisms of resistance to molecular targeted therapies based on actionable gene mutations. We underwent comprehensive genomic testing with either our original in-house CLHURC system or with OncoPrime. Samples from 24 patients with estrogen receptor-positive, human epidermal growth factor receptor 2-negative metastatic breast cancer underwent targeted sequencing between 2016 and 2018. Germline and somatic gene alterations and patients' prognosis were retrospectively analyzed according to the response to endocrine therapy. All of the patients had one or more germline and/or somatic gene alterations. Four patients with primary or secondary endocrine-resistant breast cancer harbored germline pathogenic variants of BRCA1, BRCA2, or PTEN. Among somatic gene alterations, TP53, PIK3CA, AKT1, ESR1, and MYC were the most frequently mutated genes. TP53 gene mutation was more frequently observed in patients with primary endocrine resistance compared to those with secondary endocrine resistance or endocrine-responsive breast cancer. Recurrent breast cancer patients carrying TP53-mutant tumors had significantly worse overall survival compared to those with TP53-wild type tumors. Our 160-gene cancer panel will be useful to identify clinically actionable gene alterations in breast cancer in clinical practice.

EFA6B regulates a stop signal for collective invasion in breast cancer

Cancer is initiated by somatic mutations in oncogenes or tumor suppressor genes. However, additional alterations provide selective advantages to the tumor cells to resist treatment and develop metastases. Their identification is of paramount importance. Reduced expression of EFA6B (Exchange Factor for ARF6, B) is associated with breast cancer of poor prognosis. Here, we report that loss of EFA6B triggers a transcriptional reprogramming of the cell-to-ECM interaction machinery and unleashes CDC42-dependent collective invasion in collagen. In xenograft experiments, MCF10 DCIS.com cells, a DCIS-to-IDC transition model, invades faster when knocked-out for EFA6B. In addition, invasive and metastatic tumors isolated from patients have lower expression of EFA6B and display gene ontology signatures identical to those of EFA6B knock-out cells. Thus, we reveal an EFA6B-regulated molecular mechanism that controls the invasive potential of mammary cells; this finding opens up avenues for the treatment of invasive breast cancer.

Evolutionary Trajectories and Genomic Divergence in Localized Breast Cancers after Ipsilateral Breast Tumor Recurrence

Simple Summary

Ipsilateral breast tumor relapse (IBTR) occurs in 5–10% of localized breast cancers (BCs) within 10 years of incidence, despite proper treatment of the primary lesion. However, the clinical consequences of evolutionary trajectories of BC cells and their impact on IBTR remain poorly understood. Here, we conducted a longitudinal genomic analysis of 10 matched localized BC patients with IBTR. Overall, we identified the differences in homologous recombination deficiency, chromosomal instability, and somatic mutation drivers between primary and relapsed lesions. Our analyses highlighted three clonal architectures that shape by distinct mutagenic processes and subclonal diversification during relapse progression. Finally, this study provided a framework, which integrated actionable biomarkers with clonal architectures, towards improvement of future treatment decisions.

Abstract

The evolutionary trajectories that drive clinical and therapeutic consequences in localized breast cancers (BCs) with ipsilateral breast tumor relapse (IBTR) remain largely unknown. Analyses of longitudinal paired whole-exome sequencing data from 10 localized BC patients with IBTR reveal that, compared to primary breast tumors, homologous recombination (HR) deficiency, inactivation of the HR pathway, chromosomal instability, and somatic driver mutations are more frequent. Furthermore, three major models of evolution in IBTR are summarized, through which relative contributions of mutational signatures shift, and the subclonal diversity expansions are shown. Optimal treatment regimens are suggested by the clinically relevant molecular features, such as HR deficiency (20%) or specific alterations (30%) with sensitivity to available FDA-approved drugs. Finally, a rationale for the development of the therapeutic management framework is provided. This study sheds light on the complicated evolution patterns in IBTR and has significant clinical implications for future improvement of treatment decisions.

Genomic landscape of extraordinary responses in metastatic breast cancer

Extreme responders to anticancer therapy are rare among advanced breast cancer patients. Researchers, however, have yet to investigate treatment responses therein on the whole exome level. We performed whole exome analysis to characterize the genomic landscape of extreme responders among metastatic breast cancer patients. Clinical samples were obtained from breast cancer patients who showed exceptional responses to anti-HER2 therapy or hormonal therapy and from those who did not. Matched breast tumor tissue (somatic DNA) and blood samples (germline DNA) were collected from a total of 30 responders and 15 non-responders. Whole exome sequencing using Illumina HiSeq2500 was performed for all 45 patients (90 samples). Somatic single nucleotide variants (SNVs), indels, and copy number variants (CNVs) were identified for the genomes of each patient. Group-specific somatic variants and mutational burden were statistically analyzed. Sequencing of cancer exomes for all patients revealed 1839 somatic SNVs (1661 missense, 120 nonsense, 43 splice-site, 15 start/stop-lost) and 368 insertions/deletions (273 frameshift, 95 in-frame), with a median of 0.7 mutations per megabase (range, 0.08 to 4.2 mutations per megabase). Responders harbored a significantly lower nonsynonymous mutational burden (median, 26 vs. 59, P = 0.02) and fewer CNVs (median 13.6 vs. 97.7, P = 0.05) than non-responders. Multivariate analyses of factors influencing progression-free survival showed that a high mutational burden and visceral metastases were significantly related with disease progression. Extreme responders to treatment for metastatic breast cancer are characterized by fewer nonsynonymous mutations and CNVs.

Molecular analyses of triple-negative breast cancer in the young and elderly

Background

Breast cancer in young adults has been implicated with a worse outcome. Analyses of genomic traits associated with age have been heterogenous, likely because of an incomplete accounting for underlying molecular subtypes. We aimed to resolve whether triple-negative breast cancer (TNBC) in younger versus older patients represent similar or different molecular diseases in the context of genetic and transcriptional subtypes and immune cell infiltration.

Patients and methods

In total, 237 patients from a reported population-based south Swedish TNBC cohort profiled by RNA sequencing and whole-genome sequencing (WGS) were included. Patients were binned in 10-year intervals. Complimentary PD-L1 and CD20 immunohistochemistry and estimation of tumor-infiltrating lymphocytes (TILs) were performed. Cases were analyzed for differences in patient outcome, genomic, transcriptional, and immune landscape features versus age at diagnosis. Additionally, 560 public WGS breast cancer profiles were used for validation.

Results

Median age at diagnosis was 62 years (range 26–91). Age was not associated with invasive disease-free survival or overall survival after adjuvant chemotherapy. Among the BRCA1-deficient cases (82/237), 90% were diagnosed before the age of 70 and were predominantly of the basal-like subtype. In the full TNBC cohort, reported associations of patient age with changes in Ki67 expression, PIK3CA mutations, and a luminal androgen receptor subtype were confirmed. Within DNA repair deficiency or gene expression defined molecular subgroups, age-related alterations in, e.g., overall gene expression, immune cell marker gene expression, genetic mutational and rearrangement signatures, amount of copy number alterations, and tumor mutational burden did, however, not appear distinct. Similar non-significant associations for genetic alterations with age were obtained for other breast cancer subgroups in public WGS data. Consistent with age-related immunosenescence, TIL counts decreased linearly with patient age across different genetic TNBC subtypes.

Conclusions

Age-related alterations in TNBC, as well as breast cancer in general, need to be viewed in the context of underlying genomic phenotypes. Based on this notion, age at diagnosis alone does not appear to provide an additional layer of biological complexity above that of proposed genetic and transcriptional phenotypes of TNBC. Consequently, treatment decisions should be less influenced by age and more driven by tumor biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01392-0.

Pathway-Based Drug-Repurposing Schemes in Cancer: The Role of Translational Bioinformatics

Cancer is a set of complex pathologies that has been recognized as a major public health problem worldwide for decades. A myriad of therapeutic strategies is indeed available. However, the wide variability in tumor physiology, response to therapy, added to multi-drug resistance poses enormous challenges in clinical oncology. The last years have witnessed a fast-paced development of novel experimental and translational approaches to therapeutics, that supplemented with computational and theoretical advances are opening promising avenues to cope with cancer defiances. At the core of these advances, there is a strong conceptual shift from gene-centric emphasis on driver mutations in specific oncogenes and tumor suppressors-let us call that the silver bullet approach to cancer therapeutics-to a systemic, semi-mechanistic approach based on pathway perturbations and global molecular and physiological regulatory patterns-we will call this the shrapnel approach. The silver bullet approach is still the best one to follow when clonal mutations in driver genes are present in the patient, and when there are targeted therapies to tackle those. Unfortunately, due to the heterogeneous nature of tumors this is not the common case. The wide molecular variability in the mutational level often is reduced to a much smaller set of pathway-based dysfunctions as evidenced by the well-known hallmarks of cancer. In such cases "shrapnel gunshots" may become more effective than "silver bullets". Here, we will briefly present both approaches and will abound on the discussion on the state of the art of pathway-based therapeutic designs from a translational bioinformatics and computational oncology perspective. Further development of these approaches depends on building collaborative, multidisciplinary teams to resort to the expertise of clinical oncologists, oncological surgeons, and molecular oncologists, but also of cancer cell biologists and pharmacologists, as well as bioinformaticians, computational biologists and data scientists. These teams will be capable of engaging on a cycle of analyzing high-throughput experiments, mining databases, researching on clinical data, validating the findings, and improving clinical outcomes for the benefits of the oncological patients.

Reconstructing tumor history in breast cancer: signatures of mutational processes and response to neoadjuvant chemotherapy⋆

BACKGROUND

Different endogenous and exogenous mutational processes act over the evolutionary history of a malignant tumor, driven by abnormal DNA editing, mutagens or age-related DNA alterations, among others, to generate the specific mutational landscape of each individual tumor. The signatures of these mutational processes can be identified in large genomic datasets. We investigated the hypothesis that genomic patterns of mutational signatures are associated with the clinical behavior of breast cancer, in particular chemotherapy response and survival, with a particular focus on therapy-resistant disease.

PATIENTS AND METHODS

Whole exome sequencing was carried out in 405 pretherapeutic samples from the prospective neoadjuvant multicenter GeparSepto study. We analyzed 11 mutational signatures including biological processes such as APOBEC-mutagenesis, homologous recombination deficiency (HRD), mismatch repair deficiency and also age-related or tobacco-induced alterations.

RESULTS

Different subgroups of breast carcinomas were defined mainly by differences in HRD-related and APOBEC-related mutational signatures and significant differences between hormone-receptor (HR)-negative and HR-positive tumors as well as correlations with age, Ki-67 and immunological parameters were observed. We could identify mutational processes that were linked to increased pathological complete response rates to neoadjuvant chemotherapy with high significance. In univariate analyses for HR-positive tumors signatures, S3 (HRD, P < 0.001) and S13 (APOBEC, P = 0.001) as well as exonic mutation rate (P = 0.002) were significantly correlated with increased pathological complete response rates. The signatures S3 (HRD, P = 0.006) and S4 (tobacco, P = 0.011) were prognostic for reduced disease-free survival of patients with chemotherapy-resistant tumors.

CONCLUSION

The results of this investigation suggest that the clinical behavior of a tumor, in particular, response to neoadjuvant chemotherapy and disease-free survival of therapy-resistant tumors, could be predicted by the composition of mutational signatures as an indicator of the individual genomic history of a tumor. After additional validations, mutational signatures might be used to identify tumors with an increased response rate to neoadjuvant chemotherapy and to define therapy-resistant subgroups for future therapeutic interventions.

The interesting relationship between APOBEC3 deoxycytidine deaminases and cancer: a long road ahead

Cancer is considered a group of diseases characterized by uncontrolled growth and spread of abnormal cells and is propelled by somatic mutations. Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) family of enzymes are endogenous sources of somatic mutations found in multiple human cancers. While these enzymes normally act as an intrinsic immune defence against viruses, they can also catalyse 'off-target' cytidine deamination in genomic single-stranded DNA intermediates. The deamination of cytosine forms uracil, which is promutagenic in DNA. Key factors to trigger the APOBEC 'off-target' activity are overexpression in a non-normal cell type, nuclear localization and replication stress. The resulting uracil-induced mutations contribute to genomic variation, which may result in neutral, beneficial or harmful consequences for the cancer. This review summarizes the functional and biochemical basis of the APOBEC3 enzyme activity and highlights their relationship with the most well-studied cancers in this particular context such as breast, lung, bladder, and human papillomavirus-associated cancers. We focus on APOBEC3A, APOBEC3B and APOBEC3H haplotype I because they are the leading candidates as sources of somatic mutations in these and other cancers. Also, we discuss the prognostic value of the APOBEC3 expression in drug resistance and response to therapies.

RAD51-Mediated DNA Homologous Recombination Is Independent of PTEN Mutational Status

Simple Summary

PTEN is an important tumor suppressor that is frequently mutated in malignancy. PTEN mutational loss has been associated with reduced RAD51 expression and homologous recombination deficiency (HRD), however; recent studies have failed to recapitulate these findings. Here, we show that RAD51 expression, foci formation and homologous recombination repair activity are unaltered in normal and tumorigenic PTEN-deficient cells and patient samples. Furthermore, we show that PTEN-deficient tumor cell lines do not synergize with the clinical PARP inhibitor olaparib, underscoring a need to discontinue its use in treating patients with PTEN-deficient tumors that do not otherwise exhibit HRD.

Abstract

PTEN mutation occurs in a variety of aggressive cancers and is associated with poor patient outcomes. Recent studies have linked mutational loss of PTEN to reduced RAD51 expression and function, a key factor involved in the homologous recombination (HR) pathway. However, these studies remain controversial, as they fail to establish a definitive causal link to RAD51 expression that is PTEN-dependent, while other studies have not been able to recapitulate the relationship between the PTEN expression and the RAD51/HR function. Resolution of this apparent conundrum is essential due to the clinically-significant implication that PTEN-deficient tumors may be sensitive to poly (ADP-ribose) polymerase (PARP) inhibitors (PARPi) commonly used in the clinical management of BRCA-mutated and other HR-deficient (HRD) tumors. Methods: Primary Pten-deficient (and corresponding wild-type) mouse embryonic fibroblasts (MEFs) and astrocytes and PTEN-null human tumor cell lines and primary cells were assessed for RAD51 expression (via the Western blot analysis) and DNA damage repair analyses (via alkali comet and γH2AX foci assays). RAD51 foci analysis was used to measure HR-dependent DNA repair. Xrcc2-deficient MEFs served as an HR-deficient control, while the stable knockdown of RAD51 (shRAD51) served to control for the relative RAD51/HR-mediated repair and the phospho-53BP1 foci analysis served to confirm and measure non-homologous end joining (NHEJ) activity in PTEN-deficient and shRAD51-expressing (HRD) lines. Cell proliferation studies were used to measure any potential added sensitivity of PTEN-null cells to the clinically-relevant PARPi, olaparib. RAD51 levels and DNA damage response signaling were assessed in PTEN-mutant brain tumor initiating cells (BTICs) derived from primary and recurrent glioblastoma multiforme (GBM) patients, while expression of RAD51 and its paralogs were examined as a function of the PTEN status in the RNA expression datasets isolated from primary GBM tumor specimens and BTICs. Results: Pten knockout primary murine cells display unaltered RAD51 expression, endogenous and DNA strand break-induced RAD51 foci and robust DNA repair activity. Defective HR was only observed in the cells lacking Xrcc2. Likewise, human glioblastoma multiforme (GBM) cell lines with known PTEN deficiency (U87, PTEN-mutated; U251 and U373, PTEN-null) show apparent expression of RAD51 and display efficient DNA repair activity. Only GBM cells stably expressing shRNAs against RAD51 (shRAD51) display dysfunctional DNA repair activity and reduced proliferative capacity, which is exacerbated by PARPi treatment. Furthermore, GBM patient-derived BTICs displayed robust RAD51 expression and intact DNA damage response signaling in spite of PTEN-inactivating mutations. RNA expression analysis of primary GBM tissue specimens and BTICs demonstrate stable levels of RAD51 and its paralogs (RAD51B, RAD51C, RAD51D, XRCC2, XRCC3, and DMC1), regardless of the PTEN mutational status. Conclusions: Our findings demonstrate definitively that PTEN loss does not alter the RAD51 expression, its paralogs, or the HR activity. Furthermore, deficiency in PTEN alone is not sufficient to impart enhanced sensitivity to PARPi associated with HRD. This study is the first to unequivocally demonstrate that PTEN deficiency is not linked to the RAD51 expression or the HR activity amongst primary neural and non-neural Pten-null cells, PTEN-deficient tumor cell lines, and primary PTEN-mutant GBM patient-derived tissue specimens and BTICs.

MutSignatures: an R package for extraction and analysis of cancer mutational signatures

Cancer cells accumulate somatic mutations as result of DNA damage, inaccurate repair and other mechanisms. Different genetic instability processes result in characteristic non-random patterns of DNA mutations, also known as mutational signatures. We developed mutSignatures, an integrated R-based computational framework aimed at deciphering DNA mutational signatures. Our software provides advanced functions for importing DNA variants, computing mutation types, and extracting mutational signatures via non-negative matrix factorization. Specifically, mutSignatures accepts multiple types of input data, is compatible with non-human genomes, and supports the analysis of non-standard mutation types, such as tetra-nucleotide mutation types. We applied mutSignatures to analyze somatic mutations found in smoking-related cancer datasets. We characterized mutational signatures that were consistent with those reported before in independent investigations. Our work demonstrates that selected mutational signatures correlated with specific clinical and molecular features across different cancer types, and revealed complementarity of specific mutational patterns that has not previously been identified. In conclusion, we propose mutSignatures as a powerful open-source tool for detecting the molecular determinants of cancer and gathering insights into cancer biology and treatment.

Trajectory and uniqueness of mutational signatures in yeast mutators

The acquisition of mutations plays critical roles in adaptation, evolution, senescence, and tumorigenesis. Massive genome sequencing has allowed extraction of specific features of many mutational landscapes but it remains difficult to retrospectively determine the mechanistic origin(s), selective forces, and trajectories of transient or persistent mutations and genome rearrangements. Here, we conducted a prospective reciprocal approach to inactivate 13 single or multiple evolutionary conserved genes involved in distinct genome maintenance processes and characterize de novo mutations in 274 diploid Saccharomyces cerevisiae mutation accumulation lines. This approach revealed the diversity, complexity, and ultimate uniqueness of mutational landscapes, differently composed of base substitutions, small insertions/deletions (InDels), structural variants, and/or ploidy variations. Several landscapes parallel the repertoire of mutational signatures in human cancers while others are either novel or composites of subsignatures resulting from distinct DNA damage lesions. Notably, the increase of base substitutions in the homologous recombination-deficient Rad51 mutant, specifically dependent on the Polζ translesion polymerase, yields COSMIC signature 3 observed in BRCA1/BRCA2-mutant breast cancer tumors. Furthermore, "mutome" analyses in highly polymorphic diploids and single-cell bottleneck lineages revealed a diverse spectrum of loss-of-heterozygosity (LOH) signatures characterized by interstitial and terminal chromosomal events resulting from interhomolog mitotic cross-overs. Following the appearance of heterozygous mutations, the strong stimulation of LOHs in the rad27/FEN1 and tsa1/PRDX1 backgrounds leads to fixation of homozygous mutations or their loss along the lineage. Overall, these mutomes and their trajectories provide a mechanistic framework to understand the origin and dynamics of genome variations that accumulate during clonal evolution.

Mutational signature SBS8 predominantly arises due to late replication errors in cancer

Although a majority of somatic mutations in cancer are passengers, their mutational signatures provide mechanistic insights into mutagenesis and DNA repair processes. Mutational signature SBS8 is common in most cancers, but its etiology is debated. Incorporating genomic, epigenomic, and cellular process features for multiple cell-types we develop genome-wide composite epigenomic context-maps relevant for mutagenesis and DNA repair. Analyzing somatic mutation data from multiple cancer types in their epigenomic contexts, we show that SBS8 preferentially occurs in gene-poor, lamina-proximal, late replicating heterochromatin domains. While SBS8 is uncommon among mutations in non-malignant tissues, in tumor genomes its proportions increase with replication timing and speed, and checkpoint defects further promote this signature - suggesting that SBS8 probably arises due to uncorrected late replication errors during cancer progression. Our observations offer a potential reconciliation among different perspectives in the debate about the etiology of SBS8 and its relationship with other mutational signatures.

Comprehensive molecular comparison of BRCA1 hypermethylated and BRCA1 mutated triple negative breast cancers

Homologous recombination deficiency (HRD) is a defining characteristic in BRCA-deficient breast tumors caused by genetic or epigenetic alterations in key pathway genes. We investigated the frequency of BRCA1 promoter hypermethylation in 237 triple-negative breast cancers (TNBCs) from a population-based study using reported whole genome and RNA sequencing data, complemented with analyses of genetic, epigenetic, transcriptomic and immune infiltration phenotypes. We demonstrate that BRCA1 promoter hypermethylation is twice as frequent as BRCA1 pathogenic variants in early-stage TNBC and that hypermethylated and mutated cases have similarly improved prognosis after adjuvant chemotherapy. BRCA1 hypermethylation confers an HRD, immune cell type, genome-wide DNA methylation, and transcriptional phenotype similar to TNBC tumors with BRCA1-inactivating variants, and it can be observed in matched peripheral blood of patients with tumor hypermethylation. Hypermethylation may be an early event in tumor development that progress along a common pathway with BRCA1-mutated disease, representing a promising DNA-based biomarker for early-stage TNBC.

Homologous repair deficiency score for identifying breast cancers with defective DNA damage response

Breast cancer (BC) in patients with germline mutations of BRCA1/BRCA2 are associated with benefit from drugs targeting DNA damage response (DDR), but they account for only 5-7% of overall breast cancer. To define the characteristics of these tumors and also to identify tumors without BRCA mutation but with homologous recombination deficiency (HRD) is clinically relevant. To define characteristic features of HRD tumors and analyze the correlations between BRCA1/BRCA2 and BC subtypes, we analyzed 981 breast tumors from the TCGA database using the signature analyzer. The BRCA signature was strongly associated with the HRD score top 10% (score ≥ 57) population. This population showed a high level of mutations in DDR genes, including BRCA1/BRCA2. HRD tumors were associated with high expression levels of BARD1 and BRIP1. Besides, BRCA1/2 mutations were dominantly observed in basal and luminal subtypes, respectively. A comparison of HRD features in BC revealed that BRCA1 exerts a stronger influence inducing HRD features than BRCA2 does. It reveals genetic differences between BRCA1 and BRCA2 and provides a basis for the identification of HRD and other BRCA-associated tumors.

The Genomic Landscape of Mucinous Breast Cancer

Mucinous carcinoma of the breast (MCB) is a rare histologic form of estrogen receptor (ER)-positive/HER2-negative breast cancer (BC) characterized by tumor cells floating in lakes of mucin. We assessed the genomic landscape of 32 MCBs by whole-exome sequencing and/or RNA-sequencing. GATA3 (23.8%), KMT2C (19.0%), and MAP3K1 (14.3%) were the most frequently mutated genes in pure MCBs. In addition, two recurrent but not pathognomonic fusion genes, OAZ1-CSNK1G2 and RFC4-LPP, were detected in 3/31 (9.7%) and 2/31 (6.5%) samples, respectively. Compared with ER-positive/HER2-negative common forms of BC, MCBs displayed lower PIK3CA and TP53 mutation rates and fewer concurrent 1q gains and 16q losses. Clonal decomposition analysis of the mucinous and ductal components independently microdissected from five mixed MCBs revealed that they are clonally related and evolve following clonal selection or parallel evolution. Our findings indicate that MCB represents a genetically distinct ER-positive/HER2-negative form of BC.

Modeling and analysis of site-specific mutations in cancer identifies known plus putative novel hotspots and bias due to contextual sequences

In cancer, recurrently mutated sites in DNA and proteins, called hotspots, are thought to be raised by positive selection and therefore important due to its potential functional impact. Although recent evidence for APOBEC enzymatic activity have shown that specific types of sequences are likely to be false, the identification of putative hotspots is important to confirm either its functional role or its mechanistic bias. In this work, an algorithm and a statistical model is presented to detect hotspots. The model consists of a beta-binomial component plus fixed effects that efficiently fits the distribution of mutated sites. The algorithm employs an optimal stepwise approach to find the model parameters. Simulations show that the proposed algorithmic model is highly accurate for common hotspots. The approach has been applied to TCGA mutational data from 33 cancer types. The results show that well-known cancer hotspots are easily detected. Besides, novel hotspots are also detected. An analysis of the sequence context of detected hotspots show a preference for TCG sites that may be related to APOBEC or other unknown mechanistic biases. The detected hotspots are available online in http://bioinformatica.mty.itesm.mx/HotSpotsAnnotations.

Homologous recombination DNA repair deficiency and PARP inhibition activity in primary triple negative breast cancer

Triple negative breast cancer (TNBC) encompasses molecularly different subgroups, with a subgroup harboring evidence of defective homologous recombination (HR) DNA repair. Here, within a phase 2 window clinical trial, RIO trial (EudraCT 2014-003319-12), we investigate the activity of PARP inhibitors in 43 patients with untreated TNBC. The primary end point, decreased Ki67, occured in 12% of TNBC. In secondary end point analyses, HR deficiency was identified in 69% of TNBC with the mutational-signature-based HRDetect assay. Cancers with HRDetect mutational signatures of HR deficiency had a functional defect in HR, assessed by impaired RAD51 foci formation on end of treatment biopsy. Following rucaparib treatment there was no association of Ki67 change with HR deficiency. In contrast, early circulating tumor DNA dynamics identified activity of rucaparib, with end of treatment ctDNA levels suppressed by rucaparib in mutation-signature HR-deficient cancers. In ad hoc analysis, rucaparib induced expression of interferon response genes in HR-deficient cancers. The majority of TNBCs have a defect in DNA repair, identifiable by mutational signature analysis, that may be targetable with PARP inhibitors.

SNORA71B promotes breast cancer cells across blood-brain barrier by inducing epithelial-mesenchymal transition

Background

Brain metastasis (BM) is a dreadful complication that significantly impacts the quality of life in breast cancer patients. A key process during brain metastasis is the migration of cancer cells across blood–brain barrier (BBB). However, the role of snoRNAs regulating BBB in BM is still unknown.

Methods

Here SNORic and GEO databases were used to identify differentially expressed snoRNAs between brain metastatic and non-metastatic breast cancer (BC) tissues. The effects of SNORA71B on the capacities of proliferation, migration, invasion, epithelial-mesenchymal transition (EMT), and BBB invasion of BC cells were evaluated by CCK8, transwell, western blot, and BBB model, respectively.

Results

SNORA71B was highly expressed in high BM BC tissues and cells compared to low BM BC controls. Survival analysis revealed high expression of SNORA71B was significantly associated with poor PPS and OS in breast cancer patients. ROC curve showed that SNORA71B might act as biomarker for breast cancer. Moreover, SNORA71B significantly promoted proliferation, migration, and invasion of BC cells with different BM abilities. Importantly, SNORA71B promoted the EMT process of low BM BC cells. SNORA71B knockdown inhibited the high BM BC cells across BBB, while EMT activator dramatically abrogated this inhibited effect.

Conclusions

In conclusion, SNORA71B promotes BC cells across the BBB partly via inducing EMT.

Genetic alterations and their association with clinicopathologic characteristics in advanced breast carcinomas: focusing on clinically actionable genetic alterations

Breast carcinomas (BCs) are genetically heterogeneous and associated with numerous mutations which can be used to predict outcomes and initiate targeted therapies. We investigated clinicopathologic characteristics associated with gene mutations detected using the FoundationOne CDx assay in a cohort of 223 clinically advanced BCs (66 locally recurrent and 157 metastatic) from our institution. One hundred fifty unique mutations were identified (total 1008) in the cohort, with the most prevalent (>10%) including TP53 (53.8%), PIK3CA (35%), MYC (22%), CCND1 (19.7%), FGF19 (19.7%), FGF4 (16.6%), FGF3 (16.1%), ZNF703 (14.8%), ESR1 (13.9%), FGFR1 (13.5%), PTEN (12.1%), and CDH1 (10.8%). ERBB2 genetic alteration was most common in human epidermal growth factor receptor 2 (HER2)-positive BCs, and GATA3 and ESR1 mutations were only identified in hormone receptor-positive BC. Mutations enriched in triple-negative BCs (TNBCs) included TP53, PTEN, RB1, and CDKN2A/B. CDH1 mutation was predominantly found in lobular carcinomas, and PIK3CA mutation was also enriched. Mutations enriched in metaplastic carcinomas with heterologous mesenchymal differentiation included TP53, PTEN, MCL1, CDKN2A/B, and NOTCH2. An increase in mutations of CCND1, FGF19, FGF4, FGF3, ESR1, and EMSY was identified in metastatic BCs compared with locally recurrent BCs. Overall, PIK3CA was the most frequent clinically actionable genetic alteration (35%), followed by MYC (22%), CCND1 (19.7%), and FGF3/FGF4/FGFR1 (16%). In conclusion, our study provides genetic insight into the biology of advanced BCs and summarizes their most frequent clinically actionable genetic alterations, generating useful genomic information for potential improvement of patient management.

Clinical Implications of Chromosomal Instability (CIN) and Kinetochore Abnormalities in Breast Cancers

Genetic instability is a defining property of cancer cells and is the basis of various lesions including point mutations, copy number alterations and translocations. Chromosomal instability (CIN) is part of the genetic instability of cancer and consists of copy number alterations in whole or parts of cancer cell chromosomes. CIN is observed in differing degrees in most cancers. In breast cancer, CIN is commonly part of the genomic landscape of the disease and has a higher incidence in aggressive sub-types. Tumor suppressors that are commonly mutated or disabled in cancer, such as p53 and pRB, play roles in protection against CIN, and as a result, their dysfunction contributes to the establishment or tolerance of CIN. Several structural and regulatory proteins of the centromeres and kinetochore, the complex structure that is responsible for the correct distribution of genetic material in the daughter cells during mitosis, are direct or, mostly, indirect transcription targets of p53 and pRB. Thus, despite the absence of structural defects in genes encoding for centromere and kinetochore components, dysfunction of these tumor suppressors may have profound implications for the correct function of the mitotic apparatus contributing to CIN. CIN and its prognostic and therapeutic implications in breast cancer are discussed in this article.

DNA copy number motifs are strong and independent predictors of survival in breast cancer

Somatic copy number alterations are a frequent sign of genome instability in cancer. A precise characterization of the genome architecture would reveal underlying instability mechanisms and provide an instrument for outcome prediction and treatment guidance. Here we show that the local spatial behavior of copy number profiles conveys important information about this architecture. Six filters were defined to characterize regional traits in copy number profiles, and the resulting Copy Aberration Regional Mapping Analysis (CARMA) algorithm was applied to tumors in four breast cancer cohorts (n = 2919). The derived motifs represent a layer of information that complements established molecular classifications of breast cancer. A score reflecting presence or absence of motifs provided a highly significant independent prognostic predictor. Results were consistent between cohorts. The nonsite-specific occurrence of the detected patterns suggests that CARMA captures underlying replication and repair defects and could have a future potential in treatment stratification.