- 1Department of Blood Transfusion, Ningbo Medical Center Lihuili Hospital, The Affiliated Lihuili Hospital of Ningbo University, Ningbo, China
- 2HLA Typing Laboratory, Blood Center of Zhejiang Province, Hangzhou, China
- 3Department of Radiotherapy, Ningbo Medical Center Lihuili Hospital, The Affiliated Lihuili Hospital of Ningbo University, Ningbo, China
- 4Department of Hepatopancreatobiliary Surgery, Ningbo Medical Center Lihuili Hospital, The Affiliated Lihuili Hospital of Ningbo University, Ningbo, China
Background: It has been reported that HLA was associated with susceptibility to various liver diseases. However, data on HLA in end-stage liver disease (ESLD) patients are limited. Therefore, this study aimed to explore the association of HLA with susceptibility to ESLD in patients undergoing liver transplantation (LT).
Methods: 108 ESLD patients who underwent LT from July 1, 2023, to December 31, 2024, were investigated. At the same time, 453 blood donors were randomly selected as healthy controls. Their specimens were genotyped for 11 HLA loci using next-generation sequencing (NGS).
Results: There was no significant difference in gender and ABO blood group between the two groups (P>0.05) except that ESLD patients were older (P < 0.001). The allele frequencies (AFs) of HLA-B*18:02, -C*07:02, -DRB1*12:02, -DRB1*14:54, -DQA1*01:04, -DQA1*06:01, -DPA1*02:02, -DPB1*05:01:01G and -DPB1*135:01 in ESLD were significantly higher than those of the control group (P < 0.05), while the AFs of HLA-DRB1*15:01, -DQA1*01:01, -DQA1*01:02, -DQB1*06:02, -DPA1*01:03 and -DPB1*02:01:02G were significantly lower than that of the control group (P < 0.05). However, after correction by the step-down Bonferroni method, only HLA-DQA1*06:01 and HLA-DPB1*05:01:01G were considered statistically significant (Pc < 0.05). The association between these two HLA alleles and ESLD remained significant (P < 0.05) after adjustment for age and gender by binary logistic regression analysis. Subgroup analysis further revealed that this association was primarily linked to hepatitis B virus (HBV) infection. Compared with the control group, there were statistically significant differences in 3 haplotypes with frequencies >1.00% in ESLD (P < 0.05). However, none were statistically significant after correction (Pc>0.05).
Conclusion: Our study suggested that HLA-DQA1*06:01 and HLA-DPB1*05:01:01G alleles may be associated with susceptibility to ESLD in liver transplant recipients, particularly in cases related to HBV infection.
1 Introduction
Liver disease causes two million deaths annually worldwide and accounts for 4% of all deaths (1). The concept of end-stage liver disease (ESLD) was first introduced in the 1980s, which refers to the end stage of chronic liver disease, with advanced liver injury, dysfunction, and decompensation, regardless of the etiology (2, 3). The main clinical feature of ESLD is the liver’s inability to meet the body’s physiological needs. The range of clinical disease forms of ESLD includes acute-on-chronic liver failure (ACLF), acute decompensated cirrhosis, chronic liver failure (CLF), and decompensated hepatocellular carcinoma (HCC) (4, 5). At present, liver transplantation (LT) is still the only way to cure ESLD (6). The model for end-stage liver disease (MELD) score is frequently used to assess the need for LT (7). Although LT is performed for cure ESLD in more than 100 countries worldwide, the majority of the patients with ESLD and life-threatening liver diseases are still unable to obtain a liver for transplantation (7, 8). Therefore, liver disease, especially ESLD, remains a considerable challenge.
The human leukocyte antigen (HLA) region is one of the most polymorphic genes in the human genome (9). HLA distribution differs significantly in various populations and ethnic groups (10, 11). With the rapid development of precise HLA genotyping technology, 42,193 HLA alleles have been identified in the latest IPD-IMGT/HLA database release 3.61version (2025–07) (12). HLA plays a critical role in human immunity. Class I molecules present endogenous peptides to CD8+ T cells, triggering cytotoxicity against infected cells. In contrast, class II molecules present exogenous antigens to activate T helper lymphocytes (13). In recent years, HLA has been reported to be associated with a variety of diseases (14–20), including viral hepatitis (14, 15), liver cancer (16), liver metabolic diseases (17, 18). However, the data for association of HLA with ESLD is still rare. Here, the association between HLA and ESLD in the Zhejiang Han population, China was studied.
2 Materials and methods
2.1 Study design and population
All patients with ESLD who underwent LT at Ningbo Medical Center Lihuili Hospital from July 1, 2023, to December 31, 2024 were enrolled. Inclusion criteria for the patients: (1)age ≥ 18 years; (2) diagnosis of ESLD requiring LT. In this study, the diagnosis of ESLD was primarily established by identifying at least one of the following conditions in the electronic medical records (EMR) admitting diagnoses: ACLF, decompensated cirrhosis, CLF, or decompensated HCC. These cases were subsequently confirmed as ESLD through a manual review conducted by our expert LT team. Exclusion criteria: (1) non-Zhejiang Han population; (2) a history of primary malignant tumors other than liver cancer; (3) autoimmune diseases other than autoimmune liver disease; (4) incomplete clinical data. A total of 108 patients with ESLD were finally enrolled in this study. At the same time, 453 healthy blood donors from the Zhejiang Han population in the same period were randomly selected as the control group. This study was approved by the ethics committee of Ningbo Medical Center Lihuili Hospital (Approval No: 2024 - 009). Written informed consent was obtained from all participants.
2.2 Data and specimen collection
Demographic and characteristics of ESLD patients and blood donors, including age, gender, ABO blood group were obtained through EMR or blood donation information systems, respectively. The remaining anticoagulant specimens for ABO blood group testing of ESLD patients and blood donors were collected respectively and stored below -20 °C until processing.
2.3 HLA genotyping
All specimens were genotyped for HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DRB3, HLA-DRB4, HLA-DRB5, HLA-DQA1, HLA-DQB1, HLA-DPA1 and HLA-DPB1 using commerical AllType™NGS kit (One Lambda Inc, Canoga Park, CA, USA) and HLA genotype assignment was according to our previous reports (21–23).
2.4 Statistical analysis
Categorical variables were summarized as counts and percentages, with group comparisons performed using the chi-square test or Fisher’s exact test. Continuous variables exhibiting non-normal distributions were expressed as medians (interquartile ranges), and analyzed with the Mann-Whitney U test. The allele frequencies (AFs) and haplotype frequencies (HFs) of HLA loci, as well as Hardy–Weinberg equilibrium for each HLA locus, were determined by the Arlequin software 3.5.2.2 (24). Continuous variables such as age were dichotomized at the median. A binary logistic regression analysis model was used to predict ESLD, adjusted for specific HLA alleles, age, and gender. The P values of AFs and HFs were corrected by the step-down Bonferroni method (Pc values). A P value or Pc value less than 0.05 was considered statistically significant. Statistical calculation was performed by the SPSS statistical software, version 24.0 (IBM, Armonk, NY, USA).
3 Results
3.1 High hepatitis B virus infection rate was found in the ESLD patients
According to the EMR data, the ESLD patient group had a MELD score of 22.0 (19.0-24.0). Among them, 58 (53.70%) and 50 (46.30%) patients were graded as Child-Pugh B and C, respectively. Additionally, 77 cases (71.30%) were co-infected with the hepatitis B virus (HBV), which was also the most common underlying etiology of ESLD. There was no significant difference in gender and ABO blood group between the ESLD patient group and control group (P>0.05) except that ESLD patients were older (median age 52.5 years vs. 30.0 years) (P < 0.001). The demographic and ABO blood group distribution characteristics of the ESLD patients and blood donors were shown in Table 1.
Table 1. Basic characteristics of the ESLD patients and healthy controls in the Zhejiang Han population*.
3.2 AFs distribution of the 11 HLA loci in the ESLD patients and healthy controls
The allele distribution of the 11 HLA loci was fitted Hardy-Weinberg equilibrium (P>0.05) expect for HLA-DPA1 locus in the ESLD patients (Supplementary Table S1). However the Pc value was more than 0.05 for HLA-DPA1 locus. The frequencies of HLA alleles are summarized in Supplementary Table S2. The most frequent alleles of each locus in ESLD patients and healthy controls were HLA-A*11:01 (25.46% vs 22.74%), B*40:01(15.74% vs 12.36%), C*07:02(22.22% vs 15.56%), DRB1*12:02(15.74%) vs DRB1*09:01 (14.90%), DRB3*02:02(57.61% vs 56.73%), DRB4*01:03 (89.61% vs 91.82%), DRB5*01:01(80.00% vs 78.43%), DQA1*03:02 (16.20%) vs DQA1*01:02(17.66%), DQB1*03:01(26.39% vs 21.74%), DPA1*02:02 (57.87% vs 47.57%), DPB1*05:01:01G (49.07% vs 36.31%), respectively.
3.3 DQA1*06:01 and DPB1*05:01:01G may be associated with the susceptibility of ESLD
15 HLA alleles of the HLA loci were showed uncorrected significant differences between the ESLD patient group and control group (Table 2). The AFs of HLA-B*18:02, -C*07:02, -DRB1*12:02, -DRB1*14:54, -DQA1*01:04, -DQA1*06:01, -DPA1*02:02, -DPB1*05:01:01G and -DPB1*135:01 in ESLD patient group were significantly higher than those of the control group (P < 0.05), while the AFs of HLA-DRB1*15:01, -DQA1*01:01, -DQA1*01:02, -DQB1*06:02, -DPA1*01:03 and -DPB1*02:01:02G were significantly lower than those of the control group (P < 0.05). However, only HLA-DQA1*06:01 and HLA-DPB1*05:01:01G were considered statistically significant (Pc < 0.05) after correction (Figure 1). Further it found that these associations were related to HBV infection in the subgroup analysis (Figure 2 and Supplementary Table S3-Supplementary Table S8). However, in this subgroup of HBV-infected, no significant differences were observed in the frequencies of HLA-DQA1*06:01 and HLA-DPB1*05:01:01G between those with positive and negative in the pre-operative HBV DNA detection (P>0.05; Supplementary Table S9).
Table 2. HLA alleles with significant differences between ESLD patients and healthy controls in the Zhejiang Han population.
Figure 2. Pairwise comparisons of DQA1*06:01, DPB1*05:01:01G, and DRB5 blank were performed across the HBV-related ESLD group, the non-HBV-related ESLD group, and the control group. *Pc < 0.05 corrected by the step-down Bonferroni method.
3.4 HLA-DRB5 blank may be associated with the susceptibility of ESLD
HLA-DRB5 blank frequency was higher in the ESLD patients than that in the control group (90.74% vs. 83.11%, P = 0.005, Pc = 0.035) (Supplementary Table S10). Similarly, The proportion of two copy of HLA-DRB5 absence in the ESLD patients was higher than that in the control group (82.4% vs. 68.9%, P = 0.005). However, this phenomenon was not found in HLA-DRB3 or HLA-DRB4 loci. Further subgroup analysis found that the association was related to HBV infection (Figure 2 and Supplementary Table S11-Supplementary Table S13). However, no significant differences of HLA-DRB5 blank frequencies were observed between those with positive and negative in the pre-operative HBV DNA detection (P>0.05; Supplementary Table S9).
3.5 Associations persisted after adjustment for age and gender
In the binary logistic regression analysis adjusted for age and gender, HLA-DQA1*06:01, -DPB1*05:01:01G, and the two copy of HLA-DRB5 absence remained independent risk factors for ESLD (P < 0.05; Figure 3). Additionally, among the healthy controls, no significant difference was observed in the frequencies of HLA-DQA1*06:01, -DPB1*05:01:01G, and -DRB5 blank between the age >30 and age ≤30 subgroups (Pc>0.05; Supplementary Table S14).
Figure 3. Binary logistic regression for predicting ESLD, adjusted for specific HLA alleles, age, and gender. (A) adjusted for HLA-DQA1*06:01; (B) adjusted for HLA-DPB1*05:01:01G; (C) adjusted for the absence of HLA-DRB5.
3.6 DQA1*06:01, DPB1*05:01:01G, and DRB5 blank may not be associated with ESLD severity
In patients with ESLD, the frequencies of HLA-DQA1*06:01, -DPB1*05:01:01G, and -DRB5 blank did not differ significantly between groups with MELD scores >22 and ≤22 (P>0.05; Table 3). Similarly, no significant differences in their frequencies were observed between patients with Child-Pugh grade B and grade C (P>0.05; Table 4).
Table 3. Comparison of HLA AFs between MELD score groups (>22 vs ≤22) in ESLD patients from the Zhejiang Han population.
Table 4. Comparison of HLA AFs across Child-Pugh grades B and C in ESLD patients from the Zhejiang Han population.
3.7 No HLA haplotypes association with ESLD patients
A total of 182 and 651 different HLA haplotypes were estimated in the ESLD patient group and control group. The highest frequencies in the ESLD patient group were A*11:01-C*07:02-B*40:01-DRB1*12:02-DQA1*06:01-DQB1*03:01-DPA1*02:02-DPB1*05:01:01G (2.31%) and A*30:01-C*06:02-B*13:02-DRB1*07:01-DQA1*02:01-DQB1*02:02-DPA1*02:01-DPB1*17:01:01G (2.31%), while the latter also had the highest frequency in the control group (2.76%). HLA haplotypes with frequencies >1% were listed in the Table 5. A total of 3 haplotypes were found to have uncorrected significant differences (P < 0.05), but none had substantial differences after correction (Pc>0.05).
Table 5. Comparison of the 8 HLA loci haplotypes with frequency > 1% between ESLD and healthy controls in the Zhejiang Han population.
4 Discussion
In the healthy controls of this study, the frequencies of HLA-DQA1*06:01 and -DPB1*05:01:01G were observed to be 8.72% and 36.31%, respectively, which are consistent with the previously reported frequencies of 7.87% and 39.79% in the local population (22). Meanwhile, among the HLA-DRB5 alleles, DRB5*01:01 showed the highest frequency at 78.43%, also aligning with the previously reported 77.26% in the local population (22). However, the study did not specify the frequency of HLA-DRB5 blank (22). In our study, the frequency of HLA-DRB5 blank was 83.11%, which is close to the 87.31% reported in the neighboring country of Korea (11). Furthermore, to minimize potential population stratification bias, the Hardy–Weinberg equilibrium test was performed.
There are many studies on the association between HLA and susceptibility to liver diseases, especially viral hepatitis (25–27). However, there are limited data on the association between HLA and susceptibility to ESLD. Tillmann et al. (28) found that HLA-DRB1*11 and HLA-DQB1*03 were significantly lower in hepatitis C virus (HCV) induced ESLD patients. However, it was not found in our study. This difference may be due to the different etiologies of ESLD in the two studies, which were mainly induced by HCV and HBV, respectively, and the difference in HLA distribution in the various populations.
The main etiologies of ESLD in north America are alcohol-related liver disease, non-alcoholic steatohepatitis (NASH), and chronic HCV (29). However, it is mainly associated with chronic HBV infection in Asian populations (30). Most of liver transplant recipients had chronic HBV infection in this study. It was found that HLA-DQA1*06:01 and HLA-DPB1*05:01:01G may be associated with susceptibility in ESLD in this study. However, we did not find any association between these HLA alleles and the severity of ESLD. Studies (27, 31) have reported that HLA-DQA1*06:01 may be a predisposing factor for HBV infection and progression to cirrhosis. Additional studies (25, 32) have reported HLA-DPB1*05:01 as a risk factor for HBV infection. Since about 70% of ESLD in our study were also chronically infected with HBV, therefore HLA-DQA1*06:01 and HLA-DPB1*05:01:01G as risk factors for ESLD may be original for HBV infection, which was validated in our subgroup analysis. In the HBV-infected subgroup, we further assessed the impact of the patients’ infection status by comparing those with detectable and undetectable pre-operative HBV DNA levels. However, no significant differences in HLA alleles were observed between these two groups. Furthermore, given that the vast majority of HBV-infected ESLD patients had received antiviral therapy prior to surgery (70/77) and had comorbid cirrhosis (74/77), these factors were not considered as potential sources of bias in the comparative analysis.
However, the mechanisms underlying these associations are not yet fully understood. HLA class II molecules play a critical role in shaping CD4+ T lymphocyte responses via antigen presentation. Notably, the HLA-DQ molecule is distinct from other class II antigens because most of its variable amino acid residues reside within the α-helix region of the antigen-binding groove (33). Thus, we hypothesize that the association of HLA-DQB1*06:01 with chronic HBV infection could stem from its reduced capacity to efficiently bind and present HBV antigens. Nevertheless, this hypothesis requires validation through further investigation. In addition, a proposed mechanism for the association of HLA-DPB1*05:01:01G with HBV infection involves its reduced mRNA expression, which could impair antigen presentation (25). This impairment may subsequently diminish the immune response against HBV, thereby increasing susceptibility to persistent infection. An additional theory suggests that HLA-DPB1*05:01 exhibits suboptimal binding affinity for specific peptide epitopes derived from HBV nucleocapsid proteins, leading to a reduced likelihood of viral clearance (34).
HLA-DQA1*01:02 has been suggested to be a protective factor for HBV infection (31, 35), but it was not found in this study. Previous studies (25, 26, 32) have found that HLA-DPB1*02:01 is a protective factor for HBV infection, and it was also observed in our study of the ESLD patients, but it was not statistically significant after adjustment. HLA-DPB1*09:01 has been reported (25, 26) as a risk factor for HBV infection but was not found in ESLD in this study. The reason need to further explore and one of them may be due to about 30% of ESLD patients with not HBV-infected in this study.
Interestingly, HLA-DRB5 blank was found as risk factors for ESLD patients. It has been reported that a lack of HLA-DRB5 may be associated with various diseases, including rheumatoid arthritis patients with negative antibodies against citrullinated peptides (36) and patients with multiple sclerosis (37). Therefore, the absence of HLA-DRB5 may require more attention and further research in the future.
HLA AFs may differ across age groups (38). In our study, blood donors were used as controls. While this selection may result in a “healthier” control group, the significant age difference between the younger blood donors and the older ESLD patients raised the concern that the observed differences in the frequencies of HLA-DQA1*06:01, -DPB1*05:01:01G, and -DRB5 blank might be attributable to age variation rather than disease association. To address this, we stratified the control group by median age into younger and older subgroups. However, no significant differences in the frequencies of these key HLA alleles were observed between the two subgroups. Most importantly, we performed logistic regression adjusted for age and gender as confounding factors. The results demonstrated that the associations between these HLA alleles and ESLD remained significant, thereby strengthening the reliability of our findings.
We observed significantly elevated frequencies of HLA-DQA1*06:01, -DPB1*05:01:01G, and -DRB5 blank in patients with ESLD. This finding not only carries population genetic relevance but also implies potential clinical importance for immune risk management after LT. Firstly, HLA-DQ mismatch is a well-established risk factor for de novo donor-specific antibody (DSA) development, and anti-DQ antibodies are strongly associated with antibody-mediated rejection and graft loss (39). The high frequency of HLA-DQA1*06:01 in our cohort suggests that these recipients may harbor a genetic background predisposing them to heightened immune responses against mismatched donor HLA-DQ antigens. Secondly, in hematopoietic stem cell transplantation, the immunologic impact of HLA-DPB1 is clearly defined: non-permissive mismatches potently enhance alloreactive T-cell activation, increasing the incidence of acute graft-versus-host disease and mortality (40). Although evidence in LT remains limited, the underlying principles of alloimmunity are likely shared. Thus, HLA-DPB1*05:01:01G mismatch may elevate the risk of T-cell-mediated rejection and adversely affect long-term graft survival via cytotoxic T-lymphocyte pathways. Finally, the absence of DRB5 in a recipient implies that donor DRB5 will be recognized as a foreign antigen, substantially increasing the risk of developing anti-DRB5 DSAs. In solid organ transplantation, antibodies against HLA class II antigens are powerful predictors of chronic rejection (41). Therefore, implementing high-resolution matching at these loci before transplantation and evaluating the permissibility of mismatches could improve donor-recipient selection and potentially enhance post-transplant outcomes.
There are several limitations in this study. Firstly, the collection of indicators in ESLD and healthy controls, such as body mass index (BMI), was incomplete. Secondly, the sample number was relatively small, especially for ESLD. Therefore, a multicenter study with a larger sample size is warranted to validate our findings further. Thirdly, the control group was significantly younger than the study group. These differences may lead to misinterpretation of the results because HLA alleles may differ between young and old populations. Thus, it needs to be studied further in the future.
5 Conclusion
In this study, the association between AFs and HFs of 11 HLA loci at high resolution level and ESLD was firstly analyzed in the Zhejiang Han population, China. It was found that HLA-DQA1*06:01, HLA-DPB1*05:01:01G and HLA-DRB5 absence may be associated with susceptibility to ESLD in liver transplant recipients, particularly in cases related to HBV infection. These results will help further understand the role of HLA loci in the pathogenesis of ESLD requiring LT in China.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the ethics committee of Ningbo Medical Center Lihuili Hospital (Approval No: 2024 - 009). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
GL: Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft. YH: Investigation, Visualization, Writing – review & editing. YL: Investigation, Software, Writing – review & editing. WC: Data curation, Investigation, Supervision, Writing – review & editing. CL: Data curation, Investigation, Supervision, Writing – review & editing. FZ: Data curation, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by Open Foundation of Key Laboratory of Blood Safety Research of Zhejiang Province (2023KF002), Ningbo Young Technical Backbone Talents (2021QNJSGG-LGJ), Ningbo Top Medical and Health Research Program (2024020818), and The Ningbo Major Research and Development Plan Project (2024Z179).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1684437/full#supplementary-material
References
1. Devarbhavi H, Asrani SK, Arab JP, Nartey YA, Pose E, and Kamath PS. Global burden of liver disease: 2023 update. J hepatology. (2023) 79:516–37. doi: 10.1016/j.jhep.2023.03.017
2. Houssin D, Franco D, Corlette MB, and Bismuth H. Criteria for hepatic transplantation in cirrhosis. Surgery gynecology obstetrics. (1980) 151:30–2.
3. Chen T, Chen G, Wang G, Treeprasertsuk S, Lesmana CRA, Lin HC, et al. Expert consensus on the diagnosis and treatment of end-stage liver disease complicated by infections. Hepatol Int. (2024) 18:817–32. doi: 10.1007/s12072-023-10637-3
4. Xiang Z, Li J, Zeng H, Xiang X, Gao F, Wang K, et al. Current understanding of marginal grafts in liver transplantation. Aging disease. (2024) 16:1036–58. doi: 10.14336/AD.2024.0214
5. Xiao J, Wang F, Wong NK, He J, Zhang R, Sun R, et al. Global liver disease burdens and research trends: Analysis from a Chinese perspective. J hepatology. (2019) 71:212–21. doi: 10.1016/j.jhep.2019.03.004
6. Mugaanyi J, Tong J, Lu C, Mao S, Huang J, and Lu C. Risk factors for acute rejection in liver transplantation and its impact on the outcomes of recipients. Transplant Immunol. (2023) 76:101767. doi: 10.1016/j.trim.2022.101767
7. Lucey MR, Furuya KN, and Foley DP. Liver transplantation. New Engl J Med. (2023) 389:1888–900. doi: 10.1056/NEJMra2200923
8. Feng S, Roll GR, Rouhani FJ, and Sanchez Fueyo A. The future of liver transplantation. Hepatol (Baltimore Md). (2024) 80:674–97. doi: 10.1097/HEP.0000000000000873
9. Montgomery RA, Tatapudi VS, Leffell MS, and Zachary AA. HLA in transplantation. Nat Rev Nephrology. (2018) 14:558–70. doi: 10.1038/s41581-018-0039-x
10. Mirza A, Ali Qadri MM, Zeshan B, Hafiz K, Abbas S, Ahmad N, et al. HLA class-I polymorphisms among the Punjabi population of Pakistan: A comparative analysis with country’s other ethnic groups. Hum Immunol. (2024) 85:111083. doi: 10.1016/j.humimm.2024.111083
11. Jung K, Kim JG, Shin S, Roh EY, Hong YJ, and Song EY. Allele and haplotype frequencies of 11 HLA loci in Koreans by next-generation sequencing. Hla. (2023) 101:602–12. doi: 10.1111/tan.14980
12. Robinson J, Barker DJ, and Marsh SGE. 25 years of the IPD-IMGT/HLA database. Hla. (2024) 103:e15549. doi: 10.1111/tan.15549
13. Byrne GW. Does human leukocyte antigens sensitization matter for xenotransplantation? Xenotransplantation. (2018) 25:e12411. doi: 10.1111/xen.12411
14. Naderi M, Hosseini SM, Behnampour N, Shahramian I, and Moradi A. Association of HLADQ-B1 polymorphisms in three generations of chronic hepatitis B patients. Virus Res. (2023) 325:199036. doi: 10.1016/j.virusres.2022.199036
15. Huang J, Xu R, Wang M, Liao Q, Huang K, Shan Z, et al. Association of HLA-DQB1*03:01 and DRB1*11:01 with spontaneous clearance of hepatitis C virus in Chinese Li ethnicity, an ethnic group genetically distinct from Chinese Han ethnicity and infected with unique HCV subtype. J Med virology. (2019) 91:1830–6. doi: 10.1002/jmv.25531
16. Khor SS, Ueno K, Nishida N, Kawashima M, Kawai Y, Aiba Y, et al. Novel HLA allele associations with susceptibility, staging, symptomatic state, autoimmune hepatitis and hepatocellular carcinoma events for primary biliary cholangitis in the Japanese population. Front Immunol. (2023) 14:1151502. doi: 10.3389/fimmu.2023.1151502
17. Hoofnagle JH, Bonkovsky HL, Phillips EJ, Li YJ, Ahmad J, Barnhart H, et al. HLA-B*35:01 and green tea-induced liver injury. Hepatol (Baltimore Md). (2021) 73:2484–93. doi: 10.1002/hep.31538
18. Fontana RJ, Cirulli ET, Gu J, Kleiner D, Ostrov D, Phillips E, et al. The role of HLA-A*33:01 in patients with cholestatic hepatitis attributed to terbinafine. J hepatology. (2018) 69:1317–25. doi: 10.1016/j.jhep.2018.08.004
19. Meng L, Pan Y, Liu Y, He R, Sun Y, Wang C, et al. Individuals carrying the HLA-B*15 allele exhibit favorable responses to COVID-19 vaccines but are more susceptible to Omicron BA.5.2 and XBB.1.16 infection. Front Immunol. (2024) 15:1440819. doi: 10.3389/fimmu.2024.1440819
20. Rahmouni M, De Marco L, Spadoni JL, Tison M, Medina-Santos R, Labib T, et al. The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control. Front Immunol. (.(2023) 14:1305856. doi: 10.3389/fimmu.2023.1305856
21. Wang F, Dong L, Wang W, Chen N, Zhang W, He J, et al. The polymorphism of HLA-A, -C, -B, -DRB3/4/5, -DRB1, -DQB1 loci in Zhejiang Han population, China using NGS technology. Int J immunogenetics. (2021) 48:485–9. doi: 10.1111/iji.12554
22. Tao S, You X, Chen N, Dong L, Zhao S, He Y, et al. HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 alleles in Zhejiang Han population. Immunogenetics. (2024) 76:305–14. doi: 10.1007/s00251-024-01349-5
23. Lu G, Lu Y, He Y, Chen W, and Zhu F. Impact of human leukocyte antigen mismatch between donor-recipient on acute rejection in liver transplantation using next-generation sequencing: a single-center study. Front Immunol. (2025) 16:1576815. doi: 10.3389/fimmu.2025.1576815
24. Excoffier L and Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol resources. (2010) 10:564–7. doi: 10.1111/j.1755-0998.2010.02847.x
25. Ou G, Liu X, Xu H, Ji X, Liu X, and Wang J. Variation and expression of HLA-DPB1 gene in HBV infection. Immunogenetics. (2021) 73:253–61. doi: 10.1007/s00251-021-01213-w
26. Nishida N, Sawai H, Kashiwase K, Minami M, Sugiyama M, Seto WK, et al. New susceptibility and resistance HLA-DP alleles to HBV-related diseases identified by a trans-ethnic association study in Asia. PloS One. (2014) 9:e86449. doi: 10.1371/journal.pone.0086449
27. Karra VK, Chowdhury SJ, Ruttala R, Gumma PK, Polipalli SK, Chakravarti A, et al. HLA-DQA1 & DQB1 variants associated with hepatitis B virus-related chronic hepatitis, cirrhosis & hepatocellular carcinoma. Indian J Med Res. (2018) 147:573–80. doi: 10.4103/ijmr.IJMR_1644_15
28. Huang J, Wang H, Fan ST, Zhao B, Zhang Z, Hao L, et al. The national program for deceased organ donation in China. Transplantation. (2013) 96:5–9. doi: 10.1136/gut.48.5.714
29. Cholankeril G and Ahmed A. Alcoholic liver disease replaces hepatitis C virus infection as the leading indication for liver transplantation in the United States. Clin Gastroenterol hepatology: Off Clin Pract J Am Gastroenterological Assoc. (2018) 16:1356–8. doi: 10.1016/j.cgh.2017.11.045
30. Chahal D, Marquez V, Hussaini T, Kim P, Chung SW, Segedi M, et al. End stage liver disease etiology & transplantation referral outcomes of major ethnic groups in British Columbia, Canada: A cohort study. Medicine. (2021) 100:e27436. doi: 10.1097/MD.0000000000027436
31. Liu C and Cheng B. Association of polymorphisms of human leucocyte antigen-DQA1 and DQB1 alleles with chronic hepatitis B virus infection, liver cirrhosis and hepatocellular carcinoma in Chinese. Int J immunogenetics. (2007) 34:373–8. doi: 10.1111/j.1744-313X.2007.00702.x
32. Ashouri S, Khor SS, Hitomi Y, Sawai H, Nishida N, Sugiyama M, et al. Genome-wide association study for chronic hepatitis B infection in the thai population. Front Genet. (2022) 13:887121. doi: 10.3389/fgene.2022.887121
33. Marsh SG and Bodmer JG. HLA-DR and -DQ epitopes and monoclonal antibody specificity. Immunol Today. (1989) 10:305–12. doi: 10.1016/0167-5699(89)90086-8
34. Wu TW, Chu CC, Ho TY, Chang Liao HW, Lin SK, Lin M, et al. Responses to booster hepatitis B vaccination are significantly correlated with genotypes of human leukocyte antigen (HLA)-DPB1 in neonatally vaccinated adolescents. Hum Genet. (2013) 132:1131–9. doi: 10.1007/s00439-013-1320-5
35. Huang YH, Liao SF, Khor SS, Lin YJ, Chen HY, Chang YH, et al. Large-scale genome-wide association study identifies HLA class II variants associated with chronic HBV infection: a study from Taiwan Biobank. Alimentary Pharmacol Ther. (2020) 52:682–91. doi: 10.1111/apt.15887
36. Wu X, Liu Y, Jin S, Wang M, Jiao Y, Yang B, et al. Single-cell sequencing of immune cells from anticitrullinated peptide antibody positive and negative rheumatoid arthritis. Nat Commun. (2021) 12:4977. doi: 10.1038/s41467-021-25246-7
37. Caillier SJ, Briggs F, Cree BA, Baranzini SE, Fernandez-Viña M, Ramsay PP, et al. Uncoupling the roles of HLA-DRB1 and HLA-DRB5 genes in multiple sclerosis. J Immunol (Baltimore Md: 1950). (2008) 181:5473–80. doi: 10.4049/jimmunol.181.8.5473
38. Munir S, Rahman SB, Rehman S, Saba N, Mazhar K, and Naluai ÅT. The association of HLA alleles and haplotypes with age of disease onset in Pakistani psoriatic patients. Int Arch Allergy Immunol. (2023) 184:202–10. doi: 10.1159/000527359
39. O’Leary JG, Kaneku H, Jennings LW, Bañuelos N, Susskind BM, Terasaki PI, et al. Preformed class II donor-specific antibodies are associated with an increased risk of early rejection after liver transplantation. Liver Transpl. (2013) 19:973–80. doi: 10.1002/lt.23687
40. Pidala J, Lee SJ, Ahn KW, Spellman S, Wang HL, Aljurf M, et al. Nonpermissive HLA-DPB1 mismatch increases mortality after myeloablative unrelated allogeneic hematopoietic cell transplantation. Blood. (2014) 124:2596–606. doi: 10.1182/blood-2014-05-576041
Keywords: end-stage liver disease, liver transplantation, human leukocyte antigen, next-generation sequencing, susceptibility
Citation: Lu G, He Y, Lu Y, Chen W, Lu C and Zhu F (2025) Association of HLA-DQA1*06:01 and DPB1*05:01:01G alleles with susceptibility to end-stage liver disease with HBV infection in liver transplant recipients from the Zhejiang Han population, China. Front. Immunol. 16:1684437. doi: 10.3389/fimmu.2025.1684437
Received: 12 August 2025; Accepted: 19 November 2025; Revised: 11 November 2025;
Published: 02 December 2025.
Edited by:
Mohammad Hossein Karimi, Shiraz University of Medical Sciences, IranCopyright © 2025 Lu, He, Lu, Chen, Lu and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Chen, Y3cwMDMxQDE2My5jb20=; Caide Lu, bHVjYWlkZUBuYnUuZWR1LmNu; Faming Zhu, emZtMDBAaG90bWFpbC5jb20=
Yangfang Lu3