Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

RUNX1 isoform disequilibrium promotes the development of trisomy 21-associated myeloid leukemia.

Gain of chromosome 21 (Hsa21) is among the most frequent aneuploidies in leukemia. However, it remains unclear how partial or complete amplifications of Hsa21 promote leukemogenesis and why children with Down syndrome (DS) (ie, trisomy 21) are particularly at risk of leukemia development. Here, we propose that RUNX1 isoform disequilibrium with RUNX1A bias is key to DS-associated myeloid leukemia (ML-DS). Starting with Hsa21-focused CRISPR-CRISPR-associated protein 9 screens, we uncovered a strong and specific RUNX1 dependency in ML-DS cells. Expression of the RUNX1A isoform is elevated in patients with ML-DS, and mechanistic studies using murine ML-DS models and patient-derived xenografts revealed that excess RUNX1A synergizes with the pathognomonic Gata1s mutation during leukemogenesis by displacing RUNX1C from its endogenous binding sites and inducing oncogenic programs in complex with the MYC cofactor MAX. These effects were reversed by restoring the RUNX1A:RUNX1C equilibrium in patient-derived xenografts in vitro and in vivo. Moreover, pharmacological interference with MYC:MAX dimerization using MYCi361 exerted strong antileukemic effects. Thus, our study highlights the importance of alternative splicing in leukemogenesis, even on a background of aneuploidy, and paves the way for the development of specific and targeted therapies for ML-DS, as well as for other leukemias with Hsa21 aneuploidy or RUNX1 isoform disequilibrium.

The megakaryocytic transcription factor ARID3A suppresses leukemia pathogenesis.

Given the plasticity of hematopoietic stem and progenitor cells, multiple routes of differentiation must be blocked in the the pathogenesis of acute myeloid leukemia, the molecular basis of which is incompletely understood. We report that posttranscriptional repression of the transcription factor ARID3A by miR-125b is a key event in the pathogenesis of acute megakaryoblastic leukemia (AMKL). AMKL is frequently associated with trisomy 21 and GATA1 mutations (GATA1s), and children with Down syndrome are at a high risk of developing the disease. The results of our study showed that chromosome 21-encoded miR-125b synergizes with Gata1s to drive leukemogenesis in this context. Leveraging forward and reverse genetics, we uncovered Arid3a as the main miR-125b target behind this synergy. We demonstrated that, during normal hematopoiesis, this transcription factor promotes megakaryocytic differentiation in concert with GATA1 and mediates TGFß-induced apoptosis and cell cycle arrest in complex with SMAD2/3. Although Gata1s mutations perturb erythroid differentiation and induce hyperproliferation of megakaryocytic progenitors, intact ARID3A expression assures their megakaryocytic differentiation and growth restriction. Upon knockdown, these tumor suppressive functions are revoked, causing a blockade of dual megakaryocytic/erythroid differentiation and subsequently of AMKL. Inversely, restoring ARID3A expression relieves the arrest of megakaryocytic differentiation in AMKL patient-derived xenografts. This work illustrates how mutations in lineage-determining transcription factors and perturbation of posttranscriptional gene regulation can interact to block multiple routes of hematopoietic differentiation and cause leukemia. In AMKL, surmounting this differentiation blockade through restoration of the tumor suppressor ARID3A represents a promising strategy for treating this lethal pediatric disease.

AI-driven discovery of blood xenobiotic biomarkers in neovascular age-related macular degeneration using iterative random forests.

PURPOSE: To investigate the xenobiotic profiles of patients with neovascular age-related macular degeneration (nAMD) undergoing anti-vascular endothelial growth factor (anti-VEGF) intravitreal therapy (IVT) to identify biomarkers indicative of clinical phenotypes through advanced AI methodologies. METHODS: In this cross-sectional observational study, we analyzed 156 peripheral blood xenobiotic features in a cohort of 46 nAMD patients stratified by choroidal neovascularization (CNV) control under anti-VEGF IVT. We employed Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) for measurement and leveraged an AI-driven iterative Random Forests (iRF) approach for robust pattern recognition and feature selection, aligning molecular profiles with clinical phenotypes. RESULTS: AI-augmented iRF models effectively refined the metabolite spectrum by discarding non-predictive elements. Perfluorooctanesulfonate (PFOS) and Ethyl ß-glucopyranoside were identified as significant biomarkers through this process, associated with various clinically relevant phenotypes. Unlike single metabolite classes, drug metabolites were distinctly correlated with subretinal fluid presence. CONCLUSIONS: This study underscores the enhanced capability of AI, particularly iRF, in dissecting complex metabolomic data to elucidate the xenobiotic landscape of nAMD and environmental impact on the disease. The preliminary biomarkers discovered offer promising directions for personalized treatment strategies, although further validation in broader cohorts is essential for clinical application.

A combined computational and functional approach identifies IGF2BP2 as a driver of chemoresistance in a wide array of pre-clinical models of colorectal cancer.

AIM: Chemoresistance is a major cause of treatment failure in colorectal cancer (CRC) therapy. In this study, the impact of the IGF2BP family of RNA-binding proteins on CRC chemoresistance was investigated using in silico, in vitro, and in vivo approaches. METHODS: Gene expression data from a well-characterized cohort and publicly available cross-linking immunoprecipitation sequencing (CLIP-Seq) data were collected. Resistance to chemotherapeutics was assessed in patient-derived xenografts (PDXs) and patient-derived organoids (PDOs). Functional studies were performed in 2D and 3D cell culture models, including proliferation, spheroid growth, and mitochondrial respiration analyses. RESULTS: We identified IGF2BP2 as the most abundant IGF2BP in primary and metastastatic CRC, correlating with tumor stage in patient samples and tumor growth in PDXs. IGF2BP2 expression in primary tumor tissue was significantly associated with resistance to selumetinib, gefitinib, and regorafenib in PDOs and to 5-fluorouracil and oxaliplatin in PDX in vivo. IGF2BP2 knockout (KO) HCT116 cells were more susceptible to regorafenib in 2D and to oxaliplatin, selumitinib, and nintedanib in 3D cell culture. Further, a bioinformatic analysis using CLIP data suggested stabilization of target transcripts in primary and metastatic tumors. Measurement of oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) revealed a decreased basal OCR and an increase in glycolytic ATP production rate in IGF2BP2 KO. In addition, real-time reverse transcriptase polymerase chain reaction (qPCR) analysis confirmed decreased expression of genes of the respiratory chain complex I, complex IV, and the outer mitochondrial membrane in IGF2BP2 KO cells. CONCLUSIONS: IGF2BP2 correlates with CRC tumor growth in vivo and promotes chemoresistance by altering mitochondrial respiratory chain metabolism. As a druggable target, IGF2BP2 could be used in future CRC therapy to overcome CRC chemoresistance.

Systemic Blood Proteome Patterns Reflect Disease Phenotypes in Neovascular Age-Related Macular Degeneration.

There is early evidence of extraocular systemic signals effecting function and morphology in neovascular age-related macular degeneration (nAMD). The prospective, cross-sectional BIOMAC study is an explorative investigation of peripheral blood proteome profiles and matched clinical features to uncover systemic determinacy in nAMD under anti-vascular endothelial growth factor intravitreal therapy (anti-VEGF IVT). It includes 46 nAMD patients stratified by the level of disease control under ongoing anti-VEGF treatment. Proteomic profiles in peripheral blood samples of every patient were detected with LC-MS/MS mass spectrometry. The patients underwent extensive clinical examination with a focus on macular function and morphology. In silico analysis includes unbiased dimensionality reduction and clustering, a subsequent annotation of clinical features, and non-linear models for recognition of underlying patterns. The model assessment was performed using leave-one-out cross validation. The findings provide an exploratory demonstration of the link between systemic proteomic signals and macular disease pattern using and validating non-linear classification models. Three main results were obtained: (1) Proteome-based clustering identifies two distinct patient subclusters with the smaller one (n = 10) exhibiting a strong signature for oxidative stress response. Matching the relevant meta-features on the individual patient's level identifies pulmonary dysfunction as an underlying health condition in these patients. (2) We identify biomarkers for nAMD disease features with Aldolase C as a putative factor associated with superior disease control under ongoing anti-VEGF treatment. (3) Apart from this, isolated protein markers are only weakly correlated with nAMD disease expression. In contrast, applying a non-linear classification model identifies complex molecular patterns hidden in a high number of proteomic dimensions determining macular disease expression. In conclusion, so far unconsidered systemic signals in the peripheral blood proteome contribute to the clinically observed phenotype of nAMD, which should be examined in future translational research on AMD.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

Cell type-specific novel long non-coding RNA and circular RNA in the BLUEPRINT hematopoietic transcriptomes atlas.

Transcriptional profiling of hematopoietic cell subpopulations has helped to characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases. Previously, only the genes targeted by expression microarrays could be profiled genome-wide. High-throughput RNA sequencing, however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analyzed the BLUEPRINT consortium RNA-sequencing data for mature hematopoietic cell types. The data comprised 90 total RNA-sequencing samples, each composed of one of 27 cell types, and 32 small RNA-sequencing samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type-dependent expression. We also characterized the expression of circular RNA and found that these are also cell type-specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematologic development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

CD74 and CD44 Expression on CTCs in Cancer Patients with Brain Metastasis.

Up to 40% of advance lung, melanoma and breast cancer patients suffer from brain metastases (BM) with increasing incidence. Here, we assessed whether circulating tumor cells (CTCs) in peripheral blood can serve as a disease surrogate, focusing on CD44 and CD74 expression as prognostic markers for BM. We show that a size-based microfluidic approach in combination with a semi-automated cell recognition system are well suited for CTC detection in BM patients and allow further characterization of tumor cells potentially derived from BM. CTCs were found in 50% (7/14) of breast cancer, 50% (9/18) of non-small cell lung cancer (NSCLC) and 36% (4/11) of melanoma patients. The next-generation sequencing (NGS) analysis of nine single CTCs from one breast cancer patient revealed three different CNV profile groups as well as a resistance causing ERS1 mutation. CD44 and CD74 were expressed on most CTCs and their expression was strongly correlated, whereas matched breast cancer BM tissues were much less frequently expressing CD44 and CD74 (negative in 46% and 54%, respectively). Thus, plasticity of CD44 and CD74 expression during trafficking of CTCs in the circulation might be the result of adaptation strategies.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma.

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing.

Epigenetic alterations, that is, disruption of DNA methylation and chromatin architecture, are now acknowledged as a universal feature of tumorigenesis. Medulloblastoma, a clinically challenging, malignant childhood brain tumour, is no exception. Despite much progress from recent genomics studies, with recurrent changes identified in each of the four distinct tumour subgroups (WNT-pathway-activated, SHH-pathway-activated, and the less-well-characterized Group 3 and Group 4), many cases still lack an obvious genetic driver. Here we present whole-genome bisulphite-sequencing data from thirty-four human and five murine tumours plus eight human and three murine normal controls, augmented with matched whole-genome, RNA and chromatin immunoprecipitation sequencing data. This comprehensive data set allowed us to decipher several features underlying the interplay between the genome, epigenome and transcriptome, and its effects on medulloblastoma pathophysiology. Most notable were highly prevalent regions of hypomethylation correlating with increased gene expression, extending tens of kilobases downstream of transcription start sites. Focal regions of low methylation linked to transcription-factor-binding sites shed light on differential transcriptional networks between subgroups, whereas increased methylation due to re-normalization of repressed chromatin in DNA methylation valleys was positively correlated with gene expression. Large, partially methylated domains affecting up to one-third of the genome showed increased mutation rates and gene silencing in a subgroup-specific fashion. Epigenetic alterations also affected novel medulloblastoma candidate genes (for example, LIN28B), resulting in alternative promoter usage and/or differential messenger RNA/microRNA expression. Analysis of mouse medulloblastoma and precursor-cell methylation demonstrated a somatic origin for many alterations. Our data provide insights into the epigenetic regulation of transcription and genome organization in medulloblastoma pathogenesis, which are probably also of importance in a wider developmental and disease context.

Human Lineage-Specific Transcriptional Regulation through GA-Binding Protein Transcription Factor Alpha (GABPa).

A substantial fraction of phenotypic differences between closely related species are likely caused by differences in gene regulation. While this has already been postulated over 30 years ago, only few examples of evolutionary changes in gene regulation have been verified. Here, we identified and investigated binding sites of the transcription factor GA-binding protein alpha (GABPa) aiming to discover cis-regulatory adaptations on the human lineage. By performing chromatin immunoprecipitation-sequencing experiments in a human cell line, we found 11,619 putative GABPa binding sites. Through sequence comparisons of the human GABPa binding regions with orthologous sequences from 34 mammals, we identified substitutions that have resulted in 224 putative human-specific GABPa binding sites. To experimentally assess the transcriptional impact of those substitutions, we selected four promoters for promoter-reporter gene assays using human and African green monkey cells. We compared the activities of wild-type promoters to mutated forms, where we have introduced one or more substitutions to mimic the ancestral state devoid of the GABPa consensus binding sequence. Similarly, we introduced the human-specific substitutions into chimpanzee and macaque promoter backgrounds. Our results demonstrate that the identified substitutions are functional, both in human and nonhuman promoters. In addition, we performed GABPa knock-down experiments and found 1,215 genes as strong candidates for primary targets. Further analyses of our data sets link GABPa to cognitive disorders, diabetes, KRAB zinc finger (KRAB-ZNF), and human-specific genes. Thus, we propose that differences in GABPa binding sites played important roles in the evolution of human-specific phenotypes.

ARH-seq: identification of differential splicing in RNA-seq data.

The computational prediction of alternative splicing from high-throughput sequencing data is inherently difficult and necessitates robust statistical measures because the differential splicing signal is overlaid by influencing factors such as gene expression differences and simultaneous expression of multiple isoforms amongst others. In this work we describe ARH-seq, a discovery tool for differential splicing in case-control studies that is based on the information-theoretic concept of entropy. ARH-seq works on high-throughput sequencing data and is an extension of the ARH method that was originally developed for exon microarrays. We show that the method has inherent features, such as independence of transcript exon number and independence of differential expression, what makes it particularly suited for detecting alternative splicing events from sequencing data. In order to test and validate our workflow we challenged it with publicly available sequencing data derived from human tissues and conducted a comparison with eight alternative computational methods. In order to judge the performance of the different methods we constructed a benchmark data set of true positive splicing events across different tissues agglomerated from public databases and show that ARH-seq is an accurate, computationally fast and high-performing method for detecting differential splicing events.

The direction of cross affects [corrected] obesity after puberty in male but not female offspring.

BACKGROUND: We investigated parent-of-origin and allele-specific expression effects on obesity and hepatic gene expression in reciprocal crosses between the Berlin Fat Mouse Inbred line (BFMI) and C57Bl/6NCrl (B6N). RESULTS: We found that F1-males with a BFMI mother developed 1.8 times more fat mass on a high fat diet at 10 weeks than F1-males of a BFMI father. The phenotype was detectable from six weeks on and was preserved after cross-fostering. RNA-seq data of liver provided evidence for higher biosynthesis and elongation of fatty acids (p = 0.00635) in obese male offspring of a BFMI mother versus lean offspring of a BFMI father. Furthermore, fatty acid degradation (p = 0.00198) and the peroxisome pathway were impaired (p = 0.00094). The circadian rhythm was affected as well (p = 0.00087). Among the highest up-regulated protein coding genes in obese males were Acot4 (1.82 fold, p = 0.022), Cyp4a10 (1.35 fold, p = 0.026) and Cyp4a14 (1.32 fold, p = 0.012), which hydroxylize fatty acids and which are known to be increased in liver steatosis. Obese males showed lower expression of the genetically imprinted and paternally expressed 3 (Peg3) gene (0.31 fold, p = 0.046) and higher expression of the androgen receptor (Ar) gene (2.38 fold, p = 0.068). Allelic imbalance was found for expression of ATP-binding cassette transporter gene Abca8b. Several of the differentially expressed genes contain estrogen response elements. CONCLUSIONS: Parent-of-origin effects during gametogenesis and/or fetal development in an obese mother epigenetically modify the transcription of genes that lead to enhanced fatty acid synthesis and impair ß-oxidation in the liver of male, but not female F1 offspring. Down-regulation of Peg3 could contribute to trigger this metabolic setting. At puberty, higher amounts of the androgen receptor and altered access to estrogen response elements in affected genes are likely responsible for male specific expression of genes that were epigenetically triggered. A suggestive lack of estrogen binding motifs was found for highly down-regulated genes in adult hepatocytes of obese F1 males (p = 0.074).

Network and systems biology: essential steps in virtualising drug discovery and development.

The biological processes that keep us healthy or cause disease, as well as the mechanisms of action of possible drugs are inherently complex. In the face of this complexity, attempts at discovering new drugs to treat diseases have alternated between trial-and-error (typically on experimental systems) and grand simplification, usually based on much too little information. We now have the chance to combine these strategies through establishment of 'virtual patient' models, centred on a detailed molecular characterisation of thousands or even, in the future, millions of patients. In doing so, we lay the foundations for truly personalised therapy, as well as a far-reaching virtualisation of drug discovery and development in oncology and other areas of medicine.

Influence of RNA extraction methods and library selection schemes on RNA-seq data.

BACKGROUND: Gene expression analysis by RNA sequencing is now widely used in a number of applications surveying the whole transcriptomes of cells and tissues. The recent introduction of ribosomal RNA depletion protocols, such as RiboZero, has extended the view of the polyadenylated transcriptome to the poly(A)- fraction of the RNA. However, substantial amounts of intronic transcriptional activity has been reported in RiboZero protocols, raising issues regarding their potential nuclear origin and the impact on the actual sequence depth in exonic regions. RESULTS: Using HEK293 human cells as source material, we assessed here the impact of the two commonly used RNA extraction methods and of the library construction protocols (rRNA depletion versus mRNA) on 1) the relative abundance of intronic reads and 2) on the estimation of gene expression values. We benchmarked the rRNA depletion-based sequencing with a specific analysis of the cytoplasmic and nuclear transcriptome fractions, suggesting that the large majority of the intronic reads correspond to unprocessed nuclear transcripts rather than to independent transcriptional units. We show that Qiagen or TRIzol extraction methods retain differentially nuclear RNA species, and that consequently, rRNA depletion-based RNA sequencing protocols are particularly sensitive to the extraction methods. CONCLUSIONS: We could show that the combination of Trizol-based RNA extraction with rRNA depletion sequencing protocols led to the largest fraction of intronic reads, after the sequencing of the nuclear transcriptome. We discuss here the impact of the various strategies on gene expression and alternative splicing estimation measures. Further, we propose guidelines and a double selection strategy for minimizing the expression biases, without loss of information.

Exploring the Impact of Saccharin on Neovascular Age-Related Macular Degeneration: A Comprehensive Study in Patients and Mice.

Purpose: We aimed to determine the impact of artificial sweeteners (AS), especially saccharin, on the progression and treatment efficacy of patients with neovascular age-related macular degeneration (nAMD) under anti-vascular endothelial growth factor (anti-VEGF-A) treatment. Methods: In a cross-sectional study involving 46 patients with nAMD undergoing intravitreal anti-VEGF therapy, 6 AS metabolites were detected in peripheral blood using liquid chromatography - tandem mass spectrometry (LC-MS/MS). Disease features were statistically tested against these metabolite levels. Additionally, a murine choroidal neovascularization (CNV) model, induced by laser, was used to evaluate the effects of orally administered saccharin, assessing both imaging outcomes and gene expression patterns. Polymerase chain reaction (PCR) methods were used to evaluate functional expression of sweet taste receptors in a retinal pigment epithelium (RPE) cell line. Results: Saccharin levels in blood were significantly higher in patients with well-controlled CNV activity (P = 0.004) and those without subretinal hyper-reflective material (P = 0.015). In the murine model, saccharin-treated mice exhibited fewer leaking laser scars, lesser occurrence of bleeding, smaller fibrotic areas (P < 0.05), and a 40% decrease in mononuclear phagocyte accumulation (P = 0.06). Gene analysis indicated downregulation of inflammatory and VEGFR-1 response genes in the treated animals. Human RPE cells expressed taste receptor type 1 member 3 (TAS1R3) mRNA and reacted to saccharin stimulation with changes in mRNA expression. Conclusions: Saccharin appears to play a protective role in patients with nAMD undergoing intravitreal anti-VEGF treatment, aiding in better pathological lesion control and scar reduction. The murine study supports this observation, proposing saccharin's potential in mitigating pathological VEGFR-1-induced immune responses potentially via the RPE sensing saccharin in the blood stream.

Myeloid leukemia vulnerabilities embedded in long noncoding RNA locus MYNRL15.

The noncoding genome presents a largely untapped source of new biological insights, including thousands of long noncoding RNA (lncRNA) loci. While lncRNA dysregulation has been reported in myeloid malignancies, their functional relevance remains to be systematically interrogated. We performed CRISPRi screens of lncRNA signatures from normal and malignant hematopoietic cells and identified MYNRL15 as a myeloid leukemia dependency. Functional dissection suggests an RNA-independent mechanism mediated by two regulatory elements embedded in the locus. Genetic perturbation of these elements triggered a long-range chromatin interaction and downregulation of leukemia dependency genes near the gained interaction sites, as well as overall suppression of cancer dependency pathways. Thus, this study describes a new noncoding myeloid leukemia vulnerability and mechanistic concept for myeloid leukemia. Importantly, MYNRL15 perturbation caused strong and selective impairment of leukemia cells of various genetic backgrounds over normal hematopoietic stem and progenitor cells in vitro, and depletion of patient-derived xenografts in vivo.