Fibroblast WGCNA (hdWGCNA). (A) High-dimensional WGCNA (hdWGCNA) for gene module detection in fibroblasts of this study. UMAPs are colored by eigengene of each of the seven modules. (B) Evaluation of M1 module expression in fibroblast subclusters. Bonferroni-adjusted P computed from two-sided Wilcoxon test is shown. (C) ORA by enricher on all assigned M1 module genes using MSigDB Gene Ontology Biological Processes (GOBP) gene set database. (D) Top: Bacterial burden of NHP lung granulomas by Gideon et al. (2022) grouped by the time point. Bottom: Evaluation of human TB-myofibroblast module expression in NHP TB fibroblasts on 4-wk and 10-wk samples. Two-sided Mann–Whitney U test without correction was used. Statistical annotations: P value < 0.05 (*), P value < 0.01 (**), and P value < 0.001 (***). (E) Evaluation of human TB-myofibroblast module expression in fibroblasts from granuloma vs uninvolved lungs in an independent NHP study with 4-wk post-infection (p.i) macaques (Bromley et al., 2024) (Materials and methods). Two-sided Mann–Whitney U test without correct was used. Statistical annotations: P value < 0.01 (**) and P value < 0.0001 (****).