D2型細胞質雄性不育小麥絨氈層細胞程序化死亡與活性氧代謝

劉子涵，石曉藝，閆鵬嬌，段陽，耿興俠，葉佳麗，李莎，楊雪桐，張高明，賈雨林，張玲麗，宋喜悅

?

D2型細胞質雄性不育小麥絨氈層細胞程序化死亡與活性氧代謝

劉子涵，石曉藝，閆鵬嬌，段陽，耿興俠，葉佳麗，李莎，楊雪桐，張高明，賈雨林，張玲麗，宋喜悅

（西北農林科技大學農學院，陜西楊凌 712100）

D2型細胞質是細胞質雄性不育小麥（cytoplasmic male sterility，CMS）的重要胞質來源，深入研究其敗育機理對小麥雜種優(yōu)勢利用具有重要價值。以3種同核異質D2型細胞質雄性不育小麥Va87B1-706A、C687B1-706A、Ju87B1-706A及其相應保持系706B為試材，通過體式顯微鏡和DAPI染色分別在四分體時期、單核早期、單核晚期、二核期和三核期觀察小麥花藥的外形和小孢子發(fā)育的形態(tài)；利用掃描電鏡和KI-I2染色分析三核期花藥外皮層、內皮層、烏氏體和花粉粒的表型及敗育特點；采用石蠟切片和DNA ladder的技術觀察不同發(fā)育時期D2型不育系和保持系小麥花藥絨氈層細胞形態(tài)變化特征以及絨氈層細胞程序化死亡（programmed cell death，PCD）的過程；測定花藥發(fā)育過程中活性氧（reactive oxygen species，ROS）的含量，超氧化物歧化酶（superoxide dismutase，SOD）、過氧化物酶（peroxidase，POD）、過氧化氫酶（hydrogen peroxidase，CAT）、抗壞血酸過氧化物酶（ascorbate peroxidase，APX）和谷胱甘肽過氧化物酶（glutathione peroxidase，GPX）的活性以及非酶物質還原型抗壞血酸（ascorbic acid，AsA）和谷胱甘肽（glutathione，GSH）的含量變化；同時利用RT-PCR技術對部分抗氧化酶相關基因的表達模式進行分析。與保持系706B比較，3種D2型細胞質雄性不育小麥在三核期均表現出外皮層表皮褶皺、烏氏體排列稀疏散亂以及小孢子皺縮、表皮粗糙和萌發(fā)孔閉合的不育特征，敗育類型為典敗和染敗；自單核晚期開始3種不育小麥的小孢子發(fā)生紊亂；石蠟切片觀察絨氈層細胞的變化，發(fā)現不育系均比保持系延遲一個時期在二核期發(fā)生解體，并于三核期完全降解；進一步檢測細胞凋亡DNA小片段發(fā)現3種不育系花藥絨氈層的DNA損傷均起始于單核晚期，絨氈層PCD比保持系延遲一個時期啟動，并且優(yōu)先細胞學表型檢測到凋亡；生理指標測定和定量分析的結果說明絨氈層細胞的延遲凋亡和小孢子的結構異常與不育系花藥內活性氧的過度積累、抗氧化酶活性的異常變化、非酶促抗氧化物質的嚴重下調以及部分抗氧化酶相關基因在兩系花藥中表達水平的嚴重差異密切相關。D2型細胞質雄性不育小麥敗育的關鍵時期是單核晚期，絨氈層PCD的延遲啟動誘導了活性氧的過量積累，進一步破壞了抗氧化防御系統(tǒng)的平衡，進而造成絨氈層細胞形態(tài)異常，最終導致D2型細胞質雄性不育小麥的敗育。

小麥；D2型細胞質雄性不育；絨氈層細胞程序化死亡；活性氧代謝；基因表達

0 引言

【研究意義】作物雜種優(yōu)勢是自然界普遍存在的一種生物學現象，實踐證明雜交種具有產量高，品質優(yōu)以及抗性強等優(yōu)勢[1]，許多植物雜交種的大面積推廣已取得一定的經濟效益[2]。因此，深入研究與開發(fā)利用作物雜種優(yōu)勢在提高小麥的品質和單位面積產量等方面具有重要作用，而雄性不育是實現雜種優(yōu)勢利用的一條重要途徑，對雄性不育機理的研究一直引起國內外科學家們的廣泛關注[3]。【前人研究進展】近50多年來，雜種小麥因為外源細胞質的引用而得以廣泛應用。1951年，Kihara創(chuàng)制的尾狀山羊草（）異質不育系開啟了小麥異質雄性不育的研究，之后相繼獲得130多種質核雜交小麥。其中，D2型細胞質具有一定應用潛力[4]。D2型細胞質對小麥主要農藝性狀無不利效應，且能引起許多小麥品種（系）雄性不育[5-6]。同時發(fā)現在D2型細胞質中，特別是牡山羊草（）細胞質對雜交小麥的開發(fā)利用價值最大，并且具有易恢復、增加普通小麥的生長勢、改良小麥品質以及提高小麥對白粉病的抗性等優(yōu)點，在雜交小麥育種和生產上具有廣闊的應用前景和市場價值[7-8]。盡管科學研究者們證明D2型細胞質對普通小麥的農藝性狀沒有不利效應[9-10]，但對這種細胞質雄性不育系及其敗育機理并不是很明確，且國內外研究甚少。細胞程序化死亡（programmed cell death，PCD）在多種生物的生長過程中占有重要地位，是一種由基因控制有序的選擇性細胞凋亡的過程[11]。在植物發(fā)育過程中，一些組織和器官的發(fā)育都會發(fā)生細胞程序化死亡過程，如許多特定組織細胞的分化、生殖器官的發(fā)育、葉片的自然凋亡、根冠細胞的消失以及糊粉層的瓦解等[12-14]。目前，許多學者也證實了花藥絨氈層降解的基本特征與細胞程序化死亡的過程相符[15-17]，作為花藥壁最內層的絨氈層直接和小孢子接觸，適時降解的絨氈層能為小孢子正常生長發(fā)育準確地提供脂肪、蛋白質和多糖等營養(yǎng)物質，在小孢子發(fā)育過程中起著極其重要的作用。絨氈層細胞的PCD過程提前（或延遲）均會導致植物（如水稻、玉米、向日葵、小麥等）花粉敗育[18-21]。目前，人們普遍認為，細胞程序化死亡與活性氧（reactive oxygen species，ROS）的代謝之間具有緊密的聯(lián)系，對于PCD與ROS的研究大多集中在生物脅迫和非生物脅迫對微生物、昆蟲、動物以及植物葉片的PCD影響[22-25]，而對于細胞質雄性不育小麥絨氈層的自然凋亡過程與ROS關系的研究較少。在植物體中，細胞中所有的活動所需要的能量都離不開線粒體中的呼吸作用。在正常情況下，植物線粒體產生的活性氧能夠抵御外來傷害，但一旦發(fā)生代謝紊亂，過量的ROS就會與細胞內的多種物質相互作用，導致細胞代謝的平衡。植物在面臨傷害的過程中形成一套可有效阻止或延遲氧化傷害的防御機制，即酶促抗氧化系統(tǒng)和非酶促抗氧化系統(tǒng)。其中，酶促抗氧化系統(tǒng)中的SOD酶能將超氧自由基歧化為H2O2和O2，POD協(xié)同CAT能同時清除生物體內過量的過氧化氫[26]。APX和GPX在過量H2O2的刺激下能直接參與植物的氧化反應，同時GSH被氧化成GSSG[27]。其中，馬翎健等[28]研究細胞質雄性不育小麥酶活性，發(fā)現小麥ROS含量升高、SOD和CAT活性降低，最終導致膜損傷，從而造成細胞質雄性不育小麥花粉敗育。非酶促抗氧化系統(tǒng)中，AsA-GSH循環(huán)系統(tǒng)也可以消除過多的H2O2進而調節(jié)細胞內的信號轉導途徑使植物向有利的方向發(fā)育[29]。研究表明不同的細胞質雄性不育小麥AsA和GSH代謝的異常是導致花粉敗育的根本原因[30]。目前，研究表明ROS與植物細胞的程序化死亡有關，并發(fā)現在PCD過程中ROS起到2種作用，一種是ROS作為信號分子啟動了細胞調亡的反應，一種是直接導致了細胞的PCD過程[12]?！颈狙芯壳腥朦c】小麥雄性不育的形成機理十分復雜，近年來，探究細胞質雄性不育敗育的原因已成為研究的熱點問題。而且關于細胞質雄性不育的研究一般采用核遺傳背景不同或者單一的材料進行比較，具有很大的局限性?！緮M解決的關鍵問題】本研究選用3種同核異質D2型細胞質雄性不育小麥，在消除核背景遺傳差異的前提下，從細胞學、生理生化和分子生物學等方面，利用掃描電鏡、DAPI染色、石蠟切片、DNA ladder、生理指標測定和qRT-PCR的方法，客觀詳細地分析了3種細胞質雄性不育系絨氈層PCD、ROS及抗氧化系統(tǒng)三者之間的關系，進一步揭示雄性敗育原因，為細胞質不育類型的劃分以及探索敗育機理開拓了新的方向。

1 材料與方法

1.1 植物材料

試驗材料Va87B1-706A、C687B1-706A和Ju87B1-706A，其細胞質分別源于瓦維洛夫山羊草（）、粗厚山羊草（）和牡山羊草（）。3種D2型細胞質雄性不育系以及同型保持系706B的細胞質源于普通小麥細胞質（）。所有不育系均與保持系706B于楊凌回交20代以上，不育性穩(wěn)定。所有試驗材料均在2014年10月和2015年10月播種于西北農林科技大學實驗田。2015年4月，對所有不育系隨機選10株進行自交套袋，并調查自交結實率（均為0），進一步說明不育系的穩(wěn)定性（姚盟2015）。每種材料種植2行，行長1 m，行間距0.25 m，株距6.7 cm，常規(guī)大田管理。

1.2 花藥及小孢子的表型觀察

于4月份選取3種D2型細胞質雄性不育系以及同型保持系穗中部小穗基部的花藥，用醋酸洋紅進行壓片鏡檢，確定小麥的不同發(fā)育時期（stage1，四分體時期、stage2，單核早期、stage3，單核晚期、stage4，二核期和stage5，三核期）。通過已經鑒定的小麥發(fā)育時期，記錄相應的植株形態(tài)，確保材料時期的準確性。利用體式顯微鏡（Motic K400）對不同時期的小花和花藥進行觀察和拍照。按照王書平[31]的掃描電鏡方法觀察三核期的花藥，之后用導電膠將花藥外皮、內皮和小孢子固定在銅臺上進行噴金處理，用掃描電鏡（M-6360LV）觀察其表型結構并拍照。利用FAA固定液固定不同時期的花藥，保存于4℃，用于DAPI染色和I2-K染色制片。

1.3 花藥絨氈層觀察

取出FAA固定液中各時期的花藥，參照王書平[31]的方法制作石蠟切片，切片厚度為10 μm，所有樣品制片在顯微鏡（Nikon ECLIPSE E600）下觀察并用NIS-Elements軟件拍照。

1.4 DNA的提取純化與DNA ladder的分析

取5個時期的花藥各0.25 g，參照王書平[31]的方法提取DNA并純化，將純化的DNA用核酸蛋白質檢測儀測定濃度，取10 μg DNA于2.0%的瓊脂糖凝膠進行電泳，利用凝膠成像系統(tǒng)（JY04S-3C）拍照。

1.5 生理指標的測定

在4℃環(huán)境下，選取3種D2型細胞質雄性不育系以及同型保持系5個時期的花藥進行主要生理指標的測定，重復3次，取其平均值。

1.5.1 活性氧（ROS）含量的測定 參照王書平[31]的方法測定超氧陰離子（O2-）和過氧化氫（H2O2）含量，參考Noreen等[32]方法測定丙二醛（MDA）含量，并稍有改進。

1.5.2 抗氧化酶活性的變化 參照Li等[33]方法并稍作修改測定超氧化物歧化酶（superoxide dismutase，SOD）和過氧化物酶（peroxidase，POD）活性。參照楊冬業(yè)等[34]的方法測定過氧化氫酶（hydrogen peroxidase，CAT）活性。參照張悅等[27]方法測定抗壞血酸過氧化物酶（ascorbate peroxidase，APX）以及谷胱甘肽過氧化物酶（glutathione peroxidase，GPX）活性，并加以改進。

1.5.3 非酶促抗氧化物質含量的變化 按照黃志明等[35]方法測定還原型抗壞血酸（ascorbic acid，AsA）和谷胱甘肽（glutathione，GSH），并稍作改動。

1.6 抗氧化酶相關基因表達分析

1.6.1 RNA的提取和cDNA第一鏈的合成 采用TRIzol法提取3種D2型細胞質雄性不育系及同型保持系花藥5個時期的總RNA。利用HiScriptTMQ Select RT SuperMix for qPCR（Vazyme 南京）進行cDNA反轉錄，方法參考說明書，cDNA產物用于RT-PCR反應。

1.6.2 SOD、CAT和APX酶基因的半定量分析及熒光定量PCR分析 所用引物序列（電子附表1）由上海生物工程股份有限公司合成。RT-PCR參照石曉藝[36]方法進行操作，利用凝膠成像系統(tǒng)（JY04S- 3C）拍照。

用AceQ?qPCR SYBR?Green Master Mix 試劑盒（Vazyme 南京），按照試劑盒說明書在ABI PRISM 7000上進行操作。每個模板設3次重復。目的基因相對表達量用2-ΔΔCt的方法進行計算。

1.7 數據分析

采用Microsoft Office Excel 2010和SPSS 19.0進行數據統(tǒng)計與分析，利用SSR法對數據進行方差分析。

2 結果

2.1 D2型細胞質雄性不育系敗育的表型特征

從不育系和保持系5個發(fā)育時期小花的結果可以看出，不育系的雌蕊發(fā)育均正常。然而，從單核晚期開始不育系的花藥比706B小且顏色偏淡（圖1）。尤其在三核期時，所有不育系花藥端部不開裂，底部不分叉，捻開后沒有花粉或者只可見少量花粉散出（圖1和圖2）。而保持系706B花藥飽滿且顏色亮黃，兩端分叉明顯，花藥中有大量花粉散出（圖1和圖2）。

A—E：Va87B1-706A；F—J：C687B1-706A；K—O：Ju87B1-706A；P—T：706B。A、F、K和P：四分體時期；B、G、L和Q：單核早期；C、H、M和R：單核晚期；D、I、N和S：二核期；E、J、O和T：三核期。不育系花藥從四分體時期到單核早期正常發(fā)育，但是從單核晚期開始不育系的花藥要比706B偏小且顏色偏淡。尤其在三核期時，所有不育系花藥干癟，不散粉。標尺為 0.5 mm

為進一步觀察不育系三核期花藥和小孢子的敗育特征，利用掃描電鏡觀察不育系和保持系三核期花藥的外皮層、內皮層和小孢子的形態(tài)。結果顯示，保持系的外皮層表皮光滑平整，細胞排列整齊緊密有序，而不育系的外皮層表皮褶皺粗糙，細胞排列散亂且沒有規(guī)則，細胞間相互疊加呈現出扭曲變形的不育特征（圖2）。內皮層表面是由大量烏氏體覆蓋而成，從圖中可以看出保持系內皮層表面烏氏體排列致密，形態(tài)規(guī)則且個體飽滿，而不育系內皮層表面烏氏體呈現相互黏連、排列稀疏、大小不一的不育特征（圖2）。通過掃描電鏡觀察小孢子的特征，發(fā)現保持系的小孢子形狀圓潤表面光滑，萌發(fā)孔呈現開放狀態(tài)。而不育系的小孢子表面粗糙且容易黏連微粒，形狀不規(guī)則，萌發(fā)孔閉合或皺縮（圖2）。

I2-KI染色鑒定結果顯示，3種不育系大部分花粉粒形態(tài)不規(guī)則有皺縮現象且不染色，少部分呈圓形且染色淺，因而敗育類型為典敗和染敗。保持系706B花粉粒經染色后，絕大多數花粉粒被染成棕黑色且染色充分（圖2）。

以上結果表明，D2型細胞質雄性不育系的質核互作效應嚴重影響了花藥和小孢子的正常發(fā)育，最終導致完全雄性不育。

A1—A9：Va87B1-706A；B1—B9：C687B1-706A；C1—C9：Ju87B1-706A；D1—D9：706B；A1、A6、B1、B6、C1、C6、D1和D6：花藥；A2、A3、B2、B3、C2、C3、D2和D3：外皮層細胞；A4、A5、B4、B5、C4、C5、D4和D5：內皮層細胞；Uby：烏氏體；A7-D7：碘染花粉粒；A8、A9、B8、B9、C8、C9、D8和D9：小孢子；Ap：萌發(fā)孔。保持系的外皮層表皮光滑平整，內皮層表面烏氏體排列致密形態(tài)規(guī)則且個體飽滿，小孢子形狀圓潤表面光滑；不育系外皮層表皮褶皺粗糙，細胞排列散亂，內皮層表面烏氏體現相互黏連、排列稀疏，小孢子表面粗糙且容易黏連微粒，且敗育類型為典敗和染敗。A1、A6、B1、B6、C1、C6、D1和D6的標尺為1mm；A2、B2、C2和D2的標尺為100 μm；A7—A9、B7—B9、C7—C9和D7—D9的標尺為50 μm；A3—A5、B3—B5、C—C5和D3—D5的標尺為10 μm

A1-A9: Va87B1-706A; B1-B9: C687B1-706A; C1-C9: Ju87B1-706A; D1-D9: 706B. A1, A6, B1, B6, C1, C6, D1, and D6: anther; A2, A3, B2, B3, C2, C3, D2, and D3: outer epidermal cells; A4, A5, B4, B5, C4, C5, D4, and D5: inner epidermal cells; Uby: indicates the Ubisch bodies; A7, B7, C7 and D7: microspore of I2-KI staining. A8, A9, B8, B9, C8, C9, D8 and D9: microspores; Ap: indicate the germination aperture. The outer epidermal cells of fertile anthers appeared to be slippy and smooth, and the inner epidermal cells accumulated a lot of tight Ubisch bodies. The epidermis of fertile microspores was rounded and plump with a gymnotremoid germinal aperture. The outer epidermal cells of sterile anthers appeared to be adhesive and loose, and the inner epidermal cells accumulated less and more unconsolidated Ubisch bodies. The epidermis of the CMS microspores were crenate and extremely scabrous, with a malformed germinal aperture. Scale bars represent 1mm in A1, A6, B1, B6, C1, C6, D1, and D6; 100 μm in A2, B2, C2 and D2 ; 50 μm in A7 to A9, B7 to B9, C7 to C9, and D7 to D9; 10 μm in A3 to A5, B3 to B5, C3 to C5, and D3 to D5

圖2 三核期不育系和保持系花藥及小孢子外部形態(tài)及育性的比較

Fig. 2 Comparison of scanning electron micrograph observations, I2-KI staining, and anther morphology in Va87B1-706A, C687B1-706A, Ju87B1-706A and fertile wheat plants at the trinucleate stage

2.2 D2型細胞質雄性不育系敗育的細胞學特征

為進一步探究D2型細胞質雄性不育系小孢子的發(fā)育，筆者分別對不同時期的花藥細胞進行了DAPI染色觀察（圖3）。保持系706B和不育系的花粉母細胞通過減數分裂形成四分體，單核早期小孢子細胞核于中央。與保持系單核晚期圓形的小孢子不同，不育系小孢子呈不規(guī)則形狀且表面皺縮，細胞核位于萌發(fā)孔對側；二核期保持系小孢子出現圓形的營養(yǎng)核和精核；在三核期，圓形的精核形成了2個梭形的精核，營養(yǎng)核顏色變淺。不育系Va87B1-706A、C687B1-706A和Ju87B1-706A，二核期花粉粒能正常形成精核和營養(yǎng)核，但有些細胞形狀不規(guī)則且營養(yǎng)核模糊。而三核期大部分花粉粒細胞精核呈現圓形退化變小，不能形成2個梭形的精核，且有部分細胞呈空胞化，花粉敗育。

A—E：Va87B1-706A；F—J：C687B1-706A；K—O：Ju87B1-706A；P—T：706B；A、F、K和P：四分體時期；B、G、L和Q：單核早期；C、H、M和R：單核晚期；D、I、N和S：二核期；E、J、O和T：三核期。不同于保持系小孢子的正常發(fā)育，不育系單核晚期細胞表面皺縮且呈不規(guī)則形狀；二核期花粉粒營養(yǎng)核不清晰且細胞形狀不規(guī)則；而三核期大部分花粉粒精核呈現圓形退化變小，且有部分細胞發(fā)生空胞化。標尺為 50 mm

2.3 花藥絨氈層細胞程序化死亡的觀察

為了觀察D2型細胞質雄性不育小麥花藥發(fā)育過程中絨氈層細胞程序化死亡啟動時期的異常變化，并進一步確定敗育時期和原因，利用甲苯胺藍染色法觀察不育系和保持系不同發(fā)育時期花藥絨氈層的變化（圖4）。在四分體時期，不育系和保持系花藥絨氈層體積明顯偏大，細胞質濃厚染色深，花藥壁的中層結構完整，呈帶狀分布。單核早期，不育系和保持系花藥絨氈層細胞質染色依然很深，中層開始變窄，且絨氈層沒有明顯區(qū)別。到了單核晚期，不育系和保持系的中層細胞完全退化消失，而保持系的絨氈層開始降解，即啟動了細胞程序化死亡的過程，此時期表現為絨氈層邊界模糊，并且能夠明顯地觀察到未退化的細胞核碎片及較多的細胞質殘留。而不育系在此時期并沒有觀察到絨氈層的退化，且排列緊密，結構完整，同時由于絨氈層未及時發(fā)生降解，不能為小孢子發(fā)育提供重要的營養(yǎng)物質，而導致花粉發(fā)育畸形，形態(tài)異常。二核期和三核期，保持系的絨氈層完全降解，而不育系絨氈層細胞在二核期才開始降解，三核期降解不完全，依然能夠觀察到清晰的輪廓和解體的碎片。

以上結果表明，3種D2型細胞質雄性不育小麥花藥的絨氈層細胞在花藥發(fā)育過程中均發(fā)生了細胞凋亡的延遲啟動，從而影響了小孢子的正常發(fā)育，最終導致敗育的形成。因此，推斷D2型細胞質雄性不育小麥絨氈層的延遲降解可能為敗育的細胞學原因。

A—E：Va87B1-706A；F—J：C687B1-706A；K—O：Ju87B1-706A；P—T：706B；A、F、K和P：四分體時期；B、G、L和Q：單核早期；C、H、M和R：單核晚期；D、I、N和S：二核期；E、J、O和T：三核期。E：表皮層；En：內皮層；ML：中層；T：絨氈層；Tds：四分體；Msp：小孢子。在單核晚期，保持系的絨氈層啟動了細胞程序化死亡，絨氈層邊界模糊，并且能夠明顯地觀察到未退化的細胞核碎片及較多的細胞質殘留；不育系在單核晚期時期并沒有觀察到絨氈層的退化，且排列緊密，結構完整。二核期和三核期，保持系的絨氈層完全降解；不育系絨氈層細胞在二核期開始降解，導致了絨氈層PCD的延遲啟動。標尺為50 mm

2.4 花藥絨氈層PCD檢測

為驗證D2型細胞質雄性不育小麥花藥發(fā)育過程中絨氈層PCD的異常過程，對所有材料不同時期的花藥絨氈層內細胞核的DNA片段化程度進行了DNA ladder分析。結果顯示，保持系706B在單核早期開始產生了大量以140 bp倍增的DNA片段啟動PCD，然而3種D2型細胞質雄性不育小麥均在單核晚期才發(fā)生DNA損傷，進一步證實了絨氈層PCD的延遲啟動（圖5）。因此，通過DNA ladder的結果，發(fā)現小麥花藥絨氈層呈現出典型的細胞程序化死亡特征，并且3種D2型細胞質雄性不育小麥凋亡的啟動時間明顯延遲于706B。與細胞學不同，花藥絨氈層的DNA損傷要比形態(tài)學特征早一個時期，即當DNA發(fā)生片段化后，即啟動了花藥絨氈層的PCD，最終導致花藥絨氈層細胞的表型異常。

2.5 活性氧（ROS）含量的測定

為進一步探究D2型細胞質雄性不育小麥絨氈層異常延遲細胞程序化死亡與活性氧代謝的關系，對3種D2型細胞質小麥雄性不育系和其同型保持系在不同時期花藥中O2-、H2O2和MDA含量的測定（圖6）。發(fā)現在所有不育系的花藥內作為活性氧源頭的超氧陰離子自由基O2-的含量和O2-的衍生物H2O2的含量均在單核晚期達到峰值。并在此時期，不育系的O2-和H2O2含量均顯著高于706B（＜0.05）。說明過度積累的活性氧使得花藥內細胞膜的不飽和脂肪酸發(fā)生過氧化反應降解為小分子物質MDA，因此，不育系花藥的MDA含量從單核晚期開始明顯高于706B，在后2個時期，所有不育系與保持系有顯著差異（＜0.05）。說明絨氈層延遲啟動PCD導致活性氧的過量積累，多余的ROS作為一個信號，破壞了抗氧化防御系統(tǒng)的平衡，進而可能造成絨氈層延遲降解的表型特征，最終導致了D2型細胞質雄性不育小麥的敗育。

M：Marker D2000；1：四分體時期；2：單核早期；3：單核晚期；4：二核期；5：三核期

不同小寫字母表示差異顯著（P＜0.05）。下同

2.6 抗氧化酶活性的測定

試驗測定了3種細胞質雄性不育系和相應保持系5個時期花藥SOD、CAT、POD、APX和GPX酶的活性（圖7）。所有不育系SOD、CAT和POD酶活均比保持系高，單核晚期時與保持系之間差異顯著（＜0.05）。但是，在四分體和單核早期時，保持系的APX酶活性比3種不育系低。706B的GPX酶活性整個過程持續(xù)升高，直至三核期酶活性達到最高。所有不育系的GPX酶活性均呈現先下降后升高的態(tài)勢。另外，保持系706B的GPX酶活在單核晚期分別高出3種不育系152%、261%和57%。說明抗氧化酶活性的增加，同樣反映了在小麥雄性不育系內的ROS的過度積累，破壞了抗氧化防御系統(tǒng)的平衡，最終影響了小麥的雄性不育。

2.7 非酶促抗氧化物質含量的測定

通過對小麥花藥發(fā)育過程中AsA和GSH的含量變化的檢測（圖8），3種不育系的AsA含量在單核晚期低于保持系且與保持系706B差異極顯著（＜0.01）。保持系706B的GSH含量先顯著升高，然后連續(xù)降低，單核晚期最高。不育系Va87B1-706A的GSH含量變化趨勢與保持系不同，表現為先降后升，直至三核期達到最高。而C687B1-706A和Ju87B1-706A的GSH含量先升后降再升，且在三核期達到最高。在花藥發(fā)育的單核晚期，保持系706B的GSH含極顯著高于所有細胞質雄性不育系（＜0.01）。說明花藥發(fā)育過程中AsA和GSH的含量劇烈變化，從而導致花粉敗育。

2.8 D2型細胞質雄性不育系抗氧化酶相關基因的表達分析

為研究相關抗氧化酶基因、和的表達模式以及與花藥絨氈層PCD的關系，對小麥細胞質雄性不育系花藥、和進行半定量和定量分析（圖9和圖10，電子附圖1）。其中，熒光定量以保持系706B四分體時期花藥中的基因相對表達量為對照。結果表明，3種不育系在花藥發(fā)育的各個階段的表達量均高于保持系706B。除了Ju87B1-706A，不育系Va87B1-706A、C687B1-706A和保持系706B的表達水平與酶活性的變化基本相同。隨著花藥的發(fā)育，保持系706B的表達水平持續(xù)上調，但是與酶活性的變化有很大的區(qū)別。表達水平在706B上表現出持續(xù)升高的趨勢。在四分體和單核早期，3種不育系的表達水平比保持系高。在二核和三核時期，不育系中的表達卻明顯比保持系低。

從以上結果可以推斷，由于不育系花藥內的活性氧在單核晚期的大量累積，防御系統(tǒng)為保護植物，SOD、CAT和POD酶活顯著上調，而POD的上升又會導致IAA的虧損，物質運輸嚴重受阻。APX和GPX酶的活性以及非酶抗氧化物質AsA和GSH含量的嚴重下調又使抗氧化防御系統(tǒng)更加薄弱，使活性氧進一步積累。而單核晚期ROS和抗氧化系統(tǒng)異常的變化以及相關酶基因的異常表達均與花藥絨氈層紊亂的表型特征表現一致。

3 討論

雜種優(yōu)勢被認為是提高作物產量的最有效途徑，高等植物細胞質雄性不育在作物雜種利用上具有極其重要的價值，并在經濟上為農業(yè)生產提供了巨大利益[37]。作為細胞質雄性不育小麥的重要胞質來源的D2型細胞質具有十分重大的應用價值[8]。本研究利用3種具有完全相同細胞核的不育系，消除核背景的影響，更加詳細客觀地進行比較。因此，該三種不育系和相應保持系是研究質核互作和花粉發(fā)育的一套理想材料。而雄性不育小麥敗育的過程往往伴隨著細胞結構的變化和細胞的程序性死亡[16]。由于不同植物不育的遺傳基礎不同，因而花粉敗育的時期和形式也有所不同[38-41]。對于具有（D2型細胞質）細胞質的異質小麥在長日照條件時雄蕊發(fā)生了心皮化，進而導致了花粉的敗育利用[42]。Murai等[43]也報道了具有細胞質的雄性不育小麥也表現出由雄蕊心皮化引起的完全雄性不育。然而，本研究通過對3種D2型細胞質雄性不育小麥形態(tài)學和細胞學分析并未發(fā)現雄蕊心皮化的發(fā)生，且進一步通過掃描電鏡、體式鏡和絨氈層甲苯胺藍染色詳細地觀察發(fā)現花藥的敗育形式為表皮皺縮、細胞無規(guī)則排列、花藥不開裂且沒有花粉散出的敗育特征。并由DAPI染色和石蠟切片觀察發(fā)現，3種D2型細胞質雄性不育小麥敗育的關鍵時期是單核晚期，細胞學特征表現為小孢子的異常皺縮和絨氈層的延遲降解，且敗育特征穩(wěn)定。造成這種差異的原因可能是核親本不同，質核組合不同，地理環(huán)境的影響以及溫度日照時間不同而引起的。由于植物雄性不育是一個極其復雜的過程，是多種因素共同協(xié)調導致的現象，所以造成這種差異的原因還有待更系統(tǒng)更全面的研究。同時，本研究拓寬了D2型細胞質雄性不育的敗育類型，為進一步研究D2型細胞質雄性不育小麥的細胞學特征提供了借鑒。

圖7 不育系和保持系花藥各時期抗氧化酶活性

圖8 不育系和保持系花藥各時期AsA和GSH的含量

M: Marker D2000; 1: Va87B1-706A; 2: C687B1-706A; 3: Ju87B1 i-706A; 4: 706B

大量學者表明花藥絨氈層的降解是一種典型的細胞程序化死亡過程，其PCD表現為細胞質的濃縮、畸形，最終解體[44-45]。本研究發(fā)現，不育系在單核晚期PCD的延遲啟動可能是導致敗育的細胞學原因。那么，關于細胞質雄性不育小麥花粉和小孢子的敗育及細胞程序化死亡的原因，許多研究表明與活性氧的過量積累和相關抗氧化酶活性相關[46]。Li等[47]證實了紅蓮型水稻減數分裂時期的PCD異常與ROS的過量積累有關。也有研究表明WA352在絨氈層的積累，同時與具有清除ROS作用的COXII共同作用，能夠提前啟動植物的PCD過程[48]。本研究發(fā)現在單核晚期不育系活性氧過量累積并與保持系有顯著差異。在花藥整個發(fā)育過程中，所有不育系的SOD和POD活性均比保持系高，在單核晚期與706B有顯著差異。POD在細胞的活性氧代謝中有重要作用，能夠維持細胞內源庫的平衡。POD酶活性的升高使細胞內代謝水平更加旺盛，IAA嚴重下調，造成雄性不育[49-50]。706B中CAT酶活呈現上升趨勢，而不育系中的CAT酶活卻逐漸下降。這種結果勢必使活性氧產生的速率加快，丙二醛大量積累，最終使花粉及絨氈層的發(fā)育受阻[51]。所有不育小麥的APX酶活在二核期和706B有顯著差異，而GPX酶活在單核晚期與706B有顯著差異。除Ju87B1-706A外，所有不育系AsA和GSH的含量在單核晚期均極顯著低于保持系，且與鄧明華[52]的研究結果一致。因此，進一步明確了D2型細胞質雄性不育小麥敗育的關鍵時期是單核晚期。在此關鍵期，一方面細胞質雄性不育系的ROS含量顯著升高；另一方面APX與GPX酶活性以及非酶物質含量的持續(xù)降低，說明不育系ROS的清除能力降低，單核晚期抗氧化防御系統(tǒng)薄弱，最終造成植物體內活性氧代謝異常。從而推斷延遲啟動絨氈層的導致ROS的過量積累，在抗氧化系統(tǒng)薄弱的影響下，導致了敗育的產生。

圖10 不育系和保持系花藥抗氧化酶基因的qRT-PCR分析

隨著對mtDNA探索的不斷深入，目前，越來越多的研究更注重于分子水平研究植物細胞質雄性不育。周瑋等[53]研究表明表達量高的煙草均能更好的適應各種處理。李莉等[3]發(fā)現，不同發(fā)育時期FBA及Mn-SOD酶活變化與基因表達水平相符，可能直接或間接地造成了小麥的敗育。Noreen等[32]通過提高谷胱甘肽轉移酶基因而提高棉花細胞質雄性不育系的育性。Lim等[54]將和轉入甘薯葉綠體中得到的轉基因植株APX活性與對照有顯著差異。本研究通過對3種D2型細胞質雄性不育小麥的抗氧化酶相關基因的表達分析，除了表達水平與酶活不相符外，、在不育系和保持系花藥中的表達水平與酶活性變化是一致的。差異的原因可能是CAT家族其他基因的表達對該基因有掩蓋作用[55]。本研究在分子生物學水平揭示了D2型細胞質雄性不育小麥的敗育機理，推測PCD的延遲啟動導致活性氧在敗育關鍵時期的過量積累，誘導抗氧化酶基因的異常表達，從而使抗氧化系統(tǒng)紊亂，進而造成了花藥絨氈層的延遲降解，從而引起D2型細胞質雄性不育小麥的敗育。此研究為探索詳細的分子不育機理提供了理論依據，更為植物PCD的研究增添了新的內涵。

4 結論

3種同核異質D2型小麥細胞質雄性不育系的花粉敗育關鍵期均在單核晚期，在此關鍵時期，小孢子表皮的皺縮萌發(fā)孔的閉合以及絨氈層細胞核內DNA損傷延遲至單核晚期的現象導致了活性氧的異常積累，膜脂過氧化加劇，抗氧化代謝異常，以及不育系清除活性氧能力的降低，進而造成絨氈層細胞紊亂的表型特征，使不育系花藥內花粉母細胞自單核晚期因不能滿足必要的營養(yǎng)物質和能量而出現畸形的異常表現，最終導致敗育。

[1] SCHNABLE P S, SPRINGER N M. Progress toward understanding heterosis in crop plants., 2013, 64(4): 71-88.

[2] 甘雨, 龔德華, 楊占烈, 賈先勇, 向關倫, 黃宗洪, 潘建慧. 雜交水稻新品種健優(yōu)388的選育. 貴州農業(yè)科學, 2012, 40(1): 12-14.

GAN Y, GONG D H, YANG Z L, JIA X Y, XIANG G L, HUANG Z H, PAN J H. Breeding of Jianyou 388, a new hybrid rice variety., 2012, 40(1): 12-14. (in Chinese)

[3] 李莉, 張改生, 張龍雨, 李亞鑫. 小麥生理型與遺傳型雄性不育相關基因錳超氧化物歧化酶(Mn-SOD)及果糖1, 6-二磷酸醛縮酶(FBA)的表達分析. 農業(yè)生物技術學報, 2012, 20(3): 225-234.

LI L, ZHANG G S, ZHANG L Y, LI Y XAnalysis of gene expression for manganese superoxide dismutase (Mn-SOD) and fructose-1, 6-bisphosphate Aldolase (FBA) in genetic and physiological male sterility lines of wheat., 2012, 20(3): 225-234. (in Chinese)

[4] 楊天章. 普通小麥種內一種細胞質變異. 作物學報, 1983, 9(4): 217- 219.

YANG T Z. Investigation on the intra specifically cytoplasmic variability in., 1983, 9(4): 217-219. (in Chinese)

[5] 宋喜悅, 馬翎健, 胡銀剛, 奚亞軍, 劉曙東, 何蓓如. D2型山羊草細胞質對普通小麥性狀影響的研究. 西北農林科技大學學報, 2001, 29(6): 27-29.

SONG X Y, MA L J, HU Y G, XI Y J, LIU S D, HE B RStudies on effects of the cytoplasm from D2 type cytoplasm on characters of wheat., 2001, 29(6): 27-29. (in Chinese)

[6] 姚盟, 葉佳麗, 楊智全, 段陽, 蒙立穎, 閆鵬嬌, 劉子涵, 張玲麗, 宋喜悅. 5種細胞質雄性不育小麥敗育的生物學特性及育性恢復. 麥類作物學報, 2015, 35(12): 1676-1684.

YAO M, YE J L, YANG Z Q, DUAN Y, MENG L Y, YAN P J, LIU Z H, ZHANG L L, SONG X YAbortion feature and fertility restoration of five kinds of cytoplasmic male sterile wheat lines., 2015, 35(12): 1676-1684. (in Chinese)

[7] 劉子涵, 蒙立穎, 姚盟, 宋喜悅. 牡山羊草細胞質雄性不育小麥敗育的生物學特性和細胞學研究. 中國農業(yè)大學學報, 2016, 21(5): 1-9.

LIU Z H, MENG L Y, YAO M, SONG X YBiological characteristics and cytological studies on abortion ofcytoplasmic male sterile line., 2016, 21(5): 1-9. (in Chinese)

[8] 張玲麗, 趙鳳霞, 許愛蘭. 牡山羊草()細胞質對普通小麥的效應研究. 西北農業(yè)學報, 2002, 11(4): 19-21.

ZHANG L L, ZHAO F X, XU A L. Effective ofcytoplasm on., 2002, 11(4): 19-21. (in Chinese)

[9] 劉曙東, 宋喜悅. D2型細胞質對普通小麥籽粒品質的影響. 西北農業(yè)大學學報, 1998, 26(2): 7-11.

LIU S D, SONG X Y. Influence of the D2type alloplasm on wheat quality., 1998, 26(2): 7-11. (in Chinese)

[10] 徐潔, 陳升位, 楊德, 何蓓如. D2類山羊草細胞質普通小麥核代換系育性在不同光溫條件下的差異.云南農業(yè)大學學報, 2001,16(2): 111-113.

XU J, CHEN S W, YANG D, HE B RThe reaction to phto-period and temperature of fertility the different of the wheat nucleus substitutional lines of D2genuscytoplasm., 2001, 16(2): 111-113. (in Chinese)

[11] 孫鑫博, 代小梅, 王怡杰, 韓烈保. 植物細胞程序性死亡研究進展. 生物技術通報, 2010(11): 1-6.

SUN X B, DAI X M, WANG Y J, HAN L BAdvances in programmed cell death in plant., 2010(11): 1-6. (in Chinese)

[12] WANG S P, ZHANG G S, SONG Q L, ZHANG Y X, LI Y, GUO J L, CHEN Z, NIU N, MA S C, WANG J WProgrammed cell death，antioxidant response and oxidative stress in wheat flag leaves induced by chemical hybridization agent SQ-1., 2016, 1(11): 76-86.

[13] XIE H T, WAN Z H, LI S, ZHANG Y. Spatiotemporal production of reactive oxygen species by NADPH oxidase is critical for tapetal programmed cell death and pollen development in., 2014, 26(5): 193-201.

[14] 鄭巖, 李江, 楊鵬, 陳惠萍. 水稻糊粉層 PCD 過程細胞形態(tài)觀察. 熱帶作物學報, 2016, 37(2): 298-303.

ZHENG Y, LI J, YANG P, CHEN H PThe cell morphology of rice aleurone layers during programmed cell death., 2016, 37(2): 298-303. (in Chinese)

[15] MENG L Y, LIU Z H, ZHANG L L, HU G, SONG X YCytological characterization of a thermo-sensitive cytoplasmic male-sterile wheat line having K-type cytoplasm of, 2016, 66(5): 752-761.

[16] 張虹, 梁婉琪, 張大兵. 花藥絨氈層細胞程序性死亡研究進展. 上海交通大學學報, 2008, 26(1): 86-90.

ZHANG H, LIANG W Q, ZHANG D B. Research progress on tapetum programmed cell death., 2008, 26(1): 86-90. (in Chinese)

[17] WANG S P, ZHANG G S, SONG Q L, ZHANG Y X, LI Y, GUO J L, CHEN Z, NIU N, MA S C, WANG J WAbnormal development of tapetum and microspores induced by chemical hybridization agent SQ-1 in wheat., 2015, 10: 1-15.

[18] 李建軍, 唐慧娟, 陳艷翠, 陳安國, 黃思齊, 李德芳. 質核互作型紅麻雄性不育系細胞學形態(tài)觀察. 中國麻業(yè)科學, 2015, 37(5): 229-232.

LI J J, TANG H J, CHEN Y C, CHEN A G, HUANG S Q, LI D FCytomorphological observation on cytoplasm nuclear male sterile line of Kenaf., 2015, 37(5): 229-232. (in Chinese)

[19] LI H, YUAN Z, GEMA V B, YANG C Y, LIANG W Q, ZONG J, WILSON Z A, ZHANG D BPresistent tapetal cell1 encodes a phd-finger protein that is required for tapetal cell death and pollen development in rice., 2011, 156(2): 615-630.

[20] 穆蕊, 張祖新, 張方東, 鄭用璉. 玉米CMS-S小孢子敗育過程中的細胞程序性死亡. 作物學報, 2006, 32(5): 666-670.

MU R, ZHANG Z X, ZHANG F D, ZHENG Y LProgrammed cell death during abortion of microspore in S-type cytoplasmically male-sterile maize., 2006, 32(5): 666-670. (in Chinese)

[21] SONG Y L, WANG J W, ZHANG G S, ZHAO X L, ZHANG P F, NIU N, MA S CMicrospore abortion and abnormal tapetal degeneration in a male-sterile wheat line induced by the chemical hybridizing agent SQ-1., 2015, 55(3): 1-12.

[22] GIANLUCA F, RENA B. Oxidative stress and programmed cell death in yeast., 2012, 2(64): 1-21.

[23] JYOTSNA M, SARMAN S. Miltefosine resistance ininvolves suppression of oxidative stress-induced programmed cell death., 2013, 135: 397-406.

[24] OMER M, AHINOAM E, AARON K. Dehydroascorbate: a possible surveillance molecule of oxidative stress and programmed cell death in the green alga., 2014, 202: 471-484.

[25] RUBERTI C, BARIZZA E, BODNER M, NICOLETTA L R, ROBERTO D M, FRANCESCO C, FIORELLA L S, MICHELA Z. Mitochondria change dynamics and morphology during grapevine leaf senescence., 2014, 9(7): 1-12.

[26] 劉齊元, 程元強, 朱肖文, 范東東, 李立新. 雄性不育煙草花蕾中 SOD、POD和CAT活性研究. 中國煙草學報, 2011, 17(5): 34-39.

LIU Q Y, CHENG Y Q, ZHU X W, FAN D D, LI L XStudies on SOD, POD and CAT activities in flower buds of tobacco of cytoplasmic male sterility., 2011, 17(5): 34-39. (in Chinese)

[27] 張悅, 井維霞, 劉修堂, 王濤, 杜寶貞, 幕衛(wèi), 曲愛軍. 不同濃度安打對煙草幼苗MDA、CAT、APX和GPX的影響. 生物災害科學, 2013, 36(1): 22--25.

ZHANG Y, JING W X, LIU X T, WANG T, DU B Z, MU W, QU A JEffects of avatar insecticide on the content of MDA and activities of antioxidant enzyme in tobacco seedlings., 2013, 36(1): 22--25. (in Chinese)

[28] 馬翎健, 丁勤, 王巍, 宋喜悅, 胡銀剛, 奚亞軍, 何蓓如, 劉曙東. 粗厚山羊草細胞質普通小麥核代換系SOD, CAT活性研究. 西北農林科技大學學報, 2002, 30(2): 17-20.

MA L J, DING Q, WANG W, SONG X Y, HU Y G, XI Y J, HE B R, LIU S DStudy on activities of superoxide dismutase and catalase of wheat nuclear substitutional lines incytoplasm., 2002, 30(2): 17-20. (in Chinese)

[29] ZHAO W Z, HAN J R. Effect of H2O2-induced oxidative stress on sclerotial differentiation and AsA-GSH cycle ofQ1 strain., 2016, 6(6): 1-7.

[30] CAROLINE G D M, CAMILA S, JANE M L, CLAUDRIANA L.Increase of Glutathione byscott a cmremr when it is challenged by gallic acid and nisin., 2012, 1(4): 35-43.

[31] 王書平. 小麥生理型雄性不育分子機理研究及其敗育分子模型的建拓[D]. 楊凌: 西北農林科技大學, 2016.

WANG S P. Establishment and extension of the molecular mechanism of physiological male sterility in wheat (.) [D]. Yangling: Northwest A&F University, 2016. (in Chinese)

[32] NOREEN B B, YUAN S, ZHU Y G, WANG XD.Improvements of fertility restoration in cytoplasmic male sterile cotton by enhanced expression of glutathione S-transferase (GST) gene., 2014, 33: 420-429.

[33] LI J X, DAI X M, JIAO Q H. Metabolism of reactive oxygen species in cytoplasmic male sterility of rice by marking upmost pulvinus., 2015, 175: 1263-1269. (in Chinese)

[34] 楊冬業(yè), 張麗珍, 徐淑慶. 西瓜不同生長階段的過氧化氫酶活性研究. 安徽農業(yè)科學, 2012, 40(18): 9604-9606.

YANG D Y, ZHANG L Z, XU S Q. 2012. Study on catalase activity in watermelon during its different growth periods., 2012, 40(18): 9604-9606. (in Chinese)

[35] 黃志明, 陳宇, 吳晶晶, 吳錦程. 硝普鈉對低溫脅迫下枇杷幼果線粒體AsA-GSH循環(huán)代謝的影響. 熱帶作物學報, 2011, 32(8): 1469-1474.

HUANG Z M, CHEN Y, WU J J, WU J CEffects of exogenous sodium nitroprusside on AsA-GSH circulation metabolism in mitochondria of young loquat fruits under low temperature stress., 2011, 32(8): 1469-1474. (in Chinese)

[36] 石曉藝. 5種細胞質雄性不育小麥敗育的特征與生理代謝響應[D]. 楊凌: 西北農林科技大學, 2016.

SHI XＹ.The abortive characteristics nd physiological metabolism responses of five kinds of cytoplasmic male sterile wheats[D]. Yangling: Northwest A&F University, 2016. (in Chinese)

[37] 劉玉平, 王江浩, 趙愛菊, 陳希勇, 李亞軍 , 李爾民, 柳英東. 小麥雜種優(yōu)勢利用體系研究進展. 河北農業(yè)科學, 2009, 13(4): 56-57.

LIU Y P, WANG J H, ZHAO A J, CHEN X Y, LI Y J, LI E M, LIU Y DReview of the heterosis utilization system of wheat., 2009, 13(4): 56-57. (in Chinese)

[38] 陳雪平, 羅雙霞, 軒淑欣, 張成合, 申書興. 茄子雄性不育系花粉敗育的細胞學觀察. 植物遺傳資源學報, 2013, 14(4): 754-758.

CHEN X P, LUO S X, XUAN S X, ZHANG C H, SHEN S XCytological observation of pollen abortion in eggplant male-sterile lines., 2013, 14(4): 754-758. (in Chinese)

[39] SHENG Z H, TANG L Q, SHAO G N, XIE L H, JIAO G A, TANG S Q, HU P S. The rice thermo-sensitive genic male sterility gene tms9: pollen abortion and gene isolation., 2015, 203: 145-152.

[40] CAI H, IMRAN H S, ZHAO H J. Cytological evidences of pollen abortion inAit., 2011, 10(64): 14061-14066.

[41] 張建奎, 董靜, 宗學鳳, 余國東, 李伯群, 戴秀梅, 阮仁武. 溫光敏核雄性不育小麥雄性敗育的細胞學觀察. 植物遺傳資源學報, 2009, 10(3): 411-418.

ZHANG J K, DONG J, ZONG X F, YU G D, LI B Q, DAI X M, RUAN R WCytological analysis of pollen abortion in thermo/ photoperiod-sensitive genic male sterile wheat (L.)., 2009, 10(3): 411-418. (in Chinese)

[42] 徐道欽, 徐軍生. 小麥D2型細胞質光敏雄性不育性的研究. 安徽農業(yè)通報, 2007, 13(19): 53-54.

XU D Q, XU J S. Study on D2type cytoplasm photosensitive wheat male sterility., 2007, 13(19): 53-54. (in Chinese)

[43] MURAI K, TSUTUI I, KAWANISHI Y, IKEGUCHI SDevelopment of photoperiod-sensitive cytoplasmic male sterile (PCMS) wheat lines showing high male sterility under long-day conditions and high seed fertility under short-day conditions., 2007, 159(3): 315-323.

[44] JAKYUNG Y, SUNOK M, YANG-SEOK L, ZHU L, LIANG W Q, ZHANG D B, JUNG K H, GYNHEUNG ADefective tapetum cell death 1 (DTC1) regulates ROS levels by binding to metallothionein during tapetum degeneration., 2016, 170: 1611-1623.

[45] Zhang D D, Liu D, Lü X M, WANG Y, XUN Z L, LIU Z X, LI F L, LU HThe cysteine protease CEP1, a key executor involved in tapetal programmed cell death, regulates pollen development in., 2014, 26(7): 2939-2961.

[46] 劉延忠, 王利民, 李昶. 植物細胞程序性死亡調控機制的研究進展. 山東農業(yè)科學, 2012, 44(11): 58-60.

LIU Y Z, WANG L M, LI C. Advancement of researches of plant programmed cell death regulation and mechanism., 2012, 44(11): 58-60. (in Chinese)

[47] LI S Q, WAN C X, KONG J, ZHANG Z J, LI Y S, ZHU Y G. Programmed cell death during microgenesis in a Hongllian CMS line of rice is correlated with oxidative stress in mitochondria., 2004, 31: 369-376.

[48] LUO D, XU H, LIU Z, GUO J, LI H, CHEN L, FANG C, ZHANG Q, BAI M, YAO N, WU H, WU H, JI C, ZHENG H, CHEN Y, YE S, LI X, ZHAO X, LI R, LIU Y G. A detrimental mitochondrial-nuclear interaction causes cytoplasmic male sterility in rice., 2013, 45: 573-577.

[49] ESMA H A, GULNUR A. Determination of SOD, POD, PPO and CAT enzyme activities inL.., 2016, 11(3): 1-7.

[50] 關天霞, 黨占海, 張建平, 王利民, 趙麗, 田彩萍. 亞麻溫敏型雄性不育系POD活性和內源激素含量比較分析. 甘肅農業(yè)大學學報, 2007, 42(6): 66-70.

GUAN T X, DANG Z H, ZHANG J P, WANG L M, ZHAO L, TIAN C P. Comparing analysis of POD activity and endogenous phytohormone contents in thermo-sensitivity male-sterile flax., 2007, 42(6): 66-70. (in Chinese)

[51] 劉麗華, 李保國, 齊國輝, 郭素萍, 趙志磊. 雄性不育板栗雄花序敗育與幾種酶活性及MDA含量的關系. 林業(yè)科學, 2007, 43(4): 121-124.

LIU L H, LI B G, QI G H, GUO S P, ZHAO Z LRelationships between some protective enzymes activities, MDA content and male inflorescence of male sterile chestnut., 2007, 43(4): 121-124. (in Chinese)

[52] 鄧明華. 辣椒胞質雄性不育的分子生理機制及親緣關系研究[D]. 湖南: 中南大學, 2011.

DENG M H. Studies on molecular and physiological mechanisms of pepper cytoplasmic male sterility and genetic relationship in capsicum [D]. Hunan: Central south university, 2011. (in Chinese)

[53] 周瑋, 周波, 楊雪, 侯思名, 劉明求, 劉飛虎. 轉SOD基因煙草中SOD酶活力對逆境的耐性及其遺傳學特征. 廣西植物, 2006, 26(2): 200-203.

ZHOU W, ZHOU B, YANG X, HOU S M, LIU M Q, LIU F HThe tolerance to adversity stress and its heredity characteristics of SOD activities in transgenic tobacco lines with SOD gene., 2006, 26(2): 200-203. (in Chinese)

[54] LIM S, KIM Y H, KIM S H, KWON S Y, LEE H S, KIM J S, CHO W Y, PAEK H Y, KWAK S SEnhanced tolerance of transgenic sweetpotato plants that express both Cu Zn SOD and APX in chloroplasts to methyl viologen-mediated oxidative stress and chilling., 2007, 19(3): 227-239.

[55] LUNA C M, PASTORI G M, DRISCOLL S, GROTEN K, BERNARD S, FOYER CDrought controls on H2O2accumulation, catalase (CAT)activity and CAT gene expression in wheat., 2005, 56(411): 417-423.

（責任編輯 李莉）

附表1 相關抗氧化基因熒光定量PCR引物

Table S1 Primers used for fluorescent quantitative PCR assays of antioxidant genes related to CMS

引物Primer序列Sequence (5′-3′)退火溫度Tm(℃) SOD-fSOD-rAGAACCTCAAGCCTATCAGCGACAAATCACGCAAGCACT60 CAT-fCAT-rTGCCTGTGTTTTTTATCCGAACCGTCCATGTGCCTGTAGT62 APX-fAPX-rGTTCATCCCTGGAAGACGCAGAGGGTCACGAGTCCA64 Actin-fActin-rCTCCCTCACAACAACCGCTACCAGGAACTTCCATACCAAC62

1：四分體時期；2：單核早期；3：單核晚期；4：二核期；5：三核期

M: Marker D2000; 1: Tetrad stage; 2: Early uninucleate stage; 3: Late uninucleate stage; 4: Binucleate stage; 5: Trinucleate stage

附圖1 不育系和保持系、、和的RT-PCR擴增電泳檢測結果

Supplemental Fig. 1andgenes PCR amplification of cDNA sequences in anthers fromIAMSLs during different development stages

Tapetal Programmed Cell Death, Antioxidant Response and Oxidative Stress in Wheat Anthers Associated withD2-type Cytoplasmic Male-sterility

LIU ZiHan, SHI XiaoYi, YAN PengJiao, DUAN Yang, GENG XingXia, YE JiaLi, LI Sha, YANG XueTong, ZHANG GaoMing, JIA YuLin, ZHANG LingLi, SONG XiYue

(College of Agronomy, Northwest A&F University, Yangling 712100, Shaanxi)

D2type cytoplasm is the considerable source of cytoplasm of cytoplasmic male sterility (CMS). It is of great value to further study the cytological characteristics and mechanism of physiology and molecules in heterologous cytoplasmic male sterility (CMS) wheat for the heterosis utilization.To elucidate the abortive cytological characteristics of D2-type cytoplasmic male-sterility, the phenotypes of anthers, microspores, tapetum were analyzed by scanning electron microscope, DAPI staining, and paraffin section, as well as detectingtapetal apoptotic DNA fragmentation by DNA ladder. To further clarify that the D2-type cytoplasmic male-sterility is associated with abnormal tapetal programmed cell death and reactive oxygen species metabolism, the physiological properties, including the contents of reactive oxygen species (ROS), activities of ROS-scavenging enzymes (superoxide dismutase (SOD), peroxidase, catalase (CAT), ascorbate peroxidase (APX), and glutathione peroxidase (GPX)) and nonenzymatic antioxidants (ascorbic acid (AsA) and glutathione (GSH)), and the expression levels of antioxidant-associated genes from three D2-types of cytoplasmic male sterile lines (Va87B1-706A, C687B1-706A, Ju87B1-706A) and the maintainer 706B of wheat were compared.Compared with the maintainer line, at the trinucleate stage, the inner epidermal cells of D2-type cytoplasmic male-sterility accumulated less and more unconsolidated Ubisch bodies and the epidermis of microspores were extremely scabrous. The abortive pivotal period of pollen was initiated at the later uninucleate stage, which the delayed degradation of tapetum and abnormal structure of the microspore were related to the serious oxidative stress, disorganized activities of antioxidant enzymes, down-regulated contents of non-enzymatic antioxidants and significant difference of the expression level of some related enzyme genes in the three D2-types of cytoplasmic male sterile lines. Additionally, the expression level ofandgenes was approximated to the activities of enzymes. However, there existed a distinct difference between the expression level ofgene and the activities of enzymes resulting from multiple-gene encodinggene.It was deduced that delay of start-up of the tapetal PCD (Programmed Cell Death) results from excessive ROS accumulation, and then destroyed the balance of the antioxidant defense system, and eventually led to the abortion of D2-types of cytoplasmic male sterile lines. Therefore, in order to further explore the abortive mechanism of cytoplasmic male sterility (CMS) wheat, this study provided an important theoretical basis and application significance.

wheat (L.); D2-types cytoplasmic male sterility; tapetal programmed cell death; ROS metabolism; gene expression

2017-04-10；接受日期：2017-06-02

國家自然科學基金（31271792）

劉子涵，E-mail：liuzihan643@163.com。通信作者宋喜悅，E-mail：songxiyue@126.com