Comparative proteomics analysis of plasma protein in patients with neuropsychiatric systemic lupus erythematosus.

BACKGROUND: The aim of this study was to evaluate serum biomarkers of systemic lupus erythematous (SLE) patients, with and without neuropsychiatric (NP) manifestation by high-resolution proteomic analysis.
METHODS: SLE patients with NP (NPSLE, n=9), SLE patients without NP (non-NPSLE, n=9) and healthy controls (HC, n=9) were prospectively enrolled in this study, and their plasma samples were collected and pooled into 3 NPSLE, 3 non-NPSLE and 3 HC samples for discovery profile. The TMT-LC-MS/MS-based proteomics approach was used to identify the differential proteome among the three matched groups, and the data were analyzed by bioinformatics tools, including Gene Ontology (GO) categories, Kyoto Encyclopedia of Genes and Genomes enrichment analysis, to explore canonical pathways and networks involved in the pathogenesis of NPSLE. To validation of differentially expressed proteomics results, four proteins were measured by ELISA.
RESULTS: There were altogether 223 differentially expressed proteins in NPSLE groups compared with healthy controls (HC), of which 96 proteins increased while 127 proteins decreased. Compared with non-NPSLE, there were only 49 differentially expressed proteins in NPSLE groups, of which 37 proteins increased while 12 proteins decreased. The significantly changed pathway that those proteins are involved in was complement and coagulation cascades in NPSLE group compared with health controls. However, we didn't find significantly changed pathway between NPSLE group and non-NPSLE group. Five proteins were found significantly changed in all group-comparisons with consistent tendencies using Venn analysis, including Vitamin D binding protein (VDBP), C-reactive protein (CRP), KRT16, IGHV4-4 and CTRP3. Four proteins including CTRP3, VDBP, PAPPA and TRYP2 were selected to estimate the validity of the proteomics approach by ELISA. The expression levels of CTRP3 and TRYP2 were significantly changed in NPSLE patients compared with either HC or non-NPSLE patients.
CONCLUSIONS: Our research has successfully established serum protein profiles of NPSLE and non-NPSLE patients through TMT technology and screened out five proteins significantly changed in group-comparisons with consistent tendencies. The pathway of complement and coagulation cascades may participate in pathogenesis of NPSLE and non-NPSLE. 2020 Annals of Translational Medicine. All rights reserved.

Introduction

Systemic lupus erythematosus (SLE) is a prototypic and complex inflammatory disorder featured by a variety of autoantibodies and multisystem involvement (1). Previous studies have shown that neurologic and neuropsychiatric (NP) symptoms occur in up to 75% of SLE patients, making NPSLE the most common and particularly severe form of lupus (2). The symptoms of NPSLE include mood disorders, psychosis, confusion, headache and cognitive dysfunction, which significantly degraded the quality of life and affect the survival of NPSLE patients (3,4). It was reported that NP events can adversely affect health related quality of life of SLE patients, regardless of disease activity, organ involvement and drug exposure (5). However, the mechanisms underlying NPSLE pathogenesis remain limited, so no single test combining high specificity with high sensitivity is sufficient for diagnosis of NPSLE (2). Therefore, there is an urgent need to evaluate new biomarkers of NPSLE for early diagnosis, prognosis and treatment.

As a source of new biomarkers for lots of diseases, proteomics technology has increasingly maintained attention over the past two decades (6,7). To date, tandem mass tags (TMT) based quantitative proteomics technology has been widely applied in the protein biomarker research and protein alterations quantification in many autoimmune diseases and cancers (8,9). Among the previous proteomic studies about autoimmune diseases, SLE has received some degree of attention (10-12). However, no TMT technology has been reported involved in plasma related to patients with SLE or NPSLE.

Therefore, the current study quantitatively analyzed the proteomics profiling of plasma in NPSLE, active SLE without central nervous system involvement (non-NPSLE) and healthy controls (HC). The goal of this research was to discover the protein changes possibly related with the disease progression, or a biomarker to prompt diagnosis or predict prognosis of NPSLE. So, a TMT-liquid chromatography-triple quadrupole mass spectrometry (LC–MS/MS)-based proteomics method was used to define the differential plasma proteome across plasma samples from NPSLE patients to investigate possible molecular mechanisms involved in NPSLE and pave the way for the novel potential therapeutic targets.

Methods

Patients and healthy controls

In total, eighteen patients with SLE and nine healthy controls were invited to participate in our research. Lupus patients were divided into two groups, 9 patients with NPSLE and 9 non-NPSLE patients. Detailed information of these patient samples is shown in the . All patients were diagnosed according to the American College of Rheumatology (ACR) revised SLE criteria (13). Disease activity of these patients was measured using the Systemic Lupus Erythematosus Disease Activity Index (SLEDAI) (14). NPSLE patients were fulfilled with the ACR nomenclature and case definitions for NPSLE (15).

Table S1

Clinical manifestation for each patient at baseline

	Number	Age	Sex	SLEDAI	Duration (year)	Clinical manifestations
Pool 1	1	43	F	16	5	NP, A, LN, C, V,
	2	27	F	19	0.08	NP, F, C, LN
	3	34	F	13	16	NP, C, P
Pool 2	4	59	F	13	0.25	NP, C
	5	34	F	23	6	NP, A,F
	6	31	F	14	9	NP, LN, C
Pool 3	7	34	F	8	10	NP, F
	8	47	F	5	5	NP, A, F, P
	9	29	F	12	7	NP, LN
Pool 4	10	27	F	14	3	A, LN
	11	30	F	14	0.05	F, LN
	12	37	F	17	1	F, LN
Pool 5	13	25	F	13	3	LN, P
	14	43	F	12	10	C, LN
	15	37	F	10	6	C, LN
Pool 6	16	22	F	14	14	F, LN
	17	29	F	12	5	A, C, F, LN, P
	18	25	F	10	2	C, F

A, Arthritis; C, Cytopenia; F, Febrile; LN, Lupus nephritis; NP, neuropsychiatric; P, Polyserositis; V, Vasculitis.

Altogether 27 plasma samples were collected. 2 mL of whole blood was collected separately from study subjects and K2-EDTA was used as an anticoagulant. All samples were then centrifuged at room temperature at 1,300 ×g for 10 min. The resulting plasma was collected and then kept at −80 °C before used. Within each group, every three plasmas were mixed into single samples in order to obtain three pooled samples.

Validation was performed using a second larger and independent sample set. Totally 23 NPSLE, 26 non-NPSLE and 22 age-matched healthy controls were included. Demographic features of these patient samples are shown in the . The sample cohorts were matched for age and gender. This research was approved by the Ethics Committee at The Affiliated Drum Tower Hospital of Nanjing University Medical School (ID: SC201700201) and was undertaken according to the guidelines of the Declaration of Helsinki. Written informed consent was obtained from the patient for publication of this study and any accompanying images.

Table S2

Demographic features of SLE patients

Characteristics	NPSLE (n=23)	Non-NPSLE (n=26)	HC (n=22)
Demographic
Age (years, mean ± SD)	36.8±13.9	32.3±10.8	36.2±9.7
Sex (F/M)	22/1	26/0	21/1
Disease duration (years)	5.4±6.0	7.3±5.9
SLEDAI (mean ± SD)	11.6±7.0	11.2±3.4
Clinical chart review (%)
Arthritis	52.8	69.2	–
Skin rash	39.1	73.1	–
Photosensitivity	13.0	30.4	–
Oral ulcers	13.0	7.7	–
Proteinuria	34.8	73.1	–
Cytopenia	39.1	34.6	–
Febrile	78.3	53.8	–
Current medication (%)
Glucocorticoid	100	100	–
Hydroxychloroquine	78.3	76.9	–
Cyclophosphamide	30.4	34.6	–
Mycophenolate mofetil	26.1	19.3	–
Other immunosuppressant	11.5	42.3	–

SLE, systemic lupus erythematosus.

Plasma sample preparation and SDS-PAGE

High-abundance proteins was eliminated using the Agilent Multiple Affinity Removal System (Agilent, USA) in this study, according to the manufacturer’s suggestions. The depleted plasma samples were washed three times with 0.5 mL of 25 mM Ammonium bicarbonate, concentrated to 100 µL and rapidly transferred to 1.5 mL polypropylene centrifuge tubes. Each tube was added with 0.1 mL SDT protein lysis buffer and then heated with 100 °C metal bath for 3 min, following by cooling with ultrasonic in an ice bath (50 W for 2 s, interval 8 s, total sonication time 5 min). The samples were heated with 100 °C metal bath for 3 min again. The final protein concentration was measured using a Bradford assay kit (BioRad, Hercules, CA). Subsequently, samples (20 µg) were taken for SDS-PAGE according to the quantitative data. Depletion of each set of plasma samples was completed in the same day so that all the samples were processed consistently in the same set. The protein extracts were stored at –80 °C.

Protein digestion

The protein samples were digested in standard sample buffer by FASP procedure (16). The proteins mixture was reduced for 5 min at 100 °C with 100 mM DTT. After cooling to room temperature, 200 µL of UA buffer (pH 8.5) containing 150 mM Tris-HCl and 8M Urea was added and mixed well. All samples were transferred to ultrafiltration with a 30 kDa cut-off membrane filter (Sartorius, Gottingen, Germany) and centrifuged for 30 min at 14,000 ×g, and then discarded the filtrate. This step was repeated for 3 times. Subsequently, the proteins were alkylated with 100 µL IAA (50mM IAA in UA), shaked for 1 min at 600 rpm and incubated in the dark for 30 min at 300 rpm. Then, each tube was centrifuged for 30 min at 14,000 ×g. Add 100 µL UA buffer, centrifuge for 30 min at 14,000 ×g. This step was repeated for 3 times. 100 µL dissolution buffer (100 mM/L) (Applied Biosystems, USA) was added, and centrifuged for 30 min at 14,000 ×g at room temperature (RT). This step was repeated for 3 times. Finally, the filtrate was discarded and the proteins were digested with trypsin on a thermostat mixer at 300 rpm for 18 h at 37 °C.The tryptic peptides were centrifuged (14,000 ×g, 30 min, RT), replaced with new tube, and then added with 40 µL of 25 mM DS buffer. The tubes were centrifuged (14,000 ×g, 30 min, RT) and the filtrate was collected. Finally, the tryptic peptides were quantified by OD280 (16).

TMT Labeling and high-pH reversed-phased chromatography separation

Three pooled NPSLE patient sera samples, 3 pooled non-NPSLE patient sera samples and 3 pooled healthy control sera samples, with each pool comprised of sera aliquots from 3 individuals (total of 9 healthy control and 18 SLE patients), were included in this global proteomics discovery study. Digested peptides were labeled with TMT reagents (Thermo Scientific) following procedures recommended by the manufacturer. All labeled peptides were pooled together. 1,100 Series HPLC Value System (Agilent) equipped with a Gemini-NX (Phenomemex, 00F-4453-E0) column (4.6 × 150 mm, 3 µm, 110 Å) was used for high-pH reverse-phase HPLC to fractionate peptide samples. Peptides were separated into 40 fractions using a gradient of 4.5% to 90% ACN in 10mMammonium bicarbonate (pH 10) over 130 min. The peptides were then combined into 18 fractions and dried by vacuum centrifugation. Samples were storaged at −80 °C freezer for further analysis.

LC-MS/MS analysis

Easy-nLC nanoflow HPLC system connected to Orbitrap-Fusion (Thermo Fisher Scientific, San Jose, CA, USA) were used to analyze the TMT-labeled samples. Peptides were reconstituted in 0.1% v/v trifluoroacetic acid. Each sample (1 µg) was loaded onto Thermo Scientific EASY column (two columns) using an autosampler at a flow rate of 200 nL/min. Peptides were accomplished using a segmented gradient from 5% to 28% Solvent B (0.1% formic acid in 100% ACN) for 40 min, followed by 28–90% Solvent B for 2 min and then 90% Solvent B for 18 min. The column was re-equilibrated to its initial highly aqueous solvent composition before each analysis. LC-MS/MS analysis was performed on an Obitrap Fusion mass spectrometer (Thermo Scientific). The MS scan and MS/MS scan were acquired at a resolution of 1.2×105 and 5×104 at 200 m/z, respectively. The maximum ion injection times were maintained at 105 ms for the MS/MS scans, and automatic gain control target values for Master scan modes was set to 4×105 and 1×105.

Database search and protein quantification

Using the Proteome Discoverer 2.1 software, acquired MS/MS spectra were analyzed for protein identification and quantification. Searching for the fragmentation spectra was carried out by the MASCOT search engine embedded in Proteome Discoverer against the Uniprot human database (March 8, 2017, 156914 sequences). A false discovery rate (FDR) of no more than 1% was used for data filtering.

Using Proteome Discoverer (version 2.1), the analysis of the proteins based on the ratios of TMT reporter ions from all unique peptides representing each protein was performed using Proteome Discoverer (version 2.1) in the samples. The peak intensities of the TMT reporter ions released in each of the MS/MS spectra were used, and the IS sample was used as a reference in analyzing for the TMT ratios of the reporter ions. Finally, the median of the unique peptides of the protein by calculating relative protein quantifications were first normalized.

Bioinformatics analysis

The Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) database were applied to explore the potential roles of differentially expressed proteomics by DAVID 6.8 (https://david.ncifcrf.gov/). The proteins were classified for the next three categories: molecular function, cellular component, and biological process.

ELISA experiments

Validation of differentially expressed proteomics results was performed using independent sample from the same cohort in proteomic analysis. Comprehensive considering proteomics profiles and function of proteins, we chose to verify four proteins measured by ELISA, respectively, including C1q/TNF-related protein-3 (CTRP3), vitamin D binding protein (VDBP), pregnancy-associated plasma protein-A (PAPPA) and serine protease 2 (TRYP2). Serum concentrations of CTRP3, VDBP, PAPPA and TRYP2 were all measured using commercially available ELISA kits. CTRP3 and VDBP ELISA kits were from R&D Systems. PAPPA and TRYP2 ELISA kits were from RayBiotech. We detected the serum concentrations according to the manufacturer’s instructions.

Statistics

All data were presented as mean ± SD. Differentially expressed proteins were identified by the fold of change (fold >1.5) method and independent t-test. Differences between two groups were determined by unpaired Student’s t-test if the variance was normally distributed. Comparisons among three or more groups were conducted using one-way ANOVA. Randomization and blinding strategy was used whenever possible. The correlations were analyzed by Spearman’s correlation test. Data were analyzed and visualized with SPSS16.0 software or GraphPad Prism 7.0 (GraphPad Software, La Jolla, CA, USA) and a two-tailed P value <0.05 was considered statistically significant.

Results

Protein profile changes in NPSLE patients identified by MS/MS analysis

To determine the possible role of protein profiles in NPSLE pathogenesis, we performed a TMT-labeled quantitative proteomic study to compare the protein profile differences among NPSLE patients, non-NPSLE patients and healthy subjects. Information for the samples is shown in . A total of 726 proteins were identified in our proteomic studies between SLE and healthy controls using scatter plot and hierarchical clustering analyses (). The original raw data were showed in the in total online: http://fp.amegroups.cn/cms/aced7f51ea115b5ba81dd6fdf54db336/atm.2020.04.58-1.xlsx. Hierarchical clustering demonstrated that two groups were clustered based on the protein abundance levels of the samples: proteins of NPSLE and non-NPSLE clustered together, whereas those of HC clustered independently.

Table 1

Demographic features of the 18 SLE patients and 9 healthy controls in the study

	SLE patients		Healthy controls
	NPSLE	Active without NP
Age (years, mean ± SD)	36.6±11.1	30.6±7.0	36.7±7.1
Sex (F/M)	9/0	9/0	9/0
Disease duration (years)	6.5±4.9	4.9±4.5
SLEDAI (mean ± SD)	13.7±5.4	12.9±2.2

SLE, systemic lupus erythematosus.

Figure 1

Heat map showing hierarchical clustering of proteins with expression changes greater than two-fold and P value <0.01. NPSLE: SLE1–3; non-NPSLE: SLE3–6; Healthy control group: HC1–3. Red and blue colors represent up- and downregulated proteins, respectively. NPSLE, neuropsychiatric systemic lupus erythematosus; SLE, systemic lupus erythematosus; HC, healthy control.

Further, the greatest number of proteins with different abundances occurred when analyzing NPSLE and HC samples. There were fewest proteins with significantly different abundances in NPSLE/non-NPSLE comparisons. More specifically, compared with HC group, 96 proteins were increased and 127 were decreased (cutoff ratio of >1.50 and <0.667) among proteins identified () in the NPSLE group. The details of proteins differently expressed between NPSLE and non-NPSLE groups were showed in the in total online: http://fp.amegroups.cn/cms/d75d5c8a7207877825537895b750a77a/atm.2020.04.58-2.xlsx. The volcano plot () depicted the distribution of proteins according to statistical significance (P values) and NPSLE/HC abundance ratio. Compared with non-NPSLE group, 37 proteins were increased and 12 were decreased (cutoff ratio of >1.50 and <0.667) in NPSLE/ non-NPSLE group (), so did the volcano plot show (). The details of proteins differently expressed between NPSLE and HC groups were showed in the in total online: http://fp.amegroups.cn/cms/f965d561360731ff332d3f13d6566ecc/atm.2020.04.58-3.xlsx. Similarly, 56 proteins were increased and 78 were decreased (cutoff ratio of >1.50 and <0.667) in non-NPSLE group compared to that of HC group (), so did the volcano plot show (). The details of proteins differently expressed between non-NPSLE and HC groups were showed in the in total online: http://fp.amegroups.cn/cms/6acf489810c38ddc194d3a3802d58dd9/atm.2020.04.58-4.xlsx.

Figure 2

The differentially expressed proteins in NPSLE, non-NPSLE patients and HCs. (A,C,E) Histogram showed the differentially expressed proteins between NPSLE/HC, NPSLE/non-NPSLE and non-NPSLE/HC samples; (B,D,F) volcano plot of differentially expressed proteins. The horizontal dotted lines represent a P value of 0.01 and 0.05, and vertical dotted lines represent 2.0-fold changes up and down. X axes are the fold change values (log2 scaled), and Y axes are the P values (log10 scaled). Red and blue plots represent up- and downregulated genes, respectively. NPSLE, neuropsychiatric systemic lupus erythematosus; HC, healthy control.

Verification of differentially expressed proteins by ELISA

To verify the differentially expressed proteins which were identified by TMT-labeled LC-MS/MS analysis, four aberrantly expressed proteins were selected for further analysis by ELISA in the serums from NPSLE patients, non-NPSLE patients and healthy subjects. First, we analysis the differentially expressed proteins between two different groups by Venn analysis, indicating that only five proteins were got using the following criteria: tendencies of abnormally expressed proteins were consistent between NPSLE-non-NPSLE group and between non-NPSLE-HC group (). The five proteins were showed in , including VDBP, C-reactive protein (CRP), Keratin 16 (KRT16), Immunoglobulin heavy variable 4-4 (IGHV4-4) and CTRP3. The annotative MS spectra for the unique peptides of 5 identified proteins were provided in the . In view of not all proteins have known function, we chose VDBP and CTRP3 to verify. In addition, to further analyze the difference between NPSLE and non-NPSLE, we selected two proteins (PAPPA and TRYP2) that significantly changed between NPSLE/non-NPSLE and NPSLE/HC comparisons, but with no significant difference between non-NPSLE/HC comparisons. As can be seen from , the expression of CTRP3 levels was significantly decreased in patients with NPSLE compared with patients with non-NPSLE (P<0.05) and healthy controls (P<0.001). CTRP3 levels were also significantly decreased in patients with non-NPSLE compared with healthy controls (P<0.001). The similar results could be seen in the expression of VDBP levels in three groups (). The expression of the TRYP2 levels was significantly increased in patients with NPSLE compared with patients with non-NPSLE (P<0.01, ) and healthy controls (P<0.001, ). There was no difference between patients with non-NPSLE and healthy controls (P>0.05, ). PAPPA levels were not dysregulated between these groups (P>0.05, ).

Figure 3

The differentially expressed proteins in all three groups. (A) The differentially expressed proteins were showed by Venn diagram; (B) The specific five differentially expressed proteins significantly changed in all three groups. NPSLE, neuropsychiatric systemic lupus erythematosus; HC, healthy control.

Figure 4

The differentially expressed proteins levels in plasma in of NPSLE patients, non-NPSLE patients and HCs are identified using ELISA. (A) The CTRP3 levels were significantly decreased in NPSLE patients and non-NPSLE patients. n (NPSLE) =20, n (non-NPSLE) =26, n (HC) =22; (B) the VDBP levels were significantly decreased in NPSLE patients and non-NPSLE patients. n (NPSLE) =22, n (non-NPSLE) =26, n (HC) =22; (C) the TRYP2 levels were significantly decreased in NPSLE patients instead of non-NPSLE patients. n (NPSLE) =22, n (non-NPSLE) =26, n (HC) =22; (D) there was a trend that PAPPA levels were increased in NPSLE patients. n (NPSLE) =19, n (non-NPSLE) =23, n (HC) =19. *, P<0.05; **, P<0.01; ***, P<0.001. ns, P>0.05.

Figure S1

The annotative MS spectra for the unique peptides of 5 identified proteins. (A) VDBP; (B) CRP; (C) KRT16; (D) IGHV4-4; (E) CTRP3.

Next, to define potential clinical effects of those proteins in NPSLE, we analyzed the correlations between the differentially expressed proteins levels and the clinical characteristics of NPSLE patients. As showed in , the expression levels of CTRP3 were associated with aminotransferase (ALT), uric acid (UA), glomerular filtration rate (GFR) and NK cells numbers (P<0.05). Whereas, VDBP levels were found to be associated with globin (Glo), triglyceride, serum IgA and IgG levels (P<0.05, ). Furthermore, the expression levels of TRYP2 were associated with UA, CRP, GFR and CD4+ cells numbers (P<0.05, ). The expression levels of PAPPA were associated with platelets (PLT), CRP, serum IgG levels and CD4+ cells numbers (P<0.05, ).

Figure 5

Correlation between differentially expressed proteins and clinical characteristics in NPSLE patients. (A,B,C,D) The expression levels of CTRP3 were associated with aminotransferase (ALT), uric acid (UA), glomerular filtration rate (GFR) and NK cells numbers; (E,F,G,H) the expression levels of VDBP were associated with globin (Glo), triglyceride, serum IgA and IgG levels; (I,J,K,L) the expression levels of TRYP2 were associated with UA, CRP, GFR and CD4+ cells numbers; (M,N,O,P) the expression levels of PAPPA were associated with platelets (PLT), CRP, serum IgG levels and CD4+ cells numbers.

Gene ontology (GO) analysis

The GO data were used to examine the subcellular annotation of totally expressed proteins and the function of the differentially expressed proteins.

Subcellular location annotation of totally proteins

As showed in , the totally expressed proteins were quantified according to the subcellular location annotation. Most proteins were concentrated in the following categories: cytoplasm, membrane, endoplasmic reticulum (37.8%, 34.1% and 8.8%, respectively).

Figure 6

Subcellular location annotation of totally proteins. Most proteins were concentrated in the following categories: cytoplasm, membrane, endoplasmic reticulum.

Evaluation of differentially expressed proteins by GO analysis

To further evaluate the potential biological significance of the 726 proteins identified in three groups, GO enrichment analysis was employed to analyze the proteins. As shown in , the top 10 GO items from the differentially expressed proteins between the different groups were listed. Our data demonstrated that the most significantly enriched cellular components of the differentially expressed proteins in serum of NPSLE and non-NPSLE were extracellular region, extracellular region part and membrane-bounded vesicle (). The most significantly enriched cellular components of NPSLE and healthy controls were extracellular space, extracellular region part and extracellular exosome (). The most significantly enriched cellular components of non-NPSLE and healthy controls were extracellular vesicle, extracellular organelle and extracellular region part ().

Figure 7

Enrichment analysis of GO terms for differentially expressed proteins between NPSLE, non-NPSLE and healthy controls. (A) The top 10 GO analysis which consisted of significant molecular function, cellular component and biological process of differentially expressed proteins between NPSLE and non-NPSLE; (B) the top 10 GO analysis between NPSLE and healthy controls; (C) the top 10 GO analysis between non-NPSLE and healthy controls.

As for biological processes, the most significantly enriched item of NPSLE and non-NPSLE were regulation of body fluid levels, establishment of skin barrier and regulation of water loss via skin (). The most significantly enriched biological processes of NPSLE and healthy controls were acute inflammatory response, defense response and protein activation cascade (). The most significantly enriched biological processes of non-NPSLE and healthy controls were acute inflammatory response, inflammatory response and defense response ().

As for molecular functions, the most significantly enriched item of NPSLE and non-NPSLE were endopeptidase activity, peptidase activity and serine type endopeptidase activity (). The most significantly enriched molecular functions of NPSLE and healthy controls were peptidase inhibitor activity, peptidase regulator activity and glycosaminoglycan binding (). The most significantly enriched molecular functions of non-NPSLE and healthy controls were serine-type endopeptidase activity and calcium ion binding and serine-type peptidase activity (). These together reveal the roles of significantly dysregulated proteins in NPSLE.

KEGG pathway enrichment analysis

The KEGG pathway enrichment analysis was performed to identify the statistically significant differently expressed proteins () (P<0.05). Proteins in NPSLE and controls were mainly involved in complement and coagulation cascades. The same finding was in non-NPSLE and controls. There was no significantly changed pathway in NPSLE and non-NPSLE.

Figure 8

KEGG pathway enrichment analysis of NPSLE, non-NPSLE and healthy controls. (A) The bubble chart showed the top 10 pathways enriched in NPSLE and healthy controls. (B) The bubble chart showed the top 10 pathways enriched in non-NPSLE and healthy controls.

Discussion

Proteins are the direct function executors of myriad life activities (17). Compared to genomics, proteomics is highly dynamic and regulated by multiple factors, including protein expression/degradation, subcellular localization, post-translational modifications (PTMs) and proteins-proteins interactions and other biomolecules (18). Therefore, fast, high sensitivity and high throughput proteomics technology is of great significance. In recent years, proteomics technology combined with mass spectrometry has achieved extraordinary improvements and been widely used in different clinically orientated studies (7). Chemical labeling with isobaric TMT, including TMT and iTRAQ, has been the most popular methods among the quantitative proteomic technology (19,20). Wang et al. established protein expression profiles of PBMC from patients with SLE and healthy subjects with iTRAQ-based technology and identified 67 unique proteins (21). Similarly, Wang et al. also found that STRAP, an important inhibitor of TGF-ß signaling, decrease in active SLE patients compared with healthy controls using iTRAQ-MS (11). In another study, renal tissue from patients with lupus nephritis and healthy subjects was analyzed by iTRAQ-MS and four proteins were detected, including up-regulated hnRNP-, Annexins and down-regulated ASS, aldolase (12). Therefore, we identified differently expressed proteins using TMT-labeled LC-MS/MS analysis in serum from NPSLE, non-NPSLE patients and healthy controls in this study.

Neurological and psychiatric features are a serious complication of SLE, which can result in severe neurodegenerative changes and threaten life (22). It is reported that infection, lupus nephritis, hematological abnormality, NPSLE, and interstitial pneumonia are top causes of deaths in China (23). The two recognized underlying mechanisms of the immunopathogenesis of NPSLE were ischemic and inflammatory mechanisms, including blood-brain barrier dysfunction and production of pro-inflammatory cytokines which mediated by autoantibody (5,24). However, NPSLE pathogenesis is incompletely understood. It was reported that many autoantibodies have been linked to NPSLE, including anti-NR2 and anti-ribosomal P antibodies. However, there is no enough evidence to confirm involvement in mechanisms of pathogenesis (25). S100 β were increased in serum of adults and children with NPSLE, suggesting S100 β as a biomarker (26). Nevertheless, these biomarkers lack disease specificity of NPSLE. Nevertheless, these biomarkers lack disease specificity of NPSLE. There is a desperate need for discovering reliable NPSLE biomarkers that can be used for diagnosis, disease monitoring and better prediction of response to therapy. Therefore, we investigated the abnormally expressed proteins between NPSLE patients and non-NPSLE patients. Totally 726 proteins involving different cellular localizations and biological functions were identified in our proteomic studies. Further analysis, we found higher numbers of differentially expressed proteins in NPSLE/HC comparisons than those in non-NPSLE/HC comparisons, indicating that the inflammatory damage might be progressive in NPSLE.

Then, we found only five proteins significantly changed in all three groups-comparisons by Venn analysis, including VDBP, CRP, KRT16, IGHV4-4 and CTRP3. CRP, a biomarker of inflammation, was increased in patients with SLE both in remission and active periods when compared with control subjects (27). However, the levels of CRP were not always related to SLE disease activity. The significance of CRP remains controversial and is not one of the recommended tools in SLE disease activity follow-up (28).The study on KRT16 and IGHV4-4 concerning lupus is rare. Above all, we did not verify the three proteins. VDBP (also known as Gc-globulin) is a multifunctional protein and the major plasma carrier protein of vitamin D and its metabolites (29). In our study, we found that serum VDBP levels were significantly decreased in SLE patients compared with healthy controls, and they were lower in patients with NPSLE than with non-NPSLE. It was reported that vitamin D deficiency may contribute towards the pathogenesis of SLE (30). Our result, along with other reported findings, suggests decreased levels of vitamin D, decreased levels of VDBP, or both may participate in the pathogenesis and process of SLE. However, studies on VDBP in NPSLE were rare. Our results showed that the levels of VDBP in NPSLE were significantly decreased compared with non-NPSLE, which indicate that VDBP may participate in the pathogenesis of NPSLE. Nevertheless, the detailed role it plays in NPSLE still needs further study.

CTRP3, a member of C1q TNF Related Protein (CTRP) family, plays important roles in modulating glucose and lipid metabolism (31). Recent studies have demonstrated that CTRP3 has anti-inflammatory effects in attenuating LPS-induced systemic inflammation (32). Moreover, the deficiency of CTRP3 markedly exacerbated inflammation in a mouse model of rheumatoid arthritis (33). In our study, we found that the expression of CTRP3 were extremely decreased in patients with NPSLE when compared with patients without NPSLE and healthy subjects. As cardiovascular disease is one of the leading causes of morbidity and mortality in patients with SLE and dyslipidemia in systemic lupus erythematosus has gained extensive attention (34), the role of CTRP3 that plays in SLE is worth further studying.

In addition, we also selected two proteins, PAPPA and TRYP2, which significantly changed between NPSLE/non-NPSLE and NPSLE/HC comparisons, but with no significant difference between non-NPSLE/HC comparisons. PAPPA, a member of the metzinc in metalloproteinase superfamily, was reported that exerts a proatherogenic effect in atherosclerosis (35). Besides, it was report that PAPPA was associated with the outcome of ischemic cerebrovascular disease (36). However, we only found that there was a tendency to increase of PAPPA in NPSLE. What’s more, we also detected the expression of PAPPA in cerebrospinal fluid in SLE patients and the results were consistent with in the serums. As we know, one underlying mechanisms of NPSLE progression was ischemic. We believed that more studies should be conducted to explore the mechanism by which PAPPA might be involved in the pathogenesis of NPSLE. There are few studies on trypsin 2 (TRYP2) in NPSLE, mainly focusing on tumors, and more research is needed to discover the role of TRYP2 in NPSLE.

In KEGG analysis, the pathway of complement and coagulation cascades has been emphasized between not only NPSLE-HC but also non-NPSLE-HC. The role of the complement pathway in lupus has been widely reported, and our research has shown that the complement pathway plays a critically important role not only in non-NPSLE patients, but also in NPSLE patients (37-39). Unfortunately, there was no significantly changed pathway in NPSLE and non-NPSLE in our study.

In conclusion, our study has successfully increased serum protein profiles of NPSLE, non-NPSLE patients and healthy controls through TMT technology and screened out five proteins significantly changed in NP-NPSLE patients. The pathway of complement and coagulation cascades may participate in pathogenesis of NPSLE.

The article’s supplementary files as