Abbas Shirinifard

Targeting the DNA Repair Pathway in Ewing Sarcoma

Ewing sarcoma (EWS) is a tumor of the bone and soft tissue that primarily affects adolescents and young adults. With current therapies, 70% of patients with localized disease survive, but patients with metastatic or recurrent disease have a poor outcome. We found that EWS cell lines are defective in DNA break repair and are sensitive to PARP inhibitors (PARPis). PARPi-induced cytotoxicity in EWS cells was 10- to 1,000-fold higher after administration of the DNA-damaging agents irinotecan or temozolomide. We developed an orthotopic EWS mouse model and performed pharmacokinetic and pharmacodynamic studies using three different PARPis that are in clinical development for pediatric cancer. Irinotecan administered on a low-dose, protracted schedule previously optimized for pediatric patients was an effective DNA-damaging agent when combined with PARPis; it was also better tolerated than combinations with temozolomide. Combining PARPis with irinotecan and temozolomide gave complete and durable responses in more than 80% of the mice.

The Dynamic Epigenetic Landscape of the Retina During Development, Reprogramming, and Tumorigenesis

SUMMARY In the developing retina, multipotent neural progenitors undergo unidirectional differentiation in a precise spatiotemporal order. Here we profile the epigenetic and transcriptional changes that occur during retinogenesis in mice and humans. Although some progenitor genes and cell cycle genes were epigenetically silenced during retinogenesis, the most dramatic change was derepression of cell type–specific differentiation programs. We identified developmental stage–specific super-enhancers and showed that most epigenetic changes are conserved in humans and mice. To determine how the epigenome changes during tumorigenesis and reprogramming, we performed integrated epigenetic analysis of murine and human retinoblastomas and induced pluripotent stem cells (iPSCs) derived from murine rod photoreceptors. The retinoblastoma epigenome mapped to the developmental stage when retinal progenitors switch from neurogenic to a terminal patterns of cell division. The epigenome of retinoblastomas was more similar to that of normal retina than was that of retina-derived iPSCs, and we identified retina-specific epigenetic memory.

Emergent Stratification in Solid Tumors Selects for Reduced Cohesion of Tumor Cells: A Multi-Cell, Virtual-Tissue Model of Tumor Evolution Using CompuCell3D.

Tumor cells and structure both evolve due to heritable variation of cell behaviors and selection over periods of weeks to years (somatic evolution). Micro-environmental factors exert selection pressures on tumor-cell behaviors, which influence both the rate and direction of evolution of specific behaviors, especially the development of tumor-cell aggression and resistance to chemotherapies. In this paper, we present, step-by-step, the development of a multi-cell, virtual-tissue model of tumor somatic evolution, simulated using the open-source CompuCell3D modeling environment. Our model includes essential cell behaviors, microenvironmental components and their interactions. Our model provides a platform for exploring selection pressures leading to the evolution of tumor-cell aggression, showing that emergent stratification into regions with different cell survival rates drives the evolution of less cohesive cells with lower levels of cadherins and higher levels of integrins. Such reduced cohesivity is a key hallmark in the progression of many types of solid tumors.

Development and characterization of a human orthotopic neuroblastoma xenograft

Neuroblastoma is a pediatric cancer of the developing sympathoadrenal lineage. The tumors are known to develop from the adrenal gland or paraspinal ganglia and have molecular and cellular features of sympathetic neurons such as dense core vesicles and catecholamine production. Here we present the detailed molecular, cellular, genetic and epigenetic characterization of an orthotopic xenograft derived from a high-risk stage 4 neuroblastoma patient. Overall, the xenografted tumor retained the high risk features of the primary tumor and showed aggressive growth and metastasis in the mouse. Also, the genome was preserved with no additional copy number variations, structural variations or aneuploidy. There were 13 missense mutations identified in the xenograft that were not present in the patient’s primary tumor and there were no new nonsense mutations. None of the missense mutations acquired in the xenograft were in known cancer genes. We also demonstrate the feasibility of using the orthotopic neuroblastoma xenograft to test standard of care chemotherapy and molecular targeted therapeutics. Finally, we optimized a new approach to produce primary cultures of the neuroblastoma xenografts for high-throughput drug screening which can be used to test new combinations of therapeutic agents for neuroblastoma.

3D quantitative analyses of angiogenic sprout growth dynamics

Background: Zebrafish intersegmental vessel (ISV) growth is widely used to study angiogenesis and to screen drugs and toxins that perturb angiogenesis. Most current ISV growth assays observe the presence or absence of ISVs or perturbation of ISV morphology but do not measure growth dynamics. We have developed a four‐dimensional (4D, space plus time) quantitative analysis of angiogenic sprout growth dynamics for characterization of both normal and perturbed growth. Results: We tracked the positions of the ISV base and tip for each ISV sprout in 4D. Despite immobilization, zebrafish embryos translocated globally and non‐uniformly during development. We used displacement of the ISV base and the angle between the ISV and the dorsal aorta to correct for displacement and rotation during development. From corrected tip cell coordinates, we computed average ISV trajectories. We fitted a quadratic curve to the average ISV trajectories to produce a canonical ISV trajectory for each experimental group, arsenic treated and untreated. From the canonical ISV trajectories, we computed curvature, average directed migration speed and directionality. Canonical trajectories from treated (arsenic exposed) and untreated groups differed in curvature, average directed migration speed and angle between the ISV and dorsal aorta. Conclusions: 4D analysis of angiogenic sprout growth dynamics: (1) Allows quantitative assessment of ISV growth dynamics and perturbation, and (2) provides critical inputs for computational models of angiogenesis. Developmental Dynamics 242:508–516, 2013.

Retinal Cell Type DNA Methylation and Histone Modifications Predict Reprogramming Efficiency and Retinogenesis in 3D Organoid Cultures

SUMMARY Diverse cell types can be reprogrammed into pluripotent stem cells by ectopic expression of Oct4 (Pou5f1), Klf4, Sox3, and Myc. Many of these induced pluripotent stem cells (iPSCs) retain memory, in terms of DNA methylation and histone modifications (epigenetic memory), of their cellular origins, and this may bias subsequent differentiation. Neurons are difficult to reprogram, and there has not been a systematic side-by-side characterization of reprogramming efficiency or epigenetic memory across different neuronal subtypes. Here, we compare reprogramming efficiency of five different retinal cell types at two different stages of development. Retinal differentiation from each iPSC line was measured using a quantitative standardized scoring system called STEM-RET and compared to the epigenetic memory. Neurons with the lowest reprogramming efficiency produced iPSC lines with the best retinal differentiation and were more likely to retain epigenetic memory of their cellular origins. In addition, we identified biomarkers of iPSCs that are predictive of retinal differentiation.

Detection of Phenotypic Alterations Using High-Content Analysis of Whole-Slide Images

Tumors exhibit spatial heterogeneity, as manifested in immunohistochemistry (IHC) staining patterns. Current IHC quantification methods lose information by reducing this heterogeneity in each whole-slide image (WSI) or in selective fields of view to a single staining index. The aim of this study was to investigate the sensitivity of an IHC quantification method that uses this heterogeneity to reliably compare IHC staining patterns. We virtually partitioned WSIs by a grid of square tiles, and computed the staining index distributions to quantify heterogeneities. We used samples from these distributions as inputs to non-parametric statistical comparisons. We applied our grid method to fixed tumor samples from 26 tumors obtained from a double-blind preclinical study of a patient-derived orthotopic xenograft model of pediatric neuroblastoma in CD1 nude mice. We compared the results of our grid method to the results based on whole-slide indices, the current practice. We show that our grid method reliably detects phenotypic alterations that other tests based on whole-slide indices fail to detect. Based on robustness and increased sensitivity of statistical inference, we conclude that our method of whole-slide grid quantification is superior to existing whole-slide quantification techniques.

The Nucleome of Developing Murine Rod Photoreceptors

The nuclei of rod photoreceptors in mice and other nocturnal species have an unusual inverted chromatin structure: the heterochromatin is centrally located to help focus light and improve photosensitivity. To better understand this unique nuclear organization, we performed ultra-deep Hi-C analysis on murine retina at 3 stages of development and on purified rod photoreceptors. Predicted looping interactions from the Hi-C data were validated with fluorescence in situ hybridization (FISH). We discovered that a subset of retinal genes that are important for retinal development, cancer, and stress response are localized to the facultative heterochromatin domain. We also used machine learning to develop an algorithm based on our chromatin Hidden Markov Modeling (chromHMM) of retinal development to predict heterochromatin domains and study their dynamics during retinogenesis. FISH data for 264 genomic loci were used to train and validate the algorithm. The integrated data were then used to identify a developmental stage– and cell type-specific core regulatory circuit super-enhancer (CRC-SE) upstream of the Vsx2 gene, which is required for bipolar neuron expression. Deletion of the Vsx2 CRC-SE in mice led to the loss of bipolar neurons in the retina.

Measuring Absolute Blood Perfusion in Mice Using Dynamic Contrast-Enhanced Ultrasound

We investigated the feasibility of estimating absolute tissue blood perfusion using dynamic contrast-enhanced ultrasound (CEUS) imaging in mice. We developed a novel method of microbubble administration and a model-free approach to estimate absolute kidney perfusion, and explored the kidney as a reference organ to estimate absolute perfusion of a neuroblastoma tumor. We performed CEUS on the kidneys of CD1 nude mice using the VisualSonics VEVO 2100 imaging system. We estimated individual kidney blood perfusion using the burst-replenishment (BR) technique. We repeated the kidney imaging on the mice after a week. We performed CEUS imaging of a neuroblastoma mouse xenograft tumor along with its right kidney using two sets of microbubble administration parameters to estimate absolute tumor blood perfusion. We performed statistical tests at a significance level of 0.05. Our estimated absolute kidney perfusion (425xa0±xa0123xa0mL/min/100xa0g) was within the range of previously reported values. There was no statistical difference between the estimated absolute kidney blood perfusions from the 2xa0wk of imaging (paired t-test, pxa0=xa00.09). We estimated the absolute blood perfusion in the neuroblastoma tumor to be 16.49 and 16.9xa0mL/min/100xa0g for the two sets of microbubble administration parameters (Wilcoxon rank-sum test, pxa0=xa00.6). We have established the kidney as a reliable reference organ in which to estimate absolute perfusion of other tissues. Using a neuroblastoma tumor, we have determined the feasibility of estimating absolute blood perfusion in tissues using contrast-enhanced ultrasound imaging.

Abstract 2708: Development of an individualized 3D transport model of topotecan for a patient-derived orthotopic xenograft model of pediatric neuroblastoma

Resistance to chemotherapeutics and targeted therapies in pediatric solid tumors including neuroblastoma is a common cause of poor clinical outcome. These failures in part stem from shortcomings in understanding inter- and intra-tumor heterogeneities of drug penetration due to heterogeneities in blood perfusion. Herein we propose to develop an individualized 3D transport model of topotecan (TPT) for a patient-derived orthotopic xenograft model of pediatric NB5 neuroblastoma to account for inter- and intra-tumor heterogeneities in blood perfusion. The transport model uses a 3D reaction-diffusion equation to simulate diffusion of TPT from blood vessels into the tumor tissue and its flux in and out of intracellular space. Our transport model takes three types of inputs to predict TPT exposure maps defined over the volume of an individual tumor: a) plasma concentration-time profiles from an individualized physiologically-based pharmacokinetic (PBPK) model of TPT (separate cohort), b) 3D blood perfusion map of the individual tumor from contrast enhanced ultrasound (CEUS) using VisualSonics VEVO 2100 imaging system, and c) in vitro TPT cellular uptake and efflux kinetics from two-photon imaging. We use in vitro pharmacodynamics (PD) experiments with NB5 cells exposed to TPT to derive probabilistic PD-rules for drug effects (e.g., γ-H2AX response). Based on these rules and the exposure maps, we then compute probabilities of effects for the entire tumor volume. We will validate the predicted drug effect maps by comparing them to the observed effects measured by immunohistochemistry marker for γ-H2AX from the same tumor (location matched) using spatial correlation techniques. Citation Format: Abbas Shirinifard, Suresh Thiagarajan, Yogesh T. Patel, Abigail D. Davis, Megan O. Jacus, Stacy L. Throm, Jessica Roberts, Vinay Daryani, Clinton F. Stewart, Andras Sablauer. Development of an individualized 3D transport model of topotecan for a patient-derived orthotopic xenograft model of pediatric neuroblastoma. [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 2708.