All animal experiments conducted in this study were approved by the Science and Technology Ethics Committee of Anhui University of Science and Technology (Approval No. GZ2025-048) and strictly adhered to the national standards for laboratory animal welfare (GB/T 35892-2018) of China. Figure 1 shows the study design.
Preparation of solutions
Preparation of silica suspension:
Crystalline silica (SiO₂, 80%, particle size 1–5 µm) was sterilized by autoclaving. The sterile silica powder was then dispersed in sterile saline to obtain a final concentration of 10 mg/mL, and the suspension was thoroughly homogenized using an ultrasonic homogenizer immediately before administration to ensure a uniform particle suspension.
Preparation of DMSO solution:
Under aseptic conditions, DMSO (≥99.9%) was diluted with sterile physiological saline to prepare an approximately 10% (v/v) DMSO solution, which was homogenized by ultrasonic agitation. Each mouse was injected with the DMSO solution at a dose of 0.9 g/kg, as previously described14,15.
Experimental animals
Thirty-two healthy male C57BL/6J mice aged 8–12 weeks were used in this study and acclimatized for 1 week before experimentation. Mice were randomly assigned to four experimental groups (n = 8 per group) using a computer-generated randomization method: control, DMSO, silica, and silica+DMSO groups. The sample size was determined based on previous studies using the same murine silicosis model and our preliminary experimental experience, which indicated that 8 animals per group were sufficient to detect histopathological and molecular differences between groups. Animals were housed under standard laboratory conditions with controlled temperature, a 12 h light/dark cycle, and free access to food and water. Humane endpoints were established before the study, and no animals met the predefined criteria for early euthanasia.
Silicosis model establishment and drug intervention
Silicosis was induced by a single intranasal instillation of 60 µL of sterile crystalline silica suspension (10 mg/mL). Following induction of anesthesia, mice were placed in the supine position. Once deep and slow breathing was observed, 60 µL of the silica suspension was carefully instilled dropwise into the nostrils, allowing the suspension to be inhaled spontaneously into the lungs16. Mice in the DMSO administration group received intraperitoneal injections of DMSO solution at a fixed dose twice weekly, with a 3–4 day interval between injections, for one consecutive month. Mice in the remaining groups were treated identically with an equal volume of sterile saline instead. Daily body weight monitoring was performed throughout the experiment. On the 30th day after modeling, all mice were euthanized, and lung tissues were collected under sterile conditions. Euthanasia was performed by cervical dislocation following intraperitoneal administration of tribromoethanol (0.2 mL/10 g) to induce anesthesia. Partial lung tissues were preserved in liquid nitrogen for transcriptome sequencing. A portion of tissues was fixed in 4% paraformaldehyde for 48 h prior to pathological staining, and the remaining tissues were stored at -80 °C for subsequent molecular experiments. All animal operations complied with institutional animal ethics guidelines.
Histopathological evaluation, including the Ashcroft fibrosis score, was independently performed by two investigators blinded to the experimental group allocation.
Pulmonary pathological examination
Fixed lung tissues were dehydrated, transparentized, and embedded in paraffin blocks, then sliced into 5 µm-thick sections.
hematoxylin and eosin (HE) staining: Sections were stained with hematoxylin for 3–5 min, differentiated in acidic solution, blued under running water, and counterstained with eosin for 5 min. After dehydration through graded ethanol, sections were cleared in xylene and sealed with neutral resin17.
Masson’s trichrome staining: Sections were stained sequentially with hematoxylin, Lichun Red–acid fuchsin, phosphomolybdic acid, and Bright Green. After acidic differentiation, sections were dehydrated, cleared in xylene, and mounted with neutral resin18,19.
Transcriptome sequencing
Mouse lung tissues from the four experimental groups were collected immediately after sacrifice and snap-frozen in liquid nitrogen. Total RNA was extracted using TRIzol reagent according to the manufacturer's instructions. RNA concentration and purity were assessed, and RNA integrity was evaluated using a Bioanalyzer with built-in software. Samples with RIN values greater than 8.0 were considered acceptable. Library preparation included poly(A) mRNA enrichment using Oligo(dT) beads, RNA fragmentation, cDNA synthesis, adapter ligation, and PCR amplification. Sequencing was performed on an Illumina platform, generating approximately 6.2–6.4 Gb of clean data per sample. All raw sequencing data generated in this study have been deposited in the Gene Expression Omnibus (GEO) database under accession number GSE311671.
Clean reads were aligned to the Mus musculus reference genome (mm8) using the STAR aligner. Gene expression levels were quantified as raw read counts. Differential expression analysis was performed using the limma package in R. Prior to linear modeling, raw count data were transformed using the voom function to estimate the mean-variance relationship and generate precision weights. Raw P values were adjusted for multiple testing using the Benjamini-Hochberg (BH) method, and genes with |log₂FoldChange| > 1.5 and false discovery rate (FDR) < 0.05 were considered differentially expressed. Functional enrichment analyses were performed using the clusterProfiler package in R. Gene ontology (GO) enrichment analysis was conducted to identify significantly enriched biological process (BP), cellular component (CC), and molecular function (MF) categories, whereas Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis was performed to identify significantly enriched signaling pathways. Enrichment significance was assessed using a hypergeometric test, and P values were adjusted for multiple comparisons using the Benjamini-Hochberg method. GO terms and KEGG pathways with an adjusted P value (padj) < 0.05 were considered significantly enriched20,21,22,23.
Western blot detection
Fresh mouse lung tissue was collected and homogenized with an appropriate volume of radioimmunoprecipitation assay (RIPA) lysis buffer. After centrifugation, the supernatant was collected, and protein concentration was determined using the BCA protein assay kit. The samples were then mixed with loading buffer and heated at 100 °C for 10 min. Protein separation was carried out on 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gels at a current of 200 V for 30 min. Subsequently, proteins were electrotransferred to polyvinylidene fluoride (PVDF) membranes with 300 mA constant current for 30 min. The membrane was blocked with 5% skim milk, then incubated overnight at 4 °C with primary antibodies against β-Actin (1:1500), IL-6 (1:1000), STAT3 (1:800), p-STAT3 (1:1500) and MMP12 (1:800). The following day, the membrane was washed with TBST and incubated with secondary antibody at room temperature for 1 h, followed by additional washes with TBST. Protein bands were visualized using a chemiluminescent substrate, and protein quantification was performed using ImageJ software.
Western blot analysis was performed using three independent biological replicates from each experimental group. Protein samples were prepared from individual mouse lung tissues, and each biological replicate represented one independent animal.
Statistical analysis
All experimental data were analyzed using GraphPad Prism 9.5. Normality and homogeneity of variance tests were conducted beforehand. Student’s t-test was used for two-group comparison, and one-way ANOVA combined with Tukey’s post-hoc test was adopted for multi-group comparison. P ≤ 0.05 was regarded as statistically significant.