- 1Department of Nephrology, The First Affiliated Hospital of Jinan University, Guangzhou, China
- 2Guangdong Provincial Engineering Research Center of Autoimmune Disease Precision Medicine, The First Affiliated Hospital of Southern University of Science and Technology, Shenzhen People’s Hospital, The Second Clinical Medical College of Jinan University, Shenzhen, China
- 3Fifth Department of Medicine (Nephrology/Endocrinology/Rheumatology), University Medical Centre Mannheim, University of Heidelberg, Heidelberg, Germany
Background: Ankylosing spondylitis (AS) is a chronic inflammatory disease with serious consequences and a high rate of morbidity and mortality, In our previous work, we reveal the key features of proteins in new-onset ankylosing spondylitis patients.
Material and Methods: Ankylosing spondylitis (AS) is a chronic inflammatory condition that affects the spine, and inflammation plays an essential role in AS pathogenesis. The inflammatory process in AS, however, is still poorly understood due to its intricacy. Systematic proteomic and phosphorylation analyses of peripheral blood mononuclear cells (PBMCs) were used to investigate potential pathways involved in AS pathogenesis.
Results: Liquid chromatography-tandem mass spectrometry (LC–MS/MS) analysis was performed and discovered 782 differentially expressed proteins (DEPs) and 122 differentially phosphorylated proteins (DPPs) between 9 new-onset AS patients and 9 healthy controls. The DEPs were further verified using parallel reaction monitoring (PRM) analysis. PRM analysis verified that 3 proteins (HSP90AB1, HSP90AA1 and HSPA8) in the antigen processing and presentation pathway, 6 proteins (including ITPR1, MYLK and STIM1) in the platelet activation pathway and 10 proteins (including MYL12A, MYL9 and ROCK2) in the leukocyte transendothelial migration pathway were highly expressed in the PBMCs of AS patients.
Conclusion: The key proteins involved in antigen processing and presentation, platelet activation and leukocyte transendothelial migration revealed abnormal immune regulation in patients with new-onset AS. These proteins might be used as candidate markers for AS diagnosis and new therapeutic targets, as well as elucidating the pathophysiology of AS.
Introduction
Ankylosing spondylitis (AS) is a chronic inflammatory disease that can cause significant complications and is associated with high morbidity and mortality (1–3). The pathologic mechanism of AS remains largely unclear but is thought to be immune-mediated and has a strong genetic association with HLA-B27, a class I human leukocyte antigen (4, 5). Studies on AS have revealed that there is a tendency for endothelial dysfunction (6), and inflammation plays an important role in endothelial dysfunction (7). However, the complexity of the inflammatory process in AS is still poorly understood. Platelets have a major role in the interplay among atherogenesis, thrombosis, and inflammation (8). The interaction of platelets with leucocytes, their ability to modify leucocyte functions, their adherence to the endothelium, and their ability to stimulate endothelial cells are all significant functions of platelets (9). There is evidence suggesting that platelet activation may be an indicator of AS progression (10).
Proteomics is a new technology that identifies protein networks using a high-throughput discovery approach with high sensitivity and accuracy (11, 12). Phosphorylation of proteins is one of the most fundamental, widespread, and critical methods for regulating protein activation and function (13). Abnormal phosphorylation is associated with many diseases, and it is an indicator of abnormal status of organisms (14). However, There hasn’t been a systematic proteome and phosphorylation analysis of a particular cell type in AS yet. Professional antigen-presenting cells, such as macrophages and dendritic cells, are derived from monocytes, and the role of monocytes in AS is still elusive (15).
In the present study, we used proteomic and phosphomic techniques coupled with systematic bioinformatics analysis to show that differences in protein expression and phosphorylation in AS peripheral blood mononuclear cells (PBMCs) implicate potential pathways involved in disease pathogenesis. We used PRM to evaluate target proteins to obtain a more global picture of protein expression in AS patients.
Patients and Methods
Patients
The clinical information of the participants has been reported in our previous study (16). This observational study was conducted at the Department of Nephrology, The Affiliated Hospital of Jinan University (Guangzhou, China). All procedures were approved by the local ethical committee of The Affiliated Hospital of Jinan University and followed the tenets of the Declaration of Helsinki. Nine new-onset HLA-B27-positive AS patients who fulfilled the 1984 revised New York criteria (17) and nine age-matched healthy controls were recruited [median age (interquartile range (IQR), 27 (30, 34) vs. 34 (29.5, 36), respectively]. Patients with autoimmune diseases and infectious diseases were excluded from this study. Cubital vein blood samples were collected after obtaining informed consent from the participants. An automated blood cell counter was used to determine the white blood cell count, lymphocyte, monocyte, and neutrophil differentials throughout a routine blood examination. (Beckman Coulter, Ireland).
Isolation of PBMCs
Peripheral venous blood (10 ml) was collected from each individual in ethylenediaminetetraacetic acid (EDTA) anticoagulant tubes and diluted with an equal volume of phosphate-buffered saline (PBS). Diluted blood was gently overlaid on 10 ml of Ficoll-Paque Plus solution (Amersham Biosciences, Uppsala, Sweden) and centrifuged at 500 g for 20 min at room temperature. PBMCs were collected from the interphase and washed twice with PBS. The supernatant was discarded, and the PBMCs were then resuspended at a final concentration of 1 × 106 cells/ml in RPMI 1640 medium (Gibco) containing 10% heat-inactivated fetal calf serum, 100 U/mL penicillin and 100 µg/mL streptomycin. The viability of the cells was consistently greater than 95% as measured by trypan blue solution (Sigma, Germany).
Protein Extraction and Trypsin Digestion
The protein extraction and digestion procedures were performed according to a previously published method (18). We removed samples from -80°C cryogenic storage and then sonicated them three times on ice using a high-intensity ultrasonic processor (JY92-IIN, Scientz, China) in a fourfold volume of lysis buffer (containing 8 M urea, HALT protease and phosphatase inhibitor cocktail). The cell fragments were removed by centrifugation for 10 min at 12 000 g at 4°C. Finally, the supernatant was collected and transferred to a new centrifuge tube. The protein concentration was determined using a BCA kit (Beyotime, China).
For digestion, the protein solution was reduced with 5 mM dithiothreitol (DTT) for 30 min at 56°C and alkylated with 11 mM iodoacetamide (IAA) for 15 min at room temperature in darkness. Protein samples were diluted to less than 2 M by adding 100 mM triethyl ammonium bicarbonate (TEAB), and trypsin (Promega, USA) was added at a 1:50 trypsin-to-protein mass ratio for the first digestion overnight followed by a second 1:100 trypsin-to-protein digestion for 4 h.
Modification Enrichment and Affinity Enrichment
For the enrichment of the modified peptides, peptide mixtures were incubated with IMAC microsphere suspensions with vibration in loading buffer (50% acetonitrile/6% trifluoroacetic acid). The IMAC microspheres with enriched phosphopeptides were collected by centrifugation, and the supernatant was removed. The IMAC microspheres were then sequentially washed with 50% acetonitrile/6% trifluoroacetic acid and 30% acetonitrile/0.1% trifluoroacetic acid. The enriched phosphopeptides were eluted from the IMAC microspheres with elution buffer containing 10% NH4OH. Finally, the eluted phosphopeptides were combined and vacuum dried. For LC–MS/MS analysis, the resulting peptides were desalted with C18 ZipTips (Millipore) according to the manufacturer’s instructions.
LC–MS/MS Analysis
We performed LC–MS/MS analysis on an EASY-nLC 1000 UPLC system (Thermo Fisher Scientific). Peptides were subjected to a NSI source followed by tandem mass spectrometry (MS/MS) in a Q Exactive Plus (Thermo Fisher Scientific) coupled online to the UPLC. Intact peptides were detected in the Orbitrap at a resolution of 70,000, and peptides were selected for MS/MS using a NCE setting of 28 and 30. Ion fragments were detected in the Orbitrap at a resolution of 17,500. A data-dependent procedure that alternated between one MS scan followed by 20 MS/MS scans was applied for the top 20 precursor ions above a threshold ion count of 5.0E3 in the MS survey scan with a 15.0 s dynamic exclusion. The applied electrospray voltage was 2.0 kV. Automatic gain control (AGC) was used to prevent overfilling of the orbitrap, and 5E4 ions were accumulated for the generation of MS/MS spectra. For MS scans, the m/z scan range was 350–1,800, and the fixed first mass was set as 100 m/z.
Database Search and Bioinformatics Analysis
The MS/MS data were retrieved by the MaxQuant search engine (v1.6.6.0). A human database was searched (Swiss-Prot, downloaded on May 13, 2019), retrieving 20422 protein sequences. The decoy database anti-library was used to reduce the false-positive rate (FDR). FDR was adjusted to < 1%, and the minimum score for modified peptides was set > 40. The ratio of the average protein intensity between AS patient samples and healthy control samples was considered as the relative quantitative value. Proteins with a fold-change ≥1.50 or ≤0.67 were considered expression/phosphorylation significant. Based on the protein sequence alignment method, the protein domain functions were defined by InterProScan (http://www.ebi.ac.uk/interpro/). Functional annotation enrichments of differentially expressed proteins (DEPs) and differentially phosphorylated proteins (DPPs) were performed by Gene Ontology (GO) annotation analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Significant enrichment was identified as p < 0.05 in Fisher’s exact test and q < 0.05 in Benjamini-Hochberg’s procedure. We conducted hierarchical clustering analysis for the DEPs and DPPs using the “heatmap.2” function form and “gplots” R package.
Validation of Target Proteins Using Parallel Reaction Monitoring
The differential abundance of selected proteins was verified using PRM analysis. Skyline (v.3.6) software was used to process the resulting MS data. The parameters of the peptide segment were determined as follows: the enzyme was set at trypsin [KR/P]; the maximum number of missed cleavages was set as 2; and the length of the peptide was set as 8−25. The variable modification was set as carbamidomethyl on Cys and oxidation on Met, and the max variable modifications were set as 3. The transition parameters were as follows: the precursor charges were set as 2 and 3; the ion charges were set as 1 and 2; and the ion type was set as b, y and p. Fragment ion selection started from the third to the last, and the ion match tolerance was set as 0.02 Da.
Results
The clinical and laboratory results between AS patients and healthy controls have been reported in our previous study (16).
Identification of Protein Expression and Phosphorylation Changes
A total of 2240 proteins were identified, and 782 proteins were significantly changed compared to the healthy controls. Of these altered proteins, 646 proteins were upregulated (≥1.5-fold), and 136 proteins were downregulated (≤0.67-fold) (Figure 1A). For protein phosphorylation, 210 phosphorylation sites in 122 proteins were significantly changed between AS patients and healthy controls. In addition, 50 phosphorylation sites were upregulated (≥1.5-fold) in 37 proteins, and 125 phosphorylation sites were downregulated (≤0.67-fold) in 85 proteins (Figure 1B).
Figure 1 Protein identification. (A, B) DEPs and DPPs between AS patients and healthy controls. DEPs, differentially expressed proteins; DPPs, differentially Phosphorylated proteins.
Peptide Motifs Associated With Phosphorylation
In total, 1855 sequences, containing the 12 residues, were obtained with six upstream and six downstream residues surrounding each of the phosphorylation sites. Of these sequences, 1761 (94.93%) were centered on a serine residue, and 94 (5.07%) were centered on a threonine (Figure 2A). The serine phosphorylation category included 24 overrepresented motifs as follows: the most frequent motifs were ‘sP’ (291 occurrences) and ‘Rxxs’ (215 occurrences) followed by ‘PxsP’ (136 occurrences), ‘sDxE’ (119 occurrences) and ‘sxE’ (102 occurrences) (Figure 2B). ‘tP’ was the most frequent motif in the threonine phosphorylation category (Figure 2C).
Figure 2 Enrichment analysis of phosphorylated protein motifs (A). Heatmap of the amino acid compositions of the phosphorylation sites showing the frequency of the different types of amino acids around serine (B) and threonine (C). Red indicates significant enrichment of the amino acid near the modification site, and green indicates significant reduction of the amino acid near the modification site.
Annotation and Classification of All Identified DEPs and DPPs
To investigate the functions of the identified and quantified proteins, we annotated functional category distributions according to COG categories. The differentially expressed proteins (DEPs) identified were classified into 25 COG categories. Among these, “Signal transduction mechanisms” (Group T, 107 DEPs) and “Posttranslational modification, protein turnover and chaperones “ (Group O, 105 DEPs) represented the two largest groups followed by “General function prediction only” (Group R, 85 DEPs), “Intracellular trafficking, secretion, and vesicular transport” (Group U, 72 DEPs) and “Cytoskeleton” (Group Z, 67 DEPs) (Figure 3A). Differentially phosphorylated proteins (DPPs) were classified into 16 COG categories, and the most predominant categories were “Signal transduction mechanisms” (Group T, 30 DPPs) and “Cytoskeleton” (Group Z, 19 DPPs) (Figure 3B).
Enrichment Analysis of DEPs and DPPs
In GO terms, 30 different groups of pathways were established for DEPs and DPPs, which were further divided into the following three classifications: cellular component (n=8), molecular function (n=8) and biological process (n=14). For DEPs, the most enriched GO terms were as follows: ‘actin cytoskeleton organization’ (68 proteins) in biological process; ‘secretory granule’ (193 proteins) in cellular component; and ‘glycoprotein binding’ (23 proteins) in molecular function (Figure 4A). For DPPs, the most enriched GO terms were as follows: ‘secretion’ (34 proteins) in biological process; ‘cytoskeleton’ in cellular component (51 proteins); and ‘integrin binding’ (9 proteins) in molecular function (Figure 4B).
Figure 4 GO-based enrichment analysis of DEPs (A) and DPPs (B) in terms of cellular component, molecular function and biological process.
We also performed KEGG pathway enrichment analysis by calculating Fisher’s exact test P values to characterize the predicted roles of highly altered proteins. KEGG pathway-based enrichment analysis revealed 10 pathways for DEPs, including ‘complement and coagulation cascades’ (27 proteins), ‘cGMP-PKG signaling pathway’ (27 proteins) and ‘leukocyte transendothelial migration’ (26 proteins). (Figure 5A). Regarding phosphorylation, 20 pathways were identified for DPPs, and ‘platelet activation’ (15 proteins) was the most significant pathway for enrichment (Figure 5B).
Protein domain enrichment analysis revealed that DEPs were significantly enriched in 25 terms, including small GTP-binding protein domain, thioredoxin-like fold, GroEL-like apical domain and GroEL-like equatorial domain. (Figure 6A). DPPs were significantly enriched in 10 terms, including actin-depolymerizing factor homology domain, SH3 domain and leucine-rich repeat N-terminal domain (Figure 6B).
Protein–Protein Interaction Network Analysis
To better understand the interaction of the DEPs and DPPs, we subsequently constructed a PPI network using STRING software. A total of 782 DEPs and 43 DPPs were mapped to the protein interaction database, presenting a global view of the diverse cellular functions of DEPs and DPPs in AS. As shown in Figure 7, DEPs were involved in leukocyte transendothelial migration, oxidative phosphorylation, gap junctions, bacterial invasion of epithelial cells, prion diseases, proteasome, carbon metabolism and platelet activation, and DPPs were involved in spliceosome, comprising a dense protein interaction network. The physiological interactions among these protein complexes likely contribute to the cooperation and coordination of their functions in the control of numerous intracellular signaling and regulatory pathways in AS.
Figure 7 Protein–protein interaction (PPI) network analysis of DEPs (A) and DPPs (B). The STRING database was used to annotate the functional interactions of all the identified DEPs and DPPs.
Verification Results of Target Proteins
PRM for targeted quantitative proteins by MS was used to quantitatively identify target proteins. A total of 19 proteins were quantitatively identified (Table 1). These proteins were selected for their functional significance concluded from proteome analysis, and the ratio of protein abundance varied in a wide range. Among them, 3 proteins (HSP90AB1, HSP90AA1 and HSPA8) were involved in the antigen processing and presentation pathway, of which the highest differential expression was HSP90AA1 with an AS/HC ratio of 35.40. Six proteins (including ITPR1, MYLK and STIM1) were involved in the platelet activation pathway, of which the highest differential expression was observed for STIM1 with an AS/HC ratio of 22.35. Ten proteins (including MYL12A, MYL9 and ROCK2) were involved in the leukocyte transendothelial migration pathway, of which the highest differential expression was MMP9 with an AS/HC ratio of 18.97. The fold-change values for these proteins according to the PRM results are presented in Table 1. The PRM and label-free quantification (LQ) results were highly consistent.
Discussion
AS is a chronic systemic inflammatory disorder that mainly affects the sacroiliac joint and spine (2). Due the position limitation of AS, studies of AS-related proteins have been mainly restricted to arthritogenic proteins (11, 19, 20). Additional potential proteins should be investigated to promote the understanding of the pathological mechanism of AS.
In the present study, we performed comprehensive proteomic profiling and PRM analysis of PBMCs obtained from patients with AS. We discovered 782 DEPs and 122 DPPs in AS patients compared to healthy controls. PRM analysis for 19 proteins was performed to further verify the expression levels of target proteins. As a result, 3 proteins (HSP90AB1, HSP90AA1 and HSPA8) involved in antigen processing and presentation, 6 proteins (including ITPR1, MYLK and STIM1) involved in platelet activation and 10 proteins (including MYL12A, MYL9 and ROCK2) involved in leukocyte transendothelial migration were validated to be overexpressed in the PBMCs of AS patients.
Chronic inflammation of the central axis joints is the most significant pathological alterationin in AS, and the disease might gradually extended to internal organs (21). The inflammatory process in AS is still poorly understood, making it difficult to track disease activity and choose treatment strategies. The body tightly regulates inflammation, which is associated with transient crossing of leukocytes through the blood vessel wall, a process known as transendothelial migration. (TEM) (22). The CD18 integrin family plays a crucial role in guiding leukocytes to inflammatory areas (23, 24). Kragstrup et al. (25) demonstrated that that patients with arthritis contain anti-inflammatory soluble CD18 (sCD18) complexes in their blood and synovial fluid. MMP-3 and MMP-9 are involved in CD18 shedding processes (26, 27). MMP family proteins play a role in arthritis pathogenesis (28). In particular, MMP3 is a protease that issecreted and synthesized by fibroblasts and chondrocytes in synovial joints, and it activates other MMPs, such as MMP1, MMP7 and MMP9 (29). In AS, both MMP-3 and MMP-9 serum concentrations have been linked to disease activity (30, 31). Our findings were supported by the discovery of MMP9, a previously studied AS-related protein.
In clinical practice, platelet counts generally increase during the active stage of AS (32, 33). However, there are few studies on platelet function, particularly the activation function of platelets. CD62P, a membrane glycoprotein, is expressed on the α-granule membrane of standstill platelets when they are activated, according to a previous study, while CD63, a lysosome membrane glycoprotein, is also expressed on the platelet membrane when platelets are activated. CD63 and CD62P are platelet activation indicators, with the latter serving as the gold standard (10). Another study discovered that AS patients have significantly higher levels of platelet-monocyte complex (PMC) and soluble CD40L (sCD40L) than the control group (34). In the present study, we identified 6 key proteins involved in platelet activation that were highly expressed in PBMCs of AS patients. Based on previous studies and our findings, platelet activation could indicate an aggravation of AS.
Antigen processing and presentation play a critical role in the pathophysiology of AS (35–37). Human leukocyte antigens (HLAs) are proteins encoded by genes of the human major histocompatibility complex (MHC). The expression of HLA-B27 is highly linked with the presence of AS (38, 39), and other genes related to AS have been discovered (40, 41). In our previous studies, we identified differences in the expression and phosphorylation of HLA-E, calnexin (CANX), and two heat-shock proteins (HSPs) (HSP90AA1 and HSP90AB1) in PBMCs of AS patients and healthy controls (16). In the present study, we further verified these proteins by PRM and demonstrated that HSP90AA1 and HSP90AB1 were highly expressed, indicating that HSP90AA1 and HSP90AB1 may be involved in the development of AS.
Numerous limitations should be addressed in this investigation. First, The present study was a cross-sectional investigation of a modest sample size. Second, despite the fact that thousands of proteins were detected, no more than 1000 of them have quantitative data on their expression and phosphorylation. Some key proteins may have been overlooked due to missing key information. Future research should delve deeper into these limiting factors.
In conclusion, between AS patients and healthy controls, quantitative proteomic studies revealed differently expressed and phosphorylated proteins.In the PBMCs of patients with AS, key proteins involved in antigen processing and presentation, platelet activation, and leukocyte transendothelial migration were shown to be highly elevated. These proteins could be served as candidate markers for AS diagnosis and new therapeutic targets, as well as elucidating the pathophysiology of AS.
Conclusion
Liquid chromatography-tandem mass spectrometry (LC–MS/MS) analysis was performed and discovered 782 differentially expressed proteins (DEPs) and 122 differentially phosphorylated proteins (DPPs) between 9 new-onset AS patients and 9 healthy controls. The key proteins involved in antigen processing and presentation, platelet activation and leukocyte transendothelial migration revealed abnormal immune regulation in patients with new-onset AS. These proteins may further help reveal the pathogenesis of AS and serve as candidate markers for AS diagnosis and new treatment targets.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repositories and accession numbers can be found below: http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD019667.
Ethics Statement
All procedures were approved by the local ethical committee of The Affiliated Hospital of Jinan University and followed the tenets of the Declaration of Helsinki.
Author Contributions
YD, BH, and DT designed the study. ZY and FZ generated and analysed the data. ZY and XH interpreted the data and wrote the manuscript. All the authors revised the manuscript for intellectual content and approved its final version to be published. BH is the guarantor of this work.
Funding
This work was supported by the Science and Technology Plan of Shenzhen (No. JCYJ20170307095606266); Sanming Project of Medicine in Shenzhen, the group of Rheumatology and Immunology leaded by Xiaofeng Zeng of Peking Union Medical College Hospital and Dongzhou Liu in Shenzhen People’s Hospital (SYJY201704 and SYJY201705); the Key Research and Development Program of Guangdong Province (Grant No. 2019B020229001); Shenzhen Key Medical Discipline Construction Fund (No.SZXK011); Postdoctoral Fund of the First Affiliated Hospital, Jinan University (Grant No. 809001); Young Innovative Talents Project of General Colleges and Universities in Guangdong Province (Grant No. 2018KQNCX010) and GuangDong Basic and Applied Basic Research Foundation (Grant No. 2020A1515111209).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Li Chen from Hangzhou LC-BIO Co., Ltd for help with the data analysis and programming in R.
References
1. Feltelius N, Lindh U, Venge P, Hallgren R. Ankylosing Spondylitis: A Chronic Inflammatory Disease With Iron Overload in Granulocytes and Platelets. Ann Rheum Dis (1986) 45(10):827–31. doi: 10.1136/ard.45.10.827
2. Baraliakos X, Listing J, Rudwaleit M, Sieper J, Braun J. The Relationship Between Inflammation and New Bone Formation in Patients With Ankylosing Spondylitis. Arthritis Res Ther (2008) 10(5):R104. doi: 10.1186/ar2496
3. Zhu J, Huang F. The Relationship Between Bone Inflammation and Structural Damage in Ankylosing Spondylitis. Zhonghua Nei Ke Za Zhi (2011) 50(2):89–91. doi: 10.3760/cma.j.issn.0578-1426.2011.02.001
4. Alvarez I, Lopez de Castro JA. HLA-B27 and Immunogenetics of Spondyloarthropathies. Curr Opin Rheumatol (2000) 12(4):248–53. doi: 10.1097/00002281-200007000-00003
5. Akassou A, Bakri Y. Does HLA-B27 Status Influence Ankylosing Spondylitis Phenotype? Clin Med Insights Arthritis Musculoskelet Disord (2018) 11:1179544117751627. doi: 10.1177/1179544117751627
6. Losinska K, Korkosz M, Kwasny-Krochin B. Endothelial Dysfunction in Patients With Ankylosing Spondylitis. Reumatologia (2019) 57(2):100–5. doi: 10.5114/reum.2019.84815
7. Sulicka J, Surdacki A, Korkosz M, Mikolajczyk T, Strach M, Klimek E, et al. Endothelial Dysfunction is Independent of Inflammation and Altered CCR7 T Cell Expression in Patients With Ankylosing Spondylitis. Clin Exp Rheumatol (2017) 35(5):844–9.
8. Herter JM, Rossaint J, Zarbock A. Platelets in Inflammation and Immunity. J Thromb Haemost (2014) 12(11):1764–75. doi: 10.1111/jth.12730
9. Örüm H, Pamuk GE, Pamuk ÖN, Demir M, Turgut B. Does Anti-Tnf Therapy Cause Any Change in Platelet Activation in Ankylosing Spondylitis Patients? J Thromb Thrombolysis (2011) 33(2):154–9. doi: 10.1007/s11239-011-0663-9
10. Wang F, Yan C-G, Xiang H-Y, Xing T, Wang N-S. The Significance of Platelet Activation in Ankylosing Spondylitis. Clin Rheumatol (2008) 27(6):767–9. doi: 10.1007/s10067-008-0847-7
11. Liu J, Zhu P, Peng J, Li K, Du J, Gu J, et al. Identification of Disease-Associated Proteins by Proteomic Approach in Ankylosing Spondylitis. Biochem Biophys Res Commun (2007) 357(2):531–6. doi: 10.1016/j.bbrc.2007.03.179
12. Cai A, Qi S, Su Z, Shen H, Yang Y, He L, et al. Quantitative Proteomic Analysis of Peripheral Blood Mononuclear Cells in Ankylosing Spondylitis by iTRAQ. Clin Trans Sci (2015) 8(5):579–83. doi: 10.1111/cts.12265
13. Via A, Zanzoni A. A Prismatic View of Protein Phosphorylation in Health and Disease. Front Genet (2015) 6:131. doi: 10.3389/fgene.2015.00131
14. Lenz C. Identification of Protein Phosphorylation Sites by Advanced LC-ESI-MS/MS Methods. Methods Mol Biol (2019) 1934:163–78. doi: 10.1007/978-1-4939-9055-9_11
15. Wright C, Edelmann M, diGleria K, Kollnberger S, Kramer H, McGowan S, et al. Ankylosing Spondylitis Monocytes Show Upregulation of Proteins Involved in Inflammation and the Ubiquitin Proteasome Pathway. Ann Rheumatic Dis (2008) 68(10):1626–32. doi: 10.1136/ard.2008.097204
16. Lu YP, Zhang XL, Zheng F, Yun C, Zhu C, Cai W, et al. Quantitative Proteomic Analyses To Reveal the Key Features of Proteins in New Onset Ankylosing Spondylitis Patients. ACS Omega (2020) 5(32):20153–61. doi: 10.1021/acsomega.0c01776
17. van der Linden S, Valkenburg HA, Cats A. Evaluation of Diagnostic Criteria for Ankylosing Spondylitis. A Proposal for Modification of the New York Criteria. Arthritis Rheum (1984) 27(4):361–8. doi: 10.1002/art.1780270401
18. Atallah M, Flory MR, Mallick P. A Robust Protocol for Protein Extraction and Digestion. Methods Mol Biol (2017) 1550:1–10. doi: 10.1007/978-1-4939-6747-6_1
19. Atagunduz P, Appel H, Kuon W, Wu P, Thiel A, Kloetzel PM, et al. HLA-B27-Restricted CD8+ T Cell Response to Cartilage-Derived Self Peptides in Ankylosing Spondylitis. Arthritis Rheum (2005) 52(3):892–901. doi: 10.1002/art.20948
20. Chang X, Han J, Zhao Y, Yan X, Sun S, Cui Y. Increased Expression of Carbonic Anhydrase I in the Synovium of Patients With Ankylosing Spondylitis. BMC Musculoskelet Disord (2010) 11:279. doi: 10.1186/1471-2474-11-279
21. Schett G. Independent Development of Inflammation and New Bone Formation in Ankylosing Spondylitis. Arthritis Rheum (2012). doi: 10.1002/art.34441
22. Geng JG. Directional Migration of Leukocytes: Their Pathological Roles in Inflammation and Strategies for Development of Anti-Inflammatory Therapies. Cell Res (2001) 11(2):85–8. doi: 10.1038/sj.cr.7290071
23. Gjelstrup LC, Boesen T, Kragstrup TW, Jorgensen A, Klein NJ, Thiel S, et al. Shedding of Large Functionally Active CD11/CD18 Integrin Complexes From Leukocyte Membranes During Synovial Inflammation Distinguishes Three Types of Arthritis Through Differential Epitope Exposure. J Immunol (2010) 185(7):4154–68. doi: 10.4049/jimmunol.1000952
24. El Kebir D, Jozsef L, Pan W, Filep JG. Myeloperoxidase Delays Neutrophil Apoptosis Through CD11b/CD18 Integrins and Prolongs Inflammation. Circ Res (2008) 103(4):352–9. doi: 10.1161/01.RES.0000326772.76822.7a
25. Kragstrup TW, Jalilian B, Hvid M, Kjaergaard A, Ostgard R, Schiottz-Christensen B, et al. Decreased Plasma Levels of Soluble CD18 Link Leukocyte Infiltration With Disease Activity in Spondyloarthritis. Arthritis Res Ther (2014) 16(1):R42. doi: 10.1186/ar4471
26. Gomez IG, Tang J, Wilson CL, Yan W, Heinecke JW, Harlan JM, et al. Metalloproteinase-Mediated Shedding of Integrin Beta2 Promotes Macrophage Efflux From Inflammatory Sites. J Biol Chem (2012) 287(7):4581–9. doi: 10.1074/jbc.M111.321182
27. Vaisar T, Kassim SY, Gomez IG, Green PS, Hargarten S, Gough PJ, et al. MMP-9 Sheds the Beta2 Integrin Subunit (CD18) From Macrophages. Mol Cell Proteomics (2009) 8(5):1044–60. doi: 10.1074/mcp.M800449-MCP200
28. Meng C, Bai R, Zhao Z, Huang G, Jin T, Feng W, et al. MMP-8 Single-Nucleotide Polymorphisms are Related to Ankylosing Spondylitis in Chinese Han Population. Med (Baltimore) (2018) 97(35):e12136. doi: 10.1097/MD.0000000000012136
29. Lerner A, Neidhofer S, Reuter S, Matthias T. MMP3 is a Reliable Marker for Disease Activity, Radiological Monitoring, Disease Outcome Predictability, and Therapeutic Response in Rheumatoid Arthritis. Best Pract Res Clin Rheumatol (2018) 32(4):550–62. doi: 10.1016/j.berh.2019.01.006
30. Chen CH, Lin KC, Yu DT, Yang C, Huang F, Chen HA, et al. Serum Matrix Metalloproteinases and Tissue Inhibitors of Metalloproteinases in Ankylosing Spondylitis: MMP-3 Is a Reproducibly Sensitive and Specific Biomarker of Disease Activity. Rheumatol (Oxford) (2006) 45(4):414–20. doi: 10.1093/rheumatology/kei208
31. Mattey DL, Packham JC, Nixon NB, Coates L, Creamer P, Hailwood S, et al. Association of Cytokine and Matrix Metalloproteinase Profiles With Disease Activity and Function in Ankylosing Spondylitis. Arthritis Res Ther (2012) 14(3):R127. doi: 10.1186/ar3857
32. Beyan C, Beyan E. Is Mean Platelet Volume and Inflammatory Activity Really Correlated in Patients With Ankylosing Spondylitis? Mod Rheumatol (2020) 30(2):410. doi: 10.1080/14397595.2019.1660029
33. Yazici S, Yazici M, Erer B, Erer B, Calik Y, Bulur S, et al. The Platelet Functions in Patients With Ankylosing Spondylitis: Anti-TNF-Alpha Therapy Decreases the Mean Platelet Volume and Platelet Mass. Platelets (2010) 21(2):126–31. doi: 10.3109/09537100903470306
34. Orum H, Pamuk GE, Pamuk ON, Demir M, Turgut B. Does Anti-TNF Therapy Cause Any Change in Platelet Activation in Ankylosing Spondylitis Patients? A Comparative Study. J Thromb Thrombolysis (2012) 33(2):154–9. doi: 10.1007/s11239-011-0663-9
35. Wang G, Kim TH, Li Z, Cortes A, Kim K, Bang SY, et al. MHC Associations of Ankylosing Spondylitis in East Asians are Complex and Involve non-HLA-B27 HLA Contributions. Arthritis Res Ther (2020) 22(1):74. doi: 10.1186/s13075-020-02148-5
36. Bettencourt BF, Rocha FL, Alves H, Amorim R, Caetano-Lopes J, Vieira-Sousa E, et al. Protective Effect of an ERAP1 Haplotype in Ankylosing Spondylitis: Investigating non-MHC Genes in HLA-B27-Positive Individuals. Rheumatol (Oxford) (2013) 52(12):2168–76. doi: 10.1093/rheumatology/ket269
37. Sims AM, Barnardo M, Herzberg I, Bradbury L, Calin A, Wordsworth BP, et al. Non-B27 MHC Associations of Ankylosing Spondylitis. Genes Immun (2007) 8(2):115–23. doi: 10.1038/sj.gene.6364362
38. Arevalo M, Masmitja JG, Moreno M, Calvet J, Orellana C, Ruiz D, et al. Correction to: Influence of HLA-B27 on the Ankylosing Spondylitis Phenotype: Results From the REGISPONSER Database. Arthritis Res Ther (2020) 22(1):140. doi: 10.1186/s13075-020-02237-5
39. Chen B, Li J, He C, Li D, Tong W, Zou Y, et al. Role of HLA-B27 in the Pathogenesis of Ankylosing Spondylitis (Review). Mol Med Rep (2017) 15(4):1943–51. doi: 10.3892/mmr.2017.6248
40. Abdullah H, Zhang Z, Yee K, Haroon N. KIR3DL1 Interaction With HLA-B27 is Altered by Ankylosing Spondylitis Associated ERAP1 and Enhanced by MHC Class I Cross-Linking. Discov Med (2015) 20(108):79–89.
Keywords: ankylosing spondylitis (AS), global proteomic, immune regulation, peripheral blood mononuclear cells (PBMC), phosphorylation
Citation: Yu Z, Hong X, Zhang X, Zheng F, Liu F, Xu H, Zhu C, Cai W, Liu D, Yin L, Hu B, Tang D and Dai Y (2022) Global Proteomic Analyses Reveals Abnormal Immune Regulation in Patients With New Onset Ankylosing Spondylitis. Front. Immunol. 13:838891. doi: 10.3389/fimmu.2022.838891
Received: 18 January 2022; Accepted: 23 February 2022;
Published: 17 March 2022.
Edited by:
Michael V. Volin, Midwestern University, United StatesCopyright © 2022 Yu, Hong, Zhang, Zheng, Liu, Xu, Zhu, Cai, Liu, Yin, Hu, Tang and Dai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Dai, ZGFpeW9uZzIyQGFsaXl1bi5jb20=; Bo Hu, NDIwODk1MzdAcXEuY29t; Donge Tang, ZG9uZ2U2NkAxMjYuY29t
†These authors have contributed equally to this work
‡These authors share first authorship