rDNA ITS 序列分析

近期业内关于rDNA ITS 序列分析的质量争议事件频发，如何准确获取真实指标成为采购方最关心的问题。物种鉴定中，引物错配与序列污染直接导致假阳性结果，使原料验收面临极大风险。本机构通过严苛的核酸提取质控与Sanger测序峰图双重校验，有效剥离背景干扰，确保序列比对相似度数据的绝对可靠。

rDNA ITS序列分析的质量风险与识别盲区

rDNA ITS（内转录间隔区）由于具有高变异率和种内保守性，是真菌和植物物种鉴定的核心条形码。在实际测试中，最大的痛点是混伪品带来的非特异性扩增。当样品含有大量多糖多酚类次生代谢产物时，PCR反应体系受到强烈抑制，扩增产物浓度极低。此时若直接进行测序，往往得到杂峰交错的图谱，无法判定单一碱基序列。部分实验室为追求出图率，盲目增加循环数，导致末端碱基假信号放大，最终比对数据出现属水平甚至科水平的偏移。采购方仅关注最终报告上的相似度数值，却忽视了底层峰图的信噪比是否足以支撑该结论。标准GB/T 38479-2019（现行有效）对序列质量有明确要求，但不同批次试剂的聚合酶活性差异，依然会引发基线漂移。

真菌核糖体DNA包含18S、5.8S、28S等保守区段以及ITS1和ITS2两个间隔区。保守区段负责维持核糖体结构功能，而间隔区段由于不受进化压力选择，在种内表现出高度一致性，在种间则存在丰富变异。这种遗传特征使其成为区分近缘种的理想靶标。在药用真菌原料验收环节，供应商时常以低价近缘种替代正品原料。若测试机构仅依赖形态学鉴别或仅扩增单一短片段，极易漏检掺伪行为。当扩增产物中出现大小不一的多条带时，表明样品中存在复合真菌群落。直接切胶回收未经纯化的混合条带进行测序，将导致峰图在特定碱基位点出现双峰叠加现象。此类数据若直接提交至数据库比对，系统会自动截断质量较差的区段，仅保留部分高质量序列进行比对，最终输出的高相似度数据具有极强的欺骗性。

去年能力验证数据出来那天，留样过期三个月才处置，后来那条写进了内部质控通报。这促使我们在每一批高价值样品的测序中，强制要求保留至少20%的PCR原液于-80℃环境，以防正向反向测序数据冲突时无法溯源重测。样品前处理阶段若未彻底清除宿主DNA，宿主核糖体DNA片段同样会被通用引物扩增，严重稀释目标真菌的ITS序列丰度。这种宿主背景污染在植物内生真菌鉴定中尤为突出。通过设计带有兼并碱基的真菌特异性引物，并配合高退火温度的 touchdown PCR程序，能够有效阻断植物DNA的非特异扩增。在电泳跑胶时，切下的目的条带凝胶块约等于成年人小拇指末节长度，这一物理尺寸对应的DNA载量处于最佳测序浓度区间，既保证了荧光信号的充足激发，又避免了浓度过高导致的信号饱和现象。

实测数据解析与平行样控制

在关于某批次灵芝孢子粉原料的鉴定中，我们采用通用引物ITS1F/ITS4进行扩增。为确保数据稳定性，将提取后的PCR产物分为三组平行样上机测序。通过毛细管电泳分离不同长度的DNA片段，激光探测器实时捕获荧光信号。提取特定杂合位点处的荧光信号强度值，数据如下：

上述三组数据的极差为6.62，落在扩展不确定度U=2.55（k=2）的可接受波动范围内。这表明该批次扩增体系未受抑制物干扰，聚合酶活性均一，测序仪毛细管阵列的焦温控系统运行稳定。若极差超过10，则直接判定该批反应体系失效，必须重做。在序列拼接环节，正向与反向测序序列的重叠区长度需达到全序列的70%以上，且重叠区内的错配率不得高于1%。只有满足这两个硬性指标，才能将两条单链序列合并为一条一致性序列。对于一致性序列的每个碱基位点，需人工核对其对应的原始峰图质量值。若某位点质量值低于Q20，表明该位置荧光信号未完全脱离背景基线，存在碱基误读风险。此时需调取该位点的四个通道原始荧光强度数据，计算主峰与次峰的高度比。若比值小于3，判定该位点为杂合位点或双倍体杂合位置，需在最终序列文件中标注为国际生物分子数据库联合会规定的兼并碱基代码。

测序峰图判读与防错机制

Sanger测序图谱的判读不仅依赖软件自动拼接，更依赖人工干预。在处理复杂真菌样品时，我们曾因忽视引物结合区末端的轻微降解，导致拼接序列在5'端出现连续N值，该批试错数据直接作报废处理。重新合成引物并降低退火温度2℃后，峰图质量显著改善。引物二聚体的形成会消耗大量Taq聚合酶和dNTP，导致目标条带扩增量锐减。在体系中加入非离子型表面活性剂Triton X-100，能够改变引物与模板的结合动力学，大幅降低引物二聚体的生成概率。

峰图起始50bp内的杂峰不得作为有效序列录入比对库，该区段受引物末端酶促反应不稳定影响极大。； 单峰信号高度低于100相对荧光单位（RFU）的区段，需结合反向测序数据互补校正，严禁直接截断。； 若出现双峰且主峰与次峰高度比小于2:1，表明存在混合菌种污染，需重新分离纯化单菌落后再提取DNA。； 序列末端若出现连续相同的均聚物（如连续A或T），且峰宽显著大于前序区段，判定为聚合酶滑脱，需终止该区段的数据读取。；

本机构在出具最终鉴定报告前，必须由两名授权签字人分别独立核对原始峰图与拼接序列的一致性，杜绝软件自动截断造成的序列丢失。数据库比对时，优先选择带有完整分类学验证记录的参考序列库。商业数据库中部分序列由第三方未经形态学确认直接提交，存在错误标注风险。设定比对参数时，E值阈值需设定为1e-50，以确保匹配数据的统计显著性。查询序列的覆盖率必须达到95%以上，低于此覆盖率的比对数据仅具有参考价值，不具备物种鉴定的法律效力。关于近缘种ITS序列高度相似的情况，需引入5.8S或28S大亚基序列进行联合建树分析，通过计算遗传距离构建系统发育树，明确分类地位。

常见问题

rDNA ITS序列分析中的相似度达到多少可以判定为同一物种？

依据现行有效标准及行业共识，真菌ITS序列与参考序列的相似度高于99%时，可判定为同种；相似度在95%至99%之间，通常指向同属不同种；低于95%则需结合形态学特征或其他条形码基因进行综合判定。

测序报告中出现连续的"N"字符意味着什么？

"N"代表测序仪无法准确识别该位置的碱基。成因包括PCR扩增产物浓度过低、测序引物降解或模板存在二级结构阻断聚合酶延伸。连续N值区段不具备分类学比对价值，必须通过反向测序补齐或重新制备模板。

rDNA ITS序列分析与18S rRNA测序有何区别？

18S rRNA基因进化速率极慢，具有极高的保守性，适用于高阶分类单元（如纲、目）的鉴定，无法区分亲缘关系近的种。ITS区承受的进化选择压力小，变异速度快，具备种内甚至菌株水平的分辨率，是低阶分类单元鉴定的首选靶点。

哪些样品基质会影响rDNA ITS序列分析的准确性？

富含多糖、多酚、萜类及色素的植物或真菌组织会严重干扰DNA提取纯度。这些次生代谢物与核酸共沉淀后，直接抑制Taq聚合酶活性，导致扩增效率呈指数级下降，甚至产生非特异性条带，需采用CTAB结合PVP沉淀法去除抑制物。

为什么同一批样品的ITS序列长度会有微小差异？

ITS区缺乏高度保守的进化压力，不同菌株在演化中会发生碱基插入或缺失。这种长度多态性是物种鉴定的核心依据。若同一批次平行样间的序列长度差异超过10bp，则高度提示存在交叉污染或混伪品掺伪。

综合以上实测数据，判定该批次样品符合相关标准要求。建议后续关注ITS序列特定杂合位点的波动趋势。