Linkage disequilibrium pattern and age-at-diagnosis are critical for replicating genetic associations across ethnic groups in leprosy

Machine learning approaches to identify systemic lupus erythematosus in anti-nuclear antibody-positive patients using genomic data and electronic health records

BACKGROUND

Although the 2019 EULAR/ACR classification criteria for systemic lupus erythematosus (SLE) has required at least a positive anti-nuclear antibody (ANA) titer (≥ 1:80), it remains challenging for clinicians to identify patients with SLE. This study aimed to develop a machine learning (ML) approach to assist in the detection of SLE patients using genomic data and electronic health records.

METHODS

Participants with a positive ANA (≥ 1:80) were enrolled from the Taiwan Precision Medicine Initiative cohort. The Taiwan Biobank version 2 array was used to detect single nucleotide polymorphism (SNP) data. Six ML models, Logistic Regression, Random Forest (RF), Support Vector Machine, Light Gradient Boosting Machine, Gradient Tree Boosting, and Extreme Gradient Boosting (XGB), were used to identify SLE patients. The importance of the clinical and genetic features was determined by Shapley Additive Explanation (SHAP) values. A logistic regression model was applied to identify genetic variations associated with SLE in the subset of patients with an ANA equal to or exceeding 1:640.

RESULTS

A total of 946 SLE and 1,892 non-SLE controls were included in this analysis. Among the six ML models, RF and XGB demonstrated superior performance in the differentiation of SLE from non-SLE. The leading features in the SHAP diagram were anti-double strand DNA antibodies, ANA titers, AC4 ANA pattern, polygenic risk scores, complement levels, and SNPs. Additionally, in the subgroup with a high ANA titer (≥ 1:640), six SNPs positively associated with SLE and five SNPs negatively correlated with SLE were discovered.

CONCLUSIONS

ML approaches offer the potential to assist in diagnosing SLE and uncovering novel SNPs in a group of patients with autoimmunity.

Investigation of PRKN Mutations in Levodopa-Induced Dyskinesia in Parkinson's Disease Treatment

Mitophagy is an important process that participates in mitochondrial quality control. Dysfunctions in this process can be caused by mutations in genes like PRKN and are associated with the development and progression of Parkinson's Disease (PD). The most used drug in the treatment of PD is levodopa (LD), but it can cause adverse effects, such as dyskinesia. Currently, few studies are searching for biomarkers for an effective use of lLD for this disease, especially regarding mitophagy genetics. Thus, this work investigates the association of 14 variants of the PRKN gene with LD in the treatment of PD. We recruited 70 patients with PD undergoing treatment with LD (39 without dyskinesia and 31 with dyskinesia). Genotyping was based on Sanger sequencing. Our results reinforce that age at onset of symptoms, duration of PD, and treatment and dosage of LD can influence the occurrence of dyskinesia but not the investigated PRKN variants. The perspective presented here of variants of mitophagy-related genes in the context of treatment with LD is still underexplored, although an association has been indicated in previous studies. We suggest that other variants in PRKN or in other mitophagy genes may participate in the development of levodopa-induced dyskinesia in PD treatment.

Human Genetic Variation Influences Enteric Fever Progression

In the 21st century, enteric fever is still causing a significant number of mortalities, especially in high-risk regions of the world. Genetic studies involving the genome and transcriptome have revealed a broad set of candidate genetic polymorphisms associated with susceptibility to and the severity of enteric fever. This review attempted to explain and discuss the past and the most recent findings on human genetic variants affecting the progression of Salmonella typhoidal species infection, particularly toll-like receptor (TLR) 4, TLR5, interleukin (IL-) 4, natural resistance-associated macrophage protein 1 (NRAMP1), VAC14, PARK2/PACRG, cystic fibrosis transmembrane conductance regulator (CFTR), major-histocompatibility-complex (MHC) class II and class III. These polymorphisms on disease susceptibility or progression in patients could be related to multiple mechanisms in eliminating both intracellular and extracellular Salmonella typhoidal species. Here, we also highlighted the limitations in the studies reported, which led to inconclusive results in association studies. Nevertheless, the knowledge obtained through this review may shed some light on the development of risk prediction tools, novel therapies as well as strategies towards developing a personalised typhoid vaccine.

Association of NOD2 and IFNG single nucleotide polymorphisms with leprosy in the Amazon ethnic admixed population

Leprosy is a chronic infectious disease, caused by Mycobacterium leprae, which affects skin and peripheral nerves. Polymorphisms in genes associated with autophagy, metabolism, innate and adaptive immunity confer susceptibility to leprosy. However, these associations need to be confirmed through independent replication studies in different ethnicities. The population from Amazon state (northern Brazil) is admixed and it contains the highest proportion of Native American genetic ancestry in Brazil. We conducted a case-control study for leprosy in which we tested fourteen previously associated SNPs in key immune response regulating genes: TLR1 (rs4833095), NOD2 (rs751271, rs8057341), TNF (rs1800629), IL10 (rs1800871), CCDC122/LACC1 (rs4942254), PACRG/PRKN (rs9356058, rs1040079), IFNG (rs2430561), IL6 (rs2069845), LRRK2 (rs7298930, rs3761863), IL23R (rs76418789) and TYK2 (rs55882956). Genotyping was carried out by allelic discrimination in 967 controls and 412 leprosy patients. Association with susceptibility was assessed by logistic regression analyses adjusted for the following covariates: gender, age and ancestry. Genetic ancestry was similar in case and control groups. Statistically significant results were only found for IFNG and NOD2. The rs8057341 polymorphism within NOD2 was identified as significant for the AA genotype (OR = 0.56; 95% CI, 0.37-0.84; P = 0.005) and borderline for the A allele (OR = 0.76; 95% CI, 0.58-1.00; P = 0.053) and carrier (OR = 0.76; 95% CI, 0.58-1.00; P = 0.051). The rs2430561 SNP in IFNG was associated with disease susceptibility for the AT genotype (OR = 1.40; 95% CI, 1.06-1.85; P = 0.018) and carrier (OR = 1.44; 95% CI, 1.10-1.88; P = 0.008). We confirmed that NOD2 and IFNG are major players in immunity against M.leprae in the Amazon ethnic admixed population.

The parkin-coregulated gene product PACRG promotes TNF signaling by stabilizing LUBAC

The Parkin-coregulated gene (PACRG), which encodes a protein of unknown function, shares a bidirectional promoter with Parkin (PRKN), which encodes an E3 ubiquitin ligase. Because PRKN is important in mitochondrial quality control and protection against stress, we tested whether PACRG also affected these pathways in various cultured human cell lines and in mouse embryonic fibroblasts. PACRG did not play a role in mitophagy but did play a role in tumor necrosis factor (TNF) signaling. Similarly to Parkin, PACRG promoted nuclear factor κB (NF-κB) activation in response to TNF. TNF-induced nuclear translocation of the NF-κB subunit p65 and NF-κB-dependent transcription were decreased in PACRG-deficient cells. Defective canonical NF-κB activation in the absence of PACRG was accompanied by a decrease in linear ubiquitylation mediated by the linear ubiquitin chain assembly complex (LUBAC), which is composed of the two E3 ubiquitin ligases HOIP and HOIL-1L and the adaptor protein SHARPIN. Upon TNF stimulation, PACRG was recruited to the activated TNF receptor complex and interacted with LUBAC components. PACRG functionally replaced SHARPIN in this context. In SHARPIN-deficient cells, PACRG prevented LUBAC destabilization, restored HOIP-dependent linear ubiquitylation, and protected cells from TNF-induced apoptosis. This function of PACRG in positively regulating TNF signaling may help to explain the association of PACRG and PRKN polymorphisms with an increased susceptibility to intracellular pathogens.

Human Genetic Susceptibility of Leprosy Recurrence

Host genetic susceptibility to leprosy has been intensively investigated over the last decades; however, there are no studies on the role of genetic variants in disease recurrence. A previous initiative identified three recurrent cases of leprosy for which none of the M. leprae strains, as obtained in the first and the second diagnosis, had any known genomic variants associated to resistance to Multidrug therapy; in addition, whole genome sequencing indicated that the same M. leprae was causing two out of the three recurrences. Thus, these individuals were suspected of being particularly susceptible to M. leprae infection, either as relapse or reinfection. To verify this hypothesis, 19 genetic markers distributed across 11 loci (14 genes) classically associated with leprosy were genotyped in the recurrent and in three matching non-recurrent leprosy cases. An enrichment of risk alleles was observed in the recurrent cases, suggesting the existence of a particularly high susceptibility genetic profile among leprosy patients predisposing to disease recurrence.

Genetics of leprosy: today and beyond

Leprosy is a chronic infectious disease of the skin and peripheral nerves that presents a strong link with the host genetic background. Different approaches in genetic studies have been applied to leprosy and today leprosy is among the infectious diseases with the greatest number of genetic risk variants identified. Several leprosy genes have been implicated in host immune response to pathogens and point to specific pathways that are relevant for host defense to infection. In addition, host genetic factors are also involved in the heterogeneity of leprosy clinical manifestations and in excessive inflammatory responses that occur in some leprosy patients. Finally, genetic studies in leprosy have provided strong evidence of pleiotropic effects between leprosy and other complex diseases, such as immune-mediated or neurodegenerative diseases. These findings not only impact on the field of leprosy and infectious diseases but also make leprosy a good model for the study of complex immune-mediated diseases. Here, we summarize recent genetic findings in leprosy susceptibility and discuss the overlap of the genetic control in leprosy with Parkinson's disease and inflammatory bowel disease. Moreover, some limitations, challenges, and potential new avenues for future genetics studies of leprosy are also discussed in this review.

IL8 and IL17A polymorphisms associated with multibacillary leprosy and reaction type 1 in a mixed population from southern Brazil

We evaluated the influence of the IL8 T-738A (nonidentified rs), IL8 T-353A (rs4073), IL17A G197A (rs2275913), and IL17F T7488C (rs763780) single-nucleotide polymorphisms on leprosy. The AA genotype of IL8 T-353A was observed as a risk factor for multibacillary leprosy, regardless of gender and age-of-onset of disease, considering the recessive model (OR, 3.8; 95% CI, 1.1-13.5; P, 0.023). Furthermore, the AA genotype of IL17A G197A was associated with leprosy type 1 reaction (OR, 2.4; 95% CI, 1.1-5.1; P, 0.026) when compared to the group without reaction, which was adjusted for gender and age-of-onset of disease by the model log additive. These results indicate association of IL8 and IL17A polymorphisms with the progression to multibacillary leprosy and with the type 1 reaction, respectively.

Human genetics of mycobacterial disease

Mycobacterial diseases are caused by members of the genus Mycobacterium, acid-fast bacteria characterized by the presence of mycolic acids within their cell walls. Claiming almost 2 million lives every year, tuberculosis (TB) is the most common mycobacterial disease and is caused by infection with M. tuberculosis and, in rare cases, by M. bovis or M. africanum. The second and third most common mycobacterial diseases are leprosy and buruli ulcer (BU), respectively. Both diseases affect the skin and can lead to permanent sequelae and deformities. Leprosy is caused by the uncultivable M. leprae while the etiological agent of BU is the environmental bacterium M. ulcerans. After exposure to these mycobacterial species, a majority of individuals will not progress to clinical disease and, among those who do, inter-individual variability in disease manifestation and outcome can be observed. Susceptibility to mycobacterial diseases carries a human genetic component and intense efforts have been applied over the past decades to decipher the exact nature of the genetic factors controlling disease susceptibility. While for BU this search was mostly conducted on the basis of candidate genes association studies, genome-wide approaches have been widely applied for TB and leprosy. In this review, we summarize some of the findings achieved by genome-wide linkage, association and transcriptome analyses in TB disease and leprosy and the recent genetic findings for BU susceptibility.

Genetic Susceptibility to Leprosy-From Classic Immune-Related Candidate Genes to Hypothesis-Free, Whole Genome Approaches

Genetics plays a crucial role in controlling susceptibility to infectious diseases by modulating the interplay between humans and pathogens. This is particularly evident in leprosy, since the etiological agent, Mycobacterium leprae, displays semiclonal characteristics not compatible with the wide spectrum of disease phenotypes. Over the past decades, genetic studies have unraveled several gene variants as risk factors for leprosy per se, disease clinical forms and the occurrence of leprosy reactions. As expected, several of these genes are immune-related; yet, hypothesis-free approaches have led to genes not classically linked to immune response. The PARK2, originally described as a Parkinson's disease gene, illustrates the case: Parkin-the protein coded by PARK2-was defined as an important player regulating innate and adaptive immune responses only years after its description as a leprosy susceptibility gene. Interestingly, even with the use of powerful hypothesis-free study designs such as genome-wide association studies, most of the major gene effect controlling leprosy susceptibility remains elusive. One hypothesis to explain this "hidden heritability" is that rare variants not captured by classic association studies are of critical importance. To address this question, massively parallel sequencing of large segments of the human genome-even whole exomes/genomes-is an alternative to properly identify rare, disease-causing mutations. These mutations may then be investigated through sophisticated approaches such as cell reprogramming and genome editing applied to create in vitro models for functional leprosy studies.

Genome-wide association study identifies loci associated with milk leukocyte phenotypes following experimental challenge with Streptococcus uberis

Mastitis is a detrimental disease in the dairy industry that decreases milk quality and costs upwards of $2 billion annually. Often, mastitis results from bacteria entering the gland through the teat opening. Streptococcus uberis is responsible for a high percentage of subclinical and clinical mastitis. Following an intramammary experimental challenge with S. uberis on Holstein cows (n = 40), milk samples were collected and somatic cell counts (SCC) were determined by the Dairy Herd Improvement Association Laboratory. Traditional genome-wide association studies (GWAS) have utilized test day SCC or SCC lactation averages to identify loci of interest. Our approach utilizes SCC collected following a S. uberis experimental challenge to generate three novel phenotypes: (1) area under the curve (AUC) of SCC for 0-7 days and (2) 0-28 days post-challenge; and (3) when SCC returned to below 200,000 cells/mL post-challenge (< 21 days, 21-28 days, or > 28 days). Polymorphisms were identified using Illumina's BovineSNP50 v2 DNA BeadChip. Associations were tested using Plink software and identified 16 significant (p < 1.0 × 10^-4) single-nucleotide polymorphisms (SNPs) across the phenotypes. Most significant SNPs were in genes linked to cell signaling, migration, and apoptosis. Several have been recognized in relation to infectious processes (ATF7, SGK1, and PACRG), but others less so (TRIO, GLRA1, CELSR2, TIAM2, CPE). Further investigation of these genes and their roles in inflammation (e.g., SCC) can provide potential targets that influence resolution of mammary gland infection. Likewise, further investigation of the identified SNP with mastitis and other disease phenotypes can provide greater insight to the potential of these SNP as genetic markers.

Rs2910164 in microRNA‑146a confers an elevated risk of depression in patients with coronary artery disease by modulating the expression of NOS1

Depression has been well established as an independent predictor of mortality and cardiac morbidity rates in patients with coronary artery disease (CAD). Evidence has shown that single nucleotide polymorphisms located in pre‑microRNA (miRNA) or mature miRNA may modify various biological processes and affect the process of carcinogenesis, and the downregulation of neuronal nitric oxide synthase 1 (NOS1) can induce depression. It has been shown that NOS1 is the target gene of miR‑146a, and that the rs2910164 G/C polymorphism can downregulate the expression of miR‑146a. In the present study, computational analysis was used to identify the target of miR‑146a, and a luciferase reporter assay system was used to validate NOS1 as a target gene of miR‑146a. In addition, U251 cells were treated with miR‑146a mimics/inhibitors to verify the negative regulatory association between miR‑146a and NOS1. Reverse transcription‑quantitative polymerase chain reaction analysis and western blot analysis were used to estimate the mRNA expression of NOS1 and the expression of miR‑146a. The results showed that the 'seed sequence' was located within the 3'‑untranslated region of NOS1 by searching an online miRNA database (www.mirdb.org), and the luciferase reporter assay confirmed that NOS1 was a direct target gene of miR‑146a. It was also found that the mRNA and protein expression levels of NOS1 in U251 cells treated with miR‑146a mimics and NOS1 small interfering RNA were substantially downregulated, compared with cells treated with the scramble control. The cells treated with miR‑146a inhibitors showed increased expression of NOS1. In addition, the presence of a minor allele of the rs2910164 polymorphism was significantly associated with risk of depression in patients with CAD. Taken together, the findings indicated a decreased risk of depression in the patients with CAD who were carriers of the miR‑146a rs2910164 C allele, and this association may be attributed to its ability to compromise the expression of miR‑146a, and thereby increase the expression of its target gene, NOS1.

Deciphering the genetic control of gene expression following Mycobacterium leprae antigen stimulation

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.

Each year, 200,000 new leprosy cases are reported worldwide. While there is unambiguous evidence for a role of host genetics in leprosy pathogenesis, the mechanisms by which the human host fights the infection are poorly understood. Here, we highlight the search for naturally occurring genetic variations that modulate gene expression levels following exposure to sonicate of Mycobacterium leprae, the bacterium causing the disease. Because M. leprae is not cultivable and the genuine immune cells involved in the host response during infection are still unknown, we performed a genome-wide search for such genetic variations after stimulation of whole-blood from leprosy patients with M. leprae sonicate. This design allowed to provide a general framework for the genetic control of host responses to M. leprae and outlined the contribution of host genetics to leprosy pathogenesis. Among the M. leprae-dependent genetic regulators of gene expression levels there was an enrichment of variants (i) associated with leprosy, (ii) located in transcription factor binding sites and (iii) targeted by recent positive selection.

A genome wide association study identifies a lncRna as risk factor for pathological inflammatory responses in leprosy

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.

Leprosy still affects approximately 200,000 new victims each year. A major challenge of leprosy control is the prevention of permanent disability due to nerve damage. Nerve damage occurs if leprosy remains undiagnosed for extended periods or when patients undergo pathological inflammatory responses termed Type-1 Reactions (T1R). T1R is a rare example where beneficial inflammatory responses are temporal separated from host pathological responses. There is strong experimental evidence that supports a role of host genetic factors in T1R susceptibility. Here, we employed a genome-wide association study (GWAS) to investigate susceptibility factors for T1R in Vietnamese families. We followed up the initial GWAS findings in independent population samples from Vietnam and Brazil and identified a set of cis-eQTL genetic variants for the ENSG00000235140 lncRNA as global risk factors for T1R. To test our proposal that T1R is a strong model for pathological inflammatory responses we evaluated if inflammatory bowel disease (IBD) genetic risk-factors were enriched among T1R risk factors. We observed that more than 10% of IBD-risk loci were nominally associated with risk for T1R suggesting a shared mechanism of excessive inflammatory response in the both disease etiologies.

Age-Dependent Association of TNFSF15/TNFSF8 Variants and Leprosy Type 1 Reaction

A current major challenge in leprosy control is the prevention of permanent disabilities. Host pathological inflammatory responses termed type 1 reaction (T1R) are a leading cause of nerve damage for leprosy patients. The environmental or inherited factors that predispose leprosy cases to undergo T1R are not known. However, studies have shown an important contribution of host genetics for susceptibility to T1R. We have previously identified variants encompassing the TNFSF15/TNFSF8 genes as T1R risk factors in a Vietnamese sample and replicated this association in a Brazilian sample. However, we failed to validate in Brazilian patients the strong association of TNFSF15/TNFSF8 markers rs6478108 and rs7863183 with T1R that we had observed in Vietnamese patients. Here, we investigated if the lack of validation of these variants was due to age-dependent effects on association using four independent population samples, two from Brazil and two from Vietnam. In the combined analysis across the four samples, we observed a strong association of the TNFSF15/TNFSF8 variants rs6478108, rs7863183, and rs3181348 with T1R (p_combined = 1.5E−05, p_combined = 1.8E−05, and p_combined = 6.5E−06, respectively). However, the association of rs6478108 with T1R was more pronounced in leprosy cases under 30 years of age compared to the global sample [odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.54–2.46, p_combined = 2.5E−08 versus OR = 1.46, 95% CI = 1.23–1.73, p_combined = 1.5E−05]. A multivariable analysis indicated that the association of rs6478108 with T1R was independent of either rs7863183 or rs3181348. These three variants are known regulators of the TNFSF8 gene transcription level in multiple tissues. The age dependency of association of rs6478108 and T1R suggests that the genetic control of gene expression varies across the human life span.

Soil-transmitted helminth infections and leprosy: a cross-sectional study of the association between two major neglected tropical diseases in Indonesia

Background

The clinical spectrum of leprosy is dependent on the host immune response against Mycobacterium leprae or the newly discovered Mycobacterium lepromatosis antigen. Helminth infections have been shown to affect the development of several diseases through immune regulation and thus may play a role in the clinical manifestations of leprosy and leprosy reactions. The purpose of this study is to determine the proportion of helminth infections in leprosy and its association with the type of leprosy and type 2 leprosy reaction (T2R).

Methods

History or episode of T2R was obtained and direct smear, formalin-ether sedimentation technique, and Kato-Katz smear were performed on 20 paucibacillary (PB) and 61 multibacillary (MB) leprosy participants.

Results

There are more helminth-positive participants in MB leprosy compared to PB (11/61 versus 0/20, p = 0.034) and in T2R participants compared to non-T2R (8/31 versus 3/50, p = 0.018).

Conclusions

Our results suggest that soil-transmitted helminth infections may have a role in the progression to a more severe type of leprosy, as well as the occurrence of T2R. These findings could serve as a fundamental base for clinicians to perform parasitological feces examination in patients who have MB leprosy and severe recurrent reactions to rule out the possibility of helminth infection. Further secondary confirmation of findings are needed to support these conclusions.

Pauci- and Multibacillary Leprosy: Two Distinct, Genetically Neglected Diseases

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.

Gene Association with Leprosy: A Review of Published Data

Leprosy is a chronic infectious disease caused by an obligate intracellular bacterium known as Mycobacterium leprae. Exposure to the bacillus is necessary, but this alone does not mean an individual will develop clinical symptoms of the disease. In recent years, several genes have been associated with leprosy and the innate immune response pathways converge on the main hypothesis that genes are involved in the susceptibility for the disease in two distinct steps: for leprosy per se and in the development of the different clinical forms. These genes participate in the sensing, main metabolic pathway of immune response activation and, subsequently, on the evolution of the disease into its clinical forms. The aim of this review is to highlight the role of innate immune response in the context of leprosy, stressing their participation in the signaling and targeting processes in response to bacillus infection and on the evolution to the clinical forms of the disease.

Genetics of leprosy: Expected-and unexpected-developments and perspectives

A solid body of evidence produced over decades of intense research supports the hypothesis that leprosy phenotypes are largely dependent on the genetic characteristics of the host. The early evidence of a major gene effect controlling susceptibility to leprosy came from studies of familial aggregation, twins, and complex segregation analysis. Later, linkage and association analysis, first applied to the investigation of candidate genes and chromosomal regions and more recently, to genome-wide scans, have revealed several HLA and non-HLA gene variants as risk factors for leprosy phenotypes such as disease per se, its clinical forms, and leprosy reactions. In addition, powerful, hypothesis-free strategies such as genome-wide association studies have led to an exciting, unexpected development: Leprosy susceptibility genes seem to be shared with Crohn's and Parkinson's disease. Today, a major challenge is to find the exact variants causing the biological effect underlying the genetic associations. New technologies, such as Next Generation Sequencing-that allows, for the first time, the cost- and time-effective sequencing of a complete human genome-hold the promise to reveal such variants; thus, strategies can be developed to study the functional impact of these variants in the context of infection, hopefully leading to the development of new targets for leprosy treatment and prevention.

Human Genetic Ancestral Composition Correlates with the Origin of Mycobacterium leprae Strains in a Leprosy Endemic Population

Recent reports have suggested that leprosy originated in Africa, extended to Asia and Europe, and arrived in the Americas during European colonization and the African slave trade. Due to colonization, the contemporary Colombian population is an admixture of Native-American, European and African ancestries. Because microorganisms are known to accompany humans during migrations, patterns of human migration can be traced by examining genomic changes in associated microbes. The current study analyzed 118 leprosy cases and 116 unrelated controls from two Colombian regions endemic for leprosy (Atlantic and Andean) in order to determine possible associations of leprosy with patient ancestral background (determined using 36 ancestry informative markers), Mycobacterium leprae genotype and/or patient geographical origin. We found significant differences between ancestral genetic composition. European components were predominant in Andean populations. In contrast, African components were higher in the Atlantic region. M. leprae genotypes were then analyzed for cluster associations and compared with the ancestral composition of leprosy patients. Two M. leprae principal clusters were found: haplotypes C54 and T45. Haplotype C54 associated with African origin and was more frequent in patients from the Atlantic region with a high African component. In contrast, haplotype T45 associated with European origin and was more frequent in Andean patients with a higher European component. These results suggest that the human and M. leprae genomes have co-existed since the African and European origins of the disease, with leprosy ultimately arriving in Colombia during colonization. Distinct M. leprae strains followed European and African settlement in the country and can be detected in contemporary Colombian populations.

Contemporary Colombian population is an admixture of three ancestries: Native-American, European and African. Genetic studies of human ancestry have found associations with disease, likely due to the fact that microorganisms have accompanied humans during migrations. Taking these facts into account, we studied the effect of human ancestry, Mycobacterium leprae genotype and the geographical origin of our study population, on leprosy. We found correlations between ancestral composition and M. leprae genotype: an African component is higher in the Atlantic region and a European component is higher in Andean populations (p<0.05). An interesting connection was found between the ancestral composition and two principal types of M. leprae isolates: type C54 (of African origin) was more frequent in Atlantic region populations, and type T45 (of European origin) was more frequent in the Andean region, suggesting that human and bacterial genomes have co-existed since leprosy’s origins and that leprosy has circulated with human migration.

Genetics of leprosy: expected and unexpected developments and perspectives

A solid body of evidence produced over decades of intense research supports the hypothesis that leprosy phenotypes are largely dependent on the genetic characteristics of the host. The early evidence of a major gene effect controlling susceptibility to leprosy came from studies of familial aggregation, twins, and Complex Segregation Analysis. Later, linkage and association analysis, first applied to the investigation of candidate genes and chromosomal regions and more recently, to genome-wide scans, have revealed several leukocyte antigen complex and nonleukocyte antigen complex gene variants as risk factors for leprosy phenotypes such as disease per se, its clinical forms and leprosy reactions. In addition, powerful, hypothesis-free strategies such as Genome-Wide Association Studies have led to an exciting, unexpected development: Leprosy susceptibility genes seem to be shared with Crohn's and Parkinson's diseases. Today, a major challenge is to find the exact variants causing the biological effect underlying the genetic associations. New technologies, such as Next Generation Sequencing that allows, for the first time, the cost and time-effective sequencing of a complete human genome, hold the promise to reveal such variants. Strategies can be developed to study the functional effect of these variants in the context of infection, hopefully leading to the development of new targets for leprosy treatment and prevention.

Association of TNFSF8 regulatory variants with excessive inflammatory responses but not leprosy per se

BACKGROUND

Type 1 reactions (T1R) affect a considerable proportion of patients with leprosy. In those with T1R, the host immune response pathologically overcompensates for the actual infectious threat, resulting in nerve damage and permanent disability. Based on the results of a genome-wide association study of leprosy per se, we investigated the TNFSF15 chromosomal region for a possible contribution to susceptibility to T1R.

METHODS

We performed a high-resolution association scan of the TNFSF15 locus to evaluate the association with T1R in 2 geographically and ethnically distinct populations: a family-based sample from Vietnam and a case-control sample from Brazil, comprising a total of 1768 subjects.

RESULTS

In the Vietnamese sample, 47 single-nucleotide polymorphisms (SNPs) overlapping TNFSF15 and the adjacent TNFSF8 gene were associated with T1R but not with leprosy. Of the 47 SNPs, 39 were cis-expression quantitative trait loci (cis-eQTL) for TNFSF8 including SNPs located within the TNFSF15 gene. In the Brazilian sample, 18 of these cis-eQTL SNPs overlapping the TNFSF8 gene were validated for association with T1R.

CONCLUSIONS

Taken together, these results indicate TNFSF8 and not TNFSF15 as an important T1R susceptibility gene. Our data support the need for infection genetics to go beyond genes for pathogen control to explore genes involved in a commensurate host response.

Pre-miR-146a (rs2910164 G>C) single nucleotide polymorphism is genetically and functionally associated with leprosy

Mycobacterium leprae infects macrophages and Schwann cells inducing a gene expression program to facilitate its replication and progression to disease. MicroRNAs (miRNAs) are key regulators of gene expression and could be involved during the infection. To address the genetic influence of miRNAs in leprosy, we enrolled 1,098 individuals and conducted a case-control analysis in order to study four miRNAs genes containing single nucleotide polymorphism (miRSNP). We tested miRSNP-125a (rs12975333 G>T), miRSNP-223 (rs34952329 *>T), miRSNP-196a-2 (rs11614913 C>T) and miRSNP-146a (rs2910164 G>C). Amongst them, miRSNP-146a was the unique gene associated with risk to leprosy per se (GC OR = 1.44, p = 0.04; CC OR = 2.18, p = 0.0091). We replicated this finding showing that the C-allele was over-transmitted (p = 0.003) using a transmission-disequilibrium test. A functional analysis revealed that live M. leprae (MOI 100∶1) was able to induce miR-146a expression in THP-1 (p<0.05). Furthermore, pure neural leprosy biopsies expressed augmented levels of that miRNA as compared to biopsy samples from neuropathies not related with leprosy (p = 0.001). Interestingly, carriers of the risk variant (C-allele) produce higher levels of mature miR-146a in nerves (p = 0.04). From skin biopsies, although we observed augmented levels of miR-146a, we were not able to correlate it with a particular clinical form or neither host genotype. MiR-146a is known to modulate TNF levels, thus we assessed TNF expression (nerve biopsies) and released by peripheral blood mononuclear cells infected with BCG Moreau. In both cases lower TNF levels correlates with subjects carrying the risk C-allele, (p = 0.0453 and p = 0.0352; respectively), which is consistent with an immunomodulatory role of this miRNA in leprosy.

In spite of the successful drug therapy, leprosy is still affecting people worldwide. It is well known that host genetic background influences leprosy development and that genetic variants have been associated with the disease. Therefore we conducted a study to evaluate the role of microRNAs (miRNAs) polymorphisms in leprosy. We observed that a polymorphism in miR-146a is associated with the risk to develop leprosy in Brazilians. Based on the analysis of clinical specimens, we found that the genetic variant was correlated with elevated levels of miR-146a and it is also a negative regulator of tumor necrosis factor (TNF), an important inflammatory mediator in the leprosy context. These findings provide tenable evidences that miR-146a is important in the control of gene expression during M. leprae infection and also may contribute with leprosy development by controlling TNF levels.

CUBN and NEBL common variants in the chromosome 10p13 linkage region are associated with multibacillary leprosy in Vietnam

Leprosy is caused by infection with Mycobacterium leprae and is classified clinically into paucibacillary (PB) or multibacillary (MB) subtypes based on the number of skin lesions and the bacillary index detected in skin smears. We previously identified a major PB susceptibility locus on chromosome region 10p13 in Vietnamese families by linkage analysis. In the current study, we conducted high-density association mapping of the 9.5 Mb linkage peak on chromosome region 10p13 covering 39 genes. Using leprosy per se and leprosy subtypes as phenotypes, we employed 294 nuclear families (303 leprosy cases, 63 % MB, 37 % PB) as a discovery sample and 192 nuclear families (192 cases, 55 % MB, 45 % PB) as a replication sample. Replicated significant association signals were revealed in the genes for cubilin (CUBN) and nebulette (NEBL). In the combined sample, the C allele (frequency 0.26) at CUBN SNP rs10904831 showed association [p = 1 × 10(-5); OR 0.52 (0.38-0.7)] with MB leprosy only. Likewise, allele T (frequency 0.42) at NEBL SNP rs11012461 showed association [p = 4.2 × 10(-5); OR 2.51 (1.6-4)] with MB leprosy only. These associations remained valid for the CUBN signal when taking into account the effective number of tests performed (type I error significance threshold = 2.4 × 10(-5)). We used the results of our analyses to propose a new model for the genetic control of polarization of clinical leprosy.

Mapping of PARK2 and PACRG overlapping regulatory region reveals LD structure and functional variants in association with leprosy in unrelated indian population groups

Leprosy is a chronic infectious disease caused by Mycobacterium Leprae, where the host genetic background plays an important role toward the disease pathogenesis. Various studies have identified a number of human genes in association with leprosy or its clinical forms. However, non-replication of results has hinted at the heterogeneity among associations between different population groups, which could be due to differently evolved LD structures and differential frequencies of SNPs within the studied regions of the genome. A need for systematic and saturated mapping of the associated regions with the disease is warranted to unravel the observed heterogeneity in different populations. Mapping of the PARK2 and PACRG gene regulatory region with 96 SNPs, with a resolution of 1 SNP per 1 Kb for PARK2 gene regulatory region in a North Indian population, showed an involvement of 11 SNPs in determining the susceptibility towards leprosy. The association was replicated in a geographically distinct and unrelated population from Orissa in eastern India. In vitro reporter assays revealed that the two significantly associated SNPs, located 63.8 kb upstream of PARK2 gene and represented in a single BIN of 8 SNPs, influenced the gene expression. A comparison of BINs between Indian and Vietnamese populations revealed differences in the BIN structures, explaining the heterogeneity and also the reason for non-replication of the associated genomic region in different populations.

Leprosy is a chronic granulomatous infection caused by the intracellular organism Mycobacterium leprae. The disease affects the skin and the peripheral nerves and can cause irreversible impairment of the nerve function with consequent chronic disabilities. The prevalence of leprosy has declined dramatically after the introduction of Multidrug therapy in the 1980s. However, the infection continues to survive as a major public health problem with more than 200,000 new cases reported globally every year, especially in China and India. The disease is governed by host genetic background, where several genes have been identified in association with leprosy or its clinical forms. The involvement of the PARK2 and PACRG genes with leprosy susceptibility in two distinct populations of the world, Vietnamese and Brazilian, and its non-replication in other populations suggests unravelling the reasons of heterogeneity between different population groups. The possibility of involvement of other variants and a differential LD structure for the PARK2 regulatory region in Indian populations as compared to Brazilian and Vietnamese provides an answer to the heterogeneity among associations observed previously in different population groups.

[Leprosy, a pillar of human genetics of infectious diseases]

Despite a natural reservoir of Mycobacterium leprae limited to humans and free availability of an effective antibiotic treatment, more than 200,000 people develop leprosy each year. This disease remains a major cause of disability and social stigma worldwide. The cause of this constant incidence is currently unknown and indicates that important aspects of the complex relationship between the pathogen and its human host remain to be discovered. An important contribution of host genetics to susceptibility to leprosy has long been suggested to account for the considerable variability between individuals sustainably exposed to M. leprae. Given the inability to cultivate M. leprae in vitro and in the absence of relevant animal model, genetic epidemiology is the main strategy used to identify the genes and, consequently, the immunological pathways involved in protective immunity to M. leprae. Recent genome-wide studies have identified new pathophysiological pathways which importance is only beginning to be understood. In addition, the prism of human genetics placed leprosy at the crossroads of other common diseases such as Crohn's disease, asthma or myocardial infarction. Therefore, novel lights on the pathogenesis of many common diseases could eventually emerge from the detailed understanding of a disease of the shadows.

Age-dependent association between pulmonary tuberculosis and common TOX variants in the 8q12-13 linkage region

Only a small fraction of individuals infected with Mycobacterium tuberculosis develop clinical tuberculosis (TB) in their lifetime. Genetic epidemiological evidence suggests a genetic determinism of pulmonary TB (PTB), but the molecular basis of genetic predisposition to PTB remains largely unknown. We used a positional-cloning approach to carry out ultrafine linkage-disequilibrium mapping of a previously identified susceptibility locus in chromosomal region 8q12-13 by genotyping 3,216 SNPs in a family-based Moroccan sample including 286 offspring with PTB. We observed 44 PTB-associated SNPs (p < 0.01), which were genotyped in an independent set of 317 cases and 650 controls from Morocco. A single signal, consisting of two correlated SNPs close to TOX, rs1568952 and rs2726600 (combined p = 1.1 × 10(-5) and 9.2 × 10(-5), respectively), was replicated. Stronger evidence of association was found in individuals who developed PTB before the age of 25 years (combined p for rs1568952 = 4.4 × 10(-8); odds ratio of PTB for AA versus AG/GG = 3.09 [1.99-4.78]). The association with rs2726600 (p = 0.04) was subsequently replicated in PTB-affected subjects under 25 years in a study of 243 nuclear families from Madagascar. Stronger evidence of replication in Madagascar was obtained for additional SNPs in strong linkage disequilibrium with the two initial SNPs (p = 0.003 for rs2726597), further confirming the signal. We thus identified around rs1568952 and rs2726600 a cluster of SNPs strongly associated with early-onset PTB in Morocco and Madagascar. SNP rs2726600 is located in a transcription-factor binding site in the 3' region of TOX, and further functional explorations will focus on CD4 T lymphocytes.

PARK2 mediates interleukin 6 and monocyte chemoattractant protein 1 production by human macrophages

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₃ (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 is an infectious disease with a strong host genetic component. The identification of host genetic lesions predisposing to disease is a powerful approach for mapping key junctions in the host pathogen interplay. Genetic variants located in the promoter region of the PARK2 gene are replicated leprosy susceptibility factors. To better understand a possible contribution of PARK2 to host effector mechanisms in leprosy patients, we developed a cellular model to test the contribution of the PARK2 encoded parkin protein to host responses to mycobacterial antigens. We observed that parkin was a mediator of IL-6 production in response to mycobacterial antigen in both THP-1 macrophages and human Schwann cells while human monocyte-derived macrophages needed to be pre-activated with VitD to show the same impact. Parkin also impacted on the constitutive production of MCP-1. The regulatory activity of parkin on cytokine production was found to be independent of the canonical TLR-NFκB signalling pathway. We also tested association of IL6 and CCL2 gene expression levels in whole blood assays with PARK2 polymorphisms. For both cytokines, we found significant associations with those PARK2 variants that were established leprosy susceptibility factors. Hence, our results show that genetic PARK2 variants that are correlated with leprosy susceptibility are also correlated with production of these cytokines following stimulation with M. leprae sonicate.