孫中潔,王振英,潘麗娜
(天津師范大學(xué)a.生命科學(xué)學(xué)院,b.天津市動植物抗性重點(diǎn)實(shí)驗(yàn)室天津300387)
小麥白粉病是由禾布氏白粉菌(Blumeria graminis(DC).Speer f.sp.tritici)引起的一種世界性病害,它直接影響小麥的高產(chǎn)、穩(wěn)產(chǎn),是小麥的主要病害之一[1].目前防治小麥白粉病的工作主要集中在兩方面,即新抗病基因的發(fā)掘和新型抗病品系的選育[2-4],而對白粉菌脅迫應(yīng)答的生理機(jī)制研究較少,表觀遺傳修飾尤其是DNA甲基化在白粉病抗性應(yīng)答中的作用研究更是缺乏.
表觀遺傳是指在DNA序列不變的情況下,可遺傳的基因表達(dá)的改變[5-6].可逆的表觀遺傳修飾能有效避免不必要的基因重組,在調(diào)控基因表達(dá)的同時(shí)維持基因組的穩(wěn)定性,對于植物的生長、發(fā)育、器官分化及其對環(huán)境脅迫的適應(yīng)性具有關(guān)鍵作用[7-9].其中,DNA甲基化是一種最保守的表觀遺傳修飾,指在DNA甲基轉(zhuǎn)移酶(Dnmt)的作用下,以S-腺苷甲硫氨酸(SAM)為甲基供體,將甲基基團(tuán)轉(zhuǎn)移到胞嘧啶和鳥嘌呤(CpG)二核苷酸中胞嘧啶的5位碳原子上[5,10].脊椎動物中約有3%~8%基因組DNA的胞嘧啶處于甲基化狀態(tài),而在高等植物中約有5%~30%的胞嘧啶是甲基化的[11].DNA甲基化可形成一種封閉的染色質(zhì)結(jié)構(gòu),抑制基因轉(zhuǎn)錄.對擬南芥全基因組甲基化圖譜的分析表明,DNA甲基化在轉(zhuǎn)座子和基因編碼區(qū)的分布不同,轉(zhuǎn)座子及異染色質(zhì)區(qū)域的基因處于高甲基化狀態(tài),約1/3基因的編碼區(qū)是甲基化的,但它們的啟動子區(qū)很少發(fā)生甲基化修飾[12-14].
本課題組以抗白粉病代換系6A/6V及其與京411的近等基因系(NILs)為實(shí)驗(yàn)材料,通過甲基化特異性片段長度的多態(tài)性(MSAP)分析其DNA甲基化模式的變化情況,為進(jìn)一步闡明表觀遺傳尤其是DNA甲基化修飾在小麥白粉菌脅迫應(yīng)答中的作用奠定基礎(chǔ).
感病品系京411由中國農(nóng)業(yè)大學(xué)楊作民教授惠贈;抗白粉病代換系6A/6V由南京農(nóng)業(yè)大學(xué)提供;抗病近等基因系NILs是以京411為母本,與6A/6V雜交得F1代,在F1代中選取抗白粉病植株,再以京411為父本,回交7代,再自交1代后獲得.
H/M-adapter I接頭、H/M-adapter II接頭、EcoRI-adapter I接頭、EcoRI-adapter II接頭、H/M+0引物、EcoRI+A引物、8個(gè)H/M+3引物及9個(gè)E+3引物均在上海生工生物工程有限公司合成;限制性內(nèi)切酶HpaⅡ、MspⅠ及EcoR I購自NEB公司;DNA聚合酶、T4連接酶購自寶生物工程(大連)有限公司.
1.2.1 基因組DNA的提取
取同一生長時(shí)期的0.1~0.2 g小麥葉片,采用CTAB法提取基因組DNA,利用Nanodrop微量紫外分光光度計(jì)(Nanodrop科技有限責(zé)任公司)和瓊脂糖凝膠電泳檢測DNA提取的濃度及質(zhì)量.
1.2.2 甲基化敏感限制性內(nèi)切酶酶切與連接
分別采用EcoR I/Hpa II和EcoR I/Msp I酶切基因組DNA,酶切體系如下:EcoR I 5 U,Hpa II(Msp I)5 U,10×Buffer 2.5 μL,基因組 DNA 500 ng,加水至25 μL.上述溶液混勻,37℃水浴酶切2 h后,80℃溫浴20 min,使內(nèi)切酶失活.
接頭引物各3 μL,經(jīng)95℃、5 min變性后,溫度遞降至室溫.按如下體系進(jìn)行連接:10×T4連接酶Buffer 2 μL,T4連接酶 0.1 μL,E 接頭(5 μmol/L)1 μL,H(或 M)接頭(50 μmol/L)1 μL,酶切產(chǎn)物5 μL,加水至20 μL.20℃水浴連接過夜,然后65℃溫浴20 min,使連接酶失活.
1.2.3 預(yù)擴(kuò)增
反應(yīng)體系如下:Taq DNA聚合酶(5 U/μL)0.1 μL,dNTP(10 mmol/L)0.4 μL,10×Buffer 2 μL,連接產(chǎn)物 2 μL,E0 引物(50 ng/μL)0.6 μL,M0/H0引物(50 ng/μL)0.6 μL,加水至 20 μL.將上述溶液混勻,采用以下反應(yīng)程序進(jìn)行預(yù)擴(kuò)增:94℃預(yù)變性2 min,94℃變性30 s,56℃退火1 min,72℃延伸1 min,30個(gè)循環(huán)后,72℃延伸10 min.利用質(zhì)量分?jǐn)?shù)為2%的瓊脂糖凝膠電泳檢測,選取200~1 000 bp呈彌散狀分布的產(chǎn)物,根據(jù)電泳亮度稀釋后進(jìn)行選擇性擴(kuò)增.
1.2.4 選擇性擴(kuò)增
反應(yīng)體系如下:Taq DNA聚合酶(5 U/μL)0.1 μL,10 × buffer 2.5 μL,dNTPs(10 mmol/L)0.5 μL,預(yù)擴(kuò)增稀釋樣品 2 μL,上游引物 1 μL,下游引物1 μL,加水至25 μL.混合上述溶液,采用如下體系進(jìn)行選擇性擴(kuò)增:94℃預(yù)變性2 min,94℃變性30 s,65℃每個(gè)循環(huán)遞降0.7℃,退火30 s,72℃延伸1 min,13個(gè)循環(huán),隨后56℃退火30 s,72℃延伸1 min,23個(gè)循環(huán)后,72℃延伸10 min.擴(kuò)增后,樣品中加入7 μL loading buffer(質(zhì)量分?jǐn)?shù)為98%的去離子甲酰胺;10 mmol/L EDTA,pH8.0;質(zhì)量分?jǐn)?shù)為0.25%的溴酚藍(lán);質(zhì)量分?jǐn)?shù)為0.25%的二甲苯青).95℃變性5 min,立即置于冰水混合物中準(zhǔn)備進(jìn)行聚丙烯酰胺凝膠電泳.
1.2.5 聚丙烯酰胺電泳及銀染顯帶
采用質(zhì)量分?jǐn)?shù)為6%的聚丙烯酰胺凝膠,以1×TBE為電泳緩沖液,恒功率65 W,預(yù)電泳約30 min,加選擇性擴(kuò)增樣品7 μL,待電泳至第1條指示帶到達(dá)膠板的另一前沿停止,電泳90~120 min分離選擇性擴(kuò)增產(chǎn)物.
電泳后將凝膠放入1 L質(zhì)量分?jǐn)?shù)為10%的冰乙酸溶液中,輕搖至指示劑無色,約25 min,再用1 L水洗膠板2次,每次5 min.水洗后把膠板放入銀染液(1 g AgNO3,1.5 mL質(zhì)量分?jǐn)?shù)為37%的甲醛,1 L水),輕搖約30 min.用蒸餾水快速水洗膠板5 s后放入1 L預(yù)冷的顯影液中(30 g無水碳酸鈉,1.5 mL質(zhì)量分?jǐn)?shù)為37%的甲醛,200 μL 10 mg/mL的硫代硫酸鈉)中迅速搖動,直至條帶出現(xiàn).最后,用1 L質(zhì)量分?jǐn)?shù)為10%的冰乙酸溶液終止反應(yīng),約5 min后水洗,室溫下干燥.
為研究DNA甲基化在小麥白粉病抗性應(yīng)答中的作用,課題組采用MSAP法分析抗病NILs及其感病親本京411、抗病親本6A/6V甲基化水平的變化.MSAP方法是參照經(jīng)典的分子標(biāo)記——擴(kuò)增長度多態(tài)性法(AFLP)演變而來,因此,其基本原理與AFLP大致相同.不同的是,AFLP法是將基因組DNA用2種限制性內(nèi)切酶消化,而MSAP法是用一對對甲基化敏感的同裂酶HpaⅡ和MspⅠ代替AFLP法中的一個(gè)限制性內(nèi)切酶對基因組DNA進(jìn)行消化,之后產(chǎn)生各種長度不一的具有黏性末端的酶切片段,接上人工接頭之后,作為PCR擴(kuò)增模板,利用設(shè)計(jì)好的不同引物對其進(jìn)行擴(kuò)增,得到擴(kuò)增產(chǎn)物之后進(jìn)行聚丙烯酰胺凝膠電泳,從而產(chǎn)生擴(kuò)增片段長度不同的電泳圖譜.這對同裂酶HpaⅡ和MspⅠ識別相同的四核苷酸位點(diǎn)5′-CCGG-3′序列,但是對此位點(diǎn)不同甲基化模式的敏感程度不同.HpaⅡ可以切割DNA單鏈上此位點(diǎn)的外側(cè)或內(nèi)側(cè)任意一種甲基化模式,卻不可切割DNA雙鏈上此位點(diǎn)的外側(cè)或內(nèi)側(cè)甲基化模式;而MspⅠ可以切割此序列內(nèi)側(cè)胞嘧啶甲基化模式,無論此DNA是單鏈還是雙鏈,且不管此序列外側(cè)是否甲基化,只要存在內(nèi)側(cè)胞嘧啶甲基化,MspⅠ均可將其切割成大小不一的片段[15].
因此,基于這一原理,用MSAP方法分析某DNA序列,經(jīng)這對同裂酶酶切之后可以產(chǎn)生4種情況(實(shí)驗(yàn)中如果出現(xiàn)條帶記為“+”,如果沒有出現(xiàn)條帶記為“-”):①如果通過HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切之后均產(chǎn)生條帶,即(++型),表明該序列的5′-CCGG-3′位點(diǎn)未發(fā)生甲基化;②如果通過HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切之后均未產(chǎn)生條帶,即(--型),表明雙鏈DNA的5′-CCGG-3′位點(diǎn)發(fā)生外部甲基化;③如果通過HpaⅡ/EcoRⅠ酶切產(chǎn)生條帶,而MspⅠ/EcoRⅠ酶切后沒有產(chǎn)生條帶,即(+-型),表明該序列發(fā)生了單鏈外側(cè)5′-CCGG-3′位點(diǎn)胞嘧啶的甲基化,即半甲基化;④如果通過HpaⅡ/EcoRⅠ酶切未產(chǎn)生條帶,而MspⅠ/EcoRⅠ酶切后產(chǎn)生條帶,即(-+型),表明該序列發(fā)生了雙鏈內(nèi)側(cè)5′-CCGG-3′位點(diǎn)胞嘧啶的甲基化,即全甲基化.
如圖1所示,本實(shí)驗(yàn)利用72對引物組合在NILs、京411與6A/6V基因組DNA HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切產(chǎn)物中進(jìn)行選擇性擴(kuò)增共得到4 601條帶,其中NILs、京411與6A/6V中擴(kuò)增得到條帶數(shù)分別為4 529條、4 482條和4 518條,3品系共有條帶4 411條,NILs與京411共有條帶63條,NILs與6A/6V共有條帶40條.
圖1 近等基因系NILs和親本京411、6A/6V中MSAP擴(kuò)增條帶數(shù)Fig.1 Amount of bands amplified by MSAP in NILs and its parents Jing 411,6A/6V
進(jìn)一步分析其整體甲基化水平發(fā)現(xiàn),抗病近等基因系NILs整體甲基化水平(22.9%)略低于其感病輪回親本京411(23.4%),高于其父本6A/6V(21.2%),如表1所示,說明NILs的整體基因表達(dá)水平可能高于京411,低于其抗病親本6A/6V.為排除基因重組對甲基化條帶的影響,對3品系共有的4 411條帶進(jìn)行整體甲基化水平分析,發(fā)現(xiàn)NILs的甲基化水平(21.1%)仍略低于京411(22.6%),而高于6A/6V(19.8%),如表2所示,說明抗白粉病品系的整體DNA甲基化水平低于感病品系,即相對高水平的基因表達(dá)參與了白粉病的抗性應(yīng)答.
表1 近等基因系NILs和親本京411,6A/6V中MSAP擴(kuò)增的條帶及甲基化率統(tǒng)計(jì)Tab.1 Amount of bands amplified by MSAP and the methylation rate among NILs and its parents Jing411,6A/6V
表2 近等基因系NILs和親本京411,6A/6V中MSAP擴(kuò)增的共有條帶及甲基化率統(tǒng)計(jì)Tab.2 Shared amount of bands amplified by MSAP and the methylation rate among NILs and its parents Jing411,6A/6V
MSAP分析既可在整體上估測全基因組的甲基化水平,也可用于特定位點(diǎn)不同品系基因組的甲基化模式.如圖2所示,相對親本京411,NILs大量位點(diǎn)的甲基化模式發(fā)生變異,超甲基化和去甲基化均有發(fā)生.根據(jù)HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切后擴(kuò)增泳帶的差異,可將MSAP擴(kuò)增位點(diǎn)分為A、B、C、D 4大類、12個(gè)亞類,如表3和圖2.其中A類NILs與京411的甲基化模式完全相同,即NILs相對京411無甲基化模式的變化.A類又可分3個(gè)亞類:A1代表HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切后均有條帶;A2為HpaⅡ/EcoRⅠ酶切有條帶,而MspⅠ/EcoRⅠ酶切后無條帶擴(kuò)增;A3亞類與A2相反,HpaⅡ/EcoRⅠ酶切無條帶擴(kuò)增,而MspⅠ/EcoRⅠ酶切后有條帶擴(kuò)增.B類NILs相對京411的甲基化水平降低,如B1、B2亞類,NILs基因組經(jīng)HpaⅡ/EcoRⅠ酶切后有條帶擴(kuò)增,而京411無條帶,說明該位點(diǎn)京411的DNA序列發(fā)生了外側(cè)的完全甲基化修飾;B3亞類中,NILs基因組經(jīng)MspⅠ/EcoRⅠ酶切后有條帶擴(kuò)增,而京411無條帶,說明該位點(diǎn)京411的DNA序列發(fā)生了內(nèi)側(cè)的甲基化修飾;B4亞類中,NILs經(jīng)2種酶切方式均有條帶,而京411無條帶,說明該位點(diǎn)發(fā)生了甲基化修飾,或者是基因重組位點(diǎn).C類為NILs相對京411甲基化水平提高的位點(diǎn),根據(jù)酶切擴(kuò)增條帶不同可細(xì)分為4個(gè)亞類:C1、C4為京411基因組經(jīng)MspⅠ/EcoRⅠ酶切后有條帶擴(kuò)增,而NILs無條帶;C2、C3亞類為京411基因組經(jīng)HpaⅡ/EcoRⅠ酶切后有條帶擴(kuò)增,NILs無條帶.D類為甲基化模式互換條帶,NILs經(jīng)HpaⅡ/EcoRⅠ酶切后有條帶擴(kuò)增,而京411基因組經(jīng)MspⅠ/EcoRⅠ酶切后有條帶擴(kuò)增,無法判斷其甲基化模式提高還是降低.如表3所示,A類條帶占絕大多數(shù)(96.54%),說明NILs相對京411具有極其相似的遺傳背景,基因表達(dá)差異不大;B類條帶比例(2.89%)明顯高于C類(0.55%),即NILs的甲基化水平明顯低于京411,通常高甲基化抑制基因表達(dá),說明NILs的整體基因表達(dá)水平高于京411.
圖2 部分MSAP擴(kuò)增結(jié)果Fig.2 Part of MSAP result
表3 近等基因系NILs與親本京411的甲基化模式Tab.3 Model of methylation between NILs and its parent Jing411
根據(jù)HpaⅡ/EcoRⅠ和MspⅠ/EcoRⅠ酶切后擴(kuò)增泳帶的差異,也可將NILs相對親本6A/6V的MSAP擴(kuò)增位點(diǎn)分為A、B、C、D 4大類、15個(gè)亞類,如表4所示.A類占92.52%,是NILs相對親本6A/6V甲基化模式不變的條帶;B類為NILs相對親本6A/6V甲基化水平降低的條帶;C類為NILs相對親本6A/6V甲基化水平提高的條帶;D類為NILs相對親本6A/6V發(fā)生了甲基化模式互換的條帶.如表4所示,C類(4.20%)條帶明顯高于B類(2.96%),說明NILs相對親本6A/6V甲基化水平明顯提高.
綜上,抗病NILs相對感病親本京411的甲基化水平有所降低,而相對抗病親本6A/6V甲基化水平明顯提高,說明NILs的整體基因表達(dá)水平高于京411,低于抗病親本6A/6V,即DNA甲基化水平低.基因表達(dá)相對高的品系NILs和6A/6V對白粉病有較強(qiáng)抗性,而整體基因表達(dá)相對低的品系京411對白粉病敏感,由此推測,DNA甲基化可能參與小麥白粉病的抗性應(yīng)答.
表4 近等基因系NILs與親本6A/6V的甲基化模式Tab.4 Model of methylation between NILs and its parent 6A/6V
植物在生長發(fā)育過程中,受到環(huán)境因子的多種脅迫刺激,如不適的溫度、壓力、水肥條件等物理性刺激,以及病蟲害等生物脅迫.為了生存與繁衍,生物體必須發(fā)展自身的耐受和適應(yīng)性能,進(jìn)而推動物種的進(jìn)化[16-17].脅迫誘導(dǎo)下,植物體需要發(fā)生多種代謝反應(yīng)變化以應(yīng)對和適應(yīng)脅迫,其中一種典型的反應(yīng)就是基因組變異,但通?;蛲蛔儭⒅嘏攀遣豢赡娴?,會導(dǎo)致物種的不穩(wěn)定性.而表觀遺傳修飾的改變,可在基因序列不變的情況下,引起基因表達(dá)的變化[5,18],并隨細(xì)胞分裂遺傳給子代,如此既賦予生物體耐受脅迫的能力,又不會造成基因的突變、重組,維持了物種的穩(wěn)定[3-5].小麥白粉病是一種真菌性病害,屬于典型的生物脅迫,但目前對表觀遺傳如DNA甲基化修飾在小麥白粉病脅迫應(yīng)答中的作用研究較少.
本實(shí)驗(yàn)應(yīng)用MSAP法分析比較抗病近等基因系NILs與其抗病父本6A/6V、感病輪回親本京411的甲基化差異.由于近等基因系在制備過程中經(jīng)歷了與感病親本京411的多次回交,NILs的遺傳背景與京411極其相似,而且抗病父本代換系6A/6V的V染色體(來自簇毛麥)是獨(dú)立遺傳[19-20],因此,可認(rèn)為NILs的基因組DNA除了6V染色體與代換系6A/6V相同,其他染色體與京411相同.NILs與京411對白粉病的抗性不同可能源于6V染色體上的基因,也可能在6V染色體滲入過程中引起整體基因表達(dá)的改變.相對于京411,NILs發(fā)生DNA甲基化水平降低的條帶明顯高于甲基化水平提高的條帶,且僅B4亞類是NILs獨(dú)有的擴(kuò)增條帶,可能來自6V染色體,是抗病基因.B1~B3亞類在NILs與京411中均有條帶擴(kuò)增,只是相對來講,NILs擴(kuò)增的條帶更多,說明外源染色體6V的滲入引起NILs整體基因組DNA甲基化水平降低,其原因可能是6V染色體上.外源抗病基因誘導(dǎo)下游白粉病脅迫應(yīng)答通路相關(guān)基因的DNA甲基化水平降低,基因表達(dá)提高,從而對小麥白粉菌侵染產(chǎn)生抗性.
綜上所述,NILs與其親本京411對白粉病抗性不同很可能與它們DNA甲基化水平的差異有關(guān),即表觀遺傳尤其是DNA甲基化參與了小麥白粉菌的脅迫應(yīng)答.
[1]LIU J,TAO W,DUAN X,et al.Molecular marker assisted identification of Pm genes involved in the powdery mildew resistant wheat cultivars[J].Acta Phytopathol Sin,2000,30(2):133—139.
[2]LIU J,LIU D.Progress of the study on wheat powdery mildew resistant genes[J].Acta Phytopathol Sin,2000,30(4):289—295.
[3]YU L,NIU J,MA Z,et al.Cloning,characterization and chromosome localization of two powdery mildew resistance-related gene sequences from wheat[J].Acta Bot Sin,2002,44(12):1438—1444.
[4]YAHIAOZU N,SRICHUMPA P,DUDLER R,et al.Genome analysis at different ploidy level allows cloning of the powdery mildew resistance gene Pm3b from hexaploid wheat[J].Plant J,2004,37(4):528—538.
[5]BIRD A.Perceptions of epigenetics[J].Nature,2007,447(7143):396—398.
[6]FAZZARIMJ,GREALLYJM.Epigenomics:BeyondCpGislands[J].Nat Rev Genet,2004,5(6):446—455.
[7]BENDERJ.Plantepigenetics[J].CurrBiol,2002,12(12):412—414.
[8]PARISOD C,SALMON A,ZERJAL T,et al.Rapid structural and epigenetic reorganization near transposable elements in hybrid and allopolyploid genomes in Spartina[J].New Phytol,2009,184(4):1003—1015.
[9]RAPP R A,WENDEL J F.Epigenetics and plant evolution[J].New Phytol,2005,168(1):81—91.
[10]SANTIDV,GARRETTCE,BARRPJ.Onthemechanismofinhibition of DNA-cytosine methyltrans-ferases by cytosine analogs[J].Cell,1983,33(1):9—10.
[11]CHEN T,LI E.Structure and function of eukaryotic DNA methyltransferases[J].Curr Top Dev Biol,2004,60:85—89.
[12]ZILBERMAN D,GEHRING M,TRAN R K,et al.Genomewide analysis of Arabidopsis DNA methylation uncovers an interdependence between methylation and transcription[J].Nat Genet,2007,39(1):61—69.
[13]ZHANG X,YAZAKI J,SUNDARESAN A,et al.Genomewide high-resolution mapping and functional analysis of DNA methylation in Arabidopsis[J].Cell,2006,126(6):1189—1201.
[14]GEHRING M,HENIKOFF S.DNA methylation dynamics in plant genomes[J].Biochimicaet Biophysica Acta-Gene Structure and Expression,2007,1769(5/6):276—286.
[15]MCCLELLAND M,NELSON M,RASCHKE E.Effect of site-specific modification on restriction endonucleases and DNA modificationmethyltransferases[J].NucleicAcidsRes,1994,22(17):3640—3659.
[16]ARNHOLDT-SCHMITT B.Stress-induced cell reprogramming:A role for global genome regulation[J].Plant Physiol,2004,136(1):2579—2586.
[17]DOROSZUK A,WOJEWODZIC M,KAMMENGA J.Rapid adaptive divergence of life-history traits in response to abiotic stress within a natural population of a parthenogenetic nematode[J].Proc Biol Sci,2006,273(1601):2611—2618.
[18]CHINNUSAMY V,ZHU J.Epigenetic regulation of stress responses in plants[J].Curr Opin Plant Biol,2009,12(2):133—139.
[19]段爽,王煉,甘富,等.具簇毛麥6VS的小麥抗白粉病近等基因系的培育與蛋白質(zhì)組檢測[J].中國細(xì)胞生物學(xué)學(xué)報(bào),2012,34(6):580—589.
[20]CHEN P,QI L,ZHOU B,et al.Development and molecular cytogenetic analysis of wheat-Haynaldia villosa 6VS/6AI.translocation lines specifying resistance to powdery mildew[J].Theor Appl Genet,1995,91(6/7):1125—1128.