Study design
We performed a cross-sectional study of patients treated for sJIA. The study was performed according to the Declaration of Helsinki principles and approved by the ethics review board of Chiba Children's Hospital, Chiba, Japan. (Approval number, 2020-022). Written informed consent was obtained from the study participants and/or their guardians.
Setting and participants
We recruited patients with sJIA at Chiba Children's Hospital in October 2020. Chiba Children's Hospital serves most pediatric patients with rheumatic disease in the area due to the small number of centers treating pediatric rheumatic diseases. The eligible participants were patients who had previously been diagnosed with sJIA based on the International League of Associations for Rheumatology (ILAR) classification criteria(28) and had undergone treatment at Chiba Children's Hospital between April 2013 and March 2020. The exclusion criteria were a lack of medical records, complications of other rheumatic diseases, complications of acute infection at time of serum collection, or other conditions that might induce inflammation, such as surgery, injury, and malignancy.
The patient's clinical symptoms and laboratory findings were obtained retrospectively from their medical records. The diagnosis of MAS was based on the 2016 EULAR/ACR/PRINTO classification criteria(29). The criteria for active-phase sJIA were typical symptoms, such as a fever, arthritis, hepatosplenomegaly, rash, and generalized lymphadenopathy, and increased CRP levels (>0.3 mg/dL).
Sample preparation for proteome analyses
Highly abundant serum proteins were depleted from 10 μl of serum using a High Select Top 14 Abundant Protein Depletion Mini Spin Column (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer's instructions. Depleted serum (50 μl) was diluted with 150 μl of 100 mM Tris-HCl pH 8.5 and 2% sodium dodecyl sulfate (SDS) and treated with 10 mM dithiothreitol at 50 °C for 30 min. The sample was then alkylated with 30 mM iodoacetamide in the dark at room temperature for 30 min and subjected to cleanup and digestion with single-pot solid phase-enhanced sample preparation (SP3)(30).
Two types of beads (hydrophilic and hydrophobic Sera-Mag Speed-Beads; Cytiva, Marlborough, MA, USA) were used for the SP3 method. These beads were combined at a 1:1 (v/v) ratio, rinsed with distilled water, and reconstituted in 500 mM Tris-HCl pH 7.0 at 10 μg solids/μl. The reconstituted beads (20 μl) were then added to the alkylated sample followed by ethanol to bring the final concentration to 75% (v/v), with mixing for 20 min. The beads were subsequently immobilized on a magnetic rack. The supernatant was discarded, and the pellet was rinsed with 80% ethanol and 100% acetonitrile (ACN). The beads were then resuspended in 40 μl of 50 mM Tris-HCl pH 8.0 with 1 µg trypsin/Lys-C Mix (Promega, Madison, WI, USA) and digested by gentle agitation at 37 °C overnight. The digested sample was acidified with 150 μl of 0.1% Trifluoroacetic Acid (TFA) and then desalted using GL-Tip SDB (GL Sciences Inc., Tokyo, Japan) according to the manufacturer's instructions, followed by drying with a centrifugal evaporator. The dried peptides were redissolved in 3% ACN and 0.1% formic acid and transferred to a hydrophilic-coated low-adsorption vial (ProteoSave vial; AMR Inc., Tokyo, Japan).
Proteome analyses
Peptides were directly injected onto a 75 μm × 40 cm PicoFrit emitter (New Objective, Woburn, MA, USA) packed in-house with C18 core-shell particles (CAPCELL CORE MP 2.7 μm, 160 Å material; Osaka Soda Co., Ltd., Osaka, Japan) at 60 °C and then separated with a 120-min gradient at 100 nl/min using an UltiMate 3000 RSLCnano LC system (Thermo Fisher Scientific, Waltham, MA, USA). Peptides eluting from the column were analyzed on a Q Exactive HF-X (Thermo Fisher Scientific) for overlapping window data-independent acquisition (DIA)-MS(13, 31). MS1 spectra were collected in the range of 495-745 m/z at 30,000 resolution to set an automatic gain control target of 3×106 (unit) and maximum injection time of 55 (unit). MS2 spectra were collected at >200 m/z at 45,000 resolution to set an automatic gain control target of 3×106 (unit), maximum injection time of “auto,” and stepped normalized collision energy of 22%, 26%, and 30%. The isolation width for MS2 was set to 4 m/z, and overlapping window patterns of 500-740 m/z were used for window placements optimized by Skyline.
The MS files were searched against a human spectral library using Scaffold DIA (Proteome Software, Inc., Portland, OR, USA). The human spectral library was generated from the human protein sequence database (UniProt id UP000005640, reviewed, canonical) established by Prosit. The Scaffold DIA search parameters were as follows: experimental data search enzyme, trypsin; maximum missed cleavage sites, 1; precursor mass tolerance, 8 ppm; fragment mass tolerance, 8 ppm; static modification, cysteine carbamidomethylation. The protein identification threshold was set at <1% for both peptide and protein false discovery rates. The peptide quantification was calculated by the EncyclopeDIA algorithm(32) in Scaffold DIA. For each peptide, the four highest-quality fragment ions were selected for quantitation. The protein quantitative value was estimated from the summed peptide quantitative values. Total quantitative values were calculated by normalizing the protein quantitative values between samples.
Enzyme-linked immunosorbent assays (ELISAs)
Serum IL-6 and IL-18 levels were measured using a commercial ELISA (IL-18: MBL, Nagoya, Japan; IL-6: R&D Systems, Inc., Minneapolis, MN, USA), according to the manufacturer’s instructions.
The statistical and bioinformatics analysis of the data
We performed statistical analyses and generated hierarchical clustering and heatmaps using the JMP Pro13 software program (SAS Institute Inc., Cary, NC, USA) and Qlucore omics explorer software program (Qlucore AB, Lund, Sweden).
We used proteins with at least three peptide counts that were detected in at least one subject for the subsequent analysis, with the exclusion of proteins with missing values in some subjects. After conducting log2 transformation for the quantified continuous values, an independent sample two-sided t-test or an analysis of variance (ANOVA) was used. P values of <0.05, with a confidence interval of 95%, were considered to indicate statistical significance.
We also used a principal component analysis in our study to reduce the identified proteins into a few dimensions that reduce a large amount of the variability of the original values. The first principal component (PC1) explains the highest amount of variability of the original data. Furthermore, correlations of the identified proteins were investigated using Spearman's correlation coefficient method.
An analysis of functional Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were performed with the online tool Search Tool for the Retrieval of Interacting Genes (STRING v11; http://string-db.org/) (33). An association network of the identified proteins was also established using the STRING database with the highest confidence interaction score. The proteins that showed elevated serum values were referred to the Human Protein Atlas (HPA; https://www.proteinatlas.org/) (34).