CC BY-NC-ND 4.0 · Rev Bras Ortop (Sao Paulo) 2020; 55(01): 008-016
DOI: 10.1016/j.rboe.2018.03.001
Artigo de Revisão
Sociedade Brasileira de Ortopedia e Traumatologia. Published by Thieme Revinter Publicações Ltda Rio de Janeiro, Brazil

Association of ESRα XbaI A > G, ESRα PvuII T > C and ESRβ AlwNI T > C Polymorphisms with the Risk of Developing Adolescent Idiopathic Scoliosis: A Systematic Review and Genetic Meta-analysis[*]

Article in several languages: português | English
Mohammad Reza Sobhan
1   Departamento de Ortopedia, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
,
Masoud Mahdinezhad-Yazdi*
1   Departamento de Ortopedia, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
,
2   Departamento de Genética Médica, Faculdade de Medicina, Shiraz University of Medical Sciences, Shiraz, Iran
,
Mohammadali Jafari
3   Departamento de Medicina de Emergência, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
,
Ali Raee-Ezzabadi
3   Departamento de Medicina de Emergência, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
,
Hossein Neamatzadeh
4   Departamento de Genética Médica, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
5   Centro de Pesquisa em Saúde Materno-Infantil, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
› Author Affiliations
Further Information

Endereço para correspondência

Seyed Alireza Dastgheib, PhD
Assistant Professor, Department of Medical Genetics, School of Medicine, Shiraz University of Medical Sciences
Shiraz
Iran   

Publication History

30 January 2018

13 March 2018

Publication Date:
09 January 2020 (online)

 

Abstract

Several association studies of genes polymorphisms on estrogen receptors-α and β with respect to adolescent idiopathic scoliosis (AIS) have been published in the past two decades. However, the association with AIS, especially among different ethnic subgroups, still remains controversial. Thus, we investigated these inconclusive data by performing a meta-analysis to systematically evaluate the association.

A literature search was conducted in the PubMed, ISI Web of Science, EMBASE, SCOPUS, EBSCO, Cochrane Library, China National Knowledge Infrastructure (CNKI) and Wanfang databases until January 20, 2018. The strength of relationship was assessed using odds ratios (ORs) and 95% confidence intervals (95%CIs).

A total of 12 case–control studies with 4,304 cases of AIS and 3,123 controls met our criteria. The pooled ORs indicated that the ESRα XbaI A > G, ESRα PvuII T > C and ESRβ AlwNI T > C polymorphisms were not significantly associated with the risk of developing AIS in the overall analysis. However, we found a significant association between the ESRα XbaI A > G polymorphism and AIS under the homozygote model (GG versus AA; OR = 1.448, 95%CI: 1.052–1.993; p = 0.023).

The present meta-analysis suggests that the ESRα XbaI A > G, ESRα PvuII T > C and ESRβ AlwNI T > C polymorphisms may not be associated with the risk of developing AIS in the overall analysis. However, ESRα XbaI A > G might have an influence on the susceptibility to develop AIS among Asians. Considering the limited sample size and ethnicity, further larger studies are needed to provide a more precise estimation of the associations.


#

Introduction

Adolescent idiopathic scoliosis (AIS) is a clinically significant disorder with high heritability that affects between 2% and 4% of the world population.[1] It is estimated that AIS affects up to 3% of all children, and its onset occurs after the age of 10.[2] The pathogenesis of AIS is poorly understood.[3] It seems that the cause of AIS is complex, and the possible etiology includes genetic factors, hormones and metabolic dysfunction, abnormal growth, and environmental and lifestyle factors.[4] Out of all these factors, the genetic factors are widely well documented. Dominant autosomal and dominant X-linked inheritances have been described for the etiology of AIS, but most AIS families have complex and non-Mendelian inheritance.[5] [6] Further efforts to understand the genetic basis of AIS have focused on genome wide association and candidate gene studies. Population studies of index patients and their families have shown that 11% of first-degree relatives are affected, as are 2.4% and 1.4% of second- and third-degree relatives respectively.[7]

Many previous studies have focused on the identification of genes involved in the onset and evolution of AIS, such as melatonin receptor 1B (MTNR1B), ladybird home-box[1] (LBX1), tryptophan hydroxylase[1] (TPH1), arylalkylamine N-acetyltransferase (AA-NAT) and basonuclin[2] (BNC2), and, among them, the LBX1 gene is widely investigated.[4] [8] Because AIS develops during puberty and it is more common among females, both forms of the estrogen receptors (ESRα and ESRβ) have been implicated as candidate genes.[3] [4]

Though the functional implications of the selected ESRα XbaI (A/G), ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms are still not fully elucidated, it has been suggested that intronic changes in the ERα sequence may modify the expression or affinity of this receptor by estrogen.[9] Therefore, we have conducted the current meta-analysis to better understand the association between the polymorphisms of the ESRα and ESRβ genes and AIS using human and animal-based studies.


#

Materials and Methods

Literature Search Strategy

The present meta-analysis was conducted and reported in accordance with the Preferred Reporting Items for Systematic reviews and Meta-analyses (PRISMA) guidelines. We performed a search in the Medline, ISI Web of Science, EMBASE, SCOPUS, EBSCO, Cochrane Library, China National Knowledge Infrastructure (CNKI), and Wanfang databases covering all articles published until January 20, 2018. These searches were conducted using the following keywords: scoliosis, AIS, ESR1, ESRα, ESR2, ESRβ, estrogen receptor, -351 A>G, XbaI, rs9340799, rs9340799, PvuII, rs2234693, AlwNI, and rs1256120, gene or allele, genotype, mutation, variant, variation, polymorphism. We evaluated potentially relevant publications by examining their titles and abstracts, and all of the studies matching the eligible criteria were retrieved. All eligible studies were examined carefully, and their references were checked for other relevant publications. If more than one article had been published by the same author using the same case series, we selected the study in which the most individuals were investigated. In addition, no restrictions were placed on language, and only published studies with full-text articles were included.


#

Inclusion and Exclusion Criteria

The studies included in the current meta-analysis had to meet all of the following criteria: a) evaluation of the ESRα XbaI (A/G), ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms and the risks of developing AIS; b) case–control studies; and c) sufficient published data to estimate an odds ratio (OR) with a 95% confidence interval (95%CI). The exclusion criteria were as follows: a) researches not related to AIS; b) population studies; c) summaries, comments, case reports, letters, and reviews; d) duplicates of previous publications; and e) studies with no sufficient data provided.


#

Data Extraction

Two investigators extracted the data independently, and the results were reviewed by a third investigator. From each study, the following items were noted: name of the first author, year of publication, country, number of cases and controls, gene polymorphisms, minor allele frequencies (MAFs), and deviation from the Hardy–Weinberg equilibrium (HWE) of the control group. If there were any disagreements, they were solved by discussion and consultation with another researcher.


#

Quality Assessment

The quality assessment of the included studies mostly agreed with the confirmation of the HWE for the genotype distribution of the ESRα XbaI (A/G), ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms in the controls. If the studies deviated from the HWE in the controls, they were defined as low-quality studies. Conversely, studies with the genotype distribution of the ESRα XbaI (A/G), ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms in the controls in accordance with the HWE (p > 0.05) were defined as high-quality studies.


#

Statistical Analysis

The strength of the association between the ESRα polymorphism and the risk of developing AIS was measured by ORs, whereas a sense of the precision of the estimate was provided by the 95%CIs. The ESRα XbaI (A/G) polymorphism and the susceptibility to develop AIS was evaluated using allelic (G versus A) and genotypic comparisons of codominant (GG versus AA and GA versus AA), dominant (GG + GA versus AA), and recessive (GG versus GA + AA) genetic models, in which the G allele was considered the risk allele. The ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms were evaluated according to the allele (C versus T), homozygote (CC versus TT), heterozygote (TC versus TT), dominant (TC/CC versus TT), and recessive (CC versus TC/TT) models respectively. The genetic model evaluated for the OR of the polymorphism was a dominant model. Cochran's Q statistic was used to formally test for heterogeneity. The heterogeneity was considered significant when p < 0.01 in the Q statistic. The percentage variability of the pooled OR attributable to heterogeneity among studies was quantified with the I2 metric, which is independent of the number of studies in the meta-analysis, and considers values between 0% and 100%, with higher values denoting a greater degree of heterogeneity 33 (I2 = 0–25%: no heterogeneity; I2 = 25–50%: moderate heterogeneity; I2 = 50–75%: great heterogeneity; I2 = 75–100%: extreme heterogeneity).[10] [11] The random-effects model (DerSimonian–Laird method) or fixed-effects model (Mantel–Haenszel method) were used to calculate pooled effect estimates in the presence or absence of heterogeneity. The HWE was assessed by Fisher exact test.[12] The sensitivity was assessed by changing the effect models. A statistic alteration in the significance indicated that the results were unstable.[13] In addition, one-way sensitivity analysis was also used to assess the stability of the results by omitting one of the studies each time. The publication bias was assessed by visual inspection of funnel plots, in which the standard error of log(OR) of each study was plotted against its log(OR). The publication bias was qualitatively assessed by making Begg funnel plots, and it was quantitatively evaluated by the Egger test. Values of p < 0.05 were considered representative of statistically significant publication bias. In addition, an asymmetric plot indicated a possible publication bias.[11] Subgroup analyses were performed according to sample size, ethnicity, source of control, family history status and genotyping method separately. All statistical analyses were performed with the Comprehensive Meta-Analysis (CMA, Biostat, Englewood, NJ, US) software, version 2.1. All p-values in the meta-analysis were two-sided, and p-values < 0.05 were considered significant.


#
#

Results

Characteristics of the Included Studies

The present study met the requirements of the PRISMA statement ([Fig. 1]). In the initial screening, we have identified 174 publications in the database searches, and 98 publications were excluded after we read the titles due to the high rate of repetition of articles in different databases. After eliminating the duplicates, 64 articles were excluded after we screened the titles, abstracts and full texts because they were inconsistent with our inclusion criteria. Finally, 7 eligible case–control studies with a total of 2,377 cases and 1,770 controls were included in the meta-analysis.[14] [15] [16] [17] [18] [19] [20] The characteristics of the included studies are summarized in [Table 1]. The studies were published between 2006 and 2014, and were performed in China, Japan and Poland. The distribution of the ethnic groups among these seven studies was as follows: in 2 studies with 535 cases and 425 controls, the subjects were Caucasian,[18] [20] and in 5 studies with 1,842 cases and 1,982 controls, the subjects were Asian.[14] [15] [16] [17] [18] [19] The genotype distributions among the controls of all studies followed the HWE, except for two studies on XbaI (rs9340799)[14] [16], and one study on PvuII (rs2234693).[14] Therefore, according to the quality criteria, there were five high-quality studies and two low-quality studies.

Zoom Image
Fig. 1 Flow diagram of the study selection process.
Table 1

First author

Country (ethnicity)

Case/Control

Gender

Cases

Control

MAFs

HWE

Genotypes

Alleles

Genotypes

Alleles

XbaI (rs9340799)

AA

GA

GG

A

G

AA

AG

GG

A

G

Wu et al.[14] (2006)

China (Asian)

202/174

Overall

72

76

54

220

184

82

66

26

230

118

0.339

0.042

Female

64

70

51

198

172

72

59

23

223

105

0.340

Male

8

6

3

22

12

10

7

3

7

13

0.325

Tang et al[15] (2006)

China (Asian)

540/260

Female

328

176

36

832

248

157

85

18

399

121

0.232

0.173

Zhao et al[16] (2009)

China (Asian)

100/100

Female

58

34

8

150

50

55

31

14

141

59

0.295

0.010

Takahashi et al[17] (2011)

Japan (Asian)

798/637

Female

526

248

24

1300

296

421

196

20

1038

236

0.185

0.645

Janusz et al[18] (2013)

Poland (Caucasian)

287/182

Female

96

141

50

333

241

61

92

29

214

150 0.412

0.559

PvuII (rs2234693)

TT

TC

CC

T

C

TT

TC

CC

T

C

Wu et al[14] (2006)

China (Asian)

202/174

Overall

71

92

39

234

170

64

70

40

198

150

0.431

0.017

Female

65

84

36

214

156

57

61

36

175

133

0.431

Male

6

8

3

20

14

7

9

4

23

17

0.425

Tang et al[15] (2006)

China (Asian)

540/260

Female

201

249

93

648

432

102

128

30

332

188

0.361

0.284

Zhao et al[16] (2009)

China (Asian)

100/100

Female

31

51

18

113 87

31

50

19

112 88

0.44

0.883

Janusz et al[18] (2013)

Poland (Caucasian)

287/182

Female

77

144

66

298

276

42

93

47

177

187

0.513

0.758

AlwNI (rs1256120)

CC

TC

TT

C

T

CC

TC

TT

C

T

Zhang et al[19] (2009)

China (Asian)

202/174

Overall

36

105

77

177

259

10

55

75

75

205

0.732

0.984

Female

30

85

61

145

207

5

30

45

40

120

0.750

Male

6

20

16

32

52

5

25

30

35

85

0.7083

Takahashi et al[17] (2010)

Japan (Asian)

798/637

Female

99

368

331

566

1030

79

306

252

464

810

0.6358

0.347

Kotwicki et al[20] (2014)

Poland (Caucasian)

248/243

Female

164

74

10

402

94

159

77

7

395

91

0.1872

0.521

ESRα XbaI A > G polymorphism

In total, 5 case–control studies with 1,927 cases and 1,353 controls investigated the association of ESRα XbaI A > G with the risk of developing AIS. Overall, the present meta-analysis suggested that there was no significant association between the ESRα XbaI A > G polymorphism and the risk of developing AIS (G versus A, OR = 1.071, 95%CI: 0.879–1.304, p = 0.497 ([Fig. 2A]); GA versus AA, OR = 1.037, 95%CI: 0.889–1.209, p = 0.641; GG versus AA, OR = 1.292, 95%CI: 0.979–1.705, p = 0.070; GG + GA versus AA, OR = 1.052, 95%CI: 0.910–1.217, p = 0.492; GG versus GA + AA, OR = 1.071, 95%CI: 0.922–1.243, p = 0.369). In the subgroup analyses by ethnicity, there was a significant association between ESRα XbaI A > G and the risk of developing AIS among Asians under the homozygote model (GG versus AA, OR = 1.448, 95%CI: 1.052–1.993, p = 0.023).

ESRα PvuII T > C polymorphism

In total, 4 case–control studies with 1,129 cases and 714 controls investigated the association of the ESRα PvuII T > C polymorphism with the risk of developing AIS. The pooled results based on all included studies did not show a significant association between the ESRα PvuII T > C and the risk of developing AIS under the allele (C versus T, OR = 1.018, 95%CI: 0.888–1.166, p = 0.800), the heterozygote (CT versus TT, OR = 0.996, 95%CI: 0.804–1.234, p = 0.973), the homozygote (CC versus TT, OR = 1.045, 95%CI: 0.789–1.383, p = 0.760) ([Fig. 2B]), the dominant (CC + CT versus TT, OR = 1.005, 95%CI: 0.821–1.229, p = 0.965), and the recessive models (CC versus CT + TT, OR = 1.034, 95%CI: 0.808–1.322, p = 0.792). In the subgroup analyses by ethnicity, there was no significant association between the ESRα PvuII T > C polymorphism and the risk of developing AIS under all five genetic models ([Table 2]).

Table 2

Polymorphism

Genetic model

Type of model

Heterogeneity (H)

Odds Ratio (OR)

Publication bias

I2 (%)

pH

OR

95%CI

Z test

p OR

p Begg

P Egger

XbaI (rs9340799)

Overall

G versus A

Random

61.79

0.033

1.071

0.879–1.304

0.679

0.497

0.806

0.949

GA versus AA

Fixed

0.00

0.869

1.037

0.889–1.209

0.466

0.641

0.806

0.401

GG versus AA

Fixed

46.01

0.116

1.292

0.979–1.705

1.809

0.070

0.462

0.721

GG + GA versus AA

Fixed

16.61

0.309

1.052

0.910–1.217

0.686

0.492

0.806

0.705

GG versus GA + AA

Fixed

42.61

0.137

1.071

0.922–1.243

0.897

0.369

0.086

0.233

By ethnicity

Asian

G versus A

Random

71.18

0.015

1.079

0.834–1.396

0.580

0.562

1.000

0.950

GA versus AA

Fixed

0.00

0.765

1.048

0.888–1.236

0.552

0.581

0.308

0.359

GG versus AA

Fixed

44.82

0.142

1.448

1.052–1.993

2.268

0.023

1.000

0.915

GG + GA versus AA

Fixed

36.58

0.193

1.060

0.906–1.241

0.733

0.464

1.000

0.709

GG versus GA + AA

Fixed

56.02

0.078

1.084

0.922–1.273

0.976

0.329

0.308

0.268

PvuII (rs2234693)

Overall

C versus T

Fixed

5.54

0.365

1.018

0.888–1.166

0.253

0.800

1.000

0.472

CT versus TT

Fixed

0.00

0.788

0.996

0.804–1.234

−0.033

0.973

0.734

0.903

CC versus TT

Fixed

34.62

0.204

1.045

0.789–1.383

0.305

0.760

1.000

0.608

CC + CT versus TT

Fixed

0.00

0.548

1.005

0.821–1.229

0.044

0.965

0.308

0.574

CC versus CT + TT

Fixed

45.60

0.138

1.034

0.808–1.322

0.263

0.792

1.000

0.767

By ethnicity

Asian (Chinese)

C versus T

Fixed

0.00

0.476

1.075

0.917–1.260

0.889

0.374

1.000

0.398

CT versus TT

Fixed

0.00

0.813

1.044

0.819–1.330

0.345

0.730

1.000

0.714

CC versus TT

Fixed

27.64

0.251

1.178

0.847–1.639

0.973

0.330

1.000

0.578

CC + CT versus TT

Fixed

0.00

0.949

1.085

0.864–1.364

0.704

0.481

0.296

0.041

CC versus CT + TT

Fixed

55.04

0.108

1.131

0.838–1.526

0.805

0.421

1.000

0.660

AlwNI (rs1256120)

Overall

T versus C

Fixed

0.00

0.538

1.072

0.946–1.215

1.090

0.276

1.000

0.686

TC versus CC

Fixed

0.00

0.373

0.896

0.705–1.139

−0.896

0.370

0.296

0.078

TT versus CC

Random

90.86

≤0.001

0.950

0.266–3.386

−0.080

0.937

1.000

0.489

TT + CT versus CC

Random

84.41

0.002

1.374

0.758–2.493

1.047

0.295

0.296

0.212

TT versus TC + CC

Fixed

0.00

0.391

1.032

0.857–1.242

0.333

0.739

1.000

0.970


#

ESRβ AlwNI C > T Polymorphism

A total of 3 case–control studies with 1,248 cases and 1,054 controls were selected to estimate the association of the ESRβ AlwNI C > T polymorphism and the susceptibility to develop AIS. As a result, no statistically no significant association was found in any of the genetic models (C versus T, OR = 1.072, 95%CI: 0.946–1.215, p = 0.276; CT versus TT, OR = 0.896, 95%CI: 0.705–1.139, p = 0.370; CC versus TT, OR = 0.950, 95%CI: 0.266–3.386, p = 0.937; CC + CT versus TT, OR = 1.374, 95%CI: 0.758–2.493, p = 0.295 ([Fig. 2C]); CC versus CT + TT, OR = 1.032, 95%CI: 0.857–1.242, p = 0.739); ([Fig. 2]).

Zoom Image
Fig. 2 Forest plots of the association between the ESRα XbaI A > G, ESRα PvuII T > C, and ESRβ AlwNI C > T polymorphisms under the dominant model and the risk of developing adolescent idiopathic scoliosis (AIS). (A) XbaI A > G (allele model); (B) PvuII T > C (homozygote model); and (C) AlwNI C > T (dominant model). The squares and horizontal lines correspond to the study-specific OR and 95%CI. The area of the squares reflects the study-specific weight (inverse of the variance). The diamond represents the summary OR and 95%CI.

#

Heterogeneity and Sensitivity Analyses

The heterogeneity was not significant for the ESRα XbaI A > G, ESRα PvuII T > C and ESRβ AlwNI T > C polymorphisms in most genetic models, which suggested that the polymorphisms were not the source of heterogeneity among the studies ([Table 2]). The sensitivity analyses by sequential omission of any individual study, one at a time, or by omitting studies in which the genotype distributions in the healthy controls significantly deviated from the HWE, did not materially alter the pooled ORs, indicating that the results were stable.


#

Publication Bias

The Begg funnel plot and the Egger tests were performed to estimate the publication bias of the studies regarding the association between the susceptibility to develop AIS and the ESRα XbaI A > G, ESRα PvuII T > C and ESRβ AlwNI T > C polymorphisms. The shape of the funnel plot did not reveal any evidence of obvious asymmetry for the polymorphisms under any of the genetic models. Furthermore, the p-values of the Egger tests were higher than 0.05, providing statistical evidence of the symmetry of the funnel plots. However, the results of the Egger tests showed evidence of publication bias for the ESRα PvuII T > C polymorphism in Asians under the dominant model (CC + CT versus TT: p Begg = 0.296, p Egger = 0.041); ([Fig. 3]).

Zoom Image
Fig. 3 Funnel plot to detect publication bias for association of ESRα XbaIA > G and ESRα PvuII T > C polymorphisms with risk of adolescent idiopathic scoliosis (AIS). A: XbaIA > G (allele model: G versus A); B: PvuII T > C (recessive model: CC versus CT TT).

#
#

Discussion

The pathogenesis of AIS is a complex process. It is known that genetic factors play an important role in the susceptibility to develop AIS.[14] [15] [16] However, most of the molecular mechanisms leading to AIS development are still unknown. Gene mutations in various loci have been identified by genetic studies, and a genetic basis for the AIS pathogenesis has been established.[17] Although many epidemiological studies have been conducted to assess the roles of the ESRα and ESRβ polymorphisms and the risk of developing AIS in different populations, the results have been inconclusive.[18] [20] The ESRα and ESRβ gene polymorphisms have been linked to a higher risk of developing different conditions.[21] The human ERα gene is located on chromosome 6q25, it extends for more than 140 kb, and includes 8 exons.[22] The most studied polymorphisms in this gene are the PvuII T > C and XbaI A > G in intron 1, 397 bp and 351 bp upstream of exon 2 respectively.[23] The gene that encodes ESRβ is located on chromosome 14q23,[1] and the potential association of the single-nucleotide polymorphisms (SNPs) in ESRβ (RsaI G > A, AluI A > G and AlwNI C > T) and disease have not been studied previously.[24]

The current updated meta-analysis, including 12 eligible studies with a total of 4,304 AIS cases and 3,123 controls, provides a comprehensive examination of the evidence currently available on the associations of the ESRα and ESRβ polymorphisms with the susceptibility to develop AIS. The results suggested that the polymorphisms in question had no significant association with the risk of developing AIS. The present meta-analysis is consistent with previous meta-analyses by Chen et al[25] (ESRα Xbal), Yang et al[26] (ESRα Xbal and ERα PvuII) and Cao et al[4] (ESRβ AlwNI), who found negative results for the association between the ESRα and ESRβ polymorphisms and the risk of developing AIS. However, Inoue et al[27] and Wu et al[14] reported a significant association between the ESRα XbaI A > G and ESRα PvuII T > C polymorphisms and the risk of developing AIS. In 2016, Cao et al[4] reported that there was no significant association between the ESRα XbaI A > G and ESRα PvuII T > C polymorphisms and the risk of developing AIS in a total population analysis. However, for the comparison with the study by Cao et al,[4] the subgroup analysis by ethnicity was also performed. In the present meta-analysis, we found that the ESRα XbaI A > G polymorphism was associated with the risk of developing AIS in Asians under the homozygote model (GG versus AA, OR = 1.448, 95%CI: 1.052–1.993, p = 0.023). The discrepancy between ethnicity subgroups may be due to the limited studies in Europe, since only one study was conducted among subjects of Caucasian ethnicity. Therefore, we performed a stratified analysis only on the Chinese population (Asians). Our data revealed that the XbaI G allele was a risk factor in the Chinese population. Unlike the ESRα XbaI polymorphism, no significant difference was found in the ESRα PvuII T > C polymorphism allele or genotype distribution among AIS patients and healthy individuals.

In the current meta-analysis, a fixed-effects or a random-effects models were used based on the heterogeneity testing. The differences in the studied populations with different genetic backgrounds and variations in sample selection and environmental exposures may result in heterogeneity.[28] [29] [30] Our meta-regression analysis also showed that ethnicity in case groups and control groups significantly contributed to the heterogeneity. When limiting the analysis to the studies within the HWE, the result was not altered, suggesting that the present meta-analysis is robust and credible. Moreover, we performed a sensitivity analysis according to sample size and leave-one-out cross validation to determine whether the change in the inclusion criteria by removing one study each time did not materially affected the original results.

There are some limitations to the present meta-analysis. First, the sample size was comparatively small and had insufficient statistical power for the association to be estimated. Second, we have included only studies published in English and Chinese; therefore, publication bias may have occurred. Third, the greatest proportion of statistical power was provided by Asians. There were not enough studies in Caucasians, which limited the statistical power. No study from other parts of the world was found, such as Africa and Latin America. This suggests a partial result that is only relevant to the Asian and Caucasian subgroups. Fourth, in interpreting the results, we should mention that, as with other complex traits, the risk of developing AIS may be modulated by several other genetic markers and candidate genes besides the ESRα and ESRβ genes. Therefore, further investigations on the haplotypic effect of the polymorphisms and the study of multiple polymorphisms in different genes are necessary. Finally, due to the unavailability of other detailed information, our results were based on single-factor estimates without adjustments for other risk factors. Further evaluations of the risk of developing AIS should pay more attention to the potential gene–gene, gene–environment interactions, and even the interactions of different polymorphism of the ESRα and ESRβ genes and other loci.


#

Final Consideration

In summary, the present meta-analysis suggests that the ESRα XbaI (A/G), ESRα PvuII (T/C) and ESRβ AlwNI (C/T) polymorphisms are not associated with an increased risk of developing AIS. However, ESRa XbaI A > G might have an influence on the susceptibility to develop AIS among Asians. Based on the aforementioned limitations, it is critical that large, well-designed studies are performed to re-evaluate the potential associations of the ESRα and ESRβ gene polymorphisms with other candidate gene polymorphisms and the risk of developing AIS.


#
#

Conflito de Interesses

Os autores declaram não haver conflito de interesses.

* Originally Published by Elsevier Editora Ltda.


  • Referências

  • 1 Chettier R, Nelson L, Ogilvie JW, Albertsen HM, Ward K. Haplotypes at LBX1 have distinct inheritance patterns with opposite effects in adolescent idiopathic scoliosis. PLoS One 2015; 10 (02) e0117708
  • 2 Lenke L. Idiopathic scoliosis. Philadelphia: Lippincott, Williams & Wilkins; 2004
  • 3 Buchan JG, Alvarado DM, Haller GE. , et al. Rare variants in FBN1 and FBN2 are associated with severe adolescent idiopathic scoliosis. Hum Mol Genet 2014; 23 (19) 5271-5282
  • 4 Cao Y, Min J, Zhang Q, Li H, Li H. Associations of LBX1 gene and adolescent idiopathic scoliosis susceptibility: a meta-analysis based on 34,626 subjects. BMC Musculoskelet Disord 2016; 17: 309
  • 5 Gurnett CA, Alaee F, Bowcock A. , et al. Genetic linkage localizes an adolescent idiopathic scoliosis and pectus excavatum gene to chromosome 18 q. Spine 2009; 34 (02) E94-E100
  • 6 Wise CA, Gao X, Shoemaker S, Gordon D, Herring JA. Understanding genetic factors in idiopathic scoliosis, a complex disease of childhood. Curr Genomics 2008; 9 (01) 51-59
  • 7 Cheung KM, Wang T, Qiu GX, Luk KD. Recent advances in the aetiology of adolescent idiopathic scoliosis. Int Orthop 2008; 32 (06) 729-734
  • 8 Villa MJ. Genes associated with adolescent idiopathic scoliosis: a review. Hered Genet. 2015; 4: 146
  • 9 Meng J, Mu X, Wang YM. Influence of the XbaI polymorphism in the estrogen receptor-α gene on human spermatogenic defects. Genet Mol Res 2013; 12 (02) 1808-1815
  • 10 Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003; 327 (7414): 557-560
  • 11 Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med 2002; 21 (11) 1539-1558
  • 12 Khoram-Abadi KM, Forat-Yazdi M, Kheirandish S. , et al. DNMT3B -149 C>T and -579 G>T polymorphisms and risk of gastric and colorectal cancer: a meta-analysis. Asian Pac J Cancer Prev 2016; 17 (06) 3015-3020
  • 13 Harbord RM, Egger M, Sterne JA. A modified test for small-study effects in meta-analyses of controlled trials with binary endpoints. Stat Med 2006; 25 (20) 3443-3457
  • 14 Wu J, Qiu Y, Zhang L, Sun Q, Qiu X, He Y. Association of estrogen receptor gene polymorphisms with susceptibility to adolescent idiopathic scoliosis. Spine 2006; 31 (10) 1131-1136
  • 15 Tang NL, Yeung HY, Lee KM. , et al. A relook into the association of the estrogen receptor [alpha] gene (PvuII, XbaI) and adolescent idiopathic scoliosis: a study of 540 Chinese cases. Spine 2006; 31 (21) 2463-2468
  • 16 Zhao D, Qiu GX, Wang YP, Zhang JG, Shen JX, Wu ZH. Association between adolescent idiopathic scoliosis with double curve and polymorphisms of calmodulin1 gene/estrogen receptor-α gene. Orthop Surg 2009; 1 (03) 222-230
  • 17 Takahashi Y, Matsumoto M, Karasugi T. , et al. Replication study of the association between adolescent idiopathic scoliosis and two estrogen receptor genes. J Orthop Res 2011; 29 (06) 834-837
  • 18 Janusz P, Kotwicki T, Andrusiewicz M, Kotwicka M. XbaI and PvuII polymorphisms of estrogen receptor 1 gene in females with idiopathic scoliosis: no association with occurrence or clinical form. PLoS One 2013; 8 (10) e76806
  • 19 Zhang HQ, Lu SJ, Tang MX. , et al. Association of estrogen receptor beta gene polymorphisms with susceptibility to adolescent idiopathic scoliosis. Spine 2009; 34 (08) 760-764
  • 20 Kotwicki T, Janusz P, Andrusiewicz M, Chmielewska M, Kotwicka M. Estrogen receptor 2 gene polymorphism in idiopathic scoliosis. Spine 2014; 39 (26) E1599-E1607
  • 21 Ban S, Sata F, Kurahashi N. , et al. Genetic polymorphisms of ESR1 and ESR2 that may influence estrogen activity and the risk of hypospadias. Hum Reprod 2008; 23 (06) 1466-1471
  • 22 Kos M, Reid G, Denger S, Gannon F. Minireview: genomic organization of the human ERalpha gene promoter region. Mol Endocrinol 2001; 15 (12) 2057-2063
  • 23 Yazdi MM, Jamalaldini MH, Sobhan MR. , et al. Association of ESRa Gene Pvu II T>C XbaI A>G and BtgI G>A polymorphisms with knee osteoarthritis susceptibility: a systematic review and meta-analysis based on 22 case–control studies. Arch Bone Jt Surg 2017; 5 (06) 351-362
  • 24 Fu C, Dong WQ, Wang A, Qiu G. The influence of ESR1 rs9340799 and ESR2 rs1256049 polymorphisms on prostate cancer risk. Tumour Biol 2014; 35 (08) 8319-8328
  • 25 Chen S, Zhao L, Roffey DM, Phan P, Wai EK. Association between the ESR1 -351A>G single nucleotide polymorphism (rs9340799) and adolescent idiopathic scoliosis: a systematic review and meta-analysis. Eur Spine J 2014; 23 (12) 2586-2593
  • 26 Yang M, Li C, Li M. The estrogen receptor α gene (XbaI, PvuII) polymorphisms and susceptibility to idiopathic scoliosis: a meta-analysis. J Orthop Sci 2014; 19 (05) 713-721
  • 27 Inoue M, Minami S, Nakata Y. , et al. Association between estrogen receptor gene polymorphisms and curve severity of idiopathic scoliosis. Spine 2002; 27 (21) 2357-2362
  • 28 Sobhan MR, Mehdinejad M, Jamaladini MH, Mazaheri M, Zare-Shehneh M, Neamatzadeh H. Association between aspartic acid repeat polymorphism of the asporin gene and risk of knee osteoarthritis: A systematic review and meta-analysis. Acta Orthop Traumatol Turc 2017; 51 (05) 409-415
  • 29 Mehdinejad M, Sobhan MR, Mazaheri M, Zare Shehneh M, Neamatzadeh H, Kalantar SM. Genetic association between ERCC2, NBN RAD51 gene variants and osteosarcoma risk: a systematic review and meta-analysis. Asian Pac J Cancer Prev 2017; 18 (05) 1315-1321
  • 30 Sobhan MR, Forat Yazdi M, Mazaheri M, Zare Shehneh M, Neamatzadeh H. Association between the DNA repair gene XRCC3 rs861539 polymorphism and risk of osteosarcoma: a systematic review and meta-analysis. Asian Pac J Cancer Prev 2017; 18 (02) 549-555

Endereço para correspondência

Seyed Alireza Dastgheib, PhD
Assistant Professor, Department of Medical Genetics, School of Medicine, Shiraz University of Medical Sciences
Shiraz
Iran   

  • Referências

  • 1 Chettier R, Nelson L, Ogilvie JW, Albertsen HM, Ward K. Haplotypes at LBX1 have distinct inheritance patterns with opposite effects in adolescent idiopathic scoliosis. PLoS One 2015; 10 (02) e0117708
  • 2 Lenke L. Idiopathic scoliosis. Philadelphia: Lippincott, Williams & Wilkins; 2004
  • 3 Buchan JG, Alvarado DM, Haller GE. , et al. Rare variants in FBN1 and FBN2 are associated with severe adolescent idiopathic scoliosis. Hum Mol Genet 2014; 23 (19) 5271-5282
  • 4 Cao Y, Min J, Zhang Q, Li H, Li H. Associations of LBX1 gene and adolescent idiopathic scoliosis susceptibility: a meta-analysis based on 34,626 subjects. BMC Musculoskelet Disord 2016; 17: 309
  • 5 Gurnett CA, Alaee F, Bowcock A. , et al. Genetic linkage localizes an adolescent idiopathic scoliosis and pectus excavatum gene to chromosome 18 q. Spine 2009; 34 (02) E94-E100
  • 6 Wise CA, Gao X, Shoemaker S, Gordon D, Herring JA. Understanding genetic factors in idiopathic scoliosis, a complex disease of childhood. Curr Genomics 2008; 9 (01) 51-59
  • 7 Cheung KM, Wang T, Qiu GX, Luk KD. Recent advances in the aetiology of adolescent idiopathic scoliosis. Int Orthop 2008; 32 (06) 729-734
  • 8 Villa MJ. Genes associated with adolescent idiopathic scoliosis: a review. Hered Genet. 2015; 4: 146
  • 9 Meng J, Mu X, Wang YM. Influence of the XbaI polymorphism in the estrogen receptor-α gene on human spermatogenic defects. Genet Mol Res 2013; 12 (02) 1808-1815
  • 10 Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003; 327 (7414): 557-560
  • 11 Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med 2002; 21 (11) 1539-1558
  • 12 Khoram-Abadi KM, Forat-Yazdi M, Kheirandish S. , et al. DNMT3B -149 C>T and -579 G>T polymorphisms and risk of gastric and colorectal cancer: a meta-analysis. Asian Pac J Cancer Prev 2016; 17 (06) 3015-3020
  • 13 Harbord RM, Egger M, Sterne JA. A modified test for small-study effects in meta-analyses of controlled trials with binary endpoints. Stat Med 2006; 25 (20) 3443-3457
  • 14 Wu J, Qiu Y, Zhang L, Sun Q, Qiu X, He Y. Association of estrogen receptor gene polymorphisms with susceptibility to adolescent idiopathic scoliosis. Spine 2006; 31 (10) 1131-1136
  • 15 Tang NL, Yeung HY, Lee KM. , et al. A relook into the association of the estrogen receptor [alpha] gene (PvuII, XbaI) and adolescent idiopathic scoliosis: a study of 540 Chinese cases. Spine 2006; 31 (21) 2463-2468
  • 16 Zhao D, Qiu GX, Wang YP, Zhang JG, Shen JX, Wu ZH. Association between adolescent idiopathic scoliosis with double curve and polymorphisms of calmodulin1 gene/estrogen receptor-α gene. Orthop Surg 2009; 1 (03) 222-230
  • 17 Takahashi Y, Matsumoto M, Karasugi T. , et al. Replication study of the association between adolescent idiopathic scoliosis and two estrogen receptor genes. J Orthop Res 2011; 29 (06) 834-837
  • 18 Janusz P, Kotwicki T, Andrusiewicz M, Kotwicka M. XbaI and PvuII polymorphisms of estrogen receptor 1 gene in females with idiopathic scoliosis: no association with occurrence or clinical form. PLoS One 2013; 8 (10) e76806
  • 19 Zhang HQ, Lu SJ, Tang MX. , et al. Association of estrogen receptor beta gene polymorphisms with susceptibility to adolescent idiopathic scoliosis. Spine 2009; 34 (08) 760-764
  • 20 Kotwicki T, Janusz P, Andrusiewicz M, Chmielewska M, Kotwicka M. Estrogen receptor 2 gene polymorphism in idiopathic scoliosis. Spine 2014; 39 (26) E1599-E1607
  • 21 Ban S, Sata F, Kurahashi N. , et al. Genetic polymorphisms of ESR1 and ESR2 that may influence estrogen activity and the risk of hypospadias. Hum Reprod 2008; 23 (06) 1466-1471
  • 22 Kos M, Reid G, Denger S, Gannon F. Minireview: genomic organization of the human ERalpha gene promoter region. Mol Endocrinol 2001; 15 (12) 2057-2063
  • 23 Yazdi MM, Jamalaldini MH, Sobhan MR. , et al. Association of ESRa Gene Pvu II T>C XbaI A>G and BtgI G>A polymorphisms with knee osteoarthritis susceptibility: a systematic review and meta-analysis based on 22 case–control studies. Arch Bone Jt Surg 2017; 5 (06) 351-362
  • 24 Fu C, Dong WQ, Wang A, Qiu G. The influence of ESR1 rs9340799 and ESR2 rs1256049 polymorphisms on prostate cancer risk. Tumour Biol 2014; 35 (08) 8319-8328
  • 25 Chen S, Zhao L, Roffey DM, Phan P, Wai EK. Association between the ESR1 -351A>G single nucleotide polymorphism (rs9340799) and adolescent idiopathic scoliosis: a systematic review and meta-analysis. Eur Spine J 2014; 23 (12) 2586-2593
  • 26 Yang M, Li C, Li M. The estrogen receptor α gene (XbaI, PvuII) polymorphisms and susceptibility to idiopathic scoliosis: a meta-analysis. J Orthop Sci 2014; 19 (05) 713-721
  • 27 Inoue M, Minami S, Nakata Y. , et al. Association between estrogen receptor gene polymorphisms and curve severity of idiopathic scoliosis. Spine 2002; 27 (21) 2357-2362
  • 28 Sobhan MR, Mehdinejad M, Jamaladini MH, Mazaheri M, Zare-Shehneh M, Neamatzadeh H. Association between aspartic acid repeat polymorphism of the asporin gene and risk of knee osteoarthritis: A systematic review and meta-analysis. Acta Orthop Traumatol Turc 2017; 51 (05) 409-415
  • 29 Mehdinejad M, Sobhan MR, Mazaheri M, Zare Shehneh M, Neamatzadeh H, Kalantar SM. Genetic association between ERCC2, NBN RAD51 gene variants and osteosarcoma risk: a systematic review and meta-analysis. Asian Pac J Cancer Prev 2017; 18 (05) 1315-1321
  • 30 Sobhan MR, Forat Yazdi M, Mazaheri M, Zare Shehneh M, Neamatzadeh H. Association between the DNA repair gene XRCC3 rs861539 polymorphism and risk of osteosarcoma: a systematic review and meta-analysis. Asian Pac J Cancer Prev 2017; 18 (02) 549-555

Zoom Image
Fig. 1 Fluxograma do processo de escolha de estudos.
Zoom Image
Fig. 2 Gráficos em floresta da associação entre polimorfismos em ESRα XbaI A > G, ESRα PvuII T > C e ESRβ AlwNI C > T sob modelo dominante e o risco de escoliose idiopática da adolescênci (EIA). (A) XbaI A > G (modelo alélico); (B) PvuII T > C (modelo homozigótico); e (C) AlwNI C > T (modelo dominante). Os quadrados e linhas horizontais correspondem à razão de probabilidades (RP) e ao intervalo de confiança de 95% (IC95%) específicos do estudo. A área dos quadrados reflete o peso específico do estudo (inverso da variância). O diamante representa o resumo da RP e do IC95%.
Zoom Image
Fig. 3 Análise de gráfico de funil para detecção do viés de publicação para associação de polimorfismos em ESRα XbaI A > G e ESRα PvuII T > C com risco de escoliose idiopática do adolescente (EIA). A: XbaIA > G (modelo de alelo: G versus A); B: PvuII T > C (modelo recessivo: CC versus CT TT).
Zoom Image
Fig. 1 Flow diagram of the study selection process.
Zoom Image
Fig. 2 Forest plots of the association between the ESRα XbaI A > G, ESRα PvuII T > C, and ESRβ AlwNI C > T polymorphisms under the dominant model and the risk of developing adolescent idiopathic scoliosis (AIS). (A) XbaI A > G (allele model); (B) PvuII T > C (homozygote model); and (C) AlwNI C > T (dominant model). The squares and horizontal lines correspond to the study-specific OR and 95%CI. The area of the squares reflects the study-specific weight (inverse of the variance). The diamond represents the summary OR and 95%CI.
Zoom Image
Fig. 3 Funnel plot to detect publication bias for association of ESRα XbaIA > G and ESRα PvuII T > C polymorphisms with risk of adolescent idiopathic scoliosis (AIS). A: XbaIA > G (allele model: G versus A); B: PvuII T > C (recessive model: CC versus CT TT).