Differential Expression Enrichment Tool (DEET): an interactive atlas of human differential gene expression.

Differential gene expression analysis using RNA sequencing (RNA-seq) data is a standard approach for making biological discoveries. Ongoing large-scale efforts to process and normalize publicly available gene expression data enable rapid and systematic reanalysis. While several powerful tools systematically process RNA-seq data, enabling their reanalysis, few resources systematically recompute differentially expressed genes (DEGs) generated from individual studies. We developed a robust differential expression analysis pipeline to recompute 3162 human DEG lists from The Cancer Genome Atlas, Genotype-Tissue Expression Consortium, and 142 studies within the Sequence Read Archive. After measuring the accuracy of the recomputed DEG lists, we built the Differential Expression Enrichment Tool (DEET), which enables users to interact with the recomputed DEG lists. DEET, available through CRAN and RShiny, systematically queries which of the recomputed DEG lists share similar genes, pathways, and TF targets to their own gene lists. DEET identifies relevant studies based on shared results with the user's gene lists, aiding in hypothesis generation and data-driven literature review.

Alternative polyadenylation is a determinant of oncogenic Ras function.

Alternative polyadenylation of mRNA has important but poorly understood roles in development and cancer. Activating mutations in the Ras oncogene are common drivers of many human cancers. From a screen for enhancers of activated Ras (let-60) in Caenorhabditis elegans, we identified cfim-1, a subunit of the alternative polyadenylation machinery. Ablation of cfim-1 increased penetrance of the multivulva phenotype in let-60/Ras gain-of-function (gf) mutants. Depletion of the human cfim-1 ortholog CFIm25/NUDT21 in cancer cells with KRAS mutations increased their migration and stimulated an epithelial-to-mesenchymal transition. CFIm25-depleted cells and cfim-1 mutants displayed biased placement of poly(A) tails to more proximal sites in many conserved transcripts. Functional analysis of these transcripts identified the multidrug resistance protein mrp-5/ABCC1 as a previously unidentified regulator of C. elegans vulva development and cell migration in human cells through alternative 3'UTR usage. Our observations demonstrate a conserved functional role for alternative polyadenylation in oncogenic Ras function.

Lineage-defined leiomyosarcoma subtypes emerge years before diagnosis and determine patient survival.

Leiomyosarcomas (LMS) are genetically heterogeneous tumors differentiating along smooth muscle lines. Currently, LMS treatment is not informed by molecular subtyping and is associated with highly variable survival. While disease site continues to dictate clinical management, the contribution of genetic factors to LMS subtype, origins, and timing are unknown. Here we analyze 70 genomes and 130 transcriptomes of LMS, including multiple tumor regions and paired metastases. Molecular profiling highlight the very early origins of LMS. We uncover three specific subtypes of LMS that likely develop from distinct lineages of smooth muscle cells. Of these, dedifferentiated LMS with high immune infiltration and tumors primarily of gynecological origin harbor genomic dystrophin deletions and/or loss of dystrophin expression, acquire the highest burden of genomic mutation, and are associated with worse survival. Homologous recombination defects lead to genome-wide mutational signatures, and a corresponding sensitivity to PARP trappers and other DNA damage response inhibitors, suggesting a promising therapeutic strategy for LMS. Finally, by phylogenetic reconstruction, we present evidence that clones seeding lethal metastases arise decades prior to LMS diagnosis.

The pregnant myometrium is epigenetically activated at contractility-driving gene loci prior to the onset of labor in mice.

During gestation, uterine smooth muscle cells transition from a state of quiescence to one of contractility, but the molecular mechanisms underlying this transition at a genomic level are not well-known. To better understand these events, we evaluated the epigenetic landscape of the mouse myometrium during the pregnant, laboring, and postpartum stages. We generated gestational time point-specific enrichment profiles for histone H3 acetylation on lysine residue 27 (H3K27ac), histone H3 trimethylation of lysine residue 4 (H3K4me3), and RNA polymerase II (RNAPII) occupancy by chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq), as well as gene expression profiles by total RNA-sequencing (RNA-seq). Our findings reveal that 533 genes, including known contractility-driving genes (Gap junction alpha 1 [Gja1], FBJ osteosarcoma oncogene [Fos], Fos-like antigen 2 [Fosl2], Oxytocin receptor [Oxtr], and Prostaglandin G/H synthase 2 (Ptgs2), for example), are up-regulated at day 19 during active labor because of an increase in transcription at gene bodies. Labor-associated promoters and putative intergenic enhancers, however, are epigenetically activated as early as day 15, by which point the majority of genome-wide H3K27ac or H3K4me3 peaks present in term laboring tissue is already established. Despite this early exhibited histone signature, increased noncoding enhancer RNA (eRNA) production at putative intergenic enhancers and recruitment of RNAPII to the gene bodies of labor-associated loci were detected only during labor. Our findings indicate that epigenetic activation of the myometrial genome precedes active labor by at least 4 days in the mouse model, suggesting that the myometrium is poised for rapid activation of contraction-associated genes in order to exit the state of quiescence.

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.

Correction of a splicing defect in a mouse model of congenital muscular dystrophy type 1A using a homology-directed-repair-independent mechanism.

Splice-site defects account for about 10% of pathogenic mutations that cause Mendelian diseases. Prevalence is higher in neuromuscular disorders (NMDs), owing to the unusually large size and multi-exonic nature of genes encoding muscle structural proteins. Therapeutic genome editing to correct disease-causing splice-site mutations has been accomplished only through the homology-directed repair pathway, which is extremely inefficient in postmitotic tissues such as skeletal muscle. Here we describe a strategy using nonhomologous end-joining (NHEJ) to correct a pathogenic splice-site mutation. As a proof of principle, we focus on congenital muscular dystrophy type 1A (MDC1A), which is characterized by severe muscle wasting and paralysis. Specifically, we correct a splice-site mutation that causes the exclusion of exon 2 from Lama2 mRNA and the truncation of Lama2 protein in the dy2J/dy2J mouse model of MDC1A. Through systemic delivery of adeno-associated virus (AAV) carrying clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 genome-editing components, we simultaneously excise an intronic region containing the mutation and create a functional donor splice site through NHEJ. This strategy leads to the inclusion of exon 2 in the Lama2 transcript and restoration of full-length Lama2 protein. Treated dy2J/dy2J mice display substantial improvement in muscle histopathology and function without signs of paralysis.