GammaGateR: semi-automated marker gating for single-cell multiplexed imaging.

MOTIVATION: Multiplexed immunofluorescence (mIF) is an emerging assay for multichannel protein imaging that can decipher cell-level spatial features in tissues. However, existing automated cell phenotyping methods, such as clustering, face challenges in achieving consistency across experiments and often require subjective evaluation. As a result, mIF analyses often revert to marker gating based on manual thresholding of raw imaging data. RESULTS: To address the need for an evaluable semi-automated algorithm, we developed GammaGateR, an R package for interactive marker gating designed specifically for segmented cell-level data from mIF images. Based on a novel closed-form gamma mixture model, GammaGateR provides estimates of marker-positive cell proportions and soft clustering of marker-positive cells. The model incorporates user-specified constraints that provide a consistent but slide-specific model fit. We compared GammaGateR against the newest unsupervised approach for annotating mIF data, employing two colon datasets and one ovarian cancer dataset for the evaluation. We showed that GammaGateR produces highly similar results to a silver standard established through manual annotation. Furthermore, we demonstrated its effectiveness in identifying biological signals, achieved by mapping known spatial interactions between CD68 and MUC5AC cells in the colon and by accurately predicting survival in ovarian cancer patients using the phenotype probabilities as input for machine learning methods. GammaGateR is a highly efficient tool that can improve the replicability of marker gating results, while reducing the time of manual segmentation. AVAILABILITY: The R package is available at https://github.com/JiangmeiRubyXiong/GammaGateR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Detection of Treatment Response in Triple-Negative Breast Tumors to Paclitaxel Using MRI Cell Size Imaging.

BACKGROUND: Breast cancer treatment response evaluation using the response evaluation criteria in solid tumors (RECIST) guidelines, based on tumor volume changes, has limitations, prompting interest in novel imaging markers for accurate therapeutic effect determination. PURPOSE: To use MRI-measured cell size as a new imaging biomarker for assessing chemotherapy response in breast cancer. STUDY TYPE: Longitudinal; animal model. STUDY POPULATION: Triple-negative human breast cancer cell (MDA-MB-231) pellets (4 groups, n = 7) treated with dimethyl sulfoxide (DMSO) or 10 nM of paclitaxel for 24, 48, and 96 hours, and 29 mice with MDA-MB-231 tumors in right hind limbs treated with paclitaxel (n = 16) or DMSO (n = 13) twice weekly for 3 weeks. FIELD STRENGTH/SEQUENCE: Oscillating gradient spin echo and pulsed gradient spin echo sequences at 4.7 T. ASSESSMENT: MDA-MB-231 cells were analyzed using flowcytometry and light microscopy to assess cell cycle phases and cell size distribution. MDA-MB-231 cell pellets were MR imaged. Mice were imaged weekly, with 9, 6, and 14 being sacrificed for histology after MRI at weeks 1, 2, and 3, respectively. Microstructural parameters of tumors/cell pellets were derived by fitting diffusion MRI data to a biophysical model. STATISTICAL TESTS: One-way ANOVA compared cell sizes and MR-derived parameters between treated and control samples. Repeated measures 2-way ANOVA with Bonferroni post-tests compared temporal changes in MR-derived parameters. A P-value <0.05 was considered statistically significant. RESULTS: In vitro experiments showed that the mean MR-derived cell sizes of paclitaxel-treated cells increased significantly with a 24-hours treatment and decreased (P = 0.06) with a 96-hour treatment. For in vivo xenograft experiments, the paclitaxel-treated tumors showed significant decreases in cell size at later weeks. MRI observations were supported by flowcytometry, light microscopy, and histology. DATA CONCLUSIONS: MR-derived cell size may characterize the cell shrinkage during treatment-induced apoptosis, and may potentially provide new insights into the assessment of therapeutic response. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 4.

Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors.

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

GammaGateR: semi-automated marker gating for single-cell multiplexed imaging.

Motivation: Multiplexed immunofluorescence (mIF) is an emerging assay for multichannel protein imaging that can decipher cell-level spatial features in tissues. However, existing automated cell phenotyping methods, such as clustering, face challenges in achieving consistency across experiments and often require subjective evaluation. As a result, mIF analyses often revert to marker gating based on manual thresholding of raw imaging data. Results: To address the need for an evaluable semi-automated algorithm, we developed GammaGateR, an R package for interactive marker gating designed specifically for segmented cell-level data from mIF images. Based on a novel closed-form gamma mixture model, GammaGateR provides estimates of marker-positive cell proportions and soft clustering of marker-positive cells. The model incorporates user-specified constraints that provide a consistent but slide-specific model fit. We compared GammaGateR against the newest unsupervised approach for annotating mIF data, employing two colon datasets and one ovarian cancer dataset for the evaluation. We showed that GammaGateR produces highly similar results to a silver standard established through manual annotation. Furthermore, we demonstrated its effectiveness in identifying biological signals, achieved by mapping known spatial interactions between CD68 and MUC5AC cells in the colon and by accurately predicting survival in ovarian cancer patients using the phenotype probabilities as input for machine learning methods. GammaGateR is a highly efficient tool that can improve the replicability of marker gating results, while reducing the time of manual segmentation. Availability and Implementation: The R package is available at https://github.com/JiangmeiRubyXiong/GammaGateR.

Multiplex immunohistochemistry and high-throughput image analysis for evaluation of spatial tumor immune cell markers in human breast cancer.

BACKGROUND: The clinicopathological significance of spatial tumor-infiltrating lymphocytes (TILs) subpopulations is not well studied due to lack of high-throughput scalable methodology for studies with large human sample sizes. OBJECTIVE: Establishing a cyclic fluorescent multiplex immunohistochemistry (mIHC/IF) method coupled with computer-assisted high-throughput quantitative analysis to evaluate associations of six TIL markers (CD3, CD8, CD20, CD56, FOXP3, and PD-L1) with clinicopathological factors of breast cancer. METHODS: Our 5-plex mIHC/IF staining was shown to be reliable and highly sensitive for labeling three biomarkers per tissue section. Through repetitive cycles of 5-plex mIHC/IF staining, more than 12 biomarkers could be detected per single tissue section. Using open-source software CellProfiler, the measurement pipelines were successfully developed for high-throughput multiplex evaluation of intratumoral and stromal TILs. RESULTS: In analyses of 188 breast cancer samples from the Nashville Breast Health Study, high-grade tumors showed significantly increased intratumoral CD3+CD8+ cytotoxic T lymphocyte density (P= 0.0008, false discovery rate (FDR) adjusted P= 0.0168) and intratumoral PD-L1 expression (P= 0.0061, FDR adjusted P= 0.0602) compared with low-grade tumors. CONCLUSIONS: The high- and low-grade breast cancers exhibit differential immune responses which may have clinical significance. The multiplexed imaging quantification strategies established in this study are reliable, cost-efficient and applicable in regular laboratory settings for high-throughput tissue biomarker studies, especially retrospective and population-based studies using archived paraffin tissues.

Selective Cell Size MRI Differentiates Brain Tumors from Radiation Necrosis.

Brain metastasis is a common characteristic of late-stage lung cancers. High doses of targeted radiotherapy can control tumor growth in the brain but can also result in radiotherapy-induced necrosis. Current methods are limited for distinguishing whether new parenchymal lesions following radiotherapy are recurrent tumors or radiotherapy-induced necrosis, but the clinical management of these two classes of lesions differs significantly. Here, we developed, validated, and evaluated a new MRI technique termed selective size imaging using filters via diffusion times (SSIFT) to differentiate brain tumors from radiotherapy necrosis in the brain. This approach generates a signal filter that leverages diffusion time dependence to establish a cell size-weighted map. Computer simulations in silico, cultured cancer cells in vitro, and animals with brain tumors in vivo were used to comprehensively validate the specificity of SSIFT for detecting typical large cancer cells and the ability to differentiate brain tumors from radiotherapy necrosis. SSIFT was also implemented in patients with metastatic brain cancer and radiotherapy necrosis. SSIFT showed high correlation with mean cell sizes in the relevant range of less than 20 µm. The specificity of SSIFT for brain tumors and reduced contrast in other brain etiologies allowed SSIFT to differentiate brain tumors from peritumoral edema and radiotherapy necrosis. In conclusion, this new, cell size-based MRI method provides a unique contrast to differentiate brain tumors from other pathologies in the brain. SIGNIFICANCE: This work introduces and provides preclinical validation of a new diffusion MRI method that exploits intrinsic differences in cell sizes to distinguish brain tumors and radiotherapy necrosis.

Cancer-Associated Fibroblasts and Squamous Epithelial Cells Constitute a Unique Microenvironment in a Mouse Model of Inflammation-Induced Colon Cancer.

The tumor microenvironment plays a key role in the pathogenesis of colorectal tumors and contains various cell types including epithelial, immune, and mesenchymal cells. Characterization of the interactions between these cell types is necessary for revealing the complex nature of tumors. In this study, we used single-cell RNA-seq (scRNA-seq) to compare the tumor microenvironments between a mouse model of sporadic colorectal adenoma (Lrig1CreERT2/+;Apc2lox14/+) and a mouse model of inflammation-driven colorectal cancer induced by azoxymethane and dextran sodium sulfate (AOM/DSS). While both models develop tumors in the distal colon, we found that the two tumor types have distinct microenvironments. AOM/DSS tumors have an increased abundance of two populations of cancer-associated fibroblasts (CAFs) compared with APC tumors, and we revealed their divergent spatial association with tumor cells using multiplex immunofluorescence (MxIF) imaging. We also identified a unique squamous cell population in AOM/DSS tumors, whose origins were distinct from anal squamous epithelial cells. These cells were in higher proportions upon administration of a chemotherapy regimen of 5-Fluorouracil/Irinotecan. We used computational inference algorithms to predict cell-cell communication mediated by ligand-receptor interactions and downstream pathway activation, and identified potential mechanistic connections between CAFs and tumor cells, as well as CAFs and squamous epithelial cells. This study provides important preclinical insight into the microenvironment of two distinct models of colorectal tumors and reveals unique roles for CAFs and squamous epithelial cells in the AOM/DSS model of inflammation-driven cancer.

Quantifying and correcting slide-to-slide variation in multiplexed immunofluorescence images.

MOTIVATION: Multiplexed imaging is a nascent single-cell assay with a complex data structure susceptible to technical variability that disrupts inference. These in situ methods are valuable in understanding cell-cell interactions, but few standardized processing steps or normalization techniques of multiplexed imaging data are available. RESULTS: We implement and compare data transformations and normalization algorithms in multiplexed imaging data. Our methods adapt the ComBat and functional data registration methods to remove slide effects in this domain, and we present an evaluation framework to compare the proposed approaches. We present clear slide-to-slide variation in the raw, unadjusted data and show that many of the proposed normalization methods reduce this variation while preserving and improving the biological signal. Furthermore, we find that dividing multiplexed imaging data by its slide mean, and the functional data registration methods, perform the best under our proposed evaluation framework. In summary, this approach provides a foundation for better data quality and evaluation criteria in multiplexed imaging. AVAILABILITY AND IMPLEMENTATION: Source code is provided at: https://github.com/statimagcoll/MultiplexedNormalization and an R package to implement these methods is available here: https://github.com/ColemanRHarris/mxnorm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MIRIAM: A machine and deep learning single-cell segmentation and quantification pipeline for multi-dimensional tissue images.

Increasingly, highly multiplexed tissue imaging methods are used to profile protein expression at the single-cell level. However, a critical limitation is the lack of robust cell segmentation tools for tissue sections. We present Multiplexed Image Resegmentation of Internal Aberrant Membranes (MIRIAM) that combines (a) a pipeline for cell segmentation and quantification that incorporates machine learning-based pixel classification to define cellular compartments, (b) a novel method for extending incomplete cell membranes, and (c) a deep learning-based cell shape descriptor. Using human colonic adenomas as an example, we show that MIRIAM is superior to widely utilized segmentation methods and provides a pipeline that is broadly applicable to different imaging platforms and tissue types.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Supermeres are functional extracellular nanoparticles replete with disease biomarkers and therapeutic targets.

Extracellular vesicles and exomere nanoparticles are under intense investigation as sources of clinically relevant cargo. Here we report the discovery of a distinct extracellular nanoparticle, termed supermere. Supermeres are morphologically distinct from exomeres and display a markedly greater uptake in vivo compared with small extracellular vesicles and exomeres. The protein and RNA composition of supermeres differs from small extracellular vesicles and exomeres. Supermeres are highly enriched with cargo involved in multiple cancers (glycolytic enzymes, TGFBI, miR-1246, MET, GPC1 and AGO2), Alzheimer's disease (APP) and cardiovascular disease (ACE2, ACE and PCSK9). The majority of extracellular RNA is associated with supermeres rather than small extracellular vesicles and exomeres. Cancer-derived supermeres increase lactate secretion, transfer cetuximab resistance and decrease hepatic lipids and glycogen in vivo. This study identifies a distinct functional nanoparticle replete with potential circulating biomarkers and therapeutic targets for a host of human diseases.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Induction of apically mistrafficked epiregulin disrupts epithelial polarity via aberrant EGFR signaling.

In polarized MDCK cells, disruption of the tyrosine-based YXXΦ basolateral trafficking motif (Y156A) in the epidermal growth factor receptor (EGFR) ligand epiregulin (EREG), results in its apical mistrafficking and transformation in vivo. However, the mechanisms underlying these dramatic effects are unknown. Using a doxycycline-inducible system in 3D Matrigel cultures, we now show that induction of Y156A EREG in fully formed MDCK cysts results in direct and complete delivery of mutant EREG to the apical cell surface. Within 3 days of induction, ectopic lumens were detected in mutant, but not wild-type, EREG-expressing cysts. Of note, these structures resembled histological features found in subcutaneous xenografts of mutant EREG-expressing MDCK cells. These ectopic lumens formed de novo rather than budding from the central lumen and depended on metalloprotease-mediated cleavage of EREG and subsequent EGFR activity. Moreover, the most frequent EREG mutation in human cancer (R147stop) resulted in its apical mistrafficking in engineered MDCK cells. Thus, induction of EREG apical mistrafficking is sufficient to disrupt selective aspects of polarity of a preformed polarized epithelium. This article has an associated First Person interview with the first author of the paper.

MRI-cytometry: Mapping nonparametric cell size distributions using diffusion MRI.

PURPOSE: This report introduces and validates a new diffusion MRI-based method, termed MRI-cytometry, which can noninvasively map intravoxel, nonparametric cell size distributions in tissues. METHODS: MRI was used to acquire diffusion MRI signals with a range of diffusion times and gradient factors, and a model was fit to these data to derive estimates of cell size distributions. We implemented a 2-step fitting method to avoid noise-induced artificial peaks and provide reliable estimates of tumor cell size distributions. Computer simulations in silico, experimental measurements on cultured cells in vitro, and animal xenografts in vivo were used to validate the accuracy and precision of the method. Tumors in 7 patients with breast cancer were also imaged and analyzed using this MRI-cytometry approach on a clinical 3 Tesla MRI scanner. RESULTS: Simulations and experimental results confirm that MRI-cytometry can reliably map intravoxel, nonparametric cell size distributions and has the potential to discriminate smaller and larger cells. The application in breast cancer patients demonstrates the feasibility of direct translation of MRI-cytometry to clinical applications. CONCLUSION: The proposed MRI-cytometry method can characterize nonparametric cell size distributions in human tumors, which potentially provides a practical imaging approach to derive specific histopathological information on biological tissues.

Succinate Produced by Intestinal Microbes Promotes Specification of Tuft Cells to Suppress Ileal Inflammation.

BACKGROUND & AIMS: Countries endemic for parasitic infestations have a lower incidence of Crohn's disease (CD) than nonendemic countries, and there have been anecdotal reports of the beneficial effects of helminths in CD patients. Tuft cells in the small intestine sense and direct the immune response against eukaryotic parasites. We investigated the activities of tuft cells in patients with CD and mouse models of intestinal inflammation. METHODS: We used microscopy to quantify tuft cells in intestinal specimens from patients with ileal CD (n = 19), healthy individuals (n = 14), and TNFΔARE/+ mice, which develop Crohn's-like ileitis. We performed single-cell RNA sequencing, mass spectrometry, and microbiome profiling of intestinal tissues from wild-type and Atoh1-knockout mice, which have expansion of tuft cells, to study interactions between microbes and tuft cell populations. We assessed microbe dependence of tuft cell populations using microbiome depletion, organoids, and microbe transplant experiments. We used multiplex imaging and cytokine assays to assess alterations in inflammatory response following expansion of tuft cells with succinate administration in TNFΔARE/+ and anti-CD3E CD mouse models. RESULTS: Inflamed ileal tissues from patients and mice had reduced numbers of tuft cells, compared with healthy individuals or wild-type mice. Expansion of tuft cells was associated with increased expression of genes that regulate the tricarboxylic acid cycle, which resulted from microbe production of the metabolite succinate. Experiments in which we manipulated the intestinal microbiota of mice revealed the existence of an ATOH1-independent population of tuft cells that was sensitive to metabolites produced by microbes. Administration of succinate to mice expanded tuft cells and reduced intestinal inflammation in TNFΔARE/+ mice and anti-CD3E-treated mice, increased GATA3+ cells and type 2 cytokines (IL22, IL25, IL13), and decreased RORGT+ cells and type 17 cytokines (IL23) in a tuft cell-dependent manner. CONCLUSIONS: We found that tuft cell expansion reduced chronic intestinal inflammation in mice. Strategies to expand tuft cells might be developed for treatment of CD.

Magnetic resonance imaging of mean cell size in human breast tumors.

PURPOSE: Cell size is a fundamental characteristic of all tissues, and changes in cell size in cancer reflect tumor status and response to treatments, such as apoptosis and cell-cycle arrest. Unfortunately, cell size can currently be obtained only by pathological evaluation of tumor tissue samples obtained invasively. Previous imaging approaches are limited to preclinical MRI scanners or require relatively long acquisition times that are impractical for clinical imaging. There is a need to develop cell-size imaging for clinical applications. METHODS: We propose a clinically feasible IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) approach that can characterize mean cell sizes in solid tumors. We report the use of a combination of pulse sequences, using different gradient waveforms implemented on clinical MRI scanners and analytical equations based on these waveforms to analyze diffusion-weighted MRI signals and derive specific microstructural parameters such as cell size. We also describe comprehensive validations of this approach using computer simulations, cell experiments in vitro, and animal experiments in vivo and demonstrate applications in preoperative breast cancer patients. RESULTS: With fast acquisitions (~7 minutes), IMPULSED can provide high-resolution (1.3 mm in-plane) mapping of mean cell size of human tumors in vivo on clinical 3T MRI scanners. All validations suggest that IMPULSED provides accurate and reliable measurements of mean cell size. CONCLUSION: The proposed IMPULSED method can assess cell-size variations in tumors of breast cancer patients, which may have the potential to assess early response to neoadjuvant therapy.

In vivo magnetic resonance imaging of treatment-induced apoptosis.

Imaging apoptosis could provide an early and specific means to monitor tumor responses to treatment. To date, despite numerous attempts to develop molecular imaging approaches, there is still no widely-accepted and reliable method for in vivo imaging of apoptosis. We hypothesized that the distinct cellular morphologic changes associated with treatment-induced apoptosis, such as cell shrinkage, cytoplasm condensation, and DNA fragmentation, can be detected by temporal diffusion spectroscopy imaging (TDSI). Cetuximab-induced apoptosis was assessed in vitro and in vivo with cetuximab-sensitive (DiFi) and insensitive (HCT-116) human colorectal cancer cell lines by TDSI. TDSI findings were complemented by flow cytometry and immunohistochemistry. Cell cycle analysis and flow cytometry detected apoptotic cell shrinkage in cetuximab-treated DiFi cells, and significant apoptosis was confirmed by histology. TDSI-derived parameters quantified key morphological changes including cell size decreases during apoptosis in responsive tumors that occurred earlier than gross tumor volume regression. TDSI provides a unique measurement of apoptosis by identifying cellular characteristics, particularly cell shrinkage. The method will assist in understanding the underlying biology of solid tumors and predict tumor response to therapies. TDSI is free of any exogenous agent or radiation, and hence is very suitable to be incorporated into clinical applications.

Transfer of Functional Cargo in Exomeres.

Exomeres are a recently discovered type of extracellular nanoparticle with no known biological function. Herein, we describe a simple ultracentrifugation-based method for separation of exomeres from exosomes. Exomeres are enriched in Argonaute 1-3 and amyloid precursor protein. We identify distinct functions of exomeres mediated by two of their cargo, the ß-galactoside α2,6-sialyltransferase 1 (ST6Gal-I) that α2,6- sialylates N-glycans, and the EGFR ligand, amphiregulin (AREG). Functional ST6Gal-I in exomeres can be transferred to cells, resulting in hypersialylation of recipient cell-surface proteins including ß1-integrin. AREG-containing exomeres elicit prolonged EGFR and downstream signaling in recipient cells, modulate EGFR trafficking in normal intestinal organoids, and dramatically enhance the growth of colonic tumor organoids. This study provides a simplified method of exomere isolation and demonstrates that exomeres contain and can transfer functional cargo. These findings underscore the heterogeneity of nanoparticles and should accelerate advances in determining the composition and biological functions of exomeres.