Interpreting single-cell messages in normal and aberrant hematopoiesis with the Cell Marker Accordion.

Single-cell technologies offer a unique opportunity to explore cellular heterogeneity in hematopoiesis, reveal malignant hematopoietic cells with clinically significant features and measure gene signatures linked to pathological pathways. However, reliable identification of cell types is a crucial bottleneck in single-cell analysis. Available databases contain dissimilar nomenclature and non-concurrent marker sets, leading to inconsistent annotations and poor interpretability. Furthermore, current tools focus mostly on physiological cell types, lacking extensive applicability in disease. We developed the Cell Marker Accordion, a user-friendly platform for the automatic annotation and biological interpretation of single-cell populations based on consistency weighted markers. We validated our approach on peripheral blood and bone marrow single-cell datasets, using surface markers and expert-based annotation as the ground truth. In all cases, we significantly improved the accuracy in identifying cell types with respect to any single source database. Moreover, the Cell Marker Accordion can identify disease-critical cells and pathological processes, extracting potential biomarkers in a wide variety of contexts in human and murine single-cell datasets. It characterizes leukemia stem cell subtypes, including therapy-resistant cells in acute myeloid leukemia patients; it identifies malignant plasma cells in multiple myeloma samples; it dissects cell type alterations in splicing factor-mutant cells from myelodysplastic syndrome patients; it discovers activation of innate immunity pathways in bone marrow from mice treated with METTL3 inhibitors. The breadth of these applications elevates the Cell Marker Accordion as a flexible, faithful and standardized tool to annotate and interpret hematopoietic populations in single-cell datasets focused on the study of hematopoietic development and disease.

Impact of Memory T Cells on SARS-COV-2 Vaccine Response in Hematopoietic Stem Cell Transplant.

During the COVID-19 pandemic, hematopoietic stem cell transplant (HSCT) recipients faced an elevated mortality rate from SARS-CoV-2 infection, ranging between 10-40%. The SARS-CoV-2 mRNA vaccines are important tools in preventing severe disease, yet their efficacy in the post-transplant setting remains unclear, especially in patients subjected to myeloablative chemotherapy and immunosuppression. We evaluated the humoral and adaptive immune responses to the SARS-CoV-2 mRNA vaccination series in 42 HSCT recipients and 5 healthy controls. Peripheral blood mononuclear nuclear cells and serum were prospectively collected before and after each dose of the SARS-CoV-2 vaccine. Post-vaccination responses were assessed by measuring anti-spike IgG and nucleocapsid titers, and antigen specific T cell activity, before and after vaccination. In order to examine mechanisms behind a lack of response, pre-and post-vaccine samples were selected based on humoral and cellular responses for single-cell RNA sequencing with TCR and BCR sequencing. Our observations revealed that while all participants eventually mounted a humoral response, transplant recipients had defects in memory T cell populations that were associated with an absence of T cell response, some of which could be detected pre-vaccination.

ALKBH5 modulates hematopoietic stem and progenitor cell energy metabolism through m6A modification-mediated RNA stability control.

N6-methyladenosine (m6A) RNA modification controls numerous cellular processes. To what extent these post-transcriptional regulatory mechanisms play a role in hematopoiesis has not been fully elucidated. We here show that the m6A demethylase alkB homolog 5 (ALKBH5) controls mitochondrial ATP production and modulates hematopoietic stem and progenitor cell (HSPC) fitness in an m6A-dependent manner. Loss of ALKBH5 results in increased RNA methylation and instability of oxoglutarate-dehydrogenase (Ogdh) messenger RNA and reduction of OGDH protein levels. Limited OGDH availability slows the tricarboxylic acid (TCA) cycle with accumulation of α-ketoglutarate (α-KG) and conversion of α-KG into L-2-hydroxyglutarate (L-2-HG). L-2-HG inhibits energy production in both murine and human hematopoietic cells in vitro. Impaired mitochondrial energy production confers competitive disadvantage to HSPCs and limits clonogenicity of Mll-AF9-induced leukemia. Our study uncovers a mechanism whereby the RNA m6A demethylase ALKBH5 regulates the stability of metabolic enzyme transcripts, thereby controlling energy metabolism in hematopoiesis and leukemia.

The Akt/mTOR and MNK/eIF4E pathways rewire the prostate cancer translatome to secrete HGF, SPP1 and BGN and recruit suppressive myeloid cells.

Cancer is highly infiltrated by myeloid-derived suppressor cells (MDSCs). Currently available immunotherapies do not completely eradicate MDSCs. Through a genome-wide analysis of the translatome of prostate cancers driven by different genetic alterations, we demonstrate that prostate cancer rewires its secretome at the translational level to recruit MDSCs. Among different secreted proteins released by prostate tumor cells, we identified Hgf, Spp1 and Bgn as the key factors that regulate MDSC migration. Mechanistically, we found that the coordinated loss of Pdcd4 and activation of the MNK/eIF4E pathways regulate the mRNAs translation of Hgf, Spp1 and Bgn. MDSC infiltration and tumor growth were dampened in prostate cancer treated with the MNK1/2 inhibitor eFT508 and/or the AKT inhibitor ipatasertib, either alone or in combination with a clinically available MDSC-targeting immunotherapy. This work provides a therapeutic strategy that combines translation inhibition with available immunotherapies to restore immune surveillance in prostate cancer.

Alpha-1 Adrenergic Antagonists Sensitize Neuroblastoma to Therapeutic Differentiation.

Neuroblastoma (NB) is an aggressive childhood tumor, with high-risk cases having a 5-year overall survival probability of approximately 50%. The multimodal therapeutic approach for NB includes treatment with the retinoid isotretinoin (13-cis retinoic acid; 13cRA), which is used in the post-consolidation phase as an antiproliferation and prodifferentiation agent to minimize residual disease and prevent relapse. Through small-molecule screening, we identified isorhamnetin (ISR) as a synergistic compound with 13cRA in inhibiting up to 80% of NB cell viability. The synergistic effect was accompanied by a marked increase in the expression of the adrenergic receptor α1B (ADRA1B) gene. Genetic knockout of ADRA1B or its specific blockade using α1/α1B adrenergic antagonists led to selective sensitization of MYCN-amplified NB cells to cell viability reduction and neural differentiation induced by 13cRA, thus mimicking ISR activity. Administration of doxazosin, a safe α1-antagonist used in pediatric patients, in combination with 13cRA in NB xenografted mice exerted marked control of tumor growth, whereas each drug alone was ineffective. Overall, this study identified the α1B adrenergic receptor as a pharmacologic target in NB, supporting the evaluation of adding α1-antagonists to the post-consolidation therapy of NB to more efficiently control residual disease. SIGNIFICANCE: Targeting α-adrenergic receptors synergizes with isotretinoin to suppress growth and to promote differentiation of neuroblastoma, revealing a combinatorial approach for more effective management of the disease and prevention of relapse.

Recruited monocytes/macrophages drive pulmonary neutrophilic inflammation and irreversible lung tissue remodeling in cystic fibrosis.

Persistent neutrophil-dominated lung inflammation contributes to lung damage in cystic fibrosis (CF). However, the mechanisms that drive persistent lung neutrophilia and tissue deterioration in CF are not well characterized. Starting from the observation that, in patients with CF, c-c motif chemokine receptor 2 (CCR2)+ monocytes/macrophages are abundant in the lungs, we investigate the interplay between monocytes/macrophages and neutrophils in perpetuating lung tissue damage in CF. Here we show that CCR2+ monocytes in murine CF lungs drive pathogenic transforming growth factor ß (TGF-ß) signaling and sustain a pro-inflammatory environment by facilitating neutrophil recruitment. Targeting CCR2 to lower the numbers of monocytes in CF lungs ameliorates neutrophil inflammation and pathogenic TGF-ß signaling and prevents lung tissue damage. This study identifies CCR2+ monocytes as a neglected contributor to the pathogenesis of CF lung disease and as a therapeutic target for patients with CF, for whom lung hyperinflammation and tissue damage remain an issue despite recent advances in CF transmembrane conductance regulator (CFTR)-specific therapeutic agents.

A slow-cycling/quiescent cells subpopulation is involved in glioma invasiveness.

Pediatric and adult high-grade gliomas are the most common primary malignant brain tumors, with poor prognosis due to recurrence and tumor infiltration after therapy. Quiescent cells have been implicated in tumor recurrence and treatment resistance, but their direct visualization and targeting remain challenging, precluding their mechanistic study. Here, we identify a population of malignant cells expressing Prominin-1 in a non-proliferating state in pediatric high-grade glioma patients. Using a genetic tool to visualize and ablate quiescent cells in mouse brain cancer and human cancer organoids, we reveal their localization at both the core and the edge of the tumors, and we demonstrate that quiescent cells are involved in infiltration of brain cancer cells. Finally, we find that Harmine, a DYRK1A/B inhibitor, partially decreases the number of quiescent and infiltrating cancer cells. Our data point to a subpopulation of quiescent cells as partially responsible of tumor invasiveness, one of the major causes of brain cancer morbidity.

Pre-mRNA splicing factor U2AF2 recognizes distinct conformations of nucleotide variants at the center of the pre-mRNA splice site signal.

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies.

Precision analysis of mutant U2AF1 activity reveals deployment of stress granules in myeloid malignancies.

Splicing factor mutations are common among cancers, recently emerging as drivers of myeloid malignancies. U2AF1 carries hotspot mutations in its RNA-binding motifs; however, how they affect splicing and promote cancer remain unclear. The U2AF1/U2AF2 heterodimer is critical for 3' splice site (3'SS) definition. To specifically unmask changes in U2AF1 function in vivo, we developed a crosslinking and immunoprecipitation procedure that detects contacts between U2AF1 and the 3'SS AG at single-nucleotide resolution. Our data reveal that the U2AF1 S34F and Q157R mutants establish new 3'SS contacts at -3 and +1 nucleotides, respectively. These effects compromise U2AF2-RNA interactions, resulting predominantly in intron retention and exon exclusion. Integrating RNA binding, splicing, and turnover data, we predicted that U2AF1 mutations directly affect stress granule components, which was corroborated by single-cell RNA-seq. Remarkably, U2AF1-mutant cell lines and patient-derived MDS/AML blasts displayed a heightened stress granule response, pointing to a novel role for biomolecular condensates in adaptive oncogenic strategies.

In vivo anti-tumor effect of PARP inhibition in IDH1/2 mutant MDS/AML resistant to targeted inhibitors of mutant IDH1/2.

Treatment options for patients with relapsed/refractory acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS) are scarce. Recurring mutations, such as mutations in isocitrate dehydrogenase-1 and -2 (IDH1/2) are found in subsets of AML and MDS, are therapeutically targeted by mutant enzyme-specific small molecule inhibitors (IDHmi). IDH mutations induce diverse metabolic and epigenetic changes that drive malignant transformation. IDHmi alone are not curative and resistance commonly develops, underscoring the importance of alternate therapeutic options. We were first to report that IDH1/2 mutations induce a homologous recombination (HR) defect, which confers sensitivity to poly (ADP)-ribose polymerase inhibitors (PARPi). Here, we show that the PARPi olaparib is effective against primary patient-derived IDH1/2-mutant AML/ MDS xeno-grafts (PDXs). Olaparib efficiently reduced overall engraftment and leukemia-initiating cell frequency as evident in serial transplantation assays in IDH1/2-mutant but not -wildtype AML/MDS PDXs. Importantly, we show that olaparib is effective in both IDHmi-naïve and -resistant AML PDXs, critical given the high relapse and refractoriness rates to IDHmi. Our pre-clinical studies provide a strong rationale for the translation of PARP inhibition to patients with IDH1/2-mutant AML/ MDS, providing an additional line of therapy for patients who do not respond to or relapse after targeted mutant IDH inhibition.

Influence of the microenvironment on modulation of the host response by typhoid toxin.

Bacterial genotoxins cause DNA damage in eukaryotic cells, resulting in activation of the DNA damage response (DDR) in vitro. These toxins are produced by Gram-negative bacteria, enriched in the microbiota of inflammatory bowel disease (IBD) and colorectal cancer (CRC) patients. However, their role in infection remains poorly characterized. We address the role of typhoid toxin in modulation of the host-microbial interaction in health and disease. Infection with a genotoxigenic Salmonella protects mice from intestinal inflammation. We show that the presence of an active genotoxin promotes DNA fragmentation and senescence in vivo, which is uncoupled from an inflammatory response and unexpectedly associated with induction of an anti-inflammatory environment. The anti-inflammatory response is lost when infection occurs in mice with acute colitis. These data highlight a complex context-dependent crosstalk between bacterial-genotoxin-induced DDR and the host immune response, underlining an unexpected role for bacterial genotoxins.

Transcriptome-wide quantification of double-stranded RNAs in live mouse tissues by dsRIP-Seq.

Double-stranded RNAs (dsRNAs) are abundantly present in cells, playing multiple regulatory functions. dsRNAs of viral origin activate innate immune responses. Since RNA editing and modifications affect the structure and recognition of RNAs, their alteration can result in the accumulation of aberrant endogenous dsRNAs inducing a deleterious innate immune response. Here, we present a complete protocol for the measurement of dsRNAs in a live mouse tissue using dsRNA immunoprecipitation and sequencing (dsRIP-Seq). This protocol focuses on tissue isolation, dsRNA immunoprecipitation and downstream computational analysis. For complete details on the use and execution of this protocol, please refer to Gao et al. (2020).

Combined liver-cytokine humanization comes to the rescue of circulating human red blood cells.

In vivo models that recapitulate human erythropoiesis with persistence of circulating red blood cells (RBCs) have remained elusive. We report an immunodeficient murine model in which combined human liver and cytokine humanization confer enhanced human erythropoiesis and RBC survival in the circulation. We deleted the fumarylacetoacetate hydrolase (Fah) gene in MISTRG mice expressing several human cytokines in place of their murine counterparts. Liver humanization by intrasplenic injection of human hepatocytes (huHep) eliminated murine complement C3 and reduced murine Kupffer cell density. Engraftment of human sickle cell disease (SCD)-derived hematopoietic stem cells in huHepMISTRGFah -/- mice resulted in vaso-occlusion that replicated acute SCD pathology. Combined liver-cytokine-humanized mice will facilitate the study of diseases afflicting RBCs, including bone marrow failure, hemoglobinopathies, and malaria, and also preclinical testing of therapies.

High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue.

We present deterministic barcoding in tissue for spatial omics sequencing (DBiT-seq) for co-mapping of mRNAs and proteins in a formaldehyde-fixed tissue slide via next-generation sequencing (NGS). Parallel microfluidic channels were used to deliver DNA barcodes to the surface of a tissue slide, and crossflow of two sets of barcodes, A1-50 and B1-50, followed by ligation in situ, yielded a 2D mosaic of tissue pixels, each containing a unique full barcode AB. Application to mouse embryos revealed major tissue types in early organogenesis as well as fine features like microvasculature in a brain and pigmented epithelium in an eye field. Gene expression profiles in 10-µm pixels conformed into the clusters of single-cell transcriptomes, allowing for rapid identification of cell types and spatial distributions. DBiT-seq can be adopted by researchers with no experience in microfluidics and may find applications in a range of fields including developmental biology, cancer biology, neuroscience, and clinical pathology.

The altered transcriptome of pediatric myelodysplastic syndrome revealed by RNA sequencing.

Pediatric myelodysplastic syndrome (PMDS) is a very rare and still poorly characterized disorder. In this work, we identified novel potential targets of PMDS by determining genes with aberrant expression, which can be correlated with PMDS pathogenesis. We identified 291 differentially expressed genes (DEGs) in PMDS patients, comprising genes involved in the regulation of apoptosis and the cell cycle, ribosome biogenesis, inflammation and adaptive immunity. Ten selected DEGs were then validated, confirming the sequencing data. These DEGs will potentially represent new molecular biomarkers and therapeutic targets for PMDS.

Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

We hereby provide the initial portrait of lincNORS, a spliced lincRNA generated by the MIR193BHG locus, entirely distinct from the previously described miR-193b-365a tandem. While inducible by low O2 in a variety of cells and associated with hypoxia in vivo, our studies show that lincNORS is subject to multiple regulatory inputs, including estrogen signals. Biochemically, this lincRNA fine-tunes cellular sterol/steroid biosynthesis by repressing the expression of multiple pathway components. Mechanistically, the function of lincNORS requires the presence of RALY, an RNA-binding protein recently found to be implicated in cholesterol homeostasis. We also noticed the proximity between this locus and naturally occurring genetic variations highly significant for sterol/steroid-related phenotypes, in particular the age of sexual maturation. An integrative analysis of these variants provided a more formal link between these phenotypes and lincNORS, further strengthening the case for its biological relevance.

m6A Modification Prevents Formation of Endogenous Double-Stranded RNAs and Deleterious Innate Immune Responses during Hematopoietic Development.

N6-methyladenosine (m6A) is the most abundant RNA modification, but little is known about its role in mammalian hematopoietic development. Here, we show that conditional deletion of the m6A writer METTL3 in murine fetal liver resulted in hematopoietic failure and perinatal lethality. Loss of METTL3 and m6A activated an aberrant innate immune response, mediated by the formation of endogenous double-stranded RNAs (dsRNAs). The aberrantly formed dsRNAs were long, highly m6A modified in their native state, characterized by low folding energies, and predominantly protein coding. We identified coinciding activation of pattern recognition receptor pathways normally tasked with the detection of foreign dsRNAs. Disruption of the aberrant immune response via abrogation of downstream Mavs or Rnasel signaling partially rescued the observed hematopoietic defects in METTL3-deficient cells in vitro and in vivo. Our results suggest that m6A modification protects against endogenous dsRNA formation and a deleterious innate immune response during mammalian hematopoietic development.

Low iron promotes megakaryocytic commitment of megakaryocytic-erythroid progenitors in humans and mice.

The mechanisms underlying thrombocytosis in patients with iron deficiency anemia remain unknown. Here, we present findings that support the hypothesis that low iron biases the commitment of megakaryocytic (Mk)-erythroid progenitors (MEPs) toward the Mk lineage in both human and mouse. In MEPs of transmembrane serine protease 6 knockout (Tmprss6-/-) mice, which exhibit iron deficiency anemia and thrombocytosis, we observed a Mk bias, decreased labile iron, and decreased proliferation relative to wild-type (WT) MEPs. Bone marrow transplantation assays suggest that systemic iron deficiency, rather than a local role for Tmprss6-/- in hematopoietic cells, contributes to the MEP lineage commitment bias observed in Tmprss6-/- mice. Nontransgenic mice with acquired iron deficiency anemia also show thrombocytosis and Mk-biased MEPs. Gene expression analysis reveals that messenger RNAs encoding genes involved in metabolic, vascular endothelial growth factor, and extracellular signal-regulated kinase (ERK) pathways are enriched in Tmprss6-/- vs WT MEPs. Corroborating our findings from the murine models of iron deficiency anemia, primary human MEPs exhibit decreased proliferation and Mk-biased commitment after knockdown of transferrin receptor 2, a putative iron sensor. Signal transduction analyses reveal that both human and murine MEP have lower levels of phospho-ERK1/2 in iron-deficient conditions compared with controls. These data are consistent with a model in which low iron in the marrow environment affects MEP metabolism, attenuates ERK signaling, slows proliferation, and biases MEPs toward Mk lineage commitment.

A highly efficient and faithful MDS patient-derived xenotransplantation model for pre-clinical studies.

Comprehensive preclinical studies of Myelodysplastic Syndromes (MDS) have been elusive due to limited ability of MDS stem cells to engraft current immunodeficient murine hosts. Here we report a MDS patient-derived xenotransplantation model in cytokine-humanized immunodeficient "MISTRG" mice that provides efficient and faithful disease representation across all MDS subtypes. MISTRG MDS patient-derived xenografts (PDX) reproduce patients' dysplastic morphology with multi-lineage representation, including erythro- and megakaryopoiesis. MISTRG MDS-PDX replicate the original sample's genetic complexity and can be propagated via serial transplantation. MISTRG MDS-PDX demonstrate the cytotoxic and differentiation potential of targeted therapeutics providing superior readouts of drug mechanism of action and therapeutic efficacy. Physiologic humanization of the hematopoietic stem cell niche proves critical to MDS stem cell propagation and function in vivo. The MISTRG MDS-PDX model opens novel avenues of research and long-awaited opportunities in MDS research.

HuD Is a Neural Translation Enhancer Acting on mTORC1-Responsive Genes and Counteracted by the Y3 Small Non-coding RNA.

The RNA-binding protein HuD promotes neurogenesis and favors recovery from peripheral axon injury. HuD interacts with many mRNAs, altering both stability and translation efficiency. We generated a nucleotide resolution map of the HuD RNA interactome in motor neuron-like cells, identifying HuD target sites in 1,304 mRNAs, almost exclusively in the 3' UTR. HuD binds many mRNAs encoding mTORC1-responsive ribosomal proteins and translation factors. Altered HuD expression correlates with the translation efficiency of these mRNAs and overall protein synthesis, in a mTORC1-independent fashion. The predominant HuD target is the abundant, small non-coding RNA Y3, amounting to 70% of the HuD interaction signal. Y3 functions as a molecular sponge for HuD, dynamically limiting its recruitment to polysomes and its activity as a translation and neuron differentiation enhancer. These findings uncover an alternative route to the mTORC1 pathway for translational control in motor neurons that is tunable by a small non-coding RNA.