Skip to main content

Novel bone morphogenetic protein 15 (BMP15) gene variants implicated in premature ovarian insufficiency

Abstract

Background

Bone morphogenetic protein 15 (BMP15) is expressed in oocytes and plays a crucial role in the reproduction of mono-ovulating species. In humans, BMP15 gene mutations lead to imperfect protein function and premature ovarian insufficiency. Here we investigated the BMP15 gene variants in a population of Iranian women with premature ovarian insufficiency. We conducted predictive bioinformatics analysis to further study the outcomes of BMP15 gene alterations.

Methods

Twenty-four well-diagnosed premature ovarian insufficiency cases with normal karyotype participated in this study. The entire coding sequence and exon-intron junctions of the BMP15 gene were analyzed by direct sequencing. In-silico analysis was applied using various pipelines integrated into the Ensembl Variant Effect Predictor online tool. The clinical interpretation was performed based on the approved guidelines.

Results

By gene screening of BMP15, we discovered p.N103K, p.A180T, and p.M184T heterozygous variants in 3 unrelated patients. The p.N103K and p.M184T were not annotated on gnomAD, 1000 Genome and/or dbSNP. These mutations were not identified in 800 Iranians whole-exome sequencing that is recorded on Iranom database. We identified the p.N103K variant in a patient with secondary amenorrhea at the age of 17, elevated FSH and atrophic ovaries. The p.M184T was detected in a sporadic case with atrophic ovaries and very high FSH who developed secondary amenorrhea at the age of 31.

Conclusions

Here we newly identified p.N103K and p.M184T mutation in the BMP15 gene associated with idiopathic premature ovarian insufficiency. Both mutations have occurred in the prodomain region of protein. Despite prodomain cleavage through dimerization, it is actively involved in the mature protein function. Further studies elucidating the roles of prodomain would lead to a better understanding of the disease pathogenesis.

Background

Premature Ovarian Insufficiency (POI) also called premature ovarian failure (POF) is considered by loss of ovarian activity before the age of 40. It can be distinguished as primary amenorrhea with delayed menarche, secondary amenorrhea and oligomenorrhea over 4 months. POI is diagnosed by an increase in follicle-stimulating hormone (FSH) of higher than 25mIU/ml, repeated twice over 4 weeks [1]. Antimullerian hormone (AMH) is produced by the granulosa cells of growing follicles and is of considerable value in the timely detection of POI. Ovarian reserve can be detected by transvaginal ovarian ultrasound and estimated by AMH levels in POI cases [2]. The prevalence of POI was mainly considered 1 -1.5%. However, a recent meta-analysis reported a global prevalence of 3.7% with a higher rate in medium and low developed countries [3].

In the aetiology of POI, a wide range of explanations have been considered including involvement of infectious, autoimmune, iatrogenic or genetic factors. Among them, genetic aberrations account for approximately 20–25% of cases. FMR1 permutation is the most common genetic abnormality in POI, followed by bone morphogenetic protein 15 (BMP15) gene defects [4]. BMP15 gene is located on Xp11.2 and encodes a 392- amino acid protein of transforming growth factor-β superfamily (TGF-β). The protein is composed of the signal peptide, prodomain and transforming growth factor β- like (TGFβ-like) domain [5]. BMP15 is expressed in oocytes and acts synergically with growth differentiation factor 9 (GDF9) towards the regulation of folliculogenesis [6]. Using a mouse model, it has been shown that BMP15 promote follicular development along with FSH [7]. In vitro study revealed that the simultaneous addition of GDF9 and BMP15 to the culture leads to improvement of the human primordial follicles activation [8].

BMP15 peptides either form homodimers or bind to GDF9 and make heterodimers called cumulin. Dimerization occurs via the TGFβ-like domain. A furin-like protease cleaves prodomain upon protein maturation, however, prodomain remains associated with secreted mature dimer. BMP15 and cumulin control the migration and proliferation of primordial germ cells through separate signal pathways [9, 10].

Pathogenic mutations can disrupt the transcription/translation of BMP15 protein or weaken the interaction with GDF9 [11]. Preliminary studies on knockout models have shown that GDF-9 deficiency led to impaired folliculogenesis and infertility, while a null mutation in the BMP-15 gene caused only reduced fertility in female mice [12]. Subsequent studies suggested differences in BMP15 function between mice and other mammals. A study aiming to identify factors regulating ovulation rate in sheep found that heterozygous carriers of BMP15 mutation had increased ovulation while homozygous were infertile [13]. Additional evidence from the porcine knockdown model confirmed that the role of BMP15 in single-ovulatory species is more pronounced than in multi-ovulatory species [14]. This difference is explained in the balance of GDF9 and BMP15. The dominance of GDF9 signaling leads to multiple ovulations and high fertility, whereas the production of BMP15 modulate pri-mordial germ cell sensitivity to gonadotropins and leads to single ovulation, low fertility and appropriate ovarian reserve [15].

Considering the importance of BMP15 in the pathogenesis of POI, we aimed to investigate the related mutations in a population of affected Iranian women. We further conducted predictive bioinformatics analysis in a comprehensive manner, to evaluate BMP15 gene alterations linked to POI.

Methods

Patients

This study was approved by the Ethics Committee of Tarbiat Modares University, Tehran, Iran (IR.MODARES.REC.1399.006). The clinical interpretation was performed based on the European society of human reproduction and embryology (ESHRE) guidelines [2]. All participants delivered written informed consent before participating in the study. Physical examination, blood testing and sonography were performed as part of routine clinical workup by the experienced gynecologists. Twenty four normal-karyotype women with idiopathic secondary amenorrhea before the age of 40 were entered into our study. Based on the genetic counselling, patients were further allocated into sporadic and familial groups. Having at least two serum FSH levels higher than 25mIU/ml was the main inclusion criteria. The exclusion criteria comprised any of the following: chromosomal aberrations including FMR1 premutation, personal or family history of autoimmune diseases, ovarian surgery, radiotherapy or chemotherapy. To rule out autoimmune disease, related serological tests containing thyroid peroxidase antibody, anti-tissue transglutaminase IgG and thyroglobulin antibody were checked.

DNA extraction and PCR-Sequencing

Genomic DNA was extracted from whole blood samples using the salting-out method. The purity and concentration of the extracted DNA were measured by spectrophotometry. Specific primers were designed for two exons of BMP15 using Oligo explorer and Primer Express software as follows: exon1F;GTTGTGGAGCCAGGATGCAG, exon1R;AAGCCTGACAGTAAACCCACC, exon2F;AATTTTAGGGCTGATTATAGC, exon2R;TTGGTACAGGATTACTTGCAG.

To perform the PCR reaction, 2X PCR master mix (Amplicon, Pishgam, Iran) containing an ultimate concentration of 1.5 mM MgCl2 was used in a final volume of 50 µl. The PCR program included 5 min initial priming at 95 °C and then 30 cycles with 30 s at 95 °C, 30 s at 59 °C, 30 s at 72 °C with a final extension of 72 °C for 10 min. The PCR products were electrophoresed on 1.5% agarose gel along with the 100 base pairs DNA marker. Bi-directional Sanger sequencing was performed on all samples across two coding exons and flanking intronic sequences. Sanger sequencing results were interpreted using the Chromas software version 2.01.

In silico analysis

The probable pathogenic effect of the amino acid changes was detected by SIFT, PolyPhen, etc., which were integrated into the Ensembl Variant Effect Predictor (VEP) https://asia.ensembl.org/Homo_sapiens/Tools/VEP?db=core. The prevalence of the variants was checked on the Iranome website http://www.iranome.ir/. The Iranome browser represents whole-exome sequencing data of 800 healthy individuals from eight Iranian major ethnic groups with the approximate 1:1 Female/Male ratio. All samples were selected from individuals who were > 30 years old to reduce the bias of late-onset Mendelian disorders. The mean age at blood draw was 50.61 ± 9.33, and the age range 30–84 years old in the whole project [16].

Structural modeling of BMP15 wild-type and missense variants were computed by the SWISS-MODEL servers. The protein models were assessed using the Ramachandran plot from the MolProbity program.

Statistical analyses

The collected clinical data were statistically analyzed using SPSS version 24 with a significance level chosen at 5%. Correlation studies of the quantitative and qualitative data were performed by ANOVA and Chi-square tests, respectively.

Results

Clinical status of patients

All of the patients had spontaneously entered menarche and none of the patients was exhibited primary amenorrhea. The mean age at menarche was calculated at 13.12 ± 1.56 years old. The mean age at POI diagnosis was 26.84 ± 9.07 years old ranged from 12 to 39. Median and mean levels of serum FSH were 40.8 and 75.92 mIU/ml, respectively. The analyses of AMH serum levels were available for 19 out of 24 patients and the majority of them were below 0.02ng/ml. The sonographic evaluation revealed that 79.16% of the patients had atrophic ovaries and 20.83% had normal ovarian pattern. Patients’ interviews and genetic counselling disclosed that 6 out of 24 cases had a family history of menstrual irregularities, infertility and/or secondary amenorrhea.

BMP15 genotyping

As shown in Fig. 1, by the PCR-sequencing analysis of BMP15, we discovered p.N103K, p.A180T, and p.M184T heterozygous variants in 3 unrelated patients. This is the first report of the link between p.N103K and p.M184T to the POI pathogenesis.

Fig. 1
figure 1

Electropherogram analysis of mutations and SNP alleles in unrelated POI patients

The patient carrier of p.N103K variant entered menarche at the age of 16 and secondary amenorrhea happened when she was at the age of 17. As specified in Table 1, elevated FSH levels and a low normal range of AMH were detected in this patient. The sonographic evaluation exhibited an atrophic ovaries pattern. She was born from a consanguineous background without any indication of POI or related conditions in her family members.

Table 1 Identified variants and clinical characteristics of the patients according to the ESHRE guideline

p.M184T was detected in a sporadic POI patient who was prematurely menopaused at the age of 31. High FSH levels of 100mIU/ml and very low levels of AMH were detected by hormonal profiling. Ultrasound examination showed an atrophic pattern in both ovaries. Detailed family history revealed that the parents were not relatives and there was no history of POI and/or related disorders among her relatives.

We identified p.A180T in a 42-years-old patient who was a descendant of consanguineous marriage with a familial history of POI. She had been diagnosed at age of 39 following a history of oligomenorrhea for more than 18 months. The hormonal assessment showed high FSH and very low AMH levels. The sonographic evaluation indicated normal ovarian status (Table 1).

We also detected two known variants of BMP15; c.-9 C > G and p.N103S [17]. Previous studies from Syria and India reported the p.N103S, in both cases and controls [18, 19]. This variant is mentioned in the ClinVar database as a benign variant causing ovarian dysgenesis. The annotation of its corresponding allele, c.308 A > G showed an allele frequency of 0.057 and 0.028 in gnomAD and 1000 Genome, respectively. Iranome browser specified a high c.308 A > G allele frequency of 0.075 in the Iranian population. In our study 1 out of 24 patients (1/48 alleles) was the carrier of this SNP variant.

rs3810682, c.-9 C > G

The 5`-UTR variant, c.-9 C > G, was detected in 29.16% of the patients, all in heterozygous status including 2 compound heterozygous with p.A180T and p.N103S. Statistical analysis was performed to calculate the correlation between c.-9 C > G SNP and clinicopathological parameters of the patients. Table 2 shows a statistically significant association between c.-9 C > G SNP and patients FSH levels (p value = 0.02).

Table 2 Correlation between the − 9G polymorphism and clinicopathological parameters of the patients,

In-silico analysis

Further in-silico analysis was conducted to search the effects of the newly identified p.N103K and p.M184T on protein biogenesis and function. Table 3 lists the in-silico predicted values of these mutations in comparison to previously reported nonsynonymous missense alterations of BMP15. The indicated data were extracted from the Ensembl VEP online tool.

Table 3 Variants annotation, allele frequencies and in-silico prediction of nonsynonymous missense alterations of BMP15

The newly identified p.N103K was not annotated on gnomAD, 1000 Genome or dbSNP. Investigation on the Iranome database revealed that p.N103K did not occur in 800 healthy Iranians. The SIFT and PolyPhen results showed scores of 0.35 and 0.006, respectively. In the same way, the p.M184T new mutation was not annotated in databases. Bioinformatics analysis showed scores of 0.13 and 0.003 by the SIFT and PolyPhen prediction tools, respectively.

Figure 2 shows the conserved score of the M184 and N103 residues among different species indicating a higher score for M184. Both mutations are located in the secondary structure elements of the prodomain (Fig. 3B), which are essential for regulating mature protein production, secretion, and activity. M184 and N103 are present in β6 and β1` sheet of the peptide structure, respectively [26]. As shown in Fig. 3 A, M184 and A180 are situated in the same β-sheet of the peptide structure.

Fig. 2
figure 2

The conservation scores of the first 200 amino acids of BMP15 protein among different species. Blue and red lines show the two side of the spectrum from amino acid insertions or differences with the lowest degree of conservation to sequence identity among different species, respectively. Black arrows indicate the score of N103 and M184 residues, respectively. Analyzed on https://www.ibi.vu.nl/programs/pralinewww/

Fig. 3
figure 3

In A, the location of N103, A180 and M184 are shown. B indicates the location of reported nonsynonymous missense alterations of BMP15. Most of the variants are located on prodomain except for p.R329C, p.R329H and p.C357Y which are located on the mature TGFβ-like domain

Modelling and Mapping of the wild-type and missense variants

Figure 4 A shows molecular modeling of wild-type M184 and N103 as compared with mutant T184 and K103, respectively. The results confirmed the localization of amino acids on the protein β-sheet structures. Figure 4B provides Ramachandran plots exhibiting Φ/Ψ backbone dihedral angles of the wild- type and mutant proteins. The wild-type showed 88.4% and 4.03% residues in the favored and outliers regions of the plot. The generated model for M184T was found to have fewer residues in the favored region (87.6%). The Rotamer Outliers was also increased for the M184T (0.79%) vs. wild- type (0.16%) which indicates a decrease in protein stability. The mutant N103K displayed the same configuration as that of the wild-type structure.

Fig. 4
figure 4

 A Modeled BMP15 amino acid replacement on the three-dimensional structure. B General (no proline or glycine) Ramachandran plots of the BMP15 wild-type, M184T and N103K proteins which were provided by SWISS-MODEL online server. The different colored areas correspond to favored, allowed, generously allowed and sterically disallowed conformations

Segregation analysis of the p.A180T

As mentioned earlier, p.A180T was found in a patient with a family history of POI. The c.538G˃A missense variant located in exon 2 of BMP15 is presented in databases as the rs104894767. As illustrated in the family pedigree, two older sisters of the proband were also reached menopause before the age of 40, however, their mother entered physiological menopause at the age of 55. By Sanger sequencing, the c.538G˃A variant was analyzed in the female family members. The obtained results indicated a heterozygote alteration in the 64-years-old mother and one of the affected sisters who was menopaused at 35. The other POI sister with menopause at 40 and the 43-years-old healthy sister with regular menses did not carry the p.A180T variant (Fig. 5).

Fig. 5
figure 5

Family pedigree shows the proband (arrow) and the rest of the affected family members with POI presented in colored circles. Genotyping of BMP15 is indicated below the symbols, with the A180T heterozygote variant

Discussion

In this study, we provided the molecular analysis of the BMP15 gene in a series of unrelated Iranian women with idiopathic POI. Our results revealed two novel variants; p.N103K and p.M184T as well as one known substitution; p.A180T.

With the development of next-generation sequencing, several new genes have been linked to the disease pathogenesis which supports the concept of POI heterogeneity. Investigation of copy number variations led to the identification of microdeletions/microduplications on multiple chromosomal regions. Using high-resolution SNP arrays a 4 MB deletion in Xp11.23-p11.22 was discovered in a patient with secondary amenorrhea at the age of 34. This deletion consists of 65 genes, including the entire BMP15 [28]. A more recent study reported the deletion of both paternal and maternal copies of BMP15 in a 16-year-old patient with primary amenorrhea and FSH levels of 141mIU/ml. The deletion expanded 36 KB and contained the entire BMP15 in the Xp11.22 region. In addition, the study showed intragenic duplication of the tumor protein P63 (TP63) gene in two sisters who were also the carrier of BMP15 p.Y235C [29]. Remarkably, p.Y235C is the first discovered BMP15 mutation that was reported in 2004 in the heterozygous Italian sisters [25]. The same research group further identified p.R68W and p.A180T alterations in non-familial POI patients in 2006 [21]. In 2009 they repeated the BMP15 gene screening in 300 patients and found p.R138H and p.L148P in the prodomain region. They also reported p.S5R mutation which is localized in the signal peptide of the nascent protein. By conducting a functional in vitro assay in the human granulosa cell line, they disclosed that p.R68W, p.L148P and p.R138H cause the reduction of mature BMP15 protein. While p.A180T or p.S5R displayed no deleterious impact on protein secretion or function [20]. However, in 2017, Patiño et al. estimated four-fold lower activity of p.A180T compared to wild-type BMP15. They evaluated the expression and activity of 10 BMP15 variants and showed that these nucleotide changes could reduce mature peptide production, activity or synergy with GDF9 [26].

The obtained results regarding p.A180T remained conflicting, as in several studies it was only detected in patients [11, 19], while in others, both in patients and controls [22, 30]. Based on the Iranome database, the allele frequency of the p.A180T in the Iranian healthy population is 0.006. In our study, we found the p.A180T variant in a POI patient who was born from a consanguineous marriage with a family history of POI. Segregation analysis confirmed that the p.A180T could not be the causative mutation in this family. Modifier genes and SNPs should be considered as the influential factors in the genotype/phenotype correlation.

We further identified p.N103K missense heterozygote transition. This new BMP15 mutation was detected in the offspring of first-cousin parents without any family history of POI. Our next newly identified missense heterozygote substitution; p.M184T was identified in a sporadic POI case exhibiting high FSH levels. Since p.N103K and p.M184T are not annotated in the reference databases it can be assumed that these new variants are possibly pathogenic. Their placement in the functional regions of the protein supports this conclusion.

Figure 3B indicates that p.N103K and p.M184T are localized on the prodomain region of the protein. The prodomain is cleaved during the maturation process and remains non-covalently bound to the homo or heterodimers upon secretion. The folding and dimerization of the mature protein are regulated by the prodomain [31]. It has been shown that prodomain mutations can prevent its binding to mature dimers leading to decreased BMP15 activity [32]. As mentioned previously, BMP15 reduction is correlated to increased ovulation rate and the chance of dizygotic twins’ pregnancy. This increase in ovulation can eventually lead to ovarian depletion and POI [33].

It should be noted that sheep BMP15 mutations have often been occurred in the mature protein, while most human BMP15 mutations are prodomain variants [26]. Nonsynonymous missense variants of BMP15 that are listed in Table 3 are prodomain alterations except for; p.R329C, p.R329H and p.C357Y which are located on the mature TGFβ-like domain.

In 2018, Zhang et al. reported homozygous p.C357Y mutation in the TGFβ-like domain. To our knowledge, this is the sole described homozygous nonsynonymous missense variant of the BMP15. The patient`s heterozygous mother was 48 years old with normal menstrual cycles. They assumed that the mutations in the TGFβ-like domain should be homozygous to cause the disease while the heterozygous mutations of the prodomain are sufficient for defective protein maturation and POI [27]. However, Patiño et al. reported another mutation in TGFβ-like domain; p.R329H, in a heterozygous POI patient [26]. The p.R329C mutation of the mature subunit was reported earlier from China. The patient was a 37-year-old woman who developed secondary amenorrhea at the age of 27. The patient’s mother had irregular menses and menopause at the age of 40. The data on hetero or homozygosity of these mother and daughter is not available [24]. More recently, Rossetti et al. examined p.R329C mutation, functionally. The results showed that heterozygote p.R329C mutation led to disrupted co-localization with GDF9 and SMAD pathway activation [34].

We found p.N103K and p.M184T in the heterozygous patients, like almost all previously reported BMP15 mutations [23]. This feature was often concluded by haploinsufficiency or negative dominance effects. However, the discovery of intragenic duplication of the TP63 in the carriers of p.Y235C has added to the complexity of interpreting the BMP15 mutations.

Conclusions

In this study, we found two new variants in the BMP15 gene that have not been reported in databases. These mutations were not identified in 800 Iranians whole-exome sequencing results available on Iranom’s website.

The p.N103K and p.M184T are localized on the prodomain region of the protein. Due to the cleavage of prodomain during BMP15 maturation, bioinformatics analyzes are likely to underestimate the effects of the mutations. This is while prodomain is actively involved in the homo and heterodimers secretion and function. Further studies to elucidate the roles of the prodomain could be of great value in identifying the effects of multiple mutations reported in this region.

Availability of data and materials

The authors confirm that the data supporting the findings of this study are available within the article. Raw data of are available from the corresponding author on request.

References

  1. Nelson LM. Clinical practice. Primary ovarian insufficiency. N Engl J Med. 2009;360(6):606–14.

    CAS  Article  Google Scholar 

  2. Webber L, Davies M, Anderson R, Bartlett J, Braat D, Cartwright B, et al. ESHRE Guideline: management of women with premature ovarian insufficiency. Hum Reprod (Oxford England). 2016;31(5):926–37.

    CAS  Article  Google Scholar 

  3. Golezar S, Ramezani Tehrani F, Khazaei S, Ebadi A, Keshavarz Z. The global prevalence of primary ovarian insufficiency and early menopause: a meta-analysis. Climacteric: the journal of the International Menopause Society. 2019;22(4):403–11.

    CAS  Article  Google Scholar 

  4. Persani L, Rossetti R, Cacciatore C. Genes involved in human premature ovarian failure. J Mol Endocrinol. 2010;45(5):257–79.

    CAS  Article  Google Scholar 

  5. Dube JL, Wang P, Elvin J, Lyons KM, Celeste AJ, Matzuk MM. The Bone Morphogenetic Protein 15 Gene Is X-Linked and Expressed in Oocytes. Mol Endocrinol. 1998;12(12):1809–17.

    CAS  Article  Google Scholar 

  6. Abir R, Fisch B, Johnson MH. BMP15, fertility and the ovary. Reprod Biomed Online. 2014;29(5):525–6.

    Article  Google Scholar 

  7. Salehnia M, Pajokh M, Ghorbanmehr N. Short Term Organ Culture of Mouse Ovary in the Medium Supplemented with Bone Morphogenetic Protein 15 and Follicle Stimulating Hormone: A Morphological, Hormonal and Molecular Study. J Reprod infertility. 2016;17(4):199–207.

    Google Scholar 

  8. Kedem A, Fisch B, Garor R, Ben-Zaken A, Gizunterman T, Felz C, et al. Growth differentiating factor 9 (GDF9) and bone morphogenetic protein 15 both activate development of human primordial follicles in vitro, with seemingly more beneficial effects of GDF9. J Clin Endocrinol Metab. 2011;96(8):E1246-54.

    Article  Google Scholar 

  9. Moore RK, Shimasaki S. Molecular biology and physiological role of the oocyte factor, BMP-15. Mol Cell Endocrinol. 2005;234(1–2):67–73.

    CAS  Article  Google Scholar 

  10. Peng J, Li Q, Wigglesworth K, Rangarajan A, Kattamuri C, Peterson RT, et al. Growth differentiation factor 9:bone morphogenetic protein 15 heterodimers are potent regulators of ovarian functions. 2013;110(8):E776-E85.

  11. Laissue P, Christin-Maitre S, Touraine P, Kuttenn F, Ritvos O, Aittomaki K, et al. Mutations and sequence variants in GDF9 and BMP15 in patients with premature ovarian failure. Eur J Endocrinol. 2006;154(5):739–44.

    CAS  Article  Google Scholar 

  12. Yan C, Wang P, DeMayo J, DeMayo FJ, Elvin JA, Carino C, et al. Synergistic roles of bone morphogenetic protein 15 and growth differentiation factor 9 in ovarian function. Mol Endocrinol (Baltimore Md). 2001;15(6):854–66.

    CAS  Article  Google Scholar 

  13. McNatty KP, Juengel JL, Wilson T, Galloway SM, Davis GH. Genetic mutations influencing ovulation rate in sheep. Reprod Fertil Dev. 2001;13(7–8):549–55.

    CAS  Article  Google Scholar 

  14. Qin Y, Tang T, Li W, Liu Z, Yang X, Shi X, et al. Bone Morphogenetic Protein 15 Knockdown Inhibits Porcine Ovarian Follicular Development and Ovulation. Front Cell Dev Biology. 2019;7:286.

    Article  Google Scholar 

  15. McNatty KP, Moore LG, Hudson NL, Quirke LD, Lawrence SB, Reader K, et al. The oocyte and its role in regulating ovulation rate: a new paradigm in reproductive biology. Reprod (Cambridge England). 2004;128(4):379–86.

    CAS  Article  Google Scholar 

  16. Fattahi Z, Beheshtian M, Mohseni M, Poustchi H, Sellars E, Nezhadi SH, et al. Iranome: A catalog of genomic variations in the Iranian population. Hum Mutat. 2019;40(11):1968–84.

    CAS  Article  Google Scholar 

  17. Abir R, Fisch B. Invited commentary: a single nucleotide polymorphism in BMP15 is associated with high response to controlled ovarian hyperstimulation. Reprod Biomed Online. 2011;23(1):77–80.

    CAS  Article  Google Scholar 

  18. Al-ajoury R, Kassem E, Al-halabi B, Moassess F, Al-achkar W. Investigation of some genetic variations in BMP15 accompanied with premature ovarian failure (POF) in Syrian women. Middle East Fertility Society Journal. 2015;20(2):91–6.

    Article  Google Scholar 

  19. Dixit H, Rao LK, Padmalatha VV, Kanakavalli M, Deenadayal M, Gupta N, et al. Missense mutations in the BMP15 gene are associated with ovarian failure. Hum Genet. 2006;119(4):408–15.

    CAS  Article  Google Scholar 

  20. Rossetti R, Di Pasquale E, Marozzi A, Bione S, Toniolo D, Grammatico P, et al. BMP15 mutations associated with primary ovarian insufficiency cause a defective production of bioactive protein. Hum Mutat. 2009;30(5):804–10.

    CAS  Article  Google Scholar 

  21. Di Pasquale E, Rossetti R, Marozzi A, Bodega B, Borgato S, Cavallo L, et al. Identification of new variants of human BMP15 gene in a large cohort of women with premature ovarian failure. J Clin Endocrinol Metab. 2006;91(5):1976–9.

    Article  Google Scholar 

  22. Tiotiu D, Alvaro Mercadal B, Imbert R, Verbist J, Demeestere I, De Leener A, et al. Variants of the BMP15 gene in a cohort of patients with premature ovarian failure. Hum Reprod (Oxford England). 2010;25(6):1581–7.

    CAS  Article  Google Scholar 

  23. Ferrarini E, De Marco G, Orsolini F, Gianetti E, Benelli E, Fruzzetti F, et al. Characterization of a novel mutation V136L in bone morphogenetic protein 15 identified in a woman affected by POI. J ovarian Res. 2021;14(1):85.

    CAS  Article  Google Scholar 

  24. Wang B, Wen Q, Ni F, Zhou S, Wang J, Cao Y, et al. Analyses of growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) mutation in Chinese women with premature ovarian failure. Clin Endocrinol. 2010;72(1):135–6.

    CAS  Article  Google Scholar 

  25. 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(1):106–11.

    Article  Google Scholar 

  26. Patiño LC, Walton KL, Mueller TD, Johnson KE, Stocker W, Richani D, et al. BMP15 Mutations Associated With Primary Ovarian Insufficiency Reduce Expression, Activity, or Synergy With GDF9. J Clin Endocrinol Metab. 2017;102(3):1009–19.

    PubMed  Google Scholar 

  27. Zhang W, Wang J, Wang X, Li L, Pan H, Chen B, et al. A novel homozygous mutation of bone morphogenetic protein 15 identified in a consanguineous marriage family with primary ovarian insufficiency. Reprod Biomed Online. 2018;36(2):206–9.

    CAS  Article  Google Scholar 

  28. McGuire MM, Bowden W, Engel NJ, Ahn HW, Kovanci E, Rajkovic A. Genomic analysis using high-resolution single-nucleotide polymorphism arrays reveals novel microdeletions associated with premature ovarian failure. Fertil Steril. 2011;95(5):1595–600.

    CAS  Article  Google Scholar 

  29. Bestetti I, Castronovo C, Sironi A, Caslini C, Sala C, Rossetti R, et al. High-resolution array-CGH analysis on 46,XX patients affected by early onset primary ovarian insufficiency discloses new genes involved in ovarian function. Hum Reprod (Oxford England). 2019;34(3):574–83.

    CAS  Article  Google Scholar 

  30. Ledig S, Röpke A, Haeusler G, Hinney B, Wieacker P. BMP15 mutations in XX gonadal dysgenesis and premature ovarian failure. Am J Obstet Gynecol. 2008;198(1):84.e1-5.

    Article  Google Scholar 

  31. Belli M, Shimasaki S. Molecular Aspects and Clinical Relevance of GDF9 and BMP15 in Ovarian Function. Vitam Horm. 2018;107:317–48.

    CAS  Article  Google Scholar 

  32. Inagaki K, Shimasaki S. Impaired production of BMP-15 and GDF-9 mature proteins derived from proproteins WITH mutations in the proregion. Mol Cell Endocrinol. 2010;328(1–2):1–7.

    CAS  Article  Google Scholar 

  33. Otsuka F, McTavish KJ, Shimasaki S. Integral role of GDF-9 and BMP-15 in ovarian function. Mol Reprod Dev. 2011;78(1):9–21.

    CAS  Article  Google Scholar 

  34. Rossetti R, Ferrari I, Bestetti I, Moleri S, Brancati F, Petrone L, et al. Fundamental role of BMP15 in human ovarian folliculogenesis revealed by null and missense mutations associated with primary ovarian insufficiency. Hum Mutat. 2020;41(5):983–97.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

The author would like to thank the patients who participated in this survey.

Funding

Not applicable.

Author information

Affiliations

Authors

Contributions

The author’s contributions are as follows FA: Project development, Data collection, Data analysis, Manuscript writing. SS: Project development, Data collection, Data analysis, Manuscript writing/editing. LF: Clinical Data collection and management, Manuscript editing. SD: Clinical Data collection and management, Manuscript editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Shirin Shahbazi.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. This study was approved by the Ethics Committee of Tarbiat Modares University, Tehran, Iran (IR.MODARES.REC.1399.006). https://ethics.research.ac.ir/EthicsProposalViewEn.php?id=130467.

Consent for publication

Patients signed informed consent regarding publishing their data.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Afkhami, F., Shahbazi, S., Farzadi, L. et al. Novel bone morphogenetic protein 15 (BMP15) gene variants implicated in premature ovarian insufficiency. Reprod Biol Endocrinol 20, 42 (2022). https://doi.org/10.1186/s12958-022-00913-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12958-022-00913-6

Keywords

  • POI
  • BMP15
  • Gene mutation
  • In silico analysis