Next Article in Journal
STING-Triggered CNS Inflammation in Human Neurodegenerative Diseases
Previous Article in Journal
The Cortical Silent Period in the Cricothyroid Muscle as a Neurophysiologic Feature for Dystonia Observation: E-Field-Navigated Transcranial Magnetic (TMS) Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of Polymorphisms in FSHR, INHA, ESR1, and BMP15 with Recurrent Implantation Failure

1
Department of Biomedical Science, College of Life Science, CHA University, Seongnam 13488, Republic of Korea
2
Department of Obstetrics and Gynecology, Fertility Center of CHA Bundang Medical Center, CHA University, Seongnam 13520, Republic of Korea
3
Department of Obstetrics and Gynecology, Fertility Center of CHA Gangnam Medical Center, CHA University, Seoul 06125, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Biomedicines 2023, 11(5), 1374; https://doi.org/10.3390/biomedicines11051374
Submission received: 31 March 2023 / Revised: 29 April 2023 / Accepted: 3 May 2023 / Published: 5 May 2023
(This article belongs to the Section Molecular Genetics and Genetic Diseases)

Abstract

:
Recurrent implantation failure (RIF) refers to two or more unsuccessful in vitro fertilization embryo transfers in the same individual. Embryonic characteristics, immunological factors, and coagulation factors are known to be the causes of RIF. Genetic factors have also been reported to be involved in the occurrence of RIF, and some single nucleotide polymorphisms (SNPs) may contribute to RIF. We examined SNPs in FSHR, INHA, ESR1, and BMP15, which have been associated with primary ovarian failure. A cohort of 133 RIF patients and 317 healthy controls consisting of all Korean women was included. Genotyping was performed by Taq-Man genotyping assays to determine the frequency of the following polymorphisms: FSHR rs6165, INHA rs11893842 and rs35118453, ESR1 rs9340799 and rs2234693, and BMP15 rs17003221 and rs3810682. The differences in these SNPs were compared between the patient and control groups. Our results demonstrate a decreased prevalence of RIF in subjects with the FSHR rs6165 A>G polymorphism [AA vs. AG adjusted odds ratio (AOR) = 0.432; confidence interval (CI) = 0.206–0.908; p = 0.027, AA+AG vs. GG AOR = 0.434; CI = 0.213–0.885; p = 0.022]. Based on a genotype combination analysis, the GG/AA (FSHR rs6165/ESR1 rs9340799: OR = 0.250; CI = 0.072–0.874; p = 0.030) and GG-CC (FSHR rs6165/BMP15 rs3810682: OR = 0.466; CI = 0.220–0.987; p = 0.046) alleles were also associated with a decreased RIF risk. Additionally, the FSHR rs6165GG and BMP15 rs17003221TT+TC genotype combination was associated with a decreased RIF risk (OR = 0.430; CI = 0.210–0.877; p = 0.020) and increased FSH levels, as assessed by an analysis of variance. The FSHR rs6165 polymorphism and genotype combinations are significantly associated with RIF development in Korean women.

1. Introduction

Implantation is a process by which an embryo attaches to the lumen surface of the endometrium and then migrates into the deep layer of the endometrium [1]. Successful embryo implantation requires success at each event, such as sperm and oocyte quality, development of the early embryo and endometrium, and the interaction between the blastocyst and endometrium [2,3]. Recurrent implantation failure (RIF) is defined as repeated embryo implantation failure after four high-quality embryo transfers in women. Despite advances in in vitro fertilization (IVF) technology, success rates remain stable, and approximately 10% of women undergoing IVF treatment suffer from RIF [4]. Implantation failure can be a consequence of immunological, uterine, male, embryo, or coagulation factors and genetics [5,6]. Among them, embryogenesis and development are some of the most important processes in pregnancy [7]. In the IVF procedure, only embryos classified as healthy according to specific criteria are selected and transplanted into the uterus. However, a good or “transplantable” class of embryos may increase the success rate of transplantation and pregnancy [8,9]. The morphological criteria do not necessarily correlate with implantation success rates, and approximately 50% of healthy embryos fail to implant [10]. Additionally, several factors are involved in embryogenesis and development [11,12,13,14,15,16].
Maternal hormones comprise an important factor for successful implantation during pregnancy. In particular, ovarian hormones interact with signaling molecules, including cytokines and growth factors, to assist embryo implantation [17]. Additionally, hormone receptors are associated with the implantation process and pregnancy result, and changes in the hormone receptor expression and function may affect the pregnancy outcome [18,19,20].
One of the most important hormones is follicle-stimulating hormone (FSH). FSH binds to the FSH receptor (FSHR) expressed in follicular granulosa cells to induce estrogen production and secretion, which stimulates the growth and maturation of ovarian follicles and subsequently improves the quality and recovery rate of oocytes [21]. Previously, the FSHR was thought to be expressed only in the ovaries, but a recent study showed that it was expressed in extra-gonadal tissues including the cervix, endometrium, myometrium, vascular smooth muscle, and vascular endothelium [22]. One study confirmed that FSHR expression was upregulated during decidualization and in the myometrium during pregnancy. In addition, a study in mice showed that FSHR was expressed in the placental blood vessels and was related to fetal vascular formation, and deletion of the FSHR gene affected pregnancy [23,24]. Extraovarian FSHRs play an important role in establishing and maintaining successful pregnancies in humans [22]. Previously reported studies showed that FSHR was expressed in fetal vascular endothelium [23], and it was identified that the FSH–FSHR signaling system promotes the angiogenesis of vascular endothelial cells [25]. FSHR is also expressed in the uterine myometrium and plays an important role in regulating uterine muscle contraction [26]. The INHA gene encodes inhibin α, which suppresses FSHR expression, thereby inhibiting the action of FSH [27]. Inhibin α is expressed in the luminal epithelium, glandular epithelium, stromal tissue, and vascular endothelium throughout the menstrual cycle [28]. Inhibin α is secreted by mature follicles and reflects follicular maturity. One study showed that serum inhibin α is a predictor for determining oocyte maturation [29]. Furthermore, concentrations of inhibin α in follicular fluid were higher in the pregnancy group than in the non-pregnancy group, and a positive correlation was found between the number of oocytes retrieved and the fertility rate [30]. Estrogen receptors, ERα and ERβ, mediate well-characterized effects on follicle growth, maturation, oocyte release, and endometrial preparation for implantation by binding to estrogen [31]. Among them, ERα is encoded by estrogen receptor 1 (ESR1) and plays an essential role in regulating decidualization [32]. Bone morphogenetic protein (BMP) 15 is a TGFβ family member secreted by oocytes during follicle formation and is expressed in granulocytes and follicles as well as in oocytes. BMP15 plays an important role in follicular development during primordial follicle recruitment, ovulation, and corpus luteum formation and is closely related to fertilization, embryonic quality, and pregnancy outcomes. Therefore, it can be considered a new molecular marker for predicting follicular development potential [33]. In addition, one study showed that the concentration of BMP15 in women with an intracytoplasmic sperm injection was related to the fertility level, and an increased follicular concentration of BMP15 indicated an increased fertility level [34].
Single-nucleotide polymorphisms (SNPs) have been associated with reproductive diseases [35,36,37]. Additionally, RIF has previously been associated with SNPs, and many studies have been published exploring the relevance of these associations [38,39]. Polymorphisms in regulatory regions (promoter, 5′, 3′ UTR) or gene-coding regions may alter gene expression [40,41]. Maternal hormones play an important role in maintaining pregnancy, and hormone levels may be altered by certain polymorphisms [42,43]. Therefore, four genes (FSHR, INHA, ESR1, and BMP15) related to hormones were selected, and polymorphisms located in gene regulatory or coding regions were selected. Finally, a total of seven mutations were selected: coding region (FSHR rs6165, ESR1 rs2234693, and BMP15 rs17003221), promoter region (INHA rs11893842, rs35118453, and ESR1 rs9340799), and 5′UTR region (BMP15 rs3810682). The occurrence of FSHR, INHA, ESR1, and BMP15 gene polymorphisms has also been reported to be associated with various reproductive diseases, including primary ovarian insufficiency (POI), pregnancy loss, and preeclampsia [44,45,46,47]. However, few studies have examined the associations between RIF and FSHR, INHA, ESR1, and BMP15 gene polymorphisms. In this study, we investigated whether FSHR, INHA, ESR1, and BMP15 polymorphisms are associated with RIF and whether these polymorphisms affect the levels of clinical factors in Korean women. To reveal the relationship between RIF and FSHR, INHA, ESR1, and BMP15 gene polymorphisms, we assessed the differences between RIF patients and healthy controls by examining the known FSHR (rs6165), INHA (rs11893842, rs35118453), ESR1 (rs9340799, rs2234693), and BMP15 (rs17003221, rs3810682) gene polymorphisms.

2. Materials and Methods

2.1. Study Population

Blood samples were obtained from women with RIF treated at the Department of Obstetrics and Gynecology and the Fertility Center of CHA Bundang Medical Center in Seongnam, South Korea, between March 2010 and December 2022. In total, we obtained blood samples from 133 patients with RIF and 317 control participants. 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 (reference no. CHAMC2009-12-120). All embryos were examined by an embryologist prior to transfer, and embryos that showed good quality were transferred. Implantation failure refers to cases where the level of human chorionic gonadotropin measured on the 14th day of embryo transfer is <5 U/mL [48]. In the study group, subjects diagnosed with implantation failure due to anatomical, chromosomal, hormonal, infectious, autoimmune, or thrombotic causes were excluded. Uterine anatomical abnormalities in RIF patients were confirmed by hysterosalpingography, hysteroscopy, uterine sonography, computed tomography, or magnetic resonance imaging. A karyotype analysis was conducted to confirm chromosomal abnormalities, and the karyotype analysis followed the standard protocol. Hormonal causes of RIF such as hyperproactinemia, lutein insufficiency, and thyroid disease were identified by measuring prolactin, thyroid stimulating hormone (TSH), free T4, FSH, luteinizing hormone (LH), and progesterone levels in peripheral blood. Lupus anticoagulants and anti-cardiolipin antibodies were tested to confirm the autoimmune disease lupus and antiphospholipid syndrome, respectively. A deficiency of protein C and protein S and the presence of anti-β2 glycoprotein antibodies were diagnosed as thrombosis. As a control group, women with regular menstrual cycles, pregnancy history of at least one naturally conceived pregnancy, no history of pregnancy loss, and karyotype of 46, XX were recruited from CHA Bundang Medical Center.

2.2. Estimation of Homocysteine, Folate, Total Cholesterol, Uric Acid, Blood Urea Nitrogen, Creatinine, and Blood Coagulation Status

Blood samples were collected from RIF subjects after 12 h of fasting. We performed a fluorescence polarization immunoassay using the Abbott IMx analyzer (Abbott Laboratories, Abbott Park, IL, USA) to measure the homocysteine level. Folate was measured via a competitive immunoassay using the ACS 180Plus automated chemiluminescence system (Bayer Diagnostics, Tarrytown, NY, USA). Total cholesterol, uric acid, blood urea nitrogen, and creatinine were measured using commercially available enzymatic colorimetric assays (Roche Diagnostics, GmbH, Mannheim, Germany). The platelet, white blood cell, and hemoglobin levels were obtained using the Sysmex XE 2100 automated hematology system (Sysmex Corporation, Kobe, Japan). The prothrombin time (PT) and activated partial thromboplastin time (aPTT) were measured with an ACL TOP automated photo-optical coagulometer (LSI Medience, Tokyo, Japan).

2.3. Flow Cytometry Analysis of Immune Cell Proportion

Immune cell measurement was performed by flow cytometry using CellQuest software (BD FACS Calibur; BD Biosciences, Seoul, Republic of Korea). Fluorescently-labeled [fluorescein isothiocyanate, phycoerythrin (PE), peridinin chlorophyll protein, and allophycocyanin] monoclonal antibodies specific to CD3, CD4, CD8, CD19, CD16, and CD56 were purchased from BD Biosciences. Anti-NKG2A-PE antibodies were obtained from Beckman Coulter (Fullerton, CA, USA). To determine cell surface antigen expression, peripheral blood mononuclear cells (2.5 × 105) were stained for 30 min at 4 °C in the dark, washed twice with 2% phosphate-buffered saline containing 1% bovine serum albumin and 0.01% sodium azide (FACS wash buffer), and then fixed with 1% formaldehyde solution (Sigma-Aldrich, St. Louis, MO, USA) prior to sorting, as previously described [49].

2.4. Hormone Assays

On the second to third days of the women’s menstrual cycle, blood samples were collected and measured using serum samples. Following the manufacturer’s instructions, the estradiol (E2) and TSH levels were measured using radioimmunoassays (Beckman Coulter), and the FSH and LH levels were measured using enzyme immunoassays (Siemens, Munich, Germany).

2.5. SNP Selection and Genetic Analysis

We selected FSHR, INHA, ESR1, and BMP15, which are hormone-related genes associated with pregnancy. To select polymorphisms of the FSHR, INHA, ESR1, and BMP15 genes, studies on the association between pregnancy-related diseases (recurrent pregnancy loss, recurrent implantation failure, preeclampsia, premature ovarian failure, and poor ovarian response) and polymorphisms were investigated [46,50,51,52,53,54,55]. Finally, a total of seven polymorphisms in FSHR (rs6165), INHA (rs11893842 and rs35118453), ESR1 (rs9340799 and rs2234693), and BMP15 (rs17003221 and rs3810682) were selected and studied. Genomic DNA was extracted from anticoagulated peripheral blood using a G-DEX blood extraction kit (Intron, Seongnam, Republic of Korea). All genetic polymorphisms were identified by a real-time polymerase chain reaction (PCR) using the TaqMan SNP Genotyping Assay Kit (Applied Biosystems, Foster City, CA, USA). To validate the real-time analysis, DNA sequencing was performed on approximately 10~15% of the samples by random selection using an ABI 3730XL DNA Analyzer (Applied Biosystems). The concordance of the quality control samples was 100%.

2.6. Statistical Analysis

Regarding the clinical characteristics of participants, the differences between the categorical variables were analyzed using a chi-square test and continuous variables using an independent sample t-test. The data were presented as mean and standard deviations for continuous variables, and frequency and percentage for categorical variables. The statistical normality of continuous variables was confirmed using the Kolmogorov–Smirnov test. For continuous variables showing a non-normal distribution (p < 0.05 in the Kolmogorov–Smirnov test), group difference analyses were performed using the Mann–Whitney test. Allele frequencies were determined to confirm deviations from the Hardy–Weinberg equilibrium. To select the best inheritance model for a specific polymorphism, Akaike’s information criterion was calculated. The associations between FSHR, INHA, ESR1, and BMP15 polymorphisms and RIF incidence were calculated using adjusted odds ratios (AORs) and 95% confidence intervals (95% CIs) obtained from a multivariate logistic regression analysis adjusted for age. False discovery rate correction was used to adjust multiple comparison tests and provide a measure of the expected proportion of false positives among the data. The open-source MDR software package (v.2.0, www.epistasis.org accessed on 1 April 2022) was used to perform the genetic interaction analysis. Using this MDR analysis, all possible genotype combinations for gene–gene interactions were identified and analyzed. The associations between each FSHR, INHA, ESR1, and BMP15 gene polymorphism and each clinical value (platelets, PT, aPTT, homocysteine, folate, natural killer cells, uric acid, and total cholesterol) for RIF patients were assessed using ANOVA and Kruskal–Wallis tests. For the overall statistical analysis, the level of statistical significance was set as p < 0.05.

3. Results

Table 1 presents the clinical variables of the 133 RIF patients and 317 control subjects. No significant differences in age distribution were observed between the RIF and control groups, indicating that our age frequency matching was satisfactory. The blood urea nitrogen, creatinine, PT, TSH, E2, LH, BMI, total cholesterol, and white blood cell levels were different between the patient and the control group. The blood urea nitrogen (p < 0.0001), creatinine (p < 0.0001), PT (p < 0.0001), TSH (p = 0.0001), E2 (p = 0.0002), and LH (p < 0.0001) levels increased significantly in the patient group, and on the contrary, BMI (p = 0.047), total cholesterol (p < 0.0001), and white blood cell (p = 0.005) values increased significantly in the control group.
We investigated the distribution of FSHR (rs6165), INHA (rs11893842, rs35118453), ESR1 (rs9340799, rs2234693), and BMP15 (rs17003221, rs3810682) polymorphisms in RIF patients and the control group (Table 2). The AOR with respect to age was calculated from the logistic regression analysis. The frequency of each genotype in the control group was consistent with the Hardy–Weinberg equilibrium.
We investigated the genotype frequencies according to the number of implantation failures among RIF patients (Table 2). For FSHR rs6165A>G, GG homozygous genotypes and the recessive model were found to exert a protective effect against RIF [AA vs. GG: AOR, 0.463; 95% CI, 0.218–0.980; p = 0.044 and AA+AG vs. GG (recessive model): AOR, 0.430; 95% CI, 0.210–0.877; p = 0.020]. The other gene polymorphisms (INHA rs11893842 and rs35118453, ESR1 rs9340799 and rs2234693, and BMP15 rs17003221 and rs3810682) did not show significant differences between the control and patient groups. We also identified the genotype frequencies among patients according to the number of RIFs (Table 3). FSHR rs6165A>G GG genotypes and the recessive model showed protective effects in patients with more than three implantation failures (AA vs. GG: AOR, 0.412; 95% CI, 0.183–0.931; p = 0.033 and recessive model: AOR, 0.392; 95% CI, 0.156–0.982; p = 0.046), and INHA rs3511845C>T CT genotypes and CT+TT (dominant model) showed risk effects in patients with more than four implantation failures (CC vs. CT: AOR, 3.793; 95% CI,1.171–12.285; p = 0.026 and dominant model: AOR, 3.745; 95% CI, 1.176–11.921; p = 0.025).
We performed a genotype combination analysis of RIF patients and control subjects (Table 4). In the genotype combination type of FSHR rs6165 A>G and ESR1 rs9340799 A>G, the GG/AA (OR, 0.250; 95% CI, 0.072–0.874; p = 0.030) and GG/AA+AG (OR, 0.373; 95% CI, 0.171–0.816; p = 0.014) combinations were significantly associated with RIF risk. Additionally, in the FSHR rs6165A>G/INHA rs35118453C>T combination genotype, GG/CC+CT (OR, 0.399; 95% CI, 0.190–0.841; p = 0.016) showed a significantly decreased OR. In contrast, in the genotype combination analysis of INHA rs35118453 and C>T/ESR1 rs2234693 T>C, TT/TC+CC (OR, 7.001; 95% CI, 1.298–37.776; p = 0.024) was associated with an increased RIF risk. Additionally, in the INHA rs11893842 A>G and ESR1 rs9340799 A>G genotype combination, AA+AG/GG (OR, 3.065; 95% CI, 1.081–8.690; p = 0.035) was associated with an increased RIF risk. In the genotype combination analysis of ESR1 rs9340799 A>G/ESR1 rs2234693 T>C, GG/TT+TC (OR, 12.930; 95% CI, 1.492–112.052; p = 0.020) was associated with an increased RIF risk.
We performed an analysis of the differences in clinical factors according to the genotype of the FSHR (rs6165), INHA (rs11893842, rs35118453), ESR1 (rs9340799, rs2234693), and BMP15 (rs17003221, rs3810682) polymorphisms using a one-way analysis of variance (ANOVA). We found that increased LH levels were associated with the INHA rs35118453 polymorphism among all subjects (Table 5, p < 0.05). In RIF patients, increased E2 levels were associated with INHA rs11893842, and decreased CD3 (pan T cell) levels were associated with ESR1 rs2234693 (Table S1, p < 0.05). TSH levels in the control group showed an increasing trend according to the FSHR rs6165 genotype (Table S2, p < 0.05).
We then analyzed the synergistic effects of the FSHR (rs6165), INHA (rs11893842, rs35118453), ESR1 (rs9340799, rs2234693), and BMP15 (rs17003221, rs3810682) polymorphisms and clinical factors (i.e., folate, homocysteine, E2, FSH, and LH) on RIF risk (Table S3). In Table S3, we set the reference values for homocysteine, E2, FSH, and LH levels by determining the threshold for the highest 25% of levels (homocysteine 7.94 µmol/L, E2 37.6 pg/mL, FSH 9.6 U/L, and LH 5.12 U/L), and the folate level was determined by the lowest 25% of levels (8.12 ng/mL) in RIF patients and controls. FSHR rs6165 and high E2 and FSH levels had a synergistic effect on an increased susceptibility to RIF (AA+AG with high E2 levels, AOR = 11.411, p < 0.0001; AA+AG with high FSH levels, AOR = 3.009, p = 0.025, Figure 1). Additionally, INHA rs11893842, high homocysteine levels, and low folate levels had a synergistic effect on an increasing susceptibility to RIF (AG+GG with high homocysteine levels, AOR = 6.180, p = 0.004; AG+GG with low folate levels, AOR = 8.090, p = 0.003, Figure 1). Moreover, combining the ESR1 rs9340799 AG+GG type with E2 (AOR = 6.538, p = 0.0001), FSH (AOR = 5.511, p = 0.003), and LH (AOR = 12.305, p = 0.002) resulted in an increased RIF risk (Figure 2). Furthermore, the AORs of the ESR1 rs2234693 TC+CC group when combined with each female hormone factor including ≥37.6 pg/mL E2, ≥9.6 U/L FSH, and ≥5.12 U/L LH were 6.997 (p < 0.0001), 3.371 (p = 0.003), and 6.161 (p < 0.0001), respectively.

4. Discussion

The formation and quality of gametocytes and embryo development are important factors for successful pregnancy progression, and many infertile couples have experienced recurrent failure due to the poor quality or quantity of embryos before transplantation. Currently, morphological indicators are used for embryo selection, and there is no other predictive indicator that can evaluate the quality of embryo cells. Recently, many studies have elucidated molecular and genetic factors to assess embryo quality and select developmentally competent embryos. Previous studies have reported several genes associated with the pathogenesis of oocyte maturation arrest (TUBB8, PATL2) and fertilization failures (TLE6, WEE2). [56]. In this study, seven polymorphisms were selected from four hormone-related genes (FSHR rs6165, INHA rs11893842 and rs35118453, ESR1 rs9340799 and rs2234693, BMP15 rs1700321, and rs3810682), and the association between RIF occurrence and polymorphisms was analyzed.
We examined the association between RIF occurrence and SNPs in the FSHR, INHA, ESR1, and BMP15 genes. Our results showed that the genotype frequency of the FSHR rs6165 SNP was significantly different in the control and RIF patient groups. A genotype combination analysis of the seven genetic markers revealed that the GG/AA (p = 0.030) and GG/AA+AG (p = 0.014) combinations significantly reduced the RIF risk in the FSHR rs6165 and ESR1 rs9340799 combination type. Additionally, in the genotype combination analysis of FSHR rs6165A>G/INHA rs35118453C>T, GG/CC+CT (p = 0.016) showed a significantly decreased risk of RIF. In contrast, the INHA rs35118453 TT and ESR1 rs2234693 TC+CC combination were associated with increased RIF risk (p = 0.024). Moreover, the INHA rs11893842 A>G/ESR1 rs9340799 A>G, AA+AG/GG combination (p = 0.035) was associated with an increased RIF risk. In the genotype combination of ESR1 rs9340799 A>G/ESR1 rs2234693 T>C, GG/TT+TC (p = 0.020) was associated with an increased RIF risk. Furthermore, our gene–environment combinatorial analysis data revealed statistically significant relationships between the FSHR rs6165, INHA rs11893842, ESR1 rs9340799, and rs2234693 genotypes and clinical factors (folate, homocysteine, E2, FSH, and LH) in Korean RIF patients. When combined with clinical parameters, the RIF risk increased by 2- to 11-fold.
Previously, FSRH, INHA, ESR1, and BMP15 gene polymorphisms have been studied with regard to the conditions of the female reproductive system, including polycystic ovary syndrome, poor ovarian response, POI, breast cancer, and endometriosis, as well as male infertility [50,57,58,59,60,61,62,63,64,65,66,67].
The FSHR rs6165 polymorphism is located in exon 10 of the FSHR gene and is in the transmembrane region of the FSHR protein. The rs6165 polymorphism produces a change of A to G in position 919 and changes codon 307 from threonine (ACT) to alanine (GCT) [57]. According to Ganesh et al. [50], the FSHR rs6165 polymorphism was associated with unsuccessful IVF outcomes, and a higher frequency of the heterozygous AG genotype was observed in the infertile group than in the control group. For FSHR rs6165 AA carriers, the number of oocytes retrieved was significantly higher and ovarian stimulation was significantly shorter than those in GG and AG carriers [58]. Additionally, Rod et al. [59] found that FSHR rs6165 was associated with controlled ovarian stimulation and was 3-fold higher in poor responders than in good responders. Moreover, one study showed that the rs6165 AG genotype was associated with an increased risk of male infertility [60].
INHA rs11893842 and rs35119453 are located in the promoter region of the INHA gene and have been studied with regard to POI, male infertility (sperm parameters), and adrenocortical cancer [44,61,62]. Neither variant was associated with POI. Rafaqat et al. [48] found that the GG genotype frequency was increased in male infertility patients and showed a significant association with male infertility in the Pakistani population. Furthermore, rs11893842 minor alleles showed a low frequency in adrenocortical cancer patients and the rs11893842 AA genotype was associated with decreased INHA mRNA levels [62].
ESR1 rs9340799 and rs2234693 are located in intron 1, 351 and 1397 bp upstream of exon 2, respectively. The two polymorphisms are associated with cancer in females, including breast cancer and endometrial cancer [63,64], and are associated with POI [65]. The ESR1 rs9340799 GA and rs2234693 TC genotypes were associated with a decreased risk of POI, and the ESR1 rs9340799 AA and rs2234693 TT genotypes were associated with an increased risk of POI [65]. Additionally, the ESR1 rs9340799 GG genotype was associated with a 4-fold increased risk of endometriosis and a 3-fold increased risk of IVF failure in infertile patients [66]. Furthermore, the FSHR rs6165 GG, rs6165 AA, ESR1 rs9340799 GA, and rs2234693 TC gene combination enhanced the protective effect of FSHR gene variants and was associated with a reduced risk of fibrocystic mastopathy in infertile women [67].
Several clinical factors are involved in embryo implantation and pregnancy maintenance, including folate, homocysteine, and female hormones (E2, FSH, LH). Gaskins et al. [68] found that, among women receiving assisted reproductive technology, women with high serum folate levels (>26.3 ng/mL) had a 1.62-fold higher live birth rate than women with low folate levels (<16.6 ng/mL). Additionally, Ocal et al. [69] found that high homocysteine levels in follicular fluid resulted in reduced cell division and increased fragmentation in embryo cultures, which was associated with decreased oocyte and embryo quality. Another study showed that exposure to high E2 concentrations is detrimental to blastocyst implantation and early post-implantation development. Moreover, during clinically-assisted reproductive technology, high serum E2 levels not only affected the endometrium but also directly affected the blastocyst during implantation [70]. Our study revealed that patients with several clinical factors (low folate levels and high homocysteine, E2, FSH, and LH levels) and FSHR rs6165, INHA rs11893842, ESR1 rs9340799, and rs2234693 polymorphisms had an approximately 3- to 12-fold increased risk of RIF.
This study has several limitations. First, how FSHR, INHA, ESR1, and BMP15 polymorphisms influence RIF development is still unclear. The effect of these SNPs should also be confirmed through in vitro and in vivo studies. Second, additional environmental risk factors to RIFs need to be evaluated. Finally, the size of the RIF patient and control group was small, and this study group included only Koreans. To determine whether the studied genetic polymorphism can be used as a predictor of RIF, our results should be validated using larger sample sizes and other ethnic groups.

5. Conclusions

We investigated the association between the FSHR rs6165, INHA rs11893842 and rs35118453, ESR1 rs9340799 and rs2234693, and BMP15 rs17003221 and rs3810682 polymorphisms and the risk of RIF in Korean women. We found the frequency of the GG genotype of FSHR rs6165 was lower among RIF patients than among controls, suggesting this genotype may be associated with a reduced risk of RIF. Additionally, the interaction of the FSHR rs6165, INHA rs11893842, and ESR1 rs9340799 and rs2234693 polymorphisms with some clinical factors may increase the risk of RIF.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines11051374/s1, Table S1: Clinical variables in RIF patients stratified by FSHR, INHA, ESR1, and BMP15 polymorphisms status by ANOVA and Kruskal–Wallis tests; Table S2: Clinical variables in controls stratified by FSHR, INHA, ESR1, and BMP15 polymorphisms status by ANOVA and Kruskal–Wallis tests; Table S3: Synergic effect of FSHR, INHA, ESR1, and BMP15 polymorphisms with clinical risk factor.

Author Contributions

Conceptualization, E.-J.K. and J.-E.S.; methodology, J.-Y.L.; validation, C.-S.R.; investigation, E.-J.K.; resources, Y.-R.K. and E.-H.A.; data curation, E.-J.K., J.-E.S. and J.-Y.H.; writing—original draft preparation, E.-J.K.; writing—review and editing, E.-J.K. and J.-E.S.; supervision, J.-H.K. and N.-K.K.; project administration, J.-H.K. and N.-K.K.; funding acquisition, Y.-R.K., E.-H.A., J.-H.K. and N.-K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF) grants funded by the Korean Government (MSIT; grant number 2022R1F1A1064169, 2022R1F1A1074986, and 2020R1F1A1074525). This research was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health and Welfare, Republic of Korea (grant number: HR22C1605).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of the CHA Bundang Medical Center (IRB number: 2009-12-120).

Informed Consent Statement

The institutional review board of CHA Bundang Medical Center approved the study, and all patients provided written informed consent (reference no. CHAMC2009-12-120).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Achache, H.; Revel, A. Endometrial receptivity markers, the journey to successful embryo implantation. Hum. Reprod. Update 2006, 12, 731–746. [Google Scholar] [CrossRef] [PubMed]
  2. Simon, A.; Laufer, N. Assessment and treatment of repeated implantation failure (RIF). J. Assist. Reprod. Genet. 2012, 29, 1227–1239. [Google Scholar] [CrossRef]
  3. Zhou, W.; Dimitriadis, E. Secreted microRNA to predict embryo implantation outcome: From research to clinical diagnostic application. Front. Cell. Dev. Biol. 2020, 8, 586510. [Google Scholar] [CrossRef]
  4. Busnelli, A.; Reschini, M.; Cardellicchio, L.; Vegetti, W.; Somigliana, E.; Vercellini, P. How common is real repeated implantation failure? An indirect estimate of the prevalence. Reprod. Biomed. Online 2020, 40, 91–97. [Google Scholar] [CrossRef]
  5. Bashiri, A.; Halper, K.I.; Orvieto, R. Recurrent implantation failure-update overview on etiology, diagnosis, treatment and future directions. Reprod. Biol. Endocrinol. 2018, 16, 121. [Google Scholar] [CrossRef]
  6. Tan, B.K.; Vandekerckhove, P.; Kennedy, R.; Keay, S.D. Investigation and current management of recurrent IVF treatment failure in the UK. BJOG 2005, 112, 773–780. [Google Scholar] [CrossRef] [PubMed]
  7. Zhai, J.; Xiao, Z.; Wang, Y.; Wang, H. Human embryonic development: From peri-implantation to gastrulation. Trends. Cell. Biol. 2022, 32, 18–29. [Google Scholar] [CrossRef] [PubMed]
  8. Giorgetti, C.; Terriou, P.; Auquier, P.; Hans, E.; Spach, J.L.; Salzmann, J.; Roulier, R. Embryo score to predict implantation after in-vitro fertilization: Based on 957 single embryo transfers. Hum. Reprod. 1995, 10, 2427–2431. [Google Scholar] [CrossRef] [PubMed]
  9. Gardner, D.K.; Lane, M.; Stevens, J.; Schlenker, T.; Schoolcraft, W.B. Blastocyst score affects implantation and pregnancy outcome: Towards a single blastocyst transfer. Fertil. Steril. 2000, 73, 1155–1158. [Google Scholar] [CrossRef]
  10. Gardner, D.K.; Balaban, B. Assessment of human embryo development using morphological criteria in an era of time-lapse, algorithms and ‘OMICS’: Is looking good still important? Mol. Hum. Reprod. 2016, 22, 704–718. [Google Scholar] [CrossRef]
  11. Hartshorne, G.M. Factors controlling embryo viability. Hum. Fertil. 2001, 4, 225–234. [Google Scholar] [CrossRef]
  12. Assou, S.; Boumela, I.; Haouzi, D.; Anahory, T.; Dechaud, H.; De Vos, J.; Hamamah, S. Dynamic changes in gene expression during human early embryo development: From fundamental aspects to clinical applications. Hum. Reprod. Update 2011, 17, 272–290. [Google Scholar] [CrossRef]
  13. Paonessa, M.; Borini, A.; Coticchio, G. Genetic causes of preimplantation embryo developmental failure. Mol. Reprod. Dev. 2021, 88, 338–348. [Google Scholar] [CrossRef] [PubMed]
  14. Yi, H.; Xue, L.; Guo, M.X.; Ma, J.; Zeng, Y.; Wang, W.; Cai, J.Y.; Hu, H.M.; Shu, H.B.; Shi, Y.B.; et al. Gene expression atlas for human embryogenesis. FASEB J. 2010, 24, 3341–3350. [Google Scholar] [CrossRef]
  15. Wale, P.L.; Gardner, D.K. The effects of chemical and physical factors on mammalian embryo culture and their importance for the practice of assisted human reproduction. Hum. Reprod. Update 2016, 22, 2–22. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, X.; Chen, Y.; Wang, X.; Zhang, Z.; Wang, J.; Shen, Y.; Hu, Y.; Wu, X. NINJ1 triggers extravillous trophoblast cell dysfunction through blocking the STAT3 signaling pathway. Genes Genom. 2022, 44, 1385–1397. [Google Scholar] [CrossRef]
  17. Dey, S.K.; Lim, H.; Das, S.K.; Reese, J.; Paria, B.C.; Daikoku, T.; Wang, H. Molecular cues to implantation. Endocr. Rev. 2004, 25, 341–373. [Google Scholar] [CrossRef]
  18. Robinson, D.P.; Klein, S.L. Pregnancy and pregnancy-associated hormones alter immune responses and disease pathogenesis. Horm. Behav. 2012, 62, 263–271. [Google Scholar] [CrossRef]
  19. Large, M.J.; DeMayo, F.J. The regulation of embryo implantation and endometrial decidualization by progesterone receptor signaling. Mol. Cell. Endocrinol. 2012, 358, 155–165. [Google Scholar] [CrossRef]
  20. Klonos, E.; Katopodis, P.; Karteris, E.; Papanikolaou, E.; Tarlatzis, B.; Pados, G. Endometrial changes in estrogen and progesterone receptor expression during implantation in an oocyte donation program. Exp. Ther. Med. 2020, 20, 178. [Google Scholar] [CrossRef] [PubMed]
  21. Park, S.R.; Kim, S.R.; Kim, S.K.; Park, J.R.; Hong, I.S. A novel role of follicle-stimulating hormone (FSH) in various regeneration-related functions of endometrial stem cells. Exp. Mol. Med. 2022, 54, 1524–1535. [Google Scholar] [CrossRef]
  22. Stilley, J.A.W.; Segaloff, D.L. FSH actions and pregnancy: Looking beyond ovarian FSH receptors. Endocrinology 2018, 159, 4033–4042. [Google Scholar] [CrossRef] [PubMed]
  23. Stilley, J.A.W.; Christensen, D.E.; Dahlem, K.B.; Guan, R.; Santillan, D.A.; England, S.K.; Al-Hendy, A.; Kirby, P.A.; Segaloff, D.L. FSH receptor (FSHR) expression in human extragonadal reproductive tissues and the developing placenta, and the impact of its deletion on pregnancy in mice. Biol. Reprod. 2014, 91, 74. [Google Scholar] [CrossRef]
  24. Stilley, J.A.W.; Segaloff, D.L. Deletion of fetoplacental Fshr inhibits fetal vessel angiogenesis in the mouse placenta. Mol. Cell. Endocrinol. 2018, 476, 79–83. [Google Scholar] [CrossRef]
  25. Stelmaszewska, J.; Chrusciel, M.; Doroszko, M.; Akerfelt, M.; Ponikwicka-Tyszko, D.; Nees, M.; Frentsch, M.; Li, X.; Kero, J.; Huhtaniemi, I.; et al. Revisiting the expression and function of follicle-stimulation hormone receptor in human umbilical vein endothelial cells. Sci. Rep. 2016, 6, 37095. [Google Scholar] [CrossRef]
  26. Stilley, J.A.; Guan, R.; Santillan, D.A.; Mitchell, B.F.; Lamping, K.G.; Segaloff, D.L. Differential regulation of human and mouse myometrial contractile activity by FSH as a function of FSH receptor density. Biol. Reprod. 2016, 95, 36. [Google Scholar] [CrossRef]
  27. Lu, C.; Yang, W.; Chen, M.; Liu, T.; Yang, J.; Tan, P.; Li, L.; Hu, X.; Fan, C.; Hu, Z.; et al. Inhibin A inhibits follicle-stimulating hormone (FSH) action by suppressing its receptor expression in cultured rat granulosa cells. Mol. Cell. Endocrinol. 2009, 298, 48–56. [Google Scholar] [CrossRef]
  28. Leung, P.H.; Salamonsen, L.A.; Findlay, J.K. Immunolocalization of inhibin and activin subunits in human endometrium across the menstrual cycle. Hum. Reprod. 1998, 13, 3469–3477. [Google Scholar] [CrossRef] [PubMed]
  29. Lawrenz, B.; Depret Bixio, L.; Coughlan, C.; Andersen, C.Y.; Melado, L.; Kalra, B.; Savjani, G.; Fatemi, H.M.; Kumar, A. Inhibin A—A promising predictive parameter for determination of final oocyte maturation in ovarian stimulation for IVF/ICSI. Front. Endocrinol. 2020, 11, 307. [Google Scholar] [CrossRef]
  30. Ocal, P.; Aydin, S.; Cepni, I.; Idil, S.; Idil, M.; Uzun, H.; Benian, A. Follicular fluid concentrations of vascular endothelial growth factor, inhibin A and inhibin B in IVF cycles: Are they markers for ovarian response and pregnancy outcome? Eur. J. Obstet. Gynecol. Reprod. Biol. 2004, 115, 194–199. [Google Scholar] [CrossRef] [PubMed]
  31. Goldenberg, R.L.; Vaitukaitis, J.L.; Ross, G.T. Estrogen and follicle stimulation hormone interactions on follicle growth in rats. Endocrinology 1972, 90, 1492–1498. [Google Scholar] [CrossRef] [PubMed]
  32. Pawar, S.; Laws, M.J.; Bagchi, I.C.; Bagchi, M.K. Uterine epithelial estrogen receptor-alpha controls decidualization via a paracrine mechanism. Mol. Endocrinol. 2015, 29, 1362–1374. [Google Scholar] [CrossRef]
  33. Li, Y.; Li, R.Q.; Ou, S.B.; Zhang, N.F.; Ren, L.; Wei, L.N.; Zhang, Q.X.; Yang, D.Z. Increased GDF9 and BMP15 mRNA levels in cumulus granulosa cells correlate with oocyte maturation, fertilization, and embryo quality in humans. Reprod. Biol. Endocrinol. 2014, 12, 81. [Google Scholar] [CrossRef]
  34. Wu, Y.T.; Tang, L.; Cai, J.; Lu, X.E.; Xu, J.; Zhu, X.M.; Luo, Q.; Huang, H.F. High bone morphogenetic protein-15 level in follicular fluid is associated with high quality oocyte and subsequent embryonic development. Hum. Reprod. 2007, 22, 1526–1531. [Google Scholar] [CrossRef]
  35. Barbitoff, Y.A.; Tsarev, A.A.; Vashukova, E.S.; Maksiutenko, E.M.; Kovalenko, L.V.; Belotserkovtseva, L.D.; Glotov, A.S. A data-driven review of the genetic factors of pregnancy complications. Int. J. Mol. Sci. 2020, 21, 3384. [Google Scholar] [CrossRef] [PubMed]
  36. Ryu, C.S.; Sakong, J.H.; Ahn, E.H.; Kim, J.O.; Ko, D.; Kim, J.H.; Lee, W.S.; Kim, N.K. Association study of the three functional polymorphisms (TAS2R46G>A, OR4C16G>A, and OR4X1A>T) with recurrent pregnancy loss. Genes. Genom. 2019, 41, 61–70. [Google Scholar] [CrossRef]
  37. Kwon, B.N.; Lee, N.R.; Kim, H.J.; Kang, Y.D.; Kim, J.S.; Park, J.W.; Jin, H.J. Folate metabolizing gene polymorphisms and genetic vulnerability to preterm birth in Korean women. Genes Genom. 2021, 43, 937–945. [Google Scholar] [CrossRef]
  38. Lee, J.Y.; Ahn, E.H.; Park, H.W.; Kim, J.H.; Kim, Y.R.; Lee, W.S.; Kim, N.K. Association between HOX transcript antisense RNA single-nucleotide variants and recurrent implantation failure. Int. J. Mol. Sci. 2021, 22, 3021. [Google Scholar] [CrossRef] [PubMed]
  39. Salazar Garcia, M.D.; Sung, N.; Mullenix, T.M.; Dambaeva, S.; Beaman, K.; Gilman-Sachs, A.; Kwak-Kim, J. Plasminogen activator inhibitor-1 4G/5G polymorphism is associated with reproductive failure: Metabolic, hormonal, and immune profiles. Am. J. Reprod. Immunol. 2016, 76, 70–81. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, Z.; Miteva, M.A.; Wang, L.; Alexov, E. Analyzing effects of naturally occurring missense mutations. Comput. Math Methods Med. 2012, 2012, 805827. [Google Scholar] [CrossRef]
  41. Guo, Y.; Jamison, D.C. The distribution of SNPs in human gene regulatory regions. BMC Genom. 2005, 6, 140. [Google Scholar] [CrossRef] [PubMed]
  42. Borgbo, T.; Jeppesen, J.V.; Lindgren, I.; Lundberg Giwercman, Y.; Hansen, L.L.; Yding Andersen, C. Effect of the FSH receptor single nucleotide polymorphisms (FSHR 307/680) on the follicular fluid hormone profile and the granulosa cell gene expression in human small antral follicles. Mol. Hum. Reprod. 2015, 21, 255–261. [Google Scholar] [CrossRef] [PubMed]
  43. Eldafira, E.; Prasasty, V.D.; Abinawanto, A.; Syahfirdi, L.; Pujianto, D.A. Polymorphisms of estrogen receptor-alpha and estrogen receptor-beta genes and its expression in endometriosis. Turk. J. Pharm. Sci. 2021, 18, 91–95. [Google Scholar] [CrossRef] [PubMed]
  44. Pu, D.; Xing, Y.; Gao, Y.; Gu, L.; Wu, J. Gene variation and premature ovarian failure: A meta-analysis. Eur. J. Obstet. Gynecol. Reprod. Biol. 2014, 182, 226–237. [Google Scholar] [CrossRef]
  45. Meireles, A.J.C.; Bilibio, J.P.; Lorenzzoni, P.L.; Conto, E.; Nascimento, F.C.D.; Cunha-Filho, J.S.D. Association of FSHR, LH, LHR, BMP15, GDF9, AMH, and AMHR polymorphisms with poor ovarian response in patients undergoing in vitro fertilization. JBRA Assist. Reprod. 2021, 25, 439–446. [Google Scholar] [CrossRef]
  46. Jalilvand, A.; Yari, K.; Heydarpour, F. Role of polymorphisms on the recurrent pregnancy loss: A systematic review, meta-analysis and bioinformatic analysis. Gene 2022, 844, 146804. [Google Scholar] [CrossRef] [PubMed]
  47. Ke, Y.; Bin, L.; Lin, L.; MingRong, X. ESR1 polymorphisms and risk of preeclampsia. J. Matern. Fetal. Neonatal. Med. 2022, 35, 402–409. [Google Scholar] [CrossRef]
  48. Choi, Y.; Kim, J.O.; Shim, S.H.; Lee, Y.; Kim, J.H.; Jeon, Y.J.; Ko, J.J.; Lee, W.S.; Kim, N.K. Genetic variation of methylenetetrahydrofolate reductase (MTHFR) and thymidylate synthase (TS) genes is associated with idiopathic recurrent implantation failure. PLoS ONE 2016, 11, e01608842016. [Google Scholar] [CrossRef]
  49. Kim, E.S.; Kim, J.O.; An, H.J.; Sakong, J.H.; Lee, H.A.; Kim, J.H.; Ahn, E.H.; Kim, Y.R.; Lee, W.S.; Kim, N.K. MTHFR 3′-untranslated region polymorphisms contribute to recurrent pregnancy loss risk and alterations in peripheral natural killer cell proportions. Clin. Exp. Reprod. Med. 2017, 44, 152–158. [Google Scholar] [CrossRef]
  50. Ganesh, V.; Venkatesan, V.; Koshy, T.; Reddy, S.N.; Muthumuthiah, S.; Paul, S.F.D. Association of estrogen, progesterone and follicle stimulating hormone receptor polymorphisms with in vitro fertilization outcomes. Syst. Biol. Reprod. Med. 2018, 64, 260–265. [Google Scholar] [CrossRef]
  51. Santos, M.; Cordts, E.B.; Peluso, C.; Dornas, M.; Neto, F.H.V.; Bianco, B.; Barbosa, C.P.; Christofolini, D.M. Association of BMP15 and GDF9 variants to premature ovarian insufficiency. J. Assist. Reprod. Genet. 2019, 36, 2163–2169. [Google Scholar] [CrossRef]
  52. Rah, H.; Jeon, Y.J.; Ko, J.J.; Kim, J.H.; Kim, Y.R.; Cha, S.H.; Choi, Y.; Lee, W.S.; Kim, N.K. Association of inhibin alpha gene promoter polymorphisms with risk of idiopathic primary ovarian insufficiency in Korean women. Maturitas 2014, 77, 163–167. [Google Scholar] [CrossRef]
  53. Cristina Dos Santos Lopes, A.; Perucci, L.O.; Gontijo Evangelista, F.C.; Godoi, L.C.; de Paula Sabino, A.; Gomes, K.B.; Talvani, A.; Dusse, L.M.S.; Alpoim, P.N. Association among ACE, ESR1 polymorphisms and preeclampsia in Brazilian pregnant women. Mol. Cell. Probes. 2019, 45, 43–47. [Google Scholar] [CrossRef]
  54. Neves, A.R.; Vuong, N.L.; Blockeel, C.; Garcia, S.; Alviggi, C.; Spits, C.; Ma, P.Q.M.; Ho, M.T.; Tournaye, H.; Polyzos, N.P. The effect of polymorphisms in FSHR gene on late follicular phase progesterone and estradiol serum levels in predicted normoresponders. Hum. Reprod. 2022, 37, 2646–2654. [Google Scholar] [CrossRef]
  55. Sindiani, A.M.; Batiha, O.; Al-Zoubi, E.; Khadrawi, S.; Alsoukhni, G.; Alkofahi, A.; Alahmad, N.A.; Shaaban, S.; Alshdaifat, E.; Abu-Halima, M. Association of single-nucleotide polymorphisms in the ESR2 and FSHR genes with poor ovarian response in infertile Jordanian women. Clin. Exp. Reprod. Med. 2021, 48, 69–79. [Google Scholar] [CrossRef]
  56. Solovova, O.A.; Chernykh, V.B. Genetics of oocyte maturation defects and early embryo development arrest. Genes 2022, 13, 1920. [Google Scholar] [CrossRef] [PubMed]
  57. Aittomaki, K.; Lucena, J.L.; Pakarinen, P.; Sistonen, P.; Tapanainen, J.; Gromoll, J.; Kaskikari, R.; Sankila, E.M.; Lehvaslaiho, H.; Engel, A.R.; et al. Mutation in the follicle-stimulating hormone receptor gene causes hereditary hypergonadotropic ovarian failure. Cell 1995, 82, 959–968. [Google Scholar] [CrossRef]
  58. Alviggi, C.; Conforti, A.; Santi, D.; Esteves, S.C.; Andersen, C.Y.; Humaidan, P.; Chiodini, P.; De Placido, G.; Simoni, M. Clinical relevance of genetic variants of gonadotrophins and their receptors in controlled ovarian stimulation: A systematic review and meta-analysis. Hum. Reprod. Update 2018, 24, 599–614. [Google Scholar] [CrossRef] [PubMed]
  59. Rod, A.; Jarzabek, K.; Wolczynski, S.; Benhaim, A.; Reznik, Y.; Denoual-Ziad, C.; Herlicoviez, M.; Mittre, H. ESR1 and FSHR gene polymorphisms influence ovarian response to FSH in poor responder women with normal FSH levels. Endocrinol. Metab. Synd. 2014, 3, 1–5. [Google Scholar]
  60. Wu, Q.; Zhang, J.; Zhu, P.; Jiang, W.; Liu, S.; Ni, M.; Zhang, M.; Li, W.; Zhou, Q.; Cui, Y.; et al. The susceptibility of FSHB -211G > T and FSHR G-29A, 919A > G, 2039A > G polymorphisms to men infertility: An association study and meta-analysis. BMC Med. Genet. 2017, 18, 81. [Google Scholar] [CrossRef]
  61. Rafaqat, W.; Kayani, M.R.; Fatima, T.; Shaharyar, S.; Khan, S.; Ashraf, M.; Afzal, U.; Rehman, R. Association of polymorphism c.-124G>A and c.-16 C>T in the promoter region of human INHA gene with altered sperm parameters; a pilot study. Int. J. Clin. Pract. 2020, 74, e135952020. [Google Scholar] [CrossRef] [PubMed]
  62. Hofland, J.; Steenbergen, J.; Voorsluijs, J.M.; Verbiest, M.M.; de Krijger, R.R.; Hofland, L.J.; de Herder, W.W.; Uitterlinden, A.G.; Feelders, R.A.; de Jong, F.H. Inhibin alpha-subunit (INHA) expression in adrenocortical cancer is linked to genetic and epigenetic INHA promoter variation. PLoS ONE 2014, 9, e1049442014. [Google Scholar] [CrossRef]
  63. Tan, S.C.; Low, T.Y.; Mohamad Hanif, E.A.; Sharzehan, M.A.K.; Kord-Varkaneh, H.; Islam, M.A. The rs9340799 polymorphism of the estrogen receptor alpha (ESR1) gene and its association with breast cancer susceptibility. Sci. Rep. 2021, 11, 18619. [Google Scholar] [CrossRef]
  64. Ashton, K.A.; Proietto, A.; Otton, G.; Symonds, I.; McEvoy, M.; Attia, J.; Gilbert, M.; Hamann, U.; Scott, R.J. Estrogen receptor polymorphisms and the risk of endometrial cancer. BJOG 2009, 116, 1053–1061. [Google Scholar] [CrossRef]
  65. Sadat Eshaghi, F.; Dehghan Tezerjani, M.; Ghasemi, N.; Dehghani, M. Association study of ESR1 rs9340799, rs2234693, and MMP2 rs243865 variants in Iranian women with premature ovarian insufficiency: A case-control study. Int. J. Reprod. Biomed. 2022, 20, 841–850. [Google Scholar] [CrossRef]
  66. Paskulin, D.D.; Cunha-Filho, J.S.; Paskulin, L.D.; Souza, C.A.; Ashton-Prolla, P. ESR1 rs9340799 is associated with endometriosis-related infertility and in vitro fertilization failure. Dis. Markers 2013, 35, 907–913. [Google Scholar] [CrossRef]
  67. Kornatska, A.G.; Rossokha, Z.I.; Fishchuk, L.Y.; Dubenko, O.D.; Medvedieva, N.L.; Flaksemberg, M.A.; Chubei, G.V.; Popova, O.F.; Gorovenko, N.G. ESR1 gene variants affect FSHR-depended risk of fibrocystic mastopathy in infertile women. Exp. Oncol. 2021, 43, 266–269. [Google Scholar] [PubMed]
  68. Gaskins, A.J.; Chiu, Y.H.; Williams, P.L.; Ford, J.B.; Toth, T.L.; Hauser, R.; Chavarro, J.E.; Team, E.S. Association between serum folate and vitamin B-12 and outcomes of assisted reproductive technologies. Am. J. Clin. Nutr. 2015, 102, 943–950. [Google Scholar] [CrossRef] [PubMed]
  69. Ocal, P.; Ersoylu, B.; Cepni, I.; Guralp, O.; Atakul, N.; Irez, T.; Idil, M. The association between homocysteine in the follicular fluid with embryo quality and pregnancy rate in assisted reproductive techniques. J. Assist. Reprod. Genet. 2012, 29, 299–304. [Google Scholar] [CrossRef]
  70. Chang, K.T.; Su, Y.T.; Tsai, Y.R.; Lan, K.C.; Hsuuw, Y.D.; Kang, H.Y.; Chan, W.H.; Huang, F.J. High levels estradiol affect blastocyst implantation and post-implantation development directly in mice. Biomed. J. 2022, 45, 179–189. [Google Scholar] [CrossRef]
Figure 1. Synergistic effect of FSHR rs6165 and INHA rs11893842 polymorphisms with clinical parameters. (AD) panels show AOR of FSHR rs6165 and INHA rs11893842 with clinical parameters including E2 (A), FSH (B), homocysteine (C), and folate (D).
Figure 1. Synergistic effect of FSHR rs6165 and INHA rs11893842 polymorphisms with clinical parameters. (AD) panels show AOR of FSHR rs6165 and INHA rs11893842 with clinical parameters including E2 (A), FSH (B), homocysteine (C), and folate (D).
Biomedicines 11 01374 g001
Figure 2. Synergistic effect of ESR1 rs9340799 and rs2234693 polymorphisms with clinical parameters. (AF) panels show AOR of ESR1 rs9340799 and rs2234693 with clinical parameters including E2 (A,D), FSH (B,E), and LH (C,F).
Figure 2. Synergistic effect of ESR1 rs9340799 and rs2234693 polymorphisms with clinical parameters. (AF) panels show AOR of ESR1 rs9340799 and rs2234693 with clinical parameters including E2 (A,D), FSH (B,E), and LH (C,F).
Biomedicines 11 01374 g002
Table 1. Clinical profiles between RIF patient and control subjects.
Table 1. Clinical profiles between RIF patient and control subjects.
CharacteristicsControls
(n = 317)
RIF
(n = 133)
pa
Age (years, mean ± SD)33.5 ± 3.433.6 ± 2.90.579
BMI (kg/m2, mean ± SD)22.0 ± 3.421.1 ± 3.20.047
Previous implantation failure (n, mean ± SD)-4.7 ± 2.0N/A
Live births (n, mean ± SD)1.6 ± 0.6-N/A
Mean gestational age (weeks, mean ± SD)39.3 ± 1.6-N/A
Homocysteine (μmol/L mean ± SD)6.4 ± 3.06.9 ± 1.80.402
Folate (ng/mL, mean ± SD)14.0 ± 7.515.3 ± 10.70.896
BUN (mg/dL, mean ± SD)8.8 ± 2.810.4 ± 2.9<0.0001
Creatinine (mg/dL, mean ± SD)0.6 ± 0.20.8 ± 0.1<0.0001
Uric acid (mg/dL, mean ± SD3.9 ± 1.04.0 ± 1.00.266
Total cholesterol (mg/dL)215.3 ± 57.2185.2 ± 42.6<0.0001
WBC(103/μL, mean ± SD7.9 ± 2.27.4 ± 3.00.005
Hgb (g/dL, mean ± SD)12.3 ± 1.212.4 ± 1.50.080
PLT (103/μL)230.0 ± 63.1238.7 ± 67.00.465
PT (sec)10.7 ± 1.611.3 ± 0.6<0.0001
aPTT (sec)29.1 ± 3.529.6 ± 3.40.164
CD3 (pan T) (%, mean ± SD)-66.9 ± 11.2N/A
CD4 (helper T) (%, mean ± SD)-34.4 ± 8.8N/A
CD8 (suppressor) (%, mean ± SD)-28.9 ± 7.7N/A
CD19 (B cell) (%, mean ± SD)-11.7 ± 4.7N/A
CD56 (NK cell) (%, mean ± SD)-18.1 ± 9.4N/A
TSH (mU/L, mean ± SD)1.6 ± 0.92.3 ± 1.50.0001
E2 (pg/mL, mean ± SD)26.6 ± 14.433.3 ± 15.40.0002
FSH (U/L, mean ± SD)8.1 ± 2.88.9 ± 4.50.599
LH (U/L, mean ± SD)3.7 ± 2.64.8 ± 2.1<0.0001
Note: RIF, recurrent implantation failure; BMI, body mass index; BUN, blood urea nitrogen; WBC, white blood cell; Hgb, hemoglobin; PLT, platelet count; PT, prothrombin time; aPTT, activated partial thromboplastin time; TSH, thyroid stimulating hormone; E2, estradiol; FSH, follicle stimulating hormone; LH, luteinizing hormone; N/A, not applicable. a Mann-Whitney test.
Table 2. Comparison of genotype frequencies of FSHR, INHA, ESR1, and BMP15 polymorphisms between the RIF and control subjects.
Table 2. Comparison of genotype frequencies of FSHR, INHA, ESR1, and BMP15 polymorphisms between the RIF and control subjects.
GenotypesControls
(n = 317)
RIF
(n = 133)
COR (95% CI)pFDR-pAICAOR (95% CI)pFDR-pAIC
FSHR rs6165 A>G
AA134 (42.3)62 (46.6)1.000 (reference) 1.000 (reference)
AG133 (42.0)61 (45.9)0.991 (0.647–1.520)0.9680.968430.360.996 (0.650–1.528)0.9860.986428.34
GG50 (15.8)10 (7.5)0.432 (0.206–0.908)0.0270.135282.930.463 (0.218–0.980)0.0440.165 280.03
Dominant (AA vs. AG+GG) 0.839 (0.558–1.260)0.3970.702542.200.838 (0.558–1.260)0.3960.707 540.16
Recessive (AA+AG vs. GG) 0.434 (0.213–0.885)0.0220.11540.730.430 (0.210–0.877)0.0200.100 538.69
HWE-P0.0830.340
INHA rs11893842 A>G
AA106 (33.4)42 (31.6)1.000 (reference) 1.000 (reference)
AG156 (49.2)61 (45.9)0.987 (0.621–1.570)0.9560.968518.960.988 (0.621–1.573)0.9610.986518.96
GG55 (17.4)30 (22.6)1.377 (0.778–2.436)0.2720.334363.051.358 (0.767–2.407)0.2940.350 362.07
Dominant (AA vs. AG+GG) 1.089 (0.705–1.680)0.7020.702540.031.086 (0.703–1.677)0.7100.710 539.99
Recessive (AA+AG vs. GG) 1.388 (0.842–2.287)0.1990.249536.911.386 (0.841–2.286)0.2000.250 536.87
HWE-P0.8530.386
INHA rs35118453 C>T
CC209 (65.9)84 (63.2)1.000 (reference) 1.000 (reference)
CT102 (32.2)42 (31.6)1.025 (0.660–1.590)0.9140.968480.211.025 (0.660–1.592)0.9130.986479.51
TT6 (1.9)7 (5.3)2.903 (0.948–8.892)0.0620.155292.702.865 (0.934–8.786)0.0660.165 286.68
Dominant (CC vs. CT+TT) 1.129 (0.740–1.722)0.5740.702539.631.126 (0.738–1.718)0.5830.710 539.58
Recessive (CC+CT vs. TT) 2.880 (0.949–8.737)0.0620.155534.272.875 (0.947–8.724)0.0620.155 534.09
HWE-P0.1060.564
ESR1 rs9340799 A>G
AA212 (66.9)83 (62.4)1.000 (reference) 1.000 (reference)
AG94 (29.7)41 (30.8)1.114 (0.713–1.740)0.6350.968510.371.117 (0.715–1.745)0.6270.986510.09
GG11 (3.5)9 (6.8)2.090 (0.836–5.227)0.1150.192247.082.107 (0.841–5.278)0.1120.187 240.74
Dominant (AA vs. AG+GG) 1.216 (0.798–1.854)0.3630.702539.521.218 (0.799–1.857)0.3600.710 539.46
Recessive (AA+AG vs. GG) 2.019 (0.817–4.993)0.1280.213 538.122.023 (0.818–5.003)0.1270.212 538.05
HWE-P0.8840.217
ESR1 rs2234693 T>C
TT123 (38.8)44 (33.1)1.000 (reference) 1.000 (reference)
TC144 (45.4)65 (48.9)1.262 (0.803–1.983)0.3130.968445.721.258 (0.800–1.978)0.3200.986445.09
CC50 (15.8)24 (18.0)1.342 (0.739–2.436)0.3340.334279.861.330 (0.732–2.417)0.3500.350 279.60
Dominant (TT vs. TC+CC) 1.283 (0.838–1.964)0.2520.702539.021.281 (0.836–1.962)0.2550.710 538.98
Recessive (TT+TC vs. CC) 1.176 (0.688–2.008)0.5530.553 540.001.173 (0.687–2.005)0.5590.559 539.96
HWE-P0.4700.999
BMP15 rs17003221 C>T
CC288 (90.9)124 (93.2)1.000 (reference) 1.000 (reference)
CT29 (9.1)9 (6.8)0.721 (0.331–1.568)0.4090.968539.630.718 (0.330–1.563)0.4040.986539.56
TT0 (0.0)0 (0.0)N/AN/AN/AN/AN/AN/AN/AN/A
Dominant (CC vs. CT+TT) 0.721 (0.331–1.568)0.4090.702539.630.718 (0.330–1.563)0.4040.710 539.56
Recessive (CC+CT vs. TT) N/AN/AN/AN/AN/AN/AN/AN/A
HWE-P0.3930.686
BMP15 rs3810682 C>G
CC305 (96.2)129 (97.0)1.000 (reference) 1.000 (reference)
CG12 (3.8)4 (3.0)0.788 (0.250–2.490)0.6850.968540.180.780 (0.246–2.470)0.6720.986540.11
GG0 (0.0)0 (0.0)N/AN/AN/AN/AN/AN/AN/AN/A
Dominant (CC vs. CG+GG) 0.788 (0.250–2.490)0.6850.702540.180.780 (0.246–2.470)0.6720.710 540.11
Recessive (CC+CG vs. GG) N/AN/AN/AN/AN/AN/AN/AN/A
HWE-P0.7310.860
Note: AOR was adjusted by age. RIF, recurrent implantation failure; COR, crude odds ratio; AOR, adjusted odds ratio; 95% CI, 95% confidence interval; FDR; false discovery rate; AIC, Akaike information criterion; HWE, Hardy-Weinberg equilibrium, N/A, not applicable.
Table 3. Comparison of genotype frequencies of FSHR, INHA, ESR1, and BMP15 polymorphisms between the RIF and control subjects.
Table 3. Comparison of genotype frequencies of FSHR, INHA, ESR1, and BMP15 polymorphisms between the RIF and control subjects.
GenotypesControlsRIF ≥ 3AOR (95% CI)pFDR-pAICRIF ≥ 4AOR (95% CI)pFDR-pAIC
(n = 317)(n = 119)(n = 89)
FSHR rs6165 A>G
AA106 (33.4)36 (30.3)1.000 (reference) 28 (31.5)1.000 (reference)
AG156 (49.2)54 (45.4)1.037 (0.665–1.617)0.8710.967393.8742 (47.2)0.955 (0.582–1.566)0.8540.969335.77
GG55 (17.4)29 (24.4)0.412 (0.183–0.931)0.0330.165261.4819 (21.3)0.392 (0.156–0.982)0.0460.115213.89
Dominant (AA vs. AG+GG) 0.852 (0.558–1.303)0.4610.706504.07 0.787 (0.491–1.261)0.3190.809420.25
Recessive (AA+AG vs. GG) 0.373 (0.171–0.815)0.0130.065501.78 0.372 (0.154–0.901)0.0290.073419.62
INHA rs11893842 A>G
AA209 (65.9)75 (63.0)1.000 (reference) 55 (61.8)1.000 (reference)
AG102 (32.2)38 (31.9)1.010 (0.619–1.649)0.9670.967485.2128 (31.5)1.011 (0.589–1.734)0.9690.969399.29
GG6 (1.9)6 (5.0)1.534 (0.851–2.767)0.1550.194335.636 (6.7)1.284 (0.657–2.511)0.4650.465277.85
Dominant (AA vs. AG+GG) 1.147 (0.727–1.811)0.5550.706504.15 1.082 (0.652–1.794)0.7610.809419.89
Recessive (AA+AG vs. GG) 1.536 (0.922–2.557)0.0990.165501.57 1.292 (0.720–2.319)0.3910.391415.63
INHA rs35118453 C>T
CC134 (42.3)55 (46.2)1.000 (reference) 43 (48.3)1.000 (reference)
CT133 (42.0)56 (47.1)1.029 (0.651–1.626)0.9030.967449.9140 (44.9)1.033 (0.617–1.727)0.9030.969376.03
TT50 (15.8)8 (6.7)2.771 (0.863–8.895)0.0870.194265.386 (6.7)3.793 (1.171–12.285)0.0260.115226.80
Dominant (CC vs. CT+TT) 1.123 (0.724–1.744)0.6050.706503.88 1.182 (0.725–1.925)0.5030.809419.35
Recessive (CC+CT vs. TT) 2.749 (0.868–8.704)0.0860.165497.08 3.745 (1.176–11.921)0.0250.073414.45
ESR1 rs9340799 A>G
AA212 (66.9)76 (63.9)1.000 (reference) 56 (62.9)1.000 (reference)
AG94 (29.7)35 (29.4)1.042 (0.652–1.666)0.8630.967475.8226 (29.2)1.054 (0.623–1.783)0.8440.969392.70
GG11 (3.5)8 (6.7)2.060 (0.796–5.328)0.1360.194351.207 (7.9)2.459 (0.907–6.664)0.0770.128291.35
Dominant (AA vs. AG+GG) 1.146 (0.737–1.782)0.5460.706504.06 1.196 (0.732–1.952)0.4750.809419.83
Recessive (AA+AG vs. GG) 2.025 (0.793–5.173)0.1400.175502.35 2.400 (0.900–6.397)0.0800.133417.51
ESR1 rs2234693 T>C
TT123 (38.8)38 (31.9)1.000 (reference) 29 (32.6)
TC144 (45.4)59 (49.6)1.318 (0.820–2.118)0.2540.967413.2442 (47.2)1.227 (0.720–2.089)0.4520.969338.94
CC50 (15.8)22 (18.5)1.393 (0.748–2.594)0.2970.297257.7218 (20.2)1.498 (0.762–2.945)0.2410.301219.97
Dominant (TT vs. TC+CC) 1.342 (0.858–2.100)0.1970.706502.73 1.303 (0.792–2.144)0.2980.809419.24
Recessive (TT+TC vs. CC) 1.197 (0.688–2.082)0.5240.524504.02 1.340 (0.736–2.442)0.3390.391419.45
BMP15 rs17003221 C>T
CC288 (90.9)110 (92.4)1.000 (reference) 82 (92.1)1.000 (reference)
CT29 (9.1)9 (7.6)0.804 (0.369–1.755)0.5840.967504.117 (7.9)0.834 (0.352–1.977)0.6810.969420.16
TT0 (0.0)0 (0.0)N/AN/AN/AN/A0 (0.0)N/AN/AN/AN/A
Dominant (CC vs. CT+TT) 0.804 (0.369–1.755)0.5840.706504.11 0.834 (0.352–1.977)0.6810.809420.16
Recessive (CC+CT vs. TT) N/AN/AN/AN/A N/AN/AN/AN/A
BMP15 rs3810682 C>G
CC305 (96.2)115 (96.6)1.000 (reference) 86 (96.6)1.000 (reference)
CG12 (3.8)4 (3.4)0.853 (0.269–2.709)0.7880.967504.353 (3.4)0.853 (0.234–3.104)0.8090.969420.28
GG0 (0.0)0 (0.0)N/AN/AN/AN/A0 (0.0)N/AN/AN/AN/A
Dominant (CC vs. CG+GG) 0.853 (0.269–2.709)0.7880.788504.35 0.853 (0.234–3.104)0.8090.809420.28
Recessive (CC+CG vs. GG) N/AN/AN/AN/A N/AN/AN/AN/A
Note: AOR was adjusted by age. RIF, recurrent implantation failure, AOR, adjusted odds ratio; 95% CI, 95% confidence interval; AIC, Akaike information criterion; N/A, not applicable.
Table 4. Combined genotype analysis for the FSHR, INHA, ESR1, and BMP15 polymorphisms in RIF patients and controls.
Table 4. Combined genotype analysis for the FSHR, INHA, ESR1, and BMP15 polymorphisms in RIF patients and controls.
Genotype CombinationsControlsRIF AOR (95% CI)p
(n = 317)(n = 133)
FSHR rs6165 A>G/ESR1 rs9340799 A>G
AA/AA92 (29.0)38 (28.6)1.000 (reference)
AA/AG37 (11.7)20 (15.0)1.312 (0.673–2.556)0.425
AA/GG5 (1.6)4 (3.0)2.017 (0.487–8.352)0.333
AG/AA88 (27.8)42 (31.6)1.163 (0.685–1.974)0.576
AG/AG42 (13.2)16 (12.0)1.767 (0.278–11.232)0.546
AG/GG3 (0.9)3 (2.3)0.807 (0.270–2.410)0.701
GG/AA32 (10.1)3 (2.3)0.250 (0.072–0.874)0.030
GG/AG15 (4.7)5 (3.8)0.807 (0.270–2.410)0.701
GG/GG3 (0.9)2 (1.5)1.767 (0.278–11.232)0.546
INHA rs35118453 C>T/ESR1 rs2234693 T>C
CC+CT/TT119 (37.5)42 (31.6)1.000 (reference)
CC+CT/TC+CC192 (60.6)84 (63.2)1.240 (0.802–1.917)0.333
TT/TT4 (1.3)2 (1.5)1.403 (0.247–7.976)0.703
TT/TC+CC2 (0.6)5 (3.8)7.001 (1.298–37.776)0.024
FSHR rs6165 A>G/ESR1 rs9340799 A>G
AA+AG/AA180 (56.8)80 (60.2)1.000 (reference)
AA+AG/AG+GG87 (27.4)43 (32.3)1.115 (0.710–1.750)0.636
GG/AA32 (10.1)3 (2.3)0.199 (0.059–0.670)0.009
GG/AG+GG18 (5.7)7 (5.3)0.867 (0.348–2.160)0.758
FSHR rs6165 A>G/BMP15 rs17003221 C>T
AA+AG/CC241 (76.0)114 (85.7)1.000 (reference)
AA+AG/CT+TT26 (8.2)9 (6.8)0.723 (0.328–1.594)0.421
GG/CC47 (14.8)10 (7.5)0.447 (0.218–0.918)0.028
GG/CT+TT3 (0.9)0 (0.0)N/AN/A
FSHR rs6165 A>G/BMP15 rs3810682 C>G
AA+AG/CC257 (81.1)119 (89.5)1.000 (reference)
AA+AG/CG+GG10 (3.2)4 (3.0)0.842 (0.258–2.745)0.775
GG/CC48 (15.1)10 (7.5)0.446 (0.218–0.914)0.027
GG/CG+GG2 (0.6)0 (0.0)N/AN/A
INHA rs11893842 A>G/ESR1 rs9340799 A>G
AA+AG/AA+AG255 (80.4)95 (71.4)1.000 (reference)
AA+AG/GG7 (2.2)8 (6.0)3.065 (1.081–8.690)0.035
GG/AA+AG51 (16.1)29 (21.8)1.526 (0.913–2.549)0.107
GG/GG4 (1.3)1 (0.8)0.670 (0.074–6.073)0.722
FSHR rs6165A>G/INHA rs35118453C>T
AA+AG/CC+CT261 (82.3)117 (88.0)
AA+AG/TT6 (1.9)6 (4.5)2.203 (0.695–6.979)0.180
GG/CC+CT50 (15.8)9 (6.8)0.399 (0.190–0.841)0.016
GG/TT0 (0.0)1 (0.8)N/AN/A
FSHR rs6165 A>G/ESR1 rs9340799 A>G
AA+AG/AA+AG259 (81.7)116 (87.2)1.000 (reference)
AA+AG/GG8 (2.5)7 (5.3)1.979 (0.700–5.599)0.198
GG/AA+AG47 (14.8)8 (6.0)0.373 (0.171–0.816)0.014
GG/GG3 (0.9)2 (1.5)1.450 (0.239–8.817)0.687
FSHR rs6165 A>G/ESR1 rs2234693 T>C
AA+AG/TT+TC224 (70.7)101 (75.9)1.000 (reference)
AA+AG/CC43 (13.6)22 (16.5)0.604 (0.123–2.966)0.686
GG/TT+TC43 (13.6)8 (6.0)0.401 (0.181–0.886)0.024
GG/CC7 (2.2)2 (1.5)0.604 (0.123–2.966)0.534
FSHR rs6165 A>G/BMP15 rs17003221 C>T
AA+AG/CC+CT267 (84.2)123 (92.5)1.000 (reference)
GG/CC+CT50 (15.8)10 (7.5)0.430 (0.210–0.877)0.020
FSHR rs6165 A>G/BMP15 rs3810682 C>G
AA+AG/CC+CG267 (84.2)123 (92.5)1.000 (reference)
GG/CC+CG50 (15.8)10 (7.5)0.430 (0.210–0.877)0.020
ESR1 rs9340799 A>G/ESR1 rs2234693 T>C
AA+AG/TT+TC266 (83.9)104 (78.2)1.000 (reference)
AA+AG/CC40 (12.6)20 (15.0)1.270 (0.709–2.276)0.422
GG/TT+TC1 (0.3)5 (3.8)12.930 (1.492–112.052)0.020
GG/CC10 (3.2)4 (3.0)1.024 (0.314–3.340)0.968
Note: AOR was adjusted by age. RIF, recurrent implantation failure; 95% CI, 95% confidence interval; AOR, adjusted odds ratio; NA, not applicable.
Table 5. Clinical variables in RIF patients and controls stratified by FSHR, INHA, ESR1, and BMP15 polymorphisms status by ANOVA and Kruskal–Wallis test.
Table 5. Clinical variables in RIF patients and controls stratified by FSHR, INHA, ESR1, and BMP15 polymorphisms status by ANOVA and Kruskal–Wallis test.
GenotypesBMI (kg/m2) Homocysteine (μmol/L)Folate (ng/mL)BUN (mg/dL)Creatinine (mg/dL)Uric Acid (mg/dL)Total Cholesterol (mg/dL)TSH (mU/L)
Mean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SD
FSHR rs6165 A>G
AA21.4 ± 3.67.0 ± 2.316.8 ± 11.89.2 ± 3.00.7 ± 0.13.9 ± 1.0199.1 ± 53.01.9 ± 1.4
AG21.5 ± 3.26.3 ± 1.813.3 ± 6.29.5 ± 2.90.7 ± 0.23.9 ± 1.0209.0 ± 55.61.9 ± 1.2
GG21.5 ± 2.47.0 ± 3.89.0 ± 4.09.5 ± 2.90.7 ± 0.24.0 ± 1.2202.0 ± 52.22.3 ± 1.2
pb0.6760.2730.1750.5330.5760.9870.3660.209
INHA rs11893842 A>G
AA21.4 ± 3.26.4 ± 1.813.9 ± 8.79.5 ± 3.00.7 ± 0.24.0 ± 1.0205.2 ± 57.22.1 ± 1.6
AG21.7 ± 3.66.9 ± 2.516.2 ± 12.29.3 ± 2.90.7 ± 0.13.8 ± 1.0205.6 ± 53.31.7 ± 1.0
GG20.9 ± 2.67.0 ± 2.413.9 ± 5.59.4 ± 2.80.7 ± 0.24.0 ± 1.1198.7 ± 52.11.9 ± 1.3
pb0.770.576 a0.8380.8730.7890.3810.5320.332
INHA rs35118453 C>T
CC21.3 ± 3.16.6 ± 1.914.7 ± 10.79.4 ± 3.10.7 ± 0.23.9 ± 1.1203.1 ± 56.52.0 ± 1.4
CT21.8 ± 3.97.1 ± 2.814.4 ± 7.89.2 ± 2.70.7 ± 0.24.0 ± 0.9206.6 ± 51.41.8 ± 1.0
TT19.8 ± 1.76.3 ± 1.318.3 ± 6.59.9 ± 2.20.7 ± 0.23.6 ± 1.0196.8 ± 34.72.1 ± 1.3
pb0.3170.7160.3890.5590.7460.5310.7850.789
ESR1 rs9340799 A>G
AA21.6 ± 3.66.7 ± 1.914.3 ± 7.69.4 ± 3.00.7 ± 0.23.9 ± 1.0200.7 ± 51.71.8 ± 1.2
AG21.0 ± 2.67.3 ± 3.017.9 ± 14.69.3 ± 2.80.7 ± 0.23.9 ± 1.1214.5 ± 59.52.1 ± 1.4
GG22.3 ± 4.25.6 ± 2.711.1 ± 7.79.7 ± 2.90.7 ± 0.14.1 ± 0.8184.5 ± 42.42.2 ± 1.3
pb0.6620.276 a0.5020.8470.7260.6180.1160.42
ESR1 rs2234693 T>C
TT21.8 ± 4.16.5 ± 1.612.9 ± 6.09.0 ± 3.00.7 ± 0.14.0 ± 1.2204.8 ± 51.01.9 ± 1.4
TC21.2 ± 2.76.8 ± 2.215.4 ± 11.49.5 ± 3.00.7 ± 0.23.8 ± 0.9205.5 ± 58.81.9 ± 1.3
CC21.5 ± 3.07.3 ± 3.418.3 ± 10.99.9 ± 2.50.7 ± 0.24.0 ± 1.0197.6 ± 45.91.8 ± 1.0
pb0.9740.4560.2940.0850.760.920.8720.927
BMP15 rs17003221 C>T
CC21.4 ± 3.36.7 ± 2.315.1 ± 9.79.4 ± 3.00.7 ± 0.23.9 ± 1.1204.5 ± 55.11.9 ± 1.3
CT22.4 ± 4.07.5 ± 1.08.5 ± 1.09.1 ± 2.30.7 ± 0.23.9 ± 0.8199.4 ± 44.71.8 ± 1.0
TT0.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.0
pb0.243 a0.2380.1490.8320.4440.9020.9770.94
BMP15 rs3810682 C>G
CC21.5 ± 3.36.7 ± 2.315.0 ± 9.79.4 ± 3.00.7 ± 0.23.9 ± 1.0203.8 ± 54.31.9 ± 1.3
CG20.6 ± 3.57.1 ± 1.310.9 ± 3.99.2 ± 1.80.7 ± 0.24.0 ± 0.9208.4 ± 51.41.6 ± 1.1
GG0.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.0
pb0.303 a0.6360.5370.9090.9280.6740.7850.741
GenotypesE2 (pg/mL)FSH (U/L)LH (U/L)WBC (103/μL)Hgb (g/dL)PLT (103/μL)PT (sec)aPTT (sec)
Mean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SDMean ± SD
FSHR rs6165 A>G
AA29.0 ± 14.48.5 ± 4.44.2 ± 2.47.9 ± 2.612.4 ± 1.3235.6 ± 62.311.1 ± 1.829.5 ± 3.5
AG33.0 ± 16.88.2 ± 3.24.3 ± 2.97.9 ± 2.512.4 ± 1.4231.1 ± 60.510.8 ± 0.929.3 ± 3.5
GG25.9 ± 12.98.8 ± 2.24.0 ± 1.56.8 ± 2.112.0 ± 1.5234.0 ± 89.910.7 ± 1.028.4 ± 2.8
pb0.0760.110.8430.1190.4290.5760.2070.334
INHA rs11893842 A>G
AA27.2 ± 13.38.6 ± 3.53.9 ± 2.07.6 ± 2.312.4 ± 1.2236.4 ± 73.810.8 ± 0.829.3 ± 3.9
AG30.4 ± 16.18.1 ± 4.24.1 ± 2.97.9 ± 2.712.3 ± 1.4233.5 ± 64.711.0 ± 1.829.1 ± 3.2
GG34.0 ± 15.79.0 ± 1.74.7 ± 2.27.6 ± 2.212.3 ± 1.4228.2 ± 47.510.8 ± 0.829.7 ± 3.4
pb0.0670.1130.0770.6820.8120.970.7880.521
INHA rs35118453 C > T
CC29.0 ± 15.08.5 ± 4.04.1 ± 2.77.6 ± 2.412.3 ± 1.3236.6 ± 67.610.9 ± 0.929.3 ± 3.5
CT31.0 ± 15.48.4 ± 3.04.3 ± 1.97.9 ± 2.612.4 ± 1.4229.0 ± 60.411.0 ± 2.029.3 ± 3.5
TT44.9 ± 15.07.1 ± 2.86.1 ± 1.78.6 ± 2.512.6 ± 1.6217.6 ± 50.110.7 ± 0.629.8 ± 2.8
pb0.0670.0660.0170.2320.2620.4820.8780.634
ESR1 rs9340799 A>G
AA30.5 ± 16.28.5 ± 4.04.3 ± 2.77.8 ± 2.512.4 ± 1.4234.6 ± 69.610.9 ± 0.829.2 ± 3.4
AG28.8 ± 13.48.4 ± 3.04.0 ± 2.17.8 ± 2.512.4 ± 1.3233.3 ± 55.811.0 ± 2.229.4 ± 3.6
GG30.2 ± 12.67.1 ± 2.83.3 ± 1.16.9 ± 1.912.0 ± 1.2219.8 ± 48.411.1 ± 1.029.7 ± 3.7
pb0.8820.2670.5760.4230.3810.6440.510.918
ESR1 rs2234693 T>C
TT28.7 ± 14.28.4 ± 3.64.1 ± 3.07.9 ± 2.612.4 ± 1.5232.3 ± 59.510.9 ± 0.929.3 ± 3.7
TC30.9 ± 15.58.1 ± 3.34.3 ± 2.37.7 ± 2.512.4 ± 1.2237.7 ± 71.011.0 ± 1.829.4 ± 3.4
CC30.0 ± 16.89.2 ± 4.74.0 ± 1.87.7 ± 2.312.2 ± 1.3223.1 ± 55.810.8 ± 0.829.0 ± 3.1
pb0.580.4460.6010.8610.7610.5430.7750.625
BMP15 rs17003221 C > T
CC30.1 ± 15.08.5 ± 3.74.2 ± 2.57.8 ± 2.612.4 ± 1.4233.2 ± 64.210.9 ± 1.429.4 ± 3.5
CT27.0 ± 19.87.8 ± 2.54.0 ± 0.97.7 ± 1.612.0 ± 1.0234.9 ± 69.610.8 ± 0.828.5 ± 3.1
TT0.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.0
p b0.550.50.5240.680.0390.8120.5530.146
BMP15 rs3810682 C>G
CC29.7 ± 15.28.4 ± 3.74.2 ± 2.57.7 ± 2.512.3 ± 1.3233.4 ± 65.110.9 ± 1.429.3 ± 3.4
CG35.6 ± 15.68.4 ± 3.44.4 ± 1.58.7 ± 2.212.4 ± 1.8231.9 ± 53.610.6 ± 0.630.0 ± 5.2
GG0.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.00.0 ± 0.0
Pb0.2480.9020.340.1050.680.7320.3630.761
ANOVA, Analysis of variance; BMI, body mass index; BUN, blood urea nitrogen; TSH, thyroid stimulating hormone; E2, estradiol; FSH, follicle stimulating hormone; LH, luteinizing hormone; WBC, white blood cell; Hgb, hemoglobin; PLT, platelet; PT, prothrombin time; aPTT, activated partial thromboplastin time; SD, standard deviation. a p-values were calculated by ANOVA. b p-values were calculated by Kruskal–Wallis tests.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ko, E.-J.; Shin, J.-E.; Lee, J.-Y.; Ryu, C.-S.; Hwang, J.-Y.; Kim, Y.-R.; Ahn, E.-H.; Kim, J.-H.; Kim, N.-K. Association of Polymorphisms in FSHR, INHA, ESR1, and BMP15 with Recurrent Implantation Failure. Biomedicines 2023, 11, 1374. https://doi.org/10.3390/biomedicines11051374

AMA Style

Ko E-J, Shin J-E, Lee J-Y, Ryu C-S, Hwang J-Y, Kim Y-R, Ahn E-H, Kim J-H, Kim N-K. Association of Polymorphisms in FSHR, INHA, ESR1, and BMP15 with Recurrent Implantation Failure. Biomedicines. 2023; 11(5):1374. https://doi.org/10.3390/biomedicines11051374

Chicago/Turabian Style

Ko, Eun-Ju, Ji-Eun Shin, Jung-Yong Lee, Chang-Soo Ryu, Ji-Young Hwang, Young-Ran Kim, Eun-Hee Ahn, Ji-Hyang Kim, and Nam-Keun Kim. 2023. "Association of Polymorphisms in FSHR, INHA, ESR1, and BMP15 with Recurrent Implantation Failure" Biomedicines 11, no. 5: 1374. https://doi.org/10.3390/biomedicines11051374

APA Style

Ko, E. -J., Shin, J. -E., Lee, J. -Y., Ryu, C. -S., Hwang, J. -Y., Kim, Y. -R., Ahn, E. -H., Kim, J. -H., & Kim, N. -K. (2023). Association of Polymorphisms in FSHR, INHA, ESR1, and BMP15 with Recurrent Implantation Failure. Biomedicines, 11(5), 1374. https://doi.org/10.3390/biomedicines11051374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop