Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 15 October 2025

Sec. Cardiovascular Endocrinology

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1632694

Screening of biomarkers for early diagnosis of trauma-induced coagulopathy based on untargeted metabolomics

Xu Xu,&#x;Xu Xu1,2†Qianhui Ouyang&#x;Qianhui Ouyang1†Min ShaoMin Shao1Qi YiQi Yi3Sixiang LiuSixiang Liu3Ying HuangYing Huang3Jia Wang*Jia Wang4*Xingwen Zhang*Xingwen Zhang3*Chaochao Tan*Chaochao Tan1*
  • 1Department of Clinical Laboratory, Hunan Provincial People’s Hospital (The First Affiliated Hospital of Hunan Normal University), Changsha, Hunan, China
  • 2School of Medicine, Hunan Normal University, Changsha, Hunan, China
  • 3Department of Emergency, Hunan Provincial People’s Hospital (The First Affiliated Hospital of Hunan Normal University), Changsha, Hunan, China
  • 4Department of Research, Hunan Provincial People’s Hospital (The First Affiliated Hospital of Hunan Normal University), Changsha, Hunan, China

Objectives: Trauma-induced coagulopathy (TIC) is an acute coagulation disorder characterized by massive bleeding following trauma and is a leading cause of mortality. However, current clinical methods are inadequate for predicting TIC onset, and reliable biomarkers for early diagnosis are lacking. This study aimed to identify potential biomarkers with high sensitivity and specificity for TIC using an untargeted metabolomics approach.

Methods: We analyzed serum samples from 54 trauma patients (27 with TIC and 27 without TIC) and 27 healthy controls. All samples were collected within 24 hours post-trauma. Metabolomic profiling was conducted using liquid chromatography-tandem mass spectrometry (LC-MS).

Results: Metabolite profiles differed significantly between the TIC and non-TIC groups. Two metabolites, LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) (AUC = 0.933, 95% CI: 0.849–0.995) and LysoPE(0:0/18:2(9Z,12Z)) (AUC = 0.916, 95% CI: 0.818–0.914), were identified as potential biomarkers for distinguishing TIC. The diagnostic performance of these metabolites surpassed that of both conventional coagulation tests and admission assessment scores.

Conclusion: Two LysoPE metabolites were identified as promising biomarkers for the early detection of TIC.

1 Introduction

Trauma is a leading cause of mortality worldwide (1). Uncontrolled hemorrhage, which is closely associated with trauma-induced coagulopathy (TIC), is the most common cause of early post-traumatic death (24). TIC is a coagulation disorder syndrome occurring in the early phase following trauma (5), characterized initially by a hypocoagulable state presenting as massive bleeding, followed by a late hypercoagulable phase manifesting as excessive coagulation associated with venous thrombosis and multiple organ failure (6). TIC is significantly correlated with increased mortality, greater transfusion requirements, and multiple organ dysfunction (4, 7). Epidemiological studies indicate that approximately 25%-33% of severely injured patients present with TIC upon hospital admission (8). The case fatality rate is four to six times higher than that of patients without coagulation disorders (9), with its pathogenesis involving multiple factors, including endothelial injury, coagulation factor depletion, hyperfibrinolysis, and metabolic disturbances (6).

Clinical diagnosis of TIC primarily relies on conventional coagulation tests (e.g., PT, APTT). While these methods are well-established, they have notable limitations: prolonged turnaround times and an inability to comprehensively assess thrombin generation and fibrinolytic system status (6, 10). Consequently, identifying early diagnostic biomarkers for TIC using metabolomics holds significant promise for improving patient outcomes and reducing mortality (11, 12).

Studies have shown that traumatic brain injury (TBI) is associated with enriched gene expression in coagulation/complement cascades and neuronal markers, as well as elevated levels of glycolytic metabolites and conjugated bile acids (13). In TBI patients with dural invasion, significant increases were observed in blood metabolites involved in late glycolysis, cysteine and one-carbon metabolism, as well as those related to endothelial dysfunction, arginine metabolism, and hypoxia response (14). Currently, no metabolomics studies have focused specifically on whether trauma patients develop TIC. The pathogenesis of TIC leads to alterations in peripheral blood metabolites, and metabolomics research could facilitate further investigation into TIC's pathological mechanisms, potentially uncovering novel mechanisms. Furthermore, no reliable biomarkers are currently available to identify trauma patients who develop TIC, and the associated metabolic pathways remain largely unexplored. Metabolomics, an interdisciplinary approach that systematically analyzes the dynamic changes of small molecule metabolites within biological systems, can reveal aberrant metabolic pathways underlying disease initiation and progression, and has been widely employed to explore diagnostic or prognostic biomarkers for various diseases (15). Existing detection methods fail to predict TIC onset, posing challenges for early diagnosis. The identification of robust biomarkers is crucial for early clinical recognition of TIC patients and improvement of their outcomes. This study employs metabolomics to investigate clinically meaningful biomarkers for early TIC identification. Furthermore, based on differential metabolites, we aim to elucidate significant metabolic pathways involved in TIC pathogenesis.

2 Methods

2.1 Clinical samples

This study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Hunan Provincial People's Hospital, Changsha, China ([2024]-01). Between September 2022 and December 2023, a total of 81 serum samples were collected at the Yuelushan Branch of Hunan Provincial People's Hospital, including 54 from trauma patients (27 with TIC and 27 without TIC) and 27 from healthy controls. All trauma patient samples were obtained within 24 hours after trauma.

This study employed a retrospective case-control design. The case group (TIC group) consisted of all eligible TIC patients (n = 27). To construct a comparable control group, an individual matching strategy was adopted. From trauma patients admitted during the same period who met the inclusion criteria but did not develop TIC, one non-TIC control was matched to each TIC patient based on the following key prognostic factors strongly associated with TIC occurrence: age, diabetes, hypertension, coronary heart disease, and admission scores for trauma patients. However, due to limitations in the available sample pool, matching for sex was not feasible, resulting in a difference in sex distribution between the two groups (Table 1).

Table 1
www.frontiersin.org

Table 1. Cohort characteristics of trauma-induced coagulopathy (TIC) and non-traumatic coagulopathy (Non-TIC) groups.

Diagnostic criteria for TIC (16) (1): Laboratory findings (meeting at least one of the following): PT > 18 s, APTT > 60 s, TT > 15 s, PTr > 1.6; (2) Clinical presentation: Active or potential bleeding requiring blood product transfusion or replacement therapy.

Inclusion criteria: (1) Age 18–80 years; (2) Documented trauma with hemorrhage or blood loss; (3) Time from injury to admission < 24 h; (4) Complete clinical records.

Exclusion criteria: (1) Hematologic disorders or congenital/acquired coagulation abnormalities; (2) Use of corticosteroids, immunosuppressants, or anticoagulants within the past 6 months; (3) Malignancies, severe liver cirrhosis, or other major comorbidities.

The Glasgow Coma Scale (GCS), Trauma Index (TI), Injury Severity Score (ISS) and Abbreviated Injury Scale (AIS) of trauma patients were evaluated upon admission (1719). Scoring was performed in a single-masked manner by two or more physicians, and the average score was used for analysis.

2.2 Sample preparation

Prior to extraction, a mixture of internal standards (10 μL of a cocktail containing stable isotope-labeled compounds such as LPC(17:0)-d5 and PE(17:0/17:0)-d5 at a concentration of 1 μg/mL in methanol) was added to each 100 μL serum aliquot to monitor and correct for variations in sample preparation and instrument analysis. Subsequently, the samples were processed as follows: a 100 μL aliquot, mixed with 400 μL of an 80% methanol aqueous solution, vortexed thoroughly, and incubated on ice for 5 minutes. After centrifugation at 15,000 × g for 20 minutes at 4 °C, the supernatant was collected and diluted with MS-grade water to a final methanol concentration of 53%. A second centrifugation step (15,000 × g, 4 °C, 20 minutes) was performed, and the resulting supernatant was used for LC-MS analysis.

2.3 LC-MS analysis

Metabolite profiling was performed using a Vanquish UHPLC system coupled to an Orbitrap Q Exactive™ HF-X mass spectrometer. Samples were injected onto a Hypersil Gold column (100 × 2.1 mm, 1.9 µm) at a flow rate of 0.2 mL/min. In positive ion mode, mobile phase A consisted of 0.1% formic acid in water, while mobile phase B was methanol; in negative ion mode, mobile phase A comprised 5 mM ammonium acetate (pH 9.0) and mobile phase B was methanol. The gradient elution program was as follows: 2% B (0–1.5 min), 2–85% B (1.5–3 min), 85–100% B (3–10 min), 100–2% B (10–10.1 min), and 2% B (10.1–12 min). MS parameters included: spray voltage = 3.5 kV, capillary temperature = 320 °C, sheath gas flow = 35 psi, auxiliary gas flow = 10 L/min, S-lens RF level = 60, and auxiliary gas heater temperature = 350 °C.

2.4 Data processing and quality assurance

The raw data files were processed using Compound Discoverer (CD) 3.1 software. Metabolic peaks were aligned across samples with tolerances of ≤ 0.2 minutes for retention time (RT) deviation and ≤ 5 ppm for mass deviation. Peak extraction criteria included a mass deviation of 5 ppm, signal intensity deviation of 30%, and a signal-to-noise ratio (S/N) ≥ 3. Ion peaks with more than 50% missing values within sample groups were removed from the dataset. The remaining missing values were imputed with half of the minimum detected value.

Molecular formulas were predicted based on molecular ions and fragment ions, and matched against the mzCloud (https://www.mzcloud.org/), mzVault, and a custom local database. A matching score higher than 36 out of 60 was required for the MS/MS fragment ion spectra to enhance identification confidence.

To ensure data quality and instrumental stability, a pooled quality control (QC) sample was prepared by combining equal aliquots (10 μL) from every individual serum sample in the study. This QC sample was injected repeatedly at the beginning of the analytical sequence to condition the system, and then after every 10 experimental samples throughout the run. The relative standard deviations (RSDs) of the peak areas for the detected metabolites in the QC samples were calculated. Metabolites with an RSD > 30% in the QC samples were considered unstable and were removed from subsequent data analysis, ensuring the reliability of the dataset. The repeated analysis of the same QC sample throughout the run served as a technical replicate to assess the reproducibility of the entire analytical platform.

Subsequent statistical data processing was performed on a Linux (CentOS 6.6) platform using R and Python. Finally, data from positive and negative ion modes were merged into a comprehensive data matrix containing all relevant features extracted from the raw spectra. All subsequent statistical analyses were conducted using this integrated dataset.

2.5 Statistical analysis

Clinical data were analyzed using SPSS 26. Normally distributed continuous variables were expressed as mean ± SD and compared using independent t-tests. Non-normally distributed data were presented as median (P25, P75) and analyzed using nonparametric tests. Categorical variables were reported as counts (%) and compared via chi-square tests. A two-tailed P < 0.05 was considered statistically significant.

To further evaluate the independence of the identified differential metabolites, multivariable logistic regression analysis was performed. The analysis used TIC status as the dependent variable and was adjusted for the confounding factor of sex. After adjusting for this potential confounder, both LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and LysoPE(0:0/18:2(9Z,12Z)) remained significant independent predictors of TIC (P < 0.05). (see Appendix, Tables 1, 2).

2.6 Metabolomics analysis

The data matrix was uploaded to MetaboAnalyst 6.0 (https://www.metaboanalyst.ca/). Data preprocessing on the MetaboAnalyst platform included: normalization by the sum method to adjust for systematic differences; log10 transformation to approximate a normal distribution; and Pareto scaling. Subsequently, t-tests, fold-change (FC) analysis, and volcano plot generation were conducted, along with the construction of PCA, PLS-DA, and OPLS-DA models. The data matrix was then imported into SIMCA14.1 software for a 200-permutation test of the OPLS-DA model. Finally, ROC curve analysis and metabolic pathway analysis were performed on the MetaboAnalyst 6.0 platform. Leave-one-out cross-validation (LOOCV) was performed using R software (version 4.4.3) to obtain an unbiased estimate of the AUC and avoid overoptimism.

3 Results

3.1 Study cohort and patient characteristics

The clinical characteristics of the trauma group (Tra) and control group (H) are presented in Table 2, with no significant differences observed between the groups. The clinical features of the trauma-induced coagulopathy (TIC) and non-TIC cohorts are shown in Table 1. Notably, the TIC group contained significantly fewer male patients. No significant differences were observed in the incidence of other underlying conditions or in the admission assessment scores between the groups.

Table 2
www.frontiersin.org

Table 2. Cohort characteristics of trauma (Tra) and healthy control (H) groups.

3.2 Metabolomics study

A total of 81 serum samples were analyzed, comprising 27 from TIC patients, 27 from non-TIC trauma patients, and 27 from healthy participants. Liquid chromatography–mass spectrometry (LC-MS) detected 7,433 metabolic features and identified 1,479 metabolites.

Preprocessed data were imported into MetaboAnalyst 6.0 for multivariate pattern analysis. Principal component analysis (PCA) revealed small intra-group but significant inter-group differences in serum metabolite profiles among all three groups (Figure 1A). Partial least squares-discriminant analysis (PLS-DA), a supervised method, enhanced inter-group separation while minimizing intra-group variation, demonstrating distinct metabolomic patterns across groups (Figure 1B).

Figure 1
Graphs labeled A to H display score plots and permutation tests for various groups. A to F show two-dimensional plots with different colored groups and ellipses representing confidence intervals, highlighting variations in components. G and H illustrate permutation test results for R2 and Q2 values, indicating statistical significance for different comparisons.

Figure 1. Multivariate statistical analysis of the three groups: (A) Principal Component Analysis (PCA), (B) Partial Least Squares Discriminant Analysis (PLS-DA); Multivariate statistical analysis of the Tra and H groups: (C) Principal Component Analysis (PCA), (D) Partial Least Squares Discriminant Analysis (PLS-DA), (E) Orthogonal Partial Least Squares Discriminant Analysis (OPLS-DA), (G) 200 permutation tests of OPLS-DA; Multivariate statistical analysis of the TIC and Non-TIC groups: (F) Orthogonal Partial Least Squares Discriminant Analysis (OPLS-DA), (H) 200 permutation tests of OPLS-DA.

PCA, PLS-DA, and orthogonal PLS-DA (OPLS-DA) models consistently showed minimal intra-group but marked inter-group differences between the trauma (Tra) and healthy control (H) groups (Figures 1C–E). OPLS-DA further improved model specificity by filtering out classification-unrelated noise. Application of OPLS-DA to the TIC versus non-TIC comparison revealed clear separation between the groups, indicating significant metabolomic disparities (Figure 1F). A 200-permutation test confirmed the reliability of all OPLS-DA models (Figures 1G, H). Metabolites with variable importance in projection (VIP) scores >1.5 were selected as biologically significant.

3.3 Analysis of differential metabolites

Student's t-test and fold-change (FC) analysis were performed on the metabolic data of the Tra and H groups, as well as the TIC and non-TIC groups, and volcano plots were generated. Differential metabolites between the two groups were identified. Based on FC > 2 and adjusted P < 0.05, 496 differential metabolites were screened in the Tra and H groups, of which 144 were significantly down-regulated and 352 were significantly up-regulated (Figure 2A). In the TIC and non-TIC groups, 212 differential metabolites were identified, of which 7 were significantly down-regulated and 205 were significantly up-regulated (Figure 2B).

Figure 2
Two panels labeled A and B show volcano plots. Panel A displays data points representing log2 fold change against negative log10 p-value, with colors ranging from purple to red, indicating varying fold changes. Panel B also presents similar data with a different scale, emphasizing distinct data point distributions. Both plots include legends for color and p-value significance.

Figure 2. The volcano plot of the Tra and H groups (A), and the volcano plot of the TIC and Non-TIC groups (B); The X-axis corresponds to log 2 (FC), and the Y-axis corresponds to -log 10 (P value).

By integrating multivariate and univariate analyses (FC > 2.0, adjusted P < 0.05, and VIP > 1.5), we identified 186 differential metabolites between the Tra and H groups (69 down-regulated, 117 up-regulated). For the TIC versus non-TIC comparison, 74 metabolites were selected (2 down-regulated, 72 up-regulated).

3.4 Metabolic pathway analysis

Next, receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic efficacy of the differential metabolites. The area under the curve (AUC) values of the 74 differential metabolites for distinguishing the TIC and non-TIC groups were all ≥ 0.7 (see appendix). To elucidate the potential association between these differential metabolites and traumatic coagulopathy, metabolic pathway analysis was performed on the screened differential metabolites using MetaboAnalyst 6.0. Using -log10(p) > 1.0 and Pathway Impact > 0.10 as screening criteria, the analysis revealed that the key metabolic pathways differentiating the TIC and non-TIC groups were inositol phosphate metabolism and alanine, aspartate, and glutamate metabolism (Figure 3).

Figure 3
Scatter plot with pathway impact on the x-axis and negative logarithm of p-value on the y-axis. Various colored circles represent data points, with color intensity and size varying. A cluster of larger circles appears near the top center, suggesting significant pathway impacts.

Figure 3. Metabolic pathway analysis diagram of the TIC and Non-TIC groups; The X-axis corresponds to Pathway Impact, and the Y-axis corresponds to -log 10 (P). The large circle on the left represents Alanine, aspartate and glutamate metabolism, and the one on the right represents Inositol phosphate metabolism.

3.5 Screening of potential biomarkers

Using FC > 2.0, adjusted P < 0.05, VIP > 1.5, and AUC > 0.9 as criteria, five differential metabolites were identified as candidate biomarkers for distinguishing TIC from non-TIC. After excluding three exogenous compounds, two metabolites—LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and LysoPE(0:0/18:2(9Z,12Z))—were selected as potential biomarkers with high diagnostic value (Figures 4A, B).

Figure 4
Panel A shows a Receiver Operating Characteristic (ROC) curve for LysoPE(20:4) with an AUC of 0.933, and a box plot comparing two groups with distinct red and green boxes. Panel B presents a ROC curve for LysoPE(0:0/18:2) with an AUC of 0.916, along with a similar box plot.

Figure 4. ROC curves of potential biomarkers in the TIC and Non-TIC groups (A) LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0), (B) LysoPE(0:0/18:2(9Z,12Z)); Through the ROC curve analysis, the diagnostic value of different metabolites in disease screening was evaluated. The x-axis and y-axis represent the false positive rate and the true positive rate respectively.

LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) yielded an AUC of 0.933 (95% CI: 0.849–0.995), with a sensitivity and specificity of 0.963 and 0.815, respectively. LysoPE(0:0/18:2(9Z,12Z)) had an AUC of 0.916 (95% CI: 0.818–0.914), with a sensitivity and specificity of 0.963 and 0.852, respectively (Table 3). Leave-one-out cross-validation yielded the following AUC values: 0.929 for LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and 0.885 for LysoPE(0:0/18:2(9Z,12Z)) (see Appendix; Figure 1).

Table 3
www.frontiersin.org

Table 3. Performance characteristics of biomarkers and routine coagulation indicators.

Current admission assessments for trauma patients—including the Trauma Index (TI), Glasgow Coma Scale (GCS), Injury Severity Score (ISS), and Abbreviated Injury Scale (AIS)—showed no statistical differences between TIC and non-TIC groups (P > 0.05) and demonstrated poor ability to predict TIC onset. Presently, the clinical diagnosis of TIC primarily relies on conventional coagulation tests such as PT and APTT.

Although PT, PT-INR, PT%, and APTT values differed significantly between TIC and non-TIC groups (P < 0.05), their diagnostic efficacy was markedly inferior to that of the two differential metabolites identified in this study (Figure 5). These metabolites demonstrated significantly higher sensitivity and specificity, as well as superior overall diagnostic performance (Table 3). The AUC values of LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and LysoPE(0:0/18:2(9Z,12Z)) were significantly higher than those of PT, PT-INR, PT%, and APTT (P < 0.05).

Figure 5
ROC curve displaying sensitivity versus 1-specificity for PT, PT-INR, PT%, and APTT with AUC values of 0.724, 0.725, 0.272, and 0.663, respectively. A yellow diagonal line marks the reference.

Figure 5. ROC curve of routine coagulation indexes; The X-axis corresponds to 1 - Specificity, and the Y-axis corresponds to Sensitivity.

4 Discussion

Although conventional coagulation tests such as prothrombin time (PT) and activated partial thromboplastin time (APTT) remain the standard diagnostic tools for trauma-induced coagulopathy (TIC) in clinical practice, these parameters lack the sensitivity to detect early TIC onset and fail to reliably predict subsequent coagulation deterioration. Consequently, there is an urgent unmet need for clinical methods to anticipate TIC development in trauma patients. This study aimed to address this need by employing a metabolomics approach to identify highly sensitive and specific biomarkers for the early diagnosis of TIC.

Our metabolomic analysis revealed that two lysophosphatidylethanolamines—LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and LysoPE(0:0/18:2(9Z,12Z))—serve as novel biomarkers for TIC, demonstrating exceptional diagnostic performance in early-stage detection. These metabolites represent novel quantitative indicators for TIC diagnosis, exhibiting significantly higher AUC values, sensitivity, and specificity than traditional coagulation assays (Table 3). The superior diagnostic accuracy of these biomarkers (AUC > 0.9) highlights their potential to revolutionize clinical practice by enabling earlier intervention, thereby improving patient outcomes. Furthermore, their high predictive capacity suggests utility in identifying at-risk patients prior to overt coagulopathy, facilitating preemptive therapeutic strategies. Additionally, inositol phosphate metabolism and alanine, aspartate, and glutamate metabolism were also associated with TIC.

LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) and LysoPE(0:0/18:2(9Z,12Z)) are two distinct lysophosphatidylethanolamine (LysoPE) molecules. Previous research on increased cardiovascular disease risk under hypercholesterolemia found that the oxidized low-density lipoprotein (ox-LDL) group exhibited significantly reduced APTT and PT, along with markedly elevated levels of seven LysoPEs in plasma metabolomics (20). Our observation of significantly elevated serum levels of LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) in patients with TIC provides direct evidence that an increase in LysoPE can lead to coagulation dysfunction, supporting its potential role in inducing a hypercoagulable state. However, LysoPE(0:0/18:2(9Z,12Z)) levels were significantly decreased in TIC patient plasma, a phenomenon that may be closely related to the essential pathogenesis of TIC. Specifically, TIC progresses acutely through different pathological stages (hypocoagulability, hyperfibrinolysis, and hypercoagulability), with transitions from hypocoagulability to hypercoagulability potentially occurring within minutes or hours. These dynamic changes result in a unique mixed bleeding-thrombosis phenotype throughout the disease course.

Metabolic pathway analysis revealed that the most significantly altered pathways in the TIC/non-TIC comparison were inositol phosphate metabolism and alanine, aspartate, and glutamate metabolism. The inositol phosphate metabolism pathway generates inositol trisphosphate (IP3), which promotes the release of intracellular stored calcium ions and extracellular calcium influx, thereby increasing intracellular calcium concentration (21, 22). Elevated calcium levels lead to platelet activation and thrombosis. Existing studies have demonstrated the important role of inositol phosphate metabolism in TIC pathogenesis (22, 23). Our results further support the upregulation of this pathway in TIC patients, suggesting its contribution to coagulation abnormalities.

This was a single-center study with a relatively limited sample size, which may restrict the generalizability and statistical power of our findings. For example, the gender difference between the TIC and non-TIC groups likely resulted from the limited sample size. Although the case-control design employed is suitable for early-stage biomarker screening, it does not allow for an accurate evaluation of the biomarker's actual diagnostic performance (such as sensitivity and specificity) in a consecutively enrolled prospective cohort. This study preliminarily identified two LysoPEs that exhibited significant alterations in TIC patients and demonstrated promising discriminative ability in the initial cohort. External validation in larger independent cohorts is warranted to further assess their diagnostic efficacy and clinical applicability. Furthermore, all samples were collected within 24 hours post-trauma, which precluded continuous dynamic monitoring of the two LysoPEs and thus limited the exploration of their time-dependent dynamic characteristics. The specific pathways and functional roles by which the differential metabolites influence coagulation function will need to be further validated through in vitro experiments, such as endothelial cell models. Moreover, the cost of detecting the two LysoPEs currently offers no distinct advantage over routine coagulation tests, and their use has not been widely adopted in clinical practice.

By exploring metabolomic differences in emergency severe trauma patients with or without TIC, we identified potential biomarkers with excellent diagnostic value and determined key metabolic pathways that may play crucial roles in disease progression. This study lays the foundation for future research, which should involve multiple clinical centers and larger sample sizes, as well as integrate transcriptomics and proteomics to elucidate the metabolic-coagulation interaction network and conduct more in-depth investigations.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Committee of Hunan Provincial People’s Hospital in Changsha, Hunan Province, China ([2024]-01). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

XX: Software, Formal Analysis, Writing – original draft, Investigation, Visualization, Data curation, Methodology, Validation, Conceptualization. QO: Investigation, Conceptualization, Software, Writing – original draft, Formal Analysis, Methodology. MS: Writing – review & editing, Data curation. QY: Data curation, Investigation, Writing – review & editing. SL: Writing – review & editing, Investigation. YH: Writing – review & editing. JW: Resources, Data curation, Writing – review & editing, Formal Analysis. XZ: Supervision, Conceptualization, Writing – review & editing, Resources, Data curation. CT: Formal Analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing, Conceptualization, Data curation.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was funded by grants from the Key Cultivation Project of Renshu Fund of Hunan Provincial People's Hospital (RS2022A15), the Health Commission high-level talent project of Hunan Provincial (R2023087), the Natural Science Foundation of Hunan Province (2023JJ10024), the Hunan Provincial People's Hospital Innovation Team (KCT202402), the Hunan Provincial Geriatric Research Institute (20240307-1014), the Key Research and Development Program of Hunan Province (2025JK2130), the General Project of Hunan Provincial Health Commission (20254678), the Excellent Youth Project of Hunan Provincial Department of Education (24B0073), the Changsha Natural Science Foundation (kq2502185) and the Graduate Student Research Innovation Project of Hunan Provincial Department of Education (CX20250819).

Acknowledgments

We would like to sincerely thank the First Affiliated Hospital of Hunan Normal University (Hunan Provincial People's Hospital) for its valuable suggestions. This research was approved by the Clinical Experimental Research Institute of the First Affiliated Hospital of Hunan Normal University (Hunan Provincial People's Hospital).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1632694/full#supplementary-material

Abbreviations

TIC, Trauma-induced coagulopathy; PCA, principal component analysis; OPLS-DA, orthogonal partial least-squares discrimination analysis; PLS-DA, partial least-squares discrimination analysis; VIP, variable importance in projection; FC, fold change; ROC, receiver operating characteristic; AUC, area under the curve; LysoPE, lysophosphatidylethanolamine.

References

1. Norton R and Kobusingye O. Injuries. N Engl J Med. (2013) 368:1723–30. doi: 10.1056/NEJMra1109343

PubMed Abstract | Crossref Full Text | Google Scholar

2. Kalkwarf KJ, Drake SA, Yang Y, Thetford C, Myers L, Brock M, et al. Bleeding to death in a big city: An analysis of all trauma deaths from hemorrhage in a metropolitan area during 1 year. J Trauma Acute Care Surg. (2020) 89:716–22. doi: 10.1097/TA.0000000000002833

PubMed Abstract | Crossref Full Text | Google Scholar

3. Moore HB, Moore EE, Chapman MP, McVaney K, Bryskiewicz G, Blechar R, et al. Plasma-first resuscitation to treat haemorrhagic shock during emergency ground transportation in an urban area: a randomised trial. Lancet. (2018) 392:283–91. doi: 10.1016/S0140-6736(18)31553-8

PubMed Abstract | Crossref Full Text | Google Scholar

4. Brohi K, Singh J, Heron M, and Coats T. Acute traumatic coagulopathy. J Trauma. (2003) 54:1127–30. doi: 10.1097/01.TA.0000069184.82147.06

PubMed Abstract | Crossref Full Text | Google Scholar

5. Brenner M, Stein DM, Hu PF, Aarabi B, Sheth K, and Scalea TM. Traditional systolic blood pressure targets underestimate hypotension-induced secondary brain injury. J Trauma Acute Care Surg. (2012) 72:1135–9. doi: 10.1097/TA.0b013e31824af90b

PubMed Abstract | Crossref Full Text | Google Scholar

6. Moore EE, Moore HB, Kornblith LZ, Neal MD, Hoffman M, Mutch NJ, et al. Trauma-induced coagulopathy. Nat Rev Dis Primers. (2021) 7:30. doi: 10.1038/s41572-021-00264-3

PubMed Abstract | Crossref Full Text | Google Scholar

7. Cohen MJ and Christie SA. New understandings of post injury coagulation and resuscitation. Int J Surg. (2016) 33:242–5. doi: 10.1016/j.ijsu.2016.05.037

PubMed Abstract | Crossref Full Text | Google Scholar

8. Kornblith LZ, Moore HB, and Cohen MJ. Trauma-induced coagulopathy: the past, present, and future. J Thromb Haemost. (2019) 17:852–62. doi: 10.1111/jth.14450

PubMed Abstract | Crossref Full Text | Google Scholar

9. Brohi K, Cohen MJ, Ganter MT, Schultz MJ, Levi M, Mackersie RC, et al. Acute coagulopathy of trauma: hypoperfusion induces systemic anticoagulation and hyperfibrinolysis. J Trauma. (2008) 64:1211–7; discussion 1217. doi: 10.1097/TA.0b013e318169cd3c

PubMed Abstract | Crossref Full Text | Google Scholar

10. Petros S. Trauma-induced coagulopathy. Hamostaseologie. (2019) 39:20–7. doi: 10.1055/s-0039-1677853

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ganter MT and Pittet JF. New insights into acute coagulopathy in trauma patients. Best Pract Res Clin Anaesthesiol. (2010) 24:15–25. doi: 10.1016/j.bpa.2009.09.010

PubMed Abstract | Crossref Full Text | Google Scholar

12. Ho VK, Wong J, Martinez A, and Winearls J. Trauma-induced coagulopathy: Mechanisms and clinical management. Ann Acad Med Singap. (2022) 51:40–8. doi: 10.47102/annals-acadmedsg.2020381

PubMed Abstract | Crossref Full Text | Google Scholar

13. Cralley AL, Erickson C, Schaid TR, Hallas W, Thielen O, Mitra S, et al. The proteomic and metabolomic signatures of isolated and polytrauma traumatic brain injury. Am J Surg. (2023) 226:790–7. doi: 10.1016/j.amjsurg.2023.07.040

PubMed Abstract | Crossref Full Text | Google Scholar

14. Coleman JR, D’Alessandro A, LaCroix I, Dzieciatkowska M, Lutz P, Mitra S, et al. A metabolomic and proteomic analysis of pathologic hypercoagulability in traumatic brain injury patients after dura violation. J Trauma Acute Care Surg. (2023) 95:925–34. doi: 10.1097/TA.0000000000004019

PubMed Abstract | Crossref Full Text | Google Scholar

15. Bujak R, Struck-Lewicka W, Markuszewski MJ, and Kaliszan R. Metabolomics for laboratory diagnostics. J Pharm BioMed Anal. (2015) 113:108–20. doi: 10.1016/j.jpba.2014.12.017

PubMed Abstract | Crossref Full Text | Google Scholar

16. Yue MX, Liang HP, and Du DY. Expert consensus on the diagnosis and health emergency management of acute traumatic coagulopathy and coagulopathy. Chinese Journal of Health Emergency (Electronic Edition). (2016) 2(4):197–203. Available online at: https://kns.cnki.net/kcms2/article/abstractv=MsI1aVwi7Tt_hyQh0vcCO5ve9TMpgUjz4JNGD31ZVRWSAI7SWnYlafuYporvqFzfzP0gNKU9dYim7UVqJbi0nbpizMG4ABybY5pK8omZQnxWI9t-7Lt2YqGOfD2ho_lYWZKtE4WwtAiy8f1x01UUBuzFkr7cOAbolIBS2uRgfTQCNl_Z8dsJw==&uniplatform=NZKPT&language=CHS.

Google Scholar

17. Mehta R, trainee GP, and Chinthapalli K. consultant neurologist. Glasgow coma scale explained. BMJ. (2019) 365:l1296. doi: 10.1136/bmj.l1296

PubMed Abstract | Crossref Full Text | Google Scholar

18. Teasdale G and Jennett B. Assessment of coma and impaired consciousness. A Pract Scale Lancet. (1974) 2:81–4. doi: 10.1016/S0140-6736(74)91639-0

PubMed Abstract | Crossref Full Text | Google Scholar

19. Baker SP, O’Neill B, Haddon W, and Long WB. The injury severity score: a method for describing patients with multiple injuries and evaluating emergency care. J Trauma. (1974) 14:187–96. doi: 10.1097/00005373-197403000-00001

PubMed Abstract | Crossref Full Text | Google Scholar

20. Kim M, Yoo HJ, Lee D, and Lee JH. Oxidized LDL induces procoagulant profiles by increasing lysophosphatidylcholine levels, lysophosphatidylethanolamine levels, and Lp-PLA2 activity in borderline hypercholesterolemia. Nutrit Metab Cardiovasc Dis. (2020) 30:1137–46. doi: 10.1016/j.numecd.2020.03.015

PubMed Abstract | Crossref Full Text | Google Scholar

21. Berridge MJ and Irvine RF. Inositol phosphates and cell signalling. Nature. (1989) 341:197–205. doi: 10.1038/341197a0

PubMed Abstract | Crossref Full Text | Google Scholar

22. Varga-Szabo D, Braun A, and Nieswandt B. Calcium signaling in platelets. J Thromb Haemost. (2009) 7:1057–66. doi: 10.1111/j.1538-7836.2009.03455.x

PubMed Abstract | Crossref Full Text | Google Scholar

23. Grosse J, Braun A, Varga-Szabo D, Beyersdorf N, Schneider B, Zeitlmann L, et al. An EF hand mutation in Stim1 causes premature platelet activation and bleeding in mice. J Clin Invest. (2007) 117:3540–50. doi: 10.1172/JCI32312

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: liquid chromatography-tandem mass spectrometry (LC-MS), lysophosphatidylethanolamine (LysoPE), metabolomics, trauma, trauma-induced coagulopathy (TIC)

Citation: Xu X, Ouyang Q, Shao M, Yi Q, Liu S, Huang Y, Wang J, Zhang X and Tan C (2025) Screening of biomarkers for early diagnosis of trauma-induced coagulopathy based on untargeted metabolomics. Front. Endocrinol. 16:1632694. doi: 10.3389/fendo.2025.1632694

Received: 21 May 2025; Accepted: 30 September 2025;
Published: 15 October 2025.

Edited by:

Gaetano Santulli, Albert Einstein College of Medicine, United States

Reviewed by:

Edward Wargent, University of Buckingham, United Kingdom
Bagus Gita Pranata Putra, Universitas Warmadewa, Indonesia

Copyright © 2025 Xu, Ouyang, Shao, Yi, Liu, Huang, Wang, Zhang and Tan. 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: Chaochao Tan, dzIyMzZAaHVubnUuZWR1LmNu; Jia Wang, d2FuZ2ppYUBodW5udS5lZHUuY24=; Xingwen Zhang, emhhbmd4aW5nd2VuQDE2My5jb20=

These authors have contributed equally to this work

Disclaimer: 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.