Transcription factor regulation of ribosomal RNA in hematopoiesis

PURPOSE OF REVIEW

Ribosomal RNAs (rRNAs) are transcribed within nucleoli from rDNA repeats by RNA Polymerase I (Pol I). There is variation in rRNA transcription rates across the hematopoietic tree, and leukemic blast cells have prominent nucleoli, indicating abundant ribosome biogenesis. The mechanisms underlying these variations are poorly understood. The purpose of this review is to summarize findings of rDNA binding and Pol I regulation by hematopoietic transcription factors.

RECENT FINDINGS

Our group recently used custom genome assemblies optimized for human and mouse rDNA mapping to map nearly 2200 ChIP-Seq datasets for nearly 250 factors to rDNA, allowing us to identify conserved occupancy patterns for multiple transcription factors. We confirmed known rDNA occupancy of MYC and RUNX factors, and identified new binding sites for CEBP factors, IRF factors, and SPI1 at canonical motif sequences. We also showed that CEBPA degradation rapidly leads to reduced Pol I occupancy and nascent rRNA in mouse myeloid cells.

SUMMARY

We propose that a number of hematopoietic transcription factors bind rDNA and potentially regulate rRNA transcription. Our model has implications for normal and malignant hematopoiesis. This review summarizes the literature, and outlines experimental considerations to bear in mind while dissecting transcription factor roles on rDNA.

PHF6 suppresses self-renewal of leukemic stem cells in AML

Acute myeloid leukemia is characterized by uncontrolled proliferation of self-renewing myeloid progenitors. PHF6 is a chromatin-binding protein mutated in myeloid leukemias, and its loss increases mouse HSC self-renewal without malignant transformation. We report here that Phf6 knockout increases the aggressiveness of Hoxa9-driven AML over serial transplantation, and increases the frequency of leukemia initiating cells. We define the in vivo hierarchy of Hoxa9-driven AML and identify a population that we term the 'LIC-e' (leukemia initiating cells enriched) population. We find that Phf6 loss has context-specific transcriptional effects, skewing the LIC-e transcriptome to a more stem-like state. We demonstrate that LIC-e accumulation in Phf6 knockout AML occurs not due to effects on cell cycle or apoptosis, but due to an increase in the fraction of its progeny that retain LIC-e identity. Overall, our work indicates that Phf6 loss increases AML self-renewal through context-specific effects on leukemia stem cells.

FISH-Flow to quantify nascent and mature ribosomal RNA in mouse and human cells

FISH-Flow (fluorescence in situ hybridization-flow cytometry) involves hybridizing fluorescent oligos to RNA and quantifying fluorescence at a single-cell level using flow cytometry. Here, we present a FISH-Flow protocol to quantify nascent 47S and mature 18S and 28S rRNAs in mouse and human cells, including rRNA quantification across cell cycle stages using DNA staining. We describe steps for cell preparation, hybridization of fluorescent probes against rRNA, and DNA staining. We then detail procedures for flow cytometry and data analysis. For complete details on the use and execution of this protocol, please refer to Antony et al. (2022).1.

Real-world effectiveness of intensive chemotherapy with 7&3 versus venetoclax and hypomethylating agent in acute myeloid leukemia

Intensive chemotherapy with cytarabine and anthracycline (7&3) remains the standard therapy for patients medically fit for induction, but the assessment of fitness remains controversial. Venetoclax and hypomethylating agent (ven/HMA) combination therapy has improved outcomes in unfit patients but no prospective study has assessed ven/HMA versus 7&3 as initial therapy in older, fit patients. Given no studies and expectation of ven/HMA use in patients outside of trial criteria, we evaluated retrospective outcomes among newly diagnosed patients. A nationwide electronic health record (EHR)-derived database and the University of Pennsylvania EHR identified 312 patients receiving 7&3 and 488 receiving ven/HMA who were 60-75 years old without history of organ failure. Ven/HMA patients were older and more likely to have secondary AML, adverse cytogenetics, and adverse mutations. Median overall survival (OS) for patients receiving intensive chemotherapy was 22 versus 10 months for ven/HMA (HR 0.53, 95% CI 0.40-0.60). Controlling for measured baseline characteristic imbalances reduced survival advantage by half (HR 0.71, 95% CI 0.53-0.94). A sub-group of patients with equipoise, likelihood at least 30%-70% of receiving either treatment, had similar OS outcomes (HR 1.10, 95% CI 0.75-1.6). Regarding safety outcomes, 60-day mortality was higher for ven/HMA (15% vs. 6% at 60 days) despite higher documented infections and febrile neutropenia for 7&3. In this multicenter real-word dataset, patients selected for intensive chemotherapy had superior OS but a large group had similar outcomes with ven/HMA. Prospective randomized studies, controlling for both measured and unmeasured confounders, must confirm this outcome.

Construction and validation of customized genomes for human and mouse ribosomal DNA mapping

rRNAs are transcribed from ribosomal DNA (rDNA) repeats, the most intensively transcribed loci in the genome. Due to their repetitive nature, there is a lack of genome assemblies suitable for rDNA mapping, creating a vacuum in our understanding of how the most abundant RNA in the cell is regulated. Our recent work revealed binding of numerous mammalian transcription and chromatin factors to rDNA. Several of these factors were known to play critical roles in development, tissue function, and malignancy, but their potential roles in rDNA regulation remained unexplored. This demonstrated the blind spot into which rDNA has fallen in genetic and epigenetic studies and highlighted an unmet need for public rDNA-optimized genome assemblies. Here, we customized five human and mouse assemblies-hg19 (GRCh37), hg38 (GRCh38), hs1 (T2T-CHM13), mm10 (GRCm38), and mm39 (GRCm39)-to render them suitable for rDNA mapping. The standard builds of these genomes contain numerous fragmented or repetitive rDNA loci. We identified and masked all rDNA-like regions, added a single rDNA reference sequence of the appropriate species as a ∼45 kb chromosome designated "chromosome R," and created annotation files to aid visualization of rDNA features in browser tracks. We validated these customized genomes for mapping of known rDNA-binding proteins and present a simple workflow for mapping chromatin immunoprecipitation-sequencing datasets. Customized genome assemblies, annotation files, positive and negative control tracks, and Snapgene files of standard rDNA reference sequences have been deposited to GitHub. These resources make rDNA mapping and visualization more readily accessible to a broad audience.

Rapid-kinetics degron benchmarking reveals off-target activities and mixed agonism-antagonism of MYB inhibitors

Attenuating aberrant transcriptional circuits holds great promise for the treatment of numerous diseases, including cancer. However, development of transcriptional inhibitors is hampered by the lack of a generally accepted functional cellular readout to characterize their target specificity and on-target activity. We benchmarked the direct gene-regulatory signatures of six agents reported as inhibitors of the oncogenic transcription factor MYB against targeted MYB degradation in a nascent transcriptomics assay. The inhibitors demonstrated partial specificity for MYB target genes but displayed significant off-target activity. Unexpectedly, the inhibitors displayed bimodal on-target effects, acting as mixed agonists-antagonists. Our data uncover unforeseen agonist effects of small molecules originally developed as TF inhibitors and argue that rapid-kinetics benchmarking against degron models should be used for functional characterization of transcriptional modulators.

Leukemia core transcriptional circuitry is a sparsely interconnected hierarchy stabilized by incoherent feed-forward loops

Lineage-defining transcription factors form densely interconnected circuits in chromatin occupancy assays, but the functional significance of these networks remains underexplored. We reconstructed the functional topology of a leukemia cell transcription network from the direct gene-regulatory programs of eight core transcriptional regulators established in pre-steady state assays coupling targeted protein degradation with nascent transcriptomics. The core regulators displayed narrow, largely non-overlapping direct transcriptional programs, forming a sparsely interconnected functional hierarchy stabilized by incoherent feed-forward loops. BET bromodomain and CDK7 inhibitors disrupted the core regulators' direct programs, acting as mixed agonists/antagonists. The network is predictive of dynamic gene expression behaviors in time-resolved assays and clinically relevant pathway activity in patient populations.

Control of ribosomal RNA synthesis by hematopoietic transcription factors

Ribosomal RNAs (rRNAs) are the most abundant cellular RNAs, and their synthesis from rDNA repeats by RNA polymerase I accounts for the bulk of all transcription. Despite substantial variation in rRNA transcription rates across cell types, little is known about cell-type-specific factors that bind rDNA and regulate rRNA transcription to meet tissue-specific needs. Using hematopoiesis as a model system, we mapped about 2,200 ChIP-seq datasets for 250 transcription factors (TFs) and chromatin proteins to human and mouse rDNA and identified robust binding of multiple TF families to canonical TF motifs on rDNA. Using a 47S-FISH-Flow assay developed for nascent rRNA quantification, we demonstrated that targeted degradation of C/EBP alpha (CEBPA), a critical hematopoietic TF with conserved rDNA binding, caused rapid reduction in rRNA transcription due to reduced RNA Pol I occupancy. Our work identifies numerous potential rRNA regulators and provides a template for dissection of TF roles in rRNA transcription.

Transcriptional Plasticity Drives Leukemia Immune Escape

Relapse of acute myeloid leukemia (AML) after allogeneic bone marrow transplantation has been linked to immune evasion due to reduced expression of major histocompatibility complex class II (MHCII) genes through unknown mechanisms. In this work, we developed CORENODE, a computational algorithm for genome-wide transcription network decomposition that identified a transcription factor (TF) tetrad consisting of IRF8, MYB, MEF2C, and MEIS1, regulating MHCII expression in AML cells. We show that reduced MHCII expression at relapse is transcriptionally driven by combinatorial changes in the expression of these TFs, where MYB and IRF8 play major opposing roles, acting independently of the IFNγ/CIITA pathway. Beyond the MHCII genes, MYB and IRF8 antagonistically regulate a broad genetic program responsible for cytokine signaling and T-cell stimulation that displays reduced expression at relapse. A small number of cells with altered TF abundance and silenced MHCII expression are present at the time of initial leukemia diagnosis, likely contributing to eventual relapse.

SIGNIFICANCE

Our findings point to an adaptive transcriptional mechanism of AML evolution after allogeneic transplantation whereby combinatorial fluctuations of TF expression under immune pressure result in the selection of cells with a silenced T-cell stimulation program. This article is highlighted in the In This Issue feature, p. 369.

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia

In this study, Harada et al. identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. This study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

HectD1 controls hematopoietic stem cell regeneration by coordinating ribosome assembly and protein synthesis

Impaired ribosome function is the underlying etiology in a group of bone marrow failure syndromes called ribosomopathies. However, how ribosomes are regulated remains poorly understood, as are approaches to restore hematopoietic stem cell (HSC) function loss because of defective ribosome biogenesis. Here we reveal a role of the E3 ubiquitin ligase HectD1 in regulating HSC function via ribosome assembly and protein translation. Hectd1-deficient HSCs exhibit a striking defect in transplantation ability and ex vivo maintenance concomitant with reduced protein synthesis and growth rate under stress conditions. Mechanistically, HectD1 ubiquitinates and degrades ZNF622, an assembly factor for the ribosomal 60S subunit. Hectd1 loss leads to accumulation of ZNF622 and the anti-association factor eIF6 on 60S, resulting in 60S/40S joining defects. Importantly, Znf622 depletion in Hectd1-deficient HSCs restored ribosomal subunit joining, protein synthesis, and HSC reconstitution capacity. These findings highlight the importance of ubiquitin-coordinated ribosome assembly in HSC regeneration.

YODEL: Peak calling software for HITS-CLIP data

YODEL is a peak calling software for analyzing RNA sequencing data generated by High-Throughput Sequencing of RNA isolated by Crosslinking Immunoprecipitation (HITS-CLIP; also known as CLIP-SEQ), a method to identify RNA-protein interactions genome-wide. We designed YODEL to analyze HITS-CLIP experiments, in which Argonaute proteins are immunoprecipitated, followed by sequencing of the associated RNA in order to identify bound microRNAs and their mRNA targets. The HITS-CLIP sequenced reads are mapped to the genome, and then read peaks are visualized where clustered sets of reads map to the same region. Several peak calling algorithms have been developed to define the boundaries of these peaks. In contrast to other peak callers for HITS-CLIP data, such as Piranha, YODEL does not map the starts of reads to fixed interval bins, but instead uses a heuristic approach to iteratively find the tallest point within a set clustered reads and examine bases upstream and downstream of that point until a peak has been determined. This allows the peak boundary to be defined more precisely than coordinates that are multiples of the bin size. Per-sample peak counts are also generated by YODEL, which quickly enables downstream differential representation analysis. YODEL is available athttps://github.com/LancePalmerStJude/YODEL/.

Unlinking an lncRNA from Its Associated cis Element

Long non-coding (lnc) RNAs can regulate gene expression and protein functions. However, the proportion of lncRNAs with biological activities among the thousands expressed in mammalian cells is controversial. We studied Lockd (lncRNA downstream of Cdkn1b), a 434-nt polyadenylated lncRNA originating 4 kb 3' to the Cdkn1b gene. Deletion of the 25-kb Lockd locus reduced Cdkn1b transcription by approximately 70% in an erythroid cell line. In contrast, homozygous insertion of a polyadenylation cassette 80 bp downstream of the Lockd transcription start site reduced the entire lncRNA transcript level by >90% with no effect on Cdkn1b transcription. The Lockd promoter contains a DNase-hypersensitive site, binds numerous transcription factors, and physically associates with the Cdkn1b promoter in chromosomal conformation capture studies. Therefore, the Lockd gene positively regulates Cdkn1b transcription through an enhancer-like cis element, whereas the lncRNA itself is dispensable, which may be the case for other lncRNAs.

A new 'Linc' between noncoding RNAs and blood development

Long noncoding RNAs (lncRNAs) are a large and diverse class of functional RNAs that regulate important biological processes, including cell division, survival, and differentiation. In this issue of Genes & Development, Hu and colleagues (2573-2578) report the discovery of LincRNA erythroid prosurvival (LincRNA-EPS), a murine lncRNA that facilitates red blood cell formation (erythropoiesis) by suppressing apoptosis. This finding expands the repertoire of lncRNA functions and illustrates a novel genetic pathway that potentially can be exploited for treating anemias.

Severe eosinophilia as a manifestation of acute graft-versus-host disease

Significant peripheral eosinophilia in association with acute graft-versus-host disease (GVHD) is rare. Here we report a case of eosinophilia in a 30-year-old woman with relapsed acute myelogenous leukemia after an allogeneic bone marrow transplant who was treated with donor lymphocyte infusion (DLI). A month after the DLI, she developed a pruritic erythematous rash, hepatic enzyme elevations and eosinophilia that peaked at 23,300 cells/mm(3). A diagnosis of acute GVHD was made and the patient was treated with corticosteroids with a resolution of all of the aforementioned findings. The authors suggest that eosinophilia in a bone marrow transplantation patient should raise suspicion of GVHD.