Evolution of TP53 abnormalities during CLL disease course is associated with telomere length changes

<jats:title>Abstract</jats:title><jats:sec>                 <jats:title>Background</jats:title>                 <jats:p>Telomeres are protective structures at chromosome ends which shorten gradually with increasing age. In chronic lymphocytic leukemia (CLL), short telomeres have been associated with unfavorable disease outcome, but the link between clonal evolution and telomere shortening remains unresolved.</jats:p>               </jats:sec><jats:sec>                 <jats:title>Methods</jats:title>                 <jats:p>We investigated relative telomere length (RTL) in a well-characterized cohort of 198 CLL patients by qPCR and focused in detail on a subgroup 26 patients who underwent clonal evolution of <jats:italic>TP53</jats:italic> mutations (evol<jats:italic>TP53</jats:italic>). In the evol<jats:italic>TP53</jats:italic> subgroup we explored factors influencing clonal evolution and corresponding changes in telomere length through measurements of telomerase expression, lymphocyte doubling time, and BCR signaling activity.</jats:p>               </jats:sec><jats:sec>                 <jats:title>Results</jats:title>                 <jats:p>At baseline, RTL of the evol<jats:italic>TP53</jats:italic> patients was scattered across the entire RTL spectrum observed in our CLL cohort. RTL changed in the follow-up samples of 16/26 (62%) evol<jats:italic>TP53</jats:italic> cases, inclining to reach intermediate RTL values, i.e., longer telomeres shortened compared to baseline while shorter ones prolonged. For the first time we show that <jats:italic>TP53</jats:italic> clonal shifts are linked to RTL change, including unexpected RTL prolongation. We further investigated parameters associated with RTL changes. Unstable telomeres were significantly more frequent among younger patients (<jats:italic>P</jats:italic> = 0.032). Shorter telomeres were associated with decreased activity of the B-cell receptor signaling components p-ERK1/2, p-ZAP-70/SYK, and p-NFκB (<jats:italic>P</jats:italic> = 0.04, <jats:italic>P</jats:italic> = 0.01, and <jats:italic>P</jats:italic> = 0.02, respectively).</jats:p>               </jats:sec><jats:sec>                 <jats:title>Conclusions</jats:title>                 <jats:p>Our study revealed that changes of telomere length reflect evolution in leukemic subclone proportion, and are associated with specific clinico-biological features of the explored cohort.</jats:p>               </jats:sec>

section_content: Chronic lymphocytic leukemia (CLL) is a biologically heterogeneous entity with variable clinical outcomes [1] [2] [3]. CLL tumor cell mass can be generally composed of multiple genomically distinct subclones endowed with different capacities to proliferate and survive, and their relative composition evolves during the disease course, especially under therapy pressure [4, 5]. This process, generally known as clonal evolution, crucially drives CLL disease progression. Nowadays, many established prognostic biomarkers are used to stratify CLL patients into clinically distinct subgroups. Among these, the aberrations affecting the TP53 gene have a particularly adverse effect on patient prognosis and chemotherapy response [6, 7]. TP53 defects frequently arise in CLL relapse after chemoimmunotherapy despite being undetectable at the time of diagnosis [8, 9] and are strongly associated with chemotherapy resistance [10, 11]. Therefore, the patients with aberrant TP53 should be treated with novel therapeutics based on BCR signaling or bcl-2 inhibition [12] [13] [14] [15]. 

Telomeres are complex repetitive DNA sequences bound by multiple interacting proteins, which protect ends of linear chromosomes against end-to-end chromosomal fusion and degradation. Human telomeres are composed of ∼4-11 kb long hexameric 5 '-TTA GGG -3 ' repeats in lymphocytes [16, 17], and their maintenance strategies are cell-type specific. While telomeres of somatic cells shorten with each cell division, germ cells are capable of maintaining stable telomere length using telomerase to sustain cell immortality [18]. However, even in stem cells, telomerase activity is limited and when eventually exhausted, it causes senescence and loss of ability to renew tissues [19, 20]. Of note, telomeres of normal B cells are naturally prolonged by telomerase during their activation in a germinal center of secondary lymphoid tissues [21]. 

Cancer cells are also able to activate telomerase to overcome short telomeres to ensure cell immortality [22]. In CLL, short telomeres have been associated with markers of inferior prognosis, such as unmutated IGHV genes, increased genomic complexity, 11q deletion/mutated ATM gene, and 17p deletion/mutated TP53 gene [23] [24] [25] [26] [27] [28] [29] [30], with short time to first treatment [31] and short overall survival (OS) [17, 27, 32, 33]. Telomere length is considered stable in most CLL cases; however, its changes have been observed in a subset of follow-up samples, especially in those undergoing clonal evolution [27, 34]. Furthermore, the level of telomerase activity in CLL cells is similar to its level in normal B lymphocytes, although a subgroup of IGHVmutated CLL cells had a significantly lower telomerase activity than healthy B cells [23]. In some studies, high telomerase activity was observed and correlated with poor prognosis [31, 35]. 

This present study follows up our previous efforts focused on investigating clonal evolution of TP53 aberrations [8] and on the associations of ATM aberrations with short telomere length [29] in CLL. Here we explored a well-characterized cohort of 198 CLL patients to assess associations between CLL genomic features and telomere length. Sequential samples from 26 patients, who underwent clonal evolution of TP53 mutations (evolTP53 group) over the disease course, revealed for the first time an interplay between telomere length of subclones, BCR signaling activity, and CLL clinico-biological features.
section_name: Background

section_name: Methods

section_content: Two hundred and thirty-five CLL samples from 198 CLL patients were included in the present study; clinical and laboratory data were available for all tested cases (Tables 1 and 2 ). CLL cells were separated from peripheral blood with high purity (> 98%) as we described previously [36] using Ficoll-Paque PLUS (GE Healthcare, USA) combined with RosetteSep Human B cell Enrichment Kit and CD3 + Depletion Kit (StemCell Technologies, Canada) to deplete non-B cells. 

The study was approved by the Ethical Committee of the University Hospital Brno and all patients provided informed consent to the research use of biological material and clinical data in agreement with the Helsinki declaration.
section_name: Primary CLL samples

section_content: IGHV gene somatic hypermutation status was analyzed as described previously [37, 38]. 

Exonic regions and splice sites of the TP53 and ATM genes were amplified using in-house designed primers [8, 29]. Sequencing libraries were prepared using Illumina Nextera XT library preparation kit (Illumina, USA) and sequenced using Illumina MiSeq or NextSeq machines in paired-end mode (2 × 150 bp). The results were analyzed with an in-house bioinformatic pipeline [8, 29]. The cut-off variant allelic frequency (VAF) of 10% recommended in routine clinical diagnostics [39] was used to identify TP53 and ATM mutated cases. However, for the characterization of the evolTP53 group, ultra-deep sequencing with a cut-off 0. 1% VAF was used to confirm TP53 clonal evolution [8]. 

Genomic aberrations were assessed by FISH for del(11q), del(13q), del(17p), and trisomy 12 (XL CLL Probe Kit, MetaSystems). Genome-wide copy number and heterozygosity profiles were obtained using CytoScan HD arrays (Thermo Fisher Scientific). Metaphase cytogenetics was performed on peripheral blood samples (n = 169) stimulated with IL 2 and DSP30 CpG Oligonucleotide. Complex karyotype and highly complex karyotype were defined as ≥ 3 or ≥ 5 abnormalities, respectively.
section_name: Molecular characterization of the cohort

section_content: Telomere length was measured by the quantitative polymerase chain reaction (qPCR) method described previously [29, 40] with certain modifications (see Supplementary Material) using primers specific for telomeric repetition [41]. HBG3 and HBG4 primers [42] were used for the amplification of β-globin (HBB) gene, serving as a single-copy gene control. 

Each sample was measured in triplicates for both RTL and HBB, the reactions were assembled by Freedom EVO automated system (Tecan, Switzerland) in 384well plates, and qPCR was performed by QuantStudio 12 k Flex (Applied Biosystems, USA). Each reaction contained 15 ng DNA, 100 nM F and R primers, 5 mM DTT, 10 mM Tris-HCL (pH 8), 130 mM DMSO, 1. 5 µM BSA, 50 nM ROX in the final volume 10 µl. Quantitative PCR was performed on QuantStudio 12 k Flex (Applied Biosystems) with the cycling conditions as follows: 10 min 95 C, 40x (20 s 95 °C, 1 min 56 C). The RTL was calculated as -ΔCT for CLL samples divided by -ΔCT for control human DNA sample (G3041, Promega, USA) isolated from leukocytes of healthy donors. The -ΔCT values represent the difference between average Ct values of telomere length and HBB gene. The cut-off value for RTL change in consecutive CLL cases was 5% difference of RTL between baseline and follow-up, which corresponded to the complete loss of telomeres on two chromosomes (calculated approximately as 46 chromosomes with 2 telomeric regions each = 92 telomeric regions, 5% loss corresponds to loss of ~ 4 telomeric regions).
section_name: Relative telomere length (RTL) measurement

section_content: Reverse transcription was carried out with SuperScript II Reverse Transcriptase (Thermo Fisher Scientific). One ng of reverse-transcribed cDNA was used for real-time PCR containing 2 × ABsolute QPCR Mix (low ROX, Thermo Fisher Scientific) and either hTERT or GAPDH TaqMan assay (Hs00972650_m1 for hTERT, Hs03929097_g1 for GAPDH, Thermo Fisher Scientific) following the manufacturer protocol. PCR was performed on QuantStudio 12 k Flex (Applied Biosystems) with the cycling conditions: 15 min 95 °C, 40x (15 s 95 °C, 1 min 60 °C).
section_name: Human telomerase reverse transcriptase (hTERT) expression

section_content: Fresh frozen CLL cells were thawed, fixed, and permeabilized with 1. 6% paraformaldehyde and 100% methanol. Supplementary Table 1 shows that the samples were incubated with antibodies against phosphorylated AKT, BTK/ITK, ERK1/2, IKKα/β, NF-κB, p38, PLCγ2, and ZAP-70/SYK at 20 °C for 60 min, and measured using FACSVerse flow cytometer (BD Biosciences). Median fluorescent intensity (MFI) value, representing the amount of basal phosphorylation of BCR signaling pathway components, was used for further analysis.
section_name: BCR signaling assessment

section_content: All statistical analyses were performed using software JASP (Version 0. 11. 1) [43]. Specific statistical tests used for different study variables are described in the figure legends. T-tests were Student's two-tailed, RTLs of paired samples were evaluated by Wilcoxon signed rank test, correlations were counted as Pearson's r. P values < 0. 05 were considered statistically significant.
section_name: Statistical analysis

section_content: For the original data and detailed protocols please contact the corresponding author.
section_name: Data sharing statement

section_name: Results

section_content: In order to comprehensively assess the association between genomic features and telomere length in CLL, we studied a set of 198 clinically and genetically wellcharacterized CLL cases (Table 1 ). The majority of samples (n = 153; 77. 3%) was collected prior to first treatment, and for patients with multiple available samples, only the first sample was included here. Forty samples (20. 2%) carried a TP53 mutation above 10% VAF (with concomitant del(17p) in 23/40 cases). 

The median RTL of the entire CLL cohort was 0. 81 (range 0. 46 -1. 25; Supplementary Figure 1 ), meaning that telomere length in CLL cells was generally shorter than in the control mixed human DNA sample. Additionally, we observed significantly longer RTLs in patients at Rai and Binet early disease stages at diagnosis (Supplementary Figure 2 ). To assess the impact of RTL on OS, we split the cohort into two groups with RTLs above and below median RTL (median RTL = 0. 84) and detected statistically significant difference between the groups (median OS 119 and 91 months, respectively; P = 0. 0004; Supplementary Figure 3A ). When focusing the analysis on distinct patient subgroups with genomic aberrations, we found statistically significant associations of RTL with unmutated IGHV somatic hypermutation status, recurrent cytogenetic aberrations in hierarchical order [44], complex karyotype (as defined by 3 abnormalities), and mutations in TP53 and ATM genes (Supplementary Table 2 ), which was in line with previously published results [23-26, 28, 33] thus confirming the representativeness of the patient cohort. Since 22. 7% of samples were obtained after treatment, we tested whether RTL is dependent on treatment administration. In agreement with Mansouri et al. [34], we observed that previous therapy was not significantly associated with RTL (Supplementary Table 2 ). Nonetheless, the treated and untreated cohorts were not fully biologically and clinically comparable.
section_name: Telomere length associates with clinical features of the disease

section_content: We focused in detail on the associations between telomere length and TP53 mutational status in our CLL cohort. First, we evaluated the impact of RTL on patient OS in the subgroup with TP53 mutations (n = 40). We observed that OS of this subgroup divided by the RTL value above and below median was statistically different favoring cases with longer RTL (median RTL = 0. 78; median OS 86 and 55 months, respectively; P = 0. 046; Supplementary Figure 3B ). Next, available retrospective ultra-deep NGS data [39] were used to divide the patients from the basic cohort into three groups: (I) patients with wild-type TP53 (wtTP53; n = 135; mutation status has not changed throughout the observation), (II) patients with mutated TP53 above 10% VAF (mutTP53; n = 40; median VAF = 57. 5%), and (III) patients with TP53 mutation expanding above 10% VAF later in a disease course (n = 23). All group III patients except one expanded a TP53 mutation after therapy. Median TP53 VAF of group III was 0. 2% (0% -5. 9%) at baseline and 43. 4% (12. 1%-100%) at follow-up, corresponding to median fold change of 211. Figure 1 summarizes the RTL dynamics in the context of TP53 mutation clonal evolution. When we compared RTLs between the patients with wtTP53 (group I) and mutTP53 (group II), we found that the mutTP53 samples had significantly shorter telomeres (median RTL 0. 87 and 0. 78, respectively; P < 0. 001; Fig. 1A ). Interestingly, RTLs of group III (with TP53 mutation expanding above 10% VAF later) were intermediate (median RTL 0. 83; Fig. 1A ) and they did not differ significantly either from group I (P = 0. 25) or group II (P = 0. 10). 

Apart from the above-mentioned 23 cases from group III manifesting TP53 mutation expansion, we identified another 3 cases from group II (mutTP53) who experienced replacement of a certain TP53 mutation by another one after therapy. Thus, 26 cases out of the entire cohort experienced TP53 mutation evolution (evolTP53 group) during the disease course (Table 2, Supplementary Table 3 ). In order to discover whether mutTP53 evolution showed any association with telomere length, we studied pre-and post-evolution samples. All evolTP53 cases were analyzed in at least 2 time points. The median follow-up between the pre-and post-evolution samplings was 38. 5 months (range 18 -122 months). The median RTL was 0. 82 (range 0. 52 -1. 25) before and 0. 78 (range 0. 65 -0. 98) after the mutTP53 evolution (Fig. 1B ). When focusing on individual patients, the RTL change course was inconsistent. In detail, we observed that telomeres became shorter, remained stable, or became longer in 10/26 (38. 5%), 10/26 (38. 5%), and 6/26 (23. 1%) patients, respectively (Fig. 1C ). Interestingly, short telomeres became longer over time, while in contrast, long telomeres became shorter (Pearson r = -0. 75, P < 0. 001; Fig. 1C ). Additionally, for 6 and 2 of the 26 patients, 3 and 4 time points were analyzed, respectively, to follow further RTL development. Evolutions of genomic aberrations related to changes of RTLs in two CLL cases with the shortest basal telomere length are illustrated in Figs. 2A and 2B. Importantly, the apparent correlation between RTL and genomic abnormalities indicates that telomere length was likely a feature of CLL clones, not a hallmark of the entire CLL population changing in time.
section_name: Clonal evolution of TP53 defects affects telomere length

section_content: Since the changes of telomere length could potentially be caused by deregulated telomerase activity [31, 45], we tested hTERT expression in the 26 evolTP53 patients. However, only weak or null hTERT expression was detected. Supplementary Figure 4 and Supplementary Table 3 show that we did not find any correlation with telomere length or its change.
section_name: Association of RTL dynamics with clinico-biological features in cases evolving TP53 mutations hTERT expression

section_content: Intensive cell division, which could possibly cause RTL shortening, can be estimated from LDT [24]. Supplementary Figure 5 and Supplementary Table 3 show the correlation of LTD of the evolTP53 cases with the change of RTL between the baseline and follow-up. Short doubling time was observed in available baseline samples (median 6. 5 months, range 0. 5 -39. 6; n = 21) as well as in followup samples (median 2. 6 months; range 0. 4 -81. 5; n = 25). Nevertheless, no significant association of LDT and change of RTL was found.
section_name: Lymphocyte doubling time (LDT)

section_content: To study the relation between RTL dynamics and the activity of BCR signaling, we applied flow cytometry to detect basal phosphorylation of selected intracellular markers in CLL cells obtained prior and after mutTP53 expansion (n = 11 and 10, respectively, of these 7 samples were paired; Supplementary Table 3 ). Figure 3 illustrate the correlation of RTL and BCR signaling activity. Significantly longer telomeres were detected prior as well as after mutTP53 expansion in cases with higher phosphorylation of ERK1/2 (prior mutTP53 expansion: Pearson r = 0. 63; P = 0. 04; Fig. 3A ; after mutTP53 expansion: Pearson r = 0. 66; P = 0. 04; Fig. 3C ; Supplementary Table 3 ) and ZAP-70/SYK (prior mutTP53 expansion: Pearson r = 0. 60; P = 0. 05; Fig. 3B ; after mutTP53 expansion: Pearson r = 0. 79; P = 0. 01; Fig. 3D ; Supplementary Table 3 ). Additionally, significantly longer telomeres after mutTP53 expansion were detected in cases with higher phosphorylation of NF-κB (Pearson r = 0. 70; P = 0. 02; Fig. 3E ; Supplementary Table
section_name: BCR signaling activity

section_content: All patients with changing RTLs (both shortening and extending) could be considered as one group of patients with impaired telomere maintenance resulting in unstable RTLs during the disease course. Figure 4 shows that when we examined all clinical and biological features (as listed in Supplementary Table 3 ) of the patients with unstable RTLs (n = 16), we observed that they were statistically significantly younger at diagnosis compared to the patients with stable telomeres (median 55. 5 vs. 63. 7; P = 0. 032).
section_name: Clinical and biological features

section_content: Dysregulation of telomeres is common in many cancer types and represents a means for overcoming senescence and achieving uncontrolled cell growth [46]. In CLL, short telomere length has been linked to poor prognosis [23, [31] [32] [33] 35]. Here, we used a previously described qPCR method [40, 41] to study relative telomere length (RTL) in a large cohort of 198 CLL patients. The use of pure samples containing > 98% of CLL cells enabled us to follow the RTL very accurately. We confirmed that RTL in CLL patients was generally shorter compared to the control DNA pooled from healthy donors. Short telomeres were associated with mutations in ATM gene as we reported previously [29] and with other markers of inferior prognosis such as unmutated IGHV, complex karyotype, del(17p) and del(11q), and mutations in TP53 gene in concordance with published results [23-28, 30, 31]. The shortest detected telomeres were in patients with mutated ATM gene, while the longest in patients with 13q deletion as a sole cytogenetic aberration. Additionally, our study provided an insight into the relationship between clonal evolution of CLL cells with TP53 gene aberrations and their telomere length. Within our cohort of 26 evolTP53 cases evolving TP53 mutations, we measured RTL at several time points. We demonstrated that the RTL can evolve diversely with the clonal evolution of genomic abnormalities in CLL. 

Clonal evolution is an adverse event during the course of CLL and other cancers, which gives rise to subclones of cells with novel aberrations and features supporting increased fitness. In CLL, TP53-aberrant subclones tend to expand in disease relapse frequently and have a crucial impact on CLL patients' prognosis and survival as we reported previously [8, 47]. Jebaraj et al. [27] hypothesized that telomere shortening preceded the appearance of TP53 and ATM defects, and could provide selective pressure on CLL subclones favoring those with abilities of uncontrolled cell growth [27]. In our study, patients with TP53 mutation expanding above 10% VAF later in a disease course did show a tendency towards shorter telomeres prior to TP53 clonal expansion when compared to telomeres of the wtTP53 cases. Thus, it was not possible to anticipate a forthcoming TP53 mutation expansion in individual CLL cases based on the RTL. 

In our evolTP53 cohort, we performed measurements of paired samples and showed RTL changes between baseline and follow-up samples. Mansouri et al. [34] have found a high concordance between RTLs at diagnostic and follow-up CLL samples [34] suggesting that the telomere length does not change during the disease course. Recent study of Jebaraj et al. [27] focused on telomere length in 41 CLL patients without clonal evolution and 24 CLL patients with clonal evolution between samplings (9/24 had clonal evolution of TP53 mutation detected by whole-exome sequencing). They identified significant telomere shortening between samplings only within patients with clonal evolution. In contrast, in our study specifically focused on patients with expanding TP53 mutations, we detected stable RTL only in 38. 6% of cases, whereas in the majority of evolTP53 cases RTL alterations were noted. Interestingly, 4/26 patients with mutated ATM gene had significantly shorter RTL at baseline (P = 0. 009), however, ATM mutation in 3/26 patients was not significantly associated with shorter RTL at follow-up (P = 0. 061; data not shown). 

To the best of our knowledge, our report is the first in CLL showing prolongation of RTL over time in a considerable proportion (23%) of evolTP53 patients. RTL prolongation was especially apparent in samples with very short telomeres detected at baseline. In these cases, increased telomere length was frequently found in the corresponding follow-up samples, suggesting selective pressure on subclones with longer telomeres with presumably better fitness. In this context, the prognosis of CLL cases with RTL prolongation remains to be evaluated. It seems plausible that these cases, while having longer telomeres suggesting a better prognosis, probably overcame a critically short telomere crisis and gained a wide range of new aberrations, which is typical of inferior prognosis. Nonetheless, longer telomeres in baseline samples were not necessarily maintained at follow-up. Our findings, thus, highlight an unexpected consequence of clonal evolution leading to both shorter or longer telomeres, associated with changing proportions of CLL subclones with defined genomic abnormalities. A larger study, specifically addressing the fate of individual clones carrying certain telomere lengths and distinct cytogenetic abnormalities, is needed to elucidate the causality between the telomere length, genomic features, and a clone fitness. 

To uncover the molecular mechanisms behind the RTL changes in our cohort, we studied two possible causal processes -hTERT expression and cellular proliferation. We did not confirm any clear role of hTERT in CLL telomere maintenance in concordance with others [23, 48]. Although the telomere shortening rate was previously attributed to a rapid proliferation rate [23, 31], neither the cellular proliferation, represented by lymphocyte doubling time, could explain the changing RTLs in our samples. Telomerase activity is highly restricted in normal somatic cells and was detected in normal B-lymphocytes in germinal centers [21]. CLL cells are well known to circulate between blood and lymph nodes [49]. Possibly telomerase activity in CLL could be restricted to the time when cells are staying in the supportive microenvironment of lymph nodes. Since our CLL samples are isolated from peripheral blood, we cannot anticipate the telomerase activity of CLL cells located in lymph nodes. Interestingly, we found that patients with unstable telomeres, both shortening or prolonging, with concurrent evolution of TP53 mutations, have a significantly lower age at diagnosis. We have not found any association of unstable telomeres with clinical or biological features. Some more complex factors of CLL disease nature in younger patients could be involved. 

Finally, for the first time, we studied the interplay between tonic BCR signaling and telomere length maintenance in CLL. CLL cells depend on BCR signaling, supporting their survival, proliferation [50], therapy resistance, and gene expression [50] [51] [52]. BCR pathway is successfully targeted by novel inhibitor agents introduced in clinical practice in the last decade [12] [13] [14]. Using a targeted selection of phosphorylated kinases, we aimed to analyze different branches of BCR signaling to uncover possible associations with RTL changes. Patients with shorter telomeres had lower basal phosphorylation of ZAP-70/SYK and ERK1/2 in the baseline samples and lower basal phosphorylation of ZAP-70/ SYK, ERK1/2, and p-NF-κB in the follow-up samples, indicating weakened BCR signaling at both time points. Other phosphorylated signaling components, especially BTK/ITK and PLC-γ2, showed strong, although not significant, trends towards lower basal phosphorylation in patients with shorter telomeres. This could anticipate that CLL cells with shorter telomeres are less dependent on survival and proliferative signals delivered to the cells by tonic BCR signaling.
section_name: Discussion

section_content: To summarize, we have performed an assessment of the telomere length in a well-defined cohort of CLL cases. Apart from confirming known associations between genomic features and telomere shortening, we focused on the relationship between telomere length and TP53 mutation evolution. We discovered that long telomeres can shorten, while short telomeres can, in some cases, even extend their length during the clonal evolution process. Thus, our study shows that shifts in subclone equilibriums can be reflected in telomere length changes, while telomere length per se is likely not the only factor affecting the survival capacity of a clone expanding in the evolution process.
section_name: Conclusions

section_name: Acknowledgements

section_content: The authors would like to thank Petra Sedlakova, Jiri Baloun, Boris Tichy, and Filip Pardy for technical assistance and Veronika Mancikova for editing the manuscript. H. Olbertova is a PhD candidate at Masaryk University. This work is submitted in partial fulfillment of the requirement for the PhD.
section_name: Acknowledgements

section_name: Funding

section_content: This work was supported by the projects provided by the Ministry of Health of the Czech Republic RVO 65269705, AZV NU21-08-00237 and NV19-03-00091, and by the Czech Science Foundation (GACR) 19-15737S and 19-11299S. The work was supported from European Regional Development Fund-Project " NCLG II. " (No. CZ. 02. 1. 01/0. 0/0. 0/18_046/0015515 ) and research infrastructures NCMG (ID number LM2018132 ) and EATRIS-CZ (ID number LM2018133 ), both funded by MEYS CR. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. Ministerstvo Zdravotnictví Ceské Republiky, MH-CZ RVO 65269705, AZV NU21-08-00237, AZV NV19-03-00091, Grantová Agentura České Republiky, GACR 19-15737S, GACR 19-11299S, Technologická Agentura České Republiky, LM2018132, LM2018133, Ministerstvo Školství,Mládeže a Tělovýchovy, CZ. 02. 1. 01/0. 0/0. 0/18_046/0015515
section_name: Funding

section_name: Availability of data and materials

section_content: The data used during the current study are available from the corresponding author on reasonable request.
section_name: Availability of data and materials

section_name: Supplementary Information

section_content: The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12885-022-09221-z. 

Additional file 1: Supplementary Table 1. Antibodies used for BCR signaling activity assessment. Supplementary Table 2. Results of relative telomere length (RTL) quantification in the basic CLL cohort (n = 198). Supplementary Table 3. TP53 mutation evolution cohort -overview of clinico-biological features and results in baseline and follow-up samples, PART A Abbreviations: ALZ -alemtuzumab; BR -bendamustine and rituximab; CLB -chlorambucil; CR -cyclophosphamide and rituximab; deldeletion; F -female; FCR -fludarabine, cyclophosphamide, and rituximab; FISH -fluorescence in situ hybridization; hyper CVAD -hyperfractionated cyclophosphamide, vincristine, doxorubicin, and dexamethasone; IGHV -immunoglobulin heavy chain variable region; M -male; M -mutated; mo -months; U -unmutated; yrs -years. Supplementary Table 3.
section_name: Supplementary Information

section_content: The whole study and written informed consent obtained from all patients were approved by the Ethical Committee of University Hospital Brno in concordance with the Declaration of Helsinki.
section_name: Declarations Ethics approval and consent to participate

section_content: Written consent was obtained from all patients.
section_name: Consent for publication

section_content: The authors declare that they have no competing interests.
section_name: Competing interests

section_content: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
section_name: Publisher's Note