Positive association of the DIO2 (deiodinase type 2) gene with mental retardation in the iodine-deficient areas of China
- T-W Guo1,2,*,
- F-C Zhang3,*,
- M-S Yang1,2,
- X-C Gao3,
- L Bian1,2,
- S-W Duan1,2,
- Z-J Zheng3,
- J-J Gao1,2,
- H Wang1,
- R-L Li4,
- G-Y Feng5,
- D St Clair6,
- L He1,2
- 1Institute for Nutritional Science, Chinese Academy of Sciences, 320 Yue Yang Road, Shanghai 200031, China
- 2Bio-X Life Science Research Center, Shanghai Jiao Tong University, Shanghai 200030, China
- 3Institute of Population and Health, Northwest University, Xi’an 710069, China
- 4The Second Hospital, Xi’an Jiao Tong University, Xi’an 710049, China
- 5Shanghai Institute of Mental Health, 600 South Wan Ping Road, Shanghai 200030, China
- 6Department of Mental Health, University of Aberdeen, Medical School, Foresterhill, Aberdeen, AB25 2ZD, UK
- Correspondence to: L He Bio-X Life Science Research Center, PO Box 501, Hao Ran Building, Shanghai Jiao Tong University, 1954 Hua Shan Road, Shanghai 200030, China;
- Received 18 March 2004
- Accepted 24 March 2004
Background: Iodine deficiency is the commonest cause of preventable mental retardation (MR) worldwide. However, in iodine-deficient areas not everyone is affected and familial aggregation is common. This suggests that genetic factors may also contribute. Thyroid hormone (TH) plays an important role in fetal and early postnatal brain development. The pro-hormone T4 (3,3′,5,5′-triiodothyronine) is converted in the brain to its active form, T3, or its inactive metabolite, reverse T3, mainly by the action of deiodinase type 2 (DIO2).
Methods: To investigate the potential genetic contribution of the DIO2 gene, we performed a case-control association study using three common SNPs in the gene (rs225014, rs225012, and rs225010) that were in strong linkage disequilibrium with each other.
Results: Single marker analysis showed a positive association of MR with rs225012 and rs225010. Particularly with rs225012, TT genotype frequency was significantly higher in MR cases than in controls (χ2 = 9.18, p = 0.00246). When we compared the distributions of common haplotypes, we also found significant differences between mental retardation and controls in the haplotype combination of rs225012 and rs225010 (χ2 = 15.04, df 2, global p = 0.000549). This association remained significant after Bonferroni correction (p = 0.0016470).
Conclusion: We conclude that allelic variation in the DIO2 gene may affect the amount of T3 available and in an iodine-deficient environment may partly determine overall risk of MR.
- LD, linkage disequilibrium
- MR, mental retardation
- PCR, polymerase chain reaction
- TH, thyroid hormone
Fetal iodine deficiency is the commonest cause of preventable mental retardation.1 Every year, 100 000 children are born with frank cretinism, and many times more are born with lesser mental and neurological deficits attributable to iodine deficiency,2 as a result of inadequate amounts of thyroid hormone (TH) being available to the developing fetal brain. Thyroid hormones regulate the processes of terminal brain differentiation such as dendritic and axonal growth, synaptogenesis, neuronal migration and myelination. TH also modulates the establishment of neuronal networks through regulation of the number of microglial cells producing neurotrophic factors.3
The deiodinase play a key role in the maintenance of circulating and tissue levels of thyroid hormones. There are three types of deiodinase, type 1, 2, and 3 (DIO1, DIO2, and DIO3) iodothyronine. All are seleno-enzymes characterised by a selenocysteine in the catalytic domain of the enzyme encoded by a UGA codon in the presence of a characteristic 3′ untranslated region stem loop structure, the selenocysteine insertion sequence (SECIS). DIO2 is particularly important in the brain. The pro-hormone T4 (3,3′,5,5′-triiodothyronine) is converted in the brain to its active form, T3, or its inactive metabolite, reverse T3, mainly by the action of DIO2. DIO2 is involved in an activation step converting 3,5,3′5′-tetra-iodothyronine (T4) to 3,5,3′-tri-iodothyronine (T3), and a degradation step converting 3,3′5′-triiodothyronine (reverse T3) to 3,3′-diiodothyronine. Inactivation steps are mainly regulated by deiodinase type 3 (DIO3) when converting T4 to reverse T3 and converting T3 to 3,3′-diiodothyronine.4 DIO2 appears to be a tissue-specific regulator of intracellular T3 concentrations in the brown fat, brain, and pituitary.5 The DIO2 gene maps to human chromosome 14q24.3, is about 15 kb in size, and the coding region is divided into two exons by a gap of approximately 7.4 kb.6 In light of the important action of thyroid hormone in brain development and the regulation of the active form of thyroid hormone by DIO2 in the brain, and as, up to now, there have few genetic studies performed in this specific field, we decided to evaluate whether allelic variation in the DIO2 gene might alter risk of susceptibility for MR in areas of iodine deficiency.
The study included three groups, definite mental retardation (MR; n = 96), borderline mental retardation (border; n = 116), and controls (n = 331). The mean (SD) age of groups was 9.9 (2.9) years with a 49:51 female–male ratio (table 1). All subjects were identified and recruited from Zha Shui and An Kang counties in the Qin-Ba mountain region of Shaanxi province, Western China. There is widespread soil erosion in this region (average elevation 750–1500 m) and water iodine levels are low; we recorded an iodine (SD) level in water of 1.87 (0.46) μg/l.7 There is no selenium deficiency. For several years the Shaanxi province health authorities have conducted iodination programs and popularised the use of iodine enriched common salt. This has resulted in a marked decline in the frequency of mental retardation in the Qin-Ba mountain region. However, the prevalence of mental retardation (2.78%) still remains higher than in most other areas of China (1.07%).8 Moreover, we found familial clustering in the two counties with several families displaying multiple affected members in one or more generations. It is possible that genetic factors may interact with an iodine-deficient physical environment to determine overall risk of mental retardation.
All subjects gave standard informed consent after explanation of the study. The protocol was reviewed and approved by the Ethical Committee of the National Human Genome Center. All subjects were Han Chinese in origin.
Screening for social adaptability or mental handicap
Participants were screened using the Adaptive Scale of Infant and Children revised by Zuo et al.9 Using these scales each person was given a social adaptive score or mental handicap score. Those with no disability on these scales invariably have an IQ within normal range and, for purposes of the study, were therefore classified as normal controls.
Children 4–5 years old were tested with the Chinese-Wechsler Young Children Scale of Intelligence (C-WYCSI),10 while those 6–16 years old were tested with the Chinese-Wechsler Intelligence Scale for Children (C-WISC).11 We selected an IQ of less than 70 as the cut-off for mental retardation (MR). We defined IQs of less than 70 accompanied by social disability scores of 8 or less as mental retardation (MR), and IQs of 70–79 with social disability scores of 9 as borderline MR (border).
It was usually not possible or appropriate to perform a formal IQ test in cases of frank cretinism, and we had to depend upon clinical diagnosis. A neurological examination, conducted by a physician, included tests of hearing, vision, voice and speech, reflexes, and posture and gait. We excluded from the study cases of MR if affected by trachoma, infection, trauma, toxicity, cerebral palsy, or birth complications. Controls came from the same iodine-deficient areas and were selected from families with no history of MR. If permission was granted, a blood sample was taken for routine haematology, serology, and DNA analysis. Blood samples were stored at -20°C. Genomic DNA was extracted from blood using a modified phenol/chloroform method.
We performed power calculations based on Cohen’s method12 using G*Power software and Epi info 2002. The present sample size showed over 90% power to detect significance (α<0.05) in the association with allele, genotype, and haplotype under study conditions and an effect size index of 0.2 corresponding to “weak to moderate” gene effect was used. Furthermore, the present sample sizes had a power of >73% at p = 0.05 to detect allelic association of the allele, genotype, and haplotype at a presumed odds ratio (OR) of 2.
SNPs in the DIO2 region
We selected SNPs located in the DIO2 region from information sourced from the SNP Consortium (http://snp.cshl.org/) and dbSNP (http://www.ncbi.nlm.nih.gov/SNP/) and examined their allele frequencies in 24 controls (48 chromosomes) by the direct sequencing procedure described below. After evaluation of 22 SNPs by seven pairs of sequencing primer, we selected three SNPs with minor allele frequencies of over 0.05 in the DIO2 gene. They were rs225014 (A/G) in exon 2, and rs225012 (T/C) and rs225010 (A/G) in intron 1. rs225014 and rs225012 have an interval of 1.2 kbp and rs225012 and rs225010 have an interval of 1.5 kbp.
The SNP rs225014 (A/G) in exon 2 was amplified by polymerase chain reaction (PCR) using the primers: forward: 5′-TACCTGCCATCATGCCTCTT-3′, and reverse: 5′-GGAAGTCAGCCACTGAGGAG -3′. PCRs were carried out in 96-well microtitre plates with a final reaction volume of 25 μl containing 50 mM KCl, 10 mM Tris-HCl (pH 8.0), 1.5 mM MgCl2, 200 mM dNTPs, 5 μl Q solution (Qiagen, Valencia, CA, USA), 10 pM each primer, 20 ng DNA, and 2.5 U Taq polymerase (Life Technologies, Karlsuhe, Germany). Cycle conditions were one cycle with an initial 4 min denaturation at 95°C, followed by 35 cycles of 94°C for 30 s, 57°C for 40 s, 72°C for 50 s, and a final extension period at 72°C for 10 min, using the GeneAmp PCR System 9700 (Applied BioSystems, Foster City, CA). SNPs were typed by direct sequencing. The PCR products were processed by incubation with 0.1 U of shrimp alkaline phosphatase (Roche, Basel, Switzerland) and 0.5 U of exonuclease I (New England Biolabs, Beverly, MA) at 37°C for 1 h, followed by heat inactivation at 80°C for 20 min. The PCR products were sequenced with reverse PCR primer as the sequencing primer using an ABI Prism BigDye Terminator Cycle Sequencing Kit (Applied BioSystems, Foster City, CA) on an ABI Prism 377 or 3100 sequencer.
The genotyping of SNPs rs225012 (T/C) and rs225010 (A/G) combines kinetic (real-time quantitative) PCR with allele-specific amplification, which has been described elsewhere.13 We used two separate real-time quantitative PCR reactions, each of which contains an allele-specific primer of SNP and the same common primer. Heterozygous samples have equal amounts of the two alleles, which should reach a detectable level of fluorescence at the same cycle number, but in heterogeneity the cycle number should be different for the two amplification reactions. For the rs225012 (T/C) polymorphism a 55 bp PCR product was amplified using the common primer: 5′-GTGCAAAGGGAGCACATGAA-3′, a T allele-specific primer: 5′-TAAAATATTGGGGCAGAAGAA-3′, and a C allele-specific primer: 5′-TAAAATATTGGGGCAGAAGAG-3′. Amplification was performed in a 5 μl volume of 2.5 μl 2×TaqMan universal PCR master mix (Applied Biosystems), 10 ng genomic DNA, 0.2 μM allele-specific primer, 0.2 μM common primer, and 0.2×SYBR Green I (Molecular Probe). The PCR cycles began with an UNG digestion stage at 50°C for 2 min, and an initial denaturation period at 95°C lasting for 10 min, followed by 50 cycles at 95°C for 15 s, an annealing phase conducted at 59°C for 30 s, and a dissociation stage at 95°C for 15 s. For the rs225010 (A/G) polymorphism a 55 bp PCR product was amplified under the same cycling conditions except with the following primers: the common primer: 5′-AAATTTATCTGCCTCAGTGAGCTT-3′, an allele-specific primer: 5′- GAACATAATCATATTTGGGTGAT-3′, and a G allele-specific primer: 5′-GAACATAATCATATTTGGGTGAC-3′. To check for genotyping errors, eight DNA samples were randomly selected from each 96-well plate and re-genotyped. All genotypes were identical to those obtained from the first round of genotyping.
Allele frequencies were calculated using the SPSS 10.0 software for Windows (SPSS, Chicago, IL). Deviations from Hardy-Weinberg equilibrium, differences in allele and genotype distributions, and OR with 95% confidence intervals were calculated using the method of Finetti.14 Linkage disequilibrium (LD) between two loci was measured using a two-locus LD calculator (2LD)15 and EMLD software (http://request.mdacc.tmc.edu/~qhuang/Software/pub.htm). Haplotypes were inferred by Bayesian methods,16 and implemented in the PHASE package version 1.0 (http://www.assertion.net/software/). Differences in genotype and haplotype distribution between patient and control groups were assessed by the Monte Carlo method using the CLUMP program version 1.9 with 10 000 simulations.17 Statistical significance was set at p<0.05. Odds ratios with 95% confidence intervals were estimated for the effects of high-risk haplotype and calculated by Epi info 2002 software (http://www.cdc.gov/epiinfo/). Power analysis was performed using the G*Power program18 and Epi info 2002 software.
The rs225014 SNP is an A/G polymorphism in exon 2 of the DIO2 gene, predicting a change in amino acid 92 of the protein (Thr92Ala). The minor allele G frequency of rs225014 was 0.3841 in the control population with a distribution meeting Hardy-Weinberg equilibrium (p = 0.201). The rs225012 and rs225010 SNPs are in intron 1 of the DIO2 gene. rs225012 is a T/C SNP with a minor allele T frequency of 0.2935 in our control population and an allele frequency of 0.272 in the Caucasian population. rs225010 is an A/G SNP with a minor allele G frequency of 0.337 in our control population. The distributions of genotypes of rs225012 and rs225010 were both in Hardy-Weinberg equilibrium in controls. All samples were grouped together for statistical analysis because no significant difference was found in distribution of genotype frequencies between sample from Zha Shui and An Kang counties (p>0.35).
Singular-locus association analysis
Table 2 shows the group genotypes and allele frequencies of the three SNPs. No significant difference was observed in genotypes or allele frequencies for the three SNPs between borderline mental retardation and control groups. When MR and controls were compared, rs225014 also revealed no significant difference in allele frequency (p = 0.5) or distribution of genotypes (p = 0.39). The G allele frequency of rs225010 was higher in MR (0.337) than in controls (0.2576; p = 0.045; OR = 1.43, 95% confidence interval 1.006 to 2.031). The T allele frequency of rs225012 was higher in MR (0.2935) than in controls (0.2119; p = 0.03; OR = 1.491, 95% confidence interval 1.033 to 2.152). The CC genotype frequency was significantly higher in MR (0.128) than in controls (0.043; χ2 = 9.18, p = 0.00246; OR = 3.29, 95% confidence interval 1.37 to 7.91).
To calculate the extent of LD in pairwise combinations of the three SNPs, we calculated D′, r2, and the p value, the normalised linkage disequilibrium statistic, in controls for all possible pairs of SNPs. The pairwise LD values are shown in table 3. Strong linkage disequilibrium among the three SNPs was observed (all D′>0.7, p<0.00001).
We constructed four sets of haplotypes. Three were derived from various combinations of two SNPs and one was derived from a combination of all three SNPs. All the haplotypes were estimated with the PHASE package (table 4).
We only found significant differences between controls and mental retardation in the haplotype of combination of rs225012 and rs225010 (χ2 = 15.04, df 2, global p = 0.000549). Furthermore, the data obtained from the study of haplotypes containing either rs225012C or rs225010A showed that the frequency of haplotype C-A was much lower in MR than in controls (χ2 = 19.36, df 1, p = 0.00001; OR = 0.49, 95% confidence interval 0.35 to 0.68). The three marker haplotypes GTA, GTG, GCA are all risk haplotypes for mental retardation, but haplotype GTG has the highest odds ratio (OR = 12.49, 95% confidence interval: 2.36 to 87.6).
In this study, we investigated the relationship between three polymorphisms in the DIO2 gene and mental retardation in a Chinese Han population from the Qin-Ba mountain region, a traditionally iodine-deficient area in northwest China. This study has several strengths. First, we showed that our sample size had reasonable power to detect association even when these variants had small to medium effects (effect size = 0.2) on susceptibility. Furthermore, our child samples may be less influenced by non genetic social and culture factors than adult samples. Second, all subjects came from the Qin-Ba mountain region, a relatively isolated area in northwest China. Moreover, we did not find a significant difference in allele frequencies between the two counties, which reduced the risk of stratification bias. Third, according to the definition of mental retardation,19 significantly subaverage intellectual functioning with an IQ score of 70–75 or below on a standardised individual intelligence test is classified as MR. In order to ensure the diagnosis of mental retardation, we set 70 as the IQ cut-off point of MR and classified the individuals identified with IQs of 70–79 as having borderline mental retardation. Finally we have obtained a significant result even after Bonferroni correction. rs225014 (A/G) is a coding polymorphism (Thr92Ala) predicting a change in amino acid 92 of the protein. rs225012 (C/T) and rs225010 (A/G) are located in intron 1.
Single-locus analysis showed a positive association of MR with markers rs225012 and rs225010. When we compared the distributions of common haplotypes between control and MR, we also found significant differences between controls and mental retardation in the haplotype with combined rs225012 and rs225010. This association was still significant (p<0.001647) after Bonferroni correction.20 The haplotypes of rs225012C and rs225010A showed that C-A was much more frequent in controls than in MR (χ2 = 19.36, df 1, p = 0.00001) and suggests a protective effect. We analysed whether or not two SNPs were located near an mRNA splicing site through in silico analysis by GENESPLICER (developed by the Institute for Genomic Research (TIGR), http://www.tigr.org/tdb/GeneSplicer/gene_spl.html)21 and found they were not.
Another possibility is that rs225012 and rs225010 and the haplotype combination with these two SNPs may simply be in linkage disequilibrium with a functional polymorphism elsewhere in the DIO2 gene or in a gene nearby. The nearest functional polymorphism in the DIO2 gene is rs225014 (A/G), a common non conservative variant which predicts a Thr92Ala substitution. Although the crystal structure of type 2 deiodinase is not yet known, it is worth noting that this non conservative amino acid change (aliphatic for polar group), which is not located within the conserved deiodinase catalytic domain, could potentially affect its activity.22 However, this region of the enzyme is not phylogenetically conserved. The homologous amino acid is represented by a proline in rodents and by a glycine in chick. In contrast, humans and amphibians share a threonine in this position. It is reported that this SNP is associated with obesity and insulin resistance.23 However, we found no association between this polymorphism and mental retardation.
The nearest functional candidate genes are PSEN1 (presenilin1) and TSHR (thyroid hormone receptor or thyrotropin receptor). PSEN1 is reported to be associated with early-onset Alzheimer disease type 3.24TSHR is a important gene involved both in the metabolic pathway of thyroid hormones and in a wide range of sporadic and hereditary or genetically determined changes in thyroid function such as thyroid adenomas, thyroid cancer, non autoimmune hyperthyroidism, thyrotropin resistance, and congenital hypothyroidism.25 Therefore, the effect of the combination of DIO2, PSEN1, and TSHR may play an important role in the pathogenesis of mental retardation.
It is also possible that the protective haplotype is in linkage disequilibrium with a regulatory element that affects expression of type 2 deiodinase. This in turn may influence T3 levels in the fetal brain and, if the brain is already compromised by iodine deficiency, influence the risk of mental retardation. Eighty per cent of brain T3 is formed through the enzymatic activity of type 2 deiodinase (DIO2). DIO2 is mainly found in astrocytes in vitro and vivo,26 suggesting that circulating T4 is metabolised into T3 in the glial cells and then transferred to the neurons. To examine the physiological role of DIO2, Schneider et al27 developed a DIO2 knockout mouse strain lacking DIO2 activity. Mice homozygous for the targeted deletion had no gross phenotypic abnormalities, and development and reproductive function appeared normal, except for mild growth retardation in males. It is unclear whether a similar situation pertains for DIO2 deficiency in human fetal brain whose thyroid hormone function is already compromised by iodine deficiency. We suspect that normal deiodinase type 2 may be an important protection factor from late in the first trimester and early in the second trimester of gestation, because in this period DIO2 activity in the brain is evidently increased.28 Further work is required to investigate the mechanisms by which DIO2 may affect fetal brain development in the context of iodine deficiency.
↵* The first two authors contributed equally to this work.
This work was supported by grants from the national 973 and 863 programs, the National Natural Science Foundation of China, and the Shanghai Municipal Commission for Science and Technology.
Conflict of interest: none declared.