Associations between Individual and Combined Polymorphisms of the TNF and VEGF Genes and the Embryo Implantation Rate in Patients Undergoing In Vitro Fertilization (IVF) Programs

Abstract

Background

A multiple pregnancy is now considered to be the most common adverse outcome associated with in vitro fertilization (IVF). As a consequence, the identification of women with the best chances of embryo implantation is a challenge in IVF program, in which the objective is to offer elective single-embryo transfer (eSET) without decreasing the pregnancy rate. To date, a range of hormonal and clinical parameters have been used to optimize eSET but none have significant predictive value. This variability could be due to genetic predispositions related to single-nucleotide polymorphisms (SNPs). Here, we assessed the individual and combined impacts of thirteen SNPs that reportedly influence the outcome of in vitro fertilisation (IVF) on the embryo implantation rate for patients undergoing intracytoplasmic sperm injection program (ICSI).

Materials and Methods

A 13 gene polymorphisms: FSHR(Asn680Ser), p53(Arg72Pro), AMH(Ile49Ser), ESR2(+1730G>A), ESR1(−397T>C), BMP15(−9C>G), MTHFR1(677C>T), MTHFR2(1298A>C), HLA-G(−725C>G), VEGF(+405G>C), TNFα(−308A>G), AMHR(−482A>G), PAI-1(4G/5G), multiplex PCR assay was designed to genotype women undergoing ICSI program. We analyzed the total patients population (n = 428) and a subgroup with homogeneous characteristics (n = 112).

Results

Only the VEGF(+405G>C) and TNFα(−308A>G) polymorphisms impacted fertilization, embryo implantation and pregnancy rates. Moreover, the combined VEGF+405.GG and TNFα-308.AG or AA genotype occurred significantly more frequently in women with high implantation potential. In contrast, the VEGF+405.CC and TNFα-308.GG combination was associated with a low implantation rate.

Conclusion

We identified associations between VEGF(+405G>C) and TNFα(−308A>G) polymorphisms (when considered singly or as combinations) and the embryo implantation rate. These associations may be predictive of embryo implantation and could help to define populations in which elective single-embryo transfer should be recommended (or, conversely, ruled out). However, the mechanism underlying the function of these polymorphisms in embryo implantation remains to be determined and the associations observed here must be confirmed in a larger, more heterogeneous cohort.

About the Authors:

Radia Boudjenah

* E-mail: [email protected] (RB); [email protected] (FV)

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Denise Molina-Gomes

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Antoine Torre

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Florence Boitrelle

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Stéphane Taieb

Affiliation: EA 2493, Versailles Saint Quentin University, Versailles, France

Esther Dos Santos

Affiliation: EA 2493, Versailles Saint Quentin University, Versailles, France

Robert Wainer

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Philippe de Mazancourt

Affiliation: EA 2493, Versailles Saint Quentin University, Versailles, France

Jacqueline Selva

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

François Vialard

* E-mail: [email protected] (RB); [email protected] (FV)

Affiliations Department of Reproductive Biology, Genetics, Gynaecology and Obstetrics, Poissy Saint Germain Medical Centre, Poissy, France, EA 2493, Versailles Saint Quentin University, Versailles, France

Introduction

In vitro fertilization (IVF) is a complex, multistep process. Poor control of embryo implantation is still a major obstacle to the achievement of a successful pregnancy. Only about 30% of IVF transfer cycles lead to clinical pregnancies [1]. Low implantation rates have prompted physicians to transfer multiple embryos. However, this strategy increases the likelihood of multiple pregnancy, which constitutes a risk for both the mother and the foetuses [2] and is now considered to be the most common adverse outcome associated with IVF. Indeed, multiple pregnancy increases the risk of preeclampsia, hypertension, preterm birth, preterm rupture of membranes, low birth weight, foetal death and cerebral palsy [2]. It is known that reducing the number of embryos transferred (through elective single-embryo transfer (eSET)) minimizes the incidence of multiple pregnancies. Although several strategies for selecting the right patient population (in terms of fertility criteria) for eSET have been suggested and evaluated, individual patient variability has proven to be a major drawback.

Embryo implantation is a multifactorial event that depends on interplay between the blastocyst and the receptive endometrium. Technical progress has led to improved selection of high-quality embryos, as characterized by (i) the presence of 4 to 5 blastomeres at day 2 after insemination or at least 7 to 10 blastomeres on day 3, (ii) the absence of multinucleated blastomeres, (iii) fewer than 20% of cellular fragments on days 2 and 3 after fertilization and (iv) timing of cleavage events [3]. Furthermore, many strategies have been developed for evaluating endometrial receptivity; these include the histological dating of timed endometrial biopsy tissues [4], the quantitation of endometrial cytokines in uterine flushings [5], plasma prostaglandin-E2 (PGE2) and prostaglandin-F2α levels 24 hours prior to embryo transfer [6], non-invasive ultrasound examination of the endometrium [7], transcriptomic analysis of a timed endometrial biopsy [8], and proteomic analysis of the window of implantation [9].

Separately, none of these markers has significant predictive value. However, predictive performance can be improved by combining these parameters.

Despite these advances in blastocyst and endometrium management, there is still a need to optimise individual embryo transfer protocols and thus reduce the likelihood of implantation failure or multiple pregnancy. An additional approach to this problem is based on the pharmacogenetics of embryo implantation. To date, more than 20 reports have described a variety of genetic variants associated or not with the risk of recurrent implantation failure (RIF) in patients undergoing IVF, comparing alleles frequencies between patients with RIF and control group (Table S1). These risk factors include single nucleotide polymorphisms (SNPs) in the genes for tumour suppressor protein 53 [10]–[12], follicle-stimulating hormone receptor (FSHR) [13], [14], methylenetetrahydrofolate reductase (MTHFR) [15]–[17], human leukocyte antigen-G (HLA-G) [18], [19], vascular endothelial growth factor (VEGF) [20], [21], plasminogen activator inhibitor-1 (PAI) [20], progesterone receptor (PR) [22], [23], and cyclooxygenase-2 (COX-2) [24], serotonin transporter (SERT) and serotonin receptor 1A (5-HT1A) [25]. Moreover, it was recently suggested that the implantation process depends on combinations of several genetic and environmental factors. In support of this hypothesis, an oligo-SNP model (including the p53 Arg⁷²Pro, PAI-14G/5G and VEGF+1154 A/G polymorphisms) is reportedly associated with RIF [20].

By using allele-specific polymerase chain reactions, we recently reported, for the first time, that the presence of the TNF-308.A allele was associated with higher implantation, pregnancy and multiple pregnancy rates and a lower miscarriage rate after fresh embryo transfer in patients undergoing IVF for male infertility [26].

In the present study, we used a previously described multiplex PCR [27] to evaluate the individual and combined impacts of 13 gene polymorphisms on embryo implantation rates. Variation frequency of eight of the latter (FSHR(Asn⁶⁸⁰Ser), p53(Arg72Pro), MTHFR1(677C>T), MTHFR2 (1298A>C), HLA-G (−725C>G), VEGF(+405 G>C), TNFα (−308A>G) and PAI-1 (4G/5G) have been previously reported for patients with RIF compared with control.

Our objective was to identify (after genotyping) IVF participants who could benefit from eSET and those in whom eSET should be ruled out because of the high risk of implantation failure.

Materials and Methods

Subject population

The present study analysed the associations between the participants' genotype (for 13 genes) on one hand with embryo implantation outcome (together with the ovarian response parameters reported in an earlier study) on the other [27]. Two patient populations were studied: the entire population (comprising 428 women) and a subgroup of 112 patients selected as previously described [26], [27]. The entire population group was selected based on one major criterion: women who required an intracytoplasmic sperm injection (ICSI) procedure for male fertility factors without the contribution of any female factor. Sperm were analyzed according to the World Health Organization guidelines.

More stringent inclusion criteria were used for the selected population group, including a normal karyotype, age under 38 years, serum FSH levels below 10 IU/l on Day 3 of the procedure, severe sperm alterations requiring an ICSI procedure, a unique long agonist desensitization protocol associated with recombinant FSH treatment for ovarian stimulation and Caucasian background. Embryos were cultured for 2 days according to our standard procedures. Day 2 embryo scoring was performed according to Steer's criteria (Steer et al., 1992), including morphological grade and blastomere number. At 2–3 days after oocyte retrieval, endometrium thickness was evaluated and 1–3 embryos were transferred.

The study was approved by an independent ethics committee (the Comité Consultatif de Protection des Personnes dans la Recherche Biomédicale; project reference: 01032) and was performed in accordance with the tenets of the Declaration of Helsinki. All the women provided their prior, written, informed consent to participation.

Genotyping and allele-specific PCR

Genotyping for the 13 polymorphisms was carried out as described earlier [27]. Briefly, nuclear DNA was extracted from peripheral blood leukocytes. For each of the selected SNPs, primer pairs were designed using Primer3 online software (University of Massachusetts Medical School, U.S.A.) (Table S2). Amplification with a PCR was performed for 24 cycles at an annealing temperature of 65°C [27].

Statistical analysis

Hardy-Weinberg equilibrium was assessed in a chi-squared test. Statistical analysis was performed in three steps as previously reported [27].

To summarize, for each polymorphism, a univariate analysis was performed (using the Mann–Whitney U test) for age, FSH, LH, E2, and rFSH number units, plasma E2 levels two days before oocyte retrieval, the number of mature oocytes, the number of transferred embryos and the embryo score, based on the presence or absence of the allele types. Fisher's test was used to analyze fertilization, cleavage, embryo implantation and pregnancy rates according to the presence or absence of allele types. An inter-group difference was considered to be statistically significant if two different tests gave a result with p<0.05.

A multivariate analysis was performed for all polymorphisms, and all polymorphisms found to be significant at p≤0.30 after Holm's correction were introduced as covariates into GLM models. Variables with p>0.05 and <0.10 were kept in the final models. Lastly, in order to better identify polymorphism interactions, a multilocus analysis was performed with GLM models. All polymorphisms found to be significant at p≤0.075 in the univariate analysis (without Holm correction) were introduced into the multilocus analysis; this enabled us to decrease the number of independent covariates in 1 the models. Backward selection of variables was performed as above, and variables with p>0.05 and <0.10 were kept in the final models.

Results

The PCR efficiency ranged between 91.8% and 100% (Table S3). The lowest success rate was observed for the VEGF (rs2010963; +405 G>C) genotype. The 13 polymorphisms were in Hardy-Weinberg equilibrium, with a 1% error interval. The allele frequencies in the total patient population and the selected subset were similar to those quoted on the NCBI database (http://www.ncbi.nlm.nih.gov/snp) [27].

In a univariate analysis, only two of the 13 polymorphisms (VEGF (rs2010963; +405 G>C) and TNFα (rs1800629; −308A>G)) appeared to be significantly associated with differences in the fertilization and/or embryo implantation rates. All six of the other polymorphisms previously reported in the literature as more frequent in a group of patients with recurrent embryo implantation failures (FSHR(Asn⁶⁸⁰Ser), p53(Arg⁷²Pro), MTHFR1(677C>T), MTHFR2(1298A>C), HLA-G(−725C>G), PAI-1 (4G/5G)) were not correlated with the fertilization, embryo cleavage, embryo implantation or pregnancy rates in our population (data not shown). The five other polymorphisms (AMH(Ile49Ser), AMHR(−482A>G) ESR2(+1730G>A), ESR1(−397T>C), BMP15(−9C>G), which had been selected for evaluation of their association with the ovarian response in a previous study [27]) were not associated with the embryo implantation or pregnancy rates (data not shown) in our study.

In a multivariate analysis, only the (VEGF +405 G>C and TNFα −308A>G) combination was found to be significantly associated with the embryo implantation and pregnancy rates after controlled ovarian hyperstimulation.

When comparing the various genotypes, there were no statistically significant intergroup differences in the sperm parameters.

The VEGF polymorphism (rs2010963; +405 G>C)

Given that no differences between VEGF+405.GG and VEGF+405.CG genotypes were observed in either population (Tables S4 and S5), patients with these genotype were pooled and compared with VEGF.CC patients.

Tables 1 and 2 show the data on basal hormonal patterns, ovarian stimulation treatment and response, the amount of exogenous FSH administered, the number of retrieved mature oocytes, cleaved embryo rates and the transferred embryo score as a function of the VEGF+405 polymorphism distribution in the selected subgroup (Table 1) and in the total patient population (Table 2). The VEGF polymorphism wasn't associated with any of these characteristics except the plasma E2 level on day 3, which was significantly lower for the VEGF.CC genotypes in the selected subgroup. This result was not observed in the total patient population.

[Figure omitted. See PDF.]

Table 1. VEGF alleles in the selected subgroup: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.t001

[Figure omitted. See PDF.]

Table 2. VEGF alleles in the total patient population: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.t002

In the selected subgroup, the mean fertilisation rate was found to be significantly higher (p = 0.002) in VEGF.GG+GC patients than in VEGF.CC patients (68.0% and 48.4%, respectively). This difference persisted in the total patient population (p = 0.022), with mean rates of 63.4% and 55.4%, respectively.

Similarly, VEGF.GG+GC patients in the selected subgroup had a higher mean embryo implantation rate than VEGF.CC patients (19.6% and 0%, respectively). In view of the small number of VEGF.CC patients (n = 7), this difference was not statistically significant (p = 0.058). However, a statistically significant difference (p = 0.011) was observed in the total patient population, with mean rates of 15.8% and 3.4% for the VEGF.GG+GC and VEGF.CC groups, respectively.

In the selected subgroup, VEGF.GG+GC patients had a higher pregnancy rate than VEGF.CC patients (with rates of 32.0% and 0%, respectively), although this difference was not statistically significant (p = 0.074). However, the difference was statistically significant (p = 0.033) in the total patient population, with mean rates of 24.2% and 6.9% for VEGF.GG+GC and VEGF.CC patients, respectively.

The VEGF genotype was not significantly associated with the multiple pregnancy rate in either the selected subgroup or the total patient population.

There were no significant intergroup differences in sperm characteristics (data not shown).

Considering the population as a whole, the VEGF.CC genotype's sensitivity and specificity for predicting implantation failure were respectively 8.8% and 98.2%.

The TNF polymorphism (TNFα: rs1800629; −308 A>G)

TNFα genotyping have been published previously with an allele-specific PCR [26], and similar results were obtained using our multiplex PCR. To summarize our previous data [26], in respectively the selected subgroup and the total population, a significant higher implantation (p = 0.0003 and p = 0.0146) and multiple pregnancy rates (p = 0.0037 and p = 0.0238) was observed for TNF.A allele. Pregnancy (p = 0.0239), miscarriage (p = 0.0364) and take home baby (p = 0.0014) were only impacted in the selected population.

Considering the population as a whole, the TNFα.A allele's sensitivity and specificity for predicting embryo implantation were respectively 31.6% and 78.2%.

The (TNFα −308A>G/VEGF +405 G>C) combination

Four genotype combinations were evaluated: TNF.AA+AG-VEGF.GG+GC (both of which were positive predictors), TNF.GG-VEGF.GG+GC and TNF.AA+AG-VEGF.CC (a positive and a negative predictor) and TNF.GG-VEGF.CC (both negative predictors). We did not observe any significant differences in any of the study parameters when comparing TNF.GG-VEGF.GG+GC and TNF.AA+AG-VEGF.CC (Table S6) in the entire population. Hence, the data for these genotypes were pooled. There were no patients with a TNF.AA+AG-VEGF.CC genotype in the selected subgroup and so three groups were analyzed.

The data corresponding to these three combinations (on age, basal hormonal pattern, ovarian stimulation treatment and response, the amount of exogenous FSH administered, the number of mature oocytes, cleaved embryo rates and score of transferred embryos) are shown in Table 3 (for the selected subgroup) and Table 4 (for the total patient population). There were no differences between the three groups, other than for plasma levels of E2 and FSH on day 3. In the selected subgroup, the day 3 E2 level was significantly higher (p = 0.003) for the TNF.AA+AG-VEGF.GG+GC genotype combination than for the TNF.GG-VEGF.GG+GC combination. In the total patient population, the FSH level was significantly higher (p = 0.005) for TNF.GG-VEGF.GG+GC and TNF.AA+AG-VEGF.CC than for TNF.AA+AG-VEGF.GG+GC patients.

[Figure omitted. See PDF.]

Table 3. TNFα/VEGF alleles in the selected subgroup: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.t003

[Figure omitted. See PDF.]

Table 4. TNFα/VEGF alleles in the total patient population: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.t004

In the selected subgroup, TNF.GG-VEGF.CC patients presented a lower fertilisation rate than the TNF.AA+AG-VEGF.GG+GC and TNF.GG-VEGF.GG+GC patients, with mean values of 48.4%, 64.6% and 68.8%, respectively. This was confirmed in the total patient population, with mean rates of 54.5%, 63.1% and 63.8% respectively.

After embryo transfer, patients in the selected subgroup with a TNF.AA+AG-VEGF.GG+GC genotype had a higher embryo implantation than patients with a negative predictor (the TNF.GG-VEGF.GG+GC genotype (p = 0.0003)) or two negative predictors (i.e. the TNF.GG-VEGF.CC genotype (p = 0.005)), with rates of 38.8%, 14.7% and 0%, respectively. This significant difference was confirmed in the total patient population, with rates of 21.7%, 13.7% and 1.9%, respectively.

Similarly, a higher pregnancy rate was observed for TNF.AA+AG-VEGF.GG+GC genotype in the selected subgroup, with rates of 48.0%, 26.4% and 0% for TNF.AA+AG-VEGF.GG+GC, TNF.GG-VEGF.GG+GC and TNF.GG-VEGF.CC, respectively. This difference was also observed in the total patient population, with rates of 29.2%, 22.3% and 4.2%, respectively.

These findings were consistent with a higher multiple pregnancy rate for TNF.AA+AG-VEGF.GG+GC patients (compared with TNF.GG-VEGF.GG+GC and TNF.GG-VEGF.CC patients) in the selected subgroup (with rates of 33.3%, 4.4% and 0%, respectively) and in the total patient population (with rates of 13.9%, 4.5% and 0%, respectively).

Discussion

Of the 13 polymorphisms investigated here, only two individual SNPs (VEGF+405 G/C and TNFα-308 G/A) and one SNP combination (VEGF+405/TNFα-308) appeared to be significantly associated with embryo implantation and/or pregnancy rates in an ICSI program. However, no associations were found between the other SNPs on one hand and the fertilization, embryo cleavage, embryo implantation or pregnancy rates on the other.

Patients with a VEGF.CC or TNF-308.GG genotype tended to have a lower embryo implantation rate. A combination of these genotypes was associated with a very low embryo implantation rate. To the best of our knowledge, this is the first study to have shown that the (VEGF+405/TNFα-308) combination modifies the implantation potential.

Considering the population as a whole, the TNFα.A allele had a high specificity (78.2%) for predicting embryo implantation and the VEGF.CC genotype had a high specificity (98.2%) for predicting implantation failure. However, the low sensitivities (31.6% and 8.8% for the TNFα.A allele and the VEGF.CC genotype, respectively) suggest that these markers alone cannot be used to predict embryo implantation (which is undoubtedly a multifactorial process).

The VEGF.CC genotype is associated with reduced gene expression [28]. In theory, this could lead to inadequate angiogenesis and thus lower oocyte quality, lower endometrial vascularisation and lower embryo fertilization and implantation rates (as discussed below). In our study, patients with the VEGF.CC genotype had lower oocyte fertilisation and implantation rates.

Given that pathological male factors might have reduced the fertilization, embryo quality and implantation rates, we analyzed the patients' sperm count, sperm mobility and sperm morphology as a function of genotype. No statistically significant intergroup differences were observed.

Vascular endothelial growth factor (VEGF) is a critical regulator of angiogenesis, which in turn is very important for follicular growth and selection of the dominant follicle [29], [30]. The VEGF is expressed by follicular granulosa and theca lutein cells [31]. The production of angiogenic factors by granulosa cells may help to maintain the vasculature and yield a healthy pre-ovulatory follicle [32], [33]. However, thecal rather than granulosa cells production of angiogenic factors is probably responsible for the follicular vascularisation that enables the follicle to grow and develop [30].

It has been suggested that a highly developed follicular vasculature is required to deliver hormones (and their precursors), oxygen and nutrients. The vasculature could therefore have an important role in the selection and growth of the dominant follicle [29], [30], [34], [35]. Thus, greater vascularity may be a primary determinant of follicular dominance. Conversely, low vascularity may limit the supply of nutrients, substrates and hormones to atretic follicles in vivo [30], [36]. Literature data also indicate that an active blood supply is essential for obtaining high-quality oocytes [37]. Our present data are in line with previous reports in which enhanced VEGF expression during the follicular phase was found to increase the number of predominant follicles [35], [38]–[40].

Furthermore, VEGF expression in human endometrium is depressed during the proliferative phase; it increases during the late secretory phase and reaches a peak during ovulation [41]–[43]. Given that a low embryo implantation rate could also be due to a defect in embryo invasion, the latter could lead to abnormal cytotrophoblast invasion and thus be associated with pre-eclampsia, pregnancy loss or intrauterine growth restriction. It is known that the VEGF+405.C allele is more frequent in patients with these conditions [44], [45]. Physiological changes resulting from an alteration in VEGF expression level may explain why VEGF polymorphisms are associated with differences in embryo implantation.

We have previously reported that the TNFα-308 A/G polymorphism was associated with embryo implantation rate [26]. Here, we confirmed this finding using a multiplex PCR assay and also identified an effect when the TNFα-308 A/G polymorphism is combined with the VEGF+405 G/C polymorphism.

Tumour necrosis factor and its two receptors are present in endometrium and placental trophoblasts. The factor has both inflammatory and angiogenic properties and is a potent inducer of new blood vessel growth. It also stimulates the proliferation of endometriotic stromal cells [46]. A G/A SNP in the TNFα promoter has already been linked to idiopathic miscarriage [47], [48]. However, many other researchers have failed to find significant associations in this respect [49]–[57]. In the present study, higher implantation and pregnancy rates were observed in patients with the TNFα.AG genotype, which is reportedly associated with higher TNF expression [58], [59]. Our finding is consistent with a previous report of a significant association between the TNFα concentration in endometrial secretions during the implantation window and the likelihood of a clinical pregnancy [60]. According to recent reports, the balance between pro-and anti-inflammatory cytokines is critical for implantation, placental development and pregnancy outcome. Expression of inflammation-associated T helper 1 (Th1) cytokines (such as TNF) during the pre-implantation and implantation period could favour pregnancy [61]. Moreover, TNFα is an important regulator of trophoblastic matrix metalloproteins (MMPs) (such as MMP2 [62] and MMP9 [62], [63]), which are known to mediate the invasion of an artificial extracellular matrix by trophoblastic cells [64].

Abnormal embryo–endometrial dialogue (resulting in embryo implantation failure) is reportedly related to the over- or under-expression of certain genes [65] and to the presence of polymorphic variants [66]. Our present results suggest that (i) VEGF+405 G/C influences oocyte quality (with low fertilization and implantation rates when the C allele is present) and (ii) TNF-308 influences endometrium receptivity. As discussed above, these two polymorphisms influence different steps in embryo implantation, although VEGF's association with endometrial receptivity needs to be confirmed. We also showed that a combination of polymorphisms was associated with embryo implantation. Taking VEGF.GG+GC and TNF.AA+AG as the reference for adequate implantation, the presence of a homozygous VEGF.CC or TNF.GG allele was associated with a low implantation rate and the presence of both alleles was associated with an even lower rate. The relative effect of the combination of SNPs was greater than that of a single SNP.

Our study of the association between TNF and VEGF polymorphisms and the embryo implantation rate had several limitations. Firstly, there is a lack of functional information on the manner in which polymorphisms in the TNF and VEGF genes affect the embryo implantation rate. Secondly, our study did not measure the patients' plasma VEGF and TNF levels. In the future, studies of larger, more heterogeneous cohorts will be needed to extend our understanding of the impact of TNF and VEGF on embryo implantation. If additional research identifies a definitive, causative role for TNF and VEGF polymorphisms in embryo implantation, the clinician could fine-tune his/her strategy for embryo transfer as a function of the presence of the VEGF+405 and TNFα-308 genotypes. Elective single-embryo transfer should only be performed for patients with VEGF.GG+GC and TNF.AA+AG genotypes. In contrast, eSET should not be performed for patients with VEGF.CC and TNF.GG genotypes.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Gene polymorphisms reportedly involved in embryo implantation.

https://doi.org/10.1371/journal.pone.0108287.s001

(DOC)

Table S2.

Primer design for the selected SNPs.

https://doi.org/10.1371/journal.pone.0108287.s002

(DOC)

Table S3.

PCR efficiency.

https://doi.org/10.1371/journal.pone.0108287.s003

(DOC)

Table S4.

VEGF alleles in the selected subgroup: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.s004

(DOC)

Table S5.

VEGF alleles in the total patient population: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.s005

(DOC)

Table S6.

TNFα/VEGF alleles in the total patient population: description and ART results.

https://doi.org/10.1371/journal.pone.0108287.s006

(DOC)

Acknowledgments

The authors acknowledge Christophe Fermanian for advice on the statistical analyses and David Fraser for English language revision of the manuscript.

Author Contributions

Conceived and designed the experiments: RB JS FV. Performed the experiments: RB DM ST. Analyzed the data: RB AT FV. Contributed reagents/materials/analysis tools: ED PM RW. Wrote the paper: RB AT FB FV.

Citation: Boudjenah R, Molina-Gomes D, Torre A, Boitrelle F, Taieb S, Dos Santos E, et al. (2014) Associations between Individual and Combined Polymorphisms of the TNF and VEGF Genes and the Embryo Implantation Rate in Patients Undergoing In Vitro Fertilization (IVF) Programs. PLoS ONE 9(9): e108287. https://doi.org/10.1371/journal.pone.0108287

References

1. Macklon NS, Stouffer RL, Giudice LC, Fauser BC (2006) The science behind 25 years of ovarian stimulation for in vitro fertilization. Endocr Rev 27: 170–207.

2. El-Toukhy T, Kamal A, Wharf E, Grace J, Bolton V, et al. (2009) Reduction of the multiple pregnancy rate in a preimplantation genetic diagnosis programme after introduction of single blastocyst transfer and cryopreservation of blastocysts biopsied on day 3. Hum Reprod 24: 2642–2648.

3. Meseguer M, Herrero J, Tejera A, Hilligsoe KM, Ramsing NB, et al. (2011) The use of morphokinetics as a predictor of embryo implantation. Hum Reprod 26: 2658–2671.

4. Sterzik K, Abt M, Grab D, Schneider V, Strehler E (2000) Predicting the histologic dating of an endometrial biopsy specimen with the use of Doppler ultrasonography and hormone measurements in patients undergoing spontaneous ovulatory cycles. Fertil Steril 73: 94–98.

5. Ledee-Bataille N, Lapree-Delage G, Taupin JL, Dubanchet S, Frydman R, et al. (2002) Concentration of leukaemia inhibitory factor (LIF) in uterine flushing fluid is highly predictive of embryo implantation. Hum Reprod 17: 213–218.

6. Vilella F, Ramirez L, Berlanga O, Martinez S, Alama P, et al. (2013) PGE2 and PGF2alpha concentrations in human endometrial fluid as biomarkers for embryonic implantation. J Clin Endocrinol Metab 98: 4123–4132.

7. Killick SR (2007) Ultrasound and the receptivity of the endometrium. Reprod Biomed Online 15: 63–67.

8. Horcajadas JA, Pellicer A, Simon C (2007) Wide genomic analysis of human endometrial receptivity: new times, new opportunities. Hum Reprod Update 13: 77–86.

9. Dominguez F, Pellicer A, Simon C (2009) The human embryo proteome. Reprod Sci 16: 188–190.

10. Kay C, Jeyendran RS, Coulam CB (2006) p53 tumour suppressor gene polymorphism is associated with recurrent implantation failure. Reprod Biomed Online 13: 492–496.

11. Lledo B, Turienzo A, Ortiz JA, Morales R, Ten J, et al. (2013) Negative effect of P72 polymorphism on p53 gene in IVF outcome in patients with repeated implantation failure and pregnancy loss. J Assist Reprod Genet

12. Patounakis G, Treff N, Tao X, Lonczak A, Scott RT Jr, et al. (2009) The p53 codon 72 single nucleotide polymorphism lacks a significant effect on implantation rate in fresh in vitro fertilization cycles: an analysis of 1,056 patients. Fertil Steril 92: 1290–1296.

13. Klinkert ER, te Velde ER, Weima S, van Zandvoort PM, Hanssen RG, et al. (2006) FSH receptor genotype is associated with pregnancy but not with ovarian response in IVF. Reprod Biomed Online 13: 687–695.

14. Jun JK, Yoon JS, Ku SY, Choi YM, Hwang KR, et al. (2006) Follicle-stimulating hormone receptor gene polymorphism and ovarian responses to controlled ovarian hyperstimulation for IVF-ET. J Hum Genet 51: 665–670.

15. Azem F, Many A, Ben Ami I, Yovel I, Amit A, et al. (2004) Increased rates of thrombophilia in women with repeated IVF failures. Hum Reprod 19: 368–370.

16. Haggarty P, McCallum H, McBain H, Andrews K, Duthie S, et al. (2006) Effect of B vitamins and genetics on success of in-vitro fertilisation: prospective cohort study. Lancet 367: 1513–1519.

17. Martinelli I, Taioli E, Ragni G, Levi-Setti P, Passamonti SM, et al. (2003) Embryo implantation after assisted reproductive procedures and maternal thrombophilia. Haematologica 88: 789–793.

18. Roussev RG, Coulam CB (2007) HLA-G and its role in implantation (review). J Assist Reprod Genet 24: 288–295.

19. Nardi Fda S, Slowik R, Wowk PF, da Silva JS, Gelmini GF, et al. (2012) Analysis of HLA-G polymorphisms in couples with implantation failure. Am J Reprod Immunol 68: 507–514.

20. Goodman C, Jeyendran RS, Coulam CB (2009) P53 tumor suppressor factor, plasminogen activator inhibitor, and vascular endothelial growth factor gene polymorphisms and recurrent implantation failure. Fertil Steril 92: 494–498.

21. Boudjenah R, Molina-Gomes D, Wainer R, de Mazancourt P, Selva J, et al. (2012) The vascular endothelial growth factor (VEGF) +405 G/C polymorphism and its relationship with recurrent implantation failure in women in an IVF programme with ICSI. J Assist Reprod Genet 29: 1415–1420.

22. Cramer DW, Hornstein MD, McShane P, Powers RD, Lescault PJ, et al. (2003) Human progesterone receptor polymorphisms and implantation failure during in vitro fertilization. Am J Obstet Gynecol 189: 1085–1092.

23. Coulam CB, Jeyendran RS, Roussev R (2008) Association of progesterone receptor polymorphisms with recurrent implantation failure after in vitro fertilization and embryo transfer. J Assist Reprod Genet 25: 119–122.

24. Salazar LA, Inostroza M, Jara C, Vega F, Garcia R, et al. Association of −765G>C polymorphism of the COX-2 gene with recurrent embryo implantation failure in Southern Chilean women. Clin Chim Acta 411: 1822–1824.

25. Palomares AR, Lendinez-Ramirez AM, Perez-Nevot B, Cortes-Rodriguez M, Martinez F, et al. (2013) Genetic polymorphisms of serotonin transporter and receptor 1A could influence success during embryo implantation and maintenance of pregnancy. Fertil Steril 99: 2009–e2002, 2009-2016, e2002.

26. Vialard F, El Sirkasi M, Tronchon V, Boudjenah R, Molina-Gomes D, et al. (2013) Tumor necrosis factor-308 polymorphism increases the embryo implantation rate in women undergoing in vitro fertilization. Hum Reprod 28: 2774–2783.

27. Boudjenah R, Molina-Gomes D, Torre A, Bergere M, Bailly M, et al. (2012) Genetic polymorphisms influence the ovarian response to rFSH stimulation in patients undergoing in vitro fertilization programs with ICSI. PLoS One 7: e38700.

28. Watson CJ, Webb NJ, Bottomley MJ, Brenchley PE (2000) Identification of polymorphisms within the vascular endothelial growth factor (VEGF) gene: correlation with variation in VEGF protein production. Cytokine 12: 1232–1235.

29. Geva E, Jaffe RB (2000) Role of vascular endothelial growth factor in ovarian physiology and pathology. Fertil Steril 74: 429–438.

30. Redmer DA, Reynolds LP (1996) Angiogenesis in the ovary. Rev Reprod 1: 182–192.

31. Kamat BR, Brown LF, Manseau EJ, Senger DR, Dvorak HF (1995) Expression of vascular permeability factor/vascular endothelial growth factor by human granulosa and theca lutein cells. Role in corpus luteum development. Am J Pathol 146: 157–165.

32. Redmer DA, Grazul-Bilska AT, Reynolds LP (1991) Contact-dependent intercellular communication of bovine luteal cells in culture. Endocrinology 129: 2757–2766.

33. Zheng J, Redmer DA, Reynolds LP (1993) Vascular development and heparin-binding growth factors in the bovine corpus luteum at several stages of the estrous cycle. Biol Reprod 49: 1177–1189.

34. Zeleznik AJ, Schuler HM, Reichert LE Jr (1981) Gonadotropin-binding sites in the rhesus monkey ovary: role of the vasculature in the selective distribution of human chorionic gonadotropin to the preovulatory follicle. Endocrinology 109: 356–362.

35. Wulff C, Wilson H, Wiegand SJ, Rudge JS, Fraser HM (2002) Prevention of thecal angiogenesis, antral follicular growth, and ovulation in the primate by treatment with vascular endothelial growth factor Trap R1R2. Endocrinology 143: 2797–2807.

36. Moor RM, Seamark RF (1986) Cell signaling, permeability, and microvasculatory changes during antral follicle development in mammals. J Dairy Sci 69: 927–943.

37. Fraser HM (2006) Regulation of the ovarian follicular vasculature. Reprod Biol Endocrinol 4: 18.

38. Danforth DR, Arbogast LK, Ghosh S, Dickerman A, Rofagha R, et al. (2003) Vascular endothelial growth factor stimulates preantral follicle growth in the rat ovary. Biol Reprod 68: 1736–1741.

39. Shimizu T, Jiang JY, Iijima K, Miyabayashi K, Ogawa Y, et al. (2003) Induction of follicular development by direct single injection of vascular endothelial growth factor gene fragments into the ovary of miniature gilts. Biol Reprod 69: 1388–1393.

40. Iijima K, Jiang JY, Shimizu T, Sasada H, Sato E (2005) Acceleration of follicular development by administration of vascular endothelial growth factor in cycling female rats. J Reprod Dev 51: 161–168.

41. Charnock-Jones DS, Sharkey AM, Rajput-Williams J, Burch D, Schofield JP, et al. (1993) Identification and localization of alternately spliced mRNAs for vascular endothelial growth factor in human uterus and estrogen regulation in endometrial carcinoma cell lines. Biol Reprod 48: 1120–1128.

42. Shifren JL, Tseng JF, Zaloudek CJ, Ryan IP, Meng YG, et al. (1996) Ovarian steroid regulation of vascular endothelial growth factor in the human endometrium: implications for angiogenesis during the menstrual cycle and in the pathogenesis of endometriosis. J Clin Endocrinol Metab 81: 3112–3118.

43. Gargett CE, Lederman FL, Lau TM, Taylor NH, Rogers PA (1999) Lack of correlation between vascular endothelial growth factor production and endothelial cell proliferation in the human endometrium. Hum Reprod 14: 2080–2088.

44. Banyasz I, Szabo S, Bokodi G, Vannay A, Vasarhelyi B, et al. (2006) Genetic polymorphisms of vascular endothelial growth factor in severe pre-eclampsia. Mol Hum Reprod 12: 233–236.

45. Eller AG, Branch DW, Nelson L, Porter TF, Silver RM (2011) Vascular endothelial growth factor-A gene polymorphisms in women with recurrent pregnancy loss. J Reprod Immunol 88: 48–52.

46. Iwabe T, Harada T, Tsudo T, Nagano Y, Yoshida S, et al. (2000) Tumor necrosis factor-alpha promotes proliferation of endometriotic stromal cells by inducing interleukin-8 gene and protein expression. J Clin Endocrinol Metab 85: 824–829.

47. Daher S, Shulzhenko N, Morgun A, Mattar R, Rampim GF, et al. (2003) Associations between cytokine gene polymorphisms and recurrent pregnancy loss. J Reprod Immunol 58: 69–77.

48. Costeas PA, Koumouli A, Giantsiou-Kyriakou A, Papaloizou A, Koumas L (2004) Th2/Th3 cytokine genotypes are associated with pregnancy loss. Hum Immunol 65: 135–141.

49. Quintero-Ramos A, Valdez-Velazquez LL, Hernandez G, Baltazar LM, Padilla-Gutierrez JR, et al. (2006) [Assessment of five thrombophilic genetic polymorphisms among couples with habitual abortion]. Gac Med Mex 142: 95–98.

50. Babbage SJ, Arkwright PD, Vince GS, Perrey C, Pravica V, et al. (2001) Cytokine promoter gene polymorphisms and idiopathic recurrent pregnancy loss. J Reprod Immunol 51: 21–27.

51. Baxter N, Sumiya M, Cheng S, Erlich H, Regan L, et al. (2001) Recurrent miscarriage and variant alleles of mannose binding lectin, tumour necrosis factor and lymphotoxin alpha genes. Clin Exp Immunol 126: 529–534.

52. Pietrowski D, Bettendorf H, Keck C, Burkle B, Unfried G, et al. (2004) Lack of association of TNFalpha gene polymorphisms and recurrent pregnancy loss in Caucasian women. J Reprod Immunol 61: 51–58.

53. Prigoshin N, Tambutti M, Larriba J, Gogorza S, Testa R (2004) Cytokine gene polymorphisms in recurrent pregnancy loss of unknown cause. Am J Reprod Immunol 52: 36–41.

54. Kamali-Sarvestani E, Zolghadri J, Gharesi-Fard B, Sarvari J (2005) Cytokine gene polymorphisms and susceptibility to recurrent pregnancy loss in Iranian women. J Reprod Immunol 65: 171–178.

55. Bombell S, McGuire W (2008) Tumour necrosis factor (−308A) polymorphism in pre-eclampsia: meta-analysis of 16 case-control studies. Aust N Z J Obstet Gynaecol 48: 547–551.

56. Menon R, Merialdi M, Betran AP, Dolan S, Jiang L, et al. (2006) Analysis of association between maternal tumor necrosis factor-alpha promoter polymorphism (−308), tumor necrosis factor concentration, and preterm birth. Am J Obstet Gynecol 195: 1240–1248.

57. Zammiti W, Mtiraoui N, Khairi H, Gris JC, Almawi WY, et al. (2008) Associations between tumor necrosis factor-alpha and lymphotoxin-alpha polymorphisms and idiopathic recurrent miscarriage. Reproduction 135: 397–403.

58. Liu C, Wang J, Zhou S, Wang B, Ma X (2010) Association between −238 but not −308 polymorphism of Tumor necrosis factor alpha (TNF-alpha)v and unexplained recurrent spontaneous abortion (URSA) in Chinese population. Reprod Biol Endocrinol 8: 114.

59. Zhou P, Lv GQ, Wang JZ, Li CW, Du LF, et al. (2011) The TNF-alpha-238 polymorphism and cancer risk: a meta-analysis. PLoS One 6: e22092.

60. Boomsma CM, Kavelaars A, Eijkemans MJ, Lentjes EG, Fauser BC, et al. (2009) Endometrial secretion analysis identifies a cytokine profile predictive of pregnancy in IVF. Hum Reprod 24: 1427–1435.

61. Haider S, Knofler M (2009) Human tumour necrosis factor: physiological and pathological roles in placenta and endometrium. Placenta 30: 111–123.

62. Meisser A, Chardonnens D, Campana A, Bischof P (1999) Effects of tumour necrosis factor-alpha, interleukin-1 alpha, macrophage colony stimulating factor and transforming growth factor beta on trophoblastic matrix metalloproteinases. Mol Hum Reprod 5: 252–260.

63. Cohen M, Meisser A, Haenggeli L, Bischof P (2006) Involvement of MAPK pathway in TNF-alpha-induced MMP-9 expression in human trophoblastic cells. Mol Hum Reprod 12: 225–232.

64. Staun-Ram E, Goldman S, Gabarin D, Shalev E (2004) Expression and importance of matrix metalloproteinase 2 and 9 (MMP-2 and -9) in human trophoblast invasion. Reprod Biol Endocrinol 2: 59.

65. Urman B, Yakin K, Balaban B (2005) Recurrent implantation failure in assisted reproduction: how to counsel and manage. A. General considerations and treatment options that may benefit the couple. Reprod Biomed Online 11: 371–381.

66. Madon PF, Athalye AS, Parikh FR (2005) Polymorphic variants on chromosomes probably play a significant role in infertility. Reprod Biomed Online 11: 726–732.