- Open Access
Paraoxonase single nucleotide variants show associations with polycystic ovary syndrome: a meta-analysis
Reproductive Biology and Endocrinology volume 18, Article number: 114 (2020)
Etiology of polycystic ovary syndrome (PCOS) is attributed to genetic and environmental factors. One environmental factor is oxidative stress. Paraoxonase 1 (PON1) is an antioxidant high-density lipoprotein-associated enzyme encoded by the PON1 gene. The PON1 gene has been implicated in the risk for PCOS, the influence of which appears to come from single nucleotide variants (SNVs) at multiple genetic loci. However, association study reports have been inconsistent which compels a meta-analysis to obtain more precise estimates.
From 12 publications, extracted genotype data were used in two genetic procedures. First, linkage disequilibrium (LD) was used to group eight PON SNVs into three: LD1, LD2 and LD3. Second, frequencies of the variant (var), wild-type (wt) and heterozygous (het) genotypes were used for genetic modeling (allele-genotype for LD1 and standard for LD2 and LD3). Risk associations were expressed in terms of pooled odds ratios (ORs), 95% confidence intervals (CIs) and Pa-values. Evidence was considered strong when significance was high (Pa < 0.0001) and heterogeneity absent (I2 = 0%). Pooled effects were subjected to modifier (power), subgroup (Asian/Caucasian), outlier, sensitivity and publication bias treatments. Multiple comparisons were Bonferroni-corrected.
This meta-analysis generated 11 significant outcomes, five in LD1, six in LD2 and none in LD3. All six LD2 outcomes did not survive the Bonferroni-correction but two of the five in LD1 did. These two core LD1 findings conferred greater odds of PCOS to the var allele in the highly significant (Pa < 0.0001) overall (OR 1.44, 95% CI 1.24–1.67) and Asian (OR 1.41, 95% CI 1.20–1.65) outcomes. Of these two core outcomes, the Asian effect was homogeneous (I2 = 0%) but not the overall (I2 = 29%).
Of the eight PON SNVs examined, two (rs854560 and rs662) were associated with PCOS risk. These 1.4-fold increased risk effects rendered Asians susceptible to PCOS. High statistical power, high significance, zero to low-level heterogeneity, robustness and lack of bias in the core outcomes underpinned the strong evidence for association.
Polycystic ovary syndrome (PCOS) is a multifactorial and polygenic disorder . Genetic and environmental factors have a marked influence on the progression of PCOS . Oxidative stress is an environmental factor that plays an important role in the pathogenesis of PCOS  and has become the focus of genetic association studies [4, 5] where single nucleotide variants (SNVs) in genes with anti-oxidant function, like paraoxonase 1 (PON1), have been implicated . The PON1 gene, located on chromosome 7q21.3, is composed of eight introns and nine exons spanning 26 kb and is a member of a multi-gene cluster including PON1, PON2 and PON3 . All three PON genes lie on the long arm of chromosome 7 (7q21–22). PON1 lies near the centromere, while PON2 is near the telomere with PON3 in between them . The PON2 SNV, rs7493 (Ser311Cys), is a substitution of guanine to cytosine nucleotides that results in serine to cysteine amino acid substitutions at residue 311 . Studies have also revealed other SNVs in the PON gene cluster [10, 11]. Several SNVs in PON1 have been studied in regard to its enzyme activity and risk of diseases including PCOS . Among them are five in the promoter region: (i) rs854572 (g.95325384C > G), (ii) rs705381 (g.95324637 T > C), (iii) rs705379 (g.95324583G > A), (iv) rs854571 (g.95325307 T > C) and (v) rs854573 (g.95325551C > T) and two in the coding region (exon 6 and exon 3, respectively): (i) rs662 (c.575A > G, Gln192Arg) and (ii) rs854560 (c.163 T > A, Leu55Met) [13, 14]. These PON1 variants (var) have been reported to regulate PON1 expression and affect circulating serum levels as well as catalytic activity [13, 15]. The rs705379 polymorphism partially regulates PON1 expression  by modulating the binding site for Sp1transcription factor [13, 17] while the rs705381 lies in a potential NF-1 transcription factor binding site . PON1 polymorphisms have been investigated in a variety of pathophysiological conditions that range from metabolic syndrome, cardiovascular diseases and stroke to diabetes. These reports collectively enabled better understanding of the genetic (in terms of transcription factor binding capacity and gene expression level) and physiological pathways conferring the importance of SNVs in oxidative regulation related to the etiology of PCOS. However, associations of the PON SNVs with PCOS risk have been contradictory, which gives reason to perform a meta-analysis of all eligible studies. Since the four previous meta-analyses [18,19,20,21], new primary studies have emerged. We provide a more comprehensive analysis and arrive at a reliable conclusion by reevaluating the associations of the PON SNVs with PCOS risk. To this end, we apply a number of meta-analytical tools that yields fresh insight in the PON-PCOS associations.
Materials and methods
Selection of studies
We searched MEDLINE using PubMed, Science Direct and Google Scholar for association studies as of August 03, 2019. The terms used were “PON”, “paraoxonase”, “polymorphisms”, “polycystic ovary syndrome” and “PCOS” as medical subject heading and text, unrestricted by language and time span. References cited in the retrieved articles were screened manually for additional eligible studies. Inclusion criteria were: (i) case–control studies evaluating the association between PON variants and PCOS and (ii) genotype frequency data to calculate odds ratios (ORs) and 95% confidence intervals (CIs). Exclusion criteria were: (i) animal studies, (ii) reviews, case report or case series, expert opinion and (iii) unusable genotype data.
Two investigators (AK and NP) independently extracted data and arrived at consensus. The following information was obtained from each publication: indications (yes/no) of whether each article was included in the four previous meta-analyses [18,19,20,21], first author’s name, published year, country of origin, ethnicity, diagnostic criteria, PON SNVs examined and article features needed to tally the Clark-Baudouin score. Table S1 shows the rs numbers (SNVs per study), values under cases and controls that include sample sizes and genotype frequencies as well as minor allele frequencies and P-values for the Hardy-Weinberg Equilibrium (HWE).
Statistical power, HWE and data distribution
We used the G*Power program  to evaluate statistical power, where adequacy was set at ≥80% assuming an OR of 1.5 and a genotypic risk level of α = 0.05. Control frequencies from the HWE were calculated from https://ihg.gsf.de/cgi-bin/hw/hwa1.pl with a two-tailed P < 0.05 indicating deviations. Data distribution was assessed with the Shapiro-Wilk test .
Quality assessment of the studies
Methodological quality of the included articles was assessed with the Clark-Baudouin scale , the scores of which range from 0 (worst) to 10 (best) where < 5, 5–7 and > 7 indicate low, moderate and high, respectively.
Linkage disequilibrium (LD) and genetic modeling
We performed LD analysis and adopted a suitable genetic model before examining associations of the PON SNVs with risk of PCOS. Proximity of SNVs has been posited to cause observed phenotype associations  and could merit grouping. Rationale for SNV grouping lies in the concept that SNVs in high LD are assumed to have similar association outcomes. D′ is the metric for LD, in which a value of 1 indicated complete LD . Eight PON SNVs were grouped into three (Table S1) based on D′ values of 0.97–1.00 as LD1 (rs854560 and rs662), LD2 (rs705379, rs7493 and rs854572) and LD3 (rs705381, rs854571, rs854572 and rs854573). Because multiple SNVs had different notations for each genotype, we notated variant and wild-type as var and wt, respectively. Uniformity of the minor allele frequency (all < 0.50) values across the studies in LD2 and LD3 merited use of standard genetic models: (i) recessive: var-var versus het + wt-wt, (ii) dominant: var-var + het versus wt-wt and (iii) codominant: var versus wt . Because of non-uniformity of the minor allele frequencies for LD1, we compared the following: (i) var allele with var-wt/wt-wt genotypes, (ii) wt allele with var-wt/wt-wt genotypes and (iii) wt-var (heterozygous) genotype with homozygous wt-wt and var-var genotypes.
Risk association assessments were confined to HWE-compliant studies . Using raw genotype frequency data, study-specific risks (ORs) of PCOS were estimated and pooled ORs and 95% CIs were calculated by comparing the effects on the same baseline. Associations were considered significant at a two-tailed Pa < 0.05, which were Bonferroni-corrected. Two indicators that were used to assess the strength of evidence were high significance (Pa < 0.00001) and homogeneity or zero heterogeneity (I2 = 0%). Analysis models for absence and presence of heterogeneity were fixed-effects  and random-effects , respectively. Heterogeneity was addressed in the following manner: (i) estimated with the χ2-based Q test  where significance was set at PHET < 0.10 , (ii) quantified with the I2 statistic, which measures the degree of variability between studies  and (iii) its sources examined with outlier treatment using the Galbraith plot . Outlier treatment divided the outcomes into pre-outlier and post-outlier. Robustness of the pooled ORs was assessed with sensitivity analysis, which involved serial omission of the studies followed by re-analysis. Publication bias was assessed with two criteria: (i) statistical significance and (ii) ≥ 10 studies . Normal or non-normal distribution of the operating data (ORs) warranted use of either Egger’s regression asymmetry test  or Begg-Mazumdar test of correlation . Data were analyzed using Review Manager 5.3 (Cochrane Collaboration, Oxford, England), SIGMASTAT 2.03 and SIGMAPLOT 11.0 (Systat Software, San Jose, CA).
Figure 1 outlines the study selection process in a flowchart following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines (Table S2). Twenty-one citations from the initial search were filtered to yield 12 articles for inclusion [38,39,40,41,42,43,44,45,46,47,48,49], all of which examined SNVs in PON1 except Liu et al. , which focused on rs7493 SNV in PON2.
Characteristics of the included studies
Table 1 shows which studies were and were not included in the four previous meta-analyses [18,19,20,21]. Three articles [47,48,49] were new additions to the meta-analysis literature and account for the updated associations between PON variants and PCOS. The sole African study  was included in the overall LD2 analysis, but not in LD1 (where it was HWE-deviating) nor in subgroup analysis. Asian and Caucasian subjects were in six and five articles of LD1, in five and three articles of LD2, respectively. Geography of the Asians was clearly dichotomous (China and India), which could have been sub-stratified were it not for issues of multiplicity (risk of false-positives) and reduced statistical power . Based on the mean ± standard deviation (7.17 ± 1.19) of the normally distributed (Shapiro-Wilks test: P = 0.18) Clark-Baudouin scores, methodological quality of the component studies was high. Table S1 shows the unduplicated total sample sizes of the included articles (n = 7092), five [38, 39, 46,47,48] with adequate statistical power.
This meta-analysis generated 38 comparisons where 11 significant outcomes were confined to LD1 and LD2 (Tables 2 and 3), consigning PON-PCOS associations to five SNVs (rs854560, rs662, rs705379, rs7493 and rs854572).
LD1 associations with PCOS
Table 2 shows 15 comparisons for LD1 (rs854560 and rs662), five of which were significant (Pa < 0.05). Of the five, two (Pa < 0.0001) survived the Bonferroni correction, found robust and both indicating increased risk in the var model. One was in the overall analysis seen in Fig. 2 with variable weight contributions of each study that ranged from 0.3%  to 38% , yet yielded a fixed-effects pooled effect (OR 1.44, 95% CI 1.24–1.67) with low-level heterogeneity (PHET = 0.17, I2 = 29%). This pooled effect had no evidence of publication bias (Begg-Mazumdar test: P = 0.19). The other Bonferroni-surviving outcome was homogeneous (I2 = 0%), and found in Asians (OR 1.41, 95% CI 1.20–1.65). In contrast to the Asians, Caucasian outcomes were non-significant (Pa = 0.44–0.97). The mechanism of outlier treatment is visualized in the het LD1 overall comparison (Figs. 3, 4 and 5). In the pre-outlier forest plot (Fig. 3), the pooled OR was null (OR 1.00, 95% CI 0.78–1.28, Pa = 0.97) and heterogeneous (PHET = 0.0001, I2 = 74%). Using the Galbraith plot, we identified the outlier  located above the + 2 confidence limit (Fig. 4). The post-outlier plot (Fig. 5) showed reduced risk (OR 0.85, 95% CI 0.76–0.96) with eliminated heterogeneity (PHET = 0.60, I2 = 0%), acquired significance (Pa < 0.008), found robust and no evidence of publication bias (Eggers test: P = 0.53). However, this post-outlier pooled effect did not survive the Bonferroni correction.
LD2 and LD3 associations with PCOS
Table 3 shows 18 comparisons for LD2 (rs705379, rs7493 and rs854572), six of which were significant (Pa = 0.005–0.03) with none surviving the Bonferroni correction. Of the six, two indicated reduced risks in the recessive and codominant models, both in power analysis (ORs 0.77–0.82, 95% CIs 0.62–0.97, Pa = 0.002–0.02). In contrast, the remaining four pooled ORs indicated increased risks, three of which were post-outlier derived in the overall analysis (ORs 1.14–1.34, 95% CIs 1.02–1.64, Pa = 0.005–0.03) and one in the pre-outlier Caucasian subgroup of the codominant model (OR 1.32, 95% CI 1.05–1.67, Pa = 0.02). Table 3 shows five comparisons for LD3, none of which were significant (Pa > 0.05). Of the five, four had pooled ORs that skirted the null effect (ORs 0.94–1.02).
Summary of findings
In the main findings, the two LD1 var outcomes met the criterion of high significance. However, only the Asian outcome met the other criterion of homogeneity. Nevertheless, these two core findings identified the PON rs854560 and rs662 SNVs to be associated with risk of PCOS. The LD groups differed by ethnicity, where a significant codominant risk effect in the LD2 was found in Caucasians (OR 1.32, Pa = 0.02) but not in Asians (OR 1.03, Pa = 0.59). In LD1, on the other hand, significant risk effects were found in all genetic models in Asians (Pa = 0.0001–0.03) but not in Caucasians (Pa = 0.44–0.97). Two notes regarding power outcomes (i) LD1 Asian increased risk effects were also powered and homogeneous in the var comparison; and (ii) LD2 power effects were protective (18–23%) in the recessive and codominant models. Reduced risk effects (15–16%) were also observed in the het outcomes of LD1 (overall and Asian), where the phenotypic difference (increased risk versus reduced risk) between var and het suggests a heterosis phenomenon. Heterosis occurs when subjects heterozygous for a specific genetic polymorphism show a different phenotype from homozygotes , conferring a heterozygote advantage (protection). Variations of these effects between the genetic models and ethnic subgroups suggest complex PON-PCOS associations, which is further driven by interactions between genetic and non-genetic risk factors. Gene-gene and gene-environment interactions have been reported to have roles in the associations of PON variants with PCOS. All 12 included articles acknowledged the role of environment but haplotype analysis was addressed in only four (33%) of the component studies [38, 45, 47, 49]. In addition to PON, one article  examined SNVs in another gene (superoxide dismutase-2).
The role of PON1 variants with PCOS have been addressed in four meta-analyses [18,19,20,21], which we compared with the present study in terms of general features and methodology (Table S3). The previous meta-analyses [18,19,20,21] examined PON1 only but our study included PON2 (rs7493) on account of its full LD with rs705379 and rs854572 of PON1. To our knowledge, this is the fifth meta-analysis to address the PON-PCOS associations, but the first to accomplish the following: (i) perform an umbrella review (Table S3); (ii) apply outlier treatment with the most number (n = 12) of included articles; and (iii) operate within LD parameters. These features render our study as most comprehensive, managing to accomplish two things: (i) fill the gaps and update the meta-analysis knowledge of the PON SNV-PCOS associations and (ii) minimize the methodological problems that beset primary studies including limited statistical power, unrecognized confounding factors, misleading definition of phenotypes and stratification of populations .
PCOS is a reproductive endocrinopathy  that is associated with dyslipidemia, obesity and insulin resistance . These metabolic disorders lead to disease conditions such as hypertension, cardiovascular disease and diabetes mellitus . These disease comorbidities were shown to be related with increased oxidative stress, exhibiting altered physiological conditions such as increased plasma glucose and low antioxidant reserves [55, 56]. Low antioxidant levels in PCOS patients suggest that elevated oxidative status contribute to the battery of cardiometabolic derangements . These perturbations have been attributed to reduced serum PON1 activity in PCOS patients , its genetic underpinnings  partly explained with the significant findings of our study observed in LD1, involving rs854560 (L55M) and rs662 (Q192R). These two SNVs account for the main PCOS risk associations in this meta-analysis. The rs854560 variant, not rs662, has been reported to affect enzyme concentration [15, 60]. Compared to 55LL genotype in rs854560, 55MM carriers have lower enzyme activity, which may be attributed to the correlation between the 55 M allele and reduced mRNA and protein levels . Moreover, strong LD of this variant with rs662 may partly explain the variation in PON1 catalytic activity . Depending on the assay used, the R allele of rs662 could be associated with increased or decreased PON1 activity . Reduced enzyme activity leads to elevated levels of oxidative stress altering the metabolism in PCOS patients . Oxidative stress profiles in these patients revealed the R allele carriers had impaired physiologic responses involving increased oxidization of low-density lipoprotein , high triglyceride levels  and elevated insulin resistance . Even in the absence of insulin resistance, oxidative stress levels remain high in PCOS women . This cascade of impaired physiological events, as well as inflammatory responses to cellular injuries caused by oxidative stress, primes the pathophysiology of PCOS as inflammatory mediators that have been known to regulate PON expression  and contribute to PCOS pathogenesis . Moreover, genetic combination analysis demonstrated that haplotypes containing the 192R allele was significantly associated with PCOS risk ranging from 1.6 to 8-fold [38, 49]. In sum, the R allele of rs662 appears to be a genetic determinant of PCOS susceptibility in the female population . These genetically susceptible women are the likely candidates who could benefit from the clinical application of our results. PCOS risks along with their comorbidities could be clinically reduced or delayed with modifications of environmental influences meant to reduce levels of oxidative stress.
Strengths and limitations
Limitations of our study include: (i) six (50%) of the included primary 12 articles were underpowered. However, statistical power at the aggregate level was more than adequate and (ii) credible subgrouping was suggested at no more than two levels , which left other possible subgroups (e.g. diagnostic criteria) unexamined. On the other hand, strengths of this meta-analysis are: (i) the combination of more studies, larger sample sizes and multiple meta-analysis treatments raised the level of evidence presented in this study; (ii) restricting our analysis to HW-compliant studies minimized the risk of representation and methodological bias ; (iii) potency of outlier treatment is evidenced from the LD1 and LD2 analyses, where six (50%) and eight (67%) of the 12 combined post-outlier comparisons were acquired significance and eliminated heterogeneity (Tables 2 and 3); (iv) the overall methodological quality (determined by the Clark-Baudouin Scale) of the included articles was high; (iv) all genotyping used polymerase chain reaction followed by either restriction fragment length polymorphism (11/12 articles) or direct sequencing (1/12 articles) techniques, indicating low-level heterogeneity; (v) all tissue samples were from blood, indicating source homogeneity and (vi) umbrella review of previous meta-analyses enabled comparisons of methodological treatments and findings (Table S3). This provided insight into the evolving consolidation of knowledge into the association genetics of PCOS involving the PON SNVs.
Bonferroni-corrected significance identified rs854560 and rs662 (LD1) as the PON variants associated with PCOS. Subgrouping delineated ethnic-specific effects rendering var carrier Asians susceptible. Future studies exploring other ethnic groups would substantiate conclusions on these PON-PCOS associations with sample sizes appropriate for detecting small genotypic risks.
Availability of data and materials
All data generated or analyzed during this study are included in the supplementary information files.
Goodarzi MO, Dumesic DA, Chazenbalk G, Azziz R. Polycystic ovary syndrome: etiology, pathogenesis and diagnosis. Nat Rev Endocrinol. 2011;7:219–31.
Jones MR, Goodarzi MO. Genetic determinants of polycystic ovary syndrome: progress and future directions. Fertil Steril. 2016;106:25–32.
Murri M, Luque-Ramirez M, Insenser M, Ojeda-Ojeda M, Escobar-Morreale HF. Circulating markers of oxidative stress and polycystic ovary syndrome (PCOS): a systematic review and meta-analysis. Hum Reprod Update. 2013;19:268–88.
Liu J, Zhang D. The role of oxidative stress in the pathogenesis of polycystic ovary syndrome. Sichuan Da Xue Xue Bao Yi Xue Ban. 2012;43:187–90.
Rodrigues JK, Dib LA, Ferriani RA, Jordao Junior AA, Navarro PA. Serum markers of oxidative stress and assisted reproduction procedures results in infertile patients with polycystic ovary syndrome and controls. Rev Bras Ginecol Obstet. 2010;32:118–25.
Agarwal A, Gupta S, Sharma R. Oxidative stress and its implications in female infertility - a clinician’s perspective. Reprod BioMed Online. 2005;11:641–50.
Precourt LP, Amre D, Denis MC, Lavoie JC, Delvin E, Seidman E, Levy E. The three-gene paraoxonase family: physiologic roles, actions and regulation. Atherosclerosis. 2011;214:20–36.
Ng CJ, Shih DM, Hama SY, Villa N, Navab M, Reddy ST. The paraoxonase gene family and atherosclerosis. Free Radic Biol Med. 2005;38:153–63.
Mochizuki H, Scherer SW, Xi T, Nickle DC, Majer M, Huizenga JJ, Tsui LC, Prochazka M. Human PON2 gene at 7q21.3: cloning, multiple mRNA forms, and missense polymorphisms in the coding sequence. Gene. 1998;213:149–57.
Chen Q, Reis SE, Kammerer CM, McNamara DM, Holubkov R, Sharaf BL, Sopko G, Pauly DF, Merz CN, Kamboh MI, Group WS. Association between the severity of angiographic coronary artery disease and paraoxonase gene polymorphisms in the National Heart, Lung, and Blood Institute-sponsored Women’s Ischemia Syndrome Evaluation (WISE) study. Am J Hum Genet. 2003;72:13–22.
Saeed M, Perwaiz Iqbal M, Yousuf FA, Perveen S, Shafiq M, Sajid J, Frossard PM. Interactions and associations of paraoxonase gene cluster polymorphisms with myocardial infarction in a Pakistani population. Clin Genet. 2007;71:238–44.
San Millan JL, Alvarez-Blasco F, Luque-Ramirez M, Botella-Carretero JI, Escobar-Morreale HF. The PON1-108C/T polymorphism, and not the polycystic ovary syndrome, is an important determinant of reduced serum paraoxonase activity in premenopausal women. Hum Reprod. 2006;21:3157–61.
Brophy VH, Jampsa RL, Clendenning JB, McKinstry LA, Jarvik GP, Furlong CE. Effects of 5′ regulatory-region polymorphisms on paraoxonase-gene (PON1) expression. Am J Hum Genet. 2001;68:1428–36.
Rajkovic MG, Rumora L, Barisic K. The paraoxonase 1, 2 and 3 in humans. Biochem Med (Zagreb). 2011;21:122–30.
Leviev I, James RW. Promoter polymorphisms of human paraoxonase PON1 gene and serum paraoxonase activities and concentrations. Arterioscler Thromb Vasc Biol. 2000;20:516–21.
Suehiro T, Nakamura T, Inoue M, Shiinoki T, Ikeda Y, Kumon Y, Shindo M, Tanaka H, Hashimoto K. A polymorphism upstream from the human paraoxonase (PON1) gene and its association with PON1 expression. Atherosclerosis. 2000;150:295–8.
Deakin S, Leviev I, Brulhart-Meynet MC, James RW. Paraoxonase-1 promoter haplotypes and serum paraoxonase: a predominant role for polymorphic position - 107, implicating the Sp1 transcription factor. Biochem J. 2003;372:643–9.
Chen B, Wang Y, Lin W, Jiang H, Liu C, Dong C, Gu H, Kang M, Tang W. Association between paraoxonase-1 (PON-1) polymorphisms and polycystic ovary syndrome susceptibility: a meta-analysis. Int J Clin Exp Med. 2016;9:1490–9.
Gu HF, Mou M, Liang ZG, Sun C, Ren XY, Xiao YB. The association between paraoxonase 1 gene polymorphisms and polycystic ovarian syndrome. Cell Mol Biol (Noisy-le-grand). 2016;62:44–7.
Liao D, Yu H, Han L, Zhong C, Ran X, Wang D, Mo L. Association of PON1 gene polymorphisms with polycystic ovarian syndrome risk: a meta-analysis of case-control studies. J Endocrinol Investig. 2018;41:1289–300.
Liu L, Yang J, Liu Q, Yan W, Wen Y. Risk association of polycystic ovary syndrome with paraoxonase 1 and TNF-α gene polymorphisms: a meta-analysis. Int J Clin Exp Med. 2016;9:7088–101.
Faul F, Erdfelder E, Lang A, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007;39:175–91.
Ghasemi A, Zahediasl S. Normality tests for statistical analysis: a guide for non-statisticians. Int J Endocrinol Metab. 2012;10:486–9.
Clark MF, Baudouin SV. A systematic review of the quality of genetic association studies in human sepsis. Intensive Care Med. 2006;32:1706–12.
Infante-Rivard C. Genetic association between single nucleotide polymorphisms in the paraoxonase 1 (PON1) gene and small-for-gestational-age birth in related and unrelated subjects. Am J Epidemiol. 2010;171:999–1006.
Lewontin RC. On measures of gametic disequilibrium. Genetics. 1988;120:849–52.
Pereira TV, Patsopoulos NA, Pereira AC, Krieger JE. Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication. Int J Epidemiol. 2011;40:457–69.
Trikalinos TA, Salanti G, Khoury MJ, Ioannidis JP. Impact of violations and deviations in Hardy-Weinberg equilibrium on postulated gene-disease associations. Am J Epidemiol. 2006;163:300–9.
Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. 1959;22:719–48.
DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7:177–88.
Lau J, Ioannidis JP, Schmid CH. Quantitative synthesis in systematic reviews. Ann Intern Med. 1997;127:820–6.
Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.
Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002;21:1539–58.
Galbraith RF. A note on graphical presentation of estimated odds ratios from several clinical trials. Stat Med. 1988;7:889–94.
Ioannidis JP, Trikalinos TA. The appropriateness of asymmetry tests for publication bias in meta-analyses: a large survey. CMAJ. 2007;176:1091–6.
Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.
Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50:1088–101.
Wang Y, Liu H, Fan P, Bai H, Zhang J, Zhang F. Evidence for association between paraoxonase 1 gene polymorphisms and polycystic ovarian syndrome in southwest Chinese women. Eur J Endocrinol. 2012;166:877–85.
Zhang Y, Liu H, He J, Xu K, Bai H, Wang Y, Zhang F, Zhang J, Cheng L, Fan P. Lactonase activity and status of paraoxonase 1 in Chinese women with polycystic ovarian syndrome. Eur J Endocrinol. 2015;172:391–402.
San Millan JL, Corton M, Villuendas G, Sancho J, Peral B, Escobar-Morreale HF. Association of the polycystic ovary syndrome with genomic variants related to insulin resistance, type 2 diabetes mellitus, and obesity. J Clin Endocrinol Metab. 2004;89:2640–6.
Mohamed AA, Rashed LA, Abdel Salam RF. Effect of paraoxonase gene polymorphisms on paraoxonase levels and insulin resistance index in women with polycystic ovary syndrome. Aust J Basic Appl Sci. 2009;3:3346–51.
Lenarcik A, Bidzinska-Speichert B, Tworowska-Bardzinska U. The role of chronic inflammation and Leu55Met PON1 polymorphism in the pathogenesis of polycystic ovary syndrome. Gynecol Endocrinol. 2010;26:673–83.
Zhang F, Liu HW, Fan P, Bai H, Song Q. The −108 C/T polymorphism in paraoxonase 1 gene in Chinese patients with polycystic ovary syndrome. Sichuan Da Xue Xue Bao Yi Xue Ban. 2011;42:24–8.
Paltoglou G, Tavernarakis G, Christopoulos P, Vlassi M, Gazouli M, Deligeoroglou E, Creatsas G, Mastorakos G. PON1-108 TT and PON1-192 RR genotypes are more frequently encountered in Greek PCOS than non-PCOS women, and are associated with hyperandrogenaemia. Clin Endocrinol. 2013;79:259–66.
Ferk P, Gersak K. Association of −108 C>T PON1 polymorphism with polycystic ovary syndrome. Biomed Rep. 2014;2:255–9.
Dadachanji R, Shaikh N, Khavale S, Patil A, Shah N, Mukherjee S. PON1 polymorphisms are associated with polycystic ovary syndrome susceptibility, related traits, and PON1 activity in Indian women with the syndrome. Fertil Steril. 2015;104:207–16.
Dadachanji R, Shaikh N, Patil A, Shah N, Mukherjee S. PON1 promoter polymorphisms contribute to PCOS susceptibility and phenotypic outcomes in Indian women. Gene. 2018;661:34–44.
Liu Q, Liu H, Bai H, Huang W, Zhang R, Tan J, Guan L, Fan P. Association of SOD2 A16V and PON2 S311C polymorphisms with polycystic ovary syndrome in Chinese women. J Endocrinol Investig. 2019;42:909–21.
Nalkiran HS, Sahin SB, Ayaz T, Nalkiran I, Guzel AI, Eldes T, Yildiz Y. Association of paraoxonase-1 L55M and Q192R polymorphisms with PCOS risk and potential risk factors for atherosclerosis. Biomark Med. 2019;13:279–89.
Burke JF, Sussman JB, Kent DM, Hayward RA. Three simple rules to ensure reasonably credible subgroup analyses. BMJ. 2015;351:h5651.
Comings DE, MacMurray JP. Molecular heterosis: a review. Mol Genet Metab. 2000;71:19–31.
Carmina E, Lobo RA. Polycystic ovary syndrome (PCOS): arguably the most common endocrinopathy is associated with significant morbidity in women. J Clin Endocrinol Metab. 1999;84:1897–9.
Gonzalez F, Rote NS, Minium J, Kirwan JP. Evidence of proatherogenic inflammation in polycystic ovary syndrome. Metabolism. 2009;58:954–62.
Schroeder BM, American College of O, Gynecologists. ACOG releases guidelines on diagnosis and management of polycystic ovary syndrome. Am Fam Physician. 2003;67:1619–20 1622.
Baskol G, Aygen E, Erdem F, Caniklioglu A, Narin F, Sahin Y, Kaya T. Assessment of paraoxonase 1, xanthine oxidase and glutathione peroxidase activities, nitric oxide and thiol levels in women with polycystic ovary syndrome. Acta Obstet Gynecol Scand. 2012;91:326–30.
Sabuncu T, Vural H, Harma M, Harma M. Oxidative stress in polycystic ovary syndrome and its contribution to the risk of cardiovascular disease. Clin Biochem. 2001;34:407–13.
Fenkci V, Fenkci S, Yilmazer M, Serteser M. Decreased total antioxidant status and increased oxidative stress in women with polycystic ovary syndrome may contribute to the risk of cardiovascular disease. Fertil Steril. 2003;80:123–7.
Marsillach J, Checa MA, Pedro-Botet J, Carreras R, Joven J, Camps J. Paraoxonase-1 in female infertility: a possible role against oxidative stress-induced inflammation. Fertil Steril. 2010;94:1132–4.
Aviram M, Rosenblat M, Bisgaier CL, Newton RS, Primo-Parmo SL, La Du BN. Paraoxonase inhibits high-density lipoprotein oxidation and preserves its functions. A possible peroxidative role for paraoxonase. J Clin Invest. 1998;101:1581–90.
Garin MC, James RW, Dussoix P, Blanche H, Passa P, Froguel P, Ruiz J. Paraoxonase polymorphism Met-Leu54 is associated with modified serum concentrations of the enzyme. A possible link between the paraoxonase gene and increased risk of cardiovascular disease in diabetes. J Clin Invest. 1997;99:62–6.
Brophy VH, Jarvik GP, Richter RJ, Rozek LS, Schellenberg GD, Furlong CE. Analysis of paraoxonase (PON1) L55M status requires both genotype and phenotype. Pharmacogenetics. 2000;10:453–60.
Roest M, van Himbergen TM, Barendrecht AB, Peeters PH, van der Schouw YT, Voorbij HA. Genetic and environmental determinants of the PON-1 phenotype. Eur J Clin Investig. 2007;37:187–96.
Ginsberg G, Neafsey P, Hattis D, Guyton KZ, Johns DO, Sonawane B. Genetic polymorphism in paraoxonase 1 (PON1): population distribution of PON1 activity. J Toxicol Environ Health B Crit Rev. 2009;12:473–507.
Mohamadin AM, Habib FA, Elahi TF. Serum paraoxonase 1 activity and oxidant/antioxidant status in Saudi women with polycystic ovary syndrome. Pathophysiology. 2010;17:189–96.
Torun AN, Vural M, Cece H, Camuzcuoglu H, Toy H, Aksoy N. Paraoxonase-1 is not affected in polycystic ovary syndrome without metabolic syndrome and insulin resistance, but oxidative stress is altered. Gynecol Endocrinol. 2011;27:988–92.
Fuhrman B. Regulation of hepatic paraoxonase-1 expression. J Lipids. 2012;2012:684010.
Rosenfield RL, Ehrmann DA. The pathogenesis of polycystic ovary syndrome (PCOS): the hypothesis of PCOS as functional ovarian hyperandrogenism revisited. Endocr Rev. 2016;37:467–520.
We thank Jakkree Tangthianchaichana for acquiring the key literature.
The authors declare that this study was unfunded.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Quantitative features of the included PON-PCOS studies. Table S2. Preferred Reporting Items for Systematic Reviews and Meta-Analyses checklist. Table S3. Comparisons between meta-analyses that examined the PON variants associations with PCOS.
About this article
Cite this article
Kunjantarachot, A., Pabalan, N., Jarjanazi, H. et al. Paraoxonase single nucleotide variants show associations with polycystic ovary syndrome: a meta-analysis. Reprod Biol Endocrinol 18, 114 (2020). https://doi.org/10.1186/s12958-020-00665-1
- Single nucleotide variants