1. Introduction
B-cell chronic lymphocytic leukemia (B-CLL) stands as the most common form of leukemia in the Western world [
1], with an annual incidence rate of 3.1 and 6.1 per 100,000 in females and males, respectively [
2]. Remarkable advancements in treatment and diagnosis have elevated the 5-year survival rate from 77.3% in 2009 to 87.9% today [
3,
4]. The median age at which B-CLL is diagnosed is 72 years [
5,
6]. Characterized by the accumulation of CD5+ B lymphocytes in the blood, spleen, and lymph nodes [
7]. B-CLL exhibits CD5 expression in leukemic B cells, a marker typically found in T cells but uniquely present in B-CLL cells, distinguishing them from normal B cells [
8]. Several genes, including
TP53,
ATM,
MYD88, and those involved in Notch signaling, inflammatory pathways, B cell receptor signaling, and others, have been implicated in B-CLL development [
9]. Additionally, genes regulating MAPK-ERK and MYC pathways have been connected to the disease [
10]. Recurrent deletions in chromosomes 8p, 11q, 13q, and 17p, as well as trisomy of chromosome 12, have also been reported [
9]. A 2022 study identified 82 additional putative driver mutations, each occurring in less than 2% of analyzed patients [
11]. Notably, the currently known mutated genes mainly affect fewer than 10% of B-CLL patients, suggesting other alterations contribute to the disease. Prognostic biomarkers have become essential for assessing disease prognosis, encompassing host factors such as gender and age, cell marker expression levels (CD38, ZAP70, and CD49d), serological values (β2-microglobulin, LDH), genetic alterations (deletion of chromosome arms 11q, 13q, 17p, TP53 gene mutation, and trisomy 12), and the mutational status of the IGHV gene [
12,
13,
14].
Traditionally, first-line therapy has relied on chemoimmunotherapy, such as FCR (fludarabine, cyclophosphamide, and rituximab), chlorambucil, or bendamustine with anti-CD20 antibodies (rituximab, obinutuzumab). However, long-term disease control with FCR is primarily achieved in patients with mutated IGHV genes. Recent years have witnessed a shift in treatment options, thanks to the development of kinase inhibitors targeting BCR signaling, such as ibrutinib (Bruton’s tyrosine kinase (BTK)) and idelalisib (phosphoinositide 3-kinase (PI3K)δ) [
15]. Second-generation, more specific BTK inhibitors (acalabrutinib), venetoclax (a specific inhibitor of the anti-apoptotic factor BCL2), alone or in combination with anti-CD20 antibodies, have further improved patient outcomes compared to traditional chemotherapy [
16]. For instance, venetoclax treatment induces cell death independently of
TP53 mutations and/or del(17p) [
17]. Ibrutinib has been highly successful in B-CLL treatment and is now a leading therapy, either alone or in combination with other chemotherapeutic or immunotherapeutic agents [
18].
The widely used mouse model for B-CLL is the Eμ-TCL1 model, where the T-cell leukemia oncogene TCL1 is inserted under the control of the immunoglobulin heavy chain variable region promoter and immunoglobulin heavy chain enhancer (Eμ), resulting in B-CLL-like disease development at a late age (13 to 18 months) [
19]. However, this model has drawbacks, including the extended period required for full disease development and its inability to fully represent the recurrent human B-CLL disease. Another approach involves engrafting the B-CLL cell line MEC-1 into Rag2-/-γc-/- mice, which mimics aggressive human B-CLL but falls short in expressing CD5, a canonical marker seen in human B-CLL patients [
20]. Engraftment of B-CLL patient-derived B cells into NSG mice, with or without co-injection of polyclonally-activated autologous T cells pre-stimulated in vitro, has been successful but limits the window for B-CLL biology and potential treatment testing to just 63 days [
21]. Deletion of the DLEU2/miR-15a/16-1 locus in mice, encoded on chromosome arm 13q, leads to an indolent disease mirroring human B-CLL characteristics, albeit with a disease penetrance of approximately 50% [
22]. These and other approaches to the generation of mouse models of B-CLL are thoroughly reviewed in [
23]. While these models are currently used in B-CLL research, there’s a growing need for newer models that better capture the full heterogeneity of B-CLL biology. Hence, it is crucial in mouse model development to prioritize full B-CLL penetrance, eliminating the need for xenografts, and promoting disease development in the early stages of mouse life, thus avoiding prolonged waiting times for disease study.
The RAS protein family includes numerous small guanosine triphosphate hydrolases (GTPases), some of which have been found to be mutated and implicated in various human cancers [
24]. The classical RAS members, K-RAS, H-RAS, and N-RAS, were discovered in the late 1960s, 1970s, and early 1980s [
25]. RAS-related proteins (R-RAS) share several domains and regulatory factors with classical RAS proteins, mediating their activation and inactivation cycles, including guanine exchange factors (GEF) and GTPase activating proteins (GAP) [
26]. Mutation of R-RAS2 in analogous residues to G12V and Q61L in classical RAS proteins has induced comparable cell transformation in culture and tumor growth in mice [
27,
28]. A key distinction is the high intrinsic nucleotide exchange activity of R-RAS2, allowing it to exchange GDP for GTP independently of specific GEF proteins at rates similar to those observed with H-RAS after GEF addition [
29].
R-RAS2 plays diverse roles in the body. It regulates platelet activation by interacting with the glycoprotein VI-ITAM-containing collagen receptor [
30]. R-RAS2 regulates Schwann cell migration [
31]. and mammary gland development [
32]. Overexpression of R-RAS2 protein has been observed in various cancer types, including oral squamous cell carcinoma [
33], esophageal tumors [
34], hepatocellular carcinoma, [
35] and highly aggressive skin cancer [
36]. Mutation or overexpression of
RRAS2 has led to the transformation of breast cell lines [
37] [
38]. R-RAS2 with the Q72L mutation, associated with Noonan Syndrome, has been found to trigger the formation of various tumors [
39,
40]. Recent reports have shown that R-RAS2 harboring the Q72L mutation is a potent oncogenic driver that triggers the formation of a wide variety of tumors (ovarian cystadenomas, T-ALL, etc.) [
41]. Moreover, other mutations (G23V/A/C/S, G24D/C/V, A70T, and Q72H) in human cancer have shown transforming potential when expressed in immortal cell lines [
42].
Previous research from our group has demonstrated that R-RAS2 directly interacts with both B-cell and T-cell receptors (BCR and TCR) through their immunoreceptor tyrosine-based activation motif (ITAM), predominantly in the inactive GDP-bound state of R-RAS2. This interaction provides tonic survival signals [
43]. R-RAS2 also regulates the internalization of the TCR following immune synapse formation, but only in its wild-type (WT) form [
44]. In B cells, R-RAS2 plays a crucial role in governing the proper formation of germinal centers by influencing B-cell metabolism [
45]. Through the analysis of human samples from B-CLL patients and utilizing our Rosa26-
RRAS2fl/fl-mb1-Cre (and Sox2-Cre) mouse models, we have recently demonstrated that the overexpression of RRAS2 drives the development of B-CLL [
46]. In our current study, we treated Rosa26-
RRAS2fl/fl-mb1-Cre with the clinically used drugs, ibrutinib and venetoclax [
16]. Notably, the leukemic cell population, especially in the case of ibrutinib treatment, showed significant regression. This suggests the potential utility of this mouse model as a valuable tool for exploring new therapeutic approaches for B-CLL.
4. Discussion
There are limited mouse models that fully replicate B-CLL, with the Eμ-TCL1 being the most commonly used. However, this model presents challenges, such as a long disease development period (3 to 5 months for spleen involvement and over a year for blood disease) and the need for genetic engineering to insert the TCL1 oncogene [
19]. Alternative methods include xenografts using the B-CLL cell line MEC-1 [
20], co-transplantation of patient-derived B-CLL B cells with autologous pre-stimulated T cells [
21], or deletion of the DLEU2/miR-15a/16-1 locus in mice to mimic indolent B-CLL [
22].
In contrast, the Rosa26-
RRAS2fl/flxmb1-Cre mouse model [
46] offers a novel approach. It exhibits early onset of CD19+CD5+ leukemic cells with a 100% disease penetrance, concurrent with splenomegaly and progressive blood lymphocytosis leading to reduced lifespan [
46]. In this study, we validate the Rosa26-
RRAS2fl/flxmb1-Cre model for preclinical testing of B-CLL therapies. Treatment of these mice with venetoclax and, notably, ibrutinib leads to a reduction in CD19+B220
lowCD5+ leukemic cells in the blood and spleen.
However, it’s worth noting that the treatments caused the mortality of 3 out of 9 mice in the ibrutinib group and 4 out of 9 in the venetoclax group by day 18, reflecting side effects seen in some B-CLL patients [
51]. In one case, a mouse treated with venetoclax had to be sacrificed due to severe lymphocyte loss. Venetoclax also significantly reduced the population of immature B cells in the bone marrow, which may explain the loss of B cells in some cases. Such side effects emphasize the need for therapies that balance efficacy with reduced side effects.
The splenomegaly observed in the Rosa26-RRAS2fl/flxmb1-Cre mice mirrors the human condition [
5]. After one month of ibrutinib and venetoclax treatment, splenomegaly was significantly reduced, consistent with the effects seen in B-CLL patients [
53,
54]. These observations highlight the relevance of this mouse model for studying B-CLL and its response to treatment. In future studies, gender-related effects on treatment efficacy should be evaluated, as B-CLL affects male patients approximately twice as often as females [
2].
An unexplored treatment option for B-CLL and other malignancies is the use of direct R-RAS2 inhibitors. As demonstrated previously [
46], overexpression of WT
RRAS2 drives B-CLL development. Therefore, targeting R-RAS2 directly is a promising avenue for further research. Direct inhibitors, administered alone or in combination with other therapies like ibrutinib, may enhance therapeutic efficacy and reduce the risk of ibrutinib resistance [
57,
58].
Figure 1.
Rosa26-RRAS2fl/flxmb1-Cre mice can go through a one-month long treatment with ibrutinib or venetoclax. a) Schematic representation of the treatment protocol used in Rosa26-RRAS2fl/flxmb1-Cre mice. b) Kaplan-Meier survival curve of the mice under treatment. n=9 mice per group started the experiment. c) Relative percentage of mouse weight evolution through the course of the treatment.
Figure 1.
Rosa26-RRAS2fl/flxmb1-Cre mice can go through a one-month long treatment with ibrutinib or venetoclax. a) Schematic representation of the treatment protocol used in Rosa26-RRAS2fl/flxmb1-Cre mice. b) Kaplan-Meier survival curve of the mice under treatment. n=9 mice per group started the experiment. c) Relative percentage of mouse weight evolution through the course of the treatment.
Figure 2.
Ibrutinib and venetoclax treatment reverse the splenomegaly and follicle enlargement induced by RRAS2 overexpression. a) Spleen weights of mice in the vehicle, ibrutinib and venetoclax treatment groups at the experiment endpoint, day 31. Two-tailed unpaired t test with Welch’s correction. b) Quantification of the number of CD19+ B cells in the spleens of the treated mice at the experiment endpoint. Two-tailed unpaired t test with Welch’s correction. Left, percentage of the total lymphocyte population; right, total numbers. c) Left, representative two-parameter flow cytometry plots of CD19 and CD5 expression in the spleens of mice in each of the three established groups. The CD19+CD5+ population is highlighted with a blue box. Right, quantification of the number of CD19+CD5+ B cells in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. d) Left, representative two-parameter flow cytometry plots of GFP expression vs side scatter (SSC) in the spleens of mice in each of the three established groups. The GFP+ population is highlighted with a blue box. Right, quantification of the number of GFP+ cells in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. e) Left, representative two-parameter flow cytometry plots of IgM and IgD expression within the total CD19+ B cell population. The follicular IgMlowIgD+ population is highlighted with a blue box. Right, quantification of the number of IgMlowIgD+ cells with the CD19+ gate in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total CD19+ population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. f) Representative hematoxylin and eosin stainings of the spleens of mice in the vehicle, ibrutinib and venetoclax groups at the experiment endpoint. Scale bars represent 400μm. In each of the three images, the black bars are illustrative of the diameters used to calculate follicle areas. g) Box and whiskers plot showing all points and median value of the quantification of follicle areas using the diameters illustrated in the images in f). Areas were calculated for all follicles visible in their entirety in three representative images per mouse. n=9 in the vehicle group, n=6 in the ibrutinib group, n=5 in the venetoclax group. Two-tailed unpaired t test with Welch’s correction. ns: not significant.
Figure 2.
Ibrutinib and venetoclax treatment reverse the splenomegaly and follicle enlargement induced by RRAS2 overexpression. a) Spleen weights of mice in the vehicle, ibrutinib and venetoclax treatment groups at the experiment endpoint, day 31. Two-tailed unpaired t test with Welch’s correction. b) Quantification of the number of CD19+ B cells in the spleens of the treated mice at the experiment endpoint. Two-tailed unpaired t test with Welch’s correction. Left, percentage of the total lymphocyte population; right, total numbers. c) Left, representative two-parameter flow cytometry plots of CD19 and CD5 expression in the spleens of mice in each of the three established groups. The CD19+CD5+ population is highlighted with a blue box. Right, quantification of the number of CD19+CD5+ B cells in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. d) Left, representative two-parameter flow cytometry plots of GFP expression vs side scatter (SSC) in the spleens of mice in each of the three established groups. The GFP+ population is highlighted with a blue box. Right, quantification of the number of GFP+ cells in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. e) Left, representative two-parameter flow cytometry plots of IgM and IgD expression within the total CD19+ B cell population. The follicular IgMlowIgD+ population is highlighted with a blue box. Right, quantification of the number of IgMlowIgD+ cells with the CD19+ gate in the spleens of the treated mice at the experiment endpoint. Left graph, percentage of the total CD19+ population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. f) Representative hematoxylin and eosin stainings of the spleens of mice in the vehicle, ibrutinib and venetoclax groups at the experiment endpoint. Scale bars represent 400μm. In each of the three images, the black bars are illustrative of the diameters used to calculate follicle areas. g) Box and whiskers plot showing all points and median value of the quantification of follicle areas using the diameters illustrated in the images in f). Areas were calculated for all follicles visible in their entirety in three representative images per mouse. n=9 in the vehicle group, n=6 in the ibrutinib group, n=5 in the venetoclax group. Two-tailed unpaired t test with Welch’s correction. ns: not significant.
Figure 3.
Venetoclax treatment hampers transition from pro-pre B to immature cell stage in the bone marrow. a) Left, representative two-parameter flow cytometry plots of GFP expression vs side scatter (SSC) in the bone marrow of mice in each of the three established groups. The GFP+ population is highlighted with a blue box. Right, quantification of the number of GFP+ cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. b) Left, representative two-parameter flow cytometry plots of CD19 and CD5 expression in the bone marrow of mice in each of the three established groups. The CD19+CD5+ population is highlighted with a blue box. Right, quantification of the number of CD19+CD5+ cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. c) Left, representative two-parameter flow cytometry plots of IgM and B220 expression within the total CD5- population in the bone marrow of mice in each of the three established groups. The immature B220hi IgM+ population is highlighted with a blue box. The pro-pre B cell population is highlighted with a pink box. Right, top graphs, quantification of the number of immature B cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total CD5- population; right graph, total numbers. Right, bottom graphs, quantification of the number of pro-pre B cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total CD5- population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. ns: not significant.
Figure 3.
Venetoclax treatment hampers transition from pro-pre B to immature cell stage in the bone marrow. a) Left, representative two-parameter flow cytometry plots of GFP expression vs side scatter (SSC) in the bone marrow of mice in each of the three established groups. The GFP+ population is highlighted with a blue box. Right, quantification of the number of GFP+ cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. b) Left, representative two-parameter flow cytometry plots of CD19 and CD5 expression in the bone marrow of mice in each of the three established groups. The CD19+CD5+ population is highlighted with a blue box. Right, quantification of the number of CD19+CD5+ cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total lymphocyte population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. c) Left, representative two-parameter flow cytometry plots of IgM and B220 expression within the total CD5- population in the bone marrow of mice in each of the three established groups. The immature B220hi IgM+ population is highlighted with a blue box. The pro-pre B cell population is highlighted with a pink box. Right, top graphs, quantification of the number of immature B cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total CD5- population; right graph, total numbers. Right, bottom graphs, quantification of the number of pro-pre B cells in the bone marrow of the treated mice at the experiment endpoint. Left graph, percentage of the total CD5- population; right graph, total numbers. Two-tailed unpaired t test with Welch’s correction. ns: not significant.
Figure 4.
Ibrutinib treatment reduces the circulating CD19+B220lowCD5+ leukemic population in the blood of Rosa26-RRAS2fl/flxmb1-Cre mice. a) Quantification of the total number of CD19+ B cells in the blood of mice at the start (0), halfway through the experiment (18) and the experiment endpoint (31) in the vehicle, ibrutinib and venetoclax groups. Column bars represent mean values ± SEM. b) Representative two-parameter flow cytometry plots of CD19 and CD5 expression in the blood of mice in each of the three established groups at the experiment endpoint. The CD19+CD5+ population is highlighted with a blue box. c) Quantification of the total number of CD19+CD5+ leukemic B cells in the blood of mice at the start (0), halfway through the experiment (18) and the experiment endpoint (31) in the vehicle, ibrutinib and venetoclax groups. d) Representative two-parameter flow cytometry plots of B220 and CD5 expression within the total CD19+ population at the endpoint of the experiment in the blood of the three established groups. The B220lowCD5+ population is highlighted with a blue box. e) B220lowCD5+ leukemic cells evolution, within the total CD19+ gate, in the blood of treated mice. Two-way ANOVA test. f) B220lowCD5+ population shown in e) relative evolution through the course of treatment. Two-way ANOVA test. g) B220+CD5- cells evolution, within the total CD19+ gate, in the blood of treated mice. Two-way ANOVA test. h) IgMlow IgD+ population, within the total CD19+ cells, in the blood of mice at the endpoint of the experiment. Two-tailed unpaired t test with Welch’s correction. ns: not significant.
Figure 4.
Ibrutinib treatment reduces the circulating CD19+B220lowCD5+ leukemic population in the blood of Rosa26-RRAS2fl/flxmb1-Cre mice. a) Quantification of the total number of CD19+ B cells in the blood of mice at the start (0), halfway through the experiment (18) and the experiment endpoint (31) in the vehicle, ibrutinib and venetoclax groups. Column bars represent mean values ± SEM. b) Representative two-parameter flow cytometry plots of CD19 and CD5 expression in the blood of mice in each of the three established groups at the experiment endpoint. The CD19+CD5+ population is highlighted with a blue box. c) Quantification of the total number of CD19+CD5+ leukemic B cells in the blood of mice at the start (0), halfway through the experiment (18) and the experiment endpoint (31) in the vehicle, ibrutinib and venetoclax groups. d) Representative two-parameter flow cytometry plots of B220 and CD5 expression within the total CD19+ population at the endpoint of the experiment in the blood of the three established groups. The B220lowCD5+ population is highlighted with a blue box. e) B220lowCD5+ leukemic cells evolution, within the total CD19+ gate, in the blood of treated mice. Two-way ANOVA test. f) B220lowCD5+ population shown in e) relative evolution through the course of treatment. Two-way ANOVA test. g) B220+CD5- cells evolution, within the total CD19+ gate, in the blood of treated mice. Two-way ANOVA test. h) IgMlow IgD+ population, within the total CD19+ cells, in the blood of mice at the endpoint of the experiment. Two-tailed unpaired t test with Welch’s correction. ns: not significant.