Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Meta-Analysis
. 2023 Mar;33(3):301-311.
doi: 10.1089/thy.2022.0373. Epub 2023 Mar 1.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9

Affiliations
Meta-Analysis

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9

Antoine Weihs et al. Thyroid. 2023 Mar.

Abstract

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

Keywords: DNA methylation; KLF9; Mendelian randomization; gene expression; thyroid function.

PubMed Disclaimer

Conflict of interest statement

H.J.G. has received travel grants and speakers honoraria from Fresenius Medical Care, Neuraxpharm, Servier, and Janssen Cilag as well as research funding from Fresenius Medical Care not related to this work. R.E.M. has received speaker fees from Illumina and is an advisor to the Epigenetic Clock Development Foundation. All other authors declare no conflicts of interest.

Figures

FIG. 1.
FIG. 1.
Flowchart of the analyses performed in the current study. Image created with BioRender.com
FIG. 2.
FIG. 2.
Chicago plot showing the results of the TSH discovery analysis. Red labels show the significant sites (p-value <1.1E-07) of the discovery analysis, while green labels show results of the combined analysis. The black line shows the epigenome-wide significance cutoff (1.1E-07). Points plotted in the upper panel show a positive correlation with TSH, while those plotted below have a negative association. TSH, thyrotropin.
FIG. 3.
FIG. 3.
Chicago plot showing the results of the fT4 discovery analysis. Red labels show the significant sites (p-value <1.1E-07) of the discovery analysis, while green labels show results of the combined analysis. The black line shows the epigenome-wide significance cutoff (1.1E-07). Points plotted in the upper panel show a positive correlation with fT4, while those plotted below have a negative association. fT4, free thyroxine.
FIG. 4.
FIG. 4.
Chicago plot showing the results of the fT3 discovery analysis. Red labels show the significant sites (p-value <1.1E-07) of the discovery analysis, while green labels show results of the combined analysis. The black line shows the epigenome-wide significance cutoff (1.1E-07). Points plotted in the upper panel show a positive correlation with fT3, while those plotted below have a negative association. fT3, free triiodothyronine.
FIG. 5.
FIG. 5.
Overview of the study analyses and findings. Known mechanisms and associations are colored in black, and new findings revealed by this study are marked in blue color. Image created with BioRender.com

References

    1. Cooper DS, Biondi B. Subclinical thyroid disease. Lancet 2012;379(9821):1142–1154; doi: 10.1016/S0140-6736(11)60276-6 - DOI - PubMed
    1. Wondisford, Fredric E; Radovick S, (eds). Clinical Management of Thyroid Disease, 1st ed. Elsevier: Philadelphia, USA; 2009.
    1. Wu Y, Koenig RJ. Gene regulation by thyroid hormone. Trends Endocrinol Metab 2000;11(6):207–211; doi: 10.1016/S1043-2760(00)00263-0 - DOI - PubMed
    1. Kyono Y, Sachs LM, Bilesimo P, et al. . Developmental and thyroid hormone regulation of the DNA methyltransferase 3a gene in Xenopus tadpoles. Endocrinology 2016;157(12):4961–4972; doi: 10.1210/en.2016-1465 - DOI - PMC - PubMed
    1. Lafontaine N, Campbell PJ, Castillo-Fernandez JE, et al. . Epigenome-wide association study of thyroid function traits identifies novel associations of fT3 with KLF9 and DOT1L. J Clin Endocrinol Metab 2021;106(5):e2191–e2202; doi: 10.1210/clinem/dgaa975 - DOI - PMC - PubMed

Publication types

LinkOut - more resources