Abstract

Background. Genetic studies of end-stage renal disease (ESRD), including those of human leukocyte antigen (HLA) genes, have been reported in several populations but have not yet been evaluated in Indonesia. Some studies have reported that these genes had a substantial role in ESRD. This study aims to analyze the association between HLA genes and ESRD within the Indonesian community. Method. A retrospective study to investigate HLA class I and II alleles to find out the distribution of HLA-A, -B, -C, -DPB1, -DQB1, and -DRB1 in renal transplant recipients and to ascertain their role in susceptibility to ESRD was performed on totally 149 subjects, consisting of 69 ESRD patients and 80 healthy controls. HLA typing was determined using Luminex techniques. The allele and haplotype frequencies were compared between ESRD patients and controls. Result. High-frequency alleles were HLA-A24 (43.6%), B15 (38.2%), C08 (30.8%), DRB112 (47.3%), DQB103 (50.6%), and DPB113 (22.5%). HLA-A24 and HLA-B35 were associated with a protective effect, with OR 0.537 (95%CI 0.34–0.86) and 0.316 (95%CI 0.11–0.88), respectively. There were some two-locus haplotypes associated with susceptibility to ESRD, such as B15-DRB112, B13-DRB115, A02-B15, A02-C08, and B13-DQB105. HLA-A02-B15-DRB112 and A24-B13-DRB115 appear to be associated with susceptibility to ESRD. Conclusion. The allele groups of HLA-A24 and HLA-B35 are associated with protection from ESRD. Meanwhile, HLA-B13-DRB115 and A24-B13-DRB115 are the most frequent HLAs associated with ESRD in two-locus and three-locus haplotype, respectively.

1. Introduction

The incidence of ESRD in the world is increasing as well as in Indonesia. The Indonesian Renal Registry showed that ESRD incidence rose from 9649 in 2010 to 30831 in 2017. At the same time, the prevalence of ESRD increased from 11484 in 2010 to 77892 in 2017. The causes of kidney disease are mainly due to hypertension and diabetes mellitus, with the 1-year survival rate of 83%, while the 5-year survival rate was 51.9%. The leading cause of death in 37% of cases was cardiovascular disorders [1].

End-stage renal disease (ESRD) is suspected to be influenced by genetic and nongenetic factors. Many genetic factors cause ESRD, and human leukocyte antigen (HLA) is an essential factor. HLA is encoded by the major histocompatibility complex (MHC) located on chromosome 6p21.3. The HLA molecule binds and presents peptide to T lymphocytes in cell-mediated immune response and plays a crucial role in shaping the T-cell repertoire and is also associated with allograft rejection [2, 3]. The correlation of HLA polymorphism with ESRD can be caused by HLA association with the aetiology and the progression of kidney disease. For example, in a Mexican population with type 2 diabetes mellitus (T2DM), HLA-DRB115:02 is associated with diabetic kidney disease. Meanwhile, specific HLA can also trigger a profibrogenic T-cell phenotype, possibly contributing to the progression of kidney disease severity [3]. Understanding the genetic factors that influence kidney disease is essential to help estimate the prognosis and prevention of ESRD.

Successful transplantations correlate with the HLA compatibility level between recipients and donors. Renal allotransplantation improves the quality of life and increases survival in patients with ESRD compared with long-term dialysis treatment. Primary renal allograft loss is associated with significant mortality. Therefore, the information about frequencies of HLA haplotype occurs in a population, and the association between HLA and ESRD in the Indonesian community is important. In this study, the distributions of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 in the renal transplant recipients and donors are obtained to determine HLA antigen frequencies.

2. Methods

2.1. Population Samples

This study was observational and retrospective. The data were collected from the biomolecular laboratory databases of Dr. Saiful Anwar General Hospital from 2017 to 2020. We investigated HLA class I and II alleles to determine the distribution of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 and to ascertain their role in susceptibility to ESRD in patients considering renal transplantation and healthy controls. The inclusion criteria were the subjects had complete HLA typing data from 6 loci, not family related, and Indonesian population.

2.2. DNA Extraction

DNA was isolated from whole blood containing Acid Citrate Dextrose (ACD) solution, using a Genomic DNA Isolation kit, according to the manufacturer’s instructions (Qiagen). In brief, 20 μL of Qiagen protease and 200 μL of whole blood were added into a 1.5 mL microcentrifuge tube. 200 μL of buffer (AL) was added and then incubated at 56°C for 10 minutes, and 200 μL of ethanol (100%) was added and centrifuged. The mixture was transferred into a QIAamp mini spin column, the tube was centrifuged, the QIAamp mini spin column was placed in a clean 2 mL collection tube, and the tube containing the filtrate was discarded. 500 μL of buffer (AW1) was added and centrifuged, then the QIAamp mini spin column was placed in a clean 2 mL collection tube, and the tube containing the filtrate was discarded. The procedure was repeated with 500 μL of buffer (AW2), and then, the QIAamp mini spin column was placeed in a clean 1.5 mL microcentrifuge tube. Then, 200 μL of elusion buffer (AE) was added and centrifuged, the QIAamp mini spin column was discarded, and the microcentrifuge tube containing the eluted DNA was stored at −80°C.

2.3. DNA Typing

HLA typing of loci HLA-A, HLA-B HLA -C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 was performed by Luminex Labscan 100 using the sequence-specific oligonucleotide probes (SSO) method. In summary, 2 μL of the target DNA is amplified using a group-specific primer. The PCR product is denatured and allowed to rehybridize to complementary DNA probes conjugated to fluorescently coded microspheres. A flow analyzer, the LABScan™ 100, identified the fluorescent intensity of PE (phycoerythrin) on each microsphere. The HLA Fusion software (ONE LAMBDA, INC, USA) will analyze the HLA typing for allele identification based on the reaction pattern compared to patterns associated with published HLA gene sequences.

2.4. Statistical Analysis

HLA allele, two-locus, and three-locus haplotype frequencies were performed by direct counting in percentage. Two-locus haplotype frequencies in this study will be compared with the two-locus haplotype from the other population because more studies reported the two-locus haplotype than the three-locus haplotype. However, the three-locus haplotype (HLA-ABDR) was considered more relevant in this regard, so we also performed the three-locus haplotype in this study. The association between HLA allele with ESRD was estimated with Odds Ratio (OR) and 95% confident intervals, with was considered statistically significant. Odds Ratio (OR) value greater than one was deemed positive or susceptible to ESRD, meanwhile OR values less than one were considered protective to ESRD. The exact test was used to evaluate the assumption of Hardy–Weinberg Equilibrium (HWE). The software for statistical analysis was the SPSS version 23 and OpenEpi. The institutional ethics board of Dr. Saiful Anwar General Hospital has approved the study (ethical number 400/141/K.3/302/2020).

3. Results

3.1. Characteristics of Study Participants and Hardy–Weinberg Equilibrium Tests

The study collected HLA data from 81 ESRD and 87 control participants. According to the inclusion criteria, 69 ESRD patients and 80 donors, a total of 149 subjects, were selected. There was a male predominance among the control group (63.8%) and ESRD group (51.3%). The allele number of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 in ESRD patients and donors are summarized in column charts shown in Figure 1. The HWE tests of the HLA-A, -B, -C, -DRB1, -DQB1, and -DRP1 loci showed .

3.2. Allele Frequencies at HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 Loci in Renal Transplant Recipients and Donors

The high-frequency alleles of HLA class I was HLA-A24 (35.50% in ESRD,50.60% in control group, and 43.60% in both groups), HLA-B15 (37.70% in ESRD, 38.20% in control group, and 38.20% in both groups), and HLA-C08 (31.20% in ESRD, 30% in control group, and 30.80% in both groups). Table 1 shows the class I HLA types with frequency >1% in the ESRD or control groups. If the frequency is 1% or less, the data were not performed. According to Table 1, HLA-A24 and HLA-B35 were associated with a protective effect, with OR 0.537 (95%CI 0.34–0.86) and 0.316 (95%CI 0.11–0.88), respectively.

Table 2 indicates the class II HLA types with frequency >1% in the ESRD or control groups. The high-frequency alleles of HLA class II were HLA-DRB112 (44.20% in ESRD, 50% in control group, and 47.30% in both groups), HLA-DQB103 (49.30% in ESRD, 51.90% in control group, and 50.6% in both groups), and HLA-DPB113 (22.50% in ESRD, 22.50% in control group, and 22.5% in both groups). HLA class II was not associated with ESRD.

Table 3 shows the association of two-locus haplotypes with ESRD with a significant odds ratio between the two groups ( value <0.05). There are some haplotypes associated with susceptibility to ESRD, namely, B15-DRB112, B13-DRB115, A02-B15, A02-C08, and B13- DQB105. Meanwhile, the other haplotypes are associated with protection from ESRD.

Table 4 indicates the association of the three-locus haplotypes with ESRD with a significant odds ratio between the two groups ( value <0.05). Two haplotypes are associated with susceptibility to ESRD, where the highest OR is HLA-A24-B13-DRB115 (OR = 8.171), followed by HLA-A02-B15-DRB112 (OR = 5.724). The haplotype of HLA-A24-B35-DRB112 and HLA-A24-B35-DRB115 is negatively associated with ESRD.

4. Discussion

Aside from its essential role in determining donor-recipient immune compatibility in organ transplantation, HLA genotyping is performed routinely as part of the diagnostic work-up of certain autoimmune diseases. The HLA genotyping also has contributed to understanding several diseases’ pathogenesis, including ESRD [3].

Our study observed that the ESRD participants consisted of 51.3% males and 48.7% females. It is similar to the research by Panigrahi et al. [4] and Tuladhar et al. [5], which reported that most ESRD patients were males.

The data of HLA alleles in our study showed heterogeneity in both HLA class I or class II. The number of HLA-B alleles is higher than HLA-A and HLA-C in HLA class 1. Meanwhile, HLA class II indicated that the number of HLA-DPB1 alleles is higher than HLA-DR and HLA-DQB1. The heterogeneity of HLA alleles results in amino acid substitutions that predominantly involve peptide binding sites for the effective display of a broad peptide to CD4+ and CD8+ T cells. As of July 2018, the polymorphisms of HLA-A, HLA-B, and HLA-C had 9,341 different proteins and 5,355 proteins in the polymorphisms of HLA-DR, HLA-DQ, and HLA-DP [3].

The high frequency of HLA-A24 in our study is also observed similarly in South Indian, northern Indian, and Bangladeshi Bengali populations [57]. The frequent alleles such as HLA-A02 and A11 were also observed by Ali et al., and this HLA correlated with the Bangladeshi Bengali population. The most common alleles in Nepal, such as HLA-B15 and DRB112 [5], were also found in high frequencies in our study. Pradana et al. stated that the Indonesian HLA profile had high variation at HLA-DQA106:01, DQB103:01, DRB112:02, A24:07, B15:02, B75, and B18:02 [8]. These data are consistent with our research that found the high frequency of HLA-A24, B15, DRB112, and DQB103. Polymorphism in the HLA system is also used as a tool for anthropological studies. Genetic distances and correspondence analysis demonstrated that the allele and haplotype distribution of class I and class II loci are racially and geographically restricted [9, 10].

HLA-A24 had a negative correlation with ESRD, according to our data. However, Cao et al. reported that HLA-A24 in ESRD patients is significantly higher than in the controls and associated with ESRD in the Cantonese population [11]. The different result of the study with our study is caused by that the HLA expression that affects ESRD may be influenced by not only one locus but also many loci of HLA. For example, according to our study, HLA-A24-B13-DRB115 is strongly associated with ESRD susceptibility. So, the role of HLA-A24 in chronic kidney disease and ESRD still needs more study in the future because the HLA system plays a pivotal role in the antigen presentation of intracellular and extracellular peptides and the regulation of innate and adaptive immune responses. Considering the ability of HLA to influence the thymic selection and peripheral energy of T cells, its role in the pathogenesis of kidney disease needs to be investigated [12].

HLA-B35 had a negative association with ESRD based on our study. However, in Brazil, the research by Ravazzi-Gauch et al. stated that the most common alleles for each locus among ESRD or renal transplantation patients were HLA-A02, B35, and DRB11 [13]. In China, the HLA-B locus showed HLA-B15:01, B55:02, and B39:05 emerged as susceptible alleles, whereas no protective allele was found [14]. Meanwhile, in our study, the primary HLA-B allele in the kidney transplant patient group was HLA-B15. In the HLA system, allele and haplotype distribution differ between various ethnic groups or between the members of the same ethnic group living in different geographic areas, as shown in previous studies [14, 15].

HLA-C, HLA-DR, HLA-DQ-B1, and HLA-DPB1 did not associate with ESRD, according to our data. This result is different from the study by Lowe et al. that reported HLA association with kidney function on 401,307 white British subjects. They showed that 22 types of HLA were associated with increased eGFR and 11 types of HLA were associated with decreased eGFR. As many as 7 HLA out of 11 types are associated with CKD or ESRD, namely, HLA-A01:01, B08:01, C07:01, DRB103:01, DRB301:01, DQA105:01, and DQB102:01 [16]. In the ESRD group, our study indicated that the most HLA-DPB1 locus was HLA-DPB113, but there was no significant difference with the control group. The association between HLA-DPB1 and ESRD has never been reported in Indonesia, but some studies in Indonesia stated that HLA-DPB1 is associated with hepatitis B infection susceptibility [1719].

Some two-locus haplotypes were associated with susceptibility to ESRD in our study, such as HLA-A02-B15, A02-C08, B15-DRB112, B13-DRB115, and B13-DQB105. This result is different from the research by Hamdi that reported the association of two-locus haplotypes HLA-A01-DRB113 and HLA-A30-DRBI03 with ESRD in Saudi Arabian population [20]. In Venezuela, a study by Rivera et al. stated that the haplotypes positively associated with ESRD were HLA-A02-B51, HLA-A02-B53, HLA-A23-B38, and HLA-A68-B38 [21].

According to our data, HLA-A02-B15-DRB112 and HLA-A24-B13-DRB115 were three-locus haplotypes associated with susceptibility to ESRD. HLA-A24-B13-DRB115 had a strong association with ESRD with OR 8.171. Nevertheless, the study by Cao et al. in the Cantonese community, a representative southern population of China, showed that HLA-A11-B27-DRB104 in ESRD patients is significantly higher than that in the controls [11].

Based on the data of HLA and ESRD given above, there are different results in different studies and populations. This difference may be influenced by the disease that causes kidney abnormalities. HLA has substantial risk factors in most immune-mediated renal disorders. Together with other genetic and environmental factors, HLA causes loss of tolerance and autoimmune-mediated inflammation of the kidney [3, 22]. For example, Karnes et al. reported that HLA-DRB104 and HLA-DQB103:02 had an association with diabetic kidney disease. The association of HLA-DQB103:02 and ESRD or kidney transplantation was weak (OR = 1.4), but this HLA had the risk of type 1 diabetes and diabetic kidney disease with OR 7.1 [23].

The limitations of our study were the small number of ESRD and healthy control participants. Further research with a more significant number of samples is needed. However, our research is expected to give a piece of information on the association between HLA and ESRD, especially in Indonesia, thus allowing the management of chronic kidney disease patients with HLA susceptibility more precisely and effectively.

5. Conclusions

The result of this study showed heterogeneity in both HLA class I or class II antigens, and some HLA polymorphisms have an association with ESRD in the Indonesian population. HLA-A24 and HLA-B35 are associated with protection from ESRD, whereas HLA-B13-DRB115 is the most frequent HLA related to susceptibility to ESRD in the two-locus haplotype and HLA-A24-B13-DRB115 in the three-locus haplotype.

Data Availability

The database for the study can be acquired from the principal investigator, Hani Susianti, [email protected].

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors thank Dr. Saiful Anwar General Hospital’s director and all participants involved in this study.