仲巖,張?zhí)?陳釗,程宗明

(南京農(nóng)業(yè)大學(xué)園藝學(xué)院,江蘇 南京 210095)

森林草莓(Fragariavesca)屬于薔薇科草莓屬的二倍體種,分布廣泛,具有較強(qiáng)的抗逆性和抗病性[1]。植物通過(guò)遺傳變異和自然選擇分化成不同的生態(tài)型,這些生態(tài)型的分化同生長(zhǎng)環(huán)境密切相關(guān),并且種群分布越廣,積累遺傳變異的概率也會(huì)越高[2]。森林草莓在歐洲各地均有分布,分化為不同的生態(tài)型。歐洲森林草莓不同生態(tài)型的抗病基因之間也存在差異,這些差異可能導(dǎo)致了森林草莓不同的生態(tài)型抗病能力。

抗病基因是植物抵抗自然界中病原物的重要方式,其中NBS-LRR類基因是最為重要的一類抗病基因,占已發(fā)現(xiàn)抗病基因的80%,且具有高遺傳變異度[3]。NBS-LRR基因家族的蛋白質(zhì)結(jié)構(gòu)特點(diǎn)是包含有N端結(jié)構(gòu)域、核苷酸結(jié)合位點(diǎn)(nucleotide-binding-site,NBS)和富亮氨酸重復(fù)區(qū)域(leucine-rich repeat,LRR)。根據(jù)其蛋白質(zhì)N端的TIR(toll-Interleukin-1receptor)、CC(coiled-coil)、RPW8或X(其他)結(jié)構(gòu)域,NBS-LRR基因家族可以進(jìn)一步分為TIR-NBS-LRR(TNL)基因和non-TIR-NBS-LRR(non-TNL,包括CC-NBS-LRR,CNL;RPW8-NBS-LRR,RNL和X-NBS-LRR,XNL)基因2大類型。近年來(lái),眾多園藝植物的全基因組NBS-LRR基因家族已被鑒定,如擬南芥[4]、木薯[5]、向日葵[6]、蘭花[7]、中華獼猴桃[8]、蘋果、梨、桃和梅[9-10]等,其中也包括多個(gè)草莓屬物種[11]。Zhong等[11]對(duì)森林草莓基因組中的NBS-LRR基因研究發(fā)現(xiàn),其呈現(xiàn)出快速進(jìn)化和慢速進(jìn)化2種不同的進(jìn)化模式。由于森林草莓分布較廣,分化出多種生態(tài)型,為了對(duì)這些不同生態(tài)型間的NBS-LRR基因遺傳變異情況和進(jìn)化模式有進(jìn)一步了解,根據(jù)我們構(gòu)建的森林草莓NBS-LRR基因家族的系統(tǒng)進(jìn)化樹(shù)[11],發(fā)現(xiàn)FvH4_1g15230和FvH4_5g22710基因均位于支長(zhǎng)較短的進(jìn)化支上,說(shuō)明它們可能處于快速進(jìn)化模式下,因此在森林草莓中選取了這2個(gè)具有代表性的NBS-LRR基因(其中FvH4_1g15230基因位于森林草莓1號(hào)染色體上,FvH4_5g22710基因位于5號(hào)染色體上)進(jìn)行遺傳變異分析。本研究通過(guò)分析這2個(gè)NBS-LRR基因在不同森林草莓生態(tài)型中完整核苷酸序列的理化性質(zhì)、核苷酸序列以及蛋白序列比對(duì)、蛋白質(zhì)結(jié)構(gòu)域,同時(shí)根據(jù)不同森林草莓生態(tài)型核苷酸序列差異構(gòu)建系統(tǒng)發(fā)育樹(shù)并計(jì)算Ka/Ks值及核苷酸差異度(Pi值),為草莓抗病新品種的培育提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

供試材料取自芬蘭、挪威、西班牙、英國(guó)、冰島、意大利的23種森林草莓生態(tài)型(表1)。

表1 23種森林草莓生態(tài)型的采樣國(guó)家及編號(hào)Table 1 Sampling countries and numbers of 23 woodland strawberry ecotypes

1.2 試驗(yàn)方法

1.2.1 2個(gè)NBS-LRR基因在不同森林草莓生態(tài)型中的核苷酸序列選取健康且長(zhǎng)勢(shì)一致的不同森林草莓生態(tài)型,并以FvH4_5g22710和FvH4_5g22710為目的基因。利用Plant Total DNA Isolation Kit Plus試劑盒(FOREGENE公司)提取不同生態(tài)型的基因組DNA。根據(jù)目標(biāo)基因在Phytozome v12.1(https://phytozome.jgi.doe.gov/pz/portal.html)上的全長(zhǎng)序列,采用Oligo6軟件進(jìn)行引物設(shè)計(jì)。采用PrimeSTAR?GXL DNA Polymerase(TaKaRa公司)擴(kuò)增目的基因全長(zhǎng)序列,并用10 g·L-1瓊脂糖凝膠電泳檢測(cè)目標(biāo)基因條帶,利用OMEGA膠回收試劑盒回收目的片段,送上海生工生物工程公司進(jìn)行測(cè)序。利用已公布全基因組序列的森林草莓‘Hawaii4’作為參照,使用Sequencher軟件將測(cè)序結(jié)果進(jìn)行拼接、核對(duì),最終確定不同生態(tài)型中目的基因的序列。

1.2.2 理化性質(zhì)分析及分類利用ProtParam網(wǎng)站(https://web.expasy.org/protparam/)分析目的基因編碼蛋白序列的氨基酸數(shù)量、相對(duì)分子質(zhì)量、對(duì)脂溶指數(shù)、理論等電點(diǎn)、總平均親水性等理化性質(zhì)。采用Pfam網(wǎng)站(http://pfam.xfam.org/)、SMART網(wǎng)站(http://smart.embl-heidelberg.de/)和COILS網(wǎng)站(http://embnet.vital-it.ch/software/COILS_form.html)分析其蛋白質(zhì)結(jié)構(gòu)域組成。

1.2.3 不同生態(tài)型中目標(biāo)基因及其編碼蛋白質(zhì)的序列對(duì)比及系統(tǒng)進(jìn)化樹(shù)構(gòu)建利用MEGA v7.0軟件[12]中的ClustalW程序?qū)ν换虿煌鷳B(tài)型的核苷酸及基因編碼的蛋白質(zhì)序列進(jìn)行排序,并利用DNAMEN v7.0軟件將比對(duì)結(jié)果可視化。根據(jù)上述獲得的核苷酸序列比對(duì)結(jié)果,利用MEGA v7.0軟件采用最大似然法(Bootstrap值設(shè)置為1 000,Tamura-Nei模型)對(duì)2個(gè)目的基因分別構(gòu)建系統(tǒng)進(jìn)化樹(shù)。

1.2.4 同義突變率(Ks)、非同義突變率(Ka)和核苷酸差異度(Pi值)的計(jì)算利用MEGA v7.0軟件計(jì)算不同生態(tài)型之間目的基因CDS(coding sequence)序列的Ka、Ks、Ka/Ks以及Pi值。利用PAML程序[13]通過(guò)模型0和模型7、模型8計(jì)算ω值(dN/dS,即Ka/Ks)并進(jìn)行顯著性檢測(cè)。

2 結(jié)果與分析

2.1 2個(gè)NBS-LRR基因的理化性質(zhì)分析

供試材料取自芬蘭、挪威、西班牙、英國(guó)、冰島、意大利的23種森林草莓生態(tài)型,其中FvH4_5g22710基因存在于21種森林草莓生態(tài)型中;FvH4_1g15230基因存在于20種森林草莓生態(tài)型中。

由表2可見(jiàn):在21種森林草莓生態(tài)型中,FvH4_5g22710基因所編碼蛋白的理化性質(zhì)差異較小。它們均由941個(gè)氨基酸殘基組成,但由于某幾個(gè)位點(diǎn)的氨基酸種類的差異,其蛋白相對(duì)分子質(zhì)量和理論等電點(diǎn)略有差異。其中,來(lái)自ES5、ES9、ES18、UK1等7個(gè)生態(tài)型的蛋白相對(duì)分子質(zhì)量為108 121.04,理論等電點(diǎn)為8.65;相對(duì)分子質(zhì)量為108 122.03的蛋白來(lái)自其他14種生態(tài)型,包括‘Hawaii4’(FvH4)、FIN11、FIN29、FIN5、FIN6、IT18、IT3、NOR122、NOR21、NOR13、ICE4、NOR150、NOR163和ICE16,其理論等電點(diǎn)均為8.59;因此,這21種生態(tài)型FvH4_5g22710基因編碼的蛋白均為堿性蛋白。該基因編碼蛋白的對(duì)脂溶指數(shù)和總平均親水性在這21種生態(tài)型中均沒(méi)有表現(xiàn)出差異,前者均為103.15,后者均為-0.205,說(shuō)明該基因編碼的蛋白均為親水蛋白。

表2 FvH4_5g22710基因在21種森林草莓生態(tài)型中的理化性質(zhì)Table 2 Physicochemical properties of FvH4_5g22710 gene from 21 woodland strawberry ecotypes

由表3可見(jiàn):在20種森林草莓生態(tài)型中,FvH4_1g15230基因所編碼蛋白的氨基酸數(shù)量相同,均由 817個(gè)氨基酸殘基組成,但理化性質(zhì)也表現(xiàn)出了差異?；蚓幋a的氨基酸種類差異導(dǎo)致了蛋白相對(duì)分子質(zhì)量的差異,其中來(lái)自ICE5、ICE12、ICE16、NOR163生態(tài)型的蛋白相對(duì)分子質(zhì)量最大,為93 262.11,來(lái)自生態(tài)型FIN11的蛋白相對(duì)分子質(zhì)量最小,為93 057.94。在20種生態(tài)型中,該基因編碼的蛋白理論等電點(diǎn)為5.46～5.50,表明其均為酸性蛋白;對(duì)脂溶指數(shù)為101.37和101.73;總平均親水性是-0.121～-0.109,說(shuō)明該基因編碼的蛋白均為親水蛋白。

表3 FvH4_1g15230基因在20種森林草莓生態(tài)型中的理化性質(zhì)Table 3 Physicochemical properties of FvH4_1g15230 gene from 20 woodland strawberry ecotypes

2.2 2個(gè)NBS-LRR基因的序列對(duì)比

來(lái)自21種生態(tài)型的FvH4_5g22710基因,其編碼區(qū)堿基序列的變異程度較小,僅有2處點(diǎn)突變。以FvH4的核苷酸序列為參考,來(lái)自ES5、ES9、ES18、UK1、UK9、UK10和IT24這7個(gè)生態(tài)型的FvH4_5g22710基因,在核苷酸編碼區(qū)的第157位由堿基G突變?yōu)閴A基C,該突變?yōu)榉峭x突變,導(dǎo)致其蛋白序列的第53位由谷氨酸突變?yōu)楣劝滨０?圖1);該基因在生態(tài)型ICE4、ICE16、NOR163和NOR150中,其編碼區(qū)的第207位由堿基C突變?yōu)閴A基T,該突變?yōu)橥x突變,仍舊編碼丙氨酸。該基因的內(nèi)含子區(qū)域并無(wú)突變發(fā)生。

圖1 FvH4_5g22710基因編碼區(qū)(A)及編碼的蛋白(B)在21種森林草莓生態(tài)型中突變位點(diǎn)序列比對(duì)Fig.1 Alignment of FvH4_5g22710 coding sequences(A)and protein sequences encoded by FvH4_5g22710(B) around the mutation sites in 21 ecotypes of woodland strawberry 圖中省略號(hào)代表無(wú)突變發(fā)生的核苷酸序列(A)或氨基酸序列(B)。The ellipsis means the nucleotide sequences(A)or amino acid sequences(B)without mutation.

以來(lái)自FvH4的核苷酸序列為參考,在FvH4_1g15230基因的編碼區(qū)發(fā)生了17處點(diǎn)突變,包括11處非同義突變和6處同義突變(圖2)。例如,來(lái)自ICE5生態(tài)型的該基因在其編碼區(qū)第463位發(fā)生了非同義突變,由堿基T突變?yōu)镚,因此,其蛋白序列的第155位由半胱氨酸變?yōu)楦拾彼?圖3);NOR13生態(tài)型的該基因編碼區(qū)第477位有一處同義突變,雖然堿基A突變?yōu)镚,但依舊編碼谷氨酸;其內(nèi)含子區(qū)域發(fā)生了 8處突變,包括7處點(diǎn)突變和1處缺失突變,在這處缺失突變中,來(lái)自生態(tài)型NOR13、NOR17、NOR21、NOR163、FIN11、FIN29、ICE5、ICE12和ICE16的基因在其第2個(gè)內(nèi)含子的第46～51個(gè)堿基處缺失了6個(gè)堿基(圖4)。由此可見(jiàn),基因FvH4_1g15230相較于FvH4_5g22710具有較大的遺傳變異性。

圖2 FvH4_1g15230基因編碼區(qū)在20種森林草莓生態(tài)型中突變位點(diǎn)序列比對(duì)Fig.2 Alignment of FvH4_1g15230 coding sequences around the mutation sites in 20 ecotypes of woodland strawberry圖中省略號(hào)代表無(wú)突變發(fā)生的核苷酸序列。The ellipsis means the nucleotide sequences without mutation.

圖3 FvH4_1g15230基因編碼的蛋白在20種森林草莓生態(tài)型中突變位點(diǎn)序列比對(duì)Fig.3 Alignment of the protein sequences encoded by FvH4_1g15230 gene around the mutation sites in 20 ecotypes of woodland strawberry圖中省略號(hào)代表無(wú)突變發(fā)生的氨基酸序列。The ellipsis means the amino acid sequences without mutation.

圖4 FvH4_1g15230基因內(nèi)含子區(qū)域在20種森林草莓生態(tài)型中的突變位點(diǎn)序列比對(duì)Fig.4 Alignment of FvH4_1g15230 intron regions around the mutation sites in 20 ecotypes of woodland strawberry圖中省略號(hào)代表無(wú)突變發(fā)生的核苷酸序列。The ellipsis means the nucleotide sequences without mutation.

2.3 2個(gè)NBS-LRR基因所編碼蛋白的結(jié)構(gòu)域分析

FvH4_5g22710基因編碼的蛋白中存在3種結(jié)構(gòu)域,分別為Rx_N(N-端)、NB-ARC結(jié)構(gòu)域和LRR基序;FvH4_1g15230基因編碼的蛋白中同樣存在3種結(jié)構(gòu)域,分別為RPW8(N-端)、NB-ARC結(jié)構(gòu)域和LRR基序(圖5)。兩者的N-端結(jié)構(gòu)域不同,前者屬于X-NBS-LRR類型基因,后者屬于RPW8-NBS-LRR類型基因。FvH4_5g22710基因中的2個(gè)突變位點(diǎn)均發(fā)生在Rx_N結(jié)構(gòu)域中;而在FvH4_1g15230基因編碼區(qū)的17個(gè)突變位點(diǎn)中,有7處點(diǎn)突變分布在3種結(jié)構(gòu)域中。

圖5 2個(gè)目的基因的蛋白質(zhì)結(jié)構(gòu)域組成及突變位點(diǎn)分布Fig.5 The protein domains and locations of mutation sites for the two genes 紅色箭頭為非同義突變,藍(lán)色箭頭為同義突變。The red arrows are non-synonymous mutation sites,and the blue arrows are synonymous mutation sites.

2.4 2個(gè)NBS-LRR基因系統(tǒng)進(jìn)化樹(shù)分析

由圖6可見(jiàn):FvH4_5g22710基因的系統(tǒng)進(jìn)化樹(shù)中,來(lái)自21種生態(tài)型的基因聚類為3組,其中的一組,其基因序列與FvH4的相比并無(wú)差異,包括采自芬蘭、意大利和挪威的生態(tài)型FIN5、FIN6、FIN11、FIN29、IT3、IT18、NOR13、NOR21和NOR122。其余2組的基因序列均與FvH4有一定程度的差異,第2組包括生態(tài)型ES5、ES9、ES18、UK1、UK9、UK10、IT24,大部分來(lái)自西班牙和英國(guó);第3組包括生態(tài)型ICE4、ICE16、NOR163、NOR150,分別采集自冰島和挪威。由此可見(jiàn),來(lái)源于同一國(guó)家或鄰近國(guó)家的生態(tài)型聚類在同一進(jìn)化支上,說(shuō)明它們的基因序列之間具有較小的核苷酸差異度。

FvH4_1g15230基因的系統(tǒng)進(jìn)化樹(shù)中也呈現(xiàn)出相似的聚類方式,來(lái)自20個(gè)生態(tài)型的基因分為4組,第1組由UK1、UK9、UK10、NOR122、NOR150、IT3、IT18、ES9、ES18等生態(tài)型的基因組成,大部分來(lái)自英國(guó)、挪威、意大利和西班牙;第2組僅包括來(lái)自芬蘭的生態(tài)型FIN11;第3組則大部分來(lái)自挪威的生態(tài)型NOR13、NOR17、NOR21;第4組則是多來(lái)自冰島的生態(tài)型ICE5、ICE12、ICE16。

圖6 來(lái)自森林草莓不同生態(tài)型的2個(gè)NBS-LRR基因系統(tǒng)發(fā)育樹(shù)Fig.6 Phylogenetic trees of two NBS-LRR genes in different ecotypes of woodland strawberry

2.5 選擇壓力檢驗(yàn)及核苷酸多樣性分析

FvH4_5g22710基因和FvH4_1g15230基因所受選擇壓力情況由非同義替換率(Ka)與同義替換率(Ks)的比值(Ka/Ks)和ω來(lái)反映,兩者的平均值均小于1(表4),共同說(shuō)明了2個(gè)基因均受純化選擇壓力的驅(qū)動(dòng)。此外,計(jì)算了兩者編碼的不同結(jié)構(gòu)域區(qū)段的Ka/Ks值,結(jié)果也展現(xiàn)了不同結(jié)構(gòu)域區(qū)段所受的選擇壓力情況(表5)。在FvH4_5g22710基因中,2處核苷酸突變位點(diǎn)均分布在Rx_N結(jié)構(gòu)域中,且平均Ka/Ks值為0.291 6,表明Rx_N結(jié)構(gòu)域區(qū)段仍受到純化選擇的作用。而基因FvH4_1g15230中,3種結(jié)構(gòu)域區(qū)段中均分布有突變位點(diǎn):RPW8結(jié)構(gòu)域區(qū)段中只發(fā)生非同義突變;NB-ARC和LRR結(jié)構(gòu)域中包含同義突變和非同義突變位點(diǎn),它們的平均Ka/Ks值分別為0.242 9和0.478 7,均小于1,說(shuō)明這2個(gè)結(jié)構(gòu)域區(qū)域受到了純化選擇。

為進(jìn)一步驗(yàn)證2個(gè)基因所受的選擇壓力,計(jì)算ω(dN/dS)和2Δln值。采用PAML程序的LR檢測(cè),驗(yàn)證基因編碼蛋白中受到正選擇作用的氨基酸位點(diǎn)(表4)。FvH4_5g22710的2Δln值為9.489 8,說(shuō)明其具有極顯著(2Δln>9.210,P<0.01)的正選擇作用位點(diǎn),尤其是其編碼蛋白序列的第53位的氨基酸位點(diǎn)(谷氨酸)受到正選擇的驅(qū)動(dòng)。

FvH4_5g22710的Pi值小于FvH4_1g15230(表4),說(shuō)明后者在不同生態(tài)型間的核苷酸差異度更大。在兩者不同的結(jié)構(gòu)域區(qū)段中,前者只在Rx_N結(jié)構(gòu)域中存在核苷酸差異度;而在后者中,LRR結(jié)構(gòu)域區(qū)段的Pi值最大即基因間核苷酸差異度最大,RPW8結(jié)構(gòu)域區(qū)段的Pi值最小即基因間核苷酸差異度最小(表5)。

表4 不同生態(tài)型中2個(gè)NBS-LRR基因的平均Ka/Ks值、Pi值及選擇壓力檢驗(yàn)Table 4 Average Ka/Ks value,Pi value and selection pressure test of two NBS-LRR genes in different ecotypes of woodland strawberry

表5 不同生態(tài)型中2個(gè)NBS-LRR基因編碼結(jié)構(gòu)域的平均Ka/Ks值及Pi值Table 5 Average Ka/Ks value and Pi value of domain regions encoded by two NBS-LRR genes in different ecotypes of woodland strawberry

3 討論

NBS-LRR基因是植物中存在最多的一類抗病基因,是植物細(xì)胞對(duì)抗病蟲(chóng)害侵染的重要防御組成成分。本研究對(duì)歐洲不同的森林草莓生態(tài)型進(jìn)行目標(biāo)NBS-LRR基因的遺傳變異分析。序列比對(duì)、核苷酸差異度分析結(jié)果表明FvH4_5g22710基因的遺傳變異程度小于FvH4_1g15230,推測(cè)這2個(gè)抗病基因可能擁有不同的進(jìn)化速率,以此應(yīng)對(duì)具有不同變異速率的病蟲(chóng)害[10,14]。

生態(tài)型是指來(lái)自同一物種的不同類群為了適應(yīng)不同生態(tài)環(huán)境,產(chǎn)生不同程度的遺傳變異,且能夠在類群中穩(wěn)定遺傳[2,15]。由此可見(jiàn),不同的生態(tài)型經(jīng)歷了自然選擇后,不同的基因積累了不同程度的遺傳變異。本研究中,系統(tǒng)進(jìn)化樹(shù)顯示,來(lái)自同一國(guó)家或是鄰近國(guó)家的生態(tài)型中的目標(biāo)基因,偏好聚類于同一個(gè)進(jìn)化支中,說(shuō)明相較其他基因而言,來(lái)自相同或相近地理位置的生態(tài)型,其基因序列之間具有更高的核苷酸相似度。類似的情況在不同生態(tài)型的擬南芥中也有發(fā)現(xiàn),例如孫彩玉等[16]對(duì)在意大利和瑞典不同地區(qū)取樣的擬南芥進(jìn)行春化作用相關(guān)基因(FLOWERINGLOCUSC,FLC)的序列分析發(fā)現(xiàn),采自意大利不同地區(qū)的擬南芥FLC基因序列完全相同,而來(lái)自瑞典不同地區(qū)的擬南芥FLC基因序列也完全相同。然而,來(lái)自意大利和瑞典的擬南芥在基因序列第27、461、501、638、738、884位的堿基卻不相同。

根據(jù)選擇壓力結(jié)果,目標(biāo)基因及其不同結(jié)構(gòu)域區(qū)段的平均Ka/Ks值均小于1,說(shuō)明它們?cè)谶M(jìn)化過(guò)程中均受到純化選擇的作用。Zhong等[10-11]發(fā)現(xiàn)在5個(gè)薔薇科物種和6個(gè)草莓物種中,大部分NBS-LRR基因受到純化選擇壓力的驅(qū)動(dòng),但仍有部分受到正選擇壓力作用;Jia等[17]發(fā)現(xiàn)蘋果(Malusdomestica)基因組的大部分NBS-encoding基因家族受到正選擇的壓力。豆科模式植物百脈根(Lotusjaponicus)以及白色幾內(nèi)亞山藥(Dioscorearotundata)基因組中的NBS-LRR基因基本受到純化選擇的驅(qū)動(dòng),僅有少部分基因受到正選擇的作用[18-19]。正如FvH4_5g22710基因的第53位的氨基酸位點(diǎn)(谷氨酸)受到正選擇的驅(qū)動(dòng)。植物抗病基因受到正選擇的驅(qū)動(dòng),可提高抗病基因核苷酸序列的多樣性,有效應(yīng)對(duì)病蟲(chóng)害的快速變異,提高植物的抗病性,為栽培草莓抗病新品種的培育提供理論依據(jù)和候選基因資源。

通過(guò)序列比對(duì)并結(jié)合結(jié)構(gòu)域分析發(fā)現(xiàn),一部分突變位點(diǎn)分布于目的基因編碼的蛋白結(jié)構(gòu)域中。其中,FvH4_5g22710基因中的2個(gè)突變位點(diǎn)均存在于Rx_N結(jié)構(gòu)域區(qū)段。Hao等[20]發(fā)現(xiàn)Rx_N是馬鈴薯Rx抗病蛋白N端結(jié)構(gòu)域,它與NBS和LRR共同介導(dǎo)對(duì)馬鈴薯病毒X的抗性。而FvH4_1g15230基因在其編碼的3種結(jié)構(gòu)域RPW8、NB-ARC、LRR中均發(fā)生了變異,Xiao等[21]發(fā)現(xiàn)擬南芥中的抗病基因RPW8.1和RPW8.2對(duì)白粉病病原體具有廣譜抗性;薔薇科中的RNL基因劃歸為2大類,分別具有較高和較低的差異度,呈現(xiàn)出快速和慢速2種進(jìn)化模式,用以抵御不同類型的病蟲(chóng)害[10]。LRR基序是一種富含亮氨酸的重復(fù)序列,廣泛存在于植物抗病基因編碼的蛋白中[22]。LRR結(jié)構(gòu)域與病原體效應(yīng)分子的特異性識(shí)別有關(guān),并且在抗病信號(hào)轉(zhuǎn)導(dǎo)中具有調(diào)節(jié)作用[23];由于LRR結(jié)構(gòu)域與致病效應(yīng)因子的相互識(shí)別和相互作用,LRR結(jié)構(gòu)域與N-端結(jié)構(gòu)域、NBS結(jié)構(gòu)域相比經(jīng)歷了更頻繁的正選擇壓力作用并具有更高的差異度[24-26]。本研究可為后續(xù)研究森林草莓NBS-LRR基因抗病功能以及草莓抗病分子育種提供重要的理論參考。

致謝:南京農(nóng)業(yè)大學(xué)生物信息中心高性能計(jì)算機(jī)平臺(tái)在實(shí)驗(yàn)技術(shù)方面給予了支持和幫助,謹(jǐn)致謝意。