13342 results for ""60+mice"+"6+weeks"+diet+study+2014+2017+DOI"

in 69.25 milliseconds.

• Induction of fetal meiotic oocytes from embryonic stem cells in cynomolgus monkeys

  Human in vitro oogenesis provides a framework for clarifying the mechanism of human oogenesis. To create its benchmark, it is vital to promote in vitro oogenesis using a model physiologically close to humans. Here, we establish a foundation for in vitro oogenesis in cynomolgus (cy) monkeys: cy female embryonic stem cells harboring one active and one inactive X chromosome (Xa and Xi, respectively) differentiate robustly into primordial germ cell-like cells, which in xenogeneic reconstituted ovaries develop efficiently into oogonia and, remarkably, further into meiotic oocytes at the zygotene stage. This differentiation entails a comprehensive epigenetic reprogramming, including Xi reprogramming, yet, as partly observed in vivo, Xa and Xi remain epigenetically asymmetric with incomplete Xi reactivation. In humans and monkeys, the Xi epigenome in pluripotent stem cells functions as an Xi-reprogramming determinant. We further show that developmental-pathway over-activations with suboptimal up-regulation of relevant meiotic genes impede in vitro meiotic progression. Cy in vitro oogenesis exhibits critical homology with the human system, including with respect to bottlenecks, providing a salient model for advancing human in vitro oogenesis.

  Study JGAS000573

• MATISSE WES and bulk RNA-sequencing data

  High-risk cutaneous squamous cell carcinoma is frequently seen in frail, elderly patients and often requires mutilating surgical treatment and adjuvant radiotherapy (RT). The MATISSE trial (NCT04620200) was conducted to improve clinical prospects in CSCC patients eligible for extensive curative surgery, and evaluated early treatment response and immune dynamics following ultra-short neoadjuvant immunotherapy prior to standard of care surgery with or without radiotherapy. Fifty patients were randomized to receive two courses of nivolumab 3 mg/kg (NIVO, weeks 0&2) or nivolumab 3 mg/kg (weeks 0&2) and one course of ipilimumab 1 mg/kg (week 0 only, NIVO/IPI before surgery at week 4, with the primary objective being histopathological response based on residual viable tumor cells in the surgical specimen. Ten patients choose not to undergo surgery, of whom nine reached a clinical complete response during follow up, and were disease-free at median 34 months follow up. Serial tumor and blood samples at baseline, week 1, week 2 and week 4, enabled to investigate dynamical changes through bulk WES and RNA sequencing to potentially aid future de-escalation trials.

  Study EGAS50000001003

• Non-small cell lung cancer molecular subtypes and vulnerability to immunotherapy treatment combinations

  The phase 3 IMpower150 trial in treatment-naïve patients with metastatic non-small-cell lung cancer (NSCLC) demonstrated significantly longer progression-free (PFS) and overall survival (OS) with first-line atezolizumab (anti-PD-L1)-bevacizumab (anti-VEGF)-carboplatin-paclitaxel (ABCP) than with bevacizumab-carboplatin-paclitaxel (BCP). We characterised four molecular NSCLC subtypes identified by unsupervised clustering of transcriptomes of 564 pre-treatment primary tumour samples from IMpower150 using non-negative matrix factorization (NMF1-4). Each subtype had distinct tumour PD-L1 expression levels, epithelial characteristics, immune composition, and treatment outcomes. Both NMF2 (enriched in tumour proliferation signal, macrophages, and monocytes) and NMF4 (enriched in B-cells and T cells) had elevated tumour PD-L1 expression. Of these two, only NMF4 demonstrated PFS and OS benefits with ABCP versus either BCP or atezolizumab-carboplatin-paclitaxel (ACP). Patients with NMF1 (enriched in basal and squamous-like cells) had improved outcomes on ABCP compared with ACP or BCP; those with NMF3 (enriched in adenocarcinoma signatures) showed similar outcomes among treatments. These insights could help inform individualised first-line treatment for metastatic NSCLC.

  Study EGAS50000001272

• Dynamics of checkpoint receptors in γδ T cell subsets are associated with clinical responses during anti-PD-1 immunotherapies

  Gamma delta (γδ) T cells are innate-like lymphocytes with potent anti-tumor properties. Herein, we show that immune checkpoint receptors (ICRs) display differential expression and regulation by the JAK-STAT pathway in Vδ1 and Vδ2 cells and identify constitutive (e.g. TIGIT, PD-1) and inducible (e.g. TIM-3, LAG-3, CTLA-4) ICRs. In melanoma, all γδ T cell subsets downregulated AP-1 transcription factors, but Vδ1 cells specifically expressed high levels of ICR, TOX and inhibitory killer Ig-like receptor (KIR) transcripts, reminiscent of exhaustion. However, patient-derived cells were functionally competent, although induction of LAG-3 and CTLA-4 was impaired. During anti-PD-1 monotherapy, Vδ1 cells specifically bound high levels of therapeutic antibody but only in patients who responded to treatment, revealing a potential new prognostic marker for evaluating the efficacy of IC blockade (ICB) therapy. Finally, expression of KIR genes in Vδ1 cells was downregulated in response to successful ICB therapy. Collectively, our data indicate an intricate relationship between ICRs and γδ T cells and reveal novel approaches by which these cells can be harnessed in order to discern or improve cancer immunotherapy.

  Study EGAS50000001270

• FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  CRISPR/Cas9 FLT3 knock-out (FLT3-KO) in human HSCs and LSCs followed by functional xenograft assays to test their ability to regenerate human hematopoiesis and leukemia, respectively. FLT3-KO in HSCs from human fetal liver, cord blood and adult bone marrow showed no impairment in multilineage hematopoiesis in primary and secondary xenografts. By contrast, FLT3-KO LSCs from 6 of 7 FLT3-ITD mutated AMLs were able to generate short-term engraftment but were completely exhausted by 12 weeks. Thus, FLT3 is essential for LSC long-term propagation. This dependency was unique to FLT3-ITD AML samples as non-FLT3-ITD AML samples generated leukemic grafts upon FLT3-KO. Transcriptomic analysis revealed that FLT3-KO induced downregulation of DNA repair and cell cycle checkpoints, uniquely in FLT3-ITD AML, but not in healthy HSCs or other AMLs. Our research highlights a critical distinction between healthy HSCs and LSCs: whereas healthy hematopoiesis proceeds unperturbed upon FLT3-KO, FLT3-ITD leukemogenesis is impaired through elimination of LSCs. This evidence underscores the necessity for more potent FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  Study EGAS50000000437

• Non-canonical NF-κB signaling skews B cells away from germinal center to low-affinity effector fate

  The non-canonical (NC) NF-κB pathway is a central regulator of B-cell development and survival, yet its role during B-cell mediated adaptive responses remains uncharted. Here, we uncover an unanticipated regulatory function for NC NF-κB in shaping humoral immunity. High-resolution mapping of NC NF-kB status in B cells from human lymphoid tissues revealed low activity in germinal center B-cells (GCB) than effector subsets, suggesting a selective attenuation. Consistently, mice model harboring B-cells with enforced NC NF-κB activation revealed near- complete loss of GCB and high-affinity antibodies, following vaccination. GCB loss was not driven by intrinsic defects in formation or expansion, but by an immune-niche dysregulation marked by regulatory T-cell expansion. Mechanistically, NC NF-κB activation promoted B-cell- derived IL-10, and blocking this axis normalize the niche and rescue GCB. Concurrently, NC NF- κB signaling boosted low-affinity antibody responses by transcriptionally priming follicular B- cells for plasmacytic differentiation. Together, these findings position NC NF-κB as a tunable rheostat linking B-cell signaling to immune-niche control and humoral immunity.

  Study EGAS50000001646

• Adaptive Therapy Exploits Fitness Deficits in Chemotherapy-Resistant Ovarian Cancer to Achieve Long-Term Tumor Control Open Access

  Drug resistance results in poor outcomes for cancer patients. Adaptive therapy is a potential strategy to address drug resistance that exploits competitive interactions between sensitive and resistant subclones. Here, we showed that adapting carboplatin dose according to tumor response (adaptive therapy) significantly prolonged survival of murine ovarian cancer models compared to standard carboplatin dosing, without increasing mean daily drug dose or toxicity. Platinum resistant ovarian cancer cells exhibited diminished fitness when drug was absent in vitro and in vivo, which caused selective decline of resistant populations due to reduced proliferation and increased apoptosis. Conversely fitter, sensitive cells re-grew when drug was withdrawn. Using a bioinformatics pipeline that exploits copy number changes to quantify the emergence of treatment resistance, analysis of cell‐free DNA obtained longitudinally from ovarian cancer patients during treatment showed subclonal selection through therapy, and measurements of resistant population growth correlated strongly with disease burden. These preclinical findings pave the way for future clinical testing of personalized adaptive therapy regimens tailored to the evolution of carboplatin resistance in individual ovarian cancer patients.

  Study EGAS50000001142

• Whole Genome Sequencing of Neuroblastoma

  Neuroblastoma is a childhood tumour of the peripheral sympathetic nervous system. The pathogenesis has for a long time been quite enigmatic, as only very few gene defects were identified in this often lethal tumour1. Frequently detected gene alterations are limited to MYCN amplification (20%) and ALK activations (7%)2-5. Here we present a whole genome sequence analysis of 87 neuroblastoma of all stages. Only few recurrent amino-acid changing mutations were found. In contrast, analysis of structural defects identified in 18% of high stage neuroblastoma a local shredding of chromosomes, known as chromothripsis6. These tumours are associated with a poor outcome. Structural alterations recurrently affected ODZ3, PTPRD and CSMD1, which are involved in neuronal growth cone stabilization7-9. In addition, ATRX, TIAM1 and a series of regulators of the Rac Rho pathway were mutated, further implicating defects in neuritogenesis in neuroblastoma. Most tumors with defects in these genes were aggressive high stage neuroblastoma, but did not carry MYCN amplifications. The genomic landscape of neuroblastoma therefore reveals two novel molecular defects, chromotripsis and neuritogenesis gene alterations, which frequently occur in high risk tumors.

  Study EGAS00001000222

• Novel regional age-associated DNA methylation changes within human common disease-associated loci

  Background: Advancing age progressively impacts on risk and severity of chronic disease. It also modifies the epigenome, with changes in DNA methylation, due to both random drift and variation within specific functional loci. Results: In a discovery set of 2,238 peripheral-blood genome-wide DNA methylomes (MeDIP-seq) aged 19-82 years, we identified 71 age-associated Differentially Methylated Regions (a-DMRs, p < 1.85 x 10-8) within the Linkage Disequilibrium (LD) blocks of the NIH Catalogue of published GWAS SNPs. This included 52 novel regions, 29 within loci not covered by 450k or 27k Illumina array, and with marked enrichment for Poised Promoters and Enhancers across multiple cell types. In a replication set of 2,084 DNA methylomes, 95.7% of the a-DMRs showed the same direction of ageing effect, with 80.3% and 53.3% replicated to p < 0.05 and p < 1.85 x 10-8, respectively. Conclusion: By analysing the functionally enriched disease and trait-associated regions of the human genome, we identified novel epigenetic ageing changes, which could be useful biomarkers or provide mechanistic insights into age-related common diseases.

  Study EGAS00001001910

• Comprehensive molecular profiling identifies novel genetic drivers and subtypes underlying medulloblastoma

  Current therapies for medulloblastoma (MB), a highly malignant childhood brain tumor, impose debilitating effects on the developing child, warranting the development of molecularly targeted treatments with reduced toxicities. Prior studies have failed to disclose the full spectrum of driver genes and molecular processes operative in MB or adequately explain heterogeneity among molecular subgroups. Herein, we detail the somatic landscape across 491 sequenced MBs and molecular heterogeneity amongst 1,256 epigenetically analyzed cases, identifying subgroup-specific mutational signatures, driver genes, and pathway alterations including previously unappreciated actionable therapeutic targets. Integrative approaches for assigning driver genes to subgroups explained >65-70% of Group 3 and Group 4, doubling previous knowledge. Novel subtypes underlying Group 3 and Group 4 were differentially enriched for specific driver events, including hotspot in-frame indels targeting KBTBD4 and ‘enhancer hijacking’ driving activation of PRDM6. Thus, application of integrative genomics to an unprecedented cohort of clinical samples derived from a single childhood cancer entity has disclosed a series of new cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for treating MB.

  Study EGAS00001001953

• Field_effect_of_healthy_and_diseased_livers_WGS

  Recent work in the Campbell group has revealed somatic mutations present in normal, non-cancerous human skin. A subset of the mutations conferred selective advantages to the host cells, leading to clonal expansions and raising the risk for future cancer development. Capturing such somatic mutations in normal tissue is important to advance our understanding about carcinogenesis and could provide prospective medical insights. In this project, our goal is to detect somatic mutations in normal (pre-cancerous) liver tissue. Using Laser Microdissection technology, we will dissect individual liver lobules from patient samples and submit these to sequencing. For each patient sample, we aim to sequence multiple lobules to characterise the mutagenic burden. Samples will be taken from patients with different liver disease aetiologies, including alcoholism and obesity, with a view on distinguishing the prevalent mutation types occurring in each disease context. We will perform both whole genome and targeted sequencing, initially using the WTSI cancer panel. Later we aim to use a novel bait set that captures both cancer genes as well as genes relevant to the non-cancerous samples (ie. genes implicated in hereditary disorders, immune sequences).

  Study EGAS00001002413

• Duplexseq_of_the_interstrand_crosslinks_WGS

  One of the most dangerous forms of DNA damage are interstrand crosslinks (ICLs), which covalently crosslink the two strands of the DNA double helix. The repair of these lesions is crucial for cellular survival due to their ability to block transcription and DNA replication. Initially, the major pathway that has been described in ICL repair involves a network of 22 genes that are mutated in a severe human genetic disease known as Fanconi Anemia (FA). Using synthetic lethality screens in the near-haploid human HAP1 cell line, we recently identified two potentially novel regulators of ICL repair, C1orf112 and THAP12. Loss of C1orf112 and THAP12 causes hypersensitivity to ICL-inducing DNA damaging agents, such as Mitomycin C (MMC). Additionally, C1orf112-depleted cells show elevated levels of micronuclei and accumulation of DNA damage in S-phase. To better understand how C1orf112 and THAP12 mediate the repair of ICLs, we want to perform mutational signature analysis, using the BotSeq method. Therefore, WT, C1orf112 and THAP12 knockout cells were cultured in vehicle or MMC treated conditions for 10 days and the genomic DNA was isolated. FANCA and FANCD2 knockout cells are taken along as controls in this experimental setting.

  Study EGAS00001006545

• SHANK2 mutations associated with Autism Spectrum Disorder cause hyperconnectivity of human neurons

  Heterozygous loss-of-function mutations in the synaptic scaffolding gene SHANK2 are strongly associated with autism spectrum disorder (ASD). To investigate their effect on synaptic connectivity, we generated cortical neurons from induced pluripotent stem cells (iPSC) derived from neurotypic and ASD-affected donors. We developed Sparse coculture for Connectivity (SparCon) assays where SHANK2 and control neurons were differentially labeled and sparsely seeded together on a lawn of unlabeled control neurons. We observed striking increases in dendrite length, dendrite complexity, total synapse number, and frequency of spontaneous excitatory postsynaptic currents. These findings were phenocopied in gene-edited homozygous SHANK2 knockout cells and rescued by gene correction of an ASD SHANK2 muation, supporting a role for SHANK2 as a regulator of connectivity in developing human neurons. Dendrite length increases were exacerbated by IGF1, TG003, or BDNF, and suppressed by DHPG treatment. The transcriptome in these isogenic SHANK2 neurons was deeply perturbed in synapse, plasticity, and neuronal morphogenesis gene sets and ASD gene modules, and activity-dependent dendrite extension was impaired. Our unexpected findings provide evidence for hyperconnectivity and profoundly altered transcriptome in SHANK2 neurons derived from ASD subjects.

  Study EGAS00001003436

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007602

• Distinct portrayal of lesions in synchronous multifocal lung adenocarcinoma revealed by genome sequencing

  Distinguishing multiple primary lung cancers in the synchronous multifocal intrapulmonary lesions has important significance on clinical staging and therapeutic decision. To investigate genomic aberration profiles, we applied whole genome and whole exome sequencing, and microarray-based comparative genomic hybridization on 15 intrapulmonary tumors derived from six patients with synchronous multifocal lung cancers having similar histological diagnosis. Any pair of intrapulmonary tumors in a single patient, which shared the identical genetic background and environment, showed an extinctive heterogeneity between each other. Phylogenetic relationship analysis indicated an independently branched evolution among all the tumors, suggesting they were multiple primary lung cancers. EGFR or KRAS mutations were found in 7 or 3 out of the 15 tumors, from 3 or 2 patients, respectively. Somatic mutational heterogeneity of these two genes in a single patient was also observed. Our analysis indicates genomic aberration profiling is valuable for identification of multiple primary lung cancer, especially when high histopathological concordance was observed between lesions. We also suggest a thoroughly molecular diagnosis against therapeutic target genes should be taken for each accessible nodule before making a plan for adjuvant therapy.

  Study EGAS00001001572

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

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

  Study EGAS00001005359

• Transcriptomes of human CD4+ T lymphocytes - Metabolic project

  Human CD4+ T cells are essential mediators of immune responses. By altering the mitochondrial and metabolic states, we defined metabolic requirements of human CD4+ T cells for in vitro activation, expansion, and effector function. T cell activation and proliferation were reduced by inhibiting oxidative phosphorylation, while early cytokine production was maintained by either OXPHOS or glycolytic activity. Glucose deprivation in the presence of mild mitochondrial stress markedly reduced all three T cell functions, contrasting the exposure to resveratrol, an antioxidant and sirtuin-1 activator, which specifically inhibited cytokine production and T cell proliferation, but not T cell activation. Conditions that inhibited T cell activation were associated with the downregulation of 2′,5′-oligoadenylate synthetase genes via interferon response pathways. Our findings indicate that T cell function is grossly impaired by stressors combined with nutrient deprivation, suggesting that correcting nutrient availability, metabolic stress and/or the function of T cells in these conditions will improve the efficacy of T cell-based therapies.

  Study EGAS00001005565

• Single-nucleus transcriptomic profiling of aging Down Syndrome brains

  Down syndrome (DS) – trisomy of human chromosome 21 (HSA21) – is characterized by lifelong cognitive impairments and development of neuropathological hallmarks of Alzheimer’s disease (AD). The cellular and molecular modifications responsible for these effects are not understood. Here we performed single-nucleus RNA-sequencing (snRNA-seq) employing both short (Illumina) and long-read (Pacific Biosciences) sequencing technologies on a total of 29 DS and non-DS control prefrontal cortex samples. In DS, the ratio of inhibitory-to-excitatory neurons was significantly increased, which was not observed in previous reports examining sporadic AD. DS microglial transcriptomes displayed AD-related aging and activation signatures in advance of AD neuropathology, with increased microglial expression of C1q complement genes (associated with dendritic pruning) and the HSA21 transcription factor gene RUNX1. Long-read sequencing detected vast RNA isoform diversity within and amongst specific cell types including numerous novel sequences that differed between DS and normal brains. Notably, over 8,000 genes produced RNAs containing intra-exonic junctions (IEJs), including APP that had previously been associated with somatic gene recombination. These and related results illuminate large-scale cellular and transcriptomic alterations as novel features of the aging DS brain.

  Study EGAS00001005691

• Cross-tissue transcriptomic analysis of human secondary lymphoid organ residing ILC3 reveals a default quiescent state in the absence of inflammation

  A cross tissue transcriptional comparison of human ILC3 isolated from non-inflamed lymph nodes (LN) and spleens, inflamed tonsils and the peripheral blood (PB). Using uniform gating and high-purity cell-sorting, ILC3 were sorted as Lineage(CD3,CD19,CD34,CD14,CD94)-CD45intCD127+CRTH2-CD117+. From LN, spleen, tonsils and PB, NKp44- ILC3 were isolated and, in addition, NKp44+ ILC3 were sorted from spleen and tonsils. RNA quality was measured on a RNA 6000 Pico Kit (Agilent Technologies) using a 2100 Bioanalyzer (Agilent Technologies). 100 pg high-quality RNA was isolated and cDNA was generated using the SMART-Seq (v3) Ultra Low Input RNA Kit (Clontech) with 15 cycles of amplification. Amplified and Covaris-fragmented cDNA was quality-checked on a High Sensitivity DNA chip (Agilent Technologies) and further processed according the TruSeq Nano DNA Library Preparation Kit (Illumina). Data was processed as follows: SMARTer adapters were trimmed with cutadapt and the resulting sequences were aligned to the human RefSeq transcriptome using 14 TopHat2 (Kim et al., 2013). Normalization and quantification was performed using Cufflinks (Trapnell et al., 2012). FPKM counts were determined per gene with HTSeq-count.

  Study EGAS00001002636

• The identification of genetic vulnerabilities in head and neck cancers for the development of novel treatments.

  The global cancer burden is expected to double by the year 2030 and this will disproportionately affect the Asian continent. Head and neck squamous cell carcinoma (HNSCC) is endemic in Asia. It is the sixth most common cancer worldwide and the second most common in Malaysia, with approximately 600,000 new cases per year globally. Longterm survival remains less than 50%, underscoring the need for better therapeutic strategies. In this proposal we aim to bring new technological advancements in the identification of critical genes in HNSCC and the identification of compounds for the development of new therapies in HNSCC. We propose to genomically characterise 16 unique Asian H&N cancer cell lines, to classify sub groups where identification of essential genes that can be targeted for tumor control can be applied in a more rational way. Further, we propose to use zebrafish, an established drug discovery platform to identify promising pathway-specific compounds that could be further studied using CRISPR/Cas9 lethality screens and biological models to better identify novel therapeutic strategies to combat this major health problem in Asia.

  Study EGAS00001002682

• UK10K_OBESITY_GS

  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 GS samples are part of the Obesity group and will undergo exome sequencing.

  Study EGAS00001000242

• Comprehensive characterization of cell-free tumor DNA in plasma and urine of patients with renal tumors

  Cell-free tumour-derived DNA (ctDNA) allows non-invasive monitoring of various cancers but its utility in renal cell cancer (RCC) has not been well established. A combination of untargeted and targeted sequencing methods, applied to two independent cohorts of patients (n=90) with various renal tumour subtypes, were used to determine ctDNA content in plasma and urine. Our data revealed lower plasma ctDNA levels in RCC relative to other cancers, with untargeted detection of ~33% for both cohorts. A highly sensitive personalised approach, applied to plasma and urine from select patients (n=22), improved detection to ~50%, including in patients with early stage and even benign lesions. A machine-learning based model, applied to untargeted data, predicted this detection, potentially offering a means of triaging patient samples for personalised analysis. We observed that plasma, and for the first time, urine ctDNA may better represent tumour heterogeneity than a single tissue biopsy. Longitudinal sampling of >200 plasma samples revealed that ctDNA can track disease course. These data highlight low ctDNA levels in RCC but indicate potential clinical utility provided improvement in isolation and detection approaches.

  Study EGAS00001003530

• Genome-wide cell-free DNA fragmentation in patients with cancer

  Cell-free DNA (cfDNA) in the blood provides a noninvasive diagnostic avenue for patients with cancer. However, characteristics of the origins and molecular features of cfDNA are poorly understood. We developed an approach to evaluate fragmentation patterns of cfDNA across the genome and found that cfDNA profiles of healthy individuals reflected nucleosomal patterns of white blood cells, while patients with cancer had altered fragmentation profiles. We applied this method to analyze fragmentation profiles of 236 patients with breast, colorectal, lung, ovarian, pancreatic, gastric, or bile duct cancers and 245 healthy individuals. A machine learning model incorporating genome-wide fragmentation features had sensitivities of detection ranging from 57% to >99% among the seven cancer types at 98% specificity, with an overall AUC of 0.94. Fragmentation profiles could be used to identify the tissue of origin of the cancers to a limited number of sites in 75% of cases. Combining our approach with mutation based cfDNA analyses detected 91% of cancer patients. The results of these analyses highlight important properties of cfDNA and provide a proof of principle approach for screening, early detection, and monitoring of human cancer.

  Study EGAS00001003611

• Genomic analysis of Smoothened inhibitor resistance in basal cell carcinoma

  Smoothened inhibitors are currently being investigated for the treatment of several cancers. Vismodegib is approved for the treatment of locally advanced and metastatic basal cell carcinoma (BCC). The majority of BCC patients treated with vismodegib experience significant clinical benefit, however, a small number of patients develop resistance. Knowledge of resistance mechanisms can generate predictive biomarkers and is critical for the design of additional therapeutic strategies aimed at circumventing resistance. To investigate mechanisms of resistance to vismodegib in BCC we sequenced and profiled biopsies from patients who initially responded to treatment and subsequently progressed on drug.We find that resistance is associated with re-activation of the Hedgehog (Hh) pathway through diverse mechanisms, including concurrent copy number alterations in the downstream Hh pathway components SUFU and GLI2, as well as novel SMO mutations that reduce sensitivity to vismodegib. The latter fall into two classes: those that activate SMO and those that directly affect the vismodegib-binding pocket. Additionally, we observe evidence of clonal selection of SMO mutations and intra-tumor heterogeneity of resistance mechanisms, which implies diverse and combinatorial strategies are required to overcome resistance.

  Study EGAS00001000845

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

  To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13,037 individuals enrolled in the NIHR BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor encoding gene IKZF5 and thrombocytopenia. We report five causal missense variants in or near IKZF5 zinc fingers (Znfs), of which two occurred de novo and three co-segregated in three pedigrees. A canonical DNA-Znf binding model predicts that three of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared to wild-type (WT) IKZF5 and electron microscopy revealed a reduced quantity of alpha granules in normally sized platelets. Proplatelet formation (PPF) was reduced in megakaryocytes (MKs) from seven cases relative to six controls. Comparison of RNA-seq data from platelets, monocytes, neutrophils and CD4 T-cells from three cases and 14 healthy controls showed 1,194 differentially expressed genes in platelets but only four DEGs in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

  Study EGAS00001003736

• RNA sequencing and Illumina 2.5M SNP array data collected from 675 commonly used human cancer cell lines.

  Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic response1. Although substantial effort has been directed to defining the genomic constitution of cancer cell line panels2–4, the transcriptome – which represents the active program of a cell – remains understudied. Here, we describe RNA sequencing and SNP array analysis of 675 commonly used human cancer cell lines. We explore numerous transcriptome features including coding and non-coding gene expression, transcribed mutations, gene fusion and expression of non-human sequences. Aside from many known aberrations we find new surprising characteristics, including more than 2200 unique fusion gene pairs representing a vast, testable repertoire of oncogenic fusions, many of which have analogs found in primary human tumors. We show that a combination of multiple genome and transcriptome features in a novel pathway-based approach enhances prediction of response to various targeted therapeutics. Our results provide valuable new insights into these critical pre-clinical models and provide added context for interpreting the numerous studies that employ these widely used cell lines.

  Study EGAS00001000610

• Reference Exome Data for a Northern Brazilian population

  Whole exome sequencing (WES) is a popular technology widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. However, there is not always adequate population-specific reference data available to assist clinical geneticists in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from the State of Rio Grande do Nord in Brazil. Sequence data were processed using an ‘intersect-then-combine’ (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 439 (0.755%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil.

  Study EGAS00001004112

• Oncoprint GSCCs

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007481

• Spatial_transcriptome_analysis_of_Paediatric_Thymus

  Cell atlas of Paediatric Thymus Cell atlas of Paediatric Thymus - The Human Cell Atlas initiative is an ambitious global research effort which aims to describe every cell in the human body across all developmental stages to generate a reference map to accelerate progress in science and medicine particularly relating to human development and disease. The initiative is led by investigators based at the Wellcome Sanger Institute and the Broad Institute of MIT and Harvard and funding support from the Wellcome and NIH. We aim to interrogate the cellular composition and molecular regulators underpinning a cells identity through; genetic and genomic analysis, expression analysis, non-identifiable video imaging and cell line and organoid production, as well as, computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequently validate these findings in situ using spatial genomic and imaging technologies e.g. immunohistochemistry, RNAscope, in situ sequencing and spatial transcriptomic analyses. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004281

• Gut microbiome modulates response to anti PD1 immunotherapy in metastatic melanoma patients

  There is a growing appreciation of the role of the microbiome in cancer, and evidence in pre-clinical models that the gut microbiome may modulate responses to immune checkpoint blockade though this has not been well-characterized in patients. We analyzed the oral (n=86)and gut (n=43) 16S microbiome in melanoma patients on PD-1 blockade. Significant differences were noted in the diversity and composition of the gut microbiome between responders and non-responders in patients with a fecal microbiome sample, with significantly higher alpha diversity and relative abundance of Ruminococcaceae bacteria) in R. Metagenomic studies (n=25) revealed functional differences in gut bacteria in R including enrichment of anabolic pathways. Immune profiling demonstrated enhanced systemic and anti-tumor immunity in patients with a favorable gut microbiome as well as in germ-free mice receiving fecal transplant from responding patients. Together, these data have important implications for treatment with immune checkpoint blockade in cancer. The fastq files associated with this dataset are stored at ENA under the following accession numbers: PRJEB22894 (Fecal 16S), PRJEB22874 (Oral 16S), PRJEB22895 (Murine 16S), PRJEB22893 (Fecal WGS).

  Study EGAS00001002698

• Germline HAVCR2 mutations altering TIM-3 characterize subcutaneous panniculitis-like T-cell lymphomas with hemophagocytic lymphohistiocytic syndrome

  Sub-cutaneous panniculitis-like T-cell lymphomas (SPTCL), a non-Hodgkin lymphoma, can be associated with hemophagocytic lymphohistiocytosis (HLH), a life-threatening immune activation which adversely impacts survival1,2. T-cell-immunoglobulin mucin-3 (TIM-3) is a modulator of immune responses expressed on subgroups of T- and innate immune cells. We identify in ~60% of SPTCL cases germline, loss-of-function, missense variants in highly conserved residues of TIM-3, c.245A>G(p.Y82C) and c.291A>G(p.I97M), each with specific geographic distribution. Y82C-TIM-3 occurs on a potential founder chromosome in patients with East Asian and Polynesian ancestry, while I97M-TIM-3 occurs in Caucasians. Both variants induce protein misfolding and abrogate TIM-3’s plasma-membrane expression leading to persistent immune activation, increased production of inflammatory cytokines, including TNF-α and IL-1β, promoting HLH and SPTCL. Our findings highlight HLH/SPTCL as a new genetic entity and showcase TIM-3 mutations as the first causative genetic defect in SPTCL. While TIM-3-mutant HLH/SPTCL benefit from immunomodulation, therapeutic repression of the TIM-3 checkpoint may have adverse consequences.

  Study EGAS00001002765

• Whole genome sequencing of Ewings Sarcoma

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003385

• Onco-exaptation of an Endogenous Retroviral LTR Drives IRF5 Expression in Hodgkin Lymphoma

  The transcription factor interferon regulatory factor 5 (IRF5) is upregulated in Hodgkin lymphoma (HL) and a key regulator of the aberrant transcriptome characteristic of this disease. Here we show that IRF5 upregulation in HL is driven by transcriptional activation of a normally dormant endogenous retroviral LOR1a long terminal repeat (LTR) upstream of IRF5. Specifically, through screening of RNA-seq libraries, we detected LTR-IRF5 chimeric transcripts in multiple HL cell lines but not in normal B-cell controls. In HL, the LTR was in an open and hypomethylated epigenetic state and we further show the LTR is the site of transcriptional initiation. Among HL cell lines, usage of the LTR promoter strongly correlates with overall levels of IRF5 mRNA and protein, indicating that LTR transcriptional awakening is a major contributor to IRF5 upregulation in HL. Taken together, oncogenic IRF5 overexpression in HL is the result of a specific LTR transcriptional activation. We propose that such LTR derepression is a distinct mechanism of oncogene activation (“onco-exaptation”) and that such a mechanism warrants further investigation in molecular and cancer research.

  Study EGAS00001001205

• Dissection of the molecular complexity of colorectal cancer in pre-clinical models identifies predictive signatures of sensitivity to EGFR inhibitors

  Colorectal carcinoma (CRC) is a heterogeneous tumour entity, with only a fraction of the tumours responding to currently available therapies. A detailed understanding of the biology and genetics of CRC is a prerequisite for precision medicine. To address this challenge in clinical oncology, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed an integrated pre-clinical platform analysing drug sensitivity in the established matched patient-derived in vivo and in vitro models, 46 three-dimensional cell cultures and 59 xenografts, recapitulating many of the molecular features of the donor tumours. We generated extensive omics information on patient tumours and their derived models used to generate a compendium of drug sensitivity data. We could distinguish tumour molecular sub-groups responding differently to therapy and identified a novel independent predictor of sensitivity to the EGFR inhibitor cetuximab. The findings reported here will advance our understanding of CRC and its response to therapy.

  Study EGAS00001001752

• Whole genome and transcriptome analysis of medullary thyroid cancer

  Medullary thyroid cancer is the neoplasm of the calcitonin-producing parafollicular cells of the thyroid gland. The hormone calcitonin acts to lower the blood calcium level, counteracting the effects of the parathyroid hormone; it can also serve as a biomarker for the presence or recurrence of medullary thyroid cancer. Surgery is the only curative treatment for this cancer. Targeted therapeutics vandetanib and cabozantinib are approved for the treatment of aggressive and metastatic tumors that are not amenable to surgery. The use of these multikinase inhibitors are supported by the observed overactivation of the RET oncoprotein in a large sub-population of medullary thyroid cancers; however, not all patients carry oncogenic alterations of this kinase. Hence, there is a need for molecular characterization of medullary tumors on the genomic scale with the aim of identifying targetable mutations. Here, we describe the genomic profiles of two medullary thyroid cancers and report the presence of a putative oncogenic BRAF fusion in one. Such alterations previously observed in other malignancies and known targets of available drugs, can benefit patients who currently have no therapy options.

  Study EGAS00001001473

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007604

• The demographic history and mutational load of African hunter-gatherers and farmers

  The distribution of deleterious genetic variation across human populations is a key issue in evolutionary biology and medical genetics. However, the impact of different modes of subsistence on recent changes in population size, patterns of gene flow, and deleterious mutational load remains unclear. Here, we report high-coverage exomes from various populations of rainforest hunter-gatherers and farmers from central Africa. We find that the recent demographic histories of hunter-gatherers and farmers differed considerably, with population collapses for hunter-gatherers and expansions for farmers, accompanied by increased gene flow. We show that purifying selection against newly arising deleterious alleles is of similar efficiency across African populations, in contrast with Europeans where we detect weaker purifying selection. Furthermore, the per-individual mutation load of rainforest hunter-gatherers is similar to that of farmers, under both additive and recessive models. Our results indicate that differences in the cultural practices and demographic regimes of African populations have not resulted in large differences in mutational burden, and highlight the beneficial role of gene flow in reshaping the distribution of deleterious genetic variation across human populations.

  Study EGAS00001002457

• Evolutionary analysis of primary tumors and metastatic lesions from 20 breast cancer patients (99 samples in total) using exome sequencing data.

  Metastatic breast cancers are still incurable. Characterizing their evolutionary landscape including the role of metastatic axillary lymph nodes to seed distant organ metastasis can provide rational basis for effective treatments. Here, we describe the genomic analyses of the primary tumors and metastatic lesions from 20 breast cancer patients (99 samples). Our evolutionary analyses revealed diverse spreading and seeding patterns governing tumor progression. Although linear evolution to successive metastatic sites was common, parallel evolution from the primary tumor to multiple distant sites was also evident. Metastatic spreading was frequently coupled with polyclonal seeding, where multiple metastatic subclones originated from primary tumor and/or other distant metastases. Synchronous axillary lymph node metastasis, a well-established prognosticator of breast cancer, was not involved in seeding distant metastasis, suggesting haematogenous route for cancer dissemination. Clonal evolution coincided frequently with emerging driver alterations and evolving mutational processes, notably a significant increase in apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) associated mutagenesis. Our data provide the first genomic evidence regarding the role of axillary lymph node metastasis in seeding distant organ metastasis and elucidates the evolving mutational landscape during cancer progression.

  Study EGAS00001002737

• Defective T-cell expansion in RASGRP1 deficiency

  Inherited CTPS1- CD27- and CD70-deficiencies in humans have revealed key factors of T-lymphocyte expansion, that is a prerequisite for an efficient immunity to Epstein-Barr virus (EBV) infection. RASGRP1 is a T-lymphocyte specific nucleotide exchange factor known to activate the MAP kinases pathway. A deleterious homozygous mutation in RASGRP1 leading to the loss RASGRP1 expression was identified in two siblings who both developed a persistent EBV infection leading to Hodgkin lymphoma. RASGRP1-deficient T cells exhibited defective MAPK activation and impaired proliferation that was restored by expression of wild-type RASGRP1. Similar defects were observed in T cells from healthy individuals when RASGRP1 was downregulated. RASGRP1-deficient T cells also exhibited decreased CD27-dependent proliferation toward CD70-expressing EBV-transformed B cells, a crucial pathway required for expansion of antigen-specific T cells in anti-EBV immunity. Furthermore, RASGRP1-deficient T cells failed to upregulate CTPS1, an important enzyme required for DNA synthesis. These results show that RASGRP1 deficiency leads to susceptibility to EBV infection, and demonstrate the key role of RASGRP1 at the crossroad of pathways required for the expansion of activated T lymphocytes.

  Study EGAS00001002753

• Somatic IL4R Hotspot Mutations in Primary Mediastinal Large B-cell lymphoma lead to constitutive JAK-STAT activation

  Primary mediastinal large B cell lymphoma (PMBCL) is a distinct subtype of diffuse large B cell lymphoma thought to arise from thymic medullary B cells. Gene mutations underlying the molecular pathogenesis of the disease are incompletely characterized. Here, we describe novel somatic IL4R mutations in 15 out of 62 primary cases of PMBCL (24.2%) and in all PMBCL-derived cell lines tested. The majority of mutations (11/21; 52%) were hotspot single nucleotide variants in exon 8 leading to an I242N amino acid change in the transmembrane domain. Functional analyses establish this mutation as gain-of-function leading to constitutive activation of the JAK-STAT pathway and upregulation of downstream cytokine expression profiles and B cell specific antigens. Moreover, expression of I242N mutant IL4R in a mouse xenotransplantation model conferred growth advantage in vivo. The pattern of concurrent mutations within the JAK-STAT signaling pathway suggests additive/synergistic effects of these gene mutations contributing to lymphomagenesis. Our data establish IL4R mutations as novel driver alterations and provide a strong preclinical rationale for therapeutic targeting of JAK-STAT signaling in PMBCL.

  Study EGAS00001002796

• PATL2 is a key actor of oocyte maturation whose invalidation causes infertility in women and mice

  The genetic causes of oocyte meiotic deficiency (OMD), a form of primary infertility characterised by the production of immature oocytes, remain largely unexplored. Using whole exome sequencing, we found that 26% of a cohort of 23 subjects with OMD harboured the same homozygous nonsense pathogenic mutation in PATL2, a gene encoding a putative RNA-binding protein. Using Patl2 knockout mice, we confirmed that PATL2 deficiency disturbs oocyte maturation, since oocytes and zygotes exhibit morphological and developmental defects respectively. PATL2's amphibian orthologue is involved in the regulation of oocyte mRNA as a partner of CPEB. However, Patl2's expression profile throughout oocyte development in mice, alongside colocalisation experiments with Cpeb1, Msy2 and Ddx6 (three oocyte RNA-regulators) suggest an original role for Patl2 in Mammals. Accordingly, transcriptomic analysis of oocytes from WT and Patl2-/- animals demonstrated that in the absence of Patl2, expression levels of a select number of highly relevant genes involved in oocyte maturation and early embryonic development are deregulated. In conclusion, PATL2 is a novel actor of mammalian oocyte maturation whose invalidation causes OMD in humans.

  Study EGAS00001002903

• Ewings Sarcoma RNA-Seq

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003062

• Low and variable tumor reactivity of the intratumoral TCR repertoire in human cancers

  Infiltration of human cancers by T cells is generally interpreted as a sign of immune recognition, and there is a growing effort to reactivate dysfunctional T cells at such tumor sites. However, these efforts only have value if the intratumoral T cell receptor (TCR) repertoire of such cells is intrinsically tumor-reactive, and this has not been established in an unbiased manner for most human cancers. To address this issue, we analyzed the intrinsic tumor-reactivity of the intratumoral TCR repertoire of CD8+ T cells in ovarian and colorectal cancer – two tumor types for which T cell infiltrates form a positive prognostic marker. Data obtained demonstrate that a capacity to recognize autologous tumor is limited to approximately 10% of intratumoral CD8+ T cells. Furthermore, in two out of four patient samples tested, no tumor-reactive TCRs were identified, despite infiltration of their tumors by T cells. These data indicate that the intrinsic capacity of intratumoral T cells to recognize adjacent tumor tissue can be rare and variable, and suggest that clinical efforts to reactivate intratumoral T cells will benefit from approaches that simultaneously increase the quality of the intratumoral TCR repertoire.

  Study EGAS00001003119

• Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis -Part II

  Papillary thyroid carcinoma (PTC) has a wide geographic variation in incidence, being highest in Saudi Arabia where it is only second to breast cancer as the most common cancer among females. Genomic profiling of PTC from Saudi Arabia has not been attempted previously. We performed whole-exome sequencing of 101 PTC samples and the corresponding genomic DNA to identify genes with recurrent somatic mutations, followed by sequencing of these genes using a next-generation gene panel approach in further 785 samples. In addition to BRAF, N-RAS and H-RAS, which have previously been shown to be recurrently mutated in PTC, our analysis highlights additional genes, including thyroglobulin (TG), which harbored somatic mutations in 3% of the entire cohort. Surprisingly, although TG mutations were not exclusive to mutations in the RAS-MAP kinase pathway, their presence was associated with a significantly worse clinical outcome, which suggests a pathogenic role beyond driving initial oncogenesis. Analysis of metastatic PTC tissue revealed significant enrichment for TG mutations (p<0.001), including events of apparent clonal expansion. Our results suggest a previously unknown role of TG somatic mutations in the pathogenesis of PTC and its malignant evolution.

  Study EGAS00001003357

• Clinical significance of novel subtypes of acute lymphoblastic leukemia in the context of minimal residual disease-directed therapy

  We evaluate clinical significance of recently identified subtypes of acute lymphoblastic leukemia (ALL) in 598 children treated with minimal residual disease (MRD)-directed therapy. Among the 16 B-ALL and 8 T-ALL subtypes identified by next generation sequencing, ETV6-RUNX1, high-hyperdiploid and DUX4-rearranged B-ALL had the best five-year event-free survival rates (95% to 98.4%); TCF3-PBX1, PAX5alt, T-cell, ETP, iAMP21, and hypodiploid ALL intermediate rates (80.0% to 88.2%); and BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like and KMT2A-rearranged ALL the worst rates (64.1% to 76.2%). All but three of the 142 patients with day-8 blood MRD <0.01% remained in remission. Among new subtypes, intensified therapy based on day-15 MRD≥1% improved outcome of DUX4-rearranged, BCR-ABL1-like, and ZNF384-rearranged ALL, and achievement of day-42 MRD<0.01% did not preclude relapse of PAX5alt, MEF2D-rearranged and ETV6-RUNX1-like ALL. Thus, new subtypes including DUX4-rearranged, PAX5alt, BCR-ABL1-like, ETV6-RUNX1-like, MEF2D-rearranged and ZNF384-rearranged ALL have important prognostic and therapeutic implications.

  Study EGAS00001005084

• Rucaparib in patients presenting a metastatic breast cancer with Homologous Recombination Deficiency, without germline BRCA1/2 mutation, pronostic value of high LOH genomic score and predictive value of HRDetect (microarray)

  Breast cancer may present genomic alterations leading to homologous recombination deficiency. PARP inhibitors have proved their efficacy in patients with HER2-negative metastatic breast cancer (mBC) harboring germline (g) BRCA1/2 mutations. We conducted the phase 2 RUBY trial to assess the efficacy of rucaparib in HER2-negative mBC with high genomic loss of heterozygosity (LOH) score or somatic, without gBRCA1/2 mutation. 220 of 711 patients with mBC screened for LOH presented high LOH score which was associated with a higher likelihood of death (HR = 1.39, 95% CI: 1.11-1.75, p = 0.005). The primary objective was not reached with a clinical benefit rate (objective response or SD>16 weeks) of 13.5%. Two LOH-high patients, without somatic BRCA1/2 mutation, presented a complete and durable response (14 and 32 months). HRDetect tended to be associated with response to rucaparib, whithout reaching statistical significance (median HRDetect responders versus non responders: 0.465 versus, 0.040, p = 0.2135). Our data suggests that a small subset of patients with high LOH score could derive benefit from PARP inhibitors.

  Study EGAS00001004755

• UAMS Smoldering Myeloma Myeloma Sequencing

  Identifying high risk smoldering myeloma patients and progression mechanism is a prerequisite to implement effective inception strategies and curve myeloma related morbidity and mortality. We hypothesize that genomics may help identify determinants of progression that may help predict outcome and offer effective chemo preventive targets. Eighty-two patients underwent a custom targeted panel sequencing with an additional capture for the translocation loci. These results were compared to 223 newly diagnosed patients (EGAS00001003223 and EGAD00001004117) and 17 MGUS and 10 early myeloma patients. DNA was obtained from either CD138+ cells from the bone marrow of multiple myeloma patients (tumor) or from stem cell harvests or peripheral blood cells from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using the HyperPlus kit (KAPA Biosystems). After PCR amplification the libraries were hybridized with probes against either a targeted panel consisting of 140 genes and chromosomal regions (Nimblegen) using SeqCap reagents (Nimblegen). Hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads. Overall, these data highlight the importance of dysregulation of the MAPK pathway in the progression to MM.

  Study EGAS00001003629

• Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis

  Papillary thyroid carcinoma (PTC) has a wide geographic variation in incidence, being highest in Saudi Arabia where it is only second to breast cancer as the most common cancer among females. Genomic profiling of PTC from Saudi Arabia has not been attempted previously. We performed whole-exome sequencing of 101 PTC samples and the corresponding genomic DNA to identify genes with recurrent somatic mutations, followed by sequencing of these genes using a next-generation gene panel approach in further 785 samples. In addition to BRAF, N-RAS and H-RAS, which have previously been shown to be recurrently mutated in PTC, our analysis highlights additional genes, including thyroglobulin (TG), which harbored somatic mutations in 3% of the entire cohort. Surprisingly, although TG mutations were not exclusive to mutations in the RAS-MAP kinase pathway, their presence was associated with a significantly worse clinical outcome, which suggests a pathogenic role beyond driving initial oncogenesis. Analysis of metastatic PTC tissue revealed significant enrichment for TG mutations (p<0.001), including events of apparent clonal expansion. Our results suggest a previously unknown role of TG somatic mutations in the pathogenesis of PTC and its malignant evolution.

  Study EGAS00001001788

• The proliferative history shapes the DNA methylome of B-cell tumor and predicts clinical outcome

  Here, we provide access to CLL and DLBCL DNA methylation and gene expression data in the context of a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. This approach showed that differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

  Study EGAS00001004640

• Identification of rare germline variants in familial multiple myeloma

  The risk of developing multiple myeloma (MM) in first-degree relatives of affected individuals is 2 to 4 fold higher than in the general population, implying the existence of inherited susceptibility. While over 20 common, low-penetrance variants have been associated with MM risk, less is known about rare, high-penetrance germline predisposition. We aimed to identify and characterize high/moderated penetrance germline variants in MM families using whole genome sequencing. For that purpose we identified families with at least two individuals diagnosed with MM or its precursors MGUS and smoldering MM. Also, patients with solitary plasmacytoma and AL amyloidosis were enrolled. We also recruited unaffected family members, and they were screened to exclude undetected MM or its precursor stages. After whole genome sequencing of blood DNA and filtering for rare variants (minor allele frequency <0.1%) and family segregation, variants are prioritized by a range of published in silico tools and literature search. The most promising variants are subjected to experimental validation to evaluate their functional relevance. Our data will provide guidance to worldwide efforts to unravel germline predisposition in MM and validation of the proposed variants and genes.

  Study EGAS00001004734

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

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

  Study EGAS00001004791

• Risk and modifying factors in Frontotemporal Dementia

  Understanding the molecular mechanisms underlying frontotemporal dementia (FTD) is essential for the development of successful therapies. Here, we present Phase 1 of a multi-omics, multi-model data resource for FTD research which will allows in-depth molecular research into these mechanisms. We have integrated and analysed data from the frontal lobe of FTD patients with mutations in MAPT, GRN and C9orf72 and detected common and distinct dysregulated cellular pathways. Our results highlight that excitatory neurons are the most vulnerable neuronal cell type and that vascular aberrations are a common hallmark in FTD. Via integration of multi-omics data, we detected several transcription factors and pathways which regulate the strong neuroinflammation observed in FTD-GRN. Finally, using small RNA-seq data and verification experiments in cellular models, we identified several up-regulated miRNAs that inhibit cellular trafficking pathways in FTD and lead to microglial activation. In this work we shed light on novel mechanistic and pathophysiological hallmarks of FTD. In addition, we believe that this comprehensive, multi-omics data resource will further mechanistic FTD research by the community.

  Study EGAS00001004895

• Single cell RNA sequencing of synovial B cells in early Rheumatoid Arthritis

  Rheumatoid Arthritis (RA) is a prevalent autoimmune disease characterized by inflammation of the peripheral joints and occurence of autoantibodies, i.e. Rheumatoid Factor (RF) and anti-citrullinated protein antibodies (ACPA), in circulation. Previous studies have shown that there are memory B cells and autoantibody producing plasma cells present in the joint tissue of patients with long-standing RA. However, it has remained unclear, whether these B cell populations are already present in the joint tissue of patients at the disease onset and whether these are autoreactive. Here, we used a single cell RNA sequencing approach to dissect the B cell repertoire at this early timepoint. We found evidence for B and T cell interaction and presence of memory and plasma cell pools in ACPA- and ACPA+ RA. Our results demonstrated clonal relationships between the memory and plasma cell compartments and autoreactivity within the plasma cell pool. These findings challenge our understanding of the local adaptive immune response in the joints of ACPA- and ACPA+ RA patients with direct implications for B and T cell targeting therapy in both patient subgroups.

  Study EGAS00001005144

• Epigenome wide DNA methylation assay of gingivo-buccal oral squamous cell carcinoma using single base resolution high throughput array

  Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. To achieve this, we have carried out epigenomic (n=87) profiling of paired tumor-normal tissues collected from OSCC-GB patients from India and genome-wide DNA methylation assays were performed. DNA was isolated using DNeasy Blood and Tissue Kit (Qiagen). The purity and concentration was estimated using NanoDrop 2000 (Thermo Fisher Scientific). Approximately 500 ng genomic DNA from each sample was used for sodium bisulfite conversion using the EZ DNA methylation Gold Kit (Zymo Research, USA). Genome-wide DNA methylation was assayed using iScan (Illumina), for paired tumour and adjacent normal samples of 25 patients using the Infinium Human Methylation450 BeadChip and of 62 patients using the Infinium MethylationEPIC BeadChip.

  Study EGAS00001003896

• Whole genome sequence and RNA-seq data from paired tumour and germline samples from mesothelioma patients.

  Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with RNA sequencing (RNA-seq) analysis of the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials. We analysed somatic mutation and RNA-seq data from 229 MPM samples, including 58 MPM samples that had undergone WGS from our own institutions, together with other published data. This combined analysis identified somatic driver genes, including newly identified candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed dominant signatures and showed that defects in homologous recombination repair were infrequent in our cohort. Within the tumour immune environment we identified high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immunosuppressive environment. A small subset of samples had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a ‘hot’ immune environment which is independent of the somatic mutations. We propose that our findings on genomic changes and subtypes of immune microenvironments may influence therapeutic planning in the future.

  Study EGAS00001005196

• The molecular basis of T-PLL is an actionable perturbation of TCL1/ATM- and epigenetically instructed damage responses

  T-cell prolymphocytic leukemia (T-PLL) is a rare and poor-prognostic mature T-cell malignancy. To address its incomplete molecular concept, we integrated large-scale profiling data of alterations in gene expression, allelic copy number (CN), and nucleotide sequences in 111 well-characterized patients. Besides prominent signatures of T-cell activation and prevalent clonal variants, we also identified novel hot-spots for CN variability, fusion molecules, alternative transcripts, and progression-associated dynamics. The overall lesional spectrum of T-PLL is mainly annotated to axes of DNA damage responses, T-cell receptor/cytokine signaling, and histone modulation. We formulate a multi-dimensional model of T-PLL pathogenesis centered around a unique combination of TCL1 overexpression with damaging ATM aberrations as initiating core lesions. The effects imposed by TCL1 cooperate with compromised ATM towards a leukemogenic phenotype of impaired DNA damage processing. Dysfunctional ATM appears inefficient in alleviating elevated redox burdens and telomere attrition and in evoking a p53-dependent apoptotic response to genotoxic insults. As non-genotoxic strategies, synergistic combinations of p53 reactivators and deacetylase inhibitors reinstate such cell death execution.

  Study EGAS00001002744

• MutWP4__CRUK_Grand_Challenge_Mutographs_of_Cancer__Gastric_Organoids

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

  Study EGAS00001003151

• Whole-Genome Sequencing of Salivary Gland Adenoid Cystic Carcinoma

  Adenoid cystic carcinomas (ACCs) of the salivary glands are challenging to understand, treat, and cure. To better understand the genetic alterations underlying the pathogenesis of these tumors, we performed comprehensive genome analyses of 25 fresh-frozen tumors, including whole genome sequencing, expression and pathway analyses. In addition to the well-described MYB-NFIB fusion which was found in 11 tumors (44%), we observed five different rearrangements involving the NFIB transcription factor gene in seven tumors (28%). Taken together, NFIB translocations occurred in 15 of 25 samples (60%, 95%CI=41–77%). In addition, mRNA expression analysis of 17 tumors revealed overexpression of NFIB in ACC tumors compared with normal tissues (p=0.002). There was no difference in NFIB mRNA expression in tumors with NFIB fusions compared to those without. We also report somatic mutations of genes involved in the axonal guidance and Rho family signaling pathways. Finally, we confirm previously described alterations in genes related to chromatin regulation and Notch signaling. Our findings suggest a separate role for NFIB in ACC oncogenesis and highlight important signaling pathways for future functional characterization and potential therapeutic targeting.

  Study EGAS00001002812

• Cell-free DNA analysis reveals POLR1D-mediated resistance to bevacizumab in colorectal cancer

  Bevacizumab is an approved anti-angiogenic drug for patients with metastasized colorectal cancer (mCRC) targeting VEGF. The survival benefit of anti-VEGF therapy in mCRC patients is limited to a few months and acquired resistance mechanisms are greatly unknown. Using plasma DNA, we studied the evolution of tumor genomes in a cohort of patients with mCRC (n=150) and observed a recurrent focal amplification (8.7% of cases) on chromosome 13q12.2. Analysis of TCGA data (n=619) suggested an association with later stages, which we confirmed by longitudinal plasma analyses. We defined the minimally amplified region and studied the mechanistic consequences of copy number gain of the involved genes. The amplification of one gene, POLR1D, impacted cell proliferation, resulting in upregulation of VEGFA, an important regulator of angiogenesis which has been implicated in the resistance to bevacizumab. In several patients, we observed the emergence of this 13q12.2 amplicon under bevacizumab treatment, which was invariably associated with evolution of therapy resistance. Hence, we describe a novel resistance mechanism against a widely applied treatment in mCRC patients which will impact clinical management

  Study EGAS00001003791

• Circulating tumor cells Exome sequencing from breast cancer

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

  Study EGAS00001005228

• SS18-SSX-mediated hijacking of BAF complexes drives synovial sarcoma

  Synovial sarcoma (SS) is defined by a recurrent t(x;18) chromosomal translocation, which produces the hallmark SS18-SSX oncogenic fusion. Incorporation of SS18-SSX into BAF complexes renders BAF complexes aberrant in two distinct manners: the addition of 78aa of SSX onto SS18, and concomitant loss of BAF47 assembly. However, the importance and functional contributions of each of these perturbations on BAF complex targeting and gene expression regulation remain unclear. Here we use an integrative set of genomic approaches in human cancer cell lines and primary tumor samples to define the mechanistic consequences of the SS18-SSX fusion oncoprotein. We find that SS18-SSX hijacks BAF complexes to broad polycomb domains to activate bivalent genes, driving a unique gene expression program distinct from other loss-of-function BAF complex malignancies. Importantly, restoration of BAF47 rescues enhancer activation but is dispensable for proliferative arrest in cell lines. These results demonstrate that gain-of-function SS18-SSX-mediated BAF complex targeting and gene activation is the driving event in SS, and present a mechanism by which distinct functions of BAF complexes can be co-opted to drive oncogenesis.

  Study EGAS00001002920

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

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

  Study EGAS00001005366

• Patient-derived neuroblastoma model system OHC-NB1

  Accurate disease modeling presents a bottleneck against effective translation of novel neuroblastoma therapies into the clinic. Recently, considerable intra-tumor genetic heterogeneity has been observed in neuroblastoma, which calls for preclinical models that reflect this feature. Flexibility in a testing platform is also desirable to support assessment of different endpoints and tumor cell characteristics. We present the novel neuroblastoma model OHC-NB1, propagated from a bone marrow metastasis from a patient with INRG stage M, MYCN-amplified neuroblastoma at first diagnosis. To model different aspects of the disease, several preclinical models were developed from the bone marrow aspirate: a tissue plastic-adherent cell line, a 3-dimensional (3D) spheroid culture, and cell line-derived xenografts that were subcutaneously passaged in mice. Typical aspects of MYCN-amplified neuroblastoma, such as MYCN amplification in double minutes, 1p deletion, 17q gain, and a diploid karyotype, are present in OHC-NB1 and persist in all culture types. Whole exome sequencing revealed intra-tumor heterogeneity in the bone marrow metastasis reflected in the culture models. The set of OHC-NB1 cultures with a genetically heterogeneous MYCN-amplified neuroblastoma background provide a flexible preclinical testing system that recapitulates important features of metastatic neuroblastoma.

  Study EGAS00001003031

• Germline biallelic mutation affecting the transcription factor Helios causes pleiotropic defects of immunity

  We studied a patient with recurrent respiratory infections and hypogammaglobulinemia and identified a germline homozygous missense mutation in IKZF2 encoding Helios (p.Ile325Val). We show that HeliosI325V retains DNA-binding and dimerization properties, but loses interaction with several partners, including epigenetic remodelers HDAC1, HDAC3 and the ATAC complex. Single-cell RNA-sequencing of peripheral blood mononuclear cells revealed gene expression signatures indicative of a shift towards pro-inflammatory, effector-like status in the patient’s T cells. We observed an upregulation of the anti-inflammatory gene, TSC22D3, encoding the glucocorticoid-induced leucine zipper (GILZ), whose expression was reported to increase upon IL-2 deprivation. We validated this finding in primary T cells, where we observed a pronounced reduction in IL-2 production upon their stimulation, together with a defect in their ability to proliferate, rendering them more anergic. Collectively, we identify a novel germline-encoded inborn error of immunity and define a role for Helios in conventional T cells, whereby interactions with specific binding partners is necessary to mediate the transcriptional programs that enable T-cell homeostasis in health and disease.

  Study EGAS00001005675

• A Single-cell Atlas of Progressive TKI Resistance in Chronic Myeloid Leukemia

  High-dimensional scRNA seq and CyTOF were used to identify features predictive of treatment outcome at diagnosis. An optimal TKI response was characterised by erythroid lineage skewing of the CD34+ HSPC compartment. While TKI-resistant CP was reminiscent of a BC-like state, where HSPCs were enriched for inflammation, stemness and quiescence. We designed a machine-learning approach which assessed the prognostic potential of all 32 sub-populations of our scRNA seq dataset. LSCs and NK populations harboured the highest prognostic power at diagnosis. Within the LSCs, erythroid lineage priming and MYC activation were features of an optimal and poor TKI responses, respectively. With the NK cell compartment, apart from a higher abundance of NK cells in TKI optimal responders, a memory-like HLA-DR+ adaptive NK and KLRC1+ NK populations were hallmarks of an optimal and poor TKI responses, respectively. Mechanistically, both cell intrinsic and extrinsic mechanisms contributed towards the higher sensitivity of EPs to TKI therapy. In summary, our high-dimensional single-cell atlas of TKI resistance highlights pivotal transcriptional features associated with TKI-resistance disease, some of which have the potential to be exploited as biomarkers.

  Study EGAS00001005509

• Clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancers defined by genomic analysis

  Ductal Carcinoma In Situ (DCIS) is the most common form of pre-invasive breast cancer and despite treatment a small fraction (5-10%) of DCIS patients present with invasive disease many years later. A fundamental biologic question is whether the invasive disease recurring in the same breast is established by tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial pure DCIS lesion and paired invasive recurrent tumors in 95 patients together with single cell DNA sequencing in a subset of cases. Our data shows that in 75% the invasive recurrence was clonally related to the initial DCIS, suggesting that the tumor cells were not eliminated during the initial treatment with surgery +/- radiotherapy. Surprisingly however, 18% were clonally unrelated to the DCIS, representing new independent lineages, and 7% of cases were ambiguous. Our findings show that although DCIS is often the precursor of invasive recurrence, a significant fraction of invasive recurrences are unrelated to the initial DCIS. This knowledge is essential for accurate risk evaluation of DCIS treatment de-escalation strategies and the identification of predictive biomarkers.

  Study EGAS00001005784

• De novo assembly of 150 Danish genomes reveals rich structural complexity

  Most known genetic variation in human genomes has been called from comparison of short reads to the reference genome, an approach biased against finding complex variation. We sequenced 150 individuals from 50 parent-offspring trios with multiple insert-size libraries to very high coverage. We show that each genome could be independently de novo assembled into a small number of high-quality scaffolds (median N50 > 21 Mb), each of quality comparable to long read assemblies while being very cost-effective. We show that our variant call set from comparing de novo assemblies is far more complete in terms of complex variation than previous studies. Importantly, even the complex 4-5 Mb extended MHC region was assembled and resolved into haplotypes, revealing >700kb novel sequence in this important region of the genome, and major parts of the Y chromosome including some palindromes were assembled with high accuracy. Finally, we show that our variant call-set allows for the genotyping of many more complex variants when used as a reference-panel for imputation into SNP-chip data or into previously resequenced genomes.

  Study EGAS00001002108

• Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing

  Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

  Study EGAS00001003190

• B cell receptor repertoire kinetics after SARS-CoV-2 infection and vaccination

  B cells play a central role in the immune response to both SARS-CoV-2 infection and vaccination, but the development of the B cell receptor (BCR) repertoire in both contexts has not been defined nor compared. We analysed serial samples from 171 SARS-CoV-2-infected individuals with a range of disease severities together with 63 vaccine recipients, and found marked differences in the global BCR repertoire after natural infection compared to vaccination. Following infection, the proportion of BCRs bearing IgG1/3 and IgA1 isotypes increased, somatic hypermutation (SHM) was markedly decreased and, in patients with severe disease, expansion of IgM and IgA clones was observed. In contrast, after vaccination the proportion of BCRs bearing IgD/M isotypes increased, SHM was unchanged and expansion of IgG clones was prominent. Infection generated a broad distribution of SARS-CoV-2-specific clones predicted to target the spike protein whilst vaccination produced a more focused response mainly targeting the spike’s receptor-binding domain. These findings offer insights into how different immune exposure to SARS-CoV-2 impacts upon BCR repertoire development, potentially informing vaccine strategies.

  Study EGAS00001005816

• Target gene sequencing for human normal endometrial glands

  To elucidate the timing and mechanism of the clonal expansion of somatic mutations in cancer-associated genes in the normal endometrium, we conducted target sequencing of 112 genes for 1,298 endometrial glands and matched blood samples from 36 women. By collecting endometrial glands from different parts of the endometrium, we showed that multiple glands with the same somatic mutations occupied substantial areas of the endometrium. The 112 genes are as follows: ABCC1, ACRC, ANK3, ARHGAP35, ARID1A, ARID5B, ATCAY, ATM, ATR, BARD1, BCOR, BRCA1, BRCA2, BRD4, BRIP1, CAMTA1, CDC23, CDYL, CFAP54, CHD4, CHEK1, CHEK2, CTCF, CTNNB1, CUX1, DGKA, DISP2, DYNC2H1, EMSY, FAAP24, FAM135B, FAM175A, FAM65C, FANCA, FANCB, FANCC, FANCD2, FANCE, FANCF, FANCG, FANCI, FANCL, FANCM, FAT1, FAT3, FBN2, FBXW7, FGFR2, FRG1, GPR50, HEATR1, HIST1H4B, HNRNPCL1, HOOK3, KIAA1109, KIF26A, KMT2B, KMT2C, KRAS, LAMA2, LRP1B, MLH1, MON2, MRE11A, MSH2, MSH6, MTOR, NBN, PALB2, PHEX, PIK3CA, PIK3R1, PLXNB2, PLXND1, PMS2, POLE, POLR3B, PPP2R1A, PTEN, PTPN13, RAD50, RAD51, RAD51B, RAD51C, RAD51D, RAD52, RAD54B, RAD54L, RICTOR, SACS, SIGLEC9, SLC19A1, SLX4, SPEG, STT3A, TAF1, TAF2, TAS2R31, TFAP2C, TNC, TONSL, TP53, TTC6, UBA7, VNN1, WT1, XIRP2, ZBED6, ZC3H13, ZFHX3, ZFHX4, ZMYM4.

  Study EGAS00001005914

• Identification of Germline Monoallelic Mutations in IKZF2 in Patients with Immune Dysregulation.

  Helios, encoded by IKZF2, is a member of the Ikaros family of transcription factors with pivotal roles in T-follicular helper, NK- and T-regulatory cell physiology. Somatic IKZF2 mutations are frequently found in lymphoid malignancies. Although germline mutations in IKZF1 and IKZF3, encoding Ikaros and Aiolos, have recently been identified in patients with phenotypically similar immunodeficiency syndromes, the effect of germline mutations in IKZF2 on human hematopoiesis and immunity remains enigmatic. We identified germline IKZF2 mutations (one nonsense (p.R291X)- and 4 distinct missense variants) in six patients with systemic lupus erythematosus, immune thrombocytopenia or EBV-associated hemophagocytic lymphohistiocytosis. Patients exhibited hypogammaglobulinemia, decreased number of T-follicular helper and NKcells. Single-cell RNA sequencing of PBMCs from the patient carrying the R291X variant revealed upregulation of pro-inflammatory genes associated with T-cell receptor activation and T-cell exhaustion. Functional assays revealed the inability of HeliosR291X to homodimerize and bind target DNA as dimers. Moreover, proteomic analysis by proximity-dependent Biotin Identification revealed aberrant interaction of 3/5 Helios mutants with core components of the NuRD complex conveying HELIOS mediated epigenetic and transcriptional dysregulation.

  Study EGAS00001005874

• Changes in alternative splicing and associated neo-antigens due to therapy

  Recent studies have identified alternative splicing (AS) as a novel source of neoantigens for immunotherapy. Surprisingly little is known about the AS milieu in recurrent glioblastoma (GBM), despite this being the venue for most clinical trials. We profiled 29 primary-recurrent paired human GBM specimens via RNA sequencing and re-analyzed RNA-sequencing data from non-malignant human brain tissues. From these data, we reconstructed the landscape of AS in GBM through recurrence and contrasted that to isoforms observed in non-malignant brain. The AS events we identified were cross-referenced with single-cell GBM atlases to determine cell-type specific splicing patterns. From this we identified novel splicing events in cell-surface proteins that are suitable targets for engineered T-cell therapies. We found recurrent-specific isoforms of mitogen-activated kinase pathway genes which are expressed exclusively by GBM stem-like cells that enhance invasiveness. These studies shed light on the effect of therapy on AS and identify novel targets for emerging immunotherapies.

  Study EGAS00001004524

• Patient-tailored design for selective co-inhibition of leukemic cell subpopulations

  The extensive primary and secondary drug resistance in many cancer types requires rational approaches to design personalized and selective combinatorial therapies that do not only show synergistic effect in overall cancer cell killing but also result in minimal toxic side effects on non-malignant cells. To address the combinatorial explosion in the number of relevant combinations, we implemented a machine learning approach that prioritizes patient-customized drug combinations with a desired synergy-efficacy-toxicity balance by combining single-cell RNA-sequencing with ex vivo single-agent testing in scarce patient-derived primary cells. When applied to two diagnostic and two refractory AML patient cases, each with a different genetic background, our integrated approach accurately predicted patient-specific combinations that were shown to result not only in synergistic cancer cell co-inhibition but were also capable of targeting specific AML cell subpopulations that emerge in differing stages of disease pathogenesis or treatment regimens. Our data-driven approach provides an unbiased means for systematic identification of personalized combinatorial regimens that selectively co-inhibit leukemic cells while avoiding inhibition of non-malignant cells, and highlight the relevance of considering cell heterogeneity for personalized cancer therapy.

  Study EGAS00001004614

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

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

  Study EGAS00001004684

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

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

  Study EGAS00001004786

• Accurate detection of tumor-specific fusion genes reveals strongly immunogenic personal neo-antigens

  Fusion genes arising from cancer-associated somatic mutations are a potential rich source for highly immunogenic neo-antigens. However, their exploitation as targets for personalized cancer immunotherapy is currently limited by the lack of computational tools allowing transcriptome-wide identification of unique fusion genes in an accurate and sensitive manner. Here, we present EasyFuse, a computational pipeline, to detect individual and cancer-specific fusion genes in next-generation-sequencing transcriptome data obtained from human cancer samples. Using machine learning, EasyFuse predicts personal fusion genes with high precision and sensitivity and outperforms previously described approaches as qualified by an unprecedented ground-truth dataset of >1500 verification experiments in relevant patient samples. By testing immunogenicity with autologous blood lymphocytes from cancer patients we detected pre-established CD4+ and CD8+ T cell responses for 10 of 21 (48%), and for 1 of 30 (3%) of identified fusion genes, respectively. In conclusion, we demonstrate accurate detection of cancer-specific fusion genes. The high frequency of T cell responses detected in cancer patients support the relevance of private fusion genes as neo-antigens for personalized immunotherapies, especially for tumors with low point mutation burdens.

  Study EGAS00001004877

• Heat selection enables highly scalable methylome profiling in cell-free DNA for noninvasive monitoring of cancer patients

  Genome-wide analysis of cell-free DNA (cfDNA) methylation profile has been recognized as a promising approach for sensitive and specific detection of many cancers. However, scaling such genome-wide assays for clinical translation is impractical due to the high cost of whole genome bisulfite sequencing. We have shown that the small fraction of GC-rich genome is highly enriched in CpG sites and disproportionately harbors the majority of cancer-specific methylation signature. Here, we report on the simple but effective Heat enrichment of CpG-rich regions for Bisulfite Sequencing (Heatrich-BS) platform that allows for focused methylation profiling in these highly informative regions. Our novel method and bioinformatics algorithm enable accurate tumor burden estimation with high sensitivity and quantitative tracking of colorectal cancer patient’s response to treatment, at much reduced sequencing cost suitable for frequent monitoring. We also show, for the first time, tumor epigenetic subtyping from cfDNA using Heatrich-BS, which could enable patient stratification from non-invasive liquid biopsy. As such, Heatrich-BS holds great potential for highly scalable screening and regular monitoring of cancer using liquid biopsy.

  Study EGAS00001006198

• Clonal somatic copy number altered driver events inform drug sensitivity in high-grade serous ovarian cancer

  Chromosomal instability is a major challenge to patient stratification and targeted drug development for high-grade serous ovarian carcinoma (HGSOC). Here we show that somatic copy number alterations (SCNAs) in frequently amplified HGSOC cancer genes significantly correlate with gene expression and methylation status. We identified five prevalent clonal driver SCNAs (chromosomal amplifications encompassing MYC, PIK3CA, CCNE1, KRAS and TERT) from multi-regional HGSOC data and reasoned that their strong selection should prioritise them as key biomarkers for targeted therapies. We used primary HGSOC spheroid models to test interactions between in vitro targeted therapy and SCNAs. MYC chromosomal copy number was associated with in-vitro and clinical response to paclitaxel and in-vitro response to mTORC1/2 inhibition. Activation of the mTOR survival pathway in the context of MYC-amplified HGSOC was statistically associated with increased prevalence of SCNAs in genes from the PI3K pathway. Co-occurrence of amplifications in MYC and genes from the PI3K pathway was independently observed in squamous lung cancer and triple negative breast cancer. These results suggest that identifying co-occurrence of clonal driver SCNA genes could be used to tailor therapeutics for precision medicine.

  Study EGAS00001006200

• MutWP4__CRUK_Grand_Challenge_Mutographs_of_Cancer__Pancreatic_Organoids

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

  Study EGAS00001003520

• MGP Panel: a comprehensive targeted genomics panel for molecular profiling of multiple myeloma patients

  We designed a comprehensive multiple myeloma (MM) targeted sequencing panel to identify common genomic abnormalities in a single assay and validated it against known standards. The panel comprised 228 genes/exons for mutations, 6 regions for translocations, and 56 regions for copy number abnormalities (CNAs). Toward panel validation, targeted sequencing was conducted on 233 patient samples and further validated using clinical fluorescence in situ hybridization (FISH) (translocations), multiplex ligation probe analysis (MLPA) (CNAs), whole genome sequencing (WGS) (CNAs, mutations, translocations) or droplet digital PCR (ddPCR) of known standards (mutations). Canonical IgH translocations were detected in 43.2% of patients by sequencing, and aligned with FISH except for one patient. CNAs determined by sequencing and MLPA for 22 regions were comparable in 103 samples and concordance between platforms was R2=0.969. VAFs for 74 mutations were compared between sequencing and ddPCR with concordance of R2=0.9849. In summary, we have developed a targeted sequencing panel that is as robust or superior to FISH and WGS. This molecular panel is cost effective, comprehensive, clinically actionable and can be routinely deployed to assist risk stratification at diagnosis or post-treatment to guide sequencing of therapies.

  Study EGAS00001006222

• Discovery of a highly widespread bacteriophage family and its associations to metabolic syndrome gut microbiomes

  There is significant interest in altering the course of cardiometabolic disease development via the gut microbiome. Nevertheless, the highly abundant phage members of the complex gut ecosystem -which impact gut bacteria- remain understudied. Here, we characterized gut virome changes associated with metabolic syndrome (MetS), a highly prevalent clinical condition preceding cardiometabolic disease, in 196 participants with a combination of whole genome shotgun and virus like particle sequencing. MetS gut virome populations exhibited decreased richness and diversity, but larger inter-individual variation. These populations were enriched in phages infecting Bacteroidaceae and depleted in those infecting Ruminococcaeae. Differential abundance analysis identified eighteen viral clusters (VCs) as significantly associated with either MetS or healthy viromes. Among these are a MetS-associated Roseburia VC that is related to healthy control-associated Faecalibacterium and Oscillibacter VCs. Further analysis of these VCs revealed the Candidatus Heliusviridae, a highly widespread gut phage lineage found in 90+% of the participants. The identification of the temperate Ca. Heliusviridae provides a novel starting point to a better understanding of the effect that phages have on their bacterial hosts and the role that this plays in MetS

  Study EGAS00001006260

• Clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancers defined by genomic analysis

  Ductal Carcinoma In Situ (DCIS) is the most common form of pre-invasive breast cancer and despite treatment a small fraction (5-10%) of DCIS patients present with invasive disease many years later. A fundamental biologic question is whether the invasive disease recurring in the same breast is established by tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial pure DCIS lesion and paired invasive recurrent tumors in 95 patients together with single cell DNA sequencing in a subset of cases. Our data shows that in 75% the invasive recurrence was clonally related to the initial DCIS, suggesting that the tumor cells were not eliminated during the initial treatment with surgery +/- radiotherapy. Surprisingly however, 18% were clonally unrelated to the DCIS, representing new independent lineages, and 7% of cases were ambiguous. Our findings show that although DCIS is often the precursor of invasive recurrence, a significant fraction of invasive recurrences are unrelated to the initial DCIS. This knowledge is essential for accurate risk evaluation of DCIS treatment de-escalation strategies and the identification of predictive biomarkers.

  Study EGAS00001006306

• Pharmacogenomic analysis of patient-derived tumor cells in gynecologic cancers

  Gynecologic malignancy is one of the leading causes of mortality in female adults worldwide. Comprehensive genomic analysis has revealed a list of molecular aberrations that are essential to tumorigenesis, progression, and metastasis of gynecologic tumors. However, targeting such alterations has frequently led to treatment failures due to underlying genomic complexity and activation of various tumor cell survival pathway molecules. A compilation of molecular characterization of tumors with pharmacological drug response is the next step towards clinical application of patient-tailored treatment regimens. Toward this goal, we have established a library of 139 gynecologic tumors including cervical, endometrial, and epithelial ovarian cancers (EOCs) and uterine sarcomas that were genomically and/or pharmacologically annotated and explored dynamic pharmacogenomic associations against 37 molecularly targeted drugs. We discovered lineage-specific drug sensitivities based on subcategorization of gynecologic tumors and identified TP53 mutation as a molecular determinant that elicit therapeutic response to poly (ADP-Ribose) polymerase (PARP) inhibitor. We further identified transcriptome expression of inhibitor of DNA biding 2 (ID2) as a potential predictive biomarker for treatment response to olaparib. Together, our results demonstrate potential utility of rapid drug screening combined with genomic profiling for precision treatment of gynecologic cancers.

  Study EGAS00001003556

• An Atlas of Cells in the Human Tonsil

  Palatine tonsils are secondary lymphoid organs representing the first line of immunological defense against inhaled or ingested pathogens. Here, we provide a comprehensive census of cell types forming the human tonsil by applying single cell transcriptome, epigenome, proteome and adaptive immune repertoire as well as spatial transcriptomics, resulting in an atlas of >350,000 cells. We provide a glossary of over one hundred annotated cell types and states, and disentangle gene regulatory mechanisms that drive cells through specialized lineage trajectories. Exemplarily, we stratify multiple tonsil-resident myeloid slancyte subtypes, establish a distant BCL6 superenhancer as locally active in both follicle-associated T and B cells, and describe SIX5 as a potentially novel transcriptional regulator of plasma cell maturation. Further, our atlas is a reference map to understand alterations observed in disease, here we discover immune-phenotype plasticity in tumoral cells and microenvironment shifts of mantle cell lymphomas (MCL). To facilitate such reference-based analysis, we develop HCATonsilData and SLOcatoR, a computational framework that provides programmatic and modular access to our dataset; and allows the straightforward annotation of future single-cell profiles from secondary lymphoid organs.

  Study EGAS00001006375

• Treatment-mediated selection of lethal prostate cancer clones defined by copy number architectures

  Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient relationships of metastases that evade treatment, we performed genome- wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 142 metastatic regions from 10 organs harvested post-mortem from nine men who died from prostate cancer. We identified diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number change, we confirmed a common clone of origin across metastases and diagnostic biopsies; and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. Autosome-defined clusters were characterized by cluster-specific AR gene architectures that in two index cases were topologically more congruent than by chance (p-values 0.03, 3.07x10-8). Integration with anatomical site suggested patterns of spread and points of genomic divergence. Copy number boundaries identified treatment-selected clones with putatively distinct lethal trajectories.

  Study EGAS00001006598

• DNA methylation landscape of prostate cancer

  Metastases from primary prostate cancers to rare locations, such as the brain, are becoming more common due to longer life expectancy thanks to improved treatments. Epigenetic dysregulation is a feature of primary prostate cancer, and distinct DNA methylation profiles are associated with the mutually exclusive SPOP mutant or TMPRSS2-ERG fusion genetic backgrounds. Using a cohort of 42 patients with prostate cancer brain metastases (PCBM), with matched primary samples available for 17 patients, we carried out DNA methylation analysis using Illumina EPIC array. We examined the epigenetic distinction between primary PCa and PCBM, the association between epigenetic alterations and mutational background, and particular epigenetic alterations that may be associated with PCBM in particular. Multiregion sampling of PCBM revealed epigenetic heterogeneity in metastases. We found aberrant methylation in PCBM is associated with mutational background, PRC2 complex activity and upregulation of EZH2, an effect that is particularly pronounced in SPOP mutant PCBM. While we detected a hypermethylator phenotype in PCBM, we also observed hypomethylation occurring at the promoters of genes involved in neuroactive ligand-receptor interactions or the calcium signalling pathway, such as GABRB3 and GRM5, a change not observed in non-brain prostate cancer metastases.

  Study EGAS00001006670

• Pre-neoplastic liver colonization by 11p15.5 altered mosaic cells in young children with hepatoblastoma

  In children, liver tumors are a heterogeneous group representing 1-2% of all cancers with the most common diagnosis being hepatoblastoma (HB). While HBs are predominantly sporadic, around 15% of cases develop as part of predisposition syndromes such as familial adenomatous polyposis (APC mutation) or Beckwith-Wiedemann (BWS; 11p15.5 locus altered). Here, we identified a mosaic genetic alteration of 11p15.5 locus in the liver in 17% of HB patients without a clinical diagnosis of BWS syndrome. These mosaic 11p15.5 alterations were never found in children with other types of liver tumor (hepatocellular carcinoma or adenoma, fibrolamellar carcinoma). Bulk transcriptomic analysis revealed the abundance of extracellular matrix, angiogenesis, and progenitor gene signatures in mosaic livers. In situ hybridization assay RNAscope allowed for the detection of 11p15.5 altered hepatocytes mosaicism in liver FFPE sections of HB patients presenting IGF2 overexpression and H19 downregulation together with an alteration of the zonal architecture of the liver. Spatial transcriptomics analysis revealed an overexpression of the IGF2 pathway, lipogenesis, and platelet activation in 11p15.5 altered hepatocytes. These data provide insights for 11p15.5 mosaicism detection as well as functional consequences of these alterations during the early steps of carcinogenesis.

  Study EGAS00001006692

• Targeted therapy of advanced parathyroid carcinoma guided by genomic and transcriptomic profiling (hipo_021)

  Parathyroid carcinoma (PC) is an orphan malignancy with a high risk of recurrence after surgery. Tumor-directed systemic treatments for PC are not established. We used whole-genome and RNA sequencing in four patients with advanced PC to identify molecular alterations that could guide clinical management. In two cases, the genomic and transcriptomic profiles provided targets for experimental therapies that resulted in biochemical response and prolonged disease stabilization: (i) immune checkpoint inhibition with pembrolizumab based on high tumor mutational burden and a single-base substitution signature associated with APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) overactivation; (ii) multi-receptor tyrosine kinase inhibition with lenvatinib due to overexpression of FGFR1 and RET and, later in the course of disease, PARP inhibition with olaparib prompted by signs of defective homologous recombination DNA repair. In addition, our data provided new insights into the molecular landscape of PC with respect to the genome-wide footprints of specific mutational processes and pathogenic germline alterations. These data underscore the potential of comprehensive molecular analyses to improve care for patients with ultra-rare cancers based on insight into disease biology.

  Study EGAS00001006747

• Case Report: Rare IKZF1 gene fusions identified in neonate with congenital KMT2A-rearranged Acute Lymphoblastic Leukemia

  Here we describe a rare case of congenital KMT2Ar ALL presenting with co-occurring IKZF1 gene fusions and a predictably aggressive disease trajectory. We report for the first time, the novel IKZF1::TUT1 and KDM2A::IKZF1 gene fusions. Rearrangements in-volving KMT2A are commonly retained in relapsed infant ALL, however, in this case the KMT2A::AFF1 gene fusion did not appear to be the lesion driving leukemic relapse. Instead, our data suggest that relapse was driven by IKZF1::TUT1. This gene fusion remained in all samples investigated, including the on-blinatumomab therapy sample taken immediately prior to relapse. Conversely, the KMT2A::AFF1 gene fusion was only detected in the diagnosis and refractory post-induction samples highlighting a key role for IKZF1::TUT1 in disease pathogenesis. Intriguingly, both IKZF1 gene fusions are predicted to be out-of-frame, however, our data demonstrate the IKZF1 gene is still expressed. This is not unprecedented and it has previously been observed that out-of-frame fusions can cause transcriptional activation/repression of genes involved in the fusions leading to increases or decreases of their expression and the as-sociated functional outcomes.

  Study EGAS00001006947

• HLA-DR is absent in primitive macrophages through epigenetic silencing of master regulator CIITA

  The earliest macrophages are generated during embryonic development from erythro-myeloid progenitors (EMPs) via primitive haematopoiesis. This process is poorly understood in humans and in the mouse is thought to be spatially restricted to the yolk sac. Human fetal placental macrophages, Hofbauer cells (HBC), arise during the primitive haematopoietic wave at ~18 days post conception. Here we identify a population of placental erythro-myeloid progenitors (PEMPs) in the early human placenta that give rise to HBC. PEMP are fetal CD34+CD43+ progenitors found exclusively at early gestational timepoints. Transcriptomic analyses reveal that PEMP have conserved features of primitive yolk sac EMPs, including the lack of HLF expression. Using in vitro single-cell culture experiments we show that PEMP generate HBC-like cells that completely lack HLA-DR expression  and demonstrate that this is a conserved feature of all macrophages generated by primitive haematopoiesis in humans. Finally, we demonstrate that the class II transactivator, CTIIA, the master regulator HLA-DR expression, is epigenetically silenced during primitive haematopoiesis, leading to HLA-DRneg macrophages. These findings demonstrate that HBC are derived locally from PEMP and establish the human placenta as an additional site of primitive haematopoiesis.

  Study EGAS00001006981

• RNAseq

  Familial Mediterranean Fever (FMF) is a prototypical periodic fever syndrome caused by genetic variation in MEFV. While it is known that especially IL-1β responses are dysregulated in FMF, its innate immune landscape has not been comprehensively described. Therefore, we extensively characterized the function of monocytes and neutrophils in patients with FMF in between disease attacks, as well as their circulating inflammatory proteome. We found that monocyte IL-1 and IL-6 production was enhanced following a range of stimulations, in concordance with alterations in the plasma inflammatory proteome. In contrast, neutrophil function in FMF was normal. Additionally, ATAC-seq and RNA-seq analyses revealed important epigenetic and transcriptional reprogramming arguing for monocyte function dysregulation. Interestingly, chromatin-accessibility was down-regulated in genomic regions related to cellular responses. This argues that the primary immune dysregulation in monocytes due to MEFV mutations leads to a chronic inflammatory profile that is subsequently associated with counterregulatory epigenetic and transcriptional changes reminiscent of tolerance (rather than trained immunity) induction. These data increase our understanding of the innate immune changes in FMF, that can be used in the design of novel approaches to manage chronic inflammation of these patients.

  Study EGAS00001007165

• Single-cell profiling of co-cultures of GSCCs and macrophages

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007482

• Massively parallel functional dissection of schizophrenia associated non-coding genetic variants

  Schizophrenia (SCZ) is a severe mental disorder affecting 1% of the world population. SCZ is characterized by an underlying genetic architecture that is highly polygenic. Genome wide association studies have identified thousands of genetic variants that are statistically linked to the disease. However, the translation of these associations into insights on the pathomechanisms has been challenging because the causal genetic variants, their molecular function, and their target genes remain largely unknown. To address these questions, we combined induced pluripotent stem cell technology with a massively parallel variant annotation pipeline (MVAP) to functionally characterize 35,000 SCZ associated non-coding genetic variants. This approach identified a set of 620 (1.7%) single nucleotide polymorphisms as functional on the molecular level in a highly cell type and condition specific fashion. Subsequent multi-modal integration of epigenomic data combined with CRISPR screening in human neurons enabled us to systematically translate SCZ variant associations into target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a new high-resolution map of functional variant-gene combinations and offer comprehensive biological insights into the developmental context and stimulus dependent molecular processes modulated by SCZ genetic variation beyond statistical association.

  Study EGAS00001007542

• Immune and clinicopathological features predict HER2-positive breast cancer prognosis in the neoadjuvant NeoALTTO and CALGB 40601 trials

  The identification of prognostic drivers in patients receiving neoadjuvant therapy is crucial for treatment optimization in HER2-positive breast cancer, with the immune microenvironment being a key factor. Here, we investigated the complexity of B and T cell receptor (BCR and TCR) repertoires, reflecting antigen-specific immune response, in the context of two phase III trials, namely NeoALTTO and CALGB 40601, evaluating neoadjuvant paclitaxel with trastuzumab and/or lapatinib in women affected by HER2-positive breast cancer. BCR/TCR features were heterogeneous, and BCR measures describing clonal expansion demonstrated independent prognostic value. A model developed in NeoALTTO and validated in CALGB 40601 predicted event-free survival (EFS) by integrating hormone receptor and clinical nodal status, treatment response (breast pCR)breast pathological complete response (pCR) as well as tumor-infiltrating lymphocytes and BCR evenness, measuring the equality of BCR clonal distribution. The HER2-EveNT score derived from the model allows the identification of patients with 5-year EFS > 90%, and, in those not achieving pCR, of a subgroup of immune-enriched tumors with an excellent outcome despite residual disease.

  Study EGAS00001007563

• Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression

  Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, we analyze 1,532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2,933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n=148) or the alternatively spliced isoform (n=9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal (https://cseminer.stjude.org/), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

  Study EGAS00001007766

• MutWP1: CRUK Grand Challenge Mutographs of Cancer: Oesophageal adenocarcinoma

  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 (40X and 20X depth respectively). 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. This dataset contains all the data available for this study on 2021-09-27.

  Dataset EGAD00001006083

• Genome Diversity in Africa Project: Benin (2021-02-16)

  Globally, human populations show structured genetic diversity as a result of geographical dispersion, selection and drift. Understanding this genetic variation can provide insights into the evolutionary processes that shape both human adaptation and variation in disease. Populations from SSA have the highest levels of genetic diversity. This characteristic, in addition to historical genetic admixture, can lead to complexities in the design of studies assessing the genetic determinants of disease and human variation. However, such studies of African populations are also likely to provide new opportunities to discover novel disease susceptibility loci and variants and refine gene-disease association signals. A systematic assessment of genetic diversity within SSA would facilitate genomic epidemiological studies in the region. The Genome Diversity in Africa Project (GDAP) aims to produce a comprehensive catalogue of human genetic variation in SSA, including single nucleotide polymorphisms (SNPs), structural variants, and haplotypes. This resource will make a substantial contribution to understanding patterns of genetic diversity within and among populations in SSA, as well as providing a global resource to help design, implement and interpret genomic studies in SSA populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out high depth whole genome sequencing of up to 2000 individuals across Africa (25 individuals from each ethnolinguistic group). Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from SSA accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across SSA; 3) develop a cost-efficient, next-generation genotype array for diverse populations across SSA; and 4) facilitate whole genome-sequencing association studies of complex traits and diseases by developing a reference panel for imputation and resource for enhancing fine-mapping disease susceptibility loci. These scientific objectives will be supported by cross-cutting operational activities, including network and management of the consortium, research ethics, and research capacity building in statistical genetics and bioinformatics This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2021-02-16.

  Dataset EGAD00001006970

• Development and first-in-human CAR T therapy targeting the pathognomonic MiT-fusion driven protein GPNMB

  Summary: CAR T therapy for solid tumors is constrained by the scarcity of safe, uniformly expressed cell-surface targets. Here we identify GPNMB, a MiT fusion-driven protein, as highly, homogeneously, and stably expressed in primary and relapsed translocation-positive alveolar soft part sarcoma (ASPS) and renal cell carcinoma (tRCC). We develop a GPNMB-directed CAR T cell product, GCAR1, which demonstrates potent activity against patient-matched cells, organoids and xenograft models. In a first-in-human treatment of two patients with relapsed/refractory, metastatic ASPS, GCAR1 is well tolerated and induces stable disease for up to 6 months, accompanied by resolution of many non-target lesions. GCAR1 T cells expand in peripheral blood as a polyclonal population and remain detectable for up to 6 months. Spatial transcriptomics identify multiple immunosuppressive niches adjacent to T cells infiltrating treatment-resistant lesions; biological PDL1 blockade or PD1 and TIGIT double knockout overcomes this resistance and shows synergistic activity with GCAR1 in xenograft models. Taken together, these data provide clinical proof-of-concept for treating solid tumors with CAR T cells targeting a surface antigen driven by an oncogenic gene fusion.

  Study EGAS50000001696

• Context-specific regulatory genetic variation in MTOR dampens neutrophil-T cell crosstalk in pneumonia-associated sepsis, modulating disease

  Sepsis is a heterogeneous clinical syndrome with a high mortality, requiring personalised stratification strategies. Here, we characterise genetic variation that modulates MTOR, a critical regulator of metabolism and immune responses in sepsis. The effects are highly context specific, involving a regulatory element that affects MTOR expression in activated T cells with opposite direction of effect in neutrophils. We show that the G-allele of the lead variant, rs4845987, associated with decreased risk of T2D, reduces MTOR expression in T cells and improves survival in sepsis due to pneumonia, with effects specific to sepsis endotype. Using ex vivo models, we demonstrate that activated T cells promote immunosuppressive sepsis neutrophils through released cytokines, a process dampened by hypoxia and the mTOR inhibitor rapamycin. Our work demonstrates a novel epigenetic mechanism that fine-tunes MTOR transcription and T cell activity via the variant-containing regulatory element, which further exhibits an allelic effect upon vitamin C treatment. These findings reveal how genetic variation interacts with disease state to modulate immune cell-cell communication, providing a patient stratification strategy to inform more effective sepsis treatment.

  Study EGAS50000000894

• BIOKEY: A single-cell catalogue of the dynamic changes underlying Checkpoint Immunotherapy response in Early Breast Cancer

  Checkpoint immunotherapy combined with neoadjuvant chemotherapy improves complete pathologic response in a subset of breast cancer patients. Here, we applied single-cell profiling to tumor biopsies collected before and during anti-PD1 therapy. One-third of tumors exhibited proliferative T-cells expanding along CD8+ or CD4+ lineages, which were either characterized by increased cytotoxicity and exhaustion or improved T-helper function, respectively. Lineage tracing in non-expanding tumors revealed at which point in the lineage T-cells were impaired, while gene expression modeling along these lineages revealed novel genes and underlying transcription factors involved in T-cell expansion. Interestingly, different dendritic and myeloid cell phenotypes could either stimulate or inhibit expanding T-cells, while cell-to-cell communication revealed an integrated immune context highly predictive of T-cell expansion, consisting of immune-stimulatory/-inhibitory interactions between cancer and various immune cell types. Our data yield unprecedented insights into the dynamic changes underlying checkpoint immunotherapy response in breast cancer.

  Study EGAS00001004809