杜龙岗，王美兴



玉米SLAF标记的开发及其在玉米果皮纤维素含量 BSA分析中的应用

杜龙岗，王美兴

（浙江省农业科学院作物与核技术利用研究所，杭州 310021）

【目的】降低果皮纤维素是甜玉米品质改良的重要目标，然而玉米果皮纤维素含量调控的研究甚少，相关调控基因尚未定位。利用纤维素含量差异的重组自交系（RILs）群体，通过特异位点扩增片段（specific-locus amplified fragment-sequencing，SLAF）测序和分离混池分析（bulked segregant analysis，BSA）定位控制玉米果皮纤维素含量的染色体区段，鉴定调控玉米果皮纤维素含量的候选基因。【方法】以果皮纤维素含量显著差异的E327和G5-1为亲本，构建重组自交系（RILs）。对RILs群体进行果皮纤维素含量的测定，并根据纤维素含量的结果选择纤维素含量高、低的样本进行混池，用于SLAF标签的鉴定和BSA分析。在BSA分析中，首先对两混池和2个亲本DNA用Ⅲ和166Ⅱ进行酶切，回收414—464 bp的酶切片段进行Illumina建库，并进行SLAF测序，然后根据多态性SLAF标签开发SNP标记，利用SNP标记对玉米果皮纤维素含量进行关联分析，鉴定调控甜玉米果皮纤维素含量的染色体区段。分析这些区段所包含的玉米基因，并找到它们对应的拟南芥同源基因，通过查阅拟南芥相关基因功能研究的文献，进一步鉴定控制玉米果皮纤维素含量的候选基因。【结果】两亲本和2个混池SLAF建库测序得到的SLAF符合预期，基于SLAF测序数据，鉴定了73 786个多态性SLAF标签，这些SLAF标签均匀分布在玉米的10条染色体上。在这些多态性标签中得到了523 395 SNP位点信息。通过关联分析，调控果皮纤维素变异的基因被定位到玉米基因组的6个染色区段，都位于玉米的第5染色体上。在这些区段上，一共有47个玉米基因。通过进一步的研究分析，在这些关联的染色体区段最终确定了9个候选基因。【结论】定位到调控玉米果皮纤维素的含量的基因，表明此方法可以用于基因定位。

玉米；果皮；纤维素；SLAF测序；BSA分析

0 引言

【研究意义】分子标记开发对基因定位和分子标记辅助育种有着重要意义，在玉米中采用SLAF-Seq可以快速开发出大量的多态性标记；果皮是玉米籽粒重要的组成部分之一，与玉米籽粒的抗性、脱水速度等性状相关。对甜玉米而言，果皮柔嫩度是甜玉米食用品质的决定因素之一，影响适口性这一重要的品质性状。果皮的柔嫩度作为衡量甜玉米口感品质的重要性状指标，其与果皮厚度和果皮纤维素含量有着密切的关系。鉴于果皮纤维素含量是玉米果皮柔嫩度的主要化学形成因子，对玉米纤维素含量调控基因区段的定位和候选基因的确定对玉米品质改良将具有重要意义。【前人研究进展】中国玉米种植面积和产量均居世界第二位，仅次于美国，是中国的主要粮食作物之一[1]。尽管中国玉米生产约占世界总产量的20%，但仍需进口大量玉米来满足国内所需。其中，饲料用玉米占总需求的70%以上，其余用作食用，工业等。近年来，中国特用玉米正大力发展，包括甜玉米、糯玉米、爆裂玉米等。特用玉米具有较高的经济价值、营养价值和加工价值[2]。甜玉米粮经果菜兼用，是重要的特用玉米之一[3]。其不仅含有丰富的蛋白质、维生素以及硒、镁、铁等抗癌微量元素，还含有大量的膳食纤维，具有消渴、利尿、解毒，增加肠胃蠕动，治疗便秘和防止结肠癌功效[4-5]。近年来，玉米栽培面积逐年上升，市场对优质、高产甜玉米的需求量不断增加。同时，消费者对甜玉米的品质要求也不断提高[6]。因此，甜度高、口感鲜脆、风味独特是甜玉米育种中除了丰产性好、抗病性强之外的重要目标[7-10]。目前，一些甜玉米以及超甜玉米品种如金银898[11]、沪甜１号[12]、闽甜6855[13]等品种的培育都朝着改善甜玉米的口感上发展。但是，中国的甜玉米品种与国外特别是美国相比尚有较大的差距，最大的不足是果皮柔嫩度差导致渣多、口感粗糙等问题，这也是现阶段中国甜玉米品质改良所面临的瓶颈。张士龙等[14]采用硬度计法测定了甜玉米籽粒发育过程中的果皮柔嫩度值，发现甜玉米的果皮柔嫩度受到发育时期，种植环境的影响，但果皮柔嫩度的好坏本质上仍决定于材料本身。另有研究表明，玉米果皮柔嫩度与籽粒的果皮厚度呈负相关[15]，口感与果皮粗纤维含量的负相关也达到显著水平[16]。果皮厚度是玉米果皮柔嫩度的主要物理形成因子，而果皮纤维素含量则是玉米果皮柔嫩度的主要化学形成因子[14]。一般情况下，甜玉米的粗纤维含量均较高[17]，因此，降低果皮纤维素含量对于甜玉米品质改良至关重要。浙江省农业科学院旱杂粮研究室的研究也表明甜玉米果皮细胞层数和果皮纤维素含量是影响果皮柔嫩性的关键因素，通过降低果皮细胞层数和果皮纤维素含量可以提高甜玉米的柔嫩性（王美兴等，未发表数据）。因此，甜玉米柔嫩度的改良需要在减少果皮细胞层数的同时降低果皮纤维素的含量，在果皮细胞层数无法减少的情况下，降低果皮纤维素含量则是玉米柔嫩度改良的唯一途径。【本研究切入点】尽管中国已育成了多个甜玉米品种[11-13]，但是，甜玉米柔嫩度品质改良的研究还比较少，甜玉米果皮纤维素含量调控的研究也未见报道，控制此性状的QTL位点和调控基因也尚未定位。【拟解决的关键问题】本研究通过对玉米果皮纤维素含量差异的RILs群体中纤维素含量高、低材料和亲本进行SLAF测序，通过生物信息学分析的方法开发出多态性的SLAF标记，在此基础上，通过BSA的分析方法，定位调控玉米纤维素含量的染色体区段，将此区段的基因映射到拟南芥基因组，通过对比拟南芥基因功能，确定调控玉米果皮纤维素含量的候选基因。为玉米果皮纤维素含量品质改良提供理论依据和基因资源信息。

1 材料与方法

1.1 植物材料及RILs群体的构建

以果皮纤维素含量显著差异的E327和G5-1为亲本，构建重组自交系（RILs）。E327为母本，其纤维素含量为17.62%；G5-1为父本，纤维素含量为11.34%。在138份F6代重组自交系中选择28个纤维素含量高的植株、22个纤维素含量低的植株进行关联分析。

1.2 纤维素含量测定

籽粒浸泡24 h，剥取果皮，将果皮在70℃鼓风干燥箱内烘干，然后完全粉碎。取样品测定水分含量（w）；G4玻璃砂芯漏斗于500℃灼烧至质量恒定；精确称取干燥样品1.00—1.05 g（m0），放入250 mL洁净干燥的锥形瓶中，加入25 mL硝酸-乙醇混合液[18]，装上回流冷凝管，沸水浴加热1 h，用G4玻璃砂芯漏斗抽滤去除溶剂，重复上述操作3—5次，直至纤维变白。用10 mL硝酸-乙醇混合液洗涤残渣，再用热水洗涤至洗涤液用甲基橙试验不呈酸性反应为止，最后用无水乙醇洗涤2次，抽干滤液，将盛有残渣的玻璃砂芯漏斗移入烘箱，于105℃烘干至质量恒定，称质量（m1），然后置于坩埚中500℃灼烧至质量恒定，称取质量（m2）。

1.3 SALF文库构建及高通量测序

对参与混池的各单株取等量的叶片，混合得到高纤维素混池样品（28单株）和低纤维素混池样品（22单株），采用CTAB法提取亲本E327、G5-1和两混池的DNA。对玉米参考基因组进行电子酶切，选择Ⅲ和166Ⅱ酶对DNA样品进行酶切并回收414—464 bp的酶切片段，对回收得到的酶切片段（SLAF标签）进行3′端加A处理、连接Dual-index测序接头[19]、PCR扩增、纯化、混样、切胶选取目的片段，文库质检合格后用Illumina HiSeqTM2500进行测序。

1.4 多态性SLAF标签的鉴定及分布

通过水稻对照数据评估Ⅲ+166Ⅱ的酶切效率，以此判断实验过程的准确性和有效性。根据序列相似性将每个样本测序产生的reads进行聚类，聚类到一起的reads来源于同一个SLAF片段（SLAF标签）。一个SLAF标签在不同样品间序列有差异（包括SNP和InDel），即定义为多态性SLAF标签[20]。统计玉米各染色体上SLAF标签和多态性SLAF标签的分布并作图。

1.5 BAS关联分析

根据多态性SLAF标记的信息，使用GTAK对亲本与子代混池进行SNP开发。在进行BSA分析前，过滤任一混池中read支持度小于3的位点和与相应亲本的SNP类型差异的位点，得到高质量的可信SNP位点，并在此基础上识别两混池间差异的位点。采用ED（算法[21]计算与目的基因连锁的区域。该算法利用混池间差异SNP的频率的差异，计算两混池间的SNP的ED值，ED值越大表明该SNP在两混池间的差异越大。根据连锁的原理，真实的关联区域附近的SNP位点会倾向于在两混池间表现出显著性差异，因此，目标区域内的ED值较大，通过作图，可以在目标关联区域附近观察到较明显的峰。为消除背景噪音，对原始ED值利用LOESS算法[21]进行拟合回归曲线操作，确定显著性阈值，并通过拟合曲线确定最终关联区域。公式如下：

式中，高/高/高/高：在SNP位点上A/T/C/G碱基高纤维素混池测序reads中出现频率；低/低/低/低：在SNP位点上A/T/C/G碱基低纤维素混池测序reads中出现频率。

2 结果

2.1 RILs分离群体的建立

纤维素含量受多基因共同调控。以果皮纤维素含量较高的甜玉米自交系E327和果皮纤维素含量较低的甜玉米自交系G05-1为亲本，构建了纤维素含量变化的RILs群体。在第六代RILs群体中，纤维素含量呈现连续变化，最高的纤维含量为21.3%，而最低为9.3%（电子附表1）。为了开发玉米果皮纤维素含量的分子标记和定位相关调控基因，在RILs群体中挑选28个高纤维含量（17.3%—21.3%）和22个低纤维素含量（9.3%—14.0%）的两组植株，进行BSA（bulked segregate analysis）分析。

2.2 SLAF文库的构建及评价

从JGI网站（https://jgi.doe.gov/）下载玉米参考基因组序列，其大小为2.05 GB，GC含量为46.89%。对玉米的参考基因组进行电子酶切，最终确定选用Ⅲ和166Ⅱ限制性内切酶对进行SLAF文库构建，酶切片段长度在414—464 bp的序列定义为SLAF标签，预测可得到246 254个SLAF标签（电子附表2），这些标签在玉米基因组上分布均匀（电子附表3），能满足基因定位和进一步开发分子标记的要求。

对亲本E327，G05-1和高、低纤维素混池4个样品用Ⅲ和166Ⅱ进行双酶切，Illumina建库，并筛选insert为400—500 bp的片段进行Illumnia HiSeq 2500测序。共获得39.77 M reads数据，测序平均Q30为87.57%，平均GC含量为43.59%。亲本和2个混池的read数和数据质量如表1所示。测序质量评估和碱基分布检查表明NGS（next-generation sequencing）测序获得的数据合格。双端比对效率为 95.76%，符合实验要求。根据Control测序pair-end mapped reads在基因组中的位置计算SLAF标签的实际长度，绘制Control reads插入片段的长度分布图，并估测实际片段选择范围（电子附表3），与预期待相符。通过对玉米NGS数据和Control数据的分析结果表明，所获得的玉米NGS数据符合要求，可以进一步进行SLAF标记开发和基于SLAF标记的BSA分析。

2.3 SLAF标签的确定及其在玉米染色体上的分布

通过SOAP软件[22]将SLAF定位到玉米（B73）的基因组（下载地址：https://www.ncbi.nlm.nih.gov/ assembly/GCF_000005005.1/）上，共获得268 524个SLAF标签，其中父本的测序深度为22.57×，母本的测序深度为22.64×，低纤维素混池（L）平均测序深度为26.86×，高纤维素混池（H）平均测序深度为31.05×（电子附表4）。不同染色体上的SLAF标签在染色体上均匀分布（图1-A）。通过两亲本SLAF标签的对比，找到有序列差异的SLAF标签73 786个，即为多态性SLAF标签。总体上，多态性SLAF标签占总SLAF标签数目的27.48%。SLAF标签及多态性SLAF标签的在玉米染色体上的统计见电子附表5。这些多态性标签在染色体上均匀分布（图1-B），满足作为开发分子标记的条件。

2.4 SNP标记的开发

利用多态性SLAF标签信息，采用GTAK对亲本与子代混池进行SNP标记开发，一共得到SNP标记523 395个。过滤掉低可信度的标记后，在父本E327、母本G05-1、低纤维素混池、高纤维素混池样本中分别得到414 466、411 851、459 632和464 746个SNP标记。这些SNP在4个样本中的杂合度分别位17.37%、16.62%、44.46%和40.05%（表2），杂合度的分布符合亲本及BSA混池的特征。表明基于SLAF标签开发的SNP标记可以用来进行进一步的关联分析。同时，所得到的SNP位点的变异信息（电子附表6），可以用于进一步开发SNP分子标记引物，用于玉米的基因定位或分子标记辅助育种。

图1 SLAF标签（A）和多态性SLAF（B）在染色体上的分布

表1 各样品测序数据统计表

P：父本，M：母本，L和H分别代表2个混池。下同

P: paternal parent line; M: maternal parent line; L and H represent two bulks, similarly hereinafter. The same as below

表2 本研究所检测到的SNP情况

2.5 BSA关联分析

进一步对2个亲本和2个混池的SNP进行分析，为了保证实验的准确性，首先过滤任一混池中read支持度小于3的位点265 983个，然后过滤与相应亲本的SNP类型差异的位点132 740个，最终得到高质量的可信SNP 位点344 441个。其中，280 414位点在2个混池间是一致的，64 027个SNP位点在混间具有差异。

针对这些混池间有差异的SNP位点，使用ED算法[21]计算每个SNP位点的ED值，ED值越大表明该SNP在两混池间的差异越大，差异SNP位点ED值在染色体上的分布如图2所示。为消除背景噪音，对原始ED值利用LOESS算法进行拟合回归曲线操作[21]，取所有位点ED的第998百分位数（0.69）为阈值，关联到6个区段，均位于第5染色体上，标记为Region1-Region6。表3显示Region1-Region6在玉米基因组上的位置及区段长度，这6个区段一共有47个基因，Region1-Region6所对应的基因数分别为1、1、12、5、3和25。

表3 关联区域信息统计表

A：玉米不同染色体上BSA分析得到的ED值；B：玉米第5染色体的ED值。红线显示的是ED分析的阈值

使用BLAST软件将关联区域内的47个基因分别与NR、SwissProt、GO、COG、KEGG数据库比对（电子附表7和电子附表8）。

2.6 候选基因的分析

在这6个关联区域中，Region1和Region2分别只对应一个基因，分别为和，这两个基因可以作为2个候选基因。其他区域分别对应多个基因，但是在每个区域中，可能只有少数基因真正调控纤维素含量变化。因此，需要搜索相关文献对这些基因已经报道的功能进行查阅，进一步缩小候选基因的范围。然而，在玉米中这些基因功能的研究尚未见报道，而模式植物拟南芥的基因功能研究较其他植物更为深入，因为不同植物物种中，同源基因一般具有相同或类似的功能，通过同源比对的方法把这47个基因映射到拟南芥基因组，找到它们对应的拟南芥基因及其注释信息（电子附表8），并查阅拟南芥相关基因的研究工作，推测这些玉米基因的功能。如果某一玉米基因对应的拟南芥同源基因直接或间接影响纤维素的代谢，那么这个玉米基因就确定为可能的候选基因。通过这种方法，最终在Region3、Region4、Region5、Region6分别各确定了2个基因作为候选基因（表4）。

3 讨论

SLAF-seq（specific-locus amplified fragment sequencing）技术是利用二代高通量测序发展而来的一套简化基因组测序技术。通过特定内切酶酶切基因组并筛选特定大小的片段建库测序，从而选择性地得到目标基因组的片段序列。对不同样本SLAF片段进行多样性分析可以得到多态性SLAF标签，用于BSA关联分析。SLAF-seq技术结合BSA关联分析具有通量高、准确性高、成本低、周期短等优势。另外，基于多态性SLAF标签上变异位点，可以开发出大量的SNP标记，用于后续的分子标记辅助育种[20]。目前，基于SLAF-seq的分子标记技术已经在水稻、小麦、大豆、油菜、棉花等作物中得到应用。在玉米中，SLAF标记已成功用于鉴定fea*-9LB030突变体的突变基因[23]。但是，利用SLAF-seq技术开发玉米分子标记尚未见报道。本研究利用果皮纤维素含量差异的两亲本E327、G05-1构建RILs分离群体，并对亲本和纤维素含量高、低的2个混池进行SLAF测序，得到了73 786个在亲本中存在多态性的标签。另外，与参考基因组进行比对，在4个样本中分别得到了414 466、411 851、459 632和464 746个SNP位点，这些位点的信息（电子附表6）可以用于进一步开发SNP分子标记，用于玉米基因定位和分子辅助育种。

表4 各关联区段的候选基因分析

结合BSA分析的方法，本研究一共得到6个与果皮纤维素含量差异性状关联的染色体区段，在这6个区段中一共注释到47个基因。作为模式植物拟南芥，其基因功能的研究远远领先于其他植物。鉴于此，以拟南芥基因组作为参考，根据这47个基因在拟南芥中的同源基因是否直接或间接影响纤维素的代谢来确定候选基因。在关联到的染色体区段1和区段2，虽然没有报道表明其对应的拟南芥基因与纤维素的代谢相关，但由于这两个区段各只有一个基因，因此，这两个基因（和）将作为候选基因进行后续的功能验证分析；第3区段有12个基因，其中在拟南芥中的同源基因是（）。是拟南芥根毛特异表达基因[24]，其与功能冗余，在拟南芥根毛形成中能够促进细胞壁的合成[25]，表明可能间接的影响纤维素的合成。因此，被确定为区段3的候选基因；区段4有5个基因，分别为、、、和。其中和在拟南芥中有共同的同源基因、和，它们同属于GDSL-like Lipase/ Acylhydrolase superfamily protein基因家族。这个基因家族在拟南芥中还没有参与或影响纤维素含量的报道，但最近在棉花中的一项研究表明：棉花中的同源基因位于GhMYB1转录因子的下游，参与棉花中纤维的合成[26]。因此，和被确定此区段的候选基因；染色体区段5有3个基因，其中和均与拟南芥的（）同源。PRX52属于过氧化物酶超级家族蛋白（peroxidase superfamily protein），是一个Ⅲ类过氧化物酶，在细胞壁中表达，调节次生壁的合成[27]。过氧化物酶（peroxidase，POD）可以在细胞中催化超氧化物（O2−-）和过氧化氢（H2O2）的产生，使细胞壁松弛，促进细胞壁的合成[28]。有报道表明，的表达下调，可以抑制细胞的伸长和根的生长[29-30]；在拟南芥突变体中，当植物受到铝胁迫时，纤维素合成基因和的表达显著下调[31]。此外，在玉米中也有研究发现POD的上调表达与细胞壁的合成和节间伸长相关[32]。基于这些研究，第5区段中的2个编码玉米过氧化物酶基因和被确定为候选基因；CesA是关键的纤维素合成酶，桉树中的一项研究表明，、和的启动子上含有W-box，MYC识别为点（MYC recognition site），EEC元件（EEC element）、CCA1-结合位点（CCA1-binding site）以及ARR1-结合位点（ARR1-binding element）。其中，W-box是WRKY转录因子的识别位点[33-37]，ARR1-结合原件是生长素信号途径转录因子ARR结合的位点。被子植物的树干和树枝在重力的刺激下会形成应拉木（tension wood，TW)。杨树的TW中富含纤维素，其TW和常规木（normal wood，NW）的转录组分析发现两者存在大量的差异基因。其中，属于基因家族的（对应拟南芥的WRKY20，基因位点为）在富含纤维素的TW中的表达量是NW的3.47倍，暗示WRKY转录因子与纤维素合成存在着相关性；在同样一组数据中，研究者发现12个生长素应答基因的差异表达[38]，这表明，生长素信号途径也可能间接影响纤维素的合成。也有报道表明拟南芥参与调控细胞壁的形成[39-40]这些研究暗示，纤维素的合成受到WRKY转录因子的调控，同时也受到生长素信号途径的调节。本研究利用BSA分析关联到的染色体区段6中一共有25个基因，其中，编码WRKY转录因子，与拟南芥的是同源基因；在拟南芥中的同源基因是，编码一个类SUAR生长素应答蛋白（SAUR-like auxin-responsive protein）。基于这些研究，将这两个基因确定为候选基因（表4）。

4 结论

通过对玉米果皮纤维素含量差异的亲本和RILs群体混池进行SLAF测序，共获得73 786个SLAF多态性标签，并检测到523 395个SNP位点。获得6个控制玉米果皮纤维素含量的染色体区段，共确定9个候选基因。表明SLAF测序结合BSA分析的方法可以用于快速基因定位的分析。

[1] 葛欣然. 我国玉米产业发展现状及政策调整. 新农业, 2017, 10: 7-8.

GE X R. Current situation and policy adjustment of corn industry in China., 2017, 10: 7-8. (in Chinese)

[2] 丁鑫. 我国玉米产业现状及发展趋势. 现代农村科技, 2017, 10: 10-11.

DING X. Current situation and development trend of corn industry in China.,2017, 10: 10-11. (in Chinese)

[3] 赵华, 王子明. 高产优质甜玉米新品种筛选. 玉米科学, 2017, 25(3): 38-42.

ZHAO H, WANG Z M. Screening of new sweet corn varieties with high yield and good quality.,2017, 25(3): 38-42. (in Chinese)

[4] 石明亮, 薛林, 胡加如, 黄小兰, 陈国清, 陆虎华. 玉米和特用玉米的营养保健作用及加工利用途径. 中国食物与营养, 2011, 17(2): 66-71.

SHI M L, XIE L, HU J R, HUANG X L, CHEN G Q, LU H H. Nutritional and health effects of maize and special corn and its processing and utilization., 2011, 17(2): 66-71. (in Chinese)

[5] Dewanto V, Wu X Z, Liu R H. Processed sweet corn has higher antioxidant activity.,2002, 50(17): 4959-4964.

[6] 岳武. 甜玉米育种中几个问题的探讨. 辽宁农业科学, 2017(1): 62 -63.

YUE W. Discussion on several problems in sweet corn breeding., 2017(1): 62 -63. (in Chinese)

[7] 戴惠学, 施泽平. 优质高产超甜玉米晶甜3号的选育报告. 玉米科学, 2006, 14(2): 64-65.

DAI H X, SHI Z P. Reports on the selecting & breeding of super sweet corn Jingtian No.3 with high yield and good quality.,2006, 14(2): 64-65. (in Chinese)

[8] 姚文华, 韩学莉, 汪燕芬, 谭静, 徐春霞, 陈洪梅, 番兴明. 我国甜玉米育种研究现状与发展对策. 中国农业科技导报, 2011, 13(2): 1-8.

YAO W H, HAN X L, WANG Y F, TAN J,XU C X, CHEN H M, FAN X M. Research status and development strategy for sweet corn breeding in China., 2011, 13(2): 1-8. (in Chinese)

[9] 林建新, 卢和顶, 廖长见. 甜玉米新品种闽甜4号的选育. 福建农业学报, 2012, 27(6): 611-616.

LIN J X, LU H D, LIAO C J. Breeding a new sweet corn variety, Mintian4.,2012, 27(6): 611-616. (in Chinese)

[10] 王蕴波, 杨泉女, 陈文胜, 甄畅迪. 华南地区甜玉米的育种目标初探. 佛山科学技术学院学报(自然科学版), 2010, 28(2): 68-72.

WANG Y B, YANG Q N, CHEN W S, ZHEN C D. Breeding objectives of sweet corn in South China.,2010, 28(2): 68-72. (in Chinese)

[11] 林金元, 王慧, 潘玲玲, 孙大鹏, 于典司, 施标, 卢有林, 郑洪建. 超甜玉米新品种“金银898”的选育. 上海农业学报, 2015, 31(3): 94-97.

Lin J Y, WANG H, PAN L L, SUN D P, YU D S, SHI B, LU Y L, ZHENG H J. Breeding of a new super-sweet corn cultivar ‘Jinyin 898’.2015, 31(3): 94-97. (in Chinese)

[12] 孙萍东, 潘玲玲, 徐声宇, 周美华, 王慧, 孙大鹏, 于典司, 郑洪建, 林金元. 超甜玉米新品种“沪甜1”的选育. 上海农业学报, 2016, 32 (3): 164-167.

SUN P D, PAN L L, XU S Y, ZHOU M H, WANG H, SUN D P, YU D S, ZHENG H J, LIN J Y. Breeding of a new super sweet corn cultivar ‘Hutian 1’., 2016, 32(3): 164-167. (in Chinese)

[13] 林建新, 陈山虎, 廖长见, 卢和顶. 甜玉米新品种“闽甜6855”的选育. 福建农业学报, 2016, 31(11): 1171-1174.

LIN J X, CHEN S H, LIAO C J, LU H D. Breeding a new sweet corn variety, Mintian6855., 2016, 31(11): 1171-1174. (in Chinese)

[14] 张士龙, 贺正华, 黄益勤. 超甜玉米子粒果皮柔嫩度的变化规律. 湖北农业科学, 2016, 55(5): 1105-1108.

ZHANG S L, HE Z H, HUANG Y Q. Variation law of pericarp tenderness of super sweet corn kernel.,2016, 55(5): 1105-1108. (in Chinese)

[15] Schmidt D H, Tracy W F. Effects of starchy sugary and sugary sugary endosperm on pericarp thickness in sweet corn.,1988, 23(5): 885-886.

[16] 乐素菊, 刘厚诚, 张璧, 王晓明. 超甜玉米籽粒乳熟期碳水化合物变化及食用品质. 华南农业大学学报(自然科学版), 2003, 24(2): 9-11.

YUE S J, LIU H C, ZHANG B, WANG X M. Changes of carbohydrate and taste quality in the kernels of super-sweet corn in the milky maturity stage.,2003, 24(2): 9-11. (in Chinese)

[17] 王世恒, 冯凤琴, 徐仁政. 超甜玉米营养品质分析. 玉米科学, 2004, 12(1): 61-62.

WANG S H, FENG F Q, XU R Z. Analysis on nutritional quality of super-sweet corn., 2004, 12(1): 61-62. (in Chinese)

[18] 王林风, 程远超. 硝酸乙醇法测定纤维素含量. 化学研究, 2011, 22(4): 52-55.

WANG L F, CHENG Y C. Determination the content of cellulose by nitric acid-ethanol method., 2011, 22(4): 52-55. (in Chinese)

[19] Kozich J J, Westcott S L, Baxter N T, HIGHLANDER S K, SCHLOSS P D.Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform., 2013, 79(17): 5112-5120.

[20] Sun X, Liu D, Zhang X, LI W, LIU H, HONG W, JIANG C, GUAN N, MA C, ZENG H, XU C, SONG J, HUANG L, WANG C, SHI J, WANG R, ZHENG X, LU C, WANG X, ZHEGN H. SLAF-seq: an efficient method of large-scale De novo SNP discovery and genotyping using high-throughput sequencing.,2013, 8(3): e58700.

[21] Hill J T, Demarest B L, Bisgrove B W, Gorsi B, Su Y C, Yost H J. MMAPPR: Mutation mapping analysis pipeline for pooled RNA-seq., 2013. 23: 687-697.

[22] Li R, Li Y, Kristiansen K, WANG J. SOAP: short oligonucleotide alignment program., 2008, 24(5): 713-714.

[23] Xia C, CHEN L, RONG T, LI R, XIANG Y, WANG P, LIU C, DONG X, LIU B, ZHAO D, WEI R, LAN H. Identification of a new maize inflorescence meristem mutant and association analysis using SLAF-seq method.,2015, 202: 35-44.

[24] Won S K, Lee Y J, Lee H Y, HEO Y K, CHO M, CHO H T. Cis-element-and transcriptome-based screening of root hair-specific genes and their functional characterization in., 2009, 150(3): 1459-1473.

[25] Baumberger N, Steiner M, Ryser U, KELLER B, RINGLI C. Synergistic interaction of the two paralogousgenesand, 2003, 35(1): 71-81.

[26] Yadav V K, Yadav V K, Pant P, singh s p, maurya r, sable a, sawant s v. GhMYB1 regulates SCW stage-specific expression of the GhGDSL promoter in the fibres ofL..,2017, 15(9): 1163-1174.

[27] Passardi F, Cosio C, Penel C, DUNAND C. Peroxidases have more functions than a Swiss army knife., 2005, 24(5): 255-265.

[28] Liszkay A, Kenk B, Schopfer P. Evidence for the involvement of cell wall peroxidase in the generation of hydroxyl radicals mediating extension growth.,2003, 217(4): 658-667.

[29] Van H L, Kuraishi S, Sakurai N. Aluminum-induced rapid root inhibition and changes in cell-wall components of squash seedlings., 1994, 106(3): 971-976.

[30] Goodwin S B, Sutter T R. Microarray analysis ofgenome response to aluminum stress.,2009, 53(1): 85-99.

[31] Fernandezperez F, Pomar F, Pedreno M A, NOVO-UZAL E. The suppression ofaffects fibers but not xylem lignification inby altering the proportion of syringyl units.,2015, 154(3): 395-406.

[32] Zhang Q, Cheetamun R, Dhugga K S, Rafalski J A, Tingey S V, SHIRLEY N J, TAYLOR J, HAYES K, BEATTY M, BACIC A, BURTON R A, FINCHER G B. Spatial gradients in cell wall composition and transcriptional profiles along elongating maize internodes.,2014, 14(1): 27.

[33] Eulgem T, Rushton P J, Robatzek S, SOMSSICH I E. The WRKY superfamily of plant transcription factors.,2000, 5(5): 199-206.

[34] Yu D, Chen C, Chen Z. Evidence for an important role of WRKY DNA binding proteins in the regulation of NPR1 gene expression., 2001, 13(7): 1527-1540.

[35] Dong J, Chen C, Chen Z. Expression profiles of theWRKY gene superfamily during plant defense response.,2003, 51(1): 21-37.

[36] Miao Y, Laun T M, Zimmermann P, ZENTGRAF U. Targets of the WRKY53 transcription factor and its role during leaf senescence in.,2004, 55(6): 853-867.

[37] Xu X, Chen C, Fan B, CHEN Z. Physical and Functional Interactions between Pathogen-InducedWRKY18, WRKY40, and WRKY60 Transcription Factors., 2006, 18(5): 1310-1326.

[38] Lu S, Li L, Yi X, JOSHI C P, CHIANG V L. Differential expression of three eucalyptus secondary cell wall-related cellulose synthase genes in response to tension stress., 2008, 59(3): 681-695.

[39] Wang H, Avci U, Nakashima J, Hahn MG, Chen F, Dixon RA. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants., 2010, 107: 22338-22343.

[40] Yu F, Huaxia Y, Lu W, WU C, CAO X, GUO X. GhWRKY15, a member of the WRKY transcription factor family identified from cotton (L.), is involved in disease resistance and plant development., 2012, 12(1): 144.

附表1 RILs群体纤维素含量及混池信息

Supplemental table 1 The cellulose content of RILs and the bulk information

附表2 电子酶切信息统计

Supplemental table 2 The statistical information of electrical enzyme digestion

附表3 SLAF 标签在各染色体上的数量统计

Supplemental table 3 SLAF number on different chromosome

附表4 SLAF 标签统计

Supplemental table 4 The statistical information of SLAF tags

附表5 SLAF标签和多态性SLAF标签染色体分布统计

Supplemental table 5 The statistic of SLAF tags and polymorphic SLAF tags on different chromosomes

附表6 亲本和高、低纤维素混中的多态性SNP

Supplemental Table 6 The SNPs in parent lines and high/low cellulose content bulks

附表7 关联区域内基因注释统计表

Supplemental table 7 The statistic of annotated genes in associated regions

附表8 关联区域基因的注释

Supplemental table 8 Annotation of associated gene in associated regions

附表9 关联区域玉米基因对应的拟南芥基因的注释

Supplemental table 9 Annotation of Arabidopsis homological genes of Maize genes in associated regions

请根据该文DOI：10.3864/j.issn.0578-1752.2018.08.001 或扫描右侧二维码，从相应网页“附件”处下载查阅以上信息（手机或不兼容）。

SLAF-marker Development and Its Application in BSA Analysis of Cellulose Content in Pericarp of Maize Kernel

Du Longgang, Wang Meixing

(Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021)

【Objective】Reducing cellulose content in pericarp is one of the important goals for quality improvement of sweet maize. However, the research focusing on cellulose content in pericarp of maize is limited and the related gene has not been identified. In order to map the chromosome regions and candidate genes controlling cellulose content in pericarp of sweet maize, specific-locus amplified fragment-sequencing (SLAF-seq) and bulked segregant analysis (BSA) were conducted in this study. 【Method】A recombinant inbred line (RIL) population was constructed using E327 and G5-1 as parents. In F6generation of RILs, the cellulose contents in pericarp of each line were detected. According to cellulose content results, the lines with relatively high and low cellulose contents in pericarp were selected for SLAF tags identification and BSA. For BSA, DNA was extracted from two bulked populations as well as two parent lines and digested withⅢand166Ⅱrestriction enzymes. After digestion, the 414-464 bp DNA fragments were recovered and used for Illumina library construction which were subsequently sequenced by 2nd generation sequencing technology. Based on the polymorphic SLAF tags developed in this analysis, the SNPs-cellulose content association analysis was performed to map the chromosome regions associated with cellulose content in pericarp of sweet maize. Then, the genes located on the associated regions were found. To further reduce the number of candidate genes, the genes inhomological to these genes were identified and annotated. Based on the published works related on the functional analysis of thesegenes, we further identified the candidate genes regulating cellulose content in pericarp of sweet maize. 【Result】The sequencing results of Illumina libraries were consistent with expected. As a result, 73 786 polymorphic SLAF tags were identified, which were uniformly distributed on 10 chromosomes. A total of 523 395 SNPs in those polymorphic tags were identified from the SLAF-sequencing data. Genes responsible for cellulose content in pericarp were mapped onto 6 chromosome regions of maize genome via association analysis and all these 6 chromosome regions were located on the 5th chromosome. In total, 47 gene loci in those regions and 9 genes in these associated regions were identified as candidate genes in further analysis. 【Conclusion】Through SLAF-sequencing based bulking segregated analysis, the cellulose content related genes in pericarp of sweet maize were mapped, suggesting that this method can be applied in gene mapping for other traits.

maize; pericarp; cellulose; specific-locus amplified fragment-sequencing; bulked segregant analysis

（责任编辑 李莉）

10.3864/j.issn.0578-1752.2018.08.001

2017-12-06；

2018-01-04

浙江省农业（粮食）新品种选育重大科技专项（2016C02050-9-2）

杜龙岗，E-mail：51968582@qq.com。

王美兴，E-mail：wangmeixing1975@126.com