Highlighting allelic variations at the interleukin-19 locus in term of preeclampsia predisposing factors and access to an accurate diagnostic/screening option

Abstract

Background

Preeclampsia is the main cause of preterm parturition and maternal–fetal complications. T helper 1 and T helper 2 cytokines balance is a requirement in normal pregnancy and aberrant in this immunologic balance, play an important role in the pathology of preeclampsia. In previous studies single nucleotide polymorphisms have been associated with the alteration of serum cytokine levels.

Objective

This study was aimed to discover association between interleukin-13 (rs20541, and rs56035208) and interleukin-19 (rs1028181 (T/C) and rs2243191(T/C)) polymorphisms with susceptibility to preeclampsia.

Methods

In this case-control study 300 women with and without preeclampsia (n = 150/each) who referred to Zeynabieh Hospital- Shiraz, Iran, from February 2021 to April 2022 were enrolled. For genotyping the interleukin-13 and interleukin-19 polymorphisms, the Allele-specific polymerase chain reaction and direct sequencing method was carried out.

Results

Our statistical results revealed no significant differences in allele and genotype frequencies for interleukin-13 polymorphisms compared to controls. We found that the interleukin‐13 polymorphisms are significantly associated with vulnerability to edema at rs20541 position and maternal drinking at rs56035208 position. But it was interesting to note that the differences of both the allele and genotype frequencies of interleukin-19 polymorphisms and their contribution to the risk of preeclampsia susceptibility were significant.

Conclusions

No risk of preeclampsia was found in all comparisons for interleukin-13 polymorphisms. However, the interleukin-19 polymorphisms were found to confer the risk of preeclampsia in our population.

Introduction

Preeclampsia is a common pregnancy-related disorder that affects 5–8% of all pregnancies [1]. The onset of this complication is after 20 weeks of gestation that originating at the maternal–fetal interface with obvious signs of multi-systemic involvement, from hypertensive disorders (systolic blood pressure (BP) ≥ 140 mm Hg or diastolic BP ≥ 90 mm Hg) and proteinuria (≥ 300 mg/24 h) till greater end-organ damage including hepatic alterations (hemolysis, elevated liver enzymes and low platelets (HELLP), placental structure and function impairment and subsequent increased oxidative stress and inflammation, pulmonary and neurological dysfunction including seizures [2,3,4,5,6]. Although understanding of mechanisms contributing to the pathophysiology of preeclampsia is an active area of international research, the etiology of preeclampsia remains unknown and except prophylactic use of low-dose aspirin (150 mg per day from 11 to 14 weeks of gestation until 36 weeks) in women who being identified at high risk for preterm preeclampsia, and preterm delivery that resolve signs and symptoms of preeclampsia no definitive treatment is available for preeclampsia [1, 7, 8]. Despite the exact causes of these alterations remain unclear, but it is well-established that unlike healthy pregnancy that is a controlled inflammatory process, preeclampsia is a proinflammatory state [8, 9]. Studies suggested a complex relationship between the pro-inflammatory condition of preeclampsia and immunologic aberrations. Evidence indicates that among different immune molecules, respective cytokines of T helper 1 (Th1) and T helper 2 (Th2) cells play an important role in controlling immune system function in various stages of pregnancy. In normal pregnancy, pro-inflammatory Th1 cytokines such as; interleukin 2, interferon gamma, transforming growth factor beta and angiotensin II type 1 receptor, which are involved in cellular immunity and mediate immune rejection of the fetus are decrease, by contrast, Th2 cytokines: interleukins 4, 5, 13 and 19 that involved in the maintenance of normal pregnancy occurrences, helps to neutralize pro-inflammatory cytokines, increased [10]. Imbalanced concentration of pro- and anti‐inflammatory cytokines likely contributes to the preeclampsia related complications [11]. Since the single nucleotide substitutions (SNPs) in the cytokine genes may affect cytokine transcription and influence its production [12], in our study, the association between preeclampsia and IL-13 gene rs20541 (G to A exchange) and rs56035208 (G to C exchange) polymorphisms and IL-19 rs1028181 (T to C exchange) and rs2243191(T to C exchange) have been assessed.

Methods

Study design

In this case–control study, totally 300 participants, including 150 preeclamptic and 150 healthy pregnant women who were referred to Zeynabieh Hospital- Shiraz, Iran, during February 2021 to April 2022 enrolled. It is worth noting that this study conducted in a timeframe that a population under 1,000 of pregnant women had registered for medical care under the medical assistance insurance (Medicaid). Therefore, a minimum ratio of 30% (300 individuals) ensured representativeness of the sample. Demographic characteristics are shown in Table 1. The hallmark presentations of preeclampsia, blood pressure; systole > 140 mm Hg and diastole > 90 mm Hg on two occasions at least 4 h apart (the mild form of preeclampsia) or a shorter interval timing of systole > 160 mm Hg and diastole > 110 mm Hg (the severe form of preeclampsia) after 20 weeks’ gestation with accompanying proteinuria ≥ 300 mg upon a 24 h urine collection; or UPCR ≥ 0.3 mg/mg; or a urine protein dipstick reading ≥ 2+, were considered as definitive criteria for preeclamptic women. However, inclusion criteria were expanded by additional evidence of other significant findings that may be a part of the clinical presentation including end organ damage consisting of impaired liver function, severe persistent epigastric pain, new-onset headache, pulmonary edema, or renal insufficiency and abnormal lab values. Finally, women with high blood pressure accompanied with proteinuria, who have met the criteria of other organ involvement were included. The controls were collected from the healthy pregnant women, without hypertension, with the same gestational age (≥ 20 weeks), who had been referred to the midwifery clinic of the same hospital to receive prenatal care within the same timeframe. Pregnant women with a history of chronic hypertension and current antihypertensive treatment and any underlying disease such as renal, autoimmune, metabolic, liver, diabetic neuropathy or cardiovascular disease were excluded from our study (Fig. 1).

Fig. 1

Flow chart of the selection procedure

Table 1 Demographic characteristics of participants

DNA extraction and genotype analysis

DNA was extracted from 5 ml peripheral blood, collected in EDTA tubes by Genomic DNA Isolation Kit (GeNet Bio, Daejeon, Korea), using salting out method. Allele-specific polymerase chain reaction (AS-PCR) (also known as amplification refractory mutation system (ARMS)) was designed for detection of rs20541 (G/A), rs56035208 (G/C), rs1028181 (T/C) and rs2243191(T/C) polymorphisms. For each variation, a pair of control primers specific for the normal DNA sequence (A pair of control primers which could not amplify mutant DNA at a given locus was used to confirm that the genomic DNA is, in principle, amplifiable) and allele-specific primers were designed using Oligo7 software (version 7.54, Molecular Biology Insights Inc., Cascade, CO, USA). (The designed primer sequences reported in Table 2).

Table 2 Designed primers for ARMS-PCR reactions and annealing temperature

PCR was performed in a total volume of 23.5 µl, including: 1 µl of template DNA, 12.5 µl of 2× Master Mix Red (Ampliqon), 1 µl (0/39 pmol) each primer, 5 µl H₂O. The cycling conditions was as follows: first denaturation; 5 min at 94 °C, followed by 30 cycles at 94 °C for 40s (Denaturation), annealing temperature as in Table 2 for 40s and 1 min at 72 °C (Extension) and a final extension for 7 min at 72 °C. PCR products were run by standard electrophoresis on 2% agarose gel for 10 min and visualized on UV transilluminator. Bands length presented in Fig. 2.

Fig. 2

The pattern of bands observed after agarose-gel electrophoresis. IL-13; with control primers, a 169bp amplification product was obtained in two reactions. A 223bp (rs20541) and a 221bp (rs56035208) amplification products, are corresponding to the allele -specific primers used in AS-PCR. IL-19; with control primers, a 166 bp and a 169 bp amplification product were obtained in rs2243191 and rs1028181 reactions, respectively. A 262bp (rs2243191) and a 372bp (rs1028181) amplification products, are corresponding to the allele -specific primers used in AS-PCR. The gels displayed here not cropped, and are without high-contrast (overexposure)

Direct sequencing of PCR products recovered by the GEL/PCR Purification Kit (Favorgen Biotech Corp., Taiwan, China) was performed using Genetic Analyzer 3130x (Applied Biosystems, USA). Sequences were analyzed with the CodonCode Aligner V.5.1.5 software (CodonCode Corporation, Centerville, MA, USA).

Statistical analysis

Data was analyzed by using, social sciences statistical software package for windows, version 19.0 (SPSS Inc, Chicago, Illinois, USA). The differences of the alleles and genotypes frequency were compared between preeclampsia and control groups by chi-square test. P value < 0.05 was considered statistically significant.

Results

We first evaluated the Hardy–Weinberg equilibrium by computing expected genotype values versus observed genotype values for all polymorphic loci to check whether the population was in Hardy–Weinberg equilibrium. The results showed that the deviation from Hardy–Weinberg equilibrium in the all polymorphic loci were not significant, therefore equilibrium was maintained for in question population at polymorphic rs20541 (G/A), rs56035208 (G/C), rs1028181 (T/C) and rs2243191(T/C) sites.

The accuracy and specificity of our established AS-PCR was further validated by direct sequencing of PCR products. Figure 3 shows the sequencing result for samples.

Fig. 3

Result of DNA sequencing of PCR products of IL-13 and IL-19 polymorphisms

As Table 3 displays the allele and genotype frequencies for IL-13 and IL-19 polymorphisms in the studied groups, no significant difference was observed in patients and controls groups allele (P = 0.8) and genotype (P = 0.4) frequencies for rs20541 (G/A) and allele (P = 0.2) and genotype (P = 0.6) frequencies of rs56035208 (G/C). The obtained results of the allele and genotype frequencies for the IL-19 gene SNPs showed statistically significant differences in comparisons of controls, with allele (P = 0.002), and genotype (P = 0.001) values for rs1028181, and allele (P = 0.042), and genotype (P = 0.032) values for rs2243191.

Table 3 Genotypes and alleles frequency of the IL-13 and IL-19 polymorphisms in preeclampsia and control groups

To further identify a genetic risk for preeclampsia, we tested the association between preeclampsia risk markers such as seizure, proteinuria, edema, diabetes, multipara, hypertension, gravida, hyperthyroidism, HELLP syndrome, and history of abortion, history of preeclampsia, maternal drinking, maternal smoking, socioeconomic position, maternal BMI, birth weight and maternal age and rs20541 (G/A), rs56035208 (G/C), rs1028181 (T/C) and rs2243191(T/C) polymorphisms. In the case of IL-13 SNPs, except for edema at rs20541 position (p = 0.02) and maternal drinking at rs56035208 position (p = 0.002) no significant difference between the two groups for other risk markers observed (Table 4). Yet, the IL-19 SNPs were found to be associated with the more risk markers, such as; seizure (p = 0.005 _rs1028181, 0.003 _rs2243191), proteinuria (p = 0.023 _rs1028181), edema (p = 0.017 _rs2243191), diabetes (p = 0.023 _rs2243191), hyperthyroidism (p = 0.015 _rs1028181, 0.032 _rs2243191), HELLP syndrome (p = 0.041 _rs2243191), maternal smoking (p = 0.022 _rs1028181) and maternal BMI (p = 0.045 _rs1028181, 0.012 _rs2243191). Significant risk markers are shown in Table 5.

Table 4 Association between the IL-13 (rs20541 (G/A) and rs56035208 (G/C)) polymorphisms and preeclampsia risk markers in preeclampsia and control groups

Table 5 Association between the IL-19 (rs1028181 (T/C) and rs2243191 (T/C)) polymorphisms and preeclampsia risk markers in preeclampsia and control groups

Discussion

Interleukin-13, that belongs to T helper 2 family and produced classically by T helper cells and also by eosinophils, basophils, mast cells, and NK cells, is a peptide cytokine coded by the gene IL‐13 in the chromosomal location of 5q31.1, in a cluster close to IL‐3, IL‐4, IL‐5, and CSF2 [13]. Since genetic factors are believed to be involved in the development of preeclampsia, and existing data confirm the association of IL‐13 polymorphism with inflammatory states and disorders, e.g., association of SNP IL‐13 + 1923 C/T with asthma development [14], or association of IL13 rs20541 with significantly decreased susceptibility to renal cell carcinoma (a disease with immune predisposition) [15], our study aimed to analyze the possible impact of IL-13 gene polymorphism with the susceptibility to preeclampsia. Two possible polymorphic sites of IL-13 were identified, including rs20541 (G/A) and rs56035208 (G/C). The results revealed that none of the allele or genotype frequencies of IL-13 rs20541 (G/A) and rs56035208 (G/C) polymorphisms had any significant association with preeclampsia susceptibility. Though very little is known about the IL‐13 gene polymorphism on the development of preeclampsia, the impact of several other polymorphic sites have been reported. De Lima and colleagues did not find an association between polymorphisms in genes, TNF-α (− 308 G > A), IL6 (− 174 G > C), IFN-γ (+ 874 A > T), IL10 (− 1082 A > G), (− 819 C > T) and (− 592 C > A) and TGF- (+ 869 T > C) with preeclampsia [16]. In Pinheiro article also non-association between TNF-α (− 308 G → A), IL-6 (− 174 G → C), or IL-10 (− 1082 G → A) polymorphisms and preeclampsia was reported [9]. No significant difference in the distribution of genotypes and alleles of IL-10 G-1082 A between the two groups was observed in the Sowmya study [17]. Tanaka and colleagues found no association between the rs4711998 A > G, rs8193036 C > T and rs2275913 A > G polymorphisms of the IL17A gene and the risk of developing preeclampsia [18]. In addition, findings of the current study, except edema in the position of rs20541 and maternal drinking at rs56035208 position, showed no significant association between other preeclampsia major risk factors include; seizure, proteinuria, diabetes, multipara, hypertension, gravida, hyperthyroidism, HELLP syndrome, history of abortion and history of preeclampsia, maternal drinking, maternal smoking, socioeconomic position, maternal BMI, birth weight and maternal age with rs20541 and rs56035208 polymorphisms. Similarly no significant relation in Lisi study, was found between preeclampsia risk factors, such as systolic and diastolic blood pressure, proteinuria and endothelin-1 type A receptor gene polymorphism (7231 G4A) [19]. Molvarec observed no significant relation between the TNF-alpha G-308 A polymorphism and HELLP syndrome [20]. PFAB and co-worker also reported that carrying C allele of the interleukin-6 polymorphisms (interleukin-6 S174G > C) significantly associated with new-onset edema in pregnant women [21].

On the contrary, the present study revealed association of IL-19 gene polymorphisms in the preeclampsia group compared to the controls. Interleukin-19 is a member of the IL-10 family which is known as a pregnancy compatible cytokine [22]. The association of IL-19 SNPs and preeclampsia has remained incompletely understood, but by increasing numbers of studies, strong evidence for the biological significance of IL-19 thus far has been provided which indicates its significance in asthma [23,24,25,26], psoriasis [27], rheumatoid arthritis [28] and other immune diseases. Importantly, we have identified roles for IL-19 SNPs as a genetic determiner that is associated with seizure, proteinuria, edema, diabetes, hyperthyroidism, HELLP syndrome, maternal smoking and maternal BMI during pregnancy.

Broadly speaking, new markers such as genetic characteristics, imaging, and biomarkers, are useful to improve risk prediction for patients and hold the promise of improving clinical prediction models. Based on the currently acknowledged risk factors, different professional organizations have thus far developed a number of prediction models for preeclampsia that in most models maternal serum PAPP-A and PlGF used as two biochemical predictors. Therefore, improving prediction models by maternal risk traits related markers (like maternal genetic variation considered in our study), with good discrimination and calibration, can show promising results in the early prediction of preeclampsia [29, 30]. Therefore, we conclude and recommend that IL-19 SNPs with high detection rates or low false-positive rates can be considered for improving prediction models at the time of model updating. As the available sample with the new marker in this study is small, our recommendation is using parsimonious methods for extending an existing prediction model from it that lead to the largest increase in discriminative ability of the prediction model [31].

Considering a larger sample size to evaluate association these in question polymorphisms with preeclampsia in future studies will resolve the present study limitation. Likewise, further studies on patients in a variety of ethnic populations are still required to increase our knowledge-base for these genes.

Conclusion

Focused on polymorphisms and subtly investigating their potential role in preeclampsia is biologically or clinically important and might prove a reliable tool to identify and manage the women at risk.