本发明属于分子生物学技术领域，尤其涉及一种植物DNA甲基化组织特异性获取方法。
背景技术：
DNA甲基化在植物发育过程以及基因组防御环境因子胁迫等方面扮演重要的作用，在植物中，DNA甲基化通常发生在CpG或者CpNpG序列上，并且在DNA修复的循环中可以维持不变，所有的植物发育过程都建立在精准的时空特异表达和基因表达调控的基础之上。越来越多的证据表明基因表达的时空特异性不仅由特异的DNA序列(顺式以及反式作用元件控制)调控，还由一些表观修饰因子，其中最主要的就是DNA甲基化进行调控，但目前现有技术中还未见有专门针对DNA甲基化组织特异性的获取方法的报道。
技术实现要素：
有鉴于此，本发明的目的在于针对现有技术的不足，而提供一种植物DNA甲基化组织特异性的获取方法，实现了植物组织DNA甲基化特异性的定量分析。
本发明提供了一种植物DNA甲基化组织特异性的获取方法，包括以下步骤：
1)提取植物组织基因组DNA；
2)在同一酶切体系中，用同裂酶和EcoRI对所述步骤1)得到的植物组织基因组DNA进行酶切，得到酶切产物，再对酶切产物连接限制性内切酶的接头，得到连接产物；所述同裂酶为Hpa II或Msp I；
3)以所述步骤2)接头为模板设计预扩增引物，再以所述预扩增引物对所述步骤2)得到连接产物进行PCR预扩增，得PCR预扩增产物；
4)将所述步骤3)PCR预扩增产物稀释后，加入带有选择性碱基的筛选引物进行PCR扩增，得到PCR扩增产物；
5)电泳检测所述步骤4)得到的PCR扩增产物，统计DNA条带，将(1,0)记为H，H代表半甲基化修饰的位点，将(0,1)记为M，M代表全甲基化修饰的位点，(0,0)和(1,1)记为N，N代表无效信息位点和非甲基化位点，对统计结果赋值，使H＝M＝1，N＝0；
6)根据所述步骤5)的赋值，利用公式1得到DNA甲基化片段的组织特异性：
公式1中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值；Ti等于0代表在全部组织中发生甲基化，Ti等于1代表仅在一个组织中发生了甲基化；
利用公式2得到DNA甲基化片段组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性，Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式，Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
优选的，所述步骤2)中的接头包括EcoRI接头和同裂酶接头，所述EcoRI接头的序列为SEQ ID No:1或SEQ ID No:2，所述同裂酶接头的序列为SEQ ID No:3或SEQ ID No:4。
优选的，所述步骤3)预扩增引物包括EcoRI预扩增引物和同裂酶预扩增引物，所述EcoRI预扩增引物的序列为SEQ ID No:5，所述同裂酶预扩增引物的序列为SEQ IDNo:6。
优选的，所述步骤4)筛选引物包括EcoRI筛选引物和同裂酶筛选引物，所述EcoRI筛选引物序列为SEQ ID No:7～70中的一种，所述同裂酶筛选引物序列为SEQ ID No:71～134中的一种。
优选的，所述EcoRI筛选引物序列与同裂酶筛选引物序列的组合为SEQ ID No:16+SEQ ID No:132、SEQ ID No:17+SEQ ID No:78、SEQ ID No:22+SEQ ID No:133、SEQ ID No:29+SEQ ID No:80、SEQ ID No:29+SEQ ID No:86、SEQ ID No:29+SEQ ID No:93、SEQ ID No:29+SEQ ID No:94、SEQ ID No:29+SEQ ID No:109、SEQ ID No:29+SEQ ID No:113、SEQ ID No:29+SEQ ID No:127、SEQ ID No:29+SEQ ID No:133、SEQ ID No:29+SEQ ID No:134、SEQ ID No:29+SEQ ID No:77、SEQ ID No:29+SEQ ID No:78、SEQ ID No:30+SEQ ID No:86、SEQ ID No:36+SEQ ID No:81、SEQ ID No:45+SEQ ID No:81、SEQ ID No:45+SEQ ID No:86、SEQ ID No:45+SEQ ID No:100、SEQ ID No:45+SEQ ID No:110、SEQ ID No:45+SEQ ID No:113、SEQ ID No:45+SEQ ID No:125、SEQ ID No:46+SEQ ID No:80、SEQ ID No:46+SEQ ID No:93、SEQ ID No:46+SEQ ID No:109、SEQ ID No:46+SEQ ID No:110、SEQ ID No:46+SEQ ID No:113、SEQ ID No:46+SEQ ID No:125、SEQ ID No:46+SEQ ID No:127、SEQ ID No:46+SEQ ID No:78、SEQ ID No:48+SEQ ID No:93、SEQ ID No:49+SEQ ID No:94、SEQ ID No:61+SEQ ID No:127、SEQ ID No:62+SEQ ID No:113、SEQ ID No:68+SEQ ID No:125、SEQ ID No:69+SEQ ID No:112、SEQ ID No:70+SEQ ID No:100、SEQ ID No:13+SEQ ID No:109或SEQ ID No:14+SEQ ID No:134。
优选的，所述步骤2)中酶切连接用体系为：每20μl体系包括10×扩增缓冲液2.0μl，植物组织基因组DNA4.0μl，EcoRI 0.25μL，HpaII 0.4μL或MspI0.4μL，EcoRI接头0.5μL，同裂酶接头0.5μL，T4连接酶0.5μL，ddH2O11.85μL。
优选的，所述步骤3)中预扩增用体系为：每50μL体系包括稀释100倍的所述连接产物5μL，EcoRI预扩增引物2.5μL，同裂酶预扩增引物2.5μL，10×扩增缓冲液5μL，Taq DNA聚合酶1μL，dNTPs 0.5μL，ddH2O 33.5μL；
所述预扩增的程序为：94℃预变性2min，94℃变性30s，56℃退火30s，72℃延伸80s，30～35个循环，72℃延伸5min。
优选的，所述步骤4)中扩增用体系为：每20μL扩增体系包括稀释100倍的所述PCR预扩增产物4μL，同裂酶引物2μL，EcoRI引物2μL，Taq DNA聚合酶1μL，dNTPs 0.5μL，10×扩增缓冲液2.5μL，ddH2O 8μL。
优选的，所述步骤4)扩增的程序为：94℃预变性5min，94℃变性30s，65℃退火30s，72℃延伸1min，13个循环，每个循环退火温度降0.7℃；
所述13个循环后，94℃变性30s，55℃退火30s，72℃延伸80s，23～25个循环，72℃延伸5min。
本发明提供了一种植物DNA甲基化组织特异性的获取方法，包括以下步骤：1)提取植物组织基因组DNA；2)在同一酶切体系中，用同裂酶和EcoRI对所述步骤1)得到的植物组织基因组DNA进行酶切，得到酶切产物，再对酶切产物连接限制性内切酶的接头，得到连接产物；所述同裂酶为Hpa II或MspI；3)以所述步骤2)接头为模板设计预扩增引物，再以所述预扩增引物对所述步骤2)得到连接产物进行PCR预扩增，得PCR预扩增产物；4)将所述步骤3)PCR预扩增产物稀释后，加入带有选择性碱基的筛选引物进行PCR扩增，得到PCR扩增产物；5)电泳检测所述步骤4)得到的PCR扩增产物，统计DNA条带，将(1,0)记为H，H代表半甲基化修饰的位点，将(0,1)记为M，M代表全甲基化修饰的位点，(0,0)和(1,1)记为N，N代表无效信息位点和非甲基化位点，对统计结果赋值，使H＝M＝1，N＝0；6)根据所述步骤5)的赋值，利用公式1得到DNA甲基化片段的组织特异性：
公式1中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值；Ti等于0代表在全部组织中发生甲基化，Ti等于1代表仅在一个组织中发生了甲基化；
利用公式2得到DNA甲基化片段组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性，Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式，Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。采用本发明的技术方案通过对标记位点赋值结合计算公式实现了植物组织DNA甲基化特异性的定量分析，为后续的组织特异性的DNA甲基化调控基因表达研究提供了技术支持。
附图说明
图1为植物组织DNA甲基化特异性MASP分析电泳图，其中a-j分别代表毛白杨的嫩叶、成熟叶、韧皮部、形成层、根、顶点、未成熟木质部、成熟木质部、雄花序和雌花序共10个组织。(1,0)，(0,1)，(0,0)和(1,1)分别代表半甲基化、全甲基化、无效信息位点以及非甲基化位点。
具体实施方式
本发明提供了一种植物DNA甲基化组织特异性的获取方法，包括以下步骤：
1)提取植物组织基因组DNA；
2)在同一酶切体系中，用同裂酶和EcoRI对所述步骤1)得到的植物组织基因组DNA进行酶切，得到酶切产物，再对酶切产物连接限制性内切酶的接头，得到连接产物；所述同裂酶为Hpa II或MspI；
3)以所述步骤2)接头为模板设计预扩增引物，再以所述预扩增引物对所述步骤2)得到连接产物进行PCR预扩增，得PCR预扩增产物；
4)将所述步骤3)PCR预扩增产物稀释后，加入带有选择性碱基的筛选引物进行PCR扩增，得到PCR扩增产物；
5)电泳检测所述步骤4)得到的PCR扩增产物，统计DNA条带，将(1,0)记为H，H代表半甲基化修饰的位点，将(0,1)记为M，M代表全甲基化修饰的位点，(0,0)和(1,1)记为N，N代表无效信息位点和非甲基化位点，对统计结果赋值，使H＝M＝1，N＝0；
6)根据所述步骤5)的赋值，利用公式1得到DNA甲基化片段的组织特异性：
公式I中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值；Ti等于0代表在全部组织中发生甲基化，Ti等于1代表仅在一个组织中发生了甲基化；
利用公式2得到DNA甲基化片段组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性，Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式，Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
在本发明中，所述植物组织优选包括嫩叶、成熟叶、韧皮部、形成层、根、顶点分生组织、未成熟木质部、成熟木质部、雄花序、雌花序、叶片或茎。
本发明对所述植物组织基因组DNA的提取方法没有特殊限定，采用本领域技术人员熟知的基因组DNA提取的方法即可。
得到植物组织基因组DNA后，在同一酶切体系中，用同裂酶和EcoRI对得到的植物组织基因组DNA进行酶切，得到酶切产物，再对酶切产物连接限制性内切酶的接头，得到连接产物。在本发明中，所述酶切体系的反应条件优选为14～18℃过夜，更优选为16℃过夜。
在本发明中，所述同裂酶为Hpa II或MspI。本发明对所述同裂酶和EcoRI的来源没有特殊限定，采用本领域技术人员常规选用的市售商品即可。
在本发明中，所述接头优选包括EcoRI接头和同裂酶接头，所述EcoRI接头的序列为SEQ ID No:1或SEQ ID No:2，所述同裂酶接头的序列为SEQ ID No:3或SEQ ID No:4。
在本发明中，所述酶切连接用体系优选为：每20μl体系包括10×扩增缓冲液2.0μl，植物组织基因组DNA4.0μl，EcoRI 0.25μL，HpaII 0.4μL或MspI0.4μL，EcoRI接头0.5μL，同裂酶接头0.5μL，T4连接酶0.5μL，ddH2O 11.85μL。本发明对上述试剂的来源没有特殊限定，采用本领域技术人员常用的市售商品即可。
在本发明中，以接头为模板设计预扩增引物，再以所述预扩增引物对得到连接产物进行PCR预扩增，得PCR预扩增产物。在本发明中，所述预扩增引物优选包括EcoRI预扩增引物和同裂酶预扩增引物，所述EcoRI预扩增引物的序列为SEQ ID No:5，所述同裂酶预扩增引物的序列为SEQ ID No:6。
在本发明中，所述预扩增用体系优选为：每50μL体系包括稀释100倍的所述连接产物5μL，EcoRI预扩增引物2.5μL，同裂酶预扩增引物2.5μL，10×扩增缓冲液5μL，Taq DNA聚合酶1μL，dNTPs 0.5μL，ddH2O 33.5μL。本发明对上述试剂的来源没有特殊限定，采用本领域技术人员常规选用的市售商品即可。
在本发明中，所述预扩增的程序优选为：94℃预变性2min，94℃变性30s，56℃退火30s，72℃延伸80s，30～35个循环，72℃延伸5min。
在本发明中，将所述PCR预扩增产物稀释后，加入带有选择性碱基的筛选引物进行PCR扩增，得到PCR扩增产物。
在本发明中，所述PCR预扩增产物稀释的倍数优选为80～120倍，更优选为90～110倍，最优选为100倍。
在本发明中，所述筛选引物优选包括EcoRI筛选引物和同裂酶筛选引物，所述EcoRI筛选引物序列为SEQ ID No:7～70中的一种，所述同裂酶筛选引物序列为SEQ IDNo:71～134中的一种。在本发明中，所述EcoRI筛选引物序列与同裂酶筛选引物序列的组合优选为SEQ ID No:16+SEQ ID No:132、SEQ ID No:17+SEQ ID No:78、SEQ ID No:22+SEQ ID No:133、SEQ ID No:29+SEQ ID No:80、SEQ ID No:29+SEQ ID No:86、SEQ ID No:29+SEQ ID No:93、SEQ ID No:29+SEQ ID No:94、SEQ ID No:29+SEQ ID No:109、SEQ ID No:29+SEQ ID No:113、SEQ ID No:29+SEQ ID No:127、SEQ ID No:29+SEQ ID No:133、SEQ ID No:29+SEQ ID No:134、SEQ ID No:29+SEQ ID No:77、SEQ ID No:29+SEQ ID No:78、SEQ ID No:30+SEQ ID No:86、SEQ ID No:36+SEQ ID No:81、SEQ ID No:45+SEQ ID No:81、SEQ ID No:45+SEQ ID No:86、SEQ ID No:45+SEQ ID No:100、SEQ ID No:45+SEQ ID No:110、SEQ ID No:45+SEQ ID No:113、SEQ ID No:45+SEQ ID No:125、SEQ ID No:46+SEQ ID No:80、SEQ ID No:46+SEQ ID No:93、SEQ ID No:46+SEQ ID No:109、SEQ ID No:46+SEQ ID No:110、SEQ ID No:46+SEQ ID No:113、SEQ ID No:46+SEQ ID No:125、SEQ ID No:46+SEQ ID No:127、SEQ ID No:46+SEQ ID No:78、SEQ ID No:48+SEQ ID No:93、SEQ ID No:49+SEQ ID No:94、SEQ ID No:61+SEQ ID No:127、SEQ ID No:62+SEQ ID No:113、SEQ ID No:68+SEQ ID No:125、SEQ ID No:69+SEQ ID No:112、SEQ IDNo:70+SEQ ID No:100、SEQ ID No:13+SEQ ID No:109或SEQ IDNo:14+SEQ IDNo:134。
在本发明中，所述扩增用体系优选为：每20μL扩增体系包括稀释100倍的所述PCR预扩增产物4μL，同裂酶引物2μL，EcoRI引物2μL，Taq DNA聚合酶1μL，dNTPs 0.5μL，10×扩增缓冲液2.5μL，ddH2O 8μL。本发明对上述试剂的来源没有特殊限定，采用本领域技术人员常规的市售商品即可。
在本发明中，所述扩增的程序优选为：94℃预变性5min，94℃变性30s，65℃退火30s，72℃延伸1min，13个循环，每个循环退火温度降0.7℃；所述13个循环后，94℃变性30s，55℃退火30s，72℃延伸80s，23～25个循环，72℃延伸5min。
在本发明中，将得到的PCR扩增产物进行电泳检测。在本发明中，所述电泳检测的条件优选为采用1800～2000V恒压进行5～9％聚丙烯酰胺凝胶电泳1～2h。所述电泳采用的电压优选为1800～2000V，更优选为1850～1950V，最优选为1900V。所述电泳用的聚丙烯酰胺凝胶中琼脂的质量体积百分含量优选为5～9％，更优选为6～8％，最优选为7％。所述电泳的时间优选为1～2h，更优选为1.5h。
本发明在所述电泳后，统计DNA条带，将(1,0)记为H，H代表半甲基化修饰的位点，将(0,1)记为M，M代表全甲基化修饰的位点，将(0,0)和(1,1)记为N，N代表无效信息位点和非甲基化位点，对统计结果赋值，使H＝M＝1，N＝0。本发明根据上述赋值，利用公式1得到DNA甲基化片段的组织特异性：
公式I中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值；Ti等于0代表在全部组织中发生甲基化，Ti等于1代表仅在一个组织中发生了甲基化；
利用公式2得到DNA甲基化片段组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性，Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式，Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
下面结合实施例对本发明提供的一种植物DNA甲基化组织特异性获取方法进行详细的说明，但是不能把它们理解为对本发明保护范围的限定。
实施例1
1、利用本领域常规植物组织DNA提取方法，提取毛白杨的嫩叶、成熟叶、韧皮部、形成层、根、顶点分生组织、未成熟木质部、成熟木质部、雄花序和雌花序共10个组织基因组DNA；
2、分别用Msp I和EcoRI两组酶对植物各组织基因组DNA进行酶切连接，得到连接产物；
双酶切连接反应体系的总体积为20μl，具体如下：
16℃连接过夜。
3、以接头为模板设计的预扩增引物，再以预扩增引物对连接产物进行PCR预扩增，得PCR预扩增产物；
预扩增反应体系为50μL，具体如下：
预扩增的程序为：94℃预变性2min，94℃变性30s，56℃退火30s，72℃延伸80s，30～35个循环，72℃延伸5min。
4、将PCR预扩增产物稀释100倍后，加入带有选择性碱基的筛选引物进行PCR扩增；
选择性扩增体系为20μL，具体如下：
扩增的程序为：94℃预变性5min，94℃变性30s，65℃退火30s，72℃延伸1min，13个循环，每个循环退火温度降0.7℃；
所述13个循环后，94℃变性30s，55℃退火30s，72℃延伸80s，23～25个循环，72℃延伸5min。
5、聚丙烯酰胺凝胶电泳检测PCR扩增产物，采用1900V恒压进行7％聚丙烯酰胺凝胶电泳1.5h，银染，电泳结果见图1。统计并分析DNA条带；聚丙烯酰胺凝胶电泳所涉及各种试剂在《分子克隆试验指南》中有详细配置方法。统计条带，将(1,0)记为H，(0,1)记为M，(0，0)和(1,1)记为N。对统计结果赋值，使H＝M＝1，N＝0。
6、利用公式1和公式2计算植物组织基因组DNA甲基化组织特异性和偏好性。
利用公式1计算植物组织DNA甲基化片段的组织特异性：
公式1中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织，M＝1，H＝1，M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值。Ti等于0代表在全部组织中发生甲基化。Ti等于1代表仅在一个组织中发生了甲基化。
利用公式2计算DNA甲基化片段的组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织。M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性。Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式。Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
以毛白杨的嫩叶、成熟叶、韧皮部、形成层、根、顶点、未成熟木质部、成熟木质部、雄花序和雌花序10个组织为研究材料，在4096对引物组合中筛选得到39对扩增稳定并且效率较高的引物用于毛白杨DNA甲基化组织特异性分析，共得到了642个位点，选择了28个位点进行克隆测序，共得到了26条序列，经与杨树基因组数据库比对得到22个基因，并计算了这些基因的甲基化修饰在毛白杨各组织中的特异性以及偏好性，结果见表1和表2。
表1毛白杨MSAP片段在10个组织中的甲基化模式
从表1中可以得出，毛白杨MSAP片段在10个组织中存在迥异的DNA甲基化模式。
表2毛白杨MSAP片段功能注释及组织特异性
从表2中可以得出，毛白杨MSAP片段代表了不同的候选基因，这些基因在10个组织中的甲基化模式特异性和偏好性差别显著。
实施例2
1、利用本领域常规DNA提取方法，提取桉树叶片、根、茎、花序和未成熟木质部5个组织基因组DNA；
2、分别用Hpa II和EcoRI两组酶对植物各组织基因组DNA进行酶切连接，得到连接产物；
双酶切连接反应体系的总体积为20μl，具体如下：
16℃连接过夜。
3、以接头为模板设计的预扩增引物，再以预扩增引物对连接产物进行PCR预扩增，得PCR预扩增产物；
预扩增反应体系为50μL，具体如下：
预扩增的程序为：94℃预变性2min，94℃变性30s，56℃退火30s，72℃延伸80s，30～35个循环，72℃延伸5min
4、将所述PCR预扩增产物稀释后，加入带有选择性碱基的筛选引物进行PCR扩增；
选择性扩增体系为20μL，具体如下：
扩增的程序为：94℃预变性5min，94℃变性30s，65℃退火30s，72℃延伸1min，13个循环，每个循环退火温度降0.7℃；
所述13个循环后，94℃变性30s，55℃退火30s，72℃延伸80s，23～25个循环，72℃延伸5min。
5、聚丙烯酰胺凝胶电泳检测PCR扩增产物，采用1900V恒压进行7％聚丙烯酰胺凝胶电泳1.5h，银染。统计并分析DNA条带；聚丙烯酰胺凝胶电泳所涉及各种试剂在《分子克隆试验指南》中有详细配置方法。统计条带，将(1,0)记为H，(0,1)记为M，(0，0)和(1,1)记为N(见表2)。对统计结果赋值，使H＝M＝1，N＝0。
6、利用公式1和公式2计算植物组织基因组DNA甲基化组织特异性和偏好性。
利用公式1计算植物组织DNA甲基化片段的组织特异性：
公式1中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织。M＝1，H＝1，M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值。Ti等于0代表在全部组织中发生甲基化。Ti等于1代表仅在一个组织中发生了甲基化。
利用公式2计算DNA甲基化片段的组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织。M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性。Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式。Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
以桉树叶片、根、茎、花序和未成熟木质部5个组织为研究材料，在4096对引物组合中筛选得到45对扩增稳定并且效率较高的引物用于桉树DNA甲基化组织特异性分析，共得到了831个位点，选择了22个位点进行克隆测序，共得到了18条序列，经与桉树基因组数据库比对得到18个基因，并计算了这些基因的甲基化修饰在桉树各组织中的特异性以及偏好性，结果见表3和
表3桉树MSAP片段在5个组织中的甲基化模式
从表3中可以得出，桉树MSAP片段在5个组织中存在迥异的DNA甲基化模式。
表4桉树MSAP片段功能注释及组织特异性
从表4中可以得出，桉树MSAP片段代表了不同的候选基因，这些基因在5个组织中的甲基化模式特异性和偏好性差别显著。
实施例3
1、利用本领域常规植物组织DNA提取方法，提取拟南芥叶片、茎、根、花和顶点分生组5个组织基因组DNA；
2、分别用Hpa II和EcoRI两组酶对植物各组织基因组DNA进行酶切连接，得到连接产物；
双酶切连接反应体系的总体积为20μl，具体如下：
16℃连接过夜。
3、以所述接头为模板设计的预扩增引物，再以预扩增引物对连接产物进行PCR预扩增，得PCR预扩增产物；
预扩增反应体系为50μL，具体如下：
预扩增的程序为：94℃预变性2min，94℃变性30s，56℃退火30s，72℃延伸80s，30～35个循环，72℃延伸5min。
4、将所述PCR预扩增产物稀释100倍后，加入带有选择性碱基的筛选引物进行PCR扩增；
选择性扩增体系为20μL，具体如下：
扩增的程序为：94℃预变性5min，94℃变性30s，65℃退火30s，72℃延伸1min，13个循环，每个循环退火温度降0.7℃；
所述13个循环后，94℃变性30s，55℃退火30s，72℃延伸80s，23～25个循环，72℃延伸5min。
5、聚丙烯酰胺凝胶电泳检测PCR扩增产物，采用1900V恒压进行7％聚丙烯酰胺凝胶电泳1.5h，银染。统计并分析DNA条带；聚丙烯酰胺凝胶电泳所涉及各种试剂在《分子克隆试验指南》中有详细配置方法。统计条带，将(1,0)记为H，(0,1)记为M，(0，0)和(1,1)记为N(见表2)。对统计结果赋值，使H＝M＝1，N＝0。
6、利用公式1和公式2计算植物组织基因组DNA甲基化组织特异性和偏好性。
利用公式1计算植物组织DNA甲基化片段的组织特异性：
公式1中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织。M＝1，H＝1，M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Ti代表第i个DNA甲基化片段的甲基化模式组织特异性值。Ti等于0代表在全部组织中发生甲基化。Ti等于1代表仅在一个组织中发生了甲基化。
利用公式2计算DNA甲基化片段的组织特异性的偏好性：
公式2中，n代表分析使用植物组织的数量，i代表第i个DNA甲基化片段，j代表第j个组织。M代表全甲基化修饰的位点，H代表半甲基化修饰的位点。Bi代表第i个DNA甲基化片段的甲基化模式在组织中偏好性。Bi在[-1,0)之间代表该片段在组织中甲基化模式偏好为全甲基化模式。Bi在(0,1]之间代表该片段在组织中甲基化模式偏好为半甲基化模式，Bi＝0时代表该片段在组织中甲基化模式没有偏好性。
以拟南芥叶、茎、根、花和顶点5个组织为研究材料，在4096对引物组合中筛选得到48对扩增稳定并且效率较高的引物用于拟南芥DNA甲基化组织特异性分析，共得到了712个位点，选择了29个位点进行克隆测序，共得到了24条序列，经与拟南芥基因组数据库比对得到18个基因，并计算了这些基因的甲基化修饰在拟南芥各组织中的特异性以及偏好性。结果见表5和表6。
表5拟南芥MSAP片段在5个组织中的甲基化模式
从表5中可以得出，拟南芥MSAP片段在5个组织中存在迥异的DNA甲基化模式。
表6拟南芥MSAP片段功能注释及组织特异性
从表6中可以得出，拟南芥MSAP片段代表了不同的候选基因，这些基因在5个组织中的甲基化模式特异性和偏好性差别显著。
从以上实施例可以得出，采用本发明的技术方案实现了植物组织DNA甲基化特异性的定量分析，为后续的组织特异性的DNA甲基化调控基因表达研究提供了技术支持。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
SEQUENCE LISTING
<110> 北京林业大学
<120> 一种植物DNA甲基化组织特异性获取方法
<130> 2017
<160> 134
<170> PatentIn version 3.3
<210> 1
<211> 17
<212> DNA
<213> 人工序列
<400> 1
ctcgtagact gcgtacc 17
<210> 2
<211> 15
<212> DNA
<213> 人工序列
<400> 2
aattggtacg cagtc 15
<210> 3
<211> 16
<212> DNA
<213> 人工序列
<400> 3
gatcatgagt cctgct 16
<210> 4
<211> 16
<212> DNA
<213> 人工序列
<400> 4
cgagcaggac tcatga 16
<210> 5
<211> 16
<212> DNA
<213> 人工序列
<400> 5
gactgcgtac caattc 16
<210> 6
<211> 18
<212> DNA
<213> 人工序列
<400> 6
atcatgagtc ctgctcgg 18
<210> 7
<211> 19
<212> DNA
<213> 人工序列
<400> 7
gactgcgtac caattcaaa 19
<210> 8
<211> 19
<212> DNA
<213> 人工序列
<400> 8
gactgcgtac caattcaac 19
<210> 9
<211> 19
<212> DNA
<213> 人工序列
<400> 9
gactgcgtac caattcaat 19
<210> 10
<211> 19
<212> DNA
<213> 人工序列
<400> 10
gactgcgtac caattcaag 19
<210> 11
<211> 19
<212> DNA
<213> 人工序列
<400> 11
gactgcgtac caattcaca 19
<210> 12
<211> 19
<212> DNA
<213> 人工序列
<400> 12
gactgcgtac caattcact 19
<210> 13
<211> 19
<212> DNA
<213> 人工序列
<400> 13
gactgcgtac caattcacg 19
<210> 14
<211> 19
<212> DNA
<213> 人工序列
<400> 14
gactgcgtac caattcacc 19
<210> 15
<211> 19
<212> DNA
<213> 人工序列
<400> 15
gactgcgtac caattcata 19
<210> 16
<211> 19
<212> DNA
<213> 人工序列
<400> 16
gactgcgtac caattcatc 19
<210> 17
<211> 19
<212> DNA
<213> 人工序列
<400> 17
gactgcgtac caattcatg 19
<210> 18
<211> 19
<212> DNA
<213> 人工序列
<400> 18
gactgcgtac caattcatt 19
<210> 19
<211> 19
<212> DNA
<213> 人工序列
<400> 19
gactgcgtac caattcaga 19
<210> 20
<211> 19
<212> DNA
<213> 人工序列
<400> 20
gactgcgtac caattcagc 19
<210> 21
<211> 19
<212> DNA
<213> 人工序列
<400> 21
gactgcgtac caattcagt 19
<210> 22
<211> 19
<212> DNA
<213> 人工序列
<400> 22
gactgcgtac caattcagg 19
<210> 23
<211> 19
<212> DNA
<213> 人工序列
<400> 23
gactgcgtac caattccaa 19
<210> 24
<211> 19
<212> DNA
<213> 人工序列
<400> 24
gactgcgtac caattccac 19
<210> 25
<211> 19
<212> DNA
<213> 人工序列
<400> 25
gactgcgtac caattccat 19
<210> 26
<211> 19
<212> DNA
<213> 人工序列
<400> 26
gactgcgtac caattccag 19
<210> 27
<211> 19
<212> DNA
<213> 人工序列
<400> 27
gactgcgtac caattccca 19
<210> 28
<211> 19
<212> DNA
<213> 人工序列
<400> 28
gactgcgtac caattccct 19
<210> 29
<211> 19
<212> DNA
<213> 人工序列
<400> 29
gactgcgtac caattcccg 19
<210> 30
<211> 19
<212> DNA
<213> 人工序列
<400> 30
gactgcgtac caattcccc 19
<210> 31
<211> 19
<212> DNA
<213> 人工序列
<400> 31
gactgcgtac caattccta 19
<210> 32
<211> 19
<212> DNA
<213> 人工序列
<400> 32
gactgcgtac caattcctc 19
<210> 33
<211> 19
<212> DNA
<213> 人工序列
<400> 33
gactgcgtac caattcctg 19
<210> 34
<211> 19
<212> DNA
<213> 人工序列
<400> 34
gactgcgtac caattcctt 19
<210> 35
<211> 19
<212> DNA
<213> 人工序列
<400> 35
gactgcgtac caattccga 19
<210> 36
<211> 19
<212> DNA
<213> 人工序列
<400> 36
gactgcgtac caattccgc 19
<210> 37
<211> 19
<212> DNA
<213> 人工序列
<400> 37
gactgcgtac caattccgt 19
<210> 38
<211> 19
<212> DNA
<213> 人工序列
<400> 38
gactgcgtac caattccgg 19
<210> 39
<211> 19
<212> DNA
<213> 人工序列
<400> 39
gactgcgtac caattctaa 19
<210> 40
<211> 19
<212> DNA
<213> 人工序列
<400> 40
gactgcgtac caattctac 19
<210> 41
<211> 19
<212> DNA
<213> 人工序列
<400> 41
gactgcgtac caattctat 19
<210> 42
<211> 19
<212> DNA
<213> 人工序列
<400> 42
gactgcgtac caattctag 19
<210> 43
<211> 19
<212> DNA
<213> 人工序列
<400> 43
gactgcgtac caattctca 19
<210> 44
<211> 19
<212> DNA
<213> 人工序列
<400> 44
gactgcgtac caattctct 19
<210> 45
<211> 19
<212> DNA
<213> 人工序列
<400> 45
gactgcgtac caattctcg 19
<210> 46
<211> 19
<212> DNA
<213> 人工序列
<400> 46
gactgcgtac caattctcc 19
<210> 47
<211> 19
<212> DNA
<213> 人工序列
<400> 47
gactgcgtac caattctta 19
<210> 48
<211> 19
<212> DNA
<213> 人工序列
<400> 48
gactgcgtac caattcttc 19
<210> 49
<211> 19
<212> DNA
<213> 人工序列
<400> 49
gactgcgtac caattcttg 19
<210> 50
<211> 19
<212> DNA
<213> 人工序列
<400> 50
gactgcgtac caattcttt 19
<210> 51
<211> 19
<212> DNA
<213> 人工序列
<400> 51
gactgcgtac caattctga 19
<210> 52
<211> 19
<212> DNA
<213> 人工序列
<400> 52
gactgcgtac caattctgc 19
<210> 53
<211> 19
<212> DNA
<213> 人工序列
<400> 53
gactgcgtac caattctgt 19
<210> 54
<211> 19
<212> DNA
<213> 人工序列
<400> 54
gactgcgtac caattctgg 19
<210> 55
<211> 19
<212> DNA
<213> 人工序列
<400> 55
gactgcgtac caattcgaa 19
<210> 56
<211> 19
<212> DNA
<213> 人工序列
<400> 56
gactgcgtac caattcgac 19
<210> 57
<211> 19
<212> DNA
<213> 人工序列
<400> 57
gactgcgtac caattcgat 19
<210> 58
<211> 19
<212> DNA
<213> 人工序列
<400> 58
gactgcgtac caattcgag 19
<210> 59
<211> 19
<212> DNA
<213> 人工序列
<400> 59
gactgcgtac caattcgca 19
<210> 60
<211> 19
<212> DNA
<213> 人工序列
<400> 60
gactgcgtac caattcgct 19
<210> 61
<211> 19
<212> DNA
<213> 人工序列
<400> 61
gactgcgtac caattcgcg 19
<210> 62
<211> 19
<212> DNA
<213> 人工序列
<400> 62
gactgcgtac caattcgcc 19
<210> 63
<211> 19
<212> DNA
<213> 人工序列
<400> 63
gactgcgtac caattcgta 19
<210> 64
<211> 19
<212> DNA
<213> 人工序列
<400> 64
gactgcgtac caattcgtc 19
<210> 65
<211> 19
<212> DNA
<213> 人工序列
<400> 65
gactgcgtac caattcgtg 19
<210> 66
<211> 19
<212> DNA
<213> 人工序列
<400> 66
gactgcgtac caattcgtt 19
<210> 67
<211> 19
<212> DNA
<213> 人工序列
<400> 67
gactgcgtac caattcgga 19
<210> 68
<211> 19
<212> DNA
<213> 人工序列
<400> 68
gactgcgtac caattcggc 19
<210> 69
<211> 19
<212> DNA
<213> 人工序列
<400> 69
gactgcgtac caattcggt 19
<210> 70
<211> 19
<212> DNA
<213> 人工序列
<400> 70
gactgcgtac caattcggg 19
<210> 71
<211> 21
<212> DNA
<213> 人工序列
<400> 71
atcatgagtc ctgctcggaa a 21
<210> 72
<211> 21
<212> DNA
<213> 人工序列
<400> 72
atcatgagtc ctgctcggaa c 21
<210> 73
<211> 21
<212> DNA
<213> 人工序列
<400> 73
atcatgagtc ctgctcggaa t 21
<210> 74
<211> 21
<212> DNA
<213> 人工序列
<400> 74
atcatgagtc ctgctcggaa g 21
<210> 75
<211> 21
<212> DNA
<213> 人工序列
<400> 75
atcatgagtc ctgctcggac a 21
<210> 76
<211> 21
<212> DNA
<213> 人工序列
<400> 76
atcatgagtc ctgctcggac t 21
<210> 77
<211> 21
<212> DNA
<213> 人工序列
<400> 77
atcatgagtc ctgctcggac g 21
<210> 78
<211> 21
<212> DNA
<213> 人工序列
<400> 78
atcatgagtc ctgctcggac c 21
<210> 79
<211> 21
<212> DNA
<213> 人工序列
<400> 79
atcatgagtc ctgctcggat a 21
<210> 80
<211> 21
<212> DNA
<213> 人工序列
<400> 80
atcatgagtc ctgctcggat c 21
<210> 81
<211> 21
<212> DNA
<213> 人工序列
<400> 81
atcatgagtc ctgctcggat g 21
<210> 82
<211> 21
<212> DNA
<213> 人工序列
<400> 82
atcatgagtc ctgctcggat t 21
<210> 83
<211> 21
<212> DNA
<213> 人工序列
<400> 83
atcatgagtc ctgctcggag a 21
<210> 84
<211> 21
<212> DNA
<213> 人工序列
<400> 84
atcatgagtc ctgctcggag c 21
<210> 85
<211> 21
<212> DNA
<213> 人工序列
<400> 85
atcatgagtc ctgctcggag t 21
<210> 86
<211> 21
<212> DNA
<213> 人工序列
<400> 86
atcatgagtc ctgctcggag g 21
<210> 87
<211> 21
<212> DNA
<213> 人工序列
<400> 87
atcatgagtc ctgctcggca a 21
<210> 88
<211> 21
<212> DNA
<213> 人工序列
<400> 88
atcatgagtc ctgctcggca c 21
<210> 89
<211> 21
<212> DNA
<213> 人工序列
<400> 89
atcatgagtc ctgctcggca t 21
<210> 90
<211> 21
<212> DNA
<213> 人工序列
<400> 90
atcatgagtc ctgctcggca g 21
<210> 91
<211> 21
<212> DNA
<213> 人工序列
<400> 91
atcatgagtc ctgctcggcc a 21
<210> 92
<211> 21
<212> DNA
<213> 人工序列
<400> 92
atcatgagtc ctgctcggcc t 21
<210> 93
<211> 21
<212> DNA
<213> 人工序列
<400> 93
atcatgagtc ctgctcggcc g 21
<210> 94
<211> 21
<212> DNA
<213> 人工序列
<400> 94
atcatgagtc ctgctcggcc c 21
<210> 95
<211> 21
<212> DNA
<213> 人工序列
<400> 95
atcatgagtc ctgctcggct a 21
<210> 96
<211> 21
<212> DNA
<213> 人工序列
<400> 96
atcatgagtc ctgctcggct c 21
<210> 97
<211> 21
<212> DNA
<213> 人工序列
<400> 97
atcatgagtc ctgctcggct g 21
<210> 98
<211> 21
<212> DNA
<213> 人工序列
<400> 98
atcatgagtc ctgctcggct t 21
<210> 99
<211> 21
<212> DNA
<213> 人工序列
<400> 99
atcatgagtc ctgctcggcg a 21
<210> 100
<211> 21
<212> DNA
<213> 人工序列
<400> 100
atcatgagtc ctgctcggcg c 21
<210> 101
<211> 21
<212> DNA
<213> 人工序列
<400> 101
atcatgagtc ctgctcggcg t 21
<210> 102
<211> 21
<212> DNA
<213> 人工序列
<400> 102
atcatgagtc ctgctcggcg g 21
<210> 103
<211> 21
<212> DNA
<213> 人工序列
<400> 103
atcatgagtc ctgctcggta a 21
<210> 104
<211> 21
<212> DNA
<213> 人工序列
<400> 104
atcatgagtc ctgctcggta c 21
<210> 105
<211> 21
<212> DNA
<213> 人工序列
<400> 105
atcatgagtc ctgctcggta t 21
<210> 106
<211> 21
<212> DNA
<213> 人工序列
<400> 106
atcatgagtc ctgctcggta g 21
<210> 107
<211> 21
<212> DNA
<213> 人工序列
<400> 107
atcatgagtc ctgctcggtc a 21
<210> 108
<211> 21
<212> DNA
<213> 人工序列
<400> 108
atcatgagtc ctgctcggtc t 21
<210> 109
<211> 21
<212> DNA
<213> 人工序列
<400> 109
atcatgagtc ctgctcggtc g 21
<210> 110
<211> 21
<212> DNA
<213> 人工序列
<400> 110
atcatgagtc ctgctcggtc c 21
<210> 111
<211> 21
<212> DNA
<213> 人工序列
<400> 111
atcatgagtc ctgctcggtt a 21
<210> 112
<211> 21
<212> DNA
<213> 人工序列
<400> 112
atcatgagtc ctgctcggtt c 21
<210> 113
<211> 21
<212> DNA
<213> 人工序列
<400> 113
atcatgagtc ctgctcggtt g 21
<210> 114
<211> 21
<212> DNA
<213> 人工序列
<400> 114
atcatgagtc ctgctcggtt t 21
<210> 115
<211> 21
<212> DNA
<213> 人工序列
<400> 115
atcatgagtc ctgctcggtg a 21
<210> 116
<211> 21
<212> DNA
<213> 人工序列
<400> 116
atcatgagtc ctgctcggtg c 21
<210> 117
<211> 21
<212> DNA
<213> 人工序列
<400> 117
atcatgagtc ctgctcggtg t 21
<210> 118
<211> 21
<212> DNA
<213> 人工序列
<400> 118
atcatgagtc ctgctcggtg g 21
<210> 119
<211> 21
<212> DNA
<213> 人工序列
<400> 119
atcatgagtc ctgctcggga a 21
<210> 120
<211> 21
<212> DNA
<213> 人工序列
<400> 120
atcatgagtc ctgctcggga c 21
<210> 121
<211> 21
<212> DNA
<213> 人工序列
<400> 121
atcatgagtc ctgctcggga t 21
<210> 122
<211> 21
<212> DNA
<213> 人工序列
<400> 122
atcatgagtc ctgctcggga g 21
<210> 123
<211> 21
<212> DNA
<213> 人工序列
<400> 123
atcatgagtc ctgctcgggc a 21
<210> 124
<211> 21
<212> DNA
<213> 人工序列
<400> 124
atcatgagtc ctgctcgggc t 21
<210> 125
<211> 21
<212> DNA
<213> 人工序列
<400> 125
atcatgagtc ctgctcgggc g 21
<210> 126
<211> 21
<212> DNA
<213> 人工序列
<400> 126
atcatgagtc ctgctcgggc c 21
<210> 127
<211> 21
<212> DNA
<213> 人工序列
<400> 127
atcatgagtc ctgctcgggt a 21
<210> 128
<211> 21
<212> DNA
<213> 人工序列
<400> 128
atcatgagtc ctgctcgggt c 21
<210> 129
<211> 21
<212> DNA
<213> 人工序列
<400> 129
atcatgagtc ctgctcgggt g 21
<210> 130
<211> 21
<212> DNA
<213> 人工序列
<400> 130
atcatgagtc ctgctcgggt t 21
<210> 131
<211> 21
<212> DNA
<213> 人工序列
<400> 131
atcatgagtc ctgctcgggg a 21
<210> 132
<211> 21
<212> DNA
<213> 人工序列
<400> 132
atcatgagtc ctgctcgggg c 21
<210> 133
<211> 21
<212> DNA
<213> 人工序列
<400> 133
atcatgagtc ctgctcgggg t 21
<210> 134
<211> 21
<212> DNA
<213> 人工序列
<400> 134
atcatgagtc ctgctcgggg g 21