Translational genomics research in cancers, e.g., International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), has generated large multidimensional datasets from high-throughput technologies. Data analysis at multidimensional level will greatly benefit clinical applications of genomic information in diagnosis, prognosis and therapeutics of cancers. To help, tools to effectively visualize integrated multidimensional data are important for understanding and describing the relationship between genomic variations and cancers.

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.

Diffuse intrinsic pontine glioma (DIPG) is a fatal childhood cancer. We performed a chemical screen in patient-derived DIPG cultures along with RNA-seq analyses and integrated computational modeling to identify potentially effective therapeutic strategies. The multi-histone deacetylase inhibitor panobinostat demonstrated therapeutic efficacy both in vitro and in DIPG orthotopic xenograft models. Combination testing of panobinostat and the histone demethylase inhibitor GSK-J4 revealed that the two had synergistic effects. Together, these data suggest a promising therapeutic strategy for DIPG.

Despite good initial response to chemotherapy, 30% of Ewing's sarcoma (EWS) patients with localised tumours develop recurrent disease, associated with poor prognosis. The aim of this study was to address this challenge by conducting preclinical evaluation of a death receptor targeted agent in vitro and in vivo, and by identifying predictive biomarkers.

Genetic alterations and etiology of thymic epithelial tumors (TETs) are largely unknown, hampering the development of effective targeted therapies for patients with TETs. Here TETs of advanced-stage patients enrolled in a clinical trial of molecularly-guided targeted therapies were employed for targeted sequencing of 197 cancer-associated genes. Comparative sequence analysis of 78 TET/blood paired samples obtained from 47 thymic carcinoma (TC) and 31 thymoma patients revealed a total of 86 somatic non-synonymous sequence variations across 39 different genes in 33 (42%) TETs. TCs (62%; 29/47) showed higher incidence of somatic non-synonymous mutations than thymomas (13%; 4/31; p < 0.0001). TP53 was the most frequently mutated gene in TETs (n = 13; 17%), especially in TCs (26%), and was associated with a poorer overall survival (p < 0.0001). Genes in histone modification [BAP1 (n = 6; 13%), SETD2 (n = 5; 11%), ASXL1 (n = 2; 4%)], chromatin remodeling [SMARCA4 (n = 2; 4%)], and DNA methylation [DNMT3A (n = 3; 7%), TET2 (n = 2; 4%), WT1 (n = 2; 4%)] pathways were recurrently mutated in TCs, but not in thymomas. Our results suggest a potential disruption of epigenetic homeostasis in TCs, and a substantial difference in genetic makeup between TCs and thymomas. Further investigation is warranted into the roles of epigenetic dysregulation in TC development and its potential for targeted therapy.

Papillary renal cell carcinomas (PRCC) are a histologically and genetically heterogeneous group of tumors that represent 15-20% of all kidney neoplasms and may require diverse therapeutic approaches. Alteration of the NF2 tumor suppressor gene, encoding a key regulator of the Hippo signaling pathway, is observed in 22.5% of PRCC. The Hippo signaling pathway controls cell proliferation by regulating the transcriptional activity of Yes-Associated Protein, YAP1. Loss of NF2 results in aberrant YAP1 activation. The Src family kinase member Yes also regulates YAP1 transcriptional activity. This study investigated the importance of YAP and Yes activity in three NF2-deficient PRCC cell lines. NF2-deficency correlated with increased expression of YAP1 transcriptional targets and siRNA-based knockdown of YAP1 and Yes1 downregulated this pathway and dramatically reduced cell viability. Dasatinib and saracatinib have potent inhibitory effects on Yes and treatment with either resulted in downregulation of YAP1 transcription targets, reduced cell viability, and G0-G1 cell cycle arrest. Xenograft models for NF2-deficient PRCC also demonstrated reduced tumor growth in response to dasatinib. Thus, inhibiting Yes and the subsequent transcriptional activity of YAP1 had a substantial anti-tumor cell effect both in vitro and in vivo and may provide a viable therapeutic approach for patients with NF2-deficient PRCC.

Metastasis is characterized by spreading of neoplastic cells to an organ other than where they originated and is the predominant cause of death among cancer patients. This holds true for melanoma, whose incidence is increasing more rapidly than any other cancer and once disseminated has few therapeutic options. Here we performed whole exome sequencing of two sets of matched normal and metastatic tumor DNAs.

To understand the genetic mechanisms driving variant and IGHV4-34-expressing hairy-cell leukemias, we performed whole-exome sequencing of leukemia samples from ten affected individuals, including six with matched normal samples. We identified activating mutations in the MAP2K1 gene (encoding MEK1) in 5 of these 10 samples and in 10 of 21 samples in a validation set (overall frequency of 15/31), suggesting potential new strategies for treating individuals with these diseases.

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.

Genetic mutations, metabolic dysfunction, and epigenetic misregulation are commonly considered to play distinct roles in tumor development and maintenance. However, intimate relationships between these mechanisms are now emerging. In particular, mutations in genes for the core metabolic enzymes IDH, SDH, and FH are significant drivers of diverse tumor types. In each case, the resultant accumulation of particular metabolites inhibits TET enzymes responsible for oxidizing 5-methylcytosine, leading to pervasive DNA hypermethylation.

Intrahepatic cholangiocarcinoma (ICC) and hepatocellular carcinoma (HCC) are clinically disparate primary liver cancers with etiological and biological heterogeneity. We identified common molecular subtypes linked to similar prognosis among 199 Thai ICC and HCC patients through systems integration of genomics, transcriptomics, and metabolomics. While ICC and HCC share recurrently mutated genes, including TP53, ARID1A, and ARID2, mitotic checkpoint anomalies distinguish the C1 subtype with key drivers PLK1 and ECT2, whereas the C2 subtype is linked to obesity, T cell infiltration, and bile acid metabolism. These molecular subtypes are found in 582 Asian, but less so in 265 Caucasian patients. Thus, Asian ICC and HCC, while clinically treated as separate entities, share common molecular subtypes with similar actionable drivers to improve precision therapy.

Immortal tumor cell lines are an important model system for cancer research, however, misidentification and cross-contamination of cell lines are a common problem. Seven chordoma cell lines are reported in the literature, but none has been characterized in detail. We analyzed gene expression patterns and genomic copy number variations in five putative chordoma cell lines (U-CH1, CCL3, CCL4, GB60, and CM319). We also created a new chordoma cell line, U-CH2, and provided genotypes for cell lines for identity confirmation. Our analyses revealed that CCL3, CCL4, and GB60 are not chordoma cell lines, and that CM319 is a cancer cell line possibly derived from chordoma, but lacking expression of key chordoma biomarkers. U-CH1 and U-CH2 both have gene expression profiles, copy number aberrations, and morphology consistent with chordoma tumors. These cell lines also harbor genetic changes, such as loss of p16, MTAP, or PTEN, that make them potentially useful models for studying mechanisms of chordoma pathogenesis and for evaluating targeted therapies.

Estrogen signaling plays a critical role in a number of normal physiological processes and has important implications in the treatment of breast cancer. The p160 nuclear receptor coactivator, AIB1 (amplified in breast cancer 1), is frequently amplified and overexpressed in human breast cancer and has been shown to enhance estrogen-dependent transactivation.

Skin lies at the interface between the complex physiology of the body and the external environment. This essential epidermal barrier, composed of cornified proteins encased in lipids, prevents both water loss and entry of infectious or toxic substances. We uncover that the transcription factor GATA-3 is required to establish the epidermal barrier and survive in the ex utero environment. Analysis of Gata-3 mutant transcriptional profiles at three critical developmental stages identifies a specific defect in lipid biosynthesis and a delay in differentiation. Genomic analysis identifies highly conserved GATA-3 binding sites bound in vivo by GATA-3 in the first intron of the lipid acyltransferase gene AGPAT5. Skin from both Gata-3-/- and previously characterized barrier-deficient Kruppel-like factor 4-/- newborns up-regulate antimicrobial peptides, effectors of innate immunity. Comparison of these animal models illustrates how impairment of the skin barrier by two genetically distinct mechanisms leads to innate immune responses, as observed in the common human skin disorders psoriasis and atopic dermatitis.

Thrombocytopenia-absent radii (TAR) syndrome, characterized by neonatal thrombocytopenia and bilateral radial aplasia with thumbs present, is typically caused by the inheritance of a 1q21.1 deletion and a single-nucelotide polymorphism in RBM8A on the nondeleted allele. We evaluated two siblings with TAR-like dysmorphology but lacking thrombocytopenia in infancy. Family NCI-107 participated in an IRB-approved cohort study and underwent comprehensive clinical and genomic evaluations, including aCGH, whole-exome, whole-genome, and targeted sequencing. Gene expression assays and electromobility shift assays (EMSAs) were performed to evaluate the variant of interest. The previously identified TAR-associated 1q21.1 deletion was present in the affected siblings and one healthy parent. Multiple sequencing approaches did not identify previously described TAR-associated SNPs or mutations in relevant genes. We discovered rs61746197 A > G heterozygosity in the parent without the deletion and apparent hemizygosity in both siblings. rs61746197 A > G overlaps a RelA-p65 binding motif, and EMSAs indicate the A allele has higher transcription factor binding efficiency than the G allele. Stimulation of K562 cells to induce megakaryocyte differentiation abrogated the shift of both reference and alternative probes. The 1q21.1 TAR-associated deletion in combination with the G variant of rs61746197 on the nondeleted allele is associated with a TAR-like phenotype. rs61746197 G could be a functional enhancer/repressor element, but more studies are required to identify the specific factor(s) responsible. Overall, our findings suggest a role of rs61746197 A > G and human disease in the setting of a 1q21.1 deletion on the other chromosome.

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 small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have an overall poorer outcomes, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of the tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyV-positive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyV-positive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomology-mediated end joining and coincided with focal amplifications of the tumor genome similar to other virus-mediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases, with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with an antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virus-associated.

Gastric adenocarcinoma and proximal polyposis of the stomach (GAPPS) is an autosomal-dominant cancer-predisposition syndrome with a significant risk of gastric, but not colorectal, adenocarcinoma. We mapped the gene to 5q22 and found loss of the wild-type allele on 5q in fundic gland polyps from affected individuals. Whole-exome and -genome sequencing failed to find causal mutations but, through Sanger sequencing, we identified point mutations in APC promoter 1B that co-segregated with disease in all six families. The mutations reduced binding of the YY1 transcription factor and impaired activity of the APC promoter 1B in luciferase assays. Analysis of blood and saliva from carriers showed allelic imbalance of APC, suggesting that these mutations lead to decreased allele-specific expression in vivo. Similar mutations in APC promoter 1B occur in rare families with familial adenomatous polyposis (FAP). Promoter 1A is methylated in GAPPS and sporadic FGPs and in normal stomach, which suggests that 1B transcripts are more important than 1A in gastric mucosa. This might explain why all known GAPPS-affected families carry promoter 1B point mutations but only rare FAP-affected families carry similar mutations, the colonic cells usually being protected by the expression of the 1A isoform. Gastric polyposis and cancer have been previously described in some FAP-affected individuals with large deletions around promoter 1B. Our finding that GAPPS is caused by point mutations in the same promoter suggests that families with mutations affecting the promoter 1B are at risk of gastric adenocarcinoma, regardless of whether or not colorectal polyps are present.

Osteosarcoma is the most common primary bone malignancy of adolescents and young adults. To better understand the genetic etiology of osteosarcoma, we performed a multistage genome-wide association study consisting of 941 individuals with osteosarcoma (cases) and 3,291 cancer-free adult controls of European ancestry. Two loci achieved genome-wide significance: a locus in the GRM4 gene at 6p21.3 (encoding glutamate receptor metabotropic 4; rs1906953; P = 8.1 × 10⁻⁹) and a locus in the gene desert at 2p25.2 (rs7591996 and rs10208273; P = 1.0 × 10⁻⁸ and 2.9 × 10⁻⁷, respectively). These two loci warrant further exploration to uncover the biological mechanisms underlying susceptibility to osteosarcoma.

The metazoan genome is replicated in precise cell lineage-specific temporal order. However, the mechanism controlling this orchestrated process is poorly understood as no molecular mechanisms have been identified that actively regulate the firing sequence of genome replication. Here, we develop a mechanistic model of genome replication capable of predicting, with accuracy rivaling experimental repeats, observed empirical replication timing program in humans. In our model, replication is initiated in an uncoordinated (time-stochastic) manner at well-defined sites. The model contains, in addition to the choice of the genomic landmark that localizes initiation, only a single adjustable parameter of direct biological relevance: the number of replication forks. We find that DNase-hypersensitive sites are optimal and independent determinants of DNA replication initiation. We demonstrate that the DNA replication timing program in human cells is a robust emergent phenomenon that, by its very nature, does not require a regulatory mechanism determining a proper replication initiation firing sequence.