3924 results for "RNA AND sequencing OR transcriptome"

in 57.63 milliseconds.

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

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

  Dataset EGAD00001004140

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

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

  Study EGAS50000001047

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

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

  Study EGAS00001005405

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

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

  Study EGAS00001007362

• Fecal microbiome of T2D patients undergoing semaglutide or empagliflozin treatment

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

  Dataset EGAD50000000756

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

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

  Dataset EGAD00001003984

• Genomic Analysis of Prostate Tumor Heterogeneity in Metastasis

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

  Study phs003404

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

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

  Study EGAS50000000175

• Programs, Origins, and Immunomodulatory Functions of Myeloid Cells in Gliomas

  Gliomas are challenging malignancies characterized by an immunosuppressive microenvironment filled with myeloid cells, whose roles in immune modulation are not fully understood. In this study, we utilized single-cell RNA sequencing (scRNA-seq) data from 85 primary human glioma tumors, related blood samples, and brain tumor organoid models. We found that nearly all glioma-associated myeloid cells express at least one of four distinct immunomodulatory activity programs. These programs include:1. Scavenger Immunosuppressive Program2. Complement Immunosuppressive Program3. Microglial Inflammatory Program4. Systemic Inflammatory ProgramAll four programs are found in both lower-grade and high-grade gliomas and are expressed across various myeloid cell types, whether derived from blood or resident origins. We validated the remarkable plasticity of myeloid cells by showing that the introduction of peripheral blood monocytes into glioma organoid models can induce the emergence of identities resembling microglia, macrophages, and dendritic cells, as well as all four immunomodulatory programs. Our study provides a valuable resource and a new framework for understanding the role of immunomodulatory myeloid cells in gliomas, laying the groundwork for developing effective immunotherapy strategies for glioma patients.

  Study phs003756

• Defining Cutaneous Gene Expression Signatures in Juvenile Dermatomyositis

  Skin inflammation in juvenile dermatomyositis (JDM) can signal disease onset or flare up and prevent complete disease remission. Our objective was to evaluate JDM cutaneous expression signatures through tape stripping in both lesional and non-lesional skin and as compared to paired blood. A total of 28 JDM and 20 control (CTL) patients participated in the study. All patients had tape strip sampling of non-lesional (NL) skin at enrollment visit. Within our JDM cohort, paired lesional (L) skin was additionally sampled if rash was present. 16/28 of JDM patients had at least one additional skin sample associated with a follow-up (FU) visit (total number of follow-up samples included = 26 NL and 4 L). We identified a common innate immune signature in both lesional and non-lesional JDM skin. Using unsupervised hierarchical clustering of skin expression data, we distinguished a JDM molecular subgroup characterized by a mitochondrial dysfunction signature. In comparison to blood expression data, we determined that skin may more effectively highlight the JDM subgroup characterized by a mitochondrial dysfunction expression signature. Data available include bulk RNA-sequencing data from tape strips and blood of the JDM and control patients from this study.

  Study phs003884

• Omic Studies in Children's Respiratory and Environmental Workgroup (CREW) Cohorts

  Children's Respiratory Research Workgroup (CREW) is an NIH-funded project consisting of multiple individual U.S. cohorts and three scientific centers working together to identify phenotypes and causes of childhood asthma. CREW will include data from a large number of children and their families, with broad diversity in terms of ethnicity, family characteristics, neighborhoods and geographic locations. One of the primary goals of CREW is to put together sets of data and samples of participating cohorts to identify phenotypes of childhood asthma (i.e. specific subtypes of asthma that can be distinguished by clinical features such as natural history, triggers, exacerbation frequency, concurrent allergies, lung function, sex, etc). Omics studies will include genotyping, methylation, and RNA sequencing. By integrating the DNA genotyping with the expression, immunologic and clinical data, we will identify eQTL, inflammation QTL, and lung function QTL, as well as QTL associated with other biomarkers. Analyses will identify genetic, epigenetic, transcriptional and protein/metabolite profiles that are associated with specific exposures or phenotypes. As we obtain mechanistic insights about personal and early life risk factors, we will connect asthma phenotypes with underlying causes and pathogenic mechanisms to define endotypes of childhood asthma.

  Study phs004036

• Biallelic loss of function PAX4 variants are a Novel Cause of Transient Neonatal Diabetes

  The identification of genes in which loss of function (LoF) variants cause neonatal diabetes mellitus (NDM) can provide unique insights into the development and function of human pancreatic beta cells. We describe the identification of biallelic PAX4 LoF variants in 2 unrelated individuals with transient NDM: one due to a homozygous p.(Arg126*) PAX4 stop gain variant and a second due to a homozygous c.-352_104del 2.65kb deletion affecting the first 4 exons and promoter of the PAX4 Matched Annotation from NCBI and EMBL-EBI select transcript. We confirmed reduced transcript expression due to nonsense mediated decay for the e p.(Arg126*) variant in CRISPR-edited iPSC-derived pancreatic endoderm cells. Analysis of PAX4 binding through integration of CUT&RUN and RNA-sequencing in PAX4 depleted endocrine cells identified genes directly regulated by PAX4 involved in both pancreatic islet development and glucose-sensitive insulin secretion. The remitting and relapsing diabetes in both probands demonstrates that in contrast to mice PAX4 is not essential for the development of pancreatic beta-cells and supports differences in the requirements of functional beta-cell mass across the life-course.

  Study EGAS50000000857

• Single cell atlas of large B-cell lymphoma

  Large B-cell lymphomas (LBCL) are clinically and biologically heterogeneous lymphoid malignancies, shaped by complex tumor and microenvironments that play a critical role in disease progression. There has been no direct and comprehensive single-cell characterization of tumor microenvironment in LBCL that simultaneously captures lymphoid, myeloid, and stromal populations across a sufficiently large cohort to reveal reproducible and biologically meaningful patterns. Here, we performed single-nucleus multiome profiling, coupled with bulk DNA and RNA sequencing, on 232 tumor and control biopsies to comprehensively map cell subpopulations of LBCL. Through this integrative approach, we identified Lymphoma Microenvironment Archetype Profiles (LymphoMAPs)—distinct, reproducible patterns defined by the co-occurrence of functionally specialized cell subsets. We further characterized the cell-cell communication pathways that shape these transcriptional phenotypes within archetypes. Importantly, we discovered associations between specific LymphoMAPs and clinical and molecular features of LBCL, including response to CD19 CAR T-cell therapy. Our study provides a valuable resource of lymphoid, myeloid, and non-hematopoietic cell states in LBCL and offers a framework for understanding how the abundance and functional state of cells within the tumor microenvironment may drive inter-patient variability in therapeutic outcome.

  Study EGAS50000001022

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

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

  Study EGAS00001003736

• Biased allelic expression in human primary fibroblast single cells.

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

  Study EGAS00001001009

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

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

  Study EGAS00001004614

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

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

  Study EGAS00001005675

• Oncoprint GSCCs

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

  Study EGAS00001007481

• Molecular_risk_stratification_in_patients_with_T1_colorectal_cancer_WES

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

  Study EGAS00001005734

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

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

  Study EGAS00001005874

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

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

  Study EGAS00001006747

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

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

  Study EGAS00001006818

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

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

  Study EGAS00001007482

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

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

  Dataset EGAD00001014788

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

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

  Dataset EGAD00001014790

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

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

  Study EGAS00001005039

• Whole genome sequencing data from tumor and normal samples

  Whole genome sequencing FASTQ files from tumor and normal samples from 202 patients. All experiments reported in this article were performed in the Clinical Laboratory Improvement Amendments (CLIA)-certified and College of American Pathologists (CAP)-accredited laboratories at Personalis Inc., as guided by the Association for Molecular Pathology (AMP) and CAP’s joint recommendations52. Preferred input material consisting of FFPE slides (enabling H&E and macrodissection which is part of the standard NeXT Personal process) was not available for the majority (94%) of tumor samples. Most tumor samples (72%) were processed from FFPE curls, with the remaining processed from pre-extracted DNA and fresh frozen material. Where FFPE material on slide was available, tumor sections were macrodissected to increase tumor content, with a minimum allowable cellularity threshold of 20% determined by pathological review. The tumor cellularity threshold of 20% was chosen based on prior work demonstrating that no significant correlation is observed between tumor purity and assay limit of detection (LOD) at this threshold31. Genomic DNA was isolated from tumor and normal samples and WGS libraries prepared as previously described31. Genomic DNA was isolated from matched tumor and normal samples using the Qiagen AllPrep DNA/RNA FFPE Tissue Kit or the QIAamp DNA Mini Kit (QIAGEN, Germantown, MD, USA) using internally optimized workflows. WGS sequencing libraries were prepared with 120-500 ng of acoustically sheared genomic DNA (Covaris LLC, Woburn, MA, USA) using the KAPA HyperPrep Kit (Roche Sequencing Solutions, Pleasanton, CA, USA) and customized methods. Libraries were cleaned-up using AMPure XP beads and then quantified using the KAPA Library Quantification Kit (Roche Sequencing Solutions, Pleasanton, CA, USA), before being sequenced to 30X depth of coverage using a NovaSeq 6000 instrument (Illumina, San Diego, CA, USA).

  Dataset EGAD50000001909

• Determination of cell specific regulatory enhancers in hematopoetic models

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

  Dataset EGAD00001002205

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

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

  Dataset EGAD00001010018

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

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

  Study EGAS00001003274

• Exceptional Responders Initiative

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

  Study phs001145

• Renal_habitat_WXS

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

  Study EGAS00001003703

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

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

  Study JGAS000179

• 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

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

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

  Study phs002925

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

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

  Study phs000642

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

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

  Study EGAS00001004289

• 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

• Access to "A Spatially Resolved Single-Cell Atlas of the Human Fetal Olfactory System"

  The human nasal region arises from neural crest and placodal lineages, yet its early development remains poorly understood owing to limited fetal tissue access and structural complexity. Here we present an integrated single-nucleus and spatial transcriptomic atlas of the human fetal nasal region, generated from 10 fetuses (males and females) between 7 and 12 post-conceptional weeks. Single-nucleus RNA sequencing resolved 32 distinct cell types, while integration with multiplexed error-robust fluorescence in situ hybridization (MERFISH) enabled spatial and temporal mapping of gene expression dynamics across the olfactory epithelium (OE) and adjacent tissues. We identify novel markers of olfactory sensory neuron differentiation and pathways governing epithelial patterning and OE morphogenesis. Notably, spatially organized expression of 169 olfactory receptor genes reveals molecular evidence for the “one neuron-one receptor” principle during the first trimester. Together, this work establishes the first molecular and spatial framework of early human olfactory development and provides a foundational resource for studies of sensory neurogenesis and congenital disorders.

  Dac EGAC50000000688

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

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

  Study EGAS00001007110

• Genetics of gene expression in human macrophage response to Salmonella

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

  Dataset EGAD00001003204

• The Transcriptional Landscape of SHH Medulloblastoma

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

  Dataset EGAD00001006305

• 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

• 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

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

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

  Study phs001897

• 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

• 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.&nbsp; 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.&nbsp;

  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 (&#177; 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

• 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

• 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

• 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

• 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

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

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

  Study EGAS00001005478

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

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

  Study EGAS00001006678

• 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&lt;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

• 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 &lt;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

• 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

• Common Fund (CF) Genotype-Tissue Expression Project (GTEx)

  Lay Description The aim of the Genotype-Tissue Expression (GTEx) Project is to increase our understanding of how changes in our genes affect human health and disease with the ultimate goal of improving health care for future generations. GTEx will create a database that researchers can use to study how inherited changes in genes lead to common diseases. GTEx researchers are studying genes in different tissues obtained from many different people. The GTEx project also includes a study of the GTEx donor consent process - this study will help ensure that the consent process and other aspects of the project effectively address the concerns and expectations of participants in the study. GTEx is a pioneering project that uses state-of-the-art protocols for obtaining and storing a large range of organs and tissues, and for testing them in the lab. Until now, no project has analyzed genetic variation and expression in as many tissues from the same person in such a large population as planned for GTEx. Scientific Description Understanding the role of variation in the human genome is crucial to elucidating genetic contributions to human health and disease. Despite the results of genome-wide association studies (GWAS) documenting strong statistical associations between genetic variation and human traits, the functional role for most of these variants is largely unexplained. Nearly 90% of these GWAS-implicated sites lie outside of protein-coding sequences, suggesting that these variants might regulate gene expression. The goal of the Genotype-Tissue Expression (GTEx) project is to establish a resource database and tissue biobank in which to study the relationship between genetic variation and gene expression and other molecular phenotypes in reference/non-diseased tissues. The ultimate resource will include approximately 960 post-mortem donors with several dozen tissues from each, a resource large enough to study both cis- and trans- gene expression quantitative trait loci (eQTLs). Some tissue will also be banked for additional molecular analyses. To request biosamples, please see the Biobank page for more information. GTEx was initially funded as a 2-year pilot project by the NIH Common Fund (CF) in 2010, and has been scaled up after demonstration of feasibility. The project will collect and analyze RNA levels in many different human tissues and each donor will be characterized for germline genetic variation through dense genotyping arrays and sequencing of either whole exomes or whole genomes. By treating RNA expression levels as quantitative traits, eQTLs will be identified as sites containing genetic variation that correlate with changes in RNA expression. Such eQTLs have been associated with 4%-12% of expressed human genes, and with common complex human diseases, including obesity, atherosclerosis, type 2 diabetes, Crohn's disease, and asthma. Additionally, few studies have examined the tissue specificity of eQTLs. A subset of banked tissue samples will also be analyzed for other molecular phenotypes, such as DNA methylation, DNaseI hypersensitivity sites, and proteomics. The GTEx project will thus serve as a resource database and tissue bank for many future studies, especially for understanding the functional basis of inherited susceptibility to disease. All GTEx releases since Version 5 follow the NIH Genomic Data Sharing (GDS) Policy whereby there are no restrictions on use or publication after release. Additional information about the GTEx Data Release and Publication policy can be found on the GTEx Data Portal at www.gtexportal.org.

  Study phs000424

• 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

• 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

• 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

• 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

• 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 (&lt;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 &gt; 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

• 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&#947; 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&gt;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 (&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 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 (&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 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;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

• 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&#39;s or Alzheimer&#39;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&#39; 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

• 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

• 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

• 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&#38;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, &#8220;Daratumumab-Based Response-Adaptive Therapy for Older Adults With Newly Diagnosed Multiple Myeloma.&#8221; 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 &#8216;a') with Daratumumab and Dexamethasone, and at cycle 2 day 22 (C2D22, i.e., 60 days after induction therapy start, denoted by &#8216;b'). For three NDMM patients, we also have samples at relapse (denoted by &#8216;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

• 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&#8217; 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

• 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

• 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-&#946; and IFN-&#947;), pro-inflammatory cytokines such as TNF-&#945;, 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-&#947; and TNF-&#945; 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

• 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