Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 10 February 2026

Sec. Microbial Immunology

Volume 17 - 2026 | https://doi.org/10.3389/fimmu.2026.1709981

This article is part of the Research TopicDeciphering host-pathogen interactions in tuberculosis: implications for diagnostics and therapeutics Volume IIView all 4 articles

Peripheral blood B-cell compartment dysregulation in multidrug-resistant tuberculosis is associated with reduced circulating marginal zone-like B cells

Mei Liu,,&#x;Mei Liu1,2,3†Junjie Wen,&#x;Junjie Wen4,5†Qin Gao&#x;Qin Gao1†Tao XuTao Xu6Yuanbo Lan,,*Yuanbo Lan1,2,3*Lu Meng,*Lu Meng4,6*
  • 1Department of Tuberculosis, Affiliated Hospital of Zunyi Medical University, Zunyi, Guizhou, China
  • 2Guizhou Provincial Key Laboratory of Pathogenesis and Prevention for Common Chronic Diseases, Affiliated Hospital of Zunyi Medical University, Zunyi, Guizhou, China
  • 3Guizhou Provincial 2011 Collaborative Innovation Center for Comprehensive Tuberculosis Prevention and Treatment, Affiliated Hospital of Zunyi Medical University, Zunyi, Guizhou, China
  • 4Institute of Infection and Health, Fudan University, Shanghai, China
  • 5Institute of Science, Fudan University, Shanghai, China
  • 6Department of Infectious Diseases, Huashan Hospital, Shanghai, China

Objective: Multidrug-resistant tuberculosis (MDR-TB) remains a major global health challenge. While T cell-mediated immunity in tuberculosis is well characterized, alterations in circulating B-cell subsets during chronic MDR-TB are less well defined.

Methods: Peripheral blood mononuclear cells (PBMCs) from healthy controls [interferon gamma release assay negative (IGRA)], individuals with latent tuberculosis infection (LTBI; IGRA+), and patients with active tuberculosis (ATB) were analyzed using multiparameter flow cytometry panels. Major lymphoid and myeloid populations and detailed B-cell subsets were quantified.

Results: Frequencies of major T-cell and natural killer (NK)-cell populations were broadly similar across groups. In contrast, patients with ATB showed a reduction in total CD19+ B cells. Within the B-cell compartment, ATB was characterized by an increased proportion of naïve B cells and a pronounced reduction in antibody-secreting cells (ASCs). Circulating marginal zone-like B cells (MZ B, IgD+IgM+CD27+) were also reduced in ATB compared with non-ATB groups. Receiver operating characteristic (ROC) analysis suggested that reduced MZ B-cell frequency may help discriminate individuals with ATB from those without ATB; however, this observation should be interpreted as exploratory given the cohort size and composition.

Conclusion: MDR-TB is associated with broad perturbations of the peripheral B-cell compartment, including reduced ASCs and decreased circulating MZ B cells. These findings highlight B-cell dysregulation as a feature of active disease and identify MZ B cells as a subset of interest for further investigation rather than as a stand-alone diagnostic marker.

Introduction

Tuberculosis (TB), caused by the airborne pathogen Mycobacterium tuberculosis (Mtb), remains a major global health challenge. The emergence of multidrug-resistant TB (MDR-TB), defined as resistance to at least isoniazid and rifampicin, poses a catastrophic threat to global TB control (1). According to World Health Organization reports, TB claimed approximately 1.3 million lives in 2023 alone, surpassing mortality rates of all other single infectious pathogens. This persistent burden underscores the critical need for improved diagnostic and therapeutic strategies.

Following alveolar invasion, Mtb induces a spectrum of infection states modulated by host immune competence (24). Clinically, TB status is categorized into three distinct states: active TB (ATB), characterized by symptomatic, transmissible infection; latent TB infection [LTBI; interferon gamma release assay-positive (IGRA+)], representing asymptomatic microbial persistence; and healthy controls (IGRA). Emerging evidence indicates that circulating immune cell subsets undergo differential modulation during Mtb disease progression. While cell-mediated immunity, particularly CD4+ T-cell responses, is central to mycobacterial control (5, 6), the complex interplay of multiple lymphocyte populations is increasingly recognized as critical. Th1 cells, pivotal through interferon-gamma (IFN-γ) production, act alongside functional counterparts like Th17 cells and Th1–Th17 hybrid populations, which demonstrate complementary roles in host defense (7, 8). Concurrently, innate and innate-like immune cells are significant for differentiating infection stages. Neutrophils dominate airway infiltrates in patients with ATB (9), playing dual roles in bacterial clearance and immunopathology (10). Natural killer (NK)-cell dynamics also reflect disease progression, with a progressive depletion of the CD3CD7+GZMB+ NK subset observed from healthy controls through ATB, followed by post-treatment recovery (11, 12).

Given this circulating cellular heterogeneity, comprehensive immunophenotyping through surface marker analysis offers valuable insights for biomarker discovery. Our systematic investigation evaluates multiple immune populations implicated in anti-mycobacterial responses, including CD4+/CD8+ T cells, γδ T cells, NKT cells, NK cells, and mucosal-associated invariant T (MAIT) cells (identified by TCR Vα7.2 expression) (13). Critically, we extend this characterization to B-cell subsets, whose roles in TB pathogenesis have been historically underappreciated despite emerging evidence of their immunological significance.

The conventional paradigm viewed B cells in intracellular pathogen defense primarily through the lens of humoral immunity and immunoglobulin production (14). However, contemporary research reveals multifaceted roles for B cells, including antigen presentation, cytokine modulation, and T-cell regulation (1517). B cells can internalize mycobacteria via receptor-mediated uptake and macropinocytosis, facilitating MHC class II antigen presentation. Mirroring T-cell diversity, B-cell populations demonstrate functional specialization through distinct subpopulations: naïve B cells (Bn), transitional B cells (Btr), antibody-secreting cells (ASCs), and memory B cells (MBCs) (18, 19).

Current literature, however, presents conflicting reports on peripheral B-cell dynamics during TB progression. Comparative analyses of ATB versus healthy controls describe unchanged (20), elevated (21), or diminished (22, 23) total B-cell frequencies. Similarly, contradictory patterns emerge in LTBI cohorts showing decreased B-cell counts (22), while post-treatment patients exhibit B-cell expansion (20). This ambiguity highlights the need for a more refined analysis of B-cell subset distribution and differentiation across the TB disease spectrum.

In this study, we performed multiparameter flow cytometric profiling of peripheral blood immune cells from patients with MDR-TB, individuals with LTBI, and healthy controls. We aimed to characterize broad immune alterations with particular emphasis on B-cell subset distribution, while carefully acknowledging the limitations inherent to peripheral blood-based and cross-sectional analyses.

Materials and methods

Study population and ethical considerations

This study was approved by the ethics committee of the affiliated hospital of Zunyi Medical University, China (No. KLL-2022-579). Individuals with ATB, LTBI (IGRA+), and healthy controls (IGRA) were recruited from the Affiliated Hospital of Zunyi Medical University. IGRA+ and IGRA participants were recruited from hospital staff undergoing routine health screening and community volunteers without TB-related symptoms. Asymptomatic individuals with positive IGRA results were confirmed by TB-IGRA (Beijing Wantai Biological Pharmacy Enterprise Co., Ltd.). Patients with ATB were microbiologically confirmed and, in this cohort, all ATB cases had drug-resistant TB and were receiving second-line anti-TB treatment at the time of sampling. This characteristic of the cohort is acknowledged as a limitation and may affect the generalizability of the findings to drug-sensitive TB.

The sex distribution of the enrolled cohorts, particularly the male predominance in the ATB group, reflects the established epidemiological profile of TB, where global incidence is consistently higher in male patients. This disparity is representative of the local clinical population from which participants were recruited.

PBMC isolation and processing

Peripheral blood samples (8–10 mL per participant) were collected in sodium heparin vacutainers (BD Biosciences) and processed within 2 h to ensure cell viability. Peripheral blood mononuclear cells (PBMCs) were isolated using Ficoll-Hypaque (GE Healthcare) density gradient centrifugation at 400 × g for 30 min at room temperature. The PBMC layer was carefully aspirated, washed twice with phosphate-buffered saline (PBS, Gibco), and resuspended in complete RPMI 1640 medium (Corning) supplemented with 10% heat-inactivated fetal bovine serum (FBS, Gibco).

Multicolor flow cytometry analysis

Surface staining: PBMCs were thawed and washed. The staining and gating strategy was adapted from previously published phenotypic analyses in TB (18). In brief, a total of 1 × 106 cells were stained with a pre-optimized antibody panel (details and clones provided in Supplementary Table S1) in Brilliant Stain Buffer (BD Biosciences) for 30 min at 4 °C in the dark.

Data acquisition and analysis: Stained cells were acquired on a Cytek Aurora spectral flow cytometer (Cytek Biosciences) configured with five lasers (355, 405, 488, 561, and 640 nm). Compensation was performed using single-stained UltraComp eBeads (Thermo Fisher Scientific). Gating strategies for B-cell subpopulations (based on CD27, IgD, and CD10 expression) were established using Fluorescence Minus One (FMO) controls, as depicted in Supplementary Figure S2. Data were analyzed using FlowJo v10.7.1 (BD Biosciences). For visualization, high-dimensional data were subjected to dimensionality reduction using the Uniform Manifold Approximation and Projection (tSNE) algorithm in FCS Express™ 7 Software (Denovo Software). Cell populations were quantified as percentages of parent populations (%).

Interferon gamma release assays

The IGRA test was employed to evaluate the TB infection status of all participants. IGRA was performed according to the manufacturer’s instructions (Wantai, Beijing, China). Briefly, blood samples were obtained and divided into three tubes: (1) a negative control tube (Nil tube), (2) an Mtb-specific antigen stimulation tube (Test tube), and (3) a positive control tube containing PHA (PHA tube). After incubation at 37°C for approximately 22 h, the supernatants were collected, and IFN-γ was quantitated by immunoassays. IGRA results were judged according to the manufacturer’s instruction.

Statistical analysis

Continuous variables are presented as mean ± standard deviation (SD) for normally distributed data or median [interquartile range (IQR)] for non-parametric data. Normality was assessed via Shapiro–Wilk test. The following group comparisons were employed:

Parametric: Student’s t-test (two groups) or one-way analysis of variance (ANOVA) with Tukey’s post-hoc test (≥3 groups).

Non-parametric: Mann–Whitney U test (two groups) or Kruskal–Wallis test with Dunn’s correction (≥3 groups).

Categorical variables were analyzed using Fisher’s exact test or Pearson’s χ2 test with Yates’ continuity correction. Correlation analyses utilized Spearman’s rank correlation coefficient. All tests were two-tailed, with p < 0.05 considered statistically significant. Analyses were performed using GraphPad Prism v8.4.2 (GraphPad Software).

Results

Participant characteristics

Demographic and clinical characteristics of the participants are summarized in Table 1. The mean ages were comparable across the ATB, IGRA+, and IGRA groups, supporting the validity of group-wise immunological comparisons.

Table 1
www.frontiersin.org

Table 1. Clinical characteristics of participants.

Comprehensive immune cell profiling of peripheral blood subsets via flow cytometry

A total of 24 participants were enrolled in this study and classified into three distinct groups: 8 patients with ATB, 8 individuals with LTBI (IGRA+), and 8 healthy controls (IGRA). PBMCs were isolated from each participant and cryopreserved for subsequent high-dimensional immunophenotyping by flow cytometry (Figure 1A). All samples met quality thresholds for cell count and viability to ensure robustness in downstream analyses.

Figure 1
Flowchart illustrating immune cell subsets and identification process. Panel A: Groups IGRA-, IGRA+, and ATB are linked to Panel 1 (DC, NK, Monocyte, T, B cells) and Panel 2 (Naïve B, Memory B, ASC). Panel B: Differentiation of monocytes and lymphocytes, highlighting non-classical, intermediate, and classical monocytes, along with CD3+, CD19+, CD56+ lymphocytes. Details include markers like CD14, CD16, HLA DR, and CD45, with subsets such as Treg, NKT, γδT, MAIT, and NK cells. Cell markers and differentiation pathways are indicated by arrows.

Figure 1. Study workflow and immune phenotyping of human PBMCs. (A) Schematic overview of the experimental procedure. PBMCs were isolated from all participants across the three cohorts (active tuberculosis, ATB; latent TB infection, IGRA+; and healthy controls, IGRA–). After cryopreservation, PBMCs were thawed and subjected to high-dimensional immunophenotyping using a predefined 21-marker panel by flow cytometry (including 16 markers in Panel 1 and 8 markers in Panel 2). The resulting high-parameter data were analyzed to identify immune cell populations and states associated with TB disease status. (B) Representative flow cytometry plots illustrating the sequential gating strategy used to identify major immune cell populations and subsets. The gating hierarchy began with the identification of single, live lymphocytes based on forward-scatter (FSC) and side-scatter (SSC) properties, followed by successive steps to delineate innate and adaptive immune subsets. Eighteen immune cell populations identified in Panel 1 include CD19+CD3- B cells, CD3+CD19- T cells, plasmacytoid dendritic cells (pDCs), conventional/myeloid dendritic cells (cDCs/mDCs), CD4+ T cells, CD8+ T cells, regulatory T cells (Treg), natural killer T (NKT) cells, γδ T cells, mucosal-associated invariant T (MAIT) cells, innate lymphoid cells (ILC), NK cells (including CD56bright and CD56dim subsets), and monocyte subsets (classical, intermediate, and non-classical monocytes). B-cell subsets shown in Panel 2 include memory B cells (MBCs), antibody-secreting cells (plasma/ASCs), naïve B cells (Bn), circulating marginal zone–like B cells (MZ B), and transitional B cells (Btr). This systematic approach ensured consistent and reproducible subset discrimination across all samples.

To delineate immune dynamics in ATB, we designed a multiparameter flow cytometry panel incorporating the following markers: CD45, CD3, CD19, CD56, CD11c, CD14, CD16, HLA-DR, CD123, CD303, CD127, CD25, CD4, CD8, TCRγδ, TCRVα7.2, CD10, IgD, CD27, CD38, and IgM. This panel enabled the identification of major immune lineages and subsets (Figure 1B).

Leukocyte populations were gated via CD45 expression and stratified into lymphoid (CD3+ T cells, CD19+ B cells, CD3+CD56+ NKT-like cells) and myeloid lineages. T-cell subsets were defined as CD4+/CD8+ subsets. Specialized T populations: MAIT cells (TCR Vα7.2+) and regulatory T cells (Tregs: CD25+CD127lowCD4+).

Myeloid compartment analysis from CD3CD19 cells included the following: Monocyte subsets: classical (CD14+CD16), intermediate (CD14+CD16+), and nonclassical (CD14dimCD16+). NK cell stratification: cytotoxic CD56dimCD16high vs. cytokine-secreting CD56brightCD16low subsets, with activation status assessed via HLA-DR expression. Dendritic cells (DCs): Identified as HLA-DR+CD11c+CD16.

To ensure cross-sample comparability, cryopreserved PBMCs from patients with ATB, IGRA+, and IGRA cohorts were thawed and analyzed synchronously within identical experimental runs, minimizing batch variability.

General immune cell profiling via flow cytometry

To characterize the immune landscape in ATB, we used a 16-color flow cytometry panel (Panel 1) to evaluate the frequencies of major immune cell types, including T cells, NK cells, monocytes, and B cells (Figure 2A). Gating was performed hierarchically: B cells were identified as CD19+CD3- lymphocytes. Among CD19-CD3+ T-lineage cells, NKT cells were defined as CD3+CD56+, while the remaining CD3+CD56- cells were further separated into γδ T cells (TCRγδ+) and TCRγδ- conventional T cells. From the latter, MAIT cells were gated as TCRVα7.2+, and the remaining TCRVα7.2- conventional T cells were subdivided into CD4+ T cells (CD4+CD8-) and CD8+ T cells (CD4-CD8+). Regulatory T cells (Tregs) were subsequently identified as CD25+CD127- within the CD4+ T-cell population.

Figure 2
Flow cytometry plots and bar graphs are depicted. Panel A shows gating strategies for identifying various immune cell subsets from peripheral blood mononuclear cells (PBMC), including CD3+, CD3-CD19-, and CD3-CD19- populations with detailed markers like CD56, CD16, and HLA-DR. Panel B displays bar graphs showing percentages of immune cell types, including pDC, cDC/mDC, and CD19+ cells, across IGRA- (Interferon Gamma Release Assay negative), IGRA+, and ATB (Active Tuberculosis) groups with statistical markers indicating significance.

Figure 2. Alterations in circulating immune cell populations in active tuberculosis. (A) Gating strategy for the identification of major immune cell types, including T cells (CD4+, CD8+, and Tregs), NK cells, monocytes, and B cells, using the 16-color flow cytometry panel (Panel 1). (B) Comparative frequencies of dendritic cell (DC) subsets [plasmacytoid DCs (pDCs) and conventional DCs (cDCs)] and total CD19+ B cells among the three participant groups: active tuberculosis (ATB), latent tuberculosis infection (IGRA+), and healthy controls (IGRA). Data are presented as mean ± SD. Statistical significance was determined by one-way ANOVA; *p < 0.05.

Within the CD3-CD19- compartment, innate lymphoid cells (ILCs) were defined as CD127+CD123-, and NK cells were gated as CD127-CD123-CD56+ and further stratified into CD56bright (CD16-CD56+) and CD56dim (CD16+CD56+) subsets. Plasmacytoid dendritic cells (pDCs) were identified as CD123+CD303+HLA-DR+. From the CD11c+ fraction of the CD3-CD19- lineage, conventional dendritic cells (cDCs) were defined as CD11c+CD14-, and monocyte subsets were classified as follows: classical monocytes (CD14highCD16low), intermediate monocytes (CD14highCD16highCD11chighHLADRhigh), and non-classical monocytes (CD14lowCD16highCD11chighHLADR+).

This gating strategy enabled consistent and reproducible identification of 18 distinct immune cell subsets, facilitating the analysis of their compositional changes across clinical groups. Consistent with previous reports (24, 25), we observed no significant differences in the frequencies of CD4+ T, CD8+ T, or regulatory T cells (Tregs) between patients with ATB and healthy controls (IGRA). Similarly, NK cell counts remained comparable across groups, suggesting that both circulating cellular adaptive and innate effector populations are largely maintained during active disease (Supplementary Figure S1). However, Supplementary Figure 1 reveals significant alterations in specific innate-like lymphocyte populations. Notably, both patients with ATB and IGRA+ individuals exhibited a reduction in circulating γδ T cells compared to IGRA controls, with the difference reaching statistical significance for patients with ATB. Additionally, MAIT cell frequencies were significantly lower in patients with ATB than in IGRA individuals.

In contrast, DCs—including both plasmacytoid (pDC) and conventional/myeloid (cDC/mDC) subsets—were significantly reduced in patients with ATB compared to both IGRA+ and IGRA– groups (Figure 2B). Most notably, we observed a striking reduction in total peripheral B cells (CD19+) in patients with ATB relative to both IGRA+ individuals and healthy controls (p < 0.05; Figure 2B). This finding is consistent with growing evidence of B-cell dysregulation in TB. Depletion of circulating B cells may reflect impaired humoral immune function, affecting key mechanisms such as antibody production, antigen presentation to T cells, and immunomodulatory cytokine signaling. The observed B-cell deficit suggests a broad disruption of adaptive immunity that could facilitate Mtb persistence by undermining coordinated immune responses.

B-cell maturation arrest and functional impairment in ATB

B-cell development begins in the bone marrow, where immature B cells are generated and then enter the peripheral blood. These immature B cells proceed through several transitional stages, named transitional B cells (Btr), before migrating to secondary lymphoid organs, such as the spleen. In the spleen, B cells undergo rigorous antigen selection, which influences their differentiation into mature B cells, including MBCs and plasma cells, which are both involved in the adaptive immune response.

Commonly, in the peripheral blood, B-cell ontogeny involves sequential maturation from bone marrow-derived transitional B cells (Btr: CD10+) to splenic naïve B cells (Bn: IgD+CD27), followed by antigen-driven differentiation into MBCs (IgDCD27+) and ASCs (CD19+CD38high) (18). Flow cytometry profiling (Panel 2: CD45/CD3/CD19/CD10/IgM/IgD/CD27/CD38) revealed profound perturbations in ATB: expanded naïve B cells (Bn: +40% vs. controls, p < 0.01), depleted MZ B cells (−59%, p < 0.05), and ASCs (−65%, p < 0.01) (Figures 3A-C). The significant expansion of Bn concurrent with a severe reduction in ASCs and MZ B cells indicates a profound B-cell maturation arrest, likely resulting from Mtb-induced suppression of germinal center formation or impaired antigen-presenting cell function. The paucity of ASCs further indicates compromised antibody production capacity, which may undermine opsonization and complement activation against extracellular bacilli.

Figure 3
A. Three tSNE plots compare B cell populations in IGRA-negative, IGRA-positive, and ATB samples. The percentages and density distributions differ among the groups. B. Flow cytometry gating strategy for identifying B cell subtypes, including mature, transitional, and memory B cells. Gates are applied sequentially from PBMC to specific B cell markers. C. Bar graphs display percentages of CD19-positive cells in various B cell subsets across IGRA-negative, IGRA-positive, and ATB groups, with significant differences noted for certain subsets.

Figure 3. Distribution of major circulating B-cell subsets in peripheral blood from IGRA (n = 8), IGRA+ (n = 8) and patients with ATB (n = 8). (A) t-SNE visualization of PBMCs with cells colored according to CD19 expression levels. CD19+ B cell populations are highlighted with green ellipses in the IGRA, IGRA+, and ATB groups. (B) After lymphocyte gate, the CD19+ (B cell) population gates and B-cell subsets were identified by distinct surface markers. (C) Frequencies of total B cells (CD19+) as % of live lymphocytes, transitional B cells (CD19+CD10+), and mature B cells (CD19+CD10) as % of total B cells, naïve B cells (CD19+CD10IgD+CD27), MZ B cells (CD19+CD10IgD+CD27+), and ASCs (CD19+IgDCD27+CD38+) as % of mature B cells. Statistical analyses were done by one-way ANOVA test. Significant differences are shown on top of the respective columns (p < 0.05; *p < 0.01). ns, not significant.

Reduction in circulating marginal zone B cells serves as a diagnostic indicator for active tuberculosis

Building upon our previous findings of a generalized B-cell maturation arrest and functional impairment in ATB, we sought to identify a specific and clinically measurable B-cell subset alteration. To this end, we developed a targeted five-color flow cytometry panel (CD3, CD19, CD27, IgD, and IgM) for practical immunophenotyping. After gating on lymphocytes and identifying CD19+ B cells (using CD3 for exclusion), we defined MZ B cells as the IgD+CD27+IgM+ population (Figure 4A). Strikingly, the proportion of this specific MZ B subset within total B cells was significantly contracted in patients with ATB compared to the non-ATB control group (Figure 4B). This control cohort included individuals with LTBI, closely exposed healthcare workers, and patients with other pulmonary diseases, strengthening the specificity of our observation. This finding suggests that the loss of circulating MZ B cells is a specific feature of active disease, not merely a consequence of recent exposure or general pulmonary inflammation.

Figure 4
A composite image showing multiple scientific charts and plots. A: Flow cytometry plot identifying B-cells using markers CD3, CD19, CD27, and IgD. B: Bar graph comparing percentage of B-cells between non-ATB and ATB groups with a significant difference indicated. C: Normal QQ plot displaying predicted versus actual values for non-ATB and ATB, with points along a dotted line. D: ROC curve illustrating sensitivity versus one minus specificity. E: Line chart comparing sensitivity and specificity percentages across different thresholds.

Figure 4. Decreased frequency of circulating marginal zone-like B cells in active tuberculosis and its diagnostic performance. (A) Gating strategy for the identification of human MZ B cells. Peripheral blood mononuclear cells were stained with a panel of antibodies (CD3, CD19, CD27, IgD, and IgM). After gating on lymphocytes and excluding CD3+ T cells, B cells were identified as CD19+ lymphocytes. MZ B cells were defined within the CD19+ B-cell population as IgD+CD27+IgM+. (B) Comparison of circulating MZ B-cell frequencies between ATB and non-ATB groups. MZ B-cell frequency is expressed as a percentage of total CD19+ B cells. Statistical differences between groups were analyzed using an unpaired t-test. (C) Normality of data distribution for the frequency of MZ B cells (as a percentage of total B cells) in the ATB and non-ATB groups was confirmed using Anderson–Darling, D’Agostino and Pearson, Shapiro–Wilk, and Kolmogorov–Smirnov tests (all p > 0.05). (D) Receiver operating characteristic (ROC) curve analysis evaluating the diagnostic performance of the MZ B-cell frequency for discriminating patients with ATB from non-ATB controls. The analysis resulted in an area under the curve (AUC) of 0.8025 (95% CI: 0.5915 to 1.000, p = 0.0305). (E) Sensitivity and specificity as a function of the classification threshold for MZ B-cell frequency.

The data for both groups passed all tests for normality (Anderson-Darling, D’Agostino and Pearson, Shapiro–Wilk, and Kolmogorov–Smirnov; all p > 0.05) (Figure 4C), validating the use of parametric statistical tests. Most importantly, receiver operating characteristic (ROC) curve analysis demonstrated that this simple MZ B immunophenotype held significant diagnostic power. The analysis yielded an area under the curve (AUC) of 0.8025 (95% CI: 0.5915 to 1.000, p = 0.0305) (Figure 4D), indicating good utility in discriminating ATB from non-ATB conditions. In addition, analysis of sensitivity and specificity across varying classification thresholds further illustrated the performance characteristics of MZ B-cell frequency as a diagnostic marker (Figure 4E). Therefore, the quantifiable reduction in MZ B cells provides a tangible and translatable biomarker. The ability to capture this aspect of the dysfunctional B-cell response using a minimal, clinically feasible antibody panel highlights its utility for immunodiagnostic application in ATB.

Discussion

This study offers valuable insights into the immune dysregulation observed in MDR-TB, with a particular emphasis on B-cell alterations. Our findings reveal a significant reduction in total B cells, alongside a marked depletion of ASCs in the peripheral blood of patients with MDR-TB, which may contribute to the persistence and chronicity of Mtb infection (26, 27).

Additionally, the depletion of MZ B cells, which are critical for rapid immune responses to encapsulated pathogens, further implicates B-cell dysfunction in the pathogenesis of TB. Recent studies have reported the important role of MZ B cells in the lungs, the primary organ affected by TB, in orchestrating effective immune responses (28, 29). In our study, we also observed a close relationship between MZ B cells and Mtb infection in peripheral blood, suggesting that these cells may continue to play a role throughout the course of infection. Furthermore, the use of peripheral blood as a readily accessible source for studying MZ B cells has important diagnostic implications, given the challenges in obtaining lung tissue from patients with ATB.

These findings challenge the traditional T cell-centric view of TB immunopathology. While our data confirm preserved T/NK-cell frequencies in ATB (consistent with prior reports), the concurrent collapse of humoral effector mechanisms—naïve B-cell accumulation, ASC deficiency, and MZ B-cell depletion—implies that Mtb subverts multiple immune axes. This aligns with growing recognition of B cells as active participants in TB defense through antibody-independent roles like cytokine modulation and lymphoid follicle organization (1517).

Future studies should address whether MZ B-cell depletion directly enables bacterial dissemination or merely marks advanced disease. Longitudinal tracking of this subset during treatment could clarify its biomarker utility, while spatial transcriptomics of granulomas may resolve the “missing MZ B cell” paradox.

Our study acknowledges several inherent limitations. First, the modest cohort size and the predominance of male patients in the ATB group curtail statistical power and may affect the generalizability of the results. Second, the focus exclusively on patients with MDR-TB means the findings may not be directly applicable to those with drug-sensitive TB, reflecting a more severe immunopathological state. Third, from a methodological perspective, the use of peripheral blood as a proxy for systemic immunity, though pragmatic, carries constraints. The absence of paired tissue samples from sites such as the lung, spleen, or bone marrow prevents definitive characterization of MZ B-cell trafficking and leaves unresolved the functional competence (e.g., IgM secretion) of the circulating MZ B cells we identified.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the ethics committee of the affiliated hospital of Zunyi Medical University, China. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

ML: Writing – original draft, Data curation, Formal analysis, Funding acquisition, Methodology, Validation. JW: Data curation, Writing – original draft, Formal analysis, Methodology, Validation, Software, Visualization. QG: Data curation, Methodology, Investigation, Writing – original draft. TX: Methodology, Resources, Writing – review & editing. YL: Resources, Funding acquisition, Project administration, Writing – review & editing. LM: Project administration, Conceptualization, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Guizhou Provincial Science and Technology Program (Qiankehe Foundation-ZK(2023)581), the Guizhou Provincial Scientific and Technological Program (Qiankehe Foundation-ZK(2023)Key 060), and the Zunyi Science and Technology Program (HZ(2022)254).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer BY declared a shared parent affiliation with the author(s) JW, TX, and LM to the handling editor at the time of review.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2026.1709981/full#supplementary-material

Supplementary Figure 1 | Comprehensive analysis of innate and adaptive immune cell subsets across study groups. Frequencies of various immune cell populations within PBMCs are shown for the active tuberculosis (ATB), latent tuberculosis (IGRA+), and healthy control (IGRA-) groups. Analyzed subsets include: monocytes (classical [CD14+CD16-], intermediate [CD14+CD16+], non-classical [CD14+CD16++]); innate lymphoid cells (ILCs); NK cells (CD56bright and CD56dim subsets); NKT cells; γδ T cells; and T cell subsets (Tregs, CD4+ T, CD8+ T, and MAIT cells). Data are presented as mean ± SD. Statistical significance was determined by one-way ANOVA with appropriate post-hoc tests; *p < 0.05, ns = not significant.

Supplementary Figure 2 | Flow cytometric analysis of B-cell subpopulations based on CD27, IgD, and CD10 expression across study groups. In each plot, the Fluorescence Minus One (FMO) control (left column) is shown with a defined FMO Gate (dashed rectangle or polygon). This gate sets the threshold for positive vs. negative populations for the respective y-axis marker.

Supplementary Table 1 | Antibody panel for spectral flow cytometry.

References

1. World Health Organization. Global Tuberculosis Report 2023. Geneva: World Health Organization (2023).

Google Scholar

2. Flynn JL and Chan J. Immunology of tuberculosis. Annu Rev Immunol. (2001) 19:93–129. doi: 10.1146/annurev.immunol.19.1.93

PubMed Abstract | Crossref Full Text | Google Scholar

3. Cooper AM. Cell-mediated immune responses in tuberculosis. Annu Rev Immunol. (2009) 27:393–422. doi: 10.1146/annurev.immunol.021908.132703

PubMed Abstract | Crossref Full Text | Google Scholar

4. O’Garra A, Redford PS, McNab FW, Bloom CI, Wilkinson RJ, and Berry MP. The immune response in tuberculosis. Annu Rev Immunol. (2013) 31:475–527. doi: 10.1146/annurev-immunol-032712-095939

PubMed Abstract | Crossref Full Text | Google Scholar

5. Caruso AM, Serbina N, Klein E, Triebold K, Bloom BR, and Flynn JL. Mice deficient in CD4 T cells have only transiently diminished levels of IFN-γ, yet succumb to tuberculosis. J Immunol. (1999) 162:5407–16. doi: 10.4049/jimmunol.162.9.5407

PubMed Abstract | Crossref Full Text | Google Scholar

6. Gallegos AM, Pamer EG, and Glickman MS. Delayed protection by ESAT-6-specific effector CD4+ T cells after airborne M. tuberculosis infection. J Exp Med. (2008) 205:2359–68. doi: 10.1084/jem.20080353

PubMed Abstract | Crossref Full Text | Google Scholar

7. Tong ZH and Shi HZ. Subpopulations of helper T lymphocytes in tuberculous pleurisy. Tuberculosis (Edinb). (2013) 93:279–84. doi: 10.1016/j.tube.2013.02.014

PubMed Abstract | Crossref Full Text | Google Scholar

8. Shanmugasundaram U, Bucsan AN, Ganatra SR, Ibegbu C, Quezada M, Blair RV, et al. Pulmonary Mycobacterium tuberculosis control associates with CXCR3- and CCR6-expressing antigen-specific Th1 and Th17 cell recruitment. JCI Insight. (2020) 5:e137858. doi: 10.1172/jci.insight.137858

PubMed Abstract | Crossref Full Text | Google Scholar

9. Eum SY, Kong JH, Hong MS, Lee YJ, Kim JH, Hwang SH, et al. Neutrophils are the predominant infected phagocytic cells in the airways of patients with active pulmonary TB. Chest. (2010) 137:122–8. doi: 10.1378/chest.09-0903

PubMed Abstract | Crossref Full Text | Google Scholar

10. Dallenga T, Linnemann L, Paudyal B, Repnik U, Griffiths G, and Schaible UE. Targeting neutrophils for host-directed therapy in tuberculosis. Int J Med Microbiol. (2018) 308:142–7. doi: 10.1016/j.ijmm.2017.10.001

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cai Y, Dai Y, Wang Y, Yang Q, Guo J, Wei C, et al. Single-cell transcriptomics of blood reveals a natural killer cell subset depletion in tuberculosis. EBioMedicine. (2020) 53:102686. doi: 10.1016/j.ebiom.2020.102686

PubMed Abstract | Crossref Full Text | Google Scholar

12. Roy Chowdhury R, Vallania F, Yang Q, Lopez Angel CJ, Darboe F, Penn-Nicholson A, et al. A multi-cohort study of the immune factors associated with M. tuberculosis infection outcomes. Nature. (2018) 560:644–8. doi: 10.1038/s41586-018-0439-x

PubMed Abstract | Crossref Full Text | Google Scholar

13. Gold MC, Cerri S, Smyk-Pearson S, Cansler ME, Vogt TM, Delepine J, et al. Human mucosal associated invariant T cells detect bacterially infected cells. PloS Biol. (2010) 8:e1000407. doi: 10.1371/journal.pbio.1000407

PubMed Abstract | Crossref Full Text | Google Scholar

14. Achkar JM and Casadevall A. Antibody-mediated immunity against tuberculosis: implications for vaccine development. Cell Host Microbe. (2013) 13:250–62. doi: 10.1016/j.chom.2013.02.009

PubMed Abstract | Crossref Full Text | Google Scholar

15. Lund FE. Cytokine-producing B lymphocytes—key regulators of immunity. Curr Opin Immunol. (2008) 20:332–8. doi: 10.1016/j.coi.2008.03.003

PubMed Abstract | Crossref Full Text | Google Scholar

16. Harris DP, Haynes L, Sayles PC, Duso DK, Eaton SM, Lepak NM, et al. Reciprocal regulation of polarized cytokine production by effector B and T cells. Nat Immunol. (2000) 1:475–82. doi: 10.1038/82717

PubMed Abstract | Crossref Full Text | Google Scholar

17. Mauri C and Bosma A. Immune regulatory function of B cells. Annu Rev Immunol. (2012) 30:221–41. doi: 10.1146/annurev-immunol-020711-074934

PubMed Abstract | Crossref Full Text | Google Scholar

18. du Plessis WJ, Keyser A, Walzl G, and Loxton AG. Phenotypic analysis of peripheral B cell populations during Mycobacterium tuberculosis infection and disease. J Inflammation (Lond). (2016) 13:23. doi: 10.1186/s12950-016-0133-4

PubMed Abstract | Crossref Full Text | Google Scholar

19. Tangye SG and Tarlinton DM. Memory B cells: effectors of long-lived immune responses. Eur J Immunol. (2009) 39:2065–75. doi: 10.1002/eji.200939531

PubMed Abstract | Crossref Full Text | Google Scholar

20. Barcelos W, Martins-Filho OA, Guimarães TM, Oliveira MH, Spíndola-de-Miranda S, Carvalho BN, et al. Peripheral blood mononuclear cells immunophenotyping in pulmonary tuberculosis patients before and after treatment. Microbiol Immunol. (2006) 50:597–605. doi: 10.1111/j.1348-0421.2006.tb03834.x

PubMed Abstract | Crossref Full Text | Google Scholar

21. Wu YE, Zhang SW, Peng WG, Li KS, Li K, Jiang JK, et al. Changes in lymphocyte subsets in the peripheral blood of patients with active pulmonary tuberculosis. J Int Med Res. (2009) 37:1742–9. doi: 10.1177/147323000903700610

PubMed Abstract | Crossref Full Text | Google Scholar

22. Corominas M, Cardona V, Gonzalez L, Caylà JA, Rufi G, Mestre M, et al. B lymphocytes and co-stimulatory molecules in Mycobacterium tuberculosis infection. Int J Tuberc Lung Dis. (2004) 8:98–105.

PubMed Abstract | Google Scholar

23. Kozakiewicz L, Phuah J, Flynn J, and Chan J. The role of B cells and humoral immunity in Mycobacterium tuberculosis infection. Adv Exp Med Biol. (2013) 783:225–50. doi: 10.1007/978-1-4614-6111-1_12

PubMed Abstract | Crossref Full Text | Google Scholar

24. Zhang JY, Zhang Y, Yang QT, Deng GF, Chen Q, Chen WX, et al. Dynamic changes of T cell and B cell subsets in peripheral blood of tuberculosis patients during treatment. China Trop Med. (2019) 19:678–82. doi: 10.13604/j.cnki.46-1064/r.2019.07.19

Crossref Full Text | Google Scholar

25. Estévez O, Anibarro L, Garet E, Martínez A, Pena A, Barcia L, et al. Multi-parameter flow cytometry immunophenotyping distinguishes different stages of tuberculosis infection. J Infect. (2020) 81:57–71. doi: 10.1038/s41426-018-0213-z

PubMed Abstract | Crossref Full Text | Google Scholar

26. Yang E, Yang R, Guo M, Huang D, Wang W, Zhang Z, et al. Multidrug-resistant tuberculosis (MDR-TB) strain infection in macaques results in high bacilli burdens in airways, driving broad innate/adaptive immune responses. Emerg Microbes Infect. (2018) 7:207. doi: 10.1038/s41426-018-0213-z

PubMed Abstract | Crossref Full Text | Google Scholar

27. Fonseca KL, Lozano JJ, Despuig A, Habgood-Coote D, Sidorova J, Aznar D, et al. Unravelling the transcriptome of the human tuberculosis lesion and its clinical implications. Nat Commun. (2025) 16:5028. doi: 10.1038/s41467-025-60255-w

PubMed Abstract | Crossref Full Text | Google Scholar

28. Tsai CY, Oo M, Peh JH, Yeo BCM, Aptekmann A, Lee B, et al. Splenic marginal zone B cells restrict Mycobacterium tuberculosis infection by shaping the cytokine pattern and cell-mediated immunity. Cell Rep. (2024) 43:114426. doi: 10.1016/j.celrep.2024.114426

PubMed Abstract | Crossref Full Text | Google Scholar

29. Koyuncu D, Tavolara T, Gatti DM, Gower AC, Ginese ML, Kramnik I, et al. B cells in perivascular and peribronchiolar granuloma-associated lymphoid tissue and B-cell signatures identify asymptomatic Mycobacterium tuberculosis lung infection in Diversity Outbred mice. Infection Immunity. (2024) 92:e0026323. doi: 10.1128/iai.00263-23

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: B-cell subsets, flow cytometry, immune dysregulation, marginal zone-like B cells (MZ B), multidrug-resistant tuberculosis

Citation: Liu M, Wen J, Gao Q, Xu T, Lan Y and Meng L (2026) Peripheral blood B-cell compartment dysregulation in multidrug-resistant tuberculosis is associated with reduced circulating marginal zone-like B cells. Front. Immunol. 17:1709981. doi: 10.3389/fimmu.2026.1709981

Received: 21 September 2025; Accepted: 21 January 2026; Revised: 19 January 2026;
Published: 10 February 2026.

Edited by:

Jianping Xie, Southwest University, China

Reviewed by:

Bo Yan, Fudan University, China
Socorro Miranda-Hernandez, James Cook University, Australia

Copyright © 2026 Liu, Wen, Gao, Xu, Lan and Meng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lu Meng, bWVuZ2x1QGZ1ZGFuLmR1LmNu; Yuanbo Lan, bHliaXZ5QDE2My5jb20=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.