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.

Visualizing gene expression changes in time, space, and single cells with expressyouRcell.

The last decade has witnessed massive advancements in high-throughput techniques capable of producing increasingly complex gene expression datasets across time and space and at the resolution of single cells. Yet, the large volume of big data available and the complexity of experimental designs hamper an easy understanding and effective communication of the results. We present expressyouRcell, an easy-to-use R package to map the multi-dimensional variations of transcript and protein levels in dynamic cell pictographs. expressyouRcell visualizes gene expression variations as pictographic representations of cell-type thematic maps. expressyouRcell visually reduces the complexity of displaying gene expression and protein level changes across multiple measurements (time points or single-cell trajectories) by generating dynamic representations of cellular pictographs. We applied expressyouRcell to single cell, bulk RNA sequencing (RNA-seq), and proteomics datasets, demonstrating its flexibility and usability in the visualization of complex variations in gene expression. Our approach improves the standard quantitative interpretation and communication of relevant results.

TERRA stability is regulated by RALY and polyadenylation in a telomere-specific manner.

Telomeric repeat-containing RNA (TERRA) is a long non-coding RNA transcribed from telomeres that plays key roles in telomere maintenance. A fraction of TERRA is polyadenylated, and the presence of the poly(A) tail influences TERRA localization and stability. However, the mechanisms of TERRA biogenesis remain mostly elusive. Here, we show that the stability of TERRA transcripts is regulated by the RNA-binding protein associated with lethal yellow mutation (RALY). RALY depletion results in lower TERRA levels, impaired localization of TERRA at telomeres, and ultimately telomere damage. Importantly, we show that TERRA polyadenylation is telomere specific and that RALY preferentially stabilizes non-polyadenylated TERRA transcripts. Finally, we report that TERRA interacts with the poly(A)-binding protein nuclear 1 (PABPN1). Altogether, our results indicate that TERRA stability is regulated by the interplay between RALY and PABPN1, defined by the TERRA polyadenylation state. Our findings also suggest that different telomeres may trigger distinct TERRA-mediated responses.

NOC1 is a direct MYC target, and its protein interactome dissects its activity in controlling nucleolar function.

The nucleolus is a subnuclear compartment critical in ribosome biogenesis and cellular stress responses. These mechanisms are governed by a complex interplay of proteins, including NOC1, a member of the NOC family of nucleolar proteins responsible for controlling rRNA processing and ribosomal maturation. This study reveals a novel relationship between NOC1 and MYC transcription factor, known for its crucial role in controlling ribosomal biogenesis, cell growth, and proliferation. Here, we demonstrate that NOC1 functions as a direct target of MYC, as it is transcriptionally induced through a functional MYC-binding E-box sequence in the NOC1 promoter region. Furthermore, protein interactome analysis reveals that NOC1-complex includes the nucleolar proteins NOC2 and NOC3 and other nucleolar components such as Nucleostemin1 Ns1 transporters of ribosomal subunits and components involved in rRNA processing and maturation. In response to MYC, NOC1 expression and localization within the nucleolus significantly increase, suggesting a direct functional link between MYC activity and NOC1 function. Notably, NOC1 over-expression leads to the formation of large nuclear granules and enlarged nucleoli, which co-localize with nucleolar fibrillarin and Ns1. Additionally, we demonstrate that NOC1 expression is necessary for Ns1 nucleolar localization, suggesting a role for NOC1 in maintaining nucleolar structure. Finally, the co-expression of NOC1 and MYC enhances nucleolus size and maintains their co-localization, outlining another aspect of the cooperation between NOC1 and MYC in nucleolar dynamics. This study also reveals an enrichment with NOC1 with few proteins involved in RNA processing, modification, and splicing. Moreover, proteins such as Ythdc1, Flacc, and splenito are known to mediate N6-methyladenosine (m6A) methylation of mRNAs in nuclear export, revealing NOC1's potential involvement in coordinating RNA splicing and nuclear mRNA export. In summary, we uncovered novel roles for NOC1 in nucleolar homeostasis and established its direct connection with MYC in the network governing nucleolar structure and function. These findings also highlight NOC1's interaction with proteins relevant to specific RNA functions, suggesting a broader role in addition to its control of nucleolar homeostasis and providing new insight that can be further investigated.

Hydrogen peroxide induced by nerve injury promotes axon regeneration via connective tissue growth factor.

Regeneration of the neuromuscular junction (NMJ) leverages on extensive exchange of factors released from motor axon terminals (MATs), muscle fibers and perisynaptic Schwann cells (PSCs), among which hydrogen peroxide (H2O2) is a major pro-regenerative signal. To identify critical determinants of NMJ remodeling in response to injury, we performed temporal transcriptional profiling of NMJs from 2 month-old mice during MAT degeneration/regeneration, and cross-referenced the differentially expressed genes with those elicited by H2O2 in SCs. We identified an enrichment in extracellular matrix (ECM) transcripts, including Connective Tissue Growth Factor (Ctgf), which is usually expressed during development. We discovered that Ctgf levels are increased in a Yes-associated protein (YAP)-dependent fashion in response to rapid, local H2O2 signaling generated by stressed mitochondria in the injured sciatic nerve, a finding highlighting the importance of signals triggered by mechanical force to motor nerve repair. Through sequestration of Ctgf or inactivation of H2O2, we delayed the recovery of neuromuscular function by impairing SC migration and, in turn, axon-oriented re-growth. These data indicate that H2O2 and its downstream effector Ctgf are pro-regenerative factors that enable axonal growth, and reveal a striking ECM remodeling process during nerve regeneration upon local H2O2 signaling. Our study identifies key transcriptomic changes at the regenerating NMJ, providing a rich source of pro-regenerative factors with potential for alleviating the consequences of peripheral nerve injuries.

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.

Widespread association of the Argonaute protein AGO2 with meiotic chromatin suggests a distinct nuclear function in mammalian male reproduction.

Argonaute 2 (AGO2) is a ubiquitously expressed protein critical for regulation of mRNA translation and vital to animal development. AGO2 protein is found in both cytoplasmic and nuclear compartments, and although its cytoplasmic role is well studied, the biological relevance of nuclear AGO2 is unclear. Here, we address this problem in vivo using spermatogenic cells as a model. We find that AGO2 transiently binds both chromatin and nucleus-specific mRNA transcripts of hundreds of genes required for sperm production during male meiosis in mice, and that germline conditional knockout (cKO) of Ago2 causes depletion of the encoded proteins. Correspondingly, Ago2 cKO males show abnormal sperm head morphology and reduced sperm count, along with reduced postnatal viability of offspring. Together, our data reveal an unexpected nuclear role for AGO2 in enhancing expression of developmentally important genes during mammalian male reproduction.

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.

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.

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.

Deconvolution of in vivo protein-RNA contacts using fractionated eCLIP-seq.

Thousands of RNA-binding proteins orchestrate RNA processing and altered protein-RNA interactions frequently lead to disease. Here, we present experimental and computational analysis pipelines of fractionated eCLIP-seq (freCLIP-seq), a modification of enhanced UV-crosslinking and RNA immunoprecipitation followed by sequencing. FreCLIP-seq allows transcriptome-wide analysis of protein-RNA interactions at single-nucleotide level and provides an additional level of resolution by isolating binding signals of individual RNA-binding proteins within a multicomponent complex. Binding occupancy can be inferred from read counts and crosslinking events. For complete details on the use and execution of this protocol, please refer to Biancon et al. (2022).

Natural SINEUP RNAs in Autism Spectrum Disorders: RAB11B-AS1 Dysregulation in a Neuronal CHD8 Suppression Model Leads to RAB11B Protein Increase.

CHD8 represents one of the highest confidence genetic risk factors implied in Autism Spectrum Disorders, with most mutations leading to CHD8 haploinsufficiency and the insurgence of specific phenotypes, such as macrocephaly, facial dysmorphisms, intellectual disability, and gastrointestinal complaints. While extensive studies have been conducted on the possible consequences of CHD8 suppression and protein coding RNAs dysregulation during neuronal development, the effects of transcriptional changes of long non-coding RNAs (lncRNAs) remain unclear. In this study, we focused on a peculiar class of natural antisense lncRNAs, SINEUPs, that enhance translation of a target mRNA through the activity of two RNA domains, an embedded transposable element sequence and an antisense region. By looking at dysregulated transcripts following CHD8 knock down (KD), we first identified RAB11B-AS1 as a potential SINEUP RNA for its domain configuration. Then we demonstrated that such lncRNA is able to increase endogenous RAB11B protein amounts without affecting its transcriptional levels. RAB11B has a pivotal role in vesicular trafficking, and mutations on this gene correlate with intellectual disability and microcephaly. Thus, our study discloses an additional layer of molecular regulation which is altered by CHD8 suppression. This represents the first experimental confirmation that naturally occurring SINEUP could be involved in ASD pathogenesis and underscores the importance of dysregulation of functional lncRNAs in neurodevelopment.

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).

Active Ribosome Profiling with RiboLace: From Bench to Data Analysis.

Ribosome profiling is based on the deep sequencing of RNA fragments protected by ribosomes from nuclease digestion. This technique has been extensively used to study translation, with the unique ability to provide information about ribosomes positioning along transcripts at single-nucleotide resolution. Classical ribosome profiling approaches do not distinguish between fragments protected by either actively translating or inactive ribosomes. Here we describe an original method, called active ribosome profiling or RiboLace, which is based on a unique puromycin-containing molecule capable of isolating active ribosomes by means of an antibody-free and tag-free pull-down approach. This method allows reliable estimates of the translational state of any biological system, in high concordance with protein levels. RiboLace can be applied both in vitro and in vivo and generates snapshots of active ribosome footprints at single-nucleotide resolution and genome-wide level. RiboLace data are suitable for the analysis of translated genes, codon-specific translation rates, and local changes in ribosome occupancy profiles.