A pan-cancer compendium of chromosomal instability.

Chromosomal instability (CIN) results in the accumulation of large-scale losses, gains and rearrangements of DNA1. The broad genomic complexity caused by CIN is a hallmark of cancer2; however, there is no systematic framework to measure different types of CIN and their effect on clinical phenotypes pan-cancer. Here we evaluate the extent, diversity and origin of CIN across 7,880 tumours representing 33 cancer types. We present a compendium of 17 copy number signatures that characterize specific types of CIN, with putative aetiologies supported by multiple independent data sources. The signatures predict drug response and identify new drug targets. Our framework refines the understanding of impaired homologous recombination, which is one of the most therapeutically targetable types of CIN. Our results illuminate a fundamental structure underlying genomic complexity in human cancers and provide a resource to guide future CIN research.

Signatures of copy number alterations in human cancer.

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Genomic profile analysis of leiomyomas with bizarre nuclei and fumarate hydratase deficient leiomyomas: Strengths, weaknesses, and limitations of array-CGH interpretation.

A close relationship has been demonstrated between genomic complexity and clinical outcome in uterine smooth muscle tumors. We studied the genomic profiles by array-CGH of 28 fumarate hydratase deficient leiomyomas and 37 leiomyomas with bizarre nuclei (LMBN) from 64 patients. Follow-up was available for 46 patients (from three to 249 months, mean 87.3 months). All patients were alive without evidence of disease. For 51 array-CGH interpretable tumors the mean Genomic Index (GI) was 16.4 (median: 9.8; from 1 to 57.8), significantly lower than the mean GI in LMS (mean GI 51.8, p < 0.001). We described three groups: (1) a group with FH deletion (24/58) with low GI (mean GI: 11 vs. 22,4, p = 0.02), (2) a group with TP53 deletion (17/58) with higher GI (22.4 vs. 11 p = 0.02), and (3) a group without genomic events on FH or TP53 genes (17/58) (mean GI:18.3; from 1 to 57.8). Because none of these tumors recurred and none showed morphological features of LMS we concluded that GI at the cut-off of 10 was not applicable in these subtypes of LM. By integration of all those findings, a GI <10 in LMBN remains a valuable argument for benignity. Conversely, in LMBN a GI >10 or alteration in tumor suppressor genes, should not alone warrant a diagnosis of malignancy. Nine tumors were tested with Nanocind CINSARC® signature and all were classified in low risk of recurrence. We propose, based on our observations, a diagnostic approach of these challenging lesions.

Gastrointestinal stromal tumor enhancers support a transcription factor network predictive of clinical outcome.

Activating mutations in the KIT or PDGFRA receptor tyrosine kinases are hallmarks of gastrointestinal stromal tumor (GIST). The biological underpinnings of recurrence following resection or disease progression beyond kinase mutation are poorly understood. Utilizing chromatin immunoprecipitation with sequencing of tumor samples and cell lines, we describe the enhancer landscape of GIST, highlighting genes that reinforce and extend our understanding of these neoplasms. A group of core transcription factors can be distinguished from others unique to localized and metastatic disease. The transcription factor HAND1 emerges in metastatic disease, binds to established GIST-associated enhancers, and facilitates GIST cell proliferation and KIT gene expression. The pattern of transcription factor expression in primary tumors is predictive of metastasis-free survival in GIST patients. These results provide insight into the enhancer landscape and transcription factor network underlying GIST, and define a unique strategy for predicting clinical behavior of this disease.

LIX1 regulates YAP activity and controls gastrointestinal cancer cell plasticity.

Gastrointestinal stromal tumours (GISTs), the most common mesenchymal neoplasm of the gastrointestinal tract, result from deregulated proliferation of transformed KIT-positive interstitial cells of Cajal that share mesenchymal progenitors with smooth muscle cells. Despite the identification of selective KIT inhibitors, primary resistance and relapse remain a major concern. Moreover, most patients develop resistance partly through reactivation of KIT and its downstream signalling pathways. We previously identified the Limb Expression 1 (LIX1) gene as a unique marker of digestive mesenchyme immaturity. We also demonstrated that LIX1 regulates mesenchymal progenitor proliferation and differentiation by controlling the Hippo effector YAP1, which is constitutively activated in many sarcomas. Therefore, we wanted to determine LIX1 role in GIST development. We found that LIX1 is strongly up-regulated in GIST samples and this is associated with unfavourable prognosis. Moreover, LIX1 controls GIST cell proliferation in vitro and in vivo. Upon LIX1 inactivation in GIST cells, YAP1/TAZ activity is reduced, KIT (the GIST signature) is down-regulated, and cells acquire smooth muscle lineage features. Our data highlight LIX1 role in digestive mesenchyme-derived cell-fate decisions and identify this novel regulator as a target for drug design for GIST treatment by influencing its differentiation status.

Impact of preoperative treatment on the CINSARC prognostic signature: translational research results from a phase 1 trial of the German Interdisciplinary Sarcoma Group (GISG 03).

PURPOSE: CINSARC (Complexity INdex in SARComas) is a prognostic signature for soft tissue sarcoma that determines the risk for recurrence and may serve to guide the decision for adjuvant chemotherapy. The aim of this study was to compare the CINSARC signature of pre- and posttreatment biopsies of sarcoma patients treated within a phase I trial evaluating preoperative sunitinib and irradiation. METHODS: We retrieved 14 pairs of formalin-fixed paraffin-embedded blocks from pretreatment biopsies and posttreatment resection specimens and performed expression profiling of the 67 CINSARC signature genes. RESULTS: In 5/14 patients, both probes were unsuitable for expression analysis because there was no (vital) tissue left in biopsies or resection specimens. Comparing the CINSARC risk classification before and after treatment in the remaining patients, 2/9 shifted from a high- to a low-risk classification for metastatic disease after preoperative treatment with radiation therapy plus sunitinib and 7/9 pairs of pre- and posttreatment biopsies revealed identical results. CONCLUSION: Concurrent radiation therapy and sunitinib leads to diverging results of prognostic gene array testing in a relevant proportion of sarcoma patients. These changes may reflect tumor heterogeneity, local treatment effects, or prognostic changes of the disease. Caution is advised in the selection of samples and interpretation of test results.

Value of peri-operative chemotherapy in patients with CINSARC high-risk localized grade 1 or 2 soft tissue sarcoma: study protocol of the target selection phase III CHIC-STS trial.

BACKGROUND: The value of chemotherapy in soft tissue sarcoma (STS) remains controversial. Several expert teams consider that chemotherapy provides a survival advantage and should be proposed in high-risk (HR) patients. However, the lack of accuracy in identifying HR patients with conventional risk factors (large, deep, FNCLCC grade 3, extremity STS) is an issue that cannot be neglected. For example, while the FNCLCC grading system is a powerful tool, it has several limitations. CINSARC, a 67-gene signature, has proved to be an additional independent factor for predicting metastatic spread and outperforms histological grade. Regardless of FNCLCC grade, CINSARC stratifies patients into two separate prognostic groups: one with an excellent prognosis (low-risk (LR) CINSARC) and the other with a worse outcome (HR-CINSARC) in terms of metastatic relapse. Here we evaluate the role of chemotherapy in grade 1-2 STS patients with HR-CINSARC and assess the prognostic value of CINSARC in patients treated with standard of care. METHODS: CHIC is a parallel, randomized, open-label, multicenter study evaluating the effect on metastasis-free survival of adding perioperative chemotherapy to standard of care in patients with grade ½ STS sarcoma defined as HR by CINSARC. In this target selection design, 600 patients will be screened with CINSARC to randomize 250 HR-CINSARC patients between standard of care and standard of care plus chemotherapy (4 cycles of 3 weeks of intravenous chemotherapy with doxorubicin in combination with dacarbazine or ifosfamide according to histologic subtype). LR-CINSARC patients will be treated by standard of care according to the investigator. The primary endpoint is metastasis-free survival. Secondary endpoints include overall survival, disease-free survival and safety. Furthermore, the prognostic value of CINSARC will be evaluated by comparing LR-CINSARC patients to HR-CINSARC patients randomized in standard of care. DISCUSSION: CHIC is a prospective randomized phase III trial designed to comprehensively evaluate the benefit of chemotherapy in HR-CINSARC patients and to prospectively validate the prognostic value of CINSARC in grade ½ STS sarcoma patients. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT04307277 Date of registration: 13 March 2020.

CINSARC signature as a prognostic marker for clinical outcome in sarcomas and beyond.

Prognostication is a key issue for sarcoma patients' care as it triggers the therapeutic approach including chemotherapy, which is still not standard for localized patients. Current prognostic evaluation, based on the FNCLCC grading system, has recently been improved by the CINSARC signature outperforming histology-based grading system by identifying high-risk patients in every grade, even in those considered as low. CINSARC is an expression-based signature related to mitosis and chromosome integrity with prognostic value in a wide range of cancers additional to sarcoma. First developed with frozen material, CINSARC is now coupled with NanoString technology allowing evaluation from FFPE blocks used in clinical practice. Consequently, CINSARC is currently evaluated in clinical trials with a dual objective of demonstrating the benefit of chemotherapy in sarcoma patients and testing its response prediction. Considering its overarching value in oncology, its development is welcome in any cancers where the prognostication needs to be improved.

GREB1-CTNNB1 fusion transcript detected by RNA-sequencing in a uterine tumor resembling ovarian sex cord tumor (UTROSCT): A novel CTNNB1 rearrangement.

Mutations of CTNNB1 have been implicated in tumorigenesis in many organs. However, tumors harboring a CTNNB1 translocation are extremely rare and this translocation has never been reported in a uterine mesenchymal neoplasm. We report a novel translocation t(2;3)(p25;p22) involving the GREB1 (intron 8) and CTNNB1 (exon 3) in a uterine tumor resembling ovarian sex cord tumor (UTROSCT), which exhibited extrauterine metastasis. The translocation detected by RNA-sequencing was validated by RT-PCR, and resulted in nuclear expression of ß-catenin. Juxtapositioning with GREB1, which is overexpressed in response to estrogens, resulted in overexpression of a truncated and hypophosphorylated nuclear ß-catenin in the primary and recurrent tumors. This accumulation of nuclear ß-catenin results in a constitutive activation of the Wnt/ß-catenin signaling pathway with a major oncogenic effect. The CTNNB1 gene fusion, promoted by an estrogen-responsive gene (GREB1), could be a potential driver of tumorigenesis in this case and a therapeutic target with adapted inhibitors. RT-PCR and immunohistochemistry performed on 11 additional UTROSCTs showed no CTNNB1 fusion transcript or nuclear ß-catenin immunoreactivity.

Genomic and transcriptomic comparison of post-radiation versus sporadic sarcomas.

Post-radiation sarcomas are rare secondary cancers arising from radiation therapies. To date, few genetic specificities have been described for such malignancies and the oncogenesis of sarcomas with complex genetics (both sporadic and post-radiation) remains largely misunderstood. We performed genomic and transcriptomic analyses on 77 post-radiation sarcomas using DNA-array and RNA sequencing. Consequently, we were able to investigate changes in copy number variations, transcriptome profiling, fusion gene expression, and mutational landscapes. We compare these data to a reference cohort of 93 sporadic sarcomas. At genomic level, similar chromosomal complexity was observed both in post-radiation and sporadic sarcomas with complex genetics. We found more frequent CDKN2A and CDKN2B (coding for p14/p16 and p15 proteins, respectively; at 9p21.3) losses in post-radiation (71%) than in sporadic tumors (39%; P = 6.92e-3). Among all detected fusion genes and punctual variations, few specificities were observed between these groups and such alterations are not able to drive a strong and specific oncogenesis. Recurrent MYC amplifications (96%) and KDR variants (8%) were detected in post-radiation angiosarcomas, in agreement with the literature. Transcriptomic analysis of such angiosarcomas revealed two distinct groups harboring different genomic imbalances (in particular gains of 17q24.2-17qter) with different clinical courses according to patient's vital status. Differential gene expression analysis permitted to focus on the immune response as a potential actor to tumor aggressiveness. Histochemistry validated a lower inflammation and lower immune infiltrate at tumor periphery for highly aggressive angiosarcomas. Our results provide new genomic and transcriptomic information about post-radiation sarcomas. The techniques we used (RNA-seq and DNA-arrays) did not highlight major differences in sarcomas with complex genetics depending on the radiation context, revealing similar patterns of transcriptomic profiles and chromosomal copy number variations. Additional characterizations, particularly whole genome sequencing, could measure changes in DNA following radiation therapy in such malignancies and may precise their oncogenesis.

Genome profiling is an efficient tool to avoid the STUMP classification of uterine smooth muscle lesions: a comprehensive array-genomic hybridization analysis of 77 tumors.

The diagnosis of a uterine smooth muscle lesion is, in the majority of cases, straightforward. However, in a small number of cases, the morphological criteria used in such lesions cannot differentiate with certainty a benign from a malignant lesion and a diagnosis of smooth muscle tumor with uncertain malignant potential (STUMP) is made. Uterine leiomyosarcomas are often easy to diagnose but it is difficult or even impossible to identify a prognostic factor at the moment of the diagnosis with the exception of the stage. We hypothesize, for uterine smooth muscle lesions, that there is a gradient of genomic complexity that correlates to outcome. We first tested this hypothesis on STUMP lesions in a previous study and demonstrated that this 'gray category' could be split according to genomic index into two groups. A benign group, with a low to moderate alteration rate without recurrence and a malignant group, with a highly rearranged profile akin to uterine leiomyosarcomas. Here, we analyzed a large series of 77 uterine smooth muscle lesions (from 76 patients) morphologically classified as 19 leiomyomas, 14 STUMP and 44 leiomyosarcomas with clinicopathological and genomic correlations. We confirmed that genomic index with a cut-off=10 is a predictor of recurrence (P<0.0001) and with a cut-off=35 is a marker for poor overall survival (P=0.035). For the tumors confined to the uterus, stage as a prognostic factor was not useful in survival prediction. At stage I, among the tumors reclassified as molecular leiomyosarcomas (ie, genomic index ≥10), the poor prognostic markers were: 5p gain (overall survival P=0.0008), genomic index at cut-off=35 (overall survival P=0.0193), 13p loss including RB1 (overall survival P=0.0096) and 17p gain including MYOCD gain (overall survival P=0.0425). Based on these findings (and the feasibility of genomic profiling by array-comparative genomic hybridization), genomic index, 5p and 17p gains prognostic value could be evaluated in future prospective chemotherapy trials.

Alternative PDGFD rearrangements in dermatofibrosarcomas protuberans without PDGFB fusions.

Dermatofibrosarcoma protuberans is underlined by recurrent collagen type I alpha 1 chain-platelet-derived growth factor B chain (COL1A1-PDGFB) fusions but ~ 4% of typical dermatofibrosarcoma protuberans remain negative for this translocation in routine molecular screening. We investigated a series of 21 cases not associated with the pathognomonic COL1A1-PDGFB fusion on routine fluorescence in situ hybridization (FISH) testing. All cases displayed morphological and clinical features consistent with the diagnosis of dermatofibrosarcoma protuberans. RNA-sequencing analysis was successful in 20 cases. The classical COL1A1-PDGFB fusion was present in 40% of cases (n = 8/20), and subsequently confirmed with a COL1A1 break-apart FISH probe in all but one case (n = 7/8). 55% of cases (n = 11/20) displayed novel PDGFD rearrangements; PDGFD being fused either to the 5' part of COL6A3 (2q37.3) (n = 9/11) or EMILIN2 (18p11) (n = 2/11). All rearrangements led to in-frame fusion transcripts and were confirmed at genomic level by FISH and/or array-comparative genomic hybridization. PDGFD-rearranged dermatofibrosarcoma protuberans presented clinical outcomes similar to typical dermatofibrosarcoma protuberans. Notably, the two EMILIN2-PDGFD cases displayed fibrosarcomatous transformation and homozygous deletions of CDKN2A at genomic level. We report the first recurrent molecular variant of dermatofibrosarcoma protuberans involving PDGFD, which functionally mimic bona fide COL1A1-PDGFB fusions, leading presumably to a similar autocrine loop-stimulating PDGFRB. This study also emphasizes that COL1A1-PDGFB fusions can be cytogenetically cryptic on FISH testing in a subset of cases, thereby representing a diagnostic pitfall that pathologists should be aware of.

Expression and role of TYRO3 and AXL as potential therapeutical targets in leiomyosarcoma.

BACKGROUND: Leiomyosarcoma (LMS) are 15% of adult sarcomas and remain seldom curable in metastatic phase. The TAM receptors and their ligands are overexpressed or activated in multiple malignancies, including LMS. METHODS: The TAM receptor and ligand expression was evaluated in LMS cell lines and 358 sarcoma samples by either gene expression or immunohistochemistry. TYRO3 and AXL were knocked down. Crizotinib and foretinib were investigated in vitro. RESULTS: High expression of TYRO3 and AXL was detected in LMS cell lines. TYRO3 or AXL gene knockdown reduced cell proliferation/colony formation. Crizotinib and foretinib decreased TYRO3 and AXL phosphorylation, apoptosis, G2/arrest and reduced colony formation. Immunohistochemistry performed in 107 sarcomas showed higher expression of TYRO3 and GAS6 in LMS vs other sarcomas and nuclear TYRO3 only in LMS. Microarray gene expression performed in 251 sarcomas revealed significantly higher expression of TYRO3 and GAS6 in LMS than other sarcomas. Leiomyosarcoma patients with high expression of GAS6 or PROS1 present a significantly worse PFS. CONCLUSIONS: Leiomyosarcoma patients, especially those whom develop metastasis, express higher levels of TYRO3 and GAS6. Crizotinib and foretinib showed effective antitumour activity in LMS through TYRO3 and AXL deactivation indicating that clinical trials using TYRO3 and AXL inhibitors are warranted in advanced LMS.

Differential motif enrichment analysis of paired ChIP-seq experiments.

BACKGROUND: Motif enrichment analysis of transcription factor ChIP-seq data can help identify transcription factors that cooperate or compete. Previously, little attention has been given to comparative motif enrichment analysis of pairs of ChIP-seq experiments, where the binding of the same transcription factor is assayed under different conditions. Such comparative analysis could potentially identify the distinct regulatory partners/competitors of the assayed transcription factor under different conditions or at different stages of development. RESULTS: We describe a new methodology for identifying sequence motifs that are differentially enriched in one set of DNA or RNA sequences relative to another set, and apply it to paired ChIP-seq experiments. We show that, using paired ChIP-seq data for a single transcription factor, differential motif enrichment analysis identifies all the known key transcription factors involved in the transformation of non-cancerous immortalized breast cells (MCF10A-ER-Src cells) into cancer stem cells whereas non-differential motif enrichment analysis does not. We also show that differential motif enrichment analysis identifies regulatory motifs that are significantly enriched at constrained locations within the bound promoters, and that these motifs are not identified by non-differential motif enrichment analysis. Our methodology differs from other approaches in that it leverages both comparative enrichment and positional enrichment of motifs in ChIP-seq peak regions or in the promoters of genes bound by the transcription factor. CONCLUSIONS: We show that differential motif enrichment analysis of paired ChIP-seq experiments offers biological insights not available from non-differential analysis. In contrast to previous approaches, our method detects motifs that are enriched in a constrained region in one set of sequences, but not enriched in the same region in the comparative set. We have enhanced the web-based CentriMo algorithm to allow it to perform the constrained differential motif enrichment analysis described in this paper, and CentriMo's on-line interface (http://meme.ebi.edu.au) provides dozens of databases of DNA- and RNA-binding motifs from a full range of organisms. All data and output files presented here are available at http://research.imb.uq.edu.au/t.bailey/supplementary\_data/Lesluyes2014.

Reconstructing oral cavity tumor evolution through brush biopsy.

Oral potentially malignant disorders (OPMDs) with genomic alterations have a heightened risk of evolving into oral squamous cell carcinoma (OSCC). Currently, genomic data are typically obtained through invasive tissue biopsy. However, brush biopsy is a non-invasive method that has been utilized for identifying dysplastic cells in OPMD but its effectiveness in reflecting the genomic landscape of OPMDs remains uncertain. This pilot study investigates the potential of brush biopsy samples in accurately reconstructing the genomic profile and tumor evolution in a patient with both OPMD and OSCC. We analyzed single nucleotide variants (SNVs), copy number aberrations (CNAs), and subclonal architectures in paired tissue and brush biopsy samples. The results showed that brush biopsy effectively captured 90% of SNVs and had similar CNA profiles as those seen in its paired tissue biopsies in all lesions. It was specific, as normal buccal mucosa did not share these genomic alterations. Interestingly, brush biopsy revealed shared SNVs and CNAs between the distinct OPMD and OSCC lesions from the same patient, indicating a common ancestral origin. Subclonal reconstruction confirmed this shared ancestry, followed by divergent evolution of the lesions. These findings highlight the potential of brush biopsies in accurately representing the genomic profile of OPL and OSCC, proving insight into reconstructing tumor evolution.

Age-specific induction of mutant p53 drives clonal hematopoiesis and acute myeloid leukemia in adult mice.

The investigation of the mechanisms behind p53 mutations in acute myeloid leukemia (AML) has been limited by the lack of suitable mouse models, which historically have resulted in lymphoma rather than leukemia. This study introduces two new AML mouse models. One model induces mutant p53 and Mdm2 haploinsufficiency in early development, showing the role of Mdm2 in myeloid-biased hematopoiesis and AML predisposition, independent of p53. The second model mimics clonal hematopoiesis by inducing mutant p53 in adult hematopoietic stem cells, demonstrating that the timing of p53 mutation determines AML vs. lymphoma development. In this context, age-related changes in hematopoietic stem cells (HSCs) collaborate with mutant p53 to predispose toward myeloid transformation rather than lymphoma development. Our study unveils new insights into the cooperative impact of HSC age, Trp53 mutations, and Mdm2 haploinsufficiency on clonal hematopoiesis and the development of myeloid malignancies.

The History of Chromosomal Instability in Genome-Doubled Tumors.

Tumors frequently display high chromosomal instability and contain multiple copies of genomic regions. Here, we describe Gain Route Identification and Timing In Cancer (GRITIC), a generic method for timing genomic gains leading to complex copy number states, using single-sample bulk whole-genome sequencing data. By applying GRITIC to 6,091 tumors, we found that non-parsimonious evolution is frequent in the formation of complex copy number states in genome-doubled tumors. We measured chromosomal instability before and after genome duplication in human tumors and found that late genome doubling was followed by an increase in the rate of copy number gain. Copy number gains often accumulate as punctuated bursts, commonly after genome doubling. We infer that genome duplications typically affect the landscape of copy number losses, while only minimally impacting copy number gains. In summary, GRITIC is a novel copy number gain timing framework that permits the analysis of copy number evolution in chromosomally unstable tumors. Significance: Complex genomic gains are associated with whole-genome duplications, which are frequent across tumors, span a large fraction of their genomes, and are linked to poorer outcomes. GRITIC infers when these gains occur during tumor development, which will help to identify the genetic events that drive tumor evolution. See related commentary by Taylor, p. 1766.

Pan-cancer subclonal mutation analysis of 7,827 tumors predicts clinical outcome.

Intra-tumor heterogeneity is an important driver of tumor evolution and therapy response. Advances in precision cancer treatment will require understanding of mutation clonality and subclonal architecture. Currently the slow computational speed of subclonal reconstruction hinders large cohort studies. To overcome this bottleneck, we developed Clonal structure identification through Pairwise Penalization, or CliPP, which clusters subclonal mutations using a regularized likelihood model. CliPP reliably processed whole-genome and whole-exome sequencing data from over 12,000 tumor samples within 24 hours, thus enabling large-scale downstream association analyses between subclonal structures and clinical outcomes. Through a pan-cancer investigation of 7,827 tumors from 32 cancer types, we found that high subclonal mutational load (sML), a measure of latency time in tumor evolution, was significantly associated with better patient outcomes in 16 cancer types with low to moderate tumor mutation burden (TMB). In a cohort of prostate cancer patients participating in an immunotherapy clinical trial, high sML was indicative of favorable response to immune checkpoint blockade. This comprehensive study using CliPP underscores sML as a key feature of cancer. sML may be essential for linking mutation dynamics with immunotherapy response in the large population of non-high TMB cancers.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring.

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Genome-wide in silico prediction of gene expression.

MOTIVATION: Modelling the regulation of gene expression can provide insight into the regulatory roles of individual transcription factors (TFs) and histone modifications. Recently, Ouyang et al. in 2009 modelled gene expression levels in mouse embryonic stem (mES) cells using in vivo ChIP-seq measurements of TF binding. ChIP-seq TF binding data, however, are tissue-specific and relatively difficult to obtain. This limits the applicability of gene expression models that rely on ChIP-seq TF binding data. RESULTS: In this study, we build regression-based models that relate gene expression to the binding of 12 different TFs, 7 histone modifications and chromatin accessibility (DNase I hypersensitivity) in two different tissues. We find that expression models based on computationally predicted TF binding can achieve similar accuracy to those using in vivo TF binding data and that including binding at weak sites is critical for accurate prediction of gene expression. We also find that incorporating histone modification and chromatin accessibility data results in additional accuracy. Surprisingly, we find that models that use no TF binding data at all, but only histone modification and chromatin accessibility data, can be as (or more) accurate than those based on in vivo TF binding data. AVAILABILITY AND IMPLEMENTATION: All scripts, motifs and data presented in this article are available online at http://research.imb.uq.edu.au/t.bailey/supplementary_data/McLeay2011a.