2919 results for "RNA AND sequencing OR transcriptome"

in 61.22 milliseconds.

• High-coverage whole exome sequence in non-TRU-type lung adenocarcinomas

  In order to elucidate the key genetic events in non-TRU-type lung cancer, we selected 43 non-TRU-type lung adenocarcinomas and performed a whole exome sequence and RNA-seq analysis using a next-generation sequencer. The results of the analysis identified mutations in TP53, KRAS and NKX2-1 as the top three significantly mutated genes, while the EGFR mutation was rare in this cohort. Eight NKX2-1 mutations (5 frameshift, 2 nonsense, and one missense) were identified, with one case harboring two distinct NKX2-1 mutations (missense and frameshift). All NKX2-1 mutants were transcriptionally inactive in reporter assays. Histologically, invasive mucinous adenocarcinomas accounted for most of the NKX2-1 mutations (5 cases), while 1 enteric and 1 acinar adenocarcinoma harbored the NKX2-1 mutation.

  Study JGAS000105

• Evaluation of EBUS-TBNA Aspirates from Advanced NSCLC for Comprehensive Sequencing Platforms Including Whole Genome Sequencing

  Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2–3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

  Study EGAS00001007708

• Parallel CRISPR Editing of Familial ALS Mutations Into a Healthy Control iPSC Line

  Induced pluripotent stem cells (iPSCs) provide a modeling system for human disease within the critical context of the human genome. CRISPR editing of iPSCs allows the construction of isogenic pairs where the cells have nearly identical genomes except for the desired edited change. Here, we used CRISPR editing on a healthy control line to generate three new isogenic iPSC lines that each contains a different familial amyotrophic lateral sclerosis (fALS) mutation (TARDBP[G298S/+], SOD1[G85R/+], and PFN1[G118V/+]). As controls for each mutation, we retained iPSC lines that went through the editing process but did not have the desired heterozygous mutation. We have produced high-throughput sequencing data from these isogenic lines in two different studies:(1) In “Human sensorimotor organoids derived from healthy and amyotrophic lateral sclerosis stem cells form neuromuscular junctions” (Pereira et al. 2021), we performed exome sequencing to confirm the absence of off-target CRISPR events. We also performed single-cell sequencing on sensorimotor organoids from a second healthy control iPSC line, which is deposited alongside the exome sequencing data.(2) In “iPSC Motor Neurons, But Not Other Derived Cell Types, Capture Gene Expression Changes in Postmortem Sporadic ALS Motor Neurons” (Held et al. 2023), we differentiated these isogenic lines into motor neurons, sensory neurons, cortical neurons, and astrocytes and then performed bulk RNA-sequencing. We then assessed similarities between fALS mutations, differences between cell types, and the capacity of iPSC modeling to capture differentially expressed genes in postmortem sporadic ALS motor neurons.

  Study phs002440

• Aurora US Metastatic Breast Cancer Retrospective Project

  The AURORA US Metastatic Breast Cancer project is funded by the Breast Cancer Research Foundation (BCRF) Evelyn H. Lauder Founder's Fund for Metastatic Breast Cancer Research. This multi-center effort conducted within the Translational Breast Cancer Research Consortium (TBCRC) and cancer researchers to better understand the metastatic process through the study of both the primary and metastatic tissue. In the retrospective phase, TBCRC sites submitted matched primary and metastatic tissues and blood from previous collections for piloting the process. Samples were profiled using whole genome DNA sequencing, whole exome DNA sequencing, DNA methylation arrays, and RNA sequencing. The final freeze set for samples with successful nucleic acid extraction and molecular assays included 55 patients with 31 primary tissues and 102 metastases. Twenty patients had tissue collected at autopsy and included 19 with more than one metastasis. Metastases were from 20 different tissue locations with the most common sites being liver, lung, lymph node, and brain. The median age at primary diagnosis was 49 years. The clinical subtypes of the primary tumors were 34% Triple Negative, 30% ER+HER2-, 11% ER+HER2+, and 7% ER-HER2+. Patients received an average of 3 lines of therapy in the metastatic setting. The median overall survival from primary diagnosis was 4.5 years and overall survival from metastatic diagnosis was 2 years.

  Study phs002622

• Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

  Sequencing is performed on PCR products from the amplification of gRNA segments introduced into cell lines (some are commercial, some are derived by collaborators, as described in HMDMC 13/076). Cellular DNA or RNA are not sequenced. This data is defined as "Other Biological/Biochemical Assay Data" under the WTSI datasharing guidelines and consequently will be released at the time of publication and/or through faculty websites. SEQUENCING DETAILS: in October 2016 3 DNA samples will be submitted for PCR amplification and purification through CGP pipeline using Thompson's custom primer set and PCR settings and sequencing with Hiseq using a custom protocol (with 3 forward reads to read sgRNA1, sgRNA2 and the index, optimized with Illumina Bespoke Sanger team); read length 50bp. The following set of samples will be submitted in later 2016 and ealry 2017 as cell pellet to be processed through CGP Cas9 pipeline (Thompson's oligos and protocols) and sequences on Hiseq (same custom sequencing protocol described above) We are utilizing a custom CRISPR/Cas9 library to identify synthetic lethal gene pairs. Cells will be transduced with LV expressing Cas9, then infected with LV expressign the library of 42'000 sgRNA pairs. Cell pellet will be collected from the cell population at different timepoints. The samples will be processed by CGP pipeline customized for the specific protocol and then sequenced to identify the 19bp of the sgRNA and quantify their representation.

  Study EGAS00001002117

• PsychENCODE Consortium: Epigenetic and Transcriptional Dysregulation in Autism Spectrum

  Our current understanding of autism spectrum disorders (ASD) delineates a highly heritable, yet etiologically heterogeneous disease. Forward genetic approaches to find disease associated mutations or common variation have been successful and continue to offer considerable power. Yet, given the accumulating evidence for very significant heterogeneity and environmental influences, complementary approaches to classic forward genetics become necessary. Genetic polymorphism and mutation data to date have identified dozens of causal or contributory variants, yet our preliminary data from autism brain suggest that common molecular pathways are involved in a significant subset of cases. This convergence at the tissue level suggests that other mechanisms, specifically epigenetic changes, combined with genetic background, are contributing to such final common pathways. We further tested this hypothesis by taking a comprehensive and integrative genome-wide approach to assessing brain gene-expression, miRNA levels and the related, causal epigenetic mechanisms in ASD etiology. We performed RNA-seq analyses of four cerebral cortical regions and cerebellum from ASD cases and controls, to assess mRNA, miRNA, and splicing isoform regulation. In parallel, we identified key differences in chromatin state and DNA methylation across multiple brain regions in the same ASD and control individuals used in the expression analyses using ChIP-Seq and MeDIP. We assessed the mechanisms by which changes in DNA methylation, histone modification, and DNA sequence contribute to the observed differences in gene expression. This work, which represents an unprecedented effort to unify these often disparate data (usually produced without integration in mind), delineates potential shared molecular pathways in ASD and the underlying mechanism of these differences at the level of miRNA, the chromatin regulatory apparatus, and DNA methylation. The following substudies are part of the PsychENCODE release at dbGaP and offer additional molecular data: PsychENCODE: RNA-Sequencing - SRRM4 Splicing Study phs000872 PsychENCODE: Global Changes in Patterning, Splicing and lncRNAs phs001061 PsychENCODE: Chromatin Contact Map in Fetal Cortical Laminae phs001190 PsychENCODE: Epigenetic Dysregulation in Autism Spectrum Disorder phs001220

  Study phs001022

• Single Cell RNA-Seq Reveals Malignant and Stromal Programs Associated with Metastasis in Head and Neck Cancer

  This study of head and neck cancer aimed to profile individual malignant, stromal, and immune cells in treatment-naive oral cavity squamous cell carcinoma, with or without matched metastatic cervical lymph nodes. Clinical and pathologic information was collected on all patients, in addition to a fresh biopsy sample from the primary tumor and/or metastatic lymph node at the time of initial surgery. Biopsy specimens were processed for single cell RNA-seq (scRNA-seq) and/or whole exome sequencing (WES). A total of ~6,000 single cells from 18 participants were included in the final dataset. Stromal and immune cells had consistent expression programs across patients. Conversely, malignant cells varied within and between tumors in their expression of signatures related to cell cycle, stress, hypoxia, epithelial differentiation, and partial epithelial-to-mesenchymal transition (p-EMT). One of the top scoring genes in this program was TGFβ-induced (TGFBI), implicating the classic EMT regulator TGFβ. Cells expressing the p-EMT program spatially localized to the leading edge of primary tumors. By integrating single-cell transcriptomes with bulk expression profiles for hundreds of tumors, we refined Head and Neck Squamous Cell Carcinoma (HNSCC) subtypes by their malignant and stromal composition, and established p-EMT as an independent predictor of nodal metastasis, grade, and adverse pathologic features. Our results provide insight into the HNSCC ecosystem and define stromal interactions and a p-EMT program associated with metastasis.

  Study phs001474

• Alloreactive T-cell receptor (TCR) repertoire in kidney transplantation

  We analyzed the T-cell receptor (TCR) repertoires from twelve kidney transplant recipients. Six out of the twelve kidney transplant recipients experienced a cellular rejection after kidney transplantation. TCR repertoires of CD4+ and CD8+ positive T-cells were assessed prior to transplantation and after transplantation at time of allograft biopsy using RNA based T-cell receptor beta next generation sequencing (NGS). In addition, the pre-formed alloreactive TCR repertoire for each kidney transplant recipient was identified using mixed lymphocyte reaction and donor reactive T-cells were subjected to TCR beta sequencing. In two out of the six patients with cellular rejection the TCR repertoire of graft infiltrating T-cells was additionally captured. This dataset comprises a total of 98 samples. NGS TCR beta libraries of all samples were sequenced on an Illumina NextSeq 500 and raw sequencing data (in the form of fastq files) as well assembled clonotypes and their counts (in the form of clonotype tables) are provided.

  Dataset EGAD00001007695

• Whole exome sequencing and RNA sequencing of cervical cancer

  We designed this study to identify molecular markers involved in the development and prognosis of cervical adenocarcinoma. In this study, we would like to elucidate the molecular biological characteristics of cervical adenocarcinoma in relation to the carcinogenesis of other cervical cancers (squamous cell carcinoma), vaginal and vulvar cancers. In addition, some squamous cell carcinomas of the cervix develop at a young age and have a poor prognosis. Furthermore, small cell carcinoma of the cervix also progresses rapidly and has a poor prognosis, making the establishment of diagnostic and therapeutic methods an urgent task.

  Study JGAS000582

• Whole exome sequencing and RNA sequencing of cervical cancer

  We designed this study to identify molecular markers involved in the development and prognosis of cervical adenocarcinoma. In this study, we would like to elucidate the molecular biological characteristics of cervical adenocarcinoma in relation to the carcinogenesis of other cervical cancers (squamous cell carcinoma), vaginal and vulvar cancers. In addition, some squamous cell carcinomas of the cervix develop at a young age and have a poor prognosis. Furthermore, small cell carcinoma of the cervix also progresses rapidly and has a poor prognosis, making the establishment of diagnostic and therapeutic methods an urgent task.

  Study JGAS000586

• scRNA-seq and snRNA-seq of trophoblast stem cells (TSCs) differentiation into extravillous trophoblast organoids (EVTs)

  This dataset is part of a study that aims to compare in vivo human trophoblast differentiation into EVTs to different in vitro trophoblast organoids using single-cell and single-nuclei RNA sequencing. This specific dataset includes scRNA-seq and snRNA-seq data from trophoblast stem cells (TSCs). Trophoblast stem cell (TSC) lines BTS5 and BTS11 derived by Okae and colleagues were grown as described previously (Okae et al. 2018) together with EVT differentiation media. This study shows that the main regulatory programs mediating EVT invasion in vivo are preserved in in vitro models of EVT differentiation from primary trophoblast organoids and trophoblast stem cells. Data for primary trophoblast organoids is available under E-MTAB-12650.

  Dataset EGAD00001010017

• mRNA and T cell receptor sequencing of patients with Pandemrix-associated narcolepsy type 1

  Increased risk of Narcolepsy type 1 (NT1) was reported among children and adolescents vaccinated with AS03-adjuvanted pandemic influenza A vaccine (H1N1; Pandemrix®). We hypothesized that viral T cell epitope(s) in Pandemrix mimicked self-epitope(s), and that recognition of these cross-reactive epitopes contributed to disease-specific autoimmunity. We demonstrate that pediatric, Pandemrix-associated NT1 patients had enhanced T-cell immunity against dominant T-cell epitopes of Pandemrix vaccine viral proteins neuraminidase (NA) and nucleoprotein (NP), and also responded against a self-epitope in brain-expressed protein-O-mannosyltransferase 1 (POMT1). Here, the gene expression profiles and TCR repertoire of PBMC from patients and healthy controls were compared. PBMC samples were stimulated with peptides NA-175-189 or POMT1-675-589, and bulk RNA isolated and sequenced (RNAseq and TCRseq).

  Study EGAS00001004886

• Multi-omics integration reveals only minor long-term molecular and functional sequelae in immune cells of individuals recovered from COVID-19

  The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in convalescent COVID-19 and healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays moderate changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1RA and IFN?? was higher in convalescent individuals. In conclusion, despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

  Study EGAS00001005529

• Exploring the Microbiome-Gut-Brain Axis with Respect to Psychoneurological Symptoms for Children with Solid Tumors

  The human body hosts tens of trillions of microbes, a number 10 times greater than the number of human body cells. The collection of these human-associated microbes and their genomes is called the “human microbiome”. More than 90% of these microbes live in the gastrointestinal (GI) tract, representing 500 species on average. With the development of high-throughput sequencing technology and bioinformatics tools (e.g., 16S ribosomal RNA or 16S rRNA sequencing and shotgun metagenomics analysis), recent progress has been made to elaborate the critical role that the gut microbiome plays in human health and disease. A significant disruption (dysbiosis) of the composition and function of the gut microbiome is associated with carcinogenesis, chemotherapeutic metabolism, and treatment-related symptom toxicities. The current state of science on the gut microbiome in cancer has only recently started in animal models and adults, and has not been well investigated in children diagnosed and treated for cancer, nor with respect to associations with cancer treatment-related toxicities. This proposed study will explore the hypothesis that dysbiosis of the gut microbiome is associated with cancer treatment-related GI and psychoneurological symptoms. We will collect both stool samples (for the gut microbiome) and patient-reported outcomes across chemotherapy.

  Study phs002960

• Ewing and Ewing-like-sarcoma tumoroid biobank reveals distinct drug sensitivities: translocation makes the difference

  Ewing sarcoma (ES) and Ewing-like sarcomas (ELS) are highly aggressive tumors found in soft tissues and bone of mostly children and young adults. Despite being different in many aspects, including genetics, possible cell-of-origin, and pathology, patients with any of these entities all receive the same therapeutic regimen. Currently, there is a lack of preclinical cell models that can predict drug response of these entities in vitro. Here, we describe a pediatric ES and ELS tumor organoid (tumoroid) biobank containing tumoroids with different translocations, including EWSR1::FLI1, EWSR1::ERG, CIC::DUX4, and BCOR-rearrangements. Using histology, whole genome sequencing and RNA sequencing, we demonstrate that these tumoroids retain key properties of their matching patient tumors. In addition, we traced tumor subclones in the tumoroids across patient-matched longitudinal samples, which showed that cellular heterogeneity was maintained. We subsequently performed drug screening on the tumoroid models, which unveiled entity-specific drug sensitivity to various cytotoxic compounds and targeted compounds, including MCL-1 inhibitors for CIC::DUX4 sarcomas. Taken together, this newly established ES and ELS patient-derived tumoroid biobank generates a unique source of material for future basic cancer research and drug-screening.

  Study EGAS00001007941

• RNA sequencing of undifferentiated sarcomas

  Undifferentiated sarcomas (USARCs) of adults are diverse, rare, and aggressive soft tissue cancers. Recent sequencing efforts have confirmed that USARCs exhibit one of the highest burdens of structural aberrations across human cancer. Here, we sought to unravel the molecular basis of the structural complexity in USARCs by integrating DNA sequencing, ploidy analysis, gene expression, and methylation profiling. We identified whole genome duplication as a prevalent and pernicious force in USARC tumorigenesis. Using mathematical deconvolution strategies to unravel the complex copy-number profiles and mutational timing models we infer distinct evolutionary pathways of these rare cancers. In addition, 15% of tumors exhibited raised mutational burdens that correlated with gene expression signatures of immune infiltration, and good prognosis.

  Study EGAS00001003291

• MicroRNAs, Hypertension and End Organ Damage in Humans

  Study 1) Profiling of microRNAs in plasma of patients with hypertension complicated or not by metabolic syndrome or chronic kidney disease Hypertension (HTN) or high blood pressure is associated with subclinical target organ damage such as cardiac, vascular and kidney injury, which may lead to complications and death. In this study, we investigated circulating microRNAs as biomarkers of target organ damage in HTN patients. We profiled circulating microRNAs by RNA sequencing in platelet-poor plasma of normotensive subjects and patients with HTN complicated or not by metabolic syndrome (MetS) or chronic kidney disease (CKD) (n=15 each group). Differentially expressed microRNAs were identified with a threshold of false discovery rate <0.1. Differentially expressed microRNAs were identified uniquely to associate with HTN (8), MetS (1) or CKD (13), and 8 were similarly differentially expressed in different groups. This study identified the association of differentially expressed circulating microRNAs with target organ damage in HTN patients, which could have some pathophysiological and therapeutic implications. Study 2) Profiling of microRNAs in gluteal subcutaneous small arteries of patients with hypertension complicated or not by chronic kidney disease Hypertension (HTN) is associated with vascular damage characterized by endothelial dysfunction and vascular remodeling and stiffening, which contributes to kidney injury leading to chronic kidney disease (CKD). MicroRNAs are short non-coding RNAs which repress/degrade target mRNAs. The microRNA role in vascular injury in HTN remains unclear. In this study, we aimed to identify differentially expressed microRNAs in small arteries of patients with HTN associated or not with CKD, in order to shed light on the pathophysiological molecular mechanisms. Normotensive subjects and HTN patients associated or not with CKD grades 3-4 were studied (n=15-16). Small arteries were isolated from subcutaneous gluteal biopsies, RNAs were extracted, and small and total RNA sequencing was performed by Illumina HiSeq-2500. Differentially expressed genes were identified with a P<0.05 and fold change (FC) >1.3. Differentially expressed microRNAs and mRNAs were identified uniquely to associate with HTN (microRNAs: 10, mRNAs: 68), CKD (microRNAs: 68, mRNAs: 395), and both groups (microRNAs: 2, mRNAs: 32). This study identified differentially expressed microRNAs and mRNAs in small arteries with target organ damage in HTN, which could have some pathophysiological and therapeutic implications.

  Study phs002389

• Autism Post-Mortem Brain RNA-Sequencing: SRRM4 Splicing Study

  Purpose: The goal of this study was to assess the status of splicing changes in microexons in the cortex of individuals with autism. Methods: We performed RiboZero Gold (rRNA depleted) 50bp PE RNA-seq in a larger set of case and control samples to define 12 autism and 12 control samples showing the greatest global differential gene expression change. These samples, which show differential expression of the splicing regulator SRRM4, were used to evaluate global splicing changes. Results: Within these samples, 126 of 504 (30%) detected alternative microexons display, a mean ΔPSI > 10 between ASD and control subjects of which 113 (90%) also display neural-differential regulation. By contrast, only 825 of 15,405 (5.4%) longer (i.e. > 27 nt) exons show such misregulation, of which 285 (35%) correspond to neural-regulated exons. Notably, we also observe significantly higher correlations between microexon inclusion and nSR100 mRNA expression levels across the stratified ASD samples and controls, for those microexons regulated by nSR100 relative to those microexons that are not regulated by this factor (p=1.4x10-7, Wilcoxon Sum Rank test). Conclusions: These data suggest microexon regulation is a potentially important mechanism underlying ASD and likely other neurodevelopmental disorders.

  Study phs000872

• Deep single-cell RNA sequencing data for 12346 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive NSCLC patients

  Cancer immunotherapies have shown sustained clinical responses in treating non-small cell lung cancer (NSCLC), but efficacy varies between patients and is believed to depend in part on the amount and properties of tumor infiltrating lymphocytes (TILs). To comprehensively depict and dissect the baseline landscape of the composition, lineage and functional states of TILs in lung cancer, here we generated deep single-cell RNA sequencing data for 12,346 T cells from the primary tumour, adjacent normal tissues and peripheral blood of 14 treatment-naive NSCLC patients. Combined expression and TCR-based lineage tracking revealed a significant proportion of effector T cells with common origins and similar functional states across peripheral blood and tumours pointing towards a highly migratory nature of these T cells. We also observed tumour-infiltrating CD8+ T cells undergoing extensive clonal expansion and exhaustion, with two clusters of cells exhibiting states preceding exhaustion. Survival analysis on independent datasets suggested that high ratio of pre-exhausted to exhausted T cells is associated with better prognosis of lung adenocarcinoma (LUAD). In addition, we observed further heterogeneity within the tumour regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of this group of activated tumour Tregs, which included IL1R2, correlated with poor prognosis in LUAD. The T cell clusters revealed by our single cell analyses provide a new approach for patient stratification, and the accompanying compendium of data will help the research community to gain further insight into the functional states and dynamics of T cell responses in lung cancer.

  Study EGAS00001002430

• Altered oligodendrocyte heterogeneity in Multiple sclerosis revealed by single nuclei RNA sequencing

  This study is currently hosted by the European Nucleotide Archive. To access the data contained within the Study please follow the link below: https://www.ebi.ac.uk/ena/browser/view/PRJEB39323 Oligodendrocyte (OL) pathology is increasingly implicated in neurodegenerative diseases, as they are involved in metabolic support of axons and functional cross-talk with other brain cells. Rodent OLs are heterogeneous, with developmental and biological differences, but the extent of heterogeneity in the normal human brain and its contribution to any changes to disease remains unknown. Here we performed single nuclei RNA-sequencing (snRNA-seq) from white matter (WM) areas of post mortem human brain both in control (Ctr) and multiple sclerosis (MS) patients. We identified several sub-clusters of oligodendroglia in the Ctr human WM, some similar to those in mouse, and defined new markers for these cell states. Strikingly, some of these sub-clusters were under-represented in MS tissue, while others were more prevalent than in controls. We found a lack of OL precursor cells (OPCs) and an OL subcluster in an intermediate stage of differentiation in MS lesions and in normal appearing white matter (NAWM), suggesting either depletion by the disease or by a regenerative response. The differences in mature OL sub-clusters indicate different functional states of OLs in MS tissue and, as this is similar in NAWM to lesions, that MS is a more diffuse brain disease than the focal demyelinating lesions suggest. We were also able to identify new putative markers of different MS lesion subtypes. Our findings of an altered heterogeneity of oligodendroglia in MS may have an important contribution to our understanding of disease progression and may alter therapeutic approaches to MS.

  Study EGAS00001003412

• HCA_Female_Reporductive_Adult_WSSS_RNA

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

  Study EGAS00001004210

• Upper cortical layer-driven network impairment in schizophrenia

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

  Study EGAS00001006495

• HCA_Placenta_Adult_Vento_RNA

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

  Study EGAS00001004069

• Engineering of a promoterless anti-viral RNAi hairpin into an endogenous miRNA locus

  Short hairpin RNA (shRNA) expression strategies that allow safe and persistent target mRNA knockdown are key to the success of many in vitro or in vivo RNAi applications. Here, we propose a novel solution which is expression of a promoterless miRNA-adapted shRNA (shmiRNA) from an engineered genomic miRNA locus. For proof-of-concept, we genetically “vaccinated” liver cells against a human pathogen, by using TALEns or CRISPR to integrate an anti-hepatitis C virus (HCV) shmiRNA into the liver-specific miR-122/hcr gene. Reporter assays and qRT-PCR confirmed anti-HCV shmiRNA expression as well as miR-122 integrity and functionality. Specificity and safety of shmiRNA integration were validated via PCR, cDNA and miRNA profiling, and whole genome sequencing. A subgenomic HCV replicon and a full-length reporter virus, but not a Dengue virus control, were significantly impaired in the modified cells. Our original combination of DNA engineering and RNAi expression technologies should benefit numerous applications, from basic miRNA research, to human cell and gene therapy.

  Study EGAS00001001252

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

  Despite the established role of the transcription factor MYC in cancer, little is known about the involvement of lncRNAs in mediating MYC’s function. Here we have intersected RNA-sequencing data from MYC-inducible cell lines, from a cohort of 91 mature B-cell lymphomas and from sorted germinal-center B-cells. By this approach, we identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated by MYC in the model cell lines. Among them we focused on a lncRNA that we named MINCR, showing a strong correlation with MYC expression in MYC-positive lymphomas and in pancreatic ductal adenocarcinomas. RNAi experiments showed that MINCR controls cellular viability by influencing the expression of MYC-regulated cell cycle genes. Finally we provide evidences that down-regulation of AURKA, AURKB and CTD1 may explain the reduction in cellular proliferation observed upon MINCR knock-down. We therefore suggest that MINCR acts as a modulator of MYC transcriptional control of cell cycle genes.

  Study EGAS00001001199

• CD36 defines CML cells less sensitive to imatinib

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

  Study EGAS00001002421

• Mesenchymal inflammation drives genotoxic stress in hematopoietic stem cells and predicts disease evolution in human pre-leukemia

  Mesenchymal niche cells may drive tissue failure and malignant transformation in the hematopoietic system but the molecular mechanisms and their relevance to human disease remain poorly defined. We demonstrated that perturbation of mesenchymal cells in a mouse model of the preleukemic disorder Shwachman-Diamond syndrome induces mitochondrial dysfunction, oxidative stress and activation of DNA damage responses in hematopoietic stem and progenitor cells. In this study we demonstrate, through massive parallel RNA sequencing of highly purified mesenchymal cells in a range of human preleukemic syndromes, TP53-S100A8/9-TLR4 inflammatory signaling as a common driving mechanism of genotoxic stress, which could be attenuated by TLR4 blockade. S100A8/9 expression in mesenchymal cells predicted outcome in myelodysplastic syndromes, the principal human preleukemic condition, independent of known prognostic variables. Collectively, findings reveal a concept of mesenchymal niche-induced genotoxic stress in heterotypic stem and progenitor cells through inflammatory signaling as an actionable determinant of disease outcome in human preleukemia. The data further provide novel conceptual and mechanistic insights into the intimate link between inflammation and cancer.

  Study EGAS00001001926

• Gene fusion and transcriptomic landscapes of sarcomas

  Sarcoma represents a highly heterogeneous group of tumors. We report here the first unbiased and systematic search for gene fusions and analysis of transcriptomic profiles in 100 monomorphic sarcomas cases using RNA-sequencing profiling. Fusion genes were detected in two thirds of samples. Similarly to fusion genes such as PAX3/7-FOXO1, EWSR1/FUS-ETS, SS18-SSX, BRD3/4-NUTM1 or EWSR1/FUS/TAF15-NR4A3 characterizing well-defined entities, we show that all sarcomas displaying CIC fusions- whatever their fusion partner (being DUX4, NUTM1 or FOXO4) - form a transcriptionally homogeneous group of tumors. Likewise, tumors with either a BCOR fusion (with CCNB3, MAML3 or ZC3H7B partners) or with a BCOR internal duplication form a single biological entity. We also found that EML4-ALK fusions also characterize a homogeneous infantile fibrosarcoma subgroup. In contrast, fusions like VGLL2-NCOA2/CITED2, TMP3/TRP-NTRK1 or ETV6-NTRK3 are observed in more heterogeneous tumors. Finally, we also describe a new group of bone sarcomas characterized by EWSR1- or FUS-TFCP2 fusions.

  Study EGAS00001002189

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

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

  Study EGAS00001002895

• “Castration-persistence” is a distinct state of tolerance to androgen receptor targeting therapies in prostate cancer

  Androgen receptor (AR) signalling is important in prostate cancer progression, and therapies that specifically target this pathway are the mainstay of treatment for advanced disease. Treatment however is non-curative, and resistance develops inevitably with time. Although the mechanisms that drive progressive ‘castration resistant’ disease are reasonably well characterized, how tumours survive acute pathway inhibition is unclear. We performed a neo-adjuvant study of a novel combination of AR targeting therapies, and noted that the objective response to treatment was highly variable. To determine what was driving tumour persistence in poorly responding patients, we comprehensively characterised pre- and post-treatment samples were using both whole genome and RNA-sequencing. We find that ‘castration-persistence’ is a distinct state from ‘castration-resistance’, and is mediated by global transcriptional reprogramming leading to transitional EMT state, which is shared with benign luminal epithelial cells. This appears to be AP-1 and KLF driven, and represents and integration of multiple signalling pathways, particularly IGF2, offering a number of tractable strategies to improve clinical response to AR targeting therapies.

  Study EGAS00001003172

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

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

  Study EGAS00001005786

• Immuno-genomic Profiling of Biopsy Specimens Predicts Neoadjuvant Chemotherapy Response in Esophageal Squamous Cell Carcinoma

  Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We introduce a model by machine learning and validated its potential to predict the NAC response in ESCC. To create this model, we performed whole genome sequencing and RNA-seq analysis on totally 141 ESCC biopsy specimens before NAC treatment and analyzed their immuno-genomic features and association with NAC response. Neutrophils infiltration could play an important role in ESCC response to NAC. We also demonstrated that specific copy number alterations and copy number signatures in the ESCC genome were significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a new indicator of therapeutic capability and a new therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

  Study JGAS000535

• Genomic Analyses Identify Recurrent MEF2D Fusions in Acute Lymphoblastic Leukemia

  Chromosomal rearrangements are initiating events in acute lymphoblastic leukemia (ALL). Here, using RNA-sequencing (RNAseq) of 560 ALL cases, we identify rearrangements between MEF2D (myocyte enhancer factor 2D) and five genes (BCL9, CSF1R, DAZAP1, HNRNPUL1 and SS18) in 22 B progenitor ALL (B-ALL) cases with a distinct gene expression profile, the most common of which is MEF2D-BCL9. Examination of an extended cohort of 1164 B-ALL cases identified 30 cases with MEF2D rearrangements which include an additional fusion partner, FOXJ2, thus MEF2D-rearranged cases comprise 5.3% of cases lacking recurring alterations. MEF2D-rearranged ALL is characterized by a distinct immunophenotype, DNA copy number alterations at the rearrangement sites, older diagnosis age, and poor outcome. The rearrangements result in enhanced MEF2D transcriptional activity, lymphoid transformation, activation of HDAC9 expression and sensitive to histone deacetylase inhibitors treatment. Thus, MEF2D-rearranged ALL represents a distinct form of high-risk leukemia, for which new therapeutic approaches should be considered.

  Study EGAS00001001952

• Immunogenomics of colorectal cancer response to immune checkpoint blockade

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

  Study EGAS00001004438

• Genomic Analysis of Fibrolamellar Hepatocellular Carcinoma

  Fibrolamellar hepatocellular carcinoma (FL-HCC) is a very rare and distinct subtype of hepatocellular carcinoma that occurs in the previously healthy livers of adolescents and young adults. The low incidence rate has kept FL-HCC an understudied cancer, thus the list of FL-HCC markers is limited and the genetic alterations driving its growth have remained unknown. In order to gain a better understanding of the molecular mechanisms underlying FL-HCC initiation and progression, we performed an in-depth genomic analyses of one tumor followed by immunohistochemistry validation on seven other tumors. We showed the expression of neuroendocrine markers in FL-HCC. In addition, DNA/RNA sequencing revealed that common cancer pathways are not visibly altered in FL-HCC but identified two novel structural variants, both resulting in fusion transcripts. Through in vitro studies, we demonstrated the oncogenic properties of these fusion transcripts. These experiments further highlight the tumorigenic role of gene fusions in the etiology of pediatric solid tumors and identify both candidate biomarkers and possible therapeutic targets for this lethal pediatric disease.

  Study phs000828

• Tissue and Fluid Analysis in Ocular Inflammatory Disease

  Granulomas are found in multiple infectious and non-infectious human diseases, yet the pathogenic drivers of granulomatous disease are still poorly understood. To enrich for pathologic inflammatory cells, we leveraged the intrinsic immunologic isolation of the eye by examining granulomatous uveitis, a form of ocular inflammation. We performed single-cell RNA-sequencing to obtain an unbiased gene expression survey of ocular immune cells in the anterior chamber of four live human donors with granulomatous uveitis. We paired this with antigen-receptor sequence analysis to identify clonally expanded lymphocytes. Each individual had robust clonal expansion arising from a single T or B cell lineage, suggesting distinct pathogenic processes in each patient. This variability in clonal expansion was mirrored by individualized combinations of TH1, TH1/TH17, and Tregs in ocular CD4 T cells. In contrast, ocular CD8 T cells were more transcriptionally similar across donors. Finally, clonally expanded, activated B cells displayed evidence of plasmablast differentiation and class-switching within the ocular microenvironment. Collectively, our study identified conserved and individualized features of granulomatous inflammation, illuminating parallel pathophysiologic mechanisms.

  Study phs002449

• Epigenetic Landscape of Human Parathyroids

  The goal of this study is to produce a comprehensive map of the chromatin landscape of human parathyroids. Data from 20 subjects with parathyroid adenoma causing primary hyperparathyroidism were used. For the sequencing platform, Illumina HiSeq 2000 with V3 high-output kits were used. For the data processing pipeline, reads were aligned to hg19 using Bowtie or BWA. For RNA-seq data, transcript levels were estimated per gene to RPKM using Cufflink. For ChIP-seq and ATAC-seq data, enrichment peaks were estimated compared to inputs or background using MACS2 after merging reads from replicates. For DNase-seq, significant peaks compared to background distribution were computed using Hotspot. For Hi-C data, significant interactions were detected using Juicer HiCCUPS and 5 kb resolution for the merged data from replicates. Chromatin states for parathyroids were computed using ChromHMM based on ChIP-seq profiles in parathyroids of H3K4me3, H3K4me1, H3K27ac, H3K36me3, H3K27me3 and H3K9me3 and the emission matrix from 18 states pre-trained in the Roadmap Epigenomics. Gene expression, binding positions of transcription factor, open chromatin regions, 3D-chromatin interaction, and chromatin states are submitted.

  Study phs003302

• Whole-genome sequence data for The interface of malignant and immunologic clonal dynamics in high-grade serous ovarian cancer

  Each tumor sample was cut into three pieces, yielding two end-pieces for cryovials and a middle portion placed in 10% buffered formalin. End pieces were homogenized manually and with a paddle blender (Stomacher). All paraffin-embedded blocks, including formalin-fixed tumor samples and molecular-fixed fallopian tubes, were sectioned and stained with hematoxylin and eosin prior to expert histopathological review to confirm the presence of high grade serous carcinoma. Homogenized end pieces were then flash frozen and later used for WGS. For all tumor and matched normal (peripheral blood) samples, DNA was extracted with the Qiagen AllPrep DNA/RNA kit (tumor samples from patients 25,26,28-32) or the Qiagen Blood and Tissue Extraction Kit (tumor samples from patients 1-4,7,9-17, and all blood samples). For all tumor and normal samples, DNA extraction was followed by library construction and sequencing using Illumina HiSeq2500 whole genome shotgun v4 chemistry with paired-end 125bp reads.

  Dataset EGAD00001003984

• Differential Presence of Exons in Cell-Free DNA Reveals Different Patterns in Colorectal Cancer Between Metastatic and Non-Metastatic Patients

  Tumor genotyping during disease follow-up is gaining an increasing interest for monitoring and clinical management in colorectal cancer (CRC). Growing evidences suggest that single biopsies are inefficient for molecular profiling, due to clonal evolution and intratumoral heterogeneity of primary tumors and/or metastasis. Circulating free DNA (cfDNA) has emerged as an alternative source of genetic material that seems to be representative of the continuously changing tumor molecular features. Thus, the so-called 'liquid biopsy' also offers the additional advantages of accessibility and non-invasiveness for the patient. Novel strategies based on next-generation sequencing (NGS) of cfDNA are being currently developed, with most studies focusing on targeted deep sequencing panels of potential clinically actionable genes. Here we describe an alternative approach, consisting of exome-sequencing in cfDNA at a relatively shallower depth to provide a more general overview of the circulating genome. In this study, two groups of colorectal cancer patients, with either disseminated (M) or localized disease (N), were compared using whole-exome sequencing and a RNA-seq bioinformatics analysis pipeline. Thus, a set of 379 exons present at different levels in cfDNA from both groups were identified, giving rise to a completely new concept in this field, termed as 'Differential Presence of Exons'. Differentially present exons (DPEs) were subsequently used for M and N patients clustering and classification, allowing for the design of a predictive algorithm with encouraging preliminary results in a small subset of patients that could not be initially classified attending to the selection criteria (U group, unassignable).

  Study EGAS00001002687

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

• Genomic Profiling of Peripheral T-cell Lymphoma

  Peripheral T-cell lymphomas (PTCLs) are malignant and highly aggressive hematologic tumors arising from mature post-thymic T cells. The diagnosis of PTCL includes diverse lymphoma subgroups, altogether accounting for about 15% of all non-Hodgkin lymphomas. Despite recent advances in uncovering the mutational landscape of PTCLs, the prognoses of most PTCLs remain poor. To overcome the limits of exome sequencing and identify new genomic abnormalities in PTCL, we performed RNAseq and whole genome sequencing. From our analyses, we identified a recurrent in-frame deletion (VAV1 Δ778-786) generated by a focal deletion-driven alternative splicing mechanism as well as novel VAV1 gene fusions (VAV1-THAP4, VAV1-MYO1F, and VAV1-S100A7). These genetic lesions result in increased activation of VAV1 catalytic-dependent (MAPK, JNK) and non-catalytic-dependent (nuclear factor of activated T cells, NFAT) VAV1 effector pathways. These results support an oncogenic role for VAV1 signaling in the pathogenesis of PTCL (Abate et al., 2017). Recently, we have also identified a novel chimeric oncogene, FYN-TRAF3IP2, in the RNA-seq and whole genome sequencing analyses. FYN-TRAF3IP2 leads to aberrant NF-κβ activity downstream of T-cell receptor signaling and induces NF-κβ-driven T-cell transformation when expressed in hematopoietic progenitors in mice. This novel gene fusion supports the oncogenic driver role of NF-κβ signaling in the pathogenesis of PTCL (Kim et al., in review).

  Study phs001962

• Analysis of Donor Pancreata Defines the Transcriptomic Signature and Microenvironment of Early Neoplastic Pancreatic Lesions

  The adult healthy human pancreas has been poorly studied given lack of indication to obtain tissue from the pancreas in the absence of disease and rapid postmortem degradation. We obtained pancreata from brain dead donors thus avoiding any warm ischemia time. We found that neoplastic lesions occur frequently in healthy organs, in donors as young as in their 3rd decade of life. Using the 10X Genomics Platform, we performed single cell RNA sequencing on 7 donor pancreata, 6 of which we separately sequenced tissue from the pancreas head and tail, for a total of 11 sequencing runs. Of note, because one of our runs (containing sequencing from pancreas head and tail) did not pass qc, our final analysis in our study utilized 6 donor pancreata, 5 of which contained pancreas head and tail sequenced separately. This is a robust single-cell dataset of healthy, nonpathologic pancreas tissue and includes acinar, ductal, and non-epithelial cells (myeloid cells, fibroblasts, endothelial cells, lymphocytes).Briefly, the donor pancreas was dissected by a clinically-trained pancreatobiliary surgeon. The organ was dissected into head, body, and tail, with portions of each placed into DMEM media with 1% BSA/10µM Y27632. Tissue was minced into 1mm3 pieces, then digested with 1 mg/mL collagenase P for 20-30 min at 37°C with gentle agitation. Digested tissue was rinsed three times with DMEM/1% BSA/10µM Y27632, then filtered through a 40µm mesh. Resulting cells were submitted to the University of Michigan Advanced Genomics Core for single cell sequencing using the 10x Genomics Platform. FASTQ files of all sequencing runs along with clinical annotation of donors is available through this dbGaP submission.

  Study phs003229

• The Transcriptional Landscape of SHH Medulloblastoma

  To further understand the biology of Sonic hedgehog medulloblastoma and its molecular subtypes, we studied 250 human Shh-MB using strand-specific RNA sequencing. We identified novel alterations within the cAMP dependent pathway and found that 18% of tumors have genetic events that directly target the abundance and/or stability of MYCN. We also discovered an extensive network of fusions in focally amplified regions, and several loss-of-function fusions in tumor suppressor genes PTCH, SUFU and NCOR1. Molecular convergence on a core of specific genes by nucleotide variants, copy number aberrations, and gene fusions highlights key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma.

  Dataset EGAD00001006305

• Genomic characterisation of SDH deficient renal cell carcinoma - WGS

  SDH deficient renal cell carcinomas are a rare and recently defined subtype of kidney cancer, often associated with an inherited mutation in one of the SDH gene subunits. This dataset sought to understand the genomic events that underpin tumour formation, from putative cell of origin, characterisation of the tumour microenvironment, to the genomic evolution of these rare tumours. We performed whole genome and RNA sequencing of 4 patients with SDH deficient renal cell carcinomas, including one patient who had an additional paraganglioma. An addition patient in this cohort had the initial diagnosis revised to a clear cell renal cell carcinoma.

  Dataset EGAD00001008469

• Ongoing Replication Stress Tolerance and Clonal T Cell Responses Distinguish Liver and Lung Recurrence and Patient Outcomes in Pancreatic Ductal Adenocarcinoma

  Patients with metastatic pancreatic ductal adenocarcinoma (PDAC) survive longer when disease spreads to the lung but not to the liver. We generated overlapping, multi-omic datasets to identify molecular and cellular features that distinguish patients whose disease develops liver metastasis (liver cohort) from those whose disease progression results in lung metastasis without liver metastases (lung cohort). Lung cohort patients generally survived longer than liver cohort patients, independent of tumor subtype. We developed a pORG gene signature that distinguishes primary tumors in the liver and lung cohorts. We identified ongoing replication stress (RS) response pathways in high pORG/liver cohort tumors, while low pORG/lung cohort tumors had greater densities of lymphocytes and shared T cell clonal responses. Our study demonstrates that liver-avid PDAC is associated with tolerance to ongoing RS, limited tumor immunity, and less favorable outcomes; whereas low RS, lung-avid/liver-averse tumors are associated with active tumor immunity that may account for favorable outcomes. As expected, we found high frequencies of KRAS, TP53, CDKN2A, and SMAD4 gene alterations in our tumor samples. High pORG activity (GSVA scores) in our primary tumors appear to be positively correlated with alteration frequencies in both TP53 and CDKN2A. We did not see a correlation with alterations in KRAS as almost all the tumor samples have KRAS alterations. From a de-identified dataset of 1,873 patients diagnosed with and/or treated for PDAC at our institution between 2004 and 2020, we identified 422 patients for which we had specimens with sequencing data (N=374) and/or specific evidence of disease metastasis site(s) from the OHSU cancer registry or disease-relevant computed tomography (CT) scans to allow cohort classification. Note that our study includes RNA-Seq, DNA-Seq panel, and TCR-Seq (T-cell Receptor Sequencing) data, but that this dbGaP submission just includes the patients/samples with RNA-Seq and/or DNA-Seq panel data. Therefore, this submission includes 290 samples from 278 patients. TCR sequence data is available on the Adaptive Biotechnologies platform. Clinical course timepoints, patient demographics, stage, grade, nodal involvement, resection margins, and angiolymphatic invasion were provided as deidentified data by the OHSU cancer registry with quality control data verification in a subset by pathologists (BB and TM). We reviewed all available computed tomography (CT) scans for all patients with primary tumor resection dates recorded by the cancer registrar, with tumor samples analyzed by RNA-seq, DNA-seq, or TCR-seq, and/or with additional information indicating metastatic spread (e.g., metastatic samples received for related studies). We abstracted the site of all lesions proven to be metastatic by biopsy and/or that clearly increased in size during progression or decreased in size during treatment as long as a radiologist described the lesion as “likely”, “suspicious for”, “concerning for”, or “favor” metastasis. Clinical imaging was reviewed by a radiologist (AG) to validate patient assignments to the liver, lung, and neither liver nor lung (other recurrence site) cohorts. Time to recurrence was calculated from the earliest of either the recurrence date provided by the OHSU cancer registry, or the date of earliest lesion abstracted from CT reports. The subject attributes reported in this submission were up to date at the time of our data freeze, which occurred July 2021.

  Study phs003597

• Using Genomics to Reduce Breast Cancer Disparities in the African Diaspora

  Differences in breast cancer incidence and mortality rates between North American Caucasian and African American women are well-described and transcend socioeconomic issues. Black women are diagnosed with breast cancer at a younger median age; have more clinically aggressive disease and stage-for-stage; and have higher mortality rates than age-matched Caucasian women. Black women in West Africa, the origin of the slave trade in the US in the 19th century and thus the founder population for most African Americans, have even higher rates of early-onset, poor-prognosis breast cancer than African American women. Racial difference in the distribution of intrinsic molecular subtypes has been well characterized in the US and throughout the African Diaspora as well. Despite the large efforts on characterizing racial/ethnic differences, however, the reasons women of African ancestry are disproportionately affected by breast cancer incidence and mortality remain poorly understood - largely due to paucity of data on inherent genomic differences that contribute to the disparities in incidence and progression of breast cancer across populations. West Africa Breast Cancer Study (WABCS) is an initiative that aims to comprehensively understand the genetic architecture of breast cancer in West Africans, the founder population of a large proportion of black women in the United States. The objective of the study was to provide a better understanding of the molecular genetic factors that influence prognosis in Nigerian breast cancer patients, and determine which of these alterations may be amenable to available therapy. To that end, we examined the molecular features of breast cancers of indigenous African women using a combination of whole-genome, whole-exome, and transcriptome sequencing (WGS, WES, and RNA-seq) on 194 tumors from Nigerian patients. The goal of this project was to obtain answers to two related research questions using an unscreened population without genetic admixture in Nigeria: 1) why are women of African ancestry more likely to develop aggressive young onset breast cancer? 2) What are the associated genomic and non-genomic risk factors? We hypothesize that the genomic determinants of breast cancer molecular subtypes in women of African ancestry are also molecular drivers of tumor progression and represent targets for interventions to improve clinical outcomes and close the mortality gap. By identifying causal links between genetic variants that promote aggressive tumor progression in Nigerian women in comparison to women from different population found in TCGA and ICGC, the present dataset will have significant public health impact on millions of women in the African Diaspora. The potential to identify novel pathways for interventions to reduce the increasing mortality gap between women of African and European ancestry is huge

  Study phs001687

• RNA sequencing of a glioblastoma PDX cohort

  Patient-derived xenograft (PDX) models represent an essential tool for the preclinical evaluation of novel targeted therapies and their combinations. In this project, we molecularly characterized a comprehensive panel of glioma PDX models with well conserved molecular and pathological features over multiple passages. PDX samples were subjected to next generation sequencing (NGS). Characterization of cancer-relevant features including driver mutations or cellular processes was performed using mutational and gene expression data in order to identify potential biomarker or treatment targets in glioblastoma. In summary, we reported a newly established and molecularly characterized panel of glioma PDX models with high relevance for translational preclinical research.

  Study EGAS00001007119

• Evolution and clinical impact of genetic epistasis within EGFR-mutant lung cancers: multi-timepoint exome sequencing of a single patient's disease

  This dataset is a time-series of EGFR-mutant NSCLC clinical specimens from an individual patient profiled using tumor-based whole exome sequencing and the data is in BAM format. DNA was extracted from FFPE for primary tumor and frozen tumor tissue samples and matched non-tumor tissue using the Qiagen Allprep DNA/RNA Mini Kit.  The library preparation protocol was based on the Agilent SureSelect Library Prep and Capture System. DNA was resuspended in a low TE buffer and sheared (Duty Cycle 5%; Intensity 175; Cycles/Burst: 200; Time: 300s, Corvaris S2 Utrasonicator).  Bar-coded exome libraries were prepared using the Agilent Sure Select V5 library kit per manfucaturer’s specifications. The libraries were run on the HiSeq2500. Raw paired end reads (100bp) in FastQ format generated by the Illumina pipeline were aligned to the full hg19 genomic assembly obtained from USCS, gencode 14, using bwa version 0.7.12. Picard tools version 1.117 was used to sort, remove duplicate reads and generate QC statistics. Tumor DNA was sequenced to median depth of 303X (range 114.39-383.41) and the matched germline DNA to average depth of 231.65.

  Dataset EGAD00001003769

• RNA-seq analysis of transcriptome variation with human ESC subclones

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny 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/ .

  Dataset EGAD00001004266

• Exploratory novel biomarker and resistance mechanism of milademetan, an MDM2 inhibitor, in MDM2 amplified intimal sarcoma from an open-label phase 1b/2 trial ���NCCH1806/MK004���

  Amplified Murine double minute 2 (MDM2) is found in > 70% of intimal sarcoma, known as one of the ultra-rare sarcomas. Milademetan (DS3032, RAIN-32) is a novel, specific, small-molecule MDM2 inhibitor that disrupts MDM2 and the tumor suppressor protein p53 interactions in tumor cells. We conducted a phase 1b/2 trial (Trial registration No: JMA-IIA00402) in patients with amplified MDM2 wild-type TP53 intimal sarcoma as a sub-study under the nationwide large registry for rare cancers in Japan (MASTERKEY Project). Eleven patients were enrolled, and ten were included in the efficacy analysis. Two (20%) patients had durable responses for > 15 months. Milademetan provided clinical benefits in patients with amplified MDM2 intimal sarcoma. Predictive biomarkers other than amplified MDM2 and acquired resistance mechanisms for milademetan are unknown. Whole-exome and RNA sequencing analyses of pre-treatment tissue samples were conducted to identify determinants of response. Genomic alterations were analyzed for 10 patients, and gene expression was analyzed for 9 patients using their pre-treatment tissue samples. Targeted sequencing of cell-free DNA (cfDNA) samples (liquid biopsy) was also conducted sequentially at three points [before treatment with milademetan (baseline), at Cycle 2 Day1, and at the time of disease progression] to identify determinants of response and resistance. From whole-exome and RNA sequencing analyses of pre-treatment tissue samples, we could not find any molecular pathways associated with the anti-tumor activity of milademetan. Focusing on 8 genes (CDK4, CDKN2A, CDKN2B, EGFR, ERBB3, MDM2, PDGFRA, TP53) known to be frequently affected in intimal sarcoma and 10 genes (AKT1, ATM, BBC3, CDKN1A, CDKN1C, CHEK2, MDM4, PMAIP1, PPM1D, TWIST1) reported to be associated with MDM2 inhibitor responses; we found that anti-tumor activity correlated with amplified TWIST1 (p-value = 0.028) and negatively with CDKN2A loss (p-value = 0.071). Eight of the 10 patients had their cfDNA collected sequentially at baseline, at Cycle 2 Day1 and disease progression; however, one did not consent to the exploratory analysis study using cfDNA, and one had cfDNA collected at baseline but not at disease progression due to ongoing treatment. Of the eight patients, TP53 mutations in cfDNA were detected in one and five patients at baseline and disease progression, respectively. The cfDNA allele frequency of TP53 mutations increased with disease progression. CDKN2A loss and amplified TWIST1 could be associated with the anti-tumor activity of milademetan in patients with MDM2 amplified intimal sarcoma. Acquired TP53 mutations were detected in sequential liquid biopsies as loss-of-function mutations, and these TP53 mutations might compromise the anti-tumor activity.

  Study JGAS000619

• Integrated analysis of relapsed B-cell precursor Acute Lymphoblastic Leukemia identifies subtype-specific cytokine and metabolic signatures

  Recent efforts reclassified B-Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) into more refined subtypes. Nevertheless, outcomes of relapsed BCP-ALL remain unsatisfactory, particularly in adult patients where the molecular basis of relapse is still poorly understood. To elucidate the evolution of relapse in BCP-ALL, we established a comprehensive multi-omics dataset including DNA-sequencing, RNA-sequencing, DNA methylation array and proteome MASS-spec data from matched diagnosis and relapse samples of BCP-ALL patients (n = 50) including the subtypes DUX4, Ph-like and two aneuploid subtypes. Relapse-specific alterations were enriched for chromatin modifiers, nucleotide and steroid metabolism including the novel candidates FPGS, AGBL and ZNF483. The proteome expression analysis unraveled deregulation of metabolic pathways at relapse including the key proteins G6PD, TKT, GPI and PGD. Moreover, we identified a novel relapse-specific gene signature specific for DUX4 BCP-ALL patients highlighting chemotaxis and cytokine environment as a possible driver event at relapse. This study presents novel insights at distinct molecular levels of relapsed BCP-ALL based on a comprehensive multi-omics integrated data set including a valuable proteomics data set. The relapse specific aberrations reveal metabolic signatures on genomic and proteomic levels in BCP-ALL relapse. Furthermore, the chemokine expression signature in DUX4 relapse underscores the distinct status of DUX4-fusion BCP-ALL.

  Study EGAS00001002856

• Immune Responses in Checkpoint Myocarditis Across Heart, Blood, and Tumor

  Immune checkpoint inhibitors (ICIs) are widely used in anti-cancer therapies, but they can cause morbid and potentially fatal immune-related adverse events (irAEs), such as ICI-related myocarditis (irMyocarditis). The pathogenesis of irMyocarditis and its relationship to anti-tumor immunity remain poorly understood. We sought to define immune responses in heart, tumor, and blood in irMyocarditis patients and controls by leveraging single‐cell RNA sequencing (scRNA‐seq) coupled with T-cell Receptor Sequencing (TCR-Seq). Our analysis demonstrated increased frequencies of cytotoxic T cells, inflammatory mononuclear phagocytes (MNPs), conventional dendritic cells (cDCs), and inflammatory fibroblasts in irMyocarditis heart tissue. Additionally, we revealed decreased frequencies of plasmacytoid dendritic cells, cDCs, and B lineage cells but an increased frequency of MNPs in the blood of irMyocarditis patients. Raw data for heart and blood can be found in dbGaP. All processed data and raw TCRseq data from heart tissue and tumor samples for this study can be found in GEO (Accession number GSE228597).Subject IDs with the prefix "SIC" were collected by the Severe Immunotherapy Complications (SIC) service at MGH due to suspicion of immune-related adverse events. Subject IDs with the prefix "donor" were collected by the MGH Melanoma Biobank at scheduled timepoints; this biobank systematically collects samples from patients at MGH at clinically relevant timepoints.

  Study phs003413

• Whole-Genome Sequencing Suggests Schizophrenia Risk Mechanisms in Humans with 22q11.2 Deletion Syndrome

  Chromosome 22q11.2 microdeletions impart a high but incomplete risk for schizophrenia. Possible mechanisms include genome-wide effects of DGCR8 haploinsufficiency. In a proof-of-principle study to assess the power of this model, we used high-quality, whole-genome sequencing of nine individuals with 22q11.2 deletions and extreme phenotypes (schizophrenia, or no psychotic disorder at age >50 years). The schizophrenia group had a greater burden of rare, damaging variants impacting protein-coding neurofunctional genes, including genes involved in neuron projection (nominal P = 0.02, joint burden of three variant types). Variants in the intact 22q11.2 region were not major contributors. Restricting to genes affected by a DGCR8 mechanism tended to amplify between-group differences. Damaging variants in highly conserved long intergenic noncoding RNA genes also were enriched in the schizophrenia group (nominal P = 0.04). The findings support the 22q11.2 deletion model as a threshold-lowering first hit for schizophrenia risk. If applied to a larger and thus better-powered cohort, this appears to be a promising approach to identify genome-wide rare variants in coding and noncoding sequence that perturb gene networks relevant to idiopathic schizophrenia. Similarly designed studies exploiting genetic models may prove useful to help delineate the genetic architecture of other complex phenotypes.

  Study EGAS00001002344

• Identifying transcriptional programs underlying anti-EGFR small molecule response and resistance with TraCe-seq

  Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.

  Study EGAS00001005405

• Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

  Background Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states. Methods We exposed human primary monocytes from healthy donors (n=6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2hr) and two models of endotoxin tolerance: repeated stimulations (6+6hr) and prolonged exposure to endotoxin (24hr). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing. Results We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6884 of 34616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease. Conclusion This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

  Study EGAS00001007362

• 13 WGS and 8 TCS data of ENKL patients treated with pembrolizumab.

  Whole-genome sequencing (WGS) was performed for 13 pairs of tumor-normal samples from patients diagnosed with NKTL. Genomic DNA from tumor tissue was extracted with QIAamp DNA Mini Kit. The DNA for the matching normal was obtained from blood or buccal swabs and purified by Blood and Cell Culture DNA Mini kit or E.Z.N.A. Tissue DNA Kit (Omega Bio-tek) according to manufacturer’s instructions. The quantity and quality were assessed by Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen) and agarose gel electrophoresis. All sequencing libraries were prepared using TruSeq Nano DNA Library Prep Kit (Illumina). Paired-end sequencing was performed on Illumina HiSeq 2000 or HiSeq X Ten as 2x101 bp or 2x151 bp, respectively. 8 NKTL FFPE specimens were screened for somatic mutations using deep targeted capture sequencing (TCS). FFPE rolls or slides were extracted using QIAamp DNA FFPE Tissue kit (QIAGEN). The FFPE genomic DNA was treated with NEBNext FFPE DNA Repair Mix and assessed by Quant-it PicoGreen dsDNA Assay Kit (Invitrogen). The library was generated from 10-200 ng DNA with SureSelectXT Low Input Target Enrichment System for Illumina Paired-End Sequencing Library (Agilent Technologies) according to manufacturer’s instructions. RNA based probe was designed with SureDesign (Agilent Technologies) to target-capture 140 genes. Next, the captured libraries were pooled in equimolar concentration and sequenced on Illumina Novaseq 6000 platform with SP or S1 chip. Reads aligning to 40 selected genes were isolated post-alignment for this submission. Prefix used in filenames: T - Tumor samples N - Matched-Normal samples

  Dataset EGAD00001004140

• Fecal microbiome of T2D patients undergoing semaglutide or empagliflozin treatment

  The gut microbiome has been shown to be affected by the use of many human-targeted medications, and the interaction can be bidirectional. This has been clearly demonstrated for type 2 diabetes medications that have been in clinical use for several decades. However, the bidirectional effects of novel type 2 diabetes drugs semaglutide, empagliflozin, and the gut microbiome have yet to be clearly described. We investigated the effect of semaglutide and empagliflozin initiation on the gut microbiome of type 2 diabetes patients. In addition, we analyzed whether the pre-treatment gut microbiome can predict the treatment efficacy. Twenty subjects (10 women and 10 men) were enrolled in the study between November 2019 and January 2023 from the University of Tartu Hospital, Tartu, Estonia. All participants had been diagnosed with type 2 diabetes. The study group taking semaglutide included 6 women and 4 men, and the group prescribed empagliflozin included 4 women and 6 men. Gut microbiome fecal samples donated at four timepoints (Baseline, Month 1, Month 3; Month 12) were studied using 16S ribosomal RNA gene sequencing and analysis.

  Dataset EGAD50000000756

• Determination of cell specific regulatory enhancers in hematopoetic models

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

  Dataset EGAD00001002205

• Molecular_risk_stratification_in_patients_with_T1_colorectal_cancer_WES

  The introduction of bowel cancer screening has led to a significant increase in the proportion of patients being diagnosed with asymptomatic, early-stage colorectal cancer (CRC). Although the majority of these patients are successfully treated with surgery alone, a small proportion of patients have ‘born-to-be-bad’ aggressive lesions with early dissemination leading to distant metastases. Current standard of care histological assessment is unable to distinguish between these aggressive versus non-aggressive early lesions which is essential to provide appropriate clinical management decisions. This study aims to carry out molecular and histological profiling of approximately 300 T1 CRCs in order to develop a molecular stratifier based on the risk of relapse in early-invasive CRC. This novel T1 cohort will represent the world’s largest molecularly characterised T1 cohort of samples, with digital pathology assessment alongside whole exome sequencing, copy number variation analysis and 3’ RNA-seq. This data will be used to generate a robust panel of molecular and/or histological markers applicable to formalin-fixed paraffin embedded (FFPE) archival tissue which discriminates between T1 lesions based on risk of relapse, which will ultimately be used to inform clinical management of CRC at the earliest stages of the disease.

  Study EGAS00001005734

• The impact of BNT162b2 mRNA vaccine against SARS-CoV-2 on adaptive and innate immune responses

  The mRNA-based BNT162b2 protects against severe disease and mortality caused by SARS-CoV-2 through induction of specific antibody and T-cell responses. Much less is known about its broad effects on immune responses against other pathogens. In the present study, we investigated the specific adaptive immune responses induced by BNT162b2 vaccination against various SARS-CoV-2 variants, as well as its effects on the responsiveness of human immune cells upon stimulation with heterologous viral, bacterial, and fungal pathogens. BNT162b2 vaccination induced effective humoral and cellular immunity against SARS-CoV-2 that started to wane after six months. We also observed long-term transcriptional changes in immune cells after vaccination, as assessed by RNA sequencing. Additionally, vaccination with BNT162b2 modulated innate immune responses as measured by the production of inflammatory cytokines when stimulated with various microbial stimuli other than SARS-CoV-2, including higher IL-1/IL-6 release and decreased production of IFN-α. Altogether, these data expand our knowledge regarding the overall immunological effects of this new class of vaccines and underline the need of additional studies to elucidate their effects on both innate and adaptive immune responses.

  Study EGAS00001006818

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

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

  Study EGAS00001006747

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

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

  Study EGAS00001003736

• Biased allelic expression in human primary fibroblast single cells.

  The study of gene expression in mammalian single cells using genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We have used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblast cells over 133,633 unique heterozygous single nucleotide variants (hetSNVs). We have observed that at the snapshot of analyses each cell contains mostly transcripts from one allele from the majority of genes; indeed 76.4% of the hetSNVs display stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibit haplotype consistent allelic ratio; in contrast distant sites located in two different genes are independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in majority of the single cells at a given time point are rare and enriched in highly expressed genes. The relative abundance of each allele in a cell is controlled by some regulatory mechanisms since we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

  Study EGAS00001001009

• SINGLE-CELL RNA SEQUENCING Single-cell RNA sequencing was performed on 13 ‘mild-moderate’ and 10 ‘critical’ COVID19 PBMC samples

  Summary Epidemiological and clinical reports have indicated that SARS-CoV-2 virulence hinges upon the triggering of an aberrant host immune response, more so than on direct virus-induced cellular damage. To elucidate the immunopathology underlying COVID-19 severity, we performed cytokine and multiplex immune profiling in mild-moderate and critically ill COVID-19 patients. Hypercytokinemia in COVID-19 differed from the IFN-γ-driven cytokine storm in macrophage activation syndrome, and was more pronounced in critical versus mild-moderate COVID-19. Systems modelling of cytokine levels paired with deep-immune profiling showed that classical monocytes drive this hyper-inflammatory phenotype and that a reduction in T-lymphocytes correlates with disease severity, with CD8+ cells being disproportionately affected. Expression of antigen presenting machinery was also reduced in critical disease. Furthermore, we found that neutrophils contributed to disease severity and local tissue damage by amplification of hypercytokinemia and the formation of neutrophil extracellular traps. Together these findings suggest a myeloid-driven immunopathology, in which hyperactivated neutrophils and an ineffective adaptive immune system act as mediators of COVID-19 disease severity.

  Study EGAS00001005039

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

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

  Study EGAS00001005675

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

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

  Study EGAS00001005874

• T cell landscape definition by multi-omics identifies galectin-9 as novel immunotherapy target in chronic lymphocytic leukemia (CLL)

  Failure of immunotherapy after applying checkpoint inhibitors or CAR-T cells is linked to T cell exhaustion. Here, we explored the T cell landscape in chronic lymphocytic leukemia (CLL) by single-cell omics analyses of blood, bone marrow and lymph node samples of patients and spleen samples of a CLL mouse model. By single-cell RNA- sequencing, mass cytometry (CyTOF), and multiplex image analysis of tissue microarrays, we defined the spectrum of phenotypes and transcriptional programs of T cells and and their differentiation state trajectories. We identified disease-specific accumulation of distinct regulatory T cell subsets and T cells harboring an exhausted phenotype exclusively in the CLL lymph node tissue. Integration of TCR data revealed a clonal expansion of CD8+ precursor exhausted T cells, suggesting their reactivity for CLL cells. Interactome analyses identified the TIM3 ligand Galectin-9 as novel immunoregulatory molecule in CLL. Blocking of Galectin-9 in CLL-bearing mice slowed down disease development and reduced the number of TIM3 expressing T cells. Galectin-9 expression correlated with shorter survival of CLL patients. Thus, Galectin-9 contributes to immune escape in CLL and represents a novel target for immunotherapy.

  Study EGAS00001006864

• Oncoprint GSCCs

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

  Study EGAS00001007481

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

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

  Study EGAS00001007482

• Genomic Analysis of Prostate Tumor Heterogeneity in Metastasis

  To investigate the molecular changes, mutations, and signaling pathways that characterize the initial steps of prostate cancer metastasis to the lymph nodes, we analyzed 19 discrete lymph node (LN) metastases, 97 primary tumor foci, 39 benign lymph nodes, and 10 normal adjacent prostate tissue samples from 43 patients enrolled in a clinical trial (NCT01808222) of extended lymphadenectomy and metastasis excision guided by [18F]-fluciclovine PET imaging for high risk, newly diagnosed prostate cancer conducted at Emory University. These patients fell within high or very high-risk groups (T3a, Gleason score 8-10, or PSA greater than 20 ng/ml). Of these patients, 16 had LN metastases at least 4 mm in diameter, and 20 patients had no metastases. For metastatic patients, we performed RNA and Whole Exome Sequencing (WES) on three primary tumor foci, one benign LN, and all LN metastases greater than 4mm in diameter. For non-metastatic patients, we sequenced two primary tumor foci and one benign LN. In total, RNAseq was performed on 165 tissue samples from 43 patients, and WES was performed on 137 samples from 35 patients.

  Study phs003404

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

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

  Study EGAS00001004614

• Harnessing Epigenetic Regulators to improve HSC-based lentiviral gene therapy

  Curative benefits of autologous and allogeneic transplantation of hematopoietic stem cells (HSCs) have been proven for various diseases. However, the low number of true HSCs that can be collected from patients and subsequently in vitro maintenance and expansion of true HSCs for genetic correction remain challenging. Addressing this issue, we here focused on optimizing culture conditions to improve the ex vivo expansion of true HSCs for gene therapy purposes. In particular, we explore the use of epigenetic regulators to enhance the effectiveness of HSC-based lentiviral (LV) gene therapy. The HDAC inhibitor Quisinostat and the bromodomain inhibitor CPI203 each promote ex vivo expansion of functional HSCs, as validated by xenotransplantation assays and single cell RNA-sequencing analysis. We confirmed the stealth effect of LV transduction on the loss of HSC numbers in commonly used culture protocols, while addition of Quisinostat or CPI203 improved expansion of HSCs in transduction protocols. Of note, we demonstrated that addition of Quisinostat improved LV transduction efficiency of HSCs and early progenitors. Our suggested culture conditions highlight the potential therapeutic effect of epigenetic regulators in hematopoietic stem cell biology and their clinical applications to advance HSC-based gene correction.

  Study EGAS50000000175

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - RNA-seq

  B-cell prolymphocytic leukemia (B-PLL) is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p<0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003274

• Single-cell RNA-sequencing of CD19 Chimeric Antigen Receptor (CAR) T-cells (CD19 CAR T-cells) from a Phase I clinical study in paediatric ALL: CARPALL

  Chimeric antigen receptor (CAR)-modified T-cells have become established as an effective treatment of haematological cancers. In the context of relapsed and refractory childhood pre-B cell acute lymphoblastic leukaemia (B ALL), CD19 targeting CAR T-cells often induce durable remissions. Previously, we generated a novel low-affinity CAR incorporating a CD19-specific single-chain variable fragment (scFV) called CAT, displaying a faster off-rate of interaction than the FMC63 CD19 binder used in prior clinical studies. Here, we systematically analysed CD19 CAR T-cells of ten children with relapsed or refractory B ALL enrolled in the CARPALL trial (NCT02443831). To characterize persisting CD19 CAR T-cells, we performed high throughput single-cell gene expression and T-cell receptor (TCR) sequencing of infusion products and serial blood and bone marrow samples up to five years post-infusion. We isolated CAR T-cells from peripheral blood or bone marrow by flow cytometry for CD3 and CAR expression, prior to single cell sequencing (Chromium 10X) platform.

  Dataset EGAD00001010018

• Bulk-RNA Sequencing of high-grade pancreatic and non-pancreatic Neuroendocrine Neoplasms

  Therapeutic decisions in oncology depend on a precise pathological classification of individual neoplasms. Recent years have seen an intensification of research activities aimed at the extraction of clinically relevant information from patient-derived 'omics' data based on Machine-Learning models. However, a comprehensive training of Machine-Learning models requires sufficiently large numbers of training samples, which are usually not available for rare cancer types. The problem is worsened when individual tissues segregate into different cancer subtypes, as their discrimination would require even more training samples. Methods: Here, we report on a new data-augmentation technique to support the training of Machine-Learning models on ‘omics’ data from pancreatic neuroendocrine neoplasms (panNEN). PanNENs display all properties described above: Only about 2-3% of all pancreatic neoplasms are neuroendocrine and they fall into different subtypes with distinctly different prognosis, which makes the precise classification of such samples both difficult and important for therapy decisions. The approach reconstructs a given transcriptome based on healthy pancreatic cell type signatures and creates a Machine-Learning model that integrates the observed reconstruction error and predicted cell type proportions as features. Results: A benchmark of the deconvolution model predictions with the ground-truth found that the model could efficiently predict the sample grading, disease-related patient survival time, and differentiate between different subtypes in four panNEN and one mixed panNEN and non-pancreatic NENs dataset. We compared the predictive performance of the deconvolution-trained model to that of a model trained directly on the transcriptomic data, under inclusion of the Ki-67 classification gold standard biomarker, and found the performances to be comparable. Conclusion: Our approach serves as a data-augmentation technique to facilitate the training of Machine Learning models for rare cancer types via utilization of data from healthy origin. Additionally, the deconvolution results were clinically interpretable and further research can render the deconvolution approach an effective complementary asset for the clinical classification of neoplasms.

  Study EGAS00001004861

• Genetics of gene expression in human macrophage response to Salmonella

  Understanding how cells sense and respond to their environment, and how these responses are modulated by genetic variation, are fundamental biological problems, particularly for understanding how pathogenic organisms invade and manipulate the cells of the human immune system. Macrophages recognize and respond to many important human pathogens including HIV-1, Mycobacteria tuberculosis and Salmonella. This study will focus on the cellular response of human macrophages to Salmonella infection and how this response is modulated by the genetic bacground of the individual as well as additional pro-inflammatory stimulus (interferon-gamma priming). We will acquire 100 human induced pluripotent stem cell lines from the HipSci project, differentiate the cells in vitro into macrophages and expose them to four environmental conditions: (i) no stimulation, (ii) interferon-gamma (18h), (iii) Salmonella typhimurium SL1344 (5h), (iv) interferon-gamma (18h) + Salmonella (5h).Subsequently, we will isolate RNA from the samples for sequencing.

  Dataset EGAD00001003204

• Distinct embryonic phylogenies and driver events of infant Wilms tumor - RNA

  Cancers of adults typically arise through progressive rounds of clonal diversification and intratumoral selective sweeps which generate a long mutational trunk with shorter subclonal branches. Here, we investigated whether tumors of young children exhibit the same phylogenetic configuration. We studied three infants, including two newborns, with the childhood kidney cancer, Wilms tumour, through whole genome sequencing of bulk tissues, of single cell derived organoids, and of microdissections. All three cancers exhibited unusual driver events, with tumours of newborns harbouring FOXR2 rearrangements, delineating a distinct variant of Wilms tumour. Phylogenetic analyses suggest that tumors were seeded in an early, possibly confined window of development. Unusually, following seeding there was extensive polyclonal diversification with little evidence of clonal sweeps, leading to a distinct phylogenetic configuration more reminiscent of normal tissues rather than of adult cancers. These findings indicate that some childhood cancers may diversify via unorthodox phylogenetic pathways.

  Dataset EGAD00001009813

• Single-cell and bulk RNA-sequencing of nivolumab-treated glioblastoma

  Glioblastoma is a highly aggressive type of brain tumour for which there is no curative treatment available. Immunotherapies have shown limited responses in unselected glioblastoma patients, and it is currently not well understood what separates the few responders from the non-responders at the cellular and molecular level. Here we investigate the phenotypic and transcriptional evolutionary dynamics at single-cell resolution during nivolumab immune checkpoint treatment of eight glioblastoma patients enrolled in a translational phase I/II clinical trial.

  Study EGAS00001007110

• Exceptional Responders Initiative

  The Exceptional Responders Initiative is a pilot study to investigate the underlying molecular factors driving exceptional treatment responses of cancer patients to drug therapies. Study researchers will examine molecular profiles of tumors from patients either enrolled in a clinical trial for an investigational drug(s) and who achieved an exceptional response relative to other trial participants, or who achieved an exceptional response to a non-investigational chemotherapy. An exceptional response is defined as achievement of either a complete response or a partial response for at least 6 months duration in a trial or treatment where the overall response rate is < 10%. The hope is to discover underlying molecular features that can be further investigated and may eventually predict benefit from a given drug or class of drugs for a particular patient. This pilot project will successfully characterize approximately 100 cases of tumor tissue and, when available, case-matched germline DNA. All samples will undergo whole exome sequencing, and cases with sufficient nucleic acids will undergo additional analyses (e.g. whole genome sequencing, mRNA-sequencing, mi RNA sequencing, promoter methylation analysis, SNP etc). Each case will be annotated with demographic and clinical information, along with follow-up information minimally sufficient to correlate molecular profiles with response. Both retrospective and prospective collections will be considered. The project will also accept sequencing data and clinical data from patients who have had sequencing performed outside of this project. All data will be de-identified and placed in a controlled access database so other investigators may use them for additional insights. Clinically annotated tissue specimens meeting the criteria will be provided by groups participating in the Exceptional Cases Initiative to a Biospecimen Core Resource (BCR), which will perform quality control on the tissues, and will use a standard operating procedure to isolate nucleic acids. The nucleic acids will be shipped to a sequencing center to perform whole exome sequencing and analysis. These findings will be made available to the broader cancer research community in a controlled access database.

  Study phs001145

• Renal_habitat_WXS

  This project aims at correlating genetic diversity (as well as transcriptomic, metabolic and immune diversity) observed in multiple targeted tumour biopsies and circulating tumour DNA with heterogeneity on multiparametric MRI. To this end, 20 patients with RCC (15 patients with localised/locally advanced disease and 5 with mRCC) will undergo multiparametric MRI including metabolic assessment by 13C-MRSI as well as measurements of tumour perfusion, oxygenation and water diffusion. Based on unique combinations of these imaging parameters, tumour habitats will be defined through cluster mapping. Intra-operative biopsy and immediate freezing of tumour cores from different imaging habitats will be performed. Multiregional ex-vivo tumour biopsies will be spatially registered to the MR data by using 3D printed tumour moulds generated from high-resolution volumetric sequences. These moulds allow slicing tumour specimens in identical planes to the MR acquisition and increase the accuracy of tissue sampling. The biopsies will then undergo molecular profiling such as metabolomic and RNA sequencing (including mRNA immune profiling), as well as DNA sequencing to assess driver and clonal mutations in different tumour habitats and guide design of appropriate cfDNA probes. The targeted cfDNA capture panel will be employed to determine the ability of cfDNA to capture the full genetic diversity observed in the primary tumour and to relate cfDNA to various imaging habitats in both renal and pancreatic cancer. Data integration will be undertaken between imaging habitats and molecular profiles for each region.

  Study EGAS00001003703

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

  This dataset includes 3432 paired single cell sequencing fastq files derived from synovial B cells of 5 early Rheumatoid Arthritis patients. Libraries were prepared using the SmartSeq2 protocol. All wells contained ERCC spike-ins. Libraries were sequenced on an Illumina NovaSeq instrument with 2x100bp paired-end reads yielding a median of 2.5M reads/well. Sample aliases ending on 368 and 384 represent empty control wells. Sample aliases starting with P11417_4 and P13157_1 belong to ACPA- patient A7, sample aliases starting with P11417_5 and P11417_6 belong to ACPA+ patient A1, sample aliases starting with P13157_2 belong to ACPA- patient A3, sample aliases starting with P13157_3 and P13157_6 belong to ACPA+ patient A2, sample aliases starting with P13157_4 and P13157_5 belong to ACPA- patient A4.

  Dataset EGAD00001009076

• Identification of driver oncogenes in scirrhous-type gastric cancer cell lines

  Scirrhous type gastric cancer (SGC) is one of the most intractable cancer subtypes in human, and its therapeutic targets have been rarely identified to date. Exploration of somatic mutations in the SGC genome with the next generation sequencers has been hampered by low tumor contents by the markedly increased fibrous tissues. Thus, SGC cell lines may be useful resources to search for novel oncogenes in. Here we have conducted whole exome sequencing and RNA sequencing on two SGC cell lines, OCUM-8 and OCUM-9. Interestingly, most of the mutations thus identified have not been reported. In OCUM-8 cells, a novel CD44-IGF1R fusion gene is discovered, the protein product of which ligates the amino-terminus of CD44 to the transmembrane and the tyrosine-kinase domains of IGF1R. Further, both of CD44 and IGF1R are markedly amplified in the OCUM-8 genome and abundantly expressed. CD44-IGF1R has a transforming ability, and the suppression of its kinase activity leads to rapid cell death of OCUM-8. On the other hand, OCUM-9 seems to possess multiple oncogenic events in its genome. The MET gene is highly amplified, and the treatment of OCUM-9 with crizotinib, a MET inhibitor, results in cell death. Further, a novel BORCS5-ETV6 fusion gene is identified in the OCUM-9 genome. BORCS5-ETV6 possesses oncogenic activity, and suppression of its message partially inhibits cell growth. Prevalence of these novel fusion genes among SGC awaits further investigation, but we validate the significance of cell lines as appropriate reagents for detailed genomic analyses of SGC.

  Study JGAS000179

• Coding and non-coding drivers of mantle cell lymphoma identified through exome and genome sequencing

  Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma (NHL) that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established and the features that contribute to differences in clinical course remain limited. To extend our understanding of the biological pathways involved in this malignancy, we performed a large-scale genomic analysis of MCL using data from 51 exomes and 34 genomes alongside previously published exome cohorts. To confirm our findings, we re-sequenced the genes identified in the exome cohort in 191 MCL tumors, each having clinical follow-up data. We confirmed the prognostic association of TP53 and NOTCH1 mutations. Our sequencing revealed novel recurrent non-coding mutations surrounding a single exon of the HNRNPH1 gene. In RNA-seq data from 103 of these cases, MCL tumors with these mutations had a distinct imbalance of HNRNPH1 isoforms. This altered splicing of HNRNPH1 was associated with inferior outcomes in MCL and showed a significant increase in protein expression by immunohistochemistry. We describe a functional role for these recurrent non-coding mutations in disrupting an auto-regulatory feedback mechanism, thereby deregulating HNRNPH1 protein expression. Taken together, these data strongly implicate a role for aberrant regulation of mRNA processing in MCL pathobiology.

  Study EGAS00001004289

• Genomic Analysis of Extra-Nodal Natural Killer/T-Cell Lymphoma (ENKTCL)

  In this study, we performed extensive genetic profiling on 209 extra-nodal natural killer/T-cell lymphoma (ENKTCL) nasal-type cases, including 174 Chinese from four hospitals in China, and 31 Hispanic and 4 Caucasian patients from City of Hope and Eberhard Karls University of Tübingen. All of the cases were reviewed by experienced hematopathologists to confirm diagnosis. Written informed consents from patients were reviewed and approved by the institutional review board of each institution. We deposited targeted exon sequencing data of a selected 333-gene panel on the ENKTCL cases in dbGaP. DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissues with an AllPrep DNA/RNA FFPE Kit. A custom platform (Agilent SureSelect) was used to capture target exons from genomic DNA of ENKTCL tumors. DNA was fragmented to 200bp using a Covaris S220 Focused-ultrasonicator and was end-repaired and ligated to Illumina adaptors with a KAPA Hyper Prep Kit. Ligated products were purified and amplified with a 7-cycle PCR. Exome capture from the PCR products was performed using a SureSelect Target Enrichment Kit. The captured exomes were further amplified with an 8-cycle PCR and purified for cluster generation using the cBot cluster generation system with a HiSeq PE Cluster Kit V4. 100bp paired-end reads were generated from an Illumina HiSeq 2500 system with a HiSeq SBS Kit V4.

  Study phs002925

• DNA and RNA sequencing of single human haploid germ cells

  Our goal was to look at allele-specific gene expression in single spermatids from two human donors. We first genotyped and phased each donor using a combination of linked-read sequencing (one sample per donor, 10X Chromium), and single haploid whole genome sequencing (53 spermatids for the first donor, and 19 for the second). We then sequenced mRNA from single spermatids (302 from the first donor, and 164 from the second), using the Smart-Seq2 protocol. We used these gene expression data, in combination with the haplotype data, to look for gene products that exhibit allele-specific bias reflecting the haploid genotype of the cell. We find a large class of human genes that show such allele-specific bias in spermatids, indicating that these gene products are not fully shared across the cytoplasmic bridges that connect spermatids. We term these genes “genoinformative markers” (GIMs).

  Study phs002279

• Mutational patterns and regulatory networks in epigenetic subgroups of meningioma (H033)

  DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign-1, 2, 3, intermediate-A, B and malignant. Here, we set out to identify subgroup-specific mutational patterns and pathway regulation. Whole-genome-sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic and 22 radiation-induced (Mrad) meningiomas. RNA sequencing was added for 18 cases and chromatin-immunoprecipitation for the enhancer mark H3K27ac followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural alterations were found to contribute to the spectrum of mechanisms inactivating NF2 in sporadic meningioma similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3 was detected in both sporadic meningioma and Mrad, but distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. In general, the malignant MC presented a significantly higher exposure to AC3 and higher genomic instability beyond the mutational load than the other MCs. Pathway analysis of the ChIP-seq clusters revealed that the FOXM1 is most differentially activated in high-grade cases, along with super enhancer recruitment near HOX genes and respective upregulation of expression. This data further elucidate the biological mechanisms differentiating meningiomas of different WHO grade and epigenetic subgroups and suggest leveraging the genomic instability as novel therapeutic targets.

  Study EGAS00001003481

• Characterisation of the genetic landscape of cutaneous leiomyoma and leiomyosarcoma - Leiomyoma WES

  Cutaneous leiomyoma (cLM) and leiomyosarcoma (cLMS) are rare benign and malignant soft tissue neoplasms showing smooth muscle differentiation, respectively, that arise from mesenchymal cells in the dermis and subcutis. Through whole exome sequencing of cLM and cLMS cases, we observed distinct differences between the somatic mutational profile of these tumour types. FH was identified as a driver gene in cLM with genetic alterations of FH occurring via somatic point mutation, somatic copy number loss, biallelic inactivation and germline point mutations. TP53 and RB1 were identified as driver genes in the cLMS cohort, with genetic alterations of TP53 occurring via somatic and germline point mutations, copy number loss and biallelic inactivation. Using RNA-sequencing, we identified recurrent gene fusions, including CRTC1/3-MAML2 in cLMS and a novel MYLK-MAP3K2 fusion. Analysis of the cell types present in the tumour microenvironment revealed a significantly increased presence of macrophages and decreased presence of myeloid dendritic cells in the cLMS cohort relative to the cLM cohort. Additionally, we identified common driver genes between cLMS and LMS from other sites. Thus, we provide the first in-depth profile of the genetic landscape of cLM and cLMS.

  Dataset EGAD00001014787

• Characterisation of the genetic landscape of cutaneous leiomyoma and leiomyosarcoma - Leiomyosarcoma WES

  Cutaneous leiomyoma (cLM) and leiomyosarcoma (cLMS) are rare benign and malignant soft tissue neoplasms showing smooth muscle differentiation, respectively, that arise from mesenchymal cells in the dermis and subcutis. Through whole exome sequencing of cLM and cLMS cases, we observed distinct differences between the somatic mutational profile of these tumour types. FH was identified as a driver gene in cLM with genetic alterations of FH occurring via somatic point mutation, somatic copy number loss, biallelic inactivation and germline point mutations. TP53 and RB1 were identified as driver genes in the cLMS cohort, with genetic alterations of TP53 occurring via somatic and germline point mutations, copy number loss and biallelic inactivation. Using RNA-sequencing, we identified recurrent gene fusions, including CRTC1/3-MAML2 in cLMS and a novel MYLK-MAP3K2 fusion. Analysis of the cell types present in the tumour microenvironment revealed a significantly increased presence of macrophages and decreased presence of myeloid dendritic cells in the cLMS cohort relative to the cLM cohort. Additionally, we identified common driver genes between cLMS and LMS from other sites. Thus, we provide the first in-depth profile of the genetic landscape of cLM and cLMS.

  Dataset EGAD00001014789

• Discovering the Genetic Basis of Human Neuroblastoma: A Gabriella Miller Kids First Pediatric Research Program (Kids First) Project

  The Gabriella Miller Kids First Pediatric Research Program (Kids First) is a trans-NIH effort initiated in response to the 2014 Gabriella Miller Kids First Research Act and supported by the NIH Common Fund. This program focuses on gene discovery in pediatric cancers and structural birth defects and the development of the Gabriella Miller Kids First Pediatric Data Resource (Kids First Data Resource). All of the WGS and phenotypic data from this study are accessible through dbGaP and kidsfirstdrc.org, where other Kids First datasets can also be accessed. Children with disseminated neuroblastoma have a very high risk of treatment failure and death despite receiving intensified chemotherapy, radiation therapy and immunotherapy. The long-term goal of our research program is to ultimately improve neuroblastoma cure rates by first comprehensively defining the genetic basis of the disease. The central hypothesis to be tested here is that neuroblastoma arises largely due to the epistatic interaction of common and rare heritable DNA variation. Here we will perform a comprehensive whole genome sequencing of 563 quartets of neuroblastoma patient germline and diagnostic tumor DNAs and germline DNAs from both parents. The case series was recently collected through a Children's Oncology Group epidemiology clinical trial and is robustly annotated with complete demographic (age, sex, race, ethnicity), clinical (e.g. age at diagnosis, stage, risk group), epidemiologic (parental dietary and exposure questionnaire) and biological (e.g. tumor MYCN status and multiple other tumor genomic measures) co-variates. Subjects were consented for genetic research and DNA is immediately available for shipment for sequencing. We propose Illumina-based whole genome sequencing in the 593 "trio" germline samples (Aim 1; due to missing parent: 487 full neuroblastoma triads, 106 child-single parent dyads = 1673 whole genome sequences) and matched diagnostic tumor DNA (Aim 2; N=366) at 30x sequencing depth (N=2039 whole genome sequences). Also in Aim 2 we will perform whole exome (100x) and RNA sequencing on the 366 tumor DNA and 228 tumor RNA samples from this cohort. Finally, we propose a pilot study of structural variation using long-range sequencing in 10 non-overlapping tumor samples chosen based on potentially relevant chromosomal alterations discovered with conventional NGS. Thus, a total of 2277 individual samples and 2655 sequences will be generated. We will use our established analytic pipeline that is currently being used to study the germline genomes of all cases sequenced through the NCI supported Therapeutically Applicable Research to Generate Effective Treatments program. We plan a three-stage analytic approach, first focusing on classic de novo and inherited Mendelian damaging alterations. We will next integrate our extensive epigenomic data from human neuroblastoma cell lines and genome-wide association study data (N=5,703 neuroblastoma cases to date) to guide a comprehensive assessment of noncoding variants that influence tumor initiation with a recently established analytic pipeline. Finally, we will utilize the tumor DNA analyses to inform relevance via somatic gain or loss of function effects at the sequence and/or copy number levels. All data generated in this project will be immediately placed into the Genomic Data Commons (GDC) and we will compute within this environment by importing our analytic pipelines into the GDC. These data will be fully integrated into the Kids First Data Resource and freely shared with all academically qualified petitioners. This comprehensive data set derived from a large and richly phenotyped series of neuroblastoma DNA quartets will be integrated with existing germline and/or tumor genomic data from over 6,000 neuroblastoma subjects (but none with matched patient-parent germline sequencing data) to provide an unparalleled opportunity to comprehensively discover the genetic basis of neuroblastoma.

  Study phs001436

• Regulatory and Effector T cells in Oral Immunotherapy for Food Allergy

  This study consisted of a clinical trial for peanut oral immunotherapy (OIT) and a follow-up project (see below). Patients were included from 2013-2015, and samples from the last timepoints were collected in 2017. The study focused on peanut-specific CD4+ T cell responses before, during, and after OIT with peanut protein. We analyzed T cell gene expression, cytokine production, and T cell receptor β (TCRβ) and TCRα clonotypes, both in bulk T cell populations and single cells, and sought to correlate these data with differences in clinical sensitivity between peanut-allergic patients before OIT, and with differences in clinical outcome after OIT. The goal was to better understand the role of peanut-specific effector and regulatory T cells in peanut allergy and tolerance, which may inform more effective therapies (e.g., OIT combined with biologicals, probiotics, or adjuvants) in the future. Tolerance Following Peanut Oral Immunotherapy (PNOIT2, NCT01750879): The unifying objective of this project is to determine whether peanut oral immunotherapy (PNOIT)-induced clinical tolerance in the context of food allergy is significantly associated with the expansion of a specific regulatory T cell subset (CD45RA-, CD25++, FoxP3+ +) that is thought to be inducible in the gut-associated lymphoid compartment and associated with immunological tolerance. The hypothesis of the study is that the induction of Treg cells will be associated with clinical tolerance. The investigators will measure the change from baseline of induced Treg cells as a frequency of total CD4+ T cells during active treatment and compare that between participants who achieve significant clinical tolerance and those who do not. High Threshold Peanut Challenge Study (PAID-UP, NCT02698033): This protocol was designed to better characterize a sub-population of peanut-sensitized individuals who may be non-allergic, despite significant sensitization, or who may be allergic, but at high threshold doses. By specifically targeting participants who met the initial screening criteria of an earlier study, NCT01750879, but failed to react during the pre-treatment 443 mg challenge to peanut, the investigators anticipate that they will identify individuals who have become spontaneously tolerant, despite persistent sensitization. The investigators might also find that clinical sensitivity persists but only with higher thresholds, or that sensitivity has increased (or is variable) since the previous allergen exposure. By repeating food challenges through to a full serving dose (7.4 gram), the investigators will distinguish participants who react only at higher doses from those who were not truly peanut-allergic, address whether their sensitivity has changed, and have the opportunity to further investigate their immune response to peanut allergen. Methods: Peripheral blood mononuclear cells were isolated from patient blood samples, and stimulated with peanut protein. Peanut-activated CD4+ T cells were sorted based on the activation marker CD154. For a more recent study, the activation marker CD137 was used as well. Total RNA from peanut-activated CD154+ and resting CD154- CD69- T cells was used for cDNA synthesis and amplification (SMARTer ultra low input RNA kit for sequencing - v3; Clontech Laboratories). Libraries were prepared and sequenced on the Illumina HiSeq platform, at a read depth of approximately 30 million reads per sample. Genomic DNA from CD154+ and CD154- CD69- T cells was used to amplify and sequence the TCRβ CDR3 regions (immunoSEQ assay; Adaptive Biotechnologies). The immunoSEQ approach generates an 87 base-pair fragment capable of identifying the VDJ region spanning each unique CDR3. Amplicons were sequenced using the Illumina NextSeq platform. Using a baseline developed from a suite of synthetic templates, primer concentrations and computational corrections were used to correct for the primer bias common to multiplex PCR reactions. Raw sequence data were filtered on the basis of TCRβ V, D, and J gene definitions provided by the IMGT database (www.imgt.org) and binned using a modified nearest-neighbor algorithm to merge closely related sequences and remove both PCR and sequencing errors. For single-cell analysis, we aimed to recover TCR variable sequences from high-throughput single-cell libraries generated by popular platforms such as SeqWell and DropSeq. We paired TCR information with T cell single-cell RNA sequencing results from patients undergoing peanut OIT to show clonal expansion-related T cell phenotypes.

  Study phs001897

• Clonal_haematopoiesis_in_patients_with_AAA

  Recent advances in genomics have demonstrated that clonal haemopoiesis driven by leukaemia associated somatic mutations is a relatively common phenomenon that increases in frequency with advancing age. Whilst individuals with clonal haemopoiesis have an increased risk of developing haematological malignancies, they also have an increased mortality from other causes. Additionally, certain mutations are almost exclusively seen in individuals aged 70 years or older, whilst others are seen in individuals with non-haematological cancers including breast and ovarian. Recently, clonal haemopoiesis was found to be associated with a significantly increased risk of atherosclerotic cardiovascular disease. This association is thought to be causative with clonally-derived macrophages showing elevated expression of several chemokine and cytokine genes that contribute to atherosclerosis. Another vascular pathology, abdominal aortic aneurysm (AAA), increases with age and shares risk factors with atherosclerosis (including smoking, male sex, high cholesterol). However, the impact of these risk factors and the overlap between AAA and atherosclerosis is poorly understood. To investigate a possible link between clonal haemopoiesis and AAA, we will study DNA samples from 300 patients with AAA and up to 200 controls for evidence of clonal haemopoiesis. This will be done using target DNA enrichment with biotinylated RNA baits followed by high throughput sequencing.

  Study EGAS00001002873

• Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

  Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that forms intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to show treatment resistance and conserve the tumour bulk. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogenic tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. SNP arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors, consistent with the heterogenic molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

  Study EGAS00001003438

• The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma

  Patients with metastatic uveal melanoma (UM) have an abysmal prognosis. Preclinical studies have suggested that epigenetic therapy could enhance immunogenicity of cancer cells. Here we test if epigenetic therapy would enhance PD-1 immunotherapy in patients with metastatic UM. We report the results of the PEMDAC phase 2 clinical trial (n=29; NCT02697630) where the HDAC inhibitor entinostat was combined with the PD-1 inhibitor pembrolizumab in patients with metastatic UM. The primary endpoint was objective response rate (ORR), and was met with an ORR of 14%. The clinical benefit rate at 18 weeks was 28%, median progression free survival was 2.1 months and the median overall survival was 13.4 months. Toxicities were manageable, and there were no treatment-related deaths. Extensive genomics studies were performed using DNA/RNA and single cell sequencing and flow cytometry. Objective responses and prolonged survival were seen in three patients with BAP1 wildtype tumors, and in one patient with an iris melanoma that exhibited a UV signature. Longer survival also correlated with low baseline ctDNA levels or LDH. In conclusion, HDAC inhibition and anti-PD1 immunotherapy results in durable responses in a subset of patients with metastatic UM. Further exploration of combined immunotherapy and epigenetic therapy in metastatic UM is warranted.

  Study EGAS00001005478

• Longitudinal evaluation of serum microRNAs as biomarkers for neuroblastoma burden and therapeutic p53 reactivation

  Accurate assessment of treatment response and residual disease is indispensable for the evaluation of cancer treatment efficacy. However, performing tissue biopsies for longitudinal follow-up poses a major challenge in the management of solid tumours like neuroblastoma. In the present study, we evaluated whether circulating miRNAs are suitable to monitor neuroblastoma tumour burden and whether treatment-induced changes of miRNA abundance in the tumour are detectable in serum. We performed small RNA sequencing on longitudinally collected serum samples from mice carrying orthotopic neuroblastoma xenografts that were exposed to treatment with idasanutlin or temsirolimus. We identified 57 serum miRNAs to be differentially expressed upon xenograft tumour manifestation, out of which 21 were also found specifically expressed in the serum of human high-risk neuroblastoma patients. The murine serum levels of these 57 miRNAs correlated with tumour tissue expression and tumour volume, suggesting potential utility for monitoring tumour burden. In addition, we describe serum miRNAs that dynamically respond to p53 activation following treatment of engrafted mice with idasanutlin. We identified idasanutlin-induced serum miRNA expression changes upon one day and 11 days of treatment. By limiting to miRNAs with a tumour-related induction, we put forward hsa-miR-34a-5p as a potential pharmacodynamic biomarker of p53 activation in serum.

  Study EGAS00001006678

• Transcriptomics identifies blunted immunomodulatory effects of vitamin D in people with multiple sclerosis

  Vitamin D deficiency is a risk factor for developing multiple sclerosis (MS). However, the immune effects of vitamin D in people with MS are not well understood. We analyzed transcriptomic datasets generated by RNA sequencing of immune cell subsets (CD4+, CD8+ T cells, B cells, monocytes) from 33 healthy controls and 33 untreated MS cases. We utilized a traditional bioinformatic pipeline and weighted gene co-expression network analysis (WGCNA) to determine genes and pathways correlated with endogenous vitamin D. In controls, CD4+ and CD8+ T cells had 1079 and 1188 genes, respectively, whose expressions were correlated with plasma 25-hydroxyvitamin D level (P<0.05). Functional enrichment analysis identified association with TNF-alpha and MAPK signaling. In CD4+ T cells of controls, vitamin D level was associated with expression levels of several genes proximal to multiple sclerosis risk loci (P=0.01). Genes differentially associated with endogenous vitamin D by case-control status were enriched in TNF-alpha signaling via NF-κB. WGCNA suggested a blunted response to vitamin D in cases relative to controls. Collectively, our findings provide further evidence for the immune effects of vitamin D, and demonstrate a differential immune response to vitamin D in cases relative to controls, highlighting a possible mechanism contributing to MS pathophysiology.

  Study EGAS00001007254

• Identification and Targeting of Inflammatory Macrophage-Fibroblast Crosstalk in Rheumatoid Arthritis

  Macrophages tailor their function according to the signals found in tissue microenvironments, assuming a wide spectrum of phenotypes. A detailed understanding of macrophage phenotypes in human tissues is limited. Using single-cell RNA sequencing, we defined distinct macrophage subsets in the joints of patients with the autoimmune disease rheumatoid arthritis (RA), which affects ~1% of the population. The subset we refer to as HBEGF+ inflammatory macrophages is enriched in RA tissues and is shaped by resident fibroblasts and the cytokine tumor necrosis factor (TNF). These macrophages promoted fibroblast invasiveness in an epidermal growth factor receptor-dependent manner, indicating that intercellular cross-talk in this inflamed setting reshapes both cell types and contributes to fibroblast-mediated joint destruction. In an ex vivo synovial tissue assay, most medications used to treat RA patients targeted HBEGF+ inflammatory macrophages; however, in some cases, redirecting them into a state is not expected to resolve inflammation. These data highlight how advances in our understanding of chronically inflamed human tissues and the effects of medications therein can be achieved by studies on local macrophage phenotypes and intercellular interactions. Reprinted from Kuo, Ding et al., Science Translational Medicine 2019, PMID: 31068444, with permission from American Association for the Advancement of Science.

  Study phs001340

• Generation of a Cellular Atlas of Human Adipose Tissue

  White adipose tissue (WAT), once regarded as morphologically and functionally bland, is now recognized to be dynamic, plastic, heterogenous, and involved in a wide array of biological processes including energy homeostasis, glucose and lipid handling, blood pressure control, and host defense.  High fat feeding and other metabolic stressors cause dramatic changes in adipose morphology, physiology, and cellular composition, and alterations in adiposity are associated with insulin resistance, dyslipidemia, and Type 2 diabetes (T2D). Here we provide detailed cellular atlases of human and murine subcutaneous and visceral white fat at single-cell resolution across a range of body weight. We identify subpopulations of adipocytes, adipose stem and progenitor cells (ASPCs), vascular, and immune cells and demonstrate commonalities and differences across species and dietary conditions. We link specific cell types to increased risk of metabolic disease, and we provide an initial blueprint for a comprehensive set of interactions between individual cell types in the adipose niche in leanness and obesity. We also utilize bulk RNA sequencing data of enzymatically digested human adipocytes to scale these findings to a larger cohort. These data comprise an extensive resource for the exploration of genes, traits, and cell types in the function of WAT across species, depots, and nutritional conditions. 

  Study phs002766

• Genomics and Epigenomics of the Elderly Response to Pneumococcal Vaccines

  Pneumococcal Vaccine Study: This is a study designed to develop a detailed transcriptional and cellular profile of the immune response to pneumococcal vaccination with conjugated and non-conjugated polysaccharide vaccines in the senescent immune system of older adults. In this study, 39 healthy adults ages 60 and older that have never received pneumococcal vaccination, have been randomized to receive Prevnar 13 (Pfizer), a conjugated 13-valent vaccine or Pneumovax 23 (Merck), a non-conjugated 23-valent vaccine. Four (4) study visits occurred over 70 days at Day -7 (baseline), day 1, day 10 and day 60 (± 5d) post vaccination.Bulk RNA-sequencing data are available on the 4 time points for 31 of the participants and 1 replicate. ATAC-Doublet Detector Study: We generated single nucleus ATAC-seq (snATAC-seq) data from PBMCs of 2 baseline blood samples from the pneumococcal study. Generated data was then used to benchmark a novel computational method for detecting multiplets in snATAC-seq data, namely, AMULET (https://github.com/UcarLab/AMULET). Multiplets confound downstream analyses by introducing combined epigenomic profiles from 2 or more cells/nuclei. AMULET detects multiplets by using unique snATAC-seq data features that allow it to reliably use read counts to detect multiplets for their removal from data analyses.

  Study phs002361

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC-- CA209-025

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004292

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC -- CA209-010

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004291