| ||||||||||||||
|
|
|||||||||||||



From the *Anesthesiology Service, VA Medical Center and Department of Anesthesiology, OR Health and Science University, Portland, OR;
Department of Public Health and Preventive Medicine, OR Health and Science University, Portland, OR;
Department of Mathematics and Statistics, Portland State University and Department of Public Health and Preventive Medicine, OR Health and Science University, Portland, OR;
Anesthesiology Service, VA Medical Center, Portland, OR; and ||Surgical Service, VA Medical Center and Department of Surgery, Saint Louis University Medical School, St. Louis, MO.
Address correspondence to Rose Christopherson, MD, PhD, Anesthesiology Service (P3ANES), VA Medical Center, 3710 SW US Veterans Hospital Road, Portland, OR 97229. Address e-mail to rose.christopherson{at}med.va.gov.
| Abstract |
|---|
|
|
|---|
METHODS: We studied long-term survival after resection of colon cancer in a trial of general anesthesia with and without epidural anesthesia and analgesia supplementation for resection of colon cancer in Veterans Affairs Cooperative Study No. 345. Cox and log-normal survival models were used to test the effects of pathological stage, type of anesthesia and other covariates on survival in 177 patients.
RESULTS: The presence of distant metastases had the greatest effect on survival. Thus, analyses were performed separately for patients with and without metastases. For those without metastasis, the hazard ratio for the treatment effects changed at 1.46 years. Before 1.46 years, epidural supplementation was associated with improved survival (P = 0.012), while later, the type of anesthesia did not appear to affect survival (P = 0.27). Hypertension was associated with poorer survival (P = 0.029), as was alcoholism in patients who received epidural anesthesia (P = 0.014). Survival of patients with metastases was unaffected by type of anesthesia. There was a significant age by hypertension interaction (P = 0.002). Patients survived longer if they were hypertensive, but had reduced survival if they were older than 66 years and hypertensive.
CONCLUSION: Epidural supplementation was associated with enhanced survival among patients without metastases before 1.46 years. Epidural anesthesia had no effect on survival of patients with metastases. Additional studies to confirm or refute these findings are warranted.
| Introduction |
|---|
|
|
|---|
| METHODS |
|---|
|
|
|---|
In CSP 345, patients randomized to epidural supplementation received 0.5% bupivacaine with epinephrine 1:200,000 via a lumbar or thoracic epidural catheter. A level of anesthesia of at least T-6 was established before inducing general anesthesia. General anesthesia was maintained in both groups with isoflurane, nitrous oxide, vecuronium, and fentanyl. In the group randomized to epidural supplementation, boluses of 5–10 mL of 0.5% bupivacaine with epinephrine were given as needed. The catheter was used as long as judged clinically appropriate after surgery for pain management, without a specific protocol. Epidural morphine and other drugs were used at the discretion of the local clinicians. Patients randomized to unsupplemented general anesthesia received IV opioids for postoperative pain.1
Long-Term Survival
We obtained approval from the IRBs of all the study centers willing and able to participate. This approval included authorization to perform the study without recontacting the original study patients to obtain their consents, to collect the names, social security numbers, and pathologic tumor stage of the patients enrolled in the trial, and to use this information to determine how long patients survived after surgery.
Long-term survival of CSP 345 patients who had undergone surgery for colon cancer was obtained by querying the VA Beneficiary Information and Records Locator System (BIRLS).3,4 This database is considered to have a very high degree of accuracy because it is based on the payment of burial benefits when a veteran dies. Names, social security numbers, and birth dates, used as identifiers for the BIRLS system, were obtained from patients original consent forms and records. Survival follow-up ended in December 2002. Survival was calculated as the number of days from the date of randomization in the study to the date of death or last contact with the patient and was converted to years by dividing by 365.25.
The electronic database from the original study was used to determine baseline variables and risk factors for survival, type of anesthesia, and date of surgery. Tumor-Nodes-Metastasis staging was obtained from the tumor registries at the participating hospitals. Tumor-Nodes-Metastasis staging was translated into stages 0–IV using the pathologic staging system that was in use at the time the patients had their surgery.5
Statistical Analysis
Student's t-test and
2 contingency tables were used to compare baseline characteristics between the two treatments. P values were not adjusted for multiple comparisons among these characteristics. Univariate survival analyses were performed using the Kaplan-Meier survival procedure.6 The Tarone-Ware procedure7 was used to compare the effects of treatment on survival stratified by metastasis status. Multivariate survival analyses were performed for the nonmetastasis subgroup using Cox regression.8 The log-normal regression model was used to analyze the data from the metastasis subgroup.
By way of explanation for the nonmetastasis subgroup, in survival analysis, the hazard rate, commonly called the risk and designated as h(t), is the probability of dying per unit time given survival to a given time point. This is also called the force of mortality or the time-specific death rate. A critical assumption for the Cox regression model is the proportionality of the hazard (PH) rates over time. We used the Grambsch and Therneau procedure to test this assumption.9 If the assumption is true, the survival curves will not cross over time. If the curves cross, the PH assumption is violated, as is the case for the nonmetastasis subgroup in this study. In order to analyze survival data that violate the PH assumption, a more generalized extended Cox model was used. This procedure, described in Klein and Moeschberger,10 analyzes the survival curves over two time intervals in which the required PH assumption is satisfied over each interval. The mathematical model for the hazard function at time t, excluding the other covariates that satisfy the proportional hazards assumption, is:
|
|
where h0(t) is the unspecified baseline hazard function; ET1 = TRT if t < t0, 0 otherwise, and ET2 = TRT
if t
t0, 0 otherwise; and
1 and
2 are the regression coefficients for treatment (TRT) over each of the two time intervals. Visually, t0 is the point at which the slopes of the survival curves differ because of the change in the hazard ratios (HR).
The HR is now expressed as
|
|
where TRT = 1 for unsupplemented and TRT = 0 for epidural-supplemented general anesthesia.
Estimates of
1,
2, and the regression coefficients for the covariates are obtained by maximum likelihood procedures.11 To determine the optimal change point t0 a profile likelihood procedure12,13 was used as follows: Iterate over each of the n1 plus n2 observed uncensored survival times as possible values of t0. For each observed value, compute the value of the likelihood function with the maximum likelihood estimators of
1,
2, and the regression coefficients for the covariates at that time point from the extended Cox model. Repeating this procedure for each uncensored survival time provides a graph of the likelihood function against time. The point where the graph is at the maximum is the optimal change point t0.
Since the treatment variable in this study was found to violate the PH assumption, model building was performed by initially stratifying on type of anesthesia and considering the following covariates with their two-way interactions: pathological cancer staging/ metastasis, age at randomization, presence or absence of hypertension (as indicated in the medical record), systolic blood pressure, diastolic blood pressure, hematocrit, race, presence of chronic obstructive pulmonary disease, alcoholism, hypertension, diabetes, and whether the patient had ever smoked. These baseline characteristics were taken from the CSP 345 database. Akaike's information criterion (AIC) was used to select the variables.14 This selection procedure balances the model fit against the number of covariates entered. The effect of the change in the hazard ratios at t0 is observed subsequently when the survival curves cross. The extended Cox model procedure produces estimates of the treatment effects represented by the two regression coefficients,
1 and
2, and tests for the significance of anesthesia effects, cancer staging, and the other covariate effects. A plot of the log ratio of the survival probabilities was used to visually demonstrate the difference between epidural-supplemented and unsupplemented general anesthesia over time. Pointwise 95% confidence intervals (CI) were constructed using the bootstrap resampling procedure.15
The analysis of the survival data in the metastasis subgroup consisted of fitting the log of the survival time against the candidate covariates. The procedure maximizes the likelihood function, which considers the right-censored survival times. The S-PLUS16 statistical analysis package was used to perform the statistical calculations and plot the survival graphs for the nonmetastasis and metastasis subgroups.
| RESULTS |
|---|
|
|
|---|
Baseline Characteristics and Pathological Staging
Baseline characteristics of patients who received epidural anesthesia and those who did not were comparable (Table 1). The distributions across treatment groups were comparable for Stages 0–II, but a larger proportion of patients with Stage III or IV disease were randomized to unsupplemented general anesthesia (Table 2; P = 0.029, Kruskal-Wallis test for ordered categories17). Forty-one (44.5%) of the 92 patients in the unsupplemented group had stage III or IV disease whereas 24 (28.2%) of the 85 patients in the epidural-supplemented group had stage III or IV disease. The staging categories were collapsed into nonmetastasis (stages 0–II, n = 112) and metastasis (stages III or IV, n = 65) for the survival analysis.
|
|
Kaplan-Meier Survival Analyses
Kaplan-Meier survival analysis showed a highly significant difference between patients without pathological evidence of metastasis (Stages 0–II) versus those with metastasis (Stages III–IV), combined across the two treatment arms. As expected, patients without evidence of metastasis survived significantly longer than those with metastasis; the median survival time was 6.14 years with 95% CI [5.22, 7.99] for patients without metastasis versus 2.01 years with 95% CI [1.55, 2.47] for patients with metastases (Tarone-Ware: P < 0.0001, Fig. 1).
|
Survival was further examined by treatment within the two metastasis subgroups using the Kaplan-Meier procedure stratified on metastasis status. Patients without metastasis assigned to receive epidural anesthesia exhibited better survival initially. However, the survival curves for the 2 treatments merge at approximately 2.5 years and are coincident until 4 years. At that point, the curves separate, with the group without epidural anesthesia/analgesia trending towards improved long-term survival (Fig. 1). In patients with metastasis, the survival curves for the two treatments coincide until about 1 year, after which they separate, and patients without epidural anesthesia trend towards better long-term survival. Further analyses were performed within the metastasis subgroups to determine whether the survival differences between the treatments were statistically significant.
Survival Analyses by Metastasis Subgroup
Nonmetastasis
The crossing of the survival curves in the nonmetastasis subgroup violates the PH assumption. In order to compensate for this violation and to assess the effects of a set of prognostic factors (covariates) on long-term survival, the extended Cox model described in the Methods section was used.
From the 11 covariates and their two-way interactions, the AIC stepwise procedure, performed on a Cox PH model stratified on treatment, selected hypertension and the interaction of alcoholism with type of anesthesia. These covariates were included as predictors of survival in the extended Cox model in addition to the treatment effects for the first (t < t0) and second (t > t0) time periods. The estimates and P values for the coefficients are reported in Table 3. The overall Cox model is highly significant (likelihood ratio P = 0.002) and the PH assumption is well satisfied over the two intervals, as determined by the Grambsch-Therneau test and plots of the Schoenfeld residuals (not shown). Figure 2 shows the fitted extended Cox model with the hypertension covariate set to 0.64 and 0.47, the proportions of patients with hypertension in the epidural-supplemented and unsupplemented groups, respectively. The treatment by alcoholism interaction for the group with epidural anesthesia was set to its mean of 0.26 (the proportion of these patients classified as alcoholics). Since the treatment by alcoholism interaction for the group without epidural anesthesia/analgesia was not significant and the estimated coefficients and standard errors for the other parameters were essentially the same with or without this interaction, it was removed from the model. The resulting extended Cox model fits the Kaplan-Meier survival curves.
|
|
The slope change point t0, determined by where the profile log-likelihood is maximized, is 1.46 years. This results in later crossing of the survival curves at about 3.2 years. For the reduced model (without the treatment by alcoholism interaction for patients without epidural supplementation; table not shown) the risk of dying in the first 1.46 years for patients who did not receive epidural supplementation is 4.65 times the risk for those who did (P = 0.012, 95% CI [1.40, 15.42]). For patients who survived beyond 1.46 years, not receiving epidural anesthesia appears to provide a lower risk of mortality, 0.71, but the benefit is not statistically significant (P = 0.260, 95% CI [0.39, 1.29]). Ninety-six (86%) of the 112 patients without metastasis survived beyond 1.46 years. Over both time periods, the risk of dying for patients having hypertension was 1.79 times higher than the risk for those who did not have it (P = 0.029, 95% CI [1.08, 2.96]). The risk of mortality for alcoholics who received epidurals was 2.33 times the risk for nonalcoholics receiving epidural anesthesia (P = 0.014, 95% CI [1.19, 4.56]) over both time periods.
Figure 3 shows the plot of the natural log of the ratio of the survival probabilities (equivalently, cumulative hazard rate of epidural-supplemented anesthesia minus the cumulative hazard rate of unsupplemented anesthesia) with pointwise 95% bootstrap confidence limits. The upper 95% confidence limit is below zero before 1.5 years. It is interesting to note that the time point at which this 95% CI crosses zero is essentially the same as the cut point, 1.48 years, obtained from the extended Cox model. Thus, patients who received epidural supplementation had better early survival. However the benefit was lost after 1.5–2.0 years.
|
Metastasis
The survival curves for the metastasis subgroup show satisfactory compliance with the proportional hazards assumption. Several survival models were fit, but the log-normal model fit the Kaplan-Meier curves most closely, as it predicted the natural log of the survival time, t, by a linear combination of the treatment and other covariates. In addition to epidural supplementation, the stepwise AIC procedure selected age (P = 0.059), hypertension (P = 0.003), and the age by hypertension interaction (P = 0.002) (Table 4). The resulting summary survival curves from the log-normal model, with age set to 69 (mean age) for the epidural anesthesia/analgesia group, to 68 (mean age) for the unsupplemented group, and the hypertension covariate set to 0.54 and 0.46 (the proportion of hypertensive patients for the epidural-supplemented and unsupplemented groups, respectively), fit the Kaplan-Meier survival curves very well (Fig. 4).
|
|
Type of anesthesia did not have a statistically significant effect (P = 0.22), and when it was removed from the model the coefficients and the P values for the remaining covariates change very little. The full model shown in Table 4 indicates that advanced age and the diagnosis of hypertension are associated with increased survival time among patients with metastases. However, if a patient is age 66 years or more and hypertensive, the survival time is decreased because of the negative coefficient for the age by hypertension interaction (–0.119). For example, from the log-linear model, a patient receiving epidural anesthesia/ analgesia who was 68 years old at the time of randomization without hypertension would be expected to live 2.48 years whereas the same patient with hypertension would be expected to live 2.01 years.
| DISCUSSION |
|---|
|
|
|---|
Among patients without metastases, hypertension had an adverse effect on survival, potentially due to the well-defined adverse effects of hypertension on the cardiovascular system. Alcoholism had an adverse effect on survival among patients without metastases who received epidural anesthesia, but no effect on survival among those who received unsupplemented general anesthesia. This may have been due to the presence of alcoholic liver disease, including cirrhosis, in some alcoholic patients. Likewise, cell-mediated immunity, specifically lymphocyte function, particularly CD8+ cells and natural killer cells, are believed to be important in preventing cancer metastasis but are reduced in humans after epidural-supplemented general anesthesia compared to patients receiving general anesthesia.22,23 The suppression of the natural immune response added to the hemodynamic response to epidural anesthesia could compound the effects on the liver, which may already be compromised in alcoholics.
Among patients with lymphatic or more distant metastases, the early survival benefit attributable to epidural anesthesia was not found, possibly because it was masked by the high overall death rate of these patients (Fig. 1). Hypertension was associated with improved survival regardless of type of anesthesia (P = 0.006). The most frail patients with metastatic cancer may have had low blood pressures and poor survival. The interaction of age and hypertension (P = 0.004), such that elderly hypertensive patients did not survive as long, may have been due simply to the fact that, over time, hypertension has a relentless adverse effect on the cardiovascular system. Among patients who were not hypertensive, there was a trend (P = 0.087) for older patients with metastatic disease to survive longer than younger patients. There is controversy as to whether younger patients may have more aggressive tumors and/or poorer survival than older patients.24,25
Our investigation should be viewed as a preliminary exploration of the effect of epidural anesthesia on long-term survival after cancer surgery. It is unfortunate that we do not have data on the cause of death of our patients. Cause of death might have been different in the early postoperative period, during which epidural anesthesia was associated with a significant increase in survival, from causes of death in later years, when there was no benefit associated with epidural anesthesia. Knowledge of the cause of death may have shed some light on the reduced survival probability in all time periods of alcoholic patients without metastases who received epidural anesthesia. Further studies are needed to confirm these findings. Data on other types of cancer, other patient populations, adjuvant chemotherapy or radiation, and time and location of recurrence would also be important. Accordingly, we offer our findings as an impetus for future research in this area.
| ACKNOWLEDGMENTS |
|---|
| Footnotes |
|---|
Supported by a Department of Veterans Affairs Merit Review Epidemiology Grant and by funds from the Department of Anesthesiology & Perioperative Medicine, OR Health & Science University.
None of the authors has financial interests that would be affected by the outcome of this study.
Reprints will not be available from the authors.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|