周仙莉　滕長才　張紅巖　林夕　劉玉玲　吳小燕　侯萬偉　劉玉皎

摘要： 作為一種糧食、蔬菜、飼料、綠肥兼用的食用豆類，蠶豆在人類健康、土壤改良等方面發(fā)揮了重要作用。高密度遺傳圖譜的構(gòu)建是蠶豆基因組研究的基礎(chǔ)和主要手段，至今共構(gòu)建了24張蠶豆遺傳圖譜，最為飽和的圖譜覆蓋長度為1 439 cM，平均遺傳距離為0.8 cM。同時基于圖譜開展了產(chǎn)量相關(guān)性狀、品質(zhì)相關(guān)性狀以及抗性相關(guān)性狀的QTL定位，有效地促進了蠶豆的遺傳與基因研究。

關(guān)鍵詞： 蠶豆;遺傳連鎖圖譜;QTL定位

中圖分類號： S529 文獻標識碼： A 文章編號： 1000-4440（2021）01-0237-09

Research advance of genetic linkage map and QTL location in Vicia faba L.

ZHOU Xian-li1，2， TENG Chang-cai2， ZHANG Hong-yan1，2， LIN Xi1，2， LIU Yu-ling1，2， WU Xiao-yan1，2，HOU Wan-wei2， LIU Yu-jiao1，2，3

（1.Qinghai University， Xining 810016， China;2.Qinghai Academy of Agricultural and Forestry Sciences， Xining 810016， China;3.The Co-constructing State Key Laboratory of Three Rivers Sources Ecology and Plateau Agriculture and Animal Husbandry， Qinghai University， Xining 810016， China）

Abstract： As a kind of edible legumes which can be used as food， vegetables feed and green manure， Vicia faba L. plays a crucial role in human health and soil improvement. The construction of high-density genetic linkage map is the basis and main method for genome research in V. faba L.. 24 genetic maps have been constructed so far， in which the most saturated map spanned 1 439 cM with an average genetic distance of 0.8 cM. Quantitative trait locus （QTL） mapping on yield-related traits， quality-related traits， and resistance-related traits had been carried out， which effectively promoted inheritance and gene research in V. faba L..

Key words： Vicia faba L.;genetic linkage maps;quantitative trait locus （QTL） location

蠶豆（Vicia faba L.）具有“養(yǎng)人、養(yǎng)畜、養(yǎng)地”的作用，是中國傳統(tǒng)的出口農(nóng)產(chǎn)品，也是青海省種植業(yè)結(jié)構(gòu)調(diào)整和農(nóng)民增收的優(yōu)勢作物之一[1-3]。隨著人們生活水平的提高和健康飲食需求的增加，菜用蠶豆的市場需求量越來越大。蠶豆作為一種輪作豆類，在土壤改良方面也發(fā)揮著重要作用。

遺傳連鎖圖譜是指已知基因或分子標記在連鎖群上的線性排列，標記間的相對遺傳距離通過重組率計算[4-5]。蠶豆基因組龐大（約13 Gb），缺乏可用的基因組序列數(shù)據(jù)，導致其基因組研究進展緩慢[6]。構(gòu)建高密度遺傳圖譜，進行QTL定位以及開發(fā)連鎖標記是蠶豆基因組研究的主要手段，對控制重要性狀的關(guān)鍵基因定位與克隆、關(guān)聯(lián)與連鎖分析以及分子標記輔助選擇優(yōu)良品種具有重要作用。

基于蠶豆遺傳圖譜，對產(chǎn)量、品質(zhì)、抗性等重要性狀進行基因與QTL定位，并針對部分性狀開發(fā)連鎖標記，有利于替代繁瑣且耗時的表型評估，促進分子標記輔助選擇（MAS）育種的發(fā)展，進而縮短育種進程[7]。本文對蠶豆遺傳圖譜的構(gòu)建與重要農(nóng)藝性狀、品質(zhì)性狀以及抗性的QTL定位的研究進展進行綜述，旨在為蠶豆遺傳連鎖圖譜的加密、重要性狀的基因精細定位和MAS育種等提供參考。

1 蠶豆遺傳連鎖圖譜的構(gòu)建

蠶豆遺傳連鎖圖譜的構(gòu)建起源于1991年，Van de Ven等[8]報道了蠶豆第1張遺傳圖譜。在隨后的近30年里，基于12種遺傳標記共構(gòu)建了24張遺傳圖譜（表1）。蠶豆遺傳圖譜的研究進程可大致分為以下2個階段：第1階段是利用形態(tài)學（Morphological markers）、同工酶（Isozymes）等遺傳標記與隨機擴增多態(tài)性（Random amplified polymorphic DNA，RAPD）標記、微衛(wèi)星（Microsatellite）等傳統(tǒng)標記構(gòu)建蠶豆遺傳圖譜;第2階段是利用高通量測序技術(shù)開發(fā)的表達序列標簽（Expressed sequence tag，EST）、單核苷酸多態(tài)性（Single nucleotide polymorphisms，SNP）等標記構(gòu)建蠶豆遺傳圖譜。

1.1 基于形態(tài)學、RAPD等標記的蠶豆遺傳圖譜

早期，國外學者對種皮顏色等形態(tài)學性狀與一些同工酶系統(tǒng)進行了報道，并研究了其基因在染色體上的位置[9-11]。這些形態(tài)學性狀和同工酶系統(tǒng)被作為遺傳標記廣泛應(yīng)用于蠶豆遺傳圖譜的構(gòu)建[8， 12-14]?；谕っ?、限制性片段長度多態(tài)性（Restriction fragment length polymorphism，RFLP）等傳統(tǒng)標記構(gòu)建的遺傳圖譜的標記數(shù)量較少，密度較低。Van de Ven等[8]構(gòu)建的第1張圖譜僅包含7個標記;Romn等[15]構(gòu)建的遺傳圖譜包含121個標記，平均遺傳距離為13.77 cM。構(gòu)建的蠶豆遺傳圖譜包含的連鎖群數(shù)量多于蠶豆染色體數(shù)（2n=12）。Satovic等[13]構(gòu)建的圖譜包含了48個連鎖群，而僅有6個連鎖群被錨定到了染色體上。

在作圖群體構(gòu)建上，常以Vf6[12-15]、Vf136[13，15-16]、Vf27[17-18]等品系作為親本，構(gòu)建回交群體（Backcross population，BC population）[8]、F2群體以及重組自交系（Recombinant inbred lines，RIL）等群體作為作圖群體以構(gòu)建圖譜。以蠶豆三體材料為親本構(gòu)建作圖群體，可以將連鎖群與染色體相關(guān)聯(lián)。蠶豆初級三體是以Vf6為母本與其他材料雜交所形成的，用于檢測建立遺傳圖譜中連鎖群與染色體之間的對應(yīng)關(guān)系[13-14]。Satovic等[13]將Vf6分別與不同蠶豆染色體三體材料雜交，形成了具有3、4、5、6號染色體三體的7個F2作圖群體，所構(gòu)建的圖譜包含1個形態(tài)學標記、9個同工酶標記與147個RAPD標記。由于單一群體間多態(tài)性有限，所以基于多個作圖群體構(gòu)建了復合圖譜，錨定了更多的遺傳標記，增加了圖譜密度[13-14，19]。

2002年，Poarkova等[20]從蠶豆特異性DNA文庫中開發(fā)的位于蠶豆1號染色體上的微衛(wèi)星標記也被廣泛應(yīng)用于蠶豆遺傳圖譜的構(gòu)建[16，19，21]。Ellwood等[17]則利用RIL群體構(gòu)建了包含151個內(nèi)含子靶向擴增多態(tài)性（Intron-targeted amplified polymorphic，ITAP）標記、12個連鎖群、覆蓋1 685.8 cM的遺傳圖譜，每個連鎖群上的標記數(shù)量由3個到30個不等，長度從23.6 cM到324.8 cM不等，與蒺藜苜蓿染色體組之間的同源分析結(jié)果表明兩種作物之間具有高度的共線性。

1.2 基于高通量測序技術(shù)的蠶豆遺傳圖譜

隨著分子標記和高通量測序技術(shù)的不斷發(fā)展，EST衍生標記[22-23]、SNP標記開始被應(yīng)用于蠶豆遺傳圖譜的構(gòu)建。為了對蠶豆開花與產(chǎn)量性狀進行QTL定位，Cruz-Izquierdo等[22]基于Vf6×Vf27衍生的RIL群體構(gòu)建了具有16個連鎖群、258個標記的遺傳圖譜，其中EST標記167個，其他標記則為同工酶、RAPD等標記，并將構(gòu)建的遺傳圖譜與蒺藜苜蓿和兵豆的染色體進行了共線性分析，為建立蠶豆基因組與其他豆科作物基因組之間的同源性提供了錨定位點。自2013年起，國內(nèi)開始采用EST（表達序列標簽）序列設(shè)計和已發(fā)表的簡單重復序列（Simple sequence repeats，SSR）標記構(gòu)建蠶豆遺傳連鎖圖譜，并對圖譜進行加密[5，24-25]。在Ma等[24]構(gòu)建的SSR遺傳圖譜基礎(chǔ)上，Yang等[25]繼續(xù)以91825×K1563構(gòu)建的F2群體為作圖群體進行加密，最終構(gòu)建了包含465個SSR標記、7個連鎖群、覆蓋長度為4 516.75 cM、平均遺傳距離為9.71 cM的遺傳圖譜，連鎖群上的SSR標記數(shù)量從12到136不等，連鎖群的長度范圍為129.35 cM至1 180.21 cM，圖譜覆蓋度極大提升，但平均遺傳距離仍較大。

2014年起，SNP標記開始被應(yīng)用于蠶豆遺傳圖譜構(gòu)建[26-28]。至今最飽和的蠶豆遺傳圖譜是由Sudheesh等[29]構(gòu)建的整合圖譜，Sudheesh等基于SNP標記將已發(fā)表的遺傳圖譜與其研究中新構(gòu)建的遺傳圖譜相整合，構(gòu)建了一張包含6個連鎖群、1 850個標記的整合圖譜，其覆蓋總長度為1 439 cM，平均遺傳距離為0.8 cM。最近，Carrillo-Perdomo等[30]利用在3個重組品系中表現(xiàn)多態(tài)性的1 819個SNP標記構(gòu)建了一張包含1 728個標記、圖譜長度為1 547.71 cM、平均遺傳距離為0.89 cM的高密度蠶豆遺傳圖譜。

分子標記與遺傳圖譜的構(gòu)建是蠶豆基因組研究中的有效手段。目前，RAPD、SSR、SNP等多種不同類型的遺傳標記被應(yīng)用于圖譜構(gòu)建，所構(gòu)建的蠶豆遺傳圖譜覆蓋長度增加，標記數(shù)量增多，并可對多個圖譜進行有效整合。飽和遺傳連鎖圖譜的構(gòu)建和與目標性狀緊密連鎖的分子標記的開發(fā)，對MAS育種中育種進程的縮短具有重要意義。在構(gòu)建的遺傳圖譜中，13張遺傳圖譜被應(yīng)用于蠶豆產(chǎn)量與籽粒相關(guān)性狀、開花特性、抗性等重要性狀的QTL定位。

1.3 與其他豆科作物的同源比對

將蠶豆遺傳圖譜與其他豆科作物進行同源比對，以確定蠶豆與其他作物在系統(tǒng)發(fā)育中的親緣關(guān)系，是對控制重要性狀的關(guān)鍵基因進行同源克隆與候選基因篩選的主要手段之一。將應(yīng)用來源于其他豆科作物的通用性引物構(gòu)建的蠶豆遺傳圖譜與小扁豆、蒺藜苜蓿的基因組進行同源比對，結(jié)果表明3種作物間均存在高度的共線性，蠶豆與小扁豆之間的同源水平更高;與蒺藜苜蓿的染色體比較，發(fā)現(xiàn)蠶豆和小扁豆的染色體均存在中等程度的重排，這種重排可以解釋不同作物間染色體數(shù)目的差異[17， 22]。將基于測序技術(shù)構(gòu)建的蠶豆遺傳圖譜與以上2種作物進行同源比對，發(fā)現(xiàn)蠶豆和蒺藜苜?；蚪M間存在大量的高度共線性區(qū)域，并在映射區(qū)間內(nèi)通過序列比對驗證了1個控制單寧的候選基因[28]。因此，同源比對有助于在具有高度共線性的映射區(qū)間內(nèi)預測基因含量與距離，進行基因注釋。

不同作物間的共線性程度在一定程度上反映了作物在系統(tǒng)發(fā)育中親緣關(guān)系的遠近。蠶豆與蒺藜苜蓿、鷹嘴豆、豌豆等親緣關(guān)系較近，因此在同源比對時常將蠶豆與這些豆科作物進行比較。Kaur 等[26]與Carrillo-Perdomo等[30]均將其所構(gòu)建的圖譜與多種豆科作物的基因組進行了比較。Carrillo-Perdomo等[30]發(fā)現(xiàn)蠶豆遺傳圖譜中的連鎖群與豌豆、鷹嘴豆、蒺藜苜蓿等近緣作物染色體之間的共線性區(qū)域幾乎覆蓋了整個染色體組，而蠶豆與菜豆、短豇豆、百脈根之間的親緣關(guān)系較遠，故其染色體間存在的共線性區(qū)域稍少。

2 蠶豆重要性狀相關(guān)QTL研究進展

2.1 重要農(nóng)藝性狀與產(chǎn)量相關(guān)QTL

產(chǎn)量是蠶豆育種中關(guān)注的主要問題，單株莢數(shù)、單莢粒數(shù)以及百粒質(zhì)量共同構(gòu)成了蠶豆單株產(chǎn)量[31-38]。產(chǎn)量形成過程受到多個性狀的影響，如開花時間、生長習性等。開花時間主要影響蠶豆品種的成熟特性和環(huán)境適應(yīng)性。Cruz-Izquierdo等[22]鑒定并驗證了控制開花時間、花期長度、莢長、單莢粒數(shù)、單莢胚珠數(shù)等5個性狀的QTL，在2、5、6號染色體上連續(xù)2年鑒定出12個穩(wěn)定的QTL，比較分析結(jié)果表明控制開花時間的基因組區(qū)域在其他豆科作物中具有保守性。Avila等[33]對花序特征、產(chǎn)量相關(guān)的8個性狀進行了QTL定位，檢測到大量的QTL位點，在開花節(jié)位、結(jié)莢高度、有效分枝數(shù)方面各檢測到4個QTL，在單花序花蕾數(shù)、百粒質(zhì)量方面各檢測到5個QTL，在主莖開花層數(shù)方面檢測到6個具有加性效應(yīng)的QTL，但QTL的準確性和穩(wěn)定性需進一步驗證。

百粒質(zhì)量和粒型性狀是重要的商品經(jīng)濟特征。1999年，與百粒質(zhì)量相關(guān)的QTL被首次報道[14]，QTL主要聚集在6號染色體上，與20個標記位點顯著相關(guān)，主效QTL位點（與OPM1811725標記連鎖）可解釋將近30%的表型變異率。田瑩瑩等[37]以粒型差異較大的云122與TF42為親本配置雜交組合，構(gòu)建F2群體對粒型性狀進行了QTL定位，最終檢測到4個控制籽粒長、寬、質(zhì)量的QTL（qSL-1-1、 qSW-1-2 、qSH-1-3、qSH-1-4）。

目前，產(chǎn)量方面的QTL研究較少，定位到的QTL區(qū)間較大，其穩(wěn)定性與準確性需要進一步驗證，難以開發(fā)可直接應(yīng)用于MAS育種的分子標記。對于質(zhì)量性狀，如子葉顏色和生長習性，連鎖標記的開發(fā)難度較低。沙偉超[39]采用分組混合分析（BSA）法篩選到9個與子葉顏色相連鎖的SSR標記，并初步將控制子葉顏色的基因定位到LG05，但未對連鎖標記進行驗證。Avila等[40-41]開發(fā)了用于有限生長習性選擇的酶切擴增多態(tài)性序列（Cleaved amplified polymorphic sequences，CAPS）標記Vf_TFL1，在MAS育種中具有實用性。

2.2 品質(zhì)性狀相關(guān)標記

蠶豆營養(yǎng)價值高， 富含蛋白質(zhì)、淀粉和微量元素[42]。蠶豆蛋白質(zhì)是一種優(yōu)質(zhì)的谷物蛋白質(zhì)，蛋白質(zhì)組分含量不同階段處于不斷變化中[43]。Macas等[44-45]報道了5個控制蠶豆籽粒蛋白質(zhì)的基因，并公布了相關(guān)引物序列。蠶豆籽粒中抗營養(yǎng)因子的存在降低了蛋白質(zhì)的生物學價值，在育種進程中選育抗營養(yǎng)因子含量較低的品種為品質(zhì)育種提供了新的研究方向[46]。Gutierrez等[47-49]鑒定了與單寧缺失基因（zt-1、zt-2）緊密連鎖的2個SCAR標記，報道了與蠶豆嘧啶葡糖苷和伴蠶豆嘧啶核苷含量相關(guān)基因緊密連鎖的1個CAPs標記。Hou等[50]利用596個SSR標記和100個ISSR標記篩選到1個與zt-1緊密連鎖的SSR標記（SSR84），可以準確地預測zt-1基因型。

2.3 抗性相關(guān)QTL和基因定位

病害是蠶豆產(chǎn)量的主要限制因素之一，全世界至2018年為止報道了銹病、褐斑病、霜霉病等10種蠶豆常見真菌病害[51]?？剐缘膹碗s性導致難以僅僅依據(jù)表型觀測進行抗性育種，因此有必要對抗性的遺傳機制進行研究。目前，針對銹病、褐斑病、耐寒性、抗旱性的相關(guān)QTL進行了研究，赤斑病、霜霉病等病害的遺傳研究基礎(chǔ)較薄弱，僅見病原鑒定、防治措施及遺傳機制等方面的研究報道[52-54]。

2.3.1 銹病與褐斑病抗性 蠶豆銹病和褐斑病是真菌病害，對中國秋播型蠶豆危害嚴重[55]。典型褐斑病發(fā)生時，會造成35%～40%減產(chǎn)[56]。抗病育種過程較為復雜，已育成的抗病品種由于病菌生理小種的變化也易喪失抗性，目前難以獲得完全抗性。Sillero等[57]在648份蠶豆種質(zhì)中篩選出6份（V-300、V-1271、V-1273、V-313、V-1272、V-1335）具有銹病抗性的資源?；诳逛P病品系 2N52和敏感品系VF-176雜交產(chǎn)生的F2群體， Avila等[58]采用分組混合分析法檢測到了3個與抗銹病基因（Uvf-1）緊密連鎖的RAPD標記，在標記OPI20900和Uvf-1間未檢測到重組，采用該標記對抗銹病基因型的選擇效率較高。

蠶豆褐斑病抗性的等位基因多來源于Vf6、29H、Ascot等品系，在蠶豆抗褐斑病相關(guān)QTL定位中較多地使用了這些品系。初期，基于Vf6×Vf136和29H×Vf136衍生的群體，Román等[56]檢測到的控制蠶豆褐斑病抗性的QTL （Af1）與Avila等[16]檢測到的Af3都位于3號染色體，表明二者可能位于同一基因組區(qū)域。基于Icarus× Ascot衍生的RIL群體，Kaur等[26]在2年的評估中檢測到4個相關(guān)QTL，分別位于Chr-I.A、Chr-II與Chr-VI，其中，3號位點可能在先前的研究中被報道過，而其他3個則是新檢測到的位點。最近，Sudheesh等[29]構(gòu)建了一張高密度的整合圖譜用于蠶豆褐斑病抗性相關(guān)QTL定位，作圖群體為Nura×Farah衍生的RIL，Nura是從 Icarus和Ascot雜交產(chǎn)生的系譜中選擇出來的，其抗性主要來源于Ascot，最終確定了2個基因組區(qū)域（AB_N1、AB_N2），可解釋多達49%的表型變異，2個QTL與Kaur等[26]檢測到的1號位點和4號位點為相同的位點。

2.3.2 列當抗性 在地中海沿岸，列當對于蠶豆的寄生極其嚴重，難以有效控制[46]。Díaz-Ruiz等[59-60]基于Vf6×Vf136的RIL群體檢測到了4個相關(guān)QTL區(qū)域（Oc2、Oc3、Oc4、Oc5），其中Oc4、Oc5被定位在1號染色體。Díaz-Ruiz等[32]在不同的世代和環(huán)境中對列當抗性相關(guān)QTL區(qū)域進行了驗證，最終在不同環(huán)境中的3個RIL群體中鑒定出4個控制列當抗性的QTL，而先前報道的在F2群體中檢測到的Oc1在高代群體中則不顯著。Ocaa-Moral等[36]基于SNP標記對Gutiérrez等[23]檢測到的QTL進行了驗證，發(fā)現(xiàn)在6號染色體上檢測到的Oc7在不同年份中均解釋了較大的表型變異，并在QTL區(qū)間內(nèi)增添了2個SNP標記，減小了標記間隙，而Oc8位點僅在1年中被檢測到。

2.3.3 非生物抗性 MAS育種比常規(guī)育種更快速、更有效，尤其是對于性狀表達受環(huán)境影響較大的非生物脅迫而言[61]。蠶豆的耐寒性是保證秋播型蠶豆安全越冬的前提?；蛐秃铜h(huán)境的相互作用降低了耐寒性選擇的有效性，傳統(tǒng)育種程序?qū)τ谇锊バ头N質(zhì)的遺傳改良進展緩慢[62]。Arbaoui等和Sallam等基于SNP構(gòu)建的遺傳圖譜檢測到與耐寒性和耐寒性生理相關(guān)的多個重要的推定QTL[31，63]。然而，推定的QTL無法直接應(yīng)用于分子標記輔助選擇育種。隨后，Sallam等[64]使用101個RIL和189個不同基因型的單粒傳后代作為遺傳背景，采用QTL定位和全基因組關(guān)聯(lián)分析，檢測到5個與耐寒性和脂肪酸組分相關(guān)的SNP位點，經(jīng)驗證在2個群體中有1個SNP標記（VF_Mt3g086600）與抗凍性相關(guān)基因緊密連鎖，這個標記在先前研究中被證明與產(chǎn)量性狀相關(guān)。

蠶豆在其生長過程對于干旱反應(yīng)相對敏感，開花期和結(jié)莢期的干旱脅迫會導致蠶豆的嚴重減產(chǎn)[65]。Khazaei等[66]以211個RIL為作圖群體，應(yīng)用來自蒺藜苜蓿的SNP檢測到15個與氣孔特性相關(guān)的推定QTL，并利用蠶豆與蒺藜苜蓿之間的共線性關(guān)系鑒定了位于蠶豆2號染色體上的8個QTL（qSD-2013-1、qSD-2014、qSL-2013、qSL-2014、qGS-2013、qCT-2012、qCT-2013、qCT-2014），這些QTL可能是與抗旱性相關(guān)的候選基因。

3 問題與展望

3.1 蠶豆遺傳連鎖圖譜的構(gòu)建

高質(zhì)量遺傳連鎖圖譜的構(gòu)建和與目標性狀緊密連鎖的分子標記的開發(fā)，對MAS育種中育種效率的提高具有重要意義。隨著分子標記技術(shù)的發(fā)展，大量可轉(zhuǎn)移的EST、SNP等標記被用于構(gòu)建整合圖譜，在原有圖譜的基礎(chǔ)上增添新的標記，極大地增加了圖譜的分辨率。高質(zhì)量的蠶豆遺傳圖譜也被用于與其他近緣豆科植物模型進行同源比對[26，30]。然而。與小麥、大豆等主要農(nóng)作物相比[67-68]，蠶豆遺傳連鎖圖譜仍未達到飽和。

在蠶豆遺傳圖譜構(gòu)建方面主要存在以下問題：①應(yīng)用于圖譜構(gòu)建的作圖群體較小，單一群體內(nèi)單株數(shù)量均在200株以下，影響了遺傳圖譜的精度。②已報道的蠶豆遺傳圖譜大多基于傳統(tǒng)標記，基于SNP標記構(gòu)建的遺傳圖譜僅有6張，開發(fā)的易于轉(zhuǎn)移的標記數(shù)量較少;③已構(gòu)建的高密度遺傳連鎖圖譜的總長度多在1 500 cM左右，難以覆蓋蠶豆龐大的基因組。部分基于SSR及其他標記的遺傳圖譜雖然總長度較大，但分辨率不足。

因此，應(yīng)通過擴大作圖群體，開發(fā)易轉(zhuǎn)移的EST、SNP等標記以及對已發(fā)表的圖譜進行整合等方式，進一步促進蠶豆遺傳圖譜的飽和。將所構(gòu)建的圖譜與已發(fā)表的圖譜進行整合，是增加圖譜飽和度的有效手段之一。同時由于整合圖譜包含了來自其他圖譜的標記，也有利于對不同研究中檢測到的QTL進行比較鑒定，以篩選在多個環(huán)境中穩(wěn)定表達的QTL;其次，蠶豆飽和遺傳圖譜的構(gòu)建需要更多的基因組信息，如EST測序信息、細菌人工染色體文庫等，以進一步開發(fā)更多易于轉(zhuǎn)移的標記，增加圖譜飽和程度，提高與重要性狀相關(guān)的QTL區(qū)域的分辨率。

3.2 蠶豆重要性狀相關(guān)QTL檢測

基于遺傳圖譜，大量的QTL被鑒定，涉及開花特性、產(chǎn)量、粒型、子葉顏色、籽粒蛋白質(zhì)含量、銹病抗性、褐斑病抗性、列當抗性、耐寒性、抗旱性等多個重要性狀，開發(fā)了與控制有限生長習性、單寧缺失、銹病抗性、耐寒性等基因緊密連鎖的分子標記。

在蠶豆重要性狀相關(guān)QTL定位方面，存在以下問題：①Q(mào)TL定位研究多集中于抗性領(lǐng)域，產(chǎn)量方面的QTL報道較少，其穩(wěn)定性和準確性需進一步驗證;②蠶豆遺傳圖譜的飽和度極大地限制了重要性狀的基因或QTL的精細定位和克隆[66]。對于多基因控制的數(shù)量性狀，由于無法確定其QTL的準確區(qū)域，難以找到與其緊密連鎖的分子標記，導致其實際應(yīng)用價值較低。

QTL所在基因組區(qū)域的飽和并在多個環(huán)境和遺傳背景中進行驗證是獲得可靠標記的必要前提[69]。因此，應(yīng)在QTL所在區(qū)域進一步設(shè)計引物，縮小定位區(qū)間。另外，通過蠶豆與豆科模式作物的比較基因組學分析可以有效地驗證重要性狀相關(guān)QTL位點，確定可能的候選基因[70-72]，篩選用于后代單株選擇的分子標記。同時，開發(fā)用于目標性狀選擇的分子標記，以在育種早期對含有目標性狀的個體進行快速、準確的選擇，加快育種進程。

參考文獻：

[1] 葉 茵. 中國蠶豆學 [M]. 北京： 中國農(nóng)業(yè)出版社， 2003.

[2] 劉玉皎. 調(diào)結(jié)構(gòu)，轉(zhuǎn)方式，促進青海蠶豆產(chǎn)業(yè)轉(zhuǎn)型升級 [J]. 青?？萍?， 2018， 25（1）： 35-37.

[3] 周俊玲，張蕙杰. 世界蠶豆生產(chǎn)及貿(mào)易形勢分析 [J]. 世界農(nóng)業(yè)， 2016（11）： 107-111.

[4] 馮貝貝，靳 娟，楊 磊，等. 冬棗 （Ziziphus jujube） 果實重要性狀和遺傳圖譜構(gòu)建及QTL定位研究進展 [J]. 分子植物育種， 2019， 17（24）： 8184-8190.

[5] 姜俊燁. 蠶豆微核心種質(zhì)構(gòu)建及SSR遺傳連鎖圖譜加密 [D]. 北京：中國農(nóng)業(yè)科學院， 2014.

[6] COOPER J W， WILSON M H， DERKS M F L， et al. Enhancing faba bean （Vicia faba L.） genome resources [J]. Journal of Experimental Botany， 2017， 68（8）： 1941-1953.

[7] 瞿華香，張玉燭，張 岳，等. 分子標記輔助選擇育種研究進展 [J]. 作物研究， 2008， 22 （S1）： 355-358.

[8] VAN DE VEN W T G， WAUGH R， DUNCAN N， et al. Development of a genetic linkage map in Vicia faba using molecular and biochemical techniques [J]. Aspects of Applied Biology， 1991， 27： 49-54.

[9] SUSO M J， MORENO M T. Genetic control of electrophoretic variation for glutamate oxalacetate transaminase （GOT） in Vicia faba L.[J]. Fabis Newsletter Faba Bean Information Service， 1982， 5： 14.

[10]CABRERA A， CUBERO J I， MARTIN A. Genetic mapping using trisomics in Vicia faba L. [J]. Fabis Newsletter Faba Bean Information Service， 1989， 23： 5-7.

[11]SJKDIN J. Induced translocation in Vicia faba L. [J]. Hereditas， 1971， 68： 1-34.

[12]TORRES A M， WEEDEN N F， MARTIN A. Linkage among isozyme， RFLP and RAPD markers in Vicia faba [J]. Theoretical and Applied Genetics， 1993， 85（8）： 937-945.

[13]SATOVIC Z， TORRES A M， CUBERO J I. Genetic mapping of new morphological， isozyme and RAPD markers in Vicia faba L. using trisomics [J]. Theoretical and Applied Genetics， 1996， 93： 1130-1138.

[14]PATTO M C V， TORRES A M， KOBLIZKOVA A， et al. Development of a genetic composite map of Vicia faba using F2 populations derived from trisomics plants [J]. Theoretical and Applied Genetics， 1999， 98： 736-743.

[15]ROMN B， TORRES A M， RUBIALES D， et al. Mapping of quantitative trait loci controlling broomrape （Orobanche crenata Forsk.） resistance in faba bean （Vicia faba L.） [J]. Genome， 2002， 45（6）： 1057-1063.

[16]AVILA C M， ATOVIC Z， SILLERO J C， et al. Isolate and organ-specific QTLs for ascochyta blight resistance in faba bean [J]. Theoretical and Applied Genetics， 2004， 108： 1071-1078.

[17]ELLWOOD S R， PHAN H T， JORDAN M， et al. Construction of a comparative genetic map in faba bean （Vicia faba L.）; conservation of genome structure with Lens culinaris [J]. BMC Genomics， 2008， 9： 380.

[18]DAZ-RUIZ R， SATOVIC Z， AVILA C M， et al. Conrmation of QTLs controlling Ascochyta fabae resistance in different generations of faba bean （Vicia faba L.） [J]. Crop and Pasture Science， 2009， 60： 353-361.

[19]ROMN B， SATOVIC Z， POZARKOVA D， et al. Development of a composite map in Vicia faba， breeding applications and future prospects [J]. Theoretical and Applied Genetics， 2004， 108： 1079-1088.

[20]POARKOVA D， KOBLIKOVA A， ROMN B， et al. Development and characterization of microsatellite markers from chromosome 1-specific DNA libraries of Vicia faba [J]. Biologia Plantarum， 2002， 45： 337-345.

[21]AVILA C M， ATOVIC Z， SILLERO J C， et al. QTL detection for agronomic traits in faba bean （Vicia faba L.） [J]. Agriculturae Conspectus Scientificus， 2005， 70 （3）： 65-73.

[22]CRUZ-IZQUIERDO S， AVILA C M， SATOVIC Z， et al. Comparative genomics to bridge Vicia faba with model and closely-related legume species： stability of QTLs for flowering and yield-related traits [J]. Theoretical and Applied Genetics， 2012， 125： 1767-1782.

[23]GUTIRREZ N， PALOMINO C， SATOVIC Z， et al. QTLs for Orobanche spp. resistance in faba bean： identification and validation across different environments [J]. Molecular Breeding， 2013， 32（4）： 909-922.

[24]MA Y， BAO S Y， YANG T， et al. Genetic linkage map of Chinese native variety faba bean （Vicia faba L.） based on simple sequence repeat markers [J]. Plant Breeding， 2013， 132 （4）： 397- 400.

[25]YANG T， JIANG J Y， ZHANG H Y， et al. Density enhancement of a faba bean genetic linkage map （Vicia faba） based on simple sequence repeats markers [J]. Plant Breeding， 2019， 138 （2）： 207-215.

[26]KAUR S， KIMBER R B E， COGAN N O I， et al. SNP discovery and high-density genetic mapping in faba bean （Vicia faba L.） permits identication of QTLs for ascochyta blight resistance [J]. Plant Science， 2014， 217： 47-55.

[27]SALLAM A， ARBAOUI M， ELESAWI M A， et al. Identification and verification of QTL associated with frost tolerance using linkage mapping and GWAS in winter faba bean [J]. Frontiers in Plant Science， 2016， 7： 1098.

[28]WEBB A， COTTAGE A， WOOD T， et al. A SNP based consensus genetic map for synteny-based trait targeting in faba bean （Vicia faba L.） [J]. Plant Biotechnology Journal， 2016， 14： 177-185.

[29]SUDHEESH S， KIMBER R B E， BRAICH S， et al. Construction of an integrated genetic linkage map and detection of quantitative trait loci for ascochyta blight resistance in faba bean （Vicia faba L.） [J]. Euphytica， 2019， 215 （3）： 42.

[30]CARRILLO-PERDOMO E， VIDAL A， KREPLAK J， et al. Development of new genetic resources for faba bean （Vicia faba L.） breeding through the discovery of gene-based SNP markers and the construction of a high-density consesus map [J]. Scientific Reports， 2020， 10： 6790.

[31]ARBAOUI M， LINK W， SATOVIC Z， et al. Quantitative trait loci of frost tolerance and physiologically related trait in faba bean （Vicia faba L.） [J]. Euphytica， 2008， 164： 93-104.

[32]DAZ-RUIZ R， TORRES A M， SATOVIC Z， et al. Validation of QTLs for Orobanche crenata resistance in faba bean （Vicia faba L.） across environments and generations [J]. Theoretical and Applied Genetics， 2010， 120： 909-919.

[33]AVILA C M， RUIZ-RODRGUEZ M D， CRUZ-IZQUIERDO S， et al. Identification of plant architecture and yield-related QTL in Vicia faba L. [J]. Molecular Breeding， 2017， 37：88.

[34]SATOVIC Z， AVILA C M， CRUZ-IZQUIERDO S， et al. A reference consensus genetic map for molecular markers and economically important traits in faba bean （Vicia faba L.） [J]. BMC Genomics， 2013， 14（1）： 932.

[35]CATT S C， BRAICH S， KAUR S， et al. QTL detection for flowering time in faba bean and the responses to ambient temperature and photoperiod [J]. Euphytica， 2017， 213： 125.

[36]OCAA-MORAL S， GUTIRREZ N， TORRES A M， et al. Saturation mapping of regions determining resistance to ascochyta blight and broomrape in faba bean using transcriptome-based SNP genotyping [J]. Theoretical and Applied Genetics， 2017， 130： 2271-2282.

[37]田瑩瑩，侯萬偉，劉玉皎. 蠶豆粒型性狀的遺傳分析及QTL檢測 [J]. 分子植物育種， 2018， 16 （4）： 1174-1183.

[38]楊生華，劉 榮，楊 濤，等. 蠶豆種質(zhì)資源種子表型性狀精準評價[J]. 中國蔬菜， 2016 （10）： 32-40.

[39]沙偉超. 蠶豆子葉顏色性狀的SSR標記研究 [D]. 西寧：青海大學， 2017.

[40]AVILA C M， NADAL S， MORENO M T， et al. Development of a simple PCR-based marker for the determination of growth habit in Vicia faba L. using a candidate gene approach [J]. Molecular Breeding， 2006， 17： 185-190.

[41]AVILA C M， ATIENZA S G， MORENO M T， et al. Development of a new diagnostic marker for growth habit selection in faba bean （Vicia faba L.） breeding [J]. Theoretical and Applied Genetics， 2007， 115： 1075-1082.

[42]彭 葵，李錦鴻，李育軍，等. 蠶豆的營養(yǎng)與加工研究 [J]. 長江蔬菜， 2019（12）： 42-45.

[43]劉珍珍，陳友霞，楊文藝，等. 未成熟蠶豆蛋白組分的分析 [J]. 食品研究與開發(fā)， 2019， 40（23）： 82-86.

[44]MACAS J， DOLEZEL J， LUCRETTI S， et al. Localization of seed genes on flow-sorted field bean chromosomes [J]. Chromosome Research， 1993， 1： 107-115.

[45]MACAS J， WESCHKE W， BUMLEIN H， et al. Localization of vicilin genes via polymerase chain reaction on microisolated field bean chromosomes [J]. Plant Journal， 1993， 3： 883-886.

[46]TORRES A M， AVILA C M， GUTIERREZ N， et al. Marker-assisted selection in faba bean （Vicia faba L.） [J]. Field Crops Research， 2010， 115： 243-252.

[47]GUTIERREZ N， AVILA C M， RODRIGUEZ-SUAREZ C， et al. Development of SCAR markers linked to a gene controlling absence of tannins in faba bean [J]. Molecular Breed， 2007， 19： 305-314.

[48]GUTIERREZ N， AVILA C M， MORENO M T， et al. Development of SCAR markers linked to zt-2， one of the genes controlling absence of tannins in faba bean [J]. Australian Journal of Agricultural Research， 2008， 59： 62-68.

[49]GUTIERREZ N， AVILA C M， DUC G， et al. CAPs markers to assist selection for low vicine and convicine content in faba bean （Vicia faba L.） [J]. Theoretical and Applied Genetics， 2006， 114： 59-66.

[50]HOU W W， ZHANG X J， YAN Q B， et al. Linkage map of a gene controlling zero tannins （zt-1） in faba bean （Vicia faba L.） with SSR and ISSR markers [J]. Agronomy， 2018， 8： 80.

[51]李仁慧，閆智臣，段廷玉. 蠶豆真菌病害及其研究進展 [J]. 草業(yè)科學， 2019， 36（8）： 1976-1987.

[52]孫雪梅，易紅娟，楊凌峰，等. 多種殺菌劑對蠶豆霜霉病田間防治效果比較 [J]. 農(nóng)藥科學與管理， 2019， 40（7）： 60-63.

[53]李 龍，張 蕓，郭延平，等. 8種殺菌劑對春蠶豆赤斑病的防治效果 [J]. 植物保護， 2019， 45（3）： 245-248.

[54]杜成章，龍玨臣，龔萬灼，等. 蠶豆赤斑病抗性的主基因+多基因遺傳分析 [J]. 植物保護， 2019， 45（6）： 131-137.

[55]王海飛，宗緒曉. 蠶豆種質(zhì)資源、抗病育種和 QTL 定位及抗逆性研究進展 [J]. 植物遺傳資源學報， 2011， 12 （2）： 259-270.

[56]ROMN B， SATOVIC Z， AVILA C M， et al. Locating genes associated with Ascochyta fabae resistance in Vicia faba[J]. Australian Journal of Agricultural Research， 2003， 54： 85-90.

[57]SILLERO J C， MORENO M T， RUBIALES D. Characterization of new sources of resistance to Uromyces viciae-fabae in a germplasm collection of Vicia faba [J]. Plant Pathology， 2000， 49： 389-395.

[58]AVILA C M， SILLERO J C， RUBIALES D， et al. Identification of RAPD markers linked to the Uvf-1 gene conferring hypersensitive resistance against rust （Uromyces viciae-fabae） in Vicia faba L. [J]. Theoretical and Applied Genetics， 2003， 107： 353-358.

[59]DAZ-RUIZ R， ROMN B， SATOVIC Z， et al. Validation of QTLs for broomrape resistance in an F6 population of faba bean [C]//AEP. Conference handbook-european conference on grain legumes with international conference on legume genomics & genetics. Dijon， France： European Association for Grain Legume Research， 2004： 122.

[60]DAZ-RUIZ R， SATOVIC Z， ROMN B， et al. QTL analysis of broomrape resistance in faba bean （Vicia faba L.） [C]//KOVACEVIC V， JOVANOVAC S. Proceedings of XL croatian symposium on agriculture. Opatija， Croatia： Faculty of Agriculture， 2005：181-182.

[61]SALEM S， ALGHAMDI H M， MIGDADI M H， et al. Faba bean genomics： current status and future prospects [J]. Euphytica， 2012， 186： 609-624.

[62]REDDY N R， RAGIMASALAWADA M， SABBAVARAPU M M， et al. Detection and validation of stay-green QTL in post-rainy sorghum involving widely adapted cultivar， M35-1 and a popular stay-green genotype B35 [J]. BMC Genomics， 2014， 15 （1）： 909.

[63]SALLAM A， DHANAPAL A， LIU S. Association mapping of winter hardiness and yield traits in winter faba bean （Vicia faba L.） [J]. Crop and Pasture Science， 2015， 67： 55-68.

[64]SALLAM A， ARBAOUI M， ELESAWI M A， et al. Identification and verification of QTL associated with frost tolerance using linkage mapping and GWAS in winter faba bean[J]. Frontiers in Plant Science， 2016， 7： 1098.

[65]XIA M Z. Effects of soil drought during the generative development phase of faba bean（Vicia faba）on photosynthetic characters and biomass production[J]. Journal of Agricultural Science， 1994， 122： 67-72.

[66]KHAZAEI H， OSULLIVAN D M， SILLANP M J， et al. Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean （Vicia faba L.） [J]. Theoretical and Applied Genetics， 2014， 127 （11）： 2371-2385.

[67]WANG S C， WONG D， FORREST K， et al. Characterization of polyploidy wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array [J]. Plant Biotechnology Journal， 2014， 12： 787-796.

[68]王飛飛. 大豆生育期基因TOF7的圖位克隆和功能分析[D].長春：中國科學院大學（中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所）， 2019.

[69]TORRES A M， ROMN B， AVILA C M， et al. Faba bean breeding for resistance against biotic stresses： Towards application of marker technology [J]. Euphytica， 2006， 147 （1/2）： 67-80.

[70]GNANASAMBANDAM A， PAULL J G， TORRES A， et al. Impact of molecular technologies on faba bean （Vicia faba L.） breeding strategies[J]. Agronomy， 2012， 2 （3）：132-166.

[71]BOREVITZ JO， CHORY J. Genomics tools for QTL analysis and gene discovery [J]. Current Opinion in Plant Biology， 2004， 7： 132-136.

[72]RISPAIL N， KAL P， KISS G B， et al. Model legumes contribute to faba bean breeding [J]. Field Crops Research， 2010， 115 （3）：253-269.

（責任編輯：張震林）

收稿日期：2020-07-02

基金項目：國家重點研發(fā)計劃項目（2019YFD1001300、2019YFD1001304）;國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專項（CARS-08）

作者簡介：周仙莉（1995-），女，山東濟寧人，碩士研究生，主要從事蠶豆種質(zhì)資源研究與利用。（E-mail）18797001590@163.com

通訊作者：劉玉皎，（E-mail）13997058356@163.com