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.

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.

Utility of plasma cell-free DNA for de novo detection and quantification of clonal hematopoiesis.

Although cell-free DNA (cfDNA) tests have emerged as a potential non-invasive alternative to bone marrow biopsies for monitoring clonal hematopoiesis in hematologic diseases, whether commercial cfDNA assays can be implemented for the detection and quantification of de novo clonal hematopoiesis in place of blood cells is uncertain. In this study, peripheral plasma cfDNA samples available from patients with aplastic anemia (n=25) or myelodysplastic syndromes (n=27) and a healthy cohort (n=107) were screened for somatic variants in genes related to hematologic malignancies using a Clinical Laboratory Improvement Amendments-certified panel. Results were further compared to DNA sequencing of matched blood cells. In reported results, 85% of healthy subjects, 36% of patients with aplastic anemia and 74% of patients with myelodysplastic syndromes were found to have somatic cfDNA variants, most frequently in DNMT3A, TET2, ASXL1 and SF3B1. However, concordance between cfDNA and blood cell findings was poor for the detection of clonal hematopoiesis when the allele frequency of the variants was <10%, which was mostly observed in the healthy and aplastic anemia cohorts but not in patients with myelodysplastic syndromes. After filtering data for potential artifacts due to low variant allele frequency and sequencing depth, the frequency of clonal hematopoiesis in cfDNA from healthy individuals and patients with aplastic anemia decreased to 52% and 20%, respectively. cfDNA and matched blood cells were not interchangeable for tracking changes in allele burdens as their agreement by Bland-Altman analysis was poor. A commercial cfDNA assay had good performance for de novo detection of clonal hematopoiesis in myelodysplastic syndromes, but showed no advantage over blood cells in diseases with low allele burdens or in healthy individuals.

Single-cell profiling of T lymphocytes in deficiency of adenosine deaminase 2.

Deficiency of adenosine deaminase 2 (DADA2) is a monogenic vasculitis syndrome caused by autosomal-recessive loss-of-function mutations in the ADA2 gene (previously known as CECR1). Vasculitis, vasculopathy, and inflammation are dominant clinical features of this disease; the spectrum of manifestations includes immunodeficiency and lymphoproliferation as well as hematologic manifestations. ADA2 is primarily secreted by stimulated monocytes and macrophages. Aberrant monocyte differentiation to macrophages and neutrophils are important in the pathogenesis of DADA2, but little is known about T lymphocytes in this disease. We performed combined single-cell RNA sequencing and single-cell TCR sequencing in order to profile T cell repertoires in 10 patients with DADA2. Although there were no significant alterations of T cell subsets, we observed activation of both CD8+ and CD4+ T cells. There was no clonal expansion of T cells: most TCRs were expressed at basal levels in patients and healthy donors. TCR usage was private to individual patients and not disease specific, indicating as unlikely a common pathogenic background or predisposition to a common pathogen. We recognized activation of IFN pathways as a signature of T cells and STAT1 as a hub gene in the gene network of T cell activation and cytotoxicity. Overall, T cells in DADA2 patients showed distinct cell-cell interactions with monocytes, as compared with healthy donors, and many of these ligand-receptor interactions likely drove up-regulation of STAT1 in both T cells and other immune cells in patients. Our analysis reveals previously undercharacterized cell characteristics in DADA2.

Analysis of deficiency of adenosine deaminase 2 pathogenesis based on single-cell RNA sequencing of monocytes.

Deficiency of adenosine deaminase 2 (DADA2) is a rare autosomal recessive disease caused by loss-of-function variants in the ADA2 gene. DADA2 typically presents in childhood and is characterized by vasculopathy, stroke, inflammation, immunodeficiency, as well as hematologic manifestations. ADA2 protein is predominantly present in stimulated monocytes, dendritic cells, and macrophages. To elucidate molecular mechanisms in DADA2, CD14+ monocytes from 14 patients and 6 healthy donors were analyzed using single-cell RNA sequencing (scRNA-seq). Monocytes were purified by positive selection based on CD14 expression. Subpopulations were imputed from their transcriptomes. Based on scRNA-seq, monocytes could be classified as classical, intermediate, and nonclassical. Further, we used gene pathway analytics to interpret patterns of up- and down-regulated gene transcription. In DADA2, the frequency of nonclassical monocytes was higher compared with that of healthy donors, and M1 macrophage markers were up-regulated in patients. By comparing gene expression of each monocyte subtype between patients and healthy donors, we identified upregulated immune response pathways, including IFNα/ß and IFNγ signaling, in all monocyte subtypes. Distinctively, the TNFR2 noncanonical NF-κB pathway was up-regulated only in nonclassical monocytes. Patients' plasma showed increased IFNγ and TNFα levels. Our results suggest that elevated IFNγ activates cell signaling, leading to differentiation into M1 macrophages from monocytes and release of TNFα. Immune responses and more general response to stimuli pathways were up-regulated in DADA2 monocytes, and protein synthesis pathways were down-regulated, perhaps as stress responses. Our identification of novel aberrant immune pathways has implications for therapeutic approaches in DADA2 (registered at clinicaltrials.gov NCT00071045).

Circulating exosomal microRNAs in acquired aplastic anemia and myelodysplastic syndromes.

Exosomal microRNAs modulate cancer cell metabolism and the immune response. Specific exosomal microRNAs have been reported to be reliable biomarkers of several solid and hematologic malignancies. We examined the possible diagnostic and prognostic values of exosomal microRNAs in two human bone marrow failure diseases: aplastic anemia and myelodysplastic syndromes. After screening 372 microRNAs in a discovery set (n=42) of plasma exosome samples, we constructed a customized PCR plate, including 42 microRNAs, for validation in a larger cohort (n=99). We identified 25 differentially expressed exosomal microRNAs uniquely or frequently present in aplastic anemia and/or myelodysplastic syndromes. These microRNAs could be related to intracellular functions, such as metabolism, cell survival, and proliferation. Clinical parameters and progression-free survival were correlated to microRNA expression levels in aplastic anemia and myelodysplastic syndrome patients before and after six months of immunosuppressive therapy. One microRNA, mir-126-5p, was negatively correlated with a response to therapy in aplastic anemia: patients with higher relative expression of miR-126-5p at diagnosis had the shortest progression-free survival compared to those with lower or normal levels. Our findings suggest utility of exosomal microRNAs in the differential diagnosis of bone marrow failure syndromes. (Registered at clinicaltrials.gov identifiers: 00260689, 00604201, 00378534, 01623167, 00001620, 00001397, 00217594).