Study on the mechanism of cadmium chloride pollution accelerating skin tissue metabolism disorder, aging and inhibiting hair regeneration

Drinking water contaminated by Cd2+ is one of the main pathways for Cd to enter the body. The skin barrier is destroyed when the skin is contaminated by environmental Cd2+, however, the detailed mechanism by which Cd2+ induces skin metabolic disorder, and senescence and affects hair regeneration is not completely understood. In this study, 18 C57BL/6 mice were randomly divided into a Control group, a Low-dose group, and a High-dose group with 6 mice in each group, and intragastrically administered with different concentrations of cadmium chloride once a day, respectively. After 1 month of intervention, the skin tissues on the back of mice were collected for non-targeted metabolomics analysis, and the related proteins were detected by immunofluorescence assay. Non-targeted metabolomics analysis result showed that compared with the Control group, there were 29 different metabolites, mainly including lysophospholipids, fatty acids, and bile acids, in the Low-dose group, and 39 differential metabolites in the High-dose group, in addition to the above compounds, there were more amino acid compounds, and most of the metabolites had a reduced response after administration. Immunofluorescence assay result showed that the higher the concentration of cadmium chloride led to the more obvious the proliferation inhibition and apoptosis promotion effects of skin cells, and the more significant damage to hair follicle stem cells. Thus, our findings demonstrate that cadmium chloride pollution can accelerate skin metabolism disorder, and aging and impair hair regeneration.


Introduction
Cadmium (Cd) is a kind of toxic heavy metal, which is distributed in the natural environment due to its extensive application in the industry (1), and it can also be discharged into the air through diesel exhaust. Cd 2+ is one of the major toxic substances detected in the environment. Cd 2+ pollution is still increasing in the world, and it has become a major ecological problem faced by developing countries such as China (2). Even low doses of Cd 2+ are toxic, with a long half-life of up to several decades (3). It can accumulate in the body with long-term dietary intake, resulting in and developing a range of chronic diseases, causing serious harm to human health. Chronic low doses of Cd 2+ (2 mg/kg) have been shown to induce a variety of pathological conditions, disrupting renal metabolism in rats even when exposed to very low doses of Cd 2+ (0.7 mg/kg) (4). Cd 2+ (2 mg/kg) exposure has been shown to cause fertility damage and testicular cell apoptosis in male C57BL/6J mice (5). Cd 2+ (7 mg/kg) exposure can cause liver lipid peroxidation and liver injury in male CD mice (6). In short, it affects many tissues and organs, including the liver, kidneys, brain, testis, and thymus, both in vivo and in vitro (7). Thomas et al. found that Cd concentrations in the urine of people with long-term Cd 2+ exposure were comparable to the levels of kidney and bone effects found in other populations (3). It is well known that Cd 2+ is a carcinogen, and Cd 2+ exposure can promote the development of a variety of malignant tumors, including leukemia and lung cancer.
The skin is the largest organ of the human body and is a complex human protective barrier, protecting the human body from external invasion and moisture loss through the physical and chemical barriers (8). Skin inflammation is attributed to the derangement of the epidermal barrier function (9). Studies have shown that disruption of the skin barrier can lead to the development of a variety of diseases, including atopic dermatitis, contact dermatitis, pruritus, and Sjögren syndrome, as well as systemic diseases that may be associated with aging (10,11). When the skin is contaminated by environmental Cd 2+ , the barrier is destroyed, accompanied by apoptosis, DNA damage, and lipid oxidation. At the same time, skin antioxidant enzymes such as glutathione peroxidase and methionine sulfoxide reductase are also destroyed (2). Ceramide also plays a crucial role in maintaining the skin barrier function (12). The cuticle is the skin's main mechanical barrier, formed by the multiple actions of lipids and composed of ceramide, cholesterol, and fatty acids. All three components play critical roles in skin integrity, particularly ceramide, which is essential for maintaining epidermal homeostasis, and skin ceramide levels  x Fold change (M/C) > or < . , and y P < . was the e ective variable. Blue means down, red means up. (I) Control vs. Low-dose di erential metabolite heat map. The darker the red, the higher the response, and the darker the blue, the lower the response. It can be seen that the response of most compounds decreased significantly after administration of low concentrations. decline with skin aging (11). The main mechanism by which Cd 2+ exerts its deleterious effects is through the generation of oxidative stress leading to skin senescence (13), thereby disrupting skin homeostasis. It has also been reported in the literature that disruption of skin barrier function may also be associated with autophagy (14). Although these pathways may provide useful information, the mechanisms of skin barrier disruption by Cd 2+ are still poorly understood.
Drinking water contaminated with Cd 2+ is one of the main pathways that Cd 2+ enters the body (15). Specifically, the Environmental Protection Agency (EPA) mandate that the concentration limit for Cd in drinking water should be as low as 5 ppb (16). However, emerging evidence suggests that the Cd 2+ concentration in the contaminated waters is much higher than this value (17). Some studies have shown that heavy metals have potentially harmful effects on the skin, but there are few works of literature about the effects of heavy metals on the skin, and the mechanism of Cd 2+ -induced skin tissue damage remains obscure. This study aimed to investigate how Cd 2+ induces metabolic disorder and senescence in mouse skin and to explore the effects of Cd 2+ on cell proliferation, apoptosis, and hair regeneration.

Experimental animals
Eighteen male C57BL/6 mice (5 weeks old, 10-16 g body weight) were provided by the animal experiments center of Zhejiang Chinese Medical University [Grade SPF, SCXK (Shanghai)]. All animal procedures was approved by the Zhejiang Chinese Medical University Animal Ethics Committee (No. IACUC-20211227-03) and followed the National Institutes of Health and Animal Research Committee guidelines for Animal Research.

Grouping and processing
Eighteen male C57BL/6 mice were divided into 3 groups (n = 6 per group): Control group, Low-dose group, and Highdose group. The low-dose group and the high-dose group were intragastrically administered with cadmium chloride (CdCl 2 ) once daily (2 and 7 mg/kg/time, respectively). And the mice in the Control group were intragastrically administered with equal volumes of water (100 µL). After 4 weeks treatment, the mice Du et al. .
/fpubh. . were sacrificed, and the hairs on the back were removed (first the long hairs on the back were cut off, then the remaining hairs were shaved off with a razor), and the skin on the back was taken for further examination.

LC-MS sample preparation
The skin tissue was cut into pieces, and methanol (1 mg/10 µL) was added at a ratio of 1:10 (W/V). Two magnetic beads were added, homogenized at 60 Hz for 5 min, then ultrasound was continued for 20 min, centrifugation at 13,000 r/min for 10 min, and 400 µL of the supernatant was dried with nitrogen. The supernatant was redissolved in 50 µL 70% methanol and centrifuged at 13,000 r/min for 15 min after 2 min of vortexing. The supernatant was transferred to the injection vial for testing. Some of the supernatants after centrifugation by equal absorption of homogenate, nitrogen blow-drying, re-dissolving, and preparing quality control (QC) samples. The samples were analyzed by UHPLC-QTOF/MS under liquid conditions.

UHPLC-QTOF/MS analysis Chromatography
Chromatographic separation was performed on an ExionLC system (AB Sciex, Foster City, CA, USA). A Waters Acquity HSS T3 column (2.1 × 150 mm, 1.7 µm) was applied at the temperature of 35 • C. The mobile phase A was water with 0.1% formic acid (v/v), and B was acetonitrile. The gradient was optimized as follows: 0-5 min from 3 to 8% B, 5-11 min from 8 to 30% B, 11-20 min from 30 to 80% B, 20-21 min from 80 to 95% B, 21-27 min at 95% B, then back to the initial ratio of 3% B and maintained with additional 6 min for re-equilibration. The injection volume of all samples was 2 µL.
. /fpubh. . x Fold change (M/C) > or < . , and y P < . was the e ective variable. Blue means down, red means up. (I) Control vs. High-dose di erential metabolite heat map. The darker the red, the higher the response, and the darker the blue, the lower the response. It can be seen that the response of most compounds decreased significantly after the administration of low concentrations.

Mass spectrometry
To provide high-resolution detection, a 5600 Q-TOF mass spectrometer (AB Sciex, Foster City, CA, USA) equipped with an electrospray ionization source (Turbo Ionspray) was applied. MS detection was implemented both in negative and positive ion mode with the mass rang at m/z 100-1,250. The parameters of the mass spectrometer were summarized as follows: gas 1 and gas 2, 45 psi; curtain gas, 35 psi. Heat block temperature, 550 • C; ion spray voltage, −4.5kV in negative mode and 5.5 kV in positive; declustering potential, 50 V; collision energy, ±35 V; and the collision energy spread (CES) was ±15 V. To monitor the reproducibility and stability of the acquisition system, QC samples were prepared by pooling small aliquots of each sample. The QC specimens were analyzed every five samples throughout the whole analysis procedure.

Data processing
The original map was extracted by SCIEX OS Analytics and the data matrix was transformed, including mass-to-charge ratio (m/z), retention time (RT), and intensity. All data were normalized with total peak area to generate an excel table for analysis of the subsequent metabolome. To reduce the signal interference caused by accidental error, the variables with RSD ≥ 40% in QC are eliminated in excel first.
The excel file was imported into SIMCA 14.1 (Umetrics, Umeå, Sweden) software for multivariate statistical analysis. The whole distribution of samples was observed by principal component analysis (PCA). In addition, the consistency of the samples within the group was analyzed by PCA-class analysis. Generally, under one principal component, when the samples fell outside the "2-std. dev." line, it is considered that the sample is abnormal data, and the sample data should be eliminated before the follow-up analysis.

Analysis of di erential metabolites and metabolic pathways
OPLS-DA permutation test was used to statistically analyze the validity of the OPLS-DA model. When Q 2 intersects the Yaxis, the model was validated, and then differential metabolites were screened. Based on this model, the different variables were screened according to the variable projection importance index (VIP value), and the variable with VIP > 1 was considered a meaningful variable that caused the difference. Furthermore, the partial correlation coefficient was used to screen the variables that had a great influence on the OPLS-DA model. Finally, .
/fpubh. . Mann-Whitney Test was performed on the selected variables, and P-value < 0.05 was the significant difference variable. S-plot and volcano plots are produced and visually reflect the contribution of each variable to the differential grouping. Potential markers were identified by HMDB (http://www. hmdb.ca/) and LIPID MAPS (https://www.lipidmaps.org/). The differential metabolite heatmap was made to directly reflect the response degree of the compounds after administration. Based on the results of screening and identification of significantly different metabolites, the compound name results of each group were introduced into Metabo Analyst 5.0 (http://www. metaboanalyst.ca/) for metabolic pathway analysis.
x Fold change (M/C) > or < . , and y P < . was the e ective variable. Blue means down, red means up. (I) Low-dose vs. High-dose di erential metabolite heat map. The darker the red, the higher the response, and the darker the blue, the lower the response. It can be seen that the response of most compounds decreased significantly after the administration of low concentrations.

Immunofluorescence assay
Skin tissue samples were subjected to a frozen section, rinsed in PBS, fixed in 4% PFA solution, removal of endogenous Peroxidase, BSA blocking, incubation with target primary antibody (1:500), 4 • C overnight, and PBS rinsing 3 times followed by incubation of secondary antibodies at room temperature, the nuclei were counterstained with DAPI. After sealing the slices, the images were observed and collected under fluorescence microscope, and then processed and compared with Image-J software.

Statistical analysis
All of the experimental results were expressed as the mean ± SD (standard deviation). All statistical analyses were performed using SPSS 21.0 software. The significance of differences between groups was determined by 2-tailed unpaired Student's t-test or one-way ANOVA with Dunnett's post-hoc test when samples were not distributed normally. A value of p < 0.05 was considered to be statistically significant.  Figures 2G,H is a volcano plot of Control vs. Low-dose, with red indicating variables upregulated after CdCl 2 gavage, blue indicating variables downregulated after CdCl 2 gavage, and gray indicating variables that did not differ. A total of 29 different metabolites of Control vs. Low-dose were screened and identified, which were primarily lysophospholipids, fatty acids, and bile acids. And most metabolites showed a reduced response after administration, the results are shown in Table 1. The heat map can be more intuitive to observe the difference between the two groups of variable intensity responses, the results are shown in Figure 2I.     Figures 3G,H is a volcano plot of Control vs. High-dose, with red indicating variables up-regulated after CdCl 2 gavage, blue indicating variables down-regulated after CdCl 2 gavage, and gray indicating variables that did not differ. A total of 39 different metabolites of Control vs. High-dose were screened and identified. In addition to lysophospholipids, fatty acids, and bile acids, there were also many amino acids, the results are shown in Table 2. The heat map can be more intuitive to observe the difference between the two groups of variable intensity responses, the results are shown in Figure 3I.  Table 3. The heat map can be more intuitive to observe the difference between the two groups of variable intensity responses, the results are shown in Figure 4I.

High-and low-dose Cd exposure triggers di erential metabolites in mouse skin tissues
Di erent metabolite profiles in di erent groups PLS-DA analysis was used to analyze all groups and observe the changing trend of the whole metabolic group of the three groups. Results as shown in Figures 5A-D, the three groups of Control, Low-dose, and High-dose could be distinguished, indicating the migration of the metabolic group after CdCl 2 administration. In negative ion mode, the distance between the two groups was smaller than that of the control group, suggesting that the effect of CdCl 2 on the metabolite group may be higher than that of the dose. After pooling all the differential metabolite information involved in the three groups, a total of 53 metabolites were identified, with 5 metabolites being differential metabolites shared by the three models, of which four were lysophosphatidylethanolamine; one was lysophosphatidylcholine, and 49 metabolites showed significant difference after oral administration of CdCl 2 (Low-dose or High-dose), of which 19 were common metabolites, contains phospholipids, fatty acids, acylcarnitines, bile acids compounds, the results are shown in Table 4. The number of differential metabolites is presented by the Venn diagram, as shown in Figure 5E. The response of each metabolite in three groups was compared by scattering a box plot. The results are shown in Figure 6.
Metabolomics pathway analysis of mouse skin tissues after Cd + exposure The metabolic pathway bubble map shows that the smaller the p-value, the larger the bubble, indicating that the metabolic pathway is more significant (the darker the color). According to the types of compounds involved in the metabolic pathways, the main differential compounds in Control vs. Low-dose were Fatty acid metabolism, bile acid metabolism, phospholipid metabolism, taurine metabolism, and purine metabolism. The differential metabolites of Control vs. High-dose involve more amino acid-related metabolism besides the above-mentioned metabolism. Low-dose vs. High-dose involves the metabolism of phospholipids and purines, as shown in Figure 7. The main pathways involved in the different metabolites are summarized in Table 5.
Cd + exposure a ects cell proliferation, apoptosis, and hair follicle stem cell content in skin tissue Immunofluorescence staining of skin tissue samples from three groups showed that the expression levels of cell proliferation markers Ki67 (green fluorescence) and Pcna (red fluorescence) decreased with increasing intervention concentrations. The expression of apoptosis marker Bax (red fluorescence) increased with the increase of intervention concentration, and the expression of anti-apoptosis marker Bcl-2 (red fluorescence) decreased with the increase of intervention concentration. The expression of smooth muscle marker α-SMA (green fluorescence) also decreased with the increase of intervention concentration. In the perifollicular region, the Low-dose group showed that the expression positions of the follicular stem cell markers α6 (red fluorescence), Ck15 (green fluorescence), and P63 (red fluorescence) were separated from those of the Control group, and the expression levels were decreased. In the High-dose group, the expression of hair follicle stem cell markers was further disrupted, as shown in Figures 8A,B. Based on metabolic pathway and immunofluorescence analysis, combined with KEGG and literature review results, a schematic diagram was drawn (Figure 9) to investigate the potential mechanism of the effect of CdCl 2 gavage on skin damage.

Discussion
Both endogenous and exogenous factors may contribute to skin barrier dysfunction, which may manifest as various symptoms such as skin dryness, ichthyosis, and atopic dermatitis. The lipid sheet is part of the stratum corneum, which is made up of ceramide, cholesterol and fatty acids. The stratum corneum plays a primary protective role in the barrier function, and internal water loss and invasion by microorganisms or exogenous pathogens can be prevented (18,19). The destruction of the cuticle structure may directly affect the barrier function of the skin. Fatty acids and amino acids are important nutrients in the body, which can metabolize to form acetyl CoA and participate in energy metabolism. The chemical systems of the tricarboxylic acid cycle and respiratory chains are important components of energy metabolism and occur mainly in the mitochondria. Disorders of mitochondrial function and energy metabolism can also lead to oxidative stress and ROS accumulation, which can trigger inflammation (20) and disrupt skin physiological functions. Cd 2+ is a common pollutant in the environment, more studies are needed to study the comprehensive toxic effects of Cd 2+ . Cd 2+ has toxic effects on the liver, kidney, brain, skin as well as other organs, and oxidative stress is crucial to its toxic effects (1,(21)(22)(23).
In this study, we found that both the Low-dose group (CdCl 2 2 mg/kg) and the High-dose group (CdCl 2 7 mg/kg) caused metabolic disorders in normal skin tissues of mice. Importantly, low dose Cd 2+ exposure impaired 29 metabolites in normal skin tissues of mice, while high-dose Cd 2+ exposure increased another 10 metabolites, suggesting that High-dose Cd 2+ exposure was more destructive to skin tissue metabolism. These changed metabolites were mainly concentrated in lysophospholipids, fatty acids, bile acids, and more amino acid compounds. And when the three groups of Control, Low-dose, and High-dose were compared simultaneously, five common differential metabolites were found, four of which were lysophosphatidylethanolamine and one was lysophosphatidylcholine, showing the importance of lysophospholipid metabolites in Cd 2+ exposure damage to skin tissues. Further analysis of metabolic pathways revealed that phospholipid metabolism, fatty acid metabolism, bile acid metabolism and amino acid metabolism were important metabolic pathways in Cd 2+ -exposed damaged skin, which was consistent with the combination of previous differential metabolite analysis. Besides, taurine metabolism and purine metabolism were also involved in this process. In the cell phenotype analysis, with the increase of Cd 2+ exposure dose, the cell proliferation in the skin tissue was slower, the apoptosis was faster, and the damage to hair follicle stem cells was more significant. This suggests that Cd 2+ exposure may cause accelerated aging of skin tissues and inhibit hair regeneration. Thus, our study strengthens the understanding of the damage to skin tissues by Cd 2+ exposure, provides an experimental basis for its study of potential damage mechanisms, and serves as a warning of the hazards of Cd 2+ exposure in the environment.
The main structural element of a cell membrane is phospholipids. Phospholipids (PLs) represent an important class of biomolecules, of which Glycerophospholipid (GPLs) may   play the most important role. Phospholipids are composed of the hydrophobic end of fatty acids and the hydrophilic end of substituted phosphoric acids. They are responsible not only for the organization of cell membranes but also for cell signaling, phospholipids are hydrolyzed by phospholipase into Lyso-PLs and fatty acids (24). Phospholipases are found in neutrophils, which not only secrete phospholipases but also produce ROS, a reaction that may also trigger the production of Lyso-PLs. UV radiation has been found to produce ROS causing oxidative stress, which alters cellular molecules leading to skin barrier dysfunction (25). To ward off harmful oxidation products, skin cells produce detoxifying enzymes and antioxidants. Long-wavelength ultraviolet (UVA-1) oxidation of lipids can activate NRF2-dependent antioxidant pathways (26). Sphingomyelin-containing lecithin can reduce the damage to the skin barrier and relieve skin aging by reducing the ROS level in hairless mice after UV irradiation. Our study found that phospholipids, especially lysophospholipidlike metabolites and metabolic pathways, responded less to the intervention of CdCl 2 , and the differential metabolites and metabolic pathways were more affected in the High-dose group. It is suggested that Cd 2+ may affect the balance of phospholipid metabolism (phospholipid deposition and phospholipid overhydrolysis) by inducing oxidative stress, causing dysfunction, disrupting the integrity of skin cell membranes, promoting the inflammatory process in the skin, leading to inhibition of cell proliferation (decreased Ki67 and Pcna expression) and accelerated apoptosis (increased Bax expression and decreased Bcl2 expression), which may be the exposure damages skin tissue by the first metabolic mechanism of action.
Amino acids are the basic building blocks of proteins, including the skin's most abundant fibrin, such as keratin, collagen, and elastin. Skin condition is closely related to collagen and elastic fibers, and self-repair of skin structure is essential to skin health. Amino acids are important nutrients required to promote wound healing and repair damaged skin, and also contribute to maintaining acid-base balance and moisture in cell layers, such as the stratum corneum, etc. (27). Skin is rich in extracellular matrix such as collagen, which is produced mainly by fibroblasts and protects the body from all kinds of external damage. Essential amino acids play an important role in collagen synthesis in the skin (28). Studies have shown that Branchedchain amino acids as protein components (BCAAs) can also play a role in cell signaling (29). The deficiency of leucine and L-Isoleucine can reduce the synthesis of type I and III collagen in the skin by inhibiting the action of mammalian Sirolimus target protein (mTOR) (30). Our research showed that amino acid-like metabolites and metabolic pathways were reduced in response to Cd 2+ intervention, and the differential metabolites and metabolic pathways were more affected in the High-dose group. And the expression of smooth muscle marker α-SMA also decreased with the increase of Cd 2+ intervention concentration.
. /fpubh. . It is suggested that Cd 2+ may accelerate skin aging by disrupting the skin amino acid metabolic balance, disrupting the skin collagen fiber layer (decreased α-SMA expression), and reducing skin elasticity and skin self-healing, which may be the second metabolic mechanism of action of Cd 2+ exposure damage to skin tissue. Mitochondria are one of the major producers of cellular ATP, and mitochondrial dysfunction has been associated with a variety of human diseases (31). Many of the mechanisms leading to skin aging are not fully understood, and mitochondrial function is likely to be part of a complex series of processes that lead to the decline and aging of tissue function. A growing body of research suggests that mitochondrial dysfunction and oxidative stress are key features of aging in various tissues, including the skin (32). Loss of the skin barrier during skin aging increases the risk of infection and affects wound healing. The role of mitochondria in protecting the skin barrier from microorganisms has been found to increase rapidly in glycolysis and ATP production in response to Staphylococcus aureus infection on the skin, in response to hypoxia-induced metabolic stress (33). Skin problems are directly related to mitochondrial dysfunction and the passage of complex ROS signals. Oxidative damage induced by mitochondrial ROS production is an important molecular basis for various pathophysiological conditions, including aging and cancer (34). It was found that the metabolism of fatty acids and amino acids is an important pathway of the tricarboxylic acid cycle (35). Our results revealed that fatty acid, bile acid, and amino acid metabolites and metabolic pathways were reduced in response to Cd 2+ intervention, and the differential metabolites and metabolic pathways were more affected in the High-dose group. Since fatty acid, bile acid, and amino acid metabolism are all important pathways of the tricarboxylic acid cycle, it suggests that Cd 2+ may have broken the tricarboxylic acid cycle, caused mitochondrial dysfunction, caused elevated ROS levels and the onset of oxidative stress, damaged the skin barrier, reduced the body's immune function, and thus accelerated skin aging. This may be another metabolic mechanism of action for Cd 2+ exposure damage to skin tissue.
The extensive application of Cd 2+ in industry results in serious environmental pollution, especially in the water environment, which is toxic to aquatic organisms and human beings (13,36). We call for a reduction in the discharge of heavy metal wastewater from industry and for more centralized treatment to reduce its toxic effects on humans and aquatic life, in particular the protection of endangered species. Hair growth requires ROS signaling from mitochondrial function, and Cd 2+ disrupts the antioxidant defense system and promotes .
/fpubh. . the production of lipid oxidation products by increasing ROS production and increasing DNA and protein damage (37). Certain concentrations of antioxidants AntiOxBEN2 and Antioxcin4 can induce endogenous ROS protective pathways in primary human skin fibroblasts (PHSF) (38). There are many kinds of skin cells, hair follicle stem cells are the key factor of hair regeneration, and are the foundation of the formation of functional skin. Our results showed that cell inhibition and apoptosis were most pronounced in the Highdose group, and hair follicle stem cells suffered the most severe damage (expression of α6, CK15, and P63 were all significantly decreased), suggesting that Cd 2+ may activate ROS signaling, disrupt mitochondrial function, lead to increased oxidative stress, accelerate skin aging and inhibit hair regeneration, and thus exacerbate skin function decreases. This may be one of the molecular mechanisms by which Cd 2+ exposure damages skin tissue.

Conclusion
In conclusion, our results show that lysophospholipids, fatty acids, bile acid metabolites, amino acid metabolites, and metabolic pathways changed significantly after CdCl 2 treatment, it may play a key role in the mechanism of Cd 2+ breaking the skin barrier. The disruption of metabolic homeostasis such as bile acids and purines may also affect the immune system and the body's thermogenesis. CdCl 2 inhibited the proliferation of mouse skin cells, promoted cell apoptosis, and accelerated the destruction of hair follicle stem cells. This study provides an experimental basis for predicting the risk of skin harm by Cd 2+ and exploring the potential mechanism of action and provides strong evidence for the prevention of skin harm by environmental safety. Further studies are needed to elucidate whether changes in the levels of lysophospholipids, fatty acids, bile acid metabolites, and amino acid metabolites following CdCl 2 intervention can predict skin aging risk and explore potential mechanisms.

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 authors.

Ethics statement
The animal study was reviewed and approved by the Zhejiang Chinese Medical University Animal Ethics Committee (No. IACUC-20211227-03).