李曉川 王朝海 周平 馬維 陳軍 陸燚 吳顯 王宗明 吳瑞 宋治豪 馬杰 付毅

摘要 從馬鈴薯參考基因組中篩選出了433個(gè)NBS-LRR類(lèi)基因，并繪制這些基因在基因組中的分布圖，這些基因的其中一部分成簇分布。馬鈴薯基因組測(cè)序協(xié)會(huì)（PGSC）的轉(zhuǎn)錄組數(shù)據(jù)分析顯示，只有12個(gè)NBS-LRR類(lèi)基因在PGSC的轉(zhuǎn)錄組數(shù)據(jù)中完全沒(méi)有檢測(cè)到表達(dá)，不同NBS-LRR類(lèi)基因轉(zhuǎn)錄水平不同，一部分NBS-LRR類(lèi)基因有較高的生物脅迫和激素刺激應(yīng)激性。分析了這433個(gè)基因的編碼區(qū)域的單核苷酸多態(tài)性（SNP）位點(diǎn)，SNP標(biāo)記作為在高通量測(cè)序和生物大數(shù)據(jù)分析中最常使用的遺傳標(biāo)記，有助于鑒定NBS-LRR基因在不同馬鈴薯品種中的實(shí)際抗病性。

關(guān)鍵詞 馬鈴薯；晚疫?。豢剐曰?/p>

中圖分類(lèi)號(hào) Q37? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 0517-6611（2023）10-0079-06

doi：10.3969/j.issn.0517-6611.2023.10.017

Abstract In this paper， 433 NBSLRR genes were screened from the Solanum tuberosum reference genome， and the distribution of these genes in the genome were mapped， showing that some of these genes were distributed in clusters. Only 12 NBSLRR genes were not detected in all of the transcriptome analysis data of PGSC. The expression levels of different NBSLRR genes were different， and some NBSLRR genes had higher biological and hormonal stimulation irritability. The single nucleotide polymorphism （SNP） sites in the coding regions of these 433 genes were also analyzed. As the most commonly used genetic markers in highthroughput sequencing and biological big data analysis， SNP markers are helpful for the identification of actual disease resistance of NBSLRR genes in Solanum tuberosum varieties.

Key words Solanum tuberosum;Phytophthora infestans （Mont.） DE Bary;Resistance （R） gene

馬鈴薯（Solanum tuberosum）適應(yīng)性廣，豐產(chǎn)性好，營(yíng)養(yǎng)豐富，經(jīng)濟(jì)效益高，已成為世界上繼水稻、小麥的第三大糧食作物［1］。但在生產(chǎn)中經(jīng)常受到其第一大病害——晚疫病的危害造成產(chǎn)量嚴(yán)重降低。晚疫病是由致病疫霉菌［Phytophthora infestans（Mont.）de Bary.］引起，特別是由于病原菌A2交配型的出現(xiàn)，與原A1交配型有性雜交，發(fā)生變異，導(dǎo)致最初的抗病馬鈴薯品種失去抗性［2-3］。

目前已經(jīng)發(fā)現(xiàn)多個(gè)馬鈴薯晩疫病抗病基因（R基因），包括克隆自Solanum demissum中，位于第5號(hào)染色體上的R1［4］，第4號(hào)染色體上的R2［5］，第11號(hào)染色體上的R3a［6］和R3b［7］，第9號(hào)染色體上的R8［8］；克隆自Solanum bulbocastanum，位于第8號(hào)染色體上的Rpi-blb1（RB）［9-10］，第6號(hào)染色體上的Rpi-blb2［11］，第9號(hào)染色體上的Rpi-blb3和第4號(hào)染色體上的Rpi-abpt［5］；來(lái)源于Solanum venturii位于第9號(hào)染色體上Rpi-vnt1.1、Rpi-vnt1.2和Rpi-vnt1.3［12］；來(lái)源于Solanum stoloniferum位于第8號(hào)染色體上Rpi-sto1和Rpi-pta1［13-14］；來(lái)源于Solanum americanum位于第4號(hào)染色體上Rpi-amr3i［15］。已鑒定出的R基因都擁有NBS-LRR（Nucleotidebinding site and leucinerich repeat）結(jié)構(gòu)，因此可以利用基因結(jié)構(gòu)上的保守型從參考基因組中篩選潛在的抗性基因［16-17］。目前已在多個(gè)植物的參考基因組中篩選出NBS-LRR類(lèi)基因，包括重要的作物水稻［18］、小麥［19］、玉米［20］、大豆［21］、棉花［22］、十字花科［23］、黃瓜［24］、模式生物-擬南芥［25］以及茄科植物辣椒、番茄［26］和茄子［27］等。雖然馬鈴薯晩疫病抗病基因能提供抗病性，但單個(gè)R基因?qū)霂?lái)的抗性易被變異的致病疫霉菌克服，如將RD基因?qū)攵傻脑耘喾NBiogold，在田間種植不到1年即被晚疫病菌克服［28］，但通常認(rèn)為引入多個(gè)R基因能夠增強(qiáng)晚疫病抗病表現(xiàn)［29］。因此，盡可能多地篩選出晩疫病抗病基因，有助于馬鈴薯晚疫病抗性育種。

高通量測(cè)序技術(shù)提高了人們對(duì)作物遺傳學(xué)的理解。由于普通栽培馬鈴薯基因組是高度雜合的同源4倍體的特性。馬鈴薯基因組測(cè)序協(xié)會(huì)（potato genome sequencing consortium，PGSC）首先通過(guò)對(duì)一個(gè)純合的雙單倍體（DM1-3 516 R44，DM）馬鈴薯品種測(cè)序得到了馬鈴薯參考基因組，參考基因組提供了詳細(xì)的馬鈴薯基因座信息，因此可以通過(guò)序列比對(duì)在參考基因組中查找具有保守序列和結(jié)構(gòu)域的基因［30］。并且利用參考基因組對(duì)一個(gè)（雜合二倍體RH89-039-16，RH）馬鈴薯品種重測(cè)序，基因分型得到了兩者之間的單核苷酸多態(tài)性位點(diǎn)（SNP）。PGSC還利用參考基因組對(duì)DM的16種組織包括：心皮、花瓣、萼片、雄蕊、花、果肉、未成熟果實(shí)、成熟果實(shí)、葉柄、葉、匍匐莖、未成熟塊莖、成熟塊莖、愈傷組織、根和芽的轉(zhuǎn)錄組，以及11種經(jīng)處理包括：鹽、甘露醇和熱脅迫，4種激素［IAA（吲哚乙酸）、BAP（6-芐基腺嘌呤）、GA3（赤霉酸）和ABA（脫落酸）］處理，4種生物脅迫［接種晚疫病菌、2種誘導(dǎo)植物抗病的因子苯并噻二唑（BTH）和DL-氨基丁酸（BABA）處理以及受傷葉片］后轉(zhuǎn)錄組表達(dá)水平的變化情況處理進(jìn)行了轉(zhuǎn)錄組測(cè)序。并且作為對(duì)DM轉(zhuǎn)錄組數(shù)據(jù)的補(bǔ)充和對(duì)照，PGSC還對(duì)雜合二倍體品種RH的花、雄蕊、根、葉、水分脅迫葉、葉柄、莖尖、莖、匍匐莖、未成熟塊莖、成熟塊莖、塊莖髓、塊莖皮、塊莖芽、塊莖皮質(zhì)和整株16種組織的轉(zhuǎn)錄組進(jìn)行了測(cè)序。轉(zhuǎn)錄組測(cè)序數(shù)據(jù)經(jīng)隨后的RNA定量驗(yàn)證顯示出很高的準(zhǔn)確性［30-31］。因此，可以利用R基因的NBS-LRR結(jié)構(gòu)特性，分析馬鈴薯參考基因組種的NBS-LRR類(lèi)基因，并通過(guò)轉(zhuǎn)錄組數(shù)據(jù)分析這些基因在不同組織和處理?xiàng)l件下的表達(dá)特性，以及利用基因組序列信息，分析分子標(biāo)記信息，為鑒定各基因的實(shí)際抗病性鑒定打下理論基礎(chǔ)。

1 材料與方法

1.1 馬鈴薯NBS-LRR類(lèi)基因的鑒定

使用BLASTp工具（http：//spuddb.uga.edu/blast.shtml），以期望值（E-value=1）為閾值，用已知的R基因的蛋白序列作為查詢(xún)序列，以及隱馬模型（Hidden Markov Model，HMM）［32］，以Pfam：PF00931作為查詢(xún)序號(hào)，以默認(rèn)參數(shù)從馬鈴薯參考基因組（V 4.03）中鑒定（http：//spuddb.uga.edu/integrated_searches.shtml）。利用SMART檢索基因的蛋白結(jié)構(gòu)域（http：//smart.emblheidelberg.de/）［33］；以0.9為閾值，利用 COILS 程序檢索基因的CC結(jié)構(gòu)域（http：//www.ch.embnet.org/software/COILS_form.html）［34］?；虻陌被嵝蛄欣肅lustal Omega在線服務(wù)器進(jìn)行多序列比對(duì)，利用Neighbor-Joining參數(shù)進(jìn)過(guò)計(jì)算1 000次生成樹(shù)圖（http：//www.ebi.ac.uk/Tools/msa/clustalo/）。

1.2 馬鈴薯NBS-LRR類(lèi)基因的表達(dá)分析

利用馬鈴薯基因組測(cè)序協(xié)會(huì)（potato genome sequencing consortium，PGSC）的Illumina RNA 測(cè)序數(shù)據(jù)分析NBS-LRR類(lèi)基因在不同組織和處理方式下的表達(dá)模式［30-31］。轉(zhuǎn)錄組數(shù)據(jù)的表達(dá)豐度用每百萬(wàn)個(gè)映射上的片段數(shù)中映射到外顯子中的每1 000個(gè)堿基上的片段個(gè)數(shù)（fragments per kilobase of exon model per million mapped reads，F(xiàn)PKM）值表示。對(duì)不同組織內(nèi)NBS-LRR類(lèi)基因表達(dá)水平的熱圖（Heat map）進(jìn)行構(gòu)建時(shí)，F(xiàn)PKM 值首先經(jīng)過(guò)以2為底數(shù)的對(duì)數(shù)進(jìn)行轉(zhuǎn)化。構(gòu)建生物脅迫以及激素處理時(shí)NBS-LRR類(lèi)基因表達(dá)水平變化的熱圖時(shí)，首先計(jì)算與對(duì)照的對(duì)比值，以log2比值的相對(duì)值≥1確定為差異表達(dá)基因［31］。熱圖使用 MeV 4.9.0 軟件進(jìn)行構(gòu)建［35］。相關(guān)性分析用IBM SPSS Statistics 22 軟件進(jìn)行（SPSS GmbH Software，Munich，Germany）。

1.3 馬鈴薯NBS-LRR類(lèi)基因編碼序列的SNP標(biāo)記分析

NBS-LRR類(lèi)基因編內(nèi)的SNP標(biāo)記信息來(lái)源于以DM為參考基因組，RH重測(cè)序信息，Sol-CAP對(duì)于6個(gè)常見(jiàn)栽培種的轉(zhuǎn)錄組測(cè)序信息，以及對(duì)大西洋（Atlantic）和Superior的重測(cè)序信息［6，36］。

2 結(jié)果與分析

2.1 馬鈴薯NBS-LRR類(lèi)基因的鑒定

從馬鈴薯參考基因組中鑒定得到433個(gè)NBS-LRR類(lèi)基因，其中366個(gè)基因能夠定位到12個(gè)染色體中。在這些基因中的142個(gè)基因間的距離小于1.0 mb，其中44個(gè)基因間的距離小于10 kb。21個(gè)基因定位于第4染色體上2.0 mb的染色體區(qū)域，12個(gè)基因定位于第5染色體上1.5 mb的染色體區(qū)域，22個(gè)基因定位于第11染色體上2.0 mb的染色體區(qū)域，說(shuō)明這些基因的其中一部分成簇分布（圖1）。433個(gè)NBS-LRR類(lèi)基因中的103個(gè)基因蛋白編碼序列的N端含有一個(gè)CC（coiled-coil）結(jié)構(gòu)域，這與已發(fā)現(xiàn)的晚疫病抗性基因的CC-NBS-LRR蛋白結(jié)構(gòu)域相同（圖2）。進(jìn)一步利用已知的R基因和CC-NBS-LRR基因的氨基酸序列進(jìn)行比對(duì)并繪制系統(tǒng)進(jìn)化樹(shù)。進(jìn)化樹(shù)分成2個(gè)部分（黑色和紅色），其中黑色部分沒(méi)有已克隆的R基因。紅色部分的進(jìn)化樹(shù)又可分成2個(gè)部分（基因編號(hào)突出顯示和未突出顯示的部分），這2部分內(nèi)部又可分成2～3個(gè)小組（基因編號(hào)不同顏色字體顯示）。Rpi-blb1（RB）與15個(gè)CC-NBS-LRR基因聚集成一個(gè)小組，R3a和R3b 與19個(gè)CC-NBS-LRR基因親緣關(guān)系較近，Rpi-amr3i與14個(gè)CC-NBS-LRR基因聚集成一個(gè)小組，Rpi-vnt1.1、Rpi-vnt1.2和Rpi-vnt1.3與17個(gè)CC-NBS-LRR基因聚集成一個(gè)小組，其他7個(gè)已克隆的R基因與23個(gè)CC-NBS-LRR基因親緣關(guān)系較近（圖3）。

2.2 馬鈴薯NBS-LRR類(lèi)基因的表達(dá)分析

利用PGSC在馬鈴薯不同器官和發(fā)育階段以及在脅迫處理和激素刺激時(shí)的不同轉(zhuǎn)錄組測(cè)序所產(chǎn)生的數(shù)據(jù)，包括：DM的16種組織內(nèi)基因的表達(dá)水平，DM的11種包括鹽脅迫、激素刺激、生物脅迫處理后轉(zhuǎn)錄組表達(dá)水平的變化情況，以及以DM為參考基因組，RH的16種組織內(nèi)基因的表達(dá)模式。分析NBS-LRR類(lèi)基因的表達(dá)。轉(zhuǎn)錄組數(shù)據(jù)分析顯示，只有12個(gè)NBS-LRR類(lèi)基因在所有的組織或處理中完全沒(méi)有檢測(cè)到表達(dá)（圖2）。

PGSC分別在DM和RH的16種不同組織中進(jìn)行了轉(zhuǎn)錄組分析，提取了其中易受晚疫病侵染影響的組織的轉(zhuǎn)錄組表達(dá)數(shù)據(jù)，進(jìn)一步分析了433個(gè)NBS-LRR類(lèi)基因在不同組織和樣品間的表達(dá)模式。在DM的6種組織中（根、葉、葉柄、匍匐莖、不成熟塊莖和成熟塊莖），433個(gè)NBS-LRR類(lèi)基因，表達(dá)量FPKM均值最大的10個(gè)基因分別是PGSC0003DMG400007999、PGSC0003DMG400009635、PGSC0003DMG400007469、PGSC0003DMG400003527、PGSC0003DMG402018257、PGSC0003DMG400018464、PGSC0003DMG400031405、PGSC0003DMG400020587、PGSC0003DMG400033159和PGSC0003DMG400007471，其表達(dá)量的FPKM值在52.67～11.60。有57個(gè)基因在DM的6種組織中沒(méi)有檢測(cè)到表達(dá)，208個(gè)基因在這6種組織中至少有1種組織中沒(méi)有檢測(cè)到表達(dá)。376個(gè)基因在至少DM的6種組織之一中檢測(cè)到表達(dá)，反映表達(dá)量在6種組織間離散程度的標(biāo)準(zhǔn)差（σ），最大為34.13，平均為1.76。在RH的8種組織中（根、葉、葉柄、植株、莖、匍匐莖、不成熟塊莖和成熟塊莖），433個(gè)之間NBS-LRR類(lèi)基因，表達(dá)量均值最大的10個(gè)基因分別是PGSC0003DMG400007999、PGSC0003DMG400007469、PGSC0003DMG400031405、PGSC0003DMG400013490、PGSC0003DMG400029505、PGSC0003DMG400005335、PGSC0003DMG400007462、PGSC0003DMG400018428、PGSC0003DMG400018429和PGSC0003DMG400026187，其表達(dá)量的FPKM值在60.27～12.58。有32個(gè)基因在RH的8種組織中沒(méi)有檢測(cè)到表達(dá)，179個(gè)基因在這8種組織中至少有1種組織中沒(méi)有檢測(cè)到表達(dá)。401個(gè)基因在至少RH的8種組織之一檢測(cè)到表達(dá)，反映表達(dá)量FPKM值在8種組織間離散程度的標(biāo)準(zhǔn)差（σ），最大為62.45，平均為1.57。結(jié)果顯示，相同的NBS-LRR類(lèi)基因在DM或RH的不同組織間的表達(dá)水平有一定的差異，但整體上差異并不明顯，即組織特異性不明顯。在DM和RH組織之間，相同基因表達(dá)量的FPKM值均值，其中在DM中大于RH中的基因有126個(gè)，在RH中大于DM中的基因有283個(gè)，說(shuō)明相同的NBS-LRR類(lèi)基因在DM和RH之間，整體上基因在RH中表達(dá)水平略高于在DM中。433個(gè)NBS-LRR基因在DM和RH組織間表達(dá)水平的Pearson相關(guān)系數(shù)r=0.600（P<0.01），有顯著的正相關(guān)關(guān)系，說(shuō)明NBS-LRR基因在DM和RH的組織中有類(lèi)似的表達(dá)模式。

生物脅迫處理包括接種馬鈴薯晚疫病菌和2種誘導(dǎo)植物抗病的因子［苯并噻二唑（BTH，100 mg/mL）以及DL-氨基丁酸（BABA，2 mg/mL）］處理24 h/48 h/72 h后混合樣本的表達(dá)豐度。參加分析的數(shù)據(jù)為計(jì)算每個(gè)處理相對(duì)于對(duì)照的相對(duì)表達(dá)豐度。從結(jié)果中可知，在接種晚疫病菌處理后，有32個(gè)基因表達(dá)量顯著上升（基因表達(dá)量log2比值≥1，即處理后的基因表達(dá)量是對(duì)照的基因表達(dá)量的2倍及以上）；4個(gè)基因（PGSC0003DMG400011521、PGSC0003DMG400011906、PGSC0003DMG400033161和PGSC0003DMG400009307）的基因表達(dá)量log2比值>2。在接受誘導(dǎo)植物抗病的BTH因子處理后，有36個(gè)基因表達(dá)量顯著上升，3個(gè)基因（PGSC0003DMG400011906、PGSC0003DMG400010527和PGSC0003DMG400044423）的表達(dá)量log2比值>2。在BABA因子處理后，有179個(gè)基因表達(dá)量顯著上升，32個(gè)基因的表達(dá)量log2比值>2，3個(gè)基因（PGSC0003DMG401007609、PGSC0003DMG403015877和PGSC0003DMG400028081）的表達(dá)量log2比值>4（圖2）。在3種生物脅迫方式處理后，有4個(gè)基因（PGSC0003DMG400011521、PGSC0003DMG400011906、PGSC0003DMG400033161和PGSC0003DMG400044423）至少在接受2種處理后基因表達(dá)量log2比值>2。

有研究顯示，晚疫病抗性基因的表達(dá)受到激素的影響［37］。為研究激素對(duì)NBS-LRR類(lèi)基因的表達(dá)的影響，分析整株植物經(jīng)IAA（吲哚乙酸）10 mmol /L、BAP（6-芐基腺嘌呤）10 mmol /L、GA3（赤霉酸）50 mmol /L、ABA（脫落酸）50 mmol /L 處理24 h后的相對(duì)表達(dá)豐度。在IAA刺激后，有80個(gè)基因表達(dá)量顯著上升（基因表達(dá)量log2比值≥1，即處理后的基因表達(dá)量是對(duì)照的基因表達(dá)量的2倍及以上），4個(gè)基因（PGSC0003DMG400020741、PGSC0003DMG400011521、PGSC0003DMG400029220和PGSC0003DMG400001989）的基因表達(dá)量log2比值>2。在BAP處理后，有66個(gè)基因表達(dá)量顯著上升，15個(gè)基因表達(dá)量log2比值>2，2個(gè)基因（PGSC0003DMG400010612和PGSC0003DMG400026433）的表達(dá)量log2比值>4。在GA3刺激后，有159個(gè)基因表達(dá)量顯著上升，17個(gè)基因表達(dá)量log2比值>2，2個(gè)基因（PGSC0003DMG 400020741和PGSC0003DMG400019804）的表達(dá)量log2比值>4。在ABA處理后，有51個(gè)基因表達(dá)量顯著上升，2個(gè)基因（PGSC0003DMG400020741和PGSC0003DMG400011521）的表達(dá)量log2比值>2。在4種激素刺激后，有8個(gè)基因（PGSC0003DMG400001989、PGSC0003DMG400002427、PGSC0003DMG400003527、PGSC0003DMG400011521、PGSC0003DMG400019804、PGSC0003DMG400020587、PGSC0003DMG400020741和PGSC0003DMG400026433）至少在接受2種處理后基因表達(dá)量log2比值>2（圖2）。

2.3 馬鈴薯NBS-LRR類(lèi)基因編碼序列的SNP標(biāo)記分析

SNP標(biāo)記是DNA中的單堿基變異造成的，在該研究中，以雙單倍體DM為參考基因組，利用雜合二倍體RH重測(cè)序信息，Sol-CAP對(duì)于6個(gè)常見(jiàn)栽培種的轉(zhuǎn)錄組測(cè)序信息，以及對(duì)栽培種大西洋（Atlantic）和Superior的重測(cè)序信息，分析在433個(gè)NBS-LRR類(lèi)基因編碼序列內(nèi)的SNP標(biāo)記。結(jié)果顯示，在以上9個(gè)品種的序列信息內(nèi)，只有19個(gè)NBS-LRR類(lèi)基因的編碼區(qū)域內(nèi)沒(méi)有SNP，134個(gè)NBS-LRR類(lèi)基因的編碼區(qū)域內(nèi)有1～10個(gè)SNP，67個(gè)NBS-LRR類(lèi)基因的編碼區(qū)域內(nèi)有超過(guò)50個(gè)SNP（圖2）。同時(shí)，還分析了NBS-LRR類(lèi)基因編碼序列內(nèi)SNP的數(shù)量是否影響其表達(dá)，結(jié)果顯示沒(méi)有相關(guān)性。

3 討論與結(jié)論

目前已發(fā)現(xiàn)的晚疫病抗病（R）基因都是NBS-LRR家族基因，利用已測(cè)序完成的參考基因組序列信息，可以初步篩選潛在的抗性基因。目前已在多個(gè)植物基因組中鑒定出NBS-LRR家族候選基因，包括重要的作物水稻［12］、小麥［13］、玉米［14］、大豆［15］、棉花［16］、十字花科［17］、黃瓜［18］以及模式生物-擬南芥［19］。在茄科植物中，從辣椒的參考基因組中鑒定到了305個(gè) NBS-LRR類(lèi)基因，從番茄的參考基因組中鑒定到了255個(gè)NBS-LRR類(lèi)基因［20］，從茄子的參考基因組中鑒定到了245個(gè)NBS-LRR類(lèi)基因［21］。該研究從馬鈴薯參考基因組中鑒定出433個(gè)NBS-LRR類(lèi)基因（圖1）。相較于其他茄科植物，從馬鈴薯參考基因組中鑒定出更多的NBS-LRR類(lèi)基因，可能是由于更多的基因組倍性變化，如普通栽培種馬鈴薯的基因組是四倍體，而野生馬鈴薯的基因組是二倍體至六倍體，導(dǎo)致的基因組更多變異的累加，并分化出不同的NBS-LRR類(lèi)基因。這些NBS-LRR類(lèi)基因的其中一部分成簇分布（圖1），成簇分布的基因，由于連鎖效應(yīng)，將更易于在雜交育種的過(guò)程中，一次性轉(zhuǎn)移更多的NBS-LRR類(lèi)基因進(jìn)行品種改良，而選育晚疫病抗性馬鈴薯品種需要更多的R基因累加。

基因發(fā)揮功能首先必須轉(zhuǎn)錄成mRNA，在該研究中利用PGSC的轉(zhuǎn)錄組數(shù)據(jù)，包括DM的16種組織內(nèi)基因的表達(dá)水平，DM對(duì)11種包括鹽脅迫、激素刺激、生物脅迫處理后轉(zhuǎn)錄組表達(dá)水平的變化情況，以及以DM為參考基因組，RH的16種組織內(nèi)基因的表達(dá)模式［30-31］。轉(zhuǎn)錄組數(shù)據(jù)分析顯示，只有12個(gè)NBS-LRR類(lèi)基因在所有的組織或處理中完全沒(méi)有檢測(cè)到表達(dá)。進(jìn)一步分析轉(zhuǎn)錄組數(shù)據(jù)，在易受晚疫病侵染影響的DM的6種組織和RH的8種組織中NBS-LRR基因的表達(dá)情況中，不同的NBS-LRR基因表達(dá)水平不同，但相同的NBS-LRR基因在DM或RH的不同組織間的表達(dá)水平有一定差異，但整體上差異并不明顯，即組織特異性不明顯。相同的NBS-LRR基因整體上在RH中表達(dá)水平略高于在DM中的表達(dá)水平。433個(gè)NBS-LRR基因在DM和RH組織間表達(dá)水平有顯著的正相關(guān)關(guān)系，說(shuō)明NBS-LRR基因在DM和RH中表達(dá)模式類(lèi)似。同時(shí)，部分NBS-LRR基因在生物脅迫處理和激素刺激后，表達(dá)水平顯著上升，有明顯的應(yīng)激性。

SNP標(biāo)記作為在高通量測(cè)序和生物大數(shù)據(jù)分析中最常使用的遺傳標(biāo)記，在該研究中分析了433個(gè)NBS-LRR基因的編碼區(qū)域內(nèi)的SNP位點(diǎn)（圖2），已知的SNP標(biāo)記，可以利用多態(tài)性位點(diǎn)及側(cè)翼序列，通過(guò)如KASP（Kompetitive Allele Specific PCR）高通量基因分型技術(shù)，實(shí)現(xiàn)在不同品種中對(duì)該位點(diǎn)較快速的基因分型［38］，有助于鑒定NBS-LRR基因在不同馬鈴薯品種中的實(shí)際抗病性。根據(jù)選擇性清除理論，在受選擇的基因組區(qū)域，其多態(tài)性會(huì)降低，表現(xiàn)為SNP標(biāo)記數(shù)目減少［39］，盡管研究結(jié)果顯示NBS-LRR基因編碼區(qū)域內(nèi)的SNP數(shù)量的多少與其表達(dá)不相關(guān)，但67個(gè)NBS-LRR基因的編碼區(qū)域內(nèi)存在50個(gè)以上的SNP，可能代表基因組的這部分區(qū)域重組頻繁，對(duì)基因的功能也可能產(chǎn)生一定的影響。

參考文獻(xiàn)

［1］ BHASKAR P B，VENKATESHWARAN M，WU L，et al.Agrobacteriummediated transient gene expression and silencing：A rapid tool for functional gene assay in potato［J］.PLoS One，2009，4（6）：1-8.

［2］ FRY W E，GOODWIN S B.Resurgence of the Irish potato famine fungus［J］.Bioscience，1997，47：363-371.

［3］ FRY W.Phytophthora infestans：The plant（and R gene）destroyer［J］.Molecular plant pathology，2008，9（3）：385-402.

［4］ BALLVORA A，ERCOLANO M R，WEISS J，et al.The R1 gene for potato resistance to late blight（Phytophthora infestans）belongs to the leucine zipper /NBS /LRR class of plant resistance genes［J］.Plant journal，2002，30（3）：361-371.

［5］ LOKOSSOU A A，PARK T H，VAN ARKEL G，et al.Exploiting knowledge of R/Avr genes to rapidly clone a new LZNBSLRR family of late blight resistance genes from potato linkage group IV［J］.Molecular plantmicrobe interactions，2009，22（6）：630-641.

［6］ HUANG S W，VAN DER VOSSEN E A，KUANG H H，et al.Comparative genomics enabled the isolation of the R3a late blight resistance gene in potato［J］.Plant journal，2005，42（2）：251-261.

［7］ LI G C，HUANG S W，GUO X，et al.Cloning and characterization of R3b；members of the R3 superfamily of late blight resistance genes show sequence and functional divergence［J］.Molecular plantmicrobe interactions，2011，24（10）：1132-1142.

［8］ VOSSEN J H，VAN ARKEL G，BERGERVOET M，et al.The Solanum demissum R8 late blight resistance gene is an Sw5 homologue that has been deployed worldwide in late blight resistant varieties［J］.Theoretical and applied genetics，2016，129（9）：1785-1796.

［9］ VAN DER VOSSEN E，SIKKEMA A，HEKKERT B T，et al.An ancient R gene from the wild potato species Solanum bulbocastanum confers broadspectrum resistance to Phytophthora infestans in cultivated potato and tomato［J］.Plant journal，2003，36（6）：867-882.

［10］ SONG J Q，BRADEEN J M，NAESS S K，et al.Gene RB cloned from Solanum bulbocastanum confers broad spectrum resistance to potato late blight［J］.Proc Natl Acad Sci USA，2003，100（16）：9128-9133.

［11］ VAN DER VOSSEN E A，GROS J，SIKKEMA A，et al.The Rpiblb2 gene from Solanum bulbocastanum is an Mi1 gene homolog conferring broadspectrum late blight resistance in potato［J］.Plant journal，2005，44（2）：208-222.

［12］ FOSTER S J，PARK T H，PEL M，et al.Rpivnt1.1，a Tm22 homolog from Solanum venturii，confers resistance to potato late blight［J］.Molecular plantmicrobe interactions，2009，22（5）：589-600.

［13］ VLEESHOUWERS V G A A，RIETMAN H，KRENEK P，et al.Effector genomics accelerates discovery and functional profiling of potato disease resistance and Phytophthora infestans avirulence genes［J］.PloS One，2008，3（8）：1-10.

［14］ WANG M Q，ALLEFS S，VAN DEN BERG R G，et al.Allele mining in Solanum：Conserved homologues of Rpiblb1 are identified in Solanum stoloniferum［J］.Theoretical and applied genetics，2008，116（7）：933-943.

［15］ WITEK K，JUPE F，WITEK A I，et al.Accelerated cloning of a potato late blightresistance gene using RenSeq and SMRT sequencing［J］.Nature biotechnol，2016，34（6）：656-660.

［16］ MARTIN G B，BOGDANOVE A J，SESSA G.Understanding the functions of plant disease resistance proteins［J］.Annual review of plant biology，2003，54：23-61.

［17］ VAN DER BIEZEN E A，JONES J D.The NBARC domain：A novel signalling motif shared by plant resistance gene products and regulators of cell death in animals［J］.Current biology，1998，8（7）：R226-R227.

［18］ LIU W，GHOURI F，YU H，et al.Genome wide resequencing of newly developed Rice Lines from common wild rice（Oryza rufipogon Griff.）for the identification of NBSLRR genes［J］.PLoS One，2017，12（7）：1-19.

［19］ GU L J，SI W N，ZHAO L N，et al.Dynamic evolution of NBSLRR genes in bread wheat and its progenitors［J］.Molecular genetics genomics，2015，290（2）：727-738.

［20］ CHENG Y，LI X Y，JIANG H Y，et al.Systematic analysis and comparison of nucleotidebinding site disease resistance genes in maize［J］.FEBS journal，2012，279（13）：2431-2443.

［21］ KANG Y J，KIM K H，SHIM S，et al.Genomewide mapping of NBSLRR genes and their association with disease resistance in soybean［J］.BMC plant biology，2012，12：1-13.

［22］ KHAN A M，KHAN A A，AZHAR M T，et al.Comparative analysis of resistance gene analogues encoding NBSLRR domains in cotton［J］.Journal of the science of food and agriculture，2016，96（2）：530-538.

［23］ ZHANG Y M，SHAO Z Q，WANG Q，et al.Uncovering the dynamic evolution of nucleotidebinding siteleucinerich repeat（NBSLRR）genes in Brassicaceae［J］.Journal of integrative plant biology，2016，58（2）：165-177.

［24］ WAN H J，YUAN W，BO K J，et al.Genomewide analysis of NBSencoding disease resistance genes in Cucumis sativus and phylogenetic study of NBSencoding genes in Cucurbitaceae crops［J］.BMC genomics，2013，14：1-15.

［25］ MONDRAGNPALOMINO M，MEYERS B C，MICHELMORE R W，et al.Patterns of positive selection in the complete NBSLRR gene family of Arabidopsis thaliana［J］.Genome research，2002，12（9）：1305-1315.

［26］ SHAO Z Q，XUE J Y，WU P，et al.Largescale analyses of angiosperm nucleotidebinding siteleucinerich repeat genes reveal three anciently diverged classes with distinct evolutionary patterns［J］.Plant physiology，2016，170（4）：2095-2109.

［27］ WEI H W，LIU J，GUO Q W，et al.Genomic organization and comparative phylogenic analysis of NBSLRR resistance gene family in Solanum pimpinellifolium and Arabidopsis thaliana［J/OL］.Evolutionary bioinformatics，2020，16［2022-03-17］.https：//doi.org/10.1176934320911055.

［28］ HERMSEN J G TH，RAMANNA M S.Doublebridge hybrids of Solanum bulbocastanum and cultivars of Solanum tuberosum［J］.Euphytica，1973，22（3）：457-466.

［29］ HAVERKORT A J，BOONEKAMP P M，HUTTEN R，et al.Durable late blight resistance in potato through dynamic varieties obtained by cisgenesis：Scientific and societal advances in the DuRPh project［J］.Potato research，2016，59（1）：35-66.

［30］The Potato Genome Sequencing Consortium.Genome sequence and analysis of the tuber crop potato［J］.Nature，2011，475（7355）：189-195.

［31］ MASSA A N，CHILDS K L，LIN H N，et al.The transcriptome of the reference potato genome Solanum tuberosum Group Phureja clone DM13516R44［J］.PLoS One，2011，6（10）：1-8．

［32］ ELGEBALI S，MISTRY J，BATEMAN A，et al.The Pfam protein families database in 2019［J］.Nucleic acids research，2019，47：D427-D432.

［33］ LETUNIC I，KHEDKAR S，BORK P.SMART：Recent updates，new developments and status in 2020［J］.Nucleic acids research，2021，49（D1）：D458-D460.

［34］ LUPAS A，VAN DYKE M，STOCK J.Predicting coiled coils from protein sequences［J］.Science，1991，252（5009）：1162-1164.

［35］ SAEED A I，SHAROV V，WHITE J，et al.TM4：A free，opensource system for microarray data management and analysis［J］.Biotechniques，2003，34（2）：374-378．

［36］ HAMILTON J P，HANSEY C N，WHITTY B R，et al.Single nucleotide polymorphism discovery in elite North American potato germplasm［J］.BMC genomics，2011，12：1-12.

［37］ TURNBULL D，YANG L N，NAQVI S，et al.RXLR effector AVR2 upregulates a brassinosteroidresponsive bHLH transcription factor to suppress immunity［J］.Plant physiology，2017，174（1）：356-369.

［38］ ALVAREZFERNANDEZ A，BERNAL M J，F(xiàn)RADEJAS I，et al.KASP：A genotyping method to rapid identification of resistance in Plasmodium falciparum［J］.Malar journal，2021，20（1）：1-8.

［39］ NIELSEN R，WILLIAMSON S，KIM Y，et al.Genomic scans for selective sweeps using SNP data［J］.Genome research，2005，15（11）：1566-1575.