Downloads

To download the original, full 10M cell cytokine perturbation dataset, see the Parse Biosciences website:

10 Million Human PBMCs in a Single Experiment, https://www.parsebiosciences.com/datasets/10-million-human-pbmcs-in-a-single-experiment/; Parse Biosciences.

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 columns
AIFI_L2: Level 2 (intermediate resolution) cell type used for differential expression test
fg: Foreground cytokine in the contrast
bg: Background group (PBS) in the contrast
gene: Gene symbol for the contrast test
censored: 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 test
stat: 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 NameDescriptionDownload 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 NameDescriptionDownload 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 NameDescriptionDownload 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 NameDescriptionDownload 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)