The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

Mutational processes molding the genomes of 21 breast cancers.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

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.

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.

Signatures of mutational processes in human cancer.

All cancers are caused by somatic mutations; however, understanding of the biological processes generating these mutations is limited. The catalogue of somatic mutations from a cancer genome bears the signatures of the mutational processes that have been operative. Here we analysed 4,938,362 mutations from 7,042 cancers and extracted more than 20 distinct mutational signatures. Some are present in many cancer types, notably a signature attributed to the APOBEC family of cytidine deaminases, whereas others are confined to a single cancer class. Certain signatures are associated with age of the patient at cancer diagnosis, known mutagenic exposures or defects in DNA maintenance, but many are of cryptic origin. In addition to these genome-wide mutational signatures, hypermutation localized to small genomic regions, 'kataegis', is found in many cancer types. The results reveal the diversity of mutational processes underlying the development of cancer, with potential implications for understanding of cancer aetiology, prevention and therapy.

Frequent somatic transfer of mitochondrial DNA into the nuclear genome of human cancer cells.

Mitochondrial genomes are separated from the nuclear genome for most of the cell cycle by the nuclear double membrane, intervening cytoplasm, and the mitochondrial double membrane. Despite these physical barriers, we show that somatically acquired mitochondrial-nuclear genome fusion sequences are present in cancer cells. Most occur in conjunction with intranuclear genomic rearrangements, and the features of the fusion fragments indicate that nonhomologous end joining and/or replication-dependent DNA double-strand break repair are the dominant mechanisms involved. Remarkably, mitochondrial-nuclear genome fusions occur at a similar rate per base pair of DNA as interchromosomal nuclear rearrangements, indicating the presence of a high frequency of contact between mitochondrial and nuclear DNA in some somatic cells. Transmission of mitochondrial DNA to the nuclear genome occurs in neoplastically transformed cells, but we do not exclude the possibility that some mitochondrial-nuclear DNA fusions observed in cancer occurred years earlier in normal somatic cells.

A small-cell lung cancer genome with complex signatures of tobacco exposure.

Cancer is driven by mutation. Worldwide, tobacco smoking is the principal lifestyle exposure that causes cancer, exerting carcinogenicity through >60 chemicals that bind and mutate DNA. Using massively parallel sequencing technology, we sequenced a small-cell lung cancer cell line, NCI-H209, to explore the mutational burden associated with tobacco smoking. A total of 22,910 somatic substitutions were identified, including 134 in coding exons. Multiple mutation signatures testify to the cocktail of carcinogens in tobacco smoke and their proclivities for particular bases and surrounding sequence context. Effects of transcription-coupled repair and a second, more general, expression-linked repair pathway were evident. We identified a tandem duplication that duplicates exons 3-8 of CHD7 in frame, and another two lines carrying PVT1-CHD7 fusion genes, indicating that CHD7 may be recurrently rearranged in this disease. These findings illustrate the potential for next-generation sequencing to provide unprecedented insights into mutational processes, cellular repair pathways and gene networks associated with cancer.

Large-scale analysis of whole genome sequencing data from formalin-fixed paraffin-embedded cancer specimens demonstrates preservation of clinical utility.

Whole genome sequencing (WGS) provides comprehensive, individualised cancer genomic information. However, routine tumour biopsies are formalin-fixed and paraffin-embedded (FFPE), damaging DNA, historically limiting their use in WGS. Here we analyse FFPE cancer WGS datasets from England's 100,000 Genomes Project, comparing 578 FFPE samples with 11,014 fresh frozen (FF) samples across multiple tumour types. We use an approach that characterises rather than discards artefacts. We identify three artefactual signatures, including one known (SBS57) and two previously uncharacterised (SBS FFPE, ID FFPE), and develop an "FFPEImpact" score that quantifies sample artefacts. Despite inferior sequencing quality, FFPE-derived data identifies clinically-actionable variants, mutational signatures and permits algorithmic stratification. Matched FF/FFPE validation cohorts shows good concordance while acknowledging SBS, ID and copy-number artefacts. While FF-derived WGS data remains the gold standard, FFPE-samples can be used for WGS if required, using analytical advancements developed here, potentially democratising whole cancer genomics to many.

Diagnosis of Ovarian Carcinoma Homologous Recombination DNA Repair Deficiency From Targeted Gene Capture Oncology Assays.

PURPOSE: Homologous recombination DNA repair deficiency (HRD) is a therapeutic biomarker for sensitivity to platinum and poly(ADP-ribose) polymerase inhibitor therapies in breast and ovarian cancers. Several molecular phenotypes and diagnostic strategies have been developed to assess HRD; however, their clinical implementation remains both technically challenging and methodologically unstandardized. METHODS: We developed and validated an efficient and cost-effective strategy for HRD determination on the basis of calculation of a genome-wide loss of heterozygosity (LOH) score through targeted, hybridization capture and next-generation DNA sequencing augmented with 3,000 common, polymorphic single-nucleotide polymorphism (SNP) sites distributed genome-wide. This approach requires minimal sequence reads and can be readily integrated into targeted gene capture workflows already in use for molecular oncology. We interrogated 99 ovarian neoplasm-normal pairs using this method and compared results with patient mutational genotypes and orthologous predictors of HRD derived from whole-genome mutational signatures. RESULTS: LOH scores of ≥11% had >86% sensitivity for identifying tumors with HRD-causing mutations in an independent validation set (90.6% sensitivity for all specimens). We found strong agreement of our analytic approach with genome-wide mutational signature assays for determining HRD, yielding an estimated 96.7% sensitivity and 50% specificity. We observed poor concordance with mutational signatures inferred using only mutations detected by the targeted gene capture panel, suggesting inadequacy of the latter approach. LOH score did not significantly correlate with treatment outcomes. CONCLUSION: Targeted sequencing of genome-wide polymorphic SNP sites can be used to infer LOH events and subsequently diagnose HRD in ovarian tumors. The methods presented here are readily generalizable to other targeted gene oncology assays and could be adapted for HRD diagnosis in other tumor types.

A living biobank of patient-derived ductal carcinoma in situ mouse-intraductal xenografts identifies risk factors for invasive progression.

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer (IBC). Due to a lack of biomarkers able to distinguish high- from low-risk cases, DCIS is treated similar to early IBC even though the minority of untreated cases eventually become invasive. Here, we characterized 115 patient-derived mouse-intraductal (MIND) DCIS models reflecting the full spectrum of DCIS observed in patients. Utilizing the possibility to follow the natural progression of DCIS combined with omics and imaging data, we reveal multiple prognostic factors for high-risk DCIS including high grade, HER2 amplification, expansive 3D growth, and high burden of copy number aberrations. In addition, sequential transplantation of xenografts showed minimal phenotypic and genotypic changes over time, indicating that invasive behavior is an intrinsic phenotype of DCIS and supporting a multiclonal evolution model. Moreover, this study provides a collection of 19 distributable DCIS-MIND models spanning all molecular subtypes.

A living biobank of canine mammary tumor organoids as a comparative model for human breast cancer.

Mammary tumors in dogs hold great potential as naturally occurring breast cancer models in translational oncology, as they share the same environmental risk factors, key histological features, hormone receptor expression patterns, prognostic factors, and genetic characteristics as their human counterparts. We aimed to develop in vitro tools that allow functional analysis of canine mammary tumors (CMT), as we have a poor understanding of the underlying biology that drives the growth of these heterogeneous tumors. We established the long-term culture of 24 organoid lines from 16 dogs, including organoids derived from normal mammary epithelium or benign lesions. CMT organoids recapitulated key morphological and immunohistological features of the primary tissue from which they were derived, including hormone receptor status. Furthermore, genetic characteristics (driver gene mutations, DNA copy number variations, and single-nucleotide variants) were conserved within tumor-organoid pairs. We show how CMT organoids are a suitable model for in vitro drug assays and can be used to investigate whether specific mutations predict therapy outcomes. Specifically, certain CMT subtypes, such as PIK3CA mutated, estrogen receptor-positive simple carcinomas, can be valuable in setting up a preclinical model highly relevant to human breast cancer research. In addition, we could genetically modify the CMT organoids and use them to perform pooled CRISPR/Cas9 screening, where library representation was accurately maintained. In summary, we present a robust 3D in vitro preclinical model that can be used in translational research, where organoids from normal, benign as well as malignant mammary tissues can be propagated from the same animal to study tumorigenesis.

Substitution mutational signatures in whole-genome-sequenced cancers in the UK population.

Whole-genome sequencing (WGS) permits comprehensive cancer genome analyses, revealing mutational signatures, imprints of DNA damage and repair processes that have arisen in each patient's cancer. We performed mutational signature analyses on 12,222 WGS tumor-normal matched pairs, from patients recruited via the UK National Health Service. We contrasted our results to two independent cancer WGS datasets, the International Cancer Genome Consortium (ICGC) and Hartwig Foundation, involving 18,640 WGS cancers in total. Our analyses add 40 single and 18 double substitution signatures to the current mutational signature tally. Critically, we show for each organ, that cancers have a limited number of 'common' signatures and a long tail of 'rare' signatures. We provide a practical solution for utilizing this concept of common versus rare signatures in future analyses.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer.

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

Whole-Genome Sequencing of Retinoblastoma Reveals the Diversity of Rearrangements Disrupting RB1 and Uncovers a Treatment-Related Mutational Signature.

The development of retinoblastoma is thought to require pathological genetic changes in both alleles of the RB1 gene. However, cases exist where RB1 mutations are undetectable, suggesting alternative pathways to malignancy. We used whole-genome sequencing (WGS) and transcriptomics to investigate the landscape of sporadic retinoblastomas derived from twenty patients, sought RB1 and other driver mutations and investigated mutational signatures. At least one RB1 mutation was identified in all retinoblastomas, including new mutations in addition to those previously identified by clinical screening. Ten tumours carried structural rearrangements involving RB1 ranging from relatively simple to extremely complex rearrangement patterns, including a chromothripsis-like pattern in one tumour. Bilateral tumours obtained from one patient harboured conserved germline but divergent somatic RB1 mutations, indicating independent evolution. Mutational signature analysis showed predominance of signatures associated with cell division, an absence of ultraviolet-related DNA damage and a profound platinum-related mutational signature in a chemotherapy-exposed tumour. Most RB1 mutations are identifiable by clinical screening. However, the increased resolution and ability to detect otherwise elusive rearrangements by WGS have important repercussions on clinical management and advice on recurrence risks.

Holistic cancer genome profiling for every patient.

Technological advances in the ability to read the human genome have accelerated the speed of sequencing, such that today we can perform whole genome sequencing (WGS) in one day. Until recently, genomic studies have largely been limited to seeking novel scientific discoveries. The application of new insights gained through cancer WGS into the clinical domain, have been relatively limited. Looking ahead, a vast amount of data can be generated by genomic studies. Of note, excellent organisation of genomic and clinical data permits the application of machine-learning methods which can lead to the development of clinical algorithms that could assist future clinicians and genomicists in the analysis and interpretation of individual cancer genomes. Here, we describe what can be gleaned from holistic whole cancer genome profiling and argue that we must build the infrastructure and educational frameworks to support the modern clinical genomicist to prepare for a future where WGS will be the norm.

Epigenetic modifiers DNMT3A and BCOR are recurrently mutated in CYLD cutaneous syndrome.

Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. Here, we comprehensively profile the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families and discover recurrent mutations in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumors suggest that isoform-specific DNMT3A2 mutations are associated with dysregulated methylation. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumors. These findings contribute to existing paradigms of cutaneous tumorigenesis and metastasis.

Dramatic response of metastatic cutaneous angiosarcoma to an immune checkpoint inhibitor in a patient with xeroderma pigmentosum: whole-genome sequencing aids treatment decision in end-stage disease.

"Mutational signatures" are patterns of mutations that report DNA damage and subsequent repair processes that have occurred. Whole-genome sequencing (WGS) can provide additional information to standard diagnostic techniques and can identify therapeutic targets. A 32-yr-old male with xeroderma pigmentosum developed metastatic angiosarcoma that was unresponsive to three lines of conventional sarcoma therapies. WGS was performed on his primary cancer revealing a hypermutated tumor, including clonal ultraviolet radiation-induced mutational patterns (Signature 7) and subclonal signatures of mutated DNA polymerase epsilon (POLE) (Signature 10). These signatures are associated with response to immune checkpoint blockade. Immunohistochemistry confirmed high PD-L1 expression in metastatic deposits. The anti-PD-1 monoclonal antibody pembrolizumab was commenced off-label given the POLE mutation and high mutational load. After four cycles, there was a significant reduction in his disease with almost complete resolution of the metastatic deposits. This case highlights the importance of WGS in the analysis, interpretation, and treatment of cancers. We anticipate that as WGS becomes integral to the cancer diagnostic pathway, treatments will be stratified to the individual based on their unique genomic and/or transcriptomic profile, enhancing classical approaches of histologically driven treatment decisions.

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

Whole-genome sequencing (WGS) brings comprehensive insights to cancer genome interpretation. To explore the clinical value of WGS, we sequenced 254 triple-negative breast cancers (TNBCs) for which associated treatment and outcome data were collected between 2010 and 2015 via the population-based Sweden Cancerome Analysis Network-Breast (SCAN-B) project (ClinicalTrials.gov ID:NCT02306096). Applying the HRDetect mutational-signature-based algorithm to classify tumors, 59% were predicted to have homologous-recombination-repair deficiency (HRDetect-high): 67% explained by germline/somatic mutations of BRCA1/BRCA2, BRCA1 promoter hypermethylation, RAD51C hypermethylation or biallelic loss of PALB2. A novel mechanism of BRCA1 abrogation was discovered via germline SINE-VNTR-Alu retrotransposition. HRDetect provided independent prognostic information, with HRDetect-high patients having better outcome on adjuvant chemotherapy for invasive disease-free survival (hazard ratio (HR) = 0.42; 95% confidence interval (CI) = 0.2-0.87) and distant relapse-free interval (HR = 0.31, CI = 0.13-0.76) compared to HRDetect-low, regardless of whether a genetic/epigenetic cause was identified. HRDetect-intermediate, some possessing potentially targetable biological abnormalities, had the poorest outcomes. HRDetect-low cancers also had inadequate outcomes: ~4.7% were mismatch-repair-deficient (another targetable defect, not typically sought) and they were enriched for (but not restricted to) PIK3CA/AKT1 pathway abnormalities. New treatment options need to be considered for now-discernible HRDetect-intermediate and HRDetect-low categories. This population-based study advocates for WGS of TNBC to better inform trial stratification and improve clinical decision-making.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer.

Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.