RecQ helicase expression in patients with telomeropathies.

BACKGROUND: Telomeropathies are a group of inherited disorders caused by germline pathogenic variants in genes involved in telomere maintenance, resulting in excessive telomere attrition that affects several tissues, including hematopoiesis. RecQ and RTEL1 helicases contribute to telomere maintenance by unwinding telomeric structures such as G-quadruplexes (G4), preventing replication defects. Germline RTEL1 variants also are etiologic in telomeropathies. METHODS AND RESULTS: Here we investigated the expression of RecQ (RECQL1, BLM, WRN, RECQL4, and RECQL5) and RTEL1 helicase genes in peripheral blood mononuclear cells (PBMCs) from human telomeropathy patients. The mRNA expression levels of all RecQ helicases, but not RTEL1, were significantly downregulated in patients' primary cells. Reduced RecQ expression was not attributable to cell proliferative exhaustion, as RecQ helicases were not attenuated in T cells exhausted in vitro. An additional fifteen genes involved in DNA damage repair and RecQ functional partners also were downregulated in the telomeropathy cells. CONCLUSION: These findings indicate that the expression of RecQ helicases and functional partners involved in DNA repair is downregulated in PBMCs of telomeropathy patients.

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.

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.

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.

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.

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.

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.

HLA associations, somatic loss of HLA expression, and clinical outcomes in immune aplastic anemia.

Immune aplastic anemia (AA) features somatic loss of HLA class I allele expression on bone marrow cells, consistent with a mechanism of escape from T-cell-mediated destruction of hematopoietic stem and progenitor cells. The clinical significance of HLA abnormalities has not been well characterized. We examined the somatic loss of HLA class I alleles and correlated HLA loss and mutation-associated HLA genotypes with clinical presentation and outcomes after immunosuppressive therapy in 544 AA patients. HLA class I allele loss was detected in 92 (22%) of the 412 patients tested, in whom there were 393 somatic HLA gene mutations and 40 instances of loss of heterozygosity. Most frequently affected was HLA-B*14:02, followed by HLA-A*02:01, HLA-B*40:02, HLA-B*08:01, and HLA-B*07:02. HLA-B*14:02, HLA-B*40:02, and HLA-B*07:02 were also overrepresented in AA. High-risk clonal evolution was correlated with HLA loss, HLA-B*14:02 genotype, and older age, which yielded a valid prediction model. In 2 patients, we traced monosomy 7 clonal evolution from preexisting clones harboring somatic mutations in HLA-A*02:01 and HLA-B*40:02. Loss of HLA-B*40:02 correlated with higher blood counts. HLA-B*07:02 and HLA-B*40:01 genotypes and their loss correlated with late-onset of AA. Our results suggest the presence of specific immune mechanisms of molecular pathogenesis with clinical implications. HLA genotyping and screening for HLA loss may be of value in the management of immune AA. This study was registered at clinicaltrials.gov as NCT00001964, NCT00061360, NCT00195624, NCT00260689, NCT00944749, NCT01193283, and NCT01623167.

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

Comparative Transcriptomic Analysis of the Hematopoietic System between Human and Mouse by Single Cell RNA Sequencing.

(1) Background: mouse models are fundamental to the study of hematopoiesis, but comparisons between mouse and human in single cells have been limited in depth. (2) Methods: we constructed a single-cell resolution transcriptomic atlas of hematopoietic stem and progenitor cells (HSPCs) of human and mouse, from a total of 32,805 single cells. We used Monocle to examine the trajectories of hematopoietic differentiation, and SCENIC to analyze gene networks underlying hematopoiesis. (3) Results: After alignment with Seurat 2, the cells of mouse and human could be separated by same cell type categories. Cells were grouped into 17 subpopulations; cluster-specific genes were species-conserved and shared functional themes. The clustering dendrogram indicated that cell types were highly conserved between human and mouse. A visualization of the Monocle results provided an intuitive representation of HSPC differentiation to three dominant branches (Erythroid/megakaryocytic, Myeloid, and Lymphoid), derived directly from the hematopoietic stem cell and the long-term hematopoietic stem cells in both human and mouse. Gene regulation was similarly conserved, reflected by comparable transcriptional factors and regulatory sequence motifs in subpopulations of cells. (4) Conclusions: our analysis has confirmed evolutionary conservation in the hematopoietic systems of mouse and human, extending to cell types, gene expression and regulatory elements.

Comprehensive network modeling from single cell RNA sequencing of human and mouse reveals well conserved transcription regulation of hematopoiesis.

BACKGROUND: Presently, there is no comprehensive analysis of the transcription regulation network in hematopoiesis. Comparison of networks arising from gene co-expression across species can facilitate an understanding of the conservation of functional gene modules in hematopoiesis. RESULTS: We used single-cell RNA sequencing to profile bone marrow from human and mouse, and inferred transcription regulatory networks in each species in order to characterize transcriptional programs governing hematopoietic stem cell differentiation. We designed an algorithm for network reconstruction to conduct comparative transcriptomic analysis of hematopoietic gene co-expression and transcription regulation in human and mouse bone marrow cells. Co-expression network connectivity of hematopoiesis-related genes was found to be well conserved between mouse and human. The co-expression network showed "small-world" and "scale-free" architecture. The gene regulatory network formed a hierarchical structure, and hematopoiesis transcription factors localized to the hierarchy's middle level. CONCLUSIONS: Transcriptional regulatory networks are well conserved between human and mouse. The hierarchical organization of transcription factors may provide insights into hematopoietic cell lineage commitment, and to signal processing, cell survival and disease initiation.

Eltrombopag monotherapy can improve hematopoiesis in patients with low to intermediate risk-1 myelodysplastic syndrome.

Myelodysplastic syndromes (MDS) are a group of clonal myeloid disorders characterized by cytopenia and a propensity to develop acute myeloid leukemia (AML). The management of lower-risk (LR) MDS with persistent cytopenias remains suboptimal. Eltrombopag (EPAG), a thrombopoietin receptor agonist, can improve platelet counts in LR-MDS and tri-lineage hematopoiesis in aplastic anemia (AA). We conducted a phase 2 dose modification study to investigate the safety and efficacy of EPAG in LR-MDS. EPAG dose was escalated from 50 mg/day, to a maximum of 150 mg/day over a period of 16 weeks. The primary efficacy endpoint was hematologic response at 16-20 weeks. Eleven of 25 (44%) patients responded; five and six patients had uni- or bi-lineage hematologic responses, respectively. The predictors of response were presence of a PNH clone, marrow hypocellularity, thrombocytopenia with or without other cytopenia, and elevated plasma thrombopoietin levels at study entry. The safety profile was consistent with previous EPAG studies in AA; no patients discontinued drug due to adverse events. Three patients developed reversible grade-3 liver toxicity and one patient had increased reticulin fibrosis. Ten patients discontinued EPAG after achieving a robust response (median time 16 months); four of them reinitiated EPAG due to declining counts, and all attained a second robust response. Six patients had disease progression not associated with expansion of mutated clones and no patient progressed to AML on study. In conclusion, EPAG was well-tolerated and effective in restoring hematopoiesis in patients with low to intermediate-1 risk MDS. This study was registered at clinicaltrials.gov as #NCT00932156.

Comprehensive analysis of single-cell RNA sequencing data from healthy human marrow hematopoietic cells.

OBJECTIVE: Single cell methodology enables detection and quantification of transcriptional changes and unravelling dynamic aspects of the transcriptional heterogeneity not accessible using bulk sequencing approaches. We have applied single-cell RNA-sequencing (scRNA-seq) to fresh human bone marrow CD34+ cells and profiled 391 single hematopoietic stem/progenitor cells (HSPCs) from healthy donors to characterize lineage- and stage-specific transcription during hematopoiesis. RESULTS: Cells clustered into six distinct groups, which could be assigned to known HSPC subpopulations based on lineage specific genes. Reconstruction of differentiation trajectories in single cells revealed four committed lineages derived from HSCs, as well as dynamic expression changes underlying cell fate during early erythroid-megakaryocytic, lymphoid, and granulocyte-monocyte differentiation. A similar non-hierarchical pattern of hematopoiesis could be derived from analysis of published single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), consistent with a sequential relationship between chromatin dynamics and regulation of gene expression during lineage commitment (first, altered chromatin conformation, then mRNA transcription). Computationally, we have reconstructed molecular trajectories connecting HSCs directly to four hematopoietic lineages. Integration of long noncoding RNA (lncRNA) expression from the same cells demonstrated mRNA transcriptome, lncRNA, and the epigenome were highly homologous in their pattern of gene activation and suppression during hematopoietic cell differentiation.

Sequencing of RNA in single cells reveals a distinct transcriptome signature of hematopoiesis in GATA2 deficiency.

Constitutional GATA2 deficiency caused by heterozygous germline GATA2 mutations has a broad spectrum of clinical phenotypes, including systemic infections, lymphedema, cytopenias, and myeloid neoplasms. Genotype-phenotype correlation is not well understood mechanistically in GATA2 deficiency. We performed whole transcriptome sequencing of single hematopoietic stem and progenitor cells from 8 patients, who had pathogenic GATA2 mutations and myelodysplasia. Mapping patients' cells onto normal hematopoiesis, we observed deficiency in lymphoid/myeloid progenitors, also evident from highly constrained gene correlations. HSPCs of patients exhibited distinct patterns of gene expression and coexpression compared with counterparts from healthy donors. Distinct lineages showed differently altered transcriptional profiles. Stem cells in patients had dysregulated gene expression related to apoptosis, cell cycle, and quiescence; increased expression of erythroid/megakaryocytic priming genes; and decreased lymphoid priming genes. The prominent deficiency in lympho-myeloid lineages in GATA2 deficiency appeared at least partly due to the expression of aberrant gene programs in stem cells prior to lineage commitment. We computationally imputed cells with chromosomal abnormalities and determined their gene expression; DNA repair genes were downregulated in trisomy 8 cells, potentially rendering these cells vulnerable to second-hit somatic mutations and additional chromosomal abnormalities. Cells with complex cytogenetic abnormalities showed defects in genes related to multilineage differentiation and cell cycle. Single-cell RNA sequencing is powerful in resolving transcriptomes of cell subpopulations despite a paucity of cells in marrow failure. Our study discloses previously uncharacterized transcriptome signatures of stem cells and progenitors in GATA2 deficiency, providing a broad perspective of potential mechanisms by which germline mutations modulate early hematopoiesis in a human disease. This trial was registered at www.clinicaltrials.gov as NCT01905826, NCT01861106, and NCT00001620.

Deficit of circulating CD19+ CD24hi CD38hi regulatory B cells in severe aplastic anaemia.

Immune aplastic anaemia (AA) is caused by cytotoxic T lymphocytes (CTLs) that destroy haematopoietic stem and progenitor cells. Enhanced type 1 T helper (Th1) responses and reduced regulatory T cells (Tregs) are involved in the immune pathophysiology. CD24hi CD38hi regulatory B cells (Bregs) suppress CTLs and Th1 responses, and induce Tregs via interleukin 10 (IL-10). We investigated circulating B-cell subpopulations, including CD24hi CD38hi Bregs, as well as total B cells, CD4+ T cells, CD8+ T cells and natural killer cells in 104 untreated patients with severe and very severe AA, aged ≥18 years. All patients were treated with standard immunosuppressive therapy (IST) plus eltrombopag. CD24hi CD38hi Bregs were markedly reduced in patients with AA compared to healthy individuals, especially in very severe AA, but residual Bregs remained functional, capable of producing IL-10; total B-cell counts and the other B-cell subpopulations were similar to those of healthy individuals. CD24hi CD38hi Bregs did not correlate with responses to IST, and they recovered to levels present in healthy individuals after therapy. Mature naïve B-cell counts were unexpectedly associated with IST response. Markedly reduced CD24hi CD38hi Bregs, especially in very severe AA, with recovery after IST suggest Breg deficits may contribute to the pathophysiology of immune AA.

High throughput pSTAT signaling profiling by fluorescent cell barcoding and computational analysis.

Fluorescent cell barcoding (FCB) is a multiplexing technique for high-throughput flow cytometry (FCM). Although powerful in minimizing staining variability, it remains a subjective FCM technique because of inter-operator variability and differences in data analysis. FCB was implemented by combining two-dye barcoding (DyLight 350 plus Pacific Orange) with five-color surface marker antibody and intracellular staining for phosphoprotein signaling analysis. We proposed a robust method to measure intra- and inter-assay variability of FCB in T/B cells and monocytes by combining range and ratio of variability to standard statistical analyses. Data analysis was carried out by conventional and semi-automated workflows and built with R software. Results obtained from both analyses were compared to assess feasibility and reproducibility of FCB data analysis by machine-learning methods. Our results showed efficient FCB using DyLight 350 and Pacific Orange at concentrations of 0, 15 or 30, and 250 µg/mL, and a high reproducibility of FCB in combination with surface marker and intracellular antibodies. Inter-operator variability was minimized by adding an internal control bridged across matrices used as rejection criterion if significant differences were present between runs. Computational workflows showed comparable results to conventional gating strategies. FCB can be used to study phosphoprotein signaling in T/B cells and monocytes with high reproducibility across operators, and the addition of bridge internal controls can further minimize inter-operator variability. This FCB protocol, which has high throughput analysis and low intra- and inter-assay variability, can be a powerful tool for clinical trial studies. Moreover, FCB data can be reliably analyzed using computational software.