- Research article
- Open Access
Discovery of structural alterations in solid tumor oligodendroglioma by single molecule analysis
BMC Genomics volume 14, Article number: 505 (2013)
Solid tumors present a panoply of genomic alterations, from single base changes to the gain or loss of entire chromosomes. Although aberrations at the two extremes of this spectrum are readily defined, comprehensive discernment of the complex and disperse mutational spectrum of cancer genomes remains a significant challenge for current genome analysis platforms. In this context, high throughput, single molecule platforms like Optical Mapping offer a unique perspective.
Using measurements from large ensembles of individual DNA molecules, we have discovered genomic structural alterations in the solid tumor oligodendroglioma. Over a thousand structural variants were identified in each tumor sample, without any prior hypotheses, and often in genomic regions deemed intractable by other technologies. These findings were then validated by comprehensive comparisons to variants reported in external and internal databases, and by selected experimental corroborations. Alterations range in size from under 5 kb to hundreds of kilobases, and comprise insertions, deletions, inversions and compound events. Candidate mutations were scored at sub-genic resolution and unambiguously reveal structural details at aberrant loci.
The Optical Mapping system provides a rich description of the complex genomes of solid tumors, including sequence level aberrations, structural alterations and copy number variants that power generation of functional hypotheses for oligodendroglioma genetics.
Cancer is fundamentally a disease of genomic origin. Alterations in genes and regulatory elements critical to cell cycle control lead to uncontrolled cell growth and proliferation, the common signature of all cancers. Such events can cause amplification or mutational activation of oncogenes[1, 2], deletion or mutation deactivation of tumor suppressor genes[3, 4], orientation of genes with incorrect regulatory regions, gene fusion products, etc. As cancers evolve, they accumulate a cascade of mutations, ranging in size from a single nucleotide change to the gain or loss of entire chromosomes. Coupled with the subclonal heterogeneity that is a hallmark of solid tumors[8–10], obtaining a complete portrait of the genetic landscape of human cancer remains a significant challenge.
Synergy between revolutionary genomic tools and advances in high-throughput computing has facilitated the development of a number of methods for detecting mutations. Chromosome banding and spectral karyotyping (SKY) are low-resolution techniques used to detect large-scale chromosomal features. However, obtaining metaphase spreads for performing a karyotype is often difficult, especially when working with solid tumor biopsies and paraffin embedded, formalin fixed tissue. Fluorescence in situ hybridization (FISH) and its variants are a family of molecular cytogenetic techniques developed to correlate specific sequences to cytogenetic observations. FISH offers higher resolution (compared to SKY) and has the advantage of not requiring metaphase spreads, but is limited by the fact that it requires a prior hypotheses about the locus of interest, making it unsuitable for discovery based research. Hybridization based microarray approaches, like SNP microarrays and array comparative genome hybridization (CGH), have been extensively used to detect large scale amplifications and deletions in tumor genomes[13–15], but are unable to detect changes where there is no net gain or loss of DNA, such as inversions and balanced translocations, which have been shown to be an important mechanism for oncogenic transformation[16–18]. Moreover, microarrays do not offer structural information, necessitating follow-up experiments to identify the breakpoints and sequence context of the aberration. Microarrays are also restricted to regions of the genome amenable to unique probe design, which precludes repeat-rich regions and novel insertions that are hotbeds of variation and mutation[19–22]. Most commercial microarrays (except custom designed, high-density arrays) lose sensitivity below ~50 kb, and variants, particularly insertions, in this size range have remained largely unexplored, especially in cancer genomes.
The advent of massively parallel, short read DNA sequencing- the ‘second generation’ sequencing technologies, and their application to cancer has also accelerated the pace of mutation discovery. Initially applied to targeted subsets of the genome, such as specific gene families (e.g.: all protein kinases, or ‘kinome’)[24–27], or all the coding sequences (the exome)[28–33], second-generation sequencing is increasingly being used to interrogate whole cancer genomes[34–40]. In theory, second-generation sequencing of whole genomes has the ability to discern the full range of genomic alterations. In practice, however, more than 90% of events discovered by these platforms are less than 1 kb, and are biased towards deletions rather than insertions[23, 41]. Second-generation sequencing instruments typically generate shorter reads with higher error rates from relatively short insert libraries, which present a significant computational and bioinformatic challenge in alignment and assembly. Read-pair mapping approaches have successfully identified point mutations and indels in cancer[36, 38–40], but are limited by the insert size of the DNA library to detecting base substitutions and small indels and are often confounded by repetitive regions of the genome. Further, accurate prediction of the exact breakpoints of an aberration depends on very tight size distribution of the DNA library, which can make library construction difficult. Whole genome sequencing followed by de novo assembly might mitigate some of these issues, but current assembly algorithms tend to collapse homologous sequences, and consequently dramatically under-represent repeats and segmental duplications that are known to be critical mediators of genomic rearrangement.
There remains a pressing need for discovery-based systems that can provide a scalable, comprehensive view of the cancer genome in its entirety. In this study, we present Optical Mapping as one such system. Optical Mapping creates high-resolution ordered restriction maps of whole genomes through the analysis of ensembles of single molecule restriction maps. It has previously been used to map the genomes of microbes[45–48], plants[49–52] and mammals[53–57]. However, this is the first time it has been employed to analyze the genome of a solid tumor. Optical Mapping offers several unique advantages towards assembling the complex structure of a cancer genome. Genomic DNA isolated directly from cells is analyzed, thereby obviating any bias introduced by amplification or cloning steps. Moreover, because the DNA is of high-molecular weight (300 kb - >500 kb), segmental duplications and other repeat-rich regions of the genome are revealed, and additionally, the structure and long-range context of any aberration are determined. Since the restriction maps are made from single DNA molecules, Optical Mapping effectively pieces together heterogeneous alterations, which is especially important for tumor genome analysis, as we demonstrate in oligodendroglioma.
Oligodendrogliomas are frontal lobe tumors that are thought to arise from oligodendrocytes, supporting brain cells which provide myelination for neurons[58, 59]. The concerted loss of heterozygosity (LOH) of chromosome arms 1p and 19q, observed in 50-70% of patients, is a molecular signature of this malignancy. The remarkably high prevalence of this molecular marker suggests that these regions harbor one or more tumor suppressor genes that might play an important role in the development of the tumor. Allelic losses of 1p/19q have been correlated with positive response to chemo- and radiotherapy and prolonged survival for patients with oligodendroglioma. However, it remains unclear whether LOH of 1p/19q is a prognostic biomarker for a more indolent tumor subtype that has fewer unfavorable mutations overall, rather than predictive of treatment sensitivity[62, 63]. In fact, studies have shown that 1p/19q codeleted tumors have slower growth rates and are more responsive to treatment than tumors without the codeletion[64, 65]. In order to explore each of these possibilities, Optical Mapping was used to create physical maps from two individual oligodendroglioma tumor biopsies for the purpose of identifying and characterizing structural changes on a whole genome basis.
Results and discussion
Optical map construction
We used the Optical Mapping (OM) system to explore the genomic landscape of a solid tumor. Optical Mapping creates high-resolution physical maps of genomes through the analysis of ensembles of single molecule ordered restriction maps. The tumor biopsies were disaggregated into single cells, then run through a Percoll gradient to enrich for cancer cells (methods). High molecular weight genomic DNA was extracted directly from these cells, stretched and immobilized in regular arrays on positively charged glass surfaces using a microfluidic device (Figure 1A, details in methods). After deposition, the DNA was digested with the restriction enzyme SwaI. The surface-bound restriction fragments remained in register, and were stained with a fluorescent dye and imaged by automated fluorescent microscopy (Figure 1B). Dedicated machine vision software calculated the size, in kilobase pairs (kb), of each fragment based on measurements of integrated fluorescent intensity, resulting in the high throughput, massively parallel generation of ordered restriction maps, or ‘Rmaps’, from individual genomic DNA molecules (Figure 1C). The oligodendroglioma datasets comprise close to 700,000 such Rmaps, with an average size of greater than 400 kb (Additional file1). The relative order and distance between successive restriction fragments in a single molecule optical map can be used to determine the precise location in the genome that gave rise to that molecule, by means of pair-wise alignment against an in-silico restriction map. The HF087 tumor dataset comprised 235,026 aligned Rmaps, corresponding to 36.92 fold coverage of the human genome, while the HF1551 dataset comprised 167,012 aligned maps, representing 25.36 fold coverage of the human genome (Additional file1). In the absence of a karyotype, our assessment of ploidy is based on optical map coverage and Affymetrix array analysis. Both these platforms, discussed in detail in subsequent sections, calculate chromosome copy number relative to normal, diploid genomes, and are in agreement that neither tumor sample is polyploid. They do, however, display aneuploidy, due to allelic losses of specific chromosomes/chromosome arms (1p, 19q, 13 in HF087 and 1p, 19q, 14, 21 in HF1551), so if anything the coverage is likely to be higher than what we reported.
The Rmaps that cluster together upon pair-wise alignment were then assembled into consensus optical maps and analyzed for presence of structural variants using the bioinformatics pipeline described in. The final consensus map contigs span 96.73% and 93.92% of the human genome for tumors HF087 and HF1551, respectively.
Optical map coverage analysis
Discernment of copy number variants
Copy number was inferred from aligned coverage of Rmaps, prior to assembly, in a manner analogous to read-depth based methods for detecting copy number variants from second generation sequencing data (methods). Briefly, Rmaps were aligned to the in silico reference map, and then partitioned into discrete windows spanning each chromosome. These alignments were then compared to alignments of a reference data set (comprising of a number of normal genomes) that was used to ‘normalize’ the observed coverage. This was necessary because the number of Rmaps that align to a particular region of the genome depends, in part, on the density of restriction sites in that region, which varies from chromosome to chromosome (ranging from a low of 2.5 cuts/100 kb on chromosome 22 to a high of 9.25 cuts/100 kb on chromosome 4). A Hidden Markov Model (HMM) was then fitted to this data, and copy number changes were detected (Figure 2A)[68, 69]. Optical map coverage analysis confirmed the allelic loss of chromosome arms 1p and 19q in HF087 and HF1551. The breakpoints appear to be very close to the centromere, consistent with the proposed mechanism of an unbalanced reciprocal translocation mediating the LOH event[70, 71]. Additionally, coverage analysis also detected allelic loss of chromosome 13 (HF087), 14 and 21 (HF1551), which are known to be rarer events associated with oligodendroglioma.
Solid tumor heterogeneity
The genome wide optical map of HF1551 was created using DNA from two adjacent slices of the tumor: 446,933 (~55%) Rmaps originated from slice 1 and 202,974 (~45%) Rmaps from slice 2 (Figure 2). Interestingly, when the Rmaps were partitioned according to the slice they originated from, and coverage analysis was performed separately, unique copy number profiles were obtained for each slice. In addition to allelic losses of 1p and 19q, slice 1 also had LOH of chromosomes 14 and 21, while slice 2 had evidence of LOH of 19p (Figure 2B). Solid tumors are dynamic aggregates of continually evolving subclones, resulting in spatial and temporal genetic heterogeneity. Our findings suggest that the tumor slices used for Optical Mapping evolved from distinct cancer cell clones, and is congruent with recent evidence of branched evolutionary tumor growth[72–74]. Although assembly of whole genome maps on a per slice basis was not feasible due to insufficient depth of coverage, our results establishes proof-of-principle of Optical Mapping to interrogate tumor heterogeneity.
Discovery of optical structural alterations
The optical consensus maps generated by map assembler were aligned to the in silico restriction map (generated from build 35 human reference sequence), and by comparison of the order and sizing of the 219,224 restriction fragments (fragments smaller than 0.4 kb in size were merged) between the experimental and the reference map. Such comparisons revealed structural variants in the experimental genome that were classified as four types: extra cuts (EC), where the optical consensus map displays a restriction site that was not predicted by the reference sequence; missing cuts (MC), where a cut that was predicted was not observed in the experimental map; insertions (INS), where the size of a fragment in the consensus map was significantly larger than its counterpart in the reference map; deletions (DEL), where a fragment in the experimental map was smaller than the corresponding reference fragment (or missing altogether); and finally, complex events (OTHER) involving multiple cut or size differences (methods). Approximately a third of the ECs and MCs represent small indels that are below the resolution of Optical Mapping (~3 kb). Figure 3C shows an example of each class of variant detected by Optical Mapping.
At first glance, it might appear that any one of these variants could be attributed to errors inherent in Optical Mapping. For instance, a missing cut could be due to incomplete digestion, an extra cut could result from spurious cutting by the restriction enzyme, or physical breakage of the DNA molecule, and uneven staining could lead to inaccurate estimation of fragment size. However, the high throughput advantage of Optical Mapping allows us to distinguish such random errors from legitimate genomic events. Any alteration in the optical consensus map was supported by multiple single molecule maps (Rmaps), each representing an independent observation at that locus. The Optical Mapping error models estimated the statistical significance of each structural variant, after taking into account the quality and quantity of the data.
A total of 1081 and 1085 differences were detected in HF087 and HF1551 respectively (Figure 3A and B). The distribution of structural variants across the genome is uniform and the pattern is similar for both tumors (Figure 3A). Variants range in size from single base differences to complex genomic events spanning hundreds of kilobases (Figure 4). Approximately 800 single base changes were detected in each tumor, including point mutations (such as the example depicted in Figure 4A), polymorphic SNPs where only one allele has a SwaI restriction site (referred to as snip-SNPs), and small indels that create or remove a SwaI cut site but are below the detection limit of Optical Mapping. 179 indels with a median size of 6.6 kb were detected in each sample (Figure 5). For comparison, the median size of indels reported in the Database of Genomic Variants is 2.3 kb (Figure 5, inset). ~70 complex events were found in each tumor, including known polymorphic loci such as the major histocompatibility complex (MHC), giving us confidence that these results are not spurious. Optical Mapping also discerns balanced genomic events, where there is no net gain or loss of genomic sequence. A putative inversion spanning 352 kb of chromosome 7 was observed in HF087 which appears to disrupt the ZNF92 gene (Figure 4C). Finally, the largest events detected by Optical Mapping include gains or losses of entire arms of chromosomes, for example, the allelic loss of chromosome 1 illustrated in Figure 4D, and discussed in detail previously. Intersection counts with genes, segmental duplications, published SNPs (dbSNP build 135,http://www.ncbi.nlm.nih.gov/projects/SNP/) and published structural variants (Database of Genomic Variants, November 2010 release) are shown in Additional file2. Comprehensive breakdown of the overlaps are shown in Additional file3.
Optical Mapping provides a comprehensive description of the vast and complex landscape of cancer genomes. The ability to study the genome in its entirety, including non-genic or repetitive regions using a single technology minimizes ascertainment bias. As detailed in subsequent sections, it is employed to generate a list of candidate cancer genes that is not hypothesis-limited, and elucidate their structure at sub-genic resolution.
Validation of copy number and structural variants
Experimental validation: SNP array
The Affymetrix Genome Wide Human SNP 6.0 Array, which has probes for detection of both SNPs and copy number variants (CNVs), was used to validate our findings.
Both platforms concurred on the LOH of chromosomes 1p, 19q and 13, but allelic loss of chromosome 14 in HF1551 was not detected by the Affymetrix array. The copy number profile generated by running the HMM algorithm on the maps from the first slice of tumor HF1551 was similar to that from the array, which suggests that the DNA originated from tumor sections that were in closer proximity.
Many of the SNP probes on the Affymetrix chip correspond to SwaI snip-SNPs. Hence, the array data was used to validate ECs and MCs. We observed 100% (62/62 in HF087, 44/44 in HF1551) concordance between the SNP genotype and the SwaI cut pattern at all overlapping cut differences in both tumors (Additional file4).
The copy number variants detected by the array were also compared to Optical Mapping indels. Signal intensities from the chip were normalized by global median scaling, and copy number was assessed using several different algorithms (methods), relative to a reference model file generated from the 270 HapMap samples. Though the resolution of array CGH is much lower than Optical Mapping, we were able to validate 24 structural variants in tumor HF087 and 16 in tumor HF1551 (Additional file5).
Experimental validation: PCR
The nature of many of the structural variants, being within repetitive portions of the genome, but detected by Optical Mapping unfortunately precludes their comprehensive validation by simple PCR techniques. Accordingly, we selected two variants that were amenable to PCR and overlapped genes that may offer insights into the chemo- and radio-sensitivity of oligodendroglioma. These loci were then PCR amplified, cloned and sequenced (methods).
The optical map shows an EC in the PARK2 gene in HF1551 (Figure 6A). PARK2 is a putative tumor suppressor, and mutations in this gene have been reported in multiple cancer types (detailed in ‘candidate mutations’ section of this document). An 848 bp amplicon spanning the predicted location of the EC was obtained (Figure 6B), and Sanger sequencing proved that a G to T transversion resulted in the creation of a new SwaI restriction site (Figure 6C).
We also validated an EC in tumor HF087 that occurred in the STMN2 gene (Figure 7A). As discussed in subsequent sections, STMN2 regulates microtubule dynamics and is believed to be a target of beta-catenin/TCF signalling. We amplified a 1003 bp region around the putative mutation (Figure 7B), and were able to validate the alteration via sequencing (Figure 7C).
We also validated our findings by comparing them to two sources- Optical Mapping data from several normal genomes, and publicly available SNP and structural variant data. First, oligodendroglioma structural variants were compared against structural variants found by Optical Mapping of 6 other normal human genomes by our laboratory. This internal database includes: (three lymphoblast-derived cell lines and a complete hydatiform mole (dbVar study ID nstd49), a lymphocyte-derived cell line (unpublished) and an early passage human embryonic stem cell line). 80%-90% of oligodendroglioma variants were also detected in at least one of the normal human genomes (Additional file6), suggesting that such loci are polymorphic, and affirming the veracity of our findings. Then, oligodendroglioma structural variants were compared against variants in the Database of Genomic Variants (DGV). The DGV is an extensive catalogue of structural variation in normal humans, currently holding 101,923 events detected by a variety of platforms. We observed the greatest concordance with variants found by fosmid-end sequencing (~15%) and high density oligonucleotide array CGH (~10%) (Additional file6). Finally, cut differences detected by Optical Mapping were compared to published SNPs. Detailed breakdown of these intersections are shown in Additional file3; parameters for comparisons are described in the Methods section.
Separation of mutational and polymorphic OSAs
The ultimate goal of our mapping efforts was to identify genes or genomic elements that maybe important to the biology of oligodendroglioma, with the caveat that such ‘candidates’ represent hypotheses requiring rigorous testing to establish their functional role in tumorigenesis. Distinguishing between structural polymorphisms and somatically acquired mutations is a key step towards accomplishing this goal. Unfortunately, matched normal DNA from the individuals whose tumors were optically mapped was not available. Instead, we adopted a stringent filtering scheme to remove putative polymorphisms and enrich for somatic mutations, based on comparisons to internal and publicly available data (described above). Parameters for these comparisons were determined based on the Optical Mapping error model and designed to be extremely parsimonious (methods). As a result of these operations, we arrived at a total of 21 somatic mutations (5 genes) in HF087 and 73 somatic mutations (21 genes) in HF1551.Since two mutations are seen in both tumors, 24 unique candidate cancer genes were identified in oligodendroglioma (Table 1). A few interesting candidate genes will be discussed in the next section.
Candidates common to both HF087 and HF1551
Two candidate genes, NPAS3 and OSBPL3, harbored mutations in both tumor samples (Figure 8). NPAS3 (neuronal PAS domain protein 3) shows a complex event accompanied by a ~7 kb gain in the HF087 optical map, and a missing cut in the HF1551 optical map (Figure 8A). This neuronally expressed basic helix-loop-helix transcription factor has been implicated in schizophrenia[76, 77]and bipolar disorder, and is frequently deleted or inactivated in many cancers. Recently, it has been demonstrated that NPAS3 exhibits features of a tumor suppressor which drives late progression of malignant astrocytomas, and is a negative prognostic marker for survival.
Both tumor optical maps display cut differences in the OSBPL3 (oxysterol binding protein like-3) gene (Figure 8B). This gene plays a vital role in cell adhesion, cytoskeletal organization and lipid metabolism[79–81]. It is highly expressed in B-cell associated malignancies[82, 83], where it is one of the common sites of retroviral integration. An independent study that used exon sequencing to study oligodendroglioma also found somatic mutations in OSBPL3.
Candidates observed in either HF087 or HF1551
In the HF1551 optical map, we observe a point mutation that creates a SwaI restriction site in the PARK2 gene (Figure 6). This gene encodes an E3 ubiquitin ligase, called Parkin that catalyzes the ubiquitination of a variety of target proteins for proteasome mediated degradation. Germline mutations in PARK2 have long been known to cause autosomal recessive juvenile Parkinson’s disease[85–87]. More recently, PARK2 has been identified as a tumor suppressor gene in Glioblastoma multiforme, breast, ovary, lung, colorectal and liver cancers[28, 88–94]. It encompasses most of FRA6E, the third most active common fragile site in the human genome, and shares the characteristics of other tumor suppressors such as FHIT and WWOX, that also occur in fragile sites. PARK2 is frequently deleted or inactivated in cancer cell lines and primary tumors[88, 92], and concomitantly, Parkin expression is either significantly diminished or absent[89, 92]. Unlike classical tumor suppressors where biallelic inactivation is necessary for oncogenesis, heterozygous mutations in PARK2 are sufficient to confer a growth advantage during tumor development[88, 92]. Restoring Parkin expression in Parkin-deficient cell lines reduces their profileration in vitro, while injection of Parkin-deficient cells into immunocompromised mice generate tumors in vivo. Interestingly, PARK2 also mediates chemosensitivity in breast cancer via microtubule dependent mechanism[93, 96–98].
STMN2 (stathmin-like 2) is another interesting candidate gene. We observe a point mutation in this gene in tumor HF087 (Figure 7). STMN2 is a neuron specific member of the stathmin family of small regulatory phosphoproteins which control cell profileration and differentiation. It is up-regulated in liver cancer and has been identified as a target of β-catenin/TCF-mediated transcription. STMN2 sequesters soluble tubulin, forming a ternary complex, inhibits microtubule assembly and induces their disassembly. Its highly similar, but more well-studied paralog STMN1, located on chromosome 1p, is known to sensitize cells to anti-microtubule drugs in glioma[102, 103], breast[104, 105] and prostate cancer.In light of recent studies demonstrating the synergistic epistasis between paralogous genes involved in essential cellular functions and its therapeutic implications[107, 108], we speculate that STMN1 and STMN2 might be functionally redundant, and inactivation of STMN2 might, in part, explain the treatment sensitivity of oligodendroglioma.
In the HF1551 optical map, we see an extra cut in the gene ZFYVE26 (zinc finger, FYVE domain containing 26). Spastizin, the zinc finger protein encoded by ZFYVE26, causes the neurological disorder hereditary spastic paraplegia. This gene binds to the tumor suppressor Beclin-1 and regulates cytokinesis[110, 111], and is recurrently mutated in breast cancer.
We detected a 485 kb inversion on 7q11.23 in tumor HF1551. Hemizygous deletions spanning a 1.5-1.8 MB region of this locus cause the neurodevelopmental disorder, Williams-Beuren (WB) syndrome. However, inversion of this region is polymorphic, and is present in ~6% of the general population, and in ~25% of transmitting parents in WB families[113–115]. Given the disparity in size between the aberration detected by OM and reported instances of the WB inversion, it is possible that the event we observe arose de novo and is distinct from the ‘canonical’ inversion. To test this hypothesis, we ran several targeted assemblies on the WB region. The general strategy for this approach was to modify the reference map in-silico to reflect our hypothesized structure (Figure 9A), and then use the iterative assembly framework described earlier to pull out individual restriction maps and generate an optical consensus map (methods). Since our map assembly pipeline was designed to provide the single most conservative answer, this approach is helpful in detecting large-scale aberrations that are significantly different from the reference sequence. Of the eight modified reference maps we started with, the one that reflected the canonical WB inversion/deletion did not grow, while the one that reflected the 485 kb event successfully generated a consensus map that spanned it and flanked contiguous regions on chromosome 7 (Figure 9B). The optical consensus map also closes the putative sequence gap immediately to the right of the inversion, and in fact, approximately half of the sequence gaps in the reference genome (NCBI, build 35) are spanned by optical consensus maps. The inversion encompasses the genes GTF2IRD2, PMS2P5, WBSCR16, GTF2IRD2B and NCF1, and its breakpoints appear to disrupt the genes GTF2I and STAG3L2. In the absence of matched normal DNA, it is impossible to ascertain if the inversion we detected was inherited through the germline or somatically acquired, however, this is the first report, to the best of our knowledge, of inversions in the WB region in the context of cancer.
Candidates on 1p or 19q
The concerted loss of chromosome arms 1p and 19q is a hallmark of oligodendroglioma. Seen in 50%-70% of tumors, it is believed that these regions harbor one or more tumor suppressor genes that play an important role in the development of this cancer. Hence, somatic mutations on these chromosome arms are particularly interesting. We found putative mutations on 2 genes residing on chromosome 1p (TCEB3, PAX7) and 1 gene on 19q (SIPA1L3). The roles of these genes in normal and disease states, and the structural variants we found in them are discussed briefly in the subsequent section.
We observe a 6.3 kb deletion that potentially ablates the first exon of TCEB3 in tumor HF1551. TCEB3 (transcription elongation factor B, polypeptide 3) encodes the transcriptionally active subunit of the mammalian elongin complex[116, 117]. This elongation factor stimulates the rate of transcription by suppressing the transient pausing of RNA polymerase II on the DNA template. TCEB3 is part of a multi-protein complex that functions as an elongin-based ubiquitin ligase, similar to the Von Hippel-Lindau (VHL) tumor suppressor complex, by mediating DNA damage induced ubiquitination and degradation of polymerase II.
Tumor HF1551 also has an insertion in the 1p-encoded gene PAX7 (paired box 7). The PAX genes encode a family of transcription factors that control development within the neural, myogenic and lymphoid lineages. PAX7, in particular, is essential for survival, proliferation and migration of myogenic progenitor cells, and cell fate decisions in the developing nervous system. PAX7 is the target of a recurrent gene fusion with the forkhead protein FKHR/FOXO1 that is found in ~15% of patients with alveolar rhabdomyosarcoma[121, 124]. The fusion transcript is much more abundant and transcriptionally active than wild type PAX7, suggesting that the deregulation of PAX7 downstream target genes contribute to tumorigenesis.
In the HF087 optical map, we observe a missing cut in the gene SIPA1L3 (signal induced proliferation associated 1 like 3) which is located on the long arm of chromosome 19. This gene encodes a Ras specific GTPase activating protein that is found at epithelial junctional complexes. These complexes play a crucial role in mechanical adhesion between epithelial cells to form cellular sheets and in the organization of actin cytoskeleton. Somatic mutations in SIPA1l3 have been discovered in cancers of the brain[127, 128], prostate, breast, ovary, pancreas, colon, skin and hematopoietic system, but a cohesive picture of the functional role that this gene plays in these diverse cancer types is yet to emerge.
Taken together, the candidate genes discovered by Optical Mapping point to critical roles of transcriptional control and cytoskeletal organisation in the etiology of oligodendroglioma.
Protein coding sequences comprise less than 2% of the human genome. The vast non-coding portion of the genome, once believed to be ‘junk DNA’, is rife with functional elements that orchestrate the gene expression program of cells. Recent evidence from the ENCODE (Encyclopedia of DNA Elements) consortium indicates that as much as 80.4% of the human genome encodes a defined product (for instance, a non-coding RNA) or displays a reproducible biochemical signature (for instance, a specific chromatin structure). Such signatures, either alone or in combinations, mark genomic sequences with important functions, such as promoters, enhancers, insulators and silencers. The ENCODE data sheds some light on possible functional roles of Optical Mapping candidates that are not located within genes. A number of these candidates are actively transcribed, for instance an EC on chromosome 5 of HF1551 overlaps the transcribed pseudogene GUSBP9. Several non-genic variants occur within long intergenic non-coding RNAs (lincRNA) coding regions (Additional file7). Both these classes of genomic elements provide an additional tier of gene regulation, and contribute significantly to the transcriptional landscape of human cancers[134, 135]. Several candidates also show interesting changes in their putative functions in cancer tissues. For example, we observe a MC on chromosome 2 (HF087) in a genomic region bearing a histone modification pattern characteristic of insulators in multiple different normal cell types, but the pattern changes to that of an enhancer in hepatocellular carcinoma.
Optical Mapping provides a global view of the cancer genome, free from biases introduced by cloning, amplification or hybridization, and discovers structural variation and mutation on a scale ranging from kilobases to megabases. Moreover, since the platform uses high-molecular weight DNA as analyte, the long-range context and connectivity of each variant is preserved, potentiating meaningful interpretation of candidate genes. However, Optical Mapping does not provide single-base resolution. Point mutations or indels spanning a few base pairs, such as the events frequently observed in CIC and FUBP1 genes in 1p/19q codeleted oligodendrogliomas, are below the lower limit of resolution and would remain undetected (unless they create or destroy a SwaI restriction site).
Biological significance of candidates identified by optical mapping
The aim of this study is to generate new hypotheses for oligodendroglioma genetics, and as such, functional studies are beyond the scope of this paper. However, by surveying publicly available data on the candidates discerned by Optical Mapping, we can gain some insight into the roles they might play in malignant transformation.
Moving beyond the two tumors HF087 and HF1551, we wanted to take the candidate cancer genes and analyze them in the context of other genome-wide studies. Most somatic mutations in cancer cells arise due to genomic instability and do not contribute to tumorigenesis. However, mutations in genes that promote tumor development, so-called ‘driver’ mutations tend to be recurrent. To assess the extent of recurrence of our candidates across a large number of samples, we used the Catalog of Somatic Mutations in Cancer (COSMIC). The COSMIC database is a comprehensive archive of somatic mutations in human cancer, combining manually curated data from scientific literature and the output from the Cancer Genome Project[136, 137]. 10 genes from our list of candidates had mutation frequencies (number of unique mutated samples divided by the total number of unique samples) greater than 10%, with the top hits being MYOF (19.1%), CECR2 (18.8%) and ZFYVE26 (18.7%). Mutation frequencies for all the candidate genes are listed in Table 2.
Since cells can employ a number of mechanisms to compensate for loss or mutational inactivation of genes, a more direct way of assessing the functional role of a given candidate gene is to analyze changes in its pattern of expression between normal and disease states. Array based expression profiling of tumors HF087 and HF1551 was performed by Fine et. al., and is publicly available through the NCBI GEO database, accession GSE4290. Differential expression analysis carried out using the EBarrays algorithm shows that 5 genes (NPAS3, STMN2, ZFYVE26, PHLDB2 and PLEKHM3) from our list of 24 candidates is significantly up or down-regulated (p-value 1E-03 or less). A complete list of differentially expressed genes can be found in Additional file8. To assay for functional effects in an even larger population of tumors, we queried for changes in expression of our candidate genes in REMBRANDT, a database of molecular data on brain tumors (National Cancer Institute, 2005, REMBRANDT home pagehttps://caintegrator.nci.nih.gov/rembrandt/, accessed 13thAugust 2012). The results are reported for each gene as the number of oligodendroglioma samples in the database that are differential expression by at least two-fold (Additional file9). All but 3 of our candidate genes were differentially expressed in at least 10 tumor samples. Congruent with the previous analysis, NPAS3, STMN2, ZFYVE26 and PHLDB2 are the most frequently deregulated candidate genes.
Finally, we asked what biological processes and pathways are significantly enriched or depleted in our list of candidates. This can identify fundamental cellular mechanisms that contribute to cancer development. As a whole, these candidate genes are enriched for proteins involved in cytoskeletal organization (p-value = 0.00223, after correcting for multiple testing, Ontologizer). Our candidate genes are also significantly enriched for microRNA binding targets (p-values between 0.0144-0.0198 after correcting for multiple testing, WebGestalt). Approximately half of the over-represented sites have been associated with binding of cancer-related microRNAs, underscoring the importance of post-transcriptional control of expression in oligodendroglioma.
While these results are not a direct indicator of aberrant function, this is a demonstration that Optical Mapping results can be expanded to clinical samples and used to create direct functional hypotheses.
We have applied Optical Mapping to explore the genomic landscape of solid tumor oligodendroglioma. ~2100 discrete structural variants have been discovered, ranging in size from single base changes to loss of entire chromosomes. The structure of each alteration has been elucidated at sub-genic resolution, while retaining the long-range context of the event. 94 somatic mutations have been identified, 24 of which affect genes. These novel candidate cancer genes provide focused, testable hypotheses for follow-up functional investigation. We believe that Optical Mapping provides a comprehensive, high-resolution description of the complex and disperse genomes of solid tumors.
Selection of tumors
The tumors used in this study originated from the tissue bank at the Hermelin Brain Tumor Center/Department of Neurosurgery, Henry Ford Hospital (provided by Dr. Oliver Bogler). Freshly resected tumors were snap frozen in liquid nitrogen in the operating room. Samples were sectioned in a guillotine in frozen condition, and adjacent pieces prepared for Optical Mapping and for re-review by a neuropathologist.
The tumor samples selected for Optical Mapping had to meet two criteria. First, they needed to conform to the 1p/19q paradigm of treatment sensitivity. LOH status was assessed by quantitative PCR of microsatellite markers along chromosomes 1p and 19q (data not shown). Second, they needed to have a high proportion of cancer cells as opposed to normal cells. The percentage of tumor cells present in each biopsy was estimated by MIB-1 antibody staining of an adjacent section (data not shown). The MIB-1 antibody recognizes the Ki-67 antigen, which is a cell proliferation marker. For the most part, mitotic activity is absent in the adult brain, so the measurement of the Ki-67 cell proliferation marker can be used to judge tumor aggression and composition[144, 145]. The two samples chosen for this project, HF087 and HF1551, satisfied both these criteria. Table 3 provides relevant clinical information for each tumor.
Extraction of high molecular weight DNA from solid tumor biopsies
The tumor was sectioned into 1–2 mm slices under sterile conditions in a cell culture hood. Each slice was treated with 0.8% type IV collagenase (Sigma-Aldrich, St. Louis, MO) in PBS (Phosphate buffered saline, Life Technologies, Carlsbad, CA) for 15 minutes at 37°C. The tumor tissue was mechanically disaggregated into a homogeneous suspension by repeated pipetting. The cells were pelleted by centrifugation at 1,000 RPM with a Beckman GS-6R centrifuge (Beckman Instruments, Fullerton, CA), and then resuspended in 1X HBS (Hanks Balanced Salts, Life Technologies, Carlsbad, CA) in order to lyse red blood cells. Cell debris and HBS were removed by centrifugation at 1,000 RPM. Finally, the pellet was rinsed three times with 35 mL of PBS, and resuspended in 0.5 mL of PBS.
A three layer Percoll gradient was employed to enrich for cancer cells, and minimize stromal contamination. First, a 100% solution was made by using 9 parts Percoll (Sigma-Aldrich, St. Louis, MO) and 1 part 10X HBS, which was subsequently diluted with PBS to prepare 10%, 30%, and 50% solutions. The gradient was prepared by layering 2 mL of 50% Percoll, 2 mL of 30% Percoll, and 1 mL of 10% Percoll in a 15 mL Falcon tube. The single cell suspension was then carefully layered on top, and the gradient was spun at 1,000 RPM for 10 minutes. Studies have shown that cellular debris and non-viable cells are unable to penetrate the 30% layer, while lymphocytes pelleted at the bottom of the tube. The 30% layer, containing viable cells, was carefully removed, rinsed three times with 10 mL of PBS and then resuspended in PBS at a final concentration of 1X107 cells/mL. Next, this cell suspension was mixed 1:1 (v/v) with 1.6% low gelling temperature agarose, poured into a mold and cooled to 4°C so that the agarose solidified to for gel inserts (each ~100 μL in volume).
The inserts were treated with 0.5 mg/mL proteinase K (Bioline USA, Taunton, MA), 100 mM EDTA pH 8.0 (Sigma-Aldrich, St. Louis, MO), 0.5% N-lauroylsarcosine (Sigma-Aldrich, St. Louis, MO) and incubated at 55°C overnight to lyse the tumor cells and degrade cellular proteins[45, 147–150]. Embedding cells in agarose inserts eliminates shear induced breakage of genomic DNA molecules upon lysis.
Prior to use, the gel inserts were rinsed in TE twice for 1 hour and then a third time overnight to remove the detergent and excess EDTA. DNA was electrophoretically extracted by applying a cycle of 100 V for 30 seconds and -100 V for 6 seconds.
Generation of single molecule optical maps
Optical Mapping surfaces were prepared as described earlier. Briefly, acid-cleaned glass coverslips (22 × 22 mm, Fisher's Finest, Fisher Scientific) were treated with a mixture of N-trimethoxylsilylpropyl-N,N,N-trimethylammonium chloride and vinyltrimethoxysilane (Gelest, Morrisville, PA) rendering a positive charge to the surface. Genomic DNA, mixed with a sizing standard, was elongated via c apillary flow in a microfluidic device, and immobilized by electrostatic interactions with the positively charged surface, creating arrays of stretched, biochemically accessible substrates. The surface was then washed with TE (10 mM Tris–HCl, 1 mM EDTA, pH8.0) twice, equilibrated with digestion buffer (NEB buffer 3), then incubated with the restriction endonuclease SwaI (New England Biolabs, Beverly, MA), which cleaves the genomic DNA at its cognate site. Since the elongated DNA molecule is under slight tension, upon cleavage its ends relax, creating a 1–2 micron gap, readily detected by microscopy. The resulting restriction fragments remain adsorbed to the surface, aided by a polyacrylamide overlay, and hence retain their order creating, in essence, a barcode from each genomic DNA molecule. Restriction fragments were then stained with the DNA intercalating dye YOYO-1 (0.2 μM in β-mercaptoethanol/TE, Life Technologies, Carlsbad, CA) and imaged by automated fluorescence microscopy.
The images were collected on an Optical Mapping workstation, which consists of Zeiss 135M inverted microscope (Carl Zeiss, Thornwood, NY), illuminated by 488 nm argon ion laser (Spectra Physics, Santa Clara, CA) equipped with 63X oil immersion objective. Fully automated image acquisition software, referred to as Channel Collect, takes multiple overlapping images to span the entire length of each microchannel. The images were analyzed by custom machine vision software, called Pathfinder, which identifies DNA molecules on the surface and calculates the size of each restriction fragment based on integrated fluorescence intensity measurements relative to a sizing standard. Previous studies have shown that integrated fluorescence intensity scales with fragment mass, and is independent of stretch of the DNA molecule. The end result of these operations is the high throughput, massively parallel generation of single molecule ordered restriction maps, or optical maps, containing information about both the size and order of its restriction fragments.
Pipeline for optical Map assembly and identification of structural variants
The analytical framework for assembly of optical maps is analogous to sequence assembly. First, our pipeline automatically aligned optical maps against a SwaI restriction map created in silico from the human genome reference sequence (NCBI build 35) via SOMA (Software for Optical Map Alignment) using gapped global pair wise alignment[56, 67]. SOMA uses a scoring function that assigns penalties for differences in the optical map and the reference map, including missing or extra restriction sites, or differences in the size of the fragments that could represent insertions or deletions. The parameters of SOMA were set so that we are accurately aligning the molecule to the correct location, but loose enough for allow for a small number of differences that result from the mutations or polymorphisms present in the genome. The aligned maps were then partitioned into smaller bins (1 Mb windows spanning across each chromosome, with 500 kb overlap between adjacent windows) based on their location. The optical maps in each bin were assembled into optical consensus maps by a map assembler program, using a Bayesian inference algorithm. Because some structural polymorphisms and mutations represent large-scale changes from the reference map, an iterative assembly process was used for the analysis of human data sets. The consensus map constructed in the previous step was used in place of the reference for seven more iterations of alignment and assembly, after which it was aligned to the reference sequence using SOMA. Using this strategy, Rmaps harboring major alterations that preclude alignment to the reference were gradually incorporated into the consensus map, extending it into regions that contain more complex rearrangements.
Lastly, the pipeline automatically performed analysis that tabulated structural variants using the final consensus map to the reference (derived from NCBI build 35 of the human genome) and identified five classes of differences: missing cuts, extra cuts, insertions, deletions, and ‘other’ (multiple cut and/or size differences) across each cancer genome. Each of these differences, which are largely structural variants, has to satisfy certain statistical and empirical criteria. These parameters have been detailed in Teague et al.. The only difference being the indel calling threshold, which was increased to a 13% change relative to the reference, with a 4.5 kb minimum.
Additionally, each structural variant was manually curated to ensure that the most conservative decision has been made at every locus. The genomic locations of the variants were converted to NCBI build 37 co-ordinates using the Batch Coordinate Conversion (liftover) tool from the University of California Santa Cruz Genome Browser (http://genome.ucsc.edu).
Optical map coverage analysis
Variations in depth of coverage of optical maps aligned by SOMA across the genome can be used to detect copy number alterations. Intuitively, if a region of the tumor sample has increased (or decreased) copy number relative to the ‘normal’ reference genome, more (or less) maps will originate from it on an average. This is formalized as described. Pair-wise alignments of optical maps to an in silico reference were summarized by a single number (midpoint) representing location. These locations were modeled as realizations of a non-homogeneous Poisson process. The non-homogeneity arises from the fact that the likelihood of a map aligning to a genomic region depends on the density of restriction sites, and was accounted for using alignment data from a normal genome, which are used to define random intervals with counts that follow a negative binomial distribution. These counts were then modeled by a Hidden Markov Model, incorporating spatial dependence in the data and allowing more natural estimation of certain parameters[68, 69].
Affymetrix genome wide human SNP array 6.0
DNA was prepared for hybridization using the Blood and Cell Culture Kit (Qiagen, Valencia, CA), starting from frozen cells (HF087), or tumor tissue (HF1551), disaggregated into single cells as described previously. The HF087 cells were derived from the same slice used for Optical Mapping. However, since the same was not available for HF1551, a slice adjacent to the one used for mapping was used.
The DNA was digested with NspI and StyI restriction enzymes and ligated to adaptors that recognize the 4 bp overhangs. A generic primer that anneals to the adaptor sequence was then used to amplify adaptor-ligated DNA fragments, under PCR conditions optimized to preferentially amplify fragments in the 200 to 1,100 bp size range. The amplified DNA was then fragmented, labelled, and hybridized to a Genome-Wide Human SNP 6.0 Array (experiments were performed by the DNA Facility at the Carver College of Medicine, University of Iowa). Data analysis was performed using Genotyping Console 2.0 (Affymetrix, Santa Clara, CA). CNVs were called using either the Affymetrix algorithm (with default parameters) or five different algorithms (GLAD, Circular Binary Segmentation, Fused Lasso, Gaussian Model with Adaptive Penalty, Forward-Backward Fragment Annealing Segmentation) from CGHweb (http://compbio.med.harvard.edu/CGHweb/). Only CNV calls made by two or more algorithms were considered for comparison.
Parameters for comparing oligodendroglioma structural variants
To other optical mapping datasets
Only variants of the same type were compared to each other, e.g.: MCs from HF087 were compared to MCs from lymphoblast cell line GM15510. Intersection ‘windows’ were set based on the type of OSA (100 bp for MCs, 4200 bp for ECs and 0 bp for INS, DEL and OTHER) and are reflective of the error processes inherent to each type of event.
To published SNPs and structural variants
Published SNPs were compared against Optical Mapping cut differences using 100 bp or 3000 bp windows for MCs and ECs, respectively.
Structural variants from the latest (November 2010) release of the Database of Genomic Variants were divided into two categories on the basis of their size. Events smaller than 3 kb were compared to ECs and MCs, since ~1/3rd of indels that are below the lower limit of detection for Optical Mapping manifest themselves as cut differences. Events larger than 3 kb were compared to INS, DEL and OTHER variants using a 0 bp intersection window.
Template for PCR was prepared by whole genome amplification of tumor DNA using the REPLI-g Mini kit (Qiagen Inc., Valencia, CA) as per the protocol provided by the manufacturer. Pooled normal DNA from 6 individuals (Promega Corporation, Madison, WI) was used for control reactions. Primers were designed using freely available software Primer 3 Plus. PCR reactions were performed using reagents from the Expand Long Template PCR System (Roche Applied Science, Indianapolis, IN) following the protocol supplied by the manufacturer. PCR reactions were digested with appropriate restriction enzymes to establish that the correct region had been amplified. The amplicon was then cloned in E. coli using the TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA), plasmid DNA was purified using the Qiagen Plasmid Mini Kit (Qiagen Inc., Valencia, CA), and sequenced using Sanger biochemistry.
Targeted assemblies on Williams-Beuren chromosomal region
The SwaI in-silico restriction map from the Williams-Beuren region on chromosome 14 was modified to reflect one of eight alterations: 4 possible inversions, each with unique start/end locations and spans (including the ‘canonical’ inversion), and 4 possible deletions, each with unique start/end locations and sizes (including the ‘canonical’ deletion). These modified in-silico maps were subjected to 8 rounds of iterative assembly, using the collection of HF1551 Rmaps, with the same parameters as the genome-wide assembly. The results were manually curated to rule out assembly errors.
This study was approved by the Institutional Review Board of the University of Wisconsin-Madison.
Availability of supporting data
All structural variation calls and analysis are contained within the additional files.
Single nucleotide polymorphism
Copy number variant
Hidden Markov Model
Loss of heterogeneity
Comparative genome hybridization
Major histocompatibility complex
Database of genomic variants
Catalog of somatic mutations in cancer
Encyclopedia of DNA elements.
Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL: Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987, 235 (4785): 177-182. 10.1126/science.3798106.
Bos JL: Ras oncogenes in human cancer: a review. Cancer Res. 1989, 49 (17): 4682-4689.
Murphree AL, Benedict WF: Retinoblastoma: clues to human oncogenesis. Science. 1984, 223 (4640): 1028-1033. 10.1126/science.6320372.
Orlow I, Lacombe L, Hannon GJ, Serrano M, Pellicer I, Dalbagni G, Reuter VE, Zhang ZF, Beach D, Cordon-Cardo C: Deletion of the p16 and p15 genes in human bladder tumors. J Natl Cancer Inst. 1995, 87 (20): 1524-1529. 10.1093/jnci/87.20.1524.
ar-Rushdi A, Nishikura K, Erikson J, Watt R, Rovera G, Croce CM: Differential expression of the translocated and the untranslocated c-myc oncogene in Burkitt lymphoma. Science. 1983, 222 (4622): 390-393. 10.1126/science.6414084.
Shtivelman E, Lifshitz B, Gale RP, Canaani E: Fused transcript of abl and bcr genes in chronic myelogenous leukaemia. Nature. 1985, 315 (6020): 550-554. 10.1038/315550a0.
Albertson DG, Collins C, McCormick F, Gray JW: Chromosome aberrations in solid tumors. Nat Genet. 2003, 34 (4): 369-376. 10.1038/ng1215.
Farabegoli F, Santini D, Ceccarelli C, Taffurelli M, Marrano D, Baldini N: Clone heterogeneity in diploid and aneuploid breast carcinomas as detected by FISH. Cytometry. 2001, 46 (1): 50-56. 10.1002/1097-0320(20010215)46:1<50::AID-CYTO1037>3.0.CO;2-T.
Marusyk A, Polyak K: Tumor heterogeneity: causes and consequences. Biochim Biophys Acta. 1805, 1: 105-117.
Navin N, Krasnitz A, Rodgers L, Cook K, Meth J, Kendall J, Riggs M, Eberling Y, Troge J, Grubor V: Inferring tumor progression from genomic heterogeneity. Genome Res. 2010, 20 (1): 68-80. 10.1101/gr.099622.109.
Zhang FF, Murata-Collins JL, Gaytan P, Forman SJ, Kopecky KJ, Willman CL, Appelbaum FR, Slovak ML: Twenty-four-color spectral karyotyping reveals chromosome aberrations in cytogenetically normal acute myeloid leukemia. Genes Chromosomes Cancer. 2000, 28 (3): 318-328. 10.1002/1098-2264(200007)28:3<318::AID-GCC10>3.0.CO;2-M.
Hicks DG, Tubbs RR: Assessment of the HER2 status in breast cancer by fluorescence in situ hybridization: a technical review with interpretive guidelines. Hum Pathol. 2005, 36 (3): 250-261. 10.1016/j.humpath.2004.11.010.
Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Kowbel D, Collins C, Kuo WL, Chen C, Zhai Y: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays. Nat Genet. 1998, 20 (2): 207-211. 10.1038/2524.
Ishkanian AS, Malloff CA, Watson SK, DeLeeuw RJ, Chi B, Coe BP, Snijders A, Albertson DG, Pinkel D, Marra MA: A tiling resolution DNA microarray with complete coverage of the human genome. Nat Genet. 2004, 36 (3): 299-303. 10.1038/ng1307.
Kallioniemi A, Kallioniemi OP, Sudar D, Rutovitz D, Gray JW, Waldman F, Pinkel D: Comparative genomic hybridization for molecular cytogenetic analysis of solid tumors. Science. 1992, 258 (5083): 818-821. 10.1126/science.1359641.
Pierotti MA, Santoro M, Jenkins RB, Sozzi G, Bongarzone I, Grieco M, Monzini N, Miozzo M, Herrmann MA, Fusco A: Characterization of an inversion on the long arm of chromosome 10 juxtaposing D10S170 and RET and creating the oncogenic sequence RET/PTC. Proc Natl Acad Sci USA. 1992, 89 (5): 1616-1620. 10.1073/pnas.89.5.1616.
Soda M, Choi YL, Enomoto M, Takada S, Yamashita Y, Ishikawa S, Fujiwara S, Watanabe H, Kurashina K, Hatanaka H: Identification of the transforming EML4-ALK fusion gene in non-small-cell lung cancer. Nature. 2007, 448 (7153): 561-566. 10.1038/nature05945.
Tomlins SA, Rhodes DR, Perner S, Dhanasekaran SM, Mehra R, Sun XW, Varambally S, Cao X, Tchinda J, Kuefer R: Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer. Science. 2005, 310 (5748): 644-648. 10.1126/science.1117679.
Hajirasouliha I, Hormozdiari F, Alkan C, Kidd JM, Birol I, Eichler EE, Sahinalp SC: Detection and characterization of novel sequence insertions using paired-end next-generation sequencing. Bioinformatics. 2010, 26 (10): 1277-1283. 10.1093/bioinformatics/btq152.
Liu G, Gramling S, Munoz D, Cheng D, Azad AK, Mirshams M, Chen Z, Xu W, Roberts H, Shepherd FA: Two novel BRM insertion promoter sequence variants are associated with loss of BRM expression and lung cancer risk. Oncogene. 2011, 30 (29): 3295-3304. 10.1038/onc.2011.81.
Bailey JA, Eichler EE: Primate segmental duplications: crucibles of evolution, diversity and disease. Nat Rev Genet. 2006, 7 (7): 552-564. 10.1038/nrg1895.
Mefford HC, Eichler EE: Duplication hotspots, rare genomic disorders, and common disease. Curr Opin Genet Dev. 2009, 19 (3): 196-204. 10.1016/j.gde.2009.04.003.
Alkan C, Coe BP, Eichler EE: Genome structural variation discovery and genotyping. Nat Rev Genet. 2011, 12 (5): 363-376. 10.1038/nrg2958.
Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, Shen D, Boca SM, Barber T, Ptak J: The genomic landscapes of human breast and colorectal cancers. Science. 2007, 318 (5853): 1108-1113. 10.1126/science.1145720.
Stephens P, Edkins S, Davies H, Greenman C, Cox C, Hunter C, Bignell G, Teague J, Smith R, Stevens C: A screen of the complete protein kinase gene family identifies diverse patterns of somatic mutations in human breast cancer. Nat Genet. 2005, 37 (6): 590-592. 10.1038/ng1571.
Bardelli A, Parsons DW, Silliman N, Ptak J, Szabo S, Saha S, Markowitz S, Willson JK, Parmigiani G, Kinzler KW: Mutational analysis of the tyrosine kinome in colorectal cancers. Science. 2003, 300 (5621): 949-10.1126/science.1082596.
Paez JG, Janne PA, Lee JC, Tracy S, Greulich H, Gabriel S, Herman P, Kaye FJ, Lindeman N, Boggon TJ: EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004, 304 (5676): 1497-1500. 10.1126/science.1099314.
Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486 (7403): 395-399.
Sjoblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N: The consensus coding sequences of human breast and colorectal cancers. Science. 2006, 314 (5797): 268-274. 10.1126/science.1133427.
Bettegowda C, Agrawal N, Jiao Y, Sausen M, Wood LD, Hruban RH, Rodriguez FJ, Cahill DP, McLendon R, Riggins G: Mutations in CIC and FUBP1 contribute to human oligodendroglioma. Science. 2011, 333 (6048): 1453-1455. 10.1126/science.1210557.
Agrawal N, Frederick MJ, Pickering CR, Bettegowda C, Chang K, Li RJ, Fakhry C, Xie TX, Zhang J, Wang J: Exome sequencing of head and neck squamous cell carcinoma reveals inactivating mutations in NOTCH1. Science. 2011, 333 (6046): 1154-1157. 10.1126/science.1206923.
Network CGAR: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.
Chapman MA, Lawrence MS, Keats JJ, Cibulskis K, Sougnez C, Schinzel AC, Harview CL, Brunet JP, Ahmann GJ, Adli M: Initial genome sequencing and analysis of multiple myeloma. Nature. 2011, 471 (7339): 467-472. 10.1038/nature09837.
Ding L, Ley TJ, Larson DE, Miller CA, Koboldt DC, Welch JS, Ritchey JK, Young MA, Lamprecht T, McLellan MD: Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing. Nature. 2012, 481 (7382): 506-510. 10.1038/nature10738.
Ding L, Ellis MJ, Li S, Larson DE, Chen K, Wallis JW, Harris CC, McLellan MD, Fulton RS, Fulton LL: Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010, 464 (7291): 999-1005. 10.1038/nature08989.
Lee W, Jiang Z, Liu J, Haverty PM, Guan Y, Stinson J, Yue P, Zhang Y, Pant KP, Bhatt D: The mutation spectrum revealed by paired genome sequences from a lung cancer patient. Nature. 2010, 465 (7297): 473-477. 10.1038/nature09004.
Ley TJ, Mardis ER, Ding L, Fulton B, McLellan MD, Chen K, Dooling D, Dunford-Shore BH, McGrath S, Hickenbotham M: DNA sequencing of a cytogenetically normal acute myeloid leukaemia genome. Nature. 2008, 456 (7218): 66-72. 10.1038/nature07485.
Pleasance ED, Cheetham RK, Stephens PJ, McBride DJ, Humphray SJ, Greenman CD, Varela I, Lin ML, Ordonez GR, Bignell GR: A comprehensive catalogue of somatic mutations from a human cancer genome. Nature. 2010, 463 (7278): 191-10.1038/nature08658.
Pleasance ED, Stephens PJ, O’Meara S, McBride DJ, Meynert A, Jones D, Lin ML, Beare D, Lau KW, Greenman C: A small-cell lung cancer genome with complex signatures of tobacco exposure. Nature. 2010, 463 (7278): 184-190. 10.1038/nature08629.
Shah SP, Morin RD, Khattra J, Prentice L, Pugh T, Burleigh A, Delaney A, Gelmon K, Guliany R, Senz J: Mutational evolution in a lobular breast tumour profiled at single nucleotide resolution. Nature. 2009, 461 (7265): 809-813. 10.1038/nature08489.
Mills RE, Walter K, Stewart C, Handsaker RE, Chen K, Alkan C, Abyzov A, Yoon SC, Ye K, Cheetham RK: Mapping copy number variation by population-scale genome sequencing. Nature. 2011, 470 (7332): 59-65. 10.1038/nature09708.
Alkan C, Sajjadian S, Eichler EE: Limitations of next-generation genome sequence assembly. Nat Methods. 2011, 8 (1): 61-65. 10.1038/nmeth.1527.
Consortium GP: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.
Medvedev P, Stanciu M, Brudno M: Computational methods for discovering structural variation with next-generation sequencing. Nat Methods. 2009, 6 (11 Suppl): S13-S20.
Reslewic S, Zhou S, Place M, Zhang Y, Briska A, Goldstein S, Churas C, Runnheim R, Forrest D, Lim A: Whole-genome shotgun Optical Mapping of Rhodospirillum rubrum. Appl Environ Microbiol. 2005, 71 (9): 5511-5522. 10.1128/AEM.71.9.5511-5522.2005.
Lai Z, Jing J, Aston C, Clarke V, Apodaca J, Dimalanta ET, Carucci DJ, Gardner MJ, Mishra B, Anantharaman TS: A shotgun optical map of the entire Plasmodium falciparum genome. Nat Genet. 1999, 23 (3): 309-313. 10.1038/15484.
Lin J, Qi R, Aston C, Jing J, Anantharaman TS, Mishra B, White O, Daly MJ, Minton KW, Venter JC: Whole-genome shotgun Optical Mapping of Deinococcus radiodurans. Science. 1999, 285 (5433): 1558-1562. 10.1126/science.285.5433.1558.
Schwartz DC, Li X, Hernandez LI, Ramnarain SP, Huff EJ, Wang YK: Ordered restriction maps of Saccharomyces cerevisiae chromosomes constructed by Optical Mapping. Science. 1993, 262 (5130): 110-114. 10.1126/science.8211116.
Young ND, Debelle F, Oldroyd GE, Geurts R, Cannon SB, Udvardi MK, Benedito VA, Mayer KF, Gouzy J, Schoof H: The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011, 480 (7378): 520-524.
Zhou S, Wei F, Nguyen J, Bechner M, Potamousis K, Goldstein S, Pape L, Mehan MR, Churas C, Pasternak S: A single molecule scaffold for the maize genome. PLoS Genet. 2009, 5 (11): e1000711-10.1371/journal.pgen.1000711.
Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA: The B73 maize genome: complexity, diversity, and dynamics. Science. 2009, 326 (5956): 1112-1115. 10.1126/science.1178534.
Zhou S, Bechner MC, Place M, Churas CP, Pape L, Leong SA, Runnheim R, Forrest DK, Goldstein S, Livny M: Validation of rice genome sequence by Optical Mapping. BMC Genomics. 2007, 8: 278-10.1186/1471-2164-8-278.
Kidd JM, Cooper GM, Donahue WF, Hayden HS, Sampas N, Graves T, Hansen N, Teague B, Alkan C, Antonacci F: Mapping and sequencing of structural variation from eight human genomes. Nature. 2008, 453 (7191): 56-64. 10.1038/nature06862.
Zody MC, Garber M, Adams DJ, Sharpe T, Harrow J, Lupski JR, Nicholson C, Searle SM, Wilming L, Young SK: DNA sequence of human chromosome 17 and analysis of rearrangement in the human lineage. Nature. 2006, 440 (7087): 1045-1049. 10.1038/nature04689.
Ananiev GE, Goldstein S, Runnheim R, Forrest DK, Zhou S, Potamousis K, Churas CP, Bergendahl V, Thomson JA, Schwartz DC: Optical Mapping discerns genome wide DNA methylation profiles. BMC Mol Biol. 2008, 9: 68-10.1186/1471-2199-9-68.
Teague B, Waterman MS, Goldstein S, Potamousis K, Zhou S, Reslewic S, Sarkar D, Valouev A, Churas C, Kidd JM: High-resolution human genome structure by single-molecule analysis. Proc Natl Acad Sci USA. 2010, 107 (24): 10848-10853. 10.1073/pnas.0914638107.
Church DM, Goodstadt L, Hillier LW, Zody MC, Goldstein S, She X, Bult CJ, Agarwala R, Cherry JL, DiCuccio M: Lineage-specific biology revealed by a finished genome assembly of the mouse. PLoS Biol. 2009, 7 (5): e1000112-10.1371/journal.pbio.1000112.
DeAngelis LM: Brain tumors. N Engl J Med. 2001, 344 (2): 114-123. 10.1056/NEJM200101113440207.
Reifenberger G, Louis DN: Oligodendroglioma: toward molecular definitions in diagnostic neuro-oncology. J Neuropathol Exp Neurol. 2003, 62 (2): 111-126.
Cairncross G, Jenkins R: Gliomas with 1p/19q codeletion: a.k.a. oligodendroglioma. Cancer J. 2008, 14 (6): 352-357. 10.1097/PPO.0b013e31818d8178.
Cairncross JG, Ueki K, Zlatescu MC, Lisle DK, Finkelstein DM, Hammond RR, Silver JS, Stark PC, Macdonald DR, Ino Y: Specific genetic predictors of chemotherapeutic response and survival in patients with anaplastic oligodendrogliomas. J Natl Cancer Inst. 1998, 90 (19): 1473-1479. 10.1093/jnci/90.19.1473.
Walker C, du Plessis DG, Joyce KA, Fildes D, Gee A, Haylock B, Husband D, Smith T, Broome J, Warnke PC: Molecular pathology and clinical characteristics of oligodendroglial neoplasms. Ann Neurol. 2005, 57 (6): 855-865. 10.1002/ana.20496.
Bromberg JE, van den Bent MJ: Oligodendrogliomas: molecular biology and treatment. Oncologist. 2009, 14 (2): 155-163. 10.1634/theoncologist.2008-0248.
Wharton SB, Maltby E, Jellinek DA, Levy D, Atkey N, Hibberd S, Crimmins D, Stoeber K, Williams GH: Subtypes of oligodendroglioma defined by 1p,19q deletions, differ in the proportion of apoptotic cells but not in replication-licensed non-proliferating cells. Acta Neuropathol. 2007, 113 (2): 119-127. 10.1007/s00401-006-0177-2.
van den Bent MJ, Looijenga LH, Langenberg K, Dinjens W, Graveland W, Uytdewilligen L, Sillevis Smitt PA, Jenkins RB, Kros JM: Chromosomal anomalies in oligodendroglial tumors are correlated with clinical features. Cancer. 2003, 97 (5): 1276-1284. 10.1002/cncr.11187.
Dimalanta ET, Lim A, Runnheim R, Lamers C, Churas C, Forrest DK, de Pablo JJ, Graham MD, Coppersmith SN, Goldstein S: A microfluidic system for large DNA molecule arrays. Anal Chem. 2004, 76 (18): 5293-5301. 10.1021/ac0496401.
Valouev A, Li L, Liu YC, Schwartz DC, Yang Y, Zhang Y, Waterman MS: Alignment of optical maps. J Comput Biol. 2006, 13 (2): 442-462. 10.1089/cmb.2006.13.442.
Sarkar D: On the analysis of Optical Mapping data. 2006, Madison, WI: University of Wisconsin
Sarkar D, Goldstein S, Schwartz DC, Newton MA: Statistical significance of optical map alignments. J Comput Biol. 2012, 19 (5): 478-492. 10.1089/cmb.2011.0221.
Griffin CA, Burger P, Morsberger L, Yonescu R, Swierczynski S, Weingart JD, Murphy KM: Identification of der(1;19)(q10;p10) in five oligodendrogliomas suggests mechanism of concurrent 1p and 19q loss. J Neuropathol Exp Neurol. 2006, 65 (10): 988-994. 10.1097/01.jnen.0000235122.98052.8f.
Jenkins RB, Blair H, Ballman KV, Giannini C, Arusell RM, Law M, Flynn H, Passe S, Felten S, Brown PD: A t(1;19)(q10;p10) mediates the combined deletions of 1p and 19q and predicts a better prognosis of patients with oligodendroglioma. Cancer Res. 2006, 66 (20): 9852-9861. 10.1158/0008-5472.CAN-06-1796.
Anderson K, Lutz C, van Delft FW, Bateman CM, Guo Y, Colman SM, Kempski H, Moorman AV, Titley I, Swansbury J: Genetic variegation of clonal architecture and propagating cells in leukaemia. Nature. 2011, 469 (7330): 356-361. 10.1038/nature09650.
Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, Martinez P, Matthews N, Stewart A, Tarpey P: Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012, 366 (10): 883-892. 10.1056/NEJMoa1113205.
Yap TA, Gerlinger M, Futreal PA, Pusztai L, Swanton C: Intratumor heterogeneity: seeing the wood for the trees. Sci Transl Med. 2012, 4 (127): 127-110.
Wicks SR, Yeh RT, Gish WR, Waterston RH, Plasterk RH: Rapid gene mapping in Caenorhabditis elegans using a high density polymorphism map. Nat Genet. 2001, 28 (2): 160-164. 10.1038/88878.
Macintyre G, Alford T, Xiong L, Rouleau GA, Tibbo PG, Cox DW: Association of NPAS3 exonic variation with schizophrenia. Schizophr Res. 2010, 120 (1–3): 143-149.
Pickard BS, Christoforou A, Thomson PA, Fawkes A, Evans KL, Morris SW, Porteous DJ, Blackwood DH, Muir WJ: Interacting haplotypes at the NPAS3 locus alter risk of schizophrenia and bipolar disorder. Mol Psychiatry. 2009, 14 (9): 874-884. 10.1038/mp.2008.24.
Moreira F, Kiehl TR, So K, Ajeawung NF, Honculada C, Gould P, Pieper RO, Kamnasaran D: NPAS3 demonstrates features of a tumor suppressive role in driving the progression of Astrocytomas. Am J Pathol. 2011, 179 (1): 462-476. 10.1016/j.ajpath.2011.03.044.
Lehto M, Olkkonen VM: The OSBP-related proteins: a novel protein family involved in vesicle transport, cellular lipid metabolism, and cell signalling. Biochim Biophys Acta. 2003, 1631 (1): 1-11. 10.1016/S1388-1981(02)00364-5.
Lehto M, Tienari J, Lehtonen S, Lehtonen E, Olkkonen VM: Subfamily III of mammalian oxysterol-binding protein (OSBP) homologues: the expression and intracellular localization of ORP3, ORP6, and ORP7. Cell Tissue Res. 2004, 315 (1): 39-57. 10.1007/s00441-003-0817-y.
Lehto M, Mayranpaa MI, Pellinen T, Ihalmo P, Lehtonen S, Kovanen PT, Groop PH, Ivaska J, Olkkonen VM: The R-Ras interaction partner ORP3 regulates cell adhesion. J Cell Sci. 2008, 121 (Pt 5): 695-705.
Sander B, Flygare J, Porwit-Macdonald A, Smith CI, Emanuelsson E, Kimby E, Liden J, Christensson B: Mantle cell lymphomas with low levels of cyclin D1 long mRNA transcripts are highly proliferative and can be discriminated by elevated cyclin A2 and cyclin B1. Int J Cancer. 2005, 117 (3): 418-430. 10.1002/ijc.21166.
Ek S, Andreasson U, Hober S, Kampf C, Ponten F, Uhlen M, Merz H, Borrebaeck CA: From gene expression analysis to tissue microarrays: a rational approach to identify therapeutic and diagnostic targets in lymphoid malignancies. Mol Cell Proteomics. 2006, 5 (6): 1072-1081. 10.1074/mcp.M600077-MCP200.
Bijl J, Sauvageau M, Thompson A, Sauvageau G: High incidence of proviral integrations in the Hoxa locus in a new model of E2a-PBX1-induced B-cell leukemia. Genes Dev. 2005, 19 (2): 224-233. 10.1101/gad.1268505.
Abbas N, Lucking CB, Ricard S, Durr A, Bonifati V, De Michele G, Bouley S, Vaughan JR, Gasser T, Marconi R: A wide variety of mutations in the parkin gene are responsible for autosomal recessive parkinsonism in Europe. French Parkinson’s Disease Genetics Study Group and the European Consortium on Genetic Susceptibility in Parkinson's Disease. Hum Mol Genet. 1999, 8 (4): 567-574. 10.1093/hmg/8.4.567.
Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, Yokochi M, Mizuno Y, Shimizu N: Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. 1998, 392 (6676): 605-608. 10.1038/33416.
Lucking CB, Durr A, Bonifati V, Vaughan J, De Michele G, Gasser T, Harhangi BS, Meco G, Denefle P, Wood NW: Association between early-onset Parkinson’s disease and mutations in the parkin gene. N Engl J Med. 2000, 342 (21): 1560-1567. 10.1056/NEJM200005253422103.
Veeriah S, Taylor BS, Meng S, Fang F, Yilmaz E, Vivanco I, Janakiraman M, Schultz N, Hanrahan AJ, Pao W: Somatic mutations of the Parkinson’s disease-associated gene PARK2 in glioblastoma and other human malignancies. Nat Genet. 2010, 42 (1): 77-82. 10.1038/ng.491.
Cesari R, Martin ES, Calin GA, Pentimalli F, Bichi R, McAdams H, Trapasso F, Drusco A, Shimizu M, Masciullo V: Parkin, a gene implicated in autosomal recessive juvenile parkinsonism, is a candidate tumor suppressor gene on chromosome 6q25-q27. Proc Natl Acad Sci USA. 2003, 100 (10): 5956-5961. 10.1073/pnas.0931262100.
Denison SR, Wang F, Becker NA, Schule B, Kock N, Phillips LA, Klein C, Smith DI: Alterations in the common fragile site gene Parkin in ovarian and other cancers. Oncogene. 2003, 22 (51): 8370-8378. 10.1038/sj.onc.1207072.
Fujiwara M, Marusawa H, Wang HQ, Iwai A, Ikeuchi K, Imai Y, Kataoka A, Nukina N, Takahashi R, Chiba T: Parkin as a tumor suppressor gene for hepatocellular carcinoma. Oncogene. 2008, 27 (46): 6002-6011. 10.1038/onc.2008.199.
Poulogiannis G, McIntyre RE, Dimitriadi M, Apps JR, Wilson CH, Ichimura K, Luo F, Cantley LC, Wyllie AH, Adams DJ: PARK2 deletions occur frequently in sporadic colorectal cancer and accelerate adenoma development in Apc mutant mice. Proc Natl Acad Sci USA. 2010, 107 (34): 15145-15150. 10.1073/pnas.1009941107.
Wang F, Denison S, Lai JP, Philips LA, Montoya D, Kock N, Schule B, Klein C, Shridhar V, Roberts LR: Parkin gene alterations in hepatocellular carcinoma. Genes Chromosomes Cancer. 2004, 40 (2): 85-96. 10.1002/gcc.20020.
Picchio MC, Martin ES, Cesari R, Calin GA, Yendamuri S, Kuroki T, Pentimalli F, Sarti M, Yoder K, Kaiser LR: Alterations of the tumor suppressor gene Parkin in non-small cell lung cancer. Clin Cancer Res. 2004, 10 (8): 2720-2724. 10.1158/1078-0432.CCR-03-0086.
Kahkonen M: Population cytogenetics of folate-sensitive fragile sites. I. Common fragile sites. Hum Genet. 1988, 80 (4): 344-348. 10.1007/BF00273649.
Ren Y, Jiang H, Yang F, Nakaso K, Feng J: Parkin protects dopaminergic neurons against microtubule-depolymerizing toxins by attenuating microtubule-associated protein kinase activation. J Biol Chem. 2009, 284 (6): 4009-4017.
Wang H, Liu B, Zhang C, Peng G, Liu M, Li D, Gu F, Chen Q, Dong JT, Fu L: Parkin regulates paclitaxel sensitivity in breast cancer via a microtubule-dependent mechanism. J Pathol. 2009, 218 (1): 76-85. 10.1002/path.2512.
Yang F, Jiang Q, Zhao J, Ren Y, Sutton MD, Feng J: Parkin stabilizes microtubules through strong binding mediated by three independent domains. J Biol Chem. 2005, 280 (17): 17154-17162. 10.1074/jbc.M500843200.
Manna T, Grenningloh G, Miller HP, Wilson L: Stathmin family protein SCG10 differentially regulates the plus and minus end dynamics of microtubules at steady state in vitro: implications for its role in neurite outgrowth. Biochemistry. 2007, 46 (11): 3543-3552. 10.1021/bi061819d.
Lee HS, Lee DC, Park MH, Yang SJ, Lee JJ, Kim DM, Jang Y, Lee JH, Choi JY, Kang YK: STMN2 is a novel target of beta-catenin/TCF-mediated transcription in human hepatoma cells. Biochem Biophys Res Commun. 2006, 345 (3): 1059-1067. 10.1016/j.bbrc.2006.05.017.
Riederer BM, Pellier V, Antonsson B, Di Paolo G, Stimpson SA, Lutjens R, Catsicas S, Grenningloh G: Regulation of microtubule dynamics by the neuronal growth-associated protein SCG10. Proc Natl Acad Sci USA. 1997, 94 (2): 741-745. 10.1073/pnas.94.2.741.
Ngo TT, Peng T, Liang XJ, Akeju O, Pastorino S, Zhang W, Kotliarov Y, Zenklusen JC, Fine HA, Maric D: The 1p-encoded protein stathmin and resistance of malignant gliomas to nitrosoureas. J Natl Cancer Inst. 2007, 99 (8): 639-652. 10.1093/jnci/djk135.
Liang XJ, Choi Y, Sackett DL, Park JK: Nitrosoureas inhibit the stathmin-mediated migration and invasion of malignant glioma cells. Cancer Res. 2008, 68 (13): 5267-5272. 10.1158/0008-5472.CAN-07-6482.
Alli E, Bash-Babula J, Yang JM, Hait WN: Effect of stathmin on the sensitivity to antimicrotubule drugs in human breast cancer. Cancer Res. 2002, 62 (23): 6864-6869.
Alli E, Yang JM, Hait WN: Silencing of stathmin induces tumor-suppressor function in breast cancer cell lines harboring mutant p53. Oncogene. 2007, 26 (7): 1003-1012. 10.1038/sj.onc.1209864.
Mistry SJ, Atweh GF: Therapeutic interactions between stathmin inhibition and chemotherapeutic agents in prostate cancer. Mol Cancer Ther. 2006, 5 (12): 3248-3257. 10.1158/1535-7163.MCT-06-0227.
Muller FL, Colla S, Aquilanti E, Manzo VE, Genovese G, Lee J, Eisenson D, Narurkar R, Deng P, Nezi L: Passenger deletions generate therapeutic vulnerabilities in cancer. Nature. 2012, 488 (7411): 337-342. 10.1038/nature11331.
Nijhawan D, Zack TI, Ren Y, Strickland MR, Lamothe R, Schumacher SE, Tsherniak A, Besche HC, Rosenbluh J, Shehata S: Cancer vulnerabilities unveiled by genomic loss. Cell. 2012, 150 (4): 842-854. 10.1016/j.cell.2012.07.023.
Hanein S, Martin E, Boukhris A, Byrne P, Goizet C, Hamri A, Benomar A, Lossos A, Denora P, Fernandez J: Identification of the SPG15 gene, encoding spastizin, as a frequent cause of complicated autosomal-recessive spastic paraplegia, including Kjellin syndrome. Am J Hum Genet. 2008, 82 (4): 992-1002. 10.1016/j.ajhg.2008.03.004.
Sagona AP, Nezis IP, Bache KG, Haglund K, Bakken AC, Skotheim RI, Stenmark H: A tumor-associated mutation of FYVE-CENT prevents its interaction with Beclin 1 and interferes with cytokinesis. PLoS One. 2011, 6 (3): e17086-10.1371/journal.pone.0017086.
Sagona AP, Nezis IP, Pedersen NM, Liestol K, Poulton J, Rusten TE, Skotheim RI, Raiborg C, Stenmark H: PtdIns(3)P controls cytokinesis through KIF13A-mediated recruitment of FYVE-CENT to the midbody. Nat Cell Biol. 2010, 12 (4): 362-371. 10.1038/ncb2036.
Pober BR: Williams-Beuren syndrome. N Engl J Med. 2010, 362 (3): 239-252. 10.1056/NEJMra0903074.
Hobart HH, Morris CA, Mervis CB, Pani AM, Kistler DJ, Rios CM, Kimberley KW, Gregg RG, Bray-Ward P: Inversion of the Williams syndrome region is a common polymorphism found more frequently in parents of children with Williams syndrome. Am J Med Genet C Semin Med Genet. 2010, 154C (2): 220-228. 10.1002/ajmg.c.30258.
Osborne LR, Li M, Pober B, Chitayat D, Bodurtha J, Mandel A, Costa T, Grebe T, Cox S, Tsui LC: A 1.5 million-base pair inversion polymorphism in families with Williams-Beuren syndrome. Nat Genet. 2001, 29 (3): 321-325. 10.1038/ng753.
Tam E, Young EJ, Morris CA, Marshall CR, Loo W, Scherer SW, Mervis CB, Osborne LR: The common inversion of the Williams-Beuren syndrome region at 7q11.23 does not cause clinical symptoms. Am J Med Genet A. 2008, 146A (14): 1797-1806. 10.1002/ajmg.a.32360.
Aso T, Lane WS, Conaway JW, Conaway RC: Elongin (SIII): a multisubunit regulator of elongation by RNA polymerase II. Science. 1995, 269 (5229): 1439-1443. 10.1126/science.7660129.
Conaway JW, Bradsher JN, Tan S, Conaway RC: Transcription factor SIII: a novel component of the RNA polymerase II elongation complex. Cell Mol Biol Res. 1993, 39 (4): 323-329.
Bradsher JN, Jackson KW, Conaway RC, Conaway JW: RNA polymerase II transcription factor SIII. I. Identification, purification, and properties. J Biol Chem. 1993, 268 (34): 25587-25593.
Kamura T, Burian D, Yan Q, Schmidt SL, Lane WS, Querido E, Branton PE, Shilatifard A, Conaway RC, Conaway JW: Muf1, a novel Elongin BC-interacting leucine-rich repeat protein that can assemble with Cul5 and Rbx1 to reconstitute a ubiquitin ligase. J Biol Chem. 2001, 276 (32): 29748-29753. 10.1074/jbc.M103093200.
Yasukawa T, Kamura T, Kitajima S, Conaway RC, Conaway JW, Aso T: Mammalian Elongin A complex mediates DNA-damage-induced ubiquitylation and degradation of Rpb1. EMBO J. 2008, 27 (24): 3256-3266. 10.1038/emboj.2008.249.
Barr FG: Chromosomal translocations involving paired box transcription factors in human cancer. Int J Biochem Cell Biol. 1997, 29 (12): 1449-1461. 10.1016/S1357-2725(97)00095-2.
Buckingham M, Relaix F: The role of Pax genes in the development of tissues and organs: Pax3 and Pax7 regulate muscle progenitor cell functions. Annu Rev Cell Dev Biol. 2007, 23: 645-673. 10.1146/annurev.cellbio.23.090506.123438.
Murdoch B, Delconte C, Garcia-Castro MI: Pax7 lineage contributions to the Mammalian neural crest. PLoS One. 2012, 7 (7): e41089-10.1371/journal.pone.0041089.
Davis RJ, D’Cruz CM, Lovell MA, Biegel JA, Barr FG: Fusion of PAX7 to FKHR by the variant t(1;13)(p36;q14) translocation in alveolar rhabdomyosarcoma. Cancer Res. 1994, 54 (11): 2869-2872.
Bennicelli JL, Advani S, Schafer BW, Barr FG: PAX3 and PAX7 exhibit conserved cis-acting transcription repression domains and utilize a common gain of function mechanism in alveolar rhabdomyosarcoma. Oncogene. 1999, 18 (30): 4348-4356. 10.1038/sj.onc.1202812.
Matsuda M, Kobayashi Y, Masuda S, Adachi M, Watanabe T, Yamashita JK, Nishi E, Tsukita S, Furuse M: Identification of adherens junction-associated GTPase activating proteins by the fluorescence localization-based expression cloning. Exp Cell Res. 2008, 314 (5): 939-949. 10.1016/j.yexcr.2007.11.009.
Parsons DW, Jones S, Zhang X, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Siu IM, Gallia GL: An integrated genomic analysis of human glioblastoma multiforme. Science. 2008, 321 (5897): 1807-1812. 10.1126/science.1164382.
Parsons DW, Li M, Zhang X, Jones S, Leary RJ, Lin JC, Boca SM, Carter H, Samayoa J, Bettegowda C: The genetic landscape of the childhood cancer medulloblastoma. Science. 2011, 331 (6016): 435-439. 10.1126/science.1198056.
Grasso CS, Wu YM, Robinson DR, Cao X, Dhanasekaran SM, Khan AP, Quist MJ, Jing X, Lonigro RJ, Brenner JC: The mutational landscape of lethal castration-resistant prostate cancer. Nature. 2012, 487 (7406): 239-243. 10.1038/nature11125.
Jones S, Zhang X, Parsons DW, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Kamiyama H, Jimeno A: Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science. 2008, 321 (5897): 1801-1806. 10.1126/science.1164368.
Berger MF, Hodis E, Heffernan TP, Deribe YL, Lawrence MS, Protopopov A, Ivanova E, Watson IR, Nickerson E, Ghosh P: Melanoma genome sequencing reveals frequent PREX2 mutations. Nature. 2012, 485 (7399): 502-506.
Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, Snyder M: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489 (7414): 57-74. 10.1038/nature11247.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473 (7345): 43-49. 10.1038/nature09906.
Kalyana-Sundaram S, Kumar-Sinha C, Shankar S, Robinson DR, Wu YM, Cao X, Asangani IA, Kothari V, Prensner JR, Lonigro RJ: Expressed pseudogenes in the transcriptional landscape of human cancers. Cell. 2012, 149 (7): 1622-1634. 10.1016/j.cell.2012.04.041.
Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E: Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007, 129 (7): 1311-1323. 10.1016/j.cell.2007.05.022.
Forbes SA, Bhamra G, Bamford S, Dawson E, Kok C, Clements J, Menzies A, Teague JW, Futreal PA, Stratton MR: The Catalogue of Somatic Mutations in Cancer (COSMIC). Curr Protoc Hum Genet. 2008, 10: 10-11.
Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2011, 39: 945-950. 10.1093/nar/gkq929.
Sun L, Hui AM, Su Q, Vortmeyer A, Kotliarov Y, Pastorino S, Passaniti A, Menon J, Walling J, Bailey R: Neuronal and glioma-derived stem cell factor induces angiogenesis within the brain. Cancer Cell. 2006, 9 (4): 287-300. 10.1016/j.ccr.2006.03.003.
Edgar R, Domrachev M, Lash AE: Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.
Yuan M, Kendziorski C: A unified approach for simultaneous gene clustering and differential expression identification. Biometrics. 2006, 62 (4): 1089-1098. 10.1111/j.1541-0420.2006.00611.x.
Bauer S, Grossmann S, Vingron M, Robinson PN: Ontologizer 2.0–a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics. 2008, 24 (14): 1650-1651. 10.1093/bioinformatics/btn250.
Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33: W741-W748. 10.1093/nar/gki475.
Wuchty S, Arjona D, Bozdag S, Bauer PO: Involvement of microRNA families in cancer. Nucleic Acids Res. 2012, 40 (17): 8219-8226. 10.1093/nar/gks627.
Coons SW, Johnson PC, Scheithauer BW, Yates AJ, Pearl DK: Improving diagnostic accuracy and interobserver concordance in the classification and grading of primary gliomas. Cancer. 1997, 79 (7): 1381-1393. 10.1002/(SICI)1097-0142(19970401)79:7<1381::AID-CNCR16>3.0.CO;2-W.
Reis-Filho JS, Faoro LN, Carrilho C, Bleggi-Torres LF, Schmitt FC: Evaluation of cell proliferation, epidermal growth factor receptor, and bcl-2 immunoexpression as prognostic factors for patients with World Health Organization grade 2 oligodendroglioma. Cancer. 2000, 88 (4): 862-869. 10.1002/(SICI)1097-0142(20000215)88:4<862::AID-CNCR17>3.0.CO;2-G.
VanDevanter DR, Yirdaw G, Do C, Tysseling KA, Drescher CA, Forseth BJ, Von Hoff DD, McNutt MA: Preparation of uniform, intact DNA samples from resected tumor tissues for pulsed-field gel electrophoretic analyses. Biotechniques. 1992, 13 (6): 884-887.
Lim A, Dimalanta ET, Potamousis KD, Yen G, Apodoca J, Tao C, Lin J, Qi R, Skiadas J, Ramanathan A: Shotgun optical maps of the whole Escherichia coli O157:H7 genome. Genome Res. 2001, 11 (9): 1584-1593. 10.1101/gr.172101.
Lim A: Single molecule systems: Advancements and applications to microbial and human genome analysis. 2004, PhD Thesis. University of Wisconsin-Madison
Reslewic S: The Optical Mapping of genomes: Gaining new insights on genomic structure and variation by single DNA molecule analysis. 2005, PhD Thesis. University of Wisconsin-Madison
Zhou S, Herschleb J, Schwartz DC: Optical Mapping: a single molecule system for genome analysis. New methods for DNA sequencing. 2007, 265-300.
Schwartz DC, Cantor CR: Separation of yeast chromosome-sized DNAs by pulsed field gradient gel electrophoresis. Cell. 1984, 37 (1): 67-75. 10.1016/0092-8674(84)90301-5.
Zhou S, Deng W, Anantharaman TS, Lim A, Dimalanta ET, Wang J, Wu T, Chunhong T, Creighton R, Kile A: A whole-genome shotgun optical map of Yersinia pestis strain KIM. Appl Environ Microbiol. 2002, 68 (12): 6321-6331. 10.1128/AEM.68.12.6321-6331.2002.
Anantharaman T, Mishra B, Schwartz D: Genomics via Optical Mapping: III: Contiging genomic DNA. Proc Int Conf Intell Syst Mol Biol. 1999, 1999: 18-27.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006.
Lai W, Choudhary V, Park PJ: CGHweb: a tool for comparing DNA copy number segmentations from multiple algorithms. Bioinformatics. 2008, 24 (7): 1014-1015. 10.1093/bioinformatics/btn067.
Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36 (9): 949-951. 10.1038/ng1416.
Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JA: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35: W71-W74. 10.1093/nar/gkm306.
We thank Nick Shera for his early efforts during the data acquisition phase of this project and Ezra Lyon for his bioinformatic support. We also thank NHGRI and NCI for support (D.C.S.; R01- HG000225; R33CA111933); M.R would also like to acknowledge the Morgridge Biotechnology Fellowship for funding.
The authors declare that they have no competing interests.
MR carried out the experimental studies, analysis and manuscript writing. SG developed and applied new bioinformatic tools used for the study. SZ contributed to the data analysis. KP helped with the optical mapping and data interpretation. DS and MAN contributed new statistical tools. EE and CK assisted in the statistical analysis of structural variants. OB contributed samples and overall project guidance. DCS conceived the study, its design, manuscript writing and overall management of the project. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Basic attributes of oligodendroglioma datasets. This spreadsheet describes the two optical maps in detail. Column A lists specific characteristics associated with each step of constructing the map. Columns B and C lists its values for tumor HF087 (column B) and HF1551 (column C). (XLSX 11 KB)
Additional file 2:Overlap (counts) between oligodendroglioma structural variants and genes, segmental duplications, SNPs, and variants reported by other investigators. This table lists counts of structural variants from oligodendroglioma that intersect with genes, segmental duplications, SNPs, and events from the Database of Genomic Variants (DGV). Column A specifies the genomic element; columns B and H indicate total counts for HF087 and HF1551 respectively; Columns C-G and I-M shows overlap counts by variant class. Variants from the DGV are divided based on size into those over 3 kb and those under then 3 kb, then further by study. Variants less than 3 kb are compared to optical map EC and MCs, while those over 3 kb are compared to INS, DEL and OTHER. (XLSX 14 KB)
Additional file 3:Detailed description of intersections between oligodendroglioma structural variants and genes, SNPs, variants from the Database of Genomic Variants and other normal human optical maps. This spreadsheet provides a detailed breakdown of the overlap between Each row in the spreadsheet shows an optical map difference (column A), it’s location (columns B-D), and genes (column E), variants from the Database of Genomic Variants (column F, G), snip-SNPs (column H), and structural variants from other Optical Mapping datasets (columns I-N) that overlap with it. The number is parenthesis after each column header indicates the intersection window. (XLSX 233 KB)
Additional file 4:Experimental validation of oligodendroglioma ECs and MCs by SNP array. This table lists cut differences found in HF087 and HF1551 (column A), their location (columns B-D), type (column E), the SNP genotype corresponding to it (column F) and whether it agrees with the optical map (column G). (XLSX 15 KB)
Additional file 5:Experimental validation of oligodendroglioma indels. This table lists indels found in HF087 and HF1551 (column A), their location (columns B-D) and whether it is validated by a given copy number algorithm (columns E-J). (XLSX 12 KB)
Additional file 6:Oligodendroglioma structural variants and their intersection (counts) with variants detected in six other normal human optical maps. This table lists counts of structural variants from oligodendroglioma that intersect with variants found in other normal human genomes that have been analyzed by Optical Mapping. Column A specifies the genomic element; columns B and H indicate total counts for HF087 and HF1551 respectively; Columns C-G and I-M shows overlap counts by variant class. Only variants of the same category are included in the comparison. (XLSX 11 KB)
Additional file 7:Non-genic candidates found in oligodendroglioma and functional elements from ENCODE that intersect them. Column A lists identifiers for candidate loci that do not occur within a gene, their locations (columns B-D), overlapping transcripts found by GENCODE (column E), and predicted genomic state in different cell types (columns F-N). The cell types in red font are cancer cell lines. Different genomic states are color-coded as per ENCODE website, and is detailed in the key. (XLSX 17 KB)
Additional file 8:Differentially expressed genes in GEO dataset GSE4290, analyzed by Ebarrays, p = 1E-03. Entrez gene identifer (column A), gene symbol (column B), gene name (column C), cellular location (column D) and molecular function (column E) of all differentially expressed genes in GEO dataset GSE4290. (XLSX 199 KB)
Additional file 9:Number of oligodendrogioma samples in REMBRANDT database where a given candidate gene is up or down regulated by at least 2 fold. This table lists candidate genes identified through Optical Mapping (column A), and the number of oligodendroglioma samples in the REMBRANDT database where that gene is up or down regulated by at least two fold (column B). (XLSX 10 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ray, M., Goldstein, S., Zhou, S. et al. Discovery of structural alterations in solid tumor oligodendroglioma by single molecule analysis. BMC Genomics 14, 505 (2013). https://doi.org/10.1186/1471-2164-14-505
- Cancer genomics
- Single molecule
- Optical mapping
- Structural variation