Search

We perform an immunogenomics analysis utilizing whole-transcriptome sequencing of 657 pediatric extracranial solid cancer samples representing 14 diagnoses, and additionally utilize transcriptomes of 131 pediatric cancer cell lines and 147 normal tissue samples for comparison. We describe patterns of infiltrating immune cells, T cell receptor (TCR) clonal expansion, and translationally relevant immune checkpoints. We find that tumor-infiltrating lymphocytes and TCR counts vary widely across cancer types and within each diagnosis, and notably are significantly predictive of survival in osteosarcoma patients. We identify potential cancer-specific immunotherapeutic targets for adoptive cell therapies including cell-surface proteins, tumor germline antigens, and lineage-specific transcription factors. Using an orthogonal immunopeptidomics approach, we find several potential immunotherapeutic targets in osteosarcoma and Ewing sarcoma and validated PRAME as a bona fide multi-pediatric cancer target. Importantly, this work provides a critical framework for immune targeting of extracranial solid tumors using parallel immuno-transcriptomic and -peptidomic approaches., Competing Interests: Declaration of interests A.S.B. has advisory board relationships with Bayer, EMD Serono, and Deciphera. R.J.O. receives research support from and consults for Lentigen, a Miltenyi Biotec Company, and also consults for Umoja Biopharma. Other authors declare no competing interests., (Published by Elsevier Inc.)

The RecQ DNA helicase WRN is a synthetic lethal target for cancer cells with microsatellite instability (MSI), a form of genetic hypermutability that arises from impaired mismatch repair 1-4 . Depletion of WRN induces widespread DNA double-strand breaks in MSI cells, leading to cell cycle arrest and/or apoptosis. However, the mechanism by which WRN protects MSI-associated cancers from double-strand breaks remains unclear. Here we show that TA-dinucleotide repeats are highly unstable in MSI cells and undergo large-scale expansions, distinct from previously described insertion or deletion mutations of a few nucleotides 5 . Expanded TA repeats form non-B DNA secondary structures that stall replication forks, activate the ATR checkpoint kinase, and require unwinding by the WRN helicase. In the absence of WRN, the expanded TA-dinucleotide repeats are susceptible to cleavage by the MUS81 nuclease, leading to massive chromosome shattering. These findings identify a distinct biomarker that underlies the synthetic lethal dependence on WRN, and support the development of therapeutic agents that target WRN for MSI-associated cancers.

CellMiner (http://discover.nci.nih.gov/cellminer) and CellMinerCDB (https://discover.nci.nih.gov/cellminercdb/) are web-based applications for mining publicly available genomic, molecular, and pharmacologic datasets of human cancer cell lines including the NCI-60, Cancer Cell Line Encyclopedia, Genomics of Drug Sensitivity in Cancer, Cancer Therapeutics Response Portal, NCI/DTP small cell lung cancer, and NCI Almanac cell line sets. Here, we introduce our RNA sequencing (RNA-seq) data for the NCI-60 and their access and integration with the other databases. Correlation to transcript microarray expression levels for identical genes and identical cell lines across CellMinerCDB demonstrates the high quality of these new RNA-seq data. We provide composite and isoform transcript expression data and demonstrate diversity in isoform composition for individual cancer- and pharmacologically relevant genes, including HRAS, PTEN, EGFR, RAD51, ALKBH2, BRCA1, ERBB2, TP53, FGFR2, and CTNND1. We reveal cell-specific differences in the overall levels of isoforms and show their linkage to expression of RNA processing and splicing genes as well as resultant alterations in cancer and pharmacologic gene sets. Gene-drug pairings linked by pathways or functions show specific correlations to isoforms compared with composite gene expression, including ALKBH2-benzaldehyde, AKT3-vandetanib, BCR-imatinib, CDK1 and 20-palbociclib, CASP1-imexon, and FGFR3-pazopanib. Loss of MUC1 20 amino acid variable number tandem repeats, which is used to elicit immune response, and the presence of the androgen receptor AR-V4 and -V7 isoforms in all NCI-60 tissue of origin types demonstrate translational relevance. In summary, we introduce RNA-seq data to our CellMiner and CellMinerCDB web applications, allowing their exploration for both research and translational purposes. SIGNIFICANCE: The current study provides RNA sequencing data for the NCI-60 cell lines made accessible through both CellMiner and CellMinerCDB and is an important pharmacogenomics resource for the field., (©2019 American Association for Cancer Research.)

While synthetic lethality (SL) holds promise in developing effective cancer therapies, SL candidates found via experimental screens often have limited translational value. Here we present a data-driven approach, ISLE (identification of clinically relevant synthetic lethality), that mines TCGA cohort to identify the most likely clinically relevant SL interactions (cSLi) from a given candidate set of lab-screened SLi. We first validate ISLE via a benchmark of large-scale drug response screens and by predicting drug efficacy in mouse xenograft models. We then experimentally test a select set of predicted cSLi via new screening experiments, validating their predicted context-specific sensitivity in hypoxic vs normoxic conditions and demonstrating cSLi's utility in predicting synergistic drug combinations. We show that cSLi can successfully predict patients' drug treatment response and provide patient stratification signatures. ISLE thus complements existing actionable mutation-based methods for precision cancer therapy, offering an opportunity to expand its scope to the whole genome.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

A unique resource for systems pharmacology and genomic studies is the NCI-60 cancer cell line panel, which provides data for the largest publicly available library of compounds with cytotoxic activity (∼21,000 compounds), including 108 FDA-approved and 70 clinical trial drugs as well as genomic data, including whole-exome sequencing, gene and miRNA transcripts, DNA copy number, and protein levels. Here, we provide the first readily usable genome-wide DNA methylation database for the NCI-60, including 485,577 probes from the Infinium HumanMethylation450k BeadChip array, which yielded DNA methylation signatures for 17,559 genes integrated into our open access CellMiner version 2.0 (https://discover.nci.nih.gov/cellminer). Among new insights, transcript versus DNA methylation correlations revealed the epithelial/mesenchymal gene functional category as being influenced most heavily by methylation. DNA methylation and copy number integration with transcript levels yielded an assessment of their relative influence for 15,798 genes, including tumor suppressor, mitochondrial, and mismatch repair genes. Four forms of molecular data were combined, providing rationale for microsatellite instability for 8 of the 9 cell lines in which it occurred. Individual cell line analyses showed global methylome patterns with overall methylation levels ranging from 17% to 84%. A six-gene model, including PARP1, EP300, KDM5C, SMARCB1, and UHRF1 matched this pattern. In addition, promoter methylation of two translationally relevant genes, Schlafen 11 (SLFN11) and methylguanine methyltransferase (MGMT), served as indicators of therapeutic resistance or susceptibility, respectively. Overall, our database provides a resource of pharmacologic data that can reinforce known therapeutic strategies and identify novel drugs and drug targets across multiple cancer types. Cancer Res; 77(3); 601-12. ©2016 AACR., (©2016 American Association for Cancer Research.)

Background: 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., Results: We implemented the R package, caOmicsV, to provide methods under R environment to visualize multidimensional cancer genomic data in two layouts: matrix layout and combined biological network and circular layout. Both layouts support to display sample information, gene expression (e.g., RNA and miRNA), DNA methylation, DNA copy number variations, and summarized data. A set of supplemental functions are included in the caOmicsV package to help users in generation of plot data sets from multiple genomic datasets with given gene names and sample names. Default plot methods for both layouts for easy use are also implemented., Conclusion: caOmicsV package provides an easy and flexible way to visualize integrated multidimensional cancer genomic data under R environment.

Most p53 mutations in human cancers are missense mutations resulting in a full-length mutant p53 protein. Besides losing tumor suppressor activity, some hotspot p53 mutants gain oncogenic functions. This effect is mediated in part, through gene expression changes due to inhibition of p63 and p73 by mutant p53 at their target gene promoters. Here, we report that the tumor suppressor microRNA let-7i is downregulated by mutant p53 in multiple cell lines expressing endogenous mutant p53. In breast cancer patients, significantly decreased let-7i levels were associated with missense mutations in p53. Chromatin immunoprecipitation and promoter luciferase assays established let-7i as a transcriptional target of mutant p53 through p63. Introduction of let-7i to mutant p53 cells significantly inhibited migration, invasion and metastasis by repressing a network of oncogenes including E2F5, LIN28B, MYC and NRAS. Our findings demonstrate that repression of let-7i expression by mutant p53 has a key role in enhancing migration, invasion and metastasis.

Genome-wide association studies (GWAS) have mapped risk alleles for at least 10 distinct cancers to a small region of 63 000 bp on chromosome 5p15.33. This region harbors the TERT and CLPTM1L genes; the former encodes the catalytic subunit of telomerase reverse transcriptase and the latter may play a role in apoptosis. To investigate further the genetic architecture of common susceptibility alleles in this region, we conducted an agnostic subset-based meta-analysis (association analysis based on subsets) across six distinct cancers in 34 248 cases and 45 036 controls. Based on sequential conditional analysis, we identified as many as six independent risk loci marked by common single-nucleotide polymorphisms: five in the TERT gene (Region 1: rs7726159, P = 2.10 × 10(-39); Region 3: rs2853677, P = 3.30 × 10(-36) and PConditional = 2.36 × 10(-8); Region 4: rs2736098, P = 3.87 × 10(-12) and PConditional = 5.19 × 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 × 10(-6); and Region 6: rs10069690, P = 7.49 × 10(-15) and PConditional = 5.35 × 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 × 10(-18) and PConditional = 7.06 × 10(-16)). Between three and five cancers mapped to each independent locus with both risk-enhancing and protective effects. Allele-specific effects on DNA methylation were seen for a subset of risk loci, indicating that methylation and subsequent effects on gene expression may contribute to the biology of risk variants on 5p15.33. Our results provide strong support for extensive pleiotropy across this region of 5p15.33, to an extent not previously observed in other cancer susceptibility loci., (Published by Oxford University Press 2014. This work is written by (a) US Government employee(s) and is in the public domain in the US.)

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., (Copyright © 2014. Published by Elsevier Inc.)

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 NCI-60 cell lines are the most frequently studied human tumor cell lines in cancer research. This panel has generated the most extensive cancer pharmacology database worldwide. In addition, these cell lines have been intensely investigated, providing a unique platform for hypothesis-driven research focused on enhancing our understanding of tumor biology. Here, we report a comprehensive analysis of coding variants in the NCI-60 panel of cell lines identified by whole exome sequencing, providing a list of possible cancer specific variants for the community. Furthermore, we identify pharmacogenomic correlations between specific variants in genes such as TP53, BRAF, ERBBs, and ATAD5 and anticancer agents such as nutlin, vemurafenib, erlotinib, and bleomycin showing one of many ways the data could be used to validate and generate novel hypotheses for further investigation. As new cancer genes are identified through large-scale sequencing studies, the data presented here for the NCI-60 will be an invaluable resource for identifying cell lines with mutations in such genes for hypothesis-driven research. To enhance the utility of the data for the greater research community, the genomic variants are freely available in different formats and from multiple sources including the CellMiner and Ingenuity websites., (©2013 AACR.)

CHEK2 encodes a serine/threonine kinase (Chk2) activated by ATM in response to DNA double-strand breaks. On the one hand, CHEK2 has been described as a tumor suppressor with proapoptotic, cell-cycle checkpoint and mitotic functions. On the other hand, Chk2 is also commonly activated (phosphorylated at T68) in cancers and precancerous lesions. Here, we report an extensive characterization of CHEK2 across the panel of 60 established cancer cell lines from the NCI Anticancer Screen (the NCI-60) using genomic and proteomic analyses, including exon-specific mRNA expression, DNA copy-number variation (CNV) by aCGH, exome sequencing, as well as western blot analyses for total and activated (pT68-Chk2) Chk2. We show that the high heterogeneity of Chk2 levels in cancer cells is primarily due to its inactivation (owing to low gene expression, alternative splicing, point mutations, copy-number alterations and premature truncation) or reduction of protein levels. Moreover, we observe that a significant percentage of cancer cells (12% of the NCI-60 and HeLa cells) show high endogenous Chk2 activation, which is always associated with p53 inactivation, and which is accompanied by downregulation of the Fanconi anemia and homologous recombination pathways. We also report the presence of activated Chk2 (pT68-Chk2) along with histone γ-H2AX in centrosomes.

The purpose of this study was to develop a method of classifying cancers to specific diagnostic categories based on their gene expression signatures using artificial neural networks (ANNs). We trained the ANNs using the small, round blue-cell tumors (SRBCTs) as a model. These cancers belong to four distinct diagnostic categories and often present diagnostic dilemmas in clinical practice. The ANNs correctly classified all samples and identified the genes most relevant to the classification. Expression of several of these genes has been reported in SRBCTs, but most have not been associated with these cancers. To test the ability of the trained ANN models to recognize SRBCTs, we analyzed additional blinded samples that were not previously used for the training procedure, and correctly classified them in all cases. This study demonstrates the potential applications of these methods for tumor diagnosis and the identification of candidate targets for therapy.

ASC-2 was recently discovered as a cancer-amplified transcription coactivator molecule of nuclear receptors, which interacts with multifunctional transcription integrators steroid receptor coactivator-1 (SRC-1) and CREB-binding protein (CBP)/p300. Herein, we report the identification of three mitogenic transcription factors as novel target molecules of ASC-2. First, the C-terminal transactivation domain of serum response factor (SRF) was identified among a series of ASC-2-interacting proteins from the yeast two-hybrid screening. Second, ASC-2 specifically interacted with the activating protein-1 (AP-1) components c-Jun and c-Fos as well as the nuclear factor-kappaB (NFkappaB) components p50 and p65, as demonstrated by the glutathione S-transferase pull-down assays as well as the yeast two-hybrid tests. In cotransfection of mammalian cells, ASC-2 potentiated transactivations by SRF, AP-1, and NFkappaB in a dose-dependent manner, either alone or in conjunction with SRC-1 and p300. In addition, ASC-2 efficiently relieved the previously described transrepression between nuclear receptors and either AP-1 or NFkappaB. Overall, these results suggest that the nuclear receptor coactivator ASC-2 also mediates transactivations by SRF, AP-1, and NFkappaB, which may contribute to the putative, ASC-2-mediated tumorigenesis.

Currently there are over 1,000,000 human expressed sequence tag (EST) sequences available on the public database, representing perhaps 50-90% of all human genes. The cDNA microarray technique is a recently developed tool that exploits this wealth of information for the analysis of gene expression. In this method, DNA probes representing cDNA clones are arrayed onto a glass slide and interrogated with fluorescently labeled cDNA targets. The power of the technology is the ability to perform a genome-wide expression profile of thousands of genes in one experiment. In our review we describe the principles of the microarray technology as applied to cancer research, summarize the literature on its use so far, and speculate on the future application of this powerful technique.

We used a combination of sequence analysis and exon trapping in an effort to determine the complete transcript map for a cosmid (6E5) derived from 12q13.3, a region of DNA sequence amplification in human cancers. This cosmid, previously known to contain three genes (CDK4, SAS, and OS9), was sequenced, and that information was used for computer-assisted analysis. In addition, 6E5 was subjected to both internal and 3'-terminal exon-trapping protocols, and the results of these studies were used to guide cDNA cloning experiments. These studies demonstrate that this cosmid is derived from a remarkably gene-dense region and add two new transcripts (KIAA0167 and 6E5.2) to the list of sequences that are expressed in tumors bearing amplification of this region.

We recently identified a novel tumor-suppressor gene, DPC4, at chromosome 18q21.1 and found that both alleles of DPC4 were inactivated in nearly one-half of the pancreatic carcinomas. Here, we analyzed 338 tumors, originating from 12 distinct anatomic sites, for alterations in the DPC4 gene. Sixty-four specimens were selected for the presence of the allelic loss of 18q and were further analyzed for DPC4 sequence alterations. An alteration of the DPC4 gene sequence was identified in one of eight breast carcinomas and one of eight ovarian carcinomas. These results indicate that whereas DPC4 inactivation is prevalent in pancreatic carcinoma (48%), it is distinctly uncommon (< 10%) in the other tumor types examined. The tissue restriction of alterations in DPC4, as in many other tumor-suppressor genes, emphasizes the complexity of rate-limiting checkpoints in human tumorigenesis.

It has been increasingly recognized that homogeneously staining regions (hsr) in human cancers may be complex structures composed of large amplified DNA domains containing multiple genes. It is therefore important to devise strategies for the rapid isolation of cDNAs expressed from these structures. Using a procedure we term microdissection mediated cDNA capture, we recovered hsr specific cDNAs from two different human tumors. The glioblastoma cell line TX3868 and the human sarcoma cell line OsA-CL carry hsrs containing amplified sequences from chromosome 12q13-15. We recovered 17 hsr specific cDNAs following microdissection of these hsrs which had been previously hybridized in situ with linkered cDNA. Northern blot analysis with these cDNAs revealed hybridization to distinct transcripts in OsA-CL RNA and TX3868 RNA. None of the OsA-CL cDNA clones showed cross hybridization with the TX3868 cDNAs suggesting that despite their coincident band localization on 12q, the OsA-CL and TX3868 amplification units do not completely overlap. These results significantly increase the number of amplified genes assigned to the 12q13-15 amplicon illustrating both the complexity of hsrs derived from this region and the utility of microdissection mediated cDNA capture to gain rapid access to cDNAs transcribed from amplified genes.

Chromosome-specific cosmid libraries are an extremely useful resource for positional cloning projects. Once a particular region of interest has been identified, it would be of value to have an approach for isolating chromosome band-specific cosmids that could be assembled into a sublibrary for rapid screening. We constructed a region-specific sublibrary of 700 cosmids by screening a chromosome 12-specific cosmid library with a complex probe generated by degenerate oligonucleotide-primed PCR of a microdissected homogeneously staining region containing sequences amplified from chromosome 12q13-q15. Based on fluorescence in situ hybridization, approximately 60% of the cosmids in the sublibrary were derived from the microdissected region. To demonstrate further the utility of this sublibrary, a 150-kb contig containing the SAS and CDK4 genes was constructed, as well as several additional contigs between CDK4 and MDM2. This study demonstrates the possibility of utilizing probes generated by microdissection for assembling band-specific sublibraries that are amenable to rapid screening with multiple markers.

Despite extensive data linking mutations in the p53 gene to human tumorigenesis, little is known about the cellular regulators and mediators of p53 function. MDM2 is a strong candidate for one such cellular protein; the MDM2 gene was originally identified by virtue of its amplification in a spontaneously transformed derivative of mouse BALB/c cells and the MDM2 protein subsequently shown to bind to p53 in rat cells transfected with p53 genes. To determine whether MDM2 plays a role in human cancer, we have cloned the human MDM2 gene. Here we show that recombinant-derived human MDM2 protein binds human p53 in vitro, and we use MDM2 clones to localize the human MDM2 gene to chromosome 12q13-14. Because this chromosomal position appears to be altered in many sarcomas, we looked for changes in human MDM2 in such cancers. The gene was amplified in over a third of 47 sarcomas, including common bone and soft tissue forms. These results are consistent with the hypothesis that MDM2 binds to p53, and that amplification of MDM2 in sarcomas leads to escape from p53-regulated growth control. This mechanism of tumorigenesis parallels that for virally-induced tumours, in which viral oncogene products bind to and functionally inactivate p53.

The amplification units in human tumors containing amplified myc genes were examined. The amplification unit in all cases consisted of a large genomic region coamplified with the coding region of the myc genes themselves. In eight independent neuroblastomas containing N-myc amplifications, the amplification unit was estimated to be 290 to 430 kilobases. This amplification unit was highly conserved among the different neuroblastomas, with some neuroblastomas containing almost identical units. In contrast, five tumor cell lines containing c-myc amplifications exhibited amplification units that were more variable in size (90 to 300 kilobases) and sequence content; at least three different patterns of c-myc amplification units could be discerned.