Human MCTS1-dependent translation of JAK2 is essential for IFN-γ immunity to mycobacteria.

Human inherited disorders of interferon-gamma (IFN-γ) immunity underlie severe mycobacterial diseases. We report X-linked recessive MCTS1 deficiency in men with mycobacterial disease from kindreds of different ancestries (from China, Finland, Iran, and Saudi Arabia). Complete deficiency of this translation re-initiation factor impairs the translation of a subset of proteins, including the kinase JAK2 in all cell types tested, including T lymphocytes and phagocytes. JAK2 expression is sufficiently low to impair cellular responses to interleukin-23 (IL-23) and partially IL-12, but not other JAK2-dependent cytokines. Defective responses to IL-23 preferentially impair the production of IFN-γ by innate-like adaptive mucosal-associated invariant T cells (MAIT) and γδ T lymphocytes upon mycobacterial challenge. Surprisingly, the lack of MCTS1-dependent translation re-initiation and ribosome recycling seems to be otherwise physiologically redundant in these patients. These findings suggest that X-linked recessive human MCTS1 deficiency underlies isolated mycobacterial disease by impairing JAK2 translation in innate-like adaptive T lymphocytes, thereby impairing the IL-23-dependent induction of IFN-γ.

Inborn Errors of RNA Lariat Metabolism in Humans with Brainstem Viral Infection.

Viruses that are typically benign sometimes invade the brainstem in otherwise healthy children. We report bi-allelic DBR1 mutations in unrelated patients from different ethnicities, each of whom had brainstem infection due to herpes simplex virus 1 (HSV1), influenza virus, or norovirus. DBR1 encodes the only known RNA lariat debranching enzyme. We show that DBR1 expression is ubiquitous, but strongest in the spinal cord and brainstem. We also show that all DBR1 mutant alleles are severely hypomorphic, in terms of expression and function. The fibroblasts of DBR1-mutated patients contain higher RNA lariat levels than control cells, this difference becoming even more marked during HSV1 infection. Finally, we show that the patients' fibroblasts are highly susceptible to HSV1. RNA lariat accumulation and viral susceptibility are rescued by wild-type DBR1. Autosomal recessive, partial DBR1 deficiency underlies viral infection of the brainstem in humans through the disruption of tissue-specific and cell-intrinsic immunity to viruses.

Human TMEFF1 is a restriction factor for herpes simplex virus in the brain.

Most cases of herpes simplex virus 1 (HSV-1) encephalitis (HSE) remain unexplained1,2. Here, we report on two unrelated people who had HSE as children and are homozygous for rare deleterious variants of TMEFF1, which encodes a cell membrane protein that is preferentially expressed by brain cortical neurons. TMEFF1 interacts with the cell-surface HSV-1 receptor NECTIN-1, impairing HSV-1 glycoprotein D- and NECTIN-1-mediated fusion of the virus and the cell membrane, blocking viral entry. Genetic TMEFF1 deficiency allows HSV-1 to rapidly enter cortical neurons that are either patient specific or derived from CRISPR-Cas9-engineered human pluripotent stem cells, thereby enhancing HSV-1 translocation to the nucleus and subsequent replication. This cellular phenotype can be rescued by pretreatment with type I interferon (IFN) or the expression of exogenous wild-type TMEFF1. Moreover, ectopic expression of full-length TMEFF1 or its amino-terminal extracellular domain, but not its carboxy-terminal intracellular domain, impairs HSV-1 entry into NECTIN-1-expressing cells other than neurons, increasing their resistance to HSV-1 infection. Human TMEFF1 is therefore a host restriction factor for HSV-1 entry into cortical neurons. Its constitutively high abundance in cortical neurons protects these cells from HSV-1 infection, whereas inherited TMEFF1 deficiency renders them susceptible to this virus and can therefore underlie HSE.

Human genetic and immunological determinants of critical COVID-19 pneumonia.

SARS-CoV-2 infection is benign in most individuals but, in around 10% of cases, it triggers hypoxaemic COVID-19 pneumonia, which leads to critical illness in around 3% of cases. The ensuing risk of death (approximately 1% across age and gender) doubles every five years from childhood onwards and is around 1.5 times greater in men than in women. Here we review the molecular and cellular determinants of critical COVID-19 pneumonia. Inborn errors of type I interferons (IFNs), including autosomal TLR3 and X-chromosome-linked TLR7 deficiencies, are found in around 1-5% of patients with critical pneumonia under 60 years old, and a lower proportion in older patients. Pre-existing auto-antibodies neutralizing IFNα, IFNß and/or IFNω, which are more common in men than in women, are found in approximately 15-20% of patients with critical pneumonia over 70 years old, and a lower proportion in younger patients. Thus, at least 15% of cases of critical COVID-19 pneumonia can be explained. The TLR3- and TLR7-dependent production of type I IFNs by respiratory epithelial cells and plasmacytoid dendritic cells, respectively, is essential for host defence against SARS-CoV-2. In ways that can depend on age and sex, insufficient type I IFN immunity in the respiratory tract during the first few days of infection may account for the spread of the virus, leading to pulmonary and systemic inflammation.

Human autoantibodies neutralizing type I IFNs: From 1981 to 2023.

Human autoantibodies (auto-Abs) neutralizing type I IFNs were first discovered in a woman with disseminated shingles and were described by Ion Gresser from 1981 to 1984. They have since been found in patients with diverse conditions and are even used as a diagnostic criterion in patients with autoimmune polyendocrinopathy syndrome type 1 (APS-1). However, their apparent lack of association with viral diseases, including shingles, led to wide acceptance of the conclusion that they had no pathological consequences. This perception began to change in 2020, when they were found to underlie about 15% of cases of critical COVID-19 pneumonia. They have since been shown to underlie other severe viral diseases, including 5%, 20%, and 40% of cases of critical influenza pneumonia, critical MERS pneumonia, and West Nile virus encephalitis, respectively. They also seem to be associated with shingles in various settings. These auto-Abs are present in all age groups of the general population, but their frequency increases with age to reach at least 5% in the elderly. We estimate that at least 100 million people worldwide carry auto-Abs neutralizing type I IFNs. Here, we briefly review the history of the study of these auto-Abs, focusing particularly on their known causes and consequences.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites.

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

Allele-dependent interaction of LRRK2 and NOD2 in leprosy.

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn's disease and Parkinson's disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

A computational approach for detecting physiological homogeneity in the midst of genetic heterogeneity.

The human genetic dissection of clinical phenotypes is complicated by genetic heterogeneity. Gene burden approaches that detect genetic signals in case-control studies are underpowered in genetically heterogeneous cohorts. We therefore developed a genome-wide computational method, network-based heterogeneity clustering (NHC), to detect physiological homogeneity in the midst of genetic heterogeneity. Simulation studies showed our method to be capable of systematically converging genes in biological proximity on the background biological interaction network, and capturing gene clusters harboring presumably deleterious variants, in an efficient and unbiased manner. We applied NHC to whole-exome sequencing data from a cohort of 122 individuals with herpes simplex encephalitis (HSE), including 13 individuals with previously published monogenic inborn errors of TLR3-dependent IFN-α/ß immunity. The top gene cluster identified by our approach successfully detected and prioritized all causal variants of five TLR3 pathway genes in the 13 previously reported individuals. This approach also suggested candidate variants of three reported genes and four candidate genes from the same pathway in another ten previously unstudied individuals. TLR3 responsiveness was impaired in dermal fibroblasts from four of the five individuals tested, suggesting that the variants detected were causal for HSE. NHC is, therefore, an effective and unbiased approach for unraveling genetic heterogeneity by detecting physiological homogeneity.

Neurologic complications in herpes simplex encephalitis: clinical, immunological and genetic studies.

Patients with herpes simplex virus (HSV) encephalitis (HSE) often develop neuronal autoantibody-associated encephalitis (AE) post-infection. Risk factors of AE are unknown. We tested the hypotheses that predisposition for AE post-HSE may be involved, including genetic variants at specific loci, human leucocyte (HLA) haplotypes, or the blood innate immune response against HSV, including type I interferon (IFN) immunity. Patients of all ages with HSE diagnosed between 1 January 2014 and 31 December 2021 were included in one of two cohorts depending on whether the recruitment was at HSE onset (Spanish Cohort A) or by the time of new neurological manifestations (international Cohort B). Patients were assessed for the type of neurological syndromes; HLA haplotypes; blood type I-IFN signature [RNA quantification of 6 or 28 IFN-response genes (IRG)] and toll-like receptor (TLR3)-type I IFN-related gene mutations. Overall, 190 patients (52% male) were recruited, 93 in Cohort A and 97 in Cohort B. Thirty-nine (42%) patients from Cohort A developed neuronal autoantibodies, and 21 (54%) of them developed AE. Three syndromes (choreoathetosis, anti-NMDAR-like encephalitis and behavioural-psychiatric) showed a high (≥95% cases) association with neuronal autoantibodies. Patients who developed AE post-HSE were less likely to carry the allele HLA-A*02 (4/21, 19%) than those who did not develop AE (42/65, 65%, P = 0.0003) or the Spanish general population (2005/4335, 46%, P = 0.0145). Blood IFN signatures using 6 or 28 IRG were positive in 19/21 (91%) and 18/21 (86%) patients at HSE onset, and rapidly decreased during follow-up. At Day 21 after HSE onset, patients who later developed AE had higher median IFN signature compared with those who did not develop AE [median Zs-6-IRG 1.4 (0.6; 2.0) versus 0.2 (-0.4; 0.8), P = 0.03]. However, a very high median Zs-6-IRG (>4) or persistently increased IFN signature associated with uncontrolled viral infection. Whole exome sequencing showed that the percentage of TLR3-IFN-related mutations in patients who developed AE was not different from those who did not develop AE [3/37 (8%) versus 2/57 (4%), P = 0.379]. Multivariate logistic regression showed that a moderate increase of the blood IFN signature at Day 21 (median Zs-6-IRG >1.5 but <4) was the most important predictor of AE post-HSE [odds ratio 34.8, interquartile ratio (1.7-691.9)]. Altogether, these findings show that most AE post-HSE manifest with three distinct syndromes, and HLA-A*02, but not TLR3-IFN-related mutations, confer protection from developing AE. In addition to neuronal autoantibodies, the blood IFN signature in the context of HSE may be potentially useful for the diagnosis and monitoring of HSE complications.

Genome-wide association study of resistance to Mycobacterium tuberculosis infection identifies a locus at 10q26.2 in three distinct populations.

The natural history of tuberculosis (TB) is characterized by a large inter-individual outcome variability after exposure to Mycobacterium tuberculosis. Specifically, some highly exposed individuals remain resistant to M. tuberculosis infection, as inferred by tuberculin skin test (TST) or interferon-gamma release assays (IGRAs). We performed a genome-wide association study of resistance to M. tuberculosis infection in an endemic region of Southern Vietnam. We enrolled household contacts (HHC) of pulmonary TB cases and compared subjects who were negative for both TST and IGRA (n = 185) with infected individuals (n = 353) who were either positive for both TST and IGRA or had a diagnosis of TB. We found a genome-wide significant locus on chromosome 10q26.2 with a cluster of variants associated with strong protection against M. tuberculosis infection (OR = 0.42, 95%CI 0.35-0.49, P = 3.71×10-8, for the genotyped variant rs17155120). The locus was replicated in a French multi-ethnic HHC cohort and a familial admixed cohort from a hyper-endemic area of South Africa, with an overall OR for rs17155120 estimated at 0.50 (95%CI 0.45-0.55, P = 1.26×10-9). The variants are located in intronic regions and upstream of C10orf90, a tumor suppressor gene which encodes an ubiquitin ligase activating the transcription factor p53. In silico analysis showed that the protective alleles were associated with a decreased expression in monocytes of the nearby gene ADAM12 which could lead to an enhanced response of Th17 lymphocytes. Our results reveal a novel locus controlling resistance to M. tuberculosis infection across different populations.

Clinical, Immunological, and Genetic Findings in Iranian Patients with MHC-II Deficiency: Confirmation of c.162delG RFXANK Founder Mutation in the Iranian Population.

PURPOSE: Major histocompatibility complex class II (MHC-II) deficiency is a rare inborn error of immunity (IEI). Impaired antigen presentation to CD4 + T cells results in combined immunodeficiency (CID). Patients typically present with severe respiratory and gastrointestinal tract infections at early ages. Hematopoietic stem cell transplantation (HSCT) is the only curative therapy. METHODS: We describe the clinical, immunologic, and genetic features of eighteen unrelated Iranian patients with MHC-II deficiency. RESULTS: Consanguinity was present in all affected families. The median age at the initial presentation was 5.5 months (range 7 days to 18 years). The main symptoms included failure to thrive, persistent diarrhea, and pneumonia. Autoimmune and neurologic features were also documented in about one-third of the patients, respectively. Thirteen patients carried RFXANK gene mutations, two carried RFX5 gene mutations, and three carried a RFXAP gene mutation. Six patients shared the same RFXANK founder mutation (c.162delG); limited to the Iranian population and dated to approximately 1296 years ago. Four of the patients underwent HSCT; three of them are alive. On the other hand, nine of the fourteen patients who did not undergo HSCT had a poor prognosis and died. CONCLUSION: MHC-II deficiency is not rare in Iran, with a high rate of consanguinity. It should be considered in the differential diagnosis of CID at any age. With the limited access to HSCT and its variable results in MHC-II deficiency, implementing genetic counseling and family planning for the affected families are mandatory. We are better determined to study the c.162delG RFXANK heterozygous mutation frequency in the Iranian population.

Autoantibodies Neutralizing Type I IFNs in the Bronchoalveolar Lavage of at Least 10% of Patients During Life-Threatening COVID-19 Pneumonia.

Autoantibodies (auto-Abs) neutralizing type I interferons (IFNs) are found in the blood of at least 15% of unvaccinated patients with life-threatening COVID-19 pneumonia. We report here the presence of auto-Abs neutralizing type I IFNs in the bronchoalveolar lavage (BAL) of 54 of the 415 unvaccinated patients (13%) with life-threatening COVID-19 pneumonia tested. The 54 individuals with neutralizing auto-Abs in the BAL included 45 (11%) with auto-Abs against IFN-α2, 37 (9%) with auto-Abs against IFN-ω, 54 (13%) with auto-Abs against IFN-α2 and/or ω, and five (1%) with auto-Abs against IFN-ß, including three (0.7%) with auto-Abs neutralizing IFN-α2, IFN-ω, and IFN-ß, and two (0.5%) with auto-Abs neutralizing IFN-α2 and IFN-ß. Auto-Abs against IFN-α2 also neutralize the other 12 subtypes of IFN-α. Paired plasma samples were available for 95 patients. All seven patients with paired samples who had detectable auto-Abs in BAL also had detectable auto-Abs in plasma, and one patient had auto-Abs detectable only in blood. Auto-Abs neutralizing type I IFNs are, therefore, present in the alveolar space of at least 10% of patients with life-threatening COVID-19 pneumonia. These findings suggest that these auto-Abs impair type I IFN immunity in the lower respiratory tract, thereby contributing to hypoxemic COVID-19 pneumonia.

No increased prevalence of autoantibodies neutralizing type I IFNs in idiopathic pulmonary fibrosis patients.

SARS-CoV2 infection has a poor prognosis in patients affected of idiopathic pulmonary fibrosis (IPF). Autoantibodies (auto-Abs) neutralizing type I interferons (IFNs) are found in the blood of at least 15% of patients with life-threatening COVID-19 pneumonia. Because of the elevated prevalence of some auto-Abs in IPF patients, we hypothesize that the prevalence of auto-Abs neutralizing type I IFNs might be increased in the IPF population and then explained specific poor outcome after COVID-19. We screened the plasma of 247 consecutive IPF patients for the presence of auto-Abs neutralizing type I IFNs. Three patients displayed auto-Abs neutralizing type I IFNs. Among them, the only patient with documented SARS-CoV-2 infection experienced life threatening COVID-19 pneumonia. The prevalence of auto-Abs neutralizing type I IFNs in this cohort of IPF patients was not significantly different from the one of the general population. Overall, this study did not suggest any association between auto-Abs neutralizing type I IFNs and IPF.

A genome-wide case-only test for the detection of digenic inheritance in human exomes.

Whole-exome sequencing (WES) has facilitated the discovery of genetic lesions underlying monogenic disorders. Incomplete penetrance and variable expressivity suggest a contribution of additional genetic lesions to clinical manifestations and outcome. Some monogenic disorders may therefore actually be digenic. However, only a few digenic disorders have been reported, all discovered by candidate gene approaches applied to at least one locus. We propose here a two-locus genome-wide test for detecting digenic inheritance in WES data. This approach uses the gene as the unit of analysis and tests all pairs of genes to detect pairwise gene × gene interactions underlying disease. It is a case-only method, which has several advantages over classic case-control tests, in particular by avoiding recruitment of controls. Our simulation studies based on real WES data identified two major sources of type I error inflation in this case-only test: linkage disequilibrium and population stratification. Both were corrected by specific procedures. Moreover, our case-only approach is more powerful than the corresponding case-control test for detecting digenic interactions in various population stratification scenarios. Finally, we confirmed the potential of our unbiased, genome-wide approach by successfully identifying a previously reported digenic lesion in patients with craniosynostosis. Our case-only test is a powerful and timely tool for detecting digenic inheritance in WES data from patients.

Genetic and immunologic evaluation of children with inborn errors of immunity and severe or critical COVID-19.

BACKGROUND: Most severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected individuals are asymptomatic or only exhibit mild disease. In about 10% of cases, the infection leads to hypoxemic pneumonia, although it is much more rare in children. OBJECTIVE: We evaluated 31 young patients aged 0.5 to 19 years who had preexisting inborn errors of immunity (IEI) but lacked a molecular diagnosis and were later diagnosed with coronavirus disease 2019 (COVID-19) complications. METHODS: Genetic evaluation by whole-exome sequencing was performed in all patients. SARS-CoV-2-specific antibodies, autoantibodies against type I IFN (IFN-I), and inflammatory factors in plasma were measured. We also reviewed COVID-19 disease severity/outcome in reported IEI patients. RESULTS: A potential genetic cause of the IEI was identified in 28 patients (90.3%), including mutations that may affect IFN signaling, T- and B-cell function, the inflammasome, and the complement system. From tested patients 65.5% had detectable virus-specific antibodies, and 6.8% had autoantibodies neutralizing IFN-I. Five patients (16.1%) fulfilled the diagnostic criteria of multisystem inflammatory syndrome in children. Eleven patients (35.4%) died of COVID-19 complications. All together, at least 381 IEI children with COVID-19 have been reported in the literature to date. Although many patients with asymptomatic or mild disease may not have been reported, severe presentation of COVID-19 was observed in 23.6% of the published cases, and the mortality rate was 8.7%. CONCLUSIONS: Young patients with preexisting IEI may have higher mortality than children without IEI when infected with SARS-CoV-2. Elucidating the genetic basis of IEI patients with severe/critical COVID-19 may help to develop better strategies for prevention and treatment of severe COVID-19 disease and complications in pediatric patients.

Taking population stratification into account by local permutations in rare-variant association studies on small samples.

Many methods for rare variant association studies require permutations to assess the significance of tests. Standard permutations assume that all individuals are exchangeable and do not take population stratification (PS), a known confounding factor in genetic studies, into account. We propose a novel strategy, LocPerm, in which individual phenotypes are permuted only with their closest ancestry-based neighbors. We performed a simulation study, focusing on small samples, to evaluate and compare LocPerm with standard permutations and classical adjustment on first principal components. Under the null hypothesis, LocPerm was the only method providing an acceptable type I error, regardless of sample size and level of stratification. The power of LocPerm was similar to that of standard permutation in the absence of PS, and remained stable in different PS scenarios. We conclude that LocPerm is a method of choice for taking PS and/or small sample size into account in rare variant association studies.

Autoantibodies Neutralizing Type I Interferons in 20% of COVID-19 Deaths in a French Hospital.

Recent studies reported the presence of pre-existing autoantibodies (auto-Abs) neutralizing type I interferons (IFNs) in at least 15% of patients with critical COVID-19 pneumonia. In one study, these auto-Abs were found in almost 20% of deceased patients across all ages. We aimed to assess the prevalence and clinical impact of the auto-Abs to type I IFNs in the Seine-Saint-Denis district, which was one of the most affected areas by COVID-19 in France during the first wave. We tested for the presence of auto-Abs neutralizing type I IFNs in a cohort of patients admitted for critical COVID-19 pneumonia during the first wave in the spring of 2020 in the medicine departments at Robert Ballanger Hospital, Aulnay sous Bois. We found circulating auto-Abs that neutralized 100 pg/mL IFN-α2 and/or IFN-ω in the plasma (diluted 1/10) of 7.9% (11 of 139) of the patients hospitalized for critical COVID-19. The presence of neutralizing auto-Abs was associated with an increased risk of mortality, as these auto-Abs were detected in 21% of patients who died from COVID-19 pneumonia. Deceased patients with and without auto-Abs did not present overt clinical differences. These results confirm both the importance of type I IFN immunity in host defense against SARS-CoV-2 infection and the usefulness of detection of auto-Abs neutralizing type I IFNs in the management of patients.

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.

The complex pattern of genetic associations of leprosy with HLA class I and class II alleles can be reduced to four amino acid positions.

Leprosy is a chronic disease caused by Mycobacterium leprae. Worldwide, more than 200,000 new patients are affected by leprosy annually, making it the second most common mycobacterial disease after tuberculosis. The MHC/HLA region has been consistently identified as carrying major leprosy susceptibility variants in different populations at times with inconsistent results. To establish the unambiguous molecular identity of classical HLA class I and class II leprosy susceptibility factors, we applied next-generation sequencing to genotype with high-resolution 11 HLA class I and class II genes in 1,155 individuals from a Vietnamese leprosy case-control sample. HLA alleles belonging to an extended haplotype from HLA-A to HLA-DPB1 were associated with risk to leprosy. This susceptibility signal could be reduced to the HLA-DRB1*10:01~ HLA-DQA1*01:05 alleles which were in complete linkage disequilibrium (LD). In addition, haplotypes containing HLA-DRB3~ HLA-DRB1*12:02 and HLA-C*07:06~ HLA-B*44:03~ HLA-DRB1*07:01 alleles were found as two independent protective factors for leprosy. Moreover, we replicated the previously associated HLA-DRB1*15:01 as leprosy risk factor and HLA-DRB1*04:05~HLA-DQA1*03:03 as protective alleles. When we narrowed the analysis to the single amino acid level, we found that the associations of the HLA alleles were largely captured by four independent amino acids at HLA-DRß1 positions 57 (D) and 13 (F), HLA-B position 63 (E) and HLA-A position 19 (K). Hence, analyses at the amino acid level circumvented the ambiguity caused by strong LD of leprosy susceptibility HLA alleles and identified four distinct leprosy susceptibility factors.