银屑病单细胞数据 GSE151177 分析复现

文章标题为 ***Single-cell transcriptomics applied to emigrating cells from psoriasis elucidate pathogenic versus regulatory immune cell subsets（单细胞转录组学应用于银屑病迁移细胞阐明致病性与调节性免疫细胞亚群）***，发表在 *Journal of Allergy And Clinical Immunology* 上，IF14.2（2024）

数据为 [GSE151177](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE151177)，文章提供了[分析代码](https://github.com/jaehwan79/KruegerLab_scRNAseq_emigrating_cells_human_skin)，但由于是 2021 年的文章，使用的是 Seurat V3，且代码冗余很多，因此我们使用 Seurat V5 进行复现。步骤和方法尽量保持和原文一致。

目录组织和数据下载

依据我的习惯，在项目目录 `scPsoriasis` 下新建如下目录，并将数据下载到 `data/GSE151177/`

目录结构

解压后发现是 69 个 `.gz` 文件（23 个样本 * 3），文件名格式形如 `GSM4567877_Control01_barcodes.tsv.gz`

此外，在其 Github 仓库中还提供了一个表型文件 `phenotype_data.txt`

数据预处理

首先加载包，定义数据、结果和脚本路径，并将路径设置到结果

文件整理

该数据集是标准的 10X 测序输出文件，但我们仍需要先稍做处理

自定义 `organize_files` 函数，将文件按前缀分文件夹并去除多余的前缀：

该函数会自动识别 `.gz` 文件并根据样本前缀创建相应的文件夹

执行完成后的 `data` 目录应如下：

每一个文件夹下都是三个文件：

样本分组

接下来，我们读取 `phenotype_data.txt` 文件，以获取样本的表型信息，帮助我们区分银屑病和对照样本：

根据表型信息，我们将样本分为三组：

银屑病（Psoriasis）\_试剂盒版本 2.0

银屑病（Psoriasis）\_试剂盒版本 3.0

对照（Control）\_试剂盒版本 2.0

匹配文件夹

接下来，我们匹配样本目录，以便后续读取数据

自定义 `get_phenotype_from_dir` 和 `match_sample_dirs` 函数以根据 phenotype 的样本 ID 匹配对应的文件夹

匹配对应的文件夹

命名文件夹向量，使合并后数据的 `orig.ident` 易识别

读取数据并创建 Seurat 对象

使用 `Read10X` 函数读取数据，并创建 Seurat 对象：

数据标准化

在合并样本之前，需要对每个 Seurat 对象进行标准化：

数据整合

为了消除批次效应，使用 `FindIntegrationAnchors()` 函数识别不同数据集中的细胞锚点，然后通过 `IntegrateData()` 函数整合数据：

保存整合数据

最后，我们将整合后的数据保存为 Rda 文件，便于后续分析：