Lymphocyte Doubling Time As A Key Prognostic Factor To Predict Time To First Treatment In Early-Stage Chronic Lymphocytic Leukemia

<jats:p>The prognostic role of lymphocyte doubling time (LDT) in chronic lymphocytic leukemia (CLL) was recognized more than three decades ago when the neoplastic clone’s biology was almost unknown. LDT was defined as the time needed for the peripheral blood lymphocyte count to double the of the initial observed value. Herein, the LDT prognostic value for time to first treatment (TTFT) was explored in our prospective O-CLL cohort and validated in in two additional CLL cohorts. Specifically, newly diagnosed Binet stage A CLL patients from 40 Italian Institutions, representative of the whole country, were prospectively enrolled into the O-CLL1-GISL protocol (<jats:ext-link>clinicaltrial.gov</jats:ext-link> identifier: NCT00917540). Two independent cohorts of newly diagnosed CLL patients recruited respectively at the Division of Hematology in Novara, Italy, and at the Hospital Clinic in Barcelona, Spain, were utilized as validation cohorts. In the training cohort, TTFT of patients with LDT &amp;gt;12 months was significantly longer related to those with a shorter LDT. At Cox multivariate regression model, LDT ≤ 12 months maintained a significant independent relationship with shorter TTFT along with <jats:italic>IGHV</jats:italic> unmutated (<jats:italic>IGHV</jats:italic>unmut) status, 11q and 17p deletions, elevated β2M, Rai stage I-II, and <jats:italic>NOTCH1</jats:italic> mutations. Based on these statistics, two regression models were constructed including the same prognostic factors with or without the LDT. The model with the LTD provided a significantly better data fitting (χ<jats:sup>2</jats:sup> = 8.25, P=0.0041). The risk prediction developed including LDT had better prognostic accuracy than those without LDT. Moreover, the Harrell’C index for the scores including LDT were higher than those without LDT, although the accepted 0.70 threshold exceeded in both cases. These findings were also confirmed when the same analysis was carried out according to TTFT’s explained variation. When data were further analyzed based on the combination between LDT and <jats:italic>IGHV</jats:italic> mutational status in the training and validation cohorts, <jats:italic>IGHV</jats:italic>unmut and LDT&amp;gt;12months group showed a predominant prognostic role over <jats:italic>IGHV</jats:italic>mut LTD ≤ 12 months (P=0.006) in the O-CLL validation cohort. However, this predominance was of borden-line significance (P=0.06) in the Barcelona group, while the significant prognostic impact was definitely lost in the Novara group. Overall, in this study, we demonstrated that LDT could be re-utilized together with the more sophisticated prognostic factors to manage the follow-up plans for Binet stage A CLL patients.</jats:p>

section_content: The prognostic role of lymphocyte doubling time (LDT) in chronic lymphocytic leukemia (CLL) was recognized more than three decades ago when the neoplastic clone's biology was almost unknown. LDT was defined as the time needed for the peripheral blood lymphocyte count to double the of the initial observed value. Herein, the LDT prognostic value for time to first treatment (TTFT) was explored in our prospective O-CLL cohort and validated in in two additional CLL cohorts. Specifically, newly diagnosed Binet stage A CLL patients from 40 Italian Institutions, representative of the whole country, were prospectively enrolled into the O-CLL1-GISL protocol (clinicaltrial. gov identifier: NCT00917540). Two independent cohorts of newly diagnosed CLL patients recruited respectively at the Division of Hematology in Novara, Italy, and at the Hospital Clinic in Barcelona, Spain, were utilized as validation cohorts. In the training cohort, TTFT of patients with LDT >12 months was significantly longer related to those with a shorter LDT. At Cox multivariate regression model, LDT ≤ 12 months maintained a significant independent relationship with shorter TTFT along with IGHV unmutated (IGHVunmut) status, 11q and 17p deletions, elevated b2M, Rai stage I-II, and NOTCH1 mutations. Based on these statistics, two regression models were constructed including the same

section_content: The heterogeneous course and outcome of chronic lymphocytic leukemia (CLL) are associated with clinical and laboratory parameters as well as the molecular and cytogenetic complexity of the leukemic clone can contribute to setting a prognosis. In the past, attention was centered primarily on the predictors of general outcome, intended as overall survival. In contrast, more recently, the focus also has included the definition of time to first treatment (TTFT). According to the current guidelines (32, 33), CLL is treated at progression. Many patients and clinicians often see the start of therapy as a partition between a healthy and a disease condition. Besides, an accurate prediction of the TTFT permits setting the appropriate follow-up strategy. Several methodologies have been proposed to predict the patient's course and outcome, as well as the TTFT. Although the availability of a battery of cellular/ molecular markers has opened the way to an always more refined prognostic stratification of patients, the current practice's reality indicates difficulties in carrying out several sophisticated features, often confined to a research setting (3). Thus, simplifications introducing inexpensive tests suitable for the clinical setting would be more than welcome, and such prognostic indexes would have a high likelihood of broad applicability. 

The breakthrough of novel biologic variables has led to several prognostic indexes to weigh TTFT in early-stage (Binet A) CLL patients (1, 2, (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52). Indeed, over the past few years, there has been a great effort to use novel molecular markers in prognostic modeling. Yet, questions about their usefulness to improve clinical prediction have been recently debated (53). On the other hand, as the concern for novel molecular markers detected by cutting-edge technologies has soared, performance measures for statistically quantifying their prognostic added value have risen accordingly (52). 

The German study group (53) has recently established the independent prognostic value of lymphocyte doubling time (LDT) in Binet stage A patients after more than 50 years following the recognition of a correlation between the lymphocyte proliferation pattern and clinical outcome in CLL (54). LDT, defined as the period needed for the peripheral blood lymphocyte count to reach a double value of that corresponding to the initial observation, is a simple parameter that is useful in arriving at an accurate prognosis in CLL. Whereas a high LDT (greater than 12 months) identifies a population with an excellent prognosis, a low LDT (less than or equal to 12 months) LDT predicts rapid disease progression in patients in the early clinical stages (53, 55, 56). The raising question is whether the re-introduction of LDT among the last generation prognostic factors could help determine the followup strategies for an early-stage patient, possibly improving or maintaining the prediction power of more hi-tech markers such as the IGHV gene status (51, 57). 

Herein, we investigated the LDT predictive value for TTFT in our prospective O-CLL cohort. The results of LDT prognostication power in the O-CLL training cohort were validated in two additional CLL cohorts.
section_name: INTRODUCTION

section_name: MATERIALS AND METHODS

section_content: LDT was evaluated at diagnosis, also utilizing lymphocyte counts antecedent the enrollment, if available, and defined as the time needed for the peripheral blood lymphocyte count to reach a value double that of the initial observation (52, 55, 56). LDT was calculated as reported by Hoechstetter et al. (52) through a linear regression based on four blood lymphocyte measurements, each at an interval of a minimum of four weeks from the precedent one in no more than six months before enrollment. 

O-CLL1-GISL protocol (clinicaltrial. gov identifier: NCT00917540). The ethics committees from each participating center approved this study. Informed consent was obtained from all subjects. Recruitment began in January 2007 and terminated in January 2012. According to the guidelines (32, 33), treatment was decided uniformly for all participating centers based on documented progressive and symptomatic disease. The present analysis was carried out in 498 out of 523 accrued cases where LDT was available. 

All patients from the O-CLL cohort were studied for CD38, and ZAP-70 expression, IGHV mutational status, FISH assays, and NOTCH1 and SF3B1 gene mutations as previously described (18, 28, 29, 58, 59). 

The contribution of the single institutions of the training cohort is shown in Supplementary Table 1.
section_name: Lymphocyte Doubling Time

section_content: An independent cohort of newly diagnosed and prospectively followed CLL patients recruited since 2001 at the Division of Hematology, Department of Translational Medicine, UPO, Novara, Italy, was utilized as a first validation cohort. The present analysis was restricted to 276, with LDT available, out of 283 cases included in a previous paper (5). All the prognostic factors required in this study (IGHV mutational status, Rai stage, b2M, 17(p) and 11(q) deletions, NOTCH1 coding gene mutation, and LDT were available for 257 cases. 

A further independent cohort of newly diagnosed and prospectively followed CLL patients recruited since 2001 at Hospital Clinic, Institute of Hematology and Oncology, University of Barcelona, Spain, was utilized as an additional validation cohort. The present analysis was performed in 414 cases; 355 were included in a recent paper (51). All the prognostic factors required in this study (IGHV mutational status, RAI stage, b2M, 17(p) and 11(q) deletions, NOTCH1 coding gene mutation, and LDT) were available for 247 cases.
section_name: Validation Cohorts

section_content: TTFT analyses (including the identification of risk factors for this endpoint) were performed using the Kaplan-Meier method followed by log-rank test. The prognostic impact of specific risk factors on the outcome variable was investigated by univariate and multiple Cox regression analysis. Results are expressed as hazard ratios (HR) and 95% confidence intervals (CI). The main prognostic factor in our study was the b Score. It was calculated by deriving a weight corresponding to the regression coefficients of each prognostic factor (b) (60). The regression coefficients of the independent prognostic factors were preliminarily summed up. Then, they were divided by this sum and multiplied by 100, thus deriving a weight ranging from 0 to a given percentage, These weights were summed up on an individual basis, thus deriving a score interpretable in a prognostic scale ranging from 0 to 100% (for patients exposed to all risk factors) (60). The predictive accuracy of the prognostic models was quantified by calculating the Harrell Cindex (HC-index), ranging from 0. 5 to 1. 0, the explained variation on the outcome (i. e., an index combining calibration and discrimination) (61), and the Akaike weights (AIC) (62). The integrated discrimination improvement (IDI) (53) was also calculated to assess the gain in prognostic accuracy provided by LDT. IDI is an index of risk re-classification that quantifies whether a new variable offers a clinically relevant improvement in prediction beyond and above provided by a model based on a previous risk prediction rule and not including the same variable. 

The fittings between two nested prognostic models (including and not including LDT) were compared by the -2 log likelihood statistics (63). Data analysis was performed by STATA for Windows v. 9 and SPSS Statistics v. 21.
section_name: Statistical Analysis

section_name: RESULTS

section_content: In the prospective O-CLL cohort, LDT was available in 498 Binet stage A CLL cases (median 72. 4 months). Seventy-seven cases (15. 5%) presented a LDT ≤ 12 months. Out of the 498 patients, 177 needed treatment, with a significantly higher percentage (P<0. 0001) of cases requiring therapy detected among the group with LDT ≤ 12 months (54/77, 70. 1%) as compared with that with LDT>12 months (123/431, 29. 2%). TTFT of patients with LDT ≤12 months was significantly shorter (HR=2. 9, 95% CI 2. 1-4. 0, P<0. 0001) compared to those with a longer LDT (Figure 1 ). In the same cohort, when all the correlates of outcome determined by univariate analysis (Table 1 ) were introduced into the same multiple Cox regression model, only IGHV unmutated (IGHVumut) genes, 11q and 17p deletions, elevated b2M, Rai stage I-II, NOTCH1 mutations and LDT ≤ 12 months maintained a significant independent relationship with shorter TTFT (Figure 2 ). 

Starting from this analysis, two regression models were constructed, i. e., a Cox model including all significant and independent correlates of TTFT except LDT (Model 1, Table 2 ) and a Cox model including the same set of prognostic factors and LDT (Model 2, Table 2 ). Of note, Model 2 (including LDT) provided a significantly better data fitting (c 2 = 8. 25, P=0. 0041) than Model 1 (not including LDT). Moreover, by IDI calculations, we demonstrated that LDT increased the estimated risk of +3. 8%, a result of high statistical significance (P<0. 001). These analyses were carried out in 334 cases in which all the variables were available.
section_name: O-CLL Training Cohort

section_content: Based on Model 1 and Model 2, two risk prediction rules were developed, i. e. b Score LDT and b Score no LDT (Figure 3 ). b Score which included LDT, had better prognostic accuracy than that without LDT (Table 1, training set). Moreover, the HC-index including LDT, was higher than those without LDT (75. 4 versus 74. 7), although the accepted 0. 70 threshold (25) exceeded in both cases (Table 3 ). This was also true when the same analysis was carried out according to the explained variation in TTFT, which combines the discrimination and the calibration abilities of a risk prediction rule ( b Score LDT = 47. 6% versus b Score no LDT = 45. 0%; Table 3 ).
section_name: Prediction Risk Scores in Training O-CLL Cohort

section_content: These analyses were carried out in two validation cohorts. LDT was available in 276 and in 414 Binet stage A cases respectively in Novara and Barcelona cohorts. Eighty-six and 148 patients required therapy in the Novara and Barcelona cohorts, respectively. A significantly higher percentage of treated cases were recorded among the group with LDT ≤ 12 months in both the Novara (LDT ≤ 12 70/121, 57. 9% versus LDT>12 months (16/155, 10. 3%) and the Barcelona (LDT ≤ 12 months 30/39, 76. 9% versus LDT>12 months cohort 119/375, 31. 7%). In the Novara cohort, Cox univariate analysis showed a significantly increased risk of treatment for patients with LDT ≤12 months (155 cases) (HR=6. 1, 95% CI 3. 6-10. 6, P<0. 0001) compared to those with a longer LDT (121 patients) (Figure 4A ). Similar results were detected in the Barcelona validation set, in which the 39 cases with LDT ≤12 months showed a risk of being treated 7. 4 times higher (95% CI 4. 9-11. 2, P<0. 0001) than the 375 cases with a longer LDT (HR=1) (Figure 4B ). 

Table 4 reports the clinical and biological variables significantly associated with TTFT in the multivariate models of training as well as both validation cohorts. Remarkably, the higher prognostic value, provided by risk prediction rule including LDT, found in the training cohort was fully confirmed in the two validation cohorts (Table 3 ). The HC-indexes and the explained variations in TTFT were consistently higher for bScore including LDT than for that excluding this variable (Table 3 ) in both Novara cohort (HC-Indexes 80. 4 versus 75. 2; explained variations 49. 6 versus 33. 3) and Barcelona cohort (HC-Indexes 70. 0 versus 68. 1; explained variations 30. 5 versus 25. 8). Finally, the Akaike weights coherently indicated that the risk scores including LDT had a chance to provide the best prognostic estimates ranging from 98. 1% to 100% in both the training and in the two validation cohorts (Table 3 ).
section_name: Prediction Risk Scores in Validation Novara and Barcelona Cohorts

section_content: Data were further analyzed based on the combination between LDT and IGHV mutational status in the training and validation cohorts, in cases with both variables available As expected, LDT ≤12 months significantly maintained its negative Interestingly, the probability of remaining therapy-free at 6 years was 51% in the IGHVunmut and LDT>12months group and of 76% in the IGHVmut LTD ≤ 12 months, respectively (P=0. 006) in the O-CLL validation cohort, confirming the predominant prognostic role of the IGHVumut status. However, this predominance was of borden-line significance (P=0. 06) in the Barcelona group, while the significant prognostic impact was definitely lost in the Novara group. This discrepancy could be due to a different distribution of the above-mentioned subsets (IGHVunmut and LDT>12months and IGHVmut LTD ≤ 12 months) among Novara validation cohort (Supplementary Figure 1 ).
section_name: Relationship Between LDT and IGHV Mutational Status

section_content: The clinical course of early-stage CLL is hugely heterogeneous. While some patients need treatment at the onset of the diagnosis, others remain therapy-free for many years or even do not receive any treatment livelong (1, 2). Several prognostic algorithms derived from multivariable models, nomograms and score systems have been developed to predict clinical outcomes accurately in early-stage CLL (64). The IGHV gene configuration is one of the most important single factors predicting therapy need, and it is recurrently incorporated in all prognostic models (65). 

In contrast, the prognostic role of LDT in CLL was acknowledged more than 35 years ago by Montserrat et al. (55) and soon after by Molica et al. (56), when disease biology of the neoplastic cell remained weakly recognized. The raising question is whether the re-introduction of LDT among the last generation Our results fully agree with those of the German CLL Study Group, which endorses the prognostic benefit of determining LDT in the era of new markers (52). Specifically, the German prognostic model includes genetic features, i. e., 17p and 11q deletions, as well as IGHV mutational status. Similarly, the same prognostic indicators were demonstrated by our group to be independently associated with TTFT together with LDT. Unlike CLL1-PM, NOTCH1 remained significant in our analysis, confirming previous results (3, 5, 18). Remarkably, the higher prognostic value, provided by risk prediction score including LDT, found in the training cohort was fully confirmed in the two validation cohorts. Specifically, the indicators of performance, such as the HC-indexes and the explained variations, were consistently higher for bScore including LDT than for that excluding this variable in the training and both validation cohorts. Finally, the Akaike weights coherently indicated that the risk scores including LDT had a chance to provide the best prognostic estimates ranging from 98. 1% to 100% in both the training and the two validation cohorts. More recently, the prognostic significance of LDT was demonstrated to be independent for TTFT from the CLL-IPI and the Barcelona/Brno prognostic models in a real-life cohort of 848 Binet stage A patients (66). 

Thus, the question of whether the re-introduction of LDT among the more sophisticated prognostic factors could help to determine the follow-up strategies for an early-stage patient has an affirmative answer. However, several caveats have to be considered. Although simple, LDT determinations require a precise timing and relative frequent accesses of patients to the clinic, thus preventing the setting of a definite 'watch and wait' strategy in concomitance with the work-up at diagnosis. Alternatively, more precise and possibly more rapid methodologies to measure lymphocyte proliferating potential, such as labeling with deuterated water, are too complex to be used routinely (67). Moreover, when data were further analyzed based on the combination between LDT and IGHV mutational status in the training and validation cohorts, IGHVunmut and LDT>12 months group showed a predominant prognostic role over IGHVmut LTD ≤ 12 months (P=0. 006) in the O-CLL validation cohort. However, this predominance was of borden-line significance (P=0. 06) in the Barcelona group, while the significant prognostic impact was definitely lost in the Novara group, incongruity possibly due to a different distribution of the above-mentioned subsets among Novara validation cohort. Thus, the IGHV mutation status could offer some, likely marginal, sensitivity advantage over the LDT determination. Finally, no information on the LDT stability overtime is available. Therefore, LDT values may vary over time, particularly in concomitance with the progression of the CLL towards a more aggressive form. These changes are likely to relate to a minority of cases and may not influence an entire cohort's data but may have clinical relevance for the individual patients. These variations in time have been reported for cellular markers such as CD38 or ZAP-70 and have caused their subsequent obsolescence. 

The reviewer FA declared a shared affiliation with one of the authors, AN, to the handling editor at time of review. 

Publisher's Note: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
section_name: DISCUSSION

section_name: ACKNOWLEDGMENTS

section_content: In addition to the listed Authors, the following Investigators participated in this study as part of the Gruppo Italiano Studio Linfomi (GISL) : Gianni Quintana, Divisione di Ematologia, Presidio Ospedaliero "A. Perrino", Brindisi ; Giovanni Bertoldero, Dipartimento di Oncologia, Ospedale Civile, Noale, Venezia ; Paolo Di Tonno, Dipartimento di Ematologia, Ospedale di Venere, Bari ; Robin Foà and Francesca R Mauro, Divisione di Ematologia, Università La Sapienza, Roma; Nicola Di Renzo, Unità di Ematologia, Ospedale Vito Fazzi, Lecce; Maria Cristina Cox, Ematologia, A. O. Sant'Andrea, Università La Sapienza, Roma ; Stefano Molica, Dipartimento di Oncologia ed Ematologia, Pugliese-Ciaccio Hospital, Catanzaro ; Attilio Guarini, Unità di Ematologia e Trapianto di Cellule Staminali, Istituto di Oncologia "Giovanni Paolo II", Bari ; Antonio Abbadessa, U. O. C. di Oncoematologia Ospedale "S. Anna e S. Sebastiano", Caserta ; Francesco Iuliano, U. O. C. di Oncologia, Ospedale Giannettasio, Rossano Calabro, Cosenza; Omar Racchi, Ospedale Villa Scassi Sampierdarena, Genova ; Mauro Spriano, Ematologia, A. O. San Martino, Genova; Felicetto Ferrara, Divisione di Ematologia, Ospedale Cardarelli, Napoli ; Monica Crugnola, Ematologia, CTMO, Azienda Ospedaliera Universitaria di Parma ; Alessandro Andriani, Dipartimento di Ematologia, Ospedale Nuovo Regina Margherita, Roma; Nicola Cascavilla, Unità di Ematologia e Trapianto di Cellule Staminali, IRCCS Ospedale Casa Sollievo della Sofferenza, San Giovanni Rotondo ; Lucia Ciuffreda, Unità di Ematologia, Ospedale San Nicola Pellegrino, Trani; Graziella Pinotti, U. O. Oncologia Medica, Ospedale di Circolo Fondazione Macchi, Varese ; Anna
section_name: ACKNOWLEDGMENTS

section_name: FUNDING

section_content: Associazione Italiana Ricerca sul Cancro (AIRC) Grant 5 x 1000 n. 9980 (to FM, MF, and AN); AIRC and Fondazione CaRiCal co-financed Multi-Unit Regional Grant 2014 n. 16695 to FM; AIRC, Special Program Metastases (n. 21198) 5x1000 to RF and GG; RF-2018-12365790, MoH, Rome, Italy (to GG); AIRC IG-5506 to GF, IG-14326 (to MF), IG-15426 (to FF); Compagnia S. Paolo, Turin, Italy, Project 2017. 0526 (to GF) and by the Ministry of Health (Project 5x1000, 2015 and 2016 and Current Research 2016 (to GF, GC, and FF); Swiss Cancer League, ID 3746, 4395 4660, and 4705, Bern, Switzerland; European Research Council (ERC) Consolidator Grant CLLCLONE, ID: 772051 ; Swiss National Science Foundation, ID 320030_169670/1 and 310030_192439, Berne, Switzerland; The Leukemia & Lymphoma Society, Translational Research Program, ID 6594-20, New York; Italian Ministry of Health_Ricerca Corrente 2021 to AN.
section_name: FUNDING

section_name: DATA AVAILABILITY STATEMENT

section_content: The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
section_name: DATA AVAILABILITY STATEMENT

section_name: ETHICS STATEMENT

section_content: The studies involving human participants were reviewed and approved by Comitato Etico Regionale Liguria n°OMC07. 002. The patients/participants provided their written informed consent to participate in this study.
section_name: ETHICS STATEMENT

section_content: FMo, MG, MF, and AN designed the study. FMo, GT, MG, and GD'A performed statistical analysis. FMo, MG, GT, DR, GG, EMo, MF, GC, MC, FF, FS, SB, SF, and AN analyzed and interpreted data, and wrote the manuscript. All authors contributed to the article and approved the submitted version.
section_name: AUTHOR CONTRIBUTIONS

section_content: The Supplementary Material for this article can be found online at: https://www. frontiersin. org/articles/10. 3389/fonc. 2021. 684621/full#supplementary-material Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
section_name: SUPPLEMENTARY MATERIAL