小麦PIN基因家族的鉴定及表达分析

奥鹏网院作业

小麦PIN基因家族的鉴定及表达分析
刘培勋，万洪深，郑建敏，罗江陶，蒲宗君

（四川省农业科学院作物研究所/农业农村部西南地区小麦生物学与遗传育种重点实验室，成都 610066）

摘要：【目的】PIN（puroindoline）是植物所特有的一类蛋白家族，对控制小麦籽粒硬度有重要功能。分析小麦PIN家族成员在全基因组的分布、结构及进化，研究其在不同组织的表达特异性以及在不同硬度种子的表达模式，为阐明小麦PIN基因家族的生物学功能奠定基础。【方法】根据已报道的小麦PIN基因和大麦HIN基因，利用BLASTP和HMM方法，在最新发布的小麦中国春参考序列中鉴定小麦PIN基因家族成员。利用UniProt、URGI、PFAM、CDD、expVIP等数据库，采用Clustal X、MEGA 7.0、ExPASy、MEME、GSDS、TB tools、GraphPad Prism5等软件进行生物信息学分析。采用qRT-PCR方法检测TaPIN基因家族在不同籽粒硬度小麦样品中的表达情况。【结果】共鉴定出19个小麦PIN基因，集中成簇分布于第1、5和7染色体同源群，编码148—327个氨基酸，编码蛋白相对分子量为16.39—37.19 kD，等电点为6.35—9.34。通过结构域和系统发育分析，可将19个TaPIN基因分为A和B两大类。大部分TaPINs基因仅有1个外显子，没有内含子，顺式作用元件分析发现其上游序列包含大量抗逆和种子发育相关的调控元件。转录组分析表明该基因家族在小麦籽粒中相对表达量很高，而在根茎叶等其他组织几乎不表达。实时荧光定量PCR表明，各基因间相对表达量差异显著，TaPIN9和TaPIN10表达量较高。随着小麦籽粒硬度降低，TaPIN9和TaPIN10表达上调，且表达比例增加，而TaPIN16和TaPIN6则呈现相反的趋势。【结论】小麦籽粒硬度的调节以Pina和Pinb为主，该基因家族其他成员也具有相同结构域，推测也具有相似功能，但受表达量低的限制，对籽粒硬度影响较小。从该基因的进化关系看，粗山羊草与小麦亲缘关系最近，其次是燕麦、黑麦和大麦。

关键词：小麦；籽粒硬度；PIN基因；表达分析

0 引言
【研究意义】小麦是全球重要的粮食作物，中国作为世界上最大的小麦生产国和消费国，其生产对保障国家粮食安全具有重要意义。小麦籽粒硬度是国内外小麦分类和定价的重要因素之一，是当前专用小麦育种的重要指标。不同硬度小麦面粉用途不同，硬麦适用于制作面包和面条等，软麦适用于制作饼干、糕点和酿酒等。明确小麦籽粒硬度遗传基础，对于改良小麦品质，培育具有不同硬度的小麦品种具有重大意义。【前人研究进展】前人研究发现，小麦籽粒硬度受多个遗传因子控制并受环境影响。GREENWELL等[1]首次用SDS方法从软质小麦水洗淀粉中分离获得了与籽粒硬度相关的Friabilin蛋白，其分子量为15 kD，该蛋白位于小麦5DS染色体臂的Ha位点，在软粒小麦中大量表达，在硬质小麦中含量极少或者缺失[2]。小麦籽粒硬度受PIN（puroindoline）基因控制[3]，包括Pina和Pinb，2个基因的氨基酸序列相似性达60%，共同作用影响籽粒柔软度。当二者同时存在时，小麦籽粒表现为软质，当Pina缺失或Pinb发生变异时，小麦籽粒表现为硬质。HEINZE等[4]将PIN基因转入硬粒小麦（Triticum turgidum L.var. durum），导致其籽粒品质发生改变并显著降低了磨粉所需能量。GASPARIS等[5]运用RNA介导的基因沉默技术使PIN基因沉默，发现PIN蛋白显著减少，小麦籽粒硬度增加。WILEY等[6]运用组织印记等技术，发现小麦PIN蛋白的合成和积累仅发生在籽粒胚乳细胞中。世界各地学者对不同地区小麦资源PIN基因的变异类型开展了研究并发现了大量变异类型[7-10]。随着对PIN基因研究的深入，与硬度相关的grain softness protein-1（gsp-1）也被分离出来，该基因与PIN基因高度同源，且共同存在于5DS染色体上，3个基因的物理顺序为PinB-PinA-Gsp1[11]。此外，WILKINSON等[12]在普通小麦cDNA中发现一种与puroindoline b核苷酸序列相似性高达70%的类PIN基因，称为puroindoline b-2。陈锋等[13]发现该基因有4种变异类型，分别为Pinb-v1、Pinb-v2、Pinb-v3和Pinb-v4。利用缺体材料将其分别定位为7DL、7BL、7B和7AL上。研究表明，小麦籽粒硬度不仅由基因决定，还受环境及栽培措施等的影响。王月福等[14]研究了氮肥用量对不同小麦品种籽粒蛋白质组分含量变化及加工品质的影响，结果表明，增施氮肥能够显著提高籽粒蛋白质各组分的含量，提高湿面筋含量和沉降值，从而增加籽粒硬度。李友军等[15]研究表明，氮钾配合施用可降低小麦品种豫麦50和郑麦9023的籽粒硬度。然而，其改变小麦籽粒硬度的分子遗传基础尚不明确。【本研究切入点】目前，对于小麦籽粒硬度的研究主要集中于PIN基因（Pina和Pinb），对其所属基因家族的同源基因的研究尚不全面。【拟解决的关键问题】本研究以小麦基因组参考序列（iwgsc_refseqv1.0 2018）为基础，鉴定PIN基因家族成员，分析基因分布、基因结构、保守结构域、顺式作用元件和系统进化关系等，明确小麦PIN家族基因在不同组织不同发育阶段的表达特异性和不同硬度小麦样本中的表达模式，为解析该基因功能奠定基础，为小麦遗传改良提供理论参考。

1 材料与方法
1.1 试验材料和处理方法
小麦品种蜀麦969和川麦601来自于四川省。蜀麦969是筋力最强、硬度最高的四川主推品种之一，川麦601是筋力最弱、硬度最低的品种之一。于2018年10月至2019年5月种植于四川省农业科学院德阳试验基地（31°14′N，104°24′E）。播种时各施75 kg·hm-2的磷肥（过磷酸钙）和钾肥（硫酸钾）作为基肥，常规组施150 kg·hm-2的尿素，低氮组施75 kg·hm-2的尿素，尿素的施用分为2次，一半作为基肥施用；另一半作为拔节肥施用。在开花30 d时，收取饱满的种子，液氮速冻后，-80℃保存备用，用于提取RNA。使用单粒谷物测定仪KSCS4100检测成熟小麦种子籽粒硬度指数，使用近红外谷物品质分析仪Infratec TM 1241检测粗蛋白含量。

1.2 小麦PIN基因家族的鉴定和系统发育分析
从UniProt数据库[16]下载小麦已知PIN蛋白序列[17]和大麦hordoindoline序列[18]，作为搜索小麦同源基因的探针序列（电子附表1）。小麦全基因组序列、蛋白序列和注释文件来源于URGI[19]。根据小麦PIN蛋白序列和大麦hordoindoline序列，使用HMMER3.0软件[20-21]建立HMM模型，在小麦全基因组蛋白序列中搜索同源序列。同时利用本地BLASTP程序[22]，同源比对小麦蛋白序列，进行鉴定，同时满足2种鉴定方法的被认为是小麦PIN基因家族成员。对鉴定出的小麦PIN基因，利用Clustal X 2.0进行序列比对，利用MEGA 7.0软件[23]最大似然法构建进化树。

1.3 小麦PIN基因的序列分析
使用PFAM和CDD和数据库[24-25]对鉴定出的小麦PIN基因进行保守结构域鉴定。利用ExPASy进行预测小麦PIN基因家族蛋白序列的分子量和等电点[26]。使用MEME在线程序[27]鉴定蛋白保守motif，通过GSDS在线软件[28]分析基因结构。此外，使用PlantCARE软件[29]对每个基因编码序列上游2 000 bp的DNA序列分析顺式元件，并运用TBtools工具[30]进行可视化作图。

1.4 小麦PIN基因的表达分析
从Wheat Expression Browser（powered by expVIP）[31-32]下载小麦品种中国春在不同时期不同组织的转录组表达数据，用于分析小麦PIN基因的表达模式，采用实时荧光定量PCR对小麦PIN基因的表达进行鉴定。使用Omega biotek公司植物RNA提取试剂盒（R6827）提取开花后30 d的小麦籽粒RNA，并反转录成cDNA。根据小麦PIN基因序列，利用Oligo 7软件设计定量PCR引物（表1），由上海生工生物工程有限公司合成。实时荧光定量PCR反应体系（20 μL）包括10 μL SYBR®Green Realtime PCR Master Mix、1 μL primer#1 (10 μmol·L-1)、1 μL primer #2 (10 μmol·L-1)、1 μL Template DNA和4 μL PCR grade water。反应程序为95℃ 2 min；95℃ 15 s，54℃ 30 s（同时收集信号值），40个循环，溶解曲线分析为54℃—95℃。使用β-Actin作为相对定量的内参基因，采用2-ΔΔCt方法进行数据计算，基因的相对表达量使用（平均值±标准差）表示。

表1 TaPINs实时荧光定量PCR分析所用引物
Table 1 Primer sequences used in quantitative real-time PCR analysis of TaPINs

pagenumber_ebook=9,pagenumber_book=2323
2 结果
2.1 小麦PIN基因的鉴定和序列分析
通过BLAST和Hidden Markov Model方法从小麦基因组数据库中共鉴定到19个小麦PIN家族基因，并随机命名为TaPIN1—TaPIN19，这些基因集中成簇分布于1B/1D、5A/5B/5D、7A/7B/7D染色体（电子附表2），编码148—327个氨基酸，相对分子量为16.39—37.19 kD，等电点为6.35—9.34（电子附表3）。Conserved Domains search分析结果表明，所有TaPINs基因都包含AAI_SS结构域（cl07890），该结构域是AAI_LTSS超级家族的特征，因此，PINs基因属于AAI_LTSS超级家族成员。AAI_LTSS超级家族为高等植物特有，包括谷物型阿尔法淀粉酶抑制剂（alpha amylase inhibitor）、脂质转移蛋白（lipid transfer protein）、种子储存蛋白（seed storage protein）和其他类似蛋白。Pfam数据库进一步分析表明，19个PIN蛋白结构域的Top Hit不同，其中，9个蛋白（TaPIN6—TaPIN10、TaPIN14—TaPIN16和TaPIN19）为Tryp_alpha_amyl（PF00234.22），剩下的10个（TaPIN1—TaPIN5、TaPIN11—TaPIN13、TaPIN17—TaPIN18）为Gliadin（PF13016.6）（电子附表4）。

2.2 小麦PIN基因的系统发育分析、蛋白保守结构域及基因结构分析
19个TaPINs蛋白序列的系统发育分析显示该基因家族可以分为A和B两大类，其中，作为参考的小麦PIN蛋白序列和大麦hordoindoline序列与A类聚为同一类。A类又进一步分为A1、A2和A3亚类，A1包括TaPIN9和TaPIN10，A2包括TaPIN14、TaPIN15、TaPIN16和TaPIN19，A3包括TaPIN6、TaPIN7和TaPIN8；B类可进一步分为B1和B2亚类，B1包括TaPIN11、TaPIN12、TaPIN13、TaPIN17和TaPIN18，B2包括TaPIN1、TaPIN2、TaPIN3、TaPIN4和TaPIN5（图1）。为了研究小麦PIN基因的结构特征，根据其进化关系进行保守motif分析，通过MEME数据库共鉴定12个保守motif（图2），所有序列均包含motif1、motif3、motif4和motif5，其中，属于PIN基因特征之一的富色氨酸区域位于motif3。A类中独有motif6和motif7，而B类中特有motif 8—motif12。如图2所示，除TaPIN15有1个长的内含子和TaPIN1有2个短内含子外，大部分TaPINs基因仅有1个外显子，没有内含子。19个TaPINs基因上游序列含多个顺式作用元件（图3），包括与抗逆相关的脱落酸作用元件（ABRE）、茉莉酸甲酯作用元件（CGTCA-motif和TGACG- motif）、干旱胁迫反应元件（MBS）、水杨酸作用元件（TCA-element）、抗逆元件（TC-rich repeats）、赤霉素反应元件（P-box和GARE-motif）和低温反应元件（LTR）；与种子发育相关的胚乳发育元件（GCN4_motif）、蛋白代谢调控元件（O2-site）和种子特异调控元件（RY-element）。丰度最大的3种元件为ABRE、CGTCA和TGACG，均与抗逆性有关。

2.3 小麦PIN基因在不同组织、不同时期的表达分析
利用小麦转录组数据分析了小麦品种中国春在不同组织和发育时期的表达（图4和电子附表5），19个TaPINs基因在小麦籽粒中的相对表达量显著高于其他组织部位，而在根茎叶等组织中表达量低，几乎不表达。在小麦拔节期和开花期几乎不表达，而在籽粒形成期在籽粒中大量表达表明TaPINs主要在籽粒发育过程中起作用。不同基因成员间表达差异显著，TaPIN14、TaPIN15和TaPIN12等相对表达量较低。

2.4 小麦PIN基因在不同硬度小麦中的相对表达分析
选取了3个不同籽粒硬度的小麦样品（SM969MN籽粒硬度指数为40.2，粗蛋白含量为13.4%；SM969LN硬度指数为35.3，粗蛋白含量为11.4%；CM601LN籽粒硬度指数为25.4，粗蛋白含量为10.7%），分析了TaPINs基因在不同硬度小麦中的表达模式（图6），在3个样品中，TaPIN9和TaPIN10均表达量较高，但随着施氮量降低或是在弱筋小麦品种中，其他几个基因表达量相对较低，而TaPIN9和TaPIN10表达比例非常高。比较几个表达量相对较低的基因的表达模式发现，仅有TaPIN9和TaPIN10的表达量随籽粒硬度的减低，相对表达量上调。TaPIN16和TaPIN6则呈现出相反的趋势，随籽粒硬度的降低，表达量下调。其他几个基因（TaPIN14、TaPIN15和TaPIN19）相对表达量很低。

3 讨论
根据系统发育树，可将本研究鉴定出的19个小麦PIN基因分为A和B两大类，通过同源比对和查阅小麦基因组蛋白注释可知，TaPIN9为Pina，TaPIN10为Pinb，属于A类。Pfam数据库分析表明，19个PIN蛋白结构域的Top Hit不同，A类的9个基因为Tryp_alpha_amyl（protease inhibitor/seed storage/LTPfamily），B类的10个基因为Gliadin（Cys-rich Gliadin N-terminal），推测至少A类的9个基因在影响小麦籽粒硬度中发挥相似功能，而B类则与另一类与PIN基因高度同源的基因类别Avenin-like b（AVNLB）高度相关[33]。A类又分为3个亚类，A1亚类包含了Pina和Pinb。A2亚类包括TaPIN16、TaPIN19、TaPIN14和TaPIN15，该亚类与Pinb亲缘关系较近，且与之前报道的Pinb的同源基因puroindoline b-2的变异类型数量吻合，均为4个，推测其为已报道的puroindolineb-2。前人利用小麦中国春缺体材料对该类基因进行了定位，认为Pinb-2v1位于7DL，Pinb-2v2位于7BL，Pinb-2v3位于7B，Pinb-2v4位于7AL[34]。本研究利用生物信息学方法，在小麦中国春的基因组中进行鉴定，得出TaPIN14和TaPIN15均位于7A，TaPIN16位于7B，TaPIN19位于7D染色体，与前人研究略有差异。A3亚类包括TaPIN6、TaPIN7和TaPIN8，分别位于5A、5B和5D染色体，蛋白序列比对发现，TaPIN8为已报道的gsp-1，但位于5A和5B的同源序列TaPIN6 和TaPIN7之前未见报道。A2和A3亚类均在小麦的3套染色体中存在，而A1则只存在于5D染色体，这可能是由于小麦漫长的进化过程中发生了大量基因重组（转座因子插入、基因组缺失、重复和逆序）造成[35]。

pagenumber_ebook=11,pagenumber_book=2325
图1 TaPINs蛋白序列构建的系统发育树
Fig.1 Phylogenetic tree of TaPINs using the complete protein sequences

PINA和PINB为已报道的小麦（Triticum aestivum）PIN蛋白，AINA和AINB为已报道的燕麦（Avena sativa）中同源蛋白，SINA和SINB为已报道的黑麦（Secale cereale）中同源蛋白，HINA和HINB为已报道的大麦（Hordeum vulgare）中同源蛋白，EMT21351和EMT21353为已报道的粗山羊草（Aegilops tauschii）中同源蛋白
PINA and PINB were puroindoline proteins reported in Triticum aestivum, AINA and AINB were homologous proteins in Avena sativa, SINA and SINB in Secale cereale, HINA and HINB in Hordeum vulgare, EMT21351 and EMT21353 in Aegilops tauschii, respectively

pagenumber_ebook=11,pagenumber_book=2325
图2 19个TaPINs蛋白的保守序列和基因结构图
Fig.2 The conserved motifs and gene structure analysis of 19 TaPINs

因TaPIN15太长（17 881 bp），在作图时手动删除了中间的部分内含子，最终保留了1 200 bp，以便于基因结构展示
The sectional intron of TaPIN15 gene was deleted manually, as it was too long (17881 bp), and only 1200 bp was retained for showing the gene structure

pagenumber_ebook=12,pagenumber_book=2326
图3 19个TaPINs基因顺式作用元件分析
Fig.3 Cis-acting element analysis of 19 TaPINs

ABRE：脱落酸作用元件；CGTCA-motif：茉莉酸甲酯作用元件；TGACG-motif：茉莉酸甲酯作用元件；GCN4_motif：胚乳发育元件；P-box：赤霉素反应元件；MBS：干旱胁迫反应元件；LTR：低温反应元件；O2-site：蛋白代谢调控元件；TC-rich repeats：抗逆元件；TCA-element：水杨酸作用元件；GARE-motif：赤霉素反应元件；RY-element：种子特异调控元件
ABRE: cis-acting element involved in the abscisic acid responsiveness; CGTCA-motif: cis-acting regulatory element involved in the MeJA-responsiveness;TGACG-motif: cis-acting regulatory element involved in the MeJA-responsiveness; GCN4_motif: cis-regulatory element involved in endosperm expression;P-box: Gibberellin-responsive element; MBS: MYB binding site involved in drought-inducibility; LTR: cis-acting element involved in low-temperature responsiveness; O2-site: cis-acting regulatory element involved in zein metabolism regulation; TC-rich repeats: cis-acting element involved in defense and stress responsiveness; TCA-element: cis-acting element involved in salicylic acid responsiveness; GARE-motif: Gibberellin-responsive element; RY-element:cis-acting regulatory element involved in seed-specific regulation

pagenumber_ebook=12,pagenumber_book=2326
图4 TaPINs基因在小麦品种中国春不同发育时期不同组织中的表达情况
Fig.4 Expression profiles of TaPINs in different growth stages and different tissues of Chinese Spring

A：三叶期的根；B：抽穗期的穗子；C：开花期的穗子；D：开花期的茎；E：开花后2 d的籽粒；F：五叶期的根；G：七叶期的叶子；H：开花后14 d的籽粒；I：开花后20 d的糊粉层；J：开花后10 d的胚乳；K：开花后20 d的胚乳；L：开花后30 d的胚乳
A: Root_3 leaf stage; B: Spike_emergence; C: Spike_anthesis; D:Stem_anthesis; E: Grain_2dpa; F: Shoot_5 leaf stage; G: Leaf_7 leaf stage;H: Grain_14dp; I: Aleurone layer_20dpa; J: Endosperm_10dpa; K:Endosperm_20dpa; L: Endosperm_30dpa

pagenumber_ebook=13,pagenumber_book=2327
图5 TaPINs基因家族中A类部分基因相对表达分析
Fig.5 Relative expression of TaPINs belong to class A

图中数值为3次重复平均值。SM969MN：蜀麦969，常规氮处理；SM969LN：蜀麦969，低氮处理；CM601LN：川麦601，低氮处理
The value in the figure is the average value of three replicates.SM969MN: SM969, Medium nitrogen application; SM969LN: SM969, Low nitrogen treatment;CM601LN: CM601, Low nitrogen treatment

鉴定出的19个小麦PIN基因与近缘种大麦、燕麦[36]、黑麦[37]中已报道的同源蛋白共同构建的进化树见图1中A1亚类。可见，4个麦类作物均含有该基因，其他麦类未发布完整基因组，因此，其含有PIN同源基因的数量未知。PINA与粗山羊草中的EMT21351聚在一起，PINB与EMT21353聚在一起，且蛋白序列几乎一致，再次验证了小麦D基因组来源于粗山羊草。从该基因的进化分析来看，粗山羊草与小麦亲缘关系最近，燕麦和黑麦次之，大麦与小麦亲缘最远。

通过对该基因家族编码区上游2 000 bp的DNA序列进行顺式作用元件分析，发现启动子区除了种子发育相关元件外，还存在大量抗逆相关元件，推测该类基因除影响籽粒硬度外，还参与小麦植株后期的抗性调控。DHATWALIA等[38]认为puroindoline在植物种子中具有抗菌活性，可能是一种membranotoxin膜毒素，在植物抵御微生物病原体的防御机制中发挥作用。KIM等[39]研究表明，puroindoline具有体外抗菌活性，通过转基因手段进行过表达，证明其是一种有效的抗真菌蛋白。HANEY等[40]对puroindoline抗菌活性进行了检测，表明puroindoline衍生的多肽均与DNA结合，阻断体内大分子的合成。虽然序列信息和前人研究均表明puroindoline参与抗菌活动，但基于该基因家族仅限于在种子内表达，其具体对哪些病害起作用，其调控方式等尚需进一步试验来明确。

目前，对Pina和Pinb的研究尚不能完全解释小麦籽粒硬度变异的机制。AYALA等[41]研究了西班牙南部安达鲁西亚小麦地方品种puroindoline基因的分子特征，发现13个硬质小麦样本，仅有7个可以用Pina-D1或Pinb-D1突变来解释。NIRMAL等[42]对硬籽粒小麦的基因型进行了鉴定，表明除了Pina和Pinb变异，还存在其他基因影响小麦硬度。CHEN等[43]利用电子显微镜分析的硬质和软质节节麦，不能用PIN等位基因的序列变异来解释。DARLINGTON等[44]对3个硬质大麦和3个软质大麦中的puroindoline的同源基因hordoindoline进行了序列比对，大麦中硬质与软质品种中该基因并没有差异。这些研究表明硬度是受多基因控制的数量性状。除Pina和Pinb表达量较高外，其他几个同源基因表达量非常低，对籽粒硬度的影响有限。除Pina和Pinb之外，该基因家族其他成员同样具有AAI_SS功能结构域，其在籽粒硬度调控中也可能具有共同作用，如果这些基因的表达量提升到相近水平后，能否改良籽粒硬度性状，尚需进一步试验来明确。

4 结论
在小麦全基因组范围中共鉴定出19个TaPINs基因家族成员，可分为两大类。与已报到的其他麦类作物同源基因构建的系统发育表明，从该基因来看，小麦和粗山羊草亲缘关系最近，其次是燕麦、黑麦和大麦。该基因家族均只在籽粒中表达，但部分基因表达量很低。小麦籽粒硬度的调节以Pina和Pinb为主，推测其他基因家族成员也具有相似功能，但受表达量低的限制，对籽粒硬度影响较小。

References

[1] GREENWELL P, SCHOFIELD J D.A starch granule protein associated with endosperm softness in wheat.Cereal Chemistry, 1986,63: 379.

[2] JOLLY C J, RAHMAN S, KORTT A A, HIGGINS T J V.Characterisation of the wheat Mr 15000 "grain-softness protein" and analysis of the relationship between its accumulation in the whole seed and grain softness.Theoretical and Applied Genetics, 1993,86(5): 589-597.

[3] GIROUX M J, MORRIS C F.Wheat grain hardness results from highly conserved mutations in the friabilin components puroindoline a and b. Proceedings of the National Academy of Sciences of the United States of America, 1998, 95(11): 6262-6266.

[4] HEINZE K, KISZONAS A M, MURRAY J C, MORRIS C F,LULLIEN-PELLERIN V.Puroindoline genes introduced into durum wheat reduce milling energy and change milling behavior similar to soft common wheats.Journal of Cereal Science, 2016, 71: 183-189.

[5] GASPARIS S, ORCZYK W, ZALEWSKI W, NADOLSKA-ORCZYK A.The RNA-mediated silencing of one of the Pin genes in allohexaploid wheat simultaneously decreases the expression of the other, and increases grain hardness.Journal of Experimental Botany, 2011,62(11): 4025-4036.

[6] WILEY P R, TOSI P, EVRARD A, LOVEGROVE A, JONES H D,SHEWRY P R.Promoter analysis and immunolocalisation show that puroindoline genes are exclusively expressed in starchy endosperm cells of wheat grain.Plant Molecular Biology, 2007, 64(1/2):125-136.

[7] 陈锋, 董中东, 程西永, 詹克慧, 许海霞, 崔党群.小麦puroindoline及其相关基因分子遗传基础研究进展.中国农业科学, 2010, 43(6):1108-1116.CHEN F, DONG Z D, CHENG X Y, ZHAN K H, XU H X, CUI D Q.Advances in research of molecular genetics of puroindoline and its related genes in wheat.Scientia Agricultura Sinica, 2010, 43(6):1108-1116.(in Chinese)

[8] FEIZ L, MARTIN J M, GIROUX M J.Creation and functional analysis of new Puroindoline, alleles in Triticum aestivum.Theoretical and Applied Genetics, 2009, 118(2): 247-257.

[9] WILKINSON M, WAN Y, TOSI P, LEVERINGTON M, SNAPE J,MITCHELL R A, SHEWRY P R.Identification and genetic mapping of variant forms of puroindoline b expressed in developing wheat grain.Journal of Cereal Science, 2008, 48(3): 722-728.

[10] ALI I, SARDAR Z, RASHEED A, MAHMOOD T.Molecular characterization of the puroindoline-a and b alleles in synthetic hexaploid wheats and in silico functional and structural insights into Pina-D1.Journal of Theoretical Biology, 2015, 376: 1-7.

[11] GOLLAN P, SMITH K, BHAVE M.Gsp-1 genes comprise a multigene family in wheat that exhibits a unique combination of sequence diversity yet conservation.Journal of Cereal Science, 2007,45(2): 184-198.

[12] WILKINSON M, WAN Y, TOSI P, LEVERINGTON M, SNAPE J,MITCHELL R A, SHEWRY P R.Identification and genetic mapping of variant forms of puroindoline b expressed in developing wheat grain.Journal of Cereal Science, 2008, 48(3): 722-728.

[13] CHEN F, BEECHER B S, MORRIS C F.Physical mapping and a new variant of puroindoline b-2genes in wheat.Theoretical and Applied Genetics, 2010, 120(4): 745-751.

[14] 王月福, 于振文, 李尚霞, 余松烈.施氮量对小麦籽粒蛋白质组分含量及加工品质的影响.中国农业科学, 2002, 35(9): 1071-1078.WANG Y F, YU Z W, LI S X, YU S L.Effects of nitrogen application amount on content of protein components and processing quality of wheat grain.Scientia Agricultura Sinica, 2002, 35(9): 1071-1078.(in Chinese)

[15] 李友军, 熊瑛, 骆炳山.氮、钾及其互作对两种质型小麦品质性状的影响.干旱地区农业研究, 2006, 24(2): 43-47.LI Y J, XIONG Y, LUO B S.Effects of nitrogen, potassium and their interactions on quality characteristics of two different gluten wheat cultivars.Agricultural Research in the Arid Areas, 2006, 24(2): 43-47.(in Chinese)

[16] UniProt: The universal protein knowledgebase.Nucleic Acids Research, 2016, 45(D1): D158-D169.

[17] DAY L, BHANDARI D G, GREENWELL P, LEONARD S A,SCHOFIELD J D.Characterization of wheat puroindoline proteins.FEBS Journal, 2006, 273(23): 5358-5373.

[18] BEECHER B, BOWMAN J, MARTIN J M, BETTGE A D, MORRIS C F, BLAKE T K, GIROUX M J.Hordoindolines are associated with a major endosperm-texture QTL in barley (Hordeum vulgare).Genome, 2002, 45(3): 584-591.

[19] ALAUX M, ROGERS J, LETELLIER T, FLORES R, ALFAMA F,POMMIER C, GUERCHE C.Linking the International Wheat Genome Sequencing Consortium bread wheat reference genome sequence to wheat genetic and phenomic data.Genome Biology, 2018,19(1): 111.

[20] EDDY S R, PEARSON, WILLIAM R.Accelerated profile HMM searches.PLoS Computational Biology, 2011, 7(10): e1002195.

[21] FINN R D, CLEMENTS J, EDDY S R, FINN R D, CLEMENTS J,EDDY S R.HMMER web server: Interactive sequence similarity searching.Nucleic Acids Research, 2011, 39(Web Server issue):29-37.

[22] ALTSCHUL S F.Gapped BLAST and PSI-BLAST: A new generation of protein detabase search programs.Nucleic Acids Research, 1997,25: 3389-3402.

[23] KUMAR S, STECHER G, TAMURA K.MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets.Molecular Biology & Evolution, 2016, 33(7): 1870.

[24] FINN R D, COGGILL P, EBERHARDT R Y, EDDY S R, MISTRY J,MITCHELL A L, SALAZAR G A.The Pfam protein families database: Towards a more sustainable future.Nucleic Acids Research,2015, 44(D1): D279-D285.

[25] MARCHLER-BAUER A, BO Y, HAN L, HE J, LANCZYCKI C J,LU S, GWADZ M.CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures.Nucleic Acids Research,2017, 45(D1): D200-D203.

[26] GASTEIGER E.ExPASy: The proteomics server for in-depth protein knowledge and analysis.Nucleic Acids Research, 2003, 31(13):3784-3788.

[27] BROWN P, BAXTER L, HICKMAN R, BEYNON J, MOORE J D,OTT S.MEME-LaB: Motif analysis in clusters. Bioinformatics, 2013,29(13):1696-1697.

[28] HU B, JIN J, GUO A Y, ZHANG H, LUO J, GAO G.GSDS 2.0: An upgraded gene feature visualization server.Bioinformatics, 2014,31(8): 1296.

[29] LESCOT M.Plant CARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences.Nucleic Acids Research, 2002, 30(1): 325-327.

[30] CHEN C, XIA R, CHEN H, HE Y.TBtools, a Toolkit for Biologists integrating various biological data handling tools with a user-friendly interface.BioRxiv, 2018: 289660.

[31] RAMÍREZ-GONZÁLEZ R H, BORRILL P, LANG D, HARRINGTON S A, BRINTON J, VENTURINI L, KHEDIKAR Y.The transcriptional landscape of polyploid wheat.Science, 2018, 361(6403): eaar6089.

[32] BORRILL P, RAMIREZ-GONZALEZ R, UAUY C.expVIP: A customisable RNA-seq data analysis and visualisation platform.Plant Physiology, 2016, 170(4): 2172-2186.

[33] ZHANG Y, HU X, ISLAM S, SHE M, PENG Y, YU Z, ZHANG J.New insights into the evolution of wheat avenin-like proteins in wild emmer wheat (Triticum dicoccoides).Proceedings of the National Academy of Sciences of the United States of America, 2018, 115(52):13312-13317.

[34] CHEN F, BEECHER B S, MORRIS C F.Physical mapping and a new variant of puroindoline b-2 genes in wheat.Theoretical & Applied Genetics, 2010, 120(4): 745-751.

[35] CHANTRET N, SALSE J, SABOT F, RAHMAN S, BELLEC A, LAUBIN B, GAUTIER M F.Molecular basis of evolutionary events that shaped the hardness locus in diploid and polyploid wheat species (Triticum and Aegilops).The Plant Cell, 2005, 17(4):1033-1045.

[36] GASPARIS S, ORCZYK W, NADOLSKA-ORCZYK A.Sina, and Sinb, genes in triticale do not determine grain hardness contrary to their orthologs Pina, and Pinb, in wheat.BMC Plant Biology, 2013,13(1): 190.

[37] GAZZA L, TADDEI F, CONTI S, GAZZELLONI G, MUCCILLI V,JANNI M, OVIDIO R D, ALFIERI M, REDAELLI R, POGNA N E.Biochemical and molecular characterization of Avena indolines and their role in kernel texture.Molecular Genetics and Genomics, 2015,290(1): 39-54.

[38] KIM K H, FEIZ L, DYER A T, GREY W, HOGG A C, MARTIN J M,GIROUX M J.Puroindoline: Antimicrobial wheat endosperm specific protein.Journal of Phytopathology, 2011, 7: 903-906.

[39] KIM K H, FEIZ L, DYER A T, GREY W, HOGG A C, MARTIN J M,GIROUX M J.Increased resistance to penicillium seed rot in transgenic wheat over‐expressing Puroindolines.Journal of Phytopathology,2012, 160(5): 243-247.

[40] HANEY E F, PETERSEN A P, LAU C K, JING W, STOREY D G,VOGEL H J.Mechanism of action of puroindoline derived tryptophan-rich antimicrobial peptides.Biochimica et Biophysica Acta(BBA)-Biomembranes, 2013, 1828(8): 1802-1813.

[41] AYALA M, GUZMÁN C, PEÑA R J, ALVAREZ J B.Genetic diversity and molecular characterization of puroindoline genes(Pina-D1 and Pinb-D1) in bread wheat landraces from Andalusia(Southern Spain).Journal of Cereal Science, 2016, 71: 61-65.

[42] NIRMAL R C, FURTADO A, WRIGLEY C, HENRY R J.Influence of gene expression on hardness in wheat.PLoS ONE, 2015, 11(10):e0164746.

[43] CHEN M, WILKINSON M, TOSI P, HE G, SHEWRY P.Novel puroindoline and grain softness protein alleles in Aegilops species with the C, D, S, M and U genomes.Theoretical and Applied Genetics,2005, 111(6): 1159-1166.

[44] DARLINGTON H F, ROUSTER J, HOFFMANN L, HALFORD N G,SHEWRY P R, SIMPSON D J.Identification and molecular characterisation of hordoindolines from barley grain.Plant Molecular Biology, 2001, 47(6): 785-794.

Genome-Wide Identification and Expression Analysis of PIN Genes Family in Wheat

LIU PeiXun, WAN HongShen, ZHENG JianMin, LUO JiangTao, PU ZongJun
(Crop Research Institute, Sichuan Academy of Agricultural Sciences/Key Laboratory of Wheat Biology and Genetic Improvement on Southwestern China, Ministry of Agriculture and Rural Areas, Chengdu 610066)

Abstract: 【Objective】The puroindoline (PIN) gene is a family specific to plants and plays an important role in controlling the grain hardness in wheat.In silico identification and expression analysis of PIN family genes in Triticum aestivum on whole genome level lay the foundation for elucidation the biological function of wheat PIN family.【Method】The known sequences of wheat PIN proteins and barley hordoindoline (HIN) proteins were used to query the newly released wheat peptides dataset of Chinese spring variety with HMM and BLASTP profiles.UniProt, URGI, PFAM, CDD, expVIP and other databases were used for bioinformatics analysis with Clustal X, MEGA 7.0, ExPASy, MEME, GSDS, TBtools, GraphPad Prism5 and other programs.The expression profiles of TaPIN genes in different wheat harness seeds were validated with quantitative real-time polymerase chain reaction(qRT-PCR).【Result】A total of 19 TaPINs were identified in wheat genome, which were clustered in homologous groups of chromosome 1, 5 and 7.These TaPIN proteins contained 148 to 327 amino acids, as their relative molecular weight varied from 16.39 to 37.19 kD and the isoelectric points ranged from 6.35 to 9.34.Phylogenetic and conserved domain analyses showed that the 19 TaPIN proteins were divided into A and B categories, respectively.Most TaPIN genes harbored only one exon, and there were many cis-acting regulatory elements involved in stress responsiveness and seed-specific regulation.RNA-Seq showed that this gene family expressed largely in wheat grain and hardly expressed in other tissues.The qRT-PCR results indicated that the relative expression level among TaPIN genes were significantly changed.TaPIN9 and TaPIN10 were highly expressed.The expression levels of TaPIN9 and TaPIN10, as well as its expression ratio, were up-regulated with the decrease of wheat grain hardness, while TaPIN16 and TaPIN6 showed an opposite trend.【Conclusion】Pina and Pinb genes were the main factors regulating grain hardness in wheat.It is speculated that other members of this gene family have similar functions, but have little influence on grain hardness because of low expression level.According to the evolutionary relationship of this gene, Aegilops tauschii is most closely related to wheat,followed by oats, rye and barley.

Key words: Triticum aestivum; grain hardness; PIN gene; gene expression