14658 results for "PLOS+ONE+Pennsylvania+German+ACS+language+code+pdc+open+access+study"

in 43.14 milliseconds.

• Combined clinical and gene expression score identifies high-risk individuals among follicular lymphoma patients on immunotherapy

  Follicular lymphoma (FL) is usually an indolent malignant B-cell lymphoma, but at least 30% of the patients experience transformation to Diffuse Large B-cell Lymphoma (DLBCL) with inferior survival. Early identification of patients with poor prognosis is essential for tailoring treatment. Here, we describe a network-based integrative analysis of 168 FL-biopsies from 127 patients with long clinical follow-up and transformation status. A distinct shift in expression towards a DLBCL expression phenotype is observed in BTK associated genes in cases with transformation. This shift is present years prior to transformation. An integrated score combining the BTK score and a clinical index is proposed, termed the BTK-FLIPI score. It identifies patients (16%) with very poor prognosis, with median time to progression 21.7 months (95% CI: 14.3-42.6) versus 70.9 months (95% CI: 41.5-95.0) for the remainder. The results suggest BTK as a potential target for treatment of high-risk FL-patients.

  Study EGAS00001002566

• The evolution of hematopoietic cells under cancer therapy

  Chemotherapies may increase mutagenesis of healthy cells and change the selective pressures in tissues thus influencing their evolution. However, their contributions to the mutation burden and clonal expansions of healthy somatic tissues are not clear. Here, exploiting the mutational footprint of some chemotherapies, we explore their influence on the evolution of hematopoietic cells. Cells of Acute Myeloid Leukemia (AML) secondary to treatment with platinum-based drugs showed the mutational footprint of these drugs, indicating that non-malignant blood cells received chemotherapy mutations. No trace of 5-fluorouracil (5FU) mutational signature was found in AML secondary to exposure to 5FU, suggesting that cells establishing the AML were quiescent during treatment. Using the platinum-based mutational signature as a barcode, we determined that secondary AMLs began their clonal expansion after the start of the cytotoxic treatment. Its absence in clonal hematopoiesis cases is consistent with a clonal expansion that predates the exposure to the cytotoxic agent.

  Study EGAS00001005234

• Myoepithelial progenitors as founder cells of hyperplastic human breast lesions upon PIK3CA transformation

  The myoepithelial (MEP) lineage of human breast comprises bipotent and multipotent progenitors in ducts and terminal duct lobular units (TDLUs). We here assess whether this heterogeneity impacts on oncogenic PIK3CA transformation. Single cell RNA sequencing (scRNA-seq) and multicolor imaging reveal that terminal ducts represent the most enriched source of cells with ductal MEP markers including a-smooth muscle actin (a-SMA), keratin K14, K17 and CD200. Furthermore, we find neighboring CD200high and CD200low progenitors within terminal ducts. When sorted and kept in ground state conditions, their CD200low and CD200high phenotypes are preserved. Upon differentiation, progenitors remain multipotent and bipotent, respectively. Immortalized progenitors are transduced with shp53 and mutant PIK3CA. Upon transplantation, CD200low MEP progenitors distinguish from CD200high by the formation of multilayered structures with a hyperplastic inner layer of luminal epithelial cells. We suggest a model with spatially distributed MEP progenitors as founder cells of biphasic breast lesions with implications for early detection and prevention strategies.

  Study EGAS00001005933

• Gain of function DNMT3A mutations cause microcephalic dwarfism and hypermethylation of Polycomb-regulated regions

  DNA methylation and Polycomb are key factors in the establishment of vertebrate cellular identity and fate. Here we report de novo missense mutations in DNMT3A, encoding the DNA methyltransferase DNMT3A, that cause microcephalic dwarfism, a hypocellular disorder of extreme global growth failure. Substitutions in the PWWP domain abrogate binding to the histone modifications H3K36me2/3, and alter DNA methylation in patient cells. Polycomb-associated DNA methylation canyons/valleys, hypomethylated domains encompassing developmental genes, become methylated with concomitant depletion of H3K27me3 and H3K4me3 bivalent marks. Such de novo DNA methylation occurs during differentiation of Dnmt3aW326R pluripotent cells in vitro, and is also evident in Dnmt3aW326R/+ dwarf mice. We therefore propose that the interaction of the DNMT3A PWWP domain with H3K36me2/3 normally limits DNA methylation of polycomb-marked regions. Our findings implicate the interplay between DNA methylation and polycomb at key developmental regulators as a determinant of organism size in mammals.

  Study EGAS00001003232

• Fibroblast heterogeneity drives metastatic spread in breast cancer through distinct mechanisms

  The major cause of death among breast cancer patients is metastatic dissemination, suggesting that understanding the mechanisms driving this process remains an unmet medical need. Here, by concomitantly using 5 markers, we identify 4 Carcinoma-Associated Fibroblasts (CAF) subsets in primary tumors and corresponding metastatic lymph nodes (LN). Two CAF subsets, called CAF-S1 and CAF-S4, correlate with the proportion of cancer cells in invaded LN and exert distinct functions in tumor cell spread. CAF-S1 promote cancer cell migration and initiates an epithelial to mesenchymal transition through a CXCL12 / TGFb positive loop. For their part, the highly contractile CAF-S4 foster cancer cell motility and invasion in 3-dimensional matrix via Notch pathway activation. Finally, CAF-S1 and CAF-S4 accumulation in LN is a new prognostic factor at diagnostic, independent of breast cancer subtypes and LN status, strengthening the validity of our findings in human breast cancer pathophysiology.

  Study EGAS00001003238

• Homologous recombination DNA repair deficiency and PARP inhibition activity in primary triple negative breast cancer

  Triple negative breast cancer (TNBC) encompasses molecularly different subgroups, with a subgroup harboring evidence of defective homologous recombination (HR) DNA repair. Within a phase 2 window clinical trial, RIO trial (EudraCT 2014-003319-12), we investigate the activity of PARP inhibitors in 43 patients with untreated TNBC. The primary endpoint, decreased Ki67, occurs in 12% of TNBC. Multiple different approaches identify HR deficiency in the majority of TNBC, 69% with the mutational-signature-based HRDetect assay. Cancers with mutational signatures of HR deficiency have a functional defect in HR. Following rucaparib treatment there is no association of Ki67 change with HR deficiency. In contrast, early circulating tumor DNA dynamics identify activity of rucaparib in cancers with mutation-signatures of HR deficiency, with rucaparib inducing expression of interferon response genes in HR deficient cancers. The majority of TNBCs have a defect in DNA repair, identifiable by mutational signature analysis, that may be targetable with PARP inhibitors.

  Study EGAS00001004405

• Genome-wide cell-free DNA mutational integration enables ultra-sensitive cancer monitoring

  In many areas of oncology, we lack sensitive tools to track low burden disease. While cell-free DNA (cfDNA) shows promise in detecting cancer mutations, we found that the combination of low tumor fraction (TF) and limited number of DNA fragments restricts low disease burden monitoring through the prevailing deep targeted sequencing paradigm. We reasoned that breadth may supplant depth of sequencing to overcome the barrier of cfDNA abundance. Whole genome sequencing (WGS) of cfDNA allowed ultra-sensitive detection, capitalizing on the cumulative signal of thousands of somatic mutations observed in solid malignancies, with TF detection sensitivity as low as 10-5. The WGS approach enabled dynamic tumor burden tracking and post-operative residual disease detection, associated with adverse outcome. Thus, we present an orthogonal framework for cfDNA cancer monitoring via genome wide mutational integration, enabling ultra-sensitive detection, overcoming the limitation of cfDNA abundance, and empowering treatment optimization in low-disease burden oncology care

  Study EGAS00001004406

• SARS‐CoV‐2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells

  The SARS‐CoV‐2 pandemic affecting the human respiratory system severely challenges public health and urgently demands for increasing our understanding of COVID‐19 pathogenesis, especially host factors facilitating virus infection and replication. SARS‐CoV‐2 was reported to enter cells via binding to ACE2, followed by its priming by TMPRSS2. Here, we investigate ACE2 and TMPRSS2 expression levels and their distribution across cell types in lung tissue (twelve donors, 39,778 cells) and in cells derived from subsegmental bronchial branches (four donors, 17,521 cells) by single nuclei and single cell RNA sequencing, respectively. While TMPRSS2 is strongly expressed in both tissues, in the subsegmental bronchial branches ACE2 is predominantly expressed in a transient secretory cell type. Interestingly, these transiently differentiating cells show an enrichment for pathways related to RHO GTPase function and viral processes suggesting increased vulnerability for SARS‐CoV‐2 infection. Our data provide a rich resource for future investigations of COVID‐19 infection and pathogenesis.

  Study EGAS00001004419

• Custom long non-coding RNA capture enhances detection sensitivity in different human sample types.

  Long non-coding RNAs (lncRNAs) are a heterogeneous group of transcripts that lack protein coding potential and display regulatory functions in various cellular processes. As a result of their cell- and cancer-specific expression patterns, lncRNAs have emerged as potential diagnostic and therapeutic targets. The accurate characterization of lncRNAs in bulk transcriptome data remains challenging due to their low abundance compared to protein coding genes. To tackle this issue, we describe a unique short-read custom lncRNA capture sequencing approach that relies on a comprehensive set of 565,878 capture probes for 49,372 human lncRNA genes. This custom lncRNA capture approach was evaluated on various sample types ranging from artificial high-quality RNA mixtures to more challenging formalin-fixed paraffin-embedded tissue and biofluid material. The custom enrichment approach allows the detection of a more diverse repertoire of lncRNAs, with better reproducibility and higher coverage compared to classic total RNA-sequencing.

  Study EGAS00001005418

• Detection of cell-free DNA fragmentation and copy number alterations in cerebrospinal fluid from glioma patients

  Glioma is difficult to detect or characterise using current liquid biopsy approaches. Detection of cell-free tumor DNA (cftDNA) in cerebrospinal fluid (CSF) has been proposed as an alternative to detection in plasma. We used shallow whole-genome sequencing (sWGS, at a coverage of < 0.4x) of cell-free DNA from the CSF of 13 patients with primary glioma to determine somatic copy number alterations and DNA fragmentation patterns. This allowed us to determine the presence of cftDNA in CSF without any prior knowledge of point mutations present in the tumor. We also showed that the fragmentation pattern of cell-free DNA in CSF is different from that in plasma. This low-cost screening method provides information on the tumor genome and can be used to target those patients with high levels of cftDNA for further larger-scale sequencing, such as by whole-exome and whole-genome sequencing.

  Study EGAS00001003255

• CCL22 chemokine mutations drive natural killer cell lymphoproliferative disease by biasing GPCR signaling

  Chronic lymphoproliferative disorder of natural killer cells (CLPD-NK) is characterized by clonal expansion of NK cells of incompletely understood genomic basis. Here, we report somatic mutations in the chemokine gene CCL22 as the hallmark of a distinct subset of CLPD-NK. CCL22 mutations were enriched at highly conserved residues, were mutually exclusive of STAT3 mutations, and were associated with gene expression similar to normal CD16dim/CD56bright NK cells. Mechanistically, the mutations resulted in ligand-biased chemokine receptor signaling, with decreased internalization of the G-protein-coupled receptor for CCL22, CCR4, via impaired β-arrestin recruitment. This resulted in increased cell chemotaxis in vitro, bidirectional crosstalk with the hematopoietic microenvironment, and enhanced NK cell proliferation in vivo in transgenic human IL-15 mice. Somatic CCL22 mutations illustrate a novel mechanism of tumor formation in which gain-of-function chemokine mutations promote tumorigenesis by biased G-protein-coupled receptor signaling and dysregulation of microenvironmental crosstalk.

  Study EGAS00001006009

• Molecular and functional profiling of plasmablastic lymphoma

  Plasmablastic lymphoma (PBL) represents a rare and aggressive lymphoma subtype frequently associated with immunosuppression. Clinically, patients with PBL are characterized by poor outcome. The current understanding of the molecular pathogenesis is limited. A hallmark of PBL represents its plasmacytic differentiation with loss of B-cell markers and, in 60% of cases, its association with Epstein-Barr virus (EBV). Roughly 50% of PBLs harbor a MYC translocation. Here, we provide a comprehensive integrated genomic analysis using whole exome sequencing (WES) and genome-wide copy number determination in a large cohort of 96 primary PBL samples. We identify alterations activating the RAS-RAF, JAK-STAT, and NOTCH pathways as well as frequent high-level amplifications in MCL1 and IRF4. The functional impact of these alterations is assessed using an unbiased shRNA screen in a PBL model. These analyses identify the IRF4 and JAK-STAT pathways as promising molecular targets to improve outcome of PBL patients.

  Study EGAS00001004659

• Enhanced cortical neural stem cell identity through SMAD/WNT inhibition in human cerebral organoids facilitates emergence of outer radial glial cells

  Cerebral organoids exhibit broad regional heterogeneity accompanied by limited cortical cellular diversity despite the tremendous upsurge in derivation methods, suggesting inadequate patterning of early neural stem cells (NSCs). Here we show that a short and early dual-SMAD/WNT inhibition course is necessary and sufficient for establishing robust and lasting cortical organoid NSC identity, efficiently suppressing of non-cortical NSC fates, while other widely used methods are inconsistent in their cortical NSC specification capacity. Accordingly, this method selectively enriches for outer radial glia (oRG) NSCs, which cytoarchitecturally demarcate well-defined outer sub-ventricular (oSVZ)-like regions propagating from superiorly radially organized, apical cortical rosette NSCs. Finally, this method culminates in the emergence of molecularly distinct deep and upper cortical layer neurons, and reliably uncovers cortex-specific microcephaly defects. Thus, a short SMAD/WNT inhibition is critical for establishing a rich cortical cell repertoire that enables mirroring fundamental molecular and cytoarchitectural features of cortical development and meaningful disease modeling.

  Study EGAS00001006063

• A comprehensive characterization of the cell-free transcriptome reveals tissue- and subtype-specific biomarkers for cancer detection

  Cell-free RNA (cfRNA) is a promising analyte for cancer detection. However, a comprehensive assessment of cfRNA in individuals with and without cancer has not been conducted. We performed the first transcriptome-wide characterization of cfRNA in cancer (stage III breast [n=46], lung [n=30]) and non-cancer (n=89) participants from the Circulating Cell-free Genome Atlas (NCT02889978). Of 57,820 annotated genes, 39,564 (68%) were not detected in cfRNA from non-cancer individuals. Within these low-noise regions, we identified tissue- and cancer-specific genes, defined as "dark channel biomarker" (DCB) genes, that were recurrently detected in individuals with cancer. DCB levels in plasma were correlated with tumor shedding rate and RNA expression in matched tissue, suggesting that DCBs with high expression in tumor tissue could enhance cancer detection in patients with low levels of circulating tumor DNA. Overall, cfRNA provides a unique opportunity to detect cancer, predict the tumor tissue of origin, and determine the cancer subtype.

  Study EGAS00001004704

• Whole genome sequencing of EBV Associated Nasopharyngeal Carcinoma

  Interplay between EBV infection and acquired genetic alterations during nasopharyngeal carcinoma (NPC) development remains vague. Here we report a comprehensive genomic analysis of 70 NPCs, combining whole-genome sequencing (WGS) of microdissected tumor cells with EBV oncogene expression to reveal multiple aspects of cellular-viral co-operation in tumorigenesis. Genomic aberrations along with EBV-encoded LMP1 expression underpin constitutive NF-κB activation in 90% of NPCs. A similar spectrum of somatic aberrations and viral gene expression undermine innate immunity in 79% of cases and adaptive immunity in 47% of cases; mechanisms by which NPC may evade immune surveillance despite its pro inflammatory phenotype. Additionally, genomic changes impairing TGFBR2 promote oncogenesis and stabilize EBV infection in tumor cells. Fine-mapping of CDKN2A/CDKN2B deletion breakpoints reveals homozygous MTAP deletions in 32-34% of NPCs that confer marked sensitivity to MAT2A inhibition. Our work concludes that NPC is a homogeneously NF-κB-driven and immune-protected, yet potentially druggable, cancer.

  Study EGAS00001004705

• single cell RNA-seq and ATAC-seq of human fetal forebrain tissue

  Oligodendrogenesis in the human central nervous system has been mainly observed at the second trimester of gestation, a much later developmental stage compared to mouse. Here we characterize the transcriptomic neural diversity in the human forebrain at post conceptual weeks (PCW) 8 to 10, using single-cell RNA-Seq. We find evidence of the emergence of a first wave of oligodendrocyte lineage cells as early as PCW 8, which we also confirm at the epigenomic level with single-cell ATAC-Seq. Using regulatory network inference, we predict key transcriptional events leading to the specification of oligodendrocyte precursor cells (OPCs). Moreover, by profiling the spatial expression of fifty key genes using In Situ Sequencing (ISS), we identify regions in the human ventral fetal forebrain where oligodendrogenesis first occurs. Our results indicate evolutionary conservation of the first wave of oligodendrogenesis between mouse and human and describe regulatory mechanisms involved in human OPC specification.

  Study EGAS00001006136

• Dissecting features of epigenetic variants underlying cardiometabolic risk using full-resolution epigenome profiling in regulatory elements

  Sparse profiling of CpG methylation in blood by microarrays have identified epigenetic links to common diseases. We apply methylC-capture sequencing (MCC-Seq) in a clinical population of ~200 adipose tissue and matched blood samples (Ntotal ~400), providing high-resolution methylation profiling (>1.3M CpGs) at regulatory elements. We link methylation to cardiometabolic risk through associations to circulating plasma lipid levels and identify lipid-associated CpGs with unique localization patterns in regulatory elements. We show distinct features of tissue-specific versus tissue-independent lipid-linked regulatory regions by contrasting with parallel assessments in ~800 independent adipose tissue and blood samples from the general population. We follow-up on adipose-specific regulatory regions under (1) genetic and (2) epigenetic (environmental) regulation via integrational studies. Overall, the comprehensive sequencing of regulatory element methylomes reveals a rich landscape of functional variants linked genetically as well as epigenetically to plasma lipid traits.

  Study EGAS00001003415

• Genomic analysis of a hypermutated gliosarcoma

  Gliosarcoma is a variant of glioblastoma with equally poor prognosis and characterized by mixed glial and mesenchymal pathology. Metastasis is not uncommon but involvement of the spinal cord is rare, and comprehensive genetic characterization of spinal gliosarcoma is lacking. We describe a patient initially diagnosed with a low-grade brain glioma via biopsy, followed by adjuvant radiation and temozolomide treatment. Nearly two years after diagnosis, she developed neurological deficits from an intradural, extramedullary tumor anterior to the spinal cord at T4, which was resected and diagnosed as gliosarcoma. Whole-exome sequencing (WES) of this tumor revealed a hypermutated phenotype, characterized by somatic mutations in key DNA mismatch repair (MMR) pathway genes, an abundance of C>T transitions within the identified somatic SNVs, and microsatellite stability, together consistent with temozolomide-mediated hypermutagenesis. This is the first report of a hypermutator phenotype in gliosarcoma, which may represent a novel genomic mechanism of progression from lower-grade glioma.

  Study EGAS00001004864

• The effect of pre-analytical and physiological variables on cell-free DNA fragmentation

  Background: assays that account for the biological properties and fragmentation of cell-free DNA (cfDNA) can improve the performance of liquid biopsy. However, pre-analytic and physiological differences between individuals on fragmentomic analysis are poorly defined. Methods: we analyzed the impact of collection tube, plasma processing time and physiology on the size distribution of cfDNA, their genome-wide representation and sequence diversity at the cfDNA fragment-ends using shallow Whole Genome Sequencing. Results: we observed that using different stabilizing collection tubes, or processing times does not affect the cfDNA fragment sizes, but can impact the genome-wide fragmentation patterns and fragment-end sequences of cfDNA. In addition, beyond differences depending on the gender, the physiological conditions tested between 63 individuals (age, body mass index, use of medication and chronic conditions) minimally influenced the outcome of fragmentomic methods. Conclusions: our results highlight that fragmentomic approaches have potential for implementation in the clinic, pending clear traceability of analytical and physiological factors.

  Study EGAS00001005748

• RNA-seq on patient-derived, stage II, CRC cell lines

  The gut microbiome is a key player in the immunomodulatory and pro-tumorigenic microenvironment in colorectal cancer (CRC). It is well established that Fusobacterium’s virulence factors are responsible for its pro-tumorigenic properties, however the role of its metabolic cross-talk with the tumor remains unexplored. Here we use in vitro, in vivo, and in silico methods to describe a novel microbe-host interaction keystone: formate. Omics data reveals molecular signatures, linking a metabolically driven CRC phenotype with Fusobacteria. We observed a metabolic shift towards increased F. nucleatum formate secretion in gut-on-chip model co-cultures with patient-derived CRC cells, along with an accelerated cancer glutamine metabolism. We show that high formate levels trigger AhR signaling, increase cancer stemness and cellular invasion. Finally, we observed an expansion of Th17 cells accompanying F. nucleatum-induced tumorigenesis. Moving beyond observational studies, we applied new experimental approaches for gaining ecosystem-level mechanistic understanding of F. nucleatum’s role in cancer pathogenesis.

  Study EGAS00001005948

• Loss of Epigenetic Barrier is Required for Enhancer Hijacking-Mediated Oncogenic Transcription

  Structural variants (SVs) involving enhancer hijacking can disrupt chromatin topologies to cause oncogene activation in cancer genomes, yet the molecular determinants for the transcriptional output of enhancer hijacking remain largely unknown. We developed a multimodal approach to integrate genome sequencing, chromosome conformation, and sequence-based deep learning for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified candidate pathogenic SVs including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for oncogenic activation of TLX3-dependent transcriptional programs. By engineering patient-associated t(5;14) in isogenic leukemia cells, we uncovered an uncharacterized mechanism whereby DNA methylation serves as an epigenetic barrier to enhancer hijacking and loss of epigenetic barrier is a molecular determinant for the transcriptional output of pathogenic SVs. Hence, leveraging the epigenetic barriers of SV-mediated oncogenic programs may provide new opportunities to reprogram gene regulation as epigenetic therapies in human disease.

  Study EGAS00001006140

• MutWP5__CRUK_Mutographs_of_Cancer__BRCA_Carriers___WG__Novaseq_

  Cancer is a genetic disease caused by an accumulation of mutations, however many of these mutations have been identified in pathologically normal tissue. We aim to use laser-capture microscopy (LCM) to sample individual clones from breast tissue to identify whether cancer-associated mutations appear in this normal tissue, assess the mutational burden present, and identify the mutational processes causing these mutations. We will sample from a wide age range of individuals (<20 to >60 years old) to determine whether these processes differ in pre- and post-menopausal women. We will also be comparing the tissue from healthy individuals (samples from breast reduction surgery) to those at elevated risk of breast cancer (mastectomy from BRCA1/2 patients) and those who have breast cancer (adjacent normal, distal normal, and tumour tissue from mastectomy). This will allow us to determine how these processes are different between these groups of individuals, and gain insight into the earliest stages of tumor development.

  Study EGAS00001003523

• Phase 1 CX-5461 Trial (Canadian Cancer Trials Group Trial IND.231)

  CX-5461 is a G-quadruplex stabilizer that exhibits synthetic lethality in homologous recombination (HR) deficient models. In this multicentre Phase Ib trial in patients with solid tumors, 40 patients were treated across 10 dose levels (50 - 650 mg/m2) to determine the recommended phase II dose (RP2D), and evaluate safety, tolerability, pharmacokinetics. Defective HR was explored as a predictive biomarker. CX-5461 was generally well tolerated, with a RP2D of 475 mg/m2 d1, 8 and 15 every 4 weeks, and dose limiting toxicities of phototoxicity. Responses were observed in 10% of patients, primarily in patients with defective HR (ORR 14%). Reversion mutations in PALB2 and BRCA2 were detected on progression following initial response in germline carriers, confirming the underlying synthetic lethal mechanism. In vitro characterization of UV sensitization showed this toxicity is related to the CX-5461 chemotype, independent of G4 synthetic lethality. These results establish clinical proof-of-concept for this first-in class G4 stabilizer.

  Study EGAS00001006173

• MutWP5__CRUK_Mutographs_of_Cancer__Cancer_Mastectomy__Exome__Novaseq_

  Cancer is a genetic disease caused by an accumulation of mutations, however many of these mutations have been identified in pathologically normal tissue. We aim to use laser-capture microscopy (LCM) to sample individual clones from breast tissue to identify whether cancer-associated mutations appear in this normal tissue, assess the mutational burden present, and identify the mutational processes causing these mutations. We will sample from a wide age range of individuals (<20 to >60 years old) to determine whether these processes differ in pre- and post-menopausal women. We will also be comparing the tissue from healthy individuals (samples from breast reduction surgery) to those at elevated risk of breast cancer (mastectomy from BRCA1/2 patients) and those who have breast cancer (adjacent normal, distal normal, and tumour tissue from mastectomy). This will allow us to determine how these processes are different between these groups of individuals, and gain insight into the earliest stages of tumor development.

  Study EGAS00001003528

• Arriba: accurate and efficient detection of gene fusions from RNA-Seq (H021)

  The identification of gene fusions from RNA sequencing data is a routine task in cancer research and precision oncology. However, despite the availability of many computational tools, fusion detection remains challenging. Existing methods suffer from poor prediction accuracy and are computationally demanding. We developed Arriba, a novel fusion detection algorithm with high sensitivity and short runtime. When applied to a large collection of published pancreatic cancer samples (n=803), Arriba identified a variety of driver fusions, many of which affected druggable proteins, including ALK, BRAF, FGFR2, NRG1, NTRK1, NTRK3, RET, and ROS1. The fusions were significantly associated with KRAS wild-type tumors and involved proteins stimulating the MAPK signaling pathway, suggesting that they substitute for activating mutations in KRAS. In addition, we confirmed the transforming potential of two novel fusions, RRBP1-RAF1 and RASGRP1-ATP1A1, in cellular assays. These results demonstrate Arriba’s utility in both basic cancer research and clinical translation.

  Study EGAS00001003554

• Widespread hypertranscription in aggressive human cancers

  Cancers are often defined by the dysregulation of specific transcriptional programs; however, the importance of global transcriptional changes is less understood. Hypertranscription is the genome-wide increase in RNA output. Hypertranscription’s prevalence, underlying drivers and prognostic significance are undefined in primary human cancer. This is due in part to limitations of expression profiling methods, which assume equal RNA output between samples. Here, we developed a computational method to directly measure hypertranscription in 7,494 human tumors, spanning 31 cancer types. Hypertranscription is ubiquitous across cancer, especially in aggressive disease. It defines patient subgroups with worse survival, even within well-established subtypes. Our data suggest that loss of transcriptional suppression underpins the hypertranscriptional phenotype. Single-cell analysis reveals hypertranscriptional clones, which dominate transcript production regardless of their size. Finally, patients with hypertranscribed mutations have improved response to immune checkpoint therapy. Our results provide fundamental insights into gene dysregulation across human cancers and may prove useful in identifying patients that would benefit from novel therapies.

  Study EGAS00001006365

• Resolving the immune landscape of human prostate at a single cell level in health and cancer

  The prostate gland produces prostatic fluid, high in zinc and citrate and essential for the maintenance of spermatozoa. Prostate cancer is a common condition with limited treatment efficacy in castration-resistant metastatic disease, including with immune checkpoint inhibitors. We used single-cell RNA-sequencing to perform an unbiased assessment of the cellular landscape of human prostate and identified a previously unappreciated subset of tumour-enriched androgen receptor-negative luminal epithelial cells, with increased expression of cancer-associated genes. We also found a variety of innate and adaptive immune cells in normal prostate that were transcriptionally perturbed in prostate cancer. An exception was a unique, prostate-specific, zinc transporter-expressing macrophage population (MAC-MT), that contributed to tissue zinc accumulation in homeostasis, but showed enhanced inflammatory gene expression in tumours, including T cell-recruiting chemokines. Remarkably, enrichment of the MAC-MT signature in cancer biopsies was associated with improved disease-free survival, suggesting beneficial anti-tumour functions.

  Study EGAS00001005787

• Massively parallel single-cell B-cell receptor sequencing enables rapid discovery of diverse antigen-reactive antibodies

  Using high-throughput single-cell B-cell receptor sequencing (scBCR-seq) we obtained accurately paired full-length heavy- and light-chain variable regions from thousands of individual B cells in a massively parallel fashion. We sequenced more than 250,000 B cells from rat, mouse and human repertoires to characterize their lineages and expansion. In addition, we immunized rats with chicken ovalbumin and profiled antigen-reactive B cells from lymph nodes of immunized animals. The scBCR-seq data recovered 81% (n = 56/69) of B-cell lineages identified from hybridomas generated from the same set of B cells that were subjected to scBCR-seq. Importantly, scBCR-seq identified an additional 710 candidate lineages that were not recovered as hybridomas. We synthesized, expressed and tested 93 clones from the identified lineages and found that 99% (n = 92/93) of the clones were antigen-reactive. Our results establish scBCR-seq as a powerful tool for antibody discovery.

  Study EGAS00001003663

• Molecular Profiling Reveals Unique Immune and Metabolic Features of Melanoma Brain Metastases

  There is a critical need to improve our understanding of the pathogenesis of -melanoma brain metastases (MBM). Thus, we performed RNA sequencing on 88 resected MBMs and 42 patient-matched extracranial metastases; tumors with sufficient tissue also underwent whole-exome sequencing, T-cell receptor sequencing, and IHC. MBMs demonstrated heterogeneity of immune infiltrates that correlated with prior radiation and post-craniotomy survival. Comparison with patient-matched extracranial metastases identified significant immunosuppression and enrichment of oxidative phosphorylation (OXPHOS) in MBMs. Gene-expression analysis of intracranial and subcutaneous xenografts, and a spontaneous MBM model, confirmed increased OXPHOS gene expression in MBMs, which was also detected by direct metabolite profiling and [U-13C]-glucose tracing in vivo IACS-010759, an OXPHOS inhibitor currently in early-phase clinical trials, improved survival of mice bearing MAPK inhibitor-resistant intracranial melanoma xenografts and inhibited MBM formation in the spontaneous MBM model. The results provide new insights into the pathogenesis and therapeutic resistance of MBMs.

  Study EGAS00001003672

• Renal Cell Carcinoma Tumors From S-TRAC trial (NCT00375674) Whole Transcriptome Sequencing

  We report a comprehensive genomic and transcriptomic analysis of tumor samples from 171 patients at high risk for recurrent renal cell carcinoma post nephrectomy from the S-TRAC trial (NCT00375674). We identify gene expression signatures, including STRAC11 (derived from the sunitinib-treated population). The overlap in key elements captured in these gene expression signatures, which include genes representative of the tumor stroma microenvironment, regulatory T cell, and myeloid cells, suggests they are likely to be both prognostic and predictive of the anti-angiogenic effect in the adjuvant setting. These signatures also point to the identification of potential therapeutic targets for development in adjuvant renal cell carcinoma, such as MERTK and TDO2. Finally, our findings suggest that while anti-angiogenic adjuvant therapy might be important, it may not be sufficient to prevent recurrence and other factors such as immune response and that tumor environment may be of greater importance.

  Study EGAS00001006528

• Multimodal Epigenetic Sequencing Analysis (MESA) of Cell-free DNA for Non-invasive Cancer Detection

  Multimodal characterization of cell-free DNA (cfDNA) in the blood can enable the sensitive and non-invasive detection of human cancers but remains technically challenging and costly. Here, we developed Multimodal Epigenetic Sequencing Analysis (MESA), a flexible and sensitive method of capturing and integrating multimodal epigenetic information of cfDNA using a single experimental assay, i.e., non-disruptive bisulfite-free methylation sequencing, such as Enzymatic Methyl-seq (EM-seq) and TET-assisted pyridine borane sequencing (TAPS). MESA can simultaneously infer cfDNA methylation, nucleosome occupancy, nucleosome fuzziness, and fragmentation profile for regions surrounding the promoters and polyadenylation sites (PASs). MESA’s integrated analysis of multimodal epigenetic features significantly improved the performance of cancer detection models compared to the usage of any single modality alone. MESA captures additional and highly complementary epigenetic information from cfDNA without additional experimental assays, highlighting the importance and clinical prospect of using multimodal epigenetic features for non-invasive cancer detection

  Study EGAS00001006462

• Whole Genome Sequencing on OCIAML-22

  Donor_OCIAML22_bulk: The donor (unsorted) that was used to generate OCI-AML22. OCI-AML22_bulk_1month_inVitro: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then sequenced. OCI-AML22_graft_CD34minus_WGS: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then injected in NSG-SGM3 mice. 21 weeks after injection, cells from 5 mice injected with 5 mice injected with 1475 cells per mice and 4 mice injected with 469 cells per mice were sorted. The human CD34- fraction was sequenced (CD45h+CD45m-7aad-annexinV-CD34-) OCI-AML22_graft_CD34positive_WGS: The donor (unsorted) that was used to generate OCI-AML22 was expanded for 1 month in vitro, then injected in NSG-SGM3 mice. 21 weeks after injection, cells from 5 mice injected with 5 mice injected with 1475 cells per mice and 4 mice injected with 469 cells per mice were sorted. The human CD34+ fraction was sequenced (CD45h+CD45m-7aad-annexinV-CD34+)

  Study EGAS00001006513

• Renal Cell Carcinoma Tumors From S-TRAC trial (NCT00375674) Exome Sequencing

  We report a comprehensive genomic and transcriptomic analysis of tumor samples from 171 patients at high risk for recurrent renal cell carcinoma post nephrectomy from the S-TRAC trial (NCT00375674). We identify gene expression signatures, including STRAC11 (derived from the sunitinib-treated population). The overlap in key elements captured in these gene expression signatures, which include genes representative of the tumor stroma microenvironment, regulatory T cell, and myeloid cells, suggests they are likely to be both prognostic and predictive of the anti-angiogenic effect in the adjuvant setting. These signatures also point to the identification of potential therapeutic targets for development in adjuvant renal cell carcinoma, such as MERTK and TDO2. Finally, our findings suggest that while anti-angiogenic adjuvant therapy might be important, it may not be sufficient to prevent recurrence and other factors such as immune response and that tumor environment may be of greater importance.

  Study EGAS00001006529

• Adult human kidney organoids originate from CD24+ cells and represent an advanced model for adult polycystic kidney disease

  Adult kidney organoids have been described as strictly tubular epithelial and termed tubuloids. While the cellular origin of tubuloids has remained elusive, here we report that they originate from a distinct CD24+ epithelial subpopulation. Long-term-cultured CD24+-derived tubuloids represent a functional human kidney tubule. We show that kidney tubuloids can be used to model the most common inherited kidney disease, namely autosomal dominant polycystic kidney disease (ADPKD), reconstituting the phenotypic hallmark of this disease with cyst formation. Single-cell RNA sequencing of CRISPR/Cas9 gene-edited PKD1 and PKD2 knockout tubuloids and human ADPKD and control tissue show similarities in upregulation of disease-driving genes. Furthermore, in a proof of concept, we demonstrate that tolvaptan, the only approved drug for ADPKD, shows a significant effect on cyst size in tubuloids but no effect in a pluripotent-stem-cell-derived model. Thus, tubuloids derive from a tubular epithelial subpopulation and represent an advanced model for ADPKD disease modeling.

  Study EGAS00001006551

• GRIDSS, PURPLE, LINX: Unscrambling the tumor genome via integrated analysis of structural variation and copy number

  We have developed a novel, integrated and comprehensive purity, ploidy, structural variant and copy number somatic analysis toolkit for whole genome sequencing data of paired tumor/normal samples. We show that the combination of using GRIDSS for somatic structural variant calling and PURPLE for somatic copy number alteration calling allows highly sensitive, precise and consistent copy number and structural variant determination, as well as providing novel insights for short structural variants and regions of complex local topology. LINX, an interpretation tool, leverages the integrated structural variant and copy number calling to cluster individual structural variants into higher order events and chains them together to predict local derivative chromosome structure. LINX classifies and extensively annotates genomic rearrangements including simple and reciprocal breaks, LINE, viral and pseudogene insertions, and complex events such as chromothripsis. LINX also comprehensively calls genic fusions including chained fusions. Finally, our toolkit provides novel visualisation methods providing insight into complex genomic rearrangements.

  Study EGAS00001004034

• Multi-omics identify falling LRRC15 as a COVID-19 severity marker and persistent pro-thrombotic signals in convalescence

  Patients with end-stage kidney disease (ESKD) are at high risk of severe COVID-19. Here, we perform longitudinal blood sampling of ESKD haemodialysis patients with COVID-19, collecting samples pre-infection, serially during infection, and after clinical recovery. Using plasma proteomics, and RNA-sequencing and flow cytometry of immune cells, we identify transcriptomic and proteomic signatures of COVID-19 severity, and find distinct temporal molecular profiles in patients with severe disease. Supervised learning reveals that the plasma proteome is a superior indicator of clinical severity than the PBMC transcriptome. We show that a decreasing trajectory of plasma LRRC15, a proposed co-receptor for SARS-CoV-2, is associated with a more severe clinical course. Strikingly, we observe that two months after the acute infection, patients still display dysregulated gene expression related to vascular, platelet and coagulation pathways, including PF4 (platelet factor 4), which may explain the prolonged thrombotic risk following COVID-19.

  Study EGAS00001006778

• Single-cell atlas of multiple myeloma and precursor diseases

  Single cell profiling of small cohorts of myeloma precursor disease patient samples including monoclonal gammopathy of unknown significance and smoldering multiple myeloma have shown diverse evolutionary patterns and immune changes that occur early in the disease process. Using scRNAseq and scBCRseq in a large cohort of fifty-three patients with myeloma precursor disease in comparison with myeloma and normal donors we describe the early genomic drivers of malignant transformation and describe their divergent clonal expansion in hyperdiploid compared to non-hyperdiploid samples. We describe intra-patient heterogeneity with potential therapeutic implications as well as the distinct evolution patterns (linear/branching) from myeloma precursor disease to myeloma. Finally, we describe the unique adaptations of the microenvironment as a response to distinct genomic changes in myeloma cells. These results further our knowledge to characterize myeloma precursor disease evolution, inform individual patient progression risk stratification and identify potential biomarkers that could be clinically exploited.

  Study EGAS00001006694

• Malignant progression of an ancestral bone marrow clone harboring a CIC-NUTM2A fusion in isolated myeloid sarcoma

  Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia (AML) or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis.

  Study EGAS00001006833

• 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

• 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

• 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 EGAS00001007016

• 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

• Single-cell RNA-seq of cervix and placenta

  Placental trophoblasts have been detected in cervical smears early in gestation, creating opportunities for non-invasive prenatal diagnosis. However, trophoblast isolation is limited by the lack of a cell catalogue and molecular profile of cervical smears. To establish an atlas and to explore the potential of single-cell RNA-sequencing to detect cervical trophoblasts, 10,539 single-cell transcriptomes from 12 non-invasive exocervical smears from pregnant women were profiled. In addition, we characterized 34,565 cells from six placental biopsies. We uncovered a novel extravillous trophoblast cell subtype characterized by epithelial marker genes and reduced HLA-G expression. Integration of both cell atlases demonstrated surprisingly similar expression profiles between maternal epithelial cells and placental extravillous trophoblasts, indicating that these trophoblasts retained epithelial properties without an invasive mesenchymal phenotype. Additionally, differential expression analysis identified novel markers discriminating cervical and placental cell types. Those analyses may contribute to trophoblast discovery and isolation from cervical smears.

  Study EGAS00001007044

• 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

• 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

• 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

• Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia

  In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness, and poor patient outcomes. A later arrest was associated with lineage fidelity, deep and durable leukemia remissions, and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1+ pre-leukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI-efficacy is an unexpected outcome of the differentiation stage at which this leukemia transforms.

  Study EGAS00001007167

• 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

• 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

• Repli-seq data for 'Replication timing alterations are associated with mutation acquisition during tumour evolution in breast and lung cancer'

  The essential process of DNA replication timing during each cell cycle is highly regulated ensuring the correct duplication of the genome. The extent and importance of alterations in this process during malignant transformation has not been extensively explored. Here, we evaluate the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find 6%-18% of the cancer genome displays ART. Genomic regions with a change from early to late replication exhibit an increased mutation rate and distinct mutational signatures, while regions changing from late to early contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART is a relatively early event during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

  Study EGAS00001007773

• Mutational Mechanisms in Multiply Relapsed Pediatric Acute Lymphoblastic Leukemia

  Pediatric acute lymphoblastic leukemia (ALL) is marked by low mutational load at initial diagnosis, which increases at relapse. To determine which processes are active in (relapsed) ALL and how they behave during disease progression before and after therapy, we performed whole genome sequencing on tumor samples of 29 multiply relapsed ALL patients. Mutational load increased upon relapse in 28 patients and every subsequent relapse in 22 patients, with UV-like damage, APOBEC activity, reactive oxygen species, thiopurine-associated damage and an unknown therapy component driving mutagenesis. Thiopurines were the most prominent source of new mutations in relapse, affecting over half of the studied patients and causing potential relapse-driving mutations in multiple patients. Mutational processes often affected patients over longer time periods, but could also occur in isolated events, suggesting the requirement of additional triggers. Our data pose mutational processes as prominent secondary drivers with an important role in ALL development and progression.

  Study EGAS00001007900

• Mapping the human hematopoietic stem and progenitor cell hierarchy through integrated single-cell proteomics and transcriptomics

  Single-cell transcriptomics (scRNA-seq) has enabled the characterization of cell state heterogeneity and recapitulation of differentiation trajectories. However, since proteins are the main functional entities in cells, the exclusive use of mRNA measurements comes at the risk of missing important biological information. Here we leverage recent technological advances in single-cell proteomics by Mass Spectrometry (scp-MS) to generate the first scp-MS dataset of an in vivo differentiation hierarchy encompassing over 2,500 human CD34+ hematopoietic stem and progenitor cells. Through integration with scRNA-seq, we identify proteins that are important for stem cell quiescence, which were not indicated by their mRNA transcripts, and demonstrate functional expression covariance during differentiation that is only detectable on protein level. Finally, we show that modeling translation dynamics can infer cell progression during differentiation and explain 45% more protein variation from mRNA than linear correlation. Our work serves as a framework for future single-cell multi-omics studies across biological systems.

  Study EGAS00001007930

• Novel paediatric case of a spinal high-grade astrocytoma with piloid features in a patient with Noonan Syndrome

  Noonan Syndrome (NS) is associated with an increased risk of low-grade central nervous system tumours in children but only very rarely associated with high-grade gliomas. Here we describe the first reported case of a spinal high-grade astrocytoma with piloid features (HGAP) in a child with NS. This case was a diagnostic and treatment dilemma, prior to whole-genome germline and tumour sequencing, tumour transcriptome sequencing and DNA methylation analysis. The methylation profile matched strongly with HGAP and sequencing identified somatic FGFR1 and NF1 variants and a PTPN11 germline pathogenic variant. Therapeutic targets were identified but also alterations novel to HGAP such as differential expression of VEGFA and PD-L1. The germline PTPN11 finding has not been previously described in individuals with HGAP.  This case underscores the power of precision medicine from a diagnostic, therapeutic and clinical management perspective, and describes an association between HGAP and NS which has not previously been reported.

  Study EGAS00001007937

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• Monosomy 7 delineates a primitive acute myeloid leukemia with adverse prognosis and responsiveness to epigenetic therapy

  The mechanistic basis for adverse risk associated with monosomy 7 remains undetermined in AML. In contrast to del(7q), monosomy 7 conferred adverse prognosis independent of TP53 and complex karyotype. Single-cell RNA and ATAC-seq revealed an epigenetically remodeled landscape and expansion of multipotent progenitor (MPP)-like cells. Monosomy 7 derived progenitors were distinguished by reduced chromatin accessibility and expression of STK17A and NT5C3A, which conferred resistance to cytotoxic and venetoclax-based therapy. To circumvent resistance mediated by monosomy 7 and deficient TP53, histone deacetylase inhibitors (HDACi) effectively restored apoptotic potential and augmented the activity of venetoclax. In a randomized clinical trial including older unfit patients with AML and monosomy 7, responses were observed in the majority that received low-dose cytarabine, venetoclax and the HDACi pracinostat, some with durable remissions exceeding 3 years. Longitudinal studies highlighted diminution of MPP-like cells, linking epigenetic remodeling to clinical response and therapeutic vulnerability of monosomy 7 AML.

  Study EGAS50000001641

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• Prematurity and Respiratory Outcome Program: Clinical Research Center Study of Functional and Lymphocytic Markers of Respiratory Morbidity in Hyperoxic Preemies

  Approximately 550,000 babies born prematurely each year in the United States suffer from birth at a time in development when the respiratory tract and immune system would normally be protected and maintained in a naïve state. This project is a component of the NIH Prematurity and Respiratory Outcomes Program (PROP) whose goals are the identification of disease mechanisms and biomarkers to stratify premature infants, at the time of discharge, for their risk of subsequent pulmonary morbidity. This Clinical Research Center (CRC) project will investigate prematurity-dependent alterations in cellular innate and adaptive immune systems resulting in increased susceptibility to respiratory infections and environmental irritants, and leading to respiratory morbidity in the first year of life. Prior studies have established developmental (maturity) and disease-related changes in circulating and pulmonary lymphocyte populations, but a comprehensive assessment of their relationship to disease risk/outcome has not been undertaken. We hypothesize that cellular and molecular immuno-maturity is altered due to intrinsic and extrinsic factors presented by premature birth in such a way as to reduce resistance to viral infections and to promote cytotoxic damage to the lung. We will evaluate immunologic maturity by comprehensively phenotyping lymphocyte populations in peripheral blood sampled at premature delivery, at the time of discharge from the hospital and at twelve months corrected age. The lymphocytic phenotype will be analyzed particularly in the context of gestational age and maternal-fetal stressors capable of modulating oxidative stress (oxygen exposure, infection and environmental tobacco smoke exposure). Additionally, we will assess changes in the molecular phenotype of isolated CD8 lymphocytes, a cell type preferentially recruited to the lungs of premature infants and capable of contributing to disease pathogenesis, by genome-wide expression profiling, in order to uncover novel disease pathways and define a gene expression signature associated with disease risk. We propose to build a statistical model, using cellular and molecular phenotypes and additional clinical variables, for stratifying risk of lung morbidity within the first year of life. Finally, we will assess the development of the gut microbiome in the preterm subjects to correlate with the observation of development of the adaptive immune system.

  Study phs001297

• Single Cell and Spatial Transcriptomics Studies of Fibrosis in Prospective Registry in IBD Study at MGH and GI Disease and Endoscopy Registry

  Stricturing disease or fibrosis is a common complication of inflammatory bowel disease (IBD), often requiring surgical intervention. This process is caused by the abnormal deposition of extracellular matrix in the context of inflammation, and thus involves crosstalk between multiple cell types including fibroblasts and immune cells. To understand the cellular interactions involved in fibrosis, we performed single cell and spatial transcriptomic analyses of biopsy or resection samples from IBD patients, diverticulitis patients or controls.For single cell analyses, intestinal samples were processed within three hours of collection by first stripping the epithelial layer with EDTA and then enzymatically digesting the underlying tissue. Cells from both the epithelial and tissue fractions were separately loaded in two channels on a 10x Chromium controller and libraries were generated with the 3'v3 Chemistry pipeline before sequencing on an Illumina sequencer. The resulting sequence data was demultiplexed and analyzed with the CellRanger pipeline to generate raw fastq files and processed gene-count matrices for each sample.For spatial transcriptomic analyses, intestinal resection samples were snap frozen in OCT (optimal cutting temperature embedding medium). Blocks were cut into 10μm sections in a cryostat and transferred to a 10x Genomics Visium slide v1. Spatially barcoded libraries were prepared according to manufacturer's instructions and sequencing on an Illumina sequencer. The resulting sequence data was processed with the SpaceRanger pipeline to generate fastq files and processed barcode-count matrices. Using the samples deposited here, we defined 68 fine-grained cell types across the ileum and colon of control and fibrotic intestinal tissues. We show that intestinal strictures were characterized by increased immune cells - including IgG+ plasma cells and CCR7hiCD4+ T cells - and inflammatory fibroblasts. We also define the transcriptional differences associated with fibrosis, with inflammatory fibroblasts in particular showing broad changes in gene expression. Spatial transcriptomics showed that key subsets colocalize within diseased tissues and identified additional populations such as interstitial cells of Cajal and enteric neurons. Further, we mapped gene expression onto intestinal biogeography, finding that known genetic risk loci are enriched within discrete spatial modules, defined by the presence of inflammatory fibroblasts and lymphoid follicles. Altogether, our datasets chart the key transcriptomic and cellular networks in stricturing CD and highlight the spatial organization of multicellular genetic risk factors.

  Study phs003943

• HCA_Adrenal_Foetal_WSSS_RNA_SB

  The Human Cell Atlas initiative is an ambitious global research effort which aims to describe every cell in the human body across all developmental stages to generate a reference map to accelerate progress in science and medicine particularly relating to human development and disease. The initiative is led by investigators based at the Wellcome Sanger Institute and the Broad Institute of MIT and Harvard and funding support from the Wellcome and NIH. This specific project will assess human development. We aim to interrogate the cellular composition and molecular regulators underpinning development and maturation through single cell RNA-sequencing/spatial transcriptomics and computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequent validate these findings in situ using allied imaging technologies e.g. immunohistochemistry and RNAscope and in vitro culture of human intestinal organoids (HIOs). This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004089

• Prediction of plasma ctDNA fraction and prognostic implications of liquid biopsy in advanced prostate cancer

  No consensus strategies exist for prognosticating metastatic castration-resistant prostate cancer (mCRPC). Circulating tumor DNA fraction (ctDNA%) is increasingly reported by commercial and laboratory tests but its utility for risk stratification is unclear. Here, we intersect ctDNA%, treatment outcomes, and clinical characteristics across 738 plasma samples from 491 male mCRPC patients from two randomized multicentre phase II trials and a prospective province-wide blood biobanking program. ctDNA% correlates with serum and radiographic metrics of disease burden and is highest in patients with liver metastases. ctDNA% strongly predicts overall survival, progression-free survival, and treatment response independent of therapeutic context and outperformed established prognostic clinical factors. Recognizing that ctDNA-based biomarker genotyping is limited by low ctDNA% in some patients, we leverage the relationship between clinical prognostic factors and ctDNA% to develop a clinically-interpretable machine learning tool that predicts whether a patient has sufficient ctDNA% for informative ctDNA genotyping (available online: https://www.ctDNA.org). Our results affirm ctDNA% as an actionable tool for patient risk stratification and provide a practical framework for optimized biomarker testing.

  Study EGAS50000000211

• Comprehensive spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer

  Although immunosuppressive and pro-metastatic functions of FAP+ Cancer-Associated Fibroblasts (CAF) are well-established, their plasticity and spatial distribution remain poorly understood. Here, we analyze trajectory inference, deconvolute spatial transcriptomics at single-cell level and perform functional assays to generate an integrative high-resolution map of breast cancer (BC), focusing particularly on the different subpopulations within inflammatory and myofibroblastic (iCAF/myCAF) FAP+ CAF. We identify 10 spatially-organized FAP+ CAF-related cellular niches, called EcoCellTypes, which are precisely localized within tumors. Consistent with their spatial organization, we identify DPP4- and YAP1-dependent mechanisms, by which cancer cells drive the transition of the detoxification-associated inflammatory FAP+ CAF cluster (Detox-iCAF) towards immunosuppressive extracellular matrix (ECM)-producing myofibroblasts (ECM-myCAF). In turn, ECM-myCAF polarize TREM2+ macrophages and regulatory NK cells to induce immunosuppressive EcoCellTypes. FAP+ CAF subpopulations accumulate differently depending on the invasive BC status and predict invasive recurrence of ductal carcinoma in situ (DCIS), which could help in identifying low-risk DCIS patients eligible for therapeutic de-escalation.

  Study EGAS50000000220

• Circulating Tumor Cell Heterogeneity in Neuroendocrine Prostate Cancer by Single Cell Copy Number Analysis

  Neuroendocrine prostate cancer (NEPC) is an aggressive variant of prostate cancer that may arise de novo or develop from pre-existing prostate adenocarcinoma as a mechanism of treatment resistance. Combined loss of tumor suppressors RB1, TP53, PTEN are frequent in NEPC but also present in a subset of prostate adenocarcinomas. Most clinical and preclinical studies support a trans-differentiation process, whereby NEPC arises clonally from a prostate adenocarcinoma precursor during the course of treatment resistance. Here we highlight a case of NEPC with significant intra-patient heterogeneity observed across metastases. We further demonstrate how single-cell genomic analysis of circulating tumor cells combined with phenotypic evaluation of cellular diversity can be considered as a window into tumor heterogeneity in patients with advanced prostate cancer. We performed whole genome sequencing of circulating tumor cells (CTCs) from patients with castration resistant adenocarcinoma (n=4) and neuroendocrine prostate cancer (n=8). CTCs were collected via Epic Sciences platform and single CTC genomic data (n=215) generated via Illumina NextSeq500. 

  Study phs002462

• Gut microbiome profiles according to sex, body mass index and dietary fiber intake

  Findings from recent studies suggest that the community of microbes residing in the human body is important in disease etiology; however, it remains unclear whether personal factors modulate human microbial composition. Studies based on animal models indicate that differences in composition might be attributed to sex-mediated effects. We analyzed the relationship of sex, adiposity, and dietary fiber intake with gut microbial composition using fecal samples from human subjects. We explored the associations of these factors with metrics of community composition and specific taxon abundances. We found that men and women had significantly different microbial community composition and that women had reduced abundance of a major phylum. Adiposity was associated with gut microbiome composition and specifically in women but not in men. Fiber from fruits and vegetables and fiber from beans were each associated with increased abundance of specific bacterial taxa. These findings provide initial indications that sex, adiposity, and dietary fiber might play important roles in influencing the human gut microbiome. Better understanding of these factors may have significant implications for gastrointestinal health and disease prevention.

  Study phs000884

• Genome-Wide Discovery of Novel Colon Cancer Predisposing Mutations

  Cancer risk shows clear heritability, but the inherited genetic factors remain largely unknown. This project will seek to identify new genes in which mutations confer hereditary risk for early onset colorectal cancer, with important implications for both biological understanding and clinical prediction and prevention. Identifying cancer predisposing mutations can provide new biological insights and significantly impact important clinical decision making regarding treatment, surveillance and preventive approaches. Hereditary colon cancer patients commonly present to Cancer Genetics clinics, but the majority of these patients do not have identifiable mutations in known candidate genes. We have performed whole exome sequencing to discover novel colorectal cancer susceptibility genes. The patients have been recruited from the Memorial Sloan Kettering Cancer Center (PIs Kenneth Offit, Nathan Ellis, and Zsofia Stadler), and Case Western University and University Hospitals of Cleveland (PIs Joseph Willis and Sanford Markowitz). All exome sequencing was performed at the Broad Institute of Harvard and MIT; samples sequence capture was performed using Agilent SureSelect Human All Exon Kit v2 and sequencing was performed on an Illumina HiSeq 2000 or 2500.

  Study phs000824

• The spectrum of somatic mutations in high-risk acute myeloid leukemia with -7/del(7q)

  Patients with myeloid malignancies bearing high-risk cytogenetic abnormalities lack effective therapies and have a poor overall survival. -7/del(7q) is identified in half of high-risk myeloid neoplasms. We recently identified CUX1 to be a haploinsufficient myeloid tumor suppressor gene located within the commonly deleted segment of 7q22. Here we identify the spectrum of somatic mutations that co-occur with loss of CUX1 and chromosome 7 in patients with de novo acute myeloid leukemia (AML) or a therapy-related myeloid neoplasm. -7/del(7q) leukemias have a distinct mutational profile characterized by low frequencies of alterations in major leukemogenic pathways, including genes encoding transcription factors, cohesin, and DNA-methylation-related proteins. In contrast, RAS pathway activating mutations occurred in 40% of -7/del(7q) samples, a significantly higher frequency than other AMLs and higher than previously reported. As targeted therapeutics advance, our data provide guidance for which pathways are most relevant in the treatment of adverse-risk myeloid leukemia.

  Study phs000759

• eMERGE Geisinger eGenomic Medicine (GeM) Abdominal Aortic Aneurysm Project (AAAP)

  A large research cohort of Geisinger Abdominal Aortic Aneurysm (AAA) patients was created by enrolling and consenting patients of the Geisinger Department of Vascular Surgery. Consented patients provide blood, serum and DNA samples for research and authorize use of data in their medical record for research. They also complete a data questionnaire that asks information about family history of AAA and other vascular diseases, as well as information on known or suspected AAA risk factors, including smoking history, body mass index, hypertension, type 2 diabetes, and atherosclerotic disease. The AAA cases are 78% male, with a mean age of 74 years; 39% have undergone surgical or endovascular repair; more than 85% or current or past smokers; approximately 20% of cases report a positive family history of AAA. Geisinger AAA patients undergo regular imaging studies, typically every 6 months, to monitor the progression of aneurysm expansion. This allows the growth rates of their aneurysm to be calculated. DNA samples from a total of 910 AAA patients were used for whole genome genotyping; these results are included in the dbGaP.

  Study phs000387

• Induction of HIV-1 Env-Specific Peripheral Blood Plasmablasts and Clonal Persistence as Bone Marrow Plasma Cells Following Vaccination in Humans

  Induction of persistent HIV-1 Envelope (Env) specific antibody (Ab) is a primary goal of HIV vaccine strategies, however, it is unclear if HIV Env immunization in humans induces bone marrow plasma cells, the presumed source of long-lived systemic Ab. To define the features of Env-specific plasma cells after vaccination, samples were obtained from HVTN 105, a phase I trial testing the same gp120 protein immunogen, AIDSVAX B/E, as used in RV144, along with a DNA immunogen in various prime/boost strategies. Boosting regimens that included AIDSVAX B/E induced robust peripheral blood plasmablast responses. The Env-specific immunoglobulin repertoire of the plasmablasts was dominated by VH1 gene usage and targeting of the V3 region. Numerous plasmablast-derived immunoglobulin lineages persisted in the bone marrow greater than eight months after immunization, including in the CD138+ long-lived plasma cell compartment. These findings identify a cellular linkage for the development of sustained Env-specific Ab following vaccination in humans.

  Study phs002027

• Longitudinal genome-wide analysis of progressive chronic lymphocytic leukemia under uniform front-line therapy of pentostatin, cyclophosphamide, and rituximab

  Samples from two-center prospective phase 2 clinical trail conducted at Ohio State University (Columbus, OH) and Mayo Clinic (Rochester, MN) were analyzed in 12 cases. All patients had progressive CLL as defined by National Cancer Institute (NCI) Working Group criteria. Patients provided written informed consent for correlative studies according to the Declaration of Helsinki on an institutional review board approved protocol for the collection and use of samples for research purposes from both participating institutions. Eligible patients received a regimen consisting of pentostatin (2 mg/m2), cyclophosphamide (600 mg/m2), and rituximab (375 mg/m2) provided intravenously on day 1 of a 21-day cycle for a maximum of 6 cycles. Responses were assessed by NCI Working Group criteria and included a bone marrow evaluation and two-color flow cytometry 2 months after completion of therapy. Peripheral blood samples from these patients collected longitudinally before (Pre-Baseline), at (Baseline) and after therapy (Relapse) were analyzed for genomic heterogeneity and clonal complexity by high throughput Exome sequencing.

  Study phs000794

• Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) - OncoArray Genotypes

  The data come from 40 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). CIMBA recruits individuals with pathogenic mutations in BRCA1 or BRCA2. The majority of carriers were recruited through cancer genetics clinics offering genetic testing, and were enrolled into national or regional studies. The remainder were identified by population-based sampling of cases, or community recruitment. Eligibility to participate is restricted to carriers of pathogenic BRCA1/2 mutations who were 18 years or older at recruitment. Information collected included amongst other variables: age at recruitment; ages at breast and ovarian cancer diagnosis; and estrogen receptor (ER) status. Samples were genotyped using the Illumina OncoArray beadchip 500K SNP custom array. Details of the genotyping process and sample selection are included in Phelan et al, Identification of twelve new susceptibility loci for different histotypes of epithelial ovarian cancer, Nat Genet. 2017 May;49(5):680-691 (PMID:28346442), and Milne et al, Identification of ten variants associated with risk of estrogen receptor negative breast cancer, Nat Genet (in press).

  Study phs001321

• Single-Cell DNA and RNA Sequencing over A 29-Year Period of A CLL Patient Demonstrating Evolution of Multiple Cell Clones

  Chronic lymphocytic leukemia (CLL) is a frequent B-cell malignancy characterized by recurrent somatic chromosome alterations and a low level of point mutations. Here we present single nucleotide polymorphism (SNP) microarray analyses of a single CLL patient over 29-years of observation and treatment, and transcriptome and whole genome sequencing at selected timepoints. We identified chromosome alterations of 13q14-, 6q- and 12q+ in early cell clones, elimination of clonal populations following therapy, and subsequent appearance of a clone containing trisomy 12 and a chromosome 10 copy neutral loss of heterogeneity (LOH) that marks a major population dominant at death. Serial single cell RNA sequencing revealed elevated mRNA expression of the FOS, JUN and KLF4 genes just before the appearance of acute disease requiring therapy. This gene expression pattern became undetectable following therapy, but reoccurred following relapse. Transcriptome evolution indicates that complex expression changes occur over time. In conclusion, CLL can evolve gradually during indolent phases, and undergo rapid changes following therapy.

  Study phs001181

• Adoptive Cell Therapy of Autologous T cell Receptor-Engineered T Cells Targeting the p53 Neoantigens in Human Solid Tumors

  Adoptive T cell therapy (ACT) using T cell receptor (TCR)-engineered T cells can mediate regression of advanced human cancers. How TCR-transduced T cells change post-ACT in vivo is currently not well understood. We screened 78 patients for TP53 mutations by whole exome sequencing and whole transcriptome sequencing. Of the 78 patient samples, 49 patient samples were previously reported (phs001003) and 2 additional patients were reported in phs002735. In this accession we report the sequencing data of remaining 23 patients and include 6 additional samples from previously published patients. Additionally, we treated a patient with chemo-refractory breast cancer that expressed mutant p53 R175H with autologous T cells transduced with an allogeneic TCR targeting p53 R175H with human leukocyte antigen (HLA) restriction of A*02. The patient had a partial response for 6 months with 55% tumor regression. The T cells from the infusion product and the peripheral blood T cells at 6 weeks post-ACT were analyzed using 10X single cell transcriptome analysis.

  Study phs002928

• Rare germline sequence variants, copy number variations, expression alterations, methylation variations and disease susceptibility in familial melanoma

  Although a number of high-risk melanoma genes have been identified, they account for melanoma risk in less than 40% of melanoma-prone families. The major goals of our research are to identify novel high-penetrance genes for familial melanoma in melanoma-prone families without known mutations and to identify modifier factors, including genetic and epigenetic, in melanoma-prone families with and without known major mutations. We are using whole exome sequencing to identify additional high-risk susceptibility genes and using targeted sequencing to follow up on the top genes in additional family members and in population-based melanoma cases and controls. We are also using RNASeq to investigate genome-wide allele-specific expression and eQTLs in these families to identify genes with expression alteration, which will provide additional information for gene discovery. In addition, we also conducted genome-wide copy number variation (CNV), methylation, and miRNA expression profiling analyses in our melanoma families with the goal of identifying genetic and epigenetic factors as disease modifiers.

  Study phs001177

• Tumor Fraction Guided Cell-Free DNA Profiling in Metastatic Cancer Patients

  Cell-free DNA (cfDNA) profiling is increasingly used to guide cancer care, yet mutations are not always identified. The ability to detect somatic mutations in plasma depends on both assay sensitivity and the fraction of circulating DNA in plasma that is tumor-derived (i.e cfDNA tumor fraction). We hypothesized that cfDNA tumor fraction could inform the interpretation of negative cfDNA results and guide the choice of subsequent assays of greater genomic breadth or depth. Plasma samples collected from 118 metastatic cancer patients were analyzed with cf-IMPACT, a modified version of the FDA-authorized MSK-IMPACT tumor test that can detect genomic alterations in 410 cancer-associated genes. Shallow whole genome sequencing (sWGS) was also performed in the same samples to estimate cfDNA tumor fraction based on genome-wide copy number alterations using z-score statistics. The results show that cfDNA tumor fraction can inform the interpretation of negative cfDNA results and guide the selection of subsequent sequencing platforms that are most likely to identify clinically-relevant genomic alterations.

  Study phs002290

• Pharmacogenetics of Efavirenz Discontinuation for Reported Central Nervous System Symptoms Appears to Differ by Race

  This is an evaluation of genetic associations with efavirenz (EFV) discontinuation for central nervous system (CNS) symptoms within 12 months of treatment. Patients were treated at an HIV primary care clinic in Nashville TN from 1998 to 2012. Previously known SNPs in CYP2B6 and CYP2A6 were used to define metabolizer genotypes (extensive, intermediate, slow metabolizer). Over 500,000 SNPs from genome-wide genotyping were used to define MDS (Multidimensional Scaling) coordinates to account for population stratification. Patients were defined as cases if they discontinued EFV for CNS symptoms within 12 months, otherwise they were defined as controls if they did not stop treatment. Among 563 evaluable participants, the hazard ratio for EFV discontinuation for CNS symptoms was 4.9 (95% C.I. 1.9 TO 12.4, p=0.001) in slow metabolizers compared to extensive metabolizers. This association was very significant in Whites 6.5 (95% CI: 2.3 to 18.8; p = 0.001), but not in Blacks 2.6 (95% C.I. 0.5 to 14.1; p = 0.27). The reason for this difference by race is not clear and warrants further investigation.

  Study phs001253