农学毕业论文

农学毕业论文

--WORD格式---可编辑--水稻和大麦耐盐性的比较基因组学研究郭秋林（农学，3120100118）摘要：大麦（HordeumvulgareL.）是最耐盐的禾本科作物之一，其野生种质中蕴含丰富的耐性基因资源，这些基因通过各种分子机制调控植株的耐盐性。然而，水稻是禾本科作物中耐盐差的物种之一。但大麦和水稻在基因组序列及其表达模式上具有高度相似性。这些基因往往较保守，在生长中具有功能显著性和等效性，称为共有基因，含有一些转录因子、酶家族、结构蛋白等，这些基因的作用可能与耐盐性相关。本实验选用以水稻日本晴(Oryzasativasubsp.Japonica)和西藏野生大麦XZ26（Hordeumvulgatesubsp.spontaneum）为代表材料，用转录组测序寻找盐胁迫下大麦与水稻基因组中响应盐胁迫的共有基因和物种间特异表达基因，结合对应的耐盐表型与元素含量进行分析，初步探索大麦和水稻基因组与耐盐性的关系，从基因层面揭示大麦和水稻物种间耐盐的相关机理，为水稻耐盐性改良提供理论参考。关键词：盐胁迫；差异表达基因；表现型；元素含量;转录组测序。ComparativegenomicsanalysisofsalttolerancebetweenriceandbarleyQiu-linGuo(Agronomy,3120100118)Abstract:Amonggramineouscrops,barley(HordeumvulgareL.)isoneofthemostsalt-tolerantspecies,whichcontainsrichresourcesofsalt-tolerancegenesregulatingsalttoleranceofplantsinitswildgermplasm.However,riceissaltsensitive.Barleyandricehavehighsimilarityinthegenomesequencesandtheexpressionpatternofgenes.Theseconservativegenesinbarleyandrice(cBRgenes)includingtranscriptionfactors,enzymefamilies,structuralproteinsmaintainthelawofgrowthregulation,whichplayingimportantrolesinsalttolerance.Inthisstudy,ricecultivarNipponbare(Oryzasativasubsp.Japonica)andaTibetanwildbarleyXZ26(Hordeumvulgatesubsp.Spontaneum)wereselectedforthisresearch.ThecBRsdifferentlyregulatedinriceandbarleyundersaltstresswereidentifiedthroughRNA-seq.ThenthecorrelationbetweencBRsgenesandsalttoleranceinbarleyandricewasdiscussed,revealingthemechanismofsalttolerancebetweenthetwocropsatthegenomelevel,whichcouldprovideavaluablereferenceforimprovingandbreedinginrice.Keywords:Saltstress;differentiallyexpressedgenes;phenotype;elementcontent;RNA-seq.--\n--WORD格式---可编辑--目录1.引言32.研究背景综述32.1大麦的耐盐机理32.2大麦与水稻的耐盐性比较42.3转录组及其测序技术（RNA-seq）43.材料与方法43.1实验材料53.2植株处理与取样53.3元素含量测定63.4RNA提取与RNA-Seq分析63.5数据分析84.文献综述94.1大麦和水稻耐盐性比较94.2盐胁迫对大麦及水稻元素含量的影响104.3盐胁迫下转录组分析145.讨论225.1响应基因的表达与耐盐性的关系225.2大麦和水稻耐盐性差异的基因组层面的探讨235.3研究展望及不足之处的改进23参考文献24--\n--WORD格式---可编辑--1.引言土地盐渍化是全球农业生产所面临的严峻挑战之一，盐害的产生来源多种多样，最普遍的是植物蒸腾作用引发的地下矿质资源上升到土表，这一趋势随着全球变暖和精耕细作进一步加剧。目前，全球大约有8.31亿hm2的土地已受到盐渍化的威胁，而我国盐渍土面积约为3600万hm2[3]。盐害对植物的毒害作用主要有渗透抑制、矿质营养失调和离子毒害等[2]。水稻是我国种植面积最大、单产最高粮食作物，总产目前位居我国第二。作为主要粮食，稻米在全国居民粮食消费中占65％以上。水稻的年种植面积约占全国粮食作物的30％，总产量则要占40％以上，建国以来，我国稻米生产有了较大的发展，总产量不断增加，为我国粮食生产做出了重大的贡献[3,29]。耐盐性作为水稻品种改良方向之一，备受关注[6]。水稻对于土壤盐分具有极大的敏感性，这严重威胁到了我国粮食供应，影响了农业生产，经济损失预计超过120亿[4-5]。大麦是禾本科耐盐作物的代表之一，且资源分布与适应性广[6]。因此，研究大麦的耐盐性与水稻的盐敏感性可为揭示植物耐盐机理提供重要的参考和指导；对大麦耐盐相关基因组的研究及其与水稻的同源基因比较，将对克隆耐盐基因与培育耐盐水稻新品种提供重要的理论依据和技术支持。2.研究背景综述2.1大麦的耐盐机理大麦适应性广，耐盐性强，其野生抗逆种质资源丰富大麦作为一种耐盐作物[1]。近年来，一些研究者对青藏高原一年生野生大麦开展了遗传多态性分析和干旱、酸铝及盐害等非生物胁迫耐性鉴定和优异种质的发掘，发现青藏高原一年生野生大麦具有丰富的遗传多态性和耐非生物胁迫的遗传变异[2,4]。浙江大学大麦课题组从西藏野生大麦材料中鉴定到耐盐性具有明显差异的种质资源，如XZ16、XZ26、XZ169等，其中XZ26具有优异的耐盐性[3,23]。普遍认为，渗透胁迫、离子毒害和次生胁迫是盐胁迫对植物造成危害的主要机制[1,2]。如一些植物受盐胁迫时，Cl-的毒害作用往往表现在抑制植物的光合作用[6]，而细胞中过多的Na+积累将导致离子失衡并产生特异性损伤。目前，渗透调节、离子平衡和抗氧化作用被认为是植物的主要耐盐机制[29]。如植物在胁迫下产生大量的活性氧自由基（ROS），造成一系列次生胁迫，植物通过体内的抗氧化系统清除掉一部分的ROS，转基因植物实验研究中，过量表达抗氧化酶相关基因可增强转化植株对渗透和氧化胁迫的耐性[9]。大麦耐盐机制涉及到多种生理途径与调节[2]，最主要的是离子转运平衡，包括渗透调节的改善、根部Na+的摄入与Na+回流根际土壤的调节、木质部离子运输的调控、细胞中K+、Na+的置换与保留、细胞内Na+的区隔以及对氧化胁迫的耐受等[31]。总体而言可被分为以下几方面：（1）Na+排出机制；（2）K+保留；（3）氧化胁迫耐性；（4）渗透调节[27,31]--\n--WORD格式---可编辑--。阐明各个机制的作用原理及其在大麦耐盐性上的相对作用，对于指导大麦及其它作物的耐盐性研究无疑具有重要的意义。2.2大麦与水稻的耐盐性比较相比较而言，大麦与水稻大约在5000万到1亿年发生分化，具有相当悠久的历史，但在其基因组的构成乃至RNA修饰、蛋白结构乃至功能上仍然具有较高的相似性。基因若不受功能约束，其表达形式很容易变化，长久以来使萌发调节规律得以保持的水稻、大麦基因，表达序列与形式较恒定，在萌发中具有功能显著性和等效性，有较高的研究价值[23]。此外，大麦中还存在特有的保守基因序列。鉴于大麦具有较强的耐盐性，而这从根本上又是由基因进行调控的，因而将耐盐性的差异划归为基因组的比较也是十分必要的。水稻栽培种的耐盐性普遍低于大麦，以水稻在盐处理下的生长状况作为参考指标，研究表明水稻幼苗期对盐碱较敏感，同一品种在发芽期和幼苗期的耐盐性存在一定差异，整个营养生长阶段，水稻耐盐性逐渐增强[5]。此外，相关研究[17]指出在盐碱胁迫下，水稻成熟期的株高和秆长不能很好地反映品种间耐盐性强弱，最佳取样期在分蘖期结束前，而单株分蘖数等是衡量水稻耐盐碱强弱的良好指标[4]。相关研究表明，在检测到的大麦、水稻结构与功能基因中，同源数量达到57661个，大麦特有基因为4650个，即在漫长进化史中，大麦分化出少部分具有特异性的基因[2]。而水稻与大麦在同源基因的蛋白质序列和表达形式上具有很高的相似性[23]。基于同源基因组和大麦特异基因组的研究价值，分析特有基因在表达过程中的调控机理，比较同源基因调控差异，从中找出耐盐相关基因的可能性较高。这些基因组的分析以及文库制备等也能够成为耐盐水稻开发的一项技术储备。2.3转录组及其测序技术（RNA-Seq）转录组是特定细胞或组织在特定时间或状态下转录出来的所有RNA的集合。通过对转录组的研究可以揭示生物体的基因表达、研究结构变异及发现新基因等。转录组分析的研究方法、研究平台发生着日新月异的变化，同时生物信息学分析的内容也在逐渐完善。分子生物学技术的快速发展极大地提高了人们对转录组(Transcriptome)的分析能力，特别是随着高通量测序技术的发展，RNA-Seq分析已经逐渐成为一种常用的实验手段[7]。作为一种新的转录组研究手段，RNA-Seq利用新一代测序技术能够更为快速、准确地为人们提供更多的生物体转录信息，并在生物信息学分析等相关领域得以不断应用[10]。3.　材料与方法3.1　实验材料本实验所用的水稻和大麦分别选取日本晴和西藏野生大麦XZ26作为实验材料。各取200粒左右具有活性的休眠种子。--\n--WORD格式---可编辑--3.2植株处理与取样3.2.1种子萌发水稻种子发芽：取2个500mL的锥形瓶用去离子水洗净，各加入200mL蒸馏水，供试的的日本晴种子分别放入锥形瓶中，并标记名称，放入培养箱25℃饱和湿度进行暗室培养，每隔8小时换水一次，并注意轻摇混匀。48h后，用同温去离子水冲洗一次，取2个发芽盒，标记水洗后，垫上3层用蒸馏水浸润的发芽纸，将种子均匀摆在发芽纸上，加蒸馏水至浸没。再另用滤纸加水浸湿，覆盖在发芽盒底部，发芽盒置于25℃培养箱中。恒温箱暗培养3d，出芽后补光培养。大麦种子发芽：水稻种子播种7d后进行大麦种子发芽。供试的XZ26种子用3%H2O2表面灭菌处理20min，用去离子水冲洗4次后在发芽盒中发芽，播种方法同水稻种子。发芽盒放置在25℃的生长室中暗培养3d，出芽后补光培养。3.2.2苗期盐胁迫处理大麦发芽9天后，将水稻和大麦幼苗各自移栽到6个盛有营养液的6*8孔塑料黑箱（15L）中进行水培，水稻、大麦各分为2个处理组与1个实验组。每处理（箱）各设置4´6=24个重复，每个重复取2粒发芽状况良好的种子在海绵固定下种1个孔。大麦的水培箱用气泵连续充气。所有处理在温室中培养，培养条件为光照14h，黑暗10h，昼夜温度设置为22℃/18℃，营养液pH控制在5.6左右。定期更换营养液，配方如表1，其中大麦使用的营养液编号H1到H6，水稻使用的营养液编号R1到R3及H1到H3[2]。元素种类营养液编号化学成分浓度(g/L)每10L营养液用量(mL)大量元素（大麦）H1KNO310110Ca(NO3)2·4H2O236H2MgSO4·7H2O2464H3NH4H2PO41152微量元素H4H3BO31.851MnCl2·4H2O0.99(NH4)6Mo7O24·4H2O12.36ZnSO4·7H2O1.15H5CuSO4·5H2O0.251H6Fe(Ⅲ)-EDTA42.11大量元素（水稻）R1KH2PO424.85R2MgSO4·7H2O134.85(NH4)2SO448.2R3KNO318.55Ca(NO3)2·4H2O86.4表1水稻与大麦营养液配方水培8d起进行盐胁迫处理，每个基因型的3组处理分别将营养液中NaCl的浓度设定为0mM（对照）、100mM（处理）、150mM（处理）；所有营养液每3d更换一次，加缓冲液以控制pH在6.0左右，实验组开始时NaCl浓度每次换营养液增加50mM--\n--WORD格式---可编辑--，直到达到目标浓度，使植株逐渐去适应高盐环境。3.2.3取样方法盐处理7d（播种第33d）后对材料进行取样，取样品种为XZ26、日本晴，浓度为0mM、100mM，共4组处理。每组处理各收取生长状况最好的6个重复，每个重复分别剪取地上和地下共2部分，各自作为1个样本。地下部的根用自来水冲洗3次，连续冲洗1min，洗去吸附在根表面的离子，地上部剪去除生长状况不佳的老叶片。取样时就将6个重复先后分2批收取，每批随机3个重复。第一批样品，编号地上部S1-S12，地下部R1-R12，在80℃下烘干72h称干物重，烘干后称重。烘干后的样品用于短期盐处理的金属元素含量测定。第二批样品同为3个样本，编号地上部S1-S12，地下部R1-R12，立即放于液氮冷冻烘干，-80℃保存备用于代谢含量测定（RNA-seq）[2,31]。3.3元素含量测定取样后，将地上部和地下部的第一批烘干后样品，研碎后称取重量。将样本加入试管，冷却后再加入10mlHNO3:H2O（1:1），用试管加热器以160℃消煮约2h，直至消煮液只剩试管半球部分。用滤纸过滤，收集提取液并定容至20mL，保存于试管并标记。各样本Na、K、Ca、Mg、Cu、Fe、Mn和Zn离子含量用电感耦合等离子体发射光谱仪ICP-OES（iCAP6000型号，美国赛默飞世尔科技公司）测定。3.4RNA提取与RNA-Seq分析3.4.1总RNA提取将第二批样品送样，每组取液氮冷却后的样本约0.5g，其中地上部取第三完全展开叶，采用QIAGEN的RNeasyPlantMiniKit(50)试剂盒提取总RNA，于-80℃冷库备用。3.4.2mRNA分离和纯化通过两次纯化，使用结合有poly-T寡核苷酸的磁珠从总RNA中分离纯化出含有poly-A的mRNA。步骤如下：确保总RNA起始量有1ug-5ug，加NucleasefreeWater补齐到50ul；取15ulNEBNextOligod（T）25beads到另一个新的PCR管中；取100ul2XRNABindingBuffer洗磁珠两次；往磁珠中加入50ul2XRNABindingBuffer和第一步中50ulRNA，吸打混匀；放入PCR仪，65℃5min；当PCR仪的温度达到4℃时，就拿出PCR管，室温放置5min，使mRNA充分绑定到磁珠上；将PCR管放到磁珠板上2min，使已绑定到磁珠上的mRNA分离出来；吸走上清液，加入200ulWashBuffer洗掉未绑定的RNA；置于磁珠板上2min，吸走上清液，再洗一次；往磁珠中加入50ulElutionbuffer，吸打混匀；放入到PCR仪中，80℃2min，25℃保存；当PCR仪温度达到25℃，立即拿出PCR管，往PCR管中加入50ulRNABindingBuffer，吸打混匀；室温放置5min，使mRNA充分绑定到磁珠上；按以上步骤加200ulWashBuffer洗两次；往磁珠中加入10ulNucleasefreeWater混匀，置于PCR管中80℃2min；立即将PCR管放到磁珠板上分离，等液体变澄清后，吸出6.75ul到新的PCR管中，即为--\n--WORD格式---可编辑--mRNA。3.4.3mRNA热打断加入纯化后的的6.75ulmRNA、2ulFirstStrandSynthesisReactionBuffer、0.5ulRandomPrimers，混匀后放到PCR仪中，94℃11min，时间一到立即拿出置于冰上。3.4.4mRNA反转录往上一步结束的反应中加入0.25ulMurineRNaseInhibitor、0.5ulProtoScriptIIReverseTranscriptase，混匀后按以下程序进行反应合成第一条cDNA链：25℃10min42℃15min70℃15min4℃hold往合成的第一条连中加入24ulNucleasefreewater、4ulSecondStrandSynthesisReactionBuffer、2ulSecondStrandSynthesisEnzymeMix放到PCR仪中，16℃1h，热盖温度40℃3.4.5cDNA纯化（1）往上述结束的dscDNA中加入72ul(1.8x)AMPureXPbeads（2）吸打混匀，室温放置5min（3）放置于磁珠板上，等液体变澄清（大约5min），小心吸除上清（4）往PCR管中加入200ul80%的酒精，洗30s，吸除上清，洗两次（6）往已晾干的磁珠中加入30ul的NucleasefreeWater，吸打混匀磁珠（7）将PCR管放到磁珠板上分离（8）等液体变澄清后，吸出27.75ul到新的PCR管中3.4.6末端修复通过切除或者补平3’和5’突出端，得到平末端。对5’端磷酸化后，再对3'端加‘A’，使其能与下一步含有‘T’的Adapter互补连接加入1.5ulNEBNextEndPrepEnzymeMix、3.25ulNEBNextEndRepairReactionBuffer(10x)、27.75ul纯化后dscDNA，混匀后将PCR管放入PCR仪中，按以下程序运行：20℃30min65℃30min保存于4℃3.4.7接头连接在双链cDNA片段上加上Adapter，使其能杂交到flowcell上往末端修复结束的PCR管中加入7.5ulBlunt/TALigaseMasterMix、1.25ulDilutedNEBNextAdaptor、1.25ulNucleasefreewater，混匀后将PCR管放于PCR仪中，设置为20℃15min，时间到立即加入1.5ulUSERenzyme，混匀后放置于PCR仪中，设置为37℃15min。--\n--WORD格式---可编辑--3.4.8磁珠筛选用0.6xAMPureXPbeads抓取600bp以上的DNA片段，再用0.8x的浓度对600bp以下的DNA片段进行300bp以上片段的抓取，筛选出300-600bp的范围富集PCR，步骤如下：（1）往接头连接结束的PCR管中加入1.25ulNucleasefreeWater，补齐到45ul（2）加入27ul(0.6x)AMPureXPbeads，吸打混匀，室温放置5min（3）放置于磁珠板上，等液体变澄清（大约5min），吸取上清到新的PCR管中（4）再次加入9ulAMPureXPbeads，吸打混匀，室温放置5min（5）放置于磁珠板上，等液体变澄清（大约5min），吸除上清（6）往PCR管中加入200ul80%的酒精，洗30s，吸除上清，洗两次（7）将PCR管放置于磁珠板上，打开盖子晾干，使其完全龟裂（大约5min）（8）往已晾干的磁珠中加入13ul的NucleasefreeWater，吸打混匀磁珠（9）将PCR管放到磁珠板上分离（10）等液体变澄清后，吸出11.5ul到新的PCR管中3.4.9PCR富集PCR选择性的富集两端含有接头的DNA片段，放大cDNA文库。向文库中引入检索引物以用于多样本测序。加入11.5ul筛选后的DNA片段、12.5ulNEBNextHighFidelity2xPCRMasterMix、0.5ulIndex（X）Primer、0.5ulUniversalPCRPrimer，混匀放入PCR仪中，按以下程序进行：预变性98°30s1个循环变性98°10s退火65°30s15个循环延伸72°30s终延伸72°5min1个循环保存4°∞3.4.10胶回收纯化用2%的琼脂糖凝胶进行电泳，切取350-450bp之间的片段，用QIAGEN胶回收试剂盒进行回收，剩余2ul用Qubit测RNA浓度。3.4.11文库接头效率测定使用AglientTechnologies2100Bioanalyzer仪器对每个文库进行条带大小的确定，然后再用KAPA的文库定量试剂盒进行文库接头效率的测定，以确保数据的准确性。3.5数据分析采用Excel（MicrosoftOfficeExcel2013）进行相对干物重的计算，先用t检验去除偏差值较大的样本组数，剩余数据求平均值。相对干物重=实验组平均干物重对照组平均干物重--\n--WORD格式---可编辑--对2个品种的实验组、对照组的各个元素含量的数据进行处理。按元素分类，计算均值、标准差，绘制图表。采用Blast2go3.3软件对RNA-seq检测到的各品种地上部、根部及根部的mRNA序列文库进行检索及注释。计算表达上调或下调倍数（Foldchange）。FoldChange=实验组平均表达量对照组平均表达量=X1+X2+X3Y1+Y2+Y3X1~X3是实验组三个重复的表达量，日本晴：S1~S3，R1~R3，XZ26：S7~S9，R7~R9Y1~Y3是实验组三个重复的表达量，日本晴：S4~S6，R4~R6，XZ26：S10~S12，R10~R12筛选出FDR＜0.001，相对表达FoldChange＞2或＜0.5的差异表达基因（mRNA序列），将同一基因表达出的转录组进行归纳，合并为单个基因，列出各品种各部位差异基因表。从差异基因表中选出水稻和大麦中具有高度序列相似性的共同表达基因，再从中选出在一对多或多对一中评估成绩较高的基因，列出共同差异基因表，以待进一步比较分析。4.结果与分析4.1大麦和水稻耐盐性比较图1三种盐浓度处理下大麦与水稻表现型的差异如图1，与对照组相比，在100和150mM浓度的NaCl处理条件下，大麦组、水稻组的生长情况均受到一定影响。图1为在盐处理7d取样时，各组的生长形态。对照组水稻生长状况良好，叶片正常呈现绿色、挺直；对照组大麦的生长状况良好，叶片宽大茂密，但基部老叶倒伏较严重，部分叶片有发黄，这与预期有所不符。在盐浓度100mM下，大麦总体--\n--WORD格式---可编辑--长势与对照组相比变化不大，叶片疏密适宜，呈绿色，叶宽减小；而水稻组则开始呈现一定影响，主要表现在新叶生长速率减缓，老叶开始发黄，根系生长减缓，部分叶片叶尖下垂。150mMNaCl处理下，大麦生长受到一定影响，但并不明显，仅生长速率稍减缓；与之相对，水稻所有品种的生长均受到严重影响，从株高来看，生长速率缓慢，部分植株基本停止生长，接近死亡，叶片窄小、蜷缩萎蔫。以XZ26和日本晴分别作为大麦、水稻的代表品种，第一次取样所测得干物重结果如表2。可见去除水分后与对照相比，XZ26和日本晴在100mM下干物重呈基本不变或减小趋势，由相对干物重可见变化幅度不大，大约在0.75至1.0之间，日本晴根部及XZ26地上部干物重对盐处理较为敏感。就整珠而言，高浓度盐含量时，XZ26的相对干物重随浓度增加减少更明显。表2盐处理前后大麦与水稻根和地上部干物重比较编号“S”为地上部，“R”为根部，干物重为差异显著性分析后求求平均所得部位品种处理浓度干物重/g相对干物重地上部日本晴对照0.19650.958100mM0.1883XZ26对照0.15480.815100mM0.1261根部日本晴对照0.04610.757100mM0.0349XZ26对照0.04240.892100mM0.03784.2盐胁迫对大麦及水稻元素含量的影响保持离子平衡是大麦的重要耐盐机制，也是作物正常生长的基本要素之一。在大致了解大麦与水稻主要品种的耐盐性差异后，即可测定它们对照组、100mM两个盐浓度水平的地下部和地上部的大量元素K、Ca、Na及微量元素Mg、Mn、Fe、Cu、Zn含量。从图2和图3中可以直观了解盐胁迫对大麦及水稻两个组织中元素含量的影响，其结果对该时期基因组功能及表达量变化等进一步研究具有参考价值。以日本晴和XZ26分别作为水稻和大麦的代表品种，选取盐处理后7d进行测定。盐处理后，日本晴、XZ26地上部、地下部的Na含量均显著上升，升至对照组的10-30倍，其中根部XZ26处理前为0.745mg•g-1，与日本晴的0.725mg•g-1接近，处理后上升至22.24mg•g-1较之日本晴含量更高。而地上部情况则有所不同，处理前日本晴和XZ26分别为0.278mg•g-1与0.284mg•g-1非常接近，盐处理后XZ26上升至33.52mg•g-1远低于日本晴的52.18mg•g-1，两个品种的变化趋势与根部相反，说明大麦在高盐环境下，Na自根部到地上部的转运趋向于受到抑制，体现了大麦品种的离子平衡机制。日本晴、XZ26的地上部、根部在高盐环境下，K含量均显著下降，与Na含量变化呈现显著的负相关。其中日本晴根部由39.0mg•g-1降至19.24mg•g-1，为对照组的49%，降幅最显著，XZ26由68.10mg•g-1降至52.1mg•g-1，降幅仅23%；地上部日本晴K元素含量仅降至75%，由对照组的50.83mg•g-1--\n--WORD格式---可编辑--降至38.07mg•g-1，XZ26与之相比，实验组62.76对比对照组100.78，降至62%，幅度更加明显，说明高盐环境对K的吸收产生了一定影响，并且日本晴在高盐环境下更趋向于将K元素输送到地上部，而XZ26趋向于保留在根部。两个品种的地上部、地下部的Ca含量在盐处理时发生不同程度的下降。其中对处理最不敏感的是日本晴的根部，实验组1.276mg•g-1与对照组1.152mg•g-1相比降幅10%，而XZ26根部由2.731mg•g-1降至1.802mg•g-1，降幅34%；与之相对的，地上部日本晴的实验组与对照组Ca含量分别为3.77mg•g-1与5.36mg•g-1，XZ26分别为3.84mg•g-1与6.21mg•g-1，盐处理后的Ca积累量二者相近，处理前XZ26有稍高的积累量，降幅分别为30%与38%，XZ26降幅仅仅稍大于日本晴，可见XZ26钙含量对盐处理更加敏感，而日本晴仅地上部积累量受到较大影响。微量元素方面，最值得关注的是Fe元素和Cu元素，与其他元素不同，它们的根部积累量显著高于地上部。未进行盐处理时，日本晴和XZ26根部的Fe积累量分别为0.235mg•g-1和0.822mg•g-1，大麦品种显著更高，而地上部含量分别为0.0872mg•g-1和0.0875mg•g-1，基本相同；盐处理时，日本晴根部和地上部Fe含量均一定程度的上升，分别为0.396mg•g-1和0.112mg•g-1，对照组的1.66和1.28倍，而XZ26根部和地上部Fe元素均下降，分别为0.642mg•g-1和0.0776mg•g-1，对照组的78%和89%，变化趋势相反。Cu浓度在地上部盐处理时均受到抑制，日本晴实验组0.0266mg•g-1与对照组0.0385mg•g-1相比仅为69%，XZ26实验组含量0.0097mg•g-1为对照组0.0122mg•g-1的80%；盐处理根部的Cu含量同样为抑制作用，且降幅显著，实验组0.0316mg•g-1仅为对照组0.0980mg•g-1的大约1/3，与之相对，日本晴的Cu含量小幅度上升，实验组0.465mg•g-1是对照组0.375的1.24倍。Mg含量在盐处理下变化趋势与Ca有所类似，日本晴根部基本未受到影响，100mM下Mg含量1.37mg•g-1与对照组1.36mg•g-1相比可视为不变，而XZ26也仅受到微小影响，实验组1.15mg•g-1与对照1.23mg•g-1相比达到93%；日本晴地上部分别为3.86mg•g-1与5.08mg•g-1，降幅24%，结合文献，表明Mg在由地下输送到地上部时受到一定影响，XZ26由2.44mg•g-1下降到实验组的2.03mg•g-1，对应降幅17%大于根部，但不及日本晴变化大。此外，日本晴、XZ26地上部和根部在盐胁迫下，Mn、Zn含量均发生不同幅度的下降，均为显著变化。在盐胁迫条件下，元素Na如此大幅度的上升，据猜测主要是由盐处理所引起，品种特异性的影响次之。资料表明[29]水稻地上部叶片组织代谢对于Na含量比大麦更为敏感，而本次地上部RNA-seq主要材料就是叶片，印证了取材和分析结果的正确性。除Na外，元素含量显著上升的只有日本晴地上部与根部的Fe元素，日本晴地上部的Cu元素，其他均为下降或无明显变化。可见高盐环境对于水稻和大麦的大量元素吸收具都有明显的抑制作用，不利于生长发育；微量元素在盐胁迫下可能伴随更复杂的调控方式。--\n--WORD格式---可编辑--根部Mg含量地上部Mg含量图2盐处理前后大麦与水稻元素含量差异浓度单位mg•g-1--\n--WORD格式---可编辑--地上部Fe含量根部Fe含量根部Cu含量根部Mn含量根部Zn含量地上部Cu含量地上部Mn含量地上部Zn含量根部Fe含量地上部Fe含量续图2盐处理前后大麦与水稻元素含量差异浓度单位mg•g-1--\n--WORD格式---可编辑--4.3盐胁迫下转录组分析4.3.1差异表达基因调控初步分析：通过之前所分析的生长外观表型、金属元素含量变化等，已经能大致了解水稻、大麦代表品种存在一定的耐盐性差异。而RNA-seq结果表明水稻和大麦大多数基因在盐胁迫下进行了差异表达，且差异基因的调控状况也有所不同，主要表现在不同品种、不同部位的基因上调或下调数目不同。如图3所示[30]，满足FDR＜0.001，FoldChange＜0.5或FoldChange＞2或含有表达量=0（Foldchange未知）的条件下，共筛选出地上部差异表达基因3371个，包含XZ26（大麦）187个，日本晴（水稻）3246个，地上部共同表达基因62个。有13和16个基因分别在所有品种中均上调和下调。共有基因中分别有3个和30个基因在XZ26中受到上调和下调，而在日本晴中调节方向相反。分别有36个和89个基因在日本晴中未被检测到，而在XZ26中分别上调和下调。有39个和119个基因分别仅在XZ26上调和下调，由共有差异基因和大麦特有基因的特殊性，初步推测其中含有耐盐基因。地上部根部图3大麦与水稻盐处理下差异表达基因数的比较根部差异表达基因5326个，包含XZ26（大麦）606个，日本晴（水稻）4919个，根部共同表达基因199个。有43和60个基因分别在所有品种中均上调和下调。共有基因中分别有81个和15个基因在XZ26中上调和下调，而在日本晴中调整方向相反。分别有260个和147个基因在日本晴中未被检测到，而在XZ26中分别上调和下调。有341个和162个基因分别仅在XZ26上调和下调。同理猜测共有及大麦特有基因中有应对盐胁迫的成员。差异基因越多表明对盐胁迫越敏感，根部差异表达基因显著多于地上部，水稻差异表达基因显著多于大麦。基因下调是生理活动受到抑制的表现，大麦根部基因大多上调，地上部基因大多下调，水稻则相反。而非共有基因数目显著多于共有基因说明耐盐调控方式不同。4.3.2大麦与水稻中共有的盐胁迫响应基因经过初步筛选得到地上部（表3）和地下部共有差异基因，结合Blast2go基因注释，将其归入不同的基因组，并根据其注释的功能与耐盐性之间的联系。由于软件原因，根部基因并未全部获得注释，无法进行统一，本实验仅就地上部基因进行了分析。蛋白质功能类似的基因被归为1个基因类。按作用对象筛分出A、B、C、D、E共5个基因类别。--\n--WORD格式---可编辑--A类基因中共7个成员，蛋白含有结构域，通过直接与DNA、RNA结合的蛋白而直接起到调控作用，可参与多项生理反应。1、2、48号基因的蛋白属于WRKY转录因子家族，是含高度保守的核心氨基酸序列的质膜内在蛋白，在受到病原菌、损伤、盐害等胁迫因子诱导后表达，其锌指结构域能与DNA或RNA结合，调控并参与多项生理过程包括抗病、生长发育、修复损伤等方面[28]。该基因组的这3个基因在XZ26中均下调，在日本晴均上调。16号的蛋白也在XZ26下调而在日本晴地上部显著上调，其属于VQ蛋白家族，含有VQ-motif保守序列。此外，35号锌指CCCH结构域蛋白是转录遏制物，抑制植物色素形成，通过调控GA和ABA的表达控制生长，XZ26和日本晴都上调。37号的蛋白属于螺旋-环-螺旋（HLH）DNA结构域家族，在两个品种中均发生下调。B类基因共含有6个成员，是参与激素合成或调控的酶。编号9对应蛋白为9-顺式-环氧类胡萝卜素加双氧酶，注释表明它是植物合成ABA的关键物质，在XZ26和日本晴中均上调，即促进ABA合成。编号27对应合成ACC氧化酶，主要功能是乙烯生物合成途径的限速酶，调控乙烯的生成速率，在日本晴和XZ26中表达均下调[15]。44号对应合成小麦萌发素，该类激素与逆境刺激相关，而文献[11]表明它是一类重要的胁迫响应蛋白，位于细胞外基质，使草酸分解生成H2O2，属于Cupin超家族成员，该基因在大麦中下调，水稻未知。49号参与调控FERONIA受体激酶，大麦中下调，水稻未知，该酶能正向调节生长素，抑制ABA信号转导及磷酸酶的失活，使植物生长趋向活跃；50号为脱落酸受体PYL4，与49号具有功能相似性，抑制ABA及磷酸酶失活，需在生长受抑制等胁迫环境下表达。57号为细胞分裂素-O-糖基转移酶，可以使CTK发生O-葡萄糖基化而失活，由于可通过β-葡糖苷恢复，发生这类糖基化对CTK的抑制作用影响不大，可防止CTK被降解，便于其运输、储存，该基因大麦下调，水稻显著上调[20]。C类基因共5个成员，与叶绿体的脂类及色素类调控有关的酶。30、33号基因对应α-酮戊二酸加双氧酶，与黄酮素、花青素等的生物合成相关，均在XZ26下调，日本晴上调。40号基因在XZ26上调，日本晴显著上调，参与脂氧合酶形成；46号PAP14同样存在于叶绿体并参与脂质形成，它们均在日本晴地上部显著上调[25]。62号对应A1-Igammal磷脂酶，参与叶绿体脂质代谢，可水解甘油磷脂、糖脂，在XZ26上调，日本晴下调。D类为与金属元素调控有关的酶，共有2个成员。31号在XZ26下调，日本晴上调，参与合成钙调蛋白PBP1，可结合Ca，促进其吸收，并在胁迫反应中起调节作用。47号为储铁蛋白，在叶绿体中，收集可溶性Fe2+并以Fe3+形式储存，影响光合作用，可能对Fe元素平衡有影响，在两个品种中均上调。E类为过氧化物酶，可催化过氧化氢、酚类等的分解。包括26号（XZ26显著下调、日本晴显著上调）、41号（XZ26下调、日本晴上调）、58号（XZ26下调、日本晴显著上调）。剩余编号的有标注蛋白或作用于各类蛋白质，或作用于无明确分类的小分子，由于种类繁多，对它们的功能进行进一步分类。F类为与免疫调控相关的基因组，6个成员。3、20号为病程相关蛋白，与发病原有关--\n--WORD格式---可编辑--的一种细胞外蛋白，具备多种免疫效应，XZ26下调[19]。5号在两个品种均下调，对应CBL互作激酶，与CBL蛋白结合，调节NAF域，使之以Ca依赖方式激活，激活后CBL作用于免疫反应[22]。7号对应内切几丁质酶，可以降解入侵真菌的细胞壁，XZ26下调，日本晴显著上调。17号对应GPI锚定的赖氨酸结构域蛋白，是一种细胞表面受体，日本晴下调。51号对应蛋白可引发作物过敏反应和抗病性，两个品种均下调。G类是位于细胞膜上的转运或信号蛋白，3个成员。23号的COBRA-7，离散分布，与纤维素沉积及细胞的伸长有关，在日本晴、XZ26中表达均受抑制。2号ATP与细胞膜的结合蛋白，促进细胞膜细胞骨架上ATP水解酶的催化效率，两个品种都上调。53号为丝/苏氨酸激酶受体，跨膜信号蛋白，接受胞外信号，使下游信号中丝/苏氨酸磷酸化。H类是与逆境胁迫直接相关的调控基因，9个成员。10号对应海藻糖-磷酸盐合酶，其表达在XZ26和日本晴分别上调至2.01和5.00倍。24号对应油体钙蛋白，XZ26和日本晴中分别下调和显著上调。43号对应HVA22，在XZ26和日本晴中分别上调至2.773和9.109倍。55号对应线粒体精氨酸转运蛋白，需Cl-诱导下表达，XZ26未知，日本晴显著上调至16.80倍。28、36号对应半胱氨酸受体激酶，应对胁迫、调控生长过程的保守信号组件，与病原体防御及细胞程序性死亡有关，两个基因在XZ26与日本晴中均下调[24]。29、60号对应胰蛋白酶抑制剂，能结合胰蛋白抑制底物分解，会在高盐等逆境下诱导表达，增加糖分积累，均在XZ26下调，水稻上调。此外还有之前提到的44号。除以上分类，一些对应各类参与氧化、水解、羧化等生化酶的基因被单独列出。6号α-半乳糖苷水解酶，两个品种均上调。12号枯草杆菌蛋白水解酶，XZ26下调。13号含黄素单氧酶，与含N、S、P、Se等亲和杂原子化合物的化学异物氧化代谢有关，XZ26下调，日本晴显著上调。14号对应多胺氧化酶降解多胺，均下调。15号磺基转移酶，结合磷酸盐与外源化合物或代谢产物，调控解毒、生长发育，保持黄酮化合物多样性，XZ26上调，日本晴下调[21]。18号HAD家族水解酶，XZ26下调，日本晴上调。39号对应磷酸烯醇式丙酮酸羧化酶，XZ26下调，日本晴显著上调，TCA循环中依赖Ca催化形成草酰乙酸。45号基因表达的阿魏酰辅酶A为转移酶，将阿魏酰基转移至羟基棕榈酸上形成酸酐，缺水环境下的组分存贮，两个品种均下调；而同对应酰基转移酶的56号两个品种均上调。59号对应的1-脱氧木酮糖-5-磷酸合成酶，催化3-磷酸甘油醛产生1-脱氧木酮糖-5-磷酸和丙酮酸的酮醇缩合反应生成DXP，是一种质体类异戊二烯生物合成必需的限制酶，和叶绿体的发育密切相关，该基因在XZ26中下调，而在日本晴中显著上调至70.4倍。61对应β-1-3内切葡聚糖水解酶，水解可溶性纤维素，两个品种均下调。此外还有一些功能不够明晰的基因。8号属于cupin超家族成员，在日本晴、XZ26中分别下调、上调，同属该家族的还有9、30、44号，调控方式各不相同。19号调控某种合成抑制剂，54号调控某种糖激酶，它们都在XZ26中下调，日本晴中上调。而4、11、21、25号等基因无注释，功能有待探明。--\n--WORD格式---可编辑--表3大麦（左）与水稻（右）地上部盐处理共同差异基因表达量变化及其调控蛋白No.基因MSUfoldchangeblast2go注释基因MSUfoldchangeblast2go注释1MLOC_102640.252probableWRKYtranscriptionfactor50LOC_Os05g460208.402probableWRKYtranscriptionfactor502MLOC_120790.255WRKYDNA-bindingdomainsuperfamilyLOC_Os01g402609.056WRKYtranscriptionfactor28-like3MLOC_125810.278pathogenesis-related1-17LOC_Os05g51680?pathogenesis-related1A-like4MLOC_140090.411predictedproteinLOC_Os01g486200.071---NA---5MLOC_158340.435CBL-interactingkinase21LOC_Os07g442900.293CBL-interactingkinase216MLOC_15872.478alpha-galactosidaseLOC_Os10g351107.068alpha-galactosidase7MLOC_163200.126endochitinaseLOC_Os06g5106028.679endochitinase8MLOC_1695018.027---NA---LOC_Os08g134400.179germintype19MLOC_183002.799---NA---LOC_Os03g443807.8509-cis-epoxycarotenoiddioxygenase10MLOC_187252.011probablealpha,alpha-trehalose-phosphatesynthase[UDP-forming]9LOC_Os02g548204.996probablealpha,alpha-trehalose-phosphatesynthase[UDP-forming]911MLOC_207470.492Rpr4901.2LOC_Os11g119800.371Os11g0227000,partial12MLOC_365290.414Subtilisin-likeproteaseLOC_Os04g03100?subtilisin-likeprotease13MLOC_368490.306probableflavin-containingmonooxygenase1LOC_Os04g1469065.590probableflavin-containingmonooxygenase114MLOC_373280.334probablepolyamineoxidase4LOC_Os04g575600.452probablepolyamineoxidase415MLOC_392353.716flavonolsulfotransferase-likeLOC_Os11g309100.249cytosolicsulfotransferase8-likeisoformX316MLOC_402540.407VQmotiffamilyLOC_Os06g3397030.225VQmotiffamily17MLOC_43178?lysMdomain-containingGPI-anchored1-likeLOC_Os09g278900.207lysMdomain-containingGPI-anchored1-like18MLOC_45180.440catalytichydrolaseLOC_Os03g494408.960HAD-superfamilyhydrolase,subfamilyIA,variant3,19MLOC_484290.132synthesisinhibitorIILOC_Os01g073005.361synthesisinhibitorII20MLOC_485310.277Pathogenesis-related1LOC_Os01g284505.933pathogenesis-relatedmaizeseed21MLOC_509710.303hypotheticalproteinTRIUR3_04936LOC_Os02g534105.751PREDICTED:uncharacterizedproteinLOC433088722MLOC_51830.468Pleiotropicdrugresistance4LOC_Os01g424105.957ABCtransporterGfamilymember3723MLOC_528640.391COBRA7LOC_Os07g490800.265COBRA724MLOC_529300.147caleosin2LOC_Os04g4320029.345caleosin225MLOC_536213.096predictedproteinLOC_Os04g4714011.782PREDICTED:uncharacterizedproteinLOC433664026MLOC_541290.078peroxidaseN-likeLOC_Os10g0207010.602peroxidaseN-like27MLOC_542720.361ACCoxidaseLOC_Os09g277500.2551-aminocyclopropane-1-carboxylateoxidase128MLOC_552070.175cysteine-richreceptorkinase10LOC_Os10g047200.225cysteine-richreceptorkinase1029MLOC_56460.350Bowman-BirktypetrypsininhibitorLOC_Os03g608408.734Bowman-Birktypetrypsininhibitor-like30MLOC_595960.421Naringenin,2-oxoglutarate3-dioxygenaseLOC_Os03g030345.444naringenin,2-oxoglutarate3-dioxygenase-like31MLOC_603120.258calcium-bindingPBP1-likeLOC_Os06g469505.132calcium-bindingPBP1-like32MLOC_605170.327---NA---LOC_Os11g02250?---NA---33MLOC_614970.277naringenin,2-oxoglutarate3-dioxygenase-likeLOC_Os04g491940.363naringenin,2-oxoglutarate3-dioxygenase-like34MLOC_632752.323glutamatedecarboxylaseLOC_Os03g133007.463glutamatedecarboxylase135MLOC_635252.921zincfingerCCCHdomain-containing2-likeLOC_Os05g106709.270nucleicacidbinding续表3大麦（左）与水稻（右）地上部盐处理共同差异基因表达量变化及调控蛋白--\n--WORD格式---可编辑--No.基因MSUfoldchangeblast2go注释基因MSUfoldchangeblast2go注释36MLOC_641640.451Lectin-domaincontainingreceptorkinaseLOC_Os09g169500.031---NA---37MLOC_64470.374HLHDNA-bindingdomainsuperfamilyLOC_Os02g451700.346HLHDNA-bindingdomainsuperfamily38MLOC_644750.290---NA---LOC_Os01g0673028.516receptor1239MLOC_647240.381phosphoenolpyruvatecarboxylasekinaseLOC_Os02g4158023.715phosphoenolpyruvatecarboxylasekinase40MLOC_649726.511？LOC_Os12g3726023.255lipoxygenase,chloroplastic41MLOC_652260.276rootperoxidaseLOC_Os07g480205.696peroxidase2-like42MLOC_652882.211predictedproteinLOC_Os09g0857052.887ribosomal-like43MLOC_670972.773HVA22eLOC_Os11g305009.109HVA22e44MLOC_674380.322Germin3-7LOC_Os03g59010?germin3-845MLOC_677520.399omega-hydroxypalmitateO-feruloyltransferase-likeLOC_Os11g424800.261agmatinecoumaroyltransferase46MLOC_688850.446probableplastid-lipid-associated14,chloroplasticLOC_Os01g4672028.431probableplastid-lipid-associated14,chloroplastic47MLOC_692953.315Ferritin-1,chloroplasticLOC_Os12g015304.314ferritin-1,chloroplastic48MLOC_701900.329WRKYtranscriptionfactor40LOC_Os06g4401010.132probableWRKYtranscriptionfactor6049MLOC_711250.400ReceptorkinaseFERONIALOC_Os05g25370?receptorkinaseFERONIA50MLOC_713490.337abscisicacidreceptorPYL4LOC_Os03g186000.198abscisicacidreceptorPYL4-like51MLOC_732030.332Harpin-induced1containing,expressedLOC_Os12g062200.359Harpin-induced1containing,expressed52MLOC_732082.685EHdomain-containing1LOC_Os04g5735019.957EHdomain-containing1-like53MLOC_737720.243serinethreonine-kinasereceptorLOC_Os04g013104.719---NA---54MLOC_747130.460uncharacterizedsugarkinaseslr0537LOC_Os01g016203.372uncharacterizedsugarkinaseslr053755MLOC_74725?mitochondrialargininetransporterBAC2LOC_Os01g1252016.800mitochondrialargininetransporterBAC256MLOC_751546.906acyltransferaseAt1g54570,chloroplasticLOC_Os09g335309.789acyltransferaseAt1g54570,chloroplastic57MLOC_764800.198Cytokinin-O-glucosyltransferase2LOC_Os02g5193028.456Cytokinin-O-glucosyltransferase258MLOC_791760.216Peroxidase2LOC_Os07g4801013.478peroxidase2-like59MLOC_800270.174probable1-deoxy-D-xylulose-5-phosphatesynthase2,chloroplasticLOC_Os07g0919070.248probable1-deoxy-D-xylulose-5-phosphatesynthase2,chloroplastic60MLOC_804760.393Bowman-Birktypewound-inducedproteaseinhibitorLOC_Os01g04040?wound-inducedserineproteaseinhibitor61MLOC_808490.235beta-1,3-endoglucanase,partialLOC_Os01g713500.326beta-1,3-glucanaseprecursor62MLOC_808784.456phospholipaseA1-Igamma1,chloroplastic-likeLOC_Os05g323800.234phospholipaseA1-Igamma1,chloroplastic-like4.3.3大麦中特有的盐胁迫响应基因除共有差异表达基因外，大量差异基因仅在水稻或大麦中存在。水稻和大麦分化以来不同环境的长期选择下，大麦中必然保留一些特有的基因及其表达调控，它们同样为提高大麦的耐盐性做出了贡献。所有检测出的大麦特有的响应盐胁迫的基因共125个，如表4所示。它们从离子平衡、氧化、可溶物渗透等不同途径进行调控，根据功能不同，可以筛选出一些推测的关键耐盐基因组并分类。--\n--WORD格式---可编辑--表4大麦特有的盐胁迫响应基因表达量变化及其调控蛋白基因MSUfoldchangeblast2go注释基因MSUfoldchangeblast2go注释MLOC_580640.4413-N-debenzoyl-2-deoxytaxolN-benzoyltransferase-likeMLOC_641920.260glucosyltransferaseMLOC_733800.2337-deoxyloganetinglucosyltransferase-likeMLOC_770793.383Glucosyltransferase,,expressedMLOC_526980.404ABCtransporterGfamilymember28-likeMLOC_233676.734Glutaredoxin-C1MLOC_107480.314ABCtransporterGfamilymember37MLOC_587894.996Glutaredoxin-C1MLOC_565710.102ABCtransporterGfamilymember6-likeMLOC_580410.431glycerophosphodiesterphosphodiesteraseGDPD6MLOC_60739?abscisicacidreceptorPYL4-likeMLOC_776900.337GTP-bindingSAR1AMLOC_590190.387agmatinecoumaroyltransferaseMLOC_56390.455heatshock70kDa17MLOC_145022.291aminoacidtransportMLOC_112860.356heatstresstranscriptionfactorA-4d-likeMLOC_649670.437Asparticasenepenthesin-2MLOC_74263?hypotheticalproteinF775_03829MLOC_592270.467BluecopperMLOC_787450.057hypotheticalproteinF775_03829MLOC_209013.152BTBPOZdomain-containing,partialMLOC_29100.040hypotheticalproteinF775_20921MLOC_639262.062BTBPOZdomain-containingAt3g05675-likeMLOC_617680.428hypotheticalproteinF775_32425MLOC_220130.466bZIPtranscriptionfactor60MLOC_3958511.752invertaseinhibitorMLOC_646852.813C-4methylsteroloxidaseMLOC_147870.362IQdomain-containingIQM1MLOC_442470.337caffeoylshikimateesterase-likeisoformX1MLOC_513930.258isoflavonereductaseMLOC_520190.408CBL-interactingkinase14MLOC_600220.281L-ascorbateoxidase-likeMLOC_729290.386cellnumberregulator13-likeMLOC_810030.496light-inducibleCPRF2-likeisoformX2MLOC_118750.452ceramideglucosyltransferaseMLOC_642020.349longchainacyl-synthetase4-likeMLOC_547622.119chloridechannelCLC-c-likeMLOC_666300.436L-typelectin-domaincontainingreceptorkinase-likeMLOC_77020.381cysteine-richreceptorkinase19MLOC_565070.353MLO1MLOC_232210.365cysteine-richreceptorkinase25MLOC_115480.469NADP-dependentmalicenzymeMLOC_719580.356cytochromeP450711A1MLOC_533630.211NRT1PTRFAMILY-likeMLOC_673190.434DSBA-likethioredoxindomaincontainingMLOC_593820.400NUCLEARFUSIONDEFECTIVE4-likeMLOC_650560.259E3ubiquitin-ligaseEL5-likeMLOC_61322.369OTUdomain-containingAt3g57810MLOC_662980.329endoglucanase7MLOC_631250.412pathogenesis-related1-14MLOC_136180.228ent-kaur-16-enesynthase,chloroplastic-likeMLOC_729650.265pathogenesis-related1-16MLOC_668435.464ethylene-responsivetranscriptionfactor12-likeMLOC_564710.387pentatricopeptiderepeat-containingAt1g26900,mitochondrialMLOC_432070.307Ethylene-responsivetranscriptionfactor1BMLOC_363383.327pentatricopeptiderepeat-containingAt3g29290MLOC_374922.242F-boxLRR-repeat13MLOC_662542.489pentatricopeptiderepeat-containingAt5g39980,chloroplasticMLOC_724894.935F-boxLRR-repeat3-likeMLOC_559803.303pentatricopeptiderepeat-containingAt5g46580,chloroplasticMLOC_445488.662F-boxLRR-repeatAt4g14096MLOC_70399?PeroxidaseNMLOC_572840.498ferredoxin,rootR-B1-likeMLOC_345710.198Phenylalanineammonia-lyaseMLOC_446300.369floralorganregulator1,partialMLOC_511282.766phosphoenolpyruvatecarboxykinase[ATP]-likeMLOC_45890.475GDSLesteraselipaseAt5g45910-likeMLOC_382910.134postacrosomalsheathWWdomain-binding-likeMLOC_257630.393Glucan1,3-beta-glucosidaseMLOC_119600.237PR17cprecursor--\n--WORD格式---可编辑--续表4大麦特有的盐胁迫响应基因表达量变化及其调控蛋白基因MSUfoldchangeblast2go注释基因MSUfoldchangeblast2go注释MLOC_672130.418prolineoxidaseMLOC_793350.306probablebeta-1,3-galactosyltransferase19MLOC_384000.348ReceptorkinaseHSL1MLOC_124522.669RNApolymerasesigmafactorsigCMLOC_615792.987probablecalcium-bindingCML22MLOC_748012.029RNA-bindingpno1-likeMLOC_648060.426probableglucuronosyltransferaseOs01g0926700MLOC_556630.207rootperoxidaseMLOC_719483.531probablelinoleate9S-lipoxygenase5isoformX1MLOC_652250.199rootperoxidaseMLOC_712754.435probablelinoleate9S-lipoxygenase5isoformX2MLOC_383672.105rootpimordiumdefective1MLOC_545012.577probablepolygalacturonaseAt1g80170MLOC_752230.470RPM1-interacting4-likeMLOC_678510.392probableWRKYtranscriptionfactor33MLOC_367352.305superoxidedismutase[Fe]2,chloroplasticMLOC_811310.346probableWRKYtranscriptionfactor50MLOC_615850.295sigmafactorbinding1,chloroplastic-likeMLOC_358280.353PREDICTED:uncharacterizedproteinLOC101782143MLOC_721570.436somaticembryogenesisreceptorkinase2-likeisoformX1MLOC_508392.721uncharacterizedLOC104581721isoformX1MLOC_368860.407two-componentresponseregulatorARR3-likeMLOC_509682.427uncharacterizedtransporterYBR287W-likeMLOC_750260.420UDP-glucuronicaciddecarboxylase4-likeMLOC_368252.533uncharacterizedtRNArRNAmethyltransferaseMLOC_295500.349SRG1MLOC_177790.240,partialMLOC_690912.52shockSRC2MLOC_321790.369,partialMLOC_169980.431transposon,,Mutatorsub-classMLOC_127090.329---NA---MLOC_19562.454sulfurdeficiency-induced1-likeMLOC_184510.284---NA---MLOC_216640.219subtilisin-likeproteaseMLOC_214640.377---NA---MLOC_702193.825UPF0553-likeMLOC_220002.737---NA---MLOC_718950.268vacuolar-processingenzyme-likeMLOC_24310.361---NA---MLOC_126930.316Wall-associatedreceptorkinase3MLOC_48030.379---NA---MLOC_705150.404wound-induced1-likeMLOC_56162.145---NA---MLOC_650330.488zincfinger1-likeMLOC_599140.443---NA---MLOC_159112.142predictedproteinMLOC_656930.385---NA---MLOC_405310.389predictedproteinMLOC_693250.483---NA---MLOC_739166.379predictedproteinMLOC_702580.423---NA---MLOC_77700.377predictedproteinMLOC_71480.405---NA---MLOC_816360.387predictedproteinMLOC_781310.296---NA---基因类别a与植物抗性及免疫相关。MLOC_10748（简称10748，后同）、52698、56571为ABC转运体G家族的成员，是免疫相关的防御蛋白，均下调。59019参与羟基肉桂酸酰胺的合成，这对病原体的防御作用，下调。7702、23221则与28、36号同样对应属于富半胱氨酸受体激酶，与应对胁迫和病原体防御有关。4589，GSDL脂解酶，抵抗腐生真菌黑斑病，下调。56507对应MLO1，参与病原体防御和叶细胞凋亡调控，依赖钙/钙调蛋白，对白粉病菌广有谱抗性，下调。63125、72965与3、20号同为病程相关蛋白，具有免疫效应，均下调。75223合成RPM1互作蛋白，是植物免疫体系重要调节因子，在病原菌感染下对植物抵抗侵染起到核心的作用，下调。--\n--WORD格式---可编辑--基因类别b直接参与调控DNA或RNA。20901、63926的蛋白含有BTBPOZ结构域，均上调，研究表明它们通过染色质构象的局部控制调控基因表达[33]。60739对应ABA的PYL4受体，可抑制ABA信号作用，下调。22013对应bZIP转录因子60，识别核心序列含受ABA诱导的基因的启动子元件，在植物免疫和非生物胁迫反应中起着作用，下调。67851、81131和1、2、48号基因一样属于WRKY转录因子家族。11286，热激转录因子，应对热胁迫响应，下调。43207、66843对应响应乙烯的转录因子，前者下调、后者上调。12452和74801分别对应RNA聚合酶和RNA连接酶，转录相关，均上调。基因类别c涉及蛋白质或脂质代谢。65056，泛素连接酶，将泛素分子连接到目的蛋白，使之代谢降解，下调；6132，包含OUT域，可以从蛋白质上将泛素化蛋白水解，在蛋白质的防降解及其循环中起重要调节作用，上调。37492、72489、44548均对应F-boxLRR-repeat，是SCF型泛素连接酶的底物识别组件，均上调。71958，与13号同为细胞色素P450酶，可参与外源性物质代谢，下调。13618，ent-kaur-16-ene合酶，有利于赤霉素的前体合成，下调。64202，长链酰基辅酶A合成酶，脂肪酸跨膜代谢相关，下调。66630，与36号基因同样为含植物血凝素的受体，下调。71895，液泡加工酶，可参与液泡贮藏蛋白的水解，可能与盐胁迫下的渗透调节有关，下调。基因类别d涉及离子或基团调控。54762，氯离子通道蛋白CLC，受电位调控，可能缓解盐胁迫下Cl-的毒害作用，上调。57284，铁氧还蛋白，在光合作用光反应中还原NADP，下调。61579，Ca结合蛋白，上调。77690，GTP结合蛋白，下调。14787，钙调素结合蛋白(IQM1)，与非生物胁迫相关，下调。1956号受硫元素缺乏所诱导，进行存贮硫酸盐的利用，上调，可能与离子平衡相关。36735，SOD[Fe]2，破坏盐胁迫下在细胞内产生的有毒害作用的超氧阴离子自由基，上调。基因类别e进行全局调控。72929，细胞数量调节因子，可抑制细胞数量增长，下调。56471、36338、66254、55980，三角状五肽重复区蛋白(PPR)，研究表明和植株发育有关，促进叶绿素合成，第一个在线粒体中，下调，后三个在叶绿体，上调。23367、58789，谷氧还蛋白，调节细胞发育、凋亡，均上调，可能与氧化胁迫的耐盐机制有关。44630，花器官调节因子，与花芽分化有关，下调。11286，热休克蛋白启动子，下调；5639，热休克蛋白，涉及盐胁迫下的细胞应激反应，在抗逆、抗氧化性等方面起到重要作用，下调[32]。此外有还一些彼此相关的响应基因。58064，苯甲酰基转移酶，下调；79335，半乳糖基转移酶，下调；64806，葡糖醛酸转移酶，下调；73380、11875、64192、77079，葡糖基转移酶，前三个下调，第四个上调。39585，转化酶抑制剂，显著上调。52019，与5号协同，下调。51128，与39号协同，上调。66298，与61号的协同，下调。21664与12号协同，下调。70515与60号协同，下调。12693与53号协同，下调，与细胞形态发生及发育有关。25763，β-葡糖苷酶，与57号相拮抗而使CTK恢复活性，下调。51393，异黄酮还原酶，与15号拮抗调控黄酮类化合物，下调。55263、65225与E类基因同样合成过氧化物酶，也都是下调。71275、71948对应亚油酸脂氧合酶同工酶，均上调。--\n--WORD格式---可编辑--5.讨论5.1响应基因的表达与耐盐性的关系重叠图显示地上部共有表达基因中，无论是XZ26还是日本晴，盐胁迫下根部的差异表达基因要显著多于地上部，这可能由于根际直接与高盐环境接触，需要更强的应激性，导致差异表达更加频繁。此外，大麦的下调个数要显著多于水稻，且调整幅度较小，可能想到的关联之一，是在盐胁迫下，大麦地上部基因表达活动受ABA等影响受到抑制；也可能是大麦对高盐的容忍性较高，而地上部受到Na等刺激相对较小，并未出现较强的应激反应。而水稻中，无论是与耐盐相关性较高的逆境胁迫基因，或是其他基因，其表达普遍受到了更显著的上调或下调，这表明在基因层面，水稻受到的影响更大，并出现了一定的表达紊乱现象。高盐环境下，过高的Na+浓度将会影响部分酶的活性，细胞内外的电位和水势等，从而引发胁迫环境。C类基因的脂类色素调控受到影响后，可能导致叶片颜色发生变化，且它们可能响应Cl调控光合作用[17]。F类和a类基因体现了盐胁迫与病原入侵表达模式的类似，但盐胁迫下这些基因在水稻与大麦中都普遍下调，可能为盐胁迫致使植株抗病性下降。盐胁迫下，23号基因下调可能导致细胞伸长受到抑制，植株矮小；盐胁迫影响了ATP酶的催化，相关生理合成和代谢受到抑制，52号基因上调可能是为弥补这一变化，恢复细胞膜上和骨架上的ATP代谢。43号HVA22受一些非生物逆境胁迫诱导表达，参与ABA信号转导途径，在禾本科有广泛同源性[14,25]。此外包括9号在内，多个基因的调控涉及到ABA，而ABA是一种抑制生长的植物激素，较高的ABA/IAA能使叶片发黄、脱落，延缓新陈代谢，可应对高盐逆境带来的影响。ABA还涉及基因间拮抗调节，如9号和49号。在部共同表达基因中，10号、24号、43号和55号是地上部应对盐胁迫的关键基因。相关实验表明，海藻糖-磷酸盐合酶在提高植株耐盐性中具有显著作用，它在高盐等环境下受到诱导，受胁迫环境调节，并诱导耐盐相关基因的表达有助于使植株在一段时间后缓解胁迫带来的危害[12]。24号油体钙蛋白II是水稻6个油体钙蛋白成员之一，包含3个保守的结构域，位于胞浆，植物受到盐胁迫时，一端结合Ca2+，一端磷酸化，参与逆境胁迫脂类信号途径，这些脂类中包含与胁迫信号途径重要激素茉莉酸相关的[13]。线粒体精氨酸转运蛋白受Cl-诱导，有助于脯氨酸积累以应对高渗胁迫。这些基因在高盐条件下，日本晴均表现为上调，XZ26大都上调，且日本晴的上调幅度都更大，分析认为是日本晴地上部的盐敏感性更强。有发现表明一些VQ蛋白能与WRKY转录因子(如1、2、48)相互作用，从而参与植物对抗高盐逆境的反应。49号FERONIA受体激酶经实验证明可在ABA介导进行细胞外环境碱化与防止细胞扩张的信号通路控制[26]。资料表明[11]44号germin相关基因随着盐处理时间的延长在胚芽鞘中、叶片的表达量逐渐减低，以此提高耐盐性。在XZ26地上部特有基因中，经分析判定以下4个基因组为盐胁迫关键基因上调。基因36735，分解盐胁迫下氧化反应所产生的自由基，氧化胁迫机制。基因54762，氯离子调控，缓解氯毒害，保证叶绿素功能正常等。基因组56471、36338、66254、55980，保证叶绿素--\n--WORD格式---可编辑--形成，促进光合作用。12452、74801调控DNA\RNA合成，与转录的进行乃至b类基因作用于DNA\RNA的基因调控有关。第二、三个可能与盐胁迫下XZ26叶片仍维持绿色，保证营养物质持续供应有关。5.2大麦和水稻耐盐性差异的基因组层面的探讨分析表明XZ26特有部分的基因与共有部分基因调控具有相关性，如43207、66843乙烯响应转录因子的调控可能受到27号ACC氧化酶的下调、乙烯含量增加的影响，还有资料表明它们受胁迫因子和胁迫信号转导途径的影响，并茉莉酸防御调控的关键信号集成因子。又如71275、71948亚油酸脂氧合酶同工酶的上调可能与40号脂氧合酶的上调相关。文献表明通过调控H2O2可以有效促进植物各项生理活动正常进行，如细胞壁的修饰和重建等，浓度过高将导致细胞程序性死亡[13,16,26]。此外，H2O2可以调整细胞内外pH，而消耗ATP进行质子转运形成细胞内外一定的电位梯度可排出多余的Na+[29]，H2O2的浓度将影响质子的电位差。14号的多胺氧化酶，44号的小麦萌发素在分别降解多胺和草酸时都会生成氧化产物H2O2；而E类或70399、55263、65225对应的过氧化物酶恰好能催化H2O2的分解，消除其毒性。在高盐条件下，14、44号在XZ26中下调，而共有过氧化氢酶基因均在日本晴上调、XZ26中下调，X26特有基因也均下调，体现应对胁迫时基因调节的差异和互补性[16]。47号储铁蛋白在XZ26和日本晴中地上部表达量分别上调至3.31和4.31倍，而相应的Fe元素含量分别为0.89和1.28倍，XZ26出现小幅下降。这可能与根部特性有关，也可能与元素总含量并无直接联系；此外，相关文献[18]表明盐处理会改变水稻根际氧化还原状态，吸附的Fe浓度会增加，促进铁离子被动吸收。根据31号钙调蛋白在高盐下的调控，若不考虑其他因素，这应该分别抑制和促进了XZ26和日本晴对Ca的吸收。而实际上XZ26和日本晴地上部Ca含量均下降，日本晴下降幅度更小，可见两者相对变化与猜测大致相符，但绝对值减小至少表明盐胁迫还影响了其他Ca/钙调蛋白含量变化因素[31]。5.3研究展望及不足之处的改进本结果大致符合预期，通过以XZ26和日本晴为代表的大麦和水稻外观表型、元素含量到基因组进行分析，确定了盐胁迫下的共同响应基因及部分大麦特异响应基因，为今后到耐盐基因的精细定位以及克隆，利用重点基因及基因的同源效应改善水稻耐盐性、提高产量奠定了一定的实验基础。也为禾本科作物耐盐生长机理研究和改良提供了一个有潜力的基因库资源。然而，为获得更全面、准确的实验结论，本研究仍然存在一些问题待改进，包括方案拟定、实验操作、数据处理等方面。例如：（1）由于样品设置重复数目有限，所有样品均需用于取样。而取样样品在离开营养液后需立即封装处理，取样前一刻缺少时间摆拍、观察，尤其是根部表现型的观察。因此，表现型分析缺少详细的单株对比及其图示。可以在实验前设置更多对照组。（2）用于测量干物重的第一批样品重复较少，只有3个重复，数据t检验、计算平均值的--\n--WORD格式---可编辑--结果缺乏可靠性，专用于干物重测量的重复应设置5组左右。（3）对照组大麦生长状况不佳，倒伏较严重、叶片发黄，量温室的温度、水分及营养液配置方案实施有待进一步完善，例如改善大麦根部营养液通气状况。（4）本实验没有分析根部基因表达数据，主要原因有两个，一是软件运行遇到问题，根部有注释丢失，二是若加入根部分析，会导致数据量繁多，而实验分析时间有限，篇幅也不宜过于冗长。后续有机会可进行补充实验。（5）本实验只涉及金属元素含量的测定，而若增加对于各类脂质、蛋白、激素等有机质的测定，将能更全面的阐述基因组调控差异所带来的影响极其机理，然时间有限，对于各类猜想，以后有待进一步证明。参考文献[1]李瑞芬，王雪青，王宏芝，等.野大麦耐盐适应性反应机制的研究.中国农业科学.2006,39(12):2459-2466.[2]吴德志.西藏野生大麦耐盐种质的发掘及其耐盐机制研究.浙江:浙江大学,2012.[3]李建国，濮励杰，张润森，等.土壤盐渍化研究现状及未来研究热点.地理学报.2012,67(9).[4]陈　峰，高　洁，杨敏华，等.水稻穗型的研究进展.江苏农业学报.2009,25(5):1167-1172.[5]梁正伟，杨福，王志春，等.盐碱胁迫对水稻主要生育性状的影响.生态环境.2004(1).[6]祁栋灵，郭桂珍，李明哲，等.水稻耐盐碱性生理和遗传研究进展.植物遗传资源学报.2007,8(4):486-493[7]许硕.野生大豆盐胁迫相关microRNA的功能分析.北京:中国农业科学院.2011.[8]RamC.，高安娜.谷氨酸脱羧酶活性与小麦的田间性状国外农学:麦类作物.1989,4:33-35.[9]黄有总.大麦耐盐机理及部分耐盐相关性状的QTL定位.浙江:浙江大学.2006.[10]文路，汤富酬.单细胞转录组高通量测序分析新进展.遗传(Hereditas).2014,36(11):1069-1076.[11]李莹，柳参奎.植物类萌发素蛋白研究进展.中国农学通报.2014,30(30):246-254.[12]郭蓓，胡磊，何欣，等.海藻糖-6-磷酸合成酶转基因烟草提高耐盐性的研究.ChineseBulletinofBotany.2008,25(1):41-49.[13]魏征，马红，葛晓春，等.水稻油体钙蛋白家族的进化及对干旱胁迫的响应性分析.科学通报.2011,56(20):1612-1621.[14]陈佳佳，牟巍，周树峰.玉米HVA22基因启动子的克隆与功能验证.核农学报.2014,28(4):560-568.[15]陈银华，张俊红，欧阳波，等.花椰菜ACC氧化酶基因的克隆及其RNAi对内源基因表达的抑制作用.ActaGeneticaSinica.2005,32(7):764-769.[16]刘宝宝，张士璀.多胺氧化酶在多胺代谢中的作用.鲁东大学学报(自然科学版).2014,30(4):336-341.[17]王桂君，许振文，蒋秋花，等.不同程度盐碱化土壤对两种作物种子萌发及幼苗生长的影响.安徽农业科学.2013,41(21):8857-8859，8881.[18]陈雪，刘丹青，王淑，等.不同土壤的还原状况对铁镉形态转化和水稻吸收的影响.土壤学报.2003,50(3):548-555.[19]代建丽.植物病程相关蛋白研究进展.科技探索.2008,2:123-125.[20]马新梅.拟南芥细胞分裂素O-糖基转移酶基因功能分析.山东大学出版社.2010,2:17-18.[21]董龙，罗聪件，何新华，等.芒果类黄酮磺基转移酶基因克隆及表达分析.热带作物学报.2014,35(4):706-711.--\n--WORD格式---可编辑--[1]陈奋.Cbl-b蛋白的功能及其与结构的关系.杭州师范学院学报(医学版).2007,27(1):49-52.[2]LinL.,TianSL.,ShawnK.,etc..ConservedTranscriptionalRegulatoryProgramsUnderlyingRiceandBarleyGermination.PLOSONE.2014,9(2).[3]MichaelW.,MikaelB.,JarkkoS..TranscriptionalregulationoftheCRK/DUF26groupofReceptor-likeproteinkinasesbyozoneandplanthormonesinArabidopsis.BMCPlantBiology.2010,10:10-95.[4]LiZX.,LiZ.,GaX..ROPl1GTPaseNegativelyRegulatesABASignalingbyProtectingABllPhosphataseActivityfromInhibitionbytheABAReceptorRCARl／PYL9inArabidopsis.JournalofIntegrativePlantBiology.2012,54(3):180-188.[5]YuF.,QianLC.,CandidaN.,etc..FERONIAreceptorkinasepathwaysuppressesabscisicacidsignalinginArabidopsisbyactivatingABI2phosphatase.PNAS.2012,109(36):14693-14698.[6]AyalewL.,MakiK..Insightsintothesalttolerancemechanisminbarley(Hordeumvulgare)fromcomparisonsofcultivarsthatdifferinsaltsensitivity.JournalofPlantResearch.2009,123:105–118.[7]DongHS.,TerrencePD.,DavidWB.etc..HarpininducesdiseaseresistanceinArabidopsisthroughtheeyetemicacquiredresistancepathwaymeditatedbysalicylicacidandtheNIM1gene.ThePlantJournal.1999,20(2):207-215.[8]TakujiN.,MitsuruO.,MichikoA.etc..Differencesinmechanismsofsalttolerancebetweenriceandbarleyplants.SoilScienceandPlantNutrition.1996,42(2):303-314.[9]WangXL.,WuDZ.,YangQ.,etc..IdentificationofMildFreezingShockResponsePathwaysinBarleyBasedonTranscriptomeProfiling.OriginalResearch.2016(7):106-117.[10]ZhangH.,HanB.,WangT.,etc..MechanismsofPlantSaltResponse:InsightsfromProteomics.JournalofProteomeResearch.2012,11:49-67.[11]FlowersT.,GalalH.,BromhamL..Evolutionofhalophytes:multipleoriginsofsalttoleranceinlandplants.Funct.PlantBiol.2010,37:604–612.[12]MartinoC.,DelfineS.,PizzutoR.,etc..Freeaminoacidsandglycinebetaineinleafosmoregulationofspinachrespondingtoincreasingsaltstress.NewPhytol.2003,158:455–463.--