miRNA-03.测序分析与外泌体质控

以外泌体为对象的 miRNA 测序,在生物信息学阶段需同时回答两类问题:数据本身是否可靠(技术质控),以及样品是否符合外泌体 miRNA 的预期特征(生物学质控)。本文以外泌体质量分析为主轴组织 small RNA 测序(small RNA sequencing,sRNA-seq)关注点;差异表达与生物标志物挖掘仅在边界提及。

段末注释sRNA-seq 建库通常连接 3’/5’ 接头并扩增 ~18–40 nt 片段;与 mRNA-seq 的片段化策略与读长分布截然不同。


读前说明

  • 阈值(如 miRNA 占比 >50%)为经验范围,须结合物种、组织、分离试剂盒与测序平台在本实验室基线上校准,不宜跨项目硬套绝对数。
  • 配图位于 miRNA-03-测序分析与外泌体质控/
  • 前置阅读:miRNA-01-生物学背景miRNA-02-与外泌体的关系

1. 分析流程总览

图 1 外泌体 miRNA 测序质控分析流程示意

图 1 将质控拆为实验前(样本与分离)、测序后(读段层)、注释层miRNA 谱)与判定层(与参考签名对照)。生物信息学主要覆盖后三段;但若缺少实验元数据(分离方法、血样处理时间、是否溶血),计算指标难以解释。

1.1 典型计算流程

1
2
3
4
5
6
7
原始 FASTQ
→ 接头修剪 / 长度过滤(如 15–35 nt)
→ 去低质量 / 去重复(可选 UMI)
→ 比对(miRBase + 宿主基因组 + rRNA/tRNA 库)
→ 计数(成熟体 / isomiR 层级)
→ 归一化与 QC 报表
→ 与参考谱 / 阴性对照比较

常用工具举例:cutadapt/fastp(预处理)、bowtie/STAR(比对)、miRDeep2/featureCounts/custom(定量)、R/DESeq2(统计,差异分析场景)。


2. 读段层质控(FASTQ 级)

2.1 读长分布

健康 sRNA-seq 库在读长直方图上常在 21–23 nt 出现主峰,对应成熟 miRNA;另可见 ~30 ntpiRNA/tRF 肩峰或接头二聚体污染峰。

读长模式 可能含义
主峰 21–23 nt miRNA 富集良好
主峰 <18 nt** 或 **>30 nt 降解、接头残留、非 miRNA
宽峰、无主峰 复杂混合物或文库失败

外泌体馏分若丢失 miRNA 主峰,即使总读段数充足,也不宜视为「合格外泌体 RNA」。

2.2 质量分数与接头

  • Per-base quality:3’ 端质量下降常见,需合理修剪。
  • Adapter content:接头二聚体比例过高 → 有效 miRNA 读段不足,应回溯建库。
  • PCR duplicate rate:无 UMI(unique molecular identifier)时,高重复率可能提示起始 RNA 过少或扩增循环过多。

2.3 测序深度

质控场景下的深度目标与差异分析不同:

  • 批次 QC:每样 1–5 M 有效 miRNA 读段常足以评估组成比例;
  • 低丰度标志物:需更高深度与统计重复。

深度不足时,「未检出」不等于「不存在」,但高丰度 miRNA 的相对比例仍可用于纯度判断。


3. 比对与组成比例(核心 QC 指标)

3.1 miRNA 映射率

miRNA mapping rate 定义为比对至 miRBase 成熟体(及可选发夹)的读段占总有效读段的比例。外泌体 sRNA-seq 质控中最重要的单指标之一。

映射率(经验参考) 解读
>50–80%(细胞系外泌体、优质库) 通常认为 miRNA 富集良好
30–50% 需结合 tRNA/rRNA 检查;可能仍有分析价值
<30% 警惕污染、降解、错误建库或非 EV RNA 为主

人类血浆外泌体因背景 RNA 复杂,阈值可能低于细胞培养上清,应以本实验室阴性/阳性对照建立基线。

3.2 污染 RNA 读段

应单独统计比对至以下库的读段比例:

污染类型 参考库 / 标志 升高时的含义
rRNA 45S/28S/18S 细胞裂解、核糖体共沉淀
tRNA / tRF GtRNAdb 应激、凋亡、胞质泄漏
mRNA 片段 编码区外显子 凋亡小体、MV 混合、制备不纯
线粒体 RNA mt-rRNA/mt-tRNA 细胞破损
细菌 RNA 16S(样本或试剂污染) 环境/试剂污染

外泌体质控合格的典型期望:miRNA ↑,rRNA/tRNA ↓(相对血清游离 RNA 或裂解对照)。

3.3 溶血与血小板签名

血浆/血清样本应监控文献常用标志(人源示例):

miRNA 关联
miR-451amiR-16-5p 红细胞相关;溶血时异常升高
miR-223-3pmiR-191-5p 血小板 EV;血小板活化或制备残留

这些 miRNA 升高不一定否定外泌体存在,但提示血液处理血小板污染,需在 QC 报告中标注。


4. miRNA 谱特征与来源一致性

4.1 供体细胞系参考谱

对体外培养体系,将样品 miRNA 丰度谱与历史合格批次同步制备的母细胞裂解物对比:

  • Pearson/Spearman 相关:样品间或样品与参考的相关性;
  • 主成分分析(principal component analysis,PCA):技术重复应聚类;离群点提示批次或制备问题。

合格外泌体馏分应与供体细胞已知分泌谱方向一致,而与上清可溶性游离 RNA 区分。

4.2 外泌体富集 miRNA(文献常见)

无通用「外泌体专用 miRNA 金标准」,但以下在多种体系中反复报道(细胞类型依赖,作趋势参考而非绝对清单):

miRNA 备注
miR-21-5pmiR-126-3p 多种肿瘤与正常细胞分泌
miR-148a-3pmiR-101-3p EV 相关研究常见
let-7 家族 广泛存在,丰度高,适合作内参候选
miR-92a-3p 内皮来源讨论较多

质控逻辑:关键 miRNA 是否可检出 + 相对排名是否与参考一致,而非要求每个文献 miRNA 都显著。

4.3 阴性对照设计

对照类型 用途
PBS / 缓冲液空白 试剂与环境污染
EV-depleted 上清 区分囊泡包裹 vs 可溶性 RNA
母细胞裂解物 污染时谱趋向裂解物
FBS EV(培养体系) 血清来源 miRNA 背景

阴性对照出现与样品相近的 miRNA 映射率时,说明污染或串扰,样品结论需降级。


5. 定量细节:isomiR、5p/3p 与归一化

5.1 注释粒度

策略 质控场景建议
合并为成熟体(miRBase 标准名) 批次 QC、比例稳定,推荐默认
保留 isomiR 研究加工酶与修饰;QC 报告可附 top isomiR 比例
区分 5p/3p 若参考数据分臂报告,需一致

同一 miRNA3’ 异构体isomiR)若在某批次异常增多,可能提示 3’ 端加工或建库偏好,而非生物学真实变化。

5.2 归一化

质控比较(非差异表达)常用:

  • RPM(reads per million mapped reads):简单直观,适合组成比例比较;
  • TPM-样归一化:需明确分母(仅 miRNA 读段 vs 全部读段);
  • 内参 miRNA:选稳定、高丰度且非目标信号的 miRNA(如 miR-423-5p 在血清中作内参有争议,需验证)。

禁止在未校正文库大小与组成时,直接比较原始读段数判定「外泌体更纯」。

5.3 批次效应

分离日期、操作者、试剂批号、测序 lane 均可引入系统偏移。PCA层次聚类已知内参监控是批次 QC 的最低配置;必要时用 ComBat 等校正,但校正后不宜再用于绝对质控判定,仅用于探索性比较。


6. 与蛋白/物理指标的多组学对照

miRNA-seq 质控应与实验层指标三角验证

实验指标 miRNA-seq 的对应
NTA 粒径 / 浓度 颗粒数不足 → 总 RNA 低、文库失败
Western / 流式CD63 等) 蛋白阳性但 miRNA 差 → 可能空囊泡或 RNA 降解
电镜 形态支持 EV,不反映 RNA 内容
Bioanalyzer/TapeStationRNA 片段分布应与测序读长峰一致

单一指标合格不能替代 miRNA 组成 QC;反之 miRNA 谱良好但粒径异常,亦需回溯制备。


7. QC 报告建议字段

面向外泌体质量评估,生物信息学交付建议至少包含:

  1. 原始数据:总读段、修剪后读段、长度分布图
  2. 比对摘要miRNA%、rRNA%、tRNA%、未比对%
  3. Top 20 miRNA 条形图与参考批次对比
  4. PCA / 相关性热图(含对照)
  5. 溶血/血小板 相关 miRNA 读段数
  6. 重复一致性(若有技术/生物学重复)
  7. 结论标签:如「合格 / 边缘 / 不合格」及依据条目

判定规则应在项目 SOP 中预先定义,避免事后根据结果调阈值。


8. 常见失败模式速查

现象 可能原因 建议动作
miRNA% 极低,rRNA 细胞裂解、离心条件不当 重做分离;检查超速离心转速与时间
tRF 峰突出 应激、凋亡、营养耗尽 改善培养条件;排除死细胞
FBS 对照谱相似 血清 EV 未耗尽 使用认证 EV-depleted FBS
技术重复离散 操作不一致、RNA 量低 统一 SOP;提高起始 RNA
血浆 miR-451 极高 溶血 重新采血;溶血指数检测
文库 <18 nt 为主 严重降解或错误酶切 检查 RNA 保存与建库试剂

9. 与下游分析的边界

本文聚焦质控;若样品通过 QC,下游可衔接:

  • 差异 miRNADESeq2/edgeR,需生物学重复);
  • 靶基因富集(TargetScan + GO/KEGG);
  • 蛋白组/脂质组 多组学整合。

若样品未通过 QC,不建议直接进行疾病标志物结论;应优先修复实验链路。


小结

  • 外泌体 sRNA-seq 质控的核心是读长分布 + miRNA 组成比例 + 污染 RNA 监控 + 与参考谱/对照的一致性
  • 血浆样本需额外关注溶血血小板 miRNA;培养体系需关注 FBS凋亡
  • 阈值应建立在本实验室对照数据之上;生物信息学 QC 报告应与 NTA、标志蛋白等实验指标联合判读。
  • 本系列后续可扩展:具体流水线命令、参考签名库构建与自动化 QC 脚本(待补)。

段末注释RPM(reads per million)将某 miRNA 读段数除以该样总映射读段数再乘 $10^6$,便于样品间组成比较,但不校正 miRNA 文库大小以外的组成变化。

-------------本文结束感谢您的阅读-------------