HLA-SPREAD: a natural language processing based resource for curating HLA association from PubMed abstracts

Extreme complexity in the Human Leukocyte Antigens (HLA) system and its nomenclature makes it difficult to interpret and integrate relevant information for HLA associations with diseases, Adverse Drug Reactions (ADR) and Transplantation. PubMed search displays ~ 146,000 studies on HLA reported from diverse locations. Currently, IPD-IMGT/HLA (Robinson et al., Nucleic Acids Research 48:D948–D955, 2019) database houses data on 28,320 HLA alleles. We developed an automated pipeline with a unified graphical user interface HLA-SPREAD that provides a structured information on SNPs, Populations, REsources, ADRs and Diseases information. Information on HLA was extracted from ~ 28 million PubMed abstracts extracted using Natural Language Processing (NLP). Python scripts were used to mine and curate information on diseases, filter false positives and categorize to 24 tree hierarchical groups and named Entity Recognition (NER) algorithms followed by semantic analysis to infer HLA association(s). This resource from 109 countries and 40 ethnic groups provides interesting insights on: markers associated with allelic/haplotypic association in autoimmune, cancer, viral and skin diseases, transplantation outcome and ADRs for hypersensitivity. Summary information on clinically relevant biomarkers related to HLA disease associations with mapped susceptible/risk alleles are readily retrievable from HLASPREAD. The resource is available at URL http://hla-spread.igib.res.in/. This resource is first of its kind that can help uncover novel patterns in HLA gene-disease associations.

Protective effects of HLA-DRB1*04 subtypes in Parkinson's and Alzheimer's diseases implicate acetylated Tau PHF6 sequences

Using genome-wide association data, we analyzed Human Leukocyte Antigen (HLA) associations in over 176,000 individuals with Parkinson's (PD) or Alzheimer's (AD) disease versus controls across ancestry groups. A shared genetic association was observed across diseases at rs601945 (PD: odds ratio (OR)=0.84; 95% confidence interval, [0.80; 0.88]; p=2.2x10-13; AD: OR=0.91[0.89; 0.93]; p=1.8x10-22), and with a protective HLA association recently reported in amyotrophic lateral sclerosis (ALS). Hierarchical protective effects of HLA-DRB1*04 subtypes best accounted for the association, strongest with HLA-DRB1*04:04 and HLA-DRB1*04:07, intermediary with HLA-DRB1*04:01 and HLA-DRB1*04:03, and absent for HLA-DRB1*04:05. The same signal was associated with decreased neurofibrillary tangles (but not neuritic plaque density) in postmortem brain and was more strongly associated with Tau levels than Aβ42 levels in the cerebrospinal fluid. Finally, protective HLA-DRB1*04 subtypes strongly bound aggregation-prone Tau PHF6 sequence, but only when acetylated at K311, a modification central to aggregation. A HLA-DRB1*04-mediated adaptive immune response, potentially against Tau, decreases PD, AD and ALS risk, offering the possibility of new therapeutic avenues.

Differential HLA Association of GAD65 and IA2 Autoantibodies in North Indian Type 1 Diabetes Patients

The human leucocyte antigen (HLA) association with type 1 diabetes (T1D) is well known but there are limited studies investigating the association between β-cell autoantibodies and HLA genes. We evaluated the prevalence of GAD65 and IA-2 autoantibodies (GADA and IA2A) in 252 T1D patients from North India and investigated the genetic association of GADA and IA2A with HLA class I and class II genes/haplotypes. GADA and IA2A were detected in 50.79% and 15.87% of T1D patients, respectively, while only 8.73% had both GADA and IA2A. HLA-DRB1 ∗ 03 was observed to be significantly higher in GADA+ T1D patients as compared to GADA– (91.41% vs. 66.13%, Bonferroni- corrected P P c = 1.11 × 10 − 5 ; OR = 5.45 ; 95% CI: 2.67-11.08). Similarly, HLA-DQB1 ∗ 02 was found to be significantly increased in GADA+ patients (94.53%, P c = 2.19 × 10 − 5 ; OR = 6.27 ; 95% CI: 2.7-14.49) as compared to GADA– (73.39%). The frequencies of HLA-DRB1 ∗ 04 and DQB1 ∗ 03 were increased in IA2A+ patients (45.0% and 52.5%, respectively) as compared to that in IA2A– (25.94% and 33.96%, respectively). Further, the frequency of DRB1 ∗ 03-DQB1 ∗ 02 haplotype was found to be significantly increased in GADA+ T1D patients as compared to GADA- (60.55% vs. 41.94%, P = 3.94 × 10 − 5 ; OR = 2.13 ; 95 % CI = 1.49 -3.03). Similarly, HLA-DRB1 ∗ 04-DQB1 ∗ 03 haplotype was found to be significantly increased in IA2A+ T1D patients compared to IA2A– patients (22.5% vs. 12.97%; P = 0.041 ; OR = 1.95 ; 95 % CI = 1.08 -3.52). None of the HLA class I genes (HLA-A, B, and Cw) was found to be associated with GADA or IA2A in people with T1D. Our findings suggest that HLA-DRB1 ∗ 03/DQB1 ∗ 02 and HLA-DRB1 ∗ 04/DQB1 ∗ 03 might play an important role in the development of GADA and IA2A, respectively.

HLA polymorphisms and risk of glioblastoma in Koreans

Purpose Immune responses for cancer cells can be altered according to genetic variation of human leukocyte antigen (HLA). Association of HLA polymorphism with risk of various cancer types is well known. However, the association between HLA and glioblastoma (GBM) remains uncertain. We sought to evaluate the association of HLA polymorphism with risk of GBM development in Koreans. Materials and methods A case-control study was performed to identify the odds ratios (OR) of HLA class I and II genes for GBM. The control group consisted of 142 healthy Korean volunteers, and the GBM group was 80 patients with newly diagnosed GBM at our institution. HLA class I (-A, -B, and–C) and class II (-DR, -DQ, and–DP) genotyping was performed by high-resolution polymerase chain reaction (PCR)-sequence-based typing (PCR-SBT) methods. Results There were significantly decreased frequencies of HLA-A*26:02 (OR 0.22 CI 0.05–0.98), HLA-C*08:01 (OR 0.29 CI 0.10–0.87), and HLA-DRB1*08:03 (OR 0.32 CI 0.11–0.98), while there was significantly increased frequency of HLA-C*04:01 (OR 2.29 CI 1.05–4.97). In analysis of haplotypes, the frequency of DRB1*14:05-DQB1*05:03 was significantly decreased (OR 0.22 CI 0.05–0.98). Conclusion This study suggests that genetic variations of HLA may affect GBM development in Koreans. Further investigations with larger sample sizes are needed to delineate any potential role of the HLA polymorphisms in the pathogenesis of GBM development.

Approaching Genetics Through the MHC Lens: Tools and Methods for HLA Research

The current SARS-CoV-2 pandemic era launched an immediate and broad response of the research community with studies both about the virus and host genetics. Research in genetics investigated HLA association with COVID-19 based on in silico, population, and individual data. However, they were conducted with variable scale and success; convincing results were mostly obtained with broader whole-genome association studies. Here, we propose a technical review of HLA analysis, including basic HLA knowledge as well as available tools and advice. We notably describe recent algorithms to infer and call HLA genotypes from GWAS SNPs and NGS data, respectively, which opens the possibility to investigate HLA from large datasets without a specific initial focus on this region. We thus hope this overview will empower geneticists who were unfamiliar with HLA to run MHC-focused analyses following the footsteps of the Covid-19|HLA & Immunogenetics Consortium.

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles

ObjectivesDrug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still’s disease with atypical lung disease. We sought to characterise features of patients with Still’s disease with DRESS compared with drug-tolerant Still’s controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort.MethodsIn a case/control study, we collected a multicentre series of patients with Still’s disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still’s controls (n=65). We retrospectively analysed clinical data from all Still’s subjects and typed 94/131 for HLA. European Still’s-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still’s cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still’s-DRESS cases (n=64) compared with drug-tolerant Still’s controls (n=30). KD subjects (n=19) were similarly studied.ResultsStill’s-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still’s-DRESS (64%) versus drug-tolerant Still’s (3%; p=1.2×10−14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still’s-DRESS cases versus INCHARGE Still’s controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still’s controls (p=6.3×10−10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions.ConclusionsDRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.