HLA Class Research Articles

Fibrinogen induces inflammatory responses via the immune activating receptor LILRA2.

The leukocyte immunoglobulin-like receptor (LILR) family, a group of primate-specific immunoreceptors, is widely expressed on most immune cells and regulates immune responses through interactions with various ligands. The inhibitory type, LILRB, has been extensively studied, and many ligands, such as HLA class I, have been identified. However, the activating type, LILRA, is less understood. We have previously identified microbially cleaved immunoglobulin as a non-self-ligand for LILRA2. In this study, we identified fibrinogen as an endogenous ligand for LILRA2 using mass spectrometry. Although human plasma contains fibrinogen in abundance in its soluble form, LILRA2 only recognizes solid-phase fibrinogen. In addition to the activating LILRA2, fibrinogen was also recognized by the inhibitory LILRB2 and by soluble LILRA3. In contrast, fibrin was recognized by LILRB2 and LILRA3, but not by LILRA2. Moreover, LILRA3 bound more strongly to fibrin than to fibrinogen and blocked the LILRB2-fibrinogen/fibrin interaction. These results suggest that morphological changes in fibrinogen determine whether activating or inhibitory immune responses are induced. Upon recognizing solid-phase fibrinogen, LILRA2 activated human primary monocytes and promoted the expression of various inflammation-related genes, such as chemokines, as revealed by RNA-seq analysis. A blocking antibody against LILRA2 inhibited the fibrinogen-induced inflammatory responses, indicating that LILRA2 is the primary receptor of fibrinogen. Taken together, our findings suggest that solid-phase fibrinogen is an inflammation-inducing endogenous ligand for LILRA2, and this interaction may represent a novel therapeutic target for inflammatory diseases.

Dysfunctional β-cell longevity in diabetes relies on energy conservation and positive epistasis.

Long-lived PFKFB3-expressing β-cells are dysfunctional partly because of prevailing glycolysis that compromises metabolic coupling of insulin secretion. Their accumulation in type 2 diabetes (T2D) appears to be related to the loss of apoptotic competency of cell fitness competition that maintains islet function by favoring constant selection of healthy "winner" cells. To investigate how PFKFB3 can disguise the competitive traits of dysfunctional "loser" β-cells, we analyzed the overlap between human β-cells with bona fide "loser signature" across diabetes pathologies using the HPAP scRNA-seq and spatial transcriptomics of PFKFB3-positive β-cells from nPOD T2D pancreata. The overlapping transcriptional profile of "loser" β-cells was represented by down-regulated ribosomal biosynthesis and genes encoding for mitochondrial respiration. PFKFB3-positive "loser" β-cells had the reduced expression of HLA class I and II genes. Gene-gene interaction analysis revealed that PFKFB3 rs1983890 can interact with the anti-apoptotic gene MAIP1 implicating positive epistasis as a mechanism for prolonged survival of "loser" β-cells in T2D. Inhibition of PFKFB3 resulted in the clearance of dysfunctional "loser" β-cells leading to restored glucose tolerance in the mouse model of T2D.

Correlation Between the Prevalence of Myasthenia Gravis and the Frequency of Class II Human Leucocyte Antigen Alleles in Various Geographical Locations Around the World.

Myasthenia gravis (MG) is an autoimmune condition characterised by muscle weakness due to antibodies produced against post-synaptic receptors. The impact of MG can be significant, especially with an ageing population. Human leukocyte antigens (HLA) are polymorphic genes associated with autoimmune conditions. Establishing the HLA alleles associated with MG may aid in the diagnosis, screening and early management of individuals at risk of MG. This research aims to establish the class II HLA alleles associated with the prevalence of MG in various regions of the world and identify the alleles that could predispose to the condition. A Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) flow chart and various databases including, Scopus and PubMed as well as other sources were used to find appropriate papers on HLA class II alleles associated with MG and the prevalence of MG in various countries. The frequency of selected HLA alleles in selected regions were obtained from the website, allelefrequencies.net. From this, a correlation coefficient and p-value were calculated to investigate whether the frequency of MG and the prevalence of HLA alleles had a significant association. The results highlighted two HLA alleles, DRB1*04:04 and DRB1*03, to have a significant positive association with the prevalence of MG. The frequency of the alleles showed regional variation, with European countries, particularly Northern Europe, exhibiting the highest frequencies. A significant positive correlation between HLA-DRB1*04:04 and DRB1*03 showed with the prevalence of MG, highlighting these alleles as a possible cause of the disease. Screening for these alleles, particularly in Northern Europe, may help identify individuals susceptible to MG.

HLA class II neoantigen presentation for CD4+ T cell surveillance in HLA class II-negative colorectal cancer

ABSTRACT Neoantigen-reactive CD4+ T cells play a key role in the anti-tumor immune response. However, the majority of epithelial tumors are negative for HLA class II (HLA-II) surface expression, and less is known about the processing of HLA-II antigens. Here, we directly identified naturally presented HLA-II neoantigens in HLA-II negative colorectal cancer (CRC) tissue using a proteogenomic approach. The neoantigens were immunogenic and induced patient CD4+ T cells with a Th1-like memory phenotype that produced IFN-γ, IL2 and TNF-α. Multiplex immunohistochemistry (IHC) demonstrated an interaction between Th cells and HLA-II-positive antigen-presenting cells (APCs) at the invasive margin and within the tertiary lymphoid structures (TLS). In our CRC cohort, the density of stromal APCs was associated with HLA-II antigen presentation in the tumor microenvironment (TME), and the number of TLS was positively correlated with the number of somatic mutations in the tumors. These results demonstrate the presence of neoantigen-specific CD4+ surveillance in HLA-II-negative CRC and suggest a potential role for macrophages and dendritic cells (DCs) at the invasive margin and in TLS for antigen presentation. Stromal APCs in the TME can potentially be used as a source for HLA-II neoantigen identification.

Frequencies of Human Leukocyte Antigen Alleles in Turkish Children with Kawasaki Disease

Abstract Objective Kawasaki disease (KD) is an acute systemic vasculitis that is one of the major causes of acquired heart disease especially in young children. The pathogenesis of KD is still unclear. The increased incidence of the disease among Japanese children and siblings of affected patients suggests a genetic component to KD susceptibility. Several reports have studied human leukocyte antigen (HLA) polymorphisms in different populations with KD and found various results. In the present study, we aimed to evaluate the association of HLA-A, -B, -C, -DRB1, and -DQB1 allele frequencies in Turkish children with KD. Methods The study was conducted between January 2016 and February 2018. HLA Class I (A, B, and C) and Class II (DRB1 and DQB1) alleles of patients and healthy controls were studied using the low-resolution DNA-based sequence-specific oligonucleotide method. Results Fifty children with KD and 500 healthy controls were included in this study. In the analysis of HLA-A, HLA-B, HLA-C, and HLA-DRB1 alleles, no statistical difference was found in the frequency of alleles between the patients with KD and the control group. However, a significantly lower frequency of the HLA-DQB1*03 allele was observed in the KD group than in the control group (p < 0.001, odds ratio [OR]: 0.29, 95% confidence interval [CI]: 0.15–0.55). When the patients with KD were divided into two subgroups with or without coronary artery lesions (CALs), the frequency of the HLA-DQB1*03 allele was also found lower in the KD group with CALs than the KD group without CALs (p = 0.008, OR: 0.19, 95% CI: 0.05–0.68). Conclusion The study may guide future studies on HLA-DQB1*03 whether it is a protective allele for KD and CALs in Turkish children.

HLA-C Peptide Repertoires as Predictors of Clinical Response during Early SARS-CoV-2 Infection.

The human leukocyte antigen (HLA) system plays a pivotal role in the immune response to viral infections, mediating the presentation of viral peptides to T cells and influencing both the strength and specificity of the host immune response. Variations in HLA genotypes across individuals lead to differences in susceptibility to viral infection and severity of illness. This study uses observations from the early phase of the COVID-19 pandemic to explore how specific HLA class I molecules affect clinical responses to SARS-CoV-2 infection. By analyzing paired high-resolution HLA types and viral genomic sequences from 60 patients, we assess the relationship between predicted HLA class I peptide binding repertoires and infection severity as measured by the sequential organ failure assessment score. This approach leverages functional convergence across HLA-C alleles to identify relationships that may otherwise be inaccessible due to allelic diversity and limitations in sample size. Surprisingly, our findings show that severely symptomatic infection in this cohort is associated with disproportionately abundant binding of SARS-CoV-2 structural and non-structural protein epitopes by patient HLA-C molecules. In addition, the extent of overlap between a given patient's predicted HLA-C and HLA-A peptide binding repertoires correlates with worse prognoses in this cohort. The findings highlight immunologic mechanisms linking HLA-C molecules with the human response to viral pathogens that warrant further investigation.

Variation within the non-coding genome influences genetic and epigenetic regulation of the human leukocyte antigen genes.

Variation within the non-coding genome may influence the regulation and expression of important genes involved in immune control such as the human leukocyte antigen (HLA) system. Class I and Class II HLA molecules are essential for peptide presentation which is required for T lymphocyte activation. Single nucleotide polymorphisms within non-coding regions of HLA Class I and Class II genes may influence the expression of these genes by affecting the binding of transcription factors and chromatin modeling molecules. Furthermore, an interplay between genetic and epigenetic factors may also influence HLA expression. Epigenetic factors such as DNA methylation and non-coding RNA, regulate gene expression without changing the DNA sequence. However, genetic variation may promote or allow genes to escape regulation by epigenetic factors, resulting in altered expression. The HLA system is central to most diseases, therefore, understanding the role of genetics and epigenetics on HLA regulation will tremendously impact healthcare. The knowledge gained from these studies may lead to novel and cost-effective diagnostic approaches and therapeutic interventions. This review discusses the role of non-coding variants on HLA regulation. Furthermore, we discuss the interplay between genetic and epigenetic factors on the regulation of HLA by evaluating literature based on polymorphisms within DNA methylation and miRNA regulatory sites within class I and Class II HLA genes. We also provide insight into the importance of the HLA non-coding genome on disease, discuss ethnic-specific differences across the HLA region and provide guidelines for future HLA studies.

HLA A*24:02-restricted T cell receptors cross-recognize bacterial and preproinsulin peptides in type 1 diabetes.

CD8+ T cells destroy insulin-producing pancreatic β cells in type 1 diabetes through HLA class I-restricted presentation of self-antigens. Combinatorial peptide library screening was used to produce a preferred peptide recognition landscape for a patient-derived T cell receptor (TCR) that recognized the preproinsulin-derived (PPI-derived) peptide sequence LWMRLLPLL in the context of disease risk allele HLA A*24:02. Data were used to generate a strong superagonist peptide, enabling production of an autoimmune HLA A*24:02-peptide-TCR structure by crystal seeding. TCR binding to the PPI epitope was strongly focused on peptide residues Arg4 and Leu5, with more flexibility at other positions, allowing the TCR to strongly engage many peptides derived from pathogenic bacteria. We confirmed an epitope from Klebsiella that was recognized by PPI-reactive T cells from 3 of 3 HLA A*24:02+ patients. Remarkably, the same epitope selected T cells from 7 of 8 HLA A*24+ healthy donors that cross-reacted with PPI, leading to recognition and killing of HLA A*24:02+ cells expressing PPI. These data provide a mechanism by which molecular mimicry between pathogen and self-antigens could have resulted in the breaking of self-tolerance to initiate disease.

Differential modulation of mutant CALR and JAK2 V617F-driven oncogenesis by HLA genotype in myeloproliferative neoplasms.

It has been demonstrated previously that human leukocyte antigen class I (HLA-I) and class II (HLA-II) alleles may modulate JAK2 V617F and CALR mutation (CALRmut)-associated oncogenesis in myeloproliferative neoplasms (MPNs). However, the role of immunogenetic factors in MPNs remains underexplored. We aimed to investigate the potential involvement of HLA genes in CALRmut+ MPNs. High-resolution genotyping of HLA-I and -II loci was conducted in 42 CALRmut+ and 158 JAK2 V617F+ MPN patients and 1,083 healthy controls. A global analysis of the diversity of HLA-I genotypes revealed no significant differences between CALRmut+ patients and controls. However, one HLA-I allele (C*06:02) showed an inverse correlation with presence of CALR mutation. A meta-analysis across independent cohorts and healthy individuals from the 1000 Genomes Project confirmed an inverse correlation between the presentation capabilities of the HLA-I loci for JAK2 V617F and CALRmut-derived peptides in both patients and healthy individuals. scRNA-Seq analysis revealed low expression of TAP1 and CIITA genes in CALRmut+ hematopoietic stem and progenitor cells. In conclusion, the HLA-I genotype differentially restricts JAK2 V617F and CALRmut-driven oncogenesis potentially explaining the mutual exclusivity of the two mutations and differences in their presentation latency. These findings have practical implications for the development of neoantigen-based vaccines in MPNs.

Proteasome inhibition enhances the anti-leukemic efficacy of chimeric antigen receptor (CAR) expressing NK cells against acute myeloid leukemia

BackgroundRelapsed and refractory acute myeloid leukemia (AML) carries a dismal prognosis. CAR T cells have shown limited efficacy in AML, partially due to dysfunctional autologous T cells and the extended time for generation of patient specific CAR T cells. Allogeneic NK cell therapy is a promising alternative, but strategies to enhance efficacy and persistence may be necessary. Proteasome inhibitors (PI) induce changes in the surface proteome which may render malignant cells more vulnerable to NK mediated cytotoxicity. Here, we investigated the potential benefit of combining PIs with CAR-expressing allogeneic NK cells against AML.MethodsWe established the IC50 concentrations for Bortezomib and Carfilzomib against several AML cell lines. Surface expression of class-I HLA molecules and stress-associated proteins upon treatment with proteasome inhibitors was determined by multiparameter flow cytometry. Using functional in vitro assays, we explored the therapeutic synergy between pre-treatment with PIs and the anti-leukemic efficacy of NK cells with or without expression of AML-specific CAR constructs against AML cell lines and primary patient samples. Also, we investigated the tolerability and efficacy of a single PI application strategy followed by (CAR-) NK cell infusion in two different murine xenograft models of AML.ResultsAML cell lines and primary AML patient samples were susceptible to Bortezomib and Carfilzomib mediated cytotoxicity. Conditioned resistance to Azacitidine/Venetoclax did not confer primary resistance to PIs. Treating AML cells with PIs reduced the surface expression of class-I HLA molecules on AML cells in a time-and-dose dependent manner. Stress-associated proteins were upregulated on the transcriptional level and on the cell surface. NK cell mediated killing of AML cells was enhanced in a synergistic manner. PI pre-treatment increased effector-target cell conjugate formation and Interferon-γ secretion, resulting in enhanced NK cell activity against AML cell lines and primary samples in vitro. Expression of CD33- and CD70-specific CARs further improved the antileukemic efficacy. In vivo, Bortezomib pre-treatment followed by CAR-NK cell infusion reduced AML growth, leading to prolonged overall survival.ConclusionsPIs enhance the anti-leukemic efficacy of CAR-expressing allogeneic NK cells against AML in vitro and in vivo, warranting further exploration of this combinatorial treatment within early phase clinical trials.

Hiding in plain sight: Misinterpretation of immunogenic DPB epitopes within G/P groups

The clinical impact of HLA DP antibodies is poorly understood, resulting in variable clinical strategies for transplant candidates and recipients with donor-directed HLA-DP antibodies. Complicating matters further, the DPB naming convention is not based on allelic homology and requires sequence alignments to identify potential immunogenic epitopes. Historically, G and P codes, which consolidated alleles that were identical over Exon 2, were used to simplify the reporting of HLA Class II typing as differences outside of Exon 2 have not been considered immunogenic (i.e., able to induce an antibody response). Herein, we present four cases demonstrating that polymorphisms at codons 96R/K and 170I/T, in Exon 3 of DPB, are targets for alloantibody recognition. These regions “hide in plain sight” due to the current use of G/P code-level typing, potentially leading to incorrect compatibility assessments (i.e., virtual crossmatches) and misinterpreted antibody responses. The unintentional crossing of an HLA-DPB donor-specific antibody (DSA) in a solid organ or hematopoietic stem cell transplant may lead to unforeseen deleterious clinical outcomes. Our data underscore the complexities of DPB histocompatibility assessments and highlight the need for adaptable systems that align with evolving research and clinical outcomes.

Non-coding mismatches in HLA Class I genes in the presence of DPB1 mismatches predict worse outcomes in hematopoietic cell transplantation

Non-coding mismatches in HLA Class I genes in the presence of DPB1 mismatches predict worse outcomes in hematopoietic cell transplantation

Risk of kidney graft failure (GF) associated with HLA class II amino acid mismatches (AAMM) is highest during the first year after transplant for sensitized patients

Risk of kidney graft failure (GF) associated with HLA class II amino acid mismatches (AAMM) is highest during the first year after transplant for sensitized patients

Validation of Alloseq Tx11 – Next Generation Sequencing of 11 classical HLA loci

Validation of Alloseq Tx11 – Next Generation Sequencing of 11 classical HLA loci

Immunogenetic profiles of 9 human herpes virus envelope glycoproteins

Human herpes viruses (HHV) are ubiquitous and have been implicated in numerous long-term health conditions. Since the association between viral exposure and long-term health impacts is partially influenced by variation in human leukocyte antigen (HLA) genes, we evaluated in silico the binding affinities of 9 HHV envelope glycoproteins with 127 common HLA Class I and Class II molecules. The findings show substantial variability in HHV binding affinity across viruses, HLA Class, HLA genes, and HLA alleles. Specific findings were as follows: (1) the predicted binding affinities of HHVs were characterized by four distinct groupings—[HHV1, HHV2], [HHV3, HHV4, HHV5], [HHV6A], [HHV6B, HHV7, HHV8]—with relatively lower binding affinities for HHV1, HHV2, and HHV6a compared to other HHVs; (2) significantly higher binding affinity was found for HLA Class I relative to Class II; (3) analyses within each class demonstrated that alleles of the C gene (for Class I) and DRB1 gene (for Class II) had the highest binding affinities; and (4) for each virus, predicted binding affinity to specific alleles varied, with HHV6a having the lowest affinity for HHV-HLA complexes, and HHV3, HHV4, and HHV5 having the highest. Since HLA-antigen binding is the first step in initiating an immune response to foreign antigens, these relative differences in HHV binding affinities are likely to influence long-term health impacts such that the cells infected with viruses associated with higher binding affinities across common HLA alleles may be more reduced in numbers, thereby lowering the potential for long-term sequelae of their infections.

Disruptions in antigen processing and presentation machinery on sarcoma

BackgroundThe antigen processing machinery (APM) plays a critical role in generating tumor-specific antigens that can be recognized and targeted by the immune system. Proper functioning of APM components is essential for presenting these antigens on the surface of tumor cells, enabling immune detection and destruction. In many cancers, defects in APM can lead to immune evasion, contributing to tumor progression and poor clinical outcomes. However, the status of the APM in sarcomas is not well characterized, limiting the development of effective immunotherapeutic strategies for these patients.MethodsWe investigated 126 patients with 8 types of bone and soft tissue sarcoma operated between 2001–2021. Tissue microarrays mapped 11 specific areas in each case. The presence/absence of APM protein was determined through immunohistochemistry. Bayesian networks were used.ResultsAll investigated sarcomas had some defects in APM. The least damaged component was HLA Class I subunit β2-microglobulin and HLA Class II. The proteasome LMP10 subunit was defective in leiomyosarcoma (LMS), myxoid liposarcoma (MLPS), and dedifferentiated liposarcoma (DDLPS), while MHC I transporting unit TAP2 was altered in undifferentiated pleomorphic sarcoma (UPS), gastrointestinal stromal tumor (GIST), and chordoma (CH). Among different neoplastic areas, high-grade areas showed different patterns of expression compared to high lymphocytic infiltrate areas. Heterogeneity at the patient level was also observed. Loss of any APM component was prognostic of distant metastasis (DM) for LMS and DDLPS and of overall survival (OS) for LMS.ConclusionSarcomas exhibit a high degree of defects in APM components, with differences among histotypes and tumoral areas. The most commonly altered APM components were HLA Class I subunit β2-microglobulin, HLA Class I subunit α (HC10), and MHC I transporting unit TAP2. The loss of APM components was prognostic of DM and OS and clinically relevant for LMS and DDLPS. This study explores sarcoma molecular mechanisms, enriching personalized therapeutic approaches.

SVA Regulation of Transposable Element Clustered Transcription within the Major Histocompatibility Complex Genomic Class II Region of the Parkinson's Progression Markers Initiative.

SINE-VNTR-Alu (SVA) retrotransposons can regulate expression quantitative trait loci (eQTL) of coding and noncoding genes including transposable elements (TEs) distributed throughout the human genome. Previously, we reported that expressed SVAs and human leucocyte antigen (HLA) class II genotypes on chromosome 6 were associated significantly with Parkinson's disease (PD). Here, our aim was to follow-up our previous study and evaluate the SVA associations and their regulatory effects on the transcription of TEs within the HLA class II genomic region. We reanalyzed the transcriptome data of peripheral blood cells from the Parkinson's Progression Markers Initiative (PPMI) for 1530 subjects for TE and gene RNAs with publicly available computing packages. Four structurally polymorphic SVAs regulate the transcription of 20 distinct clusters of 235 TE loci represented by LINES (37%), SINES (28%), LTR/ERVs (23%), and ancient transposon DNA elements (12%) that are located in close proximity to HLA genes. The transcribed TEs were mostly short length, with an average size of 389 nucleotides. The numbers, types and profiles of positive and negative regulation of TE transcription varied markedly between the four regulatory SVAs. The expressed SVA and TE RNAs in blood cells appear to be enhancer-like elements that are coordinated differentially in the regulation of HLA class II genes. Future work on the mechanisms underlying their regulation and potential impact is essential for elucidating their roles in normal cellular processes and disease pathogenesis.

HLA class I expression on human platelets is highly variable and correlates with distinct allele group frequencies.

Human leukocyte antigen (HLA) class I molecules are expressed on platelets and can represent a source of alloimmunization in recipients of platelet transfusions. HLA mismatch between donors and recipients may be associated with the induction of anti-HLA antibodies, which can culminate in refractoriness to platelet transfusions. In the present study we analyzed HLA allele group frequencies and HLA expression levels on human platelets from blood donors. Platelet-rich plasma was collected from 139 donors to monitor platelet HLA class I expression by flow cytometry. DNA from donors with high and low platelet HLA expression was used in the genotype studies. Frequencies of large and normal-sized platelet subpopulations were determined and HLA class I expression was studied. Mean platelet volume (MPV) and platelet large-cell ratio (P-LCR) were analyzed in both groups of donors. The analysis showed variable platelet HLA class I expression with significant differences among donors. HLA class I allele group frequencies in donors with high and low platelet HLA expression showed distinctive genotypic features strictly related to expression level. The main allele groups found in samples with high platelet HLA class I expression were HLA-A*02, -A*68, -B*15, -B*49, and -C*03. Platelet HLA class I expression did not change over time or during freezing-thawing cycles. The analysis of platelet subpopulations showed a statistically significant higher expression of HLA class I molecules on large platelets than on normal-sized platelets. Moreover, donors with high HLA class I expression showed a higher frequency of large platelets (p<0.0001). The analysis of P-LCR in both groups of donors showed a statistically significant difference (p<0.05) within high HLA-expressing donors. Our data suggest an allele-dependent expression of HLA class I molecules on human platelets with distinct HLA allele group frequencies and different platelet subpopulation frequencies among blood donors.

Identification of the novel HLA-C*08:291 allele by next-generation sequencing.

HLA-C*08:291, a novel HLA class I allele detected by next-generation sequencing.

HLA class II profile in patients with different stages of cystic echinococcosis according to the WHO ultrasound imaging classification.

The factors involving in the natural history and determinants of different features of human cystic echinococcosis (CE) are not adequately understood. Several host-related factors including the genetic structure of the host and human leukocyte antigens (HLAs) are believed to be involved in the natural history of CE in humans. The present study was conducted to investigate the association between HLA class II genes and active and inactive stages of hepatic cystic echinococcosis. Echinococcus granulosus cyst samples and patient information were collected from the biobank of the Iranian Hydatid Disease Registry from 2019 to 2022. HLA-DRB and HLA-DQB were characterized by PCR method. CE patients were categorized into three active (CE1 and CE2), inactive (CE4 and CE5), and transitional (CE3) stages according to the WHO ultrasound classification of CE. In total, 77 participants including 38 patients (36.8% men and 63.2% women) with different stages of CE as well as 39 healthy individuals (38.5% men and 61.5% women) were included in the study. Findings of the study showed that the frequency of HLA-DRB1*03 was significantly lower in the patients compared to the healthy individuals. The frequencies of HLA-DQB and HLA-DRB alleles were not differed significantly between active, inactive, and transitional stages of E. granulosus cysts. Findings of this study indicate the potential role of this allele in the susceptibility of human to cystic echinococcosis. Further large-scale studies in different endemic countries are required to document the significance of HLA-DQB and HLA-DRB as a host-related factor in the natural history of CE in human.