495 results for "cancer rna-seq"

in 41.15 milliseconds.

• Breast Cancer Single-Cell RNA-Seq Dataset

  Single-cell RNA-Sequencing of 26 primary breast cancers from Wu et al. (2021) study. Data was generated using the Chromium controller (10X Genomics) and sequenced on the NextSeq 500 platform.

  Dataset EGAD00001007495

• Exploration of coding and non-coding variants in cancer using GenomePaint.

  GenomePaint (https://proteinpaint.stjude.org/genomepaint) is a dynamic visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data, featuring a novel design that captures the inter-relatedness between DNA variations and RNA expression. Regulatory non-coding variants can be inspected and discovered along with coding variants, and their functional impact further explored by examining 3D genome and/or ChIP-seq data generated from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and can display external data such as adult cancer datasets and user-provided custom tracks. We used GenomePaint to analyze multi-omics data from 3,652 pediatric cancers representing 16 histotypes, and demonstrate the visualization features through examples, including two that led to new insights into oncogenic mechanisms in pediatric cancer. The first is the discovery of a new class of pathogenic recurrent variants that cause aberrant splicing, disrupting the RING domain of CREBBP, a driver gene frequently mutated in relapsed pediatric leukemia. The second is the cis-activation of the MYC oncogene in a subset of B-lineage acute lymphoblastic leukemia (B-ALL) via duplication of the NOTCH1-MYC enhancer (N-ME), previously discovered only in T-lineage ALL. The regulatory impact of N-ME enhancer amplification was initially confirmed by allelic imbalance in published gene expression and ChIP-seq data and verified by additional Capture-C and fluorescence in situ hybridization data generated by follow-up experiments. These examples demonstrate the power of GenomePaint in enabling not only data visualization but also integrative genomic analysis that can lead to novel biological insight for follow-up experimental validation.

  Study EGAS00001004669

• Cancer and germline exomes, and cancer RNA-seq consisiting of FASTQ paired-end reads from melanoma, lung and colon cancer samples

  Discovery of patient-specific tumor antigens usually requires in vitro-expanded autologous tumor infiltrating lymphocytes (TILs), among which tumor antigen-specific T cells are however often rare, thus limiting sensitivity. We designed in vitro culture conditions to improve the identification of rare tumor antigen-specific CD8 TILs. This innovative yet accessible pipeline allows highly-sensitive identification of tumor antigens and cognate T cell receptors (TCRs), greatly improving the selection of candidates for personalized cancer vaccines and TCR-based cellular immunotherapies.

  Study EGAS00001005513

• RNA-seq analysis of metastatic prostate cancer solid tumor biopsies

  The prevalence and features of treatment-emergent small-cell neuroendocrine prostate cancer (t-SCNC) are not well characterized in the era of modern androgen receptor (AR)–targeting therapy. We sought to characterize the clinical and genomic features of t-SCNC in a multi-institutional prospective study. Patients with progressive, metastatic castration-resistant prostate cancer (mCRPC) underwent metastatic tumor biopsy and were followed for survival. t-SCNC is present in nearly one fifth of patients with mCRPC and is associated with shortened survival.

  Study EGAS00001006369

• Evolutionary Analysis of Chronic Lymphocytic Leukemia Cells During Relapse After Allogeneic Hematopoietic Stem Cell Transplant

  We investigated the molecular dynamics driving relapse of chronic lymphocytic leukemia (CLL) following allogeneic hematopoietic stem cell transplant (allo-HSCT) through an integrated analysis of whole exome, bulk transcriptome, bulk methylome and single cell transcriptional data from paired leukemic samples representing 19 patients. We found that their clinical kinetics corresponded to distinct genetic and epigenetic evolutionary trajectories and suggest that the selection pressures of the post-transplant, immunologic bottleneck are distinct from those imposed by chemotherapy.In order to discover cancer antigens derived from annotated and unannotated protein-coding regions of the genome, we carried out ribosome profiling (Ribo-seq) and RNA-sequencing on CLL samples and isolated primary B cells from healthy donors. We discovered novel or unannotated open reading frames (nuORFs), their expression levels, as well as somatic mutations within them to predict potential neoantigens.

  Study phs001998

• Microbiome-Host Interactions in Oral Squamous Cell Carcinoma: a Metatranscriptomic Exploratory Study

  Oral cancer (predominantly squamous cell carcinoma; OSCC) is the 14th most prevalent malignancy worldwide, accounting for more than 300,000 new cases and 145,000 deaths annually. Despite advances in cancer treatment modalities, OSCC continues to have poor prognosis with 5-year survival rates less than 50%. Although use of tobacco and alcohol consumption constitute the major risk factors, there is emerging evidence for the role of the microbiome in OSCC. However, studies on the association between the microbiome and OSCC has so far focused on compositional analysis and has been largely limited to the bacterial component of the microbiome. The aim of this study is to characterize the transcriptional activity (function) as well as composition of the multi-kingdom microbiome associated with OSCC and correlate it with the gene expression of the host using ultra-deep metatranscriptome sequencing, which will provide new insight into the role of the microbiome in OSCC and the possible mechanisms by which the tumor’s microbial community may influence on its behavior. This project comprises RNA-Seq data from twenty OSCC tumor/normal tissue pairs and twenty tongue scraping samples from healthy controls. However, due to consent constraints, the human sequences from the control data are filtered out.

  Study phs002678

• Luminal progenitor cell line iHBEC(CD117)

  Basal-like breast cancer originates in luminal progenitors, frequently with an altered PI3K pathway, and focally in close association with genetically altered myoepithelial cells at the site of tumor initiation. The exact trajectory behind this bi-lineage phenomenon remains poorly understood. Here we used a breast cancer relevant transduction protocol including hTERT, shp16, shp53, and PIK3CA(H1047R) to immortalize FACS isolated luminal cells, and we identified a candidate multipotent progenitor. We found that the apparent luminal phenotype of these oncogene transduced progenitors was metastable giving rise to basal-like cells dependent on culture conditions. After culturing the cells for more than 60 passages, cells were subjected to scRNA-seq as well as bulk RNA sequencing of two subpopulations (CD271+ and CD271-).

  Dataset EGAD50000000723

• Single-cell RNA-seq data from metastatic ovarian cancer for quality control study

  Quality control is a crucial preliminary step in any single-cell RNAseq experiment, where hard thresholds are commonly used. In order to develop a methodology for a more precise cell filtering in the early steps of a scRNAseq data analysis, we collected tissue samples before chemotherapy from 4 patients.

  Study EGAS00001005066

• 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

• The Role of E2F4 in Controlling Resistance to Irinotecan (CPT-11) in Human Colon Cancer

  This study reports the results of RNA sequencing (RNA-seq) experiments performed on three isogenic clones of the human HT29 colon cancer cell line, genetically engineered to achieve the selective inactivation of the adenovirus early 2 (E2) gene transcription factor 4 (E2F4) protein, a key functional component of the multimeric transcriptional repression complex known as the dimerization partner (DP), retinoblastoma-like (RB-like), E2F and multi-vulva class B (MuvB) complex (DREAM). In previous studies, we identified E2F4 as a key regulator of colorectal cancer (CRC) resistance to irinotecan (CPT-11), a cytotoxic agent that stands at the backbone of multiple chemotherapy regimens for the treatment of metastatic CRCs. The molecular mechanisms linking E2F4 to irinotecan resistance, however, remain poorly understood. The present study was designed to identify which genes are under the transcriptional control of E2F4, either at baseline or following exposure to SN38 (i.e., the active metabolite of irinotecan, generated in vivo by biotransformation in the liver) in order to clarify the role played by E2F4 in shaping colon cancer resistance to chemotherapy. The cell line used for this study (HT29) is aneuploid and contains three copies of the E2F4 gene, because of a trisomy of chromosome 16, where the E2F4 gene is located (Kawai et al., Genes, Chromosomes & Cancer, 34:1-8, 2002; PMID: 11921276). HT29 cells were infected with lentivirus vectors encoding for CRISPR/Cas9 constructs designed to selectively inactivate the E2F4 gene, using three distinct guide-RNAs (gRNAs), each targeting a different sequence within the proximal exons of the E2F4 gene (exon 1, exon 2). Lentivirus-infected cells were sub-cloned by single-cell fluorescence-activated cell sorting (FACS), based on the differential expression of a fluorescent reporter (EGFP) expressed in tandem with the gRNAs. Isogenic clones (EGFP+) were then screened to identify those with a tri-allelic knock-out, as revealed by both genetic sequencing of the E2F4 gene (revealing the presence of three distinct frameshift mutations) and Western blot (revealing complete loss of E2F4 protein expression). The study includes twelve (n=12) RNA-seq experiments, representative of 3 experimental replicates (i.e., three distinct isogenic clones of the HT29 cell line, each displaying complete loss of E2F4 expression) and 3 negative controls (i.e., three independent preparations of the parent HT29 cell line, infected with a lentivirus vector encoding for a non-targeting gRNA) exposed to 2 distinct in vitro cell culture conditions: a) treatment with SN38 resuspended in dimethyl sulfoxide (DMSO); or b) treatment with DMSO alone (negative control). The in vitro treatment regimen used for SN38 (32 nM, 24 hours) was designed so that the exposure to SN38 experienced in vitro by HT29 cells would mirror the exposure to SN38 experienced in vivo by malignant tissues, as measured, for example, in metastatic CRC patients receiving a standard dose of irinotecan (180 mg/m2) administered by intravenous infusion over 90 minutes (Deyme et al., Cancer Chemotherapy and Pharmacology, 88:247-258, 2021; PMID: 33912999).

  Study phs003560

• UROMOL 2020 - RNA-seq data

  The molecular landscape in non-muscle-invasive bladder cancer (NMIBC) is characterized by large biological heterogeneity with variable clinical outcomes. Here, we performed a large integrative multi-omics analysis of patients diagnosed with NMIBC (n=834). Transcriptomic analysis identifies four classes (1, 2a, 2b and 3) reflecting tumor biology and disease aggressiveness. Both transcriptome-based subtyping and the level of chromosomal instability provide independent prognostic value beyond established prognostic clinicopathological parameters.

  Study EGAS00001004693

• ProstOmics - Spatial Transcriptomics

  Spatial transcriptomics data (ST) from 32 human prostate tissue samples originating from 8 prostate cancer patients (5 patients with post-surgery relapse). The ST data was acquired using the Visium Spatial Gene Expression kit which resulted in over 20 000 spatially defined spots for the 32 tissue samples. The raw transcriptomics data is RNA-seq. The individual samples have information on patient origin and sample type (cancer, cancer-adjacent field-effect normal or normal sample far from cancer). Each ST spot has metadata such as sample origin, histology class (stroma, normal epithelium, cancer of various grading etc), number of cells and estimated cell type fractions. Patient metadata include information of age at surgery, time (months) until reported relapse, total follow-up time, pre-surgery PSA, post-surgery T-stage and metastasis status.

  Dataset EGAD50000000603

• Single MII oocyte mRNA expression from women

  This dataset contains 14 paired-end FASTQ sequences from mRNA-Seq on single human M-II stage oocytes that were collected from ovarian tissue from unstimulated patients undergoing fertility preservation treatments due to cancer diagnoses, which did not influence ovarian function. Cumulus-oocyte-complexes were matured in vitro according Gruhn et al (Science 365: 1466-1469) and short term flash frozen prior to lysis, RNA extraction, full length cDNA preparation and amplification using the Ultra-low-input SMART-Seq2 v4 kit from Takara Clonetech. Further, these cDNA were used to prepare libraries for sequencing according the Nextera XT DNA library preparation kit from Illumina

  Dataset EGAD00001005971

• On-target mutations confer resistance to WRN helicase inhibitors in Microsatellite Unstable Cancer Cells

  This study investigates genetic dependencies and resistance mechanisms to Werner helicase inhibitors (WRNi) in microsatellite-unstable (MSI) cancer models. To define determinants of WRN dependency and acquired resistance, we performed genome-wide CRISPR-Cas9 knockout screens, CRISPR base-editing mutagenesis screens targeting WRN, targeted sequencing of resistance-associated variants, bulk RNA sequencing, and whole-genome sequencing. The datasets deposited in EGA comprise raw sequencing files from CRISPR-Cas9 knockout screens, CRISPR base-editing screens, targeted DNA sequencing, RNA-seq experiments, and whole-genome sequencing of parental and drug-treated cell populations, including models with spontaneously acquired resistance. These data enable investigation of WRN synthetic lethality, genetic modifiers of WRNi response, and genomic alterations emerging under selective pressure.

  Dataset EGAD00001015822

• Intratumoral plasma cells predict outcomes to PD-L1 blockade in non-small cell lung cancer

  Inhibitors of the programmed cell death-1 (PD-1/PD-L1) signaling axis are approved to treat non-small cell lung cancer (NSCLC) patients, based on their significant overall survival (OS) benefit. Using transcriptomic analysis of 891 NSCLC tumors from patients treated with either the PD-L1 inhibitor atezolizumab or chemotherapy from two large randomized clinical trials, we found a significant B cell association with extended OS with PD-L1 blockade, independent of CD8+ T cell signals. We then derived gene signatures corresponding to the dominant B cell subsets present in NSCLC from single-cell RNA-seq data. Importantly, we found increased plasma cell signatures to be predictive of OS in patients treated with atezolizumab, but not chemotherapy. B cells were also associated with the presence of tertiary lymphoid structures and organized lymphoid aggregates. Our results suggest an important contribution of B and plasma cells to PD-L1 blockade efficacy in NSCLC.

  Study EGAS00001005013

• Combined Tumor and Immune Signals From Genomes or Transcriptomes Predict Outcomes of Checkpoint Inhibition in Melanoma

  Cancer immunotherapy with checkpoint blockade (CPB) leads to improved outcomes in melanoma and other tumor types, but a majority of patients do not respond. High tumor mutation burden (TMB) and high levels of tumor-infiltrating T cells have been associated with response to immunotherapy, and integrative models using DNA and RNA have improved predictions of clinical benefit. However, identification of immune and tumor-intrinsic features in models using DNA or RNA alone have not been comprehensively explored in melanoma. In this study, we sequenced DNA and RNA of tumors from 67 melanoma patients receiving CPB, and aggregated this cohort with previously published data, yielding whole exome sequencing (WES) data for 189 patients and bulk RNA sequencing data for 178 patients. Using these datasets, we derived genomic and transcriptomic factors that predict overall survival (OS) and response to immunotherapy. Using whole exome DNA data alone, we calculated T cell burden (TCB) and B cell burden (BCB) based on rearranged TCR/Ig DNA sequences, and found that melanoma patients with high TMB and either high TCB or high BCB survived longer and had higher response rates, as compared to patients with low TMB and either low TCB or low BCB. Next, using bulk RNA-Seq data, differential expression analysis identified 83 genes which are associated with high or low OS. By combining pairs of immune-expressed genes with tumor-expressed genes, we identified three gene pairs which are associated with response and survival (Bonferroni P<0.05). All three gene pair models were validated in an independent cohort (n=180) (Bonferroni P<0.05). The best performing gene pair model included the lymphocyte-expressed MAP4K1 (Gene ID: 11184) and the transcription factor TBX3 (Gene ID: 6926) which is overexpressed in poorly differentiated melanomas. We concluded that RNA-based (MAP4K1 and TBX3) or DNA-based (TCB and TMB) models combining with immune and tumor measures improve predictions of outcome after checkpoint blockade in melanoma.

  Study phs002683

• Molecular_risk_stratification_in_patients_with_T1_colorectal_cancer_WES

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

  Study EGAS00001005734

• National Cancer Institute Cancer Genome Characterization Initiative (CGCI)

  The Office of Cancer Genomics at the National Cancer Institute sponsored a series of studies as part of the Cancer Genome Characterization Initiative (CGCI) to assess novel emerging sequencing technologies in cancer. The CGCI program included comprehensive characterization of the genetic aberrations found in different pediatric and/or adult tumors. CGCI characterized a number of B-cell non-Hodgkin lymphomas (including diffuse large B-cell lymphoma (DLBCL) from patients with and without HIV+ infection, follicular lymphoma (FL), as well as adult and pediatric Burkitt lymphomas), and additional HIV-associated tumors (including DLBCL, lung, and cervical cancers). All data from these projects is available at the NCI Genomic Data Commons. Individual project descriptions are available by disease on the substudy pages (links are included at the top of this page). Brief summaries are as follows: Non-Hodgkin Lymphoma (NHL) - CGCI investigators probed genomic alterations more deeply than had previously been possible by using state-of-the-art RNA sequencing (mRNA-seq) and whole genome shotgun sequencing (WGS) coupled with leading edge bioinformatics, data management and analysis approaches. The project sequenced tumor DNA and/or RNA from 117 NHL tumor samples and 10 cell lines. This included the genomes or exomes of 1 Follicular Lymphoma (FL) and 13 diffuse large B-cell lymphoma (DLBCL) cases, all with matched constitutional DNA sequenced to comparable depths, RNA-sequencing (mRNA-seq) of 92 DLBCL, 12 FL and 8 B-cell NHL cases with other histologies and 10 DLBCL-derived cell lines. The DLBCL cases and cell-lines are from the two major subtypes of DLBCL: germinal center B-cell (GCB) and activated B-cell (ABC). HIV+ Tumor Molecular Characterization Project (HTMCP) - This project was a joint effort of the Office of Cancer Genomics (OCG) and the Office of HIV and AIDS Malignancy (OHAM). Its goals were to characterize HIV-associated cancers (obtained from HIV-infected patients) and compare them to the same types of cancers from patients without HIV infection. Investigators performed 30X genome sequencing of 100 cases of paired tumor and germline DNA, along with transcriptome sequencing in each of 3 types of HIV+ tumors (DLBCL, lung and cervical cancers). These platforms allow discovery of mutations both in coding and non-coding genomic regions, gene expression and genomic alterations (including translocations, insertions and deletions). Comparing tumors of cancer patients both with and without HIV-infection provides insight into the potential function of this virus in certain cancers. Burkitt Lymphoma Genome Sequencing Project (BLGSP) - This project was a collaborative effort between the National Cancer Institute and the Foundation for Burkitt Lymphoma Research to develop a databank of the many alterations found in Burkitt lymphoma (BL), an uncommon type of Non-Hodgkin lymphoma that occurs most often in children and young adults. The goal of the BLGSP was to explore potential genetic changes in patients with BL that could lead to better prevention, detection and treatment of the cancer. The project characterized the alterations of the tumors' genomes (with matched normal as control) and transcriptomes by sequencing the DNA and RNA of each case. Using the data generated, the ultimate goals of the project was to discover the molecular changes that are present in BL patients and then determine how those changes correlate with treatment regimen and outcome. CGCI data is accessible at the NCI's Genomic Data Commons (GDC) via the GDC's Data Portal and from each of the CGCI Publication Pages at the GDC (Please see the "Supplemental Links" section of any CGCI publication's "Publication Information and Associated Data Files page" at the GDC.). Available datasets include raw sequencing data, datasets generated by the original CGCI research teams (fully annotated clinical information and higher-level/analyzed molecular characterization data), as well as higher level data generated by the GDC. Please see the CGCI Use and Publication Guidelines for updated details on the sharing of any CGCI substudy data, including how to cite CGCI.To learn more about the CGCI studies, please visit the CGCI Program website.

  Study phs000235

• Molecular profiling of tissue autopsies and ctDNA

  The dataset contains high-throughput sequencing data derived from a cancer autopsy series of 10 patients. As part of this study, whole-exome sequencing and RNA-seq were performed for spatially distinct tissue biopsies from the patients. In addition, plasma samples from the patients were sequenced using a custom panelt to profile ctDNA. There are 106 files containing whole-exome sequencing data, 107 files containing RNA-seq data, and 9 files containing plasma sequencing data.

  Dataset EGAD00001007040

• Ipilimumab and Decitabine in Treating Patients With Relapsed or Refractory Myelodysplastic Syndrome or Acute Myeloid Leukemia

  The NCI Experimental Therapeutics Clinical Trials Network (ETCTN)/Cancer Therapy Evaluation Program (CTEP) 10026 study (NCT02890329) tested combined ipilimumab and decitabine treatment in two cohorts: Arm A consisted of patients with relapsed acute myeloid leukemia (AML) after allogeneic hematopoietic stem cell transplantation. Arm B consisted of patients with advanced myelodysplastic syndrome (MDS) or AML that did not previously undergo stem cell transplantation. During the study, peripheral blood and bone marrow samples were obtained from study participants at study entry (screening), following the first lead-in cycle of decitabine monotherapy, during combined treatment with decitabine and ipilimumab, and at the end of treatment.This study tested the combination of the CTLA-4 blocking antibody ipilimumab and the DNA methyltransferase inhibitor decitabine in advanced myeloid disease in the transplant-naive and post-transplant context. The samples were used for bulk and single cell RNA sequencing (scRNA-seq) as well as whole exome sequencing (WES).

  Study phs003015

• Spatial Transcriptomics Uncovers Tumor Microenvironment-Based Subtypes in Invasive Lobular Carcinoma

  Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer (BC); however, little is known about its tumor microenvironment (TME). Here, we investigated the spatial architecture of 43 hormone-receptor-positive, HER2-negative (HR+/HER2-) primary ILC samples, using spatial transcriptomics and detailed morphological annotation. By integrating morphology, spatial, and single-cell data, we mapped TME-related cell types, exploring their proximity to tumor cells and their enrichment in specific compartments. Clustering of spatial transcriptomics spots revealed significant inter- and intra-tumor heterogeneity, identifying pathways and cell types linked with disease outcomes. Morphological and gene expression data were finally integrated with a multimodal approach, leading to the identification of 4 new ILC subtypes based on TME heterogeneity. By using derived, specific gene expression signatures, these subtypes were identified and validated in external bulk RNA-seq/microarray datasets (SCAN-B and METABRIC), where differences in prognosis emerged.

  Study EGAS50000001001

• 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

• Single cell RNAseq and TCRseq data from tumor and blood samples from 4 patients with muscle invasive bladder cancer

  For this dataset we performed single cell RNAseq paired with single cell TCR-seq on tumor and blood samples from 4 patients. This dataset contains 4 tumor samples as well as 4 blood samples. Each sample is made up of 2 sets of paired fastq files. The first pair contains reads corresponding to RNA transcripts (_Transcripts in file name), while the second pair contain reads corresponding to TCRs (_VDJ in file name). Sequenced on the Illumina NovaSeq6000 platform in a paired-end run using an SP flow cell (v1.5, 300 cycles).

  Dataset EGAD50000001381

• 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

• Dissociation of solid tumour tissues with cold active protease for single-cell RNA-seq minimizes conserved collagenase-associated stress responses

  Background Single-cell RNA sequencing (scRNAseq) is a powerful tool for studying complex biological systems, such as tumour heterogeneity and tissue microenvironments. However, the sources of technical and biological variation in primary solid tumour tissues and patient-derived mouse xenografts for scRNAseq, are not well understood. Here, we used low temperature (6°C) protease and collagenase (37°C) to identify the transcriptional signatures associated with tissue dissociation across a diverse scRNAseq dataset comprising 128,481 cells from patient cancer tissues, patient-derived breast cancer xenografts and cancer cell lines.Results We observe substantial variation in standard quality control (QC) metrics of cell viability across conditions and tissues. From FACS sorted populations gated for cell viability, we identify a sub-population of dead cells that would pass standard data filtering practices, and quantify the extent to which their transcriptomes differ from live cells. We identify a further subpopulation of transcriptomically “dying” cells that exhibit up-regulation of MHC class I transcripts, in contrast with live and fully dead cells. From the contrast between tissue protease dissociation at 37°C or 6°C, we observe that collagenase digestion results in a stress response. We derive a core gene set of 512 heat shock and stress response genes, including FOS and JUN, induced by collagenase (37°C), which are minimized by dissociation with a cold active protease (6°C). While induction of these genes was highly conserved across all cell types, cell type-specific responses to collagenase digestion were observed in patient tissues. We observe that the yield of cancer and non-cancer cell types varies between tissues and dissociation methods.Conclusions The method and conditions of tumour dissociation influence cell yield and transcriptome state and are both tissue and cell type dependent. Interpretation of stress pathway expression differences in cancer single cell studies, including components of surface immune recognition such as MHC class I, may be especially confounded. We define a core set of 512 genes that can assist with identification of such effects in dissociated scRNA-seq experiments.

  Study EGAS00001003753

• Genomic Profiling of Pediatric B-cell Acute Lymphoblastic Leukemia

  Pediatric leukemia is the most common type of childhood cancer and, despite tremendous progress in treatment, remains the second leading cause of cancer deaths in children. B-cell acute lymphoblastic leukemia (B-ALL) is the most common type of pediatric leukemia. Recent studies using transcriptional profiling of large cohorts of pediatric B-ALL have shown that patients can be grouped into transcriptional subtypes that are often defined by specific driver mutations. Interestingly, many of these mutational subtypes are defined by gene fusions or rearrangements that drive altered expression of transcription factors (ETV6-RUNX1, PAX5, DUX4, TCF3-PBX1, ZNF384) or in chromatin regulatory machinery (KMT2A), suggesting that transcriptional dysregulation is a major driver of B-ALL. While the transcriptional distinctions are clear, the relationship between patients in terms of their chromatin regulatory landscape remains unclear. To address this, we profiled gene expression, open chromatin, and 3D genome folding in patient derived B-ALL samples. We identified that at the level of 3D genome organization, many chromatin interactions and looping events are shared across mutational subtypes, suggesting that shared transcriptional programs, in addition to group specific driver mutations, possibly regulate the chromatin landscape in B-ALL. Using RNA-seq and ATAC-seq data, we were able to identify "latent" signatures of expression related to non-mutant transcription factors (TFs) that also contribute to expression heterogeneity across TF subtypes. These are related to patient specific patterns of 3D genome folding, and in some cases, are predictive of disease outcomes. This study helps shed light on the chromatin landscape in B-ALL and furthers our understanding of the interactions between diverse mutant TFs that drive B-ALL and non-mutant TFs that help shape the B-ALL chromatin landscape.

  Study phs003226

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

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

  Study phs003413

• Gene Expression Study of Individuals with Sex Chromosome Aneuploidies

  Males and females show dramatic differences in their vulnerability to the same diseases. For example, compared to men, lupus is six times more prevalent, thyroid cancer is three times more prevalent, and unipolar depression is twice as prevalent in women. Diseases with a strong male bias include autism (5:1), dilated cardiomyopathy (3:1), and ankylosing spondylitis (5:1). Historically, such differences have been attributed solely to extrinsic factors such as circulating sex hormones or environmental influences. We hypothesized that intrinsic factors - genetic differences between XX and XY cells - have unappreciated biological consequences throughout the body and contribute to sex differences in disease incidence and severity. This hypothesis stems from our long-term effort to sequence the sex chromosomes of diverse mammalian species, which has identified a set of homologous genes on the X and Y chromosomes that are dosage-sensitive, expressed throughout the body, and encode regulators of chromatin modification, transcription, translation, and protein stability. These X- and Y-encoded genes differ in sequence and expression pattern, which likely manifests in genome-wide differences in gene regulation between XX and XY cells and influences all aspects of human biology, including sex differences in disease susceptibility. These hard-wired molecular sex differences have been largely overlooked and understudied, representing a significant gap in our knowledge of human biology.The gene expression study of individuals with sex chromosome aneuploidies takes advantage of natural human variation in sex chromosome number, i.e. sex chromosome aneuploidy, to investigate alterations in genome-wide gene expression that correlate with changes in X- and Y-chromosome dosage. We analyzed samples from 114 individuals with a variety of sex chromosome aneuploidies, including 45,X; 47,XXY; 47,XYY; 47,XXX; 48,XXYY; and 49,XXXXY. We generated lymphoblastoid cell lines (LCLs) from blood samples and, in some cases, fibroblast cultures from skin biopsies. We supplemented our collection with previously-derived cell lines. To evaluate gene expression, we performed deep profiling of the transcriptome (RNA-seq) from these LCLs and fibroblasts. We performed parallel analyses on samples collected from 62 control 46,XX and 46,XY individuals, 6 individuals with trisomy 21, and 14 individuals with structural variations of the X and Y chromosomes. In addition, we performed CRISPRi knockdowns on 3 of the 46,XX and 3 of the 46,XY fibroblast samples for the homologous transcription factors ZFX and ZFY, encoded on the X and Y chromosomes, respectively.In the April 2024 update, we added RNA-seq datasets derived from isolated CD4+ T cells and monocytes from 76 and 72 adults, respectively, with the following sex chromosome constitutions: 45,X; 46,XX; 46,XY; 47,XXX; 47,XXY; 47,XYY. These individuals are largely a subset of the same cohort described above. In addition, we performed RNA-seq on in-vitro stimulated CD4+ T cells with the following sex chromosome constitutions: 45,X; 46,XX; 46,XY; 47,XXY.In the August 2024 update, we added RNA-seq datasets generated from the following: 1) LCLs derived from individuals with AZFa deletions of the Y chromosome, 2) DDX3X and DDX3Y knockdown (via CRISPRi) in XY fibroblasts, and 3) 5-ethyl uridine (5-EU) treatment in XY and XYYYY LCLs.

  Study phs002481

• CPTAC: Proteogenomic Studies of Ovarian Tumor Responses to Agents Targeting the DNA Damage Response

  There is a critical unmet need for predictors of chemo-refractory High-grade serous ovarian cancer (HGSOC). Despite over 3 decades of research on platinum responses in cancer, no predictive biomarker has been translated into clinical use. Predictors of refractory disease could spare these patients the unnecessary toxicity of a platinum-based regimen and provide a means to triage these patients to clinical trials to identify effective therapies for refractory disease. This study has approached the development of a predictor in a novel way enabling detection of a multi-analyte panel biomarker capable of detecting a subset of refractory HGSOCs with high specificity and implicating potential treatment vulnerabilities for refractory disease. We reported the first proteogenomic analysis of platinum refractoriness in HGSOC and made the novel dataset available as a resource to the community, including a ProTrack portal: http://ptrc.cptac-data-view.org/. Detailed proteogenomic profiles were generated for this study from a retrospective cohort of 242 treatment-naive HGSOCs with enriched representation of "exceptional non-responders" (i.e., patients with chemo-refractory disease). Whole genome sequence (WGS) and RNA seq FASTQ files are available through dbGaP. Proteomic data are available through the NCI Proteomic Data Commons (PDC): https://pdc.cancer.gov/pdc/.

  Study phs003152

• Molecular determinants of response to PD-L1 blockade across tumor types

  Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer (NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity was observed between indications, in particular within the PD-L1+ tumors. mUC and NSCLC were molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor CDKN2A was identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis revealed molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.

  Study EGAS00001004343

• Resolution of tumour cell populations enhances specificity of treatment options for precision cancer medicine

  This dataset comprises FASTQ files for single-cell RNA-seq and single-cell ATAC-seq data. These data were generated using 10x Genomics' multiome kit. In this cohort, there are 25 samples from 17 human cancer patients.

  Dataset EGAD00001015782

• Chromatin accessibility maps of chronic lymphocytic leukemia identify subtype-specific epigenome signatures and transcription regulatory networks

  Chronic lymphocytic leukemia (CLL) is characterized by substantial clinical heterogeneity, despite relatively few genetic alterations. To provide a basis for studying epigenome deregulation in CLL, we established genome-wide chromatin accessibility maps for 88 CLL samples from 55 patients using the ATAC-seq assay, and we also performed ChIPmentation and RNA-seq profiling for ten representative samples. Based on the resulting dataset, we devised and applied a bioinformatic method that links chromatin profiles to clinical annotations. Our analysis identified sample-specific variation on top of a shared core of CLL regulatory regions. IGHV mutation status – which distinguishes the two major subtypes of CLL – was accurately predicted by the chromatin profiles, and gene regulatory networks inferred for IGHV-mutated vs. IGHV-unmutated samples identified characteristic differences between these two disease subtypes. In summary, we discovered widespread heterogeneity in the chromatin landscape of CLL, established a community resource for studying epigenome deregulation in leukemia, and demonstrated the feasibility of chromatin accessibility mapping in cancer cohorts and clinical research.

  Dataset EGAD00001002110

• The Effect of the Menstrual Cycle on DNA Expression in the Normal Human Breast Epithelium

  The Susan G. Komen for the Cure® Tissue Bank at the IU Simon Cancer Center [www.komentissuebank.iu.edu] (KTB) was established expressly for the prospective collection of normal, healthy breast tissue from volunteer donors. Blood is also obtained from donors at the time of donation and is processed for serum, plasma and peripheral blood leukocyte DNA. Specimens are annotated and data includes age, race, ethnicity, personal health history, family cancer history, medication usage at the time of donation, and breast cancer risk factors. Premenopausal donors to the KTB were identified by a query of the Bank's database. Hematoxylin and eosin stained sections of the formalin-fixed paraffin-embedded tissue of the identified donors were reviewed and tissue was graded on the basis of the abundance of epithelium within the section. Only tissue containing abundant epithelium was considered for this study. Based on dates, the specimens of nine women in the follicular phase of the menstrual cycle and five in the luteal phase were chosen. Six donors using hormonal contraception at the time of donation were also included. Whole blood obtained from 19 of the 20 donors at the time of tissue donation was available and it was processed for serum. Estradiol, estriol, luteinizing hormone and progesterone concentrations were determined by the Indiana University Health Pathology Laboratory using a Beckman Unicel DxI 800 Immunoassay System. The phase of the menstrual cycle was verified by serum progesterone concentration. The epithelium of these 20 specimens was microdissected from multiple 8 micron frozen tissue sections. Total RNA extracted from the tissue was subsequently depleted of rRNA via locked nucleic acid probes. This enabled profiling of both poly-A and non-poly-A RNA species. Barcoded cDNA libraries from the 20 normal breast epithelia were prepared and sequenced on an Applied Biosystems (AB) SOLiD3 or SOLiD 4 platform. RNA-seq reads for each sample were then mapped to the human genome (hg19) using the LifeScope software version 2.5.1 (Life Technologies, Foster City, CA) and BAM (Binary Alignment/Map) files generated. Read counts for each gene were derived from the output BAM files using the RefSeq database (UCSC Genome Brower) as the gene model.

  Study phs000644

• Functional and genomic heterogeneity of long-term self-renewing compartment as the origin of treatment resistance in pancreatic tumors

  Adaptive drug-resistance mechanisms allow human tumors to evade treatment through selection and expansion of treatment-resistant clones. Here, studying clonal evolution of tumor cells derived from human pancreatic tumors, we demonstrate that in vitro cultures and in vivo tumors are maintained by a common set of tumorigenic cells that can be used to establish clonal replica tumors (CRTs), large cohorts of animals bearing human tumors with identical clonal composition. Using CRTs to conduct quantitative assessments of adaptive responses to therapeutics, we uncovered a multitude of functionally heterogeneous subpopulations of cells with differential degrees of drug sensitivity. High-throughput isolation and deep characterization of unique clonal lineages showed genetic and transcriptomic diversity underlying functionally diverse subpopulations. Molecular annotation of gemcitabine-naïve clonal lineages with distinct responses to treatment in the context of CRTs generated signatures that can predict the response to chemotherapy, representing a potential biomarker to stratify patients with pancreatic cancer. RNA-Seq data on these clones are available on EBI ENA, under accession number: PRJEB27735.

  Study EGAS00001003442

• Comprehensive Genomic Characterization of Translocation Renal Cell Carcinoma

  Xp11 translocation renal cell carcinoma (tRCC) is a rare, female-predominant cancer driven by a fusion between the TFE3 gene on chromosome Xp11.2 and a partner gene on either chrX or an autosome. It remains unknown what types of rearrangements underlie TFE3 fusions, whether fusions can arise from both the active (chrXa) and inactive X (chrXi) chromosomes, and whether TFE3 fusions from chrXi translocations account for the female predominance of tRCC. To address these questions, we performed haplotype-specific analyses of chrX rearrangements in tRCC whole genomes. We show that TFE3 fusions universally arise as reciprocal translocations and that oncogenic TFE3 fusions can arise from chrXi:autosomal translocations. Female-specific chrXi:autosomal translocations result in a 2:1 female-to-male ratio of TFE3 fusions involving autosomal partner genes and account for the female predominance of tRCC. Our results highlight how X chromosome genetics uniquely constrains somatic chrX alterations and underlies cancer sex differences.This version of the dbGaP study provides a comprehensive set of whole genome data and tables with clinical information collected for the 16 patients (15 pateints as tRCC and 1 patient as ccRCC) in the study cohort. The data also contains 12 RNA-seq samples used in the study and analyses. We anticipate that this data will considerably advance our molecular and genetic understanding of oncogenic translocations in TFE3-tRCC involving inactive chromosome X and its role in cancer sex bias.

  Study phs003008

• Childhood Cancer Data Initiative (CCDI): Pediatric In Vivo Testing Program - Neuroblastoma

  The goal of NCI's Pediatric In Vivo Testing Consortium (PIVOT) is to advance the development of effective treatments for pediatric cancers through preclinical in vivo testing of novel therapeutic agents. To design impactful experiments that can be rapidly translated to the clinic, PIVOT investigators require comprehensive genomic characterization of the patient-derived xenograft (PDX) models used across the consortium. This robust experimental design enhances collaboration with industry partners by providing proof-of-concept data for drugs in development. As part of this effort, neuroblastoma PDX models are being molecularly characterized to identify promising therapeutic candidates. The molecular profiling includes whole exome sequencing (WES), RNA-seq, MethylEPIC array, CytoSNP array, and DNA fingerprinting for quality control. Additionally, demographic and clinical data (e.g., diagnosis, disease site, disease status) are shared. These characterized PDX models, matched with patient tumors, guide model selection for preclinical drug testing.

  Study phs003163

• Single cell RNAseq of PBMC from bladder cancer patients

  Frozen PBMC samples containing at least 1 million cells were thawed for 1 minute at 37C and washed twice with RPMI complete media (10% FBS with glutamate and Pen/Strep). All of the samples had >80% viable cells. Sample processing for single-cell RNA-seq was done using Chromium Single Cell 3’ Library and Gel bead kit v2 (PN-120237) following manufacturer’s user guide (CG00052, 10x Genomics, Pleasanton, CA). The total cell density was used to impute the volume of single cell suspension needed in the reverse transcription (RT) master mix, aiming to achieve ~ 6,000 cells per sample. cDNAs and libraries were prepared following manufacturer’s user guide (10x Genomics). cDNA amplification and indexed libraries were prepared using 12 and 14 cycles of PCR, respectively. Libraries were profiled, quantified, and sequenced as 5’ single-cell gene expression libraries.

  Study EGAS00001004008

• Single-cell RNA-seq of immune cells sorted from human melanoma tumors

  Tumor immune cell compositions play a major role in response to immunotherapy but the heterogeneity and dynamics of immune infiltrates in human cancer lesions remain poorly characterized. Here we identify conserved intratumoral CD4 and CD8 T cell behaviors in scRNA-seq data from 25 melanoma patients. We discover a large population of CD8 T cells showing continuous progression from an early effector "transitional" into a dysfunctional T cell state. CD8 T cells that express a complete cytotoxic gene set are rare, and TCR sharing data suggest their independence from the transitional and dysfunctional cell states. Notably, we demonstrate that dysfunctional T cells are the major intratumoral proliferating immune cell compartment and that the intensity of the dysfunctional signature is associated with tumor-reactivity. Our data demonstrate that CD8 T cells previously defined as exhausted, are in fact a highly proliferating, clonal and dynamically differentiating cell population within the human tumor microenvironment.

  Study EGAS00001003363

• Leiomyosarcoma Cancer Transcriptome Sequencing

  RNA-Seq data (fastq raw read sequences) for French ICGC leiomyosarcoma cancer sequencing project, 78 samples representing 67 donors. Sequencing was performed on Illumina HiSeq. The libraries were then sequenced with a 2 x 75bp paired-end protocol to a minimum mean reads of 50 million.

  Dataset EGAD00001003192

• Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high risk disease

  Genomic events associated with poor outcome in chronic myeloid leukemia (CML) are poorly understood. We performed whole exome sequencing, copy number variation and/or RNA-Seq for 65 patients to discover mutations at diagnosis and blast crisis (BC). Forty-six chronic phase patients with the extremes of outcome were studied at diagnosis. Cancer gene variants were detected in 15/27 patients (56%) with subsequent BC or poor outcome and in 3/19 optimal responders (16%), P=.007. Frequently mutated genes at diagnosis were ASXL1, IKZF1 and RUNX1. The methyltransferase SETD1B was a novel recurrently mutated gene. A novel class of variant associated with the Philadelphia translocation was detected at diagnosis in 11/46 patients (24%) comprising fusions and/or rearrangement of genes on the translocated chromosomes, with evidence of fragmentation, inversion and imperfect sequence reassembly. These were more frequent at diagnosis in patients with poor outcome: 9/27 (33%) versus 2/19 optimal responders (11%), P=.07. Thirty-nine patients were tested at BC and all had cancer gene variants, including ABL1 kinase domain mutations in 58%. However, ABL1 mutations co-occurred with other mutated cancer genes in 89% of cases, and these predated ABL1 mutations in 62% of evaluable patients. Gene fusions not associated with the Philadelphia translocation occurred in 42% of patients at BC and commonly involved fusion partners that were known cancer genes (78%). Genomic analysis revealed numerous relevant variants at diagnosis in patients with poor outcome and all patients at BC. Future refined biomarker testing of specific variants will likely provide prognostic information to facilitate a risk-adapted therapeutic approach.

  Study EGAS00001003071

• A Functional Network of Gastric-Cancer-Associated Splicing Events Controlled by Dysregulated Splicing Factors

  Except for a few genes known to have oncogenic spliceforms in gastric cancer, the full scope of aberrant splicing in gastric oncogenesis remains unclear. In this study, we elucidated the alternative splicing events associated with gastric carcinogenesis and identified upstream regulators governing these events. We performed RNA-seq on 19 matched tumor/normal pairs and a panel of gastric cancer cell lines, and identified 361 tumor associated (TA) alternative splicing events (ASEs), which included known oncogenic ASEs in genes such as INSR, FGFR2, CD44 and KRAS. We further identified 8 splicing factors dysregulated in gastric tumors, the expression of which are correlated to the splice ratio of the TA ASEs. We thereby constructed a dysregulated splicing network in gastric cancer, consisting of alternative splicing changes correlated to their potential regulators. This network featured three potential regulatory modules centered around the splicing factors ESRP2, MBNL1 and PTPB1. Knockdown of each splicing factor led to the expected changes in splicing of approximately half of the ASEs assayed. Thus, dysregulation of the splicing factors indeed was the cause of many of the TA ASEs. Pathway enrichment analysis showed that genes affected by tumor associated splicing were pre-dominantly involved in cytoskeletal organization. Concordant with this, knockdown of the splicing factors regulating the TA ASEs led to the expected changes in trans-well migration activity. Thus, we have established for the first time a comprehensive splicing network that likely drives gastric tumor progression, based on which novel therapeutic targets can be identified for further evaluation.

  Study EGAS00001002256

• Multi-Modal Single-Cell, Spatial, and Genomic Analyses of Human Non-Small Cell Lung Cancer Brain Metastases

  We performed single nuclei RNA-sequencing (snRNA-seq) with matched T cell receptor sequencing (TCR-seq), and pool matched low pass whole genome sequencing (WGS) of 12 treatment-naive non-small cell lung cancer (NSCLC) primary tumors (PTs) and 31 treatment-naive NSCLC brain metastases (BMs). In total, we recovered 277,206 cell transcriptomes in 43 samples. Single-nuclei suspension from frozen NSCLC metastases were generated using methods published in Slyper et al. (Nature Medicine, 2020, PMID:32405060) with modifications. The tissue was first cut in 20 µm sections using a cryostat and stored on dry ice until the next step. Then the tissue was put on wet ice and washed once with 5 ml PBS to remove OCT. After this step nuclei were extracted by vigorous pipetting in salt-tris buffer with Tween-20 (TST) followed by a high-volume wash in salt-tris buffer (ST). Nuclei were counted and 10,000-14,000 and were loaded using a Chromium controller (10X genomics) and Chromium Single Cell 5' Reagents (V1.1 and V2). Complementary DNA libraries were generated according to manufacturer's instructions, with the modification of one additional cycle of amplification. Matched single nuclei T cell receptor (snTCR) sequencing libraries were prepared from amplified cDNA libraries using either Chromium Single Cell V(D)J Enrichment Kit for human T cells Reagent Kit, v1.1, or single cell human TCR amplification kit. Final snTCR sequencing libraries were prepared using Single Index Kit T Set A for V1.1 libraries or Library construction kit and Dual Index Kit TT set A for V2 libraries. Sequencing libraries were quantified using Tapestation D5000 reagents and a 2200 TapeStation and pooled for sequencing. Pool matched low pass whole-genome sequencing (lp-WGS) was performed in addition to snRNA-seq to improve the identification of the malignant and non-malignant cell compartments. Nuclei from snRNA-seq sample preparations were collected by centrifugation (500g, 5 min), snap frozen after removing excess ST buffer, and stored at -200C until further processing. If insufficient nuclei (0C for 9 minutes. Adapters were not diluted, and the optional PCR (5 cycles) and cleanup were performed. After library cleanup with AMPure beads (Beckman Coulter, Brea, CA), libraries were quantified with Tapestation D5000 HS tapes (Agilent) and stored at -200C until sequencing.Using multiplexed immunofluorescence in an independent cohort of treatment-naïve pairs of primary tumors and brain metastases from the same patients with NSCLC, we validated genomic and tumor-microenvironmental findings and identified a cancer cell population characterized by neural features that strongly enriched in brain metastases. This comprehensive analysis provides insights into human NSCLC brain metastasis biology and serves as an important resource for additional discovery.

  Study phs003865

• RNA-seq consisting of FASTQ paired-end reads from cancer samples

  Efforts to precisely identify tumor human leukocyte antigen presented peptides (HLAp) capable of mediating T cell based tumor rejection still face important challenges. Recent reports suggest that non-canonical cancer HLAp could be immunogenic but their identification requires highly sensitive and accurate mass-spectrometry (MS)-based proteogenomics approaches. Here, we present a MS-based analytical pipeline that can precisely characterize the non-canonical HLAp repertoire, incorporating whole exome sequencing, bulk and single cell transcriptomics, ribosome profiling, and a combination of two MS/MS search tools. This approach results in the accurate identification of hundreds of shared and tumor-specific non-canonical HLAp. Albeit often at low levels and in distinct subpopulations of cells, numerous non-canonical HLAp are shared across tumors. This analytical platform holds great promise for the discovery of novel cancer antigens for cancer immunotherapy.

  Study EGAS00001003724

• FOXM1 is a biomarker of resistance to PI3Kα inhibition in ER+ breast cancer that is detectable using metabolic imaging (RNA-seq data)

  PIK3CA, encoding the PI3Kα isoform, is the most frequently mutated oncogene in estrogen receptor-positive breast cancer. Isoform-selective PI3K inhibitors are used clinically but intrinsic and acquired resistance limits their utility. Improved selection of patients that will benefit from these drugs requires novel predictive biomarkers. We show here that persistent FOXM1 expression following drug treatment is a biomarker of resistance to PI3Kα inhibition in ER+ breast cancer. Since the transcription factor drives expression of lactate dehydrogenase (LDH) but not hexokinase 2 (HK-II), changes in FOXM1 expression can be detected using magnetic resonance imaging of LDH-catalyzed hyperpolarized 13C label exchange between pyruvate and lactate but not by PET measurements of HK-II-mediated trapping of the glucose analog 2-deoxy-2-[18F]fluorodeoxyglucose.

  Study EGAS00001004452

• Onco-exaptation of an Endogenous Retroviral LTR Drives IRF5 Expression in Hodgkin Lymphoma

  The transcription factor interferon regulatory factor 5 (IRF5) is upregulated in Hodgkin lymphoma (HL) and a key regulator of the aberrant transcriptome characteristic of this disease. Here we show that IRF5 upregulation in HL is driven by transcriptional activation of a normally dormant endogenous retroviral LOR1a long terminal repeat (LTR) upstream of IRF5. Specifically, through screening of RNA-seq libraries, we detected LTR-IRF5 chimeric transcripts in multiple HL cell lines but not in normal B-cell controls. In HL, the LTR was in an open and hypomethylated epigenetic state and we further show the LTR is the site of transcriptional initiation. Among HL cell lines, usage of the LTR promoter strongly correlates with overall levels of IRF5 mRNA and protein, indicating that LTR transcriptional awakening is a major contributor to IRF5 upregulation in HL. Taken together, oncogenic IRF5 overexpression in HL is the result of a specific LTR transcriptional activation. We propose that such LTR derepression is a distinct mechanism of oncogene activation (“onco-exaptation”) and that such a mechanism warrants further investigation in molecular and cancer research.

  Study EGAS00001001205

• WTS and WES of Renal Cell Carcinoma Tumors From a Phase III Anti-Angiogenic Adjuvant Therapy Trial

  Archival de-identified formalin-fixed paraffin-embedded RCC tumor tissue blocks from nephrectomy or tumor biopsy were processed as per below and the same sections were used for both DNA and RNA extractions. For WES (ACE version 3; Illumina NovaSeq), samples were profiled using Personalis ACE Cancer Exome (Personalis, Inc, Menlo Park, CA) Whole-transcriptome profiles were generated by RNA-seq (Accuracy and Content Enhanced (ACE) version 3; Illumina NovaSeq) using Personalis ACE Cancer Transcriptome (Personalis, Inc, Menlo Park, CA ) Of the 615 patients in the intent-to-treat population in S-TRAC trial, 193 individual specimens were available for molecular profiling, of which 171 (27.8%) (sunitinib, n = 91; placebo, n = 80) returned results for the WES analysis, and 133 (21.6%) (sunitinib, n = 72; placebo, n = 61) returned results for the GES analysis. Of the 138 WTS samples with data, replicates for two patients were summarized by median expression, and three samples were excluded from the final analysis due to low counts.

  Dataset EGAD00001009294

• 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

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

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

  Study EGAS50000001641

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

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

  Study EGAS00001002687

• Bladder cancer subtyping study across 4 atezo clinical trials

  Checkpoint inhibitors targeting the PD-1/PD-L1 interaction have revolutionized cancer therapy across many indications. In urothelial carcinoma (UC), several PD-1/PD-L1 inhibitors have been approved for treatment in localized and metastatic disease. While these agents provide durable clinical benefit in a number of patients, the majority do not respond. A better understanding of the molecular mechanisms underlying response and resistance is needed to tailor existing treatments and develop new therapeutics in UC. We profiled tumors from 2803 patients from four late-stage randomized clinical trials evaluating the PD-L1 inhibitor atezolizumab in metastatic UC and muscle-invasive bladder cancer by RNA-seq, a targeted DNA panel, immunohistochemistry, and digital pathology. A machine learning algorithm identified four UC transcriptional subtypes, representing luminal desert, stromal, immune, and basal tumors. Overall survival benefit from atezolizumab over standard-of-care was observed in immune and basal tumors (45% of patients). Digital pathology showed enrichment of human interpretable features in specific molecular subtypes, and a self-supervised approach could classify molecular subtypes from H&E slides with high accuracy, suggesting that machine-based image interpretation could accelerate tumor molecular profiling. Somatic alterations showed limited capacity to identify responders. Finally, mechanisms of benefit with PD-L1 blockade appeared different in immune and basal tumors, with immune tumors benefiting from B/plasma cell enrichment, while basal tumors benefited from granulocyte enrichment. This study represents the largest integration of UC molecular and curated clinical data in randomized clinical trials to date, paving the way for prospective clinical studies tailoring treatment regimens to specific tumor molecular subtypes in UC and other indications.

  Study EGAS50000000497

• The Haplotype-Resolved Genome and Epigenome of the Aneuploid HeLa Cancer Cell Line

  The genomic architecture of HeLa remains largely unexplored beyond its karyotype, partly because like many cancers, its extensive aneuploidy renders such analyses challenging. We carried out haplotype-resolved whole-genome sequencing of the HeLa CCL-2 strain, examined point- and indel-mutation variations, mapped copy-number variations and loss of heterozygosity regions, and phased variants across full chromosome arms. We also investigated variation and copy-number profiles for HeLa S3 and eight additional strains. We find that HeLa is relatively stable in terms of point variation, with few new mutations accumulating after early passaging. Haplotype resolution facilitated reconstruction of an amplified, highly rearranged region of chromosome 8q24.21 at which integration of the human papilloma virus type 18 (HPV-18) genome occurred and that is likely to be the event that initiated tumorigenesis. We combined these maps with RNA-seq and ENCODE Project data sets to phase the HeLa epigenome. This revealed strong, haplotype specific activation of the proto-oncogene MYC by the integrated HPV-18 genome approximately 500 kilobases upstream, and enabled global analyses of the relationship between gene dosage and expression. These data provide an extensively phased, high-quality reference genome for past and future experiments relying on HeLa, and demonstrate the value of haplotype resolution for characterizing cancer genomes and epigenomes. Reprinted from Adey et. al. 2013, Nature, with permission from Nature Publishing Group.

  Study phs000642

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

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

  Study EGAS00001005948

• Ither2 RNA-Seq dataset - Implementation of pediatric precision oncology into clinical practice: The individualized Therapies for Children with cancer program “iTHER”

  iTHER is a prospective national precision oncology program aiming to define tumor molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory pediatric tumors in order to identify relevant aberrations to inform treatment.

  Dataset EGAD00001010179

• 'KOREAN' never-smoker female adenocarcinoma RNA-seq

  Multi-layered genomic studies such as the TCGA project have greatly advanced our understanding on molecular pathogenesis of lung cancer. For Asian nonsmoker patients, however, majority of studies have been limited to mutational analyses, emphasizing ethnic differences in driver mutations (e.g. more frequent EGFR mutations and ALK fusions). In essence, a comprehensive multi-layered characterization that could lead to molecular etiology and patient stratification scheme to be translated into clinical applications for this patient group is still missing. Here we report molecular profiling of tumor and matched normal tissues from 114 non-small-cell lung adenocarcinoma patients using whole exome sequencing, transcriptome sequencing, and array comparative genomic hybridization (CGH).

  Study EGAS00001003789

• Clonal Evolution and Heterogeneity of Osimertinib Acquired Resistance Mechanisms in EGFR Mutant Lung Cancer

  Current understanding of acquired resistance to osimertinib in EGFR-mutant lung cancer is largely based on the use of targeted sequencing platforms and cell-free DNA. In this study we examined the clonal evolution and heterogeneity of acquired resistance mechanisms to osimertinib through whole exome and RNA-sequencing of prospectively collected pre-and post-osimertinib resistant tumors and metastases at rapid autopsy. We enrolled 34 patients with histologically confirmed EGFR mutant advanced lung adenocarcinoma in a prospective Phase II clinical trial to evaluate osimertinib treatment and the use of local ablative therapy (LAT) for oligoprogressive disease (NCT02759835) in EGFR mutant NSCLC from April 2016 until data cutoff in September 2019. Patients were not selected based on oligometastatic status, rather 32/34 patients had multiple organ involvement, including lung, liver, adrenal, brain and bone. Biopsies or surgeries were performed pre-osimertinib treatment and at first progression on osimertinib. A subset of patients also underwent biopsies at second progression on osimertinib, and following end of life in-patient hospice and rapid autopsy upon death. 63% (n=15/24) of first-line osimertinib and 50% (n=5/10) of second-line osimertinib treated patients had RECIST defined partial response. Overall, 21 patients had RECIST defined first progression on osimertinib. Two patients on this trial discontinued treatment and were lost to follow up. Of the 21 patients who had RECIST defined progression, 14 patients underwent LAT (8 first-line and 6 second-line). Twelve patients had paired pre- and post-osimertinib resistant tumors and three had post-osimertinib resistant tumors available for analysis. Four patients (LAT002, LAT006, LAT014 and LAT021) consented for a rapid autopsy protocol and underwent rapid autopsy after deatH. To elucidate acquired osimertinib resistance mechanisms, we performed whole-exome (WES), targeted sequencing and RNA-seq on 66 tumors and matched germline DNA from tumor tissue obtained from a total of 15 patients.

  Study phs002001

• Intratumoral Ploidy Dynamics in HCC

  Hepatocellular carcinoma (HCC) ranks as the fourth leading cause of cancer-related deaths. Effective therapies are limited, highlighting the critical need for stratification biomarkers to improve prognosis and predict treatment response. Polyploidization or hyperploidization, an early event observed in many solid cancers, may offer such potential. Here, we developed a new high throughput in situ quantitative approach to assess HCC ploidy profiles. This novel method highlights the presence of both Diploid (DP-) and Hyperploid (HP-) HCCs. RNA-seq profiling classified hyperploid HCCs as belonging to the proliferation subclass, predicting unfavourable clinical outcomes, and revealing a gene expression signature indicative of stemness. Spatially reconstructing nuclear ploidy also revealed DP/HP-HCCs with hyperploid nuclei forming discrete clusters within a predominantly diploid background. Remarkably, these emerging hyperploid foci shared a portion of their transcriptomic programme with fully developed hyperploid HCCs and displayed a signature of genomic instability, including upregulation of multiple DNA repair pathways in response to replication-associated DNA damage. Importantly, spatial transcriptomics of DP/HP HCC demonstrated that hyperploid tumor hepatocytes arise through whole-genome doubling of diploid tumor hepatocytes. Collectively, our findings underscore the prognostic value of nuclear ploidy profiling in HCC and suggest that hyperploid HCCs evolve clonally from intratumor hyperploid foci during tumor progression.

  Study EGAS00001008084

• Adult-type Granulosa Cell Tumour of the Ovary

  A recurrent mutation in FOXL2 (c.402C>G; p.C134W) is present in over 95% of adult-type granulosa cell tumours (AGCTs). In contrast, various loss-of-function mutations in FOXL2 lead to the development of blepharophimosis, ptosis and epicanthus inversus syndrome (BPES). BPES is characterized by an eyelid malformation often accompanied with primary ovarian insufficiency but not granulosa cell tumours. Two recent studies suggest that FOXL2 C402G is a gain or change-of-function mutation with altered DNA-binding specificity. Another study proposes that FOXL2 C402G is selectively targeted for degradation inducing somatic haploinsufficiency suggesting its role as a tumour suppressor. The latter study relies on data indicative of an FOXL2 allelic imbalance in AGCTs. Here we present RNA-seq data as genetic evidence that no real allelic imbalance is observed at the transcriptomic level in AGCTs. Additionally, there is no loss of protein expression in tumours harboring the mutated allele. These data and other features of this mutation as compared to other oncogenes and tumour suppressor genes argue strongly against FOXL2 being a tumour suppressor in this context. Given the likelihood that FOXL2 C402G is oncogenic, targeting the variant protein or its downstream consequences is the most viable path forward to identifying an effective treatment for this cancer.

  Study EGAS00001005414

• Deep single-cell RNA sequencing data for 11138 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive CRC patients. The DATA ACCESS AGREEMENT is provided at https://github.com/zhangyybio/single-T-cell-data-access. Applicants can request access to the data by directly downloading it or by sending an email to cancerpku@pku.edu.cn. The process that is used to approve an application includes verifying the institution, participants and research purposes of the application. In general this process will take about two weeks. In principal, any scientific research program complying with the laws and bioethic regulation policies of China will be approved.

  T cells are central players in cancer immunotherapy1, yet some of their fundamental properties such as development and migration within tumours remain elusive. The enormous T cell receptor (TCR) repertoire, required for recognising foreign and self-antigens2,3, could serve as lineage tags to track these T cells in tumours4. Here, we obtained transcriptomes of 11,138 single T cells from 12 colorectal cancer (CRC) patients and developed STARTRAC (Single T-cell Analysis by Rna-seq and Tcr TRACking) indices to quantitatively analyse dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. While both CD8+ effector and “exhausted” T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, the majority of tumour-infiltrating Tregs showed clonal exclusivity, whereas certain Treg clones were developmentally linked to multiple TH clones. Notably, we identified two IFNG+ TH1-like clusters in tumours, the GZMK+ TEM and CXCL13+ TH1-like clusters, which were associated with distinct IFN-γ-regulating transcription factors, EOMES/RUNX3 and BHLHE40, respectively. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in tumours of microsatellite-instable (MSI) patients, which might explain their favourable response rates to immune-checkpoint blockade. Furthermore, we found IGFLR1 to be highly expressed in both CXCL13+BHLHE40+ TH1-like and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provides a powerful avenue to comprehensively dissect the T cell properties in CRC, which could shed new insights into the dynamic relationships of T cells in other cancers

  Study EGAS00001002791

• Upcycling and merging data to challenge the dogma and identify new therapeutic targets for Glioblastoma

  Glioblastoma is the most common brain tumour. Characterised by a poor prognosis and its recurrence after multimodal treatments, the search for preventable risk factors has been mainly inconclusive up to date. Recently, the data merge from datasetsdeposited at the EGA allowed Aaron Diaz’s team to discover that the glioblastoma cells shift toward a mesenchymal phenotype when the tumour is recurring. Challenges in glioblastoma research As for virtually all cancer types, many efforts have been made to unveil the molecular features responsible of the disease, and several cellular pathways have indeed been identified as being frequently mutated in glioblastoma. Nonetheless, targeted therapies based on identified genes have so far failed to improve outcome, thus survival mostly relies on a standard treatment unchanged since 2005. This is a frustrating situation for both the scientific and medical communities, and above all for the patients, still facing a dreadful path. Some researchers hypothesised that this may be due to the inability to efficiently target cancer stem cells, the originators of the other cell types, thus inducing cancer relapse. In 2019 Charles P. Couturier and colleagues sequenced RNA from single cells of freshly excised glioblastomas of 16 patients, and demonstrated that glioblastoma cells replicate normal brain cell development with a conserved neural cancer cell hierarchy centered around glial progenitor-like cells. In this way, they helped identify the possible target cells to improve efficacy and durability of treatment. Data upcycling at the EGA: the glioblastoma case study Single cell RNA sequencing is generated with a laborious and expensive protocol. The quality of the starting material is crucial (cells sample freshly extracted from the patients) and often several attempts are needed before producing reliable quality sequencing results. Collecting big numbers of patients is also challenging. The sequencing data produced by Kevin Petrecca’s group in Montreal, Canada and deposited at the EGA (EGAS00001004422) was recently upcycled by Lin Wang in San Francisco, California and pooled with their freshly produced ones, and then deposited at the EGA as Dataset EGAS00001004909 The data merge allowed Aaron Diaz’s team to discover that the glioblastoma cells shift toward a mesenchymal phenotype when the tumour is recurring. They profiled 86 primary-recurrent patient-matched paired glioblastoma samples with single-nucleus RNA, among other techniques. These very comprehensive results lead the team to challenge the findings from several other cancer fields where chemotherapy standard chemo-radiation therapy selection pressure at the level of genomic alterations; this is indeed not the case for glioblastoma, where the pressure results in phenotypic transition between cellular states. Several technical controls were made to ensure that the merge of the data was not introducing a bias in the results, and an Inter-table analysis demonstrated nearly equal contribution to overall variance from each of the studies included, indicating that the findings were not due to inter-laboratory technical effect. A novel principal-component analysis showed that the largest contribution to variation in primary glioblastoma neoplastic cells was an axis between MES (mesenchymal) and proneural expression programs. In summary, treating glioblastoma often makes a MES, as commented by Lucy Stead about this remarkable work. Check Aaron Diaz’s team Nature Cancer paper, with strong technical tools, state of the art analysis and data from different sources converging to the same outcomes that made possible a significant step toward a better handling of a frightful cancer. References Couturier, C.P., Ayyadhury, S., Le, P.U. et al. Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy. Nat Commun 11, 3406 (2020). Stead, L.F. Treating glioblastoma often makes a MES. Nat Cancer 3, 1446–1448 (2022). Wang, L., Jung, J., Babikir, H. et al. A single-cell atlas of glioblastoma evolution under therapy reveals cell-intrinsic and cell-extrinsic therapeutic targets. Nat Cancer 3, 1534–1552 (2022). Related links: Kevin Petrecca’s group dataset deposited at the EGA: EGAS00001004422 Lin Wang dataset deposited at the EGA: EGAS00001004909

  Blog upcycling-and-merging-data-to-challenge-glioblastoma

• Childhood Cancer Data Initiative (CCDI): Molecular Characterization across Pediatric Brain Tumors and Other Solid and Hematologic Malignancies for Research, Diagnostic, and Precision Medicine

  This study contains tumor and germline WGS, RNA-Seq, Clinical Panel Sequencing, and other omics and molecular data for patients with pediatric brain tumors and other solid and hematologic malignancies. This data has been collected from multiple organizations including the multi-institute Children's Brain Tumor Network (CBTN), the Pacific Pediatric Neuro-Oncology Consortium (PNOC), and the Children's Hospital of Philadelphia Division for Genomic Diagnostics (CHOP DGD). In addition to being available pre-publication and without embargo via dbGaP and the NCI Cancer Research Data Commons, this data is made available via Cavatica, PedcBioPortal and the Kids First Data Resource Portal. Collectively, this availability strives to accelerate the research process and enable members of the scientific and patient communities to partner and make discoveries faster than ever before.In 2019, the CBTN launched the Pediatric Brain Tumor Atlas, which now comprises the largest collection of childhood brain tumor data in the world. This resource is helping accelerate not only brain tumor research, but is also empowering discovery for other rare childhood conditions. Brain tumors are the most common form of cancer in children aged 0-19 in the United States, and are the largest cause of cancer-related deaths. The estimated number of new cases in 2019 is nearly 3,800, and thus brain tumors are a rare disease. Despite their relative rarity, the years of potential life lost due to brain tumors in 2009 was estimated at 47,631 years for children and adolescents aged 0-19 in the United States; this is a disproportionate amount of life lost compared to adult cancers and represents an unrecognized societal threat. There is an urgent need to improve therapies for these children. This sequencing cohort defines the largest, clinically annotated pediatric brain tumor cohort study to date and seeks to define the intersection of germline and somatic underpinnings of pediatric brain tumors across a shared developmental context of cancer and structural birth defects.

  Study phs002517

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

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

  Study EGAS00001004289

• Circulating Genomic Determinants of Treatment Failure in Hodgkin Lymphoma

  In this study, we show that the plasma representation of mutations exceeds the bulk tumor representation in most cases, making classic Hodgkin Lymphoma (cHL) particularly amenable to noninvasive profiling. Leveraging single-cell transcriptional profiles of cHL tumors, we demonstrate hazard ratios (HRs) circulating tumor DNA (ctDNA) shedding to be shaped by DNASE1L3, whose increased tumor microenvironment-derived expression drives high ctDNA concentrations. Using this insight, we comprehensively profile 366 patients, revealing two distinct cHL genomic subtypes with characteristic clinical and prognostic correlates, as well as distinct transcriptional and immunological profiles. Furthermore, we identify a novel class of truncating IL4R-mutations that are dependent on IL13 signaling and therapeutically targetable with IL4R blocking antibodies. Finally, we demonstrate the clinical value of pre- and on-treatment ctDNA levels for longitudinally refining cHL risk prediction, and for detection of radiographically occult minimal residual disease. Collectively, these results support the utility of noninvasive strategies for genotyping and dynamic monitoring of cHL as well as capturing molecularly distinct subtypes with diagnostic, prognostic, and therapeutic potential.Methods utilized include a) targeted genotyping of baseline plasma samples, b) whole exome sequencing (WES) of baseline plasma samples, c) epigenetic expression inference from cell-free DNA-sequencing (EPIC-Seq) to infer gene expression from cell-free DNA (cfDNA) fragmentation profiles, d) ctDNA tracking in on-treatment samples, e) RNA-Sequencing, and f) Phased variant Enrichment and Detection Sequencing (PhasED-Seq).All patients were treated at cancer centers across Europe and North America between 2011 and 2020. The following 3 cancer centers, Stanford (CA, USA), Bellinzona [Confoederatio Helvetica (CH) or Switzerland], and Leuven (Belgium), were included. The 6 clinical trial collectives included in this study are the Phase III trial AHOD1331 trial (NCT02166463) and the Phase II cHOD17 study (NCT03755804) for enrolled pediatric patients, the Phase III AHL2011 (NCT01358747), the Phase II BREACH trial (NCT02292979), a pilot study evaluating concurrent Pembrolizumab with AVD (Doxorubicin, also known as Adriamycin, Vinblastine, and Dacarbazine; NCT03331341), and a Phase II study evaluating PVAB (Prednisone, Vinblastine, Doxorubicin and Bendamustine; NCT02414568) in older patients (61 years and older). In addition, lymph node specimens from a biorepository at St. Jude Children's Research Hospital (TN, USA) were included.

  Study phs003435

• The Role of CDX2 in Controlling ABCB1 Expression and Chemosensitivity in Human Colon Cancer

  This study reports the results of RNA sequencing (RNA-seq) experiments performed on two human colon cancer cell lines (HT-29, DLD-1) that were genetically engineered to achieve either the constitutive expression (in HT-29 cells) or the selective inactivation (in DLD-1 cells) of the transcription factor CDX2, a master regulator of embryonic gut development and a key regulator of multi-lineage differentiation in adult colon epithelial cells. Our previous studies had identified lack of CDX2 expression as a prognostic biomarker of poor prognosis and a predictive biomarker of benefit from adjuvant chemotherapy in Stage-II and Stage-III colon cancer patients (Dalerba et al., NEJM, 374:211-222, 2016; PMID: 26789870). The molecular basis of these associations, however, remained poorly understood. The present study was designed to identify which genes are under the transcriptional control of CDX2, in order to clarify their role in shaping colon cancer resistance to chemotherapy. Of the two lines included in this study, the first (HT-29) is CDX2-negative at baseline, and was infected with a lentivirus vector encoding for the cDNA of the human CDX2 gene under the control of a constitutive promoter (CMV). The second line (DLD-1) is CDX2-positive at baseline, and was infected with a series of lentivirus vectors encoding for CRISPR/Cas9 constructs designed to selectively inactivate the CDX2 gene. In both cases, infected cells were purified by fluorescence-activated cell sorting (FACS), based on the differential expression of a green fluorescent reporter (EGFP) that was also encoded by the lentivirus vectors, and transcribed in tandem with either the CDX2 cDNA (HT-29) or CRISPR/Cas9 gRNAs targeting CDX2 (DLD-1). The experiments were designed to include 3 experimental replicates per condition, and to compare the transcriptional profile of genetically engineered cells with that of autologous controls (i.e., HT-29 cells infected with an "empty" version of the same lentivirus vector used to over-express CDX2, DLD-1 cells engineered with a "non-targeting" CRISPR/Cas9 gRNA). In the specific case of DLD-1 cells engineered with CRIPSR-Cas9 constructs targeting CDX2, infected cells were sorted based on the differential expression of EGFP and then plated as single-cells, in order to sub-culture three isogenic clones, all of which were genetically sequenced to confirm bi-allelic knock-out of the CDX2 gene, and confirmed to lack CDX2 protein expression by Western blot. The results of our experiments demonstrate that, in human colon epithelial cells, CDX2 controls the expression of the ABCB1 efflux pump, by positively regulating the transcription of the corresponding gene.

  Study phs002903

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

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

  Study phs001687

• How to use the EGA search box

  Tips on how to search Tip 1: Use simple and clear keywords No matter what you are searching: a study, a dataset or a committee; the search will look recursively in all the fields of all the entities in the database. Thus, there is no need to specify which type of data you are looking for. Try to use single and specific words and avoid long and complex expressions. Try it out: cancer rna-seq instead of cancer studies using RNA-seq. Tip 2: Narrow your results using search operators Search operators give you more control over the results. They manage the way keywords are used to restrict the search. By default, the search engine applies an AND operator when a space is inserted between two words. Include OR (capitalised) between the keywords in your search in order to show pages containing either term or both terms. Try it out: Search all the studies based on transcriptomic data, including those ones that use RNA sequencing to obtain this data. Tip 3: Do not worry about details The search engine automatically checks for the most common spelling mistakes of given words and suggests the use of the correct version. It also ignores capital letters and most punctuation so none of these mistakes will hamper the accuracy of your search. Besides, when using a combination of words, the search engine will suggest similar combinations with a higher number of results. Tip 4: Use words that are likely to appear within the text Unfortunately, our search engine is not able to interpret words yet. Therefore, words which do not appear in the body text will not return any result despite them denoting a group of present terms. For instance, if we are looking for studies focused on neurological diseases, using neurological as a keyword would not be the best option. Instead, you should use the name of each neurological disease such as bipolar disorder, schizophrenia or depression. In the same fashion, incomplete words do not return any result. Please include an asterisk (*) before and/or after the incomplete word in order to allow for partial matching. Try it out: canc* instead of canc(without wildcard). Also you can try out *oma, or even *geneti* Tip 5: Narrow down your search by filtering the results You can filter the results by the following options: general, study, and dataset. Under these categories, you will find more sub-filters to refine your search. Make sure that you click on Apply in order to filter the results. You can clear the filtering anytime. Advanced Searches Additional example queries: Search for DUO:0000005 and word lung example query Or write: DUO:0000005 AND lung Search for an exome sequencing study in dbGaP example query Or write: study_type: "exome sequencing" AND repository:dbGaP Search for an exome sequencing study in EGA example query Or write: study_type: "exome sequencing" AND repository:ega Search for a dataset that is at the EGA beacon and with type "Whole genome sequencing" and with number of samples between 100 and 2000 and with technology "Illumina HiSeq 2000" example query Or write: is_in_beacon:true AND dataset_types:"Whole genome sequencing" AND samples:[100 TO 2000] AND dataset_technologies:"Illumina HiSeq 2000" Search for PMIDS example query Or write: pmids:32025007

  Documentation discovery/metadata/search-box

• Proteogenomic Characterization Unveils Biomarkers Associated With Chemoresistance in Muscle Invasive Bladder Cancer

  Cisplatin-based neoadjuvant chemotherapy (NAC) combined with radical cystectomy has been used as a standard of care for patients with clinically non-metastatic muscle invasive urothelial bladder cancer, but 40-50% patients will benefit from the treatment. Molecular subtypes of muscle invasive bladder cancer (MIBC) provide a possible biomarker roadmap that may help identify patients who may or may not respond to NAC. Recently, by integrating Mass Spectrometry-based proteomics, proteogenomic studies from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) have demonstrated the potential of this approach in prioritizing candidate driver genes and facilitating understanding of the chemotherapy resistance mechanisms. To demonstrate the proteogenomic workflow for investigating the mechanisms underlying chemotherapy resistance in MIBC, we conducted proteogenomic characterization of 46 pre- and 14 post- treatment tumors, including whole exome sequencing (WES), mRNA sequencing, Tandem Mass Tag (TMT)-proteomics, or TMT-phosphoproteomic. These samples were selected from an initial cohort of 143 MIBC tissues and filtered by tumor purity, clinical history, pathological response, and proteomic quality control metrics. The final cohort comprised samples from 53 patients, including 45 with both genomic and proteomic data and 8 with genomic data only.Our multi-omics-based clustering reveals subtypes associated with chemotherapy sensitivity. Protein abundance of a short isoform of ATAD1 and RAF family proteins were identified as biomarkers of chemosensitivity which highlights the importance of proteomics data. Integration of proteomic and phosphoproteomic data reveals Wnt signaling via GSK3B-S9 phosphorylation and the JAK/STAT pathway as potential targets to overcome chemoresistance. Correlations between PD-L1 and TROP-2/NECTIN-4 indicate an additive benefit of combination therapy targeting these proteins. Our data analysis clearly demonstrates that proteogenomic plays an important role in uncovering potential resistance mechanisms that operate only at protein and Post-Translational Modification (PTM) levels. Proteomics data can be accessed through PRIDE with ID PXD060290. RNA-seq and whole-exome sequencing (WES) FASTQ files are archived in the dbGaP database.

  Study phs004049

• Comprehensive molecular and clinicopathological profiling of lung adenocarcinoma in Japanese never or light smokers

  Studies have yet to characterize the differences in molecular profiles of lung adenocarcinoma (LUAD) among divergent ethnic groups. Herein, we conducted comprehensive molecular profiling of LUAD in never or light smokers from Asia in order to discover novel targetable mutations and prognostic biomarkers of this distinct disease entity. Methods: We analyzed 996 cases of Asian LUAD, and performed whole-exome sequencing and/or RNA-seq in 125 cases of Asian LUAD negative for the driver oncogenes defined by conventional laboratory testing. We also investigated the clinical and pathological characteristics among the 996 cases. Results: Driver oncogenes were identified in 88 cases (70.4%) with specific hotspot mutations differing from those in The Cancer Genome Atlas (TCGA) study. Clustering based on mRNA expression profiles, but not genetic mutational ones, could predict patient prognosis. The risk score generated by the expression of a three-gene set was a strong prognostic marker for overall survival and progression-free survival in our cohort, and was further validated using TCGA cohort. Among the 996 cases, each driver alteration is distributed across all histological subtypes. Adenocarcinoma in situ was identified to harbor driver mutations, suggesting that these alterations are early events in the pathogenesis of LUAD. This study indicates the value of applying gene expression profiling for predicting the prognosis after a surgical operation, and that the identification of actionable mutations is important for optimizing targeted drugs in Japanese LUAD.

  Study JGAS000215

• scRNA-seq of HGSC tumor and ascites

  Immune regulatory metabolites are key features of the tumor microenvironment (TME), yet with a few notable exceptions, their identities remain largely unknown. We uncovered the immune regulatory metabolic states and metabolomes of sorted tumor and stromal, CD4+, and CD8+ cells from the tumor and ascites of patients with high-grade serous ovarian cancer (HGSC) using high-dimensional flow cytometry and metabolomics supplemented with single cell RNA sequencing. Flow cytometry revealed that tumor cells show a consistently greater uptake of glucose than T cells, but similar mitochondrial activity. Cells within the ascites and tumor had pervasive metabolite differences, with a striking enrichment in 1-methylnicotinamide (MNA) in T cells infiltrating the tumor compared to ascites. Despite the elevated levels of MNA in T cells, the expression of nicotinamide N-methyltransferase, the gene encoding the enzyme that catalyses the transfer of a methyl group from S-adenosylmethionine to nicotinamide, was restricted to fibroblasts and tumor cells. Treatment of T cells with MNA resulted in an increase in T cell-mediated secretion of the tumor promoting cytokine tumor necrosis factor alpha. Thus, the TME-derived metabolite MNA contributes to an alternative and non-cell autonomous mechanism of immune modulation of T cells in HGSC. Collectively, uncovering the tumor-T cell metabolome may reveal metabolic vulnerabilities that can be exploited using T cell-based immunotherapies to treat human cancer.

  Study EGAS00001004829

• 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

• RNA-seq of Glioblastoma stem cells

  Chromatin accessibility discriminates stem from mature cell populations, enabling the identification of primitive stem-like cells in primary tumors, such as Glioblastoma (GBM) where self-renewing cells driving cancer progression and recurrence are prime targets for therapeutic intervention. We show, using single-cell chromatin accessibility, that primary GBMs harbor a heterogeneous self-renewing population whose diversity is captured in patient-derived glioblastoma stem cells (GSCs). In depth characterization of chromatin accessibility in GSCs identifies three GSC states: Reactive, Constructive, and Invasive, each governed by uniquely essential transcription factors and present within GBMs in varying proportions. Orthotopic xenografts reveal that GSC states associate with survival, and identify an invasive GSC signature predictive of low patient survival. Our chromatin-driven characterization of GSC states improves prognostic precision and identifies dependencies to guide combination therapies.

  Study EGAS00001003070

• Somatic Mutation Profile by Next Generation Sequencing in HER2+ Breast Cancer

  The HER pathway is the driving force behind 30% of human breast cancers. It is important to understand how targeted therapies block different cellular pathways, and mechanisms of escape from this blockage. Therapies directed at HER2 establish a successful treatment paradigm, but de novo and acquired resistance exist. The HER signaling system is a complex network with four receptors and eleven ligands, a phosphorylation signaling cascade, and many transcription factors, all complicated by both positive and negative feedback circuits. Analysis of genomes, exomes and transcriptomes by next generation sequencing is aimed at uncovering the genetic factors responsible for patient responses to HER2-directed therapies. We are sequencing HER2-overexpressing cancers, in order to identify potential somatic changes that may better select patients who will benefit from therapy, to determine new targets that may overcome resistance, and to improve outcomes with known current HER2-targeted therapies. Whole exome capture sequencing will determine somatic mutation profiles in HER2-overexpressing tumors, to comprehensively characterize the somatic alterations, with the goal of identifying those patients most likely to respond, as well as discovering new targets that may overcome resistance to HER2-directed therapy. We augment the whole exome data with RNA-seq data to determine expression levels of somatic mutations we discover.

  Study phs000770

• An oncogenic enhancer-rearrangement causes concomitant deregulation of EVI1 and GATA2 in leukemia. Targeted resequencing of chromosomal regions centered on 3q21 and 3q26 in conjunction with RNA-Seq from Acute Myeloid Leukemia patients.

  Chromosomal rearrangements without gene-fusions have been implicated in leukemogenesis by causing deregulation of proto-oncogenes via relocation of cryptic regulatory DNA elements. AML with inv(3)/t(3;3) is associated with aberrant expression of the stem-cell regulator EVI1. Applying functional genomics and genome-engineering, we demonstrate that both 3q-rearrangements reposition a distal GATA2 enhancer to ectopically activate EVI1 and simultaneously confer GATA2 functional haploinsufficiency, previously identified as the cause of sporadic familial AML/MDS and MonoMac/Emberger syndromes. Genomic excision of the ectopic enhancer restored EVI1 silencing and led to growth inhibition and differentiation of AML cells, which could be replicated by pharmacologic BET-inhibition. Our data show that structural rearrangements involving chromosomal repositioning of enhancers can cause deregulation of two unrelated distal genes, with cancer as the outcome.

  Study EGAS00001000669

• Integrated single-cell profiling dissects cell-state-specific enhancer landscapes of human tumor infiltrating T cells.

  Despite extensive studies on the chromatin landscape of exhausted T cells, the transcriptional wiring underlying functional and dysfunctional states of human tumor infiltrating lymphocytes (TILs) is incompletely understood. Here, we identify tissue-specific and general gene-regulatory landscapes in the wide breadth of CD8+ TIL functional states covering four cancer entities using single-cell chromatin profiling. We map enhancer-promoter interactions in human TILs by integrating single-cell chromatin accessibility with single-cell RNA-seq data from tumor entity-matching samples and prioritize key elements by super enhancer analysis. Our analysis reveals a human core chromatin trajectory to TIL dysfunction and identifies key enhancers, transcriptional regulators, and deregulated target genes involved in this process. Finally, we validate enhancer regulation at loci encoding PD1, TCF1, and TIM3 by targeting non-coding regulatory elements with potent CRISPR activators and repressors. In summary, our study advances the understanding of molecular regulation of TIL (dys-)function and provides a framework for modulating immunotherapeutic relevant TIL genes via their enhancers.

  Dataset EGAD00001008662

• Genetic Aberrations and Subclonal Structure Impact Chronic Lymphocytic Leukemia

  Large-scale whole-exome sequencing (WES) of primary tumor samples enables the unbiased discovery of recurrent putative driver events and patterns of clonal evolution. We report the identification of 44 recurrently mutated genes and 11 recurrent CNVs through the WES of 538 chronic lymphocytic leukemia (CLL) and matched germline DNAs. These include previously unrecognized cancer drivers (e.g., RPS15, IKZF3), and collectively identify nuclear export, MYC activity and MAPK signaling as central pathways affected by somatic mutation in CLL. A clonality analysis of this large dataset further enabled the reconstruction of temporal relationships between these driver events. Several drivers were associated with shorter progression-free survival (PFS) in 280 samples that were collected prior to uniform treatment with front line chemo-immunotherapy, with mature follow up of greater than 10 years. Direct comparison between matched pretreatment and relapse CLL from 59 samples demonstrated marked clonal evolution occurring in more than 95% of these patients. Distinct patterns of clonal evolution in relationship to specific gene alteration were observed, suggesting a hierarchy of fitness amongst mutations. Thus, large WES datasets of clinically informative samples enable the discovery of novel driver genes as well as the network of relationships between the drivers and their impact on disease relapse and clinical outcome.Additionally, we performed RNA-seq for 268 CLL samples (including 26 follow-up samples) and used them to identify expression subtypes of CLL. RRBS for 30 of these samples was also generated. In an integrative analysis of genetic, transcriptomic, and epigenetic data, incorporating known and newly identified subtypes of CLL, we built new models to improve patient prognostication.

  Study phs000922

• Super enhancers define regulatory subtypes and cell identity in neuroblastoma - RNA-seq

  Half of the children diagnosed with neuroblastoma have high-risk disease, disproportionately contributing to overall childhood cancer-related deaths. In addition to recurrent gene mutations, there is increasing evidence supporting the role of epigenetic deregulation in disease pathogenesis. Yet, comprehensive cis-regulatory network descriptions from neuroblastoma tissues are lacking. Here, using genome-wide H3K27ac profiles across 60 neuroblastomas, covering the different clinical and molecular subtypes, we identified four major super enhancer-driven epigenetic subtypes and their underlying master regulatory networks. Three of these subtypes recapitulated known clinical groups, namely MYCN amplified, MYCN non-amplified high-risk and MYCN non-amplified low-risk neuroblastomas. The fourth subtype, exhibiting mesenchymal characteristics, shared features with multipotent Schwann cell precursors, was induced by RAS activation and enriched in relapsed disease. Notably, CCND1, a disease essential gene, was regulated by both mesenchymal and adrenergic regulatory networks converging on distinct super-enhancer modules. Together, this reveals subtype-specific super-enhancer regulation in neuroblastoma.

  Study EGAS00001004552

• Molecular risk stratification in patients with T1 colorectal cancer_WES

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

  Dataset EGAD00001010890

• Integrative analysis of exome-seq, RNA-seq, ATAC-seq (bulk and single-cell), and Hi-C data generated from 3-D spatially mapped samples acquired during surgical resection from 10 patients diagnosed with IDH-WT glioblastoma

  Glioblastoma (GBM), the most common primary brain cancer in adults, remains incurable with no targeted therapies approved despite decades of investigation into its molecular landscape. Treatment failure is attributed to intratumoral heterogeneity in the GBM genome and epigenome, which foster tumor evolution and selection of resistant clones. However, tumor evolution and intratumoral heterogeneity remain poorly understood on the level of the whole tumor as most studies are based on single samples per patient and lack spatial context. Here, we have used 3-D neuronavigation during surgical resection for 10 primary IDH-WT GBM patients to collect 102 samples representing maximal tumor diversity, each mapped by 3-D spatial coordinates. We have applied a strategic set of genomic and epigenomic assays spanning multiple levels of resolution to discover, orthogonally validate, and functionally assess drivers of tumor evolution and intratumoral heterogeneity. These include extrachromosomal DNA amplifications, chromothripsis events, inversions, and translocations that disrupt both the GBM genome and epigenome while generating fusion transcripts and opportunities for therapeutic intervention. We define epigenomic programs that contribute to GBM evolution and intratumoral heterogeneity, revealing their 3-D spatial patterning within whole tumors and their cell type(s) of origin in single-cell data from the same tumor samples. Notably, we distinguish neuronal, glial, and immune programs aberrantly active in tumor cells from their counterparts in normal cells and discover NEUROD1, JUN/FOS, and NF1 transcription factors as key drivers of GBM evolution and growth. Collectively, these data provide unprecedented insight into GBM evolution and intratumoral heterogeneity from single-cell to whole-tumor resolution, redefining current understanding and providing a rich resource of targets for therapeutic investigation.

  Study EGAS00001006785

• 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

• Long-term organoid culture of a small intestinal neuroendocrine tumor rna-seq

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007108

• WXS and RNA-seq for 22 patients treated with radiation + immunotherapy

  Over 500 active clinical trials are investigating combination radiation (RT) and immune checkpoint blockade (ICB) as a cancer treatment; however, the majority of trials have found no positive interaction. We performed a comprehensive molecular analysis of a randomized phase I clinical trial of patients with non-small cell lung cancer (NSCLC) treated with concurrent or sequential ablative RT and ICB. Concurrent treatment was superior to sequential treatment in augmenting local and distant tumor responses and improving overall survival in a subset of patients with highly aneuploid tumors, but not in those with less aneuploid tumors. Our analysis suggested that concurrent treatment eliminated immunologically cold aneuploid tumors better than sequential treatment and improved patient outcomes. In addition, we report that RT alone decreased intratumoral cytotoxic effector T cell and adaptive immune gene expression signatures, in contrast to upregulation of key immune pathways after concurrent administration of RT and ICB. Our findings distinguish the differential genomic and transcriptomic effects of RT versus RT and ICB and challenge the prevailing paradigm that local ablative RT beneficially stimulates the immune response. We propose the use of tumor aneuploidy as a novel biomarker and therapeutic target in personalizing treatment approaches for patients with NSCLC treated with RT and ICB.

  Study EGAS00001006212

• 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

• B cell receptor silencing reveals origin and dependencies of high-grade B cell lymphomas with MYC and BCL2 rearrangements

  High-grade B cell lymphomas with MYC and BCL2 rearrangements (HGBCL-DH-BCL2) represent highly aggressive Germinal center (GC) derived B cell malignancies with dismal prognosis and limited therapeutic options. In the present study, we investigated the expression of BCR complex components, including immunoglobulin (IG) heavy (H) and light (L) chains, and CD79A and CD79B signaling subunits in HGBCL-DH-BCL2, scoring 65% of cases with undetectable IGH immunoreactivity (IGHUND). Transcriptomics (bulk RNA-seq) and genetic (whole exome/genome sequencing) analyses indicated that IGHUND and IGH+ HGBCL-DH-BCL2 lymphomas represent distinct biological entities. Specifically, IGHUND HGBCL-DH-BCL2 exhibit a GC DZ-like molecular program, and select subset-specific gene mutations pointing to distinct evolutionary trajectories respect to the IGH+ counterparts. IG profiling revealed that IGHUND HGBCL-DH-BCL2 preserved potentially productive, hypermutated IGH variable (V) gene rearrangements, indicating a selective advantage conferred to the malignant B cells by residual IGH molecules. Accordingly, disruption of the productively rearranged IGHG1 gene in BCRnull HGBCL-DH-BCL2 cell line model confirmed the dependency on IGH expression for optimal fitness. Using transcriptomics and in situ RNA-scope analysis, we reported that IGHUND HGBCL-DH-BCL2 originated from IG-class switched GC B cells, while IGH-expressing tumors were selected for continuous IGM expression. Moreover, while IGH+ HGBCL-DH-BCL2 favor IGM/IG-Kappa expression, IGHUND counterparts complete IG-Lambda light chain rearrangements. Notably, BCR silencing, caused by accelerated IGH turnover and reduced IGH expression, precedes HGBCL-DH-BCL2 onset, inducing RAG1/2-dependent IG light chain editing and eventually facilitating t(8;22)/IGL::MYC translocations. Finally, IGH silencing observed in HGBCL-DH-BCL2 cell models, which leads to intracellular retention/silencing of CD79B, results in reduced sensitivity to the CD79B-targeting antibody-drug conjugate Polatuzumab-vedotin. Collectively, HGBCL-DH-BCL2 originates from isotype-switched t(14;18)+ GC B cells that undergo IG light chain editing fuelling intra-clonal diversification, BCR silencing, and the acquisition of t(8;22) translocation, while preserving IGH dependency.

  Study EGAS50000001047

• Multiomic Landscape of Multiple Myeloma Precursor and Relapsed Disease

  Multiple myeloma (MM) is a treatable, though incurable, plasma cell malignancy. The molecular mechanisms driving disease onset and the emergence of drug resistance remain elusive. To better characterize the mutational, transcriptional, and epigenetic alterations that accompany MM disease evolution, we accessed molecular data from MM patients treated at Moffitt Cancer Center, with bone marrow biopsies collected between 2011 and 2023.We performed scMultiome (single-cell, paired RNA/ATAC-Seq, 10x Genomics) on 18 CD138+-selected bone marrow aspirate samples, including: 2 healthy donors, 3 patients with monoclonal gammopathy of undetermined significance (MGUS), 4 with smoldering myeloma (SMM), 3 with newly diagnosed MM (NDMM), 1 with early relapse MM (ERMM, 1 to 3 lines of therapy), and 5 with late relapse MM (LRMM, 4 or more lines of therapy). Our dataset includes a sequential sample: a patient with a premalignant condition ("scMultiome_SMM_4") that progressed to active myeloma ("scMultiome_NDMM_3"). These data allow for the assessment of transcriptional and epigenetic dysregulation at the cellular level. Our findings illustrate the dynamic epigenetic and transcriptional landscapes that accompany MM disease progression.We also conducted CUT&Tag analysis on 4 NDMM and 4 LRMM samples to map histone modifications, specifically H3K27ac, to investigate their role in transcriptional reprogramming observed in advanced stages of MM. This analysis revealed sample heterogeneity in terms of super-enhancer-regulated genes associated with active transcription.This cohort also includes samples from the clinical trial NCT04151667, “Daratumumab-Based Response-Adaptive Therapy for Older Adults With Newly Diagnosed Multiple Myeloma.” We performed scRNAseq (single-cell RNA sequencing, 10x Genomics) on 30 CD138+-selected bone marrow aspirate samples and 31 pre-sort samples featuring bone marrow mononuclear cells of various cell types from 5 healthy donors and 11 patients with newly diagnosed MM (NDMM). Our dataset includes sequential samples for all 11 NDMM patients, both prior to induction therapy (denoted by ‘a') with Daratumumab and Dexamethasone, and at cycle 2 day 22 (C2D22, i.e., 60 days after induction therapy start, denoted by ‘b'). For three NDMM patients, we also have samples at relapse (denoted by ‘c'). These data allow for the assessment of transcriptional dysregulation at the single-cell level in both the tumor and immune-tumor microenvironment, associated with Daratumumab treatment in NDMM patients.In addition to the molecular data, demographic, and phenotypic information for all samples is provided.

  Study phs003892

• 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

• 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

• A Clone's Genomic Stability as a Biomarker of Its DNA-Damage Resilience

  A main problem in the treatment of advanced cancers, including gastric cancers and glioblastoma, is the incertitude at which we predict how individual patients will respond to DNA-damaging agents, especially on the long run. Knowing the mechanism behind a patient's response, or the lack thereof, will help us depart from the oversimplified “more-is-better” and “one-size-fits-all” principles according to which DNA-damaging agents are administered. This will improve clinical outcome by allowing us to pinpoint those who would respond better and longer to lower doses than to higher doses of DNA-damaging agents. Under the assumption that the success of DNA-damaging therapy increases with the proliferation rate of a relatively homogeneous tumor population, there was little reason to assume anything other than monotonic dose-response relations. With the recent paradigm shift that most cancers are in fact DNA mosaic products of ongoing evolution, comes the urgency to reconsider these fundamental principles behind DNA-damaging therapy administration. As the developers of one of the first DNA deconvolution methods and with access to technologies to profile the transcriptomes of up to 10,000 cells simultaneously, we are equipped to embark on first personalized dose-finding strategies for DNA-damaging therapies. We will test the potential of the very long-term legacy that DNA-damage entails on a cell “genomic instability” as new biomarker of DNA-damage response. Our preliminary studies showed that, for most cancer types, DNA-damaging agents change a clone's genomic instability and that clones succumb to a limit in the amount of genomic instability they can tolerate. In particular, our results showed that patients with intermediate genomic instability have a very poor outcome and that this relation is only evident among treatment-naive patients, but not among patients treated with DNA-damaging agents. Further they show that we can measure genomic instability per clone and that clones with extreme genomic instability typically don't grow large. Our hypothesis that genomic instability, rather than proliferation rate, determines how sensitive a tumor is to DNA-damaging agents on the long-term, is founded on two unexpected findings. Patients with extremely high genomic instability per tumor clone have an exceptionally good outcome. Aim 1 will integrate exome- and single cell RNA-Seq (scRNA-Seq) data to characterize clones and to measure how much genomic instability they can tolerate. Low genomic instability is associated with reduced benefit from DNA-damaging agents. Aim 2 will use comet assays and treatment history to quantify DNA damage per clone, relating it to the clones' ability to tolerate DNA damage and to changes in the genomic instability of therapy-surviving clones.

  Study phs003762

• 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

• Targeted_sequencing_of_genes_recurrently_mutated_in_AML

  Background Massively parallel sequencing technology has transformed cancer genomics. It is now feasible, in a clinically relevant time-frame, for a clinically manageable cost, to screen DNA from patient tumours for mutations essentially genome-wide. The challenge for personalised medicine will be to increase the sample size to thousands or tens of thousands of well-characterised cases in order to attain sufficient statistical power to stratify patients accurately across the complexity and genomic heterogeneity expected for most of the common tumour types. Currently, whole genome sequencing on this scale is not feasible, and targeted sequencing of relevant portions of the genome will be required. Pilot data We have developed protocols for large-scale, multiplexed sequencing of 100-200 genes in thousands of samples. Essentially, using robotic technology, genomic DNA from the cancer specimen is processed into sequencing libraries with unique DNA barcodes, thereby allowing sequencing reads to be attributed to the sample they derive from. Currently, these sequencing libraries can be generated in a 96-well format using fully automated protocols, and we are exploring methods to expand this to a 384-well format. The sequencing libraries are pooled and hybridized to custom sets of RNA baits representing the genomic regions of interest. Sequencing of the pulled-down libraries is done in pools of 48-96 samples per lane of an Illumina Hi-Seq. This protocol is already implemented at the Sanger Institute. We have published proof that somatic mutations in novel cancer genes can be identified from exome-wide sequencing. In unpublished pilot data, we have established the feasibility of robotic library production, custom pull-down, and multiplexed sequencing of barcoded libraries for 100 known myeloid cancer genes across 760 myelodysplasia samples. Highlights of the data thus far analysed reveal that the coverage is remarkably even between samples; when 96 samples are run, average coverage per lane of sequencing is ~250, with 90-95% of targeted exons covered by >25 reads; known mutations can be discovered in the data set; and the protocol is amenable to whole genome amplified DNA. The bioinformatic algorithms for identification of substitutions and indels in pull-down data are well-established; we have pilot data proving that copy number changes, LOH and genomic rearrangements in specific regions of interest can also be identified by tiling of baits across the relevant loci. Proposal We propose to apply this methodology to 10000 samples from patients with AML enrolled in clinical trials over the last 10-20 years. Oncogenic point mutations and potentially genomic rearrangements will be identified, and linked to clinical outcome data, with a view to undertaking the following sorts of analyses: • Identification of co-occurrence, mutual exclusivity and clusters of driver mutations. • Correlation of prognosis with driver mutations and potentially gene-gene interactions • Exploration of genomic markers of drug response Ultimately, we would like to be in a position to release the mutation data together with matched clinical outcome data to genuine medical researchers via a controlled access approach, possibly within the COSMIC framework (www.sanger.ac.uk/genetics/CGP/cosmic/). The vision here is to generate a portal whereby a clinician faced with an AML patient and his / her mutational profile can obtain a ‘personalised’ prediction of outcome, together with a fair assessment of the uncertainty of the estimate. With a sufficient sample size, there would also be the potential to develop decision support algorithms for therapeutic choices based on such data.

  Study EGAS00001000408

• Neoadjuvant nivolumab or nivolumab plus ipilimumab in early-stage triple negative breast cancer: phase 2 adaptive BELLINI trial. RNA-Seq

  The adaptive BELLINI trial explored whether short combination ICI induces immune activation (primary endpoint: two-fold increase in CD8+ T cells or IFNG), providing rationale for neoadjuvant ICI without chemotherapy. In window of opportunity cohorts A (4 weeks anti-PD1) and B (4 weeks anti-PD1 + anti-CTLA4), we observed immune activation in 53% (8/15) and 60% (9/15) of patients, respectively. High tumor-infiltrating lymphocytes (TILs) correlated with response. Based on these data, we opened cohort C for patients with high TILs (≥50%) who received 6 weeks neoadjuvant anti-PD1 + anti-CTLA4 followed by surgery (primary endpoint: pathological complete response, pCR). 53% (8/15) of patients had major pathological response (< 10% viable tumor) at resection, with 33% (5/15) having pCR. All cohorts met Simon’s two-stage threshold for expansion to stage II. In conclusion, neoadjuvant immunotherapy without chemotherapy demonstrates potential efficacy and warrants further investigation in patients with early TNBC.

  Study EGAS50000000567

• Adipose Tissue Omics In Obesity

  This project used a multi-omics approach to study the adipose tissue transcriptome in humans with and without obesity. Our human cohort included males and females between the ages of 18-80 who underwent elective abdominal surgery at the University of Michigan and the Ann Arbor Veterans Affairs Healthcare System.Patients with advanced cancer or immuncompromised status were excluded. Visceral adipose tissue (VAT/omental) and subcutaneous adipose tissue (SAT/abdominal wall) were collected from subjects with obesity (BMI>=35) at the beginning of elective surgical procedures, placed on ice, and transported back to the laboratory for immediate processing. For adipose tissue stromal cells (ASC) isolation, adipose tissue was digested with Type II collagenase (2 mg/mL in PBS/2% BSA) at 37oC, 60min, centrifuged at 250 rcf, and the stromal-vascular cell pellet was retrieved, plated overnight, and the adherent cells passaged three times to enrich for ASC, which were then frozen in DMEM/F12, 15% fetal calf serum, and 10% DMSO in liquid nitrogen until use. Single nuclear sequencing was performed on preparations from intact adipose tissue. Bulk RNA sequencing and single cell sequencing were performed on ASC isolated from intact adipose tissue and cultured in 2D and 3D culture systems. We currently apply bulk, single cell, and single nuclear sequencing analyses to these tissues. One component of this project studies single nuclear sequencing in nuclei isolated from intact adipose tissues from human subjects with obesity, focusing on depot comparisons (VAT vs. SAT). This cohort also permits comparisons of adipose tissue transcriptomes between male and female subjects, controlling for age and BMI. Primary findings to date include identification of multiple ASC, adipocyte, and endothelial cell (EC) subpopulations, T-cell, and macrophage (MAC) populations that are in concordance with and expand upon other published datasets. Of two dominant ASC subpopulations, one (inflammatory-mesothelial, IM- ASC) is present in VAT, but absent in SAT, while the other (fibroadipogenic, FA-ASC) is present in both VAT and SAT. We identified multiple EC subtypes with features of lymphatic, venous, and arterial EC, and a PRDM16 expressing population with features of an EC progenitor. Immune cell populations match recent experimental validation of lipid-activated macrophage (LAM) phenotypes, TIM4 macrophages, and a prominent population of MRC1/CD206+ resident AT macrophages with expression signatures related to glucocorticoid activation. These data support depot-specific differences in ASC subpopulations and previously unrecognized EC heterogeneity within human AT, with the potential to shed light on functional differences in AT depot and disease state. Planned future studies include ATAC-Seq and other omics analyses. Another component of this project studies transcriptional responses of ASC in vitro to various stimuli, including different matrix environments. Current data includes bulk RNA Sequencing of human ASC from VAT and SAT differentiated into mature adipocytes in hydrogels of differing stiffness (3%, 5%, 10%) to define transcriptional responses to matrix environments of different stiffness. Primary findings to date include enrichment of adipogenesis and lipid and oxidative metabolic gene pathways in adipocytes in 5% relative to 3% hydrogels, while fibrotic gene pathways have been enriched in 3% hydrogels. These data demonstrate that intermediate density matrix promotes a more adipogenic, less fibrotic adipocyte phenotype geared towards increased lipid and aerobic metabolism. Planned future studies include expansion of intact tissues from our cohort studied with single nuclear sequencing, as well as further study of ASC responses in vitro to different biologic stimuli, including different matrix environments.

  Study phs003390

• RNA-seq of DF149 cells – a patient-derived xenograft model of ascites-derived, homologous recombination repair-proficient, high-grade serous ovarian carcinoma – cultured in vitro and isolated after 8 hours treatment with DMSO control (3 x biological replicates) and 2.5 µM CBL0137 (3 x biological replicates)

  Background: High-grade serous ovarian carcinomas (HGSCs) are a heterogeneous subtype of epithelial ovarian cancers and include serous cancers arising in the fallopian tube and peritoneum. These cancers are now subdivided into homologous recombination repair (HR)-deficient and proficient subgroups as this classification impacts on management and prognosis. PARP inhibitors (PARPi) have shown significant clinical efficacy, particularly as maintenance therapy following response to platinum based chemotherapy in BRCA-mutant or homologous recombination (HR)-deficient HGSCs in both the 1st and 2nd line settings. However, PARPi have limited clinical benefit in HR-proficient HGSCs which make up almost 50% of HGSC and improving outcomes in these patients is now a high priority due to the poor prognosis with current standard of care. There are a number of potential lines of investigation including efforts in sensitizing HR-proficient tumors to PARPi. Herein, we aimed to develop a novel combination therapy by targeting SSRP1 using a small molecule inhibitor CBL0137 with PARPi in HR-proficient HGSOCs. Experimental design: We tested anti-cancer activity of CBL0137 monotherapy using a panel of HGSOC cell lines and patient-derived tumor cells in vitro. RNA sequencing was used to map global transcriptomic changes in CBL0137-treated patient-derived HR-proficient HGSOC cells. We tested efficacy of CBL0137 in combination with PARPi using HGSOC cell lines and patient-derived tumor cells in vitro and in vivo. Results: We show that SSRP1 inhibition using a small molecule, CBL0137, that traps SSRP1 onto chromatin, exerts a significant anti-growth activity in vitro against HGSC cell lines and patient-derived tumor cells, and also reduces tumor burden in vivo. CBL0137 induced DNA repair deficiency via inhibition of the HR repair pathway and sensitized SSRP1-high HR-proficient HGSC cell lines and patient-derived tumor cells/xenografts to the PARPi, Olaparib in vitro and in vivo. CBL0137 also enhanced the efficacy of DNA damaging platinum-based chemotherapy in HGSC patient-derived xenografts. Conclusion: Our findings strongly suggest that combination of CBL0137 and PARP inhibition represents a novel therapeutic strategy for HR-proficient HGSCs that express high levels of SSRP1 and should be investigated in the clinic.

  Study EGAS00001006662

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

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

  Dataset EGAD00001004266

• Targeted sequencing of genes recurrently mutated in AML

  Background Massively parallel sequencing technology has transformed cancer genomics. It is now feasible, in a clinically relevant time-frame, for a clinically manageable cost, to screen DNA from patient tumours for mutations essentially genome-wide. The challenge for personalised medicine will be to increase the sample size to thousands or tens of thousands of well-characterised cases in order to attain sufficient statistical power to stratify patients accurately across the complexity and genomic heterogeneity expected for most of the common tumour types. Currently, whole genome sequencing on this scale is not feasible, and targeted sequencing of relevant portions of the genome will be required. Pilot data We have developed protocols for large-scale, multiplexed sequencing of 100-200 genes in thousands of samples. Essentially, using robotic technology, genomic DNA from the cancer specimen is processed into sequencing libraries with unique DNA barcodes, thereby allowing sequencing reads to be attributed to the sample they derive from. Currently, these sequencing libraries can be generated in a 96-well format using fully automated protocols, and we are exploring methods to expand this to a 384-well format. The sequencing libraries are pooled and hybridized to custom sets of RNA baits representing the genomic regions of interest. Sequencing of the pulled-down libraries is done in pools of 48-96 samples per lane of an Illumina Hi-Seq. This protocol is already implemented at the Sanger Institute. We have published proof that somatic mutations in novel cancer genes can be identified from exome-wide sequencing. In unpublished pilot data, we have established the feasibility of robotic library production, custom pull-down, and multiplexed sequencing of barcoded libraries for 100 known myeloid cancer genes across 760 myelodysplasia samples. Highlights of the data thus far analysed reveal that the coverage is remarkably even between samples; when 96 samples are run, average coverage per lane of sequencing is ~250, with 90-95% of targeted exons covered by &gt;25 reads; known mutations can be discovered in the data set; and the protocol is amenable to whole genome amplified DNA. The bioinformatic algorithms for identification of substitutions and indels in pull-down data are well-established; we have pilot data proving that copy number changes, LOH and genomic rearrangements in specific regions of interest can also be identified by tiling of baits across the relevant loci. Proposal We propose to apply this methodology to 10000 samples from patients with AML enrolled in clinical trials over the last 10-20 years. Oncogenic point mutations and potentially genomic rearrangements will be identified, and linked to clinical outcome data, with a view to undertaking the following sorts of analyses: ? Identification of co-occurrence, mutual exclusivity and clusters of driver mutations. ? Correlation of prognosis with driver mutations and potentially gene-gene interactions ? Exploration of genomic markers of drug response Ultimately, we would like to be in a position to release the mutation data together with matched clinical outcome data to genuine medical researchers via a controlled access approach, possibly within the COSMIC framework (www.sanger.ac.uk/genetics/CGP/cosmic/). The vision here is to generate a portal whereby a clinician faced with an AML patient and his / her mutational profile can obtain a ?personalised? prediction of outcome, together with a fair assessment of the uncertainty of the estimate. With a sufficient sample size, there would also be the potential to develop decision support algorithms for therapeutic choices based on such data.

  Dataset EGAD00001000606

• Raw sequencing files from WGS and RNA-Seq

  Acute leukemias are the most common malignancies in children with an onset of the disease among 0-6 years-old. The early development of theses malignancies has been proposed as an indicator of their prenatal origin, hypothesis that has been confirmed especially in the case of acute lymphoblastic leukemia. For acute myeloid leukemia (AML), comprising an heterogenous group of different cytogenetics subtypes, this prenatal origin has also been confirmed specially for MLL-rearranged and t(8;21) subgroups. AMLs carrying t(7;12) constitute up to 30% of cases under 2-years-old, frequently being associated with a trisomy of chromosome 19 and an overexpression of MNX1 gene. An in utero origin of this specific subgroup has been suggested but not clearly confirmed. Here, we present the case of a 5-months-old AML patient with t(7;12) in which we have fully characterized a novel breakpoint occurring among ETV6 and NOM1 genes. We demonstrate the presence of this specific ETV6-NOM1 fusion in CD34+ hematopoietic stem and progenitor cells isolated from the cord blood unit from the same patient and stored at birth. Moreover, trisomy +19 present at diagnosis was not detected at birth. Together, our results demonstrate the prenatal origin of the t(7;12) in infant AML, and establish trisomy +19 as a secondary event in this leukemia.

  Study EGAS50000000186

• A Platform Study of Combination Immunotherapy for the Neoadjuvant and Adjuvant Treatment of Patients with Surgically Resectable Adenocarcinoma of the Pancreas

  This is a single-institutional, randomized, open label platform clinical trial testing various immunotherapy combinations. Patients with newly diagnosed and surgically resectable pancreatic ductal adenocarcinoma (PDA) and scheduled to undergo the Whipple procedure at Johns Hopkins Hospital were eligible to participate in this study. Criteria for determining resectability strictly followed National Comprehensive Cancer Network (NCCN) guidelines. The study consisted of 6 parts: parts 1-5 constituted the Prime Phase (see below), and part 6 was the Extended Treatment Phase. Part 1: Participants received one cycle of immunotherapy two weeks prior to undergoing the Whipple procedure. Part 2: Subjects underwent the Whipple procedure. Part 3: Subjects received a second cycle of immunotherapy 4-10 weeks following the Whipple procedure (2-4 weeks prior to adjuvant chemoradiation). Part 4: Subjects underwent chemoradiation. Part 5: Subjects received 4 additional 28-day cycles of immunotherapy beginning 1-2 months after completing chemoradiation for a total of six cycles (two before chemoradiation and 4 following chemoradiation). Eligible subjects were randomized to Arms A and B in a 1:1 ratio, stratified up-front by age (&#8804;65, >65). Arm A participants received Cy/GVAX. Arm B participants received the combination of Cy/GVAX and nivolumab. GVAX Vaccine consists of equal numbers (2.5 x 108 each) of Panc 6.03 pcDNA1GM-CSF and Panc 10.05 pcDNA1GM-CSF combined into a single vaccination. Each vaccination consisted of six total intradermal injections. A single intravenous (IV) dose of 200 mg/m2 Cy (CytoxanR) was given one day prior to GVAX. Anti-PD-1 Therapeutic Antibody Nivolumab 480 mg every 4 weeks was administered as an intravenous infusion. The primary objective for Arms A and B was to compare IL17 expression (Th17) in vaccine-induced lymphoid aggregates between resected PDAs from patients treated with the combination of Cy/GVAX and anti-PD-1 blockade antibody vs. the treatment of Cy/GVAX alone. Secondary objectives were to assess the safety, overall survival (OS), and disease free survival (DFS) of patients treated with each of the study drug combinations, and to assess the effects of each of the immunotherapy study drug combinations on PD-L1/PD-1 associated pathways, vaccine-induced immune regulatory signatures, and peripheral and intratumoral antigen specific T cell responses. The exploratory objective was to explore the effects of therapy on tumor and peripheral blood and tumor infiltrating immune cells, and to explore potential molecular determinants of response, progression and disease stability. To these ends, tumor infiltrating immune cells (TIICs) were obtained from surgically resected tumors following the neoadjuvant treatment of the study immunotherapy. TIICs were sorted into CD3+CD4+ T cells, CD3+CD8+ T cells, CD3-CD19+ B cells, and CD3-CD19-CD11b+ myeloid cells by flow cytometry. The sorted immune cell subtype pellets were sent to the commercial vendor (Eurofins) for RNA extraction and sequencing. SMART-seq v4 Ultra Low Input mRNA Sample Prep Kit was used for library preparation with 100bp paired end reads. FASTQ files, that were returned from the vendor, are deposited here.

  Study phs003002