Impact of Rare and Multiple Concurrent Gene Fusions on Diagnostic DNA Methylation Classifier in Brain Tumors.

DNA methylation is an essential molecular assay for central nervous system (CNS) tumor diagnostics. While some fusions define specific brain tumors, others occur across many different diagnoses. We performed a retrospective analysis of 219 primary CNS tumors with whole genome DNA methylation and RNA next-generation sequencing. DNA methylation profiling results were compared with RNAseq detected gene fusions. We detected 105 rare fusions involving 31 driver genes, including 23 fusions previously not implicated in brain tumors. In addition, we identified 6 multi-fusion tumors. Rare fusions and multi-fusion events can impact the diagnostic accuracy of DNA methylation by decreasing confidence in the result, such as BRAF, RAF, or FGFR1 fusions, or result in a complete mismatch, such as NTRK, EWSR1, FGFR, and ALK fusions. IMPLICATIONS: DNA methylation signatures need to be interpreted in the context of pathology and discordant results warrant testing for novel and rare gene fusions.

Analytical demands to use whole-genome sequencing in precision oncology.

Interrogating the tumor genome in its entirety by whole-genome sequencing (WGS) offers an unprecedented insight into the biology and pathogenesis of cancer, with potential impact on diagnostics, prognostication and therapy selection. WGS is able to detect sequence as well as structural variants and thereby combines central domains of cytogenetics and molecular genetics. Given the potential of WGS in directing targeted therapeutics and clinical decision-making, we envision a gradual transition of the method from research to clinical routine. This review is one out of three within this issue aimed at facilitating this effort, by discussing in-depth analytical validation, clinical interpretation and clinical utility of WGS. The review highlights the requirements for implementing, validating and maintaining a clinical WGS pipeline to obtain high-quality patient-specific data in accordance with the local regulatory landscape. Every step of the WGS pipeline, which includes DNA extraction, library preparation, sequencing, bioinformatics analysis, and data storage, is considered with respect to its logistics, necessities, potential pitfalls, and the required quality management. WGS is likely to drive clinical diagnostics and patient care forward, if requirements and challenges of the technique are recognized and met.

Clinical interpretation of whole-genome and whole-transcriptome sequencing for precision oncology.

Whole-genome sequencing either alone or in combination with whole-transcriptome sequencing has started to be used to analyze clinical tumor samples to improve diagnosis, provide risk stratification, and select patient-specific therapies. Compared with current genomic testing strategies, largely focused on small number of genes tested individually or targeted panels, whole-genome and transcriptome sequencing (WGTS) provides novel opportunities to identify and report a potentially much larger number of actionable alterations with diagnostic, prognostic, and/or predictive impact. Such alterations include point mutations, indels, copy- number aberrations and structural variants, but also germline variants, fusion genes, noncoding alterations and mutational signatures. Nevertheless, these comprehensive tests are accompanied by many challenges ranging from the extent and diversity of sequence alterations detected by these methods to the complexity and limited existing standardization in interpreting them. We describe the challenges of WGTS interpretation and the opportunities with comprehensive genomic testing.

RNA Sequencing of Primary Cutaneous and Breast-Implant Associated Anaplastic Large Cell Lymphomas Reveals Infrequent Fusion Transcripts and Upregulation of PI3K/AKT Signaling via Neurotrophin Pathway Genes.

Cutaneous and breast implant-associated anaplastic large-cell lymphomas (cALCLs and BI-ALCLs) are two localized forms of peripheral T-cell lymphomas (PTCLs) that are recognized as distinct entities within the family of ALCL. JAK-STAT signaling is a common feature of all ALCL subtypes, whereas DUSP22/IRF4, TP63 and TYK gene rearrangements have been reported in a proportion of ALK-negative sALCLs and cALCLs. Both cALCLs and BI-ALCLs differ in their gene expression profiles compared to PTCLs; however, a direct comparison of the genomic alterations and transcriptomes of these two entities is lacking. By performing RNA sequencing of 1385 genes (TruSight RNA Pan-Cancer, Illumina) in 12 cALCLs, 10 BI-ALCLs and two anaplastic lymphoma kinase (ALK)-positive sALCLs, we identified the previously reported TYK2-NPM1 fusion in 1 cALCL (1/12, 8%), and four new intrachromosomal gene fusions in 2 BI-ALCLs (2/10, 20%) involving genes on chromosome 1 (EPS15-GNG12 and ARNT-GOLPH3L) and on chromosome 17 (MYO18A-GIT1 and NF1-GOSR1). One of the two BI-ALCL samples showed a complex karyotype, raising the possibility that genomic instability may be responsible for intra-chromosomal fusions in BI-ALCL. Moreover, transcriptional analysis revealed similar upregulation of the PI3K/Akt pathway, associated with enrichment in the expression of neurotrophin signaling genes, which was more conspicuous in BI-ALCL, as well as differences, i.e., over-expression of genes involved in the RNA polymerase II transcription program in BI-ALCL and of the RNA splicing/processing program in cALCL.

Multiregional genetic evolution of metastatic uveal melanoma.

Uveal melanoma (UM) is the most common primary intraocular malignancy in adults and leads to deadly metastases for which there is no approved treatment. Genetic events driving early tumor development are well-described, but those occurring later during metastatic progression remain poorly understood. We performed multiregional genomic sequencing on 22 tumors collected from two patients with widely metastatic UM who underwent rapid autopsy. We observed multiple seeding events from the primary tumors, metastasis-to-metastasis seeding, polyclonal seeding, and late driver variants in ATM, KRAS, and other genes previously unreported in UM. These findings reveal previously unrecognized temporal and anatomic complexity in the genetic evolution of metastatic uveal melanoma, and they highlight the distinction between early and late phases of UM genetic evolution with implications for novel therapeutic approaches.

Metastatic pediatric sclerosing epithelioid fibrosarcoma.

Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft-tissue sarcoma thought to originate in fibroblasts of the tissues comprising tendons, ligaments, and muscles. Minimally responsive to conventional cytotoxic chemotherapies, >50% of SEF patients experience local recurrence and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80%-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-yr-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. Finally, we describe the generation of an SEF primary cell line, the first such culture to be reported. The patient described here experienced persistent disease progression despite aggressive treatment including multiple resections, radiotherapy, and numerous chemotherapies and targeted therapeutics. Untreated and locally recurrent tumor and metastatic tissue were sequenced by whole-genome, whole-exome, and deep-transcriptome next-generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and discuss potential therapeutic strategies for SEF.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Identification of targetable BRAF ΔN486_P490 variant by whole-genome sequencing leading to dabrafenib-induced remission of a BRAF-mutant pancreatic adenocarcinoma.

The tumor genome of a patient with advanced pancreatic cancer was sequenced to identify potential therapeutic targetable mutations after standard of care failed to produce any significant overall response. Matched tumor-normal whole-genome sequencing revealed somatic mutations in BRAF, TP53, CDKN2A, and a focal deletion of SMAD4 The BRAF variant was an in-frame deletion mutation (ΔN486_P490), which had been previously demonstrated to be a kinase-activating alteration in the BRAF kinase domain. Working with the Novartis patient assistance program allowed us to treat the patient with the BRAF inhibitor, dabrafenib. The patient's overall clinical condition improved dramatically with dabrafenib. Levels of serum tumor marker dropped immediately after treatment, and a subsequent CT scan revealed a significant decrease in the size of both primary and metastatic lesions. The dabrafenib-induced remission lasted for 6 mo. Preclinical studies published concurrently with the patient's treatment showed that the BRAF in-frame mutation (ΔNVTAP) induces oncogenic activation by a mechanism distinct from that induced by V600E, and that this difference dictates the responsiveness to different BRAF inhibitors. This study describes a dramatic instance of how high-level genomic technology and analysis was necessary and sufficient to identify a clinically logical treatment option that was then utilized and shown to be of clinical value for this individual.

Correction to: Sequencing and curation strategies for identifying candidate glioblastoma treatments.

Following publication of the original article [1], it was reported that the given name of the fourteenth author was incorrectly published. The incorrect and the correct names are given below.

Sequencing and curation strategies for identifying candidate glioblastoma treatments.

BACKGROUND: Prompted by the revolution in high-throughput sequencing and its potential impact for treating cancer patients, we initiated a clinical research study to compare the ability of different sequencing assays and analysis methods to analyze glioblastoma tumors and generate real-time potential treatment options for physicians. METHODS: A consortium of seven institutions in New York City enrolled 30 patients with glioblastoma and performed tumor whole genome sequencing (WGS) and RNA sequencing (RNA-seq; collectively WGS/RNA-seq); 20 of these patients were also analyzed with independent targeted panel sequencing. We also compared results of expert manual annotations with those from an automated annotation system, Watson Genomic Analysis (WGA), to assess the reliability and time required to identify potentially relevant pharmacologic interventions. RESULTS: WGS/RNAseq identified more potentially actionable clinical results than targeted panels in 90% of cases, with an average of 16-fold more unique potentially actionable variants identified per individual; 84 clinically actionable calls were made using WGS/RNA-seq that were not identified by panels. Expert annotation and WGA had good agreement on identifying variants [mean sensitivity = 0.71, SD = 0.18 and positive predictive value (PPV) = 0.80, SD = 0.20] and drug targets when the same variants were called (mean sensitivity = 0.74, SD = 0.34 and PPV = 0.79, SD = 0.23) across patients. Clinicians used the information to modify their treatment plan 10% of the time. CONCLUSION: These results present the first comprehensive comparison of technical and machine augmented analysis of targeted panel and WGS/RNA-seq to identify potential cancer treatments.

A recurrent novel MGA-NUTM1 fusion identifies a new subtype of high-grade spindle cell sarcoma.

NUTM1-rearranged tumors are defined by the presence of a gene fusion between NUTM1 and various gene partners and typically follow a clinically aggressive disease course with poor outcomes despite conventional multimodality therapy. NUTM1-rearranged tumors display histologic features of a poorly differentiated carcinoma with areas of focal squamous differentiation and typically express the BRD4-NUTM1 fusion gene defining a distinct clinicopathologic entity-NUT carcinoma (NC). NCs with mesenchymal differentiation have rarely been described in the literature. In this report, we describe the characterization of two cases of high-grade spindle cell sarcoma harboring a novel MGA-NUTM1 fusion. Whole-genome sequencing identified the presence of complex rearrangements resulting in a MGA-NUTM1 fusion gene in the absence of other significant somatic mutations. Genetic rearrangement was confirmed by fluorescence in situ hybridization, and expression of the fusion gene product was confirmed by transcriptomic analysis. The fusion protein was predicted to retain nearly the entire protein sequence of both MGA (exons 1-22) and NUTM1 (exons 3-8). Histopathologically, both cases were high-grade spindle cell sarcomas without specific differentiation markers. In contrast to typical cases of NC, these cases were successfully treated with aggressive local control measures (surgery and radiation) and both patients remain alive without disease. These cases describe a new subtype of NUTM1-rearranged tumors warranting expansion of diagnostic testing to evaluate for the presence of MGA-NUTM1 or alternative NUTM1 gene fusions in the diagnostic workup of high-grade spindle cell sarcomas or small round blue cell tumors of ambiguous lineage.

Analytical Validation of Clinical Whole-Genome and Transcriptome Sequencing of Patient-Derived Tumors for Reporting Targetable Variants in Cancer.

We developed and validated a clinical whole-genome and transcriptome sequencing (WGTS) assay that provides a comprehensive genomic profile of a patient's tumor. The ability to fully capture the mappable genome with sufficient sequencing coverage to precisely call DNA somatic single nucleotide variants, insertions/deletions, copy number variants, structural variants, and RNA gene fusions was analyzed. New York State's Department of Health next-generation DNA sequencing guidelines were expanded for establishing performance validation applicable to whole-genome and transcriptome sequencing. Whole-genome sequencing laboratory protocols were validated for the Illumina HiSeq X Ten platform and RNA sequencing for Illumina HiSeq2500 platform for fresh or frozen and formalin-fixed, paraffin-embedded tumor samples. Various bioinformatics tools were also tested, and CIs for sensitivity and specificity thresholds in calling clinically significant somatic aberrations were determined. The validation was performed on a set of 125 tumor normal pairs. RNA sequencing was performed to call fusions and to confirm the DNA variants or exonic alterations. Here, we present our results and WGTS standards for variant allele frequency, reproducibility, analytical sensitivity, and present limit of detection analysis for single nucleotide variant calling, copy number identification, and structural variants. We show that The New York Genome Center WGTS clinical assay can provide a comprehensive patient variant discovery approach suitable for directed oncologic therapeutic applications.

Whole Genome Sequencing-Based Discovery of Structural Variants in Glioblastoma.

Next-generation DNA sequencing (NGS) technologies are currently being applied in both research and clinical settings for the understanding and management of disease. The goal is to use high-throughput sequencing to identify specific variants that drive tumorigenesis within each individual's tumor genomic profile. The significance of copy number and structural variants in glioblastoma makes it essential to broaden the search beyond oncogenic single nucleotide variants toward whole genome profiles of genetic aberrations that may contribute to disease progression. The heterogeneity of glioblastoma and its variability of cancer driver mutations necessitate a more robust examination of a patient's tumor genome. Here, we present patient whole genome sequencing (WGS) information to identify oncogenic structural variants that may contribute to glioblastoma pathogenesis. We provide WGS protocols and bioinformatics approaches to identify copy number and structural variations in 41 glioblastoma patient samples. We present how WGS can identify structural diversity within glioblastoma samples. We specifically show how to apply current bioinformatics tools to detect EGFR variants and other structural aberrations from DNA whole genome sequencing and how to validate those variants within the laboratory. These comprehensive WGS protocols can provide additional information directing more precise therapeutic options in the treatment of glioblastoma.

Pancreatic intraductal tubulopapillary neoplasm is genetically distinct from intraductal papillary mucinous neoplasm and ductal adenocarcinoma.

Intraductal tubulopapillary neoplasm is a relatively recently described member of the pancreatic intraductal neoplasm family. The more common member of this family, intraductal papillary mucinous neoplasm, often carries genetic alterations typical of pancreatic infiltrating ductal adenocarcinoma (KRAS, TP53, and CDKN2A) but additionally has mutations in GNAS and RNF43 genes. However, the genetic characteristics of intraductal tubulopapillary neoplasm have not been well characterized. Twenty-two intraductal tubulopapillary neoplasms were analyzed by either targeted next-generation sequencing, which enabled the identification of sequence mutations, copy number alterations, and selected structural rearrangements involving all targeted (≥300) genes, or whole-exome sequencing. Three of these intraductal tubulopapillary neoplasms were also subjected to whole-genome sequencing. All intraductal tubulopapillary neoplasms revealed the characteristic histologic (cellular intraductal nodules of back-to-back tubular glands lined by predominantly cuboidal cells with atypical nuclei and no obvious intracellular mucin) and immunohistochemical (immunolabeled with MUC1 and MUC6 but were negative for MUC2 and MUC5AC) features. By genomic analyses, there was loss of CDKN2A in 5/20 (25%) of these cases. However, the majority of the previously reported intraductal papillary mucinous neoplasm-related alterations were absent. Moreover, in contrast to most ductal neoplasms of the pancreas, MAP-kinase pathway was not involved. In fact, 2/22 (9%) of intraductal tubulopapillary neoplasms did not reveal any mutations in the tested genes. However, certain chromatin remodeling genes (MLL1, MLL2, MLL3, BAP1, PBRM1, EED, and ATRX) were found to be mutated in 7/22 (32%) of intraductal tubulopapillary neoplasms and 27% harbored phosphatidylinositol 3-kinase (PI3K) pathway (PIK3CA, PIK3CB, INPP4A, and PTEN) mutations. In addition, 4/18 (18%) of intraductal tubulopapillary neoplasms had FGFR2 fusions (FGFR2-CEP55, FGFR2-SASS6, DISP1-FGFR2, FGFR2-TXLNA, and FGFR2-VCL) and 1/18 (5.5%) had STRN-ALK fusion. Intraductal tubulopapillary neoplasm is a distinct clinicopathologic entity in the pancreas. Although its intraductal nature and some clinicopathologic features resemble those of intraductal papillary mucinous neoplasm, our results suggest that intraductal tubulopapillary neoplasm has distinguishing genetic characteristics. Some of these mutated genes are potentially targetable. Future functional studies will be needed to determine the consequences of these gene alterations.

Comparing sequencing assays and human-machine analyses in actionable genomics for glioblastoma.

OBJECTIVE: To analyze a glioblastoma tumor specimen with 3 different platforms and compare potentially actionable calls from each. METHODS: Tumor DNA was analyzed by a commercial targeted panel. In addition, tumor-normal DNA was analyzed by whole-genome sequencing (WGS) and tumor RNA was analyzed by RNA sequencing (RNA-seq). The WGS and RNA-seq data were analyzed by a team of bioinformaticians and cancer oncologists, and separately by IBM Watson Genomic Analytics (WGA), an automated system for prioritizing somatic variants and identifying drugs. RESULTS: More variants were identified by WGS/RNA analysis than by targeted panels. WGA completed a comparable analysis in a fraction of the time required by the human analysts. CONCLUSIONS: The development of an effective human-machine interface in the analysis of deep cancer genomic datasets may provide potentially clinically actionable calls for individual patients in a more timely and efficient manner than currently possible. CLINICALTRIALSGOV IDENTIFIER: NCT02725684.

The oncocytic subtype is genetically distinct from other pancreatic intraductal papillary mucinous neoplasm subtypes.

In 2010, the World Health Organization reclassified the entity originally described as intraductal oncocytic papillary neoplasm as the 'oncocytic subtype' of intraductal papillary mucinous neoplasm. Although several key molecular alterations of other intraductal papillary mucinous neoplasm subtypes have been discovered, including common mutations in KRAS, GNAS, and RNF3, those of oncocytic subtype have not been well characterized. We analyzed 11 pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms. Nine pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms uniformly exhibited typical entity-defining morphology of arborizing papillae lined by layers of cells with oncocytic cytoplasm, prominent, nucleoli, and intraepithelial lumina. The remaining two were atypical. One lacked the arborizing papilla and had flat oncocytic epithelium only; the other one had focal oncocytic epithelium in a background of predominantly intestinal subtype intraductal papillary mucinous neoplasm. Different components of this case were analyzed separately. Formalin-fixed, paraffin-embedded specimens of all cases were microdissected and subjected to high-depth-targeted next-generation sequencing for a panel of 300 key cancer-associated genes in a platform that enabled the identification of sequence mutations, copy number alterations, and select structural rearrangements involving all targeted genes. Fresh frozen specimens of two cases were also subjected to whole-genome sequencing. For the nine typical pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms, the number of mutations per case, identified by next-generation sequencing, ranged from 1 to 10 (median=4). None of these cases had KRAS or GNAS mutations and only one had both RNF43 and PIK3R1 mutations. ARHGAP26, ASXL1, EPHA8, and ERBB4 genes were somatically altered in more than one of these typical 'oncocytic subtype' of intraductal papillary mucinous neoplasms but not in the other two atypical ones. In the neoplasm with flat oncocytic epithelium, the only mutated gene was KRAS. All components of the intestinal subtype intraductal papillary mucinous neoplasms with focal oncocytic epithelium manifested TP53, GNAS, and RNF43 mutations. In conclusion, this study elucidates that 'oncocytic subtype' of intraductal papillary mucinous neoplasm is not only morphologically distinct but also genetically distinct from other intraductal papillary mucinous neoplasm subtypes. Considering that now its biologic behavior is also being found to be different than other intraductal papillary mucinous neoplasm subtypes, 'oncocytic subtype' of intraductal papillary mucinous neoplasm warrants being recognized separately.

A quantitative proteomics-based signature of platinum sensitivity in ovarian cancer cell lines.

Although DNA encodes the molecular instructions that underlie the control of cell function, it is the proteins that are primarily responsible for implementing those instructions. Therefore quantitative analyses of the proteome would be expected to yield insights into important candidates for the detection and treatment of disease. We present an iTRAQ (isobaric tag for relative and absolute quantification)-based proteomic analysis of ten ovarian cancer cell lines and two normal ovarian surface epithelial cell lines. We profiled the abundance of 2659 cellular proteins of which 1273 were common to all 12 cell lines. Of the 1273, 75 proteins exhibited elevated expression and 164 proteins had diminished expression in the cancerous cells compared with the normal cell lines. The iTRAQ expression profiles allowed us to segregate cell lines based upon sensitivity and resistance to carboplatin. Importantly, we observed no substantial correlation between protein abundance and RNA expression or epigenetic DNA methylation data. Furthermore, we could not discriminate between sensitivity and resistance to carboplatin on the basis of RNA expression and DNA methylation data alone. The present study illustrates the importance of proteomics-based discovery for defining the basis for the carboplatin response in ovarian cancer and highlights candidate proteins, particularly involved in cellular redox regulation, homologous recombination and DNA damage repair, which otherwise could not have been predicted from whole genome and expression data sources alone.

Integrative prediction of gene function and platinum-free survival from genomic and epigenetic features in ovarian cancer.

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to discover genes critical to the development, progression, and therapeutic resistance of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatics analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We identify changes in DNA methylation and expression specifically for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic gene function from integrative analysis of three modalities: copy number variation, DNA methylation, and gene expression. Our method (1) calculates the extent of genomic and epigenetic alterations of defined tumor suppressor and oncogenic features for the functional prediction of significant ovarian cancer gene candidates and (2) identifies the functional activity or inactivity of known tumor suppressors and oncogenes in ovarian cancer. We applied our protocol on 42 primary serous ovarian cancer samples using MOMA-ROMA representational array assays. Additionally, we provide the basis for incorporating epigenetic profiles of ovarian tumors for the purposes of platinum-free survival prediction in the context of TCGA data.

Loss of DOK2 induces carboplatin resistance in ovarian cancer via suppression of apoptosis.

OBJECTIVE: Ovarian cancers are highly heterogeneous and while chemotherapy is the preferred treatment many patients are intrinsically resistant or quickly develop resistance. Furthermore, all tumors that recur ultimately become resistant. Recent evidence suggests that epigenetic deregulation may be a key factor in the onset and maintenance of chemoresistance. We set out to identify epigenetically silenced genes that affect chemoresistance. METHODS: The epigenomes of a total of 45 ovarian samples were analyzed to identify epigenetically altered genes that segregate with platinum response, and further filtered with expression data to identify genes that were suppressed. A tissue culture carboplatin resistance screen was utilized to functionally validate this set of candidate platinum resistance genes. RESULTS: Our screen correctly identified 19 genes that when suppressed altered the chemoresistance of the cells in culture. Of the genes identified in the screen we further characterized one gene, docking protein 2 (DOK2), an adapter protein downstream of tyrosine kinase, to determine if we could elucidate the mechanism by which it increased resistance. The loss of DOK2 decreased the level of apoptosis in response to carboplatin. Furthermore, in cells with reduced DOK2, the level of anoikis was decreased. CONCLUSIONS: We have developed a screening methodology that analyzes the epigenome and informatically identifies candidate genes followed by in vitro culture screening of the candidate genes. To validate our screening methodology we further characterized one candidate gene, DOK2, and showed that loss of DOK2 induces chemotherapy resistance by decreasing the level of apoptosis in response to treatment.

Identification of tumor suppressors and oncogenes from genomic and epigenetic features in ovarian cancer.

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to identify genes critical to the development and progression of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatic analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We separately examined CNV and DNA methylation for 42 primary serous ovarian cancer samples using MOMA-ROMA assays and 379 tumor samples analyzed by The Cancer Genome Atlas. We have identified 346 genes with significant deletions or amplifications among the tumor samples. Utilizing associated gene expression data we predict 156 genes with altered copy number and correlated changes in expression. Among these genes CCNE1, POP4, UQCRB, PHF20L1 and C19orf2 were identified within both data sets. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We therefore identify changes in DNA methylation and expression for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic features for these modalities and perform a correlation analysis with expression. We predicted 611 potential oncogenes and tumor suppressors candidates by integrating these data types. Genes with a strong correlation for methylation dependent expression changes exhibited at varying copy number aberrations include CDCA8, ATAD2, CDKN2A, RAB25, AURKA, BOP1 and EIF2C3. We provide copy number variation and DNA methylation analysis for over 11,500 individual genes covering the genetic landscape of ovarian cancer tumors. We show the extent of genomic and epigenetic alterations for known tumor suppressors and oncogenes and also use these defined features to identify potential ovarian cancer gene candidates.