Search

Abstract Background Despite the success of antiretroviral therapy (ART), people living with HIV (PLWH) suffer from a high burden of pulmonary diseases, even after accounting for their smoking status. Cytotoxic CD8 T-cells are likely implicated in this phenomenon and may act as a double-edged sword. While being essential in viral infection control, their hyperactivation can also contribute to lung mucosal tissue damage. The effects of HIV and smoking on pulmonary mucosal CD8 T-cell dynamics has been a neglected area of research, which we address herein. Methods Bronchoalveolar lavage (BAL) fluid were obtained from ART-treated PLWH (median duration of supressed viral load: 9 years; smokers: n = 14; non-smokers: n = 21) and HIV-uninfected controls (smokers: n = 11; non-smokers: n = 20) without any respiratory symptoms or active infection. Lymphocytes were isolated and CD8 T-cell subsets and homing markers were characterized by multiparametric flow cytometry. Results Both smoking and HIV infection were independently associated with a significant increase in frequencies of total pulmonary mucosal CD8 T-cell. BAL CD8 T-cells were primarily CD69 + expressing CD103 and/or CD49a, at least one of the two granzymes (GzmA/GzmB), and little Perforin. Higher expression levels of CD103, CD69, and GzmB were observed in smokers versus non-smokers. The ex vivo phenotype of GzmA + and GzmB + cells revealed increased expression of CD103 and CXCR6 in smokers, while PLWH displayed elevated levels of CX3CR1 compared to controls. Conclusion Smoking and HIV could promote cytotoxic CD8 T-cell retention in small airways through different mechanisms. Smoking likely increases recruitment and retention of GzmB + CD8 Trm via CXCR6 and CD103. Heightened CX3CR1 expression could be associated with CD8 non-Trm recruitment from the periphery in PLWH.

Abstract Predicting COVID-19 severity is difficult, and the biological pathways involved are not fully understood. To approach this problem, we measured 4701 circulating human protein abundances in two independent cohorts totaling 986 individuals. We then trained prediction models including protein abundances and clinical risk factors to predict COVID-19 severity in 417 subjects and tested these models in a separate cohort of 569 individuals. For severe COVID-19, a baseline model including age and sex provided an area under the receiver operator curve (AUC) of 65% in the test cohort. Selecting 92 proteins from the 4701 unique protein abundances improved the AUC to 88% in the training cohort, which remained relatively stable in the testing cohort at 86%, suggesting good generalizability. Proteins selected from different COVID-19 severity were enriched for cytokine and cytokine receptors, but more than half of the enriched pathways were not immune-related. Taken together, these findings suggest that circulating proteins measured at early stages of disease progression are reasonably accurate predictors of COVID-19 severity. Further research is needed to understand how to incorporate protein measurement into clinical care.

Latent tuberculosis infection (LTBI) is common in people living with HIV (PLHIV) in high-TB-burden settings. Active TB is associated with specific stool taxa; however, little is known about the stool microbiota and LTBI in PLHIV. We characterised the stool microbiota of PLHIV with [interferon-γ release assay (IGRA)- and tuberculin skin test (TST)-positive] or without (IGRA- and TST-negative) LTBI (n = 25 per group). The 16S rRNA DNA sequences were analysed using QIIME2, Dirichlet-Multinomial Mixtures, DESeq2, and PICRUSt2. No α- or β-diversity differences occurred by LTBI status; however, LTBI-positive people were Faecalibacterium-, Blautia-, Gemmiger-, and Bacteroides-enriched and Moryella-, Atopobium-, Corynebacterium-, and Streptococcus-depleted. Inferred metagenome data showed that LTBI-negative-enriched pathways included several metabolite degradation pathways. Stool from LTBI-positive people demonstrated differential taxa abundance based on a quantitative response to antigen stimulation. In LTBI-positive people, older people had different β-diversities than younger people, whereas in LTBI-negative people, no differences occurred across age groups. Amongst female PLHIV, those with LTBI were, vs. those without LTBI, Faecalibacterium-, Blautia-, Gemmiger-, and Bacteriodes-enriched, which are producers of short-chain fatty acids. Taxonomic differences amongst people with LTBI occurred according to quantitative response to antigen stimulation and age. These data enhance our understanding of the microbiome’s potential role in LTBI.

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.

Pulmonary macrophages have two distinct ontogenies: long-lived embryonically-seeded alveolar macrophages (AM) and bone marrow-derived macrophages (BMDM). Here, we show that after infection with a virulent strain of Mycobacterium tuberculosis (H37Rv), primary murine AM exhibit a unique transcriptomic signature characterized by metabolic reprogramming distinct from conventional BMDM. In contrast to BMDM, AM failed to shift from oxidative phosphorylation (OXPHOS) to glycolysis and consequently were unable to control infection with an avirulent strain (H37Ra). Importantly, healthy human AM infected with H37Ra equally demonstrated diminished energetics, recapitulating our observation in the murine model system. However, the results from seahorse showed that the shift towards glycolysis in both AM and BMDM was inhibited by H37Rv. We further demonstrated that pharmacological (e.g. metformin or the iron chelator desferrioxamine) reprogramming of AM towards glycolysis reduced necrosis and enhanced AM capacity to control H37Rv growth. Together, our results indicate that the unique bioenergetics of AM renders these cells a perfect target for Mtb survival and that metabolic reprogramming may be a viable host targeted therapy against TB.

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.

Human genetic control is thought to affect a considerable part of the outcome of infection with Mycobacterium tuberculosis (Mtb). Most of us deal with the pathogen by containment (associated with clinical “latency”) or sterilization, but tragically millions each year do not. After decades of studies on host genetic susceptibility to Mtb infection, genetic variation has been discovered to play a role in tuberculous immunoreactivity and tuberculosis (TB) disease. Genes encoding pattern recognition receptors (PRRs) enable a consistent, molecularly direct interaction between humans and Mtb which suggests the potential for co-evolution. In this review, we explore the roles ascribed to PRRs during Mtb infection and ask whether such a longstanding and intimate interface between our immune system and this pathogen plays a critical role in determining the outcome of Mtb infection. The scientific evidence to date suggests that PRR variation is clearly implicated in altered immunity to Mtb but has a more subtle role in limiting the pathogen and pathogenesis. In contrast to ‘effectors’ like IFN-γ, IL-12, Nitric Oxide and TNF that are critical for Mtb control, ‘sensors’ like PRRs are less critical for the outcome of Mtb infection. This is potentially due to redundancy of the numerous PRRs in the innate arsenal, such that Mtb rarely goes unnoticed. Genetic association studies investigating PRRs during Mtb infection should therefore be designed to investigate endophenotypes of infection – such as immunological or clinical variation – rather than just TB disease, if we hope to understand the molecular interface between innate immunity and Mtb.

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.

N/A

Background: Mycobacterium tuberculosis (Mtb) infection is inferred from positive results of T-cell immune conversion assays measuring Mtb-specific interferon gamma production or tuberculin skin test (TST) reactivity. Certain exposed individuals do not display T-cell immune conversion in these assays and do not develop TB. Here we report a hitherto unknown form of this phenotype: HIV-1-positive persistently TB, tuberculin and IGRA negative (HITTIN). Methods: A community-based case-control design was used to systematically screen and identify adults living with HIV (HIV+), aged 35–60 years, who met stringent study criteria, and then longitudinally followed up for repeat IGRA and TST testing. Participants had no history of TB despite living in TB hyper-endemic environments in Cape Town, South Africa with a provincial incidence of 681/100,000. Mtb-specific antibodies were measured using ELISA and Luminex. Findings: We identified 48/286 (17%) individuals who tested persistently negative for Mtb-specific T-cell immunoreactivity (three negative Quantiferon results and one TST = 0mm) over 206±154 days on average. Of these, 97·2% had documented CD4 counts

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).

Tuberculosis (TB), caused by the human pathogens Mycobacterium tuberculosis (Mtb) and Mycobacterium africanum, has plagued humanity for millennia and remains the deadliest infectious disease in the modern world. Mycobacterium tuberculosis and M. africanum can be subdivided phylogenetically into seven lineages exhibiting a low but significant degree of genomic diversity and preferential geographic distributions. Human genetic variability impacts all stages of TB pathogenesis ranging from susceptibility to infection with Mtb, progression of infection to disease, and the development of distinct clinical subtypes. The genetic study of severe childhood TB identified strong inborn single-gene errors revealing crucial pathways of vulnerability to TB. However, the identification of major TB-susceptibility genes on the population level has remained elusive. In particular, the replication of findings from candidate and genome-wide association studies across distinct human populations has proven difficult, thus hampering the characterization of reliable host molecular markers of susceptibility. Among the possible confounding factors of genetic association studies is Mtb genomic variability, which generally was not taken into account by human genetic studies. In support of this possibility, Mtb lineage was found to be a contributing factor to clinical presentation of TB and epidemiological spread of Mtb in exposed populations. The confluence of pathogen and human host genetic variability to TB pathogenesis led to the consideration of a possible coadaptation of Mtb strains and their human hosts, which should reveal itself in significant interaction effects between Mtb strain and TB-susceptibility/resistance alleles. Here, we present some of the most consistent findings of genetic susceptibility factors in human TB and review studies that point to genome-to-genome interaction between humans and Mtb lineages. The limited results available so far suggest that analyses considering joint human–Mtb genomic variability may provide improved power for the discovery of pathogenic drivers of the ongoing TB epidemic.

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.

Certain individuals are able to resist Mycobacterium tuberculosis infection despite persistent and intense exposure. These persons do not exhibit adaptive immune priming as measured by tuberculin skin test (TST) and interferon-γ (IFN-γ) release assay (IGRA) responses, nor do they develop active tuberculosis (TB). Genetic investigation of individuals who are able to resist M. tuberculosis infection shows there are likely a combination of genetic variants that contribute to the phenotype. The contribution of the innate immune system and the exact cells involved in this phenotype remain incompletely elucidated. Neutrophils are prominent candidates for possible involvement as primers for microbial clearance. Significant variability is observed in neutrophil gene expression and DNA methylation. Furthermore, inter-individual variability is seen between the mycobactericidal capacities of donor neutrophils. Clearance of M. tuberculosis infection is favored by the mycobactericidal activity of neutrophils, apoptosis, effective clearance of cells by macrophages, and resolution of inflammation. In this review we will discuss the different mechanisms neutrophils utilize to clear M. tuberculosis infection. We discuss the duality between neutrophils' ability to clear infection and how increasing numbers of neutrophils contribute to active TB severity and mortality. Further investigation into the potential role of neutrophils in innate immune-mediated M. tuberculosis infection resistance is warranted since it may reveal clinically important activities for prevention as well as vaccine and treatment development.

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.

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.

A variety of host immunogenetic factors appear to influence both an individual's susceptibility to infection with Mycobacterium leprae and the pathologic course of the disease. Animal models can contribute to a better understanding of the role of immunogenetics in leprosy through comparative studies helping to confirm the significance of various identified traits and in deciphering the underlying mechanisms that may be involved in expression of different disease related phenotypes. Genetically engineered mice, with specific immune or biochemical pathway defects, are particularly useful for investigating granuloma formation and resistance to infection and are shedding new light on borderline areas of the leprosy spectrum which are clinically unstable and have a tendency toward immunological complications. Though armadillos are less developed in this regard, these animals are the only other natural hosts of M. leprae and they present a unique opportunity for comparative study of genetic markers and mechanisms associable with disease susceptibility or resistance, especially the neurological aspects of leprosy. In this paper, we review the recent contributions of genetically engineered mice and armadillos toward our understanding of the immunogenetics of leprosy.

Type-1 (T1R) and Type-2 (T2R) leprosy reactions (LR), which affect up to 50% of leprosy patients, are aggressive inflammatory episodes of sudden onset and highly variable incidence across populations. LR are often diagnosed concurrently with leprosy, but more frequently occur several months after treatment onset. It is not uncommon for leprosy patients to develop recurring reactional episodes; however, they rarely undergo both types of LR. Today, LR are the main cause of permanent disabilities associated with leprosy and represent a major challenge in the clinical management of leprosy patients. Although progress has been made in understanding the immunopathology of LR, the factors that cause a leprosy patient to suffer from LR are largely unknown. Given the impact that ethnic background has on the risk of developing LR, host genetic factors have long been suspected of contributing to LR. Indeed, polymorphisms in seven genes [Toll-like receptors (TLR)1, TLR2, nucleotide-binding oligomerisation domain containing 2, vitamin D receptor, natural resistance-associated macrophage protein 1, C4B and interleukin-6] have been found to be associated with one or more LR outcomes. The identification of host genetic markers with predictive value for LR would have a major impact on nerve damage control in leprosy. In this review, we present the recent advances achieved through genetic studies of LR.

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.