The impact of AI on research.

Large and complex datasets have made artificial intelligence (AI) an invaluable tool for discovery across biological research. We asked experts how AI has impacted their work. Their experiences and perspectives offer thoughtful insights into potential offered by AI for their fields.

Predictability of B cell clonal persistence and immunosurveillance in breast cancer.

B cells and T cells are important components of the adaptive immune system and mediate anticancer immunity. The T cell landscape in cancer is well characterized, but the contribution of B cells to anticancer immunosurveillance is less well explored. Here we show an integrative analysis of the B cell and T cell receptor repertoire from individuals with metastatic breast cancer and individuals with early breast cancer during neoadjuvant therapy. Using immune receptor, RNA and whole-exome sequencing, we show that both B cell and T cell responses seem to coevolve with the metastatic cancer genomes and mirror tumor mutational and neoantigen architecture. B cell clones associated with metastatic immunosurveillance and temporal persistence were more expanded and distinct from site-specific clones. B cell clonal immunosurveillance and temporal persistence are predictable from the clonal structure, with higher-centrality B cell antigen receptors more likely to be detected across multiple metastases or across time. This predictability was generalizable across other immune-mediated disorders. This work lays a foundation for prioritizing antibody sequences for therapeutic targeting in cancer.

Cancer Treatment in the Genomic Era.

The complexity of human cancer underlies its devastating clinical consequences. Drugs designed to target the genetic alterations that drive cancer have improved the outcome for many patients, but not the majority of them. Here, we review the genomic landscape of cancer, how genomic data can provide much more than a sum of its parts, and the approaches developed to identify and validate genomic alterations with potential therapeutic value. We highlight notable successes and pitfalls in predicting the value of potential therapeutic targets and discuss the use of multi-omic data to better understand cancer dependencies and drug sensitivity. We discuss how integrated approaches to collecting, curating, and sharing these large data sets might improve the identification and prioritization of cancer vulnerabilities as well as patient stratification within clinical trials. Finally, we outline how future approaches might improve the efficiency and speed of translating genomic data into clinically effective therapies and how the use of unbiased genome-wide information can identify novel predictive biomarkers that can be either simple or complex.

Promoter of lncRNA Gene PVT1 Is a Tumor-Suppressor DNA Boundary Element.

Noncoding mutations in cancer genomes are frequent but challenging to interpret. PVT1 encodes an oncogenic lncRNA, but recurrent translocations and deletions in human cancers suggest alternative mechanisms. Here, we show that the PVT1 promoter has a tumor-suppressor function that is independent of PVT1 lncRNA. CRISPR interference of PVT1 promoter enhances breast cancer cell competition and growth in vivo. The promoters of the PVT1 and the MYC oncogenes, located 55 kb apart on chromosome 8q24, compete for engagement with four intragenic enhancers in the PVT1 locus, thereby allowing the PVT1 promoter to regulate pause release of MYC transcription. PVT1 undergoes developmentally regulated monoallelic expression, and the PVT1 promoter inhibits MYC expression only from the same chromosome via promoter competition. Cancer genome sequencing identifies recurrent mutations encompassing the human PVT1 promoter, and genome editing verified that PVT1 promoter mutation promotes cancer cell growth. These results highlight regulatory sequences of lncRNA genes as potential disease-associated DNA elements.

A Biobank of Breast Cancer Explants with Preserved Intra-tumor Heterogeneity to Screen Anticancer Compounds.

The inter- and intra-tumor heterogeneity of breast cancer needs to be adequately captured in pre-clinical models. We have created a large collection of breast cancer patient-derived tumor xenografts (PDTXs), in which the morphological and molecular characteristics of the originating tumor are preserved through passaging in the mouse. An integrated platform combining in vivo maintenance of these PDTXs along with short-term cultures of PDTX-derived tumor cells (PDTCs) was optimized. Remarkably, the intra-tumor genomic clonal architecture present in the originating breast cancers was mostly preserved upon serial passaging in xenografts and in short-term cultured PDTCs. We assessed drug responses in PDTCs on a high-throughput platform and validated several ex vivo responses in vivo. The biobank represents a powerful resource for pre-clinical breast cancer pharmacogenomic studies (http://caldaslab.cruk.cam.ac.uk/bcape), including identification of biomarkers of response or resistance.

Multi-omic machine learning predictor of breast cancer therapy response.

Breast cancers are complex ecosystems of malignant cells and the tumour microenvironment1. The composition of these tumour ecosystems and interactions within them contribute to responses to cytotoxic therapy2. Efforts to build response predictors have not incorporated this knowledge. We collected clinical, digital pathology, genomic and transcriptomic profiles of pre-treatment biopsies of breast tumours from 168 patients treated with chemotherapy with or without HER2 (encoded by ERBB2)-targeted therapy before surgery. Pathology end points (complete response or residual disease) at surgery3 were then correlated with multi-omic features in these diagnostic biopsies. Here we show that response to treatment is modulated by the pre-treated tumour ecosystem, and its multi-omics landscape can be integrated in predictive models using machine learning. The degree of residual disease following therapy is monotonically associated with pre-therapy features, including tumour mutational and copy number landscapes, tumour proliferation, immune infiltration and T cell dysfunction and exclusion. Combining these features into a multi-omic machine learning model predicted a pathological complete response in an external validation cohort (75 patients) with an area under the curve of 0.87. In conclusion, response to therapy is determined by the baseline characteristics of the totality of the tumour ecosystem captured through data integration and machine learning. This approach could be used to develop predictors for other cancers.

Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups.

The rates and routes of lethal systemic spread in breast cancer are poorly understood owing to a lack of molecularly characterized patient cohorts with long-term, detailed follow-up data. Long-term follow-up is especially important for those with oestrogen-receptor (ER)-positive breast cancers, which can recur up to two decades after initial diagnosis1-6. It is therefore essential to identify patients who have a high risk of late relapse7-9. Here we present a statistical framework that models distinct disease stages (locoregional recurrence, distant recurrence, breast-cancer-related death and death from other causes) and competing risks of mortality from breast cancer, while yielding individual risk-of-recurrence predictions. We apply this model to 3,240 patients with breast cancer, including 1,980 for whom molecular data are available, and delineate spatiotemporal patterns of relapse across different categories of molecular information (namely immunohistochemical subtypes; PAM50 subtypes, which are based on gene-expression patterns10,11; and integrative or IntClust subtypes, which are based on patterns of genomic copy-number alterations and gene expression12,13). We identify four late-recurring integrative subtypes, comprising about one quarter (26%) of tumours that are both positive for ER and negative for human epidermal growth factor receptor 2, each with characteristic tumour-driving alterations in genomic copy number and a high risk of recurrence (mean 47-62%) up to 20 years after diagnosis. We also define a subgroup of triple-negative breast cancers in which cancer rarely recurs after five years, and a separate subgroup in which patients remain at risk. Use of the integrative subtypes improves the prediction of late, distant relapse beyond what is possible with clinical covariates (nodal status, tumour size, tumour grade and immunohistochemical subtype). These findings highlight opportunities for improved patient stratification and biomarker-driven clinical trials.

Clustering of HR + /HER2- breast cancer in an Asian cohort is driven by immune phenotypes.

Breast cancer exhibits significant heterogeneity, manifesting in various subtypes that are critical in guiding treatment decisions. This study aimed to investigate the existence of distinct subtypes of breast cancer within the Asian population, by analysing the transcriptomic profiles of 934 breast cancer patients from a Malaysian cohort. Our findings reveal that the HR + /HER2- breast cancer samples display a distinct clustering pattern based on immune phenotypes, rather than conforming to the conventional luminal A-luminal B paradigm previously reported in breast cancers from women of European descent. This suggests that the activation of the immune system may play a more important role in Asian HR + /HER2- breast cancer than has been previously recognized. Analysis of somatic mutations by whole exome sequencing showed that counter-intuitively, the cluster of HR + /HER2- samples exhibiting higher immune scores was associated with lower tumour mutational burden, lower homologous recombination deficiency scores, and fewer copy number aberrations, implicating the involvement of non-canonical tumour immune pathways. Further investigations are warranted to determine the underlying mechanisms of these pathways, with the potential to develop innovative immunotherapeutic approaches tailored to this specific patient population.

Positive correlation between transcriptomic stemness and PI3K/AKT/mTOR signaling scores in breast cancer, and a counterintuitive relationship with PIK3CA genotype.

A PI3Kα-selective inhibitor has recently been approved for use in breast tumors harboring mutations in PIK3CA, the gene encoding p110α. Preclinical studies have suggested that the PI3K/AKT/mTOR signaling pathway influences stemness, a dedifferentiation-related cellular phenotype associated with aggressive cancer. However, to date, no direct evidence for such a correlation has been demonstrated in human tumors. In two independent human breast cancer cohorts, encompassing nearly 3,000 tumor samples, transcriptional footprint-based analysis uncovered a positive linear association between transcriptionally-inferred PI3K/AKT/mTOR signaling scores and stemness scores. Unexpectedly, stratification of tumors according to PIK3CA genotype revealed a "biphasic" relationship of mutant PIK3CA allele dosage with these scores. Relative to tumor samples without PIK3CA mutations, the presence of a single copy of a hotspot PIK3CA variant was associated with lower PI3K/AKT/mTOR signaling and stemness scores, whereas the presence of multiple copies of PIK3CA hotspot mutations correlated with higher PI3K/AKT/mTOR signaling and stemness scores. This observation was recapitulated in a human cell model of heterozygous and homozygous PIK3CAH1047R expression. Collectively, our analysis (1) provides evidence for a signaling strength-dependent PI3K-stemness relationship in human breast cancer; (2) supports evaluation of the potential benefit of patient stratification based on a combination of conventional PI3K pathway genetic information with transcriptomic indices of PI3K signaling activation.

Clinical interpretation of whole-genome and whole-transcriptome sequencing for precision oncology.

Whole-genome sequencing either alone or in combination with whole-transcriptome sequencing has started to be used to analyze clinical tumor samples to improve diagnosis, provide risk stratification, and select patient-specific therapies. Compared with current genomic testing strategies, largely focused on small number of genes tested individually or targeted panels, whole-genome and transcriptome sequencing (WGTS) provides novel opportunities to identify and report a potentially much larger number of actionable alterations with diagnostic, prognostic, and/or predictive impact. Such alterations include point mutations, indels, copy- number aberrations and structural variants, but also germline variants, fusion genes, noncoding alterations and mutational signatures. Nevertheless, these comprehensive tests are accompanied by many challenges ranging from the extent and diversity of sequence alterations detected by these methods to the complexity and limited existing standardization in interpreting them. We describe the challenges of WGTS interpretation and the opportunities with comprehensive genomic testing.

Clinical utility of whole-genome sequencing in precision oncology.

Precision diagnostics is one of the two pillars of precision medicine. Sequencing efforts in the past decade have firmly established cancer as a primarily genetically driven disease. This concept is supported by therapeutic successes aimed at particular pathways that are perturbed by specific driver mutations in protein-coding domains and reflected in three recent FDA tissue agnostic cancer drug approvals. In addition, there is increasing evidence from studies that interrogate the entire genome by whole-genome sequencing that acquired global and complex genomic aberrations including those in non-coding regions of the genome might also reflect clinical outcome. After addressing technical, logistical, financial and ethical challenges, national initiatives now aim to introduce clinical whole-genome sequencing into real-world diagnostics as a rational and potentially cost-effective tool for response prediction in cancer and to identify patients who would benefit most from 'expensive' targeted therapies and recruitment into clinical trials. However, so far, this has not been accompanied by a systematic and prospective evaluation of the clinical utility of whole-genome sequencing within clinical trials of uniformly treated patients of defined clinical outcome. This approach would also greatly facilitate novel predictive biomarker discovery and validation, ultimately reducing size and duration of clinical trials and cost of drug development. This manuscript is the third in a series of three to review and critically appraise the potential and challenges of clinical whole-genome sequencing in solid tumors and hematological malignancies.

Analytical demands to use whole-genome sequencing in precision oncology.

Interrogating the tumor genome in its entirety by whole-genome sequencing (WGS) offers an unprecedented insight into the biology and pathogenesis of cancer, with potential impact on diagnostics, prognostication and therapy selection. WGS is able to detect sequence as well as structural variants and thereby combines central domains of cytogenetics and molecular genetics. Given the potential of WGS in directing targeted therapeutics and clinical decision-making, we envision a gradual transition of the method from research to clinical routine. This review is one out of three within this issue aimed at facilitating this effort, by discussing in-depth analytical validation, clinical interpretation and clinical utility of WGS. The review highlights the requirements for implementing, validating and maintaining a clinical WGS pipeline to obtain high-quality patient-specific data in accordance with the local regulatory landscape. Every step of the WGS pipeline, which includes DNA extraction, library preparation, sequencing, bioinformatics analysis, and data storage, is considered with respect to its logistics, necessities, potential pitfalls, and the required quality management. WGS is likely to drive clinical diagnostics and patient care forward, if requirements and challenges of the technique are recognized and met.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer.

BACKGROUND: Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance. METHODS: Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis. RESULTS: We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care. CONCLUSION: In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

Somatic mutations reveal asymmetric cellular dynamics in the early human embryo.

Somatic cells acquire mutations throughout the course of an individual's life. Mutations occurring early in embryogenesis are often present in a substantial proportion of, but not all, cells in postnatal humans and thus have particular characteristics and effects. Depending on their location in the genome and the proportion of cells they are present in, these mosaic mutations can cause a wide range of genetic disease syndromes and predispose carriers to cancer. They have a high chance of being transmitted to offspring as de novo germline mutations and, in principle, can provide insights into early human embryonic cell lineages and their contributions to adult tissues. Although it is known that gross chromosomal abnormalities are remarkably common in early human embryos, our understanding of early embryonic somatic mutations is very limited. Here we use whole-genome sequences of normal blood from 241 adults to identify 163 early embryonic mutations. We estimate that approximately three base substitution mutations occur per cell per cell-doubling event in early human embryogenesis and these are mainly attributable to two known mutational signatures. We used the mutations to reconstruct developmental lineages of adult cells and demonstrate that the two daughter cells of many early embryonic cell-doubling events contribute asymmetrically to adult blood at an approximately 2:1 ratio. This study therefore provides insights into the mutation rates, mutational processes and developmental outcomes of cell dynamics that operate during early human embryogenesis.

Imaging breast cancer using hyperpolarized carbon-13 MRI.

Our purpose is to investigate the feasibility of imaging tumor metabolism in breast cancer patients using 13C magnetic resonance spectroscopic imaging (MRSI) of hyperpolarized 13C label exchange between injected [1-13C]pyruvate and the endogenous tumor lactate pool. Treatment-naïve breast cancer patients were recruited: four triple-negative grade 3 cancers; two invasive ductal carcinomas that were estrogen and progesterone receptor-positive (ER/PR+) and HER2/neu-negative (HER2-), one grade 2 and one grade 3; and one grade 2 ER/PR+ HER2- invasive lobular carcinoma (ILC). Dynamic 13C MRSI was performed following injection of hyperpolarized [1-13C]pyruvate. Expression of lactate dehydrogenase A (LDHA), which catalyzes 13C label exchange between pyruvate and lactate, hypoxia-inducible factor-1 (HIF1α), and the monocarboxylate transporters MCT1 and MCT4 were quantified using immunohistochemistry and RNA sequencing. We have demonstrated the feasibility and safety of hyperpolarized 13C MRI in early breast cancer. Both intertumoral and intratumoral heterogeneity of the hyperpolarized pyruvate and lactate signals were observed. The lactate-to-pyruvate signal ratio (LAC/PYR) ranged from 0.021 to 0.473 across the tumor subtypes (mean ± SD: 0.145 ± 0.164), and a lactate signal was observed in all of the grade 3 tumors. The LAC/PYR was significantly correlated with tumor volume (R = 0.903, P = 0.005) and MCT 1 (R = 0.85, P = 0.032) and HIF1α expression (R = 0.83, P = 0.043). Imaging of hyperpolarized [1-13C]pyruvate metabolism in breast cancer is feasible and demonstrated significant intertumoral and intratumoral metabolic heterogeneity, where lactate labeling correlated with MCT1 expression and hypoxia.

Residual cancer burden after neoadjuvant chemotherapy and long-term survival outcomes in breast cancer: a multicentre pooled analysis of 5161 patients.

BACKGROUND: Previous studies have independently validated the prognostic relevance of residual cancer burden (RCB) after neoadjuvant chemotherapy. We used results from several independent cohorts in a pooled patient-level analysis to evaluate the relationship of RCB with long-term prognosis across different phenotypic subtypes of breast cancer, to assess generalisability in a broad range of practice settings. METHODS: In this pooled analysis, 12 institutes and trials in Europe and the USA were identified by personal communications with site investigators. We obtained participant-level RCB results, and data on clinical and pathological stage, tumour subtype and grade, and treatment and follow-up in November, 2019, from patients (aged ≥18 years) with primary stage I-III breast cancer treated with neoadjuvant chemotherapy followed by surgery. We assessed the association between the continuous RCB score and the primary study outcome, event-free survival, using mixed-effects Cox models with the incorporation of random RCB and cohort effects to account for between-study heterogeneity, and stratification to account for differences in baseline hazard across cancer subtypes defined by hormone receptor status and HER2 status. The association was further evaluated within each breast cancer subtype in multivariable analyses incorporating random RCB and cohort effects and adjustments for age and pretreatment clinical T category, nodal status, and tumour grade. Kaplan-Meier estimates of event-free survival at 3, 5, and 10 years were computed for each RCB class within each subtype. FINDINGS: We analysed participant-level data from 5161 patients treated with neoadjuvant chemotherapy between Sept 12, 1994, and Feb 11, 2019. Median age was 49 years (IQR 20-80). 1164 event-free survival events occurred during follow-up (median follow-up 56 months [IQR 0-186]). RCB score was prognostic within each breast cancer subtype, with higher RCB score significantly associated with worse event-free survival. The univariable hazard ratio (HR) associated with one unit increase in RCB ranged from 1·55 (95% CI 1·41-1·71) for hormone receptor-positive, HER2-negative patients to 2·16 (1·79-2·61) for the hormone receptor-negative, HER2-positive group (with or without HER2-targeted therapy; p<0·0001 for all subtypes). RCB score remained prognostic for event-free survival in multivariable models adjusted for age, grade, T category, and nodal status at baseline: the adjusted HR ranged from 1·52 (1·36-1·69) in the hormone receptor-positive, HER2-negative group to 2·09 (1·73-2·53) in the hormone receptor-negative, HER2-positive group (p<0·0001 for all subtypes). INTERPRETATION: RCB score and class were independently prognostic in all subtypes of breast cancer, and generalisable to multiple practice settings. Although variability in hormone receptor subtype definitions and treatment across patients are likely to affect prognostic performance, the association we observed between RCB and a patient's residual risk suggests that prospective evaluation of RCB could be considered to become part of standard pathology reporting after neoadjuvant therapy. FUNDING: National Cancer Institute at the US National Institutes of Health.

Androgen receptor function and targeted therapeutics across breast cancer subtypes.

Despite significant progress in breast cancer (BC) therapy, it is globally the most commonly diagnosed cancer and leads to the death of over 650,000 women annually. Androgen receptor (AR) is emerging as a potential new therapeutic target in BC. While the role of AR is well established in prostate cancer (PCa), its function in BC remains incompletely understood. Emerging data show that AR's role in BC is dependent on several factors including, but not limited to, disease subtype, tumour microenvironment, and levels of circulating oestrogens and androgens. While targeting AR in PCa is becoming increasingly effective, these advances have yet to make any significant impact on the care of BC patients. However, this approach is increasingly being evaluated in BC and it is clear that improvements in our understanding of AR's role in BC will increase the likelihood of success for AR-targeted therapies. This review summarizes our current understanding of the function of AR across BC subtypes. We highlight limitations in our current knowledge and demonstrate the importance of categorizing BC subtypes effectively, in relation to determining AR activity. Further, we describe the current state of the art regarding AR-targeted approaches for BC as monotherapy or in combination with radiotherapy.

Copy number heterogeneity identifies ER+ breast cancer patients that do not benefit from adjuvant endocrine therapy.

BACKGROUND: Endocrine therapy forms the backbone of adjuvant treatment for oestrogen-receptor-positive (ER+) breast cancer. However, it remains unclear whether adjuvant treatment improves survival rates in low-risk patients. Low intra-tumour heterogeneity (ITH) has been shown to confer low risk for recurrent disease. Here, it is studied if chromosomal copy-number ITH (CNH) can identify low-risk ER+, lymph-node-negative breast cancer patients who do not benefit from adjuvant endocrine therapy. METHODS: Lymph-node-negative ER+ patients from the observational METABRIC dataset were retrospectively analysed (n = 708). CNH was determined from a single bulk copy-number measurement for each patient. Survival rates were compared between patients that did or did not receive adjuvant endocrine therapy for CNH-low, middle and high groups with Cox proportional-hazards models, using propensity-score weights to correct for confounders. RESULTS: Adjuvant endocrine therapy improved the relapse-free survival (RFS) for CNH-high patients treatment (HR = 0.55), but not for CNH-low patients treatment (HR = 0.88). For CNH-low patients adjuvant endocrine therapy was associated with impaired OS (HR = 1.62). CONCLUSIONS: This retrospective study of lymph-node-negative, ER+ breast cancer finds that patients identified as low risk using CNH do not benefit from adjuvant endocrine therapy.

Modeling the Prognostic Impact of Circulating Tumor Cells Enumeration in Metastatic Breast Cancer for Clinical Trial Design Simulation.

Despite the strong prognostic stratification of circulating tumor cells (CTCs) enumeration in metastatic breast cancer (MBC), current clinical trials usually do not include a baseline CTCs in their design. This study aimed to generate a classifier for CTCs prognostic simulation in existing datasets for hypothesis generation in patients with MBC. A K-nearest neighbor machine learning algorithm was trained on a pooled dataset comprising 2436 individual MBC patients from the European Pooled Analysis Consortium and the MD Anderson Cancer Center to identify patients likely to have CTCs ≥ 5/7 mL blood (StageIVaggressive vs StageIVindolent). The model had a 65.1% accuracy and its prognostic impact resulted in a hazard ratio (HR) of 1.89 (Simulatedaggressive vs SimulatedindolentP < .001), similar to patients with actual CTCs enumeration (HR 2.76; P < .001). The classifier's performance was then tested on an independent retrospective database comprising 446 consecutive hormone receptor (HR)-positive HER2-negative MBC patients. The model further stratified clinical subgroups usually considered prognostically homogeneous such as patients with bone-only or liver metastases. Bone-only disease classified as Simulatedaggressive had a significantly worse overall survival (OS; P < .0001), while patients with liver metastases classified as Simulatedindolent had a significantly better prognosis (P < .0001). Consistent results were observed for patients who had undergone CTCs enumeration in the pooled population. The differential prognostic impact of endocrine- (ET) and chemotherapy (CT) was explored across the simulated subgroups. No significant differences were observed between ET and CT in the overall population, both in terms of progression-free survival (PFS) and OS. In contrast, a statistically significant difference, favoring CT over ET was observed among Simulatedaggressive patients (HR: 0.62; P = .030 and HR: 0.60; P = .037, respectively, for PFS and OS).

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.