Skip to main content

Vascular endothelial growth factor (VEGFA) gene variation in polycystic ovary syndrome in a Tunisian women population

Abstract

Background

Polycystic ovary syndrome (PCOS) is characterized by the growth of a number of small cysts on the ovaries which leads to sex hormonal imbalance. Women who are affected by this syndrome suffer from irregular menstrual cycles, decline in their fertility, excessive hair growth, obesity, acne and most importantly cardiac function problems. The vascular endothelial growth factor (VEGF) plays a pivotal role in tissue vascularization in general and in the pathogenesis of many diseases. The PCOS was found to be associated with high expression levels of VEGF. In women who undergo assisted reproductive procedures (ART), VEGF was found to be a key mediator of other factors to control ovary angiogenesis. Here, we set out to examine the association of VEGFA gene polymorphism with PCOS and its components in a population of Tunisia women to enhance our understanding of the genetic background leading angiogenesis and vascularization abnormalities in PCOS.

Methods

The association of VEGFA gene with PCOS and its components was examined in a cohort of 268 women from Tunisia involving 118 PCOS patients and 150 controls. VEGFA gene variations were assessed through the analysis of the following SNPs rs699947 (A/C), rs833061 (C/T), rs1570360 (G/A), rs833068 (G/A), rs3025020 (C/T), and rs3025039 (C/T). The linkage disequilibrium between SNPs was assessed using HAPLOVIEW software while combination of SNPs into haplotypes in the population and the reconstruction of the cladogram were carried-out by PHASE and ARLEQUIN programs, respectively. Genetic association and genotype-phenotype correlations were calculated by logistic regression and non-parametric tests (Kruskall-Wallis and Mann–Whitney tests), respectively, using StatView program.

Results

We observed 10 haplotypes in our studied cohort whereH1 (ACGG), H2 (ACAG), H7 (CTGG) and H8 (CTGA) were the most frequent. We observed the association of the genotype CT of the SNP rs30225039 with PCOS phenotype (P = 0.03; OR 95 % CI = 2.05 [1.07–3.90]) and a trend for correlation of the pair of haplotypes H2/H2 with prolactin levels in plasma (P = 0.077; 193.5 ± 94.3 vs 45.7 ± 7.2). These data are consistent with literature and highlight one more time the role of vascularization in the pathogeny of PCOS.

Conclusions

LD pattern in VEGF locus showed a similar LD pattern between the Tunisian population and the CEU. More haplotypes in the Tunisian population than in CEU was observed (22 haplotypes vs 16 haplotypes) suggesting higher recombination rate in Tunisians. The study showed that there was any advantage of using haplotypes compared with SNPs taken alone.

Background

The polycystic ovary syndrome (PCOS) is a disease affecting 6 – 12 % of women at reproductive age and characterized by abnormalities in relation to reproduction, hyperandrogenism as well as metabolic abnormalities like insulin resistance, obesity, impaired fasting glucose (IGF) and metabolic syndrome (MetS) [1, 2].

Knowing the role of blood flow and vascular pattern of an organ in the organ’s morphology and function, previous studies focused on the comparison of blood flow and vascularization between women with PCOS and healthy women showing larger and more vascularized ovaries in PCOS women compared to the control [3]. Moreover, abnormalities in ovarian angiogenesis were involved in the ovarian hyperstimulation syndrome (OHSS) and other disorders of anovulation, subfertility and pathogenic conditions as endometriosis [4]. In this context, the vascular endothelial growth factor (VEGF) system, composed of ligands and receptors plays a pivotal role in tissue vascularization and endothelial cell growth [5]. It was also reported to induce cell proliferation, promote cell migration, inhibit apoptosis and induce permeabilization of blood vessels. It is a major component of the family of angiogenic factors, which includes placental growth factor, angiopoetin, basic fibroblast growth factor and the VEGF family (A, B, C, D and E). VEGFA in particular is a protein encoded by a gene (Gene ID: 7422) located in chromosome 6 and extended on more than 16 kb (GRCh37/hg19; chr6: 43737946–43754224). The gene involves 8 exons and encodes for many VEGFA isoforms as VEGF-165 and VEGF-145.

Previous studies reported the association of VEGFA gene with PCOS through SNPs +9812, +13553, −2578 (rs699947), −460 (rs833061) and +405 [6, 7] and their haplotypes in populations of different origins [6–8]. Moreover, SNPs +405 and −460 were reported as having a trend towards higher insulin resistance [9]. It was suggested the role of some of these SNPs on the regulation of VEGFA gene expression and protein production and secretion [10].

In this context, we focused on the study of the association of VEGFA gene with PCOS and its components in a population of women from Tunisia well characterized at the phenotypic level. The aim was to better characterize the genetic background responsible for angiogenesis and vascularization abnormalities in PCOS. This will be achieved through testing the association of VEGFA gene in a population where it was previously highlighted the need for haplotyping approach in order to understand better the genetic association in populations with mosaic anthropogenetic components [11, 12].

Subjects and methods

  1. 1.

    Subjects

    Following an informed consent, a total of 268 unrelated Tunisian women, consisting of 118 PCOS patients and 150 controls. PCOS patients group (mean age: 29.8 ± 0.4 years) were recruited from the outpatient Obstetrics and Gynecology Department, CHU Farhat Hached (Sousse) in Central Tunisia (Table 2). The diagnosis of PCOS was based on the 2003 Rotterdam Criteria [13]. PCOS was then attested when two of the following three conditions are observed: anovulation, hyperandrogenism, and the presence of polycystic ovaries on ultrasound examination. Exclusion criteria were androgen-producing tumors, 21-hydroxylase-deficiency, non-classical adrenal hyperplasia, hyperprolactinemia, active thyroid disease, and Cushing’s syndrome.

    Control group consisted of 150 healthy women (mean age: 33.5 ± 0.5 years), with regular menstrual cycles and no evidence of hirsutism, acne, alopecia, or endocrinopathies. None of the controls was on hormonal therapy (including oral contraceptives) for the previous three months or longer, and none of PCOS women or controls was on medication known to affect carbohydrate metabolism or endocrine parameters, at least for the last three months before inclusion in the study. Demographic data and history of hypertension, diabetes, and hypercholesterolemia were recorded for all subjects. Obesity was defined as body-mass index (BMI) ≥30 kg/m2. The protocol of the study was approved by the local ethics committee, and written informed consent was obtained from all subjects.

  2. 2.

    SNPs genotyping

    Total genomic DNA was isolated from peripheral blood lymphocytes by the salting-out method. Six SNPs were selected because of their association with PCOS, angiogenesis and other metabolic components in literature: rs699947 (A/C), rs833061 (C/T), rs1570360 (G/A), rs833068 (G/A), rs3025020 (C/T), andrs3025039 (C/T). The 6 SNPs were genotyped by the allelic discrimination TaqMan SNP Genotyping Assays (Applied Biosystems; Foster City, CA). The PCR primers and TaqManprobes were directly available from Applied Biosystem (predesigned assays) and analyses were done according to the manufacturer’s protocol on a StepOne plus apparatus (Applied Biosystems) (Table 1).

    Table 1 SNPs probe sequence

Statistical analysis

Data were expressed as mean ± SE (continuous variables) or as percent of total (categorical variables), and intergroup significance was assessed by χ2 test. Allele frequencies were calculated by gene-counting method; each SNP was tested for Hardy–Weinberg equilibrium using Haploview 4.2 (http://www.broadinstitute.org/haploview). Linkage disequilibrium (LD) was assessed using Haploview 4.2. Genetic association and genotype-phenotype correlation were assessed by logistic regression and non-parametric tests (Mann–Whitney and Kruskall-Walli tests) using StatView program.

Calculation of the power of the study was assessed using SAMPSIZE software (http://sampsize.sourceforge.net).

Results and discussions

Phenotype data

Women with PCOS have higher BMI than controls (P < 0.0001; 28.4 ± 0.7 vs 23.1 ± 0.2 kg/m2), higher random glycemia (P <0.0001, 7.9 ± 0.2 vs 4.5 ± 0.1 mmol/L), fasting insulin (P = 0.01, 15.7 ± 1.2vs 7.7 ± 0.4 mU/L), levels of total testosterone (P < 0.0001, 2.9 ± 0.2 vs 1.0 ± 0.1 noml/L) and of plasma triglycerides (TG) (P = 0.0001, 1.7 ± 0.1 vs 1.0 ± 0.1). Prolactin levels in plasma were lower in women with PCOS compared with controls (P = 0.0001, 73.1 ± 11.7 vs 148.8 ± 9.4) (Table 2).

Table 2 Phenotypic features of our women cohort used in the study

Genotype data

All SNPs were in Hardy-Weinberg equilibrium. Allele frequencies of SNPs in PCOS and controls were respectively0.56 and 0.56 for rs699947, 0.51 and 0.54 for rs833061, 0.34 and 0.31 for rs1570360, 0.37 and 0.36 for rs833068, 0.25 and 0.22 for rs3025020 and 0.87 and 0.91 for rs3025039 (Table 3).

Table 3 Allele frequencies of SNPs in PCOS and controls

The study of LD pattern (assessed by D’ index) in the region delimited by these SNPs showed relative high LD between SNPs rs699947, rs833061, rs1570360 and rs833068 (Fig. 1) allowing then reconstruction of haplotypes with these 4 polymorphisms. We observed 10 haplotypes among which 4 haplotypes were the most frequent: H1 (ACGG)(13.4 %), H2 (ACAG) (29.3 %), H7 (CTGG) (19.6 %) and H8 (CTGA) (32.3 %)(Table 4).

Fig. 1
figure 1

Linkage disequilibrium (LD) pattern in the locus of VEGFA gene delimited by SNPs rs699947 and rs3025039. Numbers in the squares indicate D’ index (level of LD) between the corresponding SNPs

Table 4 List of haplotypes resulting from the combination of SNPs rs699947, rs833061, rs1570360 and rs833068 in the Tunisian population

Reconstruction of the cladogram in the population allowed observing that H2 arose from H1 after the mutation of the SNP r1570360 while H7 appeared after the mutation of rs1570360, rs833061 and rs699947. Finally, H8 arose after the mutation of rs1570360, rs833061, rs699947 and rs833068 (Fig. 2).

Fig. 2
figure 2

Cladistic representation of VEGFA haplotypes in the studied population. SNPs mutations responsible for the transition between haplotypes are indicated in the cladogram. H1 is the ancestral haplotype

Logistic regression showed the association of the genotype CT of rs30225039 with PCOS phenotype (P = 0.03; OR 95 % CI = 2.05 [1.07 – 3.90]) (Table 5) while non-parametric tests allowed observing a trend for correlation of the pair of haplotypes H2/H2 to prolactin levels in plasma (P = 0.077; 193.5 ± 94.3 vs 45.7 ± 7.2). According to the calculation, this study had a power of 58.38 %. PCOS is a disease that affects women at reproductive age and characterized by endocrine disorders as menses disturbance and subfertility and by metabolic abnormalities as insulin resistance, and glucose intolerance [14, 15]. Angiogenesis process has also been reported as probably implied in disorders of anovulation, subfertility and other pathogenic conditions as endometriosis. Since VEGF system plays a pivotal role in tissue vascularization, VEGFA gene is considered as candidate in reproductive abnormalities as PCOS. Many studies have been carried-out and concluded in a part of them that VEGFA levels in serum or follicular fluid are higher in women with PCOS [16–19]. In another hand, previous studies in Korean and Turkish populations reported the association of VEGFA gene with PCOS through many SNPs and one haplotype [6–8]. In that context, our investigation focused on the study of VEGFA gene association with PCOS in a Tunisian population well characterized at the phenotypic level, through 6 SNPs and their combination into haplotypes in the population. We observed the association of the gene (SNP rs30225039) with PCOS phenotype and a trend of correlation (pair of haplotypes H2/H2) with prolactin levels in plasma. Prolactin is secreted mainly by the pituitary gland and is implied in many processes as lactation and also tissue vascularization [20]. Indeed, many studies showed the role of prolactin as a pro-angiogenic factor while its enzymatically cleaved 16 kDa N-terminal fragment has a well-defined anti-angiogenic effect [21]. Moreover, it was shown the effect of prolactin on the induction of VEGF expression [22–24]. To the best of our knowledge, our study is the first work that reports the correlation, even with a trend, of VEGFA gene variations with prolactin plasma levels. These findings suggest possible regulation and/or feedback control of prolactin secretion pathway through VEGF and the role of VEGFA gene variation in this process. Deeper genetic investigations taking into account more SNPs and larger locus are necessary to demystify these molecular pathways.

Table 5 Distribution of VEGF genotype in PCOS cases and control women

Differences of VEGF plasma levels between PCOS and healthy women were observed in some populations but not in others and this could be due to genetic differences between ethnic groups [16–19] Studying VEGFA gene association in complex disorders should take into account the gene variability in the population through the LD pattern assessment and SNPs combination into haplotypes. We reported in a previous study the role of haplotypes as powerful genetic markers for better characterization of gene association in the Tunisian population and observed Fat mass and obesity associated (FTO) gene association with PCOS through haplotypes but not through SNPs taken alone [12]. Indeed, we observed LD pattern in the FTO locus less in the Tunisian population compared with the CEU Caucasian one justifying so the ability of haplotypes for detecting genetic association in such conditions. In the present study, the investigation of the LD pattern in VEGF locus showed a similar LD pattern between the Tunisian population and the CEU one (Fig. 3). However, we found more haplotypes in the Tunisian population than in CEU (22 haplotypes vs 16 haplotypes) suggesting more recombination rate operated in Tunisians. Although we were not able to observe any advantage of using haplotypes compared with SNPs taken alone in this study, these results could suggest a more complex anthropogenetics context concerning VEGF locus or the need for a better characterization of this locus by using more density of SNPs at larger genomic region/scale and using larger cohorts of patients. Therefore, better management approaches can be identified [25]. Additionally, and despite the relatively reduced size of our cohort, it remains that the powering of the study was acceptable (58.38 %) attesting then for the reliability of our findings.

Fig. 3
figure 3

Comparison of linkage disequilibrium (LD) pattern (LD is assessed by r2 index) in the VEGF locus delimited between SNPs rs699947 and rs3025039, between the studied women population and the CEU reference population. Numbers in the squares indicates levels of LD as informed by r2 index

Conclusions

Deeper investigations of VEGFA gene encompassing larger regions in the locus and using more SNPs need to be envisioned in order to better characterize the gene association in the Tunisian population and to shed light on the genetic/genomic players involved and develop better diagnostic, management and therapeutic approaches.

References

  1. Ehrmann DA. Polycystic Ovary Syndrome. N Engl J Med. 2005;352(12):1223–36.

    Article  CAS  PubMed  Google Scholar 

  2. Klein J, Craven M, Vuguin PM. Polycystic Ovarian Syndrome. Adolesc Med State Art Rev. 2015;26(2):326–42.

    PubMed  Google Scholar 

  3. Pan H-A, Wu M-H, Cheng Y-C, Li C-H, Chang F-M. Quantification of Doppler signal in polycystic ovary syndrome using three-dimensional power Doppler ultrasonography: a possible new marker for diagnosis. Human Raproduction. 2002;17(1):201–6.

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  5. Shibuya M. Vascular Endothelial Growth Factor (VEGF) and Its Receptor (VEGFR) Signaling in Angiogenesis: A Crucial Target for Anti- and Pro-Angiogenic Therapies. Genes Cancer. 2011;2(12):1097–105.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lee EJ, Oh B, Lee JY, Kimm K, Park JM, Baek KH. Association study between single nucleotide polymorphisms in the VEGF gene and polycystic ovary syndrome. Fertil Steril. 2008;89(6):1751–9.

    Article  CAS  PubMed  Google Scholar 

  7. Vural P, Kusku-Kiraz Z, Dogru-Abbasoglu S, Cil E, Karadag B, Akgul C, Uysal M. Vascular endothelial growth factor −2578 A/C, −460 T/C and +405 G/C polymorphisms in polycystic ovary syndrome. Eur J Obstet Gynecol Reprod Biol. 2009;147(1):57–60.

    Article  CAS  PubMed  Google Scholar 

  8. Guruvaiah P, Govatati S, Reddy TV, Lomada D, Deenadayal M, Shivaji S, Bhanoori M. The VEGF +405 G > C 5′ untranslated region polymorphism and risk of PCOS: a study in the South Indian Women. J Assist Reprod Genet. 2014;31(10):1383–9.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Vural P, Kusku-Kiraz Z, Dogru-Abbasoglu S, Cil E, Karadag B, Uysal M. Vascular endothelial growth factor +405G/C,-460 T/C and −2578 A/C polymorphisms are not associated with insulin resistance in polycystic ovary syndrome. Int J Immunogenet. 2010;37(4):239–43.

    Article  CAS  PubMed  Google Scholar 

  10. Choudhuri S, Chowdhury IH, Das S, Dutta D, Saha A, Sarkar R, Mandal LK, Mukherjee S, Bhattacharya B. Role of NF-kappaB activation and VEGF gene polymorphisms in VEGF up regulation in non-proliferative and proliferative diabetic retinopathy. Mol Cell Biochem. 2015;405(1–2):265–79.

    Article  CAS  PubMed  Google Scholar 

  11. Ben Salem A, Attaoua R, Mtiraoui N, Belkahla S, Ezzidi I, Ajina M, Souissi M, Poucheret P, Vintila M, Grigorescu F, et al. Common polymorphisms of calpain-10 and the risk of polycystic ovary syndrome in Tunisian population: a case–control study. Mol Biol Rep. 2014;41(10):6569–74.

    Article  CAS  PubMed  Google Scholar 

  12. Ben Salem A, Attaoua R, Mtiraoui N, Meddeb S, Kacem O, Ajina M, Souissi M, Poucheret P, Normand C, Mahjoub T, et al. Haplotyping strategy highlights the specificity of FTO gene association with polycystic ovary syndrome in Tunisian women population. Gene. 2015;565(2):166–70.

    Article  CAS  PubMed  Google Scholar 

  13. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome. Fertility and Sterility. 2004;81(1):19–25.

  14. Palioura E, Diamanti-Kandarakis E. Polycystic ovary syndrome (PCOS) and endocrine disrupting chemicals (EDCs). Rev Endocr Metab Disord. 2015;16(4):365–71.

    Article  CAS  PubMed  Google Scholar 

  15. Davies N. PCOS: Polycystic Ovarian Syndrome. Diabetes Self Manag. 2016;33(1):44–7.

    PubMed  Google Scholar 

  16. Agrawal R, Jacobs H, Payne N, Conway G. Concentration of vascular endothelial growth factor released by cultured human luteinized granulosa cells is higher in women with polycystic ovaries than in women with normal ovaries. Fertil Steril. 2002;78(6):1164–9.

    Article  PubMed  Google Scholar 

  17. Abd El Aal DE, Mohamed SA, Amine AF, Meki AR. Vascular endothelial growth factor and insulin-like growth factor-1 in polycystic ovary syndrome and their relation to ovarian blood flow. Eur J Obstet Gynecol Reprod Biol. 2005;118(2):219–24.

    Article  CAS  PubMed  Google Scholar 

  18. Artini PG, Monti M, Matteucci C, Valentino V, Cristello F, Genazzani AR. Vascular endothelial growth factor and basic fibroblast growth factor in polycystic ovary syndrome during controlled ovarian hyperstimulation. Gynecol Endocrinol. 2006;22(8):465–70.

    Article  CAS  PubMed  Google Scholar 

  19. Artini PG, Ruggiero M, Parisen Toldin MR, Monteleone P, Monti M, Cela V, Genazzani AR. Vascular endothelial growth factor and its soluble receptor in patients with polycystic ovary syndrome undergoing IVF. Hum Fertil (Camb). 2009;12(1):40–4.

    Article  CAS  Google Scholar 

  20. VanKlompenberg MK, Manjarin R, Donovan CE, Trott JF, Hovey RC. Regulation and localization of vascular endothelial growth factor within the mammary glands during the transition from late gestation to lactation. Domest Anim Endocrinol. 2016;54:37–47.

    Article  CAS  PubMed  Google Scholar 

  21. Corbacho AM, Martinez de la Escalera G, Clapp C. Roles of prolactin and related members of the prolactin/growth hormone/placental lactogen family in angiogenesis. J Endocrinol. 2002;173:219–38.

    Article  CAS  PubMed  Google Scholar 

  22. Goldhar AS, Vonderhaar BK, Trott JF, Hovey RC. Prolactin-induced expression of vascular endothelial growth factor via Egr-1. Mol Cell Endocrinol. 2005;232(1–2):9–19.

    Article  CAS  PubMed  Google Scholar 

  23. Malaguarnera L, Imbesi RM, Scuto A, D’Amico F, Licata F, Messina A, Sanfilippo S. Prolactin increases HO-1 expression and induces VEGF production in human macrophages. J Cell Biochem. 2004;93(1):197–206.

    Article  CAS  PubMed  Google Scholar 

  24. Nakajima K-I, Nakamura M, Ishisaki A, Kozakai T. Synergistic Effect of Dexamethasone and Prolactin on VEGF Expression in Bovine Mammary Epithelial Cells via p44/p42 MAP Kinase. Asian-Aust J Anim Sci. 2009;22(6):788–95.

    Article  CAS  Google Scholar 

  25. Seow KM, Lee WL, Wang PH. A challenge in the management of women with polycystic ovary syndrome. Taiwan J Obstet Gynecol. 2016;55(2):157–8.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

This work was supported by the recurrent budget from Laboratory of Human genome and multifactorial diseases (LR12ES07) at University of Monastir Tunisia.

Declaration

This article has been published as part of BMC Genomics Volume 17 Supplement 9, 2016: Proceedings of the 3rd International Genomic Medicine Conference: genomics. The full contents of the supplement are available online at http://0-bmcgenomics-biomedcentral-com.brum.beds.ac.uk/articles/supplements/volume-17-supplement-9.

Availability of data and material

The data sets supporting the results are included within the article.

Authors’ contribution

TM initiated and designed the study, and is responsible for the facility in which the study is conducted. AB wrote the first draft of the manuscript and was responsible for the genetic analysis, the statistical analysis and interpretation of the data. MA, OK, FHB contributed to the original data collection. FM, MS, SH, FJ contributed to the biological and immunoassay. MAE, MA, MHAQ have read the manuscript, revised it, agreed that the work is ready for publication.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

All patients gave written consent for participating and in the study and for publication including consequences that may result.

Ethics approval and consent to participate

The study was approved by the Medical Ethics and Research Committee of the University Hospital Farhat Hached (Sousse).

Publication fees declaration

The publication charges of this manuscript was paid by the Centre of Excellence in Genomic Medicine Research, King Abdulaziz University, Jeddah, Saudi Arabia.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Assila Ben Salem.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ben Salem, A., Megdich, F., Kacem, O. et al. Vascular endothelial growth factor (VEGFA) gene variation in polycystic ovary syndrome in a Tunisian women population. BMC Genomics 17 (Suppl 9), 748 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-3092-5

Download citation

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-3092-5

Keywords