柚木无性系指纹图谱的构建方法及其应用的制作方法
【专利摘要】本发明公开一种柚木无性系指纹图谱的构建方法及其应用。柚木无性系指纹图谱的构建方法，包括如下步骤：柚木无性系基因组DNA提取，基因组DNA的SSR扩增，PCR扩增产物上测序仪分型，获得柚木无性系指纹图谱；本发明通过试验研究，优化了柚木SSR分子标记PCR扩增体系和程序，设计和筛选出柚木SSR多态性引物。本发明用筛选出来的15对具有多态性的柚木SSR引物，构建了26个柚木无性系的DNA指纹图谱。该方法经济、精确、快速和高效，为今后快速进行无性系指纹图谱构建和鉴定、柚木种质资源遗传评价、遗传图谱构建、分子标记辅助育种提供了技术支撑。
【专利说明】柚木无性系指纹图谱的构建方法及其应用
【技术领域】
[0001]本发明涉及一种乔木无性系指纹图谱的构建方法及其应用，具体涉及一种柚木无性系指纹图谱的构建方法及其应用。
【背景技术】
[0002]袖木(Tectonagrandis L.f.)为马鞭草科(Verbenaceae)袖木属(Tectona)高大乔木，天然分布于北纬9°~26°，东经73°~104°的印度南部与中部、泰国北部、緬甸和老挝，是世界著名的热带珍贵用材树种，素有“万木之王”的美誉，其心材材质坚韧耐久，结构致密，纹理美观，是制造军舰、码头、桥梁、建筑、家具、雕刻和木地板等的高级用材。全球木材贸易网数据显示，2011年6月16~30日緬甸柚木市场的4类原木离岸价都已达到3471美元.m-3。国内市场进口柚木原木价格较高，杭州木材市场柚木价格为9000~18000元.m-3。正是具有生长快，材质优良，用途广和投资回报率高等特点，柚木在热带、南亚热带地区广为种植，是世界上人工林种植面积最大的4个树种之一，也是单位面积产值最高的造林树种。
[0003]我国引进柚木栽培已有180多年的历史，20世纪70年代初我国开始了柚木遗传改良和培育技术的系统研究，为我国发展柚木事业提供了栽培材料和技术的保障。目前，我国柚木推广种植范围遍及10个省(自治区)60多个县(市)，总面积约1.5万hm2。近年来，云南省、海南省把柚木列为主要造林树种，广东、广西、福建、贵州省采用柚木造林的私人与企业也越来越多。柚木商品林种植多以组培优良无性系为主,选育出高产、优质和多抗的柚木优良新品系并加以推广 利用，提高良种化程度，是大力发展柚木种植业的前提。我国柚木研究人员经过多年的努力，已获得了一批适合我国不同立地条件和气候的柚木优良无性系，授权指定了个别组培单位进行生产繁殖和推广。
[0004]但是，目前市场上出售的柚木苗繁殖方法多样，繁殖材料来源不清，苗木质量良莠不齐，而现有技术凭人工经验难以辨别出柚木不同无性系，导致我国柚木商品林发展良种使用率低，在很大程度上制约了我国柚木事业的发展。
【发明内容】
[0005]为了克服现有技术的缺点与不足，本发明的目的在于提供一种柚木无性系指纹图谱的构建方法。
[0006]本发明的另一目的在于提供通过上述方法构建的柚木无性系指纹图谱在柚木无性系品种鉴定上的应用。本发明通过研制一套柚木无性系鉴定方法，构建出我国柚木无性系指纹图谱，为柚木无性系的鉴定和品种权保护提供技术支撑，并可有效监督和规范我国各地柚木种植苗木的良种信息情况，提高柚木商品林发展良种化程度。
[0007]本发明的目的通过下述技术方案实现:一种柚木无性系指纹图谱的构建方法，包括如下步骤:
[0008](I)柚木无性系基因组DNA提取[0009]以柚木无性系的叶片为材料，提取柚木无性系的基因组DNA，得到基因组DNA样品;
[0010](2)柚木无性系基因组DNA的SSR扩增
[0011]以步骤(1)中得到的基因组DNA样品为模板，用柚木SSR多态性引物进行PCR扩增，得到PCR扩增产物；
[0012](3)柚木无性系PCR扩增产物上测序仪分型
[0013]取步骤(2)中获得的PCR扩增产物，加入超纯甲酰胺和LIZ-500分子量内标，得到上样液；将上样液变性，降温，然后进行柚木无性系SSR分型；
[0014](4)获得柚木无性系的指纹图谱
[0015]用GeneMapper4.0软件进行谱带分析，读取每对引物扩增出的SSR片段的大小(单位bp)，得到柚木SSR多态性引物构建出的柚木无性系指纹图谱(见表7)。
[0016]步骤(1)中所述的提取柚木无性系的基因组DNA优选采用CTAB法提取；
[0017]步骤(1)中所述的基因组DNA样品优选是放于_20°C保存备用；
[0018] 步骤(2)中所述的柚木SSR多态性引物包括15对柚木SSR多态性引物，分别是 AJ968929、AJ968930、AJ968931、AJ968932、AJ968933、AJ968934、AJ968935、AJ968936、AJ968938、AJ968940、AJ968941、AJ968942、AJ968943、AJ511753、AJ511764，其核苷酸序列依次如序列表SEQ ID N0.1~30所示；
[0019]步骤(2)中所述的PCR扩增的反应体系优选为:含浓度为25mM Mg2+的IOXbufferl.5 μ L,浓度为 IOmM 的 dNTPs0.03 μ L,酶活为 5U/μ L 的 Taq DNApolymerase(Taq DNA 酶)0.3 μ L,浓度为 IOng/ μ L 的 DNA1.5 μ L,浓度为 5 μ M 的 forward primer (正向引物)0.15 μ L,浓度为5 μ M的reverse primer (反向引物)0.15 μ L,浓度为Inmol/μ L的 fluorescent-dUTP (荧光-dUTP) 0.015 μ L,加 ddH20 至 15 μ L ;
[0020]步骤(2)中所述的PCR扩增的反应程序优选为:采用Touch-Down PCR扩增程序:94°C预变性4min ;94°C变性30s、Tm+10°C至Tm°C退火30s、72°C延伸lmin，20个循环，每个循环降低0.50C ;94°C变性30s、Tm°C退火30s、72°C延伸lmin，26个循环；72°C延伸IOmin ；4 °C forever ；
[0021]所述的Tm (退火温度)优选为51、53或56°C ；
[0022]步骤(3)中所述的上样液的体系优选为:PCR扩增产物2yL，超纯甲酰胺7.82yL，LIZ-500分子量内标0.18 μ L ;
[0023]步骤(3)中所述的变性的条件优选为95°C变性5min ；
[0024]步骤(3)中所述的降温的条件优选是放入冰浴或4°C中进行降温4分钟；降温后尽快上样分析；
[0025]步骤(3)中所述的进行柚木无性系SSR分型优选是将变性后PCR产物在ABI3130X1测序仪上进行柚木无性系SSR分型；
[0026]步骤(5)中所述的读取每对引物扩增出的SSR片段的大小，其中，不同无性系能得到不同或相同大小的SSR片段，也有的无性系不含有个别引物扩增的SSR片段，因此扩增没有产物；
[0027]—种柚木无性系指纹图谱，通过上述方法构建得到，如表7所示。
[0028]所述的柚木无性系指纹图谱在柚木无性系品种鉴定上的应用。[0029]利用柚木无性系指纹图谱进行柚木品种鉴别；依据15对柚木SSR多态性引物的检测结果进行判定:利用上述构建方法，提取未知柚木无性系的基因组DNA，进行PCR扩增，将扩增产物上测序仪分型，获得每对引物扩增出的SSR片段大小(单位bp)，若与上述柚木无性系指纹图谱中的片段完全一致，即获得该未知柚木无性系的品系名称。
[0030]本发明相对于现有技术，具有如下的优点及效果:
[0031](I)本发明通过EMBL-EBI网站找到了柚木基因组DNA微卫星序列，利用Serafer等相关软件设计出了柚木SSR引物，通过试验研究，筛选出了 15对柚木SSR多态性引物，确定了每对引物的退火温度；同时优化了柚木SSR分子标记PCR扩增体系和程序，最终获得了基于荧光-dUTP和自动化测序仪的柚木SSR分子标记技术体系。该方法经济、精确、快速和高效，为今后快速进行柚木种质资源遗传评价、遗传图谱构建、分子标记辅助育种、无性系指纹图谱构建和鉴定提供了技术支撑。
[0032](2)利用上述方法，用适宜数量的具有多态性的柚木SSR引物，对柚木无性系基因组DNA分别进行扩增和上测序仪进行SSR分型，由于不同无性系所含的微卫星片段重复数不同，从而对不同无性系扩增分析得到的微卫星片段大小不同，由此构建出无性系的指纹图谱。本发明用筛选出来的15对具有多态性的柚木SSR引物，构建了 26个柚木无性系的DNA指纹图谱，可把26个柚木无性系区分开，将为柚木无性系的鉴定和品种权保护提供技术支撑，同时可有效监督和规范我国各地柚木种植苗木的良种信息情况，提高柚木商品林发展良种化程度。
【专利附图】
【附图说明】
[0033] 图1是26个柚木无性系的遗传聚类图。
[0034]图2为5对引物29、30、31、32、33扩增柚木无性系70-13的峰图；
[0035]其中，29为 AJ968929、30 为 AJ968930、31 为 AJ96893U32 为 AJ968932、33 为AJ968933。
[0036]图3为5对引物34、35、36、38、40扩增柚木无性系70-13的峰图；
[0037]其中34 为 AJ968934、35 为 AJ968935、36 为 AJ968936、38 为 AJ968938、40 为AJ968940。
[0038]图4为5对引物41、42、43、53、64扩增柚木无性系70-13的峰图；
[0039]其中，41为 AJ96894U42 为 AJ968942、43 为 AJ968943、53 为 AJ511753、64 为AJ511764。
[0040]图5为5对引物29、30、31、32、33扩增柚木无性系70-29的峰图；
[0041]其中，29为 AJ968929、30 为 AJ968930、31 为 AJ96893U32 为 AJ968932、33 为AJ968933。
[0042]图6为5对引物34、35、36、38、40扩增柚木无性系70-29的峰图；
[0043]其中34 为 AJ968934、35 为 AJ968935、36 为 AJ968936、38 为 AJ968938、40 为AJ968940。
[0044]图7为5对引物41、42、43、53、64扩增柚木无性系70-29的峰图；
[0045]其中，41为 AJ96894U42 为 AJ968942、43 为 AJ968943、53 为 AJ511753、64 为AJ511764。【具体实施方式】
[0046]下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
[0047]以下材料在文献“抽木种源主要性状聚合遗传值的评价.林业科学研究.1996，9 (I): 7-14”公开:印度种源3070、印度种源3071、印度种源3072、印度种源3074和尼日利亚种源3078。
[0048]以下材料在文献“农杆菌介导柚木腋芽原基遗传转化的研究.中国林业科学研究院.博士学位论文.2007”公开:无性系7531、无性系7544、无性系7552、无性系8301和无性系108。 [0049]以下材料在文献“柚木无性系促萌与采穗圃营建技术.中国林业科学研究院.硕士学位论文.2007”公开:无性系7514和无性系7559。
[0050]以下材料在文献“柚木愈伤组织诱导培养的研究.中南林业科技大学学报.2012，32(3):53-58” 公开:无性系 7114。
[0051]以下材料在文献“柚木种子园无性系开花特性与结实差异分析.种子.2011，30(8):5-8”公开:无性系70-13、无性系70-29、无性系71-5、无性系71-22、无性系71-38、无性系71-46、无性系72-10、无性系72-19和无性系7549。
[0052]以下材料在文献“柚木嫩枝扦插.中南林学院学报.2005，25 (3) ”公开:无性系71-12。
[0053]无性系7137、无性系7412、无性系J731、无性系Z408、无性系7509、无性系7542、无性系7555和无性系3078-5，均购自中国林业科学研究院热带林业研究所。
[0054]实施例1
[0055](I)柚木无性系基因组DNA提取。
[0056]柚木无性系70-13 (已在“柚木种子园无性系开花特性与结实差异分析”公开)的叶片材料，采用CTAB法提取柚木无性系70-13的基因组DNA，将DNA样品放于_20°C保存备用。
[0057]所述的CTAB法的具体步骤如下:
[0058](a)取柚木无性系叶片0.2g于研钵中，液氮冷冻下迅速研磨成粉末。
[0059](b)将粉末转入1.5mL离心管中，加入600 μ L65°C预热的CTAB提取液，加入20 μ LRNase溶液(10mg/mL),充分混匀。
[0060](c)将离心管放入65°C水浴中温浴30min，每隔7~8分钟混匀I次。
[0061]Cd)取出离心管，加入600 μ L的Tris饱和酌.:氯仿(I:1),轻轻混匀5min后4~25°C、12000rpm,离心 IOmin0
[0062](e)取上清液于新离心管中，加入等体积的氯仿:异戊醇(24:1)，轻轻混匀5min后4 ~25°C、12000rpm,离心 IOmin0
[0063](f)取上清液于新离心管中，加入预冷的I倍体积的预冷异丙醇，轻轻摇匀，室温放置2~3min。
[0064](g) 4~25°C、12000rpm,离心lOmin。倒去上清液,用75%的酒精颠倒漂洗沉淀2~3次。将酒精倒出后放入通风橱，把酒精挥发干后加入50 μ LI X TE溶解DNA。[0065]所述的CTAB 提取液包括如下成分:100mmol/L Tris-HCl (pH8.0), 20mmol/L EDTA(乙二胺四乙酸二钠)，1.4mol/L NaCl,2%CTAB (w/v)，2%PVP (聚乙烯毗咯烷酮)，1%β -琉基乙醇(每次提取之前混合均匀后使用)。
[0066]所述的IXTE 包括如下成分:10mmol/L Tris-HCl、lmmol/L EDTA, ΡΗ=8.0。
[0067](2)柚木无性系基因组DNA的SSR扩增。
[0068]以步骤(1)中得 到的基因组DNA样品为模板，用筛选出的15对柚木SSR多态性引物(见表3)分别进行PCR扩增，PCR扩增体系见表2。
[0069]柚木SSR多态性引物筛选的具体步骤如下:
[0070]从EMBL网站上搜索到143条柚木基因组DNA微卫星，用SSRhunter软件分别找到其SSR序列的重复单元和重复数，再根据每条SSR序列两端的侧翼序列，用Primerf.0、01igo6.0或Serafer软件设计出对应的143对SSR引物。从26个柚木无性系(见表6)中任选5个柚木无性系基因组DNA样品混合后作为样品，用所设计的引物分别进行PCR扩增。PCR体系见表1。
[0071]表1柚木SSR引物预筛选PCR体系
[0072]
【权利要求】
1.一种柚木无性系指纹图谱的构建方法，其特征在于包括如下步骤: (1)柚木无性系基因组DNA提取 以柚木无性系的叶片为材料，提取柚木无性系的基因组DNA，得到基因组DNA样品； (2)柚木无性系基因组DNA的SSR扩增 以步骤(1)中得到的基因组DNA样品为模板，用柚木SSR多态性引物进行PCR扩增，得到PCR扩增产物； (3)柚木无性系PCR扩增产物上测序仪分型 取步骤(2)中获得的PCR扩增产物，加入超纯甲酰胺和LIZ-500分子量内标，得到上样液；将上样液变性， 降温，然后进行柚木无性系SSR分型； (4)获得柚木无性系的指纹图谱用GeneMapper4.0软件进行谱带分析,读取每对引物扩增出的SSR片段的大小，得到柚木SSR多态性引物构建出的柚木无性系指纹图谱； 步骤(2 )中所述的柚木S S R多态性引物包括15对柚木S S R多态性引物，分别是AJ968929、AJ968930、AJ968931、AJ968932、AJ968933、AJ968934、AJ968935、AJ968936、AJ968938、AJ968940、AJ968941、AJ968942、AJ968943、AJ511753、AJ511764，其核苷酸序列依次如序列表SEQ ID N0.1~30所示。
2.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(1)中所述的基因组DNA样品是放于_20°C保存备用。
3.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(2)中所述的PCR扩增的反应体系为:含浓度为25mM Mg2+的IOXbufferl.5 μ L，浓度为IOmM的dNTPs0.03 μ L,酶活为 5U/ μ L 的 Taq DNA 酶 0.3 μ L，浓度为 IOng/ μ L 的 DNAL 5 μ L,浓度为5μ M的正向引物0.15yL，浓度为5μΜ的反向引物0.15 μ L，浓度为Inmol/μ L的荧光-dUTP0.015 μ L，加 ddH20 至 15 μ L。
4.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(2)中所述的PCR扩增的反应程序为:采用Touch-Down PCR扩增程序:94°C预变性4min ;94°C变性30s、Tm+10°C至Tm°C退火30s、72°C延伸lmin，20个循环,每个循环降低0.5°C;94°C变性30s、Tm°C退火 30s、72°C延伸 lmin, 26 个循环；72°C延伸 IOmin ；4°C forever ； 所述的Tm为51、53或56°C。
5.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(3)中所述的上样液的体系为:PCR扩增产物2yL，超纯甲酰胺7.82 μ L, LIZ-500分子量内标0.18 μ L。
6.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(3)中所述的变性的条件为95°C变性5min。
7.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(3)中所述的降温是放入冰浴或4°C中进行降温4分钟。
8.根据权利要求1所述的柚木无性系指纹图谱的构建方法，其特征在于:步骤(3)中所述的进行柚木无性系SSR分型是将变性后PCR产物在ABI3130X I测序仪上进行柚木无性系SSR分型。
9.通过权利要求1~8任一项所述的构建方法得到的柚木无性系指纹图谱在柚木无性系品种鉴定上的应用。
10.根据权利要求9所述的应用，其特征在于:利用柚木无性系指纹图谱进行柚木品种鉴别；依据15对柚木SSR多态性引物的检测结果进行判定:利用权利要求1~8任一项所述的构建方法，提取未知柚木无性系的基因组DNA，进行PCR扩增，将扩增产物上测序仪分型，获得每对引物扩增出的SSR片段的大小，若与上述柚木无性系指纹图谱中的片段完全一致，即获得该未知柚木无性系的`品系名称。
【文档编号】C12Q1/68GK103725785SQ201410010871
【公开日】2014年4月16日 申请日期:2014年1月9日 优先权日:2014年1月9日
【发明者】黄桂华, 梁坤南, 周再知, 马华明 申请人:中国林业科学研究院热带林业研究所