Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 768259 10.3389/fgene.2021.768259 Genetics Original Research Fine-Mapping of the Major Histocompatibility Complex Region Linked to Leprosy in Northern China Zhang et al. HLA Imputation Related to Leprosy Zhang Ruixue 1 2 3 4 5 † Cao Lu 1 2 3 4 5 † Chen Weiwei 1 2 3 4 5 Ge Huiyao 1 2 3 4 5 Hu Xia 1 2 3 4 5 Li Zhuo 1 2 3 4 5 Wang Yirui 1 2 3 4 5 Fan Wencheng 1 2 3 4 5 Yong Liang 1 2 3 4 5 Yu Yafen 1 2 3 4 5 Mao Yiwen 1 2 3 4 5 Zhen Qi 1 2 3 4 5 Liu Hong 6 Zhang Furen 6 Sun Liangdan 1 2 3 4 5 * Department of Dermatology, No. 1 Hospital, Anhui Medical University, Hefei, China Institute of Dermatology, Anhui Medical University, Hefei, China Key Laboratory of Dermatology, Anhui Medical University, Ministry of Education, Hefei, China Anhui Provincial Institute of Translational Medicine, Hefei, China Inflammation and Immune Mediated Diseases Laboratory of Anhui Province, Hefei, China Shandong Provincial Hospital for Skin Diseases and Shandong Provincial Institute of Dermatology and Venereology, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China

Edited by: Deng-Feng Zhang, Kunming Institute of Zoology (CAS), China

Reviewed by: Erick C. Castelli, São Paulo State University, Brazil

Craig Kinnear, Stellenbosch University, South Africa

 *Correspondence: Liangdan Sun, ahmusld@163.com

These authors have contributed equally to this work

This article was submitted to Applied Genetic Epidemiology, a section of the journal Frontiers in Genetics

 16 12 2021 2021 12 768259 15 09 2021 29 11 2021 Copyright © 2021 Zhang, Cao, Chen, Ge, Hu, Li, Wang, Fan, Yong, Yu, Mao, Zhen, Liu, Zhang and Sun. 2021 Zhang, Cao, Chen, Ge, Hu, Li, Wang, Fan, Yong, Yu, Mao, Zhen, Liu, Zhang and Sun

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Background: Leprosy is a chronic infectious skin and neurological disease, and genetic background is considered to be one of the major factors of risk. The major histocompatibility complex (MHC) region not only affects susceptibility to leprosy but also its development and outcome. Given the complex traits of the MHC region, variants and the potential mechanism by which HLA influences leprosy development need to be further explored.

Methods: We extracted previous genome-wide association study data from the Northern Han Chinese population to perform HLA fine-mapping. Using the 1,000 Genome Project Phase 3 dataset as the reference panel, single-nucleotide polymorphisms (SNP), insertion and deletion (INDEL) and copy number variant (CNV) imputation were carried out. HLA classical alleles and amino acids in the MHC region were imputed using the HAN-MHC database. Further stepwise regression analysis was conducted to analyze independent signals of variants related to leprosy.

Results: We identified four independent variants: esv3608598, rs7754498, rs3130781 and rs144388449. Among them, esv3608598 is a CNV and the first HLA CNV associated with leprosy risk. SNP annotation using RegulomeDB, HaploReg, and rVarBase showed that three SNPs are likely to affect the pathogenesis of leprosy.

Conclusion: In summary, this is the first study to assess the association between HLA CNV and leprosy susceptibility in a Northern Han Chinese population. By fine mapping of the MHC region in this population, our findings provide evidence for the contribution of HLA to leprosy susceptibility.

 leprosy MHC imputation single nucleotide polymorphism copy number variant Introduction

Leprosy, or Hansen’s disease, is a chronic infectious skin and neurological disease caused by Mycobacterium leprae. In addition to the skin and peripheral nervous system, it can affect the eyes, bones, mucous membranes, testes and other body parts to a certain extent (White and Franco-Paredes, 2015). Although more than six million people have been cured by multidrug therapy since the 1980s, leprosy is still an important public health problem in many regions of the world (Moraes and Düppre, 2021). Indeed, the World Health Organization (WHO) reported 202,185 new cases of leprosy in 2019 (WHO, 2020).

Ridley and Jopling proposed in 1966 that leprosy can be divided into five main clinical forms according to the immune response, histopathology and bacterial factors (4). Tuberculoid leprosy, at one end of the clinical spectrum, involves an immune response mediated by Th1 cells, which is characterized by limited skin lesions and a small number of bacteria. At the other end is lepromatous leprosy, with weak immunity, that is characterized by a large number of lesions and intense growth of the bacillus. Between the two ends of the spectrum are borderline tuberculoid leprosy, mid-borderline leprosy and borderline lepromatous leprosy, immune-instability forms. In the early stage, indeterminate leprosy can evolve into any of the above forms or heal spontaneously (Ridley and Jopling, 1966). To standardize clinical treatment, the WHO divides leprosy into multibacillary and paucibacillary types (WHO, 1994), (Parkash, 2009).

M. leprae is an intracellular acid-fast pathogen that cannot grow in culture medium and is therefore strongly dependent on the host cell environment. Host-pathogen interactions in leprosy over thousands of years have caused M. leprae to lose part of its genome. However, the bacterium can still invade human macrophages and Schwann cells (Misch et al., 2010). In fact, only a small percentage of people infected will develop leprosy (Moet et al., 2006). Environmental factors, pathogen burden, and human genetic susceptibility may be responsible for clinical leprosy caused by exposure to M. leprae. The extremely low rate of leprosy suggests that the genetic background of the host is an important aspect of susceptibility (Cole et al., 2001).

With the development and popularization of molecular biotechnology, genome-wide association studies (GWASs) and further fine-mapping analysis in specific regions have become important methods to study leprosy susceptibility genes (Fava et al., 2020). Studies in different populations, such as Brazil and India, have found that the major histocompatibility complex (MHC) region is associated with the risk of leprosy (Covolo de Souza-Santana et al., 2021), (Alter et al., 2011). MHC, also known as human leukocyte antigen (HLA), is the most polymorphic genomic region and correlated strongly with the immune system. HLA-DRB1 is the most closely related gene to leprosy per se as well as clinical forms of the disease. HLA-A, HLA-B, HLA-DQB1 and HLA-DQA1 alleles have also been reported to be associated with leprosy (Fava et al., 2020). In China, a genome-wide association study of leprosy showed that HLA-DRB1*15 was the most significant risk allele (Wang et al., 2016). Nevertheless, due to MHC complex genetic characteristics, elucidating the exact role of MHC in leprosy is challenging.

Copy number variant (CNV) is a type of structural variant that affects DNA copy number, such as large deletions and duplications. CNVs affect approximately 12% of the human genome (Redon et al., 2006). Overall, CNVs are associated with a variety of diseases via various molecular mechanisms, such as Mendelian diseases, HIV susceptibility and autoimmunity (Zhang et al., 2009a). To our knowledge, the association between HLA CNV and leprosy susceptibility remains unclear. To identify more MHC loci, including CNVs, SNPs, HLA alleles and amino acid polymorphisms, associated with leprosy in the Chinese Northern Han population, we extracted MHC data from a previous GWAS including 1,363 leprosy patients and 1,238 controls and conducted HLA imputation based on a large Han-MHC reference panel.

 Methods Study Participants

We used northern population data from a previous genome-wide association study (Wang et al., 2016), (Liu et al., 2015). In short, leprosy samples were collected from dermatology clinics of various hospitals, and all cases were diagnosed by at least two experts. The controls were selected as individuals who did not have histories of leprosy, autoimmune, or systemic disorders or a family history of leprosy (among first-, second-, or third-degree relatives). Patients and controls were matched according to ethnic origin and geographic region of recruitment. The study subjects were homogeneous, and there was no systematic bias or potential population stratification. The study was approved by the ethics committee of the local institutions and conducted in accordance with the Declaration of Helsinki. Informed consent was obtained from all subjects or their family members. All participants were Northern Han Chinese.

 HLA Imputation and Quality Control

We extracted the SNP genotypes located in the MHC region (from 24 to 36 Mb on chromosome 6) from the GWAS data. Based on two reference panels: the 1,000 Genomes Project Phase 3 reference data and the HLA reference panel of the Chinese Han population (10,689 individuals) (Zhou et al., 2016) constructed in the previous study, HLA imputation was performed using the software Beagle 4.1. (Browning and Browning, 2007). To decrease the real computing time and ensure the imputation process was successful, we split each chromosome into chunks and imputed them in parallel in multiple computer processors. For post-imputed quality control, we obtained high-quality variants which satisfied: 1) call rate of >95%; 2) MAF >0.01; 3) Hardy-Weinberg test p value of >1.0 × 10−3; 4) imputation dosage Rsquared value of >0.5 for CNVs, imputation dosage R-squared value > 0.9 for SNPs, and imputation dosage R-squared value >0.7 for HLA amino acids and alleles.

 Association and Stepwise Regression Analyses of the MHC Region

For the association study, the biallelic variants were encoded as allele 1 and allele 2, and the multiallelic variants including multi-residue positions and HLA alleles were encoded as the presence or absence of an individual allele. Stepwise regression analysis was conducted to analyze independent signals of variants related to leprosy. In each conditional regression, the top variant was used as a covariate in the model until there were no significant variants. When the most significant variant was an SNP, HLA alleles, amino acids and CNVs with strong LD associated with it were preferred for stepwise regression analysis; otherwise, we selected the SNP as a covariate. All analyses were performed with PLINK 1.9 software (https://www.cog-genomics.org/plink2).

 In Silico Bioinformatics Analysis

To investigate the function of SNPs, we performed bioinformatics analysis using three online prediction websites: HaploReg (http://pubs.broadinstitute.org/mammals/haploreg/haploreg.php), RegulomeDB (http://regulomedb.org/) and rVarBase (http://rv.psych.ac.cn/). HaploReg was employed to explore noncoding genomic annotations for variants and to identify their potential causal relationship with the pathogenesis of the disease (Ward and Kellis, 2016). RegulomeDB was used to annotate variants with regulatory elements by providing scores (Boyle et al., 2012). The rVarBase database annotates SNPs from three aspects: chromatin states, overlapping regulatory elements and potential target genes (Guo et al., 2016).

 Results HLA Imputation and Association Analyses

After applying quality control, a total of 1,363 leprosy cases and 1,238 controls were included in this study. There were 36,006 high-quality variants after postimputation quality control. Based on Bonferroni adjustment, p < 1.39 × 10−6 (0.05/36,006) was considered the study-wide significance threshold. The top three significant variants were rs9271147, esv3608598 (p = 2.18 × 10−25; OR = 1.98; 95% CI, 1.74–2.25) and DRB1*1501 (p = 4.12 × 10−27; OR = 2.30; 95%CI, 1.98–2.68). After condition analysis, esv3608598 showed strong LD with rs9271471 and DRB1*1501.

 Stepwise Logistic Regression

To identify independent signals with variants in the MHC region, we chose the strongest associated CNV (esv3608598) to conduct stepwise logistic regression. Conditioning on esv3608598 revealed a strong second association with rs7754498 (p = 2.75 × 10−9, OR = 1.42, 95%CI 1.27–1.60). After further condition analysis of rs7754498, the most significant locus was rs3130781 (p = 9.25 × 10−8, or = 1.6, 95% CI 1.34–1.90). When rs3130781 was added to the covariate for logistic regression analysis, rs1443884 was most significant (p = 8.25 × 10−7, OR = 0.69, 95%CI 0.60–0.80) (Figure 1). However, no variant satisfied the significance threshold after conditioning of all the above loci. Thus, esv3608598, rs7754498, rs3130781 and rs144388449 may influence leprosy susceptibility independently (Table1).

Analysis of the HLA associated with leprosy. The x-axis shows genomic position and the y-axis shows–log10(p value) of the variants. The dashed horizontal line represents the significance threshold of p = 1.39 × 10−6. The dots marked red in each panel represent the sites used for condition analysis (esv3608598, rs7754498, rs3130781 and rs144388449).

Stepwise logistic regression of variants associated with leprosy.

STEP Variant Location (hg19) Variant type Raw Stepwise analysis after adjusting Gene annotation
Case Control OR (95% CI) p Value OR (95% CI) p Value
1 esv3608598 Chr6:32571334_32571609 CNV 0.35 0.22 1.98 (1.74–2.25) 2.18 × 10−25 NA NA IntergenicHLA-DRB1, HLA-DQA1
2 rs7754498 Chr6:32453500 SNP 0.55 0.43 1.62 (1.45–1.81) 1.68 × 10−17 1.42 (1.27–1.60) 2.75 × 10−9 IntergenicHLA-DRA, HLA-DRB5
3 rs3130781 Chr6:30914552 SNP 0.14 0.12 1.15 (0.98–1.35) 0.084 1.6 (1.34–1.90) 9.25 × 10−8 Intronic DPCR1
4 rs144388449 Chr6: 31454104 SNP 0.23 0.29 0.74 (0.65–0.84) 1.68 × 10−6 0.69 (0.60–0.80) 8.25 × 10−7 IntergenicHCG26, MICB
Functional Annotation

The results of biological functional annotation are shown in Table 2. The SNP rs7754498 is located 32 kb downstream of HLA-DRB5. The Regulome DB score shows that rs7754498 may be associated with transcription factor binding or DNase peak. The SNP rs3130781 is located in the intron region of DPCR1 and rs144388449 is located 8 kb upstream of MICB. The scores for rs3130781 and rs144388449 provided by Regulome DB were 1f and 2b, respectively, indicating that the two SNPs may be involved in transcription factor binding and DNase peaks. According to HaploReg, three SNPs are related to the regulation of enhancer, promoter and motif changes. Moreover, according to rVarBase results, HCG21, GTF2H4 and VARS2 are target genes of rs3130781.

Functional annotation from bioinformatics analysis.

SNP Regulome DB HaploReg v3.0 rVarBase
Score Enhancer histone marks Proteins bound Motifs changed GENCODEgenes Chromatin state a TF binding rSNPs/rCNVs
rs7754498 5 GI — — 32 kb 3′ of HLA-DRB5 1 0 329
rs3130781 1f GI — Ets DPCR1 intronic 4 0 21
rs144388449 2b — CTCF AP-1,ERalpha-a,SREBP 8.6 kb 5′ of MICB 33 0 147

Regulome DB, score: 5:TF, binding or DNase, peak; 1f: eQTL + TF, binding/DNase, peak; 2b:TF, binding + any motif + DNase, Footprint + DNase, peak.

SNP: Single Nucleotide Polymorphism.

Evidence about the chromatin state of the surrounding region.

 Discussion

Previous GWASs have reported more than 30 independent variants and genes related to leprosy, confirming that the host genetic background plays an important role in leprosy susceptibility. Most leprosy-related genes are associated with immunity, which is consistent with the fact that leprosy is caused by pathogen infection. The MHC region involved in immune responses has also been reported to be associated with leprosy in different populations (Wong et al., 2010; Jarduli et al., 2013). Fine mapping studies of the MHC region have been successfully applied to many diseases, such as systemic lupus erythematosus, psoriasis and vitiligo (Okada et al., 2014; Hanscombe et al., 2018; Yang et al., 2018). On the basis of previous GWAS data, we used Beagle 4.1, an HLA-specific imputation tool, the 1,000 Genome Project Phase 3 dataset and the Han-MHC dataset to conduct HLA imputation and identified four independent variants that may explain the associated of leprosy susceptibility with the MHC region.

In the first leprosy GWAS in Han Chinese individuals, the highest risk was conferred by an SNP at the HLA-DR-DQ locus: rs602875 (Zhang et al., 2009b). Our team previously used the Han-MHC database as a reference panel for fine mapping and identified four independent sites associated with leprosy (HLA-DRβ1 amino acid Tyr26, HLA-C*08:01, HLA-DQA1*03:03) (Zhang et al., 2019). In this study, we observed that a CNV associated with leprosy can further reveal links between the MHC region and leprosy. The CNV esv3608598 is a deletion in the intergenic region of HLA-DRB1 and HLA-DQA1, both of which have been reported to be associated with leprosy risk in many countries (Joko et al., 2000; Vanderborght et al., 2007; Zhang et al., 2009c). In addition, rs9271147 was the strongest variant and showed strong LD with esv3608598.

The SNP rs7754498 is located in the intergenic region of HLA-DRA and HLA-DRB5. A study of leprosy in Indian population showed that BTNL2-DRA intergenic SNPs confer risk for leprosy (Ali et al., 2013). HLA-DRB5 is reported to be associated with susceptibility to leprosy (Koçak et al., 2002). SNP rs3130781 is located in the intron region of DPCR1, which was first identified in 2002 as an MHC class I molecule. DPCR1 modulates NF-κB signaling and plays a role in cell growth (Yan et al., 2018). Moreover, functional annotations by three websites showed that rs3130781 may change DNA motifs and have regulatory effects on transcription. In northern China, DPCR1 rs2844695 confers a high risk of esophageal squamous cell carcinoma (Shen et al., 2014). However, the association between DPCR1 variations and risk of leprosy has not been reported. We think rs3130781 may be a significant variant in future studies to validate the function of DPCR1. The SNP (rs144388449) is located in the intergenic region of HCG26 and MICB. The latter gene, also known as MHC class I polypeptide-related sequence B, encodes a heavily glycosylated protein involved in both innate and adaptive immunity and is located in the MHC class I region. Previous studies have indicated that MICB is closely related to the occurrence of tuberculous leprosy (Tosh et al., 2006). Therefore, our findings suggest that the above three SNPs may be involved in the immune function of MHC and participate in the pathogenesis of leprosy.

There are some limitations in our study. First, the contribution of rare variants (MAF < 0.01) to leprosy was not considered due to the small sample size. Second, HLA variants were not tested in a replication set. In addition, further studies are needed to uncover the underlying mechanism of the MHC region and to confirm its functional role in leprosy.

In conclusion, this is the first study to assess the association between HLA CNV and leprosy susceptibility in a Chinese Northern Han population. In addition to a new CNV, our study reports three independent SNPs in the MHC region that are associated with leprosy risk. The variants we identified emphasize the relationship between host genetic factors and leprosy.

 Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author

 Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of the First Affiliated Hospital of Anhui Medical University. The participants provided their written informed consent to participate in this study.

 Author Contributions

LS were responsible for the study design and revised the manuscript. RZ and LC performed data analysis and wrote the manuscript. WC and HG analyzed GWAS data. XH and ZL carried out genotyping. YW and WF conducted HLA imputation. LY and YY performed stepwise regression analysis. YM and QZ performed bioinformatics analysis. HL and FZ revised the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This study was funded by Anhui Institute of Translational Medicine (ZHYX2020A005), the University Synergy Innovation Program of Anhui Province (GXXT-2020-064), Clinical medicine discipline construction project of Anhui Medical University (2021lcxk008).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

We are grateful to the patients and their family members for participating in this study.

 References Ali S. Srivastava A. K. Chopra R. Aggarwal S. Garg V. K. Bhattacharya S. N. (2013). IL12BSNPs and Copy Number Variation inIL23Rgene Associated with Susceptibility to Leprosy. J. Med. Genet. 50 (1), 34–42. 10.1136/jmedgenet-2012-101214 Alter A. Huong N. T. Singh M. Orlova M. Van Thuc N. Katoch K. (2011). Human Leukocyte Antigen Class I Region Single-Nucleotide Polymorphisms Are Associated with Leprosy Susceptibility in Vietnam and IndiaThe Journal of Infectious Diseases. J. Infect. Dis. 203 (9), 1274–1281. 10.1093/infdis/jir024%J Boyle A. P. Hong E. L. Hariharan M. Cheng Y. Schaub M. A. Kasowski M. (2012). Annotation of Functional Variation in Personal Genomes Using RegulomeDB. Genome Res. 22 (9), 1790–1797. 10.1101/gr.137323.112 Browning S. R. Browning B. L. (2007). Rapid and Accurate Haplotype Phasing and Missing-Data Inference for Whole-Genome Association Studies by Use of Localized Haplotype Clustering. Am. J. Hum. Genet. 81 (5), 1084–1097. 10.1086/521987 Cole S. T. Eiglmeier K. Parkhill J. James K. D. Thomson N. R. Wheeler P. R. (2001). Massive Gene Decay in the Leprosy bacillus. Nature 409 (6823), 1007–1011. 10.1038/35059006 Covolo de Souza-Santana F. Querino G. A. Mendes Camargo R. Nieto Brito de Souza V. Bettoni Ballallai Mangilli P. Távora Mira M. (2021). HLA-DPB1 and HLA-C Alleles Are Associated with Leprosy in a Brazilian Population. Hum. Immunol. 82 (1), 11–18. 10.1016/j.humimm.2020.10.008 Fava V. M. Dallmann-Sauer M. Schurr E. (2020). Genetics of Leprosy: Today and beyond. Hum. Genet. 139 (6-7), 835–846. 10.1007/s00439-019-02087-5 Guo L. Du Y. Qu S. Wang J. (2016). rVarBase: an Updated Database for Regulatory Features of Human Variants. Nucleic Acids Res. 44 (D1), D888–D893. 10.1093/nar/gkv1107 Hanscombe K. B. Morris D. L. Noble J. A. Dilthey A. T. Tombleson P. Kaufman K. M. (2018). Genetic fine Mapping of Systemic Lupus Erythematosus MHC Associations in Europeans and African Americans. Hum. Mol. Genet. 27 (21), 3813–3824. 10.1093/hmg/ddy280 Jarduli L. R. Sell A. M. Reis P. G. Sippert E. Â. Ayo C. M. Mazini P. S. (2013). Role ofHLA,KIR,MICA, and Cytokines Genes in Leprosy. Biomed. Res. Int. 2013, 1–17. 10.1155/2013/989837 Joko S. Numaga J. Kawashima H. Namisato M. Maeda H. (2000). Human Leukocyte Antigens in Forms of Leprosy Among Japanese Patients. Int. J. Lepr. Other Mycobact Dis. 68 (1), 49–56. Koçak M. Balcı M. Pençe B. Kundakçı N. (2002). Associations between Human Leukocyte Antigens and Leprosy in the Turkish Population. Clin. Exp. Dermatol. 27 (3), 235–239. 10.1046/j.1365-2230.2002.01004.x Liu H. Irwanto A. Fu X. a. Yu G. Yu Y. Sun Y. (2015). Discovery of Six New Susceptibility Loci and Analysis of Pleiotropic Effects in Leprosy. Nat. Genet. 47 (3), 267–271. 10.1038/ng.3212 Misch E. A. Berrington W. R. Vary J. C. Hawn T. R. (2010). Leprosy and the Human Genome. Microbiol. Mol. Biol. Rev. 74 (4), 589–620. 10.1128/MMBR.00025-10 Moet F. J. Pahan D. Schuring R. P. Oskam L. Richardus J. H. (2006). Physical Distance, Genetic Relationship, Age, and Leprosy Classification Are Independent Risk Factors for Leprosy in Contacts of Patients with Leprosy. J. Infect. Dis. 193 (3), 346–353. 10.1086/499278 Moraes M. O. Düppre N. C. (2021). Leprosy post-exposure Prophylaxis: Innovation and Precision Public Health. Lancet Glob. Health 9 (1), e8–e9. 10.1016/S2214-109X(20)30512-X Okada Y. Han B. Tsoi L. C. Stuart P. E. Ellinghaus E. Tejasvi T. (2014). Fine Mapping Major Histocompatibility Complex Associations in Psoriasis and its Clinical Subtypes. Am. J. Hum. Genet. 95 (2), 162–172. 10.1016/j.ajhg.2014.07.002 Parkash O. (2009). Classification of Leprosy into Multibacillary and Paucibacillary Groups: an Analysis. FEMS Immunol. Med. Microbiol. 55 (1), 1–5. 10.1111/j.1574-695X.2008.00491.x Redon R. Ishikawa S. Fitch K. R. Feuk L. Perry G. H. Andrews T. D. (2006). Global Variation in Copy Number in the Human Genome. Nature 444 (7118), 444–454. 10.1038/nature05329 Ridley D. S. Jopling W. H. (1966). Classification of Leprosy According to Immunity. A Five-Group System. Int. J. Lepr. Other Mycobact Dis. 34 (3), 255–273. Shen F.-F. Yue W.-B. Zhou F.-Y. Pan Y. Zhao X.-K. Jin Y. (2014). Variations in the MHC Region Confer Risk to Esophageal Squamous Cell Carcinoma on the Subjects from High-Incidence Area in Northern China. PloS one 9 (3), e90438. 10.1371/journal.pone.0090438 Tosh K. Ravikumar M. Bell J. T. Meisner S. Hill A. V. S. Pitchappan R. (2006). Variation in MICA and MICB Genes and Enhanced Susceptibility to Paucibacillary Leprosy in South India. Hum. Mol. Genet. 15 (19), 2880–2887. 10.1093/hmg/ddl229 Vanderborght P. R. Pacheco A. G. Moraes M. E. Antoni G. Romero M. Verville A. (2007). HLA-DRB1*04 and DRB1*10 Are Associated with Resistance and Susceptibility, Respectively, in Brazilian and Vietnamese Leprosy Patients. Genes Immun. 8 (4), 320–324. 10.1038/sj.gene.6364390 Wang Z. Sun Y. Fu X. a. Yu G. Wang C. Bao F. (2016). A Large-Scale Genome-wide Association and Meta-Analysis Identified Four Novel Susceptibility Loci for Leprosy. Nat. Commun. 7, 13760. 10.1038/ncomms13760 Ward L. D. Kellis M. (2016). HaploReg V4: Systematic Mining of Putative Causal Variants, Cell Types, Regulators and Target Genes for Human Complex Traits and Disease. Nucleic Acids Res. 44 (D1), D877–D881. 10.1093/nar/gkv1340 White C. Franco-Paredes C. (2015). Leprosy in the 21st century. Clin. Microbiol. Rev. 28 (1), 80–94. 10.1128/CMR.00079-13 Who (1994). Chemotherapy of Leprosy. Report of a WHO Study Group. Geneva, Switzerland: World Health Organ Tech Rep Ser, 847. Who (2020). Global Leprosy (Hansen Disease) Update, 2019 Time to Step-Up Prevention Initiatives. Geneva, Switzerland: WHO. Wong S. H. Gochhait S. Malhotra D. Pettersson F. H. Teo Y. Y. Khor C. C. (2010). Leprosy and the Adaptation of Human Toll-like Receptor 1. Plos Pathog. 6 (7), e1000979. 10.1371/journal.ppat.1000979 Yan J. Chen G. Zhao X. Chen F. Wang T. Miao F. (2018). High Expression of Diffuse Panbronchiolitis Critical Region 1 Gene Promotes Cell Proliferation, Migration and Invasion in Pancreatic Ductal Adenocarcinoma. Biochem. Biophysical Res. Commun. 495 (2), 1908–1914. 10.1016/j.bbrc.2017.12.031 Yang C. Wu J. Zhang X. Wen L. Sun J. Cheng Y. (2018). Fine-mapping Analysis of the MHC Region for Vitiligo Based on a New Han-MHC Reference Panel. Gene 648, 76–81. 10.1016/j.gene.2018.01.053 Zhang F.-R. Huang W. Chen S.-M. Sun L.-D. Liu H. Li Y. (2009). Genomewide Association Study of Leprosy. N. Engl. J. Med. 361 (27), 2609–2618. 10.1056/NEJMoa0903753 Zhang F. Gu W. Hurles M. E. Lupski J. R. (2009). Copy Number Variation in Human Health, Disease, and Evolution. Annu. Rev. Genom. Hum. Genet. 10, 451–481. 10.1146/annurev.genom.9.081307.164217 Zhang F. Liu H. Chen S. Wang C. Zhu C. Zhang L. (2009). Evidence for an Association of HLA-Drb1*15 and DRB1*09 with Leprosy and the Impact of DRB1*09 on Disease Onset in a Chinese Han Population. BMC Med. Genet. 10, 133. 10.1186/1471-2350-10-133 Zhang X. Cheng Y. Zhang Q. Wang X. Lin Y. Yang C. (2019). Meta-Analysis Identifies Major Histocompatiblity Complex Loci in or Near HLA-DRB1, HLA-DQA1, HLA-C as Associated with Leprosy in Chinese Han Population. J. Invest. Dermatol. 139 (4), 957–960. 10.1016/j.jid.2018.09.029 Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. (2016). Deep Sequencing of the MHC Region in the Chinese Population Contributes to Studies of Complex Disease. Nat. Genet. 48 (7), 740–746. 10.1038/ng.3576