
==== Front
Clin Epigenetics
Clin Epigenetics
Clinical Epigenetics
1868-7075
1868-7083
BioMed Central London

1743
10.1186/s13148-024-01743-6
Research
Causal association of epigenetic age acceleration and risk of subacute thyroiditis: a bidirectional Mendelian randomization study
Shen Bingbing 12
Pu Yusheng 34
Zheng Xiaofeng 56
Liu Yang 7
Yang Lin 8
Liu Jiaye liujiaye925@wchscu.cn

129
Li Zhihui rockoliver@vip.sina.com

12
1 grid.412901.f 0000 0004 1770 1022 Division of Thyroid Surgery, Department of General Surgery, Laboratory of Thyroid and Parathyroid Diseases, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, No 37. Guoxue Alley, Chengdu, 610000 China
2 grid.412901.f 0000 0004 1770 1022 Department of Respiratory and Critical Care Medicine, Frontiers Science Center for Disease-Related Molecular Network, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, West China Hospital, Sichuan University, Chengdu, 610000 China
3 grid.412901.f 0000 0004 1770 1022 Center of Gerontology and Geriatrics, West China Hospital, Sichuan University, Chengdu, 610000 China
4 https://ror.org/011ashp19 grid.13291.38 0000 0001 0807 1581 West China School of Nursing, Sichuan University, Chengdu, 610000 China
5 grid.412901.f 0000 0004 1770 1022 Department of Endocrinology and Metabolism, West China Hospital, Sichuan University, Chengdu, China
6 grid.412901.f 0000 0004 1770 1022 Center for Diabetes and Metabolism Research, Division of Endocrinology and Metabolism, West China Hospital, Sichuan University, Chengdu, 610000 China
7 https://ror.org/017z00e58 grid.203458.8 0000 0000 8653 0555 Department of Obstetrics and Gynecology, The Second Affiliated Hospital, Chongqing Medical University, Chongqing, 400010 China
8 https://ror.org/017z00e58 grid.203458.8 0000 0000 8653 0555 Department of Gynecology, The First Affiliated Hospital, Chongqing Medical University, Chongqing, 400042 China
9 https://ror.org/011ashp19 grid.13291.38 0000 0001 0807 1581 Frontiers Medical Center, Tianfu Jincheng Laboratory, Sichuan University, Chengdu, 610000 China
19 9 2024
19 9 2024
2024
16 13328 4 2024
13 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Background

Epigenetic age accelerations (EAAs) are a promising new avenue of research, yet their investigation in subacute thyroiditis (SAT) remains scarce. Our study endeavors to fill this void by exploring the potential causal association between EAAs and SAT.

Methods

Our study utilized publicly available genome-wide association study (GWAS) data of European ancestry to conduct a bidirectional Mendelian randomization (MR) study. Five MR methods were employed to measure causal association between EAAs and SAT multiple analyses were utilized to perform quality control.

Results

Our study evaluated causal association between SAT and four EAAs, included GrimAge acceleration (GrimAA), Hannum age acceleration (HannumAA), PhenoAge acceleration (PhenoAA), intrinsic epigenetic age acceleration (IEAA). Results showed that there is a significant causal association between PhenoAA and SAT (OR 1.109, 95% CI 1.000–1.228, p = 0.049, by IVW method). On the contrary, SAT was associated with IEAA (OR 0.933, 95% CI 0.884–0.984, p = 0.011, by IVW method; OR 0.938, 95% CI 0.881–0.998, p = 0.043, by weighted median method). Leave-one-out sensitivity analysis, heterogeneity test, pleiotropy test, and MR-PRESSO analysis provide good quality control.

Conclusion

The bidirectional MR analysis concluded that an increase in PhenoAA was correlated with a higher risk of SAT, indicating a potential causal relationship between PhenoAA and risk of SAT. Conversely, SAT was found to be closely associated with IEAA, suggesting that SAT may accelerate the aging process. Slowing down biological aging has emerged as a new research direction in curbing SAT.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-024-01743-6.

Keywords

Subacute thyroiditis
Epigenetic age acceleration
Mendelian randomization
Age
The fellowship of China Postdoctoral Science Foundation2021M702340 Liu Jiaye The Postdoctoral research grant of Sichuan University2023SCU12047 Liu Jiaye The Science and Technology Department of Sichuan Province2021ZYCD016,2022NSFSC1441 Li Zhihui issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Subacute thyroiditis (SAT) is characterized by inflammation of the thyroid gland due to follicle damage, often accompanied by neck pain and tenderness in the affected area [1]. The majority of individuals affected by SAT are typically in the middle-aged bracket, with women exhibiting a significantly higher propensity, ranging from four to seven times more likely to develop the symptoms compared to men [2, 3]. Nonsteroidal anti-inflammatory medications (NSAIDs) and corticosteroids were utilized to alleviate patient symptoms and attenuate inflammatory reactions [4, 5]. Despite the utilization of medication in the treatment of SAT, the recurrence rate remains notably high, ranging from 1.6% to over 20%. Elevated recurrence rates and prolonged treatment durations have emerged as significant challenges in managing SAT [6, 7].

Aging is intricately linked to an increased risk of cancer, degenerative diseases, metabolic disorders, and cardiovascular conditions, influenced by a combination of environmental and genetic factors [8–12]. As is well known, physiological age measured over time does not precisely reflect the true aging status of the human body, including its tissues and organs. Individual aging is influenced by complex interactions of confounding factors and genetic element [13]. Epigenetic clocks emerge as the most promising predictors of biological age, capable of forecasting both chronological age and mortality [14]. DNA methylation (DNAm) at specific cytosine-phospho-guanine (CpG) sites was found to be strongly correlated with chronological age, so DNAm can be used to predict chronological age. Epigenetic clocks are heritable indicators of biological aging based on DNAm. Different epigenetic clocks based on different DNAm levels were used to assess the characteristics of epigenetic aging [15]. Epigenetic age acceleration (EAA), reflecting the actual variance between epigenetic age and chronological age is closely intertwined with numerous age-associated diseases and life expectancy across racial groups [13, 16]. At present, a series of methods were utilized to assess EAA. Intrinsic epigenetic age acceleration (IEAA) represents aging independent of blood cell type composition [17]. DNA methylation Hannum age acceleration (HannumAA) more accurately reflects extrinsic aging, which evaluated chronological age based on 71 age-related CpGs [18]. DNA methylation GrimAge acceleration (GrimAA) predicted chronological age by 1030 CpGs associated with smoking pack-years and 7 plasma proteins and DNA methylation PhenoAge acceleration (PhenoAA) are based on 513 CpGs associated with mortality and 9 clinical biomarkers. GrimAA and PhenoAA belong to the second generation and are designed to better predict age-associated diseases and mortality [10, 19]. Recent studies have unveiled correlations between epigenetic age acceleration and kidney function [20], as well as infertility [21], atrial fibrillation [11], and schizophrenia [22]. However, there is no study that has revealed correlations between epigenetic age acceleration and SAT.

Mendelian randomization (MR) is a causal inference method linking exposure and disease, employing genetic variants as instrumental variables (IVs) for exposure [23, 24]. To reduce systematic bias and provide results similar to RCTs, the MR approach uses genetic variations associated with exposure as IVs. By using this method, confounding variables within populations are lessened. A recent genome-wide association study (GWAS) meta-analysis disclosed 137 genetic loci related to EAAs, so can be applied to MR study [25]. To far, no research has used MR analysis to investigate the causal relationship between SAT and EAA. In this work, we investigate the causal relationship between the risk of SAT development and the epigenetic clock (IEAA, HunnumAA, GrimAA, PhenoAA) using two-sample MR analysis. For forward MR analysis, EAAs were applied with exposure, single-nucleotide polymorphisms (SNPs) closely associated with EAAs were considered as IVs, and SAT was outcome. For reverse MR analysis, SAT was used as exposure, SNPs associated with SAT were regarded as IVs, and EAAs were regarded as outcome. Our study was based on three assumptions: (1) the relevance assumption with genetic variants closely related to disease; (2) the independence assumption without confounding factors of disturbing casual association between exposure and outcome; and (3) the exclusion restriction assumption without ways for IVs influence outcome [26, 27]. Figure 1 displays MR assumptions and study design.Fig. 1 Study design of the correlation between EAAs and SAT. Assumption 1, the relevance assumption; Assumption 2, the independence assumption; Assumption 3, the exclusion restriction assumption. SAT, subacute thyroiditis; EAAs, epigenetic age accelerations; SNPs, single-nucleotide polymorphisms

Materials and methods

Data sources

The MR study strictly adheres to the STROBE-MR guidelines (Table S1) [28]. To observe casual association between EAAs and subacute thyroiditis, we performed two-sample MR analyses. Two-sample MR analyses identified causal association between exposure and outcome by utilizing genetic variants of exposure as IVs, which remedied shortcoming of observational studies and utilized publicly GWAS datasets of exposure (as a risk factor) and outcome (as a disease).

To reduce bias from population stratification, GWAS meta-analyses of European ancestry were utilized to our study. GWAS data of epigenetic clocks were originated from GWAS database (https://gwas.mrcieu.ac.uk/). We selected most-recent largescale GWAS meta-analyses of EAAs, including IEAA, HannumAA, GrimAA and PhenoAA [25]. GWAS data of subacute thyroiditis were originated from FINNGEN database (https://www.finngen.fi/en/access_results). GWAS summary statistics of subacute thyroiditis were used to our study (finn-b-E4_THYROIDITSUBAC, Table S2) [29].

Selection of instrument variables

Appropriate IVs were used in order to guarantee the correctness of MR findings. First, based on the criteria (p < 5e−08), single-nucleotide polymorphisms (SNPs) connected to the exposure were selected as IVs. Nevertheless, we discovered that only four SNPs were chosen for GrimAA exposure in European ancestry, and six SNPs were chosen for exposure for subacute thyroiditis in FINNGEN. Therefore, in order to obtain a comparatively large number of IVs, we gently adjusted the threshold (p < 5e−06). Furthermore, the clump_data function is used to exclude the linkage disequilibrium (LD) and sets the kilobase pairs (kb) to 10,000. The r2 threshold is set to 0.001 (Tables S3 and S4). In addition, we calculated F value through evaluate the intensity of IVs by utilizing the formula beta2/SE2 with F value > 10 were selected [26, 30] (Tables S5 and S6). The formula for calculating the proportion of variation in traits explained by genetic instruments is R2 = 2 × MAF × (1-MAF) × beta, where MAF stands for minor allele frequency, beta for effect size, SE for standard error, N for sample size, and k for the number of IVs [31]. To verify whether selected IVs possessed the independence assumption, LDLink database (https://ldlink.nih.gov/?tab=ldtrait) was applied to exclude SNPs, which may be associated with outcome [32–34]. SNPs linked to inflammation and infection, such as immune cell counts (neutrophils, monocytes, white blood cells) and oral ulcers, are eliminated for EAA exposure based on the GWAS characteristic (Table S7). On the other hand, SNPs linked to age and smoking for SAT as exposure (Table S8). SNPs are covered in depth in Tables S9, S10, S11, S12 and S13, with GWAS trait and source citations.

Statistical analysis

In our study, we use five MR methods to perform study. The IVW method served as the primary analytical approach, while others were used as supplementary methods. If all SNPs are deemed valid IVs, then IVW can yield precise estimations [35]. MR Egger can detect and correct for pleiotropy but might have lower precision [36]. Weighted median can offer accurate estimates if over 50% of SNPs are deemed valid IVs [37]. Though simple mode might not wield the same statistical power as IVW, it retains the ability to detect pleiotropy [38]. The selection of the optimal bandwidth for mode estimation can influence weighted mode [39]. Therefore, we primarily rely on the IVW method.

Furthermore, we performed a series of tests to obtain robust and stable results. Heterogeneity test assess heterogeneity of IVs through calculating Q statistic [40]. The MR-PRESSO analysis was employed to exclude SNPs with a p-value less than 0.05. Additionally, any outliers detected were removed from the analysis. In addition, pleiotropy test was used to evaluate horizontal pleiotropy through calculating MR Egger intercept values [41]. Leave-one-out sensitivity analysis was employed to analyze each SNPs, which verify total effect of SNPs after removing one of the SNPs.

All analyses in this study were performed using R software (version 4.3.3). We mainly used “TwoSampleMR” package. All statistical tests were conducted with a two-sided significance level set at p < 0.05.

Results

Selection of instrument variables

To enhance the feasibility of obtaining robust MR results, we strictly selected appropriate SNPs as IVs to examine the relationship between EAA and SAT. For EAAs as the exposure and subacute thyroiditis as the outcome, we selected 26 SNPs for IEAA, 14 SNPs for PhenoAA, 15 SNPs for GrimAA, and 21 SNPs for HannumAA. Conversely, for subacute thyroiditis as the exposure and EAAs as the outcome, 5 SNPs for subacute thyroiditis were chosen. Tables S7 and S8 provide detailed information on these IVs.

Causal analysis of epigenetic age accelerations on subacute thyroiditis

Through MR analyses based on five methods, we found that a partial causal relationship between EAAs and SAT. IVW method revealed a significant causal association between PhenoAA and SAT (OR 1.109, 95% CI 1.000 to 1.228, p = 0.049). However, there are no significant causal association based on other four methods (OR 1.122, 95% CI 0.905–1.390, p = 0.316 by MR Egger; OR 1.141, 95% CI 0.999–1.304, p = 0.052 by weighted median; OR 1.149, 95% CI 0.943–1.400, p = 0.190 by simple mode; OR 1.147, 95% CI 0.962–1.368, p = 0.141 by weighted mode) (Table 1, Table S14). Thus, we deduced that there was a strong correlation between PhenoAA and the elevated risk of SAT. Furthermore, based on five MR methodologies, our data indicate that there is not a significant relationship between SAT and other EAAs, such as GrimAA, HannumAA, and IEAA (Table 1, Table S14). Scatter plot displayed the significant causal association of EAAs on SAT (Fig. 2). Forest plots of SNPs are displayed in Fig. S2.Table 1 The causal association of EAAs on SAT based on five MR methods

Exposure	Outcome	Method	nSNP	OR (95%CI)	p-value	
GrimAA	SAT	MR Egger	15	1.505 (0.866, 2.617)	0.171	
Weighted median	15	1.090 (0.890, 1.335)	0.401	
IVW	15	1.092 (0.944, 1.264)	0.237	
Simple mode	15	0.942 (0.662, 1.340)	0.744	
Weighted mode	15	0.958 (0.690, 1.330)	0.800	
HannumAA	SAT	MR Egger	21	0.882 (0.614, 1.267)	0.504	
Weighted median	21	1.039 (0.881, 1.225)	0.650	
IVW	21	0.963 (0.856, 1.085)	0.539	
Simple mode	21	1.051 (0.801, 1.378)	0.724	
Weighted mode	21	1.073 (0.850, 1.353)	0.559	
IEAA	SAT	MR Egger	26	1.013 (0.855, 1.200)	0.884	
Weighted median	26	1.014 (0.903, 1.138)	0.817	
IVW	26	1.035 (0.948, 1.130)	0.443	
Simple mode	26	1.057 (0.866, 1.292)	0.589	
Weighted mode	26	1.013 (0.882, 1.163)	0.861	
PhenoAA	SAT	MR Egger	14	1.122 (0.905, 1.390)	0.316	
Weighted median	14	1.141 (0.999, 1.304)	0.052	
IVW	14	1.109 (1.000, 1.228)	0.049	
Simple mode	14	1.149 (0.943, 1.400)	0.190	
Weighted mode	14	1.147 (0.962, 1.368)	0.151	
SAT, subacute thyroiditis; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; SNPs, single-nucleotide polymorphisms; MR, Mendelian randomization; IVW, inverse variance weighted; OR, odds ratio; CI, confidence interval

Fig. 2 Scatter plots of casual relationship evaluate for EAAs on SAT. A, GrimAA on SAT. B, HannumAA on SAT. C, IEAA on SAT. D, PhenoAA on SAT. Contains 15, 21, 26, and 14 IVs, respectively. SAT, subacute thyroiditis; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; IVs, instrumental variables

There was no significant heterogeneity among SNPs by heterogeneity test based on IVW and MR Egger methods (GrimAA p = 0.868, and Q = 8.386 by IVW, p = 0.902 and Q = 6.995 by MR Egger; HannumAA p = 0.337, and Q = 22.073 by IVW, p = 0.295 and Q = 21.777 by MR Egger; PhenoAA p = 0.205, and Q = 30.536 by IVW, p = 0.171 and Q = 30.428 by MR Egger; IEAA p = 0.393, and Q = 13.732 by IVW, p = 0.319 and Q = 13.715 by MR Egger, Table 2). Leave-one-out sensitivity analysis revealed that after deleting SNP one by one, there are no changes of the total effect (Fig. S1). MR-PRESSO global test was utilized to check detect the presence of horizontal pleiotropy and MR-PRESSO outlier test was used to delete SNPs with horizontal pleiotropy. Results show that there is no horizontal pleiotropy in there SNPs (Table 2, Table S15). Funnel plots further revealed that there is no heterogeneity of SNPs (Fig. S3). MR Egger intercept of pleiotropy test shows the absence of horizontal pleiotropy (Table 2, Table S16). These results imply that our findings from the MR results exhibit substantial confidence, demonstrating robustness and consistency.Table 2 Sensitivity analyses of EAAs on SAT

Exposure	Outcome	Heterogeneity	Pleiotropy
MR Egger intercept
(p-value)	MR-PRESSO
(p-value)	
MR egger	IVW	
Q (p-value)	
GrimAA	SAT	6.995 (0.902)	8.386 (0.868)	− 0.065 (0.259)		
HannumAA	SAT	21.777 (0.295)	22.073 (0.337)	0.023 (0.617)		
IEAA	SAT	30.428 (0.171)	30.536 (0.205)	0.007 (0.773)		
PhenoAA	SAT	13.715 (0.319)	13.732 (0.393)	− 0.004 (0.905)	15.591 (0.452)	
SAT, subacute thyroiditis; MR-PRESSO, MR pleiotropy residual sum and outlier; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; SNPs, single-nucleotide polymorphisms; MR, Mendelian randomization; IVW, inverse variance weighted

Casual association of subacute thyroiditis on epigenetic age accelerations

To analyze reverse casual association, we performed a series of analyses by using subacute thyroiditis as exposure. Results show that SAT was associated with IEAA (OR, 0.933, 95% CI 0.884–0.984, p = 0.011, by IVW method; OR 0.938, 95% CI 0.881–0.998, p = 0.043, by weighted median) and the OR values of other methods are in the same direction (OR 0.967, 95% CI 0.877–1.066, p = 0.544 by MR Egger; OR 0.940, 95% CI 0.849–1.040, p = 0.296 by simple mode; OR 0.40, 95% CI 0.878–1.007, p = 0.151 by weighted mode) (Table 3, Table S17). Scatter plot displayed the significant causal association of SAT on EAAs (Fig. 3). Thus, we concluded that SAT was associated with decreased IEAA. Based on five methods, we discovered no causal relationship between SAT and other EAAs (including GrimAA, HannumAA and PhenoAA). Forest plots of SNPs are displayed in Fig. S5. Leave-one-out sensitivity analysis also presents higher the total effect after removing the SNPs one by one (Fig. S4). Heterogeneity test and pleiotropy test did not find heterogeneity and pleiotropy (Table 4, Tables S18–S19). MR-PRESSO analysis did not found SNPs to be excluded in relation to the association between SAT and IEAA (Table 4). Funnel plots further revealed that there is no heterogeneity of SNPs (Fig. S6).Table 3 The causal association of SAT on EAAs based on five MR methods

Exposure	Outcome	Method	nSNP	OR (95%CI)	p-value	
SAT	GrimAA	MR Egger	5	0.924 (0.839, 1.017)	0.206	
Weighted median	5	1.006 (0.942, 1.074)	0.855	
IVW	5	1.026 (0.946, 1.112)	0.535	
Simple mode	5	1.110 (0.926, 1.320)	0.328	
Weighted mode	5	0.968 (0.908, 1.032)	0.378	
SAT	HannumAA	MR Egger	5	0.982 (0.833, 1.157)	0.840	
Weighted median	5	0.982 (0.881, 1.225)	0.560	
IVW	5	0.975 (0.901, 1.055)	0.529	
Simple mode	5	0.997 (0.898, 1.106)	0.954	
Weighted mode	5	0.979 (0.918, 1.044)	0.549	
SAT	IEAA	MR Egger	5	0.967 (0.877, 1.066)	0.544	
Weighted median	5	0.938 (0.881, 0.998)	0.042	
IVW	5	0.933 (0.884, 0.984)	0.011	
Simple mode	5	0.940 (0.849, 1.010)	0.296	
Weighted mode	5	0.940 (0.878, 1.006)	0.151	
SAT	PhenoAA	MR Egger	5	0.923 (0.816, 1.044)	0.292	
Weighted median	5	0.944 (0.870, 1.024)	0.163	
IVW	5	0.957 (0.894, 1.024)	0.207	
Simple mode	5	0.943 (0.837, 1.062)	0.387	
Weighted mode	5	0.941 (0.865, 1.023)	0.229	
SAT, subacute thyroiditis; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; SNPs, single-nucleotide polymorphisms; MR, Mendelian randomization; IVW, inverse variance weighted; OR, odds ratio; CI, confidence interval

Fig. 3 Scatter plots of casual relationship evaluate for SAT on EAAs. A, GrimAA on SAT. B, HannumAA on SAT. C, IEAA on SAT. D, PhenoAA on SAT. It contains 15, 21, 26, and 14 IVs, respectively. SAT, subacute thyroiditis; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; IVs, instrumental variables

Table 4 Sensitivity analyses of SAT on EAAs

Exposure	Outcome	Heterogeneity	Pleiotropy
MR Egger intercept
(p-value)	MR-PRESSO RSSobs
(p-value)	
MR egger	IVW	
Q (p-value)	
SAT	GrimAA	2.782 (0.427)	9.317 (0.054)	0.087 (0.083)		
SAT	HannumAA	8.913 (0.030)	8.941 (0.063)	− 0.006 (0.928)		
SAT	IEAA	2.677 (0.444)	3.425 (0.489)	− 0.030 (0.451)	6.260 (0.562)	
SAT	PhenoAA	1.297 (0.730)	1.782 (0.775)	− 0.030 (0.536)		
SAT, subacute thyroiditis; MR-PRESSO, MR pleiotropy residual sum and outlier; EAAs, epigenetic age accelerations; IEAA, intrinsic epigenetic age acceleration; HunnumAA, DNA methylation Hannum age acceleration; GrimAA, DNA methylation GrimAge acceleration; PhenoAA, DNA methylation PhenoAge acceleration; SNPs, single-nucleotide polymorphisms; MR, Mendelian randomization; IVW, inverse variance weighted

Discussion

In the current study, we utilized GWAS datasets to investigate the association between EAAs and SAT using MR analysis. Our study revealed that there is relationship between EAAs and SAT, increase in PhenoAA was correlated with a higher risk of SAT and SAT was found to be closely associated with decreased IEAA. Health aging possessed obvious heterogeneity, so it is necessary that search higher predictor to cognize and measure aging [42]. Recently, epigenetic clock and telomere length were considered as the most remarkable predictor in biological aging [14]. Previous study shows telomere length was correlated with valve stenosis (AVS) but the predictive ability is lower [43, 44]. Comparing with chronological age and other age-associated biomarkers, epigenetic clock is best predictor of biological aging status, which reflects biological aging status by measuring DNA methylation levels [13, 14, 17, 18]. Pan et al. identified that increased levels of IEAA and GrimAA were associated with elevated risk of chronic kidney disease (CKD) and a decline in kidney function by MR analysis. Furthermore, it was observed that accelerated HannumAA, GrimAA, and PhenoAA may lead to a reduction in log-transformed estimated glomerular filtration rate (log-eGFR). Specifically, for each 1-unit rise in log-eGFR, HannumAA decreased by 3.14 years, GrimAA decreased by 1.99 years, and PhenoAA decreased by 2.88 years, respectively [20]. Roberts JD found that EAA measures of GrimAA and PhenoAA were closely related to higher risk of atrial fibrillation (AF) [11]. Pan et al. concluded that PhenoAA was link to the risk of aortic valve stenosis (AVS) [45].

Nevertheless, no research has looked at the connection between SAT and EAAs. Our research is the first to show a causal association between EAAs and SAT as well as a genetic correlation. There is a strong causal relationship between SAT and EAAs. Thus, a new line of inquiry for SAT prevention is to slow down biological aging. Furthermore, our research revealed an extremely interesting finding that SAT is associated with lower IEAA. Reduced EAA was generally thought to be a beneficial protective factor. Decreased HannumAA, GrimAA, and PhenoAA are connected with increased log-eGFR, which indicated a progressive improvement in kidney function [20]. A recent MR study about obesity and EAAs revealed that EAAs was related to lower risk of obesity and obesity increased EAAs (including GrimAA and PhenoAA) [46], increasing age, and decreasing immunity [47]. The incidence of SAT is closely related to the low immunity [48], which interpreted our result that increased EAAs increase the risk of SAT. SAT is characterized by inflammation of the thyroid gland [1]. Previous study concluded that inflammation is not necessarily harmful to the body, triggering an inflammatory response by releasing inflammatory mediators to destroy pathogens and repair damaged tissues [49]. Therefore, a moderate inflammatory response is necessary to maintain human health. Our study uncovered that SAT related to decreased IEAA. We speculated that SAT may trigger inflammatory response to resist aging, which is a good aspect for SAT patients.

Our study explored causal association between EAAs and SAT. To build reliable causal inferences and overcome influences of confounding and reverse causality biases, we employed basic theory of random allocation of alleles by MR analysis [50]. In order to reduce the potential bias caused by population differences, we only used GWAS data of European ancestry in this study. In addition, rigorous quality evaluation shows that our results are reliable. Our study first reported that there are association between EAAs and SAT. The epigenetic clock may become surveillance indicator to evaluate risk of SAT for clinicians and preventive medicine practitioners, sowing down biological aging has become a new research direction for SAT.

Certainly, our study also has some limitations. Firstly, we only selected European ancestry, which led to concerns that our results would not apply to other populations. Secondly, when selecting SNPs, we used a standard bioinformatic threshold of p < 5 × 10−8, only 4 SNPs were selected for exposure for GrimAA in European-ancestry and 6 SNPs were selected for exposure for subacute thyroiditis in FINNGEN. Consequently, we chose to select SNPs under the threshold of p < 5 × 10−6, as previously employed in other studies [51–54]. Of course, this approach may introduce bias from weak instrumental variables. To mitigate this limitation, we used the F-test to seek strong evidence. Lastly, epigenetic aging is essential associated with environmental exposures rather than genetic factors, so it is significance to interpret results.

Conclusion

In conclusion, the bidirectional MR analysis revealed that an increase in PhenoAA was correlated with a higher risk of SAT, indicating a potential causal relationship between PhenoAA and risk of SAT. Conversely, SAT was found to be closely associated with IEAA, suggesting that SAT may accelerate the aging process. Slowing down biological aging has emerged as a new research direction in curbing SAT.

Supplementary Information

Additional file1 (DOCX 45 KB)

Additional file2 (DOCX 1561 KB)

Additional file3 (XLSX 248 KB)

Additional file4 (XLSX 141 KB)

Additional file5 (XLSX 129 KB)

Additional file6 (XLSX 267 KB)

Additional file7 (XLSX 76 KB)

Acknowledgements

We want to acknowledge the participants and investigators of GWAS and FinnGen study.

Author contributions

ZHL, BBS, YSP and JYL conceptualized and designed the study. BBS and JYL provided the “TwoSampleMR” package codes in R software and analyzed the data in the study. BBS drafted the manuscript. YL and LY download data of GWAS. JYL and XFZ gave constructive suggestions when writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the fellowship of China Postdoctoral Science Foundation (2021M702340), the Science and Technology Department of Sichuan Province (2021ZYCD016, 2022NSFSC1441) and the Postdoctoral research grant of Sichuan University (2023SCU12047).

Availability of data and materials

No datasets were generated or analyzed during the current study.

Declarations

Ethics approval and consent to participate

All data for the study were derived from the GWAS database, which has been reviewed and approved by the relevant ethics committee.

Consent for publication

Not applicable.

Conflict of interest

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Bingbing Shen and Yusheng Pu contributed equally to this article.
==== Refs
References

1. Desailloud R Hober D Viruses and thyroiditis: an update Virol J 2009 6 5 10.1186/1743-422X-6-5 19138419
Desailloud R, Hober D. Viruses and thyroiditis: an update. Virol J. 2009;6:5.19138419
2. Samuels MH Subacute, silent, and postpartum thyroiditis Med Clin North Am 2012 96 2 223 233 10.1016/j.mcna.2012.01.003 22443972
Samuels MH. Subacute, silent, and postpartum thyroiditis. Med Clin North Am. 2012;96(2):223–33.22443972
3. Stasiak M Michalak R Stasiak B Lewinski A Clinical characteristics of subacute thyroiditis is different than it used to be-current state based on 15 years own material Neuro Endocrinol Lett 2019 39 7 489 495 30860680
Stasiak M, Michalak R, Stasiak B, Lewinski A. Clinical characteristics of subacute thyroiditis is different than it used to be-current state based on 15 years own material. Neuro Endocrinol Lett. 2019;39(7):489–95.30860680
4. Ross DS Burch HB Cooper DS 2016 American thyroid association guidelines for diagnosis and management of hyperthyroidism and other causes of thyrotoxicosis Thyroid 2016 26 10 1343 1421 10.1089/thy.2016.0229 27521067
Ross DS, Burch HB, Cooper DS, et al. 2016 American thyroid association guidelines for diagnosis and management of hyperthyroidism and other causes of thyrotoxicosis. Thyroid. 2016;26(10):1343–421.27521067
5. Sencar ME Calapkulu M Sakiz D An evaluation of the results of the steroid and non-steroidal anti-inflammatory drug treatments in subacute thyroiditis in relation to persistent hypothyroidism and recurrence Sci Rep 2019 9 1 16899 10.1038/s41598-019-53475-w 31729433
Sencar ME, Calapkulu M, Sakiz D, et al. An evaluation of the results of the steroid and non-steroidal anti-inflammatory drug treatments in subacute thyroiditis in relation to persistent hypothyroidism and recurrence. Sci Rep. 2019;9(1):16899.31729433
6. Stasiak M Tymoniuk B Stasiak B Lewiński A The risk of recurrence of subacute thyroiditis is HLA-dependent Int J Mol Sci 2019 20 5 1089 10.3390/ijms20051089 30832406
Stasiak M, Tymoniuk B, Stasiak B, Lewiński A. The risk of recurrence of subacute thyroiditis is HLA-dependent. Int J Mol Sci. 2019;20(5):1089.30832406
7. Zhang J Ding G Li J Risk factors for subacute thyroiditis recurrence: a systematic review and meta-analysis of cohort studies Front Endocrinol (Lausanne) 2021 12 783439 10.3389/fendo.2021.783439 35002966
Zhang J, Ding G, Li J, et al. Risk factors for subacute thyroiditis recurrence: a systematic review and meta-analysis of cohort studies. Front Endocrinol (Lausanne). 2021;12:783439.35002966
8. Niccoli T Partridge L Ageing as a risk factor for disease Curr Biol 2012 22 17 R741 752 10.1016/j.cub.2012.07.024 22975005
Niccoli T, Partridge L. Ageing as a risk factor for disease. Curr Biol. 2012;22(17):R741-752.22975005
9. Dartigues JF Le Bourdonnec K Tabue-Teguo M Co-Occurrence of geriatric syndromes and diseases in the general population: assessment of the dimensions of aging J Nutr Health Aging 2022 26 1 37 45 10.1007/s12603-021-1722-3 35067701
Dartigues JF, Le Bourdonnec K, Tabue-Teguo M, et al. Co-Occurrence of geriatric syndromes and diseases in the general population: assessment of the dimensions of aging. J Nutr Health Aging. 2022;26(1):37–45.35067701
10. Levine ME Lu AT Quach A An epigenetic biomarker of aging for lifespan and healthspan Aging (Albany NY) 2018 10 4 573 591 10.18632/aging.101414 29676998
Levine ME, Lu AT, Quach A, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10(4):573–91.29676998
11. Roberts JD Vittinghoff E Lu AT Epigenetic age and the risk of incident atrial fibrillation Circulation 2021 144 24 1899 1911 10.1161/CIRCULATIONAHA.121.056456 34587750
Roberts JD, Vittinghoff E, Lu AT, et al. Epigenetic age and the risk of incident atrial fibrillation. Circulation. 2021;144(24):1899–911.34587750
12. Zhang F Deng S Zhang J Causality between heart failure and epigenetic age: a bidirectional Mendelian randomization study ESC Heart Fail 2023 10 5 2903 2913 10.1002/ehf2.14446 37452462
Zhang F, Deng S, Zhang J, et al. Causality between heart failure and epigenetic age: a bidirectional Mendelian randomization study. ESC Heart Fail. 2023;10(5):2903–13.37452462
13. Horvath S Raj K DNA methylation-based biomarkers and the epigenetic clock theory of ageing Nat Rev Genet 2018 19 6 371 384 10.1038/s41576-018-0004-3 29643443
Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat Rev Genet. 2018;19(6):371–84.29643443
14. Jylhävä J Pedersen NL Hägg S Biological age predictors EBioMedicine 2017 21 29 36 10.1016/j.ebiom.2017.03.046 28396265
Jylhävä J, Pedersen NL, Hägg S. Biological age predictors. EBioMedicine. 2017;21:29–36.28396265
15. Liu Z Leung D Thrush K Underlying features of epigenetic aging clocks in vivo and in vitro Aging Cell 2020 19 10 e13229 10.1111/acel.13229 32930491
Liu Z, Leung D, Thrush K, et al. Underlying features of epigenetic aging clocks in vivo and in vitro. Aging Cell. 2020;19(10):e13229.32930491
16. Chen BH Marioni RE Colicino E DNA methylation-based measures of biological age: meta-analysis predicting time to death Aging (Albany NY) 2016 8 9 1844 1865 10.18632/aging.101020 27690265
Chen BH, Marioni RE, Colicino E, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY). 2016;8(9):1844–65.27690265
17. Horvath S DNA methylation age of human tissues and cell types Genome Biol 2013 14 10 R115 10.1186/gb-2013-14-10-r115 24138928
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.24138928
18. Hannum G Guinney J Zhao L Genome-wide methylation profiles reveal quantitative views of human aging rates Mol Cell 2013 49 2 359 367 10.1016/j.molcel.2012.10.016 23177740
Hannum G, Guinney J, Zhao L, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67.23177740
19. Lu AT Quach A Wilson JG DNA methylation GrimAge strongly predicts lifespan and healthspan Aging (Albany NY) 2019 11 2 303 327 10.18632/aging.101684 30669119
Lu AT, Quach A, Wilson JG, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY). 2019;11(2):303–27.30669119
20. Pan Y Sun X Huang Z Effects of epigenetic age acceleration on kidney function: a Mendelian randomization study Clin Epigenetics 2023 15 1 61 10.1186/s13148-023-01476-y 37031184
Pan Y, Sun X, Huang Z, et al. Effects of epigenetic age acceleration on kidney function: a Mendelian randomization study. Clin Epigenetics. 2023;15(1):61.37031184
21. Lee Y Bohlin J Page CM Associations between epigenetic age acceleration and infertility Hum Reprod 2022 37 9 2063 2074 10.1093/humrep/deac147 35771672
Lee Y, Bohlin J, Page CM, et al. Associations between epigenetic age acceleration and infertility. Hum Reprod. 2022;37(9):2063–74.35771672
22. Wu X Ye J Wang Z Zhao C Epigenetic age acceleration was delayed in schizophrenia Schizophr Bull 2021 47 3 803 811 10.1093/schbul/sbaa164 33269797
Wu X, Ye J, Wang Z, Zhao C. Epigenetic age acceleration was delayed in schizophrenia. Schizophr Bull. 2021;47(3):803–11.33269797
23. Smith GD Ebrahim S 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol 2003 32 1 1 22 10.1093/ije/dyg070 12689998
Smith GD, Ebrahim S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.12689998
24. Boef AG Dekkers OM le Cessie S Mendelian randomization studies: a review of the approaches used and the quality of reporting Int J Epidemiol 2015 44 2 496 511 10.1093/ije/dyv071 25953784
Boef AG, Dekkers OM, le Cessie S. Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int J Epidemiol. 2015;44(2):496–511.25953784
25. McCartney DL Min JL Richmond RC Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging Genome Biol 2021 22 1 194 10.1186/s13059-021-02398-9 34187551
McCartney DL, Min JL, Richmond RC, et al. Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging. Genome Biol. 2021;22(1):194.34187551
26. Davey Smith G Hemani G Mendelian randomization: genetic anchors for causal inference in epidemiological studies Hum Mol Genet 2014 23 R1 R89 98 10.1093/hmg/ddu328 25064373
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23(R1):R89-98.25064373
27. Bowden J Davey Smith G Burgess S Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression Int J Epidemiol 2015 44 2 512 525 10.1093/ije/dyv080 26050253
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.26050253
28. Skrivankova VW Richmond RC Woolf BAR Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration BMJ 2021 375 n2233 10.1136/bmj.n2233 34702754
Skrivankova VW, Richmond RC, Woolf BAR, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ. 2021;375:n2233.34702754
29. Kurki MI Karjalainen J Palta P FinnGen provides genetic insights from a well-phenotyped isolated population Nature 2023 613 7944 508 518 10.1038/s41586-022-05473-8 36653562
Kurki MI, Karjalainen J, Palta P, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613(7944):508–18.36653562
30. Palmer TM Lawlor DA Harbord RM Using multiple genetic variants as instrumental variables for modifiable risk factors Stat Methods Med Res 2012 21 3 223 242 10.1177/0962280210394459 21216802
Palmer TM, Lawlor DA, Harbord RM, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21(3):223–42.21216802
31. Pierce BL Ahsan H Vanderweele TJ Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants Int J Epidemiol 2011 40 3 740 752 10.1093/ije/dyq151 20813862
Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. 2011;40(3):740–52.20813862
32. Machiela MJ Chanock SJ LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants Bioinformatics 2015 31 21 3555 3557 10.1093/bioinformatics/btv402 26139635
Machiela MJ, Chanock SJ. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics. 2015;31(21):3555–7.26139635
33. Lin SH Brown DW Machiela MJ LDtrait: an online tool for identifying published phenotype associations in linkage disequilibrium Cancer Res 2020 80 16 3443 3446 10.1158/0008-5472.CAN-20-0985 32606005
Lin SH, Brown DW, Machiela MJ. LDtrait: an online tool for identifying published phenotype associations in linkage disequilibrium. Cancer Res. 2020;80(16):3443–6.32606005
34. Lin SH Thakur R Machiela MJ LDexpress: an online tool for integrating population-specific linkage disequilibrium patterns with tissue-specific expression data BMC Bioinform 2021 22 1 608 10.1186/s12859-021-04531-8
Lin SH, Thakur R, Machiela MJ. LDexpress: an online tool for integrating population-specific linkage disequilibrium patterns with tissue-specific expression data. BMC Bioinform. 2021;22(1):608.
35. Yuan S Larsson S Causal associations of iron status with gout and rheumatoid arthritis, but not with inflammatory bowel disease Clin Nutr 2020 39 10 3119 3124 10.1016/j.clnu.2020.01.019 32044136
Yuan S, Larsson S. Causal associations of iron status with gout and rheumatoid arthritis, but not with inflammatory bowel disease. Clin Nutr. 2020;39(10):3119–24.32044136
36. Burgess S Thompson SG Interpreting findings from Mendelian randomization using the MR-Egger method Eur J Epidemiol 2017 32 5 377 389 10.1007/s10654-017-0255-x 28527048
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.28527048
37. Bowden J Davey Smith G Haycock PC Burgess S Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator Genet Epidemiol 2016 40 4 304 314 10.1002/gepi.21965 27061298
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.27061298
38. Milne RL Kuchenbaecker KB Michailidou K Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer Nat Genet 2017 49 12 1767 1778 10.1038/ng.3785 29058716
Milne RL, Kuchenbaecker KB, Michailidou K, et al. Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer. Nat Genet. 2017;49(12):1767–78.29058716
39. Hartwig FP Davey Smith G Bowden J Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption Int J Epidemiol 2017 46 6 1985 1998 10.1093/ije/dyx102 29040600
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.29040600
40. Greco MF Minelli C Sheehan NA Thompson JR Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome Stat Med 2015 34 21 2926 2940 10.1002/sim.6522 25950993
Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34(21):2926–40.25950993
41. Verbanck M Chen CY Neale B Do R Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases Nat Genet 2018 50 5 693 698 10.1038/s41588-018-0099-7 29686387
Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.29686387
42. Lowsky DJ Olshansky SJ Bhattacharya J Goldman DP Heterogeneity in healthy aging J Gerontol A Biol Sci Med Sci 2014 69 6 640 649 10.1093/gerona/glt162 24249734
Lowsky DJ, Olshansky SJ, Bhattacharya J, Goldman DP. Heterogeneity in healthy aging. J Gerontol A Biol Sci Med Sci. 2014;69(6):640–9.24249734
43. Kurz DJ Kloeckener-Gruissem B Akhmedov A Degenerative aortic valve stenosis, but not coronary disease, is associated with shorter telomere length in the elderly Arterioscler Thromb Vasc Biol 2006 26 6 e114 117 10.1161/atvb.26.6.1303 16627805
Kurz DJ, Kloeckener-Gruissem B, Akhmedov A, et al. Degenerative aortic valve stenosis, but not coronary disease, is associated with shorter telomere length in the elderly. Arterioscler Thromb Vasc Biol. 2006;26(6):e114-117.16627805
44. Saraieva I Benetos A Labat C Franco-Cereceda A Bäck M Toupance S Telomere length in valve tissue is shorter in individuals with aortic stenosis and in calcified valve areas Front Cell Dev Biol 2021 9 618335 10.3389/fcell.2021.618335 33777932
Saraieva I, Benetos A, Labat C, Franco-Cereceda A, Bäck M, Toupance S. Telomere length in valve tissue is shorter in individuals with aortic stenosis and in calcified valve areas. Front Cell Dev Biol. 2021;9:618335.33777932
45. Pan W Huang Q Zhou L Epigenetic age acceleration and risk of aortic valve stenosis: a bidirectional Mendelian randomization study Clin Epigenetics 2024 16 1 41 10.1186/s13148-024-01647-5 38475866
Pan W, Huang Q, Zhou L, et al. Epigenetic age acceleration and risk of aortic valve stenosis: a bidirectional Mendelian randomization study. Clin Epigenetics. 2024;16(1):41.38475866
46. Li J Wang W Yang Z Causal association of obesity with epigenetic aging and telomere length: a bidirectional Mendelian randomization study Lipids Health Dis 2024 23 1 78 10.1186/s12944-024-02042-y 38475782
Li J, Wang W, Yang Z, et al. Causal association of obesity with epigenetic aging and telomere length: a bidirectional Mendelian randomization study. Lipids Health Dis. 2024;23(1):78.38475782
47. Linterman MA Age-dependent changes in T follicular helper cells shape the humoral immune response to vaccination Semin Immunol 2023 69 101801 10.1016/j.smim.2023.101801 37379670
Linterman MA. Age-dependent changes in T follicular helper cells shape the humoral immune response to vaccination. Semin Immunol. 2023;69:101801.37379670
48. Bindra A Braunstein GD Thyroiditis Am Fam Physician 2006 73 10 1769 1776 16734054
Bindra A, Braunstein GD. Thyroiditis. Am Fam Physician. 2006;73(10):1769–76.16734054
49. Hotamisligil GS Inflammation, metaflammation and immunometabolic disorders Nature 2017 542 7640 177 185 10.1038/nature21363 28179656
Hotamisligil GS. Inflammation, metaflammation and immunometabolic disorders. Nature. 2017;542(7640):177–85.28179656
50. Nazarzadeh M Pinho-Gomes AC Bidel Z Plasma lipids and risk of aortic valve stenosis: a Mendelian randomization study Eur Heart J 2020 41 40 3913 3920 10.1093/eurheartj/ehaa070 32076698
Nazarzadeh M, Pinho-Gomes AC, Bidel Z, et al. Plasma lipids and risk of aortic valve stenosis: a Mendelian randomization study. Eur Heart J. 2020;41(40):3913–20.32076698
51. Zheng J Baird D Borges MC Recent developments in Mendelian randomization studies Curr Epidemiol Rep 2017 4 4 330 345 10.1007/s40471-017-0128-6 29226067
Zheng J, Baird D, Borges MC, et al. Recent developments in Mendelian randomization studies. Curr Epidemiol Rep. 2017;4(4):330–45.29226067
52. Ference BA Majeed F Penumetcha R Flack JM Brook RD Effect of naturally random allocation to lower low-density lipoprotein cholesterol on the risk of coronary heart disease mediated by polymorphisms in NPC1L1, HMGCR, or both: a 2 × 2 factorial Mendelian randomization study J Am Coll Cardiol 2015 65 15 1552 1561 10.1016/j.jacc.2015.02.020 25770315
Ference BA, Majeed F, Penumetcha R, Flack JM, Brook RD. Effect of naturally random allocation to lower low-density lipoprotein cholesterol on the risk of coronary heart disease mediated by polymorphisms in NPC1L1, HMGCR, or both: a 2 × 2 factorial Mendelian randomization study. J Am Coll Cardiol. 2015;65(15):1552–61.25770315
53. Zou XL Wang S Wang LY Childhood obesity and risk of stroke: a Mendelian randomisation analysis Front Genet 2021 12 727475 10.3389/fgene.2021.727475 34868204
Zou XL, Wang S, Wang LY, et al. Childhood obesity and risk of stroke: a Mendelian randomisation analysis. Front Genet. 2021;12:727475.34868204
54. Cao Z Wu Y Li Q Li Y Wu J A causal relationship between childhood obesity and risk of osteoarthritis: results from a two-sample Mendelian randomization analysis Ann Med 2022 54 1 1636 1645 10.1080/07853890.2022.2085883 35703935
Cao Z, Wu Y, Li Q, Li Y, Wu J. A causal relationship between childhood obesity and risk of osteoarthritis: results from a two-sample Mendelian randomization analysis. Ann Med. 2022;54(1):1636–45.35703935
