4051 results for "*geneti*"

in 78.80 milliseconds.

• Targeting AXL Kinase Uniquely Sensitizes Therapy-Insensitive Leukemic Stem and Progenitor Cells to Venetoclax Treatment in Acute Myeloid Leukemia

  The abundance of genetic abnormalities and phenotypic heterogeneities in AML pose significant challenges to developing improved treatments. Here we demonstrated that a key GAS6/AXL axis is highly activated in AML patient cells, particularly in leukemic stem cells. We developed a potent, selective AXL inhibitor that has excellent pharmaceutical properties and efficacy against preclinical patient-derived xenotransplantation models of AML. Importantly, inhibition of AXL sensitized AML stem/progenitor cells to venetoclax treatment, with strong synergistic effects in vitro and in vivo. Mechanistically, single-cell RNA-sequencing and functional validation studies uncovered that AXL inhibition or in combination with venetoclax potentially targets intrinsic metabolic vulnerabilities of AML stem/progenitor cells, which shows a distinct transcriptomic profile and inhibits mitochondrial oxidative phosphorylation. These findings have direct translational impact on the treatment of AML and other cancers with high AXL activity.

  Study EGAS00001004663

• International consensus definition of DNA methylation subgroups in juvenile myelomonocytic leukemia

  Known clinical and genetic markers have limitations in predicting disease course and outcome in juvenile myelomonocytic leukemia (JMML). DNA methylation (DNAme) patterns in JMML have correlated with outcome across multiple studies, suggesting it as a biomarker to improve patient stratification. However, standardized approaches to classify JMML based on DNAme patterns are lacking. We therefore sought to define an international consensus for DNAme subgroups in JMML and develop classification methods for clinical implementation.

  Study EGAS00001004682

• Allele-specific expression of GATA2 due to epigenetic dysregulation in double mutated CEBPA AML

  Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid or cancer-related gene. GATA2 ASE was strongly associated with CEBPA double mutations (CEBPA DM), with 95% of cases presenting GATA2 ASE. In CEBPA DM AML with GATA2 mutations, the mutated allele was preferentially expressed. We found that GATA2 ASE is a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPA DM AML. Acquisition of GATA2 ASE involved silencing of one allele via promoter methylation, compensated by overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission together with GATA2 ASE. In summary, we propose that GATA2 ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPA DM. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

  Study EGAS00001004684

• Chromosomal copy number heterogeneity predicts survival rates across cancers

  Survival rates of cancer patients vary widely within and between malignancies. While genetic aberrations are at the root of all cancers, individual genomic features cannot explain these distinct disease outcomes. In contrast, intra-tumour heterogeneity (ITH) has the potential to elucidate pan-cancer survival rates and the biology that drives cancer prognosis1,2. Unfortunately, a comprehensive and effective framework to measure ITH across cancers is missing3. Here, we introduce a scalable measure of chromosomal copy number heterogeneity that predicts patient survival across cancers. We show that the level of ITH can be derived from a single-sample copy number profile. Using gene-expression data we demonstrate that ongoing chromosomal instability underlies the observed heterogeneity. Analysing 11,534 primary cancer samples from 35 different malignancies, we find that copy number heterogeneity can be accurately deduced and predicts cancer survival across tissues of origin and stages of disease. Our results provide a unifying molecular explanation for the different survival rates observed between cancer types.

  Study EGAS00001004702

• Whole genome sequencing of EBV Associated Nasopharyngeal Carcinoma

  Interplay between EBV infection and acquired genetic alterations during nasopharyngeal carcinoma (NPC) development remains vague. Here we report a comprehensive genomic analysis of 70 NPCs, combining whole-genome sequencing (WGS) of microdissected tumor cells with EBV oncogene expression to reveal multiple aspects of cellular-viral co-operation in tumorigenesis. Genomic aberrations along with EBV-encoded LMP1 expression underpin constitutive NF-κB activation in 90% of NPCs. A similar spectrum of somatic aberrations and viral gene expression undermine innate immunity in 79% of cases and adaptive immunity in 47% of cases; mechanisms by which NPC may evade immune surveillance despite its pro inflammatory phenotype. Additionally, genomic changes impairing TGFBR2 promote oncogenesis and stabilize EBV infection in tumor cells. Fine-mapping of CDKN2A/CDKN2B deletion breakpoints reveals homozygous MTAP deletions in 32-34% of NPCs that confer marked sensitivity to MAT2A inhibition. Our work concludes that NPC is a homogeneously NF-κB-driven and immune-protected, yet potentially druggable, cancer.

  Study EGAS00001004705

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• STAT5 is a therapeutically targetable vulnerability in cutaneous T-cell lymphoma

  Cutaneous T-cell lymphomas represent a heterogeneous group of lymphoproliferative disorders characterized by the infiltration and expansion of mature malignant T-cells into the skin. The disease encompasses distinct variants such as mycosis fungoides a skin resident variant, in which the malignant cells are restricted to skin lesions, and leukemic variant sézary syndrome characterized by circulating malignant T-cells. Here we integrated shallow whole genome sequencing and whole exome sequencing, on samples from four sézary syndrome patients and two advanced mycosis fungoides patients with blood involvement, to characterize genetic alterations and assess their therapeutic actionability. We identified an elevated copy number of 17q, resulting in increased STAT5 expression and activation, as a key factor in the pathophysiology of the disease.

  Study EGAS00001004719

• Mutational Landscape and Tumor Burden Assessed by Cell-Free DNA in Diffuse Large B-Cell Lymphoma: a Population-based Study

  Purpose: We analyzed the utility of cfDNA in a prospective population-based cohort to determine the mutational profile, assess tumor burden, and estimate its impact in response rate and outcome in patients with DLBCL. Experimental Design: One-hundred patients were diagnosed with DLBCL during the study period. Mutational status of 112 genes was studied in cfDNA by targeted next-generation sequencing. Paired formalin-fixed paraffin-embedded samples and volumetric PET/CT were assessed when available. Results: Appropriate cfDNA to perform the analyses was obtained in 79/100 cases. At least one mutation could be detected in 69/79 cases (87%). The sensitivity of cfDNA to detect the mutations was 68% (95% CI: 56.2-78.7). The mutational landscape found in cfDNA samples was highly consistent with that shown in the tissue and allowed genetic classification in 43% of the cases. A higher amount of ctDNA significantly correlated with clinical parameters related to tumor burden (elevated LDH and β2-microglobulin serum levels, advanced stage, and high-risk IPI) and total metabolic tumor volume assessed by PET/CT. In patients treated with curative intent, high ctDNA levels (>2.5 log hGE/mL) were associated with lower complete response (65% vs. 96%, P<0.004), shorter progression-free survival (65% vs. 85%, P=0.038) and overall survival (73% vs. 100%, P=0.007) at 2 years, although it did not maintain prognostic value in multivariate analyses. Conclusions: In a population-based prospective DLBCL series, cfDNA resulted an alternative source to estimate tumor burden and to determine the tumor mutational profile and genetic classification, which have prognostic implications and may contribute to a future tailored

  Study EGAS00001004733

• Prognostic and therapeutic significance of leukemia subtypes and minimal residual disease measurements in pediatric acute lymphoblastic leukemia treated with contemporary risk-directed trial: a cohort study

  Background Genomic analyses to identify leukemia subtypes and refinement of minimal residual disease (MRD) determination are two recent advances that would have substantial impact on the management of children with acute lymphoblastic leukemia (ALL). However, the combined application of these two approaches to determine the prognostic and therapeutic significance of MRD measurement and full spectrum of leukemia subtypes in the context of contemporary risk-directed therapy has not been studied comprehensively.MethodsBetween October 29, 2007 and March 26, 2017, 598 consecutive children with newly diagnosed ALL were enrolled in a clinical trial at St. Jude Children’s Research Hospital. Patients were provisional classified according to the presenting clinical and laboratory features including immunophenotype and cytogenetics. Patients with ETV6-RUNX1 and high hyperdiploid B-progenitor ALL were provisionally included in the low-risk groups; those with TCF3-PBX1, hypodiploidy and T-cell ALL in the standard (intermediate)-risk group; and patients with BCR-ABL1, KMT2A and early T-cell precursor ALL in the high-risk group. Final risk assignment was based primarily on MRD levels measured by flow cytometry on Day 15 (≥1%) and Day 42 (<0.01% or ≥1%) of remission induction. Using genomic classification based on cytogenetics and genomic analysis (single nucleotide polymorphism array analysis and transcriptome, exome, and genome sequencing), patients were classified into 17 subtypes. The primary aim was to determine the event-free survival, overall survival and cumulative risk of relapse according to specific leukemia subtypes and MRD levels on Days 8, 15 and 42 of remission induction. This clinical trial was registered at ClinicalTrials.gov, number NCT00549848.FindingsPatients with ETV6-RUNX1, high hyperdiploidy and DUX4-rearranged B-ALL all received low- or standard-risk therapy, and had the highest 5-year event-free survival rates: 98.4% (95% CI, 95.9-100), 95.3% (91.2-99.4), and 95.0% (84.2-100), respectively. All 142 patients with Day 8 MRD<0.01%, except two with KMT2A-rearranged ALL and one with TCF3-PBX1, were in continuous complete remission. The TCF3-PBX1, PAX5alt, T-cell, early T-cell precursor, iAMP21 and hypodiploid ALL had intermediate outcome with 5-year event-free survival rates ranging from 80.0% (39.4-100) to 88.2% (71.7-100). Patients with KMT2A-rearranged, BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like ALL had the worst outcomes, with 5-year event-free survival rates between 64.1% (43.9-84.3) and 76.2% (51.9-100), despite most patients receiving standard- or high-risk therapy. Treatment intensification based on MRD and genotype improved outcome. Relapse did not occur in any of the patients with ETV6-RUNX1, DUX4-rearranged, TCF3-PBX1, iAMP21, ZNF384-rearranged, hypodiploid, or BCR-ABL1-like B-ALL who had day 15 MRD≥1% and received standard-risk therapy. However, day 42 MRD<0.01% did not preclude relapse in patients with intermediate-risk or unfavorable subtypes including TCF3-PBX1, PAX5alt, T-cell, iAMP21, BCR-ABL1, BCR-ABL1-like, ETV6-RUNX1-like, KMT2A-rearranged, and MEF2D-rearranged ALL, despite intensified treatment in the standard- or high-risk arms. InterpretationContemporary risk-directed treatment based on genotypes and MRD can cure virtually all patients with ETV6-RUNX1, hyperdiploid>50 and DUX4-rearranged ALL. While intensified treatment improve outcome for patients with high MRD during induction, substantial proportions of patients with intermediate-risk or unfavorable subtypes would still relapse even achieving negative MRD status at the end of induction and require novel therapeutics to improve outcome.

  Study EGAS00001004739

• DIPG RNA and exome sequencing

  DIPG is a paediatric high-grade glioma that diffusely infiltrates the brainstem and is the leading cause of paediatric brain tumour death. Here we performed RNA and whole exome sequencing on DIPG and control samples to analyse transcriptomic and genetic alterations associated with this devastating disease, and to understand the relationship between them.

  Study EGAS00001004749

• Genetic ancestry contributes to somatic mutations in lung cancers from admixed Latin American populations

  Inherited lung cancer risk, particularly in non-smokers, is poorly understood. Genomic and ancestry analysis of 1,153 lung cancers from Latin America revealed striking associations between Native American ancestry and their somatic landscape, including tumor mutational burden (TMB), and specific driver mutations in EGFR, KRAS, and STK11. A local Native American ancestry risk score predicted EGFR and KRAS mutation frequency more strongly than global ancestry, suggesting that germline genetics (rather than environmental exposure) underlie these disparities.

  Study EGAS00001004752

• Second hit rare genetic variants in families with seemingly GBA gene associated Parkinson’s disease

  Study EGAS00001004777

• Neuroblastoma tumor heterogeneity and cell plasticity (from PDX and cell lines)

  Two cell identities, noradrenergic and mesenchymal, have been characterized in neuroblastoma cell lines according to their epigenetic landscapes relying on specific core regulatory circuitries of transcription factors. Yet, their relationship and relative contribution in patient tumors remain to be defined. Here, we document a spontaneous plasticity between the noradrenergic and mesenchymal identities in a subset of cell lines. Strikingly, mesenchymal neuroblastoma cells revert to a noradrenergic phenotype in vivo. Consistently, tumor cells with a mesenchymal identity are detected in single-cell transcriptomic analyses of neuroblastoma tumors and PDX models. Our data also highlight neuroblastoma intra-tumor heterogeneity with the co-existence of distinct tumor populations, including sympathoblast-like and chromaffin-like cells suggesting that neuroblastoma cells arise from a common sympatho-adrenal progenitor.

  Study EGAS00001004781

• Rna-Seq Leiomyosarcoma subtypes

  Leiomyosarcomas (LMS) are genetically heterogeneous tumors differentiating along smooth muscle lines. Currently, LMS treatment is not informed by molecular subtyping and is associated with highly variable survival. While disease site continues to dictate clinical management, the contribution of genetic factors to LMS subtype, origins, and timing are unknown. Here we analyzed 70 genomes and 130 transcriptomes of LMS, including multiple tumor regions and paired metastases. Molecular profiling highlighted the extremely early origins of LMS. Three specific subtypes of LMS develop from distinct lineages of smooth muscle cells. Of these, dedifferentiated LMS with high immune infiltration and tumors primarily of gynaecological origin acquired the highest burden of genomic mutation and are associated with worse survival. Homologous recombination defects lead to genome-wide mutational signatures, and a corresponding sensitivity to PARP and other DNA damage response inhibitors, suggesting a novel therapeutic strategy for LMS. Finally, by phylogenetic reconstruction, we present evidence that clones seeding lethal metastases arise decades prior to LMS diagnosis.

  Study EGAS00001004783

• Genomics-based personalized oncology in cancer of unknown primary (CUP, project H021)

  To perform a comprehensive genomic characterization of 70 patients suffering from cancer of unknown primary (CUP) we used whole-exome, whole-genome, transcriptome and methylome analysis. We detected a substantial mutational heterogeneity with genes most commonly affected by SNVs, indels and fusions being TP53, TTN, MUC16, ABCA13, COL6A3, KRAS, LRP1B, XIRP2 and CSMD3. The most common fusion involved FGFR2, the most common focal deletion affected CDKN2A. A molecular tumor board recommended genomics-based therapies in 56/70 (80%) patients which were applied in 20/56 (35.7%) cases. Entity predictions based on transcriptome and methylome data could be made in up to 62/70 (88.6%) cases but were conclusive in only 16/48 (33.3%) cases. Germline analysis revealed 6 (likely) pathogenic mutations in 5 patients. Recommended therapies translated into a mean PFS2/1 ratio of 3.61 (median=2.25) with a median PFS1 of 89 days (n=17) compared to a median PFS2 of 182.5 days (n=20). Our data emphasize the clinical benefit of comprehensive genetic approaches in diagnostic and therapeutic management and underline the need for innovative, mechanism-based clinical trials in this heterogeneous group of diseases.

  Study EGAS00001004786

• Methylation analysis for plasma DNA of patients with organ transplantation

  We developed Genetic-Epigenetic Tissue Mapping (GETMap) to determine the tissue composition of plasma DNA carrying genetic variants not present in the constitutional genome through comparing their methylation profiles with relevant tissues. We validated this approach by showing that, in pregnant women, circulating DNA carrying fetal-specific alleles was entirely placenta-derived.

  Study EGAS00001004788

• Familial psychosis associated with a missense mutation at MACF1 gene combined with the rare duplications DUP3p26.3 and DUP16p23.3, affecting the CNTN6 and CDH13 genes

  Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant's contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide polymorphisms (SNPs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C>G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16p23.3, were also identified in the family affecting relevant genes (CNTN6 and CDH13, respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C>G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

  Study EGAS00001004791

• Comprehensive Genomic and Transcriptomic Analysis of Rare Cancers for Guiding of Therapy (H021)

  Here we analyzed the molecular profiles and clinical outcomes of 1,310 patients (proportion of rare cancers, 75%) enrolled in a prospective observational study by the German Cancer Consortium that applies whole-genome/exome and RNA sequencing to inform the care of adults with incurable cancers. Based on 472 single and six composite biomarkers, a cross-institutional molecular tumor board provided evidence-based management recommendations, including diagnostic reevaluation, genetic counseling, and experimental treatment, in 88% of cases. Recommended therapies were administered in 362 of 1,138 patients (31.8%) and resulted in significantly improved response and disease control rates (23.9% and 55.3%, respectively) compared to previous therapies, translating into a progression-free survival ratio >1.3 in 35.7% of patients.

  Study EGAS00001004813

• Reconstruction of human phylogenetic trees using single-cell genome sequencing

  The trillions of cells that constitute the human body are developed from a fertilized egg through embryogenesis. However, cellular dynamics and developmental outcomes of embryonic cells remain to be largely unknown due to the technical and ethical challenges. Here we sequenced whole-genomes single-cell-derived clones and various bulk tissues dissected from entire human bodies. Using the discovered somatic mutations as an intrinsic barcode, we reconstructed cellular phylogenetic trees that provide novel insights into early human embryogenesis.

  Study EGAS00001004824

• The genetic structure of Norway

  The aim of the present study was to describe the genetic structure of the Norwegian population using genotypes from 6369 unrelated individuals with detailed information about places of residence. Using standard single marker- and haplotype-based approaches, we report evidence of two regions with distinctive patterns of genetic variation, one in the far northeast, and another in the south of Norway, as indicated by fixation indices, haplotype sharing, homozygosity and effective population size. We detect and quantify a component of Uralic Sami ancestry that is enriched in the North. On a finer scale, we find that rates of migration have been affected by topography like mountain ridges. In the broader Scandinavian context, we detect elevated relatedness between the mid- and northern border areas towards Sweden. The main finding of this study is that despite Norway’s long maritime history and as a former Danish territory, the region closest to mainland Europe in the south appears to have been the most isolated region in Norway, highlighting the open sea as a barrier to gene flow into Norway.

  Study EGAS00001004826

• Genomic analysis of a hypermutated gliosarcoma

  Gliosarcoma is a variant of glioblastoma with equally poor prognosis and characterized by mixed glial and mesenchymal pathology. Metastasis is not uncommon but involvement of the spinal cord is rare, and comprehensive genetic characterization of spinal gliosarcoma is lacking. We describe a patient initially diagnosed with a low-grade brain glioma via biopsy, followed by adjuvant radiation and temozolomide treatment. Nearly two years after diagnosis, she developed neurological deficits from an intradural, extramedullary tumor anterior to the spinal cord at T4, which was resected and diagnosed as gliosarcoma. Whole-exome sequencing (WES) of this tumor revealed a hypermutated phenotype, characterized by somatic mutations in key DNA mismatch repair (MMR) pathway genes, an abundance of C>T transitions within the identified somatic SNVs, and microsatellite stability, together consistent with temozolomide-mediated hypermutagenesis. This is the first report of a hypermutator phenotype in gliosarcoma, which may represent a novel genomic mechanism of progression from lower-grade glioma.

  Study EGAS00001004864

• Targeted de-methylation of the FOXP3-TSDR

  CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic-editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

  Study EGAS00001004867

• Genetic makeup of agnospheres

  Agnospheres are cultures of cancer stem cells isolated from tissue biopsies from patients with Cancer of Unknown Primary (CUPs). The latter is a lethal malignancy (median survival <1 year), representing approximately the 3% of all cancers, and characterized by metastatic dissemination in the absence of a primary tumor recognizable by a standardized protocol including thorough body imaging and histopathological analysis. Agnospheres are long-term propagated in culture and, after subcutaneous injection into immunocompromised mice, they display high tumorigenicity and quickly reproduce the clinical presentation of the original tumor, featuring a multi-organ metastatic pattern and poorly differentiated histology. The goal of the project is to analyse the mutational profile of agnospheres in detail, as to verify its correspondance with that of the original tumors. Genetic information will be exploited to investigate the genetic drivers of the CUP phenotype, in particular the hypermetastatic ability and the differentiative block.

  Study EGAS00001004868

• MutWP4__CRUK_Grand_Challenge_Mutographs_of_Cancer__Mutagen_treated_organoids

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, Kings College London will characterise the mutational signatures induced by putative human carcinogens in order to identify the origins of mutational signatures found in human cancers. To achieve this human organoid cell cultures will be exposed to a representative catalogue of known or suspected human carcinogens and mutagens and, using whole genome sequencing, the patterns of mutations induced by them will be determined. Somatic mutational signatures will be subsequently extracted by non-negative matrix factorisation methods and correlated with exposure data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001004873

• Targeted sequencing DDR genes in cancer stem cells

  Colorectal cancer stem cells (CSCs) show heterogeneous levels of constitutive replication stress leading distinct dependence on the replication stress response. The aim of the projects was to find genetic, epigenetic, phenotypic factors associated with the resistance of colorectal CSCs to inhibitors of the ATR/CHK1 pathway, which regulates of the replication stress response. Different pairs of colorectal CSCs sensitive and resistant to ATR or CHK1 inhibitors were analyzed

  Study EGAS00001004892

• Single cell multi-omic data of glioblastoma evolution under therapy

  We report a single-cell multi-omics atlas of GBM longitudinal specimens that represents decades of biobanking. We identify novel cell-intrinsic and cell-extrinsic targets for therapy development which we validate in situ and in low-passage patient-derived cell lines. We identify novel genetic correlates of the mesenchymal cellular phenotype. These studies integrate disparate findings in the literature around hypermutation status and immunogenicity as well as the source of the mesenchymal shift at recurrence.

  Study EGAS00001004909

• The Genetic History of Greenlandic-European contact

  The Inuit ancestors of the Greenlandic people arrived in Greenland close to 1000 years ago. Since then, Europeans from many different countries have been present in Greenland. Consequently, the present-day Greenlandic population has ~25% of its genetic ancestry from Europe. In this study, we investigated to what extent different European countries have contributed to this genetic ancestry.

  Study EGAS00001004933

• Single-cell RNA sequencing reveals the mesangial identity and species diversity of glomerular cell transcriptomes

  Molecular characterization of the individual cell types in human kidney as well as model organisms are critical in defining organ function and understanding translational aspects of biomedical research. Previous studies have uncovered gene expression profiles of several kidney glomerular cell types, however, important cells, including mesangial (MCs) and glomerular parietal epithelial cells (PECs), were missing or incompletely described, and a systematic comparison between mouse and human kidney is lacking. To this end, we used Smart-seq2 to profile 4332 individual glomerulus-associated cells isolated from human living donor renal biopsies and mouse kidney. The analysis revealed genetic programs for all four glomerular cell types (podocytes, glomerular endothelial cells, MCs and PECs) as well as rare glomerulus-associated macula densa cells (MDCs). Importantly, we detected heterogeneity in glomerulus-associated Pdgfrb-expressing cells, including bona fide intraglomerular MCs with the functionally active phagocytic molecular machinery, as well as a unique mural cell type located in the central stalk region of the glomerulus tuft. Furthermore, we observed remarkable species differences in the individual gene expression profiles of defined glomerular cell types that highlight translational challenges in the field and provide a guide to design translational studies.

  Study EGAS00001004943

• Genetic insights into the biological mechanisms governing human ovarian ageing

  Reproductive longevity is critical for fertility and impacts healthy ageing in women, yet insights into the underlying biological mechanisms and treatments to preserve it are limited. Here, we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in ~200,000 women of European ancestry. These common alleles influence clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations. Identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increase fertility and extend reproductive life in mice. Causal inference analyses using the identified genetic variants indicates that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases risks of hormone-sensitive cancers. These findings provide insight into the mechanisms governing ovarian ageing, when they act across the life-course, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

  Study EGAS00001004947

• Comparison of the diagnostic yield of aCGH and NGS across different neurodevelopmental disorders

  Current recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). The present study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: Global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and Other NDDs. These categories were further subclassified into subcategories, based on the most frequent accompanying signs and symptoms: isolated forms, forms with epilepsy; forms with micro or macrocephaly and syndromic forms. A subset of 245 patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results support the use of NGS as the first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

  Study EGAS00001004949

• Cell-free DNA TAPS for early cancer detection

  Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection. However, due to technological challenges associated with DNA methylation sequencing in low input cfDNA samples, most studies have been limited by DNA damage caused by bisulfite sequencing, or the qualitative nature of enrichment-based sequencing. Recently, we developed TET-assisted Pyridine Borane Sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and non-cancer controls. From just 10 ng cfDNA (1-3 mL of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

  Study EGAS00001004962

• The genomic landscape of childhood acute lymphoblastic leukaemia with intrachromosomal amplification of chromosome 21 (iAMP21-ALL)

  Intrachromosomal amplification of chromosome 21 (iAMP21) results from breakage-fusion-bridge cycles and chromothripsis defines a cytogenetic subgroup of childhood B-cell precursor acute lymphoblastic leukemia (B-ALL) cases associated with a poor prognosis. To systematically characterize the mechanism of somatic alternations in iAMP21, we perform integrated genomic, cytogenetic and transcriptomic analysis of 128 cases.

  Study EGAS00001004998

• Genetic heterogeneity and dynamics of transcriptional subtypes in matched primary and recurrent head and neck squamous cell carcinomas

  To get a better understanding of different HNSCC subtypes, whole-exome (matched tumor-normal paits) and whole transcriptome sequencing to to elucidate underlying molecular processes. Here, we present the matched whole-exome sequencing data from 28 patients (normal tissue, tumor tissue, relapse).

  Study EGAS00001005005

• Primary prostate Hi-C

  Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. Despite large-scale tumour sequencing efforts, the impact of mutations on the genetic architecture in cancer remains ill-defined due to limited integration of genomics data across dimensions. We addressed this limitation by assessing the impact of structural variants on the chromatin states and the three-dimensional organization across benign and malignant primary prostate genomes. We find high concordance in the three-dimensional genome organization between malignant and benign prostate tissues, arguing for constraints to the three-dimensional genome of prostate tumours. Moreover, we identify structural variants as effectors of changes to focal chromatin interactions, guiding cis-regulatory element hijacking that imposes opposing expression changes on genes found at antipodes of a rearrangement. This leads to the repression of tumour suppressor gene expression and up-regulation of oncogenes, such as at the TMPRSS2-ERG and PMEPA1-ZNF156 loci. Collectively, our results argue that cis-regulatory element hijacking by structural variants overshadows large-scale topological changes to alter gene regulation and promote oncogenesis.

  Study EGAS00001005014

• Exome sequencing reveals pathogenic variants in known and novel candidate genes for severe sperm motility disorders

  STUDY QUESTION: What are the causative genetic variants in patients with male infertility due to severe sperm motility disorders? SUMMARY ANSWER: We identified high confidence disease-causing variants in multiple genes previously associated with severe sperm motility disorders in 9 out of 21 patients (43%). We also identified variants in novel candidate genes in 8 additional patients (38%). WHAT IS KNOWN ALREADY: Severe sperm motility disorders are a form of male infertility characterised by vital, but immotile sperm often in combination with a spectrum of structural abnormalities of the sperm flagellum. Examples of this spectrum of disorders include dysplasia of the fibrous sheath / multiple morphological abnormalities of the sperm flagella (DFS-MMAF) and severe forms of asthenozoospermia associated to non-specific structural defects. These disorders frequently recur in families and may include chronic respiratory disease. As such, genetic causes are suspected and many genes have recently been reported. Currently, depending on the clinical sub-categorisation, up to 50% of causality in patients with severe sperm motility disorders can be explained by pathogenic variants in at least 21 known genes. STUDY DESIGN, SIZE, DURATION: We performed exome sequencing in patients with severe sperm motility disorders from two different clinics, aiming to identify the underlying genetic defects and establishing a diagnostic yield in both. PARTICIPANTS/MATERIALS, SETTING, METHOD: Two groups of infertile men, one from Argentina (n= 9) and one from Australia (n=12), with clinically defined sperm motility disorders were included. All patients in the Argentine cohort were diagnosed with DFS-MMAF, based on transmission electron microscopy. Exome sequencing was performed in all 21 patients and variants with an allele frequency of <1% in the gnomAD population were prioritised and interpreted. MAIN RESULTS AND ROLE OF CHANCE: In 9/21 patients (43%), we identified pathogenic variants in known sperm motility disorder genes: CFAP43 (3 patients); CFAP44 (2 patients), QRICH2 (2 patients), DNAH1 (1 patient) and DNAH6 (1 patient). The diagnostic rate did not differ markedly between the Argentinian and the Australian cohort (44% and 42%, respectively). Furthermore, we identified patients with variants in the novel candidate sperm motility genes CFAP58, DNAH12, DRC1, MDC1, PACRG, SSPL2C and TPTE2. One patient presented with two possibly damaging variants in four candidate genes and it remains unclear which variants were responsible for the severe sperm motility defect in this patient. LIMITATIONS, REASONS FOR CAUTION: In this study, we described patients with potential bi-allelic variants in known and novel candidate genes for severe sperm motility disorders. Due to unavailability of parental DNA, we have not assessed de novo or maternally inherited dominant variants and could not phase the mutations to establish in all cases that the mutations occur on both alleles. WIDER IMPLICATIONS OF THE FINDINGS: Our results confirm an important role for five known genes for sperm motility and we demonstrate that exome sequencing is an effective method to diagnose patients with severe sperm motility defects (9/21 diagnosed; 43%). Furthermore, our analysis revealed seven novel candidate genes for severe sperm motility disorders. Genome-wide sequencing of additional patient cohorts and re-analysis of exome data of currently unsolved cases may reveal additional variants in these novel candidate genes.

  Study EGAS00001005018

• Transcriptomic and epigentic analysis of human peripheral blood NK cell subsets revealing role of Bcl11b in differentiation

  Transcriptomic analysis of seven human peripheral blood NK cell subsets and three T cell subsets purified from healthy cytomeglovirus seropositive donors. Epigenetic analysis including identification of open chromatin and active enhancer elements confirms key regulatory factors defining an NK cell differentiation axis. Chromatin-immunoprecipitation sequencing reveals contribution of key transcription factors associated with NK cell differentiation including Bcl11b.

  Study EGAS00001005025

• Summary statistics of meta-analysis using two genome-wide association study of inflammatory bowel disease in Koreans.

  Inflammatory bowel disease (IBD), a chronic inflammatory disorder of the gastrointestinal tract, is thought to develop due to dysregulated mucosal immune responses to gut flora in genetically susceptible individuals. Crohn’s disease (CD) and ulcerative colitis (UC) are the two major subtypes of IBD. To identify additional susceptibility loci for IBD in Asians, we performed meta-analyses using two genome-wide association studies.

  Study EGAS00001005026

• Test Study for EGA using data from 1000 Genomes Project - Phase 3

  The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies...

  Study EGAS00001005042

• Genetic analysis of HLA and immune escape genes in Diffuse Large B-cell Lymphoma

  Fifty percent of diffuse large B-cell lymphoma (DLBCL) cases lack cell-surface expression of the class-I major histocompatibility complex (MHC-I), thereby escaping immune recognition by cytotoxic T cells. In order to comprehensively identify the mechanisms involved in MHC-I loss, we first performed immunophenotypic analysis of both MHC class-I and -II in 657 cases across the spectrum of B-NHL, revealing that loss of MHC-I, but not MHC-II, is preferentially restricted to DLBCL. We then used whole exome and targeted deep-sequencing to examine genes involved in MHC-I expression in 74 DLBCL samples representative of MHC-I positive and negative cases. We show here that somatic biallelic or monoallelic inactivation of B2M and/or HLA-I is present in 80.9% (34/42) of MHC-I negative tumors. Furthermore, 68.8% (22/32) of MHC-I positive DLBCLs also harbored monoallelic HLA-I genetic alterations (MHC-I positivemono) that lead to allelic imbalance, suggesting allele-specific inactivation. Both MHC-I negative and MHC-I positivemono cases showed a significantly higher mutational burden as well as a higher number of inferred neo-antigens, suggesting co-selection of HLA-I loss and sustained neo-antigen production. Interestingly, the analysis of > 500.000 individuals in two databases revealed homozygosity of germline HLA-I genes in 26-38% of DLBCL patients, a frequency significantly higher than that observed in any other cancer type. In mice, germinal-center B cells lacking HLA-I expression did not progress to lymphoma and were counterselected in the context of oncogene-driven lymphomagenesis, suggesting that additional events are needed to license immune evasion. These results suggest a multi-step process of HLA-I loss including both germ-line and somatic events in DLBCL development, and have direct implications for the pathogenesis and immunotherapeutic targeting of this disease.

  Study EGAS00001005054

• Characterization of DLBCL with a PMBL gene expression signature

  Primary mediastinal large B-cell lymphoma (PMBL) is a type of aggressive B-cell lymphoma that typicallyaffects young adults, characterized by presence of a bulky anterior mediastinal mass. Lymphomas withgene expression features of PMBL have been described in non-mediastinal sites, raising questions abouthow these tumors should be classified.Here, we investigated whether these "non-mediastinal PMBLs" are indeed PMBLs or instead represent adistinct group within DLBCL. From a cohort of 325 de novo DLBCL cases, we identified tumors frompatients without evidence of anterior mediastinal involvement that expressed a PMBL expression signature(nm-PMBLsig-pos, n=16, 5%). The majority of these tumors expressed MAL and CD23 – proteins typicallyobserved in bona fide PMBL (bf-PMBL). Evaluation of clinical features of nm-PMBLsig-pos cases revealedclose associations with DLBCL, and the majority displayed a germinal center B-cell-like cell-of-origin(GCB). In contrast to bf-PMBL, nm-PMBLsig-pos patients presented at an older age, did not show pleuraldisease, and bone/bone marrow involvement was observed in three cases. However, while clinicallydistinct from bf-PMBL, nm-PMBL-sig-pos tumors resembled bf-PMBL at the molecular level with upregulationof immune response, JAK-STAT, and NF-kB signatures. Mutational analysis revealed frequent somatic genemutations in SOCS1, IL4R, ITPKB and STAT6, as well as CD83 and BIRC3, with the latter genes beingsignificantly more frequently affected than in GCB-DLBCL and bf-PMBL.Our data establish nm-PMBLsig-pos lymphomas as a group of DLBCL with distinct phenotypic and geneticfeatures, and potential implications for gene expression- and mutation-based subtyping of aggressive Bcell lymphoma and related targeted therapies.

  Study EGAS00001005057

• Mutational analysis of an oligoprogressive sarcomatoid hepatocellular carcinoma treated with an immune checkpoint inhibitor.

  In this study, we describe a patient diagnosed with metastatic sarcomatoid hepatocellular carcinoma that showed a remarkable tumor response to Nivolumab, an immune checkpoint inhibitor, followed by an oligoprogression. Molecular analysis of the primary tumor in addition to a bone and brain metastasis was performed and we were able to identify driver genetic alterations indicating a possible role in the tumorigenesis process besides some alterations that may be associated with metastasis.

  Study EGAS00001005064

• Clonal evolution and clinical implications of genetic abnormalities in blastic transformation of chronic myeloid leukaemia

  Blast crisis (BC) predicts dismal outcomes in patients with chronic myeloid leukaemia (CML). We comprehensively investigated genetic abnormalities in 136 BC and 148 chronic phase (CP) samples obtained from 216 CML patients using exome and targeted sequencing.

  Study EGAS00001005075

• Multiple migrations to the Philippines during the last 50,000 years

  A key link to understand human history in Island Southeast Asia is the Philippine archipelago and its poorly investigated genetic diversity. We collected and analyzed the most comprehensive set of population-genomic data for the Philippines, 1028 individuals covering 115 indigenous communities. We demonstrate that the Philippines were populated by at least five waves of human migration. The Cordillerans migrated into the Philippines prior to the arrival of rice agriculture, where some remain to be the least admixed East Asians carrying an ancestry shared by all Austronesian-speaking populations; thereby challenging an exclusive Out-of-Taiwan model of joint farming-language-people dispersal. Altogether, our findings portray Philippines as a crucial gateway, with a multilayered history, that ultimately changed the genetic landscape of the Asia-Pacific region.

  Study EGAS00001005083

• single-stranded DNA study

  A ssDNA library protocol was applied to cfDNA from plasma samples obtained from different DNA extraction methods and revealed significant differences in DNA fragmentation patterns in comparison to dsDNA-based protocols. In particular, a specific combination of methods revealed a population of ultrashort fragments, organized at ~50 bp. We observed significant differences in the relative abundance of these ultrashort DNA fragments in plasma from healthy individuals and cancer patients. Through shallow whole genome sequencing (sWGS, <0.5-fold coverage) and the analysis of somatic copy number aberrations (SCNA), we determined the landscape of genetic alterations in this newly identified population of cfDNA fragments. In addition, we studied their potential link with regulatory regions by investigating the genome-wide coverage patterns at transcription start sites (TSS). Furthermore, we demonstrated that the ultrashort cfDNA fragments map to regions associated with secondary DNA structures, G-quadruplexes (G4s).

  Study EGAS00001005093

• Whole genome sequencing of 76 tumor and normal samples from 11 SI-NET patients

  We have sampled 11 patients with multifocal SI-NETs. The complete genome sequences of 65 tumor and tissue samples including 42 intestinal lesions, 14 adjacent lymph node metastases, 5 distant metastases, as well as of 4 normal tissue and 11 normal blood samples were determined. This allowed elucidation of phylogenetic relationships between tumors within individual patients.

  Study EGAS00001005096

• Genetic landscape of inherited retinal dystrophies affected cases in Spain

  Inherited retinal diseases (IRDs), defined by dysfunction or progressive loss of photoreceptors, are disorders characterized by elevated heterogeneity, both at the clinical and genetic levels. Our main goal was to address the genetic landscape of IRD in a Spanish cohort of patients. A retrospective hospital-based cross-sectional study was carried out on IRD affected individuals, referred for genetic testing from all the Spanish autonomous communities. Clinical, demographic and familiar data were collected from each patient, including family pedigree, age of appearance of visual symptoms, presence of any systemic findings and geographical origin. Genetic studies were performed to all families with available DNA using different molecular techniques. Our study provides the general landscape for IRD in Spain and our results have important implications for genetic diagnosis, counselling and new therapeutic strategies to both the Spanish population and other related populations.

  Study EGAS00001005104

• Chronic myelomonocytic leukemia

  Chronic myelomonocytic leukemia is an aggressive hematological malignancy with dismal outcomes, with the pCMML subtype in particular having median OS of <2 years and high rates of AML transformation (Patnaik et al., 2014). Given limited therapeutic options for affected patients, we carried out this study to define the genetic and epigenetic landscape of pCMML and identify therapies that could modify disease biology.

  Study EGAS00001005107

• Sequencing of cancer autopsies and ctDNA

  Tumor heterogeneity is a major obstacle for precision oncology strategies. To characterize inter-lesion heterogeneity in pre-treated cancers, we performed genetic and transcriptomic profiling of synchronous lesions from ten rapid autopsies comprising six cases of breast cancer, three of cholangiocarcinoma and one lung cancer case. For nine patients, we studied mutation representation in circulating cell-free DNA (cfDNA).

  Study EGAS00001005109

• International Genetics of Parkinson Disease Progression (IGPP) Consortium GWAS Summary Results

  To characterize how genetic variation influences the progression of Parkinson’s disease to dementia, which is a major determinant for quality of life, we performed a longitudinal genome-wide survival study (GWSS). Risk profile analysis was performed with a novel polygenic hazard score in longitudinal cohorts in the International Genetics of Parkinson Disease Progression (IGPP) Consortium.

  Study EGAS00001005110

• Accurate detection and classification of pediatric sarcomas based on cell-free DNA fragmentation patterns

  Sequencing of cell-free DNA in the blood of cancer patients (“liquid biopsy”) provides attractive opportunities not only for early diagnosis, but also for minimally invasive monitoring of treatment response and disease courses. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we developed an integrated genetic/epigenetic analysis method and applied it to 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. We achieved sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. We evaluated different metrics for cell-free DNA fragmentation analysis and developed LIQUORICE, a bioinformatic tool for detecting circulating tumor DNA based on tumor-specific chromatin structure. Using machine learning methods, we combined several fragmentation-based metrics into an integrated approach for liquid biopsy analysis tailored to cancers with low mutation rates but widespread epigenetic deregulation. Clinical associations highlighted the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. Additionally, we performed low coverage whole-genome-sequencing on 43 tumor biopsy samples from patients with Ewing sarcoma, in order to compare copy number aberrations detected in cell-free DNA and biopsy samples of the same patients. For validation of the epigenetic signatures inferred from cell-free DNA, we further performed reduced representation bisulfite sequencing (RRBS) on 38 matched biopsy samples from patients with Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent somatic mutations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

  Study EGAS00001005127

• Multiregional sequencing of IDH-WT glioblastoma reveals high genetic heterogeneity and a dynamic evolutionary history

  Glioblastoma is the most common and aggressive primary brain malignancy in adults. In addi-tion to extensive inter-patient heterogeneity, glioblastoma shows intra-tumor extensive cellular and molecular heterogeneity, both spatially and temporally. This heterogeneity is one of the main reasons for the poor prognosis and overall survival. Moreover, it raises the important question of whether the molecular characterization of a single bioptic sample, as performed in standard di-agnostics, actually represents the entire lesion. In this study, we sequenced the whole-exome of nine spatially different cancer regions of three primary glioblastomas. We characterized their mutational profiles and copy number alterations, with implications for our understanding of tu-mor biology in relation to clonal architecture and evolutionary dynamics, as well as therapeuti-cally relevant alterations.

  Study EGAS00001005128

• Epigenomic profiling and transcriptomic profiling of LUAD patients tumor tissues and tumor adjacent tissues

  To provide an in-depth understanding of the epigenomic heterogeneity of LUAD, we here investigated the H3K27ac histone modification profiles of tumors and adjacent normal lung tissues from 42 LUAD patients and explored the role of epigenetic alterations in LUAD progression. We also investigated transcriptomic alterations among 36 patients tumor tissues.

  Study EGAS00001005132

• Pediatric Papillary Thyroid Carcinoma Whole Exome Sequencing

  BACKGROUND: Pediatric papillary thyroid carcinoma (PPTC) is distinct from adult-onset disease. Despite higher rates of regional and distal metastases and more frequent recurrence, there is, paradoxically, a highly favourable prognosis. Historically, only 40-50% of childhood PTC are driven by genomic variants common to adult PTC, while the oncogenic drivers in the remainder are unknown. METHODS: Fifty-two PPTCs were genotyped by candidate gene testing, followed, in a subset (n=32), by whole exome (WES) and transcriptome sequencing (RNASeq) to identify oncogenic variants in driver-unknown tumors and to explore gene expression patterns. FINDINGS: Thirty-one tumors (60%) carried variants identified by candidate gene testing, specifically: BRAFV600E(22%), NRASQ61R(2%), RET-CCDC6(19%), RET-NCOA4(15%) and PAX8-PPARG. Among the 21 remaining tumors with negative first-line genotyping, oncogenic variants were identified in 19 (90%). These were enriched with oncogenic fusions, with 11 non-recurrent fusion transcripts identified, including two that had not been previously described: STRN-RET and TG-PBF. Most fusions were associated with 3’ receptor tyrosine kinase (RTK) moieties: RET, MET, ALK and NTRK3. Gene expression analysis defined three distinct clusters, which demonstrated differing expression of genes involved in thyroid differentiation and MAPK signaling. INTERPRETATION: Exome and transcriptome-wide sequencing identify genomic drivers of PPTC in 98% of tumors, whereas conventional genotyping leaves many tumors genetically unexplained. PPTC is driven, in the majority of cases, by oncogenic fusion transcripts involving one of a limited number of RTKs. These observations have broad implications for diagnostic testing and for therapeutic options in children with advanced disease.

  Study EGAS00001005187

• mRNA capture sequencing and RT-qPCR for the detection of pathognomonic, novel and secondary fusion transcripts in formalin-fixed paraffin-embedded tissue: a sarcoma showcase

  RNA sequencing is widely used to study gene abundance patterns, but less frequently to assess structural alterations such as fusion transcripts. In this study, we assess the performance of messenger RNA (mRNA) capture sequencing to identify fusion transcripts in formalin-fixed paraffin-embedded (FFPE) tissue of different types of sarcoma, followed by orthogonal confirmation using RT-qPCR. To validate our workflow, we used the TruSight RNA Pan-Cancer Panel to analyze a first cohort of 6 positive control tumors for which the diagnostic workup demonstrated the presence of a specific chromosomal rearrangement. Fusion transcript calling by FusionCatcher confirmed all known aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole transcriptome TruSeq RNA Exome sequencing was applied to a second cohort of 17 patients, diagnosed with fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS), for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. Remarkably, for 6 out of 17 patients, mRNA capture sequencing readily detected a pathognomonic fusion transcript, i.e. PAX3-FOXO1 in 2 ARMS patients, and EWSR1-FLI1, EWSR1-ERG or EWSR1-NFATC2 in 4 URCS patients. These data show that RNA capture sequencing may enhance the detection rate of pathognomonic fusion genes in sarcoma. For the 11 patients with no known pathognomonic fusion transcript, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2 and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in two patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53 and UBE2F-ALS2). Our results provide new insights into the underlying genetic causes of these malignancies.

  Study EGAS00001005202

• Novel optineurin frameshift insertion causing familial frontotemporal dementia and parkinsonism without amyotrophic lateral sclerosis

  Frontotemporal Dementia (FTD) is a common cause of Young Onset Dementia and has diverse clinical manifestations involving behaviour, executive function, language and motor function, including parkinsonism. Up to 50% of FTD patients report a positive family history, supporting a strong genetic basis, particularly in cases with both FTD and amyotrophic lateral sclerosis (FTD-ALS). Mutations in three genes are associated with the majority of familial FTD (fFTD) cases - microtubule associated protein tau gene (MAPT), progranulin gene (PGRN), and hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72), whilst mutations in other genes such as optineurin (OPTN) have rarely been reported. Mutations in OPTN have been reported mostly in familial and sporadic cases of ALS, or in rare cases of FTD-ALS, but not in association with pure or predominant FTD and/or parkinsonian phenotype. Here, we report for the first time, a family from the Philippines with 4 affected members harbouring the same OPTN frameshift insertion, presenting with FTD-related phenotypes, including one sibling with predominant parkinsonism resembling corticobasal syndrome. Notably, none of the affected members showed any evidence of motor neuron disease or ALS at the time of writing, both clinically and on electrophysiological testing, expanding the phenotypic spectrum of OPTN mutations. Close follow-up of mutation carriers for the development of new clinical features and wider investigation of additional family members with further genetic analyses will be conducted to investigate the possibility of other genetic modifiers in this family which could explain phenotypic heterogeneity.

  Study EGAS00001005220

• Single-cell landscapes of primary glioblastomas and matched organoids and cell lines reveal variable retention of inter- and intra-tumor heterogeneity

  Glioblastomas (GBMs) are aggressive primary malignant brain tumors characterized by extensive levels of inter- and intra-tumor genetic and phenotypic heterogeneity. Patient-derived organoids (PDOs) have recently emerged as useful models to study such heterogeneity. Here, we present bulk exome as well as single-cell genome and transcriptome profiles of primary IDH wild type GBMs from ten patients, including two recurrent tumors, as well as PDOs and brain tumor-initiating cell (BTIC) lines derived from these patients. We find that PDOs are genetically similar to and variably retain gene expression characteristics of their parent tumors. At the phenotypic level, PDOs appear to exhibit similar levels of transcriptional heterogeneity as their parent tumors, whereas BTIC lines tend to be enriched for cells in a more uniform transcriptional state. The datasets introduced here will provide a valuable resource to help guide experiments using GBM-derived organoids, especially in the context of studying cellular heterogeneity.

  Study EGAS00001005227

• Circulating tumor cells Exome sequencing from breast cancer

  Single cell technologies allow interrogating the extent of tumor heterogeneity, providing novel insights into tumor evolution and treatment resistance mechanisms. Here, we investigated the level of genetic heterogeneity among circulating tumor cells (CTCs) at the single nucleotide variation and copy number aberration (CNA) levels and compared it to synchronous matched primary tumor/metastatic lesion from 3 metastatic breast cancer patients with triple-negative (TNBC), HER2-positive (HER2+) and estrogen-receptor-positive (ER+) tumors. Whole exome sequencing was performed on 3 synchronous bulk tumor tissue and 21 single/pooled CTC samples. When comparing bulk tissue with CTCs from the same patient, we found similar CNA profiles and the same patient-specific driver mutations for the HER2+ and TNBC tumors. Different CNA profiles and driver mutations were found for the ER+ tumor, with two distinct clones being identified in the CTCs, one harboring an activating ESR1 mutation, the other one bearing a TP53 mutation. These data suggest that tumor tissue and CTCs provide in a subset of patients complementary clinically relevant information to map tumor heterogeneity and monitor tumor evolution. Further validation is needed.

  Study EGAS00001005228

• Non-coding mutations reveal cancer driver cistromes in luminal breast cancer

  Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and non-coding cancer driver plexuses from somatic mutations. However, differentiating between driver and passenger events among non-coding genetic variants remains a challenge to understand the etiology of cancer and inform the delivery of personalized cancer medicine. Herein, we reveal enrichment of non-coding mutations in cis-regulatory elements that cover a subset of transcription factors linked to tumor progression in luminal breast cancers. Using a cohort of 26 primary luminal ER+PR+ breast tumors, we compiled a catalogue of ~100,000 unique cis-regulatory elements from ATAC-seq data. Integrating this catalogue with somatic mutations from 350 publicly available breast tumor whole genomes, we identified four recurrently mutated individual cis-regulatory elements. By then partitioning the non-coding genome into cistromes, defined as the sum of binding sites for a transcription factor, we uncovered cancer driver cistromes for ten transcription factors, namely CTCF, ELF1, ESR1, FOSL2, FOXA1, FOXM1 GATA3, JUND, TFAP2A, and TFAP2C in luminal breast cancer. Nine of these ten transcription factors were shown to be essential for growth in breast cancer, with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of non-coding mutations over cis-regulatory elements concatenated into a functional unit drawn from an accessible chromatin catalogue derived from primary cancer tissues.

  Study EGAS00001005235

• Targeted sequencing of follicular lymphoma tumour samples from the UK's Haematological Malignancy Research Network

  Follicular lymphoma displays morphological and clinical heterogeneity, underpinned by shared molecular pathogenesis. Mutations in epigenetic regulators alongside t(14;18) have been identified as disease initiating events. Additional mutational diversity suggests disease heterogeneity. Whether this resolves into mechanistically distinct subsets remains an open question. Targeted sequencing was applied to an unselected population-based FL cohort (n=548) diagnosed in a catchment population of ~4 million (14 centers), with full clinical follow-up (n=538) that included 96 DLBCL transformations. We investigated whether molecular subclusters of FL can be identified, and whether mutational data provide predictive information relating to transformation.

  Study EGAS00001005238

• Radiation-induced gliomas represent H3-/IDH-wild type pediatric gliomas with recurrent PDGFRA amplification and loss of CDKN2A/B

  Long-term complications such as radiation-induced second malignancies occur in a subset of patients following radiation-therapy, particularly relevant in pediatric patients due to the long follow-up period in case of survival. Radiation-induced gliomas (RIGs) have been reported in patients after treatment with cranial irradiation for various primary malignancies such as acute lymphoblastic leukemia (ALL) and medulloblastoma (MB). We perform comprehensive (epi-) genetic and expression profiling of RIGs arising after cranial irradiation for MB (n=23) and ALL (n=9). Our study reveals a unifying molecular signature for the majority of RIGs, with recurrent PDGFRA amplification and loss of CDKN2A/B and an absence of somatic hotspot mutations in genes encoding histone 3 variants or IDH1/2, uncovering diagnostic markers and potentially actionable targets.

  Study EGAS00001005243

• The_Little_Princess_Knowledge_Bank_RNAseq

  The purpose of the knowledge bank is to serve as a community resource and to address the following specific questions about Wilms tumour: (1) Do genetic, epigenetic, or transcriptional features predict outcome and thus lend themselves to treatment stratification? (2) What are the genetic drivers of Wilms tumour, considering the entire histological spectrum of Wilms tumour including the most common subtypes classified as intermediate risk?

  Study EGAS00001005244

• Host genetic determinants of HIV infection

  Broadly, the objectives of the study are to characterize genetic variation in any gene or genomic region in South African populations, in the first place HIV. The purpose is to apply this knowledge to establish assays and gene signatures to test in disease-association studies, and to study underlying mechanisms of disease causation/severity.

  Study EGAS00001005245

• Multimodal single-cell and bulk glioma analyses

  Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes, and bulk multi-omic profiles across 11 adult IDH-mutant or IDH-wild-type gliomas to delineate sources of intratumoral heterogeneity. We show that local DNA methylation disorder associates with cell-to-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption, and is altered in environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression, and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identifies an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

  Study EGAS00001005300

• COVID19 Host Genetic Initiative

  The COVID-19 host genetics initiative is a bottom-up collaborative effort that has three main goals: 1. Provide an environment to foster the sharing of resources to facilitate COVID-19 host genetics research (e.g. protocols, questionnaires). 2. Organize analytical activities across studies to identify genetic determinants of COVID-19 susceptibility and severity. 3. Provide a platform to share the results from such activities, as well as the individual-level data where possible, to benefit the broader scientific community. Here, we release the harmonized individual-level data for the subset of individuals who participated in COVID-19 HGI. Please refer to this data dictionary: https://docs.google.com/spreadsheets/d/1whu_MjShXxbdoAJMCYxZ5C-vWOcRivwfWglTcBRNZgQ/edit#gid=470747846 Please cite both Nakanishi, T el al. J Clin Invest 2021 (PMID: 34597274) and COVID-19 Host Genetics Initiative, Nature 2021 (PMID: 34237774) when you use these data for your publications.

  Study EGAS00001005304

• Bulk-tissue RNA-sequencing of anterior cingulate cortex samples derived from Lewy body disease patients

  Parkinson’s disease (PD), Parkinson’s disease with dementia (PDD) and dementia with Lewy bodies (DLB) are three clinically, genetically and neuropathologically overlapping neurodegenerative diseases collectively known as the Lewy body diseases (LBDs). A variety of molecular mechanisms have been implicated in PD pathogenesis, but the mechanisms underlying PDD and DLB remain largely unknown, a knowledge gap that presents an impediment to the discovery of disease-modifying therapies. Transcriptomic profiling can contribute to addressing this gap, but remains limited in the LBDs. Here, we applied paired bulk-tissue and single-nucleus RNA-sequencing to anterior cingulate cortex samples derived from 28 individuals, including healthy controls, PD, PDD and DLB cases (n = 7 per group), to transcriptomically profile the LBDs. Using this approach, we (i) found transcriptional alterations in multiple cell types across the LBDs; (ii) discovered evidence for widespread dysregulation of RNA splicing, particularly in PDD and DLB; (iii) identified potential splicing factors, with links to other dementia-related neurodegenerative diseases, coordinating this dysregulation; and (iv) identified transcriptomic commonalities and distinctions between the LBDs that inform understanding of the relationships between these three clinical disorders. Together, these findings have important implications for the design of RNA-targeted therapies for these diseases and highlight a potential molecular “window” of therapeutic opportunity between the initial onset of PD and subsequent development of Lewy body dementia.

  Study EGAS00001005305

• Evolutionary predictability of genetic versus non genetic resistance to anticancer drugs in melanoma

  Therapy resistance arises from heterogeneous drug-tolerant persister cells or minimal residual disease (MRD) through genetic and nongenetic mechanisms. A key question is whether specific molecular features of the MRD ecosystem determine which of these two distinct trajectories will eventually prevail. We show that, in melanoma exposed to MAPK therapeutics, emergence of a transient neural crest stem cell (NCSC) population in MRD concurs with the development of nongenetic resistance. This increase relies on a GDNF-dependent signaling cascade, which activates the AKT survival pathway in a focal adhesion kinase (FAK)-dependent manner. Ablation of the NCSC population through FAK inhibition delays relapse in patient-derived tumor xenografts. Strikingly, all tumors that ultimately escape this treatment exhibit resistance-conferring genetic alterations and increased sensitivity to ERK inhibition. These findings identify an approach that abrogates the nongenetic resistance trajectory in melanoma and demonstrate that the cellular composition of MRD deterministically imposes distinct drug resistance evolutionary paths.

  Study EGAS00001005314

• Neuroblastoma tumor heterogeneity and cell plasticity (from patients)

  Two cell identities, noradrenergic and mesenchymal, have been characterized in neuroblastoma cell lines according to their epigenetic landscapes relying on specific core regulatory circuitries of transcription factors. Yet, their relationship and relative contribution in patient tumors remain to be defined. Here, we document a spontaneous plasticity between the noradrenergic and mesenchymal identities in a subset of cell lines. Strikingly, mesenchymal neuroblastoma cells revert to a noradrenergic phenotype in vivo. Consistently, tumor cells with a mesenchymal identity are detected in single-cell transcriptomic analyses of neuroblastoma tumors and PDX models. Our data also highlight neuroblastoma intra-tumor heterogeneity with the co-existence of distinct tumor populations, including sympathoblast-like and chromaffin-like cells suggesting that neuroblastoma cells arise from a common sympatho-adrenal progenitor.

  Study EGAS00001005322

• Multi-site tumor sampling highlights molecular intra-tumor heterogeneity in malignant pleural mesothelioma

  Background:Malignant pleural mesothelioma (MPM) is a heterogeneous cancer. Better knowledge of molecular and cellular intra-tumor heterogeneity throughout the thoracic cavity is required to develop efficient therapies. This study focuses on molecular intra-tumor heterogeneity using the largest series to date in MPM and is the first to report on the multi-omics profiling of a substantial series of multi-site tumor samples. Methods:Intra-tumor heterogeneity was investigated in 16 patients from whom biopsies were taken at distinct anatomical sites. The paired biopsies collected from apex, side wall, costo-diaphragmatic or highest metabolic sites as well as 5 derived cell lines were screened using targeted sequencing. Whole exome sequencing, RNA sequencing, and DNA methylation were performed on a subset of the cohort for deep characterization. Molecular classification, recently defined histo-molecular gradients, and cell populations of the tumor microenvironment were assessed. Results:Sequencing analysis identified heterogeneous variants notably in NF2, a key tumor suppressor gene of mesothelial carcinogenesis. Subclonal tumor populations were shared among paired biopsies, suggesting a polyclonal dissemination of the tumor. Transcriptome analysis highlighted dysregulation of cell adhesion and extracellular matrix pathways, linked to changes in histo-molecular gradient proportions between anatomic sites. Methylome analysis revealed contribution of epigenetic mechanisms in two patients. Finally, significant changes in the expression of immune mediators and genes related to immunological synapse, as well as differential infiltration of immune populations in the tumor environment were observed and led to a switch from a hot to a cold immune profile in three patients. Conclusions:This comprehensive analysis reveals patient-dependent spatial intra-tumor heterogeneity at the genetic, transcriptomic and epigenetic levels and in the immune landscape of the tumor microenvironment. Results support

  Study EGAS00001005328

• SmMIP-tools:a computational toolset for processing and analysis of single-molecule molecular inversion probes derived data

  Single-molecule molecular inversion probes (smMIPs) provides a modular and cost-effective platform for high-multiplex targeted next-generation sequencing (NGS). Nevertheless, translating the raw smMIP-derived sequencing data into accurate and meaningful information currently requires proficient computational skills and a large amount of computational work, prohibiting wide-scale adoption of smMIP-based technologies. To enable easy, efficient, and accurate interrogation of smMIP-derived data, we developed SmMIP-tools, a computational toolset that combines the critical analytic steps for smMIP data interpretation into a single computational pipeline. Here, we describe in detail two major components of the software. The first is a read processing tool that performs quality control steps, generates read-smMIP linkages and retrieves molecular tags. The second is an error-aware variant caller capable of detecting single nucleotide variants (SNVs) and short insertions and deletions (indels). Using a cell-line DNA dilution series and a cohort of blood cancer patients, we benchmarked SmMIP-tools and evaluated its performance against clinical sequencing reports. We anticipate that SmMIP-tools will increase accessibility to smMIP-technology, enabling cost-effective genetic research to push personalized medicine forward.

  Study EGAS00001005359

• The genetic scenario of Mercheros: an under-represented population within the Iberian Peninsula

  In this work we analyze for the first time the Merchero population, a Spanish minority ethnic group that has been scarcely studied and historically persecuted. We examine the demographic history and the genetic scenario of Mercheros by using genome-wide array data and whole mitochondrial sequences.

  Study EGAS00001005360

• Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

  Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell-type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9-binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (>=5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

  Study EGAS00001005366

• SG Peranakan Project: Genetic admixture in the culturally unique Peranakan Chinese population in Southeast Asia

  The Peranakan Chinese are culturally unique descendants of immigrants from China who settled in the Malay Archipelago ~300-500 years ago. Today, among large communities in Southeast Asia, the Peranakans have preserved Chinese traditions with strong influence from the local indigenous Malays. Yet, whether or to what extent genetic admixture co-occurred with the cultural mixture has been a topic of ongoing debate. We performed whole-genome sequencing (WGS) on 177 Singapore (SG) Peranakans and analyzed the data jointly with WGS data of Asian and European populations. We estimated that Peranakan Chinese inherited ~5.62% (95% confidence interval [CI]: 4.75-6.46%) Malay ancestry, much higher than that in SG Chinese (1.08%, 0.69-1.53%), southern Chinese (0.86%, 0.57-1.31%), and northern Chinese (0.25%, 0.18-0.33%). A sex-biased admixture history, in which the Malay ancestry was contributed primarily by females, was supported by X chromosomal variants, and mitochondrial (MT) and Y haplogroups. Finally, we identified an ancient admixture event shared by Peranakan Chinese and SG Chinese ~1,612 (95% CI: 1,345-1,923) years ago, coinciding with the settlement history of Han Chinese in southern China, apart from the recent admixture event with Malays unique to Peranakan Chinese ~190 (159-213) years ago. These findings greatly advance our understanding of the dispersal history of Chinese and their interaction with indigenous populations in Southeast Asia.

  Study EGAS00001005379

• Papua New Guinean Genome Diversity Project

  Papua New Guinean populations have one of the highest genetic diversity of the world. Our dataset provides genomic data covering most of the Papua New Guinean territory (n=58).

  Study EGAS00001005393

• Analysis of RAD51C promoter methylation using targeted bisulfite sequencing (amplicon sequencing) in ovarian cancer pre-clinical models and patient samples.

  In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in the DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved. In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages. In a cohort of 11 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic 'scarring', indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after 6 lines of therapy (4 platinum-based), whilst another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neo-adjuvant chemotherapy. As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy.

  Study EGAS00001005395

• Interferon lambda 4 impairs viral antigen presentation and attenuates T cell responses

  Genetic variants of the interferon lambda (IFNL) gene locus are strongly associated with spontaneous and IFN treatment induced clearance of hepatitis C virus (HCV) infections. Individuals with the ancestral IFNL4-dG allele are not able to clear HCV in the acute phase and develop chronic hepatitis C (CHC) with more than 90% probability. Paradoxically, the IFNL4-dG allele encodes a fully functional IFNλ4 protein with antiviral activity against HCV and coronavirus. Here we describe a hitherto unknown effect of IFNλ4 on HCV antigen presentation. Only minor amounts of IFNλ4 are secreted, because the protein is largely retained in the endoplasmic reticulum (ER) where it induces ER stress. Stressed cells were significantly weaker activators of HCV specific CD8+ T cells. This was not due to reduced MHC I surface presentation or extracellular IFNλ4 effects, since T cell responses could be restored by exogenous loading of MHC with HCV antigens. Rather, IFNλ4 induced ER stress impaired HCV antigen processing and/or loading onto the MHC I complex. Our results provide a potential explanation for the IFNλ4 – HCV paradox.

  Study EGAS00001005396

• Genomic characterization of co-existing biliary tract intraepithelial neoplasia and carcinoma lesions reveals distinct evolutionary paths of gallbladder cancer

  Gallbladder carcinoma is the most common cancer of the biliary tract with dismal survival largely due to delayed diagnosis. Biliary tract intraepithelial neoplasia (BilIN) is the common benign tumor that is suspected to be precancerous lesions. However, the genetic and evolutionary relationships between BilIN and carcinoma remain unclear. Here we performed whole-exome sequencing of coexisting low-grade BilIN (adenoma), high-grade BilIN, and carcinoma lesions, and normal tissues from the same patients.

  Study EGAS00001005402

• Germline RUNX1 Variation and Predisposition to Childhood Acute Lymphoblastic

  Somatic mutations of RUNX1, which encodes the myeloid and lymphoid transcriptional factor RUNX1, are common in both B- and T- acute lymphoid leukemia (ALL) and are associated with poor prognosis of T-ALL. However, there has been no comprehensive investigation of the pattern or prevalence of RUNX1 germline mutation in both B- and T-ALL. Here we report germline RUNX1 variants in 1.23% of B-ALL and 2.11% of T-ALL, identifying 31 unique variants in 62 B-ALL and 18 unique variants in 26 T-ALL children. The majority of frameshift and nonsense variants affected RUNX1 function in transcriptional regulation, hematopoiesis, and cellular proliferation. We identified JAK3 as the most frequent somatic mutation in T-ALL with RUNX1 variants. These results not only identify RUNX1 as a leukemia predisposition gene but also further underline the importance of germline genetic variants to the development of ALL

  Study EGAS00001005403

• Identifying transcriptional programs underlying anti-EGFR small molecule response and resistance with TraCe-seq

  Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.

  Study EGAS00001005405

• Philippine Ayta possess the highest level of Denisovan ancestry in the world

  Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ~2.3 million genotypes from 118 ethnic groups of the Philippines including 25 diverse self-identified Negrito populations, along with high-coverage genomes of AustraloPapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world – ~30-40% greater than that of Australians and Papuans – consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans, and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

  Study EGAS00001005407

• Philippine Ayta possess the highest level of Denisovan ancestry in the world

  Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ~2.3 million genotypes from 118 ethnic groups of the Philippines including 25 diverse self-identified Negrito populations, along with high-coverage genomes of AustraloPapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world – ~30-40% greater than that of Australians and Papuans – consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans, and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

  Study EGAS00001005408

• 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

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

  In this, whole exome sequencing was performed on 185 patients suffering from Azoospermia or severe oligozoospermia and on their respective fertile parents in order to investigate de novo and inherited genetic variants likely to be the cause of infertility in these men. After bioinformatic analysis and comprehensive variant interpretation several de novo point mutations and copy number variations were identified as well as multiple homozygous and compound heterozygous mutations in recessive genes that are likely influencing the infertility phenotypes observed in the patients.

  Study EGAS00001005417

• Feasibility of Ultra-Rapid Exome Sequencing in Critically Ill Infants and Children With Suspected Monogenic Conditions in the Australian Public Health Care System

  Critically ill infants and children with suspected monogenic conditions underwent ultra-rapid whole exome genetic testing. A molecular diagnosis was established in 51% of the patients. This study suggests feasibility of ultra-rapid genomic testing in critically ill pediatric patients.

  Study EGAS00001005430

• Spatial concordance of DNA methylation classification in diffuse glioma

  Background: Intratumoral heterogeneity is a hallmark of diffuse gliomas. DNA methylation profiling is an emerging approach in the clinical classification of brain tumors. The goal of this study is to investigate the effects of intratumoral heterogeneity on classification confidence. Methods: We used neuronavigation to acquire 133 image-guided and spatially-separated stereotactic biopsy samples from 16 adult patients with a diffuse glioma (7 IDH-wildtype and 2 IDH-mutant glioblastoma, 6 diffuse astrocytoma, IDH-mutant and 1 oligodendroglioma, IDH-mutant and 1p19q codeleted), which we characterized using DNA methylation arrays. Samples were obtained from regions with and without abnormalities on contrast enhanced T1 weighted and fluid-attenuated inversion recovery MRI. Methylation profiles were analyzed to devise a three-dimensional reconstruction of (epi)genetic heterogeneity. Tumor purity was assessed from clonal methylation sites. Results: Molecular aberrations indicated that tumor was found outside imaging abnormalities, underlining the infiltrative nature of this tumor and the limitations of current routine imaging modalities. We demonstrate that tumor purity is highly variable between samples and explains a substantial part of apparent epigenetic spatial heterogeneity. We observed that DNA methylation subtypes are often, but not always, conserved in space taking tumor purity and prediction accuracy into account. Conclusion: Our results underscore the infiltrative nature of diffuse gliomas and suggest that DNA methylation subtypes are relatively concordant in this tumor type, although some heterogeneity exists.

  Study EGAS00001005434

• Whole exome sequencing (bam files) of 55 samples of myxofibrosarcoma and 44 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the whole exome sequence was performed.

  Study EGAS00001005442

• Targeted-capture sequencing (bam files) of 81 samples of myxofibrosarcoma and 44 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the Targeted-capture sequencing was performed.

  Study EGAS00001005443

• Whole genome sequencing (bam files) of 5 samples of myxofibrosarcoma and 5 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the whole genome sequence was performed.

  Study EGAS00001005444

• MutWP4__CRUK_Grand_Challenge_Mutographs_of_Cancer__organoid_control_WGS

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, Kings College London will characterise the mutational signatures induced by putative human carcinogens in order to identify the origins of mutational signatures found in human cancers. To achieve this human organoid cell cultures will be exposed to a representative catalogue of known or suspected human carcinogens and mutagens and, using whole genome sequencing, the patterns of mutations induced by them will be determined. Somatic mutational signatures will be subsequently extracted by non-negative matrix factorisation methods and correlated with exposure data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001005449

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__head___neck

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute. Non-tumour DNA is also studied for some cancer types. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001005450

• MutWP1__CRUK_Grand_Challenge__normal_kidney_Nanoseq

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute. Non-tumour DNA is also studied for some cancer types. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001005451

• Next Generation Sequencing platform for targeted Personalized Therapy of Leukemia - Acute myeloid leukemia

  The advent of high-throughput next generation sequencing (NGS) technologies that are revolutionizing genomics and transcriptomics by providing a single base resolution tool for a unified deep analysis of diseases complexity allows a fast and cost-efficient fine-scale assessment of the genetic variability hidden within cohorts of patients affected by the same leukemia. That being so, by potentially highlighting inter-individual differences that may play a role in the differential success of diverse therapeutic interventions, they promise to be crucial for selecting the most appropriate medical treatments. The project aims at the identification of novel prognostic biomarkers for acute myeloid leukemia (AML) and studied the molecular differences between aneuploid and euploid AML.

  Study EGAS00001005464

• An epigenomics time course analysis of covid19 patients from Quebec, Canada

  Investigating the cellular and epigenetic changes at play in severely ill hospitalized COVID-19 patients at admission and how it evolves with prognosis. Single-cell RNA-seq analysis done on longitudinal samples highlighted changes in the monocytic compartment, which we explored using whole-genome epigenomics methods (ATACseq abd WGBS).

  Study EGAS00001005468

• 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

• Role of cohesin/CTCF in human monocyte differentiation

  Here we study genome architecture in a postmitotic differentiation setting, the differentiation of human blood monocytes (MO). We profiled and compared epigenetic, transcriptome and 3D conformation landscapes during MO differentiation (either into dendritic cells or macrophages) across the genome, with and without siRNA-mediated depletion of RD21 or CTCF, and detected numerous architectural changes, ranging from higher level compartments

  Study EGAS00001005508

• Multi-omics integration reveals only minor long-term molecular and functional sequelae in immune cells of individuals recovered from COVID-19

  The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in convalescent COVID-19 and healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays moderate changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1RA and IFN?? was higher in convalescent individuals. In conclusion, despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

  Study EGAS00001005529

• Pancreatic islets PISA RNA-seq samples

  GWAS have identified more than 700 genetic signals associated with type 2 diabetes (T2D). To gain insight into the underlying molecular mechanisms, we created the Translational human pancreatic Islet Genotype tissue-Expression Resource (TIGER), aggregating >500 human islet RNA-seq and genotyping datasets. We imputed genotypes using 4 reference panels and meta-analyzed cohorts to improve coverage of expression quantitative trait loci (eQTL) and developed a method to combine allele-specific expression across samples (cASE). We identified >1 million islet eQTLs (56% novel), of which 53 colocalize with T2D signals (60% novel). Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identified 8 novel cASE colocalizations, among which an SLC30A8 T2D associated variant. We make all the data available through the open-access TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translate this into therapeutic insight and precision medicine for T2D.

  Study EGAS00001005535

• Germline variants in patients with rare cancers and their implications for precision cancer medicine: experiences from the Multicenter MASTER Trial by the German Cancer Consortium (HIPO_021)

  The era of precision oncology entails an increasing amount of molecular data from tumor sequencing in research and diagnostics. In this setting, germline variant analysis may add clinical benefits for patients and family members. In this study, we evaluated germline variants from genome/exome matched tumor-blood sequencing in 1.485 advanced cancer patients within a prospective observational study carried out by the German Cancer Consortium. Included patients had cancers of rare entities (79.2%) or were under the age of 51 at first diagnosis (76.6%). Ten percent of the patients harbored pathogenic germline variants in 35 tumor predisposition genes with autosomal-dominant inheritance such as BRCA1, BRCA2, TP53, RB1, CHEK2, ATM, PALB2, SDHB etc., while 4.6% of the patients were carriers of pathogenic germline variants in 18 genes with autosomal-recessive inheritance of tumor predisposition. Based on pathogenic germline variants, 14% of the patients received genetic counselling recommendations in the molecular tumor board report. Forty five percent of pathogenic germline variants supported a molecular treatment recommendation, which was implemented for 28 patients. Furthermore, we compared our approach and results to previous studies, broadening the insight into pan cancer germline findings but also revealing differences between studies on various levels, such as cohort demographics, tumor entities, gene nominations, variant detection and variant interpretation. Our data support the importance of germline variant interpretation in precision oncology studies for patients and their families, if applicable. However, guidelines for evaluation and reporting of germline variants in the context of genetic tumor profiling need to be established and further research on the relevance of germline mutations as biomarkers for treatment recommendations is needed.

  Study EGAS00001005537

• Investigation of the keratinocytic gene expression pattern in Hidradenitis suppurativa

  The underlying pathogenetic factors generating the innate immune signal necessary for T cell activation, initiation and chronification of Hidradenitis suppurativa (HS, also known as Acne inversa) are still poorly understood. Emerging evidence suggests that defective keratinocyte function critically contributes to HS disease development and progression. To elucidate the role of keratinocytes in HS lesion formation, we compared the transcriptomes of isolated lesional and perilesional HS epidermis by RNA sequencing.

  Study EGAS00001005544

• Expression-Based Subtypes Define Pathologic Response to Neoadjuvant Immune-Checkpoint Inhibitors in Muscle-Invasive Bladder Cancer

  Checkpoint immunotherapy (CPI) has increased survival for some patients with advanced stage bladder cancer (BCa). However, most patients do not respond. Here, we characterized the tumor and immune microenvironment in pre- and post-treatment tumors from the PURE01 neoadjuvant pembrolizumab immunotherapy trial, using a consolidative approach that combined transcriptional and genetic profiling with digital spatial profiling.

  Study EGAS00001005549

• Genomic characterization of retinoblastoma (targeted sequencing)

  Retinoblastoma is the most frequent intraocular malignancy in children, originating from a maturing cone precursor in the developing retina. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using multi-omics data, we demonstrate the existence of two retinoblastoma subtypes. Subtype 1, of earlier onset, includes most of the heritable forms. It harbors few genetic alterations other than the initiating RB1 inactivation and corresponds to differentiated tumors expressing mature cone markers. By contrast, subtype 2 tumors harbor frequent recurrent genetic alterations including MYCN-amplification. They express markers of less differentiated cone together with neuronal/ganglion cell markers with marked inter- and intra-tumor heterogeneity. The cone dedifferentiation in subtype 2 is associated with stemness features including low immune and interferon response, E2F and MYC/MYCN activation and a higher propensity for metastasis. The recognition of these two subtypes, one maintaining a cone-differentiated state, and the other, more aggressive, associated with cone dedifferentiation and expression of neuronal markers, opens up important biological and clinical perspectives for retinoblastomas. In this study, we performed targeted sequencing of the exonic regions of RB1, BCOR and ARID1A on 37 samples from the series of 102 retinoblastomas. Among them, 23 samples were not subjected to whole-exome sequencing, and 14 samples have the whole-exome sequencing data (deposited in EGAS00001005248).

  Study EGAS00001005550