# Metagenomic and metatranscriptomic profiling of bronchoalveolar lavage fluid identifies microbial and host biomarkers of drug-resistant tuberculosis

**Authors:** Haiqing Zhang, Limao Zhang, Bin Yang, Chunjing Gao, Huimei Liu, Yunshi Zhang, Xiaoyou Chen

PMC · DOI: 10.3389/fcimb.2025.1726935 · Frontiers in Cellular and Infection Microbiology · 2026-01-29

## TL;DR

This study uses metagenomic and transcriptomic data from lung fluid to identify microbial and immune markers that distinguish drug-resistant from drug-susceptible tuberculosis.

## Contribution

The study introduces a dual-omics approach to uncover DR-TB-specific microbial and host signatures, including an 8-gene classifier for accurate discrimination.

## Key findings

- DR-TB airways show distinct microbial diversity with enriched Streptococcus spp. and specific phages.
- DR-TB is associated with 494 differentially expressed genes linked to metabolic and immune pathways.
- An 8-gene host signature accurately differentiates DR-TB from DS-TB with an AUC of 0.837.

## Abstract

Drug-resistant tuberculosis (DR-TB) undermines global TB control, yet how resistant Mycobacterium tuberculosis strains interact with the lung microbiome, phage communities, and local host immunity remains poorly defined.

In a prospective cohort of 130 pulmonary TB patients (49 DR-TB, 81 drug-susceptible TB [DS-TB] patients), bronchoalveolar lavage fluid (BALF) was subjected to paired metagenomic and transcriptomic profiling. Microbial and bacteriophage community structures were assessed by diversity metrics and differential abundance testing, whereas host responses were characterized by gene expression, pathway enrichment, and immune cell deconvolution. A Random Forest model was trained to evaluate the diagnostic potential of host transcriptional signatures.

DR-TB airways presented distinct microbial beta diversity, with enrichment of Streptococcus spp. and streptococcal-targeting phages (e.g., Javan variants, phi-Ssu5SJ28rum). Transcriptomic analysis revealed 494 differentially expressed genes, which were associated with increased oxidative phosphorylation, suppressed ion channel and transporter activity, and enrichment of extracellular matrix remodeling pathways. Immune profiling demonstrated a significant reduction in γδ T cells in DR-TB patients (P = 0.0059). An 8-gene host-derived signature (ARHGEF5, PTGES3L, GAL3ST1, RANBP17, ACTA2_AS1, CBY3, MAMSTR, and LOC102031319) discriminated DR-TB from DS-TB with high accuracy (AUC = 0.837).

This dual-omics study defines the airway niche of DR-TB as a convergence of microbial dysbiosis, phage imbalance, and host immune–metabolic dysfunction. By uncovering DR-TB–specific microbial and transcriptional signatures, and deriving a predictive host-based classifier, our findings provide mechanistic insights and highlight novel opportunities for microbiome- and host-directed interventions in drug-resistant tuberculosis.

## Linked entities

- **Genes:** ARHGEF5 (Rho guanine nucleotide exchange factor 5) [NCBI Gene 7984], PTGES3L (prostaglandin E synthase 3 like) [NCBI Gene 100885848], GAL3ST1 (galactose-3-O-sulfotransferase 1) [NCBI Gene 9514], RANBP17 (RAN binding protein 17) [NCBI Gene 64901], ACTA2-AS1 (ACTA2 antisense RNA 1) [NCBI Gene 100132116], CBY3 (chibby family member 3) [NCBI Gene 646019], MAMSTR (MEF2 activating motif and SAP domain containing transcriptional regulator) [NCBI Gene 284358]
- **Diseases:** tuberculosis (MONDO:0018076), drug-resistant tuberculosis (MONDO:0041806)
- **Species:** Mycobacterium tuberculosis (taxon 1773)

## Full-text entities

- **Genes:** EPC1-AS1 (EPC1 antisense RNA 1) [NCBI Gene 102031319], GAL3ST1 (galactose-3-O-sulfotransferase 1) [NCBI Gene 9514] {aka CST}, RANBP17 (RAN binding protein 17) [NCBI Gene 64901], CBY3 (chibby family member 3) [NCBI Gene 646019], MAMSTR (MEF2 activating motif and SAP domain containing transcriptional regulator) [NCBI Gene 284358] {aka MASTR}, ACTA2-AS1 (ACTA2 antisense RNA 1) [NCBI Gene 100132116] {aka UC001kfo, ZXF1, uc001kfo.1}, ARHGEF5 (Rho guanine nucleotide exchange factor 5) [NCBI Gene 7984] {aka GEF5, P60, TIM, TIM1}, PTGES3L (prostaglandin E synthase 3 like) [NCBI Gene 100885848]
- **Diseases:** DR-TB (MESH:D018088), metabolic dysfunction (MESH:D008659), DS-TB (MESH:D014390)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mycobacterium tuberculosis (species) [taxon 1773], Bacteriophage sp. (species) [taxon 38018]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12894345/full.md

## Figures

6 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12894345/full.md

## References

59 references — full list in the complete paper: https://tomesphere.com/paper/PMC12894345/full.md

---
Source: https://tomesphere.com/paper/PMC12894345