Abstract
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.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
The research was supported by the Bill and Melinda Gates Foundation (grant OPP1151836 to T.R.H., W.H.B., C.M.S., and H.M.K.), the National Institutes of Health (grant R01AI124348 to W.H.B., T.R.H., C.M.S., and H.M.K.; grant U01AI115642 to W.H.B., T.R.H., C.M.S., and H.M.K.; grant K24AI137310 to T.R.H.; R33AI138272 (to TRH, WHB, HMK, CMS), NIH grants K08AI143926 and T32AI007044 (to JDS), T32NR016913 (to HH), grant U19AI162583 (to HMK, WHB, CMS, and TRH), contract no. 75N93019C00071 to T.R.H., W.H.B., C.M.S., and H.M.K.; and contract no. NO1AI70022 to W.H.B., T.R.H., C.M.S., and H.M.K.). The funders had no role in the experimental design or analysis.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The study protocols were approved by the National AIDS Research Committee, the Uganda National Council on Science and Technology, and the Institutional Review Boards at the University Hospitals Cleveland Medical Center and the University of Washington.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Footnotes
Conflict of interest: The authors have declared that no conflict of interest exists.
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/hawn-lab/RSTR_eQTL_public.
