The relationship between genetic variants associated with primary ovarian insufficiency and lipid profile in women recruited from MASHAD cohort study.

BACKGROUND AND AIM: Primary Ovarian Insufficiency (POI) is defined by the occurrence of menopause before the age of 40 years. It is often associated with cardiovascular disease (CVD). The purpose of this study was to explore the relationship between POI-associated genotypes cardiometabolic disorder risk factors.
METHODS: One hundred seventeen women with POI and one hundred eighty-three healthy women without POI were recruited in this study. DNA was extracted and analyzed using ASO-PCR or Tetra ARMS-PCR. Lipid profiles were also assessed.
RESULTS: Multivariate logistic regression analysis showed that individuals with GG vs. TT genotype of the rs1046089 SNP were more likely to have a higher serum LDL (p = 0.03) compared to the control group. There was also a significant association between low serum HDL and rs2303369 and rs4806660 SNP genotypes in the POI group. In the POI group, the percentage of those with high total cholesterol was lower in those with a CC genotype compared to those with a TT genotype (p = 0.03).
CONCLUSION: Some SNPs reported to be associated with POI appear to be independently associated with dyslipidemia. These results may be helpful to identify subjects with POI who may be susceptible to CVD.

Introduction

Menopause is defined as occurring when menstruation stops for 12 consecutive months due to loss of ovarian follicular activity; this usually occurs around the age of 45–55 years [1]. Menopause that happens before the age of 40 years is called Primary Ovarian Insufficiency (POI) which may be natural or due to or reproductive surgeries [2-4]. This condition is accompanied by amenorrhea, estrogen deficiency and an increase of gonadotrophin levels [5]. It has been reported that about 3.6% of women develop POI [6]. Various factors such as smoking [7], certain medications [8], infections and genetic and autoimmune disorders have been associated with POI [9]. Iatrogenic causes such as radiotherapy, chemotherapy, pelvic surgery are also associated with POI [5, 9–14]. About 60–90% of POI cases are idiopathic [14]. Several studies have also found that POI increases the risk of hypertension [15], cardiovascular disease [3], osteoporosis [16] cerebral infarction [17], all-cause mortality [18] type 2 diabetes mellitus [15], and other negative health consequences [16].

Genetic factors also play a significant role in POI [19]. The heritability of menopausal age is estimated about 30–85% [20]. About 20–25% of POI cases are due to genetic causes [14]. Genome-wide association studies (GWAS) have identified a polymorphism (rs16991615) of minichromosome maintenance 8 homologous recombination repair factor (MCM8) gene involved in the age of natural menopause [21]. Also, rs1046089 and rs4806660, located on Proline Rich Coiled-Coil 2A (PRRC2A) and transmembrane (TMEM) gene, respectively are associated with the age at menopause [22, 23]. GWAS identified several other variants that are associated with POI [24].

Deleterious changes in risk factors for cardiometabolic disorders often occur around the age of menopausal [18, 25–27]. Estrogen is involved in dilating blood vessels and helping blood flow [18]. Various studies have also shown that estrogen therapy in postmenopausal women reduces serum total cholesterol and low density lipoprotein (LDL) cholesterol concentrations, and increases serum high-density lipoprotein (HDL) cholesterol and triglyceride concentrations [28, 29]. Moreover, lack of ovarian function in the menopause is involved in the activation of the renin-angiotensin system, leads to immunodeficiency, inflammation and endothelial dysfunction [18, 30]. These are associated with obesity, diabetes and high blood pressure [18, 30]. Several studies have shown that age at menopause is associated with cardiovascular disease [18, 31–33]. A Japanese study found women with early menopause had a higher risk for hypercholesterolemia [34], and another study showed that early menopause is associated with hypertension [35]. Sarnowski et al. founded that genetic variants associated with early menopause are also associated with increased cardiovascular disease risk [36]. There appeared to be a need to evaluate the relationship between POI-related variants with lipid profile and susceptibility to cardiometabolic disease risk factors. No previous study has investigated association between polymorphisms related to early menopause with lipid profiles.

Therefore, we aimed to explore the associations between POI-related variants with lipid profile and susceptibility to cardiometabolic disease risk factors in Mashhad stroke and heart atherosclerotic disorders (MASHAD) cohort study population.

Methods

Totally, 117 women who had POI were included in the case group. Healthy women (n = 183) were recruited into the control group. All of participants were recruited as part of the MASHAD study. The MASHAD study was a cohort study from 2010–2020 that were included 9704 participants aged 35–65 years who will be follow-up exams every three years until 2020 [37]. The inclusion criteria were as follows: the diagnostic criterion was based on the definition of POI: (1) Women who go through menopause before the age of 40 years; (2) 12 continuous months have passed since the last bleeding; and (3) serum FSH > 40 IU/L. Exclusion criteria were: women over 40 years old, with a history of diseases and surgeries affecting menstruation (oophorectomy, hysterectomy), history of genetic abnormalities and syndromes associated with an early menopause is a part of their manifestation, history of using drugs affecting menstruation. Blood samples were collected into vacutainer® plain tubes, and were taken after 12 h fast. Blood samples were centrifuged at 4° C in 5000 rpm for 15 min and the serum part was used for lipid profile measurement. Body Mass Index (BMI) was measured using standard method [37]. Kidney, liver, and thyroid activity were normal in all participants.

DNA extraction and quality controls

Participants' DNA was extracted from 200 µl blood or buffy coat samples using a DNA extraction kit (Pars Tous, Mashhad, Iran). Qualitative and quantitative quality control was performed by agarose gel electrophoresis (Pars Tous, Mashhad, Iran) and Nano drop 2000 (Thermo Fisher Scientific, USA) in 280 and 260 nm wavelengths, respectively.

Allele-specific oligonucleotide polymerase chain reaction (ASO-PCR)

The ASO-PCR reaction volume was 15 µl which included: 1.5 µl water, 2 µl genomic DNA, 1 µl of each primer, and 7.5 µl master mix (Pars Tous, Mashhad, Iran). First, to carry out PCR, we performed one cycle of denaturation for 7 min at 95 °C. After that, 35 cycles include the following: 95 °C for 30 s, annealing for 30 s at 60 °C, 72 °C for 30 s, and eventually one cycle of 7 min was done for final extension.

Tetra amplification refractory mutation system PCR (ARMS-PCR)

Tetra ARMS was carried out by the same method and the same composition of 15 µl reaction volume that performed in ASO-PCR. Primers were designed with Primer1 software [38].

Lipid profile measurements and dyslipidemia diagnosis

Total cholesterol (TC), triglycerides (TG), and high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C) levels were measured from serum which was taken from participants 12 h fasting using standard method [39, 40]. The NCEP ATPIII criteria [41] were used to diagnose dyslipidemia. (1) If serum cholesterol levels ≥ 200 mg /dl (5.2 mmol/l) is considered hypercholesterolemia. (2) If HDL cholesterol levels < 50 mg/dl (< 1.3 mmol/l) for women is considered low HDL cholesterol. (3) If LDL cholesterol levels ≥ 130 mg/dl (≥ 3.4 mmol/l) is considered high LDL cholesterol. (4) If serum triglycerides ≥ 150 mg/dl is considered Hypertriglyceridemia.

Physical activity level assessment

The equations of James and Schofield for energy requirements, were used to assess physical activity of all participants. Questions regarding the physical activity level were based on the mentioned equations which were selected from World Health Organization MONICA project questionnaires [42]. The level of physical activity was calculated by total energy expenditure (TTE) and basal metabolism rate (BMR) during a whole day and night [43].

Ethics

All steps of the study were approved by the Mashhad University of Medical Sciences (MUMS) Ethics Committee. Informed consents were obtained from all subjects and the procedure and possible risks were explained completely, pursuant to the Declaration of Helsinki [44].

Statistical analysis

All analysis tests were performed by Statistical Package for Social Sciences (SPSS) (IBM Corp. Released 2016. IBM SPSS Statistics for Windows, Version 24.0. Armonk, NY: IBM Corp.). The values in this study have been reported frequently with percentage or mean and standard deviation. Single Nucleotide Polymorphisms' (SNPs') genotypes between POI cases and healthy controls were compared by Chi-square test. Assessment of normal distribution in quantitative data was performed by Kolmogorov-Smirnoff test. Man-Whitney test was used for comparing quantity of normally distributed values between the subgroups. In this study with having at least 117 subjects in each group and confidence interval 95% and power 80% the value of the standardized effect size about 37% reported which has a good precision. Besides, we used multivariate logistic regression to prevent confounder's factors from affecting our results. P-value < 0.05 was statistically significant.

Results

The clinical characteristics of the population have been summarized in Table 1. In our cross-sectional analysis, participants had a mean age of 55 years averagely. The mean age of the case and control group were 55.21 ± 5.56 and 54.62 ± 2.89, respectively. There was not significantly difference (P = 0.4). The control-women were matched by age. The details on the genotype distribution of these SNPs between the two groups was investigated in our previous studies [45]. All of the polymorphisms were in Hardy–Weinberg equilibrium.

Table 1

Demographic features and characteristics of the study population

Characteristics	POI cases, N = 117	Controls, N = 183	P-value
Age (Y)	55.21 ± 5.56	54.62 ± 2.89	0.4
BMI (kg/m2)	28.78 ± 5.06	29.34 ± 4.22	0.323
FBG (mg/dl)	88.74 ± 20.06	88.96 ± 31.43	0.948
PAL	1.78 ± 0.25	1.70 ± 0.23	0.008
Non smoker	86 (74.8%)	146 (79.8%)	0.390
Ex-smoker	10 (8.7%)	9 (4.9%)
Current smoker	19 (16.5%)	28 (15.3%)
TC (mg/dl)	207.8 ± 35.8	188.8 ± 33.1	< 0.001
TG (mg/dl)	111.0 (82.0–162.5)	117.0 (80.0–159.0)	0.684
LDL-C (mg/dl)	128.39 ± 33.96	110.47 ± 33.42	< 0.001
HDL-C (mg/dl)	47.45 ± 9.33	45.06 ± 11.34	0.040

Statistically significant values are shown in bold

Data are shown as Mean ± SD; Student t-test and Chi-square test were used; POI: Primary Ovarian Insufficiency; SD: Standard Deviation; BMI: Body Mass Index; FBG: Fasting Blood Glucose; PAL: Physical Activity Level; TC: Total Cholesterol; TG: Triglyceride; HDL: High Density Lipoprotein; LDL: Low Density Lipoprotein

Differences in level of lipid serum between different genotypes of POI-related polymorphisms were examined (Table 2). In the POI cases, serum total cholesterol was significantly different between various rs16991615, rs244715, rs4806660, and rs10183486 SNP genotypes; however, this difference was not observed in the healthy controls. Interestingly, three of four of the investigated factors including serum total cholesterol, triglyceride, and HDL were substantially difference in different genotypes of rs4806660 SNP in participants with POI and this association was not detected in controls group.

Table 2

Association of genotypes related to POI with lipid profile in target population divided into POI cases and healthy controls

Polymorphisms	POI cases (N = 117)					Controls (N = 183)
rs16991615	AA	GA	GG	P	AA		GA	GG	P
Serum total cholesterol (mg/dl)	222.3 ± 41.4a	195.8 ± 35.3b	209.4 ± 29.5ab	0.01	190.2 ± 33.8		191.2 ± 33.5	186.4 ± 32.7	0.64
Triglyceride (mg/dl)	108.0 (78.0–181.0)	112.0 (82.0–147.0)	107.0 (89.0–156.0)	0.96	98.5 (65.8–121.3)		129.0 (87.0–159.0)	108.0 (79.5–162.3)	0.14
HDL (mg/dl)	48.6 ± 7.7	45.9 ± 8.9	48.0 ± 10.5	0.44	45.0 (40.8–53.0)		44.0 (35.0–53.0)	43.0 (35.6–54.3)	0.58
LDL (mg/dl)	136.1 ± 37.8	120.3 ± 29.9	130.6 ± 34.1	0.15	116.7 ± 35.2		112.7 ± 34.0	107.1 ± 32.5	0.39

Statistically significant values are shown in bold

POI: Primary Ovarian Insufficiency; HDL: High-Density Lipoprotein; LDL: Low Density Lipoprotein; One-way analysis of variance (ANOVA) OR Kruskal–Wallis tests were used;  Significant P value < 0.05

140.5

(77.5–202.7)

(85.0–179.0)

(82.7–151.5)

(84.0–231.3)

122.0

(85.5–163.5)

113.0

(77.0–149.5)

43.5

(33.8–49.3)

42.0

(36.3–55.7)

(35.6–52.9)

126.8

(111.9–158.8)

134.4

(115.9–146.5)

124.1

(93.1–142.6)

201.0

(191.0–236.0)

207.5

(178.0–231.7)

208.5

(181.2–223.2)

194.5

(174.0–211.5)

181.5

(164.5–201.5)

192.0

(162.0–217.0)

120.0

(82.0–185.0)

111.5

(78.3–153.0)

107.5

(84.5–151.5)

106.0

(88.3–144.0)

124.0

(77.3–166.5)

115.0

(80.0–154.0)

45.5

(39.3–54.4)

43.5

(34.3–54.5)

42.8

(36.6–52.1)

196.0

(177.0–225.0)

207.0

(184.0–232.0)

209.0

(186.5–224.5)

(79.0–214.0)

(82.0–153.0)

109.0

(84.0–154.5)

(80.0–155.0)

118.0

(84.0–162.5)

117.0

(76.0–159.0)

52.0

(40.0–58.0)a

(34.0–55.0)ab

(35.1–50.0)b

75.5

(58.6–111.0)a

(87.1–137.0)b

122.5

(101.1–138.6)bc

201.5

(196.3–233.5)

200.5

(176.0–225.5)

213.0

(200.0–234.0)

101.5

(73.3–167.0)

(81.3–150.7)

118.0

(94.0–186.0)

145.0

(93.5–164.3)

118.5

(84.7–152.3)

(77.0–161.0)

53.5

(41.5–57.7)

(35.0–52.5)

42.8

(36.0–53.1)

93.5

(73.3–168.0)ab

122.0

(97.5–183.5)a

100.0

(73.0–140.0)b

133.0

(100.0–167.0)

111.5

(80.7–160.5)

117.0

(76.7–156.0)

46.0

(40.0–56.0)

43.7

(35.5–54.3)

43.2

(35.7–52.5)

97.5

(74.7–119.5)

107.5

(80.7–150.7)

122.0

(94.0–185.0)

143.0

(91.0–180.5)

111.0

(77.5–144.0)

120.0

(82.5–164.5)

35.5

(31.4–39.3)a

(39.0–54.5)b

43.5

(34.0–53.3)bc

98.0

(69.0–129.0)

117.0

(94.5–182.0)

106.0

(77.5–152.0)

106.0

(75.0–150.0)

120.0

(84.5–162.3)

117.0

(76.0–154.5)

56.5

(49.8–60.8)a

44.5

(34.0–55.0)b

42.8

(35.3–48.7)bc

117.6

(100.9–139.6)

134.4

(116.3–149.0)

123.6

(103.1–147.8)

Furthermore, Tables 3 and 4 show the results of multivariate logistic regression analysis before and after adjustment for age and physical activity level. These results for rs1046089 showed individuals with GG genotype were more likely to have low LDL risk (OR = 5.48, CI = 1.14–26.34, p = 0.02) than individuals with the AA genotype in control group using a multivariate logistic regression test. Also, the results demonstrate that there was a significant association with high TG risk in CC variant vs. TT in rs10183486 in POI group (OR = 4.63, CI = 1.17–18.29, p = 0.02). Furthermore, these results suggest that individuals carried the recessive homozygous genotype (CC) of rs2303369 SNP compared to individuals with dominant homozygous genotype (TT) had an increased risk of a high serum LDL, low serum HDL and the risk of high TC was decreased in control group.

Table 3

Investigation of the relationship between the risk of polymorphisms related to Primary Ovarian Insufficiency and dyslipidemia factors in the study population in a non-adjusted model

PolymorphismsLipid profile abnormalities	POI cases				Controls
	OR (95% CI)	P	OR (95% CI)	P	OR (95% CI)	P	OR (95% CI)	P
rs16991615	GA/AA		GG/AA		GA/AA		GG/AA
High LDL (mg/dl)	0.65 (0.24, 1.74)	0.39	1.15 (0.45,2.97)	0.77	0.65 (0.23, 1.84)	0.42	0.46 (0.16,1.30)	0.14
High Triglyceride (mg/dl)	0.51 (0.17, 1.50)	0.22	0.87 (0.32, 2.35)	0.79	1.93 (0.51, 7.32)	0.33	2.41 (0.64, 9.02)	0.19
Low HDL (mg/dl)	1.43 (0.52, 3.89)	0.48	1.08 (0.41, 2.80)	0.87	0.96 (0.32, 2.85)	0.94	0.88 (0.30, 2.59)	0.82
High total cholesterol (mg/dl)	0.36 (0.12, 1.02)	0.054	0.89 (0.32, 2.51)	0.83	0.64 (0.23, 1.83)	0.41	0.63 (0.22, 1.78)	0.39

Binary logistic regression models were performed; POI: Primary Ovarian Insufficiency; LDL: Low-Density Lipoprotein; HDL: High-Density Lipoprotein; High LDL > 130 and Normal LDL < 130; High Triglyceride > 150 and Normal Triglyceride < 150; Low HDL < 40 (Male) or < 50 (Female) and Normal HDL > 40 (Male) or > 50 (Female), High Total Cholesterol > 200 and Normal Total Cholesterol < 200; Normal groups of all lipid profile measures were considered as references; The common genotypes of each studied variants were considered as references; Significant P-value < 0.05

Table 4

Investigation of the relationship between the risk of polymorphisms related to Primary Ovarian Insufficiency and dyslipidemia factors in the study population in an adjusted model for age and physical activity level

PolymorphismsLipid profile abnormalities	POI cases				Control
	OR (95% CI)	P	OR (95% CI)	P	OR (95% CI)	P	OR (95% CI)	P
rs16991615	GA/AA		GG/AA		GA/AA		GG/AA
High LDL (mg/dl)	0.63 (0.23, 1.72)	0.36	1.21 (0.46, 3.17)	0.69	0.67 (0.22, 2.04)	0.48	0.37 (0.12, 1.15)	0.08
High Triglyceride (mg/dl)	0.45 (0.15, 1.38)	0.16	0.84 (0.30, 2.33)	0.74	1.38 (0.34, 5.53)	0.64	1.69 (0.43, 6.61)	0.45
Low HDL (mg/dl)	1.37 (0.50, 3.76)	0.54	1.05 (0.40, 2.76)	0.91	0.63 (0.19, 2.05)	0.45	0.59 (0.18, 1.89)	0.38
High Total Cholesterol (mg/dl)	0.34 (0.12, 0.99)	0.04	0.91 (0.32, 2.57)	0.86	0.67 (0.22, 2.02)	0.48	0.55 (0.18, 1.65)	0.29

Statistically significant values are shown in bold

Binary logistic regression models were performed; POI: Primary Ovarian Insufficiency; LDL: Low-Density Lipoprotein; HDL: High-Density Lipoprotein; High LDL > 130 and Normal LDL < 130; High Triglyceride > 150 and Normal Triglyceride < 150; Low HDL < 40 (Male) or < 50 (Female) and Normal HDL > 40 (Male) or > 50 (Female), High Total Cholesterol > 200 and Normal Total Cholesterol < 200; Normal groups of all lipid profile measures were considered as references; The common genotypes of each studied variants were considered as references; Adjusted with Age & Physical Activity Level; Significant P-value < 0.05

The risk of low HDL was increased in individuals carrying rs23303369 variant (CT) compared to non-carriers (TT) in both studied groups, (OR, 6.6; 95% CI, 1.88–19.53, P = 0.003 in POI group) and (OR = 4.21, CI = 1.17–15.18, p = 0.02 in control group). Table 5 shows the association of genotypes related to POI with lipid profile in the POI cases and healthy controls.

Table 5

Association of different genotypes related to POI with lipid profile in target population divided into POI cases and healthy controls

Polymorphisms
rs16991615	AA			GA			GG
	POI cases (N = 48)	Controls (N = 86)	P	POI cases (N = 40)	Controls (N = 79)	P	POI cases (N = 29)	Controls (N = 18)	P
Serum total cholesterol (mg/dl)	207.00 (50.00)	190.00 (57.00)	0.020	195.79 ± 35.28	191.16 ± 33.51	0.489	209.40 ± 29.54	186.42 ± 32.73	< 0.001
Triglyceride (mg/dl)	108.00 (103.00)	98.50 (56.00)	0.187	112.00 (65.00)	129.00 (72.00)	0.286	107.00 (67.00)	108.00 (83.00)	0.785
HDL (mg/dl)	48.62 ± 7.73	46.92 ± 8.39	0.489	45.94 ± 8.94	44.59 ± 12.12	0.538	48.00 (13.80)	43.00 (18.68)	0.120
LDL (mg/dl)	136.10 ± 37.79	116.74 ± 35.22	0.091	120.32 ± 29.94	112.70 ± 34.01	0.236	130.64 ± 34.13	107.10 ± 32.51	< 0.001

Statistically significant values are shown in bold

POI: Primary Ovarian Insufficiency; HDL: High-Density Lipoprotein; LDL: Low Density Lipoprotein; T-Student Test or Mann Whitney Test were used; Significant P value < 0.05

Discussion

Our findings suggest that serum levels of several parameters in the fasted serum lipid profile consisting of total cholesterol, LDL, and HDL, but not serum triglycerides, were associated with POI. Further analyses indicated that genotypes of polymorphisms, which were previously reported to be related with the incidence of POI, are substantially associated with the level of lipid profile factors in POI cases. Moreover, we found that some genotypes in specific polymorphisms including rs4806660, rs10183486, and rs2303369 SNPs were significantly related to abnormalities regarding total cholesterol, LDL and HDL level in the cases' serum.

Initially our results found significant difference between POI cases and control participants for serum HDL, LDL, and total cholesterol levels, and these factors were significantly higher in the POI group. Gulhan and his coworkers have also found that among the 4 lipid factors, only total cholesterol and LDL were substantially different between cases diagnosed with premature ovarian failure (POF) and control subjects. POF group in this study had higher levels of total cholesterol and LDL [46]. Gulhan et al. work had included only women with previous history of successful childbirth and without any hormone therapy within the last 6 months, this inconsistency in HDL serum level had happened. However, this difference might have occurred due to their use of a small sample size as it was as one third as our study or due to the role of age as a confounding factor and also the role of some genetic variants related to POI. Interestingly, a recent study on 3 Dutch university medical centers did not report any significant difference in lipid profile between previously POI-diagnosed participants and population-based controls [47]. In their study, secondary amenorrhea (cessation of menstruation for at least 3 consecutive months) was one of the criteria for including POI cases; While, this period was too short compared to our criteria (12 consecutive months) and this important thing might have affected their results. Moreover, they have not indicated whether their participants had any previous history of surgeries, diseases, or taking medications related to female reproductive tract or not; Thus, this factor might have not been considered in their study which cause this disagreement.

Knauff et al. have reported changes in lipid profile of cases with POF compared to the controls is potentially related to the rate of ovarian function [48]. Another study which included cases who enter menopause by surgical ovariectomy, clarified that this intervention on female's reproductive tract, has caused impaired lipid metabolism 6 months post-surgery and substantial increase in all four lipid indicators (Total Cholesterol, Triglyceride, LDL, and HDL) [49]. Moreover it has been indicated that Free Androgen Index (FAI) was significantly lower in POI cases than people with regular menstrual cycles; while, it had been previously suggested that higher FAI was associated with high serum triglyceride and LDL in POI cases [50]. Finally, these results suggest that impairment in sexual hormones level, as a result of decreased females' reproductive system activity, is related to weaken lipid metabolism.

Our study has for the first time investigated the association between POI-related polymorphisms' genotypes and lipid profile status. In our analysis, we observed that genotypes of 2 different SNPs including rs4806660 (TMEM150B) and rs10183486 (TLK1) are substantially associated with abnormalities of lipid profile. Regarding the first one, low level of HDL was observed by 5.26 times higher in POI cases with TC genotype rather than CC participants. This association was not found in the controls. It is possible that these results are due to hormonal disorders that occur in postmenopausal individuals and consequently the lipid profile in individuals is impaired. Furthermore, for the second mentioned SNP, hypercholesterolemia was found over 450 percent higher in cases who carried CC than those with TT genotype. Previous studies have reported that both rs4806660 [51] and rs10183486 [22] polymorphisms are associated with the incidence of POI. While it has been proved that TLK1 gene encodes nuclear serine-threonine kinases [22] in which actively transfers signals from receptors of estrogen within the cell membrane to the nucleus [52], the exact function of TMEM150B gene product, called transmembrane protein 150B, has not yet been discussed.

Several SNPs, have been found to be associated with POI, a condition in which ovaries stop releasing sexual hormones, especially estrogen, making females infertile as early as before the age of 40. This estrogen hormone deficiency affects metabolism of lipids and thus, POI cases might face some lipid profile abnormalities that could increase the risk of cardiovascular disorders in participants of our study.

Our study was mainly limited by its design as a cross-sectional study, and changes in the level of lipid profile factors were not prospectively assessed. Data regarding pattern of the target population diet for the consumption of oils and meats, as the most common lipids, were not recorded in this study. We suggest future studies consider these limitations to achieve more reliable results. Moreover, several POI-related SNPs were not investigated in our study that might be associated with the cases' metabolic status and it is highly recommended to include these polymorphisms in the analysis.

In conclusion, SNPs which are previously found to be attributed to POI, could cause impairment in cases' lipid profile status through hormonal abnormalities. Present study clarifies that TC, CC and, CT genotypes of rs4806660, rs10183486, and rs2303369 SNPs, respectively increase the rate of dyslipidemia by approximately 5 times compared to their reference genotype (rs4806660: CC; rs10183486: TT; rs2303369: TT) in POI cases. But, GG and CC genotypes of rs1046089 and rs2303369 SNPs, respectively increase the rate of dyslipidemia about 6 times compared to their reference genotype in healthy population.