水稻矮桿突變體的細胞學特征及基因定位研究

郝小花，向玉婷，曾 孟，楊友偉，謝 鵬，李 密，李東屏，田連福*

(1.湖南文理學院生命與環(huán)境科學學院，中國湖南常德415000；2.湖南師范大學生命科學學院作物不育資源創(chuàng)新與利用湖南省重點實驗室，中國湖南長沙410081)

株高是水稻最重要的性狀之一，理想的株高在水稻育種中占有重要地位，中等矮化有利于水稻品種耐肥、抗倒伏和高產。矮化突變體的分離和鑒定有助于闡明植物生長發(fā)育的分子機制[1]。研究表明，水稻的高度受株高主效基因的控制，同時受基因修飾、基因抑制等多種因素的影響[2]。這些基因的變化可以導致植物細胞水平的變化，如節(jié)間細胞縮短、數(shù)量減少等，其中一些與信號轉導有關，另一些與相關激素的生物合成有關[3]。已報道的涉及水稻株高的基因有很多，如SD1[4]、SLR1[5]、EUI1[6]、D11[7]、BRD1[8]、D2[9]、D33[10]、D53[11]等，這些基因在赤霉素(gibberellin，GA)或油菜素內酯(brassinolide，BR)的信號轉導和合成中起著非常重要的作用。

GA是調節(jié)植物生長發(fā)育不可缺少的植物激素之一，不僅在植物株高調控過程中具有重要的作用，而且在調控種子萌發(fā)、下胚軸伸長、葉片伸展及花、果實、種子發(fā)育等眾多生理過程中也起著重要的作用[12]。近年來，隨著擬南芥、水稻等模式植物的應用，GA在高等植物中的生物合成調控、信號轉導途徑等方面的研究取得突破性進展[13]。Kobayashi等[14]對水稻日本晴的營養(yǎng)組織和生殖組織中的內源GA進行了詳細分析，一共鑒定到13種GA，不同種類的GA可能參與不同的生理過程。水稻的兩個矮化品種Tan-ginbozu(dx突變體)和Waito-C(dy突變體)的地上部分均表現(xiàn)出GA缺陷的表型[15]，外源施加貝殼杉烯可以刺激Tanginbozu地上部分的伸長[16]。然而，尚未確定由Dx和Dy基因所控制的GA生物合成途徑的具體步驟[17]。d35是一個與GA合成相關的突變體，D35基因編碼GA合成初始階段的合成酶——貝殼杉烯氧化酶(kaurene oxidase，KO)，該基因的突變使GA前體物質的產生減少，最終出現(xiàn)植株矮化的現(xiàn)象[3]。sbl(shortened basal internodes)是一個半矮桿突變體，SBL基因編碼GA合成關鍵酶GA2ox(gibberellin 2-oxidase)，該基因的突變導致GA合成受阻，從而使植株出現(xiàn)矮化現(xiàn)象[18]。SBL主要在水稻莖稈的基部表達，所以矮化現(xiàn)象集中出現(xiàn)在水稻基部節(jié)間位置，從而使水稻獲得更強的抗倒伏能力。水稻EUI(elongated uppermost internode)基因編碼一個可以失活GA的P450蛋白酶，該基因突變后導致活性GA在水稻最上節(jié)間過度積累，從而引起長頸穗的植株表型[19]。

GA與其他植物生長物質之間存在相互作用，如細胞分裂素(cytokinin，CK)。Fleishon 等[20]證實在番茄花青素積累過程中，施加CK可以促進花青素的積累，施加GA則抑制花青素的積累，而同時施加GA和CK時，GA會完全抑制CK的作用。另外，在番茄下胚軸伸長過程中，外施GA可以促進下胚軸的伸長，而外施CK對下胚軸的伸長沒有顯著影響，當同時施加GA和CK時，則完全抑制由GA所引起的下胚軸伸長。此外，在番茄座果形成及發(fā)育過程中，生長素應答因子ARF7(auxin response factor 7)參與GA誘導的細胞膨大與生長素誘導的細胞分裂[21]。因此，GA所參與的植株形態(tài)的改變不僅是由于其含量的變化，也涉及GA與其他生長物質相互作用而進行的多重調節(jié)。

課題組前期從水稻Kitaake(Oryza sativa L.japonica)甲基磺酸乙酯(ethyl methane sulfonate，EMS)誘變的M2群體中分離到一株超矮桿突變體ag1(ai-gan1)，表型分析顯示突變體節(jié)間長度及其所占整株比例發(fā)生顯著改變。為了對導致植株矮化的基因進行精確定位，我們利用MBS(mapping-bysequencing)方法對超矮桿突變體與野生型(wild type，WT)雜交F1代自交產生的F2后代進行全基因組序列比對分析，共篩選到6個導致蛋白質可讀框發(fā)生改變的突變基因，并利用CRISPR-Cas9技術最終確定導致植株矮化的基因位于1號染色體，編碼赤霉素3β-羥化酶2(gibberellin 3-betahydroxylase 2，GA3ox2)，是D18的等位基因。d18突變體的發(fā)現(xiàn)較早[22]，但其細胞學研究卻較少。本文通過對突變體節(jié)間和劍葉的細胞學特性進行觀察及對部分細胞分裂和伸長相關基因的表達水平進行檢測，發(fā)現(xiàn)OsGA3ox2不僅抑制倒二節(jié)間細胞的伸長，而且影響劍葉和節(jié)間細胞縱向、橫向的分裂能力，在水稻地上部分的形態(tài)發(fā)育中起著重要作用。

1 材料與方法

1.1 實驗材料

水稻ag1矮桿突變體來源于本實驗室構建的EMS誘變的M2突變體庫，其遺傳背景為粳稻品種Kitaake。經本實驗室長沙和海南實驗基地的多代自交和選擇，突變體表型穩(wěn)定。本實驗使用的引物和 LB、N6D、AB、AAM、2N6-AS培養(yǎng)基試劑購自上海生工生物工程股份有限公司。N6D、AB、AAM、2N6-AS培養(yǎng)基配方參考文獻[23]。RNA提取試劑Trizol、反轉錄酶、高保真酶、內切酶等購自Invitrogen 公司(美國)。q-RT-PCR(real-time quantitative polymerase chain reaction)試劑為SYBR GreenⅠ，購自Promega公司(美國)。實驗所用Cas9載體、sg-RNA載體由北京大學瞿禮嘉教授提供；農桿菌EHA105、大腸桿菌TOP10為本實驗室保存。

取顆粒飽滿、活性一致(同時期種植、同時收獲)的野生型和突變體種子同期催芽、播種，比較野生型和突變體各個時期地上部分的差異，并記錄各個時期相關性狀的數(shù)據(jù)，主要包括株高、分蘗數(shù)、劍葉長、劍葉寬、千粒重、各節(jié)間長等。將ag1矮桿突變體與野生型雜交，獲得F1代，F(xiàn)1自交后獲得的F2代。材料種植成8行×8株規(guī)格的小區(qū)，行株距為20 cm×20 cm，成熟時從地面至穗頂測量株高。節(jié)間縱向細胞總數(shù)的估算方法為:節(jié)間的平均長度/細胞平均長度。

1.2 MBS全基因組測序定位

分別以野生型植株20株、自交F2后代中100株正常高度和100株矮桿植株的葉片為材料，構建3個DNA庫。利用MBS_QTL全基因組測序對3個DNA庫進行高通量測序，篩選候選基因。將過濾后的高質量reads與參考基因組序列數(shù)據(jù)庫(http://www.ricedata.cn/gene/)進行比對，對篩選出的單核苷酸多態(tài)性(single nucleotide polymorphism，SNP)位點進行鑒定與分型，根據(jù)篩選出的SNP位點在3個DNA庫中出現(xiàn)的頻率統(tǒng)計，篩選得到符合比例的SNP位點。測序和SNP位點篩選由上海歐易生物醫(yī)學科技有限公司完成。

1.3 CRISPR-Cas9基因敲除載體構建

從NCBI數(shù)據(jù)庫網(wǎng)站(https://www.ncbi.nlm.nih.gov/genome/)下載初步篩選的候選基因的編碼序列(coding sequence，CDS)，按照 CRISPR-Cas9 系統(tǒng)靶序列設計的要求為每個候選基因設計兩個不同的靶序列[24]，并委托上海生工生物工程股份有限公司合成靶序列(表1)。將正反向引物稀釋至10 μmol/L后混合，PCR退火后再稀釋成0.1 μmol/L。用BasⅠ酶切sg-RNA原始載體，體系為sg-RNA 2 μL、buffer 2 μL、BasⅠ內切酶 2 μL、ddH2O 14 μL，37 ℃酶切過夜，純化酶切產物。用T4 DNA連接酶連接sg-RNA和PCR退火產物，體系為ddH2O 6 μL、buffer 1 μL、sg-RNA 酶切產物 1 μL、PCR 退火產物1 μL、T4連接酶1 μL，22 ℃水浴 6～8 h。用42 ℃熱激法將連接產物轉化大腸桿菌感受態(tài)細胞TOP10，隨后將轉化產物涂布于含卡那霉素(kanamycin，Kan)的LB平板，篩選陽性克隆并測序，將測序正確的陽性克隆擴繁，小量提取質粒并與Cas9載體進行LR重組反應，體系包含ddH2O 7 μL、陽性質粒 1 μL、Cas9 載體 1 μL、LR mix 1 μL，25 ℃水浴 1 h，隨后加入 1 μL protein K，37 ℃水浴20 min終止反應。用42℃熱激法將LR重組反應的產物轉化大腸桿菌感受態(tài)細胞TOP10，隨后將轉化產物涂布于含壯觀霉素(spectinomycin，Spe)的LB平板，篩選陽性克隆并測序，將測序正確的陽性克隆擴繁，小量提取質粒，得到CRISPR-Cas9基因敲除載體。

表1 Os01g0177400基因CRISPR-Cas9的靶位點引物Table 1 Primers used for Os01g0177400 knockout with CRISPR-Cas9 system

1.4 農桿菌轉化

電轉杯預先用無水乙醇浸泡，再依次用超純水清洗5次、無水乙醇清洗5次，最后在超凈工作臺內置于干凈濾紙上晾干。將EHA105感受態(tài)與CRISPR-Cas9基因敲除載體置于冰上溶解3 min，同時用手套包住電轉杯置于冰上。取1 μL CRISPRCas9基因敲除載體質粒(稀釋100倍)混勻于EHA105感受態(tài)中，冰浴5～15 min。在無菌的EP管中加入1 mL LB液體培養(yǎng)基，置于37℃水浴鍋中預熱。將混有質粒的感受態(tài)轉移至電轉杯，并于電轉儀上進行轉化。電轉后加入300 μL已預熱的LB液體培養(yǎng)基，混勻，并將菌液完全吸出至無菌EP管中，37℃水浴1 h。吸取100 μL菌液涂布于具有Spe抗性的LB平板，30℃倒置培養(yǎng)30～48 h。

1.5 轉基因植株的獲得

將在AB培養(yǎng)基活化3 d的農桿菌EHA105(含轉化質粒)挑取到100 mL AAM液體培養(yǎng)基(加入1 g/L乙酰丁香酮1.5 mL)中，180 r/min振蕩培養(yǎng)30 min，使其OD600為0.1～0.15。將水稻成熟胚誘導7 d的愈傷組織挑入上述AAM菌液中，輕搖2～3 min，棄菌液，于超凈工作臺晾干；將晾干的愈傷組織轉移到2N6-AS培養(yǎng)基中(墊一張無菌濾紙)，于25℃恒溫培養(yǎng)箱中暗培養(yǎng)3 d，隨后用無菌水洗滌8次，每次3 min，再用羧芐青霉素水洗滌5 min，置于無菌濾紙上晾干；將晾干的愈傷組織接種到含羧芐青霉素和潮霉素的篩選培養(yǎng)基上，30℃連續(xù)光照培養(yǎng)至新愈傷出現(xiàn)，將新愈傷轉移到分化培養(yǎng)基上(篩選和分化過程中每7～10 d更換一次培養(yǎng)基)，待分化出芽和根即可將其移到生根瓶中培養(yǎng)[23]。7～10 d后，打開瓶蓋，加入少量無菌水繼續(xù)培養(yǎng)1周，然后將其轉移到土壤中，使其繼續(xù)生長。

1.6 石蠟切片的制作與觀察

石蠟切片的制作參照文獻[25]的方法進行，主要步驟如下:將材料于卡諾氏固定液固定，隨后用不同濃度的乙醇進行脫水；在通風櫥內依次用乙醇-二甲苯混合液與純二甲苯進行透明，每次30 min；用二甲苯、石蠟的混合液和純石蠟對材料進行逐步浸蠟；向折好的紙盒中倒入含有材料的純石蠟，靜置，待石蠟凝固即可進行常規(guī)切片；將切好的蠟帶用載玻片進行展片；將晾干后的載玻片浸泡于純二甲苯中，5 min，重復兩次；滴加中性樹膠于載玻片上，用鑷子夾上蓋玻片，制成永久裝片。用普通光學顯微鏡進行切片觀察。

1.7 實時熒光定量PCR

以野生型和突變體分蘗期未伸長的新葉為材料，利用Trizol試劑提取總RNA。根據(jù)RNA濃度計算反轉錄所需的RNA體積。cDNA第一鏈的合成采用Thermo Fisher Scientific公司的反轉錄試劑，全程配戴口罩和乳膠手套。利用Primer 5軟件根據(jù)目的基因CDS序列設計引物(表2)，實時熒光定量PCR總反應體積為30 μL，試劑為SYBR GreenⅠ(Promega)，PCR 反應在 Quant Studio 5(Life Technologies公司，美國)中完成，實驗重復3次[25]。

表2 實時熒光定量PCR引物Table 2 Primers used for real-time quantitative PCR analysis

1.8 數(shù)據(jù)處理

實驗所得統(tǒng)計數(shù)據(jù)和實時熒光定量PCR數(shù)據(jù)用Excel 2003軟件進行初步處理，采用 SPSS 19.0軟件進行單因素方差分析。采用Origin 8.0軟件制圖。

2 結果與分析

2.1 突變體地上部分的形態(tài)分析

分蘗期野生型和突變體植株的表型觀察結果顯示，野生型和突變體的劍葉長分別為(23.81±3.51)cm、(10.72±2.42)cm，突變體劍葉長只有野生型的 45.02%(圖 1A)，而劍葉寬達到(1.29±0.14)cm，是野生型(1.04±0.13)cm 的 124.04%(圖 1C)；突變體的平均分蘗數(shù)為9.31個，野生型為6.82個，突變體的分蘗數(shù)顯著增多(圖1B，D)。此外，成熟期野生型和突變體籽粒千粒重的統(tǒng)計數(shù)據(jù)顯示，突變體千粒重為(22.72±1.02)g，野生型為(26.21±1.36)g，突變體顯著低于野生型(圖1E)。以上形態(tài)分析結果表明，AG1基因的突變對水稻地上部分的生長發(fā)育有較大影響。

圖1 野生型和突變體的農藝性狀(A)成熟期劍葉；(B)分蘗期植株形態(tài)；(C)劍葉寬；(D)分蘗數(shù)；(E)千粒重。標尺=2 cm；數(shù)據(jù)為平均值±標準差，**表示差異極顯著(P＜0.01)。Fig.1 Agronomic traits of WT and mutant plants(A)Mature flag leaves；(B)Tillering plant morphology；(C)Flag leaf width；(D)Tillering number；(E)1 000-grain weight.Bar=2 cm.Data are presented as mean ± standard deviation.Asterisks indicate significant difference from the WT plants at**P＜0.01 by Student’s t test.

此外，我們分析了不同時期突變體的株高，發(fā)現(xiàn)幼苗期、分蘗期和成熟期的突變體均表現(xiàn)出極度矮化(圖2A)，其中成熟期野生型的平均株高為(52.43±5.36)cm，而突變體為(18.25±3.12)cm，僅為野生型的34.81%(圖2B)。進一步統(tǒng)計成熟期植株的穗長以及倒一節(jié)間、倒二節(jié)間和倒三節(jié)間的長度，結果顯示野生型各部分的平均長度依次為(10.91±2.32)cm、(17.13±2.12)cm、(9.91±2.36)cm和(3.64±1.27)cm，而突變體各部分的平均長度依次為(6.12±1.24)cm、(6.91±2.32)cm、(1.11±0.52)cm和(0.37±0.12)cm，均顯著縮短(圖 2C)。通過計算各部分長度所占整株的比例，我們發(fā)現(xiàn)野生型穗長以及倒一節(jié)間、倒二節(jié)間和倒三節(jié)間的長度所占整株比例分別為20.81%、32.67%、18.90%和6.94%，而突變體則為33.53%、37.86%、6.08%和2.03%，突變體穗長和倒一節(jié)間所占比例有所增加，而倒二節(jié)間和倒三節(jié)間所占整株的比例顯著降低。以上數(shù)據(jù)表明AG1基因突變影響穗長及各節(jié)間的伸長，導致植株變矮及各節(jié)間變短且占整株比例發(fā)生改變。

圖2 野生型和突變體的植株形態(tài)及各部分長度統(tǒng)計(A)幼苗期、分蘗期和成熟期的植株形態(tài)；(B)成熟期植株高度的統(tǒng)計；(C)穗長及各節(jié)間長度的統(tǒng)計。標尺=2 cm；數(shù)據(jù)為平均值±標準差，**表示差異極顯著(P＜0.01)。Fig.2 Morphology of WT and mutant plants and their length statistics(A)Mutant morphology at different stages；(B)Statistics of plant height at maturity stage；(C)Statistics of internode and panicle lengths.Bar=2 cm.Data are presented as mean±standard deviation.Asterisks indicate significant difference from the WT plants at**P＜0.01 by Student’s t test.

2.2 突變體節(jié)間細胞的特征分析

為了解突變體節(jié)間長度變化的原因，我們對倒一節(jié)間和倒二節(jié)間的中部進行切片觀察，并對細胞大小和數(shù)目進行測量與統(tǒng)計。節(jié)間縱切面的分析結果顯示:突變體和野生型倒一節(jié)間細胞的平均長度分別為(109.85±12.21)μm和(126.14±11.24)μm，突變體倒一節(jié)間的細胞長度略小于野生型，約為野生型的87.09%(圖3A，D)；突變體倒二節(jié)間的細胞長度為(28.32±5.31)μm(n=50)，顯著小于野生型的(106.23±13.42)μm(n=50)，僅為野生型細胞長度的26.66%(圖3B，E)。倒二節(jié)間縱列細胞總數(shù)目的估算結果顯示，突變體縱向細胞數(shù)目顯著減少，約為453個，但野生型約為1 178個?？v切面的細胞學分析結果說明，AG1突變不僅導致倒二節(jié)間細胞長度變短，而且縱向細胞數(shù)目也顯著減少。節(jié)間的橫切結果顯示:突變體的莖稈壁厚度顯著大于野生型(圖3G)；突變體倒二節(jié)間橫向薄壁細胞的平均層數(shù)為16層，顯著多于野生型的11層(圖3C，F(xiàn))。這些結果表明，AG1基因的突變顯著抑制倒二節(jié)間細胞的伸長和縱向分裂，促進倒二節(jié)間細胞的橫向分裂，從而導致倒二節(jié)間顯著縮短，莖稈壁顯著增厚。

圖3 野生型與突變體倒一節(jié)間和倒二節(jié)間的細胞形態(tài)(A)倒一節(jié)間縱切面細胞形態(tài)；(B)倒二節(jié)間縱切面細胞形態(tài)；(C)倒二節(jié)間橫切面細胞形態(tài)；(D)倒一節(jié)間縱向細胞的長度和寬度統(tǒng)計；(E)倒二節(jié)間縱向細胞的長度和寬度統(tǒng)計；(F)倒二節(jié)間薄壁細胞層數(shù)的統(tǒng)計；(G)倒二節(jié)間橫切面莖厚度的統(tǒng)計。標尺=100 μm；**表示差異極顯著(P＜0.01)。Fig.3 The peduncle and penultimate internode sections of the WT and mutant plants(A)Cell morphology of the peduncle internode in longitudinal section；(B)Cell morphology of the penultimate internode in longitudinal section；(C)Cell morphology of the penultimate internode in transverse section；(D)Cell length and width statistics of the peduncle internode in longitudinal section；(E)Cell length and width statistics of the penultimate internode in longitudinal section；(F)Cell layer statistics of the penultimate internode；(G)Stem thickness statistics of the penultimate internode.Bar=100 μm.Asterisks indicate significant difference from the WT plants at**P＜0.01 by Student’s t test.

2.3 突變體劍葉細胞的特征分析

由于ag1突變體和野生型的劍葉寬度具有明顯差異，因此對野生型和突變體劍葉中部的橫切面進行了分析，結果顯示:突變體劍葉中維管束數(shù)目顯著增加，其中橫切面大維管束數(shù)目平均比野生型多2.43個，小維管束數(shù)目平均增加9.52個(圖4A)，但突變體中維管束周圍的縱向薄壁細胞數(shù)目顯著減少(圖4B)。

圖4 野生型和突變體的劍葉中部橫切面細胞形態(tài)(A)劍葉中部橫切面，黑色箭頭標示兩個大維管束之間的小維管束；(B)葉片中間大葉脈的放大圖，紅色箭頭標示葉脈中的薄壁細胞。標尺=500 μm。Fig.4 Cell morphology of the transverse section at the middle part of flag leaves of the WT and mutant plants(A)The transverse section of the middle part of flag leaves.The black arrow marks the small vascular bundles between two large vascular bundles；(B)The enlarged picture of the large vascular in the middle of the leaves.The red arrow marks the parenchyma cells in veins.Bar=500 μm.

2.4 水稻AG1基因的定位與功能驗證

在ag1突變體和野生型雜交獲得的F1代中，全部植株的高度均正常，但在F1自交后獲得的F2群體中，高、矮植株數(shù)量分別為354株和116株，符合3︰1分離比，表明該突變體矮桿性狀是由單基因控制的隱性遺傳。為了確定矮桿性狀的控制基因，利用MBS全基因組測序對突變基因進行定位，共篩選到6個導致蛋白質可讀框發(fā)生改變的SNP位點，這6個SNP位點分別分布于1號染色體的Os01g0156000基因和Os01g0177400基因、8號染色體的Os08g0352100基因和Os08g-0447100基因、9號染色體的Os09g0337300基因以及11號染色體的Os11g0200000基因(表3)。參考中國水稻網(wǎng)(http://www.ricedata.cn/gene/)數(shù)據(jù)，從6個基因中篩選出與下胚軸伸長功能相關的Os01g0156000基因和與GA合成功能相關的Os01g0177400基因作為候選基因進行功能驗證。

利用CRISPR-Cas9系統(tǒng)對Os01g0156000和Os01g0177400兩個候選基因分別進行基因編輯，在Os01g0156000和Os01g0177400兩個基因的外顯子區(qū)域分別設計兩個靶序列，構建CRISPRCas9基因敲除載體，通過農桿菌介導分別得到Os01g0156000和Os01g0177400基因敲除的純合體植株。基因敲除陽性植株的表型顯示，Os01g-0156000基因敲除的陽性植株中沒有出現(xiàn)矮稈性狀的植株，而Os01g0177400基因敲除的陽性植株中存在大量的矮化植株。將Os01g0177400基因敲除的純合陽性植株分別命名為ag1-2(圖5A)、ag1-3a和 ag1-3b(圖 5B)，ag1-2、ag1-3a和 ag1-3b均表現(xiàn)出與ag1類似的矮桿性狀(圖5C)，因此，確定AG1位于1號染色體的Os01g0177400基因座，該基因編碼一個赤霉素3β-羥化酶，參與GA的合成代謝。

圖5 Os01g0177400基因靶序列位點及陽性植株形態(tài)(A)第一個靶序列位點；(B)第二個靶序列位點；(C)基因敲除陽性植株形態(tài)。標尺=2 cm。Fig.5 Target sequence loci for the gene Os01g0177400 and morphology of positive plants(A)Firsttargetsequencelocus；(B)Secondtargetsequenceloci；(C)Themorphologyofpositiveplantsthroughgeneknockout.Bar=2cm.

2.5 突變體中與細胞分裂和伸長相關基因的表達量分析

為了了解AG1基因突變后突變體節(jié)間和劍葉細胞大小與數(shù)目變化的原因，我們進一步對野生型和突變體地上部分細胞伸長和分裂相關基因的表達量進行了實時熒光定量PCR分析。根據(jù)已知基因的功能，選取影響細胞數(shù)目的細胞周期相關基因 OsCDKB2；1[26]、CYCA3；1[27]和 CYCB2；2[26]及與細胞伸長相關的 OsE2Fl、OsMCM3[28]、OsEXPA16[29]和OsXTH28[30]基因進行表達量檢測。圖6結果顯示，OsCDKB2；1、OsE2Fl、OsMCM3、OsEXPA16和OsXTH28基因的表達量在突變體中顯著降低，說明AG1的突變對細胞分裂和伸長相關基因的表達都有一定的影響。

圖6 突變體和野生型中與細胞分裂和伸長相關的部分基因的表達量變化*表示差異顯著(P＜0.05)，**表示差異極顯著(P＜0.01)。Fig.6 Expression levels of genes related to cell division and elongation in the mutant and WT plantsAsterisks indicate significant difference from the WT plants at*P＜0.05 and**P＜0.01 by Student’s t test.

3 討論

Os01g0177400編碼赤霉素3β-羥化酶，該酶是植物體內催化產生活性GA的關鍵酶，在GA的合成過程中具有重要的作用。早在2001年，人們就在水稻基因組中克隆到兩個赤霉素3β-羥化酶基因OsGA3ox1和OsGA3ox2[31]，研究表明這兩個基因具有完全不同的表達模式，分別主要在植物的營養(yǎng)器官和生殖器官中表達。這兩個基因的表達方式可能與植物營養(yǎng)器官和生殖器官中發(fā)現(xiàn)的優(yōu)勢生物活性GA的差異有關[14，32]，其通過調控營養(yǎng)器官和生殖器官中不同形式的活性GA從而調控植株不同部位的生長和發(fā)育。

表3 MBS_QTL全基因組篩選的SNP位點Table 3 SNP loci screened from MBS_QTL

在前期研究中，我們將野生型粳稻品種Kitaake用EMS進行誘變，在誘變的M2代群體中發(fā)現(xiàn)了超矮桿突變體ag1(圖1A)。將突變體ag1與野生型雜交，獲得的F1代群體全部表現(xiàn)為高桿，而F1自交后代F2群體中高桿和矮桿的數(shù)量比為3︰1，說明此矮桿性狀為單基因控制性狀。通過MBS_QTL技術篩選到6個候選的突變基因，進一步通過CRISPR-Cas9技術對6個候選基因進行基因敲除驗證，發(fā)現(xiàn)只有Os01g0177400基因敲除的后代中出現(xiàn)矮桿性狀的植株(圖5C)，由此確定突變基因為Os01g0177400。本研究的統(tǒng)計分析發(fā)現(xiàn)，突變體成熟期的高度只有野生型的34.81%，突變體倒二節(jié)間和倒三節(jié)間的長度僅占整株高度的6.08%和2.03%，所占整株比例顯著低于野生型，是突變體株高變矮的主要原因(圖2C)。同時，突變體劍葉變寬，分蘗數(shù)也顯著增加。通過對倒一節(jié)間、倒二節(jié)間和劍葉的中部進行切片分析發(fā)現(xiàn)，突變體倒一節(jié)間細胞的長度有所變短，但較野生型沒有顯著差異；而倒二節(jié)間細胞的長度較野生型顯著變短；倒二節(jié)間縱向細胞總數(shù)目和劍葉中大維管束周圍的縱向薄壁細胞數(shù)目都顯著減少，從而導致突變體倒二節(jié)間長度顯著縮短、劍葉大維管束變小。同時，突變體劍葉中小葉脈的數(shù)目和倒二節(jié)間橫向薄壁細胞層數(shù)較野生型都顯著增多，從而導致突變體劍葉變寬、倒二節(jié)間的莖壁厚度增加(圖3～4)。細胞形態(tài)的觀察結果說明，AG1基因突變不僅影響了倒二節(jié)間細胞的伸長，而且對倒二節(jié)間和劍葉細胞的縱向和橫向分裂能力也具有一定影響。

已有研究表明，GA與其他的生長物質之間存在相互作用并共同調控植株的發(fā)育，如與CK的拮抗作用等。因此，水稻ag1突變體倒二節(jié)間和劍葉的細胞數(shù)目及形態(tài)的變化可能與OsGA3ox2的表達模式有關，OsGA3ox2的突變導致不同部位的活性GA含量發(fā)生改變，從而影響相互作用的其他生長物質發(fā)生變化，進而導致細胞的形態(tài)和分裂能力都發(fā)生變化。在本研究的超矮桿突變體ag1中，與細胞分裂相關的基因OsCDKB2；1、OsE2Fl、OsMCM3以及與細胞伸長相關的基因OsEXPA16和OsXTH28的表達量都較野生型發(fā)生了明顯的變化，說明AG1基因在影響GA合成的同時，還可能與其他生長物質交叉互作從而嚴重影響植株的形態(tài)變化，因此，研究GA與其他生長物質的交叉互作，將為闡明植物形態(tài)發(fā)育提供重要的線索。