Search Results - EGA European Genome-Phenome Archive

Whole genome sequencing, SNP array and RNA-seq of uveal melanomas

We conducted whole genome sequencing and DNA SNP array of 12 uveal melanoma genomes and their matched DNA from blood. We also conducted RNA-seq of the 12 tumour samples.

Study EGAS00001000472

Osteosarcoma_whole_genome_rearrangement_screen

Genomic libraries (500 bps) will be generated from total genomic DNA derived from Osteosarcoma cancer patients and subjected to short paired end sequencing on the llumina platform. Paired reads will be mapped to build 37 of the human reference genome to facilitate the generation of genome wide copy number information, and the identification of novel rearranged cancer genes and gene fusions.

Study EGAS00001000330

Personalized RNA mutanome vaccines mobilize poly-specific therapeutic immunity against cancer RNA-Seq

The goal of the sequencing study was to determine the neoantigen profile of melanoma patients. We performed whole exome sequening of the tumor and respective patien matched PBMC samples (2x50nt) and RNA-sqeuencing of the tumor samples. The was analyzed for somatic point mutaions, using the RNA-Seq data to determine mutation expression. This study contains the patient matched RNA-Sequencing data.

Study EGAS00001003307

Whole Exome Sequencing of cohorts of Mutant Braf mouse model melanoma DNA and germline DNA.

Whole Exome Sequencing of cohorts of Mutant Braf mouse model melanoma DNA and germline DNA. The cohorts are (1) Mutant Braf mouse model melanomas, (2) Mutant Braf mouse model melanomas from UVR exposed mice and (3) Mutant Braf mouse model melanomas from UVR exposed, sunscreen protected mice.

Study EGAS00001000729

MYCN Amplification and ATRX Mutations are Incompatible in Neuroblastoma

Age and stage of disease at diagnosis is associated with outcome in neuroblastoma. A previous study identified ATRX point mutations and in-frame deletions in older patients with stage 4 disease. To validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables, we collaborated with the Children’s Oncology Group and sequenced the entire ATRX genomic locus in 475 neuroblastoma tumor samples to validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables. Older age at diagnosis, stage 4 disease, and unfavorable histology were significantly associated with mutations in the ATRX gene; MYCN amplification was mutually exclusive from ATRX mutations. To directly test if ATRX mutations and MYCN amplification are incompatible in neuroblastoma, we mutated the ATRX gene in MYCN amplified neuroblastoma cell lines and ectopically expressed MYCN in ATRX mutant neuroblastoma cell lines. Both approaches showed that ATRX and MYCN amplification are incompatible in neuroblastoma in culture and in vivo. To gain a better understanding of the underlying molecular and cellular mechanisms of this incompatibility, we performed a series of experiments including telomere analysis, epigenetic profiling, electron microscopic analysis and metabolic profiling. We discovered that induction of MYCN expression in ATRX mutant neuroblastomas cells leads to metabolic reprogramming, mitochondrial dysfunction, increased reactive oxygen species and replicative stress. We propose that this leads to synthetic lethality in ATRX mutant neuroblastomas because of the underlying replicative stress in those cells as a result of dysfunction of this essential histone chaperone.

Study EGAS00001003257

Whole genome sequencing and whole exome sequencing of mucosal melanoma

We conducted whole genome sequencing of 5 mucosal melanomas and their matched DNA from blood. We conducted whole exome sequencing of a further 5 mucosal melanomas and their matched DNA from blood.

Study EGAS00001000474

Genetic and epigenetic characterization of adenoid cystic carcinoma

Translocation events are frequent in cancer and may create chimeric fusions or regulatory rearrangements that drive oncogene overexpression. Although regulatory rearrangements are increasingly recognized in hematopoietic and even solid tumors, the underlying mechanisms remain obscure. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a unifying theme in adenoid cystic carcinoma (ACC). Whole genome sequencing data and chromatin state maps for 13 primary tumors and xenografts reveal distinct chromosomal rearrangements that juxtapose super-enhancer clusters to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers physically interact with the MYB promoter. Remarkably, the MYB product itself binds to the translocated enhancers, thus creating a positive feedback loop that sustains its own expression. MYB also binds a large number of active enhancers that drive opposing regulatory programs in alternate cell lineages in ACC. MYB cooperates with the transcription factor TP63 in the myoepithelial component of the tumors, but promotes a Notch program in the luminal epithelial component. Bromodomain inhibitors, which disrupt enhancer function, slow tumor growth in ACC primagraft models in vivo, but are ineffective against high grade tumors that harbor activating mutations in the Notch pathway. Thus, our study identifies super-enhancer translocations as a unifying feature of ACC, and provides insight into the mechanism by which sustained MYB overexpression drives alternate cell fates in this disease.Processed ChIP-seq data is available on GEO under accession number GSE76465.

Study EGAS00001001457

WGS_of_healhy_mesothelial_cells_and_primary_mesothelima_cell_lines

Mesothelioma is a rare and aggressive cancer associated with previous exposure to asbestos. Currently there are no effective treatments for mesothelioma and majority of patients will die within a year after diagnosis. Although a large number of preclinical and clinical trials assessed the efficacy of various therapuetic modes, none of the tested molecules entered the clinic. Thus, to better understand what drives mesothelioma carcinogenesis and to identify novel targets for therapy, here we aim at performing WGS of heathy mesothelial cells and a panel of mesothelioma cells lines. Primary mesothelial cell lines were established in prof. Marciniak lab from human pleural specimens, while primary mesothelioma cell lines were obtained from MesobanK. Both models constitute low passage cell lines and therefore should closely represent genetics of the original tissue.

Study EGAS00001005559

WGS of 32 paired SRCC samples

Signet-ring cell carcinoma (SRCC) has specific epidemiology and oncogenesis in gastric cancer, however with no systematical investigation for prognostic genomic features. We systematically investigated the clinical characteristics of 1,868 Chinese gastric cancer patients, and found that signet-ring cell content was significantly related to multiple clinical characteristics and treatment outcomes. Next, we performed whole genome sequencing for 32 SRCC patients with > 80% signet-ring cells in tumor, and identified frequent5 CLDN18-ARHGAP26/6 fusion (25%). With additional 797 patients for validation, we notice that prevalence of CLDN18-ARHGAP26/6 fusion is positively associated with signet-ring cell content (P = 4.1 x 10 -9 ), younger age at diagnosis (P = 4.2 x 10 -10 ), female/male ratio (P = 1.7 x 10 -9 ), and advanced TNM stage (P = 1.7 x 10 -5 ). Importantly, patients with CLDN18-ARHGAP26/6 fusion had a worse survival outcomes (P = 0.03), and got no benefit0 from platinum/fluoropyrimidines based chemotherapy (P = 0.92) compared to fusion-free patients (P = 0.001), which is consistent with the fact of oxaliplatin/fluoropyrimidine resistance acquired in CLDN18-ARHGAP26 introduced gastric cancer cell lines. Overall, this study provides new insights into the clinical and genomic features of SRCC, and highlights the importance of frequent CLDN18-ARHGAP26/6 fusions in chemotherapy response for SRCC.

Study EGAS00001002668

Ewings_Sarcoma_RNA_seq_drug_sensitivity

Illumina RNA-Seq will be performed on four Ewing’s sarcoma cell lines and two control cell lines. RNA was extracted from all the lines using a basic Trizol extraction protocol.

Study EGAS00001000419

Foundation Medicine Genomic Data Used to Identify Prognostic Markers and Fusion Genes in Multiple Myeloma

The purpose of this study is to identify prognostic markers and treatment targets using a clinically certified sequencing panel in multiple myeloma. Mutational burden was associated with maf and proliferation gene expression groups, and a high-mutational burden was associated with a poor prognosis. We identified homozygous deletions that were present in multiple myeloma within key genes, including CDKN2C, RB1, TRAF3, BIRC3 and TP53, and that bi-allelic inactivation was significantly enriched at relapse. Alterations in CDKN2C, TP53, RB1 and the t(4;14) were associated with poor prognosis. Alterations in RB1 were predominantly homozygous deletions and were associated with relapse and a poor prognosis which was independent of other genetic markers, including t(4;14), after multivariate analysis. Bi-allelic inactivation of key tumor suppressor genes in myeloma was enriched at relapse, especially in RB1, CDKN2C and TP53 where they have prognostic significance. In addition, chromosomal rearrangements that result in oncogenic kinase activation are present in many solid and hematological malignancies, but none have been reported in multiple myeloma (MM). Here we detected fusion genes in 1.5% of patients. These fusion genes were in-frame and the majority of them contained kinase domains from either receptor tyrosine kinases (ALK, ROS1, NTRK3, and FGFR1) or cytoplasmic kinases (BRAF, MAP3K14, and MAPK14) which would result in the activation of MEK/ERK, NF-κB or inflammatory signaling pathways. Fusion genes were present in smoldering MM, newly diagnosed MM and relapse patient samples indicating they are not solely late events. Most fusion genes were subclonal in nature, but one EML4-ALK fusion was clonal indicating it is a driver of disease pathogenesis. Samples with fusions of receptor tyrosine kinases were not found in conjunction with clonal Ras/Raf mutations indicating a parallel mechanism of MEK/ERK pathway activation. Fusion genes involving MAP3K14 (NIK), which regulates the NF-κB pathway, were detected as were t(14;17) rearrangements involving NIK in 2% of MM samples. Activation of kinases in myeloma through rearrangements presents an opportunity to use treatments existing in other cancers.

Study EGAS00001002874

An Unusual Genomic Variant of Pancreatic Ductal Adenocarcinoma with an Indolent Clinical Course

We describe a patient with biopsy-proven locally advanced pancreatic ductal adenocarcinoma (PDA) who underwent a modified course of gemcitabine and stereotactic body radiation therapy. Whole exome sequencing of this patient's tumor revealed a simple genome landscape with no evidence of mutations, copy number changes, or structural alterations in genes most commonly associated with PDA, i.e. KRAS, CDKN2A, TP53 or SMAD4. An analysis of germline DNA revealed no pathogenic variants of significance. Whole exome and whole genome sequencing identified a somatic mutation of RNF213 and an inversion/deletion of CTNNA2 as the genetic basis of this PDA. Although PDA is classically characterized by a predictable set of mutations, these data suggest that alternate genetic paths to PDA may exist, which can be associated with a more indolent clinical course.

Study EGAS00001002192

Genetic landscape of relapsed DLBCL

Diffuse large B-cell lymphoma is a heterogenous malignancy that can arise de novo or can result from histologic transformation from indolent lymphomas. Considerable sequencing effort has been placed on determining the mutations that are common at the time of diagnosis in DLBCL but little is known about the nature of relapsed disease. We have sequenced tumour and matched normal exomes of 38 rrDLBCL patients including 25 de novo DLBCLs and 13 transformed lymphomas to identify somatic mutations and copy number alterations associated with relapse.

Study EGAS00001001553

Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study

This data comes from a genome-wide association study in 1090 cases of scarring trachoma and 1531 controls from a non-selected multi-ethnic sampling of The Gambia. The data constitutes hard base calls at 1457295 features from the HumanOmni2.5-8 (Omni2.5) BeadChip. Three genotype-calling algorithms were used : Illuminus (Bioinformatics 2007, 23:2741–6), GenCall (Patent US7035740, 2006) and GenoSNP (Bioinformatics 2008, 24:2209–14). Data from each calling algorithm was filtered to retain only SNPs with a call rate ≥ 0.98. To obtain a merged set of SNPs, all genotypes that matched across call-sets were retained, whilst those that mismatched between call-sets were set to missing. Genotypes that were present in one call set and missing in others were also retained. The merged SNP set contained 1467876 SNPs, which were filtered for call rate ≥ 0.99 and Hardy-Weinberg equilibrium P value < 5 x 10-8. The case and control groups were approximately equivalent with respect to gender and ethnicity. 70% of cases and 63% of controls were female. Self-described ethnic composition of the cases was 29% Jola, 27% Mandinka, 21% Wolof, 10% Fula and 13% other/no data. In controls the composition was 25% Jola, 24% Mandinka, 21% Wolof, 8% Fula and 22% other/no data. Median age was 49 (range 32-60) in cases and 37 (range 12-52) in controls (t = -3.6019, P = 0.0003). This data set has substantial levels of population structure, which reflects extensive networks of pairwise kinship between participants. First pass tests of association in EMMAX (Nat Genet 2010, 42:348–54) led to genome-wide deflation of the test statistics (λ = 0.982 SE = 1.74e-05) that could be directly attributed to the differing age and gender distributions between cases and controls. Modelling the phenotype to adjust for age and gender successfully controlled for this and no genome-wide deviation from the null was subsequently observed (λ = 1.001, SE = 8.7 x 10-7). Neither principal components analysis (PCA) nor tests of proportional identity by state (IBS) variance in PLINK identified significant levels of within or between group genetic variance. The phenotypes submitted in the phenotypes file include both uncorrected (case/control status) and corrected (age and gender corrected case/control status) phenotypes, as well as age and self-described ethnicity.

Study EGAS00001001516

Whole genome, whole exome, and targeted sequencing of high-grade meningioma tumor samples.

High-grade meningiomas frequently recur and are associated with high rates of morbidity and mortality. To determine the factors that promote the development and evolution of these tumors, we analyzed the genomes of 274 high-grade meningiomas and compared this information with data from 456 low-grade meningiomas. High-grade meningiomas had a higher mutation burden than low-grade meningiomas but did not harbor any statistically significant mutated genes aside from NF2. High-grade meningiomas also possessed significantly elevated rates of chromosomal gains and losses, especially among tumors with monosomy 22. Meningiomas previously treated with adjuvant radiation had significantly more copy number alterations than radiation-induced or radiation-naïve meningiomas. Across serial recurrences, genomic disruption preceded the emergence of nearly all mutations, remained largely uniform across time, and when present in low-grade meningiomas, correlated with subsequent progression to a higher grade. In contrast to the largely stable copy number alterations, mutations were strikingly heterogeneous across tumor recurrences, likely due to extensive geographic heterogeneity in the primary tumor. While high-grade meningiomas harbored significantly fewer overtly targetable alterations than low-grade meningiomas, they contained numerous mutations that are predicted to be neoantigens, suggesting that immunologic targeting may be of therapeutic value.

Study EGAS00001002294

MED12L Gene Alterations Define Aggressive BRCA2-Mutant Prostate Cancers

Germline mutation of BRCA2 increases the lifetime risk of developing prostate cancer (PCa) by over 700%. BRCA2-mutant PCa have poorer prognosis than sporadic PCa, with rapid development of metastatic, castrate-resistant prostate cancer (mCRPC) and 5-year cancer-specific survival rates of ~50-60%1-4. Despite this, unique genomic driver events that explain the aggressiveness of localized BRCA2-mutant PCa are lacking. We used whole-genome sequencing to fully characterize 14 BRCA2-mutant PCa and demonstrate that BRCA2-mutant PCa is associated with increased genetic instability and a mutually exclusive mutational burden when compared to sporadic PCa. Importantly, BRCA2-mutated cancers are defined by unique copy number gains and hypomethylation events, including alterations in the MED12L/MED12 axis, which are found solely in mCRPC and are enriched in PCa tumours harbouring aggressive intraductal carcinoma (IDC-P) pathologic sub-types. Our findings begin to explain the clinical entity of BRCA2-mutated PCa as these tumours have a unique and aggressive genotype de novo, associated with IDC-P and mCRPC, even in the hormone-naive setting.

Study EGAS00001001615

These are from Korean HCC samples with exome sequencing

Hepatic resection is the most curative treatment option for early-stage hepatocellular carcinoma, but is associated with a high recurrence rate, which exceeds 50% at 5 years after surgery. Understanding the genetic basis of hepatocellular carcinoma at surgically curable stages may enable the identification of new molecular biomarkers that accurately identify patients in need of additional early therapeutic interventions. Whole exome sequencing and copy number analysis was performed on 231 hepatocellular carcinomas (72% with hepatitis B viral infection) that were classified as early-stage hepatocellular carcinomas, candidates for surgical resection. Recurrent mutations were validated by Sanger sequencing. Unsupervised genomic analyses identified an association between specific genetic aberrations and postoperative clinical outcomes. Recurrent somatic mutations were identified in 9 genes, including TP53, CTNNB1, AXIN1, RPS6KA3, and RB1. Recurrent homozygous deletions in FAM123A, RB1, and CDKN2A, and high-copy amplifications in MYC, RSPO2, CCND1, and FGF19 were detected. Pathway analyses of these genes revealed aberrations in the p53, Wnt, PIK3/Ras, cell cycle, and chromatin remodelling pathways. RB1 mutations were significantly associated with cancer-specific and recurrence-free survival after resection (p = 0.016 and p = 0.001, respectively). FGF19 amplifications, known to activate Wnt signalling, were mutually exclusive with CTNNB1 and AXIN1 mutations, and significantly associated with cirrhosis (p = 0.017). RB1 mutations can be used as a prognostic molecular biomarker for resectable hepatocellular carcinoma. Further study is required to investigate the potential role of FGF19 amplification in driving hepatocarcinogenesis in patients with liver cirrhosis and to investigate the potential of anti-FGF19 treatment in these patients.

Study EGAS00001000604

Automated system for scoring hematoxylin and eosin-stained ovarian cancer sections by identifying single cells uncovered that stromal cell ratio is a significant predictor for overall survival and progression-free survival.

Concerted efforts in genomic studies have revealed profound insights in prognostic ovarian cancer subtypes. On the other hand, abundant histology slides have been generated to date, yet their uses remain very limited and largely qualitative. Our goal is to develop automated histology analysis as an alternative subtyping technology for ovarian cancer that is cost-efficient and do not rely on DNA quality. We develop an automated system for scoring hematoxylin and eosin-stained (H&E) primary tumour sections of 91 late-stage ovarian cancer to identify single cells including cancer and stromal cells. We demonstrated high accuracy of our system based on expert pathologists’ scores (cancer=97.1%, stromal=89.1%) as well as compared to immunohistochemistry scoring (correlation=0.87). Quantitative stromal cell ratio is significantly associated with poor overall survival after controlling for clinical parameters including debulking status and age (multivariate analysis p=0.0021, HR=2.54, CI=1.40-4.60) and progression-free survival (multivariate analysis p=0.022, HR=1.75, CI=1.09–2.82). We demonstrate how automated image analysis enables objective quantification of microenvironmental composition of ovarian tumours. Our analysis reveals a strong effect of the tumour microenvironment on ovarian cancer progression and highlights the potential of therapeutic interventions that target the stromal compartment or cancer-stroma signalling in the stroma-high, late-stage ovarian cancer subset.

Study EGAS00001001694

Chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in paediatric high grade glioma cell lines KKNS4 and SF188 with and without a G34V mutation in H3F3A

Glioblastomas of children and young adults have a median survival of only 12-15months and are clinically and biologically distinct from histologically similar cancers in older adults1. They are defined by highly specific mutations in the gene encoding the histone H3.3 variant H3F3A2, occurring either at or close to key residues marked by methylation for regulation of transcription – K27 and G34. We performed chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in order to test the hypothesis that, rather than total H3K36me3, the G34V mutation may instead result in differential binding of the trimethyl mark throughout the genome. Compared with H3F3A wild-type SF188 paediatric GBM cells, H3K36me3 was found to be significantly differentially bound in KNS42 cells at 5130 distinct regions of the genome corresponding to 156 genes. Concurrent whole genome DNA sequencing demonstrated that this was not confounded by copy number alterations. The cerebral hemispheric-specific G34 mutation drives a distinct expression signature through differential genomic binding of the K36 trimethylation mark (H3K36me3). The transcriptional program induced recapitulates that of the developing forebrain, and involves numerous markers of stem cell maintenance, cell fate decisions and self-renewal. Critically, H3F3A G34 mutations cause profound upregulation of MYCN, a potent oncogene which is causative of glioblastomas when expressed in the correct developmental context. This driving aberration is selectively targetable in this patient population by inhibiting kinases responsible for stabilisation of the protein.

Study EGAS00001001437

Whole genome shotgun sequencing and somatic mutations data in Hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common malignancy of the liver. Genomic analysis is conducted to identify genetic alterations in driver genes which are all druggable targets for cancer therapy.

Study EGAS00001002578

Whole-Exome Sequencing of Salivary Gland Mucoepidermoid Carcinoma

Mucoepidermoid carcinoma (MEC) is the most common salivary gland malignancy. To explore the genetic origins of MEC, we performed systematic genomic analyses of these tumors.Experimental DesignWhole-exome sequencing and gene copy number analyses were performed for 18 primary cancers with matched normal tissue. Fluorescence in situ hybridization (FISH) was used to determine the presence or absence of the MECT1-MAML2 translocation in 17 tumors.ResultsTP53 was the most commonly mutated gene in MEC (28%), and mutations were found only in intermediate- and high-grade tumors. Tumors with TP53 mutations had more mutations overall than tumors without TP53 mutations (p=0.006). POU6F2 was the second most frequently mutated gene, found in three low-grade MECs with the same in-frame deletion. Somatic alterations in IRAK1, MAP3K9, ITGAL, ERBB4, OTOGL, KMT2C, and OBSCN were identified in at least two of the 18 tumors sequenced. FISH analysis confirmed the presence of the MECT1-MAML2 translocation in 15 of 17 tumors (88%).ConclusionsThrough these integrated genomic analyses, MECT1-MAML2 translocation and somatic TP53 and POU6F2 mutations appear to be the main drivers of mucoepidermoid carcinoma.

Study EGAS00001002811

RNAseq of 76 samples from Uveal Melanoma tumors

Hotspot mutations in the spliceosome gene SF3B1 are reported in 20% of uveal melanomas. SF3B1 is involved in 3'-splice site (3'ss) recognition during RNA splicing; however, the molecular mechanisms of its mutation have remained unclear. Here we show, using RNA-Seq analyses of uveal melanoma, that the SF3B1 R625/K666 mutation results in deregulated splicing at a subset of junctions, mostly by the use of alternative 3'ss

Study EGAS00001002932

HumanMethylation450K data from Purified Plasma Cells of Monoclonal gammopathy of unknown significance and Multiple myeloma patients and Healthy donors

Study EGAS00001000841

Pre-neoplastic somatic mutations including MYD88L265P in lymphoplasmacytic lymphoma

Normal-cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B-lymphocytes prior to lymphoma remains uncertain. In this study, we sequenced multiple stages of the B-lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of high-prevalence of mutated MYD88.

Study EGAS00001005656

Copy number profiling of primary samples and cell lines of retinoblastoma

While RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes. We performed SCNA profiling of 45 primary retinoblastoma samples and 8 retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N=310). Comprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N=34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades. Since focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.

Study EGAS00001001715

5524 results for "*oma"

in 25.79 milliseconds.