Identification of thrombotic biomarkers in orthopedic surgery patients by plasma proteomics
Journal of Orthopaedic Surgery and Research volume 18, Article number: 222 (2023)
Due to the poor specificity of D-dimer, more accurate thrombus biomarkers are clinically needed to improve the diagnostic power of VTE.
The plasma samples were classified into low-risk group (n = 6) and high-risk group (n = 6) according to the Caprini Thrombosis Risk Assessment Scale score. Data-independent acquisition mass spectrometry (DIA-MS) was performed to identify the proteins in the 12 plasma samples. Bioinformatics analysis including volcano plot, heatmap, KEGG pathways and chord diagram analysis were drawn to analyze the significantly differentially expressed proteins (DEPs) between the two groups. Then, another 26 plasma samples were collected to verify the key proteins as potential biomarkers of VTE in orthopedic surgery patients.
A total of 371 proteins were identified by DIA-MS in 12 plasma samples. Volcano plotting showed that there were 30 DEPs. KEGG pathway enrichment analysis revealed that the DEPs were majorly involved in the blood coagulation pathway. The chord diagram analysis demonstrated that proteins SAA1, VWF, FLNA, ACTB, VINC, F13B, F13A and IPSP in the DEPs were significantly related to blood coagulation. VWF and F13B were selected for validation experiments. ELISA test showed that, as compared with those in the low-risk group, the level of VWF in the high-risk sera was significantly increased.
The level of VWF in the high-risk group of thrombosis after orthopedic surgery was significantly higher than that in the low-risk group of preoperative thrombosis, suggesting that VWF may be used as a potential thrombus biomarker in orthopedic surgery patients.
Venous thromboembolism (VTE) is the leading reason of cardiovascular death worldwide. The annual incidence of VTE was estimated to be approximately 1–2 per 1000 population in Europe and the USA, and the incidence increases with age growth [1,2,3]. Approximately two-thirds of patients with symptomatic VTE manifest deep vein thrombosis alone (DVT), while one third has pulmonary embolism (PE) [3,4,5]. According to Virchow triad, the main factors for the formation of thrombus include stasis of blood flow, endothelial injury and hypercoagulability . In orthopedic patients, bone fracture can lead to venous injury. Immobilization and bed rest after surgery cause blood stasis in the lower extremities of the patients. In addition, fasting, anesthesia, hemostatic drugs, intraoperative blood loss and postoperative fever lead to the patients' blood hypercoagulability. These factors result in high incidence of VTE in patients after orthopedic surgery. More than 10% incidence of VTE has been reported after knee or hip joint replacement surgery . And PE contributes to a mortality rate of 70% among patients . Unfortunately, there may be no obvious signs of thrombosis, and thrombosis can occur at any age and at any time, seriously threatening the life and health of patients. Therefore, early diagnosis and prevention of thrombosis are crucial.
Currently, evaluating whether a patient suffers from VTE is mainly based on the patient's symptoms, signs, and auxiliary examinations. Computed tomography pulmonary angiography (CTPA) has been applied as the gold standard for diagnosing PE. However, radiation exposure and expensive price restrict the application of the diagnosis . Featuring with high sensitivity and accuracy, color Doppler ultrasound has been applied to the diagnosis of DVT [10, 11]. However, the guiding value for clinical doctors to prevent thrombosis is limited. As a final degradation product of cross-linked fibrin, D-dimer is often elevated in patients who are suffering from acute VTE . D-dimer is a biomarker of coagulation activation and secondary hyperfibrinolysis, and its elevated value does not necessarily indicate thrombosis . Obviously, the specificity of D-dimer for the detection of VTE is poor [14, 15]. Recent studies have attempted to identify new biomarkers that can act as a complement to D-dimer [16, 17]. However, the value of these biomarkers is susceptible to many factors such as acute coronary syndromes, local or systemic infection, trauma and malignancy . Therefore, there is a need for more accurate biomarkers for predicting VTE.
Proteomic technology is an efficient tool for identifying potential biomarkers associated with pathological states [19, 20]. With the rapid development of mass spectrometry (MS), proteomics can be used to accurately profile global proteomic alterations in various cells, tissues and body fluids. By applying proteomic techniques to serology, more potential clinical plasma biomarkers associated with venous thromboembolism will likely be discovered. With high reproducibility and sensitivity, data-independent acquisition mass spectrometry (DIA-MS)  can generate digital proteomic maps, providing highly reproducible retrospective analysis on cell and tissue samples and thus holding great promise for biomarker discovery.
In this work, DIA-MS proteomics was performed to quantify and analyze the expression of proteins in 12 plasma samples that were classified into two groups with low and high risk of VTE according to Caprini score. Bioinformatic tools including volcano plot, heatmap, KEGG pathways enrichment and chord diagram analysis were employed to sort out the key proteins involved in related pathways associated with thrombus. Then, ELISA experiment was applied to another 26 plasma samples to verify the key proteins as potential biomarkers of VTE in orthopedic surgery patients.
Materials and methods
All plasma samples came from the patients who were admitted to our department from January 2019 to January 2021. The basic characteristics of these patients were collected. And they were classified into two groups according to Caprini Thrombosis Risk Assessment Scale . A total score of 0 to 1 was defined as low risk of thrombosis, 2 as intermediate, 3 to 4 as high and ≥ 5 as very high. All the work described in this manuscript has been carried out in accordance with the Code of Ethics of the World Medical Association. Approval for this study was obtained from the institutional review board of our institution. All samples from patients were analyzed under the same procedures (Fig. 1).
Preparation of samples for DIA-MS analysis
Blood samples were collected into 5-ml anticoagulant tubes with ethylenediaminetetraacetic acid (EDTA), and, immediately, all sample tubes were homogenized upside down. Afterward, all blood specimens collected from the patients were centrifuged at 1000 g for 5 min at 4 °C, and then the respective supernatants were collected to a new corresponding centrifuge tube, which were stored at − 80 °C before further processing. Before DIA-MS analysis, the frozen plasma samples were thawed at 4 °C. The concentration of the plasma protein was detected by BCA protein assay kit (Beyotime, P0012S). Urea buffer (8 M, 500 µL) was added to 500 μL plasma samples, and the mixture was vortexed for 30 s; then, each sample was added with DTT (50 mM) and fully vortexed. After 1 h in 37 °C, IAA was added in the samples and stood for 30 min at room temperature in a dark environment. The above mixed solution was then filtered with a flat-bottom ultrafiltration tube. The filtrate was collected in a new centrifuge tube. TEAB (50 mM) and trypsin (1 μg/μL) were orderly added in the filtrate. After 16–18 h in 37 °C, the samples were centrifuged at 12,000 g for 15 min. Finally, freeze-dried protein powder was obtained in a freeze dryer. All the freeze-dried protein powders were stored at − 80 °C. For DIA-MS analysis, the protein powders were dissolved in appropriate gold water containing 0.1% FA to get the samples with the concentration of 0.5 μg/μL. Then the samples were centrifuged twice at 12,000 g for 10 min each time, and the supernatants were transferred to a new centrifuge tube. The labeled peptide (5 × iRT) was added to the centrifuge tube. Then the tubes were fully vortexed and centrifuged twice at 12,000 g for 10 min. The supernatants were collected and applied to DIA-MS analysis.
The collected supernatants in “Preparation of samples for DIA-MS analysis” secttion were taken for DIA-MS analysis. All samples were separated by a nano-flow HPLC liquid phase system Easy nLC2000. Liquid phase parameters: (a) column information: 300 µm idx 5 mm, Acclaim PepMap RSLC C18, 5 µm, 100 Å (Thermo, 160,454); Acclaim PepMap 75 µm × 150 mm C18, 3 µm, 100 Å (Thermo, 160,321). (b) Mobile phase information: Mobile phase A: 0.1% formic acid. Mobile phase B: 0.1% formic acid, 80% ACN. Flow rate = 300 nL/min. (c) Analysis time = 65 min. Effective gradient: B phase rose from 5 to 90%; 0–5 min: the linear gradient of liquid B was 0–5%; 5–45 min: the linear gradient of liquid B was 5–50%; 45–50 min, the linear gradient of liquid B was 50–90%; 50–55 min: the linear gradient of liquid B maintained at 90%; 55–65 min: the linear gradient of liquid B decreased to 5%. Mass spectrometry parameters were as follows: (a) Primary mass spectrometry parameters: resolution = 70,000; AGC target = 3e6; maximum IT = 40 ms; scan range = 350 to 1800 m/z. (b) MS parameters: resolution = 17,500; AGC target = 1e5; maximum IT = 60 ms; Top N = 20; NCE/stepped NCE = 27.
Data processing and analysis
The original data were transformed into.MzML format by Trans-Proteomic Pipeline software which also was used for database search. The theoretical protein sequence database is the Uniport human nonredundant library.fasta format file; the search database used the comet search engine and the corresponding file.params format; no decoy was added, the whole enzyme was used, and the trypsin missed cleavage site was set to 0. The fixed modified cysteine was alkylated (57.02146 Da), and the non-fixed modified methionine was oxidized (15.9949 Da). The error range of precursor ions was controlled within 10–5, and the error of secondary debris was controlled within 0.02 Da. The 12 database search results were integrated and constructed by TPP and imported into Skyline for data analysis of DIA.
DIA data acquisition was performed on 12 samples, and the mass scanning range m/z 400–1200 was equally divided into 32 consecutive 25 u windows. In each window, all product ion information of all precursor ions was selected and repeated. The Skyline software was used for quantitative analysis of label-free proteomics. DIA data were imported to screen and extract protein quantitative information. The automatic auxiliary function of Skyline was used to screen suitable chromatographic peaks. Information on protein name, peptide sequence, and product ion area was exported.
Two DIA quantification data were normalized by total ionic strength (TIC) via R (version 3.2.2) software, and the average of the two results was used for subsequent analysis. The peptides with ion quantitative correlation of more than 0.6 in the two groups of samples were retained. The product ion areas of the same peptides were summed to the peptide level, and two-tailed t-test was performed at the peptide level to screen peptides with significant difference (P < 0.05). The different peptides were normalized to the protein level to get DEPs. Then, the R language was applied to draft the volcano plot for visualizing the expression of proteins in plasma and Heatmap for cluster analysis. The data of DEPs were uploaded to the KEGG database (http://www.genome.jp/kegg/) for enrichment analysis of metabolic pathways. Chordogram analysis of GO enrichment was performed to present the DEPs in three pathways of blood coagulation, humoral immune response and acute phase response.
Another 26 plasma samples of each group were collected from patients using the method described in “Preparation of samples for DIA-MS analysis” section. Levels of VWF and F13B in the plasma samples were measured using ELISA kit according to the instruction manual.
All date were imported into GraphPad prism 8.3.0 for analysis. All results of ELISA were imported into SPSS 22.0 software (IBM, Armonk, NY, USA) for statistical analysis. Measurement data are presented as mean ± standard deviation (SD). One-way ANOVA was used for statistical analysis. A value of P < 0.05 was considered statistically difference.
Patient baseline characteristics
The basic information of the 6 patients is shown in Table 1. The information of age, gender, height, weight, BMI, hypertension, diabetes, diagnosis and Carprini score are presented in the table. Apart from the factor of orthopedic surgery, there is no any difference between the samples of two groups.
The results of DIA-MS
A total of 371 proteins were identified in 12 plasma samples and 30 DEPs were found between the two sample groups before and after orthopedic surgery. The volcano plot (Fig. 2) and the heat map analysis (Fig. 3) of the 30 DEPs showed that the expression of 14 proteins was upregulated and 16 proteins were downregulated in the high-risk group of thrombosis as compared with the low-risk group of thrombosis.
Next, we performed KEGG pathway enrichment analysis on the DEPs. The results showed that the top ten pathways participated by the DEPs were (listed in descending order), blood coagulation, hemostasis, coagulation, platelet degranulation, humoral immune response, neutrophils Cell degranulation, platelet activation, complement activation, acute phase response, and alternative pathways of complement activation (Fig. 4). It can be seen that blood coagulation pathway was mainly involved.
Then, we performed chord diagram analysis on the DEPs, by focusing on three relevant pathways of blood coagulation, humoral immune response and acute phase response. We found that the proteins SAA1, VWF, FLNA, ACTB, VINC, F13B, F13A and IPSP were clearly associated with blood coagulation (Fig. 5). Among them, the expressions of VWF and SAA1 were significantly upregulated, while the expressions of F13B, F13A and IPSP were significantly downregulated.
The results of ELISA
We then performed ELISA experiments to validate the results obtained from our proteomic analysis. Another 26 pairs of plasma samples were thus collected in our department (Additional file 1: Tables S1, S2) for this validation experiments. ELISA testing was carried out on each sample to measure the concentration of VWF and F13B, respectively, selected from the up- and downregulated DEPs. The results showed that, as compared with those in the low-risk group, the content of VWF in the high-risk sera was significantly increased, while the expression of F13B did not alter significantly. This validation suggests that VWF protein may serve as a potential orthopedic surgical thrombosis biomarker (Fig. 6).
VTE is a complication with a high incidence rate after major orthopedic surgery, and is also one of the important factors of perioperative death and unexpected death in hospital [23, 24]. It is necessary to apply a more effective biomarker characterized by high sensitivity and specificity to guide orthopedic surgeons to prevent thrombosis. Currently, color Doppler ultrasound and D-dimer are used to evaluate the risk of thrombosis before surgery in clinical practice [25, 26], and their weakness cannot be ignored, however. Recent studies reported that some proteins such as NLR, PLR, MLR and SII in blood can be biomarkers of thrombosis, but their stability and specificity are questioned . It is thus urgently needed to discover suitable biomarkers for clinicians.
Based on the well development of MS technology and proteomics, better identification, characterization and quantification of expressed proteins in plasma become a reality . Here we employed the prevalent DIA-MS method to profile the proteomic alterations in sera with different coagulation risks, aiming to globally screen for potential biomarkers of thrombosis. The identified DEPs were subjected to bioinformatic analyses including KEGG pathway enrichment and chord diagram analysis to better understand the molecular pathways involved in thrombosis. ELISA experiment was then carried out to verify the significant alterations of key proteins. Our results demonstrate that plasma VWF in patients after orthopedic surgery is significantly higher than preoperative VWF and thus may be used as a biomarker of thrombosis in patients with orthopedic surgery.
Thrombosis refers to abnormal blood coagulation in the circulating blood due to certain factors or blood deposits on the inner wall of the heart or blood vessel walls . However, the molecular pathways and key proteins involved in this process are largely unknown. Our KEGG pathway enrichment analysis on the plasma DEPs revealed that coagulation, hemostasis, blood coagulation, platelet degranulation, and platelet activation are the major pathways affected in thrombosis. These pathways are involved in different stages of the thrombosis process . And the key proteins that play functional roles in these pathways may proceed alteration in expression level and thus become direct indicators of thrombosis.
To sort out the key proteins from the plasma DEPs, we performed chord diagram analysis and found that SAA1, VWF, FLNA, ACTB, VINC, F13B, F13A and IPSP, from the pathways of blood coagulation, humoral immune response and acute phase response were clearly associated with thrombosis. Among them, the expressions of SAA1 and VWF were significantly upregulated, while the expressions of F13B, F13A and IPSP were significantly downregulated. And our ELISA experiment verified that VWF does have a higher expression in high-risk sera and thus could serve as a thrombus biomarker in orthopedic surgery patients.
VWF (von Willebrand factor) is a major determinant of hemostasis and clot formation in arteries and veins . VWF deficiency can lead to von Willebrand disease (VWD), and VWF overactivition can lead to thrombotic thrombocytopenic purpura (TTP) [31, 32]. In the process of hemostasis, VWF binds to the platelet membrane GPIb-IX complex and subendothelial collagen to mediate the adhesion of platelets at the site of vascular injury . As a carrier, VWF can stabilize factor VIII via binding to factor VIII . VWF and some hemostatic and fibrinolytic components are independent risk factors for thrombotic disease . It has been also tested that VWF may be a marker for thrombotic risk in COVID-19 . In this regard, our current finding provided new experimental evidences to support these observations. The high expression of VWF in the sera can be applied to evaluate the risk of thrombus in orthopedic patients. Due to the important function of VWF in hemostasis, further investigation on its functional mechanism may provide possibility of intervention in thrombosis and sustain its role as a potential biomarker for VTE.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Deep vein thrombosis
Computed tomography pulmonary angiography
Data-independent acquisition mass spectrometry
Differentially expressed proteins
Von Willebrand factor
Cho HJ, Lee IK, Lee YS, et al. Application of venous thromboembolism prophylaxis program in patients with colorectal cancer using the enhanced recovery after surgery protocol. Eur J Surg Oncol. 2022;48:1384–9.
Chen J, Lv M, Wu S, et al. Severe bleeding risks of direct oral anticoagulants in the prevention and treatment of venous thromboembolism: a network meta-analysis of randomised controlled trials. Eur J Vasc Endovasc Surg. 2022;63(3):465–74.
Khan F, Tritschler T, Kahn SR, et al. Venous thromboembolism. Lancet. 2021;398(10294):64–77.
Wendelboe AM, Raskob GE. Global burden of thrombosis: epidemiologic aspects. Circ Res. 2016;118(9):1340–7.
White RH. The epidemiology of venous thromboembolism. Circulation. 2003;107(23 Suppl 1):I4-8.
Hackett TR, Godin JA. Editorial commentary: should the virchow triad have been a quartet? Is high altitude a risk factor for deep venous thrombosis after knee arthroscopy? Arthroscopy. 2016;32(11):2355–6.
Falck-Ytter Y, Francis CW, Johanson NA, et al. Prevention of VTE in orthopedic surgery patients: antithrombotic therapy and prevention of thrombosis, 9th ed: American college of chest physicians evidence-based clinical practice guidelines. Chest. 2012;141(2 Suppl):e278S-e325S.
Xie J, Jiang M, Lin Y, et al. Rivaroxaban versus aspirin in prevention of venous thromboembolism: a meta-analysis of 9 randomized controlled trials comprising 7,656 patients. Thromb Haemost. 2019;119(9):1517–26.
Kearon C, De Wit K, Parpia S, et al. Diagnosis of pulmonary embolism with d-dimer adjusted to clinical probability. N Engl J Med. 2019;381(22):2125–34.
Permpikul C, Chaiyasoot W, Panitchote A. Incidence of proximal deep vein thrombosis in medical critical care patients. Thromb J. 2022;20(1):5.
Sturm L, Bettinger D, Klinger C, et al. Validation of color Doppler ultrasound and computed tomography in the radiologic assessment of non-malignant acute splanchnic vein thrombosis. PLoS ONE. 2021;16(12): e0261499.
Di Minno MND, Calcaterra I, Papa A, et al. Diagnostic accuracy of D-Dimer testing for recurrent venous thromboembolism: a systematic review with meta-analysis: VTE recurrence and D-dimer. Eur J Intern Med. 2021;89:39–47.
Kearon C, Parpia S, Spencer FA, et al. D-dimer levels and recurrence in patients with unprovoked VTE and a negative qualitative D-dimer test after treatment. Thromb Res. 2016;146:119–25.
Di Nisio M, Van Es N, Büller HR. Deep vein thrombosis and pulmonary embolism. Lancet. 2016;388(10063):3060–73.
Xue J, Ma D, Jiang J, et al. Diagnostic and prognostic value of immune/inflammation biomarkers for venous thromboembolism: Is it reliable for clinical practice? J Inflamm Res. 2021;14:5059–77.
Teperman J, Carruthers D, Guo Y, et al. Relationship between neutrophil-lymphocyte ratio and severity of lower extremity peripheral artery disease. Int J Cardiol. 2017;228:201–4.
Köse N, Yıldırım T, Akın F, et al. Prognostic role of NLR, PLR, and LMR in patients with pulmonary embolism. Bosn J Basic Med Sci. 2020;20(2):248–53.
Özkaya M, Bedel C. Can neutrophil to lymphocyte ratio and platelet to lymphocyte ratio predict the severity of sternum fractures? Current Thoracic Surgery. 2018;3:116.
Nurmohamed NS, Belo Pereira JP, Hoogeveen RM, et al. Targeted proteomics improves cardiovascular risk prediction in secondary prevention. Eur Heart J. 2022;43(16):1569–77.
Zhen J, Zhao P, Li Y, et al. The multiomics analyses of gut microbiota, urine metabolome and plasma proteome revealed significant changes in allergy featured with indole derivatives of tryptophan. J Asthma Allergy. 2022;15:117–31.
Gebreyesus ST, Siyal AA, Kitata RB, et al. Streamlined single-cell proteomics by an integrated microfluidic chip and data-independent acquisition mass spectrometry. Nat Commun. 2022;13(1):37.
Caprini JA. Risk assessment as a guide to thrombosis prophylaxis. Curr Opin Pulm Med. 2010;16(5):448–52.
Lin HY, Lin CY, Huang YC, et al. Deep vein thrombosis after major orthopedic surgery in Taiwan: a prospective cross-sectional study and literature review. J Formos Med Assoc, 2022.
Streiff MB, Haut ER. The CMS ruling on venous thromboembolism after total knee or hip arthroplasty: weighing risks and benefits. JAMA. 2009;301(10):1063–5.
Li B, Liu S, Wang Y, et al. The influence of carotid atherosclerosis on surgical outcomes of patients with cervical spondylotic myelopathy: a retrospective study. Medicine (Baltimore). 2022;101(5): e28743.
Hao L, Wu Y, Zhang Y, et al. SdrG gene activated joint surface membrane protein PTPRJ to induce periprosthetic joint infection after total hip arthroplasty. J Clin Pharm Ther. 2022;47(4):517–22.
Da Silva BF, Meng C, Helm D, et al. Towards understanding male infertility after spinal cord injury using quantitative proteomics. Mol Cell Proteomics. 2016;15(4):1424–34.
Li X, Guo T, Feng Q, et al. Progress of thrombus formation and research on the structure-activity relationship for antithrombotic drugs. Eur J Med Chem. 2022;228: 114035.
Zhu Y, Liu H, Qin S, et al. Antibacterial sericin cryogels promote hemostasis by facilitating the activation of coagulation pathway and platelets. Adv Healthc Mater. 2022: e2102717.
Mojzisch A, Brehm MA. The manifold cellular functions of von willebrand factor. Cells 2021;10(9):2351.
Xiang Y, Hwa J. Regulation of VWF expression, and secretion in health and disease. Curr Opin Hematol. 2016;23(3):288–93.
Lenting PJ, Christophe OD, Denis CV. von Willebrand factor biosynthesis, secretion, and clearance: connecting the far ends. Blood. 2015;125(13):2019–28.
Shao B, Hoover CM, Shi H, et al. Deletion of platelet CLEC-2 decreases GPIba-mediated integrin aIIbb3 activation and decreases thrombosis in TTP. Blood. 2022;139(16):2523–33.
Javitt G, Khmelnitsky L, Albert L, et al. Assembly Mechanism of Mucin and von Willebrand Factor Polymers. Cell. 2020;183(3):717-729.e716.
Sabbagh Q, André-Grégoire G, Alves-Nicolau C, et al. The von Willebrand factor stamps plasmatic extracellular vesicles from glioblastoma patients. Sci Rep. 2021;11(1):22792.
Ladikou EE, Sivaloganathan H, Milne KM, et al. Von Willebrand factor (vWF): Marker of endothelial damage and thrombotic risk in COVID-19? Clin Med (Lond). 2020;20(5):e178–82.
Ethics approval and consent to participate
All study participants provided informed consent to participate, after which peripheral blood samples were collected. The medical ethics commissioner of the affiliated hospital of Jinan University approved this study (KY-2021-034).
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, CQ., Gao, YJ., Lin, GX. et al. Identification of thrombotic biomarkers in orthopedic surgery patients by plasma proteomics. J Orthop Surg Res 18, 222 (2023). https://doi.org/10.1186/s13018-023-03672-1