22928 results for "ega"

in 70.55 milliseconds.

• RNA-seq of high grade serous ovarian tumours

  Somatic RNA for 37 samples was extracted using the Qiagen Qiasymphony RNA protcol (cat no 931636). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including DNase digestion). The resulting RNA the underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as measures of RNA purity. A260/280 had to be 2.0 and A260/230 had to be 2.0-2.2. Then RNA was quantified using LifeTechnologies Qubit RNA BR kit (cat no Q10210). RNAseq was carried out by the Edinburgh Clinical Research Facility on an Illumina NExtSeq500. Total RNA samples were assessed on the Agilent Bioanalyser (Agilent Technologies, #G2939AA) with the RNA 6000 Nano Kit (#5067-1512) for quality and integrity of total RNA, and then quantified using the Qubit 2.0 Fluorometer (Thermo Fisher Scientific Inc, #Q32866) and the Qubit RNA HS assay kit (#Q32855). Libraries were prepared from total-RNA sample using the NEBNext Ultra 2 Directional RNA library prep kit for Illumina (#E7760S) with the NEBNext rRNA Depletion kit (#E6310) according to the provided protocol. 400ng of totalRNA was then added to the ribosomal RNA (rRNA) depletion reaction using the NEBNext rRNA depletion kit (Human/mouse/rat) (#E6310). This step uses specific probes that bind to the rRNA in order to cleave it. rRNA-depleted RNA was then DNase treated and purified using Agencourt RNAClean XP beads (Beckman Coulter Inc, #66514). RNA was then fragmented using random primers before undergoing first strand and second strand synthesis to create cDNA. cDNA was end repaired before ligation of sequencing adapters, and libraries were enriched by PCR using the NEBNext Multiplex oligos for Illumina set 1 and 2 (#E7500). Final libraries had an average peak size of 271bp. Libraries were quantified by fluorometry using the Qubit dsDNA HS assay and assessed for quality and fragment size using the Agilent Bioanalyser with the DNA HS Kit (#5067-4626). Sequencing was performed using the NextSeq 500/550 High-Output v2 (150 cycle) Kit (# FC- 404-2002) on the NextSeq 550 platform (Illumina Inc, #SY-415-1002). Libraries were combined in an equimolar pool based on the library quantification results and run across 5 High-Output Flow Cell v2.5.

  Dataset EGAD00001010139

• Medulloblastoma cfDNA low-coverage WGS

  Whole genome sequencing of cell free DNA from CSF across timepoints from medulloblastoma clinical trial patients.

  Dataset EGAD00001008277

• panALL exome data set3, 650 samples

  panALL exome data set3, 650 samples

  Dataset EGAD00001008290

• RNA-Seq of 584 DLBCL Cases of Non-China Cohort from Phoenix Clinical Trial

  Standard RNA-Seq datasets. Check the associated paper for more details.

  Dataset EGAD00001008131

• Thymic epithelial transplantation for complete DiGeorge syndrome: RNA (2025-07-22)

  Single cell transcriptomics study of thymic transplant biopsies Allogeneic thymus transplantation is the only curative therapy for complete DiGeorge Syndrome (cDGS), a rare severe primary immunodeficiency characterised by athymia. GOSH is one of only two centres worldwide to offer this treatment. Despite a lack of major histocompatibility complex (MHC)-matching between donor and host, transplanted thymus becomes repopulated by recipient bone marrow derived precursor cells and supports development of functional T-cells. The mechanisms underlying thymopoiesis in this context are poorly understood, but over time we observe reconstitution of T-cell immunity, with the ability to produce host naïve T-cells showing a broad T-cell receptor (TCR) repertoire and to generate MHC-restricted T-cell proliferative responses. Although lifesaving, the achieved immunological reconstitution is typically not complete with circulating T-cell numbers usually remaining below the age related normal ranges. Additionally, we observe persistence of donor-derived T-cells of unknown clinical significance. To gain more insight into the mechanisms by which MHC-mismatched transplanted thymus supports T-cell development with self-tolerance, as well as into the basis of suboptimal T-cell immunity, we now aim to investigate immune reconstitution after thymus transplantation in further detail by using single-cell transcriptomics, applied to thymic transplant biopsies and peripheral blood samples collected during standard post-transplant patient care. By identifying which lineages of host- derived cells repopulate the thymic tissue after transplantation, we will be able to address the role of MHC in positive and negative T-cell selection during T-cell differentiation. We will also be able to clarify the exact ontogeny of the persistent donor T-cells, as well as their possible role. Understanding the mechanisms of action of HLA-mismatched transplanted thymus will contribute to treatment optimisation. Additionally, our research provides a unique opportunity to further investigate key immunological concepts, such as tolerance and autoimmunity, challenging existing paradigms in thymus immunology. 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 2025-07-22.

  Dataset EGAD00001015648

• Germline variants in patients with rare cancers - control samples WGS and WES

  The control samples (mostly blood) of 351 samples (paired WGS and WES sequencing) are in this dataset. The WGS was in nearly all cases at an Illumina HighSeq X Ten with the Illumina TruSeq Nano DNA Kit. The WES mostly on Illumina HighSeq 4000 with the Agilent SureSelect V5 plus UTRs Kit.

  Dataset EGAD00001010047

• Gut microbiota of the Pinggu cohort

  Shotgun metagenomic sequencing data of a total 2,338 fecal DNA samples from adults of the Pinggu cohort.

  Dataset EGAD00001007064

• Tumor transcriptome profiling from the OAK and POPLAR randomized clinical trials of atezolizumab in 2L NSCLC

  This dataset includes full tumor transcriptomes from 891 advanced NSCLC tumors. These data originate from pre-treatment samples from two large randomized clinical trials for second-line non-small cell lung cancer (POPLAR and OAK). The patients in these trials were treated with either the PD-L1 inhibitor atezolizumab or chemotherapy.

  Dataset EGAD00001007703

• H3Africa - Identification and characterization of novel hereditary neurological disease genes in Mali.

  Despite the vast diversity of its populations, genetic studies in Africa have been limited. African populations, Malians in particular, have a high rate of intra-ethnic and consanguineous marriage, resulting in increased prevalence of autosomal recessive diseases. Family-based genetic studies can be limited in developed countries due to small sibships. The average fertility rate in Mali is over 6 births per woman, offering a unique opportunity to find new disease genes or mutations that can then be studied in other populations. Hereditary neurological diseases are very debilitating diseases, and developing countries, particularly in Africa, pay a high price in terms of disability-adjusted life years. Although most are currently untreatable, increasing awareness about hereditary neurological disorders can reduce this burden. This study will help to identify and characterize novel hereditary neurological disease genes in Mali. These genes are likely to be important in the normal function of the nervous system and to have important pathophysiological implications for African and other populations. This study will also train physicians and students in the characterization of neurodegenerative diseases as well as in genetic technology and molecular biology, and build a suitable research environment that will retain them.

  Study EGAS00001003016

• Genome‐wide postnatal changes in immunity following fetal inflammatory response

  The fetal inflammatory response (FIR) increases the risk of perinatal brain injury, particularly in extremely low gestational age newborns (ELGANs, < 28 weeks of gestation). One of the mechanisms contributing to such a risk is a postnatal intermittent or sustained systemic inflammation (ISSI) following FIR. The link between prenatal and postnatal systemic inflammation is supported by the presence of well‐established inflammatory biomarkers in the umbilical cord and peripheral blood. However, the extent of molecular changes contributing to this association is unknown. Using RNA sequencing and mass spectrometry proteomics, we profiled the transcriptome and proteome of archived neonatal dried blood spot (DBS) specimens from 21 ELGANs. Comparing FIR‐affected and unaffected ELGANs, we identified 782 gene and 27 protein expression changes of 50% magnitude or more, and an experiment‐wide significance level below 5% false discovery rate. These expression changes confirm the robust postnatal activation of the innate immune system in FIR‐affected ELGANs and reveal for the first time an impairment of their adaptive immunity. In turn, the altered pathways provide clues about the molecular mechanisms triggering ISSI after FIR, and the onset of perinatal brain injury.

  Study EGAS00001003635

• Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

  Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes depending on age and other factors. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome sequencing data from 702 neuroblastoma samples, half of which are new data. Over 95% of samples harbor at least one recurrent driver alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years, suggesting age-specific susceptibility to genetic alterations. Mutational signature analysis shows that reactive oxygen-species (ROS)-induced mutations are the most common cause of single nucleotide variant (SNV) drivers in neuroblastoma, including most ALK and Ras pathway SNVs. ROS-induced mutagenesis is both an early event, as it causes truncal variants, and continuous as it also induces relapse-specific mutations. ROS-induced mutations are more abundant in neuroblastomas with MYCN amplification, 17q gains, and increased expression of mitochondrial ribosome genes. We observe recurrent FGFR1 N546K or internal tandem duplication variants in 6 patients and ALK N-terminal structural alterations in 5 samples, identifying additional patients amenable to precision therapy.

  Study EGAS00001003931

• SINGLE-CELL RNA SEQUENCING Single-cell RNA sequencing was performed on 13 ‘mild-moderate’ and 10 ‘critical’ COVID19 PBMC samples

  Summary Epidemiological and clinical reports have indicated that SARS-CoV-2 virulence hinges upon the triggering of an aberrant host immune response, more so than on direct virus-induced cellular damage. To elucidate the immunopathology underlying COVID-19 severity, we performed cytokine and multiplex immune profiling in mild-moderate and critically ill COVID-19 patients. Hypercytokinemia in COVID-19 differed from the IFN-γ-driven cytokine storm in macrophage activation syndrome, and was more pronounced in critical versus mild-moderate COVID-19. Systems modelling of cytokine levels paired with deep-immune profiling showed that classical monocytes drive this hyper-inflammatory phenotype and that a reduction in T-lymphocytes correlates with disease severity, with CD8+ cells being disproportionately affected. Expression of antigen presenting machinery was also reduced in critical disease. Furthermore, we found that neutrophils contributed to disease severity and local tissue damage by amplification of hypercytokinemia and the formation of neutrophil extracellular traps. Together these findings suggest a myeloid-driven immunopathology, in which hyperactivated neutrophils and an ineffective adaptive immune system act as mediators of COVID-19 disease severity.

  Study EGAS00001005039

• An alternative splicing modulator decreases mutant HTT and improves the molecular fingerprint in Huntington’s disease patient neurons

  Huntington's disease (HD) is a neurodegenerative disorder caused by poly-Q expansion in the Huntingtin (HTT) protein. Here, we delineate elevated mutant HTT (mHTT) levels in patient-derived cells including fibroblasts and iPSC derived cortical neurons using a GLP approved HTT assay. HD patients’ fibroblasts and cortical neurons recapitulate aberrant alternative splicing as a molecular fingerprint of HD. Branaplam is a splicing modulator currently tested in a phase II study in HD (NCT05111249). The drug lowers total HTT (tHTT) and mHTT levels in fibroblasts, iPSC, cortical progenitors, and neurons in a dose dependent manner at an IC50 consistently below 10nm without inducing cellular toxicity. Branaplam promotes inclusion of non-annotated novel exons. Amongst, a 115bp frameshift-inducing exon in the HTT transcript in Branaplam treated cells from Ctrl and HD patients leading to a profound reduction of HTT RNA and protein levels. Importantly, Branaplam ameliorates aberrant alternative splicing in HD patients’ fibroblasts and cortical neurons. These findings highlight the applicability of splicing modulators in the treatment of CAG repeat disorders and decipher their molecular effects associated with the pharmacokinetic and -dynamic properties in patient-derived cellular models.

  Study EGAS00001006289

• H3Africa - Kidney Disease Research Network

  Chronic kidney disease (CKD) affects an estimated 14% of adults in sub-Saharan Africa but very little research has been done on the etiology, progression and prevention of CKD on the continent. As part of the Human Heredity and Health in Africa (H3A) Consortium, the H3Africa Kidney Disease Research Network was established to study prevalent forms of kidney disease in sub-Saharan Africa and increase the capacity for genetics and genomics research. The study is performing comprehensive phenotypic characterization with analysis of environmental and genetic factors from seven clinical centers in two African countries (Ghana and Nigeria) over a five-year period. Approximately 4,000 participants with specified kidney disease diagnoses and 4000 control participants will be enrolled in the four African countries. In addition, approximately 50 families with hereditary glomerular disease will be enrolled. The study includes both pediatric and adult participants ages <1 to 74yrs across a broad spectrum of kidney diseases secondary to hypertension, diabetes, human immunodeficiency virus, sickle cell disease, biopsy-proven glomerular disease and CKD of unknown origin. Clinical and demographic data with biospecimens are collected for assessment of clinical, biochemical and genetic markers of kidney disease.

  Study EGAS00001006558

• WTSI CGP Data access committee

  Dac EGAC00001000000

• NPC genome-wide human SNP array data

  nasopharyngeal carcinoma genome-wide human SNP array data for 4083 NPC cases and 4811 controls

  Dataset EGAD00010002296

• Efficacy of two different FGFR-inhibitors in a patient with extrahepatic cholangiocarcinoma harboring an FGFR2 mutation

  Cholangiocarcinomas (CCAs) is a type of cancer with few effective systemic therapies. Elucidation of the molecular landscape of the disease from genomic studies based on next generation sequencing (NGS) has contributed to the introduction of new targeted therapies. One of these treatments consists of a class of small molecules that target members of the FGFR family of receptor tyrosine kinases. These drugs are effective and have been approved for cholangiocarcinomas with fusions or rearrangements of FGFR genes. In contrast, the role of these inhibitors in cholangiocarcinomas with mutations in FGFR genes is less well defined. We report here a patient with a cholangiocarcinoma bearing a FGFR2 p.Ser252Trp mutation. The patient was treated with two different FGFR inhibitors, as the first caused ocular toxicity. She obtained clinical benefit from both. This case illustrates the efficacy of FGFR inhibitors on cholangiocarcinoma with specific point mutations. This is the first case to report the clinical benefit of these drugs in FGFR2 p.Ser252Trp mutation. Clinical benefit can be sustained, as seen in our patient. Our case also shows that FGFR inhibitors-induced adverse effects, such as ocular toxicities, may not recur after re-challenge with an alternative drug of the same class.

  Study EGAS50000000015

• snRNA-seq analysis in multinucleated myogenic FSHD cells identifies heterogeneous FSHD transcriptome signatures associated with embryonic-like program activation and oxidative stress-induced apoptosis

  The sporadic nature of DUX4 expression in FSHD muscle challenges comparative transcriptome analyses between FSHD and control samples. A variety of DUX4 and FSHD-associated transcriptional changes have been identified, but bulk RNA-seq strategies prohibit comprehensive analysis of their spatiotemporal relation, interdependence and role in the disease process. In this study, we used single-nucleus RNA-sequencing of nuclei isolated from patient- and control-derived multinucleated primary myotubes to investigate the cellular heterogeneity in FSHD. Taking advantage of the increased resolution in snRNA-sequencing of fully differentiated myotubes, two distinct populations of DUX4-affected nuclei could be defined by their transcriptional profiles. Our data provides insights into the differences between these two populations and suggests heterogeneity in two well-known FSHD-associated transcriptional aberrations: increased oxidative stress and inhibition of myogenic differentiation. Additionally, we provide evidence that DUX4-affected nuclei share transcriptome features with early embryonic cells beyond the well-described cleavage stage, progressing into the 8-cell and blastocyst stages. Altogether, our data suggests that the FSHD transcriptional profile is defined by a mixture of individual and sometimes mutually exclusive DUX4-induced responses and cellular state-dependent downstream effects.

  Study EGAS00001007635

• Single-cell profiling reveals mechanisms of response to anti-PD-L1 versus anti-PD-L1 combined with anti-CTLA4 in head and neck squamous cell carcinoma

  Cancer patients often receive a combination of PD-(L)1 and CTLA4 inhibitors, but the mechanisms underlying their concerted actions remain unclear. We conducted a neoadjuvant study in head and neck squamous cell carcinoma (HNSCC) involving anti-PD-L1 monotherapy versus anti-PD-L1/anti-CTLA4 combination. Single-cell profiling of on- versus pre-treatment biopsies confirmed T-cell expansion as response biomarker. We identified a type 1 immune response accompanying T-cell expansion on-treatment across treatment arms, and herein characterized co-localized IgG plasma cell expansion as important contributor. In pre-treatment biopsies, similar features correlated with expansion upon anti-PD-L1, but not anti-PD-L1/anti-CTLA4. By profiling tumor-draining lymph nodes, we found the addition of anti-CTLA4 to trigger activation and subsequent trafficking of CD4+ T-cells via the blood, to become T-helper 1 cells in the tumor. anti-PD-L1/anti-CTLA4 indeed facilitated co-localized expansion of CD4+ and CD8+ T-cells in tumors versus mostly CD8+ T-cells with anti-PD-L1 alone. Hence, we identify complimentary mechanisms underlying anti-PD-L1/anti-CTLA4 therapy in HNSCC.

  Study EGAS50000000037

• Single cell RNAseq of stenotic, inflamed and non-inflamed transmural lesions from patients with Crohn's disease

  This dataset contains a gene-cell matrix derived from single-cell RNA sequencing (scRNA-seq) data of ileal tissue from Crohn's disease (CD) patients and colorectal cancer (CRC) patients. It includes: Crohn's Disease Patients: A trio of transmural lesions (stenotic, inflamed, and non-inflamed) from each patient. Colorectal Cancer Patients: Unaffected ileal tissue used as external non-inflamed control. Cell Level Metadata: The dataset includes relevant cell-level metadata such as cell type annotations used in the study. Experimental Details: Platform: 10x Genomics Chromium Single Cell 3' GEX Sequencing: Illumina NovaSeq Processing: Data processed with Cell Ranger software. Resulting count matrices were merged for downstream analysis, including integration and dimensionality reduction. Dataset Composition: Crohn's Disease Patients: 10 patients with 3 samples each (non-inflamed, inflamed, stenotic), totaling 30 samples. Colorectal Cancer Patients: 5 patients with 1 sample each of unaffected tissue, totaling 5 samples. Data Provided: Merged Raw Count Matrix: The final merged raw count matrix used for downstream analysis. Cell Metadata File: Contains details of sample, tissue, and patient for each cell in the count matrix. Barcodes File: Indicate each cell barcode which also encodes the sample, tissue, and patient details for each cell. CD.S_Inf: Stenotic Corhn's disease inflamed samples CD.S_Sten: Stenotic CD patient stenosis sample CD.S_Prox: Stenotic CD Patient - proximal non-inflamed sample CC.C_Prox: CRC Patient proximal unaffected sample eg: A barcode 'CC.C_1_Prox_AAGTCGTAGACCCTTA' indicates CRC Patient unaffected proximal sampe from CRC Patient no.1 and the nucleic acid sequence indicate a unique cell from this sample. Total Samples: Crohn's Disease (CD) Patients: 30 samples Colorectal Cancer (CRC) Patients: 5 samples Patient_no Sample Sample_type 1 CC.C_1 CC.C_1_Prox CC.C_Prox 2 CD.S_1 CD.S_1_Prox CD.S_Prox 3 CD.S_1 CD.S_1_Infl CD.S_Infl 4 CD.S_1 CD.S_1_Sten CD.S_Sten 5 CC.C_2 CC.C_2_Prox CC.C_Prox 6 CD.S_2 CD.S_2_Prox CD.S_Prox 7 CD.S_2 CD.S_2_Infl CD.S_Infl 8 CD.S_2 CD.S_2_Sten CD.S_Sten 9 CC.C_3 CC.C_3_Prox CC.C_Prox 10 CC.C_4 CC.C_4_Prox CC.C_Prox 11 CD.S_3 CD.S_3_Prox CD.S_Prox 12 CD.S_3 CD.S_3_Infl CD.S_Infl 13 CD.S_3 CD.S_3_Sten CD.S_Sten 14 CD.S_4 CD.S_4_Prox CD.S_Prox 15 CD.S_4 CD.S_4_Infl CD.S_Infl 16 CD.S_4 CD.S_4_Sten CD.S_Sten 17 CC.C_5 CC.C_5_Prox CC.C_Prox 18 CD.S_5 CD.S_5_Prox CD.S_Prox 19 CD.S_5 CD.S_5_Infl CD.S_Infl 20 CD.S_5 CD.S_5_Sten CD.S_Sten 21 CD.S_6 CD.S_6_Prox CD.S_Prox 22 CD.S_6 CD.S_6_Infl CD.S_Infl 23 CD.S_6 CD.S_6_Sten CD.S_Sten 24 CD.S_7 CD.S_7_Prox CD.S_Prox 25 CD.S_7 CD.S_7_Infl CD.S_Infl 26 CD.S_7 CD.S_7_Sten CD.S_Sten 27 CD.S_8 CD.S_8_Prox CD.S_Prox 28 CD.S_8 CD.S_8_Infl CD.S_Infl 29 CD.S_8 CD.S_8_Sten CD.S_Sten 30 CD.S_9 CD.S_9_Prox CD.S_Prox 31 CD.S_9 CD.S_9_Infl CD.S_Infl 32 CD.S_9 CD.S_9_Sten CD.S_Sten 33 CD.S_10 CD.S_10_Prox CD.S_Prox 34 CD.S_10 CD.S_10_Infl CD.S_Infl 35 CD.S_10 CD.S_10_Sten CD.S_Sten

  Dataset EGAD50000000559

• Dataset developed for use with EOSC4Cancer of synthetic colorectal cancer tumor/normal pairs.

  This dataset contains 10 tumor and normal pairs synthetic WGS data of colorectal cancer that were simulated in a standard format of Illumina paired-end reads. The NEAT read simulator (version 3.0, https://github.com/zstephens/neat-genreads) was utilized to synthetize these 10 pairs of tumor and normal WGS data. In the procedure of data generation, simulated parameters (i.e., sequencing error statistics, read fragment length distribution and GC% coverage bias) were learned from data models provided by NEAT. The average sequencing depth for tumor and normal samples aimed to reach around 110X and 60X, respectively. For generation of synthetic normal WGS data per each sample, a germline variant profile from a real patient was down-sampled randomly, representing 50% germline variants of a given patient. These were mixed with the other 50% in silico germline variants that were modelled randomly using an average mutation rate (0.001), finally constituting a full germline profile for normal synthetic WGS data. For generation of synthetic tumor WGS data per each sample, a pre-defined somatic short variant profile (SNVs+Indels) learnt from a real CRC patient was added to the germline variant profile used for creating the normal synthetic WGS data of the same patient, consisting of the variants for tumor sample. Neither copy number profile nor structural variation profile was introduced into the tumor synthetic WGS data. Tumor content and ploidy were assumed to be 100% and 2, respectively. For mapping/variant detection, the Sarek pipeline v3.1.2 (https://nf-co.re/sarek/3.1.2) was used, specifically: 1. BWA v0.7.17-r1188 for read mapping 2. GATK v4.3.0.0 for pre-processing BAM file (including markduplicates and recalibration). 2. Mutect2 (GATK v4.3.0.0) for somatic variant calling 3. Strelka2 v2.9.10 for germline and somatic variant calling Metadata information of 10 CRC patients used for the generation of synthetic normal and tumor WGS data: Patient_id Tumor_barcode Normal_barcode Age Sex Tissue Cancer SIM007 SIM007_T SIM007_N 71 F Rectal Primary CRC SIM008 SIM008_T SIM008_N 45 F Colon Neuroendocrine Metastasis CRC SIM010 SIM010_T SIM010_N 62 M Colon Metastasis CRC SIM011 SIM011_T SIM011_N 55 M Colon Neuroendocrine Metastasis CRC SIM012 SIM012_T SIM012_N 57 M Rectal Metastasis CRC SIM013 SIM013_T SIM013_N 69 M Colon Metastasis CRC SIM014 SIM014_T SIM014_N 68 M Colon Neuroendocrine primary CRC SIM015 SIM015_T SIM015_N 58 F Colon Primary CRC SIM016 SIM016_T SIM016_N 49 M Colon/Rectal Primary CRC SIM017 SIM017_T SIM017_N 78 M Colon Neuroendocrine primary CRC

  Dataset EGAD50000000564

• Cell type mapping of inflammatory muscle diseases highlights selective myofiber vulnerability in inclusion body myositis

  Inclusion body myositis (IBM) is the most prevalent inflammatory muscle disease in older adults with no effective therapy available. In contrast to other inflammatory myopathies like subacute immune-mediated necrotizing myopathy (IMNM), IBM follows a chronic disease course with both inflammatory and degenerative features of pathology. Moreover, causal factors and molecular drivers of IBM progression are largely unknown. Therefore, we paired single-nucleus RNA sequencing with spatial transcriptomics from patient muscle biopsies to map cell type-specific drivers underlying IBM pathogenesis compared to IMNM muscles and non-inflammatory skeletal muscle samples. In IBM muscles, we observed a selective loss of type 2 myonuclei paralleled by increased levels of cytotoxic T and cDC1 cells. IBM myofibers were either characterized by upregulation of cell stress markers featuring GADD45A and NORAD, or protein degradation markers including RNF7 associated with p62 aggregates. GADD45A upregulation was preferentially seen in type 2A myofibers associated with severe tissue inflammation. We also noted IBM-specific upregulation of ACHE encoding acetylcholinesterase, which can be regulated by NORAD activity and result in functional denervation of myofibers. Our results provide promising insights into possible mechanisms of myofiber degeneration in IBM and suggest a selective type 2 fiber vulnerability linked to genomic stress and denervation pathways.

  Study EGAS50000000310

• Reconstructing oral cavity tumor evolution through brush biopsy

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

  Study EGAS50000000602

• IMPRESS: Improved methylation profiling using restriction enzymes and smMIP sequencing, combined with a new biomarker panel, creating a multi-cancer detection assay

  Background Despite the worldwide progress in cancer diagnostics, more sensitive diagnostic biomarkers are needed. The methylome has been extensively investigated in the last decades, but a low-cost, bisulfite-free detection method for multiplex analysis is still lacking. Methods We developed a methylation detection technique called IMPRESS, which combines methylation-sensitive restriction enzymes and single-molecule Molecular Inversion Probes. We used this technique for the development of a multi-cancer detection assay for eight of the most lethal cancer types worldwide. We selected 1791 CpG sites that can distinguish tumor from normal tissue based on DNA methylation. These sites were analysed with IMPRESS in 35 blood, 111 tumor and 114 normal samples. Finally, a classifier model was built. Results We present the successful development of IMPRESS and validated it with ddPCR. The final classifier model discriminating tumor from normal samples was built with 358 CpG target sites and reached a sensitivity of 0.95 and a specificity of 0.91. Moreover, we provide data that highlight IMPRESS’s potential for liquid biopsies. Conclusions We successfully created an innovative DNA methylation detection technique. By combining this method with a new multi-cancer biomarker panel, we developed a sensitive and specific multi-cancer assay, with potential use in liquid biopsies.

  Study EGAS50000000624

• Full characterization of structural variation

  Structural variants (SVs) are important contributors to human disease. Their characterization remains however difficult due to their size and association with repetitive regions. Long-read sequencing (LRS) and optical genome mapping (OGM) can aid as their molecules span multiple kilobases and capture SVs in full. In this study, we selected six individuals who presented with unresolved SVs. We applied LRS onto all individuals and OGM to a subset of three complex cases. LRS detected and fully resolved the interrogated SV in all samples. This enabled a precise molecular diagnosis in two individuals. Overall, LRS identified 100% of the junctions at single-basepair level, providing valuable insights into their formation mechanisms without need for additional data sources. Application of OGM added straightforward variant phasing, aiding in the unravelment of complex rearrangements. These results highlight the potential of LRS and OGM as follow-up molecular tests for complete SV characterization. We show that they can assess clinically relevant structural variation at unprecedented resolution. Additionally, they detect (complex) cryptic rearrangements missed by conventional methods. This ultimately leads to an increased diagnostic yield, emphasizing their added benefit in a diagnostic setting. To aid their rapid adoption, we provide detailed laboratory and bioinformatics workflows in this manuscript.

  Study EGAS50000000520

• Identifying New Genetic Subtypes in Follicular Lymphoma

  Follicular lymphoma (FL) exhibits considerable variability in biological features and clinical trajectories across patients. To dissect the diversity of FL, we utilized a Bernoulli mixture model to identify genetic subtypes in 713 pre-treatment tumor tissue samples. Our analysis revealed the existence of five subtypes with unique genetic profiles that correlated with clinicopathological characteristics. The clusters were enriched in specific mutations as follows: CS (CREBBP and STAT6), TT (TNFAIP3 and TP53), GM (GNA13 and MEF2B), Q (quiescent, for low mutation burden), and AR (mutations of mTOR pathway-related genes). The subtype Q was enriched for patients with stage I disease and was associated with a lower proliferative history compared to the other subtypes. The AR subtype was unique in its enrichment for IgM-expressing FL cases and was associated with advanced-stage and more than 4 nodal sites. The existence of subtypes was validated in an independent cohort of 418 samples from the GALLIUM trial. Notably, patients assigned to the TT subtype consistently experienced inferior progression-free survival when treated with immunochemotherapy. Collectively, our findings offer insight into core pathways distinctly linked with each FL cluster and are expected to be informative in the era of targeted therapies.

  Study EGAS50000000435

• Early detection of ovarian cancer using cell-free DNA fragmentomes and protein biomarkers

  Ovarian cancer is the leading cause of death from gynecological cancers worldwide. No effective screening methods exist, and tumors are commonly detected at advanced stages where treatment is much less effective. In this study, we used whole-genome cell-free DNA (cfDNA) fragmentome and protein biomarker (CA-125 and HE4) analyses to evaluate 591 women from the European Union or the United States with ovarian cancer, benign adnexal masses, or without ovarian lesions. Using a machine learning model that incorporated multi-analyte fragmentome data and protein measurements, we detected ovarian cancer with high specificity >99% and sensitivity of 72%, 69%, 87%, and 100% for stages I–IV, respectively (AUC=0.96, 95% CI: 0.94-0.99), including 90% of high grade serous ovarian cancers. At the same specificity, CA-125 alone detected 34%, 62%, 63%, and 100% of ovarian cancers for stages I–IV (p=0.001, two-sided test of equal proportions). Additionally, our approach distinguished benign masses from ovarian cancers with high accuracy (AUC of 0.88, 95% CI=0.83-0.92). These results were validated in an independent population. These findings show that integrated cfDNA fragmentome and protein analyses detect ovarian cancers with high performance, enabling a new accessible approach for noninvasive ovarian cancer screening and diagnostic evaluation.

  Study EGAS50000000484

• A neoadjuvant, phase II trial demonstrates efficacy and tolerability of Talimogene laherparepvec in cutaneous basal cell carcinoma (NeoBCC trial)

  We present a single arm, phase II neoadjuvant trial with the oncolytic virus Talimogene laherparepvec (T-VEC) in 18 patients with difficult-to-resect basal cell carcinomas. Six cycles of T-VEC were applied intralesionally. The primary endpoint of the study, defined as the proportion of the patients, who after six cycles of T-VEC (13 weeks) at the time point of surgery, become resectable with primary wound closure without the need for plastic reconstructive surgery like skin flaps or skin grafts, was achieved (9/18 patients; 50.0%). Secondary endpoints also included the relapse free survival rate (100%), overall survival (100%) and biological read outs in the tumor samples. T-VEC led to a significant increase of cytotoxic T cells, B cells and myeloid cells, and a decrease of Tregs within the tumor microenvironment. We observed a potent induction of antitumoral and antiviral T cell clones, as well as a strong humoral immune response. The objective response rate was 55.6% and the complete pathological response rate was 33.3%. Only mild adverse events occurred. Together, neoadjuvant T-VEC represents a viable treatment option for patients with difficult-to-resect BCCs (EudraCT Number 2018-002165-19).

  Study EGAS50000000252

• snRNA: A Cultured Leukemia Stem Cell Model Enables Validation of CDK6 as a Stemness Target Against Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) represents a group of aggressive hematological malignancies, the clinical management of which is made challenging due to the persistence of rare and therapy resistant leukemia stem cells (LSCs) which serve as a source of disease relapse and poor outcomes. There are currently a paucity of methods to reliably enrich and study LSCs, hindering the development of therapies that specifically target LSCs. In this study, we deeply characterize the OCI-AML8227 culture model, which maintains a functional stemness hierarchy originating from its highly primitive CD34⁺CD38⁻ cells, to elucidate LSC biology and uncover LSC-specific therapeutic vulnerabilities. We analyzed both bulk and single-cell proteomics, transcriptomics, and epigenomics to generate a LSC protein-protein interaction network, which was then integrated with an LSC-focused small molecule screen using this model. From these findings, CDK6 was discovered as a therapeutic vulnerability specific to LSCs, which was validated in findings from the BEAT-AML cohort and a patient-derived xenograft (PDX) panel of AML samples through palbociclib treatment. Taken together, our studies validate CDK6 as a druggable vulnerability in LSCs, and authenticate OCI-AML8227 cells as a LSC target discovery engine.

  Study EGAS50000000786

• bulk ATACseq: A Cultured Leukemia Stem Cell Model Enables Validation of CDK6 as a Stemness Target Against Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) represents a group of aggressive hematological malignancies, the clinical management of which is made challenging due to the persistence of rare and therapy resistant leukemia stem cells (LSCs) which serve as a source of disease relapse and poor outcomes. There are currently a paucity of methods to reliably enrich and study LSCs, hindering the development of therapies that specifically target LSCs. In this study, we deeply characterize the OCI-AML8227 culture model, which maintains a functional stemness hierarchy originating from its highly primitive CD34⁺CD38⁻ cells, to elucidate LSC biology and uncover LSC-specific therapeutic vulnerabilities. We analyzed both bulk and single-cell proteomics, transcriptomics, and epigenomics to generate a LSC protein-protein interaction network, which was then integrated with an LSC-focused small molecule screen using this model. From these findings, CDK6 was discovered as a therapeutic vulnerability specific to LSCs, which was validated in findings from the BEAT-AML cohort and a patient-derived xenograft (PDX) panel of AML samples through palbociclib treatment. Taken together, our studies validate CDK6 as a druggable vulnerability in LSCs, and authenticate OCI-AML8227 cells as a LSC target discovery engine.

  Study EGAS50000000787

• scATAC: A Cultured Leukemia Stem Cell Model Enables Validation of CDK6 as a Stemness Target Against Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) represents a group of aggressive hematological malignancies, the clinical management of which is made challenging due to the persistence of rare and therapy resistant leukemia stem cells (LSCs) which serve as a source of disease relapse and poor outcomes. There are currently a paucity of methods to reliably enrich and study LSCs, hindering the development of therapies that specifically target LSCs. In this study, we deeply characterize the OCI-AML8227 culture model, which maintains a functional stemness hierarchy originating from its highly primitive CD34⁺CD38⁻ cells, to elucidate LSC biology and uncover LSC-specific therapeutic vulnerabilities. We analyzed both bulk and single-cell proteomics, transcriptomics, and epigenomics to generate a LSC protein-protein interaction network, which was then integrated with an LSC-focused small molecule screen using this model. From these findings, CDK6 was discovered as a therapeutic vulnerability specific to LSCs, which was validated in findings from the BEAT-AML cohort and a patient-derived xenograft (PDX) panel of AML samples through palbociclib treatment. Taken together, our studies validate CDK6 as a druggable vulnerability in LSCs, and authenticate OCI-AML8227 cells as a LSC target discovery engine.

  Study EGAS50000000788

• bulk RNA-seq: A Cultured Leukemia Stem Cell Model Enables Validation of CDK6 as a Stemness Target Against Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) represents a group of aggressive hematological malignancies, the clinical management of which is made challenging due to the persistence of rare and therapy resistant leukemia stem cells (LSCs) which serve as a source of disease relapse and poor outcomes. There are currently a paucity of methods to reliably enrich and study LSCs, hindering the development of therapies that specifically target LSCs. In this study, we deeply characterize the OCI-AML8227 culture model, which maintains a functional stemness hierarchy originating from its highly primitive CD34⁺CD38⁻ cells, to elucidate LSC biology and uncover LSC-specific therapeutic vulnerabilities. We analyzed both bulk and single-cell proteomics, transcriptomics, and epigenomics to generate a LSC protein-protein interaction network, which was then integrated with an LSC-focused small molecule screen using this model. From these findings, CDK6 was discovered as a therapeutic vulnerability specific to LSCs, which was validated in findings from the BEAT-AML cohort and a patient-derived xenograft (PDX) panel of AML samples through palbociclib treatment. Taken together, our studies validate CDK6 as a druggable vulnerability in LSCs, and authenticate OCI-AML8227 cells as a LSC target discovery engine.

  Study EGAS50000000789

• Human inflammatory cardiomyopathies following SARS-CoV2 infection and COVID-19 vaccination

  Title: The cellular and molecular responses in human inflammatory cardiomyopathies following SARS-CoV2 infection and COVID-19 vaccination. Abstract: Abstract: Myocarditis is a cardiac disorder with multiple etiologies characterised by inflammatory cell infiltration and can occur after SARS-CoV-2 infection or rarely following mRNA-based COVID-19 vaccination. The underlying cellular and molecular mechanisms driving these pathologies remain poorly understood.Here we performed single-nucleus-RNA-sequencing in left ventricular endomyocardial biopsies from patients with Non-COVID-19 myocarditis, following SARS-CoV-2 infection (Post-COVID-19) and COVID-19 vaccination (Post-Vaccination). We detected specific cytokine expression patterns highlighting a particular role of interferon-γ in Post-COVID-19 and upregulated IL16 and IL18 expression as a Post-Vaccination hallmark. While the myeloid response was similar between groups, CD4+T-cell proportions were higher in Post-Vaccination and cytotoxic CD8+T and NK cells expanded in Post-COVID-19. Endothelial cells showed gene expression changes suggestive of vascular barrier function deficiency in Post-COVID-19 and angiogenesis response to cardiac inflammation in all groups. Together, our results illuminate shared and distinct cellular and molecular architectures of Non-COVID-19, Post-COVID-19 and Post-Vaccination associated myocardial inflammation in human hearts.

  Study EGAS50000000769

• 3D tissue engineered human skeletal muscle modelling Facioscapulohumeral Muscular Dystrophy

  acioscapulohumeral muscular dystrophy (FSHD) is caused by sporadic misexpression of the transcription factor double homeobox 4 (DUX4) in skeletal muscles. So far, monolayer cultures and animal models have been used to study the FSHD disease mechanism and for FSHD therapy development, but these models do not fully recapitulate the disease and there is a lack of knowledge on how DUX4 misexpression leads to skeletal muscle dysfunction. To overcome these barriers, we have developed a 3D tissue engineered skeletal muscle (3D-TESM) model by generating genetically matched myogenic progenitors (MPs) from human induced pluripotent stem cells of three mosaic FSHD patients. 3D-TESMs derived from genetically affected MPs recapitulate pathological features including DUX4 and DUX4 target gene expression, smaller myofiber diameters, and reduced absolute forces upon electrical stimulation. RNA sequencing data illustrates increased expression of DUX4 target genes in 3D-TESMs compared to 2D myotubes, and cellular differentiation was improved by 3D culture conditions. Treatment of 3D-TESMs with small molecules identified in drug development screens in 2D muscle culture showed no improvements in contractile force and sarcomere organization, highlighting the potential benefit of 3D-TESMs for preclinical drug potency screening by eliminating false-positive compounds for FSHD research and drug development.

  Study EGAS50000000502

• Use of new methodologies to achieve a thorough molecular characterization in pediatric acute leukemia

  In childhood leukemia, a better genetic characterization is needed to define new clinic-biological entities, refine the risk-stratification of patients and help to reduce the toxicity by applying a risk-directed therapy. This project aims to obtain information that will allow us to personalize the management of pediatric patients with leukemia as much as possible, so that we can offer a complete diagnosis, a sensitive follow-up and a more targeted treatment to the real risk of each case. The use of massive genome analysis technologies such as transcriptome sequencing (RNA-Sequencing) will allow a comprehensive biological characterization of pediatric leukemia patients, especially those cases that we have not been able to characterize by conventional diagnostic techniques. We believe that the RNA-Seq technique will allow us to identify new biomarkers that will be useful for: 1) defining new biological entities with their own clinical behavior; 2) refining the stratification of patients into risk groups; 3) finding targets for a better follow-up of MRD and 4) identifying therapeutic targets for the application of targeted therapies to reduce toxicity. PI: Camós-Guijosa, Mireia; Vega-García, Nerea Funding: PI21/00213 Instituto de Salud Carlos III (ISCIII)

  Study EGAS50000000701

• APOBEC3 mutagenesis drives therapy resistance in breast cancer

  Acquired genetic alterations commonly drive resistance to endocrine and targeted therapies in metastatic breast cancer, however the underlying processes engendering these diverse alterations are largely uncharacterized. To identify the mutational processes operant in breast cancer and their impact on clinical outcomes, we utilized a well-annotated cohort of 3,880 patient samples with paired tumor-normal sequencing data. The mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) enzymes were highly prevalent and enriched in post-treatment compared to treatment-naïve hormone receptor-positive (HR+) cancers. APOBEC3 mutational signatures were independently associated with shorter progression-free survival on antiestrogen plus CDK4/6 inhibitor combination therapy in patients with HR+ metastatic breast cancer. Whole genome sequencing (WGS) of breast cancer models and selected paired primary-metastatic samples demonstrated that active APOBEC3 mutagenesis promoted resistance to both endocrine and targeted therapies through characteristic alterations such as RB1 loss-of-function mutations. Evidence of APOBEC3 activity in pre-treatment samples illustrated a pervasive role for this mutational process in breast cancer evolution. The study reveals APOBEC3 mutagenesis to be a frequent mediator of therapy resistance in breast cancer and highlights its potential as a biomarker and target for overcoming resistance.

  Study EGAS50000000875

• Hematopoietic Tet2 inactivation enhances the response to checkpoint blockade immunotherapy

  Somatic mutations inactivating TET2 are among the most common drivers of clonal hematopoiesis (CH). While TET2 inactivation is associated with monocyte-derived inflammation and improved chimeric antigen-receptor-T cell function, its impact on immunotherapy response is unknown. In murine models, hematopoietic Tet2 mutation enhanced the immune checkpoint blockade (ICB) response via the combined presence of phagocytes, CD4 and CD8 T cells. This effect was lost with myeloid- or T-cell restricted Tet2 inactivation, or in mice with 20% Tet2-mutant hematopoiesis. Mechanistically, in Tet2-mutant tumor-infiltrating leukocytes (TILs), ICB preferentially restricted cell states linked to tumor progression while inducing anti-tumor states. Tet2-mutant monocytes activated costimulatory programs, while Tet2-mutant T cells showed enhanced T cell memory signatures, alongside decreased exhaustion and regulatory phenotypes. This murine data was clinically relevant, since tumors from colorectal cancer and melanoma patients with TET2 driver mutation-CH (TET2-CH) showed enhanced immune infiltration, inflammation, and T cell activation. In melanoma patients treated with ICB, TET2-CH was associated with 6-fold greater odds of clinical benefit. Collectively, we demonstrate that hematopoietic Tet2 inactivation primes leukocytes for anti-tumor states associated with immunotherapy response and provides a potential biomarker for personalized therapy.

  Study EGAS50000001191

• Drug Development against Tumor Microtube Networks in Glioblastoma

  Glioblastomas are aggressive brain tumors that resist current therapies by hijacking neurodevelopmental programs to invade the brain and form resilient multicellular networks. Tumor Microtubes (TMs), neurite-like processes extended by cancer cells, drive this pathobiology, but TM-inhibiting drugs are lacking. Here we developed a comprehensive, combined in vitro and in vivo drug screening approach, including novel AI-based analysis tools of TM parameters. Two Protein Kinase C (PKC) activators, TPPB and PMA, exhibited robust inhibition of TM formation in cell-based assays. Since unconnected glioblastoma cells exhibit increased sensitivity to cytotoxic therapy in vivo, brain tumor-bearing mice received TPPB, the most consistent TM-inhibiting compound, in combination with radiotherapy. Long-term intravital 2-photon microscopy confirmed anti-TM and anti-tumor effects of the co-treatment. Mechanistically, TPPB treatment inhibited several key molecular and functional components of TM-related glioblastoma cell network communication and decreased the expression of TTYH1, a key driver of invasive TMs. By establishing a novel screening pipeline for anti-TM drug development, our study identifies a druggable TM master regulator pathway and a small molecule with robust anti-TM effects. It provides a proof-of-principle that disrupting TMs and disconnecting tumor networks is a pharmacologically feasible strategy to make therapies against glioblastoma more effective.

  Study EGAS50000000477

• Recurrent intra-tumour heterogeneity is a hallmark of metastatic prostate cancer

  The evolution of cancers from low grade to metastatic disease is a major determinant of cancer mortality. Cancer evolution involves a complex interplay between cell intrinsic genetics and transcriptional alterations and the extrinsic microenvironment. To define the key mechanisms underpinning metastatic development, we focused on the most lethal form of prostate cancer, metastatic castration-resistant prostate cancer, and employed single-cell multi-omics and whole-genome sequencing to deeply profile 34 metastatic lesions from 9 patients obtained from rapid autopsy. We found that intra-tumour heterogeneity is an indicator of key evolutionary processes, characterised by recurrent tumour populations acting as critical functional components of the tumour ecosystem, irrespective of clonal and microenvironmental backgrounds. Unexpectedly, microenvironments across metastatic sites played a limited role while clonal evolution primarily promoted transcriptional noise. Intra-patient functional convergence of tumour ecosystems was observed across metastatic sites, pointing towards system-level selection pressures that drive the heterogeneity landscape of mCRPC. Our findings reveal functional evolutionary convergence of metastatic disease into units of intra-tumour heterogeneity identifying critical determinants of metastatic progression for therapeutic targeting. Note: matched WGS for some tumour single-cell experiments, as well as germline WGS for patients CA0027, CA0034, CA0035, CA0043, and CA0076, can be found in EGAS00001006598.

  Study EGAS50000001312

• Bruno et al.: Interferon gamma rebalances immunopathological signatures in Chronic Granulomatous Disease through metabolic rewiring

  Chronic granulomatous disease (CGD) is a primary immunodeficiency characterized by recurrent life-threatening infections and hyperinflammatory complications. It is caused by mutations in the NADPH oxidase complex and the consequent loss of reactive oxygen species (ROS) production. Recombinant human interferon gamma (rIFN-γ) prophylaxis reduces the risk of severe infections, but the mechanisms behind its efficacy in CGD are still an open question, as it does not restore NADPH oxidase-dependent ROS production. Here, we show that myeloid innate immune cells of CGD patients are transcriptionally and functionally reprogrammed to a hyperactive inflammatory status, displaying an impaired in vitro induction of trained immunity. CGD monocytes have reduced intracellular amino acids concentrations and profound functional metabolic defects, both at the level of glycolysis and mitochondrial respiration. Ex vivo and in vivo treatment with IFN-γ restored these metabolic defects and reduced excessive IL-1β and IL-6 production in response to fungal stimuli in CGD monocytes. These data suggest that prophylactic rIFNγ modulates the metabolic status of innate immune cells in CGD. These data shed light on the effects of NADPH-derived ROS deficiency to the metabolic programs of immune cells and pose the basis for targeting this immunometabolic axis, potentially beyond CGD, with IFN-γ immunotherapy.

  Study EGAS00001005463

• Phylogenetic analysis of paired breast carcinomas identifies genetic events associated with clonal recurrence and invasive progression

  Development of ipsilateral breast carcinoma following diagnosis of breast ductal carcinoma in situ (DCIS) has been assumed to represent recurrence of the primary tumour. However, this may not be the case and it is important to know how often recurrences are new tumours. Ipsilateral primary-recurrence pairs (n=78) were sequenced to test their clonal relatedness. Shared genetic events were identified from whole exome sequencing (n=54 pairs) using haplotype-specific copy number and phylogenetic analysis. The remaining pairs were sequenced by a targeted panel or low-coverage whole genome sequencing. We included 32 non-recurrent DCIS to compare recurrent and non-recurrent disease. We found that 7% of DCIS recurrences were non-clonal by whole exome sequencing, indicative of a new breast carcinoma. Lower resolution methods detected a higher non-clonality rate (29%). Comparing primary DCIS with their recurrences found that evolution of DCIS to invasive disease was associated with increased ploidy and copy number events. TP53 mutations were enriched in DCIS with clonal recurrence compared with non-recurrent DCIS. Our results verify that de novo “recurrent tumours” of independent origin occur in patients who may be at high risk. (https://www.biorxiv.org/content/10.1101/2024.05.19.594731v2.full)

  Study EGAS50000001298

• Emergence of oncofetal plasticity is ubiquitous in early colorectal cancers

  Metastasis formation is classically considered a late-stage event in colorectal cancer (CRC) evolution. Yet, the time and patterning by which metastatic competence is acquired remain poorly understood. Here, we show that metastasis-associated oncofetal cell states already emerge at the earliest stages of CRC, concurrent with invasive front formation. However, while necessary for metastasis, we detect them ubiquitously among early non-metastatic cancers, highlighting additional bottlenecks like immune evasion. To understand how oncofetal cells first emerge, we generated multiregional organoid models that reflect successive tumor progression stages within individual early-stage CRCs. Whole-genome sequencing and growth factor-dependency assays exclude tumor cell-intrinsic acquired traits. In contrast, single-cell spatial atlases of the tumor-microenvironment before and after malignant transformation revealed stereotypic patterning of fibroblast subtypes resembling normal tissue architecture, resulting in distinct regional microenvironments. At the onset of malignant growth into the submucosa, the first cancer-associated fibroblasts (CAFs) to appear strongly resemble submucosal trophocytes and colocalize with oncofetal cell states at invasive fronts. Functionally, fibroblast-organoid co-cultures confirm that these trophocyte-like CAFs induce plastic transitioning to oncofetal states. Thus, interactions between tumor and submucosal fibroblasts directly following malignant transformation dictate the timing and location at which oncofetal plasticity first occurs during CRC progression.

  Study EGAS50000001532

• Assessing gene expression profiling from FFPE Patient Samples: A Comparison of Two Library Preparation Approaches and Recommendations

  Next-Generation Sequencing (NGS) has transformed cancer research and clinical practice, with Whole Exome Sequencing (WES) driving advances in mutational profiling and personalized oncology. Yet, transcriptomic signatures remain essential for understanding disease mechanisms, including therapy resistance pathways. RNA sequencing (RNA-seq), however, faces unique challenges when dealing with low-input or degraded RNA, as often found in archival formalin-fixed paraffin-embedded (FFPE) tissues. Although previous studies have compared library preparation protocols, rapidly evolving technologies call for updated evaluations. Here, we present a direct comparison of two FFPE-compatible stranded RNA-seq library preparation kits: TaKaRa™ SMARTer® Stranded Total RNA-Seq Kit v2 (Kit A) and Illumina Stranded Total RNA Prep Ligation with Ribo-Zero Plus (Kit B). Both kits generated high-quality RNA-seq data, yet important differences emerged. Notably, Kit A achieved comparable gene expression quantification to Kit B while requiring 20-fold less RNA input, a crucial advantage for limited samples, albeit with increased sequencing depth. We critically discuss these results in relation to RNA availability, technical performance, cost-effectiveness, processing time, and automation potential, offering practical guidance for selecting optimal RNA-seq strategies in clinical and translational research settings

  Study EGAS50000001066

• Myeloid cell networks govern re-establishment of original immune landscapes in recurrent ovarian cancer

  Immunotherapy has produced disappointing results in recurrent ovarian cancer (OC). However, the prognostic value of tumour-infiltrating lymphocytes (TILs) is largely based on the analysis of treatment-naive tumours. To understand the immunobiology of recurrent cancers, and their evolution, we profiled a total of 697 OC tumour samples (566 primary and 131 recurrent including 69 patient-matched samples) from 595 patients of five independent multi-institutional cohorts. By capturing heterogeneous TIL distributions, we identified four immune phenotypes associated with differential prognosis, TILs states and TILs:myeloid networks, which dictate malignant progression after chemotherapy and recurrence. Notably, recurrent tumours recapitulate the immunogenic landscapes of original cancers. Mirroring inflamed human OC, preclinical recurrent Brca1mut tumours maintained activated TILs:dendritic cells (DCs) niches but evaded immune control through COX/PGE2 signaling upregulation. Conversely, recurrent Brca1wt tumours displayed loss of TILs:DCs niches and accumulated immunosuppressive myeloid networks featuring Trem2/ApoEhigh TAMs and Nduf4l2high/Galectin3high malignant states. Our study highlights that persistent immunogenicity in recurrent OC is governed by the BRCA-dependent crosstalk between dissimilar myeloid cells and TILs. We propose new targets to modulate the malignant and TAM compartments which improve efficacy of standard of care chemotherapy in mouse OC models with HR deficiency.

  Study EGAS50000001069

• Finding the way towards the eradication of therapy-related myeloid neoplasms

  Therapy-related neoplasms are directly caused by the drugs used to eradicate a previous cancer. It is well known that a subset of drugs designed to kill cancer cells through DNA-damage-induced apoptosis can intrinsically, as secondary effect, cause DNA-damage in other cells, thus inducing a neoplasm. Those neoplasms have especially bad prognosis. Why are these drugs still in use? Because for many diseases we do not have anything better, at the moment. Approximately 0.5-1% of the kids treated for various types of cancers develop a therapy-related myeloid neoplasm (tMN) as a result of the treatment. Despite having the same origin, their prognostic and molecular features are different from the adult ones, and a complete understanding of the pathogenesis of tMN in children is still missing. A total of 84 paediatrics neoplasms have been profiled by either Whole Genome or Whole Exome sequencing by Jason R. Schwartz and his colleagues, in a collaborative effort guided at the St. Jude Children's hospital, in Memphis. They identified a mean of 28 somatic mutations per patient, a number which is significantly bigger than the number of mutations found in paediatrics primary neoplasms, as depicted in the original paper's figure 1b, shown here. Their sequencing data identified that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements as the most frequent drivers of the disease. Clonal evolution analysis (the comparative study of the behaviour of cells harbouring different mutations) highlighted that, unlike in adults’ tMN, there were not pre-existing leukemic clones, and virtually all mutated cells were found arising as a direct consequence of cytotoxic therapy. This represents yet another important difference between adult and paediatric tMN. This important study have been published in Nature communications in February 2021. The sequencing data have been deposited at EGA under the study number EGAS00001004850. Thanks to this, other researchers around the world could re-use it to uncover more information, integrate and validate their studies. Revealing the molecular features underlying the development of cancers is the only way towards the development of new generation targeted therapies, which will allow us to finally drop the cytotoxic drugs eliciting such unwanted and catastrophic secondary effects as tMN. It is not an overstatement to say that, while improving the treatment of primary tumours, we could entirely eradicate therapy-related neoplasms right from their origin. Looking forward.

  Blog eradication-of-therapy-related-myeloid-neoplasms

• Genome-wide DNA-methylation assessment by MethylCap-seq and Infinium HumanMethylation450 BeadChips: an independent large-scale comparison

  DNA-methylation is an important epigenetic feature in health and disease. Two cost-efficient genome-scale methodologies to assess DNA-methylation are MethylCap-seq and Illumina's Infinium HumanMethylation450 BeadChips (HM450). However, objective information regarding the best-suited methodology for a specific research question is scant. Therefore, we performed a large-scale evaluation on a set of 70 brain tissue samples obtained from 65 glioblastoma and 5 non-tumoral brain tissues, using a gold standard free Bayesian modeling procedure. While conditional specificity was adequate for both approaches, conditional sensitivity was systematically higher for HM450. Also the genome-wide characteristics were compared, revealing that the HM450 probes assess less than 10% of the regions identified as methylated by MethylCap-seq. Hence the latter method may detect more potentially relevant DNA-methylation, defined by either functional location or previously reported differentially methylated candidate regions. Our results therefore indicate that – at least for the tissue under study - both methodologies are complementary, with a higher sensitivity for HM450, but a far larger genome-wide coverage for MethylCap-seq. Note that here only the relevant MethylCap-seq data is deposited, for the HM450 data we refer to GEO (GSE60274).

  Study EGAS00001001191

• UK10K NEURO ABERDEEN

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This sample set comprises cases of schizophrenia with additional cognitive measurements, collected in Aberdeen, Scotland.For further information on the Aberdeen cohort please contact David St Clair (d.stclair@abdn.ac.uk).

  Study EGAS00001000109

• Single_Cell_Sequencing_of_Sperm__scSperm_

  The aim of this project is to genotype and sequence single spermatozoa from two men, one in his twenties and the other in his seventies. The resulting data is used to quantify the mutations that have arisen in the gametes of both individuals in order to better understand the effect of aging on mutation rates and modes.Project Outline. In order to quantify mutations, semen from two individuals are sequenced. 48 single sperm cells are isolated from each individual, and their DNA is extracted. The resulting genomes are amplified using PicoPlex, GenomiPhi MDA, Repli-G MDA, and MALBAC. QC step is applied to check the quality of WGA DNA using standard Sequenom plex (26 SNPs). A subset of 32 amplification products which pass the intiall QC, are genotyped using Affymetrix SNP6 chips. 12 of the genotyped amplification products are also sequenced. In addition, one multi-cell sample per individual is sequenced as a reference and for validation purposes.Altogether, 12 single cell sperm genomes and two multi-cell genomes are sequenced, coming to a total of 14 genomes. Of the single cell sperm genomes, 2 are sequenced to 50x coverage, and the other 10 to 25x coverage. Both multi-cell genomes are sequenced to 25x coverage.

  Study EGAS00001000935

• The Genomic Landscape of Core-Binding Factor Acute Myeloid Leukemias

  Pediatric de novo acute myeloid leukemia (AML) is a heterogeneous disease that can be divided into clinically distinct subtypes based on the presence of specific chromosomal abnormalities or gene alterations. One of the best characterized subtypes of AML involves leukemias with alterations of the core-binding factor (CBF)-complex, which comprises the FAB subtypes M2 and M4Eo and associates with a favorable outcome. Patients with the AML M2 subtype harbor a translocation between chromosomes 8 and 21 [t(8;21)] that yields the chimeric fusion gene RUNX1(AML1)-RUNX1T1(ETO), while patients with AML M4Eo express the chimeric fusion gene CBFβ-SMMHC(MYH11) as a result of an inversion/translocation event of chromosome 16 [inv(16)/t(16;16)]. In an effort to define the total complement of genetic changes in CBF-leukemia, we performed paired-end whole genome sequencing (WGS) on diagnostic leukemia blasts and matched germ line samples from 17 pediatric CBF-leukemia patients using the Illumina platform. Somatic alterations, including single nucleotide variations (SNVs) and structural variations (SVs), including insertions, deletions, inversions, and inter- and intra-chromosomal rearrangements, were detected using complementary analysis pipelines (Bambino, CREST and CONSERTING). Recurrent screening of identified mutations will be performed in a cohort of approximately 94 cases of CBF-leukemias.

  Study EGAS00001000349

• ERG ALTERATIONS DEFINE A NOVEL SUBTYPE OF ACUTE LYMPHOBLASTIC LEUKEMIA

  Structural chromosomal alterations are a hallmark of acute lymphoblastic leukemia (ALL) yet many patients lack an identifiable abnormality on conventional cytogenetic analysis. Here, using integrated analysis of genome-wide DNA copy number and gene expression data of 1764 childhood and adult ALL cases, including whole-genome, exome and mRNA-sequencing of 72 cases, we report a novel subtype of B-progenitor ALL (B-ALL) representing 4-13% of B-ALL cases characterized by a distinct gene expression profile, and focal deletions of ERG (ETS related gene) in 56% of cases. These cases are characterized by expression of a ERG transcript that utilizes a novel exon in intron 6 of ERG spliced to the canonical reading frame of ERG exons 7-10, encoding a 27 kDa C-terminal ERG protein that retains the ETS and transactivating domains but lacks the N-terminal pointed and regulatory domains. This protein is a competitive inhibitor of wild-type ERG and is leukemogenic in vivo. ERG-altered ALL cases have a favourable outcome despite the presence of alterations associated with poor prognosis in non-ERG ALL, such as deletions of IKZF1. These findings identify a new subtype of leukemia characterized by a novel mechanism of ETS transcription factor dysregulation in cancer.

  Study EGAS00001000514

• An RCOR1 loss-associated gene expression signature identifies a prognostically significant DLBCL subgroup

  Effective treatment of diffuse large B cell lymphoma (DLBCL) is plagued by heterogeneous responses to standard therapy, and molecular mechanisms underlying unfavorable outcomes in lymphoma patients remain elusive. Here, we profiled 148 genomes with 91 matching transcriptomes in an R-CHOP-treated DLBCL cohort to uncover molecular subgroups linked to treatment failure. Systematic integration of high-resolution genotyping arrays and RNA-seq data revealed novel deletions in RCOR1 to be associated with unfavorable progression free survival (P = 0.001). Integration of expression data from the clinical samples with data from RCOR1 knockdowns in the lymphoma cell lines KM-H2 and Raji yielded an RCOR1 loss-associated gene signature comprising of 233 genes. This signature identified a subgroup of patients with unfavorable overall survival (OS) (P = 0.023). The prognostic significance of the 233-gene signature for OS was reproduced in an independent cohort comprising of 195 R- CHOP treated patients (P = 0.039). Additionally, we discovered that within the IPI low risk group, the gene signature provides additional prognostic value that was independent of the cell-of-origin phenotype. Taken together, we present a novel and reproducible molecular subgroup of DLBCL, impacting risk-stratification of R- CHOP treated DLBCL patients and revealing a new avenue for therapeutic intervention strategies.

  Study EGAS00001001000

• Accurate sample assignment in a multiplexed, ultra-sensitive, high-throughput sequencing assay for minimal residual disease

  High throughput sequencing (HTS) (next generation sequencing) of the rearranged immunoglobulin and T-cell receptor genes promises to be cheaper and more sensitive than current methods for monitoring minimal residual disease (MRD) in patients with acute lymphoblastic leukemia. However, adoption of new approaches by clinical laboratories requires careful evaluation of all potential sources of error and the development of strategies to ensure the highest accuracy. Timely and efficient clinical use of HTS platforms will depend on combining multiple samples (multiplexing) in each sequencing run. Here we examine immunoglobulin heavy chain gene HTS on the Illumina MiSeq platform for MRD (HTS-MRD). We identify errors associated with multiplexing that could potentially impact on the accuracy of MRD analysis. We optimise a strategy combining high purity, sequence-optimised oligonucleotides, dual-indexing and an error-aware demultiplexing approach to minimise errors and maximise sensitivity. We present a probability-based demultiplexing pipeline, Error-Aware Demultiplexer (EAD) - that is suitable for all MiSeq sequencing strategies and accurately assigns samples to the correct identifier without excessive loss of data. Finally using controls quantified by digital PCR, we show that HTS-MRD can accurately detect as few as 1 in 106 copies of specific leukemic MRD.

  Study EGAS00001001303

• Rare SNPs in receptor tyrosine kinases are negative outcome predictors in multiple myeloma

  Multiple myeloma (MM) is a plasma cell disorder that is characterized by a great genetic heterogeneity. Recent next generation sequencing studies revealed an accumulation of tumor-associated mutations in receptor tyrosine kinases (RTKs) which may also contribute to the activation of survival pathways in MM. To investigate the clinical role of RTK-mutations in MM, we deep-sequenced the coding DNA-sequence of EGFR, EPHA2, ERBB3, IGF1R, NTRK1 and NTRK2 which were previously found to be mutated in MM, in 75 uniformly treated MM patients of the "Deutsche Studiengruppe Multiples Myelom". Subsequently, we correlated the detected mutations with common cytogenetic alterations and clinical parameters. We identified 11 novel non-synonymous SNVs or rare patient-specific SNPs, not listed in the SNP databases 1000 genomes and dbSNP, in 10 primary MM cases. The mutations predominantly affected the tyrosine-kinase and ligand-binding domains and no correlation with cytogenetic parameters was found. Interestingly, however, patients with RTK-mutations, specifically those with rare patient-specific SNPs, showed a significant lower overall, event-free and progression-free survival. This indicates that RTK SNVs and rare patient-specific RTK SNPs are of prognostic relevance and suggests that MM patients with RTK-mutations could potentially profit from treatment with RTK-inhibitors.

  Study EGAS00001001665

• Recurrent somatic JAK-STAT mutations within a novel RUNX1-mutated pedigree

  The acquisition of somatic mutations is an emerging field of investigation in familial leukemia. Currently, genetic profiles in familial MDS/AML are considered analogous to sporadic disease, although the patterns of clonal evolution within families are poorly defined. We performed whole exome profiling of tumour samples from a novel RUNX1 mutated family, to determine the stepwise evolution of MDS/AML across 4 young siblings. Three siblings developed monocytic AML/RAEB2 at 5 years of age, with hepatosplenomegaly and somatic mutations upregulating JAK-STAT signalling, the latter are typically detected in <5% of sporadic MDS/AML. Two siblings acquired the canonical JAK2 V617F mutation, while another acquired a unique missense mutation of SH2B3, a negative regulator of JAK2. Notably, 2/3 siblings demonstrated dosage amplification of these mutations due to acquired uniparental disomy of chromosomes 9p and 12q (encompassing JAK2 and SH2B3, respectively). All 4 siblings were heterozygous for the 46/1 JAK2 haplotype associated with predisposition to sporadic V617F myeloproliferative disorders, which likely influenced the acquisition of JAK2 mutations. Our findings provide further evidence that relatives with shared germline mutations may acquire somatic mutations in a non-random manner leading to convergent patterns of disease evolution.

  Study EGAS00001001862

• Comparison of EGF and PDGF driven glioblastomas.

  Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth-factor (EGF) and platelet-derived growth-factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in these tumors. Despite years of research into tyrosine-kinase inhibitors, clinical trials have shown disappointing efficacy. One of the underlying reasons is intra-tumor heterogeneity. To assess the effect of intra-tumor heterogeneity on gene expression, we derived a general approach to map single-cell expression profiles to sequentially-acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. A deletion in a PDGF-receptor domain involved in dimerization correlated to an up-regulation of the PI3K/AKT pathway in a proneural GBM, and ectopic expression of this mutant PDGFRA in vitro enhanced proliferation and invasiveness. Descending the phylogenetic tree of the PDGFR-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type that expresses growth and angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed a gradual up-regulation of pro-invasive genes.

  Study EGAS00001001900

• Genomics of enteropathy associated T cell lymphoma (EATL)

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

  Study EGAS00001001954

• Whole genome sequencing of tumour and normal paired samples of diffuse intrinsic pontine gliomas

  Diffuse intrinsic pontine glioma (DIPG) have a universally dismal prognosis (median 9-12 months), with neither chemotherapeutic nor targeted agents showing any substantial survival benefit in clinical trials in children with these tumours. DIPG are highly infiltrative malignant glial neoplasms of the ventral pons, which due to their location within the brain, make them unsuitable for surgical resection and consequently have a universally dismal clinical outcome. Recent high-throughput sequencing approaches have revealed a striking prevalence of K27M mutations in the genes encoding the histone variants H3.3 (H3F3A) or H3.1 (HIST1H3B) in the childhood brain tumour DIPG. This K-to-M substitution confers a trans-dominant ablation of global H3K27 trimethylation, which likely profoundly alters gene expression through de-repression of polycomb repressive complex 2 (PRC2) target genes. Despite these advances in our understanding of the distinct biology of these tumours, approaches for specific novel therapeutic interventions are not clear, and little has been reported of the secondary mutations accompanying these changes. We carried out whole genome sequencing on a series of 20 tumour normal pairs of pre-treatment biopsy samples of DIPG, for which the patients underwent a safe stereotactic procedure, and whole exome sequencing on a further six tumour normal pairs obtained at autopsy.

  Study EGAS00001000572

• Pediatric Non-Down Syndrome Acute Megakaryoblastic Leukemia is Characterized by Distinct Genomic Subsets with Varying Outcomes

  Acute Megakaryoblastic Leukemia (AMKL) is a subtype of acute myeloid leukemia (AML) in which cells morphologically resemble abnormal megakaryoblasts.  While extremely rare in adults, AMKL accounts for 4-15% of newly diagnosed childhood AML1-3.  A significant proportion of pediatric AMKL cases occur in children with Down syndrome (DS).  These patients have excellent outcomes and are characterized at the genomic level by a founding GATA1 mutation4-6.  In contrast, AMKL in patients without DS (non-DS-AMKL) is frequently associated with poor outcomes.  Previous efforts have identified chimeric oncogenes in a significant number of cases, including RBM15-MKL1, CBFA2T3-GLIS2, KMT2A gene rearrangements and NUP98-KDM5A7-9.  The etiology of 30-40% of cases, however, remains unknown.  To better understand the genomic landscape of non-DS-AMKL, we performed RNA and exome sequencing on specimens from 99 patients (75 pediatric and 24 adult).  We demonstrate that pediatric non-DS-AMKL is a heterogeneous malignancy that can be divided into seven subgroups with varying outcomes.  These subgroups are characterized by chimeric oncogenes with cooperating mutations in epigenetic and kinase signaling genes.  Overall, these data shed light on the etiology of AMKL and provide useful information for treatment tailoring. 

  Study EGAS00001002183

• Whole-exome and targeted sequencing of pediatric hyperdiploid B-cell precursor acute lymphoblastic leukemia

  High hyperdiploidy defines the largest genetic entity of childhood acute lymphoblastic leukemia (ALL). Despite its relatively low recurrence risk, this subgroup generates a high proportion of relapses. The cause and origin of these relapses remains obscure. We therefore explored the mutational landscape in high hyperdiploid (HD) ALL with whole-exome (n=19) and subsequent targeted deep sequencing of 60 genes in 100 relapsing and 51 non-relapsing cases. We identified multiple clones at diagnosis that were primarily defined by a variety of mutations in receptor tyrosine kinase (RTK)/Ras pathway and chromatin-modifying genes. The relapse clones consisted of reappearing as well as new mutations, and overall contained more mutations. Although RTK/Ras pathway mutations were similarly frequent between diagnosis and relapse, both intergenic and intragenic heterogeneity was essentially lost at relapse. CREBBP mutations, however, increased from initially 18-30% at relapse, then commonly co-occurred with KRAS mutations (P<0.001) and these relapses appeared primarily early (P=0.012). Our results confirm the exceptional susceptibility of HD ALL to RTK/Ras pathway and CREBBP mutations, but, more importantly, suggest that mutant KRAS and CREBBP might cooperate and equip cells with the necessary capacity to evolve into a relapse-generating clone.

  Study EGAS00001001113

• Novel genes for Intellectual Disability identified using whole genome sequence and pathway analysis

  Intellectual Disability (ID) is one of the most common global disorders. However for ~30% of patients with ID no cause is identified, despite a clinical diagnostic odyssey that includes genome wide clinical microarray (CMA). We carefully selected 8 trios, each composed of a child with ID and brain structural defect, and both normal parents for whole genome sequencing (WGS). We conducted WGS on the trio and filtered out de novo single nucleotide variants (SNVs) and then used a pathway-based refinement to select candidates. We confirmed a de novo pathogenic SNV in ARID1B, and likely pathogenic SNVs in CACNB3, SPRY4, PHF6, SQSTM1 and UPF1 in 5 of the 8 children. All genes except ARID1B and PHF6 are previously unreported for ID. We analysed our WGS data using 4 independent algorthims for copy number and structural variation including using de novo whole genome assembly. We confirmed one likely contributory 165kb de novo CNV (missed by CMA). We conducted a validation study of over 2000 exomes for our novel candidate genes and also present a strategy to analyze the non-coding sequence space. This study provides an extensive analysis of WGS in the context of ID and yielded causative variants in >60% of cases.

  Study EGAS00001001386

• Ethiopia_Genome_Project__high_coverage_

  The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

  Study EGAS00001000237

• UK10K_RARE_NEUROMUSCULAR

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The neuromuscular disorder samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Francesco Muntoni (f.muntoni@ucl.ac.uk).

  Study EGAS00001000101

• UK10K_RARE_FIND

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The Raymond samples will be part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Lucy Raymond (flr24@cam.ac.uk).

  Study EGAS00001000128

• UK10K_RARE_CILIOPATHIES

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The Ciliopathies samples will be part of the rare disease group, and will undergo exome sequencing. For further information with regard to this cohort please contact Phil Beales (p.beales@ich.ucl.ac.uk).

  Study EGAS00001000126

• UK10K_RARE_SIR

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes.Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.The SIR (Severe Insulin Resistance) samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Robert Semple (rks16@cam.ac.uk).

  Study EGAS00001000130

• Stratifying and Targeting Pediatric Medulloblastoma through Genomics

  In this project, genomic analyses of pediatric medulloblastoma samples, obtained through the international medulloblastoma consortium, will be performed. RNA and miRNA expression profiles of 1000 samples, representing all four subgroups (Wnt, Shh, Group C, and D), will be studied to identify novel subtypes within each subgroup. The resulting subtype-specific expression profiles will support the development of reliable and robust biomarkers to more accurately and reliably classify medulloblastomas for treatment in clinical trials. For that purpose, two assays will be developed: an antibody-based immunohistochemical assay and an orthogonal nucleic acid-based hybridization assay. Additional genomic DNA analysis of the 300 high risk subgroup cases will support the discovery of subgroup specific somatic mutations in order to inform current clinical trials of targeted therapies, and to identify genes and pathways already targeted in other diseases. Such therapies could be rapidly transitioned to Phase II trials in medulloblastoma. Furthermore, the discovery of somatic mutations could be used for developing as well as validating specific biomarkers. The project team will also try to identify risk factors that predispose children to this type of cancer.

  Study EGAS00001000273

• UK10K RARE CHD

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The CHD (congenital heart disease) samples will be part of the “rare disease” group, and will undergo exome sequencing.For further information with regard to this cohort please contact Shoumo Bhattacharya (shoumo@me.com).

  Study EGAS00001000125

• Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR02

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Paired DNA samples – isolated from the tumour tissue and from the blood of fifty six patients – have been analyzed to catalog germline and somatic mutations. Association between each observed genomic alteration and exposure to environmental risk factors is being explored.Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumour tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Exome capture and deep resequencing are being performed at the National Institute of Biomedical Genomics, Kalyani. We have analyzed the exomes of paired blood and tumour DNA samples using exome capture by Illumina TrueSeq Exome Enrichment kit followed by sequencing on Illumina HiSeq-2000. Each exome is being sequenced at a mean depth of 70x for Tumour and 40x for Blood. In addition, we have used Illumina Omni 2.5 SNP-chips to generate genotype data.

  Study EGAS00001001028

• P647_Targeted_resequencing_project

  In collaboration with Dr Robert Semple we have identified a family harbouring an autosomal dominant variant, which leads to severe insulin resistance (SIR), short stature and facial dysmorphism. This family is unique within the SIR cohort in having normal lipid profiles, preserved adiponectin and normal INSR expression and phosphorylation. DNA is available for 7 affected and 7 unaffected family members across 3 generations. All 14 samples have been genotyped using microsatellites and the Affymetrix 6.0 SNP chip. Linkage analysis identified an 18.8Mb haplotype on chromosome 19 as a possible location of the causative variant. However, Exome sequencing of 3 affected and 1 unaffected family members has not identified the causative variant suggesting the possibility of an intronic or intergenic variant in this region or elsewhere in the genome. We propose to conduct whole genome sequencing of 5 members of the pedigree at a depth of 20X. The chosen samples are two sets of parents plus one member of an unaffected branch of the pedigree who shares the risk haplotype on chromosome 19. Sequencing of the two sets of parents will be used along with the genome-wide SNP data to impute 4 affected children giving an effect sample size of 6 affected individuals.

  Study EGAS00001000305

• The Genomic Landscape of Childhood and Adolescent Melanoma

  Despite remarkable advances in the genomic characterization of adult melanoma, the molecular pathogenesis of pediatric melanoma remains largely unknown. We analyzed 15 conventional melanomas (CM), 3 melanomas arising in congenital nevi (CNM), and 5 spitzoid melanomas (SM), using various platforms, including whole genome or exome sequencing, molecular inversion probe assay, and/or targeted sequencing. CM demonstrated a high burden of somatic single nucleotide variations (SNV), with each case containing a TERT promoter (TERT-p) mutation, 13/15 containing an activating BRAF V600 mutation, and >80% of the identified SNV consistent with UV damage. By contrast, the 3 CNM contained an activating NRAS Q61 mutation and no TERT-p mutations. SM were characterized by chromosomal rearrangements resulting in activated kinase signaling in 40%, and an absence of TERT-p mutations except for the one SM that succumbed to hematogenous metastasis. We conclude that pediatric CM has a very similar UV-induced mutational spectrum to that found in the adult counterpart, emphasizing the need to promote sun protection practices in early life and to improve access to therapeutic agents being explored in adults in young patients. By contrast, the pathogenesis of CNM appears to be distinct. TERT-p mutations may identify the rare subset of spitzoid melanocytic lesions prone to disseminate.

  Study EGAS00001000901

• Developmental_Dysplasia_of_the_Hip__DDH_

  The National Joint Registry for England and Wales (NJR) has become established as one of the largest repositories of clinical information about patients undergoing hip and knee replacement. The primary aim of this proposal is to examine the feasibility of adding value to that clinical information by partnering it with a DNA dataset to allow its use to be extended to the study of the underlying pathogenesis of diseases of the hip and knee, and the complications of their treatment. We will do this 1) by establishing a mechanism for accessing and linking pseudo-anonymised demographic and phenotype data between the NJR database and the proposed DNA Biobank; 2) by validating the collection of postal DNA samples from consenting NJR patients; and 3) by genotyping genome-wide 900 patients with the condition developmental dysplasia of the hip (DDH) to find variants that play a role in this heritable condition. The successful demonstration that the NJR can be used to establish a DNA Biobank will provide a powerful platform to examine the genetic basis for a wide range of joint diseases, and the complications of their treatment. This proposal is to carry out the genotyping portion of this project on the current CoreExome array

  Study EGAS00001000916

• Comprehensive genomic profiles of small cell lung cancer

  We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. We found bi-allelic inactivation of TP53 and RB1 in nearly all the tumors analyzed, sometimes by complex genomic rearrangements. Two tumors with wild-type RB1 had evidence of chromothripsis leading to overexpression of Cyclin D1, revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumor suppressors TP53 and RB1 is obligatory in SCLC. We further discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In few cases, SCLC tumors exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signaling in a pre-clinical SCLC mouse model dramatically reduced the number of tumors and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

  Study EGAS00001000925

• Exploring the role of mtDNA variation in Multiple Sclerosis in a large cohort of discordant monozygotic twins

  Several lines of evidence indicate that mitochondrial DNA (mtDNA) variants might predispose to multiple sclerosis (MS). We examined this in 49 monozygotic (MZ) twin pairs clinically discordant for MS at study entry. Since the nuclear DNA of MZ twins is identical, our study provides a unique setting to functionally associate unique mtDNA variants and skewed heteroplasmy with MS development. In order to identify mtDNA variants in the twin cohort, we deeply sequenced the mitochondrial genome in blood of all 49 MZ twin pairs using next generation sequencing (Illumina, HiSeq or MiSeq) with an average coverage depth of ~25,000 sequences per base. Heteroplasmic variants were validated using targeted deep sequencing (TDS) (Illumina, MiSeq) in DNA isolated from blood and buccal swabs. In addition, whole blood was also available of 5 additional MS-affected or healthy siblings and 1 dizygotic twin pair discordant for MS, comprising in total 6 non-identical sib pairs. Moreover, all cases were screened for pathogenic mtDNA polymerase gamma (POLG) mutations. The data provides valuable insights in the possible involvement of mtDNA variants in the pathogenesis of MS, and gives information on the segregation of heteroplasmic variants within MZ twins and non-identical siblings and across different tissues.

  Study EGAS00001001240

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

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

  Study EGAS00001001810

• A somatic reference standard for cancer genome sequencing with COLO829

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

  Study EGAS00001001385

• The HLA-DQβ1 insertion is a strong achalasia risk factor and displays a geospatial north-south gradient among Europeans

  Idiopathic achalasia is a severe motility disorder of the esophagus and is characterized by a failure of the lower esophageal sphincter to relax due to a loss of neurons in the myenteric plexus. Most recently, we identified an 8 amino acid insertion in the cytoplasmic tail of HLA-DQβ1 as strong achalasia risk factor in a sample set from Central Europe, Italy, and Spain. Here, we tested whether the HLA-DQβ1 insertion also confers achalasia risk in the Polish and Swedish population. We could replicate the initial findings and the insertion shows strong achalasia association in both samples (Poland P = 1.84×10-04, Sweden P = 7.44×10-05). Combining all five European data sets – Central Europe, Italy, Spain, Poland, and Sweden – the insertion is achalasia associated with Pcombined = 1.67×10-35. In addition, we observe that the frequency of the insertion shows a geospatial north-south gradient. The insertion is less common in northern (around 6-7% in patients and 2% in controls from Sweden and Poland) compared to southern Europeans (around 16% in patients and 8% in controls from Italy) and shows a stronger attributable risk (AR) in the southern European population. Our study provides evidence that the prevalence of achalasia may differ between populations.

  Study EGAS00001001515

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

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

  Study EGAS00001001690

• Ethiopia_Genome_Project__low_coverage_

  The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

  Study EGAS00001000238

• UK10K OBESITY TWINSUK

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These TwinsUK samples are a subset of the TwinsUK cohort samples with a BMI>=40 that will be exome sequenced by the Obesity group.

  Study EGAS00001000306

• The Prediction and Prevention of Preeclampsia

  The Prediction and Prevention of Preeclampsia (PREDO) is a prospective birth cohort study of Finnish women who were pregnant between 2005 and 2010 and their children. The PREDO study cohort was set up to identify novel risk factors and biomarkers in pregnant women associated with the development of preeclampsia and intrauterine growth restriction (IUGR). Two groups of pregnant women were enrolled during the first trimester of their pregnancy: first, pregnant women with a known clinical risk factor status for preeclampsia and IUGR (N=1079; 969 of these women had at least one and 110 had none of the known risk factors for preeclampsia and IUGR), and second, pregnant women who volunteered to participate regardless of their risk factor status for preeclampsia and IUGR (N=3698). Both samples were monitored during pregnancy with four clinical visits and bi-weekly self-reports. The post-delivery follow-up has taken place at approximately 2 weeks, 6 months, and 3.5 years after the delivery. The most recent follow-up started in 2016 and is ongoing. In the high-risk sample, cord blood samples were collected and genome-wide genotyping has been conducted with Illumina OmniExpressExome 1.2 and genome-wide methylation status has been analysed with Illumina 450k HM array.

  Study EGAS00001001898

• Convergent somatic mutations in effectors of insulin signalling in chronic liver disease

  Chronic liver disease is associated with metabolic dysregulation, liver failure and hepatocellular carcinoma. We analysed somatic mutations from 1202 genomes across 32 liver samples, including normal controls, alcohol-related and non-alcoholic fatty liver disease. Five of 27 patients with liver disease carried hotspot driver mutations in FOXO1, the major transcription factor downstream of insulin signalling. FOXO1 mutations were independently acquired by up to 5 distinct clones within the same patient’s sample, and impaired insulin-mediated nuclear export of FOXO1. GPAM, which produces storage triacylglycerol from dietary calories, also had significant excess of mutations, similarly exhibiting convergent evolution within biopsies. Telomeres were shorter in diseased than normal liver, with attrition more pronounced in larger clones. Multiple independent acquisitions of drivers within one small liver sample imply that such mutations could affect hundreds of grams of tissue across the whole organ, potentially contributing to systemic metabolic dysfunction.

  Dataset EGAD00001006255

• PSCP_bisulphite analysis in hESCs (2018-08-13)

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. 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-08-13.

  Dataset EGAD00001004295

• NIPT samples for systematic evaluation of NIPT aneuploidy detection software tools

  This dataset contains two sets of samples. The reference sample set consists of a total of 669 samples that had been reported previously to be euploid by the NIPTIFY screening test. The validation sample set is based on a previously published validation study by Zilina et al. (1), consisting of 423 samples, of which 259 were high-risk pregnancies that had undergone diagnostic invasive prenatal analysis (1). All samples were sequenced with Illumina NextSeq 500 platform, producing 85 bp single-end reads with an average per-sample coverage of 0.32× at the University of Tartu, Institute of Genomics Core Facility, according to the manufacturer’s standard protocols, as described previously (1). This study was performed with the approval of the Research Ethics Committee of the University of Tartu (#315/T-13). 1. Zilina O, Rekker K, Kaplinski L, Sauk M, Paluoja P, Teder H, et al. Creating basis for introducing non‐invasive prenatal testing in the Estonian public health setting. Prenat Diagn [Internet]. 2019 Dec 6;39(13):1262-8. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1002/pd.5578

  Dataset EGAD00001007712

• Mutational landscape of normal epithelial cells in Lynch Syndrome patients

  Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole-genome sequenced 152 crypts from normal and neoplastic epithelial tissues from LS patients. In normal tissues the repertoire of mutational processes and mutation rates were similar to those found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and mutational signatures consistent with MMR deficiency was identified, which may represent a very early stage of LS pathogenesis. Phylogenetic tress of tumour crypts indicated that the most recent ancestor cell of each tumour was already MMR deficient and had experienced multiple clonal evolution cycles. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

  Dataset EGAD00001008092

• FASTQ files of total RNA-Seq data from the POPS PET (pre-eclamptic) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003508

• RNA-seq analysis of transcriptome variation with human ESC subclones

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny 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/ .

  Dataset EGAD00001004266

• Genetics causes of male infertility in 185 patient-parent trios from Netherlands and UK

  The dataset is composed by the raw and processed sequencing data generated from 185 Patients affected by azoospermia or severe oligozoospermia recruited from the Netherlands and the UK.

  Dataset EGAD00001007862

• PSCP_bisulphite analysis in hESCs

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. 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/

  Dataset EGAD00001002235

• FASTQ files of total RNA-Seq data from the POPS SGA (Small for Gestational Age) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003507

• Swarm Learning to identify COVID-19, tuberculosis and leukemia patients based on blood transcriptomes

  Identification of patients with life-threatening diseases including leukemias or infections such as tuberculosis or COVID-19 is an important goal of modern precision medicine. However, there is an increasing divide between what is technically possible and what is allowed because of privacy legislation. We have recently illustrated that classical machine learning can identify leukemia patients based on their blood transcriptomes. To facilitate integration of any omics data from any data owner world-wide without violating privacy laws, we here introduce Swarm Learning (SL), a decentralized machine learning approach uniting edge computing, artificial intelligence (AI), blockchain and privacy protection without the need for a central coordinator thereby going beyond federated learning. To illustrate its feasibility, using more than 12,000 transcriptomes from peripheral blood mononuclear cells and more than 2,000 peripheral blood transcriptomes we demonstrate that SL of omics data distributed across different individual sites leads to disease classifiers that outperform those developed at individual sites. Yet, SL completely protects local privacy regulations by design. We propose this approach to noticeably accelerate the introduction of precision medicine.

  Dataset EGAD00001006231

• PSCP_mutation analysis in hESCs

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. This data is part of a pre-publication release. For information on the proper use of pre-publication data shred by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Dataset EGAD00001002231

• Whole Exome Sequencing PPGL

  Whole exome sequencing (WES) libraries were prepared with the SureSelectXT Human All Exon V6+COSMIC target enrichment system (Agilent, 5190-9307) by following kit manufacturer`s instructions. 200 ng of genomic DNA from FFPE tumoral samples and 1 µg from FF or germline samples were used as starting material. SureSelectXT HS (Agilent, G9704A) and SureSelectXT2 (Agilent, G9621A) reagent kits were alternatively used for FFPE derived and good quality samples, respectively. DNAs were fragmented on a Covaris S2 sharing instrument based on initial sample’s integrity. Fragmented DNA samples were processed through subsequent enzymatic treatments of end-repair, dA-tailing, and ligation to Illumina's adapters. Adapter-ligated libraries were PCR amplified with Illumina PE primers and subsequently hybridized to SureSelect Oligo Capture Library Mix. Biotinylated hybrids were captured and the enriched fraction eluted. Libraries were completed by limited-cycle PCR with Illumina TruSeq primers and KAPA HiFi HotStart DNA pol (Roche KK2501). Sample-specific libraries were combined into equimolarly balanced pools and applied to Illumina flow cells for cluster generation. Sequencing was performed on either the Illumina HiSeq2500 or NovaSeq6000 in a paired-end 100bp reads mode to a median target coverage of ~200x and ~100x for tumors and germline, respectively.

  Study EGAS00001006043

• The Spatial Heterogeneity in Multiple Myeloma - from the Subclonal Architecture to the Immune Microenvironment (partly hipo_K08K)

  Multiple myeloma (MM) cells show pronounced heterogeneity not only within a given patient but also between spatially separated tumor loci in the bone marrow. Understanding this spatial heterogeneity is emerging as a critical challenge for the successful treatment of the disease. Yet, our understanding of spatial differences in the subclonal architecture, molecular signatures and interactions with the tumor microenvironment remains very limited. To address this shortcoming, we performed bulk and single-cell multi-region sequencing, including random bone marrow samples from the iliac crest and paired imaging-guided focal lesion specimens from 15 newly diagnosed MM patients. We found a median of 5 subclones per patient and unique subclones in focal lesions. Central features of spatial heterogeneity included a consistent down-regulation of the chemoattractant cytokines CXCL7 and CXCL12 in focal lesion tumor cells as well as a significant depletion of macrophages in the focal lesion microenvironment. In contrast, a site-specific expansion of T-cell clones was not detected in focal lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity with potential implications for immunotherapies and molecular studies analyzing the role of gene signatures and MM-microenvironment interactions in disease progression.

  Study EGAS00001006090

• Cardiovascular disease biomarkers derived from circulating cell-free DNA methylation

  Acute coronary syndrome (ACS) remains a major cause of worldwide mortality. The syndrome occurs when blood flow to the heart muscle is decreased or blocked, causing muscle tissues to die or malfunction. There are three main types of ACS: Non-ST-elevation myocardial infarction, ST-elevation myocardial infarction, and unstable angina. The treatment depends on the type of ACS, and this is decided by a combination of clinical findings, such as electrocardiogram and plasma biomarkers. Circulating cell-free DNA (ccfDNA) is proposed as an additional marker for ACS since the damaged tissues can release DNA to the bloodstream. We used ccfDNA methylation profiles for differentiating between the ACS types and provided computational tools to repeat similar analysis for other diseases. We leveraged cell type specificity of DNA methylation to deconvolute the ccfDNA cell types of origin and to find methylation-based biomarkers that stratify patients. We identified hundreds of methylation markers associated with ACS types and validated them in an independent cohort. Many such markers were associated with genes involved in cardiovascular conditions and inflammation. ccfDNA methylation showed promise as a non-invasive diagnostic for acute coronary events. These methods are not limited to acute events, and may be used for chronic cardiovascular diseases as well.

  Study EGAS00001007263

• Human liver NPCs single cell project

  Independent of their inflammatory phenotype, macrophages are key orchestrators of hepatic metabolism. Non-alcoholic fatty liver disease (NAFLD) often occurs in obese individuals and is among the most common causes of cirrhosis, the terminal chronic liver disease that may necessitate liver transplantation. While multiple populations of macrophages have been described in the human liver, their function and turnover in obese patients at high risk of developing NAFLD and cirrhosis is currently unknown. Herein we identified a specific human population of resident liver myeloid cells that protects against the metabolic impairment associated with obesity. By studying the turnover of liver myeloid cells in individuals undergoing liver transplantation using markers of donor-recipient mismatch, we made the novel discovery that liver myeloid cell turnover differs between humans and mice. Using single cell techniques and flow cytometry we determined that the proportion of the protective resident liver myeloid cells, denoted liver myeloid cells 2 (LM2), decreases during obesity. Functional validation approaches using human 2D and 3D cultures revealed that the presence of LM2 ameliorates the oxidative stress associated with obese conditions. Our study indicates that resident myeloid cells could be a therapeutic target to decrease the oxidative stress associated with NAFLD.

  Study EGAS00001007194

• Comparison of structural variations from 10X Genomics linked-reads and conventional Illumina short-reads sequencing

  Structural variations (SVs) are large genomic rearrangements that can drive many diseases. Conventional short-reads whole genome sequencing (cWGS) allows their identification with base-pair resolution, but suffers from high false discovery rate. cWGS taps in short-range information from short-reads while linked-reads sequencing (10XWGS) utilizes long-range information. 10XWGS allows linkage of short-reads originating from the same large DNA molecule with a unique barcode captured in a gel bead in emulsion. This mitigates alignment-based artefacts from cWGS especially in repetitive regions. However, the false discovery rate of this technology is unclear. In this study, we performed a comprehensive analysis of different type and size of SVs predicted from these two technologies. The SVs common between both technologies were found to be highly specific by PCR and Sanger sequencing while validation rate dropped for uncommon events. Further, we propose a novel enrichment approach for filtering out false positive calls from both the technologies independently. To this end, we trained a machine learning model for respective technologies and used it to characterise SVs from MCF7 cell line and a primary breast cancer tumor with high precision. This approach would be valuable in understanding true mechanisms driven by SVs in various diseases.

  Study EGAS00001004093

• Pancreatic, Small-intestinal and Pulmonary Neuroendocrine Tumors

  Purpose: The primary origin of neuroendocrine tumor metastases can be difficult to determine by histopathology alone, but is critical for therapeutic decision making. DNA methylation-based profiling is now routinely used in the diagnostic workup of brain tumors. This has been enabled by the availability of cost-efficient array-based platforms. We have extended these efforts to augment histopathological diagnosis in neuroendocrine tumors. Experimental Design and Results: We compiled data of 69 small-intestinal, pulmonary, and pancreatic neuroendocrine tumors. These data were used to build a ridge regression calibrated random forest classification algorithm (NEN-ID) that predicts the origin of tumor samples with high accuracy (> 95%). The model was validated during 3x3 nested cross validation and tested in a local (n=26) and external (n=172) cohort. In addition, we show that our diagnostic approach is robust across a range of possible confounding experimental parameters such as tumor purity and array quality. A software infrastructure and online user interface was built to make the model available to the scientific community. Conclusions: This DNA methylation-based prediction model can be used in the workup for patients with neuroendocrine tumors of unknown primary. To facilitate validation and clinical implementation, we provide a user-friendly, publicly available web-based version of NEN-ID.

  Study EGAS00001004878

• Sclerosing epithelioid fibrosarcoma case report

  Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft tissue sarcomathought to originate in fibroblasts of the tissues surrounding tendons, ligaments and muscles. Minimally responsive to conventional cytotoxic chemotherapies, greater than 50% of SEF patients experience relapse and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-year-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. The patient described in this report experienced persistent disease progression despite being heavily treated with multiple surgeries, radiation, numerous chemotherapies and targeted therapeutics. Primary tumor, metastatic and relapse sites were sequenced by whole genome, whole exome, and deep transcriptome next generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and a discussion as to potential therapeutic strategies for SEF.

  Study EGAS00001005214

• Adult-type Granulosa Cell Tumour of the Ovary

  A recurrent mutation in FOXL2 (c.402C>G; p.C134W) is present in over 95% of adult-type granulosa cell tumours (AGCTs). In contrast, various loss-of-function mutations in FOXL2 lead to the development of blepharophimosis, ptosis and epicanthus inversus syndrome (BPES). BPES is characterized by an eyelid malformation often accompanied with primary ovarian insufficiency but not granulosa cell tumours. Two recent studies suggest that FOXL2 C402G is a gain or change-of-function mutation with altered DNA-binding specificity. Another study proposes that FOXL2 C402G is selectively targeted for degradation inducing somatic haploinsufficiency suggesting its role as a tumour suppressor. The latter study relies on data indicative of an FOXL2 allelic imbalance in AGCTs. Here we present RNA-seq data as genetic evidence that no real allelic imbalance is observed at the transcriptomic level in AGCTs. Additionally, there is no loss of protein expression in tumours harboring the mutated allele. These data and other features of this mutation as compared to other oncogenes and tumour suppressor genes argue strongly against FOXL2 being a tumour suppressor in this context. Given the likelihood that FOXL2 C402G is oncogenic, targeting the variant protein or its downstream consequences is the most viable path forward to identifying an effective treatment for this cancer.

  Study EGAS00001005414

• The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma

  Patients with metastatic uveal melanoma (UM) have an abysmal prognosis. Preclinical studies have suggested that epigenetic therapy could enhance immunogenicity of cancer cells. Here we test if epigenetic therapy would enhance PD-1 immunotherapy in patients with metastatic UM. We report the results of the PEMDAC phase 2 clinical trial (n=29; NCT02697630) where the HDAC inhibitor entinostat was combined with the PD-1 inhibitor pembrolizumab in patients with metastatic UM. The primary endpoint was objective response rate (ORR), and was met with an ORR of 14%. The clinical benefit rate at 18 weeks was 28%, median progression free survival was 2.1 months and the median overall survival was 13.4 months. Toxicities were manageable, and there were no treatment-related deaths. Extensive genomics studies were performed using DNA/RNA and single cell sequencing and flow cytometry. Objective responses and prolonged survival were seen in three patients with BAP1 wildtype tumors, and in one patient with an iris melanoma that exhibited a UV signature. Longer survival also correlated with low baseline ctDNA levels or LDH. In conclusion, HDAC inhibition and anti-PD1 immunotherapy results in durable responses in a subset of patients with metastatic UM. Further exploration of combined immunotherapy and epigenetic therapy in metastatic UM is warranted.

  Study EGAS00001005478

• Unifying recovery dynamics in heterogeneous diseases exemplified by COVID-19

  Studying the recovery dynamics in different diseases poses many difficulties, as patients often display high heterogeneity in their recovery courses. Moreover, most attempts to study disease dynamics focus on disease progression rather than disease recovery mechanisms. To model recovery dynamics, using severe COVID-19 as the example, we align heterogeneous recovery trajectories via a novel computational scheme applied to longitudinally sampled blood transcriptomes. We thus generate pseudotime trajectories, which we then link to cellular and molecular mechanisms based on cell deconvolution analysis and molecular pathway prediction, thus presenting a unique framework for studying recovery processes over time. Specifically, mature neutrophils displayed a gradual decrease during recovery, allowing superior useability for outcome prediction compared to currently used clinical markers. Further, we discovered a recovery-related regulatory change in gene programs resulting in immune rebalancing between interferon and NFkB activity and the restoration of cell homeostasis. We thus propose regulatory mechanisms governing COVID-19 recovery and suggest mature neutrophils and additional gene markers, as novel clinical biomarkers for disease outcome. Overall, we present a novel clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

  Study EGAS00001005735

• Pancreatic tropism of metastatic renal cell carcinoma

  Renal cell carcinoma (RCC) is characterized by a particularly broad metastatic swath, and enigmatically, when the pancreas is a destination, the disease is associated with improved survival. Intrigued by this observation, we sought to characterize the clinical behavior, therapeutic implications, and underlying biology. While pancreatic metastases (PM) are infrequent, we identified 31 patients across two institutional cohorts and show that improved survival is independent of established prognostic variables, that these tumors are exquisitely sensitive to anti-angiogenic agents and resistant to immune checkpoint inhibitors (ICI), and that they are characterized by a distinctive biology. Primary tumors of patients with PM exhibited frequent PBRM1 mutations, 3p loss, and 5q amplification, along with a lower frequency of aggressive features such as BAP1 mutations and loss of 9p, 14q, and 4q. Gene expression analyses revealed constrained evolution with remarkable uniformity, reduced T-effector gene signatures and increased angiogenesis. Similar findings were observed histopathologically. Thus, RCC metastatic to the pancreas is characterized by indolent biology, heightened angiogenesis, and an uninflamed stroma, likely underlying its good prognosis, sensitivity to anti-angiogenic therapies, and refractoriness to ICI. These data suggest that metastatic organotropism may be an indicator of a particular biology with prognostic and treatment implications for patients.

  Study EGAS00001004208

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

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

  Study EGAS00001002211

• Host pathogen interaction long read transcriptome

  Our incomplete knowledge of the human transcriptome impairs the detection of disease-causing variants, in particular in transcripts only expressed under certain conditions. These transcripts are often lacking from reference transcript sets, such as Ensembl/GENCODE and RefSeq, and could be relevant for establishing genetic diagnoses. We present SUsPECT (Solving Unsolved Patient Exomes/gEnomes using Custom Transcriptomes), a pipeline based on the Ensembl Variant Effect Predictor (VEP) to predict variant impact on custom transcript sets, such as those generated by long-read RNA-sequencing, for downstream prioritization. Our pipeline predicts the functional consequence and likely deleteriousness scores for missense variants in the context of novel open reading frames predicted from any transcriptome. We demonstrate the utility of SUsPECT by uncovering potential mutational mechanisms of pathogenic variants in ClinVar that are predicted to be benign using the reference transcript annotation. In further support of SUsPECT’s utility, we identified an enrichment of immune-related variants predicted to have a more severe molecular consequence when annotating with a newly generated transcriptome from stimulated immune cells instead of the reference transcriptome. Our pipeline outputs crucial information for further prioritization of potentially disease-causing variants for any disease and will become increasingly useful as more long-read RNA sequencing datasets become available.

  Study EGAS00001006779

• Direct Detection of Early-stage Cancers Using Circulating Tumor DNA

  Early detection and intervention are likely to be the most effective means for reducing morbidity and mortality of human cancer. However, development of methods for noninvasive detection of early-stage tumors has remained a challenge. We have developed an approach called targeted error correction sequencing (TEC-Seq) that allows ultrasensitive direct evaluation of sequence changes in circulating cell-free DNA using massively parallel sequencing. We have used this approach to examine 58 cancer-related genes encompassing 81 kb. Analysis of plasma from 44 healthy individuals identified genomic changes related to clonal hematopoiesis in 16% of asymptomatic individuals but no alterations in driver genes related to solid cancers. Evaluation of 200 patients with colorectal, breast, lung, or ovarian cancer detected somatic mutations in the plasma of 71, 59, 59, and 68%, respectively, of patients with stage I or II disease. Analyses of mutations in the circulation revealed high concordance with alterations in the tumors of these patients. In patients with resectable colorectal cancers, higher amounts of preoperative circulating tumor DNA were associated with disease recurrence and decreased overall survival. These analyses provide a broadly applicable approach for noninvasive detection of early-stage tumors that may be useful for screening and management of patients with cancer.

  Study EGAS00001002577

• Spatiotemporal Genomic Profiling of Intestinal Metaplasia Reveals Clonal Dynamics of Gastric Cancer Progression

  Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. We analyzed 1256 gastric samples (1152 IMs) from 692 subjects from a 10-year prospective study. We identified 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9), largely occurring as subclonal events. Analysis of clonal dynamics between and within subjects, and also longitudinally across time, revealed that IM clones are likely transient but increase in size upon progression to dysplasia, with eventual transmission of genetic events to paired GCs. Single-cell and spatial profiling highlighted changes in tissue ecology and lineage heterogeneity in IM, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling revealed expression-based molecular subtypes of IM, including a body-resident “pseudoantralized” subtype associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical-genomic models outperform clinical-only models in predicting IMs likely to progress. Our results raise opportunities for GC precision prevention and interception by highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression.

  Study EGAS00001007067

• Long-term organoid culture of a small intestinal neuroendocrine tumor

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007093

• High Grade Serous Ovarian Carcinomas Originate in the Fallopian Tube

  High-grade serous ovarian carcinoma (HGSOC) is the most frequent type of ovarian cancer and has a poor outcome. It has been proposed that fallopian tube cancers may be precursors of HGSOC but evolutionary evidence for this hypothesis has been limited. To provide insights into the origins of HGSOC, we We performed whole-exome sequence and copy number analyses of laser-capture microdissected fallopian tube lesions (p53 signatures, serous tubal intraepithelial carcinomas (STICs), and fallopian tube carcinomas), ovarian cancers, and metastases from nine patients with HGSOC. The majority of tumor-specific sequence and structural alterations in ovarian cancers were present in STICs, including those affecting TP53, BRCA1, BRCA2 or PTEN genes. An eEvolutionary analyseis revealed that p53 signatures and STICs were are the precursors of ovarian carcinoma which in turn gave rise to metastatic lesionsand identify a. In one patient we identified a second STIC as a metastasis in the fallopian tube opposite from the affected ovary. These analyses revealed a window of seven years between the development of a STIC and the initiation of ovarian carcinomas, with development of metastases following rapidly thereafter. Our results provide insights into the etiology of ovarian cancer and have implications for the prevention, early detection and therapeutic intervention of this disease.

  Study EGAS00001002589

• Mapping_gene_environment_interactions_in_macrophages

  Immune cells sense and respond to external stimuli such as pathogens and initiate an inflammatory response. Genetic variants altering the behaviour and responses of immune cells may inform the results of disease association studies and provide an insight into the regulation of immune response. Studies mapping quantitative trait loci (QTL) are usually performed in quiescent or naïve cells and therefore any variants with a condition-specific manner (e.g. those with an effect only detectable in a specific time-point after cell stimulation) are likely to remain undiscovered. So far, only a handful of studies have attempted to study QTL in multiple states and even on those studies the number of states has remained relatively low (<5). Only very recently Horst and colleagues published the largest study to date investigating the effect of host variable (e.g. age, sex, contraceptive use etc.) on protein levels in macrophages stimulated with 19 different conditions (Horst et al, 2016). This project aims to develop a high-throughput experimental platform to detect QTL in a large number of states using iPSC derived macrophages. 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 EGAS00001002268

• Single cell transcriptomic landscape of pediatric B-cell acute lymphoblastic leukemia: dissection of transcriptional heterogeneity and B-cell developmental state

  Introduction: Genomic analyses of bulk ALL samples have improved our understanding of the genetic basis and risk stratification of B-ALL, but do not directly examine intratumor heterogeneity or enable inference of leukemia developmental state and cell of origin. Methods: We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) (10X Genomics 5’v2) and compared them to a scRNA-seq reference map of normal human B-cell development. Results: Analysis of heterogeneity of inferred DNA copy number alterations at single cell level showed that aneuploid ALL with near haploidy and hyperdiploid harbor chromosomal losses or gains in all blasts, consistent with early, synchronous origin rather than sequential accumulation. Intra-sample heterogeneity was driven by five gene expression signatures: cell cycle, metabolism, differentiation, inflammation, refining leukemic subtyping. An atlas of B cell development was established to map B-ALL samples to specific developmental states. Hyperdiploid samples resembled pro-B cells, while MEF2D-r and TCF3::PBX1 showed high pre-B abundance. Patients with high early lymphoid abundance had poorer outcomes. Conclusions: Understanding variation in transcriptional programs and developmental states of B-ALL blasts by sc transcriptome refines existing clinical and genomic classifications and provides novel prognostic markers.

  Study EGAS00001007512

• Ewing and Ewing-like-sarcoma tumoroid biobank reveals distinct drug sensitivities: translocation makes the difference

  Ewing sarcoma (ES) and Ewing-like sarcomas (ELS) are highly aggressive tumors found in soft tissues and bone of mostly children and young adults. Despite being different in many aspects, including genetics, possible cell-of-origin, and pathology, patients with any of these entities all receive the same therapeutic regimen. Currently, there is a lack of preclinical cell models that can predict drug response of these entities in vitro. Here, we describe a pediatric ES and ELS tumor organoid (tumoroid) biobank containing tumoroids with different translocations, including EWSR1::FLI1, EWSR1::ERG, CIC::DUX4, and BCOR-rearrangements. Using histology, whole genome sequencing and RNA sequencing, we demonstrate that these tumoroids retain key properties of their matching patient tumors. In addition, we traced tumor subclones in the tumoroids across patient-matched longitudinal samples, which showed that cellular heterogeneity was maintained. We subsequently performed drug screening on the tumoroid models, which unveiled entity-specific drug sensitivity to various cytotoxic compounds and targeted compounds, including MCL-1 inhibitors for CIC::DUX4 sarcomas. Taken together, this newly established ES and ELS patient-derived tumoroid biobank generates a unique source of material for future basic cancer research and drug-screening.

  Study EGAS00001007941

• Integrated genomic analysis for HCC

  Liver cancer is a common malignancy and is the third cause of cancer-related death in the world and hepatocellular carcinoma (HCC) is the major type of liver cancer, accounting for ~90% of cases. Only 20 ~ 30% of patients can benefit from combined immunotherapy and targeted therapy due to the high inter-tumor heterogeneity of HCC. This work aims to identify homogeneous and robust molecular subtypes in HCC based on a large, homogenous and well-annotated cohort. We performed whole genome/exome sequencing (WGS/WES) and/or RNA sequencing (RNA-seq) on 529 HCC from 461 patients collected mainly in France. Based on the great consistence between genomic and transcriptomic data, we identified 9 robust HCC subtypes mainly based on driver mutations. And we further characterized HCC subtypes using genomic, transcriptomic and clinicopathological features. Among 9 subtypes, 5 subtypes belonged to chromosome instable tumors, while 3 subtypes belonged to chromosome stable tumors. In addition, our subtypes were associated with prognosis and we also found different distributions of subtypes across various features, such as etiology and gender. Our study provided a robust molecular classification based on a large HCC dataset from Western country, it improves our understanding of the mechanisms of carcinogenesis and will facilitate development of genome-based precision medicine in HCC.

  Study EGAS00001007957

• Detection of cancers three years prior to diagnosis using plasma cell-free DNA

  A major goal of early cancer detection is to identify subclinical disease when the tumor burden is low, so that treatments are more effective. But how early can cancers be detected prior to clinical signs or symptoms? This question can be answered only through the evaluation of participants whose clinical course has not been altered by the study itself. We here describe such an evaluation, performed on prospectively collected plasma samples from the Atherosclerosis Risk in Communities (ARIC) study, including 26 participants diagnosed with cancer and 26 matched controls. At the index time point, eight of these 52 participants scored positively with a multicancer early detection (MCED) blood test. All eight of these participants were diagnosed with cancer within 4 months after blood collection. In six of these 8 participants, we were able to assess an earlier plasma sample collected 3.1 to 3.5 years prior to clinical diagnosis. In four of these six participants, the same mutations detected by the MCED test could be identified, but at 8.6 to 79-fold lower levels. These results demonstrate that it is possible to detect circulating tumor DNA (ctDNA) more than three years prior to clinical diagnosis, and provide benchmark sensitivities required for the success of ctDNA-based tests for this purpose.

  Study EGAS00001008068

• Tumor-derived exosomes modulate PD-L1 expression in monocytes

  In chronic lymphocytic leukemia (CLL), monocytes and macrophages are skewed toward protumorigenic phenotypes, including the release of tumor-supportive cytokines and the expression of immunosuppressive molecules such as programmed cell death 1 ligand 1 (PD-L1). To understand the mechanism driving protumorigenic skewing in CLL, we evaluated the role of tumor cell–derived exosomes in the cross-talk with monocytes. We carried out RNA sequencing and proteome analyses of CLL-derived exosomes and identified noncoding Y RNA hY4 as a highly abundant RNA species that is enriched in exosomes from plasma of CLL patients compared with healthy donor samples. Transfer of CLL-derived exosomes or hY4 alone to monocytes resulted in key CLL-associated phenotypes, including the release of cytokines, such as C-C motif chemokine ligand 2 (CCL2), CCL4, and interleukin-6, and the expression of PD-L1. These responses were abolished in Toll-like receptor 7 (TLR7)–deficient monocytes, suggesting exosomal hY4 as a driver of TLR7 signaling. Pharmacologic inhibition of endosomal TLR signaling resulted in a substantially reduced activation of monocytes in vitro and attenuated CLL development in vivo. Our results indicate that exosome-mediated transfer of noncoding RNAs to monocytes contributes to cancer-related inflammation and concurrent immune escape via PD-L1 expression.

  Study EGAS00001002377

• Pooled scRNA-seq of iPSC-derived neural stem cells from ADHD and control individuals

  This dataset comprises paired-end FASTQ files from single-cell RNA sequencing (scRNA-seq) of neural stem cells (NSCs) derived from induced pluripotent stem cells (iPSCs). The iPSCs were reprogrammed from somatic cells (keratinocytes, PBMCs, or lymphoblastoid cell lines) collected from individuals diagnosed with attention-deficit hyperactivity disorder (ADHD) and age-matched healthy controls (ages 6–18 years). Reprogramming was carried out using the non-integrative Sendai virus system (CytoTune-iPS 2.0, Thermo Fisher Scientific), introducing OCT3/4, SOX2, KLF4, and cMYC. iPSCs were cultured in Essential 8 medium on vitronectin and underwent rigorous quality control including immunocytochemistry, qRT-PCR, mycoplasma testing, and SNP genotyping (Illumina Global Screening Array). NSCs were generated using a standardized neural induction protocol and validated via qRT-PCR and immunostaining for markers including PAX6, SOX2, NESTIN, TUJ1, and FOXG1. Two independent scRNA-seq multiplexed experiments were performed: Experiment 1 (ADHD: n = 7; Control: n = 6) and Experiment 2 (ADHD: n = 9; Control: n = 7). Samples were pooled and demultiplexed using the Ensemblex algorithm, incorporating prior genotype data. After quality control filtering (>500 unique RNA molecules per cell), 30,433 single-cell transcriptomes were retained. Differential gene expression analysis was conducted to identify transcriptional signatures associated with ADHD compared to controls.

  Study EGAS00001008169

• Patient-derived pediatric brain tumor organoids faithfully recapitulate primary tumors

  Extensive molecular analysis by many groups has revealed a heterogeneous landscape of embryonal brain tumors, including medulloblastoma (MB) and atypical teratoid/rhabdoid tumor (ATRT). MB, the most common malignant pediatric brain tumor, is comprised of four molecularly distinct subgroups: Wingless (WNT), Sonic Hedgehog (SHH), Group 3 (G3), and Group 4 (G4). Among these, MB-SHH harboring MYCN amplification and TP53 mutation and MB-G3 with MYC overexpression or amplification constitute the most aggressive tumors, have the poorest clinical outcomes, and lack effective therapeutic options for recurrent disease. Here we describe the establishment of tumor organoid (TO) models from patient-derived orthotopic xenografts (PDOX) of MYCN-amplified and TP53-mutant MB-SHH, MYC-amplified MB-G3 and two ATRT which all grew in mice when reimplanted into the cortices of naïve recipient animals (TOX). Methylation, bulk and single-cell RNA sequencing, and whole genome sequencing demonstrated that TOs and TOXs faithfully recapitulated the molecular landscape and heterogeneity of PDOXs and primary patient tumors compared to established cell lines. The development of TOs will facilitate in vitro high-throughput drug or CRISPR screens without the need for use tumor-bearing mice. This will also accelerate the identification and validation of new vulnerabilities and in vivo preclinical testing toward clinical trials.

  Study EGAS00001008305

• Botswana 15 autosomal unlinked microsatellites

  Population structure was investigated in 990 Botswana individuals according to ethno-linguistics, Bantu and Khoisan, and geography (the nine administrative districts) using the AmpFLSTR Identifiler autosomal microsatellite markers. Genetic diversity and forensic parameters were calculated for the overall population, and according to ethno-linguistics and geography. The overall combined power of exclusion (CPE) was 0.9999965412 and the combined match probability 6,28 x 10-19. CPE was highest for the Khoisan Tuu ethnolinguistic group and the Northeast District at 0.9999582029 and 0.9999922652 respectively. CMP ranged from 6.28 x 10-19 (Khoisan Tuu) to 1,02 x 10-18 (Northwest district). Using pairwise genetic distances (FST), analysis of molecular variance (AMOVA), factorial correspondence analysis (FCA), and the unsupervised Bayesian clustering method found in STRUCTURE and landscape genetics software TESS, ethno-linguistics were found to have a greater influence on population structure than geography. FCA showed clustering between Bantu and Khoisan, and within the Bantu. This Bantu sub-structuring was not seen with STRUCTURE and TESS, which detected clustering only between Bantu and Khoisan. The patterns of population structure revealed highlight the need for regional reference databases that include ethno-linguistic and geographic location information. These markers have important potential for bio-anthropological studies as well as for forensic applications.

  Study EGAS00001002380

• Inhibiton of the GABPB1L-containing GABP tetramer is sufficient to reverse replicative immortality in TERT promoter mutant glioblastoma cells.

  Over 50 types of cancer acquire TERT promoter mutations. These single point mutations reactivate telomerase, allowing for indefinite maintenance of telomere length and enabling cellular immortalization. The transcription factor binding site created by the point mutations specifically recruit the ETS factor GABP, a multimeric transcription factor composed of the GABPα and GABPβ subunits. GABP can form two functionally independent transcription factor species – a dimer or a tetramer – depending on which of the structurally distinct GABPβ isoforms is incorporated into the complex. In this study, we show that genetic disruption of GABPβ1L, a tetramer forming isoform of GABPβ that is dispensable in normal development, results in TERT silencing in a TERT promoter mutation dependent manner. Failure to activate TERT expression by GABPβ1L culminates in telomere dysfunction, DNA damage, and mitotic cell death exclusively in TERT promoter mutant cells. Furthermore, exogenous expression of TERT is sufficient to prevent telomere degradation and loss of cell viability in GABPβ1L-reduced lines bearing the mutant TERT promoter. Orthotopic injection of tetramer-deficient mutant TERT promoter GBM cells rendered lower tumor burden and prolonged the overall survival of the tumor-bearing mice. These results highlight the potentially widespread role of GABPB1L in enabling replicative immortality of TERT promoter cancers.

  Study EGAS00001002582

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC -- CA209-009

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004290

• Allele Balance Bias Identifies Systematic Genotyping Errors and False Disease Associations

  In recent years, Next Generation Sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the allele balance distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies (RVAS).

  Study EGAS00001003027

• Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

  Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that forms intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to show treatment resistance and conserve the tumour bulk. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogenic tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. SNP arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors, consistent with the heterogenic molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

  Study EGAS00001003438