Name: Rup（RNA-seq 可用性评估管道）- 真核生物中批量 RNA-seq 实验的质量控制 (视频) | JoVE Journal
Uploaded: 2025-11-07T13:30:00.000+00:00
Duration: 5 min 7 s
Description: 学习使用Rup管道进行真核生物批量RNA-seq质量控制，通过自动化工具评估测序数据可靠性。该方法帮助生物学家在基因表达分析前识别低质量读段和标记错误，确保转录组学研究的可重复性。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

我们感谢吉森尤斯图斯-李比希大学系统生物学教授职位下属的生物信息学核心设施提供的技术协助，以及de.NBI网络内的BiGi服务中心（BMFB资助号031A533）提供的计算资源和常规支持。本研究所获经费来自德国研究基金会（DFG）资助项目BE2547/24-1（授予A.B.），同时我们亦感谢德国吉森尤斯图斯-李比希大学的支持。

<ol><li>Kivivirta, K. I., Herbert, D., Roessner, C., De Folter, S., Marsch-Martinez, N., Becker, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptome+analysis+of+gynoecium+morphogenesis+uncovers+the+chronology+of+gene+regulatory+network+activity%5BTitle%5D)+AND+%22Plant+Physiol%22%5BJournal%5D)">Transcriptome analysis of gynoecium morphogenesis uncovers the chronology of gene regulatory network activity.</a> Plant Physiol. 185 (3), 1076-1090 (2021).</li>
<li>Hohenfeld, C. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparative+analysis+of+infected+cassava+root+transcriptomics+reveals+candidate+genes+for+root+rot+disease+resistance%5BTitle%5D)+AND+%22Sci+Rep%22%5BJournal%5D)">Comparative analysis of infected cassava root transcriptomics reveals candidate genes for root rot disease resistance.</a> Sci Rep. 14 (1), 10587(2024).</li>
<li>Li, Q., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Time-course+transcriptomic+information+unravels+the+mechanisms+of+improved+drought+tolerance+by+drought-priming+in+wheat%5BTitle%5D)+AND+%22J+Integr+Agric%22%5BJournal%5D)">Time-course transcriptomic information unravels the mechanisms of improved drought tolerance by drought-priming in wheat.</a> J Integr Agric. 24 (8), 2902-2919 (2024).</li>
<li>DeLuca, D. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA-SeQC%3A+RNA-seq+metrics+for+quality+control+and+process+optimization%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">RNA-SeQC: RNA-seq metrics for quality control and process optimization.</a> Bioinformatics. 28 (11), 1530-1532 (2012).</li>
<li>Wang, L., Wang, S., Li, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RSeQC%3A+Quality+control+of+RNA-seq+experiments%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">RSeQC: Quality control of RNA-seq experiments.</a> Bioinformatics. 28 (16), 2184-2185 (2012).</li>
<li>Zhou, Q., Su, X., Jing, G., Chen, S., Ning, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA-QC-chain%3A+Comprehensive+and+fast+quality+control+for+RNA-seq+data%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">RNA-QC-chain: Comprehensive and fast quality control for RNA-seq data.</a> BMC Genomics. 19 (1), 144(2018).</li>
<li>Deshpande, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA-seq+data+science%3A+From+raw+data+to+effective+interpretation%5BTitle%5D)+AND+%22Front+Genet%22%5BJournal%5D)">RNA-seq data science: From raw data to effective interpretation.</a> Front Genet. 14, 997383(2023).</li>
<li>Li, D., Zand, M. S., Dye, T. D., Goniewicz, M. L., Rahman, I., Xie, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+evaluation+of+RNA-seq+differential+analysis+methods%5BTitle%5D)+AND+%22PLoS+One%22%5BJournal%5D)">An evaluation of RNA-seq differential analysis methods.</a> PLoS One. 17 (9), e0264246(2022).</li>
<li>FastQC: A quality control tool for high throughput sequence data. , Babraham Bioinformatics. 
 <a target="_blank" href="http://www.bioinformatics.babraham.ac.uk/projects/fastqc/">http://www.bioinformatics.babraham.ac.uk/projects/fastqc/</a> (2025).</li>
<li>Kassambara, A. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aKassambara%2C+A+%22fastqcr%3A+Quality+control+of+sequencing+data%22">fastqcr: Quality control of sequencing data</a>. , https://rpkgs.datanovia.com/fastqcr/index.html (2023).</li>
<li>Chen, S., Zhou, Y., Chen, Y., Gu, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((fastp%3A+An+ultra-fast+all-in-one+FASTQ+preprocessor%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">fastp: An ultra-fast all-in-one FASTQ preprocessor.</a> Bioinformatics. 34 (17), i884-i890 (2018).</li>
<li>Wang, W., Luo, J. D., Carroll, T. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aWang%2C+W+%22Rfastp%22">Rfastp</a>. , https://www.bioconductor.org/packages/release/bioc/html/Rfastp.html (2025).</li>
<li>Wickham, H. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aWickham%2C+H+%22ggplot2%3A+Elegant+graphics+for+data+analysis%22">ggplot2: Elegant graphics for data analysis</a>. , Springer-Verlag. New York. https://ggplot2.tidyverse.org (2016).</li>
<li>Shi, W. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aShi%2C+W+%22Rsubread%22">Rsubread</a>. , https://bioconductor.org/packages/release/bioc/html/Rsubread.html (2025).</li>
<li>Wagner, G. P., Kin, K., Lynch, V. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Measurement+of+mRNA+abundance+using+RNA-seq+data%3A+RPKM+measure+is+inconsistent+among+samples%5BTitle%5D)+AND+%22Theory+Biosci%22%5BJournal%5D)">Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples.</a> Theory Biosci. 131 (4), 281-285 (2012).</li>
<li>Liu, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Evaluating+the+impact+of+sequencing+depth+on+transcriptome+profiling+in+human+adipose%5BTitle%5D)+AND+%22PLoS+One%22%5BJournal%5D)">Evaluating the impact of sequencing depth on transcriptome profiling in human adipose.</a> PLoS One. 8, e66883(2013).</li>
<li>Schurch, N. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((How+many+biological+replicates+are+needed+in+an+RNA-seq+experiment+and+which+differential+expression+tool+should+you+use%5BTitle%5D)+AND+%22RNA%22%5BJournal%5D)">How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use.</a> RNA. 22 (6), 839-851 (2016).</li>
<li>Liu, Y., Zhou, J., White, K. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA-seq+differential+expression+studies%3A+More+sequence+or+more+replication%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">RNA-seq differential expression studies: More sequence or more replication.</a> Bioinformatics. 30 (3), 301-304 (2014).</li>
<li>Deschamps-Francoeur, G., Simoneau, J., Scott, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Handling+multi-mapped+reads+in+RNA-seq%5BTitle%5D)+AND+%22Comput+Struct+Biotechnol+J%22%5BJournal%5D)">Handling multi-mapped reads in RNA-seq.</a> Comput Struct Biotechnol J. 18, 1569-1576 (2020).</li>
<li>Seppey, M., Manni, M., Zdobnov, E. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((BUSCO%3A+Assessing+genome+assembly+and+annotation+completeness%5BTitle%5D)+AND+%22Methods+Mol+Biol%22%5BJournal%5D)">BUSCO: Assessing genome assembly and annotation completeness.</a> Methods Mol Biol. 1962, 227-245 (2019).</li>
<li>Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., Kingsford, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Salmon+provides+fast+and+bias-aware+quantification+of+transcript+expression%5BTitle%5D)+AND+%22Nat+Methods%22%5BJournal%5D)">Salmon provides fast and bias-aware quantification of transcript expression.</a> Nat Methods. 14 (4), 417-419 (2017).</li>
</ol>

作者声明无利益冲突。

差异基因表达分析的质量高度依赖于两个因素：每个样本及重复中测序的读段数量16，以及每个样本的重复次数17,18。本文介绍了用户友好的Rup分析流程，用于确定每个重复中适用于基因表达定量的读段数量。通过多种评估指标，研究人员可以了解为何某些重复的比对读段数偏低，并识别样本间相关性方面存在的问题。尽管Rup最初是为植物RNA-seq样本的质量评估而开发，但它同样适用于其他真核生物，且在用于非植物样本时无需额外调整（补充文件1）。

Rup 的第一步是确定读段修剪前后总的读段数量。测序读段的数量决定了可检测基因的数量，如果读段数量过低，许多差异表达的基因将无法被检测到。在质量修剪和过滤过程中被丢弃的读段占比较大时，可能表明 RNA 存在降解，这可能导致对高度降解基因的表达水平被低估。然而，一个样本是否可用于后续分析，取决于目标生物体以及研究目的。例如，根据我们的经验，要检测大多数植物基因的表达，需要 3000 万至 5000 万条读段，而对于真菌，可能仅需 1000 万条读段就已足够。因此，Rup 不会设定排除问题样本的质量阈值，而是提供一系列指标，用于识别可能遇到的各种问题。

第二步（比对质量）评估测序读段与参考基因组比对的准确性和可靠性，以确定其是否适用于后续分析。并非所有测序得到的读段都可用于基因丰度计算；在大多数情况下，无法比对到基因组或转录组的读段，或比对到基因组多个位置的读段将被忽略。19 且不计入总体读段计数。分析流程第二步的结果可用于理解为何某些读段未被用于丰度计算。未比对读段数量较高可能表明RNA提取过程中存在污染（例如植物病原体或植食性生物的污染）或参考基因组不完整。参考基因组中基因模型不完整可能导致大量“无特征”读段。多重比对读段占比较高可能源于文库中rRNA基因数量过多，提示在测序文库构建过程中rRNA去除不充分（若这些多重比对读段确实来自rRNA，可通过向分析流程提供rRNA注释文件以进一步分析，流程将自动计算可能的rRNA读段数量）。

第三个同样重要的质量指标是同一样本重复之间的相关性。通常情况下，同一样本重复之间的相关性应高于不同样本重复之间的相关性。重复样本之间相关性较低可能表明重复间存在较大的生物学变异、样本/条件/组织之间高度相似、其他批次效应，甚至样本混淆或标签错误。Rup 可计算所有样本间的成对相关性，并生成样本相关性的聚类热图。此外，还可对样本进行主成分分析（PCA），以识别在后续下游分析中需要关注的潜在批次效应。尽管批次效应可在后续分析中进行校正，但如果检测到的差异过大，更合适的做法可能是剔除存在问题的样本。

起始材料直接影响测序质量：组织取样、RNA提取和文库构建是减少RNA-seq质量问题的关键步骤。应尽可能保持一致的实验条件，例如使用生长培养箱。应及时采取病虫害防控措施，仅选择健康个体进行取样。建议在同一天和/或同一时间点采集样品，以降低昼夜节律引起的转录组变异。市面上有多种RNA提取试剂盒，针对目标物种选择合适的试剂盒可提高mRNA质量。文库构建方案中应包含polyA+ RNA富集步骤，以最大限度减少rRNA比例。

Rup 可作为差异基因表达分析中初步的质量控制步骤。该质量控制因多个关键原因而必需，因为它可确保输入数据的质量（识别低质量的重复样本/样品，设定测序深度和比对率的质量阈值），并发现诸如测序错误、批次效应或样本标签错误等技术问题。如果输入的 RNA 质量足够高，Rup 可通过修剪低质量的测序读段或识别标签错误的样本提供帮助。然而，Rup 无法弥补输入 RNA 质量差的问题。在 RNA 质量较低或测序数据存在错误的情况下，可能需要重新采集样本、调整 RNA 提取方案和/或重复测序。对于因污染导致大量未比对读段的样本，同样适用此处理原则。然而，根据我们的经验，并非所有 RNA 提取都能因材料不足而简单重复。在此情况下，这些样本不适用于某些下游应用：含有较少唯一比对读段的样本不应用于差异基因表达分析，但仍可用于转录本存在性分析。此时，在组织/处理/条件下转录本的缺失情况不能纳入转录本丰度的分析与定量中。

Rup 存在一些局限性。例如，它不能直接检测 RNA 降解，因为在文库构建和测序之前需要对 RNA 完整性进行直接测定。然而，本流程的代码仓库中包含一个使用 RSeQC 模块 geneBodyCoverage.py 和 tin.py 来识别 RNA 降解的脚本。此外，Rup 不检测 GC 含量和转录本长度偏差。目前参考基因组的大小限制为 4 Gb，根据我们的经验，读段比对模块在多倍体中会高估多重比对的读段数量，这一点在补充文件中通过 E. californica 单倍体与二倍体基因组版本的读段比对比较得到了证实，因此该流程更推荐使用单倍体基因组作为输入。 Rup 的输出质量在很大程度上依赖于参考基因组和注释的质量，不完整或高度碎片化的基因组序列可能导致未比对读段数量被高估。此外，基因模型注释不完整会导致未分配读段数量被高估。基因组序列和基因注释的完整性与重复性可通过 BUSCO20 等工具进行评估。

Rup 是一个独立工具，适用于 RNA-seq 实验中所有必需的初始质量控制步骤。与 RSeQC 和 RNA-SeQC 不同，Rup 无需对原始测序读段进行预处理。RNA-SeQC 无法进行测序质量分析，而 RSeQC 和 RNA-SeQC 均不计算样本相关性热图（表 1）。此外，RSeQC 的标准化输出为 FPKM 格式，且 RSeQC 和 RNA-SeQC 均未将输出可视化作为所有模块的标准功能。因此，这两种替代工具未能检测若干质量控制问题，包括读段数量过低、修剪后读段比例过高以及重复样本离群值的识别。Rup、RSeQC 和 RNA-SeQC 的比较详见 补充文件 1。此外，Rup 可与 RNA-SeQC 或 RSeQC 互补使用；例如，本流程生成的已排序 BAM 文件可作为这些工具的输入文件。

目前，Rup 已针对少量样本进行了优化，但最耗时的步骤（如质量修剪和读段比对）可以预先计算，例如在计算机集群或云基础设施上完成。对于大于 4 Gb 的基因组，这一步是必需的，因为 Rsubread 比对工具仅支持小于 4 Gb 的基因组。rRNA 注释由 barrnap 完成，该工具可识别高度保守的 rRNA 基因。根据我们的经验，rRNA 基因注释无需完全覆盖，因为即使某些异常的 rRNA 基因缺失，仅对数据集中 rRNA 存在情况的总体概览也足以用于数据集质量评估。未来版本的 Rup 可能会切换至其他比对方法，例如伪比对工具 salmon21，以减少运行时间。此外，还可向 Rup 中添加更多标准化和校正方法，如 GC 偏好性或长度偏好性校正。

总之，Rup 提供了确保 RNA-seq 数据分析可靠性和可重复性的关键信息。该流程全面报告了主要的 RNA-seq 质控指标，并生成可直接用于下游分析的输出文件。它被设计为一个独立工具，便于生物信息学知识有限的研究人员通过直观的可视化方式评估其原始测序数据的质量。

转录组学（RNA-seq）实验可全面分析塑造表型的转录特征。该方法在植物分子遗传学中已成为不可或缺的工具，用于鉴定单个或共表达的基因以及参与特定生物学过程的通路，例如发育、植物与病原体互作或非生物胁迫抗性1, 2 ,3。RNA-seq 技术的最新进展提高了检测特异性，并能够在单碱基分辨率下识别不同的转录本异构体和变异形式，从而可鉴定从较大的插入缺失到单核苷酸多态性（SNPs）的序列变异。RNA-seq 所获得的数据具有广泛的动态范围，有助于同时检测高丰度和低表达水平的转录本，但需要适当的实验条件以确保数据一致性。此外，RNA-seq 数据集规模庞大，分析和存储过程对计算资源要求较高，因此需要高效的数据管理和充足的计算能力。RNA-seq 在整个实验流程的每个步骤均需要高质量的起始材料，因为任一环节的不足都可能逐级放大并影响最终结果。RNA 完整性差或文库构建过程中的技术问题将导致偏差并降低结果准确性。获取高质量原始测序读段的参数可能因不同的高通量测序（NGS）平台而异。
根据我们的经验，我们建议仅使用RNA完整性数值（RIN）高于7的RNA作为RNA-seq文库构建的起始材料，因为该数值表明mRNA结构基本完整。为确保测序成功，标准文库构建方案通常需要约2 µg总RNA，浓度为50–200 ng/µL。应通过分光光度计检测确认RNA纯度，OD260/280比值应在1.8至2.1之间，OD260/230比值应大于1.5。测序深度不足、比对率偏低或与参考基因组错配可能会进一步扭曲基因表达和剪接的定量结果。此外，实验设计不当，例如未能区分不同组织或处理组的样本，或重复样本间变异过大，均可能引入噪声并降低可重复性。尽管许多实验室常规开展转录组分析，但对最初关键分析步骤的严格质量控制常常未在发表文献中报告，甚至完全缺失。这可能导致对转录组分析结果的过度解读，进而造成结果不可重复。
本文提供了一套用于高质量转录组分析中mRNA表达谱初始步骤的质量控制工作流程，旨在检测转录水平的变化。我们的目标是使湿实验生物学家即使在生物信息学知识有限的情况下，也能够评估其原始转录组数据。Rup（图1）适用于具备R语言基础知识的研究人员。执行本文所述的工作流程，将使研究人员能够深入理解其原始数据及其在后续分析中可能存在的局限性。据我们所知，目前尚缺乏将实用的原始转录组数据评估流程与区分高质量与低质量数据的指导原则相结合的系统方法。
<img alt="包含质量控制方法、输入文件来源和数据可视化的RNA测序流程图。" class="xfigimg" src="/files/ftp_upload/69253/69253fig1.jpg" data-altupdatedat="1783396841000" data-alt="图1"> 
图1：RNA测序 in silico 质量控制流程的工作流程图。 质量控制的输入文件来源于对植物材料进行测序所获得的RNA测序数据，以及公开可用的数据集。所提供的R语言流程通过三个主要方面评估序列质量：测序质量、比对质量和重复样本质量。计算多种质量指标，并利用常见的R语言软件包（如ggplot2和pheatmap）对统计结果进行可视化展示（例如条形图和热图）。<a href="https://www.jove.com/files/ftp_upload/69253/69253fig1large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
既往报道的评估流程需要预处理数据（例如，以 .bam 文件形式的测序读段比对结果），未能涵盖所有评估指标，或已不再维护4,5,6。本文所展示工作流程的优势在于其全面性，能够将最常见的质量控制问题整合到单一分析流程中。此外，我们提供了适用于所有下游分析应用的高质量数据示例，同时也提供了低质量数据的示例，并讨论了这些数据在进一步分析中的具体局限性。已有若干先前发表的研究阐述了RNA-seq分析工具的目的，并对其性能进行了比较评估7,8。然而，RNA-seq的质量控制和方法学报告标准尚未统一，这加剧了转录组学实验在可重复性以及生物学意义解读方面的困难。
Rup 整合了标准的高质量工具、质量控制分析以及结果可视化功能。Rup 需要以原始测序读段和注释好的基因组作为输入，可在 Mac OS、Linux 以及 Windows 系统上的“Windows 子系统 for Linux”（WSL）上运行。该工具集成了用于评估测序质量的程序，可量化测序读段在基因组上的比对情况，包含对样本中 rRNA 含量的测定，并提供样本相关性分析以评估重复样本之间的一致性。测试数据是通过激光显微切割技术获取植物物种Eschscholzia californica两个不同发育阶段的中央分生组织组织，采用超低起始量文库构建方案，并在 Novaseq 6000 平台上测序获得的。

<table><tbody><tr><td>fastqcr</td><td>R</td><td>0.1.3</td><td></td></tr><tr><td>FUJITSU ESPRIMO D958 台式机</td><td>FUJITSU</td><td></td><td>该分析流程在 Ubuntu Linux 24.02 系统（32 Gb 内存）上开发和测试</td></tr><tr><td>getopt</td><td>R</td><td>1.20.4</td><td></td></tr><tr><td>ggplot2</td><td>R</td><td>3.5.2</td><td></td></tr><tr><td>MacBook Pro&amp;nbsp;</td><td>Apple</td><td></td><td>该分析流程在 macOS 15.4.1 系统（16 Gb 内存）上测试</td></tr><tr><td>pheatmap</td><td>R</td><td>1.0.13</td><td></td></tr><tr><td>R</td><td></td><td>4.4.3</td><td></td></tr><tr><td>reshape2</td><td>R</td><td>1.4.4</td><td></td></tr><tr><td>rfastp</td><td>bioconductor</td><td>1.16.0</td><td></td></tr><tr><td>rsamtools</td><td>bioconductor</td><td>2.22.0</td><td></td></tr><tr><td>rsubread</td><td>bioconductor</td><td>2.20.0</td><td></td></tr></tbody></table>

rup（rna-seq 可用性评估流程）— 真核生物批量 rna-seq 实验的质量控制

Rup (RNA-seq Usability Assessment Pipeline) - Quality Control for Bulk RNA-seq Experiments in Eukaryotes

分子植物科学的现代研究方法通常需要进行大规模 RNA 测序（bulk RNA-seq）实验，例如追踪处理条件下转录组的全局变化，或鉴定调控通路中的关键组分。因此，植物科学的多个领域都依赖高质量且可重复的大规模 RNA-seq 数据来推动科学进步。然而，根据我们的经验，RNA-seq 数据集中质量控制措施的相关知识和实际应用往往不足。本文介绍了 Rup（RNA-seq 可用性评估流程），用于大规模 RNA-seq 数据的质量控制，以支持后续的基因表达分析。Rup 是一个独立运行的工具，适用于具备基础 R 语言知识的实验生物学研究人员。Rup 可帮助区分高质量、适合下游基因表达分析的测序数据与不适合进一步通用分析的数据。Rup 包含多项针对常见问题的检测功能，例如测序读段数量或比对率不足、污染序列的识别、总 RNA-seq 数据中 rRNA 组分的定量、技术重复间的相似性检验，并通过真实数据演示提供直观的可视化结果。Rup 提供了一整套工具，可在标准化转录组分析之前识别实验中的缺陷，从而提升个体研究人员及整个领域的数据质量。该工具增强了对大规模 RNA-seq 数据分析结果的信心，并为未来制定定义最低质量控制标准的指南提供了基础，从而提高已发表 RNA-seq 数据的可靠性和透明度。Rup 及测试数据可在 https://github.com/oliverrupp/rup 获取。

该流程完全以 R 脚本实现，并已在 Linux 和 Mac OS 操作系统上进行测试。Windows 用户可使用 Windows Subsystem for Linux (WSL)。代码和测试数据作为 GitHub 仓库提供：https://github.com/oliverrupp/rup。测序数据可在 ENA EBI 项目 PRJEB96400 下获取。
通过将两个真实样本人工构建为十个样本，以展示在使用 Rup 进行分析的 bulk RNA-seq 质量控制过程中可能遇到的各种问题。样本 s2_r1 被设计为总读段数较低，样本 s2_r2 包含大量低质量读段，将在剪切过程中被剔除。样本 s2_r3 包含大量 rRNA 读段，而样本 s2_r4 包含未能比对到参考基因组的污染读段。样本 s1_r5 和 s2_r5 的名称被互换，以说明重复样本间相关性低的情况。
第2个实验方案部分能够在修剪前或修剪后识别出读数较低的样本（图2）。柱状图显示，样本s1_r1在修剪前后读数均较低，此处初始读数原本就偏低。样本s1_r2在修剪后读数偏低，提示其中含有大量接头序列、测序错误、引物序列或被测序的poly-A/T片段，在修剪过程中被去除。输入RNA的降解也可能导致高质量读段数量减少。
第3部分的实验方案指出了读段分配中存在的问题。图3显示了样本s2_r3（绿色）中多位置比对读段数量升高，以及核糖体RNA（rRNA）读段数量较多的情况。样本s2_r4的污染情况表现为大量读段未能比对到参考基因组（粉色）。这些读段并不对应于rRNA，而是来源于非目标生物体的序列。图4展示了转录组中分配读段数量过低所引发的普遍问题。在比对到参考基因组的读段比例较低的样本（s2_r1至r4）中，仅有约三分之一的基因被分配到超过100个读段；而在其他样本中，约有一半的基因属于此类。在s2_r1至r4样本中，低表达水平基因的读段可能无法被检测到，因此使用这些样本进行的比较表达分析将极不可靠，应予以避免。
第4部分方案可用于识别重复样本中的离群值。相同样本/条件/组织的重复样本之间预期具有比其他样本/条件/组织的重复样本之间更高的相关性。图5展示了两个样本（S1和S2）各自五个重复样本的相关性热图。图上方和侧边的树状图显示了两个聚类簇，每个簇包含五个重复样本。左侧簇包含样本1的四个重复样本和样本2的一个重复样本（s2_r5），右侧簇包含样本2的四个重复样本和样本1的一个重复样本（s1_r5）。在此情况下，若将样本名称s1_r5和s2_r5重新对调，则每个簇均包含单一样本的全部重复样本，这可能提示存在重复样本标记错误。重复样本未能聚在一起的其他原因可能包括样本/条件/组织之间缺乏差异性，或重复样本仅因测序质量因素而聚类。后一种情况可能发生在所有在剪切和比对后具有异常低或异常高读段数的重复样本自行形成一个簇时。
<img alt="相关性热图；通过颜色梯度可视化数据关系以进行模式分析。" class="xfigimg" src="/files/ftp_upload/69253/69253fig5.jpg" data-altupdatedat="1783396841000" data-alt="图 5"> 
图 5：样本相关性热图。该样本相关性热图基于 log2 转换后的 TPM 值构建，显示了两个明显聚类，每个聚类包含五个样本。生物学/技术学重复样本之间的相关性应高于来自其他组织/处理组的样本。左侧聚类包含样本 1 的四个重复和样本 2 的一个重复（s2_r5），右侧聚类包含样本 2 的四个重复和样本 1 的一个重复（s1_r5）。应对这些单个重复样本进行检查，以确认是否存在可能的样本混淆或批次效应，从而解释其在热图中的聚类情况。<a href="https://www.jove.com/files/ftp_upload/69253/69253fig5large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
表1：RNA-seq质量评估流程的比较。 详细分析请参见补充文件1。<a href="https://www.jove.com/files/ftp_upload/69253/69253_Table%201.xlsx" target="_blank">请点击此处下载该表格。</a>
补充文件 1：参数和参考选择。 这些分析展示了参数和参考选择对整体质控结果的影响。<a href="https://www.jove.com/files/ftp_upload/69253/20250912_JoveManuscript_SupplementalFile.docx" target="_blank">请点击此处下载该文件。</a>

本方案可为生物湿实验研究人员（生物信息学经验有限）提供RNA测序实验的初步质量控制。

Rup（RNA-seq 可用性评估流程）— 真核生物批量 RNA-seq 实验的质量控制

现代分子植物科学的许多研究方法通常需要进行大量RNA测序实验，例如，用于追踪转录组在特定处理或发育阶段下的全局变化……

我们开发了创新的生物信息学工具，以简化、自动化并整合高通量实验的数据分析。我们利用高通量测序、先进的生物信息学软件以及强大的计算基础设施，实现系统的生物学分析。首先，使用 Bioconductor 包管理器安装所有必需的 R 软件包。创建一个源文件夹，用于整理分析所需的输入文件。将参考基因组序列以FASTA格式添加为ReferenceGenome.fa，并存入该文件夹中。将名为 ReferenceAnnotation.gtf 的基因模型注释文件添加到同一文件夹中。可选地，将 rRNA 基因注释作为名为 ReferenceRRNA.gtf 的 GTF 文件一并包含。将所有测序读段以压缩的 FASTQ 文件格式放入名为 Reads 的文件夹中。确保每个文件遵循命名格式。然后根据所使用的测序方法设置分析参数。为了评估序列质量，使用 Rsubread 软件包根据基因组 FASTA 文件构建参考基因组的索引。针对每个样本，使用 align 函数迭代并将测序读段比对至参考基因组。将生成的比对文件以 bam 格式存储在输出文件夹中。现在使用特征计数功能来统计比对到每个基因的读段。注释文件应为 GTF 格式。确保仅统计在基因组上具有唯一比对位置的读段。使用 RRNA 基因 GTF 文件，通过 feature counts 功能统计比对到 RRNA 基因的读段数量。允许将多重比对的读段包含在该统计中。获取 feature counts 功能为每个样本生成的读段分配统计结果。这些统计数据包括被分类为已分配、未比对、多重比对及其他类别的读段数量。单独收集核糖体 RNA 基因分配的统计数据。然后生成柱状图，以可视化展示前述步骤中的读段比对统计结果。根据分配给基因的读段数量对基因进行分组。将分类结果以柱状图形式展示。样本 S2R1 在修剪前后读段数量均较低。与样本 S2R2 的原始读段数量相比，其修剪后的读段数量明显减少，表明在修剪过程中去除了低质量的读段。比对分析发现了读段分配中的问题。样本 S2R3 出现了大量多重比对的读段以及较高比例的核糖体 RNA 读段。在样本 S2R4 中，有大量比例的测序读段未能比对到参考基因组，提示可能存在非目标生物序列的污染。样本 S2R1 至 S2R4 中，被分配到超过 100 条读段的基因数量较少。在相关性热图中，样本 S2R5 与样本 S1 的重复样本聚类在一起，而样本 S1R5 与样本 S2 的重复样本聚类在一起，表明可能存在重复样本标签混淆的错误。我们解决了RNA-Seq中缺失的质量控制问题，确保在下游基因表达分析之前进行可靠的数据评估。我们的工具整合了多种质量指标，为生物学家提供便捷、自动化且可重复的RNA-Seq评估。该工具能够探索RNA-Seq质量如何影响生物学解释，为透明且可重复的转录组学研究铺平道路。

rup-rna-seq-usability-assessment-pipeline-quality-control-for-bulk

Research

科研

JoVE Journal

Biology

生物学

Rup（RNA-seq 可用性评估管道）- 真核生物中批量 RNA-seq 实验的质量控制

我们介绍一种基于靶向RNA测序的方法，包括制备带索引的cDNA文库、使用定制探针进行杂交与捕获，以及数据分析，以检测特定转录本的基因表达、突变和基因融合。靶向RNA测序可在台式测序仪上对选定转录本进行经济高效且快速的评估。...

靶向RNA测序实验用于分析基因表达与基因组变异

本方案描述了染色质免疫沉淀测序（ChIP-seq）、4-硫代尿苷测序（4sU-seq）、总RNA测序（total...

转录因子结合、转录、翻译及降解的实时分析以展示细胞激活过程中的全局事件

本方案为研究人员提供了一种新工具，可用于监测多种模式生物中转录的保真度。...

真核生物中转录错误的全基因组监测

我们提供了一种偏倚低于标准方法且对2'-O-甲基化RNA具有更高灵敏度的详细小RNA文库构建方案。该方案可使用自制试剂以降低成本，也可使用现成试剂盒以提高操作便利性。...

改进小RNA测序：减少偏差并更好检测2'-O-甲基化RNA

本方法介绍了如何提高从福尔马林固定石蜡包埋（FFPE）RNA样本中获得的测序数据的质量和数量。我们描述了更准确评估FFPE-RNA样本质量的方法、测序文库的制备步骤，以及对FFPE-RNA样本数据进行分析的流程。...

降解FFPE-RNA样本测序及分析的优化

本实验方案的目的是利用RNA测序数据研究候选基因的进化与表达。...

利用RNA-seq研究分子进化与基因表达的生物信息学分析流程

基于转录组学的低成本药物筛选

绝对定量RNA测序（AQRNA-seq）是一种用于量化生物混合物中所有小RNA分布的技术。本文展示了AQRNA-seq的文库构建和数据处理步骤，用于量化在饥饿诱导的休眠期间Mycobacterium bovis...

AQRNA-seq用于小RNA定量

平面脂双层，也称为人工脂双层，可让您在明确定义的环境中研究离子导电通道。本文中，我们演示了制备脂双层室、电极所需的各个步骤，以及如何检测脂双层是否适用于单通道记录测量。...

用于电生理实验的人工双层膜制备

我们展示了利用DNA微阵列对神经系统进行表达谱分析的方法。本文描述了RNA质量控制、样品标记以及芯片杂交和扫描过程。...

DNA 微阵列：样品质量控制、芯片杂交与扫描

学习使用Rup管道进行真核生物批量RNA-seq质量控制，通过自动化工具评估测序数据可靠性。该方法帮助生物学家在基因表达分析前识别低质量读段和标记错误，确保转录组学研究的可重复性。

Rup（RNA-seq 可用性评估流程）— 真核生物批量 RNA-seq 实验的质量控制

探索更多视频