Spatial Mapping of Genetic Liability to Psychiatric Disorders in the Adult Human Hippocampus.
The 1 match
- [1] § Methods and Materials › ST Data ↔ index.Rmd, lines 19–88 · score 0.60 · Genomics Visium Spatial, Visium spots, Gene Expression, histologically, anterior, platform
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
R Markdown · 114 lines · 9.6 KB · GPL-3.0 · 1 match
- ---
- title: "VistoSeg: Visium Histology Image Segmentation and Processing Pipeline"
- author:
- - name: Madhavi Tippani
- affiliation:
- - &libd Lieber Institute for Brain Development, Johns Hopkins Medical Campus
- email: [email hidden]
- site: bookdown::bookdown_site
- apple-touch-sin: "icon_192.png"
- apple-touch-icon-size: 192
- favicon: "icon_32.png"
- github-repo: "LieberInstitute/VistoSeg"
- documentclass: book
- output:
- bookdown::gitbook: default
- #bookdown::pdf_book: default
- ---
- # Overview {-}
- Histological imaging is a critical first step of the **spatial transcriptomics** workflow, a barcoding-based transcriptome-wide technology released by [10x Genomics](https://www.10xgenomics.com).
- <center>
- 
- </center>
- <center>
- 
- </center>
- ## Why *spatial transcriptomics* or *Visium imaging*? {-}
- Methods like single-nucleus and single cell RNA-sequencing (RNA-seq) can profile single cells transcriptome-wide enabling researchers to identify cell type compositions; however, these methods necessarily destroy information about spatial positioning. On the other hand, multiplexing and in situ sequencing methods can provide spatial information, but have significant limitations on the number of genes that can be processed as well as issues with microscopy and related computational challenges. **Spatial transcriptomics**, including the 10X Genomics Visium platform, provides solutions to these limitations by allowing researchers to quantify gene expression with high spatial resolution. A critical component for the 10X Genomics platform to provide spatial context is the **Visium imaging** of the [Visium gene expression slide](https://www.10xgenomics.com/products/spatial-gene-expression). On this slide, the experimental tissue sections that are being analyzed are mounted onto the **capture areas (A1,B1,C1,D1)** located on the slide. The whole slide is then imaged, producing a large output image file that contains all of the capture areas. This image of the whole slide then has to be subsequently split into individual capture area images (necessarily JPEG or tif), which are then processed accordingly for the downstream gene expression analyses.
- This website describes the steps required to split, visualize and process the Visium images from spatial transcriptomics projects generated by the [10x Genomics Visium](https://www.10xgenomics.com/spatial-transcriptomics) commercial platform.
- \centering
- 
- The above figure describes the **`VistoSeg` pipeline**, **(A)** The data presented here is from tissue sections obtained from three levels of the human dorsolateral prefrontal cortex (DLPFC; posterior, middle and anterior). Each tissue section spans the six cortical layers plus the white matter. **(B)** Shows the original *'Visium gene expression slide'* with 4 capture areas, and the slide scanner used to image the slide. **(C)** Shows the large tif file produced by the slide scanner, which is then split into the respective capture areas using the function [splitSlide](https://github.com/LieberInstitute/VistoSeg/blob/main/code/split.m) described in [Step 1](http://research.libd.org/VistoSeg/step-1-split-visium-histology-whole-slide-image-into-individual-capture-area-images.html). **(D)** Shows the individual tif images of capture areas produced by [splitSlide](https://github.com/LieberInstitute/VistoSeg/blob/main/code/split.m), and the corresponding nuclei segmentations produced by the functions [VNS](https://github.com/LieberInstitute/VistoSeg/blob/main/code/VNS.m) (Visium Nuclei Segmentation) and [refineVNS](https://github.com/LieberInstitute/VistoSeg/blob/main/code/refineVNS.m) explained in [Step 2](http://research.libd.org/VistoSeg/step-2-nuclei-segmentation-of-individual-capture-areas-images.html). **(E)** The tif images from *'(D)'* serve as input to the [Spaceranger](https://support.10xgenomics.com/spatial-gene-expression/software/pipelines/latest/using/count) module (explained in [Step 3](https://research.libd.org/VistoSeg/step-3-spaceranger.html)), which generates [tissue_positions_list.csv](https://github.com/LieberInstitute/VistoSeg/blob/main/pipeline_outputs/spaceranger/tissue_positions_list.csv) file and [scalefactors_json.json](https://github.com/LieberInstitute/VistoSeg/blob/main/pipeline_outputs/spaceranger/scalefactors_json.json) file that contain *'Visium spot metrics'*. **(F)** The function [countNuclei](https://github.com/LieberInstitute/VistoSeg/blob/main/code/countNuclei.m) explained in [Step 4](https://research.libd.org/VistoSeg/step-4-gui-to-count-nuclei-in-a-visium-spot.html), gives the nuclei count per Visium spot info that is stored in [tissue_spot_counts.csv](https://github.com/LieberInstitute/VistoSeg/blob/main/pipeline_outputs/spotspotcheck/tissue_spot_counts.csv) file. **(G)** Finally, the pipeline provides a GUI called [spotspotcheck](https://research.libd.org/VistoSeg/step-4-gui-to-count-nuclei-in-a-visium-spot.html#spotspotcheck) that allows the user to perform visual inspection of the nuclei segmentations by allowing the user to toggle between the Visium and binary images, and also provides zoom in/out options to clearly identify cell bodies within a Visium spot.
- ## Cite `VistoSeg` {-}
- We hope that [`VistoSeg`](http://research.libd.org/VistoSeg/) will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!
- ```
- @article {Tippani2022,
- author = {Tippani, Madhavi and Divecha, Heena R. and Catallini II,
- Joseph L. and Kwon, Sang Ho and Weber, Lukas M. and Spangler, Abby and Jaffe,
- Andrew E. and Hicks, Stephanie C. and Martinowich, Keri and
- Collado-Torres, Leonardo and Page, Stephanie C. and Maynard,
- Kristen R.},
- title = {VistoSeg: processing utilities for high-resolution Visium/Visium-IF
- images for spatial transcriptomics data},
- year = {2022},
- doi = {[https://doi.org/10.1101/2021.08.04.452489](https://doi.org/10.1101/2021.08.04.452489)},
- publisher = {TODO},
- URL = {TODO},
- journal = {bioRxiv}
- }
- ```
- Project lead: [Madhavi Tippani](https://twitter.com/MadhaviTippani), Staff Scientist in the **Imaging Development Group** at the [Lieber Institute for Brain Development](https://www.libd.org/).
- <center>
- <img src="http://lcolladotor.github.io/img/LIBD_logo.jpg" width="250px">
- </center>
- Initial versions of `spotspotcheck` and `countNuclei` were developed by [Joseph L. Catallini II](https://www.linkedin.com/in/joseph-catallini-ii/).
- ## Image Acquisition {-}
- The 10X [Visium Spatial Gene Expression Imaging Guidelines](https://support.10xgenomics.com/spatial-gene-expression/imaging/doc/technical-note-visium-spatial-gene-expression-imaging-guidelines) are followed for acquiring the images. Images are acquired at 40x magnification using a Leica CS2 slide scanner and saved as *'.SVS files.'* These *'.SVS files'* are then exported as *'TIF files'* for downstream analysis. The entire Visium slide (4 capture areas with fiducial frames), is scanned in a single file (~20GB).
- ## Software Requirements {-}
- The pipeline was developed under the following software configuration.
- <center>
- 
- </center>
- `VistoSeg` has been tested on Linux, Windows and MacOS.
- 1. `MATLAB`
- MATLAB version R2019a 64-bit or later is required to run the `VistoSeg` pipeline with the [Image Processing Toolbox](https://www.mathworks.com/products/image.html) preloaded.
- 2. `Memory`
- Visium whole slide images are high resolution, and the typical size of these multiplane tif images produced in-house is ~25GB. The system RAM (we use ~75GB) should be thrice the size of the multiplane tif image to load it into MATLAB and split them into individual capture areas. The rest of the processing, on individual capture tifs can be performed on a system with as little as 16GB of RAM.
- 3. `Installation`
- The pipeline is available at https://github.com/LieberInstitute/VistoSeg, which can be download to your system from the Github website directly or the main repository can be cloned to your system using the following command on terminal/command prompt.
- ```{bash, eval = FALSE}
- git clone https://github.com/LieberInstitute/VistoSeg.git
- ```
- All the code exists in the [code](https://github.com/LieberInstitute/VistoSeg/tree/main/code) directory inside the main [`VistoSeg`](https://github.com/LieberInstitute/VistoSeg) directory.
- The user's working directory on MATLAB should be the path to the [code](https://github.com/LieberInstitute/VistoSeg/tree/main/code) directory in the downloaded repository, to run any functions this pipeline provides. Once the repository is downloaded, the user can run either of the following code to change their working directory on MATLAB to the [code](https://github.com/LieberInstitute/VistoSeg/tree/main/code) directory.
- ```{bash, eval = FALSE}
- cd /path_to_the_downloaded_repository/VistoSeg/code/
- ```
- ```MATLAB
- addpath(genpath('/path_to_the_downloaded_repository/VistoSeg/code/'))
- ```
- ## Getting Help {-}
- If you are using `VistoSeg` and are running into unexpected problems, please [report them](https://github.com/LieberInstitute/VistoSeg/issues) publicly such that other users might benefit from the answers. Thank you!
- ## Data Availability {-}
- The raw Visium `.tif` file is available at through [AWS](https://visium-libd.s3.amazonaws.com/Lieber_Institute_OTS-20-7690_rush_anterior.tif).
- **Other potential datasets**
- 1. [LIBD pilot DLPFC](https://github.com/LieberInstitute/spatialLIBD#raw-data)
- 2. [10x Genomics spatial data sets](https://support.10xgenomics.com/spatial-gene-expression/datasets)
- 3. All input files and output files for the samples images used in every step in the pipeline are provided at the end of the section.
index.Rmd at commit cd8ca31, under GPL-3.0 · at the source
Overview
- Centre for Neuropsychiatric Genetics & Genomics, Division of Psychological Medicine & Clinical Neurosciences, Cardiff University, Cardiff, United Kingdom
- Betul-Ziya Eren Genome and Stem Cell Center, Erciyes Üniversitesi, Kayseri, Türkiye
- UK Dementia Research Institute at Cardiff, Cardiff University, Cardiff, United Kingdom
Abstract
Background: Common genetic variants associated with psychiatric disorders are enriched in genes with high expression specificity for hippocampal neurons. However, to date, these studies have been based on measures of gene expression from nuclei in liquid suspensions, where information on the precise location of cells within the assayed tissue is lost.
Methods: We applied genetically informed spatial mapping of cells for complex traits (gsMAP) to test enrichment of common-variant genetic liability to schizophrenia, bipolar disorder, and major depressive disorder (MDD) in 13 hippocampal subregions according to expression specificity of associated genes using spatial transcriptomic (ST) data from 10 neurotypical adult donors. Then we used cellular deconvolution data from ST spots to perform further gsMAP analyses on cell populations within enriched subregions.
Results: Compared with other hippocampal subregions, common-variant liability to schizophrenia, bipolar disorder, and, to a lesser extent, MDD, was significantly enriched in genes with higher expression specificity for the subiculum and neuron-rich areas of the cornu ammonis subfields and dentate gyrus. Within implicated regions, enrichments for schizophrenia and bipolar disorder were pronounced in ST spots containing glutamatergic neurons. Genes with higher expression specificity for the granule cell layer and subgranular zone of the dentate gyrus were more significantly enriched for association with bipolar disorder than for schizophrenia.
Conclusions: Our findings support a role for glutamatergic neurons in all major subregions of the hippocampus in mediating common-variant liability to schizophrenia and bipolar disorder and provide evidence for the greater relative importance of the dentate gyrus in bipolar disorder.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repository
Its files are read in the Code ↔ Paper reader above, with 1 match between paragraphs and lines of code.
LieberInstitute/VistoSeg
cd8ca310eaf6ca04909e6ca3a4e36d721c5b1d1d, 3 August 2023Availability: 1 check, the latest on 30 September 2026: the link answers
- 30 September 2026: the link answers
21 files
- Step1.Rmd, R, 106 lines
- Step2.Rmd, R, 96 lines
- Step3.Rmd, R, 21 lines
- Step4.Rmd, R, 111 lines
- Step5.Rmd, R, 129 lines
- XX-Rinfo.Rmd, R, 17 lines
- code/
InFormStitch.m , MATLAB, 52 lines - code/
VNS.m , MATLAB, 39 lines - code/
clickPlot.m , MATLAB, 90 lines - code/
countNuclei.m , MATLAB, 34 lines - code/
countSpots.m , MATLAB, 56 lines - code/
countSpots_centroid.m , MATLAB, 66 lines - code/
extractMD.m , MATLAB, 20 lines - code/
refineVNS.m , MATLAB, 34 lines - code/
splitSlide.m , MATLAB, 62 lines - code/
splitSlide_IF.m , MATLAB, 37 lines - code/
spotspotcheck.m , MATLAB, 83 lines - code/
tif2mat.m , MATLAB, 12 lines - index.Rmd, R, 114 lines, 1 match
- LICENSE, License, 674 lines
- README.md, Text, 13 lines
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 19 scripts, each with its path and the digest of its content;
- 1 match between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
Datasets cited
- geo:GSE264624, at NCBI GEO; found in the text, “ST Data”
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 30 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 7 authors, 7 keywords, 2 funders, 39 references.
Cite
This paper
Baran, Y., Cameron, D., Pryce-Roberts, A., Richards, A. L., Webber, C., O’Donovan, M. C., & Bray, N. J. (2026). Spatial Mapping of Genetic Liability to Psychiatric Disorders in the Adult Human Hippocampus. Biological psychiatry global open science, 6(3), 100719. https://
BibTeX
@article{baran2026spatia
author = {Baran, Yusuf and Cameron, Darren and Pryce-Roberts, Adele and Richards, Alexander L and Webber, Caleb and O’Donovan, Michael C and Bray, Nicholas J},
title = {{Spatial Mapping of Genetic Liability to Psychiatric Disorders in the Adult Human Hippocampus}},
journal = {Biological psychiatry global open science},
year = {2026},
month = mar,
volume = {6},
number = {3},
pages = {100719},
publisher = {Elsevier},
issn = {2667-1743},
doi = {10.1016/
url = {https://
pmid = {42016326},
pmcid = {PMC13094429}
}
RIS
TY - JOUR
AU - Baran, Yusuf
AU - Cameron, Darren
AU - Pryce-Roberts, Adele
AU - Richards, Alexander L
AU - Webber, Caleb
AU - O’Donovan, Michael C
AU - Bray, Nicholas J
TI - Spatial Mapping of Genetic Liability to Psychiatric Disorders in the Adult Human Hippocampus
T2 - Biological psychiatry global open science
J2 - Biol Psychiatry Glob Open Sci
PY - 2026
DA - 2026/
VL - 6
IS - 3
SP - 100719
SN - 2667-1743
PB - Elsevier
DO - 10.1016/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1016/
"type": "article-journal",
"title": "Spatial Mapping of Genetic Liability to Psychiatric Disorders in the Adult Human Hippocampus",
"container-title": "Biological psychiatry global open science",
"author": [
{
"family": "Baran",
"given": "Yusuf"
},
{
"family": "Cameron",
"given": "Darren"
},
{
"family": "Pryce-Roberts",
"given": "Adele"
},
{
"family": "Richards",
"given": "Alexander L"
},
{
"family": "Webber",
"given": "Caleb"
},
{
"family": "O’Donovan",
"given": "Michael C"
},
{
"family": "Bray",
"given": "Nicholas J"
}
],
"container-title-short":
"volume": "6",
"issue": "3",
"page": "100719",
"DOI": "10.1016/
"PMID": "42016326",
"PMCID": "PMC13094429",
"ISSN": "2667-1743",
"publisher": "Elsevier",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
3,
5
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1016/j.celrep.2026.117500 [code]
- Spatio-molecular gene expression reflects dorsal anterior cingulate cortex structure and function in the human brain.Journal: Cell reportsIn common: genetics / omics, cellular / molecular, 7 references
- [2] doi:10.1038/s41398-026-04200-5 [code]
- Postmortem brain single-nucleus and bulk gene expression analyses identify shared and distinct abnormalities in bipolar disorder and major depressive disorder.Journal: Translational psychiatryIn common: bipolar, depression, genetics / omics, 1 other category, 4 references
- [3] doi:10.1038/s41380-026-03497-4 [code]
- Transcriptome-informed brain cartography of polygenic risk and association with brain structure in major psychiatric disorders.Journal: Molecular psychiatryIn common: bipolar, schizophrenia / psychosis, depression, 2 other categories, 3 references
- [4] doi:10.1162/imag.a.1152 [code]
- Brain patterns linked to neuropsychiatric genetic risk mirror those seen in disease.Journal: Imaging neuroscience (Cambridge, Mass.)In common: schizophrenia / psychosis, genetics / omics, 5 references
- [5] doi:10.1371/journal.pcbi.1014422 [code]
- Deciphering cell type-specific causal genetic effects on brain imaging-derived phenotypes and disorders with single-cell Mendelian randomization.Journal: PLoS computational biologyIn common: genetics / omics, cellular / molecular, 5 references
- [6] doi:10.1073/pnas.2609814123 [code]
- Genetic architectures of brain-related traits are shaped by strong selective constraints.Journal: Proceedings of the National Academy of Sciences of the United States of AmericaIn common: genetics / omics, cellular / molecular, 5 references
- [7] doi:10.1038/s41588-026-02646-3 [code]
- Co-expression-based models improve eQTL predictions for transcriptome-wide association studies and highlight new schizophrenia-associated
genes. Journal: Nature geneticsIn common: schizophrenia / psychosis, genetics / omics, cellular / molecular, 4 references - [8] doi:10.64898/2026.03.06.709690
- Genetic insights on the mechanisms of human cortical foldingJournal: bioRxiv (preprint)In common: genetics / omics, cellular / molecular, 5 references
- [9] doi:10.1038/s41380-026-03571-x [code]
- Convergent coexpression reveals shared biological mechanisms underlying common and rare variant risk in six neuropsychiatric disorders.Journal: Molecular psychiatryIn common: genetics / omics, cellular / molecular, 5 references
- [10] doi:10.1038/s41562-026-02486-5 [code]
- Genome-wide association studies of infant and toddler temperament in European and multi-ancestry populations.Journal: Nature human behaviourIn common: genetics / omics, 5 references
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 19 scripts, and 1 match between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:8820f222b18dcf7a…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
