Genome-wide association study of leprosy in Malawi and Mali.

Leprosy is a chronic infection of the skin and peripheral nerves caused by Mycobacterium leprae. Despite recent improvements in disease control, leprosy remains an important cause of infectious disability globally. Large-scale genetic association studies in Chinese, Vietnamese and Indian populations have identified over 30 susceptibility loci for leprosy. There is a significant burden of leprosy in Africa, however it is uncertain whether the findings of published genetic association studies are generalizable to African populations. To address this, we conducted a genome-wide association study (GWAS) of leprosy in Malawian (327 cases, 436 controls) and Malian (247 cases, 368 controls) individuals. In that analysis, we replicated four risk loci previously reported in China, Vietnam and India; MHC Class I and II, LACC1 and SLC29A3. We further identified a novel leprosy susceptibility locus at 10q24 (rs2015583; combined p = 8.81 × 10-9; OR = 0.51 [95% CI 0.40 - 0.64]). Using publicly-available data we characterise regulatory activity at this locus, identifying ACTR1A as a candidate mediator of leprosy risk. This locus shows evidence of recent positive selection and demonstrates pleiotropy with established risk loci for inflammatory bowel disease and childhood-onset asthma. A shared genetic architecture for leprosy and inflammatory bowel disease has been previously described. We expand on this, strengthening the hypothesis that selection pressure driven by leprosy has shaped the evolution of autoimmune and atopic disease in modern populations. More broadly, our data highlights the importance of defining the genetic architecture of disease across genetically diverse populations, and that disease insights derived from GWAS in one population may not translate to all affected populations.

The Awesome Power of Human Genetics of Infectious Disease.

Since the identification of sickle cell trait as a heritable form of resistance to malaria, candidate gene studies, linkage analysis paired with sequencing, and genome-wide association (GWA) studies have revealed many examples of genetic resistance and susceptibility to infectious diseases. GWA studies enabled the identification of many common variants associated with small shifts in susceptibility to infectious diseases. This is exemplified by multiple loci associated with leprosy, malaria, HIV, tuberculosis, and coronavirus disease 2019 (COVID-19), which illuminate genetic architecture and implicate pathways underlying pathophysiology. Despite these successes, most of the heritability of infectious diseases remains to be explained. As the field advances, current limitations may be overcome by applying methodological innovations such as cellular GWA studies and phenome-wide association (PheWA) studies as well as by improving methodological rigor with more precise case definitions, deeper phenotyping, increased cohort diversity, and functional validation of candidate loci in the laboratory or human challenge studies.

Identification of shared loci associated with both Crohn's disease and leprosy in East Asians.

Genome-wide association studies (GWAS) of Crohn's disease (CD) in European and leprosy in Chinese population have shown that CD and leprosy share genetic risk loci. As these shared loci were identified through cross-comparisons across different ethnic populations, we hypothesized that meta-analysis of GWAS on CD and leprosy in East Asian populations would increase power to identify additional shared loci. We performed a cross-disease meta-analysis of GWAS data from CD (1621 cases and 4419 controls) and leprosy (2901 cases 3801 controls) followed by replication in additional datasets comprising 738 CD cases and 488 controls and 842 leprosy cases and 925 controls. We identified one novel locus at 7p22.3, rs77992257 in intron 2 of ADAP1, shared between CD and leprosy with genome-wide significance (P = 3.80 × 10-11) and confirmed 10 previously established loci in both diseases: IL23R, IL18RAP, IL12B, RIPK2, TNFSF15, ZNF365-EGR2, CCDC88B, LACC1, IL27, NOD2. Phenotype variance explained by the polygenic risk scores derived from Chinese leprosy data explained up to 5.28% of variance of Korean CD, supporting similar genetic structures between the two diseases. Although CD and leprosy shared a substantial number of genetic susceptibility loci in East Asians, the majority of shared susceptibility loci showed allelic effects in the opposite direction. Investigation of the genetic correlation using cross-trait linkage disequilibrium score regression also showed a negative genetic correlation between CD and leprosy (rg [SE] = -0.40[0.13], P = 2.6 × 10-3). These observations implicate the possibility that CD might be caused by hyper-sensitive reactions toward pathogenic stimuli.

Deep resequencing identifies candidate functional genes in leprosy GWAS loci.

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.

Predictive nomogram for leprosy using genetic and epidemiological risk factors in Southwestern China: Case-control and prospective analyses.

BACKGROUND: There is a high incidence of leprosy among house-contacts compared with the general population. We aimed to establish a predictive model using these genetic factors along with epidemiological factors to predict leprosy risk of leprosy household contacts (HHCs). METHODS: Weighted genetic risk score (wGRS) encompassing genome wide association studies (GWAS) variants and five non-genetic factors were examined in a case-control design associated with leprosy risk including 589 cases and 647 controls from leprosy HHCs. We constructed a risk prediction nomogram and evaluated its performance by concordance index (C-index) and calibration curve. The results were validated using bootstrap resampling with 1000 resamples and a prospective design including 1100 HHCs of leprosy patients. FINDING: The C-index for the risk model was 0·792 (95% confidence interval [CI] 0·768-0·817), and was confirmed to be 0·780 through bootstrapping validation. The calibration curve for the probability of leprosy showed good agreement between the prediction of the nomogram and actual observation. HHCs were then divided into the low-risk group (nomogram score ≤ 81) and the high-risk group (nomogram score > 81). In prospective analysis, 12 of 1100 participants had leprosy during 63 months' follow-up. We generated the nomogram for leprosy in the validation cohort (C-index 0·773 [95%CI 0·658-0·888], sensitivity75·0%, specificity 66·8%). Interpretation The nomogram achieved an effective prediction of leprosy in HHCs. Using the model, the risk of an individual contact developing leprosy can be determined, which can lead to a rational preventive choice for tracing higher-risk leprosy contacts. FUNDING: The ministry of health of China, ministry of science and technology of China, Chinese academy of medical sciences, Jiangsu provincial department of science and technology, Nanjing municipal science and technology bureau.

Family-based genome-wide association study of leprosy in Vietnam.

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.

Human genetics of Buruli ulcer.

Buruli ulcer, the third most common mycobacterial disease worldwide, is caused by Mycobacterium ulcerans and characterized by devastating necrotizing skin lesions. Susceptibility to Buruli ulcer is thought to depend on host genetics, but very few genetic studies have been performed. The identification of a microdeletion on chromosome 8 in a familial form of severe Buruli ulcer suggested a monogenic basis of susceptibility. The role of common host genetic variants in Buruli ulcer development has been investigated in only three candidate-gene studies targeting genes involved in mycobacterial diseases. A recent genome-wide association study suggested a probable role for long non-coding RNAs and strengthened the contribution of autophagy as a major defense mechanism against mycobacteria. In this review, we summarize the history, epidemiological and clinical aspects of Buruli ulcer, focusing particularly on genetic findings relating to susceptibility to this disease. Finally, we discuss exciting new genetic avenues arising, in particular, from studies of mouse models, and the need for different disciplines to work together, to benefit from the extensive work on other mycobacterial diseases, mostly tuberculosis and leprosy. We are convinced that such pooling of effort will lead to the development of efficient novel strategies for combatting Buruli ulcer.

A review of HLA allele and SNP associations with highly prevalent infectious diseases in human populations.

Human leucocyte antigen (HLA) alleles and single nucleotide polymorphisms (SNPs) lying in the HLA region are known to be associated with several infectious diseases among which acquired immunodeficiency syndrome, hepatitis B, hepatitis C, tuberculosis, leprosy and malaria are highly prevalent in many human populations worldwide. Distinct approaches such as case-control comparisons, immunogenetic analyses, bioinformatic peptide-binding predictions, ancient DNA and genome-wide association studies (GWAS) have contributed to improving this knowledge during the last decade, although many results still need stronger statistical and/or functional support. The present review updates the information regarding the main HLA allele and SNP associations observed to date for six of the most widespread and some other infectious diseases, and provides a synthetic illustration of these findings on a schematic HLA genomic map. It then discusses these results by stressing the importance of integrating information on HLA population diversity in disease-association studies.

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.

Discovery of 4 exonic and 1 intergenic novel susceptibility loci for leprosy.

Seven new risk coding variants have been identified through an exome-wide association study (EWAS), which studied the contributions of protein-coding variants to leprosy susceptibility. But some potential susceptibility loci were not studied in the previous EWAS study because of the project consideration. Seventeen unstudied potential susceptibility loci of the previous EWAS were validated in 3169 cases and 9814 controls in this study. Four disease-associated exonic loci were identified: rs671 in ALDH2 (P = 2.0 × 10-20 , odds ratio [OR] = 1.35), rs13259978 in SLC7A2 (P = 1.74 × 10-8 , OR = 1.28), rs925368 in GIT2 (P = 9.18 × 10-17 , OR = 1.44), and rs75680863 in TCN2 (P = 8.37 × 10-21 , OR = 0.74). Potentially implicating ZFP36L1 as a new susceptibility gene, 1 intergenic single nucleotide polymorphism (SNP), rs1465788 (P = 7.81 × 10-6 , OR = 0.88), was also suggested to be associated with leprosy. A luciferase reporter assay showed that the rs1465788 risk allele notably decreased the transcription activity of the flanking sequence. These findings suggest the possible involvement of lipid metabolism, NF-κB homeostasis and macrophage antimicrobial pathways in leprosy pathogenesis.

Genetic and molecular aspects of androgenetic alopecia.

Androgenetic alopecia is the most common form of progressive hair loss in humans. A genetic predisposition and hormonal status are considered as major risk factors for this condition. Several recent advances in molecular biology and genetics have increased our understanding of the mechanisms of hair loss in androgenetic alopecia. We review these advances and examine the trends in the genetic and molecular aspects of androgenetic alopecia.

A pathway-based association analysis identified FMNL1-MAP3K14 as susceptibility genes for leprosy.

The nuclear transcription factor-κB (NF-κB) plays a pivotal role in controlling both innate and adaptive immunity and regulates the expressions of many immunological mediators. Abundant evidences have showed the importance of NF-κB pathway in the host immune responses against Mycobacterium leprae in the development of leprosy. However, no particular association study between leprosy and NF-κB pathway-related gene polymorphisms was reported. Here, we performed a large-scale and two-stage candidate association study to investigate the association between 94 NF-κB pathway-related genes and leprosy. Our results showed that rs58744688 was significantly associated with leprosy (P = 7.57 × 10-7 , OR = 1.12) by combining the previous genomewide association data sets and four independent validation sample series, consisting of a total of 4631 leprosy cases and 6413 healthy controls. This founding implicated that MAP3K14 and FMNL1 were susceptibility genes for leprosy, which suggested the involvement of macrophage targeting and NF-κB pathway in the development of leprosy.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behçet's disease susceptibility.

We analyzed 1,900 Turkish Behçet's disease cases and 1,779 controls genotyped with the Immunochip. The most significantly associated SNP was rs1050502, a tag SNP for HLA-B*51. In the Turkish discovery set, we identified three new risk loci, IL1A-IL1B, IRF8, and CEBPB-PTPN1, with genome-wide significance (P < 5 × 10-8) by direct genotyping and ADO-EGR2 by imputation. We replicated the ADO-EGR2, IRF8, and CEBPB-PTPN1 loci by genotyping 969 Iranian cases and 826 controls. Imputed data in 608 Japanese cases and 737 controls further replicated ADO-EGR2 and IRF8, and meta-analysis additionally identified RIPK2 and LACC1. The disease-associated allele of rs4402765, the lead marker at IL1A-IL1B, was associated with both decreased IL-1α and increased IL-1ß production. ABO non-secretor genotypes for two ancestry-specific FUT2 SNPs showed strong disease association (P = 5.89 × 10-15). Our findings extend the list of susceptibility genes shared with Crohn's disease and leprosy and implicate mucosal factors and the innate immune response to microbial exposure in Behçet's disease susceptibility.

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.

Leucine-rich repeat kinase 2 inhibitors: a patent review (2014-2016).

INTRODUCTION: Leucine-rich repeat kinase 2 (LRRK2) is a member of the Tyrosine Kinase-Like (TKL) branch of the kinome tree and is a multi-domain protein that includes GTPase and kinase activity. While genome-wide association studies (GWAS) has linked LRRK2 with Crohn's disease and leprosy, it has received the greatest attention due to it being implicated as one of the genetic loci associated with autosomal dominant inheritance in Parkinson's disease (PD). Areas covered: In this review, the small molecule patent literature from 2014-2016 with a focus on composition of matter and use patents was surveyed. Scifinder was primarily searched using 'LRRK2' as the query to identify all relevant literature and then triaged for small molecule patents. Expert opinion: The patent landscape around LRRK2 continues to develop. The early patents covered using existing kinase inhibitors for use against LRRK2. This evolved to compounds specifically designed for selectivity against LRRK2, but key exemplified compounds lacked sufficient brain exposure to affect sufficient efficacy. More recent compounds have addressed this deficiency and show greater potential for treating PD. While potency will be necessary to generate medicines with low human daily doses, brain penetration and safety will be the key differentiators for ultimately determining the most effective LRRK2 disease-modifying treatment for PD.

A genome wide association study identifies a IncRna 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.

A large-scale genome-wide association and meta-analysis identified four novel susceptibility loci for leprosy.

Leprosy, a chronic infectious disease, results from the uncultivable pathogen Mycobacterium leprae (M. leprae), and usually progresses to peripheral neuropathy and permanent progressive deformity if not treated. Previously published genetic studies have identified 18 gene/loci significantly associated with leprosy at the genome-wide significant level. However as a complex disease, only a small proportion of leprosy risk could be explained by those gene/loci. To further identify more susceptibility gene/loci, we hereby performed a three-stage GWAS comprising 8,156 leprosy patients and 15,610 controls of Chinese ancestry. Four novel loci were identified including rs6807915 on 3p25.2 (P=1.94 × 10-8, OR=0.89), rs4720118 on 7p14.3 (P=3.85 × 10-10, OR=1.16), rs55894533 on 8p23.1 (P=5.07 × 10-11, OR=1.15) and rs10100465 on 8q24.11 (P=2.85 × 10-11, OR=0.85). Altogether, these findings have provided new insight and significantly expanded our understanding of the genetic basis of leprosy.

Fine mapping of the GWAS loci identifies SLC35D1 and IL23R as potential risk genes for leprosy.

BACKGROUND: Previous genome-wide association study (GWAS) identified two new leprosy associated loci (1p31.3 [rs3762318] and 6q24.3 [rs2275606]). However, there were insufficient validations in independent populations. OBJECTIVE: To validate the association and to map the potentially causal variants/genes underlying the association between the confirmed GWAS hit and leprosy. METHODS: We genotyped 10 variants in the regions encompassing the two loci in 1110 Han Chinese subjects with and without leprosy, followed by expression quantitative trait loci (eQTL), mRNA expression profiling, and network analysis. We further sequenced the exon region of four genes that were located in the confirmed GWAS hit region in 80 leprosy patients and 99 individuals without leprosy. RESULTS: We validated the positive association of rs3762318 with multibacillary leprosy (P=7.5×10-4), whereas the association of rs2275606 could not be validated. eQTL analysis showed that both the GWAS locus rs3762318 and one surrounding positively associated SNP rs2144658 (P=1.8×10-3) significantly affected the mRNA expression of a nearby gene SLC35D1, which might be involved in metabolism. Moreover, SLC35D1 was differentially expressed in skin tissues of leprosy patients, and the differential expression pattern was consistent among leprosy subtypes. Rare damaging missense variants in IL23R were significantly enriched in leprosy patients. CONCLUSION: Our results supported the positive association between the GWAS reported rs3762318 and leprosy, and SLC35D1 and IL23R might be the causal genes.

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.