Study population
Skin lesions were derived from 21 psoriasis patients (see Supplementary Tables 1). The diagnosis of psoriasis was confirmed and quantified by a dermatologist using the Psoriasis Area Severity Index (PASI) 20, ensuring that all psoriasis patients were treated in the same manner. None of the included study participants received any systemic anti-psoriasis treatments or topical therapy within 4 weeks before the collection of biopsy skin samples, including psoriatic lesions and non-lesions. These samples were obtained by collecting 3 mm surgical biopsies from representative lesions.
Mice
Six-week-old BALB/c female mice were used to establish the IMQ-induced psoriasis-like mouse model. The mice were bred and housed in the Department of Laboratory Animals, Central South University, and maintained under specific pathogen-free conditions. Before each experiment, the mice were randomly assigned to the control or treatment group. No sample-size estimation was performed before the experimentation.
IMQ-induced psoriasis-like mouse establishment and measurement of skin scores and epidermal thickness
Six-week-old BALB/c female mice were used. Metabolomic profiling of the skin lesions of the IMQ-induced psoriasis-like mouse model was conducted as follows. A daily dose of 62.5 mg of 5% IMQ cream (Med-shine Pharmaceutical Co., Ltd.) was applied to the shaved backs of mice for 6 consecutive days. C1P treatment of the IMQ-induced psoriasis-like mouse model was as follows. A daily dose of 31.25 mg of 5% IMQ cream was applied to each of the ears of the mice, and two ears were treated. IMQ cream was applied for 6 consecutive days. The ears were intradermally injected with 1 mM C16 ceramide-1-Phosphate (d18:1/16:0) (Avanti Cat. 860533P) (10 μl) or vehicle control (normal saline) once per day for 6 consecutive days. The mice were sacrificed on Day 7. The clinical skin scores of the mice were determined beginning on Day 1 (the first day of IMQ treatment) and every other day until Day 7 using the modified psoriasis severity index score (PASI). The degree of skin erythema, induration, and scaling was classified as follows: 0, no symptoms, 1, mild, 2, moderate, 3, severe, or 4, very severe. Histological evaluation of mouse ears samples was performed by H&E staining. For each section, the epidermal thickness was measured from the stratum basal to the stratum granulosum using Image Pro-Plus. The average value of 7 random fields of view was calculated for each mouse.
LC-MS/MS based lipidomics of psoriatic skin
The lipids were extracted from skin samples following protocols as previously reported protocols 21. In short, 10 mg skin samples were weighted with 0.5 ml of prechilled solution of dichloromethane/methanol (3:1) spiked in Lipid internal standard mix (SPLASH® LIPIDOMIX® Mass Spec Standard, Avanti, USA). The mixture was lysed using a TissueLyser by 50HZ for 10 minutes. After centrifuge the mixture at 10,000g for 10 min at 4 °C, the supernatant was moved into glass vials for LC-MS analysis. Equal amount supernatant from each sample was mixed as quality control samples (QCs). Lipids were separated by a CSH C18 column (1.7 μm 2.1*100 mm, Waters, USA) and analyzed by a QExactive mass spectrometer (Thermo Fisher Scientific, USA) for qualification and quantification. Mobile phase A consisted of 10 mM ammonium formate with 0.1% formic acid in water-ACN (40:60, v/v), and mobile phase B was 10 mM ammonium formate with 0.1% formic acid in ACN-IPA (10:90, v/v). The following gradient was used for elution: 0~2 min, 40%-43% B solution, 2~2.1 min, 43%~50% B solution, 2.1~7 min, 50%~54% B solution, 7-7.1 min, 54%~70% B solution, and 7.1-13 min, 70%~99% B solution with a flow rate of 0.35 mL/min. The mass spectrometric settings for positive/negative ionization modes were as follows: spray voltage, 3.8/–3.2 kV, aux gas heater temperature, 350 °C, and capillary temperature, 320 °C. The full scan range was 200–2000 m/z with a resolution of 70,000, and the AGC target for MS acquisitions was set to 3e6 with a maximum ion injection time of 100 ms. The top three precursors were selected for subsequent MSMS fragmentation with a maximum ion injection time of 50 ms and resolution of 17,500, and the AGC was 1e5. The stepped normalized collision energy was set to 15, 30, and 45 eV. QC samples were injected after every 10 sample injections followed previous description 22 for system quality control.
LipidSearch 4.1 SP2 software (Thermo Fisher, USA) was used for lipid identification and quantitation. The processing parameters for product ion search including precursor and product tolerance as 10 ppm, adduct ions of positive mode as +H and +NH4, adduct ions of negative mode as -H, +HCOO and -2H. Data scaling and normalization were further processed using metaX 23.
Flow cytometry
Single-cell suspensions were prepared from the ear skin of mice. Skin lesions were dissected and digested with 2.0 mg/ml collagenase IV (Sigma-Aldrich Cat. V900893) and 1.0 mg/ml dispase II (Sigma-Aldrich Cat. D4693) for 60 min at 37°C. All the single-cell suspensions were filtered through 40-micron pores (BD Cat. 352340) and then were stained with fluorophore-conjugated antibodies. The cells were incubated with live/dead stain (Zombie Aqua™ Fixable Viability Kit, BioLegend Cat. 432102) and Fc block (BioLegend Cat. 101302). The cells were then washed and stained with antibodies specific for CD45 (30-F11, BioLegend Cat. 103116), CD4 (RM4-5, BioLegend Cat. 100510 or Cat. 100438), and CD25 (PC61, BioLegend Cat. 102012). For intracellular cytokine staining, the cells were stimulated for 5 hours with PMA (10 ng/ml, Sigma Cat. 79346), ionomycin (1 mg/ml, TOCRIS Cat. 56092-8-0), and GolgiPlug (1000×, BD Cat. 555029). For intracellular and intranuclear staining, the cells were fixed and permeabilized using the Foxp3/Transcription Factor Staining Buffer Kit (eBioscience Cat. 00-5523-00) according to the manufacturer’s protocol. The cells were intracellularly stained with antibodies specific for IFN-γ (XMG1.2, BioLegend Cat. 505810), or IL-17A (TC11-18H10, BD Biosciences Cat. 559502), or Foxp3 (FJK-16s, eBioscience Cat. 12-5773-82). For analysis of cell proliferation, cells were prestained with 5 mM CFSE (eBioscience Cat. 65-0850-84) before activation, according to the manufacturer’s instructions. Cell proliferation was evaluated by measuring CFSE dilution. All samples were acquired on a Cytek Dxp Athena flow cytometer. Analysis of acquired data was performed with FlowJo software. Analysis of the stained populations was performed by gating on single, live cells.
Cell culture
Human primary keratinocytes (KCs) were isolated from the foreskin of healthy donors obtained after circumcision. The skin samples were cut into small pieces and then incubated in dispase II (Sigma-Aldrich Cat. D4693, 2 mg/mL) overnight at 4 °C. The next day, the epidermis was removed from the dermis, subsequently digested with 0.05% trypsin for 10 min at 37 °C and filtered through 40-micron filters (BD Cat. 352340). After centrifugation at 1000 rpm for 5 min, the cells were cultured in the serum-free basal medium supplemented with growth factors (PromoCell Cat. C-20011). The medium was refreshed every day, and the cells were sub-cultured according to the cell fusion. Mouse naïve CD4+ T cells were obtained from total mouse splenocytes by magnetic cell sorting with the CD4+CD62L+ T Cell Isolation Kit following the manufacturer’s instructions (Miltenyi Biotech Cat. 130-106-643).
Keratinocyte proliferation
Keratinocyte proliferation was analyzed by CCK8 cell proliferation assay using the Cell Counting Kit-8 (CCK-8) (Bimake, B34302). In brief, keratinocytes were seeded into 96-well plates at 5000 cells/well and cultured for 24, 48, or 72 hours in the presence of different concentrations of C1P or vehicle control. Then CCK8 (10 µl) mixed with 100 µl medium was added to each well, and plates were incubated at 37 °C for 2 hours. Then, the absorbance was measured at 450 nm using an automated microplate reader (Beckman, USA). Five replicates per sample were analyzed.
Quantitative Real-Time PCR (RT-qPCR)
RNA was extracted from cells using TRIpure Reagent (Bioteke Cat. RP1001) according to manufacturer’s instructions. RNA was converted to cDNA using HiScript II Q RT SuperMix for qPCR (+gDNA wiper) (Vazyme Cat. R223-01), and gene expression was determined by RT-qPCR using the UltraSYBR One-Step RT-qPCR Kit (CWBIO Cat. CW0659) on a 7500 Fast thermocycler (Applied Biosystems). The relative expression of target genes was confirmed using the quantity of target gene/quantity of β-Actin. The fold change of gene expression was calculated by 2–(ΔCt experimental group–ΔCt control group), which normalized to the control group. All primer sequences used for RT-qPCR were as follows: IL-6: forward, 5′ - CACTGGTCTTTTGGAGTTTGAG -3′, and reverse, 5′ -GGACTTTTGTACTCATCTGCAC -3′, IL-1β: forward, 5′ - CTCCACCTCCAGGGACAGGATATG -3′, and reverse, 5′ - TCATCTTTCAACACGCAGGACAGG -3′, IL-17A: forward, 5′ - GAGATATCCCTCTGTGATCTGG -3′, and reverse, 5′ - GACAGAGTTCATGTGGTAGTCC -3′, β-Actin: forward, 5′ - AGAGCTACGAGCTGCCTGAC -3′, and reverse, 5′ - AGCACTGTGTTGGCGTACAG -3′. All primers were purchased from Sangon Biotech.
In vitro mouse T cell differentiation
Purified mouse naïve CD4+ T cells were stimulated with plate-bound anti-CD3 (5 μg/ml, eBioscience Cat. 145-2c11) and anti-CD28 (2 μg/ml, BD Cat. 553294) in RPMI 1640 (BI Cat. 01-100-1ACS), supplemented with 10% FBS (Gibco Cat. 16140071), 1000× 2-mercaptoethanol (Gibco Cat. 21985023), 100× HEPES Buffer (Gibco Cat. 15630080), 100× sodium pyruvate (Gibco Cat. 11360-070), 100× nonessential amino acid solution (Gibco Cat. 11140050) and 100× penicillin/streptomycin nystatin solution (Geneview Cat. B103010106). To assess Th1 cells polarization, naïve CD4+ T cells were activated with IL-12 (20 ng/ml, R&D Systems Cat. 419-ML-010) in the presence of anti-IL-4 (5 μg/ml, BioXcell Cat. BE0045). To assess Th17 cells differentiation, naïve CD4+ T cells were stimulated with IL-6 (30 ng/ml, R&D Systems Cat. 406-ML-005), IL-1β (10 ng/ml, PeproTech Cat. 211-11B-10), IL-23 (20 ng/ml, R&D Systems Cat. 1887-ML-010), and TGF-β (2 ng/ml, R&D Systems Cat. 7666-MB-005) in the presence of anti-IFN-γ (10 μg/ml, BioXcell Cat. BE0054) and anti-IL-4 (10 μg/ml, BioXcell Cat. BE0045). We cultured the cells culture in 96-well plates in a total volume of 0.2 ml/well of culture medium with 2 × 105 naïve CD4+ T cells. For C1P stimulation, C1P (Avanti Cat. 860533P) was dissolved at a concentration of 1 mM in normal saline and stored at -20 °C until use. The cells were incubated at 37 °C for 3 days under polarization conditions in the presence of C1P or vehicle control.
Statistical analysis
Metaboanalyst 24 was used for biomarker analysis. All figures were drawn using corresponding R packages. All quantified lipids were clustered into time patterns as follows. The mean abundance level for each lipid was determined at each of the three time points (N0, Q3 and Q6). Each row was then centered and clustered using k-means clustering with different cluster numbers on a Euclidean distance matrix using the R package TCseq (version 1.10.0). Appropriate cluster number was determined by assessing the stability of the clustering in terms of within-cluster sums of squares over 100 iterations of the clustering.
Clustering was performed using consensus clustering algorithm implemented in ConsensusClusterPlus R package 25 with subsampling psoriasis samples over 1000 iterations and 10 as maximum k (number of clusters) using differential lipids identified between psoriasis and control tissue samples. Before classification, samples were z-score normalized followed by lipid expression z-score normalization. Clustering results were evaluated visually in the resulting consensus matrix (one matrix per each k) as well as quantitatively using the cumulative distribution function of the area under the curve for each k.
Statistical analyses were performed on GraphPad Prism 8.0 software. Data are expressed as means ± SEM. A Student’s t test was used to compare two conditions, and an analysis of variance (ANOVA) with Bonferroni or Newman-Keuls correction was used for multiple comparisons. The level of significance was defined as p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. The error bars depict the SEM.