Missense Variants in HIF1A and LACC1 Contribute to Leprosy Risk in Han Chinese

Genome-wide association studies (GWASs) and genome-wide linkage studies (GWLSs) have identified numerous risk genes affecting the susceptibility to leprosy. However, most of the reported GWAS hits are noncoding variants and account for only part of the estimated heritability for this disease. In order to identify additional risk genes and map the potentially functional variants within the GWAS loci, we performed a three-stage study combining whole-exome sequencing (WES; discovery stage), targeted next-generation sequencing (NGS; screening stage), and refined validation of risk missense variants in 1,433 individuals with leprosy and 1,625 healthy control individuals from Yunnan Province, Southwest China. We identified and validated a rare damaging variant, rs142179458 (c.1045G>A [p.Asp349Asn]) in HIF1A, as contributing to leprosy risk (p = 4.95 × 10^-9, odds ratio [OR] = 2.266). We were able to show that affected individuals harboring the risk allele presented with multibacillary leprosy at an earlier age (p = 0.025). We also confirmed the association between missense variant rs3764147 (c.760A>G [p.Ile254Val]) in the GWAS hit LACC1 (formerly C13orf31) and leprosy (p = 6.11 × 10^-18, OR = 1.605). By using the population attributable fraction, we have shown that HIF1A and LACC1 are the major genes with missense variants contributing to leprosy risk in our study groups. Consistently, mRNA expression levels of both HIF1A and LACC1 were upregulated in the skin lesions of individuals with leprosy and in Mycobacterium leprae-stimulated cells, indicating an active role of HIF1A and LACC1 in leprosy pathogenesis.