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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

GRHL2-controlled gene expression networks in luminal breast cancer

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

Abstract 3627: The genomic landscape of breast cancers from CHEK2 c.1100delC mutation carriers

Women who carry the germline CHEK2 c.1100delC mutation have a 2.3-fold increased life-time risk to develop breast cancer (BC). Their BCs are mostly of the luminal/ER+ subtype and diagnosed at a younger age than BCs in sporadic patients. BC patients carrying a CHEK2 c.1100delC mutation are at an increased risk of developing a contralateral BC and have a worse survival compared to sporadic BC patients although resistance to either endocrine or chemotherapy does not appear to play a role herein. To unravel the biological mechanism by which CHEK2 c.1100delC exerts its tumorigenic potential, whole genome sequencing was performed of matching normal and primary BC (pBC) DNA from twenty female CHEK2 c.1100delC mutation carriers at Hartwig Medical Foundation (HMF). We compared the somatic pBC genomes of these twenty patients to 560 pBC genomes from the International Cancer Genome Consortium including 30 BRCA1 and 25 BRCA2 mutation carriers and 329 sporadic ER+ pBC patients. The validation cohort consisted of 517 metastatic BC (mBC) genomes from the Center for Personalized Cancer Treatment that were sequenced at HMF and included 24 mBC genomes from CHEK2 c.1100delC mutation carriers. We found that CHEK2 pBC and mBC genomes were dissimilar to pBC and mBC genomes of BRCA1 and BRCA2 mutation carriers and did not show evidence of a homologous recombination deficiency (HRD) phenotype. Instead, CHEK2 pBC and mBC genomes were very similar to genomes of sporadic ER+ pBC and mBC patients in terms of tumor mutational burden, single and double base substitution signatures as well as small insertion/deletion signatures. The size distribution of structural variants (SVs), however, was different for CHEK2 pBCs and mBCs compared with ER+ as well as BRCA1 and BRCA2 pBCs and mBCs. Importantly, the CHEK2 c.1100delC mutation was associated with an increased frequency of larger (> 1Mb) deletions, inversions and tandem duplications (TDs) with peaks at specific sizes. Although CHEK2 mBC genomes more frequently harbored chromothripsis compared with other mBC subtypes, this did not explain the CHEK2-specific SV size distribution pattern. Moreover, increased chromothripsis in CHEK2 mBC genomes was not associated with a higher frequency of somatic TP53 mutations. In fact, we did not observe any somatic TP53 mutations among CHEK2 pBC and mBC genomes. This mutual exclusiveness suggests CHEK2-driven tumorigenesis could act through the TP53 signaling pathway. To conclude, CHEK2 BC genomes do not harbor somatic TP53 mutations and display a unique SV size distribution profile. Moreover, the absence of an HRD phenotype among CHEK2 BC genomes suggests PARP inhibitors will not be effective for treatment of BC in CHEK2 c.1100delC mutation carriers. Increasing knowledge of the mechanisms of CHEK2-driven tumorigenesis may ultimately lead to the development of tailored strategies for treatment and prevention of BC for CHEK2 mutation carriers.

Citation Format: Marcel Smid, Marjanka K. Schmidt, Wendy J. Prager - van der Smissen, Kirsten Ruigrok - Ritstier, Sten Cornelissen, Maartje A. Schreurs, J. Margriet Collee, Muriel A. Adank, Maartje J. Hooning, John W. Martens, Antoinette Hollestelle. The genomic landscape of breast cancers from CHEK2 c.1100delC mutation carriers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3627.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

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

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

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

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

Differential quantities of immune checkpoint-expressing CD8 T cells in soft tissue sarcoma subtypes

Introduction

Local T-cell immunity is recognized for its contribution to the evolution and therapy response of various carcinomas. Here, we investigated characteristics of tumor-infiltrating lymphocytes (TILs), as well as T-cell evasive mechanisms in different soft tissue sarcoma (STS) subtypes.

Methods

Liposarcoma, gastrointestinal stromal tumor (GIST), leiomyosarcoma, myxofibrosarcoma and pleomorphic sarcomas were assessed for T-cell numbers and phenotypes using flow cytometry. Next-generation sequencing was used to analyze T-cell receptor repertoire, mutational load, immune cell frequencies, and expression of immune-related genes.

Results

GIST, myxofibrosarcoma and pleomorphic sarcoma showed high numbers of CD8+ TILs, with GIST having the lowest fraction of effector memory T cells. These TILs coexpress the immune checkpoints PD1, TIM3, and LAG3 in myxofibrosarcoma and pleomorphic sarcoma, yet TILs coexpressing these checkpoints were near negligible in GIST. Fractions of dominant T-cell clones among STS subtypes were lowest in GIST and liposarcoma, whereas mutational load was relatively low in all STS subtypes. Furthermore, myeloid-derived cells and expression of the costimulatory ligands CD86, ICOS-L and 41BB-L were lowest in GIST when compared with other STS subtypes.

Conclusion

STS subtypes differ with respect to number and phenotypical signs of antitumor responsiveness of CD8+ TILs. Notably, GIST, myxofibrosarcoma and pleomorphic sarcoma harbor high numbers of CD8+ T cells, yet in the GIST microenvironment, these T cells are less differentiated and non-exhausted, which is accompanied with a relatively low expression of costimulatory ligands.

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

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

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

3545

Background: Consensus molecular subtypes (CMSs) can guide precision treatment of colorectal cancer (CRC). Currently available assays can identify CMS1 and CMS4 cases well, while a dedicated test to distinguish CMS2 and 3 is lacking. This study aimed to identify a panel of methylation markers to distinguish between CMS2 and 3 in patients with CRC. Methods: Fresh-frozen tumor tissue of 239 patients with stage I-III CRC was included. CMS classification was performed on RNA-seq data using the single-sample-prediction parameter from the “CMSclassifier” package. Methylation profiles were obtained using the Infinium HumanMethylation450 BeadChip. We performed adaptive group-regularised logistic ridge-regression with post-hoc group-weighted elastic net marker selection to build prediction models for classification of CMS2 and CMS3 based on 15, 10 or 5 markers. Data from TCGAwas used for validation. Results: Overall methylation profiles differed between CMS2 and CMS3. Group-regularisation of the probes was done based on their location either relative to a CpG island or relative to a gene present in the CMS classifier resulting in two different prediction models and subsequently different marker panels. For both panels, even when using only 5 markers, sensitivity, specificity, and accuracy were > 90%. Validation showed comparable performances. Conclusions: Our highly sensitive and specific methylation marker panel can be used to distinguish CMS2 and 3. This enables development of a qPCR DNA methylation assay in patients with CRC to provide a specific and non-invasive classification tool.

Circular RNA in Chemonaive Lymph Node Negative Colon Cancer Patients

Circular RNAs (circRNAs) appear important in tumor progression of colon cancer (CC). We identified an extensive catalog of circRNAs in 181 chemonaive stage I/II colon tumors, who underwent curative surgery between 2007 and 2014. We identified circRNAs from RNAseq data, investigated common biology related to circRNA expression, and studied the association between circRNAs and relapse status, tumor stage, consensus molecular subtypes (CMS), tumor localization and microsatellite instability (MSI). We identified 2606 unique circRNAs. 277 circRNAs (derived from 260 genes) were repeatedly occurring in at least 20 patients of which 153 showed a poor or even negative (R < 0.3) correlation with the expression level of their linear gene. The circular junctions for circSATB2, circFGD6, circKMT2C and circPLEKHM3 were validated by Sanger sequencing. Multiple correspondence analysis showed that circRNAs were often co-expressed and that high diversity in circRNAs was associated with favorable disease-free survival (DFS), which was confirmed by Cox regression analysis (Hazard Ratio (HR) 0.60, 95% CI 0.38-0.97, p = 0.036). Considering individual circRNAs, absence of circMGA was significantly associated with relapse, whereas circSATB2, circNAB1, and circCEP192 were associated with both MSI and CMS. This study represents a showcase of the potential clinical utility of circRNAs for prognostic stratification in patients with stage I-II colon cancer and demonstrated that high diversity in circRNAs is associated with favorable DFS.

Intratumoral heterogeneity of second-harmonic generation scattering from tumor collagen and its effects on metastatic risk prediction

Background

Metastases are the leading cause of breast cancer-related deaths. The tumor microenvironment impacts cancer progression and metastatic ability. Fibrillar collagen, a major extracellular matrix component, can be studied using the light scattering phenomenon known as second-harmonic generation (SHG). The ratio of forward- to backward-scattered SHG photons (F/B) is sensitive to collagen fiber internal structure and has been shown to be an independent prognostic indicator of metastasis-free survival time (MFS). Here we assess the effects of heterogeneity in the tumor matrix on the possible use of F/B as a prognostic tool.

Methods

SHG imaging was performed on sectioned primary tumor excisions from 95 untreated, estrogen receptor-positive, lymph node negative invasive ductal carcinoma patients. We identified two distinct regions whose collagen displayed different average F/B values, indicative of spatial heterogeneity: the cellular tumor bulk and surrounding tumor-stroma interface. To evaluate the impact of heterogeneity on F/B’s prognostic ability, we performed SHG imaging in the tumor bulk and tumor-stroma interface, calculated a 21-gene recurrence score (surrogate for OncotypeDX®, or S-ODX) for each patient and evaluated their combined prognostic ability.

Results

We found that F/B measured in tumor-stroma interface, but not tumor bulk, is prognostic of MFS using three methods to select pixels for analysis: an intensity threshold selected by a blinded observer, a histogram-based thresholding method, and an adaptive thresholding method. Using both regression trees and Random Survival Forests for MFS outcome, we obtained data-driven prediction rules that show F/B from tumor-stroma interface, but not tumor bulk, and S-ODX both contribute to predicting MFS in this patient cohort. We also separated patients into low-intermediate (S-ODX < 26) and high risk (S-ODX ≥26) groups. In the low-intermediate risk group, comprised of patients not typically recommended for adjuvant chemotherapy, we find that F/B from the tumor-stroma interface is prognostic of MFS and can identify a patient cohort with poor outcomes.

Conclusions

These data demonstrate that intratumoral heterogeneity in F/B values can play an important role in its possible use as a prognostic marker, and that F/B from tumor-stroma interface of primary tumor excisions may provide useful information to stratify patients by metastatic risk.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-020-07713-4.

Interconnectivity between molecular subtypes and tumor stage in colorectal cancer

BACKGROUND

There are profound individual differences in clinical outcomes between colorectal cancers (CRCs) presenting with identical stage of disease. Molecular stratification, in conjunction with the traditional TNM staging, is a promising way to predict patient outcomes. We investigated the interconnectivity between tumor stage and tumor biology reflected by the Consensus Molecular Subtypes (CMSs) in CRC, and explored the possible value of these insights in patients with stage II colon cancer.

METHODS

We performed a retrospective analysis using clinical records and gene expression profiling in a meta-cohort of 1040 CRC patients. The interconnectivity of tumor biology and disease stage was assessed by investigating the association between CMSs and TNM classification. In order to validate the clinical applicability of our findings we employed a meta-cohort of 197 stage II colon cancers.

RESULTS

CMS4 was significantly more prevalent in advanced stages of disease (stage I 9.8% versus stage IV 38.5%, p < 0.001). The observed differential gene expression between cancer stages is at least partly explained by the biological differences as reflected by CMS subtypes. Gene signatures for stage III-IV and CMS4 were highly correlated (r = 0.77, p < 0.001). CMS4 cancers showed an increased progression rate to more advanced stages (CMS4 compared to CMS2: 1.25, 95% CI: 1.08-1.46). Patients with a CMS4 cancer had worse survival in the high-risk stage II tumors compared to the total stage II cohort (5-year DFS 41.7% versus 100.0%, p = 0.008).

CONCLUSIONS

Considerable interconnectivity between tumor biology and tumor stage in CRC exists. This implies that the TNM stage, in addition to the stage of progression, might also reflect distinct biological disease entities. These insights can potentially be utilized to optimize identification of high-risk stage II colon cancers.

Nature of fiber-coupled detector responsivity measurements at 0.1% using a primary standard

We demonstrate the capability to measure the absolute power responsivity of optical fiber-coupled detectors at an expanded uncertainty of 0.1%, by direct comparison with a cryogenic primary standard. To facilitate synchronous power measurements, commercial all-fiber beam-splitters direct laser diode light simultaneously to the device under test and the primary standard. We investigate the use of single-mode, polarisation maintaining, and photonic crystal fibers to access the cryogenic standard, and report a reduction in the temperature dependent effective refractive index of these fibers of 0.1%, 0.15% and 0.3% respectively in going from room temperature to 5 K. We also evaluate the polarisation dependent loss of the beam-splitters, the stability of the beam-splitter ratio between the cryogenic detector and the device under test and the temporal and modal stability of the Fabry-Pérot laser diode sources. It is shown that the stability of the optical fiber beam-splitters limits the overall performance of the measurement system to an expanded uncertainty of 0.1%.

Abstract P3-03-03: Inflammatory breast cancer exhibits amplification and transcriptional activation of MYC in conjunction with NOTCH and TGFβ signaling

Inflammatory breast cancer (IBC) is an aggressive and poorly characterized breast cancer subtype with few targets available for molecular therapy. Several attempts aimed at characterizing the disease at the genomic, transcriptomic and proteomic level. Together, these studies suggest that MYC is central to the biology of IBC. Here, we analyzed data sets from patient samples and preclinical models to corroborate this hypothesis.

Genomic and transcriptomic profiles were obtained using next-generation sequencing and Affymetrix GeneChip analysis on a series of 101 and 137 samples from patients with IBC. In addition, non-IBC (nIBC) control samples were included from our own repositories (N=252) as well as from publicly available data sets (i.e. TCGA and METABRIC). Finally, RNA-sequencing and peptide phosphorylation data from preclinical IBC and nIBC models, both with and without TGFβ treatment, were explored.

Apart from TP53 and ERBB2, MYC is the 3rd most frequently mutated gene with single nucleotide variants, indels and amplifications found in 22% of 101 patients. Furthermore, MYC mutations are mutually exclusive with alterations in the NOTCH genes (i.e. NOTCH1, -2, -3, and -4 and CREBBP). At the gene expression level, a MYC activation score is elevated in IBC as compare to nIBC and correlates with a NOTCH activation signature. Genes overexpressed in IBC are enriched for various sets of MYC target genes and are associated with MYC amplifications in TCGA and METABRIC samples. In preclinical models, MYC target genes are particularly responsive to TGFβ treatment in IBC cells but not in nIBC cells and analysis of peptide phosphorylation patterns suggest this interaction may be regulated via PTK2B.

The present data suggest that MYC could be an important driver of IBC biology, which corroborates with the recent observation that MYC is an important hub in the signal transduction network of SUM149 and -190 cells. Furthermore, our data reveal that MYC activation is linked to NOCTH and TGFβ signaling. In earlier studies it was demonstrated that lymphovascular tumor emboli in the Mary-X mouse model of IBC exhibit NOTCH3 addition with downstream MYC overexpression, indicating that NOTCH/MYC signaling could be crucial for metastatic progression in IBC. In this context, it should be noted that an allele associated with metastatic risk in IBC patients resides at the 8q24 locus that harbors the MYC gene.

Citation Format: Charlotte Rypens, François Bertucci, Patrice Viens, John WM Martens, Marcel Smid, Carolien Schröder, Naoto T Ueno, Massimo Cristofanilli, Piet Dirix, Peter Vermeulen, Luc Dirix, Steven Van Laere. Inflammatory breast cancer exhibits amplification and transcriptional activation of MYC in conjunction with NOTCH and TGFβ signaling [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P3-03-03.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Comprehensive evaluation of methods to assess overall and cell-specific immune infiltrates in breast cancer

Background

Breast cancer (BC) immune infiltrates play a critical role in tumor progression and response to treatment. Besides stromal tumor infiltrating lymphocytes (sTILs) which have recently reached level 1B evidence as a prognostic marker in triple negative BC, a plethora of methods to assess immune infiltration exists, and it is unclear how these compare to each other and if they can be used interchangeably.

Methods

Two experienced pathologists scored sTIL, intra-tumoral TIL (itTIL), and 6 immune cell types (CD3⁺, CD4⁺, CD8⁺, CD20⁺, CD68⁺, FOXP3⁺) in the International Cancer Genomics Consortium breast cancer cohort using hematoxylin and eosin-stained (n = 243) and immunohistochemistry-stained tissue microarrays (n = 254) and whole slides (n = 82). The same traits were evaluated using transcriptomic- and methylomic-based deconvolution methods or signatures.

Results

The concordance correlation coefficient (CCC) between pathologists for sTIL was very good (0.84) and for cell-specific immune infiltrates slightly lower (0.63–0.66). Comparison between tissue microarray and whole slide pathology scores revealed systematically higher values in whole slides (ratio 2.60–5.98). The Spearman correlations between microscopic sTIL and transcriptomic- or methylomic-based assessment of immune infiltrates were highly variable (r = 0.01–0.56). Similar observations were made for cell type-specific quantifications (r = 0.001–0.54). We observed a strong inter-method variability between the omics-derived estimations, which is further cell type dependent. Finally, we demonstrated that most methods more accurately identify highly infiltrated (sTIL ≥ 60%; area under the curve, AUC, 0.64–0.99) as compared to lowly infiltrated tumors (sTIL ≤ 10%; AUC 0.52–0.82).

Conclusions

There is a lower inter-pathologist concordance for cell-specific quantification as compared to overall infiltration quantification. Microscopic assessments are underestimated when considering small cores (tissue microarray) instead of whole slides. Results further highlight considerable differences between the microscopic-, transcriptomic-, and methylomic-based methods in the assessment of overall and cell-specific immune infiltration in BC. We therefore call for extreme caution when assessing immune infiltrates using current methods and emphasize the need for standardized immune characterization beyond TIL.

Clonality, Antigen Recognition, and Suppression of CD8+ T Cells Differentially Affect Prognosis of Breast Cancer Subtypes

PURPOSE

In breast cancer, response rates to immune therapies are generally low and differ significantly across molecular subtypes, urging a better understanding of immunogenicity and immune evasion.

EXPERIMENTAL DESIGN

We interrogated large gene-expression data sets including 867 node-negative, treatment-naïve breast cancer patients (microarray data) and 347 breast cancer patients (whole-genome sequencing and transcriptome data) according to parameters of T cells as well as immune microenvironment in relation to patient survival.

RESULTS

We developed a 109-immune gene signature that captures abundance of CD8 tumor-infiltrating lymphocytes (TIL) and is prognostic in basal-like, her2, and luminal B breast cancer, but not in luminal A or normal-like breast cancer. Basal-like and her2 are characterized by highest CD8 TIL abundance, highest T-cell clonality, highest frequencies of memory T cells, and highest antigenicity, yet only the former shows highest expression level of immune and metabolic checkpoints and highest frequency of myeloid suppressor cells. Also, luminal B shows a high antigenicity and T-cell clonality, yet a low abundance of CD8 TILs. In contrast, luminal A and normal-like both show a low antigenicity, and notably, a low and high abundance of CD8 TILs, respectively, which associates with T-cell influx parameters, such as expression of adhesion molecules.

CONCLUSIONS

Collectively, our data argue that not only CD8 T-cell presence itself, but rather T-cell clonality, T-cell subset distribution, coinhibition, and antigen presentation reflect occurrence of a CD8 T-cell response in breast cancer subtypes, which have been aborted by distinct T-cell-suppressive mechanisms, providing a rationale for subtype-specific combination immune therapies.

The genomic landscape of metastatic breast cancer highlights changes in mutation and signature frequencies

The whole-genome sequencing of prospectively collected tissue biopsies from 442 patients with metastatic breast cancer reveals that, compared to primary breast cancer, tumor mutational burden doubles, the relative contributions of mutational signatures shift and the mutation frequency of six known driver genes increases in metastatic breast cancer. Significant associations with pretreatment are also observed. The contribution of mutational signature 17 is significantly enriched in patients pretreated with fluorouracil, taxanes, platinum and/or eribulin, whereas the de novo mutational signature I identified in this study is significantly associated with pretreatment containing platinum-based chemotherapy. Clinically relevant subgroups of tumors are identified, exhibiting either homologous recombination deficiency (13%), high tumor mutational burden (11%) or specific alterations (24%) linked to sensitivity to FDA-approved drugs. This study provides insights into the biology of metastatic breast cancer and identifies clinically useful genomic features for the future improvement of patient management.

Cancer-associated fibroblast-derived Gremlin 1 promotes breast cancer progression

Background

Bone morphogenetic proteins (BMPs) have been reported to maintain epithelial integrity and to antagonize the transforming growth factor β (TGFβ)-induced epithelial to mesenchymal transition. The expression of soluble BMP antagonists is dysregulated in cancers and interrupts proper BMP signaling in breast cancer.

Methods

In this study, we mined the prognostic role of BMP antagonists GREMLIN 1 (GREM1) in primary breast cancer tissues using in-house and publicly available datasets. We determined which cells express GREM1 RNA using in situ hybridization (ISH) on a breast cancer tissue microarray. The effects of Grem1 on the properties of breast cancer cells were assessed by measuring the mesenchymal/stem cell marker expression and functional cell-based assays for stemness and invasion. The role of Grem1 in breast cancer-associated fibroblast (CAF) activation was measured by analyzing the expression of fibroblast markers, phalloidin staining, and collagen contraction assays. The role of Grem1 in CAF-induced breast cancer cell intravasation and extravasation was studied by utilizing xenograft zebrafish breast cancer (co-) injection models.

Results

Expression analysis of clinical breast cancer datasets revealed that high expression of GREM1 in breast cancer stroma is correlated with a poor prognosis regardless of the molecular subtype. The large majority of human breast cancer cell lines did not express GREM1 in vitro, but breast CAFs did express GREM1 both in vitro and in vivo. Transforming growth factor β (TGFβ) secreted by breast cancer cells, and also inflammatory cytokines, stimulated GREM1 expression in CAFs. Grem1 abrogated bone morphogenetic protein (BMP)/SMAD signaling in breast cancer cells and promoted their mesenchymal phenotype, stemness, and invasion. Moreover, Grem1 production by CAFs strongly promoted the fibrogenic activation of CAFs and promoted breast cancer cell intravasation and extravasation in co-injection xenograft zebrafish models.

Conclusions

Our results demonstrated that Grem1 is a pivotal factor in the reciprocal interplay between breast cancer cells and CAFs, which promotes cancer cell invasion. Targeting Grem1 could be beneficial in the treatment of breast cancer patients with high Grem1 expression.

Electronic supplementary material

The online version of this article (10.1186/s13058-019-1194-0) contains supplementary material, which is available to authorized users.

Circulating Tumor Cell Enumeration and Characterization in Metastatic Castration-Resistant Prostate Cancer Patients Treated with Cabazitaxel

(1) Background: Markers identifying which patients with metastatic, castration-resistant prostate cancer (mCRPC) will benefit from cabazitaxel therapy are currently lacking. Therefore, the aim of this study was to identify markers associated with outcome to cabazitaxel therapy based on counts and gene expression profiles of circulating tumor cells (CTCs). (2) Methods: From 120 mCRPC patients, CellSearch enriched CTCs were obtained at baseline and after 6 weeks of cabazitaxel therapy. Furthermore, 91 genes associated with prostate cancer were measured in mRNA of these CTCs. (3) Results: In 114 mCRPC patients with an evaluable CTC count, the CTC count was independently associated with poor progression-free survival (PFS) and overall survival (OS) in multivariable analysis with other commonly used variables associated with outcome in mCRPC (age, prostate specific antigen (PSA), alkaline phosphatase, lactate dehydrogenase (LDH), albumin, hemoglobin), together with alkaline phosphatase and hemoglobin. A five-gene expression profile was generated to predict for outcome to cabazitaxel therapy. However, even though this signature was associated with OS in univariate analysis, this was not the case in the multivariate analysis for OS nor for PFS. (4) Conclusion: The established five-gene expression profile in CTCs was not independently associated with PFS nor OS. However, along with alkaline phosphatase and hemoglobin, CTC-count is independently associated with PFS and OS in mCRPC patients who are treated with cabazitaxel.

Downregulation of matrix Gla protein is a biomarker for tamoxifen-resistant and radioresistant breast cancer

Aim: To establish downregulated biomarkers for cross-resistance for radiotherapy (RT) in tamoxifen (TAM)-resistant breast cancer cells. Materials & methods: RNA sequencing was performed on TAM and RT-resistant breast cancer cells. Breast cancer patient cohorts were queried in silico for associations between genes of interest and outcome after TAM treatment or irradiation. Results: 20 genes showed decreased expression in both TAM-resistant and RT-resistant breast cancer cells. Only matrix Gla protein in the primary tumor was associated with outcome after TAM treatment or RT in breast cancer patients. Conclusion: Matrix Gla protein is a biomarker for therapy sensitivity in breast cancer.

A kinase inhibitor screen identifies a dual cdc7/CDK9 inhibitor to sensitise triple-negative breast cancer to EGFR-targeted therapy

BACKGROUND

The effective treatment of triple-negative breast cancer (TNBC) remains a profound clinical challenge. Despite frequent epidermal growth factor receptor (EGFR) overexpression and reliance on downstream signalling pathways in TNBC, resistance to EGFR-tyrosine kinase inhibitors (TKIs) remains endemic. Therefore, the identification of targeted agents, which synergise with current therapeutic options, is paramount.

METHODS

Compound-based, high-throughput, proliferation screening was used to profile the response of TNBC cell lines to EGFR-TKIs, western blotting and siRNA transfection being used to examine the effect of inhibitors on EGFR-mediated signal transduction and cellular dependence on such pathways, respectively. A kinase inhibitor combination screen was used to identify compounds that synergised with EGFR-TKIs in TNBC, utilising sulphorhodamine B (SRB) assay as read-out for proliferation. The impact of drug combinations on cell cycle arrest, apoptosis and signal transduction was assessed using flow cytometry, automated live-cell imaging and western blotting, respectively. RNA sequencing was employed to unravel transcriptomic changes elicited by this synergistic combination and to permit identification of the signalling networks most sensitive to co-inhibition.

RESULTS

We demonstrate that a dual cdc7/CDK9 inhibitor, PHA-767491, synergises with multiple EGFR-TKIs (lapatinib, erlotinib and gefitinib) to overcome resistance to EGFR-targeted therapy in various TNBC cell lines. Combined inhibition of EGFR and cdc7/CDK9 resulted in reduced cell proliferation, accompanied by induction of apoptosis, G2-M cell cycle arrest, inhibition of DNA replication and abrogation of CDK9-mediated transcriptional elongation, in contrast to mono-inhibition. Moreover, high expression of cdc7 and RNA polymerase II Subunit A (POLR2A), the direct target of CDK9, is significantly correlated with poor metastasis-free survival in a cohort of breast cancer patients. RNA sequencing revealed marked downregulation of pathways governing proliferation, transcription and cell survival in TNBC cells treated with the combination of an EGFR-TKI and a dual cdc7/CDK9 inhibitor. A number of genes enriched in these downregulated pathways are associated with poor metastasis-free survival in TNBC.

CONCLUSIONS

Our results highlight that dual inhibition of cdc7 and CDK9 by PHA-767491 is a potential strategy for targeting TNBC resistant to EGFR-TKIs.

MicroRNAs as possible indicators of drug sensitivity in breast cancer cell lines

MicroRNAs (miRNAs) regulate gene expression post-transcriptionally. In this way they might influence whether a cell is sensitive or resistant to a certain drug. So far, only a limited number of relatively small scale studies comprising few cell lines and/or drugs have been performed. To obtain a broader view on miRNAs and their association with drug response, we investigated the expression levels of 411 miRNAs in relation to drug sensitivity in 36 breast cancer cell lines. For this purpose IC₅₀ values of a drug screen involving 34 drugs were associated with miRNA expression data of the same breast cancer cell lines. Since molecular subtype of the breast cancer cell lines is considered a confounding factor in drug association studies, multivariate analysis taking subtype into account was performed on significant miRNA-drug associations which retained 13 associations. These associations consisted of 11 different miRNAs and eight different drugs (among which Paclitaxel, Docetaxel and Veliparib). The taxanes, Paclitaxel and Docetaxel, were the only drugs having miRNAs in common: hsa-miR-187-5p and hsa-miR-106a-3p indicative of drug resistance while Paclitaxel sensitivity alone associated with hsa-miR-556-5p. Tivantinib was associated with hsa-let-7d-5p and hsa-miR-18a-5p for sensitivity and hsa-miR-637 for resistance. Drug sensitivity was associated with hsa-let-7a-5p for Bortezomib, hsa-miR-135a-3p for JNJ-707 and hsa-miR-185-3p for Panobinostat. Drug resistance was associated with hsa-miR-182-5p for Veliparib and hsa-miR-629-5p for Tipifarnib. Pathway analysis for significant miRNAs was performed to reveal biological roles, aiding to find a potential mechanistic link for the observed associations with drug response. By doing so hsa-miR-187-5p was linked to the cell cycle G2-M checkpoint in line with this checkpoint being the target of taxanes. In conclusion, our study shows that miRNAs could potentially serve as biomarkers for intrinsic drug resistance and that pathway analyses can provide additional information in this context.

Partially methylated domains are hypervariable in breast cancer and fuel widespread CpG island hypermethylation

Global loss of DNA methylation and CpG island (CGI) hypermethylation are key epigenomic aberrations in cancer. Global loss manifests itself in partially methylated domains (PMDs) which extend up to megabases. However, the distribution of PMDs within and between tumor types, and their effects on key functional genomic elements including CGIs are poorly defined. We comprehensively show that loss of methylation in PMDs occurs in a large fraction of the genome and represents the prime source of DNA methylation variation. PMDs are hypervariable in methylation level, size and distribution, and display elevated mutation rates. They impose intermediate DNA methylation levels incognizant of functional genomic elements including CGIs, underpinning a CGI methylator phenotype (CIMP). Repression effects on tumor suppressor genes are negligible as they are generally excluded from PMDs. The genomic distribution of PMDs reports tissue-of-origin and may represent tissue-specific silent regions which tolerate instability at the epigenetic, transcriptomic and genetic level.

The circular RNome of primary breast cancer

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

Laser-wakefield accelerators for high-resolution X-ray imaging of complex microstructures

Laser-wakefield accelerators (LWFAs) are high acceleration-gradient plasma-based particle accelerators capable of producing ultra-relativistic electron beams. Within the strong focusing fields of the wakefield, accelerated electrons undergo betatron oscillations, emitting a bright pulse of X-rays with a micrometer-scale source size that may be used for imaging applications. Non-destructive X-ray phase contrast imaging and tomography of heterogeneous materials can provide insight into their processing, structure, and performance. To demonstrate the imaging capability of X-rays from an LWFA we have examined an irregular eutectic in the aluminum-silicon (Al-Si) system. The lamellar spacing of the Al-Si eutectic microstructure is on the order of a few micrometers, thus requiring high spatial resolution. We present comparisons between the sharpness and spatial resolution in phase contrast images of this eutectic alloy obtained via X-ray phase contrast imaging at the Swiss Light Source (SLS) synchrotron and X-ray projection microscopy via an LWFA source. An upper bound on the resolving power of 2.7 ± 0.3 μm of the LWFA source in this experiment was measured. These results indicate that betatron X-rays from laser wakefield acceleration can provide an alternative to conventional synchrotron sources for high resolution imaging of eutectics and, more broadly, complex microstructures.

Androgen receptor expression in circulating tumor cells of patients with metastatic breast cancer

The androgen receptor (AR) has potential clinical relevance in metastatic breast cancer (mBC) since it might be a treatment target and has been associated with endocrine resistance. A minimal-invasive way to determine AR expression on metastatic tumor cells is by characterization of circulating tumor cells (CTCs). Here, we assessed AR mRNA expression in CTCs (CTC-AR) and in matched primary tumor samples from mBC patients representing different breast cancer subtypes. In addition, we explored CTC-AR-status in relation to outcome on endocrine therapy. AR, and 92 AR or estrogen receptor (ER) related genes, were measured in CellSearch-enriched CTCs from 124 mBC patients and in 52 matched FFPE primary tissues using quantitative reverse-transcriptase PCR. AR in CTCs was considered positive if the expression was 1 standard deviation higher than the expression measured in 11 healthy blood donors. A total of 31% of the mBC patients had AR-positive (AR+) CTCs. 58% of the matched CTC and primary tumor samples were discordant with respect to AR status, observing both switches from AR+ to AR-negative (AR-) and vice versa. There was no statistically significant difference in progression-free survival for patients treated with ER-targeting drugs and CTC-AR-status (13 AR+/ 37 AR- cases, p = 0.28). Thus, AR can be determined in RNA isolated from CTCs, with in our set 31% AR-positive samples. Given the discordance between AR status in CTC samples and corresponding primary tumors, determination of AR expression in CTCs might be a promising tool to select mBC patients for AR inhibiting agents.

Mitochondrial RNA Expression and Single Nucleotide Variants in Association with Clinical Parameters in Primary Breast Cancers

The human mitochondrial DNA (mtDNA) encodes 37 genes, including thirteen proteins essential for the respiratory chain, and RNAs functioning in the mitochondrial translation apparatus. The total number of mtDNA molecules per cell (mtDNA content) is variable between tissue types and also between tumors and their normal counterparts. For breast cancer, tumors tend to be depleted in their mtDNA content compared to adjacent normal mammary tissue. Various studies have shown that primary breast tumors harbor somatic mtDNA variants. A decrease in mtDNA content or the presence of somatic variants could indicate a reduced mitochondrial function within breast cancer. In this explorative study we aimed to further understand genomic changes and expression of the mitochondrial genome within breast cancer, by analyzing RNA sequencing data of primary breast tumor specimens of 344 cases. We demonstrate that somatic variants detected at the mtRNA level are representative for somatic variants in the mtDNA. Also, the number of somatic variants within the mitochondrial transcriptome is not associated with mutational processes impacting the nuclear genome, but is positively associated with age at diagnosis. Finally, we observe that mitochondrial expression is related to ER status. We conclude that there is a large heterogeneity in somatic mutations of the mitochondrial genome within primary breast tumors, and differences in mitochondrial expression among breast cancer subtypes. The exact impact on metabolic differences and clinical relevance deserves further study.

Abstract 4978: On circular RNAs in breast cancer

Introduction Circular RNA was already discovered many decades ago and circRNAs were long considered idiosyncrasies of the splicing machinery processing precursor mRNA to mature mRNA. A much more recent take on this phenomenon showed an unanticipated abundance of circRNA isoforms in human cells. In breast cancer some work has already been done on publicly available TCGA data, but this analysis has a big limitation since the RNA-seq data were prepared using a poly(A) selection step, thereby missing most circRNAs (which lack a poly(A) tail). Here we describe the identification of circRNAs in a large cohort of 348 primary breast tumors, using RNA-seq data obtained without using a poly(A) selection step in the generation of the sequence libraries. We developed a method that, in contrast to previous identification methods, does not rely on unmapped reads or known splice-junctions.

Methods RNA-seq data were generated, using total RNA and removing ribosomal RNA using Ribo Zero (Illumina, USA). Paired-end (75 bases) sequencing was performed on an Illumina HiSeq 2000. The resulting fastq files were mapped to GRCh38 using STAR (version 2.4.2a) and the resulting bam files were sorted and indexed. CircRNAs were identified directly on the bam file generated by STAR; assuming a circular RNA molecule is present, the sequence read that aligns over the crossing of the non-canonical junction would, on a linear reference, get two different mapped locations. The read-mate aligns somewhere in between these two locations with the head of the read-mate pointing toward the head of the junction read, as would be expected in a properly paired read. Finding additional read-pairs showing this configuration, but always with a breakpoint at the exact same location, strengthens the evidence for circular transcripts.

Results In our cohort of 348 primary breast cancer cases we identified 25,783 circRNAs having a junction over exons of a gene. Of these, 13,853 (54%) were identified in at least 2 samples, and 426 circRNAs were found in at least 150 samples. The top recurrent circRNA was CDR1, a known circRNA that functions as a miR-7 sponge. Other previously described circRNAs were also present in our list. Of the circRNAs with sufficient data, 668 (41%) showed a poor (R&lt;0.2) and 210 (13%) showed a negative correlation with the expression level of the linear gene. We validated several circRNAs by PCR and sequencing, including CNOT2 and CREBBP, the former showing prognostic value, the latter a known breast cancer driver gene. Additional analyses are ongoing and will be presented at the annual meeting.

Conclusion circRNAs are an abundant RNA species in breast cancer, several of which may just be a remnant of the splicing machinery. However, several circRNAs have characteristics of being regulated and thus may have a functional role in breast cancer. Due to their structural stability, candidate circRNAs may serve as good biomarkers for diagnosis and therapy-monitoring in liquid biopsies.

Citation Format: Marcel Smid, Vanja de Weerd, Saskia Wilting, Anieta Sieuwerts, Serena Nik-Zainal, John Martens. On circular RNAs in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4978.

Abstract 4345: Molecular comparison of imatinib-naïve and resistant gastrointestinal stromal tumors: Differentially expressed microRNAs and mRNAs

Background and aim: Gastrointestinal stromal tumors (GISTs) are rare mesenchymal tumors of the gastrointestinal tract. Gain-of-function mutations in KIT or, less frequently, PDGFRA underlie the pathogenesis of the majority of these tumors. The introduction of imatinib mesylate, a selective tyrosine kinase inhibitor, has dramatically improved the outcome for GIST patients. However, half of the advanced patients treated with imatinib eventually progress and acquire resistance within two years of treatment, underscoring the need to get better insight into the resistance mechanisms. In this context we examined the miRNA and mRNA expression profiles in primary (imatinib-naïve) and imatinib resistant GIST samples.

Material and Methods: Fifty-three frozen GIST samples derived from various anatomical sites (small intestine, stomach, colon) harboring a KIT mutation including exon 9 (n=11), exon 11 (n=41) and exon 17 (n=1), were analyzed. Total RNA was isolated from imatinib naïve (IM-n; n=33) and resistant (IM-r; n=20) tumors. The resistant tumors bearing secondary mutations in exon 13 (n=3), exon 17 (n=6) or no secondary mutations (n=11). The miRNA expression profiles were determined using LNA™ oligonucleotide arrays (Exiqon) capable of detecting 725 human miRNAs. Furthermore, from a subset of samples (IM-n; n=14; IM-r; n=15) the mRNA expression profile was established using Affymetrix U133A oligonucleotide arrays. Using the mRNA expression data and the Biocarta and KEGG databases, global testing identified a number of biochemical pathways containing genes that are differentially expressed between the primary and progressive tumors.

Results: Thirty-five differentially expressed miRNAs (P&lt;0.009, FDR&lt;20%) between IM-n and IM-r GIST samples were identified. MiR-30c, miR-181a and miR-144 were among the most significantly different, with p values of 0.006, 0.006 and 0.003, respectively. Analysis of the mRNA data has revealed the up-regulation of genes (P &lt;0002, FDR&lt;30%) involved in cell cycle, DNA replication and cancer related pathways (i.e. proliferation, genomic instability, resistance to chemotherapy, evading apoptosis, etc.) in IM-r samples. To further elucidate the possible involvement of the deregulated miRNAs in the regulation of the indicated genes, the most significantly differentially expressed genes were correlated to miRNA expression. Negative correlations indicate miRNA mediated gene regulation, and were further examined in target prediction programs.

Conclusion: A molecular comparison between primary, IM-n and IM-r GIST samples revealed significant differentially expressed miRNAs and mRNAs. Bioinformatic analyses provided insight into biochemical pathways, and their putative regulation by miRNAs, that may contribute to imatinib resistance as observed in GIST patients.

Citation Format: Azadeh Amirnasr, Caroline M. Gits, Patricia F. Kuijk, Marcel Smid, Maria Debiec-Rychter, Stefan Sleijfer, Erik A. Wiemer. Molecular comparison of imatinib-naïve and resistant gastrointestinal stromal tumors: Differentially expressed microRNAs and mRNAs [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4345.

Association of microRNA-7 and its binding partner CDR1-AS with the prognosis and prediction of 1st-line tamoxifen therapy in breast cancer

The large number of non-coding RNAs (ncRNAs) and their breadth of functionalities has fuelled many studies on their roles in cancer. We previously linked four microRNAs to breast cancer prognosis. One of these microRNAs, hsa-miR-7, was found to be regulated by another type of ncRNA, the circular non-coding RNA (circRNA) CDR1-AS, which contains multiple hsa-miR-7 binding sites. Based on this finding, we studied the potential clinical value of this circRNA on breast cancer prognosis in a cohort based on a cohort that was previously analysed for hsa-miR-7 and in an adjuvant hormone-naïve cohort for 1^st-line tamoxifen treatment outcomes, in which we also analysed hsa-miR-7. A negative correlation was observed between hsa-miR-7 and CDR1-AS in both cohorts. Despite associations with various clinical metrics (e.g., tumour grade, tumour size, and relapse location), CDR1-AS was neither prognostic nor predictive of relevant outcomes in our cohorts. However, we did observe stromal CDR1-AS expression, suggesting a possible cell-type specific interaction. Next to the known association of hsa-miR-7 expression with poor prognosis in primary breast cancer, we found that high hsa-miR-7 expression was predictive of an adverse response to tamoxifen therapy and poor progression-free and post-relapse overall survival in patients with recurrent disease.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons

BACKGROUND

Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data.

RESULTS

We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality.

CONCLUSIONS

We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Epithelial-Mesenchymal Transition in Human Prostate Cancer Demonstrates Enhanced Immune Evasion Marked by IDO1 Expression

Cancer invasion and metastasis are driven by epithelial-mesenchymal transition (EMT), yet the exact mechanisms that account for EMT in clinical prostate cancer are not fully understood. Expression of N-cadherin is considered a hallmark of EMT in clinical prostate cancer. In this study, we determined the molecular mechanisms associated with N-cadherin expression in patients with prostate cancer. We performed laser capture microdissection of matched N-cadherin-positive and -negative prostate cancer areas from patient samples (n = 8), followed by RNA sequencing. N-cadherin expression was significantly associated with an immune-regulatory signature including profound upregulation of indoleamine 2,3-dioxygenase (IDO1; log₂-fold change = 5.1; P = 2.98E-04). Fluorescent immunostainings of patient samples confirmed expression of IDO1 protein and also its metabolite kynurenine in primarily N-cadherin-positive areas. N-cadherin-positive areas also exhibited a local decrease of intraepithelial cytotoxic (CD8⁺) T cells and an increase of immunosuppressive regulatory T cells (CD4⁺/FOXP3⁺). In conclusion, EMT in clinical prostate cancer is accompanied by upregulated expression of IDO1 and an increased number of regulatory T cells. These data indicate that EMT, which is an important step in tumor progression, can be protected from effective immune control in patients with prostate cancer.Significance: These findings demonstrate EMT is linked to an immunosuppressive environment in clinical prostate cancer, suggesting that patients with prostate cancer can potentially benefit from combinatorial drug therapy. Cancer Res; 78(16); 4671-9. ©2018 AACR.

MicroRNA expression in pre-treatment plasma of patients with benign breast diseases and breast cancer

Background

MicroRNAs (miRs) are small RNA molecules, influencing messenger RNA (mRNA) expression and translation, and are readily detectable in blood. Some have been reported as potential breast cancer biomarkers. This study aimed to identify and validate miRs indicative of breast cancer.

Results

Based on the discovery and literature, 18 potentially informative miRs were quantified in the validation cohort. Irrespective of patient and tumour characteristics, hsa-miR-652-5p was significantly upregulated in the malignant compared to benign patients (1.26 fold, P = 0.005) and therefore validated as potential biomarker. In the validation cohort literature-based hsa-let-7b levels were higher in malignant patients as well (1.53 fold, P = 0.011). Two miRs differentiated benign wildtype from benign BRCA1 mutation carriers and an additional 8 miRs differentiated metastastic (n = 8) from non-metastatic (n = 41) cases in the validation cohort.

Methods

Pre-treatment plasma samples were collected of patients with benign breast disease and breast cancer and divided over a discovery (n = 31) and validation (n = 84) cohort. From the discovery cohort miRs differentially expressed between benign and malignant cases were identified using a 2,000-miR microarray. Literature-based miRs differentiating benign from malignant disease were added. Using RT-qPCR, their expression was investigated in a validation cohort consisting of pre-treatment benign, malignant and metastatic samples. Additionally, benign and malignant cases were compared to benign and malignant cases of BRCA1-mutation carriers.

Conclusions

Plasma microRNA levels differed between patients with and without breast cancer, between benign disease from wildtype and BRCA1-mutation carriers and between breast cancer with and without metastases. Hsa-miR-652-5p was validated as a potential biomarker for breast cancer.