Promoter-Intrinsic and Local Chromatin Features Determine Gene Repression in LADs.

It is largely unclear whether genes that are naturally embedded in lamina-associated domains (LADs) are inactive due to their chromatin environment or whether LADs are merely secondary to the lack of transcription. We show that hundreds of human promoters become active when moved from their native LAD position to a neutral context in the same cells, indicating that LADs form a repressive environment. Another set of promoters inside LADs is able to "escape" repression, although their transcription elongation is attenuated. By inserting reporters into thousands of genomic locations, we demonstrate that escaper promoters are intrinsically less sensitive to LAD repression. This is not simply explained by promoter strength but by the interplay between promoter sequence and local chromatin features that vary strongly across LADs. Enhancers also differ in their sensitivity to LAD chromatin. This work provides a general framework for the systematic understanding of gene regulation by repressive chromatin.

Systematic analysis of intrinsic enhancer-promoter compatibility in the mouse genome.

Gene expression is in part controlled by cis-regulatory elements (CREs) such as enhancers and repressive elements. Anecdotal evidence has indicated that a CRE and a promoter need to be biochemically compatible for promoter regulation to occur, but this compatibility has remained poorly characterized in mammalian cells. We used high-throughput combinatorial reporter assays to test thousands of CRE-promoter pairs from three Mb-sized genomic regions in mouse cells. This revealed that CREs vary substantially in their promoter compatibility, ranging from striking specificity to broad promiscuity. More than half of the tested CREs exhibit significant promoter selectivity. Housekeeping promoters tend to have similar CRE preferences, but other promoters exhibit a wide diversity of compatibilities. Higher-order transcription factors (TF) motif combinations may account for compatibility. CRE-promoter selectivity does not correlate with looping interactions in the native genomic context, suggesting that chromatin folding and compatibility are two orthogonal mechanisms that confer specificity to gene regulation.

Dissection of acute stimulus-inducible nucleosome remodeling in mammalian cells.

Accessibility of the genomic regulatory information is largely controlled by the nucleosome-organizing activity of transcription factors (TFs). While stimulus-induced TFs bind to genomic regions that are maintained accessible by lineage-determining TFs, they also increase accessibility of thousands of cis-regulatory elements. Nucleosome remodeling events underlying such changes and their interplay with basal positioning are unknown. Here, we devised a novel quantitative framework discriminating different types of nucleosome remodeling events in micrococcal nuclease ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) data sets and used it to analyze nucleosome dynamics at stimulus-regulated cis-regulatory elements. At enhancers, remodeling preferentially affected poorly positioned nucleosomes while sparing well-positioned nucleosomes flanking the enhancer core, indicating that inducible TFs do not suffice to overrule basal nucleosomal organization maintained by lineage-determining TFs. Remodeling events appeared to be combinatorially driven by multiple TFs, with distinct TFs showing, however, different remodeling efficiencies. Overall, these data provide a systematic view of the impact of stimulation on nucleosome organization and genome accessibility in mammalian cells.

SHMT2 inhibition disrupts the TCF3 transcriptional survival program in Burkitt lymphoma.

Burkitt lymphoma (BL) is an aggressive lymphoma type that is currently treated by intensive chemoimmunotherapy. Despite the favorable clinical outcome for most patients with BL, chemotherapy-related toxicity and disease relapse remain major clinical challenges, emphasizing the need for innovative therapies. Using genome-scale CRISPR-Cas9 screens, we identified B-cell receptor (BCR) signaling, specific transcriptional regulators, and one-carbon metabolism as vulnerabilities in BL. We focused on serine hydroxymethyltransferase 2 (SHMT2), a key enzyme in one-carbon metabolism. Inhibition of SHMT2 by either knockdown or pharmacological compounds induced anti-BL effects in vitro and in vivo. Mechanistically, SHMT2 inhibition led to a significant reduction of intracellular glycine and formate levels, which inhibited the mTOR pathway and thereby triggered autophagic degradation of the oncogenic transcription factor TCF3. Consequently, this led to a collapse of tonic BCR signaling, which is controlled by TCF3 and is essential for BL cell survival. In terms of clinical translation, we also identified drugs such as methotrexate that synergized with SHMT inhibitors. Overall, our study has uncovered the dependency landscape in BL, identified and validated SHMT2 as a drug target, and revealed a mechanistic link between SHMT2 and the transcriptional master regulator TCF3, opening up new perspectives for innovative therapies.

Rewiring of B cell receptor signaling by Epstein-Barr virus LMP2A.

Epstein-Barr virus (EBV) infects human B cells and reprograms them to allow virus replication and persistence. One key viral factor in this process is latent membrane protein 2A (LMP2A), which has been described as a B cell receptor (BCR) mimic promoting malignant transformation. However, how LMP2A signaling contributes to tumorigenesis remains elusive. By comparing LMP2A and BCR signaling in primary human B cells using phosphoproteomics and transcriptome profiling, we identified molecular mechanisms through which LMP2A affects B cell biology. Consistent with the literature, we found that LMP2A mimics a subset of BCR signaling events, including tyrosine phosphorylation of the kinase SYK, the calcium initiation complex consisting of BLNK, BTK, and PLCγ2, and its downstream transcription factor NFAT. However, the majority of LMP2A-induced signaling events markedly differed from those induced by BCR stimulation. These included differential phosphorylation of kinases, phosphatases, adaptor proteins, transcription factors such as nuclear factor κB (NF-κB) and TCF3, as well as widespread changes in the transcriptional output of LMP2A-expressing B cells. LMP2A affected apoptosis and cell-cycle checkpoints by dysregulating the expression of apoptosis regulators such as BCl-xL and the tumor suppressor retinoblastoma-associated protein 1 (RB1). LMP2A cooperated with MYC and mutant cyclin D3, two oncogenic drivers of Burkitt lymphoma, to promote proliferation and survival of primary human B cells by counteracting MYC-induced apoptosis and by inhibiting RB1 function, thereby promoting cell-cycle progression. Our results indicate that LMP2A is not a pure BCR mimic but rather rewires intracellular signaling in EBV-infected B cells that optimizes cell survival and proliferation, setting the stage for oncogenic transformation.

Thrombopoietin signaling to chromatin elicits rapid and pervasive epigenome remodeling within poised chromatin architectures.

Thrombopoietin (TPO) is a critical cytokine regulating hematopoietic stem cell maintenance and differentiation into the megakaryocytic lineage. However, the transcriptional and chromatin dynamics elicited by TPO signaling are poorly understood. Here, we study the immediate early transcriptional and cis-regulatory responses to TPO in hematopoietic stem/progenitor cells (HSPCs) and use this paradigm of cytokine signaling to chromatin to dissect the relation between cis- regulatory activity and chromatin architecture. We show that TPO profoundly alters the transcriptome of HSPCs, with key hematopoietic regulators being transcriptionally repressed within 30 minutes of TPO. By examining cis-regulatory dynamics and chromatin architectures, we demonstrate that these changes are accompanied by rapid and extensive epigenome remodeling of cis-regulatory landscapes that is spatially coordinated within topologically associating domains (TADs). Moreover, TPO-responsive enhancers are spatially clustered and engage in preferential homotypic intra- and inter-TAD interactions that are largely refractory to TPO signaling. By further examining the link between cis-regulatory dynamics and chromatin looping, we show that rapid modulation of cis-regulatory activity is largely independent of chromatin looping dynamics. Finally, we show that, although activated and repressed cis-regulatory elements share remarkably similar DNA sequence compositions, transcription factor binding patterns accurately predict rapid cis-regulatory responses to TPO.

HSP90 promotes Burkitt lymphoma cell survival by maintaining tonic B-cell receptor signaling.

Burkitt lymphoma (BL) is an aggressive B-cell neoplasm that is currently treated by intensive chemotherapy in combination with anti-CD20 antibodies. Because of their toxicity, current treatment regimens are often not suitable for elderly patients or for patients in developing countries where BL is endemic. Targeted therapies for BL are therefore needed. In this study, we performed a compound screen in 17 BL cell lines to identify small molecule inhibitors affecting cell survival. We found that inhibitors of heat shock protein 90 (HSP90) induced apoptosis in BL cells in vitro at concentrations that did not affect normal B cells. By global proteomic and phosphoproteomic profiling, we show that, in BL, HSP90 inhibition compromises the activity of the pivotal B-cell antigen receptor (BCR)-proximal effector spleen tyrosine kinase (SYK), which we identified as an HSP90 client protein. Consistently, expression of constitutively active TEL-SYK counteracted the apoptotic effect of HSP90 inhibition. Together, our results demonstrate that HSP90 inhibition impairs BL cell survival by interfering with tonic BCR signaling, thus providing a molecular rationale for the use of HSP90 inhibitors in the treatment of BL.

Elucidation of tonic and activated B-cell receptor signaling in Burkitt's lymphoma provides insights into regulation of cell survival.

Burkitt's lymphoma (BL) is a highly proliferative B-cell neoplasm and is treated with intensive chemotherapy that, because of its toxicity, is often not suitable for the elderly or for patients with endemic BL in developing countries. BL cell survival relies on signals transduced by B-cell antigen receptors (BCRs). However, tonic as well as activated BCR signaling networks and their relevance for targeted therapies in BL remain elusive. We have systematically characterized and compared tonic and activated BCR signaling in BL by quantitative phosphoproteomics to identify novel BCR effectors and potential drug targets. We identified and quantified ∼16,000 phospho-sites in BL cells. Among these sites, 909 were related to tonic BCR signaling, whereas 984 phospho-sites were regulated upon BCR engagement. The majority of the identified BCR signaling effectors have not been described in the context of B cells or lymphomas yet. Most of these newly identified BCR effectors are predicted to be involved in the regulation of kinases, transcription, and cytoskeleton dynamics. Although tonic and activated BCR signaling shared a considerable number of effector proteins, we identified distinct phosphorylation events in tonic BCR signaling. We investigated the functional relevance of some newly identified BCR effectors and show that ACTN4 and ARFGEF2, which have been described as regulators of membrane-trafficking and cytoskeleton-related processes, respectively, are crucial for BL cell survival. Thus, this study provides a comprehensive dataset for tonic and activated BCR signaling and identifies effector proteins that may be relevant for BL cell survival and thus may help to develop new BL treatments.

During Drosophila disc regeneration, JAK/STAT coordinates cell proliferation with Dilp8-mediated developmental delay.

Regeneration of fragmented Drosophila imaginal discs occurs in an epimorphic manner involving local cell proliferation at the wound site. After disc fragmentation, cells at the wound site activate a restoration program through wound healing, regenerative cell proliferation, and repatterning of the tissue. However, the interplay of signaling cascades driving these early reprogramming steps is not well-understood. Here, we profiled the transcriptome of regenerating cells in the early phase within 24 h after wounding. We found that JAK/STAT signaling becomes activated at the wound site and promotes regenerative cell proliferation in cooperation with Wingless (Wg) signaling. In addition, we showed that the expression of Drosophila insulin-like peptide 8 (dilp8), which encodes a paracrine peptide to delay the onset of pupariation, is controlled by JAK/STAT signaling in early regenerating discs. Our findings suggest that JAK/STAT signaling plays a pivotal role in coordinating regenerative disc growth with organismal developmental timing.

DNAshapeR: an R/Bioconductor package for DNA shape prediction and feature encoding.

UNLABELLED: DNAshapeR predicts DNA shape features in an ultra-fast, high-throughput manner from genomic sequencing data. The package takes either nucleotide sequence or genomic coordinates as input and generates various graphical representations for visualization and further analysis. DNAshapeR further encodes DNA sequence and shape features as user-defined combinations of k-mer and DNA shape features. The resulting feature matrices can be readily used as input of various machine learning software packages for further modeling studies. AVAILABILITY AND IMPLEMENTATION: The DNAshapeR software package was implemented in the statistical programming language R and is freely available through the Bioconductor project at https://www.bioconductor.org/packages/devel/bioc/html/DNAshapeR.html and at the GitHub developer site, http://tsupeichiu.github.io/DNAshapeR/ CONTACT: rohs@usc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PcG-mediated higher-order chromatin structures modulate replication programs at the Drosophila BX-C.

Polycomb group proteins (PcG) exert conserved epigenetic functions that convey maintenance of repressed transcriptional states, via post-translational histone modifications and high order structure formation. During S-phase, in order to preserve cell identity, in addition to DNA information, PcG-chromatin-mediated epigenetic signatures need to be duplicated requiring a tight coordination between PcG proteins and replication programs. However, the interconnection between replication timing control and PcG functions remains unknown. Using Drosophila embryonic cell lines, we find that, while presence of specific PcG complexes and underlying transcription state are not the sole determinants of cellular replication timing, PcG-mediated higher-order structures appear to dictate the timing of replication and maintenance of the silenced state. Using published datasets we show that PRC1, PRC2, and PhoRC complexes differently correlate with replication timing of their targets. In the fully repressed BX-C, loss of function experiments revealed a synergistic role for PcG proteins in the maintenance of replication programs through the mediation of higher-order structures. Accordingly, replication timing analysis performed on two Drosophila cell lines differing for BX-C gene expression states, PcG distribution, and chromatin domain conformation revealed a cell-type-specific replication program that mirrors lineage-specific BX-C higher-order structures. Our work suggests that PcG complexes, by regulating higher-order chromatin structure at their target sites, contribute to the definition and the maintenance of genomic structural domains where genes showing the same epigenetic state replicate at the same time.

Sensitive and highly resolved identification of RNA-protein interaction sites in PAR-CLIP data.

BACKGROUND: PAR-CLIP is a recently developed Next Generation Sequencing-based method enabling transcriptome-wide identification of interaction sites between RNA and RNA-binding proteins. The PAR-CLIP procedure induces specific base transitions that originate from sites of RNA-protein interactions and can therefore guide the identification of binding sites. However, additional sources of transitions, such as cell type-specific SNPs and sequencing errors, challenge the inference of binding sites and suitable statistical approaches are crucial to control false discovery rates. In addition, a highly resolved delineation of binding sites followed by an extensive downstream analysis is necessary for a comprehensive characterization of the protein binding preferences and the subsequent design of validation experiments. RESULTS: We present a statistical and computational framework for PAR-CLIP data analysis. We developed a sensitive transition-centered algorithm specifically designed to resolve protein binding sites at high resolution in PAR-CLIP data. Our method employes a Bayesian network approach to associate posterior log-odds with the observed transitions, providing an overall quantification of the confidence in RNA-protein interaction. We use published PAR-CLIP data to demonstrate the advantages of our approach, which compares favorably with alternative algorithms. Lastly, by integrating RNA-Seq data we compute conservative experimentally-based false discovery rates of our method and demonstrate the high precision of our strategy. CONCLUSIONS: Our method is implemented in the R package wavClusteR 2.0. The package is distributed under the GPL-2 license and is available from BioConductor at http://www.bioconductor.org/packages/devel/bioc/html/wavClusteR.html .

Combinatorial modeling of chromatin features quantitatively predicts DNA replication timing in Drosophila.

In metazoans, each cell type follows a characteristic, spatio-temporally regulated DNA replication program. Histone modifications (HMs) and chromatin binding proteins (CBPs) are fundamental for a faithful progression and completion of this process. However, no individual HM is strictly indispensable for origin function, suggesting that HMs may act combinatorially in analogy to the histone code hypothesis for transcriptional regulation. In contrast to gene expression however, the relationship between combinations of chromatin features and DNA replication timing has not yet been demonstrated. Here, by exploiting a comprehensive data collection consisting of 95 CBPs and HMs we investigated their combinatorial potential for the prediction of DNA replication timing in Drosophila using quantitative statistical models. We found that while combinations of CBPs exhibit moderate predictive power for replication timing, pairwise interactions between HMs lead to accurate predictions genome-wide that can be locally further improved by CBPs. Independent feature importance and model analyses led us to derive a simplified, biologically interpretable model of the relationship between chromatin landscape and replication timing reaching 80% of the full model accuracy using six model terms. Finally, we show that pairwise combinations of HMs are able to predict differential DNA replication timing across different cell types. All in all, our work provides support to the existence of combinatorial HM patterns for DNA replication and reveal cell-type independent key elements thereof, whose experimental investigation might contribute to elucidate the regulatory mode of this fundamental cellular process.

Mixture models and wavelet transforms reveal high confidence RNA-protein interaction sites in MOV10 PAR-CLIP data.

The Photo-Activatable Ribonucleoside-enhanced CrossLinking and ImmunoPrecipitation (PAR-CLIP) method was recently developed for global identification of RNAs interacting with proteins. The strength of this versatile method results from induction of specific T to C transitions at sites of interaction. However, current analytical tools do not distinguish between non-experimentally and experimentally induced transitions. Furthermore, geometric properties at potential binding sites are not taken into account. To surmount these shortcomings, we developed a two-step algorithm consisting of a non-parametric two-component mixture model and a wavelet-based peak calling procedure. Our algorithm can reduce the number of false positives up to 24% thereby identifying high confidence interaction sites. We successfully employed this approach in conjunction with a modified PAR-CLIP protocol to study the functional role of nuclear Moloney leukemia virus 10, a putative RNA helicase interacting with Argonaute2 and Polycomb. Our method, available as the R package wavClusteR, is generally applicable to any substitution-based inference problem in genomics.

Rknots: topological analysis of knotted biopolymers with R.

MOTIVATION: Rknots is a flexible R package providing tools for the detection and characterization of topological knots in biological polymers. The package is well documented and provides a simple syntax for data import and preprocessing, structure reduction, topological analyses and 2D and 3D visualization. Remarkably, Rknots is not limited to protein knots and allows researchers from interdisciplinary fields to analyze different topological structures and to develop simple yet fully custom pipelines.

PDP1 is a key metabolic gatekeeper and modulator of drug resistance in FLT3-ITD-positive acute myeloid leukemia.

High metabolic flexibility is pivotal for the persistence and therapy resistance of acute myeloid leukemia (AML). In 20-30% of AML patients, activating mutations of FLT3, specifically FLT3-ITD, are key therapeutic targets. Here, we investigated the influence of FLT3-ITD on AML metabolism. Nuclear Magnetic Resonance (NMR) profiling showed enhanced reshuffling of pyruvate towards the tricarboxylic acid (TCA) cycle, suggesting an increased activity of the pyruvate dehydrogenase complex (PDC). Consistently, FLT3-ITD-positive cells expressed high levels of PDP1, an activator of the PDC. Combining endogenous tagging of PDP1 with genome-wide CRISPR screens revealed that FLT3-ITD induces PDP1 expression through the RAS signaling axis. PDP1 knockdown resulted in reduced cellular respiration thereby impairing the proliferation of only FLT3-ITD cells. These cells continued to depend on PDP1, even in hypoxic conditions, and unlike FLT3-ITD-negative cells, they exhibited a rapid, PDP1-dependent revival of their respiratory capacity during reoxygenation. Moreover, we show that PDP1 modifies the response to FLT3 inhibition. Upon incubation with the FLT3 tyrosine kinase inhibitor quizartinib (AC220), PDP1 persisted or was upregulated, resulting in a further shift of glucose/pyruvate metabolism towards the TCA cycle. Overexpression of PDP1 enhanced, while PDP1 depletion diminished AC220 resistance in cell lines and peripheral blasts from an AC220-resistant AML patient in vivo. In conclusion, FLT3-ITD assures the expression of PDP1, a pivotal metabolic regulator that enhances oxidative glucose metabolism and drug resistance. Hence, PDP1 emerges as a potentially targetable vulnerability in the management of AML.

The proteogenomic subtypes of acute myeloid leukemia.

Acute myeloid leukemia (AML) is an aggressive blood cancer with a poor prognosis. We report a comprehensive proteogenomic analysis of bone marrow biopsies from 252 uniformly treated AML patients to elucidate the molecular pathophysiology of AML in order to inform future diagnostic and therapeutic approaches. In addition to in-depth quantitative proteomics, our analysis includes cytogenetic profiling and DNA/RNA sequencing. We identify five proteomic AML subtypes, each reflecting specific biological features spanning genomic boundaries. Two of these proteomic subtypes correlate with patient outcome, but none is exclusively associated with specific genomic aberrations. Remarkably, one subtype (Mito-AML), which is captured only in the proteome, is characterized by high expression of mitochondrial proteins and confers poor outcome, with reduced remission rate and shorter overall survival on treatment with intensive induction chemotherapy. Functional analyses reveal that Mito-AML is metabolically wired toward stronger complex I-dependent respiration and is more responsive to treatment with the BCL2 inhibitor venetoclax.

MicroRNA signature of malignant mesothelioma with potential diagnostic and prognostic implications.

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of target genes, and may behave as oncogenes or tumor suppressors. Human malignant mesothelioma is an asbestos-related cancer, with poor prognosis and low median survival. Here we report, for the first time, a cross-evaluation of miRNA expression in mesothelioma (MPP-89, REN) and human mesothelial cells (HMC-telomerase reverse transcriptase). Microarray profiling, confirmed by real-time quantitative RT-PCR, revealed a differential expression of miRNAs between mesothelioma and mesothelial cells. In addition, a computational analysis combining miRNA and gene expression profiles allowed the accurate prediction of genes potentially targeted by dysregulated miRNAs. Several predicted genes belong to terms of Gene Ontology (GO) that are associated with the development and progression of mesothelioma. This suggests that miRNAs may be key players in mesothelioma oncogenesis. We further investigated miRNA expression on a panel of 24 mesothelioma specimens, representative of the three histotypes (epithelioid, biphasic, and sarcomatoid), by quantitative RT-PCR. The expression of miR-17-5p, miR-21, miR-29a, miR-30c, miR-30e-5p, miR-106a, and miR-143 was significantly associated with the histopathological subtypes. Notably, the reduced expression of two miRNAs (miR-17-5p and miR-30c) correlated with better survival of patients with sarcomatoid subtype. Our preliminary analysis points at miRNAs as potential diagnostic and prognostic markers of mesothelioma, and suggests novel tools for the therapy of this malignancy.

High-throughput identification of human SNPs affecting regulatory element activity.

Most of the millions of SNPs in the human genome are non-coding, and many overlap with putative regulatory elements. Genome-wide association studies (GWAS) have linked many of these SNPs to human traits or to gene expression levels, but rarely with sufficient resolution to identify the causal SNPs. Functional screens based on reporter assays have previously been of insufficient throughput to test the vast space of SNPs for possible effects on regulatory element activity. Here we leveraged the throughput and resolution of the survey of regulatory elements (SuRE) reporter technology to survey the effect of 5.9 million SNPs, including 57% of the known common SNPs, on enhancer and promoter activity. We identified more than 30,000 SNPs that alter the activity of putative regulatory elements, partially in a cell-type-specific manner. Integration of this dataset with GWAS results may help to pinpoint SNPs that underlie human traits.

Selective inactivation of hypomethylating agents by SAMHD1 provides a rationale for therapeutic stratification in AML.

Hypomethylating agents decitabine and azacytidine are regarded as interchangeable in the treatment of acute myeloid leukemia (AML). However, their mechanisms of action remain incompletely understood, and predictive biomarkers for HMA efficacy are lacking. Here, we show that the bioactive metabolite decitabine triphosphate, but not azacytidine triphosphate, functions as activator and substrate of the triphosphohydrolase SAMHD1 and is subject to SAMHD1-mediated inactivation. Retrospective immunohistochemical analysis of bone marrow specimens from AML patients at diagnosis revealed that SAMHD1 expression in leukemic cells inversely correlates with clinical response to decitabine, but not to azacytidine. SAMHD1 ablation increases the antileukemic activity of decitabine in AML cell lines, primary leukemic blasts, and xenograft models. AML cells acquire resistance to decitabine partly by SAMHD1 up-regulation. Together, our data suggest that SAMHD1 is a biomarker for the stratified use of hypomethylating agents in AML patients and a potential target for the treatment of decitabine-resistant leukemia.