22929 results for "ega"

in 60.90 milliseconds.

• Spatio-temporal analysis of prostate tumors in situ suggests pre-existence of treatment-resistant clones

  The molecular mechanisms underlying lethal castration-resistant prostate cancer remain poorly understood, with intratumoral heterogeneity a likely contributing factor. To examine the temporal aspects of resistance, we analyzed tumor heterogeneity in needle biopsies collected before and after treatment with androgen deprivation therapy. By doing so, we were able to couple clinical responsiveness and morphological information such as Gleason score to transcriptome-wide data. Our data-driven analysis of transcriptomes identified several distinct intratumoral cell populations, characterized by their unique gene expression profiles. Strikingly, certain cell populations present before treatment exhibited gene expression profiles that matched those of resistant tumor cell clusters, present after treatment. We confirmed that these clusters were resistant by the localization of active androgen receptors to the nuclei in cancer cells post-treatment. Our data also demonstrated that most stromal cells adjacent to resistant clusters do not express the androgen receptor, and we identified differentially expressed genes for these cells. Altogether, this study shows the potential to increase the power in predicting resistant tumors.

  Study EGAS00001006113

• Clonal hematopoiesis is associated with adverse outcomes in patients with COVID-19

  Clonal hematopoiesis of indeterminate potential (CHIP) is defined as the occurrence of an expanded proportion of mature blood cells derived from a mutant hematopoietic precursor without evidence of hematological malignancies. The principle behind this is that the somatic mutation confers a fitness advantage to the cell in which it arose. Different clinical consequences are linked with this expansion. Early evidence of an association with higher mortality risk was provided. This was not related to higher rates of cancer but was associated in particular with increased cardiovascular mortality. Mechanistically, inflammatory processes are not only related to the development of clonal hematopoiesis, but in turn it is also a driver of inflammation. Besides pulmonary symptoms, COVID-19 evokes complex extra-pulmonary manifestations driving the pathophysiology. Among them, both inflammatory and cardiac-associated mechanisms have been deciphered. With the aim of assessing the impact of clonal hematopoiesis on the pathophysiology of COVID-19, hospitalized patients with severe or critical course were evaluated for the presence of CHIP driver mutations and, more importantly, the association with the clinical picture.

  Study EGAS00001006218

• Ovarian cancer organoid biobank

  Ovarian cancer (OC) is a heterogeneous disease usually diagnosed at a late stage. Experimental in vitro models that faithfully capture the hallmarks and tumor heterogeneity of OC are limited and hard to establish. We present a novel protocol that enables efficient derivation and long-term expansion of OC organoids. Utilizing this protocol, we have established 56 organoid lines from 32 patients, representing the spectrum of ovarian neoplasms, including non-malignant borderline tumors, as well as mucinous, clear-cell, endometrioid, low- and high-grade serous carcinomas. OC organoids recapitulate histological and genomic features of the pertinent lesion from which they were derived, illustrating intra- and inter-patient heterogeneity, and can be genetically modified. We show that OC organoids can be used for drug screening assays and capture different tumor subtype responses to the gold standard platinum-based chemotherapy, including acquisition of chemoresistance in recurrent disease. Finally, OC organoids can be xenografted, enabling in vivo drug sensitivity assays. Taken together, this demonstrates their potential application for research and personalized medicine.

  Study EGAS00001003073

• Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures

  Cancer genomes harbor a broad spectrum of structural variants (SV) driving tumorigenesis, a relevant subset of which are likely to escape discovery in short reads. We employed Oxford Nanopore Technologies (ONT) sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assemble complex rearrangements and such associated with telomeric sequences, including a 1.55 Megabasepair chromothripsis event. We uncover a complex SV pattern termed "templated insertion thread", characterized by short (mostly less than 1kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50kbp in size. Templated insertion threads occur in 3% of cancers, with a prevalence ranging to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in seven cancer-driver genes. Our study shows the potential of long-read sequencing in cancer.

  Study EGAS00001006576

• DNA methylation atlas of normal human cell types

  DNA methylation is a fundamental epigenetic mark that governs chromatin organization, cell identity, and gene expression. Here we describe a human methylome atlas, based on deep whole-genome bisulfite sequencing allowing fragment-level analysis across thousands of unique markers for 39 cell types sorted from 205 healthy tissue samples. Replicates of the same cell-type are >99.5% identical, demonstrating robustness of cell identity programs to genetic variation and environmental perturbation. Unsupervised clustering of the atlas recapitulates key elements of tissue ontogeny, and identifies methylation patterns retained since gastrulation. Loci uniquely unmethylated in an individual cell type often reside in transcriptional enhancers and contain DNA binding sites for tissue-specific transcriptional regulators. Uniquely hyper-methylated loci are rare and are enriched for CpG islands, polycomb targets, and CTCF binding sites, suggesting a novel role in shaping cell type-specific chromatin looping. The atlas provides an essential resource for interpretation of disease-associated genetic variants, and a wealth of potential tissue-specific biomarkers for use in liquid biopsies.

  Study EGAS00001006791

• Modulation of the peripheral blood immune cell transcriptome by vitamin D3 supplementation in people with a first demyelinating event: a randomized placebo-controlled trial

  Vitamin D deficiency is a risk factor for developing multiple sclerosis (MS). Both in vitro and animal studies suggest an immunomodulatory effect of vitamin D. The PrevANZ trial, a phase IIb randomized placebo-controlled trial of oral vitamin D3 supplementation in people with a first demyelinating event (FDE), was conducted to determine if supplementation can prevent recurrent disease activity in this cohort at high risk of developing definite MS. As a sub-study of this trial, we used whole blood transcriptomic analyses to investigate the effect of vitamin D3 supplementation on peripheral immune cells in people with an FDE, and to gain insight into potential mechanisms by which vitamin D3 may regulate MS risk and disease activity. The PrevANZ trial randomized participants to 1000 IU, 5000 IU or 10,000 IU daily of oral vitamin D3 or placebo. Peripheral blood was collected at baseline and 12 weeks in PAXgene Blood RNA tubes. Transcriptomic datasets were generated by RNA sequencing.

  Study EGAS00001007346

• DERMATLAS__Leiomyosarcoma_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are leiomyosarcomas.

  Study EGAS00001007631

• DERMATLAS__SG_basal_cell_adenoma_and_adenocarcinoma_WES

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are SG basal cell adenomas and adenocarcinomas.

  Study EGAS00001007745

• Multiple Myeloma follow-up sequencing study

  A total of 40 Multiple Myeloma (MM) patients at clinical relapse who progressed during Proteasome Inhibitors (PIs) or Immunomodulating Drugs (IMiDs)-based therapies and who are assigned to antiCD38-based salvage treatments, will be enrolled. We will collect bone marrow (BM) and peripheral blood (PB) samples from patients at specific timepoints: baseline (BM, PB and buccal swab) every 3 month (PB) achievement of response (≥ Very Good Partial Response (VGPR)) (BM and PB) relapse or refractory status to antiCD38-based treatments (BM and PB) Samples will be processed and stored in the "Hematological Laboratory" located in the University of Turin (Italy) for various proposed analyses: at specific time-points CD138+ (Plasma Cells-PCs) awill be immunomagnetically enriched from the BM mononuclear cells and frozen as viable cells in dimethyl sulfoxide (DMSO); PB mononuclear cells (PBMCs) will be isolated from whole blood by density-gradient centrifugation, and frozen as above; plasma fraction from PB and BM will be obtained by centrifugation and stored frozen; Neutrophils will be obtained at the time of enrollment as a source of control germline DNA and stored frozen.

  Study EGAS00001007092

• Inference of transcription factor binding from cell-free DNA enables tumor subtype prediction and early detection

  Deregulation of transcription factors (TFs) is an important driver of tumorigenesis, but non-invasive assays for assessing transcription factor activity are lacking. We Here we developed and validated a minimally invasive method for assessing TF activity based on cell-free DNA sequencing and nucleosome footprint analysis. We analyzed whole genome sequencing data for >1,000 cell-free DNA samples from cancer patients and healthy controls using a newly developed bioinformatics pipeline developed by us that infers accessibility of TF binding sites from cell-free DNA fragmentation patterns. We observe patient-specific as well as tumor-specific patterns, including accurate prediction of tumor subtypes in prostate cancer, with important clinical implications for the management of patients. Furthermore, we show that cell-free DNA TF profiling is capable of detection of early-stage colorectal carcinomas. Our approach for mapping tumor-specific transcription factor binding in vivo based on blood samples makes a key part of the noncoding genome amenable to clinical analysis.

  Study EGAS00001003206

• DERMATLAS__Leiomyosarcoma_WES

  irtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are leiomyosarcomas.

  Study EGAS00001007628

• DERMATLAS__Porocarcinoma_WES

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are porocarcinomas.

  Study EGAS00001005720

• DERMATLAS__Porocarcinoma_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are porocarcinomas.

  Study EGAS00001005721

• Extensive three-dimensional intratumor proteomic heterogeneity revealed by multiregion sampling in high-grade serous ovarian tumor specimens

  Enriched tumor epithelium, tumor-associated stroma, and whole tissue were collected by laser microdissection from thin sections across spatially separated levels of ten high-grade serous ovarian carcinomas (HGSOCs) and analyzed by mass spectrometry, reverse phase protein arrays, and RNA sequencing. Unsupervised analyses of protein abundance data revealed independent clustering of an enriched stroma and enriched tumor epithelium, with whole tumor tissue clustering driven by overall tumor “purity.” Comparing these data to previously defined prognostic HGSOC molecular subtypes revealed protein and transcript expression from tumor epithelium correlated with the differentiated subtype, whereas stromal proteins (and transcripts) correlated with the mesenchymal subtype. Protein and transcript abundance in the tumor epithelium and stroma exhibited decreased correlation in samples collected just hundreds of microns apart. These data reveal substantial tumor microenvironment protein heterogeneity that directly bears on prognostic signatures, biomarker discovery, and cancer pathophysiology and underscore the need to enrich cellular subpopulations for expression profiling. (Companion proteomic data is available in the ProteomeXchange Consortium, PRIDE: PXD026570).

  Study EGAS00001005786

• Salivary Gland Cancer TSO500 dataset

  Salivary gland cancer (SGC) is a rare cancer for which systemic treatment options are limited. Therefore, it is important to characterize its genetic landscape in search for actionable aberrations. This research aimed to identify these actionable aberrations by NGS-based analysis of DNA (single and multiple nucleotide variants, copy number variants, microsatellite instability and tumor mutational burden) in a large cohort of SGC patients. DNA was extracted from archival tissue of 119 patients with various SGC subtypes and sequenced using a targeted NGS panel encompassing 523 cancer related genes (TruSight Oncology 500, TSO500). The 119 patients included 46 AdCC patients, 43 SDC patients, 15 MEC patients, 9 AciCC patients and 6 patients with other subtypes (1 secretory carcinoma, 1 polymorphic adenocarcinoma (PAC), 1 adenocarcinoma NOS, 1 myoepithelial carcinoma, 1 epithelial/myoepithelial carcinoma and 1 mixed PAC/myoepithelial carcinoma). Median exon coverage ranged from 80-904 (median 308) unique reads and the median percentage of exon coverage with at least 100 unique reads was median 96.6% in these 119 samples.

  Study EGAS00001006232

• scRNASeq of human innate lymphoid cells from different compartments

  Group 1 innate lymphoid cells (ILCs) comprise a heterogeneous family of cytotoxic natural killer (NK) cells and ILC1s. We identifiy a population of “liver-type” ILC1s with transcriptional, phenotypic, and functional features distinct from those of conventional and liver-resident NK cells as well as from other previously described human ILC1 subsets. LT-ILC1s were CD49a+CD94+CD200R1+, expressed the transcription factor T-BET, do not express the activating receptor NKp80, or the transcription factor EOMES. Similar to NK cells, liver-type ILC1s produce IFN-γ, TNF-α, and GM-CSF; however, liver-type ILC1s also produce IL-2 and lack perforin and granzyme B. Liver- type ILC1s are expanded in cirrhotic liver tissues, and they can be produced from blood-derived ILC precursors in vitro in the presence of TGF-β1 and liver sinusoidal endothelial cells. Cells with similar signature and function can also be found in tonsil and intestinal tissues. Collectively, our study identifies and classifies a population of human cross-tissue ILC1s.

  Study EGAS00001006847

• Genotype data from 'Dense sampling of ethnic groups within African countries reveals fine-scale genetic structure and extensive historical admixture.'

  Previous studies have highlighted how African genomes have been shaped by a complex series of historical events. Despite this, genome-wide data has only been obtained from a small proportion of present-day ethnolinguistic groups. By analysing new autosomal genetic variation data of 1333 individuals from over 150 ethnic groups from Cameroon, Republic of the Congo, Ghana, Nigeria and Sudan, we demonstrate a previously underappreciated fine-scale level of genetic structure within these countries, for example correlating with historical polities in western Cameroon. By comparing genetic variation patterns among populations, we infer that many northern Cameroonian and Sudanese groups share genetic links with multiple geographically disparate populations, likely resulting from long-distance migrations. In Ghana and Nigeria, we infer signatures of intermixing dated to over 2000 years ago, corresponding to reports of environmental transformations possibly related to climate change. We also infer recent intermixing signals in multiple African populations, including Congolese, that likely relates to the expansions of Bantu language-speaking peoples.

  Study EGAS00001006944

• This study explored and validated the clinical application of targeted NGS of circulating tumor DNA in identifying tumor-specific mutations and uncovering clinical actionable targets in a variety of solid tumors in large patient cohorts.

  Cancer is a disease of genetic alterations. Comprehensive genetic diagnosis is needed to match each patient to appropriate cancer therapy. However, acquisition of representative tumor samples is invasive and often infeasible. Circulating tumor DNA (ctDNA) is a promising non-invasive biomarker for cancer mutational profiling. Here we implemented targeted next generation sequencing (NGS) with a customized pan-cancer gene panel on 605 clinical ctDNA samples in multiple cancer types. Overall, tumor-specific mutations were identified in 87% of ctDNA samples, with mutation spectra highly concordant with their matched tumor tissues. 71% of patients have at least one clinical actionable mutation with 76% of which have suggested drugs approved or in clinical trials. In particular, our study reveals a unique mutation spectrum in Chinese lung cancer patients, which could be used to guide treatment decision and monitor the appearance of drug-resistant mutations. Taken together, our study demonstrated the clinical utility of target NGS-based ctDNA mutational profiling to guide cancer treatment decision.

  Study EGAS00001002251

• Integrated single cell analysis in transformed follicular lymphoma

  Histological transformation of follicular lymphoma to aggressive form occurs 2-3% per year with poor outcome. Divergent evolution and an altered tumour-microenvironment (TME) have been implicated during transformation. However, phenotypic consequences of this evolution and its implication in TME remain unknown. To address this, we performed single cell whole genome and whole transcriptome sequencing of paired pre/post transformation patient samples. We further performed scWTS of additional samples from patients without transformation. Our analysis revealed evolutionary dynamics of transformation at unprecedented resolution. Integration of scWGS and scWTS identified pathways upregulated during evolution. scWTS analysis revealed a shifting TME landscape, with an exhausted signature emerging during transformation. Using multi-color immunofluorescence we transferred these findings to a TME-based transformation biomarker, subsequently validated in independent pretreatment cohorts. Taken together, our results provide a comprehensive view of the combined genomic and phenotypic evolution of malignant cells during transformation, and shifting cross-talk between malignant cells and TME.

  Study EGAS00001007023

• RNA sequencing data of 257 samples from 106 patients with HR+/HER2- breast cancer treated with AC plus paclitaxel or letrozole plus ribociclib (SOLTI-1402 CORALLEEN trial)

  Early-stage Luminal B breast cancer is frequent and is a major cause of breast cancer death due to its poor prognosis. Our proposal aims to study the biology behind the sensitivity and resistance of Luminal B breast cancer to chemotherapy (CHT) or a non-CHT regimen composed of hormone therapy in combination with ribociclib, a CDK4/6 inhibitor. To accomplish this, we first completed the SOLTI-1402 CORALLEEN phase II trial, a study where 106 patients with early-stage Luminal B breast cancer were randomized to standard neoadjuvant CHT for 6 months, or neoadjuvant letrozole and ribociclib for 6 months. After treatment, patients underwent surgery. The primary results of the study, which showed that the response rate to letrozole+ribociclib was similar to CHT, was reported (Prat et al; Lancet Oncol). Tumor biopsies were available at baseline, week 3 and surgery. A total of 257 samples were analyzed using the Illumina TruSeq Stranded Total RNA w/Ribo Zero Gold with MiSeq in TGL (Sequencer NovaSeq S4/PE/100x).

  Study EGAS00001007060

• DERMATLAS__SG_basal_cell_adenoma_and_adenocarcinoma_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are SG basal cell adenomas and adenocarcinomas.

  Study EGAS00001007746

• Selective targeting of TBXT with DARPins identifies regulatory networks and therapeutic vulnerabilities in chordoma

  The embryonic transcription factor TBXT (brachyury) drives chordoma, a spinal neoplasm without effective drug therapies. TBXT’s regulatory network is poorly understood, and strategies to disrupt its activity for therapeutic purposes are lacking. We developed designed ankyrin repeat proteins that block TBXT-DNA binding (T-DARPins). In chordoma cells, T-DARPins reduced cell cycle progression, spheroid formation, and tumor growth in mice and induced signs of senescence and differentiation. Transcriptomic and proteomic analyses identified gene networks involved in cell cycle regulation, embryonic cell identity, and interferon response and revealed features of regulome components, such as druggability for future translational studies and the fine-tuning of TBXT downstream effectors through IGFBP3. Finally, we found high interferon signaling in chordoma cell lines and patient tumors, which was promoted by TBXT and associated with sensitivity to JAK2 inhibitors. These findings demonstrate the potential of DARPins for probing nuclear proteins to understand the regulatory networks of transcription factor-driven cancers, including entry points for therapies that warrant testing in patients.

  Study EGAS00001008140

• Mutational signature in colorectal cancer induced by genotoxic pks+ E. coli

  Various species of the intestinal microbiota have been associated with the development of colorectal cancer (CRC), yet a direct role of bacteria in the occurrence of oncogenic mutations has not been established. Escherichia coli can carry the pathogenicity island pks, which encodes a set of enzymes that synthesize colibactin. This compound alkylates DNA on adenine residues and induces double strand breaks in cultured cells. Here, we exposed human intestinal organoids to genotoxic pks+ Escherichia coli by repeated luminal injection over a period of 5 months. Whole genome sequencing (WGS) of clonal organoids before and after this exposure reveals a distinct mutational signature, absent from organoids injected with isogenic pks-mutant bacteria. The same mutational signature is detected in a subset of 3668 human metastatic cancer genomes, predominantly in a subset of CRC cases. Our study describes a distinct mutational signature in CRC and implies that the underlying mutational process directly results from past exposure to bacteria carrying the colibactin-producing pks pathogenicity island.

  Study EGAS00001003934

• The_genetics_of_thinness_compared_to_obesity

  The variation in weight within a shared environment is largely attributable to genetic factors. Whilst many genes/loci confer susceptibility to obesity, little is known about the genetic architecture of thinness. In this study we performed a genome-wide association study of 1,622 persistently thin healthy individuals (STILTS), 1,985 severe childhood onset obesity cases (SCOOP) and 10,433 population based individuals (UKHLS) used as a common set of controls. All participants were genotyped on the Illumina Core Exome array, including 551,839 markers and imputed to the combined UK10K and 1000G (phase3) reference panel. We contrast the genetic architecture of thinness with that of severe early onset obesity and explore whether the genetic loci influencing thinness are the same as those influencing obesity pr whether there are important genetic differences between them. 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 EGAS00001002624

• Endocrine therapy reprogramming of breast cancer facilitates metastatic escape via upregulation of P-Rex1/Rac1 signalling

  The estrogen receptor (ER) drives growth in ~70% of breast cancers, and endocrine therapies prevent recurrence. However, 30% of patients recur, with many fatal cases occurring late (after >5 years), presenting a challenge due to their often non-proliferative nature and limited options. To study this, we generated slow-growing resistant cells that form small primary tumours but readily metastasize. Single-cell RNA sequencing (scRNAseq) revealed that endocrine therapy reprograms these cells, notably upregulating the Rac1 signalling component P-Rex1. We found in clinical cohorts that P-Rex1 in high in ER+ breast cancer. Intravital imaging demonstrated that Rac1 signalling is active in ER+ cells following endocrine therapy. Targeting the Rac1 pathway with small molecule inhibitors (NSC23766, R-ketorolac) reduced survival and motility in resistant cells, inhibited in vivo Rac1 activity, and reduced tumour burden when combined with tamoxifen in a drug-refractory patient-derived xenograft. This work identifies the P-Rex1/Rac1 axis, as a potential therapeutic target for late recurring ER+ breast cancer.

  Study EGAS00001008353

• Field_effect_of_healthy_and_diseased_livers

  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 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 EGAS00001002382

• Integration of human pancreatic islet genomic data refines regulatory mechanisms at Type 2 Diabetes susceptibility loci

  Human genetic studies have emphasised the dominant contribution of pancreatic islet dysfunction to development of Type 2 Diabetes (T2D). However, limited annotation of the islet epigenome has constrained efforts to define the molecular mechanisms mediating the, largely regulatory, signals revealed by Genome-Wide Association Studies (GWAS). We characterised patterns of chromatin accessibility (ATAC-seq, n=17) and DNA methylation (whole-genome bisulphite sequencing, n=10) in human islets, generating high-resolution chromatin state maps through integration with established ChIP-seq marks. We found enrichment of GWAS signals for T2D and fasting glucose was concentrated in subsets of islet enhancers characterised by open chromatin and hypomethylation, with the former annotation predominant. At several loci (including CDC123, ADCY5, KLHDC5) the combination of fine-mapping genetic data and chromatin state enrichment maps, supplemented by allelic imbalance in chromatin accessibility pinpointed likely causal variants. The combination of increasingly-precise genetic and islet epigenomic information accelerates definition of causal mechanisms implicated in T2D pathogenesis.

  Study EGAS00001002592

• Analysis of translatome, truncating mutations, lncRNA, circRNA and microproteins of 80 human DCM cases and controls

  The regulation of gene expression in healthy and diseased human tissues is incompletely understood and has primarily been studied on the transcriptional level. Here, we quantify the impact of translation on gene expression and annotate novel translational events in 80 human hearts. These include 65 tissues from patients with dilated cardiomyopathy (DCM) - the most common cause of heart transplantation – and 15 unaffected controls. We show that translational regulation modulates disease-associated molecular processes. Further, we evaluate ribosome dynamics at protein truncating variants revealing incomplete termination of translation at titin-truncating variants (TTNtv) that commonly cause DCM. Finally, we identify dozens of unknown microproteins translated from lncRNAs and circRNAs and link these to mitochondrial processes. Several microproteins are expressed from lncRNAs with previously demonstrated function, suggesting a dual coding and noncoding role. This resource will stimulate novel functional investigations into cardiac biology in health and disease and may serve as a blueprint for characterizing the translational landscape of other human tissues.

  Study EGAS00001003263

• Immunogenomics of colorectal cancer response to immune checkpoint blockade

  Molecular and cellular determinants of response to immune checkpoint inhibitors are mostly unknown. Here we use whole exome, RNA and T cell receptor sequencing, imaging mass cytometry, immunohistochemistry and PD1-PDL1 interaction detection to profile 543 regions from 16 colorectal cancers (CRCs) subsequently treated with Pembrolizumab or Nivolumab. Independent of treatment response, CRCs express low levels of PD1 and PDL1 proteins and form less PD1-PDL1 complex than other cancers. Non hypermutated CRCs have no benefit from treatment, show WNT activation and low T cell infiltrates. Among hypermutated CRCs, those with durable benefit have clonal immunogenic alterations, clonally expanded T cells and potential for immune escape through dysfunctional antigen presentation. Responsive tumours are enriched in PD1-expressing cytotoxic and proliferating CD8 T cells and PDL1-expressing antigen presenting macrophages, which form high-density clusters of interacting cells. Our study shows that immune checkpoint inhibitors are most effective in highly infiltrated CRCs where they may release the interactions between macrophages and CD8 T cells thus promoting their priming and expansion in intra-tumour niches.

  Study EGAS00001004438

• Analysis of mechanisms of CD19- relapse following novel low affinity CD19 Chimeric Antigen Receptor (CAR) T-cells (CD19 CAR T-cells) in a Phase I clinical study in paediatric ALL: CARPALL

  Phase I clinical study (CARPALL, NCT02443831) utilizing a novel low affinity CD19CAR CAR in patients age <25 with high risk CD19+ ALL. Seventeen patients were enrolled and 14 received an infusion of CAR T cells. Among those achieving complete remission, 6 subsequently relapsed, 5 with CD19- disease, 1 with CD19+ disease. Whole exome sequencing of bone marrow DNA from patients with CD19- relapse showed mutations in the CD19 gene predicted to result in loss of surface expression in 4/5 patients. These mutations were not detectable prior to CAR T cell therapy. With a median follow-up of 14 months, 5/14 patients (36%) are alive and disease-free. Overall survival was 84% at 6 months and 63% at 12 months and by the criteria used in the ELIANA study, event-free survival was 67% and 46%. Using more stringent criteria where an event is defined as molecular relapse, molecular EFS was 55% and 29 %. The median duration of remission in responding patients was 7.4 months.

  Study EGAS00001003733

• Population_dynamics_in_abnormal_haematopoiesis

  Background: Our recent study has made direct insight into the clonal dynamics of blood cell production in an unperturbed setting in a human. Using whole genome sequencing to track acquired somatic mutations, we were able to construct a phylogenetic tree of different stem and progenitor cells and used this tree to determine the dynamics of mature blood cell production in longitudinal samples. These data set the baseline for understanding population dynamics in abnormal haematopoiesis, and how malignant clones outcompete their normal counterparts to drive haematological malignancies. Aim: This study will measure population dynamics in abnormal human haematopoiesis in patients with clonal blood stem cell disorders. Methods: From patient bone marrow or peripheral blood, we have isolated colonies grown from single blood stem cells/progenitors. Colony DNA will be analysed by whole-genome sequencing to build (or populate) a phylogenetic tree to determine clonal relationships and track stem cell contribution to mature blood cell production. These data will also reveal disease-associated genomic features (mutational burden and mutational signatures) in haematopoietic cells.

  Study EGAS00001003181

• Somatic pathogenic variants in the normal mammary gland of sporadic breast cancer patients.

  The mammary gland undergoes hormonally stimulated cycles of proliferation, lactation and involution. We hypothesized that these factors increase the mutational burden in glandular tissue and may explain high cancer incidence rate in the general population and recurrent disease. Hence, we investigated the DNA sequence variants in the normal mammary gland, tumor and peripheral blood from 52 reportedly sporadic breast cancer patients, including breast-conserving surgery cases. Targeted resequencing of 542 cancer associated genes revealed mosaic somatic pathogenic variants of: PIK3CA, TP53, AKT1, MAP3K1, CDH1, RB1, NCOR1, MED12, CBFB, TBX3 and TSHR in the normal mammary gland, at considerable allelic frequencies (9x10-2 to 5.2x10-1) indicating clonal expansion. Further evaluation of the frequently damaged PIK3CA and TP53 genes by ultra-sensitive duplex sequencing demonstrated a diversified picture of multiple low level-mosaic (in 10-2 to 10-4 alleles) hotspot pathogenic variants. Our results raise a question about the oncogenic potential in non-tumor mammary gland tissue of breast-conserving surgery patients.

  Study EGAS00001005698

• ALPI deficiency causes refractory Inflammation Bowel Disease

  Herein, we report the first identification of biallelic-inherited mutations in ALPI as a Mendelian cause of inflammatory bowel disease in two unrelated patients. ALPI encodes for intestinal phosphatase alkaline, a brush border metalloenzyme that hydrolyses phosphate from the lipid A moiety of lipopolysaccharides and thereby drastically reduces Toll-like receptor 4 agonist activity. Prediction tools and structural modelling indicate that all mutations affect critical residues or inter-subunit interactions, and heterologous expression in HEK293T cells demonstrated that all ALPI mutations were loss of function. ALPI mutations impaired either stability or catalytic activity of ALPI and rendered it unable to detoxify lipopolysaccharide-dependent signalling. Furthermore, ALPI expression was reduced in patients' biopsies, and ALPI activity was undetectable in ALPI-deficient patient's stool. Our findings highlight the crucial role of ALPI in regulating host-microbiota interactions and restraining host inflammatory responses. These results indicate that ALPI mutations should be included in screening for monogenic causes of inflammatory bowel diseases and lay the groundwork for ALPI-based treatments in intestinal inflammatory disorders.

  Study EGAS00001003350

• Somatic mutations reveal lineage relationships and age-related mutagenesis in human hematopoiesis

  Mutation accumulation during human life can contribute to hematopoietic dysfunction; however, the underlying dynamics are unknown. Somatic mutations in blood progenitors can provide insight into the rate and processes underlying this accumulation, as well as the developmental lineage tree and stem cell division numbers. Here, we catalogue somatic mutations in the genomes of human bone marrow-derived and cord blood-derived hematopoietic stem and multipotent progenitor cells (HSPCs) and find that base substitutions accumulate with approximately 17 base substitutions per year in both populations, while insertions and deletions occur sporadically and at low numbers. The majority of mutations in adult HSPCs were acquired after birth and could be explained by the constant activity of various endogenous processes, which also explains the mutation load in acute myeloid leukemia (AML). We construct a developmental lineage tree revealing a polyclonal architecture of the hematopoietic progenitor compartment and providing evidence that developmental clones exhibit multipotency, though lineage biases can be detected. Our approach highlights novel features of human native hematopoiesis and its implications for leukemogenesis.

  Study EGAS00001003068

• Uveal melanoma patient with germline MBD4 nonsense mutation

  There is currently no effective treatment for metastasised uveal melanoma (UM). Recently, it was reported that a UM patient was responsive to checkpoint inhibitor (CI) treatment, due to a high tumour mutation burden correlated with a germline loss-of-function MBD4 mutation. Here, we report on another UM patient who carried an MBD4 germline nonsense variant (p.Leu563Ter) and her tumour showed a 5-fold higher than average mutation burden. We confirmed the association between germline loss-of-function variant in MBD4 and CI response. The patient experienced stable disease (10 months) and survived two years with metastatic disease, which is twice as long as median survival. Additionally, the frequency of MBD4 loss-of-function variants in reported UM cohorts was >20 times higher than in an aggregated population genome database (P < 5x10-5), implying a potential role as UM predisposition gene. These findings provide a strong basis for the inclusion of MBD4 in screening of potential UM-prone families as well as stratification of immunotherapy.

  Study EGAS00001003362

• Clinical Cancer Genomic Profiling by Three-Platform Sequencing of Whole Genome, Whole Exome and Transcriptome

  To evaluate the potential of an integrated clinical test to detect diverse classes of somatic and germline mutations relevant to pediatric oncology, we performed three-platform whole-genome (WGS), whole exome (WES) and transcriptome (RNA-Seq) sequencing of tumors and normal tissue from 78 pediatric cancer patients in a CLIA-certified, CAP-accredited laboratory. Our analysis pipeline achieves high accuracy by cross-validating variants between sequencing types, thereby removing the need for confirmatory testing, and facilitates comprehensive reporting in a clinically-relevant timeframe. Three-platform sequencing has a positive predictive value of 97-99%, 99% and 91% for somatic SNVs, indels and structural variations, respectively, based on independent experimental verification of 15,225 variants. We report 240 pathogenic variants across all cases, including 84 of 86 known from previous diagnostic testing (98% sensitivity). Combined WES and RNA-Seq, the current standard for precision oncology, achieved only 78% sensitivity. These results emphasize the critical need for incorporating WGS in pediatric oncology testing.

  Study EGAS00001002217

• Drug screening and whole genome sequencing of primary cells and cell lines from ovarian cancer patients to associate genomic aberrations with in vitro drug sensitivities

  Identifying drugs for cancers with poor prognosis like ovarian cancer (OC) is hindered by inappropriate cell models and the limitations of current viability assays. We developed an automated confocal microscopy-based live cell assay (DeathPro) that robustly de-convolves cell death and proliferation over time by employing a specific image processing strategy. Using DeathPro, we screened nine OC patient-derived cell lines and one normal cell line (HOSEpiC) with 22 drugs or combinations. Whole genome was performed to analyze if the cell lines harbor OC-typical genomic aberrations and to relate genomic features to drug sensitivities. Commonly observed alterations in OC patients like TP53 mutations and MYC or PI3KCA amplifications were absent in the normal cell line but frequently detected in the OC cell lines. Moreover, genome-derived homologous recombination deficiency (HRD) scores of the OC cell lines correlated with their responses to carboplatin, olaparib and DNA demethylation drugs. Overall, we propose an imaging-based platform to guide drug development and genome-drug sensitivity analysis.

  Study EGAS00001002239

• Analysis of four key cell types (epithelial, fbroblast, myeloid and T cells)

  Alterations of the normal prostate epithelial cell genome trigger the emergence of carcinoma, however cancer-only analyses of genomic aberrations have resulted in little clinical translation. A deep understanding of the tumor microenvironment is likely to inform on key drivers of disease progression and novel therapeutic opportunities. We investigated the role of key cell types along disease progression, profiling the RNA abundance of 52 samples from 13 patients enriched for one of four key cell types, spanning a wide CAPRA-S risk score range. Using a novel differential transcription Bayes method, TABI, we avoided a priori patient stratification, and mapped transcriptional alteration events to specific risk scores. The developmental and cell-type resolution yielded a landscape of concurrent transcriptional alterations from the tumour microenvironment, including key hallmarks. These transcriptional signatures permitted the identification of a “high survival” patient cluster in an independent prostate cancer dataset. Most importantly, together with our integrated transcriptional analysis, differential tissue composition analysis provided evidence for a key role for monocytes and macrophages in prostate cancer progression and disease recurrence.

  Study EGAS00001003579

• CD36 defines CML cells less sensitive to imatinib

  Tyrosine kinase inhibitors (TKIs) are highly effective for treatment of chronic myeloid leukemia (CML), but very few patients are cured. Major drawbacks with TKIs are their low efficacy in eradicating the leukemic stem cells responsible for disease maintenance and relapse upon TKI cessation. Here, we performed RNA-sequencing of flow-sorted primitive (CD34+CD38low) and progenitor (CD34+CD38+) chronic phase CML cells and identified transcriptional upregulation of 32 cell surface molecules relative to corresponding normal bone marrow cells. Focusing on novel markers with increased expression on primitive CML cells, we confirmed upregulation of the scavenger receptor CD36 and the leptin receptor (LEPR) by flow cytometry. We also delineate a subpopulation of primitive CML cells expressing CD36 that is less sensitive to imatinib treatment. Using CD36 targeting antibodies, we show that the CD36 positive cells can be targeted and killed by antibody dependent cellular cytotoxicity (ADCC). In summary, CD36 defines a subpopulation of primitive CML cells with decreased imatinib sensitivity that can be effectively targeted by ADCC using an anti-CD36 antibody.

  Study EGAS00001002421

• Genomic landscape of human diversity across Madagascar

  Although situated ~400 km from the east coast of Africa, Madagascar exhibits cultural, linguistic, and genetic traits from both Southeast Asia and Eastern Africa. The settlement history remains contentious; we therefore used a grid-based approach to sample at high-resolution the genomic diversity (including maternal lineages, paternal lineages, and genome-wide data) across 257 villages and 2704 Malagasy individuals. We find a common Bantu and Austronesian descent for all Malagasy individuals with a limited paternal contribution from Europe and the Middle East. Admixture and demographic growth happened recently, suggesting a rapid settlement of Madagascar during the last millennium. However, the distribution of African and Asian ancestry across the island reveals that the admixture was sex biased and happened heterogeneously across Madagascar, suggesting independent colonization of Madagascar from Africa and Asia rather than settlement by an already-admixed population. In addition, there are geographic influences on the present genomic diversity, independent of the admixture, showing that a few centuries is sufficient to produce detectable genetic structure in human populations.

  Study EGAS00001002549

• TGF-β attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells.

  We examined tumours from a large cohort of patients with metastatic urothelial bladder cancer (mUC) treated with an anti–PD-L1 agent (atezolizumab) and identified major determinants of clinical response and immune escape. Whole transcriptome profiles were generated for 368 patients using TruSeq RNA Access technology (Illumina). Somatic mutations were determined via whole exome sequencing for a subset of these patients. - Response was associated with CD8+ T effector cell phenotype and, to an even greater extent, high neoantigen or tumour mutation burden (TMB). On the other hand, lack of response was associated with a transforming growth factor β (TGF-β) signature, particularly in patients with CD8+ T cells preferentially residing in collagen-rich matrix surrounding tumours. Integration of these three independent biological features provided the best basis for understanding outcome in this setting. -

  Study EGAS00001002556

• Genetic and Functional Drivers of Diffuse Large B Cell Lymphoma

  Diffuse large B cell lymphoma (DLBCL) is the most common hematologic malignancy and is characterized by a striking degree of heterogeneity. We enrolled a total of 1001 DLBCL patients and comprehensively defined the landscape of genetic mutations, copy number alterations and expression through whole-exome and transcriptome sequencing. We identified 150 genetic drivers of DLBCL including many novel, clinically relevant genes (e.g. SPEN, SETD1B and KLHL14). Genetic drivers were highly enriched among essential genes identified by CRISPR screening in DLBCL cell lines including those of immediate potential therapeutic relevance such as MTOR, SYK, SF3B1 and XPO1. RHOA emerged as a critical driver gene from our genetic analysis and CRISPR screen; we further delineated its functional role in vitro and in vivo models to define its role as an essential oncogene in DLBCL. Our work thus identifies the functional landscape of genetic drivers and their clinical and therapeutic relevance in DLBCL patients.

  Study EGAS00001002606

• Sequencing of an organoid biobank for childhood kidney cancers that captures disease and tissue heterogeneity.

  Kidney tumors are among the most common solid tumors in children, comprising several distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first pediatric cancer organoid biobank. It contains tumor and matching normal organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumors, malignant rhabdoid tumors, renal cell carcinomas, and congenital mesoblastic nephromas. The malignant rhabdoid tumor organoids represent the first organoid model for tumors of non-epithelial origin. The tumor organoids retain key properties of native tumors, useful for revealing patient specific drug vulnerabilities. We further demonstrate that organoid cultures derived from Wilms tumors consist of multiple different cell types, including epithelial, stromal and blastemal-like. Our organoid biobank captures the cellular heterogeneity of pediatric kidney tumors, providing a representative collection of well-characterized models for basic cancer research, drug-screening, and personalized medicine.

  Study EGAS00001003853

• The whole blood of female volunteers and sperm from the male volunteer were used to extract genomic DNA to do whole genome sequencing. Distinguished SNPs between parental genomes were retained to analyze parental allele-specific DNA methylation and chromatin accessibility in scCOOL-seq data.

  DNA methylation, chromatin states, and their interrelationships represent critical epigenetic information, but those are largely unknown in human early embryos. Here, we apply scCOOL-seq (Chromatin Overall Omic-scale Landscape Sequencing) to generate a genome-wide map of DNA methylation and chromatin accessibility at single-cell resolution during human preimplantation development. Unlike in mice, the chromatin of paternal genome is already more open than that of maternal genome at the mid-zygote stage in humans, and this state is maintained until the 4-cell stage. After fertilization, genes with high variations in DNA methylation and those with high variations in chromatin accessibility tend to be two different sets. Furthermore, 1,797 (35%) out of 5,155 widely open chromatin regions in promoters closed when transcription activity was inhibited, indicating a feedback mechanism between transcription and open chromatin maintenance. Our work paves the way for dissecting the complex, yet highly coordinated, epigenetic reprogramming during human preimplantation development.

  Study EGAS00001002987

• CD8+ T-cell exhaustion induced by leukemic cells drives progression in Chronic Lymphocytic Leukemia

  In order to shed light on mechanisms that lead to CLL progression, we investigated longitudinal changes in both the genetic and immunological landscapes. For that, we analyzed paired samples obtained from CLL patients at diagnosis and progression before front-line treatment or at diagnosis and long-term asymptomatic follow-up. By whole-exome sequencing (WES), we found that CLL cells displayed limited genetic evolution at progression. Conversely, T cells from progressing patients showed significant immunophenotypical and transcriptional changes over time not observed in patients that did not progress. Specifically, we found that CD8+ T cells at progression were characterized by an increasing terminally exhausted effector phenotype (T-betdim/- Eomeshi PD1hi) with high co-expression of inhibitory receptors (PD1, CD244 and CD160). This exhausted status in CD8+ T cells was induced by progressing leukemic cells mainly through a mechanism dependent on soluble factors, including IL-10. Thus, we demonstrate that immunological changes are of paramount importance for CLL clinical progression, thereby providing a rationale for the use of early immunotherapeutic intervention in this disease.

  Study EGAS00001004116

• Dissecting intratumor heterogeneity of nodal B cell lymphoma on the transcriptional, genetic, and drug response level

  Tumor heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is well-known to affect the efficacy of anti-cancer, most personalized treatment approaches do not account for intratumor heterogeneity. We addressed this issue by studying the heterogeneity of lymph node-derived B cell non-Hodgkin lymphoma (B-NHL) by single cell RNA-sequencing (scRNA-seq) and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subclones and compared their drug response and genomic profiles. Malignant subclones of the same patient responded strikingly different to anti-cancer drugs ex vivo, which recapitulated subclone-specific drug sensitivity during in vivo treatment. Tumor infiltrating T cells represented the majority of non-malignant cells, whose gene expression signatures were similar across all donors, whereas the frequencies of T cell subsets varied significantly between the donors. Our data provide new insights into the heterogeneity of B-NHL and highlight the relevance of intratumor heterogeneity for personalized cancer therapies.

  Study EGAS00001004335

• Molecular subsets in renal cancer determine outcome to checkpoint and angiogenesis blockade

  We conducted integrated multi-omics evaluation of 823 tumors from advanced renal cell carcinoma (RCC) patients and identified molecular subsets associated with differential clinical outcomes to angiogenesis blockade alone or with a checkpoint inhibitor. Unsupervised transcriptomic analysis revealed six molecular subsets with distinct angiogenesis, immune, cell cycle, metabolism, and stromal programs. While sunitinib was effective in subsets with high angiogenesis, atezolizumab+bevacizumab improved clinical benefit in tumors with high T-effector and/or cell cycle transcription. Somatic mutations in PBRM1 and KDM5C associated with high angiogenesis and AMPK/Fatty acid oxidation gene expression, while CDKN2A/B and TP53 alterations were more prevalent in tumors with increased cell cycle and anabolic metabolism. Sarcomatoid tumors exhibited low prevalence of PBRM1 mutations and angiogenesis markers, frequent CDKN2A/B, BAP1 and TP53 alterations, and increased PD-L1 expression. These findings can be applied to molecularly stratify advanced RCC patients, explain improved outcomes of sarcomatoid tumors to checkpoint blockade, and develop personalized therapies.

  Study EGAS00001004353

• Whole-genome sequencing of rare disease patients in a national healthcare system

  Most patients with rare diseases do not receive a molecular diagnosis and the aetiologicalvariants and mediating genes for more than half such disorders remain to be discovered. Weimplemented whole-genome sequencing (WGS) in a national healthcare system to streamlinediagnosis and to discover unknown aetiological variants, in the coding and non-coding regionsof the genome. In a pilot study for the 100,000 Genomes Project, we generated WGS data for13,037 participants, of whom 9,802 had a rare disease, and provided a genetic diagnosis to1,138 of the 7,065 patients with detailed phenotypic data. We identified 95 Mendelianassociations between genes and rare diseases, of which 11 have been discovered since 2015and at least 79 are confirmed aetiological. Using WGS of UK Biobank1, we showed that rarealleles can explain the presence of some individuals in the tails of a quantitative red blood cell(RBC) trait. Finally, we reported 4 novel non-coding variants which cause disease through thedisruption of transcription of ARPC1B, GATA1, LRBA and MPL. Our study demonstrates asynergy by using WGS for diagnosis and aetiological discovery in routine healthcare.

  Study EGAS00001004364

• Hypertension delays viral clearance and exacerbates airway hyperinflammation in patients with COVID-19

  In coronavirus disease 2019 (COVID-19), hypertension and cardiovascular diseases are major risk factors for critical disease progression. However, the underlying causes and the effects of the main anti-hypertensive therapies-angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs)-remain unclear. Combining clinical data (n = 144) and single-cell sequencing data of airway samples (n = 48) with in vitro experiments, we observed a distinct inflammatory predisposition of immune cells in patients with hypertension that correlated with critical COVID-19 progression. ACEI treatment was associated with dampened COVID-19-related hyperinflammation and with increased cell intrinsic antiviral responses, whereas ARB treatment related to enhanced epithelial-immune cell interactions. Macrophages and neutrophils of patients with hypertension, in particular under ARB treatment, exhibited higher expression of the pro-inflammatory cytokines CCL3 and CCL4 and the chemokine receptor CCR1. Although the limited size of our cohort does not allow us to establish clinical efficacy, our data suggest that the clinical benefits of ACEI treatment in patients with COVID-19 who have hypertension warrant further investigation.

  Study EGAS00001004772

• RNA-seq study of a Princess Margaret Cancer Centre human acute myeloid leukemia patient cohort

  The PMCC AML RNAseq dataset consists of 81 AML patient samples (clinical data in Supplemental Table 11 of manuscript), processed in two batches. These patient samples are able to engraft in the NSG (NOD.Cg PrkdcscidIl2rgtm1Wjl /SzJ) mouse model. Five patients (90543, 598, 90240, 110484, 100500) were included in both batches. Viaably frozen material from the Leukemia Tissue Bank at Princess Margaret Cancer Centre/ University Health Network were thawed by dropwise addition of X-VIVO + 50% fetal calf serum supplemented with DNase (100μg/mL final concentration, Roche). RNA was extracted from bulk peripheral blood mononuclear cells (PBMC) using the RNeasy Micro Kit (Qiagen Inc.). A paired-end 76 base-pair flow-cell lane Illumina High seq 2000 yielded an average of 240 million sequence reads aligning to genome per sample at the Genome Sciences Centre, BC Cancer Agency for cohort 1. Cohort 2 was subjected to 125 bp, paired-end RNA-sequencing on the Illumina HiSeq 2500 with an average of 50 million reads/sample at the Centre for Applied Genomics, Sick Kids Hospital.

  Study EGAS00001004792

• The genetic structure of Norway

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

  Study EGAS00001004826

• Liverpool Preterm Birth Biomarker Study

  Preterm birth is a multifactorial condition defined as birth less than 37 weeks of gestation. This study aimed to identify early pregnancy biomarkers from maternal genome-wide and transcriptome-wide data of women experiencing two types of spontaneous PTB (sPTB): spontaneous preterm birth (SPTB) and preterm premature rupture of membranes (PPROM). Whole blood samples were obtained from women with singleton pregnancies (N=567) recruited at the Liverpool Women’s Hospital at 16 and 20 weeks of gestation. Women with a previous history of sPTB pregnancies <34 weeks were categorised on their index pregnancy outcome as spontaneous preterm birth (SPTB) and preterm premature rupture of membranes (PPROM) or high-risk term control (HTERM, ≥37 weeks). Women with no history of SPTB/PPROM and delivered at term were LTERM (≥39 weeks of gestation). DNA was genotyped on the UK Biobank Axiom™ array (Applied Biosystems™, Thermo Fisher Scientific) and RNA who processed using the Clariom™ D Human assay (Thermo Fisher Scientific). Genome-wide association analysis, differential expression analysis and expression quantitative trait loci (eQTL) were completed.

  Study EGAS00001005076

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

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

  Study EGAS00001005396

• Combined Metabolic Activators Reduces Liver Fat in Nonalcoholic Fatty Liver Disease Patients

  Nonalcoholic fatty liver disease (NAFLD) refers to excess fat accumulation in the liver. In animal experiments and human kinetic study, we found that administration of combined metabolic activators (CMA) promotes the oxidation of fat, attenuates the resulting oxidative stress, activates mitochondria and eventually removes excess fat from the liver. Here, we tested the safety and efficacy of CMA in NAFLD patients in a placebo-controlled 10-week study. We found that CMA significantly decreased hepatic steatosis and levels of aspartate aminotransferase, alanine aminotransferase, uric acid, and creatinine, whereas found no differences on these variables in the placebo group after adjustment for weight loss. By integrating clinical data with plasma metabolomics and inflammatory proteomics as well as oral and gut metagenomics data, we revealed the underlying molecular mechanisms associated with the reduced hepatic fat and inflammation in NAFLD patients and identified the key players involved in the host-microbiome interactions. In conclusion, we showed that CMA can be used to develop a pharmacological treatment strategy in NAFLD patients.

  Study EGAS00001005616

• DERMATLAS__Poroma_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are poromas.

  Study EGAS00001005759

• Combination therapies to inhibit LCK tyrosine kinase and mTOR signaling in T-cell Acute Lymphoblastic Leukemia

  Relapse and refractory T cell acute lymphoblastic leukemia (T-ALL) has a poor prognosis and new combination therapies are sorely needed. This work uncovered potent drug synergies between AKT/mTORC1 inhibitors and the general tyrosine kinase-inhibitor, dasatinib. The combination effectively curbed T-ALL growth in human cell lines and primary human T-ALL samples, and was well tolerated and effective in suppressing leukemia growth in patient-derived xenografts grown in mice. Mechanistically, dasatinib inhibited phosphorylation and activation of the lymphocyte-specific protein tyrosine kinase (LCK) to blunt the T-cell receptor (TCR) signaling pathway and when complexed with mTORC1 inhibition, induced potent T-ALL cell killing through reducing MCL-1 protein expression. Analysis of responses across T-ALL subtypes and association with specific genetic mutations or expression patterns failed to identify any obvious biomarkers or correlations in predicting therapy responses, although there appeared to be a trend toward higher TCR/LCK pathway activity in samples that exhibited synergistic therapy responses.

  Study EGAS00001005945

• Biallelic HMBS Inactivation Defines a Homogenous Clinico-Molecular Subtype of Hepatocellular Carcinoma

  Acute intermittent porphyria (AIP), caused by heterozygous germline mutations of the heme synthesis pathway enzyme HMBS (hydroxymethylbilane synthase), confers high risk of hepatocellular carcinoma (HCC) development. Yet, the role of HMBS in liver tumorigenesis remains unclear. Here we explore HMBS alterations in a large series of 758 HCC. We identify recurrent bi-allelic HMBS inactivation, both in AIP patients acquiring a second somatic HMBS mutation and in sporadic HCC with two somatic hits. HMBS alterations are enriched in truncative mutations, in particular in splice regions, leading to abnormal transcript structures. Bi-allelic HMBS inactivation results in a massive accumulation of its toxic substrate porphobilinogen and synergizes with CTNNB1 activating mutations, leading to the development of well differentiated tumors with a transcriptomic signature of Wnt/ß-catenin pathway activation and a DNA methylation signature related to ageing. HMBS-inactivated HCC mostly affect females, in absence of fibrosis and classical HCC risk factors. These data identify HMBS as a tumor suppressor gene whose bi-allelic inactivation defines a homogenous clinical and molecular HCC subtype.

  Study EGAS00001005986

• Genetic history of the Swahili population

  At the dawn of the second millennium, the expansion of the Indian Ocean trading network aligned with the emergence of an outward-oriented community along the East African coast to create a cosmopolitan cultural and trading zone known as the Swahili Corridor. Based on analyses of genome-wide genotyping data in 140 individuals from coastal Kenya and the Comoros islands, along with 3,477 individuals from the Indian Ocean rim, we reconstruct historical population dynamics showing that the Swahili Corridor is largely an East Bantu genetic continuum. Within this continuity significant gene flows from the Middle East can be seen in the Swahili and Comorians at dates corresponding to the Islamization of East Africa. However the main external gene flow into insular populations of the Swahili Corridor, such as Comorian and Malagasy groups, came from Island Southeast Asia as early as the 10th century. Remarkably, our results reflect archaeological and linguistic data in suggesting that the Comoros archipelago is the oldest place of contact between Austronesian and African populations in the Swahili Corridor.

  Study EGAS00001002569

• DERMATLAS__Hidradenoma_papilliferum_WES

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are hidradenoma papilliferums.

  Study EGAS00001005714

• DERMATLAS__Hidradenoma_papilliferum_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are hidradenoma papilliferums.

  Study EGAS00001005715

• Genome-wide mutational consequences of nucleotide excision repair-deficiency through XPC deletion in a human adult stem cell culture

  Nucleotide excision repair (NER) is one of the main DNA repair pathways that protect cells against genomic damage. Deficiency in this pathway can contribute to the development of cancer and accelerate aging. NER-deficiency is an important determinant for cancer treatment outcome, as NER-deficient tumors are selectively sensitive to cisplatin treatment. While NER-deficiency has been linked to mutational Signature 5, not all NER-deficient tumors are characterized by a high Signature 5 contribution, illustrating the importance to further characterize the mutational consequences of NER-deficiency. Here, we analyzed the mutational profile of a human adult stem cell (organoid) culture that is deficient in NER through XPC deletion by CRISPR-Cas9 gene-editing, and subsequent whole-genome sequencing analysis of a clonally derived organoid. We found that XPC deletion results in an increase in base substitution load and specifically induces Signature 8 mutations, a mutational signature with previously unknown etiology. Presence of Signature 8 may, therefore, serve as a novel biomarker for NER-deficient tumors and could improve personalized treatment strategies.

  Study EGAS00001002681

• BAMSE (Swedish abbreviation for Children, Allergy, Milieu, Stockholm, Epidemiology)

  The BAMSE (Swedish abbreviation for Children, Allergy, Milieu, Stockholm, Epidemiology) study is an ongoing longitudinal, population-based prospective birth cohort including 4,089 children born between 1994 and 1996 in Stockholm, Sweden. The cohort was initially designed to study risk factors for asthma, allergic diseases and lung function in childhood, and to study factors of importance for prognosis at already established disease.Questionnaires on respiratory symptoms and medication were answered at age of 1, 2, 4, 8 and 16 years. Response rates ranged from 96-82% at each occasion with very minor selection bias over the years. Exposure to air pollutants and other environmental factors has been mapped since birth.At the 8 and 16 year follow-up, spirometry and FeNO measurements were performed and at the latter, impulse oscillometry was measured. Blood samples (including plasma) from around 2,500 children were taken at 4, 8 and 16 years and have been analyzed for different IgE-ab and biomarkers. Genome-wide genetic, global methylation and transcriptomic data exist on a subset of the children.

  Study EGAS00001002746

• Whole-Exome Sequencing of Salivary Gland Mucoepidermoid Carcinoma

  Mucoepidermoid carcinoma (MEC) is the most common salivary gland malignancy. To explore the genetic origins of MEC, we performed systematic genomic analyses of these tumors.Experimental DesignWhole-exome sequencing and gene copy number analyses were performed for 18 primary cancers with matched normal tissue. Fluorescence in situ hybridization (FISH) was used to determine the presence or absence of the MECT1-MAML2 translocation in 17 tumors.ResultsTP53 was the most commonly mutated gene in MEC (28%), and mutations were found only in intermediate- and high-grade tumors. Tumors with TP53 mutations had more mutations overall than tumors without TP53 mutations (p=0.006). POU6F2 was the second most frequently mutated gene, found in three low-grade MECs with the same in-frame deletion. Somatic alterations in IRAK1, MAP3K9, ITGAL, ERBB4, OTOGL, KMT2C, and OBSCN were identified in at least two of the 18 tumors sequenced. FISH analysis confirmed the presence of the MECT1-MAML2 translocation in 15 of 17 tumors (88%).ConclusionsThrough these integrated genomic analyses, MECT1-MAML2 translocation and somatic TP53 and POU6F2 mutations appear to be the main drivers of mucoepidermoid carcinoma.

  Study EGAS00001002811

• Whole exome sequencing of longitudinal samples from a melanoma patient receiving MEK plus CDK4/6 inhibitor therapy.

  Combined MEK and CDK4/6 inhibition (MEKi+CDK4i) has shown promising clinical outcomes in NRAS mutant melanoma patients. Here, we interrogated longitudinal biopsies from a patient who initially responded to MEKi+CDK4i therapy but subsequently developed resistance. Whole exome sequencing and functional validation identified an acquired PIK3CAE545K mutation as conferring drug resistance. We demonstrate that PIK3CAE545K pre-existed in a rare subpopulation that was missed by both clinical and research testing, but was revealed upon multi-region sampling due to PIK3CAE545K being non-uniformly distributed. This resistant population rapidly expanded after the initiation of MEKi+CDK4i therapy and persisted in all successive samples even after immune checkpoint therapy and distant metastasis. Functional studies identified activated S6K1 as both a key marker and specific therapeutic vulnerability downstream of PIK3CAE545K-induced resistance. These results demonstrate that comprehensive analysis of pre-treatment samples can reveal rare pre-existing resistant subpopulations and also posit S6K1 as a common downstream therapeutic nexus for the MAPK, CDK4/6, and PI3K pathways.

  Study EGAS00001002846

• Using human induced pluripotent stem cells (iPSC) and iPSC-derived neurons to explore cellular phenotypes associated with schizophrenia.

  Schizophrenia and major mental illness are common devastating conditions. Whilst partially effective treatments are available, none are disease-modifying. Notwithstanding important insights provided by animal models there is a need for novel models of mental illness to develop more effective therapies. A study of sporadic schizophrenia-derived human induced pluripotent stem cells (hiPS) has demonstrated a neuronal phenotype representing a powerful proof of concept of in vitro disease modeling, even for diseases with complex causation such as schizophrenia. These techniques have been used here to examine the behaviour of neural stem cells and neurons derived from fibroblasts from subjects carrying NDE1 copy number variations (CNVs) associated with increased risk of major mental illness. Using hiPS from individuals with and without disease-associated CNVs we have investigated how disease risk is conferred at a cellular and molecular level through detailed comparative studies of neural progenitor proliferation, neuronal morphology, physiology, cytoskeletal organization, synaptic function, glutamate (NMDA) receptor expression and also RNA sequencing and proteomic studies of these patient-derived cells.

  Study EGAS00001002895

• Pheno-Seq, linking 3D phenotypes of clonal tumor spheroids to gene expression

  3D-culture systems have advanced cancer modeling by reflecting physiological characteristics of in-vivo tissues, but our understanding of functional intratumor heterogeneity including visual phenotypes and underlying gene expression is still limited. Transcriptional heterogeneity can be dissected by single-cell RNA-sequencing, but these technologies suffer from low RNA-input and fail to directly correlate gene expression with contextual cellular phenotypes. Here we present pheno-seq for integrated high-throughput imaging and transcriptomic profiling of clonal tumor spheroids derived from 3D models of breast and colorectal cancer. Specifically, we identify characteristic expression signatures that are associated with heterogeneous invasive and proliferative behavior including a rare cell subtype. Furthermore, we identify functionally relevant transcriptional regulators missed by single-cell RNA-seq, link visual phenotypes defined by heterogenous expression to inhibitor response and infer single-cell regulatory states by deconvolution. We anticipate that directly linking molecular features with patho-phenotypes of cancer cells will improve the understanding of intratumor heterogeneity and consequently prove to be useful for translational research.

  Study EGAS00001002999

• Multi-omic analysis of cell-of-origin and epigenomic state in pediatric H3K27M gliomas

  It is not known whether midline pediatric gliomas driven by Histone 3 K27M mutations (oncohistones) or EZHIP expression (oncohistone-mimics) share a common cell-of-origin, or arise from distinct cell types with unique vulnerabilities to PRC2 inhibition and partner alterations. Here, we define the etiological and oncogenic relationship between pediatric midline gliomas characterized by inhibition of K27M and K27M-like oncohistones. We assemble an extensive reference for gliogenesis from the developing mouse and human fetal brain. With bulk and single-cell transcriptomics and epigenomics, we profiled a large cohort of primary tumors comprising H3.1K27M and H3.3K27M pontine gliomas, H3.3K27M thalamic gliomas, and EZHIP+ posterior fossa ependymomas. We focus on differences across axes of location (thalamus vs. pons vs. posterior fossa), tumor type (diffuse midline gliomas vs. ependymomas), and oncohistone (H3.1/2K27M vs. H3.3K27M vs. EZHIP). We use tumor molecular features to delineate transcriptional states and cell-of-origin in each entity. Finally, we use culture models in isogenic contexts to interrogate the effect of each oncohistone or mimic.

  Study EGAS00001005773

• Primary prostate Hi-C

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

  Study EGAS00001005014

• Whole genome, transcriptome and methylome profiling enhances actionable target discovery in high-risk paediatric cancer

  Zero Childhood Cancer is a precision medicine program to benefit children with poor-outcome, rare, relapsed or refractory cancer. Utilising tumour and germline whole genome sequencing (WGS) and RNA sequencing (RNAseq) across 252 tumours from high-risk paediatric cancer patients, we identified 968 reportable molecular aberrations (39.9% in WGS and RNAseq, 35.1% WGS-only, 25.0% RNAseq-only). 93.7% of patients had at least one germline or somatic aberration, 71.4% had therapeutic targets and 5.2% had a change in diagnosis. WGS identified pathogenic cancer-predisposing variants in 16.2% of patients. In 76 CNS tumours, methylome analysis confirmed diagnosis in 71.1% and contributed to a change of diagnosis in two cases (2.6%). To date 43 patients have received a recommended therapy, 38 of whom could be evaluated. A small but encouraging proportion of 31% showed some evidence of clinical benefit. Comprehensive molecular profiling resolved the molecular basis of virtually all high-risk cancers, which in some cases has significant clinical impact.

  Study EGAS00001004572

• SNP array files, IDAT files, from 34 members of a Family with a high prevalence of psychosis

  We investigated the genetic causes of major mental disorders (MMDs) including schizophrenia, bipolar disorder I, major depressive disorder, and attention deficit hyperactive disorder, in a large family pedigree from Alpujarras, South of Spain, a region with a high prevalence of psychotic disorders. We applied a systematic genomic approach based on karyotyping (N=4), genotyping by genome-wide SNP array (N=33), whole-exome sequencing (N=4), and whole-genome sequencing (N=12). We performed genome-wide linkage analysis, family-based association analysis, and polygenic risk score estimates. Significant linkage was obtained at chromosome 9 (9q33.1-33.2, LOD score= 4.11), a suggestive region that contains five candidate genes ASTN2, BRINP1, C5, TLR4, and TRIM32, previously associated to MMDs. Comprehensive analysis associated the MMD phenotype with genes of the immune system with dual brain functions. Moreover, the psychotic phenotype was enriched for genes involved in synapsis. These results should be considered once studying the genetics of psychiatric disorders in other families, especially the ones from the same region, since founder effects may be related to the high prevalence.

  Study EGAS00001004592

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

  The UK’s Haematological Malignancy Research Network (www.HMRN.org) was established in 2004 to provide robust generalizable data to inform clinical practice and research. HMRN is a collaboration between researchers in the Epidemiology & Statistics Group (ECSG) at the University of York, a unified Clinical Network operating across 14 hospitals, and an integrated Haematological Malignancy Diagnostic Service (HMDS) in Leeds. Covering a population of around 4 million, HMRN collects detailed information about all patients diagnosed with a haematological malignancy within the HMRN region, accruing around 2,400 new diagnoses each year. The population has a similar socio-demographic profile to the country as a whole, and HMRN’s maturing data present an increasingly valuable resource to address real questions of concern to haematologists, commissioners and health service researchers – locally, nationally and internationally. This study forms part of a larger project employing targeted exome sequencing on tumour samples from patients diagnosed with a variety of both lymphoid and myeloid malignancies. The results for samples collected from 39 Burkitt lymphoma patients, diagnosed within HMRN from 2005 to 2012, are provided.

  Study EGAS00001004649

• atrial appendage RNA-seq in non-, paroxysmal and persistent atrial fibrillation

  Understanding the interconnectedness of structural remodelling processes atrial fibrillation (AF) in patients could identify targets for future therapies. These data include transcriptome sequencing of atrial tissues of patients without AF, with paroxysmal AF and persistent AF (total n=64). In AF patients, epicardial cell gene expression decreased, contrasting with an upregulation of epithelial to mesenchymal transition (EMT) and mesenchymal cell gene expression. Immunohistochemistry demonstrated thickening of the epicardium and an increased proportion of (myo)fibroblast-like cells in the myocardium, supporting enhanced EMT in AF. We furthermore report an upregulation of endothelial cell proliferation, angiogenesis and endothelial signalling. EMT and endothelial cell proliferation concurred with increased interstitial (myo)fibroblast-like cells and extracellular matrix gene expression including enhanced tenascin-c, thrombospondins, biglycan and versican. Morphological analysis discovered increased and redistributed glycosaminoglycans and collagens in the atria of AF patients. Signalling pathways, including cell-matrix interactions, PI3K-AKT and Notch signalling, that could regulate mesenchymal cell activation, were upregulated. These dynamic processes may provide targets to prevent or reverse structural remodelling.

  Study EGAS00001005295

• HCA_Placental_Infection_Atlas

  HCA Placental Infection Atlas Self-renewing human trophoblast stem cells (hTSC) can give rise to major trophoblast lineages, including syncytiotrophoblast (SCT) and extravillous trophoblast (EVT). hTSC has facilitated the study of placental development and function in vitro. Here, we aim to utilise the hTSC line (HMDMC 20/0005) to study human placental infection and inflammation. This study will provide the opportunity to better understand the pathogenesis of human placenta in the event of infection and inflammation. In brief, differentiated hTSC will be infected with relevant pathogens and stimulated with a panel of inflammatory cytokines. We will use single-cell genomic and spatial approaches to characterise cellular heterogeneity in response to infection or cytokines stimulus. In parallel, we will investigate the interaction between iPSC-derived macrophages with different trophoblast lineages during the course of infection using single-cell genomic and spatial methods. 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 EGAS00001004722

• Reference exome data for Australian Aboriginal populations to support health-based research

  Whole exome sequencing (WES) is a popular and successful technology which is widely used in both research and clinical settings. However, there is a paucity of reference data for Aboriginal Australians to underpin the translation of health-based genomic research. Here we provide a catalogue of variants called after sequencing the exomes of 50 Aboriginal individuals from the Northern Territory (NT) of Australia and compare these to 72 previously published exomes from a Western Australian (WA) population of Martu origin. Sequence data for both NT and WA samples were processed using an ‘intersect-then-combine’ (ITC) approach, using GATK and SAMtools to call variants. A total of 289,829 variants were identified in at least one individual in the NT cohort and 248,374 variants in at least one individual in the WA cohort. Of these, 166,719 variants were present in both cohorts, whilst 123,110 variants were private to the NT cohort and 81,655 were private to the WA cohort. Our data set provides a useful reference point for genomic studies on Aboriginal Australians.

  Study EGAS00001003745

• An essential role for MYB in driving oncogenic EVI1 expression in enhancer-rearranged leukemias

  In acute myeloid leukemia (AML) with inv(3)(q21;q26) or t(3;3)(q21;q26), a translocated GATA2 enhancer drives oncogenic expression of EVI1. We generated an EVI1-GFP AML model and applied an unbiased CRISPR/Cas9 enhancer scan to uncover sequence motifs essential for EVI1 transcription. Using this approach, we pinpoint a single regulatory element in the translocated GATA2 enhancer that is critically required for aberrant EVI1 expression while being dispensable for GATA2 expression from its endogenous locus. This element contains a DNA binding motif for the transcription factor MYB that heavily occupies this site specifically at the translocated allele. MYB knockout as well as peptidomimetic blockade of p300-dependent MYB function resulted in downregulation of EVI1 but not of GATA2. Targeting MYB or mutating its DNA-binding motif within the GATA2 enhancer resulted in myeloid differentiation and cell death, suggesting that interference with MYB-driven EVI1 transcription provides a potential entry point for therapy of inv(3)/t(3;3) AMLs.

  Study EGAS00001004839

• WGS of T-cell and NK-cell lymphoma for ICGC (NKTL-SG)

  Lymphomas are malignancies originating from the lymphatic system and there is a worldwide significant increase in incidences of this cancer. Lymphomas can be further classified into 2 categories namely Hodgkin’s (HL), and Non-Hodgkin’s lymphoma (NHL). The latter can be divided further into B, T and NK cell NHL. Peripheral T-cell lymphoma (PTCL) and Natural Killer/ T-cell lymphoma (NKTL) arise from mature T lymphocytes and NK cells. The more common subtypes include NKTL, PTCL not-otherwise-specified (PTCL-NOS), systemic anaplastic large cell lymphoma (ALCL), follicular PTCL (F-PTCL), angioimmunoblastic T-cell lymphoma (AITL) and enteropathy-associated T-cell lymphoma (EATL). Whole genome analyses on T-cell and NK-cell lymphoma samples from various countries to gain a better understanding of the pathogenesis and biology of these tumors. The comprehensive profiling and comparison among patients from various countries with different etiopathogenesis will hopefully facilitate the elucidation of key pathways and underlying molecular mechanisms driving these cancers, opening up novel therapeutic strategies in the future.

  Study EGAS00001002398

• Cell-free DNA TAPS for early cancer detection

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

  Study EGAS00001004962

• A combined circulating tumor DNA and protein biomarker-based liquid biopsy for the earlier detection of pancreatic cancer

  The earlier diagnosis of cancer is one of the keys to reducing future cancer deaths. Here we describe our efforts to develop a non-invasive, blood test for the detection of pancreatic ductal adenocarcinoma. We combined blood tests for KRAS gene mutations with carefully thresholded protein biomarkers to determine whether the combination of these markers was superior to any single marker. The cohort tested included 221 patients with resectable pancreatic ductal adenocarcinoma and 185 control patients without known cancers. KRAS mutations were detected in the plasma of 66 patients (30%) and every mutation found in the plasma was identical to that subsequently found in the patients' primary tumor (100% concordance). The use of KRAS in conjunction with two thresholded protein biomarkers (CA19-9 and prolactin) increased the sensitivity to 67%. Only one of the 185 plasma samples from the control cohort were positive for any of the three markers (99.5% specificity). This combinatorial approach may prove useful for the earlier detection of many cancer types.

  Study EGAS00001002444

• Mutations in SPINK2 induce azoospermia

  Azoospermia, characterized by the absence of spermatozoa in the ejaculate is a common cause of male infertility with a poorly characterized etiology. Exome sequencing analysis of two azoospermic brothers allowed the identification of a homozygous splice mutation in SPINK2, encoding a serine protease inhibitor believed to target acrosin, the main sperm acrosomal protease. In accord with these findings we observed that homozygous Spink2 KO male mice had azoospermia. Moreover, despite normal fertility, heterozygous male mice had a high rate of morphologically abnormal spermatozoa and a reduced sperm motility. Further analysis demonstrated that in the absence of Spink2, protease-induced stress initiates Golgi fragmentation and prevents acrosome biogenesis leading to spermatid differentiation arrest. We also observed a deleterious effect of acrosin overexpression in HEK cells, effect that was alleviated by SPINK2 coexpression confirming its role as acrosin inhibitor. These results demonstrate that SPINK2 is necessary to neutralize proteases during their cellular transit towards the acrosome and that its deficiency induces a pathological continuum ranging from oligoasthenoteratozoospermia in heterozygotes to azoospermia in homozygotes.

  Study EGAS00001002450

• Description and determinants of the faecal resistome and microbiome of farmers and slaughterhouse workers: a metagenome-wide cross-sectional study.

  Within this study, we extracted DNA of 194 stool samples of persons living or working on pig and broiler farms (54 pig farmers, 24 broiler farmers), 70 pig slaughterhouse workers and 46 control subjects (Lifelines population living in the Northern parts of the Netherlands, and without recent antimicrobial use and occupational livestock contact). Total DNA was sequenced and metagenomic analysis was performed. Subsequently, the faecal resistomes and microbiomes were described and analysed (doi: 10.1016/j.envint.2020.105939). Highlights study: • Food chain workers are exposed to antimicrobial resistant (AMR) bacteria in livestock. • We describe faecal resistomes and microbiomes of farmers and slaughterhouse workers. • Faecal AMR gene loads in pork exposed workers exceed those in broiler farmers/controls. • Resistome and microbiome differences are shown between and within occupational groups. • On-farm working hours and working/living on a pig farm determine the faecal resistome. • Contact with livestock may be a determinant for human AMR gene carriage. This study is part of the EFFORT project: http://www.effort-against-amr.eu/.

  Study EGAS00001003944

• Integrative genomic analysis identifies multiple subtypes and therapeutic targets in acute erythroid leukemia

  Acute erythroid leukemia (AEL) is a high risk leukemia of poorly understood genetic basis, and controversy regarding diagnosis in the spectrum of myelodysplasia and myeloid leukemia. We compared genomic and clinical features of 159 childhood and adult AEL cases to non-AEL myeloid disorders, and identified patterns of mutation that describe 5 age-related subgroups with distinct transcriptional profiles: adult, TP53-mutated; NPM1-mutated; KMT2A-mutated/rearranged; adult, DDX41-mutated; and pediatric, NUP98-rearranged. Mutational spectra were age-dependent and significantly associated with outcome, with NPM1 mutations and high HOXB9 expression being associated with favorable prognosis, and TP53, FLT3 and RB1 alterations associated with poor survival. Targetable signaling mutations were present in 45% of cases, and include recurrent mutations of ALK and NTRK1, the latter of which drive erythroid leukemogenesis sensitive to TRK inhibition. This genomic landscape of AEL provides the framework for accurate diagnosis and risk stratification of this disease, and the rationale for testing targeted therapies in this form of leukemia that has hitherto and dismal outcome.

  Study EGAS00001002537

• Genetic history of the Comorian populations.

  At the dawn of the second millennium, the expansion of the Indian Ocean trading network aligned with the emergence of an outward-oriented community along the East African coast to create a cosmopolitan cultural and trading zone known as the Swahili Corridor. Based on analyses of genome-wide genotyping data in 140 individuals from coastal Kenya and the Comoros islands, along with 3,477 individuals from the Indian Ocean rim, we reconstruct historical population dynamics showing that the Swahili Corridor is largely an East Bantu genetic continuum. Within this continuity significant gene flows from the Middle East can be seen in the Swahili and Comorians at dates corresponding to the Islamization of East Africa. However the main external gene flow into insular populations of the Swahili Corridor, such as Comorian and Malagasy groups, came from Island Southeast Asia as early as the 10th century. Remarkably, our results reflect archaeological and linguistic data in suggesting that the Comoros archipelago is the oldest place of contact between Austronesian and African populations in the Swahili Corridor.

  Study EGAS00001002565

• Genes associated with pancreas development and function maintain open chromatin in iPSCs generated from human pancreatic beta cells

  Current in-vitro islet differentiation protocols suffer from heterogeneity and low efficiency. Induced-pluripotent stem cells (iPSCs) derived from pancreatic beta-cells (BiPSCs) preferentially differentiate towards endocrine pancreas-like cells versus those from fibroblasts (FiPSCs). We interrogated genome-wide open chromatin in BiPSCs and FiPSCs via ATAC-seq, and identified ~8.3k significant, differential open chromatin sites (DOCS) between the two iPSC subtypes (FDR<0.05). DOCS where chromatin was more accessible in BiPSCs (Bi-DOCS) were significantly enriched for known regulators of endodermal development, including bivalent and weak enhancers, and FOXA2 binding sites (FDR<0.05). Bi-DOCS were associated with genes related to pancreas development and beta-cell function, including transcription factors mutated in monogenic diabetes (PDX1, NKX2-2, HNF1A; FDR<0.05). Moreover, Bi-DOCS correlated with enhanced gene expression in BiPSC-derived definitive endoderm and pancreatic progenitor cells. Bi-DOCS therefore highlight genes and pathways governing islet-lineage commitment, which can be exploited for differentiation protocol optimisation, diabetes disease modelling, and therapeutic purposes.

  Study EGAS00001002591

• Intratumoral plasma cells predict outcomes to PD-L1 blockade in non-small cell lung cancer

  Inhibitors of the programmed cell death-1 (PD-1/PD-L1) signaling axis are approved to treat non-small cell lung cancer (NSCLC) patients, based on their significant overall survival (OS) benefit. Using transcriptomic analysis of 891 NSCLC tumors from patients treated with either the PD-L1 inhibitor atezolizumab or chemotherapy from two large randomized clinical trials, we found a significant B cell association with extended OS with PD-L1 blockade, independent of CD8+ T cell signals. We then derived gene signatures corresponding to the dominant B cell subsets present in NSCLC from single-cell RNA-seq data. Importantly, we found increased plasma cell signatures to be predictive of OS in patients treated with atezolizumab, but not chemotherapy. B cells were also associated with the presence of tertiary lymphoid structures and organized lymphoid aggregates. Our results suggest an important contribution of B and plasma cells to PD-L1 blockade efficacy in NSCLC.

  Study EGAS00001005013

• Kalirin-RAC controls nucleokinetic migration in ADRN-type neuroblastoma

  The migrational propensity of neuroblastoma is affected by cell identity, but the mechanisms behind the divergence remain unknown. Using RNAi and time-lapse imaging, we show that ADRN-type NB cells exhibit RAC1- and kalirin-dependent nucleokinetic (NUC) migration that relies on several integral components of neuronal migration. Inhibition of NUC migration by RAC1 and kalirin-GEF1 inhibitors occurs without hampering cell proliferation and ADRN identity. Using three clinically-relevant expression dichotomies, we reveal that most of up-regulated mRNAs in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells are associated with low-risk characteristics. The computational analysis shows that, in a context of overall gene set poverty, the upregulomes in RAC1- and kalirin-GEF1-suppressed ADRN-type cells are a batch of AU-rich elements (ARE)-containing mRNAs, which suggests a link between NUC migration and mRNA stability. GSEA-based search for vulnerabilities reveals prospective weak points in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells, including activities of H3K27- and DNA methyltransferases. Altogether, these data support the introduction of NUC inhibitors into cancer treatment research.

  Study EGAS00001005023

• HCA_Female_Reporductive_Adult_WSSS_RNA

  Cell atlas of the human placenta 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. This specific project will assess human development. We aim to interrogate the cellular composition and molecular regulators underpinning development and maturation through single cell RNA-sequencing/spatial transcriptomics and computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequent validate these findings in situ using allied imaging technologies e.g. immunohistochemistry and RNAscope and in vitro culture of human intestinal organoids (HIOs). 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 EGAS00001004210

• Genomic and functional fidelity of small cell lung cancer patient-derived xenografts

  Small cell lung cancer (SCLC) patient-derived xenografts (PDXs) can be generated from biopsies or circulating tumor cells (CTCs), though scarcity of tissue and low efficiency of tumor growth have previously limited these approaches. Applying an established clinical-translational pipeline for tissue collection and an automated microfluidic platform for CTC-enrichment, we generated 17 biopsy-derived PDXs and 17 CTC-derived PDXs in a two-year timeframe, at efficiencies of 89% and 38%, respectively. Whole exome sequencing showed that somatic alterations are stably maintained between patient tumors and PDXs. Early-passage PDXs maintain the genomic and transcriptional profiles of the founder PDX. In vivo treatment with etoposide and cisplatin (EP) in 30 PDX models demonstrated greater sensitivity in PDXs from EP naïve patients, and resistance to EP corresponded to increased expression of a MYC gene signature. Finally, serial CTC-derived PDXs generated from an individual patient at multiple time points accurately recapitulated the evolving drug sensitivities of that patient’s disease. Collectively, this work highlights the translational potential of this strategy.

  Study EGAS00001002853

• Exome sequencing demonstrates a dual origin of relapses in retinoic-acid resistant acute promyelocytic leukemia.

  Retinoic acid (RA) and arsenic target the PML/RARA driver of Acute Promyelocytic Leukemia (APL), their combination now curing over 95% patients. Here we report exome sequencing of 64 matched diagnosis/remission/relapses patients treated with historical RA/chemotherapy protocols. A first subgroup presents a high incidence of oncogenic mutations disrupting key epigenetic or transcriptional regulators (primarily WT1) or activating MAPK signaling at diagnosis. Relapses retain these cooperating oncogenes and exhibit additional oncogenic alterations and/or mutations impeding therapy response (RARA, NT5C2). The second group primarily exhibits FLT3 activation at diagnosis, which is lost upon relapse together with most other passenger mutations, implying that relapses derive from pre-leukemic PML/RARA-expressing cells that survived RA/chemotherapy. Accordingly, clonogenic activity of PML/RARA-immortalized progenitors ex vivo is only transiently affected by RA, but selectively abrogated by arsenic. Our studies stress the role of cooperating oncogenes in direct relapses and suggest that targeting pre-leukemic cells by arsenic contributes to its clinical efficacy.

  Study EGAS00001002893

• Somatic_mutation_profiling_of_intestinal_crypts_from_IBD

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

  Study EGAS00001002896

• Exome Sequecning of MDS xenografted samples

  Pre-clinical research of Myelodysplastic Syndromes (MDS) is hampered by a lack of feasible disease models. Previously, we have established a robust patient-derived xenograft (PDX) model for MDS. Here, we demonstrate for the first time that this model is applicable as a pre-clinical platform to address pending clinical questions by interrogating the efficacy and safety of the thrombopoietin receptor agonist eltrombopag. Our pre-clinical study included n=49 xenografts generated from n=9 MDS patient samples. Substance efficacy was evidenced by FACS-based human platelet quantification, and clonal bone marrow evolution was reconstructed by serial whole exome sequencing of the PDX samples. In contrast to clinical trials in humans, this experimental setup allowed vehicle- and replicate-controlled analyses on a patient-individual level deciphering substance-specific effects from natural disease progression. We found that eltrombopag effectively stimulated thrombopoiesis in MDS PDX without adversely affecting the patients’ clonal composition. In conclusion, our MDS PDX model is a useful tool for testing new therapeutic concepts in MDS preceding clinical trials.

  Study EGAS00001005329

• An instructive role for IL7RA in the development of human B-cell precursor leukemia

  Kinase signaling fuels growth of B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Yet its role in leukemia initiation is unclear and has not been shown in primary human hematopoietic cells. We previously described activating mutations in interleukin-7 receptor alpha (IL7RA) in poor-prognosis “ph-like” BCP-ALL. Here we show that expression of activated mutant IL7RA in human CD34+ hematopoietic cells creates a preleukemic state in transplanted NOD/LtSz-scid IL2Rγnull mice, characterized by persistence of self-renewing Pro-B cells with non-productive V(D)J gene rearrangements. Preleukemic CD34+CD10highCD19+ cells evolved to a BCP-ALL with spontaneously acquired CDKN2A deletions as commonly observed in primary human BCP-ALL. CRISPR mediated silencing of CDKN2A in primary human CD34+ cells transduced with activated IL7RA resulted in robust development of BCP-ALLs in-vivo. Thus, we demonstrate for the first time that constitutive activation of IL7RA can initiate preleukemia in primary human hematopoietic progenitors and cooperates with CDKN2A silencing in progression into BCP-ALL.

  Study EGAS00001005347

• Mapping the breast cancer metastatic cascade onto circulating tumour DNA using genetic and epigenetic clonal tracking

  Circulating tumour DNA (ctDNA) allows tracking the evolution of human cancers at high resolution, overcoming many limitations of tissue biopsies. However, exploiting ctDNA to determine how a patient’s cancer is evolving in order to aid clinical decisions remains difficult. This is because ctDNA is a mix of fragmented alleles, and the contribution of different cancer cell subpopulations from distinct lesions to ctDNA remains largely unknown. Profiling of ctDNA almost invariably requires prior knowledge of what genomic alterations to track. Here we leverage on a rapid autopsy programme to demonstrate that unbiased genomic characterisation of several metastatic sites and concomitant ctDNA profiling at whole-genome resolution reveals the extent to which ctDNA is representative of widespread disease. We also present a ctDNA methylation profiling method that allows tracking evolutionary changes in the cancer cell population at single-molecule resolution without prior knowledge. These results have critical implications for the use of liquid biopsies to monitor cancer evolution in humans and potentially guide treatment.

  Study EGAS00001004014

• HCA_Placenta_Adult_Vento_RNA

  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. This specific project will assess human development. We aim to interrogate the cellular composition and molecular regulators underpinning development and maturation through single cell RNA-sequencing/spatial transcriptomics and computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequent validate these findings in situ using allied imaging technologies e.g. immunohistochemistry and RNAscope and in vitro culture of human intestinal organoids (HIOs). 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 EGAS00001004069

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

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

  Study EGAS00001005407

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

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

  Study EGAS00001005408

• Comprehensive pharmacogenomic characterization of gastric cancer

  Gastric cancer is among the most lethal human malignancies. Previous studies have identified molecular aberrations that constitute dynamic biological networks and genomic complexities of gastric tumors. However, the clinical translation of molecular-guided targeted therapy is hampered by challenges. Notably, solid tumors often harbor multiple genetic alterations, complicating the development of effective treatments. To address such challenges, we established a comprehensive dataset of molecularly annotated patient derivatives coupled with pharmacological profiles for 60 targeted agents to explore dynamic pharmacogenomic interactions in gastric cancers. We identified lineage-specific drug sensitivities based on histopathological and molecular subclassification, including substantial sensitivities toward VEGFR and EGFR inhibition therapies in diffuse- and signet ring-type gastric tumors, respectively. We identified potential therapeutic opportunities for WNT pathway inhibitors in ALK-mutant tumors, a significant association between PIK3CA-E542K mutation and AZD5363 response, and transcriptome expression of RNF11 as a potential predictor of response to gefitinib. Collectively, our results demonstrate the feasibility of drug screening combined with tumor molecular characterization to facilitate personalized therapeutic regimens for gastric tumors.

  Study EGAS00001004106

• Cancer immune control needs senescence induction by Stat1 dependent cell cycle regulator pathways in tumours

  Immune checkpoint bloackade (ICB)-based or natural cancer immune responses largely eliminate tumours. Yet, they require additional mechanisms to arrest those cancer cells that are not rejected. Cytokine-induced senescence (CIS) can stably arrest cancer cells, suggesting that interferon-dependent induction of senescence-inducing cell cycle regulators is needed to control those cancer cells that escape from killing. Here we report in two different cancers sensitive to T cell-mediated rejection, we show that deletion of the senescence-inducing cell cycle regulators p16Ink4a/p19Arf (Cdkn2a) or p21Cip1 (Cdkn1a) in the tumour cells abrogated both, the natural and the ICB-induced cancer immune control. Also in humans, melanoma metastases that progressed rapidly during ICB have losses of senescence-inducing genes and amplifications of senescence inhibitors. Metastatic cells also resist CIS. Such genetic and functional alterations are infrequent in metastatic melanomas regressing during ICB. Thus, activation of tumour-intrinsic, senescence-inducing cell cycle regulators is required to stably arrest those cancer cells that escape from eradication.

  Study EGAS00001004151

• A blood atlas of COVID-19 defines hallmarks of disease severity and specificity

  Treatment of severe COVID-19 is currently limited by clinical heterogeneity and incomplete understanding of potentially druggable immune mediators of disease. To advance this, we present a comprehensive multi-omic blood atlas in patients with varying COVID-19 severity and compare with influenza, sepsis and healthy volunteers. We identify immune signatures and correlates of host response. Hallmarks of disease severity revealed cells, their inflammatory mediators and networks as potential therapeutic targets, including progenitor cells and specific myeloid and lymphocyte subsets, features of the immune repertoire, acute phase response, metabolism and coagulation. Persisting immune activation involving AP-1/p38MAPK was a specific feature of COVID-19. The plasma proteome enabled sub-phenotyping into patient clusters, predictive of severity and outcome. Tensor and matrix decomposition of the overall dataset revealed feature groupings linked with disease severity and specificity. Our systems-based integrative approach and blood atlas will inform future drug development, clinical trial design and personalised medicine approaches for COVID-19.

  Study EGAS00001005493

• Automated image-based profiling of complex drug induced phenotypes in patient-derived organoids

  We established a platform for high-throughput confocal fluorescence microscopy and automated image analysis to perform large-scale drug response profiling with patient derived organoids (PDOs). Using PDO lines newly established from colorectal cancer patients with a standardized workflow, we performed chemical perturbations with more than 500 experimental and clinically used compounds. Subsequently, we captured approx. six million images by confocal microscopy. To analyze drug-induced PDO re-organization on a single organoid level we developed SCOPE, a software framework for automated PDO high-throughput image analysis. We observed a rich variety of divergent and reoccurring compound induced phenotypes in organoids. Perturbation of cellular processes, including signaling by MEK, GSK3β, cyclin dependent kinases (CDKs) and others, led to recurring and distinct architectural changes. Highlighting the clinical relevance, PDO viability was heterogeneous after perturbation with anticancer drugs and matched clinical response of corresponding colorectal cancer patients. Our work opens new avenues for image-based profiling for accelerated drug discovery and clinical decision-making.

  Study EGAS00001003140

• Pancreatic islets PISA RNA-seq samples

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

  Study EGAS00001005535

• Molecular origins of mpMRI visibility

  Multiparametric magnetic resonance imaging (mpMRI) has transformed the landscape of prostate cancer (PCa) management. Prostate Imaging Reporting and Data System (PI-RADSv2) scores reflect function and correlate phenotypic features with tumor aggressiveness, is an independent predictor of biochemical recurrence, and may distinguish low-risk tumors requiring no treatment from high-risk tumors requiring aggressive management. We hypothesize that differences in functional mpMRI of prostate tumors reflect the genetic and epigenetic heterogeneity of PCa. We profiled the genomes and transcriptomes of 40 Gleason score 3+4 tumors, of which 20 tumors were mpMRI invisible (not PI-RADSv2 3-5) and 20 tumors were visible (PI-RADSv2 5). Genomes of visible tumors were more unstable than those of invisible tumors. We identified 62 enriched transcripts, including numerous snoRNAs, Relaxin 1 and SCHLAP1, in visible tumors. These data suggest a confluence of aggressive pathological and microenvironmental phenomenon in PI-RADSv2 5 tumors reminiscent of nimbosus, while linking the molecular profiles of tumor aggressiveness with their functional features on mpMRI.

  Study EGAS00001003179

• MethCORR: DNA Methylation-based Characterization, Classification and Prognostication of Colorectal Cancer using Archival Formalin-fixed, Paraffin-embedded Tissue

  Transcriptional characterization and classification has potential to resolve the inter-tumor heterogeneity of colorectal cancer (CRC) and improve patient management. Yet, robust transcriptional profiling is difficult using formalin-fixed, paraffin-embedded (FFPE) samples, which complicates testing in clinical and archival material. We present MethCORR, an approach that allows uniform molecular characterization and classification of fresh-frozen and FFPE samples. MethCORR identifies genome-wide correlations between RNA expression and DNA methylation in fresh-frozen samples. This information is used to infer gene expression information in FFPE samples from their DNA methylation profiles. MethCORR was applied to methylation profiles from 877 fresh-frozen and FFPE samples and comparative analysis identified the same two major subtypes in four independent cohorts. Furthermore, subtype-specific prognostic biomarkers that better predicted relapse-free survival (HR=2.66, 95% CI [1.67-4.22], P<0.001) than TNM staging and MSI status were identified and validated using DNA methylation-specific PCR assays. The MethCORR approach is general, and may be similarly successful for other cancer types.

  Study EGAS00001004293

• Single-nuclei gene-expression analysis of pheochromocytoma and paraganglioma links tumor subtypes with tumor microenvironment

  Pheochromocytomas and paragangliomas (PCPG) are rare neuroendocrine tumors associated with autonomic nerves. We used single nuclei-RNA-seq (snRNA-seq) for analysis of 30 PCPG representing 13 known driver genes, plus two normal adrenal medullas to dissect cell composition, refine PCPG subtypes and compare PCPG and normal tissue expression. Incorporating bulk-tissue and snRNA-seq data we identified seven PCPG gene-expression subtypes with genotype and cell type associations. The tumor microenvironment reflected HIF pathway activation and neovascularization by detection of abundant tip-like endothelial cells and pericytes concordant with VEGFA expression in neoplastic cells with mutations in VHL, genes encoding SDH subunits and MAML3-fusions. Macrophages were abundant in PCPG with distinct tumor-associated features. Schwann-cell-like cells were also enriched in some subtypes. The transcriptional profile of neoplastic PCPG cells resembled normal chromaffin cells but early chromaffin and neuroblast cell markers were expressed in some subtypes including GPR139, a putative therapeutic target for a subset of pseudohypoxic PCPG.

  Study EGAS00001005861

• Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosissusceptibility

  Multiple sclerosis (MS) is characterized by a targeted attack on oligodendroglia (OLG)and associated myelin by immune cells, which are thought to be the main drivers of MS susceptibility. Assessing chromatin accessibility and the transcriptome simultaneously at the single cell level, we found that immune genes exhibit a primedchromatin state in mouse and human OLG in a non-disease context, compatible withrapid transitions to immune-competent states in MS. We identified transcription factors as BACH1 and STAT1 involved in immune gene regulation in oligodendrocyteprecursor cells (OPCs). A subset of immune genes present bivalency of H3K4me3/H3K27me3 in OPCs, with Polycomb inhibition leading to their increased activation upon interferon-gamma (IFNg) treatment. Some MS susceptibility single-nucleotide polymorphisms (SNPs) overlap with these regulatory regions in mouse andhuman OLG, and treatment of mouse OPCs with IFNg leads to chromatin architecture remodeling at these loci and altered expression of interacting genes. Our data indicates that susceptibility for MS may involve OLG, which therefore constitute novel targets for immunological-based therapies for MS.

  Study EGAS00001005911

• Molecular stratification of endometrioid ovarian carcinoma predicts clinical outcome

  Endometrioid ovarian carcinoma (EnOC) is an under-investigated type of ovarian carcinoma. Here, we report the largest genomic study of EnOCs to date, performing whole exome sequencing of 112 cases following rigorous pathological assessment. High frequencies of mutation were detected in CTNNB1(43%), PIK3CA(43%), ARID1A(36%), PTEN(29%), TP53(26%) and SOX8(19%), a novel target of recurrent mutation in EnOC. POLE and mismatch repair protein-encoding genes were mutated at lower frequency (6%, 18%) with significant co-occurrence. A molecular taxonomy was constructed using a novel algorithm (PRISTINE), identifying clinically distinct EnOC subtypes: TP53m cases demonstrated greater genomic complexity, were frequently FIGO stage III/IV at diagnosis (48%) and incompletely debulked (44%), and demonstrated inferior survival; conversely, CTNNB1m cases demonstrated low complexity and excellent clinical outcome, were predominantly stage I/II at diagnosis (89%) and completely resected (87%). Tumour complexity provides further resolution within the TP53wt/CTNNB1wt group. Moreover, we identify the WNT, MAPK/RAS and PI3K pathways as good candidate targets for molecular therapeutics in EnOC.

  Study EGAS00001004366

• Sensitive detection of tumor mutations from blood and its application to immunotherapy prognosis

  Cell-free DNA (cfDNA) is attractive for many applications, including cancer detection, locating, and monitoring. A fundamental task underlying these applications is the SNV calling from cfDNA, which, however, faces a new challenge, namely, the generally very low tumor content in cfDNA. Thus all existing callers fail to achieve satisfactory performance. Here we present cfSNV, a method incorporating multi-layer error suppression and hierarchical mutation calling, to address this important challenge. Furthermore, by leveraging cfDNA’s comprehensive coverage of clonal landscape, for the first time cfSNV can profile mutations even in subclones. In both simulated and real patient data, cfSNV vastly outperforms existing tools, showing tens of times increase in sensitivity in detecting mutations with low allele frequency while maintaining high precision. cfSNV can enhance the clinical utilities of cfDNA by dramatically reducing the required sequencing depth and therefore reduce the cost by magnitudes, and further make the Whole-Exome-Sequencing of cfDNA a viable option. As an example, we demonstrate that cfDNA-WES allows a new biomarker to effectively select patients for immunotherapy.

  Study EGAS00001004373

• Tumor-intrinsic expression of the autophagy gene Atg16l1 suppresses anti-tumor immunity in colorectal cancer

  The vast majority of colorectal cancer (CRC) patients fail to respond to immune checkpoint inhibitors (ICI). Understanding tumor-intrinsic determinants of immunotherapy resistance is critical to improve CRC patient outcomes. Here, we demonstrate that transcript levels of the core autophagy gene ATG16L1 are prognostic of poor outcome and predictive of ICI resistance in late-stage CRC. Deletion of Atg16l1 in murine CRC organoids markedly inhibited tumor growth in primary and metastatic niches in syngeneic hosts. ATG16L1 deficiency rendered CRC organoids more responsive to interferon (IFN) signaling and enhanced programmed cell death while decreasing tumor stem-like populations. IFN hypersensitivity consequently promoted cytotoxic anti-tumor immunity in vivo. In IMblaze370, a large phase III clinical trial of anti-PD-L1 in advanced metastatic CRC, low ATG16L1 levels identified patients with increased lymphocyte infiltration and improved outcome in a KRAS mutant setting. This work reveals autophagy as a clinically relevant immunosuppressive mechanism in CRC and provides a rationale for autophagy inhibition to boost immunotherapy responses in the clinic.

  Study EGAS00001005952

• Infiltrative and drug-resistant slow-cycling cells support metabolic heterogeneity in glioblastoma

  Metabolic reprogramming has been described in rapidly growing tumors, which are thought to mostly contain fast-cycling cells (FCCs) that have impaired mitochondrial function and rely on aerobic glycolysis. Here, we characterize the metabolic landscape of glioblastoma (GBM) and explore metabolic specificities as targetable vulnerabilities. Our studies highlight the metabolic heterogeneity in GBM, in which FCCs harness aerobic glycolysis, and slow-cycling cells (SCCs) preferentially utilize mitochondrial oxidative phosphorylation for their functions. SCCs display enhanced invasion and chemoresistance, suggesting their important role in tumor recurrence. SCCs also demonstrate increased lipid contents that are specifically metabolized under glucose-deprived conditions. Increased fatty acid transport in SCCs is targetable by pharmacological inhibition or genomic deletion of FABP7, both of which sensitize SCCs to metabolic stress. Furthermore, FABP7 inhibition, whether alone or in combination with glycolysis inhibition, leads to overall increased survival. Our studies reveal the existence of GBM cell subpopulations with distinct metabolic requirements and suggest that FABP7 is central to lipid metabolism in SCCs and that targeting FABP7-related metabolic pathways is a viable therapeutic strategy.

  Study EGAS00001003251

• Single-cell analysis of airway samples identifies immune cell activation correlating with COVID-19 disease severity

  To investigate the immune response and mechanisms associated with severe COVID-19, we performed single-cell RNA-seq on nasopharyngeal and bronchial samples from 19 clinically well-characterized patients with moderate or critical disease and from 5 healthy controls. We identified airway epithelial cell types and states vulnerable to SARS-CoV-2 infection. In COVID-19 patients, epithelial cells showed an average threefold increase in expression of the SARS-CoV-2 entry receptor ACE2, which correlated with interferon signals by immune cells. Compared with moderate cases, critical cases exhibited stronger interactions between epithelial and immune cells, as indicated by ligand–receptor expression profiles, and activated immune cells , including inflammatory macrophages expressing CCL2, CCL3, CCL20, CXCL1, CXCL3, CXCL10, IL8, IL1B and TNF . The transcriptional differences in critical cases compared with moderate cases likely contribute to clinical observations of heightened inflammatory tissue damage, lung injury and respiratory failure. Our data suggest that pharmacologic inhibition of the CCR1 and/or CCR5 pathways may suppress immune hyperactivation in critical COVID-19.

  Study EGAS00001004481

• Three-dimensional human alveolar stem cell culture models reveal infection response to SARS-CoV-2

  Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), which is the cause of a present global pandemic, infects human lung alveolar type 2 (hAT2) cells. Characterizing pathogenesis is crucial for developing vaccines and therapeutics. However, the lack of models mirroring the cellular physiology and pathology of hAT2 cells limits the study. Here, we develop a feeder-free, long-term three-dimensional (3D) culture technique for hAT2 cells derived from primary human lung tissue, and investigate infection response to SARS-CoV-2. By imaging-based analysis and single-cell transcriptome profiling, we reveal rapid viral replication and the increased expression of interferon-associated genes and pro-inflammatory genes in infected hAT2 cells, indicating robust endogenous innate immune response. Further tracing of viral mutations acquired during transmission identifies full infection of individual cells effectively from a single viral entry. Our study provides deep insights into the pathogenesis of SARS-CoV-2, and the application of defined 3D hAT2 cultures as models for respiratory diseases.

  Study EGAS00001004508

• Whole-exome sequencing of glioblastomas with long-term relapse interval

  Treatment options for GBM patients have not changed over the last 15 years, and the success of targeted therapy or immunotherapy has been disappointinglimited. Thus, the standard of care for GBM patients remains surgical resection followed by concomitant DNA alkylating agent therapy with temozolomide (TMZ) and radiation therapy (RT). However, the time to relapse is generally short (6.9 months on average) owing to the tumor’s ability to infiltrate normal parenchyma and its intrinsic resistance to radio-chemotherapy. Due to the lack of alternative treatment options, elucidating underlying resistance mechanism has come into major focus of research, with the hope for opening new ways of therapeutic intervention overcoming resistance. To explore therapy resistance driving events, we performed a study using a single-centric collective of matched pairs of therapy-naïve and post-therapeuticc tissue from GBM patients, selected based on a maximally long interval between resections of the initial and the recurrent tumors. This rare collective was analyzed by whole exome sequencing (WES), mRNA expression and miRNA expression profiling.

  Study EGAS00001006022

• Functional screening on patient-derived organoids identifies a therapeutic bispecific antibody that triggers EGFR degradation in LGR5+ tumor cells

  Patient-derived organoids (PDOs) have demonstrated predictive value in prospective clinical trials supporting selection of personalized treatments. Because PDOs retain the organization and physiological functions of their source tissue, PDO biobanks could also be an ideal substrate to screen for novel therapeutic interventions. Here we describe a large-scale functional screen of dual targeting bispecific antibodies (bAbs) on a colorectal cancer (CRC) PDO biobank to target their dependency on cancer stem cells. A novel drug discovery pipeline was assembled where therapeutic bAb panels generated against WNT and receptor tyrosine kinases (RTK) targets were functionally evaluated by high content imaging to capture the complexity of PDO responses across a wide range of different CRCs and paired normal colonic mucosa samples. Our strategy resulted in a bAb that specifically triggers EGFR degradation in LGR5+ CSCs and that exhibits unique therapeutic properties such as potent growth inhibition of KRAS mutant CRCs, blockade of metastasis initiation in preclinical models and minimal toxicity towards normal LGR5+ colon stem cells.

  Study EGAS00001004584

• Single-cell RNA-seq of celiac disease-specific plasma cells

  Disease-specific plasma cells (PCs) reactive with transglutaminase 2 (TG2) or deamidated gluten peptides (DGP) are abundant in celiac disease (CeD) gut lesions. Their contribution toward CeD pathogenesis is unclear. We assessed expression of markers associated with PC longevity in 15 untreated and 26 treated CeD patients in addition to 13 non-CeD controls, and performed RNA-sequencing with clonal inference and transcriptomic analysis of 3251 single PCs. We observed antigen-dependent V-gene selection and stereotypic antibodies. Generation of recombinant DGP-specific antibodies revealed a key role of a heavy-chain residue that displays polymorphism, suggesting that immunoglobulin gene polymorphisms may influence CeD-specific antibody responses. We identified transcriptional differences between CeD-specific vs non-disease-specific PCs and between short-lived vs long-lived PCs. The short-lived CD19+CD45+ phenotype dominated in untreated and short-term-treated CeD, in particular among disease-specific PCs but also in the general PC population. Thus, the disease lesion of untreated CeD is characterized by massive accumulation of short-lived PCs that are not only directed against disease-specific antigens.

  Study EGAS00001004623

• Tumour gene expression signature in primary melanoma predicts long-term outcomes: A prospective multicentre study

  Adjuvant systemic therapies are now routinely used following resection of stage III melanoma, however accurate prognostic information is needed to better stratify patients. We used differential expression analyses of primary tumours from 204 RNA-sequenced melanomas within a large adjuvant trial, identifying a 121 metastasis-associated gene signature. This signature strongly associated with progression-free (HR=1.63, p=5.24x10-5) and overall survival (HR=1.61, p=1.67x10-4), and validated in 177 regional lymph nodes metastasis as well as two externally ascertained datasets. The machine learning classification models trained using the signature genes performed significantly better in predicting metastases than models trained with clinical covariates (pAUROC=0.02), or published prognostic signatures. The signature score negatively correlated with measures of immune cell infiltration (ρ=-0.75, p<2.2x10-16), with a higher score representing reduced lymphocyte infiltration and a higher 5-year risk of death in stage II melanoma. Our expression signature identifies melanoma patients at higher risk of metastases and warrants further evaluation in adjuvant clinical trials.

  Study EGAS00001004664

• Functional and genomic heterogeneity of long-term self-renewing compartment as the origin of treatment resistance in pancreatic tumors

  Adaptive drug-resistance mechanisms allow human tumors to evade treatment through selection and expansion of treatment-resistant clones. Here, studying clonal evolution of tumor cells derived from human pancreatic tumors, we demonstrate that in vitro cultures and in vivo tumors are maintained by a common set of tumorigenic cells that can be used to establish clonal replica tumors (CRTs), large cohorts of animals bearing human tumors with identical clonal composition. Using CRTs to conduct quantitative assessments of adaptive responses to therapeutics, we uncovered a multitude of functionally heterogeneous subpopulations of cells with differential degrees of drug sensitivity. High-throughput isolation and deep characterization of unique clonal lineages showed genetic and transcriptomic diversity underlying functionally diverse subpopulations. Molecular annotation of gemcitabine-naïve clonal lineages with distinct responses to treatment in the context of CRTs generated signatures that can predict the response to chemotherapy, representing a potential biomarker to stratify patients with pancreatic cancer. RNA-Seq data on these clones are available on EBI ENA, under accession number: PRJEB27735.

  Study EGAS00001003442

• RNA-seq of M-CSF differentiated human peripheral monocyte-derived macrophages (MDMs) - Validation macIDR

  Macrophages represent multifunctional leukocytes defined by their stimulus-specific transcriptional reprogramming. As in vivo macrophages are often difficult to obtain, in vitro macrophage models are often used. We aggregated public expression data to define consensus expression profiles for eight commonly-used in vitro macrophage models and built the classifier macIDR, capable of distinguishing macrophage subsets with high accuracy (>0.95). Classification of in vivo macrophages suggested that alveolar macrophages resembled interleukin-10 activated macrophages in general whereas chronic obstructive pulmonary disease patients displayed decreased similarity to interferon-γ stimulated macrophages. Adipose tissue-derived macrophages were classified as unstimulated macrophages, but would resemble LPS-stimulated macrophages more in diabetic-obese patients. Rheumatoid arthritic synovial macrophages were similar to macrophages stimulated with interleukin-10 or interferon-γ. Altogether, our results suggest that macIDR is capable of identifying in vitro macrophages. By projecting in vivo macrophages onto the in vitro macrophages, we were capable of elucidating macrophage-specific changes as a result of tissue and disease.

  Study EGAS00001003451

• An Isogenic Cell Line Panel for Sequence-based Screening of Targeted Anti-cancer Drugs

  We describe the creation and characterization of an isogenic cell line panel representing common cancer pathways, with multiple features optimized for high-throughput screening. More than 1,800 cell lines from three normal human cells were generated using CRISPR-technologies. Surprisingly, we discovered most of these lines did not result in complete gene inactivation, despite integration of sgRNA at the desired genomic site. However, a subset of the lines harbored true, biallelic disruptions of the targeted tumor suppressor gene, yielding a final panel of 100 well-characterize lines covering 19 pathways frequently subject to loss of function in cancers. This panel included genetic markers optimized for sequence-based ratiometric assays for drug-based screening assays. To illustrate the potential utility of this panel, we developed a multiplexed high-throughput screen that identified Wee1 inhibitor MK-1775 as a selective growth inhibitor of cells with inactivation of TP53. These cell lines and screening approach should prove useful for researchers studying a variety of cellular and biochemical phenomena

  Study EGAS00001005974

• Genetic and epigenetic variation at regulatory regions contribute to cancer evolution under endocrine treatment

  Comprehensive profiling of hormone-dependent breast cancer (HDBC) has identified hundreds of protein-coding alterations contributing to cancer initiation, but only a handful have been linked to endocrine therapy resistance, potentially contributing to 40% of relapses. If other mechanisms underlie the evolution of HDBC under adjuvant therapy is currently unknown. In this work, we employ integrative functional genomics to dissect the contribution of cis-regulatory elements (CREs) to cancer evolution by focusing on 12 megabases of non-coding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Massive parallel perturbation in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance9. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies non-coding changes involved in therapy resistance. Overall, our data uncover actionable transient transcriptional programs critical for dormant persister cells and unveil new regulatory nodes driving evolutionary trajectories towards disease progression.

  Study EGAS00001006340

• Understanding_Self__Organising_Capacity_of_Stem_Cells_during_Implantation_and_Early_Post_implantation_Development_in_vitro_and_in_vivo__Implications_for_Human_Development_

  The aim of this project is to differentiate human embryonic stem cells to an extra-embryonic fate, specifically the hypoblast. This is of uttermost importance given the current lack of human hypoblast stem cells. We hypothesized that the pluripotent characteristics of the starting human embryonic stem cell population may dictate the competency for extra-embryonic cell fate specification. Based on this hypothesis and using human embryonic stem cells maintained in different naïve-like culture regimes, we have now developed conditions that allow the differentiation of human embryonic stem cells to a stable GATA6+ SOX2- population. This suggests that these cells may be putative human hypoblast stem cells. To validate this finding here we propose to perform RNA sequencing experiments of the differentiated human embryonic stem cells. By comparing their RNA expression profile to the single cell sequencing data of the human embryo that we are currently generating, we will be able to determine the identity of our GATA6+ SOX2- cells, and establish whether they represent the in vivo human hypoblast.

  Study EGAS00001003571

• Proteogenomics reveals two distinct biological pilocytic astrocytoma subgroups

  Pilocytic astrocytoma (PA) is the most common pediatric brain tumor and driven by aberrant MAPK signaling, typically mediated by BRAF alterations. While five-year overall survival rates exceed 95%, tumor recurrence constitutes a major clinical challenge in incompletely resected tumors despite chemotherapeutic or radiation based therapies. Therefore, we used proteogenomics to discern the biological heterogeneity of PA to improve classification of this tumor entity and identify novel therapeutic targets. Our proteogenomics approach integrates RNA sequencing and LC/MS-based proteomic profiling data from a cohort of 58 confirmed, primary PA samples. An integrative genomics approach was conducted to discern the biological heterogeneity of PA and to identify aberrant pathway activation in these biological subgroups. In summary, Pilocytic astrocytomas segregate into two groups where younger patients are significantly associated with Group 1. Importantly, we validate the two distinct biological subgroups in two non-overlapping cohorts. The biological heterogeneity seen here may improve biological classification and reveal novel therapeutic targets specifically useful for non-resectable tumors with high risk of recurrent or progressive disease.

  Study EGAS00001006402

• Genomic analysis and evolutionary modeling of breast and larynx cancer, based on specimens from Polish population

  We explore the genomic characteristics of the original dataset of breast and larynx cancer specimens of several cell types from the collection of the Department of Applied Radiology of the Maria Sklodowska-Curie National Research Institute of Oncology, Krakow (Poland) Branch and Institute of Human Genetics Polish Academy of Sciences. DNA from paired tissue samples from primary tumor locations and concurrent metastasis to regional lymph nodes was sequenced and matched to controls. A range of genomic analyses was carried out, such as estimation of variant allele frequencies (VAF) and copy number variation (CNV). They show that subsets of variants from primary tumors are present in metastatic lymph nodes and that the CNV pattern is similar, but not identical to that known in the literature. In addition to that, we fit the VAF data to theoretical site frequency spectra predicted by multiple models found in the literature. Analysis suggests existence of secondary advantageous subclones arising on the clonal background as well as of a number of neutral or slightly deleterious passenger mutations.

  Study EGAS00001006456