Search Results - EGA European Genome-Phenome Archive

Tumor suppressor miR-133a modulates the prostate cancer epigenome by repressing BAZ2A

Our study revealed a new subgroup of prostate tumors defined by the expression of a novel prognostic marker, BAZ2A, and the first time to show the establishment of CIMP via an microRNA-epigenetic interaction in prostate tumors.

Study EGAS00001000568

Novel Epigenetic Markers for Toxicity after Intraoperative and Conventional Radiotherapy for Breast Cancer

Ionizing radiation is a common treatment option for cancer but its use is limited by the unpredictable and highly heterogeneous onset of late side effects, especially radiation-induced fibrosis. Clinically applicable biomarkers and effective treatments for radiation fibrosis are currently unavailable. In order to identify novel markers we ran a genome-wide DNA methylation screen in primary dermal fibroblasts obtained from breast cancer patients before intraoperative radiotherapy. Cells from patients developing fibrosis within a three-year follow up were compared to those without fibrosis (12 individuals per group). Illumina Infinium HumanMethylation450 BeadChip Technology was used.

Study EGAS00001001279

ABIS_1_MeDIP-seq

We investigated changes in whole genome DNA-methylation in a group of 48 healthy five-year-old children differentiated by their hair cortisol levels. By combining this novel biomarker for chronic stress with methylated DNA immunoprecipitation sequencing (MeDIP-seq), we found that high cortisol were associated with a genome-wide decrease in DNA-methylation targeted to specific zinc finger transcription factor binding sites, repetitive elements and genes important for neurodevelopment and calcium transport. Large scale hypomethylation of repetitive regions, as well as a deficient calcium transport, have been reported in many stress related diseases, as well as in many cancers and biological aging. Thus our findings may aid in our understanding of how susceptibility to these diseases develops.

Study EGAS00001001099

RB Loss in Resistant EGFR Mutant Lung Adenocarcinomas that Transform to Small Cell Lung Cancer

Tyrosine kinase inhibitors (TKIs) are effective treatments for non-small cell lung cancers (NSCLCs) with epidermal growth factor receptor (EGFR) mutations. However, they do not lead to cures, and, on average, relapse occurs typically occurs afterafter an average of one year of continuous treatment. In a subset of patients, aA fundamental histological transformation from NSCLC to small cell lung cancer (SCLC) is observed in a subset of the resistant cancers, but the molecular changes associated with this transformation remain unknown. Analysis of a cohort of tumor samples and cell lines derived from resistant EGFR mutant patients with SCLC transformation revealed that RB is lost in 100% of these SCLC transformed cases, but rarely in those that remain NSCLC. Global changes in gene expression, including increased neuroendocrine marker expression and absence of EGFR expression, are observed in resistant EGFR mutant cancers that transformed to SCLC. Consistent Further, with their similarities to classical SCLC at the genetic and gene expression levelincreased neuroendocrine marker and decreased EGFR expression as well as greater sensitivity to BCL2 family inhibition are observed in resistant SCLC transformed cancers compared to resistant NSCLCs. , cell lines derived from resistant EGFR mutant SCLC biopsies are substantially more sensitive to Bcl-2 family inhibition with ABT-263 compared to those derived from resistant EGFR mutant NSCLCs. Together, these findings suggest that despite developing initially as EGFR mutant adenocarcinomas, this subset of resistant cancers ultimately adopt many of the molecular and phenotypic characteristics of classical SCLC.

Study EGAS00001001102

Combined hereditary and somatic mutations of replication error repair genes result in rapid onset of ultra-hypermutated cancers

DNA replication-associated mutations are repaired by two components: polymerase proofreading and mismatch repair. The mutational consequences of disruption to both repair components in humans are not well studied. We sequenced cancer genomes from children with inherited biallelic mismatch repair deficiency (bMMRD). High-grade bMMRD brain tumors exhibited massive numbers of substitution mutations (>250/Mb) – higher than all childhood and most cancers (>7000 analyzed). All ultra-hypermutated bMMRD cancers acquired early somatic driver mutations in DNA polymerases epsilon or delta. The ensuing mutation signatures and numbers are unique and diagnostic of childhood germline bMMRD (p<10e-13). Sequential tumor biopsy analysis revealed that bMMRD/polymerase mutant cancers rapidly amass an excess of simultaneous mutations (~600 mutations/cell division), reaching but not exceeding ~20,000 exonic mutations in <6 months. Thus implying a threshold compatible with cancer cell survival. We suggest a new mechanism of cancer progression in which mutations develop in a rapid burst after ablation of replication repair.

Study EGAS00001001112

Reproducibility of variant calls in replicate next generation sequencing experiments

In this study, we examined the reproducibility of nucleotide variant calls in replicate sequencing experiments of the same genomic DNA. We performed targeted sequencing of all known human protein kinase genes (kinome) (~3.3 Mb) using the SOLiD v4 platform. Seventeen breast cancer samples were sequenced in duplicate (n=14) or triplicate (n=3) to assess concordance of all calls and single nucleotide variant (SNV) calls.

Study EGAS00001000826

Whole exome sequencing reveals the mutational spectrum of testicular germ cell tumours

Testicular germ cell tumours (TGCTs) are the most common cancer in young men. Here we perform whole exome sequencing of 42 TGCTs to comprehensively study the mutational profile of TGCT. The mutation rate is uniformly low in all of the tumours (mean 0.5 mutations per megabase [Mb]) as compared to the common cancers, consistent with the embryological origin of TGCT. In addition to expected copy number gain of chromosome 12p and mutation of KIT we identify recurrent mutations in the tumour suppressor gene CDC27 (11.9%). Copy number analysis reveals recurring amplification of the spermatocyte development gene FSIP2 (15.3%) and a 0.4Mb region at Xq28 (15.3%). Two treatment-refractory patients are shown to harbour XRCC2 mutations, a gene strongly implicated in defining cisplatin resistance. Our findings provide further insights into genes involved in the development and progression of TGCT.

Study EGAS00001001084

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci

Multiple regulatory elements distant to their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we have developed a capture Hi-C (cHi-C) approach to allow an agnostic characterisation of these physical interactions on a genome-wide scale. Single nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci has allowed us to identify key long-range chromatin interactions in cis and trans involving these loci.

Study EGAS00001001085

Therapeutic Resistance to PI3K-alpha Inhibitors

We have studied the tumor genomic evolution in a patient with metastatic breast cancer bearing an activating PIK3CA mutation. The patient was treated with the PI3Kα inhibitor BYL719 and achieved a lasting clinical response on BYL719, but eventually progressed to treatment and died shortly thereafter. A rapid autopsy was performed and a total of 14 metastatic lesions were collected for further analysis. In order to identify possible genetic determinants of acquired resistance to PI3Kα inhibition, we took a three-step approach. First, we examined both the primary tumor (before BYL719 treatment) and the new lung metastasis by whole genome sequencing (DNA from the spleen was used as a normal control). Then, we analyzed the primary tumor, lung metastasis, and the peri-aortic lesion that remained stable (responding) at the time of progression to BYL719 therapy by whole exome sequencing. Finally, to confirm and expand our findings, we sequenced the primary tumor and all the metastatic lesions to >500-fold coverage using a custom targeted deep-sequencing assay, MSK-IMPACT. Using this 3-step approach, we have demonstrated that patient’s lesion underwent parallel genetic evolution at multiple metastatic sites with different PTEN genomic alterations leading to a convergent PTEN- null phenotype resistant to PI3Kα inhibition. In order to expand our observations, we analyzed paired samples (pre-treatment and at progression) from six additional patients enrolled in the BYL719 trial. Acquired bi-allelic loss of PTEN was found in one additional patient treated with BYL719 whereas in two patients PIK3CA mutations present in the primary tumor were no longer detected at the time of progression. To functionally characterize our findings, inducible PTEN knockdown in sensitive cells resulted in resistance to BYL719, while simultaneous PI3K-p110β blockade reverted this resistance phenotype, both in cell lines and in PTEN-null xenografts derived from our patient. We conclude that parallel genetic evolution of separate sites with different PTEN genomic alterations leads to a convergent PTEN-null phenotype resistant to PI3Kα inhibition.

Study EGAS00001000991

Whole exome sequencing of papillary thyroid carcinoma in the Chinese population

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Here we submitted the sequencing results for PTC using 91 tumor-normal pairs through exome sequencing in the Chinese Han population.

Study EGAS00001001268

Comprehensive analysis of atypical teratoid rhabdoid tumour (ATRT) using genomic, epigenomic and transcriptomic techniques.

ATRTs (Atypical teratoid rhabdoid tumours) represent one of the most aggressive pediatric brain cancers, but intriguingly exhibit few other recurrently mutated loci except for SMARCB1/hSNF5. We integrated whole genome (n=15), exome, copy number, gene expression and methylation analyses to comprehensively interrogate 64 ATRTs and observed that structural events were relatively frequent in the ATRT genome (~3 tumour). In addition to SMARCB1, which was targeted by structural events in a majority (49/64) of tumours, recurrent structural alterations targeting the LRP1B, CDH13, BCR and MKL1 loci were observed. We observed novel translocation events, including 2 targeting SMARCB1 which were detected only by NGS analyses. Significantly, integration of gene expression and methylation profiles with genomic analyses revealed ATRT comprise two sub-groups with distinct clinical and genetic features. Group 1 tumours were characterized by supra-tentorial brain location, focal intragenic alterations of SMARCB1 and a pro-neural gene expression signature with evidence of NOTCH pathway activation. While Group 2 ATRTs were characterized by infra-tentorial brain location, broad SMARCB1 alterations and an activation of the BMP signalling pathway. Our findings highlight epigenetic mechanisms as important determinants of ATRT tumour phenotypes, and demonstrate for the first time that ATRTs arising in different anatomical compartments comprise distinct molecular and therapeutic sub-groups.

Study EGAS00001000506

An integrated multi-omic cellular atlas of human breast cancers

Study EGAS00001005173

Copy-number signatures and mutational processes in ovarian carcinoma

This study contains the data underlying our approach for copy-number signature identification in high-grade serous ovarian cancer. The samples studied are part of the BriTROC-1 consortium, a national study into relapsed high-grade serous ovarian cancer. The dataset consists of shallow-whole genome sequencing of 300 samples from 143 patients with match tagged-amplicon sequencing of TP53. The study also contains deeper whole-genome sequencing of 48 cases.

Study EGAS00001002557

These are from Korean HCC samples with exome sequencing

Hepatic resection is the most curative treatment option for early-stage hepatocellular carcinoma, but is associated with a high recurrence rate, which exceeds 50% at 5 years after surgery. Understanding the genetic basis of hepatocellular carcinoma at surgically curable stages may enable the identification of new molecular biomarkers that accurately identify patients in need of additional early therapeutic interventions. Whole exome sequencing and copy number analysis was performed on 231 hepatocellular carcinomas (72% with hepatitis B viral infection) that were classified as early-stage hepatocellular carcinomas, candidates for surgical resection. Recurrent mutations were validated by Sanger sequencing. Unsupervised genomic analyses identified an association between specific genetic aberrations and postoperative clinical outcomes. Recurrent somatic mutations were identified in 9 genes, including TP53, CTNNB1, AXIN1, RPS6KA3, and RB1. Recurrent homozygous deletions in FAM123A, RB1, and CDKN2A, and high-copy amplifications in MYC, RSPO2, CCND1, and FGF19 were detected. Pathway analyses of these genes revealed aberrations in the p53, Wnt, PIK3/Ras, cell cycle, and chromatin remodelling pathways. RB1 mutations were significantly associated with cancer-specific and recurrence-free survival after resection (p = 0.016 and p = 0.001, respectively). FGF19 amplifications, known to activate Wnt signalling, were mutually exclusive with CTNNB1 and AXIN1 mutations, and significantly associated with cirrhosis (p = 0.017). RB1 mutations can be used as a prognostic molecular biomarker for resectable hepatocellular carcinoma. Further study is required to investigate the potential role of FGF19 amplification in driving hepatocarcinogenesis in patients with liver cirrhosis and to investigate the potential of anti-FGF19 treatment in these patients.

Study EGAS00001000604

Epigenomic profile of diverse cancer

Epigenomic profile of diverse cancer.

Study EGAS00001004352

Search for new loci and low-frequency variants influencing glioma risk by exome-array analysis

To identify protein altering variants (PAVs) for glioma we analysed Illumina HumanExome BeadChip exome array data on 1,882 glioma cases and 8,079 controls from three independent European populations. In addition to single variant tests we incorporated information on the predicted functional consequences of PAVs and analysed sets of genes with a higher likelihood of having a role in glioma on the basis of the profile of somatic mutations documented by large-scale sequencing initiatives. Globally tThere was a strong relationship between effect size and SNPs predicted to be damaging (P=2.29x10-49); however, these variants which are most likely to impact on risk, are rare (MAF<5%). While no single variant showed an association which was statistically significant at the genomewide threshold a number of represented promising Notable associations - were seen with the BRCA2:c.9976A>T, p.(Lys3326Ter)BRCA2 p.Lys332X variant, which has been shown documented to influence breast and lung cancer risks (OR=2.3, P=4.00x10-4 for glioblastoma [GBM]) and IDH2:c.782G>A, p.(Arg261His)IDH2 p.Arg261His (OR=3.21, P=7.67x10-3, for non-GBM). Additionally, gene burden tests revealed a statistically significant association for HARS2 and risk of GBM (P=2.20x10-6). Genome scans of low frequency PAVs represent a complementary strategy to identify disease-causing variants compared with scans based on tagSNPs. Strategies to lessen the multiple testing burden by restricting analysis to PAVs with higher priors affords an opportunity to maximise study power.

Study EGAS00001001258

Robust methylation-based classification of brain tumours using nanopore sequencing

Background: DNA methylation-based classification of cancer provides a comprehensive molecular approach to diagnose tumours. In fact, DNA methylation profiling of human brain tumours already profoundly impacts clinical neuro-oncology. However, current implementations using hybridization microarrays are time-consuming and costly. We recently reported on shallow nanopore whole-genome sequencing for rapid and cost-effective generation of genome-wide 5-methylcytosine profiles as input to supervised classification. Here, we demonstrate that this approach allows to discriminate a wide spectrum of primary brain tumours. Results: Using public reference data of 82 distinct tumour entities, we performed nanopore genome sequencing on N=382 tissue samples covering 46 brain tumour (sub)types. Using bootstrap sampling in a cohort of N = 56 cases, we find that a minimum set of 1,000 random CpG features is sufficient for high-confidence classification by ad hoc random forests. We implemented score recalibration as confidence measure for interpretation in a clinical context and empirically determined a platform-specific threshold in a randomly sampled discovery cohort (N = 185). Applying this cut-off to an independent validation series (N = 184) yielded 148 classifiable cases (sensitivity 80.4%) and demonstrated 100 % specificity. Cross-lab validation demonstrated robustness with concordant results across four laboratories in 10/11 (90.9%) cases. In a prospective benchmarking (N = 15), median time to results was 21.1 hours. Conclusions: In conclusion, nanopore sequencing allows robust and rapid methylation-based classification across the full spectrum of brain tumours. Platform-specific confidence scores facilitate clinical implementation for which prospective evaluation is warranted and ongoing.

Study EGAS00001006540

Clinical Implications of Genomic Alterations in the Tumor and Circulation of Pancreatic Cancer Patients

Pancreatic adenocarcinoma has the worst mortality of any solid cancer. To evaluate the clinical implications of genomic alterations in this tumor type, we performed whole-exome analyses of 24 tumors, targeted genomic analyses of 77 tumors, and used non-invasive approaches to examine tumor-specific mutations in the circulation of these patients. These analyses reveal somatic mutations in chromatin regulating genes MLL, MLL2, MLL3, and ARID1A in 20% of patients that were associated with improved survival. We observe alterations in genes with potential clinical utility in over a third of cases. Liquid biopsy analyses demonstrate that 43% of patients with localized disease have detectable circulating tumor DNA (ctDNA) at diagnosis. Detection of ctDNA after resection predicts clinical relapse and poor outcome, and recurrence by ctDNA is detected 6.5 months earlier than with CT imaging. These observations provide genetic predictors of outcome in pancreatic cancer and have implications for new avenues of therapeutic intervention.

Study EGAS00001001257

The Genomic Landscape of Response to EGFR Blockade in Colorectal Cancer

Colorectal cancer (CRC) is the third most common cancer world-wide with 1.2 million patients diagnosed yearly. In late stage CRC, the most commonly used targeted therapies are monoclonal antibodies cetuximab and panitumumab, which inactivate EGFR. Recent studies have identified alterations in KRAS and other genes as likely mechanisms of primary and secondary resistance to anti-EGFR antibody therapy. Despite these efforts, additional mechanisms of resistance to EGFR blockade are thought to be present in CRC and little is known about determinants of sensitivity to this therapy. To examine the effect of somatic genetic changes in CRC on response to anti-EGFR antibody therapy, we performed complete exome sequence and copy number analyses of 129 patient-derived tumorgrafts and targeted genomic analyses of 55 patient tumors, all of which were KRAS wild-type. We analyzed the response of tumors to anti-EGFR antibody blockade in tumorgraft models or in clinical settings. In addition to previously identified genes, we detected mutations in ERBB2, EGFR, FGFR1, PDGFRA, and MAP2K1 as potential mechanisms of primary resistance to this therapy. Novel alterations in the ectodomain of EGFR were identified in patients with acquired resistance to EGFR blockade. Amplifications and sequence changes in the tyrosine kinase receptor adaptor gene IRS2 were identified in tumors with increased sensitivity to anti-EGFR therapy. Therapeutic resistance to EGFR blockade could be overcome in tumorgraft models through combinatorial therapies targeting actionable genes. These analyses provide a systematic approach to evaluate response to targeted therapies in human cancer, highlight new mechanisms of responsiveness to anti-EGFR therapies, and provide new avenues for intervention in the management of CRC.

Study EGAS00001001305

Chromothripsis followed by circular recombination drives oncogene amplification in human cancer

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. We here identified a type of amplification, termed as “seismic amplification”, that is characterized by multiple rearrangements and discontinuous copy number levels. Seismic amplifications occurred in 9.9% (274/2,756) of cases across 38 cancer types and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification revealed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and over-expression in a substantial fraction of human malignancies.

Study EGAS00001005424

Genomic analysis of HPV positive versus HPV negative esophageal adenocarcinoma

High-risk human papillomavirus (hr-HPV) has been implicated in a subset of patients with esophageal adenocarcinoma (EAC). We therefore hypothesized that HPV associated EAC may have distinct genomic aberrations compared with viral negative esophageal cancer as shown in Head and Neck cancer and cervical cancer. Whole exome sequencing (WES) was performed to explore the mutational landscape and potential molecular signature of HPV positive versus HPV negative EAC, which includes four hr-HPV-positive and 8 HPV-negative treatment naive fresh-frozen EAC tissue specimens and matched normal tissues.

Study EGAS00001001340

Whole exome sequencing data of germline and two independent primary leukemias of five patients

The contribution of genetic predisposing factors to the development of pediatric acute lymphoblastic leukemia (ALL), the most frequently diagnosed cancer in childhood, has not been fully elucidated. Children presenting with multiple de novo leukemias are more likely to suffer from genetic predisposition. Here, we selected five of these patients and analyzed the mutational spectrum of normal and malignant tissues. In two patients, we identified germline mutations in TYK2, a member of the JAK tyrosine kinase family. These mutations were located in two adjacent codons of the pseudokinase domain (p.Pro760Leu and p.Gly761Val). In silico modeling revealed that both mutations affect the conformation of this auto-regulatory domain. Consistent with this notion, both germline mutations promote TYK2 autophosphorylation and activate downstream STAT family members, which could be blocked with the JAK kinase inhibitor I. These data indicate that germline activating TYK2 mutations predispose to the development of ALL.

Study EGAS00001001889

Genetic and epigenetic characterization of adenoid cystic carcinoma

Translocation events are frequent in cancer and may create chimeric fusions or regulatory rearrangements that drive oncogene overexpression. Although regulatory rearrangements are increasingly recognized in hematopoietic and even solid tumors, the underlying mechanisms remain obscure. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a unifying theme in adenoid cystic carcinoma (ACC). Whole genome sequencing data and chromatin state maps for 13 primary tumors and xenografts reveal distinct chromosomal rearrangements that juxtapose super-enhancer clusters to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers physically interact with the MYB promoter. Remarkably, the MYB product itself binds to the translocated enhancers, thus creating a positive feedback loop that sustains its own expression. MYB also binds a large number of active enhancers that drive opposing regulatory programs in alternate cell lineages in ACC. MYB cooperates with the transcription factor TP63 in the myoepithelial component of the tumors, but promotes a Notch program in the luminal epithelial component. Bromodomain inhibitors, which disrupt enhancer function, slow tumor growth in ACC primagraft models in vivo, but are ineffective against high grade tumors that harbor activating mutations in the Notch pathway. Thus, our study identifies super-enhancer translocations as a unifying feature of ACC, and provides insight into the mechanism by which sustained MYB overexpression drives alternate cell fates in this disease.Processed ChIP-seq data is available on GEO under accession number GSE76465.

Study EGAS00001001457

The role of the mammalian SWI/SNF chromatin remodeling complex in neuroendocrine prostate cancer

We sought to determine the effects of SMARCA4 and SMARCA2 depletion in prostate cancer cell lines. We performed siRNA-mediated knock-down of SMARCA4 and SMARCA2 in an androgen-sensitive (LNCaP) cell line and in a castration-resistant prostate cancer (CRPC)-adenocarcinoma cell line (22Rv1) and compared global transcriptional alterations using RNA-seq.

Study EGAS00001004177

MINCR is a MYC-induced lncRNA able to modulate MYC’s transcriptional network in Burkitt lymphoma cells

Despite the established role of the transcription factor MYC in cancer, little is known about the involvement of lncRNAs in mediating MYC’s function. Here we have intersected RNA-sequencing data from MYC-inducible cell lines, from a cohort of 91 mature B-cell lymphomas and from sorted germinal-center B-cells. By this approach, we identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated by MYC in the model cell lines. Among them we focused on a lncRNA that we named MINCR, showing a strong correlation with MYC expression in MYC-positive lymphomas and in pancreatic ductal adenocarcinomas. RNAi experiments showed that MINCR controls cellular viability by influencing the expression of MYC-regulated cell cycle genes. Finally we provide evidences that down-regulation of AURKA, AURKB and CTD1 may explain the reduction in cellular proliferation observed upon MINCR knock-down. We therefore suggest that MINCR acts as a modulator of MYC transcriptional control of cell cycle genes.

Study EGAS00001001199

The aim of this project is to identify, on 15 French Caucasian and 10 African-Caribbean men, through an integrative approach of DNA sequencing and transciptomic analyses, relevant genomic events that characterize or allow targeting the various phenotypes of aggressiveness of early stages of prostate cancer.

Every year about 70,000 new cases of prostate cancer are diagnosed in France and 8,700 men died from this disease. But the incidence and the rate of mortality per 100,000 men in the French West Indies are respectively 1.5x and 2x upper than the rates observed in the metropolitan France. The aim of this project is to identify, on 15 French Caucasian and 10 African-Caribbean men, through an integrative approach of DNA sequencing and transciptomic analyses, relevant genomic events that characterize or allow targeting the various phenotypes of aggressiveness of early stages of prostate cancer.

Study EGAS00001002176

Breast Cancer - Subtype defined by an amplification of the HER2 gene

Clinically positive HER2 (cHER2+) breast cancers (BC) are defined by the amplification and overexpression of the human epidermal growth factor receptor (HER2/ERBB2) gene on chromosome region 17q12 and represent 15% of all BCs. Although anti-HER2 targeted therapies have improved BC’s natural history, the magnitude of their benefit varies widely between patients and deciphering the genomic heterogeneity of cHER2+ BC may provide a basis to better understand patients’ outcome. As part of the ICGC Breast Cancer Working Group effort, we combined whole genome sequencing and transcriptomic analyses of 64 cHER2+ primary invasive carcinomas.

Study EGAS00001001431

Chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in paediatric high grade glioma cell lines KKNS4 and SF188 with and without a G34V mutation in H3F3A

Glioblastomas of children and young adults have a median survival of only 12-15months and are clinically and biologically distinct from histologically similar cancers in older adults1. They are defined by highly specific mutations in the gene encoding the histone H3.3 variant H3F3A2, occurring either at or close to key residues marked by methylation for regulation of transcription – K27 and G34. We performed chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in order to test the hypothesis that, rather than total H3K36me3, the G34V mutation may instead result in differential binding of the trimethyl mark throughout the genome. Compared with H3F3A wild-type SF188 paediatric GBM cells, H3K36me3 was found to be significantly differentially bound in KNS42 cells at 5130 distinct regions of the genome corresponding to 156 genes. Concurrent whole genome DNA sequencing demonstrated that this was not confounded by copy number alterations. The cerebral hemispheric-specific G34 mutation drives a distinct expression signature through differential genomic binding of the K36 trimethylation mark (H3K36me3). The transcriptional program induced recapitulates that of the developing forebrain, and involves numerous markers of stem cell maintenance, cell fate decisions and self-renewal. Critically, H3F3A G34 mutations cause profound upregulation of MYCN, a potent oncogene which is causative of glioblastomas when expressed in the correct developmental context. This driving aberration is selectively targetable in this patient population by inhibiting kinases responsible for stabilisation of the protein.

Study EGAS00001001437

Single Cell RNAseq of blood and tumor from renal cancer patients

Single-cell RNASeq (scRNAseq) showed that IL8 is primarily expressed in circulating and intratumoral myeloid cells and high IL8 expression was associated with the downregulation of the antigen presentation machinery in myeloid cells.

Study EGAS00001004230

Inherited damaging mutations in immune-related genes favour the development of genetically heterogeneous synchronous colorectal cancer.

Synchronous colorectal cancers are physically separated tumours that develop simultaneously in the intestine. To understand how common genetic and environmental backgrounds influence the development of multiple tumours, we conduct a genome-wide comparative analysis of somatic alterations of 20 syCRCs from 10 patients. We show that syCRCs have independent genetic origins, acquire discordant somatic alterations, and follow different clonal developments.

Study EGAS00001001461

Sequencing of serial plasma and multiregional tumor samples in a patient with metastatic breast cancer

Circulating tumor DNA analysis can be used to track tumor burden and analyze cancer genomes non-invasively but the extent to which it represents metastatic heterogeneity is unknown. Here, we follow a patient with metastatic ER and HER2 positive breast cancer receiving two lines of targeted therapy over 3 years. We characterize genomic architecture and infer clonal evolution in 8 tumor biopsies and 9 plasma samples collected over 1,193 days of clinical follow-up using exome and targeted amplicon sequencing.

Study EGAS00001001466

Fixative optimisation study for BRITROC project

The British Translational Research Ovarian Cancer Collaborative (BriTROC), spanning more than 10 centres across UK, is a prospective observational project focused on relapsed high-grade serous ovarian cancer (HGSOC). The first study, as a part of the project, is the Fixative optimisation study, which analyses the utility for next-generation sequencing of UMFIX samples compared to NBF-fixed and fresh frozen samples. Various DNA and sequencing quality analyses were performed to compare fixation methods in bulk tumour and biopsy samples.

Study EGAS00001001433

The genomic landscape of serrated lesion of the colorectum

Serrated lesions of the colorectum are the precursors of 15-30% of colorectal cancers (CRCs). These lesions have a peculiar morphological appearance, and they are more difficult to detect than conventional adenomatous polyps. In this study we sought to define the genomic landscape of these lesions using high depth targeted sequencing.

Study EGAS00001005648

Data from the paper Context-specific Effects of TGFβ/SMAD3 in Cancer Are Modulated by the Epigenome. Tufegdzic et al, Cell Reports 2015

The transforming growth factor beta (TGFβ) signalling pathway exerts opposing effects on cancer cells, acting either as a tumor-promoter or as a tumor-suppressor. Here we show these opposing effects are a result of the synergy between SMAD3, a downstream effector of TGFβ signalling, and the distinct epigenomes of breast tumor initiating cells (BTICs). These effects of TGFβ are associated with distinct gene expression programs, but genomic SMAD3 binding patterns are highly similar in the BTIC-promoting and BTIC- suppressing contexts. Our data show cell type-specific epigenetic configurations provide a modulatory layer by determining accessibility of genes to regulation by TGFβ/SMAD3. LBH, one such context-specific target gene, is regulated according to its DNA methylation status, and is crucial for TGFβ-dependent promotion of BTICs. Overall, these results reveal that the epigenome plays a central and previously overlooked role in shaping the context-specific effects of TGFβ in cancer.

Study EGAS00001001570

Genomic analysis reveals novel secondary drivers and progression pathways in skin basal cell carcinoma

Basal cell carcinoma of the skin (BCC) is the most common malignant neoplasm in humans. BCC is primarily driven by the aberrant activation of the Sonic Hedgehog (Hh) pathway. However, its extensive phenotypical variation remains to be explained. The genetic profiling of 293 BCCs revealed the highest mutation rate observed in cancer (65 Mutations/Mb), with strong prevalence of UV-light signature mutations. 85% of BCCs harbored mutations in Hh pathway genes: mutually exclusive PTCH1 (73%) and SMO (20%) (P=6.6x10-8), SUFU (8%), and in TP53 (61%). 85% of BCCs also harbored secondary driver mutations in other genes implicated in BCC tumorigenesis. Recurrent secondary driver mutations were observed in MYCN (30%), PPP6C (15%), STK19 (10%), LATS1 (8%), ERBB2 (4%), PIK3CA (2%), RAC1 (1%) and N/K/H-RAS (2%). Loss of function (LoF) and deleterious missense mutations were observed in PTPN14 (23%), RB1 (8%) and FBXW7 (5%). In line with the mutational profiles detected by RNAseq, we observed activation of the Hh pathway as well as upregulation of target genes of the Hippo-YAP pathway and activation of MYCN target genes. The functional analysis of the novel tumorigenic driver mutations in MYCN, PTPN14 and LATS1 demonstrates their relevance in BCC tumorigenesis and provides an expanded molecular understanding of BCC.

Study EGAS00001001540

A somatic reference standard for cancer genome sequencing with COLO829

Identification of somatic alterations in cancer has become feasible with the massive adoption of next generation sequencing. However, due to variability in sequencing and informatics pipelines, a common somatic reference is needed. We thus performed paired whole genome sequencing of a metastatic melanoma cell line (COLO829) and a matched lymphoblastoid line (COLO829BL) across three institutions (TGen, Illumina, Genome Sciences Centre at the British Columbia Cancer Agency). We performed a meta-analysis of all data, in combination with the originally reported analysis of these cell lines by Pleasance et al (PMID: 20016485), and report a somatic reference standard based on consensus events. DNA extractions, library preparation, sequencing, and analysis were separately performed at each site. Results were compiled with sequence data previously generated for the cell lines to identify true positive variants, collectively representing a somatic reference standard. Overall, common somatic events include point mutations, small insertion/deletions, and genes affected by concurrent copy number changes. We additionally show divergence of somatic mutations in COLO829 associated with differences in cell line lineage. We present this reference to the community as a standard for enabling interpretation and identification of somatic mutations across institutions and analytical pipelines.

Study EGAS00001001385

MYO5B mutations in Pheochromocytoma/Paraganglioma promote cancer progression

Identification of additional cancer-associated genes and secondary mutations driving the metastatic progression in pheochromocytoma and paraganglioma (PPGL) is important. We recently reported novel recurrent nonsynonymous mutations in the MYO5B gene in metastatic PPGL. Here, we explored the functional impact of these MYO5B mutations, and analyzed MYO5B expression in primary PPGL tumor cases in relation to mutation status.Immunohistochemistry and mRNA expression analysis in 30 primary PPGL tumors revealed an increased MYO5B expression in metastatic compared to non-metastatic tumor cases. In addition, subcellular localization of MYO5B protein was altered from cytoplasmatic to membranous in some metastatic tumors, and the strongest and most abnormal expression pattern was observed in a PGL case harboring a MYO5B:p.G1611S mutation. Mutation analysis of additional PPGL tumors revealed additional, and hence recurrent, mutations in the paralog MYO5A. Three missense mutations in MYO5B were functionally studied using site directed mutagenesis and stable transfection into human neuroblastoma cells (SK-N-AS) and embryonic kidney cells (HEK293 cell line). Functional analysis of the mutants verified a significant increased proliferation rate in all MYO5B mutated clones (p.L587P, p.G1611S, and p.R1641C). The two somatically derived mutations p.L587P and p.G1611S were also found to increase the migration rate, and all three mutants affected intracellular transport by increased transferrin uptake. Downstream expression analysis of MYO5B mutants demonstrated a transcriptional activation of genes involved in glucose metabolism, proliferation, homeostasis, and migration. Our study validates the functional role of novel MYO5B mutations in proliferation and migration, and suggest the MYO5-pathway to be involved the malignant progression of PPGL tumors.

Study EGAS00001003991

The spatial organization of intratumor heterogeneity andevolutionary trajectories of metastases in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) has one of the poorest survival amongst cancers. Using multi-regional sampling of nine resected HCC with different etiologies, we find that phylogenetic relationships of these sectors show diverse levels of genetic sharing spanning early to late diversification. Unlike the variegated pattern found in colorectal cancers, a large proportion of HCC display a clear Isolation-By-Distance pattern where spatially closer sectors are genetically more similar. Two resected intra-hepatic metastases showed genetic divergence occurring before and after primary tumor diversification respectively. Interestingly, metastatic tumors had much higher variability than their primary tumors, suggesting that intra-hepatic metastasis is accompanied by rapid diversification at the distant location. The presence of co-existing mutations offers the possibility of drug repositioning for HCC treatment. Taken together, these insights into intra-tumor heterogeneity allow for a comprehensive understanding of the evolutionary trajectories of HCC and suggest novel avenues for personalized therapy.

Study EGAS00001001603

Comparison of HCC cell lines and primary HCCs

There are concerns about whether cancer cell lines could faithfully represent the matched primary cancer cells. Comparison of the HCC cell lines and primary HCCs demonstrated that, during long-term in vitro culture, cell lines retain the genetic landscape of the matched primary HCCs.

Study EGAS00001001678

THE NOWAC BLOOD-BREAST CANCER COMMITEE DATA ACCESS AGREEMENT

Dac EGAC00001000743

Molecular analysis of circulating tumour cells identifies distinct profiles in chemosensitive and chemorefractory patients with small cell lung cancer

In a search for profiles that distinguish chemosensitive versus chemorefractory SCLC disease, we examined copy number aberrations (CNA) in circulating tumour cells (CTCs) from pre-treatment SCLC blood samples. A CNA classifier was generated from 88 CTCs isolated from 13 patients (training set) then verified in 18 additional patients (112 CTC samples) (testing set) and in six SCLC patient-derived CTC explants explant tumours. These data highlight the potential utility of CTCs for molecular profiling and stratification of SCLC patients in future clinical trials.

Study EGAS00001001951

Whole exome sequencing on primary retinoblastoma tissues and matching lymphocyte DNA.

Retinoblastoma is a rare childhood cancer initiated by RB1 mutation or MYCN amplification, while additional alterations may be required for tumor development. However, the view on single nucleotide variants is very limited. To better understand oncogenesis, we determined the genomic landscape of retinoblastoma. We performed exome sequencing of 71 retinoblastomas and matched blood DNA. Next, we determined the presence of single nucleotide variants, copy number alterations and viruses. Aside from RB1, recurrent gene mutations were very rare. Only a limited fraction of tumors showed BCOR (7/71, 10%) or CREBBP alterations (3/71, 4%). No evidence was found for the presence of viruses. Instead, specific somatic copy number alterations were more common, particularly in patients diagnosed at later age. Recurrent alterations of chromosomal arms often involved less than one copy, also in highly pure tumor samples, suggesting within-tumor heterogeneity. Our results show that retinoblastoma is among the least mutated cancers and signify the extreme sensitivity of the childhood retina for RB1 loss. We hypothesize that retinoblastomas arising later in retinal development benefit more from subclonal secondary alterations and therefore, these alterations are more selected for in these tumors. Targeted therapy based on these subclonal events might be insufficient for complete tumor control.

Study EGAS00001001690

Automated system for scoring hematoxylin and eosin-stained ovarian cancer sections by identifying single cells uncovered that stromal cell ratio is a significant predictor for overall survival and progression-free survival.

Concerted efforts in genomic studies have revealed profound insights in prognostic ovarian cancer subtypes. On the other hand, abundant histology slides have been generated to date, yet their uses remain very limited and largely qualitative. Our goal is to develop automated histology analysis as an alternative subtyping technology for ovarian cancer that is cost-efficient and do not rely on DNA quality. We develop an automated system for scoring hematoxylin and eosin-stained (H&E) primary tumour sections of 91 late-stage ovarian cancer to identify single cells including cancer and stromal cells. We demonstrated high accuracy of our system based on expert pathologists’ scores (cancer=97.1%, stromal=89.1%) as well as compared to immunohistochemistry scoring (correlation=0.87). Quantitative stromal cell ratio is significantly associated with poor overall survival after controlling for clinical parameters including debulking status and age (multivariate analysis p=0.0021, HR=2.54, CI=1.40-4.60) and progression-free survival (multivariate analysis p=0.022, HR=1.75, CI=1.09–2.82). We demonstrate how automated image analysis enables objective quantification of microenvironmental composition of ovarian tumours. Our analysis reveals a strong effect of the tumour microenvironment on ovarian cancer progression and highlights the potential of therapeutic interventions that target the stromal compartment or cancer-stroma signalling in the stroma-high, late-stage ovarian cancer subset.

Study EGAS00001001694

Similarity and diversity of the tumor microenvironment in multiple metastases: critical implications for overall and progression-free survival of high-grade serous ovarian cancer.

The tumor microenvironment is pivotal in influencing cancer progression and metastasis. Different cells co-exist with high spatial diversity within a patient, yet their combinatorial effects are poorly understood. We investigate the similarity of the tumor microenvironment of 192 local metastatic lesions in 61 ovarian cancer patients. An ecologically inspired measure of microenvironmental diversity derived from multiple metastasis sites is correlated with clinicopathological characteristics and prognostic outcome. We demonstrate a high accuracy of our automated analysis across multiple sites. A low level of similarity in microenvironmental composition is observed between ovary tumor and corresponding local metastases (stromal ratio r = 0.30, lymphocyte ratio r = 0.37). We identify a new measure of microenvironmental diversity derived from Shannon entropy that is highly predictive of poor overall (p = 0.002, HR = 3.18, 95% CI = 1.51-6.68) and progression-free survival (p = 0.0036, HR = 2.83, 95% CI = 1.41-5.7), independent of and stronger than clinical variables, subtype stratifications based on single cell types alone and number of sites. Although stromal influence in ovary tumors is known to have significant clinical implications, our findings reveal an even stronger impact orchestrated by diverse cell types. Quantitative histology-based measures can further enable objective selection of patients who are in urgent need of new therapeutic strategies such as combinatorial treatments targeting heterogeneous tumor microenvironment.

Study EGAS00001002065

Whole-exome sequencing of breast cancer metastasis and corresponding blood samples

This study included samples from metastatic breast cancer patients who underwent a biopsy in the context of SAFIR01 (NCT01414933), SAFIR02 (NCT02299999), SHIVA (NCT01771458) and MOSCATO (NCT01566019) prospective trials. These French multicentre trials used high throughput genome analysis on fresh frozen tumor biopsies as a therapeutic decision tool for metastatic cancer patients, with solid cancers (SHIVA, MOSCATO) or specifically with breast cancer (SAFIR01, SAFIR02). All patients gave their informed consent for translational research and genetic analyses of their somatic DNA. Overall, whole-exome sequencing for a total of 216 pairs of metastatic tumor and unmutated DNA derived from corresponding blood samples was performed using Illumina technology. The results of the analysis of the mutational profile of metastatic breast cancer was reported in Lefebvre et al. (PMID:28027327).

Study EGAS00001001695

Korean Young Age Diffuse Gastric Cancers

The incidence of diffuse gastric cancers (DGC) remains steadily high, but relatively small numbers of DGCs have been evaluated with whole exome sequencing (WES). Clinicopathological features of DGC vary according to age, but the molecular basis for these variations is unclear. To identify genomic alteration in this unique subset of gastric cancers, whole exome and SNP6 analyses were performed using frozen cancer tissue. The same analyses were also performed using blood of the same patients.

Study EGAS00001001711

Transcriptome profiling for Korean Diffuse Gastric cancers

Gastric cancers in young patients are aggressive and predominantly of diffuse histology. Also, they are more frequent in female. To identify transcriptome profile in this unique subset of gastric cancers, RNA-seq analyses were performed using frozen cancer tissue. Adjacent normal tissue of the same patients were used in differently expressed gene selection and fusion gene prediction.

Study EGAS00001001859

Genomic characterization of esophageal squamous cell carcinoma reveals critical genes underlying tumorigenesis and poor prognosis

The genetic mechanisms underlying the poor prognosis of esophageal squamous cell carcinoma (ESCC) are not well understood. In this study, we comprehensively characterized somatic mutations, copy number alterations (SCNAs)/structural variants (SVs) found in ESCC from sequencing 10 whole-genome and 57 whole-exome matched tumor-normal pairs. We identified multiple somatic mutations seen previously in known cancer pathways and identified candidate genes for ESCC including VANGL1 and MIR4707. A survival analysis based on the expression profiles of 321 ESCC individuals indicated that the somatically altered genes we found were significantly associated with ESCC poorer survival. Subsequently, we performed functional studies to validate the roles of the altered genes in tumor proliferation and metastasis.

Study EGAS00001001723

Sequencing of liver cancer cell lines

Liver cancer is a dreadful disease with limited therapeutic options. How to translate genomic findings into biomarkers and precision treatment is challenging. A large cell line-based platform was established, termed as Liver Cancer Model Repository. Heterogeneities in patients were faithfully recapitulated by 81 liver cancer cell lines at genomic and transcriptomic levels. This platform provides a rich resource and models to study cancer drivers and drug responses in liver cancer.

Study EGAS00001002237

Dissection of the molecular complexity of colorectal cancer in pre-clinical models identifies predictive signatures of sensitivity to EGFR inhibitors

Colorectal carcinoma (CRC) is a heterogeneous tumour entity, with only a fraction of the tumours responding to currently available therapies. A detailed understanding of the biology and genetics of CRC is a prerequisite for precision medicine. To address this challenge in clinical oncology, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed an integrated pre-clinical platform analysing drug sensitivity in the established matched patient-derived in vivo and in vitro models, 46 three-dimensional cell cultures and 59 xenografts, recapitulating many of the molecular features of the donor tumours. We generated extensive omics information on patient tumours and their derived models used to generate a compendium of drug sensitivity data. We could distinguish tumour molecular sub-groups responding differently to therapy and identified a novel independent predictor of sensitivity to the EGFR inhibitor cetuximab. The findings reported here will advance our understanding of CRC and its response to therapy.

Study EGAS00001001752

Identification of fusion transcripts by RNA-sequencing and Whole genome sequencing of a METABRIC patient sample

Fusion transcripts and their exact breakpoints were identified in RNA-seq data (30X) of an ER-,HER2+ breast cancer. DNA-level breakpoints associated with the fusion events were determined by analyzing whole genome sequencing (WGS) data (30X) from the same donor. The RNA-seq and WGS data for this case were generated on the Illumina platform using standard protocols. The patient (MB-0152) derives from the METABRIC cohort for which array based copy number and gene expression data were previously described (EGAS00000000083).

Study EGAS00001002475

METABRIC: Data from Batra et al (2021); DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, RRBS-based profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, the study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Study EGAS00001004327

QSEA – modelling of genome-wide DNA methylation from sequencing enrichment experiments

We describe a new method, QSEA, for analyzing methylation enrichment data. We generated a benchmark experiment consisting of MeDIP-seq enrichment data, targeted BS-seq validation data, and RNA-seq data on five pairs of tumor and adjacent normal samples of non-small cell lung cancer patients. Details of the method and its performance on the data have been described in Lienhard et al. (2016).

Study EGAS00001001822

Epigenome and transcriptome profiling of chronic lymphocytic leukemia patients

Chronic lymphocytic leukemia (CLL) is characterized by substantial clinical heterogeneity, despite relatively few genetic alterations. To provide a basis for studying epigenome deregulation in CLL, we established genome-wide chromatin accessibility maps for 88 CLL samples from 55 patients using the ATAC-seq assay, and we also performed ChIPmentation and RNA-seq profiling for ten representative samples. Based on the resulting dataset, we devised and applied a bioinformatic method that links chromatin profiles to clinical annotations. Our analysis identified sample-specific variation on top of a shared core of CLL regulatory regions. IGHV mutation status – which distinguishes the two major subtypes of CLL – was accurately predicted by the chromatin profiles, and gene regulatory networks inferred for IGHV-mutated vs. IGHV-unmutated samples identified characteristic differences between these two disease subtypes. In summary, we discovered widespread heterogeneity in the chromatin landscape of CLL, established a community resource for studying epigenome deregulation in leukemia, and demonstrated the feasibility of chromatin accessibility mapping in cancer cohorts and clinical research.

Study EGAS00001001821

Whole exome sequencing of Finnish hereditary breast cancer families

A remarkable proportion of factors causing genetic predisposition to breast cancer (BC) are unknown in non-BRCA1/2 families. Exome sequencing was performed for 13 high-risk Finnish hereditary breast and/or ovarian cancer (HBOC) families to detect variants contributing to BC susceptibility

Study EGAS00001001835

A novel TP53-KPNA3 translocation defines a de novo treatment-resistant clone in osteosarcoma

Osteosarcoma is the most common primary bone cancer. It can be cured by aggressive surgery and chemotherapy, but outcomes for metastatic or chemoresistant disease remain dismal. Cancer sequencing studies have shown that the p53 pathway is dysregulated in nearly every case, often by translocation; however, no studies of osteosarcoma evolution or intratumor heterogeneity have been done to date. We studied a patient with chemoresistant, metastatic disease over the course of three years. We performed exome sequencing on germline DNA and DNA collected from tumor at three separate timepoints. We compared variant calls and variant allele frequencies between different samples. We identified subclonal driver mutations in several different genes in the primary tumor sample and found that one particular subclone dominated subsequent tumor samples at relapse. This clone was marked by a novel TP53-KPNA3 translocation and loss of the opposite-strand wild-type TP53 allele. Future research must focus on functional significance of such clones and strategies to eliminate them.

Study EGAS00001001805

Spatial genomic heterogeneity in multiple myeloma revealed by multi- region sequencing

The bone-marrow (BM) containing skeletal system is the reservoir of Multiple Myeloma (MM), the malignant counterpart of antibody-secreting plasma cells. To gain insight into its spatial clonal architecture, we performed multi-region whole-exome sequencing of radiology-guided fine-needle aspirates from 51 MM patients. We found spatial genomic heterogeneity in a substantial proportion of cases including bi-allelic inactivation of tumor suppressors and mutations affecting cancer genes. The amount of heterogeneous mutations associated with the size of focal tumor lesions in whole body imaging consistent with focal outgrowth of highly advanced clones. In conclusion, our results support the regional initiation of transformation processes leading to high-risk disease as well as the pre-existence of key driver events in restricted areas. As such, our study does not only provide new insights in the underlying biology of MM progression but also has considerable implications for standard diagnostic approaches in the clinic.

Study EGAS00001002111

Subtype specific progression from DCIS to invasive breast cancer

Breast cancer consists of at least five main molecular “intrinsic” subtypes, which are reflected in both pre-invasive and invasive disease. Although previous studies have suggested that many of the molecular features of invasive breast cancer are established early, it is unclear what mechanisms drive progression, and whether the mechanisms of progression are dependent or independent of subtype. We have generated mRNA, miRNA and DNA copy number profiles from a total of 59 in situ lesions and 85 invasive tumors, in order to comprehensively identify those genes, signaling pathways, processes, and cell types that are involved in breast cancer progression. Our work provides evidence that there are molecular features associated with disease progression that are unique to the intrinsic subtypes. We additionally establish subtype-specific signatures that are able to identify a small proportion of pre-invasive tumors with expression profiles that resemble invasive carcinoma, indicating a higher likelihood of future disease progression.

Study EGAS00001001866

BC Cancer, part of the Provincial Health Services Authority, Technology Development Office, Data Access Committee (PHSA TDO DAC)

Dac EGAC00001003254

Genome-wide quantification of rare somatic mutations in normal human tissues using massively parallel sequencing

We present Bottleneck Sequencing System (BotSeqS), a nextgeneration sequencing method that simultaneously quantifies rare somatic point mutations across the mitochondrial and nuclear genomes. BotSeqS combines molecular barcoding with a simple dilution step immediately prior to library amplification. We use BotSeqS to show age and tissue-dependent accumulations of rare mutations and demonstrate that somatic mutational burden in normal tissues can vary by several orders of magnitude, depending on biologic and environmental factors. We further show major differences between the mutational patterns of the mitochondrial and nuclear genomes in normal tissues. Lastly, the mutation spectra of normal tissues were different from each other, but similar to those of the cancers that arose in them. This technology can provide insights into the number and nature of genetic alterations in normal tissues and can be used to address a variety of fundamental questions about the genomes of diseased tissues.

Study EGAS00001001838

Circulating Tumor DNA Analysis Detects Minimal Residual Disease and Predicts Recurrence in Patients with Stage II Colon Cancer

Detection of circulating tumor DNA (ctDNA) after resection of stage II colon cancer may identify those at the highest risk of recurrence and help inform adjuvant treatment decisions. We used massively-parallel-sequencing–based assays to evaluate the ability of ctDNA to detect minimal residual disease in 1046 plasma samples from a prospective cohort of 230 patients with resected stage II colon cancer. In patients not treated with adjuvant chemotherapy, ctDNA was detected post-operatively in 14 of 178 (7.9%) patients, 11 (79%) of whom had recurred at a median follow-up of 27 months; recurrence occurred in only 16 (9.8 %) of the 164 patients with negative ctDNA (HR, 18; 95% confidence interval [CI], 7.9 to 36; P<0.001). In patients treated with chemotherapy, the presence of ctDNA following completion of chemotherapy was also associated with an inferior recurrence-free survival (HR, 11; 95% CI, 1.8 to 68; P=0.001). ctDNA detection after stage II colon cancer resection provides direct evidence of residual disease and identifies patients at very high risk of recurrence.

Study EGAS00001001839

Whole-Exome Sequencing of Salivary Gland Mucoepidermoid Carcinoma

Mucoepidermoid carcinoma (MEC) is the most common salivary gland malignancy. To explore the genetic origins of MEC, we performed systematic genomic analyses of these tumors.Experimental DesignWhole-exome sequencing and gene copy number analyses were performed for 18 primary cancers with matched normal tissue. Fluorescence in situ hybridization (FISH) was used to determine the presence or absence of the MECT1-MAML2 translocation in 17 tumors.ResultsTP53 was the most commonly mutated gene in MEC (28%), and mutations were found only in intermediate- and high-grade tumors. Tumors with TP53 mutations had more mutations overall than tumors without TP53 mutations (p=0.006). POU6F2 was the second most frequently mutated gene, found in three low-grade MECs with the same in-frame deletion. Somatic alterations in IRAK1, MAP3K9, ITGAL, ERBB4, OTOGL, KMT2C, and OBSCN were identified in at least two of the 18 tumors sequenced. FISH analysis confirmed the presence of the MECT1-MAML2 translocation in 15 of 17 tumors (88%).ConclusionsThrough these integrated genomic analyses, MECT1-MAML2 translocation and somatic TP53 and POU6F2 mutations appear to be the main drivers of mucoepidermoid carcinoma.

Study EGAS00001002811

Cancer-Associated Mutations in Endometriosis without Cancer

Endometriosis, defined as the presence of ectopic endometrial stroma and epithelium, affects approximately 10% of reproductive-age women and can cause pelvic pain and infertility. Endometriotic lesions are considered to be benign inflammatory lesions but have cancerlike features such as local invasion and resistance to apoptosis

Study EGAS00001003576

An Isogenic Cell Line Panel for Sequence-based Screening of Targeted Anti-cancer Drugs

We describe the creation and characterization of an isogenic cell line panel representing common cancer pathways, with multiple features optimized for high-throughput screening. More than 1,800 cell lines from three normal human cells were generated using CRISPR-technologies. Surprisingly, we discovered most of these lines did not result in complete gene inactivation, despite integration of sgRNA at the desired genomic site. However, a subset of the lines harbored true, biallelic disruptions of the targeted tumor suppressor gene, yielding a final panel of 100 well-characterize lines covering 19 pathways frequently subject to loss of function in cancers. This panel included genetic markers optimized for sequence-based ratiometric assays for drug-based screening assays. To illustrate the potential utility of this panel, we developed a multiplexed high-throughput screen that identified Wee1 inhibitor MK-1775 as a selective growth inhibitor of cells with inactivation of TP53. These cell lines and screening approach should prove useful for researchers studying a variety of cellular and biochemical phenomena

Study EGAS00001005974

Multi-omics analysis of primary glioblastoma cell-lines shows recapitulation of pivotal molecular features of parental tumors

Glioblastoma (GBM) is the deadliest primary brain cancer in adults. Emerging innovative therapies hold promise for personalized cancer treatment. Improving therapeutic options depends on research relying on relevant preclinical models. In this line we have established in the setting of the GlioTeX project (GBM and experimental therapeutics), a GBM-patient derived cell line library (GBM-PDCL). Multi-OMIC approach was used to determine the molecular landscape of PDCL and the extent to which they represent GBM tumors.

Study EGAS00001001871

Same-day genomic and epigenomic diagnosis of brain tumors using realtime nanopore sequencing

Molecular classification of cancer has entered clinical routine to inform diagnosis, prognosis and treatment decisions. At the same time, new tumor entities have been identified that cannot be defined histologically. For central nervous systems tumors, the current World Health Organization classification explicitly demands molecular testing, e.g. for 1p/19q- codeletion or IDH mutations, to make an integrated histomolecular diagnosis. However, a plethora of sophisticated technologies is currently needed to assess different genomic and epigenomic alterations and turnaround times are in the range of weeks, which makes standardized and widespread implementation difficult and hinders timely decision making. Here, we explored the potential of a pocket-size nanopore sequencing device for multimodal and rapid molecular diagnostics of cancer. Low-pass whole genome sequencing was used to simultaneously generate copy number (CN) and methylation profiles from native tumor DNA in the same sequencing run. Single nucleotide variants in IDH1, IDH2, TP53, H3F3A and the TERT promoter region were identified using deep amplicon sequencing. Nanopore sequencing yielded ~0.1X genome coverage within six hours and resulting CN and epigenetic profiles correlated well with matched microarray data. Diagnostically relevant alterations, such as 1p/19q codeletion, and focal amplifications could be recapitulated. Using ad hoc random forests, we could perform supervised pan-cancer classification to distinguish gliomas, medulloblastomas and brain metastases of different primary sites. Single nucleotide variants in IDH1, IDH2 and H3F3A were identified using deep amplicon sequencing within minutes of sequencing. Detection of TP53 and TERT promoter mutations shows that sequencing of entire genes and GC-rich regions is feasible. Nanopore sequencing allows same-day detection of structural variants, point mutations and methylation profiling using a single device with negligible capital cost. It outperforms hybridization-based and current sequencing technologies with respect to time-to diagnosis and required laboratory equipment and expertise, aiming to make precision medicine possible for every cancer patient, even in resource restricted settings.

Study EGAS00001002213

Molecular profiling reclassifies adult astroblastoma into known and clinically distinct tumor entities with frequent MAPK pathway alterations

Study EGAS00001003798

The genomic landscape of cutaneous squamous cell carcinoma from immunosuppressed and immunocompetent patients reveals common drivers and a novel mutational signature associated with chronic azathioprine exposure

Cutaneous squamous cell carcinoma (cSCC) has a high tumour mutational burden (50 mutations per megabase DNA pair). Combining whole exome analyses from 40 primary cSCC tumours, comprising 20 well differentiated and 20 moderately/poorly differentiated tumours, with accompanying clinical data from a longitudinal study of immunosuppressed and immunocompetent patients, coupled with integration of independent gene expression studies has allowed the identification of commonly mutated genes and altered pathways and processes. Mutational signature analysis revealed the presence of a novel signature (signature 32), which was found exclusively in the tumours from immunosuppressed patients who had been exposed to azathioprine. Characterisation of a panel of 15 keratinocyte cancer cell lines derived from tumours from immunocompetent and immunosuppressed patients reveals that they accurately reflect the mutational signatures and genomic alterations of primary tumours and provide a resource for the validation of tumour drivers and therapeutic targets.

Study EGAS00001002612

Impact of mutational profiles on response of primary oestrogen receptor-positive breast cancers to oestrogen deprivation

Presurgical studies allow study of the relationship between mutations and response of estrogen receptor positive (ER+) breast cancer to aromatase inhibitors (AIs) but have been limited to small biopsies. Here in Phase I of this study, we perform exome sequencing on baseline, surgical core-cuts and blood from 60 patients (40 AI treated, 20 Controls). In poor responders (based on Ki67 change) we find significantly more somatic mutations than good responders. Subclones exclusive to baseline or surgical cores occur in approximately 30% of tumours. In Phase II we combine targeted sequencing on another 28 treated patients with Phase I. We find six genes frequently mutated: PIK3CA, TP53, CDH1, MLL3, ABCA13 and FLG with 71% concordance between paired cores. TP53 mutations are associated with poor response. We conclude that multiple biopsies are essential for confident mutational profiling of ER+ breast cancer and TP53 mutations are associated with resistance to oestrogen deprivation therapy.

Study EGAS00001001940

Multisample genomic analysis of solid childhood cancers using high resolution SNP-arrays, Whole Exome Sequencing and Targeted Deep Sequencing.

We mapped the prevalence of genetically distinct clones over 248 regions in 54 childhood cancers. This revealed that primary tumors can simultaneously follow up to four evolutionary patterns in different anatomic areas. The most common pattern consists of a fluctuating presence over anatomic space of subclones with very few mutations. The second most common is a surprisingly stable coexistence, over vast areas, of clones with different changes in chromosome numbers. This is contrasted by a third, less frequent, pattern where a novel clone harboring driver mutations or structural chromosome rearrangements completely replaces its progenitor. The fourth and rarest pattern is the local emergence of a myriad of clones with TP53 inactivation. Strikingly, death from disease was limited to tumors exhibiting the two latter, most dynamic patterns.

Study EGAS00001002662

Mesenchymal inflammation drives genotoxic stress in hematopoietic stem cells and predicts disease evolution in human pre-leukemia

Mesenchymal niche cells may drive tissue failure and malignant transformation in the hematopoietic system but the molecular mechanisms and their relevance to human disease remain poorly defined. We demonstrated that perturbation of mesenchymal cells in a mouse model of the preleukemic disorder Shwachman-Diamond syndrome induces mitochondrial dysfunction, oxidative stress and activation of DNA damage responses in hematopoietic stem and progenitor cells. In this study we demonstrate, through massive parallel RNA sequencing of highly purified mesenchymal cells in a range of human preleukemic syndromes, TP53-S100A8/9-TLR4 inflammatory signaling as a common driving mechanism of genotoxic stress, which could be attenuated by TLR4 blockade. S100A8/9 expression in mesenchymal cells predicted outcome in myelodysplastic syndromes, the principal human preleukemic condition, independent of known prognostic variables. Collectively, findings reveal a concept of mesenchymal niche-induced genotoxic stress in heterotypic stem and progenitor cells through inflammatory signaling as an actionable determinant of disease outcome in human preleukemia. The data further provide novel conceptual and mechanistic insights into the intimate link between inflammation and cancer.

Study EGAS00001001926

Deep single-cell RNA sequencing data for 11138 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive CRC patients. The DATA ACCESS AGREEMENT is provided at https://github.com/zhangyybio/single-T-cell-data-access. Applicants can request access to the data by directly downloading it or by sending an email to cancerpku@pku.edu.cn. The process that is used to approve an application includes verifying the institution, participants and research purposes of the application. In general this process will take about two weeks. In principal, any scientific research program complying with the laws and bioethic regulation policies of China will be approved.

T cells are central players in cancer immunotherapy1, yet some of their fundamental properties such as development and migration within tumours remain elusive. The enormous T cell receptor (TCR) repertoire, required for recognising foreign and self-antigens2,3, could serve as lineage tags to track these T cells in tumours4. Here, we obtained transcriptomes of 11,138 single T cells from 12 colorectal cancer (CRC) patients and developed STARTRAC (Single T-cell Analysis by Rna-seq and Tcr TRACking) indices to quantitatively analyse dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. While both CD8+ effector and “exhausted” T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, the majority of tumour-infiltrating Tregs showed clonal exclusivity, whereas certain Treg clones were developmentally linked to multiple TH clones. Notably, we identified two IFNG+ TH1-like clusters in tumours, the GZMK+ TEM and CXCL13+ TH1-like clusters, which were associated with distinct IFN-γ-regulating transcription factors, EOMES/RUNX3 and BHLHE40, respectively. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in tumours of microsatellite-instable (MSI) patients, which might explain their favourable response rates to immune-checkpoint blockade. Furthermore, we found IGFLR1 to be highly expressed in both CXCL13+BHLHE40+ TH1-like and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provides a powerful avenue to comprehensively dissect the T cell properties in CRC, which could shed new insights into the dynamic relationships of T cells in other cancers

Study EGAS00001002791

Deep single-cell RNA sequencing data for 12346 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive NSCLC patients

Cancer immunotherapies have shown sustained clinical responses in treating non-small cell lung cancer (NSCLC), but efficacy varies between patients and is believed to depend in part on the amount and properties of tumor infiltrating lymphocytes (TILs). To comprehensively depict and dissect the baseline landscape of the composition, lineage and functional states of TILs in lung cancer, here we generated deep single-cell RNA sequencing data for 12,346 T cells from the primary tumour, adjacent normal tissues and peripheral blood of 14 treatment-naive NSCLC patients. Combined expression and TCR-based lineage tracking revealed a significant proportion of effector T cells with common origins and similar functional states across peripheral blood and tumours pointing towards a highly migratory nature of these T cells. We also observed tumour-infiltrating CD8+ T cells undergoing extensive clonal expansion and exhaustion, with two clusters of cells exhibiting states preceding exhaustion. Survival analysis on independent datasets suggested that high ratio of pre-exhausted to exhausted T cells is associated with better prognosis of lung adenocarcinoma (LUAD). In addition, we observed further heterogeneity within the tumour regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of this group of activated tumour Tregs, which included IL1R2, correlated with poor prognosis in LUAD. The T cell clusters revealed by our single cell analyses provide a new approach for patient stratification, and the accompanying compendium of data will help the research community to gain further insight into the functional states and dynamics of T cell responses in lung cancer.

Study EGAS00001002430

Tracing the origin of disseminated tumor cells in breast cancer using single-cell sequencing

Background Single-cell micro-metastases of solid tumors often occur in the bone marrow. These disseminated tumor cells (DTCs) may resist therapy and lay dormant or progress to cause overt bone and visceral metastases. The molecular nature of DTCs remains elusive, as well as when and from where in the tumor they originate. Here, we apply single-cell sequencing to identify and trace the origin of DTCs in breast cancer. Results We sequence the genomes of 63 single cells isolated from six non-metastatic breast cancer patients. By comparing the cells DNA copy number aberration (CNA) landscapes with those of the primary tumors and lymph node metastasis, we establish that 53% of the single cells morphologically classified as tumor cells are DTCs disseminating from the observed tumor. The remaining cells represent either non-aberrant ‘normal’ cells or aberrant cells of unknown origin that have CNA landscapes discordant from the tumor. Further analyses suggest that the prevalence of aberrant cells of unknown origin is age-dependent, and that at least a subset is hematopoietic in origin. Evolutionary reconstruction analysis of bulk tumor and DTC genomes enables ordering of CNA events in molecular pseudo-time and traced the origin of the DTCs to either the main tumor clone, primary tumor subclones, or subclones in an axillary lymph node metastasis. Conclusions Single-cell sequencing of bone marrow epithelial-like cells, in parallel with intra-tumor genetic heterogeneity profiling from bulk DNA, is a powerful approach to identify and study DTCs, yielding insight into metastatic processes. A heterogeneous population of CNA-positive cells is present in the bone marrow of non-metastatic breast cancer patients, only part of which are derived from the observed tumor lineages.

Study EGAS00001002102

Fibroblast heterogeneity drives metastatic spread in breast cancer through distinct mechanisms

The major cause of death among breast cancer patients is metastatic dissemination, suggesting that understanding the mechanisms driving this process remains an unmet medical need. Here, by concomitantly using 5 markers, we identify 4 Carcinoma-Associated Fibroblasts (CAF) subsets in primary tumors and corresponding metastatic lymph nodes (LN). Two CAF subsets, called CAF-S1 and CAF-S4, correlate with the proportion of cancer cells in invaded LN and exert distinct functions in tumor cell spread. CAF-S1 promote cancer cell migration and initiates an epithelial to mesenchymal transition through a CXCL12 / TGFb positive loop. For their part, the highly contractile CAF-S4 foster cancer cell motility and invasion in 3-dimensional matrix via Notch pathway activation. Finally, CAF-S1 and CAF-S4 accumulation in LN is a new prognostic factor at diagnostic, independent of breast cancer subtypes and LN status, strengthening the validity of our findings in human breast cancer pathophysiology.

Study EGAS00001003238

Inter and intra - tumor heterogeneity in Colorectal Cancer

We used high depth massive parallel sequencing analyses to study 27 CRC (colorectal cancer) patients, with a total of 97 tumors (both primary and matched metastases) in addition to healthy control samples. We found inter-tumor concordance for coding mutations; in contrast, gene copy numbers were found to be highly discordant between primary tumors and metastases as validated by fluorescent in-situ hybridization. in a further step, we dissected a single tumor into 68 spatially defined samples and sequenced them separately. Here we identify evenly distributed coding mutations in APC and TP53 in all tumor areas, yet highly variable gene copy numbers in numerous genes.

Study EGAS00001002150

The subclonal architecture of metastatic breast cancer: Results from a prospective community-based rapid autopsy program 'CASCADE'

Background: Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. The majority of work in this area has targeted primary tumours however, and very few studies have performed comprehensive profiling of advanced disease. Evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease. Methods and Findings: Three patients with ER-positive, HER2-negative breast cancer and one patient with triple negative breast cancer underwent rapid autopsy as part of an institutional prospective community-based rapid autopsy program. Cases represented a range of management problems in breast cancer, including late relapse after early stage disease; de novo metastatic disease; discordant disease response and disease refractory to treatment. Between 5 and 12 metastatic sites were collected at autopsy together with available primary tumours and longitudinal metastatic biopsies taken during life. Samples underwent paired tumour-normal whole exome sequencing and single nucleotide polymorphism arrays. Subclonal architectures were inferred by jointly analysing all samples from each patient. Mutations were validated using high depth amplicon sequencing.Between cases, there were significant differences in mutational burden, driver mutations, mutational processes and copy number variation. Within each case, we found dramatic heterogeneity in subclonal structure from primary to metastatic disease and between metastatic sites, such that no single lesion captured the breadth of disease. Metastatic cross seeding was found in each case and treatment drove subclonal diversification. Subclones displayed parallel evolution of treatment resistance in some cases, and apparent augmentation of key oncogenic drivers as an alternative resistance mechanism. We also observed the key role of mutational processes in subclonal evolution.Limitations of this study include the potential for bias introduced by joint analysis of formalin fixed archival specimens with fresh specimens, and the difficulties in resolving subclones with whole exome sequencing. Other alterations that could define subclones such as structural variants or epigenetic modifications were not assessed. Conclusions: This study highlights the variety of mechanisms that shape the genome of metastatic breast cancer, and the value of studying advanced disease in detail. Treatment drives significant genomic heterogeneity in breast cancers which has implications for disease monitoring and treatment selection in the personalised medicine paradigm.

Study EGAS00001002153

Evolutionary Origins of Recurrent Pancreatic Cancer

Surgery is the only curative option for patients with Stage I/II pancreatic cancer, nonetheless most patients will recur after surgery and die of their disease. To identify novel opportunities for management of recurrent pancreatic cancer we performed whole exome or targeted sequencing of 10 resected primary cancers and their matched intrapancreatic recurrences or distant metastases. We identified that recurrent disease develops from a monophyletic or polyphyletic origin, and this distinction has potential clinical relevance. In all patients, treatment induced genetic bottlenecks led to a modified genetic landscape and subclonal heterogeneity for driver gene alterations in part due to intermetastatic seeding. In one patient what was believed to be recurrent disease was an independent (second) primary tumor. These findings strongly advocate for routine post-treatment sampling in the management of recurrent pancreatic cancer.

Study EGAS00001004097

miR-200-regulated CXCL12β promotes fibroblast heterogeneity and immunosuppression in ovarian cancers

High-grade serous epithelial ovarian cancer (HGSOC) has been recently subdivided into molecular subgroups, including the mesenchymal HGSOC associated with partial tumor-debulking surgery and poor patient survival. Consistent with stromal-related genes signatures defining mesenchymal HGSOC, we show here that stroma plays a key function in this HGSOC molecular subtype. We first highlight stromal heterogeneity in HGSOC by identifying 4 subsets of Carcinoma-Associated Fibroblast (CAF-S1-4). Mesenchymal HGSOC significantly accumulate the activated CAF-S1 subset, which exhibits immunosuppressive functions by promoting attraction, survival and activation of regulatory T-lymphocytes. The CXCL12β chemokine reliably characterizes mesenchymal HGSOC and specifically accumulates in the CAF-S1 subset through a regulation mediated by the miR-141/200a. CXCL12β promotes attraction of T-lymphocytes and is thus a key actor in CAF-S1-mediated immunosuppressive functions. In conclusion, we highlight here for the first time stromal heterogeneity and immunosuppression in mesenchymal HGSOC that could participate in the poor patient prognosis.

Study EGAS00001002184

CAF-S3 subset in human breast and ovarian cancers

Carcinoma-associated fibroblasts (CAF) are key players in the tumor microenvironment. By combining 6 well-known stromal markers, we identify four CAF subsets (CAF-S1 to CAF-S4) in human breast and ovarian cancers.

Study EGAS00001003344

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

Study EGAS00001006293

Comparing nodal versus bony metastatic spread using tumour phylogenies

The role of lymph node metastases in distant prostate cancer dissemination and lethality is ill defined. Patients with metastases restricted to lymph nodes have a better prognosis than those with distant metastatic spread, suggesting the possibility of distinct aetiologies. To explore this, we traced patterns of cancer dissemination using tumour phylogenies inferred from genome-wide copy-number profiling of 48 samples across 3 patients with lymph node metastatic disease and 3 patients with osseous metastatic disease. Our results show that metastatic cells in regional lymph nodes originate from evolutionary advanced extra-prostatic tumour cells rather than less advanced central tumour cell populations. In contrast, osseous metastases do not exhibit such a constrained developmental lineage, arising from either intra or extra-prostatic tumour cell populations, at early and late stages in the evolution of the primary. Collectively, this comparison suggests that lymph node metastases may not be an intermediate developmental step for distant osseous metastases, but rather represent a distinct metastatic lineage.

Study EGAS00001001801

Sequencing of pancreatic cancer primary tumors and metastases

The extent of heterogeneity of driver gene mutations present in naturally occurring metastases is largely unknown, i.e. treatment-naïve metastatic disease. To address this issue, 60x whole genome sequencing of 26 metastases from 4 patients was carried out. We found that the identical driver gene mutations were present in every metastatic lesion of each patient studied. Passenger gene mutations not known or predicted to have functional consequences accounted for all intratumoral heterogeneity. Even with respect to these passenger gene mutations, the genetic similarity among the founding cells of metastases was markedly higher than that expected for any two cells randomly taken from a normal tissue. The uniformity of driver gene mutations among metastases in the same patient has critical, encouraging implications for the success of future targeted therapies in advanced stage disease.

Study EGAS00001002186

Evolutionary analysis of pancreatic cancer and coexistent precursor lesions using whole exome sequencing data

Most adult carcinomas develop from noninvasive precursor lesions, a progression that is supported by genetic analysis. Further, the evolution of these lesions contextualizes the dynamics of advanced stage disease. We analyzed the somatic variants of co-existing pancreatic cancers and precursor lesions sampled from distinct regions of the same pancreas. After inferring evolutionary relationships, we found that the ancestral cell had initiated and clonally expanded to form one or more lesions, and that subsequent driver gene mutations eventually led to a pancreatic cancer. We also found that this multi-step progression generally spans many years. Our data suggest that independent, high-grade pancreatic lesions observed in a histologic cross section often represent a single neoplasm that can spread and colonize the ductal system, accumulating spatial and genetic divergence over time.

Study EGAS00001002778

Phylogenetic analysis of treatment-naive metastases using whole exome and genome sequencing data

Metastases are responsible for the majority of cancer related deaths and are often difficult to treat successfully. Even for metastases that occur subsequent to treatment, we do not know whether such relapsing metastases were originally heterogeneous or if the observed heterogeneity is a consequence of therapy. Given that treatment can influence evolutionary dynamics by selecting resistant clones, imposing bottlenecks on cancer cell populations, and even inducing novel somatic mutations, characterizing the prior standing variation among metastases remains an important goal to predict initial therapeutic response. To quantify the heterogeneity at clinical presentation of advanced disease, we surveyed the literature for patients in which at least two treatment-naïve metastases underwent genome/exome-wide sequencing. Across all cancer types surveyed, only 18 subjects were found to fulfill this requirement. Including previously unpublished data from two subjects, we analyzed data of 74 untreated metastases and inferred cancer phylogenies. Putative driver gene mutations were acquired at a 2-fold higher rate on the trunk of all metastases than on branches. Mutations in driver genes along the trunk were strongly enriched for predicted functional consequences. Using a stochastic mathematical model, we find that driver gene heterogeneity among metastases mostly occurs in slowly growing cancers for highly advantageous driver gene mutations. Seeding efficiency only weakly affects the heterogeneity among metastases. These results explain why functional driver gene heterogeneity is uncommon in advanced disease prior to treatment, thus providing optimism for future therapies that depend on mutations across metastases.

Study EGAS00001002777

A Unifying Paradigm for Transcriptional Heterogeneity and Squamous Features in Pancreatic Ductal Adenocarcinoma

Pancreatic cancer expression profiles largely reflect a classical or basal-like phenotype. The extent to which these profiles vary within a patient is unknown. We integrated evolutionary analysis and expression profiling in multiregion sampled metastatic pancreatic cancers, finding that squamous features are the histologic correlate of an RNA-seq defined basal-like subtype. In patients with coexisting basal/squamous and classical/glandular morphology, phylogenetic studies revealed that squamous morphology represented a subclonal population in an otherwise classical/glandular tumor. Cancers with squamous features were significantly more likely to have clonal mutations in chromatin modifiers, intercellular heterogeneity for MYC amplification, and entosis. These data provide a unifying paradigm for integrating basal-type expression profiles, squamous histology, and somatic mutations in chromatin modifier genes in the context of clonal evolution of pancreatic cancer.

Study EGAS00001003974

Application of high-throughput, high-depth, targeted single-nucleus DNA sequencing in pancreatic cancer

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimized a highly automated nuclei extraction workflow that achieved fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic ductal adenocarcinoma (PDAC) patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals novel genomic evolution patterns of PDAC as well.

Study EGAS00001006024

Comprehensive analyses of somatic TP53 mutation in tumors with variable mutant allele frequency

Somatic mutation of the tumor suppressor gene TP53 is reported in at least 50% of human malignancies. Most high-grade serous ovarian cancers (HGSC) have a mutant TP53 allele. Accurate detection of these mutants in heterogeneous tumor tissue is paramount as therapies emerge to target mutant p53. We used a Fluidigm Access Array™ System with Massively Parallel Sequencing (MPS) to analyze DNA extracted from 76 serous ovarian tumors. This dataset has been made available to researchers through the European Genome-phenome Archive (EGA). Herein, we present detailed analyses of this dataset using HaplotypeCaller and MuTect2 through the Broad Institute’s Genome Analysis Toolkit (GATK). We also present comparisons of the ability to detect TP53 mutation by MPS and Sanger sequencing in these tumors. We anticipate that this TP53 mutation dataset will be useful to researchers developing and testing new software to accurately determine high and low frequency variant alleles in heterogeneous aneuploid tumor tissue. Further, the analysis pipeline we present provides a valuable framework for determining somatic variants more broadly in tumor tissue.

Study EGAS00001002200

Drug screening and whole genome sequencing of primary cells and cell lines from ovarian cancer patients to associate genomic aberrations with in vitro drug sensitivities

Identifying drugs for cancers with poor prognosis like ovarian cancer (OC) is hindered by inappropriate cell models and the limitations of current viability assays. We developed an automated confocal microscopy-based live cell assay (DeathPro) that robustly de-convolves cell death and proliferation over time by employing a specific image processing strategy. Using DeathPro, we screened nine OC patient-derived cell lines and one normal cell line (HOSEpiC) with 22 drugs or combinations. Whole genome was performed to analyze if the cell lines harbor OC-typical genomic aberrations and to relate genomic features to drug sensitivities. Commonly observed alterations in OC patients like TP53 mutations and MYC or PI3KCA amplifications were absent in the normal cell line but frequently detected in the OC cell lines. Moreover, genome-derived homologous recombination deficiency (HRD) scores of the OC cell lines correlated with their responses to carboplatin, olaparib and DNA demethylation drugs. Overall, we propose an imaging-based platform to guide drug development and genome-drug sensitivity analysis.

Study EGAS00001002239

Pheno-Seq, linking 3D phenotypes of clonal tumor spheroids to gene expression

3D-culture systems have advanced cancer modeling by reflecting physiological characteristics of in-vivo tissues, but our understanding of functional intratumor heterogeneity including visual phenotypes and underlying gene expression is still limited. Transcriptional heterogeneity can be dissected by single-cell RNA-sequencing, but these technologies suffer from low RNA-input and fail to directly correlate gene expression with contextual cellular phenotypes. Here we present pheno-seq for integrated high-throughput imaging and transcriptomic profiling of clonal tumor spheroids derived from 3D models of breast and colorectal cancer. Specifically, we identify characteristic expression signatures that are associated with heterogeneous invasive and proliferative behavior including a rare cell subtype. Furthermore, we identify functionally relevant transcriptional regulators missed by single-cell RNA-seq, link visual phenotypes defined by heterogenous expression to inhibitor response and infer single-cell regulatory states by deconvolution. We anticipate that directly linking molecular features with patho-phenotypes of cancer cells will improve the understanding of intratumor heterogeneity and consequently prove to be useful for translational research.

Study EGAS00001002999

Heterogeneous metabolic signatures are linked to cancer cell differentiation in colorectal cancer

Colorectal cancer (CRC) is characterized by functional intratumor heterogeneity that shares many similarites with the hierarchical organization of the normal intestinal epithelium. In order to relate transcriptional subtypes to functional tumor cell heterogeneity we applied scRNA-seq to 12 patient-derived CRC spheroid cultures. We identified shared expression programs that relate to intestinal lineages and revealed metabolic signatures that are linked to cancer cell differentiation. In addition, we validated and complemented sequencing results by quantitative microscopy using live-dyes and multiplexed RNA fluorescence in situ hybridization, thereby revealing metabolic compartmentalization and potential cell-cell interactions. FInally, we demonstrate functional differences between metabolically distinct lineage subtypes that might have strong implications for future treatment strategies of CRC.

Study EGAS00001004064

This study explored and validated the clinical application of targeted NGS of circulating tumor DNA in identifying tumor-specific mutations and uncovering clinical actionable targets in a variety of solid tumors in large patient cohorts.

Cancer is a disease of genetic alterations. Comprehensive genetic diagnosis is needed to match each patient to appropriate cancer therapy. However, acquisition of representative tumor samples is invasive and often infeasible. Circulating tumor DNA (ctDNA) is a promising non-invasive biomarker for cancer mutational profiling. Here we implemented targeted next generation sequencing (NGS) with a customized pan-cancer gene panel on 605 clinical ctDNA samples in multiple cancer types. Overall, tumor-specific mutations were identified in 87% of ctDNA samples, with mutation spectra highly concordant with their matched tumor tissues. 71% of patients have at least one clinical actionable mutation with 76% of which have suggested drugs approved or in clinical trials. In particular, our study reveals a unique mutation spectrum in Chinese lung cancer patients, which could be used to guide treatment decision and monitor the appearance of drug-resistant mutations. Taken together, our study demonstrated the clinical utility of target NGS-based ctDNA mutational profiling to guide cancer treatment decision.

Study EGAS00001002251

Comparison of 133/144bp dominant phenotype vs 166bp dominant phenotype.

Background: Previous studies found that cell-free DNA (cfDNA) generated from tumors was shorter than that from healthy cells, and selecting short cfDNA could enrich for tumor cfDNA and improve its usage in early cancer diagnosis and treatment monitoring; however, the underlying mechanism of shortened tumor cfDNA was still unknown, which potentially limits its further clinical application.Results: Using targeted sequencing of cfDNA in a large cohort of solid tumor patient, sequencing reads harboring tumor-specific somatic mutations were isolated to examine the exact size distribution of tumor cfDNA. For the majority of studied cases, 166bp remained as the peak size of tumor cfDNA, with tumor cfDNA showing an increased proportion of short fragments (100-150bp). Less than 1% of cfDNA samples were found to be peaked at 134/144bp and independent of tumor cfDNA purity. Using whole-genome sequencing of cfDNA, we discovered a positive correlation between cfDNA shortening and the magnitude of chromatin inaccessibility, as measured by transcription, DNase I hypersensitivity, and histone modifications. Tumor cfDNA shortening occurred simultaneously at both 5’ and 3’ ends of the DNA wrapped around nucleosomes. Conclusions: Tumor cfDNA shortening exhibited two distinctive modes. Tumor cfDNA purity and chromatin inaccessibility were contributing factors but insufficient to trigger a global transition from 166bp dominant to 134/144bp dominant phenotype.

Study EGAS00001004135

Heterogeneous Genomic Evolution and Immune Microenvironments in Metastatic Lung Cancer

The mechanism underlying the occurrence of lung cancer metastasis to different tissues/organs remains elusive. We investigated the genomic evolution and immune microenvironments of paired primary-metastatic tumors by employing multi-region whole-exome sequencing in 179 samples of 106 tumors from 51 lung cancer patients and subsequent immunohistochemistry assays in 70 of them. Our data revealed differences in genomic landscapes, molecular determinants, evolutionary dynamics, and lymphocyte infiltration among different metastatic sites. We demonstrated commonly late arising of metastatic seeding of lung cancer with quantitative evidence. Most distant metastases originated from independent origins of earlier lymph node spreads. Immune-heterogeneity and -homogeneity were primarily driven by arm-level and focal copy number events in primary tumors, respectively. These findings implied the combinatorial role of multiple factors in shaping patterns of dissemination and advanced the clinical evaluation and intervention of lung cancer metastasis.

Study EGAS00001004228

CancerLocator: Non-Invasive Cancer Diagnosis and Tissue-of-Origin Prediction Using Methylation Profiles of Cell-Free DNA

Background: The detection and characterization of cell-free DNA in plasma is one of the most promising new areas in cancer diagnosis. Liquid biopsy, unlike traditional tissue biopsy, has the potential to diagnose a variety of different malignancies. Results: Here we propose a probabilistic method, CancerLocator, which exploits the diagnostic potential of cell-free DNA by determining not only the presence but also the location of tumors. CancerLocator simultaneously infers the proportions and the tissue-of-origin of tumor-derived cell-free DNA in a blood sample using genome-wide DNA methylation data. We comprehensively evaluate CancerLocator with simulations and real data, and compare its performance with that of two established multi-class classification methods. We show that the predicted tumor burdens are highly consistent with the true values. In addition, when the proportion of tumor-derived DNAs in the cell-free DNAs is low, the two popular machine learning methods completely fail for cancer diagnosis, while CancerLocator successfully overcomes the challenge. CancerLocator also achieves promising results on patient plasma samples, despite the fact that the DNA methylation data from these samples has very low sequencing coverage. Conclusions: CfDNA methylation may be developed as an important approach for non-invasive early cancer diagnosis.

Study EGAS00001002211

BGISEQ-500 Cancer Dataset - WGS tumour/normal pairs

This study contains whole genome sequencing data for 3 matched tumour / normal pairs from Mesothelioma patients. The sequencing is 50bp paired-end generated using the BGISEQ-500 and is provided as 6 aligned and duplicate-marked BAM files.

Study EGAS00001002298

The Ovarian Cancer Association Consortium OncoArray genome-wide association study

A custom Illumina genotyping array, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis as described in PMID:27697780.The array was used to genotype ovarian cancer cases and controls from multiple studies participating in the Ovarian Cancer Association Consortium, of which 54,747 passes QC.An initial principal components analysis was carried out using a method described in PMID:27697780 to assign subject to intercontinental ancestry groups: “European”, “African”, “Asian” and "Other". Within ancestry principal components analysis was carried using a set of ~33k unlinked markers. Subsequent analyses were adjusted for the ancestry specific PCs: 9 Eur, 10 Asian, 1 African and 6 Mixed. These PCs are included in the imputed dataset phenotypes file

Study EGAS00001002305

whole genome sequence data of multifocal hepatocellular carcinoma

Background & Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics pro?ling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual speci?city protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identi?cation of novel biomarkers and therapeutic targets in HCC.

Study EGAS00001002338

The aim of this study was to identify underlying hub genes and dysregulated pathways associated with the development of HCC using bioinformatics analysis.

In the present study, we identified differentially expressed protein-coding genes from RNA transcriptional profiling performed on 11 paired cancer tissues and adjacent non-cancerous tissues. Then, we conducted GO, KEGG, PPI network and centralities analyses to study and identify changes in pathways and hub genes. The aim of this study was to improve understanding of HCC carcinogenesis by providing information concerning the genetic changes that occur during disease progression and to uncover the expression of biomarkers with potential use for clinical diagnosis, treatment, and monitoring of disease progression.

Study EGAS00001002526

TCR β-chain repertoire characterization of regulatory and conventional T cells in peripheral blood from breast cancer patients and healthy individuals.

Identification of tumor-specific effects on peripheral TCRβ repertoire formation in humans, investigation of the clonal origin of regulatory T cells in breast cancer patients and impact analysis of the tumor-specific conversion of conventional T cells into induced regulatory T cells on the peripheral Treg repertoire in humans.

Study EGAS00001002699

Type

Repository

Study type

Dataset types

Dataset technologies

Phenotypes

Duo

Domain

6137 results for "canc*"

in 46.86 milliseconds.