用于研究宿主-病原体相互作用的高通量转录组分析

Summary

这里介绍的方案描述了一个完整的管道，用于分析从原始读数到功能分析的RNA测序转录组数据，包括质量控制和预处理步骤到高级统计分析方法。

Abstract

病原体可引起多种传染病。宿主对感染的反应诱导的生物过程决定了疾病的严重程度。为了研究这些过程，研究人员可以使用高通量测序技术（RNA-seq），以测量宿主转录组在感染的不同阶段，临床结果或疾病严重程度的动态变化。这项调查可以更好地了解疾病，并发现潜在的药物靶点和治疗方法。这里介绍的方案描述了一个完整的管道，用于分析从原始读取到功能分析的RNA测序数据。管道分为五个步骤:（1）数据的质量控制;（2）基因的作图和注释;（3）统计分析，鉴定差异表达基因和共表达基因;（4）测定样品扰动的分子程度;（5）功能分析。步骤 1 删除了可能影响下游分析质量的技术工件。在第2步中，根据标准文库协议绘制基因图谱并进行注释。步骤3中的统计分析可识别感染样本中差异表达或共表达的基因，与未感染的基因进行比较。在步骤4中使用分子扰动度方法验证样品变异性和潜在生物异常值的存在。最后，步骤5中的功能分析揭示了与疾病表型相关的途径。所提出的管道旨在通过来自宿主 - 病原体相互作用研究的RNA-seq数据分析来支持研究人员，并推动未来的体外 或 体内 实验，这对于了解感染的分子机制至关重要。

Introduction

登革热、黄热病、基孔肯雅热和寨卡等虫媒病毒与几次地方性疫情广泛相关，并已成为过去几十年中导致人类感染的主要病原体之一^1，2。感染基孔肯雅病毒（CHIKV）的个体经常出现发热、头痛、皮疹、多关节痛和关节炎^3，4，5。病毒可以破坏细胞的基因表达并影响各种宿主信号通路。最近，血液转录组研究利用RNA-seq鉴定与急性CHIKV感染相关的差异表达基因（DEGs）与康复⁶ 或健康对照组⁷进行比较。CHIKV感染的儿童具有上调的基因，这些基因参与先天免疫，例如与病毒RNA，JAK / STAT信号传导和Toll样受体信号通路的细胞传感器相关的基因⁶。急性感染CHIKV的成人还显示出与先天免疫相关的基因的诱导，例如与单核细胞和树突状细胞活化相关的基因，以及与抗病毒反应^{相关的基因7}。富含下调基因的信号通路包括与适应性免疫相关的信号通路，例如T细胞活化以及T细胞和B细胞的分化和富集⁷。

可以使用几种方法来分析宿主和病原体基因的转录组数据。通常，RNA-seq文库的制备从成熟poly-A转录本的富集开始。该步骤去除大部分核糖体RNA（rRNA），在某些情况下去除病毒/细菌RNA。然而，当生物学问题涉及病原体转录本检测并且RNA独立于先前的选择进行测序时，可以通过测序检测许多其他不同的转录本。例如，亚基因组mRNA已被证明是验证疾病严重程度的重要因素⁸。此外，对于某些病毒，如CHIKV和SARS-CoV-2，即使是富含poly-A的文库也会生成病毒读数，可用于下游分析^9，10。当专注于宿主转录组的分析时，研究人员可以研究样品之间的生物扰动，鉴定差异表达的基因和富集的途径，并生成共表达模块^7，11，12。该协议突出了使用不同生物信息学方法对CHIKV感染患者和健康个体的转录组分析（图1A）。来自先前发表的一项研究⁷ 的数据包括20名健康和39名CHIKV急性感染个体，用于产生具有代表性的结果。

Protocol

该协议中使用的样品由圣保罗大学生物医学科学研究所微生物学系和塞尔希培联邦大学伦理委员会批准（协议分别为54937216.5.0000.5467和54835916.2.0000.5546）。

1. Docker 桌面安装

注意:准备 Docker 环境的步骤因操作系统 （OS） 而异。因此，Mac 用户必须遵循列为 1.1 的步骤，Linux 用户必须遵循列为 1.2 的步骤，Windows 用户必须遵循列为 1.3 的步骤。

1. 在 MacOS 上安装。
   1. 访问 Get Docker 网站（材料表），单击 Docker Desktop for Mac ，然后单击" 从 Docker Hub 下载 "链接。
   2. 通过单击" 获取 Docker "按钮下载安装文件。
   3. 执行 Docker.dmg文件以打开安装程序，然后将图标拖到 "应用程序" 文件夹中。本地化并执行" 应用程序" 文件夹中的 Docker.app 以启动程序。
      注:顶部状态栏中的软件特定菜单指示软件正在运行，并且可从终端访问该软件。
2. 在 Linux 操作系统上安装容器程序。
   1. 访问 Get Docker Linux 网站（材料表），并按照使用 Docker Linux 存储库链接上提供的存储库部分进行安装的说明进行操作。
   2. 使用命令行更新所有 Linux 软件包:
      sudo apt-get update
   3. 将所需的软件包安装到 Docker:
      sudo apt-get install apt-transport-https ca-certificates curl gnupg lsb-release
   4. 创建软件归档密钥环文件:
      curl -fsSL https://download.docker.com/linux/ubuntu/gpg |sudo gpg --dearmor -o /usr/share/keyrings/docker-archive-keyring.gpg
   5. 在 source.list 文件中添加 Docker deb 信息:
      echo "deb [arch=amd64 signed-by=/usr/share/keyrings/docker-archive-keyring.gpg] https://download.docker.com/linux/ubuntu $（lsb_release -cs） stable" |sudo tee /etc/apt/sources.list.d/docker.list > /dev/null
   6. 再次更新所有软件包，包括最近添加的软件包:
      sudo apt-get update
   7. 安装桌面版本:
      sudo apt-get install docker-ce docker-ce-cli containerd.io
   8. 选择地理区域和时区以完成安装过程。
3. 在 Windows 操作系统上安装容器程序。
   1. 访问 Get Docker 网站（材料表），然后单击" 入门"。查找 Docker Desktop for Windows 的安装程序。下载文件并将其本地安装在计算机上。
   2. 下载后，启动安装文件（.exe）并保留默认参数。确保标记了"安装 WSL 2 所需的 Windows 组件"和"将快捷方式添加到桌面"这两个选项。
      注意:在某些情况下，当此软件尝试启动服务时，它会显示错误:WSL 安装不完整。要找出此错误，请访问网站 WSL2-Kernel（材料表）。
   3. 下载并安装最新的 WSL2 Linux 内核。
   4. 以管理员身份访问 PowerShell 终端并执行以下命令:
      dism.exe /online /enable-feature /featurename:Microsoft-Windows-Subsystem-Linux /all /norestart
   5. 确保软件 Docker 桌面已成功安装。
4. 从 Docker 中心（材料表）上的 CSBL 存储库下载映像。
   1. 打开 Docker 桌面，并验证工具栏左下角的状态是否为"正在运行"。
   2. 转到 Windows PowerShell 终端命令行。从 Docker 中心上的 CSBL 存储库下载此协议的 Linux 容器映像。执行以下命令下载镜像:
      Docker pull csblusp/transcriptome
      注意:下载映像后，可以在 Docker 桌面中看到该文件。若要创建容器，Windows 用户必须遵循步骤 1.5，而 Linux 用户必须遵循步骤 1.6。
5. 初始化 Windows 操作系统上的服务器容器。
   1. 从工具栏在桌面应用管理器中查看 Docker 映像文件，并访问"映像"页。
      注意:如果管道映像已成功下载，则将有一个 csblusp/transcriptome 映像可用。
   2. 通过单击" 运行 "按钮从 csblusp/transcriptome 图像启动容器。展开 "可选设置" 以配置容器。
   3. 定义容器名称（例如，服务器）。
   4. 将本地计算机中的文件夹与 Docker 中的文件夹相关联。为此，请确定主机路径。在本地计算机中设置一个文件夹，以存储将在结束时下载的已处理数据。设置容器路径。定义 csblusp/transcriptome 容器文件夹并将其链接到本地计算机路径（使用名称"/opt/transferdata"作为容器路径）。
   5. 之后，单击" 运行 "以创建 csblusp/transcriptome 容器。
   6. 要从 csblusp/transcriptome 容器访问 Linux 终端，请单击 CLI 按钮。
   7. 键入 bash 终端以获得更好的体验。为此，请执行以下命令:
      砰
   8. 执行 bash 命令后，确保终端显示 （root@:/#）:
      root@ac12c583b731:/#
6. 初始化 Linux 操作系统的服务器容器。
   1. 执行以下命令以基于映像创建 Docker 容器:
      docker run -d -it --rm --name server -v :/opt/transferdata csblusp/transcriptome
      注意:<主机路径>:定义本地文件夹计算机的路径。
   2. 执行以下命令以访问 Docker 容器的命令终端:
      Docker exec -it server bash
   3. 确保 Linux 终端的可用性，以便使用命令行执行任何程序/脚本。
   4. 执行 bash 命令后，确保终端显示 （root@:/#）:
      root@ac12c583b731:/#
      注意:默认情况下，root 密码为"转录组"。如果需要，可以通过执行以下命令来更改 root 密码:
      路卫
   5. 首先，执行 source 命令以 addpath.sh 以确保所有工具都可用。执行命令:
      源 /选择/添加路径.sh
7. 检查RNA测序文件夹的结构。
   1. 访问转录组管道脚本文件夹，并确保来自RNA测序的所有数据都存储在该文件夹中:/home/transcriptome-pipeline/data。
   2. 确保从分析中获得的所有结果都存储在路径/home/transcriptome-pipeline/results的文件夹中。
   3. 确保基因组和注释参考文件存储在路径/home/transcriptome-pipeline/datasets的文件夹中。这些文件将有助于支持所有分析。
   4. 确保所有脚本都存储在路径/home/transcriptome-pipeline/scripts的文件夹中，并按每个步骤分隔，如下所述。
8. 下载注释和人类基因组。
   1. 访问脚本文件夹:
      cd /home/transcriptome-pipeline/scripts
   2. 执行以下命令以下载参考人类基因组:
      巴什·downloadGenome.sh
   3. 要下载注释，请执行以下命令:
      bash downloadAnnotation.sh
9. 更改参考基因组的注释或版本。
   1. 打开 downloadAnnotation.sh 并 downloadGenome.sh 以更改每个文件的 URL。
   2. 将 downloadAnnotation.sh 和 downloadGenome.sh 文件复制到传输区域，然后在本地操作系统中进行编辑。
      cd /home/transcriptome-pipeline/scripts
      cp downloadAnnotation.sh downloadGenome.sh /opt/transferdata
   3. 打开" 主机路径" 文件夹，该文件夹在步骤 1.5.4 中选择用于在主机和 Docker 容器之间进行链接。
   4. 使用首选编辑器软件编辑文件并保存。最后，将修改后的文件放入脚本文件夹中。执行命令:
      cd /opt/transferdata
      cp downloadAnnotation.sh downloadGenome.sh /home/transcriptome-pipeline/scripts
      注意:这些文件可以使用 vim 或 nano Linux 编辑器直接编辑。
10. 接下来，使用命令行配置 fastq-dump 工具:
    vdb-config --interactive
    注意:这允许从示例数据下载排序文件。
    1. 使用 Tab 键导航 "工具" 页面，然后选择当前文件夹选项。导航到 "保存 "选项，然后单击" 确定"。然后， 退出 快速转储工具。
11. 开始下载以前发表的论文^{中的读数7}。每个样本的 SRA 加入编号是必需的。从SRA NCBI网站（材料表）获取SRA编号。
    注意:要分析公共数据库中可用的RNA-Seq数据，请遵循步骤1.12。要分析私有RNA-seq数据，请遵循步骤1.13。
12. 分析特定的公共数据。
    1. 访问国家生物技术信息中心（NCBI）网站并查找特定主题的关键字。
    2. 单击"基因组"部分中的"生物项目的结果"链接。
    3. 选择并单击特定研究。单击 "SRA 实验"。将打开一个新页面，其中显示了可用于此研究的所有样本。
    4. 单击加入号上方的 "发送到:" 。在 "选择目标" 选项中，选择" 文件 和 格式" 选项，然后选择" RunInfo"。单击 "创建文件" 以导出所有库信息。
    5. 将 SraRunInfo.csv 文件保存在 1.5.4 步骤中定义的主机路径中，然后执行下载脚本:
       cp /opt/transferdata/SraRunInfo.csv /home/transcriptome-pipeline/data
       cd /home/transcriptome-pipeline/scripts
       bash downloadAllLibraries.sh
13. 分析私有和未发布的测序数据。
    1. 将排序数据组织到名为 Reads 的文件夹中。
       注意:在 "读取" 文件夹中，为每个示例创建一个文件夹。对于每个示例，这些文件夹必须具有相同的名称。将每个示例的数据添加到其目录中。如果它是一个配对端的RNA-Seq，每个样本目录应该包含两个FASTQ文件，这些文件必须分别显示根据模式{sample}_1.fastq.gz和{sample}_2.fastq.gz，正向和反向序列结束的名称。例如，名为"Healthy_control"的示例必须具有同名的目录和名为 Healthy_control_1.fastq.gz 和 Healthy_control_2.fastq.gz 的 FASTQ 文件。但是，如果文库测序是单端策略，则只应保存一个读取文件以供下游分析。例如，同一示例"正常对照"必须具有名为 Healthy_control.fastq.gz 的唯一 FASTQ 文件。
    2. 创建包含所有样本名称的表型文件:将第一列命名为"样本"，将第二列命名为"类"。用样本名称填充"样本"列，样本名称必须与样本目录的名称相同，并在"类"列中填充每个样本的表型组（例如，对照组或感染）。最后，保存一个名为"metadata.tsv"的文件，并将其发送到 /home/transcriptome-pipeline/data/ 目录。查看现有的 metadata.tsv 以了解表型文件的格式。
       cp /opt/transferdata/metadata.tsv
       /home/transcriptome-pipeline/data/metadata.tsv
    3. 访问步骤 1.5.4 中定义的 主机路径 目录，并复制新的结构化目录示例。最后，将示例从 /opt/transferdata 移动到管道数据目录。
       cp -rf /opt/transferdata/reads/*
       /首页/转录组-管道/数据/读取/
14. 观察所有读取都存储在文件夹 /home/transcriptome-pipeline/data/reads 中。

2. 数据的质量控制

注:以图形方式评估序列读取中出错的概率。删除所有技术序列，例如适配器。

1. 使用 FastQC 工具访问文库的测序质量。
   1. 要生成质量图表，请运行 fastqc 程序。执行命令:
      巴什·FastQC.sh
      注意:结果将保存在/home/transcriptome-pipeline/results/FastQC文件夹中。由于序列适配器用于文库制备和测序，因此在某些情况下，适配器序列的片段可能会干扰映射过程。
2. 删除适配器序列和低质量读数。访问 "脚本" 文件夹并执行 Trimmomatic 工具的命令:
   cd /home/transcriptome-pipeline/scripts
   bash trimmomatic.sh
   注:用于测序滤波的参数为:去除前导低质量或3个碱基（低于质量3）（领先:3）;删除尾随的低质量或3个基数（低于质量3）（尾随:3）;使用4底座宽的滑动窗口扫描读数，当每个底座的平均质量降至20以下时进行切割（滑动窗口:4:20）;Drop的读数低于36个碱基长度（MINLEN:36）。可以通过编辑 Trimmomatic 脚本文件来更改这些参数。
   1. 确保结果保存在以下文件夹中:/home/transcriptome-pipeline/results/trimreads。执行命令:
      ls /home/transcriptome-pipeline/results/trimreads

3. 样本的映射和注释

注意:在获得高质量的读数后，需要将这些读数映射到参考基因组。对于此步骤，使用 STAR 映射器来映射示例示例。STAR 映射器工具需要 32 GB 的 RAM 内存来加载和执行读取和基因组图谱。对于没有 32 GB RAM 内存的用户，可以使用已映射的读取。在这种情况下，请跳转到步骤 3.3 或使用 Bowtie2 映射器。此部分包含 STAR（结果显示在所有图中）和 Bowtie2（低内存所需的映射器）的脚本。

1. 首先为参考基因组编制索引以进行映射过程:
   1. 使用命令行访问 "脚本" 文件夹:
      cd /home/transcriptome-pipeline/scripts
   2. 对于 STAR 映射器，请执行:
      bash indexGenome.sh
   3. 对于 Bowtie 映射器，请执行:
      巴什 indexGenomeBowtie2.sh
2. 执行以下命令将过滤后的读数（从步骤2中获得）映射到参考基因组（GRCh38版本）。STAR 和 Bowtie2 映射器都是使用默认参数执行的。
   1. 对于 STAR 映射器，请执行:
      bash mapSTAR.sh
   2. 对于 Bowtie2 映射器，执行:
      bash mapBowtie2.sh
      注意:最终结果是存储在 /home/transcriptome-pipeline/results/mapreads 中的每个样本的二进制对齐映射 （BAM） 文件。
3. 使用 FeatureCounts 工具对映射的读数进行注释，以获取每个基因的原始计数。运行为读取内容添加批注的脚本。
   注意:特征计数工具负责将映射的测序读取分配给基因组特征。在生物学问题之后可以更改的基因组注释的最重要方面包括，检测同种型，多个映射读取和外显子 - 外显子连接，对应于参数，GTF.attrType="gene_name"基因或不指定元特征水平的参数，允许MultiOverlap=TRUE和juncCounts=TRUE。
   1. 使用命令行访问脚本文件夹:
      cd /home/transcriptome-pipeline/scripts
   2. 要注释映射的读数以获取每个基因的原始计数，请执行命令行:
      Rscript 注释。R
      注:用于注释过程的参数为:返回基因短名称（GTF.attrType="gene_name"）;允许多个重叠（允许多重叠 = TRUE）;并指示库已配对端 （isPairedEnd=TRUE）。对于单端策略，请使用参数 isPairedEnd=FALSE。结果将保存在/home/transcriptome-pipeline/countreads文件夹中。
4. 使基因表达正常化。
   注意:归一化基因表达对于比较结果（例如，健康和感染样本）至关重要。还需要归一化来执行微扰分析的共表达和分子程度。
   1. 使用命令行访问 "脚本" 文件夹:
      cd /home/transcriptome-pipeline/scripts
   2. 使基因表达正常化。为此，请执行命令行:
      Rscript 规范化样本。R
      注意:在此实验中，原始计数表达式使用 M 值的修剪平均值 （TMM） 和每百万计数 （CPM） 方法进行归一化。本步骤旨在通过进行文库大小归一化来消除由于技术影响而导致的基因表达差异。结果将保存在/home/transcriptome-pipeline/countreads文件夹中。

4. 差异表达基因和共表达基因

1. 使用开源 EdgeR 包鉴定差异表达的基因。这涉及找到与对照相比表达更高或更低的基因。
   1. 使用命令行访问 "脚本" 文件夹:
      cd /home/transcriptome-pipeline/scripts
   2. 要识别差异表达的基因，请使用命令行执行DEG_edgeR R 脚本:
      Rscript DEG_edgeR.R
      注意:包含差异表达基因的结果将保存在/home/transcriptome-pipeline/results/degs文件夹中。数据可以传输到个人计算机。
2. 从 csblusp/transcriptome 容器下载数据。
   1. 将已处理的数据从 /home/transcriptome 管道传输到 /opt/transferdata 文件夹（本地计算机）。
   2. 通过执行命令行将所有文件复制到本地计算机:
      cp -rf /home/transcriptome-pipeline/results /opt/transferdata/pipeline
      cp -rf /home/transcriptome-pipeline/data /opt/transferdata/pipeline
      注意:现在，转到本地计算机，以确保所有结果、数据集和数据都可以在主机路径中下载。
3. 识别共表达模块。
   1. 访问共表达模块识别工具 （CEMiTool） 网站 （表格
      材料）。此工具从用户提供的表达式数据集中识别共表达式模块。在主页上，单击右上角的 "运行 "。这将打开一个新页面来上传表达式文件。
   2. 单击"表达文件"部分下方的"选择文件"，然后从主机路径上传归一化的基因表达矩阵"tmm_expression.tsv"。
      注意:步骤 4.4.是非强制性的。
4. 探索共表达模块的生物学意义。
   1. 单击"示例表型"部分中的"选择文件"，然后从"下载数据"步骤 4.2.2 metadata_cemitool.tsv 上传包含示例表型的文件。执行基因集富集分析 （GSEA）。
   2. 按"基因相互作用"部分中的"选择文件"以上传具有基因相互作用的文件（cemitool-interactions.tsv）。可以使用webCEMiTool提供的基因相互作用文件作为示例。相互作用可以是蛋白质 - 蛋白质相互作用，转录因子及其转录基因或代谢途径。此步骤为每个共表达模块生成一个交互网络。
   3. 单击"基因集"部分中的"选择文件"，以上传基因基质转置（GMT）格式文件中功能相关的基因列表。基因集文件使工具能够对每个共表达模块执行富集分析，即过度表示分析（ORA）。
      注意:此基因列表可以包含通路、GO 项或 miRNA 靶基因。研究人员可以使用血液转录模块（BTM）作为该分析的基因集。BTM 文件 （BTM_for_GSEA.gmt）。
5. 设置用于执行共表达分析的参数并获得其结果。
   1. 接下来，通过单击加号展开" 参数 "部分以显示默认参数。如有必要，请更改它们。选中 应用 VST 复选框。
   2. 在" 电子邮件 "部分中编写电子邮件，以电子邮件形式接收结果。此步骤是可选的。
   3. 按 "运行 CEMiTool" 按钮。
   4. 点击右上角的下载完整报告， 下载完整的分析报告 。它将cemitool_results.zip下载压缩文件。
   5. 使用 WinRAR 提取cemitool_results.zip的内容。
      注意:包含提取内容的文件夹包含多个文件，其中包含所有分析结果及其建立的参数。

5. 样品分子扰动程度的测定

1. 分子扰动度 （MDP） 网络版。
   1. 要运行 MDP，请访问 MDP 网站（材料表）。MDP计算每个样品与参考的分子距离。单击" 运行 "按钮。
   2. 在 "选择文件" 链接上，将表达式文件tmm_expression.tsv。然后，从下载数据步骤 4.2.2 上传表型数据文件 metadata.tsv。也可以提交GMT格式的通路注释文件，以计算与疾病相关的通路的扰动评分。
   3. 上传数据后，定义包含 MDP 使用的表型信息的"类"列。然后，通过选择与控件类对应的标签来定义控件类。
      注意:有一些可选参数会影响示例分数的计算方式。如有必要，用户可以更改统计平均方法，标准偏差和扰动基因的最高百分比。
   4. 之后，按 "运行MDP "按钮，将显示MDP结果。用户可以通过单击每个 图中的"下载图" 以及"下载 MDP 分数 文件"按钮上的 MDP 分数来下载 数字。
      注意:如果对如何提交文件或 MDP 的工作原理有疑问，只需浏览教程和关于网页即可。

6. 功能丰富分析

1. 创建一个下调 DG 列表和另一个上调 DEX 列表。基因名称必须根据Entrez基因符号。列表中的每个基因必须放在一行上。
2. 以 txt 或 tsv 格式保存基因列表。
3. 访问 Enrichr 网站（材料表）以执行功能分析。
4. 通过单击" 选择文件"选择基因列表。选择其中一个 DEG 列表，然后按 "提交 "按钮。
5. 单击网页顶部的 Pathways ，使用 ORA 方法执行功能富集分析。
6. 选择路径数据库。"Reactome 2016"通路数据库被广泛用于获取人类数据的生物学意义。
7. 再次单击路径数据库的名称。选择 条形图 并检查是否按 p 值排名排序。如果没有，请单击条形图，直到按 p 值排序。此条形图包括根据 p 值的前 10 个路径。
8. 按 配置 按钮，选择红色进行上调基因分析，选择蓝色进行下调基因分析。通过单击 svg、 png 和 jpg，以多种格式保存条形图。
9. 选择 表 ，然后单击条形图左下角 的将条目导出到表中 ，以获取txt文件中的功能丰富分析结果。
   注意:此功能富集结果文件在每行中包含一个通路的名称、提交的DEG列表和通路之间的重叠基因数、p值、调整后的p值、优势比、组合评分以及DEG列表中存在的参与通路的基因的基因符号。
10. 对其他 DEG 列表重复相同的步骤。
    注意:使用下调DEG的分析为下调基因提供了富集的途径，使用上调基因的分析为上调基因提供了富集的途径。

Results

转录组分析的计算环境是在 Docker 平台上创建和配置的。这种方法允许初学者Linux用户在没有先验管理知识的情况下使用Linux终端系统。Docker 平台使用主机操作系统的资源来创建包含特定用户工具的服务容器（图 1B）。创建了一个基于Linux OS Ubuntu 20.04发行版的容器，并完全配置为转录组分析， 可通过 命令行终端访问。在此容器中，有一个用于数据集和脚本的预定义文件夹结构，这是所有管道分析所必需的（图 1C）。我们的研究小组⁷ 发表的一项研究用于分析，其中包括来自健康个体的20个样本和来自CHIKV急性感染个体的39个样本（图1D）。

总RNA测序过程可能会产生读取错误，这可能是由具有两个或多个转录本的簇或试剂耗尽引起的。测序平台返回一组"FASTQ"文件，其中包含每个核苷酸碱基的序列（读取）和相关质量（

Discussion

测序文库的制备是以最佳方式回答生物学问题的关键一步。该研究感兴趣的转录本类型将指导选择哪种类型的测序文库并推动生物信息学分析。例如，根据测序的类型，从病原体和宿主相互作用的测序中，可以识别来自两者的序列或仅来自宿主转录本的序列。

下一代测序设备，例如Illumina平台，测量测序质量评分，这代表了碱基被错误调用的概率。下游分析对低质量序列非常敏感，并导致基因表达不足或误读。执行正确分析和解释的另一个障碍是适配器序列。适配器序列有助于文库的制备和测序，在大多数情况下，适配器也是测序的。最近的研究发现，制图工具对最终结果的影响很小¹³。然而，在病原体 - 宿主研究中，当测试不同的阈值以最小化多映射位点序列问题时，映射过程可以产生稍微更好的结果。

差异基因表达结果应谨慎解释，特别是当每组样本数量非常少且样本来自不同的测定并受到批次效应干扰DEG结果时。这些结果对几个因素很敏感:（一）所应用的数据过滤，例如去除低表达基因和要维持的样本数量;（ii）研究设计，仅比较样本组或每个感染患者与所...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
CEMiTool	计算系统生物学实验室	1.12.2	以全自动方式发现和分析共表达基因模块，同时提供带有高质量图表的用户友好型 HTML 报告。
EdgeR	Bioconductor（维护者：Yunshun Chen [yuchen at wehi.edu.au]）	3.30.3	具有生物复制的 RNA-seq 表达谱的差异表达分析
EnhancedVolcano	Bioconductor（维护者：Kevin Blighe [kevin at clinicalbioinformatics.co.uk]）	1.6.0	具有增强着色和标记的可发表的火山图
FastQC	Babraham Bioinformatics	0.11.9	旨在提供一种简单的方法来对来自高通量测序的原始序列数据进行一些质量控制检查
FeatureCounts	生物信息学部，Walter 和 Eliza Hall 医学研究所	2.0.0	将映射测序读数分配给指定的基因组特征
MDP	计算系统生物学实验室	1.8.0	扰动分子度 根据来自对照
R	R 核心组	4.0.3	计算统计计算和图形的编程语言和免费软件环境
STAR	生物信息学部，沃尔特和伊丽莎霍尔医学研究所	2.7.6a	Aligner 旨在专门解决 RNA-seq 数据映射的许多挑战，使用一种策略来解释剪接比对
Bowtie2	约翰霍普金斯大学	2.4.2	用于将测序读数与长参考序列比对的超快速且节省内存的工具
Trimmomatic	THE USADEL LAB	0.39	修剪 Illumina 双端和单端数据的 Trimming 接头序列任务
获取 Docker	Docker	20.10.2	创建可重复且可预测的 （https://docs.docker.com/get-docker/）
WSL2 内核	Windows	NA	https://docs.microsoft.com/en-us/windows/wsl/wsl2-kernel
获取 Docker Linux	Docker	NA	https://docs.docker.com/engine/install/ubuntu/
Docker Linux 存储库	Docker	NA	https://docs.docker.com/engine/install/ubuntu/#install-using-the-repository
MDP 网站	计算系统生物学实验室	NA	https://mdp.sysbio.tools
Enrichr 网站	MaayanLab	NA	https://maayanlab.cloud/Enrichr/
webCEMiTool	计算系统生物学实验室	NA	https://cemitool.sysbio.tools/
gProfiler	生物信息学、算法学和数据挖掘小组	NA	https://biit.cs.ut.ee/gprofiler/gost
goseq	Bioconductor（维护者：Matthew Young [my4 at sanger.ac.uk]）	NA	http://bioconductor.org/packages/release/bioc/html/goseq.html
SRA NCBI 研究	NCBI	NA	https://www.ncbi.nlm.nih.gov/bioproject/PRJNA507472/