Advances in multiplexed imaging technologies have drastically improved our ability to characterize healthy and diseased tissues at the single-cell level. Co-detection by indexing (CODEX) relies on ...DNA-conjugated antibodies and the cyclic addition and removal of complementary fluorescently labeled DNA probes and has been used so far to simultaneously visualize up to 60 markers in situ. CODEX enables a deep view into the single-cell spatial relationships in tissues and is intended to spur discovery in developmental biology, disease and therapeutic design. Herein, we provide optimized protocols for conjugating purified antibodies to DNA oligonucleotides, validating the conjugation by CODEX staining and executing the CODEX multicycle imaging procedure for both formalin-fixed, paraffin-embedded (FFPE) and fresh-frozen tissues. In addition, we describe basic image processing and data analysis procedures. We apply this approach to an FFPE human tonsil multicycle experiment. The hands-on experimental time for antibody conjugation is ~4.5 h, validation of DNA-conjugated antibodies with CODEX staining takes ~6.5 h and preparation for a CODEX multicycle experiment takes ~8 h. The multicycle imaging and data analysis time depends on the tissue size, number of markers in the panel and computational complexity.
Cutaneous T cell lymphomas (CTCL) are rare but aggressive cancers without effective treatments. While a subset of patients derive benefit from PD-1 blockade, there is a critically unmet need for ...predictive biomarkers of response. Herein, we perform CODEX multiplexed tissue imaging and RNA sequencing on 70 tumor regions from 14 advanced CTCL patients enrolled in a pembrolizumab clinical trial (NCT02243579). We find no differences in the frequencies of immune or tumor cells between responders and non-responders. Instead, we identify topographical differences between effector PD-1
CD4
T cells, tumor cells, and immunosuppressive Tregs, from which we derive a spatial biomarker, termed the SpatialScore, that correlates strongly with pembrolizumab response in CTCL. The SpatialScore coincides with differences in the functional immune state of the tumor microenvironment, T cell function, and tumor cell-specific chemokine recruitment and is validated using a simplified, clinically accessible tissue imaging platform. Collectively, these results provide a paradigm for investigating the spatial balance of effector and suppressive T cell activity and broadly leveraging this biomarker approach to inform the clinical use of immunotherapies.
Antitumoral immunity requires organized, spatially nuanced interactions between components of the immune tumor microenvironment (iTME). Understanding this coordinated behavior in effective versus ...ineffective tumor control will advance immunotherapies. We re-engineered co-detection by indexing (CODEX) for paraffin-embedded tissue microarrays, enabling simultaneous profiling of 140 tissue regions from 35 advanced-stage colorectal cancer (CRC) patients with 56 protein markers. We identified nine conserved, distinct cellular neighborhoods (CNs)—a collection of components characteristic of the CRC iTME. Enrichment of PD-1+CD4+ T cells only within a granulocyte CN positively correlated with survival in a high-risk patient subset. Coupling of tumor and immune CNs, fragmentation of T cell and macrophage CNs, and disruption of inter-CN communication was associated with inferior outcomes. This study provides a framework for interrogating how complex biological processes, such as antitumoral immunity, occur through concerted actions of cells and spatial domains.
Display omitted
•FFPE-CODEX multiplexed tissue imaging of 56 markers in 140 tissues of 35 CRC patients•Cellular neighborhoods reveal spatial organization of the tumor microenvironment•Altered organization of tumor and immune components in low- versus high-risk patients•Local enrichment of PD-1+CD4+ T cells correlates with survival in high-risk patients
A multiplexed tissue imaging and computational analysis framework applied to colorectal cancer allows interrogation of how spatial organization of the immune tumor microenvironment is linked to clinical outcomes.
Multiplexed imaging is a recently developed and powerful single-cell biology research tool. However, it presents new sources of technical noise that are distinct from other types of single-cell data, ...necessitating new practices for single-cell multiplexed imaging processing and analysis, particularly regarding cell-type identification. Here we created single-cell multiplexed imaging datasets by performing CODEX on four sections of the human colon (ascending, transverse, descending, and sigmoid) using a panel of 47 oligonucleotide-barcoded antibodies. After cell segmentation, we implemented five different normalization techniques crossed with four unsupervised clustering algorithms, resulting in 20 unique cell-type annotations for the same dataset. We generated two standard annotations: hand-gated cell types and cell types produced by over-clustering with spatial verification. We then compared these annotations at four levels of cell-type granularity. First, increasing cell-type granularity led to decreased labeling accuracy; therefore, subtle phenotype annotations should be avoided at the clustering step. Second, accuracy in cell-type identification varied more with normalization choice than with clustering algorithm. Third, unsupervised clustering better accounted for segmentation noise during cell-type annotation than hand-gating. Fourth, Z-score normalization was generally effective in mitigating the effects of noise from single-cell multiplexed imaging. Variation in cell-type identification will lead to significant differential spatial results such as cellular neighborhood analysis; consequently, we also make recommendations for accurately assigning cell-type labels to CODEX multiplexed imaging.
Recent identification of platelet/megakaryocyte-biased hematopoietic stem/repopulating cells requires revision of the intermediate pathway for megakaryopoiesis. Here, we show a unipotent ...megakaryopoietic pathway bypassing the bipotent megakaryocyte/erythroid progenitors (biEMPs). Cells purified from mouse bone marrow by CD42b (GPIbα) marking were demonstrated to be unipotent megakaryocytic progenitors (MKPs) by culture and transplantation. A subpopulation of freshly isolated CD41(+) cells in the lineage Sca1(+) cKit(+) (LSK) fraction (subCD41(+) LSK) differentiated only into MKP and mature megakaryocytes in culture. Although CD41(+) LSK cells as a whole were capable of differentiating into all myeloid and lymphoid cells in vivo, they produced unipotent MKP, mature megakaryocytes, and platelets in vitro and in vivo much more efficiently than Flt3(+) CD41(-) LSK cells, especially at the early phase after transplantation. In single cell polymerase chain reaction and thrombopoietin (TPO) signaling analyses, the MKP and a fraction of CD41(+) LSK, but not the biEMP, showed the similarities in mRNA expression profile and visible TPO-mediated phosphorylation. On increased demand of platelet production after 5-FU treatment, a part of CD41(+) LSK population expressed CD42b on the surface, and 90% of them showed unipotent megakaryopoietic capacity in single cell culture and predominantly produced platelets in vivo at the early phase after transplantation. These results suggest that the CD41(+) CD42b(+) LSK are straightforward progenies of megakaryocytes/platelet-biased stem/repopulating cells, but not progenies of biEMP. Consequently, we show a unipotent/highly biased megakaryopoietic pathway interconnecting stem/repopulating cells and mature megakaryocytes, the one that may play physiologic roles especially in emergency megakaryopoiesis.
Antigen-specific T cells traffic to, are influenced by, and create unique cellular microenvironments. Here we characterize these microenvironments over time with multiplexed imaging in a melanoma ...model of adoptive T cell therapy and human patients with melanoma treated with checkpoint inhibitor therapy. Multicellular neighborhood analysis reveals dynamic immune cell infiltration and inflamed tumor cell neighborhoods associated with CD8+ T cells. T cell-focused analysis indicates T cells are found along a continuum of neighborhoods that reflect the progressive steps coordinating the anti-tumor immune response. More effective anti-tumor immune responses are characterized by inflamed tumor-T cell neighborhoods, flanked by dense immune infiltration neighborhoods. Conversely, ineffective T cell therapies express anti-inflammatory cytokines, resulting in regulatory neighborhoods, spatially disrupting productive T cell-immune and -tumor interactions. Our study provides in situ mechanistic insights into temporal tumor microenvironment changes, cell interactions critical for response, and spatial correlates of immunotherapy outcomes, informing cellular therapy evaluation and engineering.
Display omitted
•CODEX multiplexed imaging of human and mouse, immunotherapy-treated melanomas•Tumor-specific T cells cause inflamed tumor phenotype and unique cell neighborhoods•T cell phenotype influences capacity to make effective tumor-killing neighborhoods•Potent therapy creates larger, linked areas of inflamed tumor and immune infiltrate
Hickey et al. discovered that tumor-specific T cells induce tumor inflammation and establish progressive inflammatory environments. Effective therapies exhibit compartmentalized immune infiltration, co-located with T cell and inflamed tumor zones, emphasizing the need for therapeutic T cells to be designed to curate local and distal microenvironments.
ALDH2 inactivating mutation (ALDH2*2) is the most abundant mutation leading to bone morphological aberration. Osteoporosis has long been associated with changes in bone biomaterial in elderly ...populations. Such changes can be exacerbated with elevated ethanol consumption and in subjects with impaired ethanol metabolism, such as carriers of aldehyde dehydrogenase 2 (ALDH2)-deficient gene, ALDH2*2. So far, little is known about bone compositional changes besides a decrease in mineralization. Raman spectroscopic imaging has been utilized to study the changes in overall composition of C57BL/6 female femur bone sections, as well as in compound spatial distribution. Raman maps of bone sections were analyzed using multilinear regression with these four isolated components, resulting in maps of their relative distribution. A 15-week treatment of both wild-type (WT) and ALDH2*2/*2 mice with 20% ethanol in the drinking water resulted in a significantly lower mineral content (
p
< 0.05) in the bones. There was no significant change in mineral and collagen content due to the mutation alone (
p
> 0.4). Highly localized islets of elongated adipose tissue were observed on most maps. Elevated fat content was found in ALDH2*2 knock-in mice consuming ethanol (
p
< 0.0001) and this effect appeared cumulative. This work conclusively demonstrates that that osteocytes in femurs of older female mice accumulate fat, as has been previously theorized, and that fat accumulation is likely modulated by levels of acetaldehyde, the ethanol metabolite.
Highly multiplexed, single-cell imaging has revolutionized our understanding of spatial cellular interactions associated with health and disease. With ever-increasing numbers of antigens, region ...sizes, and sample sizes, multiplexed fluorescence imaging experiments routinely produce terabytes of data. Fast and accurate processing of these large-scale, high-dimensional imaging data is essential to ensure reliable segmentation and identification of cell types and for characterization of cellular neighborhoods and inference of mechanistic insights. Here, we describe RAPID, a Real-time, GPU-Accelerated Parallelized Image processing software for large-scale multiplexed fluorescence microscopy Data. RAPID deconvolves large-scale, high-dimensional fluorescence imaging data, stitches and registers images with axial and lateral drift correction, and minimizes tissue autofluorescence such as that introduced by erythrocytes. Incorporation of an open source CUDA-driven, GPU-assisted deconvolution produced results similar to fee-based commercial software. RAPID reduces data processing time and artifacts and improves image contrast and signal-to-noise compared to our previous image processing pipeline, thus providing a useful tool for accurate and robust analysis of large-scale, multiplexed, fluorescence imaging data.
Dorsoventral (DV) patterning of the Drosophila embryo is controlled by a concentration gradient of Dorsal, a sequence-specific transcription factor related to mammalian NF-κB. The Dorsal gradient ...generates at least 3 distinct thresholds of gene activity and tissue specification by the differential regulation of target enhancers containing distinctive combinations of binding sites for Dorsal, Twist, Snail, and other DV determinants. To understand the evolution of DV patterning mechanisms, we identified and characterized Dorsal target enhancers from the mosquito Anopheles gambiae and the flour beetle Tribolium castaneum. Putative orthologous enhancers are located in similar positions relative to the target genes they control, even though they lack sequence conservation and sometimes produce divergent patterns of gene expression. The most dramatic example of this conservation is seen for the "shadow" enhancer regulating brinker: It is conserved within the intron of the neighboring Atg5 locus of both flies and mosquitoes. These results suggest that, like exons, an enhancer position might be subject to constraint. Thus, novel patterns of gene expression might arise from the modification of conserved enhancers rather than the invention of new ones. We propose that this enhancer constancy might be a general property of regulatory evolution, and should facilitate enhancer discovery in nonmodel organisms.
During the evolution of the Diptera there is a dramatic modification of the embryonic ectoderm, whereby mosquitoes contain separate amnion and serosa lineages while higher flies such as
Drosophila ...melanogaster contain a single amnioserosa. Whole-genome transcriptome assays were performed with isolated serosa from
Anopheles gambiae embryos. These assays identified a large number of genes implicated in the production of the larval cuticle. In
D. melanogaster, these genes are activated just once during embryogenesis, during late stages where they are used for the production of the larval cuticle. Evidence is presented that the serosal cells secrete a dedicated serosal cuticle, which protects
A. gambiae embryos from desiccation. Detailed temporal microarray assays of mosquito gene expression profiles revealed that the cuticular genes display biphasic expression during
A. gambiae embryogenesis, first in the serosa of early embryos and then again during late stages as seen in
D. melanogaster. We discuss how evolutionary modifications in the well-defined dorsal–ventral patterning network led to the wholesale deployment of the cuticle biosynthesis pathway in early embryos of
A. gambiae.