Search Results - EGA European Genome-Phenome Archive

Comprehensive Transcriptional Analysis of Early Stage Urothelial Carcinoma using whole transcriptome sequencing

The expression profile and sequence variants of 476 early stage urothelial carcinoma were studied using whole transcriptome sequencing. RNA-Seq libraries were prepared by Ribo-Zero treatment of total-RNA followed by library preparation using ScriptSeq (both Epicentre/Illumina). RNA-Seq libraries were paired-end sequenced (2x 101 bp) on Illumina HiSeq 2000 and the resulting fastq files were processed using tools from the Genome Analysis Toolkit (GATK and from the Tuxedo suite. Access to the sequence data (mapped and un-mapped bam and vcf files), containing person identifying information, needs signature on a controlled access form, and can be accessed at The European Genome-phenome Archive (EGA) following request. An expression matrice of FPKM values are available without restriction at ArrayExpress (E-MTAB-4321).

Study EGAS00001001236

Paired diagnostic and relapse medulloblastoma sequencing

Whole exome and targeted panel sequencing in a cohort of paired diagnostic and relapse medulloblastoma samples.

Study EGAS00001007120

Cell-free DNA TAPS for early cancer detection

Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection. However, due to technological challenges associated with DNA methylation sequencing in low input cfDNA samples, most studies have been limited by DNA damage caused by bisulfite sequencing, or the qualitative nature of enrichment-based sequencing. Recently, we developed TET-assisted Pyridine Borane Sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and non-cancer controls. From just 10 ng cfDNA (1-3 mL of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

Study EGAS00001004962

DNA methylation and somatic mutations converge on cell cycle and define similar evolutionary histories in brain tumors

The evolutionary history of tumor cell populations can be reconstructed from patterns of genetic alterations. In contrast to the stability of genetic events, epigenetic states are reversible and sensitive to the microenvironment, prompting the question whether epigenetics can be similarly used to discover tumor phylogeny. We examined the spatial and temporal dynamics of DNA methylation in a cohort of low-grade gliomas and their patient-matched recurrences. Genes upregulated through promoter hypomethylation during malignant progression to high-grade glioblastoma were enriched in cell cycle function, evolving in parallel with genetic alterations that deregulate the G1/S cell cycle checkpoint. Moreover, phyloepigenetic relationships robustly recapitulated phylogenetic patterns inferred from somatic mutations. These findings highlight widespread co-dependency of genetic and epigenetic events throughout brain tumor evolution.

Study EGAS00001001255

Institut Curie Neuroblastoma Whole Genome Sequencing Diagnosis Relapse

Neuroblastoma, a clinically heterogeneous pediatric cancer, is characterized by distinct genomic profiles but few recurrent mutations. As neuroblastoma is expected to have high degree of genetic heterogeneity, study of neuroblastoma's clonal evolution with deep coverage whole-genome sequencing of diagnosis and relapse samples will lead to a better understanding of the molecular events associated with relapse. Samples were included in this study if sufficient DNA from constitutional, diagnosis and relapse tumors was available for WGS. Whole genome sequencing was performed on trios (constitutional, diagnose and relapse DNA) from eight patients using Illumina Hi-seq2500 leading to paired-ends (PE) 90x90 for 6 of them and 100x100 for two. Expected coverage for sample NB0175 100x100bp was 30X for tumor and constitutional samples. For the seven other patients expected coverage was 80X for tumor samples with PE 100x100, 100X in the other tumor samples and 50X for all constitutional samples (see table 1). Following alignment with BWA (Li et al., Oxford J, 2009 Jul) allowing up to 4% of mismatches, bam files were cleaned up according to the Genome Analysis Toolkit (GATK) recommendations (Van der Auwera et al., Current Protocols in Bioinformatics, 2013, picard-1.45, GenomeAnalysisTK-2.2-16). Variant calling was performed in parallel using 3 variant callers: GenomeAnalysisTK-2.2-16, Samtools-0.1.18 and MuTect-1.1.4 (McKenna et al., Genome Res, 2010; Li et al., Oxford J, 2009 Aug; Cibulskis et al., Nature, 2013). Annovar-v2012-10-23 with cosmic-v64 and dbsnp-v137 were used for the annotation and RefSeq for the structural annotation. For GATK and Samtools, single nucleotide variants (SNVs) with a quality under 30, a depth of coverage under 6 or with less than 2 reads supporting the variant were filter out. MuTect with parameters following GATK and Samtools thresholds have been used to filter our irrelevant variants. .SNVs within and around exons of coding genes overlapping splice sites.. Then,variants reported in more than 1% of the population in the 1000 genomes (1000gAprl_2012) or Exome Sequencing Project (ESP6500) have been discarded in order to filter polymorphisms. Finally, synonymous variants were filtered out. MuTect focuses on somatic by filtering with constitutional sample. Mpileup comparison between constitutional and somatic DNAs allowed us to focus also on tumor specific SNVs with GATK and Samtools. Finally, every SNV called by our pipeline and also supported in any constitutional samples were filtered our in order to prevent putative constitutional DNA coverage deficiency. Then we analyzed CNVs (copy number variants) with HMMcopy-v0.1.1 (Gavin et al., Genome Res, 2012) and control-FREEC-v6.7 (Boeva et al., Bioinformatics 2011) with a respective window of 2000bp and 1000 bp, and auto-correction of normal contamination of tumor samples for Control-FREEC. Finally we explored Structural variants (SVs) including deletions, inversions, tandem duplications and translocations using DELLY-v0.5.5 with standard parameters (Rausch et al., Oxford J, 2012). In tumors, at least 10 supporting reads were required to make a call and 5 supporting reads for the sample NB0175 with a coverage of only 40X (see table 2). To predict SVs in constitutional samples for subsequent somatic filtering, only 2 supporting reads were required in order not to miss one. To identify somatic events, all the SVs in each normal sample were first flanked by 500 bp in both directions and any SVs called in a tumor sample which was in the combined flanked regions of respective normal sample was removed (see graph 1). Deletions with more than 5 genes impacted or larger than 1Mb and inversions or tandem duplications covering more than 4 genes, were removed. We focused on exonic and splicing events for deletions, inversions, and tandem duplications. For translocation, we keep all SVs that occurred in intronic, exonic, 5'UTR, upstream or splicing regions. Bioinformatics detection of variations with Deep sequencing approach Once PE reads merged and adaptors trimmed by SeqPrep with default parameters, merged reads were aligned via the BWA (Li H. and Durbin R. 2009 PMID 19451168) allowing up to 1 differences in the 22-base-long seeds and reporting only unique alignments. Only reads having a mapping quality 20 or more have been further analysed. Variant calling software was not used, since we aimed to predict variations at low frequencies, observed in less than 1% of reads. Such variants require a custom approach. Using DepthOfCoverage functions of the Genome Analysis Toolkit (GATK) v2.13.2 (McKenna A, et al., 2010 Genome Research PMID: 20644199), we focused on high quality coverage of bases A, C, G and T at the targeted variant position. Depth of coverage of each base following a mapping quality higher than 20 and a base quality higher than 10 have been taken into account in order to focus only on high quality data. Aiming to determine the background level of variability at the studied regions, 10 control samples were included in the analysis. The same approach and filtering criteria have been applied as introduced above over the entire amplicons. In order to highlight variants, for each sample the frequencies of each bases at each amplicon position were then compared to those observed in the set of controls. Statistical analyses were performed with the R statistical software (http://www.R-project.org). Fisher’s exact two-sided tests with a Bonferroni correction were performed to compare percentages of bases between the data sets, i.e. for a given base between a case and the controls. Finally, significant variations were filtered-in once (i) a significant increase in the percentage of avariant base and (ii) a significant decrease in the percentage of it's reference base following our p.values criteria was observed (p.val < 0.05).

Study EGAS00001001184

Single cell RNA sequencing of lung adenocarcinoma.

We performed single cell RNA sequencing (scRNA-seq) for 208,506 cells derived from 58 lung adenocarcinomas from 44 patients, which covers primary tumour, lymph node and brain metastases, and pleural effusion in addition to normal lung tissues and lymph nodes.

Study EGAS00001003681

Longitudinal Single-Cell Profiling Reveals Molecular Heterogeneity and Tumor-Immune Evolution in Refractory Mantle Cell Lymphoma

The mechanisms driving therapeutic resistance and poor outcomes of mantle cell lymphoma (MCL) are incompletely understood. We characterize the cellular and molecular heterogeneity within and across patients and delineate the dynamic evolution of tumor and immune cell compartments at single cell resolution in longitudinal specimens from ibrutinib-sensitive patients and non-responders. Temporal activation of multiple cancer hallmark pathways and acquisition of 17q are observed in a refractory MCL. Multi-platform validation is performed at genomic and cellular levels in PDX models and larger patient cohorts. We demonstrate that the molecular targeting of BIRC5/survivin, locates at 17q and upregulates in resistant MCL tumor cells, results in marked tumor inhibition in preclinical models. In addition, we discovere notable differences in the tumor microenvironment including progressive dampening of CD8+ T cells and aberrant cell-to-cell communication networks in refractory MCLs. This study reveales diverse and dynamic tumor and immune programs underlying therapy resistance in MCL.

Study EGAS00001005019

Search for new loci and low-frequency variants influencing glioma risk by exome-array analysis

To identify protein altering variants (PAVs) for glioma we analysed Illumina HumanExome BeadChip exome array data on 1,882 glioma cases and 8,079 controls from three independent European populations. In addition to single variant tests we incorporated information on the predicted functional consequences of PAVs and analysed sets of genes with a higher likelihood of having a role in glioma on the basis of the profile of somatic mutations documented by large-scale sequencing initiatives. Globally tThere was a strong relationship between effect size and SNPs predicted to be damaging (P=2.29x10-49); however, these variants which are most likely to impact on risk, are rare (MAF<5%). While no single variant showed an association which was statistically significant at the genomewide threshold a number of represented promising Notable associations - were seen with the BRCA2:c.9976A>T, p.(Lys3326Ter)BRCA2 p.Lys332X variant, which has been shown documented to influence breast and lung cancer risks (OR=2.3, P=4.00x10-4 for glioblastoma [GBM]) and IDH2:c.782G>A, p.(Arg261His)IDH2 p.Arg261His (OR=3.21, P=7.67x10-3, for non-GBM). Additionally, gene burden tests revealed a statistically significant association for HARS2 and risk of GBM (P=2.20x10-6). Genome scans of low frequency PAVs represent a complementary strategy to identify disease-causing variants compared with scans based on tagSNPs. Strategies to lessen the multiple testing burden by restricting analysis to PAVs with higher priors affords an opportunity to maximise study power.

Study EGAS00001001258

Kalirin-RAC controls nucleokinetic migration in ADRN-type neuroblastoma

The migrational propensity of neuroblastoma is affected by cell identity, but the mechanisms behind the divergence remain unknown. Using RNAi and time-lapse imaging, we show that ADRN-type NB cells exhibit RAC1- and kalirin-dependent nucleokinetic (NUC) migration that relies on several integral components of neuronal migration. Inhibition of NUC migration by RAC1 and kalirin-GEF1 inhibitors occurs without hampering cell proliferation and ADRN identity. Using three clinically-relevant expression dichotomies, we reveal that most of up-regulated mRNAs in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells are associated with low-risk characteristics. The computational analysis shows that, in a context of overall gene set poverty, the upregulomes in RAC1- and kalirin-GEF1-suppressed ADRN-type cells are a batch of AU-rich elements (ARE)-containing mRNAs, which suggests a link between NUC migration and mRNA stability. GSEA-based search for vulnerabilities reveals prospective weak points in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells, including activities of H3K27- and DNA methyltransferases. Altogether, these data support the introduction of NUC inhibitors into cancer treatment research.

Study EGAS00001005023

Exploring high-throughput drug sensitivity testing in neuroblastoma cell lines and patient-derived tumor organoids in the era of precision medicine

Despite druggable events in 80% of neuroblastoma patients within the Princess Máxima Center precision medicine program 'iTHER', clinical uptake of treatment recommendations has been low, and the clinical impact for individual patients remains hard to predict. This stresses the need for a method integrating genomics and transcriptomics with functional approaches into therapeutic decision making. We aimed to launch an online repository integrating genomics and transcriptomics with high-throughput drug screening (HTS) of nineteen commonly used neuroblastoma cell lines and fourteen generated neuroblastoma patient-derived organoids (NBL-PDOs) to improve identification of molecularly matched therapies and support clinical uptake. Cell lines, NBL-PDOs and their parental tumors were characterized utilizing (lc)WGS, WES and RNAseq. Cells were exposed to a ~200 compound library, and viability was assessed using CellTiterGlo or MTT. Molecular and HTS results were transferred to the R2 platform (http://r2platform.com/pmc_nb_drugs/), in which a comprehensive suite of visualizations and analysis options were implemented. As a result, a powerful reference set of cell lines is available, reflecting distinct known pharmacologic vulnerabilities, including resistance of TP53-inactivated lines to MDM2 inhibitor idasanutlin; ALK-inhibitor sensitivity in ALK-activated lines, and sensitivity of an ATM-deleted cell line to PARP- and CHEK1-inhibitors. Novel established NBL-PDOs retained molecular features of their parental tumor and again reflected established drug sensitivities. However, HTS also identified additional therapeutic vulnerabilities in vitro, such as a striking correlation between a positive mesenchymal signature and sensitivity to BCL2-inhibitor venetoclax. Finally, we explored personalized drug sensitivities within iTHER, demonstrating HTS can support genomic and transcriptomic results, thereby strengthening the rationale for clinical uptake. Furthermore, in vitro insensitivity may avoid ineffective treatments. Lastly, novel treatment options may be identified, as indicated by sensitivity to ponatinib correlating with overexpression of RET or KIT. In conclusion, we established a dynamic publicly available dataset with detailed genomic, transcriptomic, and pharmacological annotation of classical neuroblastoma cell lines as well as NBL-PDOs, representing the heterogeneous landscape of neuroblastoma. We anticipate that in vitro drug screening will be complementary to genomic-guided precision medicine by supporting clinical decision making, thereby improving prognosis for all neuroblastoma patients in the future.

Study EGAS00001007160

4774 results for "*oma"

in 31.36 milliseconds.