METHOD FOR ASSISTING WITH PROGNOSIS

Abstract

The present invention relates to a method for assisting with lymphoma prognosis. The prognosis of therapeutic response of patients with lymphoma is difficult. Based on a study of advanced stage DLBCL patients, the inventors showed that medical imaging such as 18F-FDG-PET/CT can provide a prognostic radiomic signature combining metrics reflecting tumor dissemination and tumor burden. In another aspect, the invention relates to a computer software comprising instructions to implement at least a part of a method according to the invention. In yet another aspect, the invention relates to a computer-readable non-transient recording medium on which a software is registered to implement a method according to the invention.

Claims

1– Method for assisting with lymphoma prognosis comprising: providing imaging data of a patient presenting a lymphoma, segmenting the imaging data so as to identify lesions, determining a distance between the two lesions which are furthest apart, determining the metabolic tumor volume, and providing a prognosis indicator based on the distance between the two lesions which are furthest apart and on the metabolic tumor volume.

2. Method according to claim 1, wherein the provided imaging data are PET scan data.

3. Method according to any of claim 1 or 2, wherein the segmentation is based on a concentration of a tracing agent.

4. Method according to claim 3, wherein the segmentation further comprises erasing regions with a naturally high concentration of tracing agent.

5. Method according to any of claims 1 to 4, further comprising a step consisting in providing help with the decisions between several medical treatments based on the prognosis indicator.

6. Method according to claim 5, wherein the several medical treatments comprise an RCHOP-type treatment and an RACVBP-type treatment.

7. Computer software comprising instructions to implement at least a part of a method according to one of claims 1 to 6 when the software is executed by a processor.

8. Computer-readable non-transient recording medium on which a software is registered to implement a method according to one of claims 1 to 6 when the software is executed by a processor.

Description

[0030] The invention can be better understood at the reading of the detailed examples below, which constitute non-limitative embodiments of the present invention and at the examining of the annexed drawing, on which:

[0031] - FIG. 1 is diagram showing the Kaplan-Meier estimates of progression-free survival (PFS) and overall survival (OS) according to metabolic tumor volume (MTV);

[0032] - FIG. 2 is a diagram showing the Kaplan-Meier estimates of progression-free survival (PFS) and overall survival (OS) according to Dmax.sub.patient, and;

[0033] - FIG. 3 is a diagram showing the Kaplan-Meier estimates of progression-free survival (PFS) and overall survival (OS) according to baseline metabolic tumour volume (MTV) and Dmax.sub.patient.

[0034] It is understood that the described embodiments are not restrictive and that it is possible to make improvements to the invention without leaving the framework thereof.

MATERIALS AND METHODS

Patients

[0035] DLBCL patients with an age adjusted IPI (aaIPI) score of 2 or 3 were randomly assigned to an induction immunochemotherapy with 4 cycles of either R-CHOP14 or R-ACVBP (rituximab, doxorubicin, cyclophosphamide, vindesine, bleomycin, prednisone). Consolidation treatment was driven by centrally reviewed PET assessment according to visual criteria after 2 and 4 treatment cycles. A baseline PET scan was mandatory, with at least one evaluable hypermetabolic lesion. Ethics approval was obtained for this trial, and all patients provided written informed consent to participate.

[0036] For the current analysis, only Ann Arbor stages 3 and 4 patients whose MTV could be computed from a baseline PET/CT scan and with at least two detectable lesions allowing distance measurement were included.

[0037] Baseline patient and disease characteristics, including individual components of the aaIPI score, progression-free survival (PFS) and overall survival (OS) defined according to the revised National Cancer Institute criteria were obtained.

PET/CT Scanning and Quantitative Analysis

[0038] Baseline PET image data in anonymized Digital Imaging and Communications in Medicine (DICOM) format was collected for functional parameter measurements. Quality control rejected scans with burning errors in DICOM retrieval or with a delay of > 90 minutes between 18F-FDG injection and scanning.

[0039] Analysis of PET data was performed by a nuclear medicine physician blinded to patient outcome, using the software according to the invention. MTV was calculated based on a supervised segmentation of tumour regions involving 41% SUVmax thresholding of automatically detected hypermetabolic regions. MTV was defined as the sum of every individual lesion metabolic volume. For each lesion, the tumour lesion glycolysis was calculated as the product of the lesion volume by the SUV mean within the lesion, and the total lesion glycolysis (TLG) was obtained by summing the tumour lesion glycolysis over all lesions. The highest SUVmax of the patient over all lesions and the number of lesions were also reported. Last, several features reflecting the spatial distribution of malignant foci throughout the whole body were computed, based on distance measurements between lesions. Each lesion location was defined as the position of its center and the distances between two lesions were calculated using the Euclidian distance between their centers.

[0040] Four dissemination features were calculated: the distance between the two lesions that were the furthest apart (Dmaxpatient), the distance between the largest lesion and the lesion furthest away from that bulk (Dmaxbulk), the sum of the distances of the bulky lesion from all other lesions (SPREADbulk) and the largest value, over all lesions, of the sum of the distances from a lesion to all the others (SPREADpatient).

Statistical Analysis

[0041] For each PET-derived feature, Receiver Operating Characteristic (ROC) analysis was used to define the optimal cut-off for predicting the occurrence of an event (PFS or OS) by maximizing the Youden index (sensitivity + specificity-1). Sensitivity (Se) and specificity (Sp) were calculated for that cut-off value. Only features with an area under the ROC curve (AUC) greater than 0.6 on PFS were retained for subsequent analyses. Survival functions were calculated by using Kaplan-Meier analyses and the survival distributions were compared using the log-rank test. Multivariate analyses involving MTV and dissemination features were performed using Cox proportional hazard models. Based on these results, a prognostic model combining MTV and a dissemination feature was built on which Kaplan Meier survival analysis was performed. Correlations between dissemination features and MTV were assessed using chi-squared tests. Mann-Whitney tests were used to test whether the patient size and MTV were significantly different in patients with low and high dissemination features.

[0042] Statistical significance was set to p<0.05. All statistical analyses were performed using MedCalc software (MedCalc Software, Ostend, Belgium).

RESULTS

[0043] In total, 95 patients, were included, whom clinical characteristics are summarised in Table 1.

TABLE-US-00001 Patient characteristics Patient characteristics n=95 (%) Sex Female 42 (44%) Male 53 (56%) Age (median, ranges) years 46 (18-59) Height (median, ranges) cm 173 (140-193) ≤ 170 cm 42 >170 cm 53 Ann Arbor stage III 9 (9.5%) IV 86 (90.5%) Performance status 1 27 (28.4%) 2 44 (46.3%) 3 19 (20%) 4 5 (5.3%) aaIPI 1 3 (3%) 2 69 (73%) 3 23 (24%) Treatment R-ACVBP 14 46 (48%) R-CHOP 14 49 (52%)

[0044] With a median follow-up of 44 months (range 27-63 months), the 4-year PFS and OS rates for the whole group were 77 % and 85 % respectively. Twenty-two patients had a PFS event with a median of 7 months, 12 in R-CHOP group and 10 in R-ACVBP group. Thirteen patients died with a median of 13 months, 8 in R-CHOP group and 5 in R-ACVBP group. Using log rank tests, neither performance status (0-1 versus 2-3) nor aaIPI (2 versus 3) were significantly associated with PFS (p=0.17, p=0.21) or OS (p=0.41, p=0.46). No significant prognostic impact of chemotherapy regimen (R-CHOP vs R-ACVBP) was observed for both PFS (p=0.69) and OS (p=0.48).

PET Features

[0045] Table 2 shows the descriptive statistics for the PET features and Table 3 gives the results of ROC analyses performed on each PET parameter.

TABLE-US-00002 Median, range, mean and standard of deviation (SD) of PET features PET Parameters PFS OS AUC Cut-off Se Sp AUC Cut-off Se Sp MTV (cm.sup.3) 0.64 394 68 60 0.69 468 77 71 SUVmax 0.58 15 41 85 0.53 23 46 71 TLG 0.53 4396 45 68 0.67 4550 61 73 Dmax.sub.patient (cm) 0.65 58 68 74 0.59 58 69 69 Dmax.sub.Bulk (cm) 0.63 43 54 82 0.60 43 54 80 SPREAD.sub.patient (cm) 0.65 1023 50 85 0.58 716 54 71 SPREAD.sub.Bulk (cm) 0.65 530 54 86 0.59 407 61 71 Nb of VOIs / patient 0.64 23 54 77 0.57 20 54 67 PET Parameters median range mean SD MTV 375 27-2525 469 392 SUVmax 20 4-49 21 8 TLG 3275 166-19428 4298 3323 Dmax.sub.patient (cm) 45 7-135 46 25 Dmax.sub.Bulk (cm) 32 7-101 32 17.5 SPREAD.sub.patient (cm) 367 7-11915 798 1420 SPREAD.sub.Bulk (cm) 205 7-4561 425.4 620 Nb of VOIs / patient 13 2-130 20 21

Table 3: ROC analysis of PET features, Area under the ROC curve (AUC), Sensitivity (Se), Specificity (Sp).

[0046] Using ROC optimal cut-off, MTV was highly predictive of outcome (PFS: p=0.027 and OS: p=0.0007) (Table 4). Patients with a high MTV had a significantly worse outcome with a 4-year PFS and OS of 67% and 73% versus 84% and 95% for patients with a lower MTV (FIG. 1).

TABLE-US-00003 PET parameters associated with PFS and OS in Log-rank Cox tests PFS OS HR (CI 95%) 4y-PFS (CI 95%) p HR (CI 95%) 4y-OS (CI 95%) P Low MTV High MTV 1 (ref) 2.6 (1.1-6.0) 84% (79-89) 67% (60-74) 0.027 1 (ref) 6.9 (2.1-21.9) 95% (92-98) 66% (56-76) 0.0007 Low Dmax.sub.patient High Dmax.sub.patient 1 (ref) 4.6 (1.9-11.2) 88% (84-92) 55% (47-63) 0.0003 1 (ref) 4.2 (1.3-13.1) 93% (90-96) 69% (60-78) 0.0095 Low Dmax.sub.bulk High Dmax.sub.bulk 1 (ref) 4.1 (1.5-11.3) 86% (82-90) 52% (42-62) 0.0003 1 (ref) 3.3 (1-11.3) 91% (88-94) 68% (57-79) 0.023 Low SPREAD .sub.patient High SPREAD .sub.patient 1 (ref) 3.7 (1.3-10.1) 85% (81-89) 52% (42-62) 0.0011 1 (ref) 1.9 (0.5-6.8) 86% (81-91) 78% (70-85) 0.24 Low SPREAD .sub.bulk 1 86% (82-90) <0.0001 1 (ref) 90% (87-93) 0.056 High SPREAD .sub.bulk 4.9 (1.7-13.9) 45% (35-55) 2.8 (0.8-9.9) 69% (59-79) Low nb of ROIs 1 (ref) 85% (81-89) 0.0052 1 (ref) 87% (82-92) 0.21 High nb of ROIs 3.1 (1.2-7.9) 58% (49-67) 1.9 (0.6-6.4) 79% (72-86)

[0047] Regarding the dissemination features, ROC AUC were always greater than 0.6 for PFS, and close to 0.6 for OS (Table 3). Table4 shows that Dmaxpatient > 58 cm, Dmaxbulk > 43 cm, SPREADpatient > 1020 cm and SPREADbulk > 530 cm were negative prognostic factors for PFS (p=0.0003, p=0.0003, p=0.0011, p<0.0001 respectively) and that for OS, only Dmaxpatient and Dmaxbulk were statistically significant (p=0.0095 and p=0.023 respectively).

[0048] No significant differences in height were observed between patients with low and high Dmaxpatient (p=0.96). Similarly, no significant differences in MTV were observed between patients with low and high Dmaxpatient (median of 344 cm3 and 415 cm3 respectively, p=0.14).

Combination MTV and Dissemination Features

[0049] In multivariate Cox regression analysis including MTV and Dmaxpatient, Dmaxpatient was significantly associated with PFS (p=0.0014; HR=4.3) whereas MTV was borderline significant (p=0.056; HR=2.3). For OS, both factors remained significant (p=0.037 HR=4.0 for MTV and p=0.029; HR=3.7 for Dmaxpatient).

[0050] Three risk categories could therefore be significantly distinguished on the basis of the presence or absence of high MTV (> 394 cm3) or Dmaxpatient (> 58 cm) (p=0.0003 for PFS and p=0.0011 for OS): group 1 with no risk factor (n=36), group 2 with one risk factor only (n=41), group 3 with both (n=18), with 4-year PFS rates of 94%, 73%, and 50%, respectively and 4-year OS rates of 97%, 88%, and 53%, respectively. Group 2 vs group 3 had significantly different PFS (p=0.041) and OS (p=0.019). Group 1 vs group 2 had significantly different PFS (p=0.013) whereas OS did not reach significance (p=0.13).

DISCUSSION

[0051] Lymphoma is a group of blood cancers that develop from lymphocytes. Although most cells in the body can migrate at one or more distinct steps during their development and differentiation, the trafficking propensity of lymphocytes is unrivaled among somatic cells. In case of malignant transformation, this property allows for rapid tumor dissemination irrespective of the conventional anatomic boundaries limiting early spread in most types of cancer. Thus, the disease can spread rapidly to different parts of the body, involving lymph nodes, possibly associated with extra nodal sites.

[0052] 18F-FDG-PET/CT is the current state-of-the-art imaging scan in lymphoma. Recent advances in PET imaging revealed that MTV, as a surrogate for tumour cell number, has a strong prognostic value in DLBCL, much higher than the presence of a bulk. Recently, this was confirmed in a large phase 3 study, GOYA, including more than 1100 patients (NCT0128774): MTV quartiles stratified the population in quartiles 1, 2, 3 and 4 with a three-year PFS of 86%, 84%, 78% and 66% respectively. The present invention is based on the demonstration that MTV maintained its prognostic power in a cohort of advanced stage patients. Patients stage 3 or 4 were significantly stratified in two different risk categories according to their MTV. Moreover, using ROC analysis, MTV was the only significant feature on both PFS and OS. It was superior to standard features such as aaIPI for both PFS and OS. A high MTV identify 64% of the PFS events (14/22).

[0053] In this study, new radiomic features were introduced, extracted from PET scans to quantify tumour dissemination. Several of these features based on distance measurement between lymphoma lesions were significant for PFS and OS in a group of stage 3 and stage 4 patients, suggesting that an advanced characterization of the lesion dissemination is relevant even among patients with an advanced disease. In particular, the distance between the two lesions that were the furthest apart, Dmaxpatient, had strong predictive power for PFS and OS. A high Dmaxpatient was associated with an adverse outcome, with a 4-years PFS and OS of 55% and 69% respectively. Similarly, SPREADpatient and SPREADbulk combining spatial dispersion information and the number of lesion, were very significantly associated with PFS (Table 4).

[0054] Dmaxpatient is a very simple 3D feature to calculate with an intuitive interpretation. The height did not influence Dmaxpatient, as height did not significantly differ between high or low Dmaxpatient groups.

[0055] In multivariate analysis, MTV and Dmaxpatient remained independently significant for PFS and OS. The model combining these two factors was able to separate significantly three different prognostic groups: group 1 with no risk factor, group 2 with one risk factor, group 3 with both. Specifically, this could identify a group with a poor prognosis so that clinicians might consider changing treatment. Indeed, patients with high baseline MTV (>394 cm3) and high Dmaxpatient (>58 cm) had a much worse prognosis than the other patients with 4-year PFS 50% and 4-year OS 53%. This group represented 19% of the cohort and included 41% of the PFS total number of events (9/22) and 54% of the OS total number of events, making this model useful for identifying patients with poor prognosis.

[0056] In the LNH073B trial, consolidation treatment was driven by centrally reviewed PET assessment after 2 and 4 cycles: patients who were classified as PET 2 and PET 4 negative received standard immunochemotherapy consolidation; patients classified as PET 2 positive and PET4 negative received 2 cycles of high-dose methotrexate (3 g/m2) and then a high-dose therapy (carmustine, etoposide, cytarabine, melphal [BEAM] or zevalin, carmustine, etoposide, cytarabine, melphalan[Z-BEAM]), followed by Autologous Stem Cell Transplantation (ASCT); PET 4 positive patients had salvage regimen followed by ASCT in responders to salvage. Despite this 18F-FDG-PET-driven consolidation strategy that might actually decrease the prognostic impact of baseline PET features, MTV and dissemination features remained significantly predictive of PFS and OS. Further studies are needed to more comprehensively establish the role dissemination features might play in lymphomas when measured at baseline and during patient monitoring. Imaging of CXCR4 receptors could be helpful in this regard.

CONCLUSION

[0057] 18F-FDG-PET/CT can provide a prognostic radiomic signature combining metrics reflecting tumor dissemination and tumor burden. In this study of advanced stage DLBCL patients, combining MTV and Dmaxpatient improved patient risk stratification at staging.

[0058] Unless otherwise specified, the word “or” is equivalent to “and/or”. Similarly, the word ‘one’ is equivalent to ‘at least one’ unless the contrary is specified. Unless otherwise specified, all percentages are weight percentages.