scispace - formally typeset
Open AccessPosted ContentDOI

Myeloid leukemia vulnerabilities at CTCF-enriched long noncoding RNA loci

Reads0
Chats0
TLDR
CTCF-enriched lncRNA loci are described as a subclass of functional genetic elements exemplified by MYNRL15, a pan-myeloid leukemia dependency identified by an lnc RNA-based CRISPRi screen, and it is shown that crucial DNA elements in the locus mediate its phenotype, triggering chromatin reorganization and downregulation of cancer dependency genes upon perturbation.
Abstract
The noncoding genome presents a largely untapped source of biological insights, including thousands of long noncoding RNA (lncRNA) loci. While some produce bona fide lncRNAs, others exert transcript-independent cis-regulatory effects, and the lack of predictive features renders mechanistic dissection challenging. Here, we describe CTCF-enriched lncRNA loci (C-LNC) as a subclass of functional genetic elements exemplified by MYNRL15, a pan-myeloid leukemia dependency identified by an lncRNA-based CRISPRi screen. MYNRL15 perturbation selectively impairs acute myeloid leukemia (AML) cells over hematopoietic stem / progenitor cells in vitro, and depletes AML xenografts in vivo. Mechanistically, we show that crucial DNA elements in the locus mediate its phenotype, triggering chromatin reorganization and downregulation of cancer dependency genes upon perturbation. Elevated CTCF density distinguishes MYNRL15 and 531 other lncRNA loci in K562 cells, of which 43-54% associate with clinical aspects of AML and 18.4% are functionally required for leukemia maintenance. Curated C-LNC catalogs in other cell types will help refine the search for noncoding oncogenic vulnerabilities in AML and other malignancies.

read more

Content maybe subject to copyright    Report

Myeloid leukemia vulnerabilities at CTCF-enriched
long noncoding RNA loci
Michelle Ng
Martin Luther University Halle-Wittenberg
Lonneke Verboon
Martin Luther University Halle-Wittenberg
Hasan Issa
Martin Luther University Halle-Wittenberg
Raj Bhayadia
Martin Luther University Halle-Wittenberg
Oriol Alejo-Valle
Martin Luther University Halle-Wittenberg https://orcid.org/0000-0002-0988-7416
Dorit Borchert
Hannover Medical School
Konstantin Schuschel
Martin Luther University Halle-Wittenberg
Eniko Regenyi
Free University of Berlin
Dirk Reinhardt
University Hospital Essen
Marie-Laure Yaspo
Max Planck Institute for Molecular Genetics
Dirk Heckl
Martin Luther University Halle-Wittenberg
Jan-Henning Klusmann ( jan-henning.klusmann@kgu.de )
Goethe University Frankfurt https://orcid.org/0000-0002-1070-0727
Genetics Article
Keywords: long noncoding RNA, C-LNC, MYNRL15
Posted Date: August 9th, 2021
DOI: https://doi.org/10.21203/rs.3.rs-727909/v1

License: This work is licensed under a Creative Commons Attribution 4.0 International License. 
Read Full License

1
Myeloid leukemia vulnerabilities at CTCF-enriched long noncoding RNA loci
Abstract word count: 165
Main text word count: 2656
References: 58
Running title: MYNRL15 in myeloid leukemia
Keywords: long noncoding RNA, myeloid leukemia, CRISPR, chromatin architecture
Abstract 1
The noncoding genome presents a largely untapped source of biological insights, including 2
thousands of long noncoding RNA (lncRNA) loci. While some produce bona fide lncRNAs, 3
others exert transcript-independent cis-regulatory effects, and the lack of predictive features 4
renders mechanistic dissection challenging. Here, we describe CTCF-enriched lncRNA loci 5
(C-LNC) as a subclass of functional genetic elements exemplified by MYNRL15, a pan-mye-6
loid leukemia dependency identified by an lncRNA-based CRISPRi screen. MYNRL15 per-7
turbation selectively impairs acute myeloid leukemia (AML) cells over hematopoietic stem / 8
progenitor cells in vitro, and depletes AML xenografts in vivo. Mechanistically, we show that 9
crucial DNA elements in the locus mediate its phenotype, triggering chromatin reorganization 10
and downregulation of cancer dependency genes upon perturbation. Elevated CTCF density 11
distinguishes MYNRL15 and 531 other lncRNA loci in K562 cells, of which 43-54% associate 12
with clinical aspects of AML and 18.4% are functionally required for leukemia maintenance. 13
Curated C-LNC catalogs in other cell types will help refine the search for noncoding onco-14
genic vulnerabilities in AML and other malignancies.
15

2
Main 16
It becomes increasingly clear that the 98% of the human genome that does not encode pro-17
tein nonetheless contains a wide range of functional elements that are vital for cellular home-18
ostasis
1,2
. These include cis-regulatory elements such as enhancers and promoters, insula-19
tors and other determinants of genome topology, as well as a large number and variety of 20
non-protein-coding transcripts. Long noncoding RNAs (lncRNAs) in particular comprise a 21
substantial portion of the noncoding transcriptome
3-5
and in recent years, have emerged as 22
important players in diverse cellular processes and contexts
6-8
. The hematopoietic system is 23
no exception, where lncRNAs have been described to regulate cell programming and fate
9
, 24
and where their dysregulation has been tied to malignancy
10-16
. LncRNAs present a signifi-25
cant opportunity to extend our understanding of human health and disease; however, the 26
fact remains that the vast majority of lncRNA loci lack functional characterization, and may 27
regulate cellular behaviour in ways yet unknown. Indeed, characterization is often a difficult 28
process complicated by cis-regulatory mechanisms unrelated to the transcriptional product
17-
29
23
. Improved functional classification systems are imperative for expediting investigations into 30
lncRNA determinants of pathophysiology, including the search for noncoding oncogenic vul-31
nerabilities. 32
CRISPRi screens of HSPC/AML lncRNAs identify MYNRL15 as a leukemia dependency 33
Aiming to identify lncRNAs that contribute to myeloid malignancy, we began by analyzing a 34
noncoding RNA expression atlas of the human blood system encompassing hematopoietic 35
stem cells (HSCs) and their differentiated progeny, as well as pediatric acute myeloid leuke-36
mia (AML) samples
16
. In addition to stem cell signatures reminiscent of those previously es-37
tablished for protein-coding genes
24-27
, we discovered progenitor- and AML subtype-associ-38
ated lncRNA profiles that could potentially serve as leukemia-specific targets, given their ab-39
sence in HSCs (Fig. 1a). To probe this resource for functionality and find novel AML vulnera-40
bilities, we conducted a CRISPRi-based dropout screen of 480 lncRNA genes from 8 distinct 41
signatures in 6 human leukemia cell lines (Fig. 1b). Five were selected to represent relevant 42

3
cytogenetic subgroups of AML (ML-2, NOMO-1 [KMT2A-rearreanged], SKNO-1, KASUMI-1 43
[standard risk with t(8:21)], M-07E [high risk with inv(16)(p13.3q24.3)]), and we also included 44
the well-studied erythroleukemia line K562. One candidate emerged as crucial in all six cell 45
lines – AC068831.3 (ID: ENSG00000224441 in Ensembl v91 [release 12/2017]), which we 46
renamed MYNRL15 (myeloid leukemia noncoding regulatory locus on chromosome 15; Fig.
47
1c-d, Extended Data Fig. 1). 48
MYNRL15 is a low-abundance, nuclear-enriched transcript (Extended Data Fig. 2a-b) origi-49
nating from chromosome 15, where it is flanked by two protein-coding genes: UNC45A and 50
HDDC3 (Fig. 1c). Given the local effect of the CRISPRi system on nearby genes (Extended 51
Data Fig. 2c), a range of gain- and loss-of-function approaches were necessary in order to 52
delineate the source of the MYNRL15 knockdown phenotype (Fig. 1e-f, Extended Data Fig. 53
2d-g). While CRISPR mediated excision of MYNRL15 recapitulated the effect produced by 54
CRISPRi, repression of the transcript via RNAi and LNA-gapmeRs had little impact on cell 55
viability (Fig. 1e, Extended Data Fig. 2d-f). Both protein-coding neighbors were also dispen-56
sable, as determined via CRISPR-Cas9 mediated knockout of UNC45A and HDDC3, and 57
CRISPRi mediated knockdown of HDDC3 (Fig. 1e, Extended Data Fig. 2d-g). In addition, 58
overexpression of MYNRL15 cDNAs failed to rescue the CRISPRi knockdown phenotype 59
(Fig. 1f). Altogether, these results indicate that neither of the flanking protein-coding genes, 60
nor the MYNRL15 transcript, is responsible for the function of this locus in myeloid leukemia 61
cells, and rather suggest MYNRL15 as an expressed noncoding regulatory locus. 62
Functional dissection of the MYNRL15 locus reveals crucial regulatory regions 63
Given the apparent dispensability of UNC45A, HDDC3, and the MYNRL15 transcript itself in 64
leukemia cells, we hypothesized that MYNRL15 may harbor DNA regulatory elements which 65
drive its leukemia dependency phenotype. To test this hypothesis, we functionally dissected 66
the MYNRL15 locus via complementary CRISPRi and CRISPR-Cas9 screens tiling a 15 kb 67
area centered on MYNRL15. Notably, this area exhibits features characteristic of cis-regula-68

Citations
More filters

Integrative Genomics Viewer

TL;DR: The sheer volume and scope of data posed by this flood of data pose a significant challenge to the development of efficient and intuitive visualization tools able to scale to very large data sets and to flexibly integrate multiple data types, including clinical data.

Activation of proto-oncogenes by disruption of chromosome neighborhoods

TL;DR: Wala et al. as discussed by the authors investigated whether proto-oncogenes occur within these structures and whether oncogene activation can occur via disruption of insulated neighborhood boundaries in cancer cells.
Journal ArticleDOI

Long noncoding RNAs as regulators of pediatric acute myeloid leukemia

TL;DR: In this paper , a mini review of the role of long noncoding RNAs (lncRNAs) in AML is presented, with a focus on pediatric AML.
References
More filters
Journal ArticleDOI

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

TL;DR: This work presents DESeq2, a method for differential analysis of count data, using shrinkage estimation for dispersions and fold changes to improve stability and interpretability of estimates, which enables a more quantitative analysis focused on the strength rather than the mere presence of differential expression.
Journal ArticleDOI

The Sequence Alignment/Map format and SAMtools

TL;DR: SAMtools as discussed by the authors implements various utilities for post-processing alignments in the SAM format, such as indexing, variant caller and alignment viewer, and thus provides universal tools for processing read alignments.
Journal ArticleDOI

Fast and accurate short read alignment with Burrows–Wheeler transform

TL;DR: Burrows-Wheeler Alignment tool (BWA) is implemented, a new read alignment package that is based on backward search with Burrows–Wheeler Transform (BWT), to efficiently align short sequencing reads against a large reference sequence such as the human genome, allowing mismatches and gaps.
Journal ArticleDOI

Trimmomatic: a flexible trimmer for Illumina sequence data

TL;DR: Timmomatic is developed as a more flexible and efficient preprocessing tool, which could correctly handle paired-end data and is shown to produce output that is at least competitive with, and in many cases superior to, that produced by other tools, in all scenarios tested.
Journal ArticleDOI

Fast gapped-read alignment with Bowtie 2

TL;DR: Bowtie 2 combines the strengths of the full-text minute index with the flexibility and speed of hardware-accelerated dynamic programming algorithms to achieve a combination of high speed, sensitivity and accuracy.
Related Papers (5)
Frequently Asked Questions (12)
Q1. What have the authors contributed in "Myeloid leukemia vulnerabilities at ctcf-enriched long noncoding rna loci" ?

The MYNRL15 in myeloid leukemia ( myeloid leukemia ) was published in this paper. 

Coding se-299  quences (CDS) from Ensembl v102 (release 11/2020) were used as inputs, and where pos-300  sible, sgRNAs were selected to target most, if not all, protein-coding isoforms. 

264 Lentiviral particles were produced by co-transfecting the expression vector and the packag-265  ing plasmids pMD2.G and psPAX2 (Addgene 12259 and 12260 respectively) into HEK293T 266  cells using polyethylenimine (PEI). 

For processing the 438  raw data, the authors used Trimmomatic94 to remove adapter sequences, followed by Kseq95 to trim 439  reads containing ≤6 bp of adapter sequence, which are not effectively handled by Trimmo-440  matic. 

Four days post-transduction, HSPCs were sorted 351  and plated in human methylcellulose complete medium HSC003 (R&D Systems) for colony-352  forming assays. 

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization. 

CTCF binding sites were determined 306  using ENCODE ChIP-seq peak calling data, and sgRNAs were selected to tile CTCF motifs 307  and/or point-source(s) within the peaks. 

476 For the analysis of C-LNCs in the context of AML, gene expression values were obtained 477  from the TCGA36 and NCI-TARGET35 AML patient cohorts. 

Three to nine 285  sgRNAs were selected per gene, depending on the number of different TSSs present in the 286  transcript isoforms and the distance between them. 

The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions. 

The sgRNA libraries used in this study were expressed from 261  the following backbones: SGL40C.EFS.dTomato (89395; CRISPRi lncRNA and MYNRL15 262  tiling), SGL.EFS.tBFP (173915; gained chromatin interaction region protein-coding), and 263  SGL.EFS.dTomato. 

qRT-PCR validation of WDR61 and IMP339  downregulation upon MYNRL15 perturbation using sgRNA C1.1. f, Retrieval of PAF1c lossassociated gene sets upon MYNRL15 perturbation in their RNA-seq data.