Plasma GAS6 predicts mortality risk in acute heart failure patients: insights from the DRAGON-HF trial

Abstract

Background

Growth arrest-specific 6 (GAS6) is a vitamin K-dependent protein related to inflammation, fibrosis, as well as platelet function. Genetic ablation of GAS6 in mice protects against cardiac hypertrophy and dysfunction. Nonetheless, the association between plasma GAS6 levels and acute heart failure (AHF) patients is still unknown.

Methods

We measured plasma GAS6 concentrations in 1039 patients with AHF who were enrolled in the DRAGON-HF trial (NCT03727828). Mean follow-up of the study was 889 days. The primary endpoint is all-cause death.

Results

In total, there were 195 primary endpoints of all-cause death and 135 secondary endpoints of cardiovascular death during the mean follow-up duration of 889 days. The higher levels of GAS6 were associated with higher rates of all-cause and cardiovascular death (P < 0.05). Baseline plasma GAS6 levels were still strongly correlated with clinical outcomes in different models after adjustment for clinical factors and N-terminal pro-brain natriuretic peptide (NT-proBNP, P < 0.05). GAS6 could further distinguish the risks of clinical outcomes based on NT-proBNP measurement.

Conclusion

Elevated plasma GAS6 levels were associated with an increased risk of all-cause and cardiovascular death in patients with AHF.

Trial registration NCT03727828 (DRAGON-HF trial) clinicaltrials.gov

Background

Heart failure (HF) is a devastating disease that remains a leading cause of death worldwide and is a primary socioeconomic burden for many countries. The mortality rate from acute heart failure (AHF) remains high despite current progress in HF treatment. As such, improved risk stratification and prognostic analysis tools are required to help guide the management of AHF patients. To achieve enhanced accuracy while not compromising convenience, novel biomarkers in plasma may enhance diagnostic and prognostic values in HF [1, 2].

In this regard, Growth arrest-specific 6 (GAS6) is a vitamin K-dependent protein which has been detected in many different tissues and cells [3,4,5]. Biological functions of GAS6 are regulated by interaction with the tyrosine kinase receptor family (TAM): TYRO3 Protein Tyrosine Kinase (Tyro3), AXL, MER Proto-Oncogene, Tyrosine Kinase (MERTK), with the highest affinity for AXL. The extracellular domain of these receptors can be proteolytically cleaved and released into the plasma in soluble forms to act as decoy receptors of GAS6, which reduces interaction between the ligand and receptor [6, 7]. GAS6/TAM interaction has been implicated in platelet function regulation [8], cell growth [9], phagocytosis of apoptotic cells [10], and reduction of inflammatory response [11]. GAS6 has been regarded as a risk factor for adverse effects in many diseases, such as systemic lupus erythematosus [12, 13], rheumatoid arthritis [14], liver cirrhosis [15] and venous thromboembolism [16, 17]. Although GAS6 showed relatively low expression levels in heart tissue under a normal situation, it has been demonstrated to increase significantly in hypertrophic and failing hearts [18]. Our previous study had demonstrated that Gas6 knockout would alleviate cardiac hypertrophy, fibrosis and contractile dysfunction [18]. A previous study also found GAS6 levels to be significantly increased at day 7 post-admission in ST elevation myocardial infarction (STEMI) patients [19]. Nonetheless, the association between GAS6 and clinical prognosis in patients with AHF remains unclear. In this study, we investigated the prognostic value of GAS6 in a cohort of AHF patients which enrolled in the DRAGON-HF trial.

Methods

Study population

The study enrolled 1039 patients between January 2017 and March 2020 with age > 18 yrs, symptomatic (NYHA class II-IV) heart failure. The diagnostic criteria of AHF referred to the 2016 ESC guidelines [20]. Patients with AHF was diagnosed by two clinicians based on a thorough history assessing symptoms, prior cardiovascular history and potential cardiac and non-cardiac precipitants. Assessment of signs /symptoms of congestion and /or hypoperfusion by physical examination and further confirmed by appropriate additional investigations such as ECG, chest X-ray, laboratory assessment (with specific biomarkers) and echocardiography. Patients with BNP < 100 pg/mL and/or NT-proBNP < 300 pg/mL were excluded for the diagnosis of AHF (based on sex, age, BMI and renal function). Patients with serious infectious diseases, severe hepatic and renal dysfunction, tumor and blood diseases were excluded from our study. In our study, we also identified AHF patients with renal dysfunction as evaluated by eGFR < 60 mL/min/1.73m², which constitutes a part of the observational, cohort trial: DRAGON-HF trial (Diagnostic, Risk Stratification and Prognostic Value of Novel Biomarkers in Patients with Heart Failure, NCT03727828). This study complied with the Declaration of Helsinki and was approved by the local ethics committee. Informed consent was obtained from all patients. Patients and the public were not involved in the design, or conduct, or reporting, or dissemination plans of our study.

Plasma sampling

Following provision of consent and with the patient recumbent for at least 15 min, 20 mL of venous blood was withdrawn by venipuncture of an antecubital vein and collected in evacuated tubes containing EDTA as anticoagulant. Blood samples were centrifuged at 1500g for 20 min at 4 ℃. Plasma was siphoned, aliquoted and shipped on dry ice to a central laboratory. Samples were stored at − 80 ℃ until analysis. At the time of assaying, plasma samples were defrosted at room temperature and analyzed in a single batch.

Biomarkers assessment

Circulating plasma concentration of GAS6 was measured by a commercially available ELISA assay (catalog: DY885B) as previously described [21]. According to the manufacturer, the limit of detection of this assay is 15.6 pg/mL and the upper reference limit is 1000 pg/mL. The intra-assay and inter-assay coefficients of variation were 6.5% and 8.5%, respectively. NT-proBNP levels were determined by the clinically available VITROS® NT-proBNP Assay (Ortho-Clinical Diagnostics, Raritan, NJ, USA).

Endpoints

The primary endpoint of the study is all-cause death. The secondary endpoint is cardiovascular death.

Statistical analysis

Continuous variables were presented as mean ± standard deviations or median (IQR). One way ANOVA test or Mann–Whitney U test was used to examine the difference between groups depending on the data distribution. Categorical variables were assessed by Chi-Squared test in different groups. The correlation between GAS6 and eGFR was visualized with scatterplots, and the Spearman correlation coefficient was calculated. Kaplan–Meier method was used to derive the cumulative endpoint-free survival estimates. Survival curves stratified according to the tertile distribution of GAS6 were compared with the log-rank test for trend. We also tested the association between GAS6 and long-term outcomes in different clinical models. Ln-transformed biomarker distributions were standardized to a mean of 0 and standard deviation of 1 to facilitate comparison of effect sizes between biomarkers. We first evaluated the GAS6 in models containing the following standard clinical risk factors: age, sex, regular cigarette smoking, body mass index, eGFR, presence of diabetes mellitus, hypertension, and atrial fibrillation, medical usage of diuretics and spironolactone. We then performed analyses incorporating all biomarkers together and add NT-proBNP, hs-TnT and hs-CRP to the clinical model. We also assessed the improvement of prediction value by IDI and NRI. The corresponding C-index, reflecting the discriminative ability of each biomarker up to mediate follow-up 889 days, was calculated. Proportional hazards assumption was confirmed by Schoenfeld’s test. All data were analyzed using IBM SPSS Statistics (version 20, 2011) and R statistical software (http://www.r-project.org/, version 4.1.0). The C-index, NRI, and IDI were calculated in R software, packages SurvC1 and Survival. Two-sided P-values < 0.05 were considered statistically significant.

Results

Patients characteristics

Baseline characteristics are shown in Table 1. Among 1039 patients recruited, 669 (64.4%) were men and median age was 71 yrs (IQR 64–81). Median follow-up time was 889 days. Mean plasma concentration of GAS6 in these patients was 10.34 ng/mL. Patients with higher GAS6 levels were older and more likely to exhibit lower left ventricular ejection fraction (LVEF), more severe HF (evaluated by NYHA functional class) and reduced renal function (assessed by eGFR). Furthermore, blood urea nitrogen and creatinine are directly proportionate to GAS6 levels. Meanwhile, a higher level of GAS6 in patients of AHF increases the prevalence of hypertension, alongside an incremental proportion of β-blocker, spironolactone and diuretic usage. Patients with high GAS6 levels also had enhanced hs-TnT, hs-C-reactive protein (hs-CRP) and decreased eGFR, although GAS6 levels were modestly correlated with hs-TnT (r = 0.196, P < 0.001), hs-CRP (r = 0.214, P < 0.001) and eGFR (r = − 0.137, p < 0.001). In contrast, GAS6 levels showed no correlation to body mass index, systolic blood pressure, heart rates, hemoglobin, smoking history, diabetes mellitus, atrial fibrillation, ACEI/ARBs, calcium channel blockers and statins.

Table 1 Baseline characteristics of patients

Baseline GAS6 levels and clinical outcomes

During a median follow-up of 889 days (interquartile range 728–1182 days), 195 (18.8%) patients died from all-causes and 128 died from cardiovascular causes. Endpoint-free survival stratified according to the distribution of GAS6 levels indicated that patients with lower GAS6 levels (Tertile 1; GAS6 < 8.60 ng/mL) had a significantly better clinical endpoint-free survival than patients with higher GAS6 levels (GAS6 ≥ 12.0 ng/mL) (as shown in Figs. 1 and 2). However, there’s no significant difference in HF re-hospitalization among tertiles (data not shown).

Fig. 1

Distribution of the number of patients and deaths according to GAS6 levels

Fig. 2

Kaplan–Meier curves for all-cause death and cardiovascular death according to GAS6 levels

Univariable and multivariable analysis with categorical and continuous GAS6 levels and adjustment models presented a significant relation between GAS6 and clinical endpoints (Table 2 and Additional file 1: Table S1). In model 1, GAS6 levels were associated with increasing all-cause death (HR, 1.65; 95% CI (1.41,1.93); P < 0.001) and cardiovascular death (HR, 1.71; 95% CI (1.40,2.09); P < 0.001). After additionally adjusting for clinical covariates, GAS6 levels remained associated with all-cause death (HR, 1.56; 95% CI (1.33,1.83); P < 0.001) and cardiovascular death (HR, 1.60; 95% CI (1.30,1.95); P < 0.001). After adding biomarkers into traditional clinical model, which includes age, sex, BMI, smoking, hypertension, diabetes, atrial fibrillation and eGFR, GAS6 levels were also significantly associated with all-cause death (HR, 1.44; 95% CI (1.22,1.70); P < 0.001) and cardiovascular death (HR, 1.44; 95% CI (1.17,1.79); P < 0.001).

Table 2 Associations between GAS6 and all-cause death or cardiovascular death, with adjustment for clinical characteristic

Biomarkers for prediction of all-cause death and cardiovascular death

Univariate and multivariate Cox proportional hazards analysis for clinical outcomes showed several variables were independent predictors of all-cause death and cardiovascular death in patients with AHF. For comparing the effect size of GAS6 and classic biomarkers in predicting the adverse outcome, we selected several variables as clinical risk factors and added biomarkers to the clinical model respectively. Separate incorporation of each biomarker into the clinical risk model showed that all biomarkers increased the C statistic for prediction of adverse events except CRP (C statistics for all-cause death: clinical risk alone 0.74, clinical risk + GAS6 0.77, clinical risk + NT-proBNP 0.78, clinical risk + hs-CRP 0.74 and clinical risk + hs-TNT 0.76. C statistics for cardiovascular death: clinical risk alone 0.74, clinical risk + GAS6 0.77, clinical risk + NT-proBNP 0.80, clinical risk + hs-CRP 0.75 and clinical risk + hs-TNT 0.77) (Table 3) Moreover, we examined the additive usefulness of GAS6 in multiple biomarkers strategy based on clinical risk factors and other biomarkers (C statistics for all-cause death: 0.80 and for cardiovascular death: 0.82). NRI and IDI showed the similar improvement of prediction ability by adding each or all biomarker to the clinical risk model.

Table 3 Associations between biomarkers and prognostic discrimination

Apart for the information provided by NT-proBNP, GAS6 provides beneficial prognostic information of AHF patients. Among the patients with baseline NT-proBNP less than the cohort median (< 1460.5 pg/mL), all-cause death was significantly different between the high (> 10.15 ng/mL) and low (< 10.15 ng/mL) GAS6 levels patients (as shown in Fig. 3; log-rank P < 0.001). With this subgroup, those with high GAS6 levels experienced all-cause death of 11.3% while those low GAS6 levels experienced mortality of 3.9%. Similarly, among patients with high NT-proBNP levels, mortality was 36.2% and 22.6% between GAS6 levels above and below the median, respectively. Therefore, circulating GAS6 levels may further distinguish the risks of adverse outcomes based on the NT-proBNP measurement. We found similar results when we use GAS6 and NT-proBNP to predict the prognostic information of cardiovascular death (Fig. 3; log-rank P < 0.001). Hs-TnT is another important biomarker of heart injury commonly used in clinical practice. We also found combination of GAS6 and hs-TnT could better differentiate patients at risk (Additional file 1: Figure S3).

Fig. 3

Kaplan–Meier curves for all-cause death and cardiovascular death according to GAS6 and NT-proBNP levels. Baseline median GAS6 was 10.15 ng/mL; the median level for NT-proBNP was 1460.5 pg/mL. ‘High’ refers to values above the median, while ‘low’ refers to values below the median

Analysis of renal dysfunction subgroups in patients with AHF

Previous study suggested that GAS6 was increased in CKD and kidney injury patients. We found GAS6 was negatively associated with eGFR (Additional file 1: Fig. S2). We classified AHF patients by history of renal dysfunction evaluated by eGFR < 60 mL/min/1.73m². Characteristics analysis indicated that patients with renal dysfunction who died were older, and more likely to had worse cardiac and renal function assessed by LVEF and eGFR respectively. Furthermore, these patients who died had increased GAS6 levels, as well as BUN, hs-TnT, hs-CRP and NT-proBNP, with decreased BMI and hemoglobin in comparison to live patients (Table 4). High GAS6 or/and NT-proBNP levels were also associated with higher mortality (Fig. 4), the association was similar between GAS6 and other biomarkers, including BUN, hs-TnT, and hs-CRP (Additional file 1: Fig. S4). Survival analysis also revealed that higher GAS6 levels were more likely to be associated with adverse endpoints (Fig. 5; log-rank P < 0.001).

Table 4 Characteristics of patients with AHF classified by eGFR levels

Fig. 4

Observed incidence of all-cause death and cardiovascular death. Biomarker stratifications are as follows: the cut-off value of GAS6 level was 10.15 ng/mL, the cut-off value of NT-proBNP level was 1460.5 pg/mL. ‘High’ refers to GAS6/NT-proBNP levels above the cut-off value, while ‘low’ refers to GAS6/NT-proBNP levels below the cut-off value

Fig. 5

Survival regarding the primary endpoint (all-cause death) and the secondary endpoint (cardiovascular death) stratified according to the tertile distribution of GAS6 in different renal function groups (based on eGFR)

Discussion

Our study demonstrated for the first time that elevated plasma GAS6 levels were not only associated with impaired LV function, but also with higher all-cause and cardiovascular death in patients with AHF, which was independent of the established cardiac biomarker NT-proBNP and hs-TnT. Furthermore, we also revealed that AHF patients with renal dysfunction had higher plasma GAS6 levels.

GAS6 is a vitamin K-dependent protein and a kind of ligand of TAM family of tyrosine kinase receptors. As a decoy receptor for GAS6, increased sAXL levels undermine the effects of the GAS6 /AXL ligand interaction [22]. We previously demonstrated that the expression of GAS6 was upregulated in human dilated cardiomyopathy, and Gas6 knockout in mice mitigated cardiac hypertrophy, fibrosis and contractile dysfunction. Conversely, cardiac-specific overexpression of Gas6 aggravated the adverse phenotype [18]. Recent data suggest that GAS6 may play a pivotal role as prognostic biomarker in cardiac remodeling and LV function. According to a previous study, GAS6 levels were significantly lower during admission in STEMI patients compared with controls, and increased only from day 7, but there was no prognostic analysis of GAS6 and other risk factors [19]. Our study, to the best of our knowledge, is the first to assess the role of plasma GAS6 levels in prediction of adverse outcomes in patients with AHF. Plasma GAS6 in the acute phase in AHF patients, is closely associated with several risk factors of adverse prognosis in AHF patients, including NYHA functional class, hypertension and renal dysfunction (evaluated by eGFR). Therefore, GAS6 may represent an easily accessible biomarker of unfavorable prognosis in patients with AHF, and may provide further prognostic information based on the existing indicators.

Heart failure is a heterogeneous disease and also a common end-stage syndrome of different cardiovascular disorders. Although pressure and volume overloading as well as myocardial infarction constitute the most common causes of cardiac remodeling and heart failure, the myocardium manifest as compensatory enlargement in the initial phase of injury [23]. In this regard, GAS6 has been found to be crucial in mediating the deleterious effects following cardiac injury albeit expressed at varying levels in different disorders. During volume overloading, the change in GAS6 may not be pronounced following a mild inflammation /fibrosis and eccentric hypertrophic response. However, genetic ablation of GAS6 alleviated cardiac hypertrophy and improved cardiac function induced by pressure overloading [18]. A six months follow-up of STEMI patients revealed a potential role of the GAS6-AXL system in the pathophysiology of left ventricular remodeling following STEMI [19]. With regards to AHF, the ensuing neurohormonal, cytokine, inflammatory and oxidative stress storm ultimately leads to end-organ damage. Biomarkers of end-organ damage and dysfunction, such as troponins, cystatin-C and transaminases, have been shown to be increased in AHF and associated with poor prognosis [24]. And GAS6 has been demonstrated to be related to these biomarkers. GAS6 is mainly generated by cardiac and extracardiac tissues, while being correlated to activation of inflammatory processes, hemodynamic compromise, plaque instability, loss of vessel integrity, and tissue injuries, which significantly determine the progression of HF. Previous studies have showed exposure of human endothelial cells to 10% cyclical stretch increased release of the GAS6 and GAS6 /Axl signaling in human monocytes potentiated interleukin 1-beta production [25]. GAS6 could also increase the production of reactive oxygen species (ROS) in vascular smooth muscle cells, and mediate increased migration in vascular dysfunction [26]. Though not an organ-specific protein, the association between GAS6 and myocardial injury indicates that GAS6 may potentially function as a new and additional indicator of injury severity in AHF patients.

By further classifying GAS6 and NT-proBNP with median levels, the prognostic risk-stratification of adverse outcomes may be further enhanced. As natriuretic peptide release is mainly stimulated by pressure and /or volume overload, this will imply that GAS6 may be involved in other pathophysiological pathways regulating myocardial adaptation and dysfunction, as opposed to solely relying on the BNP-related pathway. GAS6 may thus provide more clinical information in patients with AHF independent of pressure or volume overload. Although GAS6 levels were unlikely influenced by infectious inflammatory processes in our study, we also cannot preclude the possibility that GAS6 levels might be influenced by other unknown processes in our cohort patients.

AHF is often accompanied by renal insufficiency, which suggests a worse prognosis in AHF patients. However, the predictive value of plasma biomarkers in adverse cardiovascular events and renal dysfunction is poor. Biomarkers, such as blood urea nitrogen, creatinine and uric acid, cannot provide accurate information about renal injury. Patients with renal dysfunction usually deteriorate rapidly due to improper diagnosis and treatment. Due to limitations of our study, the etiology in patients with renal dysfunction in our cohort was not easily identified. Nevertheless, GAS6 was found to be negatively associated with eGFR, and significantly elevated in AHF patients with renal dysfunction, especially in dead patients. Survival analysis also revealed that higher levels of GAS6 were associated with lower survival rates. We also found that higher levels of GAS6 and clinical biomarkers, including NT-proBNP, TnT, BUN and CRP, were accompanied by highest mortality in AHF patients, especially in those with renal dysfunction. In addition, GAS6 plays a pivotal role in pulmonary, hepatic and renal inflammation and fibrosis [27,28,29]. In this regard, AHF patients with renal dysfunction whose circulating GAS6 levels may indicate a primary inflammatory/fibrotic susceptibility may benefit from optimization of anti-inflammatory/antifibrosis therapies, such as mineralcorticoid receptor antagonists.

Limitations

Our analysis has several limitations. First, there was no significant correlation between GAS6 stratification and HF re-hospitalization, which may lead to undue HF readmission. Second, there is limited information available on patient co-morbidities (especially chronic inflammatory conditions), which may potentially influence GAS6 levels. Third, serial measurement of GAS6 may provide more information about the association between GAS6 levels and prognosis of AHF patients. Fourth, we did not examine the effects of medical intervention on prognosis of GAS6. Anti-coagulation drugs, such as warfarin may influence GAS6 activation and levels. Fifth, we did not monitor continuous changes of eGFR nor determine the etiology of renal dysfunction.

Conclusion

GAS6 may prove to be a novel and additional clinical biomarker for risk stratification of AHF patients, especially in the cohort with renal dysfunction. A combination of GAS6 and other clinical biomarkers might further pave the path forward to formulate precision interventions targeting the different pathways. It may be worthwhile exploring the accuracy and efficiency of diagnosis and prognosis of GAS6 specifically in patients with other cardiovascular diseases in the future.