Prognostic value of patient-derived xenograft engraftment in pediatric sarcomas.

The goals of this work were to identify factors favoring patient-derived xenograft (PDX) engraftment and study the association between PDX engraftment and prognosis in pediatric patients with Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma. We used immunodeficient mice to establish 30 subcutaneous PDX from patient tumor biopsies, with a successful engraftment rate of 44%. Age greater than 12 years and relapsed disease were patient factors associated with higher engraftment rate. Tumor type and biopsy location did not associate with engraftment. PDX models retained histology markers and most chromosomal aberrations of patient samples during successive passages in mice. Model treatment with irinotecan resulted in significant activity in 20 of the PDXs and replicated the response of rhabdomyosarcoma patients. Successive generations of PDXs responded similarly to irinotecan, demonstrating functional stability of these models. Importantly, out of 68 tumor samples from 51 patients with a median follow-up of 21.2 months, PDX engraftment from newly diagnosed patients was a prognostic factor significantly associated with poor outcome (p = 0.040). This association was not significant for relapsed patients. In the subgroup of patients with newly diagnosed Ewing sarcoma classified as standard risk, we found higher risk of relapse or refractory disease associated with those samples that produced stable PDX models (p = 0.0357). Overall, our study shows that PDX engraftment predicts worse outcome in newly diagnosed pediatric sarcoma patients.

Introduction

Pediatric sarcomas are a rare and heterogeneous group of skeletal and soft tissue malignancies accounting for approximately 12% of all childhood solid tumors [1]. The most frequently occurring are Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma, altogether representing 13% of all malignant tumors in patients younger than 14 years, and up to 18% of all tumors among patients between 15 and 19 years [2, 3, 4]. With up‐to‐date treatment modalities, the 5‐year survival rates have improved in the last decade to 79% for Ewing sarcoma, 73% for osteosarcoma, and 70% for rhabdomyosarcoma [5, 6, 7]. Unfortunately, intensive and multimodal treatments currently used in pediatric sarcomas have topped their efficacy [1, 8, 9]. A substantial proportion of patients relapse and do not respond consistently to rescue chemotherapies because of tumor evolution into a chemoresistant phenotype [10].

Patient‐derived xenograft (PDX) tumor models in immunodeficient mice are well‐recognized tools for the study of rare tumors and their treatment [11, 12, 13]. In colorectal cancer, epithelial ovarian cancer, and lung adenocarcinoma, for instance, PDXs retain the main properties of the original patient tumor, including histology, genetic and genomic alterations, gene expression profile, and heterogeneity [6, 11, 14, 15]. Successful engraftment of a patient tumor biopsy in animals might predict the outcome of patients [16]. For breast, pancreatic, and lung cancers, engraftment in mice is associated with worse patient prognosis [15, 17, 18]. Likely, there is a selection process enabling the most aggressive clones to expand and prevail in successive passages in mice [11, 13, 19, 20, 21]. PDX models of most pediatric sarcomas are suitable for such establishment [22, 23, 24]. Whether the parameter ‘PDX engraftment’ could be useful for identifying pediatric sarcoma patients at risk of relapse or refractory disease is not well characterized. To address this question, here we establish PDX models from Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma patients and include available clinical data along with the PDX data [25]. To characterize this new preclinical platform, we address whether PDX tumors retained the main histologic, genomic, and functional properties of patient tumors during successive passages in mice. Sample characterization includes histopathology markers, chromosomal profiles (copy number alteration [CNA]), and preclinical treatment assays (efficacy of irinotecan).

Materials and methods

Xenotransplantation of patient samples

The appropriate Institutional Review Board approved the research protocol associated with this project (M‐1608‐C). Patients or their legal guardians signed informed consent at Sant Joan de Déu Hospital (SJD, Barcelona, Spain). Clinical information parameters obtained from patients are listed in supplementary material, Table S1. Patients were biopsied using either a surgical procedure (open surgery or minimally invasive surgery) or percutaneous tissue core biopsy with a 18G Tru‐Cut needle (Argon Medical Devices, Athens, TX, USA). Biopsies were used fresh or were cryopreserved. For the cryopreservation process, tissues were minced as 2 × 2 × 1 mm pieces with a scalpel on sterile petri dishes. Three to four of these pieces were added to 2 ml cryogenic vials (Corning, Glendale, AZ, USA) containing 1.5 ml of inactivated fetal bovine serum with 10% dimethyl sulfoxide (Life Technologies, Grand Island, NY, USA). Vials were stored overnight in a freezing container (Cool Cell; Corning) at −80 °C and then transferred to liquid nitrogen for long‐term storage.

Work with mice adhered to the European regulations and was approved by the animal experimental ethics committee at the Universitat de Barcelona (animal protocol numbers 135/11 and 134/18). In brief, we implanted patient tumor biopsies (freshly excised or cryopreserved) subcutaneously in 4–6‐week‐old female NOD‐SCID mice (Envigo, Barcelona, Spain) as previously described [26]. Tumors engrafting successfully were named F0 generation. When tumor volume achieved 1,000–1,500 mm3, we performed passage of freshly excised tissue to athymic nude mice (Envigo). We used NOD‐SCID mice for the passage of specific rhabdomyosarcoma models. We named subsequent generations with the number of the passage (e.g. F1 was the first filial generation after initial F0 engraftment). Information parameters obtained from successfully engrafted PDX models are listed in supplementary material, Table S2.

To interrogate whether clinical parameters of the patient disease, including the sample collection method (surgical biopsy versus Tru‐Cut needle), correlated with the likelihood of successful PDX engraftment, we used mixed‐effects logistic regression analysis and calculated the odds ratios and their corresponding 95% confidence intervals (CIs). To dichotomize age, we computed an optimum threshold by maximizing the Youden index using negative or positive engraftment [27].

Analysis of fusion genes

We amplified Ewing sarcoma fusion gene types EWSR1‐ERG and EWSR1‐FLI1, and rhabdomyosarcoma fusion genes PAX3‐FKHR and PAX7‐FKHR with the polymerase chain reaction (PCR). The sequences of the primers were EWSR1 forward, 5′‐TCCTACAGCCAAGCTCCAAGTC‐3′; FLI1 reverse, 5′‐GTGTCAGGCATGGAGGATGGA‐3′; ERG reverse, 5′‐GAGAAGGCATATGGCTGGTGG‐3′; PAX3 forward, 5′‐AGGCATGGATTTTCCAGCTATA‐3′; PAX7 forward, 5′‐TCTGCCTACGGAGCCCG‐3′; and FKHR reverse, 5′‐GGGACAGATTATGACGAATTGAATT‐3′. Ewing sarcoma fusion genes were detected in agarose gels, while rhabdomyosarcoma fusion genes were quantified with a real‐time quantitative PCR technique, using the probe FAM‐5′‐CCGGTCAGCAACGGCCTGTCT‐3′.

Histopathology

During the process of PDX establishment, the original patient tumors and PDX models at different passages must be characterized and compared [16]. First, we compared histology on hematoxylin and eosin staining, and the expression of tumor‐related proteins in the patient tumor biopsies and PDX tumors. To analyze whether protein expression (detected by immunostaining) changed upon successive mouse‐to‐mouse transplantation, we studied PDX tissues at three different filial generations: F0, as the first generation in mice, and F2 and F5, as early and late successive passages in mice. Tissues were fixed overnight in buffered formalin and embedded in paraffin. To stain human tumor cells in engrafted tumor tissues, we used anti‐human nuclei primary antibody (1:200, MAB4383; Merck Millipore, Burlington, MA, USA). To detect proliferating cells, we used anti‐Ki67 antibody (1:200, ACK02; Leica Biosystems, Wetzlar, Germany). We stained proteins highly expressed by specific tumor types, such as CD99 for Ewing sarcoma (1:20, NCL‐L‐CD99‐187; Leica Biosystems), MyoD1 for rhabdomyosarcoma (1:20, M3512; Agilent Technologies, Santa Clara, CA, USA), and the secreted protein acidic and rich in cysteine (SPARC) for osteosarcoma (1:80, 35‐5500; Invitrogen, Carlsbad, CA, USA).

PDX growth

We measured tumor growth rate as the time needed to achieve the experimental endpoint, i.e. tumor volume of 1,500 mm3. Tumor volume was calculated as (length × width2)/2, length being the longitudinal diameter and width the transverse diameter of the subcutaneous tumor.

Analysis of CNAs

We compared the CNA profiles of tumor biopsies and PDX tissue at generations F2 and F5. We performed whole‐genome analysis using the high‐density array CytoScan® HD platform (Affymetrix, Thermo Fisher Scientific, Waltham, MA, USA) as previously described [10]. Total DNA from frozen samples was digested, ligated, PCR‐amplified and purified, fragmented, biotin‐labeled, and hybridized according to manufacturer's instructions. We analyzed the CytoScan® HD Array data with the Chromosome Analysis Suite (ChAS) software (Affymetrix, Thermo Fisher Scientific).

Irinotecan activity

We used 20 of the new PDX models to evaluate the activity of irinotecan, as a model drug with established activity against pediatric solid tumors including sarcomas [10, 28, 29]. We reasoned that if the PDXs were truly representative of the patients' diseases, they would likely respond to this drug and such response should remain stable over successive PDX generations. Among the patients from which these PDX models were established, 11 Ewing sarcoma PDXs resulted from eight patients who received irinotecan as part of their rescue treatment (combined with temozolomide, with or without additional vincristine or trabectedin); 7 rhabdomyosarcoma PDXs from six patients who received irinotecan as part of the upfront treatment of the primary tumor (either alone or combined with carboplatin or vincristine); and 2 osteosarcoma PDXs from two patients who did not receive irinotecan. To compare responses to treatment of patients and PDXs, we selected newly diagnosed patients treated with irinotecan as upfront treatment and with available PDX models from their primary tumor. We evaluated treatment response in such patients using the responsive evaluation criteria in solid tumors (RECIST) protocol or the metabolic tumor volume.

For all efficacy studies in mice, we inserted freshly excised PDX tumors (obtained from one mouse of the immediate earlier generation) in both flanks of athymic nude or NOD‐SCID mice. Upon engraftment (tumor volumes ranging 100–500 mm3), mice were distributed to control or treatment groups, with care that tumor volume means and standard deviations (STDEV) in both groups were similar. Treatment groups received one cycle of 10 mg/kg/day irinotecan (Hospira, Lake City, IL, USA) in a 5‐day‐on‐2‐off regimen, intraperitoneal, for two consecutive weeks, as previously described [30]. Control groups received saline using the same regimen as treated groups. Tumor volume was measured three times a week, until day 14, in which response to treatment was evaluated as previously described [10]. We defined complete response (CR) as tumor mass <50 mm3 and >50% reduction at the end of treatment (day 14); partial response (PR) as tumor volume regression ≥50% at day 14 but tumor volume ≥50 mm3; stable disease (SD) as <50% regression and ≤25% increase in initial volume at day 14; and progressive disease as <50% regression from initial volume and >25% increase in initial volume at day 14.

To address whether the response to the drug changed upon consecutive passages of the PDX in mice, we evaluated the response to irinotecan of different filial generations in three of the models, one early after initial engraftment (F ≤ 2) and the second one later (F ≥ 6). After one cycle of irinotecan, we followed up all animals weekly until tumor regrowth to endpoint (1,500 mm3) or day 100 to build Kaplan–Meier curves and estimate median survival times. We used the log‐rank test to compare survival curves between treatments and between different generations from the same PDX model.

Prognostic value of positive tumor engraftment

The main purpose of this work was to address whether successful engraftment in mice of tumor biopsies of newly diagnosed or relapsed Ewing sarcoma, rhabdomyosarcoma, and osteosarcoma patients had a prognostic effect on overall survival (OS) and event‐free survival (EFS) of these patients.

To assess the prognostic effect of successful PDX engraftment on EFS and OS of the biopsied children, we used Cox models with a robust sandwich variance estimator to account for the correlation of the data. We calculated hazard ratios (HR) and their corresponding 95% CI. We determined median EFS and OS with the Kaplan–Meier method. For each patient, we calculated EFS as the time from biopsy until date of patient recurrence, death, or last follow‐up, and OS as the time from biopsy until patient's death or last follow‐up.

To address whether PDX engraftment of primary tumors predicted prognosis of patients classified as ‘standard risk’ according to clinical methods, Ewing sarcoma patients whose biopsies were obtained at initial diagnosis (diagnostic cohort) were divided into two groups, standard risk and high risk, according to their risk stratification following the clinical trial GEIS21 [5]. For each of such patient groups, we studied the association between the variables PDX engraftment and relapse, using the Fisher's exact test.

Statistics

We performed statistical analyses using GraphPad Prism 8 software (Graphpad, La Jolla, CA, USA) and R software version 3.5.2 (R Foundation for Statistical Computing, Vienna, Austria).

Results

Patients

We included 51 children, 37 males, and 14 females in this study, with newly diagnosed or relapsed Ewing sarcoma (N = 31), osteosarcoma (N = 10), or rhabdomyosarcoma (N = 10, five each embryonal and alveolar), from November 2010 to November 2019. Median age at inclusion was 11.4 years (range: 0.52–17.9). Most patients (73%) were nonmetastatic at diagnosis. Details on tumor types, fusion genes, anatomic location of primary tumors, risk stratification, and presence of metastases are in supplementary material, Tables S3 (Ewing sarcoma), S4 (osteosarcoma), and S5 (rhabdomyosarcoma).

Biopsies and PDXs

We obtained 68 biopsy samples. Engraftment outcomes of individual patient samples are in supplementary material, Table S6. Forty‐one tumor samples were Ewing sarcoma, 12 osteosarcoma, and 15 rhabdomyosarcoma. In total, 44% of samples engrafted efficiently in mice, with engraftment rates for each tumor type shown in Table 1. Age older than 12 years and relapse were factors associated with increased engraftment rate when pooling all diseases together (Table 1). Tumor type, biopsy location, and presence of metastases at diagnosis were not associated with engraftment rate. Procedural factors such as sample conservation (freshly excised or cryopreserved) and application of Matrigel (Corning, Glendale, AZ, USA) to the sample inserted in the mice did not affect engraftment rate. Biopsies obtained by surgery or by the Tru‐Cut method engrafted with a similar success rate (see supplementary material, Table S7).

Table 1

Association of patient factors with engraftment.

Factor	No. of Samples	Engrafted (%)	No. of Engrafted	Odds Ratio (95% CI)	P
Age (years)
<12	36	30.6	11
≥12	32	59.4	19	3.88 (1.1–14.3)	0.042
Tumor type
All	68	44.1	30
Ewing sarcoma	41	41.5	17
Osteosarcoma	12	41.7	5	0.99 (0.2–5.6) *
Rhabdomyosarcoma	15	53.3	8	1.97 (0.4–10.7) *	0.69
Alveolar	9	66.7	6
Embryonal	6	33.3	2
Timing of surgery
Diagnosis	31	22.6	7
Relapse	37	62.2	23	15.4 (1.3–182.8)	0.031
Biopsy origin
Limbs	22	31.8	7
Head and neck	14	64.3	9	21.1 (0.04–11 056.5) †
Chest wall and ribs	7	42.9	3	3.4 (0.05–215.9) †
Lung or pleura	13	61.5	8	15 (0.07–3,075.8) †
Pelvic bones	3	33.3	1	1.2 (0.01–134.4) †	0.17
Muscle	2	100	2	‡
Testes	1	0	0	‡
Vertebral spine	6	0	0	‡
Metastasis at diagnosis
No	48	37.5	18
Yes	20	60.0	12	3.62 (0.64–20.50)	0.15

Compared to Ewing sarcoma.

Compared to limbs.

Excluded from statistical analysis of biopsy origin due to lack of an event in the number of engrafted or nonengrafted samples.

Of the 51 patients, 15 had two or more biopsies included in this study. Nine of them produced two stable PDXs. We studied whether the engraftment of a first biopsy from one patient predicted engraftment of successive biopsies of the same patient at more advanced stages of the disease. We found that 9 out of 11 such second biopsies engrafted after a first biopsy with positive PDX engraftment, but the association between first and second engraftments was not statistically significant, likely due to the low number of cases in the study (see supplementary material, Table S8).

Because Ewing sarcoma patients were predominant in the study cohort, we also analyzed their data separately. We found that age older than 12 years was not determinant in this separate analysis, while relapse and the presence of metastases at diagnosis were significantly associated with engraftment (see supplementary material, Table S9).

Comparative histopathology of patient tumors and PDX samples

Histopathology of the tumors and expression of tumor markers CD99, SPARC, and MyoD1 did not change significantly upon successive mouse‐to‐mouse passaging in 12 PDX models studied by IHC (Figure 1). The typical small‐blue‐round‐cell tumor and histologic architecture of the original biopsies remained preserved in Ewing sarcoma xenografts at early and late passages. Stroma around Ewing sarcoma and rhabdomyosarcoma tumor cells, and osteoid in osteosarcoma tumors, did not change significantly upon engraftment. Anti‐human nuclear antigen staining was positive in all samples (see supplementary material, Figure S1). Tumor anatomy and histology markers were similar in PDX pairs established from the same patients (see supplementary material, Figure S2).

Figure 1

Comparative histology (hematoxylin and eosin and IHC staining) of six representative cases of original human tumor biopsies and the corresponding PDXs at early passages (F0/F2) and late passage (F5). CD99 (cell membrane), SPARC (cytoplasm), and MyoD1 (nuclear) are stained in brown. These representative samples were selected from six Ewing sarcomas, three osteosarcomas, and three rhabdomyosarcomas with complete histopathology studies. All images were obtained using a microscope at ×40 objective magnification. Scale bar represents 50 μm.

Tumor growth rate

The number of proliferating cells (% of Ki67‐positive nuclei by immunohistochemistry (IHC); representative images in supplementary material, Figure S3A and quantification in supplementary material, Figure S3B) in patient biopsies was 54 ± 26% (mean and STDEV of N = 12 cases). Their corresponding PDX samples at passage F0 showed a lower number of stained cells (37 ± 19%). PDX tumors at generations higher than F0 showed increasing number of Ki67‐positive cells (40 ± 23% at F2 and 44 ± 21% at F5; see supplementary material, Figure S3B). Although the differences in Ki67 staining were not statistically significant, tumors at F5 achieved endpoint significantly faster (53 ± 21 days; mean ± STDEV of N = 20 cases) than at the earlier stage F2 (92 ± 47 days; see supplementary material, Figure S3C).

CNA in patient biopsies and corresponding PDXs

We selected one Ewing sarcoma tumor with low CNAs and one rhabdomyosarcoma with highly aberrant chromosomic profile for CNA analyses. Karyotype alterations of Ewing sarcoma (HSJD‐ES‐012) were equivalent for patient tumor and PDX models at F2 and F5 [10]. Ewing sarcoma CNA included a copy‐neutral loss of heterogeneity in chromosomes 3p and 20q as previously reported [10] (see supplementary material, Figure S4). For rhabdomyosarcoma (HSJD‐RMS‐4), both PDX models, F2 and F5, shared 90% of the alterations observed in the patient tumor. Chromosomal aberrations shared by the biopsy and the PDX included whole gains in one or two copies of chromosomes 2, 3, 4, 6, 10, 11, 12, 14, 16, 17, 18, 19, 20, 21, and X; partial gain in chromosomes 1, 7, 9, and 13; and whole chromosome loss in chromosome 15. Nonequivalent alterations among biopsy and PDX included gain of two copies of chromosomes 5, 8, and 22 in the tumor biopsy, and partial gains of chromosomes 5 and 22 and one copy gain of chromosome 8 in both PDX models (see supplementary material, Figure S4).

Antitumor activity of irinotecan in subcutaneous PDXs

All PDX models (see identification and clinical details in supplementary material, Table S10) had a measurable response after one cycle of irinotecan (Figure 2A and supplementary material, Table S11). We found that 9 of 20 models (45%) achieved CR to treatment and all achieved at least SD (Figure 2B).

Figure 2

Antitumor activity of irinotecan in subcutaneous PDX. (A) Change in tumor volume (mean and STDEV of 3–15 tumors) at the end of irinotecan treatment (day 14). (B) Percentage of tumor models achieving each response. (C) Tumor volume (% of volume at treatment start) in three PDX pairs at passage F ≤ 2 or F ≥ 6, treated with one cycle of irinotecan (treatment) or not treated (control). Models and F were HSJD‐ES‐009 (F2 versus F6), HSJD‐ES‐017 (F1 versus F8) and HSJD‐aRMS‐2 (F2 versus F10).

We compared responses to treatment of patient and PDX of four newly diagnosed rhabdomyosarcoma patients receiving upfront irinotecan, three combined with carboplatin, and one with vincristine (see supplementary material, Table S12). The three patients treated with irinotecan plus carboplatin achieved the same response to treatment as their corresponding PDX models (HSJD‐aRMS‐7, PR; HSJD‐aRMS‐10, PR; and HSJD‐RMS‐11, CR) to irinotecan. The patient treated with irinotecan plus vincristine achieved a lower degree of response (SD) compared to the corresponding PDX (HSJD‐RMS‐9, CR).

In the PDX models for which irinotecan activity was evaluated at two different filial generations, both tumors at passages F ≥ 6 and F ≤ 2 showed similar responses to irinotecan (Figure 2C). Upon treatment cessation, the PDXs grew similarly, even those that achieved a CR, and achieved similar median survivals until endpoint (see supplementary material, Table S13).

Prognostic value of positive PDX engraftment in mice

Patients included in the analysis had a median follow‐up of 21.2 months (range: 2.40–101 months). Figure 3 shows the outcome of the study cohort according to engraftment in mice. Median EFS for patients with positive engraftment (7.56 months) was significantly shorter than that of patients with negative engraftment (22.4 months; p = 0.00021; HR: 3.34; 95% CI: 1.77–6.31; Figure 3A). Similarly, median OS was shorter for patients with positive engraftment (19.9 months) compared to that of patients with negative engraftment (42.5 months; p = 0.0069; HR: 2.51; 95% CI: 1.29–4.89; Figure 3B).

Figure 3

Kaplan–Meier estimation of (A) EFS and (B) OS among all patients with positive (n = 30) or negative (n = 38) engraftment of their PDX.

Because relapse was associated with increased engraftment rate, we stratified patients into those whose samples were obtained at initial diagnosis (diagnostic cohort) and at relapse (relapse cohort). We further separated each cohort according to positive and negative engraftment. In the diagnostic cohort, positive engraftment associated with shorter median EFS (20.6 months). In contrast, median EFS was not reached for negative engraftment (p = 0.040; HR: 2.87; 95% CI: 1.05–7.85; Figure 4A). In this cohort, patients with positive engraftment also had shorter median OS (34.2 months), although not statistically significant compared to patients with negative engraftment (p = 0.19; HR: 2.06; 95% CI: 0.70–6.0; Figure 4B).

Figure 4

Kaplan–Meier estimation of (A) EFS and (B) OS among patients of the diagnostic cohort with positive (n = 7) or negative (n = 24) engraftment of their PDX.

In the relapse cohort, we did not find a significant association of positive and negative engraftment with median EFS of 4.56 and 12.2 months, respectively (p = 0.17; HR: 1.70; 95% CI: 0.80–3.63; Figure 5A) and OS of 16.7 and 25.7 months, respectively (p = 0.30; HR: 1.51; 95% CI: 0.69–3.32; Figure 5B).

Figure 5

Kaplan–Meier estimation of (A) EFS and (B) OS among patients of the relapse cohort with positive (n = 23) or negative (n = 14) engraftment of their PDX.

We studied 18 newly diagnosed Ewing sarcoma patients, 7 of which were standard risk and 11 high risk, with follow‐up times of 15.5–101 months. Two out of the seven patients classified as standard risk relapsed and eventually died of disease. Biopsies from these two patients successfully established PDXs in mice, while biopsies from the remaining five standard risk patients did not. In fact, we found that positive PDX engraftment predicted relapse of patients classified as standard risk (p = 0.0357; Fisher's exact test). Among the 11 newly diagnosed Ewing sarcoma patients classified as high risk, two established PDX models. Both patients relapsed or were refractory to treatment and eventually died of disease. Among the remaining nine high‐risk patients, for whom PDX engraftment was unsuccessful, another two relapsed. The association of positive PDX engraftment and relapse was not significant in these high‐risk patients (p = 0.1091).

Discussion

PDXs are powerful tools for cancer biology and precision cancer medicine studies, and may be useful to predict patient prognosis [17]. Determination of prognosis is especially relevant in pediatric patients, 80% of whom will survive cancer and reach adulthood [31]. Identification of prognostic factors for these children helps select appropriate treatment for their disease, ensuring maximum clinical benefit along with minimum long‐term toxicity [32]. In the current study, we show that positive engraftment of tumor biopsies in mice is a prognostic factor related to poorer EFS for newly diagnosed pediatric patients with the most frequent being bone or soft tissue sarcomas.

All patients in the diagnostic cohort whose samples engrafted in mice relapsed and died of their disease while we were performing this study. This discovery might be especially relevant for newly diagnosed Ewing sarcoma patients, because PDX engraftment predicted relapse in patients otherwise classified as standard risk at diagnosis. Thus, our results suggest that PDX engraftment might help unveil patients at higher risk and stratify them accordingly.

The engraftment rate calculated in this study, 44%, is similar to previous studies of pediatric solid malignancies both for primary and relapsed tumors [23, 33, 34]. Patient determinants favoring tumor biopsy engraftment in mice are not yet clear. We show here that age older than 12 years or relapsed disease correlates highly with pediatric sarcoma biopsy engraftment. Age and relapsed disease are also determining factors for engraftment of other tumors such as lung adenocarcinoma and hepatoblastoma [15, 35]. In adult patients, advanced disease stage increases the likelihood of PDX engraftment from patients with lung cancer [36]. Also in adult patients, tumor mutations, such as in EGFR, BRCA1/BRCA2, and SMAD4 genes, decrease engraftment rates in adenocarcinoma, breast cancer, and pancreatic cancer, respectively [14, 15, 18].

Other studies have evaluated some of the technical factors addressed in our study, such as sample conservation. In agreement with our results, cryopreserved pancreatic cancer biopsies engrafted efficiently in mice [37]. We did not evaluate ‘engraftment site’ as a factor, but the orthotopic site is especially relevant to reproduce the complex microenvironment necessary for the growth of central nervous system tumors [38, 39, 40], and could facilitate and accelerate tumor implantation and recapitulate patient tumor characteristics [41, 42]. However, our selection of the subcutaneous site is justified by its advantages of offering a simple and fast surgical procedure and easier follow‐up of the mice, as it does not require imaging systems to visualize tumor site implantation [13].

Interestingly, we found that the Tru‐Cut technique, a less invasive method than the surgical collection of tumor samples, was successful for the purpose of PDX establishment, and resulted in more rapid patient recovery [43]. Patients for whom surgery would constitute an excessive risk or associated morbidity, such as relapsed patients with advanced disease, are more suitable for the Tru‐Cut method. On the other hand, the Tru‐Cut method is at a disadvantage, due to its resulting smaller sample size and tumor heterogeneity and necrosis in some tumor areas; therefore, the Tru‐Cut technique can result in greater challenges to pathologists and cancer biologists.

Our studies confirm the similarities of histology and genetics of the PDX and the original biopsy, consistent with previous data showing that CNA variations in patient biopsies were conserved in tumor xenografts of breast cancer and neuroblastoma [17, 44]. It was not within the scope of our project to perform a comprehensive comparison of paired PDX models coming from the same patient. In a previous study, we compared systematically Ewing sarcoma PDX pairs obtained from the same patients, finding functional differences among them that we attributed to patient treatments [10].

Other studies have shown changes in tumor architecture and chromosomal stability during successive passages over time that could differ from those arising in patients [13, 45]. Breast cancer tumors and their corresponding xenografts share mostly all single‐nucleotide variations, with the exception of some clusters that are dominant in the PDX but not in the tumor patients [17]. These results, among others, suggest that there exists clonal selection during xenograft establishment [21, 46, 47]. It is therefore not clear whether such alterations affect the capacity of the PDX to represent the original patient biopsy from the functional perspective, such as response to treatment. Our studies of irinotecan activity in PDX models at different stages of their engraftment partially addressed this question. We confirmed the stability of our models during successive engraftments and their ability to replicate patient response to treatment at the preclinical dosage of 10 mg/kg using the protracted dosing schedule [30]. It is likely that this dosage provides a higher SN‐38 systemic exposure than the one achievable in humans, due to the better conversion of the prodrug irinotecan in mice [48], which could explain the different response to irinotecan found between the PDX and one of our patients. Nevertheless, our data support that PDX might predict clinical efficacy of personalized medicine while keeping stable tumor architecture and heterogeneity [11, 12, 13]. Performing drug efficacy studies using the PDX models could aid in the design of personalized therapies for our patients. The main limiting factor is the lag time between the necessary time for the biopsy to engraft (9–17 weeks, according to our experience) and the time‐critical needs of the patients. Disease progression is usually faster than the time needed to engraft and expand tumors, and not all patients can benefit from this approach. According to our experience, we need at least 28 weeks to obtain preclinical results for candidate individualized treatments. This number is obtained by adding the median times to achieve engraftment of the F0 generation (13 weeks; range: 9–17 weeks), median times for expanding the F1 generation in a sufficient number of mice to start candidate treatments (13 weeks; range: 6–20 weeks), and times to treat and evaluate treatment activity (2 weeks typically). Although personalized treatments were not within the scope of our work, in our cohort approximately 50% (13) of patients with associated stable PDX progressed later than 28 weeks after biopsy, and could have obtained a theoretical benefit from preclinical assays using their PDX.

In conclusion, our study suggests that PDX engraftment at the time of patient diagnosis helps to identify aggressive pediatric sarcoma tumors with poorer prognosis, predicting disease progression. Such PDX models stably represent patient disease, both morphologically and functionally. Whether this systematic study of the PDXs from our patients has future clinical implications for personalized medicine, patient risk classification, and intensification of treatments should now be addressed in prospective clinical trials.

Author contributions statement

All authors were involved in data collection and manuscript editing. HC‐E, SP‐J, CL and AMC analyzed the data. JM and AMC conceived the work, acquired funding, and provided resources. HC‐E and AMC wrote the original draft and generated the figures. AMC supervised the project.

Figure S1. Staining of human cells in six representative cases of xenografts at early and late passages

Figure S2. Comparative histology of two PDX from the same rhabdomyosarcoma patient at relapse and necropsy

Figure S3. Analysis of cell proliferation in human tumor biopsies and their corresponding PDXs

Figure S4. Whole‐genome visualization of the chromosomal profiles of patients with Ewing sarcoma and rhabdomyosarcoma and their matched PDX models

Table S1. Clinical data of patients with Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma

Table S2. Data from PDX established from patients with Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma

Table S3. Clinical information from patients with Ewing sarcoma (N = 31; median age: 12.3 years; range: 0.52–17.9 years)

Table S4. Clinical information from patients with osteosarcoma (N = 10; median age: 11.1 years; range: 5.67–14.7 years)

Table S5. Clinical information from patients with rhabdomyosarcoma (N = 10; median age: 9.7 years; range: 4.03–15.4 years)

Table S6. Engraftment outcomes of patient samples

Table S7. Association of technical factors with engraftment

Table S8. Association of previous biopsy engraftment with the engraftment of the following biopsy obtained from the same patient

Table S9. Association of patient factors with engraftment in patients with Ewing sarcoma

Table S10. Clinical data from PDX models included in the efficacy study of irinotecan

Table S11. PDX response to treatment with a single cycle of irinotecan

Table S12. Comparative response to irinotecan‐based treatments of rhabdomyosarcoma patients and their corresponding PDX

Table S13. Comparative response to irinotecan treatment and median survival of PDXs at F ≤ 2 or F ≥ 6 passages

Click here for additional data file.