The promoter T-413A variant and elevated enzyme levels of heme oxygenase-1 associated with an increased risk of polycystic ovarian syndrome

Background: Oxidative stress and metabolic disorders significantly contribute to the development of polycystic ovarian syndrome (PCOS). Heme oxygenase-1 (HMOX1) plays a key role in the degradation of heme and the regulation of oxidative stress, ferroptosis, and glycolipid metabolism. This study explored the relationship between HMOX1 promoter T-413A single nucleotide polymorphism (SNP, rs2071746), (GT)n dinucleotide repeat variant (rs3074372), plasma HMOX1 levels, and the risk of PCOS in Chinese women.

Methods: This case-control study included 1092 women diagnosed with PCOS and 805 controls. The (GT)n and rs2071746 polymorphisms were identified using polymerase chain reaction amplification, followed by capillary electrophoresis or restriction fragment length polymorphism. HMOX1 levels and clinical, metabolic, hormonal, and oxidative stress indices were analyzed.

Results: The HOMX1 rs2071746T/A SNP was associated with an increased risk of PCOS based on genotype, recessive, dominant, and allele genetic models (P < 0.05). After adjusting for age, body mass index, and recruitment year of participants, the dominant model (odds ratio [OR] = 1.272, 95% confidence interval [CI]: 1.013–1.597, P = 0.039) and the TT genotype (OR = 1.395, 95% CI: 1.033–1.883, P = 0.030, with the AA genotype as the reference) remained a significant predictor of PCOS in the logistic regression models. No significant differences were observed in the (GT)n polymorphism of HMOX1 based on different genetic models. However, the TT/SS combined genotype of HMOX1 rs2071746T/A and (GT)n polymorphisms was associated with an increased risk of PCOS (OR = 1.442, 95% CI: 1.021–2.035, P = 0.037). Furthermore, elevated HMOX1 levels were related to a slight but significant increase in the risk of PCOS, and the rs2071746T/A and (GT)n genetic variants significantly affected obesity, oxidative stress, endocrine abnormalities, and metabolic disorders.

Conclusion: HMOX1 rs2071746T/A variant and elevated plasma HMOX1 levels are associated with an increased risk of PCOS.

1 Introduction

Polycystic ovarian syndrome (PCOS) is the leading endocrine-metabolic disorder among women of reproductive age, with a prevalence of 10–13% (1). Its clinical presentation is diverse and may include reproductive, dermatological, metabolic, and psychological symptoms such as irregular menstrual cycles, infertility, acne, hirsutism, obesity, insulin resistance, dyslipidemia, anxiety, depression, etc. (1, 2). PCOS can affect women throughout their lives and potentially lead to long-term complications, including type 2 diabetes (T2D), cardiocerebrovascular disease, and endometrial cancer. The exact etiology of PCOS remains unknown, making its treatment difficult (1, 3). Increasing evidence indicates that its etiology is complex, involving interactions between multiple predisposing genes, genetic epigenetics, and detrimental environmental factors (2, 4–6).

Heme oxygenase (HMOX) is the rate-limiting enzyme that catalyzes the breakdown of heme, resulting in the production of biliverdin (BV), carbon monoxide (CO), and free iron (Fe²⁺) (7, 8). BV is rapidly turned into bilirubin (BR) by BV reductase (8, 9). Heme and free Fe²⁺ are toxic owing to their oxidative properties; however, in the physiological state, free iron is swiftly sequestered by ferritin, ensuring a harmonious balance within the body (8). CO, BV, and BR possess vasodilator, antioxidant, and anti-inflammatory properties, but excessive accumulation of these products can be toxic (8, 10). The two main isomers of HMOX in humans, inducible HMOX1 (also known as HO-1) and constitutive HMOX2 (also known as HO-2), have similar structures and catalytic functions but are distributed in different tissues and exhibit distinct features (10). Low HMOX1 expression was observed in most tissues under normal conditions. However, its expression can be markedly increased in response to different pathophysiological stress conditions or stimulation factors (9, 11).

The HMOX1 is located on chromosome 22q13.1. Two genetic polymorphisms in the promoter of HMOX1, the (GT)n dinucleotide repeat variant (rs3074372) and rs2071746T/A single nucleotide polymorphism (SNP) (rs2071746), can affect the transcriptional activity of HMOX1 (11). These two polymorphisms are closely linked to certain diseases, including sensitivity to several cancers and coronary heart disease (11), cardiovascular events and mortality in patients undergoing hemodialysis (12), T2D (13), pre-eclampsia (14), chronic obstructive pulmonary disease (15), SARS-CoV-2 viremia in COVID-19 infection (16), and risk of encephalitis in HIV infection (17).

Oxidative stress and metabolic disorders significantly contribute to the pathophysiology and progression of PCOS (1, 18–21). Genetic variants of HMOX2 G554A and A-42G SNPs are associated with endocrine abnormalities and glycolipid metabolic irregularities in patients with PCOS (5). The levels of HMOX1 mRNA are higher in subcutaneous adipose tissue and granulosa cells (22, 23), but the concentrations of serum HMOX1 are lower in women with PCOS than those among control women (24). However, the association between HMOX1 polymorphisms and PCOS remains unclear. Therefore, we explored the relationship between HMOX1 (GT)n repeats and rs2071746T/A polymorphisms and the risk of PCOS. Additionally, we analyzed how these genetic variants affected plasma HMOX1 levels and various clinical and biochemical parameters in Southwest Chinese women.

2 Materials and methods

2.1 Study participants

This was a case-control study. All participants aged 17–40 years provided written informed consent and were recruited from the Reproductive Endocrinology Outpatient Department of the West China Second University Hospital between 2006 and 2024 (Figure 1). This study was approved by the Institutional Review Board of West China Second University Hospital, Sichuan University (2014–014 for P. Fan).

Figure 1

Figure 1. The process of recruitment and selction of the case and control group.

PCOS was diagnosed based on the revised 2003 Rotterdam European Society of Human Reproduction and Embryology/American Society for Reproductive Medicine consensus criteria, which include irregular menstrual cycles, oligo-ovulation, or anovulation (OA), biochemical and/or clinical hyperandrogenism (HA), and polycystic ovaries (PCOs) (25). Detailed definitions of HA, OA, and PCOs have been provided in previous articles (26–28). PCOS was diagnosed if a woman met at least two of the three criteria after ruling out other etiologies such as Cushing syndrome, congenital adrenal hyperplasia, and androgen-secreting tumors (1, 25). Both HA and OA are required in individuals aged < 20 years of age (1). Control women had regular menstrual cycles between 21–35 days, normal ovarian morphology on ultrasonography, and no biochemical or clinical signs of HA.

The participants were excluded if they had infections, cardiovascular diseases, liver or kidney diseases, autoimmune diseases, thyroid disorders, hypogonadism, hyperprolactinemia, premature ovarian insufficiency, endometriosis, or other malignancies. Additionally, the participants were excluded from the analysis when comparing oxidative stress and metabolic and hormonal parameters between groups/subgroups if they met any of the following criteria: (i) use of hormonal therapy and/or medications that influence glucose and lipid metabolism within 12 weeks before the study; (ii) smoking; (iii) being pregnant or in the luteal phase of their menstrual cycle; and/or (iv) having a fasting glucose (Glu) level of ≥7.0 mmol/L and/or a 2-hour plasma glucose after the glucose challenge (2-h Glu) level of ≥11.1 mmol/L in the control group.

Clinical indices, including waist circumference (WC), body mass index (BMI), waist-to-hip ratio, diastolic blood pressure (DBP), systolic blood pressure, severity of acne and hirsutism, and ovarian volume were assessed as previously described (27–29).

Blood samples were collected after fasting for 8–12 h. Blood cells were stored at 4°C, and plasma and serum aliquots were preserved at −80°C for later analysis. A 75g oral glucose tolerance test was conducted immediately after fasting blood sampling.

2.2 DNA purification and genotype measurements

Genomic DNA was purified from stored blood cells using a previously described method (30). (GT)n repeats in HMOX1 were determined using polymerase chain reaction (PCR) amplification and capillary electrophoresis. The forward primer with 6-carboxyfluorescein (FAM) was 5′-FAM-CCAGCTTTCTGGAACCTTCTG-3′, the reverse primer was 5′-GAAACAAAGTCTGGCCATAGGA-3′ (31). Samples were amplified using a touchdown PCR protocol (17). The products were then analyzed using a 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA) and GeneMapper 4.1 (Applied). Short repeats, with fewer than 27 GT repeats, were classified as S alleles, whereas long repeats, with at least 27 GT repeats, were classified as L alleles (12). For genotyping the rs2071746 SNP, we used PCR and restriction fragment length polymorphism method with a mismatched primer set (the forward primer: 5′-GTTCCTGATGTTGCCCACCAAGC-3′; the reverse primer: 5′-CTGCAGGCTCTGGGTGTGATTTTG-3′). The PCR products of the rs2071746T/A SNP (151 bp) were then digested with HindIII (New England Biolabs, Ipswich, MA, USA), resulting in 20 and 131 bp fragments with the T allele and a whole 151-bp product with the A allele. The results were verified by repeating the genotyping of > 30% of the randomly selected samples, yielding 100% concordance.

2.3 Analysis of HMOX1, oxidative stress, hormonal, and metabolic indices

Plasma HMOX1 concentrations were measured using ELISA kits (Elabscience Biotechnology Co., Ltd., Wuhan, China).

Estradiol, TT, luteinizing hormone (LH), follicle-stimulating hormone (FSH), sex hormone binding globulin (SHBG), plasma insulin (Ins) and Glu levels, triglycerides (TG), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), apolipoprotein (apo) A1, and apoB concentrations, total antioxidant capacity (T-AOC), glutathione (GSH), total oxidant status (TOS), and oxidative stress index (OSI) were also analyzed as previously described in other studies (20, 21, 28). The homeostatic model assessment of insulin resistance (HOMA-IR) and free androgen index (FAI) were calculated as previously described (18, 28) using the following formulas:

$$HOMA-IR=fastingGlu(mmol/L)\times fastingIns(\mu U/mL)/22.5$$

$$FAI=TT(nmol/L)/SHBG(nmol/L)\times 100$$

2.4 Statistical analysis

Data are presented as mean ± standard deviation. We used analysis of variance or independent sample t-tests for normally distributed variables and the Mann–Whitney U test for non-normally distributed variables. Analysis of covariance was used to evaluate differences in clinical and biochemical indices after adjusting for variations in age, BMI, and participant recruitment year. Chi-squared (χ²) analysis was performed to evaluate deviations in genotypic distribution from Hardy–Weinberg equilibrium and to compare the frequencies of genotypes and alleles between two groups. The Spearman’s correlation coefficient was used to assess the relationship between HMOX1 levels and other parameters. Differences were considered statistically significant if the P value was <0.05. Data were analyzed using the Statistical Program for Social Sciences (SPSS) 21.0 (IBM SPSS Statistics, IBM Corporation).

Power values were calculated according to the disease allele frequency of the rs2071746T/A SNP in HMOX1 and the sample size (prevalence = 0.12, significance level = 0.05) using the Genetic Association Study Power Calculator (http://csg.sph.umich.edu/abecasis/gas_power_calculator/index.html).

Using the online SNPStats (https://www.snpstats.net/start.htm), we assessed the linkage disequilibrium between two genetic polymorphic loci in view of the D’ parameter.

3 Results

3.1 Clinical and biochemical characteristics of the participants

Owing to the significant discrepancies in BMI and age between the PCOS and control groups (Table 1), and the relatively long recruitment period of participants between 2006 and 2024 in this study, we adjusted for these three confounding factors in our subsequent analyses.

Table 1

Table 1. Clinical and biochemical parameters in women with PCOS and controls.

Table 1 shows that the PCOS group had significantly higher acne grade scores, F-G scores, average ovarian volumes, DBP, WC, waist-to-hip ratio, TT, LH, LH/FSH ratio, FAI, fasting Ins, HOMA-IR, 2-h Glu and 2-hour insulin after the glucose challenge (2-h Ins), LDL-C, TC, TG, TG/HDL-C ratio, apoB, TOS, T-AOC, GSH, OSI, TOS/GSH ratio, and plasma HMOX1 levels, but lower serum FSH, SHBG, and HDL-C concentrations than the control group (P < 0.05).

We further compared the plasma HMOX1 levels in lean (BMI < 23 kg/m²) and overweight/obese (BMI ≥ 23 kg/m² and/or waist circumference > 80 cm) subgroups after adjusting for age, BMI, and recruitment period of participants. The results showed that the HMOX1 levels were higher in the lean PCOS subgroup (n = 339) than in the lean control subgroup (n = 382) (5.00 ± 5.00 vs. 4.48 ± 2.36 µg/L, P = 0.034), but no statistical significance in the overweight/obese PCOS subgroup (n = 387) than in the overweight/obese control subgroup (n = 134) (5.04 ± 4.29 vs. 4.48 ± 2.15 µg/L, P = 0.115).

3.2 Correlation of HMOX1 levels with clinical and biochemical indicators and risk of PCOS

The Spearman’s correlation analysis showed that plasma HMOX1 levels were positively correlated with 2-h Glu, WC, fasting Ins, HOMA-IR, TG/HDL-C ratio, BMI, 2-h Ins, TC, fasting Glu, FAI, T-AOC, apoB, and WHR in patients with PCOS (r = 0.138, 0.132, 0.132, 0.129, 0.127, 0.118, 0.105, 0.093, 0.092, 0.087, 0.085, 0.080, and 0.079, respectively; P < 0.05). Although statistically significant, the correlations between HMOX1 and PCOS traits were quantitatively modest.

Binary logistic regression analysis demonstrated that elevated HMOX1 levels were associated with an increased risk of PCOS after correcting for differences in participant recruitment year, age, and BMI (odds ratio [OR] = 1.053, 95% confidence interval [CI]: 1.008–1.100, P = 0.019).

3.3 Distributions of HOMX1 rs2071746T/A and (GT)n genotypes and alleles

Table 2 summarizes the genetic models for the rs2071746T/A and (GT)n repeat polymorphisms in HOMX1. The distribution of genotypes for both polymorphisms was consistent with Hardy–Weinberg equilibrium in women with and without PCOS (P > 0.05).

Table 2

Table 2. Association of HMOX1 T-413A (rs2071746) and (GT)n repeat polymorphisms with the risk of PCOS using different genetic models.

The frequencies of the TT genotype and T allele in the HOMX1 rs2071746T/A SNP were significantly higher in the PCOS group than those in the control group. The OR indicated that this difference was statistically significant for the dominant model, the recessive model, and the TT vs. AA genotype model, and the allele model (all P < 0.05). After adjusting for age, BMI, and recruitment year of participants, the dominant genetic model remained statistically significant in the binary logistic regression model (OR = 1.272, 95% CI: 1.013–1.597, P = 0.039) and the TT genotype remained a significant predictor for PCOS in a multinomial logistic regression model, with the AA genotype as the reference (OR = 1.395, 95% CI: 1.033–1.883, P = 0.030). The genetic association power is 0.984 for rs2071746T/A SNP. No statistically significant differences were observed between the two groups for the (GT)n repeat polymorphism of HMOX1 when analyzed using different genetic models (P > 0.05; Table 2).

The combined genotypes of HMOX1 rs2071746T/A and (GT)n polymorphisms exhibited a significant difference in frequency between patients with PCOS and controls (P = 0.031; Supplementary Table 1). The TT/SS was a risk factor for PCOS (OR = 1.442, 95% CI: 1.021–2.035, P = 0.037) in a multinomial logistic regression model using the AA/LL combined genotype as the reference, with participant recruitment year, age, and BMI as covariates. Moderate linkage disequilibrium was observed between the rs2071746T/A and (GT)n polymorphisms (D’= 0.8517, r² = 0.5539).

3.4 Effects of genotypes on clinical and biochemical indicators

We analyzed the effect of HMOX1 rs2071746T/A and (GT)n genetic variants on plasma HMOX1 levels and clinical and biochemical parameters in women with and without PCOS.

Supplementary Table 2 shows that patients with the AT genotype of HMOX1 rs2071746T/A SNP had a greater TOS/GSH ratio than those with the TT genotype (P = 0.028). The controls with the TT genotype exhibited a lower acne grade score than those with the AA genotype (P = 0.043) and lower HDL-C levels (P = 0.038) than those with the AT genotype; whereas the controls with the AT genotype exhibited lower GSH levels (P = 0.022) than those with the AA genotype.

The same parameters were analyzed for different genotypes of HMOX1 (GT)n repeat polymorphism (Supplementary Table 3). Patients with the LL genotype displayed a lower waist-to-hip ratio and TT levels (P < 0.05) than those with the SS genotype. The controls with the LL genotype exhibited a lower BMI than those with the SS and SL genotypes (P < 0.05). The controls with the SL genotype showed higher FAI (P = 0.019) than those with the SS genotype.

No statistically significant differences in plasma HMOX1 levels were observed between the different genotypes of HMOX1 rs2071746T/A and (GT)n genetic variants in the control and PCOS groups (P > 0.05; Supplementary Tables 2, 3).

4 Discussion

For the first time, we demonstrated that the TT genotype and T allele of the rs2071746T/A SNP are associated with an increased risk of PCOS in Chinese women. We also proved that the TT/SS combined genotype of the rs2071746T/A and (GT)n repeat variants is a risk factor for PCOS. Furthermore, we found that plasma HMOX1 levels were significantly higher in patients with PCOS than those in the control women, and elevated HMOX1 levels were related to a slight but significant increase in the risk of PCOS, suggesting that patients with PCOS have a compensatory increase in HMOX1 levels. HMOX1 rs2071746T/A and (GT)n repeat polymorphisms significantly affected BMI, waist-to-hip ratio, TT, FAI, acne grade score, HDL-C, GSH, and TOS/GSH ratio, but not plasma HMOX1 levels among the PCOS and/or control participants, supporting that the two variants may be involved in obesity, endocrine abnormalities, oxidative stress, and metabolic disorders.

Oxidative stress, metabolic disorders, and iron homeostasis imbalance play significant roles in the occurrence and progression of PCOS (18–21, 32). HMOX catalyzes the degradation of heme and is crucial for controlling the dynamic equilibrium of heme and its products (BV, BR, CO, and Fe²⁺) (7, 8). In addition to the recovery of Fe²⁺ from heme, HMOX participates in the regulation of multiple signaling pathways via its products, BV, CO, and Fe²⁺, as well as its substrate, heme (8). Under physiological conditions, low HMOX1 expression is found in most tissues, except in cells of the reticuloendothelial system (8, 10). Unlike the constitutive isoform HMOX2, which is barely regulated at the transcriptional and translational levels (5, 8), HMOX1 can be rapidly induced under various stress conditions (7, 8, 10, 11). BV and BR are important endogenous antioxidants and cellular signaling molecules that play significant roles in regulating immunity and glycolipid metabolism, and CO is a gaseous mediator with vasodilatory, anti-inflammatory, anti-proliferative, and anti-apoptotic properties, while heme and free Fe²⁺ can facilitate the production of reactive oxygen species (7, 8, 10). Therefore, besides its cytoprotective effects, HMOX1 induction may also be involved in the development of certain diseases. Several studies have indicated that genetic overexpression or chemical induction of HMOX1 can protect against hypertension, cardiovascular diseases, metabolic conditions, and kidney diseases (33–35). Increased HMOX1 activity may promote oxidative stress by increasing free intracellular iron and accelerating the consumption of cytosolic NADPH, thereby contributing to chronic inflammation, ferroptosis, and cell injury (8). It has been reported that elevated plasma HMOX1 levels in individuals with T2D are associated with a higher disease risk (36). However, another study showed that low serum HMOX1 levels in non-obese women are an independent risk factor for PCOS (24). In this study, we found that plasma HMOX1 levels were significantly higher in the PCOS group compared to the control group, and the lean PCOS subgroup compared to the lean control subgroup. Furthermore, elevated HMOX1 levels were related to an increased risk of PCOS. The possible reasons for the inconsistent results of HMOX1 levels in PCOS may be discrepancies in the sample size and study population. Contrary to the report of lower HMOX1 concentrations in PCOS (24), our finding of elevated HMOX1 levels could represent a protective compensatory response to various chronic unfavorable stimuli in PCOS, aligning with the canonical role of HMOX1 as an oxidative stress sensor.

Increased oxidative stress in PCOS, as shown in this study, can enhance HMOX1 transcription by activating nuclear factor erythroid 2-related factor 2 transcription factor (11, 37). Patients with PCOS have iron overload, abnormal heme metabolism, and chronic inflammation due to chronic oligomenorrhea, excessive androgen, and compensatory hyperinsulinemia (32, 38). High levels of heme and activation of the inflammatory factor nuclear factor kappa B can promote the expression of HMOX1 (11, 33). The HMOX1 and its downstream metabolites, including CO, BV, and BR, may play a protective role through their antioxidant, anti-inflammatory, and vasodilator functions (8, 10). However, the sustained induction of HMOX1 and the iron overload may paradoxically dysregulate ferroptosis through iron-mediated production of peroxidized lipids, potentially contributing to ovarian dysfunction. In addition to transcriptional regulation, HMOX1 activity is regulated by critical protein-protein interactions (PPIs) and post-translational modifications (39). PPIs affect the stability, oligomerization, subcellular localization, and function of HMOX1 (39–42). The post-translational modifications, such as phosphorylation, acetylation, and ubiquitination, also play an important role in regulating the level and activity of HMOX1 (43–45). In a word, the possible reasons and mechanisms for elevated HMOX1 levels are complex in PCOS. Further studies are needed to explain this phenomenon and its exact pathophysiological mechanisms.

Genetic variants in the promoter of HMOX1 may affect the expression of HMOX1 (11), thereby affecting the incidence and progression of diseases. The A allele of HMOX1 rs2071746T/A SNP is associated with a higher transcription activity of HMOX1 (11, 46). A meta-analysis revealed a lower susceptibility to coronary heart disease in individuals carrying the A allele (11). The AA genotype increases the occurrence of hypertension in the Japanese women (47), but decreases the risk of ischemic heart disease in the Japanese population (48). The T allele of the rs2071746T/A variant is a risk factor for the development of esophageal varices in patients with cirrhosis (46), and the TT genotype is more likely to cause proteinuria in Korean patients with T2D (49) and SARS-CoV-2 viremia in COVID-19 infection (16). Our findings indicated that women with the TT genotype and T allele have a higher risk of developing PCOS. Moreover, we found that this genetic polymorphism may contribute to oxidative stress, hyperandrogenism status, and metabolic disorders through influencing TOS/GSH ratio, GSH and HDL-C levels, and acne grade score in the study population. However, we did not observe significant differences in plasma HMOX1 levels according to different genotypes of the rs2071746T/A SNP, suggesting that this genetic variant may not be a key factor affecting HMOX1 expression in the study population.

Of the polymorphisms observed in the HMOX1 promoter region, the (GT)n repeat variant has been extensively studied (11, 13). The range of (GT)n repeat numbers is 10–50 (17), and the number of repeat lengths shows a bimodal distribution, with peaks at (GT)₂₃ and (GT)₃₀ repeats in East Asian and Caucasian populations, and a trimodal form, with crest values at (GT)₂₃, (GT)₃₀, and (GT)₃₉ among African-Americans (11). Generally, the short (S) allele is defined as the number of (GT)n repeats < 25 or 27, and the long (L) allele is defined as the number of (GT)n repeats ≥ 25 or 27 in different reports (11–13, 49). The S alleles are linked to increased transcriptional activity compared with the L alleles (11, 13), and individuals with the SS genotype of the (GT)n variant have higher levels of HMOX1 mRNA than those with the LL genotype (50). Individuals with the S allele or SS genotype have a reduced risk of coronary heart disease (11), T2D (13), rheumatoid arthritis (50), and encephalitis in HIV infection (17), but an increased risk of melanoma (51). Whereas individuals with the L allele or LL genotype have an increased risk of hypertension (52), chronic obstructive pulmonary disease (15), and preeclampsia (late-onset and non-severe form) (14). Our study revealed no significant differences were observed between the PCOS and control groups based on different genetic models. However, this polymorphism may be involved in obesity and endocrine disorders, probably by affecting waist-to-hip ratio, BMI, TT levels, and FAI in the study population. Moreover, we did not observe significant differences in plasma HMOX1 levels according to different genotypes of the (GT)n repeat variant.

Additionally, in this study, a moderate linkage disequilibrium was observed between the rs2071746T/A and (GT)n repeat polymorphisms. The TT/SS combined genotype of the two genetic polymorphisms was a risk factor for PCOS. However, we did not observe significant differences in the HMOX1 levels according to the different combined genotypes in the control and PCOS groups (P > 0.05, data not indicated). Several studies have shown that HMOX1 expression is also controlled by genetic polymorphisms (11). The (GT)n and rs2071746T/A polymorphisms may be involved in 5’-UTR alternative splicing of the HMOX1 primary transcript, which may affect translational efficiency and mRNA stability, and thus regulate the translational process of HMOX1 (53). Therefore, a more detailed and systematic investigation of the correlation between genotype and gene expression, along with in vitro studies, is required to clarify the potential mechanism.

The present study has some limitations. First, we did not measure the levels of BV, BR, and iron, which are critical downstream products of HMOX1 enzymatic activity; this could provide further evidence to reveal the relationships between HMOX1 genetic variants, HMOX1 levels, and PCOS and the underlying mechanism. Second, we did not evaluate oxidative stress and hormonal and metabolic indices because of confounding factors in some participants, which could have had an effect on the statistical effectiveness of these parameters. Third, lifestyle factors (e.g., dietary patterns, physical activity levels) and longitudinal treatment histories were not collected during enrollment, preventing us from adjusting for differences in these potential confounders in our analyses.

In conclusion, this study indicates that the HMOX1 rs2071746T/A SNP is associated with the risk of PCOS, and that the T allele, TT genotype, and its coexistence with the SS genotype of the (GT)n repeat variant are genetic risk factors for PCOS among Chinese women. We further demonstrated that patients with PCOS have higher plasma HMOX1 concentrations and that elevated HMOX1 levels are associated with an increased risk of PCOS. We found that the HMOX1 rs2071746T/A and (GT)n repeat polymorphisms may contribute to obesity, oxidative stress, endocrine abnormalities, and metabolic disorders. Our findings suggest that induction of the heme-degrading enzyme HMOX1 and its genetic polymorphisms in the promoter may be involved in the pathophysiology of PCOS.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Institutional Review Board of West China Second University Hospital, Sichuan University (2014-014 to PF). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

QW: Formal Analysis, Methodology, Writing – original draft, Investigation. JL: Resources, Data curation, Writing – original draft, Investigation. QL: Investigation, Writing – original draft, Project administration, Validation. HL: Writing – review & editing, Resources, Supervision, Data curation. HB: Writing – review & editing, Supervision. WH: Writing – review & editing, Supervision, Resources, Conceptualization. PF: Formal Analysis, Methodology, Conceptualization, Funding acquisition, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Natural Science Foundation of China (81370681) and the Program for Chang Jiang Scholars and Innovative Research Team in the University, Ministry of Education (IRT0935).

Acknowledgments

The authors thank all participants who donated blood samples for this study. We also thank You Li, Qi Song, Dehua Wan, Ying Wang, Feng Zhang, Jinxia Zhang, Yan Gong, Renjiao Zhang, Jun Tan, and Yujin Zhang for their support with sample and data collection.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1644373/full#supplementary-material

Abbreviations

ApoA1, apolipoprotein A1; ApoB, apolipoprotein B; BMI, body mass index; BR, bilirubin; BV, biliverdin; CI, confidence interval; CO, carbon monoxide; DBP, diastolic blood pressure; E₂, estradiol; FAI, free androgen index; F-G score, Ferriman–Gallwey score; FSH, follicle-stimulating hormone; Glu, glucose; GSH, glutathione; HA, hyperandrogenism; HDL-C, high-density lipoprotein cholesterol; HMOX1, heme oxygenase-1; HMOX2, heme oxygenase-2; HOMA-IR, the homeostatic model assessment of insulin resistance; Ins, insulin; LDL-C, low-density lipoprotein cholesterol; LH, luteinizing hormone; OA, oligo-ovulation or anovulation; OR, odds ratio; OSI, oxidative stress index; PCOs, polycystic ovaries; PCOS; polycystic ovarian syndrome; PCR, polymerase chain reaction; SBP, systolic blood pressure; SHBG, sex hormone-binding globulin; SNP, single nucleotide polymorphism; T2D, type 2 diabetes; T-AOC, total antioxidant capacity; TC, total cholesterol; TG, triglycerides; TOS, total oxidant status; TT, total testosterone; WC, waist circumference.

References

1. Teede HJ, Tay CT, Laven J, Dokras A, Moran LJ, Piltonen TT, et al. Recommendations from the 2023 international evidence-based guideline for the assessment and management of polycystic ovary syndrome. Fertil Steril. (2023) 120:767–93. doi: 10.1016/j.fertnstert.2023.07.025

PubMed Abstract | Crossref Full Text | Google Scholar

2. Escobar-Morreale HF. Polycystic ovary syndrome: definition, aetiology, diagnosis and treatment. Nat Rev Endocrinol. (2018) 14:270–84. doi: 10.1038/nrendo.2018.24

PubMed Abstract | Crossref Full Text | Google Scholar

3. Joshi A. PCOS stratification for precision diagnostics and treatment. Front Cell Dev Biol. (2024) 12:1358755. doi: 10.3389/fcell.2024.1358755

PubMed Abstract | Crossref Full Text | Google Scholar

4. Palumbo M, Della Corte L, Colacurci D, Ascione M, D'Angelo G, Baldini GM, et al. PCOS and the genome: is the genetic puzzle still worth solving? Biomedicines. (2025) 13:1912. doi: 10.3390/biomedicines13081912

PubMed Abstract | Crossref Full Text | Google Scholar

5. Zhang X, Li S, Liu H, Bai H, Liu Q, Yang C, et al. Heme oxygenase 2 genetic variants alter hormonal and metabolic traits in polycystic ovary syndrome. Endocr Connect. (2024) 13:e230463. doi: 10.1530/EC-23-0463

PubMed Abstract | Crossref Full Text | Google Scholar

6. Merkin SS, Phy JL, and Sites CK. Yang D Environmental determinants of polycystic ovary syndrome. Fertil Steril. (2016) 106:16–24. doi: 10.1016/j.fertnstert.2016.05.011

PubMed Abstract | Crossref Full Text | Google Scholar

7. Haines DD. Tosaki A heme degradation in pathophysiology of and countermeasures to inflammation-associated disease. Int J Mol Sci. (2020) 21:9698. doi: 10.3390/ijms21249698

PubMed Abstract | Crossref Full Text | Google Scholar

8. Duvigneau JC and Esterbauer H. Kozlov AV role of heme oxygenase as a modulator of heme-mediated pathways. Antioxidants (Basel). (2019) 8:475. doi: 10.3390/antiox8100475

PubMed Abstract | Crossref Full Text | Google Scholar

9. Kim YM, Pae HO, Park JE, Lee YC, Woo JM, Kim NH, et al. Heme oxygenase in the regulation of vascular biology: from molecular mechanisms to therapeutic opportunities. Antioxid Redox Signal. (2011) 14:137–67. doi: 10.1089/ars.2010.3153

PubMed Abstract | Crossref Full Text | Google Scholar

10. Nitti M and Furfaro AL. Mann GE heme oxygenase dependent bilirubin generation in vascular cells: A role in preventing endothelial dysfunction in local tissue microenvironment? Front Physiol. (2020) 11:23. doi: 10.3389/fphys.2020.00023

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ma LL, Sun L, Wang YX, Sun BH, and Li YF. Jin YL Association between HO−1 gene promoter polymorphisms and diseases (Review). Mol Med Rep. (2022) 25:29. doi: 10.3892/mmr.2021.12545

PubMed Abstract | Crossref Full Text | Google Scholar

12. Chen YH and Hung SC. Tarng DC Length polymorphism in heme oxygenase-1 and cardiovascular events and mortality in hemodialysis patients. Clin J Am Soc Nephrol. (2013) 8:1756–63. doi: 10.2215/CJN.01110113

PubMed Abstract | Crossref Full Text | Google Scholar

13. Rivera-Valdes JJ, Sifuentes-Franco S, Ramirez-Meza SM, Iniguez-Mosqueda O, Morales-Nunez JJ, Ramirez-Evangelista ML, et al. Association between microsatellite polymorphism in the Heme Oxygenase-1 (HMOX1) gene promoter and type 2 diabetes: an updated meta-analysis. J Diabetes Metab Disord. (2025) 24:63. doi: 10.1007/s40200-025-01575-y

PubMed Abstract | Crossref Full Text | Google Scholar

14. Kaartokallio T, Klemetti MM, Timonen A, Uotila J, Heinonen S, Kajantie E, et al. Microsatellite polymorphism in the heme oxygenase-1 promoter is associated with nonsevere and late-onset preeclampsia. Hypertension. (2014) 64:172–7. doi: 10.1161/HYPERTENSIONAHA.114.03337

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zhou H, Ying X, Liu Y, Ye S, and Yan J. Li Y Genetic polymorphism of heme oxygenase 1 promoter in the occurrence and severity of chronic obstructive pulmonary disease: a meta-analysis. J Cell Mol Med. (2017) 21:894–903. doi: 10.1111/jcmm.13028

PubMed Abstract | Crossref Full Text | Google Scholar

16. Roy-Vallejo E, Fernandez De Cordoba-Onate S, Delgado-Wicke P, Triguero-Martinez A, Montes N, Carracedo-Rodriguez R, et al. Occurrence of SARS-CoV-2 viremia is associated with genetic variants of genes related to COVID-19 pathogenesis. Front Med (Lausanne). (2023) 10:1215246. doi: 10.3389/fmed.2023.1215246

PubMed Abstract | Crossref Full Text | Google Scholar

17. Gill AJ, Garza R, Ambegaokar SS, Gelman BB, and Kolson DL. Heme oxygenase-1 promoter region (GT)n polymorphism associates with increased neuroimmune activation and risk for encephalitis in HIV infection. J Neuroinflammation. (2018) 15:70. doi: 10.1186/s12974-018-1102-z

PubMed Abstract | Crossref Full Text | Google Scholar

18. Zhang R, Hu K, Bai H, Liu H, Pu Y, Yang C, et al. Increased oxidative stress is associated with hyperandrogenemia in polycystic ovary syndrome evidenced by oxidized lipoproteins stimulating rat ovarian androgen synthesis in vitro. Endocrine. (2024) 84:1238–49. doi: 10.1007/s12020-024-03726-2

PubMed Abstract | Crossref Full Text | Google Scholar

19. Armanini D, Boscaro M, Bordin L, and Sabbadin C. Controversies in the pathogenesis, diagnosis and treatment of PCOS: focus on insulin resistance, inflammation, and hyperandrogenism. Int J Mol Sci. (2022) 23:4110. doi: 10.3390/ijms23084110

PubMed Abstract | Crossref Full Text | Google Scholar

20. Sun Y, Li S, Liu H, Bai H, Hu K, Zhang R, et al. Oxidative stress promotes hyperandrogenism by reducing sex hormone-binding globulin in polycystic ovary syndrome. Fertil Steril. (2021) 116:1641–50. doi: 10.1016/j.fertnstert.2021.07.1203

PubMed Abstract | Crossref Full Text | Google Scholar

21. Zhang R, Liu H, Bai H, Zhang Y, Liu Q, Guan L, et al. Oxidative stress status in Chinese women with different clinical phenotypes of polycystic ovary syndrome. Clin Endocrinol (Oxf). (2017) 86:88–96. doi: 10.1111/cen.13171

PubMed Abstract | Crossref Full Text | Google Scholar

22. Liu S, Zhao X, Meng Q, and Li B. Screening of potential biomarkers for polycystic ovary syndrome and identification of expression and immune characteristics. PloS One. (2023) 18:e0293447. doi: 10.1371/journal.pone.0293447

PubMed Abstract | Crossref Full Text | Google Scholar

23. Manneras-Holm L, Benrick A, and Stener-Victorin E. Gene expression in subcutaneous adipose tissue differs in women with polycystic ovary syndrome and controls matched pair-wise for age, body weight, and body mass index. Adipocyte. (2014) 3:190–6. doi: 10.4161/adip.28731

PubMed Abstract | Crossref Full Text | Google Scholar

24. Gao H, Meng J, Xing H, Nie S, Xu M, Zhang S, et al. Association of heme oxygenase-1 with the risk of polycystic ovary syndrome in non-obese women. Hum Reprod. (2014) 29:1058–66. doi: 10.1093/humrep/deu029

PubMed Abstract | Crossref Full Text | Google Scholar

25. Rotterdam EA-SPCWG. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome. Fertil Steril. (2004) 81:19–25. doi: 10.1016/j.fertnstert.2003.10.004

PubMed Abstract | Crossref Full Text | Google Scholar

26. Liu Q, Liu H, Bai H, Huang W, Zhang R, Tan J, et al. Association of SOD2 A16V and PON2 S311C polymorphisms with polycystic ovary syndrome in Chinese women. J Endocrinol Invest. (2019) 42:909–21. doi: 10.1007/s40618-018-0999-5

PubMed Abstract | Crossref Full Text | Google Scholar

27. Zhang J, Zhang Y, Liu H, Bai H, Wang Y, Jiang C, et al. Antioxidant properties of high-density lipoproteins are impaired in women with polycystic ovary syndrome. Fertil Steril. (2015) 103:1346–54. doi: 10.1016/j.fertnstert.2015.02.024

PubMed Abstract | Crossref Full Text | Google Scholar

28. Zhang J, Fan P, Liu H, Bai H, Wang Y, and Zhang F Apolipoprotein A-I. and B levels, dyslipidemia and metabolic syndrome in south-west Chinese women with PCOS. Hum Reprod. (2012) 27:2484–93. doi: 10.1093/humrep/des191

PubMed Abstract | Crossref Full Text | Google Scholar

29. Robert Y, Dubrulle F, Gaillandre L, Ardaens Y, Thomas-Desrousseaux P, Lemaitre L, et al. Ultrasound assessment of ovarian stroma hypertrophy in hyperandrogenism and ovulation disorders: visual analysis versus computerized quantification. Fertil Steril. (1995) 64:307–12. doi: 10.1016/S0015-0282(16)57728-0

PubMed Abstract | Crossref Full Text | Google Scholar

30. Higuchi R. PCR technology. In: Erlich HA, editor. Principles and applications for DNA amplification, 1st edn. Stockton Press, New York (1989). p. 36.

Google Scholar

31. Seu L, Burt TD, Witte JS, Martin JN, Deeks SG, and McCune JM. Variations in the heme oxygenase-1 microsatellite polymorphism are associated with plasma CD14 and viral load in HIV-infected African-Americans. Genes Immun. (2012) 13:258–67. doi: 10.1038/gene.2011.76

PubMed Abstract | Crossref Full Text | Google Scholar

32. Escobar-Morreale HF and Luque-Ramirez M. Role of androgen-mediated enhancement of erythropoiesis in the increased body iron stores of patients with polycystic ovary syndrome. Fertil Steril. (2011) 95:1730–1735.e1731. doi: 10.1016/j.fertnstert.2011.01.038

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ryter SW. Heme oxygenase-1: an anti-inflammatory effector in cardiovascular, lung, and related metabolic disorders. Antioxidants (Basel). (2022) 11:555. doi: 10.3390/antiox11030555

PubMed Abstract | Crossref Full Text | Google Scholar

34. McClung JA, Levy L, Garcia V, Stec DE, Peterson SJ, and Abraham NG. Heme-oxygenase and lipid mediators in obesity and associated cardiometabolic diseases: Therapeutic implications. Pharmacol Ther. (2022) 231:107975. doi: 10.1016/j.pharmthera.2021.107975

PubMed Abstract | Crossref Full Text | Google Scholar

35. Raghunandan S, Ramachandran S, Ke E, Miao Y, Lal R, Chen ZB, et al. Heme oxygenase-1 at the nexus of endothelial cell fate decision under oxidative stress. Front Cell Dev Biol. (2021) 9:702974. doi: 10.3389/fcell.2021.702974

PubMed Abstract | Crossref Full Text | Google Scholar

36. Bao W, Song F, Li X, Rong S, Yang W, Zhang M, et al. Plasma heme oxygenase-1 concentration is elevated in individuals with type 2 diabetes mellitus. PloS One. (2010) 5:e12371. doi: 10.1371/journal.pone.0012371

PubMed Abstract | Crossref Full Text | Google Scholar

37. Kerins MJ and Ooi A. The roles of NRF2 in modulating cellular iron homeostasis. Antioxid Redox Signal. (2018) 29:1756–73. doi: 10.1089/ars.2017.7176

PubMed Abstract | Crossref Full Text | Google Scholar

38. Luque-Ramirez M, Alvarez-Blasco F, Alpanes M, and Escobar-Morreale HF. Role of decreased circulating hepcidin concentrations in the iron excess of women with the polycystic ovary syndrome. J Clin Endocrinol Metab. (2011) 96:846–52. doi: 10.1210/jc.2010-2211

PubMed Abstract | Crossref Full Text | Google Scholar

39. Jagadeesh ASV, Fang X, Kim SH, Guillen-Quispe YN, Zheng J, Surh YJ, et al. Non-canonical vs. Canonical functions of heme oxygenase-1 in cancer. J Cancer Prev. (2022) 27:7–15. doi: 10.15430/JCP.2022.27.1.7

PubMed Abstract | Crossref Full Text | Google Scholar

40. Song J, Zhang X, Liao Z, Liang H, Chu L, Dong W, et al. 14-3-3zeta inhibits heme oxygenase-1 (HO-1) degradation and promotes hepatocellular carcinoma proliferation: involvement of STAT3 signaling. J Exp Clin Cancer Res. (2019) 38:3. doi: 10.1186/s13046-018-1007-9

PubMed Abstract | Crossref Full Text | Google Scholar

41. Hwang HW, Lee JR, Chou KY, Suen CS, Hwang MJ, Chen C, et al. Oligomerization is crucial for the stability and function of heme oxygenase-1 in the endoplasmic reticulum. J Biol Chem. (2009) 284:22672–9. doi: 10.1074/jbc.M109.028001

PubMed Abstract | Crossref Full Text | Google Scholar

42. Lin Q, Weis S, Yang G, Weng YH, Helston R, Rish K, et al. Heme oxygenase-1 protein localizes to the nucleus and activates transcription factors important in oxidative stress. J Biol Chem. (2007) 282:20621–33. doi: 10.1074/jbc.M607954200

PubMed Abstract | Crossref Full Text | Google Scholar

43. Hsu FF, Chiang MT, Li FA, Yeh CT, Lee WH, and Chau LY. Acetylation is essential for nuclear heme oxygenase-1-enhanced tumor growth and invasiveness. Oncogene. (2017) 36:6805–14. doi: 10.1038/onc.2017.294

PubMed Abstract | Crossref Full Text | Google Scholar

44. Lin PH, Lan WM, and Chau LY. TRC8 suppresses tumorigenesis through targeting heme oxygenase-1 for ubiquitination and degradation. Oncogene. (2013) 32:2325–34. doi: 10.1038/onc.2012.244

PubMed Abstract | Crossref Full Text | Google Scholar

45. Salinas M, Wang J, Rosa de Sagarra M, Martin D, Rojo AI, Martin-Perez J, et al. Protein kinase Akt/PKB phosphorylates heme oxygenase-1 in vitro and in vivo. FEBS Lett. (2004) 578:90–4. doi: 10.1016/j.febslet.2004.10.077

PubMed Abstract | Crossref Full Text | Google Scholar

46. Ellakany WI, Mahmoud MoheyEldin K, Invernizzi P, Mahmoud ElKady A, Eldin Fathy Abou Elkheir H, Abdel Haleem Abo Elwafa R, et al. Study of the influence of heme oxygenase 1 gene single nucleotide polymorphism (rs2071746) on esophageal varices among patients with cirrhosis. Eur J Gastroenterol Hepatol. (2018) 30:888–92. doi: 10.1097/MEG.0000000000001161

PubMed Abstract | Crossref Full Text | Google Scholar

47. Ono K, Mannami T, and Iwai N. Association of a promoter variant of the haeme oxygenase-1 gene with hypertension in women. J Hypertens. (2003) 21:1497–503. doi: 10.1097/00004872-200308000-00013

PubMed Abstract | Crossref Full Text | Google Scholar

48. Ono K, Goto Y, Takagi S, Baba S, Tago N, Nonogi H, et al. A promoter variant of the heme oxygenase-1 gene may reduce the incidence of ischemic heart disease in Japanese. Atherosclerosis. (2004) 173:315–9. doi: 10.1016/j.atherosclerosis.2003.11.021

PubMed Abstract | Crossref Full Text | Google Scholar

49. Lee EY, Lee YH, Kim SH, Jung KS, Kwon O, Kim BS, et al. Association between heme oxygenase-1 promoter polymorphisms and the development of albuminuria in type 2 diabetes: A case-control study. Med (Baltimore). (2015) 94:e1825. doi: 10.1097/MD.0000000000001825

PubMed Abstract | Crossref Full Text | Google Scholar

50. Rueda B, Oliver J, Robledo G, Lopez-Nevot MA, Balsa A, Pascual-Salcedo D, et al. HO-1 promoter polymorphism associated with rheumatoid arthritis. Arthritis Rheumatol. (2007) 56:3953–8. doi: 10.1002/art.23048

PubMed Abstract | Crossref Full Text | Google Scholar

51. Okamoto I, Krogler J, Endler G, Kaufmann S, Mustafa S, Exner M, et al. A microsatellite polymorphism in the heme oxygenase-1 gene promoter is associated with risk for melanoma. Int J Cancer. (2006) 119:1312–5. doi: 10.1002/ijc.21937

PubMed Abstract | Crossref Full Text | Google Scholar

52. Wu MM, Chiou HY, Chen CL, Hsu LI, Lien LM, Wang CH, et al. Association of heme oxygenase-1 GT-repeat polymorphism with blood pressure phenotypes and its relevance to future cardiovascular mortality risk: an observation based on arsenic-exposed individuals. Atherosclerosis. (2011) 219:704–8. doi: 10.1016/j.atherosclerosis.2011.08.047

PubMed Abstract | Crossref Full Text | Google Scholar

53. Kramer M, Sponholz C, Slaba M, Wissuwa B, Claus RA, Menzel U, et al. Alternative 5' untranslated regions are involved in expression regulation of human heme oxygenase-1. PloS One. (2013) 8:e77224. doi: 10.1371/journal.pone.0077224

PubMed Abstract | Crossref Full Text | Google Scholar