Development and application of EpitopeScan, a Python3 toolset for mutation tracking in SARS-CoV-2 immunogenic epitopes.

Introduction: Outbreaks of coronaviruses and especially the recent COVID-19 pandemic emphasize the importance of immunological research in this area to mitigate the effect of future incidents. Bioinformatics approaches are capable of providing multisided insights from virus sequencing data, although currently available software options are not entirely suitable for a specific task of mutation surveillance within immunogenic epitopes of SARS-CoV-2. Method: Here, we describe the development of a mutation tracker, EpitopeScan, a Python3 package with command line and graphical user interface tools facilitating the investigation of the mutation dynamics in SARS-CoV-2 epitopes via analysis of multiple-sequence alignments of genomes over time. We provide an application case by examining three Spike protein-derived immunodominant CD4+ T-cell epitopes restricted by HLA-DRB1*04:01, an allele strongly associated with susceptibility to rheumatoid arthritis (RA). Mutations in these peptides are relevant for immune monitoring of CD4+ T-cell responses against SARS-CoV-2 spike protein in patients with RA. The analysis focused on 2.3 million SARS-CoV-2 genomes sampled in England. Results: We detail cases of epitope conservation over time, partial loss of conservation, and complete divergence from the wild type following the emergence of the N969K Omicron-specific mutation in November 2021. The wild type and the mutated peptide represent potential candidates to monitor variant-specific CD4+ T-cell responses. EpitopeScan is available via GitHub repository https://github.com/Aleksandr-biochem/EpitopeScan.

Case report: Strong GAD antibody positivity and type 1 diabetes-HLA-susceptible haplotype-DRB1*04:05-DQB1*04:01 in a Japanese patient with immune checkpoint inhibitor-induced type 1 diabetes.

Immune checkpoint inhibitors (ICIs) are widely used in cancer treatment; however, they can lead to immune-related adverse events, including immune checkpoint inhibitor-induced type 1 diabetes mellitus (ICI-T1DM). While fulminant T1DM is common in East Asia, ICI-T1DM has predominantly been reported in Western countries. In this report, we present the case of a 66-year-old Japanese man with type 2 diabetes mellitus undergoing dialysis for diabetic nephropathy. The patient was diagnosed with left upper lobe lung cancer, and treatment with nivolumab and ipilimumab was initiated. After 48 days, the patient experienced impaired consciousness and difficulty moving. His blood glucose levels were 815 mg/dL, and metabolic acidosis was detected, leading to a diagnosis of diabetic ketoacidosis. The patient was subsequently treated with continuous intravenous insulin. However, his C-peptide levels rapidly depleted, and new-onset ICI-T1DM was diagnosed. Although most Japanese patients with ICI-T1DM test negative for glutamic acid decarboxylase (GAD) antibodies, this case exhibited a strong positivity. Thus, we reviewed the literature on 15 similar Japanese cases, revealing a mean HbA1c level at onset of 8.7% and a mean time from ICI administration to onset of 9.7 weeks, which was shorter than that in GAD-negative cases. Moreover, human leukocyte antigen typing revealed five cases of DRB1*04:05-DQB1*04:01, including the present case, and one case of DRB1*09:01-DQB1*03:03, both of which were susceptible to T1DM haplotypes. These findings suggest that GAD antibody positivity may be associated with acute onset and disease progression in some cases of Japanese patients with ICI-T1DM. Given that the prediction of new-onset ICI-T1DM is challenging, monitoring GAD antibody levels might be useful. However, further studies with large sample sizes and validation across different racial and ethnic populations are warranted.

Exploring the role of HLA variants in neuroblastoma susceptibility through whole exome sequencing.

Although a number of susceptibility loci for neuroblastoma (NB) have been identified by genome-wide association studies, it is still unclear whether variants in the HLA region contribute to NB susceptibility. In this study, we conducted a comprehensive genetic analysis of variants in the HLA region among 724 NB patients and 2863 matched controls from different cohorts. We exploited whole-exome sequencing data to accurately type HLA alleles with an ensemble approach on the results from three different typing tools, and carried out rigorous sample quality control to ensure a fine-scale ancestry matching. The frequencies of common HLA alleles were compared between cases and controls by logistic regression under additive and non-additive models. Population stratification was taken into account adjusting for ancestry-informative principal components. We detected significant HLA associations with NB. In particular, HLA-DQB1*05:02 (OR = 1.61; padj = 5.4 × 10-3) and HLA-DRB1*16:01 (OR = 1.60; padj = 2.3 × 10-2) alleles were associated to higher risk of developing NB. Conditional analysis highlighted the HLA-DQB1*05:02 allele and its residue Ser57 as key to this association. DQB1*05:02 allele was not associated to clinical features worse outcomes in the NB cohort. Nevertheless, a risk score derived from the allelic combinations of five HLA variants showed a substantial predictive value for patient survival (HR = 1.53; p = 0.032) that was independent from established NB prognostic factors. Our study leveraged powerful computational methods to explore WES data and HLA variants and to reveal complex genetic associations. Further studies are needed to validate the mechanisms of these interactions that contribute to the multifaceted pattern of factors underlying the disease initiation and progression.

Exploring potential correlations between HLA class II and the risk of microvascular complications in Japanese patients with type 1 diabetes.

Managing complications in Type 1 diabetes (T1D) remains challenging. HLA genes, particularly DR and DQ, are linked to T1D susceptibility. We studied 48 Japanese T1D inpatients and revealed associations between DRB1*04:05-DQB1*04:01 and DRB1*09:01-DQB1*03:03 haplotypes and complications, offering a new perspective for future research.

The HLA-DRB1*12:92 and HLA-DRB1*12:101 alleles were identified by next-generation sequencing.

Compared with HLA-DRB1*12:02, the alleles HLA-DRB1*12:92 and HLA-DRB1*12:101 each show one nucleotide substitution respectively.

Human Leukocyte Antigen Alleles (HLA-A, HLA-B, and HLA-DRB1) are associated with Acute Lymphoblastic Leukemia (ALL) : A Case-Control Study in a Sample of Iranian Population.

OBJECTIVE: This study seeks to elucidate the association between HLA-A, HLA-B, and HLA-DRB1 alleles and their relative risk contributions to ALL within an Iranian cohort. METHODS: Utilizing a robust case-control design, this research involved 71 ALL patients and 71 age and sex-matched healthy individuals. Genotyping of specified HLA alleles was performed using the advanced PCR-SSP technique. RESULTS: Our findings reveal a marked increase in the prevalence of the HLA-DRB1*04 allele among patients diagnosed with ALL compared to the control group (P<0.027). Conversely, the alleles HLA-A*26 (P=0.025), HLA-A*33 (P=0.020), and HLA-DRB1*03 (P=0.035) were observed at significantly reduced frequencies within the patient population. CONCLUSION: Our findings highlight HLA-DRB1*04 as a potential genetic marker for increased susceptibility to ALL, while HLA-A*26, HLA-A*33, and HLA-DRB1*03 emerge as protective factors.

Identification of the novel HLA-DRB1*11:298 allele in a Chinese cord blood donor.

HLA-DRB1*11:298 shows one single nucleotide substitution at position 397 T>G compared with HLA-DRB1*11:01:01:01.

Characterisation of two novel HLA-DRB1 alleles, HLA-DRB1*04:01:24 and HLA-DRB1*13:353.

Two novel HLA-DRB1 alleles were detected using next generation sequencing.

High-resolution African HLA resource uncovers HLA-DRB1 expression effects underlying vaccine response.

How human genetic variation contributes to vaccine effectiveness in infants is unclear, and data are limited on these relationships in populations with African ancestries. We undertook genetic analyses of vaccine antibody responses in infants from Uganda (n = 1391), Burkina Faso (n = 353) and South Africa (n = 755), identifying associations between human leukocyte antigen (HLA) and antibody response for five of eight tested antigens spanning pertussis, diphtheria and hepatitis B vaccines. In addition, through HLA typing 1,702 individuals from 11 populations of African ancestry derived predominantly from the 1000 Genomes Project, we constructed an imputation resource, fine-mapping class II HLA-DR and DQ associations explaining up to 10% of antibody response variance in our infant cohorts. We observed differences in the genetic architecture of pertussis antibody response between the cohorts with African ancestries and an independent cohort with European ancestry, but found no in silico evidence of differences in HLA peptide binding affinity or breadth. Using immune cell expression quantitative trait loci datasets derived from African-ancestry samples from the 1000 Genomes Project, we found evidence of differential HLA-DRB1 expression correlating with inferred protection from pertussis following vaccination. This work suggests that HLA-DRB1 expression may play a role in vaccine response and should be considered alongside peptide selection to improve vaccine design.

HLA molecular mismatches and induced donor-specific tolerance in combined living donor kidney and hematopoietic stem cell transplantation.

Introduction: We investigated the potential role of HLA molecular mismatches (MM) in achieving stable chimerism, allowing for donor-specific tolerance in patients undergoing combined living donor kidney and hematopoietic stem cell transplantation (HSCT). Methods: All patients with available DNA samples (N=32) who participated in a phase 2 clinical trial (NCT00498160) where they received an HLA mismatched co-transplantation of living donor kidney and facilitating cell-enriched HSCT were included in this study. High-resolution HLA genotyping data were used to calculate HLA amino acid mismatches (AAMM), Eplet MM, three-dimensional electrostatic mismatch scores (EMS-3D), PIRCHE scores, HLA-DPB1 T-cell epitope group MM, HLA-B leader sequence MM, and KIR ligands MM between the donor and recipient in both directions. HLA MM were analyzed to test for correlation with the development of chimerism, graft vs. host disease (GvHD), de novo DSA, and graft rejection. Results: Follow-up time of this cohort was 6-13.5 years. Of the 32 patients, 26 developed high-level donor or mixed stable chimerism, followed by complete withdrawal of immunosuppression (IS) in 25 patients. The remaining six of the 32 patients had transient chimerism or no engraftment and were maintained on IS (On-IS). In host versus graft direction, a trend toward higher median number of HLA-DRB1 MM scores was seen in patients On-IS compared to patients with high-level donor/mixed chimerism, using any of the HLA MM modalities; however, initial statistical significance was observed only for the EMS-3D score (0.45 [IQR, 0.30-0.61] vs. 0.24 [IQR, 0.18-0.36], respectively; p=0.036), which was lost when applying the Bonferroni correction. No statistically significant differences between the two groups were observed for AAMM, EMS-3D, Eplet MM, and PIRCHE-II scores calculated in graft versus host direction. No associations were found between development of chimerism and GvHD and non-permissive HLA-DPB1 T-cell epitope group MM, HLA-B leader sequence, and KIR ligands MM. Conclusion: Our results suggest an association between HLA-DRB1 molecular mismatches and achieving stable chimerism, particularly when electrostatic quality of the mismatch is considered. The non-permissive HLA-DPB1 T-cell epitope group, HLA-B leader sequence, and KIR ligands MM do not predict chimerism and GvHD in this combined kidney/HSCT transplant patient cohort. Further work is needed to validate our findings. Clinical trial registration: https://clinicaltrials.gov/study/NCT00498160, identifier NCT00498160.

Identification of two novel HLA-DRB1 alleles, HLA-DRB1*08:03:13 and HLA-DRB1*08:119.

Compared with HLA-DRB1*08:03:02:01, the alleles HLA-DRB1*08:03:13 and HLA-DRB1*08:119 each show one nucleotide substitution, respectively.

The novel HLA-DRB1*03:210 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-DRB1*03:210 differs from HLA-DRB1*03:01:01:01 by one non-synonymous nucleotide substitution in exon 3.

The novel HLA-DRB1*11:323 allele characterized by two different sequencing-based typing techniques.

The novel allele HLA-DRB1*11:323 differs from HLA-DRB1*11:01:02:01 by one non-synonymous nucleotide substitution in exon 2.

Family-based association of HLA-DRB1 and DQB1 alleles and haplotypes in a group of Iranian Type 1 diabetes children.

This family-based study was conducted in a group of Iranians with Type 1 diabetes (T1D) to investigate the transmission from parents of risk and non-risk HLA alleles and haplotypes, and to estimate the genetic risk score for this disease within this population. A total of 240 T1D subjects including 111 parent-child trios (111 children with T1D, 133 siblings, and 222 parents) and 330 ethnically matched healthy individuals were recruited. High-resolution HLA typing for DRB1/DQB1 loci was performed for all study subjects (n = 925) using polymerase chain reaction-sequence-specific oligonucleotide probe method. The highest predisposing effect on developing T1D was conferred by the following haplotypes both in all subjects and in probands compared to controls: DRB1*04:05-DQB1*03:02 (Pc = 2.97e-06 and Pc = 6.04e-10, respectively), DRB1*04:02-DQB1*03:02 (Pc = 5.94e-17 and Pc = 3.86e-09, respectively), and DRB1*03:01-DQB1*02:01 (Pc = 8.26e-29 and Pc = 6.56e-16, respectively). Conversely, the major protective haplotypes included DRB1*13:01-DQB1*06:03 (Pc = 6.99e-08), DRB1*15:01-DQB1*06:02 (Pc = 2.97e-06) in the cases versus controls. Also, DRB1*03:01-DQB1*02:01/DRB1*04:02|05-DQB1*03:02 and DRB1*03:01-DQB1*02:01/DRB1*03:01-DQB1*02:01 diplotypes conferred the highest predisposing effect in the cases (Pc = 8.65e-17 and Pc = 6.26e-08, respectively) and in probands (Pc = 5.4e-15 and Pc = 0.001, respectively) compared to controls. Transmission disequilibrium test showed that the highest risk was conferred by DRB1*04:02-DQB1*03:02 (Pc = 3.26e-05) and DRB1*03:01-DQB1*02:01 (Pc = 1.78e-12) haplotypes and the highest protection by DRB1*14:01-DQB1*05:03 (Pc = 8.66e-05), DRB1*15:01-DQB1*06:02 (Pc = 0.002), and DRB1*11:01-DQB1*03:01 (Pc = 0.0003) haplotypes. Based on logistic regression analysis, carriage of risk haplotypes increased the risk of T1D development 24.5 times in the Iranian population (p = 5.61e-13). Also, receiver operating characteristic curve analysis revealed a high predictive power of those risk haplotypes in discrimination of susceptible from healthy individuals (area under curve: 0.88, p = 5.5e-32). Our study highlights the potential utility of genetic risk assessment based on HLA diplotypes for predicting T1D risk in individuals, particularly among family members of affected children in our population.

Novel genetic variants of HLA gene associated with Thai Behcet's disease (BD) patients using next generation sequencing technology.

Behçet's disease (BD) manifests as an autoimmune disorder featuring recurrent ulcers and multi-organ involvement, influenced by genetic factors associated with both HLA and non-HLA genes, including TNF-α and ERAP1. The study investigated the susceptible alleles of both Class I and II molecules of the HLA gene in 56 Thai BD patients and 192 healthy controls through next-generation sequencing using a PacBio kit. The study assessed 56 BD patients, primarily females (58.9%), revealing diverse manifestations including ocular (41.1%), vascular (35.7%), skin (55.4%), CNS (5.4%), and GI system (10.7%) involvement. This study found associations between BD and HLA-A*26:01:01 (OR 3.285, 95% CI 1.135-9.504, P-value 0.028), HLA-B*39:01:01 (OR 6.176, 95% CI 1.428-26.712, P-value 0.015), HLA-B*51:01:01 (OR 3.033, 95% CI 1.135-8.103, P-value 0.027), HLA-B*51:01:02 (OR 6.176, 95% CI 1.428-26.712, P-value 0.015), HLA-C*14:02:01 (OR 3.485, 95% CI 1.339-9.065, P-value 0.01), HLA-DRB1*14:54:01 (OR 1.924, 95% CI 1.051-3.522, P-value 0.034), and HLA-DQB1*05:03:01 (OR 3.00, 95% CI 1.323-6.798, P-value 0.008). However, after Bonferroni correction none of these alleles were found to be associated with BD. In haplotype analysis, we found a strong linkage disequilibrium in HLA-B*51:01:01, HLA-C*14:02:01 (P-value 0.0, Pc-value 0.02). Regarding the phenotype, a significant association was found between HLA-DRB1*14:54:01 (OR 11.67, 95% CI 2.86-47.57, P-value 0.001) and BD with ocular involvement, apart from this, no distinct phenotype-HLA association was documented. In summary, our study identifies specific HLA associations in BD. Although limited by a small sample size, we acknowledge the need for further investigation into HLA relationships with CNS, GI, and neurological phenotypes in the Thai population.

Investigation of HLA susceptibility alleles and genotypes with hematological disease among Chinese Han population.

Several genes involved in the pathogenesis have been identified, with the human leukocyte antigen (HLA) system playing an essential role. However, the relationship between HLA and a cluster of hematological diseases has received little attention in China. Blood samples (n = 123913) from 43568 patients and 80345 individuals without known pathology were genotyped for HLA class I and II using sequencing-based typing. We discovered that HLA-A*11:01, B*40:01, C*01:02, DQB1*03:01, and DRB1*09:01 were prevalent in China. Furthermore, three high-frequency alleles (DQB1*03:01, DQB1*06:02, and DRB1*15:01) were found to be hazardous in malignant hematologic diseases when compared to controls. In addition, for benign hematologic disorders, 7 high-frequency risk alleles (A*01:01, B*46:01, C*01:02, DQB1*03:03, DQB1*05:02, DRB1*09:01, and DRB1*14:54) and 8 high-frequency susceptible genotypes (A*11:01-A*11:01, B*46:01-B*58:01, B*46:01-B*46:01, C*01:02-C*03:04, DQB1*03:01-DQB1*05:02, DQB1*03:03-DQB1*06:01, DRB1*09:01-DRB1*15:01, and DRB1*14:54-DRB1*15:01) were observed. To summarize, our findings indicate the association between HLA alleles/genotypes and a variety of hematological disorders, which is critical for disease surveillance.

Multiple Sclerosis Onset before and after COVID-19 Vaccination: Can HLA Haplotype Be Determinant?

A few cases of multiple sclerosis (MS) onset after COVID-19 vaccination have been reported, although the evidence is insufficient to establish causality. The aim of this study is to compare cases of newly diagnosed relapsing-remitting MS before and after the outbreak of the COVID-19 pandemic and the impact of COVID-19 vaccination. Potential environmental and genetic predisposing factors were also investigated, as well as clinical patterns. This is a single-centre retrospective cohort study including all patients who presented with relapsing-remitting MS onset between January 2018 and July 2022. Data on COVID-19 vaccination administration, dose, and type were collected. HLA-DRB1 genotyping was performed in three subgroups. A total of 266 patients received a new diagnosis of relapsing-remitting MS in our centre, 143 before the COVID-19 pandemic (until and including March 2020), and 123 during the COVID-19 era (from April 2020). The mean number of new MS onset cases per year was not different before and during the COVID-19 era and neither were baseline patients' characteristics, type of onset, clinical recovery, or radiological patterns. Fourteen (11.4%) patients who subsequently received a new diagnosis of MS had a history of COVID-19 vaccination within one month before symptoms onset. Patients' characteristics, type of onset, clinical recovery, and radiological patterns did not differ from those of patients with non-vaccine-related new diagnoses of MS. The allele frequencies of HLA-DRB1*15 were 17.6% and 22.2% in patients with non-vaccine-related disease onset before and during the COVID-19 era, respectively, while no case of HLA-DRB1*15 was identified among patients with a new diagnosis of MS post-COVID-19 vaccine. In contrast, HLA-DRB1*08+ or HLA-DRB1*10+ MS patients were present only in this subgroup. Although a causal link between COVID-19 vaccination and relapsing-remitting MS cannot be detected, it is interesting to note and speculate about the peculiarities and heterogeneities underlying disease mechanisms of MS, where the interactions of genetics and the environment could be crucial also for the follow-up and the evaluation of therapeutic options.

Two novel HLA class II alleles, HLA-DRB1*09:57 and HLA-DRB1*09:58 in Chinese individuals.

Compared with HLA-DRB1*09:01:02:05, the alleles HLA-DRB1*09:57 and HLA-DRB1*09:58 each show one nucleotide change, respectively.

Determination of HLA class II risk alleles and prediction of self/non-self-epitopes contributing Hashimoto's thyroiditis in a group of Iranian patients.

One of the probable hypotheses for the onset of autoimmunity is molecular mimicry. This study aimed to determine the HLA-II risk alleles for developing Hashimoto's thyroiditis (HT) in order to analyze the molecular homology between candidate pathogen-derived epitopes and potentially self-antigens (thyroid peroxidase, TPO) based on the presence of HLA risk alleles. HLA-DRB1/-DQB1 genotyping was performed in 100 HT patients and 330 ethnically matched healthy controls to determine the predisposing/protective alleles for HT disease. Then, in silico analysis was conducted to examine the sequence homology between epitopes derived from autoantigens and four potentially relevant pathogens and their binding capacities to HLA risk alleles based on peptide docking analysis. We identified HLA-DRB1*03:01, *04:02, *04:05, and *11:04 as predisposing alleles and DRB1*13:01 as a potentially predictive allele for HT disease. Also, DRB1*11:04 ~ DQB1*03:01 (Pc = 0.002; OR, 3.97) and DRB1*03:01 ~ DQB1*02:01 (Pc = 0.004; OR, 2.24) haplotypes conferred a predisposing role for HT. Based on logistic regression analysis, carrying risk alleles increased the risk of HT development 4.5 times in our population (P = 7.09E-10). Also, ROC curve analysis revealed a high predictive power of those risk alleles for discrimination of the susceptible from healthy individuals (AUC, 0.70; P = 6.6E-10). Analysis of peptide sequence homology between epitopes of TPO and epitopes derived from four candidate microorganisms revealed a homology between envelop glycoprotein D of herpes virus and sequence 151-199 of TPO with remarkable binding capacity to HLA-DRB1*03:01 allele. Our findings indicate the increased risk of developing HT in those individuals carrying HLA risk alleles which can also be related to herpes virus infection.

Validation of tag SNPs for multiple sclerosis HLA risk alleles across the 1000 genomes panel.

Currently, the genetic variants strongly associated with risk for Multiple Sclerosis (MS) are located in the Major Histocompatibility Complex. This includes DRB1*15:01 and DRB1*15:03 alleles at the HLA-DRB1 locus, the latter restricted to African populations; the DQB1*06:02 allele at the HLA-DQB1 locus which is in high linkage disequilibrium (LD) with DRB1*15:01; and protective allele A*02:01 at the HLA-A locus. HLA allele identification is facilitated by co-inherited ('tag') single nucleotide polymorphisms (SNPs); however, SNP validation is not typically done outside of the discovery population. We examined 19 SNPs reported to be in high LD with these alleles in 2,502 healthy subjects included in the 1000 Genomes panel having typed HLA data. Examination of 3 indices (LD R2 values, sensitivity and specificity, minor allele frequency) revealed few SNPs with high tagging performance. All SNPs examined that tag DRB1*15:01 were in perfect LD in the British population; three showed high tagging performance in 4 of the 5 European, and 2 of the 4 American populations. For DQB1*06:02, with no previously validated tag SNPs, we show that rs3135388 has high tagging performance in one South Asian, one American, and one European population. We identify for the first time that rs2844821 has high tagging performance for A*02:01 in 5 of 7 African populations including African Americans, and 4 of the 5 European populations. These results provide a basis for selecting SNPs with high tagging performance to assess HLA alleles across diverse populations, for MS risk as well as for other diseases and conditions.