Junctional Adhesion Molecule (JAM)-C recruitment of Pard3 and drebrin to cell contacts initiates neuron-glia recognition and layer-specific cell sorting in developing cerebella.

Sorting maturing neurons into distinct layers is critical for brain development, with disruptions leading to neurological disorders and pediatric cancers. Lamination coordinates where, when, and how cells interact, facilitating events that direct migrating neurons to their destined positions within emerging neural networks and control the wiring of connections in functional circuits. While the role of adhesion molecule expression and presentation in driving adhesive recognition during neuronal migration along glial fibers is recognized, the mechanisms by which the spatial arrangement of these molecules on the cell surface dictates adhesive specificity and translates contact-based external cues into intracellular responses like polarization and cytoskeletal organization remain largely unexplored. We used the cerebellar granule neuron (CGN) system to demonstrate that JAM-C receptor cis-binding on the same cell and trans-binding to neighboring cells controls the recruitment of the Pard3 polarity protein and drebrin microtubule-actin crosslinker at CGN to glial adhesion sites, complementing previous studies that showed Pard3 controls JAM-C exocytic surface presentation. Leveraging advanced imaging techniques, specific probes for cell recognition, and analytical methods to dissect adhesion dynamics, our findings reveal: 1) JAM-C cis or trans mutants result in reduced adhesion formation between CGNs and cerebellar glia, 2) these mutants exhibit delayed recruitment of Pard3 at the adhesion sites, and 3) CGNs with JAM-C mutations experience postponed sorting and entry into the cerebellar molecular layer (ML). By developing a conditional system to image adhesion components from two different cells simultaneously, we made it possible to investigate the dynamics of cell recognition on both sides of neuron-glial contacts and the subsequent recruitment of proteins required for CGN migration. This system and an approach that calculates local correlation based on convolution kernels at the cell adhesions site revealed that CGN to CGN JAM recognition preferentially recruits higher levels of Pard3 and drebrin than CGN to glia JAM recognition. The long latency time of CGNs in the inner external germinal layer (EGL) can be attributed to the combined strength of CGN-CGN contacts and the less efficient Pard3 recruitment by CGN-BG contacts, acting as gatekeepers to ML entry. As CGNs eventually transition to glia binding for radial migration, our research demonstrates that establishing permissive JAM-recognition sites on glia via cis and trans interactions of CGN JAM-C serves as a critical temporal checkpoint for sorting at the EGL to ML boundary. This mechanism integrates intrinsic and extrinsic cellular signals, facilitating heterotypic cell sorting into the ML and dictating the precise spatial organization within the cerebellar architecture.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

A fluorogenic complementation tool kit for interrogating lipid droplet-organelle interaction.

Contact sites between lipid droplets and other organelles are essential for cellular lipid and energy homeostasis. Detection of these contact sites at nanometer scale over time in living cells is challenging. Here, we developed a tool kit for detecting contact sites based on Fluorogen-Activated Bimolecular complementation at CONtact sites, FABCON, using a reversible, low affinity split fluorescent protein, splitFAST. FABCON labels contact sites with minimal perturbation to organelle interaction. Via FABCON, we quantitatively demonstrated that endoplasmic reticulum (ER)- and mitochondria (mito)-lipid droplet contact sites are dynamic foci in distinct metabolic conditions, such as during lipid droplet biogenesis and consumption. An automated analysis pipeline further classified individual contact sites into distinct subgroups based on size, likely reflecting differential regulation and function. Moreover, FABCON is generalizable to visualize a repertoire of organelle contact sites including ER-mito. Altogether, FABCON reveals insights into the dynamic regulation of lipid droplet-organelle contact sites and generates new hypotheses for further mechanistical interrogation during metabolic switch.

An adaptive stress response that confers cellular resilience to decreased ubiquitination.

Ubiquitination is a post-translational modification initiated by the E1 enzyme UBA1, which transfers ubiquitin to ~35 E2 ubiquitin-conjugating enzymes. While UBA1 loss is cell lethal, it remains unknown how partial reduction in UBA1 activity is endured. Here, we utilize deep-coverage mass spectrometry to define the E1-E2 interactome and to determine the proteins that are modulated by knockdown of UBA1 and of each E2 in human cells. These analyses define the UBA1/E2-sensitive proteome and the E2 specificity in protein modulation. Interestingly, profound adaptations in peroxisomes and other organelles are triggered by decreased ubiquitination. While the cargo receptor PEX5 depends on its mono-ubiquitination for binding to peroxisomal proteins and importing them into peroxisomes, we find that UBA1/E2 knockdown induces the compensatory upregulation of other PEX proteins necessary for PEX5 docking to the peroxisomal membrane. Altogether, this study defines a homeostatic mechanism that sustains peroxisomal protein import in cells with decreased ubiquitination capacity.

Latent Epigenetic Programs in Müller Glia Contribute to Stress, Injury, and Disease Response in the Retina.

Previous studies have demonstrated the dynamic changes in chromatin structure during retinal development that correlate with changes in gene expression. However, a major limitation of those prior studies was the lack of cellular resolution. Here, we integrate single-cell (sc) RNA-seq and scATAC-seq with bulk retinal data sets to identify cell type-specific changes in the chromatin structure during development. Although most genes' promoter activity is strongly correlated with chromatin accessibility, we discovered several hundred genes that were transcriptionally silent but had accessible chromatin at their promoters. Most of those silent/accessible gene promoters were in the Müller glial cells. The Müller cells are radial glia of the retina and perform a variety of essential functions to maintain retinal homeostasis and respond to stress, injury, or disease. The silent/accessible genes in Müller glia are enriched in pathways related to inflammation, angiogenesis, and other types of cell-cell signaling and were rapidly activated when we tested 15 different physiologically relevant conditions to mimic retinal stress, injury, or disease in human and murine retinae. We refer to these as "pliancy genes" because they allow the Müller glia to rapidly change their gene expression and cellular state in response to different types of retinal insults. The Müller glial cell pliancy program is established during development, and we demonstrate that pliancy genes are necessary and sufficient for regulating inflammation in the murine retina in vivo. In zebrafish, Müller glia can de-differentiate and form retinal progenitor cells that replace lost neurons. The pro-inflammatory pliancy gene cascade is not activated in zebrafish Müller glia following injury, and we propose a model in which species-specific pliancy programs underly the differential response to retinal damage in species that can regenerate retinal neurons (zebrafish) versus those that cannot (humans and mice).

Platelet-derived chemokines promote skeletal muscle regeneration by guiding neutrophil recruitment to injured muscles.

Skeletal muscle regeneration involves coordinated interactions between different cell types. Injection of platelet-rich plasma is circumstantially considered an aid to muscle repair but whether platelets promote regeneration beyond their role in hemostasis remains unexplored. Here, we find that signaling via platelet-released chemokines is an early event necessary for muscle repair in mice. Platelet depletion reduces the levels of the platelet-secreted neutrophil chemoattractants CXCL5 and CXCL7/PPBP. Consequently, early-phase neutrophil infiltration to injured muscles is impaired whereas later inflammation is exacerbated. Consistent with this model, neutrophil infiltration to injured muscles is compromised in male mice with Cxcl7-knockout platelets. Moreover, neo-angiogenesis and the re-establishment of myofiber size and muscle strength occurs optimally in control mice post-injury but not in Cxcl7ko mice and in neutrophil-depleted mice. Altogether, these findings indicate that platelet-secreted CXCL7 promotes regeneration by recruiting neutrophils to injured muscles, and that this signaling axis could be utilized therapeutically to boost muscle regeneration.

Progressive development of melanoma-induced cachexia differentially impacts organ systems in mice.

Cachexia is a systemic wasting syndrome that increases cancer-associated mortality. How cachexia progressively and differentially impacts distinct tissues is largely unknown. Here, we find that the heart and skeletal muscle undergo wasting at early stages and are the tissues transcriptionally most impacted by cachexia. We also identify general and organ-specific transcriptional changes that indicate functional derangement by cachexia even in tissues that do not undergo wasting, such as the brain. Secreted factors constitute a top category of cancer-regulated genes in host tissues, and these changes include upregulation of the angiotensin-converting enzyme (ACE). ACE inhibition with the drug lisinopril improves muscle force and partially impedes cachexia-induced transcriptional changes, although wasting is not prevented, suggesting that cancer-induced host-secreted factors can regulate tissue function during cachexia. Altogether, by defining prevalent and temporal and tissue-specific responses to cachexia, this resource highlights biomarkers and possible targets for general and tissue-tailored anti-cachexia therapies.

Modulation of protease expression by the transcription factor Ptx1/PITX regulates protein quality control during aging.

Protein quality control is important for healthy aging and is dysregulated in age-related diseases. The autophagy-lysosome and ubiquitin-proteasome are key for proteostasis, but it remains largely unknown whether other proteolytic systems also contribute to maintain proteostasis during aging. Here, we find that expression of proteolytic enzymes (proteases/peptidases) distinct from the autophagy-lysosome and ubiquitin-proteasome systems declines during skeletal muscle aging in Drosophila. Age-dependent protease downregulation undermines proteostasis, as demonstrated by the increase in detergent-insoluble poly-ubiquitinated proteins and pathogenic huntingtin-polyQ levels in response to protease knockdown. Computational analyses identify the transcription factor Ptx1 (homologous to human PITX1/2/3) as a regulator of protease expression. Consistent with this model, Ptx1 protein levels increase with aging, and Ptx1 RNAi counteracts the age-associated downregulation of protease expression. Moreover, Ptx1 RNAi improves muscle protein quality control in a protease-dependent manner and extends lifespan. These findings indicate that proteases and their transcriptional modulator Ptx1 ensure proteostasis during aging.

Innate frequency-discrimination hyperacuity in Williams-Beuren syndrome mice.

Williams-Beuren syndrome (WBS) is a rare disorder caused by hemizygous microdeletion of ∼27 contiguous genes. Despite neurodevelopmental and cognitive deficits, individuals with WBS have spared or enhanced musical and auditory abilities, potentially offering an insight into the genetic basis of auditory perception. Here, we report that the mouse models of WBS have innately enhanced frequency-discrimination acuity and improved frequency coding in the auditory cortex (ACx). Chemogenetic rescue showed frequency-discrimination hyperacuity is caused by hyperexcitable interneurons in the ACx. Haploinsufficiency of one WBS gene, Gtf2ird1, replicated WBS phenotypes by downregulating the neuropeptide receptor VIPR1. VIPR1 is reduced in the ACx of individuals with WBS and in the cerebral organoids derived from human induced pluripotent stem cells with the WBS microdeletion. Vipr1 deletion or overexpression in ACx interneurons mimicked or reversed, respectively, the cellular and behavioral phenotypes of WBS mice. Thus, the Gtf2ird1-Vipr1 mechanism in ACx interneurons may underlie the superior auditory acuity in WBS.

Matricellular Protein WISP2 Is an Endogenous Inhibitor of Collagen Linearization and Cancer Metastasis.

Collagen remodeling contributes to many physiologic and pathologic processes. In primary tumors, the linearization of collagen fibers promotes cancer cell invasion and metastasis and is indicative of poor prognosis. However, it remains unknown whether there are endogenous inhibitors of collagen linearization that could be exploited therapeutically. Here, we show that collagen linearization is controlled by two secreted matricellular proteins with antagonistic functions. Specifically, WISP1 was secreted by cancer cells, bound to type I collagen (Col I), and linearized Col I via its cysteine-rich C-terminal (CT) domain. In contrast, WISP2, which lacks a CT domain, inhibited Col I linearization by preventing WISP1-Col I binding. Analysis of patient data revealed that WISP2 expression is lower in most solid tumors, in comparison with normal tissues. Consequently, genetic or pharmacologic restoration of higher WISP2 levels impaired collagen linearization and prevented tumor cell invasion and metastasis in vivo in models of human and murine breast cancer. Thus, this study uncovers WISP2 as the first inhibitor of collagen linearization ever identified and reveals that collagen architecture can be normalized and metastasis inhibited by therapeutically restoring a high WISP2:WISP1 ratio. SIGNIFICANCE: Two secreted factors, WISP1 and WISP2, antagonistically regulate collagen linearization, and therapeutically increasing the WISP2:WISP1 ratio in tumors limits collagen linearization and inhibits metastasis.See related commentary by Barcus and Longmore, p. 5611.

Oxygen Tension and the VHL-Hif1α Pathway Determine Onset of Neuronal Polarization and Cerebellar Germinal Zone Exit.

Postnatal brain circuit assembly is driven by temporally regulated intrinsic and cell-extrinsic cues that organize neurogenesis, migration, and axo-dendritic specification in post-mitotic neurons. While cell polarity is an intrinsic organizer of morphogenic events, environmental cues in the germinal zone (GZ) instructing neuron polarization and their coupling during postnatal development are unclear. We report that oxygen tension, which rises at birth, and the von Hippel-Lindau (VHL)-hypoxia-inducible factor 1α (Hif1α) pathway regulate polarization and maturation of post-mitotic cerebellar granule neurons (CGNs). At early postnatal stages with low GZ vascularization, Hif1α restrains CGN-progenitor cell-cycle exit. Unexpectedly, cell-intrinsic VHL-Hif1α pathway activation also delays the timing of CGN differentiation, germinal zone exit, and migration initiation through transcriptional repression of the partitioning-defective (Pard) complex. As vascularization proceeds, these inhibitory mechanisms are downregulated, implicating increasing oxygen tension as a critical switch for neuronal polarization and cerebellar GZ exit.

Correlative three-dimensional super-resolution and block-face electron microscopy of whole vitreously frozen cells.

Within cells, the spatial compartmentalization of thousands of distinct proteins serves a multitude of diverse biochemical needs. Correlative super-resolution (SR) fluorescence and electron microscopy (EM) can elucidate protein spatial relationships to global ultrastructure, but has suffered from tradeoffs of structure preservation, fluorescence retention, resolution, and field of view. We developed a platform for three-dimensional cryogenic SR and focused ion beam-milled block-face EM across entire vitreously frozen cells. The approach preserves ultrastructure while enabling independent SR and EM workflow optimization. We discovered unexpected protein-ultrastructure relationships in mammalian cells including intranuclear vesicles containing endoplasmic reticulum-associated proteins, web-like adhesions between cultured neurons, and chromatin domains subclassified on the basis of transcriptional activity. Our findings illustrate the value of a comprehensive multimodal view of ultrastructural variability across whole cells.

Deep multiomics profiling of brain tumors identifies signaling networks downstream of cancer driver genes.

High throughput omics approaches provide an unprecedented opportunity for dissecting molecular mechanisms in cancer biology. Here we present deep profiling of whole proteome, phosphoproteome and transcriptome in two high-grade glioma (HGG) mouse models driven by mutated RTK oncogenes, PDGFRA and NTRK1, analyzing 13,860 proteins and 30,431 phosphosites by mass spectrometry. Systems biology approaches identify numerous master regulators, including 41 kinases and 23 transcription factors. Pathway activity computation and mouse survival indicate the NTRK1 mutation induces a higher activation of AKT downstream targets including MYC and JUN, drives a positive feedback loop to up-regulate multiple other RTKs, and confers higher oncogenic potency than the PDGFRA mutation. A mini-gRNA library CRISPR-Cas9 validation screening shows 56% of tested master regulators are important for the viability of NTRK-driven HGG cells, including TFs (Myc and Jun) and metabolic kinases (AMPKa1 and AMPKa2), confirming the validity of the multiomics integrative approaches, and providing novel tumor vulnerabilities.

Interleukin-15 Enhances Anti-GD2 Antibody-Mediated Cytotoxicity in an Orthotopic PDX Model of Neuroblastoma.

PURPOSE: Immunotherapy with IL2, GM-CSF, and an anti-disialoganglioside (GD2) antibody significantly increases event-free survival in children with high-risk neuroblastoma. However, therapy failure in one third of these patients and IL2-related toxicities pose a major challenge. We compared the immunoadjuvant effects of IL15 with those of IL2 for enhancing antibody-dependent cell-mediated cytotoxicity (ADCC) in neuroblastoma. EXPERIMENTAL DESIGN: We tested ADCC against neuroblastoma patient-derived xenografts (PDX) in vitro and in vivo and examined the functional and migratory properties of NK cells activated with IL2 and IL15. RESULTS: In cell culture, IL15-activated NK cells induced higher ADCC against two GD+ neuroblastoma PDXs than did IL2-activated NK cells (P < 0.001). This effect was dose-dependent (P < 0.001) and was maintained across several effector-to-tumor ratios. As compared with IL2, IL15 also improved chemotaxis of NK cells, leading to higher numbers of tumorsphere-infiltrating NK cells in vitro (P = 0.002). In an orthotopic PDX model, animals receiving chemoimmunotherapy with an anti-GD2 antibody, GM-CSF, and a soluble IL15/IL15Rα complex had greater tumor regression than did those receiving chemotherapy alone (P = 0.012) or combined with anti-GD2 antibody and GM-CSF with (P = 0.016) or without IL2 (P = 0.035). This was most likely due to lower numbers of immature tumor-infiltrating NK cells (DX5+CD27+) after IL15/IL15Rα administration (P = 0.029) and transcriptional upregulation of Gzmd. CONCLUSIONS: The substitution of IL15 for IL2 leads to significant tumor regression in vitro and in vivo and supports clinical testing of IL15 for immunotherapy in pediatric neuroblastoma.

The tumor cell-secreted matricellular protein WISP1 drives pro-metastatic collagen linearization.

Collagen linearization is a hallmark of aggressive tumors and a key pathogenic event that promotes cancer cell invasion and metastasis. Cell-generated mechanical tension has been proposed to contribute to collagen linearization in tumors, but it is unknown whether other mechanisms play prominent roles in this process. Here, we show that the secretome of cancer cells is by itself able to induce collagen linearization independently of cell-generated mechanical forces. Among the tumor cell-secreted factors, we find a key role in this process for the matricellular protein WISP1 (CCN4). Specifically, WISP1 directly binds to type I collagen to promote its linearization in vitro (in the absence of cells) and in vivo in tumors. Consequently, WISP1-induced type I collagen linearization facilitates tumor cell invasion and promotes spontaneous breast cancer metastasis, without significantly affecting gene expression. Furthermore, higher WISP1 expression in tumors from cancer patients correlates with faster progression to metastatic disease and poor prognosis. Altogether, these findings reveal a conceptually novel mechanism whereby pro-metastatic collagen linearization critically depends on a cancer cell-secreted factor.

Circadian gene variants and the skeletal muscle circadian clock contribute to the evolutionary divergence in longevity across Drosophila populations.

Organisms use endogenous clocks to adapt to the rhythmicity of the environment and to synchronize social activities. Although the circadian cycle is implicated in aging, it is unknown whether natural variation in its function contributes to differences in lifespan between populations and whether the circadian clock of specific tissues is key for longevity. We have sequenced the genomes of Drosophila melanogaster strains with exceptional longevity that were obtained via multiple rounds of selection from a parental strain. Comparison of genomic, transcriptomic, and proteomic data revealed that changes in gene expression due to intergenic polymorphisms are associated with longevity and preservation of skeletal muscle function with aging in these strains. Analysis of transcription factors differentially modulated in long-lived versus parental strains indicates a possible role of circadian clock core components. Specifically, there is higher period and timeless and lower cycle expression in the muscle of strains with delayed aging compared to the parental strain. These changes in the levels of circadian clock transcription factors lead to changes in the muscle circadian transcriptome, which includes genes involved in metabolism, proteolysis, and xenobiotic detoxification. Moreover, a skeletal muscle-specific increase in timeless expression extends lifespan and recapitulates some of the transcriptional and circadian changes that differentiate the long-lived from the parental strains. Altogether, these findings indicate that the muscle circadian clock is important for longevity and that circadian gene variants contribute to the evolutionary divergence in longevity across populations.

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

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.

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 (425 ± 123 mL/min/100 g) was within the range of previously reported values. There was no statistical difference between the estimated absolute kidney blood perfusions from the 2 wk of imaging (paired t-test, p = 0.09). We estimated the absolute blood perfusion in the neuroblastoma tumor to be 16.49 and 16.9 mL/min/100 g for the two sets of microbubble administration parameters (Wilcoxon rank-sum test, p = 0.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.

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

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 terminal patterns of cell division. The epigenome of retinoblastomas was more similar to that of the normal retina than that of retina-derived iPSCs, and we identified retina-specific epigenetic memory.

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.