Fine-mapping analysis of a chromosome 2 region linked to resistance to Mycobacterium tuberculosis infection in Uganda reveals potential regulatory variants

Genetic epidemiology of resistance to M. tuberculosis Infection: importance of study design and recent findings

Resistance to M. tuberculosis, often referred to as "RSTR" in the literature, is being increasingly studied because of its potential relevance as a clinical outcome in vaccine studies. This review starts by addressing the importance of epidemiological characterization of this phenotype, and ongoing challenges in that characterization. Then, this review summarizes the extant genetic and genomic studies of this phenotype, including heritability studies, candidate gene studies, and genome-wide association studies, as well as whole transcriptome studies. Findings from recent studies that used longitudinal characterization of the RSTR phenotype are compared to those using a cross-sectional definition, and the challenges of using tuberculin skin test and interferon-gamma release assay are discussed. Finally, future directions are proposed. Since this is a rapidly evolving area of public health significance, this review will help frame future research questions and study designs.

Differentially expressed transcript isoforms associate with resistance to tuberculin skin test and interferon gamma release assay conversion

BACKGROUND

A mechanistic understanding of uncommon immune outcomes such as resistance to infection has led to the development of novel therapies. Using gene level analytic methods, we previously found distinct monocyte transcriptional responses associated with resistance to Mycobacterium tuberculosis (Mtb) infection defined as persistently negative tuberculin skin test (TST) and interferon gamma release assay (IGRA) reactivity among highly exposed contacts (RSTR phenotype).

OBJECTIVE

Using transcript isoform analyses, we aimed to identify novel RSTR-associated genes hypothesizing that previous gene-level differential expression analysis obscures isoform-specific differences that contribute to phenotype.

MATERIALS AND METHODS

Monocytes from 49 RSTR versus 52 subjects with latent Mtb infection (LTBI) were infected with M. tuberculosis (H37Rv) or left unstimulated (media) prior to RNA isolation and sequencing. RSTR-associated gene expression was then identified using differential transcript isoform analysis.

RESULTS

We identified 81 differentially expressed transcripts (DETs) in 70 genes (FDR <0.05) comparing RSTR and LTBI phenotypes with the majority (n = 79 DETs) identified under Mtb-stimulated conditions. Seventeen of these genes were previously identified with gene-level bulk RNAseq analyses including genes in the IFNγ response that had increased expression among LTBI subjects, findings consistent with a clinical phenotype based on IGRA reactivity. Among the subset of 23 genes with positive differential expression among Mtb-infected RSTR monocytes, 13 were not previously identified. These novel DET genes included PDE4A and ZEB2, which each had multiple DETs with higher expression among RSTR subjects, and ACSL4 and GAPDH that each had a single transcript isoform associated with RSTR.

CONCLUSION AND LIMITATIONS

Transcript isoform-specific analyses identify transcriptional associations, such as those associated with resistance to TST/IGRA conversion, that are obscured when using gene-level approaches. These findings should be validated with additional RSTR cohorts and whether the newly identified candidate resistance genes directly influence the monocyte Mtb response requires functional study.

Monocyte Transcriptional Responses to Mycobacterium tuberculosis Associate with Resistance to Tuberculin Skin Test and Interferon Gamma Release Assay Conversion

Heavy exposure to Mycobacterium tuberculosis, the etiologic agent of tuberculosis (TB) and among the top infectious killers worldwide, results in infection that is cleared, contained, or progresses to disease. Some heavily exposed tuberculosis contacts show no evidence of infection using the tuberculin skin test (TST) and interferon gamma release assay (IGRA); yet the mechanisms underlying this "resister" (RSTR) phenotype are unclear. To identify transcriptional responses that distinguish RSTR monocytes, we performed transcriptome sequencing (RNA-seq) on monocytes isolated from heavily exposed household contacts in Uganda and gold miners in South Africa after ex vivo M. tuberculosis infection. Gene set enrichment analysis (GSEA) revealed several gene pathways that were consistently enriched in response to M. tuberculosis among RSTR subjects compared to controls with positive TST/IGRA testing (latent TB infection [LTBI]) across Uganda and South Africa. The most significantly enriched gene set in which expression was increased in RSTR relative to LTBI M. tuberculosis-infected monocytes was the tumor necrosis factor alpha (TNF-α) signaling pathway whose core enrichment (leading edge) substantially overlapped across RSTR populations. These leading-edge genes included candidate resistance genes (ABCA1 and DUSP2) with significantly increased expression among Uganda RSTRs (false-discovery rate [FDR], <0.1). The distinct monocyte transcriptional response to M. tuberculosis among RSTR subjects, including increased expression of the TNF signaling pathway, highlights genes and inflammatory pathways that may mediate resistance to TST/IGRA conversion and provides therapeutic targets to enhance host restriction of M. tuberculosis intracellular infection. IMPORTANCE After heavy M. tuberculosis exposure, the events that determine why some individuals resist TST/IGRA conversion are poorly defined. Enrichment of the TNF signaling gene set among RSTR monocytes from multiple distinct cohorts suggests an important role for the monocyte TNF response in determining this alternative immune outcome. These TNF responses to M. tuberculosis among RSTRs may contribute to antimicrobial programs that result in early clearance or the priming of alternative (gamma interferon-independent) cellular responses.

Resistance to TST/IGRA conversion in Uganda: Heritability and Genome-Wide Association Study

Background

Pulmonary tuberculosis (TB) is one of the most deadly pathogens on earth. However, the majority of people have resistance to active disease. Further, some individuals, termed resisters (RSTRs), do not develop traditional latent tuberculosis (LTBI). The RSTR phenotype is important for understanding pathogenesis and preventing TB. The host genetic underpinnings of RSTR are largely understudied.

Methods

In a cohort of 908 Ugandan subjects with genome-wide data on single nucleotide polymorphisms, we assessed the heritability of the RSTR phenotype and other TB phenotypes using restricted maximum likelihood estimation (REML). We then used a subset of 263 RSTR and LTBI subjects with high quality phenotyping and long-term follow-up to identify DNA variants genome-wide associated with the RSTR phenotype relative to LTBI subjects in a case-control GWAS design and annotated and enriched these variants to better understand their role in TB pathogenesis.

Results

The heritability of the TB outcomes was very high, at 55% for TB vs. LTBI and 50.4% for RSTR vs. LTBI among HIV- subjects, controlling for age and sex. We identified 27 loci associated with the RSTR phenotype (P<5e-05) and our annotation and enrichment analyses suggest an important regulatory role for many of them.

Interpretation

The heritability results show that the genetic contribution to variation in TB outcomes is very high and our GWAS results highlight variants that may play an important role in resistance to infection as well as TB pathogenesis as a whole.

Interaction between M. tuberculosis Lineage and Human Genetic Variants Reveals Novel Pathway Associations with Severity of TB

Tuberculosis (TB) remains a major public health threat globally, especially in sub-Saharan Africa. Both human and Mycobacterium tuberculosis (MTBC) genetic variation affect TB outcomes, but few studies have examined if and how the two genomes interact to affect disease. We hypothesize that long-term coexistence between human genomes and MTBC lineages modulates disease to affect its severity. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which we identified three MTBC lineages, of which one, L4.6-Uganda, is clearly derived and hence recent. We quantified TB severity using the Bandim TBscore and examined the interaction between MTBC lineage and human single-nucleotide polymorphisms (SNPs) genome-wide, in two independent cohorts of TB cases (n = 149 and n = 127). We found a significant interaction between an SNP in PPIAP2 and the Uganda lineage (combined p = 4 × 10⁻⁸). PPIAP2 is a pseudogene that is highly expressed in immune cells. Pathway and eQTL analyses indicated potential roles between coevolving SNPs and cellular replication and metabolism as well as platelet aggregation and coagulation. This finding provides further evidence that host–pathogen interactions affect clinical presentation differently than host and pathogen genetic variation independently, and that human–MTBC coevolution is likely to explain patterns of disease severity.

Monocyte metabolic transcriptional programs associate with resistance to tuberculin skin test/interferon-γ release assay conversion

After extensive exposure to Mycobacterium tuberculosis (Mtb), most individuals acquire latent Mtb infection (LTBI) defined by a positive tuberculin skin test (TST) or interferon-γ release assay (IGRA). To identify mechanisms of resistance to Mtb infection, we compared transcriptional profiles from highly exposed contacts who resist TST/IGRA conversion (resisters, RSTRs) and controls with LTBI using RNAseq. Gene sets related to carbon metabolism and free fatty acid (FFA) transcriptional responses enriched across 2 independent cohorts suggesting RSTR and LTBI monocytes have distinct activation states. We compared intracellular Mtb replication in macrophages treated with FFAs and found that palmitic acid (PA), but not oleic acid (OA), enhanced Mtb intracellular growth. This PA activity correlated with its inhibition of proinflammatory cytokines in Mtb-infected cells. Mtb growth restriction in PA-treated macrophages was restored by activation of AMP kinase (AMPK), a central host metabolic regulator known to be inhibited by PA. Finally, we genotyped AMPK variants and found 7 SNPs in PRKAG2, which encodes the AMPK-γ subunit, that strongly associated with RSTR status. Taken together, RSTR and LTBI phenotypes are distinguished by FFA transcriptional programs and by genetic variation in a central metabolic regulator, which suggests immunometabolic pathways regulate TST/IGRA conversion.

BCG-induced protection against Mycobacterium tuberculosis infection: Evidence, mechanisms, and implications for next-generation vaccines

The tuberculosis (TB) vaccine Bacillus Calmette-Guérin (BCG) was introduced 100 years ago, but as it provides insufficient protection against TB disease, especially in adults, new vaccines are being developed and evaluated. The discovery that BCG protects humans from becoming infected with Mycobacterium tuberculosis (Mtb) and not just from progressing to TB disease provides justification for considering Mtb infection as an endpoint in vaccine trials. Such trials would require fewer participants than those with disease as an endpoint. In this review, we first define Mtb infection and disease phenotypes that can be used for mechanistic studies and/or endpoints for vaccine trials. Secondly, we review the evidence for BCG-induced protection against Mtb infection from observational and BCG re-vaccination studies, and discuss limitations and variation of this protection. Thirdly, we review possible underlying mechanisms for BCG efficacy against Mtb infection, including alternative T cell responses, antibody-mediated protection, and innate immune mechanisms, with a specific focus on BCG-induced trained immunity, which involves epigenetic and metabolic reprogramming of innate immune cells. Finally, we discuss the implications for further studies of BCG efficacy against Mtb infection, including for mechanistic research, and their relevance to the design and evaluation of new TB vaccines.

Genome-wide association study of resistance to Mycobacterium tuberculosis infection identifies a locus at 10q26.2 in three distinct populations

The natural history of tuberculosis (TB) is characterized by a large inter-individual outcome variability after exposure to Mycobacterium tuberculosis. Specifically, some highly exposed individuals remain resistant to M. tuberculosis infection, as inferred by tuberculin skin test (TST) or interferon-gamma release assays (IGRAs). We performed a genome-wide association study of resistance to M. tuberculosis infection in an endemic region of Southern Vietnam. We enrolled household contacts (HHC) of pulmonary TB cases and compared subjects who were negative for both TST and IGRA (n = 185) with infected individuals (n = 353) who were either positive for both TST and IGRA or had a diagnosis of TB. We found a genome-wide significant locus on chromosome 10q26.2 with a cluster of variants associated with strong protection against M. tuberculosis infection (OR = 0.42, 95%CI 0.35–0.49, P = 3.71×10⁻⁸, for the genotyped variant rs17155120). The locus was replicated in a French multi-ethnic HHC cohort and a familial admixed cohort from a hyper-endemic area of South Africa, with an overall OR for rs17155120 estimated at 0.50 (95%CI 0.45–0.55, P = 1.26×10⁻⁹). The variants are located in intronic regions and upstream of C10orf90, a tumor suppressor gene which encodes an ubiquitin ligase activating the transcription factor p53. In silico analysis showed that the protective alleles were associated with a decreased expression in monocytes of the nearby gene ADAM12 which could lead to an enhanced response of Th17 lymphocytes. Our results reveal a novel locus controlling resistance to M. tuberculosis infection across different populations.

There is strong epidemiological evidence that a proportion of highly exposed individuals remain resistant to M. tuberculosis infection, as shown by a negative result for Tuberculin Skin Test (TST) or IFN-γ Release Assays (IGRAs). We performed a genome-wide association study between resistant and infected individuals, which were carefully selected employing a household contact design to maximize exposure by infectious index patients. We employed stringently defined concordant results for both TST and IGRA assays to avoid misclassifications. We discovered a locus at 10q26.2 associated with resistance to M. tuberculosis infection in a Vietnamese discovery cohort. This locus could be replicated in two independent cohorts from different epidemiological settings and of diverse ancestries enrolled in France and South Africa.

Phenotype Definition for "Resisters" to Mycobacterium tuberculosis Infection in the Literature-A Review and Recommendations

Tuberculosis (TB) remains a worldwide problem. Despite the high disease rate, not all who are infected with Mycobacterium Tuberculosis (Mtb) develop disease. Interferon-γ (IFN-γ) specific T cell immune assays such as Quantiferon and Elispot, as well as a skin hypersensitivity test, known as a tuberculin skin test, are widely used to infer infection. These assays measure immune conversion in response to Mtb. Some individuals measure persistently negative to immune conversion, despite high and prolonged exposure to Mtb. Increasing interest into this phenotype has led to multiple publications describing various aspects of these responses. However, there is a lack of a unified "resister" definition. A universal definition will improve cross study data comparisons and assist with future study design and planning. We review the current literature describing this phenotype and make recommendations for future studies.

HDAC3 inhibitor RGFP966 controls bacterial growth and modulates macrophage signaling during Mycobacterium tuberculosis infection

RATIONALE

Host-directed therapeutics for Mycobacterium tuberculosis (Mtb) offer potential strategies for combatting antibiotic resistance and for killing non-replicating bacilli. Phenylbutyrate, a partially selective histone-deacetylase (HDAC) inhibitor, was previously shown to control Mtb growth and alter macrophage inflammatory pathways at 2-4 mM concentrations.

OBJECTIVE

To identify a more potent and selective HDAC inhibitor that modulates macrophage responses to mycobacteria and has direct antibacterial effects against Mtb.

METHODS

We used cellular approaches to characterize the role of pharmacologic inhibition of HDAC3 on Mtb growth and Mtb-induced peripheral and alveolar macrophage immune functions.

MEASUREMENTS AND MAIN RESULTS

RGFP966, an HDAC3 inhibitor, controlled Mtb, BCG and M. avium growth directly in broth culture and in human peripheral blood monocyte-derived and alveolar macrophages with an MIC50 of approximately 5-10 μM. In contrast, RGFP966 did not inhibit growth of several other intracellular and extracellular bacteria. We also found that RGFP966 modulated macrophage pro-inflammatory cytokine secretion in response to Mtb infection with decreased IL6 and TNF secretion.

CONCLUSIONS

We identified a potent and selective small molecule inhibitor of HDAC3 with direct antimicrobial activity against Mtb and modulation of macrophage signaling pathways.

Genetics and evolution of tuberculosis pathogenesis: New perspectives and approaches

Tuberculosis is the most lethal infectious disease globally, but the vast majority of people who are exposed to the primary causative pathogen, Mycobacterium tuberculosis (MTB), do not develop active disease. Most people do, however, show signs of infection that remain throughout their lifetimes. In this review, we develop a framework that describes several possible transitions from pathogen exposure to TB disease and reflect on the genetics studies to address many of these. The evidence strongly supports a human genetic component for both infection and active disease, but many of the existing studies, including some of our own, do not clearly delineate what transition(s) is being explicitly examined. This can make interpretation difficult in terms of why only some people develop active disease. Nonetheless, both linkage peaks and associations with either active disease or latent infection have been identified. For transition to active disease, pathways defined as active TB altered T and B cell signaling in rheumatoid arthritis and T helper cell differentiation are significantly associated. Pathways that affect transition from exposure to infection are less clear-cut, as studies of this phenotype are less common, and a primary response, if it exists, is not yet well defined. Lastly, we discuss the role that interaction between the MTB lineage and human genetics can play in TB disease, especially severity. Severity of TB is at present the only way to study putative co-evolution between MTB and humans as it is impossible in the absence of disease to know the MTB lineage(s) to which an individual has been exposed. In addition, even though severity has been defined in multiple heterogeneous ways, it appears that MTB-human co-evolution may shape pathogenicity. Further analysis of co-evolution, requiring careful analysis of paired samples, may be the best way to completely assess the genetic basis of TB.

Interaction between host genes and Mycobacterium tuberculosis lineage can affect tuberculosis severity: Evidence for coevolution?

Genetic studies of both the human host and Mycobacterium tuberculosis (MTB) demonstrate independent association with tuberculosis (TB) risk. However, neither explains a large portion of disease risk or severity. Based on studies in other infectious diseases and animal models of TB, we hypothesized that the genomes of the two interact to modulate risk of developing active TB or increasing the severity of disease, when present. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which there were 3 MTB lineages of which L4-Ugandan (L4.6) is the most recent. TB severity, measured using the Bandim TBscore, was modeled as a function of host SNP genotype, MTB lineage, and their interaction, within two independent cohorts of TB cases, N = 113 and 121. No association was found between lineage and severity, but association between multiple polymorphisms in IL12B and TBscore was replicated in two independent cohorts (most significant rs3212227, combined p = 0.0006), supporting previous associations of IL12B with TB susceptibility. We also observed significant interaction between a single nucleotide polymorphism (SNP) in SLC11A1 and the L4-Ugandan lineage in both cohorts (rs17235409, meta p = 0.0002). Interestingly, the presence of the L4-Uganda lineage in the presence of the ancestral human allele associated with more severe disease. These findings demonstrate that IL12B is associated with severity of TB in addition to susceptibility, and that the association between TB severity and human genetics can be due to an interaction between genes in the two species, consistent with host-pathogen coevolution in TB.

Genetic Resistance to Mycobacterium tuberculosis Infection and Disease

Natural history studies of tuberculosis (TB) have revealed a spectrum of clinical outcomes after exposure to Mycobacterium tuberculosis, the cause of TB. Not all individuals exposed to the bacterium will become diseased and depending on the infection pressure, many will remain infection-free. Intriguingly, complete resistance to infection is observed in some individuals (termed resisters) after intense, continuing M. tuberculosis exposure. After successful infection, the majority of individuals will develop latent TB infection (LTBI). This infection state is currently (and perhaps imperfectly) defined by the presence of a positive tuberculin skin test (TST) and/or interferon gamma release assay (IGRA), but no detectable clinical disease symptoms. The majority of healthy individuals with LTBI are resistant to clinical TB, indicating that infection is remarkably well-contained in these non-progressors. The remaining 5-15% of LTBI positive individuals will progress to active TB. Epidemiological investigations have indicated that the host genetic component contributes to these infection and disease phenotypes, influencing both susceptibility and resistance. Elucidating these genetic correlates is therefore a priority as it may translate to new interventions to prevent, diagnose or treat TB. The most successful approaches in resistance/susceptibility investigation have focused on specific infection and disease phenotypes and the resister phenotype may hold the key to the discovery of actionable genetic variants in TB infection and disease. This review will not only discuss lessons from epidemiological studies, but will also focus on the contribution of epidemiology and functional genetics to human genetic resistance to M. tuberculosis infection and disease.