Human spleen organoid single-cell trancsriptomics, enriched for γδ T cells
Data files
Apr 02, 2026 version files 622.26 MB
-
_1_p258a_CKDL240012182-1A_225CC7LT4_L2_1_Seurat.rds
339.80 MB
-
_1_p258b_CKDL240012182-1A_225CC7LT4_L2_1_Seurat.rds
282.44 MB
-
merge_rhapsody.R
1.78 KB
-
README.md
3.39 KB
-
SampleMetadata_for_Dryad.xlsx
11.88 KB
Abstract
Vaccine effectiveness against malaria is dramatically reduced in malaria-exposed compared to malaria-naïve populations, potentially due to altered immune responses in secondary lymphoid organs following repeated infection. Newly developed human tonsil and spleen organoids, which replicate key features of B and T cell immunity, provide an exciting opportunity to overcome challenges of other models and to improve our understanding of innate-adaptive interactions in lymphoid tissue. The objectives of this study were to use these organoids to investigate the impact of malaria parasites on 1) cells within lymphoid tissues and 2) responses to a heterologous antigen.
To better characterize the impact of iRBC stimulation on human lymphoid reponses, we cocultured three human spleen organoids with iRBC and uninfected red blood cells for 10 days. We then enriched for Vd2+ cells using fluorescence-activated cell sorting, and used the BD Rhapsody platform to generate multiomic single-cell data, combining transcriptomic, T cell receptor sequence, and surface protein expression data. Data contained in this repository represents the Seurat object generated after alignment and cell calling, but before quality control. We hope these data are useful for researchers investigating the utility of organoids as a means of enabling tissue-specific immunity in humans. Our experiments were done with post-autopsy human tissues, following established protocols and institutional ethics review. We report no conflicts of interest.
Dataset DOI: 10.5061/dryad.3r2280gw3
Description of the data and file structure
BD Rhapsody v2 beads were used to capture transcriptomes, TCR sequences and oligo tagged antibodies bound to surface proteins of γδ-enriched spleen organoids. Two Seurat objects were generated from the alignments of these sequencing libaries for, one object for each lane on the sequencer that was used. The SampleMetadata.xslx file contains metadata on the samples that were used in this experiment.
Files and variables
File: _1_p258b_CKDL240012182-1A_225CC7LT4_L2_1_Seurat.rds
Description: Seurat object made from libraries sequences in lane 2
File: _1_p258a_CKDL240012182-1A_225CC7LT4_L2_1_Seurat.rds
Description: Seurat object made from libraries sequences in lane 1
Both RDS files contain Seurat objects containing transcriptomic, antigen receptor sequence, and surface protein expression data; these objects are minimally processed, meaning quality control, doublet removal, normalization, scaling, clustering etc. still have to be performed. Please refer to Press et al. 2026, PLOS Pathogens for the analysis methods.
File: SampleMetadata_for_Dryad.xlsx
Description: Metadata for each sample used in the experiment. These data mostly are also contained in the Seurat objects metadata. The sheet "Sort" contains experimental medatata in the following columns:
- sample #
- numeric sample ID
- donor #
- de-identified donor ID
- organ
- tissue of origin
- stim condition
- stimulation condition
- viability (%)
- percentage viable cells before staining
- sorted gd
- number of γδ T cells sorted
- sorted non-gd
- number of cells sorted that are not γδ T cells
- total cells added to capture
- sum of γδ T and non-γδ T cells added to capture
- % cells gd
- percentage of cells in capture that were γδ T
- Notes
- experimental notes
- HIMC tag
- numeric hash tag for multiplexing
- HIMC cells/ul
- post-sort cell count before loading
- HIMC viability (%)
- post-sort cell viability
- total cells captured
- total cells captured on BD rhapsody cartridge
- total live cells
- total live cells captured on BD rhapsody cartridge
- cartridge
- cartridge ID
The "Donor info" sheet contains age and sex of each organ donor in the study.
File: merge_rhapsody.R
Description: R script to merge the two Seurat objects and add experimental metadata.
Code/software
64GB of RAM are recommended for loading and manipulating these data. You can load these Seurat objects and merge them with the metadata by running the the code contained in merge_rhapsody.R, to do this we recommend installing R version 4.2, Seurat version 4.3, dplyr version 1.0.10 and readxl version 1.4.
After running the above code, you should be able to access all relevant experimental and metadata of this single-cell RNAseq dataset.
Human subjects data
No personally identifiable data is contained in the data published in this repository.
Written informed consent was obtained from all participants or legal representatives. Ethics approval for malaria-naïve samples was obtained from Stanford University’s Human and Animal Research Subjects Compliance Office (Stanford IRB-60741 and IRB-30387).
Spleens organoids were prepared as per Wager et al. 2021 (PMC7891554). Briefly, cryopreserved splenic cell suspensions were thawed and counted, 6 x 106 cells in 100 µl immune organoid media were added to each transwell (Millipore PICM01250) in 12-well plates. Immune organoid media containing 0.5 µl g/mL recombinant human B cell activation factor (BAFF, Biolegend 559608) was added to the outside of the well. Spleen tissue from 3 donors was used and an aliquot of cells from d0 was frozen for each. Organoids were set up and stimulated with uRBC or 3D7 iRBC (12 x 106 iRBC or uRBC/well) for 10 days with media changes every 3-4 days. Stimulated cells, as well as d0 cells, were then stained with an anti-Vδ2+ antibody (PerCP-conjugated, 1:50) and a cell viability dye (Live/Dead Aqua, 1:400) for 30 minutes, washed, and then sorted on a Sony cell sorter. For each sample, 8,000-40,000 Vδ2+ γδ T cells were sorted and depending on the number collected, 10,000-30,000 Vδ2- cells were sorted in order to reach a total of 40,000-50,000 cells per sample. Cells were sorted into cold PBS containing 2% fetal bovine serum. Fewer cells were collected for corresponding PBMC samples, which as would be expected from a blood sample following autopsy, had a lot of dead cells and debris. The percentage of total cells that were Vδ2+ T cells ranged from 20-80% per sample (80% for all iRBC-stimulated cells, lower for some uRBC-stimulated and d0). Samples were transferred to the Stanford Genomics Core, where they were then barcoded and stained with 13 BD AbSeq Antibody-Oligonucleotides Conjugates corresponding to cell markers of interest (see table below). Single cells (15,000-35,000 live cells per sample with the exception of the PBMC sample, which had 4925) were captured on the BD Rhapsody and mRNA was captured and sequenced. The BD Rhapsody Analysis Pipeline was used to generate sequencing data for transcriptomes, as well as T and B cell receptors.
