History of tuberculosis disease is associated with genetic regulatory variation in Peruvians
============================================================================================

* Sara Suliman
* Victor E. Nieto-Caballero
* Samira Asgari
* Kattya Lopez
* Sarah K. Iwany
* Yang Luo
* Aparna Nathan
* Daniela Fernandez-Salinas
* Marcos Chiñas
* Chuan-Chin Huang
* Zibiao Zhang
* Segundo R León
* Roger I Calderon
* Leonid Lecca
* Megan Murray
* Ildiko Van Rhijn
* Soumya Raychaudhuri
* D. Branch Moody
* Maria Gutierrez-Arcelus

## Abstract

A quarter of humanity is estimated to be latently infected with *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 dendritic cells (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. Five genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding 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.

Keywords
*   tuberculosis
*   expression quantitative trait loci
*   dendritic cells
*   macrophages

## Introduction

Tuberculosis disease (TB), caused by infection with *Mycobacterium tuberculosis* (*Mtb*), is a leading cause of death from infection globally1. Notably, only 5-10% of *Mtb*-infected individuals are estimated to develop active TB, and these individuals are at higher risk of recurrent TB, suggesting the existence of durable host factors that influence disease outcome2. Unbiased systems biology approaches have uncovered several host pathways associated with TB progression, including interferon signaling3, 4, metabolic dysregulation5, and depletion of Th17 cells in the peripheral blood6. Genetic association studies from our group7, 8, and others9–13, identified host genetic variants associated with higher risk of active TB. However, these studies did not systematically integrate the association between genetic variation and transcriptional profiles of myeloid cells, which are the first line of defense following *Mtb* infection.

Growing evidence supports critical roles for myeloid cells in progression to TB disease in human cohorts3, 4, 14–17, where macrophages are the main target of infection18, and dendritic cells present mycobacterial antigens to prime T cells19. Monocytes can be differentiated *in vitro* to generate monocyte-derived DCs20 and macrophages21, which have become useful tools to study innate responses to pathogens. For instance, monocyte-derived DCs from TB susceptible individuals showed elevation of autophagy associated genes, including Fasciculation And Elongation Protein Zeta 2 (*FEZ2*), a repressor of autophagosome maturation22, and group I CD1 genes, antigen presenting molecules that mark DC maturation13, 23. Hence, these *in vitro* models can be used to study the impact of genetic variants associated with differential gene expression that may underlie TB susceptibility or disease-induced regulatory effects24, 25.

Multiple studies mapping genetic variation that influence gene expression levels as expression quantitative trait loci (eQTL) have highlighted the importance of genetic regulatory variation in disease8, 26. These genetic regulatory effects are often context-dependent; that is, the same genetic variant can affect gene regulation to varying degrees depending on the cell type or cell state27–31. For example, eQTL can be present in some human tissues but not others26, 27, they can also vary between subsets of a cell type such as between naïve and memory B cells32. They can also vary between cell states of a given cell type, such as resting and activated monocytes33. It is becoming more evident that identifying these cell state-dependent genetic regulation events are crucial for understanding disease mechanisms34, 35. Multiple studies have shown that *in vitro* pathogen exposure changes cell states and consequently changes the regulatory landscape, revealing pathogen-dependent genetic regulatory effects31, 36–39. Collectively, pathogens, such as *Mtb*, could exploit the cell type-specific impacts of host genetic variation on gene expression to establish virulence or mediate long term pathogenic effects in diverse populations40. Hence, integration of genetic variation with transcriptional differences in myeloid cells may provide insights into disease mechanisms, as supported by prior studies in monocyte-derived DCs22. These genetically mediated transcriptional changes may control host immune responses in ways that change *Mtb* infection outcomes, as well as susceptibility to long term sequalae of infection24.

In this study, we re-enrolled patients in Peru that were rigorously defined as index TB patients and their *Mtb-*infected household-contacts who either progressed or did not progress to TB disease6–8, 41, 42. We profiled gene expression in monocyte-derived DCs and macrophages from these previously genotyped participants7. We identified genetic variants that differentially regulated transcription in DCs and macrophages in progressors compared to non-progressors. The top regulatory event was mediated by a variant on chromosome 15 that specifically regulated expression of fumarylacetoacetate hydrolase (*FAH*) in DCs from progressors, but not non-progressors. This cell type-specific regulation of *FAH* in recovered TB patients provides a biologically plausible metabolic determinant of TB disease outcome.

## Methods

### Human Participants

Participants were nested and re-recruited from the original prospective cohort of index TB patients and their household contacts42 **(Figure 1A)**. Briefly, household contacts in the original study were screened for symptoms of TB disease. All household contacts received an intradermal tuberculin skin test at enrollment, which was read 48-72 hours later. Cases in the current study were defined to include both former HIV-negative index TB patients, with pulmonary TB and microbiologically confirmed *Mtb* culture (primary cases), as well as household contacts of index TB patients who progressed to active tuberculosis disease within 12 months of contact (secondary cases). All cases received antibiotic treatment according to the standard of care per Peruvian national guidelines.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F1)

Figure 1: Myeloid cell transcriptomic study in a Peruvian cohort reveals genetic variants associated with *FAH* transcripts in dendritic cells for subjects with history of TB disease but not for household contacts exposed to *Mtb* that did not progress to active disease.
**(A)** Cohort scheme: human participants included former household contacts of TB patients who either progressed (cases, progressors) or did not progress to active TB disease but have converted their tuberculin skin test within 12 months of exposure (non-progressors, controls). Genotyping data were previously obtained from all participants. **(B)** Transcriptomic study scheme: Monocytes were magnetically enriched from peripheral blood mononuclear cells and differentiated for 6 days into either monocyte-derived dendritic cells (GM-CSF + IL4) or macrophages (M-CSF). Samples were analyzed by low-input RNA-sequencing, and genotyping data were integrated with transcriptional profiles in an eQTL study. Icons were generated in Biorender. (C-D) *FAH* expression levels are stratified by genotype of indel rs142312981. P-value shown is for the TB status interaction with the genetic variant on their effects on gene expression (case=TB progressor, control= non-progressor) in DCs **(C)** and macrophages **(D)**, details of the linear regression and likelihood ratio test are described in Methods. **(E-F)** Figure shows re-analysis of public data from Pacis et al. (PMID: 30886108, GSE116405). (E) *FAH* expression levels are shown from RNA-seq data from uninfected (blue), *Mtb-*infected (red), heat-killed *Mtb* treated (green) monocyte derived DCs using a mixed-effects linear model. **(F)** Figure shows average mean change in DNA methylation levels between *Mtb*-infected and non-infected monocyte-derived DCs from 6 donors is shown. Data for CpG sites within 250kb (upstream and downstream) of *FAH* transcriptional start site (TSS) are shown.

Controls were defined as former household contacts of TB patients who had a positive tuberculin skin test result during 12 months of follow up. Importantly, this group did not develop active TB disease between the initial recruitment into the original longitudinal prospective study7, 42 and re-recruitment into the current nested case-control study. Genotyping data from all individuals in this study have been previously generated for a genome-wide association study7 and integrated with transcriptomic data as outlined below. The Institutional Review Board of the Harvard Faculty of Medicine and Partners Healthcare (protocol number IRB16-1173), and the Institutional Committee of Ethics in Research of the Peruvian Institutes of Health approved this study protocol.

### Samples

Peripheral blood mononuclear cell (PBMC) samples were separated using standard Ficoll density gradient centrifugation from 50mL of venous blood and cryopreserved at 5X106 cells/vial. Samples were shipped to the Brigham and Women’s Hospital for storage in liquid nitrogen and subsequent experiments.

### Differentiation of monocyte-derived dendritic cells and macrophages *in vitro*

PBMC samples were thawed, and CD14-expressing monocytes were magnetically sorted using the Miltenyi pan monocyte isolation kit, humans (Miltenyi) following the manufacturer’s instructions in magnetic-associated cell sorting (MACS) buffer: 0.5% bovine serum albumin, 2mM ethylene-diamine-tetra-acetic acid (EDTA) in 1X phosphate-buffered saline (PBS). To generate either monocyte-derived dendritic cells or macrophages, monocytes were washed and resuspended in supplemented complete RPMI (5% fetal calf serum, 1mM 2-mercaptoethanol, penicillin-streptomycin, L-Glutamine, HEPES buffer, both essential and non-essential amino acids, and sodium hydroxide). The media was either supplemented with 300 Units/mL of granulocyte-macrophage colony stimulating factor (GM-CSF) and 200 Units/mL of interleukin (IL)-4 to generate monocyte-derived dendritic cells (DCs), or 20ng/mL of macrophage colony-stimulating factor (M-CSF) to generate monocyte-derived macrophages. For both cell types, monocytes were differentiated in 6-well tissue culture plates for 6 days, then split between analytical flow cytometry to test viability and differentiation, and cell sorting to generate pure populations for RNA extraction.

### Cell sorting and flow cytometry analysis

Magnetically isolated CD14 cells treated with cytokines for 6 days were sorted by gating on forward and side scatter profiles after differentiation to remove subcellular debris directly into the TCL RNA lysis buffer (Qiagen). Only samples with at least 1000 cells were included in the final RNA-sequencing analysis. In parallel, to confirm viability and differentiation, cultured cells were stained first with a fixable blue viability cell stain (ThermoFisher Scientific) according to manufacturer’s instructions, followed by different cocktails of fluorochrome-conjugated antibody cocktails for each cell type **(Supplementary Table 1)**, and analyzed by analytical flow cytometry on a 5-laser Fortessa (BD Biosciences) to confirm viability and differentiation. Flow cytometry panels were optimized by titrating individual antibodies to find optimal signal-to-noise ratios, and minimal spillover from spectral overlap of emissions using fluorescence minus-one. Flow cytometry (FCS3) files were analyzed using FlowJo v9.9.

### RNA sequencing

RNA-seq library preparation was performed at the Broad Technology Labs (Broad Institute) using a modified Smart-seq243 protocol for low-input RNA-sequencing, which improves detection of full-length transcripts using template switching and pre-amplification. There were 3 plates of 96 samples each **(Supplementary Table 2)** including blank samples. Samples from cases and controls and the two cell-types were randomized across all three plates. Sequencing was performed in an Illumina NextSeq500 as paired-end 2 x 25 bp reads with an additional 8 cycles for each index.

### Analysis of RNA-sequencing data

Read alignment was performed using STAR (v2.4.2)44 with the following parameters: -- twopassMode Basic, -- alignIntronMax 1000000, --alignMatesGapMax 1000000, --sjdbScore 2, \- --|quantMode TranscriptomeSAM, --sjdbOverhang 24. RSEM (v1.2.21)45 was then used for gene quantification with paired-end mode. We used the hg38 reference genome (University of California Santa Cruz Genome Browser) and GENCODE annotation version 24 (Ensembl 83). Quality summary statistics were gathered using Picard tools. Data quality was evaluated using RNA-SeQC46. RNA-seq data analysis was performed with R (version 3.6.0) and R Bioconductor. Samples with low (<80%) proportion of common genes detected (common genes defined as those present in >75% of samples) or presenting as outliers in principal component analysis (performed on top 1000 variable genes based on standard deviation) were filtered out. Five samples that did not match their genotype information (based on proportion of both alleles seen in RNA-seq data over heterozygous sites) were filtered out. Differential expression of genes by TB status was determined using a multivariable linear regression model adjusting for plate (batch), age, and gender, using the R limma package47 and P-values were adjusted using the Benjamini-Hochberg false discovery rate method.

### Genotype data

We used genotyping data from the customized Affymetrix LIMAArray genome-wide genotyping chip as described7, where QC and imputation was performed. Briefly, individuals were filtered out if they presented higher than 5% missing genotyping data, or had excess of heterozygous genotypes, or had >40 years at age of diagnosis. Variants were filtered out if they fell in either of these criteria: call rate <95%, duplicated based on coordinates, presentation of batch effect, deviation from Hardy-Weinberg equilibrium (HWE, P < 10-5), large missing rate differences between cases and controls. After these filters, the 677,232 SNPs left were used for imputation with SHAPEIT248 and IMPUTE249. After requiring imputation score r2 >= 0.4 and re-filtering based on HWE and missing rate, there were 7,756,401 genetic variants. For eQTL analyses, we removed SNPs with minor allele frequency (MAF) <0.09 to ensure presence of at least 10 minor alleles in tested genes, and avoid problems due to outliers, especially for the interaction analyses. We reported all genomic positions using Grch37 coordinates.

### Expression quantitative trait loci (eQTL) analysis

In total, 126 individuals (with 118 and 109 RNA-seq samples, corresponding to DCs and macrophages respectively) had high data quality for both genotype and gene expression for cis-eQTL mapping. We performed a separate analysis for each cell type. Genes expressed at low levels where log2(transcripts per million (TPM)+1) <1 in more than 30 samples, were filtered out in each cell subset. We restricted the analysis to SNPs that were within 1Mb of the gene start coordinate. QTLtools50 was run in the permutation pass mode (1000 permutations for top SNP per gene) to identify associations between genetic variants and quantile normalized gene expression levels expressed as log2(TPM+1) using linear regression. To control for covariates, we included 5 genotyping principal components (PCs), 15 expression PCs, age, and sex.

To screen for potential interactions of genetic variants with TB status, we followed a two-step nested approach that was recently shown to improve the discovery of gene-by-environment interactions51. First, we selected SNP-gene pairs with a suggestive eQTL (p<0.001), where only top SNP per gene was selected. Next, we applied Levene’s test (LT) to prioritize SNPs that were significantly associated with the variance in gene expression at 5% FDR using the qvalue R package52. Finally, we used these variance quantitative trait loci (vQTLs) to test for an interaction between the SNP and TB status on gene expression. Specifically, we tested for the SNP-TB status interaction effects by performing a likelihood ratio test between two nested models using the R anova function to conduct an analysis of variance (ANOVA) test. The null model estimates the effects of the SNP, TB status, and covariates evaluated using 5 genotyping PCs, 15 expression PCs, age, and sex) on a gene’s log2(TPM+1) normalized expression. The alternative model has an additional SNP by TB status interaction term: ![Formula][1]</img>  where *Ei,j* is the normalized gene expression for the *i*th sample form the *j*th subject, *θ* is the intercept, *βgeno*’ is the effect (eQTL) of the genotype for subject *j* (*gj*), *βtb* is the effect of the TB status (cases vs. controls) for subject *j* (*gj*), and *βx* is the effect of the TB status by genotype interaction (*tj* · *gj*).

As before, we used the qvalue package to calculate the FDR across tests and called all events with q<0.1 (<10% FDR) as significant. Given our limited sample size, we sought to rigorously determine if the significant events were driven by potential outlier events, so we performed permutations by shuffling the TB status for each individual 1000 times and re-computing the SNP by TB status interaction.

To identify cell-type specific eQTL interactions, we first identified the best SNP per gene for expression using both cell types together with QTLtools. For the covariates’ matrix, we removed the first expression PC because it separated the samples according to cell type. We then followed the same two-step nested approach described above with the cell type as an interaction term to identify significant SNP-cell type interactions at 10% FDR.

### Concordance between eQTLs

To verify the agreement between the mapped eQTLs in our study and the previous data, we selected the variant-gene pairs that were significant at FDR <0.05 in each cell type and compared the effect size magnitude (beta from QTLtools results) with the effect size in the other cell type. Additionally, we compared the beta of the most associated variant for each of the genes with an eQTL (FDR <0.05) in macrophages with the corresponding effect size of that variant-gene pair reported in the eQTL catalogue53 from a previously published dataset39.

### Stimulation with irradiated *Mtb* and real time PCR

We obtained de-identified leukoreduction filters (leukopak) samples from healthy blood bank donors through The Brigham and Women’s Hospital Specimen Bank, as approved by the Institutional Review Board of MassGeneralBrigham HealthCare. Monocytes were enriched by negative selection of CD14-negative cells (Miltenyi) and cultured as above in GM-CSF and IL-4 to generate monocyte-derived DCs, and M-CSF to generate macrophages. Approximately 106 monocytes were cultured in 6-well plates and stimulated on day 3 with irradiated *Mtb* (BEI resources, catalogue: NR-14819) at a multiplicity of infection of 1:1. After 24 hours of stimulation, RNA was extracted using RNeasy (Qiagen), and cDNA was extracted using the quantitect reverse transcription kit (Qiagen) following the manufacturer’s protocol. Reverse transcriptase quantitative real-time PCR as performed using Brilliant III Ultra-Fast SYBR QPCR (Agilent) with the following primers to amplify *FAH* transcripts: 5’ TTCCAGCCACCATAGGAGAC, and 3’ CAGACACCACGACAGAGGAG, normalized to GAPDH amplified using: 5’: GGTCACCAGGGCTGCTTTTA and 3’: TTCCCGTTCTCAGCCTTGAC. Differences in *FAH* expression were calculated using a paired Wilcoxon test.

### Analysis of public RNA-seq data

We obtained gene expression levels from public RNA-seq samples of *Mtb*-infected and non-infected monocyte-derived DCs in 6 individuals54. We used a mixed-effects linear model to test for differential expression between uninfected and or infected DCs, along with heat-inactivated *Mtb* treated DCs, controlling for time as a fixed effect and donor as a random effect.

### Analysis of public methylation data

We retrieved the CpG sites that showed significant differential methylation in dendritic cells in response to either live or heat-killed *Mtb* according to Pacis et al54, following the methylation difference values at four time points post-infection (2, 18, 48, 72h) within a 250 kb window around the transcription start site of *FAH*.

### Data availability

We created a public interactive website for data visualization of (1) gene expression levels in DCs and MPs with differential expression analysis results, (2) cell-type dependent eQTLs, (3) all eQTLs found in each cell-type. Current address: [https://dcmpregulation.shinyapps.io/main/](https://dcmpregulation.shinyapps.io/main/) (it will be hosted in a different domain soon). We will make the gene expression and meta data matrices available in GEO, and raw RNA-seq data available in dbGaP.

## Results

Participants for this study were nested from a larger cohort of index Peruvian TB patients and their household contacts42, from which we previously conducted a genome-wide association study7. We re-recruited a subset of these former household contacts who had been infected with *Mtb* and either progressed to active TB disease (progressors, cases) or remained TB-disease free until sampling (non-progressors, controls) after initial enrollment into the parent cohort6 **(Table 1)**. Due to sample availability, we restricted the analysis to high quality data for 63 progressors and 63 non-progressor controls, sampled a median of 6 years after exposure to TB in the household **(Figure 1A)**. We standardized the protocol to differentiate monocyte-derived dendritic cells (DCs) and macrophages using magnetic enrichment of CD14+ monocytes, and differentiation with a combination of interleukin-4 (IL4) and granulocyte macrophage-colony stimulating factor (GM-CSF) or macrophage-colony stimulating factor (M-CSF), respectively **(Figure 1B)**.

View this table:
[Table 1:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/T1)

Table 1: 
Cohort demographics of individuals included in the RNA-seq and eQTL analysis

We profiled the transcriptome of cells using low-input mRNA sequencing. Macrophages and DCs showed expression of expected genes that mark their cell type specificity **(Supplementary Figure 1)**. For example, genes characteristic of the DC lineage that are absent in macrophages, including group 1 CD1 genes (*CD1a*, *CD1b*, *CD1c* and *CD1e*) 20, 21, were significantly up-regulated in DCs. Similarly, genes characteristic of macrophage differentiation, including V-Set And Immunoglobulin Domain Containing 1 (*VSIG*), *CD14*, Fc Gamma Receptor Ia (*FCGR1A*) and IIIa (*FCGR3A*)21, were significantly upregulated in macrophages compared to DCs.

![Supplementary Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F3.medium.gif)

[Supplementary Figure 1:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F3)

Supplementary Figure 1: 
Volcano plot depicting differential gene expression between samples corresponding to monocyte-derived DCs and M2 macrophages using a multivariate linear regression analysis on log2(transcripts per million +1) adjusted for age, sex and batch (R limma package). Several genes were distinctly cell type-specific with more than 32-fold change and below an FDR-adjusted p-value threshold of 0.05.

First, we asked whether we could detect gene expression differences directly associated with TB disease status of donors. Comparing cells from former progressors and controls, we did not identify any statistically significant differentially expressed genes (FDR < 0.05) in either cell type **(Supplementary Figures 2A and 2B)**. This outcome is in part expected given that RNA was extracted from *in vitro* differentiated monocytes that were collected a median of 6 years after cases had active disease.

![Supplementary Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F4.medium.gif)

[Supplementary Figure 2:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F4)

Supplementary Figure 2: 
Differential gene expression in monocyte-derived DCs (A) or M2 macrophages (B) between former TB cases and controls, using a multivariate linear regression analysis adjusted for age, sex and batch (R limma package). There are no statistically significant TB-associated differentially expressed genes in either cell type.

However, we hypothesized that prior TB progression could be associated with changes in the regulatory landscape of differentiated myeloid cells, which could modulate the effects of genetic variants on gene expression. Therefore, we first identified genes whose expression levels were associated with genetic variants within 1Mb of the gene start (*cis* eQTLs), controlling for age, biological sex, library preparation batch, 15 gene expression and 5 genotyping principal components (PCs). We identified significant eQTLs for 257 genes in macrophages, and 330 genes in DCs (FDR <0.05). The top eQTL in DCs showed that the minor allele C at rs199659767 was associated with a significant change in the expression of Zinc Finger Protein 57 *ZPF57* (p = 1.8 × 10-31) **(Supplementary Table 3)**. The most significant eQTL (p = 3.4 × 10-30) in macrophages showed a positive association between Glutathione S-Transferase Mu 1 (*GSTM1*) expression and the allele A at rs140584594 **(Supplementary Table 4)**. Our results demonstrated a high level of concordance between the mapped eQTLs within our study (macrophages and DCs) and between our macrophage eQTLs and a previously published eQTL study in macrophages39 whose summary statistics were made available by the eQTL Catalogue53. In all comparisons we observed that most of the significant variant-gene pairs (FDR <0.05) had a beta (eQTL effect) in a consistent direction and comparable magnitude (**Supplementary Figure 3**), further supporting the strong agreement between our results and previous findings.

![Supplementary Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F5.medium.gif)

[Supplementary Figure 3:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F5)

Supplementary Figure 3: 
Concordance between effect sizes of significant DC-eQTL events (FDR<0.05) with the macrophage dataset (A), or significant macrophage (MP)-eQTL events (FDR<0.05) with the DC dataset (B). (C) Validation of significant macrophage eQTL using effect sizes of significant eQTL events from monocyte-derived macrophages from Nedelec et al (PMID: 27768889) publicly available on the eQTL catalogue. Pearson correlation analysis results are displayed for A-C.

We then sought to identify genetic regulatory variants whose effects on gene expression were differentially associated with TB disease history. Previously, it has been observed that SNPs associated with the variance of a quantitative trait (vQTL) can increase the detection of gene by environment (GxE) interactions with large effects51. We followed this statistical framework to reduce the number of eQTLs tested for interaction. First, for genes with suggestive eQTLs (p < 0.001), we identified 105 vQTL genes for DCs and 96 for macrophages (FDR < 0.05). For each of these vQTLs, we then tested for TB status interaction **(Methods)**. In the DCs, we observed five vQTLs with significant TB status-eQTL interaction effects (FDR < 0.1) **(Supplementary Table 5A)**. We conducted permutation analysis to rigorously confirm the validity of the eQTL interactions observed in DCs and macrophages by reassigning the TB status 1000 times across samples and retesting for TB status-eQTL interactions **(Supplementary Figure 4A)**. Four out of the five TB status eQTL interactions passed the test (ppermute < 0.05), further suggesting that genetic regulatory variation is associated with history of TB disease. These interactions found in DCs included two protein-coding genes: fumarylacetoacetate hydrolase (*FAH*) and N-ethylmaleimide sensitive factor (*NSF*), and two long non-coding RNA genes: *RP3-340B19.2* and *RP11-667K14.4*. We identified a single strong interaction event for macrophages at this threshold with Lysine Demethylase 6B (*KDM6B*), which did not pass the permutation p-value threshold of 0.05 **(Supplementary Table 5B, Supplementary Figure 4B)**.

![Supplementary Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F6.medium.gif)

[Supplementary Figure 4:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F6)

Supplementary Figure 4: 
Permutation analysis to confirm the validity of the eQTL interactions observed in DCs. We reshuffled TB status (cases and controls) 1000 times across samples and retesting for TB status-eQTL interactions. P-value cutoffs based on this permutation analysis are shown for top TB:genotype interaction model in DCs **(A)** and macrophages **(B)**.

![Supplementary Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F7.medium.gif)

[Supplementary Figure 5:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F7)

Supplementary Figure 5: 
Analysis of *FAH* expression normalized to GAPDH as a housekeeping gene by RT-qPCR in monocyte-derived DC (A) or macrophage (B) samples from healthy donors with or without overnight stimulation with γ-irradiated *Mtb*. P-values are calculated using Wilcoxon rank sum test.

In DCs, we found that an indel on chromosome 15 was the top genetic variant associated with expression levels of *FAH*, a gene encoding Fumarylacetoacetate Hydrolase: the last enzyme in mammalian tyrosine catabolism55. The minor allelic variant of rs142312981 (AAG insertion) was associated with lower expression of *FAH* in former TB progressors but not in non-progressors **(Figure 1C)**. The allele frequency in our cohort is 12%, compared to the estimated Peruvian allele frequency in the 1000 Genomes Project at 7%. Interestingly, this variant’s effect on gene regulation was not observed in monocyte-derived macrophages from the same donors **(Figure 1D)**, suggesting a disease and cell type-specific regulatory axis.

To validate these results, we first reanalyzed public RNA-seq and methylation dataset of *Mtb*-infected and non-infected monocyte-derived DCs in 6 individuals54. Stimulation with live (p = 1.4X10-10) or heat inactivated (p = 5.2X10-10) *Mtb* downregulated *FAH* expression **(Figure 1E)**. In order to experimentally validate these cohort-based transcriptional results, we obtained peripheral monocytes from 6 adult healthy individuals. We differentiated monocytes into dendritic cells or M2 macrophages, and we activated them overnight with γ-irradiated *Mtb*. We measured *FAH* expression levels with RT-qPCR and observed that in all individuals *FAH* was indeed down-regulated after *Mtb* stimulation **(**p = 0.03 for each cell-type, Wilcoxon-test, **Supplementary Figure 5)**. In the same public dataset54, we observed that CpG sites upstream of the *FAH* transcriptional start site and near rs142312981 were differentially methylated after *Mtb* infection **(Figure 1F)**, consistently with an association between *Mtb* infection and epigenetic remodeling of the *FAH* locus.

Several studies have demonstrated that gene expression can be regulated by genetic variants in a cell type-dependent manner27. Hence, we hypothesized that some genetic variants would regulate gene expression at varying degrees in the two cell types in an expression quantitative trait loci analysis. We identified 30 cell-type eQTL interactions at 10% FDR **(Supplementary Table 6)**. For example, rs3907022 was strongly associated with *GSDMA* expression in macrophages, but not DCs **(Figure 2A)**. In contrast, *CD1A* expression was influenced by the rs366316 SNP in DCs but not in macrophages **(Figure 2B)**. We validated this association between rs366316 and CD1A levels in DCs at the protein level with flow cytometry **(**p < 2X10-16, **Figure 2C)**. These two cell type-dependent eQTL events were driven by high expression in one cell type and low expression in the other. However, we also observed other eQTL events in which a gene has similar levels of expression in both cell types, but is regulated by a genetic variant in only one (e.g. *CTSH*, *DFNA5*, **Figure 2D and 2E**). Hence, our analysis captured different types of cell type-dependent genetic regulatory effects in two related myeloid cell types (macrophages and DCs).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/27/2023.06.20.23291558/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/06/27/2023.06.20.23291558/F2)

Figure 2: Cell type-dependent effects of genetic variation on gene expression.
Top cell type-dependent eQTL in DCs **(A)** and macrophages **(B)** were seen when transcript was expressed at higher levels in one cell type compared to the other. **(C)** Validation of CD1a transcript eQTL at the protein level is shown by flow cytometry. P-value is derived from a univariate linear regression model testing the association between rs366316 genotype and CD1a protein expression by flow cytometry. **(D and E)** Figure shows examples of cell type-dependent eQTL events for genes expressed in both DCs **(D)** and macrophages **(E)**.

## Discussion

It remains unclear how host genetic heterogeneity impacts susceptibility to different *Mtb* infection outcomes7, 9, 10, 12, 25, 41, 56. In this study, we re-recruited Peruvians with a history of TB disease or asymptomatic infection with *Mtb*. We analyzed transcriptomes of monocyte-derived DCs and macrophages, which showed the expected cell type-specific transcriptional programs21. Genetic variants in differentiated myeloid cells are known to influence gene expression levels in a cell type and state-dependent manner22, 24, 57. Here, we identified 30 cell type-dependent eQTLs, in an admixed Peruvian cohort, specifically comparing DCs with macrophages, which can serve as a resource for gene regulation in non-European populations. For instance, the regulatory effect of the rs366316 on CD1a, a protein that presents lipid antigens, was previously described in a European-ancestry cohort58, and points to shared gene and protein expression QTL between Europeans and Peruvians. Interestingly, an independent study reported an intronic variant rs411089 in the CD1a locus in a Vietnamese cohort, which was associated with higher TB susceptibility13. However, the DC-specific rs366316-mediated regulation of CD1a expression was not associated with TB history, despite reported roles of CD1a in presenting mycobacterial antigens59, 60.

Noticeably, some of the genetic regulatory events in DCs were associated with the history of TB disease. We discovered a clear interaction between TB status and a non-coding variant in their effects on expression of *FAH*, which encodes the Fumarylacetoacetate Hydrolase enzyme. The genetic regulatory effect was surprisingly restricted to monocyte-derived DC samples in individuals with TB disease history, but not the household contacts who remained healthy. We propose two possible explanations for this genetic regulatory effect associated with TB progression: (1) individuals with susceptibility to TB progression had pre-existing differences in their immune cell state which uncovered their *FAH* eQTL activity, which may have contributed to TB progression, or (2) the *FAH* eQTL appeared in progressors *after* active TB. The TB status interaction eQTL variant was not significantly associated with susceptibility to primary TB progression in the GWAS conducted on this same population (p=0.49)7, and we found evidence that the *FAH* locus changed epigenetically and transcriptionally after *Mtb* infection of DCs. Therefore, we believe the second scenario is more likely. Although durable remodeling of host immunity has classically been associated with acquired immunity in lymphocytes, mouse61 and human62 studies provide precedent for durable regulation of myeloid cells through epigenetic mechanisms, a phenomenon known as trained immunity. Although both progressors and non-progressors have been infected with *Mtb*, a high bacterial burden or strong stimulus occurring during disease can modify the cell state of the host’s immune cells, for example through chromatin architecture, with a lasting effect years after disease and treatment. To our knowledge, pathogen-associated eQTLs have been described using *in vitro* infection models, but not after natural infectious disease in humans 24, 38. However, ‘epigenetic scars’ have been reported as long-term sequalae of TB and thought to be associated with a higher long-term risk of mortality and morbidity63. Therefore, our data raises the question of whether eQTLs activated following an infectious disease could contribute to risk of post-infection sequelae such as recurrent TB, akin to reported hypermethylation in the promoter region of interferon induced protein 44 Like (*IFI44L*) gene post-COVID-1964.

Several loss-of-function mutations in *FAH* are associated with hereditary tyrosinema type I in humans, characterized by elevated blood tyrosine levels65. In contrast to a system-wide missense or antisense loss-of-function mutation, the variant we describe here is a regulatory variant that alters *FAH* expression in a cell type and state-dependent manner. In this metabolic pathway, tyrosine is converted into homogentisate, then into fumarylacetoacetate, which the FAH enzyme converts into fumarate and acetoacetate. Since infection of DCs with live or heat-killed *Mtb* results in down-regulation of *FAH* expression, the data raise the hypothesis that interference with tyrosine metabolism is associated with TB progression. Consistent with this hypothesis, a prior study of African household contacts of index TB patients reported that serum tyrosine and phenylalanine, a precursor to tyrosine, were higher in individuals who progressed to active TB than non-progressor counterparts, consistent with impaired tyrosine metabolism during TB progression5. The role of *FAH* and tyrosine metabolism in TB warrants detailed investigation in future studies.

Population level differences in TB pathogenesis, highlight the importance of including underrepresented groups in genetic studies66, 67. Despite successes in defining variants associated with TB disease susceptibility, many susceptibility loci have been reported to be population-specific7, 9, 10, 12, 56. Differences in genetic ancestry and population history effects have been shown to affect susceptibility to TB, especially in admixed populations41, 68–70. The minor allele frequency (AAG insertion) of rs142312981 allele is 7% in Peruvian (12% in our cohort) compared to 22% for the European reference genomes (NCI, LDproxy tool71). Therefore, this regulatory allele, and others, may have higher impacts in populations where the allele frequencies are higher.

Limitations of this study include profiling the transcriptomes of monocyte-derived DCs and macrophages years after TB disease, where some of the regulatory changes may have reverted to baseline. Furthermore, the *in vitro* differentiation of DCs and macrophages may alter the epigenetic landscape in ways that are different from unmanipulated *ex vivo* cells. Finally, the limited sample size may have reduced the power to detect additional significant eQTL genes and eQTL interactions. Hence, future replication studies with larger sample sizes and additional populations will be needed to validate the findings and discover new post-TB regulatory changes, as currently there are no studies with a similar experimental design with a high enough sample size for QTL analysis with TB interaction.

In conclusion, our study identifies previously unknown genetic regulatory variation associated with history of active TB and with cell type of origin. Especially notable is the DC-specific downregulation of FAH because it is a biologically plausible metabolic gene that produces metabolites, fumarate and acetoacetate, which have been previously proposed to have immune modulatory roles72, 73. The data raise the question of whether these changes could then compromise subsequent immune responses to *Mtb* infection and contribute to the known higher risk of recurrence among individuals with prior TB disease. Subsequent studies should aim at defining mechanistic links to TB susceptibility and mapping these regulatory changes in different cell types and states, to better understand the long-term regulatory sequelae of TB disease.

## Supporting information

Supplementary Tables 1-6 [[supplements/291558_file02.xlsx]](pending:yes)

## Data Availability

We created a public interactive website for data visualization of (1) gene expression levels in DCs and MPs with differential expression analysis results, (2) cell-type dependent eQTLs, (3) all eQTLs found in each cell-type. Current address: https://dcmpregulation.shinyapps.io/main/ (it will be hosted in a different domain soon). We will make the gene expression and meta data matrices available in GEO, and raw RNA-seq data available in dbGaP.

[https://dcmpregulation.shinyapps.io/main/](https://dcmpregulation.shinyapps.io/main/) 

## Acknowledgements

This study was funded by the National Institutes of Health (NIH) TB Research Unit Network (Grant U19 AI111224-01, PIs: Moody and Murray), the National Institute of Allergy and Infectious Diseases (R01 AI049313, PI: Moody, and R01 AI175614, PI: Suliman). SS is also supported by an investigator grant from the Chan Zuckerberg Biohub. None of the authors declare any competing interests.

*   Received June 20, 2023.
*   Revision received June 20, 2023.
*   Accepted June 27, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  WHO. Global Tuberculosis Report 2022. Accessed online: [https://www.who.int/publications/i/item/9789240061729](https://www.who.int/publications/i/item/9789240061729) (2022).
    
    

2.  Vega, V., Rodriguez, S., Van der Stuyft, P., Seas, C. & Otero, L. Recurrent TB: a systematic review and meta-analysis of the incidence rates and the proportions of relapses and reinfections. Thorax 76, 494–502 (2021). [https://doi.org:10.1136/thoraxjnl-2020-215449](https://doi.org:10.1136/thoraxjnl-2020-215449)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjg6Ijc2LzUvNDk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  Berry, M. P. et al. An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature 466, 973–977 (2010). [https://doi.org:10.1038/nature09247](https://doi.org:10.1038/nature09247)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09247&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20725040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281030300034&link_type=ISI) 

4.  Zak, D. E. et al. A blood RNA signature for tuberculosis disease risk: a prospective cohort study. Lancet 387, 2312–2322 (2016). [https://doi.org:10.1016/S0140-6736(15)01316-1](https://doi.org:10.1016/S0140-6736(15)01316-1)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(15)01316-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27017310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

5.  Weiner, J., 3rd. et al. Metabolite changes in blood predict the onset of tuberculosis. Nat Commun 9, 5208 (2018). [https://doi.org:10.1038/s41467-018-07635-7](https://doi.org:10.1038/s41467-018-07635-7)
    
    

6.  Nathan, A. et al. Multimodally profiling memory T cells from a tuberculosis cohort identifies cell state associations with demographics, environment and disease. Nat Immunol 22, 781–793 (2021). [https://doi.org:10.1038/s41590-021-00933-1](https://doi.org:10.1038/s41590-021-00933-1)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41590-021-00933-1&link_type=DOI) 

7.  Luo, Y. et al. Early progression to active tuberculosis is a highly heritable trait driven by 3q23 in Peruvians. Nat Commun 10, 3765 (2019). [https://doi.org:10.1038/s41467-019-11664-1](https://doi.org:10.1038/s41467-019-11664-1)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-11664-1&link_type=DOI) 

8.  Nathan, A. et al. Single-cell eQTL models reveal dynamic T cell state dependence of disease loci. Nature 606, 120–128 (2022). [https://doi.org:10.1038/s41586-022-04713-1](https://doi.org:10.1038/s41586-022-04713-1)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-04713-1&link_type=DOI) 

9.  Curtis, J. et al. Susceptibility to tuberculosis is associated with variants in the ASAP1 gene encoding a regulator of dendritic cell migration. Nat Genet 47, 523–527 (2015). [https://doi.org:10.1038/ng.3248](https://doi.org:10.1038/ng.3248)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25774636&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

10. Thye, T. et al. Genome-wide association analyses identifies a susceptibility locus for tuberculosis on chromosome 18q11.2. Nat Genet 42, 739–741 (2010). [https://doi.org:10.1038/ng.639](https://doi.org:10.1038/ng.639)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.639&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20694014&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281388400005&link_type=ISI) 

11. van Tong, H., Velavan, T. P., Thye, T. & Meyer, C. G. Human genetic factors in tuberculosis: an update. Trop Med Int Health 22, 1063–1071 (2017). [https://doi.org:10.1111/tmi.12923](https://doi.org:10.1111/tmi.12923)
    
    

12. Schurz, H. et al. A Sex-Stratified Genome-Wide Association Study of Tuberculosis Using a Multi-Ethnic Genotyping Array. Front Genet 9, 678 (2018). [https://doi.org:10.3389/fgene.2018.00678](https://doi.org:10.3389/fgene.2018.00678)
    
    

13. Seshadri, C. et al. A polymorphism in human CD1A is associated with susceptibility to tuberculosis. Genes Immun 15, 195–198 (2014). [https://doi.org:10.1038/gene.2014.5](https://doi.org:10.1038/gene.2014.5)
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24500401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

14. Khader, S. A. et al. Targeting innate immunity for tuberculosis vaccination. J Clin Invest 129, 3482–3491 (2019). [https://doi.org:10.1172/JCI128877](https://doi.org:10.1172/JCI128877)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI128877&link_type=DOI) 

15. Scriba, T. J. et al. Sequential inflammatory processes define human progression from M. tuberculosis infection to tuberculosis disease. PLoS Pathog 13, e1006687 (2017). [https://doi.org:10.1371/journal.ppat.1006687](https://doi.org:10.1371/journal.ppat.1006687)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1006687&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29145483&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

16. Suliman, S. et al. Four-gene Pan-African Blood Signature Predicts Progression to Tuberculosis. Am J Respir Crit Care Med (2018). [https://doi.org:10.1164/rccm.201711-2340OC](https://doi.org:10.1164/rccm.201711-2340OC)
    
    

17. Singhania, A. et al. A modular transcriptional signature identifies phenotypic heterogeneity of human tuberculosis infection. Nat Commun 9, 2308 (2018). [https://doi.org:10.1038/s41467-018-04579-w](https://doi.org:10.1038/s41467-018-04579-w)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-04579-w&link_type=DOI) 

18. VanderVen, B. C., Huang, L., Rohde, K. H. & Russell, D. G. The Minimal Unit of Infection: Mycobacterium tuberculosis in the Macrophage. Microbiol Spectr 4 (2016). [https://doi.org:10.1128/microbiolspec.TBTB2-0025-2016](https://doi.org:10.1128/microbiolspec.TBTB2-0025-2016)
    
    

19. Rodrigues, T. S., Conti, B. J., Fraga-Silva, T. F. C., Almeida, F. & Bonato, V. L. D. Interplay between alveolar epithelial and dendritic cells and Mycobacterium tuberculosis. J Leukoc Biol 108, 1139–1156 (2020). [https://doi.org:10.1002/JLB.4MR0520-112R](https://doi.org:10.1002/JLB.4MR0520-112R)
    
    

20. Sallusto, F. & Lanzavecchia, A. Efficient presentation of soluble antigen by cultured human dendritic cells is maintained by granulocyte/macrophage colony-stimulating factor plus interleukin 4 and downregulated by tumor necrosis factor alpha. J Exp Med 179, 1109–1118 (1994). [https://doi.org:10.1084/jem.179.4.1109](https://doi.org:10.1084/jem.179.4.1109)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjEwOiIxNzkvNC8xMTA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

21. Sander, J. et al. Cellular Differentiation of Human Monocytes Is Regulated by Time-Dependent Interleukin-4 Signaling and the Transcriptional Regulator NCOR2. Immunity 47, 1051–1066 e1012 (2017). [https://doi.org:10.1016/j.immuni.2017.11.024](https://doi.org:10.1016/j.immuni.2017.11.024)
    
    

22. Blischak, J. D. et al. Predicting susceptibility to tuberculosis based on gene expression profiling in dendritic cells. Sci Rep 7, 5702 (2017). [https://doi.org:10.1038/s41598-017-05878-w](https://doi.org:10.1038/s41598-017-05878-w)
    
    

23. Roura-Mir, C. et al. Mycobacterium tuberculosis regulates CD1 antigen presentation pathways through TLR-2. J Immunol 175, 1758–1766 (2005). [https://doi.org:10.4049/jimmunol.175.3.1758](https://doi.org:10.4049/jimmunol.175.3.1758)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE3NS8zLzE3NTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8yNy8yMDIzLjA2LjIwLjIzMjkxNTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

24. Barreiro, L. B. et al. Deciphering the genetic architecture of variation in the immune response to Mycobacterium tuberculosis infection. Proc Natl Acad Sci U S A 109, 1204–1209 (2012). [https://doi.org:10.1073/pnas.1115761109](https://doi.org:10.1073/pnas.1115761109)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTA5LzQvMTIwNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzI3LzIwMjMuMDYuMjAuMjMyOTE1NTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

25. Jiao, L. et al. A Novel Genetic Variation in NCF2, the Core Component of NADPH Oxidase, Contributes to the Susceptibility of Tuberculosis in Western Chinese Han Population. DNA Cell Biol 39, 57–62 (2020). [https://doi.org:10.1089/dna.2019.5082](https://doi.org:10.1089/dna.2019.5082)
    
    

26. Consortium, G. T. et al. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017). [https://doi.org:10.1038/nature24277](https://doi.org:10.1038/nature24277)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29022597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000412829500039&link_type=ISI) 

27. Consortium, G. T. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 369, 1318–1330 (2020). [https://doi.org:10.1126/science.aaz1776](https://doi.org:10.1126/science.aaz1776)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

28. Gutierrez-Arcelus, M. et al. Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing. PLoS Genet 11, e1004958 (2015). [https://doi.org:10.1371/journal.pgen.1004958](https://doi.org:10.1371/journal.pgen.1004958)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004958&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25634236&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

29. Soskic, B. et al. Immune disease risk variants regulate gene expression dynamics during CD4(+) T cell activation. Nat Genet 54, 817–826 (2022). [https://doi.org:10.1038/s41588-022-01066-3](https://doi.org:10.1038/s41588-022-01066-3)
    
    

30. Strober, B. J. et al. Dynamic genetic regulation of gene expression during cellular differentiation. Science 364, 1287–1290 (2019). [https://doi.org:10.1126/science.aaw0040](https://doi.org:10.1126/science.aaw0040)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjQvNjQ0Ny8xMjg3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

31. Randolph, H. E. et al. Genetic ancestry effects on the response to viral infection are pervasive but cell type specific. Science 374, 1127–1133 (2021). [https://doi.org:10.1126/science.abg0928](https://doi.org:10.1126/science.abg0928)
    
    

32. Yazar, S. et al. Single-cell eQTL mapping identifies cell type-specific genetic control of autoimmune disease. Science 376, eabf3041 (2022). [https://doi.org:10.1126/science.abf3041](https://doi.org:10.1126/science.abf3041)
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

33. Fairfax, B. P. et al. Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science 343, 1246949 (2014). [https://doi.org:10.1126/science.1246949](https://doi.org:10.1126/science.1246949)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE2OiIzNDMvNjE3NS8xMjQ2OTQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

34. Gutierrez-Arcelus, M. et al. Allele-specific expression changes dynamically during T cell activation in HLA and other autoimmune loci. Nat Genet 52, 247–253 (2020). [https://doi.org:10.1038/s41588-020-0579-4](https://doi.org:10.1038/s41588-020-0579-4)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-0579-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

35. Umans, B. D., Battle, A. & Gilad, Y. Where Are the Disease-Associated eQTLs? Trends Genet 37, 109–124 (2021). [https://doi.org:10.1016/j.tig.2020.08.009](https://doi.org:10.1016/j.tig.2020.08.009)
    
    

36. Fairfax, B. P. & Knight, J. C. Genetics of gene expression in immunity to infection. Curr Opin Immunol 30, 63–71 (2014). [https://doi.org:10.1016/j.coi.2014.07.001](https://doi.org:10.1016/j.coi.2014.07.001)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.coi.2014.07.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25078545&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

37. Klunk, J. et al. Evolution of immune genes is associated with the Black Death. Nature 611, 312–319 (2022). [https://doi.org:10.1038/s41586-022-05349-x](https://doi.org:10.1038/s41586-022-05349-x)
    
    

38. Oelen, R. et al. Single-cell RNA-sequencing of peripheral blood mononuclear cells reveals widespread, context-specific gene expression regulation upon pathogenic exposure. Nat Commun 13, 3267 (2022). [https://doi.org:10.1038/s41467-022-30893-5](https://doi.org:10.1038/s41467-022-30893-5)
    
    

39. Nedelec, Y. et al. Genetic Ancestry and Natural Selection Drive Population Differences in Immune Responses to Pathogens. Cell 167, 657–669 e621 (2016). [https://doi.org:10.1016/j.cell.2016.09.025](https://doi.org:10.1016/j.cell.2016.09.025)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2016.09.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27768889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

40. Comas, I. et al. Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nat Genet 45, 1176–1182 (2013). [https://doi.org:10.1038/ng.2744](https://doi.org:10.1038/ng.2744)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2744&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23995134&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

41. Asgari, S. et al. Higher native Peruvian genetic ancestry proportion is associated with tuberculosis progression risk. Cell Genom 2 (2022). [https://doi.org:10.1016/j.xgen.2022.100151](https://doi.org:10.1016/j.xgen.2022.100151)
    
    

42. Becerra, M. C. et al. Transmissibility and potential for disease progression of drug resistant Mycobacterium tuberculosis: prospective cohort study. BMJ 367, l5894 (2019). [https://doi.org:10.1136/bmj.l5894](https://doi.org:10.1136/bmj.l5894)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNjcvb2N0MjRfMy9sNTg5NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzI3LzIwMjMuMDYuMjAuMjMyOTE1NTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

43. Picelli, S. et al. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc 9, 171–181 (2014). [https://doi.org:10.1038/nprot.2014.006](https://doi.org:10.1038/nprot.2014.006)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nprot.2014.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24385147&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

44. Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013). [https://doi.org:10.1093/bioinformatics/bts635](https://doi.org:10.1093/bioinformatics/bts635)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts635&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23104886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312654600003&link_type=ISI) 

45. Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011). [https://doi.org:10.1186/1471-2105-12-323](https://doi.org:10.1186/1471-2105-12-323)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-12-323&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21816040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

46. DeLuca, D. S. et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 28, 1530–1532 (2012). [https://doi.org:10.1093/bioinformatics/bts196](https://doi.org:10.1093/bioinformatics/bts196)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts196&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22539670&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304537000017&link_type=ISI) 

47. Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43, e47 (2015). [https://doi.org:10.1093/nar/gkv007](https://doi.org:10.1093/nar/gkv007)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25605792&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

48. O’Connell, J. et al. A general approach for haplotype phasing across the full spectrum of relatedness. PLoS Genet 10, e1004234 (2014). [https://doi.org:10.1371/journal.pgen.1004234](https://doi.org:10.1371/journal.pgen.1004234)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004234&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24743097&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

49. Howie, B., Fuchsberger, C., Stephens, M., Marchini, J. & Abecasis, G. R. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet 44, 955–959 (2012). [https://doi.org:10.1038/ng.2354](https://doi.org:10.1038/ng.2354)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2354&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22820512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

50. Delaneau, O. et al. A complete tool set for molecular QTL discovery and analysis. Nat Commun 8, 15452 (2017). [https://doi.org:10.1038/ncomms15452](https://doi.org:10.1038/ncomms15452)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms15452&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28516912&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

51. Marderstein, A. R. et al. Leveraging phenotypic variability to identify genetic interactions in human phenotypes. Am J Hum Genet 108, 49–67 (2021). [https://doi.org:10.1016/j.ajhg.2020.11.016](https://doi.org:10.1016/j.ajhg.2020.11.016)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2020.11.016&link_type=DOI) 

52. Storey, J. D. a. B., A. J. qvalue: Q-value estimation for false discovery rate control. R package version 2.28.0 [http://github.com/jdstorey/qvalue](http://github.com/jdstorey/qvalue) (2022).
    
    

53. Kerimov, N. et al. A compendium of uniformly processed human gene expression and splicing quantitative trait loci. Nat Genet 53, 1290–1299 (2021). [https://doi.org:10.1038/s41588-021-00924-w](https://doi.org:10.1038/s41588-021-00924-w)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00924-w&link_type=DOI) 

54. Pacis, A. et al. Gene activation precedes DNA demethylation in response to infection in human dendritic cells. Proc Natl Acad Sci U S A 116, 6938–6943 (2019). [https://doi.org:10.1073/pnas.1814700116](https://doi.org:10.1073/pnas.1814700116)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE2LzE0LzY5MzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8yNy8yMDIzLjA2LjIwLjIzMjkxNTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. Weiss, A. K. H., Loeffler, J. R., Liedl, K. R., Gstach, H. & Jansen-Durr, P. The fumarylacetoacetate hydrolase (FAH) superfamily of enzymes: multifunctional enzymes from microbes to mitochondria. Biochem Soc Trans 46, 295–309 (2018). [https://doi.org:10.1042/BST20170518](https://doi.org:10.1042/BST20170518)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoicHBiaW9zdCI7czo1OiJyZXNpZCI7czo4OiI0Ni8yLzI5NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzI3LzIwMjMuMDYuMjAuMjMyOTE1NTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

56. Thye, T. et al. Common variants at 11p13 are associated with susceptibility to tuberculosis. Nat Genet 44, 257–259 (2012). [https://doi.org:10.1038/ng.1080](https://doi.org:10.1038/ng.1080)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.1080&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22306650&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

57. Blischak, J. D., Tailleux, L., Mitrano, A., Barreiro, L. B. & Gilad, Y. Mycobacterial infection induces a specific human innate immune response. Sci Rep 5, 16882 (2015). [https://doi.org:10.1038/srep16882](https://doi.org:10.1038/srep16882)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep16882&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26586179&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

58. Seshadri, C. et al. Human CD1a deficiency is common and genetically regulated. J Immunol 191, 1586–1593 (2013). [https://doi.org:10.4049/jimmunol.1300575](https://doi.org:10.4049/jimmunol.1300575)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE5MS80LzE1ODYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8yNy8yMDIzLjA2LjIwLjIzMjkxNTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

59. Kasmar, A. G. et al. Cutting Edge: CD1a tetramers and dextramers identify human lipopeptide-specific T cells ex vivo. J Immunol 191, 4499–4503 (2013). [https://doi.org:10.4049/jimmunol.1301660](https://doi.org:10.4049/jimmunol.1301660)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE5MS85LzQ0OTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8yNy8yMDIzLjA2LjIwLjIzMjkxNTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

60. Moody, D. B. et al. T cell activation by lipopeptide antigens. Science 303, 527–531 (2004). [https://doi.org:10.1126/science.1089353](https://doi.org:10.1126/science.1089353)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMDMvNTY1Ny81MjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8yNy8yMDIzLjA2LjIwLjIzMjkxNTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

61. Netea, M. G., Quintin, J. & van der Meer, J. W. Trained immunity: a memory for innate host defense. Cell Host Microbe 9, 355–361 (2011). [https://doi.org:10.1016/j.chom.2011.04.006](https://doi.org:10.1016/j.chom.2011.04.006)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2011.04.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21575907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291289500005&link_type=ISI) 

62. Kleinnijenhuis, J. et al. Bacille Calmette-Guerin induces NOD2-dependent nonspecific protection from reinfection via epigenetic reprogramming of monocytes. Proc Natl Acad Sci U S A 109, 17537–17542 (2012). [https://doi.org:10.1073/pnas.1202870109](https://doi.org:10.1073/pnas.1202870109)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzQzLzE3NTM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjcvMjAyMy4wNi4yMC4yMzI5MTU1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. DiNardo, A. R., Netea, M. G. & Musher, D. M. Postinfectious Epigenetic Immune Modifications - A Double-Edged Sword. N Engl J Med 384, 261–270 (2021). [https://doi.org:10.1056/NEJMra2028358](https://doi.org:10.1056/NEJMra2028358)
    
    

64. Lee, Y. et al. EWAS of post-COVID-19 patients shows methylation differences in the immune-response associated gene, IFI44L, three months after COVID-19 infection. Sci Rep 12, 11478 (2022). [https://doi.org:10.1038/s41598-022-15467-1](https://doi.org:10.1038/s41598-022-15467-1)
    
    

65. Morrow, G., Angileri, F. & Tanguay, R. M. Molecular Aspects of the FAH Mutations Involved in HT1 Disease. Adv Exp Med Biol 959, 25–48 (2017). [https://doi.org:10.1007/978-3-319-55780-9_3](https://doi.org:10.1007/978-3-319-55780-9_3)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-3-319-55780-9_3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28755182&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

66. Bentley, A. R., Callier, S. & Rotimi, C. N. Diversity and inclusion in genomic research: why the uneven progress? J Community Genet 8, 255–266 (2017). [https://doi.org:10.1007/s12687-017-0316-6](https://doi.org:10.1007/s12687-017-0316-6)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12687-017-0316-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 

67. Swart, Y. et al. cis-eQTL mapping of TB-T2D comorbidity elucidates the involvement of African ancestry in TB susceptibility. MedRxiv pre-print: [https://www.biorxiv.org/content/10.1101/2022.10.19.512814v1](https://www.biorxiv.org/content/10.1101/2022.10.19.512814v1) Accessed online on December 23rd, 2022 (2022).
    
    

68. Kerner, G. et al. Human ancient DNA analyses reveal the high burden of tuberculosis in Europeans over the last 2,000 years. Am J Hum Genet 108, 517–524 (2021). [https://doi.org:10.1016/j.ajhg.2021.02.009](https://doi.org:10.1016/j.ajhg.2021.02.009)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2021.02.009&link_type=DOI) 

69. Chimusa, E. R. et al. Genome-wide association study of ancestry-specific TB risk in the South African Coloured population. Hum Mol Genet 23, 796–809 (2014). [https://doi.org:10.1093/hmg/ddt462](https://doi.org:10.1093/hmg/ddt462)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddt462&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24057671&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F27%2F2023.06.20.23291558.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000330841700019&link_type=ISI) 

70. Daya, M., van der Merwe, L., van Helden, P. D., Moller, M. & Hoal, E. G. The role of ancestry in TB susceptibility of an admixed South African population. Tuberculosis (Edinb*)* 94, 413–420 (2014). [https://doi.org:10.1016/j.tube.2014.03.012](https://doi.org:10.1016/j.tube.2014.03.012)
    
    

71. Alexander, T. A. & Machiela, M. J. LDpop: an interactive online tool to calculate and visualize geographic LD patterns. BMC Bioinformatics 21, 14 (2020). [https://doi.org:10.1186/s12859-020-3340-1](https://doi.org:10.1186/s12859-020-3340-1)
    
    

72. Adam, C. et al. Acetoacetate protects macrophages from lactic acidosis-induced mitochondrial dysfunction by metabolic reprograming. Nat Commun 12, 7115 (2021). [https://doi.org:10.1038/s41467-021-27426-x](https://doi.org:10.1038/s41467-021-27426-x)
    
    

73. Mita-Mendoza, N. K. et al. Dimethyl fumarate reduces TNF and Plasmodium falciparum induced brain endothelium activation in vitro. Malar J 19, 376 (2020). [https://doi.org:10.1186/s12936-020-03447-7](https://doi.org:10.1186/s12936-020-03447-7)

 [1]: /embed/graphic-2.gif