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

in 115.90 milliseconds.

• Single cell RNAseq of PBMC from RCC patients

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

  Study EGAS00001004451

• The endometrial transcription landscape of MRKH syndrome

  The Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome (OMIM 277000) is characterized by agenesis of the uterus and upper part of the vagina in females with normal ovarian function. While genetic causes have been identified for a small subset of patients and epigenetic mechanisms presumably contribute to the pathogenic unfolding, too, the etiology of the syndrome has remained largely enigmatic. A comprehensive understanding of gene activity in the context of the disease is crucial to identify etiological components and their potential interplay. So far, this understanding is lacking, primarily due to the scarcity of samples and suitable tissue. In order to close this gap, we profiled endometrial tissue of uterus rudiments in a large cohort of MRKH patients using RNA-seq and thereby provide a genome-wide view on the altered transcription landscape of the MRKH syndrome.

  Study EGAS00001004601

• HBCC Postmortem Psychiatric Molecular Studies

  This postmortem study examines molecular, genetic and epigenetic signatures in the brains of hundreds of subjects with or without mental disorders conducted by the DIRP NIMH Human Brain Collection Core (HBCC). The brain tissues are obtained under protocols approved by the CNS IRB (NCT00001260), with the permission of the next-of-kin (NOK) through the Offices of the Chief Medical Examiners (MEOs) in the District of Columbia, Northern Virginia and Central Virginia. Additional samples were obtained from the University of Maryland Brain and Tissue Bank (contracts NO1-HD-4-3368 and NO1-HD-4-3383) (http://www.medschool.umaryland.edu/btbank/ and the Stanley Medical Research Institute: http://www.stanleyresearch.org/brain-research/). Clinical characterization, neuropathological screening, toxicological analyses, and dissections of various brain regions were performed as previously described (Lipska et al. 2006; PMID: 16997002). All patients met DSM-IV criteria for a lifetime Axis I diagnosis of psychiatric disorders including schizophrenia or schizoaffective disorder, bipolar disorder and major depression. Controls had no history of psychiatric diagnoses or addictions. SNP array: Array-based genotyping was performed on most samples published in this collection. The number of SNPs assayed via Illumina chips varied between 650,000 and 5 Million. Cerebellar tissue was generally used for genotyping studies. # Diagnosis SNP Array 1 Anxiety Disorder 1 2 Autism Spectrum Disorder 13 3 Bipolar Disorder 114 4 Control 387 5 Eating Disorder (ED) 2 6 Major Depressive Disorder (MDD) 186 7 Obsessive Compulsive Disorder (OCD) 5 8 Post-Traumatic Stress Disorder (PTSD) 0 9 Schizophrenia 220 10 Other 7 11 Tic Disorder 3 12 Undetermined 1 13 Williams Syndrome 2 Table: Numbers of samples in each diagnostic category. DNA extraction: 45-80 mg of cerebellar tissue was pulverized for DNA extractions. The QIAamp DNA mini Kit (Qiagen) method was employed for tissue DNA extraction. The tissue was initially lysed using Tissue Lyser (Qiagen) and extractions were accomplished according to manufacturer's protocol. The DNA was captured in 500uL elution buffer. The concentrations were measured using Thermo Scientific's NanoDrop 1000/NanoDrop ONE. The mean yield was 128.85 uG (+/- 79.48), the mean ratio of 260/280 was 1.87 (+/- 0.105), and the mean ratio of 260/230 was 2.48 (+/-1.75). Genotyping methods: Three types of Illumina Beadarray chips were used: HumanHap650Y, Human1M-Duo, and HumanOmni5M-Quad (San Diego, California). The genotyping was done according to the manufacturer's protocol (Illumina Proprietary, Catalog # WG-901-5003, Part # 15025910 Rev.A, June 2011). Approximately, 400ng DNA was used and each DNA sample was QC tested for 260/280 ratio by nanodrop and DNA band intactness on 2% agarose gel. Briefly, the samples were whole-genome amplified, fragmented, precipitated and resuspended in appropriate hybridization buffer. Denatured samples were hybridized on prepared Bead Array Chips. After hybridization, the Bead Chip oligonucleotides were extended by a single fluorescent labeled base, which was detected by fluorescence imaging with an Illumina Bead Array Reader, iScan. Normalized bead intensity data obtained for each sample were loaded into the Illumina Genome Studio (Illumina, v.2.0.3) with cluster position files provided by Illumina, and fluorescence intensities were converted into SNP genotypes. Microarray: We generated RNA expression data using array technology for psychiatric subjects compared to non-psychiatric subjects as controls. We used tissues from three different brain regions i.e. hippocampus, dorsolateral prefrontal cortex (DLPFC), and dura mater for a large cohort of individuals (total number 552 subjects for hippocampus, 800 for DLPFC and 146 for dura). Total RNA was extracted from ~100 mg of tissue using the RNeasy kit (Qiagen) according to the manufacturer's protocol. RNA quality and quantity were examined using the Bioanalyzer (Agilent, Inc) and NanoDrop (Thermo Scientific, Inc), respectively. Samples with RNA integrity number (RIN) # Diagnosis DLPFC Hippo Dura 1 Anxiety Disorder 1 0 0 2 Autism Spectrum Disorder 14 6 0 3 Bipolar Disorder 90 49 0 4 Control 336 270 75 5 Eating Disorder (ED) 2 1 0 6 Major Depressive Disorder (MDD) 144 87 0 7 Obsessive Compulsive Disorder (OCD) 5 3 0 8 Post-Traumatic Stress Disorder (PTSD) 6 0 0 9 Schizophrenia 192 125 71 10 Other 5 6 0 11 Tic Disorder 3 3 0 12 Undetermined 1 1 0 13 Williams Syndrome 2 1 0 Table: Numbers of samples in each diagnostic category. RNA-Seq of Dorso-lateral prefrontal cortex: All brains were collected and the dorsolateral prefrontal cortical (DLPFC) samples dissected at the HBCC, DIRP, NIMH. Dorsolateral prefrontal cortex (DLPFC) specimens were dissected from right or left hemisphere of frozen coronal slabs. The study was funded by the DIRP, NIMH under contract (#HHSN 271201400099C) with Icahn School of Medicine at Mount Sinai,1106402 One Gustave L. Levy Place, Box 3500, New York NY 10029-6574. RNA extraction, library preparation and sequencing were performed under contract at Icahn School of Medicine. The Common Mind Consortium (CMC) provided project management support. RNA isolation: Total RNA from 468 HBCC samples was isolated from approximately 100 mg homogenized tissue from each sample by TRIzol/chloroform extraction and purification with the Qiagen RNeasy kit (Cat#74106) according to manufacturer's protocol. Samples were processed in randomized batches of 12. The order of extraction for schizophrenia, bipolar, and MDD disorders and control samples was assigned randomly with respect to diagnosis and all other sample characteristics. The mean total RNA yield was 24.2 ug (+/- 9.0). The RNA Integrity Number (RIN) was determined by 4200 Agilent TapeStation System. Samples with RIN DLPFC RNA-Seq quantified expression data are provided for 364 samples. Data were generated, QC'd, processed and quantified as follows: RNA library preparation and sequencing: All samples submitted to the New York Genome Center for RNAseq were prepared for sequencing in randomized batches of 94. The sequencing libraries were prepared using the KAPA Stranded RNAseq Kit with RiboErase (KAPA Biosystems). rRNA was depleted from 1ug of RNA using the KAPA RiboErase protocol that is integrated into the KAPA Stranded RNAseq Kit. The insert size and DNA concentration of the sequencing library was determined on Fragment Analyzer Automated CE System (Advanced Analytical) and Quant-iT PicoGreen (ThermoFisher) respectively. Schizophrenia Bipolar Control 89 65 210 Table: Numbers of samples in each diagnostic category. RNA-Seq of subgenual anterior cingulate cortex (sgACC): All the 200 post-mortem brain samples (61 controls; 39 bipolar disorder; 46 schizophrenia; 54 major depressive disorder) were collected by the HBCC, DIRP, NIMH. RNA Extraction and Quality Assessment: Tissue from sgACC was pulverized and stored at -80°C. Total RNA was extracted from 50-80 mg of the tissue using QIAGEN RNeasy Lipid Tissue Mini Kit (QIAGEN, Cat. # 74804) with DNase treatment (QIAGEN, Cat. # 79254). The RNA Integrity Number (RIN) for each sample was assessed with high-resolution capillary electrophoresis on the Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, California). The concentration of RNA and their 260/280 ratio (2.1+/- 0.032 SD) were determined with NanoDrop (Thermo Scientific). RNA sequencing: Stranded RNA-Seq libraries were constructed after rRNA depletion using Ribo-Zero GOLD (Illumina). RNA sequencing was performed at National Institute of Health Intramural Sequencing Center (NISC). Schizophrenia Bipolar Control MDD 46 39 61 54 Table: Numbers of samples in each diagnostic category. Whole Genome Sequencing: All brains were collected and dissected at the HBCC, DIRP, NIMH. This study generates whole genome sequencing data using sequencing of DNA in the dorsolateral prefrontal cortex (DLPFC), anterior cingulate cortex (ACC) or cerebellum of 443 individuals with schizophrenia, bipolar disorder and major depressive disorder and non-psychiatric controls. The study was funded by the DIRP, NIMH under contract (#HHSN 271201400099C) with Icahn School of Medicine at Mount Sinai,1106402 One Gustave L. Levy Place, Box 3500, New York NY 10029-6574. DNA extraction, library preparation and sequencing were performed under contract at Icahn School of Medicine. The Common Mind Consortium (CMC) provided project management support. All specimens were dissected from right or left hemisphere of frozen coronal slabs. DNA Library Preparation and Sequencing: All samples submitted to the New York Genome Center for WGS were prepared for sequencing in randomized batches of 95. The sequencing libraries were prepared using the Illumina PCR-free DNA sample preparation Kit. The insert size and DNA concentration of the sequencing library was determined on Fragment Analyzer Automated CE System (Advanced Analytical) and Quant-iT PicoGreen (ThermoFisher) respectively. A quantitative PCR assay (KAPA), with primers specific to the adapter sequence, was used to determine the yield and efficiency of the adaptor ligation process. Performed on the Illumina HiSeqX with 30X coverage. Schizophrenia Bipolar Control 115 78 230 Table: Numbers of samples in each diagnostic category. ChIP-Seq: All brains were collected and the dorsolateral prefrontal cortical (DLPFC) samples dissected at the HBCC, DIRP, NIMH. This study generates epigenetic data using sequencing of DNA after chromatin immunoprecipitation (ChIP-Seq) for marks H3K4me3 and H3K27ac in the dorsolateral prefrontal cortex (DLPFC). Dorsolateral prefrontal cortex (DLPFC) specimens were dissected from right or left hemisphere of frozen coronal slabs. The study was funded by the DIRP, NIMH under contract (#HHSN 271201400099C) with Icahn School of Medicine at Mount Sinai,1106402 One Gustave L. Levy Place, Box 3500, New York NY 10029,6574. Chromatin precipitation, library preparation and sequencing were performed under contract at Icahn School of Medicine. The Common Mind Consortium (CMC) provided project management support. Chromatin immunoprecipitation (ChIP) assays for histone marks H3K4me3 and H3K27ac were carried out using Native ChIP. Micrococcal Nuclease (MNase) (Sigma, N3755) treatment was used to digest chromatin into mononucleosomes. The following antibodies were used for chromatin pull-down: anti-H3K4me3 (Cell Signaling, Cat# 9751BC, lot 7) and anti-H3K27ac (Active Motif, Cat# 39133, Lot # 31814008). Histone modification-enriched genomic DNA fragments were recovered using Protein A/G magnetic beads (Thermo Scientific, 88803-88938 or Millipore 16-663), and then washed, eluted, and treated with RNAse A and proteinase K. Final ChIP DNA products were isolated using phenol-chloroform extraction followed by ethanol precipitation. The efficiency of each ChIP assay was validated using Qubit concentration measurement and qPCR for positive (GRIN2B, DARPP32) and negative (HBB) control genomic regions. Only ChIP assays that passed quality control were further processed for library preparation and sequencing; this included ChIP DNA that was not detectable on Qubit but showed a good signal and expected enrichment patterns in qPCR. HISTONE_MARK H3K27ac H3K4me3 Input Bipolar 56 4 7 Control 158 11 24 Schizophrenia 79 11 12 Table: Numbers of individuals in each assay grouped by histone mark or input.Long-Read Whole-Genome Sequencing (WGS) Cohort Description: Brain specimens were obtained from the Human Brain Collection Core (HBCC), part of the NIH NeuroBioBank. Samples were collected under protocols approved by the NIH CNS Institutional Review Board (IRB) (NCT03092687), with informed consent from next-of-kin (NOK). Collection was coordinated through the Offices of the Chief Medical Examiners (MEOs) in Washington, D.C., Northern Virginia, and Central Virginia. Clinical metadata and documentation are publicly available via the NIMH Data Archive (NDA) (Collection #3151) https://nda.nih.gov/edit_collection.html?id=3151 Eligibility Criteria No clinical diagnosis of major neuropsychiatric or neurodegenerative diseaseNo diagnosis of cognitive impairment during life All individuals were confirmed to be neurologically normal at time of deathDemographics Initial cohort size: 155 individuals Ancestry: All individuals self-identified as African or African-admixed Mean age at death: 44.2 years (range: 18–85 years) Sex distribution: 36.4% femaleSample Processing: Frozen frontal cortex tissue was dissected and processed according to the public protocol: https://www.protocols.io/view/processing-human-frontal-cortex-brain-tissue-for-p-kxygxzmmov8j/v2. High-molecular-weight DNA was extracted and libraries were prepared using the Oxford Nanopore Technologies (ONT) LSK-114 kit. Sequencing was performed using ONT PromethION flow cells (R10.4.1 chemistry) Data Processing and Quality Control: Basecalling: Conducted using Guppy v6.38 Read Alignment: Reads were aligned to the GRCh38 reference genome using minimap2 Sample Identity Verification: Sample identity was validated by comparing ONT-derived SNP calls with matched short-read WGS genotypes to ensure concordance and prevent sample swaps Variant Calling and Phasing: Reads were base-called with Guppy v6.38. Reads were aligned to GRCh38 using minimap2. We verified sample identity by cross-checking ONT SNV calls with the existing short-read WGS genotypes, confirming no sample switches. The napu pipeline (https://github.com/nanoporegenomics/napu_wf) produced; haplotype-resolved assemblies, joint small-variant (SNV/indel) calls, and multi-caller structural-variant sets, all reported on GRCh38 and phased where possible. Raw signal data were basecalled to obtain 5-methyl-cytosine (5mC) status; methylation tags were added to the phased BAM files. Genome-wide methylation summaries are provided in BED format.Dataset Filtering and Exclusions: All 155 samples underwent sequencing and SNP-based ancestry inference 8 samples were excluded due to ancestry inconsistent with African or African-admixed background 1 sample was excluded due to insufficient sequencing quality Final Sample Set: 146 high-quality samples from individuals of African or African-admixed ancestry were retained for downstream analyses See PMID: 39764002 for further analysis detailsDiagnosis#SamplesControl155Table: Diagnostic Summary.Note: The data derived from HBCC resources were removed from dbGAP and are now available in the NIMH Data Archive (NDA). They include genotypes, short read whole genome sequencing (WGS), epigenetics (DNA methylation, ChIP-seq for histones), RNA expression (qPCR, microarray, RNA-seq, single nucleus RNA-seq) of various brain regions in cases with schizophrenia, bipolar disorder, major depression, substance use disorders and normative controls. Please access our NDA collection (https://nda.nih.gov/edit_collection.html?id=3151) for further detail.

  Study phs000979

• 10x genomics single cell RNA-seq and ATAC-seq from 10 pre-menopausal patients fimbriae and ampullary part of the fallopian tube

  The fallopian tube, connecting the uterus with the ovary, is a dynamic organ that undergoes cyclical changes and is the site of several diseases, including serous cancer. Here, we use single-cell technologies to construct a comprehensive cell map of healthy pre-menopausal fallopian tubes, capturing the impact of the menstrual cycle and menopause on different fallopian tube cells at the molecular level. The comparative analysis between pre- and post-menopausal fallopian tubes reveals substantial shifts in cellular abundance and gene expression patterns, highlighting the physiological changes associated with menopause. Further investigations into menstrual cycle phases illuminate distinct molecular states in secretory epithelial cells caused by hormonal fluctuations. The markers we identified characterizing secretory epithelial cells provide a valuable tool for classifying ovarian cancer subtypes.

  Study EGAS50000000628

• RNA-seq data from 224 advanced prostate tumors generated by the West Coast Dream Team

  224 pairs of FASTQ files from metastatic Castration-Resistant Prostate Cancer (mCRPC) sequenced on HiSeq 4000 instruments. Patients were enrolled in the West Coast Dream Team study. Biopsies include various tissue sites including bone, soft tissue, and lymph node.

  Dataset EGAD00001008487

• RNA-seq data from 164 advanced prostate tumors generated by the West Coast Dream Team including 42 paired samles

  164 pairs of FASTQ files from metastatic Castration-Resistant Prostate Cancer (mCRPC) sequenced on HiSeq 4000 instruments. Patients were enrolled in the West Coast Dream Team study. Biopsies include various tissue sites including bone, soft tissue, and lymph node. 42 pairs prior to enzalutamide treatment and at progression from 21 patients are included.

  Dataset EGAD00001009065

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

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

  Study phs002756

• Utility of Capillary Blood in Gene Expression Studies

  The FAA Functional Genomics Team examined the feasibility of analyzing gene expression profiles using capillary blood collected by fingerstick, with the goal of developing practical fingerstick blood collections to expand blood collection capabilities and allow large-scale or field-based collections. We compared the RNA sequencing results obtained from two different capillary blood collection and RNA extraction methods with the results obtained from a standard venous blood collection and RNA extraction protocol. We also compared capillary blood processing methods, as well as fingerstick locations, within one capillary collection/extraction method. Within each method we sought to distinguish male from female participants to assess the capacity of each collection/extraction method for gene expression profiling. The study population included 40 participants aged 18-60 who were free from any acute illness or any illness that would preclude blood collection. We determined that capillary blood possesses large-scale significant differences from venous blood and that the two capillary methods produced similar gene expression profiles. We further determined that all methods were useful for detecting gene expression differences between assessed groups (male and female participants in this case). We conclude that capillary blood is suitable for gene expression profiling in quantities as low as 100 μL whole blood, but recommend that researchers should rely on a single blood collection method if they desire to compare results between individual studies. Raw sequencing files (.fastq) from this study are available in dbGaP.

  Study phs003496

• Cell Type-Specific and Disease-Associated eQTL in the Human Lung

  This study was conducted to investigate the genetic control of gene expression in a cell-type and context-specific manner in the human lung. We performed single-cell RNA-sequencing (scRNA-seq) of lung tissue from 66 Interstitial Lung Disease (ILD) cases and 48 unaffected donors. Employing a pseudo-bulk approach, we mapped expression quantitative trait loci (eQTL) across 38 cell types. The scRNA-seq data presented here include previously published (PMID 34262047) and unpublished samples. ILD tissue samples were obtained from lungs removed at the time of lung transplantation at either Vanderbilt University Medical Center (VUMC) or the Norton Thoracic Institute (NTI). Control tissue samples were obtained from lungs declined for organ donation either at the Donor Network of Arizona (DNA) or VUMC. The self-reported ethnicities of donors were 66.7% European, 9.6% African American, 17.5% not available (N/A), and 6.1% other. The smoking history includes 46.5% ever smoker, 29.8% never smoker, 23.7% N/A. For genotyping, flash-frozen tissue in DNA/RNA shield was homogenized via a bullet blender. Genomic DNA was extracted using Zymo Quick DNA/RNA microprep plus kit. Library preparation and low-pass Whole Genome Sequencing (WGS) were carried out at TGen or by Gencove Inc. At TGen, libraries were prepared using PCR-free Watchmaker kits (Watchmaker Genomics) with 200 ng input. Genomes were sequenced on NovaSeq at low coverage (typically 0.4–1x). The resulting sequence data were processed and imputed using Gencove's imputation platform.

  Study phs003521

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

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

  Study EGAS50000000973

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

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

  Study EGAS50000000135

• 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

• RNA Transcriptomic and DNA Methylation Landscape in FTLD-TDP and Controls

  Frontotemporal lobar degeneration (FTLD) is one of the leading causes of dementia in individuals younger than 65 years, with the aggregation of TDP-43 as one of the most common forms of FTLD. FTLD-TDP is clinically, genetically and pathologically heterogeneous, with GRN and C9orf72 as the most common genetic forms (although more than 50% of the cases are genetically unexplained), and A, B and C as the most common pathological subtypes. To investigate the molecular differences between the different forms of FTLD-TDP, we performed bulk RNA sequencing (RNAseq) and reduced representation bisulfite sequencing (RRBS) from different brain regions of FTLD-TDP patients and controls. The RNAseq dataset includes samples from the frontal cortex of 149 individuals, comprising 127 FTLD-TDP cases (genetically unexplained types A (27), B (20) and C (22), 24 GRN mutation carriers, and 34 C9orf72 expansion carriers) and 22 controls. The RRBS dataset includes samples from the frontal cortex of 158 individuals —119 FTLD-TDP cases (genetically unexplained types A (25), B (24), and C (25), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 39 controls; and from the cerebellum of 151 individuals —113 FTLD-TDP cases (genetically unexplained types A (24), B (21), and C (23), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 38 controls. In total, the study consists of data from 183 unique individuals, with some overlap across the different datasets.

  Study phs004075

• Relapse CHL study

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

  Study EGAS00001008222

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

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

  Study EGAS00001003950

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

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

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

  Study EGAS00001003452

• 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

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

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

  Study EGAS00001003975

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• Exome-seq, RNA-Seq, SNP array profiling of gastric tumor samples and cell lines.

  Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we performed exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumors and 32 cell lines. Meta-analysis of the exome data and previously published data sets revealed 24 significantly mutated genes in MSS tumors, and 16 in MSI tumors. More than half of the patients in our collection could potentially benefit from targeted therapies. We identified 55 splice-site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumors. The ZAK kinase isoform TV1 is preferentially upregulated in gastric tumors and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder, and breast cancers. Overexpression of this particular isoform activated multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibited proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

  Study EGAS00001000736

• atrial appendage RNA-seq in non-, paroxysmal and persistent atrial fibrillation

  Understanding the interconnectedness of structural remodelling processes atrial fibrillation (AF) in patients could identify targets for future therapies. These data include transcriptome sequencing of atrial tissues of patients without AF, with paroxysmal AF and persistent AF (total n=64). In AF patients, epicardial cell gene expression decreased, contrasting with an upregulation of epithelial to mesenchymal transition (EMT) and mesenchymal cell gene expression. Immunohistochemistry demonstrated thickening of the epicardium and an increased proportion of (myo)fibroblast-like cells in the myocardium, supporting enhanced EMT in AF. We furthermore report an upregulation of endothelial cell proliferation, angiogenesis and endothelial signalling. EMT and endothelial cell proliferation concurred with increased interstitial (myo)fibroblast-like cells and extracellular matrix gene expression including enhanced tenascin-c, thrombospondins, biglycan and versican. Morphological analysis discovered increased and redistributed glycosaminoglycans and collagens in the atria of AF patients. Signalling pathways, including cell-matrix interactions, PI3K-AKT and Notch signalling, that could regulate mesenchymal cell activation, were upregulated. These dynamic processes may provide targets to prevent or reverse structural remodelling.

  Study EGAS00001005295

• Assessment of RNA-Seq Sample Preparation Methodology

  The goal of this study was two-fold, to determine if common purification techniques have any effect on downstream differential expression analysis and to evaluate combinations of alignment and differential expression software for reliability. To this end, blood was collected from three individuals and pooled during and after extraction. After pooling and mixing was completed, samples were divided into three aliquots for testing. The first aliquot was a control and not purified or concentrated. It was diluted to produce samples of varying concentrations for testing. The second aliquot was diluted to 20 ng/μL and used to test six different variations on the AMPure XP bead purification procedure. The last aliquot was also diluted to varying concentrations of 60, 30, and 9.6 ng/μL and purified using MinElute columns. Samples were submitted for total RNA-Seq library preparation and sequencing. Library preparation was performed using the TruSeq Stranded Total RNA with Ribo-Zero Globin kit (20020612, Illumina Inc.), and 2x150 bp PE sequencing was done on an Illumina NovaSeq 6000 with the S4 reagent kit. Comparisons were made between methods (AMPure vs. unpurified, AMPure vs. MinElute, MinElute vs. unpurified) to assess the effects of purification method on downstream differential expression. Comparisons were also made within methods using the varying concentrations tested for unpurified samples and for MinElute purified samples to assess the effects of concentration on differential expression. Variations of the AMPure procedure were also compared to assess the effectiveness of the variations tested in comparison to the unmodified procedure. A subset of samples was selected for use with alignment and differential expression package comparison. Unpurified high concentration samples eluted in RNAse-free water were compared to unpurified high concentration samples eluted in BR5, a buffer from the PAXgene blood miRNA kit, with the expectation that there should be no or very few differentially expressed genes. Unpurified high concentration samples eluted in RNAse-free water were also compared to unpurified low concentration samples also eluted in RNAse-free water, with possibly a small number of differentially expressed genes anticipated. A third dataset of simulated RNA-seq data was created with a known rate of differential expression. These files were aligned using Bowtie2, HISAT2, kallisto, RSEM, Rsubread, Salmon, and STAR. Results were then analyzed for differential expression using ALDEx2, baySeq, DEGseq, DESeq2, edgeR, limma, NOISeq, PoissonSeq, and SAMseq. Differential expression results of all three comparisons were evaluated to determine which combinations provided the most reliable results for both real and simulated data.

  Study phs003001

• Single-cell RNA-seq profiling of patient derived organoids

  The aim of this study was to profile developement of resistance to cancer treatments using organoids. Two patient derived colorectal cancer organoids were profiled using single cell RNAseq after treatment with a combination of different anticancer drugs. The organoids were profiled both after 4 weeks of treatment and after a drug holiday of 2 weeks.

  Study EGAS50000001025

• Mutation tracking in single cell RNA-Seq reveals consequences of subclonal evolution in acute myeloid leukemia

  Total bone marrow from AML patients was FACS index sorted with staining for Lineage markers (Lin), CD33, CD34, CD38, CD45RA, CD90, CD99, Tim3 and GPR56, and enrichment for Lin- and Lin-CD34+ cells. Single cells were subjected to a modified smart-seq2 protocol with targeting of mutant sites of interest ("MutaSeq”).

  Study EGAS00001003414

• March 2017 data update (bam/fastq) (containing H3K27ac for CEMT_87, RNA-Seq for CEMT_141, CEMT_142, CEMT_145, CEMT_146, H3K27me3, H3K9me3, H3K4me3, H3K4me1, ChIP-Seq Input for CEMT_88, CEMT_90 and CEMT_91) for reference epigenomes generated at Centre for Epigenome Mapping Technologies, Genome Sciences Center, B.C. Cancer Agency as part of the International Human Epigenome Consortium.

  March 2017 data update (bam/fastq) for reference epigenomes generated at Centre for Epigenome Mapping Technologies (Canadian Epigenetics, Environment and Health Research Consortium), Genome Sciences Center, B.C. Cancer Agency, Vancouver, Canada as part of the International Human Epigenome Consortium.

  Dataset EGAD00001003272

• Gene Expression in an African American Schizophrenia Dataset

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

  Study phs002842

• RNA sequencing data from glioblastoma primary cell lines treated with indisulam data access committee

  Four glioblastoma primary cell lines (BT592/pt1, BT1007/pt2, BT1009/pt3; BT1012/pt4) were derived from glioblastoma Cavitron Ultrasonic Aspirator (CUSA) material and cultured as neurospheres (GB-NS) in DMEM/F12+Glutamax medium containing B27 supplement (Thermofisher, Waltham, MA, USA) and the mitogenic factors epidermal growth factor (EGF) and fibroblast growth factor b (bFGF) (Peprotech, Rocky Hill, NJ, USA). GB-NS were confirmed mycoplasma-free by PCR Test. Cells were seeded in triplicate in 6-well plates, treated with 5 μM indisulam (Sigma Aldrich, Merk), or vehicle control, and collected after 6 and 16 hours for RNA extraction. Total RNA was extracted using TRIzol Reagent (Life Technologies, Thermo Fisher, cat n. 15596026), RNeasy Mini Kit, and RNase-Free DNase Set (Qiagen, cat n. 74104 and 79254). Sample QC, RNA library preparations, and sequencing reactions were conducted at GENEWIZ (Germany GmbH). The concentration of RNA was quantified using a Qubit Fluorometer (Life Technologies, Carlsbad, CA) and RNA integrity was assayed using a TapeStation. Sequencing libraries were prepared with RNA with PolyA selection and ERCC spike-in with standard ERCC kit (ERCC RNA Spike-In Mix,Thermo Fisher 4456740). Final libraries (Illumina, RNA with PolyA selection) were sequenced on a NovaSeq platform, producing 2×150 base paired-end reads.

  Dac EGAC50000000407

• RNA sequencing of a glioblastoma PDX cohort

  Patient-derived xenograft (PDX) models represent an essential tool for the preclinical evaluation of novel targeted therapies and their combinations. In this project, we molecularly characterized a comprehensive panel of glioma PDX models with well conserved molecular and pathological features over multiple passages. PDX samples were subjected to next generation sequencing (NGS). Characterization of cancer-relevant features including driver mutations or cellular processes was performed using mutational and gene expression data in order to identify potential biomarker or treatment targets in glioblastoma. In summary, we reported a newly established and molecularly characterized panel of glioma PDX models with high relevance for translational preclinical research.

  Study EGAS00001007119

• Understanding Self- Organising Capacity of Stem Cells during Implantation and Early Post-implantation Development in vitro and in vivo: Implications for Human Development (2020-04-20)

  The aim of this project is to differentiate human embryonic stem cells to an extra-embryonic fate, specifically the hypoblast. This is of uttermost importance given the current lack of human hypoblast stem cells. We hypothesized that the pluripotent characteristics of the starting human embryonic stem cell population may dictate the competency for extra-embryonic cell fate specification. Based on this hypothesis and using human embryonic stem cells maintained in different naïve-like culture regimes, we have now developed conditions that allow the differentiation of human embryonic stem cells to a stable GATA6+ SOX2- population. This suggests that these cells may be putative human hypoblast stem cells. To validate this finding here we propose to perform RNA sequencing experiments of the differentiated human embryonic stem cells. By comparing their RNA expression profile to the single cell sequencing data of the human embryo that we are currently generating, we will be able to determine the identity of our GATA6+ SOX2- cells, and establish whether they represent the in vivo human hypoblast. This dataset contains all the data available for this study on 2020-04-20.

  Dataset EGAD00001006056

• 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

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

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

  Dataset EGAD00001015660

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

  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 EGAD00001014787

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

  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 EGAD00001014789

• Single Cell RNA-Seq Reveals Malignant and Stromal Programs Associated with Metastasis in Head and Neck Cancer

  This study of head and neck cancer aimed to profile individual malignant, stromal, and immune cells in treatment-naive oral cavity squamous cell carcinoma, with or without matched metastatic cervical lymph nodes. Clinical and pathologic information was collected on all patients, in addition to a fresh biopsy sample from the primary tumor and/or metastatic lymph node at the time of initial surgery. Biopsy specimens were processed for single cell RNA-seq (scRNA-seq) and/or whole exome sequencing (WES). A total of ~6,000 single cells from 18 participants were included in the final dataset. Stromal and immune cells had consistent expression programs across patients. Conversely, malignant cells varied within and between tumors in their expression of signatures related to cell cycle, stress, hypoxia, epithelial differentiation, and partial epithelial-to-mesenchymal transition (p-EMT). One of the top scoring genes in this program was TGFβ-induced (TGFBI), implicating the classic EMT regulator TGFβ. Cells expressing the p-EMT program spatially localized to the leading edge of primary tumors. By integrating single-cell transcriptomes with bulk expression profiles for hundreds of tumors, we refined Head and Neck Squamous Cell Carcinoma (HNSCC) subtypes by their malignant and stromal composition, and established p-EMT as an independent predictor of nodal metastasis, grade, and adverse pathologic features. Our results provide insight into the HNSCC ecosystem and define stromal interactions and a p-EMT program associated with metastasis.

  Study phs001474

• Lineage Plasticity and Immune Cell Heterogeneity Are Coordinately Dysregulated Through Changes in FOXA1 Expression in Bladder Cancers with Squamous Differentiation

  We performed whole exome sequencing (WES) of paired, macrodissected regions of urothelial carcinoma (UC), not otherwise specified (NOS-UC) and squamous differentiation (SqD) from 21 bladder tumors in which distinct and separable NOS-UC and SqD regions were present (total 42 exomes, in addition to matched normal exomes). Mean tumor mutational burden (TMB) was higher in the SqD versus the NOS-UC regions. Consistent with the higher TMB, the SqD regions also had a higher neoantigen burden than the patient-matched NOS-UC regions. The SqD regions also exhibited higher karyotypic complexity with a higher median ploidy relative to the paired NOS-UC regions. Phylogenetic analysis of the WES data confirmed that the morphologically distinct regions of all 21 tumors arose from a shared precursor with a median of 68 (range 16 – 852) non-synonymous mutations shared between the NOS-UC and SqD regions, and the overall data were indicative of early branched evolution. As whole exome sequencing did not identify a shared mutational alteration or genomic signature unique to the SqD or NOS-UC components, we therefore performed expression profiling of separable SqD and NOS-UC regions of 12 bladder tumors from the WES cohort with RNA quality sufficient for whole transcriptome RNA-sequencing (RNAseq, total of 24 regions). Based on centroid clustering analysis of expression data, all 12 regions of SqD were basal-squamous subtype based on the TCGA molecular classification schema. Notably, 8 of the NOS-UC regions were also basal-squamous subtype, with only four pairs exhibiting discordant transcriptional subtypes (three NOS-UC regions were luminal-infiltrated, one luminal-papillary). To better quantify differences in the transcriptional profiles of the paired NOS-UC and SqD regions, we used single sample gene set enrichment analysis (ssGSEA) to assess for enrichment of luminal or basal-squamous gene expression. ssGSEA analysis revealed that the SqD regions expressed higher levels of basal-squamous genes based on the BASE47 gene set and UPK/KRT gene sets; whereas the matched NOS-UC regions expressed higher levels of luminal genes. These differences in the levels of expression of basal and luminal genes were observed even in cases in which both the NOS-UC and SqD regions were basal-squamous subtype based the TCGA molecular classification schema. Exomes from eight patients in this study have previously been deposited in dbGaP: accession number phs001783 as part of the study "Exome Recapture and Sequencing of Prospectively Characterized Clinical Specimens From Cancer Patients".

  Study phs002357

• T-ALL RNA-Seq raw data files

  Bone marrow or peripheral blood samples were collected of adult patients at first diagnosis of T-precursor acute lymphoblastic leukemia. RNA was isolated from mononuclear cells and subjected to mRNA library prep using Poly-A selection and sequencing on a NovaSeq 6000 system. Obtained gene expression profiles and gene fusion calls were used to allocate samples to molecular disease subtypes.

  Study EGAS50000000213

• RNA seq of intestinal biopsies and blood cells from patients with celiac disease and controls

  This study presents RNA sequenced from several intestinal phenotypes such as from patients with damaged intestinal mucosa as well as from those without a damaged intestine but still an autoantibody response to tissue transglutaminase (“potential celiac disease”) or those with treated celiac disease. By combining phenotypes in a much larger number of tissue samples from both peripheral blood and mucosal biopsies, we found highly differentially expressed genes, but also identify genes involved in controlling and triggering celiac disease associated changes. This study identifies molecular mechanisms involved in celiac disease and new possible targets for treatment and for identification of individuals at risk.

  Dataset EGAD50000001597

• Genomic and epigenomic study of Japanese renal cell carcinoma including WGS, RNA-seq, ATAC-seq, and methyl-seq

  We performed whole-genomic and transcriptomic sequencing of more than 100 Japanese renal cell carcinoma (RCC) cases, including clear cell RCC, papillary RCC, chromophobe RCC, and TFE3-translocated RCC. In addition, we conducted epigenome sequencing analyses: the assay for transposase-accessible chromatin sequencing (ATAC-seq) and enzymatic methyl sequencing.

  Study EGAS00001006919

• Tracking the evolution of esophageal squamous cell carcinoma under dynamic immune selection by multi-omics sequencing

  Intratumoral heterogeneity (ITH) has been linked to decreased efficacy of clinical treatments. However, although genomic ITH has been characterized in genetic, transcriptomic and epigenetic alterations are hallmarks of esophageal squamous cell carcinoma (ESCC), the extent to which these are heterogeneous in ESCC has not been explored in a unified framework. Further, the extent to which tumor-infiltrated T lymphocytes (TILs) are directed against cancer cells, but how the immune infiltration acts as a selective force to shape the clonal evolution of ESCC is unclear. In this study, we perform multi-omic sequencing on 186 samples from 36 primary ESCC patients. Through multi-omics analyses, it is discovered that genomic, epigenomic, and transcriptomic ITH are underpinned by ongoing chromosomal instability. Based on the RNA-seq data, we observe diverse levels of immune infiltrate across different tumor sites from the same tumor. We reveal genetic mechanisms of neoantigen evasion under distinct selection pressure from the diverse immune microenvironment. Overall, our work offers an avenue of dissecting the complex contribution of the multi-omics level to the ITH in ESCC and thereby enhances the development of clinical therapy.

  Dataset EGAD00001009482

• RNA, ATAC, ChIP datasets from Ankylosing Spondylitis patients and Healthy Controls

  This dataset consists of functional genomic data from 20 Ankylosing Spondylitis patients and 35 Healthy Controls taken from CD4+ T cells, CD8+ T cells and CD14 Monocytes. It contains 364 paired end fastq files consisting of 104 total RNA-seq samples and 116 ATAC samples, for ChIP there are 46 H3K4me3 samples, 46 H3K27ac samples and 3 H3K4me1 samples, along with 49 paired input samples. The samples were sequenced on Illumina HiSeq4000, Illumina NextSeq500 and Illumina NovaSeq 6000 platforms.

  Dataset EGAD00001008758

• Whole genome landscape of adult T-cell leukemia/lymphoma

  Whole-genome sequence (WGS) data of tumor-normal pairs from 139 ATL patients and RNA sequence (RNA-seq) data of tumors from 28 ATL patients.

  Dataset EGAD00001008289

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

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

  Study phs002031

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

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

  Study phs003448

• 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

• 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

• 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

• Colorectal cancer organoids expressing BRAF (fusion) genes

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

  Study EGAS00001003558

• Healthy_ageing_thymus

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

  Study EGAS00001004311