ORIGINAL RESEARCH article

Front. Immunol., 05 January 2026

Sec. Immunological Tolerance and Regulation

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1697839

Genetic variants of the HLA-G/LILRB1 ligand-receptor axis in donors or recipients are prognostic covariates for rejection after living kidney transplantation

  • 1. Institute for Transfusion Medicine, University Hospital Essen, University Duisburg-Essen, Essen, Germany

  • 2. Department of Nephrology, University Hospital Essen, University Duisburg-Essen, Essen, Germany

  • 3. HLA-Laboratory, Stefan-Morsch-Foundation, Birkenfeld, Germany

  • 4. Department of Infectious Diseases, University Hospital Essen, University Duisburg-Essen, Essen, Germany

Abstract

Background:

HLA-G is a non-classical HLA class I molecule that promotes transplant tolerance. It engages the inhibitory receptor LILRB1 on immune effector cells, suppressing cytotoxic responses and inflammation, while promoting tolerogenic and regulatory immune phenotypes. Polymorphisms in the HLA-G 3′ untranslated region (3′UTR) modulate HLA-G expression levels, and LILRB1 promoter variants influence receptor expression. The combined effect on kidney transplant (KTx) rejection has not been systematically studied.

Methods:

Living donor–recipient pairs undergoing KTx were genotyped for nine variants in the HLA-G 3′UTR region and two single nucleotide polymorphisms (SNPs) in the LILRB1 promoter (PROMO) regions. Haplotypes were arranged for both loci. Clinical endpoints were biopsy-proven T cell-mediated rejection (TCMR) within one year and antibody-mediated rejection (AMR) within five years post-transplant.

Results:

Donor positivity for HLA-G 3′UTR-1 or UTR-2 or negative for UTR-3 haplotype were associated with a significantly higher risk of TCMR in both univariate or multivariate analyses. Recipients lacking the LILRB1-PROMO CG haplotype also had an increased TCMR risk. The combination of an HLA-G 3’UTR-2 positive donor with a LILRB1-PROMO CG haplotype negative recipient was found to be an independent predictor of TCMR. In contrast, HLA-G 3′UTR variants were not associated with AMR, while the presence of the recipient LILRB1-PROMO CG haplotype emerged as an independent AMR risk factor.

Conclusions:

Donor HLA-G 3’UTR and recipient LILRB1-PROMO haplotypes define a functional immunogenetic axis that differentially influence TCMR and AMR. These results support the clinical potential of HLA-G/LILRB1 genetic profiling to improve donor selection in living KTx and to guide the development of novel rejection therapies.

1 Introduction

Human leukocyte antigen G (HLA-G) is a non-classical HLA class I molecule with potent immunomodulatory functions (). Unlike classical HLA class I antigens, HLA-G exhibits limited polymorphism, restricted tissue distribution, and multiple isoforms, including both membrane-bound and soluble variants (). Under physiological conditions, HLA-G expression is predominantly restricted to immune-privileged sites, such as the maternal–fetal interface, where it contributes to immune tolerance (, ). In pathological contexts, including transplantation, infection, autoimmunity, and malignancy, HLA-G expression can be induced and contributes to immune evasion by suppressing the activation and effector functions of various immune cells ().

The immunosuppressive activity of HLA-G is primarily mediated through its interaction with inhibitory receptors expressed on immune effector cells (). One of the key receptors is immunoglobulin-like transcript 2 (ILT2), also referred to as LILRB1 or CD85j (). LILRB1 is expressed on a broad spectrum of immune cells, including subsets of CD4+ and CD8+ T cells, natural killer (NK) cells, B cells, and myeloid antigen-presenting cells (APCs). Notably, LILRB1 expression also varies by cell type, with high levels on monocytes and B cells and intermediate to low levels on NK cells and T cells. Upon engagement with HLA-G, LILRB1 transmits inhibitory signals that suppress cytotoxic activity, reduce cytokine production, inhibit cell proliferation, and promote tolerogenic or regulatory immune phenotypes (, ). Although HLA-G is its primary ligand, LILRB1 also binds classical HLA class I molecules and HLA-F, suggesting a broader role in fine-tuning immune responses (28, 29).

The HLA-G-LILRB1 interaction serves as a critical immune checkpoint that is essential for peripheral immune tolerance (30, 31). The efficiency of this interaction depends on the expression levels of both HLA-G ligand and LILRB1 receptor. Polymorphisms in the 3′ untranslated region (3′UTR) of the HLA-G, including multiple single nucleotide polymorphisms (SNPs) and the well characterized 14 base pair insertion/deletion (14 bp Ins/Del), modulate mRNA stability, splicing, and microRNA targeting, collectively shaping HLA-G expression (3236). Functional studies have shown that the 14 bp insertion is associated with reduced mRNA stability and lower soluble HLA-G (sHLA-G) levels, whereas the +3142G allele enhances miR-148a/miR-152 binding and downregulates expression. These variants form defined haplotypes with differential expression and consequently immunoregulatory potential (37). Consequently, UTR-1 and UTR-2 haplotypes are linked to higher sHLA-G expression, while UTR-3 is associated with lower levels. Similarly, LILRB1 expression is regulated by SNPs such as rs10416697C/G in the distal promoter and rs1004443G/A in the proximal promoter regions (38, 39). These SNPs may alter transcriptional activity and surface receptor density. The rs10416697G and rs1004443A alleles are associated with higher promoter activity and increased LILRB1 surface expression in certain immune effector cells.

In kidney transplantation (KTx), alloimmune recognition of donor antigens by recipient immune system is a principal cause of rejection (40). Rejection is classified as T cell-mediated rejection (TCMR), involving direct cytotoxic and inflammatory responses, or antibody-mediated rejection (AMR), characterized by donor-specific antibodies (DSAs), complement activation, and endothelial injury. TCMR involves infiltration of CD8+ cytotoxic T cells and CD4+ helper T cells, leading to tissue destruction and inflammation. In contrast, AMR results from B cells activation and DSA production. Those antibodies bind to donor HLA molecules, activating the classical complement cascade, and inducing microvascular inflammation. NK cells and monocytes contribute to antibody-dependent cellular cytotoxicity through Fc receptor engagement, amplifying endothelial injury. C4d deposition in peritubular capillaries is a diagnostic hallmark. While TCMR can often be managed by immunosuppressive therapy, AMR remains difficult to treat and is closely associated with chronic graft dysfunction (41, 42).

Given its central role in immune regulation, the HLA-G–LILRB1 axis is of significant interest in transplantation immunology. Increased HLA-G expression has been associated with reduced rejection and improved graft survival. This effect being likely mediated through suppression of both cellular and humoral immune responses. However, the clinical relevance of this pathway may depend on the genetic variability of both ligand and receptor. While previous studies have evaluated the individual contributions of HLA-G 3′UTR polymorphisms, the role of LILRB1 promoter variants and the combined effect of donor HLA-G 3′UTR haplotypes and recipient LILRB1 haplotypes have not been systematically investigated.

This study aims to focus on the joint impact of donor and recipient HLA-G 3′UTR haplotypes and LILRB1 promoter (PROMO) SNPs (rs10416697 and rs1004443) on rejection risk following KTx. Understanding this genetic interplay may provide novel insights into the mechanisms of transplant tolerance and contribute to the development of genotype-informed strategies for risk assessment, therapeutic targeting, and long-term graft survival optimization.

2 Materials and methods

2.1 Study population and clinical data

This study included 293 living donor–recipient pairs who underwent living KTx at the University Hospital Essen, Germany, between 2005 and 2017 as part of the institutional living donor kidney transplant program. Clinical and demographic data for donors and recipients were obtained from institutional medical records and transplant databases (Table 1). Donors were predominantly under 60 years of age (82.6%) and 56% were female. Recipients were mostly male (61.1%) and under 60 years of age (89.4%). The median cold and warm ischemia times were 135 minutes (range: 45–407) and 19 minutes (range: 10–76), respectively. High immunological disparity was present, with >2 HLA-A, -B, or -DR mismatches in 87.4% of cases, sex mismatches in 66.9%, and ABO-incompatible transplantation in 14.7%. Pre-transplant panel-reactive antibodies were detected in 11.9% of recipients, while 13.0% developed post-transplant donor-specific anti-HLA antibodies (DSAs). All patients received standard post-transplant follow-up care. The primary clinical endpoints were biopsy-proven TCMR within the first 12 months and AMR, diagnosed within five years post-transplant. Rejection episodes were classified according to the Banff classification system applicable at the time of biopsy (40, 4346), based on histopathological criteria and the presence of DSAs where relevant. All biopsies were performed based on clinical indication. Induction immunosuppression consisted of basiliximab (anti-IL-2 receptor CD25 monoclonal antibody), followed by maintenance therapy with a calcineurin inhibitor (tacrolimus or cyclosporine), corticosteroids, and an antiproliferative agent (mycophenolate mofetil or azathioprine), in line with institutional and national guidelines. Antithymocyte globulin (ATG) was administered in 4.8% of recipients. No significant association was observed between the etiology of primary kidney disease and the key clinical outcomes assessed in this study.

Table 1

Donor (N=293)Number
 Sex (men/women)128/165
 Age (<60 years/>60 years)242/51
Recipient (N=293)
 Sex (men/women)179/114
 Age (<60 years/>60 years)262/31
KTx related parameters
 Cold ischemia time, Median (range); minutes135 (45 – 407)
 Warm ischemia time, Median (range); minutes19 (10 – 76)
Mismatch situation (MM)
 HLA-A,-B, or -DR MM (<2MM/>2MM)37/256
 Sex MM (no/yes)97/196
 AB0iTx (no/yes)250/43
Panel reactive antibody (PRA)
 PRA pre-KTX (yes/no)35/258
Initial Immunosuppressive therapy
 ATG/Basiliximab based induction therapy14/279
 CNI-administration yes/no288/5
 MMF co-administration yes/no253/40
 Steroid co-administration yes/no290/3
 other1/292

Demographic and clinical characteristics of recipients and donors at baseline.

HLA, Human Leukocyte Antigen; ATG, antithymocyte globulin; CNI, calcineurin inhibitor; MMF, mycophenolate mofetil; KTx, kidney transplant, AB0iTx, AB0 incompatible transplantation.

All participants provided written informed consent. The study protocol was approved by the local ethics committee of the Medical Faculty of the University of Duisburg-Essen (approval number 12-5312-BO) and conducted in accordance with the Declaration of Helsinki.

2.2 HLA-G 3′UTR and LILRB1 rs10416697/rs1004443 genotyping

Genomic DNA was extracted from peripheral blood leukocytes using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) in accordance with the manufacturer’s protocol. The 3′ untranslated region (3′UTR) of the HLA-G gene was amplified by polymerase chain reaction (PCR), and sequencing was performed as recently described (, ). Sequence chromatograms were evaluated using FinchTV software version 1.4.0 (Geospiza Inc.). The sequencing analysis covered 10 polymorphic positions in the HLA-G 3′UTR, including the 14 bp insertion/deletion (+2961) and SNPs at +3001 C/T, +3003 C/T, +3010 C/G, +3027 A/C, +3035 C/T, +3142 C/G, +3187 A/G, +3196 C/G, and +3227 A/G.

The LILRB1 promoter (LILRB1-PROMO) variants of rs10416697C/G (47) and of rs1004443G/A were defined using TaqMan SNP Genotyping Assays (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Genotyping was performed on a real-time PCR platform (Quant Studio 6 Real-Time-PCR-System, Applied Biosystems, Thermo Fisher Scientific), and allele calls were made by endpoint fluorescence analysis using standard software provided by Thermo Fisher.

Haplotype inference was conducted using PHASE version 2.1 with default parameters (). Inferred haplotypes had posterior probabilities ranging from 0.98 to 1.0 and demonstrated consistency across 10 independent runs.

2.3 Statistics

Statistical analyses were performed using SPSS Statistics version 27 (IBM Corp., Armonk, NY, USA), GraphPad Prism version 10.0 (GraphPad Software, San Diego, CA, USA), R (Rstudio), or MedCalc® Statistical Software version 23.2.1 (MedCalc Software Ltd, Ostend, Belgium; https://www.medcalc.org; 2025). Graphical abstract was performed using BioRender. Allelic and genotypic distributions were assessed for deviation from Hardy–Weinberg equilibrium using Haploview version 4.2 (Broad Institute, Cambridge, MA, USA), which was also used to construct linkage disequilibrium (LD) plots and haplotype blocks based on LILRB1 promoter (LILRB1-PROMO) SNPs. Associations between HLA-G 3′UTR haplotypes, LILRB1-PROMO haplotypes, and clinical endpoints (including biopsy-proven rejection) were evaluated using univariate analysis Fisher’s exact test. Kaplan-Meier survival analysis with log-rank tests were performed by survminer R package version 0.4.9 (https://CRAN.R-project.org/package=survminer). Multivariate analysis was conducted using Cox proportional hazards regression to identify independent predictors of rejection, adjusting for relevant clinical and transplant-related variables. A two-tailed p-value < 0.05 was considered statistically significant.

3 Results

3.1 Distribution of HLA-G 3′UTR haplotypes in living kidney transplant recipients and donors

HLA-G 3′UTR haplotypes were determined in 280 kidney transplant recipients and 279 living donors. Nine haplotypes with a frequency >1% were identified in both groups. HLA-G 3′UTR-1 and HLA-G 3′UTR-2 were the most prevalent, accounting for 36.8% and 31.9% of haplotypes in recipients and 32.8% and 30.1% in donors, respectively. Haplotype distributions did not differ significantly between groups, consistent with prior population data (see Table 2, ()). However, homozygosity for the HLA-G 3′UTR-1 haplotype was significantly more frequent in kidney transplant recipients than in donors (p = 0.002, Supplementary Table 1). In recipients, HLA-G 3′UTR-1 homozygosity was associated with underlying primary diseases such as glomerulonephritis (see Supplementary Data).

Table 2

HaplotypeFrequency recipientFrequency donorpaOR (95% CI)Frequency europe*
UTR-136.832.80.1611.19 (0.93 - 1.53)34.7
UTR-231.930.10.5021.09 (0.84 - 1.40)28.3
UTR-38.48.80.8170.95 (0.63 - 1.44)9.1
UTR-414.815.10.9130.98 (0.71 - 1.36)14.8
UTR-52.33.60.2120.63 (0.30 - 1.28)3.3
UTR-61.9160.6571.22 (0.53 - 3.07)0.5
UTR-76.15.00.4421.22 (0.73 - 2.02)5.5
UTR-131.41.80.6290.79 (0.32 – 2.04)0.2
UTR-182.12.0>0.9991.06 (0.49 - 2.17)3.3

HLA-G 3’UTR haplotype frequencies of recipients (n=280) and donors (n=279).

Haplotype phasing was assessed by PHASE 2.1 software using default parameters. Only haplotypes with frequencies >1 % were listed. UTR – untranslated region. pa-values were calculated by GraphPad Prism using two-sided Fisher’s exact test when evaluating haplotypes frequencies genotypes; alpha<0.05; OR-odds ratio; CI -confidence interval; *:Literature: https://doi.org/10.1038/s41598-021-02106-4.

3.2 Distribution of LILRB1 promoter SNP genotypes in living kidney transplant recipients and donors

Genotyping of two single nucleotide polymorphisms (SNPs) - rs10416697 and rs1004443 - located in the promoter regions of the LILRB1 gene was performed in kidney transplant recipients and living donors (Table 3). For rs10416697, there was no significant difference in C carrier frequency (CC + CG: 54.5% in recipients vs. 53.4% in donors; p = 0.8027; OR = 1.05, 95% CI: 0.76–1.44) or G carrier frequency (CG + GG: 91.3% vs. 88.4%; p = 0.2733; OR = 1.38, 95% CI: 0.81–2.33) was observed. For rs1004443 neither G carrier status (GG + AG: 53.2% vs. 53.8%; p = 0.9333; OR = 0.98, 95% CI: 0.71–1.35) nor A-carrier status (AA + AG: 91.5% vs. 87.9%; p = 0.1704; OR = 1.48, 95% CI: 0.87–2.55) differed significantly between recipients and donors. Overall, the distribution of LILRB1-PROMO SNP allelic frequencies was similar in both cohorts as well as in the European population-base data set ( (38), see Table 4).

Table 3

SNPGenotypesRecipient%Donors%pcOR (95% CI)
rs10416697CC258.73411.6
CG13145.812241.80.4135
GG13045.513646.6
C pos
CC +CG vs.15654.515653.40.80271.05 (0.76 - 1.44)
GG13045.513646.6
G pos
CG + GG vs.26191.325888.40.27331.38 (0.81 - 2.33)
CC258.73411.6
rs1004443AA13346.813446.2
AG12744.712141.70.3435
GG248.53512.1
G pos
GG + AG vs.15153.215653.80.93330.98 (0.71 - 1.35)
AA13346.813446.2
A pos
AA + AG vs.26091.525587.90.17041.48 (0.87 - 2.55)
GG248.53512.1

Distribution of LILRB1 SNP rs10416697 and rs1004443 genotypes detected at the LILRB1 promoter region in recipient (n=286/n=284) and donors (n=292/n=290).

Cpos, Positive for C allele. Gpos, Positive for G allele.Apos, Positive for A allele.

Table 4

SNPAlleleRecipient
%
Donors %pcOR (95% CI)Frequency Europe
rs10416697C31.632.50.75310.96 (0.75 - 1.23)30.0
G68.467.570.0
rs1004443A69.267.10.4481.10 (0.86 - 1.42)70.0
G30.832.930.0

Distribution of LILRB1 SNP rs10416697 and rs1004443 alleles detected at the LILRB1 promoter region in recipient (n=286/n=284) and donors (n=292/n=290).

3.3 LILRB1 promoter haplotype distribution in living kidney transplant recipients and donors

Both LILRB1-PROMO SNPs rs10416697 and rs1004443 were in strong linkage disequilibrium. Haplotype inference revealed four haplotypes (Table 5), which are referred to here as LILRB1-PROMO haplotype(s). The GA haplotype was most frequent in both groups (68.2% in recipients vs. 66.8% in donors; p = 0.611; OR = 1.07, 95% CI: 0.83–1.36), followed by CG (30.8% vs. 32.4%; p = 0.560; OR = 0.93, 95% CI: 0.72–1.19). The CA and GG haplotypes were rare (<1% in both groups) and did not differ significantly (CA: p = 0.988; GG: p = 0.186).

Table 5

rs10416697, rs1004443
Haplotype
Frequency recipientFrequency donorpaOR (95% CI)
CA0.20.20.9881.01 (0.05 - 19.43)
CG30.832.40.5600.93 (0.72 - 1.19)
GA68.266.80.6111.07 (0.83 - 1.36)
GG0.20.70.1860.25 (0.02 - 1.53)

LILRB1 gene promoter rs10416697 and rs1004443 haplotype (LILRB1-PROMO) frequencies of recipients and donors.

3.4 Donor HLA-G 3′UTR-1, -2, -3 haplotypes are associated with the incidence of TCMR in the first year post living-kidney transplantation

To evaluate the impact of donor HLA-G 3′UTR polymorphisms on early TCMR, the association between individual donor haplotypes and the incidence of biopsy-proven TCMR within the first 12 months following kidney transplantation was analyzed.

Recipients of grafts positive for HLA-G 3′UTR-1 or 3′UTR-2 showed a higher likelihood of TCMR in univariate analysis (Figures 1A–C) (p = 0.047; OR = 1.80 and (p = 0.038; OR = 1.84, respectively). In contrast, donor HLA-G 3′UTR-3 carrier status was significantly associated with a lower incidence of TCMR compared to negative status (p = 0.042; OR = 0.38). No significant associations were observed for other HLA-G 3′UTR haplotypes from donors or recipients (data not shown).

Figure 1

Taking into account the period of one year after KTx, these results were further corroborated by a Kaplan-Meier survival analysis with log-rank test (Figures 2A–C). A trend toward reduced TCMR-free survival probability was observed in HLA-G 3′UTR-1 haplotype-positive donor grafts, but statistical significance was not achieved (log-rank p = 0.052). Recipients who received a HLA-G 3′UTR-2-positive donor graft were associated with reduced TCMR-free survival (log-rank p = 0.031). In contrast, HLA-G 3′UTR-3-positive donor status was associated with prolonged TCMR-free survival probability (log-rank p = 0.047), suggesting a protective effect. Thus, donor HLA-G 3′UTR-1 and UTR-2 carrier status were related to greater TCMR risk, while HLA-G 3′UTR-3 was associated to protection in the univariate analysis.

Figure 2

3.5 Donor HLA-G 3′UTR-1 and HLA-G 3′UTR-2 haplotype status are independent predictors for TCMR post living-kidney transplantation in multivariate analysis

To assess the independent predictive value of donor HLA-G 3′UTR polymorphisms for TCMR, a multivariate Cox proportional hazards regression analysis was performed. The model was adjusted for transplant relevant clinical and immunological covariates, including ATG-based induction therapy, recipient and donor age, sex mismatch (MM), number of HLA-A, B, DR mismatch status (MM>2), ABO incompatibility (AB0i), panel-reactive antibody (PRA) status, cold and warm ischemia time. The forest plot illustrated the hazard ratios (HR) and 95% confidence intervals (CI) for each covariate (Figure 3). Among all variables tested, only the presence of donor HLA-G 3′UTR-1 (HR = 1.9, 95% CI: 1.1-3.4, p = 0.015) and HLA-G 3′UTR-2 (HR = 2.1, 95% CI: 1.2-3.6, p = 0.008) haplotypes were significantly associated with an increased risk of TCMR. In the Cox regression model, HLA-G 3′UTR-3 (HR: 0.47, 95% CI: 0.2-1.2, p = 0.130) did not reach the level of significance as an independent factor influencing TCMR risk. Furthermore, none of the clinical covariates included in the model were identified as independent prognostic covariates for TCMR. These results highlighted the independent and adverse impact of HLA-G 3′UTR haplotypes on TCMR.

Figure 3

3.6 Recipient LILRB1 promoter CG haplotype status is a covariate for TCMR post living-kidney transplantation in univariate analysis

To investigate the contribution of LILRB1-PROMO polymorphisms in modulating early alloimmune responses, we examined the association between the recipient LILRB1-PROMO haplotypes and the incidence of biopsy-proven TCMR within the first 12 months following kidney transplantation. Carriers of the LILRB1-PROMO CG haplotype exhibited a markedly lower incidence of TCMR (p = 0.01; OR: 0.48, Figure 4A). Consistently with this, Kaplan–Meier analysis with log-rank testing showed significant improved probability of 1-year TCMR-free survival in recipients carrying a LILRB1-PROMO CG haplotype compared to recipients lacking the LILRB1-PROMO CG haplotype (log-rank p = 0.011, Figure 4B). Importantly, no significant association was observed any donor LILRB1-PROMO haplotype and TCMR incidence (data not shown). Thus, these findings underscore the functional importance of LILRB1 genetic variation in the promoter regions of the recipient in shaping the cellular alloimmune response following transplantation. The presence of LILRB1-PROMO CG haplotype confers a protective effect against TCMR.

Figure 4

3.7 The donor HLA-G 3′UTR-1/2 status and recipient LILRB1-PROMO haplotype status are independent covariates for TCMR post living-kidney transplantation in multivariate analysis

To assess the individual contributions of donor HLA-G and recipient LILRB1-PROMO polymorphisms to the risk of TCMR), a multivariate Cox proportional hazards regression analysis was performed using clinical variables, consistent with those used in the evaluation of HLA-G 3′UTR haplotypes (Figure 5). The analysis demonstrated that donors positive for HLA-G 3′UTR-1 (HR = 1.9, 95% CI: 1.1-3.3, p = 0.023) or HLA-G 3′UTR-2 (HR = 2.0, 95% CI: 1.2-3.6, p = 0.011) haplotype had an significantly increased risk of TCMR within the first 12 months following transplantation. Additionally, the absence of LILRB1-PROMO CG haplotype in the recipient was independently associated with an increased risk of TCMR (HR = 1.48, 95% CI: 1.1-3.1, p = 0.021). No other genetic or clinical covariates reached statistical significance in the multivariate model.

Figure 5

3.8 The combination of an HLA-G 3′UTR-2-positive donor transplant and a LILRB1-PROMO CG haplotype-negative recipient is the strongest covariate for TCMR post living-kidney transplantation in the multivariate analysis

Since both donor and recipient haplotypes independently influence the risk of TCMR, and given the immunological significance of the HLA-G/LILRB1 checkpoint axis, an additional multivariate analysis was performed to assess whether specific donor-recipient haplotype combinations influence the occurrence of TCMR within the first year after transplantation. For this purpose, a high-risk combination was defined consisting of recipients lacking the LILRB1-PROMO-CG haplotype and donors carrying either the HLA-G 3′UTR-1 or HLA-G 3′UTR-2 haplotype (Figure 6A). The multivariate Cox regression analysis included the same transplant-related covariates as before. In this model, the only independent predictive TCMR risk factor was the immunogenetic pairing of recipients without the LILRB1-PROMO CG haplotype with donor transplants that were positive for HLA-G 3′UTR-2 (HR = 3.6, 95% CI: 1.6–8.5, p = 0.003). The pairing of LILRB1-PROMO CG negative-haplotype recipients with grafts positive for HLA-G 3′UTR-1 did not reach statistical significance in the multivariate analysis (HR = 1.7, 95% CI: 0.9–3.7, p = 0.12). None of the clinical covariates reached statistical significance. Accordingly, Kaplan-Meier survival analysis showed a significantly reduced probability of TCMR -free 1-year survival in recipients with a negative LILRB1-PROMO CG haplotype who received HLA-G 3′UTR-2-positive transplants compared to the other genetic combinations of donors and recipients (Log-rank p = 0.00079, Figure 6B).

Figure 6

3.9 Recipient LILRB1-PROMO CG positive-haplotype as a covariate for AMR post living-kidney transplantation in uni- and multivariate analysis

To evaluate the impact of genetic and clinical parameters on the risk of AMR a comprehensive analysis was performed. Neither the donor’s nor the recipient’s HLA-G 3′UTR haplotypes showed a significant association with the occurrence of AMR during a five-year observation period. In contrast to TCMR, where the presence of the LILRB1-PROMO CG haplotype in recipients has a protective effect against rejection post living-kidney transplantation (Figures 4, A, B), this haplotype was significantly more common in patients with AMR than in patients without this LILRB1-PROMO haplotype (78% vs. 50%, p = 0.007, Figure 7A). Kaplan–Meier survival analysis further supported this finding, showing significantly reduced AMR -free survival probability in recipients carrying LILRB1-PROMO CG haplotype compared to non-carriers (log-rank p = 0.0058, Figure 7B). Importantly, the multivariate Cox proportional hazards model including the clinical covariates as before, demonstrated the LILRB1-PROMO CG haplotype carrier status (HR = 2.6, 95% CI: 1.1-6.5, p = 0.003) in addition to pre-transplant panel-reactive antibody positivity as independent predictors for AMR (HR = 3.5, 95% CI: 1.4-8.8, p = 0.037, Figure 7C). The latter finding has been previously described as an independent immunological risk factor in other kidney transplant cohorts. No other clinical variables demonstrated significant associations.

Figure 7

4 Discussion

Living donor kidney transplantation remains the best available therapeutic option for patients with end-stage renal disease, offering superior long-term graft survival, reduced waiting times, and lower rates of delayed graft function compared to deceased donor transplantation (48, 49). Immune-mediated rejection, however, remains a major cause of graft loss. Identifying genetic factors that influence rejection risk is critical for improving transplant outcomes. In our study donor HLA-G 3′UTR haplotypes and recipient LILRB1 (PROMO) promoter polymorphisms (rs10416697 and rs1004443) were identified as immunogenetic relevant determinants of rejection risk following living kidney transplantation. Distinct donor–recipient haplotype constellations within the HLA-G/LILRB1 immune checkpoint axis differentially influence the incidence TCMR, whereas exclusively the recipient LILRB1-PROMO CG haplotype influences AMR. This suggests fundamentally different immune mechanisms underlying these two types of allograft injury.

HLA-G and its cognate LILRB1 receptor play a pivotal role in modulating alloimmune responses after transplantation (, 50). This pathway induces induction of tolerogenic dendritic cells and regulatory T cells, impairment the effector cell function of T-, NK- and B- cells, thus contributing to both local and systemic immune modulation (51). The expression of both receptor and ligand is tightly regulated by genetic polymorphisms (39, 52).

HLA-G, unlike classical HLA class I molecules, is selectively expressed in immune-tolerant environments such as the placenta, tumors, and transplanted tissues and is tightly regulated at the transcriptional and post-transcriptional levels (36, 53). Functional variants in the 3′UTR influence HLA-G expression by altering mRNA stability, microRNA binding, and isoform expression. Six microRNAs (miR-148a, miR-148b, miR-152, miR-133a, miR-628-5p and miR-548q) have been shown to bind sequence specifically to certain SNPs within the HLA-G 3′UTR, thereby down-regulating its expression levels (54, 55). Their regulatory activity appears tissue-dependent (, 56, 57), and many have been detected in renal cell carcinoma and other renal diseases (58, 59). Upon KTx, effect of immunosuppression on miRNA levels has so far not been clarified. The polymorphisms in the 3´UTR region arranged into nine major haplotypes, with HLA-G 3′UTR-1 and -2 being the most common (37, 60). These haplotypes have been associated with differential soluble (s)HLA-G expression. HLA-G 3′UTR-1 was associated with increased sHLA-G levels, whereas HLA-G 3′UTR-2 has low to intermediate HLA-G levels (32). Of note, only soluble HLA-G dimers - especially HLA-G5 - are reported to effectively transmit inhibitory signals via LILRB1 to suppress cytotoxic CD8+ T cell activity by downregulating Granzyme (, 6163). In line with this, increase sHLA-G expression has been associated with improved graft outcomes in different transplant settings (, 50, 6467).

The analysis of both HLA-G 3′UTR donor -recipient haplotypes in transplantation is scarce (). In contrast to our study analysing matched recipient and the donor pairs, the majority of studies focus solely on recipient HLA-G 3′UTR polymorphisms. So far in lung transplantation certain HLA-G 3′UTR recipient haplotypes have been associated with impaired graft outcomes (68). In contrast to HLA-G 3′UTR haplotypes, defined SNPs of the HLA-G 3′UTR region are discussed to be associated with graft rejection or survival: particularly the 14-bp ins/ins genotype and the +3142GG variant, have been associated with a reduced risk of rejection (6975). In contrast, the +3196G allele has been linked to an increased risk of graft loss (69). In our previous study the HLA-G 3´UTR-4 in both donor and recipient has been linked to an increased risk of BK polyomavirus (BKPyV) nephropathy and CMV replication following transplantation (). Both infections are well known to be detrimental in transplant recipients, as BKPyVAN is a major cause of allograft dysfunction and loss, while CMV contributes to morbidity, mortality, and graft injury in kidney transplant recipients (7678). Furthermore, the donor HLA-G 3′UTR-2 was already related with higher KTx rejection rates in univariate analysis (). In the current expanded transplant cohort, this association between donor HLA-G 3´UTR-2 and acute rejection could be confirmed in univariate and in addition in multivariate analyses. Moreover, in the present study, the donor HLA-G 3′UTR-1 also emerged as a risk factor for TCMR in both univariate and multivariate analyses, whereas the donor HLA-G 3′UTR-3 seemed to exhibit a protective effect against TCMR. However, this association was only significant in univariate analysis and not in multivariate analysis. No associations were observed between recipient HLA-G 3′UTR haplotypes and TCMR, and no significant associations were found between donor or recipient HLA-G 3′UTR haplotypes and AMR, emphasizing the regulatory complexity of the HLA-G 3′UTR region in different pathological situations. Thus, our findings suggest that the predominant effect likely arises from locally produced HLA-G within the graft microenvironment rather than from systemic recipient levels.

Interestingly, we found that HLA-G 3′UTR-1 homozygosity was associated with underlying primary diseases such as glomerulonephritis, a group of primary diseases for which genetic predisposition has been implicated. This finding is consistent with previous studies linking HLA-G 3′UTR-1 to other diseases with a known genetic predisposition (). However, regarding the primary outcomes, namely the occurrence of rejection episodes there was no association with the underlying primary diseases in our study cohort (data not shown).

HLA-G modulates local and systemic immune tolerance by engaging inhibitory receptors, primarily LILRB1 and LILRB2 (ILT4), on lymphoid and myeloid cells, thereby suppressing cytotoxicity and inflammation while promoting tolerogenic dendritic cells and regulatory T cells (63, 7982). Thus, impaired LILRB1 function or expression has been implicated in the pathogenesis of several autoimmune diseases (83, 84). While LILRB2 has a high-affinity receptor for HLA-G (85), LILRB1 predominantly interacts with the β2m-associated HLA-G isoforms that are more abundantly expressed in vivo. Given its broader expression on both lymphoid and myeloid cells, LILRB1 is therefore more likely than LILRB2 to mediate HLA-G–dependent immune regulation in transplantation contexts.

Emerging evidence suggests that LILRB1 expression is modulated by polymorphisms in its promoter region, particularly rs10416697 and rs1004443 variants. These variants are believed to influence transcription factor binding and promoter activity, ultimately affecting LILRB1 surface expression on immune cells. The LILRB1-PROMO CG haplotype has been associated with a reduced frequency of NK cells expressing this receptor (71). However, this variation does not affect receptor expression on myeloid effector cells, which constitutively express LILRB1. So far, data regarding LILRB1 expression on T cells was not established. Of note, although LILRB1 is often expressed on T cells, its expression is heterogeneous across CD4+ and CD8+ T cell subsets, with only a small subpopulation of especially CD8+ T cells displaying LILRB1 on their surface. This highlights the complex regulation of HLA-G–LILRB1 interactions.

In has been demonstrated that elevated sHLA-G levels or the presence of extracellular vesicles containing HLA-G increase the frequency of CD8+ LILRB1+ T cells (). Moreover, an increase of LILRB1 mRNA expression was reported in NK and T cells under the influence of elevated sHLA-G levels (86). With respect to transplantation setting, there are no data available how LILRB1 regulation and its expression are altered under the influence of immunosuppressive medication.

In our study, the analysis of donor and recipient LILRB1-PROMO haplotype revealed exclusively the importance of the presence of recipient LILRB1-PROMO CG haplotype to be protective for TCMR, whereas it was an independent risk factor of AMR besides the well-known recipient PRA-positivity. Recent findings underscore the immunological relevance of LILRB1 in the context of kidney transplant rejection. LILRB1 expression is elevated in circulating monocytes of kidney transplant recipients. Notably, myeloid cells isolated from kidney biopsy specimens show up-regulation of LILRB1, LILRB2 and LILRB3 following AMR, highlighting the involvement of this receptor family in the immune response to allografts (87). In line with this, single-cell transcriptomic profiling has demonstrated that recipient-derived monocytes expressing LILRB1 infiltrate the allograft during rejection episodes (87). These data indicate a functional role for LILRB1-expressing myeloid cells in mediating alloresponses. Collectively, these findings highlight LILRB1 as a key immune checkpoint receptor and support its potential as a target for modulating allo-immune activation and impacting transplant outcomes.

LILRB1 is not limited to binding HLA-G, it also recognizes classical HLA class I molecules (HLA-A, -B, -C) and HLA-F but also ligands that derive from various pathogens (especially Cytomegalovirus UL18 and E. coli) (88, 89). These alternative ligands, whose expression may be unregulated during inflammation or rejection, can engage LILRB1 in a conformation-dependent manner, potentially modulating immune outcomes in a context-specific fashion. Beyond ligand diversity, recent studies have highlighted that engagement of the CD47–SIRPα axis can influence the expression of leukocyte Ig-like receptors on myeloid cells, suggesting additional layers of regulation that may intersect with HLA-G/LILRB1 signaling pathways (90, 91). In addition, LILRB1 has also been linked to metabolic pathways, including cholesterol homeostasis and ferroptosis resistance in malignancies such as multiple myeloma, underscoring its functional relevance beyond transplant immunology and highlighting its potential as a target for novel therapeutic strategies (92).

While enhanced LILRB1 expression may support a more tolerogenic or inhibitory profile in antigen-presenting cells, which could suppress early T cell priming and mitigate TCMR, the same polymorphic variant may also influence B cell signaling, differentiation, or survival—thereby facilitating humoral alloimmunity (93). Indeed, LILRB1 signaling on B cells has been shown to impact their responsiveness and function, particularly in inflammatory or alloantigen-rich environments (94). In view of this context LILRB1 promoter haplotype was only important for AMR when present in the recipient, whereas the HLAG 3´UTR of the donor had no impact on AMR rate.

Taken together, our results suggest that the protective effect of HLA-G-LILRB1 interactions in TCMR is primarily driven by donor-recipient specific genetic determinants, whereas in AMR the recipient LILRB1 genetic variants have the most prominent effect. Thus, the HLA-G/LILRB1 axis is an important immune checkpoint not only in tumor biology but also in transplant settings. Interestingly in triple negative breast cancer the HLA-G 3´UTR haplotypes was not associated with disease outcome, however, high levels of sHLA-G combined with the LILRB1-PROMO rs10416697C allele being part of the here defined LILRB1-PROMO CG haplotype, were associated in adverse disease outcome in uni- and multivariate analysis (47).

These findings are clinically significant for enhancing pre-transplant risk stratification, particularly in living donor kidney transplantation cohort, which is the focus of this study. Our data suggest that specific HLA-G and LILRB1 genetic variants may differentially influence rejection risks. Specific donor HLA-G haplotypes and recipient LILRB1-PROMO variants not only influence the risk of rejection independently, but also interact functionally to define graft tolerance or injury. Our results support the integration of HLA-G and LILRB1 genotyping into routine pre-transplant immunogenetic assessment and highlight their potential utility in guiding individualized immunosuppressive management (95). Furthermore, our study is also applicable with regard to therapeutical approaches in terms of using sHLA-G or HLA-G positive extracellular vesicles as modulator of immune responses especially in the difficult to treat situation of AMR (94, 9699). In the setting of stem cell transplantation, the application of HLA-G positive extracellular vesicles drastically improved therapy-refractory graft-versus-host disease (100).

Nevertheless, this study has several limitations making further studies necessary: i) the functional impact of HLA-G 3′UTR haplotypes remains context- and microenvironment-dependent, and is still poorly understood since no mRNA or protein levels were analyzed, ii) the sample size and exclusive focus on living-donor kidney transplantation from one transplant center in Germany may limit the generalizability of the findings to other transplant settings, iii) LILRB1 expression on T cell subsets and other immune effector cells in the donor organ as well as in the recipient periphery is insufficiently characterized, yet may be critical for understanding its role in rejection and tolerance induction and iv) our study did not include direct functional validations.

In conclusion, this study underscores the importance of the donor HLA-G and recipient LILRB1 immune checkpoint axis in shaping transplant outcomes. The results of our study support the potential utility of the HLA-G/LILRB1 immune checkpoint axis. Integrating HLA-G/LILRB1 genotyping into pre-transplant assessment could enhance donor–recipient matching and guide personalized immunosuppressive strategies. Future work should explore therapeutic modulation of this pathway, including the potential use of sHLA-G or HLA-G–positive extracellular vesicles to promote immune regulation and improve long-term graft survival.

Statements

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 Ethics committee of the Medical Faculty of the University of Duisburg-Essen (approval number 12-5312-BO). 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

JH: Funding acquisition, Writing – review & editing, Methodology, Formal analysis, Writing – original draft, Investigation, Data curation, Visualization. SS: Writing – review & editing, Investigation, Methodology. FH: Writing – review & editing. AH: Writing – review & editing. AG: Writing – review & editing. BW: Resources, Writing – review & editing. JR: Visualization, Writing – review & editing. YB: Writing – review & editing, Data curation. NG: Data curation, Writing – review & editing. WP: Data curation, Writing – review & editing. PH: Resources, Writing – review & editing. OW: Funding acquisition, Resources, Writing – review & editing. HR: Methodology, Supervision, Conceptualization, Software, Funding acquisition, Project administration, Visualization, Writing – review & editing, Formal analysis, Validation, Writing – original draft, Data curation. VR: Conceptualization, Writing – original draft, Methodology, Data curation, Validation, Visualization, Software, Funding acquisition, Supervision, Resources, Project administration, Writing – review & editing, Formal analysis.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by ELAN (JH), the Jackstädt Stiftung (HR), and the Rudolf Ackermann Stiftung (OW). The funding organizations had no role in the design of the study, data collection, analysis, interpretation of data, or in writing the manuscript.

Acknowledgments

The authors acknowledge the contributions of all participating patients and the support of ELAN (Essen training programme ‘Laboratory and Science’ for medical students, Medical faculty of Essen), the Jackstädt Stiftung, and the Rudolf Ackermann Stiftung. Special thanks go to the recipients and donors kindly providing their samples. We are grateful for the technical support by the teams of the Department of Nephrology, of the Department of Infectious Diseases, and the Institute for Transfusion Medicine, all University Hospital Essen.

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.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

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.1697839/full#supplementary-material

References

  • 1

    CarosellaEDHoWangYinKYFavierBLeMaoultJ. HLA-G-dependent suppressor cells: Diverse by nature, function, and significance. Hum Immunol. (2008) 69:700–7. doi: 10.1016/j.humimm.2008.08.280

  • 2

    ArnsTAntunesDAAbellaJRRigoMMKavrakiLEGiuliattiSet al. Structural modeling and molecular dynamics of the immune checkpoint molecule HLA-G. Front Immunol. (2020) 11:575076. doi: 10.3389/fimmu.2020.575076

  • 3

    CarosellaEDRouas-FreissNTronik-Le RouxDMoreauPLeMaoultJ. HLA-G: an immune checkpoint molecule. Adv Immunol. (2015) 127:33144. doi: 10.1016/bs.ai.2015.04.001

  • 4

    CarosellaEDFavierBRouas-FreissNMoreauPLemaoultJ. Beyond the increasing complexity of the immunomodulatory HLA-G molecule. Blood. (2008) 111:4862–70. doi: 10.1182/blood-2007-12-127662

  • 5

    CastelliECde AlmeidaBSMunizYCNSilvaNSBPassosMRSSouzaASet al. HLA-G genetic diversity and evolutive aspects in worldwide populations. Sci Rep. (2021) 11:23070. doi: 10.1038/s41598-021-02106-4

  • 6

    RebmannVda Silva NardiFWagnerBHornPA. HLA-G as a tolerogenic molecule in transplantation and pregnancy. J Immunol Res. (2014) 2014:297073. doi: 10.1155/2014/297073

  • 7

    HuntJSLangatDKMcIntireRHMoralesPJ. The role of HLA-G in human pregnancy. Reprod Biol Endocrinol. (2006) 4 Suppl 1:S10. doi: 10.1186/1477-7827-4-S1-S10

  • 8

    Rouas-FreissNLeMaoultJMoreauPDaussetJCarosellaED. HLA-G in transplantation: a relevant molecule for inhibition of graft rejection? Am J Transplant. (2003) 3:11–6. doi: 10.1034/j.1600-6143.2003.30103.x

  • 9

    RohnHSchwichETomoya MichitaRSchrammSDolffSGacklerAet al. HLA-G 3’ untranslated region gene variants are promising prognostic factors for BK polyomavirus replication and acute rejection after living-donor kidney transplant. Hum Immunol. (2020) 81:141–6. doi: 10.1016/j.humimm.2019.09.011

  • 10

    RizzoRBortolottiDBolzaniSFainardiE. HLA-G molecules in autoimmune diseases and infections. Front Immunol. (2014) 5:592. doi: 10.3389/fimmu.2014.00592

  • 11

    SeligerB. Role of HLA-G in tumors and upon COVID-19 infection. Hum Immunol. (2024) 85:110792. doi: 10.1016/j.humimm.2024.110792

  • 12

    RohnHElischerFLarbigLJansenSSchrammSOtteMet al. Influence of HLA-G 3’ Untranslated region haplotypes and SNP + 3422 gene variants as host genetic factors on the outcomes of SARS-CoV-2 infection during acute and post-acute phases in a German cohort. HLA. (2024) 104:e15799. doi: 10.1111/tan.15799

  • 13

    RohnHSchrammSPansikakiKJansenSHendriksCPlatteMet al. Peripheral HLA-G/ILT-2 immune checkpoint axis in acute and convalescent COVID-19 patients. Hum Immunol. (2023) 84:393400. doi: 10.1016/j.humimm.2023.03.002

  • 14

    RebmannVSchwichEMichitaRTGruntkemeierLBittnerAKRohnHet al. Systematic evaluation of HLA-G 3’Untranslated region variants in locally advanced, non-metastatic breast cancer patients: UTR-1, 2 or UTR-4 are predictors for therapy and disease outcome. Front Immunol. (2021) 12:817132. doi: 10.3389/fimmu.2021.817132

  • 15

    SchwichERebmannVMichitaRTRohnHVonckenJWHornPAet al. HLA-G 3’ untranslated region variants +3187G/G, +3196G/G and +3035T define diametrical clinical status and disease outcome in epithelial ovarian cancer. Sci Rep. (2019) 9:5407. doi: 10.1038/s41598-019-41900-z

  • 16

    TizaouiKJalouliMBoujelbeneNHarrathAHOuzariHIRizzoRet al. The relationship of 3’UTR HLA-G14-bp insertion/deletion and +3142 C/G polymorphisms and soluble HLA-G expression with gynecological cancers: An updated meta-analysis. Immun Inflammation Dis. (2022) 10:e645. doi: 10.1002/iid3.645

  • 17

    AlmeidaRSFerreiraMLBSononPCordeiroMTSadissouIDinizGTNet al. Cytokines and soluble HLA-G levels in the acute and recovery phases of arbovirus-infected Brazilian patients exhibiting neurological complications. Front Immunol. (2021) 12:582935. doi: 10.3389/fimmu.2021.582935

  • 18

    CarosellaEDGregoriSTronik-Le RouxD. HLA-G/LILRBs: A cancer immunotherapy challenge. Trends Cancer. (2021) 7:389–92. doi: 10.1016/j.trecan.2021.01.004

  • 19

    JacquierALambertTDelattreJFDjouadouMVerineJDumontCet al. Tumor infiltrating and peripheral CD4(+)ILT2(+) T cells are a cytotoxic subset selectively inhibited by HLA-G in clear cell renal cell carcinoma patients. Cancer Lett. (2021) 519:105–16. doi: 10.1016/j.canlet.2021.06.018

  • 20

    Le RondSLe MaoultJCreputCMenierCDeschampsMLe FriecGet al. Alloreactive CD4+ and CD8+ T cells express the immunotolerant HLA-G molecule in mixed lymphocyte reactions: in vivo implications in transplanted patients. Eur J Immunol. (2004) 34:649–60. doi: 10.1002/eji.200324266

  • 21

    LeMaoultJRouas-FreissNCarosellaED. Immuno-tolerogenic functions of HLA-G: relevance in transplantation and oncology. Autoimmun Rev. (2005) 4:503–9. doi: 10.1016/j.autrev.2005.04.006

  • 22

    DumontCJacquierAVerineJNoelFGoujonAWuCLet al. CD8(+)PD-1(-)ILT2(+) T cells are an intratumoral cytotoxic population selectively inhibited by the immune-checkpoint HLA-G. Cancer Immunol Res. (2019) 7:1619–32. doi: 10.1158/2326-6066.CIR-18-0764

  • 23

    FavierBLemaoultJLesportECarosellaED. ILT2/HLA-G interaction impairs NK-cell functions through the inhibition of the late but not the early events of the NK-cell activating synapse. FASEB J. (2010) 24:689–99. doi: 10.1096/fj.09-135194

  • 24

    JacquierADumontCCarosellaEDRouas-FreissNLeMaoultJ. Cytometry-based analysis of HLA-G functions according to ILT2 expression. Hum Immunol. (2020) 81:168–77. doi: 10.1016/j.humimm.2020.02.001

  • 25

    NajiAMenierCMakiGCarosellaEDRouas-FreissN. Neoplastic B-cell growth is impaired by HLA-G/ILT2 interaction. Leukemia. (2012) 26:1889–92. doi: 10.1038/leu.2012.62

  • 26

    Khalil-DaherIRiteauBMenierCSedlikCPaulPDaussetJet al. Role of HLA-G versus HLA-E on NK function: HLA-G is able to inhibit NK cytolysis by itself. J Reprod Immunol. (1999) 43:175–82. doi: 10.1016/S0165-0378(99)00032-7

  • 27

    SchwichEHoGTLeMaoultJBade-DodingCCarosellaEDHornPAet al. Soluble HLA-G and HLA-G bearing extracellular vesicles affect ILT-2 positive and ILT-2 negative CD8 T cells complementary. Front Immunol. (2020) 11:2046. doi: 10.3389/fimmu.2020.02046

  • 28

    LepinEJBastinJMAllanDSRoncadorGBraudVMMasonDYet al. Functional characterization of HLA-F and binding of HLA-F tetramers to ILT2 and ILT4 receptors. Eur J Immunol. (2000) 30:3552–61. doi: 10.1002/1521-4141(200012)30:12<3552::AID-IMMU3552>3.0.CO;2-L

  • 29

    WillcoxBEThomasLMBjorkmanPJ. Crystal structure of HLA-A2 bound to LIR-1, a host and viral major histocompatibility complex receptor. Nat Immunol. (2003) 4:913–9. doi: 10.1038/ni961

  • 30

    AlvesCCArnsTOliveiraMLMoreauPAntunesDACastelliECet al. Computational and atomistic studies applied to the understanding of the structural and behavioral features of the immune checkpoint HLA-G molecule and gene. Hum Immunol. (2023) 84:374–83. doi: 10.1016/j.humimm.2023.01.004

  • 31

    Arnaiz-VillenaAJuarezISuarez-TrujilloFLopez-NaresAVaqueroCPalacio-GruberJet al. HLA-G: Function, polymorphisms and pathology. Int J Immunogenet. (2021) 48:172–92. doi: 10.1111/iji.12513

  • 32

    PorasIYaghiLMartelli-PalominoGMendes-JuniorCTMunizYCCagninNFet al. Haplotypes of the HLA-G 3’ Untranslated region respond to endogenous factors of HLA-G+ and HLA-G- cell lines differentially. PloS One. (2017) 12:e0169032. doi: 10.1371/journal.pone.0169032

  • 33

    RousseauPLe DiscordeMMouillotGMarcouCCarosellaEDMoreauP. The 14 bp deletion-insertion polymorphism in the 3’ UT region of the HLA-G gene influences HLA-G mRNA stability. Hum Immunol. (2003) 64:1005–10. doi: 10.1016/j.humimm.2003.08.347

  • 34

    TanZRandallGFanJCamoretti-MercadoBBrockman-SchneiderRPanLet al. Allele-specific targeting of microRNAs to HLA-G and risk of asthma. Am J Hum Genet. (2007) 81:829–34. doi: 10.1086/521200

  • 35

    Martelli-PalominoGPancottoJAMunizYCMendes-JuniorCTCastelliECMassaroJDet al. Polymorphic sites at the 3’ untranslated region of the HLA-G gene are associated with differential hla-g soluble levels in the Brazilian and French population. PloS One. (2013) 8:e71742. doi: 10.1371/journal.pone.0071742

  • 36

    CastelliECVeiga-CastelliLCYaghiLMoreauPDonadiEA. Transcriptional and posttranscriptional regulations of the HLA-G gene. J Immunol Res. (2014) 2014:734068. doi: 10.1155/2014/734068

  • 37

    CastelliECMendes-JuniorCTDeghaideNHde AlbuquerqueRSMunizYCSimoesRTet al. The genetic structure of 3’untranslated region of the HLA-G gene: polymorphisms and haplotypes. Genes Immun. (2010) 11:134–41. doi: 10.1038/gene.2009.74

  • 38

    YuKDavidsonCLWojtowiczALisboaLWangTAiroAMet al. LILRB1 polymorphisms influence posttransplant HCMV susceptibility and ligand interactions. J Clin Invest. (2018) 128:1523–37. doi: 10.1172/JCI96174

  • 39

    DavidsonCLLiNLBurshtynDN. LILRB1 polymorphism and surface phenotypes of natural killer cells. Hum Immunol. (2010) 71:942–9. doi: 10.1016/j.humimm.2010.06.015

  • 40

    RacusenLCSolezKColvinRBBonsibSMCastroMCCavalloTet al. The Banff 97 working classification of renal allograft pathology. Kidney Int. (1999) 55:713–23. doi: 10.1046/j.1523-1755.1999.00299.x

  • 41

    ChandranSMannonRB. T cell-mediated rejection in kidney transplant recipients: The end(point) is also the beginning. Am J Transpl. (2022) 22:683–4. doi: 10.1111/ajt.16964

  • 42

    MayerKASchrezenmeierEDieboldMHalloranPFSchatzlMSchranzSet al. A randomized phase 2 trial of felzartamab in antibody-mediated rejection. N Engl J Med. (2024) 391:122–32. doi: 10.1056/NEJMoa2400763

  • 43

    HaasMLoupyALefaucheurCRoufosseCGlotzDSeronDet al. The Banff 2017 Kidney Meeting Report: Revised diagnostic criteria for chronic active T cell-mediated rejection, antibody-mediated rejection, and prospects for integrative endpoints for next-generation clinical trials. Am J Transpl. (2018) 18:293307. doi: 10.1111/ajt.14625

  • 44

    LoupyAHaasMSolezKRacusenLGlotzDSeronDet al. The Banff 2015 kidney meeting report: current challenges in rejection classification and prospects for adopting molecular pathology. Am J Transpl. (2017) 17:2841. doi: 10.1111/ajt.14107

  • 45

    HaasMSisBRacusenLCSolezKGlotzDColvinRBet al. Banff 2013 meeting report: inclusion of c4d-negative antibody-mediated rejection and antibody-associated arterial lesions. Am J Transpl. (2014) 14:272–83. doi: 10.1111/ajt.12590

  • 46

    RacusenLCColvinRBSolezKMihatschMJHalloranPFCampbellPMet al. Antibody-mediated rejection criteria - an addition to the Banff 97 classification of renal allograft rejection. Am J Transpl. (2003) 3:708–14. doi: 10.1034/j.1600-6143.2003.00072.x

  • 47

    HoffmannOWormlandSBittnerAKHolzenbeinJSchwichESchrammSet al. Elevated sHLA-G plasma levels post chemotherapy combined with ILT-2 rs10416697C allele status of the sHLA-G-related receptor predict poorest disease outcome in early triple-negative breast cancer patients. Front Immunol. (2023) 14:1188030. doi: 10.3389/fimmu.2023.1188030

  • 48

    KnollGAFortinMCGillJGrimshawJMHartellDPKarnabiPet al. Measuring quality in living donation and kidney transplantation: moving beyond survival metrics. Kidney Int. (2020) 98:860–9. doi: 10.1016/j.kint.2020.07.014

  • 49

    DavisCLDelmonicoFL. Living-donor kidney transplantation: a review of the current practices for the live donor. J Am Soc Nephrol. (2005) 16:2098–110. doi: 10.1681/ASN.2004100824

  • 50

    BrugiereODreyfussDGuiletRRongSHirschiSRenaud-PicardBet al. Circulating vesicular-bound HLA-G as noninvasive predictive biomarker of CLAD after lung transplantation. Transplantation. (2025) 109:736–45. doi: 10.1097/TP.0000000000005175

  • 51

    AmodioGComiMTomasoniDGianoliniMERizzoRLeMaoultJet al. HLA-G expression levels influence the tolerogenic activity of human DC-10. Haematologica. (2015) 100:548–57. doi: 10.3324/haematol.2014.113803

  • 52

    HviidTVRizzoRChristiansenOBMelchiorriLLindhardABaricordiOR. HLA-G and IL-10 in serum in relation to HLA-G genotype and polymorphisms. Immunogenetics. (2004) 56:135–41. doi: 10.1007/s00251-004-0673-2

  • 53

    AlegreERizzoRBortolottiDFernandez-LandazuriSFainardiEGonzalezA. Some basic aspects of HLA-G biology. J Immunol Res. (2014) 2014:657625. doi: 10.1155/2014/657625

  • 54

    PortoIOMendes-JuniorCTFelicioLPGeorgRCMoreauPDonadiEAet al. MicroRNAs targeting the immunomodulatory HLA-G gene: a new survey searching for microRNAs with potential to regulate HLA-G. Mol Immunol. (2015) 65:230–41. doi: 10.1016/j.molimm.2015.01.030

  • 55

    CastelliECMoreauPOya e ChiromatzoAMendes-JuniorCTVeiga-CastelliLCYaghiLet al. In silico analysis of microRNAS targeting the HLA-G 3’ untranslated region alleles and haplotypes. Hum Immunol. (2009) 70:1020–5. doi: 10.1016/j.humimm.2009.07.028

  • 56

    SeligerB. Role of microRNAs on HLA-G expression in human tumors. Hum Immunol. (2016) 77:760–3. doi: 10.1016/j.humimm.2016.04.006

  • 57

    Jasinski-BergnerSSchmiedelDMandelboimOSeligerB. Role of HLA-G in viral infections. Front Immunol. (2022) 13:826074. doi: 10.3389/fimmu.2022.826074

  • 58

    Jasinski-BergnerSRechesAStoehrCMassaCGonschorekEHuettelmaierSet al. Identification of novel microRNAs regulating HLA-G expression and investigating their clinical relevance in renal cell carcinoma. Oncotarget. (2016) 7:26866–78. doi: 10.18632/oncotarget.8567

  • 59

    MukhadiSHullRMbitaZDlaminiZ. The role of microRNAs in kidney disease. Noncoding RNA. (2015) 1:192221. doi: 10.3390/ncrna1030192

  • 60

    CastelliECRamalhoJPortoIOLimaTHFelicioLPSabbaghAet al. Insights into HLA-G genetics provided by worldwide haplotype diversity. Front Immunol. (2014) 5:476. doi: 10.3389/fimmu.2014.00476

  • 61

    BahriRHirschFJosseARouas-FreissNBidereNVasquezAet al. Soluble HLA-G inhibits cell cycle progression in human alloreactive T lymphocytes. J Immunol. (2006) 176:1331–9. doi: 10.4049/jimmunol.176.3.1331

  • 62

    AjithAPortik-DobosVNguyen-LefebvreATCallawayCHoruzskoDDKapoorRet al. HLA-G dimer targets Granzyme B pathway to prolong human renal allograft survival. FASEB J. (2019) 33:5220–36. doi: 10.1096/fj.201802017R

  • 63

    HoWangYinKYLoustauMWuJAlegreEDaouyaMCaumartinJet al. Multimeric structures of HLA-G isoforms function through differential binding to LILRB receptors. Cell Mol Life Sci. (2012) 69:4041–9. doi: 10.1007/s00018-012-1069-3

  • 64

    CreputCLe FriecGBahriRAmiotLCharpentierBCarosellaEet al. Detection of HLA-G in serum and graft biopsy associated with fewer acute rejections following combined liver-kidney transplantation: possible implications for monitoring patients. Hum Immunol. (2003) 64:1033–8. doi: 10.1016/j.humimm.2003.08.356

  • 65

    AjithAMamouniKMusaAHoruzskoDDGaniIMulloyLLet al. IL-10-producing memory B regulatory cells as a novel target for HLA-G to prolong human kidney allograft survival. Hum Immunol. (2023) 84:366–73. doi: 10.1016/j.humimm.2023.03.003

  • 66

    BortolottiDGentiliVRotolaAPotenaLRizzoR. Soluble HLA-G pre-transplant levels to identify the risk for development of infection in heart transplant recipients. Hum Immunol. (2020) 81:147–50. doi: 10.1016/j.humimm.2019.10.003

  • 67

    BrugiereOThabutGKrawice-RadanneIRizzoRDauriatGDanelCet al. Role of HLA-G as a predictive marker of low risk of chronic rejection in lung transplant recipients: a clinical prospective study. Am J Transpl. (2015) 15:461–71. doi: 10.1111/ajt.12977

  • 68

    Di CristofaroJReynaud-GaubertMCarliniFRoubertouxPLoundouABasireAet al. HLA-G*01:04 approximately UTR3 recipient correlates with lower survival and higher frequency of chronic rejection after lung transplantation. Am J Transpl. (2015) 15:2413–20. doi: 10.1111/ajt.13305

  • 69

    Ciliao AlvesDCde Oliveira CrispimJCCastelliECMendes-JuniorCTDeghaideNHBarros SilvaGEet al. Human leukocyte antigen-G 3’ untranslated region polymorphisms are associated with better kidney allograft acceptance. Hum Immunol. (2012) 73:52–9. doi: 10.1016/j.humimm.2011.10.007

  • 70

    KangSWOhEChoWKimMParkEJKwackKHet al. HLA-G 14bp Ins/del polymorphism in the 3’UTR region and acute rejection in kidney transplant recipients: an updated meta-analysis. Med (Kaunas). (2021) 57. doi: 10.3390/medicina57101007

  • 71

    DarbasSYilmazVTKocakHKisaogluADemiryilmazIAydinliBet al. New markers for predictions of acute and chronic rejection and graft outcomes in kidney transplant recipients; HLA-G gene 3’UTR 14 bp polymorphism and sHLA-G. Gene. (2021) 790:145712. doi: 10.1016/j.gene.2021.145712

  • 72

    DurmanovaVBandzuchovaHZilinskaZTirpakovaJKubaDBucMet al. Association of HLA-G polymorphisms in the 3’UTR region and soluble HLA-G with kidney graft outcome. Immunol Invest. (2019) 48:644–58. doi: 10.1080/08820139.2019.1610888

  • 73

    GuberinaHTomoya MichitaRDolffSBienholzATrillingMHeinemannFMet al. Recipient HLA-G +3142 CC genotype and concentrations of soluble HLA-G impact on occurrence of CMV infection after living-donor kidney transplantation. Int J Mol Sci. (2017) 18. doi: 10.3390/ijms18112338

  • 74

    PiancatelliDMaccaroneDLiberatoreGParzaneseIClementeKAzzaroneRet al. HLA-G 14-bp insertion/deletion polymorphism in kidney transplant patients with metabolic complications. Transplant Proc. (2009) 41:1187–8. doi: 10.1016/j.transproceed.2009.03.028

  • 75

    HauerVRistiMMirandaBLMda SilvaJSCidralALPozziCMet al. The association of HLA-G polymorphisms and the synergistic effect of sMICA and sHLA-G with chronic kidney disease and allograft acceptance. PloS One. (2019) 14:e0212750. doi: 10.1371/journal.pone.0212750

  • 76

    FishmanJAIssaNC. Infection in organ transplantation: risk factors and evolving patterns of infection. Infect Dis Clin North Am. (2010) 24:273–83. doi: 10.1016/j.idc.2010.01.005

  • 77

    RohnHMichitaRTSchrammSDolffSGacklerAKorthJet al. HLA-E polymorphism determines susceptibility to BK virus nephropathy after living-donor kidney transplant. Cells. (2019) 8. doi: 10.3390/cells8080847

  • 78

    GuberinaHda Silva NardiFMichitaRTDolffSBienholzAHeinemannFMet al. Susceptibility of HLA-E*01:03 allele carriers to develop cytomegalovirus replication after living-donor kidney transplantation. J Infect Dis. (2018) 217:1918–22. doi: 10.1093/infdis/jix638

  • 79

    NajiALe RondSDurrbachAKrawice-RadanneICreputCDaouyaMet al. CD3+CD4low and CD3+CD8low are induced by HLA-G: novel human peripheral blood suppressor T-cell subsets involved in transplant acceptance. Blood. (2007) 110:3936–48. doi: 10.1182/blood-2007-04-083139

  • 80

    Rouas-FreissNLeMaoultJVerineJTronik-Le RouxDCulineSHennequinCet al. Intratumor heterogeneity of immune checkpoints in primary renal cell cancer: Focus on HLA-G/ILT2/ILT4. Oncoimmunology. (2017) 6:e1342023. doi: 10.1080/2162402X.2017.1342023

  • 81

    WuCLSvendsenSGRiviereADesgrandchampsFCarosellaEDLeMaoultJ. Multiplex bead-based immunoassay for the free soluble forms of the HLA-G receptors, ILT2 and ILT4. Hum Immunol. (2016) 77:720–6. doi: 10.1016/j.humimm.2016.01.017

  • 82

    GarciaMPalmaMBVerineJMiriukaSIndaAMErrecaldeALet al. The immune-checkpoint HLA-G/ILT4 is involved in the regulation of VEGF expression in clear cell renal cell carcinoma. BMC Cancer. (2020) 20:624. doi: 10.1186/s12885-020-07113-8

  • 83

    Monsivais-UrendaANino-MorenoPAbud-MendozaCBarandaLLayseca-EspinosaELopez-BotetMet al. Analysis of expression and function of the inhibitory receptor ILT2 (CD85j/LILRB1/LIR-1) in peripheral blood mononuclear cells from patients with systemic lupus erythematosus (SLE). J Autoimmun. (2007) 29:97105. doi: 10.1016/j.jaut.2007.05.003

  • 84

    VeitTDChiesJASwitalaMWagnerBHornPABusattoMet al. The paradox of high availability and low recognition of soluble HLA-G by LILRB1 receptor in rheumatoid arthritis patients. PloS One. (2015) 10:e0123838. doi: 10.1371/journal.pone.0123838

  • 85

    ShiroishiMTsumotoKAmanoKShirakiharaYColonnaMBraudVMet al. Human inhibitory receptors Ig-like transcript 2 (ILT2) and ILT4 compete with CD8 for MHC class I binding and bind preferentially to HLA-G. Proc Natl Acad Sci U S A. (2003) 100:8856–61. doi: 10.1073/pnas.1431057100

  • 86

    LeMaoultJZafaranlooKLe DanffCCarosellaED. HLA-G up-regulates ILT2, ILT3, ILT4, and KIR2DL4 in antigen presenting cells, NK cells, and T cells. FASEB J. (2005) 19:662–4. doi: 10.1096/fj.04-1617fje

  • 87

    LamartheeBGenetCCattinFDangerRGiralMBrouardSet al. Single-cell mapping of leukocyte immunoglobulin-like receptors in kidney transplant rejection. Front Transpl. (2022) 1:952785. doi: 10.3389/frtra.2022.952785

  • 88

    ChapmanTLHeikemanAPBjorkmanPJ. The inhibitory receptor LIR-1 uses a common binding interaction to recognize class I MHC molecules and the viral homolog UL18. Immunity. (1999) 11:603–13. doi: 10.1016/S1074-7613(00)80135-1

  • 89

    Redondo-GarciaSBarrittCPapagregoriouCYeboahMFrendeusBCraggMSet al. Human leukocyte immunoglobulin-like receptors in health and disease. Front Immunol. (2023) 14:1282874. doi: 10.3389/fimmu.2023.1282874

  • 90

    DaiHLanPZhaoDAbou-DayaKLiuWChenWet al. PIRs mediate innate myeloid cell memory to nonself MHC molecules. Science. (2020) 368:1122–7. doi: 10.1126/science.aax4040

  • 91

    ZhaoDDaiHMacedoCSandersSMEliasCFridayAJet al. Donor-recipient mismatch at the SIRPA locus adversely affects kidney allograft outcomes. Sci Transl Med. (2025) 17:eady1135. doi: 10.1126/scitranslmed.ady1135

  • 92

    XianMWangQXiaoLZhongLXiongWYeLet al. Leukocyte immunoglobulin-like receptor B1 (LILRB1) protects human multiple myeloma cells from ferroptosis by maintaining cholesterol homeostasis. Nat Commun. (2024) 15:5767. doi: 10.1038/s41467-024-50073-x

  • 93

    RohnHLangCSchrammSHeinemannFMTrillingMGacklerAet al. Effect of HLA-G5 immune checkpoint molecule on the expression of ILT-2, CD27, and CD38 in splenic B cells. J Immunol Res. (2022) 2022:4829227. doi: 10.1155/2022/4829227

  • 94

    NajiAMenierCMorandiFAgaugueSMakiGFerrettiEet al. Binding of HLA-G to ITIM-bearing Ig-like transcript 2 receptor suppresses B cell responses. J Immunol. (2014) 192:1536–46. doi: 10.4049/jimmunol.1300438

  • 95

    BahriRNajiAMenierCCharpentierBCarosellaEDRouas-FreissNet al. Dendritic cells secrete the immunosuppressive HLA-G molecule upon CTLA4-Ig treatment: implication in human renal transplant acceptance. J Immunol. (2009) 183:7054–62. doi: 10.4049/jimmunol.0803054

  • 96

    LeMaoultJDaouyaMWuJLoustauMHoruzskoACarosellaED. Synthetic HLA-G proteins for therapeutic use in transplantation. FASEB J. (2013) 27:3643–51. doi: 10.1096/fj.13-228247

  • 97

    SeligerB. The non-classical antigens of HLA-G and HLA-E as diagnostic and prognostic biomarkers and as therapeutic targets in transplantation and tumors. Clin Transpl. (2013) 465–72.

  • 98

    von WebskyMWKitamuraKLudwig-PortugallIKurtsCvon LaffertMLeMaoultJet al. Recombinant HLA-G as tolerogenic immunomodulant in experimental small bowel transplantation. PloS One. (2016) 11:e0158907. doi: 10.1371/journal.pone.0158907

  • 99

    RebmannVKonigLNardi FdaSWagnerBManvailerLFHornPA. The potential of HLA-G-bearing extracellular vesicles as a future element in HLA-G immune biology. Front Immunol. (2016) 7:173. doi: 10.3389/fimmu.2016.00173

  • 100

    KordelasLRebmannVLudwigAKRadtkeSRuesingJDoeppnerTRet al. MSC-derived exosomes: a novel tool to treat therapy-refractory graft-versus-host disease. Leukemia. (2014) 28:970–3. doi: 10.1038/leu.2014.41

Summary

Keywords

HLA-G 3´UTR-1, HLA-G 3´UTR-2, LILRB1, ILT-2, kidney transplantation, acute cellular rejection, humoral rejection

Citation

Hölzenbein J, Schramm S, Heinemann FM, Heinold A, Gäckler A, Reinold J, Wilde B, Busch Y, Gruenen N, Peter W, Horn PA, Witzke O, Rohn H and Rebmann V (2026) Genetic variants of the HLA-G/LILRB1 ligand-receptor axis in donors or recipients are prognostic covariates for rejection after living kidney transplantation. Front. Immunol. 16:1697839. doi: 10.3389/fimmu.2025.1697839

Received

02 September 2025

Accepted

20 October 2025

Published

05 January 2026

Volume

16 - 2025

Edited by

Jules Russick, Commissariat à l’Energie Atomique et aux Energies Alternatives (CEA), France

Reviewed by

Baptiste Lamarthée, Université de Franche-Comté, EFS, INSERM, UMR RIGHT, France

Arwa Kamoun, Hedi Chaker University Hospital and Faculty of Medicine of Sfax, Tunisia

Updates

Copyright

*Correspondence: Hana Rohn, ; Vera Rebmann,

† These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics