1. Introduction
Human Leukocyte Antigen-G (HLA-G) is a member of non-classical human major histocompatibility complex class Ib family human major histocompatibility complex exhibiting a limited polymorphism. HLA-G is a key molecule in maternal-fetal immune tolerance maintaining immune homeostasis [
1]. It has been described initially in extravillous trophoblast (EVT), but later has been identified cornea, erythroid precursors, pancreatic islets, endothelial cell precursors and erythroblasts and thymus, and displays different characteristics that of classical HLA-A, HLA-B and HLA-C (class Ia) molecules.
The HLA-G molecule constitutes of a heavy chain subunit comprised of three extracellular domains, located on chromosome 6, associated non-covalently with β2-microglobulin, encoded on chromosome 15. The HLA-G gene has 8 exons and 7 non-coding intervening sequences. Depending on the localization, membrane- bound and soluble HLA-G molecules are known due to seven different isoforms resulting from alternative splicing and differential association with 𝛽2-microglobulin. The membrane-bound subtypes are HLA-G1-4, which are expressed by ET cells in the placenta. The soluble forms released from the cell (HLA-G1, -G5, -G6, and -G7), due to the lack of the transmembrane and intracellular domains of membrane-bound HLA-G. The soluble subtypes are produced by various cells such as extravillous and villous cytotrophoblasts, chorion and amnion membranes, syncytiotrophoblasts, decidual stromal cells, fetal endothelial cells, and immune cells in the placenta [
2,
3].
Basically, HLA-G exhibits an immune tolerant function by direct inhibition of different immune-competent cells. The suppressing potential is mediated by HLA-G binding to its target receptors. The HLA-G receptors are inhibitory receptors such as killer cell immunoglobulin-like receptor 2DL4 (KIR2DL4) and immunoglobulin-like transcription receptor type 2 and type 4 (ILT2 and ILT4 respectively) [
4,
5].
Killer-cell immunoglobulin-like receptor (KIR) belongs to the family of type I transmembrane glycoproteins, which are mainly expressed on the plasma membrane of two subtypes of the natural killer (NK) cells: decidual natural killer (dNK) and peripherial blood natural killer (pbNK) cells. The ligands for several KIR proteins are subsets of HLA class I molecules, therefore KIR proteins are thought to play a significant role in down-regulation of maternal immune response. KIR2DL4 is a member of unusual killer cell immunoglobulin-like receptor (KIR) family in terms of its structure, expression, cellular localization and signaling properties. Unlike all other KIRs, which are expressed on the surface of NK cells, KIR2DL4 resides in endosomes. HLA-G, particularly the sHLA-G isoforms are solely ligands of the KIR2DL4 [
4,
5]. Importantly, HLA-G inhibits NK cell cytotoxicity and required to protect trophoblasts against NK cell-induced lysis [
6] yet the mechanism by which EVTs uniquely express HLA-G remains unknown. Since NK cells constitute the predominant lymphocyte subset in the placenta, the pro-inflammatory/pro-angiogenic outcome of the interaction between KIR2DL4 and soluble HLA-G supports the role of KIR2DL4 in extensive remodeling of maternal vasculature during early weeks of pregnancy [
7]. The ligand-receptor interaction activates the angiogenesis and vasculogenesis-related processes and inhibits the cytotoxic activity of the NK cells. sHLA-G isoforms induce trophoblast invasion and remodeling of spiral arteries by producing pro-inflammatory cytokines (IL-6, TNF-alfa) and proangiogenic factors (VEGF-A, MMP) by the NK cells [
1,
8]. sHLA-G promotes the secretion of IL-8 in NK cells and machrophages manifesting in trophoblast invasion. The activation processes encompass production of interferon- γ (IFN- γ) by the non-activated NK cells enhancing the paracrine secretion of HLA-G. The HLA-G induced NK-derived IFN-γ is essential for the decidual disorganization, vessel modification through promoting cell adhesion, smooth muscle proliferation [
9].
Moreover, HLA-G interacts with ILT2 and ILT4 receptors, thereby inhibits cytotoxicity and enhances the production of anti-inflammatory cytokines (IL-4, IL-10, TGF-beta). ILT2 is on the surface of T cells, B cells, natural killer (NK) cells and antigen presenting cells (APCs), whereas ILT4 is exclusive to dendritic cells, macrophages/monocytes and APCs. HLA-G molecule binds with high affinity to both ILT2 and ILT4 inhibiting the maternal immune response [
10]. The direct interaction of HLA-G with different immune cell subpopulations and endothel cells is accompanied by the generation and maintenance of tolerance at different stages of the immune response, e.g., differentiation, proliferation, cytolysis and cytokine secretion. On the other hand, HLA-G can be bound to ILT2 on dNK-cells and activates the PI3K-ACT signal transmission path. This interaction stimulates the production of growth promoting factors, such as pleiotrophin, osteoglycin and osteopontin, which are crucial for embryonic and fetal development [
1,
8]. Soluble forms have an antiproliferative effect on CD8+ T lymphocytes through promoting their apoptosis. Furthermore, they trigger the Th1-Th2 and Th17-regulatory T cells shift to reduce the maternal immune reaction [
1,
2].
Its major immunosuppression effected is exerted via bounding to its receptors appearing on the surface of immunocompetent cells (CD14+ monocytes NK, B, T cells and APCs) in any other organs as well apart from the placenta. HLA-G has been implicated in immune modulatory processes in organ transplantation, viral infections, cancer progression, autoimmunity, and even embryonic tissues. Both soluble and membrane-bound HLA-G have immunosuppressive effect and may contribute to tissue remodelling in both the maternal decidua and in tumors as it is protective against a graft-versus-host reaction since play a role in avoiding immune surveillance [
11,
12]. Both in cancers and in throphoblasts, the HLA-G may be important in tissue remodelling [
1,
11].
Under pathophysiological conditions HLA-G antigens are expressed on various types of malignant cells signifying that HLA-G antigen expression on tumors is a potent strategy to avoid immune surveillance. Several polymorphisms in the 3’UTR region condition changes in HLA-G expression (14bp and +3142C/G, among others), which have been associated with both the development and outcome of patients with different tumor types [
11]. Since HLA-G antigens can be expressed on immune cells (including macrophages, neutrophils and NK cells) that infiltrates the tumor mass by chemoattraction to pro-inflammatory cytokines secreted by the cancer cells. The peripheral blood from cancer patients includes also HLA-G presenting cells [
13], and HLA-G is suggested to impair the immune response to the tumor. HLA-G+ APC can inhibit the function of CD4+ T cells and induce their differentiation into regulatory T cells [
14]. These immune cells contribute to tissue remodelling within cancer, including promotion of metastasis of tumor cells and angiogenesis, furthermore NK cells exhibit citotoxicity towards tumor cells [
15]. The depletion of NK cells increases tumor growth and low pbNK cytotoxic activity is associated with increased risk of cancer [
16], whereas NK cell presence in tumor is positively associated with prognosis [
17]. Recent research describes the potential use of membrane-bound and soluble HLA-G as a prognostic biomarker to identify tumors and to assess disease stage, as well as on the use of HLA-G as a novel therapeutic target in cancer.
Some studies have found a reduced expression and polymorphism of HLA-G in pregnancy pathologies such as preeclampsia [
18] or recurrent spontaneous abortion [
19] in comparison with healthy placentas.
Soluble HLA-G can be detected both in maternal blood [
7,
8,
9] and amniotic fluid [
20]. In this study, we have investigated sHLA-G (sHLA-G1 and sHLA-G5) concentration in amniotic fluid and in maternal serum in mid-pregnancy. Considering uncomplicated pregnancies, we have investigated the correlations of sHLA-G levels and fetoplacental growth during mid-pregnancy.
2. Materials and Methods
2.1. Study Design
A prospective, cross-sectional cohort study was conducted in pregnant women undergoing amniocentesis at the Department of Obstetrics and Gynecology, University of Szeged, Hungary between January 2021 and May 2021. The demography of the study group was homogenic (Caucasian female subjects). During the study period, all singleton pregnancies with increased risk of chromosomal abnormality, where amniocentesis (AC) was performed between 16+0 and 22+0 weeks of gestation, were recruited into our study. The indications for AC were increased nuchal translucency (NT) at first trimester scan (≥2 MoM for gestational age (GA)) (n=2), chromosome aberration or gene disorder concerning the previous pregnancy (n=2), and advanced maternal age (n=37).
Exclusion criteria of the study were identified as follows: multiple pregnancies; fetal or neonatal structural or genetic anomaly, improper localization of the placenta for sonography (placenta praevia, posterior placenta), pathological placentation (placenta accreta spectrum), self-reported drug, alcohol, caffeine or nicotine abuse or exposure to circulatory medication (oxerutins, calcium dobesilate) and systemic disease (e.g., essential hypertension, any type of pregestational diabetes mellitus, autoimmune disease, vasculitis, hemophylia, thrombophylia, chronic infections).
In addition, women with complications during late pregnancy (gestational diabetes mellitus treated with diet, hypertension-related diseases, small for gestational age at delivery, large for gestational age at delivery) were excluded from the study. Therefore 41 healthy pregnancies were included in our study.
The study protocol was approved by the Clinical Research Ethics Committee of the University of Szeged (date of approval: 10 February 2017, reference number: SZTE 09/2017). The study was carried out according to the Principles of the Declaration of Helsinki. We obtained written informed consent from all participants.
2.2. Conventional 2-Dimensional (2-D) Sonographic Examinations
All pregnancies were dated by using the measurement of crown rump length (CRL) at nuchal screening. NT and anatomic assessment between 11+0 and 13+6 weeks were performed by utilizing conventional methods. Ultrasound examination took place before measuring AC to determine the number of fetuses, fetal biometry, fetal anomalies, placental location and the amount of amniotic fluid. Fetal weight was estimated according to the method of Hadlock et al.[
21]. after measuring the necessary sonographic parameters (biparietal diameter, head circumference, abdominal circumference and femur length). Estimated fetal weight percentile was calculated according to the local standards [
22]. The ultrasound investigations were conducted by J. S. and A. S.
2.3. Volume Acquisition
The images used for the determination of placental volume and 3-dimensional Power Doppler (3-DPD) indices were acquired at the time of the visit. All 3-D scans were performed by A. S. Voluson 730 Expert ultrasound machine (GE Medical Systems, Kretztechnik GmbH & Co OHG, Tiefenbach, Austria) equipped with a multifrequency probe (2–5 MHz) was used to acquire all images. Each sample was examined using 3D rendering mode, in which the color and gray value information was processed and combined to give a 3D image (mode cent; smooth: 4/5; FRQ: low; quality: 16; density: 6; enhance: 16; balance: 150; filter: 2; actual power: 2 dB; pulse repetition frequency: 0.9). We used fast low-resolution acquisition to avoid any kind of artifacts. The 3-D static volume box was placed over the highest villous vascular density zone at the umbilical cord insertion [
23]. Each image was recovered from the disc in succession for processing. We recorded one sample from each patient during gestation.
2.4. Determination of Power Doppler Indices
The stored volumes were further analyzed using the virtual organ computer-aided analysis (VOCAL) program pertaining to computer software 4DView (GE Medical Systems, Zipf, Austria, version 10.4) by the same expert in 3-D analysis (A. S.). The image used for recovering from the hard disc was captured and processed using multiplanar system. The spherical sample volume was consistently 28 mL. The VOCAL program calculated automatically the grey- and color scale values from the acquired spherical sample volume in a histogram in all cases. The combined use of power Doppler with three-dimensional ultrasound provides the possibility of quantifying blood in motion within a volume of interest. Three indices were calculated, namely the vascularization index (VI), flow index (FI), and vascularization flow index (VFI), as estimates of the percentage of the volume filled with detectably moving blood. VI (expressed as a percentage) is the proportion of color voxels in the studied volume, representing the proportion of blood vessels within the tissue. FI (expressed as a scale of 0–100) is the average value of all color voxels, representing the average power Doppler amplitude within blood vessels. VFI (expressed as a scale of 0–100) is the average color value of all grey and color voxels, a product of the number of color voxels as a percentage and the relative amplitude of these voxels [
24,
25].
The intra-observer errors were evaluated by repeated measurements of 3-DPD indices at initiation of the study. The intra-class correlation coefficients for all Doppler indices were excellent (0.99) in case of all indices.
2.5. Amniocentesis Procedure
The patients were informed about the procedure and possible complications before a consent form was signed prior to the procedure. All procedures were performed by the same operator expert (J.S.) at the outpatient unit, who followed the standard protocol. A local antiseptic was applied to the skin. A 22-gauge spinal needle was inserted under continuous ultrasound guidance, and needle insertion through the placenta was avoided. Amniotic fluid (8–10 mL) was taken, and the first 2 mL of each sample was discarded to prevent contamination with maternal cells. Blood-contaminated amniotic fluid was not utilized. Fetal heart rate was evaluated after the procedure, and no stillbirth or premature rupture was observed. Following amniocentesis, anti-D immunoglobulin was administered, when it was necessary.
2.6. Samples
Amniotic fluid and maternal venous blood were collected from each patient at the time of amniocentesis. Blood samples were centrifuged at 3400 rpm for 15 minutes. Serum and amniotic fluid samples were stored at −80 °C until assay.
2.7. Enzyme-Linked Immunosorbent Assay (ELISA)
Human sHLA-G in maternal serum and amniotic fluid were determined by ELISA. The laboratory staff members who performed the assays were blinded to pregnancy outcome, and the clinician recruiting women did not participate at analyzing the samples.
Concentration of sHLA-G was measured using the kits from Elabscience Biotechnology Corporation (Houston, TX, USA). The sensitivity of assay was 0.38 ng/ml. The intra- and inter-assay coefficients of variation were <10% according to the manufacturer.
2.8. Data and Statistical Analysis
Statistical analyses were performed using SPSS version 23 (IBM Corp., Armonk, NY, USA). Continuous variables were expressed as mean ± standard deviation (SD) and categorical variables were expressed as numbers and percentages. The relationship between the level of sHLA-G and other continuous variables was assessed using univariate and multivariate regression analyses characterized by correlation coefficient (ß) and 95% confidence interval (CI). Multiple linear regression was adjusted for well-known confounders such as maternal age, body mass index (BMI) at the time of amniocentesis, number of previous pregnancies and gestational age at the time of amniocentesis, as these factors determine the actual placental volume and fetal weight. Paired samples t-test was applied to analyze the differences between serum and amniotic levels of the analytes. The two-tailed statistical significance level was set at 5% and p-values were adjusted using Holm–Bonferroni correction for multiple comparisons.
3. Results
Descriptive statistics are shown in
Table 1. As expected by guidelines, amniocentesis has been proposed for pregnant women with advanced maternal age (above 37 years) and suspicion of aneuploidy. Mean maternal age (33.63 years) in our case group was minimally higher compared to the national average reference age at delivery, which was 29.1 years for primiparous women and 30.4 years for the total number of pregnant women in 2021 [
26]. Women who participated in our study had a mean BMI of 26.4 kg/m
2 and slightly overweighted and 29.3% of the participants were primiparous. Birth weights of our sample was also near to that of the national average of healthy pregnant. As expected, the birth weight of the offsprings very without extremities. Ultrasound characteristics are demonstrated in
Table 2. Large NT was scarcely observed in our group (n=2). Mean gestational age at the time of amniocentesis was 18.4 weeks and mean gestational age at the time of delivery was 39.0 weeks (range: 37.0-41.0 weeks). Typically, all of the fetal sonographic measurements were with no outliers and with a low standard deviation (<1) explicating that the data are clustered tightly around the mean.
Table 3 gives an overview of the levels of angiogenic factors. The mean concentration was 53.39 ng/ml for sHLA-G in amniotic fluid and 51.05 ng/ml for sHLA-G in serum.
Figure 1 presents the levels of angiogenic factors in body fluids adjusted to the gestational age. The sHLA-G levels were steady during the gestational period, which was investigated in our study
. The sHLA-G level in serum was significantly associated with sHLA-G level in amniotic fluid, both in the unadjusted (β
=0.63, 95% CI= 0.30-0.72, p<0.01) and in the adjusted model (β
=0.66, 95% CI= 0.32-0.74, p<0.01) controlled for age, BMI, previous parity and gestational age at amniocentesis.
Table 4 displays the sonographic correlates of sHLA-G levels. The sHLA-G serum concentration was significantly correlated positively with abdominal circumference percentiles in the unadjusted analysis
(β
=0.41, 95% CI= -0.08-0.75, p<0.05).
Similarly, we found a significant interaction between sHLA-G levels in serum and the estimated fetal weight percentile (β
=0.41, 95% CI= -0.02-0.84, p<0.05). An inverse correlation could be detected between sHLA-G level in serum and FI in the unadjusted linear regression model
(β
=-0.34, 95% CI= -3.58-0.46, p<0.05). Reduced nuchal translucency appears to be correlated with high sHLA-G level in amniotic fluid (β
=-0.30, 95% CI= -21,73-0.04, p<0.05). Indirect correlations were observed between sHLA-G concentrations in amniotic fluid and placental perfusion [VI:
univariate analysis: β
=-0.34, 95% CI= -2.13-0.06, p<0.05, multivariate analysis: β
=-0.38, 95% CI= -2.47-0.03, p<0.05; and VFI: univariate analysis: β
=-0.44, 95% CI= -3.28-0.63, p<0.05, multivariate analysis: β
=-0.51, 95% CI= -3.79-0.72, p<0.05].
4. Discussion
HLA-G is a molecule that was first known to confer protection to the fetus from destruction by maternal immune system, thus critically contributing to maternal-fetal immune tolerance. The appropriate vasculogenesis and angiogenesis have an important role in placental development. The trophoblast invasion and remodeling of spiral arteries are regulated by the balance and interaction between pro-angiogenic and anti-angiogenic factors. The imbalance of angiogenic factors can lead to gestational pathological conditions, such as gestational hypertension, preeclampsia and prematurity [
8,
20,
27,
28,
29]
We have found a positive correlation between sHLA-G levels in serum and amniotic fluid, even when adjusted for confounders. sHLA-G in maternal circulation is secreted by EVT cells, while sHLA-G expressed by amniotic membrane is detectable in amniotic fluid, ensuring that it is involved in both local and systemic maternal immune modulation, which is in line with previous reports [
19,
20]. The association can be explained that the local and systemic immune response is in concordance in its extent. However, why the amniotic membrane cells secret the sHLA-G into the amniotic fluid is in line with the sHLA-G serum level is not understood, but maybe it is related to the size of the fetus. A larger of amniotic epithel cells area may present a higher cell mass that is interacting with a larger amount of immune cells.
Furthermore, HLA-G is primarily a checkpoint molecule. In our research, we have realized that sHLA-G levels in maternal serum and amniotic fluid are constant during mid-pregnancy, although maternal serum level has been decreased slightly and non-significantly.
Maternal serum level exhibits a positive correlation with fetal growth in terms of abdominal circumference and estimated fetal weight, which is expressed in percentiles after sonographic measurement. HLA-G levels in serum have a wide range of variety in a relatively short time period during gestation as published by others [
30]. Furthermore, fetal growth percentiles in our sample are constant and independent of gestation between 16th and 22nd weeks. However, our study supports the idea that sHLA-G in serum is able to predict the weight percentile among fetuses with normal weight during this short time period of the study. Moreover, a previous study has revealed that sHLA-G is significantly decreased in maternal blood of patients having the diagnosis of fetal growth restriction [
31]. The possible molecular explanation can be that sHLA-G stimulates decidual NK cells that promotes fetal growth and controls trophoblast invasion during mid-pregnancy [
1,
32]. Another explanation could be that HLA-G has been found to facilitate fetal growth by stimulating secretion of growth promoting factors (GPFs) in NK cells, according to recent studies [
32,
33,
34].
One can speculate that a higher EVT mass can coordinate a larger amount of spiral arteries that can supply the higher demand of gas and nutrition exchange of a larger fetus.
In contrast with our former article where the dataset on normal and pathologic pregnancies altogether was interpreted [
35], our current results show no association between serum sHLA-G levels and placental volume. An explanation could be that our new sample was included only asymptomatic, healthy pregnancies with no extreme values which should correspond to pathologic pregnancies.
HLA-G is unusual among HLA molecules in its unique pattern of expression in healthy individuals [
36]. Similar to our former results among healthy and complicated pregnancies [
28], sHLA-G level in amniotic fluid exhibited a negative correlation with the vascular perfusion of placenta. In this study we have explored an inverse correlation between HLA-G level in maternal serum and the placental blood perfusion with regard to healthy pregnancies. We assume that sHLA-G secretion may be reactive to reduced blood flow and vascular network in the expanding placenta, since sHLA-G facilitates the process of vasculogenesis [
30]. The expansion of placental volume may outweigh the increase in capillary branching. In order to evaluate the placental perfusion, we used “Mercé-type sonobiopsy”, which is a reproducible and validated method [
37,
38], and by obtaining a representative sample of the placental tree, it is applicable throughout the entire pregnancy, unlike other methods [
24], when the entire placenta needs to be visualized [
39,
40]. It can be also hypothesized that the HLA-G specific interaction of decidual natural killer cells and macrophages provides the secretion of angiogenic and pro-inflammatory factors implicated in vascular remodeling and the extent of trophoblast invasion, which is connected to lower perfusion indices as previously detected by ultrasound examinations.
Interestingly, we have unraveled an indirect correlation between NT measurements during 11-14 weeks of gestation and sHLA-G levels in amniotic fluid during 16 and 22 weeks of gestation. Our results support the assumption that low vascularization and perfusion of the placenta with high sHLA-G level in amniotic fluid during second trimester is preceded by an impaired placental 3D-PD indices and higher NT in the first trimester. This is in line with the results of Metzenbauer et al.[
41], who found that trisomies are associated with higher NT and smaller placental volume between 11 to 14 weeks of gestation. Based on the literature we can confirm that low sHLA-G expression is commonly associated with implantation failure and miscarriage [
42,
43], which is supported by our results as well, since an obvious correlation exists between increased nuchal translucency and fetal loss [
44].
Our study is focused on a relatively large number of amniotic samples deriving from a single center, where the ultrasound measurements were performed by two highly skilled sonographers.
In conclusion, increasing amount of data are collected concerning the impact of HLA-G polymorphism on the quantity of HLA-G expression and the correlation with complications in pregnancy. sHLA-G is unchanged in serum and amniotic fluid during second trimester in high-risk, but non-complicated pregnancies. Furthermore, our results indicate the role of decreased amniotic concentration of sHLA-G in the pathogenesis of increased nuchal translucency. It is of clinical importance that serum sHLA-G might be utilized as a predictor of anomaly of fetal growth and sHLA-G in amniotic fluid and maternal serum can be a potent biomarker of placental perfusion. To establish sHLA-G as a meaningful clinical biomarker, it is of paramount importance to trace back the structural diversity to tolerance-mediating functions.