22929 results for "ega"

in 68.77 milliseconds.

• 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

• Multi-focal genomic dissection of synchronous primary and metastatic tissue from de novo metastatic prostate cancer

  De novo metastatic prostate cancer is highly aggressive, but the paucity of routinely-collected tissue in this setting has hindered genomic stratification and jeopardizes precision oncology efforts. Here, we leveraged a rare study of surgical intervention in 43 de novo metastatic prostate cancers to assess somatic genotype across 607 synchronous primary and metastatic tissue foci plus circulating tumor DNA. Intra-prostate heterogeneity was pervasive and impacted clinically-relevant genes, resulting in frequent discordance between select primary-restricted foci and synchronous metastases. Additional complexity was driven in several instances by polyclonal metastatic seeding from the reservoir of phylogenetically-related primary populations. When simulating standard clinical practice relying on a single tissue focus, genomic heterogeneity plus highly variable per-sample tumor fraction caused false genotyping of dominant disease. However, in silico modeling demonstrated that analysis of multiple biopsy cores can rescue misassigned somatic genotypes. Our results define the relationship between synchronous treatment-naive primary and metastatic lesions in de novo metastatic prostate cancer and provide a framework for implementing genomics-guided patient management.

  Study EGAS00001006466

• Upper cortical layer-driven network impairment in schizophrenia

  Schizophrenia is one of the most wide-spread and complex mental disorders. To characterize the impact of schizophrenia, we performed single nucleus RNA sequencing (snRNA-seq) of >220,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls. Additionally, >115,000 neurons were analyzed topographically by immunohistochemistry. Compositional analysis of snRNA-seq data revealed a reduction in abundance of GABAergic neurons and a concomitant increase in principal neurons, most pronounced for upper cortical layer subtypes, which was substantiated by histological analysis. Many neuronal subtypes showed extensive transcriptomic changes, the most dramatic - in upper layer GABAergic neurons, including downregulation in energy metabolism and upregulation in neurotransmission. Transcription factor network analysis demonstrated a developmental origin of transcriptomic changes. Finally, Visium spatial transcriptomics further corroborated upper layer neuron vulnerability in schizophrenia. Overall, our results point towards general network impairment within upper cortical layers as a core substrate associated with schizophrenia symptomatology. Study is available on bioRxiv (https://doi.org/10.1101/2020.11.17.386458) and in upcoming publication.

  Study EGAS00001006495

• The Precision Medicine in Liver Cancer across an Asia-pacific NETwork (PLANET) - Intratumoural immune heterogeneity

  The clinical relevance of immune landscape intratumoural heterogeneity (immune-ITH) and its role in tumour evolution remain largely unexplored. Here, we uncover significant spatial and phenotypic immune–ITH from multiple tumour sectors and decipher its relationship with tumour evolution and disease progression in hepatocellular carcinomas (HCC). Immune–ITH is associated with tumour transcriptomic-ITH, mutational burden, and distinct immune microenvironments. Tumours with low immune–ITH experience higher immunoselective pressure and escape via loss of heterozygosity in human leukocyte antigens and immunoediting. Instead, the tumours with high immune-ITH evolve to a more immunosuppressive/exhausted microenvironment. This gradient of immune pressure along with immune-ITH represents a hallmark of tumour evolution, which is closely linked to the transcriptome-immune networks that contributes to disease progression and immune inactivation. Remarkably, high immune-ITH and its transcriptomic signature are predictive for worse clinical outcome in HCC patients. This in-depth investigation of ITH provides evidence on tumour-immune co-evolution along HCC progression.

  Study EGAS00001003814

• Epigenetic and metabolomic data from type 2 diabetes adolescents

  Objective: Rates of type 2 diabetes (T2D) among adolescents are on the rise. Epigenetic changes could be associated with the metabolic alterations in adolescents with T2D. Methods: We performed a cross-sectional integrated analysis of DNA methylation data from peripheral blood mononuclear cells with serum metabolomic data from First Nation adolescents with T2D and controls participating in the Improving Renal Complications in Adolescents with type 2 diabetes through REsearch (iCARE) cohort study, to explore the molecular changes in adolescents with T2D. Results: Our analysis showed that 43 serum metabolites and 36 differentially methylated regions (DMR) were associated with T2D. Several DMRs were located near the transcriptional start site of genes with established roles in metabolic disease and associated with altered serum metabolites (e.g. glucose, leucine, and gamma-glutamylisoleucine). These included the free fatty acid receptor-1 (FFAR1), upstream transcription factor-2 (USF2), and tumor necrosis factor-related protein-9 (C1QTNF9), among others. Conclusions: We identified DMRs and metabolites that merit further investigation to determine their significance in controlling gene expression and metabolism which could define T2D risk in adolescents.

  Study EGAS00001003816

• DERMATLAS__Tubular_adenoma

  IVirtually 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 tubular adenomas.

  Study EGAS00001006686

• The miR-185/PAK6 Axis Predicts Therapy Response and Regulates Survival of Drug-Resistant Leukemic Stem Cells in Chronic Myeloid Leukemia

  Overcoming drug resistance and targeting cancer stem cells remain challenges for curative cancer treatment. To investigate the role of miRNAs in regulating drug resistance and leukemic stem cell (LSCs) fate, we performed global transcriptome profiling in treatment-naïve chronic myeloid leukemia (CML) stem/progenitor cells and identified that miR-185 levels anticipate their response to ABL tyrosine kinase inhibitors (TKIs). miR-185 functions as a tumor suppressor; its restored expression impaired survival of drug-resistant cells, sensitized them to TKIs in vitro, and markedly eliminated long-term repopulating LSCs and infiltrating blast cells, conferring a survival advantage in pre-clinical xenotransplantation models. Integrative analysis with mRNA profiles uncovered PAK6 as a crucial target of miR-185 and pharmacological inhibition of PAK6 perturbed the RAS/MAPK pathway and mitochondrial activity, sensitizing therapy-resistant cells to TKIs. Thus, miR-185 presents as a potential predictive biomarker, and dual targeting of miR-185-mediated PAK6 activity and BCR-ABL may provide a valuable strategy for overcoming drug resistance in patients.

  Study EGAS00001004153

• The clinical utility of genomics in childhood cancer extends beyond targetable mutations - WGS data

  We deeply sequenced 864 cancer-associated genes and the complete genomes and transcriptomes for 300 pediatric and adolescent/young adult (AYA) patients with poor prognosis or rare tumors. Using integrative somatic-germline analyses, we assessed the clinical utility of cancer genomics in this population. Clinically actionable variants were identified in 56% of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants, found in 54% of patients, were of unanticipated timing and type, with over 20% of these derived from the germline. Enrichment in corroborating mutational signature 3 (‘BRCAness’) in patients with germline homologous recombination defects suggests they are drivers of pediatric cancers and potential targets for PARP inhibition. Comprehensive somatic-germline cancer genomic profiling is useful at multiple points in the care trajectory for pediatric/AYA patients, supporting its integration into early clinical management. Re-biopsy at the time of relapse should be considered, and specific attention to mutational burden is indicated. Defective DNA repair genes in the germline may represent important targetable drivers.

  Study EGAS00001006610

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

  Cancer genomes harbor a broad spectrum of structural variants (SV) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read 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’ (TI thread), characterized by short (mostly <1kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50kbp in size. TI 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 cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

  Study EGAS00001006629

• The clinical utility of genomics in childhood cancer extends beyond targetable mutations - Cancer Panel data

  We deeply sequenced 864 cancer-associated genes and the complete genomes and transcriptomes for 300 pediatric and adolescent/young adult (AYA) patients with poor prognosis or rare tumors. Using integrative somatic-germline analyses, we assessed the clinical utility of cancer genomics in this population. Clinically actionable variants were identified in 56% of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants, found in 54% of patients, were of unanticipated timing and type, with over 20% of these derived from the germline. Enrichment in corroborating mutational signature 3 (‘BRCAness’) in patients with germline homologous recombination defects suggests they are drivers of pediatric cancers and potential targets for PARP inhibition. Comprehensive somatic-germline cancer genomic profiling is useful at multiple points in the care trajectory for pediatric/AYA patients, supporting its integration into early clinical management. Re-biopsy at the time of relapse should be considered, and specific attention to the mutational burden is indicated. Defective DNA repair genes in the germline may represent important targetable drivers.

  Study EGAS00001006642

• DERMATLAS__Tubular_adenoma_RNAseq

  IVirtually 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 tubular adenomas.

  Study EGAS00001006685

• Structural Expression of BMMF in tissues of colorectal, lung and pancreatic cancer patients

  Chronic inflammation, linked to the presence of bovine milk and meat factors (BMMF) and specific subsets of macrophages, results in oxygen radical synthesis and induction of mutations in DNA of actively replicating cells and replicating single stranded DNA (zur Hausen et al., 2017). Cancers arising from this process have been characterized as indirect carcinogenesis by infectious agents (without persistence of genes of the agent in premalignant or cancers cells). Here we investigate structural properties of pleomorphic vesicles, regularly identified by staining peritumour tissues of colorectal, lung and pancreatic cancer for expression of BMMF Rep. The latter represents a subgroup of BMMF1 proteins involved in replication of small single-stranded circular plasmids of BMMF, but most likely also contributing to pleomorphic vesicular structures found in the periphery of colorectal, lung and pancreatic cancers. Structurally dense regions are demonstrated in preselected areas of colorectal cancer, after staining with monoclonal antibodies against BMMF1 Rep. Similar structures were observed in human embryonic cells (HEK293TT) overexpressing Rep. These data suggest that Rep or Rep isoforms contribute to the structural formation of vesicles.

  Study EGAS00001006744

• Detection of isoforms and genomic alterations by high-throughput full-length single-cell RNA sequencing in HGSOC

  Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Long-read single-cell RNA sequencing (scRNA-seq), capturing full-length transcripts, lacked the depth to provide this information so far. Here, we increased the PacBio sequencing depth to 12,000 reads per cell, leveraging multiple strategies, including artifact removal and transcript concatenation, and applied the technology to samples from three human ovarian cancer patients. Our approach captured 152,000 isoforms, of which over 52,000 were novel, detected cell type- and cell-specific isoform usage, and revealed differential isoform expression in tumor and mesothelial cells. Furthermore, we identified gene fusions, including a novel scDNA sequencing-validated IGF2BP2::TESPA1 fusion, which was misclassified as high TESPA1 expression in matched short-read data, and called somatic and germline mutations, confirming targeted NGS cancer gene panel results. With multiple new opportunities, especially for cancer biology, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

  Study EGAS00001006807

• RNA sequencing of CCO- and CCO+ human hepatocytes

  RNAseq data from Passman et al 2023. Clonal CCO-deficient hepatocyte patches and nearby CCO-proficient hepatocytes were identified in morphologically normal human livers and sampled at varying distances along the PT-CV axis. Samples are named according to their location within the liver lobule, with "PT" denoting samples abutting the portal triad, "CV" denoting samples abutting the central hepatic vein, and "Mid" sampled acquired midway between these structures. Normal human liver is thought to be generally quiescent, however clonal hepatocyte expansions have been observed but neither their cellular source nor their expansion dynamics have been determined. Knowing the hepatocyte cell of origin, and their subsequent dynamics and trajectory within the human liver will provide an important basis to understand disease-associated dysregulation. Here we use in vivo lineage tracing and a combination of methylation sequence analysis to demonstrate normal human hepatocyte ancestry. We exploit next generation mitochondrial sequencing to determine hepatocyte clonal expansion dynamics across spatially-distinct areas of laser-capture microdissected clones, in tandem with computational modelling in morphologically-normal human liver.

  Study EGAS00001006984

• Relation between transcriptome, karyotype and age in cases of sex chromosome aneuploidies.

  Why women live longer than men is still an open question in human biology. Sex chromosomes have been proposed to play a role in the observed sex gap in longevity, and the Y male chromosome has been suspected of having a potential toxic genomic impact on male longevity. It has been hypothesized that transposable element (TE) repression declines with age, potentially leading to detrimental effects such as somatic mutations and disrupted gene expression, which may accelerate the aging process. Given that the Y chromosome is rich in repeats, age-related increases in TE expression could be more pronounced in males, likely contributing to their reduced longevity compared to females. In this work, we first studied whether TE expression is associated with the number of sex chromosomes in humans. We analyzed blood transcriptomic data obtained from individuals of different karyotype compositions: 46,XX females (normal female karyotype), 46,XY males (normal male karyotype), as well as males with abnormal karyotypes, such as 47,XXY, and 47,XYY.

  Study EGAS00001007462

• Effective reprogramming of patient-derived M2-polarized glioblastoma-associated microglia/macrophages by treatment with GW2580

  Purpose: Targeting immunosuppressive and pro-tumorigenic glioblastoma-associated macrophages and microglial cells (GAMs) has great potential to improve patient outcomes. Colony stimulating factor-1 receptor (CSF1R) has emerged as a promising target for reprogramming anti-inflammatory M2-like GAMs. However, treatment data on patient-derived, tumor-educated GAMs and their influence on the adaptive immunity are lacking. Experimental Design: CD11b+-GAMs freshly isolated from patient tumors were treated with CSF1R-targeting drugs PLX3397, BLZ945, and GW2580. Phenotypical changes upon treatment were assessed using RNAseq. Changes in GAM activation were confirmed in a complex patient-derived 3D tumor organoid model serving as a tumor avatar. Results: The most effective reprogramming of GAMs was observed upon GW2580 treatment, which led to the downregulation of M2-related markers, IL-6, and IL-10, ERK1/2 and MAPK signaling pathways, while M1-like markers and gene set enrichment indicating activated MHC-II presentation were substantially increased. Moreover, treatment of patient-derived glioblastoma organoids with GW2580 confirmed successful reprogramming.

  Study EGAS00001007466

• Human NXPE1 mediates variation in sialic acid O-acetylation in Colon Tissue

  Sialylation, and its modifications, have important consequences for development, differentiation and immune system regulation in nearly all human tissues. An accounting of all the genes responsible for modifying sialic acids is incomplete, with only one pathway for O-acetylation being well-defined in humans. Using whole genome sequencing of normal human colorectal tissue we identify a haplotype on chromosome 11 that is associated with differences in mild periodic acid Schiff staining (mPAS), a marker of acetylation of sialic acids. Of the genes in this region, only SNPs from NXPE1 correlate perfectly in a validation data set. NXPE1 protein expression by IHC also correlates with mPAS staining patterns in homozygous, as well as heterozygous colon crypts that have undergone spontaneous gain of the mPAS phenotype. Manipulation of NXPE1 expression with CRISPR-Cas9 causes changes to modified sialic acid levels in both 2D cell line and 3D normal colon organoid models. These findings suggest NXPE1 may cause O-acetylation of sialic acids in colorectal tissue, potentially defining NXPE as a new family of sialic acid-modifying genes.

  Study EGAS00001007704

• Sequential Antigen-loss and Branching Evolution in Lymphoma after Anti-CD19 and Anti-CD20 Targeted T Cell Engaging Immunotherapy

  CD19 CAR T cells and CD20 targeting T cell engaging bispecific antibodies have been approved in B-cell Non-Hodgkin lymphoma lately, heralding a new clinical setting where patients are treated with both approaches, sequentially. The aim of our study was to investigate the selective pressure of CD19 and CD20 directed therapy on the clonal architecture in lymphoma. Using a broad analytical pipeline, we identified truncating mutations in the gene encoding CD20 conferring antigen loss in 80% of patients relapsing from CD20 bispecs. Pronounced T cell exhaustion was identified in cases with progressive disease and retained CD20 expression. We also confirmed CD19 loss after CAR T cell therapy and report the case of sequential CD19 and CD20 loss. We observed branching evolution with re-emergence of CD20 positive subclones at later time points and spatial heterogeneity for CD20 expression in response to targeted therapy. Our results highlight immunotherapy as an evolutionary bottleneck selecting for antigen-loss variants but also complex evolutionary pathways underlying disease progression from these novel therapies.

  Study EGAS00001007561

• GOSH_Childhood_Leukemia_Behjati_WellcomeCore_RNA

  Paediatric single-cell archival project This project is taking sequential samples of children with B-cell acute lymphoblastic leukaemia (B-ALL) treated at Great Ormond Street Hospital (GOSH) between 2011 and the present day. Patients selected, are those that have done badly (e.g. had to escalate treatment, relapsed or died) despite having apparently ""good"" or ""intermediate"" risk cytogenetics, and often also had immunotherapy (e.g., the bi-specific T-cell engager blinatumomab) because of this. Using samples from key treatment-points, we hope to better characterise their disease and treatment response, to re-define understandings of risk for patients who did unexpectedly poorly. Additionally, we hope to define the disease/immunological factors that determine response to increasingly used immunotherapies within the field of paediatric B-ALL. Each of the samples will be viably frozen bone marrow-derived mononuclear cells. They will come in Eppendorf's, frozen and de-anonymised, from GOSH. We plan to primarily do single cell analysis (GEX, BCR +/-TCR), and then tumour-only WGS and bulk-RNA expression on left-over

  Study EGAS00001007616

• DERMATLAS__Leiomyoma_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 leiomyomas.

  Study EGAS00001007629

• DERMATLAS__Leiomyoma_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 leiomyomas.

  Study EGAS00001007630

• DERMATLAS__Poroma_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 poromas.

  Study EGAS00001007705

• MutWP5: CRUK Mutographs of Cancer: Cancer Mastectomy (WG)(Novaseq)

  Sequencing of LCM-derived microbiopsies from 40 women who underwent mastecomies due to breast cancer. LCM and sequencing will be conducted on both normal, unaffected breast, and, where possible, tumour tissue. Goal to assess the mutational burden, spectrum, and clonal dynamics within the tissue, and compare findings between the normal and associated cancer tissues. Whole-genome sequencing will be conducted on samples identified as promising from the initial targeted data. Results from this poriton of the study will be compared to women who had cosmetic breast reduction surgeries and those who are BRCA carriers. . This dataset contains all the data available for this study on 2023-03-08.

  Dataset EGAD00001010114

• The cellular immune and airway epithelial profile throughout childhood and in response to COVID-19 at multi-omic single cell level.

  Single-cell data gene expression data set (5’Chromium 10X) of healthy paediatric volunteers, and paediatric and adult COVID-19 patients. Gene expression was determined from samples of nasal, tracheal and bronchial brushings and blood (PBMCs). In addition to gene expression, PBMC’s were assayed by CITE-seq. A subset of samples have VDJ sequencing data for T cell receptors (TCR) and B cell receptors (BCR).

  Dataset EGAD00001007718

• AML Proteogenomic Landscape Whole-transcriptome RNA-Sequencing

  Whole transcriptome RNA-sequencing of purified bone marrow blasts of 136 de novo, treatment naive AML patients. For further details, we refer to the manuscript "The Proteogenomic Landscape of AML" by Jayavelu, Wolf, Buettner et al. mRNA extraction and whole transcriptome sequencing For transcriptome analysis the TruSeq Total Stranded RNA kit was used, starting with 250ng of total RNA, to generate RNA libraries following the manufacturer’s recommendations (Illumina, San Diego, CA, USA). 100bp paired-end reads were sequenced on the NovaSeq 6000 (Illumina) with a median of 57 mio. reads per sample. RNA Data Analysis Data quality control was performed with FastQC v0.11.9. Reads were aligned to the human reference genome (Ensembl GRCh38 release 82) using STAR v2.6.1. Gene count tables were generated while mapping, using Gencode v31 annotations. All downstream analyses were carried out using R v4.0 and BioConductor v3.12 (Huber et al., 2015; R Core Team, 2020). Size-factor based normalization was performed using DESeq2 v1.28.1(Love et al., 2014).

  Dataset EGAD00001008484

• Clonal dynamics after allogeneic haematopoietic cell transplantation using genome-wide somatic mutations - TGS

  Allogeneic haematopoietic cell transplantation (HCT) replaces the stem cells responsible for blood production with those harvested from a donor, and is received by 40,000 patients worldwide each year. To quantify dynamics of long-term stem cell engraftment, we sequenced whole genomes of 2,824 single-cell-derived haematopoietic colonies from blood samples of 10 donor-recipient pairs taken 9-31 years after HLA-matched sibling HCT. With younger donors, 10,000-50,000 stem cells had engrafted and were still contributing to haematopoiesis at time of sampling, but estimates were 10-fold lower with older donors. Engrafted stem cells made multilineage contributions to myeloid, B-lymphoid and T-lymphoid populations, although individual clones often showed biases towards one or other mature cell type. Recipients had lower clonal diversity than matched donors, equivalent to ~10-15 years of additional ageing, arising from up to 25-fold greater expansion of stem cell clones. An HCT-related population bottleneck alone could not explain these differences: instead, phylogenetic trees evinced two distinct modes of HCT-specific selection. In 'pruning selection', cell divisions underpinning recipient-enriched clonal expansions had occurred in the donor, preceding transplant - their selective advantage derived from preferential mobilisation, harvest, survival ex vivo or initial homing. In 'growth selection', cell divisions underpinning clonal expansion occurred through proliferative advantage in the recipient's marrow after homing - clones with multiple driver mutations especially demonstrated this pattern. Uprooting stem cells from their native environment and transplanting them to foreign soil exaggerates selective pressures, distorting and accelerating the loss of clonal diversity compared to the unperturbed haematopoiesis of donors.

  Dataset EGAD00001010874

• Deep whole genome ctDNA chronology of treatment-resistant prostate cancer

  Circulating tumor DNA (ctDNA) in blood plasma is an emerging tool for clinical cancer genotyping and longitudinal disease monitoring. We performed deep whole-genome sequencing of serial plasma and synchronous metastases in patients with aggressive prostate cancer. We comprehensively assess all classes of genomic alterations and demonstrate that ctDNA harbors multiple dominant populations whose evolutionary histories frequently indicate whole-genome doubling and shifts in mutational processes. Although tissue and ctDNA showed concordant clonally-expanded cancer driver alterations, each individual metastasis contributed only a minor share of total ctDNA. By comparing serial ctDNA before and after clinical progression on potent androgen receptor (AR) pathway inhibitors, we reveal population restructuring converging solely on AR augmentation as the dominant genomic driver of acquired treatment-resistance. Finally we leverage nucleosome footprints in ctDNA to infer mRNA abundance in synchronously biopsied metastases, including treatment-induced changes in AR pathway transcriptional activity.

  Dataset EGAD00001008460

• Genomic Advances in Sepsis (GAinS): RNA-seq

  Sepsis is defined as life-threatening organ dysfunction caused by a dysregulated host response to infection. This cohort comprises a subset of patients enrolled in the Genomic Advances in Sepsis (GAinS) study, an established biobank of adult sepsis patients. Patients with sepsis due to community acquired pneumonia or faecal peritonitis were recruited from 34 hospitals across the UK from 2005-2018, with samples for functional genomics and detailed clinical information collected on the first, third and/or fifth day following ICU admission. RNA was extracted from leukocytes isolated at the bedside using LeukoLOCK kits. We have previously identified sepsis response signatures (SRSs), transcriptomic endotypes that are associated with differential early mortality (Davenport et al, Lancet Respir Med, 2016; Burnham et al, AJRCCM, 2017) and response to treatment in a clinical trial (Antcliffe et al, AJRCCM, 2018). We generated RNA sequencing data on 903 samples, including 134 samples repeated from our previously released microarray data. Libraries were prepared using NEB Ultra II Library Prep kits (Illumina) and sequenced on a NovaSeq 6000. Reads were aligned to the reference genome (GRCh38) using STAR and gene counts quantified using featureCounts (annotation Ensembl v99). Counts were TMM-normalised and log-transformed. Following QC, processed data were available on 864 samples from 667 unique patients.

  Dataset EGAD00001008730

• Somatic mutation and selection at epidemiological scale - Sanger_NanoSeq_RandD

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

  Dataset EGAD00001015624

• cell-free DNA Whole Genome Sequencing of 3784 samples on MGI and Illumina platform

  The dataset of Integrative modeling of tumor genomes and epigenomes for enhanced cancer diagnosis by cell-free DNA includes 3784 whole genome sequencing bam files on the MGI and Illumina platform. The analyzed samples include plasma samples from normal individuals and patients with cancer.

  Dataset EGAD00001009335

• Whole-genome sequencing of gastric cancer 81 samples in the Japanese population.

  Genomic DNA of 81 cases from Japanese gastric cancer was extracted from tumor and matched normal tissues, and libraries with an insert size of 350–550 bp were prepared. The libraries were sequenced on a HiSeq 2500 instrument (Illumina) with paired-end reads of 101 bp. The read data are stored as FASTQ formatted files.

  Dataset EGAD00001008610

• Genomic Data of Pediatric MDS

  The dataset consists of - 126 whole exome sequencings (SAMD9/9Lmut: 64; GATA2mut 24, MDS wildtype 38/471) performed using SureSelect Human All Exon V6 enrichment (Agilent, cat# 5190-8863). The generated libraries were sequenced on the Illumina Hiseq 2500 with 150bp paired-end reads. FASTQ files were processed using SeqNext platform (JSI medical system, Germany), with gene-based alignment to a virtual panel of 300 genes (including 28 MDS-associated genes, SAMD9, and SAMD9L), consisting of genes relevant to bone marrow failure, MDS predisposition, and hematological cancers as per the Pan-Cancer studies with cohorts of >10,000 cancers. The respective BAM files are provided. - Custom panel targeting SAMD9, SAMD9L, and 22 single nucleotide polymorphisms (SNP) on chromosome 7q (allele frequency >35% in all ethnic sub-populations in gnomAD) (Ampliseq #IAD104171) were performed in 666/669 cases. And Custom panel targeting 28 MDS-associated genes (GATA2, RUNX1, HOXA9, CEBPA, GATA1, KRAS, NRAS, CBL, PTPN11, ASXL1, EZH2, SETBP1, FLT3, KIT, JAK2, JAK3, CSF3R, MPL, SH2, BCOR, BCORL1; RAD21, STAG2, CTCF, TP53, PTEN, CALR, VPS45) was performed in 544 cases (Ampliseq #IAD51150). Both custom panel libraries were prepared using NEBNext Ultra II DNA library prep kit (New England BioLabs, cat#E7645S/L) per manufacturer’s instruction and samples were sequenced on an Illumina Miseq 2000 with 2 x 150 bp reads. The respective BAM files are provided - 4 SAMD9/9L patients were subjected to MissionBio custom single-cell panel (CO-112) targeting 250 heterozygous gnomAD population polymorphisms on 7q arm and 69 amplicons in SAMD9/9L and other cancer genes. All libraries were sequenced on an Illumina NovaSeq6000 with 150 base-paired ending multiplexed runs. Fastq files were processed using the Tapestri Pipeline V2 and python-based Mosaic package (multi-omics analysis, data visualization). The derived BAM and loom files are provided.

  Dataset EGAD00001007856

• Somatic mutation and selection at epidemiological scale - TwinsUK_RENanoSeq_Buccal

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

  Dataset EGAD00001015621

• Single-cell profiling of neoantigen-specific T cells in lung cancers treated with neoadjuvant PD-1 blockade

  This dataset includes 406 samples from non small cell lung cancer patients treated with neoadjuvant anti-PD-1. The Single Cell 5’ V(D)J and 5’ DGE kits (10X Genomics) were used to capture immune repertoire information and gene expression from the same cell in an emulsion-based protocol at the single cell level. Libraries were generated and sequenced on an Illumina NovaSeq instrument using 2x150bp paired end sequencing.

  Dataset EGAD00001007728

• Dataset RNAseq from regions of insitu and invasive human mammary ductal disease

  BAM files from RNAseq study from regions of insitu and invasive human mammary ductal disease

  Dataset EGAD00001008586

• Whole genome sequencing of breast cancer and related clones

  Paired WGA (whole-genome amplificated) samples from 65 single cell-derived organoids and six fresh single cell-derived organoid samples both of which were established from normal mammary epithelial cells, 84 FFPE LCM (laser-capture microdissection) samples of breast cancer and related clones, 79 fresh-frozen LCM samples of non-cancer lobules of breast cancer patients, and 36 matched germline controls were subjected to whole genome sequencing using NovaSeq 6000 system (Illumina) or DNBSEQ-G400RS (MGI Tech).

  Dataset EGAD00001010833

• Dual targeting of polyamine synthesis and uptake in diffuse intrinsic pontine gliomas

  Bulk RNA-seq data of tumour data.

  Dataset EGAD00001006794

• Host factors dictate gut microbiome alterations in chronic kidney disease more strongly than to kidney function

  Despite recent progress, microbial markers in Chronic Kidney Disease (CKD) remain inconclusive. We combined quantitative faecal metagenomics (N=130 CKD cohort) and cross-study biomarker comparisons (Ntotal=4420) to study microbiome associations with glomerular filtration rate (eGFR; kidney function) and 4-year CKD progression. Transit time (TT) and medications primarily explained microbiome variation, trumping eGFR-related effects. Metagenome analysis showed increased p-cresol and indole synthesis potential and decreased plant-to-animal CAZyme ratios with lower eGFR, coupled to a Bact1 to Rum enterotype shift, suggesting saccharolytic-to-proteolytic microbiome transitions linked to dietary guidelines and slowed-down TT. Peritoneal dialysis patients showed increased (dysbiotic) Bact2 prevalence and linked intestinal inflammation. Only Escherichia coli, an unnamed Alistipes species and Bifidobacterium adolescentis were confounder-independent markers for eGFR, but these nor previous markers convincingly replicated across 11 studies. No predictors for CKD progression were found. Together, poor evidence for microbiome biomarkers but substantial eGFR-independent effects were found, suggesting potential for nutritional and TT interventions in CKD.

  Study EGAS50000000646

• Distinct genomic profiles and clinical outcomes in constitutional mismatch repair deficiency-associated high-grade gliomas: insights into mutational signatures and clonal evolution

  Constitutional mismatch repair deficiency (CMMRD) is a rare cancer-predisposition syndrome, yet the mutational processes and clonal evolution of CMMRD-associated high-grade gliomas (HGGs) remain incompletely understood. In this study, mutational signatures and tumor evolution were analyzed in 25 CMMRD-associated HGGs. Germline biallelic mutations were most frequently identified in MSH6 (56.0%), followed by PMS2 (36.0%) and MLH1 (8.0%). Affected patients showed very early disease onset (mean 5.8 years) and poor clinical outcomes, with a mean progression-free survival of 16 months. Distinct genomic and clinical differences were observed between MSH6- and PMS2-associated tumors. MSH6-deficient HGGs were enriched for SBS6, whereas PMS2-deficient tumors showed predominance of SBS21, and these mutational signatures were correlated with prognosis. Clonal evolution analyses revealed early involvement of POLE/POLD1 mutations and persistence of the founding clone at recurrence. Overall, the study highlights the importance of mutational signatures and clonal evolution in shaping tumor behavior and clinical outcomes in CMMRD-associated HGGs.

  Study EGAS50000001506

• Combined hereditary and somatic mutations of replication error repair genes result in rapid onset of ultra-hypermutated cancers

  DNA replication-associated mutations are repaired by two components: polymerase proofreading and mismatch repair. The mutational consequences of disruption to both repair components in humans are not well studied. We sequenced cancer genomes from children with inherited biallelic mismatch repair deficiency (bMMRD). High-grade bMMRD brain tumors exhibited massive numbers of substitution mutations (>250/Mb) – higher than all childhood and most cancers (>7000 analyzed). All ultra-hypermutated bMMRD cancers acquired early somatic driver mutations in DNA polymerases epsilon or delta. The ensuing mutation signatures and numbers are unique and diagnostic of childhood germline bMMRD (p<10e-13). Sequential tumor biopsy analysis revealed that bMMRD/polymerase mutant cancers rapidly amass an excess of simultaneous mutations (~600 mutations/cell division), reaching but not exceeding ~20,000 exonic mutations in <6 months. Thus implying a threshold compatible with cancer cell survival. We suggest a new mechanism of cancer progression in which mutations develop in a rapid burst after ablation of replication repair.

  Study EGAS00001001112

• TGF-β signaling mediates microglial resilience to spatiotemporally restricted myelin degeneration

  Microglia survey and regulate central nervous system myelination during embryonic development and adult homeostasis. However, whether microglia-myelin interactions are spatiotemporally regulated remains unexplored. By examining spinal cord white matter tracts, we found that myelin degeneration was particularly prominent in the dorsal column (DC) during normal aging. This was accompanied by molecular and functional changes in DC microglia as well as an upregulation of TGF-β signaling. Disrupting TGF-β signaling in microglia led to unchecked microglial responses and myelin loss in the DC, accompanied by neurological deficits exacerbated with aging. snRNA-seq revealed the emergence of a TGF-β signaling-sensitive microglial subset and a disease-associated oligodendrocyte subset, both spatially restricted to the DC. We further discovered that microglia rely on a TGFβ autocrine mechanism to prevent damage of myelin in the DC. These findings demonstrate that TGF-β signaling is indispensable for maintaining microglial resilience to myelin degeneration in the DC during the aging process. This highlights a previously unresolved checkpoint mechanism of TGF-β signaling with regional specificity and spatially restricted microglia-oligodendrocyte interactions.

  Study EGAS50000001413

• Djerroudi et al., E-cadherin inactivation shapes tumor microenvironment specificities in invasive lobular carcinoma

  Invasive lobular breast carcinoma (ILC) shows specific stromal features, T lymphocyte infiltration (TIL) being associated with poor prognosis. Here, we reveal the involved mechanism by performing single cell RNA sequencing, combined immunohistochemistry, deconvolution of bulk RNA sequencing from large retrospective ILC series, and functional assays using primary cells. We show that ILC accumulate FAP+ inflammatory cancer-associated fibroblasts (iCAF) through a previously undescribed mechanism mediated by E-cadherin/CDH1 on CAF plasticity. Indeed, CDH1 inactivation in ILC cancer cells prevents differentiation of iCAF into myofibroblastic CAF (myCAF), leading to iCAF accumulation. In turn, iCAF increase TIL infiltration and shape their spatial organization in ILC. Subsequently, CDH1-inactivated ILC cancer cells promote immune escape by lack of retention and activation of ITGAE-expressing resident memory CD8+ T lymphocytes (TRM). Hence, our study uncovers reciprocal interactions between CDH1-inactivated cancer cells, iCAF and CD8+ TRM, revealing why and how TILs have a poor prognosis in ILC patients, a mechanism extrapolated to other CDH1-mutated cancer types.

  Study EGAS50000000761