Search

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn's disease and Parkinson's disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.

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-8, 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-9). 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.

Leprosy is a chronic disease caused by Mycobacterium leprae. Worldwide, more than 200,000 new patients are affected by leprosy annually, making it the second most common mycobacterial disease after tuberculosis. The MHC/HLA region has been consistently identified as carrying major leprosy susceptibility variants in different populations at times with inconsistent results. To establish the unambiguous molecular identity of classical HLA class I and class II leprosy susceptibility factors, we applied next-generation sequencing to genotype with high-resolution 11 HLA class I and class II genes in 1,155 individuals from a Vietnamese leprosy case-control sample. HLA alleles belonging to an extended haplotype from HLA-A to HLA-DPB1 were associated with risk to leprosy. This susceptibility signal could be reduced to the HLA-DRB1*10:01~ HLA-DQA1*01:05 alleles which were in complete linkage disequilibrium (LD). In addition, haplotypes containing HLA-DRB3~ HLA-DRB1*12:02 and HLA-C*07:06~ HLA-B*44:03~ HLA-DRB1*07:01 alleles were found as two independent protective factors for leprosy. Moreover, we replicated the previously associated HLA-DRB1*15:01 as leprosy risk factor and HLA-DRB1*04:05~HLA-DQA1*03:03 as protective alleles. When we narrowed the analysis to the single amino acid level, we found that the associations of the HLA alleles were largely captured by four independent amino acids at HLA-DRβ1 positions 57 (D) and 13 (F), HLA-B position 63 (E) and HLA-A position 19 (K). Hence, analyses at the amino acid level circumvented the ambiguity caused by strong LD of leprosy susceptibility HLA alleles and identified four distinct leprosy susceptibility factors.

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

Complex traits are known to be influenced by a combination of environmental factors and rare and common genetic variants. However, detection of such multivariate associations can be compromised by low statistical power and confounding by population structure. Linear mixed effects models (LMM) can account for correlations due to relatedness but have not been applicable in high-dimensional (HD) settings where the number of fixed effect predictors greatly exceeds the number of samples. False positives or false negatives can result from two-stage approaches, where the residuals estimated from a null model adjusted for the subjects' relationship structure are subsequently used as the response in a standard penalized regression model. To overcome these challenges, we develop a general penalized LMM with a single random effect called ggmix for simultaneous SNP selection and adjustment for population structure in high dimensional prediction models. We develop a blockwise coordinate descent algorithm with automatic tuning parameter selection which is highly scalable, computationally efficient and has theoretical guarantees of convergence. Through simulations and three real data examples, we show that ggmix leads to more parsimonious models compared to the two-stage approach or principal component adjustment with better prediction accuracy. Our method performs well even in the presence of highly correlated markers, and when the causal SNPs are included in the kinship matrix. ggmix can be used to construct polygenic risk scores and select instrumental variables in Mendelian randomization studies. Our algorithms are available in an R package available on CRAN (https://cran.r-project.org/package=ggmix).

Increased prevalence of latent tuberculosis infection (LTBI) has been observed among high-risk populations such as healthcare workers (HCWs). The results may depend on the method of LTBI assessment, interferon-gamma release assay (IGRA) and/or tuberculin skin test (TST). Here, we investigated the prevalence and risk factors for LTBI assessed by both IGRAs and TST in HCWs living in Morocco, a country with intermediate tuberculosis (TB) endemicity and high BCG vaccination coverage. HCWs were recruited in two Moroccan hospitals, Rabat and Meknes. All the participants underwent testing for LTBI by both IGRA (QuantiFERON-TB Gold In-Tube, QFT-GIT) and TST. Different combinations of IGRA and TST results defined the LTBI status. Risk factors associated with LTBI were investigated using a mixed-effect logistic regression model. The prevalence of LTBI among 631 HCWs (age range 18-60 years) varied from 40.7% (95%CI 36.9-44.5%) with QFT-GIT to 52% (95%CI 48.2-56.0%) with TST using a 10 mm cut-off. The highest agreement between QFT-GIT and TST (κ = 0.50; 95%CI 0.43-0.56) was observed with the 10 mm cut-off for a positive TST. For a definition of LTBI status using a double positive result for both QFT-GIT and TST, significant associations were found with the following risk factors: being male (OR = 2.21; 95%CI 1.40-3.49; p = 0.0007), belonging to age groups 35-44 years (OR = 2.43; 95%CI 1.45-4.06; p = 0.0007) and even more 45-60 years (OR = 4.81; 95%CI 2.72-8.52; p = 7.10-8), having a family history of TB (OR = 6.62; 95%CI 2.59-16.94; p = 8.10-5), and working at a pulmonology unit (OR = 3.64; 95%CI 1.44-9.23; p = 0.006). Smoking was associated with LTBI status when defined by a positive QFT-GIT result (OR = 1.89; 95%CI 1.12-3.21; p = 0.02). A high prevalence of LTBI was observed among HCWs in two Moroccan hospitals. Male gender, increased age, family history of TB, and working at a pulmonology unit were consistent risk factors associated with LTBI.

Leprosy is a human infectious disease caused by Mycobacterium leprae. A strong host genetic contribution to leprosy susceptibility is well established. However, the modulation of the transcriptional response to infection and the mechanism(s) of disease control are poorly understood. To address this gap in knowledge of leprosy pathogenicity, we conducted a genome-wide search for expression quantitative trait loci (eQTL) that are associated with transcript variation before and after stimulation with M. leprae sonicate in whole blood cells. We show that M. leprae antigen stimulation mainly triggered the upregulation of immune related genes and that a substantial proportion of the differential gene expression is genetically controlled. Indeed, using stringent criteria, we identified 318 genes displaying cis-eQTL at an FDR of 0.01, including 66 genes displaying response-eQTL (reQTL), i.e. cis-eQTL that showed significant evidence for interaction with the M. leprae stimulus. Such reQTL correspond to regulatory variations that affect the interaction between human whole blood cells and M. leprae sonicate and, thus, likely between the human host and M. leprae bacilli. We found that reQTL were significantly enriched among binding sites of transcription factors that are activated in response to infection, and that they were enriched among single nucleotide polymorphisms (SNPs) associated with susceptibility to leprosy per se and Type-I Reaction, and seven of them have been targeted by recent positive selection. Our study suggested that natural selection shaped our genomic diversity to face pathogen exposure including M. leprae infection.

Leprosy Type-1 Reactions (T1Rs) are pathological inflammatory responses that afflict a sub-group of leprosy patients and result in peripheral nerve damage. Here, we employed a family-based GWAS in 221 families with 229 T1R-affect offspring with stepwise replication to identify risk factors for T1R. We discovered, replicated and validated T1R-specific associations with SNPs located in chromosome region 10p21.2. Combined analysis across the three independent samples resulted in strong evidence of association of rs1875147 with T1R (p = 4.5x10-8; OR = 1.54, 95% CI = 1.32-1.80). The T1R-risk locus was restricted to a lncRNA-encoding genomic interval with rs1875147 being an eQTL for the lncRNA. Since a genetic overlap between leprosy and inflammatory bowel disease (IBD) has been detected, we evaluated if the shared genetic control could be traced to the T1R endophenotype. Employing the results of a recent IBD GWAS meta-analysis we found that 10.6% of IBD SNPs available in our dataset shared a common risk-allele with T1R (p = 2.4x10-4). This finding points to a substantial overlap in the genetic control of clinically diverse inflammatory disorders.

After sustained exposure to Mycobacterium leprae, only a subset of exposed individuals develops clinical leprosy. Moreover, leprosy patients show a wide spectrum of clinical manifestations that extend from the paucibacillary (PB) to the multibacillary (MB) form of the disease. This "polarization" of leprosy has long been a major focus of investigation for immunologists because of the different immune response in these two forms. But while leprosy per se has been shown to be under tight human genetic control, few epidemiological or genetic studies have focused on leprosy subtypes. Using PubMed, we collected available data in English on the epidemiology of leprosy polarization and the possible role of human genetics in its pathophysiology until September 2015. At the genetic level, we assembled a list of 28 genes from the literature that are associated with leprosy subtypes or implicated in the polarization process. Our bibliographical search revealed that improved study designs are needed to identify genes associated with leprosy polarization. Future investigations should not be restricted to a subanalysis of leprosy per se studies but should instead contrast MB to PB individuals. We show the latter approach to be the most powerful design for the identification of genetic polarization determinants. Finally, we bring to light the important resource represented by the nine-banded armadillo model, a unique animal model for leprosy.

Depending on the epidemiological setting, a variable proportion of leprosy patients will suffer from excessive pro-inflammatory responses, termed type-1 reactions (T1R). The LRRK2 gene encodes a multi-functional protein that has been shown to modulate pro-inflammatory responses. Variants near the LRRK2 gene have been associated with leprosy in some but not in other studies. We hypothesized that LRRK2 was a T1R susceptibility gene and that inconsistent association results might reflect different proportions of patients with T1R in the different sample settings. Hence, we evaluated the association of LRRK2 variants with T1R susceptibility.An association scan of the LRRK2 locus was performed using 156 single-nucleotide polymorphisms (SNPs). Evidence of association was evaluated in two family-based samples: A set of T1R-affected and a second set of T1R-free families. Only SNPs significant for T1R-affected families with significant evidence of heterogeneity relative to T1R-free families were considered T1R-specific. An expression quantitative trait locus (eQTL) analysis was applied to evaluate the impact of T1R-specific SNPs on LRRK2 gene transcriptional levels.A total of 18 T1R-specific variants organized in four bins were detected. The core SNP capturing the T1R association was the LRRK2 missense variant M2397T (rs3761863) that affects LRRK2 protein turnover. Additionally, a bin of nine SNPs associated with T1R were eQTLs for LRRK2 in unstimulated whole blood cells but not after exposure to Mycobacterium leprae antigen.The results support a preferential association of LRRK2 variants with T1R. LRRK2 involvement in T1R is likely due to a pathological pro-inflammatory loop modulated by LRRK2 availability. Interestingly, the M2397T variant was reported in association with Crohn's disease with the same risk allele as in T1R suggesting common inflammatory mechanism in these two distinct diseases.

In an accident later known as the Lübeck disaster, 251 neonates were orally given three doses of the new Bacille Calmette-Guérin (BCG) antituberculosis (TB) vaccine contaminated with Mycobacterium tuberculosis. A total of 173 infants developed clinical or radiological signs of TB but survived the infection, while 72 died from TB. While some blamed the accident on BCG itself by postulating reversion to full virulence, such a possibility was conclusively disproven. Rather, by combining clinical, microbiological, and epidemiological data, the chief public health investigator Dr. A. Moegling concluded that the BCG vaccine had been contaminated with variable amounts of fully virulent M. tuberculosis. Here, we summarize the conclusions drawn by Moegling and point out three lessons that can be learned. First, while mortality was high (approximately 29%), the majority of neonates inoculated with M. tuberculosis eventually overcame TB disease. This shows the high constitutional resistance of humans to the bacillus. Second, four semiquantitative levels of contamination were deduced by Moegling from the available data. While at low levels of M. tuberculosis there was a large spread of clinical phenotypes reflecting a good degree of innate resistance to TB, at the highest dose, the majority of neonates were highly susceptible to TB. This shows the dominating role of dose for innate resistance to TB. Third, two infants inoculated with the lowest dose nevertheless died of TB, and their median time from inoculation to death was substantially shorter than for those who died after inoculation with higher doses. This suggests that infants who developed disease after low dose inoculation are those who are most susceptible to the disease. We discuss some implications of these lessons for current study of genetic susceptibility to TB.

Natural Killer (NK) cells contribute to the control of viral infection by directly killing target cells and mediating cytokine release. In C57BL/6 mice, the Ly49H activating NK cell receptor plays a key role in early resistance to mouse cytomegalovirus (MCMV) infection through specific recognition of the MCMV-encoded MHC class I-like molecule m157 expressed on infected cells. Here we show that transgenic expression of Ly49H failed to provide protection against MCMV infection in the naturally susceptible A/J mouse strain. Characterization of Ly49H(+) NK cells from Ly49h-A transgenic animals showed that they were able to mount a robust cytotoxic response and proliferate to high numbers during the course of infection. However, compared to NK cells from C57BL/6 mice, we observed an intrinsic defect in their ability to produce IFNγ when challenged by either m157-expressing target cells, exogenous cytokines or chemical stimulants. This effect was limited to NK cells as T cells from C57BL/6 and Ly49h-A mice produced comparable cytokine levels. Using a panel of recombinant congenic strains derived from A/J and C57BL/6 progenitors, we mapped the genetic basis of defective IFNγ production to a single 6.6 Mb genetic interval overlapping the Ifng gene on chromosome 10. Inspection of the genetic interval failed to reveal molecular differences between A/J and several mouse strains showing normal IFNγ production. The chromosome 10 locus is independent of MAPK signalling or decreased mRNA stability and linked to MCMV susceptibility. This study highlights the existence of a previously uncovered NK cell-specific cis-regulatory mechanism of Ifnγ transcript expression potentially relevant to NK cell function in health and disease.

Leprosy is a persistent infectious disease caused by Mycobacterium leprae that still affects over 200,000 new patients annually. The host genetic background is an important risk factor for leprosy susceptibility and the PARK2 gene is a replicated leprosy susceptibility candidate gene. The protein product of PARK2, Parkin, is an E3 ubiquitin ligase that is involved in the development of various forms of Parkinsonism. The human macrophage is both a natural host cell of M. leprae as well as a primary mediator of natural immune defenses, in part by secreting important pro-inflammatory cytokines and chemokines. Here, we report that down-regulation of Parkin in THP-1 macrophages, human monocyte-derived macrophages and human Schwann cells resulted in a consistent and specific decrease in interleukin-6 (IL-6) and monocyte chemoattractant protein 1 (MCP-1/CCL2) production in response to mycobacteria or LPS. Interestingly, production of IL-6 at 6 hours by THP-1 cells stimulated with live M. leprae and M. bovis BCG was dependent on pretreatment with 1,25-dihydroxyvitamin D(3) (VD). Parkin knockdown in VD-treated cells blocked IL-6 induction by mycobacteria. However, IκB-α phosphorylation and levels of IκB-ξ, a nuclear protein required for IL-6 expression, were not affected by Parkin silencing. Phosphorylation of MAPK ERK1/2 and p38 was unaffected by Parkin silencing while JNK activation was promoted but did not explain the altered cytokine production. In a final set of experiments we found that genetic risk factors of leprosy located in the PARK2 promoter region were significantly correlated with M. leprae sonicate triggered CCL2 and IL6 transcript levels in whole blood assays. These results associated genetically controlled changes in the production of MCP-1/CCL2 and IL-6 with known leprosy susceptibility factors.

Leprosy reversal reactions type 1 (T1R) are acute immune episodes that affect a subset of leprosy patients and remain a major cause of nerve damage. Little is known about the relative importance of innate versus environmental factors in the pathogenesis of T1R. In a retrospective design, we evaluated innate differences in response to Mycobacterium leprae between healthy individuals and former leprosy patients affected or free of T1R by analyzing the transcriptome response of whole blood to M. leprae sonicate. Validation of results was conducted in a subsequent prospective study. We observed the differential expression of 581 genes upon exposure of whole blood to M. leprae sonicate in the retrospective study. We defined a 44 T1R gene set signature of differentially regulated genes. The majority of the T1R set genes were represented by three functional groups: i) pro-inflammatory regulators; ii) arachidonic acid metabolism mediators; and iii) regulators of anti-inflammation. The validity of the T1R gene set signature was replicated in the prospective arm of the study. The T1R genetic signature encompasses genes encoding pro- and anti-inflammatory mediators of innate immunity. This suggests an innate defect in the regulation of the inflammatory response to M. leprae antigens. The identified T1R gene set represents a critical first step towards a genetic profile of leprosy patients who are at increased risk of T1R and concomitant nerve damage.

The host response to mycobacterial infection depends on host and pathogen genetic factors. Recent studies in human populations suggest a strain specific genetic control of tuberculosis. To test for mycobacterial-strain specific genetic control of susceptibility to infection under highly controlled experimental conditions, we performed a comparative genetic analysis using the A/J- and C57BL/6J-derived recombinant congenic (RC) mouse panel infected with the Russia and Pasteur strains of Mycobacterium bovis Bacille Calmette Guérin (BCG). Bacillary counts in the lung and spleen at weeks 1 and 6 post infection were used as a measure of susceptibility. By performing genome-wide linkage analyses of loci that impact on tissue-specific bacillary burden, we were able to show the importance of correcting for strain background effects in the RC panel. When linkage analysis was adjusted on strain background, we detected a single locus on chromosome 11 that impacted on pulmonary counts of BCG Russia but not Pasteur. The same locus also controlled the splenic counts of BCG Russia but not Pasteur. By contrast, a locus on chromosome 1 which was indistinguishable from Nramp1 impacted on splenic bacillary counts of both BCG Russia and Pasteur. Additionally, dependent upon BCG strain, tissue and time post infection, we detected 9 distinct loci associated with bacillary counts. Hence, the ensemble of genetic loci impacting on BCG infection revealed a highly dynamic picture of genetic control that reflected both the course of infection and the infecting strain. This high degree of adaptation of host genetics to strain-specific pathogenesis is expected to provide a suitable framework for the selection of specific host-mycobacteria combinations during co-evolution of mycobacteria with humans.