{"paper_id":"0d4e5189-680d-4735-bb00-021dd754341b","body_text":"Few of the extensive number of mutations compiled for human tumors, since the advent of exome sequencing, have impacted our understanding of tumor biology like those in human diffuse low grade gliomas. These tumors represent a unique cancer subgroup largely because of a distinctive subset of variants, which includes isocitrate dehydrogenase 1/2 ( IDH1/IDH2 ) [ 1 – 3 ]. These variants create a neomorphic enzyme that generates a metabolite, 2-hydroxyglutarate (2-HG), not normally detectable in the human brain [ 4 ]. The metabolite disrupts function of proteins involved in epigenetic modification and interrupts the differentiation process [ 5 ,  6 ]. Examination of other tumor types has revealed a role for  IDH1/2  variants in the development of these cancers as well [ 7 – 9 ]. Yet, there is a clear need to further delineate the significance of these mutations with regard to the biology of tumors in general and particularly in the clinic.\nDespite a less aggressive clinical course than other malignant gliomas, such as glioblastoma, the diffuse glioma subtype, oligodendroglioma, remains an incurable disease with a range of survival times after the primary diagnosis. Oligodendroglioma exhibits a high frequency of  IDH1/2  variants, particularly in cases harboring the signature genetic event, 1p/19q codeletion [ 3 ]. The putative targets of 1p/19q codeletion, tumor suppressor genes  FUBP1  and  CIC , along with  TERT  promoter and  IDH1/2  variants, generate a unique constellation of genetic aberrations which distinguish these tumors not only from other gliomas, but from most other human cancers [ 10 – 12 ]. However, our understanding of the biology of these tumors and the molecular mechanisms at play during therapy are extremely limited due until recently to the lack of appropriate models that reflect the human disease. Previously, we orthotopically engrafted a human oligodendroglioma in two successive implantations [ 13 ]. Other  in vivo  models of oligodendroglioma have been described, but only recently has a more detailed analysis of the underlying genetic alterations been reported [ 14 – 17 ].\nIn this work, a complete analysis of the exome sequence dataset from a recurrent primary oligodendroglioma (recurrent O 2010 ) and the derived xenograft (xenograft 1 ) was performed. The single nucleotide variants (SNVs) and insertions/deletions (InDels) detected were subsequently analyzed in DNA from the tumor at the primary diagnosis in 2005 (primary O 2005 ) and a second-generation xenograft (xenograft 2 ). A subset of potential driver mutations was identified, which provides insight into the initiating events in oligodendroglioma as well as a candidate tumor suppressor pathway on chromosome 4.\n\nRecurrent O 2010  and the derived xenograft (xenograft 1 ) had been previously analyzed for gross chromosomal aberrations by array comparative genomic hybridization (CGH) [ 13 ]. Array CGH revealed the characteristic genomic profiles for oligodendroglioma, which included the 1p/19q codeletion as well as chromosome 4, 6q and 14q losses that are common in higher grade oligodendrogliomas [ 18 ]. The tumor genome at a gross level appeared to be highly stable upon passage through mice, and an increase in ploidy for chromosomes 1 and 19 was not evident based on fluorescence  in situ  hybridization analysis (Miletic, H., data not shown, [ 19 ]). One major difference was a gain of the chromosomal arm 4q in xenograft 1 , which could be the result of the outgrowth of a subclone that was either present but previously undetectable in recurrent O 2010 , or it represented the generation of a new clone altogether in the mouse.\nA more cryptic form of genomic instability, copy neutral loss of heterozygosity (cnLOH), has been previously described in oligodendrogliomas. In cnLOH, the chromosomes appeared to be diploid but were homozygous, indicating loss with subsequent duplication of a remaining chromosome. Because the predominant clone in recurrent O 2010  had lost 4q, a subclone would only be distinguishable if it retained heterozygosity at 4q. Sequencing datasets for recurrent O 2010  and xenograft 1  were therefore analyzed for copy number and regions of heterozygosity and homozygosity. Analysis of SNVs revealed regions of homozygosity correlating with chromosomal losses on 1p, 19q, 4, 6q, and 14q in recurrent O 2010  as expected [ 13 ]. LOH was also observed on chromosomes with two copies, such as 9p and 12. The cnLOH at 9p was previously observed in oligodendroglioma and was associated with loss of  CDKN2A  expression in approximately 30% of the cases examined [ 20 ].\nLOH mapped similarly in the xenograft across the genome, including the cnLOH at 9p and 12. Genomic sequencing also revealed cnLOH on 4q indicating that loss of one chromosome had occurred with subsequent duplication of the remaining arm ( Figure 1 ; [ 21 ]). Thus whether the xenograft represented a subclone from the original tumor or a new variant that developed during the engraftment process cannot be differentiated. However, these regions of cnLOH were consistent with an underlying level of genomic instability that may promote tumor development through loss of wild type and subsequent duplication of chromosomal regions harboring mutated tumor suppressor genes or oncogenes.\nDeletion of chromosome 4 in recurrent O 2010  ( A ) and copy neutral LOH 4q in xenograft 1  ( B ). LOH and CNV analysis was performed using the CEQer software [ 21 ] on the whole exome sequencing data from both tumors and data from patient blood as the control. Shown is the Log2 ratio of LOH/AI data according to the formula: ratio = MaxCovAllele / (CovAllele1 + CovAllele2), whereas MaxCovAllele is the allele with the highest coverage. Small green dots indicate copy loss of CNV exons, and horizontal green bars highlight individual copy loss of CNV regions. Small yellow dots show copy neutral regions. Big blue circles indicate either a loss of heterozygosity (-1 copy loss), an allelic imbalance (copy gain) or a copy neutral imbalance. Big grey circles indicate a conserved heterozygosity.\nWe performed in-depth analysis of the exome sequencing datasets to compare the genomes of the primary tumor/xenograft pair (recurrent O 2010 /xenograft 1 ) and to identify putative genes targeted by losses of chromosomes 4q, 6q, and 14q which typically occur in higher grade oligodendrogliomas [ 18 ].\nExome sequencing revealed 80 somatic variants, 73 single nucleotide variants (SNVs) and 7 small insertions and deletions (InDels) with potential effects on protein function in recurrent O 2010 . This number of variants ranks in a higher range for oligodendroglioma, but is in overall agreement with the number observed in tumors previously treated with radiation and chemotherapy, as in this case [ 10 ,  22 ,  23 ]. Of these, 71 SNVs and all 7 InDels were also present in the xenograft. The variants found in both included  IDH1  ( NM_005896 : c.395G>A;  R132H ),  CIC  ( NM_015125 : c.4421T>G; p.Val1474Gly - missense), and  FUBP1  ( NM_003902 :c.1307_1310delTAGA; p.Ile436fs), as previously reported.  FUBP1  on chromosome 1p sustained a small deletion, which created a truncation/frameshift [ 13 ]. Variants in tumor associated genes not typically associated with oligodendroglioma were also found, such as missense variants in  TET1  ( NM_030625 :c.4399G>A; p.Glu1467Lys) and  TGFBR2  ( NM_003242 :c.426G>C; p.Glu142Asp as well as a frameshift variant in  ARID1A  ( NM_006015 :c. 2290_2291insC; p.Gln766fs), a gene that is also located at 1p36 and frequently mutated in endometriosis associated ovarian cancers [ 24 ]. Variants in several genes were observed on chromosomes exhibiting homozygosity and/or loss, including  COL25A1 ,  HS3ST1 ,  NDST4 , and  ZCCHC4  on chromosome 4 and  ATXN2 ,  NUAK1 ,  RASAL1  and  ZC3H12  on chromosome 12. However, no variant or amplification was found in  PDGFRα  on 4q, a gene involved in the development of human oligodendroglioma, or in  CDKN2A  on 9p, one of the most frequently altered genes in diverse cancer types.\nOnly 7 novel variants were found in the xenografts. These variants could potentially be newly acquired, or they were either undetectable by our methods or in the tumor sections available. All together, these results indicate that the xenograft was genomically highly stable and reflects the genetic make-up of recurrent O 2010 .\nAs recurrent O 2010  had been treated, it was of interest to identify variants present before therapy, which could help to define new driver initiating  versus  progression or passenger mutations. In addition, whether variants in the xenograft were already present either in primary O 2005  or recurrent O 2010  was of interest in terms of potential selective advantages during engraftment.\nDeep sequencing was performed for all variants with DNAs from the patient’s blood and all tumors. Deep sequencing confirmed the presence of 81/87 variants in recurrent O 2010 , while 6 variants identified originally in xenograft 1  (out of the 7 by exome sequencing exclusively in xenograft 1  detected ones) were not present in recurrent O 2010  ( Table 1 ), also not at low frequency ( Figure 2 ).\nTimeline for the acquisition of all genetic events in primary O 2005 , recurrent O 2010 , and xenografts 1&2 , including array CGH. Exome sequencing was performed on recurrent O 2010  and xenograft 1 , and these datasets were used to generate the set of variants (total = 87) to be further analyzed by deep sequencing. Deep sequencing detected a subset of these variants in primary O 2005  ( n  = 11) and was used to generate an allele frequency for each variant in each tumor. Only 3/11 (red) variants, including  IDH1 R132H  ( NM_005896 :c.395G>A), were considered to be in the majority of the cells in primary O 2005  (allele frequency ≥ 50%). The other 8 (blue) were detectable in primary O 2005  but at lower frequencies (< 50%). In addition to these 11 variants, 70 novel variants occurred in the recurrent O 2010  (total = 81, 70 new, distinguished by brown, gold, and light blue). Of these, 63 were present in the majority of the cells (brown). Included are variants in  CIC  ( NM_015125 :c.4421T>G) and  FUBP1  ( NM_003902 :c.1307_1310delTAGA). A single variant,  ARHGAP6   NM_013423 :c.455G>C (light blue), was found only in recurrent O 2010 . In the xenografts, only 6 new variants (purple, ≥ 50% and pink, < 50%, perhaps latest changes) appeared in addition to those from recurrent O 2010  ( n  = 86; 6 new, 1 lost). Additional targeted Sanger sequencing revealed the presence of a variant in  TERT  promoter in all tumors.\nOnly 11/87 somatic variants were also present in the primary O 2005  ( Table 2 ), so that the other mutations ( n  = 76) were considered to be due to progression and/or treatment. All variants found in primary O 2005  were confirmed using Sanger sequencing. With the exception of one variant (one of two identified by exome sequencing),  ARHGAP6  ( NM_013423 :c.455G>C; p.Arg152Pro) on the X chromosome, all variants observed in recurrent O 2010  were confirmed in both xenograft 1  and xenograft 2 . One small InDel,  NUP210L  ( NM_207308 :c.2072delT; p.Leu691fs), was also detected in the DNA from the patient’s blood, but with an allele frequency of only 1.1%. This result is possibly due to a sequencing error.\nOf the genes already associated with the development of oligodendroglioma ( IDH1/IDH2 ,  CIC , and  FUBP1 ), remarkably only the mutation in  IDH1  ( NM_005896 :c.395G>A) was present in the primary O 2005 .  FUBP1  and  CIC  variants were not present, even though 1p/19q codeletion had been detected by array CGH in primary O 2005  (Supplementary Figure 1). Thus,  FUBP1  and  CIC  variants appeared to occur secondary to the losses of 1p and 19q during tumor progression and were not necessary changes for tumor initiation in this case.\nIn addition to the  FUBP1  frameshift variant, we also found a loss-of-function mutation of the known tumor suppressor gene  ARID1A  on 1p in the recurrent O 2010  but not in primary O 2005 . In the TCGA dataset, alterations were detected in  ARID1A  in low grade glioma (SNV: 3.9%, 11/283; deletion: 0.4%, 1/283) and glioblastoma (SNV: 0.7%, 2/281) [ 25 ,  26 ]. As this variant emerged following treatment of the patient, the relationship to tumor development remains unclear.  CIC ,  FUBP1 , and  ARID1A  variants were all maintained through xenograft 2  and in a homozygous state.\nMost variants in the  TERT  promotor occurring in human gliomas are one of two changes:  NM_198253.2 :c.-124G>A (C228T) and  NM_198253.2 : c.-146G>A (C250T) located in a common SNV  rs2853669  [ 12 ,  27 ,  28 ]. In oligodendrogliomas, a variant in the  TERT  promoter was found in 78% of cases ( n  = 45), and of these, 43% were the C250T variant [ 12 ]. To determine whether the  TERT  promoter also harbored a variant in our system, we performed targeted Sanger sequencing on DNA from each of the tumor samples. C250T ( NM_198253.2 : c.-146G>A; p.Cys250Thr) was found in all samples and was heterozygous.\nAllele frequency based on the sequencing reads was also used to more precisely reconstruct the temporal acquisition of variants ( Figures 2  and  3 ). We identified 10 variants that displayed a significantly lower frequency where the variant first appeared, than in the tumor that followed. Eight were variants that were already detectable in primary O 2005 : the missense variants in  ATP7B  ( NM_000053 :c.1829C>T; allele frequency, AF 3%);  C3orf15  ( NM_033364 :c.1820G>A; AF, 10%);  DUSP11  ( NM_003584 :c.575A>G; AF, 10%);  EPHA1  ( NM_005232 :c.1262C>A; AF 10%);  NDST4  ( NM_022569 :c.2392C>G; AF 20%);  SYPL2  ( NM_001040709 :c.392C>G; AF 20%) and the frame shift mutations in  NUP210L  ( NM_207308 :c. 2072delT; AF 1%) and  SMAD7  ( NM_001190823 :c.484_485insT; AF 10%). The remaining two variants initially appeared in recurrent O 2010  and became predominant in xenograft 1 :  ZAK  ( NM_016653 :c.1618C>T; p.Gln540*; AF 20%;  Figure 3 ) and the InDel in  KPNA1  ( NM_002264 :c.179_181delAAG; p.lu60del; AF 20%). Only two variants,  PSKH1  ( NM_006742.2 :c.650G>A; p.Arg217Gln), and  SNX12  ( NM_001256188 :c.470G>A; p.Arg157His) were present in ≥ 50% of reads as early as primary O 2005 . Variants in these genes and  IDH1  exhibited an allele frequency that remained stable across all tumor samples. They were considered to be heterozygous and uniformly distributed amongst tumor cells, indicating that they had occurred early in tumor development.\nGraphic representation of allele frequency (% of reads for variant indicated) for variants in 4 genes in the tissue sample indicated. ( A )  PSKH1 ; ( B )  NDST4 ; ( C )  CIC ; and ( D )  ZAK .\nVariants that exhibited an allele frequency of < 50% were considered to have occurred later and in a subset of tumor cells. Any increase in allele frequency could be an indication that a subclone bearing the variant had been selected for. Examples of such variants were  NDST4  ( NM_022569 :c.2392C>G; p.Leu798Val) and  ZAK  ( NM_016653 :c.1618C>T; p.Gln540*), which were first detected in primary O 2005  and recurrent O 2010 , respectively ( Figure 3 ).  NDST4  represented an interesting case because of the changing status of chromosome arm 4q. In primary O 2005 , where 4q was diploid, the mutation was detected at a frequency of ~20%. In recurrent O 2010 , the allele frequency of the  NDST4  variant was 80%. This result was consistent with selection of tumor cells harboring reduction to homozygosity at  NDST4 /chromosome arm 4q. The allele frequency of the variant increased to 100% in xenograft 1  in the absence of contaminating non-neoplastic human cells. In xenograft 1  with cnLOH on the 4q arm, a normal variant of  NDST4  was no longer present, indicating that duplication of the remaining mutated allele had occurred.\nChromosomal losses 4q and 14q have been associated with higher grade oligodendroglioma.  NDST4  was a candidate tumor suppressor gene on chromosome 4 as the  NDST4  previous variant was already detectable in the primary tumor and the wild-type allele was lost in recurrent O 2010 , indicating a loss of function mechanism. To investigate if alterations in  NDST4  might play a role in the development of oligodendroglioma, we performed Sanger sequencing to analyze all coding exons of the gene in our own cohort of primary oligodendrogliomas ( n  = 15). We found no variants in coding regions of  NDST4  in this tumor set; only known polymorphisms and variants in intronic sequences were observed. However, only 3 cases in our cohort had sustained loss of chromosome 4. In TCGA datasets, 8 missense mutations and deletions in  NDST4  were found in lower grade gliomas, while 17 missense and 1 truncating variant were found in glioblastoma [ 25 ,  26 ]. All together the somatic mutation frequency was 0.6% [ 25 ,  26 ]. Thus, although  NDST4  remains a gene of interest in tumor development and indeed might play a role in the progression from primary O 2005  to recurrent O 2010 ,  NDST4  mutation seems to be a rare event in oligodendroglioma.\n\nExome sequencing coupled with array CGH performed on a rare oligodendroglioma xenograft revealed a highly stable genome relative to the recurrent human tumor from which it was derived. Using deep sequencing (and targeted Sanger sequencing for the  TERT  promoter), we identified a set of early and potential driver variants ( n  = 12) already present in genes in primary O 2005  which included  IDH1  and  TERT , but surprisingly not  FUBP1  and  CIC . Allele frequencies for 10 of these 12 variants was ≤ 50%, but increased for some variants indicating possible selection of tumor cell subpopulations in response to normal progression, treatment, or engraftment. Finally, the increase in allele frequency of the  NDST4  variant in the tumor series was consistent with the selection of a tumor suppressor gene.\nThe overall number of variants ( n  = 82 - including  TERT ) in recurrent O 2010  was significantly higher compared to primary O 2005  ( n  = 12 - including  TERT ) and the average for primary oligodendrogliomas, but was consistent with cases from patients who had undergone treatment [ 10 ,  11 ]. All variants were also detected in xenograft 1 , along with 6 new changes including 5 SNVs and 1 small InDel. These 6 changes in xenograft 1  might have spontaneously occurred during progression in xenograft 1  or were previously undetectable but selected for during engraftment into the mouse. Only one variant in  ARHGAP6  was exclusively found in recurrent O 2010  but at a low allele frequency (12%).\nThrough the deep sequencing analysis emerged a sequence of mutational events that may be more generally relevant to the development of a subset of human oligodendrogliomas.  FUBP1  and  CIC  variants first appeared in recurrent O 2010 . Thus, the scenario for inactivation of these tumor suppressor genes in some cases may be loss of heterozygosity before selection of cells with variants. It has been previously suggested that haploinsufficiency at 1p is tumor promoting [ 29 ]. In oligodendrogliomas, haploinsufficiency occurs with a high frequency in the context of an  IDH1/IDH2  variant [ 3 ]. In such cases, loss in combination with epigenetic modification may sufficiently knock down expression of key target genes to promote uncontrolled proliferation and/or inhibit differentiation. Variants in  FUBP1  and  CIC  may further stabilize the phenotype established by haploinsufficiency along with epigenetic changes.\nIn support of this scenario, results from several studies indicate that concurrent variants in  FUBP1  and  CIC  appear infrequently. In two early studies, only ~10% of human oligodendrogliomas with codeletion of 1p and 19q were found to harbor simultaneous variants [ 10 ,  22 ]. Our own study from 2013 corroborated these results.  FUBP1  and  CIC  were both mutated in only 12% (2/17) of oligodendrogliomas examined, and all had sustained a variant in  IDH1  [ 30 ]. In the dataset from the TCGA, simultaneous variants occurred in a higher proportion of tumors, 24% (13/53) [ 25 ,  26 ]. However, while we used deep sequencing in our analysis, there are other possibilities that might account for the low frequency of simultaneous mutations in primary tumors, such as technique and sample bias due to tumor heterogeneity or small undetectable subclones.\nAlthough ultimately it is the accumulation of genetic events that results in cancer, the putative timeline of our oligodendroglioma might provide insight for the contribution of  IDH ,  TERT  and 1p/19 codeletion to the development of oligodendroglioma. First, the results might support studies indicating that  IDH  mutation occurs early in a stem/progenitor cell type with subsequent  p53  mutations or 1p/19 codeletion defining the path to astrocytoma or oligodendroglioma.  IDH  mutation appears to be a gatekeeper, triggering uncontrolled proliferation which renders cells vulnerable to additional tumor promoting genetic events [ 31 ]. Second, 1p/19q codeletion might drive the development of oligodendroglioma through mechanisms other than loss of  FUBP1  and  CIC  [ 32 ]. Examination of expression datasets for oligodendroglioma have led to the identification of a network of genes in the codeleted region of 1p/19q that may promote the development of oligodendroglioma. Finally, we can wonder whether  TERT  or  IDH  mutations have a more direct impact on genomic instability.  TERT  promoter mutations have been shown to extend cellular lifespan but do not necessarily prevent shortening of telomeres [ 33 ]. These shorter telomeres however later in tumor development contribute to genomic instability. A second possibility is mutation of the  IDH  gene itself.  IDH  mutation promotes epigenetic modification of DNA, and methylation has been implicated in chromosomal instability, although more generally through global hypomethylation [ 34 ]. The fact that the incidence of  IDH  mutation is so high in 1p/19q codeleted oligodendrogliomas and that cnLOH is also commonly observed further supports this possible role for  IDH  mutation in cancer [ 3 ,  20 ].\nFinally, deep sequencing performed on primary O 2005  importantly also enabled us to distinguish potential driver from passenger variants in recurrent O 2010 .  NDST4  on chromosome 4q emerged as a potential tumor suppressor gene as it was one of the few genes with a variant already detectable in primary O 2005 . Unlike  IDH1  or  PSKH1 , this variant was probably present only in a subpopulation of cells in primary O 2005  which later became predominant. The protein is involved in heparan sulfate synthesis, and the distribution of heparan sulfate is altered in tumorigenesis. Heparan sulfate synthesis has also been shown to be necessary for growth factor activity. Furthermore, heparan sulfates contribute to cell-cell and cell-extracellular matrix interactions forming barriers for tumor invasion [ 35 ]. Interestingly, two additional genes mutated in recurrent O 2010 ,  EXTL3  (chromosome 8) and  HS3ST1  (chromosome 4), are also involved in heparan sulfate synthesis. These genes are altered in many cancers and like  NDST4 , frequent alterations are present in skin cutaneous melanoma and lung cancer [ 25 ,  26 ]. Therefore, alteration of heparan sulfate synthesis could be important in tumor biology.\nIn conclusion, we have used exome sequencing datasets from a unique series of a recurrent human tumor and the derived xenografts to construct an evolutionary genetic profile. Although the primary oligodendroglioma exhibited 1p/19q codeletion, the variants detected did not include changes in  FUBP1  and  CIC . Only a few variants, in addition to  IDH1  R132H and  TERT  C250T, were observed in the primary tumor that gained predominance and might represent interesting candidate genes for oligodendroglioma development. Because mutations in several of these genes are rare in oligodendrogliomas, and human cancer in general, some of these variants may simply represent hits in the so-called “hills” of the oligodendroglioma genome landscape [ 31 ]. Functional experiments are therefore required to determine whether the genes identified in this study are in fact true driver genes in human cancer.\n\nPatient material was obtained from surgeries performed at the Haukeland University Hospital (Bergen, Norway). Written consent was obtained from all patients to use tissues and clinical data for research purposes. Procedures were approved for the project (project number 013.09) by the Regional Ethics Committee (Bergen, Norway). All animal protocols were approved by authorities in an AAALAC accredited facility at the Haukeland University Hospital and in accordance with the national regulations of Norway (project numbers 2010 2658 and 2011 3079).\nTissue was obtained from a recurrent primary tumor from 2010 (oligodendroglioma WHO grade III; recurrent O 2010 ), and implanted as previously described [ 13 ]. Xenograft tissue was harvested at the first implantation to obtain the first-generation xenograft (xenograft 1 ). Cells from this tumor were reimplanted to obtain the second-generation xenograft (xenograft 2 ). Portions of all three tumors were harvested and stored at –80°C for molecular analysis. Tissue from the original tumor (oligodendroglioma WHO grade III; primary O 2005 ) was obtained from archived paraffin embedded surgical material in 2005.\nDNAs were isolated from patient material that had been immediately frozen in liquid nitrogen subsequent to surgery. Frozen sections were prepared, and every fifteenth section was stained with H&E to ensure > 60% tumor cell content. As a control, DNA was obtained from frozen patient blood (–20°C) by isolating nuclei with Buffer C1 from the Genomic DNA Buffers Kit (Qiagen; Germantown, MD, USA). Samples were treated with proteinase K overnight at 50°C in ATL buffer (Qiagen) and then DNAse free RNAseA (Fermentas/Thermo Fisher Scientific; Waltham, MA, USA) for 5 min at room temperature. Reactions were extracted with phenol:chloroform:isoamyl alcohol 25:24:1 saturated with 10 mM tris, pH 8.0, 1 mM EDTA (Sigma; St. Louis, MO, USA) and precipitated in 2.5M NH 4 OAc and 2.5 volumes of 100% ethanol. DNAs were resuspended in Nuclease-Free water (Qiagen). DNA was extracted from FFPE tissue samples using the QIAamp DNA FFPE Tissue Kit (Qiagen) according to the manufacturer’s protocol.\nAgilent SureSelect target enrichment kit (SureSelect Human All Exon 50 Mb kit, Agilent; Santa Clara, CA, USA) was used for exome library construction according to the manufacturer’s protocol. Captured DNA libraries were sequenced on the Illumina HiSeq2000 platform at a coverage of 100× (Illumina; San Diego, CA, USA).\nBriefly, the genomes derived from human tissues were aligned against the human genome hg19. Analysis of the sequencing reads from the xenografts underwent a cleaning process to eliminate contaminating mouse reads. Two alignment steps were performed. In the first step, the sequencing reads were aligned against the human (hg19) and the mouse (mm10) genomes separately. Reads with better mapping results (higher mapping quality) in the mouse genome than in the human genome were deleted. In the second step, surviving reads were aligned against a chimeric genome, which was obtained through the concatenation of the human and mouse genomes. The resulting BAM file was then separated into those reads that mapped to the human or the mouse genome using SAMtools (version 0.1.18) [ 36 ,  37 ]. The cleaning process was therefore accomplished by cross-referencing reads with a high mapping quality in the other set and eliminating them.\nThe pipeline for all the samples was the following: alignment using BWA (version 0.7.2) [ 38 ], and post-alignment processing using GATK [ 39 ] for local realignment, duplicate removal and score recalibration, and finally variant calling between each pair of samples using VarScan (version 2.3.2) [ 40 ]. SAMtools (version 0.1.18) was used to create the needed pileup file [ 36 ,  37 ]. False positives were eliminated by a false positive filtering (strandedness (0.1); maximum mismatch quality sum (100); minimum average distance to effective 3’end (0.1) [ 40 ].\nAfter a statistical filtering for somatic variants (for details see Supplementary Material), we selected for mutation type (non-synonymous substitution, start and stop codon gain or loss, and frame-shifts), to focus on variants with potential effects at the protein level. All called variants were hand-curated, which included an evaluation using the Integrative Genomics Viewer provided by the Broad Institute [ 41 ,  42 ], applying the BLAT - BLAST- like Alignment Tool [ 43 ] of the UCSC Genome Browser [ 44 ], and including additional information from the UCSC Genome Browser, such as pseudogenes and self-chain or repetitive regions. Variants were determined to be damaging/deleterious and disease causing by several prediction software programs, including PolyPhen-2, Mutation Taster, Mutation Assessor, PROVEAN and SIFT [ 45 – 51 ]. Details of the variant selection process are described in Supplementary Materials.\nAll variants were additionally verified using Sanger and/or deep sequencing. Primers were designed against regions of interest using Primer3 software (version 4.1.0) [ 52 ,  53 ]. Primers and conditions used are shown in Supplementary Data. A first/preliminary analysis of the xenograft sequencing data revealed a large number of novel variants that were suspected to be due to contamination with mouse tissue rather than novel variants occurring during engraftment. Therefore, for validation of variants, we designed all primers to bind to human sequences only. Validation of primers was performed on normal human control DNA obtained from blood and DNA from mouse tissue to verify that primers would only yield a PCR-product in human but not mouse DNA.\nFor ultra-deep sequencing, PCR was carried out using a proofreading polymerase (Phusion ®  Hot Start Flex DNA Polymerase; New England Biolabs; Ipswich, MA, USA) according to the manufacturer’s instructions under the following cycling conditions: initial denaturation step: 98°C, 30 sec; amplification step (34×): 98°C, 10 sec, 61.5 – 66.5°C, 30 sec, 72°C, 30 sec; final extension: 72°C, 10 min. The correct size of PCR products was confirmed by agarose gel electrophoresis (2.0% agarose impregnated with GelRed; 1 h, 110 V) and purified from contaminating nucleotides and primers by PEG precipitation. Libraries were sequenced to an average coverage of 36,000× on an Illumina HiSeq2000 Next Generation Sequencer (Illumina; San Diego, CA, USA). Results were analyzed using CLC Genomics Workbench Version 7.5. Reads were aligned to the human reference genome (hg19) and variants were called using the Basic Variant Detection algorithm (minimum coverage 1000×; minimum count 2, minimum frequency 1%; the base quality filter was used with default values, no direction or position filters were applied).\nVariants found in primary O 2005  were validated using Sanger’s dye terminator method on an ABI 3130x1 DNA Analyzer (Applied Biosystems/Thermo Fisher Scientific). Products were cleaned with rAPid alkaline phosphatase (Roche Diagnostics GmbH; Mannheim, Germany) and exonuclease I treatment (New England Biolabs) at 37°C for 1 h followed by inactivation at 80°C for 20 min.\nArray CGH was performed on Agilent’s SurePrint G3 Human CGH Microarray Kit 2x400K (Design ID021850, Agilent, Santa Clara, CA, USA) according to the manufacturer’s protocol. Scanning of the hybridized arrays was carried out using the Agilent High-Resolution Microarray scanner. Raw data were processed by the Feature Extraction 9.5 (Agilent) software and normalization was performed using the default settings. Feature Extraction files were imported into Genomic Workbench 7.0 (Agilent Technologies) for visualization and analysis. After diploid centralization and GC correction, aberrations were called using the ADM2 algorithm with a threshold setting of 20, centralization on with threshold of 25 and an aberration filter min Probes = 5 and minAvgAbsLogRatio = 0.35 for amplifications and deletions [ 13 ].\nTCGA data was analyzed using the cbioPortal for Cancer Genomics [ 25 ,  26 ].","source_license":"CC-BY-4.0","license_restricted":false}