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.

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.

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.

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.

A Missense LRRK2 Variant Is a Risk Factor for Excessive Inflammatory Responses in Leprosy

Background

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.

Methodology

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.

Principal Findings

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.

Significance

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.

A major challenge of current leprosy control is the management of host pathological immune responses coined Type-1 Reactions (T1R). T1R are characterized by acute inflammatory episodes whereby cellular immune responses are directed against host peripheral nerve cells. T1R affects up half of all leprosy patients and are a major cause of leprosy-associated disabilities. Since there is evidence that host genetic factors predispose leprosy patients to T1R, we have conducted a candidate gene study to test if LRRK2 gene variants are T1R risk factors. The choice of LRRK2 was motivated by the fact that LRRK2 was associated with leprosy per se in some but not in other studies. We reasoned that this may reflect different proportions of leprosy patients with T1R in the different samples and that LRRK2 may in truth be a T1R susceptibility gene. Here, we show that variants overlapping the LRRK2 gene, reported as suggestive leprosy per se susceptibility factors in a previous genome-wide association study, are preferentially associated with T1R. The main SNP carrying most of the association signal is the amino-acid change M2397T (rs3761863) which is known to impact LRRK2 turnover. Interestingly, eQTL SNPs counterbalanced the effect of the M2397T variant but this compensatory mechanism was abrogated by Mycobacterium leprae antigen stimulation.

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.

Gene set signature of reversal reaction type I in leprosy patients

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.

Leprosy type 1 reversal reactions (T1R) are an important cause of nerve damage in leprosy patients and accurate prediction of patients at increased risk of T1R is a major challenge of current leprosy control. The incidence of T1R differs widely from 6% to 67% of leprosy patients in different leprosy endemic settings. Whether or not this reflects the impact of unknown environmental triggers or differences in the genetic background across ethnicities is not known. We performed a comparative transcriptome analysis between leprosy patients affected and free of T1R in response to M. leprae antigens. As the discovery sample we enrolled cured leprosy patients who had been diagnosed with T1R at the time of leprosy diagnosis and leprosy patients who had never undergone T1R (retrospective arm). Whole genome transcriptome analysis after stimulation of blood with M. leprae antigen resulted in the definition of a T1R signature gene set. We validated the T1R gene set in RNA samples obtained from T1R-free patients at the time of leprosy diagnosis and followed for 3 years for development of T1R (prospective arm). These results confirm the role of innate factors in T1R and are a first step towards a predictive genetic T1R signature.

Recommendations for a patient-centered approach to the assessment and treatment of scalp psoriasis: a consensus statement from the Asia Scalp Psoriasis Study Group

BACKGROUND

International consensus statements on the management of scalp psoriasis are available, but no such recommendations exist for Asia.

METHODS

The Asia Scalp Psoriasis Study Group (ASPSG) met in May 2011 to review the epidemiologic pattern of scalp psoriasis in Southeast Asia and to develop Asia-specific recommendations for its management.

RESULTS

The overall prevalence of psoriasis in Asia is <0.3%, but 75-90% have scalp involvement, whether isolated or with lesions elsewhere, which can negatively impact quality of life (QoL). Treatment decisions should be based primarily on objective disease severity, but should also take account of patient QoL. Psychosocial support and more aggressive treatment should be offered to all patients with moderate to severe QoL impairment. Topical therapy is indicated first-line in all patients, with combination therapy (corticosteroid + calcipotriol), more occlusive formulations, keratolytics, and very potent corticosteroids for patients needing greater or faster efficacy. Systemic therapies, light or laser treatments should be reserved for patients with severe and recalcitrant disease.

CONCLUSIONS

The ASPSG recommends a patient-centered approach to scalp psoriasis management, consistent with the international consensus statements. Asian physicians should also consider patient QoL, prior treatment response, formulation preferences, likely adherence, cost, time available for self-management, and potential adverse events.

Crohn's disease susceptibility genes are associated with leprosy in the Vietnamese population

A genomewide association study in Chinese patients with leprosy detected association signals in 16 single-nucleotide polymorphisms (SNPs) belonging to 6 loci, of which 4 are related to the NOD2 signaling pathway and are Crohn's disease susceptibility loci. Here, we studied these 16 SNPs as potential leprosy susceptibility factors in 474 Vietnamese leprosy simplex families. We replicated SNPs at HLA-DR-DQ, RIPK2, CCDC122-LACC1, and NOD2 as leprosy susceptibility factors in Vietnam. These results validated the striking overlap in the genetic control of Crohn's disease and leprosy.

Human leukocyte antigen class I region single-nucleotide polymorphisms are associated with leprosy susceptibility in Vietnam and India

Experimental evidence suggested the existence of unidentified leprosy susceptibility loci in the human leukocyte antigen (HLA) complex. To identify such genetic risk factors, a high-density association scan of a 1.9-mega-base (Mb) region in the HLA complex was performed. Among 682 single-nucleotide polymorphisms (SNPs), 59 were associated with leprosy (P <.01) in 198 Vietnamese single-case leprosy families. Genotyping of these SNPs in an independent sample of 292 Vietnamese single-case leprosy families replicated the association of 12 SNPs (P <.01). Multivariate analysis of these 12 SNPs showed that the association information could be captured by 2 intergenic HLA class I region SNPs (P = 9.4 × 10⁻⁹)-rs2394885 and rs2922997 (marginal multivariate P = 2.1 × 10⁻⁷ and P = .0016, respectively). SNP rs2394885 tagged the HLA-C*15:05 allele in the Vietnamese population. The identical associations were validated in a third sample of 364 patients with leprosy and 371 control subjects from North India. These results implicated class I alleles in leprosy pathogenesis.

Genetic and functional analysis of common MRC1 exon 7 polymorphisms in leprosy susceptibility

The chromosomal region 10p13 has been linked to paucibacillary leprosy in two independent studies. The MRC1 gene, encoding the human mannose receptor (MR), is located in the 10p13 region and non-synonymous SNPs in exon 7 of the gene have been suggested as leprosy susceptibility factors. We determined that G396S is the only non-synonymous exon 7-encoded polymorphism in 396 unrelated Vietnamese subjects. This SNP was genotyped in 490 simplex and 90 multiplex leprosy families comprising 704 patients (47% paucibacillary; 53% multibacillary). We observed significant under-transmission of the serine allele of the G396S polymorphism with leprosy per se (P = 0.036) and multibacillary leprosy (P = 0.034). In a sample of 384 Brazilian leprosy cases (51% paucibacillary; 49% multibacillary) and 399 healthy controls, we observed significant association of the glycine allele of the G396S polymorphism with leprosy per se (P = 0.016) and multibacillary leprosy (P = 0.023). In addition, we observed a significant association of exon 7 encoded amino acid haplotypes with leprosy per se (P = 0.012) and multibacillary leprosy (P = 0.004). Next, we tested HEK293 cells over-expressing MR constructs (293-MR) with three exon 7 haplotypes of MRC1 for their ability to bind and internalize ovalbumin and zymosan, two classical MR ligands. No difference in uptake was measured between the variants. In addition, 293-MR failed to bind and internalize viable Mycobacterium leprae and BCG. We propose that the MR-M. leprae interaction is modulated by an accessory host molecule of unknown identity.

Genomewide linkage analysis of the granulomatous mitsuda reaction implicates chromosomal regions 2q35 and 17q21

The Mitsuda reaction, a delayed granulomatous skin reaction elicited by the intradermal injection of heat-killed Mycobacterium leprae, is an in vivo test reflecting the ability to generate an immune granuloma after sensitization by diverse mycobacterial infections. Accumulating evidence for the genetic control of the Mitsuda reaction has been reported. We performed a genomewide linkage scan for the quantitative Mitsuda reaction in 19 large families from Vietnam with a history of leprosy (114 offspring). Suggestive linkage was found at chromosomal regions 2q35 (P = 9 x 10(-4) at the SLC11A1 locus) and 17q21-25 (P = 8 x 10(-4)). Interestingly, these 2 regions have been previously linked to mycobacterial infection and other granulomatous diseases.

HLA-DRB1*04 and DRB1*10 are associated with resistance and susceptibility, respectively, in Brazilian and Vietnamese leprosy patients

The host genetic background has been considered one of the factors that influence leprosy outcome, a chronic infectious disease caused by Mycobacterium leprae. Genome scans demonstrated that the 6p21 region is associated with leprosy and a substantial number of population-based studies analyzing human leukocyte antigen (HLA) class II loci suggested association of HLA-DR with leprosy. However, some studies lacked robustness as they had limited power. Indeed, experimental designs require increased sample size to achieve adequate power, as well as replication studies with independent samples for confirmation of previous findings. In this work, we analyzed the influence of the HLA-DRB1 locus on leprosy susceptibility per se and disease type using a case-control design carried out in Brazilians (578 cases and 691 controls) and a replication study based on a family design in a Vietnamese population (n=194 families). The results showed that HLA-DRB1*10 is associated with susceptibility to leprosy and HLA-DRB1*04 is associated with resistance, both in the Brazilian and Vietnamese populations suggesting that these alleles play an important role in the activation of cellular immune responses against M. leprae.

A Recessive Major Gene Controls the Mitsuda Reaction in a Region Endemic for Leprosy

BACKGROUND

Leprosy is a chronic infectious disease caused by Mycobacterium leprae. The Mitsuda reaction is a delayed granulomatous skin reaction elicited by intradermal injection of heat-killed M. leprae. Interestingly, results of the Mitsuda test are positive in the majority of individuals, even in areas not endemic for M. leprae. Like leprosy, the Mitsuda reaction is thought to be genetically controlled, but its mode of inheritance is unknown, although the role of the NRAMP1 gene has previously been reported.

METHODS

We conducted a segregation analysis of quantitative Mitsuda reactivity in 168 Vietnamese nuclear families ascertained through patients with leprosy.

RESULTS

We found strong evidence (P<10-9) for a major gene controlling the Mitsuda reaction independently of leprosy clinical status. Subsequent linkage analysis showed that this major gene was distinct from NRAMP1. Under the major-gene model, approximately 12% of individuals are homozygous for the recessive predisposing allele and are predicted to display high levels of Mitsuda reactivity (mean, approximately 10 mm, versus 5 mm in other individuals).

CONCLUSION

We provide evidence that the Mitsuda reaction is controlled by a major gene. Our study paves the way for the identification of this gene and should provide novel insight into the mechanisms involved in granuloma formation, especially in M. leprae infection.

Chromosome 6q25 is linked to susceptibility to leprosy in a Vietnamese population

Leprosy, a chronic infectious disease caused by Mycobacterium leprae, affects an estimated 700,000 persons each year. Clinically, leprosy can be categorized as paucibacillary or multibacillary disease. These clinical forms develop in persons that are intrinsically susceptible to leprosy per se, that is, leprosy independent of its specific clinical manifestation. We report here on a genome-wide search for loci controlling susceptibility to leprosy per se in a panel of 86 families including 205 siblings affected with leprosy from Southern Vietnam. Using model-free linkage analysis, we found significant evidence for a susceptibility gene on chromosome region 6q25 (maximum likelihood binomial (MLB) lod score 4.31; P = 5 x 10(-6)). We confirmed this by family-based association analysis in an independent panel of 208 Vietnamese leprosy simplex families. Of seven microsatellite markers underlying the linkage peak, alleles of two markers (D6S1035 and D6S305) showed strong evidence for association with leprosy (P = 6.7 x 10(-4) and P = 5.9 x 10(-5), respectively).