Bioinformatic analysis
Transcriptome datasets were obtained from NCBI Gene Expression Omnibus libraries GSE145725 and GSE226331, where relative gene expression levels were calculated using Cufflinks v2.2.1 in fragments per kilobase per million reads (FPKM). Transcript annotation files were obtained from the Gencode v22 Comprehensive Gene Annotation gtf file (http://www.gencodegenes.org). EdgeR, a Bioconductor software package, was used to identify differentially expressed genes (DEGs) between normal skin tissue and fibrotic keloids for gene ontology (GO) pathway analysis, the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and to draw Volcano Plot and heat maps. A PPI network was created utilizing the Search Tool for Retrieval of Interacting Genes (STRING) database (https://string-db.org/), in which homo sapiens was the only species considered. The PPI network was constructed by mutual mapping of disease targets and drug targets, incorporating the mapping results into the STRING database.
Collection of hypertrophic scar and normal skin
Samples of hypertrophic scar were obtained from 3 donors (at an average age of 42 years) undergoing resection of scar in accordance with procedures approved by the Ethics Committee of Shanghai Ninth People's Hospital. Each participant was asked to sign an informed consent form before surgery. Hypertrophic scar and normal skin were identified upon characteristics of clinical diagnosis and pathology. Following harvest of samples, the tissues were processed for histology, fibroblast culture and qRT-PCR analysis, respectively.
Creation of skin fibrosis in mice
Fifteen 8-week-old male BALB/c mice (weighing 30–35g) were purchased from Shanghai SLAC Laboratory Animal Co.Ltd., and maintained under specific pathogen-free (SPF) conditions in animal center of Tongji University School of Medicine. All procedures were approved by the Animal Research Committee of Tongji University.
As reported previously[15], skin fibrosis was induced by subcutaneously injecting bleomycin (BLM, 5mg/kg) into a one-square-centimeter area of dorsal skin in mice every two days for four weeks (BLM group, n = 5). Administration of AcAc was performed by subcutaneous injection of AcAc (3mg/kg) in combination with bleomycin every other day for 4 weeks (AcAc group, n = 5). Injection of phosphate buffered saline (PBS: Invitrogen, USA) served as control (PBS group, n = 5). After 28 days, dorsal skin with subcutaneous adipose layer was harvested for histology and qRT-PCR.
Histological and immunofluorescent staining
Samples of human normal skin, hypertrophic scar and bleomycin-treated dorsal skin of mice were fixed with formalin and then embedded with paraffin, respectively. Following serial section in 4 µm, slices were deparaffined and dehydrated and subjected to H&E, Masson’s trichrome and Sirius Red staining. For immunofluorescent staining, primary antibodies against COL1A1 (1:500, PA5-29569, Thermo Fisher Scientific), BDH1 (1:100, 15417-1-AP, Proteintech), OXCT1 (1:600, 12175-1-AP, Proteintech), p-Smad2/3 (1:100, PA5-99378, Thermo Fisher Scientific) were incubated at 4°C overnight. Slides were washed and incubated for 1 h with Alexa Fluor 568-conjugated goat anti-rabbit IgG (1:500 dilution in 5% BSA). Nuclei were stained for 10 min at RT with 5 µg/ml Hoechst 33242 from Thermo Fisher Scientific. Images were obtained using a Fluoview FV10i confocal microscope or a BX53 fluorescence microscope from Olympus (Tokyo, Japan).
Cultivation of dermal fibroblasts
The freshly harvested dermal tissue was thoroughly washed in PBS after removing of epidermal and adipose tissues, and subjected to digestion with type I collagenase (Worthington, USA) at 37°C for one hour. Following digestion, cell pellets were obtained by centrifuging at 400×g for 10 minutes and resuspended in culture medium consisting of Dulbecco's modified Eagle's medium (DMEM; Gibco, USA), supplemented with 10% fetal bovine serum (FBS; Hyclone, USA), and 1% penicillin-streptomycin. The resuspended fibroblasts were cultured at 37°C with 5% CO2 for approximately 14 days. When reaching 80 ~ 90% confluence, cells were passaged and those at passage 3–5 were used in the following study.
Transition of fibroblasts to myofibroblasts was established by treating cells with TGF-β1 (10ng/ml) for 24 Hours as previously reported[16]. The effect of AcAc on cell phenotype transition was observed by adding AcAc (50µM) in culture medium, in the presence of absence of TGF-β1 (10ng/ml), respectively.
Knock-down of OXCT1 in vitro and in vivo
Adenoviruses engineered to express mouse OXCT1 shRNA (Ad-shOXCT1) were procured from Wei Nuo Sai Biology (Wuhan, China). These adenoviruses were propagated by infecting HEK-293a cells for a duration of 48 hours. Subsequently, the HEK cells were collected and lysed through a series of freeze-thaw cycles to release the adenoviruses, which were then subjected to CsCl density gradient ultracentrifugation for purification.
Fibroblasts at passage 1, reaching a confluence of 70–80%, were cultured in serum-free medium for 12 hours prior to transfection with either Ad-shControl or Ad-shOXCT1 at a multiplicity of infection (MOI) of 800 for 48 hours, thereby establishing stable fibroblasts with OXCT1 knockdown. Followed by a 24-hour culture period with or without TGF-β1 (10ng/ml) before undergoing further analysis.
Knock-down of OXCT1 in vivo was performed by subcutaneous injection of Ad-shOXCT1 (1×109 plaque forming units (PFU) in a total volume of 10 µl) once a week for 4 weeks with or without bleomycin (Ad-shOXCT1 group, Ad-shOXCT1 + BLM group, n = 5). Subcutaneous injection of empty adenovirus (Ad-shControl) was used as a control (Ad-shControl group, Ad-shControl + BLM group, n = 5). Four weeks post-initial injection, the mice were euthanized, and their skin was meticulously harvested for subsequent analysis.
Quantitative real-time PCR
Total RNA was extracted with TRIzol (R401, Vazyme) and subsequently reverse-transcribed into cDNA with HiScript II Q RT SuperMix (R222-01, Vazyme) in preparation for qPCR. Quantitative real-time polymerase chain reaction (qRT-PCR) was consequently conducted with an Applied Biosystems 7500 device with ChamQ SYBR qPCR Master Mix (Q331-02, Vazyme). 2-ΔΔCT approach was employed to calculate the relative expression levels, in turn normalized to β-actin. The sequences of all primers are illustrated in Table 1.
Table 1
Primers used in this study.
Target gene | Primer sequences (5’–3’) |
mβ-actin | 5'ATATCGCTGCGCTGGTCGTC3' 5'AGGATGGCGTGAGGGAGAGC3' |
mα-SMA | 5'GTCCCAGACATCAGGGAGTAA3' 5'TCGGATACTTCAGCGTCAGGA3' |
mCOL1A1 | 5’GCTCCTCTTAGGGGCCACT3’ 5’CCACGTCTCACCATTGGGG3’ |
mCOL3A1 | 5’CTGTAACATGGAAACTGGGGAAA3’ 5’CCATAGCTGAACTGAAAACCACC3’ |
mFN1 | 5’ATGTGGACCCCTCCTGATAGT3’ 5’GCCCAGTGATTTCAGCAAAGG3’ |
mMCT1 | 5'CCAGGCACCTGTGGTCAAC3' 5'GGTCTCGTGATTGGATCTGCT3' |
mBDH1 | 5'CGGCTAGTGGCAAAGCTATC3' 5'GTTGCAGACATTGAGCTGGA3' |
mOXCT1 | 5'CATAAGGGGTGTGTCTGCTACT3' 5'GCAAGGTTGCACCATTAGGAAT3' |
mACAT1 | 5'CAGGAAGTAAGATGCCTGGAAC3' 5'TTCACCCCCTTGGATGACATT3' |
hβ-actin | 5'CATGTACGTTGCTATCCAGGC3' 5'CTCCTTAATGTCACGCACGAT3' |
hα-SMA | 5’AAAAGACAGCTACGTGGGTGA3’ 5’GCCATGTTCTATCGGGTACTTC3’ |
hCOL1A1 | 5’GAGGGCCAAGACGAAGACATC3’ 5’ CAGATCACGTCATCGCACAAC3’ |
hCOL3A1 | 5’GGAGCTGGCTACTTCTCGC3’ 5’GGGAACATCCTCCTTCAACAG3’ |
hFN1 | 5’CGGTGGCTGTCAGTCAAAG3’ 5’AAACCTCGGCTTCCTCCATAA3’ |
hMCT1 | 5'GGCCAAGTGTGAGTTCTTCAA3' 5'GGCTCGATAATCGTGTCCCC3' |
hBDH1 | 5'GGCAGAAGTGAACCTTTGGG3' 5'GCAGTCCGAGAAAGCCTCTAC3' |
hOXCT1 | 5'GTTGGTGGTTTTGGGCTATGT3' 5'AGACCATGCGTTTTATCTGCTT3' |
hACAT1 | 5'ATGCCAGTACACTGAATGATGG3' 5'GATGCAGCATATACAGGAGCAA3' |
Western blot analysis
Cells were implanted into six-well plates and cultured until attachment. TGF-β1, Hcy, and THF were added correspondingly and incubated. The protein concentration was determined using the Pierce BCA protein assay kit (Thermo Fisher Scientific). The collected proteins were separated through SDS-PAGE and were electrotransferred to polyvinylidene fluoride (PVDF) membranes. After blocking, the membranes were reacted with the corresponding primary antibodies against BDH1(1:5000, 15417-1-AP, Proteintech), OXCT1(1:6000, 12175-1-AP, Proteintech), COL1A1 (1:1000, TA7001, Abmart), COL3A1 (1:1000, T510299, Abmart), α-SMA (1:1000, T55295, Abmart, Shanghai, China), Fibronectin (1:1000, T59537, Abmart), Smad2/3 (1:4000, PA5-17621, Thermo Fisher Scientific), p-Smad2/3 (1:1000, PA5-110155, Thermo Fisher Scientific), GAPDH(1:1000, TT0004, Abmart) overnight at 4°C. Finally, they were incubated with HRP-linked rabbit IgG antibody (1:3000,7074S/ 7076S, Cell Signaling Technology) for 1 hour. Signals were developed with SuperSignal West Dura substrate (Thermo Fisher Scientific) and imaged with a Gel Doc XR + device (Bio-Rad Laboratories, Inc.). Quantification analysis was performed with the Volume tool in ImageLab 6.0.1 (Bio-Rad), adjusted to GAPDH expression, and normalized to respective Ctrl samples to calculate fold change to Ctrl.
Quantitative and statistical analysis
Data were expressed as the mean ± standard deviation (SD). To determine the statistical significance of observed differences, we utilized one-way or two-way analysis of variance (ANOVA), followed by Tukey's post hoc test for multiple comparisons when appropriate. For direct comparisons between two groups, we employed unpaired two-tailed Student's t-tests. The levels of statistical significance were set at *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001. All calculations were performed using Prism version 9.0 (GraphPad Software, Inc.) on a Windows 11 operating system.