To download the original, full 10M cell cytokine perturbation dataset, see the Parse Biosciences website:
The Parse Biosciences dataset and downloads derived by the Allen Institute are licensed under the Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC 4.0) license. Non-commercial users are free to share and adapt the data by providing appropriate credit. See the Citation & Contributors page for citation information. Contact Parse Biosciences (support@parsebiosciences.com), for commercial use terms related to the Parse Biosciences dataset.
Differentially expressed genes
We utilized DESeq2 to compare pseudobulk scRNA-seq population data between cytokine treatments and the PBS control condition for the 12 subjects in this study. These comparisons were performed with all available genes, as well as with a subset of available genes after removing mitochondrial, ribosomal, and pseudogenes (see note in the descriptors below). The results of these pairwise comparisons are provided below, and descriptions of the columns are available by clicking on the header below:
DEG Results columnsAIFI_L2: Level 2 (intermediate resolution) cell type used for differential expression testfg: Foreground cytokine in the contrastbg: Background group (PBS) in the contrastgene: Gene symbol for the contrast testcensored: Boolean indicating whether or not the gene set was censored to remove mitochondrial, ribosomal, and non-coding genes*.log2fc: Log2(Fold Change) of average expression between groups.padj: Adjusted P-value for the contrast test (using the Benjamini and Hochberg FDR method)pvalue: Original P-value generated by DESeq2 for the contrast teststat: Test statistic generated by DESeq2 for the contrast test
*Censoring was performed by removing genes found in HGNC Gene Groups including the words "Mitochondria" or "Ribosomal" in the Group name field, as well as genes with "RNA" or "pseudogene" in the Locus type field. We also excluded genes beginning with "7SK", "RN7SK", or "MT-". If `censored` is False, all genes were included in analysis.
Cytokine vs. PBS Pseudobulk DESeq2 DEGs
| File Name | Description | Download Link | 
|---|---|---|
| parse_cytokine_deg_censored.csv | Censored DESeq2 DEGs | |
| parse_cytokine_deg_uncensored.csv | Uncensored DESeq2 DEGs | 
Gene set enrichment analysis results
After performing DEG analysis, we ranked DEGs for each cytokine by nominal p-value and direction of change relative to the PBS control, then performed Gene Set Enrichment Analysis (GSEA) using the fgsea package for R against gene sets from MSigDB (Hallmark pathways) and Reactome. These results are used to display Pathway enrichment in our DEG Explorer visualization tool.
Cytokine vs. PBS Gene Set Enrichment Results
| File Name | Description | Download Link | 
|---|---|---|
| parse_cytokine_censored_hallmark_gsea.tsv | Censored DEG Hallmark GSEA | |
| parse_cytokine_censored_reactome_gsea.tsv | Censored DEG Reactome GSEA | |
| parse_cytokine_uncensored_hallmark_gsea.tsv | Uncensored DEG Hallmark GSEA | |
| parse_cytokine_uncensored_reactome_gsea.tsv | Uncensored DEG Reactome GSEA | 
Gene sets for cytokine response enrichment analysis
After identifying differentially expressed genes, we ranked DEGs for each cytokine based on the direction of differential expression relative to PBS (up- or down-regulated), and selected the top 500 genes in each direction for each treatment condition. We provide these gene sets for use in Gene Set Enrichment Analysis tools, such as fgsea, to assist in identifying cytokine signatures in differential expression results.
Cytokine vs. PBS Gene Sets for GSEA
| File Name | Description | Download Link | 
|---|---|---|
| parse_cytokine_gsea_gene_sets.csv | Gene sets for GSEA | 
CellTypist models for Parse PBMCs
To apply cell type labels that matched the Human Immune Health Atlas, we carefully labeled the PBS control cells in the Parse GigaLab dataset to match the cell types available in the atlas. We then used the labeled data to build CellTypist labeling models for AIFI_L1 (low resolution cell classes) and AIFI_L2 (intermediate resolution cell types), and applied these models to the cytokine perturbation datasets.
To enable other researchers to label their own Parse Biosciences PBMC scRNA-seq data, we provide the model files below, suitable for annotation using the celltypist package in Python.
CellTypist models for Parse PBMC data
| File Name | Description | Download Link | 
|---|---|---|
| aifi_parse_pbmc_celltypist_model_AIFI_L1.pkl | Coarse cell class model | |
| aifi_parse_pbmc_celltypist_model_AIFI_L2.pkl | Intermediate resolution cell type model | 
scRNA-seq data with AIFI labels
scRNA-seq data from the Parse GigaLab cytokine perturbation dataset were obtained from Parse Biosciences. We then carefully labeled the PBS control sample using many of the same cell type annotations found in our Human Immune Health Atlas at Level 2 (intermediate resolution).
Below, we provide data for subsets of the Parse GigaLab cytokine perturbation dataset based on family groupings of the 90 cytokines with our additional cell type labels. Cell and sample metadata generated by Parse Biosciences are included as well. Click the header below for descriptions of these metadata:
We group the scRNA-seq data based on families of cytokines/treatment molecules. The treatment groups are shown below along with the gene symbol for the gene that encodes each molecule in parentheses if it differs from the molecule name.
Control:
PBS
gp130 Cytokines:
CT-1 (CTF1), IL-6, IL-11, IL-27, LIF, OSM
Non-gp130 Cytokines:
G-CSF (CSF3), IL-12 (IL12A), IL-23, IL-31, IL-35 (EBI3), Leptin (LEP)
RTK Cytokines:
FTL3L (FTL3LG), IL-34, M-CSF (CSF1), SCF (KITLG)
Non-RTK Cytokines:
IL-1α (IL1A), IL-1β (IL1B), IL-1Ra (IL1RN), IL-16, IL-18, IL-33, IL-36α (IL36A), IL-36Ra (IL36RN)
Rβc Family Cytokines:
GM-CSF (CSF2), IL-3, IL-5
Rɣc Family Cytokines:
IL-2, IL-4, IL-7, IL-9, IL-13, IL-15, IL-21, TSLP
IL-10 Family Cytokines:
IL-10, IL-19, IL-20, IL-22, IL-24, IL-26
IL-17 Family Cytokines:
IL-17A, IL-17B, IL-17C, IL-17D, IL-17E, IL-17F
Interferons:
IFN-α1 (IFNA1), IFN-β (IFNB1), IFN-ε (IFNE), IFN-γ (IFNG), IFN-λ1 (IFNL1), IFN-λ2 (IFNL2), IFN-λ3 (IFNL3), IFN-ω (IFNW1)
TNF Ligands, part 1:
4-1BBL (TNFSF9), APRIL (TNFSF13), BAFF (TNFSF20), CD27L (TNFSF7), CD30L (TNFSF8), CD40L (TNFSF5), FasL (TNFSF6), GITRL (TNFSF18), LIGHT (TNFSF14),
TNF Ligands, part 2:
LT-α1-β2 (TNFSF1/TNFSF3), LT-α2-β1 (TNFSF1/TNFSF3), OX40L (TNFSF4), RANKL (TNFSF11), TL1A (TNFSF15), TNF-α (TNFSF2), TRAIL (TNFSF10), TWEAK (TNFSF12)
Growth Factors:
EGF, EPO, FGF-β (FGF2), GDNF, HGF, IFG-1 (IGF1), Noggin (NOG), PSPN, TPO (THPO), VEGF (VEGFA)
Other:
ADSF (RETN),C3a (C3), C5a (C5), Decorin (DCN), IL-32-β, IL-8 (CXCL8), PRL, TGF-β1 (TGFB1)