Source code for the PBMC Ageing Atlas analysis that includes over a million cells from young and old donors and seven studies.
![image](https://private-user-images.githubusercontent.com/13019221/326982984-55eb7ce5-1772-4b24-8bf7-56e3e977d10a.png?jwt=eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJpc3MiOiJnaXRodWIuY29tIiwiYXVkIjoicmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbSIsImtleSI6ImtleTUiLCJleHAiOjE3MjUzMTA4NTMsIm5iZiI6MTcyNTMxMDU1MywicGF0aCI6Ii8xMzAxOTIyMS8zMjY5ODI5ODQtNTVlYjdjZTUtMTc3Mi00YjI0LThiZjctNTZlM2U5NzdkMTBhLnBuZz9YLUFtei1BbGdvcml0aG09QVdTNC1ITUFDLVNIQTI1NiZYLUFtei1DcmVkZW50aWFsPUFLSUFWQ09EWUxTQTUzUFFLNFpBJTJGMjAyNDA5MDIlMkZ1cy1lYXN0LTElMkZzMyUyRmF3czRfcmVxdWVzdCZYLUFtei1EYXRlPTIwMjQwOTAyVDIwNTU1M1omWC1BbXotRXhwaXJlcz0zMDAmWC1BbXotU2lnbmF0dXJlPTg0ZmIzNmRmZjQ0OGU1ZWQ3ODNmMzgxYmViNmU1N2NiMjVjMTc0MmQ3MmNhM2VjN2E1YjBmOWNlZjZmMmRiNjgmWC1BbXotU2lnbmVkSGVhZGVycz1ob3N0JmFjdG9yX2lkPTAma2V5X2lkPTAmcmVwb19pZD0wIn0.f9FKdLWWvX6feqY6RK3wXrMxx_LvqpjH7uOE3pYwyCY)
The datasets can be found under the following accessions on NCBI, GSA, and Synapse: GSE157007, GSE213516, GSE214546, HRA000203, HRA000624, HRA003766, syn22255433.
Please note that scVI
scripts run much faster on GPU.
GSE157007
, GSE213516
, GSE214546
, HRA000203
, HRA000624
, HRA003766
, syn22255433
contain scripts to combine the sample data for each study and perform the quality control. The scripts in each folder should be run in the following order:
create_adatas.py
to import the data into AnnData format.doublets.R
to perform the doublet calls for each sample.combine.py
to create a single AnnData file with all the samples and doublet calls for each dataset.qc.py
to peform the quality control for each dataset.
atlas
contains code for integrating the seven datasets and performing the downstream analyses. The scripts in each folder should be run in the following order:
combine_datasets.py
to create a single AnnData object for seven datasets.prepare_combined.py
to peform the clean up, such as removing V(D)J genes.integrate.py
to run scVI integration (preferably on a GPU).viz_integrated.py
to perform additional QC (removing doublets and RBC contamination) and PBMC annotation.
T_integrate.py
, B_integrate.py
, MALAT1_integrate.py
perform the T, B, and MALAT1+ cell re-analysis, respectively.
celltypist_run.py
and celltypist_viz.py
perform CellTypist classification for T cells.
scanpy_DE.py
performs DE test between young and old individuals.
cell_type_props_plots.R
plots of cell type proportions in the datasets.
harmony_integrate.py
, harmony_cluster.py
- alternative integration with Harmony.
abundance_heatmap.R
- cell type proportions heatmap in each sample and dataset.
T_markers
- visualization of the T cell reference for validating marker genes.
utils_py
- shared Python utilities.
This study was supported by the funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No.: 955321, as well by Estonian Research Council grant PRG2011. This publication is based upon work partially supported by the Google Cloud Research Credits program award No.: GCP19980904.