潘磊 ，李依 ，余曉露 ，李佳楠 ，陳禪友

（1.江漢大學(xué)生命科學(xué)學(xué)院，武漢，430056；2.湖北省豆類（蔬菜）植物工程技術(shù)研究中心）

豇豆[Vigna unguiculata（L.） Walp.]屬于豆科蝶形花亞科豇豆屬，為一年生草本植物，耐高溫干旱。豇豆是一種重要的豆類作物，在世界范圍內(nèi)種植廣泛，主要種植區(qū)域位于熱帶和亞熱帶的35°N和30°S之間，包括亞洲、大洋洲、中東、歐洲南部、非洲、美國南部和中南美洲（圖1）。主產(chǎn)國為尼日利亞、尼日爾、埃塞俄比亞、突尼斯、中國、印度、菲律賓等。中國的豇豆種植面積常年維持在33萬hm2以上，主要產(chǎn)區(qū)為河北、河南、江蘇、浙江、安徽、四川、重慶、湖北、湖南、廣西等?。ㄗ灾螀^(qū)）。

豇豆屬有3個(gè)栽培亞種，分別是：短豇豆（V.unguiculatassp.cylindrica）、 普 通 豇 豆 （V.unguiculatassp.unguiculata）和長豇豆（V.unguiculatassp.sesquipedalis）。普通豇豆主要分布在非洲的撒哈拉地區(qū)；長豇豆主要分布在中國和印度，而且中國是長豇豆的次生起源中心和多樣性中心。豇豆可供食用的部位為嫩莢和種子，亞洲地區(qū)主要栽培以食用嫩莢為主的長豇豆，而非洲地區(qū)主要栽培以食用種子為主的普通豇豆。根據(jù)聯(lián)合國糧農(nóng)組織（Food and Agriculture Organization，F(xiàn)AO）的統(tǒng)計(jì)，全世界的豇豆種植面積一千多萬公頃，種子常年產(chǎn)量約 600 萬 t（http://faostat3.fao.org/faostat-gateway/go/to/download/Q/QC/E）。由于豇豆富含多種人體所需的植物性蛋白質(zhì)以及膳食纖維，在人們的日常飲食和營養(yǎng)保健中發(fā)揮了重要作用。豇豆除可供鮮食外，還可以進(jìn)行干制、腌制等產(chǎn)后加工，國內(nèi)外市場需求旺盛。

與水稻、小麥、大豆等作物相比，豇豆的分子遺傳學(xué)理論與應(yīng)用基礎(chǔ)研究比較滯后，被稱為“孤兒作物”?？上驳氖牵S著現(xiàn)代生物學(xué)理論與技術(shù)的不斷發(fā)展進(jìn)步，近年來豇豆DNA分子水平上的研究愈來愈成為學(xué)者們關(guān)注的焦點(diǎn)和熱點(diǎn)。下面本文分述豇豆分子遺傳相關(guān)的研究進(jìn)展。

圖1 世界豇豆種植分布主要區(qū)域

1 豇豆的基因庫資源

1.1 豇豆種質(zhì)資源概況

豇豆是二倍體植物（2n=2x=22），基因組大小約為587 Mb[1]；其栽培亞種之一長豇豆（V.unguiculatassp.sesquipedalis）的核型為 10m+1sm[2]。

當(dāng)前，在世界范圍內(nèi)的豇豆種質(zhì)資源的搜集和保存等研究方面已經(jīng)卓有成效（表1）。位于尼日利亞的國際熱帶農(nóng)業(yè)研究所（International Institute for Tropical Agriculture，IITA；http://old.iita.org）建立了世界上最大的豇豆種質(zhì)資源庫，搜集保存了來自89個(gè)國家和地區(qū)的1 507份野生豇豆和15 003份豇豆栽培種，其中對逾12 000份進(jìn)行了28個(gè)農(nóng)藝性狀的評價(jià)，并構(gòu)建了核型種質(zhì)資源庫，共篩選出1 701份地方品種、225份改良的栽培種（包括品種、品系或株系）和130份代表性種質(zhì)材料[3]。美國農(nóng)業(yè)部（United States Department of Agriculture，USDA；http://www.ars-grin.gov）保存了世界各地搜集的6 845份豇豆種質(zhì)資源。位于我國臺灣省的亞洲蔬菜研究開發(fā)中心（The Asian Vegetable Research and Development Center，AVRDC；http://avrdc.org/）保存了1 572份豇豆資源；中國農(nóng)業(yè)科學(xué)院（Academy of Agricultural Sciences，CAAS）國家農(nóng)作物種質(zhì)資源保存中心（http://icgr.caas.net.cn/）保存有1 202份不同類型的豇豆種質(zhì)材料。

1.2 豇豆及其相關(guān)的基因組數(shù)據(jù)庫資源

已建立了豇豆分子生物學(xué)方面的專業(yè)數(shù)據(jù)庫資源，且豆科中其他種屬的植物基因組學(xué)數(shù)據(jù)庫資源也可借鑒與利用（表2）。當(dāng)前，豇豆的數(shù)據(jù)庫主要有4個(gè)，包括NCBI數(shù)據(jù)庫（http://archive-dtd.ncbi.nlm.nih.gov/），豇豆 EST 數(shù)據(jù)庫 HarvEST（http://www.harvest-web.org/）、豇豆甲基化位點(diǎn)相關(guān)的數(shù)據(jù)庫CGKB（http://cowpeagenomics.med.virginia.edu/CGKB）和豇豆物理圖譜數(shù)據(jù)庫Physical Map of Cowpea（http://phymap.ucdavis.edu/cowpea/）。此外，豆科模式植物大豆 （Glycine max）、蒺藜苜 蓿（Medicago truncatula）等建立的基因組學(xué)相關(guān)數(shù)據(jù)庫資源，亦可作為豇豆基因組學(xué)研究的重要參考。

表1 世界范圍內(nèi)豇豆種質(zhì)資源庫概況

2 豇豆的DNA分子標(biāo)記研究

2.1 豇豆DNA分子標(biāo)記的發(fā)掘

當(dāng)前，DNA分子標(biāo)記已經(jīng)廣泛應(yīng)用于植物種質(zhì)資源和作物遺傳育種研究。DNA分子標(biāo)記是遺傳標(biāo)記的一種，是基因組上存在核酸變異的特異性DNA序列，能夠反映個(gè)體或種群之間基因組中存在的某些差異。與形態(tài)標(biāo)記、細(xì)胞標(biāo)記和生化標(biāo)記相比較而言，DNA分子標(biāo)記存在諸多優(yōu)點(diǎn)：①具有較高的遺傳多態(tài)性；②直接以DNA的形式呈現(xiàn)，不受環(huán)境條件、發(fā)育時(shí)期的影響，也沒組織器官特異性；③分布于整個(gè)基因組，而且數(shù)量多；④許多分子標(biāo)記表現(xiàn)為共顯性，能夠區(qū)分基因型中的雜合子；⑤一般表現(xiàn)為選擇中性；⑥經(jīng)濟(jì)方便而且易于檢測分析。

在豇豆的分子標(biāo)記研究中，大多采用基于PCR（Polymerase chain reaction，聚合酶鏈?zhǔn)椒磻?yīng)）的DNA分子標(biāo)記。利用顯性分子標(biāo)記，如RAPD、ISSR等，已經(jīng)篩選獲得一批重復(fù)性好、穩(wěn)定可靠的分子標(biāo)記引物[4，5]。近年來，SSR和SNP等共顯性分子標(biāo)記的應(yīng)用備受青睞。Muchero等[6]2009年首次報(bào)道了基于EST序列的1 375個(gè)SNP位點(diǎn)。Li等[7]2001年率先開發(fā)出27個(gè)SSR標(biāo)記應(yīng)用于豇豆研究。Gupta等[8]2010年從NCBI數(shù)據(jù)庫（http://archive-dtd.ncbi.nlm.nih.gov/）的豇豆unigenes序列中設(shè)計(jì)篩選出102個(gè)SSR標(biāo)記。類似的研究中，Xu等[9]2010年則從豇豆基因組數(shù)據(jù)庫HarvEST（http://www.harvest-web.org/）和 CGKB（http://cowpeagenomics.med.virginia.edu/CGKB）中發(fā)掘出172個(gè)多態(tài)性豇豆SSR分子標(biāo)記（45個(gè)EST-SSR標(biāo)記和127個(gè)gSSR標(biāo)記）。

表2 豇豆及豆科植物基因組數(shù)據(jù)庫資源

2.2 豇豆種質(zhì)資源的DNA分子標(biāo)記應(yīng)用研究

探究豇豆種質(zhì)資源DNA水平的遺傳變異，特別是對野生資源和栽培種遺傳多樣性現(xiàn)狀與親緣關(guān)系的研究，是開展豇豆分子遺傳改良及其應(yīng)用研究的基礎(chǔ)。因此，世界豇豆分布區(qū)的國家和地區(qū)對豇豆種質(zhì)資源開展了不同分子標(biāo)記方面的研究。

①豇豆野生資源的分子標(biāo)記研究 一般認(rèn)為豇豆起源于非洲，因?yàn)橐吧乖诜侵薜貐^(qū)分布廣泛，但是，豇豆在非洲的起源地點(diǎn)具體位于何處尚不清楚，學(xué)者們采用不同DNA分子標(biāo)記進(jìn)行了研究，認(rèn)為非洲東部地區(qū)是野生豇豆起源地，而非洲南部地區(qū)可能是野生豇豆的多樣性中心。Coulibaly等[10]2002年采用AFLP分子標(biāo)記分析了117份豇豆材料，發(fā)現(xiàn)非洲地區(qū)野生豇豆比栽培豇豆的多樣性水平高，推測非洲東部地區(qū)可能是野生豇豆的起源地。與此類似的研究中，Ba等[11]2004年通過對來自非洲東部、西部和南部地區(qū)的栽培豇豆和野生豇豆進(jìn)行RAPD分子標(biāo)記分析，也推測非洲東部可能是野生豇豆的起源地。Ogunkanmi等[12]2008年基于SSR分子標(biāo)記分析非洲野生豇豆的遺傳多樣性，揭示出非洲南部地區(qū)野生豇豆的PIC（Polymorphic information content，多態(tài)信息含量）較高，推測非洲南部地區(qū)可能是野生豇豆的多樣性中心。

②豇豆品種資源的分子標(biāo)記研究 在豇豆品種資源的研究上，采用DNA分子標(biāo)記技術(shù)探究了中國、泰國等國家的豇豆品種間的遺傳多樣性水平和遺傳關(guān)系。徐雁鴻等[13]2007年對來自中國等亞洲國家和部分非洲地區(qū)的316份豇豆品種進(jìn)行SSR分子標(biāo)記分析，發(fā)現(xiàn)中國豇豆品種比外來品種的遺傳多樣性水平低，在中國的豇豆品種遺傳分布不均，廣西和湖北等省份的豇豆具有較高遺傳多樣性，而安徽、吉林、黑龍江和山西等省的則較低；而且基于RAPD分子標(biāo)記的分析結(jié)果表明，我國不同豇豆品種之間遺傳差異較大[4]。在泰國長豇豆遺傳關(guān)系的研究中，通過綜合農(nóng)藝表型、SSR和ISSR分子標(biāo)記等進(jìn)行聚類分析，能夠有效區(qū)分不同的豇豆材料[14]。

在豇豆地方品種的DNA分子標(biāo)記研究中，非洲和亞洲地區(qū)的地方品種比較受關(guān)注，并且在DNA分子標(biāo)記種類的選擇上進(jìn)行了比較分析。Tosti等[15]2002年發(fā)現(xiàn)在豇豆地方品種檢測效率方面，AFLP比RAPD技術(shù)揭示的多樣性指數(shù)更高，尤其是在分析遺傳背景較狹窄的材料中，AFLP技術(shù)更加有效。在AFLP分子標(biāo)記研究方面，劉永華等[16]于2007年構(gòu)建了一套優(yōu)化的豇豆AFLP分析技術(shù)體系。Fang等[17]2007年采用6個(gè)AFLP分子標(biāo)記引物組合分析了來自非洲、亞洲和南美地區(qū)的27個(gè)豇豆地方品種，揭示出這些材料具有較高的遺傳相似性，亞洲和美洲的豇豆可能有共同的起源，而與來自西非的材料不同。此外，在豇豆地方品種的RAPD分子標(biāo)記研究中，研究者發(fā)現(xiàn)馬拉維豇豆地方品種的RAPD標(biāo)記聚類結(jié)果與其形態(tài)特征之間并不一致[18]；而基于RAPD分子標(biāo)記揭示出貝寧地區(qū)豇豆之間存在較大的遺傳差異[19]。豇豆地方品種的SSR分子標(biāo)記研究中，Badiane等[20]2012年利用44個(gè)EST-SSR分子標(biāo)記將塞內(nèi)加爾的22個(gè)豇豆地方品種聚類到一起。Lee等[21]2009年則采用6個(gè)SSR分子標(biāo)記分析了492份韓國豇豆地方品種，揭示出較高的遺傳多樣性水平。

③豇豆種質(zhì)材料的分子標(biāo)記研究 除了豇豆野生種、品種、地方品種外，對不同國家地區(qū)的不同豇豆種質(zhì)材料也開展了DNA分子標(biāo)記的研究。Malviya等[22]2012年采用18個(gè)RAPD標(biāo)記將10份印度豇豆材料分為2個(gè)類群；類似的研究中，Gajera等[23]2014年發(fā)現(xiàn)，基于RAPD和ISSR這兩種分子標(biāo)記能較好地將豇豆品種材料與優(yōu)良基因型材料進(jìn)行聚類區(qū)分。Gillaspie等[24]2005年的研究表明，采用3個(gè)AFLP引物組合和10對SSR引物能有效鑒定出豇豆亞種之間遺傳關(guān)系較近的材料和雜合個(gè)體。Asare等[25]2010年采用SSR分子標(biāo)記評估了來自加納國內(nèi)9個(gè)地區(qū)的141份豇豆材料，可以聚類為5個(gè)類群，每個(gè)類群與其地理來源存在一定的相關(guān)性。Sawadogo等[26]2010年通過SSR分子標(biāo)記篩選出豇豆對寄生雜草（Striga gesnerioides）的抗性材料和敏感型材料。此外，魯忠富等[27]2010年從研究長豇豆的SSR引物中篩選出10對診斷性引物，建立了基于SSR分子標(biāo)記技術(shù)的長豇豆種子純度快速鑒定方法。

近年來，以測序?yàn)榛A(chǔ)的SNP分子標(biāo)記也逐漸應(yīng)用于豇豆研究中。Xu等[28]2012年采用1 127個(gè)基于EST序列的SNP標(biāo)記，分析了包含有長豇豆和普通豇豆的99份核心材料，結(jié)果表明，多態(tài)性SNP的比例相對較低（39%），平均每個(gè)位點(diǎn)1.33個(gè)SNP，Bayesian群體結(jié)構(gòu)分析揭示出2個(gè)亞群，大體上分別與SV亞群（Standard vegetable type）和NSV亞群（Non-standard vegetable type）一致；LD（r2）較高，且LD在亞群SV中比在亞群NSV中持續(xù)得更長，而LD衰減在2個(gè)亞群中均較快，LD衰減在不同染色體中存在差別，最長的是最短的約5倍。該研究首次進(jìn)行了豇豆群體遺傳結(jié)構(gòu)分析，并證明了基于豇豆SNP的全基因組關(guān)聯(lián)分析可以對復(fù)雜性狀進(jìn)行關(guān)聯(lián)作圖。

綜上所述，不同類型DNA分子標(biāo)記已經(jīng)在豇豆種質(zhì)資源研究中成功應(yīng)用（表3），而且以基于PCR技術(shù)的DNA分子標(biāo)記居多，近年來基于測序的SNP分子標(biāo)記也越來越受到青睞。在基于DNA分子標(biāo)記的研究中，認(rèn)為非洲東部地區(qū)是野生豇豆起源地，而非洲南部地區(qū)可能是野生豇豆的多樣性中心，為豇豆的起源和遺傳多樣性中心問題提供了分子依據(jù)。DNA分子標(biāo)記能有效揭示豇豆品種、地方品種等現(xiàn)存豇豆種質(zhì)資源中蘊(yùn)含較為豐富的遺傳變異，而且不同地區(qū)豇豆種質(zhì)資源的遺傳多樣性水平不均一，存在地區(qū)差異，通過聚類可以揭示豇豆個(gè)體之間的遺傳關(guān)系，通常情況下，基于分子標(biāo)記的聚類與性狀特征之間沒有一致性，Schut等[29]認(rèn)為其原因可能是所觀測的性狀、性狀的遺傳變異、性狀相關(guān)的基因數(shù)量等偏少，也可能是基因的上位相互作用。

表3 豇豆研究中采用的幾種主要分子標(biāo)記

表4 成功構(gòu)建的豇豆遺傳連鎖圖譜

2.3 豇豆分子遺傳圖譜研究進(jìn)展

豇豆中分子遺傳圖譜方面的研究雖然起步較晚，但是近年來發(fā)展迅速，已經(jīng)成功構(gòu)建了多個(gè)高密度連鎖遺傳圖（表4），為豇豆的分子遺傳學(xué)研究奠定了良好的基礎(chǔ)。

Ouédraogo等[30]2002年采用133個(gè) RAPD、36個(gè)RFLP和267個(gè)AFLP標(biāo)記，利用F2群體，構(gòu)建了包含11個(gè)連鎖群的豇豆遺傳圖譜，總長2 670.0 cM。2009年，豇豆中第一張高密度SNP連鎖遺傳圖譜成功構(gòu)建[6]，此遺傳圖譜將928個(gè)SNP標(biāo)記整合到11個(gè)連鎖群上，覆蓋了680.0 cM，平均0.73 cM 1個(gè)SNP標(biāo)記；同時(shí)，在與豆科模式植物的大豆和蒺藜苜蓿進(jìn)行基因組同線性分析時(shí)發(fā)現(xiàn)，長豇豆與大豆和蒺藜苜蓿在連鎖遺傳圖上的共線性比率分別為85%和82%。在國內(nèi)的研究中，Xu等[31]2011年利用長豇豆雜交組合（ZN016×ZJ282），構(gòu)建了 1個(gè)覆蓋11個(gè)連鎖群的375個(gè)標(biāo)記位點(diǎn)（其中有191個(gè)SNP標(biāo)記和184個(gè)SSR標(biāo)記），覆蓋基因組745.0 cM，每個(gè)標(biāo)記平均為1.98 cM。Lucas等[32]2011年采用基于EST序列的SNP，對13個(gè)群體的1 293個(gè)個(gè)體用Illumina 1536 GoldenGate Assay分析時(shí)，構(gòu)建了基于 1 107個(gè) SNP的 11個(gè)連鎖群（680.0 cM），并將其與大豆基因組進(jìn)行了共線性比較分析，揭示兩者基因組上存在較高的同源性。Kongjaimun等[33]2012年從豇豆屬的紅豆（V.angularis）和綠豆（V.radiata）基因組中開發(fā)SSR標(biāo)記，利用長豇豆JP81610和野豇豆（V.unguiculatasubsp.unguiculatavar.spontanea）種質(zhì)TVnu457的雜交后代，構(gòu)建了包含226個(gè)SSR標(biāo)記在11個(gè)連鎖群上覆蓋852.4 cM的連鎖圖譜，每個(gè)標(biāo)記間平均3.96 cM。

3 豇豆中功能基因的發(fā)掘與鑒定研究

功能基因的發(fā)掘與鑒定是進(jìn)行分子遺傳改良的前提和基礎(chǔ)。與水稻、玉米、小麥和棉花等農(nóng)作物相比較而言，豇豆功能基因的發(fā)掘與鑒定研究十分薄弱，包括與高溫干旱等非生物逆境抗性相關(guān)的基因，與抗炭腐病（莖枯?。⒏?、枯萎病等生物逆境相關(guān)的基因，以及與葉片、豆莢、花期等生長發(fā)育相關(guān)的基因的研究。

3.1 豇豆轉(zhuǎn)基因研究

由于當(dāng)前豇豆轉(zhuǎn)基因體系尚不成熟，對于細(xì)菌轉(zhuǎn)化體系、電轉(zhuǎn)化體系和基因槍等不同轉(zhuǎn)基因手段的選擇，目標(biāo)組織的遴選等轉(zhuǎn)基因技術(shù)體系中的關(guān)鍵問題仍需不斷深入研究（表5）。近些年，豇豆中已經(jīng)開展nptII基因、CPMV基因、gus基因、hpt基因、αAI-1基因、Cry1Ab基因和Atahas基因等轉(zhuǎn)基因研究[34～50]；最近還報(bào)道了一種利用6-磷酸甘露糖異構(gòu)酶基因作選擇標(biāo)記篩選轉(zhuǎn)基因豇豆的方法[51]。

3.2 耐旱、耐熱相關(guān)基因研究

在耐旱相關(guān)基因的研究中，F(xiàn)ranca等[52]2008年通過克隆表達(dá)豇豆耐旱相關(guān)基因PAP（Phosphatidic acid phosphatase，磷脂酸磷酸酶）的2個(gè)cDNA（Vu-PAP alpha和 VuPAP beta），對 VuPAP alpha啟動子的生物信息學(xué)分析，鑒定出了一些與干旱相關(guān)的調(diào)節(jié)元件。另外，采用深度測序分析豇豆干旱脅迫下表達(dá)的miRNA，鑒定出44個(gè)與干旱相關(guān)的miRNA，其中30個(gè)在干旱中是上調(diào)作用，14個(gè)為下調(diào)作用[53]。Mondal等[54]2014年克隆了豇豆中的β碳酸酐酶基因VuCA1，該基因的表達(dá)具有明顯的組織差異性，在花蕾、莖稈和根部表達(dá)較弱，而在葉片中表達(dá)較強(qiáng)，且在葉片中的基因表達(dá)水平因干旱和鹽脅迫而上升。

表5 豇豆中轉(zhuǎn)基因研究發(fā)展概況

在對豇豆耐旱性相關(guān)基因的QTL定位中，發(fā)掘了一些相關(guān)的QTL位點(diǎn)。采用AFLP分子標(biāo)記，對豇豆耐旱型與敏感型的重組自交系群體進(jìn)行QTL分析，揭示出與苗期干旱耐性和成熟性相關(guān)的12個(gè)QTL[55]；進(jìn)一步的研究中，發(fā)現(xiàn)了7個(gè)AFLP標(biāo)記位于干旱或非生物脅迫誘導(dǎo)相關(guān)的同源區(qū)EST[56]。Lucas等[57]基于SNP分子標(biāo)記，發(fā)掘出與豇豆耐熱相關(guān)的5個(gè)QTL位點(diǎn)Cht-1～Cht-5（表6）。

表6 豇豆中QTL作圖研究

3.3 豇豆抗蟲、抗病相關(guān)基因的研究

在豇豆抗蟲相關(guān)基因的研究方面，豇豆胰蛋白酶抑制劑（Cowpea trypsin inhibitor，CpTI）基因具有抗蟲譜廣而且昆蟲不易對其產(chǎn)生耐受性等特點(diǎn)，Hilder等[58]1987年首先獲得轉(zhuǎn)cpti基因的抗蟲煙草植株；目前，已經(jīng)應(yīng)用在水稻、棉花、紅薯、油棕等多種作物中[59～64]。

Muchero等[65]2011年在豇豆中將抗炭腐病相關(guān)的9個(gè) QTL（Mac-1～Mac-9）位點(diǎn)，定位到 5個(gè)連鎖群 LG2、LG3、LG5、LG6、LG11 上。Pottorff等[66]2012年發(fā)現(xiàn)位于連鎖群LG6上的1個(gè)QTL位點(diǎn)Fot3-1與抗鐮刀菌根腐病相關(guān)。Pottorff等[67]2014年分析出豇豆中與抗鐮刀霉菌枯萎病相關(guān)的2個(gè)QTL位點(diǎn)，F(xiàn)ot4-1和Fot4-2，分別定位到豇豆通用連鎖遺傳圖的 LG5和 LG3上（表6）。

3.4 豇豆發(fā)育相關(guān)基因的研究

豇豆發(fā)育相關(guān)的功能基因研究鮮見報(bào)道，主要涉及戟形葉片、莢長、花期、結(jié)莢、豆莢柔嫩度和可溶性固形物總量、始花時(shí)間、盛花時(shí)間、種皮顏色等相關(guān)基因的QTL研究。

Pandey等[68]2004年研究發(fā)現(xiàn)，豇豆中退化的托葉（Rudimentary stipules，RS）由顯性基因控制，隱性基因控制著葉狀葉托（FS）表型。Pottorff等[69]2012年采用SNP標(biāo)記，將1個(gè)與豇豆戟狀葉片相關(guān)的QTL位點(diǎn)（即Hls基因位點(diǎn)），定位于豇豆第4號連鎖群上，而且發(fā)現(xiàn)1個(gè)與Hls基因區(qū)域共分離的連鎖標(biāo)記SNP 1_0349，序列比對分析發(fā)現(xiàn)，Hls候選基因與擬南芥 EZA1/SWINGER（AT4G02020.1）基因同源性高（表6）。

此外，現(xiàn)有的研究表明，籽粒大小、豆莢大小、莖稈大小和葉片大小等馴化相關(guān)基因的QTL主要位于連鎖群LG7[70]。莢長相關(guān)的7個(gè)QTL[33]，始花期相關(guān)的1個(gè)QTL位點(diǎn)Qfld.zaas-11，始花節(jié)位相關(guān)的1個(gè)QTL位點(diǎn)Qnff.zaas-11，單株結(jié)莢數(shù)相關(guān)的1個(gè)QTL位點(diǎn)Qpn.zaas-3，葉片衰老相關(guān)的一個(gè)QTL位點(diǎn)Qls.zaas-11[71]，豆莢柔嫩度相關(guān)的3個(gè)QTL位點(diǎn)Psn7.1，Psn8.1，Psn11.1，可溶性固形物總量相關(guān)的2個(gè)QTL位點(diǎn)Psw1.1和Psw3.1[72]，始花時(shí)間相關(guān) 的 5 個(gè) QTL 位 點(diǎn)qfot1.1、qfot1.2、qfot1.3、qfot2、qfot10，盛花時(shí)間相關(guān)的3個(gè)QTL位點(diǎn)qdtf1、qdtf2、qdtf7[73]；豇豆種皮褐變相關(guān)的3個(gè)QTL位點(diǎn)Hbs-1、Hbs-2和Hbs-3[74]（表6），上述研究報(bào)道的 QTL位點(diǎn)中，部分位點(diǎn)已經(jīng)整合到豇豆通用遺傳圖譜[7]，包括Hls基因位點(diǎn)，F(xiàn)ot3-1基因位點(diǎn)，Mac-1～Mac-9基因位點(diǎn)，F(xiàn)ot4-1～Fot4-2基因位點(diǎn)，Hbs-1～Hbs-3（圖2）。

4 豇豆的分子進(jìn)化研究

在豇豆的起源、進(jìn)化與分類等的研究中，已經(jīng)獲得一些分子生物學(xué)證據(jù)，較多開展了基于ITS（Internal transcribed spacer，內(nèi)部轉(zhuǎn)錄間隔區(qū)）、IGS（Intergenic spacer，基因間隔區(qū)）等序列變異的系統(tǒng)發(fā)生學(xué)研究，不論是在豇豆屬的層次，或是在種屬之內(nèi)，ITS和IGS序列長度都存在較大遺傳變異，多態(tài)性位點(diǎn)較多，能有效用于豇豆分子鑒定和遺傳系譜分析。

在豇豆屬的分類水平上，Goel等[75]測序比較了豇豆屬29個(gè)物種和菜豆屬9個(gè)物種的rDNA基因中ITS序列（18S-26S rDNA重復(fù)區(qū)），發(fā)現(xiàn)在豇豆屬和菜豆屬中的ITS序列存在差別，ITS1在豇豆屬中的長度在 187～243 bp，而在菜豆屬中為217～290 bp；ITS2在豇豆屬中的長度變化為187～219 bp，而在菜豆屬中為225～243 bp；進(jìn)一步分析發(fā)現(xiàn)，ITS譜系圖與基于形態(tài)學(xué)、生物化學(xué)、細(xì)胞學(xué)和孢粉學(xué)的分類基本一致，而且這些ITS序列可以用于區(qū)分V.mungo、V.radiata、V.umbellata和V.unguiculata的野生種。

在探究豇豆亞屬各物種間的系統(tǒng)發(fā)生學(xué)關(guān)系的研究中，Vijaykumar等[76]2010年基于rRNA基因的ITS區(qū)域，發(fā)現(xiàn)356個(gè)多態(tài)性位點(diǎn)中的80%為簡約信息位點(diǎn)，通過鄰接法（Neighbor joining）和最大簡約法（Maximum parsimony）進(jìn)行遺傳發(fā)育系統(tǒng)關(guān)系重構(gòu)，可以將豇豆亞屬的57個(gè)豇豆材料（屬于15個(gè)種）聚類到5個(gè)主要分支。之后，Vijaykumar等[77]于2011年又分析了豇豆亞屬不同種的豇豆rDNA基因間隔區(qū)（IGS）序列，結(jié)果表明，5S IGS在長度（189～237 bp）和序列（58%多態(tài)性位點(diǎn)）上存在極大變異，而且大多數(shù)的豇豆種都有一種單一類型的5S rRNA重復(fù)單元，但是V.unguiculata和V.reticulata之間是例外，因?yàn)樗鼈冎g表現(xiàn)出多種5S rRNA類型。此外，Doi等[78]研究了豇豆亞屬Ceratotropis的系統(tǒng)發(fā)生學(xué)關(guān)系，比較了rDNA ITS與葉綠體DNA的atpB-rbcL基因間隔區(qū)，2種序列都約 700 bp，rDNA-ITS比atpB-rbcL的多態(tài)性位點(diǎn)多，是atpB-rbcL簡約信息位點(diǎn)（Parsimony-informative sites）的5倍，atpB-rbcL間隔區(qū)更適于對豇豆屬在物種水平上進(jìn)行分析；基于rDNA-ITS聚類的3個(gè)類群，基本與豇豆亞屬Aconitifoliae、Angulares和Ceratotropis相對應(yīng)。

5 結(jié)語

圖2 部分性狀在豇豆通用遺傳圖譜上的QTL定位

當(dāng)前，被稱為“孤兒作物”的豇豆其分子遺傳學(xué)研究發(fā)展迅速，極大地促進(jìn)了豇豆的理論與應(yīng)用基礎(chǔ)研究。在全球范圍內(nèi)已經(jīng)建立了類型豐富的豇豆種質(zhì)資源庫，總份數(shù)超過2.5萬份，而且豇豆及一些豆科模式植物中已建立了專業(yè)數(shù)據(jù)庫資源，為豇豆的基礎(chǔ)理論與應(yīng)用研究奠定了堅(jiān)實(shí)基礎(chǔ)；通過DNA分子標(biāo)記技術(shù)的應(yīng)用，世界豇豆資源的遺傳多樣性總體水平較高，但是存在地域性差異；在豇豆分子連鎖遺傳圖譜研究方面，以往多采用AFLP和SSR分子標(biāo)記，近年來基于測序的SNP分子標(biāo)記越來越受青睞，雖然標(biāo)記數(shù)量較少，不足2 000個(gè)分子標(biāo)記[67]，圖譜標(biāo)記密度不高，但是有效地加深了對目標(biāo)性狀的基因定位研究；當(dāng)前，豇豆中僅開展了少數(shù)目的基因的轉(zhuǎn)基因研究（nptII、CPMV、gus、hpt、αAI-1、Cry1Ab和Atahas等基因），獲得了轉(zhuǎn)基因的愈傷組織或者單株，但是豇豆的轉(zhuǎn)基因體系仍不完善，尚未廣泛應(yīng)用，需要構(gòu)建高效穩(wěn)定的轉(zhuǎn)基因體系；此外，豇豆的分子進(jìn)化研究為傳統(tǒng)的起源、進(jìn)化與分類等提供了分子水平的佐證。

盡管在豇豆的分子遺傳學(xué)研究方面取得了上述進(jìn)展，加深了豇豆的分子基礎(chǔ)理論與應(yīng)用研究，但是仍然有許多問題值得我們持續(xù)關(guān)注，尤其是在現(xiàn)代生物技術(shù)大發(fā)展的背景下，豇豆分子水平的研究面臨新的機(jī)遇與挑戰(zhàn)：①雖然已在世界范圍內(nèi)開展豇豆種質(zhì)資源的交流與合作研究，實(shí)現(xiàn)了豇豆資源的合理有效利用，進(jìn)行了豇豆種質(zhì)創(chuàng)新和遺傳改良，但是由于我國豇豆種質(zhì)資源的收集、鑒定和評價(jià)等研究發(fā)展緩慢，在國家種質(zhì)資源庫中的豇豆數(shù)量偏少，亟待充實(shí)更新，且需要建立健全豇豆公共數(shù)據(jù)庫；②豇豆DNA分子標(biāo)記需要大量的開發(fā)及應(yīng)用，尤其是高密度分子遺傳連鎖圖譜的構(gòu)建，與生物逆境（蟲害、病害等）、非生物逆境（高溫、干旱、鹽堿、濕澇等）以及與正常發(fā)育相關(guān)的功能基因鑒定與發(fā)掘；③完善豇豆的轉(zhuǎn)基因系統(tǒng)，通過豇豆的轉(zhuǎn)基因研究，闡明豇豆基因的功能與作用；④隨著基因組學(xué)、轉(zhuǎn)錄組學(xué)、蛋白組學(xué)和代謝組學(xué)等蓬勃發(fā)展，“組學(xué)”將在多個(gè)層次上拓寬和加深豇豆的分子遺傳機(jī)理闡釋。

[1]Iwata A，Greenland C M，Jackson S A.Cytogenetics of legumes in the Phaseoloid clade[J].Plant Genome，2013，8.

[2]She C W，Jiang X H，Ou L J，et al.Molecular cytogenetic characterisation and phylogenetic analysis of the seven cultivatedVignaspecies(Fabaceae)[J].Plant Biology，2014.DOI:10.1111/plb.12174(online).

[3]Mahalakshmi V，Ng Q，Lawson M，et al.Cowpea[Vigna unguiculata(L.) Walp.]core collection defined by geographical，agronomical and botanical descriptors[J].Plant Genet Resour，2007，5:113-119.

[4]陳禪友，潘磊，胡志輝，等.長豇豆品種資源的RAPD分析[J].江漢大學(xué)學(xué)報(bào)：自然科學(xué)版，2008（36）：77-83.

[5]Chen C Y，Peng H.Differentiation and phylogenetic relationship among different cultivars of asparagus bean[Vignaunguiculata(L.)ssp.sesquipedalis]assessed using ISSR markers[J].Nord J Bot，2010，28:251-256.

[6]Muchero W，Diop N N，Bhat P R，et al.A consensus genetic map of cowpea[Vigna unguiculata(L.)Walp.]and synteny based on EST-derived SNPs[J].Proc Natl Acad Sci USA，2009，106:18 159-18 164.

[7]Li C D，F(xiàn)atokun C A，Ubi B，et al.Determining genetic similarities and relationships among cowpea breeding lines and cultivars by microsatellite markers[J].Crop Sci，2001，41:189-197.

[8]Gupta S K，Gopalakrishna T.Development of unigenederived SSR markers in cowpea (Vigna unguiculata)and their transferability to otherVignaspecies[J].Genome，2010，53:508-523.

[9]Xu P，Wu X，Wang B，et al.Development and polymorphism ofVigna unguiculatassp.unguiculatamicrosatellite markers used for phylogenetic analysis in asparagus bean[Vigna unguiculatassp.sesquipedialis(L.)Verdc.][J].Mol Breeding，2010，25:675-684.

[10]Coulibaly S，Pasquet R S，Papa R，et al.AFLP analysis of the phenetic organization and genetic diversity ofVigna unguiculataL.Walp.reveals extensive gene flow between wild and domesticated types[J].Theor Appl Genet，2002，104:358-366.

[11]Ba F S，Pasquet R S，Gepts P.Genetic diversity in cowpea[Vigna unguiculata(L.)Walp.]as revealed by RAPD markers[J].Genet Resour Crop Ev，2004，51:539-550.

[12]Ogunkanmi L A，Ogundipe O T，Ng N Q，et al.Genetic diversity in wild relatives of cowpea(Vigna unguiculata)as revealed by simple sequence repeats(SSR)markers[J].J Food Agric Environ，2008，6:263-268.

[13]徐雁鴻，關(guān)建平，宗緒曉.豇豆種質(zhì)資源SSR標(biāo)記遺傳多樣性分析[J].作物學(xué)報(bào)，2007（33）：1 206-1 209.

[14]Tantasawat P，Trongchuen J，Prajongjai T，et al.Variety identification and comparative analysis of genetic diversity in yardlong bean (Vigna unguiculataspp.sesquipedalis)using morphological characters，SSR and ISSR analysis[J].Sci Hortic，2010，124:204-216.

[15]Tosti N，Negri V.Efficiency of three PCR-based markers in assessing genetic variation among cowpea (Vigna unguiculatasubsp.unguiculata)landraces[J].Genome，2002，45:268-275.

[16]劉永華，李國景，吳曉花，等.長豇豆AFLP技術(shù)體系的構(gòu)建與優(yōu)化[J].浙江農(nóng)業(yè)學(xué)報(bào)，2007（19）：156-159.

[17]Fang J G，Chao C T，Roberts P A，et al.Genetic diversity of cowpea[Vigna unguiculata(L.)Walp.]in four West African and USA breeding programs as determined by AFLP analysis[J].Genet Resour Crop Ev，2007，54:1 197-1 209.

[18]Nkongolo K K.Genetic characterization of Malawian cowpea[Vigna unguiculata(L.)Walp.]landraces:diversity and gene flow among accessions[J].Euphytica，2003，129:219-228.

[19]Zannou A，Kossou D K，Ahanchede A，et al.Genetic variability of cultivated cowpea in Benin assessed by random amplified polymorphic DNA[J].Afr J Biotechnol，2008，7:4 407-4 414.

[20]Badiane F A，Gowda B S，Cissé N，et al.Genetic relationship of cowpea (Vigna unguiculata)varieties from Senegal based on SSR markers[J].Genet Mol Res，2012，11:292-304.

[21]Lee J R，Back H J，Yoon M S，et al.Analysis of genetic diversity of cowpealandraces from Korea determined by simple sequence repeats and establishment of a core collection[J].Kor J Breed Sci，2009，41:369-376.

[22]Malviya N，Sarangi B K，Yadav M K，et al.Analysis of genetic diversity in cowpea (Vigna unguiculataL.Walp.)cultivars with random amplified polymorphic DNA markers[J].Plant Syst Evol，2012，298:523-526.

[23]Gajera H P，Domadiya R K，Patel S V，et al.Appraisal of RAPD and ISSR markers for genetic diversity analysis among cowpea(Vigna unguiculataL.)genotypes[J].J Crop Sci Biotech，2014，17:79-88.

[24]Gillaspie A G，Hopkins M S，Dean R E.Determining genetic diversity between lines ofVigna unguiculatasubspecies by AFLP and SSR markers[J].Genet Resour Crop Ev，2005，52:245-247.

[25]Asare A T，Gowda B S，Galyuon I K A，et al.Assessment of the geneticdiversity in cowpea (Vigna unguiculataL.Walp.)germplasm from Ghana using simple sequence repeat markers[J].Plant Genet Resour，2010，8:142-150.

[26]Sawadogo M，Ouedraogo J T，Gowda B S，et al.Genetic diversity of cowpea(Vigna unguiculataL.Walp.)cultivars in Burkina Faso resistant toStriga gesnerioide[J].Afr J Biotechnol，2010，9:8 146-8 153.

[27]魯忠富，徐沛，汪寶根，等.基于SSR分子標(biāo)記技術(shù)的長豇豆種子純度快速鑒定技術(shù)[J].浙江農(nóng)業(yè)學(xué)報(bào)，2010，22：727-730.

[28]Xu P，Wu X，Wang B，et al.Genome wide linkage disequilibrium in Chinese asparagus bean (Vigna unguiculatassp.sesquipedialis)germplasm implications for domestication history and genome wide association studies[J].Heredity，2012，109:34-40.

[29]Schut J W，Qi X，Stam P.Association between relationship measures based on AFLP markers，pedigree data and morphological traits in barley[J].Theor Appl Genet，1997，95:1 161-1 168.

[30]Ouédraogo J T，Gowda B S，Jean M，et al.An improved genetic linkage map for cowpea (Vigna unguiculataL.)combining AFLP，RFLP，RAPD，biochemical markers，and biological resistance traits[J].Genome，2002，45:175-188.

[31]Xu P，Wu X，Wang B，et al.A SNP and SSR based genetic map of asparagus bean (Vigna unguiculatassp.sesquipedialis)and comparison with the broader species[J].PLoS One，2011，6:e15952.

[32]Lucas M R，Diop N N，Wanamaker S，et al.Cowpea-soybean synteny clarified through an improved genetic map[J].Plant Genome，2011，4:218-225.

[33]Kongjaimun A，Kaga A，Tomooka N，et al.An SSR-based linkage map of yardlong bean(Vigna unguiculataL.Walp.subsp.unguiculataSesquipedalis group)and QTL analysis of pod length[J].Genome，2012，55:81-92.

[34]Garcia J A，Hille J，Goldbach R.Transformation of cowpea(Vigna unguiculata)cells with an antibiotic-resistance gene using a Ti-plasmid-derived vector[J].Plant Sci，1986，44:37-46.

[35]Garcia J A，Hille J，Vos P，et al.Transformation of cowpea(Vigna unguiculata)with a full-length DNA copy of cowpea mosaic virus M-RNA[J].Plant Sci，1987，48:89-98.

[36]Penza R，Lurquin P F，F(xiàn)ilippone E.Gene transfer by cocultivation of mature embryos withAgrobacterium tumefaciens:application to cowpea(Vigna unguiculataL.Walp.)[J].J Plant Physiol，1991，138:39-43.

[37]Penza R，Akella V，Lurquin P F.Transient expression and histological localization of aguschimeric gene after direct transfer to mature cowpea embryos[J].Biotechniques，1992，13:576-580.

[38]Akella V，Lurquin P F.Expression in cowpea seedlings of chimeric transgenes after electroporation into seed derived embryos[J].Plant Cell Rep，1993，12:110-117.

[39]Muthukumar B，Mariamma M，Veluthambi K，et al.Genetic transformation of cotyledonexplants of cowpea (Vigna unguiculataL.Walp)usingAgrobacterium tumefaciens[J].Plant Cell Rep，1996，15:980-985.

[40]Ikea J，Ingelbrecht I，Uwaifo A，et al.Stablegene transformation in cowpea(Vigna unguiculataL.Walp.)using particle gun method[J].Afr J Biotechnol，2003，2:211-218.

[41]Popelka J C，Gollasch S，Moore A，et al.Genetic transformation of cowpea(Vigna unguiculataL.)and stable transmission of the transgenes to progeny[J].Plant Cell Rep，2006，25:304-312.

[42]Chaudhury D，Madanpotra S，Jaiwal R，et al.Agrobacterium tumefaciens-mediated high frequency genetic transformation of an Indian cowpea (Vigna unguiculataL.Walp.)cultivar and transmission of transgenes into progeny[J].Plant Sci，2007，172:692-700.

[43]Solleti S K，Bakshi S，Purkayastha J，et al.Transgenic cowpea (Vigna unguiculata)seeds expressing a bean αamylase inhibitor 1 confers resistance tostorage pests，bruchid beetles[J].Plant Cell Rep，2008，27:1 841-1 850.

[44]Ivo N L，Nascimento C P，Vieira L S，et al.Biolistic-mediated genetic transformation of cowpea(Vigna unguiculata)and stable Mendelian inheritance of transgenes[J].Plant Cell Rep，2008，27:1 475-1 483.

[45]Umaharan P，Ariyanayagan R P，Haque S Q.Resistance to cowpea severe mosaic virus，determined by three dosage dependent genes inVigna unguiculataL.Walp.[J].Euphytica，1997，95:49-55.

[46]Adesoye A，Machuka J，Togun A.CRY 1ABtransgenic cowpea obtained by nodal electroporation[J].Afr J Biotechnol，2008，7:3 200-3 210.

[47]Adesoye A，Machuka J，Togun A.Transformation of cowpea(Vigna unguiculataL.Walp.)byAgrobacteriuminfiltration[J].J Appl Biosci，2010，30:1 845-1 860.

[48]Raveendar S，Ignacimuthu S.ImprovedAgrobacteriummediated transformation in cowpeaVigna unguiculataL.Walp.[J].Asian J Plant Sci，2010，9:256-263.

[49]Citadin C T，Cruz A R R，Aragao F J L.Development of transgenic imazapyr-tolerant cowpea (Vigna unguiculata)[J].Plant Cell Rep，2013，32:537-543.

[50]Citadin C T，Ibrahim A B，Aragao F J L.Genetic engineering in cowpea (Vigna unguiculata):history，status and prospects[J].GM Crops，2011，2:144-149.

[51]Bakshi S，Saha B，Roy N K，et al.Successful recovery of transgenic cowpea (Vigna unguiculata)using the 6-phosphomannose isomerase gene as the selectable marker[J].Plant Cell Rep，2012，31:1 093-1 103.

[52]Franca M G C，Matos A R，D'Arcy-Lameta A，et al.Cloning and characterization of drought-stimulated phosphatidic acid phosphatase genes fromVigna unguiculata[J].Plant Physiol Biochem，2008，6:1 093-1 100.

[53]Barrera-Figueroa B E，Gao L，Diop N N，et al.Identification and comparative analysis of drought-associated microRNAs in two cowpea genotypes[J].BMC Plant Biol，2011，11:127.

[54]Mondal T K，Ganie S A，Niraj R R K，et al.Cloning and in silico analysis of a gene encoding a putative β-carbonic anhydrase from cowpea (Vigna unguiculataL.Walp.)[J].J Plant Interact，2014，9:504-513.

[55]Muchero W，Ehlers J D，Close T J，et al.Mapping QTL for drought stress-induced premature senescence and maturity in cowpea[Vigna unguiculata(L.)Walp.][J].Theor Appl Genet，2009，118:849-863.

[56]Muchero W，Ehlers J D，Roberts P A.Restriction site polymorphism-based candidate gene mapping for seedling drought tolerance in cowpea[Vigna unguiculata(L.)Walp.][J].Theor Appl Genet，2010，120:509-518.

[57]Lucas M R，Ehlers J D，Huynh B L，et al.Markers for breeding heat-tolerant cowpea[J].Mol Breeding，2013，31:529-536.

[58]Hilder V A，Gatehouse A M，Sheerman S E，et al.A novel mechanism of insect resistance engineered into tobacco[J].Nature，1987，330:160-163.

[59]劉春明，朱禎，周兆斕，等.豇豆胰蛋白酶抑制劑抗蟲轉(zhuǎn)基因煙草的獲得[J].科學(xué)通報(bào)，1992（37）：1 694-1 697.

[60]NewellC A，LoweJM，MerryweatherA，etal.Transformation of sweet potato(Ipomoea batatasL.Lam.)withAgrobacterium tumefaciensand regeneration of plants expressing cowpea trypsin inhibitor and snowdrop lectin[J].Plant Sci，1995，107:215-227.

[61]楊朝輝，宋洪元，何鳳田.豇豆胰蛋白酶抑制劑基因轉(zhuǎn)化芥菜及抗蟲鑒定[J].中國生物化學(xué)與分子生物學(xué)報(bào)，2003（19）：731-735.

[62]徐鴻林，翟紅利，王鋒，等.豇豆胰蛋白酶抑制劑基因（cpti）及其在抗蟲轉(zhuǎn)基因作物中的應(yīng)用[J].中國農(nóng)業(yè)科技導(dǎo)報(bào)，2008（10）：18-27.

[63]Ismanizan I，F(xiàn)ong S L，Ruslan A，et al.Molecular and expression analysis of cowpea trypsin inhibitor(CpTI)gene in transgenicElaeis guineensisJacq leaves[J].Aust J Crop Sci，2010，4:37-48.

[64]Wu H，Zhang Y，Liu P，et al.Effects of transgenicCry1Ac+CpTIcotton on non-target mealybug pestFerrisia virgataand its predatorCryptolaemus montrouzieri[J].PloS One，2014，9:e95537.

[65]Muchero W，Ehlers J D，Close T J，et al.Genic SNP markers and legume synteny reveal candidate genes underlying QTL forMacrophomina phaseolinaresistance and maturity in cowpea[Vigna unguiculata(L.)Walp.][J].BMC Genomics，2011，12:8-21.

[66]Pottorff M，Wanamaker S，Ma Y Q，et al.Genetic and physical mapping of candidate genes for resistance toFusarium oxysporumf.sp.tracheiphilumrace 3 in cowpea[Vigna unguiculata(L.)Walp.][J].PloS One，2012，7:e41600.

[67]Pottorff M，Li G，Ehlers J D，et al.Genetic mapping，synteny，and physical location of two loci forFusarium oxysporumf.sp.tracheiphilumrace 4 resistance in cowpea[Vigna unguiculata(L.)Walp.][J].Mol Breeding，2014，33:779-791.

[68]Pandey R N，Dhanasekar P.Morphological features and inheritance of foliaceous stipules of primary leaves in cowpea(Vigna unguiculata)[J].Ann Bot，2004，94:469-471.

[69]Pottorff M，Ehlers J D，F(xiàn)atokun C，et al.Leaf morphology in cowpea[Vigna unguiculata(L.)Walp.]:QTL analysis，physical mapping and identifying a candidate gene using synteny with model legume species[J].BMC Genomics，2012，13:234-245.

[70]Kongjaimun A，Kaga A，Tomooka N，et al.The genetics of domestication of yardlong bean，Vigna unguiculata(L.)Walp.ssp.unguiculatacv.-gr.sesquipedalis[J].Ann Bot，2012，109:1 185-1 200.

[71]Xu P，Wu X，Wang B，et al.QTL mapping and epistatic interaction analysis in asparagus bean for several characterized and novel horticulturally important traits[J].BMC Genet，2013，14:4-13.

[72]Kongjaimun A，Somta P，Tomooka N，et al.QTL mapping of pod tenderness and total soluble solid in yardlong bean[Vigna unguiculata(L.)Walp.subsp.unguiculatacv.-gr.sesquipedalis][J].Euphytica，2013，189:217-223.

[73]Andargie M，Pasquet R S，Muluvi G M，et al.Quantitative trait loci analysis of flowering time related traits identified in recombinant inbred lines of cowpea(Vigna unguiculata)[J].Genome，2013，56:289-294.

[74]Pottorff M，Roberts P A，Timothy J，et al.Identification of candidate genes and molecular markers for heat-induced brown discoloration of seed coats in cowpea[Vigna unguiculata(L.)Walp.][J].BMC Genomics，2014，15:328-338.

[75]Goel S，Raina S N，Ogihara Y.Molecular evolution and phylogenetic implications of internal transcribed spacer sequences of nuclear ribosomal DNA in thePhaseolus-Vignacomplex[J].Mol Phylogenet Evol，2002，22:1-19.

[76]Vijaykumar A，Saini A，Jawali N.Molecular characterization of intergenic spacer region of 5S ribosomal RNA genes in subgenusVigna:extensive hybridization amongV.unguiculatasubspecies[J].Plant Syst Evol，2011，294:39-55.

[77]Vijaykumar A，Saini A，Jawali N.Phylogenetic analysis of subgenusVignaspecies using nuclear ribosomal RNA ITS:evidence of hybridization amongVigna unguiculatasubspecies[J].J Hered，2010，101:177-188.

[78]Doi K，Kaga A，Tomooka N，et al.Molecular phylogeny of genusVignasubgenusCeratotropisbased on rDNA ITS andatpB-rbcL intergenic spacer of cpDNA sequences[J].Genetica，2002，114:129-145.