Search Results - EGA European Genome-Phenome Archive

MultiOMICS study of a pair of infant monozygotic twins with concordant B-cell ALL (WGBS)

B-cell acute lymphoblastic leukemia (B-cell ALL) is the most common cancer in childhood. Studying identical twins with B-cell ALL provides a unique and tractable model for deciphering the developmental timing of pre- and post-natal mutations contributing to clonal evolution. To date, this has mainly focused on major cytogenetic subgroups of childhood B-cell ALL, including MLL fusions, ETV6-RUNX1, hyperdiploidy, and BCR-ABL1. However, formal demonstration of the prenatal origin and “backtracking” the natural history of the leukemia remains understudied in “B-other”/Normal Karyotype (NK) B-cell ALL. To characterize the epigenetic landscape of this particular leukemia subtype, we performed DNA bisulfite-sequencing on a pair of 8-month-old monozygotic twins diagnosed with concordant “B-other”/NK B-cell ALL.

Dataset EGAD00001005018

Pediatric Papillary Thyroid Carcinoma RNA-Seq

Papillary thyroid carcinoma (PTC) is the most common malignancy amongst adolescent and young adult women. Clinical and genomic features of PTC in children differ from those in adults. Historically, a substantial proportion (~50%) of pediatric PTC (PPTC) lack any of the common driver mutations typical of adult PTC. Thus, molecular diagnostics developed for adults may not be valid in children. We applied novel bioinformatic pipelines to RNASeq analysis of pediatric thyroid tumors to identify known and novel fusion oncogenes as drivers of tumorigenesis and to define gene expression patterns among tumors. Analysis of 38 samples resulted in identification of 16 previously reported oncogenic fusions and three novel candidate fusions. In conducting these analyses, we reduced the proportion of cancers with unknown drivers from approximately 50% to 6%. Gene expression analysis identified three distinct categories of tumors. These studies represent the first unbiased genomic approach to childhood PTC. Novel bioinformatic pipeline. GEC Identified 3 distinct classes. DICER1, PTEN segregated with benign tumors. >90% of PPTC carry identifiable ongognic driver variants, however, the landscape differs from that in adults and existing tools for tumor genotyping in adults may not be applicable to pediatric tumors. BACKGROUND: Papillary thyroid carcinoma demonstrates clinically distinct behaviour in pediatric patients when compared to adults, with increased regional and distal metastases, more frequent recurrence, yet overall highly favourable prognosis. Recent tumor genotyping has demonstrated that in pediatric tumors, ***-*** do not carry the most common genomic variants found in adult PTC, leaving a substantial proportion of driver-unknown or “dark matter” tumors. METHODS: We established a cohort of 52 pediatric PTC and performed quantitative-PCR genotyping followed in a subset by whole exome and whole transcriptome (RNASeq) sequencing using a novel, robust bioinformatic pipeline to identify oncogenic fusion transcripts. Gene expression analysis was used to derive distinct clusters. RESULTS: In 52 pediatric PTC, genotyping identified BRAFV600E(21%), NRASQ61R(2%), RET/PTC1 (21%), RET/PTC3 (17%) or PAX8-PPAR in 33 (63%) of tumors. Among the remaining 19 tumors, oncogenic variants were identified in /19 (%). The variants included SNVs in DICER1, ETV6-NTRK3, PTEN and 9 additional oncogenic fusions including two previously-undescribed fusion transcripts, *** and ***. CONCLUSIONS: The genomic landscape of PTC differs between children and adults. In children, there is a substantial fraction of tumors that are unexplained by variants typical of adult PTC. This discrepancy is explained largely by uncommon oncogenic fusions which can be identified by whole-transcriptome analysis.

Study EGAS00001005182

Genomic Basis of Phenotypic Variability of Complex Disorders

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

Study phs002450

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

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

Study EGAS00001001418

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

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

Study EGAS00001008073

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

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

Study EGAS00001002328

Liver biopsy derived induced pluripotent stem cells provide an unlimited supply for the generation of patient-specific hepatocyte-like cells

Background: Hepatocyte-like cells (HLCs) differentiated from induced pluripotent stem cells (iPSCs) have emerged as a promising cell culture model to study metabolism, biotransformation, viral infections and genetic liver diseases. However the derivation of iPSCs is restricted by ethical and procedural constraints. Furthermore, incomplete HLC differentiation remains a major challenge. iPSC may carry-on a tissue of origin dependent expression memory influencing iPSC differentiation into different cell types. If liver derived iPSCs (Li-iPSCs) maintain a liver specific gene expression profile is not known. It is also not known, if Li-iPSCs would allow the generation of more fully differentiated HLCs.Methods: Primary liver cells (PLCs) were expanded from liver needle biopsies and reprogrammed into Li-iPSCs using a non-integrative Sendai virus-based system. Li-iPSCs were differentiated into HLCs using established differentiation protocols. The HLC phenotype was characterized at the protein, functional and transcriptional level. RNA sequencing data were generated from the originating liver biopsies, the Li-iPSCs, fibroblast derived iPSCs, and differentiated HLCs, and used to characterize and compare their transcriptome profiles.Results: PLCs were obtained from small pieces of liver biopsies and could be reprogrammed into Li-iPSCs. Li-iPSCs indeed retain a liver specific transcriptional footprint. Li-iPSCs can be propagated to provide an unlimited supply of cells for differentiation into Li-HLCs. Similar to HLCs derived from fibroblasts, Li-HLCs could not be fully differentiated into hepatocytes. Relative to the originating liver, Li-HLCs showed lower expression of liver specific transcription factors and increased expression of genes involved in the differentiation of other tissues.Conclusions: PLCs and Li-iPSCs obtained from small pieces of human needle liver biopsies constitute a novel unlimited source for the production of HLCs. Despite the preservation of a liver specific gene expression footprint in Li-iPSCs, the generation of fully differentiated hepatocytes cannot be achieved with the current differentiation protocols.

Study EGAS00001002676

A circulating biomarker for severity of facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable skeletal myopathy. Clinical trials for FSHD are hindered by heterogeneous biomarkers poorly associated with clinical severity, requiring invasive muscle biopsy. Macroscopically FSHD presents with slow fatty replacement of muscle, rapidly accelerated by inflammation. Mis-expression of the transcription factor DUX4 is currently accepted to underlie FSHD pathogenesis and mechanisms including PAX7 target gene repression have been proposed. Here we perform RNA-sequencing on MRI guided inflamed (n=24) and isogenic non-inflamed (n=24) muscle biopsies from the same clinically characterised FSHD patients, alongside isogenic peripheral blood mononucleated cells from a subset of patients (n=13), and unaffected controls. We employ multivariate models to evaluate the clinical associations of 5 published FSHD transcriptomic biomarkers. We demonstrate that PAX7 target gene repression can discriminate control, inflamed and non-inflamed FSHD muscle independently of age and sex (p<0.013), while the discriminatory power of DUX4 target genes is limited to distinguishing FSHD muscle from control. Importantly, the level of PAX7 target gene repression in non-inflamed muscle associates with clinical assessments of FSHD severity (p=0.04). DUX4 target gene biomarkers in FSHD muscle show associations with lower limb fat fraction and D4Z4 array length, but not clinical assessment. Lastly, PAX7 target gene repression in FSHD muscle correlates with the level in isogenic peripheral blood mononucleated cells (p=0.002). A refined PAX7 target gene biomarker comprising 143/601 PAX7 target genes computed in peripheral blood (the FSHD muscle-blood biomarker) associates with clinical severity in FSHD patients (p<0.036). Our new circulating biomarker validates as a classifier of clinical severity in an independent data set of 54 FSHD patient blood samples, with improved power in older patients (p=0.03). In summary, we present the minimally invasive FSHD muscle-blood biomarker of FSHD clinical severity valid in patient muscle and blood, with potential utility in routine disease monitoring and clinical trials.

Study EGAS00001007350

A genotype-phenotype study of tumors from patients with inherited mutations in DNA repair genes

This study investigated the genomic landscape of tumors from patients with inherited mutations in DNA repair genes. The aim of this study was to determine the relationship between locus-specific loss of heterozygosity (LOH) and genomic signatures of homologous recombination-based DNA repair (HR) deficiency in tumors associated with pathogenic germline variants in DNA repair genes. Cohorts of pathogenic germline variant (PGV) carriers with different cancer types were identified via a cancer genetics registry and clinical practice at Penn Medicine and through our collaborator network. Tumors were analyzed with next generation sequencing (whole exome and/or targeted panel sequencing) with or without transcriptomic analyses. Sub-study #1 "Non-breast/ovarian tumors in germline BRCA1/2 carriers". This study investigated primary non-breast/ovarian tumors that arise in individuals with BRCA1/BRCA2 pathogenic germline variants (PGVs). Whole exome sequencing (WES) using Agilent All-Exon Kit v7 hybrid capture probes was performed on matched tumor and blood germline normal DNA on an Illumina HiSeq. Individual WES data from DNA samples from tumor and normal are available. Sub-study #2: "Breast tumors in germline CHEK2 carriers". This study investigated primary breast tumors that arise in individuals with CHEK2 pathogenic germline variants (PGVs). Whole exome sequencing (WES) using Agilent All-Exon Kit v7 hybrid capture probes was performed on matched tumor and blood germline normal DNA on an Illumina HiSeq. Whole transcriptome RNAseq was performed on tumors. Individual WES data from DNA samples from tumor and normal and RNAseq data from tumor are available. Sub-study #3: "Breast tumors in germline TP53 carriers". This study investigated primary breast tumors that arise in individuals with TP53 pathogenic germline variants (PGVs). Whole genome sequencing (WGS) and/or custom panel targeted sequencing was performed on matched tumor and blood germline normal DNA on an Illumina HiSeq. Whole transcriptome RNAseq was performed on tumors and normal breast tissue. Individual WGS or targeted data from DNA samples from tumor and normal and RNASeq data from tumor and normal are available. Sub-study #4: "Prostate tumors in germline DNA repair carriers". This study investigated prostate tumors that arise in individuals with a variety of pathogenic germline variants (PGVs). Whole genome sequencing (WGS) and/or custom panel targeted sequencing was performed on matched tumor and blood germline normal DNA on an Illumina HiSeq. Whole transcriptome RNAseq was performed on tumors and normal prostate tissue. Individual WGS or targeted data from DNA samples from tumor and normal and RNASeq data from tumor and normal are available.

Study phs003348

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

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

Study phs002780

A Multimodal Atlas of Human Brain Cell Types

"A multimodal atlas of human brain cell types" includes sample data targeting two cell types that show species differences between mouse and human. First, it includes a detailed transcriptomic, morphological, and electrophysiological characterization of cell types in layer 1 of human middle temporal gyrus, focusing primarily on different inhibitory cell types. This project also includes gene expression data collected from nuclei in layer 5 of human fronto-insula, with a goal of identifying transcriptomic signatures of Von Economo neurons. Control samples collected as part of the same experiment are also included in the data set. This study was conducted as part of a collaboration between the Allen Institute for Brain Science, the University of Szeged, and the J. Craig Venter Institute. Collaborators request that publications resulting from these data cite their original publication: Transcriptomic and morphophysiological evidence for a specialized human cortical GABAergic cell type (PMID: 30150662).

Study phs001791

DNA and RNA sequencing of single human haploid germ cells

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

Study phs002279

Whole genome sequencing data from tumor and normal samples

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

Dataset EGAD50000001909

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

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

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

Dataset EGAD00001004140

11926 results for "cancer rna-seq"

in 18.34 milliseconds.