Search Results - EGA European Genome-Phenome Archive

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

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

Study phs001340

Generation of a Cellular Atlas of Human Adipose Tissue

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

Study phs002766

Genomics and Epigenomics of the Elderly Response to Pneumococcal Vaccines

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

Study phs002361

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

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

Study EGAS00001005478

3D tissue engineered human skeletal muscle modelling Facioscapulohumeral Muscular Dystrophy

acioscapulohumeral muscular dystrophy (FSHD) is caused by sporadic misexpression of the transcription factor double homeobox 4 (DUX4) in skeletal muscles. So far, monolayer cultures and animal models have been used to study the FSHD disease mechanism and for FSHD therapy development, but these models do not fully recapitulate the disease and there is a lack of knowledge on how DUX4 misexpression leads to skeletal muscle dysfunction. To overcome these barriers, we have developed a 3D tissue engineered skeletal muscle (3D-TESM) model by generating genetically matched myogenic progenitors (MPs) from human induced pluripotent stem cells of three mosaic FSHD patients. 3D-TESMs derived from genetically affected MPs recapitulate pathological features including DUX4 and DUX4 target gene expression, smaller myofiber diameters, and reduced absolute forces upon electrical stimulation. RNA sequencing data illustrates increased expression of DUX4 target genes in 3D-TESMs compared to 2D myotubes, and cellular differentiation was improved by 3D culture conditions. Treatment of 3D-TESMs with small molecules identified in drug development screens in 2D muscle culture showed no improvements in contractile force and sarcomere organization, highlighting the potential benefit of 3D-TESMs for preclinical drug potency screening by eliminating false-positive compounds for FSHD research and drug development.

Study EGAS50000000502

Pediatric Non-Down Syndrome Acute Megakaryoblastic Leukemia is Characterized by Distinct Genomic Subsets with Varying Outcomes

Acute Megakaryoblastic Leukemia (AMKL) is a subtype of acute myeloid leukemia (AML) in which cells morphologically resemble abnormal megakaryoblasts. While extremely rare in adults, AMKL accounts for 4-15% of newly diagnosed childhood AML1-3. A significant proportion of pediatric AMKL cases occur in children with Down syndrome (DS). These patients have excellent outcomes and are characterized at the genomic level by a founding GATA1 mutation4-6. In contrast, AMKL in patients without DS (non-DS-AMKL) is frequently associated with poor outcomes. Previous efforts have identified chimeric oncogenes in a significant number of cases, including RBM15-MKL1, CBFA2T3-GLIS2, KMT2A gene rearrangements and NUP98-KDM5A7-9. The etiology of 30-40% of cases, however, remains unknown. To better understand the genomic landscape of non-DS-AMKL, we performed RNA and exome sequencing on specimens from 99 patients (75 pediatric and 24 adult). We demonstrate that pediatric non-DS-AMKL is a heterogeneous malignancy that can be divided into seven subgroups with varying outcomes. These subgroups are characterized by chimeric oncogenes with cooperating mutations in epigenetic and kinase signaling genes. Overall, these data shed light on the etiology of AMKL and provide useful information for treatment tailoring.

Study EGAS00001002183

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

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

Study EGAS00001004290

Clonal_haematopoiesis_in_patients_with_AAA

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

Study EGAS00001002873

Cell type mapping of inflammatory muscle diseases highlights selective myofiber vulnerability in inclusion body myositis

Inclusion body myositis (IBM) is the most prevalent inflammatory muscle disease in older adults with no effective therapy available. In contrast to other inflammatory myopathies like subacute immune-mediated necrotizing myopathy (IMNM), IBM follows a chronic disease course with both inflammatory and degenerative features of pathology. Moreover, causal factors and molecular drivers of IBM progression are largely unknown. Therefore, we paired single-nucleus RNA sequencing with spatial transcriptomics from patient muscle biopsies to map cell type-specific drivers underlying IBM pathogenesis compared to IMNM muscles and non-inflammatory skeletal muscle samples. In IBM muscles, we observed a selective loss of type 2 myonuclei paralleled by increased levels of cytotoxic T and cDC1 cells. IBM myofibers were either characterized by upregulation of cell stress markers featuring GADD45A and NORAD, or protein degradation markers including RNF7 associated with p62 aggregates. GADD45A upregulation was preferentially seen in type 2A myofibers associated with severe tissue inflammation. We also noted IBM-specific upregulation of ACHE encoding acetylcholinesterase, which can be regulated by NORAD activity and result in functional denervation of myofibers. Our results provide promising insights into possible mechanisms of myofiber degeneration in IBM and suggest a selective type 2 fiber vulnerability linked to genomic stress and denervation pathways.

Study EGAS50000000310

Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

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

Study EGAS00001003438

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

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

Study EGAS00001004291

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

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

Study EGAS00001004292

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

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

Study EGAS00001006678

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

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

Study EGAS00001006864

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

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

Study EGAS00001007254

Investigating Delayed-Onset Drug Hypersensitivity Reactions Prospectively

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

Study phs003344

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

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

Study EGAS00001003481

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

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

Study phs001436

Exceptional Responders Initiative

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

Study phs001145

Cell and Circuit-Specific Exploration of HIV Neurogenomics in Context of Opiate and Cocaine Misuse

Utilizing frozen postmortem brain tissue from the Manhattan HIV Brain Bank, we performed cell-type-specific, genome-wide sequencing of three different types: 1) single-nucleus RNA-sequencing (snRNA-seq), 2) in situ Hi-C, and 3) HIV integration site sequencing (IS-seq). We assayed tissue from three different donor types: 1) HIV-uninfected control (HIV-), HIV-infected (HIV+), and HIV-infected with HIV encephalitis (HIVE), who served as a positive control for active CNS viral infection. The diagnosis of HIVE was made by a board-certified neuropathologist on routine histopathology, based on the presence of a microglial nodule encephalitis with characteristic multinucleated syncytial cells and staining for HIV p24 antigen. Other HIV-infected individuals (the HIV+ group) did not have this histologic evidence of active HIV brain replication, and were variable with regard to cART status and virologic control prior to death. Donor characteristics, including ante-mortem plasma HIV loads, CD4 counts, and cART status, were documented. Quality controls were used for each assay: For snRNA-seq we applied filters to remove low quality nuclei from analysis and performed experiments in which mouse tissue was mixed with human tissue to determine a reliable HIV read count threshold. For Hi-C, we used paired HIV- and HIVE samples and processed all samples in parallel, at the same time. Correlation analysis was used to ensure biological replicates were aligned. For IS-seq, each time the assay was run, a positive and negative control were included. This allowed us to ensure a lack of contamination in IS-seq libraries and that, when a particular sample failed to yield library, this was due to low levels of HIV rather than a failure of the assay. Altogether, we built an integrative dataset from 79 sequencing files, including snRNA-seq, and cell type specific Hi-C, IS-seq, and ChIP-seq datasets. Human autopsy brain samples were collected by the Manhattan HIV Brain Bank (MHBB, U24MH100931), using protocols under the supervision of the Icahn School of Medicine at Mount Sinai (ISMMS) Institutional Review Board. Written informed consent was obtained either from decedents or their primary next of kin. Major Findings: Reorganization of open/repressive (A/B) compartment structures in HIVE microglia encompassing 6.4% of the genome was linked to transcriptional activation of interferon (IFN) signaling and cell migratory pathways and was partially recapitulated by IFN-gamma stimulation of cultured microglia. In contrast, decreased expression and repressive compartmentalization of genes regulating neuronal health and signaling was seen in both HIVE and HIV+ microglia. IS-seq recovered 1,221 integration sites in the brain that displayed distinct genomic patterns as compared to peripheral lymphocyte integration and were enriched for chromosomal domains newly mobilized into a permissive chromatin environment in HIVE microglia. Viral transcription occurred in a subset of highly activated microglia comprising 0.003% of all nuclei in HIVE brain. Our findings point to a disruption of microglia-neuronal interactions in the HIV infected brain and an interrelation of retroviral integration and expression with interferon-associated remodeling of the microglial 3D genome during progression to HIVE. Search terms: HIV, postmortem brain, microglia, substance abuse, viral integration, chromosomal conformation capture, Hi-C, chromatin, single nuclei RNA-seq, Irf8, NeuN, HIV encephalitis, neuron-microglia interaction, cytokine

Study phs003080

Effect of ETS2 modulation on chromatin accessibility and enhancer activity in inflammatory (TPP) macrophages

This dataset consists of H3K27ac ChIP-seq and ATAC-seq in primary human macrophages following modulation of ETS2 expression. There are three datasets: H3K27ac ChIP-seq in ETS2-edited and unedited inflammatory (TPP) macrophages, H3K27ac ChIP-seq in ETS2-overexpressing and control M0 macrophages, and ATAC-seq in ETS2-edited and unedited inflammatory macrophages. For editing experiments, monocytes were transfected with CRISPR-Cas9 ribonucleoproteins containing a gRNA that was validated to edit ETS2, or a non-targeting control gRNA. Inflammatory macrophage differentiation was performed using conditions that model chronic inflammation (TPP): 3 days GM-CSF (50ng/mL) followed by 3 days GM-CSF, TNFa (50ng/mL), PGE2 (1mg/mL), and Pam3CSK4 (1mg/mL). For overexpression experiments, controlled overexpression of ETS2 mRNA or control mRNA (an equivalent amount of mRNA encoding the reverse complement of ETS2 – thereby controlling for the quantity, length and purine/pyrimidine composition of the transfected RNA but with a transcript that would not be translated) was induced in resting, non-activated (M0) macrophages by transfecting 500ng mRNA. To minimise non-specific activation due to the transfected RNA, in vitro transcription was performed using co-transcriptional capping (to minimise uncapped products), and incorporating modified, minimally immunogenic nucleotides (replacing uridine with N1-methyl-pseudouridine and cytidine with methylcytidine). After 18 hours, transfected cells were were activated with low dose LPS and harvested 6 hours later. For H3K27ac ChIP-seq, cells were harvested, cross-linked, quenched, lysed, and sheared. Immunoprecipitation of histone-DNA complexes was performed overnight at 4C with rotation using an anti-H3K27ac antibody and the SimpleChIP Plus Sonication ChIP kit (Cell Signaling Technology). Following reverse cross-linking, 50ng of immunoprecipitated DNA or input DNA were used to prepare sequencing libraries using the iDeal Library Preparation kit (Diagenode), according to manufacturer instructions. 10 PCR cycles were used for the amplification step and size selection was not performed. The quality and molarity of all libraries was assessed using a BioAnalyzer 2100 (Agilent) and the libraries were sequenced on an Illumina NovaSeq 6000. Raw data are provided as raw and aligned sequencing reads from H3K27ac-bound DNA and the input chromatin. ATAC-seq in TPP macrophages was performed using the Omni-ATAC protocol (Corces et al 2017) with the following modifications: cell number was increased to 75,000 cells, cell lysis time was increased to 5 minutes; volume of Tn5 transposase in the transposition mixture was doubled; duration of the transposition step was extended to 40 minutes. Sequencing libraries were prepared by PCR with the NEBNext Ultra II Q5Master Mix (NEB) as described previously (Corces et al 2017). Amplified libraries were cleaned using AMPure XP beads (Beckman Coulter). Sequencing reads were trimmed using TrimGalore! (Phred score 24), filtered to remove reads < 20 bp and aligned to the human genome (hg19) using Burrows-Wheeler Aligner (BWA). Aligned reads were converted to BAM files, sorted and indexed using SAMtools (Li et al, 2009). PCR duplicates and unmapped reads were removed, and the resulting BAM files were re-sorted and indexed. Raw data are provided as 100 bp paired-end Illumina reads from n = 3 donors for each experiment. BAM files are also provided.

Study EGAS50000000109

DNA and RNA sequencing of single human haploid germ cells

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

Study phs002279

Multivalent State Transitions Shape the Intratumoral Composition of Small Cell Lung Carcinoma

Small cell lung carcinoma (SCLC) is an aggressive, tobacco-associated tumor characterized by rapid growth, early metastases, and initial response followed by almost invariable resistance to chemotherapy. Studies to date have not resolved the extent that diverse transcriptional programs drive SCLC and contribute to its lethality. We combined patient-derived xenograft (PDX) (n = 64) resources with multi-omic profiling, single-cell fluorescence tracking of ex vivo tumor cells, and mathematical and statistical models to study the topology of the SCLC transcriptional state space and its plasticity.

Study phs003102

Blood DNA Methylation in Post-Acute Sequelae of COVID-19 (PASC): a Prospective Cohort Study

Background: DNA methylation integrates environmental signals with transcriptional programs. COVID-19 infection induces changes in the host methylome. While post-acute sequelae of COVID-19 (PASC) is a long-term complication of acute illness, its association with DNA methylation is unknown. No universal blood marker of PASC, superseding single organ dysfunctions, has yet been identified. Methods: In this single center prospective cohort study, PASC, post-COVID without PASC, and healthy participants were enrolled to investigate their symptoms association with peripheral blood DNA methylation data generated with state-of-the-art whole genome sequencing (WGS). PASC-induced quality-of-life deterioration was scored with a validated instrument, SF-36. Analyses were conducted to identify potential functional roles of differentially methylated loci, and machine learning algorithms were used to resolve PASC severity. Findings: 103 patients with PASC (22.3% male, 77.7% female), 15 patients with previous COVID-19 infection but no PASC (40.0% male, 60.0% female), and 27 healthy volunteers (48.1% male, 51.9% female) were enrolled. Whole genome methylation sequencing revealed 39 differentially methylated regions (DMRs) specific to PASC, each harboring an average of 15 consecutive positions, that differentiate patients with PASC from the two control groups. Motif analyses of PASC-regulated DMRs identify binding domains for transcription factors regulating circadian rhythm and others. Some DMRs annotated to protein coding genes were associated with changes of RNA expression. Machine learning support vector algorithm and random forest hierarchical clustering reveal 28 unique differentially methylated positions (DMPs) in the genome discriminating patients with better and worse quality of life.

Study phs003658

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

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

Study JGAS000146

Biallelic tumor suppressor loss and DNA repair defects in de novo small cell prostate cancer

Small cell prostatic carcinoma (SCPC) is an aggressive pathology that is managed similar to small cell lung cancer. SCPC can evolve from prostatic adenocarcinoma in response to androgen deprivation therapy but in rare cases is present at initial cancer diagnosis. The molecular etiology of de novo SCPC is incompletely understood due to the scarcity of tumor tissue and the brief life expectancy of patients. Through a retrospective search of our regional oncology pharmacy database, we identified 18 patients diagnosed with de novo SCPC between 2004 and 2017. Ten patients had pure SCPC pathology and the remainder had some admixed adenocarcinoma foci, but all were treated with first line platinum-based chemotherapy. Median overall survival was 28 months. We performed targeted DNA sequencing, whole exome sequencing and mRNA profiling on formalin-fixed paraffin embedded archival tumor tissue. We observed frequent biallelic deletion and/or mutation of tumor suppressors TP53, RB1 and PTEN, similar to treatment related SCPC. Indeed, at the RNA level pure de novo SCPC closely resembled treatment-related SCPC. However, five patients had biallelic loss of DNA repair genes including BRCA1, BRCA2, ATM and MSH2/6: potentially underlying the high genomic instability of this rare disease variant. Two patients with pure de novo SCPC harbored ETS gene rearrangements to androgen-driven promoters, consistent with evolution of de novo SCPC from an androgen-driven ancestor. Overall, our results reveal a highly aggressive molecular landscape underlies this unusual pathologic variant, and propose opportunities for targeted therapy strategies in a disease with few treatment options.

Study EGAS00001003007

Pharmacogenomic Analysis Reveals New Therapeutic Options for Pleural Mesothelioma

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

Study EGAS00001007866

Dataset for Sarcoma RNA sequencing linked from study EGAS00001004813

Sequencing data of 410 sarcoma tumor runs, which were uploaded to EGAS00001004813. The sequencing was always paired.

Dataset EGAD00001010276

The Genomic and Transcriptomic Landscape of a HeLa Cell Line

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

Study phs000643

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

Analysis of CD20 loss in patients treated with Mosunetuzumab

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

Study EGAS50000000151

T-cell receptor sequencing of tumor-infiltrating lymphocytes (TIL) in breast cancer

The spatial distribution of tumor-infiltrating lymphocytes (TIL) predicts breast cancer outcome and response to systemic therapy, highlighting the importance of an intact tissue structure for characterizing patients’ tumors. Here, we present ST-FFPE, a spatial transcriptomics method for the analysis of formalin-fixed paraffin-embedded samples, which opens the possibility of interrogating widely available archival tissue. The method involves extraction, exome capture and sequencing of RNA from different tumor compartments microdissected by laser-capture, and can be used to study the cellular composition of various tumor environments. Focusing on triple-negative breast cancer (TNBC), we characterized T cells, B cells, dendritic cells, fibroblasts and endothelial cells in both stromal and intra-epithelial compartments. We found a highly variable spatial distribution of immune cell subsets among tumors. This analysis revealed that the immune repertoires of intra-epithelial T and B cells were consistently less diverse than those of stromal T and B cells. T-cell receptor (TCR) sequencing confirmed a reduced diversity and revealed higher clonality of intra-epithelial T cells relative to the corresponding stromal T cells. Analysis of the top ten dominant clonotypes in the two compartments showed a majority of shared but also some unique clonotypes both in stromal and intra-epithelial T cells. Hyperexpanded clonotypes were more abundant among intra-epithelial than stromal T cells. These findings validate the ST-FFPE method and suggest an accumulation of antigen-specific T cells within tumor core. Because ST-FFPE is readily applicable for analysis of previously collected tissue samples, it could be useful for rapid assessment of intratumoral cellular heterogeneity in multiple disease and treatment settings.

Study EGAS00001007125

Multiomic Landscape of Multiple Myeloma Precursor and Relapsed Disease

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

Study phs003892

Single-cell RNA-seq of peripheral blood mononuclear cells in classic Hodgkin lymphoma

Classic Hodgkin lymphoma (cHL) is a largely MHC class I-negative tumor with recurrent 9p24.1/ PD-L1/ PD-L2 copy gains and the highest reported response rates to PD-1 blockade. We utilized scRNA sequencing to characterize the peripheral immune response to PD-1 blockade and more broadly define non-CD8+ dependent mechanisms of immune evasion in cHL. Peripheral blood mononuclear cells were obtained from 20 patients with relapsed/refractory (R/R) cHL treated with PD-1 blockade (nivolumab) on the CheckMate 205 clinical trial (paired samples from cycle 1 day 1 [C1D1] and cycle 4 day 1 [C4D1]), 11 patients with newly diagnosed, previously untreated cHL and 13 healthy donors.

Study EGAS00001007569

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Eccrine poroma (EP) and porocarcinoma (EPC) are rare benign and malignant adnexal neoplasms of the terminal sweat gland duct, respectively. Both can arise de novo, however, EPCs can also arise from a pre-existing EP. To-date, genetic investigation of these tumors has involved studies with small sample sizes and/or limited analyses. To comprehensively compare the driver events and mutational landscape of these tumors, we performed a retrospective multi-institutional whole-exome sequencing and RNA sequencing study on the largest cohort of EPs and EPCs to-date (n=90). We uncovered novel events and delineated different pathways of tumorigenesis underlying these tumors, with EPs driven largely by fusion genes, and EPCs driven largely by somatic mutations, with rare YAP1 and frequent PAK gene novel fusions.

Dataset EGAD00001015376

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Eccrine poroma (EP) and porocarcinoma (EPC) are rare benign and malignant adnexal neoplasms of the terminal sweat gland duct, respectively. Both can arise de novo, however, EPCs can also arise from a pre-existing EP. To-date, genetic investigation of these tumors has involved studies with small sample sizes and/or limited analyses. To comprehensively compare the driver events and mutational landscape of these tumors, we performed a retrospective multi-institutional whole-exome sequencing and RNA sequencing study on the largest cohort of EPs and EPCs to-date (n=90). We uncovered novel events and delineated different pathways of tumorigenesis underlying these tumors, with EPs driven largely by fusion genes, and EPCs driven largely by somatic mutations, with rare YAP1 and frequent PAK gene novel fusions.

Dataset EGAD00001015377

Multiple Myeloma treated with anti-BCMA T-cell-redirecting therapy (1st Batch)

Bone marrow samples were collected from patients treated with either anti-BCMA CART and/or anti-BCMA TCE. All samples were sorted using the StemCell EasySep Human CD138 positive selection kit II (StemCell, Cambridge, MA). The positive CD138+ were submitted for WGS, using peripheral blood mononuclear cells as normal matched. WGS was carried out by the Molecular Genomics Core Facility at Moffitt Cancer Center contingent on the availability of high-quality pre-treatment samples. Germline DNA was obtained from peripheral blood mononuclear cells, while tumor DNA was extracted from CD138-enriched bone marrow aspirate specimens using magnetic bead-based selection. Genomic DNA was isolated using the AllPrep DNA/RNA/miRNA Universal Kit (Qiagen). For library preparation, 100 ng of genomic DNA was fragmented using the Covaris ML230 focused ultrasonicator (Covaris, Inc.), followed by library construction with the TruSeq DNA Nano Library Prep Kit (Illumina) in accordance with manufacturer instructions. Library quality and fragment size distribution were assessed via Agilent TapeStation, and quantification was performed using the Kapa Library Quantification Kit (Roche). Sequencing was conducted on the Illumina NovaSeq 6000 system using S4-300 reagent kits with a 2×150 bp paired-end configuration. Tumor samples were sequenced to a median coverage of 60×, and matched germline controls to 30×. Reference Genome: GRCh38.

Dataset EGAD50000002632

BRAIN Cell EncyclOpeDia of Transcribed Elements (BRAINcode)

BRAINCODE: How Does the Human Genome Function in Specific Brain Neurons? The human brain comprises about 86 billion neurons whose function is central to human biology. How does the human genome program high performing neurons and neural networks in response to experience? What subprograms does the genome express in physiologically and morphologically distinct brain cells? The goal of the BRAIN Cell encyclOpeDia of transcribed Elements Consortium (BRAINcode) is to provide a map of gene expression - both protein-coding and non-coding - in specific cell types, not in culture, but in situ in brains of people. Going beyond traditional mRNA sequencing, polyadenylated and non-polyadenylated transcripts were ultra deeply sequenced using ribo-depleted RNA from neurons laser-captured from human post-mortem brains. Three prototypical neuron types, dopamine neurons, pyramidal neurons, and Betz cells, were prioritized because of their key biologic roles and differential vulnerability to important neurodegenerative diseases such as Parkinson's or Alzheimer's disease. Genetic variation between individuals was examined for correlation with differences in transcribed sequences to identify regions of the genome that influence whether, how, and how much a transcript is expressed in specific cell types in human brains. Our results indicate a vast universe of annotated and novel non-coding RNAs expressed in brain cells and suggest a more diverse and much more complex transcriptional architecture than previously imagined.

Study phs001556

Integrated Metabolic Profiling and Gene Expression Analysis Reveals Therapeutic Modalities in Breast Cancer

In this study, we performed metabolite profiling and RNA-sequencing on human breast cancer tumors and normal breast tissue samples, including triple negative breast cancer (TNBC), estrogen receptor positive breast cancer (ER+) and paired normal samples, and patient-derived xenografts (PDX) samples. Here, we report two distinctive groups defined by metabolites; a TNBC-HIGH group that shows high levels of pyrimidine pathway metabolites and biosynthetic enzymes, and an ER-HIGH group that shows high levels of fatty acid and arginine biosynthesis intermediates. We identify different metabolic enrichment profiles between cell lines grown in vitro vs. in vivo; cell lines grown in vivo recapitulate patient tumors metabolic profiles. Further, we identify a subset of genes that strongly correlates with the TNBC-HIGH metabolic profile using integrated metabolic and gene expression profiling, which strongly predicts patient prognosis when tested on larger human datasets. As a proof-of-principle, when we target TNBC-HIGH metabolic dysregulation with a pyrimidine biosynthesis inhibitor (Brequinar), and/or a glutaminase inhibitor (CB-839), we observe therapeutic efficacy and decreased tumor growth in representative TNBC cell lines, murine p53-null mammary tumors, and in vivo patient-derived xenografts (PDXs) upon multi-agent drug treatment. This study highlights potential new therapeutic opportunities in breast cancers guided by a genomic biomarker, which could prove impactful for breast cancers that rapidly proliferate.

Study phs002396

Genetic and Hormonal Contributions to Gene Expression in Immune Cells

The origins of sex differences in human disease are elusive, in part because of difficulties in separating the effects of sex hormones and sex chromosomes. To isolate the effects of testosterone, estradiol, and sex chromosome complement (XX or XY) in humans, we implemented a four-state, cross-sectional study design of trans- and cisgender individuals, the former receiving exogenous testosterone or estradiol for at least a year as part of their gender-affirming medical care. This study design allowed us to compare four combinations of sex hormones and sex chromosomes in humans: XX individuals treated with exogenous testosterone (n=21), XY individuals treated with exogenous estradiol (n=13), untreated XX individuals (n=20), and untreated XY individuals (n=15). We chose to focus on the human immune system because it displays an impressively broad range of sex-differential biology. We collected blood samples from each individual and isolated PBMCs the same day. We performed 10X Genomics 3' single-cell RNA-sequencing (scRNA-seq) on PBMCs from the 69 individuals. We generated 150-bp paired end reads from a total of 358,426 cells (with an average of 48,007 reads/cell) and identified 18 distinct immune cell types. We found that testosterone, estradiol, and sex chromosome complement have unique, non-overlapping effects on PBMC cell abundances and gene expression.

Study phs003860

Investigating Host Genetic Risk Factors for Tuberculosis in Highly Endemic South African Populations

It is apparent from our recent population genetic and admixture mapping work 1–3 that a substantial part of the ancestral input into the South African Coloured (SAC) population is from the San or Khoe groups, which are not well represented in publicly available genetic databases. It is a reasonable assumption that the SAC population in the Western Cape may have derived more genetic input from the present KhoeSan population in the Northern Cape than the San residing in Namibia. It was concluded by our group that it is indeed the southern African KhoeSan group, the ≠Khomani, that best represent the KhoeSan contribution seen in the SAC 4. More recently, it has been shown that there are ancestry related increases in TB susceptibility, especially with increased Bantu-speaking African and KhoeSan ancestry 5. We aim to elucidate the epidemiological and human-host genetic risk factors for TB and the immunological pathways modulating TB infection. The proposed sub-study will enroll people evaluated for TB at Northern Cape community health clinics and their contacts from their households and/or community. We will conduct a demographic interview, collect saliva, which will be used for genetic analysis, and blood to determine latent TB infection status and capture immunological responses to mycobacterial infections at the cellular level and RNA sequencing.

Study EGAS00001007850

INSPIRE - Investigator-initiated Phase 2 Study of Pembrolizumab Immunological Response in Metastatic Solid Tumors

Whole exome and RNA sequencing source data for: Clouthier, D.L., Lien, S.C., Yang, S.Y.C. et al. An interim report on the investigator-initiated phase 2 study of pembrolizumab immunological response evaluation (INSPIRE). j. immunotherapy cancer 7, 72 (2019). https://doi.org/10.1186/s40425-019-0541-0 Bratman S*, Yang SYC*, Iafolla MAJ, Liu Z, Hansen AR, Bedard P, Lheureux S, Spreafico A, Razak AA, Shchegrova S, Louie M, Billings P, Zimmerman B, Sethi H, Aleshi A, Torti D, Marsh K, Eagles J, Cirlan I, Hanna Y, Clouthier D, Lien SC, Ohashi PS, Xu W, Siu LL, Pugh TJ. Personalized circulating tumor DNA analysis as a predictive biomarker in solid tumor patients treated with pembrolizumab. Nature Cancer. 2020:1:873-881. https://doi.org/10.1038/s43018- 020-0096-5. Yang SYC, Lien SC, Wang BX, Clouthier DL, Hanna Y, Cirlan I, Zhu K, Bruce JP, El Ghamrasni S, Iafallo MAJ, Oliva M, Hansen AR, Spreafico A, Bedard PL, Lheureux S, Razak A, Speers V, Berman HK, Aleshin A, Haibe-Kains B, Brooks DG, McGaha TL, Butler MO, Bratman SV, Ohashi PS, Siu LL, and Pugh TJ. Pan-cancer analysis of longitudinal metastatic tumors reveals genomic alterations and immune landscape dynamics associated with pembrolizumab sensitivity. Nature Comm 2021:12, 5137. https://doi.org/10.1038/s41467-021-25432-7

Study EGAS00001003280

A prognostic human brain network for diffuse midline glioma

Diffuse midline glioma (DMG) are near-universally lethal tumours of the paediatric central nervous system. In animal models, DMG form brain-wide, integrated networks through neuron-to-glioma synapses and glioma-to-glioma gap junctional coupling. This extensive connectivity robustly promotes the growth and invasion of DMG and other glial malignancies through paracrine mechanisms and direct neuron-to-glioma synapses. The organisation and clinical implications of these connections in the living human brain, however, remain to be elucidated. Here, we develop tumour network mapping to compute the brain-wide connectivity profile of DMG: defining a conserved brain network across pontine and thalamic DMG associated with patient short-term survival (DMG network). Tumour functional connectivity with the DMG network was independently predictive of patient overall survival across two external validation cohorts. Tumour growth mapped to DMG network-specific trajectories and peak, in-network neurometabolic changes across development spatiotemporally aligned with the peak age incidence of DMG. Single-nucleus RNA-sequencing confirmed diverse synaptic gene enrichment in high-connectivity DMG. Strikingly, incidental surgical resection of high-connectivity thalamic DMG tissue conferred a significant survival advantage. Collectively, these data define a conserved and prognostically important brain network in children with DMG, consistent with the hypothesis that DMG exploit otherwise healthy brain circuits to promote tumour growth.

Study EGAS50000001752

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

Single cell transcriptomics of hESC-derived midbrain dopaminergic neurons generated by a new human development-based protocol

Recent advances in stem cell technologies provide new opportunities for regenerative medicine and Parkinson’s disease modelling. However, the precise composition and quality of human embryonic stem cell (hESC)-derived preparations at the single cell level, compared to endogenous brain standards, remains to be determined. In addition, multiple developmental pathways important for midbrain dopaminergic (mDA) neuron development remain to be implemented in hESCs. Here we use mouse developmental knowledge and human single cell RNA-sequencing (scRNA-seq) data to generate a new development-based protocol to differentiate hESCs into mDA neurons. We found that Laminin 511, dual WNT activation with CHIR99021 and WNT5A, together with FGF8b, improved different aspects of midbrain patterning. Moreover, mDA neurogenesis and differentiation were enhanced by activation of liver X receptors and inhibition of fibroblast growth factor signaling. Analysis of hESC-derived midbrain cultures by scRNA-seq revealed multiple hESC-derived clusters, which resembled those in the endogenous human ventral midbrain. Moreover, high quality mDA neurons were detected and found to progressively acquire functionality in vitro. Thus, we hereby provide a hESC differentiation protocol that recapitulates several aspects of endogenous mDA neuron development, and a strategy to determine cell composition and quality of hESC-derived preparations that can be useful for disease modeling and cell therapy.

Study EGAS00001006313

BMP4 and temozolomide synergize in the majority of patient derived glioblastoma cultures

The poor prognosis of patients with glioblastoma is partly due to resistance to standard treatment, which combines chemoradiation using temozolomide (TMZ). The concept of cancer stem cells provided new insight into tumor resistance and management, and isolated Glioblastoma stem-like cells (GSCs) which contribute to this therapy resistance. Bone morphogenetic protein 4 (BMP4) is a differentiation factor that stimulates astroglial differentiation of such GSCs, thereby reducing their self-renewal. Exposing GSCs to BMP4 could therefore sensitize them to chemotherapy with TMZ, but previous studies with various treatment paradigms have produced contradictory results in this matter. Therefore, we re-defined the optimal treatment paradigm for the combination of TMZ and BMP4 and subsequently re-assessed the inter-tumoral variability in response to co-treatment with TMZ+BMP4 in vitro. First, the simultaneous treatment was found more effective than sequential therapy. Second, after application of our optimized treatment protocol, 70% of a total of 20 patient-derived glioblastoma cultures showed synergy when the treatments were combined. The synergistic effect was a result of increased cellular apoptosis and not of decreased proliferation. Comparative RNA-sequencing suggested that combination treatment results in decreased MAPK signaling. We conclude that addition of BMP4 to standard therapy with TMZ could enlarge the group of patients with good therapeutic response.

Study EGAS00001007095

Resistance to Checkpoint Blockade Therapy Through Inactivation of Antigen Presentation

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

Study phs001427

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

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

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

Study EGAS00001002201

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

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

Study EGAS00001004564

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

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

Study EGAS00001004567

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

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

Study EGAS00001004555

Addiction of Primary Cutaneous Gamma-Delta T-cell Lymphomas to JAK-STAT Signaling

Primary cutaneous γδ T-cell lymphoma (PCGDTL) is a mature T-cell lymphoma (TCL) of the skin with a median survival of 31 months and no effective therapies. PCGDTLs frequently possess mutations in the JAK/STAT pathway, including in JAK1, JAK3, STAT3, STAT5B, and SOCS1. We hypothesized that single-agent JAK inhibition may be effective in PCGDTLs, particularly those with targetable JAK/STAT pathway mutations. Herein, we demonstrate the near-complete response of two PCGDTL cases to single-agent JAK inhibitors, ruxolitinib or cerdulatinib. However, both patients eventually relapsed with acquired mutations in JAK1, JAK3, or STAT5B.

Study phs003837

Paired DNA and RNA sequencing uncovers common and rare genomic variants regulating gene expression in the human retina

In this study, we uncovered common and rare genomic variants shaping retinal expression profiles. This includes 1,483,595 significant cis-expression quantitative trait loci (eQTLs) impacting 9,959 and 3,699 genes in NSR and RPE, respectively, with associated genomic variants enriched to cis-candidate regulatory elements and notable shared eGenes between NSR and RPE. We also detected 1051 expression outliers and prioritised 299 rare non-coding single-nucleotide, structural variants or copy number variants as plausible drivers for 28% of outlier events. This study increases understanding of gene expression regulation in the human retina.

Study EGAS50000001443

RNA sequencing data of 142 samples from 142 patients with HER2+ breast cancer treated with letrozole or tamoxifen (SOLTI PAMELA trial)

PAMELA (NCT01973660) is a non-randomized, open-label, multicentric, prospective translational research study in stage I-IIIA HER2+ breast cancer designed to evaluate the ability of the PAM50 intrinsic subtypes to predict pCR in the breast (pCRB) following 18 weeks of neoadjuvant lapatinib and trastuzumab (+/- hormonal therapy if hormone receptor [HR] positive). Patients with HR+ disease received letrozole (if postmenopausal) or tamoxifen (if pre-menopausal). The primary objective was to compare the pCRB rates of the HER2-E versus the non-HER2-E subtypes in the intent-to-treat population.

Study EGAS00001006410

RNA sequencing of circulating human immune cells before and after interleukin-2 immunotherapy in systemic lupus erythematosus patients

Within the investigator-initiated clinical trial Charact-IL-2 (ClinicalTrials.gov Identifier: NCT03312335), we treated 12 female systemic lupus erythematosus patients with low-dose interleukin-2. Patients received 4 cycles of 5 daily injections of 1.5 million international units interleukin-2 over 9 weeks. This study aimed at investigating interleukin-2-induced transcriptional changes in circulating human immune cells, and to define changes in T cell receptor clonality before and after interleukin-2 immunotherapy. For this purpose, blood was collected from patients before and after each treatment cycle.

Study EGAS50000000458

TSO500 STJAN33, BRAF mutated CUP

The data set contains FASTQ files (filetype) of a NEXTSEQ550DX run (instrument). De FASTQ files are from DNA and RNA sample. The Library prep is a NGS TSO500 library prep (Illumina) (technology). A Study to Examine the Clinical Value of Comprehensive Genomic Profiling Performed by Belgian NGS Laboratories: a Belgian Precision Study of the BSMO in Collaboration With the Cancer Centre (BALLETT) This 2-year study involves the consortium of 9 cooperating Belgian NGS laboratories and will enroll 936 metastatic or locally advanced cancer patients coming from 13 different Belgian hospitals and cancer centers. Upon inclusion, all cancer patients will be offered 'comprehensive genomic profiling' (CGP) using Illumina's TSO500 NGS panel. This targeted NGS panel of 523 genes allows for the detection of single nucleotide variants, small indels, copy number variations and fusions, as well as for the determination of the 'tumor mutational burden' (TMB) and the 'microsatellite-instability' status (MSI). Both the wet lab execution of the CGP as well as the biological and clinical classification of the variants will be performed in a fully standardized way among the 9 participating Belgian local NGS laboratories.

Dataset EGAD50000000689

RNA sequencing data of pediatric KMT2A-rearranged acute lymphoblastic leukemia

Pediatric KMT2A-rearranged leukemia is an aggressive subtype of leukemia commonly seen in infants (less tha 1 year old) and young children. It involves rearrangements of the KMT2A gene (formerly MLL) on chromosome 11q23, resulting in fusion genes that drive uncontrolled cell growth. This subtype is most often found in acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML) and is associated with a poor prognosis. In infants, the 5-year event-free and overall survival rates are around 35–40%. Despite intensive chemotherapy, outcomes remain challenging, prompting ongoing research into targeted therapies and novel treatment strategies.

Study EGAS50000001799

RNA sequencing data of pediatric T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignancy accounting for approximately 15–25% of pediatric and adult ALL cases. Despite significant therapeutic improvements over the last decade, the prognosis of T-ALL remains less favorable than that of B-ALL, with 5-year EFS and OS of around 76–82% and 80–88%, respectively, depending on the treatment protocol used. T-ALL patients often suffer from failure of induction therapy or early relapse, which is associated with a dismal prognosis. Approximately 10–15% of T-ALL patients experience relapse, and survival after relapse is poor.

Study EGAS50000001800

Alignment BAM files and gene count files of the Illumina Sequencing Data (10 tumor samples)

Reads were processed with the RNA-seq workflow 1.3.0 developed by the DKFZ Omics IT and Data Management Core Facility (https://github.com/DKFZ-ODCF/RNAseqWorkflow). First, FASTQ reads were aligned via two-pass alignment using STAR 2.5.3a. The STAR index was generated from the 1000 Genomes assembly and GENCODE Version 19 gene models with a sjdbOverhang of 200. Duplicate marking of the resultant main alignment file was done with sambamba 0.6.5. Gene-specific read counting was performed using featureCounts (from Subread 1.5.1) over exon features based on GENCODE Version 19 gene models. Both reads of a paired fragment were used for counting, and the quality threshold was set to 255, indicating that STAR found a unique alignment. Strand-specific counting was also used. For RPKM and TPM calculations, all genes on chromosomes X and Y, the mitochondrial genome, as well as rRNA and tRNA genes were omitted as they are likely to introduce library size estimation biases.

Dataset EGAD00001008969

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

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

Dataset EGAD00001009813

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

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

Study EGAS00001005667

Exome sequencing

Human telomere biology disorders (TBD), a heterogeneous group of disorders characterized by telomere attrition and premature aging phenotypes, are caused by inherited loss-of-function mutations in telomere-associated genes. Here, we identify three germline heterozygous missense variants in RPA1 gene in four unrelated probands presenting with short telomeres and varying clinical features of TBD including bone marrow failure, myelodysplastic syndrome, T- and B-cell lymphopenia, pulmonary fibrosis, and skin manifestations. All variants cluster to DNA binding domain A of RPA1 protein. RPA1 is an essential single-strand DNA-binding protein essential for DNA replication and repair. RPA1 has been implicated in telomere maintenance but the mechanism remains elusive. Thus far, no genetic disorder had been linked to RPA1. Using FRET assays, we showed that RPA1E240K and RPA1V227A exhibit increased binding to single-strand and telomeric DNA, implying a gain in DNA-binding function while RPA1T270A has binding properties similar to wild type. To study the mutational effect, RPA1E240K, with the highest binding capacity, was knocked-in iPSCs using CRISPR/Cas9, which showed severe telomere shortening and impaired hematopoietic differentiation. Furthermore, in the RPA1E240K patient, stable blood counts over two decades due coincided with clonal rescue hematopoiesis with an RNA-degrading truncating RPA1 mutation and a uniparental isodisomy 17p with loss of RPA1E240K allele. Using single cell DNA sequencing, the two somatic genetic rescue events were proven to be independently acquired in hematopoietic stem cells. In summary, we describe the first human disease caused by germline RPA1 variants in individuals with TBD.

Study EGAS00001005761

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

Genomic profiling of paediatric high grade gliomas from the HERBY clinical trial

The HERBY trial was a phase-II open-label, randomised, multicentre trial evaluating bevacizumab in patients with newly-diagnosed non-brainstem high grade glioma (HGG) between the ages of 3-18yrs. 121 patients were randomised with 1yr event-free survival as the primary end-point. Confirmation of HGG diagnosis by reference pathology was mandatory before randomisation, followed by review with five independent expert neuropathologists. We collected specimens from 89 patients consenting to translational research, and performed Sanger sequencing for H3F3A, Illumina 450K BeadArray methylation profiling, and whole exome sequencing and RNA sequencing where possible. 7/89 patients (8%) harboured H3F3A G34R/V mutations, whilst 24/89 (27%) harboured H3F3A K27M, the latter reflecting a high proportion of the novel entity recognised in the 2016 WHO classification of diffuse midline glioma with H3K27M mutation. Both histone mutations conferred a significantly shorter progression-free (p=0.0104) and overall survival (p=0.00159). 450K methylation subtyping additionally identified a number of subgroups in histone wild-type cases. These included infrequent (4/86, 4.6%) cases with IDH1 mutation in older children, further categorised into astrocytic ATRX-mutant (n=3) or 1p19q co-deleted with TERT promoter mutation (n=1). 9/74, (12%) of cases were classified as biologically resembling pleomorphic xanthoastrocytoma (PXA) by methylation profiling, with 5/9 harbouring BRAF V600E mutations and epithelioid histology, and 3/9 NF1 mutation and giant cell features. Three cases had methylation profiles more closely resembling low grade gliomas, though were morphologically high grade, and harboured MAPK dysregulation in two cases, with the third part of an additional infant (<3yrs) cohort. Finally, four cases had a mutational burden several orders of magnitude higher than the rest, with (2197-5332) somatic coding variants/sample reflecting a hypermutator phenotype. These data provide an important insight into the wide biological diversity of ‘HGG’ found in the paediatric age group and allow for a profound refinement of clinical trial interpretation.

Study EGAS00001002328

The G2 gene expression signature and MYC overexpression are independent poor prognostic factors in childhood high-grade osteosarcoma

Purpose In high-grade osteosarcoma, prognostic factors at diagnosis are insufficient for stratifying patients into relevant subgroups. Recently, a transcriptomic study developed the G1/G2 gene expression signature, in which the G2 signature was associated with unfavorable survival. An orthogonal study identified MYC amplification as an unfavorable prognostic factor using targeted next-generation sequencing. The purpose of this study was to validate the independent prognostic value and to investigate the combined prognostic value of the G1/G2 signature with, MYC amplification and/or MYC expression for survival prediction. Material and methods This study included pediatric and adolescent patients with high-grade osteosarcoma. RNA sequencing was performed in 48 patients. Whole-exome sequencing was performed in 40 patients. Gene expression signature scores, MYC amplification (defined as > 7 copies) and MYC expression levels were calculated. Multivariable Cox proportional hazards analysis was performed for event-free survival (EFS) (primary endpoint) and overall survival (OS) (secondary endpoint). Results In the full cohort, the 3-year event-free survival rate was 37%. In multivariable Cox regression analysis with metastatic disease stage (n=21, 44%) as covariate, the G2 signature and MYC expression were independently associated with worse outcomes in terms of EFS (Hazard ratio (HR) 3.32, 95% confidence interval (CI) 1.34-8.21 and HR 3.38, 95%-CI 1.71-6.66, respectively) and OS (HR 4.07, 95%-CI 1.19-13.9 and HR 2.88, 95%-CI 1.22-6.76, respectively). MYC amplification was not associated with EFS or OS in univariable analysis (HR 1.88, 95% -CI 0.74-4.77 and HR 0.79, 95%-CI 0.21-3.05, respectively). Conclusion The G2 gene expression signature and MYC expression were independently associated with unfavorable outcomes in a pediatric cohort of patients with high-grade osteosarcoma. The combined prognostic value warrants further prospective validation and could potentially serve as a stratification marker for future osteosarcoma treatment protocols.

Study EGAS00001008073

RNA-seq study of longitudinal blood cell samples from children at risk of type 1 diabetes

The appearance of type 1 diabetes (T1D)-associated autoantibodies is the first and only measurable parameter to predict progression toward T1D in genetically susceptible individuals. However, autoantibodies indicate an active autoimmune reaction, wherein the immune tolerance is already broken. Therefore, there is a clear and urgent need for new biomarkers that predict the onset of the autoimmune reaction preceding auto-antibody positivity or reflect progressive b-cell destruction. Here we report the mRNA sequencing–based analysis of 306 samples including fractionated samples of CD4+ and CD8+ T cells as well as CD4, CD8 cell fractions and unfractionated PBMC samples longitudinally collected from seven children who developed beta-cell autoimmunity (Case subjects) at a young age and matched control subjects.

Study EGAS00001004071

IFN-γ and TNF-α drive a CXCL10+ CCL2+ Macrophage Phenotype Expanded in Severe COVID-19 Lungs and Inflammatory Diseases with Tissue Inflammation

We explored how the human macrophages response to different inflammatory factors, focusing particularly on the effects of antiviral interferons (IFN-β and IFN-γ), pro-inflammatory cytokines such as TNF-α, and other mediators like IL-4. Co-cultured fibroblasts were a component in some conditions to generate factors produced by resident stroma. We used a single-cell antibody-based hashing strategy to multiplex samples from different stimulatory conditions in one sequencing run. This macrophage transcriptomic data reveals distinct macrophage activation states and polarizations shaped by different tissue-related inflammatory conditions at single-cell resolution. We further performed an unbiased integration between tissue-level macrophages and this stimulated human blood-derived macrophages, which pinpoint an IFN-γ and TNF-α synergistically driven inflammatory macrophage phenotype expanded in severe COVID-19 lungs and other inflamed disease tissues.

Study phs002780

Anaplastic oligodendroglioma exome and RNA sequencing data

Anaplastic oligodendrogliomas (AOs) are rare primary brain tumors which are generally incurable, with heterogeneous prognosis and few treatment targets identified. Most oligodendrogliomas have chromosome 1p/19q co-deletion and IDH mutation. We analyzed 51 AOs by whole-exome sequencing, identifying previously reported frequent somatic mutations in CIC and FUBP1. We also identified recurrent mutations in TCF12 and in an additional series of 83 AO. Overall 7.5% of AO are mutated for TCF12, which encodes an oligodendrocyte-related transcription factor. 80% of TCF12 mutations identified were in either the bHLH domain, which is important for TCF12 function as a transcription factor, or were frame shift mutations leading to TCF12 truncated for this domain. We show that these mutations compromise TCF12 transcriptional activity and are associated with a more aggressive tumor type. Our analysis provides further insights into the unique and shared pathways driving AO.

Study EGAS00001001209

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

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

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

Study phs002756

Relapse CHL study

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

Study EGAS00001008222

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

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

Study EGAS50000000973

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

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

Study EGAS00001003950

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

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

Study EGAS50000000135

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

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

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

Study EGAS00001003975

Identification of Somatic Changes in Tumors from Fanconi Anemia Patients

This study focused on the identification of genomic signatures present in tumors from Fanconi anemia patients. Fanconi anemia is a rare hereditary disease caused by defects in DNA interstrand crosslink repair. The primary patient phenotypes are bone marrow failure and cancer predisposition. We concentrated on squamous cell carcinomas and used whole genome Illumina sequencing (22 tumor samples and 12 normal samples) , whole exome sequencing (37 tumor samples and 33 normal samples), PacBio long-read whole-genome (nine tumor and four normal samples),10X linked-read whole-genome sequencing, 10x single-cell/nuclei 3’ RNA sequencing (three tumor samples), 10x Visium spatial transcriptomics (one tumor sample), EPIC 850K methylation array (six samples), and bulk RNAseq (six tumor samples) to identify somatic changes present in tumors. The study aimed to identify changes that drive tumorigenesis in Fanconi anemia, give insight into the early development and aggressive nature of these tumors, and suggest therapeutic avenues for these patients. Another goal was to compare the genomic changes identified in tumors from Fanconi anemia patients with sporadic head and neck cancers. This study demonstrated that the tumors from Fanconi anemia individuals are generally not driven by HPV. Instead, they acquire mutations in the p53 tumor suppressor. The main characteristic of the tumors is the presence of a high number of structural variants in the form of small deletions, unbalanced translocations, and fold-back inversions, which often form complex rearrangements. These rearrangements result in amplifications and deletions of a multitude of oncogenes and tumor suppressors, respectively. The study also shows that lack of the Fanconi anemia DNA repair pathway leads to an epithelial-to-mesenchymal transition in a mouse model and in human samples that may influence the behavior of tumors in Fanconi anemia. Finally, comparing genomic changes in Fanconi anemia and sporadic tumors of the head and neck, we propose that the structural variants present in sporadic cancer result in a functional deficiency of the Fanconi anemia DNA repair pathway.

Study phs002652

Exome Sequencing to Define the Landscape of Plasma Cells in Systemic Light chain Amyloidosis

Systemic light chain amyloidosis (AL) is characterized by the deposition of immunoglobulin light chains as amyloid fibrils in different organs, where they form toxic protein aggregates. The underlying disease is a plasma cell disorder, but limited whole exome data are available. We report the findings of an exome sequencing study in AL to define a plasma cell signature and compare this to monoclonal gammopathy of undefined significance (MGUS) and myeloma (MM). Twenty-four samples from unselected newly diagnosed untreated AL patients were analysed. CD138+ cells were isolated from bone marrow cells using MACSorting (Miltenyi Biotech, Bisley, UK). Cells were lysed in RLT+ buffer and DNA/RNA extracted using the AllPrep kit (Qiagen, Manchester, UK). Peripheral blood was isolated from patients, white blood cells purified by Ficoll-Pacque and DNA extracted using the QIAamp DNA mini kit (Qiagen). There were a median of 47 acquired nonsynonymous variants (range 4-253) per patient. Thirty percent of patients yielded a mutation in a myeloma driver gene. These mutations were either clonal or sub-clonal. There was evidence of MAPK activation with NRAS mutations, NF-κB activation and DNA repair pathway alterations. Exome data were used to determine the cytogenetic groups of AL samples and identified hyperdiploidy in 33% of cases and t(11;14) in 33% of cases. All translocations were generated via class switch recombination and breakpoints on chromosome 11 fell into the expected regions. A novel translocation, unlikely to be driving disease progression was identified involving RCC1. Finally, further analysis of the APCS gene encoding for the SAP protein did not show any biased genotype in the AL population in comparison to the general population. We conclude that exome sequencing identifies a genetic make-up of AL that is similar to other plasma cell disorders such as MM and MGUS. This not only includes copy number abnormalities and translocations but also nonsynonymous mutations.

Study EGAS00001001418

CEBP-Beta/IL-1-Beta/ TNF-alpha Feedback Loop Drives Drug Resistance to BCL2 and MDM2 Inhibitors in Monocytic Leukemia Cells

To uncover the mechanisms of MDM2 inhibitor resistance that are mediated by IL-1, we performed RNAseq on three primary AML patient cells treated with idasanutlin in the presence or absence of IL-1β. These three leukemia samples were chosen because they exhibit high blast counts and are protected by IL-1β when treated with idasanutlin, as determined by apoptosis assay. We observed that a majority of the major p53 pathway targets remained upregulated in the presence of IL-1β. We next performed differential expressed gene (DEG) analysis following IL-1β treatment with and without idasanutlin. IL-1β treatment upregulated 3295, 1353, and 680 genes with ≥1 Log2FC, in each of the three patient samples, and a significant majority of these genes (82%, 62.4%, and 47.6%) remained upregulated in the presence of idasanutlin. Pathway analysis showed that IL-1β upregulated genes were associated with NF-κB, Interleukin-1 family signaling pathways, immune and inflammatory associated genes and myeloid maturation markers. To identify downstream targets that contribute to the protective effect of IL-1β during idasanutlin treatment, we performed DEG analysis on IL-1β + idasanutlin treatment compared to idasanutlin treatment alone. Interestingly, IL-1β induces its own expression and other pro-survival cytokines including TNF, CSF2 (GMCSF), IL-6, CSF3 (GCSF), as well as their respective receptors and/or coactivators, including IL1RL1 and multiple TNF receptors. We verified that IL-1β and TNF-α induce the upregulation of their own membrane-bound forms, their corresponding receptors, CSF2RB, and monocytic markers, as determined by flow cytometry immunophenotyping. Furthermore, IL-1β upregulated a number of pro-inflammatory chemokines, including CXCL13, CXCL2, etc. These findings highlight a pleiotropic role of IL-1β. The RNA expression of these three samples in the presence and absence of IL-1β and idasanutlin treatment will be deposited in dbGaP.

Study phs003479

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

Gene Expression in an African American Schizophrenia Dataset

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

Study phs002842

RNA sequencing data of pediatric BCR::ABL1 acute lymphoblastic leukemia

The BCR::ABL1 fusion results from the t(9;22)(q34;11.2) translocation, also known as the Philadelphia chromosome. In children, BCR::ABL1-positive leukemia is rare, accounting for approximately 3% of cases. The BCR::ABL1 oncoprotein activates multiple signaling pathways that drive cell proliferation. Historically, this subtype was associated with an extremely poor prognosis; however, the incorporation of tyrosine kinase inhibitors into therapy has markedly improved outcomes, with overall survival now approaching 80% in pediatric patients. Relapse most commonly occurs through the emergence of clones harboring kinase domain mutations, yet 30–40% of relapses arise in the absence of such mutations, a phenomenon that remains poorly understood.

Study EGAS50000001798

RNA sequencing data of pediatric hypodiploid acute lymphoblastic leukemia

Hypodiploid acute lymphoblastic leukemia (ALL) is a rare subtype of B-cell ALL characterized by chromosome number < 45 , encompassing near-haploid (24–29 chromosomes) and low-hypodiploid (30–39 chromosomes) subgroups. These subgroups display distinct genetic landscapes, including Ras and receptor tyrosine kinase pathway alterations in near-haploid ALL, and frequent IKZF2, RB1, and germline TP53 alterations in low-hypodiploid ALL. Clinically, hypodiploid ALL is associated with dismal prognosis, with 5-year overall survival rates often between 50-70% in children and markedly worse in adults. Despite advances in risk stratification and MRD-guided therapy, allogeneic stem cell transplantation remains the current standard recommendation for affected patients.

Study EGAS50000001806

Epigenome maps of time-resolved monocyte to macrophage differentiation and innate immune memory

Innate immune memory is the phenomenon whereby innate immune cells such as monocytes or macrophages undergo functional reprogramming after exposure to microbial components such as LPS. We apply an integrated epigenomic approach to characterize the molecular events involved in LPS-induced tolerance in a time dependent manner. ChIP-seq, RNA-seq, WGBS and ATAC-seq data were generated. This analysis identified epigenetic programs in tolerance and trained macrophages, and the potential transcription factors involved. Experimental set-up Time-course in vitro culture of human monocytes. Two innate immune memory states can be induced in culture through an initial exposure of primary human monocytes to either LPS or BG for 24 hours, followed by removal of stimulus and differentiation to macrophages for an additional 5 days. Cells were collected at baseline (day 0), 1 hour, 4 hour, 24 hour and 6 days.

Dataset EGAD00001002693

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

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

Dataset EGAD00001015660

Influence of the Microbiome on Epigenetic Mechanisms in IBD

Live CD4 T cells were sorted from inflamed and non-inflamed tissue samples of IBD patients or from healthy and IBD blood samples. ATAC-Seq libraries were generated from live CD4 T cells sorted from i) inflamed and non-inflamed tissue samples, ii) healthy and IBD blood samples, or from iii) CD4 T cell subsets polarised from healthy blood samples. After isolating crude nuclei, live CD4+ T cells were treated with Tagment DNA buffer and Tagment DNA Enzyme (Nextera DNA Library Prep Kit, Illumina), and then the DNA was purified by MinElute PCR Purification Kit (Qiagen). Transposed DNA fragments were amplified using specific adapters followed by purification with MinElute PCR Purification Kit (Qiagen). Fragments from 240-360pb were selected in the PippinHT system (Sage Science). The quality of the library and its DNA concentration were assessed by Bioanalyzer instruments (Agilent Technologies) and ultimately submitted for sequencing using Illumina HiSeq 2500 sequencer, V4 chemistry. On the other hand, single cell RNA-Seq libraries were generated exclusively from inflamed and non-inflamed tissue samples of Crohn’s disease patients. Briefly, live CD4 T cells were captured and encapsulated before cDNA amplification using the 10X Genomics Chromium Platform. Samples were prepared as outlined by 10x genomics Single Cell 3’ Reagent Kits v2 user guide. Samples were sequenced on a HiSeq 2500 with the following run parameters: Read 1 – 26 cycles, read 2 – 98 cycles, index 1 – 8 cycles.

Dataset EGAD00001011066

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Eccrine poroma (EP) and porocarcinoma (EPC) are rare benign and malignant adnexal neoplasms of the terminal sweat gland duct, respectively. Both can arise de novo, however, EPCs can also arise from a pre-existing EP. To-date, genetic investigation of these tumors has involved studies with small sample sizes and/or limited analyses. To comprehensively compare the driver events and mutational landscape of these tumors, we performed a retrospective multi-institutional whole-exome sequencing and RNA sequencing study on the largest cohort of EPs and EPCs to-date (n=90). We uncovered novel events and delineated different pathways of tumorigenesis underlying these tumors, with EPs driven largely by fusion genes, and EPCs driven largely by somatic mutations, with rare YAP1 and frequent PAK gene novel fusions.

Dataset EGAD00001015378

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Eccrine poroma (EP) and porocarcinoma (EPC) are rare benign and malignant adnexal neoplasms of the terminal sweat gland duct, respectively. Both can arise de novo, however, EPCs can also arise from a pre-existing EP. To-date, genetic investigation of these tumors has involved studies with small sample sizes and/or limited analyses. To comprehensively compare the driver events and mutational landscape of these tumors, we performed a retrospective multi-institutional whole-exome sequencing and RNA sequencing study on the largest cohort of EPs and EPCs to-date (n=90). We uncovered novel events and delineated different pathways of tumorigenesis underlying these tumors, with EPs driven largely by fusion genes, and EPCs driven largely by somatic mutations, with rare YAP1 and frequent PAK gene novel fusions.

Dataset EGAD00001015379

Single-Cell Mitochondrial Mutation Lineage Tracing of Non-Dysplastic and Dysplastic Barrett's Esophagus

This study addresses a critical knowledge gap in understanding Barrett's esophagus (BE), a common metaplastic condition and precursor to esophageal adenocarcinoma, by identifying lineage connections and cell states underlying its development and malignant transformation. We collected biopsies from 13 Barrett's esophagus patients undergoing routine surveillance endoscopy, sampling metaplastic tissue alongside normal squamous esophagus and gastric cardia tissues when clinically appropriate. Our novel approach integrated single-cell RNA sequencing of nearly 180,000 cells with mitochondrial DNA mutation analysis for lineage tracing, coupled with spatial transcriptomics and comprehensive whole-exome sequencing. For one patient with high-grade dysplasia, we performed detailed whole-exome sequencing analysis at 100x coverage on dissociated Barrett's esophagus cells alongside matched normal esophagus and gastric cardia tissues. This genomic profiling revealed critical genetic drivers of malignant transformation, including a truncating mutation in the mutation cluster region of APC (a negative WNT regulator), mutations in the tumor suppressor genes CDKN2A and TP53, and significant copy number alterations including loss of CDKN2A on chromosome 9. These genetic alterations corresponded with the transcriptional changes observed in our single-cell data, particularly WNT pathway dysregulation and the expansion of LGR5-expressing stem cells, providing molecular evidence for how single clones can initiate the dysplastic transformation. Data deposited in dbGaP consist specifically of the whole-exome sequencing raw files, processed variant calls, and copy number analyses from the single patient with high-grade dysplasia and matched normal esophagus and stomach control tissue. Whole exome sequencing was performed on an Illumina NovaSeq 6000. Sequencing data was preprocessed following the Genome Analysis Toolkit's (GATK) Best Practices. Sequencing data was analyzed with Mutect2 (GATK, v4.2.5.0) to identify short somatic mutations including single-nucleotide polymorphisms and insertions and deletions. The Barrett's esophagus sample was run with the two matched normals, as well as the publicly available germline resource somatic-hg38_af-only-gnomad.hg38.vcf.gz (GATK). Somatic variant calls were filtered using FilterMutectCalls (GATK) and loaded into the Integrative Genomics Viewer (v2.12.3) for identification of high-quality variants in known esophageal adenocarcinoma regulators. Copy number alterations in the Barrett's esophagus sample were estimated using CNVkit (v0.9.9). CNV, Mutect2, and SVC analyses will be made available through dbGAP. Copy number alterations in the Barrett's esophagus sample were estimated using CNVkit (v0.9.9) with the default run settings; a pooled reference was generated from the normal esophagus and stomach samples. Copy number results were plotted directly within CNVkit using the function scatter.

Study phs003949

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

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

Study phs002031

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

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

Study phs003448

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

Colorectal cancer organoids expressing BRAF (fusion) genes

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

Study EGAS00001003558

Healthy_ageing_thymus

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

Study EGAS00001004311

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

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

Study EGAS00001004668

Transcriptomic Analysis of HIV-Infected Cells

Defining the attributes of rare CD4 T cells that harbor latent HIV under therapy with anti-retroviral agents has been challenging due to difficulties in identifying and isolating these cells. These cells represent an important barrier to cure, and characterizing them would aid the design and development of novel HIV cure strategies. In this study, we addressed this challenge using a custom microfluidic technology named Focused Interrogation of cells by Nucleic acid Detection and Sequencing (FIND-seq) that isolates HIV-infected cell transcriptomes, based on HIV DNA detection. We applied this method to HIV DNA+ memory CD4 T cells isolated from blood from five persons with HIV who were receiving anti-retroviral agents. After curating samples based on data quality, we proceeded with bioinformatics analysis from three persons. This analysis identified six transcriptomic pathways that were inhibited in HIV DNA+ cells compared to HIV DNA- cells. These included death receptor signaling, necroptosis signaling, and anti-proliferative G-alpha12/13 signaling. Furthermore, we identified two gene clusters, including 60 and 85 genes, that were significantly associated with HIV DNA+ cells identified by co-expression network analysis. These genes included negative regulators of HIV transcription that were higher in HIV DNA+ cells, positive regulators of HIV transcription that were lower in HIV DNA+ cells, and other genes involved in the negative regulation of mRNA translation, additional RNA processing functions, and the regulation of cellular state and fate. We further examined these signatures in transcriptomic data from CD4 T cell subsets sorted by flow cytometry from nine persons with HIV. We found a partial similarity between the smaller cluster of genes and the signature of CCR6- peripheral TFH cells. Raw sequencing data from all eleven human participant samples included in the study are available. These include samples from two participants who were studied only by FIND-seq, three participants who were studied by both FIND-seq and flow cytometry and six participants who were studied only by flow cytometry.

Study phs003095

Adult Eosinophilic Esophagitis Registry Atlas

Eosinophilic esophagitis (EoE) arises in the setting of type 2 cytokine-driven inflammation, is characterized by the infiltration of eosinophils in the esophageal mucosa, and is associated with dysphagia and fibrostenotic complications. Disease development and progression involve complex interactions between epithelial cells, stromal cells and immune cells. Indeed, increased epithelial permeability is consistently found in EoE, potentially leading to increased antigen exposure and inflammation. Persistent inflammation can in turn lead to fibroblast activation and fibrosis, which, together with the mediators produced by eosinophils and mast cells, leads to dysmotility. To decipher the interactions between these cell types across distinct disease states, we obtained biopsies from 7 healthy controls and 15 EoE subjects, including both active disease and clinical remission, to perform single cell RNA-sequencing (scRNA-seq).Each biopsy sample was enzymatically digested within 3 hours of collection with collagenase/dispase to generate a single cell suspension. Suspensions were then loaded on a 10x Chromium controller, and libraries were generated with the 3' Chemistry pipeline (v2 or v3) before sequencing on an Illumina sequencer. The resulting sequence data was demultiplexed and analyzed with the CellRanger pipeline to generate raw fastq files and processed gene-count matrices for each sample. Altogether, our study generated transcriptomic profiles from 421,312 individual cells after quality control and filtering.Our analyses reveal 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unknown features of EoE. Active disease displayed enrichment of ALOX15+ macrophages, novel PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Receptor–ligand expression uncovered eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures included mast and CD4+ Tissue-resident memory T (TRM) cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors.

Study phs003574

Genomic Basis of Phenotypic Variability of Complex Disorders

The 520-kbp deletion on chromosome 16p12.1 is associated with multiple neurodevelopmental outcomes, including intellectual disability/developmental delay (ID/DD), schizophrenia, autism, and epilepsy. This variant is inherited in >95% of cases from carrier parents manifesting neuropsychiatric features, while affected children with the deletion were more likely to carry another large Copy Number Variation (CNV) or rare deleterious mutation elsewhere in the genome, or "second-hits", compared to carrier parents. Thus, while the 16p12.1 deletion confers significant risk for a range of disorders, the ultimate phenotypic trajectory is determined by "second-hit" variants in the genetic background. However, outstanding questions remain about the genetic background in these disorders, including the functional effects of "second-hit" variants towards gene expression or developmental phenotypes, the full spectrum of "second-hit" variants, and the mechanisms driving accumulation of these variants.In version 1, we performed deep clinical phenotyping, whole-genome sequencing (WGS), SNP microarray analysis, and RNA-sequencing of lymphoblastoid cell lines of 32 individuals in five large families with the 16p12.1 deletion. We found that specific classes of rare single-nucleotide and structural variants in coding and non-coding regions contributed towards expression changes in affected children based on their inheritance patterns, including a subset of variants that acted synergistically with the deletion to alter expression. We further found that rare variant-mediated expression changes contribute to specific developmental phenotypes in the carrier children, suggesting a potential genetic mechanism for variable expressivity of the deletion.In version 2, we additionally analyzed genetic and phenotypic data from individuals with 16p12.1 deletion and their spouses to assess the role of assortative mating modulating phenotypic outcomes of the 16p12.1 deletion. We identified significant within- and cross- disorder phenotypic correlations between spouses, but observed no significant genetic correlations between spouses, in line with expectations.

Study phs002450

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

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

Study EGAS00001006662

3883 results for "RNA AND sequencing OR transcriptome"

in 55.21 milliseconds.

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

Generation of a Cellular Atlas of Human Adipose Tissue

Genomics and Epigenomics of the Elderly Response to Pneumococcal Vaccines

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

3D tissue engineered human skeletal muscle modelling Facioscapulohumeral Muscular Dystrophy

Pediatric Non-Down Syndrome Acute Megakaryoblastic Leukemia is Characterized by Distinct Genomic Subsets with Varying Outcomes

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

Clonal_haematopoiesis_in_patients_with_AAA

Cell type mapping of inflammatory muscle diseases highlights selective myofiber vulnerability in inclusion body myositis

Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

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

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

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

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

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

Investigating Delayed-Onset Drug Hypersensitivity Reactions Prospectively

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

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

Exceptional Responders Initiative

Cell and Circuit-Specific Exploration of HIV Neurogenomics in Context of Opiate and Cocaine Misuse

Effect of ETS2 modulation on chromatin accessibility and enhancer activity in inflammatory (TPP) macrophages

DNA and RNA sequencing of single human haploid germ cells

Multivalent State Transitions Shape the Intratumoral Composition of Small Cell Lung Carcinoma

Blood DNA Methylation in Post-Acute Sequelae of COVID-19 (PASC): a Prospective Cohort Study

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

Biallelic tumor suppressor loss and DNA repair defects in de novo small cell prostate cancer

Pharmacogenomic Analysis Reveals New Therapeutic Options for Pleural Mesothelioma

Dataset for Sarcoma RNA sequencing linked from study EGAS00001004813

The Genomic and Transcriptomic Landscape of a HeLa Cell Line

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

Analysis of CD20 loss in patients treated with Mosunetuzumab

T-cell receptor sequencing of tumor-infiltrating lymphocytes (TIL) in breast cancer

Multiomic Landscape of Multiple Myeloma Precursor and Relapsed Disease

Single-cell RNA-seq of peripheral blood mononuclear cells in classic Hodgkin lymphoma

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Exploring the driver events of eccrine poromas and porocarcinomas: A retrospective, cross-institutional study of 90 cases

Multiple Myeloma treated with anti-BCMA T-cell-redirecting therapy (1st Batch)

BRAIN Cell EncyclOpeDia of Transcribed Elements (BRAINcode)

Integrated Metabolic Profiling and Gene Expression Analysis Reveals Therapeutic Modalities in Breast Cancer

Genetic and Hormonal Contributions to Gene Expression in Immune Cells

Investigating Host Genetic Risk Factors for Tuberculosis in Highly Endemic South African Populations

INSPIRE - Investigator-initiated Phase 2 Study of Pembrolizumab Immunological Response in Metastatic Solid Tumors

A prognostic human brain network for diffuse midline glioma

scRNA-seq of HGSC tumor and ascites

Single cell transcriptomics of hESC-derived midbrain dopaminergic neurons generated by a new human development-based protocol

BMP4 and temozolomide synergize in the majority of patient derived glioblastoma cultures

Resistance to Checkpoint Blockade Therapy Through Inactivation of Antigen Presentation

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

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

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

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

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

Addiction of Primary Cutaneous Gamma-Delta T-cell Lymphomas to JAK-STAT Signaling

Paired DNA and RNA sequencing uncovers common and rare genomic variants regulating gene expression in the human retina

RNA sequencing data of 142 samples from 142 patients with HER2+ breast cancer treated with letrozole or tamoxifen (SOLTI PAMELA trial)

RNA sequencing of circulating human immune cells before and after interleukin-2 immunotherapy in systemic lupus erythematosus patients

TSO500 STJAN33, BRAF mutated CUP

RNA sequencing data of pediatric KMT2A-rearranged acute lymphoblastic leukemia

RNA sequencing data of pediatric T-cell acute lymphoblastic leukemia

Alignment BAM files and gene count files of the Illumina Sequencing Data (10 tumor samples)

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

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

Exome sequencing

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

Genomic profiling of paediatric high grade gliomas from the HERBY clinical trial

The G2 gene expression signature and MYC overexpression are independent poor prognostic factors in childhood high-grade osteosarcoma

RNA-seq study of longitudinal blood cell samples from children at risk of type 1 diabetes

IFN-γ and TNF-α drive a CXCL10+ CCL2+ Macrophage Phenotype Expanded in Severe COVID-19 Lungs and Inflammatory Diseases with Tissue Inflammation

Anaplastic oligodendroglioma exome and RNA sequencing data

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

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

Relapse CHL study

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

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

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

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

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

Identification of Somatic Changes in Tumors from Fanconi Anemia Patients