Data and code from: Comparative transcriptomics reveals elevated TCL1A expression in human circulating immune cells across chronic pain conditions
Data files
May 26, 2026 version files 47.99 MB
-
BayesPrism_Analysis_Final.ipynb
12.82 MB
-
combined_dataset_raw_counts.csv
23.45 MB
-
combined_metadata_dryad.csv
25.87 KB
-
Meta-Analysis_Across-Sex.ipynb
1.69 MB
-
Meta-Analysis_Men.ipynb
177.27 KB
-
Meta-Analysis_Women_V2.ipynb
1.53 MB
-
Prepare_BayesPrism_Benchmark_Final.ipynb
2.69 MB
-
Prepare_BayesPrism_cell_reference_Final.ipynb
608.96 KB
-
Prepare_Bulk_Raw_Counts.ipynb
2.89 MB
-
Prepare_neutrophil_single_cell_reference.ipynb
486.63 KB
-
Prepare_PBMC_single_cell_reference_V2.ipynb
64.30 KB
-
process_bayes_prism_pseudobulk_control.R
1.45 KB
-
process_bayes_prism_V2.R
1.66 KB
-
README.md
5.22 KB
-
RNA_Seq_Meta_Analysis.R
2.81 KB
-
run_bayesprism_pseudobulk_control.R
4.77 KB
-
run_bayesprism_V2.R
5.63 KB
-
run_DESEQ2_meta-analysis_Combined.R
4.14 KB
-
run_DESEQ2_meta-analysis_Female.R
5 KB
-
run_DESEQ2_meta-analysis_Male.R
4.02 KB
-
Run_Gene_Ontology.R
1.81 KB
-
UK_Biobank_TCL1A_V12.ipynb
1.51 MB
Abstract
Background
Chronic pain affects millions globally. The complex milieu of underlying risk factors and symptom heterogeneity, coupled with the absence of biomarkers, complicates management. Many studies examine pain through the lens of the human circulating immune system, as immune dysregulation is one potential contributor to pain. Typically, studies analyze both sexes together; however, there are known sex differences in the immune response to disease, with women at higher risk of autoimmunity. Similarly, chronic pain disproportionately affects women. Analyzing both men and women together and adjusting for sex may inadvertently remove sex-specific immune signatures associated with chronic pain, further hindering our understanding of its etiology.
Methods
To circumvent this, using transcriptomic meta-analysis, we reprocessed several bulk RNA-sequencing studies of human circulating immune cells across multiple chronic pain conditions to simultaneously examine both common and previously overlooked sex-specific transcriptomic signatures.
Results
Our meta-analysis encompassed bulk RNA-sequencing samples derived from human circulating immune cells from 142 chronic pain cases and 154 controls across six chronic pain conditions. When examining both sexes together, we identified 19 genes whose expression was altered in chronic pain. Notably, when examining women alone, we identified 34 altered genes, some with roles in autoimmunity, such as TCL1A. Further, we found that TCL1A expression significantly correlated with neuropathic symptom severity. Finally, this altered expression was confirmed at the protein level in women with neuropathic pain in a separate replication analysis.
Conclusion
Through transcriptomic meta-analysis of open-access data, we identified genes conserved across pain conditions and uncovered sex-specific signatures, including TCL1A, as a potential biomarker for neuropathic pain in women.
Dataset DOI: 10.5061/dryad.37pvmcvz9
Description of the data and file structure
This repository contains scripts for the analyses undertaken in our study. Additionally, there are two .csv files-- one contains the raw count matrices from RSEM and the other has metadata. The original FASTQ files can be found in the original studies' repositories. Due to data sharing limitations of the UK Biobank, we are unable to share those files here, but anyone with UK Biobank access can re-create the analysis using our Methods section.
Files and variables
File: Prepare_BayesPrism_Benchmark_Final.ipynb
Description: Preparing files for BayesPrism (performed in R). Further details are provided within the notebook.
File: BayesPrism_Analysis_Final.ipynb
Description: Simple analysis of BP results (extracted from R). Further details are provided within the notebook.
File: combined_metadata_dryad.csv
Description: Metadata for each sample.
Variables
- batch: Study ID
- disease: disease
- Seq: sequencing
- Machine: sequencing machine
- Sex: Genetic Sex
- PBMC: Is this a PBMC dataset(Yes) or whole blood (No)
File: combined_dataset_raw_counts.csv
Description: RSEM raw count matrices. Simple concatenation. No cross-study normalization.
File: Meta-Analysis_Across-Sex.ipynb
Description: Combined sex meta analysis script. Contains code to reproduce figures. Further details are provided within the notebook via detailed commented code.
File: Meta-Analysis_Women_V2.ipynb
Description: Women only meta analysis script. Contains code to reproduce figures. Further details are provided within the notebook via detailed commented code.
File: Prepare_BayesPrism_cell_reference_Final.ipynb
Description: Subsample initial single cell RNA seq atlas and prepare reference data for BayesPrism. Contains code to reproduce figures. Further details are provided within the notebook via detailed commented code.
File: Meta-Analysis_Men.ipynb
Description: Men only meta analysis script. Contains code to reproduce figures. Further details are provided within the notebook via detailed commented code.
File: Prepare_Bulk_Raw_Counts.ipynb
Description: Process Pain cohort raw data and filter additional samples based on quality control metrics. Further details are provided within the notebook via detailed commented code.
File: process_bayes_prism_V2.R
Description: Simple export from R to universal dataframe objects.
File: Prepare_neutrophil_single_cell_reference.ipynb
Description: Identify and isolate Neutrophiles in order to be used in the Bayeprism reference. Further details are provided within the notebook via detailed commented code.
File: RNA_Seq_Meta_Analysis.R
Description: Meta analysis script. Further details are provided within the script via detailed commented code.
File: run_DESEQ2_meta-analysis_Female.R
Description: DESEQ 2 for women only. Further details are provided within the script via detailed commented code.
File: run_bayesprism_V2.R
Description: Code to run BayesPrism algorithm. Further details are provided within the script via detailed commented code.
File: process_bayes_prism_pseudobulk_control.R
Description: Simple export from R similar to process_bayes_prism_V2.R.
File: Run_Gene_Ontology.R
Description: Performs Gene Ontology analysis.
File: Prepare_PBMC_single_cell_reference_V2.ipynb
Description: Subsample large single cell reference for use in BayesPrism. Further details are provided within the notebook via detailed commented code.
File: UK_Biobank_TCL1A_V12.ipynb
Description: UK Biobank analysis. Further details are provided within the notebook via detailed commented code.
File: run_bayesprism_pseudobulk_control.R
Description: Benchmark Bayeprism against test set.
File: run_DESEQ2_meta-analysis_Male.R
Description: DESEQ 2 for men only. Further details are provided within the script via detailed commented code.
File: run_DESEQ2_meta-analysis_Combined.R
Description: DESEQ 2 for both sexes with sex as a covariate. Further details are provided within the script via detailed commented code.
Code/software
Each script/notebook contains detailed commented code which should guide the user through the analysis.
Access information
Other publicly accessible locations of the data:
- See paper for details.
Data was derived from the following sources:
- See paper for details.
Human subjects data
This study was a meta-analysis of previously published transcriptomic studies. All data used was previously de-identified by the original authors and does not contain PII. The original FASTQ files are available without a material transfer agreement and are deposited in their respective depositories outlined in the original publication. Scripts are original to this study and do not contain PII.
