Network analysis identifies regulators of lineage-specific phenotypes in Mycobacterium tuberculosis

The Mycobacterium tuberculosis (Mtb) complex includes ten phylogenetically distinct human-adapted lineages with varying geographical distributions and pathogenicity. Lineage 1 (L1) is associated with low virulence, while Lineage 2 (L2) is linked to hyper-virulence, increased transmission, and drug resistance. We performed multi-layer comparative analyses integrating whole-genome sequencing with quantitative transcriptomic and proteomic profiling of a panel of L1 and L2 clinical strains, each grown under two in vitro conditions. Our data revealed variable correlations between transcript and protein levels across strains and gene categories, indicating lineage-specific post-translational regulation. Transcriptional and translational differences scaled with phylogenetic distance, with one in three SNPs on average leading to gene expression changes. A newly developed genome-scale transcriptional regulatory model identified master transcription factors--closely linked to the sigma factor network--whose targets were differentially expressed between L1 and L2. Notably, DosR proteins showed higher basal levels and a stronger nitric oxide (NO) response in L2. Time-course validation using an Rv1985c induction and wild-type H37Rv under hypoxia and reaeration confirmed the role of Rv1985c in the differentially expressed genes between L1 and L2. These findings suggest that limited genetic variation can translate into significant phenotypic differences via differential regulation of key transcriptional networks. HighlightsO_LIProteomic and transcriptomic characterization of fully sequenced diverse L1 and L2 clinical strains of Mtb. C_LIO_LIPost translational control mechanisms for regulatory and virulence genes are mitigated in Mtb L2. C_LIO_LIA genome-scale transcriptional framework identifies DosR, Rv1985c, Lsr2 and Rv0691c as master transcription factors responsible for differential target gene expression in L2 compared to L1 strains. C_LIO_LIL1 and L2 DosR proteins respond differently to nitric oxide stress, thus determining a relevant phenotype. C_LIO_LITime-course phenotypic and transcriptomic data under hypoxia and reaeration validate Rv1985c as a key transcription factor differentiating L2 and L1 strains. C_LI