ORIGINAL RESEARCH article

Front. Genet., 31 July 2024

Sec. Immunogenetics

Volume 15 - 2024 | https://doi.org/10.3389/fgene.2024.1442938

Combined maternal KIR2DL4 and fetal HLA-G polymorphisms were associated with preeclampsia in a Han Chinese population

  • 1. College of Basic Medical Sciences, Kunming Medical University, Kunming, China

  • 2. Institute of Biomedical Engineering, Kunming Medical University, Kunming, China

  • 3. Prenatal Diagnosis Center, The first Affiliated Hospital of Kunming Medical University, Kunming, China

  • 4. Precision Medicine Center, The Affiliated Hospital of Yunnan University (The Second People’s Hospital of Yunnan Province), Kunming, China

  • 5. Clinical Laboratory, Jiangxi Province Hospital of Integrated Chinese and Western Medicine, Nanchang, China

  • 6. Obstetrics Department, Yan’an Hospital of Kunming City, Kunming, China

  • 7. Clinical Laboratory, Yan’an Hospital of Kunming City, Kunming, China

Abstract

Preeclampsia is the main cause of maternal and infant mortality and morbidity during pregnancy. Killer cell immunoglobulin-like receptor 2DL4 (KIR2DL4) and human leukocyte antigen G (HLA-G) play crucial roles in immune tolerance at the maternal-fetal interface. In this case‒control study, 154 maternal–fetal pairs were recruited, including 74 pairs with preeclampsia (56 of 74 pairs from family triads) and 80 pairs with a normal pregnancy (78 of 80 pairs from family triads). SNaPshot technology was used to detect genetic polymorphisms for 7 TagSNPs in the KIR2DL4 and HLA-G genes. Among the fetal HLA-G gene polymorphisms, rs9380142 (A vs. G: OR = 2.802, 95% CI = 1.761–4.458) and rs1063320 (G vs. C: OR = 1.807, 95% CI = 1.144–2.852) differed between the preeclampsia group and the control group. The transmission disequilibrium test (TDT) suggested that the differences in the rs9380142G/A polymorphism in foetuses between preeclampsia triads and control triads were due to differences in transmission from the parents (P = 0.001). There was no significant difference in the distribution of maternal KIR2DL4 alleles or genotype frequency between the preeclampsia group and the control group. Gene‒gene interaction analysis revealed that the combined genotypes of maternal rs649216-CC and fetal rs9380142-GG, maternal rs1051456-CG/GG and fetal rs9380142-GG, maternal rs34785252-CC and fetal rs9380142-AA/GA, and maternal rs34785252-CC/AA and fetal rs9380142-GG were associated with a significantly lower risk of preeclampsia. Therefore, this study suggested that the combination of maternal KIR2DL4 and fetal HLA-G polymorphisms was associated with preeclampsia in a Han Chinese population.

1 Introduction

Preeclampsia is a pregnancy-specific disorder characterized by hypertension and proteinuria after 20 weeks of gestation (Gestational Hypertension and Preeclampsia 2020). It is often accompanied by edema, fetal growth restriction, and other complications that pose serious health risks to both the mother and the fetus, and is a common obstetrical complication with an incidence of 5%–10% (Al-Jameil et al., 2014). According to statistics, preeclampsia causes 3.9% of maternal deaths and 10.7% of perinatal fetal deaths (Colbern, et al., 1994), seriously affecting the health of mothers and infants. Currently, most scholars believe that disruption of the immune balance at the maternal-fetal interface in the gestational state is the main cause of preeclampsia (Mukherjee, et al., 2023).

The maternal-fetal interface consists mainly of trophoblasts from the embryo and decidua from the mother. Placental trophoblasts do not express classical human leukocyte antigen (HLAIa) or HLAII but rather express nonclassical class I molecules, i.e., HLA-C, HLA-E and HLA-G (Colbern, et al., 1994; Ferreira, et al., 2017). The HLA-G gene is located on chromosome 6 and contains 8 exons and 7 introns. It has a strict tissue distribution and is expressed at high levels in extravillous trophoblasts (EVTs) (Ferreira, et al., 2017). Low expression and/or secretion of HLA-G in the placenta is associated with preeclampsia (Emmer, et al., 2004; Tang, et al., 2015). The HLA-G gene is polymorphic, particularly in the 3′-untranslated region (UTR), which contains several regulatory elements, including AU-rich elements, polyadenosine signals and signals regulating mRNA spatial and temporal expression (Reches, et al., 2020). HLA-G gene polymorphisms may primarily affect the stability of HLA-G mRNA, which in turn affects HLA-G expression. HLA-G exists in both soluble and membrane-bound forms, and its interactions with maternal immune cells play a critical role in maintaining immune tolerance during pregnancy (Mao, et al., 2023).

The decidua of early pregnancy has a large number of immune cells (Kaminski, et al., 2019), 70% of which are uterine natural killer (uNK) cells (Parham, 2004), which have low cytotoxicity (Moffett and Shreeve, 2015). uNK cells express killer cell immunoglobulin-like receptors (KIRs) on their surface. The extracellular domain of the KIR protein binds to its ligand, the HLA-I molecule. The KIR gene family is located on chromosome 19 and consists of 15 genes and 2 pseudogenes (Ritari, et al., 2022). KIR gene polymorphisms are closely related to their expression and the functional activity of their products. KIR2DL4, a member of the KIR family, is expressed on the surface of all NK cells, and its ligand is HLA-G (Colucci, 2017).

EVT cells are highly invasive uterus cells that mediate the remodeling of the uterine spiral arteries to ensure the normal exchange of material between the mother and foetus, which is essential for normal embryonic development. EVT invasion is regulated mainly by uNK cells (Kieckbusch, et al., 2014; Colucci, 2017). Maternal KIRs bind to fetal HLAs and regulate uNK cell function through the KIR-HLA receptor‒ligand pathway (Kieckbusch, et al., 2014; Colucci, 2017). When abnormal uNK cell function leads to insufficient EVT invasion into the uterus, the placenta is poorly formed, triggering preeclampsia (Hiby, et al., 2004; Nakimuli, et al., 2015). Thus, maternal and fetal factors play an equally important role in the immune imbalance between uNK cells and EVTs that triggers preeclampsia. Most of the current studies only consider genetic factors, molecular expression or associated microenvironmental changes in mothers, ignoring the influence of fetal factors. In this study, we investigated the polymorphisms of maternal KIR2DL4 and fetal HLA-G genes and assessed the correlation between different maternal-fetal genotype combinations and preeclampsia to identify preeclampsia susceptibility genes and to aid in the early diagnosis and prediction of this disease.

2 Materials and methods

2.1 Subjects

In this case‒control study, 154 maternal–fetal pairs were recruited, including 74 pairs with preeclampsia (56 of 74 pairs from family triads) and 80 pairs with a normal pregnancy (78 of 80 pairs from family triads). All samples were collected from the outpatient department of prenatal diagnosis at the First Affiliated Hospital of Kunming Medical University and Yan’an Hospital of Kunming city. These samples were used for preeclampsia evaluation. All participants were enrolled between September 2020 and June 2023. Preeclampsia was diagnosed using the criteria of the International Association for Pregnancy-Induced Hypertension (ISSHP). The diagnostic criteria were as follows: systolic blood pressure (SBP) and diastolic blood pressure (DBP) ≥140/90 mmHg present after 20 weeks of gestation, measured at least twice or at least 4 h apart, and albuminuria ≥300 mg/24 h. The control group was from the same center and consisted of pregnant volunteer women with normal blood pressure, and no history of preeclampsia. The exclusion criteria were metabolic abnormalities, chronic cardiovascular diseases, liver, kidney and endocrine diseases, Rh blood group incompatibility, autoimmune diseases and chromosome abnormalities. All participants signed an informed consent form for the collection of maternal peripheral blood and fetal umbilical vein blood samples. The study was approved by the Medical Ethics Committee of Kunming Medical University (Approval No. KMMU2020MEC027).

2.2 Sample collection and DNA extraction

During the 20th week of pregnancy, prenatal amniocentesis was performed by a professional obstetrician. In order to minimize risks to the mother and fetus, amniotic fluid is collected under sterile conditions and with ultrasound guidance. This procedure collected amniotic fluid for genetic analysis and scientific research. Peripheral blood samples were collected from parents for genotyping purposes. According to the manufacturer’s instructions, DNA was extracted from peripheral blood using the DNeasy Blood and Tissue Kit (QIAGEN, Germany). The obtained DNA was stored in a freezer at −20°C until analysis.

2.3 TagSNP selection and analysis

The necessary genotype data were obtained from the 1,000 Genomes Browser (https://www.ncbi.nlm.nih.gov/variation/tools/1000genomes/). TagSNPs were selected using HaploView 4.2 software with the following settings: upstream and downstream, 2K; minor allele frequency (MAF) ≥0.05; and linkage disequilibrium threshold (R2) ≥0.8. Three TagSNPs for the KIR2DL4 gene and four TagSNPs for the HLA-G gene were identified (Supplementary Table S1). The SNaPshot method was employed to analyze SNP polymorphisms. Details of the primers used are provided in Supplementary Table S2.

2.4 Statistical analysis

SPSS.229 (IBM) was used for the statistical analysis. For continuous variables, the Shapiro‒Wilk test was used for the normality test. Normally distributed continuous variables were expressed as the mean ± standard deviation, and nonnormally distributed variables were expressed as the median and quartile range (IQR). Classification variables were expressed as numbers or percentages. For each SNP, Pearson’s goodness-of-fit (χ2) statistic was utilized to evaluate the Hardy‒Weinberg equilibrium. The chi-square test or Fisher’s exact test was used to compare gene distributions among groups, and multivariate logistic regression was employed to calculate odds ratios (ORs) and 95% confidence intervals (CIs). The allelic associations between the preeclampsia group and the control group was evaluated. Taking advantage of the triad data, the transmission/disequilibrium test (TDT) statistic was calculated as follows: TDT χ2 = (a − b)2/(a + b). The p-value was calculated from a χ2 table with 1 df. Here, a was the number of times the interest allele is passed, and b was the number of times it has not been passed on. p < 0.05 was considered to indicate statistical significance.

3 Results

3.1 Clinical characteristics and SNP information

As shown in Table 1, no significant differences were observed between the two groups in age, BMI, previous misunderstandings, cesarean section and newborn sex (p > 0.05). However, there were statistical differences in gestational age, primiparity, parity, neonatal birth weight and neonatal length between the two groups (p < 0.05).

TABLE 1

VariablePreeclampsiaControlp
Mother
Age (years)34.19 ± 3.8633.31 ± 5.160.052
BMI (kg/m2)23.9 ± 3.724.2 ± 2.90.574
Gestational age (weeks)38.51 ± 2.5237.66 ± 1.920.019
Primiparity (%)58.1037.500.015
Parity (times)1.14 ± 0.271.97 ± 0.79<0.001
Previous miscarriages (%)29.7228.750.919
Caesarian section (%)39.1938.750.991
Offspring
weight (g)2,854.10 ± 711.923,271.62 ± 431.55<0.001
length (cm)48.04 ± 3.5749.27 ± 3.490.032
Male (%)59.4656.250.427

Clinical characteristics of preeclampsia group and control group.

In this study, we tested 3 SNPs (rs649216, rs1051456, and rs34785252) for the maternal KIR2DL4 gene and 4 SNPs (rs9380142, rs1063320, rs1130363, and rs1630185) for the HLA-G gene in all family triads. The genotype distribution of all 7 SNP loci in the control group was consistent with HWE (Supplementary Table S3).

3.2 Allele frequencies in the preeclampsia group and control group

Among the four SNP loci in the fetal HLA-G gene, the frequencies of rs9380142A (OR = 2.802, 95% CI = 1.761∼4.458, p = 0.001) and rs1063320G (OR = 1.807, 95% CI = 1.144∼2.852, p = 0.01) were greater in the preeclampsia group than in the control group. No significant differences in allele frequency for the other loci in the fetal HLA-G gene were observed between the two groups.

Among the four SNP loci in the paternal HLA-G gene, the frequencies of the rs1630185G (OR = 1.644, 95% CI = 1.005∼2.688, p = 0.046) and rs1130363A (OR = 1.644, 95% CI = 1.005∼2.688, p = 0.046) alleles were lower in the preeclampsia group than in the control group. Allele frequencies at other loci of the paternal HLA-G gene were not significantly different between the two groups.

Maternal allele frequencies for all SNP loci of the KIR2DL4 gene and HLA-G gene were not significantly different between the preeclampsia group and the control group. See Table 2.

TABLE 2

SNPsGeneAlleleMinor alleleMinor allele preeclampsiaMinor allele controlOR (95%CI)X2p
Mother
rs649216KIR2DL4C/TT23 (0.155)31 (0.193)0.765 (0.423∼1.385)0.7810.376
rs1051456KIR2DL4G/CG90 (0.608)80 (0.500)1.551 (0.986∼2.44)2.8540.056
rs34785252KIR2DL4C/AA62 (0.418)75 (0.468)0.817 (0.52∼1.282)0.7730.379
rs9380142HLA-GG/AA61 (0.412)80 (0.5)1.426 (0.908∼2.238)2.3890.122
rs1063320HLA-GC/GG53 (0.358)60 (0.375)0.929 (0.584∼1.478)0.0940.758
rs1630185HLA-GA/GG72 (0.486)66 (0.412)1.349 (0.859∼2.117)1.7010.192
rs1130363HLA-GG/AA73 (0.493)66 (0.412)1.386 (0.883∼2.175)2.0240.154
Father
rs9380142HLA-GG/AA45 (0.401)52 (0.333)1.389 (0.853∼2.262)1.3220.25
rs1063320HLA-GC/GG48 (0.428)72 (0.461)0.875 (0.536∼1.426)0.2860.592
rs1630185HLA-GA/GG56 (0.500)97 (0.621)1.644 (1.005∼2.688)3.9470.046
rs1130363HLA-GG/AA56 (0.500)97 (0.621)1.644 (1.005∼2.688)3.9470.046
Offspring
rs9380142HLA-GG/AA85 (0.574)62 (0.387)2.802 (1.761∼4.458)19.3520.001
rs1063320HLA-GC/GG74 (0.5)57 (0.356)1.807 (1.144∼2.852)6.4990.010
rs1630185HLA-GA/GG61 (0.412)75 (0.468)0.794 (0.506∼1.247)0.9980.317
rs1130363HLA-GG/AA62 (0.418)68 (0.425)0.975 (0.620∼1.533)0.0110.914

Distributions of allele frequencies in the preeclampsia group and the control group.

The meaning of bold is: p < 0.05, which is statistically significant.

3.3 Genotype frequencies in the preeclampsia group and control group

For fetal rs9380142 in the HLA-G gene, the alleles are A and G, and the A allele is the minor allele. The GG, GA and AA genotype frequencies of the preeclampsia group were 21.6%, 41.8%, 36.4%, and 0.00%, respectively, and those of the control group were 42.5%, 37.5%, and 20.0%, respectively (p = 0.01). The A allele, which has a lower frequency in the population, is assumed to be a mutant gene. In the dominant mode, mothers of fetuses carrying this mutation were more likely to develop preeclampsia (AA + GA vs. GG, OR = 2.196, 95% CI = 1.009–4.780; P = 0.047). Under the recessive mode, this difference was not significant (AA vs. GA + GG). Therefore, we suggest that if the fetal A allele is associated with the development of preeclampsia in the mother, it may function in the dominant mode.

For fetal rs1063320 in the HLA-G gene, the alleles are G and C; the G allele is the minor allele. The CC, CG and GG genotype frequencies were 25.6%, 48.6%, and 25.6%, respectively, in the preeclampsia group and 46.2%, 35.2%, and 17.5%, respectively, in the control group (p = 0.029). The G allele, which has a lower frequency in the population, is assumed to be a mutant gene. Under the dominant mode, mothers of fetuses carrying this mutation were more likely to develop preeclampsia (GG + GC vs. CC, OR = 2.491, 95% CI = 1.259–4.927; p = 0.009). Under the recessive mode, this difference was not significant (GG vs. GC + CC). Therefore, we suggest that if the fetal G allele is associated with the development of preeclampsia in the mother, it may function in the dominant mode.

The genotype frequencies of rs1630185 and rs1130363 in the paternal HLA-G gene were both significantly different between the preeclampsia group and the control group (p = 0.01 and p = 0.01). They were all inherited in a recessive model in the population. No significant differences in genotype frequencies for the remaining detected SNPs in fetuses, fathers or mothers were observed between the two groups. See Table 3.

TABLE 3

SNPsPreeclampsia (n = 74)Control (n = 80)pAB + BB vs. AABB vs. AA + AB
AAABBBTotalAAABBBTotalOR (95% CI)pOR (95% CI)p
Mathor
rs64921653 (0.716)19 (0.256)2 (0.027)7453 (0.662)23 (0.287)4 (0.050)800.7250.778 (0.392∼1.544)0.4730.528 (0.094∼2.970)0.468
rs105145613 (0.175)32 (0.432)29 (0.391)7423 (0.287)34 (0.425)23 (0.287)800.1960.626 (0.320∼1.227)0.1720.528 (0.245∼1.141)0.104
rs3478525225 (0.33736 (0.486)13 (0.175)7426 (0.325)33 (0.412)21 (0.262)800.4060.944 (0.482∼1.847)0.8660.599 (0.599∼1.305)0.197
rs938014231 (0.418)25 (0.337)18 (0.243)7424 (0.300)32 (0.400)24 (0.300)800.3040.594 (0.306∼1.156)0.1250.430 (0.367∼1.533)0.43
rs106332037 (0.5)21 (0.283)16 (0.216)7434 (0.425)32 (0.4)14 (0.175)800.3140.739 (0.391∼1.396)0.3511.300 (0.585∼2.892)0.519
rs163018518 (0.243)40 (0.54)16 (0.216)7427 (0.337)40 (0.5)13 (0.162)800.391.585 (0.783∼3.207)0.21.422 (0.631∼3.202)0.396
rs113036318 (0.243)39 (0.527)17 (0.229)7427 (0.337)40 (0.5)13 (0.162)800.3471.585 (0.783∼3.207)0.21.537 (0.688∼3.434)0.295
Father
rs938014225 (0.446)17 (0.303)14 (0.25)5637 (0.474)30 (0.384)11 (0.141)780.2541.119 (0.562∼2.229)0.7492.030 (0.843∼4.889)0.114
rs106332026 (0.464)12 (0.214)18 (0.321)5627 (0.346)30 (0.384)21 (0.269)780.1061.182 (0.571∼2.446)0.6521.464 (0.729∼2.942)0.284
rs163018518 (0.321)20 (0.357)18 (0.321)5628 (0.358)41 (0.525)9 (0.115)780.0101.182 (0.571∼2.446)0.6523.632 (1.487∼8.868)0.005
rs113036318 (0.321)20 (0.357)18 (0.321)5628 (0.358)41 (0.525)9 (0.115)780.0101.182 (0.571∼2.446)0.6523.632 (1.487∼8.868)0.005
Offspring
rs938014216 (0.216)31 (0.418)27 (0.364)7434 (0.425)30 (0.375)16 (0.200)800.0102.196 (1.009∼4.780)0.0471.633 (0.736∼3.622)0.227
rs106332019 (0.256)36 (0.486)19 (0.256)7437 (0.462)29 (0.362)14 (0.175)800.0292.491 (1.259∼4.927)0.0091.629 (0.748∼3.544)0.219
rs163018524 (0.324)39 (0.527)11 (0.148)7420 (0.250)45 (0.562)15 (0.187)800.5550.694 (0.344∼1.401)0.3090.757 (0.323∼1.773)0.521
rs113036322 (0.297)42 (0.567)10 (0.135)7427 (0.337)38 (0.475)15 (0.187)800.4771.204 (0.610∼2.378)0.5930.677 (0.283∼1.618)0.38

Distributions of genotype frequencies in the preeclampsia group and the control group.

A: major allele, B: minor allele. The meaning of bold is: p < 0.05, which is statistically significant.

3.4 Transmission of rs9380142A and rs1063320G in family triads

To further examine the preeclampsia-associated SNP loci in fetuses identified in this case‒control study and to exclude false-positive results due to factors such as population stratification, we analyzed the rs9380142 and rs1063320 loci by using the transmission disequilibrium test (TDT). It should be noted that TDT analysis is not usually performed in a control group, but since this study had a case‒control design, we believe that performing TDT analyses in the control triads, with the help of opposite trends of transmission, would strengthen the results found in the case triads. The offspring were not the patients in our study, but the mothers were. Based on these considerations, we performed TDT in both the preeclampsia and control triads.

The genotype distributions of the rs9380142 and rs1063320 loci within preeclampsia and control triads are shown in Supplementary Table S4. TDTs were performed when parent-of-origin was not considered or when paternal or maternal transmissions were considered separately. As shown in Table 4, while no significant results were observed in the preeclampsia triad, in the control triad, the rs9380142A allele was insufficiently transmitted (less frequently inherited from fathers and mothers to their offspring). Only 4 of the 18 heterozygous fathers transmitted the rs9380142A allele to their offspring (p = 0.018; X2 = 5.556, df = 1). Similarly, only 4 of 18 heterozygous mothers transmitted the rs9380142A allele (P = 0.018; X2 = 5.556, df = 1). TDT revealed that differences in the rs9380142G/A polymorphism in offspring between the preeclampsia group and the control group were associated with differential parental transmission (p = 0.001). The significantly insufficient transmission of the rs9380142A allele in the control group reinforces the results of the present study that rs9380142A is a preeclampsia-associated allele (the risk allele). There were no significant differences in the transmission of the rs1063320G allele from heterozygous fathers or mothers to their offspring in the preeclampsia and control triads.

TABLE 4

SNPAlleleΧ2, df 1p-value
rs9380142GA
Preeclampsia
Transmitted from heterozygous parent19102.7930.095
Transmitted from heterozygous mother1252.8820.09
Transmitted from heterozygous father530.50.48
Control
Transmitted from heterozygous parent28811.1110.001
Transmitted from heterozygous mother1445.5560.018
Transmitted from heterozygous father1445.5560.018
rs1063320CG
Preeclampsia
Transmitted from heterozygous parent14242.6320.105
Transmitted from heterozygous mother61220.157
Transmitted from heterozygous father61010.317
Control
Transmitted from heterozygous parent29280.0180.895
Transmitted from heterozygous mother1180.4740.491
Transmitted from heterozygous father14160.1330.715

Transmission disequilibrium test analysis of SNPs in HLA-G gene in family triads.

The meaning of bold is: p < 0.05, which is statistically significant.

3.5 Association of combined maternal KIR2DL4 and fetal HLA-G gene polymorphisms with preeclampsia

Maternal KIRs bind to fetal HLAs and regulate uNK cell function through the KIR-HLA receptor‒ligand pathway, which in turn affects placenta formation. Therefore, correlations between different maternal-fetal gene polymorphism combinations and preeclampsia were evaluated. The allele that has a lower frequency in the population is assumed to be a mutant allele. The entire SNP locus studied in the KIR2DL4 gene did not differ significantly between the two groups. The mothers were divided into mutant allele-containing and mutant allele-free groups. The TDT results indicated significant differences in allele transmission patterns between the preeclampsia group and the control group. Compared to that in preeclampsia triads, the transmission of the rs9380142A allele from parents to offspring is less frequent in the control triads. This finding suggested that the rs9380142A allele may increase the risk of preeclampsia by altering familial transmission dynamics, indicating that only rs9380142 in offspring is associated with the development of preeclampsia. rs9380142A in the HLA-G gene in offspring was inherited in a dominant mode in the study population. Therefore, the offspring were divided into two groups, i.e., minor allele-containing (AA + AG) and minor allele-free (GG) groups.

One-way ANOVA revealed that the genotype combinations of the three studied maternal KIR2DL4 gene SNP loci and fetal HLA-G gene rs9380142 were significantly different between the two groups (p < 0.05) (Table 5).

TABLE 5

MotherOffspringPreeclampsiaControlX2pMultiple logistic regression
OR95% CIp
rs649216rs9380142
CCAA/GA432711.3080.010Ref
CCGG10260.2420.101–0.5790.001
CT/TTAA/GA15190.4960.216–1.1380.098
CT/TTGG680.4710.147–1.5060.204
rs1051456rs9380142
CCAA/GA91511.8250.0080.380.148–0.9730.054
CCGG480.3160.088–1.1400.078
GC/GGAA/GA4931Ref
GC/GGGG12260.2920.120–0.6620.003
rs34785252rs9380142
CCAA/GA91612.8070.0050.3440.135–0.8770.025
CCGG480.3060.085–1.1050.071
CA/AAAA/GA4930Ref
CA/AAGG12260.2830.124–0.6420.003

The analysis of genotype combinations of maternal KIR2DL4 and fetal HLA-G.

The meaning of bold is: p < 0.05, which is statistically significant.

Multiple logistic regression analyses were performed using the maternal KIR2DL4-fetal HLA-G genotype combinations that differed between the two groups via one-way ANOVA as the independent variable, preeclampsia as the dependent variable, and the genotype combination with the highest incidence as the reference group (Table 5). The results showed that for maternal rs649216-fetal rs9380142, the CC-GG genotype combination had a significantly lower risk of morbidity than the reference (OR = 0.242, 95% CI = 0.101–0.579, p = 0.001), whereas the remaining two combinations were not significantly different from the reference. For maternal rs1051456-fetal rs9380142, the GC/GG-GG genotype combination had a significantly lower risk of morbidity than the reference (OR = 0.292, 95% CI = 0.120–0.662, p = 0.003), whereas the remaining two combinations were not significantly different from the reference. For maternal rs34785252-fetal rs9380142, the CC-AA/GA genotype combination and CA/AA-GG genotype combination both had a significantly lower risk of morbidity than the reference (OR = 0.344, 95% CI = 0.135–0.877, P = 0.025; OR = 0.283, 95% CI = 0.124–0.642, p = 0.003), whereas the risk of morbidity for the remaining combinations was not significantly different from that of the reference.

4 Discussion

KIR2DL4 is a special member of the killer cell immunoglobulin-like receptor (KIR) family that can inhibit or activate NK cells through its different N-terminal structures. It is structurally expressed in all NK cells (Colbern, et al., 1994; Goodridge, et al., 2007) and is also highly expressed in decidual NK cells and placental NK cells (Brusilovsky, et al., 2013). KIR2DL4 activates resting NK cells to produce cytokines without causing cytotoxicity (Rajagopalan, et al., 2001). The production of cytokines such as tumor necrosis factor alpha (TNF-α), interleukin-1β (IL-1β) and interferon-γ (IFN-γ) induces the production of vascular endothelial growth factor (VEGF) by trophoblast cells, which in turn affects uterine angiogenesis (Brusilovsky, et al., 2013). At present, there are few reports on the relationship between KIR2DL4 polymorphisms and preeclampsia (Tan, et al., 2009). Our study revealed no significant differences in the genotype frequencies or allele frequencies of the three SNP loci in the KIR2DL4 gene between mothers in the preeclampsia group and mothers in the control group. In case‒control studies in Malaysia and Australia, no association between maternal KIR2DL4 allele frequency and preeclampsia was observed (Witt, et al., 2002; Tan, et al., 2009). This is consistent with the results observed in our study. Therefore, we concluded that the maternal KIR2DL4 gene is not strongly associated with the development of preeclampsia if it is studied as an independent genetic factor.

During placental implantation, EVT damages the arterial wall as it approaches the uterine spiral arteries through the decidua, allowing the fetus to receive an adequate supply of blood from the mother (Pollheimer, et al., 2018). Unlike other cells, EVT cells express three HLA class I molecules, HLA-C, HLA-E, and HLA-G. A high level of HLA-G expression is unique to EVT cells (Papuchova, et al., 2020). Many scholars believe that HLA-G molecules play an important role in maternal-fetal immunotolerance. The expression of HLA-G in the placenta is associated with preeclampsia (Colbern, et al., 1994; Goldman-Wohl, et al., 2000; Luo, et al., 2018). It has been shown that HLA-G is absent or has reduced expression levels in the placenta of most patients with preeclampsia (Goldman-Wohl, et al., 2000). Our results showed that both the allele frequencies and genotype frequencies of rs9380142 and rs1063320 in fetuses were significantly different between the two groups. Both rs9380142 and rs1063320 are located in the 3′UTR of the gene and are associated with mRNA localization, mRNA stability and protein translation regulation. In the study fetal population, the rs9380142A and rs1063320G allele frequencies (MAFs) were significantly greater in the preeclampsia group than in the control group, and the ORs were greater than 1, suggesting that these loci may be risk factors for disease development in the fetus.

Further TDT analyses revealed no significant difference in the frequency of transmission of the rs1063320G allele from the heterozygote parent to the offspring between the preeclampsia group and the control group. We therefore believe that the association of this SNP with disease was likely to be a false-positive result due to factors such as population stratification and should be excluded from subsequent analysis. TDT also showed that there was a significant transmission deficit of the rs9380142A allele in the control group, with both heterozygous parents being much less likely to transmit A to their offspring. Since the fetal rs9380142A allele was a risk factor for preeclampsia, its undertransmission in controls was in line with expectations. The findings of the TDT reinforce the results of the present study that rs9380142A is a risk allele in fetuses. Therefore, we suggest that fetal HLA-G gene polymorphisms are associated with the development of preeclampsia and are risk factors for disease development. Mutations may disrupt maternal-fetal immune tolerance and ultimately increase the risk of disease development by decreasing mRNA stability and downregulating HLA-G protein expression (Michita, et al., 2018).

At the time of embryo implantation and placenta formation, maternal uNK cells synergize with fetal EVTs, allowing EVTs to invade the uterus to clear and replace the smooth muscle of maternal spiral arterioles and allowing the recasting of small uterine spiral arterioles to nourish the growing fetus (Wei, et al., 2022). Interactions between the uNK cell receptors and their ligands on EVTs must occur to strike a delicate balance between overinvasion, which damages the mother, and underinvasion, which is insufficient to supply blood and oxygen to the fetus (Pollheimer, et al., 2018). Hiby et al. (2004) reported that the maternal KIR AA genotype and fetal HLA-C2 allotype are more common in the preeclampsia group, indicating that when the fetus carries the HLA-C2 allotype, mothers who lack most or all activating KIRs have a significantly increased risk of developing preeclampsia. Another finding showed that the direct interaction between HLA-G (+) breast cancer cells and KIR2DL4 (+) tissue mast cells was positively correlated with lymph node metastasis or lymphatic invasion in clinical specimens (Riteau, et al., 2001). These findings suggest that maternal-fetal-specific gene combinations mediating the interaction of uNK cells with EVTs are strongly correlated with disease development. Our study revealed four KIR2DL4/HLA-G genotype combinations that were associated with the development of preeclampsia and significantly reduced the risk of the disease in mothers. Interestingly, in all four genotype combinations, the fetus carried the rs9380142G allele in the HLA-G gene. The role of fetal HLA-G in maternal-fetal tolerance has been well established (Lynge Nilsson, et al., 2014), and this topic has been reported in many studies (Loisel, et al., 2013; Hsu, et al., 2014). Combined with our findings, we believe that fetal HLA-G contributes more to the establishment of maternal-fetal tolerance and therefore may be more relevant to the pathological mechanisms of preeclampsia. However, we did not analyze HLA-G expression in placental tissue. This may be a limitation of this study. On the other hand, although we observed significant associations between certain maternal KIR2DL4 and fetal HLA-G polymorphisms and preeclampsia, these findings need to be validated in larger and more diverse populations. Additionally, we did not consider potential environmental and lifestyle factors that might influence the development of preeclampsia. Future studies should incorporate comprehensive genomic approaches and consider the interactions between genetic, environmental, and lifestyle factors.

In conclusion, our findings underscore significant HLA-G gene polymorphisms in fetuses within the Chinese Han population, implicating fetal HLA-G as pivotal in preeclampsia pathogenesis. Despite the lack of a direct association between maternal KIR gene polymorphisms and preeclampsia, specific maternal KIR2DL4 and fetal HLA-G combinations warrant further exploration. This study contributes crucial insights into maternal-fetal immune tolerance mechanisms, guiding future genetic and clinical studies on preeclampsia.

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 the Medical Ethics Committee of Kunming Medical University (No: kmmu2020MEC027). 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

YM: Methodology, Writing–original draft, Writing–review and editing. YQ: Data curation, Writing–original draft, Writing–review and editing. HJ: Formal Analysis, Software, Writing–original draft. HM: Data curation, Formal Analysis, Writing–original draft. YW: Data curation, Formal Analysis, Writing–original draft. YY: Writing–original draft, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. 1. Yunnan Applied Basic Research Projects-Joint Special Project , Grant/Award Number: 202401AY070001-193. 2. National Natural Science Foundation of China 81760273. 3. Health Research Project of Kunming Medical Technology Center, 2022-05-02-002. 4. Kunming Health Science and Technology Talent Training Project Medical Technology Center, 2022-SW (Tech)-16. 5. Leader Project of Medical Science and Technology Discipline in Kunming Health Science and Technology Talent Abstract, 2023-SW (Leader)-06.

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.

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

References

  • 1

    Al-JameilN.Aziz KhanF.Fareed KhanM.TabassumH. (2014). A brief overview of preeclampsia. J. Clin. Med. Res.6 (1), 1–7. 10.4021/jocmr1682w

  • 2

    BrusilovskyM.CordobaM.RosentalB.HershkovitzO.AndrakeM. D.PecherskayaA.et al (2013). Genome-wide siRNA screen reveals a new cellular partner of NK cell receptor KIR2DL4: heparan sulfate directly modulates KIR2DL4-mediated responses. J. Immunol.191 (10), 5256–5267. 10.4049/jimmunol.1302079

  • 3

    ColbernG. T.ChiangM. H.MainE. K. (1994). Expression of the nonclassic histocompatibility antigen HLA-G by preeclamptic placenta. Am. J. Obstet. Gynecol.170 (5 Pt 1), 1244–1250. 10.1016/s0002-9378(94)70134-2

  • 4

    ColucciF. (2017). The role of KIR and HLA interactions in pregnancy complications. Immunogenetics69 (8-9), 557–565. 10.1007/s00251-017-1003-9

  • 5

    EmmerP. M.JoostenI.SchutM. H.ZusterzeelP. L.HendriksJ. C.SteegersE. A. (2004). Shift in expression of HLA-G mRNA spliceforms in pregnancies complicated by preeclampsia. J. Soc. Gynecol. Investig.11 (4), 220–226. 10.1016/j.jsgi.2003.10.011

  • 6

    FerreiraL. M. R.MeissnerT. B.TilburgsT.StromingerJ. L. (2017). HLA-G: at the interface of maternal-fetal tolerance. Trends Immunol.38 (4), 272–286. 10.1016/j.it.2017.01.009

  • 7

    Gestational Hypertension and Preeclampsia (2020). Gestational Hypertension and Preeclampsia: ACOG Practice Bulletin, Number 222. Obstet. Gynecol.135 (6), e237–e260. 10.1097/aog.0000000000003891

  • 8

    Goldman-WohlD. S.ArielI.GreenfieldC.Hochner-CelnikierD.CrossJ.FisherS.et al (2000). Lack of human leukocyte antigen-G expression in extravillous trophoblasts is associated with pre-eclampsia. Mol. Hum. Reprod.6 (1), 88–95. 10.1093/molehr/6.1.88

  • 9

    GoodridgeJ. P.LathburyL. J.SteinerN. K.ShulseC. N.PullikotilP.SeidahN. G.et al (2007). Three common alleles of KIR2DL4 (CD158d) encode constitutively expressed, inducible and secreted receptors in NK cells. Eur. J. Immunol.37 (1), 199–211. 10.1002/eji.200636316

  • 10

    HibyS. E.WalkerJ. J.O'Shaughnessy KM.RedmanC. W.CarringtonM.TrowsdaleJ.et al (2004). Combinations of maternal KIR and fetal HLA-C genes influence the risk of preeclampsia and reproductive success. J. Exp. Med.200 (8), 957–965. 10.1084/jem.20041214

  • 11

    HsuP.Santner-NananB.JoungS.PeekM. J.NananR. (2014). Expansion of CD4(+) HLA-G(+) T Cell in human pregnancy is impaired in pre-eclampsia. Am. J. Reprod. Immunol.71 (3), 217–228. 10.1111/aji.12195

  • 12

    KaminskiV. L.EllwangerJ. H.ChiesJ. A. B. (2019). Extracellular vesicles in host-pathogen interactions and immune regulation - exosomes as emerging actors in the immunological theater of pregnancy. Heliyon5 (8), e02355. 10.1016/j.heliyon.2019.e02355

  • 13

    KieckbuschJ.GaynorL. M.MoffettA.ColucciF. (2014). MHC-dependent inhibition of uterine NK cells impedes fetal growth and decidual vascular remodelling. Nat. Commun.5, 3359. 10.1038/ncomms4359

  • 14

    LoiselD. A.BillstrandC.MurrayK.PattersonK.ChaiworapongsaT.RomeroR.et al (2013). The maternal HLA-G 1597ΔC null mutation is associated with increased risk of pre-eclampsia and reduced HLA-G expression during pregnancy in African-American women. Mol. Hum. Reprod.19 (3), 144–152. 10.1093/molehr/gas041

  • 15

    LuoS.PeiJ.LiX.GuW. (2018). Decreased expression of JHDMID in placenta is associated with preeclampsia through HLA-G. J. Hum. Hypertens.32 (6), 448–454. 10.1038/s41371-018-0062-1

  • 16

    Lynge NilssonL.DjurisicS.HviidT. V. (2014). Controlling the immunological crosstalk during conception and pregnancy: HLA-G in reproduction. Front. Immunol.5, 198. 10.3389/fimmu.2014.00198

  • 17

    MaoJ.FengY.ZhuX.MaF. (2023). The molecular mechanisms of HLA-G regulatory function on immune cells during early pregnancy. Biomolecules13 (8), 1213. 10.3390/biom13081213

  • 18

    MichitaR. T.KaminskiV. L.ChiesJ. A. B. (2018). Genetic variants in preeclampsia: lessons from studies in Latin-American populations. Front. Physiol.9, 1771. 10.3389/fphys.2018.01771

  • 19

    MoffettA.ShreeveN. (2015). First do no harm: uterine natural killer (NK) cells in assisted reproduction. Hum. Reprod.30 (7), 1519–1525. 10.1093/humrep/dev098

  • 20

    MukherjeeI.SinghS.KarmakarA.KashyapN.MridhaA. R.SharmaJ. B.et al (2023). New immune horizons in therapeutics and diagnostic approaches to Preeclampsia. Am. J. Reprod. Immunol.89 (2), e13670. 10.1111/aji.13670

  • 21

    NakimuliA.ChazaraO.HibyS. E.FarrellL.TukwasibweS.JayaramanJ.et al (2015). A Kir B centromeric region present in Africans but not Europeans protects pregnant women from pre-eclampsia. Proc. Natl. Acad. Sci. U. S. A.112 (3), 845–850. 10.1073/pnas.1413453112

  • 22

    PapuchovaH.KshirsagarS.XuL.Bougleux GomesH. A.LiQ.IyerV.et al (2020). Three types of HLA-G+ extravillous trophoblasts that have distinct immune regulatory properties. Proc. Natl. Acad. Sci. U. S. A.117 (27), 15772–15777. 10.1073/pnas.2000484117

  • 23

    ParhamP. (2004). NK cells and trophoblasts: partners in pregnancy. J. Exp. Med.200 (8), 951–955. 10.1084/jem.20041783

  • 24

    PollheimerJ.VondraS.BaltayevaJ.BeristainA. G.KnöflerM. (2018). Regulation of placental extravillous trophoblasts by the maternal uterine environment. Front. Immunol.9, 2597. 10.3389/fimmu.2018.02597

  • 25

    RajagopalanS.FuJ.LongE. O. (2001). Cutting edge: induction of IFN-gamma production but not cytotoxicity by the killer cell Ig-like receptor KIR2DL4 (CD158d) in resting NK cells. J. Immunol.167 (4), 1877–1881. 10.4049/jimmunol.167.4.1877

  • 26

    RechesA.BerhaniO.MandelboimO. (2020). A unique regulation region in the 3' UTR of HLA-G with a promising potential. Int. J. Mol. Sci.21 (3), 900. 10.3390/ijms21030900

  • 27

    RitariJ.HyvärinenK.PartanenJ.KoskelaS. (2022). KIR gene content imputation from single-nucleotide polymorphisms in the Finnish population. PeerJ10, e12692. 10.7717/peerj.12692

  • 28

    RiteauB.MenierC.Khalil-DaherI.MartinozziS.PlaM.DaussetJ.et al (2001). HLA-G1 co-expression boosts the HLA class I-mediated NK lysis inhibition. Int. Immunol.13 (2), 193–201. 10.1093/intimm/13.2.193

  • 29

    TanC. Y.ChongY. S.LoganathA.ChanY. H.RavichandranJ.LeeC. G.et al (2009). Possible gene-gene interaction of KIR2DL4 with its cognate ligand HLA-G in modulating risk for preeclampsia. Reprod. Sci.16 (12), 1135–1143. 10.1177/1933719109342280

  • 30

    TangY.LiuH.LiH.PengT.GuW.LiX. (2015). Hypermethylation of the HLA-G promoter is associated with preeclampsia. Mol. Hum. Reprod.21 (9), 736–744. 10.1093/molehr/gav037

  • 31

    WeiX. W.ZhangY. C.WuF.TianF. J.LinY. (2022). The role of extravillous trophoblasts and uterine NK cells in vascular remodeling during pregnancy. Front. Immunol.13, 951482. 10.3389/fimmu.2022.951482

  • 32

    WittC. S.WhitewayJ. M.WarrenH. S.BardenA.RogersM.MartinA.et al (2002). Alleles of the KIR2DL4 receptor and their lack of association with pre-eclampsia. Eur. J. Immunol.32 (1), 18–29. 10.1002/1521-4141(200201)32:1<18::AID-IMMU18>3.0.CO;2-7

Summary

Keywords

preeclampsia, polymorphism, Placenta, HLA-G, KIR2DL4

Citation

Ma Y, Qian Y, Jiang H, Meng H, Wang Y and Yang Y (2024) Combined maternal KIR2DL4 and fetal HLA-G polymorphisms were associated with preeclampsia in a Han Chinese population. Front. Genet. 15:1442938. doi: 10.3389/fgene.2024.1442938

Received

04 June 2024

Accepted

15 July 2024

Published

31 July 2024

Volume

15 - 2024

Edited by

Ali H. Hajeer, King Abdulaziz Medical City, Saudi Arabia

Reviewed by

Nagarjuna Sivaraj, Great Eastern Medical School and Hospital, India

Valéria de Lima Kaminski, Anhembi Morumbi University, Brazil

Updates

Copyright

*Correspondence: Yuling Yang,

‡ 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

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics