- Case control study
- Open Access
Lack of association between estrogen receptor β dinucleotide repeat polymorphism and autoimmune thyroid diseases in Japanese patients
BMC Medical Geneticsvolume 2, Article number: 1 (2001)
The autoimmune thyroid diseases (AITDs), such as Graves' disease (GD) and Hashimoto's thyroiditis (HT), appear to develop as a result of complex interactions between predisposing genes and environmental triggers. Susceptibility to AITDs is conferred by genes in the human leukocyte antigen (HLA) and genes unlinked to HLA, including the CTLA-4 gene. Recently, estrogen receptor (ER) β, located at human chromosome 14q23-24.1, was identifed. We analyzed a dinucleotide (CA)n repeat polymorphism located in the flanking region of ERβ gene in patients with AITDs and in normal subjects. High heterozygosity makes this polymorphism a useful marker in the genetic study of disorders affecting female endocrine systems. We also correlated a ERβ gene microsatellite polymorphism with bone mineral density (BMD) in the distal radius and biochemical markers of bone turnover in patients with GD in remission.
Fourteen different alleles were found in 133 patients with GD, 114 patients with HT, and 179 controls subjects. The various alleles were designated as allele*1 through allele*14 according to the number of the repeats, from 18 to 30. There was no significant difference in the distributions of ERβ alleles between patient groups and controls. Although recent study demonstrated a significant relation between a allele*9 in the ERβ gene and BMD in postmenopausal Japanese women, there were no statistically significant interaction between this allele and BMD in the distal radius, nor biochemical markers in patients with GD in remission.
The present results do not support an association between the ERβ microsatellite marker and AITD in the Japanese population. We also suggest that the ERβ microsatellite polymorphism has at most a minor pathogenic importance in predicting the risk of osteoporosis as a complication of GD.
As outlined in a paper we have already published, we analyzed a dinucleotide (TA)n repeat polymorphism lying upstream of the human estrogen receptor (ER) α gene in patients with autoimmune thyroid diseases (AITDs) in normal subjects . Seventeen different alleles were found in 130 patients with Graves' disease (GD), 93 patients with Hashimoto's thyroiditis (HT), and 190 control subjects . There was no significant difference in the distributions of ERα alleles between patients and controls .
The pathogenesis of AITDs involves complex interactions between genetic and environmental factors [2,3]. Susceptibility to AITDs is conferred by genes in the human leukocyte antigen (HLA) [4,5,6,7,8] and genes unlinked to HLA, including the cytotoxic T lymphocyte antigen (CTLA)-4 gene . The existence of an ERα gene polymorphism has been documented, and its association to some variant ERα genotypes found in breast cancer [10,11], hypertension , osteoporosis [13,14], generalized osteoarthritis , and some autoimmune diseases such as rheumatoid arthritis  has been reported.
A gene for a second type of human ER (ERβ) was recently identified . The ERβ is located at human chromosome 14q23-24.1, where the aberration of DNA copy number in bone disorders is frequently found . This locus was close to, but distinct from, a new susceptibility gene (GD-1) on chromosome 14q31 . Recently, an association between some ERβ genotypes and osteoporosis has been reported . ER belongs to the nuclear hormone receptor superfamily and modulates the transcription of target genes in response to estrogen , a potent immunomodulatory hormone . Estrogens appear to play a central role in the immune response and immune-mediated diseases . In view of the possible role of estrogens in the pathogenesis of AITDs, we analyzed a dinucleotide (CA)n repeat polymorphism located in the flanking region of ERβ gene  in patients with AITDs and in normal subjects. We also studied associaitons of a ERβ gene microsatellite polymorphism with distal radius bone mineral density (BMD) and biochemical markers of bone turnover to determine how this polymorphism might influence the development of osteoporosis as a complication of GD.
Results and discussion
Allele frequencies of ERβ gene microsatellite polymorphism in patients and controls
In general, estrogen action is thought to be mediated by ER, a member of the nuclear receptor superfamily . Recent studies have shown the presence of ERs on cells involved in immune response, namely thymocyte, macrophages and endothelial cells . Particular attention has been focused on the dose-dependent influence of estrogen on immune response, which appears to be related to the clinical symptoms of autoimmunity (i.e., the effects of pregnancy or oral contraceptive pills) . The influence of estrogens on cytokine production by target cells, through interference with their transcriptional activity, has also been the focus of various studies . The effect of estrogens on the expression of the proto-oncogenes and oncosuppressor genes involved in programmed cell death (apoptosis) might also be relevant to human autoimmunity; in particular, the uncontrolled synovial lining cell hyperplasia associated with rheumatoid arthritis and the prolonged T-cell survival in systemic lupus erythematosus may be important .
The ERβ gene, isolated by Mosselman et al. , has 96% conserved amino acid residues in the DNA binding domain and 58% conserved residues in the ligand binding domain compared with the respective residues for ERα. The ERβ microsatellite polymorphism we analyzed in the present study has been recently characterized by Tsukamoto et al. . High heterozygosity (0.93) makes this polymorphism a useful marker in the genetic study of disorders affecting female endocrine systems . Indeed, an association of some ERβ genotypes with osteoporosis has been reported . Thus, we investigated the relation between a dinucleotide (CA)n repeat polymorphism lying upstream of the human ER gene and AITDs.
Fourteen different alleles were found in 133 patients with GD, 114 patients with HT, and 179 controls subjects. The various alleles were designated as allele*1 through allele*14 according to their sizes, which ranged from 150 bp (containing 18 CA repeats with a 114 bp segment of amplified flanking sequences) to 176 bp (30 CA repeats). The distribution of ERβ alleles in the three groups is shown in Table 1. The allelic frequencies in our control population are in agreement with those previously published . Allele frequencies in our GD patients and our control subjects did not differ significantly (Χ2 = 11.71, 13 degrees of freedom, P = 0.55). Allele frequencies in our HT patients and our control subjects also did not differ significantly (Χ2 = 16.05, 13 degrees of freedom, P = 0.25). According to our power calculations, we could exclude both associations with relative risk > 0.5 for the groups analyzed in Table 1 (GD patients vs. control subjects) and associations with relative risk > 0.3 for the groups analyzed in Table 1 (HT patients vs. control subjects).
Associations of ERβ gene microsatellite polymorphism with distal radius BMD and biochemical markers of bone turnover
Ogawa et al.  demonstrated that, when they separated postmenopausal Japanese women into two groups bearing at least one allele*9 (26 CA repeat) and those who did not, the former subjects had significantly higher Z score of L2-4 BMD. Thus, Z score of distal radius BMD and bone turnover markers in 92 female GD patients in remission were compared between the groups (Table 2). There was no significant difference of the Z score of distal radius BMD and bone turnover markers between the groups.
Our data did not appear to indicate any association between a dinucleotide repeat polymorphism located in the flanking region of ERβ gene and the AITDs analyzed. However, this does not rule out the ERβ gene as a whole; other polymorphisms within the gene could still be associated with AITD. There are also limitations in using microsatellite markers with multiple allleles in case-control studies. Our result might be indicative of a large diversity in the genetic backgrounds of AITDs, although this observation deserves further analysis in a larger group of AITD patients. We also suggest that the ERβ microsatellite polymorphism has at most a minor pathogenic importance in predicting the risk of osteoporosis as a complication of GD, although further study is warranted to confirm the affects of this polymorphism on BMD and risk of osteoporosis.
Materials and methods
Materials and methods used in this paper follow closely on the ones described in a paper we have already published .
One hundred and thirty-three unrelated Japanese women with GD and 114 unrelated Japanese women with HT participated in this study. GD was diagnosed from clinical symptoms and biochemical confirmation of hyperthyroidism, including diffuse goiter, ophthalmopathy, elevated radioactive iodine uptake, and elevated thyroid hormone levels. HT patients had documented clinical and biochemical hypothyroidism requiring thyroid hormone replacement therapy and showed autoantibodies against thyroid peroxidase with or without antibodies against thyroglobulin. Remission of GD is defined as follows: both thyrotropin receptor antibody (TRAb) and thyroid-stimulating antibody (TSAb) in a negative range for 6 months, and a subsequent relatively stable euthyroid state over at least 6 months after cessation of antithyroid drug treatment. Remission ranged from 0.5 to 22.0 years. The effect of the ERβ microsatellite polymorphism on BMD in the distal radius and biochemical markers of bone turnover were investigated in the 92 GD patients optional in remission (age range, 47 to 78 years), none of whom had been treated with radioiodine, surgery, or drugs that could affect bone turnover or BMD prior to the study. One hundred and seventy-nine unrelated Japanese women without clinical evidence or family history of any autoimmune diseases were selected as normal control subjects. All study participants lived in the Tokyo metropolitan area. The research protocol was approved by the ethics committee of Showa University hospital, and informed consent was obtained from all patients and controls.
Measurement of biochemical markers
Blood and urine samples were collected early in the morning after an overnight fast. Both serum and urine were stored at -70°C until assayed. Serum concentration of bone alkaline phosphatase (B-ALP) was measured by densitometric scanning following polyacrylamide disc gel electrophoresis (PAGE). Osteocalcin (OC) was measured in serum with an IRMA. Materials for both assays were obtained from Mitsubishi Kagaku BCL (Tokyo, Japan). Intra- and interassay coefficients of variation were less than 10%. Urinary hydroxyproline (Hp) and deoxypyridinoline (Dpyr) were measured with high-performance liquid chromatography (HPLC; Japan Spectroscopic, Tokyo, Japan). Intra-assay coefficients of variation were less than 2.0% for Hp and 5.8% for Dpyr. Interassay coefficients of variation were less than 0.2% for Hp and 6.8% for Dpyr. Values of urinary Hp and Dpyr were corrected for urinary creatinine concentration, which was measured with an autoanalyzer.
Bone mineral density analysis
BMD was measured at the distal one third of the left radius with dual energy X-ray absorptiometry (DEXA; DTX-200, Osteometer, Copenhagen, Denmark). With this equipment, coefficients of variation in our laboratory were 1.0%. BMD was expressed as a Z-score representing the number of standard deviations (SDs) by which a value deviated from the mean BMD in 4095 healthy Japanese female controls aged 20 to 90 years.
Determination of microsatellite polymorphism by polymerase chain reaction
Genomic DNA was isolated from whole blood with a Genomix kit (Talent, Trieste, Italy). Microsatellite marker loci were typed with fluorescence-based methods [23,24]. Polymerase chain reaction (PCR) was performed with oligonucleotide primers designed to amplify a polymorphic (CA)n repeat in the flanking region of the human ERβ gene . PCR was performed in a total volume of 20 μl of the following mixture: 100 ng of human genomic DNA, 5 pmol of each primer (Cy-5'-GGTAAACCATGGTCTGTACC-3' and 5'-AACAAAATGTTGAATGAGTGGG-3'), 200 μM of each dNTPs, 2 μl of 10 x reaction buffer (Takara Shuzo Co., Kyoto, Japan); and 1 unit of Taq DNA polymerase (Takara Shuzo Co.). Thirty PCR cycles of 30 sec at 94°C, 30 sec at 61°C, and 30 sec at 72°C were performed. The amplified fragments were analyzed by electrophoresis with an automatic DNA sequencer with 8% or 6% polyacrylamide gels containing 7 M Urea (Amersham Pharmacia Biotech, Buckingamshire, UK). The raw data were converted to dinucleotide repeat polymorphic band patterns with the use of a software program (Amersham Pharmacia Biotech). The number of CA repeats in each amplified product was determined by comparison with a ladder of control DNAs.
Comparisons between the various alleles in patients with AITDs and in controls were made with the Χ2- test. Fisher's exact test was used when necessary. The relative risk was calculated by Woolf's method . Comparison of BMD in Z scores and biochemical markers between the group of individuals possessing one or two alleles of the genotype and the group without that genotype were performed using Student's t test. A P value less than 0.05 was considered statistically significant.
Ban Y, Taniyama M, Tozaki T, Tomita M, Ban Y: Estrogen receptor alpha dinucleotide repeat polymorphism in Japanese patients with autoimmune thyroid diseases. BMC Med Genet. 2000, 1: 1-10.1186/1471-2350-1-1.
Tomer Y, Barbesino G, Greenberg DA, Davies TF: The immunogenetics of autoimmune diabetes and autoimmune thyroid disease. Trends Endocrinol Metab. 1997, 8: 63-70. 10.1016/S1043-2760(96)00266-4.
Tomer Y, Davies TF: Infection, thyroid disease and autoimmunity. Endocr Rev. 1993, 14: 107-120. 10.1210/er.14.1.107.
DeGroot LJ, Quintans J: The causes of autoimmune thyroid disease. Endocr Rev. 1989, 10: 537-562.
Mckenzie JM, Zakarija M: Hyperthyroidism. In Endocrinology, 2nd ed. Edited by Degroot LJ, Besser GM, Cahill GF, et al. Philadelphia: Saunders,. 1988, 646-682.
Volpe R: Graves' disease. In Werner and Ingbar's the Thyroid: A Fundamental and Clinical Text. Edited by Braverman LE, Utiger RD. Philadelphia: Lippincott-Raven,. 1991, 648-657.
Yanagawa T, Mangklabruks A, Chang YB, Okamoto Y, Fisfalen ME, Curran PG, DeGroot LJ: Human histocompatibility leukocyte antigen DQA1 *0501 allele associated with genetic susceptibility to Graves' disease in a Caucasian population. J Clin Endocrinol Metab. 1993, 76: 1569-1574. 10.1210/jc.76.6.1569.
Badenhoop K, Walfish PG, Rau H, Fischer S, Nicolay A, Bogner U, Schleusener H, Usadel KH: Susceptibility and resistance alleles of human leukocyte antigen (HLA) DQA1 and HLA DQB1 are shared in endocrine autoimmune disease. J Clin Endocrinol Metab. 1995, 80: 2112-2117. 10.1210/jc.80.7.2112.
Yanagawa T, Hidaka Y, Guimaraes V, Soliman M, DeGroot LJ: CTLA-4 gene polymorphism associated with Graves' disease in a Caucasian population. J Clin Endocrinol Metab. 1995, 80: 41-45. 10.1210/jc.80.1.41.
Hill SM, Fuqua SAW, Chamness GC, Greene GL, McGuire WL: Estrogen receptor expression in human breast cancer associated with an estrogen receptor gene restriction fragment length polymorphism. Cancer Res. 1989, 49: 145-148.
Yaich L, Dupont WD, Cavener DR, Parl FF: Analysis of the PvuII restriction fragment-length polymorphism and exon structure of the estrogen receptor gene in breast cancer and peripheral blood. Cancer Res. 1992, 52: 77-83.
Lehrer S, Rabin J, Kalir T, Schachter BS: Estrogen receptor variant and hypertension in women. Hypertension. 1993, 21: 439-441.
Kobayashi S, Inoue S, Hosoi T, Ouchi Y, Shraki M, Orimo H: Association of bone mineral density with polymorphism of the estrogen receptor gene. J Bone Miner Res. 1996, 11: 306-311.
Sano M, Inoue S, Hosoi T, Ouchi Y, Emi M, Shiraki M, Orimo H: Association of estrogen receptor dinucleotide repeat polymorphism with osteoporosis. Biochem Biophys Res Commun. 1995, 217: 378-383. 10.1006/bbrc.1995.2787.
Ushiyama T, Ueyama H, Inoue K, Nishioka J, Ohkubo I, Hukuda S: Estrogen receptor gene polymorphism and generalized osteoarthritis. J Rheumatol. 1998, 25: 134-137.
Ushiyama T, Mori K, Inoue K, Huang J, Nishioka J, Hukuda S: Association of oestrogen receptor gene polymorphisms with age at onset of rheumatoid arthritis. Ann Rheum Dis. 1999, 58: 7-10.
Mosselman S, Polman J, Dijkema R: ER-beta: identification and characterization of a novel human estrogen receptor. FEBS Lett. 1996, 392: 49-53. 10.1016/0014-5793(96)00782-X.
Ogawa S, Hosoi T, Shiraki M, Orimo H, Emi M, Muramatsu M, Ouchi Y, Inoue S: Association of estrogen receptor β gene polymorphism with bone mineral density. Biochem Biophys Res Commun. 2000, 269: 537-541. 10.1006/bbrc.2000.2285.
Tomer Y, Barbesino G, Greenberg DA, Concepcion E, Davies TF: Linkage analysis of candidate genes in autoimmune thyroid disease. III. Detailed analysis of chromosome 14 localizes GD-1 close to MNG-1. J Clin Endocrinol Metab. 1998, 83: 4321-4327. 10.1210/jc.83.12.4321.
Evans RM: The steroid and thyroid hormone receptor superfamily. Science. 1988, 240: 889-895.
Cutolo M, Sulli A, Seriolo B, Accardo S, Masi AT: Estrogens, the immune response and autoimmunity. Clin Exp Rheumatol. 1995, 13: 217-226.
Tsukamoto K, Inoue S, Hosoi T, Orimo H, Emi M: Isolation and radiation hybrid mapping of dinucleotide repeat polymorphism at the human estrogen receptor β locus. J Hum Genet. 1998, 43: 73-74. 10.1007/s100380050043.
Reed PW, Davies JL, Copeman JB, Bennett ST, Palmer SM, Pritchard LE, Gough SCL, Kawaguchi Y, Cordell HJ, Balfour KM, et al: Chromosome-specific microsatellite sets for fluorescence-based, semi-automated genome mapping. Nat Genet. 1994, 7: 390-395.
Pritchard LE, Kawaguchi Y, Reed PW, Copeman JB, Davies JL, Barnett AH, Bain SC, Todd JA: Analysis of the CD3 region and type 1 diabetes: application of fluorescence-based technology to linkage disequilibrium mapping. Hum Mol Genet. 1995, 4: 197-202.
Woolf B: On estimating the relation between blood group and disease. Ann Hum Genet. 1955, 19: 251-255.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/content/backmatter/1471-2350-2-1-b1.pdf
This work was supported in part by the High-Technology Research Center Project from the Ministry of Education, Science, Sports and Culture of Japan.