趙輝 張麗麗 郭靜遠 霍姍姍 賀萍萍 夏啟玉 符冬妹 郭安平

摘 要 狗尾草是一種新型的模式植物，符合作為模式植物的所有要求，是C4光合作用的優(yōu)秀模型，也是研究非生物脅迫耐受性很好的模型。在狗尾草高效遺傳轉(zhuǎn)化技術(shù)的基礎(chǔ)上，以適宜轉(zhuǎn)化的狗尾草ME34品系非轉(zhuǎn)基因成熟種子為外植體，根據(jù)甘露醇高滲培養(yǎng)基培養(yǎng)法，設(shè)計了0、200、300、400、500、600 mmol/L 6個梯度，對狗尾草ME34品系種子萌發(fā)期的抗旱性進行了鑒定。結(jié)果表明：甘露醇200～300 mmol/L的梯度是鑒定狗尾草ME34品系抗旱性的關(guān)鍵梯度，種子萌發(fā)期發(fā)芽率、鮮重、根長、芽長等各項生長指標(biāo)都受到顯著或極顯著的影響，各生長指標(biāo)對甘露醇干旱脅迫的敏感性依次為鮮重>根長>芽長>發(fā)芽率，參數(shù)間兩兩顯著或極顯著正相關(guān)；鮮重對甘露醇的脅迫的反應(yīng)最敏感和直接，在200 mmol/L梯度下已表現(xiàn)為極顯著降低，可作為該品系種子萌發(fā)期抗旱性鑒定的一級指標(biāo)；根長對甘露醇梯度的反應(yīng)是低濃度促進根的生長，高濃度抑制根的生長，極顯著促進生長的梯度（200 mmol/L）適合用作鑒定依據(jù)，此外，芽長對甘露醇也敏感，200、300 mmol/L分別顯著和極顯著抑制芽長，根長和芽長可作為該品系種子萌發(fā)期抗旱性鑒定的二級指標(biāo)；種子發(fā)芽率對甘露醇的敏感性稍弱，200 mmol/L條件下發(fā)芽率變化不明顯，在300 mmol/L條件下才表現(xiàn)為極顯著下降，可作為該品系種子萌發(fā)期抗旱性鑒定的三級指標(biāo)。實驗為狗尾草作為抗旱模式植物的研究準(zhǔn)備了抗旱性評價指標(biāo)基礎(chǔ)，無論是外源抗旱基因在狗尾草內(nèi)的過表達，還是內(nèi)源抗旱基因的敲除，都需要這些具體評價指標(biāo)的支撐。

關(guān)鍵詞 狗尾草；萌發(fā)期；抗旱性鑒定；甘露醇

中圖分類號 S544.3 文獻標(biāo)識碼 A

Abstract Setaria viridis was a new model plant， which meets all the requirements of model plants. It was an excellent model of C4 photosynthesis， and for studying the tolerance of abiotic stress. Based on the efficient genetic transformation of Setaria viridis， the non-transgenic mature seeds of ME34 strain which was Setaria viridis suitable for transformation were choosen as explants to study drought tolerance. According to Mannitol medium culture method，0， 200， 300， 400， 500， 600 mmol/L 6 gradients were designed. The germination rate，bud length，root length and fresh weight of each gradient were measured and analyzed. The results showed that the 200-300 mmol/L gradient of Mannitol were the key gradients of the Mannitol tolerance for ME34 strain. The germination rate，fresh weight，root length and bud length of seed germination were significant or extremely significant affected. The susceptibility of each growth index to mannitol drought stress was fresh weight>root length>bud length>germination rate. There were significant or extremely significant correlations between the four parameters. The fresh weight was the most sensitive and direct to the mannitol stress and extremely significant reduced at 200 mmol/L gradient，which could be used as the first grade index for drought resistance identification at seed germination stage. The sensitivity of the root length to mannitol stress was also very strong，but as a drought resistance index，only the extremely significant gradient（200 mmol/L）was suitable for identification. In addition，bud length was also sensitive to mannitol，200， 300 mmol/L，was the significant and extremely significant reduced gradient，respectively. Root length and bud length could be used as the second grade index for drought resistance identification at seed germination stage. The seed germination rate was not significant at 200 mmol/L，and it showed a extremely significant decrease under the condition of 300 mmol/L. Compared with other growth indexes，the susceptibility was weak， which could be used as the third grade of drought resistance identification index. The study prepared the identification and evaluation index of drought tolerance at seed germination stage for the pattern plant of Setaria viridis under Mannitol Stress. Whether the overexpression of exogenous drought-resistant genes of Setaria viridis，or the knockout of endogenous anti-drought gene，both need the support of these specific evaluation indicators.

Key words Setaria viridis； germination stage； identification of drought tolerance； mannitol

doi 10.3969/j.issn.1000-2561.2017.11.011

狗尾草（Setaria viridis）起源于熱帶，屬單子葉植物綱禾本科黍亞科，該科植物具有熱帶植物的典型特征。狗尾草與谷子、玉米、高粱、甘蔗、薏苡以及重要能源草類親緣關(guān)系接近，具有C4光合作用系統(tǒng)，在國際上主要作為C4光合作用的優(yōu)秀模型。此外，狗尾草本身耐熱、耐旱、耐鹽，是研究非生物脅迫耐受性很好的模型；狗尾草是新型的轉(zhuǎn)基因模式植物，與雙子葉模式植物擬南芥相比具有生長周期短、植株矮小、易于種植、容易轉(zhuǎn)化、基因組小、呈二倍體、能產(chǎn)生大量自交系種子等優(yōu)點，是優(yōu)良的單子葉模式植物[1-10]，幾乎符合作為一種模式植物的所有要求。

在擬南芥模式植物實際應(yīng)用中科學(xué)家們發(fā)現(xiàn)，由于許多糧食作物為禾本科、單子葉植物，而擬南芥為雙子葉植物，較遠的親緣關(guān)系，導(dǎo)致同家族基因的功能出現(xiàn)嚴(yán)重分化，致基因功能完全相反。因此，擬南芥中基因功能的研究，離其在最終農(nóng)作物中的應(yīng)用還很遠[11-12]。于是，科學(xué)家們又發(fā)展了水稻作為禾本科的模式生物，由于其基因組小，易于轉(zhuǎn)化等諸多優(yōu)點，以水稻為模式植物開展的研究也越來越多，尤其是中國在水稻的功能基因組學(xué)方面的研究已經(jīng)走在世界的前列。但是，水稻作為模式植物有先天的缺點，它對生長對環(huán)境的要求比較嚴(yán)格，在北方每年只能種一季。

狗尾草在研究C4植物的光合作用機理和非生物脅迫耐受性方面，具有得天獨厚的優(yōu)勢，與許多禾本科農(nóng)作物如谷子、玉米等具有較近的親緣關(guān)系，它們的基因功能更具有較高的保守性。Huang等[13]以狗尾草為模式植物，通過NMU誘變，篩選到了一個花絮稀疏的突變體，通過BSA重測序方法，鑒定到了突變基因，而且玉米中的同源基因突變后也表現(xiàn)出與狗尾草突變體相同的表型，這進一步證明了以狗尾草為模式植物進行玉米基因組功能研究的優(yōu)越性。

干旱是一種世界性的災(zāi)害，因干旱而引起的糧食減產(chǎn)比其它所有逆境因素所導(dǎo)致的減產(chǎn)總和還要多。 據(jù)統(tǒng)計，世界上干旱及半干旱地區(qū)約占地球陸地面積的34.9%，中國干旱、半干旱地區(qū)占國土面積的47%，占總耕地面積的51%，干旱缺水是中國農(nóng)業(yè)面臨的最嚴(yán)重的問題之一[14-17]。因此，作為模式植物的狗尾草，在高效遺傳轉(zhuǎn)化技術(shù)的基礎(chǔ)上開展的抗旱基因功能研究在中國必將具有更重要的意義，但國內(nèi)外將狗尾草作為抗旱轉(zhuǎn)基因育種模式植物的研究還很少。本研究在狗尾草高效遺傳轉(zhuǎn)化技術(shù)的基礎(chǔ)上，以適宜轉(zhuǎn)化的狗尾草ME34品系非轉(zhuǎn)基因種子為外植體，根據(jù)甘露醇高滲培養(yǎng)基培養(yǎng)法，設(shè)計了0、200、300、400、500、600 mmol/L 6個梯度，對狗尾草ME34品系種子萌發(fā)期的抗旱性進行了鑒定，為狗尾草作為抗旱轉(zhuǎn)基因育種模式植物的抗旱性評價指標(biāo)的鑒定打下基礎(chǔ)。

1 材料與方法

1.1 材料

以適宜遺傳轉(zhuǎn)化的狗尾草ME34品系為植物材料，材料由美國Donald Danforth Plant Science Center的Thomas P. Brutnell 實驗室提供，以ME34品系經(jīng)休眠處理的成熟種子為外植體材料，研究種子萌發(fā)期的抗旱性。

1.2 方法

1.2.1 相關(guān)組培技術(shù) 甘露醇高滲培養(yǎng)基配制：1/2MS培養(yǎng)基，10 g蔗糖，pH5.8，Gelzan植物凝膠4 g，在基礎(chǔ)培養(yǎng)基內(nèi)分別添加甘露醇，添加濃度梯度為0、200、300、400、500、600 mmol/L，于120 ℃高壓滅菌20 min后倒皿備用。

種子的無菌處理和接種：選取經(jīng)休眠處理的狗尾草干燥成熟種子，去除種皮后將約1 g種子裝入14 mL離心管，加入10 mL有效濃度約1%～3%的次氯酸鈉溶液，再加50 μL吐溫20溶液，滅菌消毒處理約10 min；用無菌水清洗多次后，將種子接種于含不同濃度梯度的甘露醇的高滲培養(yǎng)基上，種子的芽朝上，每皿培養(yǎng)基約接20粒種子，每個梯度接3皿（3個重復(fù)），于24 ℃光/暗培養(yǎng)（12 000 lx，16 h/8 h，50%～60%濕度），7 d以后觀察記錄各生長指標(biāo)。

1.2.2 調(diào)查項目及分析方法

發(fā)芽率=發(fā)芽種子數(shù)/接種的種子數(shù)×100%。

芽長平均值=同一脅迫梯度下發(fā)芽種子苗地上芽長的均值。

根長平均值=同一脅迫梯度下發(fā)芽種子苗地下主根長的均值。

相對發(fā)芽率=脅迫處理的發(fā)芽率/對照處理的發(fā)芽率×100%。

芽生長抑制率=（對照芽長的均值-脅迫處理芽長的均值）/對照芽長的均值×100%。

根生長抑制率=（對照根長的均值-脅迫處理根長的均值）/對照根長的均值×100%。

鮮重抑制率=（對照鮮重的均值-脅迫處理鮮重的均值）/對照鮮重的均值×100%。

1.3 數(shù)據(jù)分析

數(shù)據(jù)、圖標(biāo)處理采用Excel 2010軟件完成， 利用Sigma Plot 10.0軟件對數(shù)據(jù)進行統(tǒng)計分析、差異顯著性檢驗和相關(guān)性分析。

2 結(jié)果與分析

2.1 甘露醇對狗尾草種子萌發(fā)的整體影響

從圖1可知，與無甘露醇的對照相比，從種子萌發(fā)第7天的整個生長情況看，200 mmol/L甘露醇已經(jīng)開始抑制狗尾草生長量（B），300 mmol/L表現(xiàn)出非常明顯的抑制作用（C），B、C、D、E、F中狗尾草中的生長受抑制情況依次加重，400 mmol/L以上非轉(zhuǎn)基因的ME34種子已很難生長（D）。根據(jù)測算，表1具體列出了各生長指標(biāo)數(shù)據(jù)。

2.2 甘露醇對狗尾草種子發(fā)芽率的影響

圖2是發(fā)芽率隨甘露醇濃度梯度的變化圖。從圖中可知，發(fā)芽率隨著甘露醇濃度的提高總體呈下降趨勢，在0～300 mmol/L的范圍內(nèi)發(fā)芽率下降較快，300～400 mmol/L的范圍內(nèi)變化平緩，400～600 mmol/L范圍內(nèi)再次急速下降。此外，從各梯度發(fā)芽率方差分析可知，0 mmol/L與200 mmol/L的發(fā)芽率差異不顯著（p=0.095>0.05），與300～600 mmol/L差異極顯著（p=0.000 54，0.000 79，0.000 24， 1.17E-05<0.01）；200 mmol/L與300～600 mmol/L差異極顯著（p=1.41E-05，8.03E-05，7.37E-05，8.42E-08<0.01）；300 mmol/L與400 mmol/L差異不顯著（p=0.4625>0.05），與500～600 mmol/L差異極顯著（p=0.003 3，2.47E-06<0.01）；400 mmol/L與500～600 mmol/L差異極顯著（p=0.0032，6.1521E-06<0.01）；500 mmol/L與600 mmol/L差異極顯著（p=0.0002<0.01）。從各梯度變化圖和方差分析的整體情況看，200 mmol/L對發(fā)芽率影響不顯著，300 mmol/L極顯著影響發(fā)芽率，如果需要通過發(fā)芽率考察ME34品系轉(zhuǎn)基因后的抗旱能力，需要用到300 mmol/L的梯度。

2.3 甘露醇對狗尾草種子萌發(fā)期鮮重的影響

圖3是鮮重隨甘露醇濃度梯度的變化圖。從圖3可知，鮮重隨著甘露醇濃度的提高總體呈下降趨勢，在0～200 mmol/L的范圍內(nèi)鮮重下降明顯，200～600 mmol/L范圍內(nèi)鮮重下降平緩。此外，從各梯度鮮重方差分析可知，0 mmol/L與200、300、400、500、600 mmol/L的鮮重差異極顯著（p=0.006 2， 0.002 3，0.001 8，0.001 4，0.000 8<0.01）；200 mmol/L與300、400 mmol/L差異不顯著（p=0.124 6， 0.051 8>0.05），與500 mmol/L差異顯著（p=0.021 1 <0.05），與600 mmol/L差異極顯著（p=0.003 6<0.01）；300 mmol/L與400 mmol/L差異不顯著（p=0.226 7>0.05），與500 mmol/L差異顯著（p=0.042 2 <0.05），與600 mmol/L差異極顯著（p=0.001 3<0.01）；400 mmol/L與500 mmol/L差異不顯著（p=0.202 5>0.05），與600 mmol/L差異極顯著（p=0.005 2 <0.01）；500 mmol/L與600 mmol/L差異顯著（p=0.031 3<0.05）。從各梯度變化圖和方差分析的整體情況看，200 mmol/L極顯著影響鮮重，如果需要通過鮮重考察ME34品系轉(zhuǎn)基因后的抗旱能力，需要用到200 mmol/L的梯度。

2.4 甘露醇對狗尾草種子胚芽生長的影響

圖4是芽長隨甘露醇濃度梯度的變化圖。從圖中可知，芽長隨著甘露醇濃度的提高總體呈縮短趨勢，在0～300 mmol/L的范圍內(nèi)芽長下降明顯，300～600 mmol/L范圍內(nèi)芽長下降平緩。此外，從各梯度芽長方差分析可知：0 mmol/L與200 mmol/L的芽長差異顯著（p=0.017 9<0.05），與300～600 mmol/L差異極顯著（p=0.000 8，0.000 4，0.000 3， 0.000 2<0.01）；200 mmol/L與300～600 mmol/L差異極顯著（p=0.002 8，0.001 0，0.000 4，0.000 2<0.01）；300 mmol/L與400 mmol/L差異顯著（p=0.026 1<0.05），與500～600 mmol/L差異極顯著（p=0.002 7，0.000 5<0.01）；400 mmol/L與500 mmol/L差異顯著（p=0.046 1<0.05），與600 mmol/L差異極顯著（p=0.002 5<0.01）；500 mmol/L與600 mmol/L差異顯著（p=0.012 9<0.05）。從各梯度變化圖和方差分析的整體情況看，200 mmol/L顯著影響芽長，300 mmol/L極顯著影響芽長。如果需要通過芽長考察ME34品系轉(zhuǎn)基因后的抗旱能力，需要用到200 mmol/L和300 mmol/L的梯度。

2.5 甘露醇對狗尾草種子胚根生長的影響

圖5是根長隨甘露醇濃度梯度的變化圖。從圖中可知，根長的變化與發(fā)芽率、鮮重、芽長的變化有所不同，隨著甘露醇濃度的提高總的趨勢是先升高后降低，在0～300 mmol/L的范圍內(nèi)根長有先顯著升高后下降的趨勢，300～600 mmol/L范圍內(nèi)根長一直下降。此外，從各梯度根長方差分析可知，0 mmol/L與200 mmol/L的根長差異顯著（p=0.012 7<0.05），與300 mmol/L差異不顯著（p=0.383 8>0.05），與400 mmol/L差異顯著（p=0.032 6<0.05），與500～600 mmol/L差異極顯著（p=0.004 5，0.000 3<0.01）；200 mmol/L與300～600 mmol/L差異極顯著（p=0.006 1， 0.001 1，0.000 4，7.19E-05<0.01）；300 mmol/L與400 mmol/L差異不顯著（p=0.1140>0.05），與500～600 mmol/L差異極顯著（p=0.009 8，0.000 4<0.01）；400 mmol/L與500 mmol/L差異顯著（p=0.011 9<0.05），與600 mmol/L差異極顯著（p=5.85E-05<0.01）；500 mmol/L與600 mmol/L差異極顯著（p=0.000 4<0.01）。從各梯度變化圖和方差分析的整體情況看，根長隨著甘露醇脅迫濃度的增加有先增長后縮短的趨勢，在200 mmol/L梯度下根長極顯著增長。如果考察ME34品系轉(zhuǎn)基因后的抗旱能力，200 mmol/L梯度下根長的增長量可以用來鑒定抗旱性的強弱；由于0 mmol/L與300 mmol/L差異不顯著，同時，300 mmol/L的梯度下非轉(zhuǎn)基因ME34品系的總體生長量已經(jīng)很小，所以300 mmol/L以上的梯度不適合用來從根長的角度進行轉(zhuǎn)基因狗尾草抗旱鑒定。

2.6 甘露醇脅迫下狗尾草種子萌發(fā)期各生長指標(biāo)的相關(guān)性分析

甘露醇脅迫下狗尾草ME34品系種子萌發(fā)期各生長指標(biāo)的相關(guān)性見表2。其中，發(fā)芽率與根長、芽長、鮮重呈極顯著正相關(guān)；根長與芽長呈極顯著正相關(guān)，根長與鮮重呈顯著正相關(guān)；芽長與鮮重呈極顯著正相關(guān)。

3 討論

從狗尾草的起源、分布和生存環(huán)境可知，狗尾草具有豐富的抗性基因資源，狗尾草是谷子（Setaria italica）的野緣近親，有研究表明谷子的很多栽培品系明顯抗旱，1 g的干物質(zhì)產(chǎn)量，玉米需要470 g水，小麥需要510 g水，谷子只需要257 g水[18-19]。狗尾草是比擬南芥更貼近主要農(nóng)作物的模式植物，目前針對狗尾草的研究主要還在C4機理方面，針對其抗旱性開展的研究還很少，有很大的研究空間[11-13]。狗尾草的抗旱機理涉及代謝、蛋白水解和信號傳遞等各種機制[20-21]。在狗尾草高效遺傳轉(zhuǎn)化技術(shù)的基礎(chǔ)上，以狗尾草為模式植物，通過轉(zhuǎn)基因技術(shù)的過表達或基因編輯等技術(shù)對狗尾草的抗旱機理進行研究，必將取得豐碩的研究成果，并指導(dǎo)作物育種工作。

植物的抗旱性是一個極其復(fù)雜的數(shù)量性狀，針對模式植物狗尾草尚無一個簡單有效、經(jīng)濟方便、實用可靠、便于實驗室參考使用的抗旱性鑒定指標(biāo)體系。本研究在參考玉米、大豆、水稻等大宗作物抗旱性鑒定方法的基礎(chǔ)上，選擇甘露醇高滲培養(yǎng)基培養(yǎng)法，在組織培養(yǎng)條件下，對狗尾草ME34品系種子萌發(fā)期的抗旱性進行了鑒定，為模式植物狗尾草抗旱轉(zhuǎn)基因研究提供了抗旱性鑒定指標(biāo)體系。

實驗結(jié)果表明，甘露醇200～300 mmol/L梯度是鑒定狗尾草ME34品系抗旱性的關(guān)鍵梯度，種子萌發(fā)期發(fā)芽率、鮮重、根長、芽長等各項生長指標(biāo)都受到顯著或極顯著的影響。其中鮮重對甘露醇的脅迫最敏感，在200 mmol/L梯度下已表現(xiàn)出極顯著降低趨勢，可作為該品系種子萌發(fā)期抗旱性鑒定的一級指標(biāo)；根長對甘露醇脅迫的敏感性也很強，在200 mmol/L梯度下根長極顯著增加，但根長對甘露醇梯度的反應(yīng)是低濃度促進根的生長，高濃度抑制根的生長，在高濃度抑制的條件下種子已難以萌發(fā)生長，所以作為抗旱指標(biāo)，極顯著促進生長的梯度（200 mmol/L）適合用作鑒定依據(jù)；此外，芽長對甘露醇也敏感，200、300 mmol/L分別顯著和極顯著抑制芽長，根長和芽長可作為該品系種子萌發(fā)期抗旱性鑒定的二級指標(biāo)；種子發(fā)芽率在200 mmol/L條件下變化不明顯，在300 mmol/L條件下才表現(xiàn)出極顯著下降趨勢，相對其它的生長指標(biāo)，其敏感性稍弱，可作為該品系種子萌發(fā)期抗旱性鑒定的三級指標(biāo)。

參考文獻

[1] Brutnell T P， Wang L， Swartwood K， et al. Setaria viridis： A model for C4 pHotosynthesis[J]. The Plant Cell， 2010， 22（8）： 2 537-2 544.

[2] Van E J， Swarwood K. Setaria viridis， agrobaterium protocols Volume[J]. Methods in Molecular Biology， 2015， 1 223： 5 767.

[3] Li P H， Brutnell T P. Setaria viridis and Setaria italica， model genetic systems for the Panicoid grasses[J]. Journal of Experimental Botany， 2011， 62（9）： 3 031-3 037.

[4] Jiang H， Barbier H， Brutnell T. Methods for performing crosses in Setaria viridis， a new model system for the grasses[J]. Journal of Visualized Experiments， 2013， 80： 50 527.

[5] 李秀蘭， 陳瑞陽. 谷子和狗尾草核型分析[J]. 武漢植物學(xué)研究， 1985， 3（4）： 409-412.

[6] 李 瑋， 姜翠茹， 陳龍俊. C4模式植物的研究進展[J]. 現(xiàn)代農(nóng)業(yè)科技， 2010（23）： 22.

[7] 龔秀秀， 王 凱， 董 雯， 等. C4光合作用調(diào)節(jié)機制[J]. 生物技術(shù)通報， 2013（3）： 19-23.

[8] 智 慧， 刁現(xiàn)民， 呂 芃， 等. 人工鹽脅迫法鑒定谷子及狗尾草物種耐鹽基因型[J]. 河北農(nóng)業(yè)科學(xué)， 2004， 8（4）： 15-18.

[9] 殷立娟， 李美榮. 中國C4植物的地理分布與生態(tài)學(xué)研究Ⅰ. 中國C4植物及其與氣候環(huán)境的關(guān)系[J]. 生態(tài)學(xué)報， 1997， 17（4）： 350-354.

[10] 孫 偉， 王德利， 王 立， 等.狗尾草蒸騰特性與水分利用效率對模擬光輻射增強和CO2濃度升高的響應(yīng)[J]. 植物生態(tài)學(xué)報， 2003， 27（4）： 448-453.

[11] Brutnell T P. Model grasses hold key to crop improvement[J]. Nature Plants， 2015， 1（5）： 15 062.

[12] Brutnell T P， Bennetzen J L， Vogel J P. Brachypodium distachyon and Setaria viridis： model genetic systems for the grasses[J]. Annual Review of PlantBiology， 2015， 66（1）： 465.

[13] Huang P， Jiang H， Zhu C， et al. Sparse panicle1 is required forinflorescence development in Setaria viridis and maize[J]. Nature Plants， 2017， 3（5）： 17 054.

[14] 湯章城. 植物對水分脅迫的反應(yīng)和適應(yīng)性——Ⅰ.抗逆性的一般概念和植物的抗?jié)承訹J]. 植物生理學(xué)報， 1983（3）： 26-31.

[15] 西北農(nóng)林大學(xué). 旱農(nóng)學(xué)[M]. 北京： 農(nóng)業(yè)出版社， 1991： 1-6.

[16] 康紹忠. 新農(nóng)業(yè)科技革命與21世紀(jì)我國節(jié)水農(nóng)業(yè)的發(fā)展[J]. 干旱地區(qū)農(nóng)業(yè)研究， 1998， 16（1）： 11-17.

[17] 李原園， 李英能， 蘇人瓊. 中國農(nóng)業(yè)水危機及其對策[R] . 北京： 中國國家科學(xué)技術(shù)委員會農(nóng)村科技司， 1997： 52-54.

[18] Diao X M. Adance in foxtial millet biotechnology and its future directions[J]. Journal of Hebei Agricultural Sciences， 2005（9）：61.

[19] Diao X. Foxtail millet production in China and its future development tendency[M]//Chai Y， Wan F S. The industrial development of China special crops. Beijing： Chinese Agricultural Science and Technology Press， 2007.

[20] Zhang J， Liu T， Fu J， et al. Construction and application of EST library from Setaria italica in response to dehydration stress[J]. Genomics， 2007， 90（1）： 121.

[21] Lata C， Sahu P P， Prasad M. Comparative transcriptome analysis of differentially expressed genes in foxtail millet （Setaria italica L.）during dehydration stress[J]. Biochemical and Biophysical Research Communications， 2010， 393（4）： 720-727.