Search Results - EGA European Genome-Phenome Archive

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

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

Study EGAS00001001940

Inferring expressed genes by whole-genome sequencing of plasma DNA

The analysis of cell-free DNA (cfDNA) in plasma represents a rapidly advancing field in medicine. cfDNA consists predominantly of nucleosome-protected DNA shed into the bloodstream by cells undergoing apoptosis. We performed whole-genome sequencing (WGS) of plasma DNA and identified two discrete regions at transcription start sites (TSS) where the nucleosome occupancy results in different read-depth coverage patterns in expressed and silent genes. By employing machine learning for gene classification we found that the plasma DNA read depth patterns from healthy donors reflected the expression signature of hematopoietic cells. In cancer patients with metastatic disease, we were able to classify expressed cancer driver genes in regions with somatic copy number gains with high accuracy. We could even determine the expressed isoform of genes with several TSSs as confirmed by RNA-Seq of the matching primary tumor. Our analyses provide functional information about the cells releasing their DNA into the circulation.

Study EGAS00001001754

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

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

Chordoma_Sequencing_Project_Whole_Genome

This is a continuation of the Chordoma Sequencing Project, prelim #G0175-122-CF. All cancers arise due to somatically acquired abnormalities in DNA sequence. Systematic sequencing of cancer genomes allows acquisition of complete catalogues of all classes of somatic mutation present in cancer. These mutation catalogues will allow identification of the somatically mutated cancer genes that are operative and characterise patterns of somatic mutation that may reflect previous exogenous and endogenous mutagenic exposures. In this application, we aim to perform whole genome sequencing on 10 chordoma matched genome pairs. The Chordoma Foundation will fund the RNA Sequencing/Methylation and SNP6 costs for this project while the whole genome and library sequencing will be covered under the CGP, 2012-2013 core WTSI internal funding. An additional sequencing of three cancer cell lines will be added to this work, financed by the Chordoma Foundation. In this project, the Chordoma foundation will pay for the additional costs of methylation, RNA Seq,

Study EGAS00001000409

APCDR Uganda GWAS: Genome-wide sequence variation and susceptibility loci for cardiometabolic traits in a sub-Saharan African population (UG2G component)

Genomic studies in African populations provide unique opportunities to understand disease aetiology, human genetic diversity and population history in a regional and a global context. To leverage the relative benefits of different strategies, we undertook a combined approach of genotyping and whole-genome sequencing (WGS) in a population-based study of 6,400 individuals from a geographically defined rural community in South-West Uganda. We present data from 4,778 individuals with genotypes for ~2.2 million SNPs from the Uganda GWAS resource (UGWAS), and sequence data on up to 1,978 individuals spanning 41.5M SNPs and 4.5M indels (UG2G); 343 individuals overlap between the two datasets. We highlight the value of the largest sequence panel from Africa to date as a global resource for variant discovery, imputation and understanding the mutational spectrum and its clinical relevance in African populations. Alongside phenotype data, we provide a rich new genomic resource for researchers in Africa and globally.

Study EGAS00001000545

CD74-NRG1 fusions in lung adenocarcinoma

We discovered a novel somatic gene fusion, CD74-NRG1, by transcriptome sequencing of 25 lung adenocarcinomas of never smokers. By screening 102 lung adenocarcinomas negative for known oncogenic alterations we found four additional fusion-positive tumors, all of which were of the invasive mucinous subtype. Mechanistically, CD74-NRG1 leads to extracellular expression of the EGF-like domain of NRG1 III-beta3, thereby providing the ligand for ERBB2-ERBB3 receptor complexes. Accordingly, ERBB2 and ERBB3 expression was high in the index case and expression of phospho-ERBB3 was specifically found in tumors bearing the fusion (p<0.0001). Ectopic expression of CD74-NRG1 in lung cancer cell lines expressing ERBB2 and ERBB3 activated ERBB3 and the PI3K-AKT pathway, and led to increased colony formation in soft agar. Thus, CD74-NRG1 gene fusions are activating genomic alterations in invasive mucinous adenocarcinomas and may offer a therapeutic opportunity for a lung tumor subtype with, so far, no effective treatment.

Study EGAS00001000653

Biallelic DICER1 mutations in sporadic pleuropulmonary blastoma

Pleuropulmonary blastoma (PPB) is an extremely rare pediatric malignancy in the lung, whose pathogenesis is poorly understood, except for recent reports of frequent germline heterozygous DICER1 mutations. To investigate the genetic basis of PPB, we performed whole-exome sequencing in 7 representative PPB cases, followed by targeted deep sequencing in 12 cases with PPB. DICER1 mutations were found in 11/12 cases. Biallelic DICER1 mutations were common in PPB, in which RNase IIIb domain mutations were found in all cases with or without nonsense/frameshift mutations and were somatic in all evaluable cases. A majority of cases had mutated DICER1 alleles in germline with or without an additional somatic mutation in the remaining allele, while other cases had exclusively somatic mutations involving the RNase IIIb domain. TP53 deletions/mutations were detected in 8/12 cases. Our results highlight a unique role of the RNase IIIb domain mutations and TP53 inactivation in PPB pathogenesis.

Study EGAS00001000662

Tracking the origins and drivers of subclonal metastatic expansion in prostate cancer

Tumour heterogeneity in primary prostate cancer is a well-established phenomenon. However, how the subclonal diversity of tumours changes during metastasis and progression to lethality is poorly understood. Here we reveal the precise direction of metastatic spread across four lethal prostate cancer patients using whole-genome and ultra-deep targeted sequencing of longitudinally collected primary and metastatic tumours. We find one case of metastatic spread to the surgical bed causing local recurrence, and another case of cross-metastatic site seeding combining with dynamic remoulding of subclonal mixtures in response to therapy. By ultra-deep sequencing end-stage blood, we detect both metastatic and primary tumour clones, even years after removal of the prostate. Analysis of mutations associated with metastasis reveals an enrichment of TP53 mutations, and additional sequencing of metastases from 19 patients demonstrates that acquisition of TP53 mutations is linked with the expansion of subclones with metastatic potential which we can detect in the blood.

Study EGAS00001000942

Origins and functional consequence of somatic mitochondrial DNA mutations

Recent sequencing studies have extensively explored the somatic alterations present in the nuclear genomes of cancers. Although mitochondria control energy metabolism and apoptosis, the origins and impact of cancer-associated mutations in mitochondrial DNA (mtDNA) are unclear. Here, we analysed somatic alterations in mtDNA from 1,675 tumors across 31 histologies. We identified 1,907 somatic substitutions, which exhibited dramatic replicative strand bias, predominantly C>T and A>G on the mitochondrial heavy strand. This strand-asymmetric signature differs from those found in nuclear cancer genomes but matches the inferred germline process shaping primate mtDNA sequence content. Numbers of mtDNA mutations showed considerable heterogeneity across tumor types. Missense mutations were selectively neutral and often gradually drifted towards homoplasmy over time. In contrast, mutations resulting in protein truncation undergo negative selection and were almost exclusively heteroplasmic. Our findings indicate that the endogenous mutational mechanism has far greater impact than any other external mutagens in mitochondria, and is fundamentally linked to mtDNA replication.

Study EGAS00001000968

Frequent alterations in cytoskeleton remodeling genes in primary and metastatic Chinese lung adenocarcinomas

The landscape of genetic alterations in lung adenocarcinoma derived from Asian patients is largely uncharacterized. Here we present an integrated genomic and transcriptomic analysis of 335 primary lung adenocarcinomas and 35 corresponding lymph node metastases from Chinese patients. Altogether 13 significantly mutated genes are identified, including themost commonly mutated gene TP53 and novel mutation targets such as RHPN2, GLI3 and MRC2. TP53 mutations are furthermore significantly enriched in tumours from patients harbouring metastases. Genes regulating cytoskeleton remodelling processes are also frequently altered, especially in metastatic samples, of which the high expression level of IQGAP3 is identified as a marker for poor prognosis. Our study represents the first large-scalesequencing effort on lung adenocarcinoma in Asian patients and provides a comprehensive mutational landscape for both primary and metastatic tumours. This may thus form a basis for personalized medical care and shed light on the molecular pathogenesis of metastatic lung adenocarcinoma.

Study EGAS00001000982

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

Trans-ethnic genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation

We carried out a trans-ethnic genome-wide association and replication study of blood pressure phenotypes amongst up to 320,251 individuals of East Asian, European and South Asian ancestry. We find genetic variants at 12 novel loci to be associated with blood pressure (P=3.9x10-11 to P=5.0x10-21). The sentinel blood pressure SNPs are enriched for association with DNA methylation at multiple nearby CpG sites, suggesting that at some of the loci identified DNA methylation may lie on the regulatory pathway linking sequence variation to blood pressure. The sentinel SNPs at the 12 novel loci point to genes involved in vascular smooth muscle (IGFBP3, KCNK3, PDE3A and PRDM6) and renal (ARHGAP24, OSR1, SLC22A7, TBX2) function. The novel and known genetic variants predict increased left ventricular mass, circulating levels of NT-proBNP, cardiovascular and all-cause mortality (P=0.04 to 8.6x10-6). Our results provide new evidence for the role of DNA methylation in blood pressure regulation.

Study EGAS00001001427

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

A Genomic History of Aboriginal Australia

The population history of Aboriginal Australians remains largely uncharacterised. We generated high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting early population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene (past c.10 kya) associated with limited gene flow from this region to the rest of Australia consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out of Africa dispersal, and subsequently admixed with different archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Study EGAS00001001766

Rare disruptive mutations in ciliary function genes contribute to testicular cancer susceptibility

Testicular germ cell tumor (TGCT) is the most common cancer in young men1,2. Here we aimed to identify novel risk factors for TGCT using whole-exome sequencing, which was performed on 328 affected individuals from 153 families, 634 sporadic cases and 1,644 controls. We searched for genes that were recurrently affected by rare variants (minor allele frequency <0.01) with potentially damaging effects and evidence of segregation in families. 8.7% of families carried rare disruptive mutations in the cilia-microtubule genes (CMG) as compared to 0.5% of controls3 (P=2.1x10-8). The most significantly mutated CMG was DNAAF1 with biallelic inactivation and loss of DNAAF1 expression shown in tumors from carriers. DNAAF1 as a cause of TGCT was supported by a DNAAF1Hu255h(+/-) zebrafish model with 94% penetrance for TGCT compared to 14% in wildtype fish. These data implicate cilia-microtubule inactivation as a cause of TGCT development and are the first evidence for CMGs as cancer susceptibility genes.

Study EGAS00001001789

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

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

Divergence between high metastatic tumor burden and low circulating tumor DNA concentration in metastasized breast cancer

Circulating tumor DNA (ctDNA) was reported to represent a highly sensitive biomarker of metastatic cancer disease directly reflecting tumor burden and dynamics. Here we investigated the role of ctDNA in patients with metastatic breast cancer. In an index patient with more than 100,000 circulating tumor cells (CTCs) in serial blood analyses, whole genome, exome, or targeted deep sequencing of the primary tumor, metastases, and 551 CTCs were consistent with a genetically homogeneous cancer. However, the allele fractions (AFs) of ctDNA were only 2-3% in each analysis, which did neither reflect the tumor burden nor the dynamics of this progressive disease by far. Indeed, plasma analyses of 71 further patients demonstrated highly variable AFs of mutant fragments, which frequently did not correspond to the tumor burden. These results provide insights into mechanisms involved in CTC and ctDNA release into the circulation and have important implications for diagnostic tests based on liquid biopsies.

Study EGAS00001000625

Transcriptional and functional profiling defines human small intestinal macrophage subsets

Macrophages (Mf) are instrumental in maintaining immune homeostasis in the intestine, yet studies on origin and heterogeneity of human intestinal Mf are scarce. Here, we identified four distinct Mf subpopulations in human small intestine (SI). Assessment of their turnover in duodenal transplants revealed that all Mf subsets were completely replaced over time; Mf1 and Mf2, phenotypically similar to peripheral blood monocytes (PBMo), were largely replaced within 3 weeks, whereas two subsets with features of mature Mf, Mf3 and Mf4, exhibited significantly slower replacement. Mf3 and Mf4 localized differently in SI; Mf3 formed a dense network in mucosal lamina propria, whereas Mf4 was enriched in submucosa. Transcriptomic analysis showed that all Mf subpopulations subsets were markedly distinct from PBMo and dendritic cells (DC). Compared to PBMo, Mf subpopulations showed reduced responsiveness to proinflammatory stimuli, but were proficient at endocytosis of particulate and soluble material. These data provide a comprehensive analysis of human SI Mf population and suggest a precursor-progeny relationship with PBMo.

Study EGAS00001002093

Myelodysplastic cells in patients re-program mesenchymal stromal cells to establish a transplantable stem cell-niche disease unit.

Myelodysplastic syndromes (MDS) are a heterogeneous group of myeloid neoplasms with defects in hematopoietic stem/progenitor cells (HSPCs) and possibly the HSPC niche. Here we show that patient-derived mesenchymal stromal cells (MDS MSCs) display a disturbed differentiation program and are essential for the propagation of MDS-initiating lin-CD34+CD38- stem cells in orthotopic xenografts. Overproduction of niche factors such as N-Cadherin, IGFBP2, VEGFA and LIF is associated with the ability of MDS MSCs to enhance MDS expansion. These factors represent putative therapeutic targets to disrupt critical hematopoietic-stromal interactions in MDS. Finally, healthy MSCs adopt "MDS-MSC like" molecular features when exposed to hematopoietic MDS cells, indicative of an instructive remodeling of their microenvironment. This patient-derived xenograft model therefore provides functional and molecular evidence that MDS is a complex disease involving both the hematopoietic and stromal compartments. The resulting deregulated expression of niche factors may well also be a feature of other hematopoietic malignancies.

Study EGAS00001000716

Resequencing_candidate_genes_for_male_spermatogenic_impairment

This study aims to identify novel candidate variants from human Y-chromosomal genes DAZ, BPY2 and CDY1 by resequencing the coding regions of these genes from male patients with spermatogenic impairment. The coding regions of the genes have been amplified by standard PCR, amplicon lengths range from 244 to 486 bp. A total of 61 amplicons were amplified for each of the 96 patients, totalling to approx. 25 kb per sample. Amplicons were quantified by gel electrophoresis and pooled in approx. equimolar concentrations per patient. For each of the 96 submitted samples, approx. 1 microgram of amplified DNA pool is provided in a total volume of 120 microlitres. The samples need to be indexed and libraries prepared for a PE250bp Illumina MiSeq run. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001002157

Whole exome sequencing for gallbladder cancer in Xinhua Hospital Affiliated to Shanghai Jiao Tong University, School of Medicine

Patients with gallbladder carcinoma (GBC), the most aggressive malignancy of the biliary tract, have a poor prognosis. Here, we report our identification of somatic mutations of GBCs in 57 tumor-normal pairs by use of a combination of exome sequencing and ultra-deep sequencing of cancer-related genes. The mutation pattern is defined by a dominative prevalence of C>T mutations at TCN sites. Genes with a significant frequency of non-silent mutations include TP53 (47.1%), KRAS (7.8%), and ERBB3 (11.8%). Moreover, ErbB signaling (including EGFR, ERBB2, ERBB3, ERBB4 and their downstream genes) is the most extensively mutated pathway, affecting 36.8% (21 of 57) of the GBC samples. Multivariate analyses further reveal that patients with ErbB pathway mutations have a worse outcome (P = 0.001). These findings provide insight into the somatic mutational landscape in GBC and highlight the key role of the ErbB signaling pathway in the pathogenesis of GBCs.

Study EGAS00001000853

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

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

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

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

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

Study EGAS00001001570

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

Comparing nodal versus bony metastatic spread using tumour phylogenies

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

Study EGAS00001001801

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

Genetic landscape of Early T-cell precursor acute lymphoblastic leukaemia

Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL.

Study EGAS00001000348

Epigenomic alterations define lethal CIMP-positive ependymomas of infancy

Ependymomas are common childhood brain tumours that occur throughout the nervous system, but are most common in the paediatric hindbrain. Current standard therapy comprises surgery and radiation, but not cytotoxic chemotherapy as it does not further increase survival. Whole-genome and whole-exome sequencing of 47 hindbrain ependymomas reveals an extremely low mutation rate, and zero significant recurrent somatic single nucleotide variants. Although devoid of recurrent single nucleotide variants and focal copy number aberrations, poor-prognosis hindbrain ependymomas exhibit a CpG island methylator phenotype. Transcriptional silencing driven by CpG methylation converges exclusively on targets of the Polycomb repressive complex 2 which represses expression of differentiation genes through trimethylation of H3K27. CpG island methylator phenotype-positive hindbrain ependymomas are responsive to clinical drugs that target either DNA or H3K27 methylation both in vitro and in vivo. We conclude that epigenetic modifiers are the first rational therapeutic candidates for this deadly malignancy, which is epigenetically deregulated but genetically bland.

Study EGAS00001000443

Reference epigenomes generated as part of the International Human Epigenomics Consortium (IHEC)

The Centre for Epigenome Mapping Technologies (CEMT) is an epigenome sequencing platform funded by the Canadian Institutes of Health Research (CIHR) and Genome BC as part of the Canadian Epigenetics, Environment and Health Research Consortium (CEEHRC). CEMT aims to produce 100 human reference epigenomes, each comprising the following datasets: whole genome; whole genome bisulfite and oxidative bisulfite; transcriptome; microRNA; and chromatin immunoprecipitation (ChIP)-seq for the core chromatin marks (H3K4me1, H3K4me3, H3K9me3, H3K27me3, H3K36me3, and H3K27ac).CEEHRC is a full member of the International Human Epigenomics Consortium (IHEC), which aims to coordinate the production of 1,000 reference maps of human epigenomes for key cellular states relevant to health and diseases. As such, all datasets submitted by the CEMT platform conform to IHEC assay standards and metadata guidelines.

Study EGAS00001000552

Genomic catastrophes frequently arise in esophageal adenocarcinoma and drive tumorigenesis.

Oesophageal adenocarcinoma (EAC) incidence is rapidly increasing in Western countries. A better understanding of EAC underpins efforts to improve early detection and treatment outcomes. While large EAC exome sequencing efforts to date have found recurrent loss-of-function mutations, oncogenic driving events have been underrepresented. Here we use a combination of whole-genome sequencing (WGS) and single-nucleotide polymorphism-array profiling to show that genomic catastrophes are frequent in EAC, with almost a third (32%, n=40/123) undergoing chromothriptic events. WGS of 22 EAC cases show that catastrophes may lead to oncogene amplification through chromothripsis-derived double-minute chromosome formation (MYC and MDM2) or breakage-fusion-bridge (KRAS, MDM2 and RFC3). Telomere shortening is more prominent in EACs bearing localized complex rearrangements. Mutational signature analysis also confirms that extreme genomic instability in EAC can be driven by somatic BRCA2 mutations. These findings suggest that genomic catastrophes have a significant role in the malignant transformation of EAC.

Study EGAS00001000750

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

RNA_seq_of_Toxoplasma_gondii_response_in_human_macrophages

Toxoplasmosis is a zoonotic disease caused by a ubiquitous protozoan parasite called Toxoplasma gondii, which caninfect all mammal and bird species throughout the world. seroprevalence varies widely between countries. Studieshave estimated that between 7-34% of people in the UK have been infected with T. gondii. The vast majority of thesepeople will not have noticed any symptoms, however about 10% of people develop a mild to moderate self limitingflu-like illness. Following the acute active stage of the infection the parasite persists in the body in the form of cysts,particularly in heart and skeletal muscle and nervous system tissues, for many years, and usually for life. Inimmunocompetent persons these cysts do not pose a health risk. We will use RNA-seq to quantify the transcriptionalresponse of macrophages to T gondii infection.This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001001708

Novel mutational mechanisms and drivers in Pancreatic Neuroendocrine Tumours

Pancreatic neuroendocrine tumours (PanNETs) are increasing in prevalence due to earlier detection, consequently creating challenges in clinical care, as current classification permits unsatisfactory therapeutic stratification. We performed whole genome sequencing of 102 primary PanNETs and defined the genomic events underlying their pathogenesis. We describe the mutational signatures they harbour, including a novel G:C>T:A Base-Excision-Repair-deficiency signature due to MUTYH inactivation. We uncover a larger than expected proportion of germline mutations in clinically sporadic PanNETs. These include previously unreported mutations in DNA repair genes MUTYH, CHEK2, BRCA2, and PALB2, which together with MEN1 and VHL occur in 16% of patients. Somatic mutations (point mutations and structural rearrangements) commonly occurred in genes involved in 4 main pathways: chromatin remodelling, DNA damage repair, activation of mTOR signalling (including novel EWSR1 gene fusions) and telomere maintenance. Aberrations in these mechanisms are associated with subtypes of PanNET with potential clinical application for prognostic and therapeutic stratification.

Study EGAS00001001732

Exome Sequencing of Gastric Cancer

GASTRIC cancer (GC) is a heterogeneous disease with multiple environmental etiologies and alternative pathways of carcinogenesis. Beyond mutations in TP53, alterations in other genes or pathways account only for small subsets of the disease. We performed exome sequencing of 22 GCs and identified mutated genes and pathway alterations previously unreported: genes involved in chromatin modification were commonly mutated. A downstream validation study confirmed frequent inactivating mutations and/or protein deficiency of the SWI/SNF chromatin remodeling complex gene ARID1A in 83% of GCs with microsatellite instability, 73% of GCs with Epstein-Barr virus (EBV) infection, and 11% of non-EBV microsatellite stable GCs. The mutation spectrum of ARID1A differs between molecular subtypes and the mutation prevalence is negatively associated with mutations in TP53. Clinically, ARID1A alterations were associated with better prognosis in a stage-independent manner. These results revealed the genomic landscape and highlighted the importance of chromatin remodeling in the molecular taxonomy of GCs.

Study EGAS00001000153

Complete Genomics paired end sequencing; Ovarian cancer

We sequenced the genomes from a monozygotic twin discordant for schizophrenia and a tumor-normal pair of an ovarian cancer patient. Using whole-genome twin data to discriminate between correctly identified single nucleotide variants (SNVs) and errors a strategy for the accurate detection of SNVs was developed. By applying stringent sequencing quality measures, excluding error-prone regions and selecting SNVs identified by different mapping and variation calling algorithms, error rates were ~37-fold reduced. This enabled us to identify the first discordant SNVs in monozygotic twins using whole-genome sequencing. In addition, by showing that novel SNVs are highly enriched in errors, accurate estimates of the number of novel and rare SNVs occurring in unrelated Caucasian individuals were obtained. Finally, somatic mutations in coding and regulatory sequences were reliably identified in the highly rearranged ovarian tumor. Overall, our data demonstrate that strategies to reduce error rates in whole-genomes are required for disease gene discovery.

Study EGAS00001000158

Exome_sequencing_of_blastic_plasmacytoid_dendritic_cell_neoplasms

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare and aggressive haematological malignancy derived from precursors of plasmacytoid dendritic cells. Due to the rarity of BPDCNs our knowledge of their molecular pathogenesis was until recently confined to observations describing reccurent chromosomal deletions involving chromosomes 5q, 12p, 13q, 6q, 15q and 9. A recent publication went on to delineate the common deleted regions using aCGH and demonstrated that these centred around known tumour suppressor genes including CDKN2A/B (9p21.3), RB1 (12p13.2-14.3), CDKN1B (13q11-q12) and IKZF1 (7p12.2).These mutations are found recurrently in several different cancers and in most cases are thought to be involved in tumour progression rather than initiation. However, the well-defined nature and cellular ontogeny of these neoplasms suggests strongly that they share one or a few characteristic mutations as has been demonstrated for other uncommon but well-defined neoplasms such as Hairy Cell Leukemia (BRAF) and ovarian Granulosa Cell tumours (FOXL2).

Study EGAS00001000171

Association of Age at Diagnosis and Genetic Mutations in Patients with Neuroblastoma

Neuroblastoma is diagnosed over a wide age range from birth through young adulthood, and older age at diagnosis is associated with a decline in survivability. We performed whole genome sequencing of DNA from diagnostic tumors and their matched germlines from 40 patients with metastatic neuroblastoma obtained between 1987 and 2009. To identify genetic lesions that are associated with age at diagnosis in patients with metastatic neuroblastoma. Age groups at diagnosis included infants (0-18 months), children (18 months-12 years), and adolescents and young adults (>12 years). To confirm the findings from this discovery cohort, validation testing using tumors from an additional 64 patients obtained between 1985 and 2009 was also performed. Formalin-fixed paraffin-embedded tumor tissue was used for immunohistochemistry and fluorescent in situ hybridization. Telomere lengths were analyzed using the whole genome sequencing data, quantitative polymerase chain reaction and fluorescent in situ hybridization

Study EGAS00001000213

Exome-sequencing identifies new oncogenes and tumor suppressor genes recurrently altered in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy. High-resolution copy number analysis of 125 tumors of which 24 were subjected to whole-exome sequencing identified 135 homozygous deletions and 994 somatic gene mutations with predicted functional consequences. We identified new recurrent alterations in 4 genes (ARID1A, RPS6KA3, NFE2L2 and IRF2) not previously described in HCC. Functional analyses demonstrated tumor suppressor properties for IRF2 whose inactivation, exclusively found in hepatitis B virus related tumors, leads to impaired TP53 function. Alternatively, inactivation of proteins involved in chromatin remodeling was frequent and predominant in alcohol related tumors. Moreover, activation of the oxidative stress metabolism and inactivation of RPS6KA3 were new pathways associated with WNT/beta-catenin activation, thereby suggesting a cooperative effect in tumorigenesis. This study shows the dramatic somatic genetic diversity in HCC, it reveals interactions between oncogene and tumor suppressor gene mutations markedly related to specific risk factors.

Study EGAS00001000217

Genome Landscape of High-Grade Serous Ovarian Cancer

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.

Study EGAS00001000397

A Genomics-Based Classification of Human Lung Tumors

We characterized genome alterations in 1255 clinically annotated lung tumors of all histological subgroups to identify genetically defined and clinically relevant subtypes. More than 55% of all cases had at least one oncogenic genome alteration potentially amenable to specific therapeutic intervention, including several personalized treatment approaches that are already in clinical evaluation. Marked differences in the pattern of genomic alterations existed between and within histological subtypes, thus challenging the original histomorphological diagnosis. Immunohistochemical studies confirmed many of these reassigned subtypes. The reassignment eliminated almost all cases of large cell carcinomas, some of which had therapeutically relevant alterations. Prospective testing of our genomics-based diagnostic algorithm in 5145 lung cancer patients enabled a genome-based diagnosis in 3863 (75%) patients, confirmed the feasibility of rational reassignments of large cell lung cancer, and led to improvement in overall survival in patients with EGFR-mutant or ALK-rearranged cancers. Thus, our findings provide support for broad implementation of genome-based diagnosis of lung cancer.

Study EGAS00001000647

High number of somatic mutations found in the healthy blood compartment of a 115-year-old woman reveals oligoclonal hematopoiesis

The somatic mutation burden in healthy white blood cells (WBCs) is not well known. Based on deep whole genome sequencing, we estimate that ~450 somatic mutations accumulated in the non-repetitive genome within the healthy blood compartment of a 115-year old woman. The detected mutations appear to have been harmless passenger mutations: they were enriched in non-coding, AT-rich regions that are not evolutionarily conserved, and they were depleted for genomic elements where mutations might have favorable or adverse effects on cellular fitness, such as regions with actively transcribed genes.The distribution of variant allele frequencies of these mutations suggests that the majority of the peripheral white blood cells were offspring of two related hematopoietic stem cell (HSC) clones. Moreover, telomere lengths of the WBCs were significantly shorter than telomere lengths from other tissues. Together, this suggests that the finite lifespan of HSCs, rather than somatic mutation effects, may lead to hematopoietic clonal evolution at extreme ages.

Study EGAS00001000660

Germline and somatic SMARCA4 mutations characterize small-cell carcinoma of the ovary, hypercalcemic type.

Small-cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is the most common undifferentiated ovarian malignancy in women under 40 years1. We sequenced the exomes of six individuals from three SCCOHT families. After discovering segregating deleterious germline mutations in SMARCA4 in all three families, we tested DNA from a fourth family, which also possessed a segregating SMARCA4 germline mutation. All familial tumors sequenced harboured either a somatic mutation or loss of the wild-type allele. Immunohistochemical (IHC) analysis of these and additional familial and non-familial cases showed loss of the SMARCA4 protein (BRG1) in 38/40 tumors overall. Sequencing of cases with available DNA revealed at least one germline or somatic deleterious SMARCA4 mutation in 30/32 cases. Additionally, the SCCOHT cell line BIN-67 contained bi-allelic deleterious mutations in SMARCA4. Our findings identify alterations in SMARCA4 as the major cause of SCCOHT, which could lead to genetic counseling and new treatment approaches.

Study EGAS00001000721

Integrated genetic and epigenetic analysis defines novel molecular clusters in rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in childhood. To unravel the genetic/epigenetic basis of RMS, we studied 60 RMS cases using whole exome/transcriptome sequencing, copy number (CN) profiling, and DNA methylome analysis. Based on DNA methylation patterns, RMS was clustered into 4 distinct subtypes, which exhibited remarkable correlation with mutation/CN profiles, histological phenotypes, and clinical behaviors. A1 and A2 subtypes, especially A1, largely corresponded to alveolar histology with frequent PAX3/7 fusions and alterations in cell cycle regulators. In contrast, mostly showing embryonal histology, both E1 and E2 subtypes were characterized by high frequency of CN alterations and/or allelic imbalances, FGFR4/RAS/AKT pathway mutations and PTEN mutations/methylation and in E2, also by p53 inactivation. Despite generally poor prognosis in RMS, patients in E1 cluster exceptionally showed excellent prognosis. Our results highlight the close relationship between DNA methylation status and gene mutations and biological behavior in RMS.

Study EGAS00001000884

65 prostate cancer cases WGS and transcriptome sequencing project

We obtained whole-genome and transcriptome sequencing of tumor-benign pairs from 65 Asian prostate cancer (PCa) patients. Taken together with the TCGA PCa cohort, our integration analysis provides the first comprehensive mutational landscape of PCas across major ethnics worldwide and reveals potential novel role of frequent altered chromatin remodeler genes in PCas without androgen receptor (AR) alteration and complex rearrangements derived from Chromoplexy. In addition, we identified 5 putative tumor suppressor genes (TSGs) chains including one on chromosome 13q with novel findings of PCDH9 (deletions in 15/65 tumors) which was implicated in aggressive progression of prostate cancer. Finally, the integrative pathway and co-expression network analysis followed by the functional study indicated the frequent altered genes, particularly focal amplifications of PLXNA1 (11/65), in semaphoring signaling of axon guidance pathway may interact with PI3K/AKT and AR signaling and also contribute to the aggressiveness and progression of prostate cancer.

Study EGAS00001000888

SeqControl: Process Control for DNA Sequencing

As next-generation sequencing (NGS) continues to increase in speed and throughput, routine clinical and industrial application draws steadily closer. These “production” uses of NGS will require enhanced quality-monitoring and quality-control to optimize output and reduce costs. We therefore developed a framework called SeqControl for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, base-wise coverage and base-wise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references we derive multivariate models that predict sequencing quality and depth. SeqControl robustly predicts how much sequencing is required to reach a given coverage depth (AUC = 0.993), accurately classifies clinically relevant formalin-fixed paraffin-embedded samples and makes predictions from as little as 1/8 of a lane of sequencing data (AUC = 0.967). These techniques can be immediately incorporated into existing NGS pipelines to monitor data quality in real-time. SeqControl represents a first step towards statistical process-control for NGS.

Study EGAS00001000899

Genetic_profiling_of_mucosal_melanoma

Primary mucosal melanomas (MMs) arise from melanocytes located in mucosal membranes lining the respiratory, gastrointestinal and urogenital tracts. MMs frequently present late and have a poor prognosis; the 5-year survival rate is only 14%. MM makes up only ~1.4% of all melanomas and it is this rarity that makes knowledge of the genetic changes that contribute to its pathogenesis limited to a small number of exome/genome studies and other targeted studies. Thus to investigate the somatic alterations and mutation spectra in MM genomes, we have extracted genomic DNA from formalin-fixed, paraffin-embedded (FFPE) human MMs, and subjected them to whole exome sequencing. Given the propensity of MM to metastasize, we will also be sequencing metastatic MM lesions; primary and metastatic lesions from the same individual represent an excellent opportunity to identify potential drivers of metastasis in MM. Finally we will sequence ‘normal’ DNA from the same individual, where possible, to exclude germline variations.

Study EGAS00001001115

Exome sequencing of hepatocellular carcinomas identifies new mutational signatures and potential therapeutic targets

Genomic analyses promise to improve tumor characterization in order to optimize personalized treatment for patients with hepatocellular carcinoma (HCC). Exome sequencing analysis of 243 liver tumors revealed mutational signatures associated with specific risk factors, mainly combined alcohol/tobacco consumption, aflatoxin B1 and HBV infection. We identified 161 putative driver genes associated with 11 recurrent pathways. Associations of mutations defined 3 groups of genes related to risk factors and centered on CTNNB1 (alcohol), TP53 (HBV) and AXIN1. Analyses according to tumor stage progression revealed TERT promoter mutation as an early event whereas FGF/CCND1 amplification, TP53 and CDKN2A alterations appeared at more advanced stages in aggressive tumors. In 24% of the tumors, we identified genetic alterations potentially targetable by FDA-approved drugs. In conclusion, we identified risk factor–specific mutational signatures and defined the extensive landscape of altered genes and pathways in HCC, which will be useful to design clinical trials for targeted therapy.

Study EGAS00001001002

Whole-genome plasma sequencing reveals focal amplifications as a driving force in metastatic prostate cancer

Genomic alterations in metastatic prostate cancer remain incompletely characterized. Here we analyze 493 prostate cancer cases from the TCGA database and perform whole-genome plasma sequencing on 95 plasma samples derived from 43 patients with metastatic prostate cancer. From these samples, we identify established driver aberrations in a cancer-related gene in nearly all cases (97.7%), including driver gene fusions (TMPRSS2:ERG), driver focal deletions (PTEN, RYBP, SHQ1), and driver amplifications (AR, MYC). In serial plasma analyses, we observe changes in focal amplifications in 40% of cases. The mean time interval between new amplifications was 26.4 weeks (range: 5-52 weeks), suggesting that they represent rapid adaptations to selection pressure. An increase in neuron-specific enolase is accompanied by clonal pattern changes in the tumor genome, most consistent with subclonal diversification of the tumor. Our findings suggest a high plasticity of prostate cancer genomes with newly occurring focal amplifications as a driving force in progression

Study EGAS00001001018

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

Recurrent mTORC1-activating RRAGC mutations in follicular lymphoma

Follicular lymphoma (FL) is an incurable B-cell malignancy characterized by the t(14;18) and mutations in one or more components of the epigenome. Whilst frequent gene mutations in signaling pathways, including JAK-STAT, NOTCH and NF-κB, have also been defined, the spectrum of these mutations typically overlap with the closely-related diffuse large B cell lymphoma (DLBCL). A combination of discovery exome and extended targeted sequencing revealed recurrent somatic mutations in RRAGC uniquely enriched in FL patients (17%). More than half of the mutations preferentially co-occurred with ATP6V1B2 and ATP6AP1 mutations, components of the vacuolar H+-adenosine triphosphate ATPase (v-ATPase) known to be necessary for amino acid-dependent mTORC1 activation. RagC mutants increased raptor binding whilst rendering mTORC1 signaling resistant to amino acid deprivation through nucleotide-independent mechanisms. Collectively, the activating nature of the RRAGC mutations, their existence within the dominant clone and stability during disease progression supports their potential as an excellent candidate to be therapeutically exploited.

Study EGAS00001001190

Preeclampsia InterPregGen Consortium: Large-scale GWAS meta-analysis of maternal preclampsia cases and controls from Europe and Central Asia. Plus extension of earlier European fetal preeclampsia GWAS meta-analysis (see EGAS00001001048) by adding Central Asian fetal preeclampsia cases and controls. Datasets provided under this study are GWAS meta-analysis summary statistics and individual level GWAS genotype data. Related InterPregGen Consortium data are also provided under studies EGAS00001000416 and EGAS00001000417 (whole genome sequence data for 100 unrelated Uzbeks and 100 unrelated Kazakhs) and EGAS00001001048 (European fetal preeclampsia GWAS summary statistics and genotype data

Preeclampsia (PE) is a syndrome affecting pregnant mothers and fetus/babies characterised by hypertension and proteinuria, and is a leading cause of maternal and fetal death and of premature births worldwide. The InterPregGen Consortium was funded by a European Framework 7 (FP7) grant and grew out of the WTCCC3 GWAS comparing ~2000 UK PE mothers with ~6000 common UK controls. The EGAS00001001266 study includes two components the FP7 InterPregGen project: (a) GWAS case-control meta-analysis of maternal PE cases from Europe (UK, Iceland, Norway, Denmark, Finland) and Central Asia (Kazakhstan, Uzbekistan); (b) GWAS case-control meta-analysis of PE fetal (baby) cases from Europe (UK, Iceland, Norway, Denmark) and Central Asia (Kazakhstan, Uzbekistan)

Study EGAS00001001266

Recurrent adeno-associated virus 2-related insertional mutagenesis in human hepatocellular carcinomas

Hepatocellular carcinomas (HCC) are liver tumors related to various etiologies including alcohol intake, hepatitis B (HBV) or C (HCV) virus infection. Additional risk factors remain to be identified, particularly in patients who develop HCC without cirrhosis. We identified clonal integration of adeno-associated virus type 2 (AAV2) in 11 out of 193 HCC. These AAV2 integrations occurred within known cancer driver genes, namely, CCNA2 (Cyclin A2, 4 cases), TERT (Telomerase Reverse-Transcriptase, 1 case), CCNE1 (Cyclin E1, 3 cases), TNFSF10 (Tumor Necrosis Factor member 10, 2 cases) or KMT2B (Lysine (K)-Specific Methyltransferase 2B, 1 case) leading to over-expression of the target genes. Tumors with viral integration mainly developed in non-cirrhotic liver (9 out of 11 cases) and without known risk factors (6 out of 11) suggesting a pathogenic role of AAV2 in these patients. In conclusion, AAV2 is a DNA virus associated with oncogenic insertional mutagenesis in human HCC.

Study EGAS00001001284

Barcoding reveals complex clonal dynamics of de novo transformed human mammary cells

Most human breast cancers have already diversified genomically when they first become clinically evident, by which time extensive heterogeneous histopathologies, transcriptomes and growth patterns are also apparent. Accordingly, important initial events and the cellular context in which they occur have been difficult to characterize. Using DNA barcoding, we now demonstrate the high efficiency with which both purified basal and luminal cells isolated directly from normal adult human mammary tissue can be rapidly transformed by a single oncogene (KRASG12D) resulting in the production within 8 weeks in vivo of serially transplantable, polyclonal, invasive ductal carcinomas that are phenotypically heterogeneous and transcriptionally distinct from the initial cells transduced. Barcoding also revealed a consistent dramatic change in the clonal content of passaged tumours. This system thus provides a powerful new platform for examining early events in the genesis, evolution and treatment response of malignant human mammary cells generated using defined mutations.

Study EGAS00001001310

A comprehensive assessment of somatic mutation detection in cancer using whole genome sequencing

As whole genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Using tumor-normal sample pairs from two different types of cancer, chronic lymphocytic leukemia and medulloblastoma, we conducted a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines, and even validation methods. Here we show that using PCR-free methods and increasing sequencing depth to ~100x showed benefits, as long as the tumor:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artifact-prone nature of the raw data and lack of standards for dealing with the artifacts. However, armed with the benchmark mutation set we have created, we show that many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Study EGAS00001001539

Integrated genomic characterization of IDH1 mutant Glioma malignant progression

Gliomas represent approximately 30% of all central nervous system tumors, and 80% of malignant brain tumors. To understand the molecular mechanisms underlying the malignant progression of low-grade Isocitrate Dehydrogenase 1 (IDH1) mutant gliomas, we studied paired tumor samples from 41 patients, comparing higher-grade, progressed samples to their lower-grade counterparts. Integrated genomic analyses, including whole-exome sequencing, copy number, gene expression and DNA methylation profiling, demonstrated non-linear clonal expansion of the original tumors and identified oncogenic pathways driving progression. These include activation of the MYC and RTK-RAS-PI3K pathways, up-regulation of the FOXM1- and E2F2-mediated cell cycle transition genes, as well as epigenetic silencing of developmental transcriptional factors bound by Polycomb repressive complex 2 in human embryonic stem cells. Our results not only provide mechanistic insight into the genetic and epigenetic mechanisms driving glioma progression but also identify inhibition of the bromodomain and extra-terminal (BET) family as a potential therapeutic approach.

Study EGAS00001001588

Spatial and Temporal Homogeneity of Driver Mutations in Diffuse Intrinsic Pointine Glioma

Diffuse Intrinsic Pontine Gliomas (DIPG) are deadly pediatric brain tumors where needle biopsies help guide diagnosis and targeted therapies. To address spatial heterogeneity, we analyzed 134 specimens from various neuroanatomical structures of whole autopsy brains from nine DIPG patients. Evolutionary reconstruction indicates histone 3 (H3)K27M – including novel H3.2K27M - mutations potentially arise first and are invariably associated with specific, high-fidelity obligate partners throughout the tumor and its spread, from diagnosis to end-stage disease, suggesting mutual need for tumorigenesis. These H3K27M ubiquitously-associated mutations involve alterations in TP53 cell-cycle (TP53/PPM1D) or specific growth factor pathways (ACVR1/PIK3R1). Later oncogenic alterations arise in sub-clones and often affect the PI3K pathway. Our findings are consistent with early tumor spread outside the brainstem including the cerebrum. The spatial and temporal homogeneity of driver mutations in DIPG implies they will be captured by limited biopsies and emphasizes the need to develop therapies specifically targeting obligate oncohistone partnerships.

Study EGAS00001001654

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

Multi-layered population structure in Island Southeast Asians

The history of human settlement in Southeast Asia has been complex and involved several distinct dispersal events. Here we report the analyses of 1825 individuals from Southeast Asia including new genome-wide genotype data for 146 individuals from three Mainland Southeast Asian (Burmese, Malay and Vietnamese) and four Island Southeast Asian (Dusun, Filipino, Kankanaey and Murut) populations. While confirming the presence of previously recognized major ancestry components in the Southeast Asian population structure, we highlight the Kankanaey Igorots from the highlands of the Philippine Mountain Province as likely the closest living representatives of the source population that may have given rise to the Austronesianexpansion. This conclusion rests on independent evidence from various analyses of autosomal data and uniparental markers. Given the extensive presence of trade goods, cultural and linguistic evidence of Indian influence in Southeast Asia starting from 2.5kya we also detect traces of a South Asian signature in different populations in the region dating to the last couple of thousand years.

Study EGAS00001001738

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

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

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

Study EGAS00001001838

whole-genome sequencing in 17 ESCC cases and whole-exome sequencing in 71 cases

Esophageal cancer is one of the most aggressive cancers and the sixth leading cause of cancer death worldwide1. Approximately 70% of the global esophageal cancers occur in China and over 90% histopathological forms of this disease are esophageal squamous cell carcinoma (ESCC)2-3. Currently, there are limited clinical approaches for early diagnosis and treatment for ESCC, resulting in a 10% 5-year survival rate for the patients. Meanwhile, the full repertoire of genomic events leading to the pathogenesis of ESCC remains unclear. Here we show a comprehensive genomic analysis in 158 ESCC cases, as part of the International Cancer Genome Consortium (ICGC) Research Projects (http://icgc.org/icgc/cgp/72/371/1001734). We conducted whole-genome sequencing in 17 ESCC cases and whole-exome sequencing in 71 cases, of which 53 cases and additional 70 ESCC cases were subjected to array comparative genomic hybridization (a-CGH) analysis.

Study EGAS00001000709

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

Variation in the Glucose Transporter gene SLC2A2 is associated with glycaemic response to metformin

Metformin is the first-line antidiabetic drug with over 100 million users worldwide, yet its mechanism of action remains unclear1. Here the Metformin Genetics (MetGen) Consortium reports a three-stage genome wide association study (GWAS), consisting of 13,123 participants of different ancestries. The C-allele of rs8192675 in the intron of SLC2A2, which encodes the facilitated glucose transporter GLUT2, was associated with a 0.17% (p=6.6x10-14) greater metformin induced HbA1c reduction in 10,577 participants of European ancestry. rs8192675 is the top cis-eQTL for SLC2A2 in 1,226 human liver samples, suggesting a key role for hepatic GLUT2 in regulation of metformin action. In obese individuals C-allele homozygotes at rs8192675 had a 0.33% (3.6mmol/mol) greater absolute HbA1c reduction than T-allele homozygotes.This is about half the effect seen with the addition of a DPP-4 inhibitor, and equates to a dose difference of 550mg of metformin, suggesting rs8192675 as a potential biomarker for stratified medicine.

Study EGAS00001001875

Hi-C dataset for testicular germ cell tumour GWAS risk loci, as described in the Oncoarray Litchfield et al. 2016 paper.

Genome-wide association studies (GWAS) have transformed our understanding of testicular germ cell tumour (TGCT) susceptibility but much of the heritability remains unexplained. Here we report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 7,319 TGCT cases and 23,082 controls. We identify 19 new TGCT risk loci, approximately doubling the number of known TGCT risk loci to 44. By performing in-situ Hi-C in TGCT cells, we establish a network of physical interactions between all 44 TGCT risk SNPs and candidate causal genes. Our findings reveal widespread disruption of developmental transcriptional regulators as a basis of disease susceptibility, consistent with failed primordial germ cell differentiation as an initiating step in TGCT oncogenesis1. Defective microtubule assembly and dysregulation of KIT-MAPK signalling also feature as recurrently disrupted pathways. Our findings support a polygenic model of disease risk and provide insight into the biological basis of TGCT.

Study EGAS00001001930

New Brain Tumor Entities Emerge from Molecular Classification of CNS-PNETs

Primitive neuroectodermal tumors of the central nervous system (CNS-PNETs) are highly aggressive, poorly differentiated embryonal tumors occurring predominantly in young children but also affecting adolescents and adults. Herein we demonstrate that a significant proportion of institutionally diagnosed CNS-PNETs display molecular profiles indistinguishable from those of various other well defined CNS tumor entities, facilitating diagnosis and appropriate therapy for patients with these tumors. From the remaining fraction of CNS-PNETs we identify four new CNS tumor entities, each associated with a recurrent genetic alteration and distinct histopathological and clinical features. These new molecular entities, designated “CNS neuroblastoma with FOXR2 activation (CNS NB-FOXR2)”, “CNS Ewing sarcoma family tumor with CIC alteration (CNS EFT-CIC)”, “CNS high grade neuroepithelial tumor with MN1 alteration (CNS HGNET-MN1)”, and “CNS high grade neuroepithelial tumor with BCOR alteration (CNS HGNET-BCOR)”, will enable meaningful clinical trials and the development of therapeutic strategies for patients affected by poorly differentiated CNS tumors.

Study EGAS00001001632

A living biobank of breast cancer organoids captures disease heterogeneity

Breast Cancer (BC) comprises multiple distinct subtypes that differ genetically, pathologically, and clinically. Xeno-transplantation and direct culturing of tumor tissue are increasingly used in drug development and personalized medicine strategies. Here, we describe a robust protocol for long-term culturing of human mammary epithelial organoids. Using this protocol, 101 BC organoid lines were generated, broadly recapitulating the diversity of the disease. BC organoid morphologies typically matched the histopathology, hormone receptor-, and HER2 status of the original tumor. DNA copy number variations as well as sequence changes were consistent within tumor-organoid pairs and largely retained even after extended passaging. BC organoids furthermore populated all major gene expression-based classification groups and allowed drug screens in vitro and upon xenotransplantation. This study describes a representative collection of well-characterized BC organoids available for cancer research and drug development, as well as a strategy to assess in vitro drug response in a personalized fashion.

Study EGAS00001002158

The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

Dataset EGAD00001007033

Fetal origins of malignant germ cell tumours

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells (PGCs). They are uniquely able to generate embryonic and extraembryonic tissues, which in malignant GCTs carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we studied the phylogenetic and transcriptional diversity of 15 malignant gonadal GCTs and four normal testis biopsies by sequencing 131 whole genomes and 416 transcriptomes from 14 gonadal histologies, excised by laser capture microdissection. Our findings demonstrate that tumours were initiated by whole genome duplication likely in embryogenesis, within ~5-8 cell divisions post-PGC specification, followed by chromosome 12p gains associated with invasive disease. Of note, 12p imbalances were not only generated through GCT-typical isochromosomes, but also through non-isochromosomic configurations. Whilst tumours developed along homogenous phylogenetic pathways, they spawned manifold tissues independent of genetic subclonal diversification. A key feature of GCT tissues was the expression of fetal-specific genes. The transcriptional diversity notwithstanding, we found universal transcriptional elements correlated with hallmark 12p gains. Overall, our study reveals stereotyped phylogenies and transcriptomes underpinning the development of GCT that originate in fetal life and may lend themselves to therapeutic manipulation.

Dataset EGAD00001007037

Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - TGS

Mutations in cancer-associated genes drive tumour outgrowth. However, the timing of driver mutations and dynamics of clonal expansion that lead to human cancers are largely unknown. We used 580,133 somatic mutations from whole-genome sequencing of 1013 clonal haematopoietic colonies to reconstruct the phylogeny of haematopoiesis, from embryogenesis to clinical disease, in 12 patients with myeloproliferative neoplasms which are blood cancers more common in older age. JAK2V617F, the pathognomonic mutation driving the majority of these cancers, was acquired in utero or childhood, with upper estimates of age of acquisition from 33 weeks gestation to 10.8 years, in all 5 patients in whom JAK2V617F was either the only or the first driver event. Driver mutations associated with age-related clonal haematopoiesis occurred prior to or following JAK2V617F, as independent clonal expansions in JAK2V617F-mutated patients, and as large clonal expansions in JAK2V617F-unmutated patients . These mutations were also acquired in utero or childhood, with DNMT3A mutations occurring by 8 weeks of gestation to 7.6 years across 4 patients, and PPM1D mutation occurring by age 5.8yrs in a patient with MPN lacking phenotypic driver mutations. Sequential driver mutation acquisition was common, separated by decades across life, and often outcompeted ancestral clones. The mean latency between JAK2V617F acquisition and clinical presentation was 30 years (range 11-54 years). Rates of clonal expansion were inferred from phylogenetic trees and varied substantially (3% to 190% expansion/year), were affected by additional driver mutations, and were predictive of latency to clinical presentation. Driver mutations and rates of expansion would have been detectable in blood one to four decades before clinical presentation. This study reveals how driver mutation acquisition early in life with life-long growth and evolution underlie adult myeloproliferative neoplasms, providing opportunities for early detection and intervention and a new paradigm for cancer development.

Dataset EGAD00001007715

Distribution of ctDNA levels in plasma of early-stage non-small cell lung cancer patients measured using personalised ctDNA analysis

Blood-based assays have shown increasing ability to detect circulating tumour DNA (ctDNA) in patients with early-stage cancer. However, detection of ctDNA in patients with non-small cell lung cancer (NSCLC) has continued to prove challenging. We performed retrospective analysis to quantify ctDNA levels in a cohort of 100 patients with early-stage NSCLC prior to treatment with curative intent. Where tumour tissue was available for whole exome sequencing, mutations identified were used to define patient-specific sequencing assays. For those 90 patients, plasma cell-free DNA was sequenced to high depth across capture panels targeting a median of 328 mutations specific to each patient. Data was analysed using Integration of Variant Reads (INVAR), detecting ctDNA in 66.7% of patients, including 52.7% (29 of 55) patients with stage I disease and >88% detection for patients with stage II and III disease (16/18 and 15/17). ctDNA was detected in plasma at fractional concentrations as low as 9.1x10-6, and in patients with tumour volumes as low as 0.23 cm3. A 36-gene sequencing panel (InVisionFirst-LungTM) was used to analyse plasma DNA in 27 samples including the 10 cases without tumour exome data, and detected ctDNA in 59% of samples tested (16 of 27). Across the entire cohort, detection rates were higher in squamous cell carcinoma patients compared to adenocarcinoma patients (81% vs. 59%). Detection of ctDNA prior to treatment was associated with significantly shorter time free from relapse, across all patients and in patient subgroups, with Hazard Ratios ranging from 2.25 to >11. Our analysis indicates that for patients with stage I NSCLC, the median ctDNA fraction in plasma is approx. 12 parts per million (0.0012%). This indicates the limits of detection that would be required for ctDNA-based liquid biopsies to detect ctDNA in the majority of patients with early-stage NSCLC.

Dataset EGAD00001006230

pancreatic ductal adenocarcinoma exomes in study: Pro-immunogenic Impact of MEK inhibition combined with an anti-CD40 immunostimulatory antibody

This dataset contains 3 pairs of exomes, germline (from whole blood) and patient-derived xenograft (PDX), from human pancreatic durctal adenocarcinoma patients. The data is referred to in the publication: "Pro-immunogenic impact of MEK inhibition synergizes with agonist anti-CD40 immunostimulatory antibodies in tumor therapy" (Nature Communications, 2020) Abstract: Cancer types with lower mutational load and a non-permissive tumor microenvironment are intrinsically resistant to immune checkpoint blockade. While the combination of cytostatic drugs and immunostimulatory antibodies constitutes an attractive concept for overcoming this refractoriness, suppression of immune cell function by cytostatic drugs may limit therapeutic efficacy. Here we show that targeted inhibition of mitogen-activated protein kinase (MAPK) kinase (MEK) does not impair dendritic cell-mediated T-cell priming and activation. Accordingly, combining MEK inhibitors (MEKi) with agonist antibodies (Abs) targeting the immunostimulatory CD40 receptor resulted in potent synergistic anti-tumor efficacy. Detailed analysis of the mechanism of action of MEKi GDC-0623 by means of flow cytometric analysis of the tumor immune infiltrate and whole tumor transcriptomics showed that, in addition to its cytostatic impact on tumor cells, this drug exerts multiple pro-immunogenic effects, including the suppression of M2-type macrophages, myeloid derived suppressor cells and CD4+ T-regulatory cells. In addition, MEKi was found to induce tumor-cell intrinsic interferon signaling, which contributed to antigen presentation by tumor cells. Finally, the tumoridical impact of MEKi involves the activation of multiple pro-inflammatory pathways involved in immune cell effector function in the tumor microenvironment. Our data therefore indicate that the combination of MEK inhibition with agonist anti-CD40 Ab is a promising therapeutic concept, especially for the treatment of mutant Kras-driven tumors such as pancreatic ductal adenocarcinoma.

Dataset EGAD00001005938

Preferential infiltration of unique Vγ9Jγ2‐Vδ2 T cells into glioblastoma multiforme

Glioblastoma multiforme (GBM) is clinically highly aggressive as a result of evolutionary dynamics induced by cross-talk between cancer cells and a heterogeneous group of immune cells in tumor microenvironment. The brain harbors limited numbers of immune cells with few lymphocytes and macrophages; thus, innate‐like lymphocytes, such as γδ T cells, have important roles in antitumor immunity. Here, we characterized GBM‐infiltrating γδ T cells, which may have roles in regulating the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor‐infiltrating and blood‐circulating γδ T cells from four patients were analyzed by next-generation sequencing-based T-cell receptor (TCR) sequencing in addition to mutation and immune profiles in four GBM cases. In all tumor tissues, abundant innate and effector/memory lymphocytes were detected, accompanied by large numbers of tumor‐associated macrophages and closely located tumor‐infiltrating γδ T cells, which appear to have anti-tumor activity. The immune-related gene expression analysis using the TCGA database showed that the signature gene expression extent of γδ T cells were more associated with those of cytotoxic T and Th1 cells and M1 macrophages than those of Th2 cells and M2 macrophages. Although the most abundant γδ T cells were Vγ9Vδ2 T cells in both tumor tissues and blood, the repertoire of intratumoral Vγ9Vδ2 T cells was distinct from that of peripheral blood Vγ9Vδ2 T cells and was dominated by Vγ9Jγ2 sequences, not by canonical Vγ9JγP sequences that are mostly commonly found in blood γδ T cells. Collectively, unique GBM‐specific TCR clonotypes were identified by comparing TCR repertoires of peripheral blood and intra‐tumoral γδ T cells. These findings will be helpful for the elucidation of tumor-specific antigens and development of anticancer immunotherapies using tumor-infiltrating γδ T cells.

Dataset EGAD00001004862

Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - WGS

Mutations in cancer-associated genes drive tumour outgrowth. However, the timing of driver mutations and dynamics of clonal expansion that lead to human cancers are largely unknown. We used 580,133 somatic mutations from whole-genome sequencing of 1013 clonal haematopoietic colonies to reconstruct the phylogeny of haematopoiesis, from embryogenesis to clinical disease, in 12 patients with myeloproliferative neoplasms which are blood cancers more common in older age. JAK2V617F, the pathognomonic mutation driving the majority of these cancers, was acquired in utero or childhood, with upper estimates of age of acquisition from 33 weeks gestation to 10.8 years, in all 5 patients in whom JAK2V617F was either the only or the first driver event. Driver mutations associated with age-related clonal haematopoiesis occurred prior to or following JAK2V617F, as independent clonal expansions in JAK2V617F-mutated patients, and as large clonal expansions in JAK2V617F-unmutated patients . These mutations were also acquired in utero or childhood, with DNMT3A mutations occurring by 8 weeks of gestation to 7.6 years across 4 patients, and PPM1D mutation occurring by age 5.8yrs in a patient with MPN lacking phenotypic driver mutations. Sequential driver mutation acquisition was common, separated by decades across life, and often outcompeted ancestral clones. The mean latency between JAK2V617F acquisition and clinical presentation was 30 years (range 11-54 years). Rates of clonal expansion were inferred from phylogenetic trees and varied substantially (3% to 190% expansion/year), were affected by additional driver mutations, and were predictive of latency to clinical presentation. Driver mutations and rates of expansion would have been detectable in blood one to four decades before clinical presentation. This study reveals how driver mutation acquisition early in life with life-long growth and evolution underlie adult myeloproliferative neoplasms, providing opportunities for early detection and intervention and a new paradigm for cancer development.

Dataset EGAD00001007714

The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

Dataset EGAD00001007034

Transcriptional profiling of tauopathies in human IPS-derived neurons (2019-08-21)

This project is a pilot study, in collaboration with Maria Grazia Spillantini and Mariangela Iovino (Cambridge Centre for Brain Repair), to investigate the utility of IPS-derived neurons for the study of neurodegenerative disorders. Our aim is to characterise the transcriptional consequences of tauopathies using neurons derived from differentiated IPSCs as a model system. We will use IPS cells derived from six individuals, four with known mutations in the tau protein, 2 without. RNA will be extracted at Day 0 and Day 65 of differentiation by which time the neuronal tauopathy is apparent. RNA will be extracted and the transcriptome of each line characterised using RNAseq. We will then search for genes that are differentially expressed between the transcriptomes of individuals with tau mutations versus those in controls. My lab will analyse the RNAseq data, comparing both affected and controls and both time-points, to establish candidate genes. Darren Logan’s lab, along with our collaborators, will experimentally verify and further investigate these genes in additional lines and animal models. From this analysis we will generate a list of candidate genes that are differentially expressed between cases and controls. This study will not only help us understand the molecular basis of tauopathies, but also identify gene candidates for biomarkers of neurodegenerative disease. It will serve as a proof of principle for future planned studies into generating and transcriptomically analysing an allelic series of Tau mutations in IPSCs with a controlled genetic background. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2019-08-21.

Dataset EGAD00001005277

Molecular diagnosis of albinism (2018-03-14)

Albinism is genetically heterogeneous rare genetic condition affecting 1:17000 in the Western world (but more frequent in Africa) whose main feature is a profound visual impairment, characterised by foveal hypoplasia, abnormal chiasmatic connections, nystagmus and photofobia. All these features result in severly altered visual acuity (<0,1), absent depth perception and poor night vision. People with albinism are primarily visually handicapped. In addition, for some types of albinism, the visual phenotype can be presented with partial or total hypopigmentation, hence resulting in a secondary phenotype which can lead to skin cancer if skin is not adequately protected. Recently a new syndrome has been described, FHONDA, with the same visual abnormalities of albinism but without pigment alteration. The traditional classification differentiates Oculoculatenous albinism (OCA), where hypopigmentation involves hair, skin and eyes versus Ocular Albinism (OA), where hypopigmentation only affects the eyes. These are non-sydrimic types of albinism. Some syndromic forms (Hermansky-Pudlak=HPS, Chediak-Higashi=CHS) affect cells beyond pigment cells, present in the lungs, immune system, platelets and intestines, resulting in more severe phenotypes that can be fatal. Mutations in at least 19 genes are assocaited with the corresponding types of albinism. Most hospitals will only diagnose the most frequent cases using traditional Sanger, MLPA approaches. Some will use CGH arrays. We aim to diagnose all cases of albinism through the Albinochip proposal, which combines a Sequenom first step of known mutations combined with subsequent NGS approaches. In some cases we fail to find a second mutation, these are good candidates for further full exome analyses. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2018-03-14.

Dataset EGAD00001004039

Natural history of clonal haematopoiesis (2017-09-04)

The incidence of acute myeloid leukemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 60. Only 10-15% of cases evolve from a pre-existing myeloproliferative or myelodysplastic disorder; the remaining cases arise de novo without a detectable prodrome and are diagnosed upon development of bone marrow failure. Analysis of diagnostic blood samples has demonstrated that de novo AML is preceded by the accumulation of somatic mutations in pre-leukemic hematopoietic stem and progenitor cells (preL-HSPCs) that subsequently undergo clonal expansion. If individuals in this pre-leukemic phase could be identified, methods for determination of risk and monitoring for progression to overt AML could be developed. However recurrent AML mutations also accumulate during aging in healthy individuals who never develop AML, referred to as age related clonal hematopoiesis (ARCH). To distinguish individuals with preL-HSPCs at high risk of developing AML from those with ARCH, we undertook deep targeted sequencing of genes recurrently mutated in AML in blood samples from 133 individuals in the European Prospective Investigation into Cancer and Nutrition (EPIC) study taken on average 6 years before they developed AML (pre-AML group), together with 683 matched healthy individuals (Control group). Pre-AML cases displayed accelerated age-correlated accumulation of somatic mutations.The identity, number and variant allele frequency (VAF) of mutations differed between the two groups, and were incorporated into a computational model of AML risk prediction that accurately distinguished pre-AML cases from controls on average 7 years prior to AML development. Our findings provide proof of concept that early prediction of AML development is feasible in high-risk populations, paving the way for early disease detection, monitoring, and potentially prevention.

Dataset EGAD00001003703

Genetic Heterogeneity of the familial gastric neuroendocrine tumors (2018-06-06)

Gastric neuroendocrine tumors (gNETs) occur with an estimated frequency of 2 per 100,000 in the general population. Type I gastric neuroendocrine tumors (NETs) represent the 75% of gNTEs and arise from gastric enterochromaffin-like (ECL) cells. They have late age of onset and usually benigh course. Classically, hypergastrinemia in patients who have autoimmune atrophic gastritis, causes hyperplasia of gastric ECL cells that progresses into type I gastric NETs and parietal cell (PC) destruction. The genetic bases in families with this disease are unknown. We performed an exome sequencing study of an atypical aggressive familial gNETs case (with early age onset, nodal infiltrations and gastric adenocarcinomas) that followed a recessive model. We identified a deleterious mutation in homozygosis in the ATP4A gene, which encodes the proton pump responsible for acid secretion by gastric parietal cells. This mutation lead to achlorhydria first, and hypergastrinemia and gNET developing as consequence (Calvete et al. 2014). Recently, two more families with gNETs, classical clinical traits and recessive model have been studies by WES but we didn't find any mutation in the ATP4a gene. However, putative mutations affecting genes that contribute to the development and the integrity of PC have been found suggesting that genetic alterations associated to this disorder target to a unique cell type (parietal cells). In order to cinfirm this hypothesis, it is necessary the search for new genes implicated in the gNETs, more familial cases are needed to be studied. We have identified four more new familial gNETs cases. Here, we propose their study by WES. The first family is formed by thress siblings with gNETs. The other families include two siblings with gNETs. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2018-06-06.

Dataset EGAD00001004153

Spiradenocarcinoma (2018-10-29)

Spiradenocarcinoma is a rare cutaneous sweat gland adnexal cancer with potential for aggressive behaviour. They are classified histologically into low- and high-grade tumours, with morphologically low-grade tumours thought to behave more favourably. However, limited information is available, with only 18 published cases. We have collected morphologically low-grade spiroadenocarcinomas (one with a lung metastasis) and high-grade spiroadenocarcinomas, as well as some spiradenomas (benign lesions), cylindromas (another type of malignant cutaneous sweat gland adnexal tumour) and hybrid spiradenoma-cylindromas. H&E-stained sections were reviewed, follow-up was obtained, and immunohistochemistry for Ki-67, p53 and, MYB has been performed. The tumours were solitary, measuring 0.8-7?cm (median: 2.7?cm), with a predilection for the head and neck of elderly patients (median age: 72 years; range 53-92) without gender bias. Histologically, the tumours were multinodular and located in deep dermis and subcutis. A pre-existing spiradenoma was present in all cases. The malignant component was characterized by expansile growth with loss of the dual cell population, up to moderate cytological atypia and increased mitotic activity (median: 10/10 HPF; range 1-28). Additional findings included squamoid differentiation (n=9), necrosis (n=7), and ulceration (n=5). P53 expression was variable and no significant differences were noted in the benign compared with the malignant parts of the tumours. In contrast, in the malignant components the Ki-67 proliferative index was slightly increased, and MYB expression was lost. Follow-up (median: 67 months; range: 13-132) available for 16 patients (84%) revealed a local recurrence rate of 19% but no metastases or disease-related mortality. Here we wish to exome sequence these cases to define the first genomic landscape for this malignancy. This dataset contains all the data available for this study on 2018-10-29.

Dataset EGAD00001004426

Constitutional and somatic genomic rearrangements coherently restructure chromosome 21 in acute lymphoblastic leukaemia

Changes in gene dosage are a major driver of cancer1, engineered from a finite, but increasingly well annotated, repertoire of mutational mechanisms2-6. These processes operate over levels ranging from individual exons to whole chromosomes, often generating correlated copy number alterations across hundreds of linked genes. An example of the latter is the 2% of childhood acute lymphoblastic leukemia (ALL) characterized by recurrent intrachromosomal amplification of megabase regions of chromosome 21 (iAMP21)7,8 To dissect the interplay between mutational processes and selection on this scale, we used genomic, cytogenetic and transcriptional analysis, coupled with novel bioinformatic approaches, to reconstruct the evolution of iAMP21 ALL. We find that individuals born with the rare constitutional Robertsonian translocation between chromosomes 15 and 21, rob(15;21)(q10;q10)c, have ~2700-fold increased risk of developing iAMP21 ALL compared to the general population. In such cases, amplification is initiated by chromothripsis involving both sister chromatids of the dicentric Robertsonian chromosome. In contrast, sporadic iAMP21 is typically initiated by breakage-fusion-bridge (BFB) events, often followed by chromothripsis or other rearrangements. In both sporadic and iAMP21 in rob(15;21)c individuals, the final stages of amplification frequently involve large-scale duplications of the abnormal chromosome. The end-product is a derivative chromosome 21 or a derivative originating from the rob(15;21)c chromosome, der(15;21), respectively, with gene dosage optimised for leukemic potential, showing constrained copy number levels over multiple linked genes. In summary, the constitutional translocation, rob(15;21)c, predisposes to leukemia through a novel mechanism, namely a propensity to undergo chromothripsis, likely related to its dicentric nature. More generally, our data illustrate that several cancer-specific mutational processes, applied sequentially, can co-ordinate to fashion copy number profiles over large genomic scales, incrementally refining the fitness benefits of aggregated gene dosage changes.

Dataset EGAD00001000658

Determination of cell specific regulatory enhancers in hematopoetic models

The BLUEPRINT project is a large-scale project investigating epigenetic mechanisms involved in blood formation, in health and disease. The human variation workpackage (WP10, led by NS) of the project seeks to characterize the effect of common sequence variation on the epigenome status of a cell. To do this, the project will use highly purified blood cells to minimise "experimental noise" and therefore enhance the power to discover modest effects. Two peripheral blood cell types, the CD14+CD16- monocyte (an important central orchestrator of adaptive immunity and a bridge between innate and adaptive immunity) and the CD65+CD9- neutrophilic granulocyte (the frontline cell for innate immunity) have been selected for this purpose. The two types of cells will be obtained at high purity from adult blood (AB) of 200 healthy males and females, respectively. Cells will be purified by using already validated and fully operational protocols that are based on density gradient centrifugation of the buffy coat obtained from whole blood, followed by magnetic bead-based purification using monoclonal antibodies against Cluster of Differentiation (CD) lineage-specific cell surface markers. Units of 475 ml of AB will be obtained from consenting volunteers of the Cambridge BioResource (CBR), a panel of 10,000 healthy volunteers local to Cambridge who have already consented to participate in biomedical research and of whom biological samples (DNA, plasma, serum) and lifestyle data have been deposited in a repository and database, respectively. We are requesting funding from the Human Diversity project to sequence the genomes of the 200 CBR volunteers at low pass (6x coverage). Nuclei, DNA and RNA will be recovered from the purified cells and made available for RNA-seq, DNA-seq and ChIP-seq and genomic DNA for entire genome sequencing will be recovered from the DNA repository.

Dataset EGAD00001002205

WGS BAM files fromxa0Genomic Features and Classification of Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma

Background and Aims: Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline(g) alterations in BRCA1, BRCA2 and PALB2. Methods: We interrogated whole genome sequencing (WGS) data on 391 patients including 49 carriers of pathogenic variants (PVs) in g_BRCA_ and PALB2. HRD classifiers were applied to the dataset and included: 1) the genomic instability score (GIS) used by Myriad MyChoice HRD assay; 2) substitution base signature 3 (SBS3); 3) HRDetect; and, 4) Structural Variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status. Results: Biallelic tumour inactivation of g_BRCA_ or PALB2 was evident in 43/49 germline carriers identifying HRD-PDAC. HRDetect (score ?0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumour classifiers suggested that 7-10% of PDAC that do not harbor g_BRCA/PALB2_ have features of HRD. Of the somatic HRDetect cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared to BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC. However in advanced disease the GIS (p=0.02), SBS3 (p=0.03) and HRDetect score (p=0.005) were predictive of platinum response and superior survival. PVs in g_ATM_ (n=6) or g_CHEK2_ (n=2) did not result in HRD-PDAC by any of the classifiers. In four patients, BRCA2 reversion mutations associated with platinum resistance. Conclusions: Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7-10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.

Dataset EGAD00001007571

Immune variation leads to diverse outcomes in human malaria (2020-01-15)

Falciparum malaria is clinically heterogeneous and yet in most cases the risk of life-threatening disease dramatically declines after the first few infections of life because children rapidly acquire disease tolerance (resistance to severe malaria without improved control of parasite burden). Identifying the factors that determine clinical outcome in a malaria-naive host is therefore paramount to reduce malaria mortality. However, the relative contribution of disease-causing variants of the Plasmodium var gene family versus pathogenic inflammatory cytokine cascades remains fiercely debated - we sought to reconcile these conflicting arguments by studying their interaction in vivo. To this end, two human challenge models were used to reveal the parasite-host interactions that underpin variation in falciparum malaria. To capture the diversity of human immune responses, each individual was analysed independently by tracking dynamic changes in their whole blood transcriptome through time. And to uncover evidence of preferential expansion of disease-causing variants, var gene expression was tracked in vivo from the start to end of infection. In this way, we could show that group A var genes are always expressed upon liver egress but in a minority population that does not increase over 10-days of blood cycling; there is no selection of disease-causing variants in the naive host. In fact, parasites do not respond in any way to differences or changes in host environment. On the other hand, host-intrinsic variation determines the intensity of inflammation and progression to clinical malaria. And furthermore, regulation of the interferon signaling network controls host fate. These data emphasise the role of human immune decision-making in shaping course & outcome of infection. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2020-01-15.

Dataset EGAD00001005790

GATA6 expression distinguishes classical and basal-like subtypes in advanced pancreatic cancer.

PURPOSE: To determine the impact of basal-like and classical subtypes in advanced PDAC and to explore GATA6 expression as a surrogate biomarker. EXPERIMENTAL DESIGN: Within the COMPASS trial patients proceeding to chemotherapy for advanced PDAC undergo tumour biopsy for RNA sequencing. Overall response rate (ORR) and overall survival (OS) were stratified by subtypes and according to chemotherapy received. Correlation of GATA6 with the subtypes using gene expression profiling, in situ hybridization (ISH) were explored. RESULTS: Between December 2015-May 2019, 195 patients (95%) had enough tissue for RNA sequencing; 39 (20%) were classified as basal-like and 156 (80%) as classical. RECIST response data were available for 157 patients; 29 basal-like and 128 classical where the ORR was 10% vs. 33% respectively (p=0.02). In patients with basal-like tumours treated with modified FOLFIRINOX (mFFX) (n=22) the progression rate was 60% compared to 15% in classical PDAC (p= 0.0002). Median OS in the intention to treat population (n=195) was 9.3 months for classical vs. 5.9 months for basal-like PDAC (HR 0.47 95% CI 0.32-0.69, p=0.0001). GATA6 expression by RNAseq highly correlated with the classifier (p<0.001) and ISH predicted the subtypes with sensitivity of 89% and specificity of 83%. In a multivariable analysis, GATA6 expression was prognostic (p=0.02). In exploratory analyses, basal-like tumours, could be identified by keratin 5, were more hypoxic and enriched for a T cell inflamed gene expression signature. CONCLUSIONS: The basal-like subtype is chemoresistant and can be distinguished from classical PDAC by GATA6 expression.

Dataset EGAD00001006081

Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma_CM67-RNASeq

Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (> median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

Dataset EGAD00001006306

Comprehensive cancer predisposition gene testing in an adult multiple primary tumor series shows a broad range of deleterious variants and atypical tumor phenotypes

Multiple primary tumors (MPT) affect a substantial proportion of cancer survivors and may result from various causes including inherited predisposition. Currently, germline genetic testing of MPT cases for cancer predisposition gene (CPG) variants is mostly targeted by tumor type. We ascertained pre-assessed MPT cases from genetics centers (defined as ≥2 primaries by age 60 years or ≥3 by 70) and performed whole genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment/molecular investigations, pathogenic variants in moderate and high-risk CPGs were detected in 67/440 (15.2%) of probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies including structural variants at low frequency (6/440 (1.4%) of probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant) the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice site CPG variants and at least one discordant tumor type was significantly higher than a control population (χ2=43.642 P=<0.0001). 2/67 (3%) of probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Summing together variant detection rates from a similarly ascertained previous MPT case series, the present results suggest that first-line comprehensive CPG analysis in a clinical genetics referral-based MPT cohort would detect a deleterious variant in about a third of cases.

Dataset EGAD00001004088

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 impact of a new class of transcriptional regulators, the long non-coding RNAs (lncRNAs), on the way MYC is able to influence cellular transcriptome. To this aim we have intersected RNA-sequencing data from two MYC-inducible cell lines and from a cohort of 91 mature B-cell lymphomas carrying, or not carrying, genetic variants resulting in MYC over-expression. By this approach, we identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them we focused on a lncRNA that we named MINCR, for MYC-Induced long Non-Coding RNA, showing a strong correlation with MYC expression in MYC-positive lymphomas and also in pancreatic ductal adenocarcinomas. To understand its cellular role we performed RNA interference (RNAi) experiments and found that MINCR knock-down is associated with a reduction in cellular viability, due to an impairment in cell cycle progression. Differential gene expression analysis following RNAi showed a strongly significant enrichment of cell cycle genes among the genes down-regulate following MINCR knock-down. Interestingly these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of MYC transcriptional program. Accordingly, following MINCR knock-down, we observed a reduction in the binding of MYC to the promoters of selected 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 is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Dataset EGAD00001001441

Whole-exome ultra-high throughput sequencing in brain samples of suicide victims who had suffered from major depressive disorder and control subjects who had died from other causes

We carried out whole-exome ultra-high throughput sequencing in brain samples of suicide victims who had suffered from major depressive disorder and control subjects who had died from other causes. This study aimed to reveal the selective accumulation of rare variants in the coding and the UTR sequences within the genes of suicide victims. We also analysed the potential effect of STR and CNV variations, as well as the infection of the brain with neurovirulent viruses in this behavioural disorder. As a result, we have identified several candidate genes, among others three calcium channel genes that may potentially contribute to completed suicide. We also explored the potential implication of the TGF-β signalling pathway in the pathogenesis of suicidal behaviour. To our best knowledge, this is the first study that uses whole-exome sequencing for the investigation of suicide.

Study EGAS00001002818

Adult granulosa cell tumor WGS data cohort with corresponding reference germline WGS data

Adult granulosa cell tumors (AGCTs) harbor a somatic FOXL2 c.402C>G mutation in ~95% of cases and are mainly surgically removed due to limited systemic treatment effect. In this study, potentially targetable genomic alterations in AGCTs were investigated by whole genome sequencing on 46 tumor samples and matched normal DNA. Copy number variant (CNV) analysis confirmed gain of chromosome 12 and 14, and loss of 22. Pathogenic TP53 mutations were identified in three patients with highest tumor mutational burden and mitotic activity, defining a high-grade AGCT subgroup. Within-patient tumor comparisons showed 29-80% unique somatic mutations per sample, suggesting tumor heterogeneity. FOXL2-wildtype AGCTs had DICER1, TERT(C228T) and TP53 mutations and similar CNV profiles as FOXL2-mutant tumors. Our study confirms that absence of the FOXL2 c.402C>G mutation does not exclude AGCT diagnosis. The lack of overlapping variants in targetable cancer genes, indicates the need for personalized treatment for AGCT patients

Study EGAS00001004249

Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities

Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically defined by mutations in SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be subdivided in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subtype-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT subgroups. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing representative pre-clinical models allowing for the development of subtype-specific therapies.

Study EGAS00001006422

Early on-treatment changes in circulating tumor DNA fraction and response to enzalutamide or abiraterone in metastatic castration-resistant prostate cancer

On-treatment circulating tumor DNA (ctDNA) quantification is touted as a surrogate biomarker for understanding metastatic cancer response to systemic therapy, but there is relatively little supporting data from standardized patient cohorts. Here, in a prospective observational cohort, we demonstrate that 4-week on-treatment changes in plasma ctDNA levels can accurately identify metastatic castration-resistant prostate cancer that will progress within 6 months of initiating first-line abiraterone or enzalutamide treatment. On-treatment ctDNA measurements outperformed existing clinical variables for predicting disease response and were also closely associated with differential overall survival. Although optimal timing of blood collections and thresholds for ctDNA detection will require further testing, our data highlights the potential for sequential ctDNA measurements to provide an early and reliable read-out for therapy response. The high accuracy of this information suggests that ctDNA could be useful for guiding clinical trials of therapy changes in initially non-responsive metastatic cancers.

Study EGAS00001006856

Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA

These data relate to Abbosh et al 2023 'Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA'. Some patients were included in this work that were excluded from the other TRACERx 421 papers (Frankell et al and Al-bakir et al) often due to exclusion criteria in the TRACERx clinical trial. Reasons include partial contamination with C>A artefact mutational signature (CRUK0480 ,0490) which were included for ctDNA detection analyses using high confidence clonal mutations but not clonality analyses, incomplete resection at primary surgery (CRUK0291, 0234, 0230, 0387, 0622) and a concurrent oesophageal second primary cancer (CRUK0498). The primary tissue exome sequencing for these patients that were used in this study are included in this repository. Deep targeted sequencing of (median 200) tumour specific mutations in plasma for 1069 samples from 198 TRACERx patients are also included. Additional details related to these patients and samples can be found in the supplementary or data repositories from Frankell et al and Abbosh et al.

Study EGAS00001006923

Study of complex rearrangements and mutational signatures in neuroblastoma heterogeneous risk groups.

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although certain complex genomic alterations, such as extrachromosomal DNA amplifications (ecDNA), are associated with adverse outcome and have been recurrently detected in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography of complex rearrangements along with mutational signatures derived from all variant classes, we identify previously unrecognized co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is linked to differences in the processes driving these mutational scenarios. Whereas high-risk MYCN-amplified neuroblastoma genomes were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas, on the other hand, were marked by footprints of chromosome missegregation and TOP1 mutational activity. This analysis provides a systematic perspective on the repertoire of mutational patterns that contribute to clinical neuroblastoma heterogeneity.

Study EGAS00001006983

MutWP5__CRUK_Mutographs_of_CancerBreastReduction_Mammoplasty__WG_

Cancer is a genetic disease caused by an accumulation of mutations, however many of these mutations have been identified in pathologically normal tissue. We aim to use laser-capture microscopy (LCM) to sample individual clones from breast tissue to identify whether cancer-associated mutations appear in this normal tissue, assess the mutational burden present, and identify the mutational processes causing these mutations. We will sample from a wide age range of individuals (<20 to >60 years old) to determine whether these processes differ in pre- and post-menopausal women. We will also be comparing the tissue from healthy individuals (samples from breast reduction surgery) to those at elevated risk of breast cancer (mastectomy from BRCA1/2 patients) and those who have breast cancer (adjacent normal, distal normal, and tumour tissue from mastectomy). This will allow us to determine how these processes are different between these groups of individuals, and gain insight into the earliest stages of tumor development.

Study EGAS00001003311

Methylome sequencing of cell-free DNA and RRBS of solid tissue

Early cancer detection by cell-free DNA (cfDNA) faces multiple challenges: low fraction of tumor cfDNA, molecular heterogeneity of cancer, and sample sizes not sufficient to reflect diverse patient population. We develop a cancer detection approach to address these challenges. It consists of an assay, cfMethyl-Seq, for cost-effective sequencing of cfDNA methylome (with >12-fold enrichment over whole genome bisulfite sequencing in CpG islands), and a computational method to extract methylation information and diagnose patients. Applying our system to 408 colon, liver, lung, and stomach cancer patients and controls, at 97.9% specificity we achieve 80.7% and 74.5% sensitivity in detecting all-stage and early-stage cancer, and 89.1% and 85.0% accuracy for locating tissue-of-origin of all-stage and early-stage cancer, respectively. Our approach cost-effectively retains methylome profiles of cancer abnormalities, thereby permitting the classification models to learn new features as training cohorts grow, and expanding their scope to other cancer types.

Study EGAS00001006020

Orthotopic Patient-Derived Xenografts of Pediatric Solid Tumors

Pediatric solid tumors arise from endodermal, ectodermal, or mesodermal lineages. Although the overall survival of children with solid tumors is 75%-80%, that of children with recurrent disease is below 30%. To capture the complexity and diversity of pediatric solid tumors and establish new models of recurrent disease, we developed a protocol to produce orthotopic patient-derived xenografts (O-PDXs) at diagnosis, recurrence, and autopsy. Tumor specimens were received from 168 patients, and 64 O-PDXs were established for 12 types of pediatric solid tumors. The origins of the O-PDX tumors were reflected in their gene-expression profiles and epigenomes. Genomic profiling of the tumors, including detailed clonal analysis, was performed to determine whether the clonal population in the xenograft recapitulated the patient’s tumor. We identified several drug vulnerabilities and showed that the combination of a WEE1 inhibitor (AZD1775), irinotecan, and vincristine can lead to complete response in multiple rhabdomyosarcoma O-PDX tumors in vivo.

Study EGAS00001002528

Serial TERT rearrangement breakpoint quantification in circulating tumor DNA enables minimal residual disease monitoring in patients with neuroblastoma

Approximately 30% of high-risk neuroblastomas harbor genomic TERT rearrangements at diagnosis, resulting in high TERT transcript levels and active telomerase. Combined with mutations in the RAS/MAPK/ALK signaling transduction network, these patients have a dismal prognosis. Our objective was to establish a liquid biopsy-based monitoring strategy for this particularly vulnerable patient subgroup, for whom no in-time molecular diagnostic tools exist to date. Detection and quantification of TERT rearrangement breakpoints in cell-free tumor DNA circulating in peripheral blood and bone marrow plasma from patients with TERT-driven neuroblastoma improved minimal residual disease monitoring and early relapse detection in individual patients. TERT rearrangement breakpoints, as a single marker or combined with mutations in the RAS/MAPK/ALK signaling transduction network, can be reliably analyzed from 1ng DNA using droplet digital PCR assays fit for clinical routine use. TERT rearrangement breakpoints are reflected in cell-free tumor DNA and can serve as robust biomarkers for disease activity.

Study EGAS00001007365

Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames. To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a step-wise approach to employ multiple CRISPR-Cas9 screens to elucidate functional non-canonical ORFs implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream open reading frames (uORFs) exhibited selective functionality independent of the main coding sequence. One of these, ASNSD1-uORF or ASDURF, was upregulated, associated with the MYC family oncogenes, and was required for medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future cancer genomics studies seeking to define new cancer targets.

Study EGAS00001007426

Unraveling mutagenic processes influencing the tumor mutational patterns of individuals with Constitutional Mismatch Repair Deficiency

Constitutional mismatch repair deficiency (CMMRD), caused by bi-allelic germline variants in one of the mismatch repair (MMR) genes, is a childhood cancer predisposition syndrome that often results in the development of multiple tumors early in life. A better understanding of mutational processes driving subsequent CMMRD tumors could advance optimal treatment. Therefore, we performed a genomic characterization of 41 tumors from 17 individuals with CMMRD. Mutational patterns in these tumors were found to be influenced by multiple factors, including the affected MMR genes, tumor types and somatic polymerase proofreading mutations. Temozolomide treatment left prominent mutational signatures in two second primary hematologic malignancies. Furthermore, a novel indel signature was found in 54% of the tumors, characterized by one base pair cytosine insertions in cytosine homopolymers. In conclusion, the analysis of sequential tumors in individuals with CMMRD reveals diverse mutational patterns influenced by the underlying mutated MMR gene, tumor type and treatment history.

Study EGAS00001007660

Northern Ireland COhort for the Longitudinal study of Ageing

The Northern Ireland COhort for the Longitudinal study of Ageing (NICOLA) is a representative sample of ~8,500 people from across Northern Ireland. The study, which was set up in 2012, aims to understand what it is like to grow older in Northern Ireland. • NICOLA has a strong focus on molecular biomarkers so there is complementary genetic, epigenetic and transcriptomic data available for a subset of individuals. • We inherit much of our DNA from our parents while a small amount of this material changes as we get older. In NICOLA we have 551,839 directly genotyped and 18,148,478 imputed Single Nucleotide Polymorphisms (SNPs) currently available for 2969 participants. • Summary statistics for the association of these gene polymorphisms with ~30 phenotypes were generated. • Epigenetics provides a link between our inherited DNA and environmental influences from a person’s diet, medication and lifestyle. NICOLA has the epigenetic quality controlled profiles of 1984 individuals arising from variations in DNA methylation at 862,927 genetic sites.

Study EGAS00001007915

Whole genome sequencing of patients with or at risk for HCC

Liver cancer is a major cause of cancer mortality worldwide. Screening individuals at high risk, including those with cirrhosis and viral hepatitis, provides an avenue for improved survival, but current screening methods are inadequate. In this study, we used whole-genome cell-free DNA (cfDNA) fragmentome analyses to evaluate 724 individuals from the United States, the European Union, or Hong Kong with hepatocellular carcinoma (HCC) or who were at average or high-risk for HCC. Using a machine learning model that incorporated multifeature fragmentome data, the sensitivity for detecting cancer was 88% in an average-risk population at 98% specificity and 85% among high-risk individuals at 80% specificity. We validated these results in an independent population. cfDNA fragmentation changes reflected genomic and chromatin changes in liver cancer, including from transcription factor binding sites. These findings provide a biological basis for changes in cfDNA fragmentation in patients with liver cancer and provide an accessible approach for noninvasive cancer detection.

Study EGAS00001007249

Human islet 3D chromatin maps provide insights into type 2 diabetes

The transcriptional output of each cell type is controlled by thousands of enhancers, many of which contain genetic risk variants for common diseases such as type 2 diabetes (T2D). To gain insight into how enhancer variation influences diabetes risk, we created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated islet enhancers with their target genes, often located hundreds of kilobases away. We identified hubs that show spatial and functional connections between enhancersand target genes related to islet function and diabetes. We demonstrate that genetic variants distributed across hub enhancers have a major impact on T2D heritability, and use this knowledge to identify individuals in whom islet regulatory variation has a prominent role in T2D risk. Our results demonstrate the importance of 3D chromatin architecture for molecular interpretation of T2D genetic association signals, and define genomic spaces that harbor a distinct component of the T2D polygenic burden.

Study EGAS00001002917

Engineering large chromosomal deletions by CRISPR-Cas9

Arm-level chromosomal deletions are a prevalent and defining feature of cancer. A high degree of tumor-type and subtype specific recurrencies suggest a selective oncogenic advantage. However, due to their large size it has been difficult to pinpoint the oncogenic drivers that confer this advantage. Suitable functional genomics approaches to study the oncogenic driving capacity of arm-level deletions are limited. Here we present an effective technique to engineer arm-level deletions by CRISPR-Cas9 and create isogenic cell line models. We simultaneously induce double-strand breaks (DSBs) at two ends of a chromosomal arm and select the cells that have lost the intermittent region. Using this technique, we induced arm-level deletions on chromosome 11q (65 MB) and chromosome 6q (53 MB) in neuroblastoma cell lines. Such isogenic models enable further research on the role of arm-level deletions in tumor development and growth, and their possible therapeutic potential.

Study EGAS00001005134

Robust detection of translocations in lymphoma FFPE samples using Targeted Locus Capture-based sequencing - NGS targeted capture control cohort

In routine diagnostic pathology, cancer biopsies are preserved by Formalin-Fixed, Paraffin-Embedding (FFPE) procedures for examination of (intra-) cellular morphology. Such procedures inadvertently induce DNA fragmentation, which compromises sequencing-based analyses of chromosomal rearrangements. Yet, rearrangements drive many types of hematolymphoid malignancies and solid tumors, and their manifestation is instructive for diagnosis, prognosis and treatment. Here, we present FFPE-Targeted Locus Capture (FFPE-TLC) for targeted sequencing of proximity-ligation products formed in FFPE tissue blocks, and PLIER, a computational framework that allows automated identification and characterization of rearrangements involving selected, clinically-relevant, loci. FFPE-TLC, blindly applied to 149 lymphoma and control FFPE samples, identifies the known and novel rearrangement partners. It outperforms fluorescence in situ hybridization (FISH) in sensitivity and specificity, and shows clear advantages over standard capture-NGS methods, finding rearrangements involving repetitive sequences which they typically miss. FFPE-TLC is therefore a powerful clinical diagnostics tool for accurate targeted rearrangement detection in FFPE specimens.

Study EGAS00001005325

Establishment and characterization of circulating tumor cells-derived organoids from metastatic breast cancer patients.

Over the last decade, molecular stratification and associated therapies have greatly improved the treatment options and prognosis of patients diagnosed with breast cancer. However, metastatic disease remains incurable. Breast cancer cells mainly spread to the lung, liver, bone and brain and circulating tumor cells (CTCs) that migrate in the lymphatics or blood have been implicated as the cellular source of metastasis. Detailed molecular and functional characterization of patient CTCs remains poor due to their low number in patient’s bloodstream. In addition, CTCs are difficult to culture and expand ex vivo. Here we have established a novel method for culturing and unlimitedly expanding breast cancer patient-derived CTCs in CTC-derived organoids (CDOs). Our platform can serve as a foundation for serial longitudinal isolation and expansion of breast cancer cells from peripheral blood. This may enable the identification and functional validation of therapy resistance mechanisms during disease progression, which can be used for new approaches in precision medicine.

Study EGAS00001007582

The genomic landscape of relapsed infant and childhood KMT2A-rearranged acute leukemia

To study the mechanisms of relapse in KMT2A-rearranged (KMT2A-r) acute lymphoblastic (ALL) and acute myeloid leukemia (AML), we performed whole-genome and exome sequencing of infants and children with relapsed ALL/AML (n=36), and longitudinal deep-sequencing of 257 samples in 30 patients. Somatic alterations in drug-response genes, most commonly in TP53 and IKZF1 (64%), were highly enriched in early relapse ALL (79%, 9-36 months after diagnosis), but rare in very early relapse ALL (<9 months, 9%). A marked chemotherapy-exposure signature was detected for mutations in early relapse ALL but not in very early ALL or AML relapse, in line with different mechanisms of relapse. Longitudinal analyses could track residual leukemia cells, clonal drug responses, and the upcoming relapse. These results highlight that KMT2A-r ALL and AML evade therapy differently and provide new insights into the mechanisms of relapse in this highly lethal form of pediatric acute leukemia.

Study EGAS00001008197

Human induced pluripotent stem cells display a similar mutation burden as embryonic pluripotent cells in vivo

Induced pluripotent stem cells (iPSCs) hold great promise for regenerative medicine, but genetic instability is a major concern. Embryonic pluripotent cells also accumulate mutations during early development, but how this relates to the mutation burden in iPSCs remains unknown. Here, we directly compared the mutation burden of cultured iPSCs with their isogenic embryonic cells during human embryogenesis. We generated developmental lineage trees of human fetuses by phylogenetic inference from somatic mutations in the genomes of multiple stem cells, which were derived from different germ layers. Using this approach, we characterized the mutations acquired pre-gastrulation and found a rate of 1.65 mutations per cell division. When cultured in hypoxic conditions, iPSCs generated from fetal stem cells of the assessed fetuses displayed a similar mutation rate and spectrum. Our results show that iPSCs maintain a genomic integrity during culture at a similar degree as their pluripotent counterparts do in vivo.

Study EGAS00001005939

Ductal keratin 15+ luminal progenitors in normal breast exhibit a basal-like breast cancer transcriptomic signature

Normal breast luminal epithelial progenitors have been implicated as cell of origin in basal-like breast cancer, but their anatomical localization remains understudied. Here, we combine micro-collection of uncultured organoids and single cell mRNA sequencing (scRNA-seq) of FACS-sorted luminal epithelial cells with multicolor imaging to profile ducts and terminal ductal lobular units (TDLUs) and compare them with breast cancer subtypes. Unsupervised clustering reveals eleven distinct clusters and a differentiation trajectory starting with keratin 15+ (K15+) progenitors enriched in ducts. Spatial mapping of luminal progenitors is confirmed at the protein level by staining with critical duct markers. Comparison of the gene expression profiles of normal luminal cells with those of breast cancer subtypes reveals a strong correlation between normal breast K15+ ductal progenitors and basal-like breast cancer. We propose that K15+ basal-like breast cancers originate in ductal progenitors, which emphasizes the importance of not only lineages but also cellular position within the ductal-lobular tree.

Study EGAS00001005963

Immune-awakening revealed by peripheral T cell dynamics after one cycle of immunotherapy

BACKGROUND: Our understanding of T cell evolution under checkpoint inhibitors (CPI) is still incomplete, restraining the achievement of full benefit from CPI.OBJECTIVES: We studied peripheral T cell turnover and evolution and their prognostic value after 3weeks on treatment with CPI (one cycle) analysing T cell receptor-β (TCR) sequences in plasma cell-free DNA (cfDNA) and PBMC, and performing a phenotypic analysis of peripheral T cell subsets.RESULTS: Peripheral T cell turnover and TCR repertoire dynamics correlated with response. Additionally, the cfDNA TCR repertoire reorganisation fingerprint correlated with the expansion of an immune-effector subset of peripheral T cells that predicted treatment response and identified the patients with longer overall survival.CONCLUSIONS: Here we show that prognostic changes in peripheral T cells occur within 3 weeks of commencing CPI treatment. This dynamic immune-awakening informs on the immune-system reorganization can be monitored using minimally invasive liquid biopsies.

Study EGAS00001004043

MGA-NUTM1 fusion in high grade spindle cell sarcoma

NUT carcinoma (NC) with mesenchymal differentiation have rarely been described in the literature. In this report, we describe the characterization of two cases of high grade spindle cell sarcoma harboring a novel MGA-NUTM1 fusion. Whole genome sequencing identified the presence of complex rearrangements resulting in a MGA-NUTM1 fusion gene in the absence of other significant somatic mutations. The fusion protein was predicted to retain nearly the entire protein sequence of both MGA (exons 1-22) and NUTM1 (exons 3-8). In contrast to typical cases of NC, these cases were successfully treated with aggressive local control measures (surgery and radiation) and both patients remain alive without disease. These cases describe a new subtype of NUTM1-rearranged tumors warranting expansion of diagnostic testing to evaluate for the presence of MGA-NUTM1 or alternative NUTM1 gene fusions in the diagnostic workup of high grade spindle cell sarcomas or small round blue cell tumors of ambiguous lineage.

Study EGAS00001003341

GENOMIC MUTATION LANDSCAPE OF SKIN CANCERS FROM DNA REPAIR-DEFICIENT XERODERMA PIGMENTOSUM PATIENTS

Xeroderma pigmentosum (XP) is a genetic disorder caused by mutations in genes of the Nucleotide Excision Repair (NER) pathway (groups A-G) or in Translesion Synthesis (TLS) DNA polymerase η (group V). XP is associated with an increased skin cancer risk, reaching, for some groups, several thousand-fold compared to the general population. Here, we analyzed 38 skin cancer genomes from five XP groups. We found that the activity of NER determines heterogeneity of the mutation rates across skin cancer genomes and that transcription-coupled NER extends beyond the gene boundaries reducing the intergenic mutation rate. Mutational profile in XP-V tumors revealed the role of polymerase η in the error-free bypass of (i) rare TpG and TpA DNA lesions, (ii) 3’ nucleotides in pyrimidine dimers, and (iii) TpT photodimers. Our study unravels the genetic basis of skin cancer risk in XP and provides insights into the mechanisms reducing UV-induced mutagenesis in the general population.

Study EGAS00001006732

HCA_Gonads_Foetal_EU_H2020_HUGODECA_RNA

Cell Atlas of human fetal gonads The HUGODECA consortium brings together leading European academic and industrial experts (from 5 different EU countries). It includes active contributors and executives of the Human Cell Atlas (HCA) which will ensure complementarity with other ongoing HCA efforts. The overall HUGODECA concept is grounded on the integration of multiple synergistic expertise and technologies: Single cell profiling, spatial transcriptomic, 2D Mass cytometry and cyclic immunofluorescence and 3D imaging of optically cleared gonads. HUGODECA will implement novel tools, analytical and computational methods to process and integrate multidimensional OMICS and image data across different platforms. It will evaluate the accuracy of ex vivo culture models of human gonadal development and assess consequences of altering key signalling pathways. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001004723

A biobank of patient-derived pediatric brain tumor models

Brain tumors are the leading cause of cancer-related death in children. Genomic studies have provided insights into molecular subgroups and oncogenic drivers of pediatric brain tumors that may lead to novel therapeutic strategies. To evaluate new treatments, better preclinical models adequately reflecting the biological heterogeneity are needed. Through the Children’s Oncology Group ACNS02B3 study, we have generated and comprehensively characterized 30 patient-derived orthotopic xenograft models and 7 cell lines representing 14 molecular subgroups of pediatric brain tumors. Patient-derived orthotopic xenograft models were found to be representative of the human tumors they were derived from in terms of histology, immunohistochemistry, gene expression, DNA methylation, copy number, and mutational profiles. In vivo drug sensitivity of targeted therapeutics was associated with distinct molecular tumor subgroups and specific genetic alterations. These models and their molecular characterization provide an unprecedented resource for the cancer community to study key oncogenic drivers and to evaluate novel treatment strategies.

Study EGAS00001002536

Epigenomics and Single-cell Sequencing Define a Developmental Hierarchy in Langerhans Cell Histiocytosis

Langerhans cell histiocytosis (LCH) is a rare neoplasm predominantly affecting children. It occupies a hybrid position between cancers and inflammatory diseases, and it provides an attractive model for studying cancer development. To explore the molecular mechanisms underlying the pathophysiology of LCH and its characteristic clinical heterogeneity, we investigated the transcriptomic and epigenomic diversity in primary LCH lesions. Using single-cell RNA sequencing, we identified multiple recurrent types of LCH cells within these biopsies, including putative LCH progenitor cells and several subsets of differentiated LCH cells. We confirmed the presence of proliferative LCH cells in all analysed biopsies using immunohistochemistry, and we defined an epigenomic and gene-regulatory basis of the different LCH cell subsets by chromatin accessibility profiling. In summary, our single-cell analysis of LCH uncovered an unexpected degree of cellular, transcriptomic, and epigenomic heterogeneity among LCH cells, indicative of complex developmental hierarchies in LCH lesions.

Study EGAS00001003822

A cycling, progenitor-like cell population at the root of atypical teratoid rhabdoid tumor subtype differentiation trajectories.

Background: Atypical teratoid rhabdoid tumors (ATRTs), characterized by the loss of SMARCB1, are among the most lethal pediatric central nervous system (CNS) tumors. Three molecular subtypes have been identified, each defined by distinct molecular and clinical features. No subtype-specific treatments are available, highlighting the necessity to better understand inter- and intra-subtype heterogeneity. Methods: We generated a single-nucleus transcriptome atlas of ATRTs, validated by single-cell ATAC-seq and spatial transcriptomics, to study subtype-specific differentiation trajectories. Results: We identified brain progenitor-like expression profiles within unique subtype-specific differentiation lineages. A shared cycling, progenitor-like cell population, interspersed throughout tumors, was observed across all ATRT samples. Moreover, we demonstrate that subtype-specific differentiation trajectories can be induced pharmacologically in ATRT tumoroids directing tumor cells towards a non-proliferative, mature phenotype. Conclusions: Our study reveals that ATRT subtypes exhibit signaling trajectories mirroring normal fetal brain development, enabling the development of maturation therapies tailored towards ATRTs.

Study EGAS00001008108

Whole-Genome Sequencing of a Healthy Aging Cohort.

Studies of long-lived individuals have revealed few genetic mechanisms for protection against age-associated disease. Therefore, we pursued genome sequencing of a related phenotype-healthy aging-to understand the genetics of disease-free aging without medical intervention. In contrast with studies of exceptional longevity, usually focused on centenarians, healthy aging is not associated with known longevity variants, but is associated with reduced genetic susceptibility to Alzheimer and coronary artery disease. Additionally, healthy aging is not associated with a decreased rate of rare pathogenic variants, potentially indicating the presence of disease-resistance factors. In keeping with this possibility, we identify suggestive common and rare variant genetic associations implying that protection against cognitive decline is a genetic component of healthy aging. These findings, based on a relatively small cohort, require independent replication. Overall, our results suggest healthy aging is an overlapping but distinct phenotype from exceptional longevity that may be enriched with disease-protective genetic factors.

Study EGAS00001002306

White blood cell and cell-free DNA analyses for detection of residual disease in gastric cancer

Liquid biopsies are providing new opportunities for detection of residual disease in cell-free DNA (cfDNA) after surgery but may be confounded through identification of alterations arising from clonal hematopoiesis. Here, we identify circulating tumor-derived DNA alterations (ctDNA) through ultrasensitive targeted sequencing analyses of matched cfDNA and white blood cells from the same patient. We apply this approach to analyze samples from patients in the CRITICS trial, a phase III randomized controlled study of perioperative chemotherapy in patients with operable gastric cancer. After filtering alterations derived from matched white-blood cells, the presence of ctDNA predicts recurrence when analyzed within nine weeks after preoperative treatment and after surgery in patients eligible for multimodal treatment. These analyses provide a facile method for distinguishing ctDNA from other cfDNA alterations and highlight the utility of ctDNA as a predictive biomarker of patient outcome to perioperative cancer therapy and surgical resection in patients with gastric cancer.

Study EGAS00001004114

22929 results for "ega"

in 73.88 milliseconds.

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

Inferring expressed genes by whole-genome sequencing of plasma DNA

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

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

Chordoma_Sequencing_Project_Whole_Genome

APCDR Uganda GWAS: Genome-wide sequence variation and susceptibility loci for cardiometabolic traits in a sub-Saharan African population (UG2G component)

CD74-NRG1 fusions in lung adenocarcinoma

Biallelic DICER1 mutations in sporadic pleuropulmonary blastoma

Tracking the origins and drivers of subclonal metastatic expansion in prostate cancer

Origins and functional consequence of somatic mitochondrial DNA mutations

Frequent alterations in cytoskeleton remodeling genes in primary and metastatic Chinese lung adenocarcinomas

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

Trans-ethnic genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation

Searching for variants associated with endometriosis

A Genomic History of Aboriginal Australia

Rare disruptive mutations in ciliary function genes contribute to testicular cancer susceptibility

EATL-II STUDY

Recurrent somatic mutations in POLR2A define a distinct subset of meningiomas

Divergence between high metastatic tumor burden and low circulating tumor DNA concentration in metastasized breast cancer

Transcriptional and functional profiling defines human small intestinal macrophage subsets

Myelodysplastic cells in patients re-program mesenchymal stromal cells to establish a transplantable stem cell-niche disease unit.

Resequencing_candidate_genes_for_male_spermatogenic_impairment

Whole exome sequencing for gallbladder cancer in Xinhua Hospital Affiliated to Shanghai Jiao Tong University, School of Medicine

Genetic Basis of Hepatosplenic T Cell Lymphoma (HSTL)

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

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

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

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

Comparing nodal versus bony metastatic spread using tumour phylogenies

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

Genetic landscape of Early T-cell precursor acute lymphoblastic leukaemia

Epigenomic alterations define lethal CIMP-positive ependymomas of infancy

Reference epigenomes generated as part of the International Human Epigenomics Consortium (IHEC)

Genomic catastrophes frequently arise in esophageal adenocarcinoma and drive tumorigenesis.

Tracking the genomic evolution of esophageal adenocarcinoma through neoadjuvant chemotherapy

RNA_seq_of_Toxoplasma_gondii_response_in_human_macrophages

Novel mutational mechanisms and drivers in Pancreatic Neuroendocrine Tumours

Exome Sequencing of Gastric Cancer

Complete Genomics paired end sequencing; Ovarian cancer

Exome_sequencing_of_blastic_plasmacytoid_dendritic_cell_neoplasms

Association of Age at Diagnosis and Genetic Mutations in Patients with Neuroblastoma

Exome-sequencing identifies new oncogenes and tumor suppressor genes recurrently altered in hepatocellular carcinoma

Genome Landscape of High-Grade Serous Ovarian Cancer

A Genomics-Based Classification of Human Lung Tumors

High number of somatic mutations found in the healthy blood compartment of a 115-year-old woman reveals oligoclonal hematopoiesis

Germline and somatic SMARCA4 mutations characterize small-cell carcinoma of the ovary, hypercalcemic type.

Integrated genetic and epigenetic analysis defines novel molecular clusters in rhabdomyosarcoma

65 prostate cancer cases WGS and transcriptome sequencing project

SeqControl: Process Control for DNA Sequencing

Genetic_profiling_of_mucosal_melanoma

Exome sequencing of hepatocellular carcinomas identifies new mutational signatures and potential therapeutic targets

Whole-genome plasma sequencing reveals focal amplifications as a driving force in metastatic prostate cancer

Neuroblastoma relapse trio series from the AMC

Recurrent mTORC1-activating RRAGC mutations in follicular lymphoma

Recurrent adeno-associated virus 2-related insertional mutagenesis in human hepatocellular carcinomas

Barcoding reveals complex clonal dynamics of de novo transformed human mammary cells

A comprehensive assessment of somatic mutation detection in cancer using whole genome sequencing

Integrated genomic characterization of IDH1 mutant Glioma malignant progression

Spatial and Temporal Homogeneity of Driver Mutations in Diffuse Intrinsic Pointine Glioma

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

Multi-layered population structure in Island Southeast Asians

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

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

whole-genome sequencing in 17 ESCC cases and whole-exome sequencing in 71 cases

Subtype specific progression from DCIS to invasive breast cancer

Variation in the Glucose Transporter gene SLC2A2 is associated with glycaemic response to metformin

Hi-C dataset for testicular germ cell tumour GWAS risk loci, as described in the Oncoarray Litchfield et al. 2016 paper.

New Brain Tumor Entities Emerge from Molecular Classification of CNS-PNETs

A living biobank of breast cancer organoids captures disease heterogeneity

The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

Fetal origins of malignant germ cell tumours

Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - TGS

Distribution of ctDNA levels in plasma of early-stage non-small cell lung cancer patients measured using personalised ctDNA analysis

pancreatic ductal adenocarcinoma exomes in study: Pro-immunogenic Impact of MEK inhibition combined with an anti-CD40 immunostimulatory antibody

Preferential infiltration of unique Vγ9Jγ2‐Vδ2 T cells into glioblastoma multiforme

Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - WGS

The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

Transcriptional profiling of tauopathies in human IPS-derived neurons (2019-08-21)

MutWP5__CRUK_Mutographs_of_CancerBreastReduction_Mammoplasty__WG_