绿竹（Dendrocalamopsis oldhami（Munro））属禾本科（Gramineae）竹亚科（Bambusoideae）绿竹属（Dendrocalamopsis）丛生竹种，是绿竹属中分布最广、栽培最多的竹种^[1]。绿竹具有发笋期长、生长迅速、产量丰富、笋味鲜美，质地柔软的优点，是我国南方优良速生的笋材两用丛生竹种之一，具有较高的商品效益，是主产地的经济来源之一^[2]。绿竹原产我国台湾省淡水，我国主要分布在浙江南部、福建、台湾、广东、广西和海南等省区^[3]。

目前，国内关于绿竹的研究报道，主要集中在良种繁育技术研究、生物量研究、生理特性研究，蛋白组学研究等方面^[4-10]，而关于分子遗传学研究^[11-13]较少，且绿竹分子标记的研究中主要是运用随机扩增多态性DNA（Random amplified polymorphic DNA，RAPD）技术，余学军等^[11]利用RAPD标记对22个不同栽培类型的绿竹进行多态性分析。吴益民等^[12]构建孝顺竹、凤尾竹、绿竹、白绿竹4个竹种RAPD指纹图谱。但是RAPD在实验过程中容易产生假带、共迁移等现象，所以结果不够稳定，重复性较低。简单重复序列区间（Inter simpce sequence repeat，ISSR）技术是1994年Zietkiewicz等^[14]在简单重复序列（Simple sequence repeat，SS-R）分子标记技术的基础上发展起来的一种新的分子标记方法，具有简单快捷易操作，DNA用量较少，试验成本低等优点，由于其多态性好、稳定性高和重复性高的特性被广泛应用于多种物种的品种鉴定和遗传多样性分析等方面^{[15, 16]}。而利用ISSR分子标记技术对绿竹ISSR-PCR体系建立尚未有报道。本实验采用正交试验设计的方法对dNTPs浓度、Mg²⁺浓度、Taq DNA聚合酶浓度、引物浓度和模板DNA浓度设计5因素4水平试验，建立绿竹ISSR-PCR最佳反应体系，并优化扩增程序、筛选出高多态性引物，旨在为绿竹的种质资源的鉴定、基因定位、遗传多样性分析、分子育种、DNA指纹图谱构建等研究奠定了基础。

1 材料与方法 1.1 材料

供试的材料于2016年8月采自于福建省东山国有赤山林场的绿竹地理种源试验样地。对不同的植株随机选取3-5片嫩叶，放置于塑料自封袋中，并用变色硅胶进行迅速干燥，将干燥的材料速带回实验室，于-80℃超低温冰箱保存备用。

1.2 方法 1.2.1 绿竹基因组总DNA的提取与检测

绿竹总DNA的提取采用的是杭州博日公司Biospin植物基因组DNA提取试剂盒（Cat#BSC13S1）法，并对该种方法略作改良。提取完后用微量紫外分光光度计测定DNA的浓度和纯度。并用1%的琼脂糖凝胶电泳检测DNA的完整性与清晰性，看是否有降解或断裂现象，将满足试验条件的DNA稀释到20 ng/μL，并置于-20℃冰箱保存。

1.2.2 绿竹ISSR-PCR扩增反应体系的优化

用L16（4⁵）正交试验设计方法对dNTPs浓度、Mg²⁺浓度、Taq DNA聚合酶浓度、引物浓度和模板DNA浓度设计5因素4水平试验（表 1），共16个处理，每个处理2个重复，来研究各因素对ISSR-PCR扩增的影响。反应体系的总体积为20 μL，引物选择UBC815。扩增程序为：94℃预变性5 min；94℃变性45 s，（根据引物的退火温度）复性30 s，72℃延伸90 s，循环35次，72℃延伸10 min，4℃保存。

将PCR产物于1.5%的琼脂糖凝胶上，电泳45 min，电压设为5 v/cm。并用紫外凝胶成像分析仪拍照记录电泳情况。利用直观分析法对每个处理结果进行评分，一共划分为16个等级，以条带的清晰度、明亮度、特异性高低为评分标准，最高的记为16分，最低的记为1分^[17]。评分结果用DPS7.05统计软件进行方差分析。

1.2.3 绿竹ISSR-PCR扩增反应程序的优化

在最优反应体系的基础上，采用单因素试验的方法对退火温度和循环次数进行筛选，根据引物U815的Tm值52.9℃，在45-60℃之间设置了12个温度梯度，系统自动生成的12个温度梯度分别为44.9℃，45.2℃，46.1℃，47.5℃，49.3℃，51.2℃，53.3℃，55.2℃，57.1℃，58.5℃，59.6℃，60.1℃。对循环次数设置了30次、32次、35次、38次、40次5个梯度，以筛选出最佳的ISSR-PCR扩增反应程序。

1.2.4 绿竹ISSR扩增体系稳定性验证

以不同种源的绿竹的基因组总DNA为模板，在最佳反应体系和最佳扩增程序下，用引物UBC815进行ISSR-PCR扩增，检验建立的扩增体系的稳定性。

1.2.5 ISSR有效引物的筛选

以优化的反应体系及扩增程序体系为基础进行引物筛选，在100条ISSR引物中筛选出14条扩增条带清晰、多态性较高、重复性好的引物。

2 结果 2.1 绿竹基因组总DNA的提取结果

部分绿竹基因组总DNA凝胶电泳结果如图 1所示，我们可以看出DNA条带清晰明亮，无明显拖尾现象，说明所提取的DNA基因组完整性较好，无明显降解现象。利用微量紫外分光光度计测定DNA的A₂₆₀/A₂₈₀的比值均在1.8-2.0之间，说明DNA的纯度较高。

2.2 绿竹ISSR-PCR正交试验结果

绿竹正交试验的电泳结果如图 2所示，在16个处理组合中，由于dNTPs、Mg²⁺、Taq DNA聚合酶、引物和模板DNA浓度5个影响因素各水平组合的不同，扩增结果存在明显的差异。对正交试验的各处理的扩增结果进行评分，评分结果见表 1。根据评分计算每个因素同一水平下的总和T和平均值X，以及极差R，计算结果见表 2。

极差R值越大，说明该因素对绿竹ISSR-PCR扩增的影响越大。由表 2可知，各因素对ISSR-PCR扩增结果的影响从大到小依次为：Mg²⁺ > dNTPs > 模板DNA > Taq DNA聚合酶 > 引物。平均值X最大值所对应的水平即为该因素的最佳浓度。所以dNTPs水平3、Mg²⁺水平1、Taq DNA聚合酶水平3、引物水平2、模板DNA水平4最好。因此绿竹ISSR-PCR正交试验得出的最佳水平为：dNTPs浓度为0.2 mmol/L，Mg²⁺浓度为2.0 mmol/L，TaqDNA聚合酶浓度为1.5 U，引物浓度为0.5 μmol/L，DNA浓度为为60 ng。

为了进一步判断各个因素对绿竹ISSR-PCR扩增影响的显著程度，对正交试验的结果进行方差分析，结果见表 3。其中F值越大，说明该因素对绿竹ISSR-PCR扩增的影响越大。由表 3可知，各因素对ISSR-PCR扩增结果的影响从大到小依次为：Mg²⁺ > dNTPs > 模板DNA > Taq DNA聚合酶 > 引物。与极差分析的结果一致。其中dNTPs、Mg²⁺和模板DNA对绿竹ISSR-PCR扩增影响达到极显著的水平（P < 0.01）。Taq DNA聚合酶对绿竹ISSR-PCR扩增影响达到显著的水平（P < 0.05）。引物浓度对绿竹ISSR-PCR扩增影响不显著。

为了进一步确定每个因素的最适宜的浓度水平，对5个因素进行Duncan多重比较，结果见表 4。由之前的极差分析和方差分析可知，dNTPs对绿竹ISSR-PCR扩增的影响仅次于Mg²⁺，从多重比较可以看出，dNTPs水平1、水平2、水平3和水平4差异显著，而水平1、水平2、水平3之间的差异并不显著，从经济角度最终选择dNTPs最佳浓度为0.2 mmol/L。Mg²⁺对绿竹ISSR-PCR扩增的影响最大，Mg²⁺水平1、水平2和水平3、水平4差异显著，水平1和水平2两者差异不显著，综合而言，最终选择Mg²⁺最佳浓度为2.0 mmol/L。当TaqDNA聚合酶浓度为1.5 U时，极差分析的均值最大，且与其他3个水平差异显著，所以最终选择TaqDNA聚合酶最佳浓度为1.5 U。在本试验中，引物浓度对绿竹ISSR-PCR扩增的影响最小，且4个水平间差异不显著，从经济角度最终选择引物最佳浓度为0.4 μmol/L。当模板DNA浓度为60 ng，极差分析的均值最大，且与其他3个水平差异显著，所以最终选模板DNA最佳用量为60 ng。

综上所述，最终确定的绿竹ISSR-PCR扩增最佳反应体积为：dNTPs浓度为0.2 mmol/L，Mg²⁺浓度为2.0 mmol/L，TaqDNA聚合酶浓度为1.5 U，引物浓度为0.4 μmol/L，DNA浓度为60 ng。

2.3 绿竹ISSR-PCR扩增反应程序的优化结果

不同的循环次数对ISSR-PCR扩增有一定的影响。本实验对循环次数设置了30次、32次、35次、38次和40次5个梯度，实验结果如图 3所示。当循环次数为38次时，扩增条带最为清晰明亮且稳定。所以，选择38次循环为ISSR-PCR反应的最佳循环次数。

不同的退火温度对ISSR-PCR扩增也有一定的影响。根据引物U815的Tm值52.9℃，在45-60℃之间设置了12个温度梯度，结果（图 4）显示退火温度对扩增影响明显，当退火温度为47.5℃时，扩增的条带数最为清晰稳定，当退火温度大于47.5℃时，高分子量区域和中等分子量区域条带的扩增效率下降，直至无扩增条带，且在低分子量区域也逐渐出现弥散现象，整体电泳条带数逐渐减少，清晰度也逐渐减弱，因此，最终确定引物U815的最佳退火温度为47.5℃。

2.4 绿竹ISSR扩增体系稳定性验证

利用引物UBC818对不同种源的绿竹个体的基因组总DNA样品进行进行ISSR-PCR扩增，扩增结果见图 5。由图 5可知，在该体系下均能扩增出条带清晰、稳定，多态性丰富的条带。说明该体系稳定可靠，适合绿竹ISSR-PCR扩增反应。

2.5 ISSR有效引物的筛选

根据绿竹优化后的最佳反应体系和扩增程序，对100条随机引物进行扩增。部分引物筛选结果如图 6所示。结果显示利用建立的体系筛选不同的引物时，大部分引物都能扩增出清晰明亮的条带，只有少部分引物扩增出的条带数较少或不能扩增出条带。说明本实验建立的绿竹最佳反应体系具有较高的重复性和稳定性。最终在100条引物中筛选出14条扩增条带清晰不拖尾，稳定性较好的引物（表 5）。

3 讨论

ISSR分子标记技术已经广泛应用于种质资源的鉴定、基因定位、遗传多样性分析、指纹图谱的构建等方面的研究，是目前使用最广泛的分子标记技术之一^[18-20]。ISSR-PCR扩增结果受到物种类型、反应体系、扩增程序、引物等的综合影响。本试验采用正交试验设计的方法对dNTPs浓度、Mg²⁺浓度、Taq DNA聚合酶浓度、引物浓度、模板DNA用量设计5因素4水平试验，并对扩增程序进行优化，利用优化后的体系对100条ISSR引物进行筛选。正交试验设计可以克服单因素试验不能考虑到各影响因素之间的相互作用的缺点，还可以解决完全组合设计需要大量试验、耗时耗费的不足，具有设计均衡分散、齐整可比、效应明显、节约时间和成本等优点^[17]，能够通过较少的处理组合快速找到最佳的ISSR-PCR反应体系，并能分析出不同因素对扩增结果影响的显著性。但是正交试验直观分析又有一定的主观随机性，所以需要设置重复试验来减少打分误差，增加数据的可靠性^[21]。

在本实验中，不同的处理组合对绿竹ISSR-PCR扩增结果的影响差异明显。极差分析和方差分析的结果均表明各因素对绿竹ISSR-PCR扩增结果的影响从大到小依次为：Mg²⁺ > dNTPs > 模板DNA > Taq DNA聚合酶 > 引物。Mg²⁺是影响绿竹ISSR-PCR扩增最显著的因素，Mg²⁺浓度变化直接影响Taq DNA聚合酶的活性，并间接影响引物与模板DNA的双联杂交体的解链温度、引物的退火温度、扩增产物的特异性以及引物二聚体的生成^[22-23]。靳晓丽等^[24-26]研究均显示Mg²⁺对ISSR-PCR扩增影响最大，所以正确的把握好Mg²⁺的用量很关键。dNTPs对绿竹ISSR-PCR扩增的影响仅次于Mg²⁺，dNTPs是PCR扩增的底物，dNTPs浓度过低会降低DNA的合成速率，导致扩增产量的下降；dNTPs浓度过高又会与Mg²⁺相鳌合从而间接影响Taq DNA聚合酶的活性，导致碱基的错误渗入和非特异性扩增^{[27, 28]}。模板DNA浓度对绿竹ISSR-PCR扩增的影响也是极显著的，在本实验中DNA模板浓度越高，扩增的条带越亮产量越多。Taq DNA聚合酶对绿竹ISSR-PCR扩增的影响是显著的，Taq DNA聚合酶浓度过低会导致扩增效率较低，条带较弱；浓度过高又会产生非特异性条带^[29]。在本实验中引物浓度对绿竹ISSR-PCR扩增的影响不显著，这可能与本实验中引物浓度设置的水平范围有关，在本实验中，引物浓度在0.4-0.7 μmol/L之间设置了4个水平，4个水平间差异不显著，与包蕊等^[30]研究的在ISSR-PCR反应中引物浓度影响最小的研究结果一致。

循环次数和退火温度对绿竹ISSR-PCR扩增也有一定影响，循环次数较少，扩增的条带数也较少且清晰度较低；循环次数过多又会增加非特异性扩增。退火温度较低会导致模板DNA与引物的非特异性结合且背景深；退火温度较高又会使得模板DNA与引物的结合率下降，从而使得扩增出来的条带数较少^[31]。在一定的温度范围内，随着退火温度的增加，扩增出的产物的特异性越好，一般最佳的退火温度低于引物的真实Tm值5℃左右^[32]。

4 结论

本研究采用L16（4⁵）正交试验设计方法对dNTPs浓度、Mg²⁺浓度、Taq DNA聚合酶浓度、引物浓度和模板DNA浓度设计5因素4水平试验，并通过直观分析和正交极差分析对比，确定了这5个因素对ISSR-PCR扩增的影响力大小以及最佳用量。最后确定适合绿竹ISSR-PCR的最佳反应体系为：dNTPs浓度为0.2 mmol/L，Mg²⁺浓度为2.0 mmol/L，TaqDNA聚合酶浓度为1.5 U，引物浓度为0.4 μmol/L，DNA浓度为60 ng。验证实验表明，该体系稳定可靠，适合于绿竹ISSR-PCR扩增。