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.

Cohesin-dependent regulation of gene expression during differentiation is lost in cohesin-mutated myeloid malignancies

Cohesin complex disruption alters gene expression, and cohesin mutations are common in myeloid neoplasia, suggesting a critical role in hematopoiesis. Here, we explore cohesin dynamics and regulation of hematopoietic stem cell homeostasis and differentiation. Cohesin binding increases at active regulatory elements only during erythroid differentiation. Prior binding of the repressive Ets transcription factor Etv6 predicts cohesin binding at these elements and Etv6 interacts with cohesin at chromatin. Depletion of cohesin severely impairs erythroid differentiation, particularly at Etv6-prebound loci, but augments self-renewal programs. Together with corroborative findings in acute myeloid leukemia and myelodysplastic syndrome patient samples, these data suggest cohesin-mediated alleviation of Etv6 repression is required for dynamic expression at critical erythroid genes during differentiation and how this may be perturbed in myeloid malignancies.

The Use of Figurative Language to Describe Frailty in Older Adults

INTRODUCTION

Frailty affects 15% of non-institutionalized older adults in the United States, yet confusion remains in defining and, in turn, assessing frailty. Figurative language, such as metaphor, can help to explain difficult scientific concepts and to form new theories. We aimed to examine the use of figurative language to describe frailty and to identify themes in the way figurative expressions are used. Understanding how frailty is described figuratively may offer insights for developing useful communication approaches in research settings.

METHODS

We performed a comprehensive review of editorials in the scientific literature to explore figurative language used to describe frailty in older adults. We categorized themes among the figurative expressions, which may help to inform how to effectively communicate about frailty.

RESULTS

We found 24 editorials containing 32 figurative expressions. The figurative expressions conceptualized frailty in six ways: 1) a complex, multifaceted concept; 2) an important issue in health and medicine; 3) indicative of something that is failing or faulty; 4) indicative of fragility; 5) representative of vulnerable, ignored persons; and 6) an opportunity for self-awareness and reflection.

DISCUSSION

Our review highlights the heterogeneity in depictions of frailty, which is consonant with the lack of a standardized definition of frailty. We also found a novel aspect to the concept of frailty, which merits attention: frailty characterized as an opportunity for self-awareness and reflection. Figurative language, which often juxtaposes familiar with challenging, complex concepts, can offer insights on issues in frailty research and holds potential as a tool for researchers to improve communication about this important and debated medical condition.

Notch Signaling Mediates Secondary Senescence

Oncogene-induced senescence (OIS) is a tumor suppressive response to oncogene activation that can be transmitted to neighboring cells through secreted factors of the senescence-associated secretory phenotype (SASP). Currently, primary and secondary senescent cells are not considered functionally distinct endpoints. Using single-cell analysis, we observed two distinct transcriptional endpoints, a primary endpoint marked by Ras and a secondary endpoint marked by Notch activation. We find that secondary oncogene-induced senescence in vitro and in vivo requires Notch, rather than SASP alone, as previously thought. Moreover, Notch signaling weakens, but does not abolish, SASP in secondary senescence. Global transcriptomic differences, a blunted SASP response, and the induction of fibrillar collagens in secondary senescence point toward a functional diversification between secondary and primary senescence.

Clonal approaches to understanding the impact of mutations on hematologic disease development

Interrogation of hematopoietic tissue at the clonal level has a rich history spanning over 50 years, and has provided critical insights into both normal and malignant hematopoiesis. Characterization of chromosomes identified some of the first genetic links to cancer with the discovery of chromosomal translocations in association with many hematological neoplasms. The unique accessibility of hematopoietic tissue and the ability to clonally expand hematopoietic progenitors in vitro has provided fundamental insights into the cellular hierarchy of normal hematopoiesis, as well as the functional impact of driver mutations in disease. Transplantation assays in murine models have enabled cellular assessment of the functional consequences of somatic mutations in vivo. Most recently, next-generation sequencing-based assays have shown great promise in allowing multi-"omic" characterization of single cells. Here, we review how clonal approaches have advanced our understanding of disease development, focusing on the acquisition of somatic mutations, clonal selection, driver mutation cooperation, and tumor evolution.

Hydroxycarbamide Plus Aspirin Versus Aspirin Alone in Patients With Essential Thrombocythemia Age 40 to 59 Years Without High-Risk Features

PURPOSE

Cytoreductive therapy is beneficial in patients with essential thrombocythemia (ET) at high risk of thrombosis. However, its value in those lacking high-risk features remains unknown. This open-label, randomized trial compared hydroxycarbamide plus aspirin with aspirin alone in patients with ET age 40 to 59 years and without high-risk factors or extreme thrombocytosis.

PATIENTS AND METHODS

Patients were age 40 to 59 years and lacked a history of ischemia, thrombosis, embolism, hemorrhage, extreme thrombocytosis (platelet count ≥ 1,500 × 109/L), hypertension, or diabetes requiring therapy. In all, 382 patients were randomly assigned 1:1 to hydroxycarbamide plus aspirin or aspirin alone. The composite primary end point was time to arterial or venous thrombosis, serious hemorrhage, or death from vascular causes. Secondary end points were time to first arterial or venous thrombosis, first serious hemorrhage, death, incidence of transformation, and patient-reported quality of life.

RESULTS

After a median follow-up of 73 months and a total follow-up of 2,373 patient-years, there was no significant difference between the arms in the likelihood of patients reaching the primary end point (hazard ratio, 0.98; 95% CI, 0.42 to 2.25; P = 1.0). The incidence of significant vascular events was low, at 0.93 per 100 patient-years (95% CI, 0.61 to 1.41). There were also no differences in overall survival; in the composite end point of transformation to myelofibrosis, acute myeloid leukemia, or myelodysplasia; in adverse events; or in patient-reported quality of life.

CONCLUSION

In patients with ET age 40 to 59 years and lacking high-risk factors for thrombosis or extreme thrombocytosis, preemptive addition of hydroxycarbamide to aspirin did not reduce vascular events, myelofibrotic transformation, or leukemic transformation. Patients age 40 to 59 years without other clinical indications for treatment (such as previous thrombosis or hemorrhage) who have a platelet count < 1,500 × 109/L should not receive cytoreductive therapy.

Classification and Personalized Prognosis in Myeloproliferative Neoplasms

BACKGROUND

Myeloproliferative neoplasms, such as polycythemia vera, essential thrombocythemia, and myelofibrosis, are chronic hematologic cancers with varied progression rates. The genomic characterization of patients with myeloproliferative neoplasms offers the potential for personalized diagnosis, risk stratification, and treatment.

METHODS

We sequenced coding exons from 69 myeloid cancer genes in patients with myeloproliferative neoplasms, comprehensively annotating driver mutations and copy-number changes. We developed a genomic classification for myeloproliferative neoplasms and multistage prognostic models for predicting outcomes in individual patients. Classification and prognostic models were validated in an external cohort.

RESULTS

A total of 2035 patients were included in the analysis. A total of 33 genes had driver mutations in at least 5 patients, with mutations in JAK2, CALR, or MPL being the sole abnormality in 45% of the patients. The numbers of driver mutations increased with age and advanced disease. Driver mutations, germline polymorphisms, and demographic variables independently predicted whether patients received a diagnosis of essential thrombocythemia as compared with polycythemia vera or a diagnosis of chronic-phase disease as compared with myelofibrosis. We defined eight genomic subgroups that showed distinct clinical phenotypes, including blood counts, risk of leukemic transformation, and event-free survival. Integrating 63 clinical and genomic variables, we created prognostic models capable of generating personally tailored predictions of clinical outcomes in patients with chronic-phase myeloproliferative neoplasms and myelofibrosis. The predicted and observed outcomes correlated well in internal cross-validation of a training cohort and in an independent external cohort. Even within individual categories of existing prognostic schemas, our models substantially improved predictive accuracy.

CONCLUSIONS

Comprehensive genomic characterization identified distinct genetic subgroups and provided a classification of myeloproliferative neoplasms on the basis of causal biologic mechanisms. Integration of genomic data with clinical variables enabled the personalized predictions of patients' outcomes and may support the treatment of patients with myeloproliferative neoplasms. (Funded by the Wellcome Trust and others.).

Population dynamics of normal human blood inferred from somatic mutations

Haematopoietic stem cells drive blood production, but their population size and lifetime dynamics have not been quantified directly in humans. Here we identified 129,582 spontaneous, genome-wide somatic mutations in 140 single-cell-derived haematopoietic stem and progenitor colonies from a healthy 59-year-old man and applied population-genetics approaches to reconstruct clonal dynamics. Cell divisions from early embryogenesis were evident in the phylogenetic tree; all blood cells were derived from a common ancestor that preceded gastrulation. The size of the stem cell population grew steadily in early life, reaching a stable plateau by adolescence. We estimate the numbers of haematopoietic stem cells that are actively making white blood cells at any one time to be in the range of 50,000-200,000. We observed adult haematopoietic stem cell clones that generate multilineage outputs, including granulocytes and B lymphocytes. Harnessing naturally occurring mutations to report the clonal architecture of an organ enables the high-resolution reconstruction of somatic cell dynamics in humans.

Single-cell approaches identify the molecular network driving malignant hematopoietic stem cell self-renewal

Recent advances in single-cell technologies have permitted the investigation of heterogeneous cell populations at previously unattainable resolution. Here we apply such approaches to resolve the molecular mechanisms driving disease in mouse hematopoietic stem cells (HSCs), using JAK2V617F mutant myeloproliferative neoplasms (MPNs) as a model. Single-cell gene expression and functional assays identified a subset of JAK2V617F mutant HSCs that display defective self-renewal. This defect is rescued at the single HSC level by crossing JAK2V617F mice with mice lacking TET2, the most commonly comutated gene in patients with MPN. Single-cell gene expression profiling of JAK2V617F-mutant HSCs revealed a loss of specific regulator genes, some of which were restored to normal levels in single TET2/JAK2 mutant HSCs. Of these, Bmi1 and, to a lesser extent, Pbx1 and Meis1 overexpression in JAK2-mutant HSCs could drive a disease phenotype and retain durable stem cell self-renewal in functional assays. Together, these single-cell approaches refine the molecules involved in clonal expansion of MPNs and have broad implications for deconstructing the molecular network of normal and malignant stem cells.

Thrombopoietin signaling to chromatin elicits rapid and pervasive epigenome remodeling within poised chromatin architectures

Thrombopoietin (TPO) is a critical cytokine regulating hematopoietic stem cell maintenance and differentiation into the megakaryocytic lineage. However, the transcriptional and chromatin dynamics elicited by TPO signaling are poorly understood. Here, we study the immediate early transcriptional and cis-regulatory responses to TPO in hematopoietic stem/progenitor cells (HSPCs) and use this paradigm of cytokine signaling to chromatin to dissect the relationship between cis-regulatory activity and chromatin architecture. We show that TPO profoundly alters the transcriptome of HSPCs, with key hematopoietic regulators being transcriptionally repressed within 30 min of TPO. By examining cis-regulatory dynamics and chromatin architectures, we demonstrate that these changes are accompanied by rapid and extensive epigenome remodeling of cis-regulatory landscapes that is spatially coordinated within topologically associating domains (TADs). Moreover, TPO-responsive enhancers are spatially clustered and engage in preferential homotypic intra- and inter-TAD interactions that are largely refractory to TPO signaling. By further examining the link between cis-regulatory dynamics and chromatin looping, we show that rapid modulation of cis-regulatory activity is largely independent of chromatin looping dynamics. Finally, we show that, although activated and repressed cis-regulatory elements share remarkably similar DNA sequence compositions, transcription factor binding patterns accurately predict rapid cis-regulatory responses to TPO.

Mutant calreticulin knockin mice develop thrombocytosis and myelofibrosis without a stem cell self-renewal advantage

Somatic mutations in the endoplasmic reticulum chaperone calreticulin (CALR) are detected in approximately 40% of patients with essential thrombocythemia (ET) and primary myelofibrosis (PMF). Multiple different mutations have been reported, but all result in a +1-bp frameshift and generate a novel protein C terminus. In this study, we generated a conditional mouse knockin model of the most common CALR mutation, a 52-bp deletion. The mutant novel human C-terminal sequence is integrated into the otherwise intact mouse CALR gene and results in mutant CALR expression under the control of the endogenous mouse locus. CALRdel/+ mice develop a transplantable ET-like disease with marked thrombocytosis, which is associated with increased and morphologically abnormal megakaryocytes and increased numbers of phenotypically defined hematopoietic stem cells (HSCs). Homozygous CALRdel/del mice developed extreme thrombocytosis accompanied by features of MF, including leukocytosis, reduced hematocrit, splenomegaly, and increased bone marrow reticulin. CALRdel/+ HSCs were more proliferative in vitro, but neither CALRdel/+ nor CALRdel/del displayed a competitive transplantation advantage in primary or secondary recipient mice. These results demonstrate the consequences of heterozygous and homozygous CALR mutations and provide a powerful model for dissecting the pathogenesis of CALR-mutant ET and PMF.

The significance of tumour microarchitectural features in breast cancer prognosis: a digital image analysis

Background

As only a minor portion of the information present in histological sections is accessible by eye, recognition and quantification of complex patterns and relationships among constituents relies on digital image analysis. In this study, our working hypothesis was that, with the application of digital image analysis technology, visually unquantifiable breast cancer microarchitectural features can be rigorously assessed and tested as prognostic parameters for invasive breast carcinoma of no special type.

Methods

Digital image analysis was performed using public domain software (ImageJ) on tissue microarrays from a cohort of 696 patients, and validated with a commercial platform (Visiopharm). Quantified features included elements defining tumour microarchitecture, with emphasis on the extent of tumour-stroma interface. The differential prognostic impact of tumour nest microarchitecture in the four immunohistochemical surrogates for molecular classification was analysed. Prognostic parameters included axillary lymph node status, breast cancer-specific survival, and time to distant metastasis. Associations of each feature with prognostic parameters were assessed using logistic regression and Cox proportional models adjusting for age at diagnosis, grade, and tumour size.

Results

An arrangement in numerous small nests was associated with axillary lymph node involvement. The association was stronger in luminal tumours (odds ratio (OR) = 1.39, p = 0.003 for a 1-SD increase in nest number, OR = 0.75, p = 0.006 for mean nest area). Nest number was also associated with survival (hazard ratio (HR) = 1.15, p = 0.027), but total nest perimeter was the parameter most significantly associated with survival in luminal tumours (HR = 1.26, p = 0.005). In the relatively small cohort of triple-negative tumours, mean circularity showed association with time to distant metastasis (HR = 1.71, p = 0.027) and survival (HR = 1.8, p = 0.02).

Conclusions

We propose that tumour arrangement in few large nests indicates a decreased metastatic potential. By contrast, organisation in numerous small nests provides the tumour with increased metastatic potential to regional lymph nodes. An outstretched pattern in small nests bestows tumours with a tendency for decreased breast cancer-specific survival. Although further validation studies are required before the argument for routine quantification of microarchitectural features is established, our approach is consistent with the demand for cost-effective methods for triaging breast cancer patients that are more likely to benefit from chemotherapy.

Determination of complex subclonal structures of hematological malignancies by multiplexed genotyping of blood progenitor colonies

Current next-generation sequencing (NGS) technologies allow unprecedented insights into the mutational profiles of tumors. Recent studies in myeloproliferative neoplasms have further demonstrated that, not only the mutational profile, but also the order in which these mutations are acquired is relevant for our understanding of the disease. Our ability to assign mutation order from NGS data alone is, however, limited. Here, we present a strategy of highly multiplexed genotyping of burst forming unit-erythroid colonies based on NGS results to assess subclonal tumor structure. This allowed for the generation of complex clonal hierarchies and determination of order of mutation acquisition far more accurately than was possible from NGS data alone.

Multiplexing for Oxidative Bisulfite Sequencing (oxBS-seq)

DNA modifications, especially methylation, are known to play a crucial part in many regulatory processes in the cell. Recently, 5-hydroxymethylcytosine (5hmC) was discovered, a DNA modification derived as an intermediate of 5-methylcytosine (5mC) oxidation. Efforts to gain insights into function of this DNA modification are underway and several methods were recently described to assess 5hmC levels using sequencing approaches. Here we integrate adaptation based multiplexing and high-efficiency library prep into the oxidative Bisulfite Sequencing (oxBS-seq) workflow reducing the starting amount and cost per sample to identify 5hmC levels genome-wide.

Myeloproliferative neoplasms: from origins to outcomes

Substantial progress has been made in our understanding of the pathogenetic basis of myeloproliferative neoplasms. The discovery of mutations in JAK2 over a decade ago heralded a new age for patient care as a consequence of improved diagnosis and the development of therapeutic JAK inhibitors. The more recent identification of mutations in calreticulin brought with it a sense of completeness, with most patients with myeloproliferative neoplasm now having a biological basis for their excessive myeloproliferation. We are also beginning to understand the processes that lead to acquisition of somatic mutations and the factors that influence subsequent clonal expansion and emergence of disease. Extended genomic profiling has established a multitude of additional acquired mutations, particularly prevalent in myelofibrosis, where their presence carries prognostic implications. A major goal is to integrate genetic, clinical, and laboratory features to identify patients who share disease biology and clinical outcome, such that therapies, both existing and novel, can be better targeted.

Myeloproliferative neoplasms: from origins to outcomes

Substantial progress has been made in our understanding of the pathogenetic basis of myeloproliferative neoplasms. The discovery of mutations in JAK2 over a decade ago heralded a new age for patient care as a consequence of improved diagnosis and the development of therapeutic JAK inhibitors. The more recent identification of mutations in calreticulin brought with it a sense of completeness, with most patients with myeloproliferative neoplasm now having a biological basis for their excessive myeloproliferation. We are also beginning to understand the processes that lead to acquisition of somatic mutations and the factors that influence subsequent clonal expansion and emergence of disease. Extended genomic profiling has established a multitude of additional acquired mutations, particularly prevalent in myelofibrosis, where their presence carries prognostic implications. A major goal is to integrate genetic, clinical, and laboratory features to identify patients who share disease biology and clinical outcome, such that therapies, both existing and novel, can be better targeted.

Mbd3/NuRD controls lymphoid cell fate and inhibits tumorigenesis by repressing a B cell transcriptional program

Differentiation of lineage-committed cells from multipotent progenitors requires the establishment of accessible chromatin at lineage-specific transcriptional enhancers and promoters, which is mediated by pioneer transcription factors that recruit activating chromatin remodeling complexes. Here we show that the Mbd3/nucleosome remodeling and deacetylation (NuRD) chromatin remodeling complex opposes this transcriptional pioneering during B cell programming of multipotent lymphoid progenitors by restricting chromatin accessibility at B cell enhancers and promoters. Mbd3/NuRD-deficient lymphoid progenitors therefore prematurely activate a B cell transcriptional program and are biased toward overproduction of pro-B cells at the expense of T cell progenitors. The striking reduction in early thymic T cell progenitors results in compensatory hyperproliferation of immature thymocytes and development of T cell lymphoma. Our results reveal that Mbd3/NuRD can regulate multilineage differentiation by constraining the activation of dormant lineage-specific enhancers and promoters. In this way, Mbd3/NuRD protects the multipotency of lymphoid progenitors, preventing B cell-programming transcription factors from prematurely enacting lineage commitment. Mbd3/NuRD therefore controls the fate of lymphoid progenitors, ensuring appropriate production of lineage-committed progeny and suppressing tumor formation.

Pathogenesis of Myeloproliferative Disorders

Myeloproliferative neoplasms (MPNs) are a set of chronic hematopoietic neoplasms with overlapping clinical and molecular features. Recent years have witnessed considerable advances in our understanding of their pathogenetic basis. Due to their protracted clinical course, the evolution to advanced hematological malignancies, and the accessibility of neoplastic tissue, the study of MPNs has provided a window into the earliest stages of tumorigenesis. With the discovery of mutations in CALR, the majority of MPN patients now bear an identifiable marker of clonal disease; however, the mechanism by which mutated CALR perturbs megakaryopoiesis is currently unresolved. We are beginning to understand better the role of JAK2(V617F) homozygosity, the function of comutations in epigenetic regulators and spliceosome components, and how these mutations cooperate with JAK2(V617F) to modulate MPN phenotype.

Proliferation Drives Aging-Related Functional Decline in a Subpopulation of the Hematopoietic Stem Cell Compartment

Aging of the hematopoietic stem cell (HSC) compartment is characterized by lineage bias and reduced stem cell function, the molecular basis of which is largely unknown. Using single-cell transcriptomics, we identified a distinct subpopulation of old HSCs carrying a p53 signature indicative of stem cell decline alongside pro-proliferative JAK/STAT signaling. To investigate the relationship between JAK/STAT and p53 signaling, we challenged HSCs with a constitutively active form of JAK2 (V617F) and observed an expansion of the p53-positive subpopulation in old mice. Our results reveal cellular heterogeneity in the onset of HSC aging and implicate a role for JAK2V617F-driven proliferation in the p53-mediated functional decline of old HSCs.

SC3: consensus clustering of single-cell RNA-seq data

Single-cell RNA-seq enables the quantitative characterization of cell types based on global transcriptome profiles. We present single-cell consensus clustering (SC3), a user-friendly tool for unsupervised clustering, which achieves high accuracy and robustness by combining multiple clustering solutions through a consensus approach (http://bioconductor.org/packages/SC3). We demonstrate that SC3 is capable of identifying subclones from the transcriptomes of neoplastic cells collected from patients.

Hoxa9 and Meis1 Cooperatively Induce Addiction to Syk Signaling by Suppressing miR-146a in Acute Myeloid Leukemia

The transcription factor Meis1 drives myeloid leukemogenesis in the context of Hox gene overexpression but is currently considered undruggable. We therefore investigated whether myeloid progenitor cells transformed by Hoxa9 and Meis1 become addicted to targetable signaling pathways. A comprehensive (phospho)proteomic analysis revealed that Meis1 increased Syk protein expression and activity. Syk upregulation occurs through a Meis1-dependent feedback loop. By dissecting this loop, we show that Syk is a direct target of miR-146a, whose expression is indirectly regulated by Meis1 through the transcription factor PU.1. In the context of Hoxa9 overexpression, Syk signaling induces Meis1, recapitulating several leukemogenic features of Hoxa9/Meis1-driven leukemia. Finally, Syk inhibition disrupts the identified regulatory loop, prolonging survival of mice with Hoxa9/Meis1-driven leukemia.

Order Matters: The Order of Somatic Mutations Influences Cancer Evolution

Cancers evolve as a consequence of multiple somatic lesions, with competition between subclones and sequential subclonal evolution. Some driver mutations arise either early or late in the evolution of different individual tumors, suggesting that the final malignant properties of a subclone reflect the sum of mutations acquired rather than the order in which they arose. However, very little is known about the cellular consequences of altering the order in which mutations are acquired. Recent studies of human myeloproliferative neoplasms show that the order in which individual mutations are acquired has a dramatic impact on the cell biological and molecular properties of tumor-initiating cells. Differences in clinical presentation, complications, and response to targeted therapy were all observed and implicate mutation order as an important player in cancer biology. These observations represent the first demonstration that the order of mutation acquisition influences stem and progenitor cell behavior and clonal evolution in any cancer. Thus far, the impact of different mutation orders has only been studied in hematological malignancies, and analogous studies of solid cancers are now required.

Abstract P6-09-16: Identification of proliferation related derivers and their roles in precision medicine for breast cancers: A retrospective multidimensional comparative, integrated genomic, transcriptomic, and protein analysis

Backgound and Aim: The best test to guide the choice of systemic therapy for breast cancer (BC) has not yet been identified. We did this study to identify factors that drive proliferation features in BC and assess their association with clinical outcomes after systemic therapy.

Methods: We applied an artificial neural network-based integrative data mining approach to three cohorts of patients with untreated lymph node (LN)-negative BC (Wang et al; n=286, Desmedt et al; n=198 and Schmidt et al; n=200). The results were validated in four cohorts of BC patients (the Nottingham discovery cohort (n=171), Uppsala cohort (n=249), The Cancer Genome Atlas-Breast Cancer project [TCGA-BRCA; n= 970] and Molecular Taxonomy of Breast Cancer International Consortium [METABRIC cohort; n=1980]. Genes that featured prominently in our interactome map of proliferation have been chosen to take them forward to investigate their clinicopathological relevance of their gene copy number aberrations (CNAs), mRNA transcript expression, and protein expression and their associations with breast cancer-specific survival (BCSS), distant relapse-free survival (DRFS) and pathological complete response (pCR) in ten international cohorts of BC (n&gt;12000 patients).

Findings: ESR1, SPAG5, EGFR, BCL2, and FOXA1 were among the 39 common gene probes that were predictive across most proliferation features and datasets. In TCGA-BRCA cohort, SPAG5 gene mutation, gain/amplification and loss at the Ch17q11.2 locus were detected in 43 (4.4%), 177 (18.2%) and 180 (18.8%) of 970 patients, respectively and 65 (31%) of 479 ER-positive /HER-positive patients showed gain/amplification of SPAG5 gene. In multivariable analysis, high SPAG5 transcript and SPAG5 protein expression were associated with reduced BCSS compared with lower expression (METABRIC: HR 1·27, 95% CI 1·02–1·58, p=0·034; untreated LN-negative cohort: 2·34, 1·24–4·42, p=0·0090; and Nottingham-cohort: 1·73, 1·23–2·46, p=0·0020). In patients with ER-negative/HER2-negative or ER-positive/HER2-negative BC, high SPAG5 transcript expression was associated with an increased pCR compared with low SPAG5 transcript expression after receiving anthracycline neoadjuvant chemotherapy (AC-NeoACT) [(Multicentre phase 2 clinical trial cohort; n=136; OR 2·47, 95% CI 1·17–5·21, p=0.016) and (MD Anderson- taxane+AC-NeoACT cohort; n=287; OR 3·16, 95% CI 1·46–6·84, p=0.003); respectively]. In patients with ER-positive/HER2-negative BC who received taxane+AC-NeoACT followed by adjuvant tamoxifen (Adj-Tam) for 5 years (MD Anderson- taxane+AC-NeoACT cohort; n=287), high and low SPAG5 transcript expression had similar DRFS (HR 1·40, 95% CI 0.76–2·58, p=0.282). Whereas in ER-positive/HER2-negative BC patients who received only adj-Tam (n=298), high SPAG5 transcript expression was associated with reduced DRF at 5 years compared with lower expression (HR 1.98, 95% CI 1.19–3.27, p=0.008).

Interpretation: The transcript and protein products of SPAG5 are independent prognostic and predictive biomarkers that might have clinical utility as biomarkers for combination cytotoxic chemotherapy sensitivity in ER-positive/HER-negative BC.

Citation Format: Abdel-Fatah TMA, Agarwal D, Zafeiris D, Pongor L, Györffy B, Rueda OM, Moseley PM, Green AR, Liu D-X, Pockley AG, Rees RC, Caldas C, Ellis IO, Ball GR, Chan SYT. Identification of proliferation related derivers and their roles in precision medicine for breast cancers: A retrospective multidimensional comparative, integrated genomic, transcriptomic, and protein analysis [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P6-09-16.