516 results for "cancer rna-seq"

in 47.80 milliseconds.

• A Phase 2 Study of Lamivudine in Patients with p53 Mutant Metastatic Colorectal Cancer

  This was an open-label, single-arm Phase II study of lamivudine in patients who had progressed on systemic therapy for advanced colorectal cancer with TP53 mutations. The phase II study had a two-stage design, with a target accrual of 20 evaluable patients for the first stage and a total of 32 patients for the whole study. The first 9 patients were treated with lamivudine at 150 mg PO bid continuously for 28-day cycles. Subsequent patients (10-32) received a higher dose of 600 mg PO bid continuously for 28-day cycles. Tumor assessments were performed every 8 weeks until documented disease progression by Response Evaluation Criteria in Solid Tumors (RECIST) or drug intolerance. Whole genome sequencing (WGS) and total RNA-seq were performed on pre- and post- treatment biopsies on this trial.

  Study phs002833

• Epigenetic profiling of colorectal cancer initiating cells (CC-ICs) to identify bivalently marked genes (H3K4me3 and H3K27me3 ChIP-seq), and investigation of changes in transcriptome following EZH2 inhibition using RNA-seq.

  Epigenetic regulation of transcription plays a crucial role in lineage commitment of embryonic stem cells. Promoters of key lineage-specific differentiation genes are found in a repressed bivalent state, having both activating H3K4me3 and repressive H3K27me3 histone marks, making them poised for transcription upon loss of H3K27me3 in response to environmental cues. Whether the tumour-initiating, self-renewing, cancer-initiating cells (C-ICs) have similar epigenetic regulatory mechanism that prevent lineage commitment is unknown. In order to investigate bivalently marked and repressed promoters, we used a patient-derived CC-IC enriched model to identify the changes in transcriptome following inhibition of EZH2, the H3K27 methyltransferase. We also performed ChIP-seq for H3K27me3 and H3K4me3 at baseline in order to identify repressed and bivalently marked promoters.

  Study EGAS00001003003

• Genomic and Immune Profiling of Breast Cancer Brain Metastases

  Previous work has shown that breast cancer brain metastases (BCBrM) can have unique molecular profiles compared to extracranial tumors (ECTs), resulting in some patients demonstrating subtype switching between primary/extracranial tumors and brain metastases. Description of the biology of BCBrM and associations between biological features and clinical outcomes, particularly between molecular intrinsic subtypes, will facilitate the design of rational therapies for patients with BCBrM. To better characterize the molecular profiles of BCBrM, we conducted whole exome sequencing (WES) and total RNA sequencing (RNA-Seq) on BCBrM tissues, regardless of clinical or molecular subtype, obtained from a cohort of 42 patients from Duke University Medical Center. This cohort adds valuable patient samples to the field and enables strengthened investigation of somatic mutations, copy number alterations, gene expression profiles, and immune populations within and between BCBrM patients and other similar cohorts. The Duke BCBrM cohort includes clinical and molecular annotation of the 42 patient specimens and their associated WES and RNA-Seq data. Clinical and molecular (PAM50) intrinsic subtypes in ECTs and brain metastases were largely concordant within patients. There were observed differences in intrinsic subtype calls between frozen and FFPE specimens of the same tissue, with discordance mostly observed with LumB calls. Copy number alterations in chromosomal regions included clinically relevant genes, and somatic alterations were largely conserved within patients across matched samples, with few genes, including TP53, being recurrently altered in the majority of patient specimens. In matched FFPE brain metastasis and ECT samples, we observed differentially expressed genes and pathways that had increased or decreased expression in brain metastasis, including in oxidative phosphorylation and DNA repair. Survival outcomes for patients differed by intrinsic molecular subtypes, with some association between expression of select gene pathways and survival. Data available through dbGaP will be the sequencing data and sample attributes from BCBrM tissues obtained from 42 patients from Duke University Medical Center.

  Study phs003673

• Exploring_the_heterogeneity_of_sarcoma_using_single_cell_sequencing_

  Multi region samples are collected from patients, with consent, immediately after resection of the tumour. Samples are digested and sorted using FACS as single cells into lysis buffer. Cells are then stored until further processing for G&T-seq. After sequencing, we will explore intra-tumour heterogeneity using computational approaches to integrate RNA and DNA data onto the tumour phylogeny This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002866

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

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

  Study EGAS50000001401

• Transcriptional Profiling of Macrophages In Situ in Metastatic Melanoma Reveals Localization-Dependent Phenotypes and Function

  Therapeutic interventions modulating immune function at the tumor site could improve outcomes in cancer. Here we analyzed metastatic melanoma, a tumor type highly responsive to T-cell based therapies, and found that the tumor-infiltrating T cells are localized closely to CD14+ monocytes/macrophages rather than to melanoma cells. Using novel immunofluorescence-guided laser capture micro-dissection, we analyzed transcriptomes of CD3+ T cells, CD14+ monocytes/macrophages and melanoma cells in non-dissociated tissue. CD14+ cells localized in tumors beds displayed specific transcriptional signatures distinct from CD14+ cells that were localized in tumor stroma.The 8 distinct metastatic melanoma samples are from the Cooperative Human Tissue Network (CHTN). 64 RNA-seq files are available from 16 intraCD14+, 16 stromaCD14+, 16 cancer cells, 8 intraCD3+ and 8 stromaCD3+ cells. 2 LCM harvests were completed for the CD14+ and cancer cells on two sections of the samples. One harvest was done for the CD3+ samples.The processed data (TPM/raw counts tables) are available on GEO under GSE180124. 

  Study phs002564

• Dissecting intratumor heterogeneity of nodal B cell lymphoma on the transcriptional, genetic, and drug response level

  Tumor heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is well-known to affect the efficacy of anti-cancer, most personalized treatment approaches do not account for intratumor heterogeneity. We addressed this issue by studying the heterogeneity of lymph node-derived B cell non-Hodgkin lymphoma (B-NHL) by single cell RNA-sequencing (scRNA-seq) and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subclones and compared their drug response and genomic profiles. Malignant subclones of the same patient responded strikingly different to anti-cancer drugs ex vivo, which recapitulated subclone-specific drug sensitivity during in vivo treatment. Tumor infiltrating T cells represented the majority of non-malignant cells, whose gene expression signatures were similar across all donors, whereas the frequencies of T cell subsets varied significantly between the donors. Our data provide new insights into the heterogeneity of B-NHL and highlight the relevance of intratumor heterogeneity for personalized cancer therapies.

  Study EGAS00001004335

• Development of the prevention and therapy of CRC using patient derived culture tissues.

  Tumor tissues consist of functionally heterogeneous cells including those with high clonogenic potential called cancer stem cells (CSCs) which produce differentiated non-tumorigenic progenies. Transplantation assays have been performed to identify cell surface markers in solid tumors including colorectal cancers, which identified several CSC markers such as CD133 and CD44. In addition to the transplantation, the lineage tracing of genetically labelling cells in animal models of cancer qhas been employed to explore the cellular heterogeneity of tumors by identifying the CSCs and its progenies. This approach is particularly useful to investigate solid tumors, where cellular proliferation and differentiation are controlled by niche factors and cell-cell contact with adjacent cells. Lineage tracing of mouse model of cancer unravel the mode of growth of autochthonous tumor in their native environment. In this study, we used patient derived organoids (PDO) established from surgically resected CRC patients to identify the CSC. Single cell RNA-seq analysis identified potential CRC markers, which may exhibit self-renewal activity and produce progenitors

  Study JGAS000139

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

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

  Study EGAS00001003452

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

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

  Study EGAS00001001754

• Epigenome and transcriptome profiling of chronic lymphocytic leukemia patients

  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.

  Study EGAS00001001821

• Integrative analysis reveals a macrophage-predominant, immunosuppressive immune microenvironment and subtype-specific therapeutic vulnerabilities in advanced salivary gland cancer

  Background Salivary gland cancers (SGC) are rare and heterogeneous malignant tumors. Advanced SGC lack established treatment options and show poor response to immunotherapy. Here, an integrative multi-omics analysis in a large cohort of advanced SGC revealed insights into the tumor immune microenvironment (TIM) and distinct mechanisms of immune evasion. Methods A total of 104 patients with recurrent/metastatic SGC from the DKTK MASTER program were included in this study. Whole-exome or whole-genome sequencing and RNA-sequencing was performed on fresh frozen tumor tissue. The tumor immune microenvironment was analyzed using CIBERSORT deconvolution analysis and immune gene expression scores in bulk RNA-sequencing data. Single-nuclei sequencing and immunohistochemistry analyses were performed in selected samples. Results were validated in bulk RNA-sequencing data of a previously published independent dataset. Results Bulk transcriptome analysis revealed an immune-deserted TIM in the majority of advanced SGC samples. Immune exclusion was most prominent in adenoid cystic carcinoma (ACC) subgroup 1 exhibiting a downregulation of the antigen processing machinery. Only a small subset of advanced SGC, including few adenoid cystic carcinoma, exhibited T-cell inflammation, which was correlated with tumor mutational burden in Non-ACC samples. Subtype specific expression of immune checkpoints as well as cancer testis antigens were identified with prominent expression of VTCN1 in luminal cells within ACC. Single-cell RNA-sequencing and bulk RNA-seq deconvolution analysis validated immune cell exclusion and revealed a TIM that was dominated by M2 macrophages across SGC subtypes. Among evaluable patients treated with immune checkpoint inhibitors, a high T-cell to macrophage ratio was associated with clinical benefit. Conclusions These data support biomarker-based development of immune-checkpoint inhibition and the development of novel immune-checkpoint inhibitors and cellular therapies in SGC.

  Study EGAS50000000809

• Human Pilocytic Astrocytoma Single Cell RNA Sequencing

  Pilocytic astrocytoma (PA), the most common childhood brain tumor, is a low-grade glioma with a single driver BRAF rearrangement. Here, we perform scRNAseq in six PAs using methods that enabled detection of the rearrangement. When compared to higher-grade gliomas, a strikingly higher proportion of the PA cancer cells exhibit a differentiated, astrocyte-like phenotype. A smaller proportion of cells exhibit a progenitor-like phenotype with evidence of proliferation. These express a mitogen-activated protein kinase (MAPK) program that was absent from higher-grade gliomas. Immune cells, especially microglia, comprise 40% of all cells in the PAs and account for differences in bulk expression profiles between tumor locations and subtypes. These data indicate that MAPK signaling is restricted to relatively undifferentiated cancer cells in PA, with implications for investigational therapies directed at this pathway. Note that 931 out of 1234 samples in the Sequence Read Archive (SRA) passed QC and have transcript count data available from the Broad Institute Single Cell Portal. The transcript count data can be found at https://singlecell.broadinstitute.org/single_cell/study/SCP271/pilocytic-astrocytoma-single-cell-rna-seq#study-download

  Study phs001854

• International Cancer Proteogenome Consortium (ICPC): Proteogenomics of Oral Squamous Cell Carcinoma in Taiwan

  Oral cavity squamous cell carcinoma (OSCC) is the major subtype (~90%) of oral cancer, which is leading cancer with high mortality worldwide. The environmental risk factors associated with OSCC are quite different in Taiwan as compared to other populations. In most non-Taiwanese populations, OSCC is often associated with tobacco chewing, cigarette smoking, alcohol drinking, or human papillomavirus (HPV) infection. In Taiwan, in contrast, the OSCC occurs mainly in males aged 30 years or older who have habitual betel nut chewing, cigarette smoking, and alcohol drinking. To uncover molecular mechanisms and identify new therapeutic strategies of oral squamous cell carcinoma, this study recruited 112 treatment-naïve OSCC patients and generated multi-omics data, including WES, RNA-Seq, TMT-labeled proteomes, and phosphoproteomes from paired normal/tumor tissues. We have investigated our multi-omics data, covering molecular characterizations such as mutation signatures, impacts of mutations on RNA/protein/phosphoprotein, RNA-protein correlation, somatic copy number alteration, and the relationship with RNA & protein, proteome, and phosphoproteome analysis (T vs NAT), novel peptides (somatic variants, neoantigens, CT antigens, T cell receptors, etc.), as well as molecular subtyping (& pathway analysis). Notably, we have discovered a higher frequency of APOBEC3B deletion polymorphism in the Taiwanese population, and APOBEC3A is a potential prognostic biomarker. Our proteogenomic data enable us to further stratify patients into four groups based on APOBE3A expression level, which is found inversely correlated with EGFR expression. Therefore, the patients with high EGFR and low APOBEC3A expression may benefit from anti-EGFR therapy, while those with low EGFR and high APOBEC3A expression may benefit from immunotherapy. To clarify if there is an interaction between genetic polymorphisms and betel nut chewing habit in Taiwanese OSCC, we plan to concurrently evaluate the effect of habitual betel nut chewing on the gene mutation pattern, the dysregulated protein/phosphoproteins as novel biomarkers and/or therapeutic targets, and the molecular subtyping of patients for effective treatments. 

  Study phs002580

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

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

  Study EGAS00001003554

• Case Report: Pre-Clinical Combination Targeting VEGF and PI3K in a Rare, Aggressive Mixed Endometrial Carcinoma

  We report a rare case of a young patient (VENUS 167) initially diagnosed with grade 1 endometrioid endometrial cancer, which, following endocrine treatment, presented with mixed aggressive carcinoma with three distinct histologic patterns: grade 1 endometrioid, large cell neuroendocrine, and undifferentiated carcinoma. The surgical specimen at the time of disease progression was used to establish OPTO.85, a patient-derived organoid (PDO), followed by a corresponding organoid-derived xenograft (ODX). Multi-omic analyses confirmed that OPTO.85 accurately reflected the patient’s tumor characteristics. WES analysis identified oncogenic alterations in PIK3CA, ARID1A, and CTNNB1. Further RNA-seq and ATAC-seq analyses revealed enrichment in VEGF and Wnt signaling pathways, suggesting therapeutic vulnerabilities. A high-throughput drug screen was conducted using ApexBio-Approved and Epigenetic Drug Libraries, along with Kinase Inhibitor and Tool Compound Libraries developed at the Ontario Institute of Cancer Research. The OPTO.85 PDO exhibited sensitivity to PI3K inhibitors and responsiveness to VEGF inhibition. Cediranib demonstrated synergy with BKM120, significantly reducing organoid growth. This combination also showed in vivo efficacy in the ODX model, where dual inhibitors significantly suppressed tumor growth compared to single compounds. This case exemplifies the impact of genomic profiling and patient-derived models in identifying actionable molecular changes in rare cancers with limited therapeutic options and poor prognosis. It highlights that high-throughput sequencing for individual patient tumors and generation of patient-derived models are feasible in endometrial cancer. This preclinical model may assist clinical decision and personalized therapy requiring validation in prospective studies.

  Study EGAS50000001665

• Breast Cancer FRT RNA seq

  We propose to definitively characterise the somatic genetics of ER+ve, HER2-ve breast cancer through generation of comprehensive catalogues of somatic mutations in breast cancer cases by high coverage genome sequencing coupled with integrated transcriptomic and methylation analyses.

  Dataset EGAD00001000338

• Data Access Committee for the DNA sequencing data included in the study "Multifocal cohort analysis unveils cell types associated with regional lymph node seeding in prostate cancer".

  A data transfer agreement (DTA), based on the example DTA provided by UZ Gent, is to be filled out by the requester and Multifocal RNA-seq and tNGS data from primary and metastatic tissues of locally advanced and de novo metastatic cohorts - Data Access Committee. https://www.uzgent.be/nl/over-uz-gent/organisatie/data-access-comite/wat-zijn-de-opdrachten-van-het-data-access-comite Coordinator: Muriel Fouquet E-mail: dataaccesscommittee@uzgent.be

  Dac EGAC50000000930

• Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma - Ribo-Seq

  This study aims to investigate the dysregulation of RNA translation and identify functional non-canonical open reading frames (ORFs) as potential targets for medulloblastoma treatment. The study involves ribosome profiling and RNAseq of medulloblastoma tissues and cell lines to observe the translation of non-canonical ORFs. Multiple CRISPR-Cas9 screens will be used to identify functional non-canonical ORFs implicated in medulloblastoma cell survival.

  Dataset EGAD00001011193

• Comprehensive single cell study of lung adenocarcinoma from early to metastatic stages

  We performed single cell RNA sequencing (scRNA-seq) for 208,506 cells derived from 58 lung adenocarcinomas from 44 patients, which covers primary tumour, lymph node and brain metastases, and pleural effusion in addition to normal lung tissues and lymph nodes. The extensive single cell profiles depicted a complex cellular atlas and dynamics during lung adenocarcinoma progression which includes cancer, stromal, and immune cells in the surrounding tumor microenvironments.

  Dataset EGAD00001005054

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

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

  Study JGAS000105

• 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

• 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

• Childhood Cancer Data Initiative (CCDI): Pediatric In Vivo Testing Program –Sarcoma, Kidney, and Liver Cancers

  The goal was to complete whole exome sequencing (WES) and RNA sequencing (RNA-Seq) on patient derived xenograft samples. Cancer in children is rare with approximately 15,700 new cases diagnosed annually in children 21 years or younger in the U.S. Through use of multimodality therapy (surgery, radiation therapy, and aggressive chemotherapy), 70% of patients will be 'cured' of their disease, and 5-year Event-Free Survival (EFS) exceeds 80%. Consequently, the number of patients that can be enrolled in phase I/II clinical trials is small, and most patients will have been extensively treated, hence drug/radiation resistant. Thus, preclinical studies that accurately translate into effective clinical therapy are an essential component of pediatric drug development. Our group has contributed to studies in the Pediatric Preclinical Testing Program (PPTP) and Pediatric Preclinical Testing Consortium (PPTC) that have led to clinical studies through Children's Oncology Group (COG). Of importance, we have developed and characterized over 330 Patient Derived Xenografts (PDX), developed from tumors both at diagnosis and relapse, that can be used to facilitate pediatric drug development as directed by FDA under the Research to Accelerate Cures and Equity for Children Act (RACE for Children Act).

  Study phs003160

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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