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

in 89.21 milliseconds.

• 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 sequencing CYLD cutaneous syndrome

  Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. We comprehensively profiled the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families. Novel recurrent mutations were found in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumours suggest that isoform specific DNMT3A2 mutations are associated with dysregulated Wnt/ß-catenin pathway signaling. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumours. These findings add novel dimensions to existing paradigms of cutaneous tumorigenesis and metastasis.

  Study EGAS00001003841

• Single cell characterization of T-cell lymphoma: RNA (2025-07-31)

  Lymphomas are types of blood cancer derived from lymphocytes, and can be classified into two main categories, Hodgkin's lymphomas (HL) and non-Hodgkin lymphomas (NHL). T-cell lymphoma is a type of NHL that develops from T lymphocytes and accounts for 10-15% of total NHL cases. As a rare disease, T-cell lymphoma is functionally, pathologically, and clinically complex and heterogeneous. The 2016 revision of the WHO classification has listed 29 subtypes of mature T/NK-cell neoplasms, among which over 25 are T lymphocyte in origin. These subtypes of T-cell lymphomas have different biological behaviours and clinical prognosis. Clinically, T-cell lymphoma ranges from indolent (slow-growing) to aggressive (fast-growing and spreading) disease. According to the guideline booklet provided by the Lymphoma Research Foundation (https://lymphoma.org/aboutlymphoma/nhl/), indolent T-cell lymphoma subtypes include adult T-cell leukaemia/lymphoma (ATLL), cutaneous T-cell lymphoma (CTCL), and mycosis fungoides (MF), while the aggressive subtypes include anaplastic large cell lymphoma (ALCL), angioimmunoblastic T-cell lymphoma (AITL), peripheral T-cell lymphoma (PTCL), and T lymphoblastic lymphoma (TLBL). PTCL and ALCL represent the most common subtypes of aggressive T-cell lymphoma whereas TLBL is a rare and aggressive subtype mostly diagnosed in children. In this study, we will use single cell RNA sequencing and single cell multiomic methods three aggressive T-cell lymphoma subtypes (ALCL, PTCL and TLBL) and try to gain a deeper insight into the genotypic and phenotypic features of these T-cell lymphomas. . This dataset contains all the data available for this study on 2025-07-31.

  Dataset EGAD00001015669

• Title: Divergent levels of CD112 and INKA1 define a distinct subset of human long-term hematopoietic stem cells

  The data provided here was critical in establishing that human long-term hematopoietic stem cells (LT-HSC), previously described as the most primitive HSC population, is actually composed of distinct subsets that can be prospectively isolated. Via mechanistic studies centering around the Rho-GTPase effector kinase PAK4 and its inhibitor INKA1, we identified the immune checkpoint ligand CD112 as a marker for hematopoietic stem and progenitor cells, that is highest expressed on LT-HSC. More importantly, CD112 can be used to stratify functionally distinct subsets within LT-HSC: In response to regeneration-mediated stress, the CD112low subset exhibits a transient restraint (termed latency) before contributing to hematopoietic reconstitution, while the CD112high subset is primed to respond rapidly. High resolution RNA-seq of the CD112 surface expression spectrum within rare LT-HSC subsets (human umbilical cord blood) demonstrated that more genes are differentially upregulated in the deeper quiescent and less metabolic active subset. Genes enriched in this subset centre around cell adhesion and Rho-GTPase signaling. This is in agreement with the scRNAseq data from human G-CSF mobilized peripheral blood (mPB) generated here that was used as an model of in vivo activation/priming revealing via RNA-velocity and pseudo-time analysis that INKA1high versus PAK4high, CDK6high and CD112high enrichment are either detected early or late in diffusion pseudotime indicative of quiescent versus primed cell status, respectively. RNAseq following INKA1 overexpression in LT-HSC and ST-HSC revealed by GSEA an overall stemness preserving phenotype and particularly in LT-HSC, but not in short-term HSC (ST-HSC), suppression of transcriptional programs linked to activation. Collectively, our data decipher the molecular intricacies underlying HSC heterogeneity and self-renewal regulation and point to latency as an orchestrated physiological response that integrates quiescence control with HSC fate choices to preserve a stem cell reservoir.

  Dataset EGAD00001006541

• Dual targeting of polyamine synthesis and uptake in diffuse intrinsic pontine gliomas

  Bulk RNA-seq data of tumour data.

  Dataset EGAD00001006794

• A Phase I Study of the Treatment of Recurrent Malignant Glioma with CAN-3110 (AKA rQNestin34.5v.2), a Genetically Engineered HSV-1 Virus

  This study is a "first-in-human" phase 1 trial of recurrent, IDH wild-type glioblastoma (rGBM) patients treated with CAN-3110 (aka rQNestin34.5v.2), a genetically modified oncolytic herpes simplex virus 1 (HSV-1). The study was a 3+3 dose-escalation study starting at a single intra-tumoral injection of 106 PFU CAN-3110 going up to 1010 PFU in half-log increments (cohorts 1-9). An expansion cohort (cohort 10) was treated with 109 PFU injected into up to 5 sites within a single tumor. Principal findings of the study conclude that:CAN-3110 was safe at all tested doses (no dose-limiting toxicities reached)HSV-1 seropositive patients survived significantly longer than HSV-1 seronegative patients in this study - longer than expected based on historical rGBM controlsHSV-1 seropositive patients were significantly more likely to clear CAN-3110 (HSV-1) antigen from their tumors based on immunohistochemical stainingCAN-3110 treatment-induced CD8+/CD4+ T-cell recruitment to injected tumor sites altered TCR beta diversity in both tumor infiltrating leukocytes (TILs) and PBMCs, and initiated changes in specific TCR beta clonotypes and immune gene signatures which were associated with post-treatment survival.Data provided in dbGaP for this study include:Bulk RNA-Seq data for 13 pre/post-treatment pairs of rGBM tumors corresponding to 12 unique patients (patient 28 was treated with CAN-3110 at two different timepoints). This data is provided as raw sequencing reads (.fastq), kallisto pseudoalignment transcript-level read counts (abundance.tsv files), gene-level abundance matrix (Gene_Level_Expression_Matrix.tsv), and immune gene signature scores matrix (Immune_Gene_Signatures_Matrix.tsv). Sequencing was performed as paired 150bp sequencing on NovaSeq6000 S4 Illumina flow cells. Library preparation varied by batch as described in the included RNAseq_Metadata.xlsx table.TCR beta DNA sequencing performed by Adaptive Biotechnologies for 21 pre/post-treatment pairs of rGBM tumors and PBMCs representing 20 unique patients (patient 28 was treated with CAN-3110 at two different timepoints). Note that all tumor pairs with TCR beta sequencing data also have available bulk RNA-Seq data. The TCR beta sequencing data is provided as an excel spreadsheet (TCRbeta_Sample_Overview.xslx), which provides a broad level summary of the sequencing on a per-sample basis (# of TCRbeta templates, T cell fraction, TCR beta diversity metrics, etc.) and as a large tab separated value file (Full_TCRbeta_Rearrangement_Details.tsv), which provides information per sample at the level of each unique TCR beta sequence (amino/nucleic acid sequence, VDJ gene usage, etc.).An excel table of de-identified metadata for the patients involved in this study.

  Study phs003378

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

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

  Study phs003152

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

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

  Study phs000770

• An RCOR1 loss-associated gene expression signature identifies a prognostically significant DLBCL subgroup

  Effective treatment of diffuse large B cell lymphoma (DLBCL) is plagued by heterogeneous responses to standard therapy, and molecular mechanisms underlying unfavorable outcomes in lymphoma patients remain elusive. Here, we profiled 148 genomes with 91 matching transcriptomes in an R-CHOP-treated DLBCL cohort to uncover molecular subgroups linked to treatment failure. Systematic integration of high-resolution genotyping arrays and RNA-seq data revealed novel deletions in RCOR1 to be associated with unfavorable progression free survival (P = 0.001). Integration of expression data from the clinical samples with data from RCOR1 knockdowns in the lymphoma cell lines KM-H2 and Raji yielded an RCOR1 loss-associated gene signature comprising of 233 genes. This signature identified a subgroup of patients with unfavorable overall survival (OS) (P = 0.023). The prognostic significance of the 233-gene signature for OS was reproduced in an independent cohort comprising of 195 R- CHOP treated patients (P = 0.039). Additionally, we discovered that within the IPI low risk group, the gene signature provides additional prognostic value that was independent of the cell-of-origin phenotype. Taken together, we present a novel and reproducible molecular subgroup of DLBCL, impacting risk-stratification of R- CHOP treated DLBCL patients and revealing a new avenue for therapeutic intervention strategies.

  Study EGAS00001001000

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

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

  Study EGAS00001004343

• Adult-type Granulosa Cell Tumour of the Ovary

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

  Study EGAS00001005414

• Stereotyped B-cell responses are linked to IgG constant region polymorphisms in multiple sclerosis

  Clonally related B cells infiltrate the brain, meninges, and cerebrospinal fluid of multiple sclerosis (MS) patients, but the mechanisms driving the B-cell response and shaping the immunoglobulin repertoires remain unclear. Here, we used single-cell full-length RNA-seq and B-cell receptor reconstruction to simultaneously assess the phenotypes, isotypes, constant region polymorphisms, and the paired heavy- and light-chain repertoires in intrathecal B cells. We detected extensive clonal connections between the memory B cell and antibody-secreting cell compartments and observed clonally related cells of different isotypes, including IgM/IgG1, IgG1/IgA1, IgG1/IgG2, and IgM/IgA1. There was a strong dominance of the G1m1 allotype constant region polymorphisms in antibody-secreting cells, but not in memory B cells. Tightly linked to the G1m1 allotype, we found a preferential pairing of the immunoglobulin heavy-chain variable (IGHV)4 gene family with the κ variable (IGKV)1 gene family. The IGHV4-39 gene was most used and showed the highest frequency of pairing with IGKV1-5 and IGKV1(D)-33. These results link IgG constant region polymorphisms to stereotyped B-cell responses in MS and indicate that the intrathecal B-cell response in these patients could be directed against structurally similar epitopes.

  Study EGAS00001005745

• Characterization of a New Case of XMLV (Bxv1) Contamination in the Human Cell Line Hep2 (Clone 2B) Using High-Throughput DNA- and RNA-Seq

  We have sequenced DNA and RNA from the Hep2 (clone 2B) cell line (Sigma-Aldrich catalog number: 85011412-1VL) using short read Illumina sequencing. We show that this cell line is indistinguishable from the HeLa cell line (as mentioned on the Sigma-Aldrich website) by mapping integrations of the human papillomavirus 18. We also show that the cell line is contaminated by the Bxv1 xenotropic murine leukemia virus. This contamination has to be taken into account when using the cell line in future experiments.

  Study phs001944

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

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

  Study EGAS00001005513

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

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

  Study EGAS00001006369

• Single Cell RNAseq at various stages of HiPSCs differentiating toward definitive endoderm and endoderm derived lineages

  When comparing the differentiation capacities of pluripotent stem cell lines that have different genetic backgrounds, batch to batch experimental variablility poses a significant challenge, especially when trying to identify smaller effects. One way to address this issue is to differentiate several different lines in the same culture dish, thereby elimating experimental variation. In addition, it allows researchers to analyze many more lines with less experiments. Parallel single cell RNA-Seq exploits that individual cells are tagged and hence each cell can be reliably assigned to the donor of origin based on the genetic variants it contains. In addition, analyzing the genetic signature of single cells within a differentiating population can reveal differentation stages that are not easily detected in bulk RNAseq data. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Dataset EGAD00001005741

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

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

  Dataset EGAD00001008969

• NHGRI GREGoR Consortium: Genomics Research to Elucidate the Genetics of Rare Disease

  The NHGRI GREGoR (Genomics Research to Elucidate the Genetics of Rare Disease) Consortium was established in June 2021 with the goal of developing novel tools and approaches to advance the discovery of the genetic basis of rare conditions. Numerous types of molecular data are generated in GREGoR and available on the AnVIL cloud platform via dbGaP application, including short- and long-read genome and exome sequencing, transcriptomics, metabolomics, methylomics, and proteomics. De-identified clinical and demographic data is obtained, with a focus on standardized ontologies.Visit the GREGoR Consortium data webpage for summary information about the GREGoR Dataset, including numbers of participants and data types, methods documentation, and Release Notes. The Consortium comprises five Research Centers (RCs - see below), a Data Coordinating Center (DCC), and various partner members and external collaborators.Baylor College of Medicine Research Center (BCM-GREGoR) The Baylor College of Medicine GREGoR program, which is part of the GREGoR consortium, enrolls individuals, families, and cohorts with suspected rare disease across a range of syndromic and non-syndromic phenotypes. Subjects are enrolled from national and international collaborating physician referrals. Subjects provide written informed consent for future re-contact. Data generated and shared include family structure, detailed phenotypes, exome or short-read genome data, and in some cases long-read genome or RNA-sequencing, and these are shared upon completion of standard quality control checks and annotation. Broad Institute (Broad) The Broad Center for Mendelian Genomics, part of the GREGoR consortium uses next-generation sequencing (exome, genome, transcriptome, and long read sequencing), computational approaches, and functional studies to discover the variants and genes that underlie Mendelian conditions with a particularly focus on neuromuscular, neurodevelopmental, and syndromic phenotypes. Samples come from collaborators and direct enrollment through the Rare Genomes Project and we are committed to rapid data sharing without an embargo period. University of California, Irvine (UCI-GREGoR) To accelerate the pace of Mendelian disease gene discovery and clinical implementation, we propose a Mendelian Genomics Research Center, part of the GREGoR Consortium, leveraging the broad pediatric and adult clinical and research expertise at Children's National Hospital and University of California, Irvine. Our goal is to develop best practices to increase the diagnostic yield of rare diseases, engage the community to reduce health disparities for complex diagnoses, while creating a dataset accessible to all. Our Center will unite world class experts combining basic and translational research with innovative approaches to phenotyping, variant identification and functional investigation of both coding and non-coding sequence changes with the goals of discovering novel Mendelian gene variations and identifying variants not detected on current sequencing pipelines, disambiguating uncertain variants into disease-causing versus benign categorizations, and sharing information by working collaboratively with the GREGoR community.GREGoR Stanford Site (GSS) The goal of the GREGoR Stanford Site (GSS) is to provide a platform for functional genomics research and validation to improve diagnosis in Mendelian disease. Participants included individuals with undiagnosed suspected Mendelian disease who had non-diagnostic exome sequencing and their immediate family members. Participants and their family members provided written, informed consent and biological samples from which DNA, RNA, plasma, fibroblasts, PBMCs and other cell types were generated and stored. Samples from research participants and their immediate family members may have undergone short and long-read genome sequencing, transcriptome sequencing, metabolomics and/or lipidomics profiling, methyl-capture-sequencing and ATAC-sequencing. De-identified clinical data extracted from participant medical records are linked to the samples. University of Washington Center for Rare Disease Research (UW-CRDR) The goals of the University of Washington Center for Rare Disease Research are to: (1) maximize novel gene discovery for Mendelian conditions by recruitment, short- and long-read whole genome sequencing, transcriptome sequencing and analysis of families with rare conditions for which the gene is either unknown or the gene is known but no pathogenic variant can be identified via clinical testing; (2) develop new strategies for gene discovery for Mendelian conditions caused by variants that are difficult to detect using conventional testing strategies, variants of unknown function effect (e.g., regulatory, structural variants) or have unusual modes of inheritance; and (3) implement high-throughput screening and targeted follow-up functional studies to prioritize and validate candidate non-coding variants. De-identified data and phenotypic information are shared via MyGene2, ClinVar, and AnVIL.

  Study phs003047

• Bulk-RNA Sequencing of high-grade pancreatic and non-pancreatic Neuroendocrine Neoplasms

  Therapeutic decisions in oncology depend on a precise pathological classification of individual neoplasms. Recent years have seen an intensification of research activities aimed at the extraction of clinically relevant information from patient-derived 'omics' data based on Machine-Learning models. However, a comprehensive training of Machine-Learning models requires sufficiently large numbers of training samples, which are usually not available for rare cancer types. The problem is worsened when individual tissues segregate into different cancer subtypes, as their discrimination would require even more training samples. Methods: Here, we report on a new data-augmentation technique to support the training of Machine-Learning models on ‘omics’ data from pancreatic neuroendocrine neoplasms (panNEN). PanNENs display all properties described above: Only about 2-3% of all pancreatic neoplasms are neuroendocrine and they fall into different subtypes with distinctly different prognosis, which makes the precise classification of such samples both difficult and important for therapy decisions. The approach reconstructs a given transcriptome based on healthy pancreatic cell type signatures and creates a Machine-Learning model that integrates the observed reconstruction error and predicted cell type proportions as features. Results: A benchmark of the deconvolution model predictions with the ground-truth found that the model could efficiently predict the sample grading, disease-related patient survival time, and differentiate between different subtypes in four panNEN and one mixed panNEN and non-pancreatic NENs dataset. We compared the predictive performance of the deconvolution-trained model to that of a model trained directly on the transcriptomic data, under inclusion of the Ki-67 classification gold standard biomarker, and found the performances to be comparable. Conclusion: Our approach serves as a data-augmentation technique to facilitate the training of Machine Learning models for rare cancer types via utilization of data from healthy origin. Additionally, the deconvolution results were clinically interpretable and further research can render the deconvolution approach an effective complementary asset for the clinical classification of neoplasms.

  Study EGAS00001004861

• 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

• Identification of Somatic Changes in Tumors from Fanconi Anemia Patients

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

  Study phs002652

• Bacterial Vaginosis, Cervical Immune Cells and HIV Susceptibility

  This was a cross-sectional study to characterize cervical immune cell composition and activation status in the presence versus absence of bacterial vaginosis (BV) using bulk messenger RNA sequencing (mRNA-Seq). Women enrolled in a longitudinal, open cohort study and engaged in a sex work in Mombasa, Kenya, were eligible. On the day of enrollment, participants self-collected vaginal swabs for generation of Gram-stained slide. The slides were scored from 0-10 using the method of Nugent and Hillier; in this system, women with scores of 0-3 are considered to have normal vaginal flora, those with scores of 4-6 are considered to have intermediate flora, and those with scores of 7-10 are considered to have BV. For this study, only women with normal flora (Nugent scores 0-3) or BV (Nugent scores 7-10) were eligible. Participants underwent speculum-assisted pelvic examination and collection of cervical biopsy specimens using Tischler forceps. Biopsies were placed immediately in ice cold medium and transferred to the laboratory where they were weighed and transferred to freezing medium (10% DMSO + fetal bovine serum) and cooled slowly overnight. Samples were shipped from Kenya to the United States in liquid nitrogen dry shippers, collagenase digested into single cell suspensions, and sorted on a FACS aria to enrich for a population of live, CD45+ cells. Sorted samples (N=17 with BV and N=13 without BV) were immediately lysed and processed using the Takara SMART-Seq v4 Ultra Low Input RNA Kit for sequencing. Library generation was completed using the manufacturer's instructions and sequencing performed on an Illumina NextSeq500 platform using Illumina NextSeq 500/550 High Output v2.5 (150 Cycles). The FASTQ files generated and BV status (N=9 with BV and N=12 without BV that passed quality control) of the donor sample will be available on dbGaP. We identified over 1000 differentially expressed genes between women with versus without BV at an adjusted p-value of <0.05 (listed under "gene names" below).

  Study phs002329

• Parallel Genomic Alterations of Antigen and Payload Targets Mediate Polyclonal Acquired Clinical Resistance to the Antibody Drug Conjugate Sacituzumab Govitecan

  Sacituzumab Govitecan (SG), a novel antibody drug conjugate (ADC) incorporating the anti-TROP2 antibody hRS7 conjugated to a topoisomerase-1 inhibitor (SN-38) payload, is the first ADC to be approved for advanced triple negative breast cancer (TNBC). However, mechanisms governing therapeutic resistance to SG are not known. We sought to identify mechanisms of de novo and acquired resistance to SG through unbiased whole-exome sequencing (WES) and RNA sequencing analysis of pre-treatment and multi-site post-progression (autopsy) tumor specimens. We examined three metastatic TNBC cases exhibiting (1) de novo progression, (2) stable disease, and (3) a deep response followed by progression, then mapped the temporal and spatial genomic evolution of acquired resistance in the responding patient. We then conducted additional pre-clinical experiments to validate the observed resistance mechanisms. TROP2 RNA and gene copy number were associated with de novo resistance, as case (1) was found to have absent TROP2 expression in all specimens, case (2) expressed TROP2, while case (3) exhibited both expression and focal genomic amplification of the TACSTD2/TROP2 locus. The genomic phylogeny tree inferred from case (3) post-progression specimens revealed one branch harboring an acquired canonical E418K resistance mutation in TOP1 (encoding topoisomerase 1) and a subsequent sub-clonal TOP1 inactivating frameshift mutation, while a distinct branch exhibited acquisition of a novel T256R missense mutation in TACSTD2 (encoding TROP2). Both the TOP1- and TACSTD2-mutant clones seeded multiple distinct metastatic sites. Through reconstitution experiments in TROP2 negative cells we found that TROP2 T256R is a stable protein with defective cell membrane localization and reduced cell surface binding by RS7 compared to wild-type TROP2. Collectively, these findings underscore TROP2 as a determinant of initial response to SG, and they reveal parallel and mutually exclusive polyclonal molecular mechanisms of acquired resistance involving the direct antibody target and drug payload target in distinct metastatic subclones of a single patient. While further research is needed to extend these novel findings, this study highlights the specificity of SG and illustrates how identifying such mechanisms will inform rational therapeutic strategies to overcome ADC resistance.&nbsp;&nbsp;

  Study phs002555

• 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

• Epigenomics of Human CD8 T cell Differentiation and Aging

  The efficacy of the adaptive immune response declines dramatically with age, but the cell-intrinsic mechanisms driving the changes characteristic of immune aging in humans remain poorly understood. One hallmark of immune aging is the loss of self-renewing naive cells and the accumulation of differentiated but dysfunctional cells within the CD8 T cell compartment. Using ATAC-seq, we first inferred the transcription factor binding activities that maintain the naive and central and effector memory CD8 T cell states in young adults. Integrating our results with RNA-seq, we determined that BATF, ETS1, Eomes, and Sp1 govern transcription networks associated with specific CD8 T cell subset properties, including activation and proliferative potential. Extending our analysis to aged humans, we found that the differences between memory and naive CD8 T cells were largely preserved across age, but that naive and central memory cells from older individuals exhibited a shift toward a more differentiated pattern of chromatin openness. Additionally, aged naive cells displayed a loss in chromatin openness at gene promoters, a phenomenon that appears to be due largely to a loss in binding by NRF1, leading to a marked drop-off in the ability of the naive cell to initiate transcription of mitochondrial genes. Our findings identify BATF- and NRF1-driven gene regulation as targets for delaying CD8 T cell aging and restoring T cell function.

  Study phs001187

• Human CD4 Memory T Cell Activation Time Course

  We performed a high-resolution time course RNA-seq experiment on CD4+ memory T cells during activation with anti-CD3/CD28 beads. We performed this in 24 adult individuals of European ancestry with no autoimmune diseases. We used 1 million cells as starting material, per time point (up to 8 time points per individual from 0 to 72 hours). We identified widespread dynamic allele-specific expression (dynASE) across the genome, i.e. ASE that changes with time after activation, which reflects cell state-dependent genetic regulatory effects. DynASE events were enriched in risk loci for multiple autoimmune diseases. This highlights the importance of studying multiple cellular states in order to elucidate the mechanisms of autoimmune disease. It also provides a successful framework to do it in other cellular systems relevant for other traits. We also showed an unprecedented amount of regulatory variation for HLA genes. HLA-DQB1 has 3 distinct regulatory haplotypes with expression profiles that could come from 3 different genes. We performed ATAC-seq in replicate for CD4+ memory T cells stimulated for 72 hours from another healthy individual. These data, together with genetic fine-mapping, CRISPR editing, EMSA and luciferase assays, helped us find a SNP in the promoter of HLA-DQB1 that causes at least part of the large cell state-dependent change in expression, at the RNA and protein levels.

  Study phs002259

• 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

• Human Pancreatic Islet RNAseq - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (&lt;73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004042

• Phenotyping data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (&amp;lt;73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004056

• GWAS data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (&lt;73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004044

• An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors

  We undertook an empirical approach to dissect the tumor microenvironment by developing a bulk tumor dissection algorithm, DisHet, and leveraging RNA-Seq data of tumorgrafts (patient-derived tumors implanted in mice), in which only the tumor cell component expands. We found that approximately 65% of previously defined immune signature genes are not abundantly expressed in the Renal Cell Carcinoma (RCC) microenvironment, and we identified more than two times as many novel immune/stromal transcripts. By using refined immune/stroma-specific genes and genomics, electronic medical record data and imaging data of 1084 RCC patients, we discovered a highly-inflamed pan-RCC subtype enriched for Treg cells, NK cells, Th1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. This inflamed subtype (IS) is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. Interestingly, IS is correlated with systemic manifestations of inflammation in patients such as thrombocytosis and anemia, whose pathogenesis were poorly defined and have been predictors of poor prognosis in RCC. Indeed, we discovered that IS is a strong predictor of poor survival. Lastly, our analyses show that tumor cells may drive stromal immune response. Overall, these data provide a missing link between the tumor cells, the tumor microenvironment, and systemic factors.

  Study EGAS00001002786

• 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

• Control of Focal Adhesion Kinase Activation by RUNX1-regulated miRNAs in high-risk AML

  We recently described a 16-gene expression signature for improved risk stratification of acute myeloid leukemia (AML) patients called the AML Prognostic Score (APS). A subset of APS high-risk AML patients showed increased levels of focal adhesion kinase (FAK), encoded by the Protein Tyrosine Kinase 2 (PTK2) gene, which was correlated with RUNX1 mutations. RUNX1 mutant cells are more sensitive to PTK2 inhibitors. As we were not able to detect RUNX1 binding sites in the PTK2 promoter, we hypothesized that RUNX1 might regulate micro(mi)RNAs that repress PTK2, such that loss-of-function RUNX1 mutations would result in reduced miRNA expression and derepression of PTK2. Examination of paired RNA-seq and miRNA-seq data from 301 AML cases revealed two miRNAs that positively correlated with RUNX1 expression, contained RUNX1 binding sites in their promoters and were predicted to target PTK2. We show that the hsa-let7a-2-3p and hsa-miR-135a-5p promoters are regulated by RUNX1, and that PTK2 is a direct target of both miRNAs. Even in the absence of RUNX1 mutations, hsa-let7a-2-3p and hsa-miR-135a-5p regulate PTK2 expression, and reduced expression of these two miRNAs sensitizes AML cells to PTK2 inhibition. These data explain how RUNX1 regulates PTK2, and identify potential miRNA biomarkers for targeting AML with PTK2 inhibitors.

  Study EGAS00001006491

• RNA-seq of high grade serous ovarian tumours

  Somatic RNA for 40 samples matched to the WGS was extracted using the Qiagen Qiasymphony RNA protcol (cat no 931636). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including DNase digestion). The resulting RNA the underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as measures of RNA purity. A260/280 had to be 2.0 and A260/230 had to be 2.0-2.2. Then RNA was quantified using LifeTechnologies Qubit RNA BR kit (cat no Q10210). RNAseq was carried out by the Edinburgh Clinical Research Facility on an Illumina NExtSeq500. Total RNA samples were assessed on the Agilent Bioanalyser (Agilent Technologies, #G2939AA) with the RNA 6000 Nano Kit (#5067-1512) for quality and integrity of total RNA, and then quantified using the Qubit 2.0 Fluorometer (Thermo Fisher Scientific Inc, #Q32866) and the Qubit RNA HS assay kit (#Q32855). Libraries were prepared from total-RNA sample using the NEBNext Ultra 2 Directional RNA library prep kit for Illumina (#E7760S) with the NEBNext rRNA Depletion kit (#E6310) according to the provided protocol. 400ng of totalRNA was then added to the ribosomal RNA (rRNA) depletion reaction using the NEBNext rRNA depletion kit (Human/mouse/rat) (#E6310). This step uses specific probes that bind to the rRNA in order to cleave it. rRNA-depleted RNA was then DNase treated and purified using Agencourt RNAClean XP beads (Beckman Coulter Inc, #66514). RNA was then fragmented using random primers before undergoing first strand and second strand synthesis to create cDNA. cDNA was end repaired before ligation of sequencing adapters, and libraries were enriched by PCR using the NEBNext Multiplex oligos for Illumina set 1 and 2 (#E7500). Final libraries had an average peak size of 271bp. Libraries were quantified by fluorometry using the Qubit dsDNA HS assay and assessed for quality and fragment size using the Agilent Bioanalyser with the DNA HS Kit (#5067-4626). Sequencing was performed using the NextSeq 500/550 High-Output v2 (150 cycle) Kit (# FC- 404-2002) on the NextSeq 550 platform (Illumina Inc, #SY-415-1002). Libraries were combined in an equimolar pool based on the library quantification results and run across 5 High-Output Flow Cell v2.5.

  Dataset EGAD00001009048

• RNA-seq of high grade serous ovarian tumours

  Somatic RNA for 37 samples was extracted using the Qiagen Qiasymphony RNA protcol (cat no 931636). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including DNase digestion). The resulting RNA the underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as measures of RNA purity. A260/280 had to be 2.0 and A260/230 had to be 2.0-2.2. Then RNA was quantified using LifeTechnologies Qubit RNA BR kit (cat no Q10210). RNAseq was carried out by the Edinburgh Clinical Research Facility on an Illumina NExtSeq500. Total RNA samples were assessed on the Agilent Bioanalyser (Agilent Technologies, #G2939AA) with the RNA 6000 Nano Kit (#5067-1512) for quality and integrity of total RNA, and then quantified using the Qubit 2.0 Fluorometer (Thermo Fisher Scientific Inc, #Q32866) and the Qubit RNA HS assay kit (#Q32855). Libraries were prepared from total-RNA sample using the NEBNext Ultra 2 Directional RNA library prep kit for Illumina (#E7760S) with the NEBNext rRNA Depletion kit (#E6310) according to the provided protocol. 400ng of totalRNA was then added to the ribosomal RNA (rRNA) depletion reaction using the NEBNext rRNA depletion kit (Human/mouse/rat) (#E6310). This step uses specific probes that bind to the rRNA in order to cleave it. rRNA-depleted RNA was then DNase treated and purified using Agencourt RNAClean XP beads (Beckman Coulter Inc, #66514). RNA was then fragmented using random primers before undergoing first strand and second strand synthesis to create cDNA. cDNA was end repaired before ligation of sequencing adapters, and libraries were enriched by PCR using the NEBNext Multiplex oligos for Illumina set 1 and 2 (#E7500). Final libraries had an average peak size of 271bp. Libraries were quantified by fluorometry using the Qubit dsDNA HS assay and assessed for quality and fragment size using the Agilent Bioanalyser with the DNA HS Kit (#5067-4626). Sequencing was performed using the NextSeq 500/550 High-Output v2 (150 cycle) Kit (# FC- 404-2002) on the NextSeq 550 platform (Illumina Inc, #SY-415-1002). Libraries were combined in an equimolar pool based on the library quantification results and run across 5 High-Output Flow Cell v2.5.

  Dataset EGAD00001010139

• FASTQ files of total RNA-Seq data of POPS control samples

  Placental biopsies (n = 64 female placentas, n = 67 male placentas) were selected from healthy pregnancies from the POPs cohort. These patients had no evidence of hypertension at booking and during pregnancy, did not experience pre-eclampsia, Hemolysis, Elevated Liver enzymes, and Low Platelets (HELLP) syndrome, gestational diabetes, or diabetes mellitus type I or type II and other obstetric complications. They delivered live babies with a birth weight percentile in the normal range (20-80th percentile), with no evidence of slowing in fetal growth trajectory. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003457

• Type I Interferon Exacerbates Mycobacterium Tuberculosis Induced Human Macrophage Death

  In this study, we tested whether type I interferons (IFN-I), namely IFN-alpha and IFN-beta, may promote Mycobacterium tuberculosis (Mtb) pathogenesis by exacerbating cell death of human monocyte-derived macrophages (MDM) infected with Mtb. We differentiated the MDM under two different published methods (see Document Method.pdf for Table 1 of Methods). To determine if MDM infected with Mtb upregulates IFN-I transcripts, signaling factors, or IFN-I stimulated genes (ISGs), the MDM were infected with Mtb or treated with an activator of stimulator of IFN genes (STING), diamidobenzimidazole (di-ABZI). STING activation acted as a positive control for IFN-I induction and ISG upregulation. Indeed, di-ABZI treated MDM differentiated under both methods upregulated ISGs 24 hpi. At 24 hour post treatment, MDM differentiated under one method, but not under another method, upregulated ISGs, such as ISG15, IFITM3, and CXCL10, upon Mtb infection. MDM differentiated under both methods upregulated ISGs upon di-ABZI treatment.

  Study phs003607

• Exome and RNA sequencing data from 32 ocular and extraocular sebaceous carcinomas

  Sebaceous carcinomas (SeC) are cutaneous malignancies that, in rare cases, metastasize and prove fatal. Here we report whole exome sequencing on 32 SeC, revealing distinct mutational classes that explain both cancer ontogeny and clinical course. A UV-damage signature predominated 10/32 samples, while 9 were instead defined by microsatellite instability (MSI) mutations. UV-damage SeC exhibited poorly differentiated, infiltrative histopathologycompared to MSI signature SeC (p = 0.003), features previously associated with dissemination. Strikingly, UV-damage SeC transcriptomes and anatomic distributionclosely resembling those of cutaneous squamous cell carcinomas (SCC), implicating sun-exposed keratinocytes as a cell of origin. Like SCC, this UV-damage subclass harbors a high somatic mutation burden with &gt;50 mutations/Mb, predicting immunotherapeutic response. In contrast, ocular SeC acquire far fewer mutations without a dominant signature, but show frequent truncating mutations in the ZNF750 epidermal differentiation regulator. Our data exemplify how different mutational processes convergently drive histopathologically related but clinically distinct cancers.

  Study EGAS00001002869

• 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.

  Dataset EGAD00001006309

• Integrated Transcriptomic and Regulatory RNA Profiling Reflects Complex Pathophysiology and Uncovers a Conserved Gene Signature in End Stage Heart Failure RNA-Seq data

  Heart failure (HF) is a complex syndrome. Despite availability of multiple treatment options, the mortality remains high and the quality of life poor. Better understanding of the underlying pathophysiological processes can lead to development of novel therapies. Multiple comparative transcriptomics studies, which revealed gene level changes in the key pathophysiological pathways in failing hearts, point towards heterogeneity from interplay of disease stage, etiologies and ethnicity. Transcriptomic characterization of HF in patients from different ethnicities can potentially help in understanding the heterogeneity imparted by various factors and the core elements in heart failure.

  Study EGAS50000000810

• Whole genome bisufite sequencing of smoking and non-smoking mother-child pairsBisufite sequencing, RNA-seq and ChIP-Seq data of whole blood samples from smoking and non-smoking mothers and their children at gestation/birth and follow-up years.

  Whole-blood samples of 16 mother-child pairs with differing maternal smoking behaviour were bisulfite treated and sequenced with deep coverage to identify smoking-associated differentially methylated regions.RNA-sequencing and ChIP-sequencing of 4 chromatin marks complete the data set and allow for integrative analyses, functional region annotation and to study the impact of epigenetic changes on gene expression.Longitudinal data spanning several years provides the means to investigate the stability of observed differences of all data types.

  Study EGAS00001000455

• BLUEPRINT: Epigenetic programming during monocyte to macrophage differentiation and trained innate immunity

  Monocyte differentiation into macrophages represents a cornerstone process for host defense. Concomitantly, immunological imprinting of either tolerance or trained immunity determines the functional fate of macrophages and susceptibility to secondary infections. Transcriptomes (RNA-Seq) and epigenomes (ChIP-Seq H3K4me1,H3K4me3,H3K27ac) in four primary cell types: monocytes, in vitro differentiated naive, tolerized and trained macrophages were characterized. Inflammatory and metabolic pathways were modulated in macrophages, including decreased inflammasome activation, and pathways functionally implicated in trained immunity were identified. Strikingly, B-glucan training elicits an exclusive epigenetic signature, revealing a complex network of enhancers and promoters. Analysis of transcription factor motifs in DNase I hypersensitive sites at cell-type specific epigenetic loci unveiled differentiation and treatment specific repertoires. Altogether, this study provides a resource to understand the epigenetic changes that underlie innate immunity in humans.

  Dataset EGAD00001001011

• Transcriptome analysis of leiomyosarcoma

  We performed transcriptome analysis of leiomyosarcoma

  Study EGAS50000000594

• Gene Expression and Regulatory Networks in Human Leukocytes - Immunological Variation Consortium

  Profiling of gene expression with microarrays holds great potential for human health, for illuminating disease pathways or providing biomarkers to monitor disease or its resolution. Using high-throughput approaches for genotyping, immunophenotyping and gene expression analysis, the project will examine the basis for the control of gene expression in human immune cells, and how it is influenced by natural genetic variation or aging. In practice, the project will combine several cross-informative approaches: 1) Building on the established sample/data pipelines and robust protocols of the Immunological Genome (ImmGen) project, and on established cohorts of ethnically diverse healthy volunteers at hand, microarray techniques will be used to generate whole-genome expression profiles from purified naive CD4+ lymphocytes and monocytes from 600 healthy volunteers. A dense genetic map will be established for all donors. The results will elucidate how variation in the human genome affects the expression of immune genes, of key importance in understanding gene variants that bring susceptibility to immune or inflammatory disease. Computational analysis of these rich data will allow the reconstruction of regulatory connections between genes, helping to establish general modules and those specific of a given immune cell type. These data will be complemented by an orthogonal data group, generated from a restricted subset of 10 individuals, in which we will profile a larger set of 28 carefully delineated cell populations that exist in human blood. This work will also benefit from powerful interspecies comparison with similar experiments being performed in mice by ImmGen. 2) In addition, RNA from the same set of 28 defined cell populations will be probed with microarrays that explore other aspects of the transcriptome: i) microRNAs and other non-coding RNAs; ii) exon or splice junction arrays that will establish a map of differential splicing in human blood leukocyte. 3) The composition and reactivity of blood cells from the same donors will be established at the time of sample collection using high-throughput flow cytometry, correlating immune phenotypes with gene expression and genetic variation. 4) Genetic variability conditions the baseline levels of gene expression, but also the responsiveness to activating challenges. With samples from the same donors, NanoString technology for transcript counting will be used to analyze the transcriptional response of defined gene sets, representing response signatures of T or dendritic cells, for a fine-grained dissection of responses to different triggers (different bacterial ligands for dendritic cells, different cytokine environment for T cells). In keeping with the resource aspect of this project, all data and interpretations will be made publicly available rapidly upon curation, allowing public querying and browsing of the data, by using and evolving the existing web architectures of the ImmGen project, of the Broad Institute, and of international repositories. RELEVANCE: Exploration with DNA chips of the genome&#39;s expression microarrays holds great potential for human health, to better understand disease and to serve as diagnostic tools. Using a combination of high-throughput genomic techniques and computational biology, we will perform a broad exploration of gene expression in human blood cells across groups of African-American, Asian and European ancestry, asking how these profiles are affected by genetic variation or by age. These results will provide an invaluable reference benchmark for the interpretation of genetic and immunological studies. Additionally, over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, in a second study we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4+ T cells in up to 105 healthy donors. We found regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution. 15% of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (ATAC-QTLs). ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression, and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression. In another study we analyzed RNA-seq profiles from primary CD4+ T cells of two healthy donors at six time points from 0 to 72 hours after stimulation with bead-bound antibodies against CD3 and CD28.

  Study phs000815

• Multimodal Immune Profiling to Determine Mechanisms of COVID-19 Clinical Trajectory in Uganda

  We conducted a prospective cohort study of adults (RESERVE-U-C19) admitted to Entebbe Regional Referral Hospital (ERRH), a national COVID-19 referral center, across three viral variant-driven phases of the pandemic in Uganda. In this cohort, which reflects the entire spectrum of SARS-CoV-2 infection, we integrated clinical data with serum proteomics (N=306) and whole-blood transcriptomics (N=100) to determine a range of host responses associated with COVID-19 severity, SARS-CoV-2/HIV co-infection, and circulating viral variants. Beyond the prognostic importance of myeloid cell-driven immune activation and lymphopenia, we show that multifaceted impairment of host protein synthesis and extensive redox imbalance define core metabolic signatures of severe COVID-19, with central roles for IL-7, IL-15, and lymphotoxin-&#945; in COVID-19 respiratory failure. While prognostic signatures were generally consistent in SARS-CoV-2/HIV-coinfection, type I interferon responses uniquely scaled with COVID-19 severity in persons living with HIV. Throughout the pandemic, COVID-19 severity peaked during phases dominated by A.23/A.23.1 and Delta B.1.617.2/AY variants, with Delta phase COVID-19 distinguished by exaggerated pro-inflammatory myeloid cell and inflammasome activation, NK and CD8+ T-cell depletion, impaired host protein synthesis, and upregulation of viral cell entry and trafficking pathways. Data from this study available through dbGaP include raw RNA-Seq data from whole-blood of adult participants (N=100).

  Study phs003246

• MAITS in HCC

  Hepatocellular Carcinoma (HCC) is a leading cause of cancer-related death and can be considered a prototype of inflammation-derived cancer arising from chronic liver injury. The cell composition of the HCC tumor immune microenvironment (TiME) has a major impact on cancer biology as the TiME can have divergent capacities on tumor initiation, progress, and response to therapy. Recent development of multi-omics and single-cell technologies help us to comprehensively quantify the cellular heterogeneity and spatial organization of the TiME and to further our understanding of antitumor immunity. We investigated the cellular composition of liver cancer patient samples (n=8) using single-cell RNA-seq. For this purpose, mononuclear cells were isolated from human liver cancer samples. Three locations within the tumor-bearing liver were used: adjacent liver, rim and tumor core. Then, cells were FACS sorted (either CD45+ cells or MAIT cells) and subjected to 10X Chromium-based scRNA-seq. We investigated immune cell changes and the cellular heterogeneity of the HCC TiME with a focus on MAIT cells. Annotated H5 files are provided. Clinical metadata including TMN stage, sex, gender, ethnicity, pretreatment, and histopathological reports are available for all patient samples. Further details on the study can be obtained in our paper once it&#8217;s published.A complementary dataset derived from ultrahighplex CO-Detection by indEXing (CODEX) from paired patient samples can be accessed through The Cancer Imaging Archives (TCIA) under DOI: https://doi.org/10.7937/bh0r-y074.

  Study phs003279

• Genomic Profiling of Peripheral T-cell Lymphoma

  Peripheral T-cell lymphomas (PTCLs) are malignant and highly aggressive hematologic tumors arising from mature post-thymic T cells. The diagnosis of PTCL includes diverse lymphoma subgroups, altogether accounting for about 15% of all non-Hodgkin lymphomas. Despite recent advances in uncovering the mutational landscape of PTCLs, the prognoses of most PTCLs remain poor. To overcome the limits of exome sequencing and identify new genomic abnormalities in PTCL, we performed RNAseq and whole genome sequencing. From our analyses, we identified a recurrent in-frame deletion (VAV1 &#916;778-786) generated by a focal deletion-driven alternative splicing mechanism as well as novel VAV1 gene fusions (VAV1-THAP4, VAV1-MYO1F, and VAV1-S100A7). These genetic lesions result in increased activation of VAV1 catalytic-dependent (MAPK, JNK) and non-catalytic-dependent (nuclear factor of activated T cells, NFAT) VAV1 effector pathways. These results support an oncogenic role for VAV1 signaling in the pathogenesis of PTCL (Abate et al., 2017). Recently, we have also identified a novel chimeric oncogene, FYN-TRAF3IP2, in the RNA-seq and whole genome sequencing analyses. FYN-TRAF3IP2 leads to aberrant NF-&#954;&#946; activity downstream of T-cell receptor signaling and induces NF-&#954;&#946;-driven T-cell transformation when expressed in hematopoietic progenitors in mice. This novel gene fusion supports the oncogenic driver role of NF-&#954;&#946; signaling in the pathogenesis of PTCL (Kim et al., in review).

  Study phs001962

• The Library of Integrated Network-Based Cellular Signatures (LINCS) - NeuroLINCS

  The NeuroLINCS Center is part of the NIH Common Fund&#39;s Library of Integrated Network-based Cellular Signatures (LINCS) program, which aims to characterize how a variety of human cells, tissues and the entire organism respond to perturbations by drugs and other molecular factors. As Part of the LINCS program, the NeuroLINCS study concentrates on human brain cells, which are far less understood than other cells in the body. Our initial focus is to produce diseased motor neurons from patients by utilizing high-quality induced pluripotent stem cell (iPSC) lines from Amyotrophic Lateral Sclerosis (ALS) and Spinal Muscular Atrophy (SMA) patients in addition to unaffected normal healthy controls. Using state-of-the-art OMICS methods (genomics, epigenomics, transcriptomics, and proteomics), we intend to create a wealth of cellular data that is patient-specific in the context of their baseline genetic perturbations and in the presence of other genetic and environmental perturbagens (e.g. endoplasmic reticulum stress). The primary data will be used to build cell signatures that convey the key features that distinguish the state of a cell and determine its behavior. Ultimately, the analysis of these datasets will lead to the identification of a network of unique signatures relevant to each of these motor neuron diseases. The datasets represented in this study are generated from assays interrogating RNA expression (RNA-seq), chromatin accessibility (ATAC-seq) and whole genome sequencing.

  Study phs001231

• Detection and Targeting of Splicing Deregulation in Pediatric Acute Myeloid Leukemia Stem Cells

  This study evaluates the gene expression and splicing changes between pediatric patients with and without Acute Myeloid Leukemia (AML) and pediatric and adult patients with AML. We obtained bone marrow or peripheral blood from patients and isolated CD34+ stem and/or progenitor cells from pediatric patients with AML (10) or without AML (6) and adult patients with AML (5). In parallel, a targeted analysis for known, deleterious, genetic mutations in the KRAS, NRAS, and KMT2C genes identified 2, 5, and 6 AML patients with alterations, respectively. Additionally, Tapestri analysis further identified a mutation in the intronic region of the SF3B1 gene in multiple pediatric AML patients. RNA-Seq was performed on the purified cell populations at The Scripps Research Institute Next Generation Core on an Illumina NextSeq 500 sequencer with 150bp paired-end reads.These datasets revealed that, similar to adult AML, pediatric AML samples were enriched for Leukemia Stem Cells (LSC) markers. A deeper analysis showed widespread dysregulation of splicing in pediatric AML samples compared to age-matched non-leukemic samples. In addition, the splicing regulator RBFOX2 was downregulated with an accompanying increase in specific CD47 splice isoforms in the pediatric AML samples. Thus, the increased presence of LSCs in pediatric AML is possibly driven by global splicing deregulation, and the treatment with splicing modulators like Rebecsinib has a potential therapeutic value for eradicating LSCs in these patients.

  Study phs003196

• Region-Specific Neural Stem Cell Lineages Revealed by Single-Cell RNA-Sequences from Human Embryonic Stem Cells

  During development of the human brain, multiple cell types with diverse regional identities are generated. Here we report a system to generate early human brain forebrain and mid/hindbrain cell types from human embryonic stem cells (hESCs), and infer and experimentally confirm a lineage tree for the generation of these types based on single-cell RNA-Seq analysis. We engineered SOX2Cit/+ and DCXCit/Y hESC lines to target progenitors and neurons throughout neural differentiation for single-cell transcriptomic profiling, then identified discrete cell types consisting of both rostral (cortical) and caudal (mid/hindbrain) identities. Direct comparison of the cell types were made to primary tissues using gene expression atlases and fetal human brain single-cell gene expression data, and this established that the cell types resembled early human brain cell types, including preplate cells. From the single-cell transcriptomic data a Bayesian algorithm generated a unified lineage tree, and predicted novel regulatory transcription factors. The lineage tree highlighted a prominent bifurcation between cortical and mid/hindbrain cell types, confirmed by clonal analysis experiments. We demonstrated that cell types from either branch could preferentially be generated by manipulation of the canonical Wnt/beta-catenin pathway. In summary, we present an experimentally validated lineage tree that encompasses multiple brain regions, and our work sheds light on the molecular regulation of region-specific neural lineages during human brain development.

  Study phs001205