Single nucleotide polymorphisms in binding site of miRNA-135a and targeted gene IRS2 are correlated with multiple clinical features of PCOS: A study in Chinese women
Abstract
BACKGROUND:
The etiology of polycystic ovary syndrome (PCOS) remains unclear with highly heterogeneous clinical manifestations, recently growing evidence revealing genetic variants play a crucial part in its pathogenesis.
OBJECTIVE:
This study aimed to examine the correlation between SNPs in miRNA-135a’s binding site of targeted gene IRS2 and clinical manifestations of PCOS in Chinese females.
METHOD:
A total of 126 Chinese women with PCOS and 109 healthy women were enrolled, divided into 4 groups based on different clinical features of hyperandrogenemia (HA), insulin resistance (IR), polycystic ovary morphology (PCOM) and obesity. We analyzed 2 single nucleotide polymorphisms (SNPs) of the IRS2 gene (rs2289046 and rs1865434) and clinical features’ laboratory measurements such as sex hormone, fasting plasma glucose (FPG), fasting plasma insulin (FINS).
RESULTS:
Located in miRNA-135a binding site of IRS2 gene, the rs2289046’s triple genotypes distribution showed a significant difference between PCOS/control group and PCOM/non-PCOM group (
CONCLUSION:
The results revealed the two SNPs as rs2289046 and rs1865434 in the IRS-2 binding region of miRNA-135a have correlations with the clinical features of PCOS in Chinese population.
1.Introduction
Polycystic ovary syndrome (PCOS) is a common female endocrine and metabolic disorder that affects 5–10% women in childbearing age [1] and is characterized by menstrual irregularities, chronic anovulation, obesity, polycystic ovary morphology (PCOM), hyperandrogenemia, insulin resistance (IR) or hyperinsulinemia (HA) [2, 3, 4, 5]. PCOS is considered to be a polygenic feature under the influence of environmental factors, which may be ascribed with the interaction of protective and susceptible genomic variants [6]. Although the etiology of PCOS remains unknown, growing evidences showed the point of genetic defects played a crucial role in PCOS pathogenesis.
Insulin resistance (IR) is reported as one of the main pathophysiological characteristics of reproductive and metabolic disorders in PCOS [7]. The IRS-2 gene is located on chromosome 13q34, encoding the insulin receptor substrate protein 2 (IRS-2), which is a cytoplasmic signaling molecule composed of 1354 amino acids. As a molecular adaptor, IRS-2 regulates peripheral glucose metabolism and pancreatic islet
Single nucleotide polymorphisms (SNPs) in the miRNA-binding sites of target gene are correlated with susceptibility to multiple diseases as multiple miRNA expression profiles are altered in pathophysiology of PCOS [10]. Studies showed that the expression of miR-135a was up-regulated in PCOS’s follicular fluid, while the potential target gene IRS-2 was down-regulated [10, 11]; miR-135a down-regulates the expression of the IRS2 product by binding to the 3’UTR site to play an crucial role in the mechanism of IR. In this study, we explored the relevance of two SNPs rs2289046 and rs1865434 in miRNA-135a and its target gene IRS-2 binding region with the clinical features of PCOS.
2.Materials and methods
2.1Study subjects
Study participants were recruited from the Guangdong Provincial Family Planning Specialist Hospital from June 2014 to December 2016. This study was approved by the Medical Ethics Committee of the Guangdong Provincial Family Planning Specialist Hospital. All participants obtained the informed consent before taken part in this study.
All participants were assigned to the case group and the control group. A total of 126 patients with a diagnosis of PCOS were enrolled in the study. Diagnosis of PCOS was based on the “2003 ESHRE/ASRM guidelines” if more than two of the following criteria were fulfilled: oligo-/anovulation (AO), polycystic ovaries on ultrasonography (USG), and clinical and/or biochemical signs of hyperandrogenism [12]. Disorders with similar clinical features, e.g., congenital adrenal hyperplasia, hypercortisolism, or androgen-secreting tumors, were not considered. Subjects who had received drugs which have known to interfere with hormonal levels within 3 months were also excluded. The control group consisted of 109 healthy women between 20–35 years old with a regular menstrual cycle range of 21–35 days and at least one successful pregnancy history, while no family history of diabetes, endocrine disorders and sonographic signs of PCOS were shown.
2.2Parameters analyzed
Hormone level and metabolic examinations follicle-stimulating hormone (FSH), luteinizing hormone (LH), the LH/FSH ratio, estradiol (E2), prolactin (PRL) and total testosterone were measured by chemiluminescence immunoassay. Before ovarian stimulation, all participants were given a 75 g oral glucose tolerance test and measured the insulin levels. Plasma glucose and insulin were measured by the glucose oxidase method.
In addition, the general information of patients were also collected, i.e., age, body height (m), weight (kg), waist circumference, body mass index (BMIe) calculated as weight (kg)/body height
2.3SNP genotyping
Peripheral blood (2 ml) was collected from all participants. Genomic DNA was isolated by using a TIANamp Genomic DNA kit (TIANGEN Biotech (Beijing) Co., Ltd.). The primer design synthesis and SNP typing were performed with the assistance of Shanghai Jereh Bioengineering Co. The multiple single-base extension (SNapShot) method was used to analyze the polymorphisms of the SNPs. Two SNP sites for the IRS-2 gene, rs2289046 and rs1865434 were analyzed. Total volume of the amplification reaction was set for 15
Table 1
SNPs | Nucleotide change | Primers | Product size (bp) |
---|---|---|---|
rs2289046 | G | Forward: TACCTGCGATGTTTACGTCCAC | 198 |
Reverse: TATTCATCCCCTTCCCAAAGC | |||
rs1865434 | A | Forward: ACTCCAGAGATTGCTCTGTTC | 191 |
Reverse: ACACAGTCATTGCTCAGATCC |
2.4Statistical analysis
The database was created using EpiData 3.0 software. Analyses were performed using SPSS Statistics version 22.0. The SNPs were analyzed for deviation from Hardy-Weinberg equilibrium (HWE). We calculated the genotype and allelic frequencies of target SNPs. The
3.Results
3.1Basic features of the participants in this research
Clinical, hormonal and demographic features of women with PCOS and controls are summarized in Table 2. The demographic features of cases and controls were comparable, showing a statistically significant difference (
Table 2
Parameters | Controls ( | PCOS ( |
| |||
---|---|---|---|---|---|---|
Age (years) | 30.97 | 28.61 | 3.77 | |||
Weight (kg) | 54.73 | 59.34 | 0.01 | |||
BMI (kg/m | 19.18 | 23.50 | 0.06 | |||
FSH (mIU/ml) | 4.70 | 5.15 | ||||
LH (mIU/ml) | 3.83 | 7.88 | ||||
LH/FSH | 0.85 | 1.74 | ||||
E2 (pmmol/L) | 144.95 | 167.39 | 0.36 | |||
PRL (mIU/L) | 328.18 | 339.54 | 0.45 | |||
T (nmol/L) | 0.95 | 2.38 | ||||
HOMA_IR | 1.60 | 4.27 |
3.2Genotype frequencies of IRS-2 gene
The distribution of genotypes at the genetic loci IRS-2 rs2289046 and rs1865434 are in accordance with HWE in both control and PCOS samples (
Table 3
SNP | Genotype | PCOS | Control |
|
| OR (95% CI) | ||
rs2289046 | G/A | |||||||
GG | 29 | (23) | 10 | (9.2) | 1 | |||
AG | 78 | (61.9) | 77 | (70.6) | 6.339 | 0.012 | 0.349 (0.159–0.766) | |
AA | 19 | (15.1) | 22 | (20.2) | 5.422 | 0.02 | 0.298 (0.116–0.766) | |
rs1865434 | A/G | |||||||
AA | 76 | (60.3) | 64 | (58.7) | 1 | |||
AG | 45 | (35.7) | 45 | (41.3) | 0.404 | 0.417 | 0.842 (0.495–1.431) | |
GG | 5 | (4) | 0 | 4.092 | 0 | 1.066 (1.008–1.127) |
OR: odds ratio; CI: confidence intervals;
3.3Genotype frequencies distribution of IRS-2 gene in accordance with clinical manifestations of PCOS
Based on different phenotypes and features of PCOS, a total of 235 females were divided into several subgroups: PCO morphology (PCOM) and non-PCOM on pelvic ultrasound; HA and non-HA; IR and non-IR; obesity and non-obesity on BMI.
3.3.1. Polycystic ovaries were diagnosed according to the Rotterdam criteria with presenting follicles diameter and ovarian volume in sonography. A total of 235 subjects were divided into PCOM group and non-PCOM group. The two SNPs (rs2289046 and rs1865434) genotype distributions of IRS-2 gene in the PCOM and non-PCOM groups are shown in Table 4. In the comparison of genotype distributions, the rs2289046 genotypes (AG and AA) and rs1865434AG genotype showed significantly lower frequencies in the PCOM compared to the non-PCOM group (58.9% vs. 71.9%,
Table 4
SNP | Genotype | Non-PCOM | PCOM |
|
| OR (95% CI) |
---|---|---|---|---|---|---|
rs2289046 | G/A | |||||
GG | 11 (8.6) | 28 (26.2) | 1 | |||
AG | 92 (71.9) | 63 (58.9) | 10.922 | 0.001 | 0.269 (0.125–0.580) | |
AA | 25 (19.5) | 16 (14.9) | 7.399 | 0.007 | 0.251 (0.098–0.642) | |
rs1865434 | A/G | |||||
AA | 72 (56.3) | 68 (63.5) | 1 | |||
AG | 56 (43.7) | 34 (31.8) | 8.391 | 0.004 | 0.643 (0.375–1.103) | |
GG | 5 (4.7) | 5.108 | 0.024 | 1.074 (1.009–1.142) |
OR: odds ratio; CI: confidence intervals;
3.3.2. All 235 subjects were divided into the HA group (T
Table 5
SNP | Genotype | Non-HA | HA |
|
| OR (95% CI) |
---|---|---|---|---|---|---|
rs2289046 | G/A | |||||
GG | 29 (18.6) | 12 (15.2) | 1 | |||
AG | 106 (67.9) | 49 (62) | 2.306 | 0.129 | 0.539 (0.264–1.102) | |
AA | 21 (13.5) | 18 (22.8) | 1.764 | 0.184 | 0.483 (0.192–1.213) | |
rs1865434 | A/G | |||||
AA | 90 (57.7) | 50 (63.3) | 1 | |||
AG | 63 (40.4) | 27 (34.2) | 0.803 | 0.067 | 0.771 (0.437–1.362) | |
GG | 3 (1.9) | 2 (2.5) | 0.039 | 0.738 | 1.200 (0.194–7.423) |
OR: odds ratio; CI: confidence intervals;
3.3.3. A total of 235 subjects were divided into IR (HOMA
Table 6
SNP | Genotype | Non-IR | IR |
|
| OR (95% CI) |
---|---|---|---|---|---|---|
rs2289046 | G/A | |||||
GG | 20 (13.9) | 19 (20.9) | 1 | |||
AG | 93 (64.6) | 62 (68.1) | 0.974 | 0.324 | 0.702 (0.347–1.421) | |
AA | 31 (21.5) | 10 (11) | 4.120 | 0.042 | 0.340 (0.131–0.878) | |
rs1865434 | A/G | |||||
AA | 85 (59) | 55 (60.4) | 1 | |||
AG | 58 (40.3) | 32 (35.2) | 0.324 | 0.246 | 0.853 (0.492–1.476) | |
GG | 1 (1) | 4 (4.4) | 3.316 | 0.063 | 6.182 (0.673–56.770) |
OR: odds ratio; CI: confidence intervals;
Table 7
SNP | Genotype | Non-obesity | Obesity |
|
| OR (95% CI) |
---|---|---|---|---|---|---|
rs2289046 | G/A | |||||
GG | 32 (17) | 7 (14.9) | 1 | |||
AG | 120 (63.8) | 35 (74.5) | 0.394 | 0.53 | 1.333 (0.542–3.281) | |
AA | 36 (19.2) | 5 (10.6) | 0.519 | 0.471 | 0.635 (0.183–2.200) | |
rs1865434 | A/G | |||||
AA | 110 (58.5) | 30 (63.8) | 1 | |||
AG | 76 (40.4) | 14 (29.8) | 1.221 | 0.024 | 0.675 (0.336–1.358) | |
GG | 2 (1.1) | 3 (6.4) | 4.086 | 0.044 | 5.500 (0.879–34.430) |
OR: odds ratio; CI: confidence intervals;
3.3.4. According to BMI, all 235 subjects were divided into the obesity (BMI
4.Discussion
PCOS is considered with highly heterogeneous clinical manifestations and unclear etiology syndrome of ovarian dysfunction [16]. Hyperandrogenism (HA) and polycystic ovary morphology (PCOM) are known as cardinal features [17], while other clinical manifestations like menstrual irregularities, insulin resistance (IR), and obesity, appear additionally. A meta-analysis showed that Asian and American women with PCOS presented a higher prevalence of T2DM with 4.4- and 4.7- fold increasing versus to those women without PCOS [7]. Insulin resistance occurs in PCOS patients with a variable range of 40–70% [18, 19]. Chen et al. showed that dysregulated glucose metabolism coexists among 23.7% of the Chinese patients with PCOS [20]. Li et al. indicated impaired
In this survey, we detected the two SNPs in IRS-2 gene SNP among women with PCOS and performed a parallel with the genotype frequencies with healthy women without PCOS. The rs2289046 genotype and the rs186534 genotype showed a significant association with PCOS and other laboratory measurements suggested that the IRS-2 gene polymorphism was correlated with the evolvement of PCOS and/or metabolic disorders; however, the genotypes were not associated with HA. The rs2289046 genotype (AA and AG) displayed a negative correlation with the development of PCOS and PCOM; meanwhile, the rs2289046 AA genotype was negatively associated with IR; and the rs1865434AG genotype was negatively associated with PCOM and obesity, make a consequence to a strong protective effect. However, the rs1865434 GG genotype showed a positive correlation to PCOS, PCOM and obesity and was a strong risk factor with a higher odds ratio. This series of findings indicated that the rs2289046 and rs186534 genotypes were associated with the development and clinical manifestations of PCOS, except for HA. Among the genotypes being analyzed, rs2289046 genotype (AA and AG) and rs1865434 genotype (AG and GG) showed significant associations with laboratory measurements, including FSH, LH, LH/FSH and HOMA-IR. Different genetic variants may contribute to the differences in serum T concentrations between the PCOS and controls.
Insulin resistance (IR) is a collaborative metabolic character of both obesity and PCOS. It was widely documented that women with PCOS, independent of obesity, are insulin resisted and have compensatory hyperinsulinemia [24]. Obesity sensitizes thecal cells to LH stimulation and amplifies functional ovarian hyperandrogenism by upregulating ovarian androgen production [25]. Furthermore, insulin might enhance androgenicity through the increased levels of free (biologically available) testosterone accordingly by suppressing production of sex hormone-binding globulin (SHBG) within the liver [26]. These factors lead to disturbances in follicle development and maturation, formatting small follicles and inducing follicular atresia, resulting in clinical manifestations of PCOS (menstrual irregularity or amenorrhea, polycystic ovarian enlargement, obesity, infertility).
Polycystic ovaries are the morphological ovarian phenotype in women with the PCOS. Diamanti et al. reported the prevalence rates of PCOM in the order of 60–80% and a positive relationship was found between the androgens levels, insulin resistance index, and PCOM [24]. Carmina et al. reported a positive correlation between insulin resistance indices and ovarian volume in 326 women with PCOS [27]. In multiple cell types, IRS-2 is involved in preserving insulin action whereas reduction in IRS expression and/or function may cause development of
In women genetically predisposed to the development of PCOS, weight-gain and obesity often make consequence to clinical and biochemical characters. Obesity is associated with PCOS in between 38% or 88% of cases [30]. A case-control association study showed that IRS-2 SNP (rs1865434) was significantly higher in 273 Caucasian with PCOS paralleled to 173 controls (
5.Conclusion
rs2289046 genotype (AA and AG) and rs1865434 AG genotype may have protective effects on PCOS development, while rs1865434 GG genotype may be a risk factor for the evolvement of PCOS. Our data-depending correlation between genotype, disease, biochemical and hormonal variables may support a role for IRS-2 gene polymorphisms in the pathogenesis of PCOS. Further investigations are still required to be performed in a larger number of subjects to clarify the roles of IRS-2 genetic variants in patients with PCOS.
Acknowledgments
The authors would like to thank all participants in this study. This work was supported by the Guangzhou Science and Technology Innovation Project (Grant number: 201804010003) and Guangdong Medical Research foundation (Grant number: A2020467).
Conflict of interest
None to report.
References
[1] | Ong M, Peng J, Jin XL, Qu XQ. Chinese herbal medicine for the optimal management of polycystic ovary syndrome. Am J Chin Med. (2017) ; 45: (3): 405-422. |
[2] | Wijeyaratne CN, Udayangani SAD, Balen AH. Ethnic-specific polycystic ovary syndrome: Epidemiology, significance and implications. Expert Rev Endocrinol Metab. (2013) ; 8: (1): 71-79. |
[3] | Burks HR, Wild RA. Diagnostic criteria and epidemiology of PCOS. In: Pal L (ed) Polycystic ovary syndrome: current and emerging concepts. Springer, New York, (2014) , pp. 3-10. |
[4] | Hutchison SK, Stepto NK, Harrison CL, Moran LJ, Strauss BJ, Teede HJ. Effects of exercise on insulin resistance and body composition in overweight and obese women with and without polycystic ovary syndrome. J Clin Endocrinol Metabol. (2011) ; 96: (1): 48-56. |
[5] | Kakoly NS, Earnest A, Teede HJ, Moran LJ, Joham AE. The impact of obesity on the incidence of type 2 diabetes among women with polycystic ovary syndrome. Diabetes Care. (2019) ; 42: (4): 560-567. |
[6] | Diamanti-Kandarakis E, Piperi C. Genetics of polycystic ovary syndrome: Searching for the way out of the labyrinth. Hum Reprod Update. (2005) ; 11: : 631-643. |
[7] | Kakoly NS, Khomami MB, Joham AE, Cooray SD, Misso ML, Norman RJ, Harrison CL, Ranasinha S, Teede HJ, Moran LJ. Ethnicity, obesity and the prevalence of impaired glucose tolerance and type 2 diabetes in PCOS: A systematic review and meta-regression. Human Reproduction Update. (2018) ; 24: : 455-467. |
[8] | El Mkadem SA, Lautier C, Macari F, et al. Role of allelic variants gly972arg of irs-1 and gly1057asp of irs-2 in moderate-to-severe insulin resistance of women with polycystic ovary syndrome. Diabetes. (2001) ; 50: (9): 2164-2168. |
[9] | Anagnostis P, Tarlatzis BC, Kauffman RP. Polycystic ovarian syndrome (PCOS): Long-term metabolic consequences. Metabolism. (2018) ; 86: : 33-43. |
[10] | Sørensen AE, Wissing ML, Salö S, et al. MicroRNAs related to polycystic ovary syndrome (PCOS). Genes (Basel). (2014) ; 5: (3): 684-708. |
[11] | Roth LW, Mccallie B, Alvero R, et al. Altered microRNA and gene expression in the follicular fluid of women with polycystic ovary syndrome. J Assist Reprod Genet. (2014) ; 31: (3): 355-362. |
[12] | Rotterdam E-A. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome (PCOS). Hum Reprod. (2004) ; 19: : 41-47. |
[13] | Toscani M, Migliavacca R, Sisson de Castro JA, Spritzer PM. Estimation of truncal adiposity using waist circumference or the sum of trunk skinfolds: A pilot study for insulin resistance screening in hirsute patients with or without polycystic ovary syndrome. Metabolism. (2007) ; 56: (7): 992-997. |
[14] | Wiltgen D, Spritzer PM. Variation in metabolic and cardiovascular risk in women with different polycystic ovary syndrome phenotypes. Fertil Steril. (2010) ; 94: (6): 2493-2496. |
[15] | Xing X, Yang ZJ. The diagnostic significance of homeostasis model assessment of insulin resistance in metabolic syndrome among subjects with different glucose tolerance. Chin J Diab. (2004) ; 1: : 182-186. |
[16] | Long W, Zhao C, Ji C, Ding H, Cui Y, Guo X, et al. Characterization of serum microRNAs profile of PCOS and identification of novel non-invasive biomarkers. Cell Physiol Biochem. (2014) ; 33: (5): 1304-1315. |
[17] | Laven JS, Imani B, Eijkemans MJ, Fauser BC. New approach to polycystic ovary syndrome and other forms of anovulatory infertility. Obstet Gynecol Surv. (2002) ; 57: (11): 755-767. |
[18] | Nawrocka-Rutkowska J, Ciecwiez S, Marciniak A, et al. Insulin resistance assessment in patients with polycystic ovary syndrome using different diagnostic criteria-impact of metformin treatment. Ann Agric Environ Med. (2013) ; 20: : 528-532. |
[19] | Tabassum R, Imtiaz F, Sharafat S, Shukar-Ud-Din S, Nusrat U. Prevalence and clinical profile of insulin resistance in young women of poly cystic ovary syndrome: A study from Pakistan. Pak J Med Sci. (2013) ; 29: (2): 593-596. |
[20] | Chen C, Jing G, Li Z, Juan S, Bin C, Jie H. Insulin resistance and polycystic ovary syndrome in a Chinese population. 278 Endocr Pract. (2017) . doi: 10.4158/EP171849.OR. |
[21] | Li W, Li Q. Dysregulation of glucose metabolism even in Chinese PCOS women with normal glucose tolerance. Endocr J. (2012) ; 59: : 765-770. |
[22] | Shi X, Xie X, Jia Y, Li S. Associations of insulin receptor and insulin receptor substrates genetic polymorphisms with polycystic ovary syndrome: A systematic review and meta-analysis. J. Obstet. Gynaecol. Res. (2016) ; 42: (7): 844-854. |
[23] | Feigelson HS, Teras LR, Diver WR, Tang W, et al. Genetic variation in candidate obesity genes ADRB2, ADRB3, GHRL, HSD11B1, IRS1, IRS2, and SHC1 and risk for breast cancer in the Cancer Prevention Study II. Breast Cancer Res. (2008) ; 10: (4): R57. |
[24] | Diamanti-Kandarakis E, Panidis D. Unravelling the phenotypic map of polycystic ovary syndrome (PCOS): A prospective study of 634 women with PCOS. Clin Endocrinol (Oxf). (2007) ; 67: (5): 735-742. |
[25] | Charles JG, Naila G. Characteristics of obesity in polycystic ovary syndrome: Etiology, treatment, and genetics. Metabolism. (2019) ; 92: : 108-120. |
[26] | Lim SS, Norman RJ, Davies MJ, Moran LJ. The effect of obesity on polycystic ovary syndrome: a systematic review and meta-analysis. Obes Rev. (2013) ; 14: (2): 95-109. |
[27] | Carmina E, Orio F, Palomba S, Longo RA, Lombardi G, Lobo RA. Ovarian size and blood flow in women with polycystic ovary syndrome and their correlations with endocrine parameters. Fertil Steril. (2005) ; 84: (2): 413-419. |
[28] | Brady MJ. IRS2 takes center stage in the development of type 2 diabetes. J Clin Invest. (2004) ; 114: (7): 886-888. |
[29] | Perez-Martinez P, Delgado-Lista J, Garcia-Rios A, et al. Insulin receptor substrate-2 gene variants in subjects with metabolic syndrome: Association with plasma monounsaturated and n-3 polyunsaturated fatty acid levels and insulin resistance. Mol Nutr Food Res. (2012) ; 56: (2): 309-315. |
[30] | Barber TM, Franks S. Genetics of polycystic ovary syndrome. Polycystic Ovary Syndrome. 40: Karger Publishers; (2013) . 28-39. |
[31] | Goodarzi MO, Louwers YV, Taylor KD, et al. Replication of association of a novel insulin receptor gene polymorphism with polycystic ovary syndromeFertil. Steril. (2011) ; 95: (5): 1736-1741. |
[32] | Butte NF, Voruganti VS, Cole SA, Haack K, et al. Resequencing of IRS2 reveals rare variants for obesity but not fasting glucose homeostasis in Hispanic children. Physiol Genomics. (2011) ; 43: (18): 1029-1037. |