A structural equation model for imaging genetics using spatial transcriptomics.

Imaging genetics deals with relationships between genetic variation and imaging variables, often in a disease context. The complex relationships between brain volumes and genetic variants have been explored with both dimension reduction methods and model-based approaches. However, these models usually do not make use of the extensive knowledge of the spatio-anatomical patterns of gene activity. We present a method for integrating genetic markers (single nucleotide polymorphisms) and imaging features, which is based on a causal model and, at the same time, uses the power of dimension reduction. We use structural equation models to find latent variables that explain brain volume changes in a disease context, and which are in turn affected by genetic variants. We make use of publicly available spatial transcriptome data from the Allen Human Brain Atlas to specify the model structure, which reduces noise and improves interpretability. The model is tested in a simulation setting and applied on a case study of the Alzheimer’s Disease Neuroimaging Initiative.FSW – Publicaties zonder aanstelling Universiteit Leide

Stochastic neighbor embedding as a tool for visualizing the encoding capability of magnetic resonance fingerprinting dictionaries

Objective To visualize the encoding capability of magnetic resonance fingerprinting (MRF) dictionaries. Materials and methods High-dimensional MRF dictionaries were simulated and embedded into a lower-dimensional space using t-distributed stochastic neighbor embedding (t-SNE). The embeddings were visualized via colors as a surrogate for location in low-dimensional space. First, we illustrate this technique on three different MRF sequences. We then compare the resulting embeddings and the color-coded dictionary maps to these obtained with a singular value decomposition (SVD) dimensionality reduction technique. We validate the t-SNE approach with measures based on existing quantitative measures of encoding capability using the Euclidean distance. Finally, we use t-SNE to visualize MRF sequences resulting from an MRF sequence optimization algorithm. Results t-SNE was able to show clear differences between the color-coded dictionary maps of three MRF sequences. SVD showed smaller differences between different sequences. These findings were confirmed by quantitative measures of encoding. t-SNE was also able to visualize differences in encoding capability between subsequent iterations of an MRF sequence optimization algorithm. Discussion This visualization approach enables comparison of the encoding capability of different MRF sequences. This technique can be used as a confirmation tool in MRF sequence optimization.Radiolog

Hierarchical prediction of registration misalignment using a convolutional LSTM: application to chest CT scans

In this paper we propose a supervised method to predict registration misalignment using convolutional neural networks (CNNs). This task is casted to a classification problem with multiple classes of misalignment: "correct" 0-3 mm, "poor" 3-6 mm and "wrong" over 6 mm. Rather than a direct prediction, we propose a hierarchical approach, where the prediction is gradually refined from coarse to fine. Our solution is based on a convolutional Long Short-Term Memory (LSTM), using hierarchical misalignment predictions on three resolutions of the image pair, leveraging the intrinsic strengths of an LSTM for this problem. The convolutional LSTM is trained on a set of artificially generated image pairs obtained from artificial displacement vector fields (DVFs). Results on chest CT scans show that incorporating multi-resolution information, and the hierarchical use via an LSTM for this, leads to overall better F1 scores, with fewer misclassifications in a well-tuned registration setup. The final system yields an accuracy of 87.1%, and an average F1 score of 66.4% aggregated in two independent chest CT scan studies.Radiolog

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.Imaging- and therapeutic targets in neoplastic and musculoskeletal inflammatory diseas

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.Imaging- and therapeutic targets in neoplastic and musculoskeletal inflammatory diseas

BrainScope: interactive visual exploration of the spatial and temporal human brain transcriptome

Spatial and temporal brain transcriptomics has recently emerged as an invaluable data source for molecular neuroscience. The complexity of such data poses considerable challenges for analysis and visualization. We present BrainScope: a web portal for fast, interactive visual exploration of the Allen Atlases of the adult and developing human brain transcriptome. Through a novel methodology to explore high-dimensional data (dual t-SNE), BrainScope enables the linked, all-in-one visualization of genes and samples across the whole brain and genome, and across developmental stages. We show that densities in t-SNE scatter plots of the spatial samples coincide with anatomical regions, and that densities in t-SNE scatter plots of the genes represent gene co-expression modules that are significantly enriched for biological functions. We also show that the topography of the gene t-SNE maps reflect brain region-specific gene functions, enabling hypothesis and data driven research. We demonstrate the discovery potential of BrainScope through three examples: (i) analysis of cell type specific gene sets, (ii) analysis of a set of stable gene co-expression modules across the adult human donors and (iii) analysis of the evolution of co-expression of oligodendrocyte specific genes over developmental stages. BrainScope is publicly accessible at www.brainscope.nl.FSW – Publicaties zonder aanstelling Universiteit Leide

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.Nephrolog

A novel software tool for semi-automatic quantification of thoracic aorta dilatation on baseline and follow-up computed tomography angiography

A dedicated software package that could semi-automatically assess differences in aortic maximal cross-sectional diameters from consecutive CT scans would most likely reduce the post-processing time and effort by the physicians. The aim of this study was to present and assess the quality of a new tool for the semi-automatic quantification of thoracic aorta dilation dimensions. Twenty-nine patients with two CTA scans of the thoracic aorta for which the official clinical report indicated an increase in aortic diameters were included in the study. Aortic maximal cross-sectional diameters of baseline and follow-up studies generated semi-automatically by the software were compared with corresponding manual measurements. The semi-automatic measurements were performed at seven landmarks defined on the baseline scan by two operators. Bias, Bland–Altman plots and intraclass correlation coefficients were calculated between the two methods and, for the semi-automatic software, also between two observers. The average time difference between the two scans of a single patient was 1188 ± 622 days. For the semi-automatic software, in 2 out of 29 patients, manual interaction was necessary; in the remaining 27 patients (93.1%), semi-automatic results were generated, demonstrating excellent intraclass correlation coefficients (all values ≥ 0.91) and small differences, especially for the proximal aortic arch (baseline: 0.19 ± 1.30 mm; follow-up: 0.44 ± 2.21 mm), the mid descending aorta (0.37 ± 1.64 mm; 0.37 ± 2.06 mm), and the diaphragm (0.30 ± 1.14 mm; 0.37 ± 1.80 mm). The inter-observer variability was low with all errors in diameters ≤ 1 mm, and intraclass correlation coefficients all ≥ 0.95. The semi-automatic tool decreased the processing time by 40% (13 vs. 22 min). In this work, a semi-automatic software package that allows the assessment of thoracic aorta diameters from baseline and follow-up CTs (and their differences), was presented, and demonstrated high accuracy and low inter-observer variability

Characterization and Evaluation of the Artemis Camera for Fluorescence-Guided Cancer Surgery

Purpose: Near-infrared (NIR) fluorescence imaging can provide the surgeon with real-time visualization of, e.g., tumor margins and lymph nodes. We describe and evaluate the Artemis, a novel, handheld NIR fluorescence camera.Procedures: We evaluated minimal detectable cell numbers (FaDu-luc2, 7D12-IRDye 800CW), preclinical intraoperative detection of sentinel lymph nodes (SLN) using indocyanine green (ICG), and of orthotopic tongue tumors using 7D12-800CW. Results were compared with the Pearl imager. Clinically, three patients with liver metastases were imaged using ICG.Results: Minimum detectable cell counts for Artemis and Pearl were 2 × 105 and 4 × 104 cells, respectively. In vivo, seven SLNs were detected in four mice with both cameras. Orthotopic OSC-19-luc2-cGFP tongue tumors were clearly identifiable, and a minimum FaDu-luc2 tumor size of 1 mm3 could be identified. Six human malignant lesions were identified during three liver surgery procedures.Conclusions: Based on this study, the Artemis system has demonstrated its utility in fluorescence-guided cancer surgery

Genome-wide coexpression of steroid receptors in the mouse brain: Identifying signaling pathways and functionally coordinated regions

Steroid receptors are pleiotropic transcription factors that coordinate adaptation to different physiological states. An important target organ is the brain, but even though their effects are well studied in specific regions, brain-wide steroid receptor targets and mediators remain largely unknown due to the complexity of the brain. Here, we tested the idea that novel aspects of steroid action can be identified through spatial correlation of steroid receptors with genome-wide mRNA expression across different regions in the mouse brain. First, we observed significant coexpression of six nuclear receptors (NRs) [androgen receptor (Ar), estrogen receptor alpha (Esr1), estrogen receptor beta (Esr2), glucocorticoid receptor (Gr), mineralocorticoid receptor (Mr), and progesterone recep