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.

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.

MON-725 Transcriptome Profiling in Postnatal Pituitary Gland Identifies Cell Type-Driven Sex-Specific Changes

The pituitary gland is integral to the regulation of growth, metabolism, puberty, reproduction, and stress responses. Previously, we found that many genes associated with age-at-menarche in genome-wide association studies (GWAS) displayed increasingly sex-biased expression across the pubertal transition in the mouse pituitary. However, whether this trend exists beyond puberty-related genes was not known. In addition, the regulatory mechanisms underlying these gene expression changes remained to be explored. To answer these questions, we profiled the transcriptome, including microRNAs, of mouse pituitary in both sexes across pubertal transition in an unbiased manner and leveraged a recently published pituitary single cell transcriptome to explore cellular composition changes. We found that the most dynamic temporal changes in both mRNA and miRNA expression occur prior to puberty, underscoring a role for regulation of early pituitary postnatal development. We also observed ~900 genes displaying sex-biased expression patterns, arising during early development and becoming increasingly biased across puberty, including known sex-biased genes such as Fshb and Lhb. However, sex differences in miRNA expression are less pronounced, only 13 miRNAs were found to be sex-biased, suggesting lesser contribution of miRNAs to sex-biased gene expression relative to other forms of regulation. To assess whether pituitary cellular composition could underlie changes in gene expression across pubertal transition, we performed single cell deconvolution of our bulk pituitary gland gene expression. Interestingly, we found that sex differences in cell proportions were estimated to emerge across puberty: a greater proportion of lactotropes was found among females, and greater proportions of gonadotropes and somatotropes were found among males. We observed sex-biased expression patterns of marker genes for these cell types, including Prl, Fshb, and Gh. This finding suggests that cell proportion differences between sexes likely contribute to whole pituitary transcriptome changes we observed, however, to what extent remains to be studied. Together our study indicates that miRNAs play a substantial role in regulation of pituitary postnatal development but that differences in cellular composition may contribute more robustly to sex-biased gene expression.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-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 promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes

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.

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.

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.

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.

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.

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.

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.