ESR1 polymorphism (rs2234693) influences femoral bone mass in patients with Turner syndrome

in Endocrine Connections
Authors:
Renata C Scalco Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil
Disciplina de Endocrinologia, Faculdade de Ciências Médicas da Santa Casa de São Paulo, São Paulo, São Paulo, Brazil
Unidade de Endocrinologia do Desenvolvimento, Laboratorio de Hormonios e Genetica Molecular LIM/42, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Renata C Scalco in
Current site
Google Scholar
PubMed
Close
,
Ericka B Trarbach Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Ericka B Trarbach in
Current site
Google Scholar
PubMed
Close
,
Edoarda V A Albuquerque Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Edoarda V A Albuquerque in
Current site
Google Scholar
PubMed
Close
,
Thais K Homma Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Thais K Homma in
Current site
Google Scholar
PubMed
Close
,
Thais H Inoue-Lima Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Thais H Inoue-Lima in
Current site
Google Scholar
PubMed
Close
,
Mirian Y Nishi Unidade de Endocrinologia do Desenvolvimento, Laboratorio de Hormonios e Genetica Molecular LIM/42, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Mirian Y Nishi in
Current site
Google Scholar
PubMed
Close
,
Berenice B Mendonca Unidade de Endocrinologia do Desenvolvimento, Laboratorio de Hormonios e Genetica Molecular LIM/42, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Berenice B Mendonca in
Current site
Google Scholar
PubMed
Close
, and
Alexander A L Jorge Unidade de Endocrinologia Genetica, Laboratorio de Endocrinologia Celular e Molecular LIM/25, Disciplina de Endocrinologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, São Paulo, Brazil

Search for other papers by Alexander A L Jorge in
Current site
Google Scholar
PubMed
Close

Correspondence should be addressed to R C Scalco: rescalco@gmail.com
Open access

Sign up for journal news

Most patients with Turner syndrome (TS) need hormone replacement therapy because of hypergonadotropic hypogonadism; individual outcomes, however, are highly variable. Our objective was to assess the influence of five estrogen receptor 1 gene (ESR1) polymorphisms (rs543650, rs1038304, rs2046210, rs2234693 and rs9340799) on adult height, breast development, uterine volume and bone mineral density (BMD). We studied 91 TS patients from a tertiary hospital using adult estrogen dose. In our group, ESR1 rs2234693 was associated with femoral neck and total hip BMD, and it accounted for around 10% of BMD variability in both sites (P < 0.01). Patients homozygous for C allele in this polymorphism had significantly lower femoral neck BMD (0.699 ± 0.065 g/cm2 vs 0.822 ± 0.113 g/cm2, P = 0.008) and total hip BMD (0.777 ± 0.118 g/cm2 vs 0.903 ± 0.098 g/cm2, P = 0.009) than patients homozygous for T allele. The other four ESR1 polymorphisms were not able to predict any of the above estrogen therapy outcomes in an isolated manner. Patients homozygous for the haplotype GCG formed by polymorphisms rs543650, rs2234693 and rs9340799 had an even more significantly lower femoral neck BMD (0.666 ± 0.049 vs 0.820 ± 0.105 g/cm2, P = 0.0047) and total hip BMD (0.752 ± 0.093 vs 0.908 ± 0.097 g/cm2, P = 0.0029) than patients homozygous for haplotypes with a T allele in rs2234693. In conclusion, homozygosity for C allele in ESR1 rs2234693 and/or for GCG haplotype appears to be associated with lower femoral neck and total hip BMD. We believe that the identification of polymorphisms related to estrogen outcomes may contribute to individualization of treatment in TS.

Abstract

Most patients with Turner syndrome (TS) need hormone replacement therapy because of hypergonadotropic hypogonadism; individual outcomes, however, are highly variable. Our objective was to assess the influence of five estrogen receptor 1 gene (ESR1) polymorphisms (rs543650, rs1038304, rs2046210, rs2234693 and rs9340799) on adult height, breast development, uterine volume and bone mineral density (BMD). We studied 91 TS patients from a tertiary hospital using adult estrogen dose. In our group, ESR1 rs2234693 was associated with femoral neck and total hip BMD, and it accounted for around 10% of BMD variability in both sites (P < 0.01). Patients homozygous for C allele in this polymorphism had significantly lower femoral neck BMD (0.699 ± 0.065 g/cm2 vs 0.822 ± 0.113 g/cm2, P = 0.008) and total hip BMD (0.777 ± 0.118 g/cm2 vs 0.903 ± 0.098 g/cm2, P = 0.009) than patients homozygous for T allele. The other four ESR1 polymorphisms were not able to predict any of the above estrogen therapy outcomes in an isolated manner. Patients homozygous for the haplotype GCG formed by polymorphisms rs543650, rs2234693 and rs9340799 had an even more significantly lower femoral neck BMD (0.666 ± 0.049 vs 0.820 ± 0.105 g/cm2, P = 0.0047) and total hip BMD (0.752 ± 0.093 vs 0.908 ± 0.097 g/cm2, P = 0.0029) than patients homozygous for haplotypes with a T allele in rs2234693. In conclusion, homozygosity for C allele in ESR1 rs2234693 and/or for GCG haplotype appears to be associated with lower femoral neck and total hip BMD. We believe that the identification of polymorphisms related to estrogen outcomes may contribute to individualization of treatment in TS.

Introduction

Turner syndrome (TS) is characterized by complete or partial absence of the second sex chromosome in phenotypic females with one or more clinical manifestations (1). It is one of the most common chromosomal abnormalities, affecting approximately 25–50 per 100,000 females. Hypergonadotropic hypogonadism is a frequent manifestation, present in more than 95% of TS patients. Therefore, most patients will need hormone replacement therapy (HRT) for puberty induction, attaining of peak bone mass and maintenance of female secondary sex characteristics (1).

However, individual outcomes of estrogen treatment are highly variable among these patients. For instance, uterine development descriptions range from lower than average even with adequate HRT (2) to normal in all treated patients (3). This feature was previously associated with diverse factors such as the presence of spontaneous puberty (4), duration of HRT (4, 5), dose of HRT (5, 6), type of HRT (4) and karyotype (7). Nevertheless, there are controversies among studies and none of these factors can thoroughly explain the interindividual variation.

Genome-wide association studies (GWAS) showed that polymorphisms in estrogen receptor 1 gene (ESR1) participate in the variability of many characteristics, such as height, bone density, breast size and risk of cancer (8, 9, 10, 11). There is only one study involving TS patients that showed an association between two ESR1 polymorphisms and bone mass gain with estrogen therapy (12). To our knowledge, there are no data on the impact of ESR1 polymorphisms on height variability or on the development of secondary sex characteristics in these patients. Thus, the aim of this study was to assess the influence of ESR1 polymorphisms on estrogen therapy outcomes in TS patients.

Patients and methods

Subjects

The study protocol was approved by the Research Ethics Committee of the Hospital das Clinicas da Faculdade de Medicina da Universidade de Sao Paulo. Written informed consent was obtained from all patients before starting the molecular studies. Ninety-one TS patients from our endocrinology clinic were selected using the following criteria: (1) regular use of adult estrogen dose (i.e. at least 1 mg of estradiol valerate daily or equivalent) for a minimum period of one year; (2) absence of spontaneous menarche; and (3) absence of Y chromosome material in karyotype. These patients were submitted to clinical evaluation for height and breast development (according to Tanner staging), pelvic ultrasound for assessment of uterine volume, and dual-energy X-ray absorptiometry (DEXA) for determination of bone mineral density at lumbar spine and proximal femur. Since DEXA has a known tendency to underestimate the bone density of people with short stature, we calculated bone mineral apparent density (BMAD) by using bone mineral content and projected bone area (13).

Cytogenetic study

Karyotype data were collected from medical records. Twenty to fifty cells were evaluated. Metaphase spreads prepared from peripheral lymphocyte cultures from each patient were studied by conventional staining and G and C banding techniques and were analyzed through an optical microscope (Carl Zeiss). Patients with chromosome fragments or rings of unidentified origin were evaluated for the presence of Y-chromosomal material with real-time PCR with multiple Y-specific probes, according to a previously described protocol (5).

ESR1 polymorphisms

ESR1 polymorphisms were selected through active search in PubMed and in GWAS catalogs. We chose five polymorphisms that were recurrent in studies or that were Tag-single nucleotide polymorphisms (Tag-SNPs), provided that they had a minor allele frequency (MAF) which allowed the present study with the available sample. Two of them are well-known ESR1 polymorphisms (rs2234693 and rs9340799) and the others were previously found to influence height (rs543650) (8), bone mineral density (rs1038304) (9) or risk of breast cancer (rs2046210) (11) in GWAS.

Molecular studies

Genomic DNA was isolated from peripheral blood leukocytes by standard methods from all patients. The ESR1 polymorphisms were genotyped by allelic discrimination using specific TaqMan probes and primers in the Step One Plus Real-Time PCR system (Applied Biosystems) according to manufacturer’s instructions. Ten percent of all samples were randomly selected and sequenced by Sanger for quality control. Briefly, PCR amplification was performed using GoTaq Green PCR systems (Promega) with 25 ng of DNA, 1× buffer, 10 mM dNTP, 0.2 mM of each forward and reverse primers and 1 U of enzyme, in a final 15 µL reaction mixture (Table 1). PCR cycling conditions were as follows: 98°C for 5 min, followed by 30 cycles of 95°C for 30 min, 56°C of annealing temperature for 45 min and 72°C for 30 min; and 10 min of final elongation step at 72°C. Subsequently, the PCR products were sequenced using BigDye Terminator v3.1 cycle sequencing protocol in the ABI PRISM 310 automatic sequencer (Thermo Fisher). The agreement of genotyping between both methods was 100%.

Table 1

List of TaqMan assays and PCR primers used in this study.

ESR1 SNPs TaqMan assay ID PCR primer sequences (5′ to 3′) PCR amplicon (pb)
Forward Reverse
rs2234693a C_3163590-10 CATGAACCACCATGCTCAGT AAAAACATACTACCTGCACCAGAA 388
rs9340799a C_3163591-10
rs543650 C_551076-10 TCCCTGGAGGAAACGTAAAA TCCCAACACCTTCCAGACTC 360
rs1038304 C_219719-10 CCTGAGTAGCCGGGAGTACA AGTCTGAGGAATGGGAGCAA 359
rs2046210 C12034236-10 GTGGGTCAAGACCAGCATTT CCATCGTCCACATCTCACAC 337

aBoth SNPs were PCR genotyped in the same amplicon.

Statistical analysis

In this study, we evaluated the association among ESR1 polymorphisms and the following outcomes: adult height, breast development, uterine volume and bone mineral density at lumbar spine and proximal femur. Qualitative variables are listed as percentages, whereas quantitative variables are shown as mean ± s.d. Nominal variables were compared with chi-square test or Fisher exact test, as appropriate. Numerical variables were analyzed using t-test, Mann–Whitney rank-sum test or Kruskal–Wallis one-way ANOVA on ranks, as appropriate. To assess whether genetic variants had independent prognostic significance for outcomes, we performed linear and multiple-regression analyses adjusting for established influential factors. Since we tested associations with clinical parameters for multiple SNPs (n = 5), we applied Bonferroni’s correction for the P values and only associations with P value <0.01 were considered significant. Haplotype blocks were constructed by using the Solid Spine method from the software Haploview (14) and their association with estrogen outcomes was also assessed. All statistical analyses were performed using SigmaStat for Windows (version 3.5).

Results

We evaluated 91 TS patients aged 25.1 ± 6.8 years old (range 16 to 42 years old) (Table 2). 45,X karyotype was the most frequently found (n = 55, 60%), followed by 46,XiXq with or without mosaicism (n = 12), other structural abnormalities (n = 7), 45,X/46,XX (n = 6), ring X chromosome mosaicisms (n = 6), 45,X/47,XXX (n = 3), marker chromosome mosaicisms (n = 2). Eighty-three percent of them received recombinant human growth hormone (rhGH) treatment for a mean duration of 6.3 ± 3.1 years, starting at 9.9 ± 4.0 years old. Their mean adult height was 148.4 ± 6.9 cm. Their age at puberty onset or induction was 14.3 ± 1.9 years old and 9% of them had spontaneous thelarche. Uterine volume had a wide range, from 5.9 to 107 cm3 (mean 31.3 ± 18.4 cm3). Ninety percent of patients reached Tanner stage 5 breast development. Eighty percent of patients had normal bone mass after calculation of BMAD based on bone mineral content and projected bone area.

Table 2

Clinical and radiological data from 91 Turner syndrome patients.

Basal characteristics
 45,X karyotype 60%
 Target height 158.6 ± 4.7 cm
 Use of rhGH 83%
 Age at rhGH onset 9.9 ± 4.0 years
 Duration of rhGH use 6.3 ± 3.1 years
 Age at puberty onset or induction 14.3 ± 1.9 years old
 Spontaneous thelarche 9%
Current data
 Age 25.1 ± 6.8 years old
 Adult height 148.4 ± 6.9 cm
 Breast development (Tanner staging) B5 90%
 Uterine volume 31.3 ± 18.4 cm3
 Normal bone mass in lumbar spine (BMD) 68%
 Normal bone mass in lumbar spine (BMAD) 80%
 Normal bone mass in femoral neck (BMD) 80%
 Normal bone mass in femoral neck (BMAD) 96.5%

Qualitative variables are listed as percentages, whereas quantitative variables are shown as mean ± s.d.

B5, Tanner 5 breast development; BMAD, bone mineral apparent density; BMD, bone mineral density; rhGH, recombinant human growth hormone.

All patients were genotyped for the five selected ESR1 polymorphisms: rs543650, rs1038304, rs2046210, rs2234693 and rs9340799. We did not detect deviations from the Hardy-Weinberg equilibrium in their genotype distributions. We first evaluated the association among these polymorphisms and adult height, uterine volume, breast development and bone density by using linear regression. ESR1 rs2234693, also known as PvuII (after the name of the restriction enzyme initially used to study it) was associated with femoral neck and total hip BMD, and it accounted for around 10% of BMD variability in these sites (Table 3). This polymorphism was not related to adult height, uterine volume, breast development and lumbar spine BMD or BMAD in our patients. The other four ESR1 polymorphisms were not able to predict any of the above estrogen therapy outcomes in an isolated manner.

Table 3

Linear regressions among five ESR1 polymorphisms and selected estrogen replacement outcomes (P values).

ESR1 polymorphisms rs543650 rs1038304 rs2046210 rs2234693 rs9340799
Adult height 0.318 0.89 0.831 0.129 0.218
Breast development 0.719 0.1 0.154 0.495 0.687
Femoral neck BMD 0.558 0.789 0.395 0.007 (R2 = 0.107) 0.16
Femoral neck BMAD 0.768 0.129 0.266 0.025 0.352
Lumbar spine BMD 0.936 0.063 0.074 0.364 0.158
Lumbar spine BMAD 0.289 0.375 0.301 0.280 0.437
Total hip BMD 0.958 0.908 0.606 0.009 (R2 = 0.108) 0.08
Uterine volume 0.998 0.643 0.445 0.171 0.851

Only associations with P < 0.01 were considered significant after application of Bonferroni’s correction.

BMAD, bone mineral apparent density; BMD, bone mineral density; ESR1, estrogen receptor 1 gene; ns, not significant.

ESR1 rs2234693 distribution was C/C 14%, C/T 57% and T/T 29%. Genotypic groups were similar concerning karyotype distribution, target height, use of rhGH treatment, age at rhGH onset and duration of treatment, age at puberty onset/induction and frequency of spontaneous thelarche (Table 4). There were also no significant differences concerning estrogen dose, adult height, uterine volume, breast development according to Tanner staging or lumbar spine BMD and BMAD. However, femoral neck and total hip BMD were significantly different among the groups, and these differences were due to a significantly lower BMD in both sites in patients homozygous for C allele in comparison to patients homozygous for T allele.

Table 4

Clinical and radiological data from 91 Turner syndrome patients grouped according to ESR1 rs2234693 genotype.

C/C C/T T/T C/C vs C/T vs T/Ta CC vs TTb
P value P value
Number of patients 13 52 26
45,X karyotype 61% 65% 50% 0.164 0.094
Target height (cm) 159.9 ± 4.2 158.9 ± 5.1 157.2 ± 4.0 0.216 0.068
Use of rhGH 92% 83% 81% 0.641 0.365
Age at rhGH onset (years) 11.1 ± 4.1 9.7 ± 4.0 9.2 ± 3.9 0.463 0.221
Duration of rhGH use (years) 5.7 ± 3.3 6.3 ± 3.0 6.9 ± 3.4 0.50 0.311
Age at puberty onset or induction (years) 14.7 ± 1.6 14.2 ± 2.0 14.1 ± 1.9 0.609 0.323
Spontaneous thelarche 7.7% 9.6% 7.7% 0.96 0.908
Adult height (cm) 151 ± 5.4 148.3 ± 6.8 147.2 ± 7.8 0.279 0.126
Breast development (Tanner) M5 85% M5 86% M5 100% 0.151 0.05
Estradiol valerate or equivalent dose (mg/day) 1.1 ± 0.6 1.4 ± 0.5 1.3 ± 0.4 0.118 0.15
Femoral neck BMAD (g/cm3) 0.150 ± 0.017 0.167 ± 0.028 0.181 ± 0.037 0.08 0.025
Femoral neck BMD (g/cm2) 0.699 ± 0.065 0.760 ± 0.119 0.820 ± 0.110 0.028c 0.007
Lumbar spine BMD (g/cm2) 0.867 ± 0.122 0.898 ± 0.142 0.910 ± 0.124 0.64 0.305
Lumbar spine BMAD (g/cm3) 0.154 ± 0.040 0.146 ± 0.020 0.144 ± 0.014 0.474 0.907
Total hip BMD (g/cm2) 0.777 ± 0.118 0.857 ± 0.111 0.901 ± 0.096 0.027d 0.007
Uterine volume (cc3)e 29.9 ± 17.8

29.5 (8.0–76.2)
29.0 ± 16.5

25.2 (5.9–82.9)
35.6 ± 22.6

22.6 (10.6–107)
0.329 0.35

Only associations with P < 0.01 after application of Bonferroni’s correction were considered significant.

aOne-way ANOVA or chi-square/Fisher. bt test or chi-square/Fisher. cHolm–Sidak method: T/T vs C/C P = 0.0111; T/T vs C/T ns; C/T vs C/C ns. dHolm–Sidak method: T/T vs C/C P = 0.00786; T/T vs C/T ns; C/T vs C/C ns. eFor its wide range, uterine volume is also presented as median (minimum-maximum).

BMAD, bone mineral apparent density; BMD, bone mineral density; ns, not significant; rhGH, recombinant human growth hormone.

ESR1 rs2234693 was located in a linkage disequilibrium block with rs543650 and rs9340799 (Fig. 1). The most common haplotypes were TTA (32.5%, formed by a T allele in ESR1 rs543650, a T allele in rs2234693 and an A allele in rs9340799), followed by GCG (31.1%), GTA (22.9%) and GCA (11.2%). ESR1 rs1038304 and rs2046211 were located in another linkage disequilibrium block nearby. HRT outcomes were compared among haplotypes, and patients homozygous for the haplotype GCG had significantly lower femoral neck BMD and total hip BMD than patients homozygous for haplotypes with a T allele in rs2234693 (Table 5). It was not possible to compare patients homozygous for GCG haplotype with patients homozygous for other haplotypes containing a C allele in rs2234693 because there was only one patient with this latter haplotype. None of the haplotypes was linked to adult height, uterine volume or L1-L4 BMD.

Figure 1
Figure 1

Haploview linkage disequilibrium plot for five selected SNPs within ESR1. Haplotype blocks were constructed by using the Solid Spine method. ESR1 rs2234693 was located in a linkage disequilibrium block with rs543650 and rs9340799. ESR1 rs1038304 and rs2046211 were located in another linkage disequilibrium block nearby.

Citation: Endocrine Connections 8, 11; 10.1530/EC-19-0398

Table 5

Comparison of bone mass data among Turner syndrome patients grouped according to ESR1 haplotypes (involving rs543650, rs2234693 and rs9340799).

Homozygous GCG (n = 8) Heterozygous GCG (n = 38) xCx (except GCG) (n = 20) Homozygous xTx (n = 25) P value all groups P value homozygous GCG vs homozygous xTx
Femoral neck BMD (g/cm2) 0.666 ± 0.049 0.766 ± 0.117 0.775 ± 0.123 0.820 ± 0.105 0.039 0.0047
Lumbar spine BMD (g/cm2) 0.846 ± 0.103 0.884 ± 0.131 0.947 ± 0.147 0.889 ± 0.134 0.224 0.417
Total hip BMD (g/cm2) 0.752 ± 0.093 0.852 ± 0.099 0.869 ± 0.132 0.908 ± 0.097 0.025 0.0029

BMD, bone mineral density; heterozygous GCG, patients heterozygous for the haplotype formed by a G allele in ESR1 rs543650, a C allele in rs2234693 and a G allele in rs9340799; homozygous GCG, patients homozygous for the haplotype formed by a G allele in ESR1 rs543650, a C allele in rs2234693 and a G allele in rs9340799; homozygous xTx, patients homozygous for haplotypes formed by any allele in ESR1 rs543650, a T allele in rs2234693 and any allele in rs9340799; ns, not significant; xCx, patients heterozygous (and 1 homozygous) for haplotypes formed by any allele in ESR1 rs543650, a C allele in rs2234693 and any allele in rs9340799.

Discussion

Polymorphisms influence physiological, morphological and pathological variation in the human population (15). In TS, polymorphisms in the growth hormone receptor gene (GHR) (16, 17) and the insulin-like growth factor-binding protein 3 gene (IGFBP3) (16) were correlated with increased responsiveness to rhGH treatment (16, 17) and adult height (16). One polymorphism in the suppressor of cytokine signaling 2 gene (SOCS2) was also linked to adult height after rhGH treatment (18). Furthermore, polymorphisms in the protein tyrosine phosphatase non-receptor type 2 gene (PTPN22) (19) and the myosin IXB gene (MYO9B) (20) were associated with a higher risk of autoimmune diseases; polymorphisms in the vitamin D receptor gene (VDR) were related to lower BMD and higher levels of bone markers (21); and a polymorphism in the cytotoxic T lymphocyte-associated protein 4 gene (CTLA-4) was linked to obesity in these patients (22). Finally, ESR1 rs2234693 and rs9340799 were associated with bone mass gain in lumbar spine after the onset of estrogen replacement in adult TS patients (12).

In the present study, we found that TS patients homozygous for C allele in ESR1 rs2234693 had lower femoral neck and total hip BMD than TS patients homozygous for T allele. Additionally, patients homozygous for the haplotype GCG formed by polymorphisms rs543650, rs2234693 and rs9340799 had an even more significantly lower femoral neck BMD and total hip BMD than patients homozygous for haplotypes with a T allele in rs2234693. Similar results were found in two other studies involving young premenopausal healthy patients (23, 24). Massart et al. found that, in a group of Italian women aged between 20 and 30 years old, the presence of at least one C allele in rs2234693 was associated with lower lumbar spine and proximal femur BMD (23). Furthermore, McGuigan et al. showed that, in a group of Irish women around 22 years old, homozygosity for T allele in ESR1 rs2234693 was related to higher hip BMD (24). Our results are also aligned to those of the previous study on bone mass gain after the onset of estrogen replacement in TS patients (12), which showed significant bone gain in lumbar spine in patients who had at least one T allele in ESR1 rs2234693 and/or one A allele in ESR1 rs9340799 but not in patients homozygous for C allele in ESR1 rs2234693 and/or homozygous for G allele in ESR1 rs9340799. Unfortunately, proximal femur BMD was not evaluated in that study.

Briefly, estrogens increase bone formation at the periosteum through its action in osteoblast progenitors and suppress bone resorption in trabecular and endocortical bone surfaces by decreasing osteoclast numbers. These actions are mainly mediated by estrogen receptor 1 (ER1) nuclear and non-nuclear initiated signaling pathways. However, the molecular details by which ER1 signaling influences bone mass remain unclear (25). In a in vitro model using renal cells, Herrington et al. showed that the presence of a C allele in ESR1 rs2234693 produces a potential binding site for myb transcription factors in ER1, and might produce ER1 isoforms that have different properties than the full-length gene product (26). If the presence of a C allele in ESR1 rs2234693 reduces ER1 function in bone, this could explain why TS patients who were homozygous for this allele had lower bone mass gain with estrogen replacement (12). Moreover, it could also explain the lower femoral neck and total hip BMD seen in young women with this genotype. The reason for the lack of association between this allele and lumbar spine BMD in our study and in McGuigan’s (24) is still unclear but could be probably explained by the relatively small size of both studies.

Other factors that could have influenced bone mass, such as use and duration of rhGH treatment, age at puberty onset/induction or estrogen dose, were not significantly different among genotypic groups. Karyotype is considered not to be associated with BMD or fracture risk by most groups (1), and was not linked to bone mass in our patients either. Limitations of our study include its small size and absence of a replication analysis in young Brazilian women without TS. Consequently, further studies, preferably prospective and multicenter, with a larger number of TS patients and healthy young controls are necessary to confirm our findings. On the other hand, this was the first study to evaluate the association of an ESR1 polymorphism to bone mass in TS patients already receiving adult-dose estrogen replacement therapy. Moreover, the presence of a similar association in other groups of young women with and without TS strengthens our result and suggests a potential contribution of ESR1 rs2234693 to bone mass acquisition. Finally, although we could not find associations among ESR1 polymorphisms and adult height, breast development and uterine volume, the high variability especially of the latter suggests that genetic variants play an important role. More efforts to identify polymorphisms associated with estrogen outcomes may allow for individualization of treatment in TS.

Declaration of interest

Merck has performed a scientific review of the publication, but the views and opinions described in the publication do not necessarily reflect those of Merck.

Funding

This study was supported by Grant for Growth Innovation from Merck KGaA, Darmstadt, Germany and by Grants 2016/03318-0 (to R C S) and 2013/03236-5 (to A A L J) from São Paulo Research Foundation (FAPESP).

References

  • 1

    Gravholt CH, Andersen NH, Conway GS, Dekkers OM, Geffner ME, Klein KO, Lin AE, Mauras N, Quigley CA, Rubin K, et al.Clinical practice guidelines for the care of girls and women with Turner syndrome: proceedings from the 2016 Cincinnati International Turner Syndrome Meeting. European Journal of Endocrinology 2017 G1G70. (https://doi.org/10.1530/EJE-17-0430)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 2

    Cleemann L, Holm K, Fallentin E, Skouby SO, Smedegaard H, Møller N, Borch-Christensen H, Jeppesen EM, Wieslander SB, Andersson AM, et al.Uterus and ovaries in girls and young women with Turner syndrome evaluated by ultrasound and magnetic resonance imaging. Clinical Endocrinology 2011 756761. (https://doi.org/10.1111/j.1365-2265.2011.03995.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 3

    McDonnell CM, Coleman L, Zacharin MR. A 3-year prospective study to assess uterine growth in girls with Turner’s syndrome by pelvic ultrasound. Clinical Endocrinology 2003 446450. (https://doi.org/10.1046/j.1365-2265.2003.01737.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 4

    Bakalov VK, Shawker T, Ceniceros I, Bondy CA. Uterine development in Turner syndrome. Journal of Pediatrics 2007 528.e1531.e1. (https://doi.org/10.1016/j.jpeds.2007.04.031)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 5

    Elsedfy HH, Hamza RT, Farghaly MH, Ghazy MS. Uterine development in patients with Turner syndrome: relation to hormone replacement therapy and karyotype. Journal of Pediatric Endocrinology and Metabolism 2012 441445. (https://doi.org/10.1515/jpem-2012-0040)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 6

    Snajderova M, Mardesic T, Lebl J, Gerzova H, Teslik L, Zapletalova J & Czech National Study Group for HRT Optimization in Paediatric and Adolescent Endocrinology and Gynaecology. The uterine length in women with Turner syndrome reflects the postmenarcheal daily estrogen dose. Hormone Research 2003 198204. (https://doi.org/10.1159/000073233)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 7

    Doerr HG, Bettendorf M, Hauffa BP, Mehls O, Partsch CJ, Said E, Sander S, Schwarz HP, Stahnke N, Steinkamp H, et al.Uterine size in women with Turner syndrome after induction of puberty with estrogens and long-term growth hormone therapy: results of the German IGLU Follow-up Study 2001. Human Reproduction 2005 14181421. (https://doi.org/10.1093/humrep/deh764)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 8

    Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Willer CJ, Jackson AU, Vedantam S, Raychaudhuri S, et al.Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 2010 832838. (https://doi.org/10.1038/nature09410)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 9

    Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, Ingvarsson T, Jonsdottir T, Saemundsdottir J, Center JR, Nguyen TV, et al.Multiple genetic loci for bone mineral density and fractures. New England Journal of Medicine 2008 23552365. (https://doi.org/10.1056/NEJMoa0801197)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 10

    Eriksson N, Benton GM, Do CB, Kiefer AK, Mountain JL, Hinds DA, Francke U, Tung JY. Genetic variants associated with breast size also influence breast cancer risk. BMC Medical Genetics 2012 53. (https://doi.org/10.1186/1471-2350-13-53)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 11

    Garcia-Closas M, Couch FJ, Lindstrom S, Michailidou K, Schmidt MK, Brook MN, Orr N, Rhie SK, Riboli E, Feigelson HS, et al.Genome-wide association studies identify four ER negative-specific breast cancer risk loci. Nature Genetics 2013 392.e1398.e1. (https://doi.org/10.1038/ng.2561)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 12

    Sowińska-Przepiera E, Andrysiak-Mamos E, Chełstowski K, Adler G, Friebe Z, Syrenicz A. Association between ER-α polymorphisms and bone mineral density in patients with Turner syndrome subjected to estroprogestagen treatment – a pilot study. Journal of Bone and Mineral Metabolism 2011 484492. (https://doi.org/10.1007/s00774-010-0247-3)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 13

    Bakalov VK, Bondy CA. Fracture risk and bone mineral density in Turner syndrome. Reviews in Endocrine and Metabolic Disorders 2008 145151. (https://doi.org/10.1007/s11154-008-9076-2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 14

    Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005 263265. (https://doi.org/10.1093/bioinformatics/bth457)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 15

    Wright AF. Genetic Variation: Polymorphisms and Mutations. Encyclopedia Of Life Sciences. Hoboken NJ, USA: John Wiley & Sons, 2005. (https://doi.org/10.1038/npg.els.0005005)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 16

    Braz AF, Costalonga EF, Montenegro LR, Trarbach EB, Antonini SR, Malaquias AC, Ramos ES, Mendonca BB, Arnhold IJ, Jorge AAL. The interactive effect of GHR-exon 3 and -202 A/C IGFBP3 polymorphisms on rhGH responsiveness and treatment outcomes in patients with Turner syndrome. Journal of Clinical Endocrinology and Metabolism 2012 E671E677. (https://doi.org/10.1210/jc.2011-2521)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 17

    Binder G, Baur F, Schweizer R, Ranke MB. The d3-growth hormone (GH) receptor polymorphism is associated with increased responsiveness to GH in Turner syndrome and short small-for-gestational-age children. Journal of Clinical Endocrinology and Metabolism 2006 659664. (https://doi.org/10.1210/jc.2005-1581)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 18

    Braz AF, Costalonga EF, Trarbach EB, Scalco RC, Malaquias AC, Guerra-Junior G, Antonini SR, Mendonca BB, Arnhold IJ, Jorge AA. Genetic predictors of long-term response to growth hormone (GH) therapy in children with GH deficiency and Turner syndrome: the influence of a SOCS2 polymorphism. Journal of Clinical Endocrinology and Metabolism 2014 E1808E1813. (https://doi.org/10.1210/jc.2014-1744)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 19

    Bianco B, Verreschi IT, Oliveira KC, Guedes AD, Galera BB, Galera MF, Barbosa CP, Lipay MV. PTPN22 polymorphism is related to autoimmune disease risk in patients with Turner syndrome. Scandinavian Journal of Immunology 2010 256259. (https://doi.org/10.1111/j.1365-3083.2010.02438.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 20

    Villanueva-Ortega E, Ahedo B, Fonseca-Sánchez MA, Pérez-Durán J, Garibay-Nieto N, Macías-Galavíz MT, Trujillo-Cabrera Y, García-Latorre E, Queipo G. Analysis of PTPN22, ZFAT and MYO9B polymorphisms in Turner syndrome and risk of autoimmune disease. International Journal of Immunogenetics 2017 153157. (https://doi.org/10.1111/iji.12323)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 21

    Peralta López M, Centeno V, Miras M, Silvano L, Pérez A, Muñoz L, Sobrero G, Ulla M, Tolosa de Talamoni N. Association of vitamin D receptor gene Cdx2 polymorphism with bone markers in Turner syndrome patients. Journal of Pediatric Endocrinology and Metabolism 2012 669671. (https://doi.org/10.1515/jpem-2012-0098)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 22

    Santos LOD, Bispo AVS, Barros JV, Laranjeira RSM, Pinto RDN, Silva JA, Duarte AR, Araújo J, Sandrin-Garcia P, Crovella S, et al.CTLA-4 gene polymorphisms are associated with obesity in Turner syndrome. Genetics and Molecular Biology 2018 727734. (https://doi.org/10.1590/1678-4685-GMB-2017-0312)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 23

    Massart F, Marini F, Bianchi G, Minisola S, Luisetto G, Pirazzoli A, Salvi S, Micheli D, Masi L, Brandi ML. Age-specific effects of estrogen receptors’ polymorphisms on the bone traits in healthy fertile women: the BONTURNO study. Reproductive Biology and Endocrinology 2009 32. (https://doi.org/10.1186/1477-7827-7-32)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 24

    McGuigan FE, Murray L, Gallagher A, Davey-Smith G, Neville CE, Van’t Hof R, Boreham C, Ralston SH. Genetic and environmental determinants of peak bone mass in young men and women. Journal of Bone and Mineral Research 2002 12731279. (https://doi.org/10.1359/jbmr.2002.17.7.1273)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 25

    Almeida M, Laurent MR, Dubois V, Claessens F, O'Brien CA, Bouillon R, Vanderschueren D, Manolagas SC. Estrogens and androgens in skeletal physiology and pathophysiology. Physiological Reviews 2017 135187. (https://doi.org/10.1152/physrev.00033.2015)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 26

    Herrington DM, Howard TD, Brosnihan KB, McDonnell DP, Li X, Hawkins GA, Reboussin DM, Xu J, Zheng SL, Meyers DA, et al.Common estrogen receptor polymorphism augments effects of hormone replacement therapy on E-selectin but not C-reactive protein. Circulation 2002 18791882. (https://doi.org/10.1161/01.cir.0000016173.98826.88)

    • PubMed
    • Search Google Scholar
    • Export Citation

 

  • Collapse
  • Expand
  • Figure 1

    Haploview linkage disequilibrium plot for five selected SNPs within ESR1. Haplotype blocks were constructed by using the Solid Spine method. ESR1 rs2234693 was located in a linkage disequilibrium block with rs543650 and rs9340799. ESR1 rs1038304 and rs2046211 were located in another linkage disequilibrium block nearby.

  • 1

    Gravholt CH, Andersen NH, Conway GS, Dekkers OM, Geffner ME, Klein KO, Lin AE, Mauras N, Quigley CA, Rubin K, et al.Clinical practice guidelines for the care of girls and women with Turner syndrome: proceedings from the 2016 Cincinnati International Turner Syndrome Meeting. European Journal of Endocrinology 2017 G1G70. (https://doi.org/10.1530/EJE-17-0430)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 2

    Cleemann L, Holm K, Fallentin E, Skouby SO, Smedegaard H, Møller N, Borch-Christensen H, Jeppesen EM, Wieslander SB, Andersson AM, et al.Uterus and ovaries in girls and young women with Turner syndrome evaluated by ultrasound and magnetic resonance imaging. Clinical Endocrinology 2011 756761. (https://doi.org/10.1111/j.1365-2265.2011.03995.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 3

    McDonnell CM, Coleman L, Zacharin MR. A 3-year prospective study to assess uterine growth in girls with Turner’s syndrome by pelvic ultrasound. Clinical Endocrinology 2003 446450. (https://doi.org/10.1046/j.1365-2265.2003.01737.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 4

    Bakalov VK, Shawker T, Ceniceros I, Bondy CA. Uterine development in Turner syndrome. Journal of Pediatrics 2007 528.e1531.e1. (https://doi.org/10.1016/j.jpeds.2007.04.031)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 5

    Elsedfy HH, Hamza RT, Farghaly MH, Ghazy MS. Uterine development in patients with Turner syndrome: relation to hormone replacement therapy and karyotype. Journal of Pediatric Endocrinology and Metabolism 2012 441445. (https://doi.org/10.1515/jpem-2012-0040)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 6

    Snajderova M, Mardesic T, Lebl J, Gerzova H, Teslik L, Zapletalova J & Czech National Study Group for HRT Optimization in Paediatric and Adolescent Endocrinology and Gynaecology. The uterine length in women with Turner syndrome reflects the postmenarcheal daily estrogen dose. Hormone Research 2003 198204. (https://doi.org/10.1159/000073233)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 7

    Doerr HG, Bettendorf M, Hauffa BP, Mehls O, Partsch CJ, Said E, Sander S, Schwarz HP, Stahnke N, Steinkamp H, et al.Uterine size in women with Turner syndrome after induction of puberty with estrogens and long-term growth hormone therapy: results of the German IGLU Follow-up Study 2001. Human Reproduction 2005 14181421. (https://doi.org/10.1093/humrep/deh764)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 8

    Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Willer CJ, Jackson AU, Vedantam S, Raychaudhuri S, et al.Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 2010 832838. (https://doi.org/10.1038/nature09410)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 9

    Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, Ingvarsson T, Jonsdottir T, Saemundsdottir J, Center JR, Nguyen TV, et al.Multiple genetic loci for bone mineral density and fractures. New England Journal of Medicine 2008 23552365. (https://doi.org/10.1056/NEJMoa0801197)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 10

    Eriksson N, Benton GM, Do CB, Kiefer AK, Mountain JL, Hinds DA, Francke U, Tung JY. Genetic variants associated with breast size also influence breast cancer risk. BMC Medical Genetics 2012 53. (https://doi.org/10.1186/1471-2350-13-53)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 11

    Garcia-Closas M, Couch FJ, Lindstrom S, Michailidou K, Schmidt MK, Brook MN, Orr N, Rhie SK, Riboli E, Feigelson HS, et al.Genome-wide association studies identify four ER negative-specific breast cancer risk loci. Nature Genetics 2013 392.e1398.e1. (https://doi.org/10.1038/ng.2561)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 12

    Sowińska-Przepiera E, Andrysiak-Mamos E, Chełstowski K, Adler G, Friebe Z, Syrenicz A. Association between ER-α polymorphisms and bone mineral density in patients with Turner syndrome subjected to estroprogestagen treatment – a pilot study. Journal of Bone and Mineral Metabolism 2011 484492. (https://doi.org/10.1007/s00774-010-0247-3)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 13

    Bakalov VK, Bondy CA. Fracture risk and bone mineral density in Turner syndrome. Reviews in Endocrine and Metabolic Disorders 2008 145151. (https://doi.org/10.1007/s11154-008-9076-2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 14

    Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005 263265. (https://doi.org/10.1093/bioinformatics/bth457)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 15

    Wright AF. Genetic Variation: Polymorphisms and Mutations. Encyclopedia Of Life Sciences. Hoboken NJ, USA: John Wiley & Sons, 2005. (https://doi.org/10.1038/npg.els.0005005)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 16

    Braz AF, Costalonga EF, Montenegro LR, Trarbach EB, Antonini SR, Malaquias AC, Ramos ES, Mendonca BB, Arnhold IJ, Jorge AAL. The interactive effect of GHR-exon 3 and -202 A/C IGFBP3 polymorphisms on rhGH responsiveness and treatment outcomes in patients with Turner syndrome. Journal of Clinical Endocrinology and Metabolism 2012 E671E677. (https://doi.org/10.1210/jc.2011-2521)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 17

    Binder G, Baur F, Schweizer R, Ranke MB. The d3-growth hormone (GH) receptor polymorphism is associated with increased responsiveness to GH in Turner syndrome and short small-for-gestational-age children. Journal of Clinical Endocrinology and Metabolism 2006 659664. (https://doi.org/10.1210/jc.2005-1581)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 18

    Braz AF, Costalonga EF, Trarbach EB, Scalco RC, Malaquias AC, Guerra-Junior G, Antonini SR, Mendonca BB, Arnhold IJ, Jorge AA. Genetic predictors of long-term response to growth hormone (GH) therapy in children with GH deficiency and Turner syndrome: the influence of a SOCS2 polymorphism. Journal of Clinical Endocrinology and Metabolism 2014 E1808E1813. (https://doi.org/10.1210/jc.2014-1744)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 19

    Bianco B, Verreschi IT, Oliveira KC, Guedes AD, Galera BB, Galera MF, Barbosa CP, Lipay MV. PTPN22 polymorphism is related to autoimmune disease risk in patients with Turner syndrome. Scandinavian Journal of Immunology 2010 256259. (https://doi.org/10.1111/j.1365-3083.2010.02438.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 20

    Villanueva-Ortega E, Ahedo B, Fonseca-Sánchez MA, Pérez-Durán J, Garibay-Nieto N, Macías-Galavíz MT, Trujillo-Cabrera Y, García-Latorre E, Queipo G. Analysis of PTPN22, ZFAT and MYO9B polymorphisms in Turner syndrome and risk of autoimmune disease. International Journal of Immunogenetics 2017 153157. (https://doi.org/10.1111/iji.12323)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 21

    Peralta López M, Centeno V, Miras M, Silvano L, Pérez A, Muñoz L, Sobrero G, Ulla M, Tolosa de Talamoni N. Association of vitamin D receptor gene Cdx2 polymorphism with bone markers in Turner syndrome patients. Journal of Pediatric Endocrinology and Metabolism 2012 669671. (https://doi.org/10.1515/jpem-2012-0098)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 22

    Santos LOD, Bispo AVS, Barros JV, Laranjeira RSM, Pinto RDN, Silva JA, Duarte AR, Araújo J, Sandrin-Garcia P, Crovella S, et al.CTLA-4 gene polymorphisms are associated with obesity in Turner syndrome. Genetics and Molecular Biology 2018 727734. (https://doi.org/10.1590/1678-4685-GMB-2017-0312)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 23

    Massart F, Marini F, Bianchi G, Minisola S, Luisetto G, Pirazzoli A, Salvi S, Micheli D, Masi L, Brandi ML. Age-specific effects of estrogen receptors’ polymorphisms on the bone traits in healthy fertile women: the BONTURNO study. Reproductive Biology and Endocrinology 2009 32. (https://doi.org/10.1186/1477-7827-7-32)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 24

    McGuigan FE, Murray L, Gallagher A, Davey-Smith G, Neville CE, Van’t Hof R, Boreham C, Ralston SH. Genetic and environmental determinants of peak bone mass in young men and women. Journal of Bone and Mineral Research 2002 12731279. (https://doi.org/10.1359/jbmr.2002.17.7.1273)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 25

    Almeida M, Laurent MR, Dubois V, Claessens F, O'Brien CA, Bouillon R, Vanderschueren D, Manolagas SC. Estrogens and androgens in skeletal physiology and pathophysiology. Physiological Reviews 2017 135187. (https://doi.org/10.1152/physrev.00033.2015)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 26

    Herrington DM, Howard TD, Brosnihan KB, McDonnell DP, Li X, Hawkins GA, Reboussin DM, Xu J, Zheng SL, Meyers DA, et al.Common estrogen receptor polymorphism augments effects of hormone replacement therapy on E-selectin but not C-reactive protein. Circulation 2002 18791882. (https://doi.org/10.1161/01.cir.0000016173.98826.88)

    • PubMed
    • Search Google Scholar
    • Export Citation