Urinary peptides provide information about the risk of mortality across a spectrum of diseases and scenarios

Background There is evidence of pre-established vulnerability in individuals that increases the risk of their progression to severe disease or death, although the mechanisms causing this are still not fully understood. Previous research has demonstrated that a urinary peptide classifier (COV50) predicts disease progression and death from SARS-CoV-2 at an early stage, indicating that the outcome prediction may be partly due to vulnerabilities that are already present. The aim of this study is to examine the ability of COV50 to predict future non-COVID-19-related mortality, and evaluate whether the pre-established vulnerability can be generic and explained on a molecular level by urinary peptides. Methods Urinary proteomic data from 9193 patients (1719 patients sampled at intensive care unit (ICU) admission and 7474 patients with other diseases (non-ICU)) were extracted from the Human Urinary Proteome Database. The previously developed COV50 classifier, a urinary proteomics biomarker panel consisting of 50 peptides, was applied to all datasets. The association of COV50 scoring with mortality was evaluated. Results In the ICU group, an increase in the COV50 score of one unit resulted in a 20% higher relative risk of death [adjusted HR 1.2 (95% CI 1.17–1.24)]. The same increase in COV50 in non-ICU patients resulted in a higher relative risk of 61% [adjusted HR 1.61 (95% CI 1.47–1.76)], consistent with adjusted meta-analytic HR estimate of 1.55 [95% CI 1.39–1.73]. The most notable and significant changes associated with future fatal events were reductions of specific collagen fragments, most of collagen alpha I (I). Conclusion The COV50 classifier is predictive of death in the absence of SARS-CoV-2 infection, suggesting that it detects pre-existing vulnerability. This prediction is mainly based on collagen fragments, possibly reflecting disturbances in the integrity of the extracellular matrix. These data may serve as a basis for proteomics-guided intervention aiming towards manipulating/ improving collagen turnover, thereby reducing the risk of death. Supplementary Information The online version contains supplementary material available at 10.1186/s12967-023-04508-6.


Background
Pre-existing vulnerabilities play a key role in determining an individual's risk for disease progression or death [1], highlighting the importance of considering these factors when managing diseases.Given the complexity of the disease-associated molecular mechanisms [2][3][4][5] and the factors impacting the outcome, recognizing and understanding the pre-established vulnerabilities can help identify high-risk individuals and tailor treatment strategies, ultimately improving outcome.
Among theories attempting to explain this phenomenon (e.g., in the context of trauma) [6], the "two-hit" model has evolved.According to this scenario, the stress response encompasses the physiological reaction to the initial injury (referred to as the "first hit"), followed by a reaction to the secondary insult/ intervention (known as the "second hit") [6,7].The model is rooted in the fundamental idea that consecutive insults, which may not have significant effects individually, can result in a profound physiological response.This reaction can manifest in various biological systems and can be evaluated by measuring multiple parameters [8].However, the molecular mechanisms responsible for the "two-hit" model are complex and not fully understood [6,8].In general, the "first hit" acts as a priming event that predisposes the patient to develop a systemic inflammatory syndrome, with a key feature being a leak of the endothelium.This initially manifests in a specific body region, but eventually affects multiple organs.Subsequently, a second insult can trigger an exaggerated inflammatory response, responsible for potentially life-threatening conditions such as multiple organ failure and multiple organ dysfunction syndrome [6,8].Understanding the vulnerability to the "second hit" can support minimising the impact of complications, potentially leading to a better outcome.
Recently, it has been suggested that SARS-CoV-2 infection could act as a "second hit".SARS-CoV-2 is among the main conditions associated with collapsing glomerulopathy, acting as a "second hit" in susceptible patients with APOL1 risk alleles, similar to human immunodeficiency virus and other viruses [9].Another example involves complement-mediated disorder, which seems to be a predominant form of thrombotic microangiopathy associated with COVID-19.Considering the development of thrombotic microangiopathy following SARS-CoV-2 infection, it was suggested that the virus acted as a "secondary trigger", revealing an underlying complement defect [10].
Previous biomarker research demonstrated the capability of a urinary peptide-based classifier (COV50) to predict disease progression and death from SARS-CoV-2 at the earliest possible date, i.e., upon the first positive indication of a SARS-CoV-2 infection [11,12].This assessment was based on the measurement of 50 specific urinary peptides, with the most prominent changes involving the reduction of peptides derived from collagen alpha 1(I), polymeric immunoglobulin receptor and CD99 antigen, and an increase in peptides derived from alpha-1-antitrypsin [13].The ability to predict outcome very soon after infection suggests that the prediction may not be solely based on molecular events associated with SARS-CoV-2 infection, but at least in part due to pre-established vulnerability, resulting from a "first hit".This would indicate that prediction of severe disease course may be feasible even before the infection.We hypothesised that this pre-established vulnerability can be generic (expanded to other indications) and COV50 could serve as a biomarker for detecting vulnerable subjects who are adversely affected by other clinical insults.
Hence, the present study aimed to examine the ability of COV50 to predict future non-COVID-19-related mortality in patients admitted to the intensive care unit (ICU) or having other diseases (non-ICU).If the hypothesis is confirmed, a significantly higher number of individuals in this vulnerable population (defined by a high COV50 score) should experience death compared to the population with a lower score.

Patients
ICU: Patients from the medical, surgical, or mixed ICUs at 14 university hospitals from the FROG-ICU study were included [14].Inclusion criteria were mechanical ventilation or administration of vasoactive agents for at least 24 h.Exclusion criteria were age under 18, severe head injury with a Glasgow Coma Scale below 8, brain death or persistent vegetative state, pregnancy or breastfeeding, transplantation in the past 12 months, moribund status, and lack of social security coverage.All capillary electrophoresis coupled to mass spectrometry (CE-MS) datasets with a 1-year follow-up and information on relevant co-variables (age, body mass index (BMI), sex, blood pressure, estimated glomerular filtration rate (eGFR), presence of diabetes, kidney, cardiovascular disease, hypertension) were included in the present study without pre-selection.
Non-ICU: The assessment of COV50 in the non-ICU population was based on 7474 datasets from the Human Urinary Proteome Database [15,16] with available information on age, sex, eGFR, blood pressure, BMI, presence of diabetes, kidney disease, cardiovascular disease, hypertension, and a follow-up data.
All datasets were obtained from previously published studies and fully anonymized.Ethical review and approval were waived for this study by the ethics committee of the Hannover Medical School, Germany (no. 3116-2016), as all data were fully anonymized.The number of subjects per study and patient characteristics are listed in Table 1 and Additional file 1.

Urinary proteome/peptidome data
The urinary proteome is well characterized and reference standards are available [17].Urinary proteome analysis was conducted on urine samples collected at study inclusion and subsequently bio-banked until assayed.Detailed information on urine sample preparation, proteome analysis by CE-MS, data processing, and sequencing of the urinary peptides allowing for the identification of parental proteins is available in previous publications [11,[18][19][20] and described in detail in Additional file 2.

Statistics
As descriptive statistics for the ICU and non-ICU samples, shown in Table 1, median and 1 st and 3 rd quartile (IQR) were used for continuous variables and absolute (N) and relative frequencies (%) for categorical variables.
Hypotheses of no differences in scale or distribution of patient characteristics between the death and non-death groups were tested with Wilcoxon-Mann-Whitney tests for continuous and with χ 2 -homogeneity tests for categorical variables.Kernel density estimates of the distribution of COV50 scores divided by ICU and mortality groups are depicted in Fig. 1A.Mortality per person-time, stratified by age and COV50 groups, as shown in Fig. 1B, was estimated as the ratio of the number of the deceased to the sum of all patients` observation times within each group, scaled to 100 person-years.Corresponding mortality probabilities with their 95% confidence intervals (CI) for each group, presented in Fig. 1C, were estimated through a logistic regression involving all 9193 patients.
For each study, separate unadjusted Cox regressions were conducted to analyse the effect of the COV50 score on experiencing death, as listed in Fig. 2A.In Fig. 2B, these models were additionally adjusted for age, female, log(BMI), mean arterial pressure (MAP) and log(eGFR).All regressors besides female and COV50 were normalized (mean 0, sd 1).The natural logarithms of the estimated hazard ratios (logHR) and their standard errors were combined in meta-analyses to determine the effect of the COV50 score on mortality.A random effects model was estimated after the assumption that all included studies are heterogeneous, i.e., coming from different populations.Study weights are based on the logHR estimates` uncertainty, namely their standard errors.Studies were categorized into more homogenous subgroups, and estimates for each subgroup were displayed in Fig. 2. Overall and group-wise between-study heterogeneity is presented with τ 2 and assessed by Higgins & Thompson's I 2 statistic.χ 2 -Tests for heterogeneity and subgroup differences were based on Cochran's Q. Random effects meta-analysis estimates were presented with 95% CIs and a 95% prediction interval for the overall effect.One Cox regression, stratified by study pooling all 9193 patients, was used as a benchmark to the meta-analytic approach.As displayed in Table 2, the model`s adjustment specification matched to the adjusted separate study regressions (Fig. 2B).To be comparable to the adjusted meta-analytic estimate, HRs for COV50 interacted with ICU and non-ICU, as well as for the above-mentioned non-ICU subgroups, were estimated.Standard errors were clustered on the study level for more robust inference and due to unobserved heterogeneity between studies.The models log-likelihood, associated Wald test and concordance were reported in Table 2.We allowed for a type 1 error of 5%, all hypotheses were two-sided.All analyses were carried out using R 4.2.2.

Results
First, we assessed the hypothesis that the COV50 classifier defines a vulnerable population at the molecular level, irrespective of SARS-CoV-2 infection.For that purpose, we examined datasets from subjects from the FROG-ICU study [21], as this study was more comparable to the CRIT-COV study (patients in ICU) and had available a large number of endpoints.We identified 1719 datasets to be included in this study, for which follow-up and information on relevant co-variables were available [38].
To further support our analysis, we also investigated whether COV50 could predict mortality in subjects outside the ICU.Studies with more than 50 individuals and available follow-up and information on relevant co-variables were selected from the Human Urinary Proteome Database [15,16].
Demographic information on the subjects included in the study is presented in Table 1, separated into ICU and non-ICU groups, as well as by death status.More detailed information is provided in Additional file 1.Among the risk factors for death, we found significant differences at the aggregate level in both, the ICU and the non-ICU subgroups, as expected.The median COV50 score is significantly higher (p < 0.001) in patients who experienced death during the observation period, as also displayed in its distributions in Fig. 1A.
Considering that age is a crucial risk factor for mortality, we investigated the relationship between COV50 and mortality across different age groups.The results are illustrated in Fig. 1B, C. Panel B depicts mortality in person-time in COV50 groups, whereas panel C relates mortality as a percentage along the continuum of COV50 scores.In both subgroups, an increase in COV50 accompanies higher mortality, with this effect being more pronounced among older individuals.
The crude HRs in Fig. 2A for all studies generally indicate an association of a higher relative risk of death with increasing COV50 scores, with all but 5 studies showing a significantly elevated relative risk.In the meta-analysis, the combined HR estimate for all subgroups, except the diabetes-related studies, significantly differs from 1, as evident from the 95% CIs.Overall, the adjustment for risk factors reduced the COV50 HR estimates, in line with expectations, as adjustment typically improves comparability by accounting for observed between study heterogeneity at the patient level.However, in studies with low numbers of events (particularly PersTIgAN, STOPI-gAN), variance increased substantially with the adjustment.The estimates from the meta-analysis resulted in an unadjusted HR of 1.77 [95% CI 1.58-1.97]and an adjusted HR of 1.55 [95% CI 1.39-1.73].Although displaying a trend within and between the subgroups, neither heterogeneity nor subgroup differences were found to be statistically significant (Fig. 2).
The subgroup HR estimates derived from the adjusted meta-analysis in Fig. 2B are robust, as they closely align with the corresponding estimates from the pooled adjusted Cox regression shown in Table 2. Within the ICU group, an increase in the COV50 score of one unit results (on average) in a 20% higher relative risk of death [adj.HR 1.2 (95% CI 1.17-1.24)].As the absolute risk of death is considerably lower in non-ICU patients, the same increase in COV50 in non-ICU patients results in a higher relative risk of 61% [adj.HR 1.61 (95% CI 1.47-1.76)].These findings align well with the adjusted HR estimate of 1.55 [95% CI 1.39-1.73]obtained from the meta-analysis.
COV50 is a composite score based on 50 distinct urinary peptides.To examine which of these 50 peptides served as individual predictors of death in the cohorts investigated (ICU, non-ICU), we compared the distribution of the 50 peptides in the datasets from survivors with those from subjects that died.The results of this   List of 50 urinary peptides included in the COV50 classifier and their respective regulation trend in investigated cohorts (ICU and non-ICU), and corresponding HRs for predicting mortality.P-values below 0.05 are marked in bold.Peptides with an increased abundance in the case vs. the control group (> 1.1) are marked in red, while those with decreased abundance (< 0.91) are marked in green.The regulation trend was calculated by dividing average abundances in the individual case vs. the control group BH Benjamini-Hohberg, CI confidence interval, HR hazard ratio, ICU intensive care unit analysis are shown in Table 3.A high degree of concordance was observed when comparing the peptides regulation trend in the context of COVID-19, death in or after ICU, or death without ICU stay.The association of single peptides with the mortality was also supported by the Cox regression analysis.The most notable and significant changes associated with future fatal events are the reductions in specific collagen fragments, with most of them derived from collagen alpha I(I).
In the ICU subjects, 33 out of the 50 peptides were found to be significantly associated with future death, with 29 having a regulation trend > 1.1 or < 0.91.Among the latter, 28 exhibited a regulation trend in a similar direction as for critical/lethal COVID-19.However, one peptide from alpha-1-acid glycoprotein 1 had an opposing regulation.Upon investigating the most prominent peptides derived from collagen, all significant changes were concordant between death in COVID-19 or ICU.In the non-ICU subjects, 34 of the 50 peptides were significantly associated with future death, with 29 having a regulation trend > 1.1 or < 0.91.Of these 29 peptides, 22 showed a regulation concordant with the one in critical/ lethal COVID-19, while 7 peptides changed in an opposing direction.The latter group includes peptides derived from polymeric immunoglobulin receptor, neurosecretory protein VGF, keratin, type II, hornerin, collagen alpha-1(I), and collagen alpha-1(III).The most prominent difference in comparison to the distribution in COVID-19 patients was observed for peptides derived from CD99 antigen and polymeric immunoglobulin receptor.While a consistent trend and significant reduction of multiple CD99 antigen and polymeric immunoglobulin receptor peptides was associated with severe disease and mortality in critical COVID-19 patients, this distribution was less evident or not observed for polymeric immunoglobulin receptor in the non-ICU population and for CD99 antigen in all cohorts not infected with SARS-CoV-2.

Discussion
This study is the first to investigate a peptide-based classifier, COV50, and specific urinary peptides in a large and diverse population of patients both inside and outside the ICU.The data demonstrate that COV50 not only predicts an unfavourable outcome of a COVID-19 episode but also appears to identify "vulnerable subjects" who are likely at substantially higher risk of severe or lethal COVID-19.This vulnerability also seems relevant in other clinical situations (e.g., non-SARS-CoV-2 infections), including those that lead to ICU admission, thereby increasing the risk of death in various pathological conditions.This is likely related to the fact that urinary peptides reflect local and systemic changes.It has been proposed that approximately 70% of urinary proteins under normal physiological circumstances are derived from the kidney and the urinary tract.The remaining 30% originates from other organs and is released into the bloodstream [39].Although some of the peptides contained in COV50 have been previously identified in plasma [40,41], the origin of specific naturally occurring urinary peptides cannot be predicted with high certainty.The most prominent and consistent findings are the reduction of several specific urinary collagen fragments, most from collagen alpha-I(I).This decrease in collagen fragments may indicate reduced collagen degradation within the extracellular matrix, which is expected to result in increased fibrosis.Fibrosis has been associated with various diseases affecting different organs, including the liver, kidney, lungs, and heart [42].Previous studies have demonstrated an association between fibrosis and poor outcome in patients with various pathologies [43][44][45].Fibrosis may constitute the "first hit" and induce vulnerability to "second hit" events either in e.g., infectious or general (cardiovascular) scenarios.In this context, a pre-existing fibrotic condition may render an organ/ tissue more vulnerable to further damage or insults from a second event or trigger.Fibrosis alters the normal structure and function of the affected tissue, compromising its capacity to respond and recover from subsequent insults.Consequently, when a second hit, such as infection or inflammation occurs, it can lead to more severe complications and worsen the overall outcome.
The concordance of significant changes in individual peptides observed due to critical/lethal COVID-19 appears to be higher in the context of ICU than in non-ICU subjects.While an objective measure to assess significant differences does not seem to exist, a concordance (based on up-or down-regulation) of 97% (in the case of ICU) compared to 76% (in the case of non-ICU) is at least indicative.
As expected, there are similarities in changes in biomarkers in patients developing the critical condition, irrespective of the underlying pathology and disease aetiology.At the same time, it becomes evident that specific changes, a decrease of peptides from CD99 antigen, are associated more specifically with critical COVID-19, and cannot be consistently associated with all-cause death, neither in nor outside ICU.This suggests that the "second hit" in the context of a SARS-CoV-2 infection is depicted via peptides deregulated in severe COVID-19 only, like CD99 antigen.This study's findings align with prior research reporting an association of urinary peptides (or classifiers based on theme) with unfavourable outcome.A Pubmed search using the keywords (urine OR urinary) AND (peptidom* OR proteom*) AND (death OR mortality) in the title or abstract resulted in 96 publications.After a manual assessment by three authors, 11 manuscripts were found to be relevant.These studies explored the association of urinary peptides with mortality in humans and include those describing the development of COV50 [11,12].Currie et al. described a significant value of CKD273, a classifier based on 273 urinary peptides, in predicting mortality in 155 microalbuminuric type 2 diabetic patients [46].Similar results were presented by Verbeke et al., linking CKD273 to mortality in 451 chronic kidney disease patients [36].Nkuipou-Kenfack et al. reported an association between urinary peptides and death, developing a classifier to predict mortality after ICU stay in 1243 patients [38].In 2021, Martens et al. described the connection between multiple urinary peptides, including many collagen fragments, biological age, and mortality [18].Batra et al. presented a proteomics-based mortality signature in COVID-19 and acute respiratory distress syndrome patients [47].In the context of hepatocellular carcinoma, Bannaga et al. identified several urine peptides being significantly associated with death [48].Recently, Wei et al. reported on the detection of urinary peptides related to pulsewave velocity also linked to mortality [49].In a robust study involving 1170 patients that underwent cardiac surgery, Piedrafita and colleagues identified 204 urinary peptides associated with acute kidney injury [50].A classifier based on these 204 peptides was validated in an independent cohort of 1569 ICU patients, demonstrating good performance and significant association with mortality.
Collagen peptides were consistently prominent biomarkers across many of these studies, with reduced abundance being associated with an increased risk of death, as also demonstrated by He and colleagues in the context of heart failure [25].Data from large cohorts in ICU and subjects not in critical condition at the time of sampling indicate that urinary peptides and classifiers derived thereof hold significant predictive value for a patient-relevant endpoint: death.In line with previous studies, the prediction of death appears to predominantly rely on collagen fragments, potentially reflecting attenuation of collagen degradation, and consequently progressing fibrotic processes.Evidently, the COV50 classifier was not designed to predict death in the general population.Additionally, considering the observation in this study that several peptides contained in this classifier show opposite regulation on predicting critical COVID-19 or death from any cause, it is to be expected that a classifier designed exclusively for death prediction, based solely on peptides significantly associated with death, could be of substantial value in guiding death-preventing interventions.Such a classifier is likely to be based mainly on collagen fragments.
The study has limitations.It relies on previously generated datasets; however, the large number of datasets, the high number of endpoints assessed, and the very high significance level of the findings strongly support the generalizability of the results.In fact, a strength of the study is the inclusion of datasets from various studies, underscoring the robust basis for this assessment.

Conclusions
Collectively, this study demonstrates that the urinary COV50 classifier is significantly associated with future death in both ICU and non-ICU patients, allowing for the identification of "vulnerable" subjects, irrespective of the underlying conditions.Further research is necessary to assess whether specific, personalized intervention guided by urinary collagen fragments can significantly improve outcomes, ultimately reducing the risk of future mortality.
• fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year

•
At BMC, research is always in progress.

Learn more biomedcentral.com/submissions
Ready to submit your research Ready to submit your research ?Choose BMC and benefit from: ? Choose BMC and benefit from:

Fig. 1
Fig. 1 The distribution of COV50 scores.A Density of the COV50 distribution in ICU and non-ICU subjects.B Mortality per person-years for FROG and non-ICU cohorts given age and COV50.C Mortality as share [0-1] from a logistic regression for FROG and non-ICU cohorts given age and COV50 score

Fig. 2
Fig. 2 Random effects meta-analyses based on the log-HR and the standard errors from the separate cox regressions.A Unadjusted, B adjusted for sex, age, kidney function and BMI.The size of dot symbols is proportional to weight and weight is inverse proportional to HR standard error

Table 1
Descriptive statistics for the ICU and non-ICU samples analysed within this study Categorical variables are described with absolute (N) and group-wise relative frequencies (%), continuous variables with median (IQR).P-values for group differences result from chi-squared homogeneity tests for categorical and for Wilcoxon rank sum test for continuous variables BMI body mass index, BP blood pressure, eGFR estimated glomerular filtration rate, FU follow-up, ICU intensive care unit, yrs years

Table 2
Estimates from the pooled adjusted Cox regression CI Confidence Interval, CKD chronic kidney disease, eGFR estimated glomerular filtration rate, HR Hazard Ratio, ICU intensive care unit; MAP mean arterial pressure EffectGroup (Non-ICU Subgroup) HR 95% CI p Age 1.83 1.43, 2.35 < 0.001

Table 3
Urinary peptides included in the COV50 classifier