Research | Open | Published:
The prognostic role of systemic inflammatory markers on HIV-infected patients with non-Hodgkin lymphoma, a multicenter cohort study
Journal of Translational Medicinevolume 13, Article number: 89 (2015)
The systemic inflammatory response has been postulated as having prognostic significance in a wide range of different cancer types. We aimed to assess the prognostic role of inflammatory markers on survival in HIV-infected patients with Non-Hodgkin Lymphoma (NHL), and to compute a prognostic score based on inflammatory biomarkers.
We evaluated data on HIV patients with NLH diagnosis between 1998 and 2012 in a HIV Italian Cohort. Using Cox proportional regression model, we assessed the prognostic role of Neutrophil-Lymphocyte Ratio (NLR), Platelet-Lymphocyte Ratio (PLR), Glasgow Prognostic Score (GPS), modified Glasgow Prognostic Score (mGPS), Prognostic Index (PI), and Prognostic Nutritional Index (PNI). We also computed a risk score equation, assigning patients to a derivation and a validation sample. The area under the curve (AUC) was use to evaluate the predictive ability of this score.
215 non-Hodgkin lymphoma cases (80.0% males) with a mean age of 43.2 years were included. Deaths were observed in 98 (45.6%) patients during a median follow up of 5 years. GPS, mGPS, PI and PNI were independently associated with risk of death. We also computed a mortality risk score which included PNI and occurrence of an AIDS event within six months from NHL diagnosis. The AUCs were 0.69 (95% CI 0.58 to 0.81) and 0.69 (95% CI 0.57 to 0.81) at 3 and 5 years of the follow-up, respectively.
GPS, mGPS, PI and PNI are independent prognostic factors for survival of HIV patients with NHL.
Non-Hodgkin Lymphoma (NHL) is one of the most common cancer and cause of death among HIV-infected patients [1-3]. Overall survival is poor and more than half subjects die within five years from NHL diagnosis, although it has improved after the introduction of the combined Antiretroviral Therapy (cART) [2-9].
The International Prognostic Index (IPI) was developed at the beginning of 1990s by the International Non-Hodgkin Lymphoma Prognostic Factors Project to create a prognostic tool for patients with aggressive NHL treated with doxorubicin-containing chemotherapy in the general population . The IPI categorizes patients into low, intermediate and high risk groups based on the baseline characteristics including age, Eastern Cooperative Oncology Group performance status, Lactate Dehydrogenase (LDH) level, Ann Arbor stage, and extranodal involvement. The IPI has also been demonstrated prognostic for lymphomas associated with HIV infection .
Tumor microenvironment and inflammatory response have an important role at different stages of tumor development, including initiation, promotion, malignant conversion, invasion and metastasis . Recent studies showed that some blood parameters associated with systemic inflammation, i.e. C-Reactive Protein (CRP), white cell, neutrophil and platelet count, and hypoalbuminemia are predictors of survival in patients with primary solid cancers . Various combinations of these factors have been used to derive inflammation-based scores to predict survival of subjects with solid cancers [13-16], such as Glasgow Prognostic Score (GPS), Neutrophil/Lymphocyte Ratio (NLR) and Platelet/Lymphocyte Ratio (PLR).
The aims of this study were: (i) to evaluate the prognostic role of inflammatory markers on survival in HIV-infected patients with NHL and, (ii) to compute a prognostic score based on demographical and clinical variables and inflammatory biomarkers.
The Italian MASTER cohort is a hospital-based multicenter, open, dynamic cohort established in the mid-1990s with retrospective patients’ enrolment since 1986. At present, the cohort includes about 24 500 HIV-infected subjects, aged 18 years or older, at eight Infectious Diseases Units throughout Italy. In this work we included HIV-infected patients with a first NHL diagnosis between January 1998 and December 2012. The study design and cancer collection data have been described previously .
We retrieved gender, age, country of origin, HIV exposure risk, date of enrolment in the cohort, Hepatitis B Virus (HBV) or/and Hepatitis C Virus (HCV) co-infection and cART at cancer diagnosis from the MASTER electronic database. The following parameters, measured within 6 months from the diagnosis of NHL, were also retrieved: AIDS event occurrence, HIV-RNA, CD4 cell count, CD8 cell count, CRP, albumin, white blood cell, neutrophil, lymphocyte and platelet counts. The inflammatory based prognostic variables GPS, modified Glasgow Prognostic Score (mGPS), NLR, PLR, Prognostic Index (PI), and Prognostic Nutritional Index (PNI) were defined as shown in Table 1.
The study was conducted in accordance with the guidelines of the Declaration of Helsinki and the principles of Good Clinical Practice. Informed consent was obtained according to the standards of the local ethics committees.
The primary outcome of this study was all-cause mortality. The cumulative risk of death was determined from the data of NHL diagnosis to the end of the observation period. The censor date for survival analysis was 31st December 2012. The observation period ended either on 31st December 2012, or at last follow-up visit or death, whichever occurred first. Patients lost to follow-up contributed to the time at risk until last visit.
At first, the associations of each prognostic variable with all-cause mortality were tested by univariate and multivariate analysis using Cox proportional hazard models. Gender, age at cancer diagnosis, intravenous drug use, AIDS defining event, HBV and/or HCV co-infection, cART, CD4 cell count, HIV-RNA undetectable were included as covariates. The results were expressed as Hazards Ratios (HRs), their 95% Confidence Intervals (95% CIs), and p-values according to Wald test. Due a few missing data in GPS, mGPS, PLR, PI, and PNI, the Cox proportional hazard regression models were adjusted by stabilized inverse probability weights with missing data in prognostic variables considered missing at random . In order to compare the fitness of the model with and without the inflammatory variables, we used the log-likelihood ratio test.
The Area Under the Receiver Operating Characteristic (ROC) Curve (AUC) was used to assess the predictive accuracy of each prognostic variable. As sensitivity analyses, we assessed the consistency of the prognostic role of inflammatory variables in subjects who underwent cART at NHL diagnosis.
Non-linear relationships between continuous prognostic variables and risk of death were assessed by Cox regression models with spline terms of the prognostic variables . We used the Akaike’s information criterion  to assess fitting of models with linear and non-linear terms and to choose the number of spline knots.
Subsequently, we derived a risk score of mortality on the basis of the inflammatory, demographic and clinical variables. To this end, the patients were randomly divided into two subgroups by split-sample technique: (i) derivation subgroup including 70% of subjects, and (ii) validation subgroup, including 30% of subjects. In the derivation sample, a Cox proportional regression model was used to evaluate the association between baseline factors and risk of death, using a cut-off of minimal 10 events per variable for first selection of variables. Gender, age at cancer diagnosis, intravenous drug use, AIDS defining event, HBV and/or HCV co-infection, cART, CD4 cell count, HIV-RNA undetectable and PNI were included. Bootstrapping with 1000 replications was performed to check the stability of variables included in the final models. Linear prediction equations, representing the risk scores, were derived from the final model. For visualization purposes, three risk groups were created according to tertiles of the risk scores distribution in the derivation sample. Cumulative incidence rates of mortality among these three groups were compared by the Kaplan–Meier methods and tested with the log-rank test for trend in the derivation sample. In the derivation and the validation sample the AUC was used to assess the predictive accuracy of the risk scores.
The selection of variables for fitting the most parsimonious models was performed using a backward stepwise procedure, with p = 0.20 for retaining each variable in the model. A graphical check on each regressor did not detect major departures from the proportional hazard assumption of the model.
All the statistical tests were two-sided, assumed a level of significance of 0.05 and were performed using Stata version 12.0 (StataCorp, College Station, TX, USA).
A total of 215 non-Hodgkin lymphoma cases were included. They had a mean age of 43.2 years at time of NHL diagnosis (Standard Deviation, SD = 9.0) and most of them were males (80.0%). The characteristics of the patient population at NHL diagnosis are shown in Table 2. Most of them were assuming cART and plasmatic HIV-RNA was positive (>37 cp/ml). More than half patients had severe immunodeficiency with CD4 cell count less than <200 cells/mm3.
After a median follow up of 2 years (ranging from 2 days to 13.9 years), 18 patients (8.4%) were lost and 98 (45.6%) died. The median survival was 9.3 years. At 1, 3 and 5 years, cumulative risk of death was 34.1%, 44.3% and 45.7%, respectively.
Prognostic role of inflammation-based prognostic variables
Table 3 describes the distribution of inflammatory markers and the association between each of them and risk of death using Cox regression models. Median time from NHL diagnosis to blood testing was 18 days (interquartile range: 3–53 days).
Score 0 was attributed to the majority of patients for NLR, while for PI and PNI most subjects were categorized as score 1. About half of patients were scored with 1 for GPS and mGPS, and with 0 for PLR. No associations were found between prognostic variables and HIV-RNA copy number in subjects with positive HIV-RNA at NHL diagnosis. GPS, mGPS, PI and PNI were independently associated with risk of death using both univariate and multivariate Cox models. These results were also confirmed when we limited the analyses to patients under cART (data not shown).
The likelihood ratio tests for the global fit of the model, before and after addition of inflammatory variables to the full model, were significant for GPS, mGPS, PI and PNI (p < 0.001).
Non-linear relationships between NLR, PLR, PNI and risk of death were evaluated by using multivariate Cox regression models with a restricted cubic-spline for each prognostic score (Figure 1). A trend of increasing risk of death with increasing PNI was observed, although not statistical significant, whereas no trend with PLR and NLR was found.
The ROC analysis showed that mGPS, GPS, PNI and PI had statistically significant discrimination ability, with AUC ranging from 0.60 to 0.69, with follow-up truncated at both 3 and 5 years (Table 4).
Mortality risk score
Finally, we computed a risk score for death on the basis of the final Cox model. In the derivation sample (93 subjects with PNI available, 34 of which died), age, occurrence of an AIDS defining event within six months from date of NHL diagnosis, no cART assumption, CD4 cell count lower than 200 cells/mm3 and PNI equal to 1 were significantly associated with fatal outcome in univariate Cox regression analysis (Table 5). However, in multivariate analyses, only occurrence of AIDS defining event within 6 months from NHL diagnosis and PNI were independently associated with mortality, also using bootstrapping technique.
The linear prediction equation for risk of mortality (Rm) derived from the final Cox model was: Rm = 0. 703 *(occurrence of AIDS defining event within six months from NHL diagnosis [yes =1/no = 0]) + 0.794* (PNI). The predictive accuracies of this risk score for mortality, as measured by the AUC, were 0.69 (95% CI 0.58 to 0.81) and 0.69 (0.57 to 0.81) in derivation sample, and 0.69 (0.49 to 0.90) and 0.73 (0.53 to 0.92) in validation sample (47 subjects with PNI available, 12 of which died), at 3 and 5 years of follow-up, respectively. As an application of the score, patients were categorized into three risk groups: patients without AIDS defining event within six months from the date of NHL diagnosis and with PNI > 45, patients with AIDS defining event within six months from the date of NHL diagnosis or with PNI < 45 and patients with AIDS defining event within six months from the date of NHL diagnosis and with PNI < 45. According to the 3 risk score categories, cumulative mortality rates were 20.8% (10.4 to 38.9), 31.8% (18.8 to 50.3) and 64.5% (43.3 to 84.8) in the low-, intermediate- and high-risk group (log rank test p < 0.001) (Figure 2).
As a sensitive analysis, we have also computed a prognostic score including PNI and AIDS events on the whole cohort, finding AUCs for this score of 67.0 (56.3-77.6) and 67.9 (57.8-78.1) at 3 and 5 years, respectively.
In this study, we investigated the relationship between some biomarkers of systemic inflammation and overall survival in HIV-infected patients with NHL for the first time. We found that some inflammatory-based scores including serum levels of CRP and albumin (GPS, mGPS, PI and PNI), but not those based on cellular components of the white cell count (NLR and PLR), had an independent prognostic value for risk of death. We also computed a mortality risk score including PNI and occurring of an AIDS event within six months from NHL diagnosis, which showed a rather fair predictive ability, with an AUC of 0.69 at 3 years of follow-up.
Various prognostic parameters have been used for NHL, so far. Among them, the IPI is probably the most commonly used and it has been revised recently [10,23,24]. New inflammatory indexes, based on routinely examined blood parameters, have been proposed, including NLR and GPS: NLR has been found useful in patients with diffuse large B-cell lymphoma at low risk of death , whereas an elevated GPS was associated with poorer survival in patients with extranodal natural killer/T-cell lymphoma  and in those with diffuse large cell lymphoma treated with R-CHOP .
In HIV-infected patients with NHL, however, also HIV-infection related factors, such as cART, concurrent infections and baseline immunosuppression, have been found to influence prognosis [4,7,25]. However, the studies carried out so far produced inconsistent results. A low CD4 cell count at NHL diagnosis found to be an independent risk factor for death in pre-cART and early cART eras in two large cohorts [3-5,26], whereas history of AIDS was associated with worse survival only in the pre-cART era . A protective effect of cART has been described [27,28]. Other studies, however, found that only patient- and lymphoma-, but not HIV-related factors were predictive of poor outcome in the cART era [7-9]. For this reason, a composite score for HIV-infected patients with NHL has been recently developed which involves both HIV infection factors and the variables included in IPI . Nevertheless, the IPI is not of common use for HIV patients, because of lack of some necessary data, and other inflammatory scores based on routine data, such as NLR and GPS, may be considered.
In our study, we have included all HIV-infected patients with NHL independently of histology, Ann Arbor stage or chemotherapy regimen. The 2-year cumulative risk of death from NHL diagnosis in our cohort was 45% , slightly lower than the 50-70% value reported in previous studies on HIV-infected subjects [4,5,26,30]. This is possibly due to the fact that previous studies included patients with NHL diagnosis between the late 1990s and early 2000s, when both cART and chemotherapy regimens were not so effective as they are todays.
In these subjects, we found that some indexes including serum CRP and albumin (GPS and mGPS, PNI and PI), as representative of systemic inflammation and nutrition status, were independent predictors of survival in HIV patients with NHL. Therefore, we developed a simple score which integrates clinical and laboratory variables for assessing risk of death in HIV-infected patients with NHL. The score includes one of the HIV-related variables, i.e. an AIDS defining event occurring within six months from NHL diagnosis, and one inflammatory score, the PNI, based on common blood parameters such as lymphocytes and albumin, which therefore can be easily applied to any patient with newly diagnosed NHL. We found that this score had a rather fair predictive ability with an AUC of 0.69 at 3 years of follow-up.
There are several strengths in our study, including a large sample size, multicenter national cohort and inclusion of patients receiving state of the art lymphoma and HIV care. This study has some limits, too. First, we could not analyze the prognostic factors according to type of NHL, due to lack of power of the study for subgroup analyses. Second, the derivation and validation samples used to compute and test the prognostic score were unbalanced (70% and 30%, respectively). In fact, because of the relatively small number of NHL cases, we choose to include more subjects (70%) in the derivation cohort to get better precision of the estimates. Anyway, we also computed a prognostic score including PNI and AIDS events within six months from NHL diagnosis on the whole cohort, finding an AUC for this score similar to those computed using validation and derivation cohorts.
In summary, our study extends previous findings which were focused on solid cancers and shows that some markers of systemic inflammatory response are associated with poor outcome also in HIV-infected patients with NHL.
Combined antiretroviral therapy
International prognostic index
Glasgow prognostic score
Hepatitis B virus; HCV: Hepatitis C virus
Viral hepatitis C
Modified glasgow prognostic score
Prognostic nutritional index
- 95% CI:
95% confidence interval
Receiver operating characteristic
Area under the receiver operating characteristic curve
Risk of mortality
Gotti D, Raffetti E, Albini L, Sighinolfi L, Maggiolo F, Di Filippo E, et al. Survival in HIV-infected patients after a cancer diagnosis in the cART Era: results of an Italian multicenter study. PLoS One. 2014;23:9(4).
Matthews GV, Bower M, Mandalia S, Powles T, Nelson MR, Gazzard BG. Changes in acquired immunodeficiency syndrome-related lymphoma since the introduction of highly active antiretroviral therapy. Blood. 2000;96(8):2730–4.
Robotin MC, Law MG, Milliken S, Goldstein D, Garsia RJ, Dolan GM, et al. Clinical features and predictors of survival of AIDS-related non-Hodgkin’s lymphoma in a population-based case series in Sydney, Australia. HIV Med. 2004;5(5):377–84.
Bower M, Gazzard B, Mandalia S, Newsom-Davis T, Thirlwell C, Dhillon T, et al. A Prognostic Index for Systemic AIDS-Related Non-Hodgkin Lymphoma Treated in the Era of Highly Active Antiretroviral Therapy. Ann Intern Med. 2005;143(4):265–73.
Bohlius J, Schmidlin K, Costagliola D, Fatkenheuer G, May M, Caro Murillo AM, et al. Prognosis of HIV-associated non-Hodgkin lymphoma in patients starting combination antiretroviral therapy. AIDS. 2009;23(15):2029–37.
Hoffmann C, Hentrich M, Gillor D, Behrens G, Jensen B, Stoehr A, et al. Hodgkin lymphoma is as common as non-Hodgkin lymphoma in HIV-positive patients with sustained viral suppression and limited immune deficiency: a prospective cohort study. HIV Med 2014 Sep 23.
Lim ST, Karim R, Tulpule A, Nathwani BN, Levine AM. Prognostic factors in HIV-related diffuse large-cell lymphoma: before versus after highly active antiretroviral therapy. J Clin Oncol. 2005;23(33):8477–82.
Miralles P, Berenguer J, Ribera JM, Rubio R, Mahillo B, Tellez MJ, et al. Prognosis of AIDS-Related Systemic Non-Hodgkin Lymphoma Treated With Chemotherapy and Highly Active Antiretroviral Therapy Depends Exclusively on Tumor-Related Factors. JAIDS. 2007;44(2):167–73.
Lim ST, Karim R, Nathwani BN, Tulpule A, Espina B, Levine AM. AIDS-related Burkitt’s lymphoma versus diffuse large-cell lymphoma in the pre-highly active antiretroviral therapy (HAART) and HAART eras: significant differences in survival with standard chemotherapy. J Clin Oncol. 2005;23(19):4430–8.
A predictive model for aggressive non-Hodgkin’s lymphoma. The International Non-Hodgkin’s Lymphoma Prognostic Factors Project. N Engl J Med. 1993; 329(14):987–94.
Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140(6):883–99.
Roxburgh CS, McMillan DC. Role of systemic inflammatory response in predicting survival in patients with primary operable cancer. Future Oncol. 2010;6(1):149–63.
McMillan DC. The systemic inflammation-based Glasgow Prognostic Score: a decade of experience in patients with cancer. Cancer Treat Rev. 2013;39(5):534–40.
Guthrie GJ, Charles KA, Roxburgh CS, Horgan PG, McMillan DC, Clarke SJ. The systemic inflammation-based neutrophil-lymphocyte ratio: experience in patients with cancer. Crit Rev Oncol Hematol. 2013;88(1):218–30.
Templeton AJ, McNamara MG, Seruga B, Vera-Badillo FE, Aneja P, Ocaña A, et al. Prognostic role of neutrophil-to-lymphocyte ratio in solid tumors: a systematic review and meta-analysis. J Natl Cancer Inst. 2014;106(6):dju124.
Sun K, Chen S, Xu J, Li G, He Y. The prognostic significance of the prognostic nutritional index in cancer: a systematic review and meta-analysis. J Cancer Res Clin Oncol. 2014;140(9):1537–49.
Li X, Zhang Y, Zhao W, Liu Z, Shen Y, Li J, et al. The Glasgow Prognostic Score as a significant predictor of diffuse large B cell lymphoma treated with R-CHOP in China. Ann Hematol. 2015;94(1):57–63.
Troppan K, Deutsch A, Gerger A, Stojakovic T, Beham-Schmid C, Wenzl K, et al. The derived neutrophil to lymphocyte ratio is an independent prognostic factor in patients with diffuse large B-cell lymphoma. Br J Cancer. 2014;110(2):369–74.
Li YJ, Jiang WQ, Huang JJ, Xia ZJ, Huang HQ, Li ZM. The Glasgow Prognostic Score (GPS) as a novel and significant predictor of extranodal natural killer/T-cell lymphoma, nasal type. Am J Hematol. 2013;88(5):394–9.
Seaman SR, White IR. Review of inverse probability weighting for dealing with missing data. Stat Methods Med Res. 2013;22(3):278–95.
Harrell F. General Aspects of Fitting Regression Models. In: Regression Modeling Strategies - With Applications to Linear Models, Logistic Regression and Survival Analysis. New York: Springer-Verlag; 2011. p. 23–4.
Akaike H. A New Look at the Statistical Model Identification. IEEE Trans Autom Control. 1974;AC-19(6):716–23.
Sehn LH, Berry B, Chhanabhai M, Fitzgerald C, Gill K, Hoskins P, et al. The revised International Prognostic Index (R-IPI) is a better predictor of outcome than the standard IPI for patients with diffuse large B-cell lymphoma treated with R-CHOP. Blood. 2007;109(5):1857–61.
Zhou Z, Sehn LH, Rademaker AW, Gordon LI, Lacasce AS, Crosby-Thompson A, et al. An enhanced International Prognostic Index (NCCN-IPI) for patients with diffuse large B-cell lymphoma treated in the rituximab era. Blood. 2014;123(6):837–42.
Simcock M, Blasko M, Karrer U, Bertisch B, Pless M, Blumer L, et al. Treatment and prognosis of AIDS-related lymphoma in the era of highly active antiretroviral therapy: findings from the Swiss HIV Cohort Study. Antivir Ther. 2007;12(6):931–9.
Chao C, Xu L, Abrams D, Leyden W, Horberg M, Towner W, et al. Survival of non-Hodgkin lymphoma patients with and without HIV infection in the era of combined antiretroviral therapy. AIDS. 2010;24(11):1765–70.
Barta SK, Xue X, Wang D, Tamari R, Lee JY, Mounier N, et al. Treatment factors affecting outcomes in HIV-associated non-Hodgkin lymphomas: a pooled analysis of 1546 patients. Blood. 2013;122(19):3251–62.
Gopal S, Patel MR, Yanik EL, Cole SR, Achenbach CJ, Napravnik S, et al. Temporal trends in presentation and survival for HIV-associated lymphoma in the antiretroviral therapy era. J Natl Cancer Inst. 2013;105(16):1221–9.
Barta SK, Xue X, Wang D, Lee JY, Kaplan LD, Ribera JM, et al. A new prognostic score for AIDS-related lymphomas in the Rituximab-era. Haematologica. 2014;123(6):837–42.
Chow KU, Mitrou PS, Geduldig K, Helm EB, Hoelzer D, Brodt HR. Changing incidence and survival in patients with aids-related non-Hodgkin’s lymphomas in the era of highly active antiretroviral therapy (HAART). Leuk Lymphoma. 2001;41:105–16.
We thank Prof. Giampiero Carosi for his thoughtful review of the manuscript, the editor and the reviewers for their constructive comments, which helped us to improve the manuscript, and all patients, doctors, nurses and the dedicated staff of the local centers: Institute of Clinical Infectious Diseases of Spedali Civili Hospital, University of Brescia, Brescia; the Department of Infectious Diseases of San Gerardo de’ Tintori Hospital, Monza; the Department of Infectious Diseases, Papa Giovanni XXIII Hospital, Bergamo; the Department of Infectious Diseases of Istituti Ospedalieri of Cremona, Cremona; Department of Infectious Diseases of Nuovo Polo Ospedaliero S.Anna di Cona, Ferrara; the Department of Infectious Diseases of Santa Maria Annunziata Hospital, Florence; the Institute of Clinical Infectious Diseases of Polyclinic A. Gemelli, University of Sacred Heart, Roma; the Institute of Clinical Infectious Diseases of Polyclinic of Bari, University of Bari, Bari.
The MASTER Cohort and this study were funded by M.I.S.I. foundation (Fondazione Malattie Infettive e Salute Internazionale, http://www.fondazionemisi.it/). The authors declare that they have no competing interests.
The authors participated meaningfully in the study. Specific author contributions were: Study concept and design: RE, QRE and DF; acquisition of data: CF, LN, PG, DFE, SD, CG, BA, ZF, DS, DPM, CF and QRE; analysis and interpretation of data: all author; drafting of the first version of the manuscript: RE, QRE and DF; critical revision of the manuscript for important intellectual content and approval of final version: all authors; statistical analysis: RE; study supervision: DF. All authors read and approved the final manuscript.