{"paper_id":"efb69805-211e-4b8d-9365-ad3179ec5151","body_text":"Recurrent pregnancy loss (RPL) or recurrent spontaneous abortion has been defined as the occurrence of at least two consecutive pregnancy losses prior to the 20th week of gestation [ 1 ,  2 ]. RPL occurs in approximately 1% of all pregnancies; however, the etiology for more than half of the RPLs remains undetermined [ 3 ]. Genetic variation has been suggested one of the contributing factors leading to RPL and a number of single nucleotide polymorphisms (SNPs) have been reported to be associated with RPL [ 4 ]. MicroRNAs (miRNAs) are short (approximately 22 nt) noncoding RNA molecules regulating expression of target genes at the post-transcriptional level by translational repression or messenger RNA degradation [ 5 ]. Several studies recently reported the associations between miRNA polymorphisms and RPL [ 6 – 9 ]. One study identified two SNPs in miR-125a altering the production of miR-125a which was subsequently associated with an elevated risk for RPL in the Han Chinese women [ 8 ]. Another study reported an association between two pre-miRNA polymorphisms (miR-196a2 and miR-499) and the occurrence of RPL in Korean females [ 9 ], which was supported in Iranian women [ 6 ]. The most recent study identified a polymorphism in the coding region of  miR-423  contributing to an increase in the expression of mature  miR-423  associated with RPL in the Han Chinese population [ 7 ]. Several miRNAs that are considered important during pregnancy were chosen for this study because of their elevated expression ( miR-27a ), decreased expression in the endometrium and in trophoblasts ( miR-423 ), lower expression during endometriosis ( miR-449b ), and involvement in pregnancy loss via the p53 network ( miR-605 ) [ 10 – 13 ]. In this study, we determined the susceptibility to RPL associated with genetic variants of miRNAs associated with placental or fetal development.\n\nThe study group consisted of 387 females, 33.21 ± 4.55 years of age [mean age ± standard deviation (SD)] and a body mass index (BMI) of 21.49 ± 3.84 (± SD) who were diagnosed as idiopathic RPL patients with at least two consecutive pregnancy losses prior to the 20th week of gestation according to the definitions of infertility and recurrent pregnancy loss by American Society for Reproductive Medicine [ 1 ]. These patients were enrolled in a study at the Infertility Medical Center of CHA Bundang Medical Center from March 1999 to February 2012. Among the RPL patients, none had a history of smoking or alcohol use. RPL patients due to anatomical, hormonal, chromosomal (patients or their spouses), infectious, autoimmune, or thrombotic causes are excluded from the study. The age-matched control group consisted of 225 females, 33.43 ± 5.89 years of age (mean age ± SD) and with a BMI of 21.68 ± 3.451 (± SD), each of whom had regular menstrual cycles, had a history of naturally conceived pregnancy at least once, had no history of pregnancy loss or karyotype 46,XX, and who were recruited from the CHA Bundang Medical Center. All patients and controls were Korean. The institutional review board of CHA Bundang Medical Center approved the study, and all patients provided written informed consent.\nGenomic DNA was extracted from non-coagulated peripheral blood using the G-DEX blood extraction kit (Intron, Seongnam, Korea). The nucleotide changes were determined by the polymerase chain reaction (PCR)-restriction fragment length polymorphism analyses using the isolated genomic DNA as a template. Primer sequences for PCR amplification of each polymorphism were as follows:  miR-27a A>G [rs895819], forward  5'-GAA CTT AGC CAC TGT GAA CAC  C AC TTG G-3'  and reverse  5'-TTG CTT CCT GTC ACA AAT CAC ATT G-3'  (the mismatch sequence is underlined) [ 14 ];  miR-423 C>A [rs6505162], forward  5'-GTA CAT TTT CCC GGA TGG AA-3'  and reverse  5'-GGG AGA AAC TCA AGC GC C  G-3' ;  miR-449b A>G [rs10061113], forward  5'-GGT ATC CAG AGC ACT TCA TTG ACA-3'  and reverse  5'-ACC TGA ATC AGG TAG GCA GTG T C T-3' ; and  miR-605 A>G [rs2043556], forward  5'-AGA GCA GTT ACG CCA CAT GA-3'  and reverse  5'-GCC TTC TCC TTG GGA AAA AC-3' . We performed a restriction enzyme digestion at 37°C for 16 hours using  Dra III (New England BioLabs, Ipswich, MA, USA) for the  miR-27a  polymorphism,  Bsr FI for the  miR-423  polymorphism,  Bsm AI for the  miR-449b  polymorphism, and  Hin fI for the  miR-605  polymorphism. We confirmed the genotyping of the four sites by sequencing 10% of the samples.\nBlood samples from RPL patients were collected during pregnancy. Plasma homocysteine, folate, total cholesterol, and urate concentrations, and blood coagulation factors were measured in RPL patients after fasting for 12 hours. Homocysteine levels (6.98 ± 2.10 μM) were measured using a fluorescence polarization immunoassay and the Abbott IMx analyzer (Abbott Laboratories, Abbott Park, IL, USA). Folate levels (14.21 ± 11.94 ng/mL) were determined using a competitive immunoassay with ACS:180 (Bayer Diagnostics, Tarrytown, NY, USA). Total cholesterol (187.73 ± 49.42 mg/dL) and urate levels (3.80 ± 0.84 mg/dL) were determined using commercially available enzymatic colorimetric tests (Roche Diagnostics, Mannheim, Germany). Platelet (PLT) counts, prothrombin time (PT), and activated partial thromboplastin time (aPTT) were measured to assess blood coagulation. PLT counts (255.43 ± 59.22 10 3  cells/μL) were measured using a Sysmex XE2100 automated hematology analyzer (Sysmex, Kobe, Japan). PT (11.58 ± 0.85 seconds) and aPTT (32.24 ± 4.33 seconds) were measured using an automated photo-optical coagulometer (ACL TOP; Mitsubishi Chemical Medience, Tokyo, Japan).\nPeripheral blood mononuclear cells (PBMCs) were isolated from whole blood using a cell preparation tube containing sodium citrate (Becton-Dickinson, Franklin Lakes, NJ, USA). To obtain monocytes, viable PBMCs were frozen in 80% fetal bovine serum (FBS; Lonza, Cologne, Germany), 10% dimethyl sulfoxide (Sigma-Aldrich, St. Louis, MO, USA), and 10% RPMI 1640 media (Life Technologies, Carlsbad, CA, USA) in liquid nitrogen. After thawing, the PBMCs were cultured in RPMI 1640 media supplemented with 10% FBS, 50 mg/mL gentamicin sulfate (Lonza), 50 μM 2-mercaptoethanol (Sigma-Aldrich), and 2 mM glutamine (Life Technologies). The cells were washed twice with phosphate buffered saline (Welgene, Seoul, Korea) and then resuspended in RPMI 1640 media containing 10% FBS, 1% Minimal Essential Media with nonessential amino acids (Life Technologies), and 1% sodium pyruvate (Life Technologies) at a density of 1 × 10 6  cells/mL, and incubated overnight as described above. All NK cell assays were performed after 16–20 hours of incubation.\nTo determine the absolute number of NK cells, 200 mL of diluted blood was incubated for 20 minutes on ice with phycoerythrin-conjugated anti-CD56 and peridinin chlorophyll protein-conjugated anti-CD3 monoclonal antibodies (BD Biosciences, San Jose, CA, USA). Then, 20,000 fluorescein isothiocyanate conjugated beads were added, and the blood sample was subjected to erythrocyte lysis using FACS Lysing Solution (BD Biosciences). The samples were analyzed on a flow cytometer using the BD FACSCalibur (BD Biosciences). The NK cell number of the diluted blood samples was calculated as (NK cells/mL sample) = [(CD56 + /CD3 -  cell count)/(bead count)] × 100. This FACS-based NK cell count required 0.1 mL of whole blood per tested condition.\nThe differences in four microRNAs ( miR-27a A>G,  miR-423 C>A,  miR-499b A>G, and  miR-605 A>G) genotype and haplotype frequencies between patients and normal controls were compared using Fisher’s exact test and logistic regression analyses. Allele frequencies were estimated to identify deviations from the Hardy-Weinberg equilibrium (HWE). Adjusted odds ratios (AOR) and 95% confidence interval (CI) were estimated as a measure of the strength of association between genotypes and RPL risk. Multiple comparison tests were adjusted by using the false discovery rate (FDR) correction, and associations with an FDR-adjusted  P  value < .05 were considered statistically significant [ 15 ].\nGene-gene interactions among SNP loci were analyzed with multifactor dimensionality reduction (MDR) using MDR software, version 2.0 ( www.epistasis.org ) [ 16 – 18 ]. Based on the MDR identification of the most significant models using the best maximized cross validation value, the best multilocus combinations were determined. HAPSTAT software, version 3.0 ( www.bios.unc.edu/Èlin/hapstat/ ) was used to estimate haplotype frequencies for polymorphisms that were determined by MDR analyses to have strong synergistic effects. Statistical analyses were performed using GraphPad Prism software, version 4.0 (GraphPad, San Diego, CA, USA) and StatsDirect software, version 2.4.4 (StatsDirect, Altrincham, UK). The statistical significance was set at  P  <0.05. The relationship between each of the four microRNA genotypes and each of the six clinical variables of the RPL patients (plasma homocysteine, folate, NK cell, PLT, PT, and aPTT which may contribute to an even balance of coagulation and fibrinolysis during pregnancy) was analyzed using multiple linear regression analyses. Regression models were examined for the six clinical variables with grades (10 levels) as dependent variables and the microRNA polymorphisms as independent variables. The differences in plasma homocysteine, folate, PLT, PT, aPTT, and NK cells as a function of the four microRNA genotypes and combination genotypes were evaluated by one-way analysis of variance and independent two-sample  t -tests.\n\nThe demographic characteristics and clinical profiles of RPL patients and control subjects are shown in  Table 1 . The two groups were matched for age and BMI. Platelet numbers were significantly higher in the patient group than in the control group. The genotype and allele frequencies of the four miRNA SNPs in females with RPL and controls are shown in  Table 2 . All genotypes in the study cases were in HWE. In  Table 2 , the  miR-27a  A>G polymorphism was significantly associated with a risk of RPL (AA vs. AG: AOR = 0.654; 95% CI = 0.456–0.937; AA vs. AG+GG: AOR = 0.682; 95% CI = 0.484–0.960); however, there was no association after adjustment for multiple tests using the FDR correction. When RPL patients were stratified according to the occurrence of consecutive recurrent pregnancy losses (RPL = 2 and ≥ 3 vs. all RPL patients with RPL ≥ 2), there was an association between the  miR -27a  A>G polymorphism and RPL risk in the RPL ≥ 3 subgroup alone (AA vs. AG: AOR = 0.611; 95% CI = 0.404–0.923; AA vs. AG+GG: AOR = 0.639; 95% CI = 0.432–0.945). However, the associations were not statistically significant after adjustment for multiple tests using the FDR correction. The  miR-449b  A>G polymorphism was significantly associated with RPL risk (AA vs. GG: AOR = 2.069; 95% CI = 1.033–.4.146; AA vs AG+GG: AOR = 1.406; 95% CI = 1.011–1.955). However, the association was not significant after adjustment for multiple tests using the FDR correction. Using combination analyses ( Table 3 ), the AG/CC (AOR = 0.579; 95% CI = 0.366–0.917) combined genotype for  miR-27a / miR-423  was associated with a lower RPL risk compared with reference genotypes when variant genotypes were located in the  miR-27a  loci. In addition, CC/GG (AOR = 2.888; 95% CI = 1.116–7.470), CA/AG (AOR = 1.925; 95% CI = 1.110–3.338) for the  miR-423 / miR-449b , and AG/AG (AOR = 1.804; 95% CI = 1.067–3.052) for the  miR-449b / miR-605  were associated with an increased RPL risk compared with reference genotypes when variant genotypes were located in the loci of  miR-423 ,  miR-449b , and  miR-605 . These results were consistent with associations between a RPL risk and individual microRNA genotypes; however, each association was not significant after the FDR correction for multiple comparisons, suggesting a weak association. Haplotype-based analyses of the four microRNA polymorphisms for gene-gene interactions are shown in  S1 Table  (all possible allele combinations) and  Table 4  (allele combinations suggesting associations with RPL). Interaction models suggested by the MDR were evaluated using haplotype-based analyses. Among the models of the four polymorphic loci, the G-C-A-G haplotype (OR = 0.525; 95% CI = 0.321–0.859) was associated with a decreased RPL risk whereas three haplotypes, A-C-G-A (OR = 1.870; 95% CI = 1.178–2.968), A-A-A-G (OR = 2.429; 95% = 1.153–5.114), and G-C-G-G (OR = 3.214; 95% CI = 1.441–7.172), were associated with an increased RPL risk. Among the models of the three polymorphic loci, one haplotype, G-A-G of  miR-27a/miR-449b/miR-605  (OR = 0.625; 95% CI = 0.402–0.972), was associated with a reduced RPL risk, whereas haplotypes A-C-G of  miR-27a/miR-423/miR-449b  (OR = 1.498; 95% CI = 1.030–2.179); A-A-G, A-G-A, and G-G-G of  miR-27a/miR-449b/miR-605  (OR = 1.526; 95% CI = 1.037–2.244; OR = 1.649; 95% CI = 1.086–2.504; OR = 3.089; 95% CI = 1.401–6.809), C-G-G, and A-A-G of  miR-423/miR-449b/miR-605  (OR = 1.773; 95% CI = 1.110–2.833; OR = 2.253; 95% CI = 1.203–4.220) were associated with an increased RPL risk. Among the models of the two polymorphic loci, two haplotypes, C-G and A-A of  miR-423/miR-449b  (OR = 1.518; 95% CI = 1.134–2.031), were associated with a higher RPL risk. The G-G-G of  miR-27a/miR-449b/miR-605  and C-G and A-A of  miR-423/miR-449b  haplotypes remained significant after adjustment for multiple tests using the FDR correction. Multiple linear regression analyses of clinical variables in Korean RPL patients according to the quintiles of clinical variables are shown in  Table 5 . The AG and AG+GG genotypes of the  miR-27a  polymorphism showed significant contributions to the prediction of folate levels in RPL patients, with regression coefficients of 1.069 and 0.788, respectively. Although the differences in plasma homocysteine, folate, PLT, PT, aPTT, and NK cell number in relation to the four microRNA genotypes and haplotypes were evaluated by one-way analysis of variance and independent two sample  t -tests, no significant difference was found ( S2  and  S3  Tables).\nNote: RPL = recurrent pregnancy loss; BMI = body mass index; NK = natural killer; NA = not applicable; aPTT = activated partial thromboplastin time; values are mean ± standard deviation unless otherwise noted.\na  Platelet count of 178 control subjects;\nb  Activated partial thromboplastin time of 63 control subjects;\nc  Platelet count of 205 RPL patients;\nd  Activated partial thromboplastin time of 210 RPL patients;\ne \n P  values were calculated using the two-sided  t -test for continuous variables.\nNote: RPL = recurrent pregnancy loss; AOR = adjusted odds ratio by the age of participants; OR = odds ratio; CI = confidence interval; FDR = false discovery rate; HWE = Hardy-Weinberg equilibrium.\na  Fisher’s exact test;\nb  FDR-adjusted  P  value\nNote: RPL = recurrent pregnancy loss; AOR = adjusted odds ratio; CI = confidence interval.\na  Fisher’s exact test;\nb  FDR-adjusted  P  value.\nNote: RPL = recurrent pregnancy loss; OR = odds ratio; CI = confidence interval; ORs and 95% CIs of each haplotype combination were calculated with reference to frequencies of all others using Fisher’s exact test.\na  Fisher’s exact test;\nb  FDR-adjusted  P  value.\nNote: RPL = recurrent pregnancy loss; Hcy = homocysteine; NK cells = natural killer cells; PLT = platelet count; PT = prothrombin time; aPTT = activated partial thromboplastin time; R 2  = coefficient of determination; Coef = regression coefficients; Ref = reference.\na  Homocysteine 10 quintiles: Hcy≤4.69, 4.69<Hcy≤5.46, 5.46<Hcy≤5.9, 5.9<Hcy≤6.31, 6.31<Hcy≤6.73, 6.73<Hcy≤7.2, 7.2<Hcy≤7.55, 7.55<Hcy≤8.2, 8.2<Hcy≤9.22, Hcy>9.22;\nb  Folate 10 quintiles: folate≤5.34, 5.34<folate≤6.70, 6.70<folate≤7.79, 7.79<folate≤9.62, 9.62<folate≤11.16, 11.16<folate≤13.25, 13.25<folate≤16.59, 16.59<folate≤19.65, 19.65<folate≤22.54, folate>22.54;\nc  NK cells 10 quintiles: NK cells≤9, 9<NK cells≤11, 11<NK cells≤13, 13<NK cells≤15, 15<NK cells≤17, 17<NK cells≤19.5, 19.5<NK cells≤21.5, 21.5<NK cells≤25, 25<NK cells≤30, NK cells>30;\nd  Platelet count 10 quintiles: PLT≤187, 187<PLT≤207, 207<PLT≤224, 224<PLT≤237, 237<PLT≤251, 251<PLT≤266, 266<PLT≤279, 279<PLT≤302, 302<PLT≤329, PLT>329; prothrombin time 10 quintiles:\ne  PT≤10.5, 10.5<PT≤11, 11<PT≤11.2, 11.2<PT≤11.4, 11.4<PT≤11.6, 11.6<PT≤11.75, 11.75<PT≤11.9, 11.9<PT≤12.25, 12.25<PT≤12.6, PT>12.6;\nf  Activated partial thromboplastin time: aPTT≤26.8, 26.8<aPTT≤28.4, 28.4<aPTT≤29.6, 29.6<aPTT≤30.75, 30.75<aPTT≤31.9, 31.9<aPTT≤33.15, 33.15<aPTT≤34.4, 34.4<aPTT≤36, 36<aPTT≤37.5, aPTT>37.5\n\nThe effects of miRNA polymorphisms on pregnancy loss have been reported in a limited number of studies [ 6 – 9 ]. These studies investigated the roles of  miR-125 ,  miR-196a2 ,  miR-499  and  miR-423  polymorphisms that were supported by functional evidence showing disruption of mature microRNA production and its downstream target gene [ 6 – 9 ,  19 ]. Because more evidence in support of the functional importance of miRNAs in pregnancy has been reported, we determined whether the four miRNA SNPs associated with placental or fetal development played a role in pregnancy loss [ 14 ,  20 – 27 ]. Our results showed an association between the  miR-27a  variant G allele and a lower RPL risk, and an association between the  miR-449b  variant G allele and a higher RPL risk. The effects of variant alleles of  miR-27a  and  miR-449b  were also suggested using combination and haplotype-based analyses. Multiple linear regression analyses of clinical variables in Korean RPL patients revealed statistically significant relationships between the  miR-27a  genotypes and plasma folate levels. Functional analyses indicated that the variant genotypes of  miR-27a , AG, and GG might be responsible for the elevated  miR-27a  levels [ 14 ], and it has been speculated that the elevated  miR-27a  levels contribute to the increased folate concentration that is protective against RPL [ 20 ]. However, the possible effects of  mir-27a  on RPL mediated by folate require further investigation to confirm this hypothesis. The  miR-449b  was chosen for the study because it is one of the four miRNAs that was downregulated in hatched blastocysts [ 21 ]. However, the effect of rs10061133 A>G of  miR-449b  has rarely been reported, and there have been contradictory reports on its effects, depending on the cell type and gene expression pattern [ 24 ]. In our study, we found that the GG and AG+GG genotypes of  miR-449b  were associated with an increased risk of RPL. We, therefore, hypothesize that the GG and AG+GG genotypes of  miR-449b  rs10061133 affect the risk of RPL risk by modulating the expression of mature  miR-449b  [ 24 ,  27 ].\nThe  miR-423  rs6505162 polymorphism was chosen for the present study because it was reported to be differentially regulated in placental injury [ 25 ], and rs6505162 in pre- miR-423  affects the mature miRNA expression by increasing miRNA expression with the variant A allele [ 26 ]. However, in our study, we did not find a statistically significant association between  miR-423  genotypes and RPL risk. Notably, a recent study in the Han Chinese population also failed to identify an association between the  miR-423  genotypes and RPL risk; however, the study found an association between  miR-423  alleles and RPL, with functional evidence showing that the minor A allele contributed to an increased expression of mature  miR-423  [ 7 ]. An explanation for the different results may be related to the different allele frequencies. The  miR-605  rs2043556 A>G polymorphism was included in the study because the variant G allele of  miR-605  caused a 2.6-fold reduction in the processing levels of  miR-605  [ 23 ], and  miR-605  was significantly dysregulated in placentas after exposure to Bisphenol A, which may disrupt endocrine levels [ 22 ]. However, we did not find a significant association between the  miR-605  polymorphism and a risk of RPL. A miR-125a polymorphism involving one nucleotide mutation in the pri-miR-125a coding region related to RPL was not included in the study because the location and nomenclature of the polymorphism were not clear, although the A>G mutation reduced the expression of mature miR-125a [ 8 ].\nThere were some limitations in our study. First, this case-control study identified an association between SNPs in miRNAs in RPL development, but it did not establish a cause-and-effect relationship. What we identified, therefore, warrants functional studies to elucidate the pathogenesis related to RPL. Second, we reported a significant relationship between  miR-27a  genotypes and plasma folate levels in Korean RPL patients, and we speculated that elevated  miR-27a  levels contributed to the increased folate levels. However, further functional studies are necessary to directly investigate the effect of  mir-27a  on RPL when mediated by folate levels. Third, we did not determine the effect of miRNA polymorphisms on miRNA expression in the placenta, which might have contributed to the potential roles of miRNA during the peri-implantation and early pregnancy periods.\nAccumulation of findings on post-transcriptional regulation of gene expression by miRNAs and their involvements in trophoblast differentiation, proliferation, and angiogenesis during the developments of placenta or fetus have broaden our perspectives on the roles of miRNA and microRNA machinery gene polymorphisms during the pathogenesis of idiopathic RPL since three recent reports on the association of miRNA polymorphisms with RPL in Han-Chinese and Korean women [ 6 – 9 ,  12 ,  19 ,  28 – 32 ]. In the present study, we identified associations between miRNA polymorphisms ( miR-27a A>G and  miR-449b A>G) and RPL risk in Korean females, and identified a relationship between miRNA polymorphism ( miR-27a A>G) and plasma folate levels. However, what we identified in our study warrants further investigation to elucidate the underlying mechanism in which miRNA polymorphisms modulate folate levels and RPL development.\n\n(DOCX)\nClick here for additional data file.\n(DOCX)\nClick here for additional data file.\n(DOCX)\nClick here for additional data file.","source_license":"CC-BY-4.0","license_restricted":false}