DNA对比图中的相似度,通常表示两段被拿来比较的序列有多接近,但它不是脱离比较条件的固定结论。同一个“98%”,可能代表一小段片段高度一致,也可能代表较大范围内的大部分碱基一致。读图时要把比较对象、相似度口径和覆盖范围放在一起看,再按用途判断哪个指标更重要。
相似度不等于“两个生物有多像”
DNA由碱基序列构成。比对程序会先把两条序列中可能对应的位置排列起来,再统计匹配、差异和缺口。图表上的相似度一般是对这些比对结果的概括,具体算法可能因软件、数据库或图表设置而不同。
在核苷酸序列比对中,常见指标是一致率,即比对区域里相同碱基所占的比例。比如100个对应位置中有96个相同,一致率就是96%。有些图表把它简称为“相似度”;另一些工具则可能采用更宽泛的相似性计算。因此,看到百分数时,应先看图例或说明采用什么口径。
这个数字也不能直接等同于亲缘关系、整体基因组相同程度或外形相似程度。若比较的是单个基因、一段保守区域或特定标记位点,结果只描述该比较范围,不能自动推及整条染色体或整个生物体。
一致率、覆盖度和差异数:三项指标各回答什么
| 指标 | 它回答的问题 | 读图时的重点 |
|---|---|---|
| 一致率 | 对齐位置中有多少碱基相同? | 确认统计的是碱基一致,还是其他相似性口径。 |
| 覆盖度 | 原始序列有多大比例参与了比对? | 看分母是查询序列、参考序列,还是两者中的某一方。 |
| 差异数 | 比对区域内有多少不匹配或缺口? | 留意替换、插入、缺失和无法判读的位置是否分开统计。 |
| 比对区域长度 | 百分数是建立在多长的序列上? | 短片段的高一致率,不一定代表广泛区域也相同。 |
覆盖度不是相似度。覆盖度描述“比上了多少”,一致率描述“比上的部分有多一致”。两者不能互相替代;只看一个百分数,容易把局部匹配误读成整体相同。
同样是高百分比,为什么结论会不同
下面是便于说明读法的示例数字,并非某个物种或检测项目的固定结果:
- 结果甲:一致率98%,覆盖度5%。图中只比对了目标序列的一小部分,这段区域很接近,但不足以说明整条序列普遍一致。
- 结果乙:一致率90%,覆盖度95%。匹配范围更广,整体差异较多;若要了解广泛序列差别,覆盖度高的信息更完整。
- 结果丙:一致率98%,覆盖度95%。大范围序列中一致比例较高,但仍需结合比对区域、样本质量和分析目标理解,不能只凭这一项推断个体关系。
这些例子说明,百分比必须和分母、长度一起读。若图上只标出“相似度”,没有说明比对区域和统计口径,这个数值表达的是有限的序列比较结果,而不是完整的生物学结论。
按比较目的选择该看的指标
| 比较目的 | 优先关注 | 适合的判断方式 |
|---|---|---|
| 观察两段序列是否接近 | 一致率、差异数、比对长度 | 先确认两段序列对应的是同一区域,再比较匹配与不匹配位置。 |
| 判断某序列是否覆盖充分 | 覆盖度、缺口和未比对部分 | 把覆盖范围与一致率并看,避免把局部高匹配当成全段结果。 |
| 比较不同物种的基因片段 | 同源区域、覆盖度、序列差异 | 优先比较具有对应关系的区域;不同基因的结果不能简单混为一谈。 |
| 分析个体差异或亲缘关系 | 多个标记位点及专门的分析结果 | 单个片段的相似度只能说明该片段,不能单独替代整体关系判断。 |
读DNA对比图时的选择顺序
- 先认比较对象:确认是两个人的样本、不同物种,还是两条数据库序列;对象不同,百分比所能回答的问题也不同。
- 再认统计口径:查看图例中的“identity”“similarity”等说明,分清精确碱基匹配率与更宽泛的相似性评分。
- 随后看覆盖范围:结合覆盖度、比对长度和缺口,判断百分数代表局部片段还是较大范围。
- 最后对应使用目的:看局部基因差异时侧重该区域的一致率;看序列是否完整匹配时同时关注覆盖度;分析个体或亲缘差异时,则需要符合该目的的数据和方法。
简而言之,DNA对比图中的相似度是“被比较序列在特定比对口径下有多接近”的数值。选读指标时,一致率看匹配程度,覆盖度看比较范围,差异数看具体变化;三者合起来,才比单独一个百分数更有解释力。