A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function

Genetic Susceptibility to Joint Occurrence of Polycystic Ovary Syndrome and Hashimoto's Thyroiditis: How Far Is Our Understanding?

Polycystic ovary syndrome (PCOS) and Hashimoto’s thyroiditis (HT) are endocrine disorders that commonly occur among young women. A higher prevalence of HT in women with PCOS, relative to healthy individuals, is observed consistently. Combined occurrence of both diseases is associated with a higher risk of severe metabolic and reproductive complications. Genetic factors strongly impact the pathogenesis of both PCOS and HT and several susceptibility loci associated with a higher risk of both disorders have been identified. Furthermore, some candidate gene polymorphisms are thought to be functionally relevant; however, few genetic variants are proposed to be causally associated with the incidence of both disorders together.

Predictive Value of Thyroid Hormones for Incident Hyperuricemia in Euthyroid Subjects: The Tianjin Chronic Low-Grade Systemic Inflammation and Health Cohort Study

OBJECTIVE

Thyroid hormones (THs) play an important role in both serum uric acid (SUA) excretion and purine nucleotide metabolism. Past research mainly focused on the relationships between thyroid dysfunction and hyperuricemia. Although most subjects at risk for hyperuricemia are euthyroid, few studies have investigated the predictive values of THs on incident hyperuricemia in euthyroid adults. This study aimed to examine how free triiodothyronine, thyroxine, and thyroid-stimulating hormone are related to incident hyperuricemia in euthyroid subjects.

METHODS

Participants without baseline hyperuricemia were recruited from Tianjin Chronic Low-grade Systemic Inflammation and Health Cohort Study and followed up for ∼6 years. Thyroid function was determined by chemiluminescence immunoassay methods. Hyperuricemia was defined as SUA ≥416.5 μmol/L for males and ≥357.0 μmol/L for females. Thyroid function and SUA were assessed yearly during follow-up. Adjusted Cox proportional hazards regression models were used to assess the relationship between thyroid function and hyperuricemia.

RESULTS

The incidence rates of hyperuricemia were 109 and 50 per 1000 person-years in males and females, respectively. In males, compared with the lowest quartile, the multivariable-adjusted hazards ratios for hyperuricemia in the highest quartiles of triiodothyronine, thyroxine, and thyroid-stimulating hormone were 0.57 (0.50-0.66), 0.63 (0.54-0.73), and 1.03 (0.90-1.19) (P for trend < .0001, < .0001, and .51), respectively. However, no statistically significant correlations between thyroid function and incident hyperuricemia in females were found.

CONCLUSION

This cohort study is the first to demonstrate that higher THs are related to lower risk of incident hyperuricemia in a male population with euthyroid status.

An update of cyclic nucleotide phosphodiesterase as a target for cardiac diseases

INTRODUCTION

Cyclic nucleotides, cAMP, and cGMP, are important second messengers of intracellular signaling and play crucial roles in cardiovascular biology and diseases. Cyclic nucleotide phosphodiesterases (PDEs) control the duration, magnitude, and compartmentalization of cyclic nucleotide signaling by catalyzing the hydrolysis of cyclic nucleotides. Individual PDEs modulate distinct signaling pathways and biological functions in the cell, making it a potential therapeutic target for the treatment of different cardiovascular disorders. The clinical success of several PDE inhibitors has ignited continued interest in PDE inhibitors and in PDE-target therapeutic strategies.

AREAS COVERED

This review concentrates on recent research advances of different PDE isoforms with regard to their expression patterns and biological functions in the heart. The limitations of current research and future directions are then discussed. The current and future development of PDE inhibitors is also covered.

EXPERT OPINION

Despite the therapeutic success of several marketed PDE inhibitors, the use of PDE inhibitors can be limited by their side effects, lack of efficacy, and lack of isoform selectivity. Advances in our understanding of the mechanisms by which cellular functions are changed through PDEs may enable the development of new approaches to achieve effective and specific PDE inhibition for various cardiac therapies.

The association between longitudinal trends of thyroid hormones levels and incident hypertension in a euthyroid population

Thyroid hormones, including free triiodothyronine (FT3), free thyroxine (FT4), have well-recognized effects on the cardiovascular system. However, the evidence is lacking regarding the relationship between repeated FT3, FT4, and thyroid-stimulating hormone (TSH) measurements and incident hypertension. The aim of this cohort study was to examine how longitudinal trends of serum FT3, FT4, and TSH levels are related to the development of hypertension in a euthyroid population. A prospective study (n = 5926) was performed in Tianjin, China. Participants without a history of hypertension were followed up for ~4 years (median: 3 years). Hypertension was defined according to the criteria of JNC7. FT3, FT4, and TSH were determined by chemiluminescence immunoassay methods. FT3, FT4, TSH, and blood pressure were assessed yearly during follow-up. Adjusted Cox proportional hazards regression models were used to assess the relationships between baseline, means, and annual changes in FT3, FT4, TSH, and hypertension. The incidence rate of hypertension per 1000 person-years was 73. Compared with the lowest quartile, the multivariable-adjusted hazards ratios (95% confidence interval) for hypertension in the highest quartiles of changes in FT3, FT4, and TSH were 1.51 (1.23-1.84), 2.04 (1.67-2.48), and 1.20 (0.99-1.45), respectively. Similar relationships were observed between the means of FT3, FT4, TSH, and hypertension. However, we found no correlations between baseline FT3, FT4, TSH, and incident hypertension. The present study is the first to demonstrate that the annual changes and means, but not baseline FT3 and FT4 values are independently related to the risk of incident hypertension in the euthyroid general population.

Epidemiology of thyroid disorders in the Lifelines Cohort Study (the Netherlands)

BACKGROUND

Thyroid hormone plays a pivotal role in human metabolism. In epidemiologic studies, adequate registration of thyroid disorders is warranted. We examined the prevalence of thyroid disorders, reported thyroid medication use, thyroid hormone levels, and validity of thyroid data obtained from questionnaires in the Lifelines Cohort Study.

METHODS

We evaluated baseline data of all 152180 subjects (aged 18-93 years) of the Lifelines Cohort Study. At baseline, participants were asked about previous thyroid surgery and current and previous thyroid hormone use. At follow-up (n = 136776, after median 43 months), incident thyroid disorders could be reported in an open, non-structured question. Data on baseline thyroid hormone measurements (TSH, FT4 and FT3) were available in a subset of 39935 participants.

RESULTS

Of the 152180 participants, mean (±SD) age was 44.6±13.1 years and 58.5% were female. Thyroid medication was used by 4790 participants (3.1%); the majority (98.2%) used levothyroxine, and 88% were females. 59.3% of levothyroxine users had normal TSH levels. The prevalence of abnormal TSH levels in those not using thyroid medication was 10.8%; 9.4% had a mildly elevated (4.01-10.0 mIU/L), 0.7% had suppressed (<0.40 mIU/L), while 0.7% had elevated (>10.0 mIU/L) TSH levels. Over 98% of subjects with TSH between 4 and 10 mIU/L had normal FT4. Open text questions allowing to report previous thyroid surgery and incident thyroid disorders proved not to be reliable and severely underestimated the true incidence and prevalence of thyroid disorders.

CONCLUSIONS

Undetected thyroid disorders were prevalent in the general population, whereas the prevalence of thyroid medication use was 3.1%. Less than 60% of individuals using levothyroxine had a normal TSH level. The large group of individuals with subclinical hypothyroidism (9.4%) offers an excellent possibility to prospectively follow the natural course of this disorder. Both structured questions as well as linking to G.P.'s and pharmacists' data are necessary to improve the completeness and reliability of Lifelines' data on thyroid disorders.

Transcription factor GLIS3: Critical roles in thyroid hormone biosynthesis, hypothyroidism, pancreatic beta cells and diabetes

GLI-Similar 3 (GLIS3) is a member of the GLIS subfamily of Krüppel-like zinc finger transcription factors that functions as an activator or repressor of gene expression. Study of GLIS3-deficiency in mice and humans revealed that GLIS3 plays a critical role in the regulation of several biological processes and is implicated in the development of various diseases, including hypothyroidism and diabetes. This was supported by genome-wide association studies that identified significant associations of common variants in GLIS3 with increased risk of these pathologies. To obtain insights into the causal mechanisms underlying these diseases, it is imperative to understand the mechanisms by which this protein regulates the development of these pathologies. Recent studies of genes regulated by GLIS3 led to the identification of a number of target genes and have provided important molecular insights by which GLIS3 controls cellular processes. These studies revealed that GLIS3 is essential for thyroid hormone biosynthesis and identified a critical function for GLIS3 in the generation of pancreatic β cells and insulin gene transcription. These observations raised the possibility that the GLIS3 signaling pathway might provide a potential therapeutic target in the management of diabetes, hypothyroidism, and other diseases. To develop such strategies, it will be critical to understand the upstream signaling pathways that regulate the activity, expression and function of GLIS3. Here, we review the recent progress on the molecular mechanisms by which GLIS3 controls key functions in thyroid follicular and pancreatic β cells and how this causally relates to the development of hypothyroidism and diabetes.

Genetic variation in the body mass index of adult survivors of childhood acute lymphoblastic leukemia: A report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort

BACKGROUND

Treatment characteristics such as cranial radiation therapy (CRT) do not fully explain adiposity risk in childhood acute lymphoblastic leukemia (ALL) survivors. This study was aimed at characterizing genetic variation related to adult body mass index (BMI) among survivors of childhood ALL.

METHODS

Genetic associations of BMI among 1458 adult survivors of childhood ALL (median time from diagnosis, 20 years) were analyzed by multiple approaches. A 2-stage genome-wide association study in the Childhood Cancer Survivor Study (CCSS) and the St. Jude Lifetime Cohort Study (SJLIFE) was performed. BMI was a highly polygenic trait in the general population. Within the known loci, the BMI percent variance explained was estimated, and additive interactions (chi-square test) with CRT in the CCSS were evaluated. The role of DNA methylation in CRT interaction was further evaluated in a subsample of ALL survivors.

RESULTS

In a meta-analysis of the CCSS and SJLIFE, 2 novel loci associated with adult BMI among survivors of childhood ALL (LINC00856 rs575792008 and EMR1 rs62123082; P_Meta < 5E-8) were identified. It was estimated that the more than 700 known loci explained 6.2% of the variation in adult BMI in childhood ALL survivors. Within the known loci, significant main effects for 23 loci and statistical interactions with CRT at 9 loci (P < 7.0E-5) were further identified. At 2 CRT-interacting loci, DNA methylation patterns may have differed by age.

CONCLUSIONS

Adult survivors of childhood ALL have genetic heritability for BMI similar to that observed in the general population. This study provides evidence that treatment with CRT can modify the effect of genetic variants on adult BMI in childhood ALL survivors.

Mendelian randomization supports a causative effect of TSH on thyroid carcinoma

Evidence from observational studies suggest a positive association between serum thyroid-stimulating hormone (TSH) levels and differentiated thyroid carcinoma. However, the cause-effect relationship is poorly understood and these studies are susceptible to bias and confounding. This study aimed to investigate the causal role of TSH in both benign thyroid nodules and thyroid cancer in up to 451,025 UK Biobank participants, using a genetic technique, known as Mendelian randomization (MR). Hospital Episode Statistics and Cancer Registry databases were used to identify 462 patients with differentiated thyroid carcinoma and 2031 patients with benign nodular thyroid disease. MR methods using genetic variants associated with serum TSH were used to test causal relationships between TSH and the two disease outcomes. Mendelian randomization provided evidence of a causal link between TSH and both thyroid cancer and benign nodular thyroid disease. Two-sample MR suggested that a 1 s.d. higher genetically instrumented TSH (approximately 0.8 mIU/L) resulted in 4.96-fold higher odds of benign nodular disease (95% CI 2.46-9.99) and 2.00-fold higher odds of thyroid cancer (95% CI 1.09-3.70). Our results thus support a causal role for TSH in both benign nodular thyroid disease and thyroid cancer.

VEGFA and NFE2L2 Gene Expression and Regulation by MicroRNAs in Thyroid Papillary Cancer and Colloid Goiter

Deregulation of VEGFA (Vascular Endothelial Growth Factor A) and NFE2L2 (Nuclear Factor (Erythroid-derived 2)-Like 2), involved in angiogenesis and oxidative stress, can lead to thyroid cancer progression. MiR-17-5p and miR-612 are possible regulators of these genes and may promote thyroid disorders. In order to evaluate the involvement of VEGFA, NFE2L2, hsa-miR-17-5p, and hsa-miR-612 in thyroid pathology, we examined tissue samples from colloid goiter, papillary thyroid cancer (PTC), and a normal thyroid. We found higher levels of VEGFA and NFE2L2 transcripts and the VEGFA protein in goiter and PTC samples than in normal tissue. In the goiter, miR-612 and miR-17-5p levels were lower than those in PTC. Tumors, despite showing lower VEGFA mRNA expression, presented higher VEGFA protein levels compared to goiter tissue. In addition, NRF2 (Nuclear Related Transcription Factor 2) protein levels in tumors were higher than those in goiter and normal tissues. Inhibition of miR-17-5p resulted in reduced NFE2L2 expression. Overall, both transcript and protein levels of NFE2L2 and VEGFA were elevated in PTC and colloid goiter. Hsa-miR-612 showed differential expression in PTC and colloid goiter, while hsa-miR-17-5p showed differential expression only in colloid goiter, suggesting that hsa-miR-17-5p may be a positive regulator of NFE2L2 expression in PTC.

Targeting the right population for T3 + T4 combined therapy: where are we now and where to next?

The universal applicability of levothyroxine (LT4) monotherapy for the treatment of hypothyroidism has been questioned in recent years. Indeed, it is now clear that about 10-15% of LT4-treated hypothyroid patients are dissatisfied with their treatment. It is plausible that this subset of hypothyroid patients may need T3 + T4 combined therapy to restore peripheral euthyroidism. To address this issue, many clinical trials have investigated the effect of T3 + T4 combinations versus standard LT4-based therapy. However, to date, results have been inconclusive, mainly due to the lack of markers that identify candidates for combination therapy. A breakthrough in this field came with the recent finding that several single-nucleotide polymorphisms in the deiodinase genes are associated with the persistence of hypothyroid symptoms in biochemically euthyroid LT4-treated patients, and are thus markers of candidates for combination therapy. In addition, whole-genome association studies are expanding our knowledge of other genes of the thyroid hormone (TH) pathway that affect serum TH levels. To target the right population for the T3 + T4 combined therapy, the next step is to translate these new findings into prospective trials. Hopefully, this will pave the way to personalized therapy for each hypothyroid patient.

Establishment of reference interval for thyroid-stimulating hormone using electrochemiluminescence assay in a healthy adult population from Fortaleza, Brazil

Objective This study aimed to determine the thyroid-stimulating hormone (TSH) reference interval (RI) and to assess the influence of the use of thyroid ultrasonography (TUS) on reference individual selection from a healthy adult population in Fortaleza, Brazil. Subjects and methods This cross-sectional study recruited patients (N = 272; age = 18-50 years) with normal thyroid function (NTF) and placed them in three groups according to their test results: NTF (n = 272; all participants), TUS (n = 170; participants who underwent thyroid US), RI (n = 124; reference individuals with normal TSH levels). TSH, FT4, TT3, TgAb, and TPOAb concentrations were determined by electrochemiluminescence assay. TUS was performed using a 7-12 MHz multifrequency linear transducer by two radiologists. The 2.5th and 97.5th percentiles of the distribution curve corresponded to lower and upper TSH RI levels, respectively. Results The mean TSH level was 1.74 ± 0.96 mIU/L, and TSH range was 0.56-4.45 mIU/L. There was no difference in the TSH concentrations between men and women nor between the groups. TUS did not appear to be an essential tool for the reference group selection. Conclusion The upper limit of TSH was comparable to the reference interval provided by the assay manufacturer (4.45 vs. 4.20 mIU/L) but the lower limit was not (0.56 vs. 0.27 mIU/L). This finding may have a clinical impact since these values may lead to the misdiagnosis of euthyroid patients with subclinical hyperthyroidism.

Thyroid Function Affects the Risk of Stroke via Atrial Fibrillation: A Mendelian Randomization Study

CONTEXT

Observational studies suggest that variations in normal range thyroid function are associated with cardiovascular diseases. However, it remains to be determined whether these associations are causal or not.

OBJECTIVE

To test whether genetically determined variation in normal range thyroid function is causally associated with the risk of stroke and coronary artery disease (CAD) and investigate via which pathways these relations may be mediated.

DESIGN, SETTING, AND PARTICIPANTS

Mendelian randomization analyses for stroke and CAD using genetic instruments associated with normal range thyrotropin (TSH) and free thyroxine levels or Hashimoto's thyroiditis and Graves' disease. The potential mediating role of known stroke and CAD risk factors was examined. Publicly available summary statistics data were used.

MAIN OUTCOME MEASURES

Stroke or CAD risk per genetically predicted increase in TSH or FT4 levels.

RESULTS

A 1 standard deviation increase in TSH was associated with a 5% decrease in the risk of stroke (odds ratio [OR], 0.95; 95% confidence interval [CI], 0.91-0.99; P = 0.008). Multivariable MR analyses indicated that this effect is mainly mediated via atrial fibrillation. MR analyses did not show a causal association between normal range thyroid function and CAD. Secondary analyses showed a causal relationship between Hashimoto's thyroiditis and a 7% increased risk of CAD (OR, 1.07; 95% CI, 1.01-1.13; P = 0.026), which was mainly mediated via body mass index.

CONCLUSION

These results provide important new insights into the causal relationships and mediating pathways between thyroid function, stroke, and CAD. We identify variation in normal range thyroid function and Hashimoto's thyroiditis as risk factors for stroke and CAD, respectively.

What do the structures of GCN5-containing complexes teach us about their function?

Transcription initiation is a major regulatory step in eukaryotic gene expression. It involves the assembly of general transcription factors and RNA polymerase II into a functional pre-initiation complex at core promoters. The degree of chromatin compaction controls the accessibility of the transcription machinery to template DNA. Co-activators have critical roles in this process by actively regulating chromatin accessibility. Many transcriptional coactivators are multisubunit complexes, organized into distinct structural and functional modules and carrying multiple regulatory activities. The first nuclear histone acetyltransferase (HAT) characterized was General Control Non-derepressible 5 (Gcn5). Gcn5 was subsequently identified as a subunit of the HAT module of the Spt-Ada-Gcn5-acetyltransferase (SAGA) complex, which is an experimental paradigm for multifunctional co-activators. We know today that Gcn5 is the catalytic subunit of multiple distinct co-activator complexes with specific functions. In this review, we summarize recent advances in the structure of Gcn5-containing co-activator complexes, most notably SAGA, and discuss how these new structural insights contribute to better understand their functions.

Evaluation of genetic susceptibility of common variants in SOX9 in patients with congenital talipes equinovarus in the Han Chinese population

Background

Congenital talipes equinovarus (CTEV) is a common birth defect that causes severe deformities of one or both feet. Genetics have been proven to play a key role in the risk of CTEV. Our study aimed to evaluate the genetic susceptibility of common variants in the SOX9 gene to CTEV in a Han Chinese population.

Methods

In this study, we recruited 2,205 study participants, including 692 CTEV patients and 1513 healthy controls. A total of seven selected single-nucleotide polymorphisms (SNPs) within the SOX9 gene were genotyped, and environmental variables, including maternal smoking and alcoholic drinking habits, were assessed. In addition, bioinformatics analyses were performed to explore the potential biological functions of the associated SNPs.

Results

The SNP rs73354570 was identified to be significantly associated with the risk of CTEV (OR = 1.53, P = 2.11 × 10⁻⁵), and the C allele was associated with an increased risk of CTEV. A dose-dependent pattern could be observed in genotypic analyses. The OR for individuals with AC genotypes was 1.37 (95% CI 1.09–1.71), and the OR for individuals with CC homozygotes was 1.47 (95% CI 1.18–1.82). Further analyses identified that rs73354570 is located within a region of multiple binding proteins, including CEBPB and POLR2A, which suggested that this SNP was also part of genetic motifs that are found within several cell types.

Conclusion

Our results provide evidence supporting the important role of the SOX9 gene in the contribution to the risk of CTEV.

Sex-Specific Role for the Long Non-coding RNA LINC00473 in Depression

Depression is a common disorder that affects women at twice the rate of men. Here, we report that long non-coding RNAs (lncRNAs), a recently discovered class of regulatory transcripts, represent about one-third of the differentially expressed genes in the brains of depressed humans and display complex region- and sex-specific patterns of regulation. We identified the primate-specific, neuronal-enriched gene LINC00473 as downregulated in prefrontal cortex (PFC) of depressed females but not males. Using viral-mediated gene transfer to express LINC00473 in adult mouse PFC neurons, we mirrored the human sex-specific phenotype by inducing stress resilience solely in female mice. This sex-specific phenotype was accompanied by changes in synaptic function and gene expression selectively in female mice and, along with studies of human neuron-like cells in culture, implicates LINC00473 as a CREB effector. Together, our studies identify LINC00473 as a female-specific driver of stress resilience that is aberrant in female depression.

Importance of Thyroid Hormone level and Genetic Variations in Deiodinases for Patients after Acute Myocardial Infarction: A Longitudinal Observational Study

This study aimed to examine the influence of thyroid hormone (TH) levels and genetic polymorphisms of deiodinases on long-term outcomes after acute myocardial infarction (AMI). In total, 290 patients who have experienced AMI were evaluated for demographic, clinical characteristics, risk factors, TH and NT-pro-BNP. Polymorphisms of TH related genes were included deiodinase 1 (DIO1) (rs11206244-C/T, rs12095080-A/G, rs2235544-A/C), deiodinase 2 (DIO2) (rs225015-G/A, rs225014-T/C) and deiodinase 3 (DIO3) (rs945006-T/G). Both all-cause and cardiac mortality was considered key outcomes. Cox regression model showed that NT-pro-BNP (HR = 2.11; 95% CI = 1.18– 3.78; p = 0.012), the first quartile of fT3, and DIO1 gene rs12095080 were independent predictors of cardiac-related mortality (HR = 1.74; 95% CI = 1.04–2.91; p = 0.034). The DIO1 gene rs12095080 AG genotype (OR = 3.97; 95% CI = 1.45–10.89; p = 0.005) increased the risk for cardiac mortality. Lower fT3 levels and the DIO1 gene rs12095080 are both associated with cardiac-related mortality after AMI.

The Genetic Basis of Thyroid Function: Novel Findings and New Approaches

CONTEXT

Genetic factors are major determinants of thyroid function. Over the last two decades, multiple genetic variants have been associated with variations in normal range thyroid function tests. Most recently, a large-scale genome-wide association study (GWAS) doubled the number of known variants associated with normal range thyrotropin (TSH) and free thyroxine (FT4) levels.

EVIDENCE ACQUISITION

This review summarizes the results of genetic association studies on normal range thyroid function and explores how these genetic variants can be used in future studies to improve our understanding of thyroid hormone regulation and disease.

EVIDENCE SYNTHESIS

Serum TSH and FT4 levels are determined by multiple genetic variants on virtually all levels of the hypothalamus-pituitary-thyroid (HPT) axis. Functional follow-up studies on top of GWAS hits has the potential to discover new key players in thyroid hormone regulation, as exemplified by the identification of the thyroid hormone transporter SLC17A4 and the metabolizing enzyme AADAT. Translational studies may use these genetic variants to investigate causal associations between thyroid function and various outcomes in Mendelian Randomization (MR) studies, to identify individuals with an increased risk of thyroid dysfunction, and to predict the individual HPT axis setpoint.

CONCLUSIONS

Recent genetic studies have greatly improved our understanding of the genetic basis of thyroid function, and have revealed novel pathways involved in its regulation. In addition, these findings have paved the way for various lines of research that can improve our understanding of thyroid hormone regulation and thyroid diseases, as well as the potential use of these markers in future clinical practice.

Thyroid Hormone Transporters

Thyroid hormone transporters at the plasma membrane govern intracellular bioavailability of thyroid hormone. Monocarboxylate transporter (MCT) 8 and MCT10, organic anion transporting polypeptide (OATP) 1C1, and SLC17A4 are currently known as transporters displaying the highest specificity toward thyroid hormones. Structure-function studies using homology modeling and mutational screens have led to better understanding of the molecular basis of thyroid hormone transport. Mutations in MCT8 and in OATP1C1 have been associated with clinical disorders. Different animal models have provided insight into the functional role of thyroid hormone transporters, in particular MCT8. Different treatment strategies for MCT8 deficiency have been explored, of which thyroid hormone analogue therapy is currently applied in patients. Future studies may reveal the identity of as-yet-undiscovered thyroid hormone transporters. Complementary studies employing animal and human models will provide further insight into the role of transporters in health and disease. (Endocrine Reviews 41: 1 - 55, 2020).

Hypothyroidism and Kidney Function: A Mendelian Randomization Study

Background: Uncertainty in the mechanism and directionality of observational associations between thyroid function and kidney function may be addressed by genetic analysis with an instrumental variable method termed bidirectional Mendelian randomization (MR). Methods: In the Women's Genome Health Study (WGHS), observational associations between thyroid measures and kidney function were evaluated. Genetic instruments for MR were from recent genome-wide association studies (GWAS) of hypothyroidism, thyrotropin (TSH), and free thyroxine (fT4) concentrations within the reference range, thyroid peroxidase antibodies (TPOAb), estimated glomerular filtration rate from creatinine (eGFR_crea), eGFR from cystatin C (eGFR_cys), and chronic kidney disease (CKD). In WGHS individual-level data, these instruments were used for bidirectional MR between thyroid (N = 3336) and kidney (N = 23,186) functions. To increase power, MR was also performed using GWAS summary statistics from the Chronic Kidney Disease Genetics Consortium (CKDGen) for eGFR_crea (N = 567,460), eGFR_cys (N = 24,063), CKD [N(total) = 480,698, N(cases) = 41,395], and urinary albumin/creatinine ratio (UACR/N = 54,450). Results: In the WGHS, hypothyroidism was observationally associated with decreased eGFR_crea [beta (standard error, SE): -0.024 (0.009) ln(mL/min/1.73 m²), p = 0.01]. By MR, hypothyroidism was associated with decreased eGFR_crea in the WGHS [beta (SE): -0.007 (0.002) per doubled odds hypothyroidism, p = 1.7 × 10^-3] and in CKDGen [beta (SE): -0.004 (0.0005), p = 2.0 × 10^-22], and robust to sensitivity analysis. Hypothyroidism was also associated by MR with increased CKD in CKDGen (odds ratio, OR [confidence interval, CI]: 1.05 [1.03-1.08], p = 3.3 × 10^-5), but not in the WGHS (OR [CI]: 1.02 [0.95-1.10], p = 0.57). Increased TSH within the reference range had an MR association with increased eGFR_crea in the WGHS [beta (SE): -0.018 (0.007) ln(mL/min/1.73 m²)/standard deviation, SD, p = 6.5 × 10^-3] and CKDGen [beta (SE): -0.008 (0.001) ln(mL/min/1.73 m²)/SD, p = 6.8 × 10^-17], and with CKD in CKDGen (OR [CI]: 1.10 [1.04-1.15], p = 3.1 × 10^-4). There were no MR associations of hypothyroidism or TSH with eGFR_cys or UACR, and MR associations of fT4 in the reference range with kidney function were inconsistent in both the WGHS and CKDGen. However, by MR in CKDGen, TPOAb were robustly associated with decreased eGFR_crea [beta (SE): -0.041 (0.009), p = 6.2 × 10^-6] and decreased eGFR_cys [beta (SE): -0.294 (0.065), p = 6.2 × 10^-6]. TPOAb were less robustly associated with CKD but not associated with UACR. In reverse MR in the WGHS, kidney function was not consistently associated with thyroid function. Conclusions: Bidirectional MR supports a directional association from hypothyroidism, increased TSH, and TPOAb, but not fT4, to decreased eGFR_crea and increased CKD.

Association of Thyroid Function Genetic Predictors With Atrial Fibrillation: A Phenome-Wide Association Study and Inverse-Variance Weighted Average Meta-analysis

Importance

Thyroid hormone levels are tightly regulated through feedback inhibition by thyrotropin, produced by the pituitary gland. Hyperthyroidism is overwhelmingly due to thyroid disorders and is well recognized to contribute to a wide spectrum of cardiovascular morbidity, particularly the increasingly common arrhythmia atrial fibrillation (AF).

Objective

To determine the association between genetically determined thyrotropin levels and AF.

Design, Setting, and Participants

This phenome-wide association study scanned 1318 phenotypes associated with a polygenic predictor of thyrotropin levels identified by a previously published genome-wide association study that included participants of European ancestry. North American individuals of European ancestry with longitudinal electronic health records were analyzed from May 2008 to November 2016. Analysis began March 2018.

Main Outcomes and Measures

Clinical diagnoses associated with a polygenic predictor of thyrotropin levels.

Exposures

Genetically determined thyrotropin levels.

Results

Of 37 154 individuals, 19 330 (52%) were men. The thyrotropin polygenic predictor was positively associated with hypothyroidism (odds ratio [OR], 1.10; 95% CI, 1.07-1.14; P = 5 × 10-11) and inversely associated with diagnoses related to hyperthyroidism (OR, 0.64; 95% CI, 0.54-0.74; P = 2 × 10-8 for toxic multinodular goiter). Among nonthyroid associations, the top association was AF/flutter (OR, 0.93; 95% CI, 0.9-0.95; P = 9 × 10-7). When the analyses were repeated excluding 9801 individuals with any diagnoses of a thyroid-related disease, the AF association persisted (OR, 0.91; 95% CI, 0.88-0.95; P = 2.9 × 10-6). To replicate this association, we conducted an inverse-variance weighted average meta-analysis using AF single-nucleotide variant weights from a genome-wide association study of 17 931 AF cases and 115 142 controls. As in the discovery analyses, each SD increase in predicted thyrotropin was associated with a decreased risk of AF (OR, 0.86; 95% CI, 0.79-0.93; P = 4.7 × 10-4). In a set of AF cases (n = 745) and controls (n = 1680) older than 55 years, directly measured thyrotropin levels that fell within the normal range were inversely associated with AF risk (OR, 0.91; 95% CI, 0.83-0.99; P = .04).

Conclusions and Relevance

This study suggests a role for genetically determined variation in thyroid function within a physiologically accepted normal range as a risk factor for AF. The clinical decision to treat subclinical thyroid disease should incorporate the risk for AF as antithyroid medications to treat hyperthyroidism may reduce AF risk and thyroid hormone replacement for hypothyroidism may increase AF risk.

Assessment of the Relationship Between Genetic Determinants of Thyroid Function and Atrial Fibrillation: A Mendelian Randomization Study

Importance

Increased free thyroxine (FT4) and decreased thyrotropin are associated with increased risk of atrial fibrillation (AF) in observational studies, but direct involvement is unclear.

Objective

To evaluate the potential direct involvement of thyroid traits on AF.

Design, Setting, and Participants

Study-level mendelian randomization (MR) included 11 studies, and summary-level MR included 55 114 AF cases and 482 295 referents, all of European ancestry.

Exposures

Genomewide significant variants were used as instruments for standardized FT4 and thyrotropin levels within the reference range, standardized triiodothyronine (FT3):FT4 ratio, hypothyroidism, standardized thyroid peroxidase antibody levels, and hyperthyroidism. Mendelian randomization used genetic risk scores in study-level analysis or individual single-nucleotide polymorphisms in 2-sample MR for the summary-level data.

Main Outcomes and Measures

Prevalent and incident AF.

Results

The study-level analysis included 7679 individuals with AF and 49 233 referents (mean age [standard error], 62 [3] years; 15 859 men [29.7%]). In study-level random-effects meta-analysis, the pooled hazard ratio of FT4 levels (nanograms per deciliter) for incident AF was 1.55 (95% CI, 1.09-2.20; P = .02; I2 = 76%) and the pooled odds ratio (OR) for prevalent AF was 2.80 (95% CI, 1.41-5.54; P = .003; I2 = 64%) in multivariable-adjusted analyses. The FT4 genetic risk score was associated with an increase in FT4 by 0.082 SD (standard error, 0.007; P < .001) but not with incident AF (risk ratio, 0.84; 95% CI, 0.62-1.14; P = .27) or prevalent AF (OR, 1.32; 95% CI, 0.64-2.73; P = .46). Similarly, in summary-level inverse-variance weighted random-effects MR, gene-based FT4 within the reference range was not associated with AF (OR, 1.01; 95% CI, 0.89-1.14; P = .88). However, gene-based increased FT3:FT4 ratio, increased thyrotropin within the reference range, and hypothyroidism were associated with AF with inverse-variance weighted random-effects OR of 1.33 (95% CI, 1.08-1.63; P = .006), 0.88 (95% CI, 0.84-0.92; P < .001), and 0.94 (95% CI, 0.90-0.99; P = .009), respectively, and robust to tests of horizontal pleiotropy. However, the subset of hypothyroidism single-nucleotide polymorphisms involved in autoimmunity and thyroid peroxidase antibodies levels were not associated with AF. Gene-based hyperthyroidism was associated with AF with MR-Egger OR of 1.31 (95% CI, 1.05-1.63; P = .02) with evidence of horizontal pleiotropy (P = .045).

Conclusions and Relevance

Genetically increased FT3:FT4 ratio and hyperthyroidism, but not FT4 within the reference range, were associated with increased AF, and increased thyrotropin within the reference range and hypothyroidism were associated with decreased AF, supporting a pathway involving the pituitary-thyroid-cardiac axis.

Common genetic variations of deiodinase genes and prognosis of brain tumor patients

BACKGROUND

Thyroid hormone (TH) metabolism can have prognostic significance in brain tumors. We studied the association of common variations in three deiodinase gene single-nucleotide polymorphisms (SNPs) with circulating TH concentrations and prognosis of brain tumor patients.

METHODS

Patients admitted for glioma and meningioma surgery between January, 2010 and September, 2011 were evaluated for functional status (Barthel Index or BI) and circulating free tri-iodothyronine (FT3), free thyroxine (FT4), and thyroid-stimulating hormone (TSH) concentrations. Ten common SNPs in the DIO1 gene; five SNPs in the DIO2 gene; and one SNP in the DIO3 gene were genotyped. Follow-up continued until November, 2017.

RESULTS

In glioblastoma patients, the DIO1 SNP rs2235544 CC genotype was associated with significantly lower risk of death at 2 years when compared to AA + CA genotypes after adjusting for patient gender, age, pre-operative functional status, adjuvant therapy, and extent of resection (HR = 0.34, 95% CI: 0.13-0.84, p = 0.019). The TT genotype vs. CC + TC genotypes of the DI02 SNP rs12885300 was associated with increased mortality risk after adjusting for patient gender, age, pre-operative functional status, adjuvant therapy, extent of resection, and FT3/FT4 (HR = 3.13, 95% CI: 1.20-8.16, p < 0.019). The C-allele of the DI01 SNP rs2235544 was related to increased circulating free T3/ free T4 ratio in glioma and meningioma patients, indicating greater T4 to T3 conversion.

CONCLUSIONS

SNPs of DIO1 gene (rs2235544) and DIO2 gene (rs12885300) have independent prognostic significance in glioblastoma patients. The C-allele of the DIO1 (rs2235544) is associated with greater T4 to T3 conversion.

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.

Conservation, acquisition, and functional impact of sex-biased gene expression in mammals

Sex differences abound in human health and disease, as they do in other mammals used as models. The extent to which sex differences are conserved at the molecular level across species and tissues is unknown. We surveyed sex differences in gene expression in human, macaque, mouse, rat, and dog, across 12 tissues. In each tissue, we identified hundreds of genes with conserved sex-biased expression-findings that, combined with genomic analyses of human height, explain ~12% of the difference in height between females and males. We surmise that conserved sex biases in expression of genes otherwise operating equivalently in females and males contribute to sex differences in traits. However, most sex-biased expression arose during the mammalian radiation, which suggests that careful attention to interspecies divergence is needed when modeling human sex differences.

Thyroid and Cardiovascular Disease: Research Agenda for Enhancing Knowledge, Prevention, and Treatment

Thyroid hormones have long been known to have a range of effects on the cardiovascular system. However, significant knowledge gaps exist concerning the precise molecular and biochemical mechanisms governing these effects and the optimal strategies for management of abnormalities in thyroid function in patients with and without preexisting cardiovascular disease. In September 2017, the National Heart, Lung, and Blood Institute convened a Working Group with the goal of developing priorities for future scientific research relating thyroid dysfunction to the progression of cardiovascular disease. The Working Group reviewed and discussed the roles of normal thyroid physiology, the consequences of thyroid dysfunction, and the effects of therapy in 3 cardiovascular areas: cardiac electrophysiology and arrhythmias, the vasculature and atherosclerosis, and the myocardium and heart failure. This report describes the current state of the field, outlines barriers and challenges to progress, and proposes research opportunities to advance the field, including strategies for leveraging novel approaches using omics and big data. The Working Group recommended research in 3 broad areas: (1) investigation into the fundamental biology relating thyroid dysfunction to the development of cardiovascular disease and into the identification of novel biomarkers of thyroid hormone action in cardiovascular tissues; (2) studies that define subgroups of patients with thyroid dysfunction amenable to specific preventive strategies and interventional therapies related to cardiovascular disease; and (3) clinical trials focused on improvement in cardiovascular performance and cardiovascular outcomes through treatment with thyroid hormone or thyromimetic drugs.

Thyroid and Cardiovascular Disease: Research Agenda for Enhancing Knowledge, Prevention, and Treatment

Thyroid hormones have long been known to have a range of effects on the cardiovascular system. However, significant knowledge gaps exist concerning the precise molecular and biochemical mechanisms governing these effects and the optimal strategies for management of abnormalities in thyroid function in patients with and without preexisting cardiovascular disease. In September 2017, The National Heart, Lung, and Blood Institute convened a Working Group with the goal of developing priorities for future scientific research relating thyroid dysfunction to the progression of cardiovascular disease. The Working Group reviewed and discussed the roles of normal thyroid physiology, the consequences of thyroid dysfunction, and the effects of therapy in three cardiovascular areas: cardiac electrophysiology and arrhythmias, the vasculature and atherosclerosis, and the myocardium and heart failure. This report describes the current state of the field, outlines barriers and challenges to progress, and proposes research opportunities to advance the field, including strategies for leveraging novel approaches using omics and big data. The Working Group recommended research in three broad areas: 1) investigation into the fundamental biology relating thyroid dysfunction to the development of cardiovascular disease and into the identification of novel biomarkers of thyroid hormone action in cardiovascular tissues; 2) studies that define subgroups of patients with thyroid dysfunction amenable to specific preventive strategies and interventional therapies related to cardiovascular disease; and 3) clinical trials focused on improvement in cardiovascular performance and cardiovascular outcomes through treatment with thyroid hormone or thyromimetic drugs.

Reference intervals in the diagnosis of thyroid dysfunction: treating patients not numbers

Although assigning a diagnosis of thyroid dysfunction appears quite simple, this is often not the case. Issues that make it unclear whether thyroid function is normal include transient changes in thyroid parameters, inter-individual and intra-individual differences in thyroid parameters, age-related differences, and ethnic variations. In addition, a statistically calculated distribution of thyroid analytes does not necessarily coincide with intervals or cutoffs that have predictive value for beneficial or adverse health outcomes. Based on current clincial trial data, it is unclear which individuals with mild thyroid-stimulating hormone elevations will benefit from levothyroxine treatment. For example, only a small number of patients with thyroid-stimulating hormone values of more than 10 mIU/L have been studied in a randomised manner. Even if therapy is initiated for abnormal thyroid function, not all treated individuals are maintained at the desired treatment target, and therefore might still be at risk. The consequence of this is that each patient's thyroid function needs to be assessed on an individual basis with the entire clinical picture in mind. Monitoring also needs to be vigilant, and the targets for treatment reassessed continually.

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.

Predicting adult height from DNA variants in a European-Asian admixed population

Accurate genomic profiling for adult height is of high practical relevance in forensics genetics. Adult height is a classical reference trait in the field of human complex trait genetics characterized by highly polygenic nature and relatively high heritability. A meta-analysis of genome-wide association studies by the Genetic Investigation of Anthropocentric Traits (GIANT) consortium has identified 697 DNA variants associated with adult height in Europeans; however, whether these variants will still be informative in non-Europeans is still in question. The present study investigated the predictive power of these 697 height-associated SNPs in 687 Uyghurs of European-Asian admixed origin. Among all GIANT SNPs, 11% showed nominally significant association (6.78 × 10^-4 < p < 0.05) with adult height in the Uyghur population and among the significant SNPs 77% of allele effects were in the same direction as those in Europeans reported in the GIANT study. Fitting linear and logistic models using a polygenic score consisting of all GIANT SNPs resulted in an 80-20 cross-validated mean R² of 10.08% (95% CI 3.16-18.40%) for quantitative height prediction and a mean AUC value of 0.65 (95% CI 0.57-0.72%) for qualitative "above average" prediction. Fine-tuning the SNP set using their association p values considerably improved the prediction results (number of SNPs = 62, R² = 15.59%, 95% CI 6.80-25.71%; AUC = 0.70, 95% CI 62-0.77) in the Uyghurs. Overall, our findings demonstrate substantial differences between the European and Asian populations in the genetics of adult height, emphasizing the importance of population heterogeneity underlying the genetic architecture of adult height.

Genome-wide association meta-analysis for total thyroid hormone levels in Croatian population

Thyroid hormones (THs) are key regulators of cellular growth, development, and metabolism. The thyroid gland secretes two THs, thyroxine (T4) and triiodothyronine (T3), into the plasma where they are almost all bound reversibly to plasma proteins. Free forms of THs are metabolically active, however, they represent a very small fraction of total TH levels. No genome-wide studies have been performed to date on total TH levels, comprising of protein-bound and free forms of THs. To detect genetic variants associated with total TH levels, we carried out the first GWAS meta-analysis of total T4 levels in 1121 individuals from two Croatian cohorts (Split and Korcula). We also performed GWAS analyses of total T3 levels in 577 individuals and T3/T4 ratio in 571 individuals from the Split cohort. The top association in GWAS meta-analysis of total T4 was detected for an intronic variant within SLC22A9 gene (rs12282281, P = 4.00 × 10^-7). Within the same region, a genome-wide significant variant (rs11822642, P = 2.50 × 10^-8) for the T3/T4 ratio was identified. SLC22A9 encodes for an organic anion transporter protein expressed predominantly in the liver and belongs to the superfamily of solute carriers (SLC), a large group of transport membrane proteins. The transport of THs across the plasma membrane in peripheral tissues is facilitated by the membrane proteins, and all TH transport proteins known to date belong to the same SLC superfamily as SLC22A9. These results suggest a potential role for SLC22A9 as a novel transporter protein of THs.

Interpreting Elevated TSH in Older Adults

Thyroid function is most often monitored clinically through measurement of the regulatory hormone, thyrotropin (TSH). Subclinical hypothyroidism is generally defined as a TSH level above the reference range while thyroid hormone levels remain within the reference range. Elevated TSH is more common among older adults, leading to high rates of treatment, and over-treatment, in this population. However, the use of levothyroxine in older adults with mild TSH elevations has begun to be called into question by observations that demonstrate a lack of harm from not treating and a lack of benefit from treating. Importantly, these findings suggest that the existing diagnostic algorithm for subclinical hypothyroidism, based on isolated TSH elevation, may be inappropriate for older adults. Age-specific reference ranges have been suggested as a way to avoid inappropriate treatment, but that strategy continues to rely on population norms rather than disease definitions to drive clinical decisions. Recent insight into age-related variability in the underlying pathophysiology that impacts on thyroid function tests demonstrates the need for new clinical tools to allow the targeted use of therapy where it will have benefit.

Functional Urate-Associated Genetic Variants Influence Expression of lincRNAs LINC01229 and MAFTRR

Genetic variation in the genomic regulatory landscape likely plays a crucial role in the pathology of disease. Non-coding variants associated with disease can influence the expression of long intergenic non-coding RNAs (lincRNAs), which in turn function in the control of protein-coding gene expression. Here, we investigate the function of two independent serum urate-associated signals (SUA1 and SUA2) in close proximity to lincRNAs and an enhancer that reside ∼60 kb and ∼300 kb upstream of MAF, respectively. Variants within SUA1 are expression quantitative trait loci (eQTL) for LINC01229 and MAFTRR, both co-expressed with MAF. We have also identified that variants within SUA1 are trans-eQTL for genes that are active in kidney- and serum urate-relevant pathways. Serum urate-associated variants rs4077450 and rs4077451 within SUA2 lie within an enhancer that recruits the transcription factor HNF4α and forms long range interactions with LINC01229 and MAFTRR. The urate-raising alleles of rs4077450 and rs4077451 increase enhancer activity and associate with increased expression of LINC01229. We show that the SUA2 enhancer region drives expression in the zebrafish pronephros, recapitulating endogenous MAF expression. Depletion of MAFTRR and LINC01229 in HEK293 cells in turn lead to increased MAF expression. Collectively, our results are consistent with serum urate variants mediating long-range transcriptional regulation of the lincRNAs LINC01229 and MAFTRR and urate relevant genes (e.g., SLC5A8 and EHHADH) in trans.

Functional and Symptomatic Individuality in the Response to Levothyroxine Treatment

Background: For significant numbers of patients dissatisfied on standard levothyroxine (LT4) treatment for hypothyroidism, patient-specific responses to T4 could play a significant role. Aim: To assess response heterogeneity to LT4 treatment, identifying confounders and hidden clusters within a patient panel, we performed a secondary analysis using data from a prospective cross-sectional and retrospective longitudinal study. Methods: Multivariate and multivariable linear models adjusted for covariates (gender, age, and BMI) were stratified by disease-specific treatment indication. During follow-up, pooled observations were compared from the same patient presenting either with or without self-reported symptoms. Statistical analysis was extended to multilevel models to derive intra-class correlation coefficients and reliability measures during follow-up. Results: Equilibria between TSH, FT4, and FT3 serum concentrations in 342 patients were examined by treatment indication (benign goiter, autoimmune thyroiditis, thyroid carcinoma), consequently displaying complex interactive response patterns. Seventy-seven patients treated with LT4 and monitored for thyroid carcinoma presented, in association with changes in LT4 dose, either with hypothyroid symptoms or symptom-free. Significant biochemical differences appeared between the different presentations. Leveled trajectories by subject to relief from hypothyroid symptoms differed significantly, indicating distinct responses, and denying a single shared outcome. These were formally defined by a high coefficient of the intraclass correlation (ICC1, exceeding 0.60 in all thyroid parameters) during follow-up on multiple visits at the same LT4 dose, when lacking symptoms. The intra-personal clusters were clearly differentiated from random variability by random group resampling. Symptomatic change in these patients was strongly associated with serum FT3, but not with FT4 or TSH concentrations. In 25 patients transitioning from asymptomatic to symptomatically hyperthyroid, FT3 concentrations remained within the reference limits, whilst at the same time marked biochemical differences were apparent between the presentations. Conclusions: Considerable intra-individual clustering occurred in the biochemical and symptomatic responses to LT4 treatment, implying statistically multileveled response groups. Unmasking individual differences in the averaged treatment response hereby highlights clinically distinguishable subgroups within an indiscriminate patient panel. This, through well-designed larger clinical trials will better target the different therapeutic needs of individual patients.

A Guide to Applying the Sex-Gender Perspective to Nutritional Genomics

Precision nutrition aims to make dietary recommendations of a more personalized nature possible, to optimize the prevention or delay of a disease and to improve health. Therefore, the characteristics (including sex) of an individual have to be taken into account as well as a series of omics markers. The results of nutritional genomics studies are crucial to generate the evidence needed so that precision nutrition can be applied. Although sex is one of the fundamental variables for making recommendations, at present, the nutritional genomics studies undertaken have not analyzed, systematically and with a gender perspective, the heterogeneity/homogeneity in gene-diet interactions on the different phenotypes studied, thus there is little information available on this issue and needs to be improved. Here we argue for the need to incorporate the gender perspective in nutritional genomics studies, present the general context, analyze the differences between sex and gender, as well as the limitations to measuring them and to detecting specific sex-gene or sex-phenotype associations, both at the specific gene level or in genome-wide-association studies. We analyzed the main sex-specific gene-diet interactions published to date and their main limitations and present guidelines with recommendations to be followed when undertaking new nutritional genomics studies incorporating the gender perspective.

Fine mapping of 14q13 reveals novel variants associated with different histological subtypes of papillary thyroid carcinoma

The first two genome wide association studies (GWAS) of papillary thyroid carcinoma (PTC) detected five variants associated with PTC. Two of them (rs944289 and rs116909374) are located at 14q13 making that locus an important target of research into the genetic predisposition to PTC. We aimed at uncovering other variants at 14q13 associated with PTC independently from the GWAS variants. We performed next generation sequencing of the 14q13 region and analyzed the allele frequencies of single nucleotide polymorphisms (SNPs) in n = 90 PTC cases vs. n = 379 EUR controls from the 1,000 Genome Project. The variants associated with PTC were validated in an Ohio cohort of n = 1,216 PTC cases and n = 1,416 controls. Next, we analyzed the association between SNPs and expression of nearby genes and clinical parameters. We showed that rs368187 was associated with PTC (OR = 1.31, p = 2.20 × 10^-6 ). Rs1632250, Rs1863347 and rs1755787 showed association with classical PTC (cPTC) (n = 891; OR = 1.24, 2.22 × 10^-3 , OR = 1.31, p = 2.15 × 10^-4 and OR = 1.24, p = 2.06 × 10^-3 , respectively) while variant rs28397092 showed association with follicular variant (n = 243; OR = 1.51, p = 1.36 × 10^-3 ). Rs1863347 was associated with suppression of PTCSC3 in unaffected thyroid tissue (p = 0.026). Rs1632250, rs1863347 and rs1755787 showed association with multifocality (OR = 1.85, p = 0.001, OR = 1.98, p = 0.001 and OR = 1.76, p = 0.003 respectively) and N stage (OR = 1.79, p = 0.014, OR = 1.73, p = 0.023 and OR = 1.81, p = 0.013, respectively) in microPTC (n = 328) while rs368187 was associated with M stage (OR = 0.56, p = 0.034) in cPTC. Our results disclose multiple variants associated with PTC and clinical features in the 14q13 superlocus. We suggest that translational genotype/phenotype studies should take into account not only somatic mutations but also germline variants.

Genome-Wide Association Study Reveals Distinct Genetic Susceptibility of Thyroid Nodules From Thyroid Cancer

CONTEXT

Thyroid nodules are very common, and 7% to 15% of them are diagnosed as thyroid cancer. However, the inherited genetic risk factors for thyroid nodules and their associations with thyroid cancer remain unknown.

OBJECTIVE

To identify the genetic variants associated with susceptibility to thyroid nodules in comparison with thyroid cancer.

DESIGN AND SETTING

We performed a three-stage genome-wide association study for thyroid nodules. The discovery stage involved a genome-wide scan of 811 subjects with thyroid nodules and 691 subjects with a normal thyroid from a population-based cohort. Replication studies were conducted in an additional 1981 cases and 3100 controls from the participants of a health checkup. We also performed expression quantitative trait loci analysis of public data.

RESULTS

The most robust association was observed in TRPM3 (rs4745021) in the joint analysis (OR, 1.26; P = 6.12 × 10-8) and meta-analysis (OR, 1.28; P = 2.11 × 10-8). Signals at MBIP/NKX2-1 were replicated but did not reach genome-wide significance in the joint analysis (rs2415317, P = 4.62 × 10-5; rs944289, P = 8.68 × 10-5). The expression quantitative trait loci analysis showed that TRPM3 expression was associated with the rs4745021 genotype in thyroid tissues.

CONCLUSIONS

To the best of our knowledge, we have performed the first genome-wide association study of thyroid nodules and identified a susceptibility locus associated with thyroid nodules, suggesting that thyroid nodules have a genetic predisposition distinct from that of thyroid cancer.

Sex-specific genetic influence on thyroid-stimulating hormone and free thyroxine levels, and interactions between measurements: KNHANES 2013-2015

BACKGROUND

Although a wide range of genetic influences on thyroid-stimulating hormone (TSH) and free thyroxine (fT4) levels have been reported, sex differences in the genetic influences have not been well described.

METHODS

We assessed TSH and fT4 levels in 2,250 subjects without thyroid peroxidase antibody, with data obtained from the Korea National Health and Nutrition Examination Surveys (KNHANES) conducted from 2013 to 2015. Using variance decomposition methods, the variation of TSH and fT4 levels was divided into genetic and environmental components common to both sexes, and to males and females separately. The genetic correlation between TSH and fT4 levels was also assessed in both sexes, and in males and females separately.

RESULTS

Narrow-sense heritability for TSH and fT4 were 54% and 56%, respectively. Sex-specific heritability for TSH levels was significantly higher in females than in males (75% and 41%, respectively; p = 0.037). Heritability for fT4 levels was not significantly different between males and females (62% and 52%, respectively; p = 0.335). TSH and fT4 levels showed a negative genetic correlation in females (ρg = -0.347, p = 0.040) after regressing out the influences of environmental covariates, but this correlation was not present in males (ρg = -0.160, p = 0.391).

CONCLUSIONS

The genetic influences on individual TSH levels were more prominent in females than in males. In addition, female-specific pleiotropy between TSH and fT4 might be a clue that this stronger genetic influences in females would mainly affect thyroid function per se, rather than other TSH-related factors that do not primarily trigger the negative feedback loop between TSH and fT4.

The Association of Thyroid Function With Bone Density During Childhood

CONTEXT

Although the skeleton is a well-known thyroid hormone target organ, very little data are available on the association of thyroid function with bone outcomes during childhood.

OBJECTIVE

To study the association of thyroid function with bone mass during childhood.

DESIGN, SETTING, AND PARTICIPANTS

Population-based prospective cohort including 4204 children with TSH and free T4 (FT4) measured at the age of 6 years.

MAIN OUTCOME MEASURES

Bone density was assessed by a total body dual-energy X-ray absorptiometry scan at the median age of 6 years (95% range, 5.6 to 7.9) and at the age of 10 years (95% range, 9.0 to 10.9) in 4204 and 3404 participants, respectively.

RESULTS

There was an inverse association of TSH with bone mineral density (BMD) at the age of 6 (β -0.028 ± 0.011, P = 0.009) and with follow-up measurements at the age of 10 (β -0.027 ± 0.011, P = 0.014), but not with bone mineral content (BMC) at the age of 6 (β -0.028 ± 0.015, P = 0.06) or for follow-up measurements of BMC at the age of 10 (β -0.011 ± 0.015, P = 0.47). There was an inverse association of FT4 with BMD (β -0.016 ± 0.006, P = 0.014) and BMC (β -0.023 ± 0.009, P = 0.009) cross-sectionally, and also at the age of 10 years (BMD: β -0.018 ± 0.007, P = 0.007; BMC: β -0.021 ± 0.009, P = 0.020).

CONCLUSION

A higher FT4 concentration is associated with lower bone mass at the age of 6 and at the age of 10 years. These data provide insights into the effects of thyroid function on bone physiology during childhood.

Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves’ disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets.

Thyroid dysfunction is a common public health problem and associated with cardiovascular co-morbidities. Here, the authors carry out genome-wide meta-analysis for thyroid hormone (TH) levels, hyper- and hypothyroidism and identify SLC17A4 as a TH transporter and AADAT as a TH metabolizing enzyme.

Phenomic Impact of Genetically-Determined Euthyroid Function and Molecular Differences between Thyroid Disorders

BACKGROUND

The thyroid plays a key role in development and homeostasis, but it has been difficult to establish causality with diseases and phenotypic traits because of several potential confounders.

METHODS

To determine the causal effect of euthyroid function, we conducted a two-sample Mendelian randomization study of euthyroid thyrotropin (TSH) and free thyroxine (FT4) levels with respect to 2419 traits assessed in 337,199 individuals from UK Biobank. Additionally, we investigated the molecular differences between hypothyroidism and hyperthyroidism using genome-wide data.

RESULTS

After multiple testing correction, sixteen traits appear to be affected by genetically-determined euthyroid TSH, including multiple thyroid-related traits, e.g., hypothyroidism (p = 2.39 × 10-17), height (p = 2.76 × 10-10), body fat distribution (impedance of whole body, p = 4.43 × 10-8), pulse rate (p = 2.84 × 10-8), female infertility (p = 4.91 × 10-6), and hearing aid use (p = 7.10 × 10-5). Moreover, we found a consistent genetic correlation between hypothyroidism and hyperthyroidism (rg = 0.45, p = 5.45 × 10-6) with several immune pathways shared between these diseases. Two molecular pathways survived multiple testing correction for specificity to hyperthyroidism, JAK/STAT signaling (p = 1.02 × 10-6) and Rac guanyl-nucleotide exchange factor activity (p = 4.39 × 10-6).

CONCLUSION

Our data shed new light on the inter-individual variability of euthyroid function and the molecular mechanisms of the two thyroid disorders investigated.

Assessment of the genetic and clinical determinants of fracture risk: genome wide association and mendelian randomisation study

OBJECTIVES

To identify the genetic determinants of fracture risk and assess the role of 15 clinical risk factors on osteoporotic fracture risk.

DESIGN

Meta-analysis of genome wide association studies (GWAS) and a two-sample mendelian randomisation approach.

SETTING

25 cohorts from Europe, United States, east Asia, and Australia with genome wide genotyping and fracture data.

PARTICIPANTS

A discovery set of 37 857 fracture cases and 227 116 controls; with replication in up to 147 200 fracture cases and 150 085 controls. Fracture cases were defined as individuals (>18 years old) who had fractures at any skeletal site confirmed by medical, radiological, or questionnaire reports. Instrumental variable analyses were performed to estimate effects of 15 selected clinical risk factors for fracture in a two-sample mendelian randomisation framework, using the largest previously published GWAS meta-analysis of each risk factor.

RESULTS

Of 15 fracture associated loci identified, all were also associated with bone mineral density and mapped to genes clustering in pathways known to be critical to bone biology (eg, SOST, WNT16, and ESR1) or novel pathways (FAM210A, GRB10, and ETS2). Mendelian randomisation analyses showed a clear effect of bone mineral density on fracture risk. One standard deviation decrease in genetically determined bone mineral density of the femoral neck was associated with a 55% increase in fracture risk (odds ratio 1.55 (95% confidence interval 1.48 to 1.63; P=1.5×10-68). Hand grip strength was inversely associated with fracture risk, but this result was not significant after multiple testing correction. The remaining clinical risk factors (including vitamin D levels) showed no evidence for an effect on fracture.

CONCLUSIONS

This large scale GWAS meta-analysis for fracture identified 15 genetic determinants of fracture, all of which also influenced bone mineral density. Among the clinical risk factors for fracture assessed, only bone mineral density showed a major causal effect on fracture. Genetic predisposition to lower levels of vitamin D and estimated calcium intake from dairy sources were not associated with fracture risk.

Thyroid Stimulating Hormone and Bone Mineral Density: Evidence From a Two-Sample Mendelian Randomization Study and a Candidate Gene Association Study

With population aging, prevalence of low bone mineral density (BMD) and associated fracture risk are increased. To determine whether low circulating thyroid stimulating hormone (TSH) levels within the normal range are causally related to BMD, we conducted a two-sample Mendelian randomization (MR) study. Furthermore, we tested whether common genetic variants in the TSH receptor (TSHR) gene and genetic variants influencing expression of TSHR (expression quantitative trait loci [eQTLs]) are associated with BMD. For both analyses, we used summary-level data of genomewide association studies (GWASs) investigating BMD of the femoral neck (n = 32,735) and the lumbar spine (n = 28,498) in cohorts of European ancestry from the Genetic Factors of Osteoporosis (GEFOS) Consortium. For the MR study, we selected 20 genetic variants that were previously identified for circulating TSH levels in a GWAS meta-analysis (n = 26,420). All independent genetic instruments for TSH were combined in analyses for both femoral neck and lumbar spine BMD. In these studies, we found no evidence that a genetically determined 1-standard deviation (SD) decrease in circulating TSH concentration was associated with femoral neck BMD (0.003 SD decrease in BMD per SD decrease of TSH; 95% CI, -0.053 to 0.048; p = 0.92) or lumbar spine BMD (0.010 SD decrease in BMD per SD decrease of TSH; 95% CI, -0.069 to 0.049; p = 0.73). A total of 706 common genetic variants have been mapped to the TSHR locus and expression loci for TSHR. However, none of these genetic variants were associated with BMD at the femoral neck or lumbar spine. In conclusion, we found no evidence for a causal effect of circulating TSH on BMD, nor did we find any association between genetic variation at the TSHR locus or expression thereof and BMD. © 2018 The Authors. Journal of Bone and Mineral Research Published by WileyPeriodicals, Inc.

Empowering thyroid hormone research in human subjects using OMICs technologies

OMICs subsume different physiological layers including the genome, transcriptome, proteome and metabolome. Recent advances in analytical techniques allow for the exhaustive determination of biomolecules in all OMICs levels from less invasive human specimens such as blood and urine. Investigating OMICs in deeply characterized population-based or experimental studies has led to seminal improvement of our understanding of genetic determinants of thyroid function, identified putative thyroid hormone target genes and thyroid hormone-induced shifts in the plasma protein and metabolite content. Consequently, plasma biomolecules have been suggested as surrogates of tissue-specific action of thyroid hormones. This review provides a brief introduction to OMICs in thyroid research with a particular focus on metabolomics studies in humans elucidating the important role of thyroid hormones for whole body metabolism in adults.

Genome-Wide Association Studies of Autoimmune Thyroid Diseases, Thyroid Function, and Thyroid Cancer

Thyroid diseases, including autoimmune thyroid diseases and thyroid cancer, are known to have high heritability. Family and twin studies have indicated that genetics plays a major role in the development of thyroid diseases. Thyroid function, represented by thyroid stimulating hormone (TSH) and free thyroxine (T4), is also known to be partly genetically determined. Before the era of genome-wide association studies (GWAS), the ability to identify genes responsible for susceptibility to thyroid disease was limited. Over the past decade, GWAS have been used to identify genes involved in many complex diseases, including various phenotypes of the thyroid gland. In GWAS of autoimmune thyroid diseases, many susceptibility loci associated with autoimmunity (human leukocyte antigen [HLA], protein tyrosine phosphatase, non-receptor type 22 [PTPN22], cytotoxic T-lymphocyte associated protein 4 [CTLA4], and interleukin 2 receptor subunit alpha [IL2RA]) or thyroid-specific genes (thyroid stimulating hormone receptor [TSHR] and forkhead box E1 [FOXE1]) have been identified. Regarding thyroid function, many susceptibility loci for levels of TSH and free T4 have been identified through genome-wide analyses. In GWAS of differentiated thyroid cancer, associations at FOXE1, MAP3K12 binding inhibitory protein 1 (MBIP)-NK2 homeobox 1 (NKX2-1), disrupted in renal carcinoma 3 (DIRC3), neuregulin 1 (NRG1), and pecanex-like 2 (PCNXL2) have been commonly identified in people of European and Korean ancestry, and many other susceptibility loci have been found in specific populations. Through GWAS of various thyroid-related phenotypes, many susceptibility loci have been found, providing insights into the pathogenesis of thyroid diseases and disease co-clustering within families and individuals.

An Association Mapping Framework To Account for Potential Sex Difference in Genetic Architectures

Over the past few years, genome-wide association studies have identified many trait-associated loci that have different effects on females and males, which increased attention to the genetic architecture differences between the sexes. The between-sex differences in genetic architectures can cause a variety of phenomena such as differences in the effect sizes at trait-associated loci, differences in the magnitudes of polygenic background effects, and differences in the phenotypic variances. However, current association testing approaches for dealing with sex, such as including sex as a covariate, cannot fully account for these phenomena and can be suboptimal in statistical power. We present a novel association mapping framework, MetaSex, that can comprehensively account for the genetic architecture differences between the sexes. Through simulations and applications to real data, we show that our framework has superior performance than previous approaches in association mapping.

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.