- 1Department of Nephrology and Endocrinology, Dongzhimen Hospital, Beijing University of Chinese Medicine, Beijing, China
- 2Department of Respiratory Disease, Dongzhimen Hospital, Beijing University of Chinese Medicine, Beijing, China
Background: Carotid Intima-Media Thickness (CIMT) is a key marker for atherosclerosis, with its modulation being crucial for cardiovascular disease (CVD) risk assessment. While thyroid function’s impact on cardiovascular health is recognized, the causal relationship and underlying mechanisms influencing CIMT remain to be elucidated.
Methods: In this study, Mendelian Randomization (MR) was employed to assess the causal relationship between thyroid function and CIMT. Thyroid hormone data were sourced from the Thyroidomics Consortium, while lipid traits and CIMT measurements were obtained from the UK Biobank. The primary analysis method was a two-sample MR using multiplicative random effects inverse variance weighting (IVW-MRE). Additionally, the study explored the influence of thyroid hormones on lipid profiles and assessed their potential mediating role in the thyroid function-CIMT relationship through multivariate MR analysis.
Results: The study revealed that lower levels of Free Thyroxine (FT4) within the normal range are significantly associated with increased CIMT. This association was not observed with free triiodothyronine (FT3), thyroid-stimulating hormone (TSH), or TPOAb. Additionally, mediation analysis suggested that apolipoprotein A-I and B are involved in the relationship between thyroid function and CIMT. The findings indicate a potential U-shaped curve relationship between FT4 levels and CIMT, with thyroid hormone supplementation in hypothyroid patients showing benefits in reducing CIMT.
Conclusion: This research establishes a causal link between thyroid function and CIMT using MR methods, underscoring the importance of monitoring thyroid function for early cardiovascular risk assessment. The results advocate for the consideration of thyroid hormone supplementation in hypothyroid patients as a strategy to mitigate the risk of carotid atherosclerosis. These insights pave the way for more targeted approaches in managing patients with thyroid dysfunction to prevent cardiovascular complications.
1 Introduction
Atherosclerosis, the primary cause of global mortality, constitutes a chronic inflammatory disease affecting large and medium-sized arteries, leading to conditions such as ischemic heart disease, strokes, and peripheral vascular disease collectively termed cardiovascular disease (CVD) (1). The prevalence of total CVD cases nearly doubled from 271 million in 1990 to 523 million in 2019, resulting in 18.6 million CVD-related deaths in 2019, as reported by the latest Global Burden of Disease (GBD) study (2). Carotid intima-media thickness (CIMT) is a widely used surrogate marker for atherosclerosis, providing a simple, reproducible, sensitive, and noninvasive measure, and a robust predictor of future cardiovascular events (3). Identifying modifiable risk factors for CIMT is paramount.
Thyroid function is recognized for its multifaceted impact on the cardiovascular system (4). Elevated circulating thyroid hormone levels have been linked to hypertension and hypercoagulation, while decreased thyroid hormone levels may lead to hyperlipidemia and inflammation. Observational studies have consistently associated variations in thyroid function within the normal range with an increased risk of atherosclerosis (5), atrial fibrillation (6), stroke (7), heart failure, and mortality (8). The underlying mechanisms involve factors such as endothelial dysfunction (9), disorders of hemostasis (10), hemodynamic changes (11), and direct effects of thyroid hormones on the myocardium (12). Despite these associations, previous research on the link between thyroid dysfunction and CIMT has yielded inconsistent conclusions (13–16). There is a pressing need to explore the relationship between thyroid function and CIMT, enabling the early assessment of atherosclerosis and CVD.
Mendelian randomization (MR) is an innovative approach employing genetic variations, typically single-nucleotide polymorphisms (SNPs), associated with modifiable exposures (or risk factors) to assess their potential causal relationship with clinically relevant outcomes. MR aims to mitigate biases arising from confounding and reverse causation, especially in situations where randomized controlled trials are impractical and observational studies may introduce biases (17–19). Over the past decades, MR has emerged as a time-efficient and cost-effective tool for prioritizing potential factors influencing human biology and disease. In this study, we investigate the existence of a causal relationship between thyroid function and CIMT. Additionally, we employ multivariate MR to analyze potential mediating mechanisms.
2 Methods
2.1 Study design
The study employed a two-sample univariate MR design, as depicted in Figure 1, to elucidate the causal relationships between a comprehensive array of thyroid function indicators and CIMT. These indicators included thyrotropin (thyroid-stimulating hormone [TSH]), thyroxine (free tetraiodothyronine [FT4]), free triiodothyronine (FT3), total triiodothyronine (TT3) concentrations, FT3/FT4, TT3/FT4, thyroid peroxidase antibodies (TPOAB), hyperthyroidism, hypothyroidism. The research framework aimed to provide a detailed examination of the impact of these thyroid function indicators on CIMT, which included assessing the effects of levothyroxine treatment in hypothyroid patients. Recognizing the established association between lipids, thyroid function, and CIMT, the study also incorporated a multivariate MR analysis to evaluate the influence of serum lipid profiles, such as apolipoprotein B (apoB), apolipoprotein A-I (apoA-I), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG), thereby exploring their potential role as mediators in this relationship.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1345267/fendo-15-1345267-HTML/image_m/fendo-15-1345267-g001.jpg)
Figure 1 Overview of Study Design and Methodology. This figure illustrates the two-sample univariate Mendelian Randomization (MR) design utilized in our study to investigate the causal relationships between various thyroid function indicators and Carotid Intima-Media Thickness (CIMT). The analysis encompassed a wide range of thyroid function indicators, including thyrotropin (TSH), free tetraiodothyronine (FT4), free triiodothyronine (FT3), total triiodothyronine (TT3), FT3/FT4 ratio, TT3/FT4 ratio, and thyroid peroxidase antibodies (TPOAB), as well as conditions like hyperthyroidism and hypothyroidism. The study also evaluated the effects of levothyroxine treatment in hypothyroid patients and incorporated a multivariate MR analysis to assess the role of serum lipid profiles (apoB, apoA-I, HDL-C, LDL-C, and triglycerides) as potential mediators in the thyroid function-CIMT relationship. This figure was created by Biorender (https://www.biorender.com/).
Adhering to the three core assumptions of MR studies, the research ensured that the selected instrumental variables (SNPs) were associated with the corresponding phenotype, not linked to confounding factors affecting exposure and outcomes, and influenced outcomes exclusively through their impact on exposure. A meticulous screening process was implemented to select SNPs, with a particular focus on eliminating those with potential pleiotropic effects. Multiple sensitivity analyses were conducted to affirm the robustness of the study’s conclusions. The utilization of published abstract data from human participant studies negated the requirement for additional ethical approval and patient consent, streamlining the research process.
2.2 Source of data
The study integrated genome-wide association study (GWAS) datasets for TSH, FT4, T3, and the T3/FT4 ratio, obtained from 46 independent cohorts within the Thyroidomics Consortium. Participants under the age of 18, of non-European ancestry, those on thyroid medication, or with a history of thyroid surgery were excluded from all analyses. We collected information on the gender distribution, average age, and thyroid hormone parameter measurements across all cohorts (refer to Figure 1; Supplementary Table 1) (20, 21). SNPs representing TPOAB were derived from a meta-analysis of TPOAb GWAS involving 18,297 individual (22). Data on patients with hyperthyroidism (3,557 cases/456,942 controls) and hypothyroidism (30,155 cases/379,986 controls) were obtained from a meta-analysis of the UK Biobank and FinnGen (23). Data representing levothyroxine sodium users were sourced from the UK Biobank (18,947 cases/443,986 controls, id: ukb-b-17918). CIMT measurements were derived from two-dimensional carotid scans of 45,185 UK Biobank participants, encompassing a diverse ethnic background (24). Lipid-related traits, including LDL cholesterol, HDL cholesterol, triglycerides, apolipoprotein B, and apolipoprotein A-I, were sourced from a GWAS of circulating lipoprotein lipids by the UK Biobank (25).
2.3 Selection of genetic instruments
The process of selecting genetic instruments encompassed several steps: (1) screening for SNPs associated with exposure factors at genome-wide significance levels (P < 5 × 10−8) for TSH, FT4, hyperthyroidism, hypothyroidism, and thyroid hormone users, and a threshold of P < 5 × 10−6 for other thyroid function indicators to increase the availability of instrumental variables. For TPOAb, significant loci reported in the original literature were used; (2) Exclusion of SNPs exhibiting linkage disequilibrium (r2 = 0.001, kb = 10,000, 1000G EUR population); (3) alignment of the respective exposure and outcome datasets based on effect allele frequencies, excluding SNPs for intermediate alleles (EAF > 0.42); (4) assessment of instrument strength to mitigate bias from weak instrumental variables, evaluated using F-statistics. F-statistics for all extracted SNPs were above 10, signifying the absence of weak instrumental bias. SNPs with potential pleiotropic effects were identified and excluded using phenoscanner (http://www.phenoscanner.medschl.cam.ac.uk/) with a threshold of P < 5 × 10-6 (refer to Supplementary Material for details).
2.4 Mendelian randomization analysis
Univariate MR analyses were conducted to examine the causal effect of each exposure factor on outcomes. The primary analysis utilized multiplicative random effects inverse variance weighting (IVW-MRE). Additional MR methods, including MR-Egger regression and weighted median, were employed to corroborate IVW estimates and ensure robustness in various scenarios. MR-Egger regression, useful for detecting directional pleiotropy, and the weighted median method, offering consistent estimates in the presence of some invalid instrumental variables (26), were utilized. Multivariate MR analyses, incorporating lipid profiles, were performed to assess whether lipid factors mediated the effect of thyroid function on CIMT. To account for the linear relationship between LDL and apoB, HDL and apoA-I, the relationship between lipid profile, FT4, and TSH was corrected using the MVMR-Lasso method (27).
2.5 Sample independence
Sample independence is crucial to avoid weak instrument bias in MR analyses (28). We ensured no overlap between subjects in samples estimating genetic associations between exposure and results. Although some sample overlap existed between thyroid disease and carotid intima thickness, the strong statistical strength of subsequent MR studies (F > 10, see Supplementary Information) mitigated concerns of weak instrument bias.
2.6 Sensitivity analysis
MR analyses are potentially vulnerable to the influences of pleiotropy and heterogeneity, which can significantly skew results and lead to erroneous interpretations. In this study, we employed MR-Egger regression to detect and adjust for pleiotropy, a method that identifies and corrects for bias caused by gene variants affecting the outcome through pathways other than the exposure of interest (29). This approach is crucial for ensuring the validity of MR findings, as it addresses the concern that genetic variants may have multiple effects that confound the estimated exposure-outcome relationship. Additionally, Cochrane’s Q-test (30) was utilized to estimate and quantify heterogeneity among the genetic instruments. Assessing heterogeneity is essential to evaluate the consistency of the instrumental variable effects and to ensure the robustness and reliability of our findings. By addressing these potential sources of bias, our study enhances the credibility of the MR approach in elucidating causal relationships.
2.7 Statistical methods
In the primary analysis, Bonferroni correction was applied to account for multiple testing in two distinct parts of the study. For the analysis examining the relationship between thyroid function and CIMT, a Bonferroni-corrected threshold of P < 0.05/9 was employed. Similarly, for the analysis exploring the association between thyroid function and lipid profiles, a corrected threshold of P < 0.05/35 was used. Results achieving these Bonferroni-corrected thresholds were deemed significant, while correlations with P values between the corrected threshold and 0.05 were considered suggestive. MR results were expressed as effects (β) with corresponding 95% confidence intervals (CI), indicating the impact of each standard deviation increase in exposure factors on outcomes. Statistical analyses employed the TwoSampleMR software package (version 0.5.7) in R (version 4.3.1). Reporting adhered to STROBE-MR guidelines (18). For detailed information on instrumental variable screening and MR methods, refer to the Supplementary Materials.
3 Results
3.1 Instrument variable selection results
Through rigorous instrumental variable screening procedures, 10 SNPs were identified as instrumental variables for FT4, 33 SNPs for TSH, 2 SNPs for TPOAb concentration, 7 SNPs for hyperthyroidism, and 50 SNPs for hypothyroidism. Additionally, 27 SNPs were identified for FT3, 11 SNPs for TT3, 37 SNPs for the FT3/FT4 ratio, and 24 SNPs for the TT3/FT4 ratio (Supplementary Table 2).
3.2 Thyroid function and CIMT
The IVW-MRE analysis demonstrated a significant association between elevated FT4 levels and reduced CIMT (beta = -0.069, 95% CI: -0.120, -0.018, P = 0.008). Conversely, FT3(beta = 0.049, 95% CI: -0.014 to 0.112, P = 0.130), TSH (beta = 0.025, 95% CI: -0.015 to 0.065, P = 0.221), and TPOAb concentrations (beta = -0.048, 95% CI: -0.505 to 0.409, P = 0.837) exhibited no significant effect on CIMT. Other thyroid function indicators did not impact CIMT (Figure 2). Hyperthyroidism was correlated with an increase in CIMT (beta = 0.053, 95% CI: 0.028 to 0.077, P = 2.8 × 10−5), whereas hypothyroidism did not demonstrate a significant effect on CIMT (beta = 0.010, 95% CI: -0.011 to 0.032, P = 0.347). Utilizing multivariate MR methods that accounted for hypothyroidism and thyroxin treatment, it was observed that hypothyroidism contributed to an increase in CIMT (beta = 0.106, 95% CI: 0.014 to 0.197, P = 0.023), while thyroxin treatment was associated with a decrease in CIMT (beta = -2.189, 95% CI: -4.287 to -0.091, P = 0.041) (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1345267/fendo-15-1345267-HTML/image_m/fendo-15-1345267-g002.jpg)
Figure 2 Causal Relationship Between Thyroid Function and Carotid Intima-Media Thickness. This figure presents the results of the Inverse Variance Weighted Multiplicative Random Effects (IVW-MRE) analysis, demonstrating a significant association between elevated FT4 levels and reduced CIMT. In contrast, FT3, TSH, and TPOAb concentrations showed no significant effect on CIMT. The figure also highlights the observed correlation between hyperthyroidism and increased CIMT, while hypothyroidism did not show a significant effect.
3.3 Thyroid function and lipid traits
An elevation in FT4 was observed to increase levels of apoA-I and HDL (apoA-I: beta = 0.037, 95% CI: 0.003, 0.072, P = 0.035; HDL: beta = 0.032, 95% CI: 0.009, 0.055, P = 0.006). Additionally, an increase in TSH was associated with a rise in apoB levels (beta = 0.022, 95% CI: 0.006, 0.038, P = 0.007). However, these findings did not meet the Bonferroni-corrected threshold. No impact of other thyroid function traits on lipid profiles was observed (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1345267/fendo-15-1345267-HTML/image_m/fendo-15-1345267-g003.jpg)
Figure 3 Causal Influence of Thyroid Function on Lipid Profiles. This figure demonstrates the effects of thyroid function on lipid profiles. An increase in FT4 is associated with higher levels of apoA-I and HDL, while a rise in TSH correlates with increased apoB levels. Notably, in this figure, asterisks indicate the level of statistical significance: a single asterisk (*) denotes a P-value less than 0.05, and double asterisks (**) denote a P-value less than 0.01.
3.4 Results of multivariate MR
After adjusting for lipid traits, FT4, and TSH using multivariate MR, FT4 was no longer associated with CIMT. Only apoA-I and apoB were causally associated with CIMT (apoA-I: beta = -0.048, 95% CI: -0.094, -0.002, P = 0.039, apoB: beta = 0.133, 95% CI: 0.086, 0.180, P = 2.62× 10−8) (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1345267/fendo-15-1345267-HTML/image_m/fendo-15-1345267-g004.jpg)
Figure 4 Multivariate MR Analysis Incorporating Lipid Factors. Image depicts the results of the multivariate MR analysis after adjusting for lipid traits, FT4, and TSH. The analysis revealed that only apoA-I and apoB were causally associated with CIMT, indicating their significant mediating role in the relationship between thyroid function and CIMT.
3.5 Sensitivity analysis
Cochran’s Q test revealed heterogeneity in a subset of the results, particularly in the associations between TSH, FT3/FT4, and TT3/FT4 with CIMT. To mitigate the impact of this heterogeneity on the study outcomes, the IVW-MRE was employed. Additionally, the MR-Egger test, utilized for detecting pleiotropy, demonstrated intercepts near zero, suggesting an absence of pleiotropy in these associations (refer to Supplementary Material for details).
4 Discussion
This study represents the first application of MR methods to decipher the causal relationship between thyroid function and CIMT. Our two-sample MR analysis delved into the impact of thyroid disorders on CIMT. We discovered that lower FT4 levels within the normal range correlate with an elevation in CIMT. Conversely, FT3, TSH, and TPOAb exhibited no significant influence on CIMT. The administration of thyroid hormone supplements in hypothyroid patients emerged as a beneficial strategy to diminish CIMT, in contrast to the CIMT augmentation observed in hyperthyroidism. Further, our exploration into the interplay between thyroid function and lipid traits revealed a positive association of FT4 with apoA-I and HDL-C, and of TSH with apoB. Intriguingly, only apoA-I and apoB demonstrated a causal relationship with CIMT. This finding suggests a potential U-shaped curve relationship between FT4 levels and CIMT, mediated by apoA-I, underscoring the potential benefits of thyroid hormone supplementation in hypothyroid patients. Our research sheds light on the need to consider the risk of carotid atherosclerosis in the early routine evaluation of individuals with thyroid dysfunction.
Thyroid dysfunction, whether overt or subclinical, has been associated with increased CIMT (31); however, the precise mechanisms underlying this association remain elusive. It is plausible that thyroid-related biomarkers such as TSH, FT3, FT4, or TPOAb may contribute to endothelial dysfunction (32, 33). Moreover, thyroid dysfunction is often accompanied by metabolic disturbances, such as hyperlipidemia, obesity, and insulin resistance (34). These factors are established as clear risk factors for carotid atherosclerosis, leading clinical studies to primarily focus on patients’ metabolic abnormalities rather than changes in thyroid hormone levels.
Our findings highlight FT4 as a pivotal independent risk factor for increased CIMT. In cases of hypothyroidism, the administration of levothyroxine, a common thyroid hormone replacement therapy, has been shown to effectively reduce CIMT. We also observed that decreased FT4 within the normal range is associated with increased CIMT, aligning with previous clinical research (35, 36). Alessandro P. Delitala et al. found that higher FT4 levels increase arterial stiffness in the common carotid artery, possibly related to its effect on heart rate (37). In our study, we also noted increased CIMT in hyperthyroid patients. In addition, we discovered that the relationship between FT4 and CIMT is mediated by apoA-I, potentially elucidating the mechanism through which FT4 influences CIMT. However, Ayse S Cikim et al. found no significant difference in CIMT between subclinical hyperthyroid patients and the general population (38). Our study proposes a U-shaped relationship between FT4 levels and CIMT, underscoring the importance of routinely assessing thyroid function in patients to monitor the risk of carotid atherosclerosis. While our findings suggest a link between lower FT4 levels within the normal range and increased CIMT, current guidelines for subclinical hypothyroidism caution against the use of thyroxine therapy in individuals with TSH levels below 10 mIU/ml (39).
Previous studies have showed that FT3 is negatively correlated with coronary artery disease (40, 41), and higher FT3 is cross-sectionally associated with higher total HDL particle concentration and with lower HDL particle size in euthyroid individuals (42). Furthermore, T3 has been found to influence atherosclerosis through thyroid hormone receptors, which regulate the atherosclerotic process independently of the lipid profile (43). These findings underscore the importance of thyroid function in the context of atherosclerotic cardiovascular disease. Although our study found that FT3 has no significant effect on CIMT, a finding that may have resulted from data source limitations, the significance of FT3 for carotid atherosclerosis risk is still something that warrants careful consideration.
Serum TSH levels, sensitive to subtle changes in thyroid function, may promote atherosclerosis by stimulating the proliferation of vascular smooth muscle cells (39). Clinical studies have reported an association between TSH levels and CIMT in postmenopausal women (33), but no such association was found before menopause. Our study found no correlation between TSH levels and CIMT, however, an elevation in TSH was found to be positively associated with increased apoB levels. Therefore, it is important to acknowledge the potential impact of TSH on CIMT within specific subgroups.
TPOAb was associated with autoimmune diseases, and TPOAb titers correlate positively with high-sensitivity C-reactive protein levels (44). Prolonged inflammation may damage the process of endothelial repair, ultimately leading to atherosclerotic lesions. Shimizu, Y. et al. conducted a prospective study on 1069 Japanese subjects with thyroid hormones within the normal range, finding that an increase in TPOAb titers leads to progressive CIMT increase (45). However, TPOAb elevation often accompanies the onset of autoimmune diseases, making it challenging to analyze its direct effect on CIMT clinically. Our study did not find a correlation between TPOAb and CIMT, but caution is needed in interpreting this result, possibly due to the limited statistical power resulting from the inclusion of a small number of SNPs in the instrumental variable.
Apart from the influence of age and gender, CIMT is determined by traditional risk factors such as blood pressure and lipid abnormalities. There is an interaction between thyroid function and lipid metabolism, with MR studies revealing a causal relationship between changes in normal-range thyroid function and the diagnosis of metabolic syndrome and lipid status (46). The apoB/apoA-I ratio in hypothyroid patients is significantly higher than in those with normal thyroid function (47, 48). Clinical studies have shown that thyroid hormone supplementation can reduce lipids in patients with subclinical hypothyroidism (49), consequently lowering CIMT (35, 50).
To clarify the impact of the interaction between thyroid function and lipids on CIMT, we conducted a multivariable MR analysis. The results indicate that apoA-I, apo B, mediate the effect of FT4 on CIMT, emphasizing the crucial role of lipids in the process of CIMT thickening. Lipid-lowering drugs, as a frontline therapy for clinical prevention of carotid plaque formation and progression, are equally applicable to hypothyroid patients with increased CIMT. Our study also suggests that combined supplementation of thyroid hormones and lipid-lowering treatment may achieve a better reduction in CIMT for patients with hypothyroidism.
The primary strength of this study lies in its utilization of a MR design, where SNP assignment is randomized and unaffected by reverse causality, thereby minimizing potential bias from confounding factors. Moreover, the extensive sample size derived from GWAS summary statistics surpasses that of a typical RCT, enhancing the reliability and, to some extent, the predictiveness of the study’s outcomes regarding a clinical trial. However, several limitations should be acknowledged. Firstly, while MR provides insights into the long-term effects of thyroid hormones on CIMT, the results may not fully align with the outcomes of short-term drug interventions. The dynamics of physiological responses to exogenous hormones might differ from the prolonged impact of endogenous hormone levels. Secondly, limitations arise from the availability of instrumental variables. The scarcity of suitable SNPs for TPOAb may compromise the statistical power of the analysis. This constraint emphasizes the need for future research to identify additional instrumental variables for a more comprehensive understanding. Thirdly, the absence of individual-level data hindered a detailed exploration of the potential nonlinear relationship between FT4 and CIMT. Understanding such nuances could aid in establishing optimal hormone supplementation regimens. Fourthly, the study did not delve into the underlying pathophysiological mechanisms linking FT4 deficiency to CIMT. A deeper exploration of these mechanisms could enhance our understanding of the causal pathways involved. Finally, the generalizability of our findings is constrained by the study’s exclusive focus on populations of European ancestry. The applicability of our results to other ethnic groups remains uncertain and warrants consideration in future studies.
In conclusion, our findings suggest that lower FT4 levels may contribute to increased CIMT through their impact on lipid levels. While this study provides insights into the potential association between thyroid function and cardiovascular health, further investigations are warranted to ascertain the effectiveness of thyroid hormone supplementation in mitigating cardiovascular disease risk, particularly among individuals with subclinical hypothyroidism. These future studies should encompass diverse populations to ensure the broader applicability of the findings.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
MZ: Conceptualization, Data curation, Writing – original draft. CZ: Data curation, Investigation, Software, Writing – original draft. XX: Data curation, Formal analysis, Writing – original draft. JL: Conceptualization, Project administration, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by Beijing University of Traditional Chinese Medicine Leader Youth Project (2023-JYB-JBQN-020) and Chinese Society of Traditional Chinese Medicine Joint Key Project Basic Research (2023DYPLHGG-11). The corresponding author JL obtained the fundings. The funders had no role in this work. The views expressed in the paper are those of the authors and not necessarily of the listed funders.
Acknowledgments
Thanks to the thyroidomics consortium, the UK Biobank and the FinnGen for providing summary-level data resources. Their invaluable contribution significantly enriched the quality and depth of our research, enabling a more comprehensive exploration of the intricate relationship between thyroid hormones and carotid intima thickness.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1345267/full#supplementary-material
References
1. Kobiyama K, Ley K. Atherosclerosis. Circ Res. (2018) 123:1118–20. doi: 10.1161/circresaha.118.313816
2. Roth GA, Mensah GA, Johnson CO, Addolorato G, Ammirati E, Baddour LM, et al. Global burden of cardiovascular diseases and risk factors, 1990-2019: update from the GBD 2019 study. J Am Coll Cardiol. (2020) 76:2982–3021. doi: 10.1016/j.jacc.2020.11.010
3. Nezu T, Hosomi N, Aoki S, Matsumoto M. Carotid intima-media thickness for atherosclerosis. J Atheroscler Thromb. (2016) 23:18–31. doi: 10.5551/jat.31989
4. Gu Y, Meng G, Zhang Q, Liu L, Wu H, Zhang S, et al. Association of longitudinal trends in thyroid function with incident carotid atherosclerosis in middle-aged and older euthyroid subjects: the tianjin chronic low-grade systemic inflammation and health (Tclsih) cohort study. Age Ageing. (2022) 51. doi: 10.1093/ageing/afab276
5. Bano A, Chaker L, Mattace-Raso FUS, van der Lugt A, Ikram MA, Franco OH, et al. Thyroid function and the risk of atherosclerotic cardiovascular morbidity and mortality: the rotterdam study. Circ Res. (2017) 121:1392–400. doi: 10.1161/circresaha.117.311603
6. Baumgartner C, da Costa BR, Collet TH, Feller M, Floriani C, Bauer DC, et al. Thyroid function within the normal range, subclinical hypothyroidism, and the risk of atrial fibrillation. Circulation. (2017) 136:2100–16. doi: 10.1161/circulationaha.117.028753
7. Chaker L, Baumgartner C, den Elzen WP, Collet TH, Ikram MA, Blum MR, et al. Thyroid function within the reference range and the risk of stroke: an individual participant data analysis. J Clin Endocrinol Metab. (2016) 101:4270–82. doi: 10.1210/jc.2016-2255
8. Cappola AR, Arnold AM, Wulczyn K, Carlson M, Robbins J, Psaty BM. Thyroid function in the euthyroid range and adverse outcomes in older adults. J Clin Endocrinol Metab. (2015) 100:1088–96. doi: 10.1210/jc.2014-3586
9. Dörr M, Empen K, Robinson DM, Wallaschofski H, Felix SB, Völzke H. The association of thyroid function with carotid artery plaque burden and strokes in a population-based sample from a previously iodine-deficient area. Eur J Endocrinol. (2008) 159:145–52. doi: 10.1530/eje-08-0140
10. Ellervik C, Mora S, Kuś A, Åsvold B, Marouli E, Deloukas P, et al. Effects of thyroid function on hemostasis, coagulation, and fibrinolysis: A mendelian randomization study. Thyroid. (2021) 31:1305–15. doi: 10.1089/thy.2021.0055
11. Ichiki T. Thyroid hormone and vascular remodeling. J Atheroscler Thromb. (2016) 23:266–75. doi: 10.5551/jat.32755
12. Jabbar A, Pingitore A, Pearce SH, Zaman A, Iervasi G, Razvi S. Thyroid hormones and cardiovascular disease. Nat Rev Cardiol. (2017) 14:39–55. doi: 10.1038/nrcardio.2016.174
13. Kim H, Kim TH, Kim HI, Park SY, Kim YN, Kim S, et al. Subclinical thyroid dysfunction and risk of carotid atherosclerosis. PloS One. (2017) 12:e0182090. doi: 10.1371/journal.pone.0182090
14. Peixoto de Miranda ÉJ, Bittencourt MS, Pereira AC, Goulart AC, Santos IS, Lotufo PA, et al. Subclinical hypothyroidism is associated with higher carotid intima-media thickness in cross-sectional analysis of the Brazilian longitudinal study of adult health (Elsa-brasil). Nutr Metab Cardiovasc Dis. (2016) 26:915–21. doi: 10.1016/j.numecd.2016.06.005
15. Zhou Y, Zhao L, Wang T, Hong J, Zhang J, Xu B, et al. Free triiodothyronine concentrations are inversely associated with elevated carotid intima-media thickness in middle-aged and elderly chinese population. J Atheroscler Thromb. (2016) 23:216–24. doi: 10.5551/jat.30338
16. Yazıcı D, Özben B, Toprak A, Yavuz D, Aydın H, Tarçın Ö, et al. Effects of restoration of the euthyroid state on epicardial adipose tissue and carotid intima media thickness in subclinical hypothyroid patients. Endocrine. (2015) 48:909–15. doi: 10.1007/s12020-014-0372-1
17. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the strobe-mr statement. Jama. (2021) 326:1614–21. doi: 10.1001/jama.2021.18236
18. Bell KJL, Loy C, Cust AE, Teixeira-Pinto A. Mendelian randomization in cardiovascular research: establishing causality when there are unmeasured confounders. Circ Cardiovasc Qual Outcomes. (2021) 14:e005623. doi: 10.1161/circoutcomes.119.005623
19. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: A guide, glossary, and checklist for clinicians. Bmj. (2018) 362:k601. doi: 10.1136/bmj.k601
20. Teumer A, Chaker L, Groeneweg S, Li Y, Di Munno C, Barbieri C, et al. Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation. Nat Commun. (2018) 9:4455. doi: 10.1038/s41467-018-06356-1
21. Sterenborg R, Steinbrenner I, Li Y, Bujnis MN, Naito T, Marouli E, et al. Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications. Nat Commun. (2024) 15:888. doi: 10.1038/s41467-024-44701-9
22. Medici M, Porcu E, Pistis G, Teumer A, Brown SJ, Jensen RA, et al. Identification of novel genetic loci associated with thyroid peroxidase antibodies and clinical thyroid disease. PloS Genet. (2014) 10:e1004123. doi: 10.1371/journal.pgen.1004123
23. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x
24. Yeung MW, Wang S, van de Vegte YJ, Borisov O, van Setten J, Snieder H, et al. Twenty-five novel loci for carotid intima-media thickness: A genome-wide association study in >45 000 individuals and meta-analysis of >100 000 individuals. Arterioscler Thromb Vasc Biol. (2022) 42:484–501. doi: 10.1161/atvbaha.121.317007
25. Richardson TG, Sanderson E, Palmer TM, Ala-Korpela M, Ference BA, Davey Smith G, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: A multivariable mendelian randomisation analysis. PloS Med. (2020) 17:e1003062. doi: 10.1371/journal.pmed.1003062
26. 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:693–8. doi: 10.1038/s41588-018-0099-7
27. Grant AJ, Burgess S. Pleiotropy robust methods for multivariable mendelian randomization. Stat Med. (2021) 40:5813–30. doi: 10.1002/sim.9156
28. Burgess S, Davies NM, Thompson SG. Bias due to participant overlap in two-sample mendelian randomization. Genet Epidemiol. (2016) 40:597–608. doi: 10.1002/gepi.21998
29. 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:512–25. doi: 10.1093/ije/dyv080
30. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
31. Saif A, Mousa S, Assem M, Tharwat N, Abdelhamid A. Endothelial dysfunction and the risk of atherosclerosis in overt and subclinical hypothyroidism. Endocr Connect. (2018) 7:1075–80. doi: 10.1530/ec-18-0194
32. Raterman HG, Voskuyl AE, Simsek S, Schreurs MW, van Hoogstraten IM, Peters MJ, et al. Increased progression of carotid intima media thickness in thyroid peroxidase antibodies-positive rheumatoid arthritis patients. Eur J Endocrinol. (2013) 169:751–7. doi: 10.1530/eje-13-0394
33. Sakamaki K, Tsunekawa K, Ishiyama N, Kudo M, Ando K, Akuzawa M, et al. Association between high normal-range thyrotropin concentration and carotid intima-media thickness in euthyroid premenopausal, perimenopausal and postmenopausal women. Maturitas. (2021) 144:29–36. doi: 10.1016/j.maturitas.2020.10.022
34. Yang W, Jin C, Wang H, Lai Y, Li J, Shan Z. Subclinical hypothyroidism increases insulin resistance in normoglycemic people. Front Endocrinol (Lausanne). (2023) 14:1106968. doi: 10.3389/fendo.2023.1106968
35. Cakal E, Turgut AT, Demirbas B, Ozkaya M, Cakal B, Serter R, et al. Effects of L-thyroxine replacement therapy on carotid intima-media thickness in patients with primary hypothyroidism. Exp Clin Endocrinol Diabetes. (2009) 117:294–300. doi: 10.1055/s-0028-1085998
36. Dullaart RP, de Vries R, Roozendaal C, Kobold AC, Sluiter WJ. Carotid artery intima media thickness is inversely related to serum free thyroxine in euthyroid subjects. Clin Endocrinol (Oxf). (2007) 67:668–73. doi: 10.1111/j.1365-2265.2007.02943.x
37. Delitala AP, Scuteri A, Fiorillo E, Orrù V, Lakatta EG, Schlessinger D, et al. Carotid beta stiffness association with thyroid function. J Clin Med. (2021) 10. doi: 10.3390/jcm10030420
38. Cikim AS, Oflaz H, Ozbey N, Cikim K, Umman S, Meric M, et al. Evaluation of endothelial function in subclinical hypothyroidism and subclinical hyperthyroidism. Thyroid. (2004) 14:605–9. doi: 10.1089/1050725041692891
39. Bekkering GE, Agoritsas T, Lytvyn L, Heen AF, Feller M, Moutzouri E, et al. Thyroid hormones treatment for subclinical hypothyroidism: A clinical practice guideline. Bmj. (2019) 365:l2006. doi: 10.1136/bmj.l2006
40. Li L, Cai G, Lu W, Li F, Yu L, Xiao J. Interaction between triglyceride-glucose index and thyroid hormones on coronary artery disease risk in patient with euthyroid. Front Endocrinol (Lausanne). (2023) 14:1255656. doi: 10.3389/fendo.2023.1255656
41. Wang P, Du R, Lin L, Ding L, Peng K, Xu Y, et al. Association between free triiodothyronine levels and peripheral arterial disease in euthyroid participants. BioMed Environ Sci. (2017) 30:128–33. doi: 10.3967/bes2017.016
42. Post A, Garcia E, Gruppen EG, Kremer D, Connelly MA, Bakker SJL, et al. Higher free triiodothyronine is associated with higher Hdl particle concentration and smaller Hdl particle size. J Clin Endocrinol Metab. (2022) 107:e1807–e15. doi: 10.1210/clinem/dgac044
43. Mörk LM, Rehnmark S, Davoodpour P, Norata GD, Larsson L, Witt MR, et al. The thyroid receptor modulator kb3495 reduces atherosclerosis independently of total cholesterol in the circulation in apoe deficient mice. PloS One. (2013) 8:e78534. doi: 10.1371/journal.pone.0078534
44. Liu J, Duan Y, Fu J, Wang G. Association between thyroid hormones, thyroid antibodies, and cardiometabolic factors in non-obese individuals with normal thyroid function. Front Endocrinol (Lausanne). (2018) 9:130. doi: 10.3389/fendo.2018.00130
45. Shimizu Y, Kawashiri SY, Noguchi Y, Nakamichi S, Nagata Y, Hayashida N, et al. Normal anti-thyroid peroxidase antibody (Tpo-Ab) titers and active arterial wall thickening among euthyroid individuals: A prospective study. J Clin Med. (2022) 11. doi: 10.3390/jcm11030521
46. Pleić N, Gunjača I, Babić Leko M, Zemunik T. Thyroid function and metabolic syndrome: A two-sample bidirectional mendelian randomization study. J Clin Endocrinol Metab. (2023) 108:3190–200. doi: 10.1210/clinem/dgad371
47. Liu FH, Hwang JS, Kuo CF, Ko YS, Chen ST, Lin JD. Subclinical hypothyroidism and metabolic risk factors association: A health examination-based study in northern Taiwan. BioMed J. (2018) 41:52–8. doi: 10.1016/j.bj.2018.02.002
48. Unal E, Akın A, Yıldırım R, Demir V, Yildiz İ, Haspolat YK. Association of subclinical hypothyroidism with dyslipidemia and increased carotid intima-media thickness in children. J Clin Res Pediatr Endocrinol. (2017) 9:144–9. doi: 10.4274/jcrpe.3719
49. Li X, Wang Y, Guan Q, Zhao J, Gao L. The lipid-lowering effect of levothyroxine in patients with subclinical hypothyroidism: A systematic review and meta-analysis of randomized controlled trials. Clin Endocrinol (Oxf). (2017) 87:1–9. doi: 10.1111/cen.13338
Keywords: thyroid hormones, Mendelian randomization, carotid intima-media thickness, apolipoprotein A-I, apolipoprotein B
Citation: Zhang M-z, Zhao C, Xing X-m and Lv J (2024) Deciphering thyroid function and CIMT: a Mendelian randomization study of the U-shaped influence mediated by apolipoproteins. Front. Endocrinol. 15:1345267. doi: 10.3389/fendo.2024.1345267
Received: 27 November 2023; Accepted: 01 March 2024;
Published: 22 March 2024.
Edited by:
Fausto Bogazzi, University of Pisa, ItalyReviewed by:
Iordanis Mourouzis, National and Kapodistrian University of Athens, GreeceViktor Kravchenko, National Academy of Sciences of Ukraine, Ukraine
Copyright © 2024 Zhang, Zhao, Xing and Lv. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Lv, lvjiebucm@163.com
†These authors have contributed equally to this work and share first authorship