Mycobacterium tuberculosis-dependent Monocyte Expression Quantitative Trait Loci and Tuberculosis Pathogenesis

The heterogeneity of outcomes after Mycobacterium tuberculosis (Mtb) exposure is a conundrum associated with millennia of host-pathogen co-evolution. We hypothesized that human myeloid cells contain genetically encoded, Mtb-specific responses that regulate critical steps in tuberculosis (TB) pathogenesis. We mapped genome-wide expression quantitative trait loci (eQTLs) in Mtb-infected monocytes with RNAseq from 80 Ugandan household contacts of pulmonary TB cases to identify monocyte-specific, Mtb-dependent eQTLs and their association with cytokine expression and clinical resistance to tuberculin skin test (TST) and interferon-γ release assay (IGRA) conversion. cis-eQTLs (n=1,567) were identified in Mtb-infected monocytes (FDR<0.01), including 29 eQTLs in 16 genes which were Mtb-dependent (significant for Mtb:genotype interaction [FDR<0.1], but not classified as eQTL in media condition [FDR≥0.01]). A subset of eQTLs were associated with Mtb-induced cytokine expression (n=8) and/or clinical resistance to TST/IGRA conversion (n=1). Expression of BMP6, an Mtb-dependent eQTL gene, was associated with IFNB1 induction in Mtb-infected and DNA ligand-induced cells. Network and enrichment analyses identified fatty acid metabolism as a pathway associated with eQTL genes. These findings suggest that monocyte genes contain Mtb-dependent eQTLs, including a subset associated with cytokine expression and/or clinical resistance to TST/IGRA conversion, providing insight into immunogenetic pathways regulating susceptibility to Mtb infection and TB pathogenesis.


INTRODUCTION
Despite ongoing efforts to eradicate tuberculosis (TB), it remains one of the leading infectious diseases worldwide [1]. The wide spectrum of clinical states, ranging from asymptomatic latent TB infection (LTBI) to active TB disease, presents a challenge in identifying immune events providing protection early after exposure. While various inter-individual and environmental factors can influence TB susceptibility, such as bacillary load, proximity, duration of contact, age, malnutrition, and hygienic conditions, host genetics and immunologic factors are believed to play a significant role in increasing risk for TB disease [2,3]. Therefore, it is critical to identify immunogenetic determinants of TB susceptibility to develop effective prevention strategies and provide insights for new TB treatments.
Mycobacterium tuberculosis (Mtb), the primary causative pathogen of TB, is an ancient infectious agent and has co-evolved with humans for over 10,000 years [4]. Host and pathogen genetic pressure over this long time period selects genes and variants that may regulate protective immune responses. The primary niche of this evolutionary battle occurs within myeloid cells where Mtb resides in a phagosome. Such pressures can contribute to inter-individual variation in susceptibility to Mtb infection and TB disease [5]; however, the specific major susceptibility genes and variants remain largely unknown. The search for the genetic determinants of clinical TB susceptibility has included case-control studies of candidate genes and genome-wide association studies (GWAS) [6][7][8][9][10][11]. Despite extensive efforts, the major genes and specific functional effects of the polymorphisms associated with TB remain poorly understood. One approach to address this knowledge gap is expression quantitative trait loci (eQTL) mapping [12]. By linking genetic variants with Mtb-dependent changes in intermediate traits within myeloid cells, such as cytokine . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) 5 gene expression, eQTL mapping can provide insights into underlying immunogenetic mechanisms of disease.
To investigate the functional consequences of genetic variants within macrophages responding to Mtb infection, we examine our hypothesis that a specific subset of monocyte genes has eQTLs that regulate pro-inflammatory responses to Mtb or that are associated with clinical TB phenotypes. Through genome-wide genotyping and transcriptional profiling of Mtb-infected and uninfected CD14+ monocytes from 80 household contacts (HHCs) of pulmonary TB in Uganda, we identified 29 eQTLs associated with the expression of 16 genes under Mtb-stimulated and not unstimulated conditions. Some of these 'Mtb-dependent' eQTLs were also associated with Mtb DNA sensor-dependent interferon- (IFN-) expression, resistance to tuberculin skin test (TST) and interferon- (IFN-) release assay (IGRA) conversion after Mtb exposure, and modulation of immunometabolic pathways. These findings shed light on the immunogenetic factors underlying TB susceptibility and host-pathogen interactions in the context of Mtb infection.

Monocytes express 16 genes with Mtb-dependent eQTLs.
To determine whether monocyte responses to Mtb include stimulation-specific eQTLs, we examined genome-wide genotyping (MEGA EX genotyping array) and RNAseq transcriptional profiles of cells infected with H37Rv for 6h or media controls in 101 participants in a previously described Ugandan TB household contact study [13][14][15]. We linked the genotypes and transcriptional profiles to define eQTLs among the 80 subjects from whom these datasets were available ( Figure 1). Although this primary study compared individuals who did (LTBI) or did not . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) 6 ("resister" or RSTR) convert their TST/IGRA after heavy exposure, we did not include these phenotypes in our primary analysis. Both the RSTR and LTBI groups were relatively young (mean age = 23 years), with no significant differences observed in sex, body mass index (BMI), Bacille Calmette-Guérin (BCG) vaccination history, or epidemiologic exposure risk score between RSTR and LTBI groups (Table 1) [16]. To discover cis-eQTLs, we considered single nucleotide polymorphisms (SNPs) within 1 megabase (Mb) of the transcription start site (TSS) and their association with gene expression in the media and Mtb conditions. We observed 1,567 cis-eQTLs associated with the expression of 159 genes in Mtb-infected monocytes, whereas in uninfected monocytes, we identified 2,106 cis-eQTLs associated with the expression of 261 genes (false discovery rate [FDR] < 0.01, Table S1). There was a significant overlap of eQTLs identified in infected and uninfected monocytes (n = 1,269; Wilcoxon signed rank test, p < 0.001). To discover "Mtb-dependent" cis-eQTLs, we next evaluated the interaction effect of genotypes between Mtbinfected and uninfected monocytes, while controlling for age and sex (FDR < 0.1). Among the 1,567 eQTLs identified in Mtb-infected monocytes, we found 32 eQTLs associated with the expression of 17 genes that were significant for an Mtb-infection:genotype interaction (FDR <0.1).
Three of these eQTL were also identified under unstimulated conditions (FDR <0.01), and thus were excluded (Figure 2a). For further analysis, we focused on the remaining 29 eQTLs, referred to as "Mtb-dependent eQTLs", which were associated with the expression of 16 genes (Table S2, Figure S1).
Among the 16 Mtb-dependent eQTL genes (Table 2, Table S2, Figure 2b, Figure S1), 3 were associated with multiple eQTLs. Specifically, 12 eQTLs were associated with SIRPB1 expression (Figure 3a), 2 eQTLs with PRKAG2, and 2 eQTLs with NDUFAF4 ( Figure S2). Notably, these . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) 7 eQTLs showed a high linkage disequilibrium (LD, r 2 > 0.8), indicating a potential single causative locus for each gene. Mtb infection resulted in clear genotype-dependent effects as compared to the media condition. To further categorize the Mtb effects, we compared the slope change across genotypes that was observed under unstimulated and Mtb-infected conditions. Some genes (CPSF1 and SIRPB1) had similar directionality of genotype effects in the media and Mtb condition that included a low FDR in the media condition (0.01 to 0.1), indicating that genetic effects are present even without Mtb infection, but they intensify following infection (Figure 3b To investigate the functional roles of 16 Mtb-dependent eQTL genes, we constructed a network of their protein-level interactions using the STRING database [17] and identified two highconfidence edges (STRING score > 700, Figure 4). To gain further insight into the functional roles of these genes, we performed hypergeometric mean enrichment pathway analysis of the 16 Mtbdependent eQTL genes using the Human Molecular Signatures Database (MsigDB v2023.1, Table   . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) 8 S3) [18]. In the C2 canonical pathway, we observed significant enrichment of gene sets related to the fatty acid and glutathione metabolism for GPX2, GSTO2, MLYCD, and PRKAG2 (FDR < 0.1).
Additionally, within the C5 gene ontology biological pathways, we found associations of BMP6, CHMP1B, FANCA, MAS1, MLYCD, PRKAG2, and PIBF1 with 52 gene sets related to fatty acid, lipid, and organic hydroxy metabolism as well as mitotic spindle (Table S3, FDR < 0.1). In summary, our network and enrichment analyses of the 16 Mtb-dependent eQTL genes revealed their involvement in multiple pathways associated with fatty acid metabolism, glutathione metabolism, and other cellular metabolic processes.

Polymorphisms in Mtb-dependent eQTL genes and association with resistance to TST/IGRA conversion in TB household contacts.
We examined the clinical significance of Mtb-dependent eQTL genes and variants by testing their association with resistance to TST/IGRA conversion in highly exposed HHCs of pulmonary TB subjects. We conducted a candidate gene association study with a previously collected extended sample size (RSTR [case], n=74) versus (LTBI [control], n=189) for each Mtb-dependent eQTL gene, adjusting for age, sex, and kinship using a generalized linear mixed model (GENESIS) [11].
One of the Mtb-dependent eQTLs, rs6599528 (associated with the cis-expression of CPSF1), was significantly associated with the RSTR phenotype where the adjusted odds of being RSTR increased by 1.63 times for each copy of the minor allele at this locus (p= 0.04, Table 3). We also examined whether SNPs within the up-and down-stream 5-kilobase (kb) flanking sequences of each of the 16 Mtb-dependent eQTL genes were associated with the RSTR outcome. We identified 26 SNPs in 7 genes associated with clinical RSTR status (p < 0.05; Table S4), the majority of which were found in PRKAG2 (n= 20), with 6 having previously been reported [13,19]. None of . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) 9 these findings remained significant after correction for multiple comparisons. Taken together, our results suggest that Mtb-dependent eQTLs and SNPs in their gene region may be associated with TB clinical outcomes.

Mtb-dependent eQTLs associate with monocyte cytokine expression.
Under inflammatory conditions, genetic polymorphisms of pro-inflammatory cytokines, such as tumor necrosis factor (TNF), interleukin (IL)-1β, IL-6, or type I interferons (IFN-β), have been associated with susceptibility to TB in humans [20-23]. We examined whether Mtb-dependent eQTLs are associated with the expression levels of these cytokines. We identified 8 of the 29 Mtb-dependent eQTLs that were associated with Mtb-induced cytokine expression in monocytes (additive regression model adjusted for age and sex, p < 0.05, Figure 2a). Among these, 4 eQTLs were associated with altered IFNB1 gene expression in response to Mtb infection (p < 0.05, and/or IL1B expression levels between Mtb-infected and uninfected monocytes (p < 0.05, Figure   5d-h, Table S5). We also observed that 6 of the 8 eQTL genes additionally had variants that were associated with clinical RSTR outcomes (Figure 2a). These findings indicate that Mtb-dependent eQTL genes identify host genes that may regulate Mtb-induced cytokine induction in trans and potentially contribute to clinical resistance.
BMP6 expression is associated with Mtb and DNA-ligand induced IFNB expression in monocytes.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)

10
To further test our hypothesis that Mtb-dependent eQTL genes regulate cytokine pathways in monocytes, we used gene silencing methods to investigate potential mechanisms. Among 16 Mtbdependent eQTL genes, we prioritized genes based on three criteria: (i) more pronounced slope changes between Mtb-infected and uninfected conditions, (ii) significant associations with cytokine gene expression, and (iii) significant associations with the clinical RSTR outcomes. From this analysis, BMP6, CHMP1B, KCNK5, PIBF1, PRKAG2, and TIMM44 were initially selected for gene silencing experiments. After assessing siRNA knockdown efficiency, we ultimately focused on BMP6, PRKAG2, and TIMM44 to investigate their role in cytokine signaling experiments.
After silencing expression of BMP6, PRKAG2, and TIMM44, we stimulated PMA-differentiated THP-1 cells with DNA/RNA, toll-like receptor (TLR), or inflammasome ligands. We observed a decrease in IFNB1 expression after 4h of stimulation with 4 μg/mL sheared calf thymus DNA in BMP6-knockdown cells compared to the control group ( Figure 6a). This expression pattern aligns with the HHC primary monocyte data, where the minor allele (G) of rs55966428 was associated with a significant reduction in the expression of both BMP6 ( Figure 6b) and IFNB1 (Figure 6c) after 6h of Mtb infection. Although the other DNA/RNA ligands showed a similar trend of lower IFNB1 expression, the results were not statistically significant. In contrast, we did not find significant differences in IFNB expression for PRKAG2 or TIMM44 ( Figure S3a). There were no significant differences in TNF, IL-6, or IL-1β secretion observed between the BMP6, PRKAG2, and TIMM44-silenced cells and the siRNA controls (Figure 6d-g, Figure S3b-e). Collectively, our data suggests that BMP6 regulates the expression of IFNB1 following Mtb and DNA ligand stimulation, but not TLR-or inflammasome-induced cytokine expression.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted August 29, 2023.

DISCUSSION
Myeloid cells play a critical role to restrict Mtb replication through intrinsic microbicidal pathways and by stimulating downstream cellular responses, but the genetic factors that determine Mtb infection and disease outcomes remains largely unknown. While we previously identified 260 differentially expressed genes in monocytes in response to Mtb infection and clinical TB phenotypes [14], this analysis was limited in fully capturing the genetically encoded monocyte responses following Mtb infection. In the current study, using genome-wide eQTL mapping in Mtb-infected monocytes, we identified 29 Mtb-dependent eQTLs in 16 genes. A subset of these eQTLs is also associated with Mtb-induced cytokine gene expression and/or clinical resistance to TST/IGRA conversion (RSTR phenotype). The expression profiles of these Mtb-dependent eQTL genes and their associated pathways indicate the involvement of diverse inflammatory and cellular metabolic processes, potentially influencing the mechanisms underlying resistance or susceptibility to Mtb infection.
Our primary findings suggest that BMP6 (bone morphogenetic protein 6), a member of the transforming growth factor (TGF)-β superfamily, is involved in the genetic regulation of IFN-β and may influence antimicrobial responses. We also found evidence suggesting Mtb may regulate monocyte gene expression through mRNA processing. CPSF1 (cleavage and polyadenylation specific factor 1) is crucial for mRNA maturation and alternative 3' untranslated region isoform generation [41]. The involvement of CPSF1 in various human diseases, including cancers [42-44] and ocular disorders [45, 46] has been studied, but its role in TB susceptibility remains unexplored. Importantly, our genetic analysis . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted August 29, 2023. ; https://doi.org/10.1101/2023.08.28.23294698 doi: medRxiv preprint 13 revealed that the minor allele (C) at rs6599528 was associated with a 1.63-fold increase in the odds of being RSTR, potentially indicating a protective property of the C allele against TST/IGRA conversion. Upon Mtb infection, the majority of individuals who carry the major allele (A) has decreased expression of CPSF1; while no changes in its expression among individuals with recessive genotype (C/C). These findings suggest that the genotypic regulation of CPSF1 via rs6599528 may potentially mediate alternative splicing of host transcripts that influence susceptibility to Mtb infection. Among potential host pathways that mediate this protection, recent attention has focused on the alternative role of CPSF1 as an E3 ubiquitin ligase that degrades the transcription factor hypoxia-inducible factor (HIF)-1α. HIF-1α is essential for the IFN--driven macrophage immunometabolic response [47, 48] and its activation results in increased IL6 expression [49]. Our findings indicate that individuals with the homozygous recessive genotype (C/C) at this locus, which is linked to increased CPSF1 expression, are predominantly observed in the RSTR group, which is characterized by a lack of response to IFN-γ activation signals (e.g., IGRA positivity). Increased CPSF1 expression among RSTRs, which is predicted to result in HIF-1a degradation, is consistent with IFN-γ-independent protective mechanisms [50, 51].
Interestingly, in our data, IL6 expression is elevated in resting monocytes carrying the minor allele ( Figure S4h). However, the genetic effects of rs6599528 on IL6 expression become non-significant following Mtb infection, indicating that IL-6 regulation is influenced by more complex and multifaceted factors. Therefore, further investigations are needed to elucidate the complex influence of rs6599528 on pro-inflammatory signaling, antimicrobial mechanisms against Mtb, and clinical resistance outcomes.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; https://doi.org/10.1101/2023.08.28.23294698 doi: medRxiv preprint production through oxidative phosphorylation [85]. The convergence of these findings adds weight to the results and encourages further exploration of shared genes or pathways.
This study has several limitations. First, the 'Mtb-dependent eQTLs' identified may not be specific to Mtb infection since our stimulation condition did not include other pathogens or ligands, an area for future investigation. Additionally, our findings using peripheral blood monocytes may not extend to other tissues and cell types at the site of Mtb pathogenesis such as alveolar macrophages where variants or Mtb-dependent eQTL genes may be distinct. Finally, although our siRNA silencing experiments and cytokine gene association analysis suggest a plausible causal effect that Mtb-dependent eQTL gene expression has on cytokine responses following Mtb-infection, future advanced functional assays will explore the regulatory roles of eQTLs' promoter/enhancer elements that mediate these responses.
In summary, this study demonstrates the robust application of genome-wide eQTL mapping to investigate genetic regulation of the host response against Mtb. Our findings highlight significant associations between Mtb-dependent eQTL genes and variants and alterations in the immune response to Mtb infection, as well as potential regulatory interactions between specific eQTLs and Mtb-induced cytokines. These insights into immunogenetic mechanisms provide valuable targets for host-directed therapies and TB control.

Overview of study design and sampling
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. instructions (Qiagen). The RNA samples were quantified using Nanodrop 8000 instrument (Thermo Scientific), and their quality was measured by TapeStation (Agilent) (RNA Integrity Number ≥ 8.0). Next, cDNA libraries were prepared with random hexanucleotide primers and . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; rRNA depletion using SMARTer RNAseq Kit (Takara), followed by sequencing on Illumina HiSeq 2500 and Novaseq 6000 platforms as previously described [14].

Sequence data processing
Sequences were processed as described previously [14]. Briefly, sequences were aligned to the GRCh38 reference genome using STAR 2.6.0 [87] and quantified read counts using RSEM 1.3.0 [88]. We excluded 4,691 genes ( Figure 1) that were non-protein coding, not reliably annotated, or located on a sex chromosome. The latter would have complicated assessment of deviation from Hardy-Weinberg equilibrium (HWE) due to males being hemizygous. This left 13,439 genes, which we then trimmed-mean of M-values (TMM) normalized and converted to log2 counts per million (CPM) using limma [89]. In parallel, genome-wide genotyping of 263 HHCs was performed using the Illumina MEGA EX array. We successfully genotyped 1,002,430 SNPs, of which we filtered out 731,662 SNPs that deviated from HWE (P < 10 -6 ) and had a minor allele frequency (MAF) of less than 0.05.

eQTL analysis
To identify genetic variants that regulate mRNA expression in a cis-acting manner, we tested for associations between transcript expression levels and genotypes at SNPs located within a 1Mb window centered on the genes' transcription start sites (TSS). As we lacked substantial power for trans analysis due to the small sample size, we focused solely on cis-eQTLs. Initially, we performed an additive linear regression of genotypes on log-transformed gene expression level, adjusted for age and sex in either Mtb-infected or uninfected monocytes. The R package MatrixEQTL was used to perform all regressions [90]. We estimated FDR by using a Fisher's . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; https://doi.org/10.1101/2023.08.28.23294698 doi: medRxiv preprint 20 exact test to correct for multiple testing by the Benjamini-Hochberg method to control for false positives. Significant eQTLs were defined as those with an FDR < 0.01 in Mtb-infected or uninfected monocytes.
For the 1,567 eQTLs identified in Mtb-infected monocytes, we further introduced an interaction term (Mtb-infection:genotype) in the main additive regression model to uncover whether the genetic effect of eQTLs on the level of gene expression differed by Mtb infection status. We applied a loose cut-off of FDR < 0.1 in the interaction model, given the small number of eQTLs included in this analysis [91]. From the 32 eQTLs that showed significance in Mtb-infected condition and the interaction model, we removed 3 eQTLs that overlapped with those identified in uninfected condition to capture only those that regulate target genes in response to Mtb infection, which we named "Mtb-dependent eQTLs". Finally, we defined the lead eQTLs per gene as the SNP that was most significantly associated (with the lowest FDR value) with the expression of that gene (Figure 2).

Mapping cytokine expression levels
The associations between genotype and cytokine expression were assessed for each of the 16 lead Mtb-dependent eQTLs. Specifically, we examined the expression levels of TNF, IL6, IL1B, and IFNB1 in Mtb-infected and uninfected monocytes, as genetic variation in these cytokine genes has previously been linked to TB susceptibility [20-23]. We used a linear regression model to regress the number of minor alleles on the differences in cytokine expression between Mtb-infected and uninfected monocytes (normalized log2 [Mtb-infected -uninfected relative expression]), while adjusting for age and sex, at a significance level of 0.05.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ;

Mtb-dependent eQTL gene association study
We conducted association analyses to determine whether a particular eQTL genotype co-occurs with a RSTR clinical phenotype more often than would be expected by chance. We calculated odds ratios (ORs), adjusted for age and sex as well as kinship using a quasi-likelihood

siRNA knockdown of selected Mtb-dependent eQTL genes
We used siRNA to silence Mtb-dependent eQTL genes, and 300 nM of two oligos per gene (siRNA identification s2032 and s2033 for BMP6; s3275a and s32752 for CHMP1B; s16450 and s16451 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ;

23
After 24h of treatment with BMP6 siRNA, cells were washed and replated with fresh RPMI-10 to each well, and treatment with DNA/RNA ligands was performed as described below. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; https://doi.org/10.1101/2023.08.28.23294698 doi: medRxiv preprint 24 For repeated BMP6 knockdown experiments, RNA was isolated from lysates using the same protocol with RNeasy Plus Mini Kit (Qiagen). Synthesis of the first strand cDNA was performed using SuperScript II reverse transcriptase and oligo (dT) primer (Invitrogen). qPCR was performed with the CFX96 real-time system (Bio-Rad) using the SsoFast EvaGreen Supermix with the LOW ROX kit (Bio-Rad). The following primers designed from PrimerBank were used. The PrimerBank identifications are BMP6 (133930782c1) and IFNB1 (50593016c1). CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Study Approval
The current prospective cohort is part of the larger Kawempe Community Health Study conducted in Kampala, Uganda. Detailed information on the original study's setting, recruitment procedure,

Data Availability
Access to raw transcriptomic data is available through the NCBI database of Genotypes and Phenotypes (dbGaP) Data Browser (https://www.ncbi.nlm.nih.gov/gap/) under accession 002445.v1.p1, but first must be approved by data access committees (DACs) for each study site (see Supplemental Methods [14]). All R code is available at https://github.com/hawnlab/RSTR_eQTL_public. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023.

8(5).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 14574.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ;

Figure 1. Study Design and eQTL Mapping Data Processing in the Uganda Household
Contact Study. Household contacts (HHCs) of pulmonary TB index cases were studied, enrolled, and followed longitudinally for assessment using serial TST and IGRA testing. Laboratory studies included RNASeq measurements in Mtb-infected monocytes and genomewide genotyping with the MEGA EX genotyping array.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; https://doi.org/10.1101/2023.08.28.23294698 doi: medRxiv preprint Using an additive regression model (FDR< 0.01), we identified 2,106 eQTLs in uninfected monocytes and 1,567 eQTLs in Mtb-infected monocytes after adjusting for age and sex. Mtb-dependent eQTLs were specifically defined as cis-eQTLs that showed strong statistical evidence only in Mtb-infected monocytes. We performed additive linear regression with an interaction term for the Mtb stimulation status to evaluate the interaction effect at an FDR of 10%. Out of the 2,106 eQTLs in Mtb-infected monocytes, 32 eQTLs showed significant linear associations with gene expression levels after adjusting for age and sex (FDR <0.1). Three eQTLs were excluded as they were identified in both Mtb-infected and uninfected monocytes. We identified 29 Mtb-dependent eQTLs regulating the expression levels of 16 genes. Additionally, 8 Mtb-dependent eQTLs were associated with altered levels of TNF, IL6, IL1B, or IFNB1 expression in trans (orange circle). Finally, 7 Mtb-dependent eQTL genes contained SNPs associated with the clinical RSTR phenotype (blue circle). Importantly, 6 Mtb-dependent eQTL genes and variants were associated with both clinical RSTR phenotypes and the regulation of cytokine gene expression. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ;  is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ;

Figure 6. BMP6 expression is associated with host genotype in Mtb-infected cells and is required for maximal IFNB1 responses in monocytes following DNA-ligand stimulation.
(a) The fold change expression of IFNB1 was significantly reduced in BMP6-silenced THP1 cells after 4h of stimulation with 4 μg/mL calf thymus DNA ligand, compared to siRNA control cells. BMP6-silenced cells stimulated with 1 μg/mL of super coiled plasmid DNA, 20 μM of cGAMP, or 10 μg/mL of poly(I:C) exhibited statistically non-significant decreases in IFNB1 induction compared to siRNA control cells. The fold change expression of IFNB1 was measured by qPCR and normalized against the IFNB1 expression from lipofectamine without DNA/RNA ligand to control background IFNB1 induction from lipofectamine. For TLR-specific stimulation, BMP6-silenced and siRNA control THP-1 cells were stimulated with LPS, PAM2/PAM3, Mtb whole cell lysate, or media for 24h. (b) In the Ugandan HCC cohort, the minor allele (G) of rs55966428 was associated with decreased expression of BMP6 following 6h of Mtb infection in primary monocyte data from HHCs. Host genotype (x-axis) and normalized log2 expression (y-axis) of BMP6 for the rs55966428 Mtb-dependent eQTL with false discovery rate (FDR) reported separately for the media and Mtb-stimulated conditions. (c) In the same cohort, the minor allele (G) of rs55966428 was associated with decreased IFNB1 expression following 6h of Mtb infection in primary monocyte data. This plot shows genotype on the x-axis and differences in IFNB1 expression on the y-axis (normalized log2 [Mtb-infected -uninfected relative IFNB1 expression]), adjusting for age and sex. (d) TNF, (e) IL-6, and (f) IL-1β supernatant levels were measured using ELISA. To assess (g) inflammasome-mediated IL-1β response, nucleofected cells were primed with LPS for 2h and treated with nigericin for NLRP3specific stimulation, Burkholderia thailandensis needle protein for 4h for NLRC4-specific stimulation. Needle protein was administered at 8 ng/ml in conjunction with 16 ng/ml Bacillus anthracis protective antigen. The IL-1β supernatant levels were measured by ELISA. **p < 0.01 tested by paired t-test, ns: not significant.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; Age at enrollment was used to classify individuals into two groups: children, defined as those under 15 years old, and adults, defined as those aged 15 and above. The epidemiologic exposure risk score was calculated based on the intensity of exposure to the index case during enrollment (phase 1, between 2002 and 2012). The score ranged from 0 to 10 for adults and 0 to 9 for children based on their age at enrollment. Additionally, PBMCs were collected, and RSTR status was measured during the first visit of the retracing study (phase 2, between 2014 and 2017) for each participant.
*Mean and median estimates were compared using the parametric t test and nonparametric Wilcoxon-Mann-Whitney test, respectively. The χ 2 was used to test for association between RSTR status and sex. BCG = Bacille Calmette-Guérin, BMI = body mass index, LTBI = latent tuberculosis infection, p = p-value, PBMC = peripheral blood mononuclear cell, RSTR = resister or resistant to TST/IGRA conversion after high TB exposure, SD = standard deviation.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted August 29, 2023. ; a. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.