2020年8月31日，山西農(nóng)業(yè)大學(xué)雜糧分子育種團(tuán)隊(duì)與中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所等單位合作研究論文“A mini-foxtail millet with an Arabidopsis-like life cycle as a C4 model system”在線(xiàn)發(fā)表于國(guó)際期刊Nature Plants（IF：13.256），報(bào)道了C4禾谷類(lèi)研究的理想模式植物“小米”的研究成果。山西農(nóng)業(yè)大學(xué)的楊致榮教授、中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所的張皓珊博士和山西農(nóng)業(yè)大學(xué)的李旭凱副教授為論文的共同第一作者；山西農(nóng)業(yè)大學(xué)的王興春教授、韓淵懷教授和中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所的隋毅副研究員為論文的共同通訊作者。美國(guó)普渡大學(xué)和英國(guó)諾丁漢大學(xué)等單位參與了部分研究。百邁客有幸參與了其中基因組、轉(zhuǎn)錄組及數(shù)據(jù)庫(kù)構(gòu)建工作。目前，數(shù)據(jù)庫(kù)由百邁客負(fù)責(zé)維護(hù)和運(yùn)行，暫時(shí)托管在公司服務(wù)器（http://foxtail-millet.biocloud.net/home）。
研究背景
眾所周知，與C3植物相比，C4植物具有較高的光合效率和水肥利用率。因此，若能將C4光合途徑整合到水稻等C3作物，將有望大幅提高水稻產(chǎn)量，解決由于人口增長(zhǎng)帶來(lái)的溫飽問(wèn)題。然而，目前常用的模式植物擬南芥和水稻都是C3植物，作為C4作物的模式植物時(shí)有很大的局限性，無(wú)法解決諸如C4光合代謝以及許多黍亞科特殊的基礎(chǔ)問(wèn)題。
谷子（Setaria italica）起源于中國(guó)，是我國(guó)幾千年來(lái)的主栽作物和中華民族的哺育作物。此外，具有抗旱、耐瘠薄和高光效等突出優(yōu)勢(shì)，恰恰彌補(bǔ)了擬南芥和水稻作為模式植物的不足，是極具發(fā)展?jié)摿Φ暮坦阮?lèi)模式植物。但谷子生育期較長(zhǎng)、株高較高、遺傳轉(zhuǎn)化困難，極大地限制了其作為模式植物在功能基因組學(xué)研究中的應(yīng)用。本次研究篩選到一個(gè)超早熟谷子突變體“小米”，其生命周期短，植株小，與擬南芥相似?！靶∶住蹦Ｊ街参矬w系的建立將極大促進(jìn)C4高光效、氮素高效吸收利用機(jī)制、抗旱機(jī)制、遺傳馴化和優(yōu)異品質(zhì)形成的分子基礎(chǔ)等研究，使“小米”成為C4植物功能研究的理想模式系統(tǒng)。
研究材料與方法
每株系約50株T1轉(zhuǎn)基因幼苗進(jìn)行重測(cè)序，每株系約12 Gb的數(shù)據(jù)（~28×）。晉谷21、“小米”、“小米”2進(jìn)行重測(cè)序比較分析。
11個(gè)不同的組織：3日齡種子（種子）、2周齡整株苗（苗）、根、莖、2周齡苗的第一片全展葉（葉1），30日齡植株頂?shù)诙~（葉2）、旗葉（葉3）、第四葉（葉4）、未成熟的圓錐花序（圓錐花序1）、授粉期圓錐花序（圓錐花序2）、灌漿期圓錐花序（圓錐花序3）；葉2有5個(gè)生物學(xué)重復(fù)，其他有3個(gè)生物學(xué)重復(fù)。
主要研究?jī)?nèi)容
1.?“小米”的創(chuàng)造與表型表征
2.?“小米”的特性由PHYC基因突變引起
3.?“小米”基因組組裝注釋
4.?谷子品種間基因組序列比較
5.?“小米”動(dòng)態(tài)基因表達(dá)圖譜的構(gòu)建
6.?農(nóng)桿菌介導(dǎo)高效遺傳轉(zhuǎn)化體系的建立
研究結(jié)果
1.?“小米”的創(chuàng)造與表型表征
利用EMS誘變技術(shù)對(duì)山西省名優(yōu)谷子晉谷21進(jìn)行了誘變，從中篩選到一個(gè)超早熟突變體“小米”。
“小米”主要有以下特征：（1）表現(xiàn)出極早的開(kāi)花表型：抽穗期39天（播種后天數(shù)，DAS）左右，約70天內(nèi)完成了生命周期。而野生型（晉谷21，WT）的平均抽穗期約為82 天（DAS），成熟期為130 DAS（圖1b、c）；（2）“小米”比野生型矮（圖1b），但結(jié)實(shí)率高12.83%，且種子大小沒(méi)有顯著差異（圖1d）；（3）“小米”在短日（SD，10h光照/14h暗）條件下比長(zhǎng)日（LD，16h光/8h暗）條件下的抽穗晚了約一個(gè)月，表明“小米”的早期抽穗依賴(lài)于LD條件（圖1e）。通過(guò)延長(zhǎng)日長(zhǎng)和優(yōu)化其他條件，可以將“小米”的生命周期縮短到65天，株高約為29厘米（圖1f），每年可以在成長(zhǎng)室里完成5到6代。（4）僅需較小的種植面積，類(lèi)似于種植等量的擬南芥植株。
圖1 “小米”的表型特征
2.?“小米”的特性由PHYC基因突變引起
為了找出導(dǎo)致“小米”早熟的突變，研究者將其與抽穗期為75 DAS的G1雜交。所有9株F1植株均表現(xiàn)出類(lèi)似G1的晚抽穗表型，由268株組成的F2群體表現(xiàn)出G1型晚抽穗期與“小米”型早抽穗期的分離比為3:1（?χ2=0.318<χ2?0.05（1）=3.841），說(shuō)明“小米”早熟抽穗期是由單基因座隱性突變引起的。利用106個(gè)早花F2個(gè)體，將該位點(diǎn)定位到9號(hào)染色體上的212kb區(qū)域，根據(jù)“小米”參考基因組的注釋?zhuān)搮^(qū)域含有27個(gè)基因（圖2a）。
將“小米”基因組序列與經(jīng)基因組重測(cè)序生成的晉谷21基因組序列進(jìn)行比較，發(fā)現(xiàn)在定位區(qū)域只存在一個(gè)突變——“小米”基因Si9g09200的編碼區(qū)從“G”到“T”的突變，該基因編碼一種假定的PHYC蛋白（圖2a）。在所檢測(cè)的77個(gè)早期抽穗F2個(gè)體中，均為T(mén)/T純合突變體。相比之下，在256個(gè)晚抽穗個(gè)體中，92個(gè)是G/G WT純合子，另外164個(gè)是G/T雜合子，反映了基因型和表型之間的完美關(guān)聯(lián)。這種突變產(chǎn)生了一個(gè)終止密碼子，形成一個(gè)截短的蛋白質(zhì)，占該基因約71%的轉(zhuǎn)錄本（轉(zhuǎn)錄本1）（圖2a-d）。這一突變還導(dǎo)致了113 bp的移碼缺失的剪接，形成了占該基因約29%轉(zhuǎn)錄本（轉(zhuǎn)錄本2）的截短蛋白（圖2a-d）。根據(jù)預(yù)測(cè)，被截?cái)嗟牡鞍踪|(zhì)缺少第二個(gè)PAS結(jié)構(gòu)域的三分之二（肽1）或三分之一（肽2），以及整個(gè)HK和HD域（圖2d）。
圖2 “小米”的分子特征
研究者對(duì)另一個(gè)來(lái)源于晉谷21的突變株——“小米”2號(hào)進(jìn)行了測(cè)序?！靶∶住?號(hào)突變體表現(xiàn)出與“小米”相似的早期抽穗表型。PHYC基因座序列比較顯示“小米”2號(hào)PHYC第一外顯子發(fā)生單點(diǎn)突變（T674A），導(dǎo)致從保守的亮氨酸變?yōu)榻M蛋白（圖3）。這種SNP與M2雜合子突變體的82株早抽穗（A/A基因型）和84株晚抽穗（49 T/A基因型和35 T/T基因型）M3植株的表型分離密切相關(guān)。綜上所述，這些觀(guān)察結(jié)果證實(shí)了早期抽穗表型是由PHYC位點(diǎn)突變引起的。
圖3 “小米”2號(hào)突變體的表型和分子生物學(xué)特性
為了了解PHYC位點(diǎn)突變對(duì)LD條件下開(kāi)花時(shí)間的影響，研究者利用30-DAS（“小米”抽穗前10天）植株的第二片葉片進(jìn)行了RNA-seq，并以30-DAS采集的WT葉片為對(duì)照。發(fā)現(xiàn)與擬南芥振蕩器基因、偽應(yīng)答調(diào)節(jié)蛋白（PRRs）、PHYTOCLOCK 1?(PCL1) 和GIGANTEA?(GI)同源的幾個(gè)基因在“小米”中的表達(dá)受到顯著影響（圖4）。
與Ghd7同源的下游光周期基因的表達(dá)水平下降了約95倍（長(zhǎng)日照條件下,增強(qiáng)Ghd7的表達(dá)可延遲抽穗期），而與HEADING DATE 1 (Ehd1), HEADING DATE?3a (Hd3a),APETALA1?(AP1)/FRUITFULL?(FUL)-like?MADS box基因同源的開(kāi)花基因在“小米”中的表達(dá)水平顯著增加。Ehd1是水稻光周期開(kāi)花途徑中的一個(gè)關(guān)鍵轉(zhuǎn)錄調(diào)控因子，可以激活葉片中成花素Hd3a的轉(zhuǎn)錄。在葉片中產(chǎn)生的Hd3a蛋白被轉(zhuǎn)運(yùn)到莖尖分生組織，在那里誘導(dǎo)AP1/FUL-like?MADS box基因的表達(dá)。這些觀(guān)察結(jié)果表明，“小米”早熟的表型是由于光周期通路的中斷引起的。
圖4?光周期途徑RNA-seq結(jié)果的qRT-PCR驗(yàn)證
（注：這些蛋白質(zhì)編碼基因命名：Si（代表“小米”），其次是染色體數(shù)目和染色體上的基因數(shù)，例如，1號(hào)染色體上的第一個(gè)和第二個(gè)基因分別命名為Si1G00010和Si1G00020）
3.?“小米”基因組組裝注釋
為了便于“小米”作為模型植物的使用，使用PacBio平臺(tái)組裝出429.94 Mb基因組 v1.0（k-mer預(yù)估438.26 Mb），contig N50高達(dá)18.8?Mb，Scaffold N50為42.41 Mb，并且僅包含48個(gè)Gap。其中399.4 Mb通過(guò)Hi-C錨定在9條染色體上（掛載率~93%；圖5）。根據(jù)Illumina DNA比對(duì)，組裝的錯(cuò)誤率約為0.001%（每100kb有一個(gè)錯(cuò)誤）。BUSCO評(píng)估結(jié)果為97.78%。這些結(jié)果表明“小米”基因組v1.0可以作為研究界的金標(biāo)準(zhǔn)參考。
圖5 “小米”基因組Hi-C熱圖
“小米”基因組序列中有237.28 Mb（55.19%）為重復(fù)序列?；趶念^預(yù)測(cè)、同源基因預(yù)測(cè)并結(jié)合Iso-seq和RNA-seq，對(duì)34436個(gè)蛋白質(zhì)編碼基因進(jìn)行了注釋?zhuān)渲?2743個(gè)（95.08%）位于9個(gè)假染色體上。在GO、KEGG、KOG、TrEMBL、NR數(shù)據(jù)庫(kù)中檢索這些基因，并與擬南芥和水稻的注釋進(jìn)行比較，檢索出具有已知功能的同源基因，共注釋33789個(gè)基因（98.12%）。此外還注釋了919個(gè)rRNA基因、3516個(gè)tRNA基因、2631個(gè)假基因、340個(gè)miRNA前體、28260個(gè)lncRNA前體和1318個(gè)circRNA前體（表1）。
所有的基因組和轉(zhuǎn)錄組數(shù)據(jù)可在數(shù)據(jù)庫(kù)(http://sky.sxau.edu.cn/MDSi.htm)公開(kāi)訪(fǎng)問(wèn)。在這個(gè)數(shù)據(jù)庫(kù)中，研究人員可以通過(guò)染色體坐標(biāo)、基因或轉(zhuǎn)錄符號(hào)，或者通過(guò)BLAST搜索“小米”基因組、CDS或肽序列來(lái)尋找基因組位置。
表1 “小米”基因組組裝與注釋統(tǒng)計(jì)?
4.?谷子品種間基因組序列比較
與 豫谷1、張谷和 TT8的3個(gè)基因組序列相比，“小米”基因組在基因組覆蓋率、contig N50、contig和gap數(shù)方面表現(xiàn)出*高的質(zhì)量?！靶∶住毙蛄兄?14.58 Mb（96.44%）對(duì)應(yīng)于豫谷1號(hào) 383.52 Mb（95.67%），共有1577935個(gè)SNPs。相應(yīng)區(qū)域的大小差異主要由259731個(gè)InDel，2804個(gè)（總計(jì)15.32 Mb）存在變異（>1000bp）和2722個(gè)（總計(jì)17.38 Mb）缺失變異（>1000bp）引起（圖6）。
在“小米”基因組34436個(gè)蛋白編碼基因中，有32112個(gè)基因（93.25%）與豫谷1號(hào)基因組（v2.2）共有。2280個(gè)“小米”預(yù)測(cè)基因在張谷基因組中未被發(fā)現(xiàn)。在豫谷1號(hào)和張谷中均未發(fā)現(xiàn)的1030個(gè)基因，為“小米”特有基因。
圖6?“小米”與豫谷1號(hào)基因組比較圈圖
5.?“小米”動(dòng)態(tài)基因表達(dá)圖譜的構(gòu)建
為了建立一個(gè)參考的基因表達(dá)圖譜，研究者測(cè)量了“小米”不同發(fā)育階段的11種不同組織的轉(zhuǎn)錄水平?？偣伯a(chǎn)生并分析了1054.51 M reads（每個(gè)樣本約30 M reads）。從頂部第二葉（葉2）到授粉期圓錐花序（圓錐花序2）檢測(cè)到的基因表達(dá)比例為74.26%~82.95%。11個(gè)“小米”組織共表達(dá)22202個(gè)基因。這些基因中，85個(gè)（0.25%）在所有被測(cè)組織中都有組成性表達(dá)，包括一個(gè)轉(zhuǎn)錄起始因子（Si3G07600）和兩個(gè)泛素結(jié)合酶編碼基因（Si1G37980和Si2G05250）。此外，還鑒定了1218個(gè)器官/組織特異基因和1226個(gè)器官/組織優(yōu)先表達(dá)基因。
圖7?Si9g04830基因表達(dá)模式圖
6.?農(nóng)桿菌介導(dǎo)高效遺傳轉(zhuǎn)化體系的建立
研究者對(duì)各種因素進(jìn)行了測(cè)試，以開(kāi)發(fā)出適合“小米”的農(nóng)桿菌介導(dǎo)的轉(zhuǎn)化方案。并挑戰(zhàn)使用成熟種子作為愈傷組織誘導(dǎo)的起始材料，以避免在使用新鮮組織（如幼花序或未成熟胚）誘導(dǎo)愈傷組織時(shí)的昂貴需求。在改良的愈傷組織誘導(dǎo)培養(yǎng)基（CIM）上進(jìn)行三輪繼代培養(yǎng)后，獲得了致密的胚性愈傷組織（圖8a）。通過(guò)使用綠色熒光蛋白（GFP）報(bào)告基因來(lái)監(jiān)測(cè)農(nóng)桿菌感染效率，如圖8b、c所示，以及篩選轉(zhuǎn)基因愈傷組織的有效性?！靶∶住钡脑偕芰υ贑IM培養(yǎng)基上保持良好（圖8f）。表達(dá)GFP的轉(zhuǎn)基因植株在生根培養(yǎng)基上容易誘導(dǎo)生根，移栽到土壤后成活良好（圖8g-i）。
研究者比較了兩個(gè)常用的選擇標(biāo)記NPTII（新霉素磷酸轉(zhuǎn)移酶II）和HPT（潮霉素磷酸轉(zhuǎn)移酶），NPTII的轉(zhuǎn)化效率為8.05%-38.75%，平均轉(zhuǎn)化效率為23.28%；HPT的轉(zhuǎn)化效率為3.08%-16.67%，平均為8.72%。使用分別擴(kuò)增GFP基因、UBI啟動(dòng)子、HPT或NPTII選擇標(biāo)記基因的引物，通過(guò)PCR來(lái)確認(rèn)轉(zhuǎn)基因原植物（T0）中是否存在轉(zhuǎn)基因（圖8j）。觀(guān)察到綠色熒光蛋白在干種子和發(fā)芽種子中的表達(dá)，表明轉(zhuǎn)基因已傳遞給后代（圖8k，l）。
圖8?農(nóng)桿菌介導(dǎo)的“小米”轉(zhuǎn)化研究
用HPT或NPTII標(biāo)記的13個(gè)獨(dú)立轉(zhuǎn)基因系的基因組測(cè)序證實(shí)了轉(zhuǎn)基因插入“小米”基因組的正確性，并通過(guò)PCR進(jìn)一步確定了3個(gè)檢測(cè)系的T-DNA插入位點(diǎn)。在盆栽中培育了代表8個(gè)轉(zhuǎn)基因事件的T1植株，觀(guān)察到與未轉(zhuǎn)化的“小米”植株沒(méi)有明顯的表型差異。因此，研究者建立了一個(gè)有效的方案，允許在2-3個(gè)月內(nèi)（從農(nóng)桿菌感染算起）或4-5個(gè)月內(nèi)（從成熟種子開(kāi)始愈傷組織算起），完成生產(chǎn)準(zhǔn)備移植到土壤中的轉(zhuǎn)基因植物（圖9）。
圖9?農(nóng)桿菌介導(dǎo)“小米”轉(zhuǎn)化工藝流程圖
小結(jié)
本次研究利用EMS誘變技術(shù)對(duì)山西省名優(yōu)谷子晉谷21進(jìn)行了誘變，從中篩選到一個(gè)超早熟突變體“小米”。該突變體生育期僅僅兩個(gè)月左右，株高僅30 cm左右?！靶∶住钡纳诤椭旮吲c模式植物擬南芥相當(dāng)，從而解決了谷子作為C4禾谷類(lèi)模式植物無(wú)法在室內(nèi)大規(guī)模培養(yǎng)的難題。在此基礎(chǔ)上，研究團(tuán)隊(duì)組裝了高質(zhì)量參考基因組，構(gòu)建了全生育期基因表達(dá)圖譜和谷子多組學(xué)數(shù)據(jù)庫(kù)（http://sky.sxau.edu.cn/MDSi.htm），極大地方便了“小米”功能基因組學(xué)研究。經(jīng)過(guò)大量摸索和嘗試，建立了一套方便快捷、高效穩(wěn)定的農(nóng)桿菌介導(dǎo)的“小米”遺傳轉(zhuǎn)化體系（圖7）。該體系遺傳轉(zhuǎn)化效率高達(dá)23.28%，已經(jīng)接近模式植物擬南芥和水稻的遺傳轉(zhuǎn)化水平，解決了“小米”作為模式植物遺傳轉(zhuǎn)化效率低下的難題。
京公網(wǎng)安備 11011302003368號(hào)