潘旭浩 程立銳 陳小翠 蔣彩虹 任民 楊愛國

摘 要：煙草花葉病（CMV）是煙草主要病害之一，篩選與CMV抗病QTL緊密連鎖的分子標(biāo)記是煙草抗病分子育種的重要部分。本研究以抗病材料臺煙8號和感病材料NC82為親本，獲得BC1群體，構(gòu)建抗、感池。采用BSA和SLAF-seq技術(shù)開發(fā)CMV抗病相關(guān)SNP分子標(biāo)記，共獲得180 370個SLAF標(biāo)簽，經(jīng)分析后獲得6156個多態(tài)性SNP標(biāo)記。利用SNP-index分析方法，發(fā)現(xiàn)兩個與CMV抗病高度關(guān)聯(lián)的區(qū)域。利用這些分子標(biāo)記成功定位到了兩個CMV抗病相關(guān)QTL，其中一個QTL定位于Chr01上的55.72 ～60.44 Mb區(qū)間內(nèi)，區(qū)間大小為4.72 Mb；另一個則定位于煙草Chr18上的44.21～48.70 Mb區(qū)間內(nèi)，區(qū)間大小為4.49 Mb。在關(guān)聯(lián)到的QTL區(qū)間內(nèi)設(shè)計分子標(biāo)記可以用于煙草抗CMV育種的輔助選擇。

關(guān)鍵詞：煙草；CMV抗性；QTL；SALF-seq

中圖分類號：S572.03 文章編號：1007-5119（2018）05-0001-08 DOI：10.13496/j.issn.1007-5119.2018.05.001

Abstract： Cucumber mosaic virus （CMV） is one of most important diseases in tobacco， and selection of molecular markers closely linked to CMV resistant QTLs is a key for tobacco CMV-resistance breeding. In order to identify genes closely linked to CMV resistance and develop resistant varieties， a BC1 population derived from resistant-parent Taiyan-8 and susceptible-parent NC82 was used to construct resistant and susceptible pools. By application of BSA and SLAF-seq technique， a total of 180，370 SALF labels were obtained and 6，156 polymorphic SNP markers were developed. As a result， two CMV-associated regions were identified by analysis of SNP-index data. Our results showed that abundant polymorphic SNPs can be developed between tobacco varieties using SALF-seq. By utilization of these SNPs， two CMV-resistance QTLs were successfully identified in Chr01 and Chr18， with 4.72 Mb and 4.49 Mb interval， respectively. The two identified QTLs are useful for in molecular marker assisted selection in tobacco breeding programs for CMV resistance.

Keywords： tobacco； CMV resistance； QTL； SLAF-seq

黃瓜花葉病毒?。–ucumber mosaic virus，CMV）是一種病毒性病害，在煙草中主要經(jīng)由蚜蟲或汁液摩擦進行傳播，感病后主要表現(xiàn)為煙葉花斑、葉片卷曲，嚴重時直接影響煙株發(fā)育，導(dǎo)致煙葉產(chǎn)量、質(zhì)量下降。其流行速度快，發(fā)生范圍廣，是我國發(fā)生最為普遍、危害最大的煙草病害之一。目前生產(chǎn)上對煙草CMV防治主要還是以化學(xué)防治為主，不僅浪費人力，增加種植成本，而且對于環(huán)境也會造成不利的影響，因此，培育CMV抗病品種一直被認為是防治CMV的有效手段。但是由于煙草CMV抗性由多基因控制[1]，高抗病的種質(zhì)資源缺乏，通過常規(guī)育種手段選育抗病品種的耗時長，效果差。

因此開發(fā)與抗病基因緊密連鎖的分子標(biāo)記，通過分子標(biāo)記輔助育種可以加快培育煙草抗CMV品種的進程，提高育種效率。

前人關(guān)于CMV抗病的研究多集中在辣椒[2-3]、黃瓜[4-5]等作物上。由于煙草遺傳背景狹窄[6]，可利用的傳統(tǒng)分子標(biāo)記數(shù)量少，多態(tài)性差，關(guān)于煙草CMV抗病QTL（Quantitative trait loci，數(shù)量性狀基因座）的研究較為滯后。目前，僅有少量與煙草CMV抗病QTL相關(guān)的文獻發(fā)表，范靜苑等[7]構(gòu)建了煙草BC1群體，并基于混合群體分組分析法（BSA）對CMV抗病基因進行了初步定位，獲得了兩個和CMV抗性基因連鎖的標(biāo)記SM1350和SM2270；文柯等[1]在煙草上對BC1群體進行QTL基因定位，找到了一個與抗性QTL緊密連鎖的標(biāo)記53 735，遺傳距離為0.1 cM，該QTL位于10號連鎖群，貢獻率為13%。

目前發(fā)表的煙草抗CMV分子標(biāo)記無論從數(shù)量還是準(zhǔn)確性上，離實際應(yīng)用分子標(biāo)記進行輔助育種仍有很大距離。簡化基因組測序（Specific Length Amplified Fragment Sequencing，SLAF-seq）技術(shù)是以生物信息學(xué)為基礎(chǔ)的高度自動化高通量測序技術(shù)[8]，其具有可重復(fù)性高，測序時間短，獲取信息量大，可利用多態(tài)性SNP標(biāo)記多等重要特點，可以有效克服煙草遺傳背景狹窄，傳統(tǒng)多態(tài)性標(biāo)記缺乏的缺點，利用該技術(shù)有望對煙草CMV抗病基因進行精細定位。

本研究在前人研究基礎(chǔ)上，基于390個BC1單株接種CMV后鑒定結(jié)果，構(gòu)建抗病池和感病池，利用SLAF-seq開發(fā)大量SNP多態(tài)性標(biāo)記，獲得CMV主效QTL位點，定位了這些位點在煙草染色體上的位置，為煙草抗CMV分子標(biāo)記開發(fā)以及抗病育種提供理論依據(jù)。

1 材料與方法

1.1 試驗材料

供試親本及群體：抗病親本為臺煙8號（P1），感病親本為NC82（P2），兩個親本雜交后F1代與臺煙8號回交，獲得包含390個單株BC1群體作為供試群體。

供試材料無性擴繁：利用P1、P2、F1、BC1群體材料單株的葉片進行組織擴繁、無性繁殖，每個單株均獲得與之基因組完全一樣的家系材料，這些家系材料用于接種CMV，調(diào)查病情，家系的病情即代表擴繁前的單株病情。

毒源：用于接種鑒定煙草CMV的毒源為中國農(nóng)業(yè)科學(xué)院煙草研究所病蟲害綜合防治中心提供的CMV-C（普通株系），在使用前15 d進行復(fù)壯以防止保存期病毒致病性發(fā)生變異。

1.2 群體接種、病情指數(shù)鑒定及群體病級方差分析

對供試BC1群體CMV抗性鑒定在中國農(nóng)業(yè)科學(xué)院煙草研究所青島即墨基地溫室進行。取經(jīng)過復(fù)壯的攜帶CMV病毒的枯斑三生煙葉于研缽，加少許石英砂和磷酸緩沖液，研磨充分，雙層紗布過濾，然后按終體積比1：5加入磷酸緩沖液定容。在煙苗真葉期時，選擇充分展開的第4～5片真葉，在葉面上均勻撒400～500目的石英砂，用脫脂棉球蘸取病毒水溶液輕輕摩擦1～2遍進行接種，盡量保證每片葉接種力度、深度的一致性，接種后立即用水沖洗葉面。在第1次接種5 d后進行第2次摩擦接種，兩次接種力度保持一致。

分別在第2次接種后14～20 d以及20～27 d進行CMV病害調(diào)查，并計算病情指數(shù)。依據(jù)煙草病蟲害分級及調(diào)查方法GB/T 23222—2008進行，病毒病分級標(biāo)準(zhǔn)，以株為單位，分0級，1級，3級，5級，7級，9級共六個級別，分別為：

0級：全株無病。

1級：心葉明脈或輕微花葉，病株無明顯矮化。

3級：1/3葉片花葉但不變形，或病株矮化為正常株高的3/4以上。

5級：1/3至1/2葉片花葉，或少數(shù)葉片變性，或主脈壞死，植株矮化為正常株高的2/3至3/4。

7級：1/2至2/3葉片花葉，或變形，或主側(cè)脈壞死，或植株矮化為正常株高的1/2到2/3。

9級：全株葉片花葉，嚴重變形或壞死，病株矮化為正常株高的1/2以上。

病情指數(shù)=100×∑（各級病葉數(shù)×各級代表值）/（調(diào)查總?cè)~數(shù)×最高級代表值）。群體抗性分類標(biāo)準(zhǔn)：免疫（病情指數(shù)等于0，無侵染）；高抗（0<病情指數(shù)≤15）；中抗（15<病情指數(shù)≤30）；中感（30<病情指數(shù)≤50）；高感（50<病情指數(shù)≤100）[9]。

利用SPSS軟件對BC1群體病級進行單因素ANOVA分析，以確定群體內(nèi)單株抗病性差異是否顯著。

1.3 群體DNA提取及抗感基因池構(gòu)建

對所有材料臺煙8號、NC82、F1及BC1群體用植物基因組DNA提取試劑盒進行單株提取，為了保證提取的濃度和純度，本法略微有所改進：根據(jù)人工CMV接種及鑒定的結(jié)果，從BC1分離群體中選取極端抗病株與極端感病株各52個建立抗感基因池，每個家系樣品根據(jù)所測的DNA原始濃度進行等質(zhì)量混樣，保證所加入的DNA樣品總量相等，抗池標(biāo)為R-bulk，感病標(biāo)為S-bulk。

1.4 基于SLAF-seq技術(shù)的QTL定位

根據(jù)酶切預(yù)測軟件（SLAF_Predict）進行煙草基因組酶切預(yù)測，本研究采用SspI和HaeIII進行酶切消化；37 ℃酶切15 h，反應(yīng)結(jié)束后用QIAGEN試劑盒純化，50 ?L EB 回溶；在20 ℃反應(yīng)體系下，加入總體系為100 ?L的酶切DNA片段、ddH2O、dNTPs mix、T4 DNA polymerase、Klenow DNA polymerase、T4 PNK，反應(yīng)30 min進行5'末端磷酸化修飾；隨后在20 ℃反應(yīng)條件下，加入Klenow 片段和dATP，在酶切片段的3'端加上單核苷酸A，利用T4連接酶將Duplex測序接頭連接在3'端加單核苷酸A的酶切片段上，反應(yīng)結(jié)束后QIAGEN試劑盒純化，10 ?L EB回溶；在20 ℃反應(yīng)條件下，將上一步純化產(chǎn)物同時加入DNA ligase buffer 2X、Adapter oligo mix以及DNA ligase反應(yīng)15 min，反應(yīng)結(jié)束后用QIAGEN試劑盒純化，30 ?L EB回溶；上一步純化的產(chǎn)物用2%瓊脂糖凝膠電泳60 min，電壓為120 V； 切取464～494 bp目的片段利用QIAGEN試劑盒純化后將產(chǎn)物稀釋，利用Illumina high-throughput sequencing platform（Illumina， Inc； San Diego， CA， U S）進行雙端測序。通過評估Control日本晴測序數(shù)據(jù)監(jiān)控實驗過程是否正常，并對序列的reads雙端比對效率、酶切效率、CG含量以及Q20進行分析，評估測序質(zhì)量。

基因組經(jīng)酶切處理打斷為無數(shù)小片段，每個片段相當(dāng)于一個標(biāo)記位點，同一位點reads序列通過相似性聚類，形成一個group。一個group中一般存在1～4條高深度片段，其余全為低深度片段。對測序得到的各樣品的reads數(shù)據(jù)進行評估，通過相似性聚類、基因型糾錯后，再定位到基因組上，得到SLAF標(biāo)簽，對SLAF標(biāo)簽進行多態(tài)性分析，得到樣品間的多態(tài)性SNP標(biāo)記。根據(jù)得到的多態(tài)性標(biāo)記以及抗池、感池的CMV表型鑒定結(jié)果進行關(guān)聯(lián)分析，將CMV主效QTL位點定位在煙草染色體上。

2 結(jié) 果

2.1 供試群體無性繁殖

經(jīng)過無性繁殖培養(yǎng)得到無菌苗臺煙8號40棵、NC82 40棵、F1 40棵以及BC1群體800棵。在無性擴繁過程中，排除掉家系內(nèi)新苗大小不一致，或因污染、分化數(shù)不夠等情況，共篩選出煙苗長勢一致，家系內(nèi)植株數(shù)5株以上的BC1家系452個，共7362個單株，平均每個BC1家系16株。

2.2 BC1無性繁殖家系CMV抗病鑒定與遺傳分析

在植株5～6葉期時進行兩次幼苗摩擦接毒，在接種后第20～25天進行抗病鑒定，結(jié)果發(fā)現(xiàn)（圖1）：抗性親本臺煙8號，接種CMV后無明顯感病癥狀；感病親本NC82在完全感病后表現(xiàn)為葉面花葉、皰斑嚴重，葉片畸形呈柳條狀、鐮刀狀，葉緣底部橢圓而葉尖細長，植株矮化、纖弱；F1大部分表現(xiàn)出花葉、畸形等感病癥狀（圖1）。BC1群體材料則表現(xiàn)出分離，部分植株表現(xiàn)為抗病，部分植株表現(xiàn)為感病，而多數(shù)植株發(fā)病表現(xiàn)介于兩個親本之間（圖2）。為減小CMV抗病鑒定中的人為、環(huán)境誤差，BC1群體每個單株均構(gòu)建了無性繁殖家系，通過群體單株家系的平均抗、感病表現(xiàn)，計算群體單株的病情指數(shù)，最終對390個BC1群體家系進行了鑒定分析，結(jié)果表明，臺煙8號病情指數(shù)介于10～15之間；NC82病情指數(shù)在30～50。BC1群體中，約78個家系病情指數(shù)小于15，鑒定為抗病，206個家系病情指數(shù)在15～30，鑒定為中抗，約106個家系病情指數(shù)大于30，鑒定為感病。BC1群體CMV抗病性狀整體呈現(xiàn)連續(xù)分布，因此說明抗病親本的CMV抗性由多基因控制，為數(shù)量遺傳性狀（圖2）。BC1群體病級方差分析結(jié)果表明（表1），群體內(nèi)單株間病級指數(shù)存在極顯著差異。

2.3 BC1群體抗、感池構(gòu)建

根據(jù)BC1群體獲得CMV病情指數(shù)進行抗、感

池構(gòu)建：選取的原則是單株家系數(shù)量大于等于5，且抗池病情指數(shù)小于15，感池病情指數(shù)大于33.33。根據(jù)單株DNA的檢測濃度進行等質(zhì)量混樣，混樣數(shù)為抗、感池各52個樣品。經(jīng)檢測P1、P2、R-bulk和S-bulk 的濃度均大于50 ng/?L，OD260/280在1.8～2.0之間，OD260/230大于2.0，表現(xiàn)單峰，無色素和蛋白質(zhì)污染，檢測結(jié)果合格。瓊脂糖電泳檢測DNA完整性高（圖3）。

2.4 酶切效率統(tǒng)計與建庫評估

序列的reads雙端比對結(jié)果顯示，Control日本晴水稻在本次實驗雙端比對效率為74.94%，比對效率正常；Control數(shù)據(jù)的酶切效率E-ratio為96.54%，酶切反應(yīng)良好[10]；根據(jù)酶切體系，通過Control評估所有測序reads的插入片段，繪制所有reads的插入片段分布圖，并估測實際的切膠范圍，本實驗Control測序數(shù)據(jù)的插入片段大部分都分布在450～550 bp，插入片段的中心值為482 bp（圖4）。綜上所述，文庫構(gòu)建正常。

2.5 測序數(shù)據(jù)統(tǒng)計與評估

對已構(gòu)建的BC1群體抗、感池以及兩個親本進行SLAF-seq全基因組測序，采用Q20檢測和堿基分布檢查的方法對測序結(jié)果進行評估和過濾以保證后期的分析質(zhì)量。測序結(jié)果顯示，試驗總共得到54.4 Mb測序reads，總數(shù)據(jù)量為8.7 Gb；其中臺煙8號作為抗性親本有效reads 13 854 896個，NC82作為感性親本有效reads為9 139 744個，感池S-bulk有效reads為16 321 266個，抗池R-bulk有效reads為15 023 916個。試驗的堿基質(zhì)量主要集中在20～40之間，即測序錯誤率在0.01%～1%，能夠保證后面的多態(tài)性分析和關(guān)聯(lián)分析。測序獲得平均GC 含量為32.93%（表2），GC含量均較低，說明達到測序要求。

2.6 SNP標(biāo)記開發(fā)及抗CMV主效QTL位點定位

本試驗獲得有效SLAF標(biāo)簽數(shù)為180 370個。對得到的180 370個SLAF標(biāo)簽，根據(jù)等位基因數(shù)和基因序列之間的差異進行多態(tài)性分析，共得到SNP、INDEL、No Polymorphism和Repeat 4種SLAF標(biāo)簽，其中SNP和INDEL統(tǒng)稱為多態(tài)性標(biāo)記。本次分析共得到6156個多態(tài)性標(biāo)記且均為SNP標(biāo)記，多態(tài)性比例為3.41%（表3）。

SNP-index方法計算關(guān)聯(lián)值，并采用SNPNUM方法對△SNP-index 進行擬合，取每個SNP附近200個SNP 的△SNP-index的中值作為該位點擬合后的關(guān)聯(lián)值。根據(jù)計算機模擬實驗計算[11]結(jié)果，發(fā)現(xiàn)與CMV抗性性狀關(guān)聯(lián)的SNP標(biāo)記主要集中在兩個區(qū)間之內(nèi)：其中一個QTL定位于煙草Chr01上的55.72～60.44 Mb區(qū)間內(nèi)，區(qū)間大小為4.72 Mb；另一個則定位于煙草Chr18上的44.21～48.70 Mb區(qū)間內(nèi)，區(qū)間大小為4.49 Mb。

利用茄科基因組數(shù)據(jù)庫（https：//solgenomics. net/organism/Nicotiana_tabacum/genome）對兩個區(qū)間進行基因預(yù)測。在Chr01的4.72 Mb區(qū)間內(nèi)有34個基因，其中有20個基因功能已被注釋；在Chr18的4.49 Mb區(qū)間內(nèi)有41個基因，其中有31個基因功能已被注釋。經(jīng)初步分析，兩個區(qū)間內(nèi)均存在鋅指蛋白家族基因，該基因家族成員可能與煙草CMV抗性相關(guān)。

3 討 論

分子標(biāo)記輔助育種技術(shù)可以大大提高作物育種效率[12]。在早期的研究中，SSR標(biāo)記因其穩(wěn)定、可靠、易操作的特點，成為在連鎖作圖、QTL定位中應(yīng)用最為廣泛的標(biāo)記[13]。但是，由于煙草本身遺傳背景狹窄，標(biāo)記多態(tài)性差，使得SSR標(biāo)記在煙草抗病基因QTL定位研究中利用難度大、效率低、效果差。目前，SLAF-seq技術(shù)在單體型圖譜繪制、遺傳圖譜繪制、關(guān)聯(lián)性圖譜繪制和多態(tài)性圖譜繪制等方面的應(yīng)用已有大量成功案例，利用SLAF技術(shù)和BSA進行SNP開發(fā)及QTL定位的研究已在擬南芥[14]、黃瓜[15]、甘藍型油菜[16]等不同植物上有相關(guān)報道。利用SLAF技術(shù)可以在短時間內(nèi)獲得遍布整個基因組的海量多態(tài)性SNP標(biāo)記，從而實現(xiàn)對目標(biāo)QTL的精細定位。

本研究結(jié)合BSA方法，構(gòu)建了BC1群體以及其群體單株無性繁殖家系，通過對BC1群體單株家系CMV抗病性的鑒定獲得了較為準(zhǔn)確的抗病表型數(shù)據(jù)，并構(gòu)建抗、感病極端混池。利用SLAF-seq技術(shù)獲得了覆蓋全基因組SLAF標(biāo)簽數(shù)為180 370個，經(jīng)篩選其中多態(tài)性SNP標(biāo)記共6156個，多態(tài)性比率為3.41%。

與SSR標(biāo)記進行比較，牟建英[17]以白粉病抗感材料為親本篩選了2317對SSR引物，得到61對多態(tài)性標(biāo)記，多態(tài)性比率僅為2.6%；高亭亭等[18]利用2653個SSR標(biāo)記對煙草黑脛病抗感親本進行全基因組篩選，共得到184對多態(tài)性引物，多態(tài)比率為6.9%；WU等[19]利用1376對SSR引物對白稈突變體及HD之間進行全基因組篩選，最終僅獲得96個在全基因組均勻分布的SSR多態(tài)性標(biāo)記，多態(tài)比率為6.9%；而目前密度最大的烤煙SSR標(biāo)記遺傳連鎖圖僅含覆蓋24個連鎖群的611個SSR標(biāo)記[20]。從多態(tài)性比率來看，本研究中SNP標(biāo)記多態(tài)性比率與SSR標(biāo)記差別不大，這進一步驗證了煙草遺傳背景狹窄、標(biāo)記多態(tài)性差的結(jié)論。但是由于開發(fā)的SNP標(biāo)記基數(shù)大，覆蓋密度高，與前人研究進行對比，本研究中通過SLAF-seq技術(shù)獲得的煙草多態(tài)性SNP標(biāo)記數(shù)量是常規(guī)SSR標(biāo)記的幾十倍甚至上百倍，有效地克服了SSR標(biāo)記在煙草研究上的缺點，為煙草CMV抗病QTL定位提供了有力的技術(shù)保證。

一般來說，在構(gòu)建連鎖圖譜進行基因定位的研究中，標(biāo)記覆蓋的連鎖群越全，連鎖群上標(biāo)記密度越大，定位效果越好，反之則定位QTL可靠性、可用性就會降低。在煙草CMV抗病QTL研究中，文柯[1]利用SSR標(biāo)記對煙草CMV抗性進行定位，在該研究中63個標(biāo)記僅覆蓋了11個連鎖群，其中標(biāo)記最多的連鎖群僅有8個標(biāo)記，而有5個連鎖群僅由2個標(biāo)記構(gòu)成；在煙草其他病害QTL定位中，TONG等[21]利用2000多對SSR標(biāo)記對煙草赤星病進行定位，從中篩選出了覆蓋24個連鎖群的181個多態(tài)性標(biāo)記，平均每個連鎖群標(biāo)記數(shù)為7.54個；LAN等[22]在對煙草青枯病的研究中，利用201個SSR標(biāo)記構(gòu)建了24個鎖鏈群，平均每個連鎖群標(biāo)記數(shù)8.37個。而在本研究中，利用分布于全部24個連鎖群的6156個SNP多態(tài)性用于CMV抗病QTL定位，平均每條染色體多態(tài)性標(biāo)記達到256.5個。

從CMV抗病QTL定位結(jié)果來看，文柯利用SSR標(biāo)記找到了一個與CMV抗病相關(guān)的QTL并將其定位在10號連鎖群上，10號連鎖群僅含有2個SSR標(biāo)記53 735和54 169，標(biāo)記之間間距為13.41 cM；在本研究中，我們發(fā)現(xiàn)了兩個CMV抗病QTL，其中一個QTL被定位于Chr01染色體的4.72 Mb區(qū)間內(nèi)，另一個被定位于Chr18染色體的4.49 Mb區(qū)間內(nèi)。本研究相較前人研究定位了兩個抗CMV主效QTL，這也與煙草CMV受多基因控制的結(jié)論更吻合，同時這兩個QTL定位區(qū)間也更小，為進一步基因精細定位與克隆打好了基礎(chǔ)。本研究結(jié)果充分說明在煙草上利用SLAF-seq技術(shù)開發(fā)SNP標(biāo)記進行基因定位比SSR標(biāo)記更為可靠、高效。

本研究進一步對Chr01和Chr18染色體上定位QTL區(qū)間內(nèi)包含的基因進行分析，發(fā)現(xiàn)在Chr01的4.72 Mb區(qū)間內(nèi)共有34個基因，發(fā)現(xiàn)其中有14個基因功能未知，其余20個基因功能已被注釋；在Chr18的4.49 Mb區(qū)間內(nèi)共有41個基因，其中有10個基因功能未知，另外31個基因功能已被注釋。對已注釋基因的功能進行分析發(fā)現(xiàn)，在兩個區(qū)間內(nèi)均發(fā)現(xiàn)鋅指蛋白（Zinc finger）家族基因，其中在Chr01區(qū)間內(nèi)的有兩個鋅指蛋白家族基因，它們分別為C3HC4類型和C2C2類型；而在Chr18區(qū)間內(nèi)也有兩個鋅指蛋白基因，均為C2H2類型。鋅指蛋白是一個龐大的轉(zhuǎn)錄因子家族，其功能主要是調(diào)控基因表達[23]，根據(jù)鋅指蛋白保守基序特點，可以將其分成9大類，而C2H2、C3HC4、C2C2便是其中主要的3種類型[24]。前人關(guān)于鋅指蛋白的研究主要集中在植物發(fā)育和非生物脅迫上，而抗病方面的研究目前還處于起步階段，但是根據(jù)前人有限的報道，這3種類型的鋅指蛋白均有家族成員與植物抗病功能相關(guān)。其中，C2H2和C3HC4家族基因中已有其與CMV抗性相關(guān)的報道，如煙草ZFT1基因在感染煙草花葉病毒后也會受到不同程度地誘導(dǎo)，激起敏感應(yīng)答[25]；而擬南芥STZ基因感染黃瓜花葉病毒，同樣誘導(dǎo)過敏反應(yīng)[26]；水稻鋅指蛋白基因OsBIRF1過表達載體轉(zhuǎn)入煙草后，轉(zhuǎn)基因煙草增強了對煙草花葉病毒、野火病菌和煙粉虱的抗性，防御相關(guān)基因的表達水平也升高了[27]；除此之外C2C2類型基因如OsLOL2有報道稱過表達該基因可提高水稻對稻瘟病的抗病性[28]。因此，本研究初步預(yù)測上述4個鋅指蛋白基因可能就是參與煙草CMV抗病調(diào)控的候選基因。在本研究基礎(chǔ)上，進一步開展抗病基因克隆及基因功能驗證并開發(fā)分子標(biāo)記，對于闡明煙草與CMV病毒間的互作機制和CMV抗病分子標(biāo)記輔助育種具有重要意義。

4 結(jié) 論

利用SLAF-seq技術(shù)對CMV抗、感親本及后代群體進行簡化測序，開發(fā)了覆蓋24個連鎖群的6156個多態(tài)性標(biāo)記，對BC1 群體抗病表型及基因型進行連鎖分析，找到了兩個與CMV抗性緊密連鎖的QTL，分別定位于染色體Chr01的4.72 Mb以及Chr18的4.49 Mb區(qū)間內(nèi)，本研究獲得的兩個抗病QTL為煙草CMV抗病分子育種提供了重要的理論依據(jù)。

參考文獻

[1]文軻，張志明，任民，等. 烤煙CMV抗性基因QTL定位[J]. 中國煙草科學(xué)，2013，34（3）：55-59.

WEN K， ZHANG Z M， REN M et al. QTL Analysis of the resistance gene to cmv in flue-cured tobacco[J]. Chinese Tobacco Science， 2013， 34（3）： 55-59.

[2]YAO M， LI N， WANG F， et al. Genetic analysis and identification of QTLs for resistance to cucumber mosaic virus in chili pepper （Capsicum annuum L.）[J]. Euphytica， 2013， 193（2）： 135-145.

[3]趙娟，王立浩，毛勝利，等. 辣椒抗黃瓜花葉病毒QTL分析[J]. 中國蔬菜，2011，1（18）：30-35.

ZHAO J， WANG L H， MAO S L， et al. QTL Analysis of cucumber mosaic virus resistance in pepper （capsicum annuum L.）[J]. China Vegetables， 2011， 1（18）： 30-35.

[4]張桂華，黃煥煥，韓毅科，等. 與黃瓜感花葉病毒病基因連鎖的分子標(biāo)記[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報，2010，18（2）：280-284.

ZHANG G H， HUANG H H， HAN Y K et al. Molecular marker linked to cucumber mosaic virus resistant gene in cucumber[J]. Journal of Agricultural Biotechnology， 2010， 18（2）： 280-284.

[5]黃煥煥. 黃瓜CMV抗性相關(guān)基因的分子標(biāo)記研究[D]. 咸陽：西北農(nóng)林科技大學(xué)，2007.

HUANG H H. Study on molecular marker of cucumber mosaic virus resistant-related gene in cucumber （Cucumis Sativus L.）[D]. Xianyang： Northwest Agriculture and Forest University， 2007.

[6]楊本超，肖炳光，陳學(xué)軍，等. 基于ISSR標(biāo)記的烤煙種質(zhì)遺傳多樣性研究[J]. 遺傳，2005，27（5）：753-758.

YANG B C， XIAO B G， CHEN X J， et al. Genetic fiversity of flue-cured tobacco varieties based on ISSR markers[J]. HEREDITAS， 2005， 27（5）： 753-758.

[7]范靜苑，王元英，蔣彩虹，等. 煙草CMV抗性鑒定及抗性基因的SSR標(biāo)記研究[J]. 分子植物育種，2009，7（2）：355-359.

FAN J Y， WANG Y Y， JIANG C H， et al. Resistance identification and SSR marking of resistance gene of tobacco cucumber mosaic vanirus[J]. Molecular Plant Breeding， 2009， 7（2）： 355-359.

[8]SUN X， LIU D， ZHANG X， et al. SLAF-seq： an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing[J]. Plos One， 2013， 8（3）： e58700[2018-05-02]. https：//www.ncbi. nlm.nih.gov/pmc/articles/PMC3602454.

[9]裘維蕃，梁訓(xùn)生. 調(diào)查煙草花葉病病情指數(shù)的計算方法[J]. 植物保護，1963，1（3）：142.

QIU W F， LIANG X S. The calculation method for disease index investigation of tobacco mosaic disease [J]. Plant Protection， 1963， 1（3）： 142.

[10]INITIATIVE T A G. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana[J]. Nature， 2000， 408（6814）： 796-815.

[11]TAKAGI H， ABE A， YOSHIDA K， et al. QTL-seq： rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations[J]. Plant Journal for Cell & Molecular Biology， 2013， 74（1）： 174-183.

[12]VOS P， HOGERS R， BLEEKER M， et al. AFLP： a new technique for DNA fingerprinting[J]. Nucleic Acids

Research， 1995， 23（21）： 4407.

[13]CAVAGNARO P F， SENALIK D A， YANG L， et al. Genome-wide characterization of simple sequence repeats in cucumber （Cucumis sativus L.）[J]. Bmc Genomics， 2010， 11（1）： 569[2018-05-02]. https：//bmcgenomics. biomedcentral.com/articles/10.1186/1471-2164-11-569.

[14]MENG Y， HUANG Y， WANG Q， et al. Phenotypic and genetic characterization of resistance in Arabidopsis thaliana to the oomycete pathogen Phytophthora parasitica[J]. Frontiers in Plant Science， 2015（6）： 378[2018-05-02]. https：//www.ncbi.nlm.nih.gov/pmc/ articles/PMC4445315/.

[15]XU X， LU L， ZHU B， et al. QTL mapping of cucumber fruit flesh thickness by SLAF-seq[J]. Scientific Reports， 2015， 5： 15829[2018-05-02]. https：//www.ncbi.nlm.nih. gov/pmc/articles/PMC4623748.

[16]王偉，劉凡，任莉，等. 采用SLAF-seq技術(shù)開發(fā)甘藍型油菜霜霉病抗性SNP位點[J]. 中國油料作物學(xué)報，2016，38（5）：555-562.

WANG W， LIU F， REN L， et al. Resistance SNP development to downy mildew in Brassica napus using SLAF-seq technique[J]. Chinese Journal of Oil Crop Sciences， 2016， 38（5）： 555-562.

[17]牟建英，錢玉梅，任民等. 煙草白粉病抗性基因的QTL定位[J]. 中國煙草學(xué)報，2013（4）：105-108.

MU J Y， QIAN Y M， REN M， et al. QTL analysis of resistance gene to powdery mildew in tobacco [J]. Chinese Tobacco Science， 2013（4）： 105-108.

[18]高亭亭，蔣彩虹，羅成剛，等. 煙草品種Beinhart1000-1抗黑脛病基因的QTL定位[J]. 植物遺傳資源學(xué)報，2015，16（2）：330-335.

GAO T T， JIANG C H， LUO C G， et al. Mapping of quantitative trait loci affecting resistance to black shank in tobacco line Beinhart1000-1[J]. Journal of Plant Genetic Resources， 2015， 16（2）： 330-335.

[19]WU Q， WU X， ZHANG X， et al. Mapping of two white stem genes in tetraploid common tobacco （Nicotiana tabacum L.）[J]. Molecular Breeding， 2014， 34（3）： 1065-1074.

[20]肖炳光. 煙草基因組計劃進展篇：3. 煙草分子標(biāo)記遺傳連鎖圖構(gòu)建和重要抗病基因定位[J]. 中國煙草科學(xué)，2013，34（3）：118-119.

XIAO B H. Progress in the tobacco genome project： 3. Construction of linkage map using molecular markers and mapping of important resistance genes for tobacco[J]. Chinese Tobacco Science， 2013， 34 （3）： 118-119.

[21]TONG Z， JIAO T， WANG F， et al. Mapping of quantitative trait loci conferring resistance to brown spot in flue-cured tobacco （Nicotiana tabacum L.）[J]. Plant Breeding， 2012， 131（2）： 335-339.

[22]LAN T， ZHENG S， YANG L， et al. Mapping of quantitative trait loci conferring resistance to bacterial wilt in tobacco （Nicotiana tabacum L.）[J]. Plant Breeding， 2015， 133（5）： 672-677.

[23]黃驥，王建飛，張紅生. 植物C2H2型鋅指蛋白的結(jié)構(gòu)與功能[J]. 遺傳，2004，26（3）：414-418.

HUANG J， WANG J F， ZHANG H S. Structure and Function of Plant C2H2 Zinc finger protein[J]. HEREDITAS （Beijing）， 2004， 26（3）： 414-418.

[24]TAKATSUJI H. Zinc-finger transcription factors in plants[J]. Cellular & Molecular Life Sciences Cmls， 1998， 54（6）： 582-596.

[25]UEHARA Y， TAKAHASHI Y， BERBERICH T， et al. Tobacco ZFT1， a transcriptional repressor with a Cys2/His2 type zinc finger motif that functions in spermine-signaling pathway[J]. Plant Molecular Biology， 2005， 59（3）： 435-448.

[26]MARATHE R， GUAN Z， ANANDALAKSHMI R， et al. Study of Arabidopsis thaliana resistome in response to cucumber mosaic virus infection using whole genome microarray[J]. Plant Molecular Biology， 2004， 55（4）： 501-520.

[27]LIU H， ZHANG H， YANG Y， et al. Functional analysis reveals pleiotropic effects of rice RING-H2 finger protein gene OsBIRF1， on regulation of growth and defense responses against abiotic and biotic stresses[J]. Plant Molecular Biology， 2008， 68（1-2）： 17-30.

[28]XU C， HE C. The rice OsLOL2， gene encodes a zinc finger protein involved in rice growth and disease resistance[J]. Molecular Genetics & Genomics， 2007， 278（1）： 85-94.