小麥籽粒蛋白含量相關基因NAM新等位變異的挖掘

朱 娉，馮 波，徐智斌，樊小莉，王 濤

(1.中國科學院成都生物研究所，四川成都 610041； 2.中國科學院大學，北京100049)

?

小麥籽粒蛋白含量相關基因NAM新等位變異的挖掘

朱 娉1,2，馮 波1，徐智斌1，樊小莉1，王 濤1

(1.中國科學院成都生物研究所，四川成都 610041； 2.中國科學院大學，北京100049)

NAM基因是NAC基因家族的成員之一，該基因能加速植株葉片衰老，增加籽粒蛋白質積累，并促進鋅、鐵等微量元素吸收。為選育高蛋白和高微量元素含量的小麥品種，用NAM共顯性標記Xuhw89對858份小麥材料進行了篩選，并進一步擴增全長序列。結果表明，有3個栽培二粒小麥含有功能性NAM基因，其中，F37與F40來源的NAM基因分別與已公布的 TaNAC2基因和 TtNAM-A1基因(GenBank登錄號分別為HQ872050和DQ869672)相同；F34來源的NAM基因(命名為TdNAM，GenBank登錄號為KU529317)與已公布的 TtNAM-B2基因(GenBank登錄號為DQ869676)相似度最高，共有5處核苷酸差異，1處發生在內含子，另外4處發生在外顯子區，導致了3個氨基酸殘基的變化，預測分析發現這種差異不會對蛋白功能產生影響。因此，本研究篩選獲得的NAM基因新等位變異可為進一步選育高蛋白、高微量元素含量的小麥品種提供材料和基因資源。

小麥；NAM基因；新等位變異；分子標記；序列分析

谷物是人類攝取蛋白質和微量元素的主要來源。世界衛生組織報告顯示，超過20億的人口缺乏鋅、鐵等主要微量元素，超過1/3的世界人口處于蛋白質缺乏的亞健康狀態[1-2]。小麥作為世界主要的糧食作物之一，提高其蛋白質及鋅、鐵等微量元素含量對改善人類健康至關重要。

前人的研究表明，控制籽粒蛋白質含量(Grain protein contents,GPC)的數量性狀位點 Gpc-B1位于野生二粒小麥(Triticumturgidumvar.dicoccoides)6BS染色體上，該位點具有多效性，在開花期，增加可溶性蛋白和氨基酸在旗葉中的積累，提高氮同化效率，并最終影響籽粒大小與蛋白質含量[3]。Uauy等[4]進一步將 Gpc-B1位點上控制籽粒蛋白質含量的關鍵基因圖位克隆，命名為 NAM-B1，并采用RNA干擾技術降低 NAM-B1基因表達水平，植株衰老時間延遲3周，籽粒中的蛋白質、鋅、鐵元素含量下降30%。然而，六倍體小麥中 Gpc-B1位點上的 NAM-B1基因中1個堿基的插入導致移碼框突變，NAC結構域丟失，蛋白無功能，籽粒蛋白及鋅、鐵微量元素含量明顯偏低[4]。因此，篩選具有功能性NAM基因的材料是小麥品質改良的另一途徑。

NAM轉錄因子是NAC轉錄因子家族成員之一。此類轉錄蛋白的N端為高度保守的NAC結構域(含A、B、C、D和E五個亞基)，含有蛋白質和DNA結合位點；C端為轉錄激活域，具有高度多樣性[5-8]。Jamar等[9]發現不同大麥品種的 HvNAM-1多態性明顯，其編碼蛋白的C亞基中脯氨酸代替丙氨酸可能影響蛋白質構像或轉錄過程中的DNA結合能力[10]，結合籽粒蛋白含量的測定發現，新等位變異的品種籽粒蛋白含量偏高。Wang等[11]在大麥中發現HvNAM- 1基因編碼區544位堿基的變異(脯氨酸→丙氨酸)，造成籽粒蛋白含量降低。可見，NAM基因的多態性會影響籽粒性狀。然而，目前小麥中還沒有報道影響籽粒蛋白質含量的NAM等位變異。

Mesfin等[12]和Khan等[13]將野生二粒小麥 Gpc-B1基因轉化進入面包小麥中，獲得了Glupro穩定高蛋白品種。隨后，Khan等[13]開發出限制性長度多態性標記用于篩選 Gpc-B1位點，然而這些標記的PCR產物需要酶切處理才能判斷基因型。Distelfeld等[14]開發的共顯性標記Xuhw89為特異性標記，可直接檢測到低蛋白含量的四倍體硬粒小麥LDN和高蛋白含量的野生二粒小麥DIC的 Gpc-B1位點中4 bp大小的差異，從而高效篩選出功能性NAM基因。因此，本研究利用標記Xuhw89對858份小麥材料進行篩選，以期檢測含有 Gpc-B1位點的籽粒高蛋白含量品種，并通過同源克隆獲得NAM基因，以期獲得新的等位變異，為進一步培育高蛋白和高微量元素的小麥新品種提供種質材料。

1 材料與方法

1.1 材料與試劑

本研究所用材料包括含 Gpc-B1基因的高蛋白小麥品種Glupro(陽性對照)、LDN(陰性對照)及供試的858份小麥材料(人工六倍體小麥345份、四倍體小麥100份、中國小麥微核心種質222份及地方品種小麥191份)。Glupro由美國農業部提供，人工六倍體材料由SIMMTY提供，四倍體由中國科學院遺傳與發育生物學研究所王道文課題組提供，微核心種質由中國農業科學院賈繼增課題組提供，地方品種來自本實驗室。

2×TaqPCR MasterMix和DNA純化回收試劑盒購自天根生化科技有限公司(北京)；KOD-Plus-Neo高保真酶、10×Buffer、dNTPs和MgSO4均購自TaKaRa公司(大連)；pEASY-Blunt Cloning kit、Trans-T1 Phage Resistant Chemically Competent Cell購自北京全式金生物技術有限公司；引物由成都擎科梓熙生物有限公司合成。

1.2 Xuhw89標記的檢測

采用CTAB方法提取基因組DNA[19]。利用共顯性標記Xuhw89[14](F:5′-TCTCCAAG AGGGGAGAGACA-3′; R:5′-TTCCTCTCCC ATGAATCTAGCA-3′)對陽性對照、陰性對照及供試的858份小麥材料DNA進行PCR擴增。擴增體系(10 μL)：2×TaqPCR Master Mix 5 μL，上下游引物各0.1 μL，DNA模板(100 ng·μL-1)1 μL，加ddH2O至10 μL。擴增程序：95 ℃預變性4 min；94 ℃變性30 s，57 ℃退火30 s，72 ℃延伸30 s，32個循環；72 ℃延伸10 min。PCR產物經6%聚丙烯酰胺凝膠電泳分離，銀染顯影，拍照分析。

1.3 NAM基因的同源克隆

以1.2中篩選到的與陽性對照Glupro擴增到的片段大小一致的材料DNA為模板，采用用于擴增山羊草NAM基因的引物P1/P2[5](P1：5′-AGATCTGATGAGGTCCATGGG-3′；P2：5′-ATGCCCGTATGTGGTGTTCATAT-3′)進行NAM基因全長擴增。擴增體系(20 μL)：10×Buffer 2 μL，2 mmol·L-1dNTPs 2 μL，25 mmol·L-1MgSO41 μL，上下游引物各0.2 μL，KOD-plus-Neo酶0.4 μL，100 ng·μL-1模板DNA 1 μL，加ddH2O至20 μL。擴增程序：95 ℃預變性4 min；95 ℃變性40 s，55.3 ℃退火30 s，68 ℃延伸1 min 20 s，32個循環；72 ℃延伸10 min。PCR產物用1%瓊脂糖凝膠電泳檢測。將目的片段回收后克隆，每個樣品3次重復，最后交由北京諾賽基因組研究中心有限公司進行DNA測序。

1.4 NAM基因的生物信息學分析

用DNAMAN進行序列比對分析；用BLASTn(http://blast.ncbi.nlm.nih.gov/)檢索同源序列；用New GENSCAN(http://genes.mit.edu/GENSCAN.html)進行預測分析基因結構及氨基酸組成；用NCBI的CD-search(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)進行結構域預測；用PredictProtein(https://www.predictprotein.org/) 預測蛋白功能位點。

2 結果與分析

2.1 Xuhw89標記的檢測結果

利用共顯性分子標記Xuhw89對供試的858份小麥材料進行篩選，結果(圖1)表明，有3個栽培二粒小麥材料F34、F37和F40的擴增帶型與Glupro一致，獲得大小為122 bp的擴增條帶，其余材料帶型與LDN一致，獲得大小為126 bp的擴增條帶。

2.2 3個二粒小麥NAM基因的克隆及其核酸序列分析

利用引物P1/P2對3個栽培二粒小麥材料 F34、F37和F40的基因組DNA進行PCR擴增，得到大小約1.5 kb的擴增產物(圖2)。測序后分析發現，這3個材料的NAM基因序列均含有一個完整的開放閱讀框。NCBI在線BLASTn表明，F37來源的NAM基因序列與已公布的 TaNAC2基因序列(GenBank登錄號為HQ872050)一致；F40來源的NAM基因序列與已公布的 TtNAM-A1基因序列(GenBank登錄號為DQ869672)一致；F34來源的NAM基因序列與已公布的 TtNAM-B2基因序列(GenBank登錄號為DQ869676)相似度最高。因此，推測F34來源的NAM基因序列為小麥NAM基因的一個新等位變異類型，將其暫命名TdNAM，同時將其提交至GenBank，得到序列登錄號為KU529317。

M:pUC18 DNA/MspⅠ Marker; G:Glupro; 1:F34; 2:F37; 3:F40; L:LDN; 4:F38; 5:F42; 6:S64; 7:S81; 8:h45.

圖1 標記Xuhw89擴增產物的聚丙烯酰胺凝膠電泳圖

Fig.1 Polyacrylamide gel of the amplified

products of molecular marker Xuhw89

進一步對TdNAM基因核酸序列進行分析發現，TdNAM基因序列包括3個外顯子和2個內含子。與其相似度最高的 TtNAM-B2基因序列相比發現，TdNAM基因核酸序列共有5處核苷酸變異，其中，1處發生在第一個內含子第311位(C→T)；另外4處發生在外顯子區，分別為第一個外顯子第181位(A→G)及第三個外顯子第897位(C→G)、第1 027位(G→A)和第1 322位(G→C)(圖3)。

M:DL5000 DNA Marker; G:Glupro; 1:F34; 2:F37; 3:F40.

圖3 TtNAM-B2基因與 TdNAM基因的結構比對圖

A～E:NAM亞結構域；△表示蛋白結合位點；*表示核苷酸結合位點。

A-E:NAM sub-domains; △ indicates the protein-binding region; * indicates the polynucleotide-binding region.

圖4NAM基因編碼氨基酸序列的比對分析圖

Fig.4 Sequence analysis ofNAMtranscription factors

2.3 推導的蛋白質序列分析

TaNAC2基因和 TtNAM-A1基因編碼405個氨基酸殘基，TdNAM基因編碼396個氨基酸殘基。通過序列比對分析發現，其蛋白序列的N端均有一個典型的NAC結構域，A、B、C、D和E五個亞基高度保守；C端為轉錄激活域。野生二粒小麥中的NAM等位基因在功能區域均高度保守，以 TtNAM-B2(GenBank登錄號為ABI94355)作為參考序列，TdNAM有3處氨基酸殘基差異，包括位于B亞基的第61位(Thr→Ala)及位于TAR結構域的第197位(Ala→Gly)和第339位(Gly→Arg)。通過PredictProtein軟件分析發現，本研究獲得的TdNAM基因編碼的氨基酸殘基變異沒有發生在轉錄因子的核酸結合位點及蛋白結合位點，且通過SNAP2預測，該3個堿基差異不會對蛋白功能產生影響。

3 討 論

NAM轉錄因子作為NAC家族成員，能調控植物的生長及籽粒品質性狀。野生型 NAM-B1在栽培條件下容易發生突變，因此培育品種中含有功能性的 NAM-B1基因較少[15]。Asplund等[20]在63份瑞典小麥中僅篩選到4份六倍體材料含有野生型 NAM-B1。Hagenblad等[15]利用SSR標記從367份來自世界各地的面包小麥品種中僅篩選到5份小麥材料(4份春小麥和1份冬小麥)含有野生型 NAM-B1基因。Hu等[16]研究表明，能擴增出Xuhw89特征條帶的小麥材料籽粒蛋白含量明顯高于其他品種。本研究通過共顯性標記Xuhw89在858份小麥材料中篩選，得到3個栽培二粒小麥含有控制籽粒高蛋白含量的NAM基因，可作為籽粒高蛋白含量候選材料進一步研究。

Ooka等[6]發現NAC轉錄因子的A、C和D亞基高度保守，與轉錄因子結合核苷酸或蛋白質相關； B、E亞基則與生長發育或組織特異表達相關。本研究獲得的3個NAM基因均含有3個外顯子和2個內含子，與山羊草中NAM基因結構一致[5]。通過氨基酸序列分析發現，3個NAM轉錄因子N端均含有一個NAC結構域，保守性較高；C端為激活結構域，多態性較明顯。TdNAM基因與已公布的 TtNAM-B2基因同源性最高，其編碼蛋白僅有3處氨基酸差異，B亞基中精氨酸代替蘇氨酸，其A、C和D亞基保守，可能影響植物生長。在蛋白序列中丙氨酸的替換容易影響蛋白質構像，大麥中HvNAM1、HvNAM2在102位丙氨酸到脯氨酸的差異，其籽粒蛋白含量上升[17]，本實驗獲得的TdNAM也有丙氨酸的替換現象，可能影響蛋白結構。史紅飛等[18]克隆得到的 TaNAC2基因編碼的蛋白預測定位于細胞核中，可能調節籽粒性狀。Uauy等[4]克隆得到的 TtNAM-B1、 TtNAM-B2、 TtNAM-A1位于染色體不同的位置，三者表達模式相近，均能不同程度地促進籽粒中蛋白和鋅、鐵等微量元素的積累。通過PredictProtein預測發現，TdNAM與TtNAM-B2在功能位點無差異，TdNAM應具有促進籽粒蛋白含量積累的功能。下一步將通過構建遺傳群體來驗證其在提高籽粒中蛋白含量及Zn、Fe等微量元素含量方面的功能。

本研究從858份小麥材料中篩選獲得3份含有功能性NAM基因的栽培二粒小麥，其中1個為新等位變異，為下一步選育高蛋白、微量元素含量的小麥品種提供材料和基因資源。

[1]SHEWRY P R.Improving the protein content and composition of cereal grain [J].JournalofCerealScience,2007,46(3):239-250.

[2]WHO.The World Health Report:Reducing risks and promoting healthy life [R].Geneva,Switzerland:World Health Organization,2001.

[3]JOPPA L R,DU C,HART G E,etal.Mapping gene(s) for grain protein in tetraploid wheat(TriticumturgidumL.) using a population of recombinant inbred chromosome lines [J].CropScience,1997,37:1586-1589.

[4]UAUY C,DISTELFELD A,FAHIMA T,etal.ANACgene regulating senescence improves grain protein,zinc,and iron content in wheat [J].Science,2006,314:1298-1301.

[5]HU X G,WU B H,YAN Z H,etal.Characteristics and polymorphism ofNAMgene fromAegilopssectionSitopsisspecies [J].AfricanJournalofAgriculturalResearch,2012,7(37):1586-1589.

[6]OOKA H,SATOH K,DOI K,etal.Comprehensive analysis of NAC family genes inOryzasativaandArabidopsisthaliana[J].DNAResearch,2003,10:239-247.

[7]ERNST H A,OLSEN A N,SKRIVER K,etal.Structure of the conserved domain of ANAC,a member of the NAC family of transcription factors [J].EMBORepoters,2004,5(3):297-303.

[8]OLSEN A N,ERNST H A,LEGGIO L L,etal.NAC transcription factors:structurally distinct,functionally diverse [J].TRENDSinPlantScience,2005,10(2):1360-1385.

[9]JAMAR C,LOFFET F,FRETTINGER P,etal.NAM- 1 gene polymorphism and grain protein content inHordeum[J].JournalofPlantPhysiolgy,2010,167(6):497-501.

[10]XUE G P,BOWER N I,MCINTYRE C L,etal.TaNAC69 from the NAC superfamily of transcription factors is up-regulated by abiotic stresses in wheat and recognises two consensus DNA-binding sequences [J].FunctionalPlantBiology,2005,33(1):43-57.

[11]WANG Y G,REN X F,SUN D F,etal.Origin of worldwide cultivated barley revealed byNAM- 1 gene and grain protein content [J].FrontPlantScience,2015,6:1-12.

[12]MESFIN A,FROHBERG R C,KHAN K,etal.Increased grain protein content and its association with agronomic and end-use quality in two hard red spring wheat populations derived fromTriticumtugidumL.var.dicoccoides[J].Euphytica,2000,116:237-242.

[13]KHAN K,FROHBERG R,OLSEN T,etal.Inheritance of gluten protein components of high protein hard red spring wheat lines derived fromTriticumturgidumvar.dicoccoides[J].CerealChemistry,1989,166:397-401.

[14]DISTELFELD A,UAUY C,FAHIMA T,etal.Physical map of the wheat high-grain protein content gene Gpc-B1 and development of a high-throughput molecular marker [J].NewPhytologist,2006,168:753-763.

[15]HAGENBLAD J,ASPLUND L,BSLFOURIER F,etal.Strong presence of the high grain protein content allele of NAM-B1 in Fennoscandian wheat [J].TheoreticalandAppliedGenetics,2012,125(8):1677-1686.

[16]胡 云,徐如宏.小麥高蛋白質含量基因 Gpc-B1的分子標記研究[J].天津農業科學,2014,20(9):6-10.

HU Y,XU R H.Markers linked to high protein content genes Gpc-B1 in wheats [J].TianjinAgriculturalSciences,2014,20(9):6-10.

[17]CAI S G,YU G,CHEN X H,etal.Grain protein content variation and its association analysis in barley [J].BMCPlantBiology,2013,13(35):1-11.

[18]史紅飛,高 翔,陳其皎,等.小麥NAC轉錄因子的基因克隆與序列分析[J].麥類作物學報,2011,31(3):395-401.

SHI H F,GAO X,CHEN Q J,etal.Genes cloning and sequences analysis of NAC transcripter factor in wheat[J].JournalofTriticeaeCrops,2011,31(3):395-401.

[19]ROGERS S O,BENDICH A J.Extraction of DNA from milligram amounts pf fresh,herbarium and mummified plant tissues [J].PlantMolecularBiology,1985,5:69-76.

[20]ASPLUNDA L,HAGENBLADA J,LEINO M W.Re-evaluating the history of the wheat domestication gene NAM-B1 using historical plant material [J].JournalofArchaeologicalScience,2010,37(9):2303-2307.

Exploring New Alleles ofNAMGene Related with Grain Protein Content in Wheat

ZHU Ping1,2,FENG Bo1,XU Zhibin1,FAN Xiaoli1,WANG Tao1

(1.Chengdu Institute of Biology,Chinese Academy of Science,Chengdu,Sichuan 610041,China;2.University of Chinese Academy of Science,Beijing 100049,China)

NAMgene is a member of the NAC family.It can accelerate leaf senescence,and increase the protein remobilization to grains and promote absorption of zinc,iron and other trace elements.To select the excellent quality wheat with the high protein content and high trace elements,858 wheat lines were screened by co-dominant marker Xuhw89 and the full-length ofNAMgene were amplified.Our results showed there were threeTriticumdicoccumvarieties possess the functionalNAM.The nucleotides sequences analysis showed thatNAMof F37 and F40 were respectively consistent with the TaNAC2(GenBank No.HQ872050) and TtNAM-A1(GenBank No.DQ869672).TheNAMof F34 had the closest homology with the TtNAM-B2(GenBank No.DQ869676).Compared to TtNAM-B2,five variation sites were found in F34,including one difference in intron and four in exons,resulting in three amino acid variations.The prediction analysis showed that the variation had no effect on the protein function.This new allelic variation ofNAMgene was named asTdNAM(GenBank No.KU529317).In this study,the new allele ofNAMcould provide the materials and gene resources for breeding wheat varieties with high protein and trace elements content.

Wheat;NAMgene; New allele; Molecular marker; Sequence analysis

時間：2016-07-07

2016-02-25

2016-04-12

國家轉基因專項(2016ZX08009003)

E-mail：849429995@qq.com

馮 波(E-mail:fengb2600@163.com)

S512.1；S330

A

1009-1041(2016)07-0866-06

網絡出版地址：http://www.cnki.net/kcms/detail/61.1359.S.20160707.1530.012.html