Why Link Biomarker Response to Clinical Efficacy?
A treatment may change a biomarker before a clinical endpoint can be measured, or the endpoint may be costly, slow, or difficult to observe. A biomarker can therefore help describe the treatment’s biological response and may contain information about later clinical outcomes. The development question is not simply whether exposure changes the biomarker. It is whether biomarker change is related to the clinical endpoint in a way that is useful, reliable, and appropriately interpreted.
A linked analysis usually considers at least two relationships: exposure with biomarker response, and biomarker response with clinical outcome. These relationships can support prediction and help organize evidence about a treatment’s effects. They do not, on their own, show that changing the biomarker causes clinical benefit or establish that the biomarker can replace the clinical endpoint.
Specify the Link Before Interpreting It
Let \(X\) denote a selected exposure metric, \(M\) a biomarker change, and \(Y\) a clinical outcome. For a continuous biomarker and continuous clinical endpoint, a simple pair of models is:
Here, \(\alpha_1\) describes the exposure–biomarker association. In the clinical model, \(\beta_1\) describes the expected difference in \(Y\) associated with a one-unit difference in \(M\), comparing individuals with the same modeled exposure \(X\). The coefficient \(\beta_2\) describes the exposure association with \(Y\) conditional on \(M\). These are model-based associations unless additional causal assumptions are justified.
The form of the clinical outcome model must suit the endpoint. A continuous endpoint may use a linear model; a binary endpoint may use a logistic model; and a time-to-event endpoint may use an appropriate survival model. The coefficient’s scale then depends on the model: for example, a logistic regression coefficient is on the log-odds scale, not the probability scale. The same interpretive caution applies across model types.
Define the analysis population, exposure metric, biomarker summary, clinical endpoint, and timing before fitting the models. The biomarker should be measured in a window that is meaningful for the clinical outcome; if effects are delayed, a contemporaneous biomarker value may not represent the relevant response history. Baseline values, background therapy, disease severity, and other common causes of biomarker and outcome should be considered based on the question and causal assumptions. Earlier material on causal assumptions in individual-level exposure–response is directly relevant: conditioning on a post-treatment biomarker can change the estimand and can introduce bias.
A useful prediction target is the expected clinical outcome for a specified exposure after using the exposure model to predict the biomarker. But uncertainty in both fitted models matters. Plugging a predicted biomarker into a clinical model gives a point prediction; it does not by itself provide a valid confidence interval for the complete prediction.
Worked Example: Translating Exposure Through a Biomarker Model
Suppose a hypothetical analysis uses exposure \(X\) in mg/L, biomarker change \(M\) in biomarker units, and clinical improvement \(Y\) in points, with higher \(Y\) indicating greater improvement. The fitted models are \(M = 0.2 + 1.2X\) and \(Y = 1.0 + 1.5M + 0.8X\). These coefficients are illustrative estimates, not results from a real study.
At \(X=0\) mg/L, the predicted biomarker change is \(0.2 + 1.2(0)=0.2\) biomarker units. The predicted clinical improvement is \(1.0 + 1.5(0.2) + 0.8(0)=1.3\) points.
At \(X=2\) mg/L, the predicted biomarker change is \(0.2 + 1.2(2)=2.6\) biomarker units. The predicted clinical improvement is \(1.0 + 1.5(2.6) + 0.8(2)=6.5\) points.
The model-implied difference between these exposure values is \(6.5-1.3=5.2\) points. Substituting the biomarker model into the clinical model gives a biomarker-linked component of \(1.5(1.2)(2)=3.6\) points and an exposure term of \(0.8(2)=1.6\) points; these sum to 5.2 points. This algebraic decomposition is descriptive of the specified models. Without stronger assumptions, the 3.6-point component is not a causal indirect effect mediated by the biomarker.
Quantifying Uncertainty in the Biomarker–Outcome Relationship
For a linear model with independent observations, a correctly specified mean relationship, and errors with mean zero and constant variance, the estimated slope is uncertain because the observed outcomes vary around the fitted line. A confidence interval expresses this sampling uncertainty. With normally distributed errors, the usual small-sample \(t\) interval also requires that the errors be normally distributed; sufficiently large samples may permit approximate inference under weaker distributional conditions.
For a simple linear regression of \(Y\) on \(M\), define \(S_{MM}=\sum_i(M_i-\bar{M})^2\), where \(\bar{M}\) is the sample mean of \(M\). The fitted slope is \(\hat{\beta}_1=S_{MY}/S_{MM}\), where \(S_{MY}=\sum_i(M_i-\bar{M})(Y_i-\bar{Y})\). If the residual sum of squares is \(SSE\), then the residual variance estimate is \(s^2=SSE/(n-2)\), and \(SE(\hat{\beta}_1)=s/\sqrt{S_{MM}}\). The \(100(1-\alpha)\%\) confidence interval is \(\hat{\beta}_1 \pm t_{1-\alpha/2,\,n-2}SE(\hat{\beta}_1)\).
The slope formula follows by choosing the intercept and slope that minimize the residual sum of squares. For any fixed slope, the minimizing intercept makes the residuals sum to zero, so the fitted line passes through \((\bar{M},\bar{Y})\). Write the fitted value as \(\bar{Y}+b(M_i-\bar{M})\). Differentiating the sum of squared residuals with respect to \(b\), setting the derivative to zero, and solving gives:
Worked Example: Confidence Interval for a Biomarker–Outcome Slope
In a hypothetical study of 16 participants, \(Y\) is clinical improvement in points and \(M\) is biomarker change in units. A simple regression of \(Y\) on \(M\) gives \(\hat{\beta}_1=1.20\) points per biomarker unit, \(S_{MM}=40\) biomarker units squared, and \(SSE=56\) points squared. Assume independent participants, a linear mean relationship, constant residual variance, and normally distributed residuals for the small-sample interval.
There are \(df=n-2=14\) degrees of freedom. The residual variance estimate is \(s^2=56/14=4\) points squared, so \(s=2\) points. The slope standard error is:
Using \(t_{0.975,14}=2.145\), the 95% confidence interval is:
Under this model, the data are compatible with a positive association: the interval ranges from 0.522 to 1.878 additional clinical improvement points per biomarker unit. It is not a 95% probability statement about the fixed slope. It also does not establish that increasing the biomarker would cause the clinical improvement.
Association Is Not Surrogate Validation
A biomarker can respond strongly to exposure and still be a poor indicator of clinical benefit. The exposure–biomarker relationship answers whether the biomarker changes with exposure. The biomarker–outcome relationship asks whether biomarker values and clinical outcomes co-vary. Neither question alone establishes that treatment effects on the biomarker reliably predict treatment effects on the clinical endpoint.
Individual-level association and trial-level surrogacy are different evidentiary questions. An individual-level association means that, within the analyzed population and model, people with different biomarker values tend to have different outcomes. Trial-level surrogacy concerns whether treatment effects on the biomarker across relevant trials predict treatment effects on the clinical endpoint. Even a strong individual-level association may arise because disease severity affects both measurements, or because treatment exposure affects both through separate pathways.
Worked Example: A Biomarker Response Without a Clinical Trend
Consider three hypothetical exposure groups with group-mean exposure values of 1, 2, and 3 mg/L. Their mean biomarker changes are 1, 3, and 5 units, while their mean clinical improvements are 4.0, 4.1, and 4.0 points.
The biomarker rises by 2 units for each 1 mg/L increase in mean exposure, so the group means show a clear exposure-related biomarker response. But the clinical improvement is nearly unchanged across the same groups. For these three group means, the mean clinical improvement is \(4.033\) points. The centered biomarker values are \(-2,0,2\), and the centered clinical values are approximately \(-0.033,0.067,-0.033\). Their cross-product sum is approximately \(0\), so a simple regression of the three clinical group means on the three biomarker group means has a slope of approximately \(0\) points per biomarker unit.
This small illustration is not an inferential analysis and cannot establish that the biomarker and clinical endpoint are unrelated in individuals. It shows why an exposure-linked biomarker change alone is insufficient to claim clinical benefit or surrogate validity. Proper assessment requires suitable individual-level and, when relevant, across-trial evidence, with uncertainty and study design taken into account.
Diagnostics, Alternatives, and Limitations
Before interpreting a linked model, check whether its assumptions and measurements are plausible. Useful checks include plots of clinical outcome against biomarker, residuals against fitted values and exposure, and observations with high influence. Look for curvature, unequal residual spread, outliers that dominate the slope, and differences in data availability across exposure or biomarker ranges. A narrow biomarker range can make the slope imprecise even when the fitted line appears steep. If measurement error in the biomarker is substantial, a naive regression can attenuate the estimated association; estimated exposure or biomarker uncertainty should not be ignored when it is material.
Timing and selection also matter. A biomarker measured after clinical deterioration may reflect the outcome rather than predict it. Restricting analysis to participants with both observed biomarker and clinical data can bias results if missingness is related to disease status or response. A model that omits a common cause of biomarker and outcome may report an association that is not attributable to the biomarker itself. These issues are reasons to state the target and assumptions, not reasons to relabel an association as causal.
The main alternatives answer different questions. Modeling the clinical endpoint directly as a function of exposure is often the clearest approach when the goal is exposure–efficacy characterization; it does not require the biomarker to explain the clinical response. A joint model can represent correlated biomarker and clinical outcomes and propagate uncertainty, but its conclusions still depend on its structure and assumptions. A mediation analysis targets a causal decomposition into pathways and requires additional causal assumptions beyond fitting a biomarker coefficient. Comparing treatment effects across trials is relevant to trial-level prediction but requires multiple informative trials and cannot be replaced by one individual-level regression.
A linked biomarker analysis is most useful as one part of an evidence chain: it can show whether the biomarker responds to exposure, whether it is associated with the clinical endpoint, and how uncertain those patterns are. It has limitations when the biomarker is measured imprecisely, the relationship changes across populations or time, the clinical endpoint is poorly represented by the model, or evidence is sparse at the exposures of interest. Keep the conclusion proportional to the analysis: prediction or association may be supported even when causal mediation or surrogate validity is not.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A fitted clinical model has a biomarker coefficient of 0.8 points per biomarker unit, conditional on exposure. State what this coefficient means and one conclusion it does not justify.
- In the confidence-interval example, why are there 14 degrees of freedom rather than 15?
- A biomarker changes with exposure, but the clinical endpoint does not vary across exposure groups. What can and cannot be concluded?
- Why can a strong individual-level biomarker–outcome association fail to establish trial-level surrogacy?
- Name one diagnostic pattern that would make a simple linear biomarker–outcome model questionable.
Solutions. (1) At the same modeled exposure, a one-unit difference in biomarker is associated with an average 0.8-point difference in clinical outcome; it does not show that changing the biomarker causes that difference. (2) A simple linear regression estimates an intercept and a slope, leaving \(n-2=16-2=14\) residual degrees of freedom. (3) The biomarker is exposure-responsive in the observed data, but that alone does not show clinical efficacy or surrogate validity; the analysis also does not establish that no individual-level relationship exists. (4) Individuals may differ in common causes of biomarker and outcome, while trial-level surrogacy asks whether treatment effects on the biomarker predict treatment effects across trials. (5) Curvature in the outcome-versus-biomarker pattern, changing residual spread, or a single highly influential observation could undermine the simple linear model.