22929 results for "ega"

in 54.07 milliseconds.

• The evolution of hematopoietic cells under cancer therapy

  Chemotherapies may increase mutagenesis of healthy cells and change the selective pressures in tissues thus influencing their evolution. However, their contributions to the mutation burden and clonal expansions of healthy somatic tissues are not clear. Here, exploiting the mutational footprint of some chemotherapies, we explore their influence on the evolution of hematopoietic cells. Cells of Acute Myeloid Leukemia (AML) secondary to treatment with platinum-based drugs showed the mutational footprint of these drugs, indicating that non-malignant blood cells received chemotherapy mutations. No trace of 5-fluorouracil (5FU) mutational signature was found in AML secondary to exposure to 5FU, suggesting that cells establishing the AML were quiescent during treatment. Using the platinum-based mutational signature as a barcode, we determined that secondary AMLs began their clonal expansion after the start of the cytotoxic treatment. Its absence in clonal hematopoiesis cases is consistent with a clonal expansion that predates the exposure to the cytotoxic agent.

  Study EGAS00001005234

• Myoepithelial progenitors as founder cells of hyperplastic human breast lesions upon PIK3CA transformation

  The myoepithelial (MEP) lineage of human breast comprises bipotent and multipotent progenitors in ducts and terminal duct lobular units (TDLUs). We here assess whether this heterogeneity impacts on oncogenic PIK3CA transformation. Single cell RNA sequencing (scRNA-seq) and multicolor imaging reveal that terminal ducts represent the most enriched source of cells with ductal MEP markers including a-smooth muscle actin (a-SMA), keratin K14, K17 and CD200. Furthermore, we find neighboring CD200high and CD200low progenitors within terminal ducts. When sorted and kept in ground state conditions, their CD200low and CD200high phenotypes are preserved. Upon differentiation, progenitors remain multipotent and bipotent, respectively. Immortalized progenitors are transduced with shp53 and mutant PIK3CA. Upon transplantation, CD200low MEP progenitors distinguish from CD200high by the formation of multilayered structures with a hyperplastic inner layer of luminal epithelial cells. We suggest a model with spatially distributed MEP progenitors as founder cells of biphasic breast lesions with implications for early detection and prevention strategies.

  Study EGAS00001005933

• Gain of function DNMT3A mutations cause microcephalic dwarfism and hypermethylation of Polycomb-regulated regions

  DNA methylation and Polycomb are key factors in the establishment of vertebrate cellular identity and fate. Here we report de novo missense mutations in DNMT3A, encoding the DNA methyltransferase DNMT3A, that cause microcephalic dwarfism, a hypocellular disorder of extreme global growth failure. Substitutions in the PWWP domain abrogate binding to the histone modifications H3K36me2/3, and alter DNA methylation in patient cells. Polycomb-associated DNA methylation canyons/valleys, hypomethylated domains encompassing developmental genes, become methylated with concomitant depletion of H3K27me3 and H3K4me3 bivalent marks. Such de novo DNA methylation occurs during differentiation of Dnmt3aW326R pluripotent cells in vitro, and is also evident in Dnmt3aW326R/+ dwarf mice. We therefore propose that the interaction of the DNMT3A PWWP domain with H3K36me2/3 normally limits DNA methylation of polycomb-marked regions. Our findings implicate the interplay between DNA methylation and polycomb at key developmental regulators as a determinant of organism size in mammals.

  Study EGAS00001003232

• Fibroblast heterogeneity drives metastatic spread in breast cancer through distinct mechanisms

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

  Study EGAS00001003238

• Homologous recombination DNA repair deficiency and PARP inhibition activity in primary triple negative breast cancer

  Triple negative breast cancer (TNBC) encompasses molecularly different subgroups, with a subgroup harboring evidence of defective homologous recombination (HR) DNA repair. Within a phase 2 window clinical trial, RIO trial (EudraCT 2014-003319-12), we investigate the activity of PARP inhibitors in 43 patients with untreated TNBC. The primary endpoint, decreased Ki67, occurs in 12% of TNBC. Multiple different approaches identify HR deficiency in the majority of TNBC, 69% with the mutational-signature-based HRDetect assay. Cancers with mutational signatures of HR deficiency have a functional defect in HR. Following rucaparib treatment there is no association of Ki67 change with HR deficiency. In contrast, early circulating tumor DNA dynamics identify activity of rucaparib in cancers with mutation-signatures of HR deficiency, with rucaparib inducing expression of interferon response genes in HR deficient cancers. The majority of TNBCs have a defect in DNA repair, identifiable by mutational signature analysis, that may be targetable with PARP inhibitors.

  Study EGAS00001004405

• Genome-wide cell-free DNA mutational integration enables ultra-sensitive cancer monitoring

  In many areas of oncology, we lack sensitive tools to track low burden disease. While cell-free DNA (cfDNA) shows promise in detecting cancer mutations, we found that the combination of low tumor fraction (TF) and limited number of DNA fragments restricts low disease burden monitoring through the prevailing deep targeted sequencing paradigm. We reasoned that breadth may supplant depth of sequencing to overcome the barrier of cfDNA abundance. Whole genome sequencing (WGS) of cfDNA allowed ultra-sensitive detection, capitalizing on the cumulative signal of thousands of somatic mutations observed in solid malignancies, with TF detection sensitivity as low as 10-5. The WGS approach enabled dynamic tumor burden tracking and post-operative residual disease detection, associated with adverse outcome. Thus, we present an orthogonal framework for cfDNA cancer monitoring via genome wide mutational integration, enabling ultra-sensitive detection, overcoming the limitation of cfDNA abundance, and empowering treatment optimization in low-disease burden oncology care

  Study EGAS00001004406

• SARS‐CoV‐2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells

  The SARS‐CoV‐2 pandemic affecting the human respiratory system severely challenges public health and urgently demands for increasing our understanding of COVID‐19 pathogenesis, especially host factors facilitating virus infection and replication. SARS‐CoV‐2 was reported to enter cells via binding to ACE2, followed by its priming by TMPRSS2. Here, we investigate ACE2 and TMPRSS2 expression levels and their distribution across cell types in lung tissue (twelve donors, 39,778 cells) and in cells derived from subsegmental bronchial branches (four donors, 17,521 cells) by single nuclei and single cell RNA sequencing, respectively. While TMPRSS2 is strongly expressed in both tissues, in the subsegmental bronchial branches ACE2 is predominantly expressed in a transient secretory cell type. Interestingly, these transiently differentiating cells show an enrichment for pathways related to RHO GTPase function and viral processes suggesting increased vulnerability for SARS‐CoV‐2 infection. Our data provide a rich resource for future investigations of COVID‐19 infection and pathogenesis.

  Study EGAS00001004419

• Somatic_mutation_profiling_of_intestinal_crypts_from_IBD_Full_STDY

  Following on from our pilot study we will investigate further the role that somatic mutations in intestinal stem cells play in IBD. We will use laser capture microdissection to ascertain epithelial cells from an individual crypts/vilus. We will extract DNA from these and use a new low-input material library prep methods to make accurate libraries for sequencing. We will perform 15X whole-genome sequencing of crypts biopsied from patients. Biopsies have been taken from inflammed and noninflammed patches to compare mutation rates and mutational signatures. We will also test if crypts from inflammed patches of a single patient are clonal (either within a given inflammed patch or across inflammed patches). We will quantify the rate of cancer driver mutations in inflammed and non-inflammed crypts. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001003249

• Custom long non-coding RNA capture enhances detection sensitivity in different human sample types.

  Long non-coding RNAs (lncRNAs) are a heterogeneous group of transcripts that lack protein coding potential and display regulatory functions in various cellular processes. As a result of their cell- and cancer-specific expression patterns, lncRNAs have emerged as potential diagnostic and therapeutic targets. The accurate characterization of lncRNAs in bulk transcriptome data remains challenging due to their low abundance compared to protein coding genes. To tackle this issue, we describe a unique short-read custom lncRNA capture sequencing approach that relies on a comprehensive set of 565,878 capture probes for 49,372 human lncRNA genes. This custom lncRNA capture approach was evaluated on various sample types ranging from artificial high-quality RNA mixtures to more challenging formalin-fixed paraffin-embedded tissue and biofluid material. The custom enrichment approach allows the detection of a more diverse repertoire of lncRNAs, with better reproducibility and higher coverage compared to classic total RNA-sequencing.

  Study EGAS00001005418

• Detection of cell-free DNA fragmentation and copy number alterations in cerebrospinal fluid from glioma patients

  Glioma is difficult to detect or characterise using current liquid biopsy approaches. Detection of cell-free tumor DNA (cftDNA) in cerebrospinal fluid (CSF) has been proposed as an alternative to detection in plasma. We used shallow whole-genome sequencing (sWGS, at a coverage of < 0.4x) of cell-free DNA from the CSF of 13 patients with primary glioma to determine somatic copy number alterations and DNA fragmentation patterns. This allowed us to determine the presence of cftDNA in CSF without any prior knowledge of point mutations present in the tumor. We also showed that the fragmentation pattern of cell-free DNA in CSF is different from that in plasma. This low-cost screening method provides information on the tumor genome and can be used to target those patients with high levels of cftDNA for further larger-scale sequencing, such as by whole-exome and whole-genome sequencing.

  Study EGAS00001003255

• CCL22 chemokine mutations drive natural killer cell lymphoproliferative disease by biasing GPCR signaling

  Chronic lymphoproliferative disorder of natural killer cells (CLPD-NK) is characterized by clonal expansion of NK cells of incompletely understood genomic basis. Here, we report somatic mutations in the chemokine gene CCL22 as the hallmark of a distinct subset of CLPD-NK. CCL22 mutations were enriched at highly conserved residues, were mutually exclusive of STAT3 mutations, and were associated with gene expression similar to normal CD16dim/CD56bright NK cells. Mechanistically, the mutations resulted in ligand-biased chemokine receptor signaling, with decreased internalization of the G-protein-coupled receptor for CCL22, CCR4, via impaired β-arrestin recruitment. This resulted in increased cell chemotaxis in vitro, bidirectional crosstalk with the hematopoietic microenvironment, and enhanced NK cell proliferation in vivo in transgenic human IL-15 mice. Somatic CCL22 mutations illustrate a novel mechanism of tumor formation in which gain-of-function chemokine mutations promote tumorigenesis by biased G-protein-coupled receptor signaling and dysregulation of microenvironmental crosstalk.

  Study EGAS00001006009

• ENCORAFENIB COMBINED WITH BINIMETINIB FOR BRAFV600E-MUTATED RELAPSED/REFRACTORY MULTIPLE MYELOMA: THE PHASE II GMMG-BIRMA TRIAL (Hipo_K08K)

  Activating BRAF mutations are rare events in multiple myeloma but have been shown to be promising therapeutic targets in small case studies. This multicenter phase II trial assessed the efficacy and safety of the BRAF/MEK inhibitors, encorafenib and binimetinib, in patients with relapsed/refractory multiple myeloma (RRMM) carrying an activating BRAFV600E mutation. Twelve RRMM patients with a median of five prior lines of therapy were enrolled. The study met its primary endpoint with 10 patients achieving at least a partial response according to IMWG criteria (overall response rate 83.3%). Responses occurred rapidly within the first cycle and deepened over time. The median progression-free survival (PFS) was 5.6 months and overall survival was 55% at 24 months. Genomic profiling revealed RAS mutations and structural variants involving the BRAF locus to be drivers of resistance to BRAF/MEK inhibition in RRMM. This trial demonstrates the value of translationally driven clinical trials in selected patient populations.

  Study EGAS00001005973

• Germline TP53 mutations undergo copy number gain years prior to tumor diagnosis

  Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 mutations. We performed whole-genome sequence (WGS) analysis of tumors from 22 TP53 mutation carriers. We observed somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identified near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurred earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, revealed that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising very early in life or in utero.years prior to tumor diagnosis.

  Study EGAS00001005982

• Multimodal immunogenomic biomarker analysis of tumors from pediatric patients enrolled to a phase 1-2 study of single-agent atezolizumab

  We report herein an extensive exploratory biomarker analysis of refractory tumors taken from pediatric patients prior to receiving atezolizumab monotherapy in the phase 1-2 iMATRIX-atezolizumab trial (NCT02541604). A high percentage of CD8+ T cells and elevated protein levels of programmed cell death ligand 1 (PD-L1) were associated with progression-free survival (PFS). T-cell receptor (TCR) sequencing revealed that diverse infiltrating TCR repertoire at baseline was prognostic. We found no associations between panel-based tumor mutation burden (TMB) or specific genetic aberrations with PFS in this study. Through a pan-cancer gene co-expression network analysis, we developed a novel tumor-agnostic Pediatric Cytotoxicity and Antigen Presentation (PedCAP) signature that was associated with improved PFS in the iMATRIX-atezo study. Our study highlights features of immune response in pediatric cancers when treated with immune checkpoint inhibitors and provides a multi-biomarker pediatric immunogram framework to guide prospective clinical trials in pediatric cancers.

  Study EGAS00001006004

• A Proteogenomic Analysis of Clear Cell Renal Cell Carcinoma in a Chinese Population

  Clear cell renal cell carcinoma (ccRCC) is a common and aggressive subtype of renal cancer. Here we conduct a comprehensive proteogenomic analysis of 232 tumor and adjacent non-tumor tissue pairs from Chinese ccRCC patients. By comparing with tumor adjacent tissues, we find that ccRCC shows extensive metabolic dysregulation and an enhanced immune response. Molecular subtyping classifies ccRCC tumors into three subtypes (GP1–3), among which the most aggressive GP1 exhibits the strongest immune phenotype, increased metastasis, and metabolic imbalance, linking the multi-omics-derived phenotypes to clinical outcomes of ccRCC. Nicotinamide N-methyltransferase (NNMT), a one-carbon 38 metabolic enzyme, is identified as a potential marker of ccRCC and a drug target for GP1. We demonstrate that NNMT induces DNA-dependent protein kinase catalytic subunit (DNA-PKcs) homocysteinylation, increases DNA repair, and promotes ccRCC tumor growth. This study provides insights into the biological underpinnings and prognosis assessment of ccRCC, revealing targetable metabolic vulnerabilities.

  Study EGAS00001006005

• HSC_colony_many_years_post_allogeneic_bone_marrow_transplant___Bait_set_2

  This study looks at human donors and recipients of allogeneic blood & marrow marrow transplantation. Each transplant "pair" (donor & recipient) underwent the transplant procedure 10 -30 years prior to sampling - the majority for an underlying diagnosis of leukemia. We will perform whole-genome sequencing on tens to hundreds of single HSPC (haematopoietic stem and progenitor cell)-derived colonies per individual. We will infer the phylogeny of these cells (those from both donor & recipient within a pair) by analysing the somatically acquired mutations within each cell. The structure of the phylogeny will allow estimates of hitherto unknown parameters such as the number of HSCs that engraft long-term, the dynamics of the engraftment process on the HSC level, and the additional mutation burden acquired through the arduous transplant procedure. This TGS part of the study is used for assessing the contribution of clones to mature blood subsets, and better timing mutations relative to transplantation.

  Study EGAS00001005534

• Oncogenic cooperation in a human de novo T-ALL model

  We previously modelled T-cell acute lymphoblastic leukemia (T-ALL) by over-expressing four oncogenes (activated NOTCH1, LMO2, TAL1, and BMI1) in healthy human umbilical cord blood cells through lentiviral transduction (Kusakabe et al, Nat Commun 10, 2913 (2019). https://doi.org/10.1038/s41467-019-10510-8). We have now tested individual effects of the four oncogenes by transducing human cells with all possible combinations of one, two, three, and four oncogenes. We also tested the effects of additional known T-ALL oncogenes including HOXA9, LYL1, MEF2C, NXK2.1, TLX1, and TLX3. In addition to culture in vitro on OP9-DL1 feeders, we also injected transduced cells into immunodeficient NSG mice. We observed differing abilities of the various oncogene combinations to support growth in vitro and frank leukemia in vivo. This study includes RNA-seq data from samples of transduced cells cultured in vitro as well as cells harvested from clinically moribund mice after development of leukemia in vivo.

  Study EGAS00001006055

• Molecular and functional profiling of plasmablastic lymphoma

  Plasmablastic lymphoma (PBL) represents a rare and aggressive lymphoma subtype frequently associated with immunosuppression. Clinically, patients with PBL are characterized by poor outcome. The current understanding of the molecular pathogenesis is limited. A hallmark of PBL represents its plasmacytic differentiation with loss of B-cell markers and, in 60% of cases, its association with Epstein-Barr virus (EBV). Roughly 50% of PBLs harbor a MYC translocation. Here, we provide a comprehensive integrated genomic analysis using whole exome sequencing (WES) and genome-wide copy number determination in a large cohort of 96 primary PBL samples. We identify alterations activating the RAS-RAF, JAK-STAT, and NOTCH pathways as well as frequent high-level amplifications in MCL1 and IRF4. The functional impact of these alterations is assessed using an unbiased shRNA screen in a PBL model. These analyses identify the IRF4 and JAK-STAT pathways as promising molecular targets to improve outcome of PBL patients.

  Study EGAS00001004659

• Enhanced cortical neural stem cell identity through SMAD/WNT inhibition in human cerebral organoids facilitates emergence of outer radial glial cells

  Cerebral organoids exhibit broad regional heterogeneity accompanied by limited cortical cellular diversity despite the tremendous upsurge in derivation methods, suggesting inadequate patterning of early neural stem cells (NSCs). Here we show that a short and early dual-SMAD/WNT inhibition course is necessary and sufficient for establishing robust and lasting cortical organoid NSC identity, efficiently suppressing of non-cortical NSC fates, while other widely used methods are inconsistent in their cortical NSC specification capacity. Accordingly, this method selectively enriches for outer radial glia (oRG) NSCs, which cytoarchitecturally demarcate well-defined outer sub-ventricular (oSVZ)-like regions propagating from superiorly radially organized, apical cortical rosette NSCs. Finally, this method culminates in the emergence of molecularly distinct deep and upper cortical layer neurons, and reliably uncovers cortex-specific microcephaly defects. Thus, a short SMAD/WNT inhibition is critical for establishing a rich cortical cell repertoire that enables mirroring fundamental molecular and cytoarchitectural features of cortical development and meaningful disease modeling.

  Study EGAS00001006063

• A comprehensive characterization of the cell-free transcriptome reveals tissue- and subtype-specific biomarkers for cancer detection

  Cell-free RNA (cfRNA) is a promising analyte for cancer detection. However, a comprehensive assessment of cfRNA in individuals with and without cancer has not been conducted. We performed the first transcriptome-wide characterization of cfRNA in cancer (stage III breast [n=46], lung [n=30]) and non-cancer (n=89) participants from the Circulating Cell-free Genome Atlas (NCT02889978). Of 57,820 annotated genes, 39,564 (68%) were not detected in cfRNA from non-cancer individuals. Within these low-noise regions, we identified tissue- and cancer-specific genes, defined as "dark channel biomarker" (DCB) genes, that were recurrently detected in individuals with cancer. DCB levels in plasma were correlated with tumor shedding rate and RNA expression in matched tissue, suggesting that DCBs with high expression in tumor tissue could enhance cancer detection in patients with low levels of circulating tumor DNA. Overall, cfRNA provides a unique opportunity to detect cancer, predict the tumor tissue of origin, and determine the cancer subtype.

  Study EGAS00001004704

• 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

• single cell RNA-seq and ATAC-seq of human fetal forebrain tissue

  Oligodendrogenesis in the human central nervous system has been mainly observed at the second trimester of gestation, a much later developmental stage compared to mouse. Here we characterize the transcriptomic neural diversity in the human forebrain at post conceptual weeks (PCW) 8 to 10, using single-cell RNA-Seq. We find evidence of the emergence of a first wave of oligodendrocyte lineage cells as early as PCW 8, which we also confirm at the epigenomic level with single-cell ATAC-Seq. Using regulatory network inference, we predict key transcriptional events leading to the specification of oligodendrocyte precursor cells (OPCs). Moreover, by profiling the spatial expression of fifty key genes using In Situ Sequencing (ISS), we identify regions in the human ventral fetal forebrain where oligodendrogenesis first occurs. Our results indicate evolutionary conservation of the first wave of oligodendrogenesis between mouse and human and describe regulatory mechanisms involved in human OPC specification.

  Study EGAS00001006136

• Dissecting features of epigenetic variants underlying cardiometabolic risk using full-resolution epigenome profiling in regulatory elements

  Sparse profiling of CpG methylation in blood by microarrays have identified epigenetic links to common diseases. We apply methylC-capture sequencing (MCC-Seq) in a clinical population of ~200 adipose tissue and matched blood samples (Ntotal ~400), providing high-resolution methylation profiling (>1.3M CpGs) at regulatory elements. We link methylation to cardiometabolic risk through associations to circulating plasma lipid levels and identify lipid-associated CpGs with unique localization patterns in regulatory elements. We show distinct features of tissue-specific versus tissue-independent lipid-linked regulatory regions by contrasting with parallel assessments in ~800 independent adipose tissue and blood samples from the general population. We follow-up on adipose-specific regulatory regions under (1) genetic and (2) epigenetic (environmental) regulation via integrational studies. Overall, the comprehensive sequencing of regulatory element methylomes reveals a rich landscape of functional variants linked genetically as well as epigenetically to plasma lipid traits.

  Study EGAS00001003415

• Comprehensive genomic and transcriptomic analysis of three synchronous primary tumours and a recurrence from a head and neck cancer patient

  Synchronous primary malignancies occur in a small proportion of head and neck squamous cell carcinoma patients. Here, we analysed three synchronous primaries and a recurrence from one patient, by comparing the genomic and transcriptomic profiles among the tumour samples and determining the recurrence origin. We found remarkable levels of heterogeneity among the primary tumours, and through the patterns of shared mutations, we traced the origin of the recurrence. Interestingly, the patient carried germline variants which might have predisposed him to carcinogenesis, together with a history of alcohol and tobacco consumption. Mutational signature analysis confirmed the impact of alcohol exposure, with Signature 16 present in all tumour samples. Characterisation of immune cell infiltration highlighted an immunosuppressive environment in all samples, which exceeded the potential activity of T cells. Studies such as the one described here have important clinical value and contribute to personalised treatment decisions for patients with synchronous primaries and matched recurrences.

  Study EGAS00001004857

• 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

• The effect of pre-analytical and physiological variables on cell-free DNA fragmentation

  Background: assays that account for the biological properties and fragmentation of cell-free DNA (cfDNA) can improve the performance of liquid biopsy. However, pre-analytic and physiological differences between individuals on fragmentomic analysis are poorly defined. Methods: we analyzed the impact of collection tube, plasma processing time and physiology on the size distribution of cfDNA, their genome-wide representation and sequence diversity at the cfDNA fragment-ends using shallow Whole Genome Sequencing. Results: we observed that using different stabilizing collection tubes, or processing times does not affect the cfDNA fragment sizes, but can impact the genome-wide fragmentation patterns and fragment-end sequences of cfDNA. In addition, beyond differences depending on the gender, the physiological conditions tested between 63 individuals (age, body mass index, use of medication and chronic conditions) minimally influenced the outcome of fragmentomic methods. Conclusions: our results highlight that fragmentomic approaches have potential for implementation in the clinic, pending clear traceability of analytical and physiological factors.

  Study EGAS00001005748

• The European MAPPYACTS trial: Precision Medicine Program in Pediatric and Adolescent Patients with Recurrent Malignancies

  MAPPYACTS (NCT02613962) is an international prospective precision medicine trial aiming to define tumor molecular profiles in pediatric patients with recurrent/refractory malignancies, in order to suggest the most adapted salvage treatment. From February 2016 to July 2020, 787 patients were included in France, Italy, Ireland and Spain. At least one genetic alteration leading to targeted treatment suggestion was identified in 436 patients (69%) with successful sequencing; 10% of these were considered "ready for routine use". Of 356 patients with follow-up beyond 12 months, 107 (30%) received one or more matched targeted therapies, 56% of them within early clinical trials, mainly in the AcSé-ESMART platform trial (NCT02813135). Overall, matched treatment resulted in a 17% objective response rate, those of patients with "ready for routine use" alterations was 38%. In patients with extra-cerebral tumors, 76% of actionable alterations detected in tumor tissue were also identified in circulating cell free DNA (cfDNA).

  Study EGAS00001005935

• Exome sequencing of synchronous colorectal cancers

  Approximately 4% of colorectal cancer (CRC) patients have at least two simultaneous cancers in the colon. Due to the shared environment, these synchronous CRCs (SCRCs) provide a unique setting to study colorectal carcinogenesis. Understanding whether these tumors are genetically similar or distinct is essential when designing therapeutic approaches. We performed exome sequencing of 47 primary cancers and corresponding normal samples from 23 patients. Additionally, we carried out a comprehensive mutational signature analysis to assess whether tumors had undergone similar mutational processes and the first immune cell score analysis (IS) of SCRC to analyze the interplay between immune cell invasion and mutation profile in both lesions of an individual. The tumor pairs shared only few mutations, favoring different mutations in known CRC genes and signaling pathways, and displayed variation in their signature content. Two tumor pairs had discordant mismatch repair statuses. In majority of the pairs, IS varied between primaries. Differences were not explained by any clinicopathological variable or mutation burden.

  Study EGAS00001003474

• RNA-seq on patient-derived, stage II, CRC cell lines

  The gut microbiome is a key player in the immunomodulatory and pro-tumorigenic microenvironment in colorectal cancer (CRC). It is well established that Fusobacterium’s virulence factors are responsible for its pro-tumorigenic properties, however the role of its metabolic cross-talk with the tumor remains unexplored. Here we use in vitro, in vivo, and in silico methods to describe a novel microbe-host interaction keystone: formate. Omics data reveals molecular signatures, linking a metabolically driven CRC phenotype with Fusobacteria. We observed a metabolic shift towards increased F. nucleatum formate secretion in gut-on-chip model co-cultures with patient-derived CRC cells, along with an accelerated cancer glutamine metabolism. We show that high formate levels trigger AhR signaling, increase cancer stemness and cellular invasion. Finally, we observed an expansion of Th17 cells accompanying F. nucleatum-induced tumorigenesis. Moving beyond observational studies, we applied new experimental approaches for gaining ecosystem-level mechanistic understanding of F. nucleatum’s role in cancer pathogenesis.

  Study EGAS00001005948

• The Proteogenomic Subtypes of Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) is an aggressive blood cancer with poor prognosis. We report a comprehensive proteogenomic analysis of bone-marrow biopsies from 252 uniformly treated AML patients to elucidate the molecular pathophysiology of AML in order to inform future diagnostic and therapeutic approaches. In addition to in-depth quantitative proteomics, our analysis includes cytogenetic profiling and DNA/RNA sequencing. We identify five proteomic AML subtypes, each reflecting specific biological features spanning genomic boundaries. Two of these proteomic subtypes correlate with patient outcome, but none are exclusively associated with specific genomic aberrations. Remarkably, one subtype (Mito-AML), which is only captured in the proteome, is characterized by high expression of mitochondrial proteins and confers poor outcome, with reduced remission rate and shorter overall survival upon treatment with intensive induction chemotherapy. Functional analyses reveal that Mito-AML is metabolically wired towards stronger complex I dependent respiration and is more responsive to treatment with the BCL2 inhibitor venetoclax.

  Study EGAS00001005950

• Combined gene expression and digital pathology identifies molecular mediators of T cell exclusion and immune suppression in ovarian cancer

  Close proximity between cytotoxic T lymphocytes and tumour cells is required for effective immunotherapy. However, what determines the spatial distribution of T cells in the tumour microenvironment is not well understood. Coupling digital pathology and transcriptome analysis on a large ovarian tumour cohort, we develop a machine learning approach to molecularly classify and characterize tumour-immune phenotypes. Our study identifies two important hallmarks characterizing T cell excluded tumours: 1) loss of antigen presentation on tumour cells and 2) upregulation of TGFb and activated stroma. Furthermore, we identify TGFb as a key mediator of T cell exclusion. TGFb reduces MHC-I expression in ovarian cancer cells in vitro; TGFb also activates fibroblasts and induced extracellular matrix (ECM) production as a potential physical barrier to hinder T cell infiltration. Our findings indicate that targeting TGFb may represent a promising strategy to overcome T cell exclusion and improve clinical benefits of cancer immunotherapy.

  Study EGAS00001003487

• Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs drug target discovery in ankylosing spondylitis

  Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterised by enthesitis of the spine and sacroiliac joints. Genome-wide association studies have revealed >100 genetic associations but their functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We report differential disease signatures for specific genomic regions from individual -omic modalities and multi-omic integration, showing enrichment at genetically associated loci and evidence for altered monocyte function in AS. We link putative functional SNPs with cognate genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1. We also show how disease-specific functional genomic data can be integrated with GWAS evidence to enhance therapeutic target discovery. This study combines multiple epigenetic and transcriptional analysis with GWAS to identify cell types, candidate drug target genes and transcriptional regulation of likely pathogenic relevance.

  Study EGAS00001006233

• Human genome-wide variations in the Massim region

  The Massim, a cultural region encompassing the southeastern tip of mainland Papua New Guinea (PNG) and nearby islands, is famous for the Kula trading network, in which different valuables circulate in different directions among the islands. To explore Massim genetic history, we generated genome-wide data from across the region and found variable levels of Papuan-related (indigenous) ancestry, including a novel ancestry associated with Rossel Island, and Austronesian-related ancestry that arrived later. We find genetic evidence for different patterns of contact across PNG, the Massim, and the Bismarck and Solomon Archipelagoes for Austronesian-related vs. Papuan-related ancestry, with more geographic restriction for the latter. Moreover, Kula-practicing groups share more genetic similarity than do other groups, and this likely predates the origin of Kula, suggesting that high between-group contact facilitated the formation of Kula. Our study provides the first comprehensive genome-wide assessment of Massim inhabitants and new insights into the fascinating Kula system.

  Study EGAS00001006010

• Loss of Epigenetic Barrier is Required for Enhancer Hijacking-Mediated Oncogenic Transcription

  Structural variants (SVs) involving enhancer hijacking can disrupt chromatin topologies to cause oncogene activation in cancer genomes, yet the molecular determinants for the transcriptional output of enhancer hijacking remain largely unknown. We developed a multimodal approach to integrate genome sequencing, chromosome conformation, and sequence-based deep learning for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified candidate pathogenic SVs including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for oncogenic activation of TLX3-dependent transcriptional programs. By engineering patient-associated t(5;14) in isogenic leukemia cells, we uncovered an uncharacterized mechanism whereby DNA methylation serves as an epigenetic barrier to enhancer hijacking and loss of epigenetic barrier is a molecular determinant for the transcriptional output of pathogenic SVs. Hence, leveraging the epigenetic barriers of SV-mediated oncogenic programs may provide new opportunities to reprogram gene regulation as epigenetic therapies in human disease.

  Study EGAS00001006140

• MutWP5__CRUK_Mutographs_of_Cancer__BRCA_Carriers___WG__Novaseq_

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

  Study EGAS00001003523

• Phase 1 CX-5461 Trial (Canadian Cancer Trials Group Trial IND.231)

  CX-5461 is a G-quadruplex stabilizer that exhibits synthetic lethality in homologous recombination (HR) deficient models. In this multicentre Phase Ib trial in patients with solid tumors, 40 patients were treated across 10 dose levels (50 - 650 mg/m2) to determine the recommended phase II dose (RP2D), and evaluate safety, tolerability, pharmacokinetics. Defective HR was explored as a predictive biomarker. CX-5461 was generally well tolerated, with a RP2D of 475 mg/m2 d1, 8 and 15 every 4 weeks, and dose limiting toxicities of phototoxicity. Responses were observed in 10% of patients, primarily in patients with defective HR (ORR 14%). Reversion mutations in PALB2 and BRCA2 were detected on progression following initial response in germline carriers, confirming the underlying synthetic lethal mechanism. In vitro characterization of UV sensitization showed this toxicity is related to the CX-5461 chemotype, independent of G4 synthetic lethality. These results establish clinical proof-of-concept for this first-in class G4 stabilizer.

  Study EGAS00001006173

• MutWP5__CRUK_Mutographs_of_Cancer__Cancer_Mastectomy__Exome__Novaseq_

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

  Study EGAS00001003528

• The immunopeptidome landscape associated with T cell infiltration, inflammation and immune-editing in lung cancer

  One key barrier for improving efficacy of cancer immunotherapy remains patient stratification. While patients with CD3+CD8+ T cell inflamed ‘hot’ tumors typically show better response to immune checkpoint inhibitors, it is still unknown if the repertoire of HLA bound peptides (HLAp) presented in ‘hot’ and ‘cold’ tumors is substantially different. We surveyed 61 tumor regions and adjacent non-malignant lung tissues from eight lung cancer patients and performed deep antigen discovery combining immunopeptidomics, genomics, bulk and spatial transcriptomics and explored the heterogenic expression and presentation of tumor (neo)antigens. We associated diverse immune cell populations with the immunopeptidome in CD3+CD8+ T cell excluded, infiltrated, hot and cold tumors and found evidence for lack of immune-editing and higher presentation efficiency of tumor antigens in cold and CD3+CD8+ T cell excluded tumors. This could have implications for the choice of combination therapies tailored to the patient’s mutanome and microenvironment.

  Study EGAS00001006298

• Arriba: accurate and efficient detection of gene fusions from RNA-Seq (H021)

  The identification of gene fusions from RNA sequencing data is a routine task in cancer research and precision oncology. However, despite the availability of many computational tools, fusion detection remains challenging. Existing methods suffer from poor prediction accuracy and are computationally demanding. We developed Arriba, a novel fusion detection algorithm with high sensitivity and short runtime. When applied to a large collection of published pancreatic cancer samples (n=803), Arriba identified a variety of driver fusions, many of which affected druggable proteins, including ALK, BRAF, FGFR2, NRG1, NTRK1, NTRK3, RET, and ROS1. The fusions were significantly associated with KRAS wild-type tumors and involved proteins stimulating the MAPK signaling pathway, suggesting that they substitute for activating mutations in KRAS. In addition, we confirmed the transforming potential of two novel fusions, RRBP1-RAF1 and RASGRP1-ATP1A1, in cellular assays. These results demonstrate Arriba’s utility in both basic cancer research and clinical translation.

  Study EGAS00001003554

• Widespread hypertranscription in aggressive human cancers

  Cancers are often defined by the dysregulation of specific transcriptional programs; however, the importance of global transcriptional changes is less understood. Hypertranscription is the genome-wide increase in RNA output. Hypertranscription’s prevalence, underlying drivers and prognostic significance are undefined in primary human cancer. This is due in part to limitations of expression profiling methods, which assume equal RNA output between samples. Here, we developed a computational method to directly measure hypertranscription in 7,494 human tumors, spanning 31 cancer types. Hypertranscription is ubiquitous across cancer, especially in aggressive disease. It defines patient subgroups with worse survival, even within well-established subtypes. Our data suggest that loss of transcriptional suppression underpins the hypertranscriptional phenotype. Single-cell analysis reveals hypertranscriptional clones, which dominate transcript production regardless of their size. Finally, patients with hypertranscribed mutations have improved response to immune checkpoint therapy. Our results provide fundamental insights into gene dysregulation across human cancers and may prove useful in identifying patients that would benefit from novel therapies.

  Study EGAS00001006365

• Exome & MiSeq sequencing of individuals with Huntington's disease

  The age at onset of motor symptoms in Huntington's disease (HD) is driven by HTT CAG repeat length but modified by other genes. In this study, we used exome sequencing of 683 patients with HD with extremes of onset or phenotype relative to CAG length to identify rare variants associated with clinical effect. We discovered damaging coding variants in candidate modifier genes identified in previous genome-wide association studies associated with altered HD onset or severity. Variants in FAN1 clustered in its DNA-binding and nuclease domains and were associated predominantly with earlier-onset HD. Nuclease activities of purified variants in vitro correlated with residual age at motor onset of HD. Mutating endogenous FAN1 to a nuclease-inactive form in an induced pluripotent stem cell model of HD led to rates of CAG expansion similar to those observed with complete FAN1 knockout. Together, these data implicate FAN1 nuclease activity in slowing somatic repeat expansion and hence onset of HD.

  Study EGAS00001006383

• Resolving the immune landscape of human prostate at a single cell level in health and cancer

  The prostate gland produces prostatic fluid, high in zinc and citrate and essential for the maintenance of spermatozoa. Prostate cancer is a common condition with limited treatment efficacy in castration-resistant metastatic disease, including with immune checkpoint inhibitors. We used single-cell RNA-sequencing to perform an unbiased assessment of the cellular landscape of human prostate and identified a previously unappreciated subset of tumour-enriched androgen receptor-negative luminal epithelial cells, with increased expression of cancer-associated genes. We also found a variety of innate and adaptive immune cells in normal prostate that were transcriptionally perturbed in prostate cancer. An exception was a unique, prostate-specific, zinc transporter-expressing macrophage population (MAC-MT), that contributed to tissue zinc accumulation in homeostasis, but showed enhanced inflammatory gene expression in tumours, including T cell-recruiting chemokines. Remarkably, enrichment of the MAC-MT signature in cancer biopsies was associated with improved disease-free survival, suggesting beneficial anti-tumour functions.

  Study EGAS00001005787

• Massively parallel single-cell B-cell receptor sequencing enables rapid discovery of diverse antigen-reactive antibodies

  Using high-throughput single-cell B-cell receptor sequencing (scBCR-seq) we obtained accurately paired full-length heavy- and light-chain variable regions from thousands of individual B cells in a massively parallel fashion. We sequenced more than 250,000 B cells from rat, mouse and human repertoires to characterize their lineages and expansion. In addition, we immunized rats with chicken ovalbumin and profiled antigen-reactive B cells from lymph nodes of immunized animals. The scBCR-seq data recovered 81% (n = 56/69) of B-cell lineages identified from hybridomas generated from the same set of B cells that were subjected to scBCR-seq. Importantly, scBCR-seq identified an additional 710 candidate lineages that were not recovered as hybridomas. We synthesized, expressed and tested 93 clones from the identified lineages and found that 99% (n = 92/93) of the clones were antigen-reactive. Our results establish scBCR-seq as a powerful tool for antibody discovery.

  Study EGAS00001003663

• Molecular Profiling Reveals Unique Immune and Metabolic Features of Melanoma Brain Metastases

  There is a critical need to improve our understanding of the pathogenesis of -melanoma brain metastases (MBM). Thus, we performed RNA sequencing on 88 resected MBMs and 42 patient-matched extracranial metastases; tumors with sufficient tissue also underwent whole-exome sequencing, T-cell receptor sequencing, and IHC. MBMs demonstrated heterogeneity of immune infiltrates that correlated with prior radiation and post-craniotomy survival. Comparison with patient-matched extracranial metastases identified significant immunosuppression and enrichment of oxidative phosphorylation (OXPHOS) in MBMs. Gene-expression analysis of intracranial and subcutaneous xenografts, and a spontaneous MBM model, confirmed increased OXPHOS gene expression in MBMs, which was also detected by direct metabolite profiling and [U-13C]-glucose tracing in vivo IACS-010759, an OXPHOS inhibitor currently in early-phase clinical trials, improved survival of mice bearing MAPK inhibitor-resistant intracranial melanoma xenografts and inhibited MBM formation in the spontaneous MBM model. The results provide new insights into the pathogenesis and therapeutic resistance of MBMs.

  Study EGAS00001003672

• Renal Cell Carcinoma Tumors From S-TRAC trial (NCT00375674) Whole Transcriptome Sequencing

  We report a comprehensive genomic and transcriptomic analysis of tumor samples from 171 patients at high risk for recurrent renal cell carcinoma post nephrectomy from the S-TRAC trial (NCT00375674). We identify gene expression signatures, including STRAC11 (derived from the sunitinib-treated population). The overlap in key elements captured in these gene expression signatures, which include genes representative of the tumor stroma microenvironment, regulatory T cell, and myeloid cells, suggests they are likely to be both prognostic and predictive of the anti-angiogenic effect in the adjuvant setting. These signatures also point to the identification of potential therapeutic targets for development in adjuvant renal cell carcinoma, such as MERTK and TDO2. Finally, our findings suggest that while anti-angiogenic adjuvant therapy might be important, it may not be sufficient to prevent recurrence and other factors such as immune response and that tumor environment may be of greater importance.

  Study EGAS00001006528

• Multimodal Epigenetic Sequencing Analysis (MESA) of Cell-free DNA for Non-invasive Cancer Detection

  Multimodal characterization of cell-free DNA (cfDNA) in the blood can enable the sensitive and non-invasive detection of human cancers but remains technically challenging and costly. Here, we developed Multimodal Epigenetic Sequencing Analysis (MESA), a flexible and sensitive method of capturing and integrating multimodal epigenetic information of cfDNA using a single experimental assay, i.e., non-disruptive bisulfite-free methylation sequencing, such as Enzymatic Methyl-seq (EM-seq) and TET-assisted pyridine borane sequencing (TAPS). MESA can simultaneously infer cfDNA methylation, nucleosome occupancy, nucleosome fuzziness, and fragmentation profile for regions surrounding the promoters and polyadenylation sites (PASs). MESA’s integrated analysis of multimodal epigenetic features significantly improved the performance of cancer detection models compared to the usage of any single modality alone. MESA captures additional and highly complementary epigenetic information from cfDNA without additional experimental assays, highlighting the importance and clinical prospect of using multimodal epigenetic features for non-invasive cancer detection

  Study EGAS00001006462

• HSC_colony_many_years_post_allogeneic_bone_marrow_transplant___Bait_set_3

  This study looks at human donors and recipients of allogeneic blood & marrow marrow transplantation. Each transplant "pair" (donor & recipient) underwent the transplant procedure 10 -30 years prior to sampling - the majority for an underlying diagnosis of leukemia. We will perform whole-genome sequencing on tens to hundreds of single HSPC (haematopoietic stem and progenitor cell)-derived colonies per individual. We will infer the phylogeny of these cells (those from both donor & recipient within a pair) by analysing the somatically acquired mutations within each cell. The structure of the phylogeny will allow estimates of hitherto unknown parameters such as the number of HSCs that engraft long-term, the dynamics of the engraftment process on the HSC level, and the additional mutation burden acquired through the arduous transplant procedure. This TGS part of the study is used for assessing the contribution of clones to mature blood subsets, and better timing mutations relative to transplantation.

  Study EGAS00001006572

• Whole Genome Sequencing on OCIAML-22

  Donor_OCIAML22_bulk: The donor (unsorted) that was used to generate OCI-AML22. OCI-AML22_bulk_1month_inVitro: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then sequenced. OCI-AML22_graft_CD34minus_WGS: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then injected in NSG-SGM3 mice. 21 weeks after injection, cells from 5 mice injected with 5 mice injected with 1475 cells per mice and 4 mice injected with 469 cells per mice were sorted. The human CD34- fraction was sequenced (CD45h+CD45m-7aad-annexinV-CD34-) OCI-AML22_graft_CD34positive_WGS: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then injected in NSG-SGM3 mice. 21 weeks after injection, cells from 5 mice injected with 5 mice injected with 1475 cells per mice and 4 mice injected with 469 cells per mice were sorted. The human CD34+ fraction was sequenced (CD45h+CD45m-7aad-annexinV-CD34+)

  Study EGAS00001006513

• Renal Cell Carcinoma Tumors From S-TRAC trial (NCT00375674) Exome Sequencing

  We report a comprehensive genomic and transcriptomic analysis of tumor samples from 171 patients at high risk for recurrent renal cell carcinoma post nephrectomy from the S-TRAC trial (NCT00375674). We identify gene expression signatures, including STRAC11 (derived from the sunitinib-treated population). The overlap in key elements captured in these gene expression signatures, which include genes representative of the tumor stroma microenvironment, regulatory T cell, and myeloid cells, suggests they are likely to be both prognostic and predictive of the anti-angiogenic effect in the adjuvant setting. These signatures also point to the identification of potential therapeutic targets for development in adjuvant renal cell carcinoma, such as MERTK and TDO2. Finally, our findings suggest that while anti-angiogenic adjuvant therapy might be important, it may not be sufficient to prevent recurrence and other factors such as immune response and that tumor environment may be of greater importance.

  Study EGAS00001006529

• Adult human kidney organoids originate from CD24+ cells and represent an advanced model for adult polycystic kidney disease

  Adult kidney organoids have been described as strictly tubular epithelial and termed tubuloids. While the cellular origin of tubuloids has remained elusive, here we report that they originate from a distinct CD24+ epithelial subpopulation. Long-term-cultured CD24+-derived tubuloids represent a functional human kidney tubule. We show that kidney tubuloids can be used to model the most common inherited kidney disease, namely autosomal dominant polycystic kidney disease (ADPKD), reconstituting the phenotypic hallmark of this disease with cyst formation. Single-cell RNA sequencing of CRISPR/Cas9 gene-edited PKD1 and PKD2 knockout tubuloids and human ADPKD and control tissue show similarities in upregulation of disease-driving genes. Furthermore, in a proof of concept, we demonstrate that tolvaptan, the only approved drug for ADPKD, shows a significant effect on cyst size in tubuloids but no effect in a pluripotent-stem-cell-derived model. Thus, tubuloids derive from a tubular epithelial subpopulation and represent an advanced model for ADPKD disease modeling.

  Study EGAS00001006551

• Tocilizumab treatment leads to early resolution of myeloid dysfunction and lymphopenia in patients hospitalized with COVID-19

  Tocilizumab, an anti-interleukin-6 receptor inhibitor, is recommended in global treatment guidelines for patients hospitalized with severe COVID-19. Using proteomic and transcriptomic analysis, we characterized the immune profile and identified cellular and molecular pathways directly modified by tocilizumab in peripheral blood samples collected from patients enrolled in the COVACTA study, a phase 3, randomized, double-blind, placebo-controlled trial, to assess the efficacy and safety of tocilizumab in hospitalized patients with severe COVID-19 pneumonia. We identified factors predicting disease severity and clinical outcomes, including markers of inflammation, lymphopenia, myeloid dysfunction, and organ injury. Proteomic analysis confirmed a pharmacodynamic effect for tocilizumab. Transcriptomic analysis revealed that tocilizumab treatment leads to faster resolution of lymphopenia and myeloid dysfunction associated with severe COVID-19 thus supporting an underlying anti-inflammatory mechanism of action for the beneficial effects of tocilizumab in patients hospitalized with COVID-19.

  Study EGAS00001006688

• GRIDSS, PURPLE, LINX: Unscrambling the tumor genome via integrated analysis of structural variation and copy number

  We have developed a novel, integrated and comprehensive purity, ploidy, structural variant and copy number somatic analysis toolkit for whole genome sequencing data of paired tumor/normal samples. We show that the combination of using GRIDSS for somatic structural variant calling and PURPLE for somatic copy number alteration calling allows highly sensitive, precise and consistent copy number and structural variant determination, as well as providing novel insights for short structural variants and regions of complex local topology. LINX, an interpretation tool, leverages the integrated structural variant and copy number calling to cluster individual structural variants into higher order events and chains them together to predict local derivative chromosome structure. LINX classifies and extensively annotates genomic rearrangements including simple and reciprocal breaks, LINE, viral and pseudogene insertions, and complex events such as chromothripsis. LINX also comprehensively calls genic fusions including chained fusions. Finally, our toolkit provides novel visualisation methods providing insight into complex genomic rearrangements.

  Study EGAS00001004034

• The landscape of chromothripsis across adult cancer types

  Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. We analysed chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals novel insights in the clonal evolution of tumours with chromothripsis. (H021)

  Study EGAS00001004250

• The Genetic Landscape of Ocular Adnexa MALT Lymphoma Reveals Frequent Aberrations in NFAT and MEF2B Signaling Pathways

  A comprehensive constellation of somatic non-silent mutations and copy number (CN) variations in ocular adnexa marginal zone lymphoma (OAMZL) is unknown. By utilizing whole-exome sequencing in 69 tumors we define the genetic landscape of OAMZL. Mutations and CN changes in CABIN1 (30%), RHOA (26%), TBL1XR1 (22%), and CREBBP (17%) and inactivation of TNFAIP3 (26%) were among the most common aberrations. Candidate cancer driver genes cluster in the B-cell receptor (BCR), NFkB, NOTCH and NFAT signaling pathways. One of the most commonly altered genes is CABIN1, a calcineurin inhibitor acting as a negative regulator of the NFAT and MEF2B transcriptional activity. CABIN1 deletions enhance BCR-stimulated NFAT and MEF2B transcriptional activity, while CABIN1 mutations enhance only MEF2B transcriptional activity by impairing binding of mSin3a to CABIN1. Our data provide an unbiased identification of genetically altered genes that may play a role in the molecular pathogenesis of OAMZL and serve as therapeutic targets.

  Study EGAS00001006631

• Multi-omics identify falling LRRC15 as a COVID-19 severity marker and persistent pro-thrombotic signals in convalescence

  Patients with end-stage kidney disease (ESKD) are at high risk of severe COVID-19. Here, we perform longitudinal blood sampling of ESKD haemodialysis patients with COVID-19, collecting samples pre-infection, serially during infection, and after clinical recovery. Using plasma proteomics, and RNA-sequencing and flow cytometry of immune cells, we identify transcriptomic and proteomic signatures of COVID-19 severity, and find distinct temporal molecular profiles in patients with severe disease. Supervised learning reveals that the plasma proteome is a superior indicator of clinical severity than the PBMC transcriptome. We show that a decreasing trajectory of plasma LRRC15, a proposed co-receptor for SARS-CoV-2, is associated with a more severe clinical course. Strikingly, we observe that two months after the acute infection, patients still display dysregulated gene expression related to vascular, platelet and coagulation pathways, including PF4 (platelet factor 4), which may explain the prolonged thrombotic risk following COVID-19.

  Study EGAS00001006778

• Single-cell atlas of multiple myeloma and precursor diseases

  Single cell profiling of small cohorts of myeloma precursor disease patient samples including monoclonal gammopathy of unknown significance and smoldering multiple myeloma have shown diverse evolutionary patterns and immune changes that occur early in the disease process. Using scRNAseq and scBCRseq in a large cohort of fifty-three patients with myeloma precursor disease in comparison with myeloma and normal donors we describe the early genomic drivers of malignant transformation and describe their divergent clonal expansion in hyperdiploid compared to non-hyperdiploid samples. We describe intra-patient heterogeneity with potential therapeutic implications as well as the distinct evolution patterns (linear/branching) from myeloma precursor disease to myeloma. Finally, we describe the unique adaptations of the microenvironment as a response to distinct genomic changes in myeloma cells. These results further our knowledge to characterize myeloma precursor disease evolution, inform individual patient progression risk stratification and identify potential biomarkers that could be clinically exploited.

  Study EGAS00001006694

• Malignant progression of an ancestral bone marrow clone harboring a CIC-NUTM2A fusion in isolated myeloid sarcoma

  Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia (AML) or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis.

  Study EGAS00001006833

• WGS of MAPKi acquired resistant samples from patients and PDX models

  Blocking cancer genomic instability, which generates diverse variants that enable rapid adaptations, may prevent tumor escape from therapies. We show that, after MAPK inhibitor (MAPKi) therapy in patients and patient-derived xenografts, acquired-resistant genomes of metastatic cutaneous melanoma accumulate amplicons of resistance-driver, non-homologous end-joining (NHEJ), and homologous recombination repair (HRR) genes within complex genomic rearrangements (CGRs) and extrachromosomal circular DNAs (ecDNAs). Also, chromothriptic genomic regions enrich for mutations and may engender ecDNAs and CGRs. Overall, somatic mutations within ecDNA- and CGR-amplicons enrich for HRR signatures, and, in acquired resistance, ecDNA (versus CGR) mutations enrich for APOBEC3 signatures. Breakpoint-sequence analysis suggests NHEJ as critical to double-stranded DNA break repair underlying CGR and ecDNA formation. In vivo, NHEJ-targeting by a DNA-PKCS inhibitor prevented acquired MAPKi-resistance by reducing the size of ecDNAs and CGRs early on combination treatment. Thus, targeting the cause(s), as opposed to the consequences, of genomic instability prevents acquired resistance.

  Study EGAS00001006874

• Blood transcriptome profiling links immunity to disease severity in myotonic dystrophy type 1 (DM1).

  The blood transcriptome was examined in relation to disease severity in type I myotonic dystrophy (DM1) patients who participated in the Observational Prolonged Trial In DM1 to Improve QoL- Standards (OPTIMISTIC) study. This sought to a) ascertain if transcriptome changes were associated with increasing disease severity, as measured by the muscle impairment rating scale (MIRS) and b) es-tablish if these changes in mRNA expression and associated biological pathways were also observed in the Dystrophia Myotonica Biomarker Discovery Initiative (DMBDI) microarray dataset in blood (with equivalent MIRS/DMPK repeat length). The changes in gene expression were compared using a number of complementary pathways, gene ontology and upstream regulator analyses, which suggested that symptom severity in DM1 was linked to transcriptomic alterations in innate and adaptive im-munity associated with muscle wasting. Future studies should explore the role of immunity in DM1 in more detail to assess its relevance to DM1.

  Study EGAS00001006926

• Multiome_HB

  Hepatoblastoma (HB) cells display strong phenotypic heterogeneity with a major impact on drug response, but the underlying mechanisms are poorly understood. Here, we use a single-cell multi-omic strategy to unravel the molecular determinants of this plasticity. HB display a continuum of single-cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activation underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation poles. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

  Study EGAS00001006932

• Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs drug target discovery in ankylosing spondylitis

  Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterised by enthesitis of the spine and sacroiliac joints. Genome-wide association studies have revealed >100 genetic associations but their functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We report differential disease signatures for specific genomic regions from individual -omic modalities and multi-omic integration, showing enrichment at genetically associated loci and evidence for altered monocyte function in AS. We link putative functional SNPs with cognate genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1. We also show how disease-specific functional genomic data can be integrated with GWAS evidence to enhance therapeutic target discovery. This study combines multiple epigenetic and transcriptional analysis with GWAS to identify cell types, candidate drug target genes and transcriptional regulation of likely pathogenic relevance.

  Study EGAS00001006945

• Single-cell RNA-seq and spatial transcriptomics data for the sarcoidosis baseline project

  Granulomas are lumps of immune cells that can form in various organs. Most granulomas appear unstructured, yet they have some resemblance to lymphoid organ formation. To better understand granuloma formation, we performed single-cell sequencing and spatial transcriptomics on granulomas from patients with sarcoidosis and bioinformatically reconstructed the underlying gene-regulatory networks. We discovered an immune stimulatory environment in granulomas that repurposed transcriptional programs associated with lymphoid organ development. Granuloma formation followed characteristic spatial patterns and involved genes linked to immunometabolism, cytokine and chemokine signaling, and extracellular matrix remodeling. Three cell types emerged as key players in granuloma formation: metabolically reprogrammed macrophages, cytokine-producing Th17.1 cells, and fibroblasts with inflammatory and tissue remodeling phenotypes. Pharmacological inhibition of one of the identified processes attenuated granuloma formation in a sarcoidosis mouse model. We show that human granulomas adopt characteristic aspects of normal lymphoid organ development in aberrant combinations, indicating that granulomas constitute aberrant lymphoid organs.

  Study EGAS00001006970

• Neuroblastoma heterogeneity

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

  Study EGAS00001007016

• Neuroblastoma heterogeneity

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

  Study EGAS00001007019

• Single-cell RNA-seq of cervix and placenta

  Placental trophoblasts have been detected in cervical smears early in gestation, creating opportunities for non-invasive prenatal diagnosis. However, trophoblast isolation is limited by the lack of a cell catalogue and molecular profile of cervical smears. To establish an atlas and to explore the potential of single-cell RNA-sequencing to detect cervical trophoblasts, 10,539 single-cell transcriptomes from 12 non-invasive exocervical smears from pregnant women were profiled. In addition, we characterized 34,565 cells from six placental biopsies. We uncovered a novel extravillous trophoblast cell subtype characterized by epithelial marker genes and reduced HLA-G expression. Integration of both cell atlases demonstrated surprisingly similar expression profiles between maternal epithelial cells and placental extravillous trophoblasts, indicating that these trophoblasts retained epithelial properties without an invasive mesenchymal phenotype. Additionally, differential expression analysis identified novel markers discriminating cervical and placental cell types. Those analyses may contribute to trophoblast discovery and isolation from cervical smears.

  Study EGAS00001007044

• Epigenetic landscape of mixed phenotype leukemias

  Leukemias with ambiguous lineage comprise a number of loosely defined entities, often without a clear mechanistic basis. Here, we investigated a group of such leukemias with a CpG Island Methylator Phenotype (CIMP), previously identified as CEBPA-silenced AML. Transcriptomics and epigenomics analyses revealed a hybrid myeloid/lymphoid epigenetic landscape, whereas genetic alterations were heterogenous. This suggests that CIMP leukemias are defined by their shared epigenetic profile rather than a common genetic lesion. Gene expression enrichment showed strong similarity with ETP-ALL and an early lymphoid progenitor cell of origin. Accordingly, integration of differential methylation and expression revealed widespread silencing of myeloid transcription factors (TFs), among which CEBPA was key for differentiation arrest. Hypermethylation also resulted in loss of CTCF binding, accompanied by a few changes in chromatin interactions involving critical TFs like KLF4. In conclusion, epigenetic dysregulation, and not genetic lesions, explain the mixed phenotype of a group of CIMP leukemias resembling ETP-ALL.

  Study EGAS00001007094

• Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes

  DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N=190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P<1x10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

  Study EGAS00001007118

• Spatial and temporal transcriptomics of medulloblastoma with chromothripsis identifies multiple genetic clones that resist to treatment and lead to relapse

  From all paediatric brain tumours, medulloblastomas with chromothripsis, a form of genome instability leading to massive genome rearrangements, are particularly aggressive. To dissect the complex biology underlying the aggressiveness of chromothriptic medulloblastomas, we applied spatial transcriptomics to i) chromothriptic and non-chromothriptic medulloblastomas from the same molecular subgroup (n=13) and ii) patient-derived xenografts from chromothriptic medulloblastoma, from the minimal residual disease stage to regrown tumours (n=11). Chromothriptic medulloblastomas showed higher spatial intra-tumour heterogeneity, proliferation and stemness, but lower immune infiltration and differentiation, as compared with non-chromothriptic medulloblastomas. Spatial mapping of genetic clones identified phenotypic features such as degree of differentiation, proliferation and immune infiltration enriched in specific clones. Cells from multiple genetic clones resisted treatment and gave rise to relapse in patient-derived xenografts of chromothriptic medulloblastoma. We identified a potential role of tumour microtubes in treatment resistance in chromothriptic medulloblastoma.

  Study EGAS00001007128

• Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia

  In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness, and poor patient outcomes. A later arrest was associated with lineage fidelity, deep and durable leukemia remissions, and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1+ pre-leukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI-efficacy is an unexpected outcome of the differentiation stage at which this leukemia transforms.

  Study EGAS00001007167

• The genetic history of the southern Andes from present-day Mapuche ancestry

  We generate genome-wide data from 64 participants from three Mapuche populations in Southern Chile: Pehuenche, Lafkenche, and Huilliche. Broadly, we describe three main ancestry blocks with a common origin, which characterize the Southern Cone, the Central Andes, and Amazonia. Within the Southern Cone, ancestors of the Mapuche lineages differentiated from those of the Far South during the Middle Holocene, and did not experience further migration waves from the north. We find that the deep genetic split between the Central and Southern Andes is followed by instances of gene flow, which may have accompanied the southward spread of cultural traits from the Central Andes, including crops and loanwords from Quechua into Mapudungun (the language of the Mapuche). Our findings add new perspectives on the genetic (pre)history of South America, from first settlement through to the present-day indigenous presence. Follow-up fieldwork took these results back to the indigenous communities to contextualize the genetic narrative alongside indigenous knowledge and perspectives.

  Study EGAS00001007200

• Single molecule genome-wide mutation profiles of cell-free DNA for non-invasive detection of cancer

  Somatic mutations are a hallmark of tumorigenesis and may be useful for non-invasive diagnosis of cancer. We analyzed whole-genome sequencing (WGS) data from 2,511 individuals in the Pan-Cancer Analysis of Whole Genomes (PCAWG) study as well as 489 individuals from four prospective cohorts and found distinct regional and mutation type specific frequencies in tissue and cell-free DNA (cfDNA) of cancer patients that were associated with replication timing and other chromatin features. A machine learning model using genome-wide mutational profiles combined with other features and followed by CT imaging detected >90% of lung cancer patients, including those with stage I and II disease. The fixed model was validated in an independent cohort, detected patients with cancer earlier than standard approaches, and could be used to monitor response to therapy. This approach lays the groundwork for non-invasive cancer detection using genome-wide mutation features that may facilitate cancer screening and monitoring.

  Study EGAS00001007248

• An isogenic human iPSC model unravels neurodevelopmental abnormalities in SMA

  Whether neurodevelopmental defects underlie the selective neuronal death that characterizes neurodegenerative diseases is becoming an intriguing question. To address it, we focused on the motor neuron (MN) disease Spinal Muscular Atrophy (SMA), caused by reduced levels of the ubiquitous protein SMN. Taking advantage of the first isogenic human induced pluripotent stem cell-derived SMA model that we have generated and a spinal cord organoid system, here we report that the relative and temporal expression of early neural progenitor and MN markers is altered in SMA. Furthermore, the corrected isogenic controls only partially reverse these abnormalities. These findings raise the relevant clinical implication that SMN-increasing treatments might not fully amend SMA pathological phenotypes. The approach we have taken demonstrates that the discovery of new disease mechanisms is greatly improved by using human isogenic models. Moreover, our study implies that SMA has a developmental component that might trigger the MN degeneration.

  Study EGAS00001007259

• RNAseq of Colorectal cancer organoid-stroma biobank cohort

  In colorectal cancers (CRC) the tumor microenvironment plays a key role for prognosis and therapy efficacy. Patient-derived tumor organoids (PDTOs) show enormous potential for preclinical testing, however, in purely epithelial cultures features including the ‘consensus molecular subtypes’ (CMS) are largely eradicated. To better reflect cell heterogeneity, we established the CRC organoid-stroma biobank of matched PDTOs and cancer-associated fibroblasts (CAFs) from 30 patients. Context-specific phenotyping showed that xenotransplantation or co-culture with CAFs restores the individual transcriptomic status and instructs subtype-specific stromal gene expression. Furthermore, co-culture exposed CMS4-specific resistance to Gefitinib and SN-38 and revealed prognostic gene expression signatures. Chemogenomic library screening identified patient- and therapy-dependent mechanisms of stromal resistance including MET as common target. Our results demonstrate that CRC phenotypes are encrypted in the cancer epithelium in a plastic fashion that strongly depends on the context. Consequently, CAFs are essential for faithful representation of molecular subtypes and therapy responses ex vivo.

  Study EGAS00001007300

• WES of Colorectal cancer organoid-stroma biobank cohort

  In colorectal cancers (CRC) the tumor microenvironment plays a key role for prognosis and therapy efficacy. Patient-derived tumor organoids (PDTOs) show enormous potential for preclinical testing, however, in purely epithelial cultures features including the ‘consensus molecular subtypes’ (CMS) are largely eradicated. To better reflect cell heterogeneity, we established the CRC organoid-stroma biobank of matched PDTOs and cancer-associated fibroblasts (CAFs) from 30 patients. Context-specific phenotyping showed that xenotransplantation or co-culture with CAFs restores the individual transcriptomic status and instructs subtype-specific stromal gene expression. Furthermore, co-culture exposed CMS4-specific resistance to Gefitinib and SN-38 and revealed prognostic gene expression signatures. Chemogenomic library screening identified patient- and therapy-dependent mechanisms of stromal resistance including MET as common target. Our results demonstrate that CRC phenotypes are encrypted in the cancer epithelium in a plastic fashion that strongly depends on the context. Consequently, CAFs are essential for faithful representation of molecular subtypes and therapy responses ex vivo.

  Study EGAS00001007301

• A longitudinal single-cell atlas of treatment response in pediatric AML

  Pediatric acute myeloid leukemia (pAML) is a disease characterized by heterogeneous cellular composition, driver alterations and prognosis. Characterization of this heterogeneity and how it affects treatment response remains relatively understudied in pediatric patients. We therefore used single-cell RNA sequencing and single-cell ATAC sequencing to profile 684,031 diagnosis, remission and relapse cells from 28 patients representing different pAML subtypes. Analysis revealed that at diagnosis, cellular compositions differed between distinct subgroups. Upon relapse, cellular hierarchies transitioned towards a more primitive state regardless of subtype. These primitive cells were distinct compared to cells at diagnosis, having underrepresentation of transcriptional programs involved in myeloid differentiation while programs commonly found in other lineages were overrepresented. In a subset of patients, this appeared to accompany the appearance of a B-lymphoid-like hierarchy. Together, our data reveal the emergence of apparent subtype-specific plasticity upon treatment and inform on potential novel therapeutic angles to target these cell populations.

  Study EGAS00001007323

• High-resolution, patient-level dissection of IL-23 blockade in cutaneous psoriasis

  Psoriasis vulgaris and other chronic inflammatory diseases improve markedly with therapeutic blockade of IL-23 signaling in T cells, but the genetic mechanisms of response remain poorly understood. We single-cell transcriptomically profiled CD45+ immune cells isolated from lesional psoriatic skin before and after IL-23 blockade. In clinically responsive patients, a psoriatic transcriptional signature in skin-resident memory T cells was sharply attenuated. In contrast, poorly responsive patients were distinguished by enduring T17 cell activation, a mechanism distinct from alternative cytokine signaling or downstream resistance. Spatial transcriptomic analysis suggests that successful IL-23 blockade requires dampening of > 90% of IL-17-induced signaling in lymphocyte-adjacent keratinocytes, an unexpectedly high threshold. We also detected a subset of persistent, disease-specific T17 abnormalities in responsive patients, revealing an irreversible cell identity that may necessitate ongoing IL-23 inhibition. Collectively, our data establishes a patient-level paradigm for dissecting response to immunomodulatory treatments.

  Study EGAS00001007373

• Immune landscape of oncohistone-mutant gliomas reveals diverse myeloid populations and tumor-promoting behavior

  Histone H3-mutant gliomas are deadly brain tumors characterized by a dysregulated epigenome and stalled differentiation. In contrast to the extensive datasets available on tumor cells, limited information exists on their tumor microenvironment (TME), including the immune infiltrate. Here, we characterize the TME of H3.3K27M and G34R/V-mutant gliomas, and multiple H3.3K27M mouse models, using transcriptomic and proteomic, including spatial single-cell approaches. Resolution of immune lineages indicates uniform high-infiltration of H3-mutant gliomas with diverse myeloid populations, high-level expression of immune checkpoints, and scarce lymphoid cells, findings uniformly reproduced in all H3.3K27M models tested. We show these myeloid populations communicate with H3-mutant cells mediating immunosuppression and sustaining tumor formation and maintenance. Dual inhibition of myeloid cells and immune checkpoints showed significant therapeutic benefit in pre-clinical syngeneic models. Our findings provide valuable characterization of the TME of oncohistone-mutant gliomas, and insights into means to modulate the myeloid infiltrate for the benefit of patients.

  Study EGAS00001007510

• A clinically applicable connectivity signature for glioblastoma includes the tumor network driver CHI3L1

  Tumor microtubes (TMs) connect glioma cells to a network with considerable relevance for tumor progression and therapy resistance. The determination of TM-interconnectivity in individual tumors has been challenging and the impact on patient survival unresolved. Here, a connectivity signature from single-cell RNA-sequenced (scRNA-Seq) xenografted primary glioblastoma (GB) cells has been established using a dye uptake methodology, confirmed with recording of cellular calcium epochs and validated with clinical correlations. Astrocyte-like and mesenchymal-like GB cells have the highest connectivity signature scores in scRNA-sequenced patient-derived xenografts and patient samples. In large GB cohorts, network connectivity correlated with the mesenchymal subtype and dismal patient survival. CHI3L1 has been identified and validated as a robust molecular marker of connectivity with functional relevance. The connectivity signature allows novel insights into brain tumor biology, provides a proof-of-principle that tumor cell TM-connectivity is relevant for patients’ prognosis, and serves as a robust prognostic biomarker.

  Study EGAS00001007611

• Early evolutionary branching across spatial domains predisposes to clonal replacement under chemotherapy in neuroblastoma

  Neuroblastoma (NB) is one of the most lethal childhood cancers due its propensity to treatment resistance. By spatial mapping of subclone geographies before and after chemotherapy across 89 tumor regions from 12 NBs, we find that densely packed territories of closely related subclones present at diagnosis are replaced under effective treatment by islands of distantly related survivor subclones, originating from a different most recent ancestor compared to lineages dominating before treatment. Conversely, in tumors that progressed under treatment, ancestors of subclones dominating later in disease are present already at diagnosis. Chemotherapy treated xenografts and cell culture models replicates these two contrasting scenarios and shows branching evolution to be a constant feature of proliferating NB cells. Phylogenies based on whole genome sequencing of 505 individual NB cells indicate that a rich repertoire of parallel subclones, emerges already with the first oncogenic mutations and lays the foundation for clonal replacement under treatment.

  Study EGAS00001007650

• Repli-seq data for 'Replication timing alterations are associated with mutation acquisition during tumour evolution in breast and lung cancer'

  The essential process of DNA replication timing during each cell cycle is highly regulated ensuring the correct duplication of the genome. The extent and importance of alterations in this process during malignant transformation has not been extensively explored. Here, we evaluate the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find 6%-18% of the cancer genome displays ART. Genomic regions with a change from early to late replication exhibit an increased mutation rate and distinct mutational signatures, while regions changing from late to early contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART is a relatively early event during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

  Study EGAS00001007773

• Myeloid-specific KDM6B inhibition sensitizes Glioblastoma to PD1 blockade

  Glioblastoma (GBM) tumors are enriched in immune-suppressive myeloid cells and are refractory to immune checkpoint therapy (ICT). Targeting epigenetic pathways to reprogram the functional phenotype of immune-suppressive myeloid cells to overcome resistance to ICT remains unexplored. Single-cell and spatial transcriptomic analyses of human GBM tumors demonstrated high expression of an epigenetic enzyme - histone 3 lysine 27 demethylase (KDM6B) in intra-tumoral immune-suppressive myeloid cell subsets. Importantly, myeloid-cell specific Kdm6b deletion enhanced pro-inflammatory pathways and improved survival in GBM tumor-bearing mice. Mechanistic studies elucidated that the absence of Kdm6b enhances antigen-presentation, interferon response and phagocytosis in myeloid cells by inhibiting mediators of immune suppression including Mafb, Socs3 and Sirpa. Further, pharmacological inhibition of KDM6B mirrored the functional phenotype of Kdm6b deleted myeloid cells and enhanced anti-PD1 efficacy. Thus, this study identified KDM6B as an epigenetic regulator of the functional phenotype of myeloid cell subsets and a potential therapeutic target to improve response to ICT.

  Study EGAS00001007002

• Mutational Mechanisms in Multiply Relapsed Pediatric Acute Lymphoblastic Leukemia

  Pediatric acute lymphoblastic leukemia (ALL) is marked by low mutational load at initial diagnosis, which increases at relapse. To determine which processes are active in (relapsed) ALL and how they behave during disease progression before and after therapy, we performed whole genome sequencing on tumor samples of 29 multiply relapsed ALL patients. Mutational load increased upon relapse in 28 patients and every subsequent relapse in 22 patients, with UV-like damage, APOBEC activity, reactive oxygen species, thiopurine-associated damage and an unknown therapy component driving mutagenesis. Thiopurines were the most prominent source of new mutations in relapse, affecting over half of the studied patients and causing potential relapse-driving mutations in multiple patients. Mutational processes often affected patients over longer time periods, but could also occur in isolated events, suggesting the requirement of additional triggers. Our data pose mutational processes as prominent secondary drivers with an important role in ALL development and progression.

  Study EGAS00001007900

• Mapping the human hematopoietic stem and progenitor cell hierarchy through integrated single-cell proteomics and transcriptomics

  Single-cell transcriptomics (scRNA-seq) has enabled the characterization of cell state heterogeneity and recapitulation of differentiation trajectories. However, since proteins are the main functional entities in cells, the exclusive use of mRNA measurements comes at the risk of missing important biological information. Here we leverage recent technological advances in single-cell proteomics by Mass Spectrometry (scp-MS) to generate the first scp-MS dataset of an in vivo differentiation hierarchy encompassing over 2,500 human CD34+ hematopoietic stem and progenitor cells. Through integration with scRNA-seq, we identify proteins that are important for stem cell quiescence, which were not indicated by their mRNA transcripts, and demonstrate functional expression covariance during differentiation that is only detectable on protein level. Finally, we show that modeling translation dynamics can infer cell progression during differentiation and explain 45% more protein variation from mRNA than linear correlation. Our work serves as a framework for future single-cell multi-omics studies across biological systems.

  Study EGAS00001007930

• Novel paediatric case of a spinal high-grade astrocytoma with piloid features in a patient with Noonan Syndrome

  Noonan Syndrome (NS) is associated with an increased risk of low-grade central nervous system tumours in children but only very rarely associated with high-grade gliomas. Here we describe the first reported case of a spinal high-grade astrocytoma with piloid features (HGAP) in a child with NS. This case was a diagnostic and treatment dilemma, prior to whole-genome germline and tumour sequencing, tumour transcriptome sequencing and DNA methylation analysis. The methylation profile matched strongly with HGAP and sequencing identified somatic FGFR1 and NF1 variants and a PTPN11 germline pathogenic variant. Therapeutic targets were identified but also alterations novel to HGAP such as differential expression of VEGFA and PD-L1. The germline PTPN11 finding has not been previously described in individuals with HGAP.  This case underscores the power of precision medicine from a diagnostic, therapeutic and clinical management perspective, and describes an association between HGAP and NS which has not previously been reported.

  Study EGAS00001007937

• Cross-species Comparison Reveals Therapeutic Vulnerabilities Halting Glioblastoma Progression

  The growth of a tumor is tightly linked to the distribution of its cells along a continuum of activation states. Here, we systematically decode the activation state architecture (ASA) in a glioblastoma (GBM) patient cohort through comparison to adult murine neural stem cells. Modelling of these data forecast how tumor cells organize to sustain growth and identifies rate of activation as the main predictor of growth. Accordingly, patients with a higher quiescence fraction exhibit improved outcomes. Further, DNA methylation arrays enable ASA- related patient stratification. Comparison of healthy and malignant gene expression dynamics reveals dysregulation of the Wnt- antagonist SFRP1 at the quiescence to activation transition. SFRP1 overexpression renders GBM quiescent and increases overall survival of tumor-bearing mice. Surprisingly, it does so through reprogramming the tumor’s stem-like methylome into an astrocyte- like one. Our findings offer a framework for patient stratification with prognostic value, biomarker identification, and therapeutic avenues to halt GBM progression.

  Study EGAS00001008155

• Molecular and clinical diversity in primary central nervous system lymphoma: a LOC network multi-omic PCNSL study

  Primary central nervous system lymphoma (PCNSL) is a distinct extranodal lymphoma presenting with limited stage disease but variable response rates to treatment despite homogenous pathological presentation. The likely underlying molecular heterogeneity and its clinical impact is poorly understood. The present dataset contents paired-ended whole-exome sequencing information (n=115; HyperExome Kapa hyperprep), paired-ended RNA-seq information (n=123; KAPA mRNA HyperPrep Kits), and paired-ended bisulfite sequencing (n=64; TruSeq Methyl Capture EPIC) from fresh-frozen tumor tissue immunocompetent, treatment naïve PCNSL patients. Additionally, the dataset includes single-ended RNA-seq (n=93; QuantSeq 3’ mRNA-Seq Library Prep Kit) from formalin-fixed, paraffin-embedded tissue of immunocompetent, treatment naïve PCNSL patients. All samples were sequenced in an Illumina NovaSeq 6000 instrument.

  Dataset EGAD00001008706

• RaScALL: Rapid (Ra) screening (Sc) of RNA-seq data for prognostically significant genomic alterations in acute lymphoblastic leukaemia (ALL)

  RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, a full RNA-seq Bioinformatics workflow is time-consuming and costly in a clinical setting where rapid detection and accurate reporting of clinically relevant alterations are essential. To accelerate the identification of ALL-specific alterations (including gene fusions, single nucleotide variants and focal gene deletions), we developed the rapid screening tool RaScALL, capable of identifying more than 100 prognostically significant lesions directly from raw sequencing reads. RaScALL uses the k-mer based targeted detection tool km and known ALL variant information to achieve a high degree of accuracy for reporting subtype defining genomic alterations compared to standard alignment-based pipelines. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% (164 samples) of reported cases in a 180-patient Australian study cohort and 95% (n=63) of samples in a North American validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n=12) and IKZF1 p.N159Y (n=4). Accurate detection of intragenic IKZF1 deletions resulting in aberrant transcript isoforms was also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches, ensuring accelerated risk-stratification and therapeutic triage.

  Dataset EGAD00001009087

• Diagnostic utility of whole genome sequencing in adults with B-other acute lymphoblastic leukemia

  Genomic profiling at diagnosis of B-cell precursor Acute Lymphoblastic Leukemia (BCP-ALL) in adults is used to guide disease classification, risk stratification and treatment decisions. Patients for which diagnostic screening fails to identify disease defining or risk stratifying lesions are classified as B-other ALL. We screened a cohort of 652 BCP-ALL cases enrolled in UKALL14 to identify and perform whole genome sequencing (WGS) on paired tumor-normal samples. For 52 B-other patients we compared WGS findings to data from clinical and research cytogenetics. WGS identifies a cancer associated event in 51/52 cases, this includes an established subtype defining genetic alteration in 5/52 that were previously missed by standard-of-care genetics. Of the 47 true B-other ALL we identified a recurrent driver in 87% (41).  Complex karyotype by cytogenetics emerges as a heterogeneous group, underlied by distinct genetic alterations associated with either favorable (DUX4-r) or poor outcomes (MEF2D-r, IGK::BCL2). For a subset of 31 cases, we integrate findings from RNA-sequencing (RNA-seq) analysis to include fusion gene detection, and classification by gene expression. Compared to RNA-seq, WGS was sufficient to detect and resolve recurrent genetic subtypes, however RNA-seq can provide orthogonal validation of findings. In conclusion, we demonstrate that WGS can identify clinically relevant genetic abnormalities missed by standard-of-care testing and identify leukemia driver events in virtually all cases of B-other ALL.

  Dataset EGAD00001009304

• Recurrent DNMT3B gene rearrangements are associated with unfavorable outcome in dicentric (9;20)-positive pediatric BCP-ALL

  We provide clinical data sets of Array CGH, targeted RNA-seq, total RNA-seq, whole genome bisulfite (WGBS) and whole genome DNA sequencing (WGS) obtained from bone marrow (BM) or peripheral blood (PB) mononuclear cells of 57 pediatric patients with dicentric chromosome dic(9;20) positive Acute lymphocytic leukemia (ALL), from which in 6 cases DNMT3B gene rearrangement was identified. This data is complemented by total RNA-seq and WGBS of samples from 4 additional ALL patients with a t(12;21) translocation and ETV6-RUNX1 gene fusion. DNA was isolated from BM or PB B-lymphocytes using Qiagen QIAamp DNA Blood Midi Kit to perform i) Array CGH of 58 dic(9;20) positive samples by hybridizing 500ng DNA using a Agilent 400K SurePrint G3 Custom CGH Human Genome Microarray (e-Array design 84704) ii) WGBS of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples using Tecan TrueMethyl oxBS-Seq module for library preparation and Illumina NovaSeq 6000 platform to run 2x151 cycles iii) WGS of DNMT3B rearrangement positive samples using Illumina Lotus DNA Library Prep Kit followed by sequencing running 2x160 cycles on an Illumina NovaSeq 6000 platform. RNA was isolated from PB B-lymphocytes using the PerkinElmer Chemagic 360 instrument, followed by i) targeted RNA-seq of 56 dic(9;20) positive samples prepared using Illumina TruSight RNA Pan-Cancer Panel and sequenced on an Illumina MiSeq platform running 2x75 cycles ii) total RNA-seq of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples utilizing TruSeq Stranded Total RNA Library Prep Gold kit and running 2x100 cycles on an Illumina NovaSeq 6000 platform.

  Dataset EGAD00001011122

• Comparison of sequencing assays for sensitive detection of circulating tumour DNA in stage IA-IV breast cancer

  This study compared different assays for the detection of circulating tumour DNA (ctDNA) in serial plasma from stage IA-IV breast cancer patients, targeting structural variants (SVs), single nucleotide variants (SNVs) and/or somatic copy-number aberrations (SCNAs). SV-multiplex PCR, SNV-/SV-hybrid capture, and different depths of whole-genome sequencing (WGS) were used to evaluate ctDNA levels, demonstrating concordant results. SNV-hybrid capture targeting 1,347-7,491 mutations was the most sensitive assay, detecting 67% (36/54) of samples down to an allele fraction (AF) of 0.00024%. SV-multiplex PCR, targeting 21-47 mutations, detected 63% (34/54) of samples down to 0.00047% AF and has potential as a clinical assay.

  Dataset EGAD00001008589

• Mutation analysis of core biopsies from localized breast cancer patients

  Targeted next-generation sequencing (NGS) of 93 frequently mutated genes in breast cancer using the QIAseq Human Breast Cancer Targeted Panel (QIAGEN), which uses digital sequencing by incorporating unique molecular barcodes (UMI).

  Dataset EGAD00001008396

• ANGIOPREDICT - an FP7-funded project enabling personalised medicine for patients with metastatic colorectal cancer.

  Exome (*_{N,T}{1,2}) RNAseq (polyA - *_PolyA, and RiboZero - *_RibZ) Methylation (SeqCapEpi - MAPD*).

  Dataset EGAD00001008463

• WGS of high grade serous ovarian tumours and matched blood normals

  FASTQ reads for 34 matched tumour-normal WGS pairs for high grade serous ovarian cancer patients. Scottish HGSOC samples were collected via local Bioresource facilities at Edinburgh, Glasgow, Dundee and Aberdeen and stored in liquid Nitrogen until required. HGSOC patients were determined from pathology records and were included in the study where there was matched tumour and whole blood samples. Tumour samples were divided into two for DNA and RNA extraction and slivers of tissue were taken, fixed in formalin and embedded in paraffin wax (FFPE). Samples were only included if they were confirmed as HGSOC and there was greater than 40% tumour cellularity throughout the tumour, determined using H&sE staining of the FFPE sections and pathology review. Somatic DNA was extracted from the tumour and germline DNA was extracted from whole blood. Somatic DNA was extracted using the Qiagen DNeasy Blood and tissue kit (cat no 69504). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including RNase digestion step). Germline DNA was extracted from 1-3ml whole blood using the Qiagen FlexiGene kit (cat no 51206) following the manufacturers recommended protocol. The resulting DNA underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as surrogate measures of DNA purity. A260/280 had to be 1.8 or greater and A260/230 had to be 2.0 or greater. Then, DNA was quantified using LifeTechnologies Qubit dsDNA BR kit (cat no Q32850) and we required a minimum of 50ul at 25ng/ul for WGS. Thirdly, DNA was diluted to 25ng/ul and a representative sample was loaded onto a 0.8% TAE gel, ran at 100v for 60mins and then imaged using a BioRad ChemiDoc imaging system to visualise the DNA quality.

  Dataset EGAD00001010135

• Whole exome and Transcriptome sequencing of treatment-naïve esophageal adenocarcinoma biopsies and matched peripheral blood mononuclear cells

  Exome sequencing was performed on n=28 treatment-naïve esophageal adenocarcinoma (EACs). Three to four biopsies sampling different areas of each tumor were pooled before nucleic acid extractions to mitigate the elevated heterogeneity described for EAC. WES was performed on EAC biopsies at 120X average coverage, with autologous PBMCs used as germline controls at 80X average coverage. Libraries were prepared from 30 ng of input DNA using the SureSelect QXT Human All Exon V7 kit (Agilent Technologies) and sequenced on the NextSeq 550 (Illumina), 2x150 bp. BCL files were demultiplexed to FastQ files using bcl2fastq2 software (Illumina). Three paired end sequencing batches were analyzed independently (Batch1: samples 8, 10, 11, 12, 15, 17, 18; Batch2: samples 20, 24, 25, 26, 27, 29, 30, 31, 33, 34 ; Batch3: samples 35, 37, 39, 40, 41, 43, 45, 48, 54, 55, 57). RNA sequencing was performed on n=26 treatment-naïve esophageal adenocarcinoma (EACs). Three to four biopsies sampling different areas of each tumor were pooled before nucleic acid extractions to mitigate the elevated heterogeneity described for EAC. RNAseq libraries were prepared on 50 ng of total RNA (with RNA integrity index RIN >=7) with the TruSeq Stranded mRNA library preparation kit (Illumina) in accordance with low-throughput protocol. After PCR enrichment (15 cycles) and purification of adapter-ligated fragments, the concentration and length of DNA fragments were measured using D1000 Screen Tape System (Agilent), obtaining a median insert size of 311 nucleotides. Then, RNAseq libraries were sequenced using the Illumina NovaSeq platform, 1x100 bp, obtaining on average 100 million single reads per sample.

  Dataset EGAD00001010876

• 1M-scBloodNL

  Here we present the 1M-scBloodNL study, in which we performed single-cell RNA-seq on 120 individuals of the Northern Netherlands population cohort Lifelines. For each individual peripheral blood mononuclear cells (PBMCs) were sequenced in an unstimulated condition, and after 3 and 24 hour in vitro stimulation with C. albicans (CA), M. tuberculosis (MTB) and P. aeruginosa (PA), totalling approximately 1.3 million cells. scRNA-seq was conducted with the 10X Genomics 3'-end v2 (72 libraries) and v3 (33 libraries) technology. In general, each library contains PBMCs from 8 donors and 2 different stimulation-timepoint combinations. Donors were demultiplexed using a combination of SoupOrCell (https://www.nature.com/articles/s41592-020-0820-1) and genotype information to assign the correct donor to a donor-specific cell cluster.

  Dataset EGAD00001007764

• Identification of the Global miR-130a Targetome Reveals a Novel Role for TBL1XR1 in Hematopoietic Stem Cell Self-Renewal and t(8;21) AML

  In this study, we identified miR-130a as a regulator of HSC self-renewal and differentiation. To characterize gene expression changes following enforced expression of miR-130a OE, we performed RNA-seq in CD34+ cord blood (CB) cells transduced with control and miR-130a OE lentiviruses. To capture miRNA targets in an unbiased, transcriptome-wide manner, we perfomed enhanced CLIPseq procol in 2 replicates of CD34+ CB cells and Kasumi-1 cell line, which represent a model system for t(8;21) AML. We chose this cell line, as we found miR-130a to be highly expressed in this AML subtype where it is critical for maintaining the oncogenic molecular program mediated by AML1-ETO. Chimeric Ago2 eCLIPseq in CD34+ CB cells combined with Mass Spectrometry data analysis identified TBL1XR1 as a principal target of miR-130a. To elucidate gene expression changes associated with TBL1XR1 loss of function, we performed RNA-seq in CD34+CD38- CB cells transduced with control and shRNA targeted against TBL1XR1. To determine the functional significance of high miR-130a expression levels in Kasumi-1 cells on the molecular network controlled by AML1-ETO, we performed CUT&RUN assay and RNA-seq in Kasumi-1 cells following miR-130a knock-down (KD). Collectively, our findings reveal a unique role of miR-130a in regulating normal hematopoietic stem cell self-renewal and how elevated levels of miR-130a in t(8;21) AML contribute to the leukemogenesis of this AML subtype.

  Dataset EGAD00001008412

• Genomics of Malignant Peripheral Nerve Sheath Tumor (MPNST)

  The Genomics of MPNST (GeM) Consortium dataset includes de-identified whole genome sequencing data (.bam) for germline samples (DNA primarily derived from blood) sequenced at standard (30x) coverage (n=88) and for tumor samples (DNA derived from fresh frozen tissue) sequenced at 90x coverage (n=105). This dataset also includes transcriptome profiling data (.fastq) for paired normal nerve samples (n=7) and for tumor samples (n=132).

  Dataset EGAD00001008608

• WGS of high grade serous ovarian tumours and matched blood normals

  FASTQ reads for 81 matched tumour-normal WGS pairs for high grade serous ovarian cancer patients. Scottish HGSOC samples were collected via local Bioresource facilities at Edinburgh, Glasgow, Dundee and Aberdeen and stored in liquid Nitrogen until required. HGSOC patients were determined from pathology records and were included in the study where there was matched tumour and whole blood samples. Tumour samples were divided into two for DNA and RNA extraction and slivers of tissue were taken, fixed in formalin and embedded in paraffin wax (FFPE). Samples were only included if they were confirmed as HGSOC and there was greater than 40% tumour cellularity throughout the tumour, determined using H&E staining of the FFPE sections and pathology review. Somatic DNA was extracted from the tumour and germline DNA was extracted from whole blood. Somatic DNA was extracted using the Qiagen DNeasy Blood and tissue kit (cat no 69504). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including RNase digestion step). Germline DNA was extracted from 1-3ml whole blood using the Qiagen FlexiGene kit (cat no 51206) following the manufacturers recommended protocol. The resulting DNA underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as surrogate measures of DNA purity. A260/280 had to be 1.8 or greater and A260/230 had to be 2.0 or greater. Then, DNA was quantified using LifeTechnologies Qubit dsDNA BR kit (cat no Q32850) and we required a minimum of 50ul at 25ng/ul for WGS. Thirdly, DNA was diluted to 25ng/ul and a representative sample was loaded onto a 0.8% TAE gel, ran at 100v for 60mins and then imaged using a BioRad ChemiDoc imaging system to visualise the DNA quality.

  Dataset EGAD00001009049

• RNA sequencing data for primary and recurrent ovarian granulosa cell tumors

  RNA sequencing dataset for primary and recurrent ovarian granulosa cell tumors consists of 24 .bam files with aligned reads, including 8 primary and 16 recurrent tumors. Total RNA was extracted from cryopreserved tissue of adult-type granulosa cell tumor samples. Libraries were prepared from cDNA using the NuGEN Ovation Ultralow Library System V2 (San Carlos, CA). Paired-end bulk RNA sequencing was performed on the Illumina HiSeq 2000 platform. RNA sequencing reads were aligned to the hg19/GRCh37 reference human genome using the STAR software (version 2.6.0b) with default parameters. Samples description: GCT001 Adult-type Granulosa Cell Tumor: recurrent tumor GCT002 Adult-type Granulosa Cell Tumor: primary tumor GCT003 Adult-type Granulosa Cell Tumor: recurrent tumor GCT004 Adult-type Granulosa Cell Tumor: primary tumor GCT005 Adult-type Granulosa Cell Tumor: recurrent tumor GCT006 Adult-type Granulosa Cell Tumor: recurrent tumor GCT007 Adult-type Granulosa Cell Tumor: primary tumor GCT008 Adult-type Granulosa Cell Tumor: recurrent tumor GCT009 Adult-type Granulosa Cell Tumor: recurrent tumor GCT010 Adult-type Granulosa Cell Tumor: recurrent tumor GCT011 Adult-type Granulosa Cell Tumor: primary tumor GCT012 Adult-type Granulosa Cell Tumor: recurrent tumor GCT013 Adult-type Granulosa Cell Tumor: recurrent tumor GCT014 Adult-type Granulosa Cell Tumor: recurrent tumor GCT015 Adult-type Granulosa Cell Tumor: recurrent tumor GCT016 Adult-type Granulosa Cell Tumor: recurrent tumor GCT017 Adult-type Granulosa Cell Tumor: recurrent tumor GCT018 Adult-type Granulosa Cell Tumor: primary tumor GCT019 Adult-type Granulosa Cell Tumor: recurrent tumor GCT020 Adult-type Granulosa Cell Tumor: primary tumor GCT021 Adult-type Granulosa Cell Tumor: primary tumor GCT022 Adult-type Granulosa Cell Tumor: recurrent tumor GCT023 Adult-type Granulosa Cell Tumor: recurrent tumor GCT024 Adult-type Granulosa Cell Tumor: primary tumor

  Dataset EGAD00001009108

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedNanoSeq_Buccal

  Targeted NanoSeq in buccal epithelium obtained longitudinally from a number of donors. A custom-designed panel of ~200 genes will be used to achieve high depth coverage (aiming for ~1000X at panel sites).As we age, many tissues become colonised by microscopic clones carrying somatic driver mutations. Some of these clones represent a first step towards cancer whereas others may contribute to ageing and other diseases. However, our understanding of the clonal landscapes of human tissues, and their impact on cancer risk, ageing and disease, remains limited due to the challenge of detecting somatic mutations present in small numbers of cells. Here, we introduce a new version of nanorate sequencing (NanoSeq), a duplex sequencing method with error rates of less than 55 per billion base pairs, which is compatible with whole-exome and targeted gene sequencing. Deep sequencing of polyclonal samples with single-molecule sensitivity enables the simultaneous detection of mutations in large numbers of clones, yielding accurate somatic mutation rates, mutational signatures and driver mutation frequencies in any tissue. Applying targeted NanoSeq to 1,042 non-invasive samples of oral epithelium and 371 samples of blood from a twin cohort, we found an unprecedentedly rich landscape of selection, with 46 genes under positive selection driving clonal expansions in the oral epithelium, over 62,000 driver mutations, and evidence of negative selection in some genes. The high number of positively selected mutations in multiple genes provides high-resolution maps of selection across coding and non-coding sites, a form of in vivo saturation mutagenesis. Multivariate regression models enable mutational epidemiology studies on how carcinogenic exposures and cancer risk factors, such as age, tobacco or alcohol, alter the acquisition and selection of somatic mutations. Accurate single-molecule sequencing has the potential to unveil the polyclonal landscape of any tissue, providing a powerful tool to study early carcinogenesis, cancer prevention and the role of somatic mutations in ageing and disease.

  Dataset EGAD00001015618