14459 results for ""Pennsylvania+German"+ACS+study+'Socius'+open+access"

in 147.11 milliseconds.

• Microbiomic and immunogenic biomarkers of adjuvant chemotherapy efficacy in stage III colorectal cancer

  The standard treatment for stage III colorectal cancer (CRC) is surgery, followed by adjuvant chemotherapy (AC). However, the efficacy of AC is limited, and approximately 30% of patients experience recurrence. Therefore, easily assessable molecular biomarkers that can predict response to AC are needed. Here, we aimed to identify predictive biomarkers of response to AC using tumor transcriptomic data. Tumor specimens from 253 patients who underwent surgery for stage III CRC between 1997 and 2019 were analyzed using RNA sequencing. Transcriptional subtyping and analyses of gut microbiome, immune cell fractions, and T-cell receptor (TCR) were performed to identify factors associated with recurrence-free survival (RFS) in patients treated with or without AC. Of 253 patients, 118 (46.6%) received AC. Twenty phyla and 799 species of the microbiomes were identified, among which Fusobacteria were identified in 247 patients (97.6%). Multivariable Cox proportional hazards regression analysis of patients who received AC showed that intratumoral Fusobacterium and M1 macrophage abundance and diversity of TCR alpha and delta (TRAD) were independently associated with RFS. Fusobacterium and M1 macrophage abundance, and TRAD diversity were associated with the response to AC, suggesting that these factors could serve as biomarkers for personalized treatment in stage III CRC.

  Study JGAS000875

• LySeqST: A targeted sequencing assay for robust genomic classification of diffuse large B-cell lymphoma

  Diffuse large B-cell lymphoma (DLBCL) is a biologically heterogeneous disease. Two genomic classification systems, LymphGen and DLBclass, are capable of reproducibly classifying single tumors into subtypes with prognostic relevance in the context of standard-of-care R-CHOP chemoimmunotherapy. Preliminary data from subgroup analyses of clinical trials suggest that the distinct biological features are targetable with drugs that exploit subtype-specific vulnerabilities, and large clinical trials are being designed to test this hypothesis. Although whole exome sequencing (WES) remains the gold standard for each of these classification systems, many laboratories may opt for smaller targeted sequencing panels that reduce costs associated with sequencing, data storage, and processing. Here, we describe the design and validation of a targeted sequencing panel (LySeqST) that captures the genomic features required specifically for LymphGen classification. We perform in silico and real-world validation vs. WES data and demonstrate that LySeqST can be used to accurately classify DLBCL tumors. In an unselected population-based cohort, we determine the real-world proportions and validate the prognostic relevance of LymphGen subtypes, confirming that only tumors expressing the dark zone gene expression signature are consistently associated with inferior outcomes. Our results support the use of LySeqST for accurate genomic classification of DLBCL.

  Study EGAS50000001601

• Exome sequencing identifies mutation of the ribosome in T-cell acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric malignancy and responds well to therapy in children, but not in adults. We used exome sequencing of 67 T-cell ALL (T- ALL) samples to gain insight in the mutational spectrum and age related differences of these leukemias. We detected protein-altering mutations in 508 genes, with an average of 8.2 mutations in pediatric and 21.0 in adult T-ALL. Based on stringent filtering, we predict 15 of these genes to be drivers, including 7 genes that were not previously implicated in T-ALL. We identify CNOT3 as a new tumor suppressor, as it is mutated in 12.2 % of adult T-ALLs and its knock-down causes de novo tumors in a sensitized drosophila eye tumor model. In addition, we identify mutations in RPL5 and RPL10, proteins of the 60S ribosomal complex2, in 18.2 % of pediatric T-ALLs. Mutation of arginine 98 in the X-linked RPL10 gene was most frequent. Yeast cells expressing the Rpl10 R98S mutant showed a clear ribosome biogenesis defect, excluding that this mutation is a passenger mutation. Our data provide insights in the mutational landscape of pediatric and adult T-ALL and identify the ribosome as a potential oncogenic factor in cancer.

  Study EGAS00001000296

• Single cell sequencing of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis

  The quantity of tumor infiltrating lymphocytes (TILs) in breast cancer (BC) is a robust prognostic factor for improved patient survival particularly in triple negative and HER2-overexpressing BC subtypes. Whilst T cells are the predominant TIL population, the relationship between quantitative and qualitative differences in T cell subpopulations and patient prognosis remains unknown. We performed single-cell RNA sequencing (scRNA-seq) of 6311 T cells isolated from human BCs, and show for the first time that significant heterogeneity exists in the infiltrating T cell population. We demonstrated that high TIL BCs contained significant numbers of CD8+ T cells with features of tissue-resident memory T (TRM) cell differentiation and that these CD8+ TRM cells expressed high levels of immune checkpoint molecules and effector proteins. A CD8+ TRM gene signature developed from the scRNA-seq data was significantly associated with improved patient survival in early-stage triple negative breast cancer (TNBC) and provided better prognostication compared with CD8 expression alone. Our data suggests that CD8+ TRM contribute to BC immunosurveillance and are the key targets of modulation through immune checkpoint inhibition. Further understanding of the development, maintenance and regulation of TRM cells will be crucial for successful immunotherapeutic development in BC.

  Study EGAS00001002845

• Shedding light over COVID-19 susceptibility and severity

  “One of the most remarkable features of SARS-CoV-2 infection is the variation in consequence ranging from asymptomatic to life-threatening” state the authors of the majestic Nature paper released in accelerated pre-view on July 2021. This variation in consequences represent not only an issue of burning medical relevance, but also a scientific and intellectual challenge for scientists around the globe. Researchers are known to be easily triggered by a challenge, and they recently demonstrated to be even better at double their efforts to focus on an urgent societal need. In this work led and coordinated by the COVID-19 Host Genetics Initiative (COVID-19 HGI), researchers gathered data from 49,562 cases and two million controls from 46 studies performed worldwide. The majority of the data has been deposited at the EGA, where 29 datasets have been collected under the study number EGAS00001005304. Data gathered from diverse sources conferred a huge robustness to the results. Genetic data from different ancestry origin were included in the analyses, even though data from patients with European origin was highly predominant (approx. 75%). Separated analysis of the different ancestries gave rise to relevant differences in the results, highlighting how important it is to take this aspect in consideration for future studies. Data from the individual studies was compiled into three categories according to COVID-19 severity for a meta-analysis: (a) critically ill cases (b) moderate or severe cases (c) all cases. They were paired with non-COVID-19-reported controls. All three analysis resulted in a total of 13 genetic loci associated with COVID-19 severity. Out of those, 9 loci associated specifically with critical illness, while 6 (a few overlapped) associated with moderated illness. Some genomic position matched previous findings. Interestingly, working with such an enormous amount of data allowed the depth to discriminate association with disease susceptibility vs severity, an issue of large medical relevance. Several loci mapped into biological relevant genes, like TYK2, PPP1R15A and FOXP4, previously implicated in immune response and/or lung functions. Researchers participating in the COVID-19 Host Genetics Initiative have worked in an unprecedented time-tight collaboration, prioritizing data sharing and collective results to tackle the pandemic affecting the world. The scientific community won’t miss any of the lessons to learn from this, and the EGA is ready to keep working to enable safe data sharing and fruitful collaborations.

  Blog covid-19-susceptibility-and-severity

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

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

  Study EGAS00001001862

• Comparison of EGF and PDGF driven glioblastomas.

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

  Study EGAS00001001900

• The evolutionary steps from primary to metastatic prostate cancer are largely uncharted, and the ability to use DNA present in body fluids as correlates of aggregate metastatic status is under-examined. We reconstructed phylogenies in ten prostate cancer patients with fatal disease using deep targeted sequencing of the prostate, adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. A total of 163 samples are studied.

  The evolutionary steps from primary tumor to metastasis in prostate cancer are largely uncharted, and the ability to use serum, plasma, and cerebrospinal fluid as a correlate of aggregate metastatic tumor genomic status has not been tested. We used deep targeted sequencing to reconstruct tumor evolution in ten prostate cancer patients with fatal disease encompassing examination of the prostate and adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. We show that there is substantial evolution from a common ancestor within the prostate that results in branching to multiple lineages which form an intermixed multi-clonal primary tumor mass. After the occurrence of key driver aberrations, one of these lineages will metastasize to multiple sites in a sequential fashion. These metastatic sites are then susceptible to being populated by cells from other intra-prostatic lineages or from other metastases. Genomic representation of metastases in body fluids is not uniform. Cerebrospinal fluid analysis can detect lineages not detected in circulating DNA, suggesting possible clinical utility.

  Study EGAS00001003848

• Pancreatic, Small-intestinal and Pulmonary Neuroendocrine Tumors

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

  Study EGAS00001004878

• Sclerosing epithelioid fibrosarcoma case report

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

  Study EGAS00001005214

• IDENTIFICATION AND TARGETED MANAGEMENT OF A PATIENT WITH A NEURODEGENERATIVE DISORDER CAUSED BY BIALLELIC MUTATIONS IN SLC5A6

  We describe a sibling pair displaying an early infantile-onset, progressive neurodegenerative phenotype, with symptoms of developmental delay and epileptic encephalopathy developing from 12 to 14 months of age. Using whole exome sequencing, compound heterozygous variants were identified in SLC5A6, which encodes the sodium-dependent multivitamin transporter (SMVT) protein. SMVT is an important transporter of the B-group vitamins biotin, pantothenate and lipoate. The protein is ubiquitously expressed and has major roles in vitamin uptake in the digestive system, as well as transport of these vitamins across the blood-brain barrier. Pathogenicity of the identified variants was demonstrated by impaired biotin uptake of mutant SMVT. Identification of this vitamin transporter as the genetic basis of this disorder guided targeted therapeutic intervention, resulting clinically in improvement of the patient’s neurocognitive and neuromotor function. This is the second report of biallelic mutations in SLC5A6 leading to a neurodegenerative disorder due to impaired biotin, pantothenate and lipoate uptake. The genetic and phenotypic overlap of these cases confirms mutations in SLC5A6 as the genetic cause of this disease phenotype. Recognition of the genetic disorder caused by SLC5A6 mutations is essential for early diagnosis and to facilitate timely intervention by triple vitamin (biotin, pantothenate and lipoate) replacement therapy.

  Study EGAS00001003861

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

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

  Study EGAS00001008305

• Characterizing microbiome-directed fibre snacks in gnotobiotic mice and humans

  Knowledge of the interrelationships between what we eat and the configurations of our gut microbial communities is providing important insights into how food components that are not directly metabolized by human enzymes are linked to our physiology and health status. Changing food preferences brought about by Westernization that have deleterious health effects1,2, plus rapid population expansion, ongoing challenges to sustainable agriculture, and other forces contributing to increased food insecurity, are catalyzing efforts to identify more nutritious and affordable foods3. The gut microbial community is complex, dynamic, and exhibits considerable intra- and interpersonal variation in its composition and functions. The massive number of potential interactions between its components makes it challenging to define the mechanisms by which food ingredients affect community properties. There is also a paucity of information about the ‘bioactive’ ingredients of foods that influence the fitness and expressed functions of community members. Here, plant fibres, from different sustainable sources and targeting distinct features of obese human gut microbiomes in gnotobiotic mice, were formulated into snack prototypes and used to supplement controlled diets consumed by overweight and obese adults; the results revealed fibre-specific changes in their microbiomes that were linked to changes in their plasma proteomes indicative of altered physiologic state.

  Study EGAS00001005268

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

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

  Study EGAS00001002377

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

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

  Study EGAS00001006090

• UK10K NEURO ABERDEEN

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

  Study EGAS00001000109

• Ethiopia_Genome_Project__high_coverage_

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

  Study EGAS00001000237

• Exome_trios_in_patients_with_gastroschisis

  Gastroschisis (MIM 230750) is a herniation of the intestines through a defect of the abdominal wall lateral to the umbilicus (usually on the right side), and it is not covered by a membrane [Ledbetter, 2012]. Gastroschisis is a congenital anomaly with increasing incidence, easy prenatal diagnosis and extremely variable postnatal outcomes. On the basis of clinical manifestations, epidemiologic charateristics, and the presence and type of additional malformations, gastroschisis could be considered a heterogeneous condition with no gene/s discovered yet. This congenital anomaly affects approximately 1-3 infancts per 10,000 live births [Calzolari et al.1995;Parker et al.,2010] Current knowledge about causative mutations/variants. To date, no single gene has been linked to gastroschisis. Some publications have tried to link this malformation to variants in genes (such as AEBP1 (adipocyte enhancer binding protein) gene [Feldkamp et al,. 2012] or the VEGF-NOS3 pathway [Lammer et al., 2008]. Previously, a Scribble mutant mouse model (circletail) was reported to exhibit gastroschisis, however recent studies demonstrated that the Scribble knockout fetus exhibits exomphalos phenotype of gastroschisis [Carnagham et al., 2013]. 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 EGAS00001002664

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

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

  Study EGAS00001003027

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

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

  Study EGAS00001000349

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

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

  Study EGAS00001007108

• ERG ALTERATIONS DEFINE A NOVEL SUBTYPE OF ACUTE LYMPHOBLASTIC LEUKEMIA

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

  Study EGAS00001000514

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

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

  Study EGAS00001001665

• Whole-Genome sequencing of hepatocellular carcinomas

  The French ICGC project on liver tumors is coordinated by Pr Jessica Zucman-Rossi and funded by Inca (French Institute for Cancer). The aim of the present project is to identify the catalog of somatic and germline mutations in liver tumors using whole genome and whole exome sequencing together with CGH-SNP, methylome and transcriptomic profiling. For this purpose, a series of 500 liver tumors are collected through the French National Liver Collection and these samples will be analyzed using the different omics technologies. The data will be deposited in the ICGC and EGA database to be publically available for the scientific community. Hepatocellular carcinoma (HCC) accounts for more than 90% of liver cancers, and is a major health problem. It is the 3rd cause of cancer-related mortality. Advances in genomic analyses have formed a comprehensive understanding of different underlying pathobiological layers resulting in hepatocarcinogenesis. Thus, the development of next-generation sequencing technologies has made it possible to generate more comprehensive catalogues of somatic alteration events (single nucleotide substitutions, structural variations, and epigenetic changes) in liver cancer genome than ever before. The dataset will include 50 whole genome sequencing tumor/germline pairs, of which 6 are deposited in February 2014.

  Study EGAS00001000706

• Single_Cell_Sequencing_of_Sperm__scSperm_

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

  Study EGAS00001000935

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

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

  Study EGAS00001001000

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal

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

  Study EGAS00001003774

• Rapid identification of somatic genome rearrangements as personalized biomarkers for blood-based cancer monitoring

  Increasing evidence shows the value of circulating tumour DNA (ctDNA) to detect cancer and monitor its progression. Somatic genomic structural variations (SVs) are promising personalized biomarkers for sensitive and specific detection of ctDNA in liquid biopsies. However, accurate, affordable, and fast identification of such SV biomarkers is challenging, which hinders routine use in the clinic. Here, we demonstrate a novel approach - termed SHARC - for rapid discovery of somatic SV breakpoints as personalized tumour biomarkers. SHARC combines low-coverage cancer genome sketching by using Oxford Nanopore portable sequencing with a random forest classification and a dedicated filtering pipeline to enrich for somatic SVs. Our method leverages the real-time and long-read capabilities of Nanopore sequencing to identify somatic SV breakpoints at nucleotide resolution from a tumour biopsy within two days. We applied SHARC to tumour samples of high-grade ovarian and prostate cancer and validated on average 10 somatic SVs per sample with the use of PCR mini-amplicons. Finally, we demonstrate that these somatic SV biomarkers can be used to detect tumour presence from liquid biopsies in a quantitative manner and we retrospectively monitored treatment response in patients with prostate cancer, demonstrating its potential benefit for clinical practice.

  Study EGAS00001003963

• Native American gene flow into Polynesia predating Easter Island settlement

  The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas while critics have argued that these botanical dispersals need not have been human mediated. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui). Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

  Study EGAS00001004209

• Pancreatic tropism of metastatic renal cell carcinoma

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

  Study EGAS00001004208

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

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

  Study EGAS00001002183

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

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

  Study EGAS00001002211

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

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

  Study EGAS00001000572

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

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

  Study EGAS00001001303

• Biallelic mutations in the ubiquitin ligase RFWD3 cause Fanconi anemia

  The WD40-containing E3 ubiquitin ligase RFWD3 has been recently linked to the repair of DNA damage by Homologous Recombination (HR). Here we show that an RFWD3 mutation within the WD40 domain is connected to the genetic disease Fanconi anemia (FA). An individual revealed congenital abnormalities characteristic for FA. Cells from the patient, carrying the compound heterozygous mutations c.205_206dupCC and c.1916T>A in RFWD3, show increased sensitivity to DNA interstrand cross-linking agents in terms of increased chromosomal breakage, reduced survival and cell cycle arrest in G2 phase. The cellular phenotype is mirrored in genetically engineered human and avian cells by inactivation of RFWD3 or by introduction of the patient-derived missense mutation, and is rescued by expression of wildtype RFWD3 protein. HR is disrupted in RFWD3 mutant cells, caused by impaired relocation of mutant RFWD3 to chromatin and defective physical interaction with RPA. Rfwd3 knockout mice exhibit increased embryonic lethality, are sub-fertile, show ovarian and testicular atrophy and have a reduced life span in which they resemble other FA mouse models. Although mutation was detected in a single child with FA yet, we propose RFWD3 as a novel FA gene, FANCW, supported by cellular paradigm systems and an animal model.

  Study EGAS00001002440

• High Grade Serous Ovarian Carcinomas Originate in the Fallopian Tube

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

  Study EGAS00001002589

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

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

  Study EGAS00001001113

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

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

  Study EGAS00001007512

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

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

  Study EGAS00001007941

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

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

  Study EGAS00001008169

• Genomic profiling of Acute Lymphoblastic Leukemia in Ataxia Telangiectasia patients reveals tight link between ATM mutations and chromothripsis

  Recent developments in sequencing technologies led to the discovery of a novel form of genomic instability, termed chromothripsis. This catastrophic genomic event, involved in tumorigenesis, is characterized by tens to hundreds of simultaneously acquired locally clustered rearrangements on one chromosome. We hypothesized that leukemias developing in individuals with Ataxia Telangiectasia, who are born with two mutated copies of the ATM gene, an essential guardian of genome stability, would show a higher prevalence of chromothripsis due to the associated defect in DNA double-strand break repair. Using whole-genome sequencing, fluorescence in situ hybridization and RNA sequencing, we characterized the genomic landscape of Acute Lymphoblastic Leukemia (ALL) arising in patients with Ataxia Telangiectasia. We detected a high frequency of chromothriptic events in these tumors, specifically on acrocentric chromosomes, as compared to tumors from individuals with other types of DNA repair syndromes (27 cases total, 10 with Ataxia Telangiectasia). Our data suggest that the genomic landscape of Ataxia Telangiectasia ALL is clearly distinct from that of sporadic ALL. Mechanistically, short telomeres and compromised DNA damage response in cells of Ataxia Telangiectasia patients may be linked with frequent chromothripsis. Furthermore, we show that ATM loss is associated with increased chromothripsis prevalence in additional tumor entities.

  Study EGAS00001002270

• Whole genome, whole exome, and targeted sequencing of high-grade meningioma tumor samples.

  High-grade meningiomas frequently recur and are associated with high rates of morbidity and mortality. To determine the factors that promote the development and evolution of these tumors, we analyzed the genomes of 274 high-grade meningiomas and compared this information with data from 456 low-grade meningiomas. High-grade meningiomas had a higher mutation burden than low-grade meningiomas but did not harbor any statistically significant mutated genes aside from NF2. High-grade meningiomas also possessed significantly elevated rates of chromosomal gains and losses, especially among tumors with monosomy 22. Meningiomas previously treated with adjuvant radiation had significantly more copy number alterations than radiation-induced or radiation-naïve meningiomas. Across serial recurrences, genomic disruption preceded the emergence of nearly all mutations, remained largely uniform across time, and when present in low-grade meningiomas, correlated with subsequent progression to a higher grade. In contrast to the largely stable copy number alterations, mutations were strikingly heterogeneous across tumor recurrences, likely due to extensive geographic heterogeneity in the primary tumor. While high-grade meningiomas harbored significantly fewer overtly targetable alterations than low-grade meningiomas, they contained numerous mutations that are predicted to be neoantigens, suggesting that immunologic targeting may be of therapeutic value.

  Study EGAS00001002294

• Botswana 15 autosomal unlinked microsatellites

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

  Study EGAS00001002380

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

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

  Study EGAS00001002432

• Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy

  Background: Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits including cardiovascular diseases.Results: Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 149 (97)* patients with dilated cardiomyopathy and 113 (108)* non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage and allele specific expression. Systematic annotation of genome wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome wide association signals in two independent cohorts.Conclusions: RNA transcription, splicing and allele specific expression are each important determinants of the DCM phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics. (*) numbers in parentheses represent samples that the publication is based on.

  Study EGAS00001002454

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

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

  Study EGAS00001002577

• UK10K_RARE_NEUROMUSCULAR

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

  Study EGAS00001000101

• CancerDetector: Ultrasensitive and Non-Invasive Cancer Detection at the Resolution of Individual Reads using Cell-free DNA Methylation Sequencing Data

  The detection of tumor-derived cell-free DNA in plasma is one of the most promising directions in cancer diagnosis. The major challenge in such approach is how to identify the tiny amount of tumor DNAs out of total cell-free DNAs in blood. Here we propose an ultrasensitive cancer detection method, termed “CancerDetector”, using the DNA methylation profiles of cell-free DNAs. The key of our method is to probabilistically model the joint methylation patterns of multiple adjacent CpG sites on an individual sequencing read, in order to exploit the pervasive nature of DNA methylation for signal amplification. Therefore, CancerDetector can sensitively identify a trace amount of tumor cfDNAs in plasma, at the level of individual reads. We evaluated CancerDetector on the simulated data, and showed a high concordance of the predicted and true tumor burden. Testing CancerDetector on real plasma data demonstrated its high sensitivity and specificity in detecting tumor DNAs. In addition, the predicted tumor burden showed great consistency with tumor size and survival outcome. Note that all of those testing were performed on sequencing data at low to medium coverage (1X to 10X). Therefore, CancerDetector holds the great potential to detect cancer early and cost-effectively.

  Study EGAS00001002728

• UK10K_RARE_FIND

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

  Study EGAS00001000128

• UK10K_RARE_CILIOPATHIES

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

  Study EGAS00001000126

• UK10K_RARE_SIR

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

  Study EGAS00001000130

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

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

  Study EGAS00001004290

• Stratifying and Targeting Pediatric Medulloblastoma through Genomics

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

  Study EGAS00001000273

• UK10K RARE CHD

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

  Study EGAS00001000125

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

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

  Study EGAS00001001028

• P647_Targeted_resequencing_project

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

  Study EGAS00001000305

• Exome sequencing of a novel cervical cancer cell line

  Human cancer cell lines are largely used in the searching for new antineoplastic agents. However, due to the artifacts of a long-term in culture, cell lines do not always represent the realistic tumor cell behavior. This has motivated the development of models that better mimetics the tumor tissue, among them, the establishment of primary cell cultures. In this work, we establish and characterized a low-passage cervix cancer cell line from a Brazilian patient with squamous cell carcinoma. The phenotype confirms the epithelial and tumor origin, through cytokeratins, EpCAM, and p16 staining. Whole exome sequencing showed relevant somatic mutations in several genes including BRCA2, TGFBR1, and IRX2 genes. CNV analysis by nanostring and WGS revealed amplification in genes related mainly with kinases proteins, involved in proliferation, migration and cell differentiation, such as EGFR, PIK3CA, and MAPK7. Overexpression of EGFR was confirmed by phospho RTK-array and western blot analysis. Furthermore, the cell was sensitivity to cisplatin, with IC50 13 times lower than SiHa cell line. In conclusion, this cervical cancer cell line presents molecular alterations that are an important tool for leading pre-clinical studies of new drugs that target one or more of the altered pathways.

  Study EGAS00001003343

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

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

  Study EGAS00001003438

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal_LCM

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

  Study EGAS00001003455

• Reconstituting the transcriptome and DNA methylome landscapes of human implantation at single-cell resolution

  Implantation is a milestone event during mammalian embryogenesis. Due to the extreme difficulty of obtaining in vivo human early post-implantation embryos, the gene regulatory network and epigenetics controlling human embryo implantation remains elusive. Here, combining an in vitro culture system for human post-implantation development and single-cell omics sequencing technologies, over 10,000 single cells at five representative stages of pre/post-implantation development were systematically analyzed. Unsupervised dimensionality reduction and clustering algorithm of the transcriptome data show stepwise implantation routes for the epiblast, primitive endoderm, and trophectoderm lineages, suggesting preparation for the establishment of a mother-to-offspring connection after implantation. Female embryos showed asynchronous progress of dosage loss of X chromosomes during implantation. Furthermore, using the single cell trio-Seq (scTrio-Seq) strategy, re-methylation of the genomes of all the three lineages was unambiguously revealed. Surprisingly, the genome re-methylation of PE lineage were much slower than both EPI and TE lineages during the implantation process, indicating distinct methylome features between EPI and PE although both of which were derived from ICM. Collectively, our work paves the way for understanding the complex molecular mechanisms that regulate human embryo implantation, informing new insights and future efforts in early embryonic development and reproductive medicine.

  Study EGAS00001003443

• Molecular determinants of response to PD-L1 blockade across tumor types

  Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer (NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity was observed between indications, in particular within the PD-L1+ tumors. mUC and NSCLC were molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor CDKN2A was identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis revealed molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.

  Study EGAS00001004343

• A somatic reference standard for cancer genome sequencing with COLO829

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

  Study EGAS00001001385

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

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

  Study EGAS00001001690

• Ethiopia_Genome_Project__low_coverage_

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

  Study EGAS00001000238

• UK10K OBESITY TWINSUK

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

  Study EGAS00001000306

• A system-wide approach to monitor responses to synergistic BRAF and EGFR inhibition in colorectal cancer cells

  Molecular heterogeneity of tumors, epigenetic changes and a diverse range of molecular mechanisms are main contributors to drug resistance, which represents one of the great challenges in cancer treatment. A deeper understanding of the molecular biology of cancer has resulted in better targeted therapies, where one or multiple drugs are adopted in novel therapies to tackle resistance. This beneficial effect of using combination treatments has to some extent also been observed in colorectal cancer (CRC) patients harboring the BRAF(V600E) mutation, whereby dual inhibition of BRAF(V600E) and EGFR increases antitumor activity. Notwithstanding this success, it is not clear whether this combination treatment is the only or most effective treatment to block resistance. Here, we investigate molecular responses upon single and multi-target treatments, over time, using BRAF(V600E) mutant CRC cells as a well-define model system. Through integration of transcriptomic, proteomic and phosphoproteomics data we obtain a comprehensive overview, revealing both well-known and novel responses. We primarily observe widespread upregulation of tyrosine kinases (RTKs) and and metabolic pathways. This points to by which the treated cells switch energy sources as a defensive response entering a quiescent like state, while activating signalling to re-activate the MAPK pathway.

  Study EGAS00001002654

• Bisulfite-converted duplexes for the strand-specific detection and quantification of rare mutations

  The identification of mutations that are present at low frequencies in clinical samples is an essential component of precision medicine. The development of molecular barcoding for next generation sequencing has greatly enhanced the sensitivity of detecting such mutations by massively parallel sequencing. However, further improvements in specificity would be useful for a variety of applications. We herein describe a technology (BiSeqS) that can increase the specificity of sequencing by at least two orders of magnitude over and above that achieved with molecular barcoding and can be applied to any massively parallel sequencing instrument. BiSeqS employs bisulfite treatment to distinguish the two strands of molecularly barcoded DNA; its specificity arises from the requirement for the same mutation to be identified in both strands. Because no library preparation is required, the technology permits very efficient use of the template DNA as well as sequence reads, which are nearly all confined to the amplicons of interest. Such efficiency is critical for clinical samples, such as plasma, in which only tiny amounts of DNA are often available. We show here that BiSeqS can be applied to evaluate transversions, as well as small insertions or deletions, and can reliably detect one mutation among >10,000 wild type molecules.

  Study EGAS00001002406

• Hermansky-Pudlak syndrome type 1 causes impaired anti-microbial immunity through a pathogenic lipid metabolism-mTOR circuit

  Mendelian diseases that present with immune-mediated disorders can provide insights into the molecular mechanisms that drive inflammation. Hermansky-Pudlak syndrome (HPS) types 1 and 4 are caused by defective vesicle trafficking involving the BLOC-3 complex. The presence of inflammatory complications such as Crohn’s disease-like inflammation and lung fibrosis in these patients remains enigmatic. Using mass cytometry we observe an augmented inflammatory monocyte compartment in HPS1 patient peripheral blood that may be associated with a TNF - and IL-1α-dominated cytokine dysregulation. HPS1 patient monocyte-derived macrophages express an inflammatory TNF-OSM mRNA gene signature and changes in lipid metabolism. Using stimulation experiments and lysosomal proteomics we show that defective lipid metabolism drives RAB32-dependent mTOR signaling, facilitated by the accumulation of mTOR on lysosomes. This pathogenic circuit translates into aberrant bacterial clearance, which can be rescued with mTORC1 inhibition. We reveal that a pathogenic lipid-mTOR signaling circuit acts as a metabolic checkpoint for defective anti-microbial activity. This mechanism may be relevant to the complex pathology of HPS1 patients featuring macrophage lipid accumulation, granuloma formation, defective anti-microbial activity and tissue inflammation. Lastly, this circuit may be present in a wider group of disorders with defective lipid metabolism and cholesterol accumulation.

  Study EGAS00001005053

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

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

  Study EGAS00001004291

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

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

  Study EGAS00001004018

• Hypoxia acts as an environmental cue for the human TRM differentiation program

  Tissue-resident memory T-cells (TRM) provide frontline defense against infectious diseases and contribute to anti-tumor immunity; however, aside from the necessity of TGF-β, knowledge regarding TRM-inductive cues remains incomplete, particularly for human cells. Oxygen tension is an environmental cue that distinguishes peripheral tissues from the circulation and here, we demonstrate that differentiation of human CD8+ T-cells in the presence of hypoxia and TGF-β1 led to the development of a TRM phenotype, characterized by a greater than five-fold increase in CD69+CD103+ cells expressing human TRM hallmarks and enrichment for endogenous human TRM gene signatures, including increased adhesion molecule expression and decreased expression of genes involved in recirculation. Hypoxia and TGF-β1 synergized to produce a significantly larger population of TRM phenotype cells than either condition alone, and comparison of these cells from the individual and combination conditions revealed distinct phenotypic and transcriptional profiles, indicating a programming response to milieu rather than a mere expansion. Our findings identify a previously unreported cue for the TRM differentiation program and can enable facile generation of human TRM phenotype cells in vitro for basic studies and translational applications such as adoptive cellular therapy.

  Study EGAS00001005286

• To determine the mutational impact of the in vitro culture, clonal human adult and pluripotent stem cell lines were subjected to a second clonal step after 3 months of culture. These subclones were whole genome sequenced to identify all the mutations that accumulated during the 3 month culture period.

  Genetic changes acquired during culture pose a potential risk for the successful application of stem cells. To assess the risk of in vitro expansion on mutation accumulation we have performed whole genome sequencing of clonally expanded human induced pluripotent stem cells (iPS cells) and adult stem cells (ASCs) to identify all mutations that accumulated over a fixed culture period. We find that ASCs acquire more single nucleotide variants (SNVs) and indels per population doubling than the iPS cells. When compared with ASCs, iPS cells are more vulnerable to mutations in genes and promoters. Mutational analysis revealed a clear in vitro induced mutational signature that is irrespective of stem cell type. This in vitro signature is characterized by C to G transversions that are probably caused by oxidative stress. Additionally, we observed stem cell specific mutational signatures and differences in transcriptional strand bias, indicating differential activity of DNA repair mechanisms between stem cell types in culture. In conclusion, in vitro culture induces mutation accumulation in iPS cells and ASCs. Culture under low-oxygen tension may help to reduce the number of culture-induced mutations.

  Study EGAS00001002955

• The Molecular Landscape of Asian Breast Cancers Reveals Clinically Relevant Population-Specific Differences

  Molecular profiles of breast cancer have contributed to an improved understanding of the disease, enabled development of molecular prognostic signatures to guide treatment decisions and unveiled new or more accurate therapeutic options for breast cancer patients. However, the extent to which differences in genetic, environmental and lifestyle factors influence molecular profiles in different populations remains poorly characterised, as relatively few large-scale molecular studies of breast tumours in non-Caucasian populations have hitherto been reported. Here, we present the molecular profiles of 560 Asian breast tumours and a comparative analysis of breast cancers arising in Asian and Caucasian women. Compared to the breast tumours in predominantly Caucasian women reported in TCGA and METABRIC, we show an increased prevalence of Her2-enriched molecular subtypes and higher prevalence of TP53 somatic mutations in ER+ Asian breast tumours. Using gene expression and immunohistochemistry, we observed elevated immune scores in Asian breast tumours, suggesting potential clinical response to immune checkpoint inhibitors. Whilst Her2-subtype and enriched immune score are associated with improved survival, presence of TP53 somatic mutations is associated with poorer survival in ER+ tumours. Taken together, these population differences unveil new opportunities to improve the understanding of this disease and lay the foundation for precision medicine in different populations.

  Study EGAS00001004518

• Hermansky-Pudlak syndrome type 1 causes impaired anti-microbial immunity through a pathogenic lipid metabolism-mTOR circuit - 10x Genomics scRNAseq

  Mendelian diseases that present with immune-mediated disorders can provide insights into the molecular mechanisms that drive inflammation. Hermansky-Pudlak syndrome (HPS) types 1 and 4 are caused by defective vesicle trafficking involving the BLOC-3 complex. The presence of inflammatory complications such as Crohn’s disease-like inflammation and lung fibrosis in these patients remains enigmatic. Using mass cytometry we observe an augmented inflammatory monocyte compartment in HPS1 patient peripheral blood that may be associated with a TNF - and IL-1α-dominated cytokine dysregulation. HPS1 patient monocyte-derived macrophages express an inflammatory TNF-OSM mRNA gene signature and changes in lipid metabolism. Using stimulation experiments and lysosomal proteomics we show that defective lipid metabolism drives RAB32-dependent mTOR signaling, facilitated by the accumulation of mTOR on lysosomes. This pathogenic circuit translates into aberrant bacterial clearance, which can be rescued with mTORC1 inhibition. We reveal that a pathogenic lipid-mTOR signaling circuit acts as a metabolic checkpoint for defective anti-microbial activity. This mechanism may be relevant to the complex pathology of HPS1 patients featuring macrophage lipid accumulation, granuloma formation, defective anti-microbial activity and tissue inflammation. Lastly, this circuit may be present in a wider group of disorders with defective lipid metabolism and cholesterol accumulation.

  Study EGAS00001005098

• Variation and transmission of the human gut microbiota across generations - shotgun data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005649

• Variation and transmission of the human gut microbiota across generations - 16S data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005651

• A recessive PRDM13 mutation results in congenital hypogonadotropic hypogonadism and cerebellar hypoplasia

  The PRDM13 (PR Domain containing 13) putative chromatin modifier and transcriptional regulator functions downstream of the transcription factor PTF1A, which controls GABAergic fate in the spinal cord and neurogenesis in the hypothalamus. Here, we report a novel, recessive syndrome associated with PRDM13 mutation. Patients exhibited intellectual disability, ataxia with cerebellar hypoplasia, scoliosis and delayed puberty with congenital hypogonadotropic hypogonadism (CHH). Expression studies revealed Prdm13/PRDM13 transcripts in the developing hypothalamus and cerebellum in mouse and human. An analysis of hypothalamus and cerebellum development in mice homozygous for a Prdm13 mutant allele revealed a significant reduction in the number of Kisspeptin (Kiss1) neurons in the hypothalamus and PAX2+ progenitors emerging from the cerebellar ventricular zone. The latter was accompanied by ectopic expression of the glutamatergic lineage marker TLX3. Prdm13-deficient mice displayed cerebellar hypoplasia, normal gonadal structure, but delayed pubertal onset. Together, these findings identify PRDM13 as a critical regulator of GABAergic cell fate in the cerebellum and of hypothalamic kisspeptin neuron development, providing a mechanistic explanation for the co-occurrence of CHH and cerebellar hypoplasia in this syndrome. To our knowledge, this is the first evidence linking disrupted PRDM13-mediated regulation of Kiss1 neurons to CHH in humans.

  Study EGAS00001005878

• Sclerosing epithelioid fibrosarcoma case report

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

  Study EGAS00001005215

• Comprehensive Genomic Characterization of Refractory Multiple Myeloma Reveals a Complex Mutational and Structural Landscape Associated with Drug Resistance (H067)

  Outcome of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) is still very poor. To elucidate the genetic landscape of relapsed/refractory MM (RRMM) we performed whole genome and transcriptome sequencing of 39 heavily pretreated RRMM patients with at least double refractoriness revealing complex structural changes and a high mutational load (H067/K08K). We found a median of 67 structural variants and 116 nonsynonymous exonic mutations per patient. Assessment of oncogenic drivers showed NRAS, BRAF, TP53, KRAS, SLC4A7, FAM46C, DIS3, RB1, KMT2C, MLLT4, EWSR1, HCFC2, and COPS3 as significantly mutated. We found multiple genomic regions with bi-allelic events affecting known tumor suppressor genes, most frequently TP53, RB1, and TRAF3. Supervised analysis of mutational signatures showed a major contribution of Cosmic signature 6 associated with deficiency in DNA mismatch repair and the ‘BRCAness’ signature which coincided with mutations in genes known to confer PARP inhibitor sensitivity. Recurrently mutated gene networks included genes with relevance for PI and IMiD activity, the latter including in particular members of the Cereblon and the COP9 signalosome complex. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

  Study EGAS00001004363

• tFL with a PMBL GE signature

  We investigated the clinicopathologic features of 5 follicular lymphomas (FL) that transformed to morphologic diffuse large B-cell lymphomas (DLBCL) and had a primary mediastinal large B-cell lymphoma (PMBL)-like gene expression profile. None of the transformed FL (tFL) arose in the mediastinum, all cases tested had a germinal center B-cell immunophenotype, 20% were CD30+, 60% CD23+, 80% MAL+, 20% CD200+, and 0% CD273/PDL2+. Whole exome sequencing detected alterations in genes associated with both FL/DLBCL (CREBBP, KMT2C, KMT2D, ARID1A, HIST1 members, and TNFRSF14) and PMBL (JAK-STAT pathway genes, B2M, and CD58). Copy number (CN) analysis detected gains/amplification of REL in 60% and STAT6 in 40%, gains of chromosome 16, including IL4R, in 40%, and both deletions and gains of 11q in 20%. This cohort, although limited in size, supports a subset of tFL that has blended features between FL/DLBCL and PMBL. Despite having some features that are less common in DLBCL (MAL and CD23 expression, JAK-STAT activation), these cases lack the most characteristic CN alteration seen in PMBL (9p24.1 gain/amplification). This cohort expands the biologic heterogeneity of tFL, illustrating a subset with blended FL/DLBCL and PMBL features, with potential treatment implications that include the use of novel PMBL-targeted therapies.

  Study EGAS00001005870

• The malaria-protective human glycophorin structural variant DUP4 shows somatic mosaicism and association with hemoglobin levels

  Glycophorin A and glycophorin B are red blood cell surface proteins that are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A - glycophorin B fusion gene, and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40% Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant, and reveal somatic variation in copy number of the glycophorin A- glycophorin B fusion gene. By developing a simple, specific, PCR-based assay for DUP4 we show the DUP4 variant reaches a frequency of 13% in a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation, and show that it is associated with a malarial-related phenotype in a non-hospitalizedlongitudinally-followed population.

  Study EGAS00001003239

• CRISPR/Cas9-mediated genome editing of Schistosoma mansoni acetylcholinesterase

  Here we report successful gene knock-in (KI) in the eggs of Schistosoma mansoni by combining CRISPR/Cas9 with single-stranded oligodeoxynucleotides (ssODNs). We targeted the acetylcholinesterase (AChE) gene of S. mansoni using two synthetic guide RNAs (gRNAs), X5 and X7, respectively. Liver eggs of S. mansoni were exposed to CRISPR-vector containing X5 or X7 by electroporation. Simultaneously, eggs were transfected with a ssODN donor encoding a stop codon in all six frames. Next generation sequencing analysis revealed that CRISPR/Cas9-mediated editing in S. mansoni eggs resulted in Homology-Directed Repair (HDR) when template DNA ssODN provided. Furthermore, soluble egg antigen (SEA) from AChE-modified eggs exhibited markedly reduced AChE activity compared with controls, indicative that programmed Cas9 cleavage mutated the AChE gene. Following injection of modified schistosome eggs into the tail veins of mice, a significant decrease in granuloma size in the lungs of these animals. Notably, an enhanced Th2 response induced by eggs in lung, and splenocytes small intestine-draining mesenteric lymph node cells was also generated in mice injected with X5-KI eggs in different methods. These findings further demonstrate the power and utility of CRISPR/Cas9-based genome editing for undertaking functional genomics studies in schistosomes.

  Study EGAS00001004455

• Whole Exome Sequencing PPGL

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

  Study EGAS00001006043

• Stereotyped B-cell responses are linked to IgG constant region polymorphisms in multiple sclerosis

  Clonally related B cells infiltrate the brain, meninges, and cerebrospinal fluid of multiple sclerosis (MS) patients, but the mechanisms driving the B-cell response and shaping the immunoglobulin repertoires remain unclear. Here, we used single-cell full-length RNA-seq and B-cell receptor reconstruction to simultaneously assess the phenotypes, isotypes, constant region polymorphisms, and the paired heavy- and light-chain repertoires in intrathecal B cells. We detected extensive clonal connections between the memory B cell and antibody-secreting cell compartments and observed clonally related cells of different isotypes, including IgM/IgG1, IgG1/IgA1, IgG1/IgG2, and IgM/IgA1. There was a strong dominance of the G1m1 allotype constant region polymorphisms in antibody-secreting cells, but not in memory B cells. Tightly linked to the G1m1 allotype, we found a preferential pairing of the immunoglobulin heavy-chain variable (IGHV)4 gene family with the κ variable (IGKV)1 gene family. The IGHV4-39 gene was most used and showed the highest frequency of pairing with IGKV1-5 and IGKV1(D)-33. These results link IgG constant region polymorphisms to stereotyped B-cell responses in MS and indicate that the intrathecal B-cell response in these patients could be directed against structurally similar epitopes.

  Study EGAS00001005745

• Regulatory T cell transcriptomic reprogramming characterizes adverse events by checkpoint inhibitors in solid tumors

  Immune checkpoint inhibitors (ICI), which target immune regulatory pathways to unleash antitumor responses, have revolutionized cancer immunotherapy. Despite the remarkable success of ICI immunotherapy, a significant proportion of patients whose tumors respond to these treatments develop immune-related adverse events (irAEs) resembling autoimmune diseases. Although the clinical spectrum of irAEs is well characterized, their successful management remains empiric. This is in part because the pathogenic mechanisms involved in the break-down of peripheral tolerance and induction of irAEs remain elusive. Herein, we focused on regulatory T cells (Tregs) in individuals with irAEs because these cells are vital for maintenance of peripheral tolerance, appear expanded in the peripheral blood of individuals with cancer and abundantly express checkpoint molecules, hence representing direct targets of ICI immunotherapy. Our data demonstrate an intense transcriptomic reprogramming of CD4+CD25+CD127– Tregs in the blood of individuals with advanced metastatic melanoma who develop irAEs following ICI immunotherapy, with a characteristic inflammatory, apoptotic and metabolic signature. This inflammatory signature was shared by Tregs from individuals with different types of cancer developing irAEs and individuals with autoimmune diseases. Our findings propose an inflammatory Treg reprogramming as a feature of immunotherapy-induced irAEs, that may facilitate translational approaches aiming to induce robust antitumor immunity without disturbing peripheral tolerance.

  Study EGAS00001004694

• Proteogenomics of chronic lymphocytic leukemia

  Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterized the proteome and transcriptome in-depth alongside genetic and ex-vivo drug response profiling in a clinically well annotated CLL discovery cohort (n= 68). Unsupervised clustering of the proteome data revealed six subgroups. Five of these proteomic groups were associated with genetic features, while one group was only detectable at the proteome level. This new group was characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). We developed classifiers to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n= 165, n= 169) and confirmed that ASB-CLL comprises about 20 % of CLL patients. The inferior overall survival observed in ASB-CLL was independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling. This is the transcriptomics data of the discovery cohort.

  Study EGAS00001005746

• Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas

  Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents and ranges on the spectrum of malignancy from low-grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF. However, in approximately 50% of cases, no genetic driver has been established. Clinical whole genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine−threonine kinase that activates MEK. The patient who had exhausted all other therapeutic options was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements - uncovered by comprehensive clinical sequencing of a single case - are the most common genetic event in spitzoid melanoma, are present in adult melanomas, and could be amenable to MEK inhibition.

  Study EGAS00001003430

• Androgen receptor blockade promotes response to BRAF/MEK-targeted therapy

  We studied a group of melanoma patients treated with neoadjuvant BRAF/MEK-targeted therapy (NCT02231775, n=51), and observed significantly higher rates of major pathologic response (MPR= <10% viable tumor at resection) and improved recurrence-free survival (RFS) in females versus males (MPR-66% versus 14%, p=0.001; RFS-64% versus 32% at 2 years, p=0.021). Findings were validated in a several additional cohorts 2-4 patients with unresectable metastatic melanoma treated with BRAF and/or MEK-targeted therapy (n=664 patients in total), demonstrating improved progression-free survival (PFS) and overall survival (OS) in females versus males in several of these studies. Studies in pre-clinical models demonstrated significantly impaired anti-tumor activity in male versus female BRAF/MEK-treated mice (p=0.006), with significantly higher expression of androgen receptor (AR) in tumors of male and female BRAF/MEK-treated mice versus control (p=0.0006 and 0.0025). Pharmacologic inhibition of AR signaling improved responses to BRAF/MEK-targeted therapy in male and female mice (p=0.018 and p=0.003), whereas induction of AR signaling (via testosterone administration) was associated with significantly impaired response to BRAF/MEK-targeted therapy in males and females (p=0.021 and p<0.0001).

  Study EGAS00001006196

• WXS and RNA-seq for 22 patients treated with radiation + immunotherapy

  Over 500 active clinical trials are investigating combination radiation (RT) and immune checkpoint blockade (ICB) as a cancer treatment; however, the majority of trials have found no positive interaction. We performed a comprehensive molecular analysis of a randomized phase I clinical trial of patients with non-small cell lung cancer (NSCLC) treated with concurrent or sequential ablative RT and ICB. Concurrent treatment was superior to sequential treatment in augmenting local and distant tumor responses and improving overall survival in a subset of patients with highly aneuploid tumors, but not in those with less aneuploid tumors. Our analysis suggested that concurrent treatment eliminated immunologically cold aneuploid tumors better than sequential treatment and improved patient outcomes. In addition, we report that RT alone decreased intratumoral cytotoxic effector T cell and adaptive immune gene expression signatures, in contrast to upregulation of key immune pathways after concurrent administration of RT and ICB. Our findings distinguish the differential genomic and transcriptomic effects of RT versus RT and ICB and challenge the prevailing paradigm that local ablative RT beneficially stimulates the immune response. We propose the use of tumor aneuploidy as a novel biomarker and therapeutic target in personalizing treatment approaches for patients with NSCLC treated with RT and ICB.

  Study EGAS00001006212

• Sex chromosome aneuploidies give rise to changes in the circular RNA profile

  Purpose The landscape of circular RNAs (circRNAs), an important class of non-coding RNAs that regulate gene expression, has never been described in human disorders of sex chromosome aneuploidies. We profiled circRNAs in Turner syndrome females (45,X;TS) and Klinefelter syndrome males (47,XXY; KS) to investigate how circRNAs respond to a missing or an extra X chromosome. Methods Samples of blood, muscle and fat were collected from individuals with TS (n = 33) and KS (n = 22) and from male (n = 16) and female (n = 44) controls. CircRNAs were identified using a combination of circRNA identification pipelines (CIRI2, CIRCexplorer2 and circRNA_finder). Results Differential expression of circRNAs was observed throughout the genome in TS and KS, in all tissues. The host-genes from which several of these circRNAs were derived, were associated with known phenotypic traits. Furthermore, several differentially expressed circRNAs had the potential to capture micro RNAs that targeted protein-coding genes with altered expression in TS and KS. Conclusion Sex chromosome aneuploidies introduce pervasive changes in the circRNA transcriptome, demonstrating that the genomic changes in these syndromes are more complex than hitherto thought. CircRNAs may help explain some of the genomic and phenotypic traits observed in these syndromes.

  Study EGAS00001006404

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_cell_carcinoma___2

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

  Study EGAS00001003542

• Whole genome sequencing delineates regulatory, structural, and cryptic splice variants in early onset cardiomyopathy

  Cardiomyopathy (CMP) is a heritable disorder. Over 50% cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1,953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine % cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen % harbored high-risk regulatory variants in promoters and enhancers of CMP genes (Odds ratio 2.25, p=6.70×10-7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (Odds ratio 6.7-58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early-onset CMP.

  Study EGAS00001004929

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

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

  Study EGAS00001004943

• Significant and pervasive effects of RNA degradation on Nanopore direct RNA sequencing

  Abstract Oxford Nanopore direct RNA sequencing (DRS) is capable of sequencing complete RNA molecules and accurately measuring gene and isoform expression. However, as DRS is designed to profile intact RNA, expression quantification may be more heavily dependent upon RNA integrity than alternate RNA-seq methodologies. It is currently unclear how RNA degradation impacts DRS or if it can be corrected for. To assess the impact of RNA integrity on DRS we performed a degradation time-series using SH-SY5Y neuroblastoma cells. Our results demonstrate that degradation is a significant and pervasive factor that can bias DRS measurements, including a reduction in library complexity resulting in an overrepresentation of short genes and isoforms. Degradation also biases differential expression analyses; however, we find explicit correction can almost fully recover meaningful biological signal. In addition, DRS provided less biased profiling of partially degraded samples than nanopore cDNA-PCR sequencing. Overall, we find samples with RIN > 9.5 can be treated as undegraded and samples with RIN > 7 can be utilised for DRS with appropriate correction. These results establish the suitability of DRS for a wide range of samples, including partially degraded in-vivo clinical and post-mortem samples, whilst limiting the confounding effect of degradation on expression quantification.

  Study EGAS00001006542

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

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

  Study EGAS00001004292

• T cell responses of ALS patients

  In this project, we aimed to understand how T cell responses contribute to the disease progression of amyotrophic lateral sclerosis (ALS). The present data is on single-cell sequencing isolated from human cerebrospinal fluid (CSF) cells from both ALS patients (n=5) and controls (n=4). This analysis was conducted as part of a bigger project which is summarized in the section below. We used flow cytometry to define T cell subsets and phenotypes in blood and CSF samples collected at the time of diagnosis on a cohort of 89 newly diagnosed ALS patients in Stockholm, Sweden. High frequency of CD4+FOXP3- effector T cells in blood and CSF was associated with poor survival whereas high frequency of activated regulatory T (Treg) cells and high ratio between activated and resting Treg cells in blood was associated with better survival. T cell profiles also predicted disease progression rate. On an independent cohort of cases and controls, we used single cell transcriptomics data to demonstrate that ALS patients had altered T cell gene expression patterns and clonally expanded CD4+ and CD8+ T cells in CSF. In summary, T cell responses contribute to disease progression of ALS, supporting modulation of adaptive immunity as a viable therapeutic option.

  Study EGAS00001006675

• Hexanucleotide repeat expansions in C9orf72 alter microglial responses and prevent a coordinated glial reaction in ALS

  Neuroinflammation is an important hallmark in amyotrophic lateral sclerosis (ALS). Experimental evidence has highlighted a role of microglia in the modulation of motor neuron degeneration. However, the exact contribution of microglia to both sporadic and genetic forms of ALS is still unclear. We generated single nuclei profiles of spinal cord and motor cortex from sporadic and C9orf72 ALS patients, as well as controls. We particularly focused on the transcriptomic responses of both microglia and astrocytes. We confirmed that C9orf72 is highly expressed in microglia and shows a diminished expression in carriers of the hexanucleotide repeat expansion (HRE). This resulted in an impaired response to disease, with specific deficits in phagocytic and lysosomal transcriptional pathways. Astrocytes also displayed a dysregulated response in C9orf72 ALS patients, remaining in a homeostatic state. This suggests that C9orf72 HRE alters a coordinated glial response, which ultimately would increase the risk for developing ALS. Our results indicate that C9orf72 HRE results in a selective microglial loss-of-function, likely impairing microglial-astrocyte communication and preventing a global glial response. This is relevant as it indicates that sporadic and familial forms of ALS may present a different cellular substrate, which is of great importance for patient stratification and treatment.

  Study EGAS00001006711

• T cell landscape definition by multi-omics identifies galectin-9 as novel immunotherapy target in chronic lymphocytic leukemia (CLL)

  Failure of immunotherapy after applying checkpoint inhibitors or CAR-T cells is linked to T cell exhaustion. Here, we explored the T cell landscape in chronic lymphocytic leukemia (CLL) by single-cell omics analyses of blood, bone marrow and lymph node samples of patients and spleen samples of a CLL mouse model. By single-cell RNA- sequencing, mass cytometry (CyTOF), and multiplex image analysis of tissue microarrays, we defined the spectrum of phenotypes and transcriptional programs of T cells and and their differentiation state trajectories. We identified disease-specific accumulation of distinct regulatory T cell subsets and T cells harboring an exhausted phenotype exclusively in the CLL lymph node tissue. Integration of TCR data revealed a clonal expansion of CD8+ precursor exhausted T cells, suggesting their reactivity for CLL cells. Interactome analyses identified the TIM3 ligand Galectin-9 as novel immunoregulatory molecule in CLL. Blocking of Galectin-9 in CLL-bearing mice slowed down disease development and reduced the number of TIM3 expressing T cells. Galectin-9 expression correlated with shorter survival of CLL patients. Thus, Galectin-9 contributes to immune escape in CLL and represents a novel target for immunotherapy.

  Study EGAS00001006864

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities WGS

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006865

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities RNA-Seq

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006866