Interred mechanisms of resistance and host immune evasion revealed through network-connectivity analysis of <i>M. tuberculosis</i> complex graph pangenome

<i>Mycobacterium tuberculosis</i> complex successfully adapts to environmental pressures through mechanisms of rapid adaptation which remain poorly understood despite knowledge gained through decades of research. In this study, we used 110 reference-quality, complete <i>de novo</i> assembled, long-read sequenced clinical genomes to study patterns of structural adaptation through a graph-based pangenome analysis, elucidating rarely studied mechanisms that enable enhanced clinical phenotypes offering a novel perspective to the species' adaptation. Across isolates, we identified a pangenome of 4,325 genes (3,767 core and 558 accessory), revealing 290 novel genes, and a substantially more complete account of difficult-to-sequence <i>esx/pe/pgrs/ppe</i> genes. Seventy-four percent of core genes were deemed non-essential <i>in vitro</i>, 38% of which support the pathogen's survival <i>in vivo</i>, suggesting a need to broaden current perspectives on essentiality. Through information-theoretic analysis, we reveal the <i>ppe</i> genes that contribute most to the species' diversity-several with known consequences for antigenic variation and immune evasion. Construction of a graph pangenome revealed topological variations that implicate genes known to modulate host immunity (<i>Rv0071-73, Rv2817c, cas2</i>), defense against phages/viruses (<i>cas2, csm6,</i> and <i>Rv2817c-2821c</i>), and others associated with host tissue colonization. Here, the prominent trehalose transport pathway stands out for its involvement in caseous granuloma catabolism and the development of post-primary disease. We show paralogous duplications of genes implicated in bedaquiline (<i>mmpL5</i> in all L1 isolates) and ethambutol (<i>embC-A</i>) resistance, with a paralogous duplication of its regulator (<i>embR</i>) in 96 isolates. We provide hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can escape detection by molecular diagnostics.IMPORTANCE<i>M. tuberculosis</i> complex (MTBC) has killed over a billion people in the past 200 years alone and continues to kill nearly 1.5 million annually. The pathogen has a versatile ability to diversify under immune and drug pressure and survive, even becoming antibiotic persistent or resistant in the face of harsh chemotherapy. For proper diagnosis and design of an appropriate treatment regimen, a full understanding of this diversification and its clinical consequences is desperately needed. A mechanism of diversification that is rarely studied systematically is MTBC's ability to structurally change its genome. In this article, we have <i>de novo</i> assembled 110 clinical genomes (the largest <i>de novo</i> assembled set to date) and performed a pangenomic analysis. Our pangenome provides structural variation-based hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can compromise molecular diagnostics and lead to further emergence of antibiotic resistance.