PLCG1 R707Q mutation is counter selected under targeted therapy in a patient with a hepatic angiosarcoma
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
Hepatocellular-cholangiocarcinoma (H-ChC) is a rare subtype of liver cancer with clinicopathological features of both hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). To date, molecular mechanisms underlying the co-existence of HCC and iCCA components in a single tumor remain elusive. Here we show that H-ChC samples contain substantial private mutations from WES analyses, ranging from 33.1% to 86.4%, indicative of substantive intratumor heterogeneity (ITH). However, on the other hand, numerous ubiquitous mutations shared by HCC and iCCA suggest the monoclonal origin of H-ChC. Mutated genes identify herein e.g. VCAN, ACVR2A and FCGBP are speculated to contribute to distinct differentiation of HCC and iCCA within H-ChC. Moreover, immunohistochemistry demonstrate that EpCAM is highly expressed in 80% of H-ChC implying the stemness of such liver cancer. In summary, our data highlight the monoclonal origin and stemness of H-ChC, as well as substantial intratumoral heterogeneity.
Massively parallel sequencing has revolutionized research in cancer genetics and genomics and enhanced our understanding of natural human genetic variation. Recently, Lam et al. have performed a detailed comparison of two next-generation sequencing technologies with respect to their sensitivity to call single nucleotide variants (SNV) and indels. Here, we sequenced two tumor/normal pairs obtained from two paedriatic medulloblastoma patients with Life Technologies’ SOLiD 4 and 5500xl SOLiD, Illumina’s HiSeq2000, and Complete Genomics’ technology. We then compared their ability to call SNVs with high confidence. As gold standard for SNV calling, we used genotypes determined by an Affymetrix SNP array. Additionally, we performed a detailed analysis of how evenly each technology covers the genome and how the reads are distributed across functional genomic regions. Finally, we studied how a combination of data from different technologies might help to overcome the limitations in SNV calling by any of the four technologies alone.
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.
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.
Over 50 types of cancer acquire TERT promoter mutations. These single point mutations reactivate telomerase, allowing for indefinite maintenance of telomere length and enabling cellular immortalization. The transcription factor binding site created by the point mutations specifically recruit the ETS factor GABP, a multimeric transcription factor composed of the GABPα and GABPβ subunits. GABP can form two functionally independent transcription factor species – a dimer or a tetramer – depending on which of the structurally distinct GABPβ isoforms is incorporated into the complex. In this study, we show that genetic disruption of GABPβ1L, a tetramer forming isoform of GABPβ that is dispensable in normal development, results in TERT silencing in a TERT promoter mutation dependent manner. Failure to activate TERT expression by GABPβ1L culminates in telomere dysfunction, DNA damage, and mitotic cell death exclusively in TERT promoter mutant cells. Furthermore, exogenous expression of TERT is sufficient to prevent telomere degradation and loss of cell viability in GABPβ1L-reduced lines bearing the mutant TERT promoter. Orthotopic injection of tetramer-deficient mutant TERT promoter GBM cells rendered lower tumor burden and prolonged the overall survival of the tumor-bearing mice. These results highlight the potentially widespread role of GABPB1L in enabling replicative immortality of TERT promoter cancers.
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).
