Spatial mapping of protein composition and tissue organization: a primer for multiplexed antibody-based imaging.

Tissues and organs are composed of distinct cell types that must operate in concert to perform physiological functions. Efforts to create high-dimensional biomarker catalogs of these cells have been largely based on single-cell sequencing approaches, which lack the spatial context required to understand critical cellular communication and correlated structural organization. To probe in situ biology with sufficient depth, several multiplexed protein imaging methods have been recently developed. Though these technologies differ in strategy and mode of immunolabeling and detection tags, they commonly utilize antibodies directed against protein biomarkers to provide detailed spatial and functional maps of complex tissues. As these promising antibody-based multiplexing approaches become more widely adopted, new frameworks and considerations are critical for training future users, generating molecular tools, validating antibody panels, and harmonizing datasets. In this Perspective, we provide essential resources, key considerations for obtaining robust and reproducible imaging data, and specialized knowledge from domain experts and technology developers.

FLINO: a new method for immunofluorescence bioimage normalization.

MOTIVATION: Multiplexed immunofluorescence bioimaging of single-cells and their spatial organization in tissue holds great promise to the development of future precision diagnostics and therapeutics. Current multiplexing pipelines typically involve multiple rounds of immunofluorescence staining across multiple tissue slides. This introduces experimental batch effects that can hide underlying biological signal. It is important to have robust algorithms that can correct for the batch effects while not introducing biases into the data. Performance of data normalization methods can vary among different assay pipelines. To evaluate differences, it is critical to have a ground truth dataset that is representative of the assay. RESULTS: A new immunoFLuorescence Image NOrmalization method is presented and evaluated against alternative methods and workflows. Multiround immunofluorescence staining of the same tissue with the nuclear dye DAPI was used to represent virtual slides and a ground truth. DAPI was restained on a given tissue slide producing multiple images of the same underlying structure but undergoing multiple representative tissue handling steps. This ground truth dataset was used to evaluate and compare multiple normalization methods including median, quantile, smooth quantile, median ratio normalization and trimmed mean of the M-values. These methods were applied in both an unbiased grid object and segmented cell object workflow to 24 multiplexed biomarkers. An upper quartile normalization of grid objects in log space was found to obtain almost equivalent performance to directly normalizing segmented cell objects by the middle quantile. The developed grid-based technique was then applied with on-slide controls for evaluation. Using five or fewer controls per slide can introduce biases into the data. Ten or more on-slide controls were able to robustly correct for batch effects. AVAILABILITY AND IMPLEMENTATION: The data underlying this article along with the FLINO R-scripts used to perform the evaluation of image normalizations methods and workflows can be downloaded from https://github.com/GE-Bio/FLINO. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Stratification of chemotherapy-treated stage III colorectal cancer patients using multiplexed imaging and single-cell analysis of T-cell populations.

Colorectal cancer (CRC) has one of the highest cancer incidences and mortality rates. In stage III, postoperative chemotherapy benefits <20% of patients, while more than 50% will develop distant metastases. Biomarkers for identification of patients at increased risk of disease recurrence following adjuvant chemotherapy are currently lacking. In this study, we assessed immune signatures in the tumor and tumor microenvironment (TME) using an in situ multiplexed immunofluorescence imaging and single-cell analysis technology (Cell DIVETM) and evaluated their correlations with patient outcomes. Tissue microarrays (TMAs) with up to three 1 mm diameter cores per patient were prepared from 117 stage III CRC patients treated with adjuvant fluoropyrimidine/oxaliplatin (FOLFOX) chemotherapy. Single sections underwent multiplexed immunofluorescence staining for immune cell markers (CD45, CD3, CD4, CD8, FOXP3, PD1) and tumor/cell segmentation markers (DAPI, pan-cytokeratin, AE1, NaKATPase, and S6). We used annotations and a probabilistic classification algorithm to build statistical models of immune cell types. Images were also qualitatively assessed independently by a Pathologist as 'high', 'moderate' or 'low', for stromal and total immune cell content. Excellent agreement was found between manual assessment and total automated scores (p < 0.0001). Moreover, compared to single markers, a multi-marker classification of regulatory T cells (Tregs: CD3+/CD4+FOXP3+/PD1-) was significantly associated with disease-free survival (DFS) and overall survival (OS) (p = 0.049 and 0.032) of FOLFOX-treated patients. Our results also showed that PD1- Tregs rather than PD1+ Tregs were associated with improved survival. These findings were supported by results from an independent FOLFOX-treated cohort of 191 stage III CRC patients, where higher PD1- Tregs were associated with an increase overall survival (p = 0.015) for CD3+/CD4+/FOXP3+/PD1-. Overall, compared to single markers, multi-marker classification provided more accurate quantitation of immune cell types with stronger correlations with outcomes.

Multi-protein spatial signatures in ductal carcinoma in situ (DCIS) of breast.

BACKGROUND: There is limited knowledge about DCIS cellular composition and relationship with breast cancer events (BCE). METHODS: Immunofluorescence multiplexing (MxIF) was used to image and quantify 32 cellular biomarkers in FFPE DCIS tissue microarrays. Over 75,000 DCIS cells from 51 patients (median 9 years follow-up for non-BCE cases) were analysed for profiles predictive of BCE. K-means clustering was used to evaluate cellular co-expression of epithelial markers with ER and HER2. RESULTS: Only ER, PR and HER2 significantly correlated with BCE. Cluster analysis identified 6 distinct cell groups with different levels of ER, Her2, cMET and SLC7A5. Clusters 1 and 3 were not significant. Clusters 2 and 4 (high ER/low HER2 and SLC7A5/mixed cMET) significantly correlated with low BCE risk (P = 0.001 and P = 0.034), while cluster 6 (high HER2/low ER, cMET and SLC7A5) correlated with increased risk (P = 0.018). Cluster 5 (similar to cluster 6, except high SLC7A5) trended towards significance (P = 0.072). A continuous expression score (Escore) based on these 4 clusters predicted likelihood of BCE (AUC = 0.79, log-rank test P = 5E-05; LOOCV AUC = 0.74, log-rank test P = 0.006). CONCLUSION: Multiplexed spatial analysis of limited tissue is a novel method for biomarker analysis and predicting BCEs. Further validation of Escore is needed in a larger cohort.

Highly multiplexed single-cell analysis of formalin-fixed, paraffin-embedded cancer tissue.

Limitations on the number of unique protein and DNA molecules that can be characterized microscopically in a single tissue specimen impede advances in understanding the biological basis of health and disease. Here we present a multiplexed fluorescence microscopy method (MxIF) for quantitative, single-cell, and subcellular characterization of multiple analytes in formalin-fixed paraffin-embedded tissue. Chemical inactivation of fluorescent dyes after each image acquisition round allows reuse of common dyes in iterative staining and imaging cycles. The mild inactivation chemistry is compatible with total and phosphoprotein detection, as well as DNA FISH. Accurate computational registration of sequential images is achieved by aligning nuclear counterstain-derived fiducial points. Individual cells, plasma membrane, cytoplasm, nucleus, tumor, and stromal regions are segmented to achieve cellular and subcellular quantification of multiplexed targets. In a comparison of pathologist scoring of diaminobenzidine staining of serial sections and automated MxIF scoring of a single section, human epidermal growth factor receptor 2, estrogen receptor, p53, and androgen receptor staining by diaminobenzidine and MxIF methods yielded similar results. Single-cell staining patterns of 61 protein antigens by MxIF in 747 colorectal cancer subjects reveals extensive tumor heterogeneity, and cluster analysis of divergent signaling through ERK1/2, S6 kinase 1, and 4E binding protein 1 provides insights into the spatial organization of mechanistic target of rapamycin and MAPK signal transduction. Our results suggest MxIF should be broadly applicable to problems in the fields of basic biological research, drug discovery and development, and clinical diagnostics.

3D reconstruction of skin and spatial mapping of immune cell density, vascular distance and effects of sun exposure and aging.

Mapping the human body at single cell resolution in three dimensions (3D) is important for understanding cellular interactions in context of tissue and organ organization. 2D spatial cell analysis in a single tissue section may be limited by cell numbers and histology. Here we show a workflow for 3D reconstruction of multiplexed sequential tissue sections: MATRICS-A (Multiplexed Image Three-D Reconstruction and Integrated Cell Spatial - Analysis). We demonstrate MATRICS-A in 26 serial sections of fixed skin (stained with 18 biomarkers) from 12 donors aged between 32-72 years. Comparing the 3D reconstructed cellular data with the 2D data, we show significantly shorter distances between immune cells and vascular endothelial cells (56 µm in 3D vs 108 µm in 2D). We also show 10-70% more T cells (total) within 30 µm of a neighboring T helper cell in 3D vs 2D. Distances of p53, DDB2 and Ki67 positive cells to the skin surface were consistent across all ages/sun exposure and largely localized to the lower stratum basale layer of the epidermis. MATRICS-A provides a framework for analysis of 3D spatial cell relationships in healthy and aging organs and could be further extended to diseased organs.

An atlas of inter- and intra-tumor heterogeneity of apoptosis competency in colorectal cancer tissue at single-cell resolution.

Cancer cells' ability to inhibit apoptosis is key to malignant transformation and limits response to therapy. Here, we performed multiplexed immunofluorescence analysis on tissue microarrays with 373 cores from 168 patients, segmentation of 2.4 million individual cells, and quantification of 18 cell lineage and apoptosis proteins. We identified an enrichment for BCL2 in immune, and BAK, SMAC, and XIAP in cancer cells. Ordinary differential equation-based modeling of apoptosis sensitivity at single-cell resolution was conducted and an atlas of inter- and intra-tumor heterogeneity in apoptosis susceptibility generated. Systems modeling at single-cell resolution identified an enhanced sensitivity of cancer cells to mitochondrial permeabilization and executioner caspase activation compared to immune and stromal cells, but showed significant inter- and intra-tumor heterogeneity.

Novel DNA-based T-Cell Activator Promotes Rapid T-Cell Activation and Expansion.

Autologous chimeric antigen receptor engineered T-cell therapies are beginning to dramatically change the outlook for patients with several hematological malignancies. Yet methods to activate and expand these cells are limited, often pose challenges to automation, and have biological limitations impacting the output of the injectable dose. This study describes the development of a novel, highly flexible, soluble DNA-based T-cell activation and expansion platform which alleviates the limitations of current technologies and provides rapid T-cell activation and expansion.

Comparison of Multiplexed Immunofluorescence Imaging to Chromogenic Immunohistochemistry of Skin Biomarkers in Response to Monkeypox Virus Infection.

Over the last 15 years, advances in immunofluorescence-imaging based cycling methods, antibody conjugation methods, and automated image processing have facilitated the development of a high-resolution, multiplexed tissue immunofluorescence (MxIF) method with single cell-level quantitation termed Cell DIVETM. Originally developed for fixed oncology samples, here it was evaluated in highly fixed (up to 30 days), archived monkeypox virus-induced inflammatory skin lesions from a retrospective study in 11 rhesus monkeys to determine whether MxIF was comparable to manual H-scoring of chromogenic stains. Six protein markers related to immune and cellular response (CD68, CD3, Hsp70, Hsp90, ERK1/2, ERK1/2 pT202_pY204) were manually quantified (H-scores) by a pathologist from chromogenic IHC double stains on serial sections and compared to MxIF automated single cell quantification of the same markers that were multiplexed on a single tissue section. Overall, there was directional consistency between the H-score and the MxIF results for all markers except phosphorylated ERK1/2 (ERK1/2 pT202_pY204), which showed a decrease in the lesion compared to the adjacent non-lesioned skin by MxIF vs an increase via H-score. Improvements to automated segmentation using machine learning and adding additional cell markers for cell viability are future options for improvement. This method could be useful in infectious disease research as it conserves tissue, provides marker colocalization data on thousands of cells, allowing further cell level data mining as well as a reduction in user bias.

Understanding heterogeneous tumor microenvironment in metastatic melanoma.

A systemic analysis of the tumor-immune interactions within the heterogeneous tumor microenvironment is of particular importance for understanding the antitumor immune response. We used multiplexed immunofluorescence to elucidate cellular spatial interactions and T-cell infiltrations in metastatic melanoma tumor microenvironment. We developed two novel computational approaches that enable infiltration clustering and single cell analysis-cell aggregate algorithm and cell neighborhood analysis algorithm-to reveal and to compare the spatial distribution of various immune cells relative to tumor cell in sub-anatomic tumor microenvironment areas. We showed that the heterogeneous tumor human leukocyte antigen-1 expressions differently affect the magnitude of cytotoxic T-cell infiltration and the distributions of CD20+ B cells and CD4+FOXP3+ regulatory T cells within and outside of T-cell infiltrated tumor areas. In a cohort of 166 stage III melanoma samples, high tumor human leukocyte antigen-1 expression is required but not sufficient for high T-cell infiltration, with significantly improved overall survival. Our results demonstrate that tumor cells with heterogeneous properties are associated with differential but predictable distributions of immune cells within heterogeneous tumor microenvironment with various biological features and impacts on clinical outcomes. It establishes tools necessary for systematic analysis of the tumor microenvironment, allowing the elucidation of the "homogeneous patterns" within the heterogeneous tumor microenvironment.

Efficacy and tolerability of anti-programmed death-ligand 1 (PD-L1) antibody (Avelumab) treatment in advanced thymoma.

BACKGROUND: Thymic epithelial tumors are PD-L1-expressing tumors of thymic epithelial origin characterized by varying degrees of lymphocytic infiltration and a predisposition towards development of paraneoplastic autoimmunity. PD-1-targeting antibodies have been evaluated, largely in patients with thymic carcinoma. We sought to evaluate the efficacy and safety of the anti-PD-L1 antibody, avelumab (MSB0010718C), in patients with relapsed, advanced thymic epithelial tumors and conduct correlative immunological studies. METHODS: Seven patients with thymoma and one patient with thymic carcinoma were enrolled in a phase I, dose-escalation trial of avelumab (MSB0010718C), and treated with avelumab at doses of 10 mg/kg to 20 mg/kg every 2 weeks until disease progression or development of intolerable side effects. Tissue and blood immunological analyses were conducted. RESULTS: Two of seven (29%) patients with thymoma had a confirmed Response Evaluation Criteria in Solid Tumors-defined partial response, two (29%) had an unconfirmed partial response and three patients (two thymoma; one thymic carcinoma) had stable disease (43%). Three of four responses were observed after a single dose of avelumab. All responders developed immune-related adverse events that resolved with immunosuppressive therapy. Only one of four patients without a clinical response developed immune-related adverse events. Responders had a higher absolute lymphocyte count, lower frequencies of B cells, regulatory T cells, conventional dendritic cells, and natural killer cells prior to therapy. CONCLUSION: These results demonstrate anti-tumor activity of PD-L1 inhibition in patients with relapsed thymoma accompanied by a high frequency of immune-related adverse events. Pre-treatment immune cell subset populations differ between responders and non-responders. TRIAL REGISTRATION: ClinicalTrials.gov - NCT01772004 . Date of registration - January 21, 2013.

Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas.

Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations.

A multimodal contrast agent for preoperative MR Imaging and intraoperative tumor margin delineation.

We have constructed a multimodal contrast agent suitable for near-infrared, NIR, fluorescent imaging as well as magnetic resonance imaging, MRI. This class of agents may be useful for preoperative tumor localization and tumor functional evaluation and for intraoperative delineation of tumor margins. We have covalently attached dyes of the cyanine family to a previously described polymeric contrast agent, Gd-DTPA-polylysine, of an extended, uncoiled conformation. The dual modality agent is as effective in imaging tumors by MRI as the parent compound provided that the dye loading on the polymer is such that it does not eliminate all the available free-lysine groups on the parent Gd-DTPA-polylysine polymers. NIR fluorescence from preclinical subcutaneous and orthotopic mammary gland tumors could be detected with a signal to background ratio of as high as 4.5 at 12 hours post agent injection at a dye dose of 125 nmole/kg. For intraoperative delineation of tumor margins, a wide-field illumination camera system was devised giving high signal to background NIR fluorescent images of surgically exposed orthotopic mammary gland tumors. Histologic microscopy confirmed the location of the dual modality agent at the boundary of the tumor with a margin distance of about 0.3 mm from labeled tumor cells.

Multiplexed immunofluorescence delineates proteomic cancer cell states associated with metabolism.

The phenotypic diversity of cancer results from genetic and nongenetic factors. Most studies of cancer heterogeneity have focused on DNA alterations, as technologies for proteomic measurements in clinical specimen are currently less advanced. Here, we used a multiplexed immunofluorescence staining platform to measure the expression of 27 proteins at the single-cell level in formalin-fixed and paraffin-embedded samples from treatment-naive stage II/III human breast cancer. Unsupervised clustering of protein expression data from 638,577 tumor cells in 26 breast cancers identified 8 clusters of protein coexpression. In about one-third of breast cancers, over 95% of all neoplastic cells expressed a single protein coexpression cluster. The remaining tumors harbored tumor cells representing multiple protein coexpression clusters, either in a regional distribution or intermingled throughout the tumor. Tumor uptake of the radiotracer 18F-fluorodeoxyglucose was associated with protein expression clusters characterized by hormone receptor loss, PTEN alteration, and HER2 gene amplification. Our study demonstrates an approach to generate cellular heterogeneity metrics in routinely collected solid tumor specimens and integrate them with in vivo cancer phenotypes.

Terminal phosphate labeled nucleotides: synthesis, applications, and linker effect on incorporation by DNA polymerases.

A number of terminal phosphate-labeled nucleotides with three or more phosphates and with varied length linkers attached between the terminal phosphate and the dye have been synthesized. These nucleotides have been tested as substrates for different DNA and RNA polymerases. We have also explored their utility in DNA sequencing, SNP analysis, nucleic acid amplification, quantitative PCR, and other biochemical assays.

Emerging understanding of multiscale tumor heterogeneity.

Cancer is a multifaceted disease characterized by heterogeneous genetic alterations and cellular metabolism, at the organ, tissue, and cellular level. Key features of cancer heterogeneity are summarized by 10 acquired capabilities, which govern malignant transformation and progression of invasive tumors. The relative contribution of these hallmark features to the disease process varies between cancers. At the DNA and cellular level, germ-line and somatic gene mutations are found across all cancer types, causing abnormal protein production, cell behavior, and growth. The tumor microenvironment and its individual components (immune cells, fibroblasts, collagen, and blood vessels) can also facilitate or restrict tumor growth and metastasis. Oncology research is currently in the midst of a tremendous surge of comprehension of these disease mechanisms. This will lead not only to novel drug targets but also to new challenges in drug discovery. Integrated, multi-omic, multiplexed technologies are essential tools in the quest to understand all of the various cellular changes involved in tumorigenesis. This review examines features of cancer heterogeneity and discusses how multiplexed technologies can facilitate a more comprehensive understanding of these features.

Terminal phosphate-labeled nucleotides with improved substrate properties for homogeneous nucleic acid assays.

Nucleotides with a dye attached to the terminal phosphate with four or more phosphates (tetra- or pentaphosphates) are superior substrates than the corresponding triphosphates for DNA and RNA polymerases. When fluorogenic dyes are directly attached to the terminal phosphate, they can be released by the action of polymerase and alkaline phosphatase. The released dye changes color and fluorescence properties. The fluorescent signal can also be amplified by using multiple labeled nucleotides to detect small amounts of template. We have explored the utility of these nucleotides in a variety of applications including homogeneous SNP detection methods, DNA sequencing, and quantitation of PCR and RCA.