秦 娜 朱燦燦 代書桃 宋迎輝 李君霞 王春義

（河南省農(nóng)業(yè)科學院糧食作物研究所，450002，河南鄭州）

葉片是作物進行光合作用的主要器官，葉色與葉綠素含量、光合作用強度及作物產(chǎn)量密切相關(guān)[1-4]，谷子因其基因組小、高效轉(zhuǎn)化、生育期短、抗旱耐瘠，且易于在實驗室培養(yǎng)等優(yōu)點，正在快速發(fā)展成為禾本科黍亞科功能基因組和C4光合作用的模式植物。因此，谷子葉色突變體已成為解析C4作物高效光能利用效率和探索葉綠體發(fā)育機制的重要遺傳材料[5-7]。作物葉色突變基因的鑒定與功能分析，對揭示葉綠素代謝分子機制、創(chuàng)制高光效育種新種質(zhì)奠定技術(shù)基礎(chǔ)，同時作為標記性狀在農(nóng)作物雜交種選育方面發(fā)揮重要作用[8]。

葉色突變體在水稻、玉米等作物中均有較多報道。Wang等[9]研究發(fā)現(xiàn)，水稻中200多個葉色突變體中已克隆出50多個突變基因；Guan等[10]對200多個玉米葉色突變體進行相關(guān)研究，并成功克隆出黃葉色突變基因ygl1-1。目前，谷子相關(guān)葉色突變體則鮮有報道，李傳宗[11]鑒定出苗期黃葉色突變基因Seita.7G305300，該基因編碼磷酰乙醇胺胞苷轉(zhuǎn)移酶（phosphorylethanolamine cytidylyltransferase，PECT），通過線粒體途徑影響植物葉綠體發(fā)育，并最終導致苗期黃葉色表型；Li等[12]克隆出了黃綠色葉片突變基因SiYGL1，該基因編碼葉綠素合成途徑關(guān)鍵酶D亞基（CHLD）蛋白，其突變后抑制葉綠素合成，導致葉色突變；條紋葉突變體則是由SiSTL1[13]和SiSTL2[14]突變致使葉綠體發(fā)育異常，2個基因分別編碼核糖核苷酸還原酶（ribonucleotide reductase，RNR）大亞基蛋白和脫氧胞苷單磷酸脫氨酶（deoxy-cytidine monophosphate deaminase，DCD）蛋白。目前，隨著谷子成為模式植物的研究熱點，其葉色突變相關(guān)分子機制研究不能滿足對C4光合作用和葉綠體發(fā)育機制的研究，因此迫切需要創(chuàng)制更多的葉色突變體并發(fā)掘其作用基因和功能。

隨著測序技術(shù)的快速發(fā)展，Abe等[15]利用突變位點圖譜（MutMap）技術(shù)對混池進行測序和克隆突變體基因，該技術(shù)被稱為一種克隆突變體基因的快速高效技術(shù)，已在擬南芥[16]、水稻[17]、玉米[18]和谷子[19]等作物中成功克隆多個突變基因。MutMap適合對甲基磺酸乙酯（EMS）誘變的隱性突變基因進行鑒定與克隆。植物通過EMS誘變和自交得到純合體后，將突變體和親本雜交構(gòu)建F2分離群體，分別得到野生型表型和突變體表型群體，對這2個群體DNA分別進行等量混合，得到野生型DNA混池和突變體DNA混池，對2個混池分別進行DNA測序，利用MutMap pepline進行數(shù)據(jù)分析，計算SNP在突變體混池和野生型混池出現(xiàn)的頻率，即可得到和突變表型連鎖的染色體區(qū)段和可能的突變位點，從而篩選鑒定出突變基因。

在課題組前期研究工作中，利用EMS對谷子主推品種豫谷1號進行誘變，篩選鑒定到1個黃葉色突變體ylm-1，且該突變體表現(xiàn)出較高的光合速率和籽粒產(chǎn)量[20]。為進一步深入研究該黃葉色突變體的遺傳機制，本研究以黃葉色突變體ylm-1和野生型豫谷1號雜交構(gòu)建的F2分離群體為定位群體，采用MutMap技術(shù)對突變基因ylm-1進行精細定位，開發(fā)新的分子標記，鑒定和克隆候選基因，揭示ylm-1基因在葉綠素代謝中的功能，為解析C4植物葉綠體發(fā)育機制奠定理論和技術(shù)基礎(chǔ)。

1 材料與方法

1.1 供試材料

供試材料分別為野生型親本豫谷1號（安陽市農(nóng)業(yè)科學院選育）、黃葉色突變體ylm-1及ylm-1與野生型親本雜交構(gòu)建的F2分離群體（河南省農(nóng)業(yè)科學院糧食作物研究所創(chuàng)制）。檢測突變體遺傳背景所用材料分別為冀谷41（河北省農(nóng)林科學院谷子研究所選育）、金苗紅酒谷（山西省農(nóng)業(yè)科學院谷子研究所選育）和豫谷17和豫谷28（河南省農(nóng)業(yè)科學院糧食作物研究所選育）。于2019和2020年種植在河南省現(xiàn)代農(nóng)業(yè)研發(fā)基地（河南省新鄉(xiāng)市）。

1.2 試驗方法

1.2.1 突變體ylm-1表型及遺傳背景分析 于谷子苗期分別對突變體ylm-1和野生型葉片顏色進行調(diào)查，同時選用已有的50對SSR標記（表1）對ylm-1和野生型及4個遺傳背景差異的品種豫谷17、冀谷41、金苗紅酒谷和豫谷28進行分子標記檢測，進一步驗證該突變體的遺傳背景。突變體ylm-1和野生型豫谷1號雜交獲得F1，F(xiàn)1自交獲得F2分離群體，分別于3～5葉期調(diào)查F1葉色情況和F2葉色分離情況。計算F2分離群體中黃葉色表型和野生表型的分離比，并進行χ2檢驗，分析突變體ylm-1的遺傳規(guī)律。

1.2.2 MutMap技術(shù)定位候選基因 采用MutMap技術(shù)和DNA混池測序（bulked segregant analysis，BSA）法定位和克隆候選基因[17-18]。于5葉期從F2分離群體中選取黃葉色表型植株30株葉片，提取DNA后進行等量混合，構(gòu)建突變型DNA池，野生型DNA池為對照，進行全基因組測序（百邁克生物科技有限公司）。完成測序后，將ylm-1突變型混池數(shù)據(jù)與野生型數(shù)據(jù)分別同谷子參考基因組豫谷1號比對，獲得突變型DNA池和野生型基因組的一致性序列。并對染色體上得到SNPs數(shù)據(jù)集和SNP位點進行過濾，去除豫谷1號和混池DNA共同的SNPs，篩選出混池DNA特異性的SNPs。最后計算SNP-index，將SNP-index數(shù)值接近1的區(qū)域定為關(guān)聯(lián)區(qū)間，根據(jù)關(guān)聯(lián)區(qū)間內(nèi)基因注釋結(jié)果進一步鑒定候選基因。

利用Phytozome和Plant GDB數(shù)據(jù)庫查找定位區(qū)間內(nèi)的基因，對區(qū)間內(nèi)基因進行功能注釋。同時結(jié)合SIFGD數(shù)據(jù)庫進行候選基因的轉(zhuǎn)錄模式分析[21]。

1.2.3 dCAPS分子標記開發(fā)與檢測 根據(jù)候選SNPs區(qū)域內(nèi)的DNA序列和突變位點，采用軟件Primer 5.0設(shè)計dCAPS分子標記（表1），分別以ylm-1和野生型基因組DNA及F2分離群體為模板，檢測ylm-1和野生型基因型差異，鑒定黃葉色表型是由該突變位點導致。

表1 引物名稱和序列信息Table 1 Name and sequence of primers

PCR擴增總反應(yīng)體系：模板（10ng）2.00μL，10×Buffer 2.00μL，dNTPs（2.50mmol/L）1.60μL，正反向特異性引物（10μmol/L）各 0.50μL，Taq DNA聚合酶（0.50U）0.10μL，加 ddH2O 至 20.00μL。PCR擴增程序為94℃5min；94℃40s，57℃40s，72℃30s，35個循環(huán)；72℃5min，4℃保存。

PCR產(chǎn)物酶切反應(yīng)體系：PCR產(chǎn)物6.00μL，限制性內(nèi)切酶10U，10×Buffer 1.50μL，加ddH2O至16.00μL，37℃酶切12h。取酶切產(chǎn)物8.00μL于5%瓊脂糖凝膠電泳檢測。

2 結(jié)果與分析

2.1 突變體ylm-1表型及遺傳學分析

突變體ylm-1與野生型苗期葉色表現(xiàn)出明顯的差異（圖1），隨著植株的生長，開花期葉片呈現(xiàn)黃綠色（圖2），整個生長發(fā)育過程未恢復(fù)為正常葉色，而野生型全生育期葉片均表現(xiàn)為正常綠色。在苗期和開花期，F(xiàn)1葉色表型與野生型一致。

圖1 苗期野生型（WT）、突變體ylm-1和F1葉片顏色Fig.1 Leaf color of the wild type(WT),mutant ylm-1 and F1at seedling stage

圖2 開花期野生型（WT）、突變體ylm-1和F1葉片顏色Fig.2 Leaf color of the wild type（WT）,mutant ylm-1 and F1at flowering stage

在110株BC1F2分離群體中，通過谷子苗期葉色表型調(diào)查，正常葉色植株89株，黃葉色植株21株，卡方檢驗（χ2=0.54，P=0.38），符合3:1分離比，表明該黃葉色突變性狀由1對單隱性基因控制。

選用50對分子標記對ylm-1、野生型及豫谷17、冀谷41、豫谷28和金苗紅酒谷基因型進行檢測，突變體與野生型對照的條帶在凝膠上的相對位置一致，DNA片段大小相同，說明突變體與野生型遺傳背景相同，不是其他品種混雜。而其余4個谷子品種的條帶與野生型在部分標記間顯示出差異性（圖3）。

圖3 突變體ylm-1遺傳背景檢測Fig.3 Genetic background detection of mutant ylm-1

2.2 基于MutMap技術(shù)的候選基因鑒定

為進一步探究ylm-1黃葉色形成的分子機制，利用MutMap技術(shù)對突變體ylm-1的候選基因進行鑒定[21-22]。將ylm-1與野生型豫谷1號雜交的F2群體中30株黃葉色表型植株進行混池測序，野生型豫谷1號作為對照，根據(jù)測序結(jié)果進行SNP位點分析。根據(jù)參考基因組信息，開發(fā)并過濾SNP位點，計算親本與混池中每個SNP的SNP-index值，其中SNP-index值接近1的為候選SNPs。最終得到高質(zhì)量可信的突變體ylm-1混池特異SNP 158 149個，突變體混池SNP-index分布見圖4。

圖4 SNP-index值在染色體上的分布Fig.4 Distribution map of SNP-index values on chromosomes

取所有位點擬合值的99百分位數(shù)，計算得0.97，最終得到3個區(qū)域，總長度0.45Mb，共包含26個基因（表2）。據(jù)統(tǒng)計，親本與混池間存在非同義突變的SNP共32個，分別位于13個基因上，即非同義突變的基因共13個，以上非同義突變基因與黃葉色突變性狀可能直接相關(guān)（表3）。

應(yīng)用BLAST軟件對候選區(qū)間內(nèi)13個非同義突變編碼基因進行多個數(shù)據(jù)庫（NR，Phytozome，PlantGDB，Swiss-Pro）的深度注釋，其中Seita.9G249900是一個編碼紅葉綠素代謝還原酶（RCCR）的相關(guān)基因（表3），該基因位于第9號染色體19 937 988與19 940 620bp之間。利用SIFGD數(shù)據(jù)庫對13個非同義突變基因在谷子不同組織器官的表達模式進行分析[23]，基于mRNA-seq FPKM（fragments per kilobase million）的轉(zhuǎn)錄水平，結(jié)果表明Seita.9G249900在葉片的表達量顯著高于其他部位（圖5）。因此，推測該基因為ylm-1的重要候選基因。

圖5 13個非同義突變基因在不同部位中的FPKM值Fig.5 FPKM values of thirteen non-synonymous mutation genes in different parts

2.3 突變位點驗證

為進一步驗證Seita.9G249900為黃葉色突變體ylm-1的候選基因，根據(jù)該基因突變位點序列設(shè)計了dCAPS標記（突變位點造成限制性內(nèi)切酶AluI的差異）（表1）。并利用該分子標記檢測野生型豫谷1號、突變體ylm-1及野生型×ylm-1 F2單株群體帶型，凝膠電泳結(jié)果（圖6）顯示，突變體ylm-1與野生型豫谷1號存在明顯差異，野生型×ylm-1 F2群體中Seita.9G249900基因型和表型呈共分離。Seita.9G249900位點為純合野生型或者雜合型時，葉片顏色表現(xiàn)正常；Seita.9G249900位點為純合突變型時，葉片為黃色。綜上所述，ylm-1中確實存在基因位點突變，其黃葉色表型很可能是由基因Seita.9G249900的突變導致。

圖6 野生型×ylm-1 F2分離群體中dCAPS標記凝膠分析Fig.6 Analysis of electrophoretic images of dCAPS marker for wild type×ylm-1 of F2populations

2.4 基因位點的功能分析

結(jié)合NCBI和Phytozome數(shù)據(jù)庫基因組分析結(jié)果發(fā)現(xiàn)，Seita.9G249900含有2個外顯子和1個內(nèi)含子（圖7a）。其中ylm-1突變發(fā)生在第2外顯子361bp處（A/T），導致1個谷氨酸（E）變成天冬氨酸（D）（圖7b）。Seita.9G249900編碼PCCR，主要參與葉綠素的分解代謝、生物合成及光合系統(tǒng)Ⅱ相關(guān)光復(fù)合物的分解代謝過程，與擬南芥AT4G37000基因具有同源性，其編碼的紅葉綠素代謝還原酶（RCCR）在分解葉綠素的卟啉成分中起著重要的催化作用；同時也是水稻LOC-OS10g0389100的同源基因，且編碼葉綠體中RCCR，研究表明RCCR1基因沉默或表達下調(diào)時，植株葉片呈現(xiàn)黃色且生育期增加，葉片色素含量顯著降低。綜上所述，結(jié)合候選基因的功能注釋與位點分析，ylm-1應(yīng)為谷子1個新的黃葉色突變體，初步推斷其黃葉色表型是由一個編碼RCCR基因Seita.9G249900突變導致，其基因功能與作用機制需進一步驗證與分析。

圖7 Seita.9G249900編碼區(qū)結(jié)構(gòu)圖（a）及氨基酸序列分析（b）Fig.7 The structure diagram of Seita.9G249900 coding region(a)and analysis of peptide sequence(b)

3 討論

3.1 MutMap技術(shù)在突變基因定位與克隆中的應(yīng)用

MutMap是Abe等[15]于2012年提出的一種新的基因克隆技術(shù)，較圖位克隆更快捷、更高效、成本更低，較適合于由EMS誘變產(chǎn)生的隱性突變基因分析。植物通過EMS誘變和自交得到純合體后，將其與親本雜交得到F1，F(xiàn)1自交得到的F2群體會出現(xiàn)野生型和突變體的分離。對這2個表型群體的DNA分別進行等量混合，構(gòu)建野生型DNA混池和突變體DNA混池。將2個混池分別進行DNA測序，利用MutMap pepline進行數(shù)據(jù)分析，計算SNP在突變體混池和野生型混池出現(xiàn)的頻率，即可得到和突變表型連鎖的染色體區(qū)段和可能的突變位點。另外MutMap技術(shù)在定位質(zhì)量性狀隱性基因和數(shù)量性狀QTL位點方面均有較廣泛和高效的應(yīng)用[24-26]。本研究以野生型豫谷1號和黃葉色突變體ylm-1雜交構(gòu)建的F2群體為定位群體，利用MutMap技術(shù)對ylm-1突變基因進行精細定位，通過提取30株黃葉色表型植株葉片DNA并等量混合作為突變型DNA混池，野生型豫谷1號DNA作為對照，二者進行基因組重測序，根據(jù)測序和分析結(jié)果篩選SNP-index接近1的染色體區(qū)域作為與目標性狀關(guān)聯(lián)的候選區(qū)域，結(jié)果顯示該區(qū)間內(nèi)共包含13個非同義突變基因，根據(jù)基因功能注釋結(jié)果和表達模式分析結(jié)果，進一步推測候選基因。

3.2 紅葉綠素代謝還原酶（RCCR）的功能分析

植物葉綠素代謝和信號調(diào)控途徑較為復(fù)雜，作用基因眾多，其中任何一個基因發(fā)生突變都有可能引起葉綠素代謝途徑的變化，從而表現(xiàn)為葉色異常甚至植株死亡的現(xiàn)象。RCCR是葉綠體中一種可溶性蛋白，屬于鐵氧還蛋白依賴的膽紅素還原酶家族成員之一[27]，可以將紅葉綠素代謝產(chǎn)物（RCC）的C20/C1雙鍵還原形成初級熒光葉綠素分解代謝產(chǎn)物（pFCC）[28]。葉綠素代謝相關(guān)基因的研究對解析葉綠素代謝、葉色異常及光合效率變化的機制具有重要的理論價值[29]。

隨著分子生物學的發(fā)展，有關(guān)葉綠素代謝相關(guān)基因的研究不斷深入，并取得一定的突破性進展。目前研究普遍認為葉綠素的代謝存在2種途徑，即葉綠素酶（CLH）代謝途徑和脫鎂葉綠素酶（PPH）途徑[30]。在由脫鎂葉綠酸a（pheide a）轉(zhuǎn)化成pFCC的過程中，RCCR在葉綠素代謝中起重要的調(diào)控作用，其活性與葉綠素代謝呈明顯的正相關(guān)，RCCR發(fā)揮活性的過程對葉片保持綠色起著重要調(diào)控作用[31-33]。Li等[34]研究發(fā)現(xiàn)RCCR基因表達與活性在葉片衰老過程中大幅度下調(diào)。Pruzinská等[35]報道了擬南芥中缺乏RCCR的acd2-2突變體表現(xiàn)出葉色變淺、葉綠素含量降低的表型。Fukasawa等[36]研究指出，BoRCCR表達量在青花菜黃化葉片中顯著降低。Gomez-Lobato等[37]研究發(fā)現(xiàn)細胞分裂素、1-MCP以及熱空氣處理均能增加BoRCCR的表達量。

本研究中，利用MutMap技術(shù)克隆到一個與紅葉綠素代謝相關(guān)的RCCR基因Seita.9G249900，其是擬南芥ACD2（AT4G37000）基因的同源基因，編碼的RCCR在分解葉綠素的卟啉成分中起著重要的催化作用；同時也是水稻RCCR1（LOC-OS 10g0389100）基因的同源基因，在RCCR1基因沉默或表達下調(diào)時，植株葉片呈現(xiàn)黃色且生育期增長，葉片色素含量顯著降低。根據(jù)候選基因功能注釋和表達分析，結(jié)果表明編碼RCCR蛋白的基因Seita.9G249900發(fā)生突變可能是導致ylm-1葉色發(fā)生突變的原因，該基因的功能與作用機制尚需進一步驗證與分析。

4 結(jié)論

黃葉色突變性狀受1對隱性核基因控制，定位于第9號染色體，利用MutMap技術(shù)在關(guān)聯(lián)區(qū)間內(nèi)精細定位到13個非同義突變基因，根據(jù)候選基因功能注釋和表達分析，鑒定和克隆到一個與紅葉綠素代謝相關(guān)的RCCR基因Seita.9G249900，該基因為導致ylm-1葉色發(fā)生突變的候選基因，其物理位置位于第9號染色體19 937 988與19 940 620bp之間，推測Seita.9G249900基因在葉綠素代謝和光合作用過程中起重要作用。