O'Grady, J. F.; Leonard, A. S.; Li, H.; Fang, L.; Pausch, H.; Gormley, I. C.; Gordon, S. V.; MacHugh, D. E.

Understanding the functional impact of genomic sequence variants is critical for evaluating the role of genetic variation in the host response during tuberculosis (TB) disease and anti-TB treatment (ATT). Hitherto, there have been no genome-wide in vivo response expression quantitative trait loci (reQTL) studies conducted for active TB and ATT. Here, using longitudinal peripheral blood RNA-seq data from n = 48 patients with active TB who underwent ATT, we call sequence variants directly from these transcriptomes and impute them with a multi-ancestry reference panel. Associating our variants with the expression of nearby genes, we characterise thousands of cis-eQTL and hundreds of reQTL. We further show significant changes in cell type proportions during ATT through deconvolution of the bulk RNA-seq data and identify the putative cell type specific nature of cis-eQTL. Our work sheds light on the immunogenetics of TB disease and treatment, while providing a framework for studies using only RNA-seq data.