Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

Break Check: Transcription-Driven Topoisomerase II Collisions near Chromatin Loop Anchors Are Hotspots for DNA Damage and Translocations.

In this issue of Molecular Cell, Gothe, Maman et al. and Canela et al. demonstrate that type II topoisomerase (TOP2B/TOP2A)-mediated DNA breaks, and the oncogenic translocations they generate, are dependent on transcription and governed by chromatin architecture.

Interactome Rewiring Following Pharmacological Targeting of BET Bromodomains.

Targeting bromodomains (BRDs) of the bromo-and-extra-terminal (BET) family offers opportunities for therapeutic intervention in cancer and other diseases. Here, we profile the interactomes of BRD2, BRD3, BRD4, and BRDT following treatment with the pan-BET BRD inhibitor JQ1, revealing broad rewiring of the interaction landscape, with three distinct classes of behavior for the 603 unique interactors identified. A group of proteins associate in a JQ1-sensitive manner with BET BRDs through canonical and new binding modes, while two classes of extra-terminal (ET)-domain binding motifs mediate acetylation-independent interactions. Last, we identify an unexpected increase in several interactions following JQ1 treatment that define negative functions for BRD3 in the regulation of rRNA synthesis and potentially RNAPII-dependent gene expression that result in decreased cell proliferation. Together, our data highlight the contributions of BET protein modules to their interactomes allowing for a better understanding of pharmacological rewiring in response to JQ1.

Enhancer-gene rewiring in the pathogenesis of Quebec platelet disorder.

Quebec platelet disorder (QPD) is an autosomal dominant bleeding disorder with a unique, platelet-dependent, gain-of-function defect in fibrinolysis, without systemic fibrinolysis. The hallmark feature of QPD is a >100-fold overexpression of PLAU, specifically in megakaryocytes. This overexpression leads to a >100-fold increase in platelet stores of urokinase plasminogen activator (PLAU/uPA); subsequent plasmin-mediated degradation of diverse α-granule proteins; and platelet-dependent, accelerated fibrinolysis. The causative mutation is a 78-kb tandem duplication of PLAU. How this duplication causes megakaryocyte-specific PLAU overexpression is unknown. To investigate the mechanism that causes QPD, we used epigenomic profiling, comparative genomics, and chromatin conformation capture approaches to study PLAU regulation in cultured megakaryocytes from participants with QPD and unaffected controls. QPD duplication led to ectopic interactions between PLAU and a conserved megakaryocyte enhancer found within the same topologically associating domain (TAD). Our results support a unique disease mechanism whereby the reorganization of sub-TAD genome architecture results in a dramatic, cell-type-specific blood disorder phenotype.

X-linked myotubular myopathy is associated with epigenetic alterations and is ameliorated by HDAC inhibition.

X-linked myotubular myopathy (XLMTM) is a fatal neuromuscular disorder caused by loss of function mutations in MTM1. At present, there are no directed therapies for XLMTM, and incomplete understanding of disease pathomechanisms. To address these knowledge gaps, we performed a drug screen in mtm1 mutant zebrafish and identified four positive hits, including valproic acid, which functions as a potent suppressor of the mtm1 zebrafish phenotype via HDAC inhibition. We translated these findings to a mouse XLMTM model, and showed that valproic acid ameliorates the murine phenotype. These observations led us to interrogate the epigenome in Mtm1 knockout mice; we found increased DNA methylation, which is normalized with valproic acid, and likely mediated through aberrant 1-carbon metabolism. Finally, we made the unexpected observation that XLMTM patients share a distinct DNA methylation signature, suggesting that epigenetic alteration is a conserved disease feature amenable to therapeutic intervention.

Gene expression profiling of puberty-associated genes reveals abundant tissue and sex-specific changes across postnatal development.

The timing of human puberty is highly variable, sexually dimorphic, and associated with adverse health outcomes. Over 20 genes carrying rare mutations have been identified in known pubertal disorders, many of which encode critical components of the hypothalamic-pituitary-gonadal (HPG) axis. Recent genome-wide association studies (GWAS) have identified more than 100 candidate genes at loci associated with age at menarche or voice breaking in males. We know little about the spatial, temporal or postnatal expression patterns of the majority of these puberty-associated genes. Using a high-throughput and sensitive microfluidic quantitative PCR strategy, we profiled the gene expression patterns of the mouse orthologs of 178 puberty-associated genes in male and female mouse HPG axis tissues, the pineal gland, and the liver at five postnatal ages spanning the pubertal transition. The most dynamic gene expression changes were observed prior to puberty in all tissues. We detected known and novel tissue-enhanced gene expression patterns, with the hypothalamus expressing the largest number of the puberty-associated genes. Notably, over 40 puberty-associated genes in the pituitary gland showed sex-biased gene expression, most of which occurred peri-puberty. These sex-biased genes included the orthologs of candidate genes at GWAS loci that show sex-discordant effects on pubertal timing. Our findings provide new insight into the expression of puberty-associated genes and support the possibility that the pituitary plays a role in determining sex differences in the timing of puberty.

MACE: model based analysis of ChIP-exo.

Understanding the role of a given transcription factor (TF) in regulating gene expression requires precise mapping of its binding sites in the genome. Chromatin immunoprecipitation-exo, an emerging technique using λ exonuclease to digest TF unbound DNA after ChIP, is designed to reveal transcription factor binding site (TFBS) boundaries with near-single nucleotide resolution. Although ChIP-exo promises deeper insights into transcription regulation, no dedicated bioinformatics tool exists to leverage its advantages. Most ChIP-seq and ChIP-chip analytic methods are not tailored for ChIP-exo, and thus cannot take full advantage of high-resolution ChIP-exo data. Here we describe a novel analysis framework, termed MACE (model-based analysis of ChIP-exo) dedicated to ChIP-exo data analysis. The MACE workflow consists of four steps: (i) sequencing data normalization and bias correction; (ii) signal consolidation and noise reduction; (iii) single-nucleotide resolution border peak detection using the Chebyshev Inequality and (iv) border matching using the Gale-Shapley stable matching algorithm. When applied to published human CTCF, yeast Reb1 and our own mouse ONECUT1/HNF6 ChIP-exo data, MACE is able to define TFBSs with high sensitivity, specificity and spatial resolution, as evidenced by multiple criteria including motif enrichment, sequence conservation, direct sequence pileup, nucleosome positioning and open chromatin states. In addition, we show that the fundamental advance of MACE is the identification of two boundaries of a TFBS with high resolution, whereas other methods only report a single location of the same event. The two boundaries help elucidate the in vivo binding structure of a given TF, e.g. whether the TF may bind as dimers or in a complex with other co-factors.

Age, sex, and cell type-resolved hypothalamic gene expression across the pubertal transition in mice.

BACKGROUND: The hypothalamus plays a central role in regulating puberty. However, our knowledge of the postnatal gene regulatory networks that control the pubertal transition in males and females is incomplete. Here, we investigate the age-, sex- and cell-type-specific gene regulation in the hypothalamus across the pubertal transition. METHODS: We used RNA-seq to profile hypothalamic gene expression in male and female mice at five time points spanning the onset of puberty (postnatal days (PD) 12, 22, 27, 32, and 37). By combining this data with hypothalamic single nuclei RNA-seq data from pre- and postpubertal mice, we assigned gene expression changes to their most likely cell types of origin. In our colony, pubertal onset occurs earlier in male mice, allowing us to focus on genes whose expression is dynamic across ages and offset between sexes, and to explore the bases of sex effects. RESULTS: Our age-by-sex pattern of expression enriched for biological pathways involved hormone production, neuronal activation, and glial maturation. Additionally, we inferred a robust expansion of oligodendrocytes precursor cells into mature oligodendrocytes spanning the prepubertal (PD12) to peri-pubertal (PD27) timepoints. Using spatial transcriptomic data from postpubertal mice, we observed the lateral hypothalamic area and zona incerta were the most oligodendrocyte-rich regions and that these cells expressed genes known to be involved in pubertal regulation. CONCLUSION: Together, by incorporating multiple biological timepoints and using sex as a variable, we identified gene and cell-type changes that may participate in orchestrating the pubertal transition and provided a resource for future studies of postnatal hypothalamic gene regulation.

Untangling the roles of TOP2A and TOP2B in transcription and cancer.

Type II topoisomerases (TOP2) are conserved regulators of chromatin topology that catalyze reversible DNA double-strand breaks (DSBs) and are essential for maintaining genomic integrity in diverse dynamic processes such as transcription, replication, and cell division. While controlled TOP2-mediated DSBs are an elegant solution to topological constraints of DNA, DSBs also contribute to the emergence of chromosomal translocations and mutations that drive cancer. The central importance of TOP2 enzymes as frontline chemotherapeutic targets is well known; however, their precise biological functions and impact in cancer development are still poorly understood. In this review, we provide an updated overview of TOP2A and TOP2B in the regulation of chromatin topology and transcription, and discuss the recent discoveries linking TOP2 activities with cancer pathogenesis.

Postnatal developmental trajectory of sex-biased gene expression in the mouse pituitary gland.

BACKGROUND: The pituitary gland regulates essential physiological processes such as growth, pubertal onset, stress response, metabolism, reproduction, and lactation. While sex biases in these functions and hormone production have been described, the underlying identity, temporal deployment, and cell-type specificity of sex-biased pituitary gene regulatory networks are not fully understood. METHODS: To capture sex differences in pituitary gene regulation dynamics during postnatal development, we performed 3' untranslated region sequencing and small RNA sequencing to ascertain gene and microRNA expression, respectively, across five postnatal ages (postnatal days 12, 22, 27, 32, 37) that span the pubertal transition in female and male C57BL/6J mouse pituitaries (n = 5-6 biological replicates for each sex at each age). RESULTS: We observed over 900 instances of sex-biased gene expression and 17 sex-biased microRNAs, with the majority of sex differences occurring with puberty. Using miRNA-gene target interaction databases, we identified 18 sex-biased genes that were putative targets of 5 sex-biased microRNAs. In addition, by combining our bulk RNA-seq with publicly available male and female mouse pituitary single-nuclei RNA-seq data, we obtained evidence that cell-type proportion sex differences exist prior to puberty and persist post-puberty for three major hormone-producing cell types: somatotropes, lactotropes, and gonadotropes. Finally, we identified sex-biased genes in these three pituitary cell types after accounting for cell-type proportion differences between sexes. CONCLUSION: Our study reveals the identity and postnatal developmental trajectory of sex-biased gene expression in the mouse pituitary. This work also highlights the importance of considering sex biases in cell-type composition when understanding sex differences in the processes regulated by the pituitary gland.

Loss of Epigenetic Regulation Disrupts Lineage Integrity, Induces Aberrant Alveogenesis, and Promotes Breast Cancer.

Systematically investigating the scores of genes mutated in cancer and discerning disease drivers from inconsequential bystanders is a prerequisite for precision medicine but remains challenging. Here, we developed a somatic CRISPR/Cas9 mutagenesis screen to study 215 recurrent "long-tail" breast cancer genes, which revealed epigenetic regulation as a major tumor-suppressive mechanism. We report that components of the BAP1 and COMPASS-like complexes, including KMT2C/D, KDM6A, BAP1, and ASXL1/2 ("EpiDrivers"), cooperate with PIK3CAH1047R to transform mouse and human breast epithelial cells. Mechanistically, we find that activation of PIK3CAH1047R and concomitant EpiDriver loss triggered an alveolar-like lineage conversion of basal mammary epithelial cells and accelerated formation of luminal-like tumors, suggesting a basal origin for luminal tumors. EpiDriver mutations are found in ∼39% of human breast cancers, and ∼50% of ductal carcinoma in situ express casein, suggesting that lineage infidelity and alveogenic mimicry may significantly contribute to early steps of breast cancer etiology. SIGNIFICANCE: Infrequently mutated genes comprise most of the mutational burden in breast tumors but are poorly understood. In vivo CRISPR screening identified functional tumor suppressors that converged on epigenetic regulation. Loss of epigenetic regulators accelerated tumorigenesis and revealed lineage infidelity and aberrant expression of alveogenesis genes as potential early events in tumorigenesis. This article is highlighted in the In This Issue feature, p. 2711.

Dual Regulatory Functions of SUFU and Targetome of GLI2 in SHH Subgroup Medulloblastoma.

SUFU alterations are common in human Sonic Hedgehog (SHH) subgroup medulloblastoma (MB). However, its tumorigenic mechanisms have remained elusive. Here, we report that loss of Sufu alone is unable to induce MB formation in mice, due to insufficient Gli2 activation. Simultaneous loss of Spop, an E3 ubiquitin ligase targeting Gli2, restores robust Gli2 activation and induces rapid MB formation in Sufu knockout background. We also demonstrated a tumor-promoting role of Sufu in Smo-activated MB (∼60% of human SHH MB) by maintaining robust Gli activity. Having established Gli2 activation as a key driver of SHH MB, we report a comprehensive analysis of its targetome. Furthermore, we identified Atoh1 as a target and molecular accomplice of Gli2 that activates core SHH MB signature genes in a synergistic manner. Overall, our work establishes the dual role of SUFU in SHH MB and provides mechanistic insights into transcriptional regulation underlying Gli2-mediated SHH MB tumorigenesis.

Topoisomerase II beta interacts with cohesin and CTCF at topological domain borders.

BACKGROUND: Type II DNA topoisomerases (TOP2) regulate DNA topology by generating transient double stranded breaks during replication and transcription. Topoisomerase II beta (TOP2B) facilitates rapid gene expression and functions at the later stages of development and differentiation. To gain new insight into the genome biology of TOP2B, we used proteomics (BioID), chromatin immunoprecipitation, and high-throughput chromosome conformation capture (Hi-C) to identify novel proximal TOP2B protein interactions and characterize the genomic landscape of TOP2B binding at base pair resolution. RESULTS: Our human TOP2B proximal protein interaction network included members of the cohesin complex and nucleolar proteins associated with rDNA biology. TOP2B associates with DNase I hypersensitivity sites, allele-specific transcription factor (TF) binding, and evolutionarily conserved TF binding sites on the mouse genome. Approximately half of all CTCF/cohesion-bound regions coincided with TOP2B binding. Base pair resolution ChIP-exo mapping of TOP2B, CTCF, and cohesin sites revealed a striking structural ordering of these proteins along the genome relative to the CTCF motif. These ordered TOP2B-CTCF-cohesin sites flank the boundaries of topologically associating domains (TADs) with TOP2B positioned externally and cohesin internally to the domain loop. CONCLUSIONS: TOP2B is positioned to solve topological problems at diverse cis-regulatory elements and its occupancy is a highly ordered and prevalent feature of CTCF/cohesin binding sites that flank TADs.