# Estrogen receptor 1 gene polymorphisms (PvuII and XbaI) are associated with type 2 diabetes in Palestinian women

Suheir Ereqat, Stephane Cauchi, Khaled Eweidat, Muawiyah Elqadi,
Abedelmajeed Nasereddin

► To cite this version:

Suheir Ereqat, Stephane Cauchi, Khaled Eweidat, Muawiyah Elqadi, Abedelmajeed Nasereddin. Estrogen receptor 1 gene polymorphisms (PvuII and XbaI) are associated with type 2 diabetes in Palestinian women. PeerJ, 2019, 7, pp.e7164. 10.7717/peerj.7164. hal-02346725

HAL Id: hal-02346725

https://hal.science/hal-02346725

Submitted on 18 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Estrogen receptor 1 gene polymorphisms (PvuII and XbaI) are associated with type 2 diabetes in Palestinian women

Suheir Ereqat<fnref n="1" />, Stéphane Cauchi<fnref n="2" />,<fnref n="3" />,<fnref n="4" />,<fnref n="5" />, Khaled Eweidat<fnref n="1" />, Muawiyah Elqadi<fnref n="1" /> and Abedelmajeed Nasereddin<fnref n="6" />

## ABSTRACT

**Background:** Type 2 diabetes mellitus (T2DM) is a multifactorial disease where both genetic and environmental factors contribute to its pathogenesis. The PvuII and XbaI polymorphisms of the estrogen receptor 1 (ESR1) gene have been variably associated with T2DM in several populations. This association has not been studied in the Palestinian population. Therefore, the aim of this study was to investigate the association between the PvuII and XbaI variants in the *ESR1* and T2DM and its related metabolic traits among Palestinian women.

**Methods:** This case-control study included 102 T2DM and 112 controls in which PvuII and XbaI variants of the *ESR1* gene were genotyped using amplicon based next generation sequencing (NGS).

**Results:** Allele frequencies of both PvuII and XbaI variants were not significantly different between patients and control subjects (*P* > 0.05). In logistic regression analysis adjusted for age and BMI, the *ESR1* PvuII variant was associated with risk of T2DM in three genotypic models (*P* < 0.025) but the strongest association was observed under over-dominant model (TT+CC vs. TC) (OR = 2.32, CI [1.18–4.55] adjusted *P* = 0.013). A similar but non-significant trend was also observed for the *ESR1* XbaI variant under the over-dominant model (AA+GG vs. AG) (OR = 2.03, CI [1.05–3.95]; adjusted *P* = 0.035). The frequencies of the four haplotypes (TA, CG, CA, TG) were not significantly different in the T2DM patients compared with control group (*P* > 0.025). Among diabetic group, an inverse trend with risk of cardio vascular diseases was shown in carriers of CG haplotype compared to those with TA haplotype (OR = 0.28, CI [0.09–0.90]; adjusted *P* = 0.035). Further, stratified analyses based on *ESR1* PvuII and XbaI genotypes revealed no evidence for association with lipid levels (TC, TG, HDL, LDL).

**Conclusions:** This is the first Palestinian study to conclude that *ESR1* PvuII and XbaI variants may contribute to diabetes susceptibility in Palestinian women. Identification of genetic risk markers can be used in defining high risk subjects and in prevention trials.

Submitted 13 February 2019

Accepted 21 May 2019

Published 28 June 2019

Corresponding author

Suheir Ereqat,
sereqat@staff.alquds.edu

Academic editor

Juan Riesgo-Escovar

Additional Information and
Declarations can be found on
page 11

DOI 10.7717/peerj.7164

cc Copyright

2019 Ereqat et al.

Distributed under

Creative Commons CC-BY 4.0

OPEN ACCESS

<footnote n="1">¹ Biochemistry and Molecular Biology Department-Faculty of Medicine, Al-Quds University, East Jerusalem, Palestine</footnote>

<footnote n="2">² CNRS, UMR8204, Lille, France</footnote>

<footnote n="3">³ INSERM, U1019, Lille, France</footnote>

<footnote n="4">⁴ Université de Lille, Lille, France</footnote>

<footnote n="5">⁵ Institut Pasteur de Lille, Centre d'Infection et d'Immunité de Lille, Lille, France</footnote>

<footnote n="6">⁶ Al-Quds Public Health Society, East Jerusalem, Palestine</footnote>
**Subjects** Bioinformatics, Molecular Biology, Diabetes and Endocrinology, Medical Genetics

**Keywords** T2DM, ESR1, PvuII and XbaI polymorphisms, NGS, CVD, Palestine

## INTRODUCTION

Type 2 diabetes mellitus (T2DM) is a multifactorial disease that caused by a complex combination of genetic and environmental factors. Identification of genetic polymorphisms associated with diabetes may lead to prediction of disease development and prevention of its vascular complications (Slominski et al., 2018). Recently, several reports have revealed the role of estradiol in regulating energy metabolism (Hevener et al., 2018). It was evident that estrogens increase hepatic insulin sensitivity, stimulate insulin synthesis in islets of Langerhans, prevent β-cell apoptosis and improve insulin action in skeletal muscles (Meyer et al., 2011).

Estrogen exerts its physiological functions through the estrogen receptors (ESR1, ESR2, and a G-protein coupled cell surface receptor) and might prevent menopause syndrome, cardio vascular diseases (CVD) and diabetes (Gupte, Pownall & Hamilton, 2015). Estrogen receptor 1 (ESR1) is broadly expressed in adipose tissue, skeletal muscle, liver, and immune cells. It is a ligand-activated transcription factor that regulates a large number of genes in diverse target tissues (Hevener et al., 2018). Animal studies reported that male and female ESR1 gene knockout mice developed features of the metabolic syndrome (MetS) including obesity caused by impaired fatty acid oxidation, glucose intolerance, and impaired insulin sensitivity, thus revealed the critical role of ESR1 in metabolic homeostasis (Heine et al., 2000; Ribas et al., 2010).

The ESR1 gene, which encompasses 140 kb of DNA, is found on chromosomes 6q25.1. It is a highly polymorphic gene containing more than 1,600 single nucleotide polymorphisms (SNPs). Two SNPs in ESR1, both located in the first intron, PvuII (rs2234693) and XbaI (rs9340799), are the most extensively investigated variants and reported to be associated with T2DM (Huang et al., 2006; Mohammadi et al., 2013), MetS (Zhao et al., 2018), and other diseases (Onland-Moret et al., 2005; Silva et al., 2010; Weng et al., 2015). These polymorphisms might interfere with the estrogen effect by altering ESR1 gene expression via altering the binding of its own transcription factors (Gomes-Rochette et al., 2017). Diabetes and obesity had reached an alarming rate among Palestinians especially among refugees (El Kishawi et al., 2014; Shahin, Kapur & Seita, 2015). In 2015, the Palestinian Ministry of Health reported that cardiovascular diseases were the leading cause of death and diabetes mellitus was the fourth cause of death among Palestinians (Palestinian Ministry of Health, 2015). A recent study conducted on three refugee camps in the West Bank, revealed that the overall prevalence of obesity and overweight was 63.1% (for obesity: 42% in women and 29.2% in men; overweight 25.8% and 28.9% in women and men; respectively). The prevalence of MetS among obese and overweight adults was 69.4% with no gender-based differences (Damiri et al., 2018). Obese and overweight individuals are at high risk of both cardiovascular diseases and T2DM (Alberti & Zimmet, 1998). Estrogen deficiency is an important obesity-triggering factor which enhances metabolic dysfunction and thus predisposing to T2DM and cardiovascular diseases among menopausal women (Lizcano & Guzman, 2014). A meta-analysis study
investigating the new-onset of T2DM in post-menopausal women following estrogen replacement therapy revealed a 30% lower relative risk [RR 0.7 (CI [0.6–0.9])] of diabetes compared with placebo (Salpeter et al., 2006). However, the use of estrogen to prevent chronic diseases is still challenging and controversial (Bolton, 2016).

As the *ESR1* gene is a potential candidate gene for susceptibility to T2DM, we hypothesized that the *ESR1* variants (PvuII and XbaI) might be associated with T2DM, diabetic complications and related metabolic traits in Palestinian diabetic women.

# MATERIALS AND METHODS

## Study design and participants

This case control study includes 214 Palestinian women (102 with T2DM and 112 non-diabetic controls). All participants were unrelated, aged >40 years and selected within the period of 2016–2017 from the United Nations Relief and Works Agency for Palestine Refugees clinics (Hebron and Ramallah, Palestine). Diagnosis of T2DM was based on World Health Organization criteria (fasting plasma glucose ≥ 126 mg/dl and/or currently on treatment for diabetes) (Kumar et al., 2016). The diagnosis of T2DM was confirmed based on patients’ medical records reporting initial diagnosis of diabetes mellitus after age of 40, receiving oral hypoglycemic agents without insulin for at least 1 year after initial diagnosis, and currently on treatment for diabetes. Patients with probable type 1 diabetes who received continuous insulin therapy since diagnosis were excluded. The anthropometric measurements were collected from their medical records using a standard questionnaire that included age, age at diagnosis, sex, family history, diabetic complications, and medication for diabetes. BMI was calculated as kilograms divided by the square of height in meters. Blood pressure was measured in sitting position, on the left arm, after a 5-min rest by a health clinic worker, with a mercury sphygmomanometer. Blood samples (five ml) after a 12-h minimum fast was collected in EDTA tubes for biochemical tests as described in our previous work (Sabarneh et al., 2018). Plasma glucose, cholesterol, HDL cholesterol, and triglyceride were determined using standard methods of commercial kits (Human, Wiesbaden, Germany). Low-density lipoprotein (LDL) cholesterol was calculated based on the Friedewald formula. The control group was selected from individuals who came to the same clinic for an annual health check-up, they were eligible to be included if they had no prior diagnosis for T2DM, no family history in first-degree relatives, their fasting glucose < 126 mg/dl and BMI > 25. All the participants provided written informed consent to participate in the study, the study protocol was approved by Al-Quds University Research Ethics Committee (Rf no. 2/SRC/4).

## Amplicon based next generation genotyping and bioinformatics

Genomic DNA was extracted from whole blood (300 µl) using genomic QIAamp DNA purification kit according to the manufacturer instructions (Qiagen, Hilden, Germany). All DNA samples were genotyped for the T/C, rs2234693, and A/G, rs9340799, also known as PvuII and XbaI polymorphisms, respectively using amplicon based next generation sequencing (NGS). Briefly, two primers (forward and reverse) were used to target the two SNPs as previously described (Motawi et al., 2015). Both primers were
modified with over hanged Illumina adaptor sequences at the 5' ends (italic bolded, *Table S1*) to target a partial sequence of 119 bp in length with a final product of 186 bp using conventional thermocycler.

The PCR product was visualized and captured on a 1.5% agarose gel, cleaned by Agencourt AMPure XP system (X1, A63881; Beckman Coulter Genomics, Indianapolis, IN, USA) and eluted in 25 µl elution buffer. All purified products were subjected to a second round of amplification to assign unique index sequences (barcode) for each sample using Nextera XT Index Kit (Illumina, San Diego, CA, USA). Five µl from each barcoded sample were pooled together, mixed and spin down. Then, 100 µl of the pooled product was cleaned by Agencourt AMPure XP system (X1) (A63881; Beckman Coulter Genomics), and eluted in 50 µl elution buffer. Library purity and quantity were evaluated by 4200 TapeStation System (Agilent Technologies, Inc., Santa Clara, CA, USA) using D1000 ScreenTape kit (Agilent Technologies, Inc., Santa Clara, CA, USA) and by Qubit® Fluorometer (Invitrogen, Carlsbad, CA, USA) using Qubit dsDNA high-sensitivity assay (Invitrogen, Carlsbad, CA, USA). Concentration of four nM was prepared. 20 K reads for each sample was targeted. Samples were deep sequenced on NextSeq 500/550 machine using the 150-cycle Mid Output Kit (Illumina, San Diego, CA, USA).

The obtained DNA sequences were uploaded on the Galaxy program (https://usegalaxy.org/). Workflow of filtration included Illumina adaptor trim, quality selection of Q > 20 with minimal read length of 100 bp. Four virtual probe sequences were used to identify the PvuII and XbaI variants (*Table S1*). The genotypes were determined based on the calculated ratio between the read counts for wild type and mutant alleles, for both SNPs (PvuII and XbaI) in each individual sample.

## Statistical analysis

Statistical analysis was performed with the SPSS package, version 19.0 (SPSS, Inc., Chicago, IL, USA). All tests were two-tailed and *P* < 0.05 was considered significant. The mean values and standard errors were reported in tables. Genotype and allele frequencies in T2DM and control subjects were tested by multivariable logistic regression analysis using five genetic models: co-dominant, dominant, over-dominant, recessive and additive with adjustment for age and BMI using R statistics (V 3.5.1; SNPassoc package) (Gonzalez et al., 2007). The genotype frequencies were tested for Hardy-Weinberg equilibrium by calculating a χ² statistic and corresponding *P*-value as previously described by Mayo (2008). The risk to T2DM was estimated by computing odds ratios (ORs) and 95% confidence intervals (CIs). Linkage disequilibrium (LD), coefficient (D') for haplotypes and their frequencies were performed using R statistics (V 3.5.1; SNPassoc package) (Gonzalez et al., 2007). Multivariate linear regression models taking into account age and BMI were performed to assess associations between SNPs and T2DM complications.

All quantitative parameters were normalized before analysis using the “bestNormalize” R package (https://github.com/petersonR/bestNormalize) (Peterson, 2017). To minimize the chance of obtaining type I error, Bonferroni's adjustment was applied. *P* < 0.025 was adopted as the significant threshold unless otherwise specified. Quanto 1.2.4 software
<table><caption>**Table 1** Demographic and biochemical characteristics of study subjects.</caption><thead><tr><th></th><th>Control (n = 112)</th><th>T2DM (n = 102)</th><th>P-value</th></tr></thead><tbody><tr><td>Age at sampling (years)</td><td>48.9 (0.68)</td><td>59.9 (0.95)</td><td>0.0001</td></tr><tr><td>Age at diagnosis (years)</td><td>NA</td><td>49.6 (0.83)</td><td></td></tr><tr><td>BMI (Kg/m<sup>2</sup>)</td><td>32.5 (0.52)</td><td>34 (0.6)</td><td>0.054</td></tr><tr><td>SBP (mmHg)</td><td>122 (0.93)</td><td>136.2 (1.6)</td><td>0.0001</td></tr><tr><td>DBP (mmHg)</td><td>77.9 (0.85)</td><td>78.9 (1.02)</td><td>0.4</td></tr><tr><td>FBS (mg/dl)</td><td>90.1 (0.97)</td><td>161.8 (4.9)</td><td>0.0001</td></tr><tr><td>TC (mg/dl)</td><td>182 (3.4)</td><td>192 (4.02)</td><td>0.03</td></tr></tbody></table>

**Notes:**

BMI, body mass index; SBP, systolic blood pressure; DBP, diastolic blood pressure; FBS, fasting blood sugar.

Data are presented as mean (standard error); *P* < 0.05 was considered significant (in bold). NA, Not applicable.

was used to estimate the statistical power (http://biostats.usc.edu/Quanto.html). The power of any test of statistical significance was <50%.

## RESULTS

### Biochemical and genetic analysis

Anthropometric and biochemical characteristics of study individuals are presented in Table 1. In comparison with the control group, all clinical parameters showed statistically significant increase in T2DM patients (*P* < 0.05). However, no significant difference in diastolic blood pressure was observed among the two group (*P* = 0.4). The T2DM subjects were significantly older (*P* = 0.0001) but the control group were at the same mean age at T2DM diagnosis in cases (*P* > 0.05) (Table 1). The cases had higher mean BMI than control subjects but was at borderline significance (*P* = 0.054). Among T2DM patients, 10.8, 8.8, 7.8, and 5.9% had CVD, nephropathy, diabetic foot, and retinopathy, respectively. The genotyping distribution for both variants was consistent with Hardy Weinberg equilibrium in cases and control groups and in all subjects (*P* > 0.05).

Our results showed that the heterozygous genotypes TC and AG of the *ESR1* PvuII and XbaI polymorphisms were significantly elevated in T2DM patients compared with control individuals (*P* = 0.04, *P* = 0.01, respectively), while, the frequencies of C and G alleles of the *ESR1* PvuII and XbaI polymorphisms, respectively, were comparable among cases and controls (*P* > 0.05) (Table 2). However, upon Bonferroni's adjustment, only the heterozygous genotype AG of the *ESR1* XbaI remained significantly higher in diabetic group (*P* = 0.01).

### PvuII and XbaI Polymorphisms in *ESR1* gene and risk to type 2 diabetes

To estimate the effect of the genotypes on the disease, logistic regression analysis was performed using five genetic models adjusted to age and BMI: additive, dominant, codominant, recessive, and over-dominant. In the case of the PvuII T/C rs2234693 polymorphism, T2DM patients showed significantly different genotypes distribution compared to control group in three models, carriers of TC genotype had higher risk to T2DM compared to those of TT genotype (OR = 2.84, CI [1.32–6.11]; *P* = 0.024)
<table><caption>**Table 2** Allele frequencies and genotypes distribution of *ESR1* gene PvuII and XbaI polymorphisms in control and T2DM patients.</caption><thead><tr><th>Allele/genotype</th><th>All subjects (%)</th><th>Control n (%)</th><th>T2DM n (%)</th><th>P-value</th></tr></thead><tbody><tr><td>PvuII</td><td></td><td></td><td></td><td></td></tr><tr><td>T</td><td>244 (57)</td><td>133 (59)</td><td>111 (54)</td><td>0.46</td></tr><tr><td>C</td><td>184 (43)</td><td>91 (41)</td><td>93 (46)</td><td>0.46</td></tr><tr><td>TT</td><td>70 (33)</td><td>43 (38)</td><td>27 (26)</td><td>0.06</td></tr><tr><td>TC</td><td>104 (49)</td><td>47 (42)</td><td>57 (56)</td><td>0.04</td></tr><tr><td>CC</td><td>40 (19)</td><td>22 (20)</td><td>18 (18)</td><td>0.7</td></tr><tr><td>XbaI</td><td></td><td></td><td></td><td></td></tr><tr><td>A</td><td>247 (58)</td><td>132 (59)</td><td>115 (56)</td><td>0.6</td></tr><tr><td>G</td><td>118 (42)</td><td>92 (41)</td><td>89 (44)</td><td>0.7</td></tr><tr><td>AA</td><td>73 (34)</td><td>44 (39)</td><td>29 (28)</td><td>0.09</td></tr><tr><td>AG</td><td>101 (47)</td><td>44 (39)</td><td>57 (56)</td><td><strong>0.01</strong></td></tr><tr><td>GG</td><td>40 (19)</td><td>24 (21)</td><td>16 (16)</td><td>0.35</td></tr></tbody></table>

**Notes:**

Bold numbers showed the significant correlation.

P-values lower than the Bonferroni threshold (P = 0.025) were considered statistically significant.

<table><caption>**Table 3** Association of *ESR1* PvuI and XbaI variants with T2DM.</caption><thead><tr><th>Model</th><th>Genotype</th><th>OR (95% CI)</th><th>*P-value</th></tr></thead><tbody><tr><td>PvuII SNP</td><td></td><td></td><td></td></tr><tr><td>Codominant</td><td>T/T</td><td>1</td><td></td></tr><tr><td></td><td>T/C</td><td>2.84 [1.32–6.11]</td><td><strong>0.024</strong></td></tr><tr><td></td><td>C/C</td><td>1.77 [0.66–4.73]</td><td></td></tr><tr><td>Dominant</td><td>T/T</td><td>1</td><td></td></tr><tr><td></td><td>T/C–C/C</td><td>2.50 [1.21–5.14]</td><td><strong>0.011</strong></td></tr><tr><td>Recessive</td><td>T/T–T/C</td><td>1</td><td></td></tr><tr><td></td><td>C/C</td><td>0.96 [0.41–2.25]</td><td>0.92</td></tr><tr><td>Overdominant</td><td>T/T–C/C</td><td>1</td><td></td></tr><tr><td></td><td>T/C</td><td>2.32 [1.18–4.55]</td><td><strong>0.013</strong></td></tr><tr><td>Log-additive</td><td>–</td><td>1.48 [0.92–2.37]</td><td>0.1</td></tr><tr><td>Xba SNP</td><td></td><td></td><td></td></tr><tr><td>Codominant</td><td>A/A</td><td>1</td><td></td></tr><tr><td></td><td>A/G</td><td>2.15 [1.02–4.51]</td><td>0.1</td></tr><tr><td></td><td>G/G</td><td>1.18 [0.45–3.10]</td><td></td></tr><tr><td>Dominant</td><td>A/A</td><td>1</td><td></td></tr><tr><td></td><td>A/G–G/G</td><td>1.82 [0.91–3.66]</td><td>0.089</td></tr><tr><td>Recessive</td><td>A/A–A/G</td><td>1</td><td></td></tr><tr><td></td><td>G/G</td><td>0.76 [0.32–1.80]</td><td>0.53</td></tr><tr><td>Overdominant</td><td>A/A–G/G</td><td>1</td><td></td></tr><tr><td></td><td>A/G</td><td>2.03 [1.05–3.95]</td><td>0.035</td></tr><tr><td>Log-additive</td><td>–</td><td>1.21 [0.76–1.92]</td><td>0.42</td></tr></tbody></table>

**Note:**

*P-values were from logistic regression models adjusted for age and BMI, P < 0.025 was considered significant, bold numbers showed the significant correlation.
<table><caption>**Table 4** Demographic characteristics and biochemical measurements based on ESR1 PvuII and XbaI genotypes in diabetic women.</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="4">PvuII genotypes</th><th colspan="4">XbaI genotypes</th></tr><tr><th>CC</th><th>TC</th><th>TT</th><th>P-value</th><th>AA</th><th>AG</th><th>GG</th><th>P-value</th></tr></thead><tbody><tr><td>BMI (Kg/m²)</td><td>34.5 (1.8)</td><td>33.7 (0.69)</td><td>34.1 (1.3)</td><td>0.91</td><td>34.5 (1.24)</td><td>33.5 (0.69)</td><td>34.8 (1.9)</td><td>0.67</td></tr><tr><td>SBP (mmHg)</td><td>143.5 (3.50)</td><td>135.5 (2.2)</td><td>133 (2.9)</td><td>0.09</td><td>134.6 (2.75)</td><td>135.3 (2.27)</td><td>142.4 (3.8)</td><td>0.25</td></tr><tr><td>DBP (mmHg)</td><td>79 (3.0)</td><td>78.8 (1.4)</td><td>78.7 (1.6)</td><td>0.93</td><td>79.9 (1.65)</td><td>78.5 (1.4)</td><td>78.1 (3.1)</td><td>0.79</td></tr><tr><td>FBS (mg/dl)</td><td>169.85 (13.2)</td><td>162.82 (6.84)</td><td>154.44 (8.1)</td><td>0.59</td><td>159.24 (9.39)</td><td>159.14 (6.36)</td><td>176.21 (13.98)</td><td>0.46</td></tr><tr><td>HbA1C</td><td>8.28 (0.47)</td><td>7.86 (0.26)</td><td>7.72 (0.38)</td><td>0.62</td><td>7.83 (0.4)</td><td>7.78 (0.24)</td><td>8.46 (0.52)</td><td>0.46</td></tr><tr><td>TC (mg/dl)</td><td>202.42 (11.68)</td><td>188.29 (5.3)</td><td>192.33 (7.39)</td><td>0.45</td><td>191.52 (6.84)</td><td>189.57 (5.59)</td><td>200.6 (11.78)</td><td>0.65</td></tr><tr><td>TG (mg/dl)</td><td>250.07 (62.75)</td><td>192.87 (19.25)</td><td>164.45 (6.84)</td><td>0.2</td><td>164.9 (6.37)</td><td>212.01 (27.22)</td><td>191.75 (15.81)</td><td>0.42</td></tr><tr><td>LDL (mg/dl)</td><td>118.3 (7.69)</td><td>105.95 (3.29)</td><td>105.59 (4.78)</td><td>0.19</td><td>106.1 (4.47)</td><td>106.57 (3.53)</td><td>116.78 (7.48)</td><td>0.36</td></tr><tr><td>HDL (mg/dl)</td><td>38.56 (2.8)</td><td>41.58 (1.19)</td><td>44.55 (1.62)</td><td>0.11</td><td>44.38 (1.59)</td><td>41.14 (1.2)</td><td>39.69 (3.04)</td><td>0.2</td></tr></tbody></table>

**Notes:**

TC, total cholesterol; TG, triglyceride; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol.

Values are presented as mean (standard error); *P*-value was obtained by ANOVA, *P* < 0.05 was considered significant.

<table><caption>**Table 5** Estimated haplotype distributions in cases (T2DM) and control groups.</caption><thead><tr><th>PvuII</th><th>XbaI</th><th>Frequency (case)</th><th>Frequency (control)</th><th>Frequency (total)</th><th>OR (95% CI)</th><th>P-value</th><th>Global P</th></tr></thead><tbody><tr><td>T</td><td>A</td><td>0.5391</td><td>0.557</td><td>0.5485</td><td>1</td><td>-</td><td>0.14</td></tr><tr><td>C</td><td>G</td><td>0.4312</td><td>0.374</td><td>0.4013</td><td>1.36 [0.84-2.21]</td><td>0.21</td><td></td></tr><tr><td>C</td><td>A</td><td>0.0247</td><td>0.0323</td><td>0.0286</td><td>2.28 [0.62-8.37]</td><td>0.22</td><td></td></tr><tr><td>T</td><td>G</td><td>0.005</td><td>0.0367</td><td>0.0216</td><td>0.21 [0.02-2.66]</td><td>0.23</td><td></td></tr></tbody></table>

(Table 3). In the dominant model, the TC+CC genotype carriers was significantly at higher risk to T2DM than TT genotypes (OR = 2.50, CI [1.21–5.14]; *P* = 0.011). In the over-dominant model (TT+CC vs. TC), the association was also significant (OR = 2.32, CI [1.18–4.55]; *P* = 0.013) (Table 3). The analysis of the Xba rs9340799 A/G polymorphism showed a trend for association with T2DM in the over dominant model (AA+GG vs. AG) (OR = 2.03, CI [1.05–3.95]; *P* = 0.035) (Table 3).

Biochemical characteristics of T2DM individuals according to ESR1 PvuII and XbaI genotypes are shown in Table 4. It demonstrated insignificant association of the investigated parameters across genotypes of both polymorphisms (*P* > 0.05). A stratified analysis of combined genotypes with respect to different genetic models also revealed no significant association with any biochemical or clinical parameter (data not shown).

## LD estimation between ESR1 SNPs and haplotype analysis

The polymorphisms tested were in LD (normalized Lewontin's *D'* = 0.91). Among all subjects, the most common haplotype, TA, had a frequency of 55%, and the CG haplotype had a frequency of 40% while two haplotypes CA and TG had frequencies of 2.9% and 2.2%, respectively. The frequencies of the four haplotypes were not significantly different in the T2DM patients compared with control group as shown in Table 5. In logistic
<table><caption>**Table 6** Association of *ESR1* PvuII and XbaI variants with CVD among T2DM cases.</caption><thead><tr><th colspan="2">PvuII SNP</th><th colspan="2">CVD</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">P-value</th></tr><tr><th>Model</th><th>Genotype</th><th>Yes (%)</th><th>No (%)</th></tr></thead><tbody><tr><td rowspan="2">Dominant</td><td>T/T</td><td>54.5</td><td>23.1</td><td>1</td><td></td></tr><tr><td>T/C-C/C</td><td>45.5</td><td>76.9</td><td>0.25 [0.07-0.93]</td><td>0.039</td></tr><tr><td rowspan="2">Recessive</td><td>T/T-T/C</td><td>90.9</td><td>81.3</td><td>1</td><td></td></tr><tr><td>C/C</td><td>9.1</td><td>18.7</td><td>0.41 [0.05-3.57]</td><td>0.36</td></tr><tr><td rowspan="2">Overdominant</td><td>T/T-C/C</td><td>63.6</td><td>41.8</td><td>1</td><td></td></tr><tr><td>T/C</td><td>36.4</td><td>58.2</td><td>0.42 [0.11-1.57]</td><td>0.19</td></tr></tbody></table><table><thead><tr><th colspan="2">Xba SNP</th><th colspan="2">CVD</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">P-value</th></tr><tr><th>Model</th><th>Genotype</th><th>Yes (%)</th><th>No (%)</th></tr></thead><tbody><tr><td rowspan="2">Dominant</td><td>A/A</td><td>54.5</td><td>25.3</td><td>1</td><td></td></tr><tr><td>A/G-G/G</td><td>45.5</td><td>74.7</td><td>0.26 [0.07-0.99]</td><td>0.049</td></tr><tr><td rowspan="2">Recessive</td><td>A/A-A/G</td><td>100</td><td>82.4</td><td>1</td><td></td></tr><tr><td>G/G</td><td>0</td><td>17.6</td><td>0.00 [0.00-NA]</td><td>0.036</td></tr><tr><td rowspan="2">Overdominant</td><td>A/A-G/G</td><td>54.5</td><td>42.9</td><td>1</td><td></td></tr><tr><td>A/G</td><td>45.5</td><td>57.1</td><td>0.64 [0.18-2.31]</td><td>0.49</td></tr></tbody></table><table><thead><tr><th colspan="2">PvuII/XbaI</th><th colspan="2">CVD</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">P-value</th></tr><tr><th>Model</th><th>Haplotype</th><th>Yes (%)</th><th>No (%)</th></tr></thead><tbody><tr><td rowspan="3">T/A</td><td>73</td><td>52</td><td>1</td><td></td><td></td></tr><tr><td>C/G</td><td>23</td><td>46</td><td>0.28 [0.09-0.90]</td><td>0.035</td></tr><tr><td>C/A</td><td>4</td><td>2</td><td>2.75 [0.22-33.96]</td><td>0.43</td></tr></tbody></table>

**Note:**

NA, not applicable.

regression analyses, adjusting for the same potential confounders (age and BMI) used in the genotype models, none of the four possible haplotypes were associated with increased risk of diabetes (global *P* = 0.14) (Table 5).

## Association of *ESR1* PvuII and XbaI Polymorphisms with CVD in diabetic group

The association of PvuII and XbaI SNPs within the *ESR1* gene and CVD risk was evaluated among diabetic group (n = 102) using logistic regression after adjustment of age and BMI. Under the dominant model, the combined PvuII TC+CC genotypes have a trend to lower risk of CVD compared with TT genotype (OR = 0.25, CI [0.07-0.93]; *P* = 0.039). A similar trend was observed in the combined genotypes (AG+GG) of the *ESR1* XbaI compared to AA (OR = 0.26, CI [0.07-0.99]; *P* = 0.049) as shown in Table 6. The frequency of CG haplotype was statistically higher in diabetic patients without CVD (46%) compared to those with CVD (23%) (OR = 0.28, CI [0.09-0.90]; adjusted *P* = 0.035). On the other hand, we did not find any association of *ESR1* PvuII and XbaI polymorphisms with the other T2DM complications including nephropathy, retinopathy, and diabetic foot (data not shown).
DISCUSSION

The prevalence of diabetes has increased worldwide, most markedly in the world's middle-income countries which reflects an increase in associated risk factors such as being overweight or obese (World Health Organization (WHO), 2016). At menopause, women experience health challenges due to estrogen deficiency. The risk of T2DM is increasing with post-menopause status interacting with other risk factors, that is, hypertension, dyslipidemia, and obesity (Ren et al., 2018). Estrogen may regulate insulin action directly via actions on insulin-sensitive tissues. In skeletal muscle, ESR1 is thought to have a positive effect on insulin signaling and GLUT4 expression. It has been shown that stimulation of estrogen receptor with its agonist propylpyrazoletriyl increased insulin-stimulated glucose uptake in skeletal muscles (Gorres et al., 2011). Although not all studies are in agreement, PvuII and XbaI and other polymorphisms across the ESR1 gene have been associated with risk of T2DM as reported in the Chinese (Huang et al., 2006), African-Americans, European-Americans (Sale et al., 2004), Hungarians (Speer et al., 2001), and Egyptian women (Motawi et al., 2015). In this study, the frequency for ESR1 PvuII and XbaI alleles and genotypes in Palestinian women (n = 214) was evaluated. Our results reported that the ESR1 PvuII C and XbaI G allele frequencies were similar to other studies among Arab population (El-Beshbishy et al., 2015; Motawi et al., 2015). However, several studies showed distinct prevalence for these two polymorphisms, particularly the XbaI alleles (Ganasyam et al., 2012). Discrepancies in allelic frequencies between reports may be attributed to the ethnic variability of the studied populations, heterogeneity of the analyzed diseases and sample size.

Most often, deviation from the Hardy-Weinberg equilibrium (HWE) is caused by a small sample size and poor genotyping quality. In this study, minimization of genotyping errors was achieved by amplicon based NGS. It provides a lower cost per sample alternative to restriction fragment length polymorphism (RFLP), probe based Taq-man PCR and Sanger sequencing. RFLP, a traditional genotyping method, requires a high quality, and quantity of DNA sample and numerous steps including amplification, digestion with restriction enzymes and gel electrophoresis which can be a laborious and time-consuming process. Moreover, confirmatory sequencing must be performed when a sample shows confused banding pattern particularly the bands of low molecular size. In our analysis, DNA sequences with a quality score of more than 20 (>Q20), Q20 represents an error rate of 1 in 100 with a corresponding call accuracy of 99%, were selected. Amplicon based NGS requires advanced instrumentation, however and despite this limitation, we believe that it is the method of choice for genotyping studies and for rapid genetic screening of several hundred to thousands of samples.

To our knowledge, this is the first report to investigate whether *ESR1* PvuII and *XbaI* alleles and their genotypes and haplotypes are associated with T2DM among Palestinian women.

In this study, the additive, dominant, co dominant, and recessive genotypic models of association were tested, all were adjusted for age and BMI. The *ESR1* PvuII variant, supported evidence for association with risk of T2DM in three genotypic models
(*P* < 0.025) but the strongest association was observed in the overdominant model, in which heterozygous carriers confer a higher risk compared to homozygotes (Table 3). Although such association toward the heterozygotes could result from genotyping error, this is unlikely to be the case in our study as we used variant specific probes as described above.

One possible explanation is that the variant allele may have a very strong dominant effect so that there is little difference between the effects of the variant homozygotes and heterozygotes. Given that the majority of participants in both cases and controls were heterozygotes, the real effect among the variant homozygotes can only be assessed by much larger study in the future. On the other hand, the over dominant model (upon Bonferroni correction) showed a trend for association of XbaI SNP with T2DM (OR = 2.03, CI [1.05–3.95], *P* = 0.035) which was unexpected given the tight LD between these two polymorphisms. These results could also be attributed to the poor statistical power due to small sample size. Considering the associations of both SNPs found in other studies, we believe that further study involving a larger number of participants is needed to elucidate the exact role of these polymorphisms in susceptibility to T2DM. However, in 2018, a meta-analysis including eight studies indicated that PvuII, rather than XbaI polymorphism, was associated with T2DM (OR = 0.673, 95% CI [0.550–0.823]). In that study, the C allele of PvuII polymorphism showed a protective role in T2DM in Chinese people while the G allele of XbaI polymorphism is related to a reduced risk for T2DM in Caucasian population (Yang et al., 2018). Thus, we believe that physiological pathway to diabetes may vary among different population and those differences are reflected in part by genetic differences. Therefore, the results of our study should be considered exploratory and confirmed by additional studies including other polymorphisms/genes as one single gene polymorphisms are not enough to answer the pathogenesis of a complex disease such as T2DM. Further, it is well known that estrogens regulate the cardiovascular system via their direct effects on the vessel wall and indirect effects on total cholesterol and triglycerides metabolism (Liping et al., 2013). Among diabetic group, stratified analyses based on *ESR1* PvuII and XbaI genotypes revealed no evidence for association with lipid levels (TC, TG, HDL, LDL). Moreover, we did not find any association between the two polymorphisms (PvuII and XbaI) and TC in control group (n = 114). One limitation of this study is the lack of information about subjects who received hormone replacement therapy and/or lipid lowering agent which may affect their serum lipid profile, therefore, we were unable to study the interaction between these variants and hormone replacement therapy on lipid profile which was, however, beyond the scope of this study. It is reported that the effect of PvuII and XbaI polymorphisms on serum lipids depend on several other genes not on *ESR1* gene alone. Similarly, *Matsubara et al.* (1997) and *Almeida et al.* (2006) did not find any association between the two polymorphisms and serum lipid profile in post-menopausal women (*Almeida et al.*, 2006; *Matsubara et al.*, 1997). A recent study conducted on a cohort of post-menopausal Brazilian women reported no effects of PvuII SNP of *ESR1* gene on patient's serum TC, LDL, HDL, and TG, while XbaI SNP was associated to changes in TG and total lipids mainly in obese and overweight woman (*Gomes-Rochette et al.*, 2017). In contrast,
Egyptian study showed that both PvuII and XbaI SNPs have been associated with increased levels of triglycerides, total cholesterol, and LDL (Motawi et al., 2015).

The mechanism by which intronic polymorphisms of *ESR1* gene might confer increased risk of T2DM is not fully understood. However, it is reported that the PvuII and XbaI polymorphisms change the *ESR1* gene expression by altering the binding of its own transcription factors. Moreover, LD between these two polymorphisms with other polymorphisms in the *ESR1* gene, such as TA tandem polymorphism in the promoter region could affect gene expression or function (Herrington et al., 2002). Polymorphisms in *ESR1* were first thought to be potential risk factors for the development of CVD, but a meta-analysis including 10 case-control studies demonstrated that the *ESR1* SNPs PvuII and XbaI are not associated with risk of CVD (Morselli et al., 2017). In this study we have noted an inverse trend between these polymorphisms and risk of CVD in diabetic group. Our findings cannot be a definitive conclusion to identify *ESR1* genotype/haplotype that could be useful in identification of diabetic women that are more prone to develop CVD due to small sample size with limited statistical power to detect interactions or perform subgroups analyses which was the most obvious limitation of this study.

## CONCLUSIONS

The present study suggests that the PvuII polymorphism in the *ESR1* gene is associated with increased risk of type 2 diabetes but not with lipid profile. The interaction between XbaI genotype and T2DM still need to be clarified. Identification of genetic risk markers can be used in defining high risk subjects and in prevention trials. The mechanism by which the *ESR1* PvuII and XbaI polymorphisms might be related to CVD-among diabetic group-needs further investigation with larger sample size. Also, future studies are warranted to replicate our results and to explain the influence of these intronic polymorphisms on diabetes susceptibility in Palestine.

## ACKNOWLEDGEMENTS

The authors gratefully acknowledge the UNRWA outpatient clinic members who helped to the patients' recruitment, and all study participants.

## ADDITIONAL INFORMATION AND DECLARATIONS

### Funding

The authors received no funding for this work.

### Competing Interests

The authors declare that they have no competing interests.

### Author Contributions

* Suheir Ereqat conceived and designed the experiments, analyzed the data, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.
* Stéphane Cauchi conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, approved the final draft.

* Khaled Eweidat performed the experiments, approved the final draft, collect the blood samples and extract the DNA.

* Muawiyah Elqadi performed the experiments, approved the final draft, collect the data from patients medical records.

* Abedelmajeed Nasereddin conceived and designed the experiments, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.

## Human Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

The study protocol was approved by Al-Quds University Research Ethics Committee (Rf no. 2/SRC/4).

## Data Availability

The following information was supplied regarding data availability:

Anonymized raw data are available in a [Supplemental File](#).

## Supplemental Information

Supplemental information for this article can be found online at [http://dx.doi.org/10.7717/peerj.7164#supplemental-information](http://dx.doi.org/10.7717/peerj.7164#supplemental-information).

## REFERENCES

**Alberti KG, Zimmet PZ. 1998.** Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. *Diabetic Medicine* **15**(7):539-553
DOI 10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S.

**Almeida S, Fiegenbaum M, De Andrade FM, Osorio-Wender MC, Hutz MH. 2006.** ESR1 and APOE gene polymorphisms, serum lipids, and hormonal replacement therapy. *Maturitas* **54**(2):119-126 DOI 10.1016/j.maturitas.2005.09.009.

**Bolton JL. 2016.** Menopausal hormone therapy, age, and chronic diseases: perspectives on statistical trends. *Chemical Research in Toxicology* **29**(10):1583-1590
DOI 10.1021/acs.chemrestox.6b00272.

**Damiri B, Abualsoud MS, Samara AM, Salameh SK. 2018.** Metabolic syndrome among overweight and obese adults in Palestinian refugee camps. *Diabetology & Metabolic Syndrome* **10**(1):34 DOI 10.1186/s13098-018-0337-2.

**El-Beshbishy HA, Tawfeek MA, Al-Azhary NM, Mariah RA, Habib FA, Aljayar L, Alahmadi AF. 2015.** Estrogen receptor alpha (ESR1) gene polymorphisms in pre-eclamptic saudi patients. *Pakistan Journal of Medical Sciences* **31**(4):880-885 DOI 10.12669/pjms.314.7541.

**El Kishawi RR, Soo KL, Abed YA, Muda WA. 2014.** Obesity and overweight: prevalence and associated socio demographic factors among mothers in three different areas in the Gaza Strip-Palestine: a cross-sectional study. *BMC Obesity* **1**(1):7 DOI 10.1186/2052-9538-1-7.
**Ganasyam SR, Rao TB, Murthy YS, Jyothy A, Sujatha M. 2012.** Association of estrogen receptor-α gene & metallothionein-1 gene polymorphisms in type 2 diabetic women of Andhra Pradesh. *Indian Journal of Clinical Biochemistry* **27**(1):69-73 DOI 10.1007/s12291-011-0179-2.

**Gomes-Rochette NF, Souza LS, Tommasi BO, Pedrosa DF, Eis SR, Fin ID, Vieira FL, Graceli JB, Rangel LB, Silva IV. 2017.** Association of PvuII and XbaI polymorphisms on estrogen receptor alpha (ESR1) gene to changes into serum lipid profile of post-menopausal women: Effects of aging, body mass index and breast cancer incidence. *PLOS ONE* **12**(2):e0169266 DOI 10.1371/journal.pone.0169266.

**Gonzalez JR, Armengol L, Sole X, Guino E, Mercader JM, Estivill X, Moreno V. 2007.** SNPassoc: an R package to perform whole genome association studies. *Bioinformatics* **23**(5):644-645 DOI 10.1093/bioinformatics/btm025.

**Gorres BK, Bomhoff GL, Morris JK, Geiger PC. 2011.** In vivo stimulation of oestrogen receptor alpha increases insulin-stimulated skeletal muscle glucose uptake. *Journal of Physiology* **589**(8):2041-2054 DOI 10.1113/jphysiol.2010.199018.

**Gupte AA, Pownall HJ, Hamilton DJ. 2015.** Estrogen: an emerging regulator of insulin action and mitochondrial function. *Journal of Diabetes Research* **2015**(9):916585 DOI 10.1155/2015/916585.

**Heine PA, Taylor JA, Iwamoto GA, Lubahn DB, Cooke PS. 2000.** Increased adipose tissue in male and female estrogen receptor-alpha knockout mice. *Proceedings of the National Academy of Sciences of the United States of America* **97**(23):12729-12734 DOI 10.1073/pnas.97.23.12729.

**Herrington DM, Howard TD, Brosnihan KB, McDonnell DP, Li X, Hawkins GA, Reboussin DM, Xu J, Zheng SL, Meyers DA, Bleecker ER. 2002.** Common estrogen receptor polymorphism augments effects of hormone replacement therapy on E-selectin but not C-reactive protein. *Circulation* **105**(16):1879-1882 DOI 10.1161/01.cir.0000016173.98826.88.

**Hevener AL, Zhou Z, Moore TM, Drew BG, Ribas V. 2018.** The impact of ERalpha action on muscle metabolism and insulin sensitivity—Strong enough for a man, made for a woman. *Molecular Metabolism* **15**:20-34 DOI 10.1016/j.molmet.2018.06.013.

**Huang Q, Wang TH, Lu WS, Mu PW, Yang YF, Liang WW, Li CX, Lin GP. 2006.** Estrogen receptor alpha gene polymorphism associated with type 2 diabetes mellitus and the serum lipid concentration in Chinese women in Guangzhou. *Chinese Medical Journal* **119**(21):1794-1801 DOI 10.1097/00029330-200611010-00006.

**Kumar R, Nandhini LP, Kamalanathan S, Sahoo J, Vivekanadan M. 2016.** Evidence for current diagnostic criteria of diabetes mellitus. *World Journal of Diabetes* **7**(17):396-405 DOI 10.4239/wjd.v7.i17.396.

**Liping D, Lihua H, Zhitao J, Taohong H, Huili M, Lina Z, Xincheng Q, Caiyi L. 2013.** A meta-analysis of correlation of ER gene polymorphisms and risk in Chinese population with coronary heart disease. *Life Science Journal* **10**(4):595-598.

**Lizcano F, Guzman G. 2014.** Estrogen deficiency and the origin of obesity during Menopause. *BioMed Research International* **2014**:757461 DOI 10.1155/2014/757461.

**Matsubara Y, Murata M, Kawano K, Zama T, Aoki N, Yoshino H, Watanabe G, Ishikawa K, Ikeda Y. 1997.** Genotype distribution of estrogen receptor polymorphisms in men and postmenopausal women from healthy and coronary populations and its relation to serum lipid levels. *Arteriosclerosis, Thrombosis, and Vascular Biology* **17**(11):3006-3012 DOI 10.1161/01.atv.17.11.3006.

**Mayo O. 2008.** A century of Hardy-Weinberg equilibrium. *Twin Research and Human Genetics* **11**(3):249-256 DOI 10.1375/twin.11.3.249.
**Meyer MR, Clegg DJ, Prossnitz ER, Barton M. 2011.** Obesity, insulin resistance and diabetes: sex differences and role of oestrogen receptors. *Acta Physiologica 203*(1):259-269
DOI 10.1111/j.1748-1716.2010.02237.x.

**Mohammadi F, Pourahmadi M, Mosalanejad M, Jamali H, Ghobadifar MA, Erfanian S. 2013.** Association of estrogen receptor α genes PvuII and XbaI polymorphisms with type 2 diabetes mellitus in the inpatient population of a hospital in Southern Iran. *Diabetes & Metabolism Journal* 37(4):270-277 DOI 10.4093/dmj.2013.37.4.270.

**Morselli E, Santos RS, Criollo A, Nelson MD, Palmer BF, Clegg DJ. 2017.** The effects of oestrogens and their receptors on cardiometabolic health. *Nature Reviews Endocrinology* 13(6):352-364 DOI 10.1038/nrendo.2017.12.

**Motawi TM, El-Rehany MA, Rizk SM, Ramzy MM, El-Roby DM. 2015.** Genetic polymorphism of estrogen receptor alpha gene in Egyptian women with type II diabetes mellitus. *Meta Gene* 6:36-41 DOI 10.1016/j.mgene.2015.08.001.

**Onland-Moret NC, Van Gils CH, Roest M, Grobbee DE, Peeters PH. 2005.** The estrogen receptor α gene and breast cancer risk (The Netherlands). *Cancer Causes & Control* 16(10):1195-1202 DOI 10.1007/s10552-005-0307-5.

**Palestinian Ministry of Health. 2015.** Health annual report: Palestine. Available at http://www.site.moh.ps/Content/Books/ZxC cynmiUofNqt66u4CrHRgmJR6Uv7z77srjjIEAho6xnz5V3rgLTu_RhO7xf2j2VusNiIvWkjwp84yXHLdGleB97gKrHHI5iZ9oPJ25owGEN.pdf (accessed 5 January 2019).

**Peterson RA. 2017.** Estimating normalization transformations with bestNormalize. GitHub. Available at https://github.com/petersonR/bestNormalize.

**Ren Y, Zhang M, Liu Y, Sun X, Wang B, Zhao Y, Liu D, Liu X, Zhang D, Liu F, Cheng C, Liu L, Chen X, Zhou Q, Hu D. 2018.** Association of menopause and type 2 diabetes mellitus. *Menopause* 26(3):325-330 DOI 10.1097/GME.0000000000001200.

**Ribas V, Nguyen MT, Henstridge DC, Nguyen AK, Beaven SW, Watt MJ, Hevener AL. 2010.** Impaired oxidative metabolism and inflammation are associated with insulin resistance in ERα-deficient mice. *American Journal of Physiology-Endocrinology and Metabolism* 298(2):E304-E319 DOI 10.1152/ajpendo.00504.2009.

**Sabarneh A, Ereqat S, Cauchi S, AbuShamma O, Abdelhafez M, Ibrahim M, Nasereddin A. 2018.** Common FTO rs9939609 variant and risk of type 2 diabetes in Palestine. *BMC Medical Genetics* 19(1):156 DOI 10.1186/s12881-018-0668-8.

**Sale MM, Freedman BI, Langefeld CD, Williams AH, Hicks PJ, Colicigno CJ, Beck SR, Brown WM, Rich SS, Bowden DW. 2004.** A genome-wide scan for type 2 diabetes in African-American families reveals evidence for a locus on chromosome 6q. *Diabetes* 53(3):830-837 DOI 10.2337/diabetes.53.3.830.

**Salpeter SR, Walsh JM, Ormiston TM, Greyber E, Buckley NS, Salpeter EE. 2006.** Meta-analysis: effect of hormone-replacement therapy on components of the metabolic syndrome in postmenopausal women. *Diabetes, Obesity and Metabolism* 8(5):538-554
DOI 10.1111/j.1463-1326.2005.00545.x.

**Shahin Y, Kapur A, Seita A. 2015.** Diabetes care in refugee camps: the experience of UNRWA. *Diabetes Research and Clinical Practice* 108(1):1-6 DOI 10.1016/j.diabres.2015.01.035.

**Silva IV, Rezende LCD, Lanes SP, Souza LS, Madeira KP, Cerri MF, Paes MF, Daltoe RD, Chambo-Filho A, Guimaraes MC, Graceli JB, Rangel LBA. 2010.** Evaluation of PvuII and XbaI polymorphisms in the estrogen receptor alpha gene (ESR1) in relation to menstrual cycle timing and reproductive parameters in post-menopausal women. *Maturitas* 67(4):363-367
DOI 10.1016/j.maturitas.2010.08.006.
Slominski B, Mysliwska J, Ryba-Stanislawowska M, Skrzypkowska M, Mysliwiec M. 2018. Estrogen receptor α gene polymorphism and vascular complications in girls with type 1 diabetes mellitus. *Molecular and Cellular Biochemistry* **437**(1-2):153-161
DOI 10.1007/s11010-017-3103-0.

Speer G, Cseh K, Winkler G, Vargha P, Braun E, Takacs I, Lakatos P. 2001. Vitamin D and estrogen receptor gene polymorphisms in type 2 diabetes mellitus and in android type obesity. *European Journal of Endocrinology* **144**:385-389 DOI 10.1530/eje.0.1440385.

Weng H, Zhang C, Hu YY, Yuan RX, Zuo HX, Yan JZ, Niu YM. 2015. Association between estrogen receptor-α gene XbaI and PvuII polymorphisms and periodontitis susceptibility: a meta-analysis. *Disease Markers* **2015**:741972 DOI 10.1155/2015/741972.

World Health Organization (WHO). 2016. Global report on diabetes. Available at https://apps.who.int/iris/bitstream/handle/10665/204871/9789241565257_eng.pdf?sequence=1.

Yang J, Han R, Chen M, Yuan Y, Hu X, Ma Y, Wu M, Zhang X, Wang M, Jiang S, Pan F. 2018. Associations of estrogen receptor alpha gene polymorphisms with type 2 diabetes mellitus and metabolic syndrome: a systematic review and meta-analysis. *Hormone and Metabolic Research* **50**(6):469-477 DOI 10.1055/a-0620-8553.

Zhao L, Fan X, Zuo L, Guo Q, Su X, Xi G, Zhang Z, Zhang J, Zheng G. 2018. Estrogen receptor 1 gene polymorphisms are associated with metabolic syndrome in postmenopausal women in China. *BMC Endocrine Disorders* **18**(1):65 DOI 10.1186/s12902-018-0289-4.