Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Association between Thyroid Function and Ocular Parameters

During development, thyroid hormones play an important role in eye development, while in adults, some pathological thyroid conditions can affect the normal functioning of the eyes. Thyroid eye disease is the most well-known eye pathology caused by a pathological thyroid condition. Few studies have investigated the association between ocular parameters and thyroid function. Thus, in this study, we aimed to examine whether thyroid activity affects ocular parameters. This cross-sectional study included 4633 healthy adults recruited within the 10,001 Dalmatians project of the Croatian Biobank. The plasma levels of thyroid-stimulating hormone (TSH), free triiodothyronine (fT3), free thyroxine (fT4), thyroglobulin (Tg), thyroglobulin antibodies (TgAb), and thyroid peroxidase antibodies (TPOAb) were measured by an immunoassay. We determined 20 ocular parameters for each participant (10 for each eye, including corneal radius, corneal thickness, anterior chamber depth, anterior chamber angle, lens thickness, posterior chamber length, axial length, intraocular lens power (IOL), spherical power, and cylinder power). Patients with hyperthyroidism had thicker corneas compared to euthyroid individuals. Corneal thickness was also negatively associated with plasma TSH levels. Intra-ocular lens power was higher in patients with clinical hypothyroidism, while spherical power was higher in euthyroid individuals with positive antibodies compared to euthyroid individuals. Intra-ocular lens power negatively correlated with fT4 levels, while spherical power positively correlated with TgAb, TPOAb, and Tg levels and negatively correlated with TSH levels. The anterior chamber angle was positively associated with plasma TSH levels and TPOAb levels and negatively associated with plasma fT4 levels. These findings suggest an interesting interplay between ophthalmic measures and thyroid status, detectable even in the general adult population.

Thyroid Hormones Are Not Associated with Plasma Osteocalcin Levels in Adult Population with Normal Thyroid Function

Thyroid hormones (THs) play an indispensable role in skeletal development and bone remodeling. Some studies have reported associations of THs with serum osteocalcin (OC) levels, but the results are quite inconsistent and the molecular mechanism of their simultaneous or interdependent activity on bone is almost unknown. Therefore, the aim of this study was to determine the possible associations of plasma THs with plasma OC levels and the possible mediating effect of OC on the relationship between THs and bone mineral density (BMD). For this purpose, out of the initial 1981 participants, we selected healthy euthyroid participants controlled for available confounding factors that can affect thyroid function and bone metabolism (N = 694). Given our results, we could not confirm any associations of THs with plasma OC levels nor the mediating effect of OC on the relationship between THs and BMD in euthyroid population. In the group of women controlled for menopause status (N = 396), we found a significant negative association of body mass index (BMI) with OC levels (β = −0.14, p = 0.03). We also found a negative association of free triiodothyronine (fT3) (β = −0.01, p = 0.02) and age (β = −0.003, p < 0.001) with BMD, and a positive association of BMI (β = 0.004, p < 0.001) and male gender (β = 0.1, p < 0.001) with BMD. In addition, we found significantly higher plasma OC levels and lower values of BMD in postmenopausal euthyroid women compared with premenopausal euthyroid women. In our opinion, the results of previous studies suggesting an association between circulating THs and serum OC levels may be influenced by an inconsistent selection of participants and the influence of confounding factors.

Epidemiology of Hypothyroidism, Hyperthyroidism and Positive Thyroid Antibodies in the Croatian Population

Simple Summary

The aim of this paper was to analyse the prevalence of hypothyroidism and hyperthyroidism in the Croatian population. This is the first epidemiological study of its kind conducted in our population. It is important to correctly diagnose thyroid dysfunction due to the detrimental effects of thyroid dysfunction on human health (especially in pregnant women and patients with cardiovascular diseases). The results of our study conducted on 4402 individuals, showed a higher prevalence of hypothyroidism in our country (10.5%) compared to other countries, while the prevalence of hyperthyroidism was quite similar (1.3%). We also observed that a high proportion of thyroid disorders remained undiagnosed (as many as 92.6% subclinical hypothyroid, 93.9% clinical hypothyroid, 83% subclinical hyperthyroid and 71.4% clinical hyperthyroid participants). Therefore, our study indicates that investing in prevention programs is crucial.

Abstract

Thyroid dysfunction appears to be the leading endocrine disorder. We conducted a cross-sectional study on 4402 individuals from three Croatian cohorts. The aim of this study was to analyse the prevalence of diagnosed and undiagnosed hypothyroidism, hyperthyroidism (subclinical and clinical) and positive thyroid antibodies in the Croatian population. The results of the study indicated that 17.6% of participants were euthyroid with positive antibodies. The prevalence of clinical and subclinical hypothyroidism was 3% and 7.4%, respectively, while the prevalence of clinical and subclinical hyperthyroidism was 0.2% and 1.1%, respectively. Among them, 92.6% subclinical hypothyroid, 93.9% clinical hypothyroid, 83% subclinical hyperthyroid and 71.4% clinical hyperthyroid participants were undiagnosed. Finally, the prevalence of undiagnosed subclinical and clinical hypothyroidism in our population was 6.9% and 2.8%, respectively, while the prevalence of undiagnosed subclinical and clinical hyperthyroidism was 0.9% and 0.1%, respectively. Women showed a higher prevalence of thyroid disorders; 1.57 times higher odds of euthyroidism with positive antibodies, 2.1 times higher odds of subclinical hyperthyroidism, 2.37 times higher odds of clinical hypothyroidism and 1.58 times higher odds of subclinical hypothyroidism than men. These results indicate an extremely high proportion of undiagnosed cases, and therefore require investments in a prevention programme.

Genome-Wide Association Analysis and Genomic Prediction of Thyroglobulin Plasma Levels

Thyroglobulin (Tg) is an iodoglycoprotein produced by thyroid follicular cells which acts as an essential substrate for thyroid hormone synthesis. To date, only one genome-wide association study (GWAS) of plasma Tg levels has been performed by our research group. Utilizing recent advancements in computation and modeling, we apply a Bayesian approach to the probabilistic inference of the genetic architecture of Tg. We fitted a Bayesian sparse linear mixed model (BSLMM) and a frequentist linear mixed model (LMM) of 7,289,083 variants in 1096 healthy European-ancestry participants of the Croatian Biobank. Meta-analysis with two independent cohorts (total n = 2109) identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) within the ST6GAL1 gene (p<5×10-8). BSLMM revealed additional association signals on chromosomes 1, 8, 10, and 14. For ST6GAL1 and the newly uncovered genes, we provide physiological and pathophysiological explanations of how their expression could be associated with variations in plasma Tg levels. We found that the SNP-heritability of Tg is 17% and that 52% of this variation is due to a small number of 16 variants that have a major effect on Tg levels. Our results suggest that the genetic architecture of plasma Tg is not polygenic, but influenced by a few genes with major effects.

Thyroid hormone levels are associated with metabolic components: a cross-sectional study

Aim

Methods

Results

Conclusions

Genome-Wide Analysis Identifies Two Susceptibility Loci for Positive Thyroid Peroxidase and Thyroglobulin Antibodies

INTRODUCTION

Thyroid peroxidase (TPO) and thyroglobulin (Tg) are main components of the thyroid gland and play an essential role in thyroid hormone synthesis. The development of antibodies to thyroid peroxidase (TPOAb) and thyroglobulin (TgAb) is the major diagnostic hallmark and early indicator of autoimmune thyroid disease. TPOAb and TgAb are under strong genetic influence; however, genetic factors that determine thyroid antibody positivity are largely unknown.

MATERIALS AND METHODS

To identify novel loci associated with TPOAb and/or TgAb positivity, we performed a genome-wide meta-analysis in a total of 2613 individuals from Croatia. Participants with elevated plasma TPOAb and/or TgAb were defined as cases (N = 619) and those with TPOAb and TgAb within reference values were defined as controls (N = 1994).

RESULTS

We identified 2 novel loci, of which 1 is located within the YES1 gene (rs77284350, P = 1.50 × 10-8), and the other resides within the IRF8 gene (rs16939945, P = 5.04 × 10-8).

CONCLUSIONS

Although the observed variants were associated with TPOAb and TgAb positivity for the first time, both YES1 and IRF8 were previously linked to susceptibility to other autoimmune diseases, and represent plausible biological candidates. This study adds to the knowledge of genetics underlying thyroid antibodies and provides a good basis for further research.

Genome-wide meta-analysis identifies novel loci associated with free triiodothyronine and thyroid-stimulating hormone

PURPOSE

Thyroid hormones are essential for the normal function of almost all human tissues, and have critical roles in metabolism, differentiation and growth. Free triiodothyronine (fT3), free thyroxine (fT4) and thyroid-stimulating hormone (TSH) levels are under strong genetic influence; however, most of the heritability is yet unexplained.

METHODS

In order to identify novel loci associated with fT3, fT4 and TSH serum levels we performed a genome-wide meta-analysis of 7 411 206 polymorphisms in up to 1731 euthyroid individuals from three Croatian cohorts from Dalmatia region: two genetically isolated island populations and one mainland population. Additionally, we also performed a bivariate analysis of fT3 and fT4 levels.

RESULTS

The EPHB2 gene variant rs67142165 reached genome-wide significance for association with fT3 plasma levels (P = 9.27 × 10^-9) and its significance was confirmed in bivariate analysis (P = 9.72 × 10^-9). We also found a genome-wide significant association for variant rs13037502 upstream of the PTPN1 gene and TSH plasma levels (P = 1.67 × 10^-8).

CONCLUSION

We identified a first genome-wide significant variant associated with fT3 plasma levels, as well as a novel locus associated with TSH plasma levels. These findings are biologically relevant and enrich our knowledge about the genetic basis of pituitary-thyroid axis function.

Genetic Variants in the ST6GAL1 Gene Are Associated with Thyroglobulin Plasma Level in Healthy Individuals

Background: Thyroglobulin (Tg) is a 660 kDa iodoglycoprotein that serves as a scaffold for thyroid hormone synthesis. Although a twin study showed that variability of serum Tg levels has a substantial genetic basis, no genome-wide association study (GWAS) of serum/plasma Tg levels has been performed to date. The aim of this study was to identify genetic variants associated with plasma Tg levels among healthy individuals. Methods: A GWAS was conducted on two Croatian cohorts, and a combined analysis was performed. The analyses included 1094 individuals. A total of 7,597,379 variants, imputed using the 1000 Genomes reference panel, were analyzed for association. GWAS was performed under an additive model, controlling for age, sex, and relatedness within each data set. Combined analysis was conducted using the inverse-variance fixed-effects method. Results: Sixteen variants located on chromosome 3, within the ST6GAL1 gene, reached genome-wide significance. The lead SNP was rs4012172 ( \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$p = 1.29 \times {10^{ - 10}}$$ \end{document} ), which explained 3.19% of the variance in Tg levels. ST6GAL1 belongs to the sialyltransferase protein family, which has a fundamental role in the synthesis of specific sialylated structures on various glycoproteins, including Tg. It is known that only immature Tg (poorly sialylated or desialylated) can be transferred to the bloodstream. Conclusions: A highly biologically plausible locus was identified that could have a role in the regulation of plasma Tg levels in healthy individuals.

Genome-wide association analysis suggests novel loci underlying thyroid antibodies in Hashimoto's thyroiditis

Thyroid antibodies against thyroglobulin (TgAb) and thyroid peroxidase (TPOAb) are key markers of Hashimoto's thyroiditis (HT), the most common autoimmune thyroid disorder. Genetic determinants of thyroid antibodies are still poorly known, especially as they were not studied in patients with thyroid diseases. We performed the first genome-wide association analysis of thyroid antibodies in 430 HT patients that may be considered as population extremes for thyroid antibodies distribution. We detected two suggestively associated genetic variants with TgAb, rs6972286 close to ANKRD7 and LSM8 (P = 2.34 × 10^-7) and rs756763 inside CA10 (P = 6.05 × 10^-7), and one with TPOAb, rs12507813 positioned between TRIM61 and TRIM60 (P = 4.95 × 10^-7). Bivariate analysis resulted with three suggestively associated genetic variants that predispose to both antibodies: rs13190616 inside RP11-138J23.1 (P = 2.01 × 10^-6), rs561030786 close to DUBR (P = 7.33 × 10^-6) and rs12713034 inside FSHR (P = 7.66 × 10^-6). All identified genomic regions have a substantial literature record of involvement with female-related traits, immune-mediated diseases and personality traits that are all characterized by increased thyroid antibody levels. Our findings demonstrate the existence of genetic overlap between thyroid autoimmunity in HT and different non-thyroid diseases characterized by the presence of thyroid antibodies. We also suggest that genetic variants that regulate antibody levels may differ between HT patients and individuals with normal thyroid function.

The effect of multiple nutrients on plasma parathyroid hormone level in healthy individuals

Although the effect of isolated nutrients on plasma parathyroid hormone (PTH) is somewhat familiar, the effect of multiple nutrients on plasma PTH level has not yet been studied. The aim of this study was to identify groups of food items that are associated with the plasma PTH level in healthy individuals. This cross-sectional study enrolled 1180 healthy individuals from Croatia with plasma PTH levels inside the referent values. A food frequency questionnaire containing 58 food items was completed to evaluate the dietary intake. We used principal component analysis to reduce food items into dietary groups, followed by linear regression analysis to test the association between dietary groups and the level of PTH. The results indicate that different sorts of vegetables (p = .006), sausages, salami, mushrooms, eggs (p = .033), as well as white bread (p = .009) are associated with the increase, while bran bread (p = .009) is associated with the decreased plasma PTH level.

Effects of genetic variants on serum parathyroid hormone in hyperparathyroidism and end-stage renal disease patients: A systematic review and meta-analysis

BACKGROUND

Parathyroid hormone (PTH) is one of the principal regulators of calcium homeostasis, crucial for normal functioning of the kidneys, bones, heart, and nervous system. Different pathologic conditions can affect serum PTH level resulting in hyperparathyroidism or hypoparathyroidism. Our study assessed the association of previously reported polymorphisms with the level of PTH (expressed in pg/mL) among individuals with different pathologic conditions affecting PTH level.

METHODS

We searched Web of Science, MEDLINE, and Scopus to identify relevant articles published up to July 2017. The search yielded 6967 publications of which 44 fulfilled the inclusion criteria. We conducted meta-analyses for calcium-sensing receptor gene (CaSR) rs1801725 polymorphism in patients with primary hyperparathyroidism and vitamin D receptor gene (VDR) rs1544410 polymorphism in patients with end-stage renal disease (ESRD).

RESULTS

None of the polymorphisms were significantly associated with PTH levels in the overall population. In subgroup analysis by ethnicity for VDR rs1544410 gene polymorphism, we found significant differences under dominant model (SMD: -0.18 [-0.32, -0.05], P < .01) and AA versus GG comparison (SMD: -0.29 [-0.52, -0.06], P < .01) in Asian patients with ESRD, while nominally significant results (P < .05) were observed for AG versus GG and AA versus GG comparisons in European individuals with ESRD.

CONCLUSION

Scientific evidence of genetic association of serum PTH level among individuals with different pathologic conditions remains deficient and published results provide weak evidence. Further well-conducted studies on larger sample sets designed according to evidence-based principles are warranted to assure clinically applicable findings.

Genome-wide meta-analysis identifies novel loci associated with parathyroid hormone level

Background

Parathyroid hormone (PTH) is one of the principal regulators of calcium homeostasis. Although serum PTH level is mostly accounted by genetic factors, genetic background underlying PTH level is insufficiently known. Therefore, the aim of this study was to identify novel genetic variants associated with PTH levels.

Methods

We performed GWAS meta-analysis within two genetically isolated Croatian populations followed by replication analysis in a Croatian mainland population and we also combined results across all three analyzed populations. The analyses included 2596 individuals. A total of 7,411,206 variants, imputed using the 1000 Genomes reference panel, were analysed for the association. In addition, a sex-specific GWAS meta-analyses were performed.

Results

Polymorphisms with the lowest P-values were located on chromosome 4 approximately 84 kb of the 5′ of RASGEF1B gene. The most significant SNP was rs11099476 (P = 1.15 × 10⁻⁸). Sex-specific analysis identified genome-wide significant association of the variant rs77178854, located within DPP10 gene in females only (P = 2.21 × 10^− 9). There were no genome-wide significant findings in the meta-analysis of males.

Conclusions

We identified two biologically plausible novel loci associated with PTH levels, providing us with further insights into the genetics of this complex trait.

Electronic supplementary material

The online version of this article (10.1186/s10020-018-0018-5) contains supplementary material, which is available to authorized users.

Discovery and fine mapping of serum protein loci through transethnic meta-analysis

Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p < 5 × 10(-8)) for serum albumin (HPN-SCN1B, GCKR-FNDC4, SERPINF2-WDR81, TNFRSF11A-ZCCHC2, FRMD5-WDR76, and RPS11-FCGRT, in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein (TNFRS13B, 6q21.3, and ELL2, in up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

Historic, demographic, and genetic evidence for increased population frequencies of CCR5Delta32 mutation in Croatian Island isolates after lethal 15th century epidemics

AIM

To assess the frequency of 32 base pair deletion in CCR5 (CCR5Delta32), which has been shown to confer resistance to HIV infection in a homozygous form, in 10 isolated island communities of Dalmatia, Croatia, with different histories of exposure to epidemics during and since the medieval period.

METHODS

In 2002, DNA analysis of 100 randomly selected individuals from each of the 10 isolated communities of 5 Croatian islands (Susak, Rab, Vis, Lastovo, and Mljet) showed high levels of 3-generational endogamy, indicating limited gene flow. Five of the communities were decimated by epidemics of unknown cause between 1449-1456, while the other 5 villages remained unaffected. Genotyping of the CCR5 gene was performed using the polymerase chain reaction method with primers flanking the region containing 32-bp deletion.

RESULTS

The frequency of CCR5Delta32 in the 5 villages affected by the epidemic was 6.1-10.0%, and 1.0-3.8% in the 5 unaffected villages. The Delta32 mutation was found in 71 of 916 alleles among the individuals from the affected villages (7.5%), and in 24 of 968 alleles in unaffected villages (2.5%, chi(2)=27.3, P<10-6). A previous study in 303 random Croatian blood donors showed the frequency of the CCR5 Delta32 of 7.1% in the general population. The difference remained significant after correcting for population structure using both STRAT and STRUCTURE software and the genomic control test, to ensure results do not arise from the background genetic differences.

CONCLUSION

Our results and historical evidence, suggest that the mid-15th century epidemic could have acted as a selection pressure for the CCR5Delta32 mutation.

Genome-wide association study of biochemical traits in Korcula Island, Croatia

AIM

To identify genetic variants underlying biochemical traits--total cholesterol, low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, triglycerides, uric acid, albumin, and fibrinogen, in a genome-wide association study in an isolated population where rare variants of larger effect may be more easily identified.

METHODS

The study included 944 adult inhabitants of the island of Korcula, as a part of larger DNA-based genetic epidemiological study in 2007. Biochemical measurements were performed in a single laboratory with stringent internal and external quality control procedures. Examinees were genotyped using Human Hap370CNV chip by Illumina, with a genome-wide scan containing 346027 single nucleotide polymorphisms (SNP).

RESULTS

A total of 31 SNPs were associated with 7 investigated traits at the level of P<1.00 x 10(-5). Nine of SNPs implicated the role of SLC2A9 in uric acid regulation (P=4.10 x 10(-6)-2.58 x 10(-12)), as previously found in other populations. All 22 remaining associations fell into the P=1.00 x 10(-5)-1.00 x 10(-6) significance range. One of them replicated the association between cholesteryl ester transfer protein (CETP) and HDL, and 7 associations were more than 100 kilobases away from the closest known gene. Nearby SNPs, rs4767631 and rs10444502, in gene kinase suppressor of ras 2 (KSR2) on chromosome 12 were associated with LDL cholesterol levels, and rs10444502 in the same gene with total cholesterol levels. Similarly, rs2839619 in gene PBX/knotted 1 homeobox 1 (PKNOX1) on chromosome 21 was associated with total and LDL cholesterol levels. The remaining 9 findings implied possible associations between phosphatidylethanolamine N-methyltransferase (PEMT) gene and total cholesterol; USP46, RAP1GDS1, and ZCCHC16 genes and triglycerides; BCAT1 and SLC14A2 genes and albumin; and NR3C2, GRIK2, and PCSK2 genes and fibrinogen.

CONCLUSION

Although this study was underpowered for most of the reported associations to reach formal threshold of genome-wide significance under the assumption of independent multiple testing, replications of previous findings and consistency of association between the identified variants and more than one studied trait make such findings interesting for further functional follow-up studies. Changed allele frequencies in isolate population may contribute to identifying variants that would not be easily identified in much larger samples in outbred populations.

131 I-induced changes in rat thyroid gland function

Therapeutic doses of (131)I administered to thyrotoxic patients may cause thyroid failure. The present study used a rat model to determine thyroid function after the administration of different doses of (131)I (64-277 microCi). Thirty male Fisher rats in the experimental group and 30 in the control group (untreated) were followed for 6 months. The animals were 4 months old at the beginning of the experiment and were sacrificed at an age of 9 months. Hormone concentration was determined before (131)I administration (4-month-old animals) and three times following (131)I administration, when the animals were 7, 8, and 9 months old. The thyroid glands were removed and weighed, their volume was determined and histopathological examination was performed at the end of the experiment. Significant differences in serum triiodothyronine and thyroid-stimulating hormone concentration, measured at the age of 7, 8, and 9 months, were found in the experimental group. During aging of the animals, the concentration of thyroxin fell from 64.8 +/- 8.16 to 55.0 +/- 6.1 nM in the control group and from 69.4 +/- 6.9 to 25.4 +/- 3.2 nM in the experimental group. Thyroid gland volume and weight were significantly lower in the experimental than in the control group. Thyroid glands from the experimental group showed hyaline thickness of the blood vessel wall, necrotic follicles, a strong inflammatory reaction, and peeling of necrotic cells in the follicles. In conclusion, significant differences in hormone levels and histopathological findings indicated prolonged hypothyroidism after (131)I administration to rats, which was not (131)I dose dependent.