23249 results for "ega"

in 59.03 milliseconds.

• Early on-treatment changes in circulating tumor DNA fraction and response to enzalutamide or abiraterone in metastatic castration-resistant prostate cancer

  On-treatment circulating tumor DNA (ctDNA) quantification is touted as a surrogate biomarker for understanding metastatic cancer response to systemic therapy, but there is relatively little supporting data from standardized patient cohorts. Here, in a prospective observational cohort, we demonstrate that 4-week on-treatment changes in plasma ctDNA levels can accurately identify metastatic castration-resistant prostate cancer that will progress within 6 months of initiating first-line abiraterone or enzalutamide treatment. On-treatment ctDNA measurements outperformed existing clinical variables for predicting disease response and were also closely associated with differential overall survival. Although optimal timing of blood collections and thresholds for ctDNA detection will require further testing, our data highlights the potential for sequential ctDNA measurements to provide an early and reliable read-out for therapy response. The high accuracy of this information suggests that ctDNA could be useful for guiding clinical trials of therapy changes in initially non-responsive metastatic cancers.

  Study EGAS00001006856

• Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA

  These data relate to Abbosh et al 2023 'Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA'. Some patients were included in this work that were excluded from the other TRACERx 421 papers (Frankell et al and Al-bakir et al) often due to exclusion criteria in the TRACERx clinical trial. Reasons include partial contamination with C>A artefact mutational signature (CRUK0480 ,0490) which were included for ctDNA detection analyses using high confidence clonal mutations but not clonality analyses, incomplete resection at primary surgery (CRUK0291, 0234, 0230, 0387, 0622) and a concurrent oesophageal second primary cancer (CRUK0498). The primary tissue exome sequencing for these patients that were used in this study are included in this repository. Deep targeted sequencing of (median 200) tumour specific mutations in plasma for 1069 samples from 198 TRACERx patients are also included. Additional details related to these patients and samples can be found in the supplementary or data repositories from Frankell et al and Abbosh et al.

  Study EGAS00001006923

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

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

  Study EGAS00001006874

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

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

  Study EGAS00001006926

• Multiome_HB

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

  Study EGAS00001006932

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

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

  Study EGAS00001006945

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

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

  Study EGAS00001006970

• Study of complex rearrangements and mutational signatures in neuroblastoma heterogeneous risk groups.

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

  Study EGAS00001006983

• Neuroblastoma heterogeneity

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

  Study EGAS00001007019

• Epigenetic landscape of mixed phenotype leukemias

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

  Study EGAS00001007094

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

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

  Study EGAS00001007128

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

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

  Study EGAS00001007118

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

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

  Study EGAS00001007200

• An isogenic human iPSC model unravels neurodevelopmental abnormalities in SMA

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

  Study EGAS00001007259

• RNAseq of Colorectal cancer organoid-stroma biobank cohort

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

  Study EGAS00001007300

• WES of Colorectal cancer organoid-stroma biobank cohort

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

  Study EGAS00001007301

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

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

  Study EGAS00001007323

• Serial TERT rearrangement breakpoint quantification in circulating tumor DNA enables minimal residual disease monitoring in patients with neuroblastoma

  Approximately 30% of high-risk neuroblastomas harbor genomic TERT rearrangements at diagnosis, resulting in high TERT transcript levels and active telomerase. Combined with mutations in the RAS/MAPK/ALK signaling transduction network, these patients have a dismal prognosis. Our objective was to establish a liquid biopsy-based monitoring strategy for this particularly vulnerable patient subgroup, for whom no in-time molecular diagnostic tools exist to date. Detection and quantification of TERT rearrangement breakpoints in cell-free tumor DNA circulating in peripheral blood and bone marrow plasma from patients with TERT-driven neuroblastoma improved minimal residual disease monitoring and early relapse detection in individual patients. TERT rearrangement breakpoints, as a single marker or combined with mutations in the RAS/MAPK/ALK signaling transduction network, can be reliably analyzed from 1ng DNA using droplet digital PCR assays fit for clinical routine use. TERT rearrangement breakpoints are reflected in cell-free tumor DNA and can serve as robust biomarkers for disease activity.

  Study EGAS00001007365

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

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

  Study EGAS00001007373

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

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

  Study EGAS00001007510

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

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

  Study EGAS00001007611

• Unraveling mutagenic processes influencing the tumor mutational patterns of individuals with Constitutional Mismatch Repair Deficiency

  Constitutional mismatch repair deficiency (CMMRD), caused by bi-allelic germline variants in one of the mismatch repair (MMR) genes, is a childhood cancer predisposition syndrome that often results in the development of multiple tumors early in life. A better understanding of mutational processes driving subsequent CMMRD tumors could advance optimal treatment. Therefore, we performed a genomic characterization of 41 tumors from 17 individuals with CMMRD. Mutational patterns in these tumors were found to be influenced by multiple factors, including the affected MMR genes, tumor types and somatic polymerase proofreading mutations. Temozolomide treatment left prominent mutational signatures in two second primary hematologic malignancies. Furthermore, a novel indel signature was found in 54% of the tumors, characterized by one base pair cytosine insertions in cytosine homopolymers. In conclusion, the analysis of sequential tumors in individuals with CMMRD reveals diverse mutational patterns influenced by the underlying mutated MMR gene, tumor type and treatment history.

  Study EGAS00001007660

• Myeloid-specific KDM6B inhibition sensitizes Glioblastoma to PD1 blockade

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

  Study EGAS00001007002

• Northern Ireland COhort for the Longitudinal study of Ageing

  The Northern Ireland COhort for the Longitudinal study of Ageing (NICOLA) is a representative sample of ~8,500 people from across Northern Ireland. The study, which was set up in 2012, aims to understand what it is like to grow older in Northern Ireland. • NICOLA has a strong focus on molecular biomarkers so there is complementary genetic, epigenetic and transcriptomic data available for a subset of individuals. • We inherit much of our DNA from our parents while a small amount of this material changes as we get older. In NICOLA we have 551,839 directly genotyped and 18,148,478 imputed Single Nucleotide Polymorphisms (SNPs) currently available for 2969 participants. • Summary statistics for the association of these gene polymorphisms with ~30 phenotypes were generated. • Epigenetics provides a link between our inherited DNA and environmental influences from a person’s diet, medication and lifestyle. NICOLA has the epigenetic quality controlled profiles of 1984 individuals arising from variations in DNA methylation at 862,927 genetic sites.

  Study EGAS00001007915

• A cycling, progenitor-like cell population at the root of atypical teratoid rhabdoid tumor subtype differentiation trajectories.

  Background: Atypical teratoid rhabdoid tumors (ATRTs), characterized by the loss of SMARCB1, are among the most lethal pediatric central nervous system (CNS) tumors. Three molecular subtypes have been identified, each defined by distinct molecular and clinical features. No subtype-specific treatments are available, highlighting the necessity to better understand inter- and intra-subtype heterogeneity. Methods: We generated a single-nucleus transcriptome atlas of ATRTs, validated by single-cell ATAC-seq and spatial transcriptomics, to study subtype-specific differentiation trajectories. Results: We identified brain progenitor-like expression profiles within unique subtype-specific differentiation lineages. A shared cycling, progenitor-like cell population, interspersed throughout tumors, was observed across all ATRT samples. Moreover, we demonstrate that subtype-specific differentiation trajectories can be induced pharmacologically in ATRT tumoroids directing tumor cells towards a non-proliferative, mature phenotype. Conclusions: Our study reveals that ATRT subtypes exhibit signaling trajectories mirroring normal fetal brain development, enabling the development of maturation therapies tailored towards ATRTs.

  Study EGAS00001008108

• Cross-species Comparison Reveals Therapeutic Vulnerabilities Halting Glioblastoma Progression

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

  Study EGAS00001008155

• Reference epigenomes generated as part of the International Human Epigenomics Consortium (IHEC)

  The Centre for Epigenome Mapping Technologies (CEMT) is an epigenome sequencing platform funded by the Canadian Institutes of Health Research (CIHR) and Genome BC as part of the Canadian Epigenetics, Environment and Health Research Consortium (CEEHRC). CEMT aims to produce 100 human reference epigenomes, each comprising the following datasets: whole genome; whole genome bisulfite and oxidative bisulfite; transcriptome; microRNA; and chromatin immunoprecipitation (ChIP)-seq for the core chromatin marks (H3K4me1, H3K4me3, H3K9me3, H3K27me3, H3K36me3, and H3K27ac).CEEHRC is a full member of the International Human Epigenomics Consortium (IHEC), which aims to coordinate the production of 1,000 reference maps of human epigenomes for key cellular states relevant to health and diseases. As such, all datasets submitted by the CEMT platform conform to IHEC assay standards and metadata guidelines.

  Study EGAS00001000552

• Single-cell Kinnex sequencing of Alzheimer's disease isoform profile

  Alzheimer’s disease (AD) is the most common cause of dementia, yet its molecular causes remain incompletely understood. RNA diversity in part arising from dysregulated splicing may contribute to AD pathogenesis; however, the ability to interrogate the resulting full-length isoforms in single brain cells has been technologically limited. Here we report the use of PacBio Kinnex long-read sequencing combined with 10X Genomics single-cell preparations to identify both annotated and unannotated RNA isoforms. Eight AD and seven non-diseased post-mortem human brains yielded ~70,000 single nuclei showing diverse, differentially expressed, and switched transcripts in multiple genes. Cell-type-specific isoform expression and variants with intra-exonic junctions associated with reverse transcriptase-mediated somatic gene recombination were also detected. Novel isoforms, including CHI3L1 and SEPTIN4, were altered in AD. Kinnex sequencing of RNA isoforms from single nuclei detected vast isoform diversity amongst brain cell types, representing an under-explored element in AD and other brain disorders.

  Study EGAS50000001476

• Human islet 3D chromatin maps provide insights into type 2 diabetes

  The transcriptional output of each cell type is controlled by thousands of enhancers, many of which contain genetic risk variants for common diseases such as type 2 diabetes (T2D). To gain insight into how enhancer variation influences diabetes risk, we created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated islet enhancers with their target genes, often located hundreds of kilobases away. We identified hubs that show spatial and functional connections between enhancersand target genes related to islet function and diabetes. We demonstrate that genetic variants distributed across hub enhancers have a major impact on T2D heritability, and use this knowledge to identify individuals in whom islet regulatory variation has a prominent role in T2D risk. Our results demonstrate the importance of 3D chromatin architecture for molecular interpretation of T2D genetic association signals, and define genomic spaces that harbor a distinct component of the T2D polygenic burden.

  Study EGAS00001002917

• Motif-directed chromatin repression by BCL11B shapes ectopic targeting and lineage commitment

  BCL11B encodes a critical T-lineage transcription factor that is implicated as both an oncogene and a tumor suppressor in acute leukemia. Despite its central role in development and malignancy, a basic understanding of the sequence determinants of BCL11B activity remains obscure. Here, we define two functionally distinct modes of BCL11B engagement with chromatin: Most detectable BCL11B binding events are determined by pre-existing patterns of chromatin accessibility and RUNX1 occupancy. In contrast, BCL11B directs chromatin closing at genomic regions enriched for DNA motifs containing the core “TGACC” sequence, which are most prevalent in non-T lineages. These BCL11B-closed target sites cannot be identified from conventional BCL11B chromatin profiling approaches yet represent the primary mechanism by which BCL11B regulates gene expression. We provide evidence that cell of origin impacts the ability of BCL11B to repress lineage-inappropriate gene programs, thereby promoting lineage skewing. Together, these findings define a motif-centered determinant of BCL11B repression and explain how BCL11B controls lineage specification.

  Study EGAS50000001719

• Immune cell atlas of environmental and ancestral diversity in Indonesia [scRNAseq]

  Island Southeast Asia (ISEA) remains consistently underrepresented in large-scale human genomic resources despite its rich ancestral and lifestyle diversity. The region's complex population history and environmental/pathogen pressures could have shaped its populations' unique immune responses. Here, we generated single-cell PBMC profiles from 199 Indonesians sampled across four populations in the islands of Bali and New Guinea. The dataset consists of low-pass whole genome sequencing and pooled single-cell PBMC transcriptome sequencing. These groups capture diversity in regional genetic ancestries (West ISEA-like and Papuan-like) and lifestyle contrasts (urban versus rural communities in Bali; highland versus lowland communities in New Guinea), providing a window into how ancestry and environment jointly shape immune variation across ISEA. Together, this study provides the first large-scale single-cell immune atlas of ISEA, highlighting how ancestral and environmental diversity jointly shape human immunity in this globally important yet understudied region.

  Study EGAS50000001656

• Immune cell atlas of environmental and ancestral diversity in Indonesia [WGS]

  Island Southeast Asia (ISEA) remains consistently underrepresented in large-scale human genomic resources despite its rich ancestral and lifestyle diversity. The region's complex population history and environmental/pathogen pressures could have shaped its populations' unique immune responses. Here, we generated single-cell PBMC profiles from 199 Indonesians sampled across four populations in the islands of Bali and New Guinea. The dataset consists of low-pass whole genome sequencing and pooled single-cell PBMC transcriptome sequencing, including T-cell receptor (TCR). These groups capture diversity in regional genetic ancestries (West ISEA-like and Papuan-like) and lifestyle contrasts (urban versus rural communities in Bali; highland versus lowland communities in New Guinea), providing a window into how ancestry and environment jointly shape immune variation across ISEA. Together, this study provides the first large-scale single-cell immune atlas of ISEA, highlighting how ancestral and environmental diversity jointly shape human immunity in this globally important yet understudied region.

  Study EGAS50000001655

• Evolving epigenomics of immune cells in type 1 diabetes at single nuclei resolution

  The appearance of diabetes-associated autoantibodies is the first detectable sign of the disease process leading to type 1 diabetes (T1D). Evidence suggests that T1D is a heterogenous disease, where the type of antibodies first formed imply subtypes. Here, we followed 49 children, who subsequently presented with T1D and 49 matched controls, profiling single-cell epigenomics at different time points of disease development. Quantitation of cell and nuclei populations as well as transcriptome and open-chromatin states indicated robust, early, replicable monocyte lineage differences between cases and controls, suggesting heightened pro-inflammatory cytokine secretion early among cases. The order of autoantibody emergence in cases showed variation across lymphoid and myeloid cells, potentially indicating cellular immune response divergence. The strong monocytic lineage representation in peripheral blood immune cells before seroconversion and the weaker differential coordination of these gene networks close to clinical diagnosis emphasizes the importance of early life as a critical phase in T1D development.

  Study EGAS50000000863

• Single Cell RNA Transcriptomics of Mantle Cell Lymphoma Reveals the Presence of Treatment-Resistant Subclones at the Time of Diagnosis

  Mantle cell lymphoma is a chronically relapsing B cell neoplasm characterized by clonal heterogeneity. We implemented single cell RNA sequencing of paired tumor samples from four patients with an aggressive nodal MCL at diagnosis and after failure of standard immunochemotherapy. Clustering the cells based on the inferred CNVs showed small subclones present at diagnosis with a CNV profile similar to the major subclone at relapse. Comparison of these subclones to major subclones present at diagnosis revealed 119 upregulated pathways common for the four patient with an aggressive MCL, including NF-kappaB and PTEN signaling pathways. Additionally, we sequenced a patient with an indolent to blastoid MCL transformation with intestinal tissue involvement. In contrast to other compartments at relapse, the cells in the intestine showed strikingly different CNV profile and presence of new driver mutations including NOTCH1. In conclusion, our study demonstrates that scRNAseq can potentially detect therapy resistant subclones at diagnosis.

  Study EGAS50000000825

• scRNA-seq of LN and lymphoma stroma

  The activities of non-haematopoietic cells (NHCs), including mesenchymal stromal cells and endothelial cells, in lymphomas, reportedly underlie lymphomagenesis. However, understanding of lymphoma NHCs has been hampered by unexplained NHC heterogeneity, even in normal human lymph nodes (LNs). Here, we constructed a single-cell transcriptome atlas of >100,000 NHCs collected from 27 human samples, including LNs and various nodal lymphomas, revealing 30 distinct subclusters, including some that were previously unrecognized. Notably, this atlas was highly useful for comparative analysis with lymphoma NHCs, which revealed an unanticipated landscape of subcluster-specific changes in gene expression and interaction with malignant cells in follicular lymphoma NHCs, facilitating our understanding of stromal remodelling in lymphoma and potentially suggesting novel biomarkers. Our study largely updates NHC taxonomy in human LNs and pertains to analysis of disease status, providing a rich resource and deeper insights into LN and lymphoma biology and offering to advance lymphoma management and therapy.

  Study EGAS00001005732

• Immune variation leads to diverse outcomes in human malaria (2020-01-15)

  Falciparum malaria is clinically heterogeneous and yet in most cases the risk of life-threatening disease dramatically declines after the first few infections of life because children rapidly acquire disease tolerance (resistance to severe malaria without improved control of parasite burden). Identifying the factors that determine clinical outcome in a malaria-naive host is therefore paramount to reduce malaria mortality. However, the relative contribution of disease-causing variants of the Plasmodium var gene family versus pathogenic inflammatory cytokine cascades remains fiercely debated - we sought to reconcile these conflicting arguments by studying their interaction in vivo. To this end, two human challenge models were used to reveal the parasite-host interactions that underpin variation in falciparum malaria. To capture the diversity of human immune responses, each individual was analysed independently by tracking dynamic changes in their whole blood transcriptome through time. And to uncover evidence of preferential expansion of disease-causing variants, var gene expression was tracked in vivo from the start to end of infection. In this way, we could show that group A var genes are always expressed upon liver egress but in a minority population that does not increase over 10-days of blood cycling; there is no selection of disease-causing variants in the naive host. In fact, parasites do not respond in any way to differences or changes in host environment. On the other hand, host-intrinsic variation determines the intensity of inflammation and progression to clinical malaria. And furthermore, regulation of the interferon signaling network controls host fate. These data emphasise the role of human immune decision-making in shaping course & outcome of infection. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2020-01-15.

  Dataset EGAD00001005790

• GATA6 expression distinguishes classical and basal-like subtypes in advanced pancreatic cancer.

  PURPOSE: To determine the impact of basal-like and classical subtypes in advanced PDAC and to explore GATA6 expression as a surrogate biomarker. EXPERIMENTAL DESIGN: Within the COMPASS trial patients proceeding to chemotherapy for advanced PDAC undergo tumour biopsy for RNA sequencing. Overall response rate (ORR) and overall survival (OS) were stratified by subtypes and according to chemotherapy received. Correlation of GATA6 with the subtypes using gene expression profiling, in situ hybridization (ISH) were explored. RESULTS: Between December 2015-May 2019, 195 patients (95%) had enough tissue for RNA sequencing; 39 (20%) were classified as basal-like and 156 (80%) as classical. RECIST response data were available for 157 patients; 29 basal-like and 128 classical where the ORR was 10% vs. 33% respectively (p=0.02). In patients with basal-like tumours treated with modified FOLFIRINOX (mFFX) (n=22) the progression rate was 60% compared to 15% in classical PDAC (p= 0.0002). Median OS in the intention to treat population (n=195) was 9.3 months for classical vs. 5.9 months for basal-like PDAC (HR 0.47 95% CI 0.32-0.69, p=0.0001). GATA6 expression by RNAseq highly correlated with the classifier (p<0.001) and ISH predicted the subtypes with sensitivity of 89% and specificity of 83%. In a multivariable analysis, GATA6 expression was prognostic (p=0.02). In exploratory analyses, basal-like tumours, could be identified by keratin 5, were more hypoxic and enriched for a T cell inflamed gene expression signature. CONCLUSIONS: The basal-like subtype is chemoresistant and can be distinguished from classical PDAC by GATA6 expression.

  Dataset EGAD00001006081

• Preferential infiltration of unique Vγ9Jγ2‐Vδ2 T cells into glioblastoma multiforme

  Glioblastoma multiforme (GBM) is clinically highly aggressive as a result of evolutionary dynamics induced by cross-talk between cancer cells and a heterogeneous group of immune cells in tumor microenvironment. The brain harbors limited numbers of immune cells with few lymphocytes and macrophages; thus, innate‐like lymphocytes, such as γδ T cells, have important roles in antitumor immunity. Here, we characterized GBM‐infiltrating γδ T cells, which may have roles in regulating the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor‐infiltrating and blood‐circulating γδ T cells from four patients were analyzed by next-generation sequencing-based T-cell receptor (TCR) sequencing in addition to mutation and immune profiles in four GBM cases. In all tumor tissues, abundant innate and effector/memory lymphocytes were detected, accompanied by large numbers of tumor‐associated macrophages and closely located tumor‐infiltrating γδ T cells, which appear to have anti-tumor activity. The immune-related gene expression analysis using the TCGA database showed that the signature gene expression extent of γδ T cells were more associated with those of cytotoxic T and Th1 cells and M1 macrophages than those of Th2 cells and M2 macrophages. Although the most abundant γδ T cells were Vγ9Vδ2 T cells in both tumor tissues and blood, the repertoire of intratumoral Vγ9Vδ2 T cells was distinct from that of peripheral blood Vγ9Vδ2 T cells and was dominated by Vγ9Jγ2 sequences, not by canonical Vγ9JγP sequences that are mostly commonly found in blood γδ T cells. Collectively, unique GBM‐specific TCR clonotypes were identified by comparing TCR repertoires of peripheral blood and intra‐tumoral γδ T cells. These findings will be helpful for the elucidation of tumor-specific antigens and development of anticancer immunotherapies using tumor-infiltrating γδ T cells.

  Dataset EGAD00001004862

• The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

  De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

  Dataset EGAD00001007034

• Fetal origins of malignant germ cell tumours

  Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells (PGCs). They are uniquely able to generate embryonic and extraembryonic tissues, which in malignant GCTs carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we studied the phylogenetic and transcriptional diversity of 15 malignant gonadal GCTs and four normal testis biopsies by sequencing 131 whole genomes and 416 transcriptomes from 14 gonadal histologies, excised by laser capture microdissection. Our findings demonstrate that tumours were initiated by whole genome duplication likely in embryogenesis, within ~5-8 cell divisions post-PGC specification, followed by chromosome 12p gains associated with invasive disease. Of note, 12p imbalances were not only generated through GCT-typical isochromosomes, but also through non-isochromosomic configurations. Whilst tumours developed along homogenous phylogenetic pathways, they spawned manifold tissues independent of genetic subclonal diversification. A key feature of GCT tissues was the expression of fetal-specific genes. The transcriptional diversity notwithstanding, we found universal transcriptional elements correlated with hallmark 12p gains. Overall, our study reveals stereotyped phylogenies and transcriptomes underpinning the development of GCT that originate in fetal life and may lend themselves to therapeutic manipulation.

  Dataset EGAD00001007037

• Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - WGS

  Mutations in cancer-associated genes drive tumour outgrowth. However, the timing of driver mutations and dynamics of clonal expansion that lead to human cancers are largely unknown. We used 580,133 somatic mutations from whole-genome sequencing of 1013 clonal haematopoietic colonies to reconstruct the phylogeny of haematopoiesis, from embryogenesis to clinical disease, in 12 patients with myeloproliferative neoplasms which are blood cancers more common in older age. JAK2V617F, the pathognomonic mutation driving the majority of these cancers, was acquired in utero or childhood, with upper estimates of age of acquisition from 33 weeks gestation to 10.8 years, in all 5 patients in whom JAK2V617F was either the only or the first driver event.  Driver mutations associated with age-related clonal haematopoiesis occurred prior to or following JAK2V617F,  as independent clonal expansions in JAK2V617F-mutated patients, and as large clonal expansions in JAK2V617F-unmutated patients . These mutations were also acquired in utero or childhood, with DNMT3A mutations occurring by 8 weeks of gestation to 7.6 years across 4 patients, and PPM1D mutation occurring by age 5.8yrs in a patient with MPN lacking phenotypic driver mutations. Sequential driver mutation acquisition was common, separated by decades across life, and often outcompeted ancestral clones. The mean latency between JAK2V617F acquisition and clinical presentation was 30 years (range 11-54 years). Rates of clonal expansion were inferred from phylogenetic trees and varied substantially (3% to 190% expansion/year), were affected by additional driver mutations, and were predictive of latency to clinical presentation. Driver mutations and rates of expansion would have been detectable in blood one to four decades before clinical presentation. This study reveals how driver mutation acquisition early in life with life-long growth and evolution underlie adult myeloproliferative neoplasms, providing opportunities for early detection and intervention and a new paradigm for cancer development.

  Dataset EGAD00001007714

• Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - TGS

  Mutations in cancer-associated genes drive tumour outgrowth. However, the timing of driver mutations and dynamics of clonal expansion that lead to human cancers are largely unknown. We used 580,133 somatic mutations from whole-genome sequencing of 1013 clonal haematopoietic colonies to reconstruct the phylogeny of haematopoiesis, from embryogenesis to clinical disease, in 12 patients with myeloproliferative neoplasms which are blood cancers more common in older age. JAK2V617F, the pathognomonic mutation driving the majority of these cancers, was acquired in utero or childhood, with upper estimates of age of acquisition from 33 weeks gestation to 10.8 years, in all 5 patients in whom JAK2V617F was either the only or the first driver event.  Driver mutations associated with age-related clonal haematopoiesis occurred prior to or following JAK2V617F,  as independent clonal expansions in JAK2V617F-mutated patients, and as large clonal expansions in JAK2V617F-unmutated patients . These mutations were also acquired in utero or childhood, with DNMT3A mutations occurring by 8 weeks of gestation to 7.6 years across 4 patients, and PPM1D mutation occurring by age 5.8yrs in a patient with MPN lacking phenotypic driver mutations. Sequential driver mutation acquisition was common, separated by decades across life, and often outcompeted ancestral clones. The mean latency between JAK2V617F acquisition and clinical presentation was 30 years (range 11-54 years). Rates of clonal expansion were inferred from phylogenetic trees and varied substantially (3% to 190% expansion/year), were affected by additional driver mutations, and were predictive of latency to clinical presentation. Driver mutations and rates of expansion would have been detectable in blood one to four decades before clinical presentation. This study reveals how driver mutation acquisition early in life with life-long growth and evolution underlie adult myeloproliferative neoplasms, providing opportunities for early detection and intervention and a new paradigm for cancer development.

  Dataset EGAD00001007715

• 1M-scBloodNL

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

  Dataset EGAD00001007764

• pancreatic ductal adenocarcinoma exomes in study: Pro-immunogenic Impact of MEK inhibition combined with an anti-CD40 immunostimulatory antibody

  This dataset contains 3 pairs of exomes, germline (from whole blood) and patient-derived xenograft (PDX), from human pancreatic durctal adenocarcinoma patients. The data is referred to in the publication: "Pro-immunogenic impact of MEK inhibition synergizes with agonist anti-CD40 immunostimulatory antibodies in tumor therapy" (Nature Communications, 2020) Abstract: Cancer types with lower mutational load and a non-permissive tumor microenvironment are intrinsically resistant to immune checkpoint blockade. While the combination of cytostatic drugs and immunostimulatory antibodies constitutes an attractive concept for overcoming this refractoriness, suppression of immune cell function by cytostatic drugs may limit therapeutic efficacy. Here we show that targeted inhibition of mitogen-activated protein kinase (MAPK) kinase (MEK) does not impair dendritic cell-mediated T-cell priming and activation. Accordingly, combining MEK inhibitors (MEKi) with agonist antibodies (Abs) targeting the immunostimulatory CD40 receptor resulted in potent synergistic anti-tumor efficacy. Detailed analysis of the mechanism of action of MEKi GDC-0623 by means of flow cytometric analysis of the tumor immune infiltrate and whole tumor transcriptomics showed that, in addition to its cytostatic impact on tumor cells, this drug exerts multiple pro-immunogenic effects, including the suppression of M2-type macrophages, myeloid derived suppressor cells and CD4+ T-regulatory cells. In addition, MEKi was found to induce tumor-cell intrinsic interferon signaling, which contributed to antigen presentation by tumor cells. Finally, the tumoridical impact of MEKi involves the activation of multiple pro-inflammatory pathways involved in immune cell effector function in the tumor microenvironment. Our data therefore indicate that the combination of MEK inhibition with agonist anti-CD40 Ab is a promising therapeutic concept, especially for the treatment of mutant Kras-driven tumors such as pancreatic ductal adenocarcinoma.

  Dataset EGAD00001005938

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

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

  Dataset EGAD00001009049

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

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

  Dataset EGAD00001009087

• Distribution of ctDNA levels in plasma of early-stage non-small cell lung cancer patients measured using personalised ctDNA analysis

  Blood-based assays have shown increasing ability to detect circulating tumour DNA (ctDNA) in patients with early-stage cancer. However, detection of ctDNA in patients with non-small cell lung cancer (NSCLC) has continued to prove challenging. We performed retrospective analysis to quantify ctDNA levels in a cohort of 100 patients with early-stage NSCLC prior to treatment with curative intent. Where tumour tissue was available for whole exome sequencing, mutations identified were used to define patient-specific sequencing assays. For those 90 patients, plasma cell-free DNA was sequenced to high depth across capture panels targeting a median of 328 mutations specific to each patient. Data was analysed using Integration of Variant Reads (INVAR), detecting ctDNA in 66.7% of patients, including 52.7% (29 of 55) patients with stage I disease and >88% detection for patients with stage II and III disease (16/18 and 15/17). ctDNA was detected in plasma at fractional concentrations as low as 9.1x10-6, and in patients with tumour volumes as low as 0.23 cm3. A 36-gene sequencing panel (InVisionFirst-LungTM) was used to analyse plasma DNA in 27 samples including the 10 cases without tumour exome data, and detected ctDNA in 59% of samples tested (16 of 27). Across the entire cohort, detection rates were higher in squamous cell carcinoma patients compared to adenocarcinoma patients (81% vs. 59%). Detection of ctDNA prior to treatment was associated with significantly shorter time free from relapse, across all patients and in patient subgroups, with Hazard Ratios ranging from 2.25 to >11. Our analysis indicates that for patients with stage I NSCLC, the median ctDNA fraction in plasma is approx. 12 parts per million (0.0012%). This indicates the limits of detection that would be required for ctDNA-based liquid biopsies to detect ctDNA in the majority of patients with early-stage NSCLC.

  Dataset EGAD00001006230

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

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

  Dataset EGAD00001010876

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

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

  Dataset EGAD00001008412

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

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

  Dataset EGAD00001008589

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

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

  Dataset EGAD00001008706

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

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

  Dataset EGAD00001009304

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

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

  Dataset EGAD00001010135

• Mutation analysis of core biopsies from localized breast cancer patients

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

  Dataset EGAD00001008396

• The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

  De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

  Dataset EGAD00001007033

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

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

  Dataset EGAD00001009108

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

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

  Dataset EGAD00001008463

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedNanoSeq_Buccal

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

  Dataset EGAD00001015618

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedNanoSeq_Blood

  Utilising Targeted NanoSeq, we hope to explore the driver landscape in unprecedented detail in donors from TwinsUK. These donors have paired buccal swabs and data from this will help contextualise that work. 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 EGAD00001015619

• Somatic mutation and selection at epidemiological scale - TwinsUK_ExomeNanoSeq_Buccal

  Exome NanoSeq in buccal swab samples obtained from patients with a wide range of clinical phenotypes. 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 EGAD00001015620

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedEMSeq_Blood

  Development of a targeted methylation assay to determine the cell-type composition of a sample. 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 EGAD00001015623

• Genomics of Malignant Peripheral Nerve Sheath Tumor (MPNST)

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

  Dataset EGAD00001008608

• Ancestry and somatic profile predict acral melanoma origin and prognosis – WXS

  Acral melanoma, which is not ultraviolet (UV)-associated, is the most common type of melanoma in several low- and middle-income countries including Mexico. Latin American samples are significantly underrepresented in global cancer genomics studies, which directly affects patients in these regions as it is known that cancer risk and incidence may be influenced by ancestry and environmental exposures. To address this, we characterise the genome and transcriptome of 123 acral melanoma tumours from 92 Mexican patients, a population notable because of its genetic admixture. Compared with other studies of melanoma, we found fewer frequent mutations in classical driver genes such as BRAF, NRAS or NF1. While most patients had predominantly Amerindian genetic ancestry, those with higher European ancestry had increased frequency of BRAF mutations and a lower median number of structural variants. The tumours with activating BRAF mutations have a transcriptional profile more similar to cutaneous non-volar melanocytes, suggesting that acral melanomas in these patients may arise from a distinct cell of origin compared to other tumours arising in these locations. KIT mutations were found in a subset of these tumours, and quadruple wild-type samples (non BRAF/NRAS/NF1/KIT) differed from mutated samples in their structural genomic profile and overall and recurrence-free survival patterns. Transcriptional profiling defined three expression clusters; these characteristics were associated with recurrence-free and overall survival. We highlight potential novel low-frequency drivers, such as PTPRJ, NF2 and RDH5. Our study enhances knowledge of this understudied disease and underscores the importance of including samples from diverse ancestries in cancer genomics studies.

  Dataset EGAD00001015755

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - ExomeNanoSeq

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016058

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - TargetedNanoSeq

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016059

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - TargetedEMSeq

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016060

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - PTA_RNA

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016063

• Solve-RD_EuroNMD_cohort-1_DF1+2_V1

  Solve-RD data submitted to the ERN-EuroNMD cohort for re-analysis (Data freeze 1+2) v1

  Dataset EGAD00001009768

• Whole genome and Whole exome sequencing of patient-derived xenograft models of endometrial cancer

  Whole genome and Whole exome sequencing of patient-derived xenograft models of endometrial cancer

  Dataset EGAD00001008326

• WGS BAM files fromxa0Genomic Features and Classification of Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma

  Background and Aims: Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline(g) alterations in BRCA1, BRCA2 and PALB2. Methods: We interrogated whole genome sequencing (WGS) data on 391 patients including 49 carriers of pathogenic variants (PVs) in g_BRCA_ and PALB2. HRD classifiers were applied to the dataset and included: 1) the genomic instability score (GIS) used by Myriad MyChoice HRD assay; 2) substitution base signature 3 (SBS3); 3) HRDetect; and, 4) Structural Variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status. Results: Biallelic tumour inactivation of g_BRCA_ or PALB2 was evident in 43/49 germline carriers identifying HRD-PDAC. HRDetect (score ?0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumour classifiers suggested that 7-10% of PDAC that do not harbor g_BRCA/PALB2_ have features of HRD. Of the somatic HRDetect cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared to BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC. However in advanced disease the GIS (p=0.02), SBS3 (p=0.03) and HRDetect score (p=0.005) were predictive of platinum response and superior survival. PVs in g_ATM_ (n=6) or g_CHEK2_ (n=2) did not result in HRD-PDAC by any of the classifiers. In four patients, BRCA2 reversion mutations associated with platinum resistance. Conclusions: Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7-10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.

  Dataset EGAD00001007571

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

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

  Dataset EGAD00001011122

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - PTA_DNAHyb

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016062

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - LCMB_WGS

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016064

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - PTA_WGS

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016061

• Polyclonal selection of immune checkpoint mutations in thyroid autoimmunity - LCMB_WES

  Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.

  Dataset EGAD00001016065

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma_CM67-RNASeq

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (> median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Dataset EGAD00001006306

• Natural history of clonal haematopoiesis (2017-09-04)

  The incidence of acute myeloid leukemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 60. Only 10-15% of cases evolve from a pre-existing myeloproliferative or myelodysplastic disorder; the remaining cases arise de novo without a detectable prodrome and are diagnosed upon development of bone marrow failure. Analysis of diagnostic blood samples has demonstrated that de novo AML is preceded by the accumulation of somatic mutations in pre-leukemic hematopoietic stem and progenitor cells (preL-HSPCs) that subsequently undergo clonal expansion. If individuals in this pre-leukemic phase could be identified, methods for determination of risk and monitoring for progression to overt AML could be developed. However recurrent AML mutations also accumulate during aging in healthy individuals who never develop AML, referred to as age related clonal hematopoiesis (ARCH). To distinguish individuals with preL-HSPCs at high risk of developing AML from those with ARCH, we undertook deep targeted sequencing of genes recurrently mutated in AML in blood samples from 133 individuals in the European Prospective Investigation into Cancer and Nutrition (EPIC) study taken on average 6 years before they developed AML (pre-AML group), together with 683 matched healthy individuals (Control group). Pre-AML cases displayed accelerated age-correlated accumulation of somatic mutations.The identity, number and variant allele frequency (VAF) of mutations differed between the two groups, and were incorporated into a computational model of AML risk prediction that accurately distinguished pre-AML cases from controls on average 7 years prior to AML development. Our findings provide proof of concept that early prediction of AML development is feasible in high-risk populations, paving the way for early disease detection, monitoring, and potentially prevention.

  Dataset EGAD00001003703

• Comprehensive cancer predisposition gene testing in an adult multiple primary tumor series shows a broad range of deleterious variants and atypical tumor phenotypes

  Multiple primary tumors (MPT) affect a substantial proportion of cancer survivors and may result from various causes including inherited predisposition. Currently, germline genetic testing of MPT cases for cancer predisposition gene (CPG) variants is mostly targeted by tumor type. We ascertained pre-assessed MPT cases from genetics centers (defined as ≥2 primaries by age 60 years or ≥3 by 70) and performed whole genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment/molecular investigations, pathogenic variants in moderate and high-risk CPGs were detected in 67/440 (15.2%) of probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies including structural variants at low frequency (6/440 (1.4%) of probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant) the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice site CPG variants and at least one discordant tumor type was significantly higher than a control population (χ2=43.642 P=<0.0001). 2/67 (3%) of probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Summing together variant detection rates from a similarly ascertained previous MPT case series, the present results suggest that first-line comprehensive CPG analysis in a clinical genetics referral-based MPT cohort would detect a deleterious variant in about a third of cases.

  Dataset EGAD00001004088

• Spiradenocarcinoma (2018-10-29)

  Spiradenocarcinoma is a rare cutaneous sweat gland adnexal cancer with potential for aggressive behaviour. They are classified histologically into low- and high-grade tumours, with morphologically low-grade tumours thought to behave more favourably. However, limited information is available, with only 18 published cases. We have collected morphologically low-grade spiroadenocarcinomas (one with a lung metastasis) and high-grade spiroadenocarcinomas, as well as some spiradenomas (benign lesions), cylindromas (another type of malignant cutaneous sweat gland adnexal tumour) and hybrid spiradenoma-cylindromas. H&E-stained sections were reviewed, follow-up was obtained, and immunohistochemistry for Ki-67, p53 and, MYB has been performed. The tumours were solitary, measuring 0.8-7?cm (median: 2.7?cm), with a predilection for the head and neck of elderly patients (median age: 72 years; range 53-92) without gender bias. Histologically, the tumours were multinodular and located in deep dermis and subcutis. A pre-existing spiradenoma was present in all cases. The malignant component was characterized by expansile growth with loss of the dual cell population, up to moderate cytological atypia and increased mitotic activity (median: 10/10 HPF; range 1-28). Additional findings included squamoid differentiation (n=9), necrosis (n=7), and ulceration (n=5). P53 expression was variable and no significant differences were noted in the benign compared with the malignant parts of the tumours. In contrast, in the malignant components the Ki-67 proliferative index was slightly increased, and MYB expression was lost. Follow-up (median: 67 months; range: 13-132) available for 16 patients (84%) revealed a local recurrence rate of 19% but no metastases or disease-related mortality. Here we wish to exome sequence these cases to define the first genomic landscape for this malignancy. This dataset contains all the data available for this study on 2018-10-29.

  Dataset EGAD00001004426

• Molecular diagnosis of albinism (2018-03-14)

  Albinism is genetically heterogeneous rare genetic condition affecting 1:17000 in the Western world (but more frequent in Africa) whose main feature is a profound visual impairment, characterised by foveal hypoplasia, abnormal chiasmatic connections, nystagmus and photofobia. All these features result in severly altered visual acuity (<0,1), absent depth perception and poor night vision. People with albinism are primarily visually handicapped. In addition, for some types of albinism, the visual phenotype can be presented with partial or total hypopigmentation, hence resulting in a secondary phenotype which can lead to skin cancer if skin is not adequately protected. Recently a new syndrome has been described, FHONDA, with the same visual abnormalities of albinism but without pigment alteration. The traditional classification differentiates Oculoculatenous albinism (OCA), where hypopigmentation involves hair, skin and eyes versus Ocular Albinism (OA), where hypopigmentation only affects the eyes. These are non-sydrimic types of albinism. Some syndromic forms (Hermansky-Pudlak=HPS, Chediak-Higashi=CHS) affect cells beyond pigment cells, present in the lungs, immune system, platelets and intestines, resulting in more severe phenotypes that can be fatal. Mutations in at least 19 genes are assocaited with the corresponding types of albinism. Most hospitals will only diagnose the most frequent cases using traditional Sanger, MLPA approaches. Some will use CGH arrays. We aim to diagnose all cases of albinism through the Albinochip proposal, which combines a Sequenom first step of known mutations combined with subsequent NGS approaches. In some cases we fail to find a second mutation, these are good candidates for further full exome analyses. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2018-03-14.

  Dataset EGAD00001004039

• Genetic Heterogeneity of the familial gastric neuroendocrine tumors (2018-06-06)

  Gastric neuroendocrine tumors (gNETs) occur with an estimated frequency of 2 per 100,000 in the general population. Type I gastric neuroendocrine tumors (NETs) represent the 75% of gNTEs and arise from gastric enterochromaffin-like (ECL) cells. They have late age of onset and usually benigh course. Classically, hypergastrinemia in patients who have autoimmune atrophic gastritis, causes hyperplasia of gastric ECL cells that progresses into type I gastric NETs and parietal cell (PC) destruction. The genetic bases in families with this disease are unknown. We performed an exome sequencing study of an atypical aggressive familial gNETs case (with early age onset, nodal infiltrations and gastric adenocarcinomas) that followed a recessive model. We identified a deleterious mutation in homozygosis in the ATP4A gene, which encodes the proton pump responsible for acid secretion by gastric parietal cells. This mutation lead to achlorhydria first, and hypergastrinemia and gNET developing as consequence (Calvete et al. 2014). Recently, two more families with gNETs, classical clinical traits and recessive model have been studies by WES but we didn't find any mutation in the ATP4a gene. However, putative mutations affecting genes that contribute to the development and the integrity of PC have been found suggesting that genetic alterations associated to this disorder target to a unique cell type (parietal cells). In order to cinfirm this hypothesis, it is necessary the search for new genes implicated in the gNETs, more familial cases are needed to be studied. We have identified four more new familial gNETs cases. Here, we propose their study by WES. The first family is formed by thress siblings with gNETs. The other families include two siblings with gNETs. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2018-06-06.

  Dataset EGAD00001004153

• MINCR is a MYC-induced lncRNA able to modulate MYC’s transcriptional network in Burkitt lymphoma cells

  Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long non-coding RNAs (lncRNAs), on the way MYC is able to influence cellular transcriptome. To this aim we have intersected RNA-sequencing data from two MYC-inducible cell lines and from a cohort of 91 mature B-cell lymphomas carrying, or not carrying, genetic variants resulting in MYC over-expression. By this approach, we identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them we focused on a lncRNA that we named MINCR, for MYC-Induced long Non-Coding RNA, showing a strong correlation with MYC expression in MYC-positive lymphomas and also in pancreatic ductal adenocarcinomas. To understand its cellular role we performed RNA interference (RNAi) experiments and found that MINCR knock-down is associated with a reduction in cellular viability, due to an impairment in cell cycle progression. Differential gene expression analysis following RNAi showed a strongly significant enrichment of cell cycle genes among the genes down-regulate following MINCR knock-down. Interestingly these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of MYC transcriptional program. Accordingly, following MINCR knock-down, we observed a reduction in the binding of MYC to the promoters of selected cell cycle genes. Finally we provide evidences that down-regulation of AURKA, AURKB and CTD1 may explain the reduction in cellular proliferation observed upon MINCR knock-down. We therefore suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

  Dataset EGAD00001001441

• Determination of cell specific regulatory enhancers in hematopoetic models

  The BLUEPRINT project is a large-scale project investigating epigenetic mechanisms involved in blood formation, in health and disease. The human variation workpackage (WP10, led by NS) of the project seeks to characterize the effect of common sequence variation on the epigenome status of a cell. To do this, the project will use highly purified blood cells to minimise "experimental noise" and therefore enhance the power to discover modest effects. Two peripheral blood cell types, the CD14+CD16- monocyte (an important central orchestrator of adaptive immunity and a bridge between innate and adaptive immunity) and the CD65+CD9- neutrophilic granulocyte (the frontline cell for innate immunity) have been selected for this purpose. The two types of cells will be obtained at high purity from adult blood (AB) of 200 healthy males and females, respectively. Cells will be purified by using already validated and fully operational protocols that are based on density gradient centrifugation of the buffy coat obtained from whole blood, followed by magnetic bead-based purification using monoclonal antibodies against Cluster of Differentiation (CD) lineage-specific cell surface markers. Units of 475 ml of AB will be obtained from consenting volunteers of the Cambridge BioResource (CBR), a panel of 10,000 healthy volunteers local to Cambridge who have already consented to participate in biomedical research and of whom biological samples (DNA, plasma, serum) and lifestyle data have been deposited in a repository and database, respectively. We are requesting funding from the Human Diversity project to sequence the genomes of the 200 CBR volunteers at low pass (6x coverage). Nuclei, DNA and RNA will be recovered from the purified cells and made available for RNA-seq, DNA-seq and ChIP-seq and genomic DNA for entire genome sequencing will be recovered from the DNA repository.

  Dataset EGAD00001002205

• Constitutional and somatic genomic rearrangements coherently restructure chromosome 21 in acute lymphoblastic leukaemia

  Changes in gene dosage are a major driver of cancer1, engineered from a finite, but increasingly well annotated, repertoire of mutational mechanisms2-6. These processes operate over levels ranging from individual exons to whole chromosomes, often generating correlated copy number alterations across hundreds of linked genes. An example of the latter is the 2% of childhood acute lymphoblastic leukemia (ALL) characterized by recurrent intrachromosomal amplification of megabase regions of chromosome 21 (iAMP21)7,8 To dissect the interplay between mutational processes and selection on this scale, we used genomic, cytogenetic and transcriptional analysis, coupled with novel bioinformatic approaches, to reconstruct the evolution of iAMP21 ALL. We find that individuals born with the rare constitutional Robertsonian translocation between chromosomes 15 and 21, rob(15;21)(q10;q10)c, have ~2700-fold increased risk of developing iAMP21 ALL compared to the general population. In such cases, amplification is initiated by chromothripsis involving both sister chromatids of the dicentric Robertsonian chromosome. In contrast, sporadic iAMP21 is typically initiated by breakage-fusion-bridge (BFB) events, often followed by chromothripsis or other rearrangements. In both sporadic and iAMP21 in rob(15;21)c individuals, the final stages of amplification frequently involve large-scale duplications of the abnormal chromosome. The end-product is a derivative chromosome 21 or a derivative originating from the rob(15;21)c chromosome, der(15;21), respectively, with gene dosage optimised for leukemic potential, showing constrained copy number levels over multiple linked genes. In summary, the constitutional translocation, rob(15;21)c, predisposes to leukemia through a novel mechanism, namely a propensity to undergo chromothripsis, likely related to its dicentric nature. More generally, our data illustrate that several cancer-specific mutational processes, applied sequentially, can co-ordinate to fashion copy number profiles over large genomic scales, incrementally refining the fitness benefits of aggregated gene dosage changes.

  Dataset EGAD00001000658

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedEMSeq_Buccal

  Development of a targeted methylation assay to determine the cell-type composition of a sample. 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 EGAD00001015622

• Transcriptional profiling of tauopathies in human IPS-derived neurons (2019-08-21)

  This project is a pilot study, in collaboration with Maria Grazia Spillantini and Mariangela Iovino (Cambridge Centre for Brain Repair), to investigate the utility of IPS-derived neurons for the study of neurodegenerative disorders. Our aim is to characterise the transcriptional consequences of tauopathies using neurons derived from differentiated IPSCs as a model system. We will use IPS cells derived from six individuals, four with known mutations in the tau protein, 2 without. RNA will be extracted at Day 0 and Day 65 of differentiation by which time the neuronal tauopathy is apparent. RNA will be extracted and the transcriptome of each line characterised using RNAseq. We will then search for genes that are differentially expressed between the transcriptomes of individuals with tau mutations versus those in controls. My lab will analyse the RNAseq data, comparing both affected and controls and both time-points, to establish candidate genes. Darren Logan’s lab, along with our collaborators, will experimentally verify and further investigate these genes in additional lines and animal models. From this analysis we will generate a list of candidate genes that are differentially expressed between cases and controls. This study will not only help us understand the molecular basis of tauopathies, but also identify gene candidates for biomarkers of neurodegenerative disease. It will serve as a proof of principle for future planned studies into generating and transcriptomically analysing an allelic series of Tau mutations in IPSCs with a controlled genetic background. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2019-08-21.

  Dataset EGAD00001005277

• Identification of gene mutations and fusion genes in patients with Sézary Syndrome

  Sezary syndrome (SS) is a leukemic form of cutaneous T-cell lymphoma (CTCL) with an aggressive clinical course. The goal of our study is to understand the genetic basis of the disease by looking for driver gene mutations and fusion genes in 15 erythrodermic patients with circulating Sezary cels, fourteen of them fulfilling the diagnostic criteria of SS. Peripheral blood samples were collected for each patient and CD4+ T-lymphocytes and Granulocytes were isolated. DNA and RNA was extracted for matched normal and tumor materials. Whole exome sequencing was performed on matched normal and tumor tissues of 12 patients and RNA sequencing was performed on 10 patients. For 7 of these patients, we have both exome as well as RNA sequencing data. In this study, we present a complex genomic landscape of Sezary syndrome with several point mutations, copy number variations and fusion events which could contribute to the pathogenesis of SS.

  Study EGAS00001001706

• Sensitive and Frequent Identification of High Avidity Neo-epitope Specific CD8 + T-cells in Immunotherapy-naïve Ovarian Cancer

  Immunotherapy directed against private tumor neo-antigens derived from non-synonymous somatic mutations is a promising strategy of personalized cancer immunotherapy. However, feasibility in low mutational load tumor types remains unknown. Comprehensive and deep analysis of circulating and tumor-infiltrating lymphocytes (TILs) for neo-epitope specific CD8 + T cells allowed prompt identification of oligoclonal and polyfunctional such cells from most immunotherapy-naïve patients with advanced epithelial ovarian cancer studied. Neo-epitope recognition was discordant between circulating T cells and TILs, and was more likely to be found among TILs, which displayed higher functionalavidity and unique TCRs with higher predicted affinity than their blood counterparts. Our results imply that identification of neo-epitope specific CD8 + T cells is achievable even in tumors with relatively low number of somatic mutations, and neo-epitope validation in TILs extends opportunities for mutanome-based personalized immunotherapies to such tumors.

  Study EGAS00001002803

• Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia

  Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We therefore investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population, with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of the cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

  Study EGAS00001004407

• Assessment of cannabidiol and Δ9-tetrahydrocannabiol in mouse models of medulloblastoma

  Phytocannabinoids Δ9-tetrahydrocannabinol (THC) and cannabidiol (CBD) have been demonstrated to exhibit anti-cancer activity in preclinical models of brain cancer leading to new clinical trials for adults with glioblastoma. We describe here the first report that has investigated a role for THC and CBD in paediatric brain cancer. Cannabinoids had cytotoxic activity against medulloblastoma and ependymoma cells in vitro, functioning in part through the inhibition of cell cycle progression and the induction of autophagy. Despite these effects in vitro, when tested in orthotopic mouse models of medulloblastoma or ependymoma, no impact on animal survival was observed. Furthermore, cannabinoids neither enhanced nor impaired conventional chemotherapy in a medulloblastoma mouse model. These data show that while THC and CBD do have some effects on medulloblastoma and ependymoma cells, are well tolerated and have minimal adverse effects, they do not appear to elicit any survival benefit in preclinical models of paediatric brain cancer.

  Study EGAS00001004963

• Mult-omics Palbociclib Resistance Study in HR+/HER2– Metastatic Breast Cancer

  Resistance to CDK4/6 inhibitors and endocrine therapy (ET) are common and poorly understood since they have been real game changer in patients with HR-positive and HER2-negative MBC. A comprehensive genomic and transcriptomic analysis of pre-treatment and post-treatment tumors from patients treated with palbociclib plus ET identified novel markers associated with poor prognosis such as genomic scar features caused by homologous repair deficiency (HRD), estrogen response signatures, and four unique prognostic clusters with distinct molecular features. Tumors with TP53 mutations co-occurring with a unique HRD-high cluster responded poorly. By comparing genomic and transciptomic profiles of paired samples, tumors were found to be further enriched in HRD genomic scars and many had switched to aggressive molecular subtypes. Furthermore, we identified high frequencies of acquired genomic alterations upon disease progression in RB1, ESR1, PTEN, and KMT2C. Our findings provide new insight into potential predictive biomarkers that could be targeted to overcome resistance. (NCT03401359)

  Study EGAS00001005736

• Genome-wide Rare Variant Score Associates With Morphological Subtypes of Autism Spectrum Disorder

  Defining different genetic subtypes of Autism Spectrum Disorder (ASD) can help predict developmental outcomes. Based on minor physical and major congenital anomalies, we categorized 325 Canadian children with ASD into dysmorphic and non-dysmorphic subgroupings. We developed a method to calculate a genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data for each ASD proband. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD had a significantly higher GRVS compared to those with non-dysmorphic ASD (P= 0.027). Using the polygenic transmission disequilibrium test, we observed an over-transmission of ASD-associated common variants in non-dysmorphic ASD probands (P= 2.9×10-3). We replicated our findings using 442 ASD probands that had accompanying morphology data in the Simons Simplex Collection. Our findings provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

  Study EGAS00001005753

• Accumulation of copy number alterations and clinical progression across advanced prostate cancer

  The burden of copy number alterations positively associated with radiologically-evident distant metastases at diagnosis (P=0.00006) and showed a non-linear relationship with clinical outcome on univariable and multivariable analysis, characterised by a sharp increase in the relative risk of progression (P=0.003) and death (P=0.045) for each unit increase, stabilising into more modest increases with higher burdens. This association between copy number burden and outcome was similar in each of the four metastatic states. Copy number loss occurred significantly more frequently than gain at the lowest copy number burden quartile (q=4.1X10-6). Loss of segments in chromosome 5q21-22 and gains at 8q21-24, respectively including CHD1 and cMYC, occurred more frequently in cases with higher copy number alteration (for either region: Kolmogorov–Smirnov distance, 0.5; adjusted P<0.0001). Copy number alterations showed variability across tumor regions in the same prostate. This variance associated with increased risk of distant metastases (Kruskal-Wallis test P=0.037).

  Study EGAS00001006251

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

  The structure, fragmentation patterns and terminal sequences of cell-free DNA (cfDNA) are altered by nucleases and biological mechanisms in the blood of cancer patients. The cfDNA fragment-end composition recovered from low coverage WGS (<1 fold coverage) using a bespoke software (FrEIA) is aberrant in the plasma from cancer patient (n = 418, 655 samples) compared to controls (n = 117). As a standalone test FrEIA allows detection down to ~0.2% tumor fraction in vitro and in silico at 95% specificity, leading to a sensitivity of ~71% for detecting lung cancer (14/22 stage I-II, 27/38 stage III, 92/127 stage IV) and ~68% for detecting esophageal adenocarcinoma (26/44 stage II, 46/62 stage III). Additional cfDNA biological patterns can be combined with FrEIA increasing the diagnostic potential of low coverage WGS at minimal cost (mean AUROC = 0.96). Integrating multiple cfDNA biological signal augments the diagnostic performance of liquid biopsy.

  Study EGAS00001006142

• Disease specific alterations in the olfactory mucosa of patients with Alzheimer’s disease

  OOlfactory dysfunction manifests early in several neurodegenerative disorders. Olfaction is orchestrated by olfactory mucosal cells located in the upper nasal cavity. However, it is unclear how this tissue reflects key neurodegenerative features in Alzheimer’s disease. Here we report that Alzheimer´s disease olfactory mucosal cells secrete toxic amyloid-beta. We detail cell-type-specific gene expression patterns, unveiling 147 differentially expressed disease-associated genes compared to the cognitively healthy controls, and 5 distinct populations in globose basal cell -, myofibroblast-, and fibroblast/ stromal – like cells in vitro. Overall, coordinated alteration of RNA and protein metabolism, inflammatory processes and signal transduction were observed in multiple cell populations, suggesting a key role in pathophysiology. Our results demonstrate the potential of olfactory cell cultures in modelling Alzheimer´s disease. Moreover, for the first time we provide single cell data on olfactory mucosa in Alzheimer´s disease for investigating molecular and cellular mechanisms associated with the disease.

  Study EGAS00001006019

• An IL-1β driven neutrophil-stromal cell axis fosters a BAFF-rich microenvironment in multiple myeloma

  The bone marrow is continuously occupied by high numbers of neutrophils, and a tumor-supportive bias of these cells could significantly impact bone marrow-confined malignancies. In multiple myeloma, the bone marrow is characterized by inflammatory stromal cells with the potential to influence neutrophils. Here, we investigated myeloma-associated alterations in marrow neutrophils and the impact of stromal inflammation on neutrophil function. Mature neutrophils in myeloma marrow are activated and tumor-supportive, transcribing increased levels of pro-inflammatory cytokines, including IL-1β, and myeloma cell survival factors, such as the BCMA-ligand BAFF. Neutrophils were re-activated after first-line treatment, while this regimen reduced, but did not normalize, stromal inflammation. Interactions with inflammatory stroma induced neutrophil activation, including BAFF secretion, in a STAT3-dependent manner and once activated, neutrophils gained the ability to reciprocally induce stromal activation. Combined, our data define the presence of a neutrophil-stromal cell feed-forward loop driving tumor-supportive inflammation that could impact disease recurrence.

  Study EGAS00001007038

• Pediatric Patient-Derived-Xenograft development in MAPPYACTS – international pediatric cancer precision medicine trial in relapsed and refractory tumors

  Pediatric patients with recurrent and refractory cancers are in most need for new treatments. This study developed patient-derived-xenograft (PDX) models within the European MAPPYACTS cancer precision medicine trial (NCT02613962). To date, 131 PDX models were established following heterotopical and/or orthotopical implantation in immunocompromised mice: 76 sarcomas, 25 other solid tumors, 12 central nervous system tumors, 15 acute leukemias, and 3 lymphomas. PDX establishment rate was 43%. Histology, whole exome and RNA sequencing revealed a high concordance with the primary patient’s tumor profile, human leukocyte-antigen characteristics and specific metabolic pathway signatures. A detailed patient molecular characterization, including specific mutations prioritized in the clinical molecular tumor boards are provided. Ninety models were shared with the IMI2 ITCC Paediatric Preclinical Proof-of-concept Platform (IMI2 ITCC-P4) for further exploitation. This new PDX biobank of unique recurrent childhood cancers provides an essential support for basic and translational research and new treatments development in advanced pediatric malignancies.

  Study EGAS00001007327

• SMARCB1 loss activates patient-specific distal oncogenic enhancers in malignant rhabdoid tumors

  Malignant rhabdoid tumor (MRT) is a highly malignant and often lethal childhood cancer. MRTs are genetically defined by bi-allelic inactivating mutations in SMARCB1, a member of the BRG1/BRM-associated factors (BAF) chromatin remodeling complex. Mutations in BAF complex members are common in human cancer, yet their contribution to tumorigenesis remains in many cases poorly understood. Here, we studied derailed regulatory landscapes as a consequence of SMARCB1 loss in the context of MRT. Our multi-omics approach on patient-derived MRT organoids revealed a dramatic reshaping of the regulatory landscape upon SMARCB1 reconstitution. Chromosome conformation capture experiments subsequently revealed patient-specific looping of distal enhancer regions with the promoter of the MYC oncogene. This intertumoral heterogeneity in MYC enhancer utilization is also present in patient MRT tissues as shown by combined single-cell RNA-seq and ATAC-seq. We show that loss of SMARCB1 drives patient-specific epigenetic reprogramming underlying MRT tumorigenesis.

  Study EGAS00001007590

• Sequencing data from MethylScan assay on plasma

  Machine learning models in biomedical research are often hindered by demographic imbalances in clinical datasets, leading to biased predictions that disadvantage minority populations. Existing bias correction methods face limitations in handling the heterogeneity of biomedical data and the complexity of demographic influences. Here, we present DeBias, a computational framework for mitigating demographic biases in high-dimensional biomedical datasets. DeBias eliminates bias-associated subspaces within the entire feature space using control samples and enables global bias correction across all features, effectively preserving disease-specific signals while minimizing demographic distortions. To demonstrate its effectiveness, we applied DeBias to cfDNA methylation data for cancer detection. It achieved a >6-fold reduction in the number of features containing biases and outperformed existing methods in improving the cancer detection performance for minority populations. Our method represented a promising step toward developing machine learning models that are more accurate, reliable, and equitable for minority populations.

  Study EGAS00001008127

• ARHGAP11A Maintains Cortical Progenitor Identity Through RHOA–ROCK Signalling During Human Brain Development

  Mammalian brain development involves coordinated progenitor amplification, migration, and differentiation. Apical progenitors (APs) establish the ventricular zone (VZ) and are key determinants of brain size and complexity. Rho-family guanosine triphosphatases (Rho GTPases) regulate cytoskeletal dynamics and AP behaviour, but how specific Rho GTPase-activating proteins (Rho GAPs) influence progenitor identity and VZ architecture remains unclear. Using human forebrain organoids, we investigated the Rho GAP ARHGAP11A in human corticogenesis. CRISPR/Cas9-mediated ARHGAP11A knockout revealed its essential role in maintaining VZ integrity. Loss of ARHGAP11A impaired neuroepithelial organisation and randomised mitotic cleavage plane orientation via the RHOA–ROCK–actin axis. This led to premature AP delamination, AP depletion, and reduced cell density and glial numbers. Pharmacological inhibition of RHOA or ROCK rescued these defects, highlighting ARHGAP11A’s role in cytoskeletal remodelling to maintain cortical progenitors. Our findings establish ARHGAP11A as a critical regulator of AP identity and VZ integrity, with broad implications for human corticogenesis.

  Study EGAS00001008322