Causal Relationship Between Blood Triglyceride Levels and Age Spots: A Mendelian Randomization Analysis

Abstract

Objective

This study examined the association between blood triglyceride (TG) levels (TLIB) and age spots (AS).

Methods

We acquired data from the Mendelian randomization (MR) Base database and evaluated the causal association between TLIB and AS.

Results

From genome-wide association studies, we selected 33 single nucleotide polymorphisms (SNPs) that were significantly associated with TLIB and AS. The inverse variance-weighted (IVW) and weighted median estimation methods showed that TLIB had a protective effect on AS (IVW: β=−0.214, P=0.019, odds ratio [OR]=0.807, 95% confidence interval [CI]=0.674–0.966; weighted median: β=−0.277, P=0.032, OR=0.758, 95% CI=0.589–0.977). However, the MR‐Egger analysis suggested no causal association (β=−0.234, P=0.085, OR=0.792, 95% CI=0.612–1.024). The greater precision of the weighted median estimation and IVW suggests that our results support a potential causal association between TLIB and AS.

Conclusion

The MR analysis proved that TLIB has a protective effect against AS and that triglycerides have potential preventive and therapeutic effects against AS. However, the specific dose-effect relationship requires further study.

Keywords:

• age spots
• triglycerides
• Mendelian randomization

Introduction

Age spots (AS) are uneven patches caused by pigmentation due to skin aging.¹ AS not only affects the skin’s beauty, it negatively impacts a patient’s psychological well-being and self-esteem.² The main mechanisms underlying AS development include ultraviolet (UV) radiation, natural aging, and endocrine changes.³ The current clinical treatment of AS mainly involves destroying the skin-forming spots and then stimulating skin regeneration. The main treatment methods include laser treatment, cryotherapy, chemical stripping, and microneedle treatment;^4–7 these invasive therapies often lead to skin dryness, dysplasia, and other discomfort, so we must continue exploring the mechanism of AS to develop superior treatment strategies.

Triglycerides are important lipid compounds composed of three molecules of glycerol and three molecules of fatty acids through ester bonds. Studies have shown that triglycerides are associated with skin health and can reduce UV-induced metalloproteinase-1 expression in cultured human epidermal keratinocytes.^8,9 UV irradiation is also an important factor in skin aging and AS formation.^10,11 However, the relationship between triglyceride levels in the blood (TLIB) and AS remains controversial; thus, more evidence is required.

Mendelian randomization (MR) has been used to evaluate causality in observational studies. This is based on the Mendelian principle of inheritance, which states that genetic variations are randomly inherited. Therefore, genetic variations associated with other diseases may also affect the variables of interest. MR evaluates the association between exposure and outcomes by comparing populations with different genetic variations; moreover, it enables better control of potential confounders than traditional observational studies, thereby improving the reliability of the conclusion.¹² Therefore, our study aimed to examine the association between TLIB and AS using a two-sample MR analysis.

Materials and Methods

Acquired Genetic Variant Data

We acquired data from the MR Base (http://app.mrbase.org), a platform for MR research that compiles and summarizes data from GWAS.^13,14

We used publicly available summary statistics datasets from GWAS meta-analyses of TLIB in Europe (ID: ebi-a-GCST002216, n=94,595, variants=2,410,057)¹⁵ as the exposure. We set the threshold of statistical significance as “P<5.00E⁻⁸; LD Rsq<0.1” to identify single nucleotide polymorphisms (SNPs). A total of 54 SNPs associated with triglycerides were identified. We used the same method to obtain data regarding AS in East Asian individuals (ID: ebi-a-GCST006096, n=11,253, variants= 559,843)¹⁶ as the outcome (Table 1).

Table 1 Study Exposure and Outcome

	ID	Trait	First Author	Sample Size	Variants, n	Year	Sex	Population Ethnicity	Build
Exposure	ebi-a-GCST002216	Triglycerides	Willer CJ	94595	2,410,057	2013	NA	European	HG19/GRCh37
Outcome	ebi-a-GCST006096	Age spots	Endo C	11253	559,843	2018	NA	East Asian	HG19/GRCh37

MR Analysis

Upon acquiring the data, we evaluated the association between the SNPs and triglyceride. We then separately examined the association between each SNP and AS. Finally, we used MR analysis to evaluate the causal association between triglyceride and AS. Three statistical methods were used to investigate the relationship between the TLIB and AS: MR-Egger regression, weighted median estimation, and inverse variance-weighted (IVW). MR-Egger regression is a statistical technique applied in MR studies to examine the causal associations between exposure and outcomes. This method is commonly used to identify and adjust for pleiotropic effects in which a genetic variant influences outcomes through various paths. The MR-Egger approach uses instrumental variables to estimate the causal effect of exposure on outcomes by considering any potential effects resulting from pleiotropy.¹⁷

Weighted median estimation is a statistical method that estimates the typical value of a population by computing the median score of a weighted sample. In this technique, each data point in the sample is assigned a weight that denotes its relative significance or impact on the overall estimation of population parameters. Weighted median estimation is particularly useful when outliers are present in the data distribution since it provides a more robust measure of the central tendency than traditional techniques such as the arithmetic mean.¹⁸

IVW is a statistical analysis technique commonly employed in meta-analyses to integrate the effect size estimates derived from numerous studies. In IVW, the weight assigned to each study was proportional to its inverse variance, which signified the accuracy of the effect estimate. The IVW approach constructs a pooled effect size by considering both sample size and variation in effect sizes across studies.¹⁷

Heterogeneity Test and Sensitivity Analysis

The IVW method was used to calculate the Q value¹⁹ to test for heterogeneity between SNPs. The Q_P-value was calculated to determine the heterogeneity. We also conducted a leave-one-out analysis to explore the potential influence of an individual SNP on the causal association.

Statistical Analysis

The MR analysis was completed using the “TwoSampleMR” package of the R Project (V4.2.3).²⁰ Calculations of the causal association between TLIB and AS are reported as odds ratios (OR) and corresponding 95% confidence intervals (CI). Values of P<0.05 were considered statistically significant.

Results

Detail of SNPs

We selected 33 independent SNPs from GWAS that were significantly associated with both triglyceride and AS (Table 2, Figure 1).

Table 2 Details of Each SNP Identified in the Study

SNP	effect_allele.exposure	beta.exposure	pval.exposure
rs10501321	C	−0.0216	1.41^E−08
rs10761762	C	−0.027	1.06^E−17
rs10790162	G	−0.2305	1.00^E−200
rs11057408	T	−0.0258	2.05^E−12
rs11613352	T	−0.028	9.40^E−14
rs11974409	G	−0.0899	1.36^E−100
rs1260326	C	−0.1148	1.00^E−200
rs12676857	C	0.0332	7.29^E−12
rs12678919	G	−0.1702	1.82^E−199
rs174535	C	0.047	1.73^E−41
rs17513135	T	0.022	1.63^E−08
rs2043085	C	−0.0327	7.81^E−20
rs2250802	A	0.023	1.21^E−10
rs247616	T	−0.0393	1.12^E−25
rs2665357	C	0.0212	8.33^E−10
rs287621	C	−0.0222	7.67^E−09
rs2954022	A	−0.078	2.23^E−113
rs2972146	T	0.0281	2.97^E−15
rs3760627	C	0.0189	5.29^E−09
rs3761445	A	0.0232	8.06^E−12
rs38855	G	−0.0187	2.11^E−08
rs439401	C	0.0659	1.42^E−66
rs4587594	A	−0.0694	3.50^E−82
rs4738684	G	−0.0205	8.82^E−09
rs634869	C	−0.0272	1.78^E−14
rs645040	T	0.0293	1.83^E−12
rs676210	A	−0.0733	3.28^E−71
rs6831256	G	0.0258	1.60^E−12
rs6882076	C	0.0286	1.51^E−15
rs6995541	G	0.0265	1.34^E−12
rs731839	A	−0.0224	2.65^E−09
rs749671	A	−0.0211	6.11^E−10
rs998584	A	0.0293	3.42^E−15

Abbreviation: SNP, single nucleotide protein.

Figure 1 Forest plot of SNPs associated with triglyceride and AS.

Abbreviations: AS, age spots; SNPs, single nucleotide proteins.

Mendelian Randomization

The MR-Egger regression, weighted median estimation, and IVW results are presented in Table 3 and Figure 2. The IVW and weighted median estimation methods showed that triglyceride had a protective effect on AS (IVW: β=−0.214, P=0.019, OR=0.807, 95% CI=0.674–0.966; weighted median estimation: β=−0.277, P=0.032, OR=0.758, 95% CI=0.589–0.977). On the other hand, the MR‐Egger analysis suggested no causal association between triglyceride and AS (β=−0.234, P=0.085, OR=0.792, 95% CI=0.612–1.024). However, the weighted median estimation and IVW feature greater estimation precision than the MR-Egger analysis.²¹ Therefore, our results support the potential causal association between TLIB and AS.

Table 3 Causal Associations Between Blood Triglyceride Levels and Age Spots

Method	SNP, n	b	se	pval		OR (95% CI)
MR-Egger	33	−0.234	0.132	0.085		0.792 (0.612–1.024)
Weighted median estimation	33	−0.277	0.129	0.032		0.758 (0.589–0.977)
Inverse variance-weighted	33	−0.214	0.092	0.019		0.807 (0.674–0.966)

Abbreviations: CI, confidence interval; MR, Mendelian randomization; OR, odds ratio.

Figure 2 Scatter plot of SNPs associated with triglyceride and AS. The plot presents the effect sizes of the SNP–triglyceride association (x-axis, SD units) and the SNP–AS association (y-axis, log (OR)) with 95% confidence intervals. The regression slopes of the lines correspond to causal estimates using the three MR methods (weighted median estimation, inverse variance-weighted, and Mendelian regression–Egger).

Abbreviations: AS, age spots; MR, Mendelian randomization; OR, odds ratio; SNPs, single nucleotide proteins.

Heterogeneity Test and Sensitivity Analysis

Cochran’s Q_P=0.513 (IVW) suggested no heterogeneity. Low heterogeneity indicates the reliability of the MR results.²² The horizontal pleiotropy test showed no pleiotropy in our MR analysis (MR-Egger regression intercept=0.002, standard error=0.008, P=0.840; values of P>0.05 indicate no pleiotropy and conformation to the MR). The leave-one-out analysis showed that no individual SNP significantly impacted the causal inference (Figure 3).

Figure 3 Leave-one-out analysis of SNPs associated with triglyceride and AS. Each black point represents the result of the IVW MR method applied to estimate the causal effect of triglyceride and AS.

Abbreviations: AS, age spots; IVW, inverse variance-weighted; MR, Mendelian randomization; SNPs, single nucleotide proteins.

Discussion

In the wake of enhanced living standards, the focus on aesthetic appeal and preservation of skin vitality has surged. The treatment of AS has emerged as a predominant aspect of anti-aging skincare, sparking the interest of dermatologists, cosmetologists, and skin beauty researchers worldwide. This has stimulated significant investment in research to elucidate the mechanism of AS formation.^6,23

Our study used MR as a powerful tool to measure variation in genes of known function to examine the causal effect of modifiable exposure on disease²⁴ and determine whether blood triglycerides potentially exert preventive and therapeutic effects on AS (Figure 2). Previous studies demonstrated that triglycerides retain water in the skin, and low triglyceride levels indicate insufficient or impaired lipid metabolism, which plays a crucial role in skin homeostasis.

Triglycerides are an integral component of sebum that are essential for forming the hydrophobic lipid film on the skin’s surface.²⁵ This film plays a protective role by preventing excessive water evaporation and maintaining appropriate skin hydration.²⁶ In recent years, advances in molecular technology and the theory that free radicals cause aging have been proposed; additional studies have shown that triglycerides can regulate the esterification reaction in the cell, thus reducing the accumulation of lipids in the cell and counteracting the effect of free radicals, which are among the important factors leading to AS formation.²⁷ In addition, triglycerides can regulate cholesterol metabolism, which in turn regulates cellular changes and makes the skin healthier.^9,28

These are all possible reasons why triglycerides can fight AS. Moreover, low triglyceride levels can reduce adipose tissue, which is known to exert a photoprotective effect against UV radiation. As UV radiation is a critical etiological factor in AS,²⁹ the promotion of adipogenesis by triglycerides plays an indirect role in protecting the skin. Furthermore, triglycerides also influence cholesterol metabolism, an important regulatory pathway that controls cellular changes and skin health.³⁰ Because of these multifaceted roles, triglycerides have been extensively used in topical skin care product formulations.³¹

TLIB was the exposure in this study,¹⁵ Few studies to date have examined the relationship between blood triglyceride levels and AS, which is also a highlight of this study; however, it is well known that high blood triglyceride levels may cause hyperlipidemia and do not benefit the skin’s health.³² Although we used a relatively large sample size, MR did not show a specific dose-effect relationship, which is also a study limitation. Therefore, maintaining steady-state triglyceride levels in a way that does not compromise health and can prevent senile plaques requires future investigation. Meanwhile, we found no data related to exposure and sex outcomes in the public database, which is also a disadvantage of its use. Therefore, we will conduct clinical observational studies in the future to determine the dose-effect relationship between AS and blood triglycerides to guide clinical treatment.

Conclusion

In this study, MR analysis was used to prove that blood triglycerides have a protective effect against AS and have potential preventive and therapeutic effects against it. However, the specific dose-effect relationship requires further study.

Abbreviations

AS, age spots; CI, confidence interval; GWAS, genome-wide association study; IVW, inverse variance-weighted; MR, Mendelian randomization; OR, odds ratio; SNP, single nucleotide polymorphism; TLIB, triglyceride levels in the blood.

Ethics Approval

This study was approved by the Medical Ethics Committee of Second Affiliated Hospital of Guangxi Medical University (approval no. 2023-KY (0380)).

Author Contributions

All authors made a significant contribution to the work reported, whether in the study conception, design, and execution; data acquisition, analysis, and interpretation; participated in drafting, revising, or critically reviewing the article; gave final approval of the version to be published; agree about the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Disclosure

The authors report no conflicts of interest in this work.

Acknowledgments

We thank Editage for English language editing.

Additional information

Funding

This study was supported by the 2020 Guangxi Medical and Health Appropriate Technology Development and Application Program, a study of the protective effects of dexmedetomidine in perioperative patients undergoing kidney transplantation (grant no. S2020014), and the Guangxi Medical and Health Key Discipline Construction Project.