- Research article
- Open Access
A case–control study of selenoprotein genes polymorphisms and autoimmune thyroid diseases in a Chinese population
BMC Medical Genetics volume 18, Article number: 54 (2017)
Selenium is an essential trace and there is a high selenium concentration in the thyroid gland. Selenium deficiency may impair the thyroid function. The aim of this study was to investigate the association between three selenoprotein genes polymorphisms and autoimmune thyroid diseases.
We genotyped six single-nucleotide polymorphisms (SNPs), rs6865453 in selenoprotein P gene (SELENOP), rs713041 rs2074451 rs3746165 in glutathione peroxidase 4 gene (GPX4) and rs28665122 and rs7178239 in selenoprotein S gene (SELENOS) by MassARRAY system using the chip-based matrix-assisted laser desorption ionization time-of-flight mass spectrometry technology in 1060 patients with autoimmune thyroid diseases and 938 healthy controls.
Major alleles in rs6865453 of SELENOP, rs713041, rs2074451, rs3746165 of GPX4 decreased while the major allele C in rs28665122 of SELENOS increased in AITD patients than in the control. The allele C and genotype CC in rs7178239 of SELENOS showed different trend in GD and HT patients when compared with the control. All the distribution difference showed nonsignificant. Analysis according to clinical features including ophthalmopathy, hypothyroidism and family history came out to be negative either.
Our findings suggest non-association between three selenoprotein genes and AITD, conflicting to the positive result in another population. Different selenium nutrition status in different populations may contribute to conflicting results, the contribution of genetic variants in AITD mechanism may be another reason.
Selenium is an essential trace element and its distribution in nature varies widely. The soil of central and southern European countries is generally poor in selenium content, resulting in mild regional selenopenia, whereas in most North American regions it is abundant. Severe selenopenia was observed in large areas of central Asia. Selenium is present in the active site of selenoproteins, which are crucial to immune function, inflammation, redox processes and important for thyroid hormone synthesis and metabolism. There is a high selenium concentration in the thyroid gland, and selenium deficiency may impair the thyroid function. Reactive oxygen species (ROS), formed from thyroid hormone metabolism, may cause damage to different cellular structures and trigger tissue inflammation. Selenoproteins, which are powerful antioxidant enzymes, mitigate the effects of oxidative stress by elimination of ROS. In addition, selenoproteins play a vital role in the regulation of human immune system and Se deficiency is accompanied by dysregulation of both cell-mediated immunity and B cell function.
Autoimmune thyroid diseases (AITD), mainly including Graves’ disease (GD) and Hashimoto’s thyroiditis (HT), are the most frequent autoimmune disorders. GD is characterized clinically by hyperthyroidism, diffuse goiter and the presence of thyrotropin receptor antibodies (TRAb), while HT is characterized by apoptosis of thyrocytes leading to hypothyroidism and the presence of thyroid peroxidase antibodies (TPOAb) or antibodies against thyroglobulin (TGAb). Autoimmune thyroid disease (AITD) is a multifactorial disease in which autoimmunity against thyroid antigens develops against a particular genetic background facilitated by exposure to environmental factors. Environmental factors include smoking, vitamin D, iodine, selenium and so on. GD and HT also share some but not all susceptible genes.
The value of selenium supplementation in autoimmune thyroid disorders had been investigated and most studies confirmed the beneficial effect of selenium supplementation in Hashimoto’s and Graves’s diseases, though some studies showed conflicting results as well [1,2,3,4,5,6,7,8,9,10,11]. Recently, selenium was proved to be effective in mild inflammatory thyroid associated orbitopathy (TAO) too . Different selenium nutrition state in different studies may partly explain different reactivity to the selenium supplementation, different genetics factors in different populations may be another reason. Few studies were done about the association of selenoprotein genetic variations in AITD.
The role for selenium supplementation in the treatment of autoimmune thyroid diseases is still unresolved. Clinic studies provide evidence for a preventive effect of selenium on autoimmune thyroid disease, which may be mediated by the antioxidative and anti-inflammatory properties of selenoproteins, we therefore investigated whether genetic variants in selenoproteins are associated with GD and HT. Among 25 selenoproteins known to date, various functions are characterized for each of them. These include Se transport (selenoprotein P), antioxidant/redox properties (glutathione peroxidases (GPxs), thioredoxin reductases and selenoprotein P) and anti-inflammatory properties (selenoprotein S and GPx4). We chosed to study SNPs rs6865453 in selenoprotein P gene (SELENOP), rs713041 rs2074451 rs3746165 in glutathione peroxidase 4 gene (GPX4), rs28665122 and rs7178239 in selenoprotein S gene (SELENOS) in a Han Chinese population.
Our study population (n = 1998) consisted of 1060 Chinese AITD patients including 701 GD patients, 359 HT patients and 938 healthy unrelated controls. All AITD patients in the present case–control study were recruited from the Department of Endocrinology, Jinshan Hospital Affiliated to Fudan University. The diagnostic criteria for GD were mainly determined by clinical manifestation, laboratory biochemical proof of hyperthyroidism, the presence of diffuse goiter and the positive circulating TSH receptor antibody (TRAb). HT was defined on the basis of enlarged thyroid and the high level of either TPOAb or TgAb, with or without clinical and biochemical hypothyroidism. All healthy controls were unrelated individuals with neither thyroid diseases nor other autoimmune diseases, recruited from the Medical Examination Center of the same hospital. All the subjects, including AITD patients and controls, were Han Chinese and signed the informed consent. Both patients and controls have similar diet habits and individuals with distinct diet style such as vegetarians were excluded. The research project was approved by the Ethics Committee of Jinshan Hospital Affiliated to Fudan University.
Peripheral venous blood of 2 ml from the subjects was collected in an EDTA tube. The genomic DNA was extracted by salting-out method, using RelaxGene Blood DNA System (TIANGEN BIOTECH, Beijing, China), according to the manufacturers’ protocol. Genotyping of SNPs (rs6865453 in SELENOP, rs713041, rs2074451 and rs3746165 in GPX4, rs28665122 and rs7178239 in SELENOS) was performed by MassARRAY system (Sequenom, San Diego, CA, USA) using the chip-based matrix-assisted laser desorption ionization time-of-flight mass spectrometry technology. This is a well-selected platform by many researchers during studies including SNP genotyping. It has high specificity and sensitivity. Primers were obtained from Sangon Biotech (Shanghai, China). Briefly, multiplex reaction was designed using Assay Designer software version 3.0 (Sequenom) and was processed following standard protocols for iPLEX chemistry. The reaction products were then cleaned and dispensed onto a Spectro-CHIP bioarray. The chip was scanned using MassARRAY workstation version 3.3, and the resulting spectra were analyzed using the Sequenom TYPER software.
We evaluated the Hardy-Weinberg equilibrium of SNPs within our population using the HWE program . Allele and genotype frequencies between cases and controls were computed by chi-square test or Fisher’s exact test. Differences between groups were determined by the odds ratio (OR) and 95% confidence interval (95%CI). These statistical analyses were performed using the software GraphPad Prism 6. Rectified analysis between cases and controls for ages were conducted by Binary Logistic regression analysis using the software SPSS22.0. A p value less than 0.05 was considered significant.
Clinical data analysis
As shown in Table 1, our study investigated 1060 AITD patients, who comprised 701 GD (214 men and 487 women) and 359 HT patients (44 men and 315 women). In GD patients, the average onset age was 33.86 ± 14.36, 145 individuals had family history and 154 had ophthalmopathy. In HT patients, the average onset age was 32.63 ± 13.44, 72 individuals had family history and 6 had ophthalmopathy. In addition, 184 patients suffered from hypothyroidism while other 175 were of euthyroidism in the HT group.
Allele and genotype results
The distribution of genotype frequencies for each SNP (rs6865453/rs713041/rs2074451/rs3746165/rs28665122/rs7178239) was in Hardy–Weinberg equilibrium in both the patients and the controls (Table 2). The allele distributions of the 6 SNPs of three genes, their genotype and case–control association analysis were shown in (Table 3). The major allele A and AA genotype in SELENOP rs6865453, the major allele C and CC genotype in rs713041, the major allele G and GG genotype in rs2074451, and the major allele C and CC genotype in rs3746165 of GPX4 gene decreased in both GD and HT patients than those in the control. The major allele C and genotype CC in SELENOS rs7178239 showed different trend in GD and HT patients when compared to the control. While the major allele C and CC genotype in rs28665122 of SELENOS increased in both GD and HT patients than the control, different to the result for another study, in which these allele and genotype decreased significantly in their HT patients . Our association analysis showed negative and the distribution difference showed nonsignificant.
Genotype, allele distribution and clinical phenotype association
Etiology for Graves’ ophthalmopathy remains unclear, and genetic factors contribute to it too. Here we compared the genotype, allele distribution of Graves’ ophthalmopathy patients or non-ophthalmopathy ones to the control and also did the association analysis between ophthalmopathy and non-ophthalmopathy patients, but no association was found (Table 4). When we did similar analysis between HT patients with hypothyroidism and the euthyroidism HT patients, no significant difference was found either (Table 5). We did the association analysis according to the patients’ family history too. The statistic results were negative in GD patients grouped on the basis of their family history (Table 6). When HT patients were divided into family history-positive and negative groups, the patients number was too small to do effective statistical analysis. So this part of analysis was ignored.
Selenium is incorporated into selenoproteins implicated in antioxidative defence, thyroid hormone metabolism and immunoregulation. Selenoproteins are important for the thyroid function and the thyroid health. Selenium metabolism is controlled by hepatocytes synthesizing and secreting the Se transporter selenoprotein P (SePP). SePP, encoded by SELENOP gene, is a liver-derived plasma protein and the major selenium containing protein in serum. Serum levels of SePP in humans are considered a biomarker of selenium status in the entire body. Its promoter activity has been shown to be inhibited by cytokines including interleukin 1β, tumor necrosis factor α, interferon γ and transforming growth factor β1 [15, 16], suggesting a role of SELENOP in the process of inflammation. SePP is able to bind to the endothelium and be recruited to the site of the inflammation. It also has the ability to bind epithelial cells and displays phospholipid hydroperoxide thiol peroxidase activity . Overexpression of SELENOP suppressed hydrogenperoxide-induced activation of JNK and p38 whereas the gene’s knocked-down led to increased production of lipid hydroperoxide . Studies showed that genetic variants in SELENOP were associated with the risk of breast cancer, prostate cancer, advanced colorectal adenoma and colorectal cancer as well as first phase insulin response and the occurrence of abdominal aortic aneurysm [19,20,21,22,23,24]. No studies about their roles in AITD were reported. We detected the association of rs6865453 in SELENOP with AITD, but it was negative.
Phospholipid hydroperoxide glutathione peroxidase (PHGPx, GPx4) is an intracellular selenoprotein, a major antioxidant enzyme, which plays unique roles in the protection of cells against oxidative stress by catalysing reduction of lipid hydroperoxides. GPX4 may also play a role in regulation of leukotriene biosynthesis and cytokine signaling pathways [25, 26]. GPx4 rs713041 polymorphisms could modulate the synthesis and expression of GPx4 by altering the selenocysteine insertion [27, 28]. GPx4 rs713041 were found to be associated with breast cancer mortality, the risk of colorectal cancer and cerebral stroke [23, 29, 30]. GPx4 rs2074451 variant was also found to be associated with breast cancer survival . And rs3746165 variants in GPx4 may be associated with risk of lethal prostate cancer . No relevance of these SNPs to AITD was reported. We detected these three GPx4 SNPs in this case–control study, but no association with AITD was found.
Selenoprotein S (SELENOS, gene aliases: SELS, VIMP, TANIS) has been classified as a new endoplasmic reticulum (ER) membrane protein that moves misfolded proteins from the ER to the cytosol and prevents stress responses that activate the inflammatory cascade . It is a membrane-bound selenoprotein and located in ER and other plasma membrane. It protects cells from oxidative damage and it impacts the immune and inflammatory signal pathways. It is expected that genetic polymorphisms affecting SELENOS gene transcription and subsequent SELENOS expression levels might contribute to the development and progression of inflammatory disorders. A previous study discovered that a functional single nucleotide polymorphism (SNP) in the SELENOS promoter region (G-105A, rs28665122) impairs SELENOS expression and amplifies the production of inflammatory cytokines such as tumor necrosis factor-a (TNF-a), interleukin-6 (IL-6) and interleukin-1b (IL-1b) . No significant differences were observed in the genotypic distribution of the SELENOS polymorphisms analyzed in Spanish cohorts between patients suffering from autoimmune diseases with a clear inflammatory component, like T1D, RA or IBD and controls . While the SELENOS rs28665122 polymorphism had been found to be associated with the risk of spontaneous preterm birth in a Chinese population , the existence of association between this genetic variation and HT risk was also found in a Portuguese population with a higher allele A proportion in the patient group . Here we detected the link between rs28665122 and rs7178239 polymorphisms in SELENOS gene and AITD in a Chinese population. Unfortunately, the association between AITD risk and these two SELENOS SNPs in our study was negative. Opposing to the result from the Portuguese population, the allele C and CC genotype in SELENOS rs28665122 increased in HT patients in our study. Inconsistent result suggests that different genetic variants may be involved in the etiology of AITD in different populations. Different pathogenic mechanisms may underlie autoimmune thyroid disease in different populations. This may partly explain the conflicting effects of selenium administration on AITD antibodies in studies from different populations. Different selenium nutrition status of the population from different geographical location may be another factor contributing to this inconsistent result. Clearly, further in-depth studies and evaluation are required concerning the mechanism of selenium action and the etiology of AITD.
Our findings suggested non-association between six SNPs of three selenprotein genes and AITD, non-association between SNPs and clinical genotypes. Our results were conflicted to the result from the Portuguese population study. In that study, the higher allele A in SELENOS rs28665122 had been found to be associated with the HT risk . Different selenium nutrition state in different populations may contribute to conflicting results, the contribution of genetic variants in AITD mechanism may be another reason.
Autoimmune thyroid disease
Glutathione peroxidase 4 gene
Phospholipid hydroperoxide glutathione peroxidase
Inflammatory bowel disease
Reactive oxygen species
Selenoprotein P gene
Selenoprotein S gene
Se transporter selenoprotein P
Type 1 diabetes
Thyroid associated orbitopathy
Antibodies against thyroglobulin
Tumor necrosis factor-a
Thyroid peroxidase antibodies
Thyrotropin receptor antibodies
Gartner R, Gasnier BC, Dietrich JW, Krebs B, Angstwurm MW. Selenium supplementation in patients with autoimmune thyroiditis decreases thyroid peroxidase antibodies concentrations. J Clin Endocrinol Metab. 2002;87(4):1687–91.
Bonfig W, Gartner R, Schmidt H. Selenium supplementation does not decrease thyroid peroxidase antibody concentration in children and adolescents with autoimmune thyroiditis. TheScientificWorldJOURNAL. 2010;10:990–6.
Zagrodzki P, Ratajczak R. Selenium supplementation in autoimmune thyroiditis female patient--effects on thyroid and ovarian functions (case study). Biol Trace Elem Res. 2008;126(1–3):76–82.
Moncayo R, Moncayo H, Kapelari K. Nutritional treatment of incipient thyroid autoimmune disease. Influence of selenium supplementation on thyroid function and morphology in children and young adults. Clin Nutr. 2005;24(4):530–1.
Vrca VB, Skreb F, Cepelak I, Romic Z, Mayer L. Supplementation with antioxidants in the treatment of Graves’ disease; the effect on glutathione peroxidase activity and concentration of selenium. Clin Chim Acta. 2004;341(1–2):55–63.
Calissendorff J, Mikulski E, Larsen EH, Moller M. A prospective investigation of Graves’ disease and selenium: thyroid hormones, auto-antibodies and self-rated symptoms. Eur Thyroid J. 2015;4(2):93–8.
Duntas LH, Mantzou E, Koutras DA. Effects of a six month treatment with selenomethionine in patients with autoimmune thyroiditis. Eur J Endocrinol. 2003;148(4):389–93.
Turker O, Kumanlioglu K, Karapolat I, Dogan I. Selenium treatment in autoimmune thyroiditis: 9-month follow-up with variable doses. J Endocrinol. 2006;190(1):151–6.
Nacamulli D, Mian C, Petricca D, Lazzarotto F, Barollo S, Pozza D, Masiero S, Faggian D, Plebani M, Girelli ME, et al. Influence of physiological dietary selenium supplementation on the natural course of autoimmune thyroiditis. Clin Endocrinol. 2010;73(4):535–9.
Bacic Vrca V, Skreb F, Cepelak I, Mayer L. Supplementation with antioxidants in the treatment of Graves’ disease: the effect on the extracellular antioxidative parameters. Acta Pharma. 2004;54(2):79–89.
Karanikas G, Schuetz M, Kontur S, Duan H, Kommata S, Schoen R, Antoni A, Kletter K, Dudczak R, Willheim M. No immunological benefit of selenium in consecutive patients with autoimmune thyroiditis. Thyroid. 2008;18(1):7–12.
Marcocci C, Kahaly GJ, Krassas GE, Bartalena L, Prummel M, Stahl M, Altea MA, Nardi M, Pitz S, Boboridis K, et al. Selenium and the course of mild Graves’ orbitopathy. N Engl J Med. 2011;364(20):1920–31.
Hardy-Weinberg equilibrium. http://ihg.gsf.de/cgi-bin/hw/hwa1.pl. Accessed 28 Apr 2017.
Santos LR, Duraes C, Mendes A, Prazeres H, Alvelos MI, Moreira CS, Canedo P, Esteves C, Neves C, Carvalho D, et al. A polymorphism in the promoter region of the selenoprotein S gene (SEPS1) contributes to Hashimoto’s thyroiditis susceptibility. J Clin Endocrinol Metab. 2014;99(4):E719–723.
Dreher I, Jakobs TC, Kohrle J. Cloning and characterization of the human selenoprotein P promoter. Response of selenoprotein P expression to cytokines in liver cells. J Biol Chem. 1997;272(46):29364–71.
Mostert V, Dreher I, Kohrle J, Wolff S, Abel J. Modulation of selenoprotein P expression by TGF-beta (1) is mediated by Smad proteins. BioFactors. 2001;14(1–4):135–42.
Burk RF, Hill KE. Selenoprotein P: an extracellular protein with unique physical characteristics and a role in selenium homeostasis. Annu Rev Nutr. 2005;25:215–35.
Kabuyama Y, Oshima K, Kitamura T, Homma M, Yamaki J, Munakata M, Homma Y. Involvement of selenoprotein P in the regulation of redox balance and myofibroblast viability in idiopathic pulmonary fibrosis. Genes Cells. 2007;12(11):1235–44.
Peters U, Chatterjee N, Hayes RB, Schoen RE, Wang Y, Chanock SJ, Foster CB. Variation in the selenoenzyme genes and risk of advanced distal colorectal adenoma. Cancer Epidemiol Biomark Prev. 2008;17(5):1144–54.
Strauss E, Oszkinis G, Staniszewski R. SEPP1 gene variants and abdominal aortic aneurysm: gene association in relation to metabolic risk factors and peripheral arterial disease coexistence. Sci Rep. 2014;4:7061.
Hellwege JN, Palmer ND, Ziegler JT, Langefeld CD, Lorenzo C, Norris JM, Takamura T, Bowden DW. Genetic variants in selenoprotein P plasma 1 gene (SEPP1) are associated with fasting insulin and first phase insulin response in Hispanics. Gene. 2014;534(1):33–9.
Steinbrecher A, Meplan C, Hesketh J, Schomburg L, Endermann T, Jansen E, Akesson B, Rohrmann S, Linseisen J. Effects of selenium status and polymorphisms in selenoprotein genes on prostate cancer risk in a prospective study of European men. Cancer Epidemiol Biomarkers Prev. 2010;19(11):2958–68.
Meplan C, Hughes DJ, Pardini B, Naccarati A, Soucek P, Vodickova L, Hlavata I, Vrana D, Vodicka P, Hesketh JE. Genetic variants in selenoprotein genes increase risk of colorectal cancer. Carcinogenesis. 2010;31(6):1074–9.
Pellatt AJ, Wolff RK, John EM, Torres-Mejia G, Hines LM, Baumgartner KB, Giuliano AR, Lundgreen A, Slattery ML. SEPP1 influences breast cancer risk among women with greater native american ancestry: the breast cancer health disparities study. PLoS One. 2013;8(11):e80554.
Brigelius-Flohe R. Tissue-specific functions of individual glutathione peroxidases. Free Radic Biol Med. 1999;27(9–10):951–65.
Villette S, Kyle JA, Brown KM, Pickard K, Milne JS, Nicol F, Arthur JR, Hesketh JE. A novel single nucleotide polymorphism in the 3′ untranslated region of human glutathione peroxidase 4 influences lipoxygenase metabolism. Blood Cells Mol Dis. 2002;29(2):174–8.
Meplan C, Crosley LK, Nicol F, Horgan GW, Mathers JC, Arthur JR, Hesketh JE. Functional effects of a common single-nucleotide polymorphism (GPX4c718t) in the glutathione peroxidase 4 gene: interaction with sex. Am J Clin Nutr. 2008;87(4):1019–27.
Gautrey H, Nicol F, Sneddon AA, Hall J, Hesketh J. A T/C polymorphism in the GPX4 3′UTR affects the selenoprotein expression pattern and cell viability in transfected Caco-2 cells. Biochim Biophys Acta. 2011;1810(6):584–91.
Udler M, Maia AT, Cebrian A, Brown C, Greenberg D, Shah M, Caldas C, Dunning A, Easton D, Ponder B, et al. Common germline genetic variation in antioxidant defense genes and survival after diagnosis of breast cancer. J Clin Oncol Off J Am Soc Clin Oncol. 2007;25(21):3015–23.
Polonikov AV, Vialykh EK, Churnosov MI, Illig T, Freidin MB, Vasil’eva OV, Bushueva OY, Ryzhaeva VN, Bulgakova IV, Solodilova MA. The C718T polymorphism in the 3′-untranslated region of glutathione peroxidase-4 gene is a predictor of cerebral stroke in patients with essential hypertension. Hypertens Res. 2012;35(5):507–12.
Van Blarigan EL, Ma J, Kenfield SA, Stampfer MJ, Sesso HD, Giovannucci EL, Witte JS, Erdman Jr JW, Chan JM, Penney KL. Plasma antioxidants, genetic variation in SOD2, CAT, GPX1, GPX4, and prostate cancer survival. Cancer Epidemiol Biomarkers Prev. 2014;23(6):1037–46.
Ye Y, Shibata Y, Yun C, Ron D, Rapoport TA. A membrane protein complex mediates retro-translocation from the ER lumen into the cytosol. Nature. 2004;429(6994):841–7.
Curran JE, Jowett JB, Elliott KS, Gao Y, Gluschenko K, Wang J, Abel Azim DM, Cai G, Mahaney MC, Comuzzie AG, et al. Genetic variation in selenoprotein S influences inflammatory response. Nat Genet. 2005;37(11):1234–41.
Suzuki Y, Schmitt MJ. Redox diversity in ERAD-mediated protein retrotranslocation from the endoplasmic reticulum: a complex puzzle. Biol Chem. 2015;396(5):539–54.
Wang Y, Yang X, Zheng Y, Wu ZH, Zhang XA, Li QP, He XY, Wang CZ, Feng ZC. The SEPS1 G-105A polymorphism is associated with risk of spontaneous preterm birth in a Chinese population. PLoS One. 2013;8(6):e65657.
Acknowledgements to Jian Zhang from the department of Clinical Laboratory in Jinshan Hospital affiliated to Fudan University for her help to collect patients’ samples. Acknowledgements to the fund from Shanghai Jinshan District Health and Family Planning Commission (JSKJ-KTMS-2013-02).
Shanghai Jinshan District Health and Family Planning Commission (JSKJ-KTMS-2013-02). The funding body played no role in the design of the study and the collection, analysis and interpretation of data and the manuscript writing.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
LX carried out the molecular genetic studies, participated in the sequence alignment and drafted the manuscript. JY and RS carried out the molecular genetic studies. NY and WJ participated in the sequence alignment. DL and LS participated in the design of the study and performed the statistical analysis. JZ conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All the subjects were Han Chinese and signed the informed consent. The research project was approved by the Ethics Committee of Jinshan Hospital Affiliated to Fudan University.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xiao, L., Yuan, J., Yao, Q. et al. A case–control study of selenoprotein genes polymorphisms and autoimmune thyroid diseases in a Chinese population. BMC Med Genet 18, 54 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-017-0415-6
- Selenoprotein P gene (SELENOP)
- Glutathione peroxidase 4 gene (GPX4)
- Selenoprotein S gene (SELENOS)
- Single-nucleotide polymorphisms (SNPs)
- Autoimmune thyroid disease (AITD)