Abstract

Incorporating biomarkers into cardiovascular studies, including nursing research, is a common approach when identifying underlying mechanisms and providing targets for intervention. However, effective utilization of biomarker data demands careful consideration. In the analysis, interpretation, and reporting phase, there are many facets to consider, including non-normality of the data, normalization procedures, and potential confounding influences of other clinical data. Furthermore, as many studies focus on patient-reported outcomes (PROs), it is important that the analysis and interpretation of biomarkers in relation to PROs is rigorous and reproducible. In this article, Part 2 of 2, we provide an overview of considerations for the analysis, interpretation, and reporting phases of biomarker studies. We also provide an example of these steps.

Learning objectives
  • Develop the steps for initial biomarker data analysis and reporting.

  • Identify best approaches for analysing biomarker data, especially in relation to patient-reported outcomes.

  • Present the biomarker findings in a rigorous, reproducible manner.

Introduction

Incorporating biomarkers into cardiovascular nursing research is a common approach when identifying underlying mechanisms and providing targets for intervention. Moreover, the inclusion of biomarkers into nursing research studies may help advance the movement towards precision health related to patient-reported outcomes (PROs).1 The inclusion of biomarkers can be as simple as a single plasma protein [e.g. troponin or N-terminal pro-B-type natriuretic peptide (NT-proBNP)] that is measured in a clinical core lab, or as complex as running multiple assays for various genetic, transcriptomic, metabolomic, or proteomic markers. As described in Part 1, many considerations are needed to ensure a rigorous and thorough methodology for identifying a biomarker (or multiple biomarkers), designing experiments, selecting an assay, and processing samples. After the biomarker data are resulted, there are additional considerations in analysing, interpreting, and reporting the data, particularly in the context of other clinical data. Furthermore, as many cardiovascular nursing research studies focus on PROs, it is important the analysis and interpretation of biomarkers in relation to PROs are rigorous and reproducible. In this second part, we describe how to analyse, interpret, and report biomarker data that are resulted from sample processing (Central Illustration; Box 1 for common definitions; Table 1 for examples). We also present an example to illustrate the steps and considerations.

Overview of the eight steps for incorporating biomarkers in cardiovascular nursing research.
Central Illustration

Overview of the eight steps for incorporating biomarkers in cardiovascular nursing research.

Box 1
Box of definitions
Intra-assay coefficient of variationMeasurement of the variance between data points within an assay. It determines the precision within a single run.
Inter-assay coefficient of variationMeasurement of the variance between different runs or experiments. It determines the precision across multiple runs or batches.
Natural log transformationThe natural log (ln) transformation can be used to transform data to make skewed data more normally distributed.
Intra-assay coefficient of variationMeasurement of the variance between data points within an assay. It determines the precision within a single run.
Inter-assay coefficient of variationMeasurement of the variance between different runs or experiments. It determines the precision across multiple runs or batches.
Natural log transformationThe natural log (ln) transformation can be used to transform data to make skewed data more normally distributed.
Intra-assay coefficient of variationMeasurement of the variance between data points within an assay. It determines the precision within a single run.
Inter-assay coefficient of variationMeasurement of the variance between different runs or experiments. It determines the precision across multiple runs or batches.
Natural log transformationThe natural log (ln) transformation can be used to transform data to make skewed data more normally distributed.
Intra-assay coefficient of variationMeasurement of the variance between data points within an assay. It determines the precision within a single run.
Inter-assay coefficient of variationMeasurement of the variance between different runs or experiments. It determines the precision across multiple runs or batches.
Natural log transformationThe natural log (ln) transformation can be used to transform data to make skewed data more normally distributed.
Table 1

Examples of applying Steps 5–8a

Step 5: perform the initial analysis and reportStep 6: analyse biomarker dataStep 7: consider covariatesStep 8: interpret and present results
Example 1 Butts et al. (2023)2Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersANCOVA was used to examine between-group differences at 3 months. Paired t tests were used to examine within-group differences over time. Pearson correlations were used to examine the relationship between telomere length and cytokines. Effect sizes were calculated using Hedges’ g for paired t tests and η2 for analysis of covariance modelsAdjusted for baseline values with ANCOVA. As this was a pilot study with a small sample (n = 32), additional covariates were not includedTotal telomere length increased and plasma IL-1β levels decreased in the exercise group from baseline to 3 months. Total telomere length was negatively associated with IL-1β at baseline
Example 2 Butts et al. (2023)3Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersZero-truncated Poisson regression models were constructed to estimate the effect of a predictor variableSociety of Thoracic Surgeons-Predicted Risk of Morbidity pre-operative risk score was included in multi-variate modelsMultivariable Poisson regression models with 4-h PCF chymase activity plus risk score as predictor variables demonstrate best model fitting for prediction of ICU and total hospital length of stay
Example 3 Denfeld et al. (2021)4sST2 was interpolated from a standard curve. Descriptive statistics for all variables, and biomarker data were natural log-transformedLatent growth curve modelling was used to estimate change in NT-proBNP and sST2 from pre- to 1, 3, and 6 months post-LVAD implantation. Quantified the effect of gender on intercepts and each shape or phase of changeThere were very few significant differences between women and men pre-LVAD, and thus, did not adjust for covariatesNT-proBNP: women and men had similar values of NT-proBNP at baseline and through 1 month post-implantation, but NT-proBNP decreased at a steeper rate from 1 to 6 months post-implantation compared with women.
sST2: women and men had similar values of sST2 at baseline, but the trajectories of change thereafter were significantly different as sST2 increased and then decreased for women, whereas for men, it just decreased
Example 4 Butts et al. (2019)5Descriptive statistics were calculated for all study variables and data were reviewed for normality assumptions and outliers in preparation for analysisMultiple linear regression was used to examine linear relationships. Student’s t-tests were used to compare groups (normotensive vs. resistant hypertension)Adjusted for age, sex, and body mass indexXanthine oxidase activity was increased two-fold in resistant hypertension vs. normal and was positively associated with left ventricular mass, left ventricular diastolic function, and 24-h urinary sodium
Worked examplePerformed descriptive analyses; transformed data to general normal distributions; calculated HOMA-IRUsing generalized linear modelling, examined associations between biomarkers and HRQOLAdjusted for severity of heart failure (Seattle Heart Failure Model Score) and comorbidities (Charlson Comorbidity Index)Worse HRQOL may be associated with insulin resistance or related pathways
Step 5: perform the initial analysis and reportStep 6: analyse biomarker dataStep 7: consider covariatesStep 8: interpret and present results
Example 1 Butts et al. (2023)2Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersANCOVA was used to examine between-group differences at 3 months. Paired t tests were used to examine within-group differences over time. Pearson correlations were used to examine the relationship between telomere length and cytokines. Effect sizes were calculated using Hedges’ g for paired t tests and η2 for analysis of covariance modelsAdjusted for baseline values with ANCOVA. As this was a pilot study with a small sample (n = 32), additional covariates were not includedTotal telomere length increased and plasma IL-1β levels decreased in the exercise group from baseline to 3 months. Total telomere length was negatively associated with IL-1β at baseline
Example 2 Butts et al. (2023)3Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersZero-truncated Poisson regression models were constructed to estimate the effect of a predictor variableSociety of Thoracic Surgeons-Predicted Risk of Morbidity pre-operative risk score was included in multi-variate modelsMultivariable Poisson regression models with 4-h PCF chymase activity plus risk score as predictor variables demonstrate best model fitting for prediction of ICU and total hospital length of stay
Example 3 Denfeld et al. (2021)4sST2 was interpolated from a standard curve. Descriptive statistics for all variables, and biomarker data were natural log-transformedLatent growth curve modelling was used to estimate change in NT-proBNP and sST2 from pre- to 1, 3, and 6 months post-LVAD implantation. Quantified the effect of gender on intercepts and each shape or phase of changeThere were very few significant differences between women and men pre-LVAD, and thus, did not adjust for covariatesNT-proBNP: women and men had similar values of NT-proBNP at baseline and through 1 month post-implantation, but NT-proBNP decreased at a steeper rate from 1 to 6 months post-implantation compared with women.
sST2: women and men had similar values of sST2 at baseline, but the trajectories of change thereafter were significantly different as sST2 increased and then decreased for women, whereas for men, it just decreased
Example 4 Butts et al. (2019)5Descriptive statistics were calculated for all study variables and data were reviewed for normality assumptions and outliers in preparation for analysisMultiple linear regression was used to examine linear relationships. Student’s t-tests were used to compare groups (normotensive vs. resistant hypertension)Adjusted for age, sex, and body mass indexXanthine oxidase activity was increased two-fold in resistant hypertension vs. normal and was positively associated with left ventricular mass, left ventricular diastolic function, and 24-h urinary sodium
Worked examplePerformed descriptive analyses; transformed data to general normal distributions; calculated HOMA-IRUsing generalized linear modelling, examined associations between biomarkers and HRQOLAdjusted for severity of heart failure (Seattle Heart Failure Model Score) and comorbidities (Charlson Comorbidity Index)Worse HRQOL may be associated with insulin resistance or related pathways

ANCOVA, analysis of covariance; HF, heart failure; HOMA-IR, homeostatic model assessment for insulin resistance (HOMA-IR); HRQOL, health-related quality of life; ICU, intensive care unit; IL-1β, interleukin 1 beta; LVAD, left ventricular assist device; NT-proBNP, N-terminal B-type natriuretic peptide; sST2, soluble suppressor of tumorgenicity.

aSteps 1–4 are presented in Part 1.

Table 1

Examples of applying Steps 5–8a

Step 5: perform the initial analysis and reportStep 6: analyse biomarker dataStep 7: consider covariatesStep 8: interpret and present results
Example 1 Butts et al. (2023)2Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersANCOVA was used to examine between-group differences at 3 months. Paired t tests were used to examine within-group differences over time. Pearson correlations were used to examine the relationship between telomere length and cytokines. Effect sizes were calculated using Hedges’ g for paired t tests and η2 for analysis of covariance modelsAdjusted for baseline values with ANCOVA. As this was a pilot study with a small sample (n = 32), additional covariates were not includedTotal telomere length increased and plasma IL-1β levels decreased in the exercise group from baseline to 3 months. Total telomere length was negatively associated with IL-1β at baseline
Example 2 Butts et al. (2023)3Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersZero-truncated Poisson regression models were constructed to estimate the effect of a predictor variableSociety of Thoracic Surgeons-Predicted Risk of Morbidity pre-operative risk score was included in multi-variate modelsMultivariable Poisson regression models with 4-h PCF chymase activity plus risk score as predictor variables demonstrate best model fitting for prediction of ICU and total hospital length of stay
Example 3 Denfeld et al. (2021)4sST2 was interpolated from a standard curve. Descriptive statistics for all variables, and biomarker data were natural log-transformedLatent growth curve modelling was used to estimate change in NT-proBNP and sST2 from pre- to 1, 3, and 6 months post-LVAD implantation. Quantified the effect of gender on intercepts and each shape or phase of changeThere were very few significant differences between women and men pre-LVAD, and thus, did not adjust for covariatesNT-proBNP: women and men had similar values of NT-proBNP at baseline and through 1 month post-implantation, but NT-proBNP decreased at a steeper rate from 1 to 6 months post-implantation compared with women.
sST2: women and men had similar values of sST2 at baseline, but the trajectories of change thereafter were significantly different as sST2 increased and then decreased for women, whereas for men, it just decreased
Example 4 Butts et al. (2019)5Descriptive statistics were calculated for all study variables and data were reviewed for normality assumptions and outliers in preparation for analysisMultiple linear regression was used to examine linear relationships. Student’s t-tests were used to compare groups (normotensive vs. resistant hypertension)Adjusted for age, sex, and body mass indexXanthine oxidase activity was increased two-fold in resistant hypertension vs. normal and was positively associated with left ventricular mass, left ventricular diastolic function, and 24-h urinary sodium
Worked examplePerformed descriptive analyses; transformed data to general normal distributions; calculated HOMA-IRUsing generalized linear modelling, examined associations between biomarkers and HRQOLAdjusted for severity of heart failure (Seattle Heart Failure Model Score) and comorbidities (Charlson Comorbidity Index)Worse HRQOL may be associated with insulin resistance or related pathways
Step 5: perform the initial analysis and reportStep 6: analyse biomarker dataStep 7: consider covariatesStep 8: interpret and present results
Example 1 Butts et al. (2023)2Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersANCOVA was used to examine between-group differences at 3 months. Paired t tests were used to examine within-group differences over time. Pearson correlations were used to examine the relationship between telomere length and cytokines. Effect sizes were calculated using Hedges’ g for paired t tests and η2 for analysis of covariance modelsAdjusted for baseline values with ANCOVA. As this was a pilot study with a small sample (n = 32), additional covariates were not includedTotal telomere length increased and plasma IL-1β levels decreased in the exercise group from baseline to 3 months. Total telomere length was negatively associated with IL-1β at baseline
Example 2 Butts et al. (2023)3Descriptive statistics were analysed for the study variables, and the data were reviewed for normality assumptions and outliersZero-truncated Poisson regression models were constructed to estimate the effect of a predictor variableSociety of Thoracic Surgeons-Predicted Risk of Morbidity pre-operative risk score was included in multi-variate modelsMultivariable Poisson regression models with 4-h PCF chymase activity plus risk score as predictor variables demonstrate best model fitting for prediction of ICU and total hospital length of stay
Example 3 Denfeld et al. (2021)4sST2 was interpolated from a standard curve. Descriptive statistics for all variables, and biomarker data were natural log-transformedLatent growth curve modelling was used to estimate change in NT-proBNP and sST2 from pre- to 1, 3, and 6 months post-LVAD implantation. Quantified the effect of gender on intercepts and each shape or phase of changeThere were very few significant differences between women and men pre-LVAD, and thus, did not adjust for covariatesNT-proBNP: women and men had similar values of NT-proBNP at baseline and through 1 month post-implantation, but NT-proBNP decreased at a steeper rate from 1 to 6 months post-implantation compared with women.
sST2: women and men had similar values of sST2 at baseline, but the trajectories of change thereafter were significantly different as sST2 increased and then decreased for women, whereas for men, it just decreased
Example 4 Butts et al. (2019)5Descriptive statistics were calculated for all study variables and data were reviewed for normality assumptions and outliers in preparation for analysisMultiple linear regression was used to examine linear relationships. Student’s t-tests were used to compare groups (normotensive vs. resistant hypertension)Adjusted for age, sex, and body mass indexXanthine oxidase activity was increased two-fold in resistant hypertension vs. normal and was positively associated with left ventricular mass, left ventricular diastolic function, and 24-h urinary sodium
Worked examplePerformed descriptive analyses; transformed data to general normal distributions; calculated HOMA-IRUsing generalized linear modelling, examined associations between biomarkers and HRQOLAdjusted for severity of heart failure (Seattle Heart Failure Model Score) and comorbidities (Charlson Comorbidity Index)Worse HRQOL may be associated with insulin resistance or related pathways

ANCOVA, analysis of covariance; HF, heart failure; HOMA-IR, homeostatic model assessment for insulin resistance (HOMA-IR); HRQOL, health-related quality of life; ICU, intensive care unit; IL-1β, interleukin 1 beta; LVAD, left ventricular assist device; NT-proBNP, N-terminal B-type natriuretic peptide; sST2, soluble suppressor of tumorgenicity.

aSteps 1–4 are presented in Part 1.

Step 5: perform the initial analysis and report

How a biomarker is initially analysed and reported will depend on a number of factors, including the biological material from which it was derived (e.g. plasma or tissue), the type of assay that was used, and any modifications needed after raw values are obtained (e.g. normalization, log transformation). First, there are several pre-analytic steps to produce the absolute value of the biomarker, depending on the experimental procedure. Oftentimes this involves using a standard curve to quantify the biomarker. For example, with an enzyme-linked immunosorbent assay (ELISA), a standard curve is generated based on standards of known concentrations tested alongside the samples. The standard curve then allows you to quantify the concentration of each sample based on readouts from the plate reader. At this point, you have ‘raw values’ of biomarkers. Biomarkers are typically presented as a certain amount per volume (fluid) or mass (tissue), also called a concentration (e.g. pg/mL). It is critical that the units are always reported alongside the biomarker values.

After obtaining the raw values of the biomarker, the next step is to perform simple descriptive statistics to understand range and variability, including any outliers. In order to minimize any potential confounding variability introduced across multiple assays, sometimes it is necessary to normalize the biomarker value across all assays. Normalizing across assays is possible when quality controls are included across assays (e.g. including the same pooled plasma across assays). One simple approach to normalization is to divide the biomarker value by the average of the quality control value in that assay. Of note, some assays, including the newer multiplex assays and mRNA sequencing, require normalization procedures to produce relative (and not absolute) values of biomarkers (e.g. normalized protein expression for multiplex assays6). Also, when quality controls are used across assays, intra- and inter-assay coefficients of variability can be calculated (see Box 1 for definition and calculations). These are typically reported in the methods section or supplementary material.

Oftentimes biomarker data are severely skewed, necessitating additional data manipulation, such as log transformation, to achieve normally distributed data suitable for common regression techniques. Skewness can be evaluated using visualization with histograms and skewness tests available in many common statistical software programs. Additionally, it may be helpful to create ratios of two related biomarkers (e.g. blood urea nitrogen-to-creatinine ratio) or calculate additional equations [e.g. homeostatic model assessment for insulin resistance (HOMA-IR) calculation7] to better understand the biological process. Throughout all of these steps, it is imperative to have rigorous and transparent methodology to ensure results are reliable and reproducible. Key steps include retaining all raw data files, tracking analytic steps and decisions, and having two investigators separately analyse the data to compare and validate findings.

Step 6: analyse biomarker data in relation to patient-reported outcomes

While biomarker data in and of themselves are interesting data points, and may in fact be the final step of analysis, understanding the biomarker data in the context of clinical information is very helpful in translating the findings. Furthermore, many readers may be able to understand common clinical biomarkers (e.g. troponin), but it may be difficult when novel biomarkers are measured; thus, it is important to answer the question ‘What do the data mean?’ Moreover, analysing biomarker data, common or novel, in the context of clinical information helps advance science towards precision health.

First, biomarker data (both raw and transformed data) are merged with other available data such as sociodemographic and clinical characteristics, PRO data, and so on. The type of analysis will depend on the research question. Common cross-sectional techniques involve generalized linear modelling techniques (e.g. linear or logistic regression, mediation analysis, moderation analysis). When many biomarkers are measured, especially those from the same or similar pathways (e.g. inflammation), multiple regression models may need to be evaluated to generate parsimonious models. For example, Kazukauskiene et al.8 examined associations between biomarkers of inflammation, thyroid function, and NT-proBNP (all natural log transformed) and health-related quality of life (HRQOL) among patients with coronary heart disease and heart failure. They first used univariate regression models to identify significant associations and then multivariate models including those significant in bivariate testing and adjusting for clinical characteristics.

Notably, there are multiple approaches within the family of generalized linear models that extend linear and logistic regression to data that have distributions other than normal.9 If longitudinal data are available, multiple techniques can be used such as mixed effects and latent growth curve modelling. There are many benefits to using latent growth curve modelling to analyse biomarker data, such as the ability to accommodate non-normally distributed data, missing data, and unevenly spaced time points (see example in Denfeld et al.4). Biomarker data can also be incorporated into mixture modelling approaches, such as latent class10 and growth mixture modelling,11 as well as survival analysis12 (see example in Willinger et al.13). Another consideration is the benefits and drawbacks to analysing biomarker data within cross-sectional vs. longitudinal designs. While cross-sectional biomarker studies can provide insight into associations between biomarkers (especially novel biomarkers) and outcomes of interest, they are limited in their ability to provide insights into the natural progression of a biomarker, its variability, and its relationship to disease development or progression. However, longitudinal data are needed in biomarker research due to their ability to capture patterns, trends, and individual variations over time, offering valuable insights into disease progression, treatment effectiveness, and personalized healthcare.

Step 7: consider other important covariates

When working with human models of disease, especially involving PROs, it is necessary to consider and incorporate other factors that may influence or confound the relationship being examined. When examining the association between a biomarker and PROs, it is likely necessary to adjust for other relevant clinical characteristics that may influence the relationship; for example, using a composite score of relevant clinical variables when examining associations between biomarkers and symptoms14 or frailty.15 Furthermore, many biomarker research questions focus on comparing two known groups, such as males vs. females.4,16 While simply comparing groups is one approach, a common criticism is not considering confounding variables that would explain the group differences. Similar to above, one approach is to include these confounding variables (single or a composite) in a regression model; however, with small sample sizes (a common occurrence in biomarker studies), you are limited in the number of variables that can be entered into a regression model. One way to approach this is through matching participants on key characteristics (e.g. age, sex, severity of disease) and then comparing biomarker values. The matching procedures can be as simple as matching on age or as complex as generating propensity scores17 for each individual and then matching on scores.

One important consideration when analysing biomarker data is related to sex as a biological variable (SABV). It is important to understand SABV from multiple angles, such as understanding if the biomarker is fundamentally different between males and females either physiologically or pathophysiologically (i.e. sexual dimorphisms) or if a certain value of the biomarker indicates different processes or manifestations between sexes. While there are a number of ways to approach SABV, the ‘must do’s’ include (i) including both sexes in the study unless there is a strong rationale not to include one sex (rare occasions), (ii) collecting relevant data related to SABV such as hormonal therapies or cycle phase, and (iii) communicating sex differences and similarities.18

Step 8: interpret and present the final results

The amount of time it takes to interpret the final results should not be underestimated, and this is where it is important to include collaborators. A first question to ask is: ‘what does the association mean?’ Sometimes the associations are intuitive at first glance; for example, higher beta-adrenergic receptor kinase-1, an enzyme that downregulates and desensitizes beta-adrenergic receptors, is associated with worse physical heart failure symptoms.14 But sometimes the associations are counterintuitive; for example, lower myostatin, a skeletal muscle cytokine, is associated with higher levels of physical frailty.15 With either scenario, it is critical to be well-versed in the literature to tease apart potential explanations for your findings. One key question to ask when reviewing the literature in your topic area is how did they measure, analyse, and report the biomarkers, and how is this similar to or different from your biomarker findings?

The presentation of the final results depends on the research question and can be tailored to the target audience using tables and figures. For biomarker data, you can consider presenting both raw and transformed data in order to have both clinical applicability and analytic transparency. Graphs of the data are very helpful in displaying the findings overall or by group, either at one time point or over time. For effective scientific data presentation using figures, prioritize clarity and simplicity: (i) use clear labels and titles, employing legible fonts and appropriate sizing; (ii) ensure axes are labelled with units and have clear scales; (iii) match the right type of graph (e.g. bar for group comparisons, scatter plot for linear associations) with the data; (iv) minimize clutter by avoiding unnecessary elements and emphasizing key findings using colours or annotations; (v) maintain consistency in style and formatting across figures for coherence; and (iv) provide detailed legends and captions to aid understanding. Remember, the goal is to convey complex information succinctly, fostering comprehension, and engagement for your audience.

In tandem with the presentation of the findings, the methods will need to be described in sufficient detail, including biomarker selection rationale, sample collection, processing (including type of assay used), and analytic steps. As word count limitations may influence the description, we suggest citing method papers and utilizing supplemental materials to ensure clarity and transparency.

Example application

To illustrate the analysis and interpretation of biomarker data in relation to PROs, we utilized a dataset from a previous study15 and ask the research question: ‘Is plasma glucose, insulin, or insulin resistance (i.e. HOMA-IR) associated with HRQOL among patients with heart failure?’. We assayed both fasting glucose and insulin concentrations in the plasma of patients with heart failure. As previously described, glucose was measured using a colorimetric assay (BioAssay Systems; Hayward, CA, USA), and insulin was measured using an ELISA (Mercodia; Uppsala, Sweden). To generate HOMA-IR, we used the equation: [fasting glucose (mmol/L) × fasting insulin (μU/mL)/22.5].7

The sample has been previously described.15 In the initial examination of biomarker data, insulin concentrations ranged from 1.5 to 66.0 mU/L [median 9.4 (interquartile range 5.0–16.7)] with significant skewness (Figure 1A). Glucose concentrations ranged from 71.3 to 319 mg/dL [median 105.1 (interquartile range 95–126)] also with significant skewness (Figure 1B). Thus, insulin and glucose were natural log transformed to generate more normal distributions (Figure 1C and D). HOMA-IR was calculated (after converting glucose from mg/dL to mmol/L), and this ranged from 0.32 to 30.60 [median 2.39 (interquartile range 1.30–5.48)]. Generalized linear modelling was used to analyse the association between the biomarkers and Minnesota Living with Heart Failure Questionnaire (MLHFQ) scores: total, physical, and emotional scale.

Biomarker distribution before (panels A and B) and after (panels C and D) logarithmic transformation with normal curves superimposed.
Figure 1

Biomarker distribution before (panels A and B) and after (panels C and D) logarithmic transformation with normal curves superimposed.

After adjusting for Seattle Heart Failure Model Score and Charlson Comorbidity Index (Table 2), the natural log of insulin was significantly associated with both total and emotional MLHFQ scores; that is, higher plasma insulin was associated with worse total and emotional HRQOL. Similarly, HOMA-IR was significantly associated with emotional MLHFQ scores; that is higher insulin resistance was associated with worse emotional HRQOL. There were no significant associations between the natural log of glucose and MLHFQ scores. These findings may indicate that worse HRQOL in HF, especially emotional HRQOL, may be associated with insulin resistance or related pathways (Figure 2). However, these are preliminary findings, and further analyses would need to investigate and pinpoint specific mechanisms.

Example application of analytical model and findings. HOMA-IR, homeostatic model assessment for insulin resistance; MLHFQ, Minnesota Living with Heart Failure Questionnaire. Created with BioRender.com
Figure 2

Example application of analytical model and findings. HOMA-IR, homeostatic model assessment for insulin resistance; MLHFQ, Minnesota Living with Heart Failure Questionnaire. Created with BioRender.com

Table 2

Associations between insulin and glucose and health-related quality of life

Total HRQOLPhysical HRQOLEmotional HRQOL
β Coefficient ± SEaP valueβ Coefficient ± SEaP valueβ Coefficient ± SEaP value
ln insulin (mU/L)0.23 ± 0.100.0260.21 ± 0.110.0500.37 ± 0.130.008
ln glucose (mg/dL)0.49 ± 0.290.1550.35 ± 0.300.2430.66 ± 0.380.078
HOMA-IR0.03 ± 0.010.0720.02 ± 0.010.1600.05 ± 0.020.010
Total HRQOLPhysical HRQOLEmotional HRQOL
β Coefficient ± SEaP valueβ Coefficient ± SEaP valueβ Coefficient ± SEaP value
ln insulin (mU/L)0.23 ± 0.100.0260.21 ± 0.110.0500.37 ± 0.130.008
ln glucose (mg/dL)0.49 ± 0.290.1550.35 ± 0.300.2430.66 ± 0.380.078
HOMA-IR0.03 ± 0.010.0720.02 ± 0.010.1600.05 ± 0.020.010

HOMA-IR, homeostatic model assessment for insulin resistance; HRQOL, health-related quality-of-life; ln, natural log; SE, Standard Error.

aAdjusting for Seattle Heart Failure Model projected 1-year survival and Charlson Comorbidity Index.

Table 2

Associations between insulin and glucose and health-related quality of life

Total HRQOLPhysical HRQOLEmotional HRQOL
β Coefficient ± SEaP valueβ Coefficient ± SEaP valueβ Coefficient ± SEaP value
ln insulin (mU/L)0.23 ± 0.100.0260.21 ± 0.110.0500.37 ± 0.130.008
ln glucose (mg/dL)0.49 ± 0.290.1550.35 ± 0.300.2430.66 ± 0.380.078
HOMA-IR0.03 ± 0.010.0720.02 ± 0.010.1600.05 ± 0.020.010
Total HRQOLPhysical HRQOLEmotional HRQOL
β Coefficient ± SEaP valueβ Coefficient ± SEaP valueβ Coefficient ± SEaP value
ln insulin (mU/L)0.23 ± 0.100.0260.21 ± 0.110.0500.37 ± 0.130.008
ln glucose (mg/dL)0.49 ± 0.290.1550.35 ± 0.300.2430.66 ± 0.380.078
HOMA-IR0.03 ± 0.010.0720.02 ± 0.010.1600.05 ± 0.020.010

HOMA-IR, homeostatic model assessment for insulin resistance; HRQOL, health-related quality-of-life; ln, natural log; SE, Standard Error.

aAdjusting for Seattle Heart Failure Model projected 1-year survival and Charlson Comorbidity Index.

When not to use a biomarker and alternative approaches

The decision to use a biomarker should be made carefully. For example, some biomarkers may be expensive to measure, and the cost-benefit ratio should be evaluated. Additionally, the decision to incorporate a biomarker into a study should be substantiated by ample evidence, with due regard to ethical implications. Moreover, when collecting samples in non-routine practice, the burden on participants must be considered. This involves being mindful of the potential impact on individuals contributing to research. Lastly, it is imperative to acknowledge that research, particularly involving biomarkers, results in significant waste. Therefore, a comprehensive examination of the consequences of biomarker usage on both the environment and the economy is essential.

Various alternative approaches exist for utilizing biomarkers in research without utilizing invasive procedures like venipuncture for sample collection. One option is to obtain stored samples from prior or ongoing studies, either from a collaborator or a biobank, for secondary analyses. Additionally, biospecimens can be collected from clinical waste products, such as discarded surgical tissue [e.g. cardiac tissue plug from left ventricular assist device (LVAD) surgery] or fluid from drains. Another avenue involves conducting secondary analyses using laboratory values from electronic medical records to explore the relationship between these biomarkers and patient outcomes. These methods offer valuable insights into biomarkers and their association with clinical or PROs. Given their often lower costs, these approaches are excellent for generating pilot data that can inform larger biomarker studies.

Conclusions

This Part 2 article emphasizes the important role of biomarkers in cardiovascular nursing research, offering a detailed roadmap for their incorporation and analysis. This paper, Part 2 of 2, focuses on crucial aspects of analysing biomarker data in relation to PROs, highlighting normalization and data manipulation. Furthermore, different statistical techniques and considerations for covariates were discussed to enable the translation of biomarker data into a clinical context. Lastly, it is essential to understand the literature and suggest potential explanations for your findings, while transparently presenting the results in figures or tables to facilitate comprehension. The two papers provide a comprehensive guide to equip researchers with a systematic approach, fostering precision health insights and enhancing the understanding of patient outcomes in cardiovascular nursing research.

Funding

This work is supported, in part, by Research Foundation Flanders (G072022N) to B.D., the National Institute of Nursing Research of the National Institutes of Health (R01NR019054) to Q.D., the Office of Research on Women’s Health and the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health (K12HD043488) to Q.D., and the National Institute on Aging of the National Institutes of Health (K23AG076977) to B.B. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Data availability

The data underlying this article will be shared on reasonable request to the corresponding author.

References

1

Hickey
 
KT
,
Bakken
 
S
,
Byrne
 
MW
,
Bailey
 
DCE
,
Demiris
 
G
,
Docherty
 
SL
, et al.  
Precision health: advancing symptom and self-management science
.
Nurs Outlook
 
2019
;
67
:
462
475
.

2

Butts
 
B
,
Hope
 
C
,
Herring
 
C
,
Mueller
 
K
,
Gary
 
RA
.
The effects of exercise on telomere length in persons with heart failure
.
J Cardiovasc Nurs
 
2023
. Published online ahead of print.

3

Butts
 
B
,
Goeddel
 
LA
,
Zheng
 
J
,
Pat
 
B
,
Powell
 
P
,
Mobley
 
J
, et al.  
Impact of early pericardial fluid chymase activation after cardiac surgery
.
Front Cardiovasc Med
 
2023
;
10
:
1132786
.

4

Denfeld
 
QE
,
Faulkner
 
KM
,
Davis
 
MR
,
Habecker
 
BA
,
Chien
 
CV
,
Gelow
 
JM
, et al.  
Exploring gender differences in trajectories of clinical markers and symptoms after left ventricular assist device implantation
.
Eur J Cardiovasc Nurs
 
2021
;
20
:
648
656
.

5

Butts
 
B
,
Calhoun
 
DA
,
Denney
 
TS
 Jr
,
Lloyd
 
SG
,
Gupta
 
H
,
Gaddam
 
KK
, et al.  
Plasma xanthine oxidase activity is related to increased sodium and left ventricular hypertrophy in resistant hypertension
.
Free Radic Biol Med
 
2019
;
134
:
343
349
.

6

Mitchell
 
A
,
Malmgren
 
L
,
Bartosch
 
P
,
McGuigan
 
FE
,
Akesson
 
KE
.
Pro-inflammatory proteins associated with frailty and its progression—a longitudinal study in community-dwelling women
.
J Bone Miner Res
 
2023
;
38
:
1076
1091
.

7

Matthews
 
DR
,
Hosker
 
JP
,
Rudenski
 
AS
,
Naylor
 
BA
,
Treacher
 
DF
,
Turner
 
RC
.
Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man
.
Diabetologia
 
1985
;
28
:
412
419
.

8

Kazukauskiene
 
N
,
Burkauskas
 
J
,
Macijauskiene
 
J
,
Mickuviene
 
N
,
Brozaitiene
 
J
.
Exploring potential biomarkers associated with health-related quality of life in patients with coronary artery disease and heart failure
.
Eur J Cardiovasc Nurs
 
2018
;
17
:
645
651
.

9

Lee
 
CS
,
Conway
 
C
.
The role of generalized linear models in handling cost and count data
.
Eur J Cardiovasc Nurs
 
2022
;
21
:
392
398
.

10

Lee
 
CS
,
Faulkner
 
KM
,
Thompson
 
JH
.
Identifying subgroups: part 1: patterns among cross-sectional data
.
Eur J Cardiovasc Nurs
 
2020
;
19
:
359
365
.

11

Lee
 
CS
,
Faulkner
 
KM
,
Thompson
 
JH
.
Identifying subgroups: part 2: trajectories of change over time
.
Eur J Cardiovasc Nurs
 
2020
;
19
:
444
450
.

12

Denfeld
 
QE
,
Burger
 
D
,
Lee
 
CS
.
Survival analysis 101: an easy start guide to analysing time-to-event data
.
Eur J Cardiovasc Nurs
 
2023
;
22
:
332
337
.

13

Willinger
 
L
,
Brudy
 
L
,
Häcker
 
A-L
,
Meyer
 
M
,
Hager
 
A
,
Oberhoffer-Fritz
 
R
, et al.  
High-sensitive troponin T and N-terminal pro-B-type natriuretic peptide independently predict survival and cardiac-related events in adults with congenital heart disease
.
Eur J Cardiovasc Nurs
 
2024
;
23
:
55
61
.

14

Denfeld
 
QE
,
Mudd
 
JO
,
Hasan
 
W
,
Gelow
 
JM
,
Hiatt
 
SO
,
Winters-Stone
 
K
, et al.  
Exploring the relationship between beta-adrenergic receptor kinase-1 and physical symptoms in heart failure
.
Heart Lung
 
2018
;
47
:
281
284
.

15

Denfeld
 
QE
,
Purnell
 
JQ
,
Lee
 
CS
,
Orwoll
 
ES
,
Camacho
 
SA
,
Hiatt
 
SO
, et al.  
Candidate biomarkers of physical frailty in heart failure: an exploratory cross-sectional study
.
Eur J Cardiovasc Nurs
 
2023
;
22
:
149
157
.

16

Lee
 
CS
,
Hiatt
 
SO
,
Denfeld
 
QE
,
Chien
 
CV
,
Mudd
 
JO
,
Gelow
 
JM
.
Gender-specific physical symptom biology in heart failure
.
J Cardiovasc Nurs
 
2015
;
30
:
517
521
.

17

Moons
 
P
.
Propensity weighting: how to minimise comparative bias in non-randomised studies?
 
Eur J Cardiovasc Nurs
 
2020
;
19
:
83
88
.

18

Usselman
 
CW
,
Lindsey
 
ML
,
Robinson
 
AT
,
Habecker
 
BA
,
Taylor
 
CE
,
Merryman
 
WD
, et al.  
Guidelines on the use of sex and gender in cardiovascular research
.
Am J Physiol Heart Circ Physiol
 
2024
;
326
:
H238
H255
.

Author notes

Conflict of interest: none declared.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/pages/standard-publication-reuse-rights)

Comments

0 Comments
Submit a comment
You have entered an invalid code
Thank you for submitting a comment on this article. Your comment will be reviewed and published at the journal's discretion. Please check for further notifications by email.