Survival After Liver Transplantation in the United States: A Disease-Specific Analysis of the UNOS Database Mark S. Roberts,1,2,3 Derek C. Angus,2,3,4 Cindy L. Bryce,1,2,3 Zdenek Valenta,5 and Lisa Weissfeld2,5 Our goal was to describe disease-specific survival and the clinical variables that predict survival in a large national cohort of adult liver transplant recipients. Data on 17,044 adult patients who received an initial orthotopic liver transplant between 1990 and 1996 with follow-up through 1999 was obtained from the United Network for Organ Sharing (UNOS). Disease-specific Kaplan-Meier survival plots and Cox Proportional Hazards models were estimated, and differences in the clinical characteristics of patients at the time of transplantation by disease were examined. Overall posttransplant survival currently exceeds 85% in the first year and is approaching 75% at 5 years. Unadjusted KaplanMeier survival is improved for recipients who are younger, female, and in better clinical condition. Survival is a function of disease and level of illness: cancer, fulminant liver failure, alcoholic liver disease, and the hepatitidies have the poorest prognosis, while primary billiary cirrohsis and sclerosing cholangitis have the best. Recipients who were outpatients before transplantation have longer survival than those transplanted from the hospital or intensive care unit. Although the model for end-stage liver disease (MELD) score was designed to predict pretransplant survival, patients with higher MELD scores have poorer posttransplant survival, but the MELD score is less predictive than the specific disease. Differences in disease-specific survival are partially explained by differences in disease severity at the time of transplantation. In conclusion, Disease-specific survival models indicate that there remains tremendous variability in survival as a function of underlying liver disease. However, a significant portion of the difference in survival between dis-
Abbreviations: UNOS, United Network for Organ Sharing; PSC, primary sclerosing cholangitis; MELD, model for end-stage liver disease; HCV, hepatitis C. From the 1Section of Decision Sciences and Clinical Systems Modeling, Division of General Medicine, 2Center for Research on Health Care, 3 Department of Medicine, 4Department of Critical Care Medicine, University of Pittsburgh School of Medicine, Pittsburgh, PA; and 5Department of Biostatistics, University of Pittsburgh Graduate School of Public Health, Pittsburgh, PA. Supported by grant (R01 HS09694) from the NIDDK and the AHRQ. Address reprint requests to Mark S. Roberts, MD, MPP, University of Pittsburgh, Center for Research on Health Care, 230 McKee Place, Suite 600, Pittsburgh, PA 15213. Telephone: 412-692-4849; FAX: 412692-4838; E-mail:
[email protected] Copyright © 2004 by the American Association for the Study of Liver Diseases Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/lt.20137
886
eases arises from differences in clinical characteristics at the time of transplantation. (Liver Transpl 2004;10:886–897.)
S
ince its first report in 1963, orthotopic liver transplantation has evolved to become a standard treatment for end-stage liver disease. There are now over 4,000 liver transplant procedures performed each year in the United States.1 However, liver transplantation is expensive and available only to a small portion of all patients who might benefit. This has fueled significant national debate regarding which patients should receive donor organs.2,3 Many factors influence this debate, such as competing moral viewpoints,4 state’s rights, and the efficient use of a scarce resource.5 Uncertainty regarding differences in patient outcomes furthers this concern. Therefore, key to an informed discussion is an updated, nationally representative, and accurate prediction of the likely survival given a patient’s set of clinical characteristics at the time of transplantation. Such estimates do not currently exist at a level of detail sufficient enough to impact selection or allocation decisions. Current national survival estimates, which indicate 1- and 5-year survival ratess of 87.5 and 73.9%, respectively, are not disease-specific, and do not provide estimates of the effect of important clinical covariates on outcome.1 Disease-specific risk-adjusted survival estimates represent single-center experiences,6–8 are based on older data,9,10 or consider single diseases.11 Recent analyses of national databases do not predict long-term or diseasespecific survival.12 The United Network for Organ Sharing (UNOS) maintains a registry of all candidates for organ transplantation listed at transplant centers in the United States. As part of a larger effort to determine the optimal timing of liver transplantation, we developed diseasespecific survival models based on the UNOS registry. These statistical models provide estimates of expected patient survival by disease and clinical characteristics.
Methods Patient Selection Detailed descriptions of the UNOS registry have been published elsewhere.13 Briefly, the registry follows all prospective
Liver Transplantation, Vol 10, No 7 ( July), 2004: pp 886 – 897
Survival After Liver Transplantation
887
Table 1. End-Stage Liver Disease Classification Scheme Ten-Group Classification 1 2 3 4 5 6 7
8 9
10
Category Primary biliary cirrhosis: primary biliary cirrhosis. Primary sclerosing cholangitis: Crohn’s disease; ulcerative colitis; cholangitis with no bowel disease. Alcoholic liver disease: Lae¨nnec’s cirrhosis. Autoimmune disorders: cirrhosis (drug or industrial exposure; cryptogenic; idiopathic); chronic autoimmune hepatitis (etiology unknown; postnecrotic). Hepatitis C and similar infections: postnecrotic cirrhosis (non-A, non-B; type C; type D; types B and C; types B and D); Lae¨nnec’s cirrhosis (postnecrotic; other). Hepatitis B: postnecrotic cirrhosis (HBsAg-positive). Acute hepatic failure: acute hepatic necrosis (drug exposure; hepatitis A; hepatitis B; HBsAg-positive; non-A, non-B hepatitis; hepatitis C; hepatitis D; hepatitis B and C; hepatitis B and D; other acute viral infection); etiology unknown. Cancers: primary liver malignancy (hepatoma; fibrolamellar hepatocellular carcinoma; cholangiocarcinoma; hepatoblastoma; hemangioendothelioma; hemangiosarcoma; angiosarcoma). Metabolic disorders: alpha-1-antitrypsin deficiency; glycogen storage disease type 1; glycogen storage disease type II; hemochromatosis; hemosiderosis; hyperlipidemia type II; homozygous hypercholesterolemia; primary oxalosis or oxaluria; hyperoxaluria; tyrosinemia; Wilson’s disease or other copper disorder; urea cycle disorder; Crigler-Najjar syndrome; Wolman’s disease; protoporphyria; Niemann-Pick disease; abetalipoproteinemia; Gaucher’s disease; Osler-Rendu-Weber disease; carbamoylphosphate synthase deficiency; amyloidosis; Wiskott-Aldrich syndrome. Other liver diseases: cirrhosis (postnecrotic hepatitis A); secondary biliary cirrhosis (Caroli’s disease; choledochal cyst; other); familial cholestasis (Byler disease; other); cholestatic liver disease not listed above; neonatal hepatitis; biliary atresia (extrahepatic biliary atresia; hypoplasia; Alagille syndrome; other); congenital hepatic fibrosis; cystic fibrosis; Budd-Chiari syndrome; benign tumor (hepatic adenoma; polycystic liver disease; other); liver disease induced by total parenteral nutrition or hyperalimentation; graft-versus-host disease; trauma; biliary stricture or stenosis; idiopathic adult ductopenia; unknown.
candidates listed for organ transplantation, documenting any change in status and date of any transplant. The registry records additional clinical information at the time of transplant, as well as information on the donor, and continues to follow the recipient posttransplant. Registry data include standard demographic, clinical, and laboratory information available at the time of listing, as well as the patient’s priority to receive an organ, known as the UNOS status. The UNOS status determines the order in which patients receive donor organs. We received permission from the University of Pittsburgh Internal Review Board to obtain data from UNOS for liver transplant recipients who were transplanted between 1990 and 1996 (n ⫽ 23,791). From this file, we selected all adults (more than 16 years of age) who underwent their first liver transplant. We excluded those who had combination transplant procedures (e.g., combined heart-liver procedure), providing a sample of adult, initial solitary liver transplants (n ⫽ 17,044). The clinical data collected by UNOS changed in April 1994 when several new clinical variables were added to the data collection forms. For several analyses, therefore, only data in the smaller more clinically-rich later dataset are included (n ⫽ 8,172).
Disease Classification The UNOS classification scheme for the etiology of end-stage liver disease includes more than 60 codes, many of which included too few patients for a disease-specific analysis. We asked our National Clinical Oversight Committee, an expert panel of 7 liver transplant physicians and surgeons from transplant centers throughout the country, to collapse the UNOS scheme to a list of 15 or fewer diagnoses. After discussion at a 1-day meeting and subsequent e-mail distribution, we produced a final list of 10 categories into which all UNOS disease categories were mapped (Table 1). We used the diagnosis at time of transplant and assigned each patient to one of the categories. For those without a diagnosis at transplant, we used diagnosis at time of listing (n ⫽ 151 patients, 0.9%). If no diagnosis was available (n ⫽ 36 patients, 0.2%), we assigned the patients to the miscellaneous group, “Other.” For all diagnoses, we considered the UNOS stated diagnosis as accurate: it is not possible to re-evaluate the accuracy of the diagnosis codes in the UNOS dataset.
Selection of Survival Prediction Variables The list of potential predictor variables is long: there are over 100 variables contained in the UNOS registry alone. Many
888
Roberts et al.
other factors not recorded in UNOS, such as Childs-Turcotte-Pugh or acetylketone body ratios, have also been associated with outcome. In addition, for some clinical characteristics, such as human leukocyte antigens-matching and race, there is conflicting evidence of an effect on survival.14 – 16 We therefore employed a 2-stage procedure to select variables for inclusion in the statistical models. First, we asked the National Clinical Oversight Committee to select all variables for which they believed there was either strong evidence or strong theoretical rationale for association with survival. If available in the registry, each of these recommended variables were considered candidate variables. Second, we tested each candidate variable in single variable analyses, using a log-rank test or Cox Proportional Hazards model including only the single candidate variable. All of the variables were tested for significance, and we retained those variables whose association with outcome had a significance level of less than 0.1. This final set of variables that were available, clinically important, and statistically significant in single variable analyses were used to construct the multivariable survival models. Although the model for end-stage liver disease (MELD) score17 was being developed at the time of our variable selection process, it was not included in the analysis as it was designed to predict pretransplant survival. The current analysis is interested in predicting posttransplant survival based on clinical characteristics unrelated to a particular priority scheme. However, because our dataset contained sufficient clinical material to calculate the MELD score at the time of transplant, we were able to stratify posttransplant survival by a calculated MELD score.
Statistical Analyses We described survival overall and by disease, gender, UNOS status, age, year of transplant, and MELD score using standard Kaplan-Meier methods.18 We constructed disease-specific survival estimates using the Cox Proportional Hazards regression model19 and fit the models using S-plus.20 To take advantage of the full sample size, we also built a single stratified model and controlled for disease by including disease category as a stratum indicator. The basic difference between these 2 techniques is that disease-specific models allow the effect of a specific covariate to vary across diseases, whereas a stratified model requires the magnitude of the effect of a covariate to be constant across all diseases and only allows the baseline (intercept) risk of that covariate to change. We selected variables for inclusion in all disease-specific models using a stepwise approach, dropping variables with an insignificant coefficient (P ⱖ .05). We fit the stratified model using the same approach, although the disease stratum indicator always remained in the model. We then used the models to generate survival curves for a “standardized” patient across all disease groups. As KaplanMeier estimates do not control for differences in the severity of illness across diseases at the time of transplantation, this step allows for the comparison across diseases of patients with equal severity of illness.
Results Patient Characteristics The majority of liver transplant recipients were Caucasian, and patients ranged in age from 16 to 77 years (Table 2). The most prevalent diagnoses were hepatitis C (HCV), alcoholic liver disease, and autoimmune disorders, and these accounted for nearly two-thirds of all cases. We restricted our study to patients receiving their first liver transplant and also eliminated patients who had already received a transplant for another organ (less than 1% of the sample). The clinical characteristics of transplanted patients have changed over time. From 1990 to 1996, the average recipient age increased (46.8 vs. 50.0 years, P ⫽ .003), the proportion of women decreased (45.1 vs. 41.5%, P ⫽ .011), and the proportion of Caucasians (83.4 vs. 87.1% P ⫽ .080) and African-Americans (5.5 vs. 7.2%, P ⫽ .003) increased. (Some of this difference may be related to interval changes in the details of racial coding in the UNOS database.) The distribution of disease has changed as well. Consistent with the rising incidence of HCV, that disease now represents the largest group of transplant recipients (16.8 vs. 31.8%, P ⬍ .001), with declines in alcoholic liver disease (22.2 vs. 15.7%) and cancer (6.4 vs. 2.4%) as indications for liver transplantation. Patient Survival The average survival was 83.0% at 1 year, 70.2% at 5 years, and 61.9% at 8 years (figures not shown). Survival was higher in women (83.1, 71.6, and 64.5%) than in men (82.9, 69.1, and 59.8%), and in patients transplanted under the age of 60 (Fig. 1a and b). Location at time of transplant (at home, in hospital, and in intensive care unit), which is an element of UNOS status, is a significant predictor of transplant success, with 1-year survival of 87, 81, and 79% for patients transplanted from home, hospital, or intensive care unit, espectively. Survival has improved over time: 1-year survival rose from 74.8% for those transplanted during 1990 to 86.2% for those transplanted in 1996 (P ⬍ .001). However, the improvement in survival over time has been slowing: the 1-year survival for patients transplanted in 1996 is not statistically significant from the survival of patients transplanted in 1994. Figure 2 displays the Kaplan-Meier plots by disease group. As expected, there is substantial variability in survival by disease, with 5-year survival rates ranging from ⬍40% for patients with cancer to ⬎80% for patients with primary biliary cirrhosis or primary
889
Survival After Liver Transplantation
Table 2. Sample and Subsample Characteristics
Gender Male Female Race White Black Asian Middle Eastern American Indian/Alaskan Native Indian subcontinent Other Unknown Ethnicity non-Hispanic Hispanic Not specified Age at time of transplant (years) Average Minimum Maximum Primary diagnosis Hepatitis C virus Alcoholic liver disease Autoimmune disorders, crypt. Primary biliary cirrhosis Primary sclerosing cholangitis Acute failure (fulminants) Hepatitis B virus Metabolic disorders Cancer Other Prior transplant (not liver)
Entire Sample (n ⫽ 17044)
Subsample 1990 – 1992 (n ⫽ 5857)
Subsample 1993 – 1994 (n ⫽ 5128)
Subsample 1995 – 1996 (n ⫽ 6059)
56.8% 43.2%
55.6% 44.4%
56.8% 43.2%
57.9% 42.1%
83.9% 6.3% 3.0% 1.3%
82.1% 5.5% 3.2% 1.5%
82.2% 6.5% 2.7% 1.4%
87.1% 6.8% 3.2% 1.1%
0.6% 0.2% 0.6% 4.2%
0.4% 0.2% 1.0% 6.2%
0.6% 0.1% 0.8% 5.8%
0.7% 0.2% 0.1% 0.8%
59.8% 9.2% 31.1%
40.6% 7.6% 51.8%
52.7% 8.6% 38.6%
84.3% 11.2% 4.6%
48.8 16.0 77.5
47.9 16.0 76.8
48.8 16.0 77.5
49.7 16.0 77.3
26.2% 18.5% 16.7% 9.4% 8.7% 5.6% 4.6% 3.4% 3.2% 3.6% 0.4%
20.0% 21.3% 16.1% 10.3% 9.2% 5.7% 5.7% 3.6% 4.9% 3.5% 0.5%
28.3% 18.2% 16.8% 10.0% 8.7% 5.9% 3.8% 3.3% 2.4% 2.7% 0.6%
30.8% 16.1% 17.2% 7.9% 8.3% 5.3% 4.3% 3.4% 2.2% 4.5% 0.3%
sclerosing cholangitis. Patients with fulminant liver failure; hepatitis B, HCV, and alcoholic liver disease have intermediate survivals. As expected, patients with fulminant liver disease have the highest initial posttransplant mortality, but their survival parallels the survival of primary biliary cirrohsis and PSC after the first year. Figure 3 illustrates survival based on MELD score at the time of transplantation, with MELD scores grouped into 4 categories of increasing levels of illness. Although the Kaplan-Meier curves are significantly different (logrank tests P ⬍ .000 for all 4 curves), 5-year survival in the highest MELD category (MELD Score ⬎24) is only 7% lower than those patients in the lowest MELD category (MELD score ⬍10). The range of 5-year survival is much greater when stratified by disease, indicat-
ing disease is potentially a stronger predictor of posttransplant survival than MELD score. Predictors of Survival This analysis confirms many of the expected effects of various clinical characteristics on posttransplant survival. Table 3 presents the results of the overall diseasestratified models as well as disease-specific models for each of the 10 disease groups. Survival models for 2 time periods are presented: a “reduced model”, which contains only those clinical variables that were present in the UNOS database throughout the entire study period, and a “full model”, which is estimated on the smaller set of patients after 1994 that contains the more
890
Roberts et al.
Figure 1. Survival after liver transplantation in the United States, 1990 – 1996. Kaplan-Meier plots with confidence intervals of posttransplant survival by gender, age, location, and year of transplant are shown.
complete set of clinical variables. For the disease-specific models, only the direction and strength of significant relationships are shown: the full table with coefficients and standard errors is available from the authors. It is important to note that in the full model, for the years in which all clinical data were available, all of the variables listed in Table 3 were tested for significance, and the variable was retained in the model only if the variable was significant at P :ⱕ.10. Across all disease categories, survival is adversely affected by increasing age of the recipient and the
donor, by the presence of poor renal function, and by declining hepatic function (as measured by increasing bilirubin or declining albumin). Some of these relationships are non-linear: the best fitting model includes age of the donor and recipient as both a standard and squared term, indicating that the deleterious effect of age increases with advancing age of donor and recipient. Controlling for blood type compatibility (ABO) and other disease factors, Caucasian recipients survived longer than minorities, as did recipients who received organs from Caucasian donors.
Survival After Liver Transplantation
891
Figure 2. Disease-specific survival after liver transplantation. The graph depicts Kaplan-Meier actuarial posttransplant survival by individual disease group, which are aggregations of the more than 60 diseases that form the reasons for transplantation in the UNOS database.
In the models estimated from the more clinically complete data (1994-1996), clinical variables such as the presence of encephalopathy and ventilator support adversely affect survival. Several variables that the National Clinical Oversight Committee indicated should be expected to affect outcomes were not found to be significant predictors of survival in the UNOS dataset. For example, although in overall Kaplan-Meier analysis women live longer than men, multivariable analysis found that female gender was not protective when controlling for disease, severity of illness, and other clinical factors.. However, matching of donor and recipient gender increased survival. Body mass index (BMI), prothrombin time, and blood type of the recipient were not significant in either stratified model, although these variables were significant in several specific diseases and were included in the appropriate disease-specific models. As expected, many of the variables related to organ matching and the quality of the donor organ were found to be significant predictors of mortality. Evidence of cytomegalovirus in the donor and the donor race were significant predictors of survival in the
overall models; cytomegalovirus and blood type matching were significant in several individual diseases. Several covariates suggested by the National Clinical Oversight Committee as predictive of survival (ChildPugh, presence of edema, history of ulcerative colitis, and length of stay in the intensive care unit) were not available in the UNOS database for either time period. Predicted Survival by Disease Figure 4 presents estimated predicted survival by a disease group using disease-specific Cox Proportional Hazards models for a cohort of patients with clinical characteristics at the median of our sample (Table 4). This figure compares survival across disease groups by controlling for disease severity at time of transplant: each curve is estimated assuming patients in the different disease categories had the same set of clinical characteristics (age, race, gender, bilirubin, prothrombin time, creatinine, albumin, and donor characteristics) at the time of transplantation. Although overall, the predicted survival for an average patient in the cohort was
892
Roberts et al.
Figure 3. Posttransplant survival by model for end-stage liver disease (MELD) Score. Kaplan-Meier actuarial posttransplant survival is stratified by MELD score. Although overall, the survival by MELD score groups are significantly different (log-rank P < .000), the differentiation arises from the lowest scores compared to the highest. There is no statistical difference between MELD scores of 24.
similar to that found in the Kaplan-Meier plots, there are several diseases that demonstrate differences in predicted survival. The major difference between actual (KaplanMeier) versus predicted (Cox Proportional Hazards) survival is the effect of different clinical characteristics at the time of transplantation detailed in Table 5, which describes the average value of clinical covariates by disease group at the time of transplantation. For example, the predicted 5-year survival for a patient with an underlying diagnosis of cancer and clinical characteristics similar to all patients transplanted is 34%. However,the Kaplan-Meier estimate of survival is 40%, indicating that, controlling for disease severity (as measured by equivalent laboratory values), patients with cancer have an even poorer survival compared to other diagnoses. A similar but opposite effect is seen in acute or fulminant liver failure where the 5-year Kaplan Meier survival is 67%. The predicted survival for patients with “average” disease severity is 72% indicating that, on
average, patients with fulminant liver failure are sicker at the time of transplantation than those patients with chronic liver disease. Fulminants have higher bilirubin, worse renal function, and are more likely to be in the intensive care unit at the time of transplant. The comparison between primary biliary cirrhosis and sclerosing cholangitis is also interesting. The Kaplan Meier estimates of survival for these 2 diseases are virtually identical, with a 5-year survival of 80%. However, predicted survival between the 2 diseases with patients of identical clinical characteristics indicates a substantial survival benefit to patients with primary biliary cirrhosis (5-year predicted survival of 86 vs. 80% of PSC). The only major difference in clinical characteristics is gender: women make up the majority of patients with PSC.
Discussion There are several observations from this study. First, survival after liver transplantation in the United States is
Table 3. Predictors of Posttransplant Survival: 1990 – 1996 Stratified Models (Relative Risk and 95% Limits) Covariates Suggested by Expert Panel Patient and clinical variables Age Age**2 (age squared) Gender (male indicator) Race (white indicator)
RM (1990 – 96)
1.018 1.01 (1.015,1.021) 6 (1.010, 1.021)
Donor age**2 (squared) Age interaction (patient-donor) Donor gender (male indicator) Gender interaction Donor race (white indicator) Donor Hispanic ethnicity Race interaction Donor blood type ABO match Donor anti-CMV Ischemia time (natural log) Cause of death Use of pressors History of diabetes Other variables CMV match Prior transplant (not liver) Donor organ procured locally
HCV
ALD
HBV
Auto
Cancer
Met dis
Acute
Other
RM FM RM FM RM FM RM FM RM FM RM FM RM FM RM FM RM FM RM FM
佢
1
佢
佣
佦
佢
佦
2
0.793 (0.716,0.879)
佦
⫹
1
–
1
佣
佧
0.712 (0.610,0.831) Bilirubin (natural log) 1.078 (1.041,1.117) Creatinine 1.364 (natural log) (1.279,1.455) Prethrombin time (natural log) Risk factors Anti-CMV Ascites Not available Diabetes Not available Encephalopathy Not available
Donor variables Donor age
PSC
佢 佧
Weight (BMI) ABO (blood type) Albumin (natural log)
Hemodialysis Prior RUQ surgery Spontaneous bacterial peritonitis Variceal bleeding (history of) Ventilator Support
FM (1994 – 96)
Disease-Specific Survival Models PBC
0.659 (0.511, 0.849) 1.069 (1.007, 1.136) 1.326 (1.189, 1.479)
1 2
佣
佧 佦
佦
佢
佢
佣 佦
佢
佦 佦
2
佧
1 1
1
1
佦
佦
佦
佦
佢
佢
1
佦
佦 1
1
佢
佦
1
佧 佧 佦 佦
1.225 (1.063, 1.413)
佦
1 1
Not available Not available Not available
佦 1
Not available Not available
1.512 (1.216, 1.879)
1.007 (1.005, 1.010)
1.009 (1.005, 1.012)
佢
1
佦
佢
1
1
佢
佢
佦
–
1
佢
–
佢
佢
1
佢 佦 0.912 (0.846, 0.983) 0.926 (0.862, 0.995) 0.843 (0.762, 0.933)
佧
佧
佧
佧
0.865 (0.770, 0.972) 0.738 (0.631, 0.863)
佦
2
佣
1
佣
佣 佦
佦
佣 佦
佦 佧
2 佦
1.179 (1.043, 1.332)
佧
佢
佦
2
Not Available Not Available 0.923 (0.856, 0.995) 1.598 (1.074, 2.376)
佦
佧 佢
0.879 (0.777, 0.995) 2.075 (1.001, 4.302)
1 佦
2 佦
佢
1 2
2
Abbreviations: RM, reduced model; FM, full model. *The reduced model includes the set of covariates identified by our expert panel that were available through UNOS during the entire period studied (1990 – 1996). In April 1994, UNOS began collecting additional information; this larger set of covariates is included in the full model (1994 – 1996). 佦, coefficient is positive and P ⬍ .05; 1, coefficient is positive and P ⬍ .01; 佢, coefficient is positive and P ⬍ .001; 佧, coefficient is negative and P ⬍ .05; 2, coefficient is negative and P ⬍ .01; 佣, coefficient is negative and P ⬍ .001; ⫹, coefficient is positive but not significant (variable is retained because higher order terms proved significant); – , coefficient is negative but not significant (variable is retained because higher order terms proved significant).
894
Roberts et al.
Figure 4. Cox Proportional Hazards prediction of posttransplant survival after liver transplantation by disease group. This figure depicts predicted posttransplant survival for all 10 disease groups assuming that the individuals are transplanted with identical pretransplant clinical characteristics and identical donor characteristics: the only difference is the underlying etiology of liver disease. This emphasizes the disease-specific differences in survival, controlling for the empiric observation that patients with different diseases present at different levels of illness and have different rates of progression.
excellent and has been improving over time, even with a changing distribution of the diseases leading to transplantation. Survival rates of 1 year for all but a few diseases exceed 90%, and long-term survival (8 years) approaches 70%. The epidemiology of liver disease as a cause of transplantation is rapidly changing and is becoming more weighted by HCV, which has a poorer survival relative to many other diseases. However, this relationship is complex as patients with HCV are different, from say, patients with primary billiary cirrohsis/PSC. Their clinical presentation at time of transplant is different, and they are less likely to be in the hospital at time of transplant. This work emphasizes the importance of considering not only the predictors of transplant success and failure, but also the constellation of characteristics and severity of illness with which patients present for transplant evaluation across different diseases.
Table 4. Clinical Characteristics Used to Compare Survival by Disease Using Cox Proportional Hazards Models Recipient data Age Gender Race Bilirubin (mg/dL) Albumin (g/dL) Creatinine (mg/dL) Prothrombin time (sec) Donor data Age Gender Race
49.3 years Male Caucasian 3.4 2.9 1.0 14.8 31.2 years Male Caucasian
895
Survival After Liver Transplantation
Table 5. Mean Values of Clinical Covariates at Time of Transplantation by Disease Group
Patient characteristics Age (years) Male (%) Caucasian (%) African American (%) Hispanic (%) Other race (%) Weight (BMI)
ABO (%) Type A Type B Type AB Type O Albumin (g/dL)
Total Bilirubin (mg/dL)
Creatinine (mg/dL)
Prothrombin time (sec)
Anti-CMV (% positive) In hospital (%) In ICU (%) Donor characteristics Age (yrs)
Male (%) Anti-CMV (% positive) Cold ischemia time (in minutes)
Overall (17,044)
PBC (1,596)
PSC (1,487)
HCV (4,460)
ALD (3,150)
HBV (784)
Auto (2,839)
Cancer (538)
Met dis (583)
Acute (954)
Other* (653)
48.8 (11.6) 56.8 79.1 6.1 9.2 5.5 26.4 (5.4) [15,328]
54.2 (9.4) 11.0 83.6 3.1 9.4 3.8 24.3 (4.8) [1,447]
46.7 (11.9) 65.6 85.7 8.9 3.4 2.0 24.2 (4.4) [1,364]
49.5 (9.7) 65.9 77.6 5.3 11.2 5.8 27.6 (5.3) [4,006]
50.5 (9.0) 77.6 86.2 3.4 7.9 2.5 26.5 (5.0) [2,872]
48.3 (10.2) 81.6 63.5 4.0 6.3 26.3 25.6 (4.8) [696]
48.8 (13.4) 39.9 77.8 8.2 10.3 3.8 27.8 (6.0) [2,576]
51.2 (11.8) 60.6 67.8 4.8 9.7 17.7 25.0 (4.6) [486]
44.7 (14.9) 65.5 86.1 1.7 8.7 3.4 26.6 (5.6) [534]
41.0 (13.8) 41.2 66.2 15.1 12.1 6.6 26.6 (5.6) [834]
41.1 (14.2) 40.6 75.8 11.3 9.2 3.7 24.3 (5.0) [513]
41.2 11.8 4.6 42.4 3.0 (0.7) [15,939]
41.6 11.2 3.9 43.2 3.1 (0.7) [1,498]
42.0 11.6 3.8 42.6 3.1 (0.7) [1,403]
41.3 11.7 5.1 41.9 2.9 (0.6) [4,162]
43.0 11.3 4.6 41.2 3.0 (0.7) [2,965]
38.2 12.9 6.0 42.9 2.9 (0.6) [750]
41.2 12.0 4.3 42.5 2.9 (0.7) [2,667]
38.7 15.2 4.8 41.3 3.3 (0.7) [508]
44.9 7.9 4.6 42.5 2.9 (0.7) [552]
36.6 13.5 4.6 45.3 2.9 (0.6) [881]
39.5 13.2 3.4 44.0 3.2 (0.8) [553]
7.6 (10.3) [16,709] 1.3 (1.1) [16,749] 16.2 (5.6) [16,165]
9.0 (9.8) [1,575] 1.0 (0.7) [1,576] 14.1 (4.0) [1,529]
9.5 (10.3) [1,465] 1.1 (1.0) [1,469] 14.2 (4.2) [1,424]
5.1 (8.0) [4,357] 1.2 (0.9) [4,377] 16.0 (4.5) [4,233]
6.6 (8.9) [3,107] 1.5 (1.3) [3,111] 16.3 (4.7) [2,993]
7.6 (10.4) [777] 1.4 (1.2) [779] 17.2 (5.0) [753]
6.8 (9.5) [2,797] 1.3 (1.0) [2,797] 16.4 (5.1) [2,693]
4.3 (7.6) [534] 1.2 (1.2) [535] 14.3 (3.8) [514]
10.1 (14.7) [576] 1.3 (1.1) [576] 17.8 (6.9) [557]
19.5 (14.0) [940] 1.9 (1.9) [943] 23.1 (10.4) [911]
8.3 (10.8) [586] 1.3 (1.3) [581] 16.0 (5.6) [558]
67.9 [14,717] 22.6 16.7
71.5 [1,409] 18.3 10.5
56.4 [1,319] 20.6 8.5
71.9 [3,843] 23.4 10.6
67.3 [2,702] 23.5 17.8
79.2 [664] 27.3 19.4
66.9 [2,490] 24.0 13.5
61.9 [451] 26.0 13.2
61.3 [514] 24.4 22.5
69.0 [813] 13.0 66.9
60.7 [512] 25.9 21.2
33.8 (15.9) [16,485] 63.3
32.7 (16.4) [1,558] 54.2
32.5 (15.3) [1,429] 65.8
34.9 (15.9) [4,280] 64.4
34.1 (15.7) [3,071] 66.3
33.7 (15.8) [756] 66.7
33.2 (15.8) [2,749] 61.3
34.2 (16.4) [527] 66.2
33.7 (16.0) [569] 65.9
33.8 (16.4) [920] 61.9
32.2 (16.5) [626] 59.6
57.8 615.2 (248.7) [15,481]
58.3 590.0 (250.0) [1,447]
54.6 608.9 (254.4) [1,355]
60.4 611.4 (244.6) [4,023]
57.5 637.3 (247.3) [2,886]
55.0 631.0 (263.5) [725]
57.1 614.9 (245.6) [2,606]
55.9 620.0 (247.1) 482
58.1 602.8 (241.9) [541]
56.4 591.8 (257.5) [876]
57.5 634.1 (246.7) [540]
Standard deviations are in parentheses. The number of non-missing observations (if different from the total) is listed in square brackets.
Second, this work represents a comprehensive set of disease-specific models estimated from a large representative national database. The results underscore the need for disease-specific models because we have found that survival is strongly linked to disease and that different clinical variables affect future survival in different diseases, which may be important in timing and treatment decisions. Third, our analysis confirms the findings of prior analyses indicating that the distribution of the timing of deaths is quite different between fulminant disease and chronic disease, with many deaths occurring early (within the first 6 months posttransplant) in acute liver disease. By extension, in chronic liver disease, most
deaths across all diseases occur late, which indicates efforts focusing on events in this period may therefore be most likely to improve survival in chronic liver disease and efforts to improve early survival may have more impact in acute liver disease. Fourth, the MELD score does not differentiate posttransplant survival well and is a poorer predictor than the underlying etiology of liver disease. Although this is not surprising (the MELD score was designed to predict pretransplant survival), it emphasizes that the characteristics that are predictive of survival while waiting for transplantation are not identical to those that predict posttransplant survival. It is interesting to note that this disease-specific
896
Roberts et al.
multivariate analysis indicates that a different set of clinical variables are significant predictors of posttransplant survival than has been previously reported for certain diseases. For example, in the classic descriptions of posttransplant survival for primary billiary cirrohsis and PSC described by Wiesner,8 and Grambsch21 age, bilirubin, albumin, prothrombin time, and edema are found to be significant predictors of survival. However, in the larger national UNOS sample, albumin and bilirubin are not found to be statistical predictors of posttransplant survival. It seems likely that outcome has improved since the early Mayo Clinic work with FK-506,22 maturation of surgical techniques, and widespread use of a standardized approach to cytomegalovirus prophylaxis, rejection, etc. However, these differences indicate that our understanding of the predictors of transplant success and failure remains dynamic, and that the re-evaluation of risk factors and survival predictors is an important ongoing activity. These results are also interesting when compared to international results from similar time periods. In an analysis of 1,656 patients with follow-up through 1995, Smits and colleagues found the same increasing survival over time and found donor and recipient age, blood type compatibility, and several other similar variables to be predictors of survival.23 However, overall, survival in Europe was slightly poorer, e.g., 5-year survival for patients with cirrhosis in the European experience was only 54 vs. 80% in the US experience. However, data from a more recent review from the European Liver Transplant Registry (which includes 1,821 patients transplanted up to1994) found outcomes to be closer to the US experience.24 The 5-year posttransplant survival for patients with cirrhosis was 70% and 55% for those with acute hepatic failure The US experience is 80 and 67%, respectively. Another important observation is that the current clinical variables contained within the UNOS dataset do not contain all of the variables that a national group of experts agreed were clinically important in the prediction of future survival. Further work should include collection of these variables in a sample large enough to permit an appraising of their predictive value. Such work is crucial as the national registry is currently under overhaul and identification of key variables will be essential if improved clinically acceptable models are to be built in the future. Although many variables that are clinically assumed to be important in survival may not reach statistical significance in multivariable models, the ability to test the effect of clinically relevant variables will increase the face value of statistical prediction.
Finally, the analysis presented here shows important findings for the use of disease specific models in the prediction of survival. Our work indicates that diseasespecific models provide different estimates of survival than overall survival models with disease strata indicators. As the distribution of liver disease etiology changes over time (such as the dramatic rise in HCV as the etiology of end-state liver disease), it will become more important to have accurate disease-specific models of survival to predict the overall benefit from transplantation programs. In addition, the results presented here indicate that a significant portion of the difference in survival by disease demonstrated in standard KaplanMeier plots is dependent upon differences in the clinical characteristics of the patients who present with those diseases, rather than an effect of the disease itself. This may have implications on long-term survival if the timing of transplant evaluation and listing is changed for different diseases. An understanding of the effect of clinical covariates on posttransplant survival is also directly applicable to the problem of the timing of living-related donor transplants. In this case, more complete information is available before transplantation (the characteristics of the donor are known, for example), and the major decision to be made is one of timing. Given the differences found in the clinical variables that predict future survival, it is likely that the optimal timing of living related donor transplants is different in different diseases. In summary, analysis of the UNOS database provides the most current and comprehensive analysis of disease-specific survival after transplantation that is possible in the United States today. This analysis indicates that there are significant differences in survival by type of disease, and that the distribution of diseases is rapidly changing over time. The continuing analysis of the factors that predict transplant success and failure is necessary to predict the overall effect of changes in evaluation, timing, and allocation criteria on overall survival from transplantation. Although survival has been improving, this does not imply that current transplantation times are optimal for patients: it remains possible that transplantation at different times with more favorable clinical characteristics could improve overall survival.
References 1. Annual report of the U.S. scientific registry for transplant recipients and the organ procurement and transplant network: transplant data: 1989 – 1998. U.S. Department of Health and Human Services, Health Resources and Services Administration,
Survival After Liver Transplantation
2.
3. 4. 5.
6.
7.
8.
9.
10.
11.
12.
13.
Office of Special Programs, Division of Transplantation, Rockville, MD; UNOS, Richmond, VA. Ubel PA, Loewenstein G. Public perceptions of the importance of prognosis in allocating transplantable livers to children. Med Decis Making 1996;16:234 – 241. McMaster P. Transplantation for alcoholic liver disease in an era of organ shortage. Lancet 2000;355:424 – 425. Gibbons RD, Duan N, Meltzer D. Inequities in liver transplant allocation. Science 2000;289:549 – 550. Ubel PA, Arnold RM, Caplan AL. Rationing failure. The ethical lessons of the retransplantation of scarce vital organs. JAMA 1993;270:2469 – 2474. Jain A, Reyes J, Kashyap R, Dodson SF, Demetris AJ, Ruppert K, et al. Long-term survival after liver transplantation in 4,000 consecutive patients at a single center. Ann Surg 2000;232:490 – 500. Shmueli D, Mor E, Sharabani E, Bar NN, Ben Ari Z, Yussim A, et al. A 4-year small-center experience in liver transplantation. Transplant Proc 1997;29:2868 – 2869. Wiesner RH. Liver transplantation for primary biliary cirrhosis and primary sclerosing cholangitis: predicting outcomes with natural history models. Mayo Clin Proc 1998;73:575 – 588. Jain AB, Fung JJ, Todo S, Reyes J, Selby R, Irish W, et al. One thousand consecutive primary orthotopic liver transplants under FK 506: survival and adverse events. Transplant Proc 1995;27: 1099 – 1104. Buckels JA. Liver transplantation: current status, complications and prevention. J Antimicrob Chemother 1995;36(suppl B):39 – 49. Forman LM, Lewis JD, Berlin JA, Feldman HI, Lucey MR. The association between hepatitis C infection and survival after orthotopic liver transplantation. Gastroenterology 2002;122:889–896. Smith CM, Davies DB, McBride MA. Liver transplantation in the United States: a report from the UNOS liver transplant registry. Clin Transpl 1999;23 – 34. Harper AM, Baker AS. The UNOS OPTN waiting list: 1988 – 1995. Clin Transpl 1995;69 – 84.
897
14. Donaldson PT, Thomson LJ, Heads A, Underhill JA, Vaughan RW, Rolando N, et al. IgG donor-specific crossmatches are not associated with graft rejection or poor graft survival after liver transplantation. An assessment by cytoxicity and flow cytometry. Transplantation 1995;60:1016 – 1023. 15. Pillay P, Van Thiel DH, Gavaler JS, Starzl TE. Effect of race upon organ donation and recipient survival in liver transplantation 18. Dig Dis Sci 1990;35:1391 – 1396. 16. Nikaein A, Backman L, Jennings L, Levy MF, Goldstein R, Gonwa T, et al. HLA compatibility and liver transplant outcome. Improved patient survival by HLA and cross-matching. Transplantation 1994;58:786 – 792. 17. Wiesner RH, McDiarmid SV, Kamath PS, Edwards EB, Malinchoc M, Kremers WK, et al. MELD and PELD: application of survival models to liver allocation. Liver Transpl 2001;7:567 – 580. 18. Harris EK, Abert A. Volume 114: Statistics: textbook and monographs series. Cohen, DB, ed. New York: Marcel Dekker, Inc. 1991. 19. Cox DR. Regression models and life-tables (with discussion). J Roya Stat Soc Ser B 1972;34:187 – 202. 20. S-PLUS 2000 Guide to Statistics, Volume 1, Data Analysis Products Division, MathSoft, Seattle, WA, 1999. 21. Grambsch PM, Dickson ER, Wiesner RH, Langworthy A. Application of the Mayo primary biliary cirrhosis survival model to Mayo liver transplant patients. Mayo Clin Proc 1989;64: 699 – 704. 22. Fung JJ, Starzl TE. FK506 in solid organ transplantation. Ther Drug Monit 1995;17:592 – 595. 23. Smits J, Meester JD, Persijn GG, Claas FHJ, Vanrenterghem Y. Long term results of solid organ transplantation. In: Cecka JM, Terasaki PI, editors. Clinical Transplant 1996. Report from the Eurotransplant international foundation. Los Angeles, CA: UCLA Tissue Typing Laboratory. 1996:109 – 127. 24. Bismuth H, Samuel D, Castaing D, Williams R, Pereira SP. Liver transplantation in Europe for patients with acute liver failure. Semin Liver Dis 1996;16:415 – 425.