2.1 Animal welfare statement and generation of a zymosan A-induced interstitial pneumonia arthritic mouse model.
All mouse protocols were approved by the Use of Laboratory Animal Committee of Ningxia University in accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals (KYLL-2024-1033). All mice were purchased from CLEA Japan Inc (Tokyo, Japan), and housed in a special pathogen-free facility with a 12/12h light/dark cycle with food and water ad libitum in the animal facility at Ningxia Medical University (Yinchuan, China). To generate an arthritic mouse model of interstitial pneumonia, 8-week-old SKG/jcl mice (n = 12) were intraperitioneally administered 7.5 mg of zymosan A (ZYM) (Alfa Aesar, Lancashire, UK) dissolved in 0.5 mL of phosphate buffered saline (PBS). The control mice (n = 12) were administered 0.5 mL of PBS. We examined body weight and paw size once a week. The arthritis score for each of the four paws was recorded, and the scoring system was as follows: 0 = normal joints; 1 = slight swelling or erythema of the ankle or midfoot; 2 = slight swelling of the ankle and foot; 3 = modest swelling and erythema; and 4 = severe swelling and erythema[20]. The maximum score was 16 for each mouse. The mice were euthanized via carbon dioxide (CO2) inhalation at the end of the experiment. The blood, knee joint, brain, heart, kidney, liver, spleen, colon and lung were harvested for pathological and molecular analysis at 8 and 16 weeks post-ZYM challenge.
2.2 Lung and joint histopathological analysis.
Mouse lung tissues were fixed in 4% paraformaldehyde (PFA) solution before they were either embedded in optimal cutting temperature (OCT) compound or dehydrated and processed for paraffin embedding. Five-micrometer-thick sections were prepared and stained using hematoxylin and eosin (H&E) staining, Masson’s trichrome staining, or other methods as described elsewhere[21]. The severity of interstitial injury and fibrosis in lung tissues was assessed by evaluating inflammatory cell infiltration and using the Ashcroft scale, with scores ranging from 0 to 4 for interstitial injury and 0 to 8 for fibrosis, on the basis of histological images of the lung as previously described[22].
Joint tissues were fixed in 10% formalin, decalcified with EDTA decalcification solution, embedded in paraffin, and cut into five-micrometer-thick slices. H&E and toluidine blue-O (TBO) staining were performed. To assess the degree of joint inflammation and joint destruction, inflammatory cell infiltration and cartilage damage were scored by two independent and blinded pathologists as previously described. Both joint inflammation and joint destruction were evaluated on a scale of 0–4, as previously described[23].
2.3 Measuring anti-CCP, RF, histone H3, cell-free DNA and MPO-DNA concentrations.
The concentrations of anti-citrullinated protein (anti-CCP) antibody, RF, MPO-DNA, citrullinated histone H3 (Cit-H3), and cell-free DNA were measured by an enzyme-linked immunosorbent assay (ELISA) with appropriate commercially available kits, following the manufactures’ instructions. The ELISA kits for detecting mouse anti-CCP and RF were obtained from Shanghai JiangLai Biotechnology Co.,Ltd (Shanghai, China), while the kits for detecting MPO-DNA, Cit-H3 and cell-free DNA were obtained from Shanghai HengYuan Biotechnology Co.,Ltd (Shanghai, China). The protein concentration in each sample was determined by comparison with the standard curve.
2.4 Evaluations with micro-CT scanning.
At 16 weeks postinjections, the knee joints of the mice were harvested and fixed in 4% paraformaldehyde. The samples were scanned via micro- CT (µCT 40; Scanco, Zurich, Switzerland) as described previously[24]. The scanner was set to a resolution of 10 µm with a voltage of 70 kV and an electric current of 114 µA. The region of interest (ROI) was defined to cover the entire subchondral bone in the tibial plateaus. The three-dimensional structural parameters analyzed included the bone volume (BV), bone volume/total tissue volume (BV/TV), trabecular thickness (Tb.Th), and trabecular separation (Tb.Sp).
2.5 RNA isolation and RT-PCR
Total RNA was isolated using the MiniBEST Universal RNA Extraction Kit (TaKaRa) according to the manufacturer’s instructions. Quantitative RT-PCR with SYBR Green Master Mix (TaKaRa, Dalian, China) was performed using the StepOnePlus RT-PCR System. The relative quantity of target mRNAs was calculated using the delta-delta Ct method. The primer sequences are displayed in Supplementary Table S1.
2.6 Western blotting.
Proteins were extracted from both cell lines and tissues using RIPA lysis buffer (Thermo Fisher Science, USA) under cold conditions, After centrifugation at 12000g for 10 min, the protein samples were subjected to immunoblotting. Protein concentration was determined using a BCA protein assay kit (Thermo Fisher Science, USA). Subsequently, the protein samples were loaded onto PAGE gels (Epizyme Biomedical Technology, Shanghai, China) and transferred to 0.22 µm Immobilon PVDF membranes (Millipore Sigma, USA). After blocking with 5% milk, the membranes were incubated in primary antibodies at appropriate dilutions overnight at 4℃. Next, secondary antibodies were applied at room temperature for approximately 1h, and the immunoreactivity was visualized using an ECL system (Thermo Fisher Scientific, Waltham, MA). Relative expression of target protein was measured semi-quantitatively by densometric analysis of blots. The intensity of blot area was calculated by optical densitometry using ImageJ Software version 2.0.0 (http://rsb.info.nib.gov/ij). The ration of the net intensity of each sample to that of housekeeping gene GAPDH served as an internal loading control. The values are reported as densitometric arbitrary units (A.U.). Relative target protein expression in the experimental group was divided by that in the control group to calculate the fold change in target protein abundance. Primary and secondary antibodies used for immunoblotting are listed in Supplementary Table S2.
2.7 Immunoflurescence staining.
Fixed lung tissues that were embedded in OCT compound were cryosectioned at a thickness of 8 µm for immunofluorescence (IF). Cryosections were air-dried at RT for 30 min, re-fixed in 4% PFA for 10 min, and finally permeabilized in 0.2% Triton X-100/PBS for 20 min at RT. Sections were then blocked by incubating them in 5% donkey serum in PBS for 1h before being probed with primary antibodies in dilutant buffer (1% donkey serum, 0.03% Triton X-100, and 1mM Cacl2 in PBS) overnight at 4℃. The sections were washed and then incubated with fluorescent dye-conjugated secondary antibodies at RT for 2h. The slides were washed in PBS three times for 5min, mounted with VECTASHIELD Antifade Mounting Medium with DAPI (Vector Laboratories), imaged under a Leica TCS SP2 A0BS Confocal System, and processed using Leica Confocal Software v.2.6.1 (Leica). Primary antibodies and secondary antibodies used for IF are listed in Supplementary Table S2.
2.8 Immunohistochemistry.
Paraffin-embedded tissue specimens were sectioned, deparaffinized with xylene, and rehydrated in decreasing concentrations of ethanol followed by 3 washes with 1×PBS. Antigen retrieval was performed in boiling sodium citrate buffer (10 mmol/L sodium citrate, 0.05% Tween-20, PH 6.0) or using Proteinase-K (20 µg/mL). The sections were then blocked with 3% fetal bovine serum in 1×PBS for 30 minutes, followed by overnight incubation at 4 ℃. After 3 washes with 1×PBS, the sections were further incubated with biotin or HRP-conjugated secondary antibodies for 30min at RT. Antibody staining was then developed with DAB substrate (Zsbio, Beijing, China) for HRP-conjugated secondary antibodies or incubated with an Avidin-Biotin Complex Staining Kit (ABC Kit; Vector Laboratories, Burlingame, CA) before DAB development as per the manufacturer’s instructions. Primary antibodies and secondary antibodies used for IHC are listed in Supplementary Table S2.
2.9 Human subjects.
The study and protocol were approved by the Ethics Committee for Conducting Human Research at General Hospital of Ningxia Medical University (KYLL-2024-1033). All patient participants provided written consent for the collection and analysis of their blood samples for publication, in accordance with the protocol (KYLL-2024-1033) outlined by the Ethics Committee for the Conduct of Human Research. The investigator of this study maintains human research records, including signed and dated consent documents for 10 years. All cases who fulfilled the 2010 American Collage of Rheumatology (ACR)/ European League Against Rheumatism (EULAR) and/or 1987 ACR revised criteria were collected from the outpatient rheumatology and respiratory clinic of General Hospital of Ningxia Medical University between November 2022 to July 2024[25]. Pulmonary involvement was evaluated in all patients by high-resolution computed tomography (HRCT). The assessments and subsequent reclassification of HRCT images according to the 2013 IIP classification[26] were independently analyzed by two senior radiologists at the General Hospital of Ningxia Medical University, with 8 and 15 years of experience, respectively. Ultimately, 25 RA-nonspecific interstitial pneumonia (NSIP), 24 RA-usual interstitial pneumonia (UIP), 12 RA-organizing pneumonia (OP), and 12 RA-other patterns patients were qualified and included in this study (Figure S1). Clinical characteristics and laboratory data were abstracted from medical records. The demographics of individuals involved in this study are outlined in Supplementary Figure S2.
2.11 Single-cell suspension preparation.
The preparation of single-cell suspension was followed by previous protocols[27]. Briefly, Lung lobes from the control and ZYM-treated mice were digested with collagenase-type IA and deoxyribonuclease I (Thermo Fisher Scientific Inc., Waltham, MA, USA), followed by trypsinization. The cell suspension was then filtered using a 40 µm cell strainer, washed, centrifuged, and finally resuspended in magnetic-activated cell sorting buffer (Miltenyi Biotec, Bergisch Gladbach, North Rhine-Westphalia, Germany) for single-cell RNA sequencing (ScRNA-Seq) analysis.
2.12 Construction of single-cell libraries and sequencing.
In brief, single-cell suspensions and beads were mixed with each other by adjusting co-encapsulation occupancy to 0.05. After individual droplets were collected, messenger RNA was reverse transcribed into complementary DNA, followed by cDNA amplification. Finally, a 3’gene expression library was prepared using the Chromium Next GEM Single Cell 3’Kit v3.1 (10× Genomics, Pleas-anton, CA, USA). Sequencing was conducted by using a NovaSeq 6,000 (Illumina, San Diego, CA, USA) by OE Biotech Co., Ltd. (Shanghai, China).
2.13 ScRNA-seq data quality control analysis.
Raw sequencing reads of mouse lung tissues were aligned to the mouse genome reference (GENCODE, mm10) and processed into a single-cell matrix via Cell Ranger (version 7.0.0) with the default parameters. Considering the presence of double droplets, empty droplets, dead cells, low-quality cells with the nFeature_RNA > 4000, percent.mt > 10 and nFeature_RNA < 1000 were excluded. A total of 37,735 cells were finally identified after filtration in the dataset. Subsequently, the data was normalized, the dimensionality was reduced, and cell clustering was performed using Seurat (v.4.1.1)[28]. Additionally, the clusters with relatively low gene numbers and absence of specific marker genes were also removed.
2.14 Integration and cell type annotation.
To integrate each sample and correct for batch effects, the top 2000 highly variable genes (HVGs) were identified using the Seurat function FindVariableGenes (mean.function = FastExpMean, dispersion.function = FastLogVMR). To remove the batch effects from the ScRNA-Seq data, the mutual nearest neighbors (MNN) method presented by Haghverdi et al., was performed with the R package batchelor (version 1.6.3)[29]. Graph-based clustering was performed to cluster cells according to their gene expression profile with the FindClusters function. Cells were visualized using a 2-dimensional uniform manifold approximation and projection (UMAP) algorithm with the RunUMAP function. Finally, the FindAllMarkers function (test.use = presto) was used to identify marker genes of each cluster.
2.15 Differentially expressed genes analysis.
Differentially expressed genes (DEGs) were selected using the function FindMarkers (test.use = presto). A P value < 0.05 and |log2fold change|>0.58 were set as the threshold for significantly differential expression.
2.16 Gene ontology enrichment analysis.
Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of differentially expressed genes (DEGs) were conducted using R using R (version 4.0.3), utilizing the hypergeometric distribution as the statistical method.
2.17 Microarray data and identification of hub genes.
The transcription profile dataset of lung tissues in RA-UIP, IPF-UIP and non-UIP control groups (GSE199152) was obtained from NCBI GEO database (http://www.ncbi.nlm.nih.gov/geo/). The platform is GPL16791, which included 27 subjects, of which 3 were RA-UIP patients, 20 were IPF patients, and 4 were non-UIP controls. DEGs were identified by comparing the expression values in lung tissues between RA-UIP and non-UIP control using the limma R package (http://www.bioconductor.org/packages/release/bioc/html/limma.html) in Bioconductor. The screening criteria set as adjusted P < 0.05 and FC > 1.
2.18 Cell culture and siRNA transfection.
The human promyelocytic (HL-60) and fetal lung (MRC-5) cell lines were purchased from Pricella Bioscience Inc, (Wuhan, China) and cultured in Iscoves’s Modified Dulbecco’s Medium (IMDM) (Gibco, USA) supplemented with 20% fetal bovine serum (FBS) (Gibco, USA) and 1% antibiotic-antimycotic solution at 37℃, 5% CO2. Transfection of HL-60 cells with either si-Control or si-Nr4a3 (Sangon, Shanghai, China) was performed using Lipofectamine 3000 regent (L3000015, Invitrogen) according to the instruction manual.
2.19 Statistics.
Except for the scRNA-seq data, experimental data were analyzed using two-tailed unpaired Student’s t-tests, one-way or two-way ANOVA by GraphPad Prism software (version 9.5, GraphPad Software, Inc., San Diego, CA) and “ggplot2”, “ggtext”, stats” and “car” packages in R (version 4.2.1). Data are expressed as mean ± standard errors, and differences were considered statistically significant when P < 0.05. All experiments, except for the scRNA-seq, were performed independently at least three times.