4055 results for "*oma"

in 91.67 milliseconds.

• Alternative splicing isoforms in patient-derived hepatocellular carcinoma cells

  The aim of this study was to obtain a landscape of alternative splicing (AS) in human hepatocellular carcinoma (HCC). We applied hybrid sequencing by integrating Pacbio long-read and Illumina short-read RNA sequencing data from 8 patient-derived HCC cell lines, and used an immortalized hepatocyte cell line as reference. This study revealed features of annotated and novel splice variants in HCC cells, including transcript length, AS type and protein changes. We also identified novel isoforms and tumor specific isoforms, which may represent potential biomarkers and targets for therapeutic development.

  Study EGAS00001002697

• Coverage bias and sensitivity of variant calling for four whole-genome sequencing technologies

  Massively parallel sequencing has revolutionized research in cancer genetics and genomics and enhanced our understanding of natural human genetic variation. Recently, Lam et al. have performed a detailed comparison of two next-generation sequencing technologies with respect to their sensitivity to call single nucleotide variants (SNV) and indels. Here, we sequenced two tumor/normal pairs obtained from two paedriatic medulloblastoma patients with Life Technologies’ SOLiD 4 and 5500xl SOLiD, Illumina’s HiSeq2000, and Complete Genomics’ technology. We then compared their ability to call SNVs with high confidence. As gold standard for SNV calling, we used genotypes determined by an Affymetrix SNP array. Additionally, we performed a detailed analysis of how evenly each technology covers the genome and how the reads are distributed across functional genomic regions. Finally, we studied how a combination of data from different technologies might help to overcome the limitations in SNV calling by any of the four technologies alone.

  Study EGAS00001000274

• Therapeutic Targeting of Ependymoma as Informed by Oncogenic Enhancer Profiling

  Genomic sequencing has driven precision-based oncology therapy; however, genetic drivers remain unknown or non-targetable for many malignancies, demanding alternative approaches to identify therapeutic leads. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Here, we mapped active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts. Enhancer regions revealed novel oncogenes, molecular targets, and pathways, which when subjected to small molecule inhibitor or shRNA treatment, increased survival and slowed proliferation in mouse and neurosphere patient-derived models of ependymomas. This study represents one of the largest enhancer mapping study of any cancer type to date, and provides a framework for target and drug discovery for other cancers recalcitrant to therapeutic development because of their lack of known genetic drivers.

  Study EGAS00001002696

• Mutational patterns and regulatory networks in epigenetic subgroups of meningioma (H033)

  DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign-1, 2, 3, intermediate-A, B and malignant. Here, we set out to identify subgroup-specific mutational patterns and pathway regulation. Whole-genome-sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic and 22 radiation-induced (Mrad) meningiomas. RNA sequencing was added for 18 cases and chromatin-immunoprecipitation for the enhancer mark H3K27ac followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural alterations were found to contribute to the spectrum of mechanisms inactivating NF2 in sporadic meningioma similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3 was detected in both sporadic meningioma and Mrad, but distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. In general, the malignant MC presented a significantly higher exposure to AC3 and higher genomic instability beyond the mutational load than the other MCs. Pathway analysis of the ChIP-seq clusters revealed that the FOXM1 is most differentially activated in high-grade cases, along with super enhancer recruitment near HOX genes and respective upregulation of expression. This data further elucidate the biological mechanisms differentiating meningiomas of different WHO grade and epigenetic subgroups and suggest leveraging the genomic instability as novel therapeutic targets.

  Study EGAS00001003481

• Genomic analysis of seminomas

  This study looks at the genomic landscape of seminomas, based on exome sequencing data and SNP array data. Eight seminoma tumors and their matched normal DNAs were exome sequenced on Illumina platforms. The same samples were also processed on Affymetrix SNP 6.0 arrays.

  Study EGAS00001000943

• Genomic analyses of gynecologic carcinosarcomas reveal frequent mutations in chromatin remodeling genes

  Malignant Mixed Müllerian Tumors, also known as carcinosarcomas, are rare tumors of gynecological origin. Here we perform whole exome analyses of 22 tumors using massively parallel sequencing to determine the mutational landscape of this tumor type. On average, we identify 43 mutations per tumor, excluding four cases with a mutator phenotype that harbored inactivating mutations in mismatch repair genes. In addition to mutations in TP53 and KRAS, we identify genetic alterations in chromatin remodeling genes, ARID1A and ARID1B, in histone methyltransferase MLL3, in histone deacetylase modifier SPOP, and in chromatin assembly factor BAZ1A, in nearly two thirds of cases. Alterations in genes with potential clinical utility are observed in more than three quarters of the cases and included members of the PI 3-kinase and homologous DNA repair pathways. These findings highlight the importance of the dysregulation of chromatin remodeling in carcinosarcoma tumorigenesis and suggest new avenues for personalized therapy.

  Study EGAS00001000941

• Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma

  Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analyzed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined 4 CCA clusters-fluke-positive CCAs (clusters 1/2) are enriched in ERBB2 amplifications and TP53 mutations; conversely, fluke-negative CCAs (clusters 3/4) exhibit high copy-number alterations and PD-1/PD-L2 expression, or epigenetic mutations (IDH1/2, BAP1) and FGFR/PRKA-related gene rearrangements. Whole-genome analysis highlighted FGFR2 3' untranslated region deletion as a mechanism of FGFR2 upregulation. Integration of noncoding promoter mutations with protein-DNA binding profiles demonstrates pervasive modulation of H3K27me3-associated sites in CCA. Clusters 1 and 4 exhibit distinct DNA hypermethylation patterns targeting either CpG islands or shores-mutation signature and subclonality analysis suggests that these reflect different mutational pathways. Our results exemplify how genetics, epigenetics, and environmental carcinogens can interplay across different geographies to generate distinct molecular subtypes of cancer.SIGNIFICANCE: Integrated whole-genome and epigenomic analysis of CCA on an international scale identifies new CCA driver genes, noncoding promoter mutations, and structural variants. CCA molecular landscapes differ radically by etiology, underscoring how distinct cancer subtypes in the same organ may arise through different extrinsic and intrinsic carcinogenic processes.

  Study EGAS00001001653

• Non-Mendalian inheritance of extrachromosal DNA elements can drive disease evolution in glioblastoma

  Genomic heterogeneity of glioblastoma (GBM) is suspected to contribute to the poor response to therapy of this disease. We compared molecular characteristics between primary GBM, neurospheres and orthotopic xenograft models derived from the same parental tumor. Driver alterations were in majority propagated from tumor to model systems. Extrachromosomal amplifications of MET, a proto-oncogene coding for a receptor tyrosine kinase, were detected in three primary GBM, largely discarded in neurospheres cultures, but resurfaced in xenografts. The clonal dynamics inferred by somatic single nucleotide variants (sSNVs) in MET-amplified samples diverged from the pattern delineated by the MET amplification event suggesting that the MET event and sSNVs were inherited in different manners. Our analysis shows that extrachromosomal elements are able to drive tumor progression.

  Study EGAS00001001878

• Pharmacogenomic landscape of patient-derived cells informs precision oncology therapy

  Outcomes of anticancer therapy vary dramatically among patients, which may be caused by the specific molecular alterations in each patient’s tumor. Precision oncology aims to apply optimal therapies for each tumor based on its molecular characteristics. We have established a resource reporting the genomic and transcriptomic profiles of 462 patient tumor-derived cells (PDCs) across 14 cancer types, together with responses to 60 targeted agents. Compared with long-term cultured cancer cell lines, PDCs better recapitulate the molecular profiles of the original tumors. Among other unreported associations, we identify molecular factors inducing resistance to EGFR inhibitors in glioblastoma, and we suggest repurposing ibrutinib (currently used in hematological malignancies) for EGFR-specific therapy in gliomas. Lastly, using a retrospective clinical study, we find that PDC-derived sensitivities can be used to predict patient responses.

  Study EGAS00001002515

• Spatiotemporal genomic architecture informs precision oncology in glioblsatoma

  We characterized genomic and expression profiles across 127 multi-sector or longitudinal specimens from 52 glioma patients. We showed that, while samples from nearby geographical location shared similar genomic alterations (clonal or subclonal) and expression signatures, multifocal tumors were seeded from individual clones.

  Study EGAS00001001880

• Preferential infiltration of distinct Vγ9δ2 T cells into Glioblastoma multiforme

  Glioblastoma multiforme (GBM), the most fatal CNS cancer, is highly infiltrative as a result of phenotypic alteration of cancer cells mainly driven by a mutual interaction among cancer cells, tumor-associated macrophages (TAMs) and other stromal cells. This work aims at characterizing GBM-infiltrating γδ T cells that may regulate the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor-infiltrating and blood circulating γδ T cells from 4 patients were analyzed by NGS-based TCR sequencing. RNA gene expression profiling and immunostaining of cancer tissues were performed as well.

  Study EGAS00001002790

• A Protein Neddylation Inhibitor MLN4924 Suppresses Patient-Derived Glioblastoma Cells via Inhibition of ERK and AKT Signaling

  Glioblastoma is a highly aggressive and lethal brain tumor, with limited treatment options. Abnormal activation of the neddylation pathway is observed in glioblastoma, and the NEDD8-activating enzyme (NAE) inhibitor, MLN4924 was previously shown to be effective in glioblastoma cell line models. However, its effect has not been tested in patient-derived glioblastoma stem cells. Immunohistochemistry and public data were analyzed to determine whether NEDD8 pathway proteins are important in glioblastoma development and patient survival. The effects of MLN4924 were evaluated in 4 glioblastoma cell lines and 15 patient-derived glioblastoma stem cells using high content analysis. Immunoblot analysis assessed the NEDD8 pathway proteins and biomarkers. Whole transcriptome analysis identified pathways associated with MLN4924 sensitivity. NAE1 and UBA3 levels increased in glioblastoma patients; high NEDD8 levels were associated with poor clinical outcomes. Glioblastoma cell lines and patient-derived stem cells were highly susceptible to MLN4924, while normal human astrocytes were resistant. Genomic analyses indicated that MLN4924 sensitive cells exhibited enrichment of ERK and AKT signaling. We verified that MLN4924 inhibits ERK and AKT phosphorylation in MLN4924 sensitive cells. Our findings suggest that patient-derived glioblastoma stem cells in the context of ERK and AKT activation are sensitive and highly regulated by neddylation inhibition.

  Study EGAS00001004018

• Mutation-specific non-canonical pathway of PTEN as a distinct therapeutic target for glioblastoma

  PTEN is one of most frequently altered tumor suppressor genes in malignant tumors. The dominant negative effect of PTEN alteration suggests that aberrant function of PTEN mutation might be more disastrous than deletion, the most frequent genomic event in glioblastoma (GBM). This study aimed to understand the functional properties of various PTEN missense mutations and to investigate their clinical relevance. The genomic landscape of PTEN alteration was analyzed using the Samsung Medical Center GBM cohort and validated via The Cancer Genome Atlas dataset. Several hot-spot mutations were identified, and their subcellular distributions and phenotypes were evaluated. We established a library of cancer cell lines that over-express these mutant proteins using the U87MG and patient-derived cell models lacking functional PTEN. PTEN mutations were categorized into two major subsets: missense mutations in the phosphatase domain and truncal mutations in the C2 domain. We determined the subcellular compartmentalization of four mutant proteins (H93Y, C124S, R130Q, and R173C) from the former group and found that they had distinct localizations; those associated with invasive phenotypes (‘edge mutations’) localized to the cell periphery, while the R173C mutant localized to the nucleus. Invasive phenotypes derived from edge substitutions were unaffected by an anti-PI3K/Akt agent but were disrupted by microtubule inhibitors. PTEN mutations exhibit distinct functional properties with regard to their subcellular localization. Further, some missense mutations (‘edge mutations’) in the phosphatase domain caused enhanced invasiveness associated with dysfunctional cytoskeletal assembly, thus suggesting it to be a potent therapeutic target.

  Study EGAS00001004753

• 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

• The Berlin (BLN) panel of glioblastoma cell lines: RNAseq of human gliomasphere cell lines and matched parental tumors

  We generated a panel of nine glioma stem-like cell lines under serum-free cell culture conditions from human glioblastoma samples. In addition to clinical annotation, MR imaging and histological characterization of the parental tumor, as well as xenograft studies of tumorigenecity, transcriptomic profiling of tumor/cell line pairs by RNA sequencing has been performed.

  Study EGAS00001001167

• Whole exome sequencing of cell-free DNA reveals temporo-spatial heterogeneity and identifies treatment-resistant clones in neuroblastoma

  Background: Neuroblastoma (NB) displays important clinical and genetic heterogeneity, with emergence of new mutations at tumor progression. Patients and Methods: To study clonal evolution during treatment and follow-up, an innovative method based on circulating cell-free DNA (cfDNA) analysis by whole exome sequencing (WES) paired with target sequencing was realized in sequential liquid biopsy samples of 19 NB patients. Results: WES of the primary tumor and cfDNA at diagnosis showed overlap of Single Nucleotide Variants (SNVs) and Copy Number Alterations (CNAs), with 41% and 93% of all detected alterations common to the primary NB and cfDNA. CfDNA WES at a second time point indicated a mean of 22 new SNVs for patients with progressive disease. Relapse-specific alterations included genes of the MAPK pathway and targeted the protein kinase A signaling pathway. Deep coverage target sequencing of intermediate timepoints during treatment and follow-up identified distinct subclones. For 17 seemingly relapse-specific SNVs detected by sequencing detected these alterations in minor subclones, with relapse-emerging SNVs targeting genes of neuritogenesis and cell cycle. Furthermore a persisting, resistant clone with concomitant disappearance of other clones was identified by a mutation in the ubiquitin protein ligase HERC2. Conclusion: Modelisation of mutated allele fractions in cfDNA indicated distinct patterns of clonal evolution, with either a minor, treatment-resistant clone expanding to a major clone at relapse, or minor clones collaborating towards tumor progression. Identification of treatment-resistant clones will enable development of more efficient treatment strategies.cfDNA WES at relapse but not tumor or cfDNA WES at diagnosis, deep coverage target.

  Study EGAS00001002705

• Enhanced reduced representation bisulfite sequencing (eRRBS) on 45 multiple myeloma samples and 3 normal plasma cell

  Enhanced reduced representation bisulfite sequencing (eRRBS) on 45 multiple myeloma samples and 3 normal plasma cell. Libraries were sequenced on a HiSeqTM4000 Illumina machine using 75bp paired-end reads

  Study EGAS00001004348

• 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

• Neuroblastoma relapse trio series from the AMC

  The majority of neuroblastoma patients have tumors that initially respond to chemotherapy, but a large proportion of patients will experience therapy-resistant relapses. The molecular basis of this aggressive phenotype is unknown. Whole genome sequencing of 23 paired diagnostic and relapsed neuroblastomas showed clonal evolution from the diagnostic tumor with a median of 29 somatic mutations unique to the relapse sample. Eighteen of the 23 relapse tumors (78%) showed RAS-MAPK pathway mutations. Seven events were detected only in the relapse tumor while the others showed clonal enrichment. In neuroblastoma cell lines we also detected a high frequency of activating mutations in the RAS-MAPK pathway (11/18, 61%) and these lesions predicted for sensitivity to MEK inhibition in vitro and in vivo. Our findings provide a rationale for genetic characterization of relapse neuroblastoma and show that RAS-MAPK pathway mutations may function as a biomarker for new therapeutic approaches to refractory disease.

  Study EGAS00001001183

• Human melanoma samples with and without resistance to BRAF inhibitor therapy

  Human melanoma samples were collected pre, on, and progression on BRAF inhibitor therapy. RNA was extracted and run on RNA-seq. This has provided insights into different categories of BRAF inhibitor resistance mechanisms.

  Study EGAS00001000992

• Whole Genome sequencing of adult T-cell leukemia/lymphoma

  Whole Genome sequencing of a single adult T-cell leukemia/lymphoma case

  Study EGAS00001001210

• Anaplastic oligodendroglioma exome and RNA sequencing data

  Anaplastic oligodendrogliomas (AOs) are rare primary brain tumors which are generally incurable, with heterogeneous prognosis and few treatment targets identified. Most oligodendrogliomas have chromosome 1p/19q co-deletion and IDH mutation. We analyzed 51 AOs by whole-exome sequencing, identifying previously reported frequent somatic mutations in CIC and FUBP1. We also identified recurrent mutations in TCF12 and in an additional series of 83 AO. Overall 7.5% of AO are mutated for TCF12, which encodes an oligodendrocyte-related transcription factor. 80% of TCF12 mutations identified were in either the bHLH domain, which is important for TCF12 function as a transcription factor, or were frame shift mutations leading to TCF12 truncated for this domain. We show that these mutations compromise TCF12 transcriptional activity and are associated with a more aggressive tumor type. Our analysis provides further insights into the unique and shared pathways driving AO.

  Study EGAS00001001209

• 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

• Tracking the genomic evolution of esophageal adenocarcinoma through neoadjuvant chemotherapy

  Esophageal adenocarcinomas (EACs) are associated with dismal prognosis. Deciphering the evolutionary histories of this disease may shed light on therapeutically tractable targets and reveal changing mutational processes during the disease course and following neoadjuvant chemotherapy (NAC). We exome sequenced 40 tumor regions from 8 patients with operable EACs both before and after platinum-containing NAC. This revealed the evolutionary genomic landscape of EACs with the presence of heterogeneous driver mutations, parallel evolution, early genome doubling events and an association between high intratumor heterogeneity and poor response to NAC. Multi-region sequencing demonstrated a significant reduction in T>G mutations within a CTT context when comparing early and late mutational processes and the presence of a platinum signature with enrichment of C>A mutations within a CpC context following NAC. EACs are characterized by the presence of amplifications containing targetable oncogenes present before and after chemotherapy that may provide a rationale for future therapeutic approaches.

  Study EGAS00001001254

• 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

• Integrated molecular analysis of adult T-cell leukemia/lymphoma

  This study is an integrated molecular study of adult T-cell leukemia/lymphoma which includes whole-exome (n = 81), whole-genome (n = 48), and transcriptome sequencing data (n = 57) as well as methylation (n = 109) and SNP array data (n = 426)

  Study EGAS00001001296

• 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

• Recurrent somatic mutations in POLR2A define a distinct subset of meningiomas

  RNA polymerase II mediates the transcription of all protein-coding genes in eukaryotic cells, a process that is fundamental to life. Genomic mutations in this enzyme have not been previously linked to any pathology in humans, a testament to its indispensable role in cell biology. Based on a combination of next-generation genomic analyses of 775 meningiomas, we report that recurrent somatic p.Gln403Lys or p.Leu438_His439del mutations in POLR2A, which encodes the catalytic subunit of RNA polymerase II1, are sufficient to hijack this essential enzyme and drive neoplasia. POLR2A mutant tumors reveal dysregulation of key meningeal identity genes, including WNT6 and ZIC1/ZIC42,3. In addition to POLR2A, NF2, SMARCB1, TRAF7, KLF4, AKT1, PIK3CA and SMO4-8, we also report somatic mutations in AKT3, PIK3R1, PRKAR1A, and SUFU in meningiomas. Our results identify a role for essential transcriptional machinery in driving tumorigenesis and define mutually exclusive meningioma subgroups with distinct clinical and pathological features.

  Study EGAS00001001916

• Reliable assessment of telomere maintenance mechanisms in neuroblastoma

  Telomere maintenance mechanisms (TMM) are a hallmark of high-risk neuroblastoma, and are conferred by activation of telomerase or alternative lengthening of telomeres (ALT). However, detection of TMM is not yet part of the clinical routine, and consensus on TMM detection, especially on ALT assessment, remains to be achieved. We here propose a workflow to reliably detect TMM in neuroblastoma. We show that unambiguous classification is feasible following a stepwise approach that determines both, activation of telomerase and ALT. The workflow proposed in this study can be used in clinical routine and provides a framework to systematically and reliably determine telomere maintenance mechanisms for risk stratification and treatment allocation of neuroblastoma patients.

  Study EGAS00001006510

• 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

• Independent development of lymphoid and histiocytic malignancies from a shared early precursor

  Recent reports have described lineage flexibility in patients with recurrent hematologic malignancies. We encountered a boy with T-lineage acute lymphoblastic leukemia, non-Langerhans cell histiocytosis and histiocytic sarcoma. All tumors showed identical T-cell receptor rearrangements, but SNP arrays and whole exome sequencing revealed unique aberrations in each tumor. From the genetic data we conclude that the tumors were derived from a precursor cell with lineage plasticity for both lymphoid and myeloid development. From this precursor cell, the tumors evolved with independent genetic progression in a non-linear manner.

  Study EGAS00001001353

• Genomes of Relapsing Neuroblastoma

  Neuroblastoma is a malignancy of the developing sympathetic nervous system that is often lethal when relapse occurs, but the molecular mechanisms behind this process are poorly defined. We here used whole-exome sequencing, mRNA expression, array CGH and DNA methylation analysis to holistically characterize 16 paired samples from neuroblastoma patients at diagnosis and relapse. The mutational burden significantly increased in relapsing tumors accompanied by altered mutational signatures and reduced subclonal heterogeneity. Global allele frequencies at relapse indicated clonal mutation selection during disease progression. Promoter methylation patterns were consistent over disease course and patient-specific. Recurrent alterations at relapse included mutations in the putative CHD5 neuroblastoma tumor suppressor, chromosome 9p losses, DOCK8 mutations, inactivating mutations in PTPN14 and a relapse-specific activity pattern for the PTPN14 target gene, YAP. Recurrent new mutations in HRAS, KRAS, DOCK8 and genes mediating cell-cell interaction in 13/16 relapse tumors also indicate disturbances in signaling pathways mediating mesenchymal transition. Our data shed first light on genetic alteration frequency, identity and evolution in neuroblastoma.

  Study EGAS00001001387

• Carcinoma of the oral tongue (OTSCC) genomic landscape characterisation

  Carcinoma of the oral tongue (OTSCC) is the most common malignancy of the oral cavity, characterized by frequent recurrence and poor survival. The last three decades has witnessed a change in the OTSCC epidemiological profile, with increasing incidence in younger patients, females and never-smokers. Here, we sought to characterize the OTSCC genomic landscape and to determine factors that may delineate the genetic basis of this disease, inform prognosis and identify targets for therapeutic intervention.

  Study EGAS00001001329

• 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

• Exome Sequencing to Define the Landscape of Plasma Cells in Systemic Light chain Amyloidosis

  Systemic light chain amyloidosis (AL) is characterized by the deposition of immunoglobulin light chains as amyloid fibrils in different organs, where they form toxic protein aggregates. The underlying disease is a plasma cell disorder, but limited whole exome data are available. We report the findings of an exome sequencing study in AL to define a plasma cell signature and compare this to monoclonal gammopathy of undefined significance (MGUS) and myeloma (MM). Twenty-four samples from unselected newly diagnosed untreated AL patients were analysed. CD138+ cells were isolated from bone marrow cells using MACSorting (Miltenyi Biotech, Bisley, UK). Cells were lysed in RLT+ buffer and DNA/RNA extracted using the AllPrep kit (Qiagen, Manchester, UK). Peripheral blood was isolated from patients, white blood cells purified by Ficoll-Pacque and DNA extracted using the QIAamp DNA mini kit (Qiagen). There were a median of 47 acquired nonsynonymous variants (range 4-253) per patient. Thirty percent of patients yielded a mutation in a myeloma driver gene. These mutations were either clonal or sub-clonal. There was evidence of MAPK activation with NRAS mutations, NF-κB activation and DNA repair pathway alterations. Exome data were used to determine the cytogenetic groups of AL samples and identified hyperdiploidy in 33% of cases and t(11;14) in 33% of cases. All translocations were generated via class switch recombination and breakpoints on chromosome 11 fell into the expected regions. A novel translocation, unlikely to be driving disease progression was identified involving RCC1. Finally, further analysis of the APCS gene encoding for the SAP protein did not show any biased genotype in the AL population in comparison to the general population. We conclude that exome sequencing identifies a genetic make-up of AL that is similar to other plasma cell disorders such as MM and MGUS. This not only includes copy number abnormalities and translocations but also nonsynonymous mutations.

  Study EGAS00001001418

• 33 patients with Monoclonal Gammopathy of Undetermined Significance (MGUS) had whole exome sequencing performed. Cells were sorted by FACSAria using CD138, CD19 and CD56 to obtain a pure abnormal plasma cell population. Generated somatic variants were compared to a previous study of 463 multiple myeloma patients.

  Monoclonal gammopathy of undetermined significance (MGUS) is a premalignant precursor of multiple myeloma (MM) with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS genes mutations, little is known about molecular mechanism of malignant transformation. We have performed whole exome sequencing together with SNP array analysis in 33 flow-cytometry separated abnormal PC samples of MGUS patients to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations (NS-SNVs) and copy number alterations (CNAs) were present in 97.0% and in 63.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in MGUS compared to MM (p<10-4) and we have identified 6 myeloma significantly mutated genes which are KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB in the MGUS dataset. We also found a positive correlation with increasing chromosome changes and somatic mutations. IGH translocations were present in 27.3% of cases comprising t(4;14), t(11;14), t(14;16) or t(14;20) and were in a similar frequency to MM, which corresponded with primary lesion hypothesis. Data from this study showed MGUS is a genetically comprehensive disease, however overall genetic instability is significantly lower compared to MM.

  Study EGAS00001001658

• Genome-wide identification of distinct miRNA-mRNA target regulation pairs in Non-Hodgkin lymphomas: a report from the ICGC MMML-Seq consortium

  We obtained miRNA profiles by miRnome sequencing from ICGC MMML-Seq patients diagnosed with Burkitt lymphoma, diffuse large B-cell lymphoma and follicular lymphoma and provide evidence of subtype-specific miRNA expression differences. We describe differentially expressed, mutated, edited and not previously annotated miRNAs. Addditionally, by performing argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) experiments, we obtained a set of biochemically-validated miRNA binding sites and identified miRNA-mRNA interaction pairs with a negative correlation in patient RNASeq data.

  Study EGAS00001001394

• 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

• Whole-genome bisulfite sequencing analysis of low-grade astrocytomas in the ICGC PedBrain Tumor Project

  Pilocytic Astrocytoma (PA) is the most common pediatric brain tumor. While genome and transcriptome landscapes are well-studied, data of the complete methylome, tumor cell composition and immune infiltration are scarce. We generated whole genome bisulfite sequence data (WGBS) of nine PAs to find evidence for a link of focal methylation differences and differential gene expression to immune infiltration.

  Study EGAS00001004019

• Genome-wide expression profiles to compare array-based, RNAseq and NanoString technologies in patients with resected pancreatic ductal carcinoma.

  Using formalin-fixed paraffin embedded (FFPE) tissues collected prospectively in patients with resected pancreatic ductal carcinoma (PDC) we generated genome-wide expression profiles to compare array-based, RNAseq and NanoString technologies. Samples included in the study were collected 17 patients. After careful selection by the pathologist of 17 malignant cells-enriched, 7 stroma-enriched, and 6 normal (normal pancreas in 4 cases and normal duodenum in 2 cases) regions, sections were macrodissected, and total RNA extracted. Two additional FFPE-frozen matched-paired PDAC from our local Biobank were included. Our study represents the most extensive work comparing three technologies to generate gene expression profiles from FFPE samples, and provides technical guidance for future biomarker studies in samples collected prospectively in the context of clinical trials.

  Study EGAS00001002501

• Whole exome sequencing of an invasive diffuse intrinsic pontine glioma sampled from different sites

  The lateral ventricle subventricular zone (SVZ) is a frequent site of high-grade glioma (HGG) spread and tumor invasion of the SVZ is a predictor of a worse clinical prognosis. A range of HGG types invade the stem cell niche of the SVZ, including both adult glioblastoma and pediatric high-grade gliomas such as diffuse intrinsic pontine glioma (DIPG). The cellular and molecular mechanisms mediating this frequent invasion of the SVZ are poorly understood. Here we demonstrate that neural precursor cell (NPC):glioma cell communication underpins the propensity of glioma to colonize the SVZ. SVZ NPCs secrete chemoattractant signals toward which glioma cells home. Biochemical and proteomic analyses of factors secreted by SVZ NPCs revealed a discrete list of candidate proteins. Necessity and sufficiency testing implicated the neurite outgrowth-promoting factor pleiotrophin, along with three required binding partners (secreted protein acidic and rich in cysteine (SPARC), SPARC-like protein 1, and heat shock protein 90B), as key mediators of this chemoattractant effect. Pleiotrophin protein expression is strongly enriched in the SVZ, and knockdown of pleiotrophin expression starkly reduced glioma invasion of the SVZ in the adult murine brain. Pleiotrophin, in complex with the three binding partners, activated the Rho/ROCK pathway in DIPG cells, and inhibition of ROCK resulted in decreased DIPG invasion toward SVZ neural precursor cell-secreted factors. These findings demonstrate a pathogenic role for neural precursor cell – glioma interactions and potential therapeutic targets to limit glioma invasion.

  Study EGAS00001002326

• Sequencing of paediatric High Grade Gliomas and DIPG

  Paediatric high grade glioma (pHGG) and diffuse intrinsic pontine glioma (DIPG) are a phenotypically diverse collection of tumours which molecular profiling studies have refined into age- and location-based subgroups driven by unique genetic and epigenetic alterations. Despite this, the rarity of these tumours means individual studies are underpowered to investigate subgroup-specific changes and to identify low frequency events recurrently targeting specific pathways. We have collated genomic data from 143 unpublished cases up to 30 years of age along with those from 18 previously published datasets in an integrated analysis comprising 1067 unique patients across all anatomical compartments of the CNS, with histone mutation status available for 893 cases (n=67 H3.3G34R/V, n=316 H3.3K27M, n=68 H3.1/H3.2K27M). Distinct clinicopathological and molecular subgroups are defined by these histone mutations, with several novel co-segregating mutations identified through these aggregated analyses, including loss of FBXW7 in H3.3G34R/V, TOP3A rearrangements in H3.3K27M, and BCOR mutations in H3.1K27M. Histone wild-type subgroups can be further refined by the presence of key oncogenic events (PDGFRA or EGFR amplification/mutations) or methylation profiles which molecularly more closely resemble pleomorphic xanthoastrocytoma or other low grade gliomas. The burden of genomic aberrations decreases with age, highlighting the infant population as biologically and clinically distinct. Across the whole cohort, we identify novel previously unrecognised pathway dysregulation in small subsets of tumours (e.g. splicing, WNT, immune response) which make up the umbrella classification of paediatric diffusely infiltrating gliomas. The integrated dataset helps to further define the molecular diversity of the disease, opening up novel avenues for biological study and providing a basis for functionally defined future treatment stratification.

  Study EGAS00001002314

• Genomic profiling of paediatric high grade gliomas from the HERBY clinical trial

  The HERBY trial was a phase-II open-label, randomised, multicentre trial evaluating bevacizumab in patients with newly-diagnosed non-brainstem high grade glioma (HGG) between the ages of 3-18yrs. 121 patients were randomised with 1yr event-free survival as the primary end-point. Confirmation of HGG diagnosis by reference pathology was mandatory before randomisation, followed by review with five independent expert neuropathologists. We collected specimens from 89 patients consenting to translational research, and performed Sanger sequencing for H3F3A, Illumina 450K BeadArray methylation profiling, and whole exome sequencing and RNA sequencing where possible. 7/89 patients (8%) harboured H3F3A G34R/V mutations, whilst 24/89 (27%) harboured H3F3A K27M, the latter reflecting a high proportion of the novel entity recognised in the 2016 WHO classification of diffuse midline glioma with H3K27M mutation. Both histone mutations conferred a significantly shorter progression-free (p=0.0104) and overall survival (p=0.00159). 450K methylation subtyping additionally identified a number of subgroups in histone wild-type cases. These included infrequent (4/86, 4.6%) cases with IDH1 mutation in older children, further categorised into astrocytic ATRX-mutant (n=3) or 1p19q co-deleted with TERT promoter mutation (n=1). 9/74, (12%) of cases were classified as biologically resembling pleomorphic xanthoastrocytoma (PXA) by methylation profiling, with 5/9 harbouring BRAF V600E mutations and epithelioid histology, and 3/9 NF1 mutation and giant cell features. Three cases had methylation profiles more closely resembling low grade gliomas, though were morphologically high grade, and harboured MAPK dysregulation in two cases, with the third part of an additional infant (<3yrs) cohort. Finally, four cases had a mutational burden several orders of magnitude higher than the rest, with (2197-5332) somatic coding variants/sample reflecting a hypermutator phenotype. These data provide an important insight into the wide biological diversity of ‘HGG’ found in the paediatric age group and allow for a profound refinement of clinical trial interpretation.

  Study EGAS00001002328

• The evolutionary trajectory of a highly recurrent paediatric high grade neuroepithelial tumour with MN1:BEND2 fusion

  Astroblastomas are rare brain tumours which predominate in children and young adults, and have a controversial claim as a distinct entity, with no established WHO grade. Reports suggest a better outcome than high grade gliomas, though they frequently recur. Recently, they have been described to overlap with a newly-discovered group of tumours described as ’high grade neuroepithelial tumour with MN1 alteration’ (CNS HGNET-MN1), defined by global methylation patterns and strongly associated with gene fusions targeting MN1. We have studied a unique case of astroblastoma arising in a 6 year old girl, with multiple recurrences over a period of 10 years, with the pathognomonic MN1:BEND2 fusion. Exome sequencing allowed for a phylogenetic reconstruction of tumour evolution, which when integrated with clinical, pathological and radiological data provide for a thorough understanding of disease progression, with initial treatment driving tumour dissemination along four distinct trajectories. Infiltration of distant sites was associated with a later genome doubling, whilst there was evidence of convergent evolution of different lesions acquiring distinct alterations targeting NF-κB. These data represent an usual opportunity to understand the evolutionary history of a highly recurrent childhood brain tumour, and provide novel therapeutic targets for astroblastoma / CNS HGNET-MN1.

  Study EGAS00001002432

• 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

• Melanoma RNA Seq Data commitee

  Dac EGAC00001002928

• Genomic profiling of paediatric glioma cell lines

  Although paediatric gliomas resemble their adult counterparts in many ways, there appear to be distinct clinical and biological differences. One important factor hampering the development of new targeted therapies is the relative lack of cell lines derived from childhood glioma patients, as it is unclear whether the well-established adult lines commonly used are representative of the underlying molecular genetics of childhood tumours. We have carried out a detailed molecular profiling of paediatric glioma cell lines KNS42, SF188, UW479, RES259 and RES186.

  Study EGAS00001003006

• Targets of MEK inhibition in DIPG

  The survival of children with diffuse intrinsic pontine glioma (DIPG) remains dismal, with new treatments desperately needed. In the era of precision medicine, targeted therapies represent an exciting treatment opportunity, yet resistance can rapidly emerge, playing an important role in treatment failure. In a prospective biopsy-stratified clinical trial, we combined detailed molecular profiling (methylation BeadArray, exome, RNAseq, phospho-proteomics) linked to drug screening in newly-established patient-derived models of DIPG in vitro and in vivo. We identified a high degree of in vitro sensitivity to the MEK inhibitor trametinib in samples which harboured genetic alterations targeting the MAPK pathway, including the non-canonical BRAF_G469V mutation, and those affecting PIK3R1 and NF1. However, treatment of PDX models and of a patient with trametinib at relapse failed to elicit a significant response. We generated trametinib-resistant clones in the BRAF_G469V model through continuous drug exposure, and identified acquired mutations in MEK1/2 (MEK1_K57N, MEK1_I141S and MEK2_I115N) with sustained pathway up-regulation. These cells showed the hallmarks of mesenchymal transition, and expression signatures overlapping with inherently trametinib-insensitive primary patient-derived cells that predicted an observed sensitivity to dasatinib. Combinations of trametinib with dasatinib and the downstream ERK inhibitor ulixertinib showed highly synergistic effects in vitro. These data highlight the MAPK pathway as a therapeutic target in DIPG, and show the importance of parallel resistance modelling and rational combinatorial treatments likely to be required for meaningful clinical translation.

  Study EGAS00001004495

• 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

• Malignant pheochromocytomas/paragangliomas harbor mutations in transport and cell adhesion genes

  One out of ten patients with pheochromocytoma (PCC) and paraganglioma (PGL) develop malignant disease. Today there are no reliable pathological methods to predict malignancy at the time of diagnosis. Tumors harboring mutations in the succinate dehydrogenase subunit B (SDHB) gene often metastasize but the sequential genetic events resulting in malignant progression are not fully understood. The aim of this study was to identify somatic mutations that contribute to the malignant transformation of PCC/PGL. We performed pair-wise (tumor-normal) whole-exome sequencing to analyze the somatic mutational landscape in five malignant and four benign primary PCC/sympathetic PGL (sPGL), including two biological replicates from each specimen. In total, 225 unique somatic mutations were identified in 215 genes, with an average mutation rate of 0.54 mutations/megabase. Malignant tumors had a significantly higher number of mutations compared to benign tumors (p<0.001). Three novel genes were identified as recurrently mutated; MYCN, MYO5B and VCL, and mutations in these genes were exclusively found in malignant sPGL tumors. Mutations in the MYO5B gene could be verified in two publicly available data sets. A gene ontology analysis of mutated genes showed enrichment of cellular functions related to cytoskeletal protein binding, myosin complex, and motor activity, many of which had functions in Rab and Rac/Rho GTPase pathways. Conclusively, we have identified recurrent mutations in genes related to intracellular transport and cell adhesion, and we have confirmed recurrent mutations of MYO5B in PCC/PGL cases with malignant potential. Our study suggests that deregulated Rab and Rac/Rho pathways may be important in PCC/PGL tumorigenesis.

  Study EGAS00001001601

• 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

• Exome-sequencing of human B cell lymphoma cell lines

  26 cell lines derived from human Diffuse Large B Cell lymphomas (DLBCL) or Burkit Lymphomas (BL) were subjected to whole exome sequencing. Exome capture was carried out using the SeqCap EZ Exome Library 2.0 kit (Roche/Nimblegen) and 100 bp single-read sequencing was performed on a HiSeq2500 (Illumina). 82% of the coding region was covered at least 30x.

  Study EGAS00001001463

• 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

• Copy number profiling of primary samples and cell lines of retinoblastoma

  While RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes. We performed SCNA profiling of 45 primary retinoblastoma samples and 8 retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N=310). Comprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N=34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades. Since focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.

  Study EGAS00001001715

• 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

• Hodgkin's Lymphoma WGS Commitee

  Dac EGAC00001003015

• Transcriptional regulatory networks of tumor-associated macrophages that drive malignancy in mesenchymal glioblastoma

  Background: Glioblastoma (GBM) is a complex disease with extensive molecular and transcriptional heterogeneity. GBM can be subcategorized into four distinct subtypes; tumors that shift towards the mesenchymal phenotype upon recurrence are generally associated with treatment resistance, unfavorable prognosis, and the infiltration of pro-tumorigenic macrophages. Results: We explore the transcriptional regulatory networks of mesenchymal-associated tumor-associated macrophages (MA-TAMs), which drive the malignant phenotypic state of GBM, and identify macrophage receptor with collagenous structure (MARCO) as the most highly differentially expressed gene. MARCOhigh TAMs induce a phenotypic shift towards mesenchymal cellular state of glioma stem cells, promoting both invasive and proliferative activities, as well as therapeutic resistance to irradiation. MARCOhigh TAMs also significantly accelerate tumor engraftment and growth in vivo. Moreover, both MA-TAM master regulators and their target genes are significantly correlated with poor clinical outcomes, and are often associated with genomic aberrations in neurofibromin 1 (NF1) and phosphoinositide 3-kinases/mammalian target of rapamycin/Akt pathway (PI3K-mTOR-AKT)-related genes. We further demonstrate origination of MA-TAMs from peripheral blood, as well as their potential association with tumor-induced polarization states and immunosuppressive environments. Conclusions: Collectively, our study characterizes the global transcriptional profile of TAMs driving mesenchymal GBM pathogenesis, providing potential therapeutic targets for improving the effectiveness of GBM immunotherapy.

  Study EGAS00001002443

• 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

• Searching for variants associated with endometriosis

  We explored regulatory mechanism of a SNP on chromosome 9p21 associated with endometriosis by leveraging “allele-specific” functional genomic approaches. By re-sequencing 1.29 Mb of 9p21 region and scrutinizing DNase-seq data from the ENCODE project, we prioritized rs17761446 as a candidate functional variant that was in perfect linkage disequilibrium with the original GWAS SNP (rs10965235) and located on DNase I hypersensitive site. Chromosome conformation capture followed by high-throughput sequencing revealed that the protective G allele of rs17761446 exerted stronger chromatin interaction with ANRIL promoter. We demonstrated that the protective allele exhibited preferential binding affinities to TCF7L2 and EP300 by bioinformatics and chromatin immunoprecipitation (ChIP) analyses. ChIP assays for histone H3 lysine 27 acetylation and RNA polymerase II reinforced the enhancer activity of the SNP site. The allele specific expression analysis for eutopic endometrial tissues and endometrial carcinoma cell lines showed that rs17761446 was a cis-regulatory variant where G allele was associated with increased ANRIL expression.

  Study EGAS00001001741

• Elucidating the genomic architecture of Asian EGFR-mutant lung adenocarcinoma through multi-region exome sequencing

  EGFR-mutant lung adenocarcinomas (LUAD) display diverse clinical trajectories and are characterized by rapid but short-lived responses to EGFR tyrosine kinase inhibitors (TKIs). Through sequencing of 79 spatially distinct regions from 16 early stage tumors, we show that despite low mutation burdens, EGFR-mutant Asian LUADs unexpectedly exhibit a complex genomic landscape with frequent and early whole genome doubling, aneuploidy and high clonal diversity. Multiple truncal alterations, including TP53 mutations and loss of CDKN2A and RB1, converge on cell cycle dysregulation, with late sector-specific high amplitude amplifications and deletions that potentially beget drug resistant clones. We highlight the association between genomic architecture and clinical phenotypes, such as co-occurring truncal drivers and primary TKI resistance. Through comparative analysis with published smoking-related LUAD, we postulate that the high intra-tumor heterogeneity observed in Asian EGFR-mutant LUAD may be contributed by an early dominant driver, genomic instability and low background mutation rates.

  Study EGAS00001001736

• 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

• Analysis of tumor periphery and center-specific mutations in renal cell carcinoma

  Renal cell carcinoma (RCC) is a genomically heterogeneous tumor. In the present project, the question whether intratumoral heterogeneity follows a zonal pattern indicating spatial niches was addressed. Whole exome sequencing of 16 paired samples from tumor periphery and center revealed a number of region-specific functional SNVs and Indels. Therefore, RCCs are not composed of evenly admixed tumor cells but show topological differences in their clonal composition.

  Study EGAS00001001784

• Mutational landscape of renal cell carcinoma with venous tumor thrombus

  Renal cell carcinoma with venous tumor thrombus opens the unique opportunity to investigate the clonal evolution of a tumor in two fundamentally different compartments i.e., within the confinement of an organ and inside a large blood vessel. We performed multiregion whole exome sequencing of a total of 37 samples from five consecutive patients (normal tissue, n=5; primary tumors, n=16; tumor thrombi, n=16) to >35-fold target coverage. Four patients had a clear cell RCC, one patient had a poorly differentiated type II papillary RCC (RCC-VTT-04). The latter patient had a friable thrombus, the others were of solid consistency.

  Study EGAS00001001950

• Exome sequencing data from two myelosarcomas

  Isolated myeloid sarcoma (MS, i.e. acute myeloid leukemia without detectable bone marrow infiltration) is rare and most patients subsequently develop overt AML. We prospectively performed whole exome sequencing of uterine MS, concomitant morphologically normal bone marrow (BM) and subsequently developing AML from two women. Both patients harbored preleukemic clones in the morphologically normal BM, including DNMT3A mutations. In addition, recurrent mutations were identified in the NFE2 gene, a gene not previously associated with AML or isolated MS. Retrospective studies of further MS samples, twelve with and four without preceding or concomitant bone marrow AML, revealed NFE2 mutations in two additional cases, both with isolated MS and prior myeloid malignancies. We conclude that DNMT3A mutated preleukemia may predispose for the frequent bone marrow relapses following isolated MS and that recurrent NFE2 mutations may have a role in the development of isolated MS.

  Study EGAS00001002562

• Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study

  This data comes from a genome-wide association study in 1090 cases of scarring trachoma and 1531 controls from a non-selected multi-ethnic sampling of The Gambia. The data constitutes hard base calls at 1457295 features from the HumanOmni2.5-8 (Omni2.5) BeadChip. Three genotype-calling algorithms were used : Illuminus (Bioinformatics 2007, 23:2741–6), GenCall (Patent US7035740, 2006) and GenoSNP (Bioinformatics 2008, 24:2209–14). Data from each calling algorithm was filtered to retain only SNPs with a call rate ≥ 0.98. To obtain a merged set of SNPs, all genotypes that matched across call-sets were retained, whilst those that mismatched between call-sets were set to missing. Genotypes that were present in one call set and missing in others were also retained. The merged SNP set contained 1467876 SNPs, which were filtered for call rate ≥ 0.99 and Hardy-Weinberg equilibrium P value < 5 x 10-8. The case and control groups were approximately equivalent with respect to gender and ethnicity. 70% of cases and 63% of controls were female. Self-described ethnic composition of the cases was 29% Jola, 27% Mandinka, 21% Wolof, 10% Fula and 13% other/no data. In controls the composition was 25% Jola, 24% Mandinka, 21% Wolof, 8% Fula and 22% other/no data. Median age was 49 (range 32-60) in cases and 37 (range 12-52) in controls (t = -3.6019, P = 0.0003). This data set has substantial levels of population structure, which reflects extensive networks of pairwise kinship between participants. First pass tests of association in EMMAX (Nat Genet 2010, 42:348–54) led to genome-wide deflation of the test statistics (λ = 0.982 SE = 1.74e-05) that could be directly attributed to the differing age and gender distributions between cases and controls. Modelling the phenotype to adjust for age and gender successfully controlled for this and no genome-wide deviation from the null was subsequently observed (λ = 1.001, SE = 8.7 x 10-7). Neither principal components analysis (PCA) nor tests of proportional identity by state (IBS) variance in PLINK identified significant levels of within or between group genetic variance. The phenotypes submitted in the phenotypes file include both uncorrected (case/control status) and corrected (age and gender corrected case/control status) phenotypes, as well as age and self-described ethnicity.

  Study EGAS00001001516

• 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

• Clonal selection and double hit events involving tumor suppressor genes underlie relapse from total therapy

  In order to gain insights into the mechanisms of disease progression in multiple myeloma (MM) patients treated with multi-agent and high-dose chemotherapy we have performed a longitudinal study of 33 patients entered into Total Therapy protocols using gene expression profiling, high resolution copy number analyses and whole exome sequencing. The results of this study show the importance of acquired mutations in MM driver genes, bi-allelic events affecting tumor suppressor genes and increased proliferation rates in driving relapse, which is enhanced by intra-clonal heterogeneity and Darwinian type clonal evolution. Branching was the most frequent evolution pattern, characterized by clonal loss and appearance of new clones. We found a higher number of lost and acquired copy number aberrations and an excess of bi-allelic inactivation of tumor suppressor genes in GEP70 high risk (HiR) cases, consistent with genomic instability being a key feature of HiR MM. In conclusion, our study further stresses the impact of known prognostic markers and driver genes on the fitness level of MM cells to survive multi-agent chemotherapy. The frequent bi-allelic loss of tumor suppressor genes in HiR cases and the potentially resulting failure of treatment with DNA damaging agents highlight the need for new therapies for HiR MM.

  Study EGAS00001001810

• 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

• EATL-II STUDY

  Enteropathy-associated T-cell lymphoma (EATL), a rare and aggressive intestinal malignancy of intraepithelial T lymphocytes, comprises two disease variants (EATL-I and EATL-II) differing in clinical characteristics and pathological features. Here we report findings derived from whole exome sequencing of 15 EATL-II tumor-normal tissue pairs. The tumor suppressor gene SETD2 encoding a non-redundant H3K36-specific trimethyltransferase is altered in 14/15 cases (93%), mainly by loss-of-function mutations and/or loss of the corresponding locus (3p21.31). These alterations consistently correlate with defective H3K36 trimethylation. The JAK/STAT pathway comprises recurrent STAT5B (60%), JAK3 (46%) and SH2B3 (20%) mutations, including a STAT5B V712E activating variant. Additionally, frequent mutations in TP53, BRAF and KRAS are observed. Conversely, in EATL-I, no SETD2, STAT5B or JAK3 mutations are found, and H3K36 trimethylation is preserved. This study describes SETD2 inactivation as EATL-II molecular hallmark, supports EATL-I and -II being two distinct entities, and defines potential new targets for therapeutic intervention.

  Study EGAS00001001879

• 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

• Genomic profiling of esthesioneuroblastoma

  We performed whole exome and whole genome sequencing on a cohort of esthesioneurblastoma to understand the genetic underpinnings of this neoplasm.

  Study EGAS00001003225

• 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

• 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

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

  Study EGAS00001003798

• Genetic Basis of Hepatosplenic T Cell Lymphoma (HSTL)

  Hepatosplenic T cell lymphoma (HSTL) is a rare and lethal lymphoma; the genetic drivers of this disease are unknown.  Through whole exome sequencing of 68 HSTLs, we define recurrently mutated driver genes and copy number alterations in the disease. Chromatin modifying genes including SETD2, INO80 and ARID1B were commonly mutated in HSTL, affecting 62% of cases. HSTLs manifest frequent mutations in STAT5B (31%), STAT3 (9%), and PIK3CD (9%) for which there currently exist potential targeted therapies. In addition, we noted less frequent events in EZH2, KRAS and TP53. SETD2 was the most frequently silenced gene in HSTL. We experimentally demonstrated that SETD2 acts as a tumor suppressor gene. In addition, we found that mutations in STAT5B and PIK3CD activate critical signaling pathways important to cell survival in HSTL. Our work thus defines the genetic landscape of HSTL and implicates novel gene mutations linked to HSTL pathogenesis and potential treatment targets.

  Study EGAS00001002182

• Changes in alternative splicing and associated neo-antigens due to therapy

  Recent studies have identified alternative splicing (AS) as a novel source of neoantigens for immunotherapy. Surprisingly little is known about the AS milieu in recurrent glioblastoma (GBM), despite this being the venue for most clinical trials. We profiled 29 primary-recurrent paired human GBM specimens via RNA sequencing and re-analyzed RNA-sequencing data from non-malignant human brain tissues. From these data, we reconstructed the landscape of AS in GBM through recurrence and contrasted that to isoforms observed in non-malignant brain. The AS events we identified were cross-referenced with single-cell GBM atlases to determine cell-type specific splicing patterns. From this we identified novel splicing events in cell-surface proteins that are suitable targets for engineered T-cell therapies. We found recurrent-specific isoforms of mitogen-activated kinase pathway genes which are expressed exclusively by GBM stem-like cells that enhance invasiveness. These studies shed light on the effect of therapy on AS and identify novel targets for emerging immunotherapies.

  Study EGAS00001004524

• Inactivation of TGFβ receptors in stem cells drives cutaneous squamous cell carcinoma - 30 whole exomes

  Melanoma patients treated with oncogenic BRAF inhibitors can develop cutaneous squamous cell carcinoma (cSCC) within weeks of treatment, driven by paradoxical RAS/RAF/MAPK pathway activation. Here, we identify frequent TGFBR1 and TGFBR2 mutations in human vemurafenib-induced skin lesions and in sporadic cSCC. Functional analysis reveals these mutations ablate canonical TGFβ Smad signaling which is localised to bulge stem cells in both normal human and murine skin. MAPK pathway hyperactivation (through BrafV600E or KrasG12D knockin) and TGFβ signaling ablation (through Tgfbr1 deletion) in LGR5+ve stem cells enables rapid cSCC development in the mouse. Mutation of Tp53 (which is commonly mutated in sporadic cSCC) coupled with Tgfbr1 deletion in LGR5+ve cells also results in cSCC development. These findings indicate that LGR5+ve stem cells can act as cells of origin for cSCC and that RAS/RAF/MAPK pathway hyperactivation or Tp53 mutation, coupled with loss of TGFβ signaling, are driving events of skin tumorigenesis.

  Study EGAS00001001892

• 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

• Exome sequence data for germline, primary tumor, and relapse tumor of a transformed non-Hodgkins lymphoma patient with unexpected long-time remission

  The patient was diagnosed with a transformed high-grade non-Hodgkin lymphoma stage III, and treated with conventional R-CHOP (Rituximab (R), Cyclophosphamide (C), Doxorubicin (H), Vincristine (O) and Prednisone (P)). Unfortunately, she suffered from severe toxicity but recovered during the following six months remission until a biopsy verified relapse. The patient refused second line combination chemotherapy, but accepted a 3 months palliation with R and Chlorambucil. Unexpectedly, she obtained continuous complete remission and is at present >9 years after primary diagnosis.Molecular studies and data evaluation by principal component analysis, mutation screening, and copy number variations of the primary and relapsed tumor, identified patterns of branched lymphoma evolution, most likely diverging from an unknown in situ follicular lymphoma.

  Study EGAS00001001973

• Genomics of enteropathy associated T cell lymphoma (EATL)

  Enteropathy-associated T cell lymphoma (EATL) is a lethal and the most common neoplastic complication of celiac disease. In this study, we applied whole exome sequencing in 69 EATL tumors to define the genetic landscape of EATLs. The most frequently mutated gene in EATL was SETD2 (32% of cases), with loss of function mutations comprising a large proportion (41%) of those mutations. The JAK-STAT pathway was the most frequently mutated signaling pathway, with mutations in STAT5B being the most common, followed by JAK1, JAK3, STAT3, and the negative regulator of the pathway, SOCS1. We also identified mutations in NRAS and KRAS, as well as TP53 and TERT. Interestingly, the two types of EATL had highly overlapping genetic alterations indicating shared mechanisms underlying their pathogenesis. Gene expression analysis of the two types of EATL identify robust differences in transcriptional profiles. We further modeled the effects of SETD2 loss in vivo by developing a T cell specific knockout mouse model. These mice manifested an expansion of γδ T cells that identify novel roles for SETD2 in T cell development and lymphomagenesis. Together, our data render the most comprehensive genetic portrait yet of this uncommon but lethal disease, and may inform future classification schemes.

  Study EGAS00001001954

• Single-cell RNA sequencing on 5063 single T cells isolated from peripheral blood, tumour and adjacent normal tissues from six hepatocellular carcinoma patients.

  Single-cell RNA sequencing on 5063 single T cells isolated from peripheral blood, tumour and adjacent normal tissues from six hepatocellular carcinoma patients. The transcriptional properties of these individual cells, coupled with the assembled T cell receptor (TCR) sequences, were used to delineate the molecular subtypes, lineage, and functional properties of each T cell group. For five HCC patient, the exome-seq data of matched normal/tumor and bulk RNA-seq data of tumor are also available.

  Study EGAS00001002072

• 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

• Melanoma Institute Australia Data Access Committee

  Dac EGAC00001003075

• 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

• RNA-sequencing of six Pilocytic astrocytoma tumors

  Pilocytic astrocytoma (PA) is the most common pediatric brain tumor. A recurrent feature of PA is deregulation of the mitogen activated protein kinase (MAPK) pathway most often through KIAA1549-BRAF fusion, but also by other BRAF- or RAF1-gene fusions and point mutations (e.g. BRAFV600E). These features may serve as diagnostic and prognostic markers, and also facilitate development of targeted therapy. The aim of this study was to characterize the genetic alterations underlying the development of PA tumor in six cases, and evaluate methods for fusion oncogene detection. Using a combined analysis of RNA sequencing and copy number variation data we identified a new BRAF fusion involving the 5’ gene fusion partner GTF2I (7q11.23), not previously described in PA. The novel GTF2I-BRAF 19-10 fusion was found in one case, while the other five cases harbored the frequent KIAA1549-BRAF 16-9 fusion gene. Comparing fusion detection methods, Fluorescence in situ hybridization with BRAF break apart probe was the most sensitive method for detection of the two different BRAF rearrangements (GTF2I-BRAF and KIAA1549-BRAF). Our finding of a novel BRAF fusion in PA further emphasis the important role of B-Raf in tumorigenesis of these tumor types. Also, the growing list of BRAF/RAF gene fusions suggests these to be informative tumor markers in molecular diagnostics, which could guide future treatment strategies.

  Study EGAS00001002149

• Genomic and transcriptomic profiling of signalling networks in follicular lymphoma

  This dataset is composed of sequencing data from 16 patients diagnosed with follicular lymphoma, who developed transformed disease and for whom representative frozen tumor biopsy samples were available. For 14 patients there is paired material from both disease stages. DNA and RNA libraries were constructed and captured with a kinome bait library (Agilent Technologies), before sequencing on an Illumina platform.

  Study EGAS00001002175

• 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

• Phylogenetic evolution of metastatic melanoma.

  To explore molecular factors associated with evolution and resistance in melanoma evolution, we multi-lesion sequenced primary and metastatic sites collected from autopsy in seven patients.

  Study EGAS00001003582

• Gene fusion and transcriptomic landscapes of sarcomas

  Sarcoma represents a highly heterogeneous group of tumors. We report here the first unbiased and systematic search for gene fusions and analysis of transcriptomic profiles in 100 monomorphic sarcomas cases using RNA-sequencing profiling. Fusion genes were detected in two thirds of samples. Similarly to fusion genes such as PAX3/7-FOXO1, EWSR1/FUS-ETS, SS18-SSX, BRD3/4-NUTM1 or EWSR1/FUS/TAF15-NR4A3 characterizing well-defined entities, we show that all sarcomas displaying CIC fusions- whatever their fusion partner (being DUX4, NUTM1 or FOXO4) - form a transcriptionally homogeneous group of tumors. Likewise, tumors with either a BCOR fusion (with CCNB3, MAML3 or ZC3H7B partners) or with a BCOR internal duplication form a single biological entity. We also found that EML4-ALK fusions also characterize a homogeneous infantile fibrosarcoma subgroup. In contrast, fusions like VGLL2-NCOA2/CITED2, TMP3/TRP-NTRK1 or ETV6-NTRK3 are observed in more heterogeneous tumors. Finally, we also describe a new group of bone sarcomas characterized by EWSR1- or FUS-TFCP2 fusions.

  Study EGAS00001002189

• Exome sequencing of pseudomyxoma peritonei

  Pseudomyxoma peritonei (PMP) is a subtype of mucinous adenocarcinoma mainly restricted to the peritoneal cavity and most commonly originating from the appendix. In this study, we sequenced whole coding genome of nine PMP tumors and paired normal tissues in order to identify commonly mutated genes and signaling pathways affected in PMP.

  Study EGAS00001002418

• 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

• Genomic characterization of NUT midline carcinoma

  In this study, we sequenced two NUT midline carcinoma genomes. Analysis of rearrangement breakpoint provided information on the generating event of BRD-NUT rearrangements. Mutational signature analysis also revealed biological process during tumorigenesis of NUT midline carcinomas.

  Study EGAS00001001934

• Evolutionary Genome Analysis of Transformation into Small Cell Carcinomas from Lung Adenocarcinomas

  Anecdotal case have been reported lung adenocarcinoma could occasionally transform into small cell carcinoma. Most of the reported cases were EGFR-mutant lung adenocarcinomas that transform into small cell carcinoma during EGFR tyrosine kinase inhibitor treatment. However, some cases also indicate lung adenocarcinomas without EGFR mutation could also transform into small cell carcinomas. Here we investigate this transformation process by whole genome sequencing of transformed small cell carcinomas as well as serially acquired samples of lung adenocarcinoma into small cell carcinomas.

  Study EGAS00001001757

• Transcriptional mechanisms of resistance to anti-PD-1 therapy

  Purpose: To explore factors associated with response and resistance to anti-PD-1 therapy, we analyzed multiple disease sites at autopsy in a patient with widely metastatic melanoma who had a heterogeneous response.Materials and Methods: Twenty-six melanoma specimens (four pre-mortem, 22 post-mortem) were subjected to whole-exome sequencing. Candidate immunologic markers and gene expression were assessed in ten cutaneous metastases showing response or progression during therapy.Results: The melanoma was driven by biallelic inactivation of NF1. All lesions had highly concordant mutational profiles and copy number alterations, indicating linear clonal evolution. Expression of candidate immunologic markers was similar in responding and progressing lesions. However, progressing cutaneous metastases were associated with over-expression of genes associated with extracellular matrix and neutrophil function.Conclusions: Although mutational and immunologic differences have been proposed as the primary determinants of heterogeneous response/resistance to targeted therapies and immunotherapies, respectively, differential lesional gene expression profiles may also dictate anti-PD-1 outcomes.

  Study EGAS00001002195

• An Unusual Genomic Variant of Pancreatic Ductal Adenocarcinoma with an Indolent Clinical Course

  We describe a patient with biopsy-proven locally advanced pancreatic ductal adenocarcinoma (PDA) who underwent a modified course of gemcitabine and stereotactic body radiation therapy. Whole exome sequencing of this patient's tumor revealed a simple genome landscape with no evidence of mutations, copy number changes, or structural alterations in genes most commonly associated with PDA, i.e. KRAS, CDKN2A, TP53 or SMAD4. An analysis of germline DNA revealed no pathogenic variants of significance. Whole exome and whole genome sequencing identified a somatic mutation of RNF213 and an inversion/deletion of CTNNA2 as the genetic basis of this PDA. Although PDA is classically characterized by a predictable set of mutations, these data suggest that alternate genetic paths to PDA may exist, which can be associated with a more indolent clinical course.

  Study EGAS00001002192