Mice
Female DBA/1J mice (8 weeks old, 20–24 g) were obtained from Jiangsu GemPharmatech Co. Ltd. (Nanjing, China). The mice were housed in specific-pathogen-free (SPF) facilities at the Laboratory Animal Center of Sun Yat-sen University and maintained at 25°C, 75% humidity, and a 12-hour light/dark cycle with ad libitum access to food and water. The animals had ad libitum access to standard chow and water. Prior to the experiment, all the mice were acclimatized for at least 7 days. All animal procedures complied with the Guide for the Care and Use of Laboratory Animals (8th Edition, National Academies Press, 2011, ISBN: 0309154006) and were approved by the Sun Yat-sen University Committee on Laboratory Animal Management and Ethics (No. SYSU-IACUC-2023-000203).
Human synovial samples
Synovial samples from the knee joint were collected from patients with osteoarthritis (OA) or rheumatoid arthritis (RA) during total knee arthroplasty (TKA) surgery. The study protocol was approved by the Ethics Committee of The Eighth Affiliated Hospital, Sun Yat-sen University (Shenzhen, China; Ethics No. 2022-035-03), and written informed consent was obtained from all participants.
Cells
Six healthy donors (HDs) and six rheumatoid arthritis (RA) patients were recruited for this study. The protocol was approved by the Ethics Committee of The Eighth Affiliated Hospital, Sun Yat-sen University (Ethics No. 2022-035-03), and written informed consent was obtained from all participants. Human whole blood samples were collected from HDs and RA patients.
Single-cell RNA sequencing data processing and analysis
Synovial macrophage single-cell transcriptome datasets and clinical metadata from 17 rheumatoid arthritis (RA) patients, 4 undifferentiated arthritis patients, 6 osteoarthritis (OA) patients, and 7 healthy control subjects were obtained from the European Bioinformatics Institute (EBI, E-MTAB-8322) 34 and the Gene Expression Omnibus (GEO, GSE216651 35, GSE152805 36). The raw gene expression matrices were processed via Seurat v5.0.3 package 37 in R v4.3.3, which involves quality control, normalization, and batch correction. The cells were filtered on the basis of gene count (> 300 genes), feature complexity (< 3500 unique molecular identifiers [UMIs]), and mitochondrial gene expression (< 20% of total UMIs) to exclude dead cells and multiplets.
Data from GSE216651 and GSE152805 were first integrated via the IntegrateLayers function with canonical correlation analysis (CCA), followed by clustering with FindClusters (resolution = 0.4) to identify major cell populations. Macrophages were isolated via marker genes (FCGR1A, ITGAM, CD14, LYZ) and further integrated with E-MTAB-8322 data via IntegrateLayers with reciprocal principal component analysis (RPCA), followed by reclustering at a lower resolution (0.3) to refine the subpopulations. Differential gene expression across macrophage subgroups was identified via FindAllMarkers (min.pct = 0.3, logFC threshold = 0.3), with subgroup distributions visualized via ggplot2 v3.4.2. The functional enrichment of DEGs in CD48highS100A12+ macrophages was performed via clusterProfiler v4.8.1 38 and GO.db v3.17.0, with the results visualized in ggplot2. To infer transcription factor (TF) gene regulatory networks at single-cell resolution, we employed pySCENIC v0.12.1, a high-performance Python implementation of the SCENIC pipeline 39, within a Python 3.7 environment. An 8,000-cell submatrix was first extracted from the integrated single-cell transcriptome to identify coexpression modules between TFs and their putative target genes.
Regulon (gene-TF regulatory network) construction utilized the GRNBoost2 algorithm, which models coexpression relationships across cells to predict TF-target interactions. Direct DNA-binding targets of TFs were validated via motif enrichment analysis via cisTarget databases (resources.aertslab.org/cistarget).
Regulon activity scores (RASs) were calculated for each cell to quantify the enrichment of TF-target gene coexpression, with each regulon comprising a TF and its motif-confirmed direct targets. Cell type-specific regulons were identified via the regulon specificity score (RSS), which was computed via the Jensen–Shannon divergence (JSD) algorithm, to measure regulon activity uniqueness across subpopulations.
To mitigate sampling bias, the Seurat-integrated dataset was randomly subsampled three times, with SCENIC analysis repeated for each subsample. The results were intersected to derive consensus regulons, ensuring robustness. Finally, UMAP dimensionality reduction (umap-learn v0.5.5) was performed on the basis of TF regulon activity scores, with visualizations generated via ggplot2 v3.4.2 to map regulatory heterogeneity across macrophage subclusters.
ChIP-seq data processing and analysis
IRF7 ChIP-seq raw reads from mouse bone marrow-derived macrophages stimulated with 100 ng/mL lipopolysaccharide (LPS) for 24 hours (GSE62697 40) were obtained from the GEO database. The raw reads were trimmed via TrimGalore v0.6.1 with the following parameters: -q 25 --phred33 --length 36 -e 0.1 to remove low-quality sequences and adapter contaminants. Clean reads were aligned to the mouse genome (mm10) via Bowtie2 v2.4.4 41, generating SAM files that were subsequently converted to BAM files via SAMtools v1.17 for downstream analysis. 42
Peak calling was performed via MACS2 v2.2.7.143, with unstimulated (t0) samples serving as the negative control to identify differential IRF7 binding sites. The enrichment profiles of the ChIP-seq peaks were visualized via deepTools v3.5.2 44, specifically the computeMatrix and plotHeatmap functions, to assess the binding intensity around transcription start sites (TSSs). Genomic annotation of peaks, including promoter and enhancer regions, was conducted via the Chipseeker v1.36.0 R package. 45
To identify direct IRF7 target genes, we intersected ChIP-seq-derived binding genes with IRF7 regulon-targeted genes from pySCENIC analysis, defining high-confidence regulatory interactions. The functional enrichment of these candidate genes was performed via the clusterProfiler package, which focuses on Reactome pathway analysis to characterize the biological processes regulated by IRF7 in inflammatory macrophages. 46
Bulk RNA-seq data analysis
Bulk transcriptome datasets (GSE154346 21, GSE130011 23) from a macrophage line under unstimulated (M0) and M1-polarized conditions were downloaded from the GEO database. To validate IRF7-mediated gene regulation, we generated a heatmap of IRF7 and its predicted downstream targets (e.g., PTGS2, CXCL10, NFKB1, and IL1B) via the pheatmap v1.0.12 package.47
Immunohistochemical staining of clinical knee synovial samples
Formalin-fixed, paraffin-embedded (FFPE) 4-µm-thick knee synovial tissue sections were deparaffinized in xylene followed by rehydration through a series of graded ethanol solutions. Antigen retrieval was executed via microwave-mediated heating at 95°C for 20 minutes in 10 mM citrate buffer (pH 6.0) to expose the target epitopes. Immunohistochemical staining was performed via an SP Rabbit & Mouse HRP Kit (DAB; Catalog No. CW2069S, CWBIO, China) following the manufacturer’s standardized protocol. Endogenous peroxidase activity was quenched with 3% hydrogen peroxide in methanol for 15 minutes, followed by blocking with 5% bovine serum albumin (BSA) in Tris-buffered saline (TBS) for 1 hour at room temperature (RT). Primary antibodies against CD68 (macrophage marker; 25747-1-AP, Proteintech; 1:200 dilution) and S100A12 (inflammatory marker; 16630-1-AP, Proteintech; 1:100 dilution) were applied overnight at 4°C. After being washed with TBS-T (0.1% Tween-20 in TBS), the sections were incubated with species-matched horseradish peroxidase (HRP)-conjugated secondary antibodies (CW2069S, CWBIO) for 45 minutes at RT. Signal development was performed using 3,3'-diaminobenzidine (DAB) substrate solution (10 minutes), and the nuclei were counterstained with hematoxylin for 2 minutes. Images were acquired via a Leica microscope (Leica Microsystems, Germany) at 400× magnification.
Cell culture and transfection
Human CD14 + monocytes were isolated from whole blood through density gradient centrifugation via Ficoll-Paque PLUS medium (GE Healthcare Life Sciences), followed by magnetic-activated cell sorting with anti-human CD14 monoclonal antibody-conjugated magnetic beads (Miltenyi Biotech, Germany), which was performed strictly according to the manufacturers’ protocols. The isolated cells were maintained in RPMI 1640 culture medium supplemented with 10% heat-inactivated fetal bovine serum (FBS), 2 mM L-glutamine, and a sterile antibiotic mixture (100 U/mL penicillin and 100 µg/mL streptomycin). Cultures were incubated at 37°C in a humidified atmosphere containing 5% CO₂, with half-medium replacement every 3 days to maintain cell viability and proliferation. For macrophage differentiation, CD14⁺ monocytes were seeded into culture wells and incubated with 25 ng/mL macrophage colony-stimulating factor (M-CSF; Sigma‒Aldrich, USA) for 5 days.
For M1 polarization, differentiated macrophages were stimulated with 50 ng/mL lipopolysaccharide (LPS) and 20 ng/mL interferon-γ (IFN-γ; Sigma‒Aldrich, USA) for 24 hours. IRF7 knockdown was achieved by transfecting M1-polarized macrophages with a siRNA targeting IRF7 (siRNA-IRF7; iGEBIO, Guangzhou, China) via Lipofectamine™ RNAiMAX (Invitrogen) according to the manufacturer’s protocol. A nontargeting siRNA served as the negative control (NC). One day before transfection, M1 macrophages were seeded in 6-well plates at 3 × 10⁵ cells/well. At 60–80% confluence, 5 µL of siRNA-IRF7 or siRNA-NC was diluted in 45 µL OPTI-MEM (Sigma‒Aldrich, USA), mixed with Lipofectamine™ RNAiMAX reagent, and incubated for 20 minutes at 25°C before being added to the cells for 24 hours.
RNA isolation
Total RNA was isolated from M1-polarized macrophages via TRIzol reagent (Invitrogen, USA) according to the manufacturer’s protocol. The RNA concentration and purity were quantified and assessed via spectrophotometric analysis. The integrity of the total RNA was assessed via agarose gel electrophoresis. The qualified total RNA was stored at -80°C and utilized for subsequent analysis.
Real-time quantitative polymerase chain reaction (RT‒qPCR)
Complementary DNA (cDNA) was synthesized from total RNA via the PrimeScript RT Kit (TaKaRa) according to the manufacturer’s protocol with a Bio-Rad T100 Thermal Cycler. RT‒qPCR was performed on an Applied Biosystems 7500 Real-Time PCR System using SYBR Premix Ex Taq (TaKaRa). The thermal cycling protocol comprised an initial denaturation step at 95°C for 30 seconds, followed by 40 cycles of denaturation at 95°C for 5 seconds and annealing/extension at 60°C for 20 seconds. Each sample was analyzed in technical triplicate, and the mean mRNA expression levels were calculated.
Melting curve analysis was conducted to validate specific target amplification, ensuring single-product formation without primer dimers. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) served as the internal reference gene, and relative gene expression was quantified via the 2^(-ΔΔCt) method. The forward and reverse primers for the target genes are listed in Table 1.
Table 1
RT‒qPCR real-time primer sequences.
Gene | Amplicon Size (bp) | Forward primer (5’→3’) | Reverse primer (5’→3’) |
|---|
IRF7 | 84 | GCTGGACGTGACCATCATGTA | GGGCCGTATAGGAACGTGC |
NFKB1 | 104 | AACAGAGAGGATTTCGTTTCCG | TTTGACCTGAGGGTAAGACTTCT |
PTGS2 | 94 | CTGGCGCTCAGCCATACAG | CGCACTTATACTGGTCAAATCCC |
IL1B | 132 | ATGATGGCTTATTACAGTGGCAA | GTCGGAGATTCGTAGCTGGA |
TNF | 220 | CCTCTCTCTAATCAGCCCTCTG | GAGGACCTGGGAGTAGATGAG |
CXCL10 | 198 | GTGGCATTCAAGGAGTACCTC | TGATGGCCTTCGATTCTGGATT |
Western blotting
The cells were lysed on ice for 30 minutes in radioimmunoprecipitation assay (RIPA) buffer (Sigma‒Aldrich, USA) supplemented with 1% protease and phosphatase inhibitor cocktail (Roche). The cell lysate was then centrifuged at 12,000 × g at 4°C for 30 minutes. The protein concentration in the supernatant was determined via a Bicinchoninic Acid (BCA) Protein Assay Kit (Sigma‒Aldrich, USA) following the manufacturer's protocol.
Equal amounts of protein were mixed with sodium dodecyl sulfate (SDS) loading buffer. The samples were then separated via 10% SDS-polyacrylamide gel electrophoresis (SDS‒PAGE) and electrotransferred onto a polyvinylidene fluoride (PVDF) membrane (Merck Millipore).
The membrane was blocked with Tris-buffered saline with Tween 20 (TBST) containing 5% skim milk powder at room temperature for 1 hour. The membrane was subsequently incubated overnight at 4°C with primary antibodies specific to glyceraldehyde 3-phosphate dehydrogenase (GAPDH, CST 5174), interferon regulatory factor 7 (IRF7, Proteintech Group, 22392–1 - AP), C - X - C motif chemokine ligand 10 (CXCL10, 10937-1-AP, Proteintech Group), prostaglandin-endoperoxide synthase 2 (PTGS2, 27308-1-AP, Proteintech Group), interleukin 1 beta (IL-1β, 26048-1-AP, Proteintech Group), and nuclear factor kappa B subunit 1 (NFKB1, 14220-1-AP, Proteintech Group).
After incubation with primary antibodies, the membrane was washed three times with TBST. Then, the samples were incubated at room temperature for 1 hour with an appropriate horseradish peroxidase (HRP)-conjugated secondary antibody (diluted 1:2000, Santa Cruz Biotechnology). The membrane was subsequently washed three more times with TBST. The chemiluminescent signal was detected via Immobilon Western Chemiluminescence HRP Substrate (Merck Millipore). Semiquantitative analysis of protein expression was performed via ImageJ software (National Institutes of Health, USA).
Collagen-induced arthritis (CIA) mouse model
The mice were randomly assigned to four experimental groups (n = 6 per group, for a total of 24) via a random number generator in R v4.3.3. Allocation was stratified to ensure balanced group distribution across cage positions and housing racks. Treatments and measurements were conducted in a randomized order to avoid systematic bias. The sample size was based on previous studies showing sufficient power (80%) to detect differences in clinical arthritis scores and histological inflammation, with an effect size of 1.2 and a significance level of 0.05. The group allocations were as follows:
In the si-IRF7 group, the mice received intra-articular injections of IRF7-specific siRNA (3 nmol/mouse/dose; iGEBIO, China) into the ankle joint once weekly from day 18 to day 39 after initial immunization.
Si-mock group: Mice were administered mock nontargeting siRNA via the same dosing and administration schedule.
CIA model group (positive control): Mice underwent collagen-induced arthritis (CIA) induction but received no therapeutic intervention.
Healthy control group (negative control): Mice were neither immunized nor treated.
All interventions and assessments were performed by researchers blinded to group allocation to reduce bias.
CIA was induced in female DBA/1J mice (8 weeks old, 20–24 g; GemPharmatech, China) following a validated protocol. Briefly, chicken type II collagen (Chondrex, USA) was dissolved in 0.05 M acetic acid to a final concentration of 2 mg/mL and emulsified 1:1 (v/v) with complete Freund’s adjuvant (CFA; containing 4 mg/mL Mycobacterium tuberculosis; Chondrex, USA). The mice received an initial intradermal immunization with 100 µg of collagen emulsion at the base of the tail. On day 21, a booster injection of 100 µg of chicken collagen in incomplete Freund’s adjuvant (Chondrex, USA) was administered intraperitoneally to elicit arthritis, which typically developed 7–10 days post-booster.
From day 18 post-initial immunization, the mice underwent examinations every three days to evaluate the onset of arthritis. Paw thickness was measured via a dial caliper. The primary outcome was joint inflammation, quantified by paw swelling, whereas the secondary outcomes included histopathological changes and bone erosion.
On day 43, all the mice were humanely euthanized via CO₂ inhalation followed by cervical dislocation. Ankle joints were harvested and fixed for microcomputed tomography (µCT) analysis via a Siemens Inveon system to evaluate synovial inflammation and bone erosion.
The inclusion and exclusion criteria were established a priori. The mice were included in the study if they were healthy, age-matched (8–10 weeks old), and free of signs of illness at the start of the experiment. Animals were excluded if they exhibited preexisting health issues, failed to recover from anesthesia, or showed abnormal behavior unrelated to experimental treatment. During analysis, data points were excluded only if technical errors occurred (e.g., sample loss, poor tissue quality for histology) or if the animals were euthanized for humane reasons unrelated to the experimental outcomes. No animals were excluded from the study. All the mice completed the experimental protocol. No outliers were removed.
Flow cytometric analysis of mouse ankle tissues
Fresh ankle joint samples were surface disinfected with 75% ethanol and processed under sterile conditions on ice. The tissues were dissected into small fragments via a sterile scalpel, washed twice with ice-cold sterile PBS, and digested in 1.5 U/mL Dispase II solution (Sigma‒Aldrich, USA) at 37°C for 60 minutes with gentle agitation. The digested tissues were strained through a 70-µm sterile cell strainer (BD Falcon) to obtain single-cell suspensions, which were subsequently centrifuged at 300 × g for 5 minutes at 4°C. The cell pellets were subsequently resuspended in RPMI 1640 medium (Gibco) supplemented with 10% FBS.
For M2 macrophage quantification, the cells were first stained with a Zombie Violet™ Fixable Viability Kit (DAPI conjugate; Biolegend, USA) to exclude dead cells, followed by surface staining with a FITC-conjugated anti-mouse F4/80 antibody (123107, Biolegend) for 30 minutes at room temperature in the dark. After being washed, the cells were fixed and permeabilized via a fixation/permeabilization solution (Invitrogen, USA) according to the manufacturer’s protocol and then incubated with an Alexa Fluor 647-conjugated anti-mouse CD206 antibody (565250, BD Biosciences) for 60 minutes at 4°C in the dark.
For regulatory T (Treg) cell phenotyping, the cells were stained with the Zombie Violet™ Fixable Viability Kit (BV421 conjugate; Biolegend, USA) and the surface markers APC-conjugated anti-mouse CD3 antibody (100235, Biolegend, USA) and FITC-conjugated anti-mouse CD4 antibody (100405, Biolegend, USA) for 15 minutes at room temperature. Following fixation/permeabilization, intracellular staining for FOXP3 was performed with a PE-conjugated anti-mouse FOXP3 antibody (320007, Biolegend, USA) for 60 minutes at 4°C in the dark.
All the stained samples were washed twice with PBS, resuspended in 200 µL of 1% BSA-PBS, and transferred to FACS tubes. Flow cytometry data were acquired on a BD FACS Celesta cytometer (BD Biosciences) and analyzed via FlowJo software (version 10; BD Biosciences), with gating strategies validated by isotype controls and fluorescence-minus-one (FMO) samples.
Histological and immunohistochemical analysis of mouse ankles
After being euthanized, the ankle joints of the mice were dissected and fixed in 10% neutral-buffered formalin (pH 7.4) for 48 hours at room temperature. The samples were decalcified in 10% ethylenediaminetetraacetic acid (EDTA) decalcifying solution (pH 7.2) for 72 hours with gentle agitation, processed into 5-µm paraffin sections and stained with hematoxylin‒eosin (H&E) for morphological evaluation. Synovial inflammation and hyperplasia were assessed by a blinded observer via a validated scoring system and evaluated on a 0–3 scale: synovial inflammation (0 = thin synovium (1–2 cell layers), no inflammatory cell infiltration; 3 = severe thickening (polypoid hyperplasia), dense inflammatory cell aggregation, with villi formation or fibrin exudation). 48
For IF analysis, the sections were blocked in 5% bovine serum albumin (BSA) in PBS for 1 hour at room temperature and then incubated overnight at 4°C with the following primary antibodies: a rat recombinant anti-CD68 antibody (ab53444, Abcam, USA; 1:200) and a rabbit polyclonal anti-S100A12 antibody (16630-1-AP, Proteintech, China; 1:100). After washing, the sections were incubated with the following fluorescently conjugated secondary antibodies: goat anti-rabbit IgG H&L (Alexa Fluor® 488; ab150077; Abcam; 1:500) and goat anti-rat IgG H&L (Alexa Fluor® 555; ab150158; Abcam; 1:500) for 1 hour at room temperature in the dark. The slides were mounted with mounting medium and antifading agent (with DAPI) (S2110, SolarBio, China) and imaged via a fluorescence microscope (Leica Microsystems, Germany).
IHC was performed via a streptavidin‒biotin peroxidase kit (CW2069S, Cowin Bio, China) according to the manufacturer’s protocol. The sections were pretreated with citrate buffer (pH 6.0) for antigen retrieval and then incubated overnight at 4°C with the following primary antibodies: rabbit anti-CD68 antibody (28058-1-AP, Proteintech; 1:200), rabbit anti-S100A12 antibody (16630-1-AP, Proteintech; 1:100), rabbit anti-IRF7 antibody (22392-1-AP, Proteintech; 1:100), and rabbit anti-PTGS2 antibody (27308-1-AP, Proteintech; 1:100). Following washes, the sections were incubated with HRP-conjugated secondary antibodies and developed with 3,3’-diaminobenzidine (DAB) substrate. The slides were counterstained with hematoxylin, dehydrated, and mounted with DPX mounting medium (Sigma‒Aldrich, USA).
IHC slides were imaged via a microscope (Leica Microsystems, Germany) at ×200 magnification. The average optical density (AOD) of positive staining was quantified via ImageJ software (NIH, USA) with the IHC Toolbox plugin, normalized to background levels and expressed as the mean ± SD across three independent fields per sample.
Statistical analysis
All in vitro experiments (RT‒qPCR and Western blotting) were performed in triplicate, and the data are presented as the means ± standard deviations (SDs). Statistical significance was determined via Student's t test for two groups and one-way analysis of variance (ANOVA) followed by a Bonferroni correction for three or more groups. Statistical analyses were conducted via R and GraphPad Prism 9. A P value < 0.05 was considered statistically significant.