Association of FSHR, LH, LHR, BMP15, GDF9, AMH, and AMHR polymorphisms with poor ovarian response in patients undergoing in vitro fertilization

Arivaldo José Conceição Meireles; João Paolo Bilibio; Pânila Longhi Lorenzzoni; Emily De Conto; Fábio Costa do Nascimento; João Sabino da Cunha-Filho

JBRA Assist. Reprod. 2021;25(3):439-446
ORIGINAL ARTICLE
doi: 10.5935/1518-0557.20210004

Association of FSHR, LH, LHR, BMP15, GDF9, AMH, and AMHR polymorphisms with poor ovarian response in patients undergoing in vitro fertilization

Arivaldo José Conceição Meireles^1,2, João Paolo Bilibio^1,3, Pânila Longhi Lorenzzoni^1,2, Emily De Conto^2,4, Fábio Costa do Nascimento¹, João Sabino da Cunha-Filho^2,4

¹Clinica Pronatus Centro de Reprodução Humana, Belém, Pará, Brazil
²Programa de Pós Graduação de Ciências Médicas da Universidade Federal do Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil
³Department of Obstetrics and Gynecology, Universidade Federal do Pará, Belém, Pará, Brazil
⁴Clínica Insemine Centro de Reprodução Humana, Porto Alegre, Brazil

Received July 24, 2020
Accepted January 13, 2020

Corresponding Author:
João Paolo Bilibio
Department of Obstetrics and Gynecology
Universidade Federal do Pará
Belém - PA, Brazil
E-mail: joaopaolobilibio@yahoo.com.br
ORCID: J.P. Bilibio - 0000-0003-4755-170X

CONFLICT OF INTEREST
The authors have no conflict of interest to declare.

ABSTRACT
Objective: This paper aimed to assess the correlation between LH, LHR, GDF9, FSHR, AMH, AMHR2, and BMP15 polymorphisms, which are related to follicular development, and decreased ovarian response in women undergoing controlled ovarian hyperstimulation (COH) for IVF.
Methods: This age-matched case-control study included three or four controls per woman undergoing COH. Controls were women with normal ovarian response (NOR) and cases were women with poor ovarian response (POR) in oocyte retrieval (three or fewer oocytes). DNA was extracted from peripheral blood and potential associations with gene polymorphisms related to follicular development (LH, LHR, GDF9, FSHR, AMH, AMHR2, and BMP15) were analyzed.
Results: Sixty-six patients were included, 52 in the NOR and 14 in the POR group. Two GDF9 polymorphisms were associated with follicular response after COH, one associated with POR - the presence of a mutant polymorphism in heterozygosis and homozygosis of the GDF9 398-39 (C to G) [23% NOR versus 68% POR (OR 4.01, CI 1.52-10.6, p=0.005)] - and another associated with protective response - the presence of normal homozygosis of GDF9 (C447T) [19.2% NOR versus 50% POR (OR 0.34, IC 0.14-0.84, p=0.019)]. No additional associations were found between the other analyzed polymorphisms and POR.
Conclusions: This study found that GDF9 appears to play an important role in follicular development, whereas polymorphisms in its DNA chain may negatively affect ovarian reserve, such as 398-39 (C to G), or positively, as seen in C447T.

Keywords: poor ovarian response, polymorphisms, ovarian reserve, IVF, SNPs

INTRODUCTION
Approximately 10% of women seeking fertility treatment have diminished ovarian reserve (DOR), defined as a decreased number or quality oocytes (Ferraretti et al., 2011). Poor ovarian response (POR) is defined as meeting two or three of the following criteria: age greater than 39 years; prior POR to conventional stimulation protocols (less than three oocytes retrieved); and abnormal ovarian reserve testing (Ferraretti et al., 2011). DOR can have multiple etiologies, including autoimmunity, idiopathic, iatrogenic, and genetic causes, and single nucleotide polymorphisms (SNPs) (Goswami & Conway, 2005; Greene et al., 2014). The impact of mutations and polymorphic variant genes involved in folliculogenesis is still uncertain.
Polymorphisms exist within genes and are a source of variation between individuals. Gene association studies have identified a number of SNPs that affect gonadotropin, steroid, and transforming growth factor beta (TGFβ) pathways, which include TGF-β, anti-Müllerian hormone (AMH), activins, inhibins, bone morphogenetic proteins (BMPs), and growth differentiation factors (GDFs) involved in ovarian response (Greene et al., 2014). Most of them affect mRNA levels or the protein sequence, and thus lead to quantitative functional protein variations that may account for the observed inter-individual variability in controlled ovarian hyperstimulation (COH) (Levi et al., 2001; Greene et al., 2014).
Growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP-15), members of the TGFβ superfamily, are potent regulators of folliculogenesis and ovulation and are both expressed in oocytes from early stage follicles (Aaltonen et al., 1999; Chang et al., 2016; Sanfins et al., 2018). With respect to COH phenotypes, GDF9 and BMP15 alleles have been associated with stimulation outcome (Moron et al., 2006; Wang et al. 2010a; Hanevik et al., 2011; Bilibio et al., 2020). Several genetic variants of GDF9 have been identified, and their correlation with POR has been noted, suggesting that these variants contribute to aberrant follicular development and oocyte loss (Di Pasquale et al., 2004; Shimizu et al., 2004; Abir et al., 2008; Wang et al., 2010a;b). Other variant alleles of BMP15 were associated with increased follicle production in COH, such as the 9C>G polymorphism, which was associated with high response to ovarian stimulation (Moron et al., 2006; Hanevik et al., 2011).
The follicle stimulating hormone receptor (FSHR) carries more than two thousand SNPs, and one polymorphism, Asn608Ser, was found to have an important impact on ovarian response (Livshyts et al., 2009; Sheikhha et al., 2011; Pabalan et al., 2014). Homozygous mutant women with this polymorphism have higher baseline FSH levels, require higher FSH doses during COS, and present lower estradiol levels than heterozygous and normal homozygous women. Another study showed that women homozygous for the Ser680 variant had a greater number of mature oocytes than women homozygous for the Asn680 variant (Lawson et al., 2003; Ferrarini et al., 2013).
AMH (Ile49Ser) and AMHR polymorphisms (482A>G) have been associated with variations in estradiol levels and may modulate FSH sensitivity (de Boer et al., 2002). In regard to LH polymorphisms, variant 8Arg-15Thr in women undergoing IVF is more frequently seen among poor responders to rFSH and women with ovarian resistance to rFSH who, despite requiring higher rFSH dosages, have fewer oocytes retrieved (Alviggi et al., 2009).
This paper aimed to assess the correlation between LH, LHR, GDF9, FSHR, AMH, AMHR2, and BMP15 polymorphisms, which are related to follicular development, and decreased ovarian response in women undergoing COH for IVF.

MATERIALS AND METHODS
This age-matched case-control study included three or four controls per case involving women undergoing COH for IVF. Controls were women with normal ovarian response (NOR) and cases were women with POR in oocyte retrieval (three or fewer oocytes). The study was carried out at the Pronatus Assisted Reproduction Center and at the Federal Univesity of Rio Grande do Sul. The study was approved by the National Ethics and Research Committee on Human Beings and granted certificate no. 25525413.0.0000.5327 by the Ethics Committee of the Hospital de Clínicas de Porto Alegre (an Institutional Review Board equivalent).
Sixty-six patients were included in the study. The inclusion criteria were as follows: age between 30 and 39 years; presence of both ovaries; diagnosis of infertility (more than one year) by male factor, tubal factor (confirmed by hysterosalpingography or videolaparoscopy), or unexplained causes; and no prior IVF treatment. Three or four controls were matched by age to each case. The exclusion criteria were as follows: endocrine disorders; polycystic ovary syndrome; prior chemotherapy; abnormal karyotype; and other factors affecting ovarian function, such as ovarian surgery and endometrioma.
The patients underwent antral follicle counts on Day 2 or 3 of the menstrual cycle and were tested for FSH, LH, estradiol, and anti-Müllerian hormone (AMH) levels. Ovarian stimulation was initiated on Day 3 of the menstrual cycle with recombinant FSH (Puregon, Organon) and included r-FSH 200 IU/day for the first three days and then 150 IU/day; after Day 6, a GnRH antagonist, Orgalutran (Organon), was administered daily by S.C. injection (0.25 mg/d) until the administration of recombinant human chorionic gonadotropin (rhCG) (Ovidrel 250 mg/0.5 ml, Serono). When three or more follicles reached 17 mm, oocyte maturation was stimulated with rhCG, and on this day, all follicles and serum FSH, LH, and P4 levels were measured. Ultrasound-guided oocyte retrieval was performed after 35 h and, following denudation, the oocytes were categorized as metaphase II (MII), metaphase I (MI), or prophase I (PI) stage. After oocyte retrieval, the patients were split into two groups: NOR (four or more oocytes) or POR (three or fewer oocytes).
Whole blood samples were collected and aliquots of 350 µL from each sample were used for genomic DNA extraction using the Easy DNA kit according to manufacturer’s instructions (Invitrogen, UK). Concentration and purity were determined by spectrophotometry NanoDrop ND1000 (NanoDrop Technologies Inc, USA). DNA was used for analysis of polymorphisms in the LHβ, LHR, GDF9, FSHR, AMH, AMHR2, and BMP15 genes (Table 1).

Table 1. Studied gene polymorphisms

Polymerase chain reaction (PCR) was performed to amplify the regions of interest in the LHβ, LHR, FSHR and GDF9 genes (Table 1). Forward and reverse primers delimited the region where the polymorphisms were located, as described in Table 1. Amplification of DNA at a concentration of 200 ng was performed in a thermocycler (Veriti^® 96-Well Thermal Cycler, Applied Biosystems, USA) with reagent Invitrogen (UK). The following annealing temperatures were applied: -63°C, LHR; -64°C, GDF9; -63°C, FSHR; and LH, 60°C. PCR products were stained with SYBR^® Gold (Life Technologies, USA) and verified by electrophoresis on 1.5% agarose gel, and the DNA fragments obtained had the following sizes: LHβ - 662 base pairs, LHR - 350 pairs bases, GDF9 - 491 base pairs and FSHR - 520 base pairs.
The PCR products (detection of polymorphisms of LHβ, LHR and GDF9 - direct sequencing) were purified using the purification protocol with PEG8000 and 2.5% NaCl and subjected to direct sequencing by the Sanger method on an automated sequencer ABI 3500 Genetic Analyzer (Applied Biosystems, USA). Assays were performed at the Unit of Molecular and Protein Analysis of Hospital de Clinicas de Porto Alegre (UAMP/HCPA). Primers for sequencing were used at a concentration of 4 pmol/µL. The results were compared with NCBI reference sequences LHβ - NM_000894.2, LHR - NM_000233.3, and GDF9 - NM_005260.3 (Table 1).
The product of the PCR for the FSHR gene underwent RFLP using restriction enzyme BSRI for four hours at 65 ° C. After digestion, the samples were stained with SYBR ^® Gold (Life Technologies, USA) and the results obtained were determined by electrophoresis on 2.5% agarose gel: NN (680Asn/Asn) fragments of 520 base pairs; SS (680Ser/Ser) fragments of 413 and 107 base pairs; and NS (680Asn/Ser) fragments of 520, 413 and 107 base pairs. The results were compared with the NCBI reference sequence NM_000145.3 (Table 1).
The polymorphisms in the AMH, AMHR2, and BMP15 genes were determined by TaqMan Allelic Discrimination - Real-Time PCR from DNA at a concentration of 10-20 ng. The tests shown in Table 1 (Applied Biosystems, USA) were performed on a StepOne ™ Real-Time PCR System (Applied Biosystems, USA) at the Unit of Analysis Molecular and Proteins Hospital de Clinicas de Porto Alegre (UAMP/HCPA).
Observed numbers for each genotype were compared against expected values to test whether the samples were in Hardy-Weinberg equilibrium using the Chi-Squared test with one degree of freedom. Data were tested for normality using the Kolmogorov-Smirnov test, and data with a normal distribution, as determined by the t-test, were used to compare means. The Kruskal-Wallis test was used for nonparametric data. Qualitative variables were analyzed with the Chi-squared test. Statistical analyses to determine the association between the polymorphisms were performed using the Chi-squared test and the Monte Carlo method for multiple simulations assuming that the null hypothesis was correct with fixed marginal values. Marginal values were determined from experimental data. This method determines an empirical P value for the observed data. Moreover, we utilized advanced statistical analysis to evaluate the effect of polymorphisms using logistic regression to predict the effect on POR. Statistical significance was set at p<0.05. Statistical tests were performed on the Statistical Package for the Social Sciences 23 (SPSS Inc., Chicago, IL).

RESULTS
The clinical characteristics of patients undergoing COH for IVF are presented in Table 2. Mean age was similar (34.6 NOR versus 35.5 POR, p=0.271), and the women with NOR were heavier than the individuals with POR (65.9 NOR versus 56.4 POR, p=0.001).

Table 2. Clinical characteristics of patients submitted to controlled ovarian hyperstimulation for in vitro fertilization (mean±SD)

Ovarian and hormonal characteristics of patients submitted to COH for IVF are presented in Table 3. Antral follicle count (12.3 NOR versus 4.2 POR, p<0.001), number of follicles greater than 17 mm at rhCG day (7.1 NOR versus 2.0 POR, p<0.001), MII oocyte count (9.8 NOR versus 1.7 POR, p<0.001), and AMH (2.0 NOR versus 0.6 POR, p=0.017) were lower in the POR group.

Table 3. Ovarian and hormonal characteristics of patients submitted to controlled ovarian hyperstimulation for in vitro fertilization (mean±SD)

The test for Hardy-Weinberg equilibrium was performed via the Chi-squared test, and the following results were found: p=0.51 for GFD-9 398; p=0.33 for GFD-9 C447T; p=0.74 for GFD-9 546; p=0.51 for FSHR (Asn680Ser); p=0.60 for LH (Trp8Arg); p=0.60 for LH (Ile15Thr); p=0.99 for LHR (18isnLQ); p=0.16 for AMH (Ile49Ser); p=0.91 for AMHR2 (482A>G); and p=0.98 for BMP15 (9C>G) polymorphisms.
The genotype and allele frequencies of polymorphisms in the NOR and POR groups are demonstrated in Table 4. Two polymorphisms of GDF9 were associated with follicular response after COH. One was associated with POR: the presence of a mutant polymorphism in heterozygosis and homozygosis of the GDF9 398-39 (C to G) (23% NOR versus 68% POR; OR 4.01, IC 1.52-10.6, p=0.005); and the presence of normal homozygosis of GDF9 (C447T) were associated with protective response (19.2% NOR versus 50% POR; OR 0.34, IC 0.14-0.84, p=0.019). There was no association between other studied polymorphisms and POR.

Table 4. Genotype and allele frequencies of polymorphisms in women submitted controlled ovarian hyperstimulation for IVF with normal ovarian response (NOR) and poor ovarian response (POR)

When we evaluated our data with logistic regression to investigate the influence of polymorphisms, only the GDF9 398 polymorphism was associated with POR after controlling for bias from other polymorphisms (Table 5). The logistic regression with other parameters of poor responders showed that AMH, gonadotropin, and the GDF9 C447 polymorphism were associated with POR.

Table 5. Logistic regression - influence of polymorphisms on ovarian response after controlled ovarian hyperstimulation for IVF

DISCUSSION
The present study revealed that GDF9 polymorphisms were associated with follicular response after COH for IVF. These polymorphisms can have a negative impact, i.e., the presence the allele mutant genotype in heterozygosis or homozygosis of GDF9 398-39 (C to G) was found in 69% of the individuals in the POR group and in only 23% of the women in the NOR group (OR 4.01). They may also have a protective response, i.e., the presence of normal homozygosis of GDF9 (C447T) was associated with POR (19.2% NOR versus 50% POR).
Multivariate analysis showed the possible influence of polymorphisms of genes regulating follicular growth and oocyte development (FSHR, LH, LHR, AMH, AMH2, BMP15 and GDF9). Response to administration of gonadotropins during ovarian hyperstimulation for IVF control showed that only the GDF9 398-39 (C to G) polymorphism was associated with POR, suggesting that this polymorphism plays an important role in the status of ovarian reserve and function.
We did not find correlations between the other studied polymorphisms, including GDF9 (G546A and G646A), FSHR (Asn680Ser), LH (Trp8Arg and Ile15Thr), LHR (18isnLQ), AMH (Ile49Ser), AMHR2 (482A>G), or BMP15 (9C>G), and POR. Unlike our study, other authors described associations between the FSHR (Livshyts et al., 2009; Sheikhha et al., 2011) and LHCGR polymorphisms (Skiadas et al., 2012) and POR. A possible cause for this difference is the variety of ethnic groups within the populations included in each study.
GDF9 plays a pivotal role during early folliculogenesis, and deletion of GDF9 in mice causes follicular arrest at the primary stage and infertility (Dong et al., 1996). GDF9 also stimulates granulosa cell proliferation (Vitt et al., 2000), cumulus cell expansion (Yan et al., 2001), follicular apoptosis inhibition (Orisaka et al., 2006), and enhancement of oocyte and embryo development (Hussein et al., 2006; Yeo et al., 2008). In vitro studies using recombinant GDF9 protein have clarified the biological roles and the importance of GDF9 in follicular growth and development in all stages of folliculogenesis. In the preantral stage, GDF9 effectively stimulated the growth of in vitro cultured preantral follicles (Hayashi et al., 1999). GDF9 also promotes early preantral follicular growth in human ovaries (Hreinsson et al., 2002). In the transition to the antral stage, it appears that GDF9 promotes follicular survival by suppressing granulosa cell apoptosis and follicular atresia (Orisaka et al., 2006). This may be achieved in part by GDF9 stimulation of FSHR expression, since adequate FSRH levels in granulosa cells are essential for FSH-dependent antral follicle growth.
GDF9 also plays an important role during the final stages of follicular growth prior to ovulation. Prior to LH surge, cumulus cells require GDF9 to support metabolic cascades such as glycolysis and sterol biosynthesis (Sugiura et al., 2005). GDF9 also regulates diverse processes and gene expression during the preovulatory stage (Elvin et al., 2000) and enhances cumulus cell expansion in the presence of FSH (Elvin et al., 1999), but not without FSH (Dragovic et al., 2005). Although animal models, in vitro studies, and humans studies have revealed the role of GDF9 in regulating follicular development, little is known about it in human ovarian function, since studies have demonstrated a strong correlation between GDF9 polymorphism in women with DOR and poor ovarian response followed by poor IVF outcomes, indicating that GDF9 plays an important role in determining ovarian reserve status and function (Wang et al., 2010b; 2013; Sanfins et al., 2018; Bilibio et al., 2020). In this study, we saw that the number of antral follicles, total follicles, oocytes retrieved, AMH, and estradiol day rhCG decreased in poor responders, suggesting that the GDF9 polymorphism affected oocyte development in all stages of folliculogenesis.
This study found a valuable association of polymorphisms of GDF9398 and ovarian response in patients undergoing ovarian stimulation for IVF in a relatively young and ethnically heterogeneous population. However, the total number of patients participating in this study must be considered. This is an element of great importance, since despite the number of included patients, we found significant results in an ethnically heterogeneous population. Our findings suggest that this polymorphism may be a predictor of ovarian response in other populations.

CONCLUSION
The GDF9 polymorphism (398-39; C to G and C447T) exerts an important influence on oocyte development. GDF-9 seems to influence all stages of folliculogenesis by decreasing follicular recruitment, decreasing ovarian response, or decreasing the number of MII specimens after ovarian hyperstimulation for IVF. Functional studies are required to understand the real function and a possible correction and treatment to improve the number of oocytes retrieved in patients with this polymorphism.

Funding
Fundo de Incentivo a Pesquisa (FIPE), Hospital de Clínicas de Porto Alegre, Brazil.

ACKNOWLEDGEMENTS
We are grateful for the funding provided by the Fundo de Incentivo à Pesquisa (FIPE), Hospital de Clínicas de Porto Alegre, Brazil.

REFERENCES
Aaltonen J, Laitinen MP, Vuojolainen K, Jaatinen R, Horelli-Kuitunen N, Seppä L, Louhio H, Tuuri T, Sjöberg J, Bützow R, Hovata O, Dale L, Ritvos O. Human growth differentiation factor 9 (GDF-9) and its novel homolog GDF-9B are expressed in oocytes during early folliculogenesis. J Clin Endocrinol Metab. 1999;84:2744-50. PMID: 10443672 DOI: 10.1210/jc.84.8.2744
Medline Crossref

Abir R, Ben-Haroush A, Melamed N, Felz C, Krissi H, Fisch B. Expression of bone morphogenetic proteins 4 and 7 and their receptors IA, IB, and II in human ovaries from fetuses and adults. Fertil Steril. 2008;89:1430-40. PMID: 17624341 DOI: 10.1016/j.fertnstert.2007.04.064
Medline Crossref

Alviggi C, Clarizia R, Pettersson K, Mollo A, Humaidan P, Strina I, Coppola M, Ranieri A, D’Uva M, De Placido G. Suboptimal response to GnRHa long protocol is associated with a common LH polymorphism. Reprod Biomed Online. 2009;18:9-14. PMID: 19146763 DOI: 10.1016/S1472-6483(10)60418-X
Medline Crossref

Bilibio JP, Meireles AJC, Conto E, Lorenzzoni PL, Nascimento FCD, Cunha-Filho JSD. GDF9 polymorphisms: influence on ovarian response in women undergoing controlled ovarian hyperstimulation. JBRA Assist Reprod. 2020;24:447-53. PMID: 32510893 DOI: 10.5935/1518-0557.20200027
Medline Crossref

Chang HM, Qiao J, Leung PC. Oocyte-somatic cell interactions in the human ovary-novel role of bone morphogenetic proteins and growth differentiation factors. Hum Reprod Update. 2016;23:1-18. PMID: 27797914 DOI: 10.1093/humupd/dmw039
Medline Crossref

de Boer EJ, den Tonkelaar I, te Velde ER, Burger CW, Klip H, van Leeuwen FE; OMEGA-project group. A low number of retrieved oocytes at in vitro fertilization treatment is predictive of early menopause. Fertil Steril. 2002;77:978-85. PMID: 12009354 DOI: 10.1016/S0015-0282(02)02972-2
Medline Crossref

Di Pasquale E, Beck-Peccoz P, Persani L. Hypergonadotropic ovarian failure associated with an inherited mutation of human bone morphogenetic protein-15 (BMP15) gene. Am J Hum Genet. 2004;75:106-11. PMID: 15136966 DOI: 10.1086/422103
Medline Crossref

Dong J, Albertini DF, Nishimori K, Kumar TR, Lu N, Matzuk MM. Growth differentiation factor-9 is required during early ovarian folliculogenesis. Nature. 1996; 383:531-5. PMID: 8849725 DOI: 10.1038/383531a0
Medline Crossref

Dragovic RA, Ritter LJ, Schulz SJ, Amato F, Armstrong DT, Gilchrist RB. Role of oocyte-secreted growth differentiation factor 9 in the regulation of mouse cumulus expansion. Endocrinology. 2005;146:2798-806. PMID: 15761035 DOI: 10.1210/en.2005-0098
Medline Crossref

Elvin JA, Yan C, Matzuk MM. Growth differentiation factor-9 stimulates progesterone synthesis in granulosa cells via a prostaglandin E2/EP2 receptor pathway. Proc Natl Acad Sci U S A. 2000;97:10288-93. PMID: 10944203 DOI: 10.1073/pnas.180295197
Medline Crossref

Elvin JA, Yan C, Wang P, Nishimori K, Matzuk MM. Molecular characterization of the follicle defects in the growth differentiation factor 9-deficient ovary. Mol Endocrinol. 1999; 13: 1018-34. PMID: 10379899 DOI: 10.1210/mend.13.6.0309
Medline Crossref

Ferraretti AP, La Marca A, Fauser BC, Tarlatzis B, Nargund G, Gianaroli L; ESHRE working group on Poor Ovarian Response Definition. ESHRE consensus on the definition of ‘poor response’ to ovarian stimulation for in vitro fertilization: the Bologna criteria. Hum Reprod. 2011;26:1616-24. PMID: 21505041 DOI: 10.1093/humrep/der092
Medline Crossref