Exome sequencing provides unprecedented insights into cancer biology and pharmacological response. Here we assess these two parameters for the NCI-60, which is among the richest genomic and pharmacological publicly available cancer cell line databases. Homozygous genetic variants that putatively affect protein function were identified in 1,199 genes (approximately 6% of all genes). Variants that are either enriched or depleted compared to non-cancerous genomes, and thus may be influential in cancer progression and differential drug response were identified for 2,546 genes. Potential gene knockouts are made available. Assessment of cell line response to 19,940 compounds, including 110 FDA-approved drugs, reveals ≈80-fold range in resistance versus sensitivity response across cell lines. 103,422 gene variants were significantly correlated with at least one compound (at p<0.0002). These include genes of known pharmacological importance such as IGF1R, BRAF, RAD52, MTOR, STAT2 and TSC2 as well as a large number of candidate genes such as NOM1, TLL2, and XDH. We introduce two new web-based CellMiner applications that enable exploration of variant-to-compound relationships for a broad range of researchers, especially those without bioinformatics support. The first tool, "Genetic variant versus drug visualization", provides a visualization of significant correlations between drug activity-gene variant combinations. Examples are given for the known vemurafenib-BRAF, and novel ifosfamide-RAD52 pairings. The second, "Genetic variant summation" allows an assessment of cumulative genetic variations for up to 150 combined genes together; and is designed to identify the variant burden for molecular pathways or functional grouping of genes. An example of its use is provided for the EGFR-ERBB2 pathway gene variant data and the identification of correlated EGFR, ERBB2, MTOR, BRAF, MEK and ERK inhibitors. The new tools are implemented as an updated web-based CellMiner version, for which the present publication serves as a compendium.

The genetic profile of human pancreatic cancers harbors considerable heterogeneity, which suggests a possible explanation for the pronounced inefficacy of single therapies in this disease. This observation has led to a belief that custom therapies based on individual tumor profiles are necessary to more effectively treat pancreatic cancer. It has recently been discovered that axon guidance genes are affected by somatic structural variants in up to 25% of human pancreatic cancers. Thus far, however, some of these mutations have only been correlated to survival probability and no function has been assigned to these observed axon guidance gene mutations in pancreatic cancer. In this study we established three novel pancreatic cancer cell lines and performed whole genome sequencing to discover novel mutations in axon guidance genes that may contribute to the cancer phenotype of these cells. We discovered, among other novel somatic variants in axon guidance pathway genes, a novel mutation in the PLXNA1 receptor (c.2587G>A) in newly established cell line SB.06 that mediates oncogenic cues of increased invasion and proliferation in SB.06 cells and increased invasion in 293T cells upon stimulation with the receptor's natural ligand semaphorin 3A compared to wild type PLXNA1 cells. Mutant PLXNA1 signaling was associated with increased Rho-GTPase and p42/p44 MAPK signaling activity and cytoskeletal expansion, but not changes in E-cadherin, vimentin, or metalloproteinase 9 expression levels. Pharmacologic inhibition of the Rho-GTPase family member CDC42 selectively abrogated PLXNA1 c.2587G>A-mediated increased invasion. These findings provide in-vitro confirmation that somatic mutations in axon guidance genes can provide oncogenic gain-of-function signals and may contribute to pancreatic cancer progression.

Hairy cell leukemia (HCL) is a purine analog-responsive B-cell malignancy containing the BRAF V600E mutation, expressing CD22, CD11c, CD103, tartrate resistant acid phosphatase (TRAP) CD25, CD123, and annexin 1A. BRAF V600E and the latter 4 markers are usually absent in the more aggressive and chemoresistant variant HCLv. To evaluate differences between HCL and HCLv, expression microarrays comparing HCL with HCLv were performed for 24694 genes using 47323 probes. Microarray data from 35 HCL and 27 HCLv purified samples showed the greatest HCL-HCLv difference in the muscle-associated gene MYF6, expressed by its 2 probes 18.5- and 10.8-fold higher in HCL than HCLv (p<0.0001). By real-time quantitative PCR (RQ-PCR), 100% of 152 classic HCL samples were MYF6-positive, vs 5 (6%) of 90 blood donors. MYF6-expression was also detected in 18 (35%) of 51 with HCLv, 11 (92%) of 12 with HCL expressing unmutated IGHV4-34, 35 (73%) of 48 with chronic lymphocytic leukemia (CLL), and 1 (8%) of 12 with mantle cell lymphoma. Hypomethylation status of MYF6 supported expression in HCL more than HCLv. Posttreatment blood samples becoming negative by flow cytometry remained MYF6+ by RQ-PCR in 42 (48%) of 87 HCL patients, and MYF6 RQ-PCR could detect 1 HCL in 105 normal cells. MYF6, universally expressed in HCL and in most CLL samples, may be a useful biomarker for these leukemias. Further studies are underway to determine the role of MYF6 in HCL.

A new ovarian near-diploid cell line, OVDM1, was derived from a highly aneuploid serous ovarian metastatic adenocarcinoma. A metastatic tumor was obtained from a 47-year-old Ashkenazi Jewish patient three years after the first surgery removed the primary tumor, both ovaries, and the remaining reproductive organs. OVDM1 was characterized by cell morphology, genotyping, tumorigenic assay, mycoplasma testing, spectral karyotyping (SKY), and molecular profiling of the whole genome by aCGH and gene expression microarray. Targeted sequencing of a panel of cancer-related genes was also performed. Hierarchical clustering of gene expression data clearly confirmed the ovarian origin of the cell line. OVDM1 has a near-diploid karyotype with a low-level aneuploidy, but samples of the original metastatic tumor were grossly aneuploid. A number of single nucleotide variations (SNVs)/mutations were detected in OVDM1 and the metastatic tumor samples. Some of them were cancer-related according to COSMIC and HGMD databases (no founder mutations in BRCA1 and BRCA2 have been found). A large number of focal copy number alterations (FCNAs) were detected, including homozygous deletions (HDs) targeting WWOX and GATA4. Progression of OVDM1 from early to late passages was accompanied by preservation of the near-diploid status, acquisition of only few additional large chromosomal rearrangements and more than 100 new small FCNAs. Most of newly acquired FCNAs seem to be related to localized but massive DNA fragmentation (chromothripsis-like rearrangements). Newly developed near-diploid OVDM1 cell line offers an opportunity to evaluate tumorigenesis pathways/events in a minor clone of metastatic ovarian adenocarcinoma as well as mechanisms of chromothripsis.

In contrast to adult mutant gastrointestinal stromal tumors [GISTs], pediatric/wild-type GISTs remain poorly understood overall, given their lack of oncogenic activating tyrosine kinase mutations. These GISTs, with a predilection for gastric origin in female patients, show limited response to therapy with tyrosine kinase inhibitors and generally pursue a more indolent course, but still may prove fatal. Defective cellular respiration appears to underpin tumor development in these wild-type cases, which as a group lack expression of succinate dehydrogenase [SDH] B, a surrogate marker for respiratory chain metabolism. Yet, only a small subset of the wild-type tumors show mutations in the genes coding for the SDH subunits [SDHx]. To explore additional pathogenetic mechanisms in these wild-type GISTs, we elected to investigate post-transcriptional regulation of these tumors by conducting microRNA (miRNA) profiling of a mixed cohort of 73 cases including 18 gastric pediatric wild-type, 25 (20 gastric, 4 small bowel and 1 retroperitoneal) adult wild-type GISTs and 30 gastric adult mutant GISTs. By this approach we have identified distinct signatures for GIST subtypes which correlate tightly with clinico-pathological parameters. A cluster of miRNAs on 14q32 show strikingly different expression patterns amongst GISTs, a finding which appears to be explained at least in part by differential allelic methylation of this imprinted region. Small bowel and retroperitoneal wild-type GISTs segregate with adult mutant GISTs and express SDHB, while adult wild-type gastric GISTs are dispersed amongst adult mutant and pediatric wild-type cases, clustering in this situation on the basis of SDHB expression. Interestingly, global methylation analysis has recently similarly demonstrated that these wild-type, SDHB-immunonegative tumors show a distinct pattern compared with KIT and PDGFRA mutant tumors, which as a rule do express SDHB. All cases with Carney triad within our cohort cluster together tightly.

Chordoma is a rare primary bone malignancy that arises in the skull base, spine and sacrum and originates from remnants of the notochord. These tumors are typically resistant to conventional chemotherapy, and to date there are no FDA-approved agents to treat chordoma. The lack of in vivo models of chordoma has impeded the development of new therapies for this tumor. Primary tumor from a sacral chordoma was xenografted into NOD/SCID/IL-2R γ-null mice. The xenograft is serially transplantable and was characterized by both gene expression analysis and whole genome SNP genotyping. The NIH Chemical Genomics Center performed high-throughput screening of 2,816 compounds using two established chordoma cell lines, U-CH1 and U-CH2B. The screen yielded several compounds that showed activity and two, sunitinib and bortezomib, were tested in the xenograft. Both agents slowed the growth of the xenograft tumor. Sensitivity to an inhibitor of IκB, as well as inhibition of an NF-κB gene expression signature demonstrated the importance of NF-κB signaling for chordoma growth. This serially transplantable chordoma xenograft is thus a practical model to study chordomas and perform in vivo preclinical drug testing.

Replication of mammalian genomes requires the activation of thousands of origins which are both spatially and temporally regulated by as yet unknown mechanisms. At the most fundamental level, our knowledge about the distribution pattern of origins in each of the chromosomes, among different cell types, and whether the physiological state of the cells alters this distribution is at present very limited.