Multiple myeloma (MM) is an incurable disease, whose clinical heterogeneity makes its management challenging, highlighting the need for biological features to guide improved therapies. Deregulation of specific long non-coding RNAs (lncRNAs) has been shown in MM, nevertheless, the complete lncRNA transcriptome has not yet been elucidated. In this work, we identified 40,511 novel lncRNAs in MM samples. lncRNAs accounted for 82% of the MM transcriptome and were more heterogeneously expressed than coding genes. A total of 10,351 overexpressed and 9,535 downregulated lncRNAs were identified in MM patients when compared with normal bone-marrow plasma cells. Transcriptional dynamics study of lncRNAs in the context of normal B-cell maturation revealed 989 lncRNAs with exclusive expression in MM, among which 89 showed de novo epigenomic activation. Knockdown studies on one of these lncRNAs, SMILO (specific myeloma intergenic long non-coding RNA), resulted in reduced proliferation and induction of apoptosis of MM cells, and activation of the interferon pathway. We also showed that the expression of lncRNAs, together with clinical and genetic risk alterations, stratified MM patients into several progression-free survival and overall survival groups. In summary, our global analysis of the lncRNAs transcriptome reveals the presence of specific lncRNAs associated with the biological and clinical behavior of the disease.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in older individuals. Here we analyze the transcriptome of human HSCs purified from young and older healthy adults, as well as MDS patients, identifying transcriptional alterations following different patterns of expression. While aging-associated lesions seem to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detect the upregulation of the transcription factor DNA Damage Inducible Transcript 3 (DDIT3). Overexpression of DDIT3 in human healthy HSCs induces an MDS-like transcriptional state, and dyserythropoiesis, an effect associated with a failure in the activation of transcriptional programs required for normal erythroid differentiation. Moreover, DDIT3 knockdown in CD34⁺ cells from MDS patients with anemia is able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythroid differentiation characterizing MDS patients.

Hematologic toxicity is a common side effect of chimeric antigen receptor T-cell (CAR-T) therapies, being particularly severe among patients with relapsed or refractory multiple myeloma (MM). In this study, we characterized 48 patients treated with B-cell maturation antigen (BCMA) CAR-T cells to understand kinetics of cytopenia, identify predictive factors, and determine potential mechanisms underlying these toxicities. We observed that overall incidence of cytopenia was 95.7%, and grade >3 thrombocytopenia and neutropenia, 1 month after infusion, was observed in 57% and 53% of the patients, respectively, being still present after 1 year in 4 and 3 patients, respectively. Baseline cytopenia and high peak inflammatory markers were highly correlated with cytopenia that persisted up to 3 months. To determine potential mechanisms underlying cytopenias, we evaluated the paracrine effect of BCMA CAR-T cells on hematopoietic stem and progenitor cell (HSPC) differentiation using an ex vivo myeloid differentiation model. Phenotypic analysis showed that supernatants from activated CAR-T cells (spCAR) halted HSPC differentiation, promoting more immature phenotypes, which could be prevented with a combination of interferon γ, tumor necrosis factor α/β, transforming growth factor β, interleukin-6 (IL-6) and IL-17 inhibitors. Single-cell RNA sequencing demonstrated upregulation of transcription factors associated with early stages of hematopoietic differentiation in the presence of spCAR (GATA2, RUNX1, CEBPA) and a decrease in the activity of key regulons involved in neutrophil and monocytic maturation (ID2 and MAFB). These results suggest that CAR-T activation induces HSPC maturation arrest through paracrine effects and provides potential treatments to mitigate the severity of this toxicity.

The histone H3 at lysine 27 (H3K27) demethylase lysine demethylase 6A (KDM6A) is a tumor suppressor in multiple cancers, including multiple myeloma (MM). We created isogenic MM cells disrupted for KDM6A and tagged the endogenous protein to facilitate genome-wide studies. KDM6A binds genes associated with immune recognition and cytokine signaling. Most importantly, KDM6A binds and activates NLRC5 and CIITA, which encode regulators of major histocompatibility complex genes. Patient data indicate that NLRC5 and CIITA are downregulated in MM with low KDM6A expression. Chromatin analysis shows that KDM6A binds poised and active enhancers and KDM6A loss led to decreased H3K27ac at enhancers, increased H3K27me3 levels in body of genes bound by KDM6A, and decreased gene expression. Reestablishing histone acetylation with an HDAC3 inhibitor leads to upregulation of major histocompatibility complex expression, offering a strategy to restore immunogenicity of KDM6A-deficient tumors. Loss of Kdm6a in Kirsten rat sarcoma virus (K-RAS)-transformed murine fibroblasts led to increased growth in vivo associated with decreased T-cell infiltration.

Single-cell RNA-Sequencing has the potential to provide deep biological insights by revealing complex regulatory interactions across diverse cell phenotypes at single-cell resolution. However, current single-cell gene regulatory network inference methods produce a single regulatory network per input dataset, limiting their capability to uncover complex regulatory relationships across related cell phenotypes. We present SimiC, a single-cell gene regulatory inference framework that overcomes this limitation by jointly inferring distinct, but related, gene regulatory dynamics per phenotype. We show that SimiC uncovers key regulatory dynamics missed by previously proposed methods across a range of systems, both model and non-model alike. In particular, SimiC was able to uncover CAR T cell dynamics after tumor recognition and key regulatory patterns on a regenerating liver, and was able to implicate glial cells in the generation of distinct behavioral states in honeybees. SimiC hence establishes a new approach to quantitating regulatory architectures between distinct cellular phenotypes, with far-reaching implications for systems biology.

While myelodysplastic syndromes with del(5q) (del(5q) MDS) comprises a well-defined hematological subgroup, the molecular basis underlying its origin remains unknown. Using single cell RNA-seq (scRNA-seq) on CD34⁺ progenitors from del(5q) MDS patients, we have identified cells harboring the deletion, characterizing the transcriptional impact of this genetic insult on disease pathogenesis and treatment response. Interestingly, both del(5q) and non-del(5q) cells present similar transcriptional lesions, indicating that all cells, and not only those harboring the deletion, may contribute to aberrant hematopoietic differentiation. However, gene regulatory network (GRN) analyses reveal a group of regulons showing aberrant activity that could trigger altered hematopoiesis exclusively in del(5q) cells, pointing to a more prominent role of these cells in disease phenotype. In del(5q) MDS patients achieving hematological response upon lenalidomide treatment, the drug reverts several transcriptional alterations in both del(5q) and non-del(5q) cells, but other lesions remain, which may be responsible for potential future relapses. Moreover, lack of hematological response is associated with the inability of lenalidomide to reverse transcriptional alterations. Collectively, this study reveals transcriptional alterations that could contribute to the pathogenesis and treatment response of del(5q) MDS.