Expression of the p53-inducible antiproliferative gene BTG2 is suppressed in many cancers in the absence of inactivating gene mutations, suggesting alternative mechanisms of silencing. Using a shRNA screen targeting 43 histone lysine methyltransferases (KMTs), we show that SETD1A suppresses BTG2 expression through its induction of several BTG2-targeting miRNAs. This indirect but highly specific mechanism, by which a chromatin regulator that mediates transcriptional activating marks can lead to the downregulation of a critical effector gene, is shared with multiple genes in the p53 pathway. Through such miRNA-dependent effects, SETD1A regulates cell cycle progression in vitro and modulates tumorigenesis in mouse xenograft models. Together, these observations help explain the remarkably specific genetic consequences associated with alterations in generic chromatin modulators in cancer.

Acquired chromosomal DNA amplifications are features of many tumors. Although overexpression and stabilization of the histone H3 lysine 9/36 (H3K9/36) tri-demethylase KDM4A generates transient site-specific copy number gains (TSSGs), additional mechanisms directly controlling site-specific DNA copy gains are not well defined. In this study, we uncover a collection of H3K4-modifying chromatin regulators that function with H3K9 and H3K36 regulators to orchestrate TSSGs. Specifically, the H3K4 tri-demethylase KDM5A and specific COMPASS/KMT2 H3K4 methyltransferases modulate different TSSG loci through H3K4 methylation states and KDM4A recruitment. Furthermore, a distinct chromatin modifier network, MLL1-KDM4B-KDM5B, controls copy number regulation at a specific genomic locus in a KDM4A-independent manner. These pathways comprise an epigenetic addressing system for defining site-specific DNA rereplication and amplifications.

Chromatin-modifying enzymes play a fundamental role in regulating chromatin structure so that DNA replication is spatially and temporally coordinated. For example, the lysine demethylase 4A/Jumonji domain-containing 2A (KDM4A/JMJD2A) is tightly regulated during the cell cycle. Overexpression of JMJD2A leads to altered replication timing and faster S phase progression. In this study, we demonstrate that degradation of JMJD2A is regulated by the proteasome. JMJD2A turnover is coordinated through the SKP1-Cul1-F-box ubiquitin ligase complex that contains cullin 1 and the F-box and leucine-rich repeat protein 4 (FbxL4). This complex interacted with JMJD2A. Ubiquitin overexpression restored turnover and blocked the JMJD2A-dependent faster S phase progression in a cullin 1-dependent manner. Furthermore, increased ubiquitin levels decreased JMJD2A occupancy and BrdU incorporation at target sites. This study highlights a finely tuned mechanism for regulating histone demethylase levels and emphasizes the need to tightly regulate chromatin modifiers so that the cell cycle occurs properly.

The KDM4/JMJD2 family of histone demethylases is amplified in human cancers. However, little is known about their physiologic or tumorigenic roles. We have identified a conserved and unappreciated role for the JMJD2A/KDM4A H3K9/36 tridemethylase in cell cycle progression. We demonstrate that JMJD2A protein levels are regulated in a cell cycle-dependent manner and that JMJD2A overexpression increased chromatin accessibility, S phase progression, and altered replication timing of specific genomic loci. These phenotypes depended on JMJD2A enzymatic activity. Strikingly, depletion of the only C. elegans homolog, JMJD-2, slowed DNA replication and increased ATR/p53-dependent apoptosis. Importantly, overexpression of HP1γ antagonized JMJD2A-dependent progression through S phase, and depletion of HPL-2 rescued the DNA replication-related phenotypes in jmjd-2(-/-) animals. Our findings describe a highly conserved model whereby JMJD2A regulates DNA replication by antagonizing HP1γ and controlling chromatin accessibility.

DNA double-strand break (DSB) repair involves complex interactions between chromatin and repair proteins, including Tip60, a tumour suppressor. Tip60 is an acetyltransferase that acetylates both histones and ATM (ataxia telangiectasia mutated) kinase. Inactivation of Tip60 leads to defective DNA repair and increased cancer risk. However, how DNA damage activates the acetyltransferase activity of Tip60 is not known. Here, we show that direct interaction between the chromodomain of Tip60 and histone H3 trimethylated on lysine 9 (H3K9me3) at DSBs activates the acetyltransferase activity of Tip60. Depletion of intracellular H3K9me3 blocks activation of the acetyltransferase activity of Tip60, resulting in defective ATM activation and widespread defects in DSB repair. In addition, the ability of Tip60 to access H3K9me3 is dependent on the DNA damage-induced displacement of HP1beta (heterochromatin protein 1beta) from H3K9me3. Finally, we demonstrate that the Mre11-Rad50-Nbs1 (MRN) complex targets Tip60 to H3K9me3, and is required to activate the acetyltransferase activity of Tip60. These results reveal a new function for H3K9me3 in coordinating activation of Tip60-dependent DNA repair pathways, and imply that aberrant patterns of histone methylation may contribute to cancer by altering the efficiency of DSB repair.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

Despite continued widespread use, the genomic effects of cisplatin-based chemotherapy and implications for subsequent treatment are incompletely characterized. Here, we analyze whole exome sequencing of matched pre- and post-neoadjuvant cisplatin-based chemotherapy primary bladder tumor samples from 30 muscle-invasive bladder cancer patients. We observe no overall increase in tumor mutational burden post-chemotherapy, though a significant proportion of subclonal mutations are unique to the matched pre- or post-treatment tumor, suggesting chemotherapy-induced and/or spatial heterogeneity. We subsequently identify and validate a novel mutational signature in post-treatment tumors consistent with known characteristics of cisplatin damage and repair. We find that post-treatment tumor heterogeneity predicts worse overall survival, and further observe alterations in cell-cycle and immune checkpoint regulation genes in post-treatment tumors. These results provide insight into the clinical and genomic dynamics of tumor evolution with cisplatin-based chemotherapy, suggest mechanisms of clinical resistance, and inform development of clinically relevant biomarkers and trials of combination therapies.

Current genomics methods are designed to handle tens to thousands of samples but will need to scale to millions to match the pace of data and hypothesis generation in biomedical science. Here, we show that high efficiency at low cost can be achieved by leveraging general-purpose libraries for computing using graphics processing units (GPUs), such as PyTorch and TensorFlow. We demonstrate > 200-fold decreases in runtime and ~ 5-10-fold reductions in cost relative to CPUs. We anticipate that the accessibility of these libraries will lead to a widespread adoption of GPUs in computational genomics.

Chromatin modulation provides a key checkpoint for controlling cell cycle regulated gene networks. The replicative canonical histone genes are one such gene family under tight regulation during cell division. These genes are most highly expressed during S phase when histones are needed to chromatinize the new DNA template. While this fact has been known for a while, limited knowledge exists about the specific chromatin regulators controlling their temporal expression during cell cycle. Since histones and their associated mutations are emerging as major players in diseases such as cancer, identifying the chromatin factors modulating their expression is critical. The histone lysine tri-demethylase KDM4A is regulated over cell cycle and plays a direct role in DNA replication timing, site-specific rereplication, and DNA amplifications during S phase. Here, we establish an unappreciated role for the catalytically active KDM4A in directly regulating canonical replicative histone gene networks during cell cycle. Of interest, we further demonstrate that KDM4A interacts with proteins controlling histone expression and RNA processing (i.e., hnRNPUL1 and FUS/TLS). Together, this study provides a new function for KDM4A in modulating canonical histone gene expression.

Totipotency emerges in early embryogenesis, but its molecular underpinnings remain poorly characterized. In the present study, we employed DNA fiber analysis to investigate how pluripotent stem cells are reprogrammed into totipotent-like 2-cell-like cells (2CLCs). We show that totipotent cells of the early mouse embryo have slow DNA replication fork speed and that 2CLCs recapitulate this feature, suggesting that fork speed underlies the transition to a totipotent-like state. 2CLCs emerge concomitant with DNA replication and display changes in replication timing (RT), particularly during the early S-phase. RT changes occur prior to 2CLC emergence, suggesting that RT may predispose to gene expression changes and consequent reprogramming of cell fate. Slowing down replication fork speed experimentally induces 2CLCs. In vivo, slowing fork speed improves the reprogramming efficiency of somatic cell nuclear transfer. Our data suggest that fork speed regulates cellular plasticity and that remodeling of replication features leads to changes in cell fate and reprogramming.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 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 analyses3-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.

Alterations in DNA repair pathways are common in tumors and can result in characteristic mutational signatures; however, a specific mutational signature associated with somatic alterations in the nucleotide- excision repair (NER) pathway has not yet been identified. Here we examine the mutational processes operating in urothelial cancer, a tumor type in which the core NER gene ERCC2 is significantly mutated. Analysis of three independent urothelial tumor cohorts demonstrates a strong association between somatic ERCC2 mutations and the activity of a mutational signature characterized by a broad spectrum of base changes. In addition, we note an association between the activity of this signature and smoking that is independent of ERCC2 mutation status, providing genomic evidence of tobacco-related mutagenesis in urothelial cancer. Together, these analyses identify an NER-related mutational signature and highlight the related roles of DNA damage and subsequent DNA repair in shaping tumor mutational landscape.

Tumour expression of selected microRNAs (miRs) correlates with cisplatin efficacy in multiple cancers. We investigated the role of selected miRs in patients receiving cisplatin-based therapy for advanced urothelial carcinoma (UC).

The retinoblastoma (RB) family of proteins regulate transcription. These proteins lack intrinsic DNA-binding activity but are recruited to specific genomic locations through interactions with sequence-specific DNA-binding factors. The best-known target of RB protein (pRB) is the E2F transcription factor; however, many other chromatin-associated proteins have been described that may allow RB family members to act at additional sites. To gain a perspective on the scale of E2F-dependent and E2F-independent functions, we generated genome-wide binding profiles of RBF1 and dE2F proteins in Drosophila larvae. RBF1 and dE2F2 associate with a large number of binding sites at genes with diverse biological functions. In contrast, dE2F1 was detected at a smaller set of promoters, suggesting that it overrides repression by RBF1/dE2F2 at a specific subset of targets. Approximately 15% of RBF1-bound regions lacked consensus E2F-binding motifs. To test whether RBF1 action at these sites is E2F independent, we examined dDP mutant larvae that lack any functional dE2F/dDP heterodimers. As measured by chromatin immunoprecipitation-microarray analysis (ChIP-chip), ChIP-quantitative PCR (qPCR), and cell fractionation, the stable association of RBF1 with chromatin was eliminated in dDP mutants. This requirement for dDP was seen at classic E2F-regulated promoters and at promoters that lacked canonical E2F-binding sites. These results suggest that E2F/DP complexes are essential for all genomic targeting of RBF1.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Daytime napping is a common, heritable behavior, but its genetic basis and causal relationship with cardiometabolic health remain unclear. Here, we perform a genome-wide association study of self-reported daytime napping in the UK Biobank (n = 452,633) and identify 123 loci of which 61 replicate in the 23andMe research cohort (n = 541,333). Findings include missense variants in established drug targets for sleep disorders (HCRTR1, HCRTR2), genes with roles in arousal (TRPC6, PNOC), and genes suggesting an obesity-hypersomnolence pathway (PNOC, PATJ). Association signals are concordant with accelerometer-measured daytime inactivity duration and 33 loci colocalize with loci for other sleep phenotypes. Cluster analysis identifies three distinct clusters of nap-promoting mechanisms with heterogeneous associations with cardiometabolic outcomes. Mendelian randomization shows potential causal links between more frequent daytime napping and higher blood pressure and waist circumference.

Aberrant regulation of microRNA (miRNA) machinery components is associated with various human cancers, including papillary thyroid carcinoma (PTC), which is the most common type of thyroid cancer, and a higher prevalent female malignancy. The purpose of this study is to investigate racial differences in mRNA expression levels of four miRNA machinery components, Dicer, Drosha, DGCR8, and AGO2, and their correlations with clinicopathological characteristics. Forty PTC samples from female Asian Korean PTC patients were enrolled. Using qPCR, we examined mRNA expression levels of the components and next validated our results by comparison with results of female white American in the TCGA PTC project. Interestingly, mRNA expression levels of the selected factors were altered in the TCGA PTC samples. However, only Drosha showed a significantly lower expression level in Asian Korean PTC samples. Furthermore, the mRNA expression levels of the four components showed no association with clinicopathological characteristics in both groups. On the other hand, positive correlations were observed between altered mRNA expression levels of Dicer and Drosha and DGCR8 and Drosha in TCGA PTC samples. These findings collectively revealed that altered mRNA expression levels of miRNA machinery components might be responsible for racial differences in the carcinogenesis of PTC.

Anti-PD-1/PD-L1 agents have transformed the treatment landscape of advanced non-small cell lung cancer (NSCLC). To expand our understanding of the molecular features underlying response to checkpoint inhibitors in NSCLC, we describe here the first joint analysis of the Stand Up To Cancer-Mark Foundation cohort, a resource of whole exome and/or RNA sequencing from 393 patients with NSCLC treated with anti-PD-(L)1 therapy, along with matched clinical response annotation. We identify a number of associations between molecular features and outcome, including (1) favorable (for example, ATM altered) and unfavorable (for example, TERT amplified) genomic subgroups, (2) a prominent association between expression of inducible components of the immunoproteasome and response and (3) a dedifferentiated tumor-intrinsic subtype with enhanced response to checkpoint blockade. Taken together, results from this cohort demonstrate the complexity of biological determinants underlying immunotherapy outcomes and reinforce the discovery potential of integrative analysis within large, well-curated, cancer-specific cohorts.

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations1,2. Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models3-8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

We analyzed 921 adenocarcinomas of the esophagus, stomach, colon, and rectum to examine shared and distinguishing molecular characteristics of gastrointestinal tract adenocarcinomas (GIACs). Hypermutated tumors were distinct regardless of cancer type and comprised those enriched for insertions/deletions, representing microsatellite instability cases with epigenetic silencing of MLH1 in the context of CpG island methylator phenotype, plus tumors with elevated single-nucleotide variants associated with mutations in POLE. Tumors with chromosomal instability were diverse, with gastroesophageal adenocarcinomas harboring fragmented genomes associated with genomic doubling and distinct mutational signatures. We identified a group of tumors in the colon and rectum lacking hypermutation and aneuploidy termed genome stable and enriched in DNA hypermethylation and mutations in KRAS, SOX9, and PCBP1.