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.

Tumor MHC Class I Expression Associates with Intralesional IL2 Response in Melanoma.

Cancer immunotherapy can result in lasting tumor regression, but predictive biomarkers of treatment response remain ill-defined. Here, we performed single-cell proteomics, transcriptomics, and genomics on matched untreated and IL2 injected metastases from patients with melanoma. Lesions that completely regressed following intralesional IL2 harbored increased fractions and densities of nonproliferating CD8+ T cells lacking expression of PD-1, LAG-3, and TIM-3 (PD-1-LAG-3-TIM-3-). Untreated lesions from patients who subsequently responded with complete eradication of all tumor cells in all injected lesions (individuals referred to herein as "extreme responders") were characterized by proliferating CD8+ T cells with an exhausted phenotype (PD-1+LAG-3+TIM-3+), stromal B-cell aggregates, and expression of IFNγ and IL2 response genes. Loss of membranous MHC class I expression in tumor cells of untreated lesions was associated with resistance to IL2 therapy. We validated this finding in an independent cohort of metastatic melanoma patients treated with intralesional or systemic IL2. Our study suggests that intact tumor-cell antigen presentation is required for melanoma response to IL2 and describes a multidimensional and spatial approach to develop immuno-oncology biomarker hypotheses using routinely collected clinical biospecimens.

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.

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.

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.

Microfluidic Tissue Mesodissection in Molecular Cancer Diagnostics.

We present a mesodissection platform that retains the advantages of laser-based dissection instrumentation with the speed and ease of manual dissection. Tissue dissection in clinical laboratories is often performed by manually scraping a physician-selected region from standard glass slide mounts. In this manner, costs associated with dissection remain low, but spatial resolution is compromised. In contrast, laser microdissection methods maintain spatial resolution that matches the requirements for analysis of important tissue heterogeneity but remains costly and labor intensive. We demonstrate a microfluidic tool for rapid extraction of histological regions of interest from formalin-fixed paraffin-embedded tissue, which uses a simple and automated method that is compatible with most downstream enzymatic reactions, including protocols used for next-generation DNA sequencing.

A novel, automated technology for multiplex biomarker imaging and application to breast cancer.

AIMS: Multiplexed immunofluorescence is a powerful tool for validating multigene assays and understanding the complex interplay of proteins implicated in breast cancer within a morphological context. We describe a novel technology for imaging an extended panel of biomarkers on a single, formalin-fixed paraffin-embedded breast sample and evaluating biomarker interaction at a single-cell level, and demonstrate proof-of-concept on a small set of breast tumours, including those which co-express hormone receptors with Her2/neu and Ki-67. METHODS AND RESULTS: Using a microfluidic flow cell, reagent exchange was automated and consisted of serial rounds of staining with dye-conjugated antibodies, imaging and chemical deactivation. A two-step antigen retrieval process was developed to satisfy all epitopes simultaneously, and key parameters were optimized. The imaging sequence was applied to seven breast tumours, and compared with conventional immunohistochemistry. Single-cell correlation analysis was performed with automated image processing. CONCLUSIONS: We have described a novel platform for evaluating biomarker co-localization. Expression in multiplexed images is consistent with conventional immunohistochemistry. Automation reduces inconsistencies in staining and positional shifts, while the fluorescent dye cycling approach dramatically expands the number of biomarkers which can be visualized and quantified on a single tissue section.

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.

Quantitative single cell analysis of cell population dynamics during submandibular salivary gland development and differentiation.

Epithelial organ morphogenesis involves reciprocal interactions between epithelial and mesenchymal cell types to balance progenitor cell retention and expansion with cell differentiation for evolution of tissue architecture. Underlying submandibular salivary gland branching morphogenesis is the regulated proliferation and differentiation of perhaps several progenitor cell populations, which have not been characterized throughout development, and yet are critical for understanding organ development, regeneration, and disease. Here we applied a serial multiplexed fluorescent immunohistochemistry technology to map the progressive refinement of the epithelial and mesenchymal cell populations throughout development from embryonic day 14 through postnatal day 20. Using computational single cell analysis methods, we simultaneously mapped the evolving temporal and spatial location of epithelial cells expressing subsets of differentiation and progenitor markers throughout salivary gland development. We mapped epithelial cell differentiation markers, including aquaporin 5, PSP, SABPA, and mucin 10 (acinar cells); cytokeratin 7 (ductal cells); and smooth muscle α-actin (myoepithelial cells) and epithelial progenitor cell markers, cytokeratin 5 and c-kit. We used pairwise correlation and visual mapping of the cells in multiplexed images to quantify the number of single- and double-positive cells expressing these differentiation and progenitor markers at each developmental stage. We identified smooth muscle α-actin as a putative early myoepithelial progenitor marker that is expressed in cytokeratin 5-negative cells. Additionally, our results reveal dynamic expansion and redistributions of c-kit- and K5-positive progenitor cell populations throughout development and in postnatal glands. The data suggest that there are temporally and spatially discreet progenitor populations that contribute to salivary gland development and homeostasis.

Simple and robust image-based autofocusing for digital microscopy.

A simple image-based autofocusing scheme for digital microscopy is demonstrated that uses as few as two intermediate images to bring the sample into focus. The algorithm is adapted to a commercial inverted microscope and used to automate brightfield and fluorescence imaging of histopathology tissue sections.

Long-working-distance incoherent-light interference microscope.

We describe the design and operation of a long-working-distance, incoherent light interference microscope that has been developed to address the growing demand for new microsystem characterization tools. The design of the new microscope is similar to that of a Linnik interference microscope and thus preserves the full working distance of the long-working-distance objectives utilized. However, in contrast to a traditional Linnik microscope, the new microscope does not rely on the use of matched objectives in the sample and the reference arms of the interferometer. An adjustable optical configuration has been devised that allows the total optical path length, wavefront curvature, and dispersion of the reference arm to be matched to the sample arm of the interferometer. The reference arm configuration can be adjusted to provide matching for 5x, 10x, and 20x long-working-distance objectives in the sample arm. In addition to retaining the full working distance of the sample arm objectives, the new design allows interference images to be acquired in situations in which intervening windows are necessary, such as occur with packaged microsystems, microfluidic devices, and cryogenic, vacuum, or environmental chamber studies of microsystem performance. The interference microscope is compatible with phase-shifting interferometry, vertical scanning interferometry, and stroboscopic measurement of dynamic processes.