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.

Thyroid Function and Metabolic Syndrome: A Two-Sample Bidirectional Mendelian Randomization Study

CONTEXT

Thyroid function has been associated with metabolic syndrome (MetS) in a number of observational studies but the direction of effects and the exact causal mechanism of this relationship is still unknown.

OBJECTIVE

To examine genetically predicted effects of thyroid function on MetS risk and its components, and vice versa, using large-scale summary genetic association data.

METHODS

We performed a two-sample bidirectional Mendelian randomization (MR) study using summary statistics from the most comprehensive genome-wide association studies (GWAS) of thyroid-stimulating hormone (TSH, n = 119 715), free thyroxine (fT4, n = 49 269), MetS (n = 291 107), and components of MetS: waist circumference (n = 462 166), fasting blood glucose (n = 281 416), hypertension (n = 463 010), triglycerides (TG, n = 441 016) and high-density lipoprotein cholesterol (HDL-C, n = 403 943). We chose the multiplicative random effects inverse variance weighted (IVW) method as the main analysis. Sensitivity analysis included weighted median and mode analysis, as well as MR-Egger and Causal Analysis Using Summary Effect estimates (CAUSE).

RESULTS

Our results suggest that higher fT4 levels lower the risk of developing MetS (OR = 0.96, P = .037). Genetically predicted fT4 was also positively associated with HDL-C (β = 0.02, P = .008), while genetically predicted TSH was positively associated with TG (β = 0.01, P = .044). These effects were consistent across different MR analyses and confirmed with the CAUSE analysis. In the reverse direction MR analysis, genetically predicted HDL-C was negatively associated with TSH (β = -0.03, P = .046) in the main IVW analysis.

CONCLUSION

Our study suggests that variations in normal-range thyroid function are causally associated with the diagnosis of MetS and with lipid profile, while in the reverse direction, HDL-C has a plausible causal effect on reference-range TSH levels.

Exploring Neurophysiological Mechanisms and Treatment Efficacies in Laryngeal Dystonia: A Transcranial Magnetic Stimulation Approach

Laryngeal dystonia (LD), known or termed as spasmodic dysphonia, is a rare movement disorder with an unknown cause affecting the intrinsic laryngeal muscles. Neurophysiological studies point to perturbed inhibitory processes, while conventional genetic studies reveal fragments of genetic architecture in LD. The study's aims are to (1) describe transcranial magnetic stimulation (TMS) methodology for studying the functional integrity of the corticospinal tract by stimulating the primary motor cortex (M1) for laryngeal muscle representation and recording motor evoked potentials (MEPs) from laryngeal muscles; (2) evaluate the results of TMS studies investigating the cortical silent period (cSP) in LD; and (3) present the standard treatments of LD, as well as the results of new theoretical views and treatment approaches like repetitive TMS and laryngeal vibration over the laryngeal muscles as the recent research attempts in treatment of LD. Neurophysiological findings point to a shortened duration of cSP in adductor LD and altered cSP duration in abductor LD individuals. Future TMS studies could further investigate the role of cSP in relation to standard laryngological measures and treatment options. A better understanding of the neurophysiological mechanisms might give new perspectives for the treatment of LD.

Role of ST6GAL1 in Thyroid Cancers: Insights from Tissue Analysis and Genomic Datasets

Thyroid cancer is the predominant endocrine-related malignancy. ST6 β-galactoside α2,6-sialyltransferase 1 (ST6GAL1) has been studied in various types of cancers; however, the expression and function of ST6GAL1 in thyroid cancer has not been investigated so far. Previously, we conducted two genome-wide association studies and have identified the association of the ST6GAL1 gene with plasma thyroglobulin (Tg) levels. Since Tg levels are altered in thyroid pathologies, in the current study, we wanted to evaluate the expression of ST6GAL1 in thyroid cancer tissues. We performed an immunohistochemical analysis using human thyroid tissue from 89 patients and analyzed ST6GAL1 protein expression in papillary thyroid cancer (including follicular variant and microcarcinoma) and follicular thyroid cancer in comparison to normal thyroid tissue. Additionally, ST6GAL1 mRNA levels from The Cancer Genome Atlas (TCGA, n = 572) and the Genotype-Tissue Expression (GTEx) project (n = 279) were examined. The immunohistochemical analysis revealed higher ST6GAL1 protein expression in all thyroid tumors compared to normal thyroid tissue. TCGA data revealed increased ST6GAL1 mRNA levels in both primary and metastatic tumors versus controls. Notably, the follicular variant of papillary thyroid cancer exhibited significantly higher ST6GAL1 mRNA levels than classic papillary thyroid cancer. High ST6GAL1 mRNA levels significantly correlated with lymph node metastasis status, clinical stage, and reduced survival rate. ST6GAL1 emerges as a potential cancer-associated glycosyltransferase in thyroid malignancies, offering valuable insights into its diagnostic and prognostic significance.

Vitamin D and the Thyroid: A Critical Review of the Current Evidence

Vitamin D is necessary for the normal functioning of many organs, including the thyroid gland. It is, therefore, not surprising that vitamin D deficiency is considered a risk factor for the development of many thyroid disorders, including autoimmune thyroid diseases and thyroid cancer. However, the interaction between vitamin D and thyroid function is still not fully understood. This review discusses studies involving human subjects that (1) compared vitamin D status (primarily determined by serum calcidiol (25-hydroxyvitamin D [25(OH)D]) levels) with thyroid function assessed by thyroid stimulating hormone (TSH), thyroid hormones, and anti-thyroid antibody levels; and (2) evaluated the effect of vitamin D supplementation on thyroid function. Due to the many inconsistencies in the results between the studies, it is still difficult to draw a definite conclusion on how vitamin D status affects thyroid function. Studies in healthy participants observed either a negative correlation or no association between TSH and 25(OH)D levels, while the results for thyroid hormones showed high variability. Many studies have observed a negative association between anti-thyroid antibodies and 25(OH)D levels, but equally many studies have failed to observe such an association. Regarding the studies that examined the effect of vitamin D supplementation on thyroid function, almost all observed a decrease in anti-thyroid antibody levels after vitamin D supplementation. Factors that could contribute to the high variability between the studies are the use of different assays for the measurement of serum 25(OH)D levels and the confounding effects of sex, age, body-mass index, dietary habits, smoking, and the time of year when the samples were collected. In conclusion, additional studies with larger numbers of participants are needed to fully understand the effect of vitamin D on thyroid function.

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.

Serotonin Receptor Gene Polymorphisms Are Associated with Cerebrospinal Fluid, Genetic, and Neuropsychological Biomarkers of Alzheimer's Disease

A decrease in serotonergic transmission throughout the brain is among the earliest pathological changes in Alzheimer's disease (AD). Serotonergic receptors are also affected in AD. Polymorphisms in genes of serotonin (5HT) receptors have been mostly associated with behavioral and psychological symptoms of dementia (BPSD). In this study, we examined if AD patients carrying different genotypes in 5HTR1B rs13212041, 5HTR2A rs6313 (T102C), 5HTR2C rs3813929 (-759C/T), and 5HTR6 rs1805054 (C267T) polymorphisms have a higher risk of faster disease progression (assessed by neuropsychological testing), are more prone to develop AD-related pathology (reflected by levels of cerebrospinal fluid [CSF] AD biomarkers), or have an association with an apolipoprotein E (APOE) haplotype. This study included 115 patients with AD, 53 patients with mild cognitive impairment (MCI), and 2701 healthy controls. AD biomarkers were determined in the CSF of AD and MCI patients using enzyme-linked immunosorbent assays (ELISA), while polymorphisms were determined using either TaqMan SNP Genotyping Assays or Illumina genotyping platforms. We detected a significant decrease in the CSF amyloid β_1-42 (Aβ_1-42) and an increase in p-tau₁₈₁/Aβ_1-42 ratio in carriers of the T allele in the 5HTR2C rs3813929 (-759C/T) polymorphism. A significantly higher number of APOE ε4 allele carriers was observed among individuals carrying a TT genotype within the 5HTR2A T102C polymorphism, a C allele within the 5HTR1B rs13212041 polymorphism, and a T allele within the 5HTR6 rs1805054 (C267T) polymorphism. Additionally, individuals carrying the C allele within the 5HTR1B rs13212041 polymorphism were significantly more represented among AD patients and had poorer performances on the Rey-Osterrieth test. Carriers of the T allele within the 5HTR6 rs1805054 had poorer performances on the MMSE and ADAS-Cog. As all four analyzed polymorphisms of serotonin receptor genes showed an association with either genetic, CSF, or neuropsychological biomarkers of AD, they deserve further investigation as potential early genetic biomarkers of AD.

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.

Presence or severity of Hashimoto's thyroiditis does not influence basal calcitonin levels: observations from CROHT biobank

PURPOSE

The influence of Hashimoto's thyroiditis (HT) on calcitonin (Ct) production is unresolved question. The aim of this study was to explore if basal Ct levels are influenced by the presence/severity of HT or correlated with clinical phenotypes of HT patients.

METHODS

We included 467 HT patients and 184 control participants, from Croatian Biobank of HT patients (CROHT), in this retrospective study. Calcitonin levels between HT patients and controls were compared using Mann-Whitney test. Ct levels between two subgroups of HT patients, divided by intake of levothyroxine (LT4) therapy, were additionally tested to take into account the illness severity. Spearman rank correlation test was used to analyze correlations between Ct levels and 14 relevant phenotypes.

RESULTS

We have not detected significant differences in median Ct levels between HT patients and controls (2.2 vs 2.35 pg/mL, respectively, P = 0.717) nor in-between two subgroups of HT patients (P = 0.347). We have not detected statistically significant correlations between Ct levels and clinical phenotypes, although we identified three weak nominal correlations: negative correlation of Ct with TgAb in all HT patients (r = - 0.1, P = 0.04); negative correlation of Ct with age in subgroup of HT patients without LT4 therapy (r = - 0.13, P = 0.04); positive correlation of Ct with BSA in subgroup of HT patients on LT4 therapy (r = 0.16, P = 0.042).

CONCLUSION

Our results suggest that HT patients of all disease stages preserve Ct production as healthy individuals and there is no need for Ct measurements in the absence of a nodule. Additional confirmation and clarification of observed nominal correlations are needed due to potential clinical relevance of TgAb and age-dependent Ct decrease in HT women.

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.

Environmental Factors Affecting Thyroid-Stimulating Hormone and Thyroid Hormone Levels

Thyroid hormones are necessary for the normal functioning of physiological systems. Therefore, knowledge of any factor (whether genetic, environmental or intrinsic) that alters the levels of thyroid-stimulating hormone (TSH) and thyroid hormones is crucial. Genetic factors contribute up to 65% of interindividual variations in TSH and thyroid hormone levels, but many environmental factors can also affect thyroid function. This review discusses studies that have analyzed the impact of environmental factors on TSH and thyroid hormone levels in healthy adults. We included lifestyle factors (smoking, alcohol consumption, diet and exercise) and pollutants (chemicals and heavy metals). Many inconsistencies in the results have been observed between studies, making it difficult to draw a general conclusion about how a particular environmental factor influences TSH and thyroid hormone levels. However, lifestyle factors that showed the clearest association with TSH and thyroid hormones were smoking, body mass index (BMI) and iodine (micronutrient taken from the diet). Smoking mainly led to a decrease in TSH levels and an increase in triiodothyronine (T3) and thyroxine (T4) levels, while BMI levels were positively correlated with TSH and free T3 levels. Excess iodine led to an increase in TSH levels and a decrease in thyroid hormone levels. Among the pollutants analyzed, most studies observed a decrease in thyroid hormone levels after exposure to perchlorate. Future studies should continue to analyze the impact of environmental factors on thyroid function as they could contribute to understanding the complex background of gene-environment interactions underlying the pathology of thyroid diseases.

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 for Hashimoto's thyroiditis

PURPOSE

Hashimoto's thyroiditis (HT) is the most common form of autoimmune thyroid diseases. Current knowledge of HT genetics is limited, and not a single genome-wide association study (GWAS) focusing exclusively on HT has been performed to date. In order to decipher genetic determinants of HT, we performed the first GWAS followed by replication in a total of 1443 individuals from Croatia.

METHODS

We performed association analysis in a discovery cohort comprising 405 cases and 433 controls. We followed up 13 independent signals (P < 10^-5) in 303 cases and 302 controls from two replication cohorts and then meta-analyzed results across discovery and replication datasets.

RESULTS

We identified three variants suggestively associated with HT: rs12944194 located 206 kb from SDK2 (P = 1.8 × 10^-6), rs75201096 inside GNA14 (P = 2.41 × 10^-5) and rs791903 inside IP6K3 (P = 3.16 × 10^-5). Genetic risk score (GRS), calculated using risk alleles of these loci, accounted for 4.82% of the total HT variance, and individuals from the top GRS quartile had 2.76 times higher odds for HT than individuals from the lowest GRS quartile.

CONCLUSIONS

Although discovered loci are implicated with susceptibility to HT for the first time, genomic regions harboring these loci exhibit good biological candidacy due to involvement in the regulation of the thyroid function and autoimmunity. Additionally, we observe genetic overlap between HT and several related traits, such as hypothyroidism, Graves' disease and TPOAb. Our study adds a new knowledge of underlying HT genetics and sets a firm basis for further research.

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.

Association of Established Thyroid-stimulating Hormone and Free Thyroxine Genetic Variants with Hashimoto's Thyroiditis

Hashimoto's thyroiditis (HT), the most frequent autoimmune thyroid disease (AITD), is characterized by chronic inflammation of the thyroid gland that usually results in hypothyroidism. Thyroid-stimulating hormone (TSH) and free thyroxine (FT4) levels are used as clinical determinants of thyroid function. The main aim of this study was to explore the association of established TSH and FT4 genetic variants with HT. We performed a case-control analysis using 23 genetic markers in 200 HT patients and 304 controls. Additionally, we tested the association of selected variants with several thyroid-related quantitative traits in HT cases only. Two genetic variants showed nominal association with HT: rs11935941 near NR3C2 gene (p = 0.0034, OR = 0.57, 95% CI = 0.39-0.83) and rs1537424 near MBIP gene (p = 0.0169, OR = 0.72, 95% CI = 0.55-0.94). Additionally, three SNPs showed nominal association with thyroglobulin antibody (TgAb) levels: rs4804416 in INSR gene (p = 0.0073, β = -0.51), rs6435953 near IGFBP5 gene (p = 0.0081, β = 0.75), and rs1537424 near MBIP gene (p = 0.0117, β = 0.49). GLIS3 genetic variant rs10974423 showed nominal association with thyroid peroxidase antibody (TPOAb) levels (p = 0.0465, β = -0.56) and NRG1 genetic variant rs7825175 was nominally associated with thyroid gland volume (p = 0.0272, β = -0.18). All detected loci were previously related to thyroid function or pathology. Findings from our study suggest biological relevance of NR3C2 and MBIP with HT, although these loci require additional confirmation in a larger replication study.

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.

Association of established hypothyroidism-associated genetic variants with Hashimoto's thyroiditis

PURPOSE

Hashimoto's thyroiditis (HT) as a chronic autoimmune disease of the thyroid gland is the most common cause of hypothyroidism. Since HT and hypothyroidism are closely related, the main aim of this study was to explore the association of established hypothyroidism single-nucleotide polymorphisms (SNPs) with HT.

METHODS

The case-control dataset included 200 HT cases and 304 controls. Diagnosis of HT cases was based on clinical examination, measurement of thyroid antibodies (TgAb, TPOAb), hormones (TSH and FT4) and ultrasound examination. We genotyped and analysed 11 known hypothyroidism-associated genetic variants. Case-control association analysis was performed in order to test each SNP for the association with HT using logistic regression model. Additionally, each SNP was tested for the association with thyroid-related quantitative traits (TPOAb levels, TgAb levels and thyroid volume) in HT cases only using linear regression.

RESULTS

We identified two genetic variants nominally associated with HT rs3184504 in SH2B3 gene (P = 0.0135, OR = 0.74, 95% CI = 0.57-0.95) and rs4704397 in PDE8B gene (P = 0.0383, OR = 1.32, 95% CI = 1.01-1.74). The SH2B3 genetic variant also showed nominal association with TPOAb levels (P = 0.0163, β = -0.46) and rs4979402 inside DFNB31 gene was nominally associated with TgAb levels (P = 0.0443, β = 0.41).

CONCLUSIONS

SH2B3 gene has previously been associated with susceptibility to several autoimmune diseases, whereas PDE8B has been associated with TSH levels and suggested to modulate thyroid physiology that may influence the manifestation of thyroid disease. Identified loci are novel and biologically plausible candidates for HT development and represent good basis for further exploration of HT susceptibility.

Environmental Risk Factors for Type 1 Diabetes Mellitus Development

Background Although environmental factors induce development of type 1 diabetes mellitus (T1DM) in genetically susceptible individuals, many of those factors have been uncovered. Therefore, the aim of the present study was to analyze associations of T1DM with a wide range of environmental factors. Material and Methods A case-control study was conducted on 249 diabetic and 255 healthy individuals from the Dalmatian region of South Croatia. Data regarding risk factors during pregnancy and early life period of the child were evaluated. Results History of antihypertensive intake (p=0.04) and frequency of stressful life events during pregnancy (p=0.01) were associated with higher risk of T1DM, while hypertension was associated with lower risk of T1DM (p=0.01). Maternal age<25 years at delivery was associated with a higher risk of T1DM (p=0.01).Diabetic patients had a positive family history of T1DM or T2DM (p=0.002) more frequently than controls, while history of infectious diseases was inversely associated with the risk of T1DM (p=0.03). A higher risk of T1DM was significantly associated with earlier introduction of cow's milk (p=0.001), higher number of meals consumed per day (p=0.02), higher frequency of carbohydrate (p=0.001) and meat (p=0.01) consumption and stressful life events during childhood (p=0.02) while earlier introduction of fruit was associated with a lower risk of T1DM (p=0.03) Conclusion This case-control study confirmed associations of a large number of environmental factors with development of T1DM with emphasis on the association of mother's antihypertensive intake during pregnancy, which extends our knowledge about environmental factors related with development of T1DM.