GeneSurfer Enables Transcriptome-wide Exploration and Functional Annotation of Gene Co-expression Modules in 3D Spatial Transcriptomics Data.

Gene co-expression provides crucial insights into biological functions, however, there is a lack of exploratory analysis tools for localized gene co-expression in large-scale datasets. We present GeneSurfer, an interactive interface designed to explore localized transcriptome-wide gene co-expression patterns in the 3D spatial domain. Key features of GeneSurfer include transcriptome-wide gene filtering and gene clustering based on spatial local co-expression within transcriptomically similar cells, multi-slice 3D rendering of average expression of gene clusters, and on-the-fly Gene Ontology term annotation of co-expressed gene sets. Additionally, GeneSurfer offers multiple linked views for investigating individual genes or gene co-expression in the spatial domain at each exploration stage. Demonstrating its utility with both spatial transcriptomics and single-cell RNA sequencing data from the Allen Brain Cell Atlas, GeneSurfer effectively identifies and annotates localized transcriptome-wide co-expression, providing biological insights and facilitating hypothesis generation and validation.

Accelerating Hyperbolic t-SNE.

The need to understand the structure of hierarchical or high-dimensional data is present in a variety of fields. Hyperbolic spaces have proven to be an important tool for embedding computations and analysis tasks as their non-linear nature lends itself well to tree or graph data. Subsequently, they have also been used in the visualization of high-dimensional data, where they exhibit increased embedding performance. However, none of the existing dimensionality reduction methods for embedding into hyperbolic spaces scale well with the size of the input data. That is because the embeddings are computed via iterative optimization schemes and the computation cost of every iteration is quadratic in the size of the input. Furthermore, due to the non-linear nature of hyperbolic spaces, euclidean acceleration structures cannot directly be translated to the hyperbolic setting. This article introduces the first acceleration structure for hyperbolic embeddings, building upon a polar quadtree. We compare our approach with existing methods and demonstrate that it computes embeddings of similar quality in significantly less time. Implementation and scripts for the experiments can be found at https://graphics.tudelft.nl/accelerating-hyperbolic-tsne.

ManiVault: A Flexible and Extensible Visual Analytics Framework for High-Dimensional Data.

Exploration and analysis of high-dimensional data are important tasks in many fields that produce large and complex data, like the financial sector, systems biology, or cultural heritage. Tailor-made visual analytics software is developed for each specific application, limiting their applicability in other fields. However, as diverse as these fields are, their characteristics and requirements for data analysis are conceptually similar. Many applications share abstract tasks and data types and are often constructed with similar building blocks. Developing such applications, even when based mostly on existing building blocks, requires significant engineering efforts. We developed ManiVault, a flexible and extensible open-source visual analytics framework for analyzing high-dimensional data. The primary objective of ManiVault is to facilitate rapid prototyping of visual analytics workflows for visualization software developers and practitioners alike. ManiVault is built using a plugin-based architecture that offers easy extensibility. While our architecture deliberately keeps plugins self-contained, to guarantee maximum flexibility and re-usability, we have designed and implemented a messaging API for tight integration and linking of modules to support common visual analytics design patterns. We provide several visualization and analytics plugins, and ManiVault's API makes the integration of new plugins easy for developers. ManiVault facilitates the distribution of visualization and analysis pipelines and results for practitioners through saving and reproducing complete application states. As such, ManiVault can be used as a communication tool among researchers to discuss workflows and results. A copy of this paper and all supplemental material is available at osf.io/9k6jw, and source code at github.com/ManiVaultStudio.

SpaceWalker enables interactive gradient exploration for spatial transcriptomics data.

In spatial transcriptomics (ST) data, biologically relevant features such as tissue compartments or cell-state transitions are reflected by gene expression gradients. Here, we present SpaceWalker, a visual analytics tool for exploring the local gradient structure of 2D and 3D ST data. The user can be guided by the local intrinsic dimensionality of the high-dimensional data to define seed locations, from which a flood-fill algorithm identifies transcriptomically similar cells on the fly, based on the high-dimensional data topology. In several use cases, we demonstrate that the spatial projection of these flooded cells highlights tissue architectural features and that interactive retrieval of gene expression gradients in the spatial and transcriptomic domains confirms known biology. We also show that SpaceWalker generalizes to several different ST protocols and scales well to large, multi-slice, 3D whole-brain ST data while maintaining real-time interaction performance.

Comparative transcriptomics reveals human-specific cortical features.

The cognitive abilities of humans are distinctive among primates, but their molecular and cellular substrates are poorly understood. We used comparative single-nucleus transcriptomics to analyze samples of the middle temporal gyrus (MTG) from adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets to understand human-specific features of the neocortex. Human, chimpanzee, and gorilla MTG showed highly similar cell-type composition and laminar organization as well as a large shift in proportions of deep-layer intratelencephalic-projecting neurons compared with macaque and marmoset MTG. Microglia, astrocytes, and oligodendrocytes had more-divergent expression across species compared with neurons or oligodendrocyte precursor cells, and neuronal expression diverged more rapidly on the human lineage. Only a few hundred genes showed human-specific patterning, suggesting that relatively few cellular and molecular changes distinctively define adult human cortical structure.

Identification of a Disease-Associated Network of Intestinal Immune Cells in Treatment-Naive Inflammatory Bowel Disease.

Chronic intestinal inflammation underlies inflammatory bowel disease (IBD). Previous studies indicated alterations in the cellular immune system; however, it has been challenging to interrogate the role of all immune cell subsets simultaneously. Therefore, we aimed to identify immune cell types associated with inflammation in IBD using high-dimensional mass cytometry. We analyzed 188 intestinal biopsies and paired blood samples of newly-diagnosed, treatment-naive patients (n=42) and controls (n=26) in two independent cohorts. We applied mass cytometry (36-antibody panel) to resolve single cells and analyzed the data with unbiased Hierarchical-SNE. In addition, imaging-mass cytometry (IMC) was performed to reveal the spatial distribution of the immune subsets in the tissue. We identified 44 distinct immune subsets. Correlation network analysis identified a network of inflammation-associated subsets, including HLA-DR+CD38+ EM CD4+ T cells, T regulatory-like cells, PD1+ EM CD8+ T cells, neutrophils, CD27+ TCRγδ cells and NK cells. All disease-associated subsets were validated in a second cohort. This network was abundant in a subset of patients, independent of IBD subtype, severity or intestinal location. Putative disease-associated CD4+ T cells were detectable in blood. Finally, imaging-mass cytometry revealed the spatial colocalization of neutrophils, memory CD4+ T cells and myeloid cells in the inflamed intestine. Our study indicates that a cellular network of both innate and adaptive immune cells colocalizes in inflamed biopsies from a subset of patients. These results contribute to dissecting disease heterogeneity and may guide the development of targeted therapeutics in IBD.

Co-expression patterns of microglia markers Iba1, TMEM119 and P2RY12 in Alzheimer's disease.

Microglia have been identified as key players in Alzheimer's disease pathogenesis, and other neurodegenerative diseases. Iba1, and more specifically TMEM119 and P2RY12 are gaining ground as presumedly more specific microglia markers, but comprehensive characterization of the expression of these three markers individually as well as combined is currently missing. Here we used a multispectral immunofluorescence dataset, in which over seventy thousand microglia from both aged controls and Alzheimer patients have been analysed for expression of Iba1, TMEM119 and P2RY12 on a single-cell level. For all markers, we studied the overlap and differences in expression patterns and the effect of proximity to ß-amyloid plaques. We found no difference in absolute microglia numbers between control and Alzheimer subjects, but the prevalence of specific combinations of markers (phenotypes) differed greatly. In controls, the majority of microglia expressed all three markers. In Alzheimer patients, a significant loss of TMEM119+-phenotypes was observed, independent of the presence of ß-amyloid plaques in its proximity. Contrary, phenotypes showing loss of P2RY12, but consistent Iba1 expression were increasingly prevalent around ß-amyloid plaques. No morphological features were conclusively associated with loss or gain of any of the markers or any of the identified phenotypes. All in all, none of the three markers were expressed by all microglia, nor can be wholly regarded as a pan- or homeostatic marker, and preferential phenotypes were observed depending on the surrounding pathological or homeostatic environment. This work could help select and interpret microglia markers in previous and future studies.

Tumor-specific T cells support chemokine-driven spatial organization of intratumoral immune microaggregates needed for long survival.

BACKGROUND: The composition of the tumor immune microenvironment (TIME) associated with good prognosis generally also predicts the success of immunotherapy, and both entail the presence of pre-existing tumor-specific T cells. Here, the blueprint of the TIME associated with such an ongoing tumor-specific T-cell response was dissected in a unique prospective oropharyngeal squamous cell carcinoma (OPSCC) cohort, in which tumor-specific tumor-infiltrating T cells were detected (immune responsiveness (IR+)) or not (lack of immune responsiveness (IR-)). METHODS: A comprehensive multimodal, high-dimensional strategy was applied to dissect the TIME of treatment-naive IR+ and IR- OPSCC tissue, including bulk RNA sequencing (NanoString), imaging mass cytometry (Hyperion) for phenotyping and spatial interaction analyses of immune cells, and combined single-cell gene expression profiling and T-cell receptor (TCR) sequencing (single-cell RNA sequencing (scRNAseq)) to characterize the transcriptional states of clonally expanded tumor-infiltrating T cells. RESULTS: IR+ patients had an excellent survival during >10 years follow-up. The tumors of IR+ patients expressed higher levels of genes strongly related to interferon gamma signaling, T-cell activation, TCR signaling, and mononuclear cell differentiation, as well as genes involved in several immune signaling pathways, than IR- patients. The top differently overexpressed genes included CXCL12 and LTB, involved in ectopic lymphoid structure development. Moreover, scRNAseq not only revealed that CD4+ T cells were the main producers of LTB but also identified a subset of clonally expanded CD8+ T cells, dominantly present in IR+ tumors, which secreted the T cell and dendritic cell (DC) attracting chemokine CCL4. Indeed, immune cell infiltration in IR+ tumors is stronger, highly coordinated, and has a distinct spatial phenotypical signature characterized by intratumoral microaggregates of CD8+CD103+ and CD4+ T cells with DCs. In contrast, the IR- TIME comprised spatial interactions between lymphocytes and various immunosuppressive myeloid cell populations. The impact of these chemokines on local immunity and clinical outcome was confirmed in an independent The Cancer Genome Atlas OPSCC cohort. CONCLUSION: The production of lymphoid cell attracting and organizing chemokines by tumor-specific T cells in IR+ tumors constitutes a positive feedback loop to sustain the formation of the DC-T-cell microaggregates and identifies patients with excellent survival after standard therapy.

Semi-automated background removal limits data loss and normalizes imaging mass cytometry data.

Imaging mass cytometry (IMC) allows the detection of multiple antigens (approximately 40 markers) combined with spatial information, making it a unique tool for the evaluation of complex biological systems. Due to its widespread availability and retained tissue morphology, formalin-fixed, paraffin-embedded (FFPE) tissues are often a material of choice for IMC studies. However, antibody performance and signal to noise ratios can differ considerably between FFPE tissues as a consequence of variations in tissue processing, including fixation. In contrast to batch effects caused by differences in the immunodetection procedure, variations in tissue processing are difficult to control. We investigated the effect of immunodetection-related signal intensity fluctuations on IMC analysis and phenotype identification, in a cohort of 12 colorectal cancer tissues. Furthermore, we explored different normalization strategies and propose a workflow to normalize IMC data by semi-automated background removal, using publicly available tools. This workflow can be directly applied to previously acquired datasets and considerably improves the quality of IMC data, thereby supporting the analysis and comparison of multiple samples.

Systems analysis and controlled malaria infection in Europeans and Africans elucidate naturally acquired immunity.

Controlled human infections provide opportunities to study the interaction between the immune system and malaria parasites, which is essential for vaccine development. Here, we compared immune signatures of malaria-naive Europeans and of Africans with lifelong malaria exposure using mass cytometry, RNA sequencing and data integration, before and 5 and 11 days after venous inoculation with Plasmodium falciparum sporozoites. We observed differences in immune cell populations, antigen-specific responses and gene expression profiles between Europeans and Africans and among Africans with differing degrees of immunity. Before inoculation, an activated/differentiated state of both innate and adaptive cells, including elevated CD161+CD4+ T cells and interferon-γ production, predicted Africans capable of controlling parasitemia. After inoculation, the rapidity of the transcriptional response and clusters of CD4+ T cells, plasmacytoid dendritic cells and innate T cells were among the features distinguishing Africans capable of controlling parasitemia from susceptible individuals. These findings can guide the development of a vaccine effective in malaria-endemic regions.

Iron loading is a prominent feature of activated microglia in Alzheimer's disease patients.

Brain iron accumulation has been found to accelerate disease progression in amyloid-ß(Aß) positive Alzheimer patients, though the mechanism is still unknown. Microglia have been identified as key players in the disease pathogenesis, and are highly reactive cells responding to aberrations such as increased iron levels. Therefore, using histological methods, multispectral immunofluorescence and an automated in-house developed microglia segmentation and analysis pipeline, we studied the occurrence of iron-accumulating microglia and the effect on its activation state in human Alzheimer brains. We identified a subset of microglia with increased expression of the iron storage protein ferritin light chain (FTL), together with increased Iba1 expression, decreased TMEM119 and P2RY12 expression. This activated microglia subset represented iron-accumulating microglia and appeared morphologically dystrophic. Multispectral immunofluorescence allowed for spatial analysis of FTL+Iba1+-microglia, which were found to be the predominant Aß-plaque infiltrating microglia. Finally, an increase of FTL+Iba1+-microglia was seen in patients with high Aß load and Tau load. These findings suggest iron to be taken up by microglia and to influence the functional phenotype of these cells, especially in conjunction with Aß.

ImaCytE: Visual Exploration of Cellular Micro-Environments for Imaging Mass Cytometry Data.

Tissue functionality is determined by the characteristics of tissue-resident cells and their interactions within their microenvironment. Imaging Mass Cytometry offers the opportunity to distinguish cell types with high precision and link them to their spatial location in intact tissues at sub-cellular resolution. This technology produces large amounts of spatially-resolved high-dimensional data, which constitutes a serious challenge for the data analysis. We present an interactive visual analysis workflow for the end-to-end analysis of Imaging Mass Cytometry data that was developed in close collaboration with domain expert partners. We implemented the presented workflow in an interactive visual analysis tool; ImaCytE. Our workflow is designed to allow the user to discriminate cell types according to their protein expression profiles and analyze their cellular microenvironments, aiding in the formulation or verification of hypotheses on tissue architecture and function. Finally, we show the effectiveness of our workflow and ImaCytE through a case study performed by a collaborating specialist.

InCorr: Interactive Data-Driven Correlation Panels for Digital Outcrop Analysis.

Geological analysis of 3D Digital Outcrop Models (DOMs) for reconstruction of ancient habitable environments is a key aspect of the upcoming ESA ExoMars 2022 Rosalind Franklin Rover and the NASA 2020 Rover Perseverance missions in seeking signs of past life on Mars. Geologists measure and interpret 3D DOMs, create sedimentary logs and combine them in 'correlation panels' to map the extents of key geological horizons, and build a stratigraphic model to understand their position in the ancient landscape. Currently, the creation of correlation panels is completely manual and therefore time-consuming, and inflexible. With InCorr we present a visualization solution that encompasses a 3D logging tool and an interactive data-driven correlation panel that evolves with the stratigraphic analysis. For the creation of InCorr we closely cooperated with leading planetary geologists in the form of a design study. We verify our results by recreating an existing correlation analysis with InCorr and validate our correlation panel against a manually created illustration. Further, we conducted a user-study with a wider circle of geologists. Our evaluation shows that InCorr efficiently supports the domain experts in tackling their research questions and that it has the potential to significantly impact how geologists work with digital outcrop representations in general.

Visual cohort comparison for spatial single-cell omics-data.

Spatially-resolved omics-data enable researchers to precisely distinguish cell types in tissue and explore their spatial interactions, enabling deep understanding of tissue functionality. To understand what causes or deteriorates a disease and identify related biomarkers, clinical researchers regularly perform large-scale cohort studies, requiring the comparison of such data at cellular level. In such studies, with little a-priori knowledge of what to expect in the data, explorative data analysis is a necessity. Here, we present an interactive visual analysis workflow for the comparison of cohorts of spatially-resolved omics-data. Our workflow allows the comparative analysis of two cohorts based on multiple levels-of-detail, from simple abundance of contained cell types over complex co-localization patterns to individual comparison of complete tissue images. As a result, the workflow enables the identification of cohort-differentiating features, as well as outlier samples at any stage of the workflow. During the development of the workflow, we continuously consulted with domain experts. To show the effectiveness of the workflow, we conducted multiple case studies with domain experts from different application areas and with different data modalities.

Visualizing Dynamic Changes at the Maternal-Fetal Interface Throughout Human Pregnancy by Mass Cytometry.

During healthy pregnancy, a balanced microenvironment at the maternal-fetal interface with coordinated interaction between various immune cells is necessary to maintain immunological tolerance. While specific decidual immune cell subsets have been investigated, a system-wide unbiased approach is lacking. Here, mass cytometry was applied for data-driven, in-depth immune profiling of the total leukocyte population isolated from first, second, and third trimester decidua, as well as maternal peripheral blood at time of delivery. The maternal-fetal interface showed a unique composition of immune cells, different from peripheral blood, with significant differences between early and term pregnancy samples. Profiling revealed substantial heterogeneity in the decidual lymphoid and myeloid cell lineages that shape gestational-specific immune networks and putative differentiation trajectories over time during gestation. Uncovering the overall complexity at the maternal-fetal interface throughout pregnancy resulted in a human atlas that may serve as a foundation upon which comprehension of the immune microenvironment and alterations thereof in pregnancy complications can be built.

Helminth infections drive heterogeneity in human type 2 and regulatory cells.

Helminth infections induce strong type 2 and regulatory responses, but the degree of heterogeneity of such cells is not well characterized. Using mass cytometry, we profiled these cells in Europeans and Indonesians not exposed to helminths and in Indonesians residing in rural areas infected with soil-transmitted helminths. To assign immune alteration to helminth infection, the profiling was performed before and 1 year after deworming. Very distinct signatures were found in Europeans and Indonesians, showing expanded frequencies of T helper 2 cells, particularly CD161+ cells and ILC2s in helminth-infected Indonesians, which was confirmed functionally through analysis of cytokine-producing cells. Besides ILC2s and CD4+ T cells, CD8+ T cells and γδ T cells in Indonesians produced type 2 cytokines. Regulatory T cells were also expanded in Indonesians, but only those expressing CTLA-4, and some coexpressed CD38, HLA-DR, ICOS, or CD161. CD11c+ B cells were found to be the main IL-10 producers among B cells in Indonesians, a subset that was almost absent in Europeans. A number of the distinct immune profiles were driven by helminths as the profiles reverted after clearance of helminth infections. Moreover, Indonesians with no helminth infections residing in an urban area showed immune profiles that resembled Europeans rather than rural Indonesians, which excludes a major role for ethnicity. Detailed insight into the human type 2 and regulatory networks could provide opportunities to target these cells for more precise interventions.

High-dimensional cytometric analysis of colorectal cancer reveals novel mediators of antitumour immunity.

OBJECTIVE: A comprehensive understanding of anticancer immune responses is paramount for the optimal application and development of cancer immunotherapies. We unravelled local and systemic immune profiles in patients with colorectal cancer (CRC) by high-dimensional analysis to provide an unbiased characterisation of the immune contexture of CRC. DESIGN: Thirty-six immune cell markers were simultaneously assessed at the single-cell level by mass cytometry in 35 CRC tissues, 26 tumour-associated lymph nodes, 17 colorectal healthy mucosa and 19 peripheral blood samples from 31 patients with CRC. Additionally, functional, transcriptional and spatial analyses of tumour-infiltrating lymphocytes were performed by flow cytometry, single-cell RNA-sequencing and multispectral immunofluorescence. RESULTS: We discovered that a previously unappreciated innate lymphocyte population (Lin-CD7+CD127-CD56+CD45RO+) was enriched in CRC tissues and displayed cytotoxic activity. This subset demonstrated a tissue-resident (CD103+CD69+) phenotype and was most abundant in immunogenic mismatch repair (MMR)-deficient CRCs. Their presence in tumours was correlated with the infiltration of tumour-resident cytotoxic, helper and γδ T cells with highly similar activated (HLA-DR+CD38+PD-1+) phenotypes. Remarkably, activated γδ T cells were almost exclusively found in MMR-deficient cancers. Non-activated counterparts of tumour-resident cytotoxic and γδ T cells were present in CRC and healthy mucosa tissues, but not in lymph nodes, with the exception of tumour-positive lymph nodes. CONCLUSION: This work provides a blueprint for the understanding of the heterogeneous and intricate immune landscape of CRC, including the identification of previously unappreciated immune cell subsets. The concomitant presence of tumour-resident innate and adaptive immune cell populations suggests a multitargeted exploitation of their antitumour properties in a therapeutic setting.

GPGPU Linear Complexity t-SNE Optimization.

In recent years the t-distributed Stochastic Neighbor Embedding (t-SNE) algorithm has become one of the most used and insightful techniques for exploratory data analysis of high-dimensional data. It reveals clusters of high-dimensional data points at different scales while only requiring minimal tuning of its parameters. However, the computational complexity of the algorithm limits its application to relatively small datasets. To address this problem, several evolutions of t-SNE have been developed in recent years, mainly focusing on the scalability of the similarity computations between data points. However, these contributions are insufficient to achieve interactive rates when visualizing the evolution of the t-SNE embedding for large datasets. In this work, we present a novel approach to the minimization of the t-SNE objective function that heavily relies on graphics hardware and has linear computational complexity. Our technique decreases the computational cost of running t-SNE on datasets by orders of magnitude and retains or improves on the accuracy of past approximated techniques. We propose to approximate the repulsive forces between data points by splatting kernel textures for each data point. This approximation allows us to reformulate the t-SNE minimization problem as a series of tensor operations that can be efficiently executed on the graphics card. An efficient implementation of our technique is integrated and available for use in the widely used Google TensorFlow.js, and an open-source C++ library.

Multidimensional analyses of proinsulin peptide-specific regulatory T cells induced by tolerogenic dendritic cells.

Induction of antigen-specific regulatory T cells (Tregs) in vivo is the holy grail of current immune-regulating therapies in autoimmune diseases, such as type 1 diabetes. Tolerogenic dendritic cells (tolDCs) generated from monocytes by a combined treatment with vitamin D and dexamethasone (marked by CD52hi and CD86lo expression) induce antigen-specific Tregs. We evaluated the phenotypes of these Tregs using high-dimensional mass cytometry to identify a surface-based T cell signature of tolerogenic modulation. Naïve CD4+ T cells were stimulated with tolDCs or mature inflammatory DCs pulsed with proinsulin peptide, after which the suppressive capacity, cytokine production and phenotype of stimulated T cells were analysed. TolDCs induced suppressive T cell lines that were dominated by a naïve phenotype (CD45RA+CCR7+). These naïve T cells, however, did not show suppressive capacity, but were arrested in their naïve status. T cell cultures stimulated by tolDC further contained memory-like (CD45RA-CCR7-) T cells expressing regulatory markers Lag-3, CD161 and ICOS. T cells expressing CD25lo or CD25hi were most prominent and suppressed CD4+ proliferation, while CD25hi Tregs also effectively supressed effector CD8+ T cells. We conclude that tolDCs induce antigen-specific Tregs with various phenotypes. This extends our earlier findings pointing to a functionally diverse pool of antigen-induced and specific Tregs and provides the basis for immune-monitoring in clinical trials with tolDC.

Conserved cell types with divergent features in human versus mouse cortex.

Elucidating the cellular architecture of the human cerebral cortex is central to understanding our cognitive abilities and susceptibility to disease. Here we used single-nucleus RNA-sequencing analysis to perform a comprehensive study of cell types in the middle temporal gyrus of human cortex. We identified a highly diverse set of excitatory and inhibitory neuron types that are mostly sparse, with excitatory types being less layer-restricted than expected. Comparison to similar mouse cortex single-cell RNA-sequencing datasets revealed a surprisingly well-conserved cellular architecture that enables matching of homologous types and predictions of properties of human cell types. Despite this general conservation, we also found extensive differences between homologous human and mouse cell types, including marked alterations in proportions, laminar distributions, gene expression and morphology. These species-specific features emphasize the importance of directly studying human brain.