中华人民共和国农业行业标准
NY/T 4912—2025
小麦品种及其实质性派生品种鉴定
MNP标记法
Identification of wheat varieties and their essentially derived varieties—
MNP marker method
2025-12-09 发布
2026-05-01 实施
中华人民共和国农业农村部 发 布
前 言
本文件按照 GB/T 1 .1—2020« 标准化工作导则 第 1 部分 : 标准化文件的结构和起草规则»的规定起草 .
请注意本文件的某些内容可能涉及专利 . 本文件的发布机构不承担识别专利的责任 .
本文件由农业农村部种业管理司提出 .
本文件由全国植物新品种测试标准化技术委员会(SAC/TC 277)归 口 .
本文件起草单位 :江汉大学 、农业农村部科技发展中心 、武汉明了生物科技有限公司 、山东省农业科学院作物研究所 、山东省农业科学院 、湖北省农业科学院粮食作物研究所 、深圳华大智造科技股份有限公司 .
本文件主要起草人 :李甜甜 、陈红 、张秀杰 、韩瑞玺 、彭海 、方治伟 、周俊飞 、张晗 、郑永胜 、李论 、张全芳 、高利芬 、陈利红 、肖华锋 、张静 、董静 、刘易科 、朱展望 、陈雪燕 、余进文 .
小麦品种及其实质性派生品种鉴定 MNP标记法
1 范围
本文件规定 了 利 用 多 核 苷 酸 多 态 性 (Multiple Nucleotide Polymorphism , MNP)标 记 法 进 行 小 麦(Triticum aestivum L.) 品种及其实质性派生品种鉴定的术语和定义 、原理 、试剂和材料 、仪器设备 、操作程序 、质量控制 、数据分析 、结果判定与表述 .
本文件适用于小麦品种及其实质性派生品种的鉴定 .
2 规范性引用文件
下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款 . 其中 ,注日期的引用文件 ,仅该日期对应的版 本 适 用 于 本 文 件 ; 不 注 日 期 的 引 用 文 件 , 其 最 新 版 本(包 括 所 有 的 修 改 单)适 用 于 本文件 .
GB/T 3543 .2 农作物种子检测规程 第 2 部分 :扦样
GB/T 6682 分析实验室用水规格和试验方法
GB/T 38551 植物品种鉴定 MNP标记法
3 术语和定义
下列术语和定义适用于本文件 .
3 .1
多核苷酸多态性 multiplenucleotidepolymorphism ,MNP
在一段核苷酸序列中 , 由一个或多个核苷酸变异引起的序列多态性 .
[来源 :GB/T 38551—2020 ,3 .1 ,有修改]
3 .2
实质性派生品种 essentiallyderived variety,EDV
由原始品种实质性派生 ,或者由该原始品种的实质性派生品种派生出来的品种 ,与原始品种有明显区别 ,并且除派生引起的性状差异外 ,在表达由原始品种基因型或者基因型组合产生的基本性状方面与原始品种相同 .
3 .3
平均覆盖倍数 averagecoverage
在基因组上比对到所有标记位点上的所有测序片段的总数与标记位点总数的比值 .
3 .4
检出标记位点 detected markers
至少有一个等位基因型且有 20 条及以上测序片段支持的标记位点 .
4 原理
小麦品种的基因组中存在着能够世代稳定遗传的 MNP位点 ,不同小麦品种间同一组 MNP位点的序列上可能存在差异 ,这种差异可通过基于 PCR、二代高通量测序以及生物信息学方法进行检测 ,进而区分不同的品种 .
5 试剂和材料
除非另有规定 ,仅使用分析纯试剂 ,实验用水符合 GB/T 6682 中规定的一级水的要求 .
5 .1 多重 PCR扩增与文库构建试剂盒
应匹配 MNP标记和标记检测引物 , 以及高通量测序试剂盒 .
5 .2 高通量测序试剂盒
应匹配高通量测序仪 .
5 .3 MNP标记和标记检测引物
应符合附录 A 的要求 .
6 仪器设备
6 .1 离心机
最大转速度不小于 12 000 r/min.
6 .2 电泳仪
6 .3 PCR扩增仪
6 .4 高通量测序仪
测序读长不低于 300 bp.
6 .5 计算机服务器
7 操作程序
7 .1 样品准备
送检样品为种子 、幼苗 、叶片等组织或器官 . 试验样品抽取的样本数量宜为 30 个以上 . 样品需要扦样时 ,应符合 GB/T 3543 .2 的规定 . 样品可以混合检测 ,必要时进行个体检测 .
7 .2 DNA提取
7 .2 .1 CTAB法
将混合样品在液氮中研磨至粉末 ,迅速取约 200 mg粉末置于 2 .0 mL离心管中 ;每管加入 600μL 65 ℃预热的 CTAB提取液 ,充分混合 ,65 ℃水浴 45 min~ 60 min,每间隔 15 min轻缓颠倒混匀 ,水浴后 12 000 r/min离心 10 min. 吸取上清液转移至新的离心管中 , 每管加入与提取液等体积的三氯甲烷和异戊醇 ,轻缓混匀后静置 10 min;12 000 r/min离心 10 min后吸取上清液至新的离心管中 ,加入等体积预冷的
,液,0, i,,iin溶. 上D,
度和1量:以,0推.取 ,其.他 DNA提取方法均适用于本文件 . DNA
溶液的紫外吸光度 OD260 与 OD280 的比值宜介于 1 .7 ~ 2 .0 .
注 2 :三氯甲烷和异戊醇混合液中三氯甲烷与异戊醇的体积比为 24 ∶ 1 .
7 .2 .2 试剂盒提取法(仲裁法)
利用核酸提取试剂盒提取待测样品 DNA. 提取的 DNA 在 260 nm 与 230 nm 处的吸光度值的比值大于 2 .0 ,260 nm 与 280 nm 处的吸光度比值介于 1 . 7 ~ 1 .9 . 方法见附录 C 中的 C.2 .
7 .3 多重 PCR扩增与文库构建
按多重 PCR扩增与文库构建试剂盒的说明书进行 DNA 质控 、多重 PCR 扩增 、文库构建与纯化 . 其中 ,多重 PCR 的扩增循环数不高于 20 次 . 方法见 C.3 .
7 .4 高通量测序
按高通量测序试剂盒和高通量测序仪的操作说明 ,对 7 .3 中获得的高通量测序文库进行高通量测序 .高通量测序时标记位点的平 均 覆 盖 倍 数 宜 设 置 为 700 倍 以 上 , 测 序 片 段 总 读 长 不 小 于 300 bp. 方 法 见C.4 .
8 质量控制
8 .1 环境
样品准备 、DNA提取 、多重 PCR扩增 、文库构建和高通量测序宜在规定的区域或相互隔离的区域按单一方向进行操作 ,不同区域的仪器设备需专用。
8 .2 测序数据
8 .2 .1 高通量测序原始数据质量应满足所采用的高通量测序仪的操作手册中所规定的测序质量要求。
8 .2 .2 将样品的测序数据比对到参考基因组的标记位点上 ,统计第一次检测的标记位点的平均覆盖倍数C1 。
8 .2 .3 当 C1 小于 500 时 ,判定样品的测序数据量不足 ,从 7 .4 或之前的步骤开始重新实验至检测标记位点的平均覆盖倍数 C1 大于或等于 500 。
8 .2 .4 当 C1 大于或等于 500 时 ,进一步计算检出标记位点的比例 R1 。
R1 按公式(1)计算。
R1 = × 100 … … … … … … … … … … … … … … … … … (1)
式中 :
R1 — 样品检出标记位点的比例 ,单位为百分号(%) ;
T1 — 样品检出标记位点的数目 ;
T — 样品检测标记位点的数目。
8 .2 .5 当 R1 大于或等于 95%时 ,判定测序数据合格 ;否则 ,从 7 .2 或之前的步骤重新实验至第二次检出的标记位点的平均覆盖倍数 C2 大于或等于 500 。
8 .2 .6 当 C2 大于或等于 500 时 ,进一步计算第一次和第二次共同检出的标记位点的比例 R2 。
R2 按公式(2)计算。
R2 = × 100 … … … … … … … … … … … … … … … (2)
式中 :
R2 — 第一次和第二次共同检出的标记位点的比例 ,单位为百分号(%) ;
T12 — 第一次和第二次共同检出标记位点的数目 ;
T11 — 第一次检出标记位点的数目 ;
T2 — 第二次检出标记位点的数目。
8 .2 .7 当 R2 大于或等于 95%时 ,判定测序数据合格。
注 : 附录 C提供了一个操作案例 ,其他满足 DNA提取 、文库制备和高通量测序的试剂 、仪器设备均可。
9 数据分析
9 .1 测序数据比对与记录
9 .1 .1 将测序数据比对到参考基因组上的每个 MNP标记位点上。
9 .1 .2 位点的基因型指该位点的所有检出等位基因型 ,其中 ,检出等位基因型指从该标记第一个到最后一个碱基构成的检出 DNA 片段 ,不同检出等位基因型用 “/”隔开。 位点的基因型记录实例见 C8 。
9 .2 遗传相似度计算
遗传相似度按公式(3)计算。
GS= × 100 … … … … … … … … … … … … … … (3)式中 :
GS — 待测品种与对照品种的遗传相似度 ;
nij — 待测品种与对照品种中均检出的但基因型无任何差异的标记位点的数目 ; Nij — 待测品种与对照品种中均检出标记位点的数目。
10 结果判定与表述
10 .1 品种鉴定
10 .1 .1 当待测品种与对照品种的 GS小于 96 .00%时 ,判定为“不同品种 ”.
10 .1 .2 当待测品种与对照品种的 GS大于或等于 96 .00%时 ,判定为“疑同品种 ”.
10 .2 实质性派生品种判定
10 .2 .1 当待测品种与对照品种的 GS小于 90 .00%时 ,判定为“非实质性派生品种 ”.
10 .2 .2 当待测品种与对照品种的 GS大于或等于 90 .00%时 ,判定为“疑似实质性派生品种 ”.
10 .3 结果表述
待测品种 与对照品种 比较位点数为 ,差异位点数为 ,遗传相似度为 ,判定为 .
示例 1 :待测品种 A与对照品种 B 比较位点数为 780 ,差异位点数为 5 ,遗传相似度为 99 .36% ,判定为疑同品种 .
示例 2 :待测品种 A与对照品种 B 比较位点数为 778 ,差异位点数为 150 ,遗传相似度为 80 .72% ,判定为非实质性派生品种 .
注 :结果表述也可采用表格形式 ,但需完整包含上述内容 .
附 录 A
(规范性)
MNP标记和标记检测引物
MNP标记编号 、染色体位置 、引物序列 、变异类型 、参照品种等信息见表 A.1 。
表 A.1 MNP标记和标记检测引物信息
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
表 A.1 (续)
附 录 B (资料性)
DNA提取溶液配制
试剂配制用水应符合 GB/T 6682 的要求。
B.1 0.5 mol/L EDTA 溶液
称取 186 .1 g Na2 EDTA . 2H2 O 溶于 800 mL水中 ,充分搅拌溶解 ,加 NaOH 调 pH 至 8 .0 ,加水定容至 1 000 mL,121 ℃高压灭菌 20 min。
B.2 0.5 mol/L HCl溶液
量取 25 mL浓盐酸(质量分数为 36% ~ 38%) ,加水定容至 500 mL。
B.3 1 mol/L NaOH 溶液
称取 40 .0 g NaOH 溶于 800 mL水中 ,充分搅拌溶解 ,加水定容至 1 000 mL。
B.4 1 mol/L Tris-HCl溶液
称 取 60 .55 g三羟甲基氨基甲烷溶于 400 mL水中 ,加 HCl调 pH 至 8 .0 ,加水定容至 500 mL,121 ℃高压灭菌 20 min。
B.5 CTAB提取液
称取 20 .0 g CTAB、81 .7 g NaCl置 于 1 000 mL 烧 杯 中 , 量 取 100 mL 1 mol/L Tris_HCl溶 液 和40 mL 0 .5 mol/LEDTA溶液倒入烧杯中 ,加 700 mL水 ,充分搅拌溶解 ,加水定容至 1 000 mL,121 ℃高压灭菌 20 min。
B.6 TE 缓冲液
量取 5 mL 1 mol/LTris-HCl溶液和 1 mL 0 .5 mol/LEDTA溶液 ,加水定容至 500 mL,121 ℃高
附 录 C
(资料性)
品种鉴定流程示例
C.1 样品准备
小麦样品来源于农业农村部科技发展中心 . 待测样品名称 ‘郑麦 101 ’, 实验编号 XMI220424200;对照样品名称‘徐麦 31 ’,实验编号 XMI220424070 .
从待测样品与对照样品中分别抽取 30 粒以上的种子 ,发芽至幼苗长度为 5 cm~ 10 cm. 每棵幼苗的叶片取约 10 mm 圆片 ,装入离心管中 .
C.2 DNA提取
采用某公司生产的新型植物基因组 DNA提取试剂盒提取小麦叶片 DNA,具体步骤如下 :
a) 向装 有 叶 片 的 离 心 管 中 加 入 液 氮 充 分 碾 磨 后 , 加 入 400 μL缓 冲 液 LP1 和 6 μL RNase A (10 mg/mL) ,涡旋振荡 1 min,室温放置 10 min.
b) 加入 130 μL缓冲液 LP2 ,充分混匀 ,涡旋振荡 1 min.
c) 12 000 r/min离心 5 min,将上清液移至新的离心管中 .
d) 加入 1 .5 倍体积的缓冲液 LP3 (使用前 ,请先检查是否已加入无水乙醇) ,立即充分振荡混匀 15 s,此时可能会出现絮状沉淀 .
e) 将上一步所得溶液和絮状沉淀加入一个吸附柱 CB3 中 (吸附柱放入收集管中) , 12 000 r/min离心 30 s,倒掉废液 , 吸附柱 CB3 放入收集管中 .
f) 向吸附柱 CB3 中加入 600 μL漂洗液 PW(使用前 ,请先检查是否已加入无水乙醇) ,12 000 r/min离心 30 s,倒掉废液 ,将吸附柱 CB3 放入收集管中 .
g) 重复操作步骤 f) .
h) 将吸附柱 CB3 放回收集管中 ,12 000 r/min离心 2 min,倒掉废液 . 将吸附柱 CB3 室温放置数分钟 ,彻底晾干吸附柱中残余的漂洗液 .
i) 将吸附柱 CB3 转入一个干净的离心管中 , 向吸附膜的中间部分悬空滴加 50 μL洗脱缓冲液 TE,
温,Ni .min,12 000 r/min离心 2 min,将溶液收集到离心管中 . 产物放于- 20℃冰
j) DNA 质 检 . 用 分 光 光 度 计 测 定 并 计 算 XMI220424200 和 XMI220424070 的 DNA 溶 液 在 260 nm 与 280 nm 处的吸光度比值 ,分别为 1 .86 和 1 .95 . 在 260 nm 与 230 nm 处的吸光度比值分别为 2 .04 和 2 .09;取 4 μL DNA在 1%的琼脂糖凝胶上电泳 ,检测 DNA 条带是否完整 ,结果如图 C.1 所示 ;取 1 μL DNA用 Qubit荧光定量仪测定 XMI220424200 和 XMI220424070 的 DNA浓度 ,分别为 30 .4 ng/μL和 39 .6 ng/μL.
注 :最左边条带为某公司生产的分子 量 标 准 ,其 由 上 至 下 片 段 大 小 分 别 为 23 130 bp、9 416 bp、6 557 bp、4 361 bp、2 322 bp.
C.3 多重 PCR扩增与文库构建
C.3 多重 PCR扩增
采用某公司生产的多重扩增试剂盒进行多重 PCR扩增与文库构建 ,该试剂盒匹配后续步骤中某公司的测序平台 ,需自备 80%乙醇 ,其余试剂均为试剂盒提供 .
图 C.1 XMI220424200 和 XMI220424070 的 DNA 电泳图
由于 XMI220424200 与 XMI220424070 的多重 PCR扩增与文库构建实验流程完全一样 , 因此 , 除特别说明 ,下面仅就 XMI220424200 的实验流程进行说明。
a) 配制多重 PCR扩增体系。 在 PCR 管中加入 4 μL表 A.1 中的多重 PCR 引物混合物(每条引物浓度为 0 .2 μmol/L) 、4 μL样品 XMI220424200 的基因组 DNA、10 μLGenoPlexs 3 × T Master Mix和 12 μL ddH2O,振荡混匀。
b) 多重 PCR扩增反应。 多重 PCR扩增程序 :95 ℃ ,3 min; (95 ℃ ,20 s,60 ℃ ,4 min) × 15 个循环 ;
C.3 .2 多重 PCR扩增产物纯化
a) 向 C.3 .1 中获得的多重 PCR扩增产物中加入 12 μL (C.3 .1 中获得的 30 μL多重 PCR 扩增产物的 0 .4倍体积) GenoPrep DNA Clean Beads,振荡混匀后 ,室温静置 5 min。
b) 将 PCR管置于磁力架上吸附磁珠 ,直至溶液澄清。
c) 用移液器吸取上清液至新的 1 .5 mL离心管中 ,避免吸到磁珠。
d) 向上清液中加入 18μL (C.3 .1 中获得的 30 μL多重 PCR 扩增产物的 0 .6 倍体积)的 GenoPrep DNA Clean Beads,振荡混匀后 ,室温静置 5 min。
e) 用磁力架吸附磁珠 ,直至溶液澄清。 用移液器小心吸取上清液 ,弃上清液 , 留磁珠。
澄清。 用移液器小心吸取上清液 ,弃上清液 , 留磁珠。
g) 加入 100 μL80%乙醇(现用现配) ,用移液器小心去除上清液 ,避免吸到磁珠。
h) 室温放置 ,直至乙醇挥发干净 ,避免过干 ,获得纯化的多重 PCR扩增产物。
C.3 .3 高通量测序文库构建
a) 向 C.3 .2 中获得的纯化的多重 PCR扩增产物中 ,加入 10 μLGenoPlexs3 ×T MasterMix、2 μL浓度为 5 μmol/L的 P5 primer、2 μL浓度为 5 μmol/L的 P7 barcodeprimer(引物中包含样品条形码 ,样品 XMI220424200 的条形码序列为 GGTTGTCTAG)和 16 μL ddH2 O。
b) 将配制好的反应体系振荡混匀并短暂离心 ,按如下程序进行 PCR反应 :95 ℃ ,3 min; (95 ℃ , 15 s;58 ℃ ,15 s;70 ℃ ,30 s) × 8 个循环 ;72 ℃ ,5min。
c) 反应结束后 , 即构建好 30 μL样品 XMI220424200 的高通量测序文库。
C.3 .4 高通量测序文库纯化
a) 向 C.3 .3 中获得的高通量测序文库中加 24 μL (C.3 .3 中获得的 30 μL高通量测序文库的 0 .8倍体积) GenoPrep DNA Clean Beads,振荡混匀 ,室温静置 5 min。
b) 用磁力架吸附磁珠 ,直至溶液澄清。 用移液器小心吸取上清液 ,弃上清 , 留磁珠。
c) 加入 40 μLGenoPlexsBW07 Buffer,涡旋均匀。
d) 用磁力架吸附磁珠 ,直至溶液澄清。 用移液器小心去除上清液 ,避免吸到磁珠。
f) 加入 35 μL 10 mmol/LTris_HCl(pH = 8 .0) ,充分悬浮磁珠 ,室温静置 5 min。将离心管置于
获得的纯化高通量测序文库直接用于后续实验或置于- 20 ℃ 保存。
图 C.2 样品 XMI220424200 和样品 XMI220424070 的高通量测序文库电泳图
注 :最左边条带为某公司生产的分子量标准 ,其由下至上分子量大小分别为 100 bp、200 bp、300 bp、400 bp、500 bp和600 bp。
h) 高通量测序文库混合。 分别取样品 XMI220424200 和样品 XMI220424070 各 100 ng,振荡混匀
并离心 ,获得上机混合测序文库。
C.4 高通量测序
C.4 .1 上机混合测序文库变性
对 C.3 中获得的上机混合测序文库进行变性。
a) 采用 Qubit荧光定量仪测定上机混合测序文库的浓度为 38 .2 ng/μL。该试剂盒推荐文库投入
文库的质量 X 660ng ,其中 ,DNA 主片段长度按 400 bp计算 , 1
b) l的,为B,C.增。仪上 95 ℃变性
3 min,获得文库变性产物。
C.4 .2 文库变性产物单链环化
按某公司生产的环化试剂盒的操作手册对文库进行环化。
a) 在冰上配制 12 .1 μL的单链环化反应液 ,其中包含 11 .6 μL的 SplintBuffer和 0 .5 μL的 DNA
Rapid Ligase。
b) 向上述单链环化反应液中加入 C.4 .1 中获得的文库变性产物 ,涡旋振荡 6 次 ,每次 3 s,瞬时离心将反应液收集至管底后 ,37 ℃保温 30 min,获得单链环化产物。
C.4 .3 酶切消化
对环化产物进行酶切 ,执行以下步骤。
a) 在冰上配制 4 .0 μL的酶切消化反应液 ,其中包含 1 .4 μL 的 Digestion Buffer和 2 .6 μL 的 Di_ gestion Enzyme。
b) 向酶切消化反应液中加入 C.4 .2 中获得的单链环化产物后 ,涡旋振荡 6 次 ,每次 3 s,瞬时离心将反应液收集至管底 ,37℃反应 30 min。
c) 反应结束后 ,加入 7 .5 μLDigestion StopBuffer,涡旋振荡 6 次 ,每次 3 s,瞬时离心将反应液收集至管底 ,获得酶切消化后的单链环化产物。
d) 将酶切消化后的单链环化产物转移到新的 1 .5 mL离心管中。
C.4 .4 单链环化产物纯化
对单链环化产物进行纯化 ,执行以下步骤。
a) 提前 30 min取出 DNA Clean Beads置于室温 ,使用前充分振荡混匀。
b) 吸取 170 μL DNA Clean Beads至 C.4 .3 获得的酶切消化后的单链环化产物中 ,用移液器轻轻吹
打至少 10 次至完全混匀 ,最后一次应确保将吸头中所有液体及磁珠打入 1 .5 mL离心管中。
c) 室温孵育 10 min。
d) 瞬时离心 ,将 1 .5 mL离心管置于磁力架 ,静置 2 min~ 5 min至液体澄清 ,用移液器小心吸取并丢弃上清液。
e) 保持 1 .5 mL离心管置于磁力架上 ,加入 500 μL新鲜配制的 80%乙醇漂洗磁珠及管壁 ,小心吸取并弃上清液。
f) 重复步骤 e) ,尽量吸干管内液体。
g) 保持 1 .5 mL离心管置于磁力架上 ,打开 1 .5 mL离心管管盖 ,室温干燥 ,直至磁珠表面无反光 、无开裂。
h) 将 1 .5 mL离心管从磁力架上取下 ,加入 22 μLTE Buffer进行 DNA洗脱 ,用移液器轻轻吹打至少 10 次至完全混匀。
i) 室温下溶解 10 min。
j) 瞬时离心 ,将 1 .5 mL EP管置于磁力架上 ,静置 2 min~ 5 min至液体澄清 ,将 20 μL上清液转移
k) ., 化its。Dsit量。
浓度(D)为 2 .3 ng/μL,按公式(C.1)对环化产物质量浓度进行换算 ,计算环化产物的摩尔浓度
(C)为 17 .4 fmol/μL。
式中 :
C — 环化产物的摩尔浓度 ,单位为飞摩尔每微升( fmol/μL) ;
D — 环化产物质量浓度 ,单位为 ng/μL ;
3030/400 — 本公式中的换算系数。
C.4 .5 高通量测序
采用某公司生产的高通量测序试剂盒 ,对环化产物进行高通量测序。
a) 从某公司生产的高通量测序试剂盒中取出 DNB制备缓冲液 、DNB聚合酶混合液 Ⅰ 、TE缓冲液和 DNB终止缓冲液 ,置于冰盒上约 0 .5 h,待试剂融化后 ,使用涡旋振荡器振荡混匀 5 s后 ,短暂离心置于冰盒上备用。
b) 取 0 .2 mL PCR管 ,在冰上配制 40 μL的反应混合液体系 ,其中包括 2 .3 μL的环化文库 ssDNA
( 环 化文库 ssDNAL ) 、17 .7 μL的 TE缓冲液和 20 μL的 DNB 制备缓冲液 .
c) 将 (b)中获得的混合液体系涡旋振荡器振荡混匀 ,离心 5 s,置于 PCR仪中 95 ℃ , 1 min; 65 ℃ , 1 min;40 ℃ ,1 min和 4 ℃ ,10 min.
d) 取出 DNB聚合酶混合液 Ⅱ (LC)置于冰盒上 ,短暂离心 5 s,置于冰盒上备用[请勿将 DNB聚合酶混合液 Ⅱ (LC)置于室温 ,请勿长时间触碰管壁] .
e) 当 c) 中 PCR仪达到 4 ℃后取出 PCR管 ,离心 5 s;在冰上加入 40 μLDNB聚合酶混合液 Ⅰ 和 4 μLDNB聚合酶混合液 Ⅱ (LC) ,振荡混匀 ,离心 5 s,立即置于 PCR仪中 ,在 35℃热盖条件下 30 ℃ ,25 min和 4 ℃ ,10 min;PCR仪降温至 4 ℃后 ,立即加入 20 μLDNB终止缓冲液 ,用阔口吸头缓慢地吹打混匀 5 次 ~ 8 次 (切勿振荡或剧烈吹打) ,获得制备好的 DNA纳米球(DNB) .
f) DNB浓度测定 . 采用 Qubit® ssDNA AssayKit在 Qubit荧光定量仪上检测制备好的 DNB 的浓度为 19 .2 ng/μL. 制备好的 DNB可置于 4 ℃保存 48 h.
g) DNB加载 . 取出 DNB加载缓冲液 Ⅰ 和 DNB加载缓冲液 Ⅱ (如发现 DNB加载缓冲液 Ⅱ 中有结晶 ,使用涡旋振荡器持续剧烈振荡 1 min~ 2 min至沉淀重新溶解 ,短暂离心后方可使用) ,置于冰盒上融化后 ,漩涡震荡 5 s混匀 ,短暂离心后置于冰盒上备用 .
h) 在 0 .5 mL冻存管中配制 DNB加载体系 :50 μL的 DNB加载缓冲液 Ⅰ 、50 μL的 DNB加载缓冲液 Ⅱ 、1 μL的 DNB聚合酶混合液 II(LC)和 100 μL(浓度大于 8 ng/μL)制备好的 DNA 纳米球(DNB) ;用阔口吸头缓慢混匀 5 次 ~ 8 次 (切勿离心 、振荡及剧烈吹打) ,4 ℃保存备用 .
i) 准备测序试剂槽 . 取出测序试剂槽 ,常温水浴解冻 3 h ~ 4 h(或者提前 1 d将其置于 2 ℃ ~ 8 ℃冰箱解冻)后 ,置于 2 ℃ ~ 8 ℃冰箱备用 ;使用前颠倒混匀试剂槽 3 次 ,然后将试剂槽置于正前方 ,前后左右剧烈晃动 10 次 ~ 20 次 ,直至试剂中无肉眼可见的分层 ,尤其是该测序试剂槽的 17号试剂和 18 号试剂 ;打开试剂槽盖板 ,使用无尘纸擦净冷凝水 .
j) 提前 1 h取出 dNTPs混合液 Ⅲ和 dNTPs混合液 Ⅱ ,室温融化后置于冰上或 4 ℃备用 ;加样前需使用涡旋振荡器振荡 5 s混匀 ,短暂离心(离心到底部即可)后使用 ;使用前取出 DNA 聚合酶混合液 ,置于冰上或 4 ℃备用 ,加样前需颠倒混匀 4 次 ~ 6 次 .
k) 使用洁净的 1 mL枪头在测序试剂槽的 1 号和 2 号孔边缘位置轻轻戳出一个直径约 2 cm 的加样孔位 ;在 1 号孔位中加入 0 .96 mL dNTPs混合液 Ⅲ和 0 .96 mL DNA 聚合酶混合液 ;在 2 号孔位中加入 2 .04 mL dNTPs混合液 Ⅱ 和 1 .02 mL DNA 聚合酶混合液 ;使用配套的透明封口膜将 1 号和 2 号加样孔封住(切勿盖住孔位中心位置 ,避免影响试剂针下降) .
l) 测序试剂槽水平放置在桌面上 ,双手握住两侧 ,顺时针摇晃 10 次 ~ 20 次 ,再逆时针摇晃 10 次 ~ 20 次 ,其间要确保肉眼可见旋涡 ,直至该试剂槽中 1 号孔位中的试剂上下层颜色均匀一致 , 以保证试剂充分混匀 .
m) 使用枪头戳破 15 号孔的封口膜 ;用 200 μL移液器移取 200 μL MDA 聚合酶混合液加入 MDA试剂的试剂管中 ;颠倒混匀 4 次 ~ 6 次 ,使其充分混匀 ,再将混匀液加入 15 号孔中 ,加入时确保管底部无气泡 ,至此即完成测序试剂槽上机前的准备工作 .
n) 上机测序 . 从 - 20 ℃冰箱中取出载片包装彩盒 ,将载片从中取出 ,拆开真空包装袋 ; 打开载片舱
出 ;用空气罐吹净载片平台和载片背面的灰尘(如果平台表面有可见结晶 ,需要用润湿的无尘纸
轻轻擦拭) ,按下载片吸附按钮 ,取出新的载片 ,两孔位置在左侧 ,一孔位置在右侧 ,标签位置靠
右 ,双手握住载片两端 ;载片孔位对应定位柱放置 ,保持载片空位内壁与定位柱贴合 ,将载片边框
左右两边同时按下 ,使载片吸附在平台上 ,确保载片可以牢固吸附 ,关闭载片舱门 .
o) 打开试剂舱舱门 ,按照试剂槽盖板指示方向 ,把准备好的测序试剂槽轻轻推进试剂舱 ,直到推到底部并确认测序试剂槽完全放入 ;放入要加载的 DNB冻存管 ,关闭试剂舱舱门 .
p) 在计算机软件界面上点击 “测序 ”,进入测序参数设置界面 ,输入 DNB编号 ,选择测序方案 “FCL
PE150 ”,点击 “下一步 ”,将光标放置在 “试剂槽 ID”文本框 ,打开试剂舱舱门 ,使用条码扫描枪扫描测序试剂槽条码录入试剂槽信息 ,关闭试剂舱舱门 ;把光标移至 “载片 ID”后面的文本框 ,打开载片舱舱门 ,扫描载片上的二维码录入载片信息 ;各项信息确认无误后 ,点击 “开始 ”;待测序完成后 ,点击 “完成 ”,并将测序数据拷贝至移动硬盘 .
C.5 测序数据拆分
高通量 测 序 仪 根 据 样 品 条 形 码 的 序 列 , 自 动 将 测 序 数 据 拆 分 到 样 品 XMI220424200 和 样 品XMI220424070 .
由于采用双末端测序模式 , 因此 ,每个样品的每个测序片段均包括正向和反向测序序列 . 其中 ,样品XMI220424200 的 正 向 和 反 向 测 序 序 列 存 放 文 件 的 名 称 分 别 为 XMI220424200 _ 1 .fq.gz 和XMI220424200_ 2 .fq.gz; 样 品 XMI220424070 的 正 向 和 反 向 测 序 序 列 存 放 文 件 的 名 称 分 别 为XMI220424070_1 .fq.gz和 XMI220424070_2 .fq.gz.
C.6 测序数据比对
小麦参考基因组版本为 iwgsc_refseqv1 .fasta.
iuoie(il:..对iilw:c:_/-f.e_f.t)f,a,
..scI2上2口z””–式,文;cs–._序0线1I220数程fz;
对结果输入 XMI220424200 .sam 文件中 .
比对结果采用 SAM(The Sequence Alignment/ Map format,序列比对格式)格式保存 . SAM 文件格式详细说明见 http://samtools.github.io/hts-specs/SAMv1 .pdf.
C.7 测序数据质量控制
若测序片段比对到参考基因组上的位置与标记位点在参考基因组上的位置重合 ,则判定该测序片段属于该标记位点 .
统计每个标记位点的测序片段的数目 ,作为该标记位点的覆盖倍数 . 例如 ,样品 XMI220424200 的第1 个标记位点覆盖倍数为 568 倍 .
按上述方法 ,获得 XMI220424200 所有标记位点的覆盖倍数 ,计算 XMI220424200 标记位点的平均覆
盖倍1≥(5,数 ) /780= 960 .24 倍 .
每个标记位点所有相同测序片段归为同一个等位基因型 ,统计该标记位点的每个等位基因型的测序片段的数目 . 例如 ,表 C.1 为样品 XMI220424200 在表 A.1 中的第 2 个 MNP标记位点的所有等位基因型序列及其测序片段数目 .
评论