Abstract
Aims: Th1 and Th2 balance is crucial for maintenance of pregnancy, and intracellular JAK and STAT proteins significantly contribute to it. In view of evidence linking JAK2 and STAT3 variants with recurrent spontaneous miscarriage (RSM), here we investigated the association of JAK2 (rs2230724) and STAT3 (rs1053023 and rs1053004) to RSM susceptibility in Tunisians. Subjects and Methods: A retrospective case-control study. Subjects comprised 235 RSM cases and 235 control subjects. JAK2 and STAT3 were genotyped by the allelic discrimination method. Results: STAT3 rs1053023 and, to a lower extent, rs1053004 were significantly associated with RSM under the additive and dominant, but not recessive models. This remained significant after adjustment for the covariates age, smoking, and gravida. In contrast to STAT3 variants, JAK2 rs2230724 was not associated with RSM under any of the genetic models tested. Two-locus STAT3 (rs1053023/rs1053004) haplotype analysis revealed increased frequency of the C/G haplotype in patients with RSM. Multivariate regression analysis confirmed the association of C/G haplotype with RSM (p=0.001; odds ratio=2.01; 95% confidence interval=1.32-3.07), thus conferring RSM susceptibility nature. These differences remained significant after applying the Bonferroni correction for multiple testing (Pc=0.004). Conclusions: STAT3 rs1053023, more so than the STAT3 rs1053004 or JAK2 rs2230724 polymorphisms, is associated with RSM risk.
Introduction
I
Antigen-specific T cell activation is a complex process, which converges at engagement of JAK kinases and the STAT family of intracellular proteins (Hauguel-de Mouzon and Guerre-Millo, 2006), of which the JAK2-STAT3 signaling pathway was reported to play an important role in trophoblast invasion and differentiation (Mulla et al., 2010), and is required for establishment and maintenance of pregnancy (Bao et al., 2006; Fitzgerald et al., 2008), as shown by the finding that dysregulated STAT3 expression and activation were associated with defective placentation (Tarrade et al., 2002; Catalano et al., 2005; Mulla et al., 2010). Since JAK2 and STAT3 expression are genetically determined, previous studies on the association of the JAK2 and STAT3 polymorphisms with RSM was explored, but with inconclusive findings (Mercier et al., 2007; Finan et al., 2010; Kaandorp et al., 2010). We evaluated the association of the JAK2 (rs2230724) and STAT3 (rs1053023 and rs1053004) gene variants with RSM in Tunisian subjects, and found that of these, STAT3 (rs1053023) is the only variant identified in association with RSM in our population. These three variants are single-nucleotide polymorphisms (SNPs). JAK2 rs2230724 is an exonic polymorphism, localized in exon 19. This variant was associated with susceptibility to acute leukemia and its subtypes, thus suggested to be a genetic determinant for acute leukemia and acute myeloid leukemia (Zhong et al., 2012). Both STAT3 variants rs1053023 and rs1053004 are localized in exon 24. STAT3 rs1053023 is associated with 28% risk reduction for B-cell non-Hodgkin's lymphoma (Butterbach et al., 2011), while STAT3 rs1053004 is significantly associated with idiopathic RSM (Finan et al., 2010).
Subjects and Methods
RSM cases comprised 235 unrelated women (mean age 32.1±6.1), who were recruited from the outpatient OB/GYN clinics in for assessment of idiopathic RSM, as per the Royal College of Obstetricians and Gynecologists (RCOG) Guidelines (Table 1). The inclusion criteria included ≥3 consecutive pregnancy losses with no known risk factors, and with the same partner, which occurred during the first trimester of gestation. Exclusion criteria were Rh incompatibility, older age at first miscarriage (≥40 years), preclinical miscarriages and/or biochemical pregnancy, and pre-eclampsia (elevated systolic and diastolic blood pressure (BP) >145/95 mm Hg, or rise in systolic/diastolic BP>30/15 mmHg on at least two occasions). Additional exclusion criteria included systemic autoimmune disease, thyroid dysfunction, diabetes mellitus, anatomical disorders, infections (toxoplasmosis, human cytomegalovirus, rubella, HIV, Group B streptococci, Chlamydia trachomatis, hepatitis B and C, and bacterial vaginosis), and liver function abnormalities.
Tunisian subjects consisted of 235 cases and 235 control women.
Student's t-test for continuous variables; Pearson's chi-square tests for categorical variables.
Number of subjects (percent total).
NA, not applicable.
Controls comprised 235 multiparous ethnically matched women who had at least two live births, and no miscarriages, and consisted of university and hospital employees, or volunteers. All participants were Tunisian Arabs, and non-Arab Tunisians (Berber) were excluded. On all subjects, demographic details were recorded, and were verified from clinic records where available. The Research and Ethics Committee of the University of Monastir (Tunisia) approved the study protocol. Blood samples were taken from all participants in ethylenediaminetetraacetic acid-containing tubes by venipuncture for genomic DNA extraction, after signing a consent form agreeing to be included in the study.
Genotyping
JAK2 (rs2230724) and STAT3 (rs1053023 and rs1053004) polymorphisms were genotyped using the allelic (VIC- and FAM-labeled) discrimination method. TaqMan assays were ordered from Applied Biosystems (Table 2). The reaction was performed in a 6-μL volume on StepOne real-time PCR system, as recommended by the manufacturer (Applied Biosystems). Replicate-blinded quality control samples were included to assess reproducibility of the genotyping procedure; concordance was >99%.
Location on chromosome based on dbSNP-build 125.
p determined using Haploview 4.2.
SNP, single-nucleotide polymorphism; HWE, Hardy-Weinberg equilibrium.
Statistical analysis
Data were expressed as mean±SD (continuous variables), or as percentages of total (categorical variables), and intergroup significance was assessed by Student's t-test (continuous variables), and the χ2 test (categorical variables). Allele frequencies were calculated by the gene-counting method; each SNP was tested for the Hardy-Weinberg equilibrium (HWE) using χ2 goodness-of-fit test by Haploview 4.2 (www.broad.mit.edu/mpg/haploview/). Genotypic association of each SNP with RSM susceptibility was conducted under additive, dominant, and recessive models using SNPstats (http://bioinfo.iconcologia.net/SNPstats), taking homozygous wild-type genotype as the reference (odds ratio [OR]=1.00), after controlling for age, smoking, and body-mass index (BMI).
Results
Study subjects
The demographics and clinical characteristics of study subjects are shown in Table 1. BMI (25.2±3.9 kg/m2 vs. 25.6±4.1 kg/m2; p=0.322) and the prevalence of obesity (>30 kg/m2) (p=0.325) were comparable between RSM cases and control women. In contrast, the prevalence of smokers (p=0.002), age (p=0.012), and gravida (p<0.001) was significantly different between RSM cases and control women. Accordingly, the latter were the covariates that were controlled for in subsequent analysis. In total, RSM cases reported a mean of 3.5 pregnancy losses, while control women reported an average of 3.2 live births/woman and no pregnancy losses.
JAK2 and STAT3 genotype distribution
The distribution of JAK2 rs2230724 (p=1.00), STAT3 rs1053023 (p=0.066), more so than rs1053004 (p=0.023) genotypes was in the HWE among study subjects (Table 2). Table 3 summarizes the minor allele frequencies (MAF) of the JAK2 and STAT3 variants in RSM cases and control subjects. MAF of both JAK2 rs2230724 (p=0.44) and STAT3 rs1053004 (p=0.10) were not significantly different between RSM cases and control women. In contrast, significantly higher MAF of STAT3 rs1053023 were seen in cases (p=0.006), and translated into increased RSM risk [OR (95%)=1.50 (1.13-1.97)].
Study subjects comprised 235 RSM cases and 235 control women.
Crude p value.
aOR, adjusted OR, adjusted for gender and BMI.
Bold numbers indicate statistically significant values.
MAF, minor allele frequency; OR, odds ratio; CI, confidence interval; RSM, recurrent spontaneous miscarriage.
Table 4 summarizes the association analyses between JAK2 rs2230724 and STAT3 rs1053023 and rs1053004 variants, under additive, dominant, and recessive genetic models. STAT3 rs1053023 and, to a lower extent, rs1053004 showed a significant association with RSM under the additive and dominant, but not recessive models, which remained significant after adjustment for the covariates age, smoking, and gravida. In contrast to the STAT3 variants, JAK2 rs2230724 did not associate significantly with RSM under any of the genetic models tested.
p values adjusted for age, gender, and body-mass index.
Number of subjects (frequency).
Bold numbers indicate statistically significant values.
STAT3 haplotype analysis
Two-locus haplotypes were constructed based on the prevalence of individual SNPs, and linkage disequilibrium between them. Taking the common rs1053023/rs1053004 T/A haplotype as reference (OR=1.00), multivariate regression analysis confirmed the association of the C/G haplotype with RSM (p=0.001; OR=2.01; 95% confidence interval [CI]=1.32-3.07), thus conferring RSM susceptibility nature (Table 5). These differences remained significant for C/G (Pc=0.004) after applying the Bonferroni correction for multiple testing (Table 5).
STAT3 haplotype (rs1053004 (T/C)/rs1053023 (A/G)) frequency determined by the maximum likelihood method.
Pearson's chi-square test.
Haplotype frequency.
Bold numbers indicate significant associations.
Discussion
Several lines of reasoning implicated altered JAK2/STAT3 activity with adverse pregnancy outcomes. The beneficial effect of Interleukin-11 for normal placentation in pregnancy involves upregulation of the α2-macroglobulin gene through the JAK2/STAT3 pathway (Bao et al., 2006), and JAK2 is a key factor for pregnancy-associated mammary epithelium proliferation and differentiation (Shillingford et al., 2002). In addition, STAT3 knockouts display embryonic lethality, at the gastrulation phase, further establishing a central role for STAT3 in successful pregnancy (Takeda et al., 1997), as evidenced by the findings that decidual cytokines utilize STAT3 as a common intracellular signaling protein, thereby allowing for effective invasion, proliferation, and differentiation of trophoblast cells, hence implantation (Catalano et al., 2005; Fitzgerald et al., 2008; Pu et al., 2012). Accordingly, dysregulation of STAT3 expression and activation during placental development are likely to provoke defective placentation, along with shallow invasion in pre-eclampsia or neoplasms such as human choriocarcinoma (Tarrade et al., 2002).
We analyzed the possible association between the JAK2 and STAT3 variants and RSM risk among Tunisian women. Among the acquired (modifiable) risk factors, RSM risk was linked with the prevalence of smoking, which was in apparent disagreement with a recent study on Chinese women reporting no such association (Zhang et al., 2010). This discrepancy may be related to the coffee consumption in our population compared to the Chinese cohort, which was reported to contribute to increased RSM risk in women smokers (Derbyshire and Abdula, 2008). Significant association was noted between RSM and STAT3 rs1053023, but not with JAK2 rs2230724 or STAT3 rs1053004, in apparent disagreement with the recent study of Finan, in which STAT3 rs1053004 more so than STAT3 rs1053023 was associated with RSM risk in Bahraini Arab women (Finan et al., 2010). This apparent contradiction is likely reconciled by differences in ethnic background of these distinct communities, as was earlier reported (Almawi et al., 2005).
We analyzed the contribution of the JAK2 and STAT3 variants to RSM risk at the allele and genotype levels, the latter, including the additive, dominant, and recessive models. JAK2 rs2230724, STAT3 rs1053023 A>G, and rs1053004 C>T minor allele frequencies among control subjects were similar to frequencies established for Europeans/Caucasians and Asians, were higher that the rates seen for Asians (Japanese and Chinese), but were lower than frequencies established for African-Americans, or Sub-Saharan Africans. These findings confirm racial heterogeneity in the distribution of the JAK2 and STAT3 variants, and hence their disease association.
Since both STAT3 variants (rs1053004 and rs1053023) were in LD, 2-locus haplotypes were constructed based on the prevalence of the individual SNP and the LD between them. Of the three variants tested, slight deviation from the HWE was noted for rs1053004, which is not likely to be due to genotyping errors, but rather to population stratification, as was suggested elsewhere (Hosking et al., 2004; Vine and Curtis, 2009).
We found no significant association between JAK2 rs2230724, previously linked and altered cell proliferation and differentiation (Zhong et al., 2011), and RSM risk in both populations. Mercier reported that the JAK2 V617F mutation, which occurs in 0.20% of normal Southern European primigravidas 21-36 years of age, is linked with increased risk of pregnancy loss (Mercier et al., 2007), while Kaandorp reported that the V617F variant is not a clinically relevant etiologic factor in unexplained RSM (Kaandorp et al., 2010). Melillo et al. (2009) show that the fetal loss rate in our study was still high (24.6%). In particular, the rate of fetal loss due to spontaneous abortion in the first trimester was 15.6% in this study, while it ranged between 28.1% and 31.7% in other reports (Passamonti et al., 2007; Gangat et al., 2009). The risk of spontaneous abortion was 2.5 times higher in patients with ET than in the age-matched general population in the study of Melillo et al. (2009). Coupled to other studies demonstrating lack of association of the related variant V617F (rs77375493) with pregnancy complications (Dahabreh et al., 2009; Grandone et al., 2011), this questions the contribution of altered JAK2 activity in pregnancy complications.
PSTAT3 has been previously shown to be involved in fertility, implantation, and decidualization, hence the maintenance of pregnancy (Garcia et al., 2007). Altered STAT3 activity was linked with poor pregnancy outcome, evidenced by the reduced endometrial pSTAT3 activity in women with unexplained infertility (Dimitriadis et al., 2007). We focused on the 3′UTR region containing regulatory elements for STAT3 expression, and reasoned that altered STAT3 activity induced by specific STAT3 variants influences RSM susceptibility, thus contributing to pregnancy loss by setting in local inflammation at placental membranes (Dimitriadis et al., 2006), leading to accelerated miscarriage. At the haplotype level, the double-mutant C/G haplotype was associated with increased RSM risk (OR=2.01, 95% CI (1.32-3.07), after controlling for several covariates.
In conclusion, our study demonstrates that STAT3 rs1053023, but not STAT3 rs1053004 or JAK2 polymorphisms, is associated with an RSM risk. Our study has several strengths, namely the homogeneity of the population tested (only Tunisian Arabs from central Tunisia), and limiting the RSM cases to early miscarriages, thus minimizing the possibility of patient misclassification. Our study has some shortcomings, namely that we did not correlate the STAT3 genotypes with STAT3 activity, and to the selection of the SNPs, thus raising the possibility of the contribution of other variants in neighboring and distant genes in modifying RSM risk. A larger study is needed to confirm the current findings to reveal the mechanism by which altered STAT3 activity induced by specific genotypes influences overall RSM risk.
Footnotes
Author Disclosure Statement
No competing financial interests exist.
