Abstract Background There is resurgence within drug and biomarker development communities for the use of primary tumorgraft models as improved predictors of patient tumor response to novel therapeutic strategies. Despite perceived advantages over cell line derived xenograft models, there is limited data comparing the genotype and phenotype of tumorgrafts to the donor patient tumor, limiting the determination of molecular relevance of the tumorgraft model. This report directly compares the genomic characteristics of patient tumors and the derived tumorgraft models, including gene expression, and oncogenic mutation status. Methods Fresh tumor tissues from 182 cancer patients were implanted subcutaneously into immune-compromised mice for the development of primary patient tumorgraft models. Histological assessment was performed on both patient tumors and the resulting tumorgraft models. Somatic mutations in key oncogenes and gene expression levels of resulting tumorgrafts were compared to the matched patient tumors using the OncoCarta (Sequenom, San Diego, CA) and human gene microarray (Affymetrix, Santa Clara, CA) platforms respectively. The genomic stability of the established tumorgrafts was assessed across serial in vivo generations in a representative subset of models. The genomes of patient tumors that formed tumorgrafts were compared to those that did not to identify the possible molecular basis to successful engraftment or rejection. Results Fresh tumor tissues from 182 cancer patients were implanted into immune-compromised mice with forty-nine tumorgraft models that have been successfully established, exhibiting strong histological and genomic fidelity to the originating patient tumors. Comparison of the transcriptomes and oncogenic mutations between the tumorgrafts and the matched patient tumors were found to be stable across four tumorgraft generations. Not only did the various tumors retain the differentiation pattern, but supporting stromal elements were preserved. Those genes down-regulated specifically in tumorgrafts were enriched in biological pathways involved in host immune response, consistent with the immune deficiency status of the host. Patient tumors that successfully formed tumorgrafts were enriched for cell signaling, cell cycle, and cytoskeleton pathways and exhibited evidence of reduced immunogenicity. Conclusions The preservation of the patient’s tumor genomic profile and tumor microenvironment supports the view that primary patient tumorgrafts provide a relevant model to support the translation of new therapeutic strategies and personalized medicine approaches in oncology. Keywords: Tumorgrafts, Genomics, Translational models Background Although direct transfer xenografts of fresh human tumors, or primary “tumorgrafts”, have been reported as early as the 1970s for testing of new pharmaceutical agents [[54]1-[55]4], only recently has there been a resurgence of interest in this alternative to the more traditional cell line xenograft models. This has been due, in part, to the growing realization that drugs which work in the traditional cell line xenograft models rarely exhibit comparable efficacy in patients with the anatomically/pathologically equivalent tumor [[56]5-[57]7]. Indeed, it has been highlighted that the failure of drugs during clinical trials is linked to a lack of testing in clinically relevant preclinical models [[58]8]. The advent of “precision medicine” (PMed), utilizing novel molecular techniques to design patient tumor-specific therapeutic regimens, requires preclinical cancer models that closely reflect the originating human disease to allow evaluation and rapid translation of these molecular based approaches to therapy selection in the clinic [[59]9]. While there have been a growing number of studies reporting comparative chemotherapeutic responses between tumorgraft models and patients [[60]10-[61]13], there is limited data addressing genomic similarities between tumorgraft models and the originating patient tumors [[62]14-[63]16]. Previous studies have reported excellent concordance in expression profiles between patient tumors and tumorgraft models of similar histotypes [[64]17,[65]18]; however, the comparisons were not made in a pairwise fashion between each tumorgraft and the originating tumor from the donor patient. Given the added cost and time required to develop primary patient tumorgrafts relative to the more classical cell line xenografts, a detailed cellular and molecular characterization of the tumorgraft models showing a high degree of equivalence to the originating patient’s tumor in the human host could provide the justification for the routine use of this model, even at the ectopic subcutaneous site, in translational studies and personalized medicine applications. Our goal was to compare the genomic profiles of a panel of diverse tumorgraft models to the patient tumor tissue from which the models were derived. To address inter-disease variability, the panel of tumorgraft models selected for this study contains a spectrum of neoplastic diseases, rather than a limited set of disease types [[66]14-[67]16,[68]19]. Our long term objective is to develop this panel of mixed-type tumorgraft models to evaluate new treatment strategies based on the molecular characterization of an individual disease and the subsequent treatment with targeted therapy, independent of the tumor histology and anatomical location. Methods Human subjects Ethical approval for this study was granted by the Institutional Review Boards at both the Van Andel Research Institute (VARI) and the collaborating clinical institutions: Spectrum Health Hospitals, Grand Rapids, MI; Saint Mary's Health Care, Grand Rapids, MI; Oncology Care Associates, St. Joseph, MI; Mary Crowley Cancer Research Centers, Dallas, TX. Written, informed consent was obtained from each patient or guardian prior to study enrollment. Eligibility for tumorgraft development was not limited by cancer type, cancer staging, or anatomical location. The tumor tissue used for tumorgraft development was deemed excess to that required for the patient’s diagnosis and standard of care and treatment. Tumor processing The tumor tissue was subdivided into a portion for genomic analysis and a portion for implantation into immune-compromised mice. Tissue for genomic analysis was snap frozen (30–60 s) in chilled isopentane (dry ice) and stored at −80° C. Tumor tissue for implantation into mice was maintained in 6 ml transport media, [RPMI 1640 media (Invitrogen, Carlsbad, CA), 1% penicillin/streptomycin (Invitrogen), and 50 units of heparin (Sigma, St. Louis, MO)/ml media] at 4 °C until implantation. Samples from distant sites were shipped on wet ice and were implanted within 24 hours of surgical resection. Tumorgraft development All animal studies were approved by the VARI Institutional Animal Care and Use Committee (IACUC) and used 6–8 week old athymic nu/nu mice from the VARI breeding colony. Food and water was available ad libitum for the duration of the studies. Mice for each tumorgraft model were gender matched to the donor patient. Body weights of the mice were recorded weekly during tumorgraft development. Tumorgraft volumes (½ x length x depth x height) were measured 1x/week when volumes ≤50 mm^3 and 3x/week at tumor volume >50 mm^3. Mice were euthanized and subcutaneous tumorgrafts harvested following IACUC guidelines. Upon receipt, the tumor tissue for implantation was placed into a sterile dish containing sterile phosphate buffered saline (Invitrogen) and carefully teased into ≤3 millimeters (longest axis) tumor fragments. Dependent on tumor tissue availability, tumor fragments were implanted in a maximum of five mice (1^st generation). Following administration of general anaesthesia (isoflurane), the right flank was cleaned with 70% ethyl alcohol, a small incision made, and a subcutaneous pocket created by blunt dissection. The tumor fragment was inserted into the pocket and the incision closed using a surgical staple. Immediately following surgery, the mouse received a single dose of the analgesic Ketoprofen (5 mg/kg body weight). Mice were monitored for health and tumor growth for the duration of the study. A tumorgraft model that failed to develop within 6 months in the 1^st generation mice was discontinued and the mice euthanized. When a 1^st generation tumorgraft reached a volume of ≥1500 mm^3 the mouse was euthanized and the tumorgraft was aseptically harvested. The tumorgraft was subsequently divided into pieces for four applications: 1) 5 pieces (≤3 millimeters (longest axis) were directly transplanted into 2^nd generation mice; 2) 18 pieces (≤3 millimeters (longest axis) were cryopreserved (RPMI 1640 media (Invitrogen), 10% fetal bovine serum (Invitrogen), 10% dimethyl sulfoxide (Sigma, St. Louis, MO), and 50 units heparin (Sigma)/ml RPMI 1640 media), using a “Mr. Frosty” Freezing Container (Thermo Scientific, Waltham, MA) at −1°/min cooling rate, in a −80 °C freezer, to allow for subsequent reestablishment of the tumorgraft model; 3) a single piece (3 mm x 5 mm x 5 mm) was snap frozen for future genomic/proteomic analysis; and 4) remaining tissue was formalin fixed. The 2^nd generation mice were monitored for health and tumorgraft growth characteristics as with the 1^st generation mice. Once a total of ~50 tumorgraft fragments, from either 1^st or 2^nd generation mice, had been cryopreserved, all remaining mice were euthanized and any tumorgrafts harvested (snap frozen and formalin fixed) and stored for subsequent analyses. Two osteosarcoma models were grown out to a 4^th generation and used to examine longitudinal changes such as genetic drift and mutation status. For a tumorgraft model to be considered fully established the following three criteria needed to be met: 1) 1^st generation growth, 2) 2^nd generation tumor growth following direct (fresh) transplant of a tumor fragment from a 1^st generation mouse, and 3) re-growth of the cryopreserved tumor fragment. To test cryo-viability, three 1^st generation tumor fragments, frozen in liquid nitrogen for a minimum of four weeks, were thawed rapidly, rinsed in sterile phosphate buffered saline (Invitrogen) containing 1% penicillin/streptomycin (Invitrogen), and implanted into immune-compromised mice, as described above. In the event that the cryopreserved tissue failed to develop into a tumorgraft within 6 months, the model was deemed unsuccessful and the mice euthanized. While the preservation of the histopathological features of the primary patient tumor in the resulting tumorgraft models was assessed, it was not a requirement for a successful model development. Data storage and management For this study, all clinical, preclinical, histological, and genomic data was stored and analyzed using the XenoBase-BioIntegration Suite (XB-BIS) ([69]http://xbtransmed.com/). Data analysis is accomplished by integrating molecular, preclinical and clinical data, facilitating integrated analysis within this single platform as described in the pertinent sections below. Histological analysis of tumorgrafts and patient tumors Histopathological evaluation was performed on five micron hematoxylin and eosin (H&E) stained sections of formalin-fixed paraffin-embedded (FFPE) pieces of patient tumors and of subsequent 1^st and 2^nd generation tumorgrafts. Tissue sections were reviewed by two independent pathologists (BJN & RS) to compare the pathological features of the tumorgrafts and compared directly with those of the corresponding patient tumor. Mutational analysis using the sequenom OncoCarta assay For mutation detection, the OncoCarta mutation panel was used according to manufacturer’s protocol (Sequenom, San Diego, CA). The OncoCarta™ Panel v1.0 (Sequenom, San Diego, CA) has the capacity to detect 238 mutations in 19 oncogenes. Ten representative pairs of 1^st generation tumorgrafts and originating patient tumors and 1^st - 4^th generation series of an osteosarcoma tumorgraft model were screened using the MassARRAY system and OncoCarta Panel v1.0 (Sequenom, San Diego, CA) [[70]20,[71]21]. Tumor DNA was isolated from either snap frozen of formalin-fixed paraffin embedded tissue using an initial Proteinase K digestion followed by salt ethanol precipitation. The DNA was further purified using the QIAamp DNA micro Kit (Qiagen, Valencia, CA). Only samples in which the H&E sections, from above and below the processed tissue, contained ≥50% of viable malignant cell content (assessed by nuclei vs. surface area) were submitted for analysis of oncogenic mutations. The spectra from the MassARRAY Compact Analyzer were loaded into MassARRAY Typer Analyzer software 4.0.4.20 (Sequenom). The analysis software automates the identification of mutants by comparing ratios of the wild type peak to that of all suspected mutants and adjusting these peaks for adducts it detects in the spectrum [[72]20,[73]21]. All mutations detected by this automated process were visually reviewed to identify "real" mutant peaks and to remove all artifact peaks due to salt peaks or other background peaks. Gene expression analysis Transcriptome profiling of the genomes of 24 representative 1^st generation tumorgraft/patient tumor pairs, encompassing gastrointestinal, musculoskeletal, skin, head and neck, and genitourinary cancers, and genomes of two musculoskeletal models across four tumorgraft generations was performed using the Human U133 Plus 2.0 GeneChip® (Affymetrix, Inc, Santa Clara, CA). The resulting gene expression data were stored and analyzed in XB-BIS. RNA was isolated from snap-frozen patient tumors and tumorgraft tissue that was assessed for viable malignant cell content, as described above. RNA was isolated by chloroform extraction, precipitation in isopropyl alcohol, a 70% ethanol wash and suspension in nuclease free water. A final RNA purification was performed using RNeasy Mini-/Micro-Kit^TM (Qiagen) with an on-column DNase 1 digestion. RNA quantity and purity was determined using a NanoDrop Spectrophotometer (NanoDrop, Wilmington, DE); RNA quality was determined with an Agilent 2100 Bioanalyzer (Agilent, Foster City, CA). Labeled one round cDNA was generated either by the Affymetrix method: Gene-Chip One-Cycle Target Labeling and Control Reagent (Affymetrix) or by the NuGEN method: WT-Ovation Pico RNA Amplification system and FL-Ovation cDNA Biotin Module (NuGEN Technologies, Inc., San Carlos, CA). The resulting amplification products were quantified and hybridized to GeneChip® Human Genome U133 Plus 2.0 Arrays (Affymetrix, Inc). Hybridization, washing, staining and scanning of the arrays were performed according to manufacturer’s instructions (Affymetrix, Inc.). The resulting .cel files were loaded into the Affymetrix Expression Console Version 1.1, batch processed using the Affymetrix MAS5 algorithm with a target value of 500. The resulting pivot table of normalized expression intensities was then exported into XB-BIS for subsequent analysis. Comparative analysis of genomes of tumorgrafts and matched patient tumors Genomic similarity between the tumorgraft and its matched patient tumor were assessed by determination of Pearson correlation coefficients and by unsupervised cluster analysis. For determination of Pearson correlation coefficients, gene expression data was filtered to remove those Affymetrix probes deemed absent and/or with normalized intensities less than 100 in ≥50% of the sample set of 24 matched tumor/tumorgraft pairs (n = 25,806 probe sets); the gene expression data was not censored for Affymetrix x_at probes. For unsupervised cluster analysis, the selection of the significant variance probes was made to emphasize the differences in the matched tumor pairs, which is diminished by large numbers of low intensity, non-significant probes that are similar both between pairs and across the tumor panel. The probe sets were filtered to include only probe sets with an intensity ≥150, called present in ≥20% samples, and exhibited significant variance (coefficient of variance >1). Affymetrix x_at probes were excluded from the final data used for the cluster analysis. Initial unsupervised clustering of the 24 tumorgraft/patient tumors, identified that the choice of RNA amplification methods, Affymetrix and NuGEN, was a significant confounding variable (data not shown). Subsequently, the 48 data sets were clustered independently, based on method of RNA amplification (Affymetrix, 28 genomes; NuGEN, 20 genomes). The resulting Affymetrix amplification probe set (n = 6,845) and NuGEN amplification probe set (n = 6,265) were loaded into the programming language R version 2.11.1. The distance matrix was created using Euclidean distance default parameters and complete linkage with the default parameters were utilized for the hclust function. The unrooted dendrogram was constructed using the phylo object and plotted with default parameters. To determine whether the differential expression in the genomes between the tumorgrafts and the originating patient tumors was random or represent a gain or loss of specific pathways, a paired t-test (p < 0.05) was employed on the set of 25,806 probes, generated above, to identify those genes that were differentially expressed, 2-fold or greater, between tumorgrafts and their matched tumors. No multiple testing correction was applied because the sequential pathway enrichment analysis is expected to remove or diminish the effects of random significant gene selection. The resulting genes up-regulated (n = 17) or down-regulated (n = 395) in tumorgrafts, as compared to the originating patient tumors (Additional file [74]1: Table S1), were analyzed for canonical pathway enrichment and ontology enrichment using the MetaCore™ pathway analysis suite (GeneGo-Thomson Reuters, St. Joseph, MI) using p-values of hypergeometric distribution as enrichment and ranking method [[75]22]. GeneGo’s data import methods automatically resolve duplicate gene mapping and unmapped Affymetrix probes. The attrition in the mapping of the Affymetrix probes in Additional file [76]1: Table S1 ranged from 30% to 35%. In addition, the gene lists were analyzed for enrichment in broadly defined protein functions, i.e. transcription factors, secreted proteins, kinases, receptors, proteases, etc. using MetaCore™’s Interactome workflow. Analysis of genomes of tumors that successfully formed tumorgrafts compared to tumors that failed to form tumorgrafts To begin to understand why all donor tumors do not form tumorgrafts following transplant into immune-compromised mice, the genomes of 16 patient tumors that did form tumorgrafts were compared to 16 tumors that did not form tumorgrafts. These 32 genomes were selected since they were all amplified by the Affymetrix method, thus removing the confounding effect of the RNA amplification method. For the analysis of the tumorgraft-forming tumor genomes against the non-tumorgraft forming tumor genomes, a student’s two-tailed t-test (p < 0.05) was employed to identify genes that were differentially expressed, 2-fold or greater, between the two classes. No multiple testing correction was applied because the sequential pathway enrichment analysis is expected to remove or diminish the effects of random significant gene selection. The resulting list of up-regulated (n = 491) and down-regulated (n = 691) genes in the tumorgraft-forming tumors, compared to the non-tumorgraft forming tumors (Additional file [77]2: Table S2), were submitted to MetaCore™ for enrichment analysis and interaction analysis, as described above. The attrition in the mapping of the Affymetrix probes in Additional file [78]2: Table S2 ranged from 28% to 36%. Results Forty-nine tumorgraft models, representing 18 distinct cancer pathologies have been developed by direct implantation of patient tumor tissue into immune-compromised nude mice giving the overall take rate of 27% (Table [79]1). These 49 models were distributed across eight different cancer location categories, based on the NCI Cancer by Body Location/System (Table [80]2). Tumorgraft models of the gastrointestinal tract and skin account for ~50% of successfully established tumorgraft models. There were an additional three models that formed 1^st generation tumorgrafts but that failed to re-establish tumorgrafts following cryopreservation of 1^st generation tumorgraft fragments, giving a 94% success rate in tumorgraft model re-establishment following cryopreservation. No association was identified between tumorgraft development and the time between surgical resection and implantation (p = 0.86 by simple t-test on representative sample set, data not shown). Table 1. Summary of Tumorgraft Development Location Tumor type # Models Success Failure Success (%) Breast __________________________________________________________________ Breast Cancer __________________________________________________________________ 55 __________________________________________________________________ 3 __________________________________________________________________ 52 __________________________________________________________________ 5 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 55 __________________________________________________________________ 3 __________________________________________________________________ 52 __________________________________________________________________ 5 __________________________________________________________________ Endocrine __________________________________________________________________ Neuroendocrine Tumor __________________________________________________________________ 3 __________________________________________________________________ 1 __________________________________________________________________ 2 __________________________________________________________________ 33 __________________________________________________________________ __________________________________________________________________ Papillary Thyroid Carcinoma __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 5 __________________________________________________________________ 1 __________________________________________________________________ 4 __________________________________________________________________ 20 __________________________________________________________________ Gastrointestinal __________________________________________________________________ Cecal Adenocarcinoma __________________________________________________________________ 1 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Cholangiocarcinoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Colorectal Cancer __________________________________________________________________ 17 __________________________________________________________________ 8 __________________________________________________________________ 9 __________________________________________________________________ 47 __________________________________________________________________ __________________________________________________________________ Esophageal Adenocarcinoma __________________________________________________________________ 5 __________________________________________________________________ 2 __________________________________________________________________ 3 __________________________________________________________________ 40 __________________________________________________________________ __________________________________________________________________ Gallbladder Carcinoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Gastric Carcinoma __________________________________________________________________ 3 __________________________________________________________________ 1 __________________________________________________________________ 2 __________________________________________________________________ 33 __________________________________________________________________ __________________________________________________________________ Hepatoblastoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Pancreatic Cancer __________________________________________________________________ 6 __________________________________________________________________ 3 __________________________________________________________________ 3 __________________________________________________________________ 50 __________________________________________________________________ __________________________________________________________________ Small Bowel Cancer __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 37 __________________________________________________________________ 15 __________________________________________________________________ 22 __________________________________________________________________ 41 __________________________________________________________________ Genitourinary __________________________________________________________________ Bladder Carcinoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Clear Cell Sarcoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Prostate Adenocarcinoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Renal Cell Carcinoma __________________________________________________________________ 4 __________________________________________________________________ 0 __________________________________________________________________ 4 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Urachal Adenocarcinoma __________________________________________________________________ 1 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Wilms Tumor __________________________________________________________________ 1 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 9 __________________________________________________________________ 2 __________________________________________________________________ 7 __________________________________________________________________ 22 __________________________________________________________________ Gynecologic __________________________________________________________________ Ovarian Carcinoma __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ Head and Neck __________________________________________________________________ Head and Neck Cancer __________________________________________________________________ 2 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Nasopharyngeal Cancer __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Squamous Cell Carcinoma of the mouth __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 4 __________________________________________________________________ 2 __________________________________________________________________ 2 __________________________________________________________________ 50 __________________________________________________________________ Leukemia __________________________________________________________________ Acute Lymphoblastic Leukemia __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ Lung __________________________________________________________________ Squamous Cell Carcinoma of the lung __________________________________________________________________ 5 __________________________________________________________________ 0 __________________________________________________________________ 5 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Lung Adenocarcinoma __________________________________________________________________ 4 __________________________________________________________________ 1 __________________________________________________________________ 3 __________________________________________________________________ 25 __________________________________________________________________ __________________________________________________________________ Small Cell Lung Cancer __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ NSCLC __________________________________________________________________ 13 __________________________________________________________________ 5 __________________________________________________________________ 8 __________________________________________________________________ 38 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 23 __________________________________________________________________ 6 __________________________________________________________________ 17 __________________________________________________________________ 26 __________________________________________________________________ Lymphoma __________________________________________________________________ Burkitts Lymphoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Follicular Lymphoma __________________________________________________________________ 3 __________________________________________________________________ 0 __________________________________________________________________ 3 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Hodgkins Lymphoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Lymphoangioma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Lymphoma __________________________________________________________________ 6 __________________________________________________________________ 0 __________________________________________________________________ 6 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 12 __________________________________________________________________ 0 __________________________________________________________________ 12 __________________________________________________________________ 0 __________________________________________________________________ Musculoskeletal __________________________________________________________________ Desmoid Tumor __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Ewings Sarcoma __________________________________________________________________ 3 __________________________________________________________________ 2 __________________________________________________________________ 1 __________________________________________________________________ 67 __________________________________________________________________ __________________________________________________________________ Hemangiopericytoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Leiomysarcoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Osteosarcoma __________________________________________________________________ 3 __________________________________________________________________ 3 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Rhabdomyosarcoma __________________________________________________________________ 3 __________________________________________________________________ 1 __________________________________________________________________ 2 __________________________________________________________________ 33 __________________________________________________________________ __________________________________________________________________ Synovial Sarcoma __________________________________________________________________ 1 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 13 __________________________________________________________________ 7 __________________________________________________________________ 6 __________________________________________________________________ 54 __________________________________________________________________ Neurologic __________________________________________________________________ Neuroblastoma __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ 1 __________________________________________________________________ 0 __________________________________________________________________ Skin __________________________________________________________________ Melanoma __________________________________________________________________ 18 __________________________________________________________________ 11 __________________________________________________________________ 7 __________________________________________________________________ 61 __________________________________________________________________ __________________________________________________________________ Merkel Cell Carcinoma __________________________________________________________________ 2 __________________________________________________________________ 2 __________________________________________________________________ 0 __________________________________________________________________ 100 __________________________________________________________________ __________________________________________________________________ Total __________________________________________________________________ 20 __________________________________________________________________ 13 __________________________________________________________________ 7 __________________________________________________________________ 65 __________________________________________________________________ Final Total 182 49 133 27 [81]Open in a new tab Table 2. Distribution of Tumorgraft Models across Cancer Location Categories and Staging Location Tumor type # Models Staging __________________________________________________________________ I II III IV ND Breast __________________________________________________________________ Breast Cancer __________________________________________________________________ 3 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 2 __________________________________________________________________ 1 __________________________________________________________________ Endocrine __________________________________________________________________ Neuroendocrine Tumor __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Gastrointestinal __________________________________________________________________ Cecal Adenocarcinoma __________________________________________________________________ 1 __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Colorectal Cancer __________________________________________________________________ 9 __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 7 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Esophageal Adenocarcinoma __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Gastric Carcinoma __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Pancreatic Cancer __________________________________________________________________ 3 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 3 __________________________________________________________________ __________________________________________________________________ Genitourinary __________________________________________________________________ Wilms Tumor __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ Urachal Adenocarcinoma __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ Head and Neck __________________________________________________________________ Head and Neck Cancer __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Lung __________________________________________________________________ Lung Carcinoma __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ NSCLC __________________________________________________________________ 5 __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ 2 __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ Musculoskeletal __________________________________________________________________ Ewings Sarcoma __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 2 __________________________________________________________________ __________________________________________________________________ Osteosarcoma __________________________________________________________________ 3 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 3 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Rhabdomyosarcoma __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ Synovial Sarcoma __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ Skin __________________________________________________________________ Melanoma __________________________________________________________________ 11 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 4 __________________________________________________________________ 6 __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ Merkel Cell Carcinoma __________________________________________________________________ 2 __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ 1 __________________________________________________________________ __________________________________________________________________ Total 49 [82]Open in a new tab The take rate (the ability of donor tumor to successfully be propagated in mice) was closely related to tumor stage at time of resection. Forty-three percent (35/82) of tumors with a pathological staging of Stage III and Stage IV tumors were successfully propagated in mice, accounting for 73% of fully developed models (Table [83]2). Only 10% (6/60) of the donor tumors from Stage I or II cancers successfully propagated formed tumorgrafts in mice, accounting for only 12% of the fully developed models. Pathological staging information for the remaining 22% of the tumor tissues implanted into mice (40/182) which formed tumorgrafts was not obtained. A contingency table analysis with stage III and IV as one condition and stage I and II the second condition supported the hypothesis with a p-value = 1.2 E-5 and an odds ratio of 6.7 for this data set (data not shown). Tumorgraft development was also found to be related to cancer type. Of the human tumor samples used to derive the panel of 49 tumorgraft models, a success rate of >40% was observed for skin (65%), musculoskeletal (54%), head and neck cancers (50%), and gastrointestinal cancers (41%) (Table [84]1). Lung (26%), genitourinary (22%), and endocrine (20%) cancers exhibited lower take rates for tumorgraft formation. Cancers of the breast (5%), predominantly early stage (data not shown), and the lymphatic system (0%) proved the most challenging to form tumorgrafts using the athymic nude mouse strain and methodologies adopted here. No tumorgrafts were established from tumors of gynecologic, leukemic, or neurologic origins; however, these finding are not necessarily representative of the true take rate for these tumor types as ≤2 attempts for each of these cancer types were attempted. Hematoxylin and eosin stained sections from matching patient tumors and 1^st and 2^nd generation tumorgrafts from pancreatic cancer, breast cancer, lung cancer and cervical cancer models were examined for histological similarities. In all 4 different tumorgraft models, the initial patient tumors were characterized by histological features consistent with the anatomical site of origin (Figure [85]1 A-D). The 1^st generation tumor grafts resembled the original patient tumor, including in most instances retention of peritumoral stroma which became progressively replaced in the 2^nd generation tumorgrafts by malignant cell overgrowth. In addition, the relative state of differentiation of the respective cancer type was retained following engraftment, including production of mucin and retention of atypical glandular structures for the adenocarcinomas, even after 302 and 177 days for the pancreatic and breast cancers, respectively. In the lung cancer and cervical cancer models, the 1^st and 2^nd generation tumors revealed viable malignant appearing cells. The 2^nd generation tumorgrafts retained essential diagnostic features at 286 and 345 days, respectively, for the lung and cervical cancer models. Figure 1. [86]Figure 1 [87]Open in a new tab Histological resemblance between patient tumor and resulting tumorgraft following sequential grafting into immunocompromised mice. Representative sections stained with hematoxylin and eosin of indicated tumor samples (Scale bars indicate 100 μm). A) Adenocarcinoma of pancreas with desmoplastic stromal reaction in original patient tumor. 1^st generation tumorgraft (total time in vivo, 217 days) with preservation of atypical glandular structures with inspissated mucinous material surrounded by fibrous stroma. 2^nd generation tumorgraft (303 days) with similar histology. B) Adenocarcinoma of colon with numerous atypical glands containing mucinous material and necrotic debris with surrounding fibrovascular connective tissue. 1^st generation tumorgraft after 101 days reveals dilated atypical ducts with mucin and necrotic debris in lumens comprised of multiple columnar layers of atypical epithelium surrounded by fibrovascular stroma. 2^nd generation graft after 177 days with similar histology. C) Lung carcinoma with prominent fibrovascular connective tissue in original patient neoplasm that becomes progressively overgrown by malignant cells in 1^st (216 days) and 2^nd (286 days) generation tumor grafts. D) Cervical carcinoma with fibromucinous stroma in original patient neoplasm, with expansion of large islands of malignant cells surrounded by fibrovascular stroma in 1^st (282 days) and 2^nd (345 days) generation tumorgrafts. Of the 20 samples (10 patient samples, 10 1^st generation tumorgraft samples) submitted for oncogene mutation analysis, only 15 mutations were identified. Four tumorgraft models each had a single point mutation that was also observed in the originating tumor tissue. Models VARI-LTM-044 colorectal cancer and VARI-LTM-034 pancreatic cancer both had a single KRAS mutation; melanoma models VARI-LTM-086 and VARI-LTM-027 had BRAF and NRAS mutations, respectively (Table [88]3). Model VARI-LTM-041 colorectal cancer contained three mutations (KRAS and two PIK3CA mutations) that were also present in the matching tumor. The colorectal cancer model VARI-LTM-026, did not exhibit mutations in the patient tumor, but showed a single KRAS mutation in the tumorgraft. In the osteosarcoma model VARI-LTM-021, no mutations were identified in the patient tumor with any genetic drift identified in the 1^st - 4^th generation tumorgrafts (data not shown). Table 3. Oncogene mutations detected in tumorgrafts and originating patient tumors Model Tumor Type Gene Mutation Primary Tumor Mutation Status Tumorgraft Mutation Status VARI-LTM-041 __________________________________________________________________ Colorectal Cancer __________________________________________________________________ KRAS __________________________________________________________________ G12V __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-041 __________________________________________________________________ Colorectal Cancer __________________________________________________________________ PIK3CA __________________________________________________________________ R88Q __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-041 __________________________________________________________________ Colorectal Cancer __________________________________________________________________ PIK3CA __________________________________________________________________ E545K __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-044 __________________________________________________________________ Cecal Adenocarcinoma __________________________________________________________________ KRAS __________________________________________________________________ G13D __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-026 __________________________________________________________________ Colorectal Cancer __________________________________________________________________ KRAS __________________________________________________________________ G12D __________________________________________________________________ FALSE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-086 __________________________________________________________________ Melanoma __________________________________________________________________ BRAF __________________________________________________________________ V600E __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-027 __________________________________________________________________ Melanoma __________________________________________________________________ NRAS __________________________________________________________________ Q61R __________________________________________________________________ TRUE __________________________________________________________________ TRUE __________________________________________________________________ VARI-LTM-034 Pancreatic Cancer KRAS Q61H TRUE TRUE [89]Open in a new tab Paired analysis of the 24 tumorgraft/tumor pairs exhibited a high degree of similarity, with Pearson correlation coefficients ranging from r = 0.75 to 0.99 with a median coefficient of 0.93 (Table [90]4). There were only four pairs, three gastrointestinal and one lung, that had Pearson values r < 0.90; the NSCLC model VARI-LTM-023 had the Pearson correlation coefficient of r = 0.75. For this model, the originating patient tumor RNA had a RIN (RNA Integrity Number) of 4.9 implying degraded RNA, which was likely to be the primary factor driving the low correlation and would typically be excluded from the molecular analysis. Table 4. Pearson correlation coefficient of each individual tumorgraft and originating patient tumor based on Affymetrix gene expression data Cancer Location Category Tumor Type Subject ID Pearson correlation (R) Gastrointestinal __________________________________________________________________ Cecal Adenocarcinoma __________________________________________________________________ VARI-LTM-044 __________________________________________________________________ 0.87 __________________________________________________________________ __________________________________________________________________ Colorectal Cancer __________________________________________________________________ VARI-LTM-026 __________________________________________________________________ 0.86 __________________________________________________________________ __________________________________________________________________ Colorectal Cancer __________________________________________________________________ VARI-LTM-034 __________________________________________________________________ 0.97 __________________________________________________________________ __________________________________________________________________ Colorectal Cancer __________________________________________________________________ VARI-LTM-041 __________________________________________________________________ 0.99 __________________________________________________________________ __________________________________________________________________ Gastric Carcinoma __________________________________________________________________ VARI-LTM-024 __________________________________________________________________ 0.86 __________________________________________________________________ __________________________________________________________________ Pancreatic Cancer __________________________________________________________________ VARI-LTM-035 __________________________________________________________________ 0.92 __________________________________________________________________ Genitourinary __________________________________________________________________ Wilms Tumor __________________________________________________________________ VARI-LTM-022 __________________________________________________________________ 0.98 __________________________________________________________________ Head and Neck __________________________________________________________________ Head and Neck Cancer __________________________________________________________________ VARI-LTM-042 __________________________________________________________________ 0.94 __________________________________________________________________ __________________________________________________________________ Head and Neck Cancer __________________________________________________________________ VARI-LTM-045 __________________________________________________________________ 0.90 __________________________________________________________________ Lung __________________________________________________________________ Lung Cancer __________________________________________________________________ VARI-LTM-029 __________________________________________________________________ 0.95 __________________________________________________________________ __________________________________________________________________ NSCLC __________________________________________________________________ VARI-LTM-023 __________________________________________________________________ 0.75 __________________________________________________________________ __________________________________________________________________ NSCLC __________________________________________________________________ VARI-LTM-046 __________________________________________________________________ 0.93 __________________________________________________________________ __________________________________________________________________ NSCLC __________________________________________________________________ VARI-LTM-038 __________________________________________________________________ 0.93 __________________________________________________________________ Musculoskeletal __________________________________________________________________ Ewings Sarcoma __________________________________________________________________ VARI-LTM-028 __________________________________________________________________ 0.97 __________________________________________________________________ __________________________________________________________________ Osteosarcoma __________________________________________________________________ VARI-LTM-020 __________________________________________________________________ 0.93 __________________________________________________________________ __________________________________________________________________ Osteosarcoma __________________________________________________________________ VARI-LTM-021 __________________________________________________________________ 0.95 __________________________________________________________________ __________________________________________________________________ Osteosarcoma __________________________________________________________________ VARI-LTM-031 __________________________________________________________________ 0.93 __________________________________________________________________ __________________________________________________________________ Rhabdomyosarcoma __________________________________________________________________ VARI-LTM-058 __________________________________________________________________ 0.94 __________________________________________________________________ __________________________________________________________________ Synovial Sarcoma __________________________________________________________________ VARI-LTM-019 __________________________________________________________________ 0.91 __________________________________________________________________ Skin __________________________________________________________________ Melanoma __________________________________________________________________ VARI-LTM-025 __________________________________________________________________ 0.96 __________________________________________________________________ __________________________________________________________________ Melanoma __________________________________________________________________ VARI-LTM-027 __________________________________________________________________ 0.92 __________________________________________________________________ __________________________________________________________________ Melanoma __________________________________________________________________ VARI-LTM-043 __________________________________________________________________ 0.97 __________________________________________________________________ __________________________________________________________________ Merkel Cell Carcinoma __________________________________________________________________ VARI-LTM-033 __________________________________________________________________ 0.97 __________________________________________________________________ Merkel Cell Carcinoma VARI-LTM-039 0.93 [91]Open in a new tab Of the 28 genomes amplified by the Affymetrix method, 12 of the 14 tumorgraft/donor tumor pairs co-clustered in an unsupervised analysis (Figure [92]2A). Of the 20 genomes amplified by the NuGEN method, 7 of 10 tumorgraft/ tumor pairs coclustered (Figure [93]2B). 4 of the 5 tumorgraft/tumor pairs that did not co-cluster in an unsupervised fashion exhibited Pearson correlation coefficients r ≤ 0.90 (Table [94]4). The tumorgraft and originating patient tumor from the NSCLC model VARI-LTM-023 co-clustered despite a low Pearson correlation coefficient of r = 0.75, suggesting significant variance in respect to other tumors in the cohort. Figure 2. [95]Figure 2 [96]Open in a new tab Unsupervised cluster analysis of genomes of tumorgrafts and originating patient tumors. Unsupervised clustering was performed on genomes amplified by A) Affymetrix method, 28 genomes total; and B) NuGEN method, 20 genomes total. Clustering was performed using programming language R version 2.11.1 on normalized gene expression data that was pre-processed using the Affymetrix MAS5 algorithm to an average intensity target of 500. Other than a few highlighted outliers discussed in the text, patient tumors and their derived tumorgrafts are closely related at the gene expression level. The genomic profiles of two representative osteosarcoma models VARI-LTM-020 and VARI-LTM-021 were observed over 4 generations to assess the potential for genomic drift over time in vivo. The Pearson correlation coefficients were very high, ranging from r = 0.98-0.99 between all 4 generations in both models (data not shown). The largest decrease in Pearson correlation coefficient occurred between the human tumor to 1^st generation mouse tumorgraft (VARI-LTM-020, r = 0.91; VARI-LTM-021, r = 0.94). Genomic fidelity was then maintained in subsequent mouse to mouse generation comparisons. This stability in the genomes across mouse tumorgraft generations, along with the absence of additional somatic oncogene mutations suggests the maintenance of genomic stability during serial in vivo passages and time in vivo for at least up to 186 days (VARI-LTM-021). Paired analysis of all 24 tumorgrafts to their matching tumors identified 25,806 gene probes that were present in at least 50% of the 48 samples. 395 genes were decreased and 17 genes increased two-fold or greater in the tumorgrafts (Additional file [97]1: Table S1) compared to the patient tumors (Figure [98]3A, green and red circles respectively). Figure 3. [99]Figure 3 [100]Open in a new tab Enrichment of canonical pathway maps in genes downregulated in tumorgrafts compared to their donor tumors.A) Scatter plot of expression of gene probes (25,806) graphically depicting the up-regulated (red circle) and down-regulated (green circle) 2-fold or greater (t-test, p ≤ 0.05) in tumorgrafts relative to the originating patient tumors. B) Ontology enrichment analysis was performed using the MetaCore™ pathway analysis suite (GeneGo-Thomson Reuters, St. Joseph, MI). Eight of the top ten canonical pathways enriched in the genes down-regulated in the tumorgrafts, relative to the originating patient tumors, were immune-related, consistent with immunodeficient status of the murine host. Of the 395 genes down-regulated in the murine tumorgrafts, there was a significant enrichment in pathways related to immune response (Figure [101]3B). Enrichment scores were extremely high, with the top-scored classical complement pathway showing a –log p-value of 25. Of the 17 genes up-regulated in the tumorgrafts, there were only 1–2 genes per map, suggesting no enrichment for any canonical pathway (data not shown). An independent GO biological pathway enrichment analysis was conducted using the down regulated probe list and 14 of the most significantly enriched top 100 pathways p < 4.0E-9 were associated with immune response [[102]23]. Classification of the down-regulated genes with regards to enrichment in protein function, identified a borderline significance (p = 0.05) in receptors and ligands in the tumorgrafts compared to the patient tumors. For the up-regulated genes, the only significant enrichment in protein function was for proteases, accounting for only two of the 17 genes. Comparative analysis of the patient tumors that formed tumorgrafts to patient tumors that did not form tumorgrafts identified 491 genes that were up-regulated and 691 genes down-regulated 2-fold or greater in those tumors that formed tumorgrafts compared to non-tumorgraft forming tumors (Additional file [103]2: Table S2). The up-regulated genes in the tumors that formed tumorgrafts were enriched for cell cycle, cell signaling pathways and cytoskeleton remodeling (Figure [104]4A). One key canonical pathway map, enriched with up-regulated genes that were elevated in those tumors that formed tumorgrafts, is for WNT signaling (Figure [105]5); genes up-regulated in this pathway map included VEGF-A, Frizzled, β-Catenin, c-Myc, and FAK-1. The down-regulated genes in the tumorgraft-forming tumors were enriched for immune response pathways (8/10 pathways, Figure [106]4B). One key pathway map enriched with down-regulated genes is that of the role of integrins in NK cell cytotoxicity (Figure [107]6); note that in the context of this study, the upper Target Cell in this pathway corresponds to the relevant tumorgraft-forming cells. Figure 4. [108]Figure 4 [109]Open in a new tab Enrichment of canonical pathway maps in genes up and downregulated in tumorgraft forming tumors compared to non-tumorgraft forming tumors. Ontology enrichment analysis was performed using the MetaCore™ pathway analysis suite (GeneGo-Thomson Reuters, St. Joseph, MI). A) The top ten canonical pathway maps enriched in the genes up-regulated (n = 491) in the tumors that formed tumorgrafts included cell signalling and cell cycle-related pathways. B) Eight of top ten canonical pathway maps enriched in the genes down-regulated (n = 691) in the tumorgrafts that formed tumorgrafts were immune-related pathways. Figure 5. [110]Figure 5 [111]Open in a new tab WNT signalling pathway map enriched with genes up-regulated in tumorgraft-forming tumors. The MetaCore™ pathway analysis suite (GeneGo-Thomson Reuters, St. Joseph, MI) identified a pathway map for WNT signalling. VEGF-A, Frizzled, β-catenin, FAK-1, c-Myc, Laminin 1, XIAP, p21, MELP, and 4E-BP1 present in this pathway were all up-regulated in tumorgraft forming tumors. Figure 6. [112]Figure 6 [113]Open in a new tab Immune response pathway map enriched with genes down-regulated in the tumorgraft forming patient tumors. The MetaCore™ pathway analysis suite (GeneGo-Thomson Reuters, St. Joseph, MI) identified a pathway map for an immune response involving integrins in NK cell cytotoxicity that was enriched for several ligands down-regulated in tumors that did not form tumorgrafts (Target cell denotes tumor cells). Discussion The goal of this study was to create a panel of tumorgraft models, developed directly from patient tumor tissue from a wide range of heterogeneous tumor types of various stages. Prior investigations have focused on specific cancers, e.g. Secondary liver cancer [[114]15] lung [[115]24], pancreatic cancer [[116]16], pediatric osteosarcomas [[117]14], or pediatric rhabdomyosarcoma [[118]25], rather than a spectrum of malignancies. To examine the suitability of the panel of tumorgraft models generated in our lab as molecularly-relevant preclinical models, we directly compared a number of phenotypic and genotypic features of the tumorgrafts to the originating patient tumors from which they were derived. We have fully developed 49 models, spanning 18 different cancer types. A critical component of our development process is the ability of the model to re-establish following cryopreservation, alleviating the need for continual propagation, reported by some investigators [[119]16]. Cryopreservation significantly reduces the cost and necessary resources, extending the life of the tumor model, and minimizing genetic drift that could occur following long-term continuous in vivo propagation. Our overall tumor take rate is similar to previous reports [[120]3,[121]5]. A number of strategies have been tested to improve tumorgraft take rates, including the severe combined immunodeficient mouse strains CB17SC-M-F scid^−/−[[122]9,[123]25,[124]26] and CB17/Icr scid[[125]12], suppression of the immune system of recipient CBA/CaJ mice by thymectomy, whole body irradiation, administration of 1-β-D-arabinofuranosylcytosine [[126]14,[127]25], and use of orthotopic models [[128]27]. However, these strategies could increase the cost of model development with minimal improvement in tumorgraft development over the use of subcutaneous implantation into naive athymic nu/nu mice reported in this study and by other investigators [[129]2,[130]16]. It should be noted use of the athymic nu/nu mouse for use in the development of such tumorgraft does have limitations. The use of more immunodeficient mouse strains, e.g. NOD/scid or NOD/scid gamma (NSG), are necessary for successful engraftment of certain tumor types e.g. leukemia [[131]28,[132]29], may better resemble the human microenvironment, and could provide better consistency in growth rates following implantation of human tumor tissue into mice. However, in the context of developing a panel of tumorgraft models across a wide range of tumor types, the athymic nu/nu mouse has proven a robust model. The stage of the cancer at the time of acquisition has an impact on a successful tumorgraft development. Only 10% of Stage I and II tumors successfully formed tumorgrafts following implantation into mice; this is in contrast to the 43% take rate of tumor tissue from patients with Stage III or IV cancer at time of tissue acquisition which was successfully propagated in the murine host. Overall, ~75% of the 49 tumors that formed tumorgrafts were from patients that were diagnosed either Stage III or IV. In this study, breast cancers had a take rate of only 6% following subcutaneous implantation into mice. Dobrolecki et al [[133]27] reported success in developing breast cancer models by orthotopically implanting the breast cancer tissue into the mammary fat pads of mice. While more labor intensive than the subcutaneous model, an orthotopic approach should be considered if there is need for a tumorgraft model of a specific cancer that cannot be established using the subcutaneous model. An important conclusion from this study is that these subcutaneous (predominantly ectopic with the exception of melanomas) tumorgraft models recapitulate the expression profiles of the original patient’s tumor within its natural site of site of metastatic spread. This provides some evidence that refutes the lack of relevance of the subcutaneous model for drug discovery and development, which may rather reflect the use of genetically drifted in vitro cell lines in the classical xenograft context. The primary value for any tumor model is to recapitulate the in-vivo malignant tissue state, and accurately predict the efficacy of therapeutic agents in a timely and cost-effective manner. In contrast with commonly used animal models in which established human cancer cell lines are injected subcutaneously and fail to accurately predict human clinical trial results [[134]8], this human tumorgraft model system is more likely to be a successful surrogate preclinical model for several reasons. First, tumor bearing tissue fragments are transplanted containing not only malignant cells, but also supporting stromal tissue in an anatomically correct hierarchy. Second, the transferred tumors include cells post engraftment retaining their essential cell to cell interactions and relative state of differentiation as assessed by their histological appearance. Third, because of the high success rate in establishing not only 1^st, but also subsequent generation tumorgrafts, there is expansion of the number of tumorgrafts available for testing and the ability to verify retention of histological hallmarks in a serial fashion. Fourth, assessment of tumor growth is relatively straightforward given the subcutaneous location of the xenografts. Finally, since portions of tumors can be cryopreserved and successfully engrafted post thawing, testing of therapeutic agents, or combination therapies not originally considered, can be studied with this experimental and scalable in vivo model system. In the era of molecularly targeted oncology agents, the effective translation of results in the mouse model to the clinic requires the genome of the tumorgraft model to be high comparable to the original tumor [[135]30]. Significant changes in DNA fidelity or RNA expression levels significantly diminish the translational value of a model. Entire genome characterization of RNA expression rather, than a focused RT-PCR analysis of a select cohort of genes [[136]16], allows for the assessment of transcriptome wide alterations and associated biological systems that occur during the growth of the transplanted human tumors. The high Pearson correlations seen between the genomes of the 24 pairs of tumorgraft and originating patient tumors implies that the donor tumor genome is largely maintained in the tumorgraft model. Similarly, genotypic fidelity has also recently been reported in a panel of 25 human breast cancer tumorgrafts [[137]27] and on a smaller scale in secondary liver cancers using quantitative PCR analysis of 21 genes related to oncogenesis and cell cycle [[138]15]. The lower Pearson correlations observed in the gastrointestinal cancers could be an artifact of the tissue harvest and the high levels of digestive enzymes present in the patient tumors within the gastrointestinal tract. High levels of RNases present in pancreatic cancer are reported to complicate the extraction of high quality RNA [[139]31]. This finding highlights the importance of an optimized, rapid tissue collection protocol, especially for genomic profiling. It should be noted that poor RNA quality, for genomic profiling does not necessarily impact the ability of the tissue to form a viable tumorgraft. The high degree of clustering by matched tumor/tumorgraft pairs (79%) supports the high degree of similarity of the genomes of the patient tumor and resulting mouse tumorgraft even when the opportunity for RNA degradation may occur. Genomic stability is an important characteristic to consider prior to the long term use of a tumorgraft model. The high correlation in gene expression across four tumorgraft generations observed in this study and similar findings of others [[140]3,[141]16,[142]27,[143]30] and absence of drift in somatic mutations in known oncogenes, suggest that the genomes of tumorgraft models are stable. Histological integrity of patient-derived tumorgrafts has also been demonstrated for up to ten generations in immune-compromised mice [[144]4] and up to 30 passages in mice without significant changes in growth and morphological characteristics [[145]15]. One study reported tumorgraft models cultured in vivo for 10–12 generations before regenerating the models from earlier cryopreserved generations [[146]3]. The role of somatic oncogene mutations in tumorigenesis, pathogenesis and disease progression can have a profound influence on therapy [[147]32]; mutations in oncogenes RAS [[148]33] and B-RAF [[149]34] are commonly found in a variety of cancers. Therefore it is imperative that such mutations are maintained in the tumorgraft models. Although only a few oncogenic mutations were observed, those identified in the patient tumors were conserved in the matching tumorgraft, consistent with other reports [[150]16] and consistent with known mutations in specific tumor types [[151]20,[152]21]. There are a number of possible explanations for the enrichment of canonical immune pathways within those genes down-regulated in the tumorgrafts compared to patient tumors. First, a loss of tumor cell immune pathways and proteins not necessary for growth and development when transplanted from the immunecompetent patient into the immunocompromised host. Second, there is an intentional tumor tissue response to transplantation to down regulate those genes responsible for its immunogenic signature presented to the vestigial immune defenses of the immunocompromised mouse. A third explanation, and most likely, is that the gene expression changes are due to a loss of circulating human immune cells which infiltrate and support tumor development, when the tumor was transplanted from the human into the mouse as earlier reported by Neale et al [[153]17]. In those patient tumors that formed tumorgrafts, up-regulated genes show enrichment in pathways considered “usual suspects” in highly proliferative, late stage tumors; these canonical pathways include cell cycle, cytoskeletal remodeling, and those known drivers of tumorigenesis WNT and AKT signaling [[154]35-[155]37]. These same tumors were elevated for vascular endothelial growth factor (VEGF). This increased level of VEGF in the tumor could support the neovascularization of the tumor once implanted into the mouse; Nisolle et al [[156]38] reported that the survival and growth of human endometrium transplanted into nude mice was correlated with a high VEGF content. Transgenic expression of VEGF into human islet cells followed by transplant into the livers increased the rate of tissue revascularization, survival and function [[157]39]. A recent study demonstrated the pre-treatment of human ovarian tissue with VEGF and vitamin E prior to implantation into the back muscle of immunodeficient mice demonstrated enhanced growth, compared no pre-treatment of the ovarian tissue [[158]40]. The downregulation of immunological pathways in the tumors that formed tumors would suggest that loss of immune pathways is a mark of late stage tumors and provides a selective advantage supporting the establishment of the tumors in a foreign host or location. This downregulation of immunological pathways could be a reflection of the loss of expression signature of patient immune infiltrating cells [[159]17]. Of particular note is the downregulation of natural killer (NK) cell ligands, as identified as a pathway enriched for genes down-regulated in tumors that formed tumorgrafts compared to those tumors that failed to form tumorgrafts (Figure [160]4B). Although athymic nu/nu mice cannot generate mature T lymphocytes, they maintain the ability to mount a response to T-independent antigens [[161]41] and do produce NK cells as a component of their vestigial immune system. The downregulation of the NK ligands may therefore promote human tumor development in immune compromised mice. While there is little evidence demonstrating the interaction of human integrins and murine NK receptors, Kanwar et al reports interspecies cross-talk via the demonstration that the ectopic expression of human intercellular cell adhesion molecule-1 (ICAM-1) in mouse EL-4 tumors inhibited in vivo tumor growth [[162]42]. The mechanisms of NK cell recognition of self and immunosurveillance in both mouse and humans are complex and involve numerous receptors that recognize many cell surface ligands [[163]43]. Experimental evidence clearly supports the view that murine NK cells can inhibit human tumor cell development and metastasis both in vitro and in vivo[[164]44-[165]50]. This study reports the development of a series of tumorgraft models, representing a broad spectrum of cancer types, derived from molecularly characterized fresh patient tumor tissue. The genomes of these tumorgraft models exhibit a high degree of correlation to their originating patient tumors which is maintained over multiple generations of in vivo passages. The loss of immune function pathways and gain of cell proliferation pathways suggests a logical selective pressure for certain tumors to form tumorgrafts, which could be prescreened in future studies to identify those most likely to develop into useful preclinical models for translational studies. Conclusions The maintenance of the genome in the tumorgraft, from the human to murine host, supports the concept that tumorgraft models can provide an invaluable and robust tool in the evaluation of precision medicine methodologies and novel treatment strategies prior to use in clinical trials. An important point is that a subcutaneous primary patient tumorgraft models can recapitulate the expression profiles of the original “orthotopic” disease in the human and allow for reasonably high-throughput preclinical screens for drug efficacy. These models may therefore be useful in reducing the rate of attrition in drug development that has been in part been attributed to the use of the classical cell line xenografts. Competing interests CPW has disclosed a potential conflict of interest in relation to his invention of the XB-BIS informatics system used to analyze some of the data presented in this manuscript. He is also a consultant for the company TransMed Systems who acquired licensing rights to XB-BIS. Some third parties have acquired some of the described tumorgraft models for commercial projects, but these remain available for academic (noncommercial) purposes. Authors’ contributions DJM supervised the animal work, analyzed the data, interpreted the results, and wrote the manuscript. NRM interpreted the results and revised the manuscript. DMC, PJR, AI, and YN analyzed the data. DD and SBS carried out the animal work. EE, HF, SS, and REE carried out the laboratory experiments. JR, RS, and BN analyzed the histopathological images. CPW defined the research project, interpreted the results, and revised the manuscript. All authors read and approved the final manuscript. Supplementary Material Additional file 1 Table S1. The spreadsheet contains the results of the differential gene expression analysis comparing each patient’s tumor with the paired tumorgraft. This file contains two spreadsheet tabs. The first tab displays the 17 gene probesets up-regulated in tumorgrafts (Group 2, column C) relative to the originating patient tumors (Group 1, column B), and the second tab displays the 395 gene probes down-regulated. The official gene name, mean normalized intensity in Groups 1 and 2, the associated fold-change in gene expression (Group 2/Group 1), the calculated pairwise t-test statistic, the corresponding p-value, the False Discovery Rate, and the Affymetrix probeset identifiers are shown. (XLSX 59 kb) [166]Click here for file^ (59.9KB, xlsx) Additional file 2 Table S2. The spreadsheet contains the results of the differential gene expression analysis comparing the patient tumors that successfully developed tumorgrafts with those that did not. This file contains two spreadsheet tabs. The first sheet displays the 491 gene probesets up regulated in the tumors that successfully formed tumorgrafts (Group 2, column C) relative to those that did not (Group 1, column B) and the second sheet displays the 691 gene probes down regulated. The official gene name, mean normalized intensity in Groups 1 and 2, the associated fold-change in gene expression (Group 2/Group 1), the calculated student’s t-test statistic, the corresponding p-value, the False Discovery Rate, and the Affymetrix probeset identifiers are shown. (XLSX 137 kb) [167]Click here for file^ (137.2KB, xlsx) Contributor Information David J Monsma, Email: david.monsma@vai.org. Noel R Monks, Email: noel.monks@vai.org. David M Cherba, Email: david.cherba@vai.org. Dawna Dylewski, Email: dawna.dylewski@vai.org. Emily Eugster, Email: emily.eugster@vai.org. Hailey Jahn, Email: Hailey.Jahn@gmail.com. Sujata Srikanth, Email: sujatasrikanth@yahoo.com. Stephanie B Scott, Email: stephanie.scott@vai.org. Patrick J Richardson, Email: pjrichardson1@comcast.net. Robin E Everts, Email: reverts@sequenom.com. Aleksandr Ishkin, Email: aleksandr.ishkin@thomsonreuters.com. Yuri Nikolsky, Email: yuri.nikolsky@thomsonreuters.com. James H Resau, Email: james.resau@vai.org. Robert Sigler, Email: rsigler3667@yahoo.com. Brian J Nickoloff, Email: brian.nickoloff@vai.org. Craig P Webb, Email: craig.webb@vai.org. Acknowledgements