利用CRISPR/Cas9敲除OsNramp5基因創(chuàng)制低鎘秈稻

龍起樟 黃永蘭 唐秀英 王會(huì)民 蘆明 袁林峰 萬建林

利用CRISPR/Cas9敲除基因創(chuàng)制低鎘秈稻

龍起樟 黃永蘭 唐秀英 王會(huì)民 蘆明 袁林峰 萬建林*

[江西省超級(jí)水稻研究發(fā)展中心/水稻國(guó)家工程實(shí)驗(yàn)室（南昌）/江西省綠色水稻種質(zhì)重點(diǎn)實(shí)驗(yàn)室，南昌 330200；*通訊聯(lián)系人，E-mail: ncwanjl@163.com]

為了嘗試快速培育低鎘秈稻，選取廣泛應(yīng)用的雜交水稻親本華占和五豐B以及常規(guī)品種五山絲苗和中早35為材料，通過CRISPR/Cas9技術(shù)創(chuàng)制基因突變株系，在鎘污染及正常土壤中種植并測(cè)定突變株系籽粒（糙米）鎘含量，其他相關(guān)元素含量亦同時(shí)在鎘污染土壤種植條件下測(cè)定，在非鎘污染土壤種植條件下考查基因敲除對(duì)農(nóng)藝性狀的影響。成功獲得了前述品種的基因敲除株系；非鎘污染條件下種植的4個(gè)品種基因敲除株系籽粒鎘含量低于0.02 mg/kg，平均較野生型降低85.5%；而在鎘污染土壤種植時(shí)，不同品種基因敲除株系籽粒鎘含量低于0.1 mg/kg，平均比野生型降低94.8%；錳含量也降低52.7%，鉻含量增加59.5%，鉛含量在華占中增加79.1%，而在其他品種中無變化；銅、鐵、鋅、鈣、硒和砷含量（后4種元素只在華占及衍生品系中檢測(cè)）受影響較小或不受影響；基因敲除株系株高、結(jié)實(shí)率和千粒重較野生型小幅降低，而有效分蘗略微增加，產(chǎn)量平均減少6.9%。通過基因敲除，可以顯著降低鎘積累，但在某些種植條件下，代價(jià)為小幅產(chǎn)量損失；通過本研究獲得的低鎘基因敲除品系在鎘污染地區(qū)具有較好利用潛力。

秈稻；鎘；基因編輯；育種

水稻是一種易于富集鎘的作物[1]，而耕地酸化或污染導(dǎo)致土壤中活性鎘含量增加，又加劇了水稻中鎘的積累[2]，導(dǎo)致稻米鎘超標(biāo)事件頻發(fā)。過量攝入鎘嚴(yán)重危害身體健康[3]，在我國(guó)南方部分地區(qū)，稻米鎘超標(biāo)已經(jīng)成為引起社會(huì)廣泛關(guān)注的食品安全問題。耕作方式的變革以及污染土壤的修復(fù)是解決稻米鎘污染問題的根本途徑，而選擇種植鎘低積累水稻品種可起到立竿見影的效果，是目前低成本快速解決稻米鎘污染問題的最佳選擇。研究者從現(xiàn)有水稻品種資源中篩選獲得了一些鎘低積累水稻品種[4-7]，這可以解決部分稻米鎘污染的問題，然而受水稻品種的生態(tài)適應(yīng)性限制，獲得的低鎘品種資源不適合在所有鎘污染區(qū)推廣，或者獲得的資源降鎘幅度有限，不適合在鎘污染程度較重的地區(qū)種植，因此生產(chǎn)中迫切需要適應(yīng)不同地區(qū)生產(chǎn)需求的低鎘水稻新品種。

水稻中鎘吸收轉(zhuǎn)運(yùn)相關(guān)蛋白及調(diào)控因子的鑒定為育種家通過分子手段快速培育低鎘水稻新品種奠定了基礎(chǔ)[2, 8-9]。研究顯示一些關(guān)鍵的鎘吸收轉(zhuǎn)運(yùn)蛋白編碼基因可作為分子操作對(duì)象以創(chuàng)制低鎘水稻品種。OsHMA2[10-12]是一種鋅的轉(zhuǎn)運(yùn)蛋白，參與鎘從根部向莖部的轉(zhuǎn)移，同時(shí)在莖節(jié)中通過韌皮部向發(fā)育器官轉(zhuǎn)運(yùn)鎘，其功能缺失可使鎘和鋅同時(shí)減少，而且嚴(yán)重影響生長(zhǎng)和產(chǎn)量，但過表達(dá)可增加籽粒中鋅含量并一定程度降低鎘積累；OsHMA3[13]是一種鎘的轉(zhuǎn)運(yùn)蛋白，位于液泡膜上向液泡內(nèi)轉(zhuǎn)運(yùn)鎘，其功能是限制鎘從根向莖葉轉(zhuǎn)移從而將其截留在根內(nèi)液泡中，此基因的超量表達(dá)可大幅降低籽粒鎘含量，而功能缺失將使地上部分鎘猛增；OsLCT1[14]亦是一種鎘的細(xì)胞膜離子外排轉(zhuǎn)運(yùn)蛋白，同時(shí)可轉(zhuǎn)運(yùn)鈣、鎂、鉀和錳，其參與鎘從莖葉向籽粒轉(zhuǎn)移過程，通過RNA干涉抑制其功能可使籽粒中鎘含量減少約50%，同時(shí)鐵含量增加，而且不影響植物生長(zhǎng)及產(chǎn)量，亦不影響其他有益礦物元素的積累；OsLCD[15]的具體功能未知，其功能缺失可降低籽粒中約50%的鎘含量，不影響生長(zhǎng)和產(chǎn)量；OsNramp5[16-18]是一種錳和鐵的轉(zhuǎn)運(yùn)蛋白，是根中鎘吸收的主要蛋白，其功能缺失可使籽粒中鎘含量大幅降低（含量不到對(duì)照的3%），同時(shí)錳含量也降低一半以上，因品種或種植條件不同，OsNramp5功能缺失突變體可表現(xiàn)出大幅減產(chǎn)或產(chǎn)量不受影響兩種截然不同的結(jié)果。

CRISPR/Cas9等基因編輯技術(shù)的快速發(fā)展使得植物基因功能缺失突變體的創(chuàng)制變得十分便利[19]，因此，可以通過基因編輯技術(shù)創(chuàng)制鎘吸收蛋白編碼基因的突變體以快速培育低鎘水稻品種。由前所述，、和基因功能缺失可使水稻籽粒鎘積累減少，而且在一定條件下對(duì)其他有益元素的積累和產(chǎn)量無影響或影響較小，因而可以作為基因敲除靶標(biāo)以創(chuàng)制低鎘水稻品種。本研究旨在解決南方稻米嚴(yán)重鎘超標(biāo)（如鎘含量1 mg/kg以上）地區(qū)的水稻生產(chǎn)問題，故選擇缺失可大幅降鎘的基因作為靶標(biāo)，通過CRISPR/Cas9方法對(duì)其進(jìn)行基因敲除，檢測(cè)基因失活對(duì)鎘和其他元素積累的影響，并考查突變體株高、產(chǎn)量和品質(zhì)性狀的變化，評(píng)價(jià)其利用價(jià)值。

1 材料與方法

1.1 水稻品種

Kasalath由南京農(nóng)業(yè)大學(xué)水稻研究所提供；華占和五豐B分別源于中國(guó)水稻研究所和廣東農(nóng)業(yè)科學(xué)院水稻研究所；五山絲苗和中早35購自種子公司。

1.2 載體構(gòu)建

CRISPR/Cas9載體為自行構(gòu)建，載體以雙元載體pCUbi1390(由pCAMBIA1390骨架上添加玉米Ubiquitin 3啟動(dòng)子構(gòu)建而成)[20]為基礎(chǔ)，在I和HⅠ位點(diǎn)處插入水稻密碼子優(yōu)化過的Cas9基因，繼續(xù)在Ubiquitin 3啟動(dòng)子前的d Ⅲ位點(diǎn)插入方向與基因同向的水稻U6或U3啟動(dòng)子和sgRNA骨架構(gòu)件（構(gòu)件中啟動(dòng)子和sgRNA框架之間插入一段兩端為Ⅰ位點(diǎn)的R-基因表達(dá)盒以便于克隆sgRNA向?qū)蛄衅?；?gòu)件中包含U6或U3啟動(dòng)子的轉(zhuǎn)錄起始位點(diǎn)堿基G或A），最終載體命名為pCUbi1390Cas9-U6和pCUbi1390Cas9-U3（圖1），載體在大腸桿菌ccdB Survival? 2 T1R菌株（Invitrogen公司）中保存。其中，含有水稻密碼子優(yōu)化過的Cas9基因原始載體由瞿禮嘉實(shí)驗(yàn)室惠贈(zèng)[21]；水稻U6和U3啟動(dòng)子擴(kuò)增自日本晴基因組DNA，分別參考Feng等[22]和Miao等[21]的數(shù)據(jù)；sgRNA骨架來自張峰實(shí)驗(yàn)室[23]；R-基因表達(dá)盒來自從Invitrogen公司購置的gateway載體pENTR 1A。

利用在線工具CRISPR-P（http://cbi.hzau.edu.cn/ crispr/）[24, 25]設(shè)計(jì)用于敲除基因（RAP ID: Os07g0257200）的sgRNA向?qū)蛄?。選擇兩個(gè)靶位點(diǎn)，均位于基因第10外顯子的正義鏈上[基因組ORF位置：3673?3692（位點(diǎn)1）和3685?3704（位點(diǎn)2），mRNA CDS位置：1006?1025（位點(diǎn)1）和1018?1037（位點(diǎn)2）]。構(gòu)建載體時(shí)，先根據(jù)載體構(gòu)建要求合成針對(duì)兩個(gè)靶位點(diǎn)的向?qū)蛄蠨NA接頭寡核酸，針對(duì)位點(diǎn)1為cttgGGCAGAGCTCC ACTATTAC（5′-3′；正義鏈）和aaacGTAATAGTGGA GCTCTGCC（5′-3′；反義鏈）；針對(duì)位點(diǎn)2為ctagACTATTACCGGCACATACGCTGG（5′-3′；正義鏈）和tcgaCCAGCGTATGTGCCGGTAATAGT（5′-3′；反義鏈）。將兩條寡核酸退火后即為接頭，將接頭連接到用Ⅰ酶（Thermo scientific?；貨號(hào)ER1581）線性化的CRISPR/Cas9載體（第一個(gè)位點(diǎn)使用pCUbi1390Cas9-U6，第二個(gè)位點(diǎn)使用pCUbi1390Cas9-U3）后轉(zhuǎn)化大腸桿菌，挑單克隆測(cè)序，測(cè)序引物為OsU6-F（TTGAGCGATTACAGGC GAAAGTG；用于pCUbi1390Cas9-U6載體，檢測(cè)向?qū)蛄械恼_性），OsU3-F（GGCATGCATG GATCTTGGAGGAAT；用于pCUbi1390Cas9-U3載體，檢測(cè)向?qū)蛄姓_性），35S-F（TGACGCA CAATCCCACTATCCTTC；檢測(cè)載體完整性），Cas9-R-1（TCGAGCCTGCGGGACTTAGAG；Cas9 5′端引物，檢測(cè)載體完整性）；C126（TCGTGA AGAAGACCGAGGTT；Cas9 3′端引物，檢測(cè)載體完整性，此引物由Miao等設(shè)計(jì)[21]），選取正確克隆擴(kuò)大繁殖并提取質(zhì)粒備用。

LB、RB分別為T-DNA的左、右邊界；HPTⅡ?潮霉素抗性標(biāo)記；2×35S?串聯(lián)的35S啟動(dòng)子；CmR-ccdBcassette?氯霉素抗性基因及自殺基因表達(dá)盒；sgRNA scaffold & Poly (T)7?sgRNA骨架及7個(gè)T堿基組成的終止子；ZmUbi1 promoter?玉米Ubi1基因啟動(dòng)子；Rice codon optimized Cas9?水稻密碼子優(yōu)化的Cas9基因；AarⅠ?供向?qū)蛄锌寺〉腡ypeⅡ限制性酶識(shí)別位點(diǎn)（黑色三角示剪切點(diǎn)）；PstⅠ-Sna BⅠ-MluⅠ?三個(gè)單一限制性內(nèi)切酶識(shí)別位點(diǎn)；圖中省略了HPT和Cas9基因各自表達(dá)所需的終止子CaMV 3′ UTR和Nos polyA信號(hào)序列；pCUbi1390Cas9-U6和pCUbi1390Cas9-U3載體的區(qū)別僅在于sgRNA啟動(dòng)子及轉(zhuǎn)錄起始?jí)A基不同。

Fig. 1. Structure of the T-DNA region of the CRISPR/Cas9 Ti binary vectors pCUbi1390Cas9-U6/U3.

1.3 遺傳轉(zhuǎn)化

將構(gòu)建好的Ti雙元質(zhì)粒載體轉(zhuǎn)化至EHA105農(nóng)桿菌中，農(nóng)桿菌感受態(tài)細(xì)胞制備方法采用CaCl2法，轉(zhuǎn)化采用反復(fù)凍融法[26]。對(duì)于品種Kasalath，遺傳轉(zhuǎn)化采用Toki等的方法進(jìn)行[27]；華占、五豐B、五山絲苗和中早35的轉(zhuǎn)化委托武漢伯遠(yuǎn)生物科技有限公司完成。

1.4 水稻DNA提取

DNA提取用CTAB法[28]，省略部分純化步驟。

1.5 轉(zhuǎn)基因陽性株篩選

用葉片切段潮霉素水溶液浸泡法[29]檢測(cè)T-DNA是否插入T0再生苗中，接著檢測(cè)篩除T1代轉(zhuǎn)基因陽性植株。

1.6 CRISPR/Cas9靶位點(diǎn)突變檢測(cè)

通過PCR產(chǎn)物直接測(cè)序的方法對(duì)轉(zhuǎn)基因陽性植株中的基因敲除靶位點(diǎn)進(jìn)行初步檢測(cè)，若測(cè)序結(jié)果峰圖中有套峰則認(rèn)為存在突變，繼續(xù)將片段克隆至質(zhì)粒載體中進(jìn)行測(cè)序，每個(gè)樣品至少5個(gè)陽性克隆，從中分析出兩個(gè)等位基因各自的基因型，若測(cè)序結(jié)果無套峰，則與對(duì)照進(jìn)行序列比對(duì)，獲悉基因型。PCR引物為TTCGTGGCGCTGCTGATAAAC（正向）和AGAGCGGAGAAATATGGACGAAAGT（反向）。PCR體系包括5 μL 10×PCR緩沖液（含Mg2+；TaKaRa），4 μL dNTP混合液（各2.5 mmol/L），5 μL引物（正反向混合，各5 pmol/μL），0.5 μL TaKaRa（5 U/μL），2 μL二甲基亞砜（濃度不低于99.5%），2 μL DNA，31.5 μL H2O，總體積50 μL。PCR程序：95℃下2 min；94℃下30 s，58℃下30 s，72℃下60 s，35個(gè)循環(huán)；72℃下2 min。PCR產(chǎn)物直接測(cè)序所用引物為（擴(kuò)增片段內(nèi)，位于反義鏈）GTGCACCCCTACAATTCGTCAGT。測(cè)序委托杭州擎科生物技術(shù)有限公司完成。

1.7 元素含量分析

鎘、錳、鐵、銅、鉛、鉻、汞、砷、硒、鋅和鈣元素含量測(cè)定委托具有資質(zhì)的專業(yè)第三方檢測(cè)機(jī)構(gòu)青島科創(chuàng)質(zhì)量檢測(cè)有限公司完成。田間按株系及姊妹系各自混收稻谷，脫殼成糙米（曬干，水分12%左右）后檢測(cè)，每份樣品不少于20 g，將樣品進(jìn)行粉碎（過100目篩）后送檢。Cd測(cè)定樣品中部分抽取2個(gè)重復(fù)樣檢測(cè)，最后取其均值作測(cè)定值，其余樣品只測(cè)定一份樣品。元素檢測(cè)方法采用ICP-MS（電感耦合等離子體質(zhì)譜）法。

1.8 材料繁殖及農(nóng)藝性狀測(cè)定

T0代植株溫室內(nèi)盆缽種植；T1分離群體在有防護(hù)措施的網(wǎng)室內(nèi)種植，開花前篩選后拔除轉(zhuǎn)基因植株；已去除轉(zhuǎn)基因成分的T2或T3植株在大田種植，并用塑料薄膜墻圍住隔離；供鎘測(cè)定材料在晚季種植（2017年6?11月），供農(nóng)藝性狀考查的中早35和五豐B品系作早稻（2018年4?7月）種植，華占和五山絲苗作中稻種植（2018年5?10月）。所有世代種植行距20 cm、株距17 cm，每行8株。元素含量測(cè)定所用材料種植時(shí)野生型對(duì)照設(shè)3次重復(fù)，在田間均勻排布，基因敲除株系按株系順序種植，同一轉(zhuǎn)化事件后代攜相同變異姊妹系視為重復(fù)；農(nóng)藝性狀調(diào)查所用材料種植時(shí)不設(shè)重復(fù)。所有材料均按大田常規(guī)方法進(jìn)行肥水管理。株高測(cè)定取莖稈基部至最長(zhǎng)穗頂部的長(zhǎng)度，有效分蘗測(cè)定取穗子實(shí)粒不少于5粒的分蘗，株高和有效分蘗數(shù)同時(shí)測(cè)定，連續(xù)測(cè)定一定株數(shù)（一般大于40株），邊株不測(cè)定；每穗總粒數(shù)和實(shí)粒數(shù)測(cè)定取樣10株，取樣時(shí)連續(xù)選10行，每行取一株，逐行錯(cuò)位選株（即第1，2，3，4，5，6，7，8，9和10行分取第2，3，4，5，6，7，6，5，4和3株），邊株不取，穗子混合，每穗分別統(tǒng)計(jì)；千粒重測(cè)定取曬干稻谷去除空癟粒和雜質(zhì)后用自動(dòng)數(shù)粒儀數(shù)1000粒，重復(fù)2次；米質(zhì)測(cè)定采用株系內(nèi)所有植株的混合稻谷，經(jīng)研磨成精米后測(cè)定，外觀品質(zhì)指標(biāo)通過萬深SC-E大米外觀品質(zhì)檢測(cè)分析儀進(jìn)行測(cè)定，其中，堊白度指堊白米粒中堊白面積與整粒米面積的比值，水分、蛋白含量、直鏈淀粉含量和食味分用近紅外分析儀檢測(cè)，儀器為InfratecTM1241（Foss），定量模型為廠商提供并經(jīng)過用戶矯正，數(shù)值仍為估計(jì)值。

1.9 HRM分析

高分辨率溶解曲線分析利用Lightscanner 96系統(tǒng)(Idaho科技公司)進(jìn)行，采用小片段擴(kuò)增法。PCR引物及擴(kuò)增子為CTGGGCAAGTCGAGTGCGAT（正向）和ATTACCTGCATGATGTACTGTCC（反向），擴(kuò)增產(chǎn)物104 bp。PCR體系包括1 μL 10×PCR 緩沖液（含Mg2+；TaKaRa），0.8 μL dNTP混合液（各2.5 mmol/L），1 μL引物（正反向各5 pmol/μL），0.1 μL TaKaRa（5 U/μL），0.4 μL 二甲基亞砜（濃度不低于99.5%），1 μL DNA，1 μL LC-Green熒光染料，4.7 μL H2O，總體積10 μL，上加20 μL礦物油。PCR程序：94℃下2 min；95℃下5 s，60℃下20 s，72℃下20 s，共50個(gè)循環(huán)；72℃下1 min；95℃下3 min；16℃下2 min。將PCR體系連同石蠟油全部轉(zhuǎn)入專用PCR板（Bio-rad，貨號(hào)hsp9665）中進(jìn)行HRM分析，分析溫度范圍為65℃~95℃。用系統(tǒng)自帶軟件分析數(shù)據(jù)。

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

數(shù)據(jù)分析利用Microsoft Excel 2010及其共享Real Statistics統(tǒng)計(jì)插件（http://www.real-statistics. com/）進(jìn)行。統(tǒng)計(jì)分析時(shí)，除特殊注明外，若比較分析中只含2組樣本數(shù)據(jù)，直接進(jìn)行測(cè)驗(yàn)檢測(cè)方差同質(zhì)性后利用Excel自帶函數(shù)進(jìn)行檢測(cè)（根據(jù)測(cè)驗(yàn)結(jié)果選等方差或異方差模式）；若包含3組以上樣本數(shù)據(jù)，則進(jìn)行多重比較，利用Real Statistics統(tǒng)計(jì)插件行中位數(shù)版Levene檢驗(yàn)檢測(cè)方差同質(zhì)性，若同質(zhì)后行方差分析（ANOVA）檢測(cè)組間差異，若結(jié)果顯著，采用Dunnett檢驗(yàn)檢測(cè)突變體與對(duì)照差異；若方差不同質(zhì)，采用Welch ANOVA檢測(cè)組間差異，若顯著,則采用Games Howell test檢測(cè)突變體與對(duì)照的差異顯著性。

不同品種突變體與對(duì)照相比的不同指標(biāo)增減率數(shù)據(jù)為品種內(nèi)各株系的平均值與各自野生型對(duì)照的平均值比較而得，4個(gè)品種總的增減率為4個(gè)品種各自平均增減率的均值。

2 結(jié)果與分析

2.1 OsNramp5基因編輯水稻品系的獲得

先以品種Kasalath為材料試驗(yàn)CRISPR/Cas9系統(tǒng)對(duì)基因兩個(gè)靶位點(diǎn)的編輯效果。結(jié)果顯示在兩個(gè)靶位點(diǎn)中均成功創(chuàng)制突變，選擇位點(diǎn)1（CRISPR/Cas9系統(tǒng)中使用U6啟動(dòng)子）作為進(jìn)一步試驗(yàn)的基因編輯靶位點(diǎn)，擴(kuò)大試驗(yàn)對(duì)象至純系品種五山絲苗和中早35以及雜種親本華占和五豐B。在針對(duì)靶位點(diǎn)1的試驗(yàn)中，5個(gè)品種一共獲得97個(gè)獨(dú)立轉(zhuǎn)基因株系，對(duì)T0代（轉(zhuǎn)基因當(dāng)代定義為T0代）植株的測(cè)序顯示，不同品種基因突變效率不同，突變率（突變轉(zhuǎn)化子數(shù)/總轉(zhuǎn)化子數(shù)×100%）為39%~90%，平均為58%（表1）。對(duì)這5個(gè)品種一共56個(gè)獨(dú)立基因編輯株系中的106個(gè)突變事件進(jìn)行了分析，結(jié)果顯示一共出現(xiàn)了32種突變基因型（圖2-A），理論上其中25種突變類型為移碼突變，導(dǎo)致蛋白功能失活，7種突變類型（即當(dāng)缺失或插入的堿基造成的總的堿基變化數(shù)為3的整倍數(shù)時(shí)，包括#05，#07，#12，#13，#17，#22和#24）導(dǎo)致部分氨基酸殘基發(fā)生改變，蛋白功能是否受影響未知。每個(gè)基因編輯株系所攜突變基因型見表2。針對(duì)靶位點(diǎn)2的基因編輯結(jié)果也在表1、表2和圖2-B中顯示。

進(jìn)一步對(duì)靶位點(diǎn)1所有突變基因型的發(fā)生頻率進(jìn)行了統(tǒng)計(jì)（圖2-C），發(fā)現(xiàn)主要發(fā)生的突變類型為缺失1~3個(gè)堿基和插入1個(gè)堿基的突變，占所有突變的66%，而且在這3種缺失突變中缺失的堿基位置并不具隨機(jī)性，而是缺失特定位置的堿基，而在插入1個(gè)堿基的突變類型中又以插入堿基A和T最為常見，未發(fā)現(xiàn)插入堿基C的突變類型；缺失5個(gè)堿基及以上的缺失突變以及插入缺失型的突變中每一種突變類型一般只發(fā)生1次，極少數(shù)發(fā)生2次。

在轉(zhuǎn)基因T1代通過葉片潮霉素抗性篩選獲得去除轉(zhuǎn)基因成分的植株，對(duì)這些植株中的靶位點(diǎn)進(jìn)行測(cè)序檢測(cè)，篩選獲得純合突變株，利用這些純合株的后代（T2或T3）進(jìn)行后續(xù)試驗(yàn)。目前為止除品種Kasalath外，所有品種的多數(shù)突變基因型都獲得了純合后代，詳細(xì)信息見表2。

2.2 OsNramp5基因敲除降低水稻籽粒中鎘積累

在獲得去除轉(zhuǎn)基因成分的基因突變株系后，將品種華占、五豐B、五山絲苗和中早35的部分基因突變株系（T2代，除華占品系L01-#05/#05外余均為移碼突變）及其野生型對(duì)照分別在江西省A、B兩地種植（A地為鎘中等污染土壤，而B地為非污染土壤）。分別對(duì)從兩地種植獲得的種子（糙米）進(jìn)行鎘含量測(cè)定，測(cè)試結(jié)果見表3。由表中數(shù)據(jù)可知，A、B兩地種植的基因編輯株系種子鎘含量較對(duì)照均大幅下降。在A地種植時(shí)，在對(duì)照品種鎘含量嚴(yán)重超標(biāo)——約5~10倍于國(guó)標(biāo)限定值（0.2 mg/kg；GB2762?2017）的情況下，不同品種敲除株系的鎘含量不及限定值的1/2，比對(duì)照降低91.7%~96.9%（4個(gè)品種的均值為94.8%）；在B地種植時(shí)，對(duì)照和基因敲除株系鎘含量均低于國(guó)標(biāo)限定值，其中不同品種基因敲除株系比對(duì)照降低76.2%~91.1%（均值85.5%），不到國(guó)標(biāo)限定值的1/10。值得一提的是，華占品系L01-#05/#05鎘含量的降低與在其他品系相同，因此，雖然此株系內(nèi)氨基酸序列只發(fā)生幾個(gè)氨基酸殘基的改變，同樣造成了基因功能的失活。總而言之，對(duì)于不同類型的水稻品種，敲除基因均可有效降低籽粒鎘積累。

表3數(shù)據(jù)亦顯示，無論在何地種植，不同品種野生型本身鎘積累存在一定差異，尤其是中早35，其鎘含量明顯比其他3個(gè)品種低，而敲除基因后，不同品種的鎘積累量變得相對(duì)一致。這一結(jié)果暗示了基因的功能差異可能是引起這些品種鎘積累差異的關(guān)鍵原因。

表1 不同品種陽性轉(zhuǎn)化株基因突變情況統(tǒng)計(jì)

A和B，靶位點(diǎn)1和2的序列突變情況。WT?野生型；#**或##**，突變類型編號(hào)；野生型序列中Cas9靶位點(diǎn)、Cas9系統(tǒng)中U6和U3啟動(dòng)子A圖和B圖分別為靶位點(diǎn)1和靶位點(diǎn)2的突變情況。起始轉(zhuǎn)錄堿基以及PAM位點(diǎn)分別用灰色陰影、加下劃線及加框標(biāo)示，黑色倒三角指示Cas9切割位點(diǎn)，所有序列方向均為5′-3′；突變體序列中缺失堿基用短線代替，插入堿基加下劃線顯示，其余為與野生型相同者；圖右邊數(shù)字表示缺失（用“?”表示）或插入（用“+”表示）的堿基數(shù)。C?靶位點(diǎn)1不同突變類型發(fā)生頻率統(tǒng)計(jì)。突變體類型分3類，即包含缺失1 bp、2 bp、3 bp、4 bp和≥6 bp六種類型的缺失突變，包含插入1個(gè)堿基A、T或G及插入2個(gè)堿基TT四種類型的插入突變，以及同時(shí)含有缺失和插入即InDel；一共統(tǒng)計(jì)56個(gè)獨(dú)立轉(zhuǎn)化株系，106個(gè)突變等位基因。

A and B, The sequence changes in Target sites 1 and 2 of thegene, respectively. WT, Wild type; #** or ##**, Mutation type numbers; In the WT sequence, the target site, the cut site and the PAM site are indicated with the bases highlighted in grey (The underlined base “G” or “A” indicates the transcription start site of the rice U6 or U3 promoter in the CRISPR/Cas9 system), reverse black triangle and the bases in underlined bold font, respectively; In the mutation sequences, the deleted bases are indicated by dashed lines, the inserted ones are underlined, and those identical with the wild type are shown normally; The numbers on the right indicate the counts of the deleted or inserted bases. C, The occurrence frequency of various mutation types in Target site 1. The mutations were summarized as three types, which are deletions, insertions and InDels. The deletion type includes five subtypes, which are those with 1 bp, 2 bp, 3 bp, 4 bp and ≥6 bp deletion respectively. The insertion type includes four subtypes, which are those with one A, one T, one G and double T insertion, respectively. The data are derived from the analysis of 106 mutation alleles from 56 independent transgenic lines.

圖2 基因編輯株系中基因靶位點(diǎn)突變情況

Fig. 2. Target site mutations ofin the genetically edited lines.

表2 五個(gè)品種基因編輯株系中OsNramp5靶位點(diǎn)的基因型

灰色底紋標(biāo)識(shí)的單元格表示該基因型已經(jīng)分離到去除標(biāo)記且純合的植株；其他基因型正在分離純合植株；?，未確定；?，不適用。

The grey shading contents show that marker-free homozygous plants are available for the genotypes. For the other genotypes, the work to seek such plants was in progress; ?, Not resolved; ?, Not applicable.

2.3 OsNramp5基因敲除對(duì)籽粒中其他元素積累的影響

對(duì)鎘污染土壤（A地）種植的基因敲除株系及對(duì)照進(jìn)行籽粒鎘含量測(cè)定的同時(shí)，也對(duì)錳、鐵、銅、鉛和鉻元素含量進(jìn)行了測(cè)定。由表3可知，基因敲除后，所有品種錳含量均顯著降低，減幅34.2%~63.1%（均值52.7%）；與此相反，鐵和鉻含量一定程度增加，不同品種鐵含量增幅11.4%~34.5%（均值21.4%），而鉻含量增幅更大，達(dá)到45.6%~75.5%（均值59.5%），但是，差異顯著性測(cè)驗(yàn)結(jié)果顯示，無論是鐵還是鉻，這種增加只在部分品種中顯著；鉛元素含量在華占中顯著增加，增幅達(dá)79.1%，由于其值已略微超出國(guó)家限定值（0.2 mg/kg；GB2762?2017），這種變化值得關(guān)注，而其他品種鉛含量并未因基因敲除而顯著改變；銅元素含量變化較?。▋H在中早35基因敲除株系中減少了13.3%）。此外，對(duì)華占基因敲除株系及野生型的汞、砷、硒、鋅和鈣元素含量進(jìn)行檢測(cè)（數(shù)據(jù)未顯示），其中汞元素含量在所有樣品中均低于所用方法的檢測(cè)下限（30 μg/kg），因而基因敲除對(duì)汞積累的影響未知，而砷、硒、鋅和鈣的含量也無顯著變化。

表3 OsNramp5基因敲除株系籽粒的礦物元素含量變化

1在鎘污染土壤種植；2在正常土壤種植；3數(shù)據(jù)為相同品種不同突變基因型觀測(cè)值的均值。樣品編號(hào)格式為“品種名縮寫-株系號(hào)-基因型”，具體信息見表2和圖2。HZ?華占；WF?五豐B；WS?五山絲苗；ZZ?中早35；WT?野生型；MUT?OsNramp5基因敲除株系。

數(shù)據(jù)格式為“均值±標(biāo)準(zhǔn)差(n)”或“均值(n)”；*，P<0.05；**，P<0.01；差異顯著性通過t測(cè)驗(yàn)（雙樣本等方差或雙樣本異方差）進(jìn)行檢測(cè)（先進(jìn)行F測(cè)驗(yàn)檢測(cè)方差齊性），統(tǒng)計(jì)分析時(shí)假設(shè)相同品種不同株系相互之間等同。

1Grown in a Cd-contaminated field;2Grown in a non-Cd-contaminated field;3Means of all the observed values of the genetic editing lines. The sample name is “cultivar name abbreviation-line number-genotype”, and one may refer to Table 2 and Fig. 2 for the information of the line numbers and genotypes. HZ, Huazhan; WF, Wufeng B; WS, Wushansimiao; ZZ, Zhongzao 35; WT, Wild type; MUT, Mutants(All the OsNramp5 mutation lines).

Values are “Mean±SD(n)” or “Mean(n)”;*, P<0.05; **, P<0.01; The significance was determined by Student’s t test (types: equal variances or unequal variances) after F test for variance homogeneity. In the analysis of statistical significance, all the genetic editing lines within one cultivar were considered to be identical.

被敲除后，水稻籽粒中鎘含量的大幅下降，錳含量亦減少一半左右，而鉻含量增加一半以上，鐵、鋅、鈣、硒、砷和銅積累受影響較小或不受影響，鉛可能在個(gè)別品種中有一定程度增加。

2.4 OsNramp5基因敲除對(duì)農(nóng)藝性狀的影響

為了檢測(cè)敲除對(duì)農(nóng)藝性狀的影響，對(duì)大田（非鎘污染土壤）種植的部分低鎘基因敲除株系（T3代）及野生型對(duì)照的株高和產(chǎn)量性狀進(jìn)行了測(cè)定（表4）。4個(gè)品種基因敲除株系（同一品種內(nèi)所有株系的指標(biāo)均值）與對(duì)照相比發(fā)生了如下變化：株高稍降低，降幅3.5%~6.8%（4個(gè)品種的降幅均值為4.9%），所有株系與對(duì)照相比均差異顯著；有效分蘗呈增加趨勢(shì)，增幅1.9%~9.0%（均值5.5%），但由于分蘗性狀誤差較大，可能差異為誤差所掩蓋，差異并不顯著，但由于調(diào)查的樣本較大，我們認(rèn)為這種增加是真實(shí)存在的；每穗總粒數(shù)變化無明顯規(guī)律；每穗實(shí)粒數(shù)、結(jié)實(shí)率和千粒重均呈下降趨勢(shì)，降幅分別為6.5%~10.3%（均值7.9%）、3.2%~9.0%（中早35基因敲除株系除外；均值5.3%）和1.4%~5.5%（均值3.4%），不過這種下降并非所有株系都差異顯著。通過計(jì)算理論產(chǎn)量，發(fā)現(xiàn)4個(gè)品種減產(chǎn)5.4%~9.4%（均值6.9%）。概而言之，敲除可引起株高、實(shí)粒數(shù)、結(jié)實(shí)率和千粒重輕度下降，而有效分蘗小幅增加，最終導(dǎo)致不容忽視的小幅減產(chǎn)。

表4 OsNramp5基因敲除對(duì)株高、產(chǎn)量性狀的影響

數(shù)據(jù)格式為均值±標(biāo)準(zhǔn)差（）；*，<0.05；**，<0.01；?，未檢測(cè)或不適用；突變體與野生型之間顯著性測(cè)驗(yàn)視適用情況采用測(cè)驗(yàn)（含兩組數(shù)據(jù)時(shí)用；行測(cè)驗(yàn)比較方差同質(zhì)性后進(jìn)行）或方差分析加多重比較（含三組以上數(shù)據(jù)時(shí)用；等方差時(shí)用ANOVA+Dunnett檢驗(yàn)，異方差時(shí)用Welch ANOVA +Games-Howell test，組間方差齊性用中位數(shù)版Levene檢驗(yàn)進(jìn)行檢測(cè)，千粒重?cái)?shù)據(jù)因<3未檢測(cè)方差齊性，視為等方差）；結(jié)實(shí)率數(shù)據(jù)經(jīng)平方根反正弦轉(zhuǎn)換后進(jìn)行方差分析；產(chǎn)量依據(jù)“實(shí)粒數(shù)×有效分蘗數(shù)×千粒重/1000×270 000株/1000”估算；樣品編號(hào)格式為“品種名縮寫-株系號(hào)-基因型”，基因型WS-L01-#28/#28測(cè)試兩個(gè)姊妹系。品種名縮寫：HZ?華占；WF?五豐B；WS?五山絲苗；ZZ?中早35；WT?野生型。株系號(hào)見表2，基因型編號(hào)見圖2。

Values are Mean± SD ();*,<0.05;**,<0.01; ?，Not applicable or not determined; the significances between mutants and wild types were determined bytests aftertests for variance homogeneity (for comparisons with two groups of samples and for comparisons of the 1000-grain weight trait) or by ANOVA with post hoc multiple comparison analysis (for comparisons with three or more groups of samples; ANOVA + Dunnett’stests for groups with equal variances; Welch’s ANOVA + Games-Howell tests for groups with unequal variance; Median-edition Levene’s test was used to determine the homogeneity of group variances expect those of 1000-grain weight because of<3; groups of 1000-grain weight data were viewed as having equal variances); arcsine square root transformation for the data of seed setting rate was performed before ANOVA; the theoretical yield was estimated by “filled grain number per panicle × effective tiller number × 1000-grain weight / 1000 × 270 000 / 1000”. The sample name is “cultivar name abbreviation-line number-genotype”, and the information of the line numbers and genotypes is available in Table 2 and Fig. 2. Two sister lines were analyzed for the genotype WS-L01-#28/#28. Abbreviations of cultivar names: HZ, Huazhan; WF, Wufeng B; WS, Wushansimiao; ZZ, Zhongzao 35; WT, Wild type. The sample name is “cultivar name abbreviation-line number-genotype”. One may refer to Table 2 and Fig. 2 for the information of the line numbers and genotypes.

為了檢測(cè)敲除對(duì)大米品質(zhì)的影響，我們亦對(duì)敲除突變體及對(duì)照米樣的外觀品質(zhì)進(jìn)行了測(cè)定，同時(shí)利用近紅外分析對(duì)蛋白含量、直鏈淀粉含量及食味進(jìn)行了估計(jì)（表5）。4個(gè)品種敲除株系與對(duì)照相比，粒長(zhǎng)輕微增加（增幅0.4%~2.7%，均值1.5%），粒寬略微變窄（減幅0.9%~3.9%，均值2.5%），最終所有株系均表現(xiàn)出長(zhǎng)寬比輕微增加（增幅1.8%~5.0%，均值4.0%），而這種長(zhǎng)寬比的增加引起堊白粒率和堊白度的小幅下降（降幅分別為6.7%~22.6%和0.6%~17.6%，均值分別為13.3%和10.0%）。食味指標(biāo)測(cè)定結(jié)果顯示，與對(duì)照相比，4個(gè)品種敲除株系中蛋白質(zhì)含量輕微增加（增幅1.8%~4.3%，均值3.2%），直鏈淀粉含量除品種五山絲苗外，其余品種均不同程度下降（降幅3.5%~12.6%，均值9.5%），而食味分幾乎不受敲除的影響（只有五山絲苗品種數(shù)據(jù)）?？傊贸芤鹱蚜ｉL(zhǎng)寬比小幅增加并一定程度降低堊白發(fā)生，同時(shí)引起蛋白含量輕微增加，而且在一些品種中導(dǎo)致直鏈淀粉含量小幅下降，對(duì)食味分幾乎無影。從應(yīng)用角度來看，敲除對(duì)米質(zhì)的影響是有利的。

表5 OsNramp5基因敲除株系的米質(zhì)變化

數(shù)據(jù)格式為均值±標(biāo)準(zhǔn)誤，重復(fù)3~4次，外觀品質(zhì)每次分析整精米500粒以上，食味品質(zhì)每次測(cè)定取5個(gè)次級(jí)樣品的平均值；?，含量超出模型檢測(cè)范圍；*，<0.05；**，<0.01。數(shù)據(jù)分析及樣品編號(hào)格式同表4，其中堊白粒率數(shù)據(jù)經(jīng)平方根反正弦轉(zhuǎn)換后進(jìn)行方差分析，同一品種樣品水分含量基本相同。品種名縮寫：HZ?華占；WF?五豐B；WS?五山絲苗；ZZ?中早35；WT?野生型。株系號(hào)見表2，基因型編號(hào)見圖2。

Values are Mean±SE,=3 or 4, for appearance quality more than 500 head rice grains were analyzed in each determination, and for eating quality five subsamples were included in each analysis; ?，Data out of the model range;*，<0.05;**,<0.01. Methods for data analysis are the same as those described in Table 4 and arcsine square root transformation for the data of seed setting rate was performed before ANOVA; The samples within the same cultivar had nearly identical water content. Abbreviations: HZ, Huazhan; WF, Wufeng B; WS, Wushansimiao; ZZ, Zhongzao 35; WT, Wild type. One may refer to Table 2 and Fig. 2 for the information of the line numbers and genotypes.

其中基因型WT（上圖）、WT（下圖）、#03、#05、#06、#09、#14、#26和#27的樣品數(shù)（n）分別為9、13、11、2、4、2、2、1和1。

Fig. 3. Melting curves of various mutation alleles and the wild type in HRM analysis.

2.5 通過HRM法對(duì)OsNramp5基因突變體和野生型進(jìn)行基因分型

高分辨率溶解曲線（HRM）分析是對(duì)SNP和短片段InDel變異進(jìn)行基因分型的高通量低成本方法。本研究獲得的基因突變體可用作低鎘基因供體以進(jìn)行其他低鎘新品種的培育。為了在進(jìn)一步的低鎘性狀轉(zhuǎn)育工作中采用HRM方法進(jìn)行基因分型以提高分子標(biāo)記輔助選擇的效率，本研究進(jìn)行試驗(yàn)，以篩選出易于分型的基因突變類型。結(jié)果顯示，在粗放檢測(cè)條件下（DNA粗放提取、PCR模板濃度不調(diào)為一致以及HRM分析時(shí)不加溫度內(nèi)標(biāo)），編號(hào)為#3, #05, #06, #09, #14, #26和#27（圖3）的突變基因型溶解曲線仍然能與野生型較好區(qū)分，解鏈溫度（m值）差異為0.5℃~1℃。因此在將來的低鎘性狀轉(zhuǎn)育工作中若使用HRM進(jìn)行基因分型，可選擇含有前述幾種基因型的株系作為供體親本。

3 討論

在低鎘水稻創(chuàng)制中，理想的情形是，獲得的低鎘水稻籽粒鎘含量降低的同時(shí)其他有害元素含量不增加，有益元素含量不減少，同時(shí)農(nóng)藝性狀較少或不受負(fù)面影響，而功能失活突變體似乎符合此種要求。

是水稻根部Mn和Cd吸收的關(guān)鍵離子轉(zhuǎn)運(yùn)蛋白，亦在Fe吸收中發(fā)揮作用，這一發(fā)現(xiàn)分別由日本3個(gè)獨(dú)立研究團(tuán)隊(duì)[16-18]所揭示。其中Ishikawa等[16]的研究顯示突變體在低鎘品種培育中具有重要應(yīng)用價(jià)值，作者利用從粳稻品種越光離子束誘變突變體庫中鑒定獲得的基因功能失活突變體進(jìn)行研究，發(fā)現(xiàn)突變能大幅度降低籽粒及稻草鎘積累（含量<0.1 mg/kg）的同時(shí)不改變Fe、Cu和Zn含量，亦不影響包含株高、產(chǎn)量和食味等在內(nèi)的農(nóng)藝性狀。而Sasaki等[18]利用粳稻品種中花11 T-DNA插入突變體進(jìn)行的研究顯示，雖然功能失活使得籽粒及稻草鎘含量大幅度降低，但對(duì)農(nóng)藝性狀的影響與Ishikawa等[16]報(bào)道的截然不同，其植株生長(zhǎng)受到嚴(yán)重抑制并大幅減產(chǎn)（最高減產(chǎn)89%）。突變導(dǎo)致農(nóng)藝性狀嚴(yán)重劣變的現(xiàn)象也為Ishimaru等[17]所報(bào)道。作者利用品種Tsukinohikari的RNA干涉株系進(jìn)行研究，發(fā)現(xiàn)Mn缺乏時(shí)RNA干涉植株與野生型對(duì)照相比生長(zhǎng)受到抑制，但Mn正常時(shí)無差異，揭示Mn的缺乏是導(dǎo)致農(nóng)藝性狀劣變的原因。同樣，Sasaki等[18]通過補(bǔ)充Mn部分逆轉(zhuǎn)苗期生長(zhǎng)抑制表型，進(jìn)一步確認(rèn)了Mn的缺乏導(dǎo)致了農(nóng)藝性狀劣變。Ishikawa等[16]和Sasaki等[18]的研究都是在非Cd污染土壤大田條件下評(píng)價(jià)基因敲除對(duì)農(nóng)藝性狀的影響，顯然兩個(gè)團(tuán)隊(duì)的研究結(jié)果相互矛盾，這一點(diǎn)雖然Ishikawa等在其論文中已經(jīng)提及，但并沒有給出解釋。我們認(rèn)為，這種不一致可能是不同品種對(duì)Mn的需求不同所致，提示了當(dāng)利用功能缺失突變基因?qū)Σ煌贩N進(jìn)行低鎘化改造時(shí)應(yīng)該逐例評(píng)價(jià)其對(duì)農(nóng)藝性狀的影響。

Tang等[30]通過CRISPR/Cas9技術(shù)創(chuàng)制了雜交稻品種親本華占及638S的基因敲除株系，并利用華占基因敲除品系及其與638S基因敲除品系配制獲得的雜種評(píng)價(jià)了突變對(duì)Cd、Mn、Fe、Zn和Cu元素積累以及株高、產(chǎn)量等農(nóng)藝性狀的影響。在本研究中，我們亦在不同類型的秈稻（即早秈中早35、晚秈五山絲苗、三系恢復(fù)系華占和三系保持系五豐B；其中華占品種與Tang等研究中使用的相同）中展示了敲除后水稻籽粒中的Cd含量大幅下降，同時(shí)較大幅度下降的還有Mn含量，而Fe、Zn和Cu含量因品種不同而略微改變或不受影響，這與前述日本研究者在粳稻中的研究結(jié)果以及Tang等[30]在秈稻中的結(jié)果基本一致。與前人的研究不同，我們還測(cè)定了基因敲除株系中元素Pb、Cr、As、Se和Ca的含量變化，發(fā)現(xiàn)Pb在個(gè)別品種（華占）中含量呈增加趨勢(shì)，Cr在4個(gè)品種中亦不同程度增加，As、Se和Ca含量在檢測(cè)的品種華占中并無變化，由于Pb和Cr的轉(zhuǎn)運(yùn)機(jī)制研究較少，尚不清楚其含量增加的具體原因。本研究顯示個(gè)別品種基因敲除株系中Pb含量的增加已經(jīng)略微超過國(guó)標(biāo)限定值（0.2 mg/kg；GB2762?2017），因而將來利用功能突變基因進(jìn)行品種低鎘化改造時(shí)應(yīng)該關(guān)注Pb含量的變化。Cr是人體必需的元素，由于本研究中籽粒本底Cr水平較低，我們認(rèn)為基因敲除引起的Cr含量增加并不足以對(duì)健康造成危害（中國(guó)營(yíng)養(yǎng)學(xué)會(huì)發(fā)布的2013版《中國(guó)居民膳食營(yíng)養(yǎng)素參考攝入量》中推薦我國(guó)成年人總鉻適宜攝入量為30 μg/d；GB2762?2017食品中污染物限量標(biāo)準(zhǔn)規(guī)定稻米鉻含量須低于1 mg/kg），但是在Cr污染地區(qū)，若本底Cr含量比較高，這種效應(yīng)就不能忽視，可能導(dǎo)致稻米Cr含量超標(biāo)。

本研究進(jìn)一步展示了基因敲除對(duì)水稻農(nóng)藝性狀的負(fù)面影響，發(fā)現(xiàn)不同品種中株高和產(chǎn)量都小幅度降低，這一結(jié)果與Ishikawa等[16]以及Tang等[30]的結(jié)果矛盾，而與Sasaki等[18]和Ishimaru等[17]的結(jié)果一定程度吻合，只是本研究中產(chǎn)量減少的程度比Sasaki等[18]報(bào)道的小得多。值得注意的是，我們和Tang等[30]的研究使用了相同的品種華占為試驗(yàn)材料，只是敲除靶位點(diǎn)不同，顯然敲除靶位點(diǎn)的差異不大可能是產(chǎn)生這種差異的原因，我們認(rèn)為可能是種植條件的不同導(dǎo)致了這種矛盾結(jié)果。在本研究中，檢測(cè)對(duì)農(nóng)藝性狀影響的植株種植在非Cd污染土壤（與Sasaki等[18]和Ishikawa等[16]的研究相同），而Tang等[30]的試驗(yàn)在嚴(yán)重Cd污染土壤上進(jìn)行，在鎘污染土壤種植時(shí)，一方面，由于Cd超量積累對(duì)水稻生長(zhǎng)不利[31]，野生型植株農(nóng)藝性狀可能受到負(fù)面影響，而突變體因鎘積累較低而得以避免，另一方面，突變體株系中存在基因敲除導(dǎo)致的生長(zhǎng)抑制效應(yīng)，而野生型不存在這種效應(yīng)，因而突變體與野生型中各自遭受的生長(zhǎng)抑制程度相似，導(dǎo)致兩者農(nóng)藝性狀表現(xiàn)無差別，而在一般土壤種植時(shí)，由于沒有Cd毒害，突變體與野生型對(duì)照的差別就只剩下基因敲除株系因功能失活而引起的生長(zhǎng)抑制，因此差異得以顯現(xiàn)。此外，其他種植條件的不同（如Mn的供應(yīng)）也可能是產(chǎn)生差異結(jié)果的原因。事實(shí)上，我們將華占和五山絲苗的各自兩個(gè)基因敲除株系在地點(diǎn)C種植時(shí)（正文提及的地點(diǎn)A、B之外的另一地點(diǎn)；該地稻米鎘含量只是輕微超標(biāo)），小區(qū)測(cè)產(chǎn)方產(chǎn)量測(cè)定結(jié)果顯示突變體與野生型對(duì)照幾乎無差別（數(shù)據(jù)未顯示）。當(dāng)然，以上猜想需要進(jìn)一步試驗(yàn)驗(yàn)證?？傊煌贩N及不同種植條件下，功能突變株系可能較野生型減產(chǎn)或無變化，提示生產(chǎn)中推廣通過基因敲除創(chuàng)制低鎘水稻品種時(shí)須謹(jǐn)慎，需要按照新品種的標(biāo)準(zhǔn)重新通過多年多點(diǎn)試驗(yàn)確定推廣地域。需要指出的是，對(duì)于Cd嚴(yán)重超標(biāo)稻米產(chǎn)區(qū)，一定程度的減產(chǎn)以換取稻米食用安全也是值得的。

本研究還揭示了當(dāng)基因敲除導(dǎo)致減產(chǎn)時(shí)，產(chǎn)量的降低主要由實(shí)粒數(shù)、結(jié)實(shí)率和千粒重降低所致，并為有效分蘗數(shù)的增加所部分補(bǔ)償，這是Sasaki等[18]和Ishikawa等[17]的研究中未曾提及的。吳平等[32]利用F2分離群體進(jìn)行分析的結(jié)果顯示，水稻株高和分蘗呈負(fù)相關(guān)關(guān)系。本研究發(fā)現(xiàn)基因敲除材料中伴隨株高小幅度降低分蘗略微增多，與前人報(bào)道吻合。

Allen等[33]通過在人體細(xì)胞中對(duì)40 000多個(gè)靶位點(diǎn)切割產(chǎn)生的大于1×109個(gè)突變序列進(jìn)行分析揭示了CRISPR/Cas9切割基因組DNA產(chǎn)生雙鏈斷裂后細(xì)胞修復(fù)產(chǎn)生突變的規(guī)律，發(fā)現(xiàn)突變結(jié)果與切割位點(diǎn)鄰近序列特征息息相關(guān)，并可依據(jù)靶位點(diǎn)序列信息準(zhǔn)確預(yù)測(cè)突變結(jié)果。我們通過該報(bào)道作者提供的網(wǎng)絡(luò)工具FORECasT（https://partslab.sanger. ac. uk/FORECasT）對(duì)本研究采用的基因靶位點(diǎn)突變結(jié)果進(jìn)行了預(yù)測(cè)，預(yù)測(cè)結(jié)果中概率最高的前三種突變類型與本研究的實(shí)際結(jié)果完全吻合（數(shù)據(jù)未顯示）。因此，本研究結(jié)果還驗(yàn)證了FORECasT在植物基因敲除突變結(jié)果預(yù)測(cè)中的準(zhǔn)確性。

總之，雖然水稻中功能缺失可能導(dǎo)致不同程度減產(chǎn)，但是到目前為止，若想通過改良品種的辦法使嚴(yán)重鎘超標(biāo)（如>1 mg/kg）稻米產(chǎn)區(qū)的大米鎘含量降至較低水平（如<0.2 mg/kg），除非通過轉(zhuǎn)基因過表達(dá)[13]、或基因片段[34]，否則利用功能缺失突變體仍然是唯一的選擇。而對(duì)于稻米Cd輕微污染（如含量低于0.5 mg/kg）產(chǎn)區(qū)低鎘品種的培育，更優(yōu)的選擇可能是利用[14, 35]、[15]和[36]等基因的突變體，因?yàn)檫@些基因功能缺失后能使鎘含量減少約50%的同時(shí)對(duì)其他元素含量和產(chǎn)量影響甚微。

謝辭:感謝江西省超級(jí)水稻研究發(fā)展中心蔡耀輝研究員、江西省農(nóng)業(yè)科學(xué)院土壤肥料與資源環(huán)境研究所徐昌旭研究員和謝杰博士在鎘污染耕地信息提供及材料種植中的幫助。

[1] Grant C A, Clarke J M, Duguid S, Chaney R L. Selection and breeding of plant cultivars to minimize cadmium accumulation., 2008, 390: 301-310.

[2] Sebastian A, Prasad M N V. Cadmium minimization in rice. A review., 2014, 34(1): 155-173.

[3] J?rup L, ?kesson A. Current status of cadmium as an environmental health problem., 2009, 238: 201-208.

[4] Arao T, Ae N. Genotypic variations in cadmium levels of rice grain., 2003, 49(4): 473-479

[5] Liu J, Zhu Q, Zhang Z, Xu J, Yang J, Wong M H. Variations in cadmium accumulation among rice cultivars and types and the selection of cultivars for reducing cadmium in the diet., 2005, 85: 147-153.

[6] 王林友, 竺朝娜, 王建軍, 張禮霞, 金慶生, 石春海. 水稻鎘、鉛、砷低含量基因型的篩選. 浙江農(nóng)業(yè)學(xué)報(bào), 2012, 24(1): 133-138.

Wang L Y, Zhu C N, Wang J J, Zhang L X, Jin Q S, Shi C H. Screening for rice (L.) genotyeps with lower Cd, Pb and As contents., 2012, 24(1): 133-138. (in Chinese with English abstract)

[7] 葉新新, 周艷麗, 孫波. 適于輕度Cd、As污染土壤種植的水稻品種篩選. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2012, 31(6): 1082-1088.

Ye X X, Zhou Y L, Sun B. Screening of suitable rice cultivars for the adaptation to lightly contaminated soil with Cd and As., 2012, 31(6): 1082-1088. (in Chinese with English abstract)

[8] Uraguchi S, Fujiwara T. Rice breaks ground for cadmium-free cereals., 2013, 16(3): 328-334.

[9] Uraguchi S, Fujiwara T. Cadmium transport and tolerance in rice: Perspectives for reducing grain cadmium accumulation., 2012, 5(1): 5.

[10] Satoh-Nagasawa N, Mori M, Nakazawa N, Kawamoto T, Nagato Y, Sakurai K, Takahashi H, Watanabe A, Akagi H. Mutations in rice () heavy metal ATPase 2 (OsHMA2) restrict the translocation of zinc and cadmium., 2012, 53(1): 213-224.

[11] Takahashi R, Ishimaru Y, Shimo H, Ogo Y, Senoura T, Nishizawa N K, Nakanishi H. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice., 2012, 35(11): 1948-1957.

[12] Yamaji N, Xia J, Mitani-Ueno N, Yokosho K, Feng M J. Preferential delivery of zinc to developing tissues in rice is mediated by P-type heavy metal ATPase OsHMA2., 2013, 162(2): 927-939.

[13] Ueno D, Yamaji N, Kono I, Huang C F, Ando T, Yano M, Ma J F. Gene limiting cadmium accumulation in rice., 2010, 107(38): 16 500-16 505.

[14] Uraguchi S, Kamiya T, Sakamoto T, Kasai K, Sato Y, Nagamura Y, Yoshida A, Kyozuka J, Ishikawa S, Fujiwara T. Low-affinity cation transporter (OsLCT1) regulates cadmium transport into rice grains., 2011, 108(52): 20 959-20 964.

[15] Shimo H, Ishimaru Y, An G, Yamakawa T, Nakanishi H, Nishizawa N K. Low cadmium (), a novel gene related to cadmium tolerance and accumulation in rice., 2011, 62(15): 5727-5734.

[16] Ishikawa S, Ishimaru Y, Igura M, Kuramata M, Abe T, Senoura T, Hase Y, Arao T, Nishizawa N K, Nakanishi H. Ion-beam irradiation, gene identification, and marker-assisted breeding in the development of low-cadmium rice., 2012, 109(47): 19 166-19 171.

[17] Ishimaru Y, Takahashi R, Bashir K, Shimo H, Senoura T, Sugimoto K, Ono K, Yano M, Ishikawa S, Arao T, Nakanishi H, Nishizawa N K. Characterizing the role of rice NRAMP5 in manganese, iron and cadmium transport., 2012, 2: 286.

[18] Sasaki A, Yamaji N, Yokosho K, Ma J F. Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice., 2012, 24(5): 2155-2167.

[19] Belhaj K, Chaparro-Garcia A, Kamoun S, Nekrasov V. Plant genome editing made easy: Targeted mutagenesis in model and crop plants using the CRISPR/Cas system., 2013, 9(1): 39.

[20] Peng H, Zhang Q, Li Y, Lei C, Zhai Y, Sun X, Sun D, Sun Y, Lu T. A putative leucine-rich repeat receptor kinase, OsBRR1, is involved in rice blast resistance., 2009, 230(2): 377-385.

[21] Miao J, Guo D, Zhang J, Huang Q, Qin G, Zhang X, Wan J, Gu H, Qu L J. Targeted mutagenesis in rice using CRISPR-Cas system., 2013, 23(10): 1233-1236.

[22] Feng Z, Zhang B, Ding W, Liu X, Yang D L, Wei P, Cao F, Zhu S, Zhang F, Mao Y, Zhu J K. Efficient genome editing in plants using a CRISPR/Cas system., 2013, 23(10): 1229-1232.

[23] Cong L, Ran F A, Cox D, Lin S, Barretto R, Habib N, Hsu P D, Wu X, Jiang W, Marraffini L A, Zhang F. Multiplex genome engineering using CRISPR/Cas systems., 2013, 339(6121): 819-823.

[24] Liu H, Ding Y, Zhou Y, Jin W, Xie K, Chen L L. CRISPR-P 2.0: An improved CRISPR-Cas9 tool for genome editing in plants., 2017, 10(3): 530-532.

[25] Lei Y, Lu L, Liu H Y, Li S, Xing F, Chen L L. CRISPR-P: A web tool for synthetic single-guide RNA design of CRISPR-system in plants., 2014, 7(9): 1494-1496.

[26] Holsters M, de Waele D, Depicker A, Messens E, van Montagu M, Schell J. Transfection and transformation of., 1978, 163(2): 181-187.

[27] Toki S, Hara N, Ono K, Onodera H, Tagiri A, Oka S, Tanaka H. Early infection of scutellum tissue with Agrobacterium allows high-speed transformation of rice., 2006, 47(6): 969-976.

[28] Allen G C, Flores-Vergara M A, Krasynanski S, Kumar S, Thompson W F. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethy- lammonium bromide., 2006, 1(5): 2320-2325.

[29] 劉巧泉, 陳秀花, 王興穩(wěn), 彭凌濤, 顧銘洪. 一種快速檢測(cè)轉(zhuǎn)基因水稻中潮霉素抗性的簡(jiǎn)易方法. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2001, 9(3): 264-268.

Liu Q Q, Chen X H, Wang X Y, Peng T L, Gu M H. A rapid simple method of assaying hygromysin resistance in transgenic rice., 2001, 9(3): 264-268. (in Chinese with English abstract)

[30] Tang L, Mao B, Li Y, Lv Q, Zhang L, Chen C, He H, Wang W, Zeng X, Shao Y, Pan Y, Hu Y, Peng Y, Fu X, Li H, Xia S, Zhao B. Knockout ofusing the CRISPR/Cas9 system produces low Cd-accumulating indica rice without compromising yield., 2017, 7(1): 14 438.

[31] Rizwan M, Ali S, Adrees M, Rizvi H, Zia-Ur-Rehman M, Hannan F, Qayyum M F, Hafeez F, Ok Y S. Cadmium stress in rice: Toxic effects, tolerance mechanisms, and management: A critical review., 2016, 23(18): 17 859-17 879.

[32] 吳平. 應(yīng)用RFLP標(biāo)記分析水稻株高與分蘗的遺傳相關(guān)性. 中國(guó)科學(xué)C輯: 生命科學(xué), 1996, 26(3): 264-270.

Wu P. Genetic correlation analysis of rice plant height and tiller number using RFLP markers., 1996, 26(3): 264-270. (in Chinese; the title was translated into English by us)

[33] Allen F, Crepaldi L, Alsinet C, Strong A J, Kleshchevnikov V, De Angeli P, Páleníková P, Khodak A, Kiselev V, Kosicki M, Bassett A R, Harding H, Galanty Y, Mu?oz-Martínez F, Metzakopian E, Jackson S P, Parts L. Predicting the mutations generated by repair of Cas9-induced double-strand breaks., 2018, 37: 64-72.

[34] Wang C, Guo W, Ye S, Wei P, Ow D W. Reduction of Cd in rice through expression of-like gene fragments., 2016, 9(2): 301-304

[35] Uraguchi S, Kamiya T, Clemens S, Fujiwara T. Characterization of OsLCT1, a cadmium transporter from indica rice ()., 2014, 151(3): 339-347.

[36] Hao X, Zeng M, Wang J, Zeng Z, Dai J, Xie Z, Yang Y, Tian L, Chen L, Li D. A node-expressed transporter OsCCX2 is involved in grain cadmium accumulation of rice., 2018, 9: 476.

Creation of Low-Cd-accumulatingRice by Disruption ofGene via CRISPR/Cas9

LONG Qizhang, HUANG Yonglan, TANG Xiuying, WANG Huimin, LU Ming, YUAN Linfeng, WAN Jianlin*

[,,;Corresponding author,:]

【】 To accelerate the breeding process of low-Cd-accumulating rice,【】the widely appliedpure-line cultivars, Wushansimiao and Zhongzao 35, and the popular parents of hybrids, Huazhan and Wufeng B were used as materials to createknockout lines via CRISPR/Cas9. The marker-freemutants together with the wild types were grown either in a Cd-contaminated field (Field A) or in a non-Cd-contaminated field (Field B) for subsequent experiments. The changes of the grain Cd content were investigated for plants grown in both Field A and Field B. The changes of the contents of other related minerals were also examined but only for those plants grown in Field A. The agronomic changes of the mutants were further assessed by examining those plants grown in Field B. 【】 Rice lines with mutations atwere successfully generated for all the four cultivars. On one hand, the grain Cd level inmutants grown in Field B was reduced by 85.5% on average (<0.02 mg/kg); as a result ofmutation, the Cdlevel in plants grown in Field A, was reduced by 94.8% (<0.1 mg/kg), meanwhile, the Mn level was decreased by 52.7%; on the contrary, the Cr level was elevated by 59.5% on average for the four cultivars and an increase of 79.1% in the Pb level was also observed but only inmutants of Huazhan; the levels of other mineral elements including Fe, Cu , Zn, Ca, Se and As were slightly or hardly affected. On the other hand, a slight decrease in plant height and an average yield loss of 6.9% were observed due tomutation for all cultivars; the yield loss could be ascribed to a reduction in seed setting rate and seed weight and was compensated by a slight increase in effective tiller number. 【】 Low-Cd-accumulating rice plants can be generated bydisruption at the cost of an acceptable yield loss and the low-Cd-accumulating rice lines produced in the present study may serve as promising cultivars to eliminate rice Cd contamination in many areas.

rice; cadmium; genetic editing; breeding

Q755; S511.034

A

1001-7216(2019)05-0407-14

10.16819/j.1001-7216.2019.9026

2019-03-05；

2019-05-24。

國(guó)家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目（2016YFD0102102）；江西現(xiàn)代農(nóng)業(yè)科研協(xié)同創(chuàng)新專項(xiàng)（JXXTCX2015001-003）；江西省科技計(jì)劃項(xiàng)目農(nóng)業(yè)領(lǐng)域重點(diǎn)項(xiàng)目（20141BBF60001）。