Immunosuppressive Activity of Exosomes from Granulocytic Myeloid-Derived Suppressor Cells in a Murine Model of Immune Bone Marrow Failure.

We previously reported that granulocytic myeloid-derived suppressor cells (G-MDSCs) suppressed T-cell activation and attenuated bone marrow failure (BMF) in a minor histocompatibility (minor-H) antigen mismatched murine aplastic anemia (AA) model. In the current study, we tested the hypothesis that exosomes, a subset of extracellular vesicles, are responsible at least partially for G-MDSCs' therapeutic efficacy. Indeed, exosomes isolated from GMDSCs (G-MDSC-exos) suppressed CD4+ and CD8+ T-cell proliferation in vitro and mildly attenuated immune BMF in the minor-H mismatched AA model. G-MDSC-exos treatment significantly increased red blood cells, hemoglobin, and total bone marrow (BM) cells, and moderately reduced BM CD8+ T cells. G-MDSC-exos' effects were associated with upregulations in an array of lymphocyte-suppression-related miRNAs such as hsa-miR-142-5p, miR-19a-3p, and miR-19b-3p in both BM CD4+ and CD8+ T cells. We concluded that G-MDSC-exos attenuate immune BMF via modulating the delivery of immunosuppressive miRNAs into activated T lymphocytes.

Early activation of inflammatory pathways in UBA1-mutated hematopoietic stem and progenitor cells in VEXAS.

VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome is a pleiotropic, severe autoinflammatory disease caused by somatic mutations in the ubiquitin-like modifier activating enzyme 1 (UBA1) gene. To elucidate VEXAS pathophysiology, we performed transcriptome sequencing of single bone marrow mononuclear cells and hematopoietic stem and progenitor cells (HSPCs) from VEXAS patients. HSPCs are biased toward myeloid (granulocytic) differentiation, and against lymphoid differentiation in VEXAS. Activation of multiple inflammatory pathways (interferons and tumor necrosis factor alpha) occurs ontogenically early in primitive hematopoietic cells and particularly in the myeloid lineage in VEXAS, and inflammation is prominent in UBA1-mutated cells. Dysregulation in protein degradation likely leads to higher stress response in VEXAS HSPCs, which positively correlates with inflammation. TCR usage is restricted and there are increased cytotoxicity and IFN-γ signaling in T cells. In VEXAS syndrome, both aberrant inflammation and myeloid predominance appear intrinsic to hematopoietic stem cells mutated in UBA1.

Single-cell genomics in acquired bone marrow failure syndromes.

Mechanistic studies of immune bone marrow failure are difficult because of the scarcity of residual cells, the involvement of multiple cell types, and the inherent complexities of hematopoiesis and immunity. Single-cell genomic technologies and bioinformatics allow extensive, multidimensional analysis of a very limited number of cells. We review emerging applications of single-cell techniques, and early results related to disease pathogenesis: effector and target cell populations and relationships, cell-autonomous and nonautonomous phenotypes in clonal hematopoiesis, transcript splicing, chromosomal abnormalities, and T-cell receptor usage and clonality. Dense and complex data from single-cell techniques provide insights into pathophysiology, natural history, and therapeutic drug effects.

Conditional survival and standardized mortality ratios of patients with severe aplastic anemia surviving at least one year after hematopoietic cell transplantation or immunosuppressive therapy.

Immunosuppressive treatment (IST) and hematopoietic cell transplant (HCT) are standard therapies for severe aplastic anemia (SAA). We report on conditional survival and standardized mortality ratios (SMR), which compare the mortality risk with the general population adjusted for age, gender, and race/ethnicity, in patients with SAA alive for at least 12 months after treatment with IST or HCT between 2000 and 2018. Given changes to treatment regimens and differences in length of follow-up, two treatment periods were defined a priori: 2000-2010 and 2011-2018. The SMR of patients treated during the period 2000-2010 and who survived one year were 3.50 (95% confidence interval [CI]: 2.62-4.58), 4.12 (95% CI: 3.20-5.21), and 8.62 (95% CI: 6.88-10.67) after IST, matched related donor HCT, and alternative donor HCT, respectively. For the period 2011-2018, the corresponding SMR were 2.89 (95% CI: 1.54-4.94), 3.12 (95% CI: 1.90-4.82), and 4.75 (95% CI: 3.45-6.38), respectively. For IST patients, their mortality risk decreased over time, and became comparable to the general population by five years. For patients who underwent HCT during 2000-2010 and 2011-2018, their mortality risk became comparable to the general population after ten years and after five years, respectively. Thus, 1-year survivors after IST or HCT can expect their longevity beyond five years to be comparable to that of the general US population.

Spectrum of clonal hematopoiesis in VEXAS syndrome.

Vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic (VEXAS) syndrome is caused by somatic mutations in UBA1 (UBA1mut) and characterized by heterogenous systemic autoinflammation and progressive hematologic manifestations, meeting criteria for myelodysplastic syndrome (MDS) and plasma cell dyscrasias. The landscape of myeloid-related gene mutations leading to typical clonal hematopoiesis (CH) in these patients is unknown. Retrospectively, we screened 80 patients with VEXAS for CH in their peripheral blood (PB) and correlated the findings with clinical outcomes in 77 of them. UBA1mut were most common at hot spot p.M41 (median variant allele frequency [VAF] = 75%). Typical CH mutations cooccurred with UBA1mut in 60% of patients, mostly in DNMT3A and TET2, and were not associated with inflammatory or hematologic manifestations. In prospective single-cell proteogenomic sequencing (scDNA), UBA1mut was the dominant clone, present mostly in branched clonal trajectories. Based on integrated bulk and scDNA analyses, clonality in VEXAS followed 2 major patterns: with either typical CH preceding UBA1mut selection in a clone (pattern 1) or occurring as an UBA1mut subclone or in independent clones (pattern 2). VAF in the PB differed markedly between DNMT3A and TET2 clones (median VAF of 25% vs 1%). DNMT3A and TET2 clones associated with hierarchies representing patterns 1 and 2, respectively. Overall survival for all patients was 60% at 10 years. Transfusion-dependent anemia, moderate thrombocytopenia, and typical CH mutations, each correlated with poor outcome. In VEXAS, UBA1mut cells are the primary cause of systemic inflammation and marrow failure, being a new molecularly defined somatic entity associated with MDS. VEXAS-associated MDS is distinct from classical MDS in its presentation and clinical course.

Alemtuzumab in relapsed immune severe aplastic anemia: Long-term results of a phase II study.

Immune severe aplastic anemia (SAA) is characterized by pancytopenia and immune-mediated bone marrow destruction. SAA may be treated with hematopoietic stem cell transplantation (HSCT) or immunosuppressive therapy (IST). However, 30% of patients treated with IST relapse. We previously reported a clinical trial of alemtuzumab in which more than half of 25 relapsed SAA patients (56%) responded hematologically. Here, we present long-term results of a total of 42 patients. Participants with SAA who had previously completed antithymocyte globulin (ATG)-based IST, but had relapsed, were enrolled on this study. Alemtuzumab was administered intravenously (IV) (n = 28) or subcutaneously (SC) (n = 14). The primary endpoint was hematologic response at 6 months. Secondary endpoints included relapse, clonal evolution, and survival. This trial was registered at clinicaltrials.gov (NCT00195624). Patients were enrolled over 9 years, with median follow-up of 6 years. Median age was 32 years, with 57% being female. At 6 months, 18 patients (43%) achieved response; 15 (54%) of those who received IV compared with 3 (21%) who received SC therapy. Six patients (14%) had durable long-term response without need for subsequent AA-directed therapy or HSCT at last follow-up. Nine patients had clonal evolution, with high-risk evolution occurring in 6. Overall survival was 67% at median follow-up of 6 years. Prolonged iatrogenic immunosuppression was observed as long as 2 years after alemtuzumab administration. Alemtuzumab induces responses in relapsed SAA, some of which are durable long-term. However, immunosuppression can persist for years, requiring long-term monitoring.

Differential diagnosis of bone marrow failure syndromes guided by machine learning.

The choice to postpone treatment while awaiting genetic testing can result in significant delay in definitive therapies in patients with severe pancytopenia. Conversely, the misdiagnosis of inherited bone marrow failure (BMF) can expose patients to ineffectual and expensive therapies, toxic transplant conditioning regimens, and inappropriate use of an affected family member as a stem cell donor. To predict the likelihood of patients having acquired or inherited BMF, we developed a 2-step data-driven machine-learning model using 25 clinical and laboratory variables typically recorded at the initial clinical encounter. For model development, patients were labeled as having acquired or inherited BMF depending on their genomic data. Data sets were unbiasedly clustered, and an ensemble model was trained with cases from the largest cluster of a training cohort (n = 359) and validated with an independent cohort (n = 127). Cluster A, the largest group, was mostly immune or inherited aplastic anemia, whereas cluster B comprised underrepresented BMF phenotypes and was not included in the next step of data modeling because of a small sample size. The ensemble cluster A-specific model was accurate (89%) to predict BMF etiology, correctly predicting inherited and likely immune BMF in 79% and 92% of cases, respectively. Our model represents a practical guide for BMF diagnosis and highlights the importance of clinical and laboratory variables in the initial evaluation, particularly telomere length. Our tool can be potentially used by general hematologists and health care providers not specialized in BMF, and in under-resourced centers, to prioritize patients for genetic testing or for expeditious treatment.

Granulocytic myeloid-derived suppressor cells to prevent and treat murine immune-mediated bone marrow failure.

Myeloid-derived suppressor cells (MDSCs) are immature myeloid cells that originate in the bone marrow (BM) and have immunoregulatory functions. MDSCs have been implicated in the pathogenesis of several autoimmune diseases but have not been investigated in immune aplastic anemia (AA). We examined the roles of granulocytic-MDSCs (G-MDSCs) in murine models of human AA and BM failure (BMF). As both prophylaxis and therapy, BM-derived G-MDSCs improved pancytopenia and BM cellularity and suppressed BM T-cell infiltration in major histocompatibility complex (MHC)-matched C.B10 BMF mice. These effects were not obtained in the MHC-mismatched CByB6F1 AA model, likely because of MHC disparity between G-MDSCs and donor T cells. Single-cell RNA sequencing demonstrated that G-MDSCs downregulated cell cycle-related genes in BM-infiltrated T cells, consistent with suppression of T-cell proliferation by G-MDSCs through reactive oxygen species pathways. Clearance of G-MDSCs in the MHC-mismatched CByB6F1 model using anti-Ly6G antibody facilitated T cell-mediated BM destruction, suggesting an intrinsic immunosuppressive property of G-MDSCs. However, the same anti-Ly6G antibody in the MHC-matched C.B10 AA model mildly mitigated BMF, associated with expansion of an intermediate Ly6G population. Our results demonstrate that G-MDSC eradication and therapeutic efficacy are immune context-dependent.

Effects of nandrolone decanoate on telomere length and clinical outcome in patients with telomeropathies: a prospective trial.

Androgens have been reported to elongate telomeres in retrospective and prospective trials with patients with telomeropathies, mainly with bone marrow failure. In our single-arm prospective clinical trial (clinicaltrials gov. Identifier: NCT02055456), 17 patients with short telomeres and/or germline pathogenic variants in telomere biology genes associated with at least one cytopenia and/or radiologic diagnosis of interstitial lung disease were treated with 5 mg/kg of intramuscular nandrolone decanoate every 15 days for 2 years. Ten of 13 evaluable patients (77%) showed telomere elongation at 12 months by flow-fluorescence in situ hybridization (average increase, 0.87 kb; 95% confidence interval: 0.20-1.55 kb; P=0.01). At 24 months, all ten evaluable patients showed telomere elongation (average increase, 0.49 kb; 95% confidence interval: 0.24-1.23 kb; P=0.18). Hematologic response was achieved in eight of 16 patients (50%) with marrow failure at 12 months, and in ten of 16 patients (63%) at 24 months. Seven patients had interstitial lung disease at baseline, and two and three had pulmonary response at 12 and 24 months, respectively. Two patients died due to pulmonary failure during treatment. In the remaining evaluable patients, the pulmonary function remained stable or improved, but showed consistent decline after cessation of treatment. Somatic mutations in myeloid neoplasm-related genes were present in a minority of patients and were mostly stable during drug treatment. The most common adverse events were elevations in liver function test levels in 88%, acne in 59%, and virilization in 59%. No adverse events grade ≥4 was observed. Our findings indicate that nandrolone decanoate elongates telomeres in patients with telomeropathies, which correlated with clinical improvement in some cases and tolerable adverse events.

Efficacy of JAK1/2 inhibition in murine immune bone marrow failure.

Immune aplastic anemia (AA) is a severe blood disease characterized by T-lymphocyte- mediated stem cell destruction. Hematopoietic stem cell transplantation and immunosuppression are effective, but they entail costs and risks, and are not always successful. The Janus kinase (JAK) 1/2 inhibitor ruxolitinib (RUX) suppresses cytotoxic T-cell activation and inhibits cytokine production in models of graft-versus-host disease. We tested RUX in murine immune AA for potential therapeutic benefit. After infusion of lymph node (LN) cells mismatched at the major histocompatibility complex [C67BL/6 (B6)⇒CByB6F1], RUX, administered as a food additive (Rux-chow), attenuated bone marrow hypoplasia, ameliorated peripheral blood pancytopenia, preserved hematopoietic progenitors, and prevented mortality, when used either prophylactically or therapeutically. RUX suppressed the infiltration, proliferation, and activation of effector T cells in the bone marrow and mitigated Fas-mediated apoptotic destruction of target hematopoietic cells. Similar effects were obtained when Rux-chow was fed to C.B10 mice in a minor histocompatibility antigen mismatched (B6⇒C.B10) AA model. RUX only modestly suppressed lymphoid and erythroid hematopoiesis in normal and irradiated CByB6F1 mice. Our data support clinical trials of JAK/STAT inhibitors in human AA and other immune bone marrow failure syndromes.

Time-Varying Gene Expression Network Analysis Reveals Conserved Transition States in Hematopoietic Differentiation between Human and Mouse.

(1) Background: analyses of gene networks can elucidate hematopoietic differentiation from single-cell gene expression data, but most algorithms generate only a single, static network. Because gene interactions change over time, it is biologically meaningful to examine time-varying structures and to capture dynamic, even transient states, and cell-cell relationships. (2) Methods: a transcriptomic atlas of hematopoietic stem and progenitor cells was used for network analysis. After pseudo-time ordering with Monocle 2, LOGGLE was used to infer time-varying networks and to explore changes of differentiation gene networks over time. A range of network analysis tools were used to examine properties and genes in the inferred networks. (3) Results: shared characteristics of attributes during the evolution of differentiation gene networks showed a "U" shape of network density over time for all three branches for human and mouse. Differentiation appeared as a continuous process, originating from stem cells, through a brief transition state marked by fewer gene interactions, before stabilizing in a progenitor state. Human and mouse shared hub genes in evolutionary networks. (4) Conclusions: the conservation of network dynamics in the hematopoietic systems of mouse and human was reflected by shared hub genes and network topological changes during differentiation.

CellCallEXT: Analysis of Ligand-Receptor and Transcription Factor Activities in Cell-Cell Communication of Tumor Immune Microenvironment.

(1) Background: Single-cell RNA sequencing (scRNA-seq) data are useful for decoding cell-cell communication. CellCall is a tool that is used to infer inter- and intracellular communication pathways by integrating paired ligand-receptor (L-R) and transcription factor (TF) activities from steady-state data and thus cannot directly handle two-condition comparisons. For tumor and healthy status, it can only individually analyze cells from tumor or healthy tissue and examine L-R pairs only identified in either tumor or healthy controls, but not both together. Furthermore, CellCall is highly affected by gene expression specificity in tissues. (2) Methods: CellCallEXT is an extension of CellCall that deconvolutes intercellular communication and related internal regulatory signals based on scRNA-seq. Information on Reactome was retrieved and integrated with prior knowledge of L-R-TF signaling and gene regulation datasets of CellCall. (3) Results: CellCallEXT was successfully applied to examine tumors and immune cell microenvironments and to identify the altered L-R pairs and downstream gene regulatory networks among immune cells. Application of CellCallEXT to scRNA-seq data from patients with deficiency of adenosine deaminase 2 demonstrated its ability to impute dysfunctional intercellular communication and related transcriptional factor activities. (4) Conclusions: CellCallEXT provides a practical tool to examine intercellular communication in disease based on scRNA-seq data.

Anti-COX-2 autoantibody is a novel biomarker of immune aplastic anemia.

In immune aplastic anemia (IAA), severe pancytopenia results from the immune-mediated destruction of hematopoietic stem cells. Several autoantibodies have been reported, but no clinically applicable autoantibody tests are available for IAA. We screened autoantibodies using a microarray containing >9000 proteins and validated the findings in a large international cohort of IAA patients (n = 405) and controls (n = 815). We identified a novel autoantibody that binds to the C-terminal end of cyclooxygenase 2 (COX-2, aCOX-2 Ab). In total, 37% of all adult IAA patients tested positive for aCOX-2 Ab, while only 1.7% of the controls were aCOX-2 Ab positive. Sporadic non-IAA aCOX-2 Ab positive cases were observed among patients with related bone marrow failure diseases, multiple sclerosis, and type I diabetes, whereas no aCOX-2 Ab seropositivity was detected in the healthy controls, in patients with non-autoinflammatory diseases or rheumatoid arthritis. In IAA, anti-COX-2 Ab positivity correlated with age and the HLA-DRB1*15:01 genotype. 83% of the >40 years old IAA patients with HLA-DRB1*15:01 were anti-COX-2 Ab positive, indicating an excellent sensitivity in this group. aCOX-2 Ab positive IAA patients also presented lower platelet counts. Our results suggest that aCOX-2 Ab defines a distinct subgroup of IAA and may serve as a valuable disease biomarker.

Predictors of clonal evolution and myeloid neoplasia following immunosuppressive therapy in severe aplastic anemia.

Predictors, genetic characteristics, and long-term outcomes of patients with SAA who clonally evolved after immunosuppressive therapy (IST) were assessed. SAA patients were treated with IST from 1989-2020. Clonal evolution was categorized as "high-risk" (overt myeloid neoplasm [meeting WHO criteria for dysplasia, MPN or acute leukemia] or isolated chromosome-7 abnormality/complex karyotype without dysplasia or overt myeloid neoplasia) or "low-risk" (non-7 or non-complex chromosome abnormalities without morphological evidence of dysplasia or myeloid neoplasia). Univariate and multivariate analysis using Fine-Gray competing risk regression model determined predictors. Long-term outcomes included relapse, overall survival (OS) and hematopoietic stem cell transplant (HSCT). Somatic mutations in myeloid cancer genes were assessed in evolvers and in 407 patients 6 months after IST. Of 663 SAA patients, 95 developed clonal evolution. Pre-treatment age >48 years and ANC > 0.87 × 109/L were strong predictors of high-risk evolution. OS was 37% in high-risk clonal evolution by 5 years compared to 94% in low-risk. High-risk patients who underwent HSCT had improved OS. Eltrombopag did not increase high-risk evolution. Splicing factors and RUNX1 somatic variants were detected exclusively at high-risk evolution; DNMT3A, BCOR/L1 and ASXL1 were present in both. RUNX1, splicing factors and ASXL1 somatic mutations detected at 6 months after IST predicted high-risk evolution.

Single-cell RNA sequencing coupled to TCR profiling of large granular lymphocyte leukemia T cells.

T-cell large granular lymphocyte leukemia (T-LGLL) is a lymphoproliferative disease and bone marrow failure syndrome which responds to immunosuppressive therapies. We show single-cell TCR coupled with RNA sequencing of CD3+ T cells from 13 patients, sampled before and after alemtuzumab treatments. Effector memory T cells and loss of T cell receptor (TCR) repertoire diversity are prevalent in T-LGLL. Shared TCRA and TCRB clonotypes are absent. Deregulation of cell survival and apoptosis gene programs, and marked downregulation of apoptosis genes in CD8+ clones, are prominent features of T-LGLL cells. Apoptosis genes are upregulated after alemtuzumab treatment, especially in responders than non-responders; baseline expression levels of apoptosis genes are predictive of hematologic response. Alemtuzumab does not attenuate TCR clonality, and TCR diversity is further skewed after treatment. Inferences made from analysis of single cell data inform understanding of the pathophysiologic mechanisms of clonal expansion and persistence in T-LGLL.

Long-term eltrombopag for bone marrow failure depletes iron.

Eltrombopag (EPAG) has been approved for the treatment of aplastic anemia and for immune thrombocytopenia, and a subset of patients require long-term therapy. Due to polyvalent cation chelation, prolonged therapy leads to previously underappreciated iron depletion. We conducted a retrospective review of patients treated at the NIH for aplastic anemia, myelodysplastic syndrome, and unilineage cytopenias, comparing those treated with EPAG to a historical cohort treated with immunosuppression without EPAG. We examined iron parameters, duration of therapy, response assessment, relapse rates, and common demographic parameters. We included 521 subjects treated with (n = 315) or without EPAG (n = 206) across 11 studies with multiyear follow-up (3.6 vs. 8.5 years, respectively). Duration of EPAG exposure correlated with ferritin reduction (p = 4 × 10-14 ) regardless of response, maximum dose, or degree of initial iron overload. Clearance followed first-order kinetics with faster clearance (half-life 15.3 months) compared with historical responders (47.5 months, p = 8 × 10-10 ). Risk of iron depletion was dependent upon baseline ferritin and duration of therapy. Baseline ferritin did not correlate with response of marrow failure to EPAG or to relapse risk, and timing of iron clearance did not correlate with disease response. In conclusion, EPAG efficiently chelates total body iron comparable to clinically available chelators. Prolonged use can deplete iron and ultimately lead to iron-deficiency anemia mimicking relapse, responsive to iron supplementation.

CAISC: A software to integrate copy number variations and single nucleotide mutations for genetic heterogeneity profiling and subclone detection by single-cell RNA sequencing.

BACKGROUND: Although both copy number variations (CNVs) and single nucleotide variations (SNVs) detected by single-cell RNA sequencing (scRNA-seq) are used to study intratumor heterogeneity and detect clonal groups, a software that integrates these two types of data in the same cells is unavailable. RESULTS: We developed Clonal Architecture with Integration of SNV and CNV (CAISC), an R package for scRNA-seq data analysis that clusters single cells into distinct subclones by integrating CNV and SNV genotype matrices using an entropy weighted approach. The performance of CAISC was tested on simulation data and four real datasets, which confirmed its high accuracy in sub-clonal identification and assignment, including subclones which cannot be identified using one type of data alone. Furthermore, integration of SNV and CNV allowed for accurate examination of expression changes between subclones, as demonstrated by the results from trisomy 8 clones of the myelodysplastic syndromes (MDS) dataset. CONCLUSIONS: CAISC is a powerful tool for integration of CNV and SNV data from scRNA-seq to identify clonal clusters with better accuracy than obtained from a single type of data. CAISC allows users to interactively examine clonal assignments.