Factor V Leiden 1691G > A mutation and the risk of recurrent pregnancy loss (RPL): systematic review and meta-analysis

Background Although numerous replication case-control studies have attempted to determine the association between Factor V Leiden (FVL) 1691G > A mutation and susceptibility to Recurrent pregnancy loss (RPL), there have been confliction among the results of various ethnic groups. To address this limitation, here we implemented first meta-analysis to provide with consistent conclusion of the association between FVL 1691G > A mutation and RPL risk. Methods After a systematic literature search, pooled odds ratio (OR) and their corresponding 95% confidence interval (CI) were used to evaluate the strength of the association. Additionally, meta-regression analyses were performed to find potential source of heterogeneity. Results In this meta-analysis, 62 studies, containing 10,410 cases and 9406 controls, were included in quantitative analysis. Overall population analysis revealed a significant positive association in the dominant (OR = 2.15, 95% CI = 1.84–2.50, P < 0.001), over-dominant (OR = 1.88, 95% CI = 1.61–2.19, P < 0.001), allelic (OR = 2.05, 95% CI = 1.79–2.35, P < 0.001), and heterozygote (OR = 1.97, 95% CI = 1.68–2.30, P < 0.001) models. Moreover, a significant association of dominant (OR = 3.04, 95% CI = 2.04–4.54, P < 0.001), over-dominant (OR = 2.65, 95% CI = 1.74–4.05, P < 0.001), and heterozygote (OR = 2.67, 95% CI = 1.81–4.22, P < 0.001) models was found in the Iranian population. The subgroup analysis indicated strong significant association in Asian, European, Africa population, and case-control studies but not in South Americans and cohort studies. Conclusion The FVL 1691G > A mutation and the risk of RPL confers a genetic contributing factor in increasing the risk of RPL, particularly in Iranians, except for South Americans.


Introduction
Recurrent pregnancy loss (RPL) is a heterogeneous disorder which affects women of reproductive age. Recently, The American Society of Reproductive Medicine has defined RPL as two or more than two failed pregnancies before the 20th week of pregnancy [1][2][3]. Overall, 1-5% of women during reproductive ages could be affected [4].
From pathophysiological point of view, RLP might be influenced by various items, such as genetic factors (chromosomal aberrations, genetic polymorphisms), infectious diseases, structural abnormalities of the uterus, coagulative disorders (thrombophilia), endocrinological problems (thyroid disease and diabetes), and immunological disease (autoimmune disorder and inflammatory diseases) [5][6][7]. With considering these factors, still approximately 40 to 50% of cases remained idiopathic [8].
Although pregnancy as a physiological condition is associated with a hypercoagulable state, and the contact between placenta and maternal circulation is crucial for the establishment of a successful pregnancy, but any abnormality in this circulation, especially abnormal blood clotting in the small placental blood vessels, may results in RPL [9,10]. During last decades, thrombophilia attracted a lot of attention as a risk factor for RLP. Thrombophilia is characterized as a hemostatic disorder which leads to an increased tendency of thromboembolic processes. Classically, thrombophilia could be classified into acquired and inherited forms [11,12]. In this regards, antiphospholipid syndrome is an established acquired thrombophilia factor which increase the risk of RPL. Among inherited factors, mutation in Factor V Leiden (FVL) of the FV gene, G20210A of the FII (prothrombin) gene, and C677T of the methylenetetrahydrofolate reductase (MTHFR) gene are believed to play a key role in pathogenesis of RPL [13,14].
FVL mutation shows an autosomal dominant pattern which occurs by substitution of guanine by adenine (CGA---> CAA) at the nucleotide 1691 in the exon 10. As a result of this missense mutation, arginine (Arg) at amino acid 506 is substituted with glutamine (Gln), leading to generation of FVL resistant to the activated protein C (APC). APC is a natural anticoagulant which in normal situation cleaves activated factor V at amino acid 506 and makes it inactive [15][16][17][18][19][20].
Studies have shown that FVL mutation increases the risk of venous thrombosis 7 times in heterozygote and 80 times in homozygote carriers. In addition, it has been reported that this mutation increases the risk of preeclampsia in FVL carriers [21,22]. The exact mechanism that FVL mutation influence the etiology of RPL is a controversial issue and has not yet been divulged thoroughly, but several studies suggested that production of micro thrombosis could sediment in delicate placental blood vessels and cause placental infarction and subsequent maternal and fetal complications [23,24].
In spite of all findings, still the exact association between FVL mutation and the risk RPL is unclear and several investigators worldwide try to clarify this question. Therefore, here we conducted the first and the most comprehensive meta-analysis on the association between FVL 1691G > A mutation and risk of RPL by exerting 62 studies encompassing 10,410 cases and 9406 health control to achieve more reliable conclusion.

Methods
Ethical approval is not necessary for this meta-analysis. The current meta-analysis was conducted according to the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) statement [25], including publication search, study selection, inclusion and exclusion criteria, data extraction, quality assessment, and statistical analysis.

Publication search
A comprehensive systematic search in the ISI Web of Science, Scopus, and PubMed/Medline databases was conducted to retrieve all publications evaluating the associations between FVL 1691G > A mutation and susceptibility to RPL prior to May 2020. The following combinations of key words were used: ("Miscarriage" OR "abortion" OR "pregnancy loss" OR "habitual abortion" OR "fetal loss" OR "Recurrent Pregnancy Loss") AND ("Factor V Leiden" OR "FV Leiden" OR "1691G > A" OR "rs6025") AND ("polymorphism*" OR "variant" OR "mutation" OR "genotype" OR "allele" OR "single nucleotide polymorphism" OR "SNP"). In spite of detailed search, a manual cross-check of eligible studies and reviews was carried out to include other potential studies. Original data in English language and human population studies were collected.

Study selection
Primary search strategy generates 1266 studies that were exported into Endnote X8 software. The duplicated studies were removed and title & abstract of remaining studies were reviewed by two investigators and irrelevant studies were excluded. Full-text verification was performed if we could not classify studies based on title & abstract. Any disagreements during study selection were discussed and resolved by consensus.

Inclusion and exclusion criteria
Studies considered eligible if they met the following inclusion criteria: a) Studies concerning the association between FVL 1691G > A mutation and susceptibility to recurrent pregnancy loss as the main outcome; b) Studies that their case group have recurrent pregnancy loss (two or more times of abortion); c) Studies with case-control and cohort design; d) Studies reporting sufficient data of genotype or allele frequency that could confer feasibility of calculating the odds ratios (ORs) and 95% confidence intervals (CIs). On the other hand, duplicates, case reports, book chapters, reviews, letter to editor, studies with insufficient data, and abstracts were all excluded.

Data extraction and quality assessment
According to a standardized extraction form, the following data were independently extracted by two investigators: the first author's last name, journal and year of publication, country of origin, ethnicity, allele and genotype frequency in cases and controls, mean or range of age, genotyping method, and total sample size of cases and controls. The third investigator finalized the extracted data, and potential discrepancies were resolved by consensus. For quality assessment of the included publications, the Newcastle-Ottawa Scale (NOS) was applied [26]. In this respect, studies with 0-3, 4-6 or 7-9 scores were of, respectively, low, moderate, and high-quality.

Statistical analysis
Deviation from Hardy-Weinberg equilibrium (HWE) for distribution of the genotype frequencies was analyzed by χ2-test in the control group. The strength of the association between FVL 1691G > A mutation and RPL risk was evaluated by the pooled OR and its corresponding 95% CI. Different comparison models for FVL 1691G > A mutation were as follow: dominant model (AA+GA vs. GG), overdominant model (GA vs. GG + AA), allelic model (A vs. G), and heterozygote (GA vs. GG). It should be noted that due to the AA genotype frequency of zero in both cases and controls, the recessive and homozygote models were not calculable. Presence of heterogeneity between included studies was estimated by Cochran's Q-statistic (P value< 0.10 was considered as statistically significant) [27]. Besides, to report quantitative heterogeneity I-squared (I 2 ) tests was used. The fixed-effected model (FEM) was used if P Q-statistic > 0.10 or I 2 was< 50%; otherwise, the random-effected model (REM) was applied. In order to assessed the predefined sources of heterogeneity among included studies, subgroup analysis and meta-regression analysis based on year of population, the continent of the study population, and genotyping method were performed. Additionally, sensitivity analysis was conducted in presence of heterogeneity [28,29]. Publication bias was estimated by Begg's funnel plots and Egger's regression test (P value< 0.05 was considered as statistically significant) [30,31]. The funnel plot asymmetry was assessed with the Egger's test. Practically, in case of no evidence of publication bias, studies with high precision (large study effects) will be located near the average line, and studies with low precision (small study effects) will be spread equally on both sides of the average line; any deviation from this shape can indicate publication bias. The data analyses were carried out using STATA (version 14.0; Stata Corporation, College Station, TX) and SPSS (version 23.0; SPSS, Inc. Chicago, IL) software.

Study characteristics
The four-phase search and screening process of the literatures based on the PRISMA statement is depicted in the Fig. 1. According to the aforementioned keywords, a    (Table 3).

Subgroup analysis by continent
The included studies were performed in Asia (25 studies), Europe (26 studies), South America (6 studies), Africa (4 studies) and Oceania (1 article). Since there was only one study for Oceania, we exclude it from the subgroup analysis. The final results revealed strong significant association between FVL 1691G > A mutation and the risk of RPL in Asian, European, and Africa population, but not in South Americans (Fig. 3).
The results of pooled ORs, heterogeneity tests, and publication bias tests in different analysis models are shown in the Table 3.  (Table 3).

Heterogeneity and publication bias
To check existence of publication bias, Egger's linear regression and Begg's funnel plot test were used. The shape of the funnel plots did not disclose obvious asymmetry under all the genotype model of the FVL 1691G > A mutation (Fig. 4). Additionally, some degree of heterogeneity was detected in overall population. Therefore, we stratified study by continent and study design to find its potential source.

Meta-regression analyses
Meta-regression analyses were performed to explore potential sources of heterogeneity among included studies ( Table 4). The findings indicated that none of the expected heterogeneity parameter were the source of heterogeneity (Fig. 5).

Sensitivity analysis
The impact of individual study on pooled OR was evaluated by sequential omission of each studies. The analysis results showed that no individual study significantly affected the pooled ORs under any genotype models of the FVL 1691G > A mutation (Fig. 6).

Discussion
RPL has been one of the most prevalent obstetric complications, that affect more than 30% of gestations. A remarkable amount of pregnancy losses has been attributed to genetic variations, of which over 50% have been related to chromosomal abnormalities. Several investigations have reported the association of FVL 1691G > A mutation with RPL; that notwithstanding, there have been conflicting results among various ethnicities. The inconsistent results have been attributed to variety in the race of included subjects, different diagnostic criteria of patients, little statistical power, small sample sizes, and the linkage disequilibrium (LD) between various genes and variations [88]. However, metaanalysis strategy provides a pertinent tool to settle the problem of confliction by resolving the limitations of single replication studies, such as limited statistical power and little sample size. Thus, here we conducted the first meta-analysis to find a valid estimation of the association between FVL 1691G > A mutation and risk of RPL. The FVL 1691G > A mutation is a G-to-A point mutation at nucleotide 1691 in the factor V gene, that results in the single amino-acid replacement Arg506Gln, leading to resistance to be cleaved and, therefore, inactivation by APC and promoted susceptibility to clotting [89,90]. This mutation enhances the risk of venous thrombosis up to 50-100 times in homozygote carriers [22].
In this meta-analysis, 62 studies, containing 10,410 cases and 9406 controls, were included in quantitative analysis. The analysis of overall population indicated that all genetic comparisons of the FVL 1691G > A mutation, including dominant model (OR = 2.15), over-dominant model (OR = 1.88), allelic model (OR = 2.05), and heterozygote model (OR = 1.97) significantly increased the risk of RPL susceptibility. In 2015, Sergi et al. [91] by including nine studies, containing a total of 2147 women for the FVL mutation, 1305 women with early RPL, and 842 women with no gestational complications, indicated higher carrier frequency of FVL mutation in women with early RPL (OR = 1.68). Moreover, Marcelo and colleagues [92] in 2019 revealed that there was no association between recurrent miscarriage and inherited thrombophilias in patients with polycystic ovarian syndrome, with respect to FVL (OR = 0.74; 95% CI = 0.38 -1.45; P = 0.38), among others. On the other hand, a comprehensive systematic review and meta-analysis in 2016 [93], by exerting 369 articles evaluating 124 polymorphisms of 73 genes, to explore the potential genetic biomarkers for recurrent miscarriage identified increased risk of the disease in the recessive and over-dominant models, but a decreased risk in the dominant and allelic models for FVL 1691G > A mutation, both in overall analysis and subgroup analysis in Caucasians. Our analysis is unique of its type, as it included only patients having RPL diagnosis. Moreover, our subgroup analysis based on the continent of the study population divulged a strong association between FVL 1691G > A mutation and the risk of RPL in Asian, European, and Africa populations, but not in South Americans. It should be noted that among the 62 case-control studies included, 25 studies were in Asia, 26 studies in Europe, 6 studies in South America, 4 studies in Africa, and 1 study in Oceania. Although the subgroup analysis of 6 studies in South America indicated an OR < 1 (which was not significant across all genetic models), all other populations (which made large portion of the studies included) had OR > 1, imply that the South America data had little effect on the pooled effect estimation. The other parameter for subgroup analysis was study design. In this regard, a significant positive association between FVL 1691G > A mutation and the risk of  On the other side, the analysis was also performed in the Iranian population, containing 9 publications with 1409 cases and 1160 controls. The previous metaanalysis in Iranian population by Kamali et al. [94] in 2018, by employing 7 studies, indicated significant increased risk of RPL only in the allelic (OR = 2.252) and dominant models (OR = 2.217). However, our analysis indicated that the measured genetic models, including dominant model (OR = 2.97), over-dominant model (OR = 2.58), and heterozygote model (OR = 2.67, 95%) increased the risk of RPL. The difference between our analysis and the previous one was that we included two more study with higher sample size.
There was a degree of heterogeneity during the overall analysis. From statistical perspective, this heterogeneity describes the variability between included studies and may originate from clinical or methodological heterogeneity, from other unreported, unknown study characteristics, or may be due to chance. Therefore, for finding any sources of heterogeneity and attenuating their effects, we conducted subgroup analysis and weighted meta-regression. Collectively, the results of metaregression showed that none of the parameters, including publication year, the continent of the study population, and genotyping methods were the expected source of heterogeneity. However, subgroup analysis reduced heterogeneity in all groups and explained part of the   observed heterogeneity expect Asians and studies with cohort design. Furthermore, the other way of dealing with statistical heterogeneity, which we used in our analysis, was to incorporate "Random" term to account for it in a random-effects. Random effect model typically produces more conservative estimates of the significance of a result (a wider confidence interval). As it gives proportionately higher weights to smaller studies and lower weights to larger studies than fixed effect analysis.
To address the limitations in the current metaanalysis, it should be stated that, first our literature search was limited to only studies published in English language. Second, there was a degree of heterogeneity during the overall analysis. But not in all subgroup analyses, indicating the role of genetic diversity and other confounders in susceptibility to RPL. Third, as this meta-analysis a crude estimation of the association between FVL 1691G > A mutation and the risk of RPL, Fig. 6 Sensitivity analysis in the present meta-analysis investigates the association of FVL 1691G > A mutation an risk of RPL; a: overall population, b: Iranian studies thus the roles of age, paternal genetic impression, environmental factors, and the effect of gene-gene interactions in conferring the susceptibility risk to RPL were neglected.
Considering all the facts, this meta-analysis, the first one of its type to our best knowledge, retrieved 62 studies, encompassing 10,410 cases and 9406 health controls, to find a consistent result of the association between FVL 1691G > A mutation and risk of RPL. Our results indicated statistically significant increased risk of RPL in the overall analysis. The increased susceptibility to RPL was also observed in Iranian, Asian, European, Africa populations, and studies with case-control design, but not in South Americans and studies with cohort design. Further experiments, alongside with inclusion of additional studies with large sample sizes, should consider the role cofounders in susceptibility to RPL.