ORIGINAL RESEARCH article

Front. Public Health, 19 October 2022

Sec. Environmental Health and Exposome

Volume 10 - 2022 | https://doi.org/10.3389/fpubh.2022.1035301

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

  • 1. Research Institute of Orthopedics, The Affiliated Jiangnan Hospital of Zhejiang Chinese Medical University, Hangzhou, China

  • 2. Hangzhou Xiaoshan Hospital of Traditional Chinese Medicine, Hangzhou, China

  • 3. Shaoxing Traditional Chinese Medicine Hospital Affiliated to Zhejiang Chinese Medical University, Shaoxing, China

  • 4. Institute of Orthopedics and Traumatology, The First Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou, China

  • 5. The First Clinical College of Zhejiang Chinese Medical University, Hangzhou, China

  • 6. School of Basic Medicine Sciences, Zhejiang Chinese Medical University, Hangzhou, China

Article metrics

View details

8

Citations

3,5k

Views

1,3k

Downloads

Abstract

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. Cd2+ is one of the major toxic substances detected in the environment. Cd2+ 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 Cd2+ 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 Cd2+ (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 Cd2+ (0.7 mg/kg) (4). Cd2+ (2 mg/kg) exposure has been shown to cause fertility damage and testicular cell apoptosis in male C57BL/6J mice (5). Cd2+ (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 Cd2+ exposure were comparable to the levels of kidney and bone effects found in other populations (3). It is well known that Cd2+ is a carcinogen, and Cd2+ 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 Cd2+, 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 decline with skin aging (11). The main mechanism by which Cd2+ 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 Cd2+ are still poorly understood.

Drinking water contaminated with Cd2+ is one of the main pathways that Cd2+ 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 Cd2+ 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 Cd2+-induced skin tissue damage remains obscure. This study aimed to investigate how Cd2+ induces metabolic disorder and senescence in mouse skin and to explore the effects of Cd2+ on cell proliferation, apoptosis, and hair regeneration.

Methods

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 High-dose group. The low-dose group and the high-dose group were intragastrically administered with cadmium chloride (CdCl2) 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 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.

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 differential metabolites and metabolic pathways

OPLS-DA permutation test was used to statistically analyze the validity of the OPLS-DA model. When Q2 intersects the Y-axis, 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, 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.

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.

Results

Repeatability and stability of the UHPLC-QTOF/MS method

The data in Figures 1A,B are QC Base Peak Chromatograms (BPC) in positive and negative ion modes. The base peak diagram is a continuous representation of the strongest ionic strength at each time point, which includes the ionic strength and the retention time of the ions in the chromatography. The PCA plots of all samples in positive and negative ion mode are shown in Figures 1C,G, respectively, showing that the QC samples clustered more closely, indicating good stability and reproducibility of this experiment. The consistency of the within-group samples was analyzed by PCA-Class analysis. From the results of Figures 1CJ, it can be seen that although some points negative ions D4, G6, positive ions C1, D4 fall on the 2std line, they do not exceed it, therefore, all data will be retained.

Figure 1

Figure 1

The BPC diagram of QC samples in positive and negative ion mode. (A) Negative ion mode, (B) Positive ion mode. (C–F) PCA analysis of negative ions. (C) PCA analysis of all samples (R2X 0.756, Q2 0.396), (D) PCA-class analysis of Control group, (E) PCA-class analysis of Low-dose group, (F) PCA-class analysis of High-dose group. (G–J) PCA analysis of positive ions. (G) PCA analysis of all samples (R2X 0.566, Q2 0.0992), (H) PCA-class analysis of Control group, (I) PCA-class analysis of Low-dose group, (J) PCA-class analysis of High-dose group.

Cd2+ exposure causes difference in metabolites in mouse skin tissue

Low-dose Cd2+ exposure results in 29 differential metabolites in mouse skin tissues

Figures 2AF shows the OPLS-DA of Control vs. Low-dose, showing that the two groups are distinguished and the model is valid. Figures 2G,H is a volcano plot of Control vs. Low-dose, with red indicating variables upregulated after CdCl2 gavage, blue indicating variables downregulated after CdCl2 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.

Figure 2

Figure 2

OPLS-DA: Control vs. Low-dose. (A) Negative ions (R2X 0.915, R2Y 1.000, Q2 0.992), (B) Positive ions (R2X 0.857, R2Y 0.999, Q2 0.738), (C) Negative ion replacement test, (D) Positive ion replacement test. (E) Negative ion S-plot, (F) Positive ion S-plot. (G,H) Volcano diagram. (G) Negative ions, (H) Positive ions. ① Fold change (M/C) > 2 or <0.5, and ② P <0.05 was the effective variable. Blue means down, red means up. (I) Control vs. Low-dose differential 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.

Table 1

MetabolitesFormulam/zRt minVIPp (corr)p-valueFold
N-Glycolylneuraminic acidC11H19NO10324.09361.33.00−0.910.00220.27
InosineC10H12N4O5267.07353.762.87−0.840.00430.32
ValerylcarnitineC12H23NO4246.17009.661.08−0.760.00220.44
Taurallocholic acidC26H45NO7S516.299014.011.06−0.390.04110.13
Taurocholic acidC26H45NO7S533.325515.193.86−0.470.01520.07
LysoPE(18:2/0:0)C23H44NO7P476.278319.772.04−0.630.01520.66
LysoPE(20:4/0:0)C25H44NO7P502.292819.824.52−0.790.00220.73
LysoPC(18:2/0:0)C26H50NO7P542.321720.141.19−0.770.00870.65
LysoPE(18:1/0:0)C23H46NO7P478.293920.872.230.730.00871.19
LysoPC(18:1/0:0)C26H52NO7P566.346320.981.59−0.630.00220.72
LysoPE(22:4/0:0)C27H48NO7P530.324121.132.050.650.04112.06
1-O-Hexadecyl-sn-glycero-3-phosphocholineC24H52NO6P482.360521.302.37−0.810.00220.47
Hydroxyeicosatetraenoic acidC20H32O3319.227921.532.01−0.600.01520.66
Monoethylhexyl phthalic acidC16H22O4279.159122.021.740.590.01521.50
LysoPE(18:0/0:0)C23H48NO7P480.309622.521.94−0.830.00220.54
MG(18:2/0:0/0:0)C21H38O4355.284323.041.16−0.580.04110.13
Hydroxyhexadecanoic acidC16H32O3271.227923.261.35−0.780.01520.44
FA 22:6C22H32O2327.233023.42.53−0.640.00870.67
FA 16:1C16H30O2253.217323.583.25−0.710.00430.62
FA 20:4C20H32O2303.233023.663.28−0.660.01520.64
FA 22:5C22H34O2329.248623.841.57−0.790.00220.54
FA 18:2C18H32O2279.233023.918.28−0.810.00430.51
MG(16:0/0:0/0:0)C19H38O4331.284324.041.150.810.00222.28
FA 20:3C20H34O2305.248624.321.14−0.760.00430.57
Hydroxyoctadecanoic acidC18H36O3299.259224.861.03−0.870.00220.42
FA 18:1C18H34O2281.248625.139.23−0.840.00220.54
FA 20:2C20H36O2307.264325.411.13−0.790.00430.54
MG(18:0/0:0/0:0)C21H42O4381.297525.871.980.760.02602.88
FA 20:1C20H38O2309.279927.272.41−0.860.00220.49

Basic information of differential metabolites in Control vs. Low-dose.

High-dose Cd2+ exposure results in 39 differential metabolites in mouse skin tissues

Figures 3AF shows the OPLS-DA of Control vs. High-dose, showing that the two groups are distinguished and the model is valid. Figures 3G,H is a volcano plot of Control vs. High-dose, with red indicating variables up-regulated after CdCl2 gavage, blue indicating variables down-regulated after CdCl2 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.

Figure 3

Figure 3

OPLS-DA: Control vs. High-dose. (A) Negative ions (R2X 0.652, R2Y 0.989, Q2 0.898), (B) Positive ions (R2X 0.787, R2Y 1.000, Q2 0.730), (C) Negative ion replacement test, (D) Positive ion replacement tests, (E) Negative ion S-plot, (F) Positive ion S-plot. (G,H) Volcano diagram. (G) Negative ions, (H) Positive ions. ① Fold change (M/C) > 2 or <0.5, and ② P < 0.05 was the effective variable. Blue means down, red means up. (I) Control vs. High-dose differential 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.

Table 2

MetabolitesFormulam/zRt minVIPp (corr)p-valueFold
L-LysineC6H14N2O2147.11281.111.27−0.780.01520.76
L-ArginineC6H14N4O2175.11901.241.76−0.860.00220.57
L-GlutamineC5H10N2O115.08661.251.25−0.840.00430.60
AnserineC10H16N4O3241.12951.261.27−0.640.02600.37
GlycerophosphocholineC8H20NO6P258.11011.272.68−0.700.00220.33
L-ValineC5H11NO2118.08631.422.15−0.870.00220.74
XanthosineC10H12N4O6285.08304.671.61−0.760.00220.41
ValerylcarnitineC12H23NO4246.17009.661.30−0.690.01520.50
Taurocholic acidC26H45NO7S533.325515.194.81−0.430.04110.15
9,10,13-TriHOMEC18H34O5329.233315.961.670.610.00872.17
Allocholic acidC24H40O5426.321417.661.81−0.520.04110.11
LysoPE(18:2/0:0)C23H44NO7P476.278319.771.87−0.560.04110.79
LysoPE(22:6/0:0)C27H44NO7P526.292819.816.190.740.00431.35
LysoPE(20:4/0:0)C25H44NO7P502.292819.825.280.710.01521.24
LysoPE(18:2/0:0)C23H44NO7P476.278320.031.63−0.720.01520.71
LysoPC(18:2/0:0)C26H50NO7P564.330720.132.08−0.700.04110.80
LysoPE(20:3/0:0)C25H46NO7P504.308520.431.080.620.04111.50
LysoPE(16:0/0:0)C21H44NO7P452.278320.72.21−0.690.01520.80
LysoPE(18:1/0:0)C23H46NO7P480.308520.905.720.830.00221.48
LysoPC(18:1/0:0)C26H52NO7P544.337420.981.41−0.760.01520.74
LysoPC(16:0/0:0)C24H50NO7P496.339821.043.970.910.00222.00
LysoPE(22:4/0:0)C27H48NO7P530.324121.133.590.830.00222.60
LysoPE(P-16:0/0:0)C21H44NO6P438.297921.243.34−0.670.00870.70
1-O-Hexadecyl-sn-glycero-3-phosphocholineC24H52NO6P482.360521.303.06−0.800.00430.49
Hydroxyeicosatetraenoic acidC20H32O3319.227921.531.68−0.590.04110.79
LysoPE(P-18:1(9Z)/0:0)C23H46NO6P462.299021.691.68−0.890.00430.63
LysoPC(19:1/0:0)C27H54NO7P536.371121.773.020.900.00222.22
LysoPC(P-18:0/0:0)C26H54NO6P508.376221.781.76−0.580.04110.61
LysoPC(19:1/0:0)C27H54NO7P536.371122.093.220.890.00222.61
LysoPE(18:0/0:0)C23H48NO7P482.324122.492.39−0.770.00220.59
Hydroxyhexadecanoic acidC16H32O3271.227923.261.64−0.810.00430.45
FA 22:6C22H32O2327.233023.42.70−0.760.00870.73
FA 20:4C20H32O2303.233023.663.49−0.790.00870.71
FA 22:5C22H34O2329.248623.841.81−0.890.00220.61
FA 18:2C18H32O2279.233023.918.29−0.870.00220.69
FA 20:3C20H34O2305.248624.321.17−0.800.01520.70
Hydroxyoctadecanoic acidC18H36O3299.259224.861.32−0.940.00220.40
FA 20:2C20H36O2307.264325.411.12−0.730.02600.71
13-DocosenamideC22H43NO338.341728.422.80−0.700.00870.64
L-LysineC6H14N2O2147.11281.111.27−0.780.01520.76

Basic information of differential metabolites in Control vs. High-dose.

High- and low-dose Cd exposure triggers differential metabolites in mouse skin tissues

Figures 4AF shows the OPLS-DA plot of Low-dose vs. High-dose, with OPLS-DA at negative ions showing relatively small differences between the two groups. Figures 4G,H is a volcano plot of Low-dose vs. High-dose, with red indicating variables up-regulated after CdCl2 gavage, blue indicating variables down-regulated after CdCl2 gavage, and gray indicating variables that did not differ. Thirteen different metabolites of Low-dose vs. High-dose were screened and identified. Most of them were phospholipids, and the High-dose group had a higher response. The results are shown in 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.

Figure 4

Figure 4

OPLS-DA: Low-dose vs. High-dose. (A) Negative ions (R2X 0.562, R2Y 0.860, Q2 0.536), (B) Positive ions (R2X 0.874, R2Y 0.999, Q2 0.901), (C) Negative ion replacement test, (D) Positive ion replacement tests, (E) Negative ion S-plot, (F) Positive ion S-plot. (G,H) Volcano diagram. (G) Negative ions, (H) Positive ions. ① Fold change (M/C) > 2 or <0.5, and ② P <0.05 was the effective variable. Blue means down, red means up. (I) Low-dose vs. High-dose differential 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.

Table 3

MetabolitesFormulam/zRt minVIPp (corr)p-valueFold
GlycerophosphocholineC8H20NO6P258.11011.271.76−0.680.00220.42
HypoxanthineC5H4N4O137.04584.181.18−0.610.04110.59
LysoPE(22:6/0:0)C27H44NO7P526.292819.815.650.730.00431.57
LysoPE(20:4/0:0)C25H44NO7P502.292819.826.660.910.00221.69
LysoPC(18:2/0:0)C26H50NO7P520.339820.141.010.700.01521.41
LysoPC(20:4/0:0)C28H50NO7P544.339820.152.160.580.04111.46
LysoPE(18:1/0:0)C23H46NO7P480.308520.905.270.880.00221.75
LysoPE(22:4/0:0)C27H48NO7P530.324121.131.530.620.01521.26
LysoPE(18:0/0:0)C23H48NO7P482.324122.493.11−0.740.00870.35
OleoylethanolamideC20H39NO2348.287323.401.36−0.630.04110.36
MG(16:0/0:0/0:0)C19H38O4331.284324.041.33−0.870.00220.30
StearoylethanolamideC20H41NO2328.321024.761.89−0.670.01520.43
MG(18:0/0:0/0:0)C21H42O4381.297525.881.08−0.820.00220.40

Basic information of differential metabolites in Low-dose vs. High-dose.

Different metabolite profiles in different 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 5AD, the three groups of Control, Low-dose, and High-dose could be distinguished, indicating the migration of the metabolic group after CdCl2 administration. In negative ion mode, the distance between the two groups was smaller than that of the control group, suggesting that the effect of CdCl2 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 CdCl2 (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.

Figure 5

Figure 5

PLS-DA diagram of all samples. (A) PLS-DA plot in negative ion mode (R2X 0.861, R2Y 0.996, Q2 0.938), (B) PLS-DA plot in positive ion mode (R2X 0.636, R2Y 0.801, Q2 0.559), (C) Negative ion replacement test, (D) Positive ion replacement test. (E) Venn diagram of the number of differential metabolites.

Table 4

No.NameFormulaRtIon modecalm/zm/zppmMS/MSClassControl vs. Low-doseControl vs. High-doseLow-dose vs. High-dose
1L-LysineC6H14N2O21.11M+H147.1128147.11290.8130, 84, 56Amino acid--
2L-ArginineC6H14N4O21.24M+H175.1190175.11910.8175, 130, 116, 70Amino acid--
3L-GlutamineC5H10N2O1.25M+H115.0866115.08725.2130, 84, 56Amino acid--
4AnserineC10H16N4O31.26M+H241.1295241.12981.1241, 170, 153, 126, 109Amino acid--
5GlycerophosphocholineC8H20NO6P1.27M+H258.1101258.11103.4258, 184, 124, 104, 86Glycerophospholipid-
6N-Glycolylneuraminic acidC11H19NO101.30M-H324.0936324.09493.9324, 236, 186, 116, 87Sialic acid--
7L-ValineC5H11NO21.42M+H118.0863118.0862−0.372, 55Amino acid-
8InosineC10H12N4O53.76M-H267.0735267.07485.0267, 135Purine-
9HypoxanthineC5H4N4O4.18M+H137.0458137.04590.8137, 119, 110, 91, 81, 77Purine--
10XanthosineC10H12N4O64.67M+H285.0830285.08362.2153, 136Purine-
11ValerylcarnitineC12H23NO49.66M+H246.1700246.17031.4246, 187, 85Carnitine-
12Taurallocholic acidC26H45NO7S14.01M+H516.2990516.2989−0.1533, 516, 498, 480, 462, 337Bile acid--
13Taurocholic acidC26H45NO7S15.19M+NH4533.3255533.32600.9533, 516, 498, 480, 462, 337Bile acid-
149,10,13-TriHOMEC18H34O515.96M-H329.2333329.23453.6329, 229, 211, 171, 139Fatty acid--
15Allocholic acidC24H40O517.66M+NH4426.3214426.32181.0426, 373, 355Bile acid--
16LysoPE(18:2/0:0)C23H44NO7P19.77M-H476.2783476.27850.5476, 279, 214Glycerophospholipid-
17LysoPE(22:6/0:0)C27H44NO7P19.81M+H526.2928526.29371.8526, 385, 354Glycerophospholipid-
18LysoPE(20:4/0:0)C25H44NO7P19.82M+H502.2928502.29361.7502, 361, 330, 287, 269, 203Glycerophospholipid
19LysoPE(18:2/0:0)C23H44NO7P20.03M-H476.2783476.27901.6476, 279, 196Glycerophospholipid--
20LysoPC(18:2/0:0)C26H50NO7P20.14M+Na520.3398520.34051.4542, 483, 337, 184, 146, 104Glycerophospholipid
21LysoPC(20:4/0:0)C28H50NO7P20.15M+H544.3398544.34010.7544, 184Glycerophospholipid--
22LysoPE(20:3/0:0)C25H46NO7P20.43M+H504.3085504.3084−0.1504, 363Glycerophospholipid--
23LysoPE(16:0/0:0)C21H44NO7P20.70M-H452.2783452.27891.4452, 255, 196Glycerophospholipid--
24LysoPE(18:1/0:0)C23H46NO7P20.90M+H480.3085480.30890.9482, 339, 308Glycerophospholipid
25LysoPC(18:1/0:0)C26H52NO7P20.98M+FA-H566.3463566.34741.9566, 506, 281Glycerophospholipid-
26LysoPC(16:0/0:0)C24H50NO7P21.04M+H496.3398496.34020.9496, 478, 313, 184, 104Glycerophospholipid--
27LysoPE(22:4/0:0)C27H48NO7P21.13M+H530.3241530.32460.9530, 389Glycerophospholipid
28LysoPE(P-16:0/0:0)C21H44NO6P21.24M+H438.2979438.29820.7438, 420, 364, 284, 266, 155Glycerophospholipid--
291-O-Hexadecyl-sn-glycero-3-phosphocholineC24H52NO6P21.30M+H482.3605482.3601−0.8482, 184, 104Glycerophospholipid-
30Hydroxyeicosatetraenoic acidC20H32O321.53M-H319.2279319.22934.4319, 301Fatty acid-
31LysoPE(P-18:1(9Z)/0:0)C23H46NO6P21.69M-H462.2990462.30002.2462, 265, 196, 140, 78Glycerophospholipid--
32LysoPC(19:1/0:0)C27H54NO7P21.77M+H536.3711536.37150.8536, 518, 184Glycerophospholipid--
33LysoPC(P-18:0/0:0)C26H54NO6P21.78M+H508.3762508.3761−0.1508, 367, 184, 104Glycerophospholipid--
34Monoethylhexyl phthalic acidC16H22O422.02M+H279.1591279.15982.7279, 149, 57Phthalate--
35LysoPC(19:1/0:0)C27H54NO7P22.09M+H536.3711536.37150.8536, 518, 184, 104Glycerophospholipid--
36LysoPE(18:0/0:0)C23H48NO7P22.49M+H482.3241482.32461.0482, 464, 421, 341, 310Glycerophospholipid
37MG(18:2/0:0/0:0)C21H38O423.04M+H355.2843355.28481.6355, 337, 266, 263, 245Glycerolipid--
38Hydroxyhexadecanoic acidC16H32O323.26M-H271.2279271.22904.3271, 225Fatty acid-
39OleoylethanolamideC20H39NO223.4M+Na348.2873348.28781.4326, 309, 308, 62Ethanol amine--
40FA 22:6C22H32O223.4M-H327.2330327.23434.0327, 283Fatty acid, Docosahexaenoic acid-
41FA 16:1C16H30O223.58M-H253.2173253.21875.7253, 126Fatty acid, Palmitoleic acid--
42FA 20:4C20H32O223.66M-H303.2330303.23475.9303, 259Fatty acid, Arachidonic acid-
43FA 22:5C22H34O223.84M-H329.2486329.24973.3329, 297, 281, 149Fatty acid, Docosapentaenoic acid-
44FA 18:2C18H32O223.91M-H279.2330279.23434.9279, 261Fatty acid, Linoleic acid-
45MG(16:0/0:0/0:0)C19H38O424.04M+H331.2843331.28471.3313, 257, 239, 95, 71, 57Glycerolipid-
46FA 20:3C20H34O224.32M-H305.2486305.24994.4305Fatty acid, Eicosatrienoic acid-
47StearoylethanolamideC20H41NO224.76M+H328.3210328.32192.6328, 62Ethanol amine--
48Hydroxyoctadecanoic acidC18H36O324.86M-H299.2592299.25971.6299, 253Fatty acid-
49FA 18:1C18H34O225.13M-H281.2486281.25015.3281, 263Fatty acid, Oleic acid--
50FA 20:2C20H36O225.41M-H307.2643307.26533.6307Fatty acid, Eicosadienoic acid-
51MG(18:0/0:0/0:0)C21H42O425.88M+H381.2975381.29811.6341, 281, 267, 109, 95, 57Glycerolipid-
52FA 20:1C20H38O227.27M-H309.2799309.28113.7309, 291Fatty acid, Eicosenoic acid--
5313-DocosenamideC22H43NO28.42M+H338.3417338.34231.7338, 321, 303, 177, 149, 135, 83, 69Fatty amide--
Total293913

Classification of differential metabolites.

Figure 6

Figure 6

Comparison of three groups of differential metabolites with scatter box plot.

Metabolomics pathway analysis of mouse skin tissues after Cd2+ 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.

Figure 7

Figure 7

Bubble diagram of the metabolic pathway. (A) Control vs. Low-dose; (B) Control vs. High-dose; (C) Low-dose vs. High-dose. (1) Biosynthesis of unsaturated fatty acids; (2) Aminoacyl-tRNA biosynthesis; (3) Arginine biosynthesis; (4) Linoleic acid metabolism; (5) Glycerophospholipid metabolism; (6) D-Glutamine and D-glutamate metabolism; (7) Nitrogen metabolism; (8) Valine, leucine, and isoleucine biosynthesis; (9) Taurine and hypotaurine metabolism; (10) Biotin metabolism; (11) Histidine metabolism; (12) Purine metabolism; (13) Pantothenate and CoA biosynthesis; (14) Ether lipid metabolism; (15) beta-Alanine metabolisms; (16) Lysine degradation; (17) Alanine, aspartate, and glutamate metabolism; (18) Glyoxylate and dicarboxylate metabolism; (19) Arachidonic acid metabolism; (20) Arginine and proline metabolism; (21) Pyrimidine metabolism; (22) Valine, leucine, and isoleucine degradation; (23) Primary bile acid biosynthesis; (24) Amino sugar and nucleotide sugar metabolism.

Table 5

Pathwaysp–log(p)Holm pFDRImpactDetailsComments
Control vs. Low-dose
Fatty acidsBiosynthesis of unsaturated fatty acids0.00004.82860.00120.00120.0000KEGGUnsaturated fatty acids
Linoleic acid metabolism0.04891.31081.00001.00001.0000KEGGFA 18:2
Arachidonic acid metabolism0.30560.51491.00001.00000.3329KEGG SMPFA 20:4
GlycerophospholipidsGlycerophospholipid metabolism0.30560.51491.00001.00000.0174KEGGLPC; LPE
Bile acidsTaurine and hypotaurine metabolism0.07711.11281.00001.00000.0000KEGG SMPTaurocholic acid
Primary bile acid biosynthesis0.37350.42781.00001.00000.0229KEGG SMPTaurocholic acid; Taurallocholic acid
Sialic acidAmino sugar and nucleotide sugar metabolism0.31260.50501.00001.00000.0106KEGG SMP SMPN-Glycolylneuraminic acid
PurinesPurine metabolism0.49110.30891.00001.00000.0025KEGG SMPInosine
Control vs. High-dose
Amino acidsAminoacyl-tRNA biosynthesis0.00202.69830.16630.08410.0000KEGGL-Arginine; L-Glutamine; L-Valine; L-Lysine
Arginine biosynthesis0.01131.94760.92510.31590.0761KEGGL-Arginine; L-Glutamine;
D-Glutamine and D-glutamate metabolism0.06971.15671.00000.83660.0000KEGG SMPL-Glutamine
Nitrogen metabolism0.06971.15671.00000.83660.0000KEGGL-Glutamine
Valine, leucine and isoleucine biosynthesis0.09191.03661.00000.85790.0000KEGGL-Valine
Biotin metabolism0.11360.94461.00000.95440.0000KEGG SMPL-Lysine
Histidine metabolism0.17580.75501.00001.00000.0492KEGG SMPL-Anserine
Lysine degradation0.26140.58261.00001.00000.0000KEGG SMPL-Lysine
Arginine and proline metabolism0.37040.43131.00001.00000.0579KEGG SMPL-Arginine
Valine, leucine and isoleucine degradation0.38580.41371.00001.00000.0000KEGG SMPL-Valine
Glyoxylate and dicarboxylate metabolism0.32210.49201.00001.00000.0000KEGGL-Glutamine
beta-Alanine metabolism0.22450.64881.00001.00000.0000KEGG SMPL-Anserine
Alanine, aspartate and glutamate metabolism0.28810.54051.00001.00000.1138KEGG SMPL-Glutamine
Pantothenate and CoA biosynthesis0.20530.68751.00001.00000.0000KEGG SMPL-Valine
Pyrimidine metabolism0.37810.42241.00001.00000.0000KEGG SMPL-Glutamine
Fatty acidsBiosynthesis of unsaturated fatty acids0.00073.17790.05580.05580.0000KEGGUnsaturated fatty acids
Linoleic acid metabolism0.05841.23341.00000.83661.0000KEGGFA 18:2
Arachidonic acid metabolism0.35470.45021.00001.00000.3329KEGG SMPFA 20:4
GlycerophospholipidsGlycerophospholipid metabolism0.06691.17471.00000.83660.0655KEGGLPC; LPE
Ether lipid metabolism0.21500.66761.00001.00000.0000KEGGGlycerophosphocholine
Bile acidsTaurine and hypotaurine metabolism0.09191.03661.00000.85790.0000KEGG SMPTaurocholic acid
Primary bile acid biosynthesis0.42970.36681.00001.00000.0229KEGG SMPTaurocholic acid; Allocholic acid
PurinesPurine metabolism0.18490.73321.00001.00000.0000KEGG SMPXanthosine; L-Glutamine
Low-dose vs. High-dose
GlycerophospholipidsGlycerophospholipid metabolism0.00532.27470.44620.44620.0655KEGGLPC; LPE
Ether lipid metabolism0.06471.18881.00001.00000.0000KEGGGlycerophosphocholine
PurinesPurine metabolism0.20100.69681.00001.00000.0217KEGG SMPHypoxanthine

Metabolic pathway information.

Cd2+ exposure affects 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 CdCl2 gavage on skin damage.

Figure 8

Figure 8

(A) Immunofluorescence detection of markers in different groups (X200). (B) Comparison of immunofluorescence gray values of markers in different groups. Data are represented as means ± SD (*P < 0.05 vs. control).

Figure 9

Figure 9

Schematic diagram of skin tissue damage caused by CdCl2 gavage. From left to right: CdCl2 enters the body and acts on all layers of skin tissue, inducing inflammation and impairment of skin barrier function. This action may occur through the following pathways: (1) affecting the balance of phospholipid metabolism, disrupting the integrity of cell membranes and promoting the inflammatory process, leading to inhibition of cell proliferation (reduced Ki67 and Pcna expression) and accelerated apoptosis (increased Bax expression and reduced Bcl2 expression). (2) Affects amino acid metabolic balance, disrupts collagen fibers (decreased α-SMA expression) and skin repair functions. (3) Affects the metabolic balance of fatty acids, amino acids, bile acids and purines, leading to disruption of energy metabolism, increased oxidative stress, impaired skin immune function and destruction of hair follicle structure (decreased expression of α6, CK15, and P63).

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. Cd2+ is a common pollutant in the environment, more studies are needed to study the comprehensive toxic effects of Cd2+. Cd2+ has toxic effects on the liver, kidney, brain, skin as well as other organs, and oxidative stress is crucial to its toxic effects (1, 2123).

In this study, we found that both the Low-dose group (CdCl2 2 mg/kg) and the High-dose group (CdCl2 7 mg/kg) caused metabolic disorders in normal skin tissues of mice. Importantly, low dose Cd2+ exposure impaired 29 metabolites in normal skin tissues of mice, while high-dose Cd2+ exposure increased another 10 metabolites, suggesting that High-dose Cd2+ 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 Cd2+ 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 Cd2+-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 Cd2+ 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 Cd2+ 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 Cd2+ exposure, provides an experimental basis for its study of potential damage mechanisms, and serves as a warning of the hazards of Cd2+ 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 lysophospholipid-like metabolites and metabolic pathways, responded less to the intervention of CdCl2, and the differential metabolites and metabolic pathways were more affected in the High-dose group. It is suggested that Cd2+ may affect the balance of phospholipid metabolism (phospholipid deposition and phospholipid over-hydrolysis) 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 Branched-chain 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 Cd2+ 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 Cd2+ intervention concentration. It is suggested that Cd2+ 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 Cd2+ 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 Cd2+ 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 Cd2+ 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 Cd2+ exposure damage to skin tissue.

The extensive application of Cd2+ 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 Cd2+ disrupts the antioxidant defense system and promotes 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 High-dose group, and hair follicle stem cells suffered the most severe damage (expression of α6, CK15, and P63 were all significantly decreased), suggesting that Cd2+ 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 Cd2+ 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 CdCl2 treatment, it may play a key role in the mechanism of Cd2+ 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. CdCl2 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 Cd2+ 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 CdCl2 intervention can predict skin aging risk and explore potential mechanisms.

Funding

This work was supported by National Natural Science Foundation of China (No. 81904053), Special Research Project of the Affiliated Hospital of Zhejiang Chinese Medical University (No. 2021FSYYZY43), Hangzhou Medical and Health Technology Planning Project (No. B20220021), Hangzhou Science and Technology Planning Project (No. 2020ZDSJ0042), and Hangzhou Xiaoshan District Science and Technology Planning Project (No. 2019216).

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.

Statements

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

Author contributions

WD, ZW, MW, JJ, and SY did the data collection and writing. WD, YD, and ZW were major contributors to writing the manuscript. HR and RQ contributed to the conception and design of the study. All authors read and approved the final manuscript.

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.

References

  • 1.

    MeiHYaoPWangSSLiNZhuTFChenXFet al. Chronic Low-dose cadmium exposure impairs cutaneous wound healing with defective early inflammatory responses after skin injury. Toxicol Sci. (2017) 159:32738. 10.1093/toxsci/kfx137

  • 2.

    ChavatteLJuanMMounicouSNoblesseELPaysKNizardCet al. Elemental and molecular imaging of human full thickness skin after exposure to heavy metals. Metallomics. (2020) 12:155562. 10.1039/d0mt00121j

  • 3.

    ThomasLDKHodgsonSNieuwenhuijsenMJarupL. Early kidney damage in a population exposed to cadmium and other heavy metals. Environ Health Perspect. (2009) 117:18184. 10.1289/ehp.11641

  • 4.

    XuMYWangPSunYJWuYJ. Disruption of kidney metabolism in rats after subchronic combined exposure to Low-dose cadmium and chlorpyrifos. Chem Res Toxicol. (2019) 32:12229. 10.1021/acs.chemrestox.8b00219

  • 5.

    LiMLiuCYangLLZhangLChenCHHeMDet al. G9a-mediated histone methylation regulates cadmium-induced male fertility damage in pubertal mice. Toxicol Lett. (2016) 252:1121. 10.1016/j.toxlet.2016.04.004

  • 6.

    EyblVKotyzovaDKoutenskyJ. Comparative study of natural antioxidants-curcumin, resveratrol and melatonin-in cadmium-induced oxidative damage in mice. Toxicology. (2006) 225:1506. 10.1016/j.tox.2006.05.011

  • 7.

    SonekatsuMYamadaHNishioNGuJGG. Effects on low threshold mechanoreceptors in whisker hair follicles by 5-HT, Cd2+, tetraethylammonium, 4-aminopyridine, and Ba2. Mol Pain. (2022) 18:17448069221076606. 10.1177/17448069221076606

  • 8.

    KauthMTrusovaOV. Topical ectoine application in children and adults to treat inflammatory diseases associated with an impaired skin barrier: a systematic review. Dermatol Ther. (2022) 12:295313. 10.1007/s13555-021-00676-9

  • 9.

    KuboANagaoKAmagaiM. Epidermal barrier dysfunction and cutaneous sensitization in atopic diseases. J Clin Invest. (2012) 122:44047. 10.1172/JCI57416

  • 10.

    KimBELeungDYM. Significance of skin barrier dysfunction in atopic dermatitis. Allergy Asthma Immunol Res. (2018) 10:20715. 10.4168/aair.2018.10.3.207

  • 11.

    WangZManMQLiTEliasPMMauroTM. Aging-associated alterations in epidermal function and their clinical significance. Aging. (2020) 12:555165. 10.18632/aging.102946

  • 12.

    LeoTKTanESSAminiFRehmanNNgESCTanCK. Effect of rice (oryza sativa L) ceramides supplementation on improving skin barrier functions and depigmentation: an open-label prospective study. Clin Trial. (2022) 14:2737. 10.3390/nu14132737

  • 13.

    KrutmannJLiuWLiLPanXCCrawfordMSoreGet al. Pollution and skin: from epidemiological and mechanistic studies to clinical implications. J Dermatol Sci. (2014) 76:1638. 10.1016/j.jdermsci.2014.08.008

  • 14.

    YoonSJLimCJChungHJKimJHHuhYHParkKet al. Autophagy activation by crepidiastrum denticulatum extract attenuates environmental pollutant-induced damage in dermal fibroblasts. Int J Mol Sci. (2019) 20:517. 10.3390/ijms20030517

  • 15.

    JaywantSAArifKM. A comprehensive review of microfluidic water quality monitoring sensors. Sensors. (2019) 19:4781. 10.3390/s19214781

  • 16.

    BazylewskiPMiddelkoopSVDivigalpitiyaRFanchiniG. Solid-state chemiresistors from two-dimensional MOS2 nanosheets functionalized with l-cysteine for in-line sensing of part-per-billion Cd2+ ions in drinking water. ACS Omega. (2019) 5:6439. 10.1021/acsomega.9b03246

  • 17.

    KhoshbinZMoeenfardMZahraeeHDavoodianN. A fluorescence imaging-supported aptasensor for sensitive monitoring of cadmium pollutant in diverse samples: A critical role of metal organic frameworks. Talanta. (2022) 246:123514. 10.1016/j.talanta.2022.123514

  • 18.

    KimYLimKM. Skin barrier dysfunction and filaggrin. Arch Pharm Res. (2021) 44:3648. 10.1007/s12272-021-01305-x

  • 19.

    DanzbergerJDonovanMRanklCZhuRVicicSBaltenneckCet al. Glycan distribution and density in native skin's stratum corneum. Skin Res Technol. (2018) 24:4508. 10.1111/srt.12453

  • 20.

    ChenCHouGZengCRenYChenXPengC. Metabolomic profiling reveals amino acid and carnitine alterations as metabolic signatures in psoriasis. Theranostics. (2021) 11:75467. 10.7150/thno.51154

  • 21.

    MouroVGSLadeiraLCMLoziAAde MedeirosTSSilvaMRde OliveiraELet al. Different routes of administration lead to different oxidative damage and tissue disorganization levels on the subacute cadmium toxicity in the liver. Biol Trace Elem Res. (2021) 199:462434. 10.1007/s12011-020-02570-5

  • 22.

    ZhouJZengLZhangYWangMLiYJiaYet al. Cadmium exposure induces pyroptosis in testicular tissue by increasing oxidative stress and activating the AIM2 inflammasome pathway. Sci Total Environ. (2022) 847:157500. 10.1016/j.scitotenv.2022.157500

  • 23.

    MohamedHRH. Alleviation of cadmium chloride-induced acute genotoxicity, mitochondrial DNA disruption, and ROS generation by chocolate coadministration in mice liver and kidney tissues. Biol Trace Elem Res. (2022) 200:375061. 10.1007/s12011-021-02981-y

  • 24.

    FuchsB. Mass spectrometry and inflammation–MS methods to study oxidation and enzyme-induced changes of phospholipids. Anal Bioanal Chem. (2014) 406:1291306. 10.1007/s00216-013-7534-5

  • 25.

    AhnYKimMGJoKHongKBSuhHJ. Effects of sphingomyelin-containing milk phospholipids on skin hydration in UVB-exposed hairless mice. Molecules. (2022) 27:2545. 10.3390/molecules27082545

  • 26.

    GruberFMayerHLengauerBMlitzVSandersJMKadlAet al. NF-E2-related factor 2 regulates the stress response to UVA-1-oxidized phospholipids in skin cells. FASEB J. (2010) 24:3948. 10.1096/fj.09-133520

  • 27.

    SolanoF. Metabolism and functions of amino acids in the skin. Adv Exp Med Biol. (2020) 1265:18799. 10.1007/978-3-030-45328-2_11

  • 28.

    SzokaLKarnaEHlebowicz-SaratKKaraszewskiJPalkaJA. Exogenous proline stimulates type I collagen and HIF-1I expression and the process is attenuated by glutamine in human skin fibroblasts. Mol Cell Biochem. (2017) 435:197206. 10.1007/s11010-017-3069-y

  • 29.

    KimballSRJeffersonLS. Signaling pathways and molecular mechanisms through which branched-chain amino acids mediate translational control of protein synthesis. J Nutr. (2006) 136:227S−31S. 10.1093/jn/136.1.227S

  • 30.

    YamaneTMoriokaYKitauraYIwatsukiKShimomuraYOishiY. Branched-chain amino acids regulate type I tropocollagen and type III tropocollagen syntheses via modulation of mTOR in the skin. Biosci Biotechnol Biochem. (2018) 82:61115. 10.1080/09168451.2017.1386084

  • 31.

    BulthuisEPEinerCDistelmaierFGrohLVriesSEVEWesterloEVDet al. The decylTPP mitochondria-targeting moiety lowers electron transport chain supercomplex levels in primary human skin fibroblasts. Free Radic Biol Med. (2022) 188:43446. 10.1016/j.freeradbiomed.2022.06.011

  • 32.

    StoutRBirch-MachinM. Mitochondria's role in skin ageing. Biology. (2019) 8:29. 10.3390/biology8020029

  • 33.

    WickershamMWachtelSLungTWFSoongGJacquetRRichardsonAet al. Metabolic stress drives keratinocyte defenses against staphylococcus aureus infection. Cell Rep. (2017) 18:274251. 10.1016/j.celrep.2017.02.055

  • 34.

    SreedharAAguilera-AguirreLSinghKK. Mitochondria in skin health, aging, and disease. Cell Death Dis. (2020) 11:444. 10.1038/s41419-020-2649-z

  • 35.

    LiaoWTJinQWLiuJNRuan YL LiXRShenYYet al. Mahuang decoction antagonizes acute liver failure via modulating tricarboxylic acid cycle and amino acids metabolism. Front Pharmacol. (2021) 12:599180. 10.3389/fphar.2021.599180

  • 36.

    PrinsJMFuLJGuoLWangYS. Cd2+ -induced alteration of the global proteome of human skin fibroblast cells. J Proteome Res. (2014) 13:167787. 10.1021/pr401159f

  • 37.

    YuXJYuRQGuiDZhangXYZhanFPSunXet al. Hexavalent chromium induces oxidative stress and mitochondria-mediated apoptosis in isolated skin fibroblasts of indo-pacific humpback dolphin. Aquat Toxicol. (2018) 203:17986. 10.1016/j.aquatox.2018.08.012

  • 38.

    TeixeiraJBasitFWillemsPHGMWagenaarsJAWesterloEVDAmorimRet al. Mitochondria-targeted phenolic antioxidants induce ROS-protective pathways in primary human skin fibroblasts. Free Radic Biol Med. (2021) 163:31424. 10.1016/j.freeradbiomed.2020.12.023

Summary

Keywords

skin, cadmium pollution, metabolism, aging, hair regeneration

Citation

Du W, Dong Y, Wang Z, Yao S, Wang M, Ji J, Ruan H and Quan R (2022) Study on the mechanism of cadmium chloride pollution accelerating skin tissue metabolism disorder, aging and inhibiting hair regeneration. Front. Public Health 10:1035301. doi: 10.3389/fpubh.2022.1035301

Received

02 September 2022

Accepted

04 October 2022

Published

19 October 2022

Volume

10 - 2022

Edited by

Zhiyuan Meng, Yangzhou University, China

Reviewed by

Jin Yan, Suzhou University of Science and Technology, China; Feng Zhao, Guangxi University, China

Updates

Copyright

*Correspondence: Weibin Du Hongfeng Ruan

†These authors have contributed equally to this work

This article was submitted to Environmental health and Exposome, a section of the journal Frontiers in Public Health

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics