Center for Cancer Research

We present Model-based Analysis of ChIP-Seq data, MACS, which analyzes data generated by short read sequencers such as Solexa's Genome Analyzer. MACS empirically models the shift size of ChIP-Seq tags, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome, allowing for more robust predictions. MACS compares favorably to existing ChIP-Seq peak-finding algorithms, and is freely available.

The last decade has seen a sharp increase in the number of scientific publications describing physiological and pathological functions of extracellular vesicles (EVs), a collective term covering various subtypes of cell-released, membranous structures, called exosomes, microvesicles, microparticles, ectosomes, oncosomes, apoptotic bodies, and many other names. However, specific issues arise when working with these entities, whose size and amount often make them difficult to obtain as relatively pure preparations, and to characterize properly. The International Society for Extracellular Vesicles (ISEV) proposed Minimal Information for Studies of Extracellular Vesicles ("MISEV") guidelines for the field in 2014. We now update these "MISEV2014" guidelines based on evolution of the collective knowledge in the last four years. An important point to consider is that ascribing a specific function to EVs in general, or to subtypes of EVs, requires reporting of specific information beyond mere description of function in a crude, potentially contaminated, and heterogeneous preparation. For example, claims that exosomes are endowed with exquisite and specific activities remain difficult to support experimentally, given our still limited knowledge of their specific molecular machineries of biogenesis and release, as compared with other biophysically similar EVs. The MISEV2018 guidelines include tables and outlines of suggested protocols and steps to follow to document specific EV-associated functional activities. Finally, a checklist is provided with summaries of key points.

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human 'knockout' variants in protein-coding genes.

The genomes of cancers deficient in mismatch repair contain exceptionally high numbers of somatic mutations. In a proof-of-concept study, we previously showed that colorectal cancers with mismatch repair deficiency were sensitive to immune checkpoint blockade with antibodies to programmed death receptor-1 (PD-1). We have now expanded this study to evaluate the efficacy of PD-1 blockade in patients with advanced mismatch repair-deficient cancers across 12 different tumor types. Objective radiographic responses were observed in 53% of patients, and complete responses were achieved in 21% of patients. Responses were durable, with median progression-free survival and overall survival still not reached. Functional analysis in a responding patient demonstrated rapid in vivo expansion of neoantigen-specific T cell clones that were reactive to mutant neopeptides found in the tumor. These data support the hypothesis that the large proportion of mutant neoantigens in mismatch repair-deficient cancers make them sensitive to immune checkpoint blockade, regardless of the cancers' tissue of origin.

Small noncoding microRNAs (miRNAs) can contribute to cancer development and progression and are differentially expressed in normal tissues and cancers. From a large-scale miRnome analysis on 540 samples including lung, breast, stomach, prostate, colon, and pancreatic tumors, we identified a solid cancer miRNA signature composed by a large portion of overexpressed miRNAs. Among these miRNAs are some with well characterized cancer association, such as miR-17-5p, miR-20a, miR-21, miR-92, miR-106a, and miR-155. The predicted targets for the differentially expressed miRNAs are significantly enriched for protein-coding tumor suppressors and oncogenes (P < 0.0001). A number of the predicted targets, including the tumor suppressors RB1 (Retinoblastoma 1) and TGFBR2 (transforming growth factor, beta receptor II) genes were confirmed experimentally. Our results indicate that miRNAs are extensively involved in cancer pathogenesis of solid tumors and support their function as either dominant or recessive cancer genes.

The 2005 National Institutes of Health (NIH) Consensus Conference proposed new criteria for diagnosing and scoring the severity of chronic graft-versus-host disease (GVHD). The 2014 NIH consensus maintains the framework of the prior consensus with further refinement based on new evidence. Revisions have been made to address areas of controversy or confusion, such as the overlap chronic GVHD subcategory and the distinction between active disease and past tissue damage. Diagnostic criteria for involvement of mouth, eyes, genitalia, and lungs have been revised. Categories of chronic GVHD should be defined in ways that indicate prognosis, guide treatment, and define eligibility for clinical trials. Revisions have been made to focus attention on the causes of organ-specific abnormalities. Attribution of organ-specific abnormalities to chronic GVHD has been addressed. This paradigm shift provides greater specificity and more accurately measures the global burden of disease attributed to GVHD, and it will facilitate biomarker association studies.

Diffuse large B-cell lymphoma (DLBCL) can be divided into prognostically important subgroups with germinal center B-cell-like (GCB), activated B-cell-like (ABC), and type 3 gene expression profiles using a cDNA microarray. Tissue microarray (TMA) blocks were created from 152 cases of DLBCL, 142 of which had been successfully evaluated by cDNA microarray (75 GCB, 41 ABC, and 26 type 3). Sections were stained with antibodies to CD10, bcl-6, MUM1, FOXP1, cyclin D2, and bcl-2. Expression of bcl-6 (P <.001) or CD10 (P =.019) was associated with better overall survival (OS), whereas expression of MUM1 (P =.009) or cyclin D2 (P <.001) was associated with worse OS. Cases were subclassified using CD10, bcl-6, and MUM1 expression, and 64 cases (42%) were considered GCB and 88 cases (58%) non-GCB. The 5-year OS for the GCB group was 76% compared with only 34% for the non-GCB group (P <.001), which is similar to that reported using the cDNA microarray. Bcl-2 and cyclin D2 were adverse predictors in the non-GCB group. In multivariate analysis, a high International Prognostic Index score (3-5) and the non-GCB phenotype were independent adverse predictors (P <.0001). In summary, immunostains can be used to determine the GCB and non-GCB subtypes of DLBCL and predict survival similar to the cDNA microarray.

BACKGROUND: The survival of patients with diffuse large-B-cell lymphoma after chemotherapy is influenced by molecular features of the tumors. We used the gene-expression profiles of these lymphomas to develop a molecular predictor of survival. METHODS: Biopsy samples of diffuse large-B-cell lymphoma from 240 patients were examined for gene expression with the use of DNA microarrays and analyzed for genomic abnormalities. Subgroups with distinctive gene-expression profiles were defined on the basis of hierarchical clustering. A molecular predictor of risk was constructed with the use of genes with expression patterns that were associated with survival in a preliminary group of 160 patients and was then tested in a validation group of 80 patients. The accuracy of this predictor was compared with that of the international prognostic index. RESULTS: Three gene-expression subgroups--germinal-center B-cell-like, activated B-cell-like, and type 3 diffuse large-B-cell lymphoma--were identified. Two common oncogenic events in diffuse large-B-cell lymphoma, bcl-2 translocation and c-rel amplification, were detected only in the germinal-center B-cell-like subgroup. Patients in this subgroup had the highest five-year survival rate. To identify other molecular determinants of outcome, we searched for individual genes with expression patterns that correlated with survival in the preliminary group of patients. Most of these genes fell within four gene-expression signatures characteristic of germinal-center B cells, proliferating cells, reactive stromal and immune cells in the lymph node, or major-histocompatibility-complex class II complex. We used 17 genes to construct a predictor of overall survival after chemotherapy. This gene-based predictor and the international prognostic index were independent prognostic indicators. CONCLUSIONS: DNA microarrays can be used to formulate a molecular predictor of survival after chemotherapy for diffuse large-B-cell lymphoma.

Familial cancer syndromes have helped to define the role of tumor suppressor genes in the development of cancer. The dominantly inherited Li-Fraumeni syndrome (LFS) is of particular interest because of the diversity of childhood and adult tumors that occur in affected individuals. The rarity and high mortality of LFS precluded formal linkage analysis. The alternative approach was to select the most plausible candidate gene. The tumor suppressor gene, p53, was studied because of previous indications that this gene is inactivated in the sporadic (nonfamilial) forms of most cancers that are associated with LFS. Germ line p53 mutations have been detected in all five LFS families analyzed. These mutations do not produce amounts of mutant p53 protein expected to exert a trans-dominant loss of function effect on wild-type p53 protein. The frequency of germ line p53 mutations can now be examined in additional families with LFS, and in other cancer patients and families with clinical features that might be attributed to the mutation.

Abstract Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature 1 . Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium 2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses 3–15 , enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated—but distinct—DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Genetic changes underlying clear cell renal cell carcinoma (ccRCC) include alterations in genes controlling cellular oxygen sensing (for example, VHL) and the maintenance of chromatin states (for example, PBRM1). We surveyed more than 400 tumours using different genomic platforms and identified 19 significantly mutated genes. The PI(3)K/AKT pathway was recurrently mutated, suggesting this pathway as a potential therapeutic target. Widespread DNA hypomethylation was associated with mutation of the H3K36 methyltransferase SETD2, and integrative analysis suggested that mutations involving the SWI/SNF chromatin remodelling complex (PBRM1, ARID1A, SMARCA4) could have far-reaching effects on other pathways. Aggressive cancers demonstrated evidence of a metabolic shift, involving downregulation of genes involved in the TCA cycle, decreased AMPK and PTEN protein levels, upregulation of the pentose phosphate pathway and the glutamine transporter genes, increased acetyl-CoA carboxylase protein, and altered promoter methylation of miR-21 (also known as MIR21) and GRB10. Remodelling cellular metabolism thus constitutes a recurrent pattern in ccRCC that correlates with tumour stage and severity and offers new views on the opportunities for disease treatment. The Cancer Genome Atlas Research Network reports an integrative analysis of more than 400 samples of clear cell renal cell carcinoma based on genomic, DNA methylation, RNA and proteomic characterisation; frequent mutations were identified in the PI(3)K/AKT pathway, suggesting this pathway might be a potential therapeutic target, among the findings is also a demonstration of metabolic remodelling which correlates with tumour stage and severity. The Cancer Genome Atlas consortium reports an integrative analysis of more than 400 samples of clear cell renal carcinoma on the basis of genomic, DNA methylation, RNA and proteomic characterization. The data reveal frequent mutations in the PI(3)K/AKT pathway, suggesting that this pathway might be a potential therapeutic target, in addition to an array of epigenetic alterations that are linked to specific mutations in chromatin-associated proteins. One notable finding is the presence of a metabolic shift in aggressive cancers, correlating with tumour stage and severity.

Chronic lymphocytic leukemia (CLL) is the most common human leukemia and is characterized by predominantly nondividing malignant B cells overexpressing the antiapoptotic B cell lymphoma 2 (Bcl2) protein. miR-15a and miR-16-1 are deleted or down-regulated in the majority of CLLs. Here, we demonstrate that miR-15a and miR-16-1 expression is inversely correlated to Bcl2 expression in CLL and that both microRNAs negatively regulate Bcl2 at a posttranscriptional level. BCL2 repression by these microRNAs induces apoptopsis in a leukemic cell line model. Therefore, miR-15 and miR-16 are natural antisense Bcl2 interactors that could be used for therapy of Bcl2-overexpressing tumors.

The objective of this review article is to summarize current knowledge of blood flow and perfusion-related parameters, which usually go hand in hand and in turn define the cellular metabolic microenvironment of human malignancies. A compilation of available data from the literature on blood flow, oxygen and nutrient supply, and tissue oxygen and pH distribution in human tumors is presented. Whenever possible, data obtained for human tumors are compared with the respective parameters in normal tissues, isotransplanted or spontaneous rodent tumors, and xenografted human tumors. Although data on human tumors in situ are scarce and there may be significant errors associated with the techniques used for measurements, experimental evidence is provided for the existence of a compromised and anisotropic blood supply to many tumors. As a result, O2-depleted areas develop in human malignancies which coincide with nutrient and energy deprivation and with a hostile metabolic microenvironment (e.g., existence of severe tissue acidosis). Significant variations in these relevant parameters must be expected between different locations within the same tumor, at the same location at different times, and between individual tumors of the same grading and staging. Furthermore, this synopsis will attempt to identify relevant pathophysiological parameters and other related areas future research of which might be most beneficial for designing individually tailored treatment protocols with the goal of predicting the acute and/or long-term response of tumors to therapy.

Hybridization of lambda-mRNAs to excess liver DNA yielded results compatible with gene reiteration frequencies of three or less. Purified mRNA from tumors producing structurally different lambda chains were used in competition hybridization experiments. An unlabeled lambda-mRNA competed with another, labelled mRNA to the same extent as homologous unlabelled lambda-mRNA. Mouse DNA was digested with Eco R-I restriction endonuclease and fractionated by gel electrophoresis. A DNA fragment carrying the V lambda-gene(s) was indentified in this digest. This fragment hybridized with lambda-mRNAs coding for two different lambdaV regions equally well. These results indicate that base sequence homology among lambda-mRNAs is so high that any lambda-mRNA should cross-hybridize with all or most of germ line V lambda genes. From amino acid sequence data, it is argued that there are probably more than 25 different lambdaV regions. Hence it is concluded that the number of germ line genes is too small to account for the diversity of lambda chains.

Abstract Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale 1–3 . Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4–5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter 4 ; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation 5,6 ; analyses timings and patterns of tumour evolution 7 ; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity 8,9 ; and evaluates a range of more-specialized features of cancer genomes 8,10–18 .

We describe here the development and implementation of a pilot-scale, in vitro, anticancer drug screen utilizing a panel of 60 human tumor cell lines organized into subpanels representing leukemia, melanoma, and cancers of the lung, colon, kidney, ovary, and central nervous system. The ultimate goal of this disease-oriented screen is to facilitate the discovery of new compounds with potential cell line-specific and/or subpanel-specific antitumor activity. In the current screening protocol, each cell line is inoculated onto microtiter plates, then preincubated for 24-28 hours. Subsequently, test agents are added in five 10-fold dilutions and the culture is incubated for an additional 48 hours. For each test agent, a dose-response profile is generated. End-point determinations of the cell viability or cell growth are performed by in situ fixation of cells, followed by staining with a protein-binding dye, sulforhodamine B (SRB). The SRB binds to the basic amino acids of cellular macromolecules; the solubilized stain is measured spectrophotometrically to determine relative cell growth or viability in treated and untreated cells. Following the pilot screening studies, a screening rate of 400 compounds per week has been consistently achieved.

For decades, studies of endocrine-disrupting chemicals (EDCs) have challenged traditional concepts in toxicology, in particular the dogma of "the dose makes the poison," because EDCs can have effects at low doses that are not predicted by effects at higher doses. Here, we review two major concepts in EDC studies: low dose and nonmonotonicity. Low-dose effects were defined by the National Toxicology Program as those that occur in the range of human exposures or effects observed at doses below those used for traditional toxicological studies. We review the mechanistic data for low-dose effects and use a weight-of-evidence approach to analyze five examples from the EDC literature. Additionally, we explore nonmonotonic dose-response curves, defined as a nonlinear relationship between dose and effect where the slope of the curve changes sign somewhere within the range of doses examined. We provide a detailed discussion of the mechanisms responsible for generating these phenomena, plus hundreds of examples from the cell culture, animal, and epidemiology literature. We illustrate that nonmonotonic responses and low-dose effects are remarkably common in studies of natural hormones and EDCs. Whether low doses of EDCs influence certain human disorders is no longer conjecture, because epidemiological studies show that environmental exposures to EDCs are associated with human diseases and disabilities. We conclude that when nonmonotonic dose-response curves occur, the effects of low doses cannot be predicted by the effects observed at high doses. Thus, fundamental changes in chemical testing and safety determination are needed to protect human health.

Urothelial carcinoma of the bladder is a common malignancy that causes approximately 150,000 deaths per year worldwide. So far, no molecularly targeted agents have been approved for treatment of the disease. As part of The Cancer Genome Atlas project, we report here an integrated analysis of 131 urothelial carcinomas to provide a comprehensive landscape of molecular alterations. There were statistically significant recurrent mutations in 32 genes, including multiple genes involved in cell-cycle regulation, chromatin regulation, and kinase signalling pathways, as well as 9 genes not previously reported as significantly mutated in any cancer. RNA sequencing revealed four expression subtypes, two of which (papillary-like and basal/squamous-like) were also evident in microRNA sequencing and protein data. Whole-genome and RNA sequencing identified recurrent in-frame activating FGFR3–TACC3 fusions and expression or integration of several viruses (including HPV16) that are associated with gene inactivation. Our analyses identified potential therapeutic targets in 69% of the tumours, including 42% with targets in the phosphatidylinositol-3-OH kinase/AKT/mTOR pathway and 45% with targets (including ERBB2) in the RTK/MAPK pathway. Chromatin regulatory genes were more frequently mutated in urothelial carcinoma than in any other common cancer studied so far, indicating the future possibility of targeted therapy for chromatin abnormalities. This paper reports integrative molecular analyses of urothelial bladder carcinoma at the DNA, RNA, and protein levels performed as part of The Cancer Genome Atlas project; recurrent mutations were found in 32 genes, including those involved in cell-cycle regulation, chromatin regulation and kinase signalling pathways; chromatin regulatory genes were more frequently mutated in urothelial carcinoma than in any other common cancer studied so far. This study of 131 high-grade muscle-invasive urothelial bladder carcinomas, part of The Cancer Genome Atlas (TCGA) project, reports recurrent mutations in 32 genes, including those involved in cell-cycle regulation, chromatin regulation and kinase signalling pathways. Chromatin regulatory genes were more frequently mutated in urothelial carcinoma than in any common cancer studied to date. Recurrent in-frame activating FGFR3–TACC3 fusions and expression or integration of viruses associated with gene inactivation are also identified. Importantly, potential therapeutic targets are identified in 69% of the tumours.

A gene discovered by positional cloning has been identified as the von Hippel-Lindau (VHL) disease tumor suppressor gene. A restriction fragment encompassing the gene showed rearrangements in 28 of 221 VHL kindreds. Eighteen of these rearrangements were due to deletions in the candidate gene, including three large nonoverlapping deletions. Intragenic mutations were detected in cell lines derived from VHL patients and from sporadic renal cell carcinomas. The VHL gene is evolutionarily conserved and encodes two widely expressed transcripts of approximately 6 and 6.5 kilobases. The partial sequence of the inferred gene product shows no homology to other proteins, except for an acidic repeat domain found in the procyclic surface membrane glycoprotein of Trypanosoma brucei.

We recently showed that the mammalian genome encodes >1,000 large intergenic noncoding (linc)RNAs that are clearly conserved across mammals and, thus, functional. Gene expression patterns have implicated these lincRNAs in diverse biological processes, including cell-cycle regulation, immune surveillance, and embryonic stem cell pluripotency. However, the mechanism by which these lincRNAs function is unknown. Here, we expand the catalog of human lincRNAs to approximately 3,300 by analyzing chromatin-state maps of various human cell types. Inspired by the observation that the well-characterized lincRNA HOTAIR binds the polycomb repressive complex (PRC)2, we tested whether many lincRNAs are physically associated with PRC2. Remarkably, we observe that approximately 20% of lincRNAs expressed in various cell types are bound by PRC2, and that additional lincRNAs are bound by other chromatin-modifying complexes. Also, we show that siRNA-mediated depletion of certain lincRNAs associated with PRC2 leads to changes in gene expression, and that the up-regulated genes are enriched for those normally silenced by PRC2. We propose a model in which some lincRNAs guide chromatin-modifying complexes to specific genomic loci to regulate gene expression.