Search

One of the persistent challenges of genetic association studies is the replication of genetic marker-disease associations across ethnic groups. Here, we conducted high-density association mapping of PARK2/PACRG SNPs with leprosy and identified 69 SNPs significantly associated with leprosy in 198 single-case Vietnamese leprosy families. A total of 56 associated SNPs localized to the overlapping promoter regions of PARK2/PACRG. For this region, multivariate analysis identified four SNPs belonging to two major SNP bins (rs1333955, rs7744433) and two single SNP bins (rs2023004, rs6936895) that capture the combined statistical evidence (P = 1.1 × 10(-5)) for association among Vietnamese patients. Next, we enrolled a case-control sample of 364 leprosy cases and 370 controls from Northern India. We genotyped all subjects for 149 SNPs that capture >80 % of the genetic variation in the Vietnamese sample and found 24 SNPs significantly associated with leprosy. Multivariate analysis identified three SNPs (rs1333955, rs9356058 and rs2023004) that capture the association with leprosy (P < 10(-8)). Hence, two SNPs (rs1333955 and rs2023004) were replicated by multivariate analysis between both ethnic groups. Marked differences in the linkage disequilibrium pattern explained some of the differences in univariate analysis between the two ethnic groups. In addition, the strength of association for two promoter region SNP bins was significantly stronger among young leprosy patients in the Vietnamese sample. The same trend was observed in the Indian sample, but due to the higher age-at-diagnosis of the patients the age effect was less pronounced.

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.

Host genetics has an important role in leprosy, and variants in the shared promoter region of PARK2 and PACRG were the first major susceptibility factors identified by positional cloning. Here we report the linkage disequilibrium mapping of the second linkage peak of our previous genome-wide scan, located close to the HLA complex. In both a Vietnamese familial sample and an Indian case-control sample, the low-producing lymphotoxin-alpha (LTA)+80 A allele was significantly associated with an increase in leprosy risk (P = 0.007 and P = 0.01, respectively). Analysis of an additional case-control sample from Brazil and an additional familial sample from Vietnam showed that the LTA+80 effect was much stronger in young individuals. In the combined sample of 298 Vietnamese familial trios, the odds ratio of leprosy for LTA+80 AA/AC versus CC subjects was 2.11 (P = 0.000024), which increased to 5.63 (P = 0.0000004) in the subsample of 121 trios of affected individuals diagnosed before 16 years of age. In addition to identifying LTA as a major gene associated with early-onset leprosy, our study highlights the critical role of case- and population-specific factors in the dissection of susceptibility variants in complex diseases.