Implementing subtype-specific pre-clinical models of breast cancer to study pre-treatment aspirin effects

Backgorund

Prior data suggest pre‐diagnostic aspirin use impacts breast tumour biology and patient outcome. Here, we employed faithful surgical resection models of HER2+ and triple‐negative breast cancer (TNBC), to study outcome and response mechanisms across breast cancer subtypes.

Method

NOD/SCID mice were implanted with HER2+ MDA‐MB‐231/LN/2‐4/H2N, trastuzumab‐resistant HER2+ HCC1954 or a TNBC patient‐derived xenograft (PDX). A daily low‐dose aspirin regimen commenced until primary tumours reached ~250 mm³ and subsequently resected. MDA‐MB‐231/LN/2‐4/H2N mice were monitored for metastasis utilising imaging. To interrogate the survival benefit of pre‐treatment aspirin, 3 weeks post‐resection, HCC1954/TNBC animals received standard‐of‐care (SOC) chemotherapy for 6 weeks. Primary tumour response to aspirin was interrogated using immunohistochemistry.

Results

Aspirin delayed time to metastasis in MDA‐MB‐231/LN/2‐4/H2N xenografts and decreased growth of HER2⁺/TNBC primary tumours. Lymphangiogenic factors and lymph vessels number were decreased in HER2⁺ tumours. However, no survival benefit was seen in aspirin pre‐treated animals (HCC1954/TNBC) that further received adjuvant SOC, compared with animals treated with SOC alone. In an effort to study mechanisms responsible for the observed reduction in lymphangiogenesis in HER2⁺ BC we utilised an in vitro co‐culture system of HCC1954 tumour cells and mesenchymal stromal cells (MSC). Aspirin abrogated the secretion of VEGF‐C in MSCs and also decreased the lymph/angiogenic potential of the MSCs and HCC1954 by tubule formation assay. Furthermore, aspirin decreased the secretion of uPA in HCC1954 cells potentially diminishing its metastatic capability.

Conclusion

Our data employing clinically relevant models demonstrate that aspirin alters breast tumour biology. However, aspirin may not represent a robust chemo‐preventative agent in the HER2⁺ or TNBC setting.

Molecular Subtyping Combined with Biological Pathway Analyses to Study Regorafenib Response in Clinically Relevant Mouse Models of Colorectal Cancer

PURPOSE

Regorafenib (REG) is approved for the treatment of metastatic colorectal cancer, but has modest survival benefit and associated toxicities. Robust predictive/early response biomarkers to aid patient stratification are outstanding. We have exploited biological pathway analyses in a patient-derived xenograft (PDX) trial to study REG response mechanisms and elucidate putative biomarkers.

EXPERIMENTAL DESIGN

Molecularly subtyped PDXs were annotated for REG response. Subtyping was based on gene expression (CMS, consensus molecular subtype) and copy-number alteration (CNA). Baseline tumor vascularization, apoptosis, and proliferation signatures were studied to identify predictive biomarkers within subtypes. Phospho-proteomic analysis was used to identify novel classifiers. Supervised RNA sequencing analysis was performed on PDXs that progressed, or did not progress, following REG treatment.

RESULTS

Improved REG response was observed in CMS4, although intra-subtype response was variable. Tumor vascularity did not correlate with outcome. In CMS4 tumors, reduced proliferation and higher sensitivity to apoptosis at baseline correlated with response. Reverse phase protein array (RPPA) analysis revealed 4 phospho-proteomic clusters, one of which was enriched with non-progressor models. A classification decision tree trained on RPPA- and CMS-based assignments discriminated non-progressors from progressors with 92% overall accuracy (97% sensitivity, 67% specificity). Supervised RNA sequencing revealed that higher basal EPHA2 expression is associated with REG resistance.

CONCLUSIONS

Subtype classification systems represent canonical "termini a quo" (starting points) to support REG biomarker identification, and provide a platform to identify resistance mechanisms and novel contexts of vulnerability. Incorporating functional characterization of biological systems may optimize the biomarker identification process for multitargeted kinase inhibitors.

Conservation of copy number profiles during engraftment and passaging of patient-derived cancer xenografts

Patient-derived xenografts (PDXs) are resected human tumors engrafted into mice for preclinical studies and therapeutic testing. It has been proposed that the mouse host affects tumor evolution during PDX engraftment and propagation, affecting the accuracy of PDX modeling of human cancer. Here, we exhaustively analyze copy number alterations (CNAs) in 1,451 PDX and matched patient tumor (PT) samples from 509 PDX models. CNA inferences based on DNA sequencing and microarray data displayed substantially higher resolution and dynamic range than gene expression-based inferences, and they also showed strong CNA conservation from PTs through late-passage PDXs. CNA recurrence analysis of 130 colorectal and breast PT/PDX-early/PDX-late trios confirmed high-resolution CNA retention. We observed no significant enrichment of cancer-related genes in PDX-specific CNAs across models. Moreover, CNA differences between patient and PDX tumors were comparable to variations in multiregion samples within patients. Our study demonstrates the lack of systematic copy number evolution driven by the PDX mouse host.

Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models

Resistance to chemotherapy often results from dysfunctional apoptosis, however multiple proteins with overlapping functions regulate this pathway. We sought to determine whether an extensively validated, deterministic apoptosis systems model, 'DR_MOMP', could be used as a stratification tool for the apoptosis sensitiser and BCL-2 antagonist, ABT-199 in patient-derived xenograft (PDX) models of colorectal cancer (CRC). Through quantitative profiling of BCL-2 family proteins, we identified two PDX models which were predicted by DR_MOMP to be sufficiently sensitive to 5-fluorouracil (5-FU)-based chemotherapy (CRC0344), or less responsive to chemotherapy but sensitised by ABT-199 (CRC0076). Treatment with ABT-199 significantly improved responses of CRC0076 PDXs to 5-FU-based chemotherapy, but showed no sensitisation in CRC0344 PDXs, as predicted from systems modelling. ¹⁸F-Fluorodeoxyglucose positron emission tomography/computed tomography (¹⁸F-FDG-PET/CT) scans were performed to investigate possible early biomarkers of response. In CRC0076, a significant post-treatment decrease in mean standard uptake value was indeed evident only in the combination treatment group. Radiomic CT feature analysis of pre-treatment images in CRC0076 and CRC0344 PDXs identified features which could phenotypically discriminate between models, but were not predictive of treatment responses. Collectively our data indicate that systems modelling may identify metastatic (m)CRC patients benefitting from ABT-199, and that ¹⁸F-FDG-PET could independently support such predictions.

Abstract 1118: Absence of mouse-specific tumor evolution in patient-derived cancer xenografts

Patient-Derived Xenografts (PDXs) are preclinical models largely used to study tumor biology and drug response. Recent literature highlighted the possibility that growth of human tumors in a mouse microenvironment imposes a selection driving mouse-specific genetic evolution of PDXs, which may compromise their reliability as human cancer models. Conversely, independent studies observed a conservation of the genomic landscape during PDX engraftment and passaging.

We noticed that PDX genetic evolution was particularly evident in studies based on copy number aberration (CNA) inferred from gene expression data, while it was negligible when DNA-based CNA profiles were employed. Therefore, in a joint international effort of the EurOPDX and PDXNet consortia, we assembled a dataset of 37 hepatocellular and 54 gastric carcinoma tumor or PDX samples with matched RNA-based and DNA-based CNA profiles. We found that DNA-based CNA profiles invariably yield higher concordance between patient's tumor and derived PDXs than those inferred from RNA. RNA-based profiles displayed poor concordance with matched DNA-based profiles, and much lower resolution, so that they missed many focal copy number events detected by DNA-based methods. These results revealed that CNA measurements cannot be accurately estimated by expression data and that a systematic reassessment of CNA dynamics in PDXs based on DNA data is required.

To this aim, we generated CNA profiles by low-pass whole genome sequencing (WGS) of 87 colorectal and 43 breast cancer triplets, each composed of matched patient's tumor (PT) and PDX at early (PDX-early) and later (PDX-late) passage. In this way, for each tumor type, we generated three perfectly matched PT, PDX-early and PDX-late cohorts and performed CNA recurrence analysis by GISTIC in each cohort. The hypothesis was that if the mouse host induces a selective pressure capable of shaping the CNA landscape during PDX engraftment and propagation, GISTIC analysis would highlight systematic and progressive changes, from the PT to the PDX-early cohort, and then to the PDX-late cohort. Notably instead, the CNA profiles of the PT and PDX-early/late cohorts were virtually indistinguishable, with no progressive accumulation or loss of CNA during PDX passage and only minor changes not functionally related or associated to cancer-driver or actionable genes. These results were not consequence of insufficient capture of the CNA repertoire, since the GISTIC profiles recapitulated those generated by TCGA for colorectal and breast cancer. In summary, our analyses highlighted that while RNA-based CNA inferences have inadequate resolution and accuracy to study genomic evolution in PDXs, DNA-based CNA profiles confirm retention of CNAs in PTs and PDXs, excluding a systematic mouse driven selection via copy number changes. Ultimately, these results support the robustness of PDXs as preclinical models for predicting drug response.

Citation Format: Jessica Giordano, Xing Yi Woo, Anuj Srivastava, Zi-Ming Zhao, Michael W. Lloyd, Roebi de Bruijn, Yun-Suhk Suh, Francesco Galimi, Andrea Bertotti, Adam Lafferty, Alice C. O'Farrell, Elodie Modave, Diether Lambrechts, Petra ter Brugge, Violeta Serra, Elisabetta Marangoni, Rania El Botty, Jong-Il Kim, Han-Kwang Yang, Charles Lee, Dennis A. Dean, Brandi Davis-Dusenbery, Yvonne A. Evrard, James H. Doroshow, Alana L. Welm, Bryan E. Welm, Michael T. Lewis, Bingliang Fang, Jack Roth, Funda Meric-Bernstam, Meenhard Herlyn, Michael Davies, Li Ding, Shunqiang Li, Ramaswamy Govindan, Jeffrey A. Moscow, Carol J. Bult, Claudio Isella, Livio Trusolino, Annette T. Byrne, Jos Jonkers, Jeffrey H. Chuang, Enzo Medico, EurOPDX consortium & PDXNET consortium. Absence of mouse-specific tumor evolution in patient-derived cancer xenografts [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1118.

Implementing Reverse Phase Protein Array Profiling as a Sensitive Method for the Early Pre-Clinical Detection of Off-Target Toxicities Associated with Sunitinib Malate

PURPOSE

The tyrosine kinase inhibitor (TKI) sunitinib is a multi-targeted agent approved across multiple cancer indications. Nevertheless, since approval, data has emerged to describe a worrisome side effect profile including hypertension, hand-foot syndrome, fatigue, diarrhea, mucositis, proteinuria, and (rarely) congestive heart failure. It has been hypothesized that the observed multi-parameter toxicity profile is related to "on-target" kinase inhibition in "off-target" tissues.

EXPERIMENTAL DESIGN

To interrogate off-target effects in pre-clinical studies, a reverse phase protein array (RPPA) approach is employed. Mice are treated with sunitinib (40 mg kg^-1 ) for 4 weeks, following which critical organs are removed. The Zeptosens RPPA platform is employed for protein expression analysis.

RESULTS

Differentially expressed proteins associated with damage and/or stress are found in the majority of organs from treated animals. Proteins differentially expressed in the heart are associated with myocardial hypertrophy, ischaemia/reperfusion, and hypoxia. However, hypertrophy is not evidenced on histology. Mild proteinuria is observed; however, no changes in renal glomerular structure are visible via electron microscopy. In skin, proteins associated with cutaneous inflammation, keratinocyte hyper-proliferation, and increased inflammatory response are differentially expressed.

CONCLUSIONS AND CLINICAL RELEVANCE

It is posited that pre-clinical implementation of a combined histopathological/RPPA approach provides a sensitive method to mechanistically elucidate the early manifestation of TKI on-target/organ off-target toxicities.

A Type 1 Diabetes Genetic Risk Score Predicts Progression of Islet Autoimmunity and Development of Type 1 Diabetes in Individuals at Risk

OBJECTIVE

We tested the ability of a type 1 diabetes (T1D) genetic risk score (GRS) to predict progression of islet autoimmunity and T1D in at-risk individuals.

RESEARCH DESIGN AND METHODS

We studied the 1,244 TrialNet Pathway to Prevention study participants (T1D patients' relatives without diabetes and with one or more positive autoantibodies) who were genotyped with Illumina ImmunoChip (median [range] age at initial autoantibody determination 11.1 years [1.2-51.8], 48% male, 80.5% non-Hispanic white, median follow-up 5.4 years). Of 291 participants with a single positive autoantibody at screening, 157 converted to multiple autoantibody positivity and 55 developed diabetes. Of 953 participants with multiple positive autoantibodies at screening, 419 developed diabetes. We calculated the T1D GRS from 30 T1D-associated single nucleotide polymorphisms. We used multivariable Cox regression models, time-dependent receiver operating characteristic curves, and area under the curve (AUC) measures to evaluate prognostic utility of T1D GRS, age, sex, Diabetes Prevention Trial-Type 1 (DPT-1) Risk Score, positive autoantibody number or type, HLA DR3/DR4-DQ8 status, and race/ethnicity. We used recursive partitioning analyses to identify cut points in continuous variables.

RESULTS

Higher T1D GRS significantly increased the rate of progression to T1D adjusting for DPT-1 Risk Score, age, number of positive autoantibodies, sex, and ethnicity (hazard ratio [HR] 1.29 for a 0.05 increase, 95% CI 1.06-1.6; P = 0.011). Progression to T1D was best predicted by a combined model with GRS, number of positive autoantibodies, DPT-1 Risk Score, and age (7-year time-integrated AUC = 0.79, 5-year AUC = 0.73). Higher GRS was significantly associated with increased progression rate from single to multiple positive autoantibodies after adjusting for age, autoantibody type, ethnicity, and sex (HR 2.27 for GRS >0.295, 95% CI 1.47-3.51; P = 0.0002).

CONCLUSIONS

The T1D GRS independently predicts progression to T1D and improves prediction along T1D stages in autoantibody-positive relatives.

Abstract 3010: Proteomic analysis of off-target toxicities following treatment with the tyrosine kinase inhibitor sunitinib

Background: Sunitinib is a multi-targeted agent approved for treatment of a number of cancers. However, since approval, data has continually emerged relating to toxicity profiles including hypertension, hand-foot syndrome and fatigue. Underlying mechanisms are unresolved. It has been hypothesised that the multi-parameter toxicity profile is related to ‘on-target’ kinase inhibition in ‘off-target’ tissues.

Methods: To interrogate off-target effects, the Zeptosens Reverse Phase protein Array (RPA) platform was used to assess tissues obtained from mice treated with sunitinib (40 mg/kg) for 4 weeks. Additional tissue was collected for histological analyses and pathophysiologic changes assessed.

Results: RPA data analyses on all organs collected (heart, kidney, liver, brain, intestine and skin) revealed the presence of differentially expressed proteins associated with damage and/or stress. Of note, Cyclin D1, MEK and phosphorylated-p21(CIP/WAF1) expression increased by 1.8-fold (p < 0.05), 2.2-fold (p < 0.1) and 1.7-fold (p < 0.1) respectively in the kidneys of sunitinib treated mice. These proteins are associated with kidney damage after chemotherapy, and induction of genes associated with renal immune response, inflammation and apoptosis. Mild proteinuria was observed in sunitinib treated animals, however no gross change in renal glomerular ultrastructure was visible via electron microscopy. Phosphorylated-FGFR1, cleaved Caspase-7, and phosphorylated-cMYC expression decreased by 0.45-fold (p < 0.05), 0.4-fold (p < 0.1) and 0.8-fold (p < 0.1) respectively in the skin of sunitinib treated mice. Down-regulation of these proteins has been associated with a perturbation of cutaneous wound healing (p-FGFR1/cleaved Caspase-7). Mice treated with high dose sunitinib (80 mg/kg) showed overt signs of skin toxicity. Histopathologic studies of the thyroid gland revealed decreased thyroid epithelial cell height (p < 0.05). Nevertheless, serum levels of triiodothyronine and thyroxine remained unchanged. An orthogonal validation study (Western blotting) is ongoing based on gene ontology, pathway analysis and observed fold changes in critical signalling pathway proteins involved in cellular stress.

Conclusions: Implementation of a combined histopathologic/ RPA approach may provide a sensitive method to mechanistically elucidate the off-target toxicity sequelae associated with TKIs approved in the oncology setting.

Work was supported by an EU funded Industry Academia Pathways and Partnerships Marie Curie Award (AngioTox) Grant Number 251528. RPAs were performed at Bayer Technology Services, Leverkusen, Germany.

Citation Format: Alice C. O’Farrell, Adam Lafferty, Ian Miller, Rhys Evans, Maurice Cary, David Murray, Marina Alamanou, Girish Mallya Udupi, Liam Shiels, William M. Gallagher, Annette T. Byrne. Proteomic analysis of off-target toxicities following treatment with the tyrosine kinase inhibitor sunitinib. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3010.

Relative abundance of growth hormone receptor isoforms in rhesus monkey tissues and human transformed lymphocytes

The growth hormone receptor (GHR) is expressed as one active, full-sequence isoform and one truncated, inactive one that lacks the intracellular signaling domain. The aim of this study was to investigate the variation in the tissue expression of the full and truncated mRNA and protein. Epstein-Barr virus-transformed human B lymphocyte lines were established from 9 normal individuals with a height standard deviation score (SDS) of - 0.1 +/- 1.1 (mean +/- SD). Tissues were also collected from 3 Rhesus monkeys, whose GHR has 94.1 % homology with the human molecule. mRNA quantitation was determined by Real Time Quantitative PCR. Growth hormone receptor expression in transformed lymphocytes was also studied by fluorescence-activated cell sorter analysis. Both isoforms were expressed in transformed lymphocytes, but individual variation in the relative mRNA expression was small (truncated isoform percentage of total receptor mRNA: 17.1 +/- 4.4, mean +/- SD). There was no correlation between donors' height SDS and the expression of either isoform or the ratio between them. Protein expression by FACS analysis showed wider variation among the subjects; however, the relative ratio was similar in all the subjects. In monkey tissues, the truncated receptor showed a tissue-specific distribution. In conclusion, the expression of both isoforms in transformed lymphocytes from normal subjects shows small differences at the RNA or protein levels, and does not correlate with height SDS. Growth hormone splice isoforms show tissue specificity, suggesting local regulation of splicing. Tissues with relatively high expression of the truncated isoform are likely to be more resistant to the effects of GH due to the dominant negative effect of this isoform. In addition, the differential tissue expression might influence the levels of growth hormone binding protein in the immediate milieu of each tissue.

Adrenocortical-pituitary hybrid tumor causing Cushing's syndrome

We describe the first case of an adrenocortical-pituitary hybrid tumor causing Cushing's syndrome in a 17-yr-old boy. Adrenal vein sampling confirmed elevated secretion of both cortisol and ACTH precursors from a right adrenal mass, whereas pituitary ACTH levels, as determined by bilateral inferior petrosal sinus samples (IPSS), were unresponsive to CRH and equal to peripheral levels. There was no biochemical or histological evidence for a pheochromocytoma, but, rather, the tumor demonstrated lipid-rich clear cells characteristic of an adrenocortical adenoma. Immunohistochemical analysis revealed ACTH immunoreactivity and synaptophysin proteins in the tumor. Isolation of tumor cells by the novel technique of laser capture microdissection and subsequent RT-PCR showed expression of POMC messenger ribonucleic acid and cytochrome p450 enzyme messenger ribonucleic acid within the same cells. Finally, ultrastructural analysis provided ultimate proof for adrenocortical-pituitary hybrid cells exhibiting the characteristic vesicular mitochondria and abundant smooth endoplasmic reticulum of steroid cells and the typical secretory granules of corticotrophs within the cytoplasm of the same cells. The adrenocortical tumor expressed the pituitary transcription factor pituitary homeobox factor 1 and the steroidogenic factor 1. The intermingling of the centrally located ectodermally derived pituitary tissue with the mesodermally derived adrenocortical tissue in this adenoma suggests a hitherto unrecognized genetic and phenotypic plasticity within the hypothalamic-pituitary-adrenal axis.

Anisomastia associated with interstitial duplication of chromosome 16, mental retardation, obesity, dysmorphic facies, and digital anomalies: molecular mapping of a new syndrome by fluorescent in situ hybridization and microsatellites to 16q13 (D16S419-D16S503)

Anisomastia is a common problem among developing adolescent girls. We recently evaluated a 22-yr-old female patient who had severe anisomastia (which had been repaired by surgery), associated with moderate to severe mental retardation, a stocky body habitus with mild obesity, dysmorphic facies (prominent, upslanting palpebral fissures, beaked nose, and a prominent philtrum), webbed neck, low hairline, and severe bilateral clinodactyly of the third, fourth, and fifth fingers with acral (but not large joint) flexion contractures. A peripheral blood high resolution karyotype revealed additional chromosomal material within the long arm of chromosome 16. Densitometric analysis of amplified polymorphic sequence-tagged sites (STS) mapping to 16q suggested that the duplication is defined by the noninvolved markers D16S419 [16q12-cen, 66 centimorgan (cM) from 16p terminus] and D16S421 (16q13-q21, 84.4 cM), encompassing a maximum of 18.4 cM of genetic distance. The STS analysis showed that the duplication was on the maternally derived chromosome 16, resulting in two maternal (and one paternal) copies of that region of chromosome 16. The location was further confirmed by bacterial artificial chromosomes (BACs) that were obtained from a commercially available library, labeled, and used for fluorescence in situ hybridization. The BACs containing STSs D16S408, D16S3137, and D16S3032 (markers that correspond to 16q13) showed two regions of hybridization, indicating that these sites were duplicated, whereas a BAC containing the STS D16S512 (which corresponds to 16q21-q22) revealed one hybridization signal per 16q, indicating that the corresponding region was not involved in the duplication. The distance between the probe signals suggested a tandem duplication. We conclude that even though trisomy 16 is the most common autosomal trisomy in spontaneous abortions, few patients with unbalanced chromosome 16 abnormalities survive to adulthood; in this report we describe one such patient with an interstitial chromosome 16 duplication (at 16q13), who had a specific phenotype associated with abnormal breast size. There are clinical similarities between this patient and patients with other 16q abnormalities, although the breast findings were unique. Molecular cytogenetics, including fluorescence in situ hybridization and densitometric analysis of amplified STSs, provided useful tools for the precise mapping of the syndrome to 16q13, where the gene(s) responsible for this phenotype might be localized.

Rapid diagnosis of germline p53 mutation using the enzyme mismatch cleavage method

The p53 tumor suppressor gene is the most commonly altered gene in human cancers. Germline mutations in p53 are the genetic alteration underlying predisposition to multiple cancers in Li-Fraumeni syndrome and Li-Fraumeni-like syndrome. We describe a patient who presented with developed adrenocortical carcinoma at age 19 months and a cerebral primitive neuroectodermal tumor at age 5 years. The patient did not have a family history of cancer. We used the enzyme mismatch cleavage (EMC) method to screen for mutations in the p53 gene and found a germline mutation in exon 7 (codon 248). Loss of heterozygosity analysis in one tumor revealed loss of the wild-type p53 allele. In our report we demonstrate the EMC method to be a rapid and sensitive method for mutation detection.