Next Article in Journal
Determination of a Tumor-Promoting Microenvironment in Recurrent Medulloblastoma: A Multi-Omics Study of Cerebrospinal Fluid
Next Article in Special Issue
Applying Tissue Slice Culture in Cancer Research—Insights from Preclinical Proton Radiotherapy
Previous Article in Journal
DFIQ, a Novel Quinoline Derivative, Shows Anticancer Potential by Inducing Apoptosis and Autophagy in NSCLC Cell and In Vivo Zebrafish Xenograft Models
Previous Article in Special Issue
Patient-Derived Xenograft Models of Pancreatic Cancer: Overview and Comparison with Other Types of Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of the Injection Site on Growth Characteristics, Phenotype and Sensitivity towards Cytarabine of Twenty Acute Leukaemia Patient-Derived Xenograft Models

1
Charles River Discovery Research Services Germany GmbH, Am Flughafen 12-14, 79108 Freiburg, Germany
2
Department of Medicine I, Faculty of Medicine, University of Freiburg, 79106 Freiburg, Germany
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(5), 1349; https://doi.org/10.3390/cancers12051349
Submission received: 28 April 2020 / Revised: 22 May 2020 / Accepted: 23 May 2020 / Published: 25 May 2020
(This article belongs to the Special Issue Patient-Derived Xenograft-Models in Cancer Research)

Abstract

:
Rodent models have contributed significantly to the understanding of haematological malignancies. One important model system in this context are patient-derived xenografts (PDX). In the current study, we examined 20 acute leukaemia PDX models for growth behaviour, infiltration in haemopoietic organs and sensitivity towards cytarabine. PDX were injected intratibially (i.t.), intrasplenicaly (i.s.) or subcutaneously (s.c.) into immune compromised mice. For 18/20 models the engraftment capacity was independent of the implantation site. Two models could exclusively be propagated in one or two specific settings. The implantation site did influence tumour growth kinetics as median overall survival differed within one model depending on the injection route. The infiltration pattern was similar in i.t. and i.s. models. In contrast to the s.c. implantation, only one model displayed circulating leukaemic cells outside of the locally growing tumour mass. Cytarabine was active in all four tested models. Nevertheless, the degree of sensitivity was specific for an individual model and implantation site. In summary, all three application routes turned out to be feasible for the propagation of PDX. Nevertheless, the distinct differences between the settings highlight the need for well characterized platforms to ensure the meaningful interpretation of data generated using those powerful tools.

1. Introduction

Myeloid malignancies in general are defined as clonal haemopoietic malfunction characterized by excessive proliferation, atypical self-renewal, and/or differentiation defects of haemopoietic stem cells (HSCs) and myeloid progenitor cells [1]. Leukaemia as a cancer category accounted for 440,000 new cases and more than 300,000 deaths worldwide in 2018 [2,3]. Acute myeloid leukaemia (AML) as one of the most devastating diseases in this category represents 1.2% of all new cancer cases in the U.S. in 2019 [4]. Our current understanding is that the onset of the disease is caused by different genetic and epigenetic changes in HSCs as well as functional changes in the bone marrow niche in the close vicinity of the emerging leukaemic cells. As a consequence of these developments a phenotypically distinct disease evolves that requires a specific treatment regimen depending on a multitude of risk factors [5]. Although some progress has been made specifically in the treatment of paediatric AML, the prognosis of adult AML specifically in the elderly patient has not improved significantly [6]. The high medical need is reflected in the increasing number of compounds in different stages of the development process: 1016 different drugs are currently investigated for their ability to tackle AML [7]. One important tool in this process are animal models as they are crucial for the selection of the most promising candidates at a late preclinical stage [8,9]. Furthermore, animal models are imminent for the understanding of the biology of the disease and thereby support the identification of potential new drug targets. Since the advent of the triple immunodeficient NOD-SCID-IL2rγnull (NSG) and related mouse strains like NOD/Shi-scid-IL-2Rγnull (NOG) the development of patient-derived leukaemia xenograft models was significantly facilitated [10,11]. Our group established a panel of 20 well-characterized acute leukaemia patient-derived xenograft models covering a broad range of different subtypes representing the molecular landscape of the disease. In the current study, we investigated the influence of the tumour microenvironment on different aspects of the leukaemic cell biology: the PDX were implanted via three different application routes and take rate, growth behaviour, expression of leukaemia-specific surface markers and sensitivity towards standard of care cytarabine was investigated. The implantation route had an impact explicitly on the tumour growth kinetics of an individual model. This was specifically interesting, as we identified this feature as a model intrinsic across different passages within one implantation scheme. Other model characteristics like the expression of three distinct surface markers and the take rate were model specific and stable across different engraftment sites. The sensitivity towards cytarabine was investigated in four models in a disseminated and the subcutaneous setting. In general, the sensitivity was not influenced by the implantation site. Nevertheless, using overall survival as a read-out, the sensitivity was influenced by the location of the tumour cells. In summary, the characteristics of an individual leukaemia PDX model are composed of the genetic make-up of the model itself and the tumour microenvironment modulated via the implantation route.

2. Results

2.1. Leukaemia PDX Retain the Cytogenetic Features of the Donor Patient and Mimic the Genomic Landscape of the Disease

Our group established and characterized 20 leukaemia PDX lines with the aim to develop a preclinical cancer model platform that recapitulates the human disease as close as possible to clinical reality. The collection included 18 acute myeloid leukaemia (AML), one acute promyelocytic leukaemia (APL) and one acute lymphoblastic leukaemia (ALL). The median age of the donor patients was 63 years (range 43–80 years), where 13 received no pre-treatment, 3 underwent peripheral blood stem cell transplantation (PBSCT) and 4 received other treatments mostly with a cytotoxic regimen. In the course of the establishment phase, we collected standard of care data for 19 out of 20 models. A molecular characterization was done on 16 models in total: eight models were analysed by whole exome sequencing (WES) and RNAsequencing (RNAseq), six models exclusively by WES and two models solely by RNAseq (see Table 1 and Tables S1 and S2).
During the model establishment phase the patient-derived blasts were injected in the bone marrow of recipient mice via intratibial (i.t.) injection. Mice were monitored at this time of first passage for up to 18 months daily for clinical signs of leukaemia and every two weeks for detection of human leukaemic cells in the peripheral blood by flow cytometry (FC) analysis. When disease onset was perceived, mice were sacrificed and cells from the bone marrow and spleen harvested. After the depletion of mouse cells from the suspension, human leukaemic cells were re-implanted into the bone marrow of recipient mice. This procedure was repeated until stable growth was observed. A model was defined as established when stable growth over at least four passages, regrowth from liquid nitrogen and accordance between genetic data from the PDX with the data from the patient were obtained. To investigate the influence of the tumour microenvironment growth behaviour, take rate and expression of distinct surface markers were determined in twenty different models applied by three different implantation routes. Beside the i.t. approach, which was used during the establishment phase, the intrasplenic (i.s.) route and the subcutaneous (s.c.) implantation were validated in more detail. The detailed implantation scheme is described in Figure 1a.
The molecular analyses of the individual models were not only used to verify the origin of the PDX from a specific patient but also to determine how representative the leukaemia PDX panel was for a respective patient population. The comparison of the frequency of AML specific mutations in the PDX panel and data of an AML patient population from TCGA revealed a highly significant correlation of both data sets (Pearson correlation coefficient 0.81, p < 0.00004). Thus, our PDX panel largely represents the molecular landscape of the human disease (Figure 1b).

2.2. The Overall Survival Time is a Model-Specific Feature which Evolves during the Establishment Phase of the Respective PDX.

Over the course of the development phase the individual PDX displayed differing overall survival times between the passages (Figure 2). Animals were sacrificed following stringent termination criteria enabling the comparison of those data across different models and passages. Within one model the first four passages differed significantly in terms of overall survival (OS, (Log-rank (Mantel–Cox) test, p < 0.005). From passage four on the OS time stabilized at a model specific value. The median OS time was characteristic for a specific established model and differed significantly between models. It ranged from 22 days (LEXF 2799) to 168 days (LEXF 2918) with a median OS time for the complete panel of 64.75 days (Table S3).

2.3. The Implantation Site Has an Impact on the Overall Survival Time in a Number of PDX Models but not the Complete Panel

The influence of the implantation site on OS was evaluated in 20 leukaemia PDX models. For all models the i.t. injection was compared with the s.c. approach, and for twelve models the i.s. implantation technique was additionally compared to the two other methods. In 14 out of 20 models the OS time was significantly influenced by the implantation site (Figure 3 and Table 2), whereas in six models no differences could be detected. We compared the median OS time across all models for which overall survival data in all tested settings was available (n = 12 models). In that case, the median OS of the i.s. implanted animals was 45 days (with three models showing no growth until end of the observation period). The median OS was 70 days for the i.t. injected mice and 67 days for the s.c. implanted animals. The differences were statistically significant for i.s. vs. s.c. and i.s. vs. i.t. In the case of the i.s. injection the outcome was influenced by the fact that three out of twelve tested models did not show any tumour growth during the observation period. Taking only models into account which showed successful engraftment of leukaemic cells the i.s. implanted cells grew significantly faster than the s.c growing cells (all Log-rank (Mantel–Cox) test).

2.4. The Implantation Site Has No Significant Impact on Take Rate but Does Influence the Infiltration Capacity in Haemopoietic Organs

In the course of the study we compared take rate and infiltration rate depending on the implantation techniques i.t., i.s. and s.c. In general, the take rate was very high in all models with individual models engrafting only in two out of three injection routes (Figure 4a). After i.s. injection three out of 12 models did not engraft. Nevertheless, if the engraftment was successful, the models grew significantly faster as in the i.t. or s.c. setting. One model, LEXF 2713, only engrafted reliably in a disseminated setting. This was vice versa for LEXF 2918 a reliable growth was specifically possible using the s.c. approach. Of note, this was the model with the longest OS time during the establishment phase. Overall, the differences in take rates were not statistically significant (two-way ANOVA). When mice reached termination criteria, BM, PB and spleen were harvested and analysed for infiltration of human leukaemic blasts by FC. Figure 4b plots the percentage of human CD45+ cells over all models in the three different compartments over three implantation techniques. With the exception of one model, LEXF 2431, the s.c. implanted leukaemia cells did not circulate in PB or infiltrated the investigated haemopoietic organs. As a consequence, the infiltration rates of the disseminated approaches in BM, PB and spleen were significantly higher compared to the s.c. method (p < 0.0001, Kruskal–Wallis test). For the two disseminated approaches the infiltration pattern was similar in BM, PB and spleen. Comparing the infiltration rate at the respective injection sites between the three settings also revealed no statistical differences (Kruskal–Wallis test).

2.5. Individual Leukaemia PDX Express a Specific Surface Marker Pattern which is Stable across Different Engraftment Sites

During the development phase of the PDX models, it was already observed that individual models expressed specific surface marker patterns, which were stable across the different passages. In the current study we investigated whether this was also the case when the PDX were implanted via different routes. Figure 5 gives an overview of the percentage of positive cells determined by FC for human CD45, CD34 and CD38 as the most ubiquitously expressed surface markers in our leukaemia PDX panel. The expression pattern was stable and specific across the different compartments and injection routes. A cluster analysis was performed and results were plotted as a heatmap (data not shown) and as a principle component analysis plot (Figure S1). Neither the injection route (i.t vs. i.s. vs. s.c) nor the engraftment site (BM vs. spleen vs. PB) displayed a strong correlation coefficient. Data from individual models showed correlation coefficients > 0.5 (p < 0.05) whereas between different models the data sets did not correlate (Pearson correlation)).

2.6. Cytarabine is Active in Four Leukaemia PDX Models with Varying Degree of Efficacy Depending on the Model and the Implantation Site

Cytarabine as a first line treatment in acute leukaemia was tested in four PDX models implanted i.t. and s.c., respectively. The antitumoral activity was determined by OS and percentage of human CD45+ cells in PB (i.t.) and absolute tumour volume in mm ³ (s.c.). Overall, cytarabine showed statistically significant activity in all four models independent of the injection route. The optimal test/control (T/C) values and their statistical significance of an individual model were similar between the two experimental set-ups (Figure 6). Nevertheless, when comparing the influence of the treatment on OS, only LEXF 2799 showed analogous growth kinetics. Cytarabine induced a 1.6fold increase in OS in the i.t. as well as the s.c. experiment. In LEXF 2665 (3.29fold vs 1.94fold) and LEXF 2531 (3.75fold vs 1.72fold) the prolongation of OS was more pronounced in the s.c. setting. In contrast, animal bearing LEXF 2734 in a disseminated mode benefit more from the treatment with cytarabine as mice bearing the same tumour model s.c. (1.84 fold vs complete remission). In the respective donor patients of the leukaemia PDX models cytarabine achieved a complete remission in the course of the first line treatment. Three of four patients underwent a relapse and received additional therapies during the course of their disease. Thus, despite the minor differences in the duration of the antitumoral effect, the models reflected the clinical response of the donor patient regardless of the implantation site.

3. Discussion

Rodent tumour models have been generated since the 1960s, with xenografts of human tumour models emerging in the 1980s [12]. The development of the nu/nu mutant mouse lacking functional thymus tissue and therefore T cells enabled the engraftment and serial transplantation of human tumour tissue into mice [13]. With the development of more immunosuppressive mouse strains like SCID, NOD/SCID or most recently NOD/SCID/IL2Rγnull (NOG and NSG) mice the engraftment rate of human tumour tissue improved significantly over the last 30 years. Notably the latter enabled not only an improved engraftment of solid tumours but the establishment of PDX from haematological malignancies namely acute leukaemia and different types of lymphomas [14,15].
Recently, there have been multiple efforts in industry as well as academia to establish large panels of well characterized PDX models covering a wide range of different tumour types [16]. Currently, these collections represent a predictive preclinical tool for drug development as well as tumour biology research displaying most consistently the complexity of tumour heterogeneity and molecular diversity of human cancers.
We established a panel of 20 acute leukaemia PDX models recapitulating the diversity of the disease. The models were established by injection into the bone marrow of recipient mice mimicking the bone marrow niche in the patient as close as possible. The i.t. injection was previously described as a successful implantation technique for other haematological malignancies using established cell lines and patient-derived material [17,18,19]. To the best of our knowledge the current study describes for the first time this method in a systematic approach to establish a panel of patient-derived acute leukaemia models.
The comparison with TCGA data revealed a high correlation with the patient cohort indicating that the PDX panel mirrored clinical reality to a great extent. The comparison with published data in this context led to the assumption that besides the engraftment technique and the extended observation period in the establishment phase [20], the choice of the mouse strain is crucial [21,22,23]. Over the course of the establishment phase the overall survival time developed into a model specific feature, which was used in the following experiments as an additional read-out besides the infiltration rate in haemopoietic organs.
To understand the influence of the microenvironment more in detail and to further evaluate different implantation techniques, we compared three implantation routes: the already established i.t. vs. i.s. and s.c. Take rate, infiltration rate, OS and sensitivity towards standard of care cytarabine were compared depending on the injection route. OS was significantly influenced by the engraftment site in 14 out of 20 models. The models in the two groups (OS dependent on injection site vs. non-dependent) displayed no obvious discrimination criteria. Phenotypical and molecular characteristics as well as clinical features from the donor patient were evenly distributed across both groups. Thus, the mechanistic background for this phenomenon remains to be elucidated in future studies. Of note, it has to be taken into account, that due to the nature of the different growth patterns the termination criteria had to be adapted to the tumour growth and might have an impact on differences in OS across different implantation techniques as well. Other groups reported a correlation between patient related characteristics and engraftment latency in the stage of establishment [24]. Meyer et al. [25] systematically investigated the different factors affecting human leukaemia engraftment in mice: the transplantation procedure was one relevant factor they identified as key for the engraftment kinetics. Interestingly, the implantation site had no significant impact on the take rate, although after i.s. implantation three models displayed no engraftment over the course of the experiment. Wang et al. [26] reported the successful engraftment of five out of eleven cases of human all after i.s. injection. This was in line with our observation of nine out of twelve successful engraftments. However, the major difference between the two studies is that Wang et al. [26] compared the engraftment capacity of patient material while our group used PDX models established by i.t. injection.
In our hands, the human leukaemic cells specifically infiltrated the haemopoietic organs when injected i.t. or i.s. Only one s.c. model exhibited a reproducible engraftment of human leukaemic cells in the PB. Of note, no leukaemic cells could be detected in the bone marrow and spleen of those animals. Yan et al. [27] reported the infiltration of haemopoietic organs in animals bearing s.c. implanted leukaemia PDX depending on the growth behaviour of the respective model. The differences in the outcome of those two studies might be mainly related to the fact that different mouse strains were used. The SCID mice used by Yan et al. [27] allowed, due to their higher immunological competence, engraftment of only very aggressive forms of leukaemia. In contrast the NSG used in the current study enabled the engraftment of less aggressive disease forms that do not tend to disseminate in the unfavourable environment of a s.c. implantation.
The leukaemia PDX models were characterized by a definite expression pattern of surface markers determined by FC: CD45, CD34 and CD38 were expressed in a model specific percentage on the tumour cells. The implantation site did not alter the expression of those markers significantly as shown by the correlation and principle component analysis. To the best of our knowledge, similar analyses were not published until now. In a slightly different context, Vick et al. [14] reported that the immunophenotype was stable across different passages, which is completely in line with our observations that the immunophenotype is a model inherent feature.
In four of the leukaemia models we evaluated the antitumoral activity of cytarabine in relation to the application route of the tumour cells. Cytarabine achieved a complete eradication of the tumour cells from peripheral blood and the subcutaneous injection site (if applicable) in all experiments independent of the model or the growth pattern. This reflected very well the patient’s outcome, as all four donor patients underwent complete remission under cytarabine treatment. Tumour growth kinetics after the treatment cycle has ended was impacted by the implantation site in three out of four models. Two models displayed an elongated remission phase in the s.c. setting whereas in one model the opposite was observed. Despite the fact that the database with only one compound in four models was rather small, the results do underline the predictive power of leukaemia PDX models. Our results were in line with studies from other groups investigating the predictivity of PDX for future treatment options in the clinic [28,29,30].
Taken together, all three investigated implantation techniques delivered reliable results and are feasible for preclinical studies in the drug development as well as tumour biology context. Nevertheless, the distinct differences between individual data sets highlight the need for well characterized preclinical models to ensure the meaningful interpretation of data generated using those powerful tools.

4. Materials and Methods

4.1. PDX Establishment

T-and B-cell depleted peripheral blood or bone marrow cells (3 × 10 6 cells/mouse) from leukaemia patients were injected intratibially (i.t.) in the bone marrow of 4–6-week old NSG (NOD/Shi-scid/IL-2Rγnull; Charles River, France) mice. No pre-conditioning with irradiation or chemotherapeutic based treatment was performed prior to injection. Human T-and B-cells were depleted by MACS selection using CD3 and CD19 microbeads (#130-050-101 and #130-050-301, Miltenyi Biotec, Germany) according to the manufacturer’s instructions. Leukaemic cell engraftment was determined by flow cytometry (FC) in bone marrow (BM), peripheral blood (PB) and spleen during the course of the experiment and at the end of the study. Overall survival (OS) served as an additional read-out. The PDX model was defined as established when stable growth over at least four passages and regrowth from liquid nitrogen was observed. All surgery was performed under isoflurane anaesthesia, and all efforts were made to minimize suffering.

4.2. Comparison of Different Injection Routes

Mouse-cell depleted human leukaemic blasts (3 × 10 6 cells/mouse) from established leukaemia PDX were injected into the spleen (i.s.), i.t. or subcutaneous (s.c.) in 4–6-week old NSG (Charles River, France) mice. No pre-conditioning with irradiation or chemotherapeutic based treatment was performed prior to injection. Mouse cell depletion was performed using the mouse cell depletion kit (#130-104-694, Miltenyi Biotec, Bergisch Gladbach, Germany) [31]. The i.s. injection was done via surgical intervention. Under isoflurane anaesthesia and analgesic treatment a small (3 mm) incision was made in the skin and the peritoneum. The leukaemia cells (3 × 106 cells) were injected in a total volume of approx. 30 µL. The puncture canal was closed and haemostasis was carried out using sterile cotton swabs. Skin and peritoneum were closed using resorbable suture material and sore clips. Leukaemic cell engraftment was determined by FC in BM, PB, spleen and subcutaneously growing tumour mass at the end of a study. OS served as an additional read-out. For subcutaneously implanted mice tumour volume was determined twice weekly using electronic calipers and engrafted mice were sacrificed when a tumour volume of 1.800 mm³ was gained. For the disseminated growing tumours, the animals were sacrificed when disease onset was observed by detection of human cells in peripheral blood by FC (> 2%) and deterioration of overall condition indicated by a body weight loss of > 15% over the course of 2 days. All surgery was performed under isoflurane anaesthesia, and all efforts were made to minimize suffering.

4.3. Flow Cytometry Analyses

BM, spleen, tumour or PB cells (5 × 105 to 1 × 106) from tumour bearing animals were harvested and incubated with the primary antibody: hCD45 anti human PB (MHCD4528; Invitrogen, Carlsbad, CA, USA), mCD45.1 anti-mouse PE (553776; BD Biosciences, Franklin Lakes, NJ, USA), hCD34 anti-human FITC (555821; BD Biosciences), hCD38 anti-human BV605 (562665; BD Biosciences), hCD33 anti-human APC (551378; BD Biosciences), hHLA-ABC-PECy5 (555554; BD Biosciences), hCD3-FITC anti-human FITC (555916; BD Biosciences) or isotype control and the mean fluorescence intensity was analysed by flow cytometry. The tumour load was assessed by measurement of the percentages of cells positive for different human surface markers after red blood cell lysis with ammonium chloride. Staining was performed in the presence of CD16/CD32 Abs to block non-specific staining (rat anti-mouse CD16/CD32 (FcγIII/II) receptor IgG2b, BD Pharmingen, San Diego, CA, USA). Samples were analysed on an Attune Acoustic Focusing Cytometer NXT (Applied Biosystems, Waltham, MA, USA) which recorded 50,000 events.

4.4. Standard of Care Testing

Mice were treated with cytarabine (15 mg/kg/day given intravenously; STADApharm, Bad Vilbel, Germany) or vehicle (PBS) for five consecutive days. Treatment commenced upon disease onset: for i.t. implanted animals this was determined as > 2% human CD45+ cells in PB. For s.c. implanted animals the starting median absolute tumour volume was 100 mm³. Beside those read-outs OS served as a read-out for both implantation techniques. The relative volume of an individual tumour on day X (RTVx) was calculated by dividing the absolute volume (mm3) of the respective tumour on day X (Tx) by the absolute volume of the same tumour on the day of randomization, i.e., on day 0 (T0), multiplied by 100, as shown by the following equation: RTVx [%] = Tx/T10 × 100. Tumour inhibition on a particular day (T/Cx) was calculated from the median RTV of a test group and the median RTV of a control group multiplied by 100, as shown by the following equation: T/Cx [%] = median RTVx treated group/median RTVx control group*100. The minimum T/Cx [%] value recorded for a particular group during an experiment represented the maximum anti-tumour activity for the respective compound and were thereby defined as optimal T/C. All experiments were performed according to the relevant animal welfare guidelines published by FELASA and GV-SOLAS in an AAALAC accredited animal facility. Tumour volume and body weight were determined twice per week. Percentage of human CD45+ cells in PB once a week. The calculation of the optimal T/C value was performed as described above for tumour volume.

4.5. Ethics Statement

This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Society of Laboratory Animals (GV SOLAS) in an AAALAC accredited animal facility. All animal experiments were approved by the Committee on the Ethics of Animal Experiments of the regional council (Regierungspräsidium Freiburg, Abt. Landwirtschaft, Ländlicher Raum, Veterinär- und Lebensmittelwesen - Ref. 35, permit-#: G-09/59, G-12/86, G-17/138). The analyses were performed according to the guidelines of the Declaration of Helsinki and good clinical practice. All patients gave their written informed consent for institutional-initiated research studies and analyses of clinical outcome studies conforming to our institutional review board guidelines. All experiments were approved by the ethical commission of the Albert Ludwig University Freiburg (permit-# EK Freiburg: 279/10, 07.09.2010).

4.6. Statistical Analysis

Student’s t test, two-tailed and one-way ANOVA followed by Tukey’s multiple comparisons test were used to calculate all reported p-values. Log-rank [Mantel–Cox]-test was used to calculate overall survival. Descriptive analyses were assessed whenever appropriate and were obtained using GraphPrism software [32]

4.7. Molecular Analysis and Comparison with TCGA Data

The mutational data of the PDX models were retrieved via the online compendium from Charles River [33] and compared to the mutational data from TCGA via the TCGA data portal [34]. The frequency of mutation was compared between both data sets.

4.8. Cluster Analysis

Heat maps and PCA plots were generated using the ClustVis web tool [35]. Centering and unit variance scaling was applied to the sample ID. Both surface markers and sample IDs were clustered using correlation distance and average linkage. Hierarchical clustering of the heatmap started with calculating all pairwise distances. Objects with the smallest distance were merged in each step. The clustering method defined how to go from object level to cluster level when calculating distance between two clusters. The clustering distance was defined as Pearson correlation subtracted from 1.

5. Conclusions

We took advantage of our well characterized leukaemia PDX panel and evaluated the influence of the injection site on different parameters of the PDX growth behaviour. Taken together, the features of an individual PDX model are dictated by the genetic make-up of the model itself and the tumour microenvironment modified in the current study via the engraftment site. Finally, the question of which in vivo model will most precisely reflect the patient situation ultimately may depend on the focus of the scientific question.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/5/1349/s1, Figure S1: Principle component analysis of the expression pattern of three different surface markers in 20 leukaemia PDX models depending on the application route and engraftment site, Table S1: patient characteristics of 20 leukaemia PDX models, Table S2: Mutational landscape of acute leukemia PDX models. Table S3: Statistical analysis of overall survival time of 20 established leukaemia PDX implanted intratibially in NSG mice.

Author Contributions

Conceptualization, J.S., G.G., E.O. and M.L.; data curation, G.G., D.L., M.P., A.E. and E.O.; formal analysis, J.S., A.E. and E.O.; funding acquisition, M.L.; methodology, J.S., G.G., D.L., A.E. and E.O.; project administration, J.S. and M.L.; resources, J.S. and M.L.; supervision, M.L.; validation, J.S., G.G. and E.O.; visualization, J.S. and G.G.; writing–original draft, J.S.; writing–review and editing, G.G., E.O. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Research Foundation Grant agreement CRC 992 MEDEP C04, FOR 2674, LU 429/16-1), JCLS (R14/25) for ML.

Acknowledgments

We thank Anne Löhr and Martina Klotter for excellent technical support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Korn, C.; Méndez-Ferrer, S. Myeloid malignancies and the microenvironment. Blood 2017, 129, 811–822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ferlay, J.; Colombet, M.; Soerjomataram, I.; Mathers, C.; Parkin, D.M.; Piñeros, M.; Znaor, A.; Bray, F. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int. J. Cancer 2019, 144, 1941–1953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Cronin, K.A.; Ries, L.A.; Edwards, B.K. The Surveillance, Epidemiology, and End Results (SEER) Program of the National Cancer Institute. Cancer 2014, 120 (Suppl. 23), 3755–3757. [Google Scholar] [CrossRef]
  5. Deschler, B.; Lübbert, M. Acute myeloid leukemia: Epidemiology and etiology. Cancer 2006, 107, 2099–2107. [Google Scholar] [CrossRef]
  6. Appelbaum, F.R.; Gundacker, H.; Head, D.R.; Slovak, M.L.; Willman, C.L.; Godwin, J.E.; Anderson, J.E.; Petersdorf, S.H. Age and acute myeloid leukemia. Blood 2006, 107, 3481–3485. [Google Scholar] [CrossRef]
  7. Riordan, L. Pharmaprojects; Informa: London, UK, 2020. [Google Scholar]
  8. Bhatia, S.; Daschkey, S.; Lang, F.; Borkhardt, A.; Hauer, J. Mouse models for pre-clinical drug testing in leukemia. Expert Opin. Drug Discov. 2016, 11, 1081–1091. [Google Scholar] [CrossRef]
  9. Keating, M.J. Leukemia: A model for drug development. Clin. Cancer Res. 1997, 3 12 Pt 2, 2598–2604. [Google Scholar]
  10. Ito, M.; Hiramatsu, H.; Kobayashi, K.; Suzue, K.; Kawahata, M.; Hioki, K.; Ueyama, Y.; Koyanagi, Y.; Sugamura, K.; Tsuji, K.; et al. NOD/SCID/gamma(c)(null) mouse: An excellent recipient mouse model for engraftment of human cells. Blood 2002, 100, 3175–3182. [Google Scholar] [CrossRef]
  11. Shultz, L.D.; Lyons, B.L.; Burzenski, L.M.; Gott, B.; Chen, X.; Chaleff, S.; Kotb, M.; Gillies, S.D.; King, M.; Mangada, J.; et al. Human lymphoid and myeloid cell development in NOD/LtSz-scid IL2R gamma null mice engrafted with mobilized human hemopoietic stem cells. J. Immunol. 2005, 174, 6477–6489. [Google Scholar] [CrossRef] [Green Version]
  12. Langdon, S.P.; Hendriks, H.R.; Braakhuis, B.J.; Pratesi, G.; Berger, D.P.; Fodstad, O.; Fiebig, H.H.; Boven, E. Preclinical phase II studies in human tumor xenografts: A European multicenter follow-up study. Ann. Oncol. Off. J. Eur. Soc. Med Oncol. 1994, 5, 415–422. [Google Scholar] [CrossRef] [PubMed]
  13. Pantelouris, E.M. Absence of thymus in a mouse mutant. Nature 1968, 217, 370–371. [Google Scholar] [CrossRef] [PubMed]
  14. Vick, B.; Rothenberg, M.; Sandhofer, N.; Carlet, M.; Finkenzeller, C.; Krupka, C.; Grunert, M.; Trumpp, A.; Corbacioglu, S.; Ebinger, M.; et al. An advanced preclinical mouse model for acute myeloid leukemia using patients‘ cells of various genetic subgroups and in vivo bioluminescence imaging. PLoS ONE 2015, 10, e0120925. [Google Scholar] [CrossRef] [PubMed]
  15. Sugimoto, K.; Hayakawa, F.; Shimada, S.; Morishita, T.; Shimada, K.; Katakai, T.; Tomita, A.; Kiyoi, H.; Naoe, T. Discovery of a drug targeting microenvironmental support for lymphoma cells by screening using patient-derived xenograft cells. Sci. Rep. 2015, 5, 13054. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Conte, N.; Mason, J.C.; Halmagyi, C.; Neuhauser, S.; Mosaku, A.; Yordanova, G.; Chatzipli, A.; Begley, D.A.; Krupke, D.M.; Parkinson, H.; et al. PDX Finder: A portal for patient-derived tumor xenograft model discovery. Nucleic Acids Res. 2019, 47, D1073–D1079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kohart, N.A.; Elshafae, S.M.; Supsahvad, W.; Alasonyalilar-Demirer, A.; Panfil, A.R.; Xiang, J.; Dirksen, W.P.; Veis, D.J.; Green, P.L.; Weilbaecher, K.N.; et al. Mouse model recapitulates the phenotypic heterogeneity of human adult T-cell leukemia/lymphoma in bone. J. Bone Oncol. 2019, 19, 100257. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, M.W.; Kim, H.J.; Yoo, K.H.; Kim, D.S.; Yang, J.M.; Kim, H.R.; Noh, Y.H.; Baek, H.; Kwon, H.; Son, M.H.; et al. Establishment of a bioluminescent imaging-based in vivo leukemia model by intra-bone marrow injection. Int. J. Oncol. 2012, 41, 2047–2056. [Google Scholar] [CrossRef]
  19. Schueler, J.; Wider, D.; Klingner, K.; Siegers, G.M.; May, A.M.; Wäsch, R.; Fiebig, H.H.; Engelhardt, M. Intratibial injection of human multiple myeloma cells in NOD/SCID IL-2Rγ(null) mice mimics human myeloma and serves as a valuable tool for the development of anticancer strategies. PLoS ONE 2013, 8, e79939. [Google Scholar] [CrossRef] [Green Version]
  20. Paczulla, A.M.; Dirnhofer, S.; Konantz, M.; Medinger, M.; Salih, H.R.; Rothfelder, K.; Tsakiris, D.A.; Passweg, J.R.; Lundberg, P.; Lengerke, C. Long-term observation reveals high-frequency engraftment of human acute myeloid leukemia in immunodeficient mice. Haematologica 2017, 102, 854–864. [Google Scholar] [CrossRef] [Green Version]
  21. Feuring-Buske, M.; Gerhard, B.; Cashman, J.; Humphries, R.K.; Eaves, C.J.; Hogge, D.E. Improved engraftment of human acute myeloid leukemia progenitor cells in beta 2-microglobulin-deficient NOD/SCID mice and in NOD/SCID mice transgenic for human growth factors. Leukemia 2003, 17, 760–763. [Google Scholar] [CrossRef] [Green Version]
  22. Rombouts, W.J.; Martens, A.C.; Ploemacher, R.E. Identification of variables determining the engraftment potential of human acute myeloid leukemia in the immunodeficient NOD/SCID human chimera model. Leukemia 2000, 14, 889–897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sanchez, P.V.; Perry, R.L.; Sarry, J.E.; Perl, A.E.; Murphy, K.; Swider, C.R.; Bagg, A.; Choi, J.K.; Biegel, J.A.; Danet-Desnoyers, G.; et al. A robust xenotransplantation model for acute myeloid leukemia. Leukemia 2009, 23, 2109–2117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Culen, M.; Kosarova, Z.; Jeziskova, I.; Folta, A.; Chovancova, J.; Loja, T.; Tom, N.; Bystry, V.; Janeckova, V.; Dvorakova, D.; et al. The influence of mutational status and biological characteristics of acute myeloid leukemia on xenotransplantation outcomes in NOD SCID gamma mice. J. Cancer Res. Clin. Oncol. 2018, 144, 1239–1251. [Google Scholar] [CrossRef] [PubMed]
  25. Meyer, L.H.; Debatin, K.M. Diversity of human leukemia xenograft mouse models: Implications for disease biology. Cancer Res. 2011, 71, 7141–7144. [Google Scholar] [CrossRef] [Green Version]
  26. Wang, N.; Huang, L.; Wang, D.; Wang, J.; Jiang, L.; Zhou, K.; Yang, Y.; Xu, D.; Zhou, J. Successful engraftment of human acute lymphoblastic leukemia cells in NOD/SCID mice via intrasplenic inoculation. Cancer Biol. Ther. 2012, 13, 1158–1164. [Google Scholar] [CrossRef] [Green Version]
  27. Yan, Y.; Wieman, E.A.; Guan, X.; Jakubowski, A.A.; Steinherz, P.G.; O’Reilly, R.J. Autonomous growth potential of leukemia blast cells is associated with poor prognosis in human acute leukemias. J. Hematol. Oncol. 2009, 2, 51. [Google Scholar] [CrossRef] [Green Version]
  28. Townsend, E.C.; Murakami, M.A.; Christodoulou, A.; Christie, A.L.; Köster, J.; DeSouza, T.A.; Morgan, E.A.; Kallgren, S.P.; Liu, H.; Wu, S.C.; et al. The Public Repository of Xenografts Enables Discovery and Randomized Phase II-like Trials in Mice. Cancer Cell 2016, 29, 574–586. [Google Scholar] [CrossRef]
  29. Wang, K.; Sanchez-Martin, M.; Wang, X.; Knapp, K.M.; Koche, R.; Vu, L.; Nahas, M.K.; He, J.; Hadler, M.; Stein, E.M.; et al. Patient-derived xenotransplants can recapitulate the genetic driver landscape of acute leukemias. Leukemia 2017, 31, 151–158. [Google Scholar] [CrossRef] [Green Version]
  30. Wunderlich, M.; Mizukawa, B.; Chou, F.S.; Sexton, C.; Shrestha, M.; Saunthararajah, Y.; Mulloy, J.C. AML cells are differentially sensitive to chemotherapy treatment in a human xenograft model. Blood 2013, 121, e90–e97. [Google Scholar] [CrossRef] [Green Version]
  31. Agorku, D.J.; Tomiuk, S.; Klingner, K.; Wild, S.; Ruberg, S.; Zatrieb, L.; Bosio, A.; Schueler, J.; Hardt, O. Depletion of Mouse Cells from Human Tumor Xenografts Significantly Improves Downstream Analysis of Target Cells. J. Vis. Exp. JoVE 2016, 113, 1–14. [Google Scholar] [CrossRef] [Green Version]
  32. Graphpad.com scientific-software. Available online: https://www.graphpad.com (accessed on 25 May 2020).
  33. Charles River tumor model comendium. Available online: https://compendium.criver.com/ (accessed on 25 May 2020).
  34. The Cancer Genome Atlas Program. Available online: https://portal.gdc.cancer.gov/ (accessed on 25 May 2020).
  35. Metsalu, T.; Vilo, J. ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015, 43, W566–W570. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Establishment of a patient-derived xenograft panel of 20 acute leukaemia models. (a) Implantation scheme over several passages. Patient blasts were implanted into the bone marrow of NSG mice. When disease onset was observed by detection of human cells in peripheral blood and deterioration of overall condition, mice were sacrificed and cells from bone marrow and spleen harvested. After a mouse cell depletion human leukaemic cells were re-implanted into the bone marrow of recipient mice. This procedure was repeated until stable growth was observed. Subsequently, the PDX leukaemic cells were implanted into the bone marrow, in the spleen or subcutaneously into NSG mice. The monitoring of disease onset was performed as described before. Additionally, the antitumoral activity of cytarabine was evaluated in all three settings in four different PDX models. (b) For 16 models, mutational data were available. The comparison with TCGA data revealed a highly significant correlation between the two data sets. Mutational data were either acquired based on WES or RNAseq.
Figure 1. Establishment of a patient-derived xenograft panel of 20 acute leukaemia models. (a) Implantation scheme over several passages. Patient blasts were implanted into the bone marrow of NSG mice. When disease onset was observed by detection of human cells in peripheral blood and deterioration of overall condition, mice were sacrificed and cells from bone marrow and spleen harvested. After a mouse cell depletion human leukaemic cells were re-implanted into the bone marrow of recipient mice. This procedure was repeated until stable growth was observed. Subsequently, the PDX leukaemic cells were implanted into the bone marrow, in the spleen or subcutaneously into NSG mice. The monitoring of disease onset was performed as described before. Additionally, the antitumoral activity of cytarabine was evaluated in all three settings in four different PDX models. (b) For 16 models, mutational data were available. The comparison with TCGA data revealed a highly significant correlation between the two data sets. Mutational data were either acquired based on WES or RNAseq.
Cancers 12 01349 g001
Figure 2. Overall survival in 20 leukaemia PDX models over different passages. Leukaemic blasts from donor patients were injected intratibially and mice monitored as described in detail in Figure 1 over the course of different passages. Until passage four the overall survival time changed between passages and stabilized at a model specific value afterwards.
Figure 2. Overall survival in 20 leukaemia PDX models over different passages. Leukaemic blasts from donor patients were injected intratibially and mice monitored as described in detail in Figure 1 over the course of different passages. Until passage four the overall survival time changed between passages and stabilized at a model specific value afterwards.
Cancers 12 01349 g002
Figure 3. (a) Overall survival of 12 different leukaemia PDX dependent on the implantation route. In the case of the intrasplenic injection the outcome was influenced by the fact that three out of twelve tested models did not show any tumour growth during the observation period. In the left panel the OS is plotted for all models. The differences were statistically significant for i.s. vs. s.c.. In case only the models were taken into account which did show successful engraftment of leukaemic cells the i.s. implanted cells grew significantly faster than the s.c growing cells (right panel, all Log-rank (Mantel–Cox) test). (b) Overall survival determined in 20 leukaemia PDX models for up to three different injection sites. For 20 PDX models the overall survival rate was determined over time depending on the application route of the cancer cells. For 14 out of the 20 models a statistically significant difference between the implantation techniques was determined (Log-rank (Mantel–Cox) test).
Figure 3. (a) Overall survival of 12 different leukaemia PDX dependent on the implantation route. In the case of the intrasplenic injection the outcome was influenced by the fact that three out of twelve tested models did not show any tumour growth during the observation period. In the left panel the OS is plotted for all models. The differences were statistically significant for i.s. vs. s.c.. In case only the models were taken into account which did show successful engraftment of leukaemic cells the i.s. implanted cells grew significantly faster than the s.c growing cells (right panel, all Log-rank (Mantel–Cox) test). (b) Overall survival determined in 20 leukaemia PDX models for up to three different injection sites. For 20 PDX models the overall survival rate was determined over time depending on the application route of the cancer cells. For 14 out of the 20 models a statistically significant difference between the implantation techniques was determined (Log-rank (Mantel–Cox) test).
Cancers 12 01349 g003
Figure 4. Influence of the injection site on take rate (a) and infiltration rate in different compartments (b) in 20 leukaemia PDX models. In general, the take rate was not influenced by the implantation site. Nevertheless, a few models did only engraft in two out of the three injection sites. A no take was defined by no growth after three individual experiments with n > 3 mice and an observation time >2-fold of the passaging time of the individual model. **** p < 0.0001.
Figure 4. Influence of the injection site on take rate (a) and infiltration rate in different compartments (b) in 20 leukaemia PDX models. In general, the take rate was not influenced by the implantation site. Nevertheless, a few models did only engraft in two out of the three injection sites. A no take was defined by no growth after three individual experiments with n > 3 mice and an observation time >2-fold of the passaging time of the individual model. **** p < 0.0001.
Cancers 12 01349 g004
Figure 5. Infiltration pattern of 20 leukaemia PDX in bone marrow, spleen and peripheral blood depending on the injection site of the leukaemic cells. Each PDX model displayed its specific expression pattern of the three investigated surface markers CD45, CD38 and CD34 in all tested settings. Despite for one model, LEXF 2431, leukaemic cells did not infiltrate the haemopoietic organs when injected subcutaneously. A no take was defined by no growth after three individual experiments with n > 3 mice and an observation time > 2-fold of the passaging time of the individual model. n.d., not determined.
Figure 5. Infiltration pattern of 20 leukaemia PDX in bone marrow, spleen and peripheral blood depending on the injection site of the leukaemic cells. Each PDX model displayed its specific expression pattern of the three investigated surface markers CD45, CD38 and CD34 in all tested settings. Despite for one model, LEXF 2431, leukaemic cells did not infiltrate the haemopoietic organs when injected subcutaneously. A no take was defined by no growth after three individual experiments with n > 3 mice and an observation time > 2-fold of the passaging time of the individual model. n.d., not determined.
Cancers 12 01349 g005
Figure 6. Determination of sensitivity towards cytarabine in four leukaemia PDX models depending on the implantation route of the tumour cells. Antitumor activity was determined by percentage of leukaemic blasts in peripheral blood (i.t. injected animals), tumour volume (s.c. implanted animals) and overall survival (i.t. and s.c.). Cytarabine was active an all tested models and settings although with varying degree depending on the model and the implantation site. Significance was calculated using the Log-rank (Mantel–Cox) test for OS and Mann–Whitney for tumour volume and percentage of human CD45+ cells.
Figure 6. Determination of sensitivity towards cytarabine in four leukaemia PDX models depending on the implantation route of the tumour cells. Antitumor activity was determined by percentage of leukaemic blasts in peripheral blood (i.t. injected animals), tumour volume (s.c. implanted animals) and overall survival (i.t. and s.c.). Cytarabine was active an all tested models and settings although with varying degree depending on the model and the implantation site. Significance was calculated using the Log-rank (Mantel–Cox) test for OS and Mann–Whitney for tumour volume and percentage of human CD45+ cells.
Cancers 12 01349 g006
Table 1. Characteristics of the leukaemia patient-derived xenografts (PDX) panel.
Table 1. Characteristics of the leukaemia patient-derived xenografts (PDX) panel.
Nr of Modelspts DataPDX Data
Age, Median (Range)Gender (f:m)Pre-TreatmentSubtypeMutational StatusSoC Data AvailableWES Data AvailableRNAseq Data Available
2063
(43–80) years
11:9none13AML18FLT3-ITD6191410
PBSCT3ALL1FLT3-TKD1
other4APL1NPM1-A6
Abbreviations: SoC, standard of care; WES, whole exome sequencing; RNAseq, RNA sequencing.
Table 2. Comparison of overall survival (OS) times within one model depending on the implantation route (Log-rank (Mantel–Cox) test).
Table 2. Comparison of overall survival (OS) times within one model depending on the implantation route (Log-rank (Mantel–Cox) test).
Model NameIntrasplenic vs. IntratibialIntrasplenic vs. SubcutaneousIntratibial vs. Subcutaneous
LEXF 2412n.s.n.s.n.s.
LEXF 2431<0.0001<0.00010.0091
LEXF 2531n.d.n.d.n.s.
LEXF 2665<0.0001n.s.0.0002
LEXF 2713<0.0001<0.0001<0.0001
LEXF 27340.0334n.s.0.0208
LEXF 2799n.s.n.s.n.s.
LEXF 2824<0.0001<0.00010.0001
LEXF 2848n.s.n.s.n.s.
LEXF 28970.0002<0.0001n.s.
LEXF 2918n.s.<0.0001<0.0001
LEXF 2943n.s.n.s.n.s.
LEXF 29570.0009n.s.<0.0001
LEXF 2964n.d.n.d.0.001
LEXF 2966n.d.n.d.0.0003
LEXF 2997n.d.n.d.<0.0001
LEXF 4010n.d.n.d.n.s.
LEXF 4052n.d.n.d.0.0006
LEXF 4096n.d.n.d.<0.0001
LEXF 4128n.d.n.d.<0.0001
all modelsn.s.0.049n.s
all models w/o no takes<0.0001<0.0001n.s.
Abbreviations: n.d: not determined; n.s: not significant.

Share and Cite

MDPI and ACS Style

Schueler, J.; Greve, G.; Lenhard, D.; Pantic, M.; Edinger, A.; Oswald, E.; Lübbert, M. Impact of the Injection Site on Growth Characteristics, Phenotype and Sensitivity towards Cytarabine of Twenty Acute Leukaemia Patient-Derived Xenograft Models. Cancers 2020, 12, 1349. https://doi.org/10.3390/cancers12051349

AMA Style

Schueler J, Greve G, Lenhard D, Pantic M, Edinger A, Oswald E, Lübbert M. Impact of the Injection Site on Growth Characteristics, Phenotype and Sensitivity towards Cytarabine of Twenty Acute Leukaemia Patient-Derived Xenograft Models. Cancers. 2020; 12(5):1349. https://doi.org/10.3390/cancers12051349

Chicago/Turabian Style

Schueler, Julia, Gabriele Greve, Dorothée Lenhard, Milena Pantic, Anna Edinger, Eva Oswald, and Michael Lübbert. 2020. "Impact of the Injection Site on Growth Characteristics, Phenotype and Sensitivity towards Cytarabine of Twenty Acute Leukaemia Patient-Derived Xenograft Models" Cancers 12, no. 5: 1349. https://doi.org/10.3390/cancers12051349

APA Style

Schueler, J., Greve, G., Lenhard, D., Pantic, M., Edinger, A., Oswald, E., & Lübbert, M. (2020). Impact of the Injection Site on Growth Characteristics, Phenotype and Sensitivity towards Cytarabine of Twenty Acute Leukaemia Patient-Derived Xenograft Models. Cancers, 12(5), 1349. https://doi.org/10.3390/cancers12051349

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop