7186 results for "({!join from=dataset_studies to=id v=$q1}) OR ("transcriptome" OR "RNA-seq" OR "rna sequencing" OR "rna-sequencing")"

in 52.79 milliseconds.

• 10x genomics single cell RNA-seq and ATAC-seq from 10 pre-menopausal patients fimbriae and ampullary part of the fallopian tube

  The fallopian tube, connecting the uterus with the ovary, is a dynamic organ that undergoes cyclical changes and is the site of several diseases, including serous cancer. Here, we use single-cell technologies to construct a comprehensive cell map of healthy pre-menopausal fallopian tubes, capturing the impact of the menstrual cycle and menopause on different fallopian tube cells at the molecular level. The comparative analysis between pre- and post-menopausal fallopian tubes reveals substantial shifts in cellular abundance and gene expression patterns, highlighting the physiological changes associated with menopause. Further investigations into menstrual cycle phases illuminate distinct molecular states in secretory epithelial cells caused by hormonal fluctuations. The markers we identified characterizing secretory epithelial cells provide a valuable tool for classifying ovarian cancer subtypes.

  Study EGAS50000000628

• RNA-seq data from 224 advanced prostate tumors generated by the West Coast Dream Team

  224 pairs of FASTQ files from metastatic Castration-Resistant Prostate Cancer (mCRPC) sequenced on HiSeq 4000 instruments. Patients were enrolled in the West Coast Dream Team study. Biopsies include various tissue sites including bone, soft tissue, and lymph node.

  Dataset EGAD00001008487

• RNA-seq data from 164 advanced prostate tumors generated by the West Coast Dream Team including 42 paired samles

  164 pairs of FASTQ files from metastatic Castration-Resistant Prostate Cancer (mCRPC) sequenced on HiSeq 4000 instruments. Patients were enrolled in the West Coast Dream Team study. Biopsies include various tissue sites including bone, soft tissue, and lymph node. 42 pairs prior to enzalutamide treatment and at progression from 21 patients are included.

  Dataset EGAD00001009065

• Exome-seq, RNA-Seq, SNP array profiling of gastric tumor samples and cell lines.

  Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we performed exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumors and 32 cell lines. Meta-analysis of the exome data and previously published data sets revealed 24 significantly mutated genes in MSS tumors, and 16 in MSI tumors. More than half of the patients in our collection could potentially benefit from targeted therapies. We identified 55 splice-site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumors. The ZAK kinase isoform TV1 is preferentially upregulated in gastric tumors and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder, and breast cancers. Overexpression of this particular isoform activated multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibited proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

  Study EGAS00001000736

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

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

  Study EGAS00001005295

• Characterizing the Role of the Immune Microenvironment in Multiple Myeloma Progression at a Single Cell Level

  The plasma cell disease multiple myeloma (MM) is highly dependent on its bone marrow microenvironment (BMME). Early alterations within the bone marrow microenvironment that contribute to the progression of multiple myeloma (MM) from its precursor stages could hold the key to identifying novel therapeutic approaches, yet the intrinsic variability in cellular populations between patients together with differences in sample processing and analysis methods have made it difficult to identify consistent changes between datasets. Here, we used single-cell RNA sequencing of BM cells from precursor stages, monoclonal gammopathy of unknown significance (MGUS), smoldering MM (SMM) and newly diagnosed MM and analyzed our data in combination with a previously published dataset that used a similar patient population and sample processing. Despite vast interpatient heterogeneity, some alterations were seen consistently in both datasets. We identified changes in immune cell populations as disease progressed that were characterized by a substantial decrease in memory and naïve CD4 T cells and an increase in CD8+ effector T cells and T-regulatory cells. These alterations were further accompanied by an enrichment of non-clonal memory B cells and an increase of CD14 and CD16 monocytes in MM compared to its precursor stages. These results provide crucial information of the immune changes associated with the progression to clinical MM and can help to develop immune based strategies for patient stratification and early therapeutic intervention. 

  Study phs002756

• Utility of Capillary Blood in Gene Expression Studies

  The FAA Functional Genomics Team examined the feasibility of analyzing gene expression profiles using capillary blood collected by fingerstick, with the goal of developing practical fingerstick blood collections to expand blood collection capabilities and allow large-scale or field-based collections. We compared the RNA sequencing results obtained from two different capillary blood collection and RNA extraction methods with the results obtained from a standard venous blood collection and RNA extraction protocol. We also compared capillary blood processing methods, as well as fingerstick locations, within one capillary collection/extraction method. Within each method we sought to distinguish male from female participants to assess the capacity of each collection/extraction method for gene expression profiling. The study population included 40 participants aged 18-60 who were free from any acute illness or any illness that would preclude blood collection. We determined that capillary blood possesses large-scale significant differences from venous blood and that the two capillary methods produced similar gene expression profiles. We further determined that all methods were useful for detecting gene expression differences between assessed groups (male and female participants in this case). We conclude that capillary blood is suitable for gene expression profiling in quantities as low as 100 μL whole blood, but recommend that researchers should rely on a single blood collection method if they desire to compare results between individual studies. Raw sequencing files (.fastq) from this study are available in dbGaP.

  Study phs003496

• Cell Type-Specific and Disease-Associated eQTL in the Human Lung

  This study was conducted to investigate the genetic control of gene expression in a cell-type and context-specific manner in the human lung. We performed single-cell RNA-sequencing (scRNA-seq) of lung tissue from 66 Interstitial Lung Disease (ILD) cases and 48 unaffected donors. Employing a pseudo-bulk approach, we mapped expression quantitative trait loci (eQTL) across 38 cell types. The scRNA-seq data presented here include previously published (PMID 34262047) and unpublished samples. ILD tissue samples were obtained from lungs removed at the time of lung transplantation at either Vanderbilt University Medical Center (VUMC) or the Norton Thoracic Institute (NTI). Control tissue samples were obtained from lungs declined for organ donation either at the Donor Network of Arizona (DNA) or VUMC. The self-reported ethnicities of donors were 66.7% European, 9.6% African American, 17.5% not available (N/A), and 6.1% other. The smoking history includes 46.5% ever smoker, 29.8% never smoker, 23.7% N/A. For genotyping, flash-frozen tissue in DNA/RNA shield was homogenized via a bullet blender. Genomic DNA was extracted using Zymo Quick DNA/RNA microprep plus kit. Library preparation and low-pass Whole Genome Sequencing (WGS) were carried out at TGen or by Gencove Inc. At TGen, libraries were prepared using PCR-free Watchmaker kits (Watchmaker Genomics) with 200 ng input. Genomes were sequenced on NovaSeq at low coverage (typically 0.4–1x). The resulting sequence data were processed and imputed using Gencove's imputation platform.

  Study phs003521

• Multi-omics bulk and single-cell profiling of epithelioid sarcoma

  Epithelioid sarcoma (EpS) is an aggressive sarcoma, characterized by the loss of SMARCB1 expression. EpS is traditionally classified as distal or proximal according to clinicopathological features, but its molecular characteristics remain largely unknown. To establish an EpS molecular classification and uncover determinants of inter- and intra-patient heterogeneity in EpS, we used multi-omics profiling and integrated the genomic, transcriptional and methylome landscapes with single-cell RNA sequencing on fresh samples as well as spatial transcriptomics. We identified two molecular subtypes of EpS: “distal-like” and “proximal-like”, which were distinct from the histological subtypes. Distal-like tumors include all morphological distal subtypes and a subset of proximal and hybrid EpS. They express a specific molecular single-cell derived epithelial-to-mesenchymal transition signature, which associates with improved patient survival and includes Desmoglein 2 (DSG2), which we identify as a potential routine immunohistochemical biomarker to improve diagnosis. Distal-like EpS also display increased peri-tumoral CD8+ T cell infiltrates and specific tumor-immune cell interactions. Conversely, proximal-like EpS harbor a higher inter-tumoral molecular heterogeneity with some cases resembling other SMARCB1-deficient tumors by DNA methylation profiling, and increased intra-tumoral pro-tumoral macrophage infiltrates. Our study introduces a novel molecular classification of EpS with prognostic value, associated diagnostic biomarkers and distinct tumor immune microenvironments, which extends beyond traditional clinicopathological classification. This paves the way for precision medicine-based therapeutic strategies for patients with EpS.

  Study EGAS50000000973

• Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

  Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

  Study EGAS50000000135

• The Genomic and Transcriptomic Landscape of a HeLa Cell Line

  HeLa is the most widely used model cell line for studying human cellular and molecular biology. To date, no genomic reference for this cell line has been released, and experiments have relied on the human reference genome. Effective design and interpretation of molecular genetic studies performed using HeLa cells require accurate genomic information. Here we present a detailed genomic and transcriptomic characterization of a HeLa cell line. We performed DNA and RNA sequencing of a HeLa Kyoto cell line and analyzed its mutational portfolio and gene expression profile. Segmentation of the genome according to copy number revealed a remarkably high level of aneuploidy and numerous large structural variants at unprecedented resolution. Some of the extensive genomic rearrangements are indicative of catastrophic chromosome shattering, known as chromothripsis. Our analysis of the HeLa gene expression profile revealed that several pathways, including cell cycle and DNA repair, exhibit significantly different expression patterns from those in normal human tissues. Our results provide the first detailed account of genomic variants in the HeLa genome, yielding insight into their impact on gene expression and cellular function as well as their origins. This study underscores the importance of accounting for the strikingly aberrant characteristics of HeLa cells when designing and interpreting experiments, and has implications for the use of HeLa as a model of human biology. Copyright Landry et al., "The Genomic and Transcriptomic Landscape of a HeLa Cell Line," G3: Genes | Genomes | Genetics (2013).

  Study phs000643

• RNA Transcriptomic and DNA Methylation Landscape in FTLD-TDP and Controls

  Frontotemporal lobar degeneration (FTLD) is one of the leading causes of dementia in individuals younger than 65 years, with the aggregation of TDP-43 as one of the most common forms of FTLD. FTLD-TDP is clinically, genetically and pathologically heterogeneous, with GRN and C9orf72 as the most common genetic forms (although more than 50% of the cases are genetically unexplained), and A, B and C as the most common pathological subtypes. To investigate the molecular differences between the different forms of FTLD-TDP, we performed bulk RNA sequencing (RNAseq) and reduced representation bisulfite sequencing (RRBS) from different brain regions of FTLD-TDP patients and controls. The RNAseq dataset includes samples from the frontal cortex of 149 individuals, comprising 127 FTLD-TDP cases (genetically unexplained types A (27), B (20) and C (22), 24 GRN mutation carriers, and 34 C9orf72 expansion carriers) and 22 controls. The RRBS dataset includes samples from the frontal cortex of 158 individuals —119 FTLD-TDP cases (genetically unexplained types A (25), B (24), and C (25), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 39 controls; and from the cerebellum of 151 individuals —113 FTLD-TDP cases (genetically unexplained types A (24), B (21), and C (23), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 38 controls. In total, the study consists of data from 183 unique individuals, with some overlap across the different datasets.

  Study phs004075

• Relapse CHL study

  Classical Hodgkin Lymphoma (CHL) is characterized by a complex tumor microenvironment (TME) that supports disease progression. While immune cell recruitment by Hodgkin and Reed–Sternberg (HRS) cells is well-documented, the role of non-malignant B cells in relapse remains unclear. Using single-cell RNA sequencing (scRNA-seq) on paired diagnostic and relapsed CHL samples, we identified distinct shifts in B cell populations, particularly an enrichment of naïve B cells and a reduction of memory B cells in early-relapse compared to late-relapse and newly diagnosed CHL. The enrichment of naïve B cells in early relapse biopsies was confirmed in independent validation cohorts using scRNA-seq and immunohistochemistry. Notably, naïve B cells in early-relapse samples exhibited high expression of LGALS9, an immunosuppressive molecule known to interact with T cell immunoglobulin and mucin domain 3 (TIM-3) on regulatory T cells (Tregs). Cell-cell interaction analysis revealed the importance of interactions between LGALS9+ naïve B cells and TIM-3+ Tregs in the early-relapse setting. Spatial analysis by imaging mass cytometry confirmed close proximity of LGALS9+ naïve B cells with TIM-3⁺ CD4⁺ T cells and HRS cells, pointing to their role in shaping an immunosuppressive niche. Our findings highlight a previously unrecognized population of LGALS9⁺ naïve B cells with immunoregulatory potential in early-relapse CHL and provide new insights into the spatial and transcriptional architecture of the relapsed TME in CHL.

  Study EGAS00001008222

• Characterization of HCV-specific CD4 T cells during DAA-Therapy

  Background: Chronic HCV-infection is characterized by a severe impairment of HCV-specific CD4 T cell help that is driven by chronic antigen stimulation. We aimed to study the fate of HCV-specific CD4 T cells after viral elimination. Methods: HCV-specific CD4 T cell responses were longitudinally analyzed using MHC class II tetramer-technology, multicolor flow cytometry and RNA sequencing in a cohort of chronically HCV-infected patients undergoing therapy with direct-acting antivirals. In addition, HCV-specific neutralizing antibodies and CXCL13 levels were analyzed. Results: We observed that the frequency of HCV-specific CD4 T cells increased within two weeks after initiation of DAA therapy. Multicolor flow cytometry revealed a downregulation of exhaustion and activation markers and an upregulation of memory-associated markers. While cells with a Th1 phenotype were the predominant subset at baseline, cells with phenotypic and transcriptional characteristics of follicular T helper cells increasingly shaped the circulating HCV-specific CD4 T cell repertoire, suggesting antigen-independent survival of this subset. These changes were accompanied by a decline of HCV-specific neutralizing antibodies and the germinal center activity. Conclusion: We identified a population of HCV-specific CD4 T cells with a follicular T helper cell signature that is maintained after therapy-induced elimination of persistent infection and may constitute an important target population for vaccination efforts to prevent re-infection and immunotherapeutic approaches for persistent viral infections.

  Study EGAS00001003950

• 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

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

• 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

• Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

  Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

  Study EGAS00001003975

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• Assessment of RNA-Seq Sample Preparation Methodology

  The goal of this study was two-fold, to determine if common purification techniques have any effect on downstream differential expression analysis and to evaluate combinations of alignment and differential expression software for reliability. To this end, blood was collected from three individuals and pooled during and after extraction. After pooling and mixing was completed, samples were divided into three aliquots for testing. The first aliquot was a control and not purified or concentrated. It was diluted to produce samples of varying concentrations for testing. The second aliquot was diluted to 20 ng/μL and used to test six different variations on the AMPure XP bead purification procedure. The last aliquot was also diluted to varying concentrations of 60, 30, and 9.6 ng/μL and purified using MinElute columns. Samples were submitted for total RNA-Seq library preparation and sequencing. Library preparation was performed using the TruSeq Stranded Total RNA with Ribo-Zero Globin kit (20020612, Illumina Inc.), and 2x150 bp PE sequencing was done on an Illumina NovaSeq 6000 with the S4 reagent kit. Comparisons were made between methods (AMPure vs. unpurified, AMPure vs. MinElute, MinElute vs. unpurified) to assess the effects of purification method on downstream differential expression. Comparisons were also made within methods using the varying concentrations tested for unpurified samples and for MinElute purified samples to assess the effects of concentration on differential expression. Variations of the AMPure procedure were also compared to assess the effectiveness of the variations tested in comparison to the unmodified procedure. A subset of samples was selected for use with alignment and differential expression package comparison. Unpurified high concentration samples eluted in RNAse-free water were compared to unpurified high concentration samples eluted in BR5, a buffer from the PAXgene blood miRNA kit, with the expectation that there should be no or very few differentially expressed genes. Unpurified high concentration samples eluted in RNAse-free water were also compared to unpurified low concentration samples also eluted in RNAse-free water, with possibly a small number of differentially expressed genes anticipated. A third dataset of simulated RNA-seq data was created with a known rate of differential expression. These files were aligned using Bowtie2, HISAT2, kallisto, RSEM, Rsubread, Salmon, and STAR. Results were then analyzed for differential expression using ALDEx2, baySeq, DEGseq, DESeq2, edgeR, limma, NOISeq, PoissonSeq, and SAMseq. Differential expression results of all three comparisons were evaluated to determine which combinations provided the most reliable results for both real and simulated data.

  Study phs003001

• Single-cell RNA-seq profiling of patient derived organoids

  The aim of this study was to profile developement of resistance to cancer treatments using organoids. Two patient derived colorectal cancer organoids were profiled using single cell RNAseq after treatment with a combination of different anticancer drugs. The organoids were profiled both after 4 weeks of treatment and after a drug holiday of 2 weeks.

  Study EGAS50000001025

• Mutation tracking in single cell RNA-Seq reveals consequences of subclonal evolution in acute myeloid leukemia

  Total bone marrow from AML patients was FACS index sorted with staining for Lineage markers (Lin), CD33, CD34, CD38, CD45RA, CD90, CD99, Tim3 and GPR56, and enrichment for Lin- and Lin-CD34+ cells. Single cells were subjected to a modified smart-seq2 protocol with targeting of mutant sites of interest ("MutaSeq”).

  Study EGAS00001003414

• March 2017 data update (bam/fastq) (containing H3K27ac for CEMT_87, RNA-Seq for CEMT_141, CEMT_142, CEMT_145, CEMT_146, H3K27me3, H3K9me3, H3K4me3, H3K4me1, ChIP-Seq Input for CEMT_88, CEMT_90 and CEMT_91) for reference epigenomes generated at Centre for Epigenome Mapping Technologies, Genome Sciences Center, B.C. Cancer Agency as part of the International Human Epigenome Consortium.

  March 2017 data update (bam/fastq) for reference epigenomes generated at Centre for Epigenome Mapping Technologies (Canadian Epigenetics, Environment and Health Research Consortium), Genome Sciences Center, B.C. Cancer Agency, Vancouver, Canada as part of the International Human Epigenome Consortium.

  Dataset EGAD00001003272

• Gene Expression in an African American Schizophrenia Dataset

  Schizophrenia is a common and severe psychotic disorder. While some common SNPs and rare copy number variants have been identified as being significantly associated with disease risk, the biological mechanisms remain largely undefined. To identify gene expression abnormalities in schizophrenia, we generated whole-genome gene expression profiles using RNAseq on lymphoblastoid cell lines from schizophrenia cases and controls of African American ancestry. Genes related to immune processes were enriched in the genes differentially expressed in the schizophrenia cases. Approximately half of the differentially expressed genes in this African American dataset were also differentially expressed by affection status in our previously studied European ancestry dataset (PMID: 28418402).

  Study phs002842

• RNA sequencing data from glioblastoma primary cell lines treated with indisulam data access committee

  Four glioblastoma primary cell lines (BT592/pt1, BT1007/pt2, BT1009/pt3; BT1012/pt4) were derived from glioblastoma Cavitron Ultrasonic Aspirator (CUSA) material and cultured as neurospheres (GB-NS) in DMEM/F12+Glutamax medium containing B27 supplement (Thermofisher, Waltham, MA, USA) and the mitogenic factors epidermal growth factor (EGF) and fibroblast growth factor b (bFGF) (Peprotech, Rocky Hill, NJ, USA). GB-NS were confirmed mycoplasma-free by PCR Test. Cells were seeded in triplicate in 6-well plates, treated with 5 μM indisulam (Sigma Aldrich, Merk), or vehicle control, and collected after 6 and 16 hours for RNA extraction. Total RNA was extracted using TRIzol Reagent (Life Technologies, Thermo Fisher, cat n. 15596026), RNeasy Mini Kit, and RNase-Free DNase Set (Qiagen, cat n. 74104 and 79254). Sample QC, RNA library preparations, and sequencing reactions were conducted at GENEWIZ (Germany GmbH). The concentration of RNA was quantified using a Qubit Fluorometer (Life Technologies, Carlsbad, CA) and RNA integrity was assayed using a TapeStation. Sequencing libraries were prepared with RNA with PolyA selection and ERCC spike-in with standard ERCC kit (ERCC RNA Spike-In Mix,Thermo Fisher 4456740). Final libraries (Illumina, RNA with PolyA selection) were sequenced on a NovaSeq platform, producing 2×150 base paired-end reads.

  Dac EGAC50000000407

• 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

• Understanding Self- Organising Capacity of Stem Cells during Implantation and Early Post-implantation Development in vitro and in vivo: Implications for Human Development (2020-04-20)

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

  Dataset EGAD00001006056

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

  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 EGAD00001014788

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

  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 EGAD00001014790

• Mapping Cells and Interactions in the Thymus Across Development and Aging (2025-07-28)

  Single cell RNA sequencing on thymic cells from children with syndromic diseases such as Trisomy 21 and inflammatory/autoimmune diseases such as myasthenia gravis. Our group has previously performed atlasing of the healthy human thymus using single cell transcriptomics and obtained an unprecedented understanding of human thymus throughout development. In this study, we intend to survey diseased human thymus and compare this to healthy to obtain a better understanding of thymic dysfunction in human diseases. The diseases we will survey include congenital conditions such as Down's syndrome, CHARGE syndrome, DiGeorge syndrome as well as acquired diseases such as myasthenia gravis and thymic hyperplasia. We intend to perform genetic and genomic study on the samples, including expression analysis. This will not only help us to understand the cells and molecular pathways affected in dysfunctional thymus; but may also lead to generalisable understandings in human genetic diseases and inflammatory/autoimmune diseases. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2025-07-28.

  Dataset EGAD00001015660

• 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

• 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

• Lineage Plasticity and Immune Cell Heterogeneity Are Coordinately Dysregulated Through Changes in FOXA1 Expression in Bladder Cancers with Squamous Differentiation

  We performed whole exome sequencing (WES) of paired, macrodissected regions of urothelial carcinoma (UC), not otherwise specified (NOS-UC) and squamous differentiation (SqD) from 21 bladder tumors in which distinct and separable NOS-UC and SqD regions were present (total 42 exomes, in addition to matched normal exomes). Mean tumor mutational burden (TMB) was higher in the SqD versus the NOS-UC regions. Consistent with the higher TMB, the SqD regions also had a higher neoantigen burden than the patient-matched NOS-UC regions. The SqD regions also exhibited higher karyotypic complexity with a higher median ploidy relative to the paired NOS-UC regions. Phylogenetic analysis of the WES data confirmed that the morphologically distinct regions of all 21 tumors arose from a shared precursor with a median of 68 (range 16 – 852) non-synonymous mutations shared between the NOS-UC and SqD regions, and the overall data were indicative of early branched evolution. As whole exome sequencing did not identify a shared mutational alteration or genomic signature unique to the SqD or NOS-UC components, we therefore performed expression profiling of separable SqD and NOS-UC regions of 12 bladder tumors from the WES cohort with RNA quality sufficient for whole transcriptome RNA-sequencing (RNAseq, total of 24 regions). Based on centroid clustering analysis of expression data, all 12 regions of SqD were basal-squamous subtype based on the TCGA molecular classification schema. Notably, 8 of the NOS-UC regions were also basal-squamous subtype, with only four pairs exhibiting discordant transcriptional subtypes (three NOS-UC regions were luminal-infiltrated, one luminal-papillary). To better quantify differences in the transcriptional profiles of the paired NOS-UC and SqD regions, we used single sample gene set enrichment analysis (ssGSEA) to assess for enrichment of luminal or basal-squamous gene expression. ssGSEA analysis revealed that the SqD regions expressed higher levels of basal-squamous genes based on the BASE47 gene set and UPK/KRT gene sets; whereas the matched NOS-UC regions expressed higher levels of luminal genes. These differences in the levels of expression of basal and luminal genes were observed even in cases in which both the NOS-UC and SqD regions were basal-squamous subtype based the TCGA molecular classification schema. Exomes from eight patients in this study have previously been deposited in dbGaP: accession number phs001783 as part of the study "Exome Recapture and Sequencing of Prospectively Characterized Clinical Specimens From Cancer Patients".

  Study phs002357

• Tracking the evolution of esophageal squamous cell carcinoma under dynamic immune selection by multi-omics sequencing

  Intratumoral heterogeneity (ITH) has been linked to decreased efficacy of clinical treatments. However, although genomic ITH has been characterized in genetic, transcriptomic and epigenetic alterations are hallmarks of esophageal squamous cell carcinoma (ESCC), the extent to which these are heterogeneous in ESCC has not been explored in a unified framework. Further, the extent to which tumor-infiltrated T lymphocytes (TILs) are directed against cancer cells, but how the immune infiltration acts as a selective force to shape the clonal evolution of ESCC is unclear. In this study, we perform multi-omic sequencing on 186 samples from 36 primary ESCC patients. Through multi-omics analyses, it is discovered that genomic, epigenomic, and transcriptomic ITH are underpinned by ongoing chromosomal instability. Based on the RNA-seq data, we observe diverse levels of immune infiltrate across different tumor sites from the same tumor. We reveal genetic mechanisms of neoantigen evasion under distinct selection pressure from the diverse immune microenvironment. Overall, our work offers an avenue of dissecting the complex contribution of the multi-omics level to the ITH in ESCC and thereby enhances the development of clinical therapy.

  Dataset EGAD00001009482

• RNA, ATAC, ChIP datasets from Ankylosing Spondylitis patients and Healthy Controls

  This dataset consists of functional genomic data from 20 Ankylosing Spondylitis patients and 35 Healthy Controls taken from CD4+ T cells, CD8+ T cells and CD14 Monocytes. It contains 364 paired end fastq files consisting of 104 total RNA-seq samples and 116 ATAC samples, for ChIP there are 46 H3K4me3 samples, 46 H3K27ac samples and 3 H3K4me1 samples, along with 49 paired input samples. The samples were sequenced on Illumina HiSeq4000, Illumina NextSeq500 and Illumina NovaSeq 6000 platforms.

  Dataset EGAD00001008758

• Whole genome landscape of adult T-cell leukemia/lymphoma

  Whole-genome sequence (WGS) data of tumor-normal pairs from 139 ATL patients and RNA sequence (RNA-seq) data of tumors from 28 ATL patients.

  Dataset EGAD00001008289

• T-ALL RNA-Seq raw data files

  Bone marrow or peripheral blood samples were collected of adult patients at first diagnosis of T-precursor acute lymphoblastic leukemia. RNA was isolated from mononuclear cells and subjected to mRNA library prep using Poly-A selection and sequencing on a NovaSeq 6000 system. Obtained gene expression profiles and gene fusion calls were used to allocate samples to molecular disease subtypes.

  Study EGAS50000000213

• RNA seq of intestinal biopsies and blood cells from patients with celiac disease and controls

  This study presents RNA sequenced from several intestinal phenotypes such as from patients with damaged intestinal mucosa as well as from those without a damaged intestine but still an autoantibody response to tissue transglutaminase (“potential celiac disease”) or those with treated celiac disease. By combining phenotypes in a much larger number of tissue samples from both peripheral blood and mucosal biopsies, we found highly differentially expressed genes, but also identify genes involved in controlling and triggering celiac disease associated changes. This study identifies molecular mechanisms involved in celiac disease and new possible targets for treatment and for identification of individuals at risk.

  Dataset EGAD50000001597

• Genomic and epigenomic study of Japanese renal cell carcinoma including WGS, RNA-seq, ATAC-seq, and methyl-seq

  We performed whole-genomic and transcriptomic sequencing of more than 100 Japanese renal cell carcinoma (RCC) cases, including clear cell RCC, papillary RCC, chromophobe RCC, and TFE3-translocated RCC. In addition, we conducted epigenome sequencing analyses: the assay for transposase-accessible chromatin sequencing (ATAC-seq) and enzymatic methyl sequencing.

  Study EGAS00001006919

• Decoding Human Heart Morphogenesis through Single-cell Multi-modal Analyses

  The heart, which is the first organ to develop, is highly dependent on its form to function. However, how diverse cardiac cell types spatially coordinate to create complex morphological structures critical for heart function remains to be elucidated. Here, we show that integration of single cell RNA-sequencing with high-resolution multiplexed error-robust fluorescent in situ hybridization (MERFISH) not only resolves the identity of cardiac cell types developing the human heart but also provides a spatial mapping of individual cells that enables illumination of their organization into cellular communities forming distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, supporting their specialization according to cellular ecosystem and anatomic region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities revealed signaling pathways orchestrating the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. In vivo conditional genetic mouse models and in vitro human pluripotent stem cell studies confirmed an intricate multicellular PLXN-SEMA crosstalk among specific ventricular cardiomyocyte, fibroblast and endothelial cell subpopulations that directs the compaction of the ventricular wall layers. Thus, these detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodeling the human heart offer new insights into structural heart diseases as well as engineering complex multi-cellular tissues for human heart repair. Reprinted from Farah et al. (In press), with permission from Nature.

  Study phs002031

• Non-Coding Autoimmune Risk Variant Defines Role for ICOS in T Peripheral Helper Cell Development

  Our fine-mapping of 76 non-MHC loci associated with risk for rheumatoid arthritis (RA, 11,475 cases, 15,870 controls) has recently identified rs117701653, a non-coding single nucleotide polymorphism (SNP) in the CTLA4/CD28/ICOS locus, as the variant most likely to modulate RA risk within this region, with the minor allele C showing disease protection compared to the major allele A. We found that this SNP exhibits allele-specific protein binding and transcriptional modulation, further supporting its regulatory nature. Here, we aimed to determine whether the functional non-coding variant rs117701653 modulates the expression of genes within the CD28/CTLA4/ICOS locus. To address this question, we employed low-input RNA sequencing using the Illumina SmartSeq2 platform on resting total CD4+ T cells from the blood of 24 healthy subjects (8 AA, 8 AC, 8 CC) recruited based on genotype from the Mass General Brigham biobank. We performed a targeted cis-eQTL mapping analysis for protein coding genes with transcription start sites (TSS) within a 1MB window of SNP rs117701653. These genes included WDR12, NBEAL1, CYP20A1, ABI2, RAPH1, CD28, CTLA4, ICOS, and PARD3B. Surprisingly, allelic variation at rs117701653 correlated strongly with expression of ICOS. This suggests that individuals carrying the protective allele C at rs117701653 express lower levels of ICOS, implicating ICOS as a plausible target gene influenced by this protective allele. The raw sequence data from the 24 healthy individuals are available through dbGAP. Gene expression matrices have been deposited to GEO accession GSE235868.

  Study phs003448

• Surgery in combination with immune checkpoint therapy as an effective treatment for patients with metastatic cancer.

  Abstract Surgical removal of primary tumors was shown to reverse tumor-mediated immune suppression in pre-clinical models with metastatic disease. However, how cytoreductive surgery in the metastatic setting modulates the immune responses in patients, especially in the context of immune checkpoint therapy (ICT)-containing treatments is not understood. Here, we report the first prospective, non-comparative clinical trial (N=104) using three different ICT-containing strategies plus cytoreductive or “debulking” surgery to remove the primary tumor-bearing kidney or a metastatic lesion as a treatment for patients (N=43) with metastatic clear cell renal cell carcinoma (mccRCC). For those patients (N=61) who were not candidates for cytoreductive surgery, a biopsy was obtained instead for correlative biological studies. Our data demonstrated that the combination of ICT with cytoreductive surgery was safe and feasible in patients with mccRCC. The 2-year overall survival was 84% with a median OS of 54.7 months for patients who received ICT containing regimens plus surgery. Immune-monitoring studies with co-detection by indexing (CODEX) identified distinct tumor spatial conformation of cellular subsets as a novel and improved predictor of response to ICT. Importantly, single-cell RNA-sequencing data demonstrated that surgical removal of the tumor increased antigen-presenting dendritic cell population with a concurrent reduction in KDM6B expressing immune-suppressive myeloid cells in the peripheral blood. Together, this study highlighted the feasibility of combining ICT with cytoreductive surgery in a metastatic setting and demonstrated potential enhancement of immune responses following ICT plus cytoreductive surgery in patients with metastatic disease.

  Study EGAS00001005667

• DNA demethylation is associated with malignant progression of low-grade gliomas

  To elucidate the mechanisms of tumor progression is key to discovering strategies for treating patients with lower-grade glioma. Thus to characterize the changes in the molecular profile during malignant progression, we performed methylation array analysis (122 tumors), whole-exome sequencing (36 tumors and matched normal samples), and RNA sequencing (31 tumors). This cohort included 24 matched pairs of initial lower-grade gliomas and recurrent tumors, most of which showed malignant progression when they recurred. Nearly half of IDH-mutant glioblastomas that had progressed from lower-grade gliomas exhibited characteristic partial DNA demethylation in previously methylated genomic regions of their corresponding initial tumors, which had the glioma CpG island methylator phenotype (G-CIMP). In these glioblastomas, cell cycle-related genes were upregulated, RB and PI3K-AKT pathway genes were frequently altered. Notably, bioinformatics analysis revealed that late-replicating domain was significantly enriched in the demethylated regions that were mostly located in non-regulatory regions, suggesting that the loss of DNA methylation during malignant transformation may involve mainly passive demethylation due to a delay in maintenance of methylation during accelerated cell division. Nonetheless, a limited number of genes including IGF2BP3, which potentially drives cell proliferation, were presumed to be upregulated due to demethylation of their promoter. Thus, cell proliferation was further promoted by passive demethylation in the promoter region of these oncogenic genes. Our data indicated that demethylation of the G-CIMP profile found in a subset of recurrent gliomas reflects accelerated cell divisions accompanied by malignant transformation. Oncogenic genes activated by such epigenetic change represent potential therapeutic targets.

  Study JGAS000146

• 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

• Colorectal cancer organoids expressing BRAF (fusion) genes

  Fusion genes may be oncogenic drivers and potential targets for personalized therapies in a variety of cancer types. The BRAF gene is frequently involved in oncogenic fusions, with fusion frequencies of 0.2-3% throughout different cancers. BRAF can be fused to a wide variety of genes, however, BRAF fusions rarely occur in the exact same gene configuration, making assessment of the fusion relevance and decision for drug treatment challenging. Here, we devised a colorectal cancer organoid-based platform to functionally characterize diverse BRAF fusions, containing various partner genes (AGAP3, DLG1 and TRIM24). We compared these BRAF fusions with respect to cellular behaviour, downstream signaling activation and response to targeted therapies. We found that 5’ partner choice of BRAF affects cellular localization and intracellular signaling capacities of the fusion genes. The DLG1-BRAF fusion gene showed distinct localization to the plasma membrane and exhibited increased levels of MAPK pathway activation under unperturbed conditions. Furthermore, the different BRAF fusions showed varying sensitivities to the targeted inhibition of ERK and BRAF, with the DLG1-BRAF fusion being the most sensitive. Additionally, RNA-sequencing identified distinct subsets of genes affected by the DLG1-BRAF fusion gene. Importantly, all fusion genes conveyed resistance to the clinically relevant EGFR/HER2/HER4-inhibitor afatinib, suggesting that BRAF fusions should be screened alongside other MAPK pathway alterations to identify mCRC patients amenable to cetuximab treatment. In summary, we developed a platform to efficiently assess molecular and cellular effects of fusion genes, and revealed that differential drug responses and distinct gene expression profiles can be induced by different BRAF fusions .

  Study EGAS00001003558

• Healthy_ageing_thymus

  Single cell RNA sequencing of immune and non-immune cells from healthy ageing thymus Thymus is a primary lymphoid organ that creates an environment for the T lymphocyte precursors differentiation into the naive T cells, which involves TCR gene rearrangement, negative and positive selection. Despite thymus critical function for the recognition of the pathogens it starts to atrophy very early in life. Thymic Epithelial Space (TES), where thymocyte education occurs, starts to decline rapidly just after birth, while with the onset of puberty thymus is also progressively replaced with adipose tissue. This leads to a reduced naive T cell output,which decreases an organism's ability to recognize pathogens. In our lab, we have already collected and characterized cell type composition and gene regulatory networks of the fetal tand adult human thymus (Park et al, Science 2020). In the current project we plan to expand the number of paediatric and adult samples in the study with the aim to understand how the process of thymic involution happens. We will use single-cell transcriptomics to zoom on changes that occur in immune and non-immune cell types in different phases of thymic involution with the hope to understand which gene expression changes are causal for the thymic involution process and how naïve T cell production can be increased in elderly. 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 EGAS00001004311

• Transcriptomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum

  Background; Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH.MethodsMultiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping were performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors).ResultsOverall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high ITH and were classified discordantly. Proximal tumors showed ITH related to differences in the microenvironment, while ITH within distal tumors were cancer-cell related. This subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05).ConclusionOur transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.

  Study EGAS00001004668

• A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer

  Background Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. Methods We combined immunohistochemistry (IHC) with RNA-sequencing and whole genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Further, we explored transcriptional features of immune archetypes, T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions, and assessed common adult biomarkers such as neoantigen load and TMB. Results A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. Conclusions Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

  Study EGAS00001007029

• Pharmacogenomic Analysis Reveals New Therapeutic Options for Pleural Mesothelioma

  Background: Pleural mesothelioma (PM) urgently requires novel effective treatments. Drug repurposing is a relevant strategy to discover new therapeutic options among well-known drugs used in clinical practice. This study aimed to identify new therapies in the context of molecular heterogeneity of PM. Methods: We performed a multi-omics study of the largest cohort of patient-derived primary PM cell lines (n=58) to date, using whole exome sequencing, RNA sequencing and methylome analysis. We also conducted a large-scale drug screening of over 1300 compounds on a subset of cell lines (n=11), as well as an extension screen (n=48) for selected drugs of interest. The best candidates were validated using 3D culture and mouse models. Results: We demonstrated that the molecular landscape of the cell lines cohort recapitulates the main findings in tumours and revealed new specific features of PM. The large-scale drug screening and subsequent validations highlighted notably strong sensitivity to entinostat and fluvastatin. Importantly, both compounds displayed same efficacy in 2D and 3D culture models and single treatment prolonged survival in vivo similarly or even better than standard chemotherapy. Strikingly, combining entinostat treatment with anti-PD1 led to cures in some mice and immunized them against tumour cells. Conclusions: Among several candidates, we identified entinostat, notably in combination with immune checkpoint inhibitor, and to a lesser extent fluvastatin, as promising therapeutic options for patients. We believe our study provides a valuable resource that could contribute to advances in the field and warrants further clinical investigation.

  Study EGAS00001007866

• Molecular Analysis of Short- Versus Long-Term Survivors of High-Grade Serous Ovarian Carcinoma

  The underlying determinants for long-term survival (LTS, ≥10 years) versus short-term survival (STS, <3 years) of patients diagnosed with high-grade serous ovarian carcinoma (HGSC) are largely unknown. The present study sought to identify molecular predictors of LTS for women with HGSC. A cohort of 24 frozen HGSC samples (12 LTS and 12 STS) were collected and analyzed at DNA and RNA levels. Copy number variation and mutation did not indicate large differences between LTS and STS. From RNA-Seq analysis, 11 genes were found to be differentially expressed between the STS and LTS groups (fold change >2; false discovery rate <0.01). In the validation cohort, transmembrane protein 62 (TMEM62) was found to be related to LTS. Our in vitro and in vivo studies demonstrate that restoring TMEM62 may be a novel approach for treatment of HGSC. These findings may have implications for biomarker and intervention strategies to help improve patient outcomes.

  Study phs003020

• Comprehensive Omics Analysis of Pediatric Solid Tumors and Establishment of a Repository for Related Biologic Studies (10C0086)

  This study (10-C-0086) collected and analyzed tumor and circulating tumor omics of pediatric and young adult patients with relapsed rhabdomyosarcoma who were co-enrolled in and receiving treatment on the interventional study 17-C-0049. We sought to describe the summary genomic findings of tumors, and report on the detection of circulating tumor DNA (ctDNA) in serial samples. Tumor tissue and blood were collected and analyzed. Tumor samples were subjected to whole exome and/or whole genome sequencing paired with RNA-Seq. Cell free DNA (cfDNA) was extracted from serially collected blood samples and sequenced. Fusion status, mutations, and IGF1R and YES1 expression from tumor samples revealed the presence of a PAX3 fusion in most samples, a number of common mutations, and universal but heterogeneous expression of IGF-1R and YES1. ctDNA was detected above the 3% threshold in a majority of patients and analysis of cfDNA demonstrated an ability to monitor tumor clonal evolution.

  Study phs003243

• Single-cell Profiling of an In Vitro Model of Human Interneuron Development Reveals Temporal Dynamics of Cell Type Production and Maturation

  GABAergic interneurons are essential for neural circuit function and their loss or dysfunction is implicated in human neuropsychiatric disease. In vitro methods for interneuron generation hold promise for studying human cellular and functional properties and ultimately therapeutic cell replacement. We used a protocol for generating cortical interneurons from hESCs and analyzed the properties and maturation timecourse of cell types using single-cell RNAseq (data available at GEO under: GSE93802 on March 10 2017). Transcriptomic profiles of the hESC-derived interneurons were compared to several different populations of cells from mid-gestation human neocortex that showed differing levels of PAX6 and SOX2 expression. For this study, 104 samples of 100 human neocortical cells each, have been sorted based on SOX2 and PAX6 expression, mRNA recovered from the fixed cells by FRISCR, and transcriptomic profiles generated by SmartSeq2. The RNA-seq data from the 104 100-cell samples is included in this dbGaP study.

  Study phs001276

• Functional Dynamics of the Elderly Gut Microbiome During Probiotic Consumption

  A mechanistic understanding of the health benefits conferred by consumption of probiotic bacteria has been limited by our knowledge of the resident gut microbiota and its interaction with the host. We used fecal samples from a study of twelve healthy individuals aged 65-80 to characterize the structure and functional dynamics of the gut microbiota associated with consumption of the single-organism probiotic, Lactobacillus rhamnosus GG ATCC 53103 (LGG). Samples were collected prior to probiotic consumption (day 0), on day 28 immediately after consuming 10^10 CFU of LGG twice daily for 28 days and day 56, one month after stopping LGG consumption. Our integrative approach incorporated bacterial 16S rRNA gene sequencing, whole-community expression profiling using RNA-seq, and metagenomic sequencing. We highlight the value of combinatorial 'omics methods and concomitant high-resolution informatics to probe the role that probiotics may play on the structure and function of the resident microbiota.

  Study phs000896

• Integrative Age-Related Changes in Genome and Epigenome in Human Lung in Relation to Smoking

  The Integrative Age-Related Changes in Genome and Epigenome in Human Lung in Relation to Smoking study was designed to assess DNA methylation and gene expression profiles in cytologically normal bronchial progenitor "basal" cells. The cells were collected from 54 consented bronchoscopy specimen donors at standard airway locations during clinically indicated bronchoscopy, and then outgrown in primary cultures. DNA methylation profiles were assessed using Enzymatic-Methyl sequencing (EM-seq: 96% genomic coverage, 7.5x depth) in all 54 donors' samples. The gene expression profile was evaluated by bulk RNA-Seq in 39 donors' samples. The study participants were at the ages of 33 - 81 years old at enrollment. Of the 54 donors, 31 were newly diagnosed with lung cancer. As for tobacco smoke exposure, there were 19 current, 18, former, and 17 never smokers. There are 26 subjects whose specimens were used for single cell whole-genome sequencing study phs002758.

  Study phs003317

• Multi-Ethnic Study of Atherosclerosis (BioLINCC)

  The Multi-Ethnic Study of Atherosclerosis (MESA) is a study of the characteristics of subclinical cardiovascular disease (disease detected non-invasively before it has produced clinical signs and symptoms) and the risk factors that predict progression to clinically overt cardiovascular disease or progression of the subclinical disease. MESA researchers study a diverse, population-based sample of 6,814 asymptomatic men and women aged 45-84. Thirty-eight percent of the recruited participants are white, 28 percent African-American, 22 percent Hispanic, and 12 percent Asian, predominantly of Chinese descent. Data available for request include phenotypic data previously released on BioLINCC as well as raw images (E.g. Echocardiogram) and raw tracing data files (e.g. Electrocardiogram). These new data are intended to be made available to interested researchers via BioData Catalyst, will include participants irrespective of participation in genetics (i.e. the full MESA cohort), and may be merged with currently posted (to dbGaP) phenotype and molecular datasets (WGS, GWA, RNA Seq, Methylation, Metabolomics, and Proteomics).

  Study phs003288

• Integrative Analysis for Multi-Omics Data in Non-Small-Cell Lung Cancer

  We conducted epigenetic and transcriptomic mapping of primary non-small cell lung tumors and non-neoplastic tissues. 20 non-neoplastic lung tissues and 18 non-small cell lung tumors were obtained from Biorepository and Tissue Technology Shared Resource at Moores Cancer Center, UC San Diego. All tumor samples were in pathological stages I and II. 8 pairs of samples were "matched", meaning that the tumor and normal tissues came from the same patient. Microfluidic oscillatory washing–based chromatin immunoprecipitation followed by sequencing (MOWChIP-seq) was used to profile the binding of H3K4me1, H3K4me3, H3K9me3, H3K27ac and H3K27me3. SmartSeq-2 was used to map mRNA expression. We identified a large number of differentially modified regions for both epigenetic and transcriptomic markers between tumors and normal tissues. Deep learning model showed extensive transcription factor network rewiring. Raw sequencing data of the ChIP-seq and RNA-seq libraries are available through this dbGaP submission.

  Study phs003113

• Genome sequencing identifies splice-disrupting variants in childhood heart disease

  We used genome sequencing (GS) to identify cardiac specific splice disrupting variants that were subsequently validated in a subset of participants. In turn, we developed a heart-specific model for canonical and non-canonical splice variants, which was applied to patients with congenital heart disease (CHD). These included patients with two of the most common forms of cyanotic CHD i.e. tetralogy of Fallot (TOF) and dextro-transposition of the great arteries (TGA). In addition to identifying canonical splice-disrupting variants in known CHD-related genes in 1% cases, this approach identified putatively damaging non-canonical splice-disrupting variants in 11% of isolated CHD, with deeply intronic variants representing 53% of non-canonical splice-disrupting variants in CHD genes. GS was critical for the identification of variants that would not be captured by routine clinical genetic tests including exome sequencing, while cardiac RNA-Seq allowed for high specificity in the interpretation of splice-disrupting effects in the heart.

  Study EGAS50000000586

• Dual inhibition of FLT3 and BCL-2 is effective in preclinical models of BCL11B-activated lineage ambiguous leukemia

  Aberrant activation of BCL11B (“BCL11B-a”) defines a subtype of lineage ambiguous leukemias with T-lymphoid and myeloid features, co-occurring activating FLT3 mutations, and a stem/progenitor immunophenotype and gene expression profile. As with other lineage ambiguous leukemias, optimal treatment is unclear and there are limited targeted therapeutic options. Here, we investigated the efficacy of BCL-2 and FLT3 inhibition with venetoclax and gilteritinib, respectively, in preclinical models of BCL11B-a leukemia. Despite variation in response to single agent therapies, the combination of venetoclax plus gilteritinib was highly effective in all models evaluated. BH3 profiling and single cell RNA-seq analysis suggests that resistance to venetoclax and/or gilteritinib monotherapies is due to the tumor-intrinsic activity of BCL-xL and/or MCL-1 prior to drug treatment. These data support clinical evaluation of venetoclax in combination with gilteritinib in BCL11B-a lineage ambiguous leukemias.

  Study EGAS50000000978

• Inferring expressed genes by whole-genome sequencing of plasma DNA

  The analysis of cell-free DNA (cfDNA) in plasma represents a rapidly advancing field in medicine. cfDNA consists predominantly of nucleosome-protected DNA shed into the bloodstream by cells undergoing apoptosis. We performed whole-genome sequencing (WGS) of plasma DNA and identified two discrete regions at transcription start sites (TSS) where the nucleosome occupancy results in different read-depth coverage patterns in expressed and silent genes. By employing machine learning for gene classification we found that the plasma DNA read depth patterns from healthy donors reflected the expression signature of hematopoietic cells. In cancer patients with metastatic disease, we were able to classify expressed cancer driver genes in regions with somatic copy number gains with high accuracy. We could even determine the expressed isoform of genes with several TSSs as confirmed by RNA-Seq of the matching primary tumor. Our analyses provide functional information about the cells releasing their DNA into the circulation.

  Study EGAS00001001754

• Systematic kinase inhibitor profiling identifies CDK9 as a synthetic lethal target in NUT midline carcinoma

  Kinase inhibitors represent the backbone of targeted cancer therapy, yet only a limited number of oncogenic drivers are directly druggable. By interrogating the activity of 1,505 kinase inhibitors we found that BRD4-NUT-rearranged NUT midline carcinoma (NMC) cells are specifically killed by CDK9 inhibition (CDK9i) and exhibit a unique dependency on CDK9 and Cyclin-T1 expression. We show that CDK9i leads to robust induction of apoptosis and of markers of DNA damage response in NMC cells. While both CDK9i and bromodomain inhibition over time result in reduced Myc protein expression, only bromodomain inhibition induces cell differentiation and a p21-induced cell cycle arrest in these cells. Finally, RNA-Seq and ChIP-based analyses reveal a BRD4-NUT specific CDK9i-induced perturbation of transcriptional elongation. Thus, our data provide a mechanistic basis for the genotype-dependent vulnerability of NUT midline carcinoma cells to CDK9i that may be of relevance for the development of targeted therapies for NMC patients.

  Study EGAS00001002588

• Tissue and plasma RNA from esophageal cancer and precursor lesions

  In the past decades, the incidence of esophageal adenocarcinoma has increased dramatically in Western populations. Better understanding of disease etiology along with the identification of novel prognostic and predictive biomarkers are urgently needed to improve the dismal survival probabilities. Here, we performed comprehensive RNA (both coding and non-coding) profiling in various samples from 17 patients diagnosed with esophageal adenocarcinoma, high-grade dysplastic or non-dysplastic Barrett’s esophagus. Per patient, a blood plasma sample, and a healthy esophageal and disease tissue sample were included. In total, this comprehensive dataset consists of 102 RNA-seq libraries from 51 samples. This unique resource allows for discovery of novel biomarkers and disease mechanisms, comparison of tissue and liquid biopsy profiles, integration of coding and non-coding RNA patterns, and can serve as a validation dataset in other RNA landscaping studies. Moreover, structural RNA differences can be identified in this dataset, including protein coding mutations fusion genes, and circular RNAs.

  Study EGAS00001004939

• Oncogenic cooperation in a human de novo T-ALL model

  We previously modelled T-cell acute lymphoblastic leukemia (T-ALL) by over-expressing four oncogenes (activated NOTCH1, LMO2, TAL1, and BMI1) in healthy human umbilical cord blood cells through lentiviral transduction (Kusakabe et al, Nat Commun 10, 2913 (2019). https://doi.org/10.1038/s41467-019-10510-8). We have now tested individual effects of the four oncogenes by transducing human cells with all possible combinations of one, two, three, and four oncogenes. We also tested the effects of additional known T-ALL oncogenes including HOXA9, LYL1, MEF2C, NXK2.1, TLX1, and TLX3. In addition to culture in vitro on OP9-DL1 feeders, we also injected transduced cells into immunodeficient NSG mice. We observed differing abilities of the various oncogene combinations to support growth in vitro and frank leukemia in vivo. This study includes RNA-seq data from samples of transduced cells cultured in vitro as well as cells harvested from clinically moribund mice after development of leukemia in vivo.

  Study EGAS00001006055

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

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

  Study EGAS00001007590

• 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

• RNA-Sequencing of cervical cancers

  Question: Can we identify cervical cancer patients who are at risk for distant metastatic (DM) recurrence following treatment with radiotherapy, concurrent weekly cisplatin and brachytherapy (RTCT). Findings: An immune-based 55 gene risk score was developed using a cohort of 81 patients treated with RTCT that was strongly predictive of DM and cause-specific survival (CSS). The risk score was validated in two independent patient cohorts. A high immune metastatic risk score was associated with a high tumor mutational burden and a ‘cold’, immune-excluded tumor microenvironment at diagnosis. Meaning: The immune gene expression risk score may help to identify patients at risk of DM and potential targets for mitigating this risk.

  Study EGAS50000000087

• Developing Biomarkers Incorporating High Throughput RNA, DNA, Small RNA Sequencing and Protein Expression in Inflammatory Bowel Disease Using Formalin-Fixed, Paraffin-Embedded (FFPE) Tissue Samples

  This study included pediatric patients with suspected inflammatory bowel disease (IBD) who underwent endoscopy between 2008 and 2012. This included 101 colonic tissue samples and 100 ileal tissue samples. Biopsies were obtained at the time of diagnosis, and patients were followed for a mean of 6 years to record clinical characteristics. These included baseline characteristics, such as age of diagnosis and initial Montreal classification, as well as outcomes, such as stricturing or fistulizing complications, disease remission, or progression to surgery.Genes with significant results in the non-IBD vs CD comparison in colonic tissue are listed in the "Links to Related Genes" section of this page.

  Study phs003156

• Single-cell RNA-sequencing of CD34+ bone marrow cells from patients with SLE, healthy individuals and umbilical cord blood samples

  All immune cells that contribute to the pathogenesis of systemic lupus erythematosus (SLE) derive from adult haematopoietic stem and progenitor cells (HSPCs) within the bone marrow (BM). We reasoned that the fundamental abnormalities in the disease can be traced back to HSPCs. The aim of this study was to delineate human haematopoietic CD34+ progenitor subpopulations and to identify deregulated pathways in each subpopulation in SLE as compared to healthy controls and umbilical cord blood samples. This study revealed both quantitative-as evidenced by decreased numbers of non-proliferating early progenitors and qualitative differences- characterized by an IFN signature, which is known to drive loss of function and depletion of HSPCs.

  Study EGAS00001007317

• 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

• 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

• Ewings Sarcoma RNA-Seq

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003062

• Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence (RNA-seq)

  Understanding how prostate cancer cells adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenetic plasticity towards pro-survival signaling, and uncovered circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel clinical lead for therapeutic development.

  Study EGAS00001006016

• Mapping regulatory variation in sensory neurons using IPS lines from the HIPSCI project

  Sensory neurons are nerve cells that are activated by sensory input such as heat, light and convey information to the brain. Although a key cell type in complex organisms, human sensory neurons are challenging to study because they are impossible to obtain from living donors. We have collaborated with the Neucentis Pharmaceutical Research Unit to differentiate sensory neuron like cells from human induced pluripotent stem cells derived as part of the Human Induced Pluripotent Stem Cells Initiative. We will sequence RNA from 100 IPS lines derived from healthy individuals and perform RNA-seq on the differentiated cells to identify noncoding variants that alter gene expression in human sensory neurons.

  Dataset EGAD00001003145

• Whole genome sequencing of multiple myeloma patient samples.

  Tumor (CD138+ plasma cells) and non-tumor (peripheral blood white cells or stem cell harvest) DNA from patients with a plasma cell dyscrasias were sequenced. Whole genome sequencing using high molecular weight DNA was performed using the 10X Genomics Chromium platform on either a HiSeq4000 or NovaSeq (Illumina) using 100 to 150 bp paired-end reads. The dataset consists of 111 patients, and 223 samples in total (matched tumor and control per patient; one patient had 2 tumor samples sequenced). Diseases consisted of 2 MGUS patients, 8 SMM patients, 91 newly diagnosed myeloma patients, 1 previously treated patient, 4 relapsed MM patients, and 5 PCL patients. Paired RNA-seq data are available for 81 of the samples under study EGAS00001003411.

  Dataset EGAD00001004452

• 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

• Resistance to Checkpoint Blockade Therapy Through Inactivation of Antigen Presentation

  Although immune checkpoint blockade (CPB) leads to prolonged responses in 15-40% of patients with metastatic melanoma, treatment refractory disease and progression after initial response remain major causes of mortality. While predictors of response have been reported, the common mechanisms of both primary and acquired resistance are poorly understood. To identify mechanisms of resistance and examine the evolving landscape in response to CPB, we performed whole exome sequencing (WES), immunohistochemistry (IHC), and RNA-sequencing (RNAseq) of longitudinal tumor biopsies from 17 metastatic melanoma patients treated with various CPB therapies. We found no significant changes in both mutational and neoantigen loads over time between responders and nonresponders. However, we identified abnormalities in one gene, beta-2-microglobulin (B2M), an essential component of MHC Class I antigen presentation, that were present in samples during disease progression but not regression. In total, we identified B2M aberrations in 29.4% of patients, including multiple early frameshift mutations, loss of heterozygosity (LOH) overlapping B2M, and absence of tumor-specific B2M protein expression. Additional defects in the antigen presentation and IFNγ pathways were identified but were not restricted to progressing lesions in our cohort. In two independent cohorts of 105 and 38 melanoma patients treated with ipilimumab (anti-CTLA4) and pembrolizumab (anti-PD1) respectively, we found that B2M LOH was enriched 3-fold in nonresponders (~30%) vs. responders (~10%) and associated with poorer overall survival (log-rank p=0.01, p=0.006). Loss of both copies of B2M was found only in nonresponders. We also found evidence for association of LOH overlapping IFNGR1 with poorer overall survival exclusively in the anti-PD1 cohort. Thus, B2M loss is likely a common mechanism of primary and acquired resistance to therapies targeting CTLA4 or PD-1.

  Study phs001427

• AC-ICAM: An Atlas and Compass of Immune-CAncer-Microbiome Interactions in Colon Cancer

  Fresh frozen tumor samples and matched normal colon specimens from chemotherapy naïve, consecutive patients with a histological diagnosis of colon adenocarcinoma at Leiden University Medical Center (LUMC) were subjected to multi-omics profiling at Sidra Integrated Genomics Services, part of Sidra Medicine, Doha, Qatar. This resulted in a high-quality dataset of samples collected between 2001 and 2015. Multi-step QC was performed, including assessment of the quality of the nucleic acid extraction, sequencing, and genetic concordance between DNA/RNA sequences from the same patient. After QC, this resulted in RNA sequencing on 348 tumor samples, whole exome sequencing (WES) on 281 tumor-normal sample pairs (n=562 samples), TCR sequencing of the TCR β locus (TRB) for 114 tumor and 10 matched normal samples, and 16S rDNA sequencing on 246 tumor-normal sample pairs (n=492 samples) and 42 tumor only samples, and 167 WGS non human-aligned reads of the tumor samples. Extensive clinicopathological and survival data of the cohort are available. The median follow-up time for this cohort is 4.6 years, and the distribution of gender and pathological stage are within expected ranges for this patient population with a mean age of 68. We refer to this resource as the Sidra-LUMC (SILU) colon adenocarcinoma (COAD) Atlas and Compass of Immune-Colon cancer-Microbiome interactions (AC-ICAM)."Reprinted from Nature Medicine,Roelands, J., Kuppen, P.J.K., Ahmed, E.I. et al. An integrated tumor, immune and microbiome atlas of colon cancer. Nat Med (2023). https://doi.org/10.1038/s41591-023-02324-5, with permission from Nature."

  Study phs002978

• A Phase I Study With a Personalized Neoantigen Cancer Vaccine in Glioblastoma Multiforme

  Neoantigens, which are derived from tumor-specific protein-coding gene mutations, are exempt from central tolerance, can generate robust immune responses and function as bona fide antigens that facilitate tumor rejection. A strategy of using multi-epitope and personalized neoantigen vaccination has previously been tested in patients with high-risk melanoma. We demonstrated that this strategy is feasible for glioblastoma (GBM) which typically has a relatively low mutation load and an immunologically 'cold' tumor microenvironment. We conducted whole exome sequencing (WES) of tumor and normal cells from individual GBM patients to identify tumor-specific mutations. We assessed the expression of mutated alleles by RNA-sequencing of tumor, and used personalized neoantigen-targeting vaccines to immunize newly diagnosed GBM patients following surgical resection and conventional radiotherapy in a phase I/Ib study. Patients who did not receive dexamethasone, a highly potent corticosteroid that is frequently prescribed to treat cerebral oedema in GBM patients, generated circulating polyfunctional neoantigen-specific CD4+ and CD8+ T cell responses that were enriched in a memory phenotype and showed an increase in the number of tumor-infiltrating T cells. Using single-cell T cell receptor analysis, we provided an evidence that neoantigen-specific T cells from the peripheral blood can migrate into an intracranial GBM tumor. Neoantigen-targeting vaccines thus have the potential to favorably alter the immune environment of GBM.In order to discover cancer antigens derived from annotated and unannotated protein-coding regions of the genome, we carried out matched ribosome profiling (Ribo-Seq) on two GBM samples. We discovered novel or unannotated open reading frames (nuORFs) and their expression levels as well as somatic mutations within the nuORFs, which could be used to predict potential neoantigens.

  Study phs001519

• Investigating Delayed-Onset Drug Hypersensitivity Reactions Prospectively

  Delayed-type drug hypersensitivity reactions (dtDHR) are significant causes of morbidity and mortality, but the origin, phenotype and function of pathogenic T cells across the spectrum of severity are unknown. Research has been hampered by several barriers. We leveraged recent technical advancements to study skin-resident memory T cells (TRM) versus recruited T cell subsets in the pathogenesis of severe systemic forms of disease, Stevens-Johnson syndrome/toxic epidermal necrolysis (SJS/TEN) and drug reaction with eosinophilia and systemic symptoms (DRESS), and skin-limited disease, morbilliform drug eruption (MDE). As part of the overarching study, we used single cell RNA sequencing (scRNA-seq) with CITE sequencing and TCR sequencing on prospectively collected skin and blood samples from patients with SJS/TEN and MDE patients. We compared this to healthy controls obtained from discarded skin from surgeries and discarded blood collars from a blood donation center. These data are included here in dbGaP database. The overarching study also performed bulk transcriptional profiling and microscopic analysis of clinical skin samples, as well as a mouse model of MDE (PMID: 39042477). Together, the data supported clonal expansion and recruitment of cytotoxic CD8+ T cells (identified using the mean normalized expression of NKG7, GNLY, GZMA, GZMB, and PRF1) from circulation into skin, along with both expanded and non-expanded cytotoxic CD8+ skin TRM in SJS/TEN whereas MDE displayed a cytotoxic T cell profile without expansion or recruitment of cytotoxic CD8+ T cells from circulation implicating TRM as the potential cause of disease. Concurrently, MDE displayed a pro-Treg signature, highlighted by CCR8 and TNFRSF4 (OX40), which promotes Treg survival compared to SJS/TEN.

  Study phs003344

• Analysis of CD20 loss in patients treated with Mosunetuzumab

  CD20 is an established therapeutic target in B-cell malignancies. The CD20xCD3 bispecific antibody mosunetuzumab has significant efficacy in B-cell non-Hodgkin lymphomas (NHL). Since target antigen loss is a recognized mechanism of resistance, we evaluated the contribution of CD20 expression to clinical response in patients with relapsed and/or refractory NHL in the phase 1/2 GO29781 trial (ClinicalTrials.gov: NCT02500407) investigating mosunetuzumab monotherapy. CD20 was studied using immunohistochemistry (IHC), RNA sequencing, and whole-exome sequencing performed centrally in biopsies collected prior to treatment, or as paired samples at pre-dose and on-treatment, or at-progression. Prior to treatment, most patients exhibited a high proportion of tumor cells expressing CD20; however, the proportion of CD20+ tumor cells in 14/293 (6%) patients was <10%. Analyses of paired biopsies from patients during treatment revealed that CD20 levels were maintained in 29/30 (97%), vs patients with samples at-progression, in which 11/32 (34%) showed CD20 loss. Reduced transcription or acquisition of truncating mutations explained most but not all cases of CD20 loss. In vitro modeling confirmed the effect of CD20 variants identified in clinical samples on reduction of CD20 expression and confirmed missense mutations in the extracellular domain that could block mosunetuzumab binding. This study expands knowledge about the occurrence of target antigen loss after anti-CD20 therapeutics to CD20-targeting bispecific antibodies and elucidates the underlying mechanisms for reduced CD20 expression at disease progression that may be generalizable to other anti-CD20 targeting agents. These results also establish the utility of readily-available IHC staining for CD20 as a tool to inform clinical decisions.

  Study EGAS50000000151

• ATAC-seq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001438

• 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

• Stromal-induced epithelial-mesenchymal transition induces drug resistance in acute lymphoblastic leukemia

  Bone marrow (BM) microenvironment-induced drug resistance in acute lymphoblastic leukemia (ALL) is thought to be mediated by intercellular interaction of leukemic cells with BM niches, but the mechanistic basis of resistance is poorly understood. Here, we genomic and functional analyses, we show that ALL-BM mesenchymal stromal cells (MSC) interactions encompass direct cell-cell adhesion via integrin α4β1, uni-directional secretion of soluble factors from MSC to ALL, and bidirectional transfer of cellular vesicles through tunneling nanotubes. We show that this crosstalk drives an epithelial-mesenchymal transition (EMT)-like program in ALL cells adhering to MSC, thereby triggering stem cell-like features including drug resistance and survival advantage. Moreover, single-cell RNA sequencing of cell line and primary ALL samples identified a hybrid cell state derived from B-ALL adhered to MSC, that exhibits both B-ALL and MSC gene signatures that are orchestrated by an EMT-like program. Multiple EMT-activating pathways including WNT/β-catenin, cytokine/STAT3, TGFβ, NOTCH1, and NF-κB-dependent signaling pathways are activated in MSC-adherent ALL cells; however, only the inhibition of either WNT/β-catenin or cytokine/STAT3 pathway could attenuate the survival and drug resistance of ALL cells adherent to MSC. Notably, we identified that this hybrid cell state exhibited elevated WNT/β-catenin-mediated transcriptional activity. Finally, we show that blocking of the interaction of β-catenin and CBP interaction abrogates the survival and drug resistance of ALL cells adhered to MSC, indicating that targeting of an EMT-like program is potential therapeutic regimen in targeting cell extrinsically acquired drug resistance in ALL.

  Study EGAS00001006120

• THAP11 mutations in a patient with a cblX-like phenotype implicates THAP11 in the regulation of cobalamin metabolism and early vertebrate development

  Mutations in HCFC1 are associated with cblX (MIM309541), an X-linked recessive disorder, characterized by defects in cobalamin metabolism and other developmental defects. HCFC1 is a transcriptional co-regulator, which interacts with transcription factors to regulate the expression of a myriad of genes. Notably, HCFC1 regulates cobalamin metabolism via the regulation of MMACHC expression through its interaction with THAP11, a THAP domain-containing transcription factor. In addition, the HCFC1/THAP11 complex potentially regulates genes involved in diverse cellular functions including cell cycle, proliferation, and transcription. Thus, it is likely that mutations in THAP11 may also result in biochemical and other phenotypes similar to those observed in patients with cblX. We report a patient who presented with phenotypic features that overlap cblX, but did not have any mutations in either MMACHC or HCFC1. We sequenced THAP11 by Sanger sequencing and discovered a potentially pathogenic, homozygous variant in THAP11, c.240C>G (p.Phe80Leu). Functional analysis in the developing zebrafish embryo, demonstrate that both THAP11 and HCFC1 regulate the proliferation and differentiation of neural precursors, suggesting important roles in normal brain development. Further, the loss of THAP11 in zebrafish embryos, results in craniofacial abnormalities including the complete loss of Meckel’s cartilage, the ceratohyal, and all of the ceratobranchial cartilages. These data are consistent with our previous work which demonstrated a role for HCFC1 in vertebrate craniofacial development. Furthermore, our high throughput RNA-sequencing analysis reveals several overlapping gene targets of HCFC1 and THAP11 that likely contribute to the phenotypes associated with cblX and related disorders. Thus, both HCFC1 and THAP11 play important roles in the regulation of cobalamin metabolism as well as other pathways involved in early vertebrate development.

  Study EGAS00001002201

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

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

  Study EGAS00001004564

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

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

  Study EGAS00001004567

• Targeting FGFR1 for treatment of soft-tissue sarcoma (H021)

  Purpose: Altered FGFR1 signaling has emerged as therapeutic target in various epithelial malignancies. In contrast, the role of FGFR1 in soft-tissue sarcoma (STS) has not been established. Prompted by the detection and subsequent therapeutic inhibition of amplified FGFR1 in a patient with metastatic leiomyosarcoma, we investigated the oncogenic properties and potential as drug target of FGFR1 in STS.Experimental Design: The frequency of FGFR1 amplification and overexpression, as assessed by fluorescence in situ hybridization, microarray-based comparative genomic hybridization and mRNA expression profiling, SNP array profiling, and RNA sequencing, was determined in three independent patient cohorts. The sensitivity of STS cell lines with or without FGFR1 alterations to genetic and pharmacologic FGFR1 inhibition and the signaling pathways engaged by FGFR1 were investigated using short-term viability assays, long-term colony formation assays, and biochemical analysis.Results: Increased FGFR1 copy number was detected in 74 of 190 (38.9%; Cohort 1), 13 of 79 (16.5%; Cohort 2), and 80 of 254 (31.5%; Cohort 3) patients. FGFR1 overexpression occurred in 16 of 79 (20.2%, Cohort 2) and 39 of 254 (15.4%; Cohort 3) patients. Targeting of FGFR1 by short hairpin RNA-mediated knockdown and different small-molecule inhibitors (PD173074, AZD4547, BGJ398) revealed that the requirement for FGFR1 signaling in STS cells is primarily dictated by FGFR1 expression levels, and identified the MAPK-ERK1/2 axis as critical FGFR1 effector pathway.Conclusion: These data identify FGFR1 as driver gene in multiple STS subtypes and support FGFR1 inhibition, guided by patient selection according to FGFR1 expression and monitoring of MAPK-ERK1/2 signaling, as therapeutic option in this challenging group of diseases.H021

  Study EGAS00001001844

• Profiling the unique protective properties of intracranial arterial endothelial cells

  Cardiovascular disorders, like atherosclerosis and hypertension, are increasingly known to be associated with vascular cognitive impairment (VCI). In particular, intracranial atherosclerosis is one of the main causes of VCI, although plaque development occurs later in time and is structurally different compared to atherosclerosis in extracranial arteries. Recent data suggest that endothelial cells (ECs) that line the intracranial arteries may exert anti-atherosclerotic effects due to yet unidentified pathways. To gain insights into underlying mechanisms, we isolated post-mortem endothelial cells from both the intracranial basilar artery (BA) and the extracranial common carotid artery (CCA) from the same individual (total of 15 individuals) with laser capture microdissection. RNA sequencing revealed a distinct molecular signature of the two endothelial cell populations of which the most prominent ones were validated by means of qPCR. Our data reveal for the first time that intracranial artery ECs exert an immune quiescent phenotype. Secondly, genes known to be involved in the response of ECs to damage (inflammation, differentiation, adhesion, proliferation, permeability and oxidative stress) are differentially expressed in intracranial ECs compared to extracranial ECs. Finally, Desmoplakin (DSP) and Hop Homeobox (HOPX), two genes expressed at a higher level in intracranial ECs, and Sodium Voltage-Gated Channel Beta Subunit 3 (SCN3B), a gene expressed at a lower level in intracranial ECs compared to extracranial ECs, were shown to be responsive to shear stress and/or hypoxia. With our data we present a set of intracranial-specific endothelial genes that may contribute to its protective phenotype, thereby supporting proper perfusion and consequently may preserve cognitive function. Deciphering the molecular regulation of the vascular bed in the brain may lead to the identification of novel potential intervention strategies to halt vascular associated disorders, such as atherosclerosis and vascular cognitive dysfunction.

  Study EGAS00001004479

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

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

  Study EGAS00001004555

• Comprehensive characterization of pre- and post-treatment samples of breast cancer reveal potential mechanisms of chemotherapy resistance

  When locally advanced breast cancer is treated with neoadjuvant chemotherapy, the recurrence risk is significantly higher if no complete pathologic response is achieved. Identification of the underlying resistance mechanisms is essential to select treatments with maximal efficacy and minimal toxicity. Here we employed gene expression profiles derived from 317 HER2-negative treatment-naïve breast cancer biopsies of patients who underwent neoadjuvant chemotherapy, deep whole exome and RNA-sequencing profiles of 22 matched pre- and post-treatment tumors, and treatment outcome data to identify biomarkers of response and resistance mechanisms. Molecular profiling of treatment-naïve breast cancer samples revealed that expression levels of proliferation, immune response and extracellular matrix (ECM) organization combined predict response to chemotherapy. Triple negative patients with high proliferation, high immune response and low ECM expression had a significantly better treatment response and survival benefit (HR 0.29, 95% CI 0.10-0.85; p=0.02), while in ER+ patients the opposite was seen (HR 4.73, 95% CI 1.51-14.8; 0=0.008). The characterization of paired pre-and post-treatment samples revealed that aberrations of known cancer genes were either only present in the pre-treatment sample (CDKN1B) or in the post-treatment sample(TP53, APC, CTNNB1). Proliferation-associated genes were frequently down-regulated in post-treatment ER+ tumors, but not in triple negative tumors. Genes involved in ECM were upregulated in the majority of post-chemotherapy samples. Genomic and transcriptomic differences between pre- and post-chemotherapy samples are common and may reveal potential mechanisms of therapy resistance. Our results show a wide range of distinct, but related mechanisms, with a prominent role for proliferation- and ECM-related genes.

  Study EGAS00001005876

• DNA methylation-based prognostic subtypes of chordoma tumors in tissue

  Chordomas are rare malignant bone cancers of the skull-base and spine. Patient survival is variable and not reliably predicted using clinical factors or molecular features. This study identifies prognostic epigenetic chordoma subtypes that are detected noninvasively using plasma methylomes. Methylation profiles of 68 chordoma surgical samples were obtained between 1996 and 2018 across three international centers along with matched plasma methylomes where available. Consensus clustering identified two stable tissue clusters with a disease-specific survival difference that was independent of clinical factors in a multivariate Cox analysis (HR = 14.2, 95%CI: 2.1–94.8, P = 0.0063). Immune-related pathways with genes hypomethylated at promoters and increased immune cell abundance were observed in the poorperforming “Immune-infiltrated” subtype. Cell-to-cell interaction plus extracellular matrix pathway hypomethylation and higher tumor purity were observed in the better-performing “Cellular” subtype. The findings were validated in additional DNA methylation and RNA sequencing datasets as well as with immunohistochemical staining. Plasma methylomes distinguished chordomas from other clinical differential diagnoses by applying fifty chordoma-versusother binomial generalized linear models in random 20% testing sets (mean AUROC = 0.84, 95%CI: 0.52–1.00). Tissue-based and plasma-based methylation signals were highly correlated in both prognostic clusters. Additionally, leave-one-out models accurately classified all tumors into their correct cluster based on plasma methylation data. Here, we show the first identification of prognostic epigenetic chordoma subtypes and first use of plasma methylome based biomarkers to noninvasively diagnose and subtype chordomas. These results may transform patient management by allowing treatment aggressiveness to be balanced with patient risk according to prognosis.

  Study EGAS00001006406

• Multi-omic analysis of the tumor microenvironment shows clinical correlations in Ph1 study of atezolizumab +/- SoC in MM

  Multiple myeloma (MM) remains incurable, and treatment of relapsed/refractory (R/R) disease is challenging. There is an unmet need for more targeted therapies in this setting; deep cellular and molecular phenotyping of the tumor and microenvironment in MM could help guide such therapies. This phase 1b study (NCT02431208) evaluated the safety and efficacy of the anti-programmed death-ligand 1 monoclonal antibody atezolizumab (Atezo) alone or in combination with the standard of care (SoC) treatments lenalidomide (Len) or pomalidomide (Pom) and/or daratumumab (Dara) in patients with R/R MM. Study endpoints included incidence of adverse events (AEs) and overall response rate (ORR). A novel unsupervised integrative multi-omic analysis was performed using RNA sequencing, mass cytometry immunophenotyping, and proteomic profiling of baseline and on-treatment bone marrow samples from patients receiving Atezo monotherapy or Atezo+Dara. A similarity network fusion (SNF) algorithm was applied to preprocessed data. Eighty-five patients were enrolled. Treatment-emergent deaths occurred in 2 patients; both deaths were considered unrelated to study treatment. ORRs ranged from 11.1% (Atezo+Len cohorts, n=18) to 83.3% (Atezo+Dara+Pom cohort, n=6). High-dimensional multi-omic profiling of the tumor microenvironment and integrative SNF analysis revealed novel correlations between cellular and molecular features of the tumor and immune microenvironment, patient selection criteria, and clinical outcome. Atezo monotherapy and SoC combinations were safe in this patient population and demonstrated some evidence of clinical efficacy. Integrative analysis of high dimensional genomics and immune data identified novel clinical correlations that may inform patient selection criteria and outcome assessment in future immunotherapy studies for myeloma.

  Study EGAS00001007286

• RNA-seq data in hMPC, MSC_Pat and EWIma

  Ewing sarcoma is characterized by pathognomonic translocations fusing most frequently EWSR1 with FLI1 (EF1). In addition, Ewing sarcoma can also display alterations in STAG2, TP53 and CDKN2A (SPC). Starting from Ewing sarcoma derived human mesenchymal stem cells (MSCpat), we recapitulated this translocation and SPC alterations using a CRISPR/cas9 approach and generated a bona fide Ewing sarcoma model (EWIma1) displaying transcriptomic and epigenetic hallmarks of EwS.

  Study EGAS50000001173

• Oropharynx cancers sequencing study (RNA-seq from FFPE)

  This is a study of transcriptomics in HPV-associated oropharynx tumors. The goals of this study were to characterize HPV+ oropharynx tumor subtypes and to examine HPV integration events and their downstream effects. HPV-host fusion transcripts and recurrent HPV integration genes were identified. Various alternative isoforms and splicing events of HPV genes were also studied. Multiple potential biomarkers of prognosis were associated with clinical outcomes.

  Study EGAS50000000893

• PCA Atlas Chromium scRNA-seq and demultiplexing support (FASTQs, BAMs, capture-level VCFs and mapping tables)

  This dataset contains Chromium single-cell RNA-seq data and demultiplexing support for the PCA Atlas study EGAS00001008332. Included data objects: - 16 Chromium scRNA-seq FASTQ runs (EGAR accessions) for captures PCa1–PCa16. - 16 aligned BAM analyses (one Cell Ranger possorted_genome_bam.bam per capture). - 16 capture-level genotype VCF analyses derived from the Axiom UK Biobank array, used as demultiplexing panels. - demux_map.csv (EGAZ00001945520): donor–capture–sample mapping for genotype-based demultiplexing. - chromium.csv (EGAZ00001945537): run-level mapping linking captures, EGA samples, runs (EGAR), experiments (EGAX) and BAM analyses. Together these objects provide a complete view of the Chromium scRNA-seq data and the genotype-based demultiplexing support needed to reproduce donor assignment and downstream analyses.See also https://zenodo.org/records/17372603

  Dataset EGAD00001015795

• single cell RNA-seq of small cell lung cancer circulating tumor cells

  This study investigated the transcriptomic profile of circulating tumor cells from three small cell lung cancer patients. CTCs were enriched from the blood of three SCLC patients using the CTC-iChip followed by magnetic depletion of RBCs. The enriched samples were processed with the 10x Genomics Chromium platform (Chromium GEM-X Single Cell 3' Kit v4) and sequenced on a NextSeq 2000 system.

  Study EGAS50000001401

• Biallelic variants in the non-protein coding minor spliceosome components RNU6ATAC and RNU6ATAC cause syndromic monogenic autoimmune diabetes

  Non-protein coding genes are emerging as critical contributors to the aetiology of rare diseases, providing key insights into human biology and uncovering novel disease mechanisms. We identified 7 individuals from 4 families with early-onset diabetes (diagnosed <5 years) and immune dysregulatory features caused by biallelic variants in RNU6ATAC. RNU6ATAC encodes a small nuclear RNA (snRNA) acting as a catalytic component of the minor spliceosome, a protein-RNA complex mediating splicing of ~700 genes containing U12/minor-type introns. Variant screening of the other 64 minor spliceosome genes in 276 infants with diabetes identified 12 unrelated individuals with biallelic disease-causing variants in RNU4ATAC. Biallelic pathogenic RNU4ATAC variants are known to cause a variable spectrum of clinical features, which until now did not include diabetes. Clinically, 12/19 RNU6ATAC/RNU4ATAC patients had additional immune dysregulatory features, and 50% of patients tested were islet-autoantibody positive, strongly supporting an autoimmune aetiology for their diabetes. RNA-seq in 3 individuals with biallelic RNU6ATAC variants showed a pattern of intron retention in U12 intron-containing genes similar to that seen in RNU4ATAC-individuals (n=3), supporting a shared disease mechanism. Analysis of patients’ transcriptomic, methylation and immune data revealed impaired B cell development and maturation. We conclude that biallelic RNU6ATAC variants cause a novel syndrome of early-onset autoimmune diabetes and immune dysregulation. We further show that infancy-onset diabetes is a previously unrecognised feature of RNU4ATACopathy. Our work highlights the important role of two snRNAs critical to minor spliceosome function in immune system regulation, providing novel insights into the pathogenesis of autoimmune diabetes.

  Study EGAS50000001565

• Exploiting immune cell receptor information to quantify index switching in single cell transcriptome sequencing experiment

  By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms has been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as “index switching”, which causes certain amount of reads to be assigned to an incorrect library. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-unique expression of immune cell receptors to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3500 and HiSeq 4000. By utilizing the unique T-cell receptor (TCR) expression, we could quantify the spread-of-signal from many different wells (n=51 from total three batches) due to index switching. We used TraCer to reconstruct full-length TCR from all samples, and then used Kallisto to quantify TCR gene expression. We found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 4.2% (interquartile range (IQR): 2.0%-8.7%). We did not detect any consistent pattern of certain indices to be more prone for switching than others, suggesting that index switching is a stochastic process. We confirm that index switching is a problem that affects all samples run in multiplexed libraries on Illumina HiSeq 3500 and HiSeq 4000 platforms.

  Study EGAS00001002911

• Single-cell RNA-Seq of human lymphoma reveals malignant B cell diversity and patterns of T cell immune checkpoint co-expression

  Follicular lymphoma (FL) is a generally incurable B-cell malignancy which has the potential to transform into highly aggressive lymphomas. Genomic studies indicate it is often a small subpopulation rather than the dominant population in the FL that gives rise to the more aggressive subtype. To resolve the underlying transcriptional networks of follicular B-cell lymphomas at single molecule and cell resolution, we leveraged droplet-based barcoding technology for highly parallel single cell RNA-Seq. We analyzed the transcriptomes from tens of thousands of cells derived from five primary FL tumors. Simultaneously, we conducted multi-dimensional flow cell sorting to validate our characterizing of cellular lineages and critical expressed proteins. For each tumor, we identified multiple cellular subpopulations, matching known hematopoietic lineages. Comparison of gene expression by matched malignant and normal B cells from the same patient revealed tumor-specific features. Malignant B cells exhibited restricted immunoglobulin light chain expression (either Ig Kappa or Ig Lambda), as well the expected upregulation of the BCL2 gene, but also down-regulation of the FCER2, CD52 and MHC class II genes. By leveraging the single-cell resolution on large numbers of cells per patient, we were able to examine tumor-resident T cells. We identified pairs of immune checkpoint molecules that were co-expressed, providing a potentially useful strategy for selection of patient-tailored combination immunotherapies. In summary, massively parallel measurement of single-cell expression in thousands of tumor cells and tumor-resident lymphocytes can be used to obtain a systems-level view of the tumor microenvironment and identify new avenues for therapeutic development.

  Study phs001378

• Single cell RNAseq of stenotic, inflamed and non-inflamed transmural lesions from patients with Crohn's disease

  This dataset contains a gene-cell matrix derived from single-cell RNA sequencing (scRNA-seq) data of ileal tissue from Crohn's disease (CD) patients and colorectal cancer (CRC) patients. It includes: Crohn's Disease Patients: A trio of transmural lesions (stenotic, inflamed, and non-inflamed) from each patient. Colorectal Cancer Patients: Unaffected ileal tissue used as external non-inflamed control. Cell Level Metadata: The dataset includes relevant cell-level metadata such as cell type annotations used in the study. Experimental Details: Platform: 10x Genomics Chromium Single Cell 3' GEX Sequencing: Illumina NovaSeq Processing: Data processed with Cell Ranger software. Resulting count matrices were merged for downstream analysis, including integration and dimensionality reduction. Dataset Composition: Crohn's Disease Patients: 10 patients with 3 samples each (non-inflamed, inflamed, stenotic), totaling 30 samples. Colorectal Cancer Patients: 5 patients with 1 sample each of unaffected tissue, totaling 5 samples. Data Provided: Merged Raw Count Matrix: The final merged raw count matrix used for downstream analysis. Cell Metadata File: Contains details of sample, tissue, and patient for each cell in the count matrix. Barcodes File: Indicate each cell barcode which also encodes the sample, tissue, and patient details for each cell. CD.S_Inf: Stenotic Corhn's disease inflamed samples CD.S_Sten: Stenotic CD patient stenosis sample CD.S_Prox: Stenotic CD Patient - proximal non-inflamed sample CC.C_Prox: CRC Patient proximal unaffected sample eg: A barcode 'CC.C_1_Prox_AAGTCGTAGACCCTTA' indicates CRC Patient unaffected proximal sampe from CRC Patient no.1 and the nucleic acid sequence indicate a unique cell from this sample. Total Samples: Crohn's Disease (CD) Patients: 30 samples Colorectal Cancer (CRC) Patients: 5 samples Patient_no Sample Sample_type 1 CC.C_1 CC.C_1_Prox CC.C_Prox 2 CD.S_1 CD.S_1_Prox CD.S_Prox 3 CD.S_1 CD.S_1_Infl CD.S_Infl 4 CD.S_1 CD.S_1_Sten CD.S_Sten 5 CC.C_2 CC.C_2_Prox CC.C_Prox 6 CD.S_2 CD.S_2_Prox CD.S_Prox 7 CD.S_2 CD.S_2_Infl CD.S_Infl 8 CD.S_2 CD.S_2_Sten CD.S_Sten 9 CC.C_3 CC.C_3_Prox CC.C_Prox 10 CC.C_4 CC.C_4_Prox CC.C_Prox 11 CD.S_3 CD.S_3_Prox CD.S_Prox 12 CD.S_3 CD.S_3_Infl CD.S_Infl 13 CD.S_3 CD.S_3_Sten CD.S_Sten 14 CD.S_4 CD.S_4_Prox CD.S_Prox 15 CD.S_4 CD.S_4_Infl CD.S_Infl 16 CD.S_4 CD.S_4_Sten CD.S_Sten 17 CC.C_5 CC.C_5_Prox CC.C_Prox 18 CD.S_5 CD.S_5_Prox CD.S_Prox 19 CD.S_5 CD.S_5_Infl CD.S_Infl 20 CD.S_5 CD.S_5_Sten CD.S_Sten 21 CD.S_6 CD.S_6_Prox CD.S_Prox 22 CD.S_6 CD.S_6_Infl CD.S_Infl 23 CD.S_6 CD.S_6_Sten CD.S_Sten 24 CD.S_7 CD.S_7_Prox CD.S_Prox 25 CD.S_7 CD.S_7_Infl CD.S_Infl 26 CD.S_7 CD.S_7_Sten CD.S_Sten 27 CD.S_8 CD.S_8_Prox CD.S_Prox 28 CD.S_8 CD.S_8_Infl CD.S_Infl 29 CD.S_8 CD.S_8_Sten CD.S_Sten 30 CD.S_9 CD.S_9_Prox CD.S_Prox 31 CD.S_9 CD.S_9_Infl CD.S_Infl 32 CD.S_9 CD.S_9_Sten CD.S_Sten 33 CD.S_10 CD.S_10_Prox CD.S_Prox 34 CD.S_10 CD.S_10_Infl CD.S_Infl 35 CD.S_10 CD.S_10_Sten CD.S_Sten

  Dataset EGAD50000000559