Search

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte- derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a proteincoding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Mycobacterium tuberculosis (M.tb) results in 10 million active tuberculosis (TB) cases and 1.5 million deaths each year, making it the world's leading infectious cause of death. Infection leads to either an asymptomatic latent state or TB disease. Memory T cells have been implicated in TB disease progression, but the specific cell states involved have not yet been delineated because of the limited scope of traditional profiling strategies. Furthermore, immune activation during infection confounds underlying differences in T cell state distributions that influence risk of progression. Here, we used a multimodal single-cell approach to integrate measurements of transcripts and 30 functionally relevant surface proteins to comprehensively define the memory T cell landscape at steady state (i.e., outside of active infection). We profiled 500,000 memory T cells from 259 Peruvians > 4.7 years after they had either latent M.tb infection or active disease and defined 31 distinct memory T cell states, including a CD4+CD26+CD161+CCR6+ effector memory state that was significantly reduced in patients who had developed active TB (OR = 0.80, 95% CI: 0.73-0.87, p = 1.21 x 10-6). This state was also polyfunctional; in ex vivo stimulation, it was enriched for IL-17 and IL-22 production, consistent with a Th17-skewed phenotype, but also had more capacity to produce IFNgamma than other CD161+CCR6+ Th17 cells. Additionally, in progressors, IL-17 and IL-22 production in this cell state was significantly lower than in non-progressors. Reduced abundance and function of this state may be an important factor in failure to control M.tb infection. ### Competing Interest Statement The authors have declared no competing interest.