多年生黑麦草温度胁迫差异表达转录因子的比较转录组分析

多年生黑麦草是禾本科冷季型草种,温度胁迫严重影响其分布和产量。转录因子调控基因表达在植物响应逆境胁迫中发挥重要作用。该研究以多年生黑麦草品种雅晴为材料,采用高通量RNA-seq技术,对热(40℃)、冷冻(-10℃)和对照(22℃)处理的样本进行转录因子应答分析,以探索多年生黑麦草转录因子对温度胁迫的响应规律,并筛选抗逆转录因子候选基因。结果显示:(1)共鉴定了694个转录因子unigenes,分属于AP2/ERF、GTF、HSF、MYB、NAC、WRKY、bHLH和bZIP等32个家族。(2)在热和冷冻胁迫下,ERF(AP2/ERF)、MYB、NAC和bZIP家族基因均以上调为主,而WRKY家族则多数下调。HSF、GTF和DREB(AP2/ERF)家族成员大多数受热诱导上调,而多数bHLH家族成员受冷冻胁迫上调。(3)功能富集结果表明,多年生黑麦草差异表达的转录因子基因主要参与植物激素信号转导、昼夜节律和病原体互作等通路。研究表明,多年生黑麦草中与胁迫适应相关的转录因子基因普遍上调表达,而与生长和抗病相关的基因多数下调。该研究筛选到大量抗逆转录因子候选基因,为分子育种提供抗逆基因资源。     

植物AP2/ERF转录因子家族的研究进展

作为植物最大转录因子家族之一,AP2/ERF转录因子广泛存在于植物中,并因其在基因育种等方面具有重要作用而倍受关注。AP2/ERF转录因子家族成员至少含有一段60个左右氨基酸构成的AP2保守结构域,根据结构域的数量和识别序列的不同可以将其分为AP2、ERF、DREB、RAV和Soloist 5个亚家族,且在拟南芥、水稻、玉米和番茄等植物中,AP2/ERF转录因子及其各个亚家族成员的数量各不相同。研究发现,AP2/ERF转录因子可通过响应乙烯、细胞分裂素和生长素的调节从而直接或间接参与种子发育过程、花和果实等器官的形态建成等植物发育的多个进程;除了初生代谢,AP2/ERF转录因子还在植物次生代谢尤其是在调控药用植物主要药用活性成分(如青蒿素、紫杉醇和木质素)合成方面效果显著。同时,有报道称拟南芥AP2/ERF基因具有正向调节抗灰霉病的功能,一些AP2/ERF基因也被报道在植物应对高盐、干旱、缺氧、低温等非生物胁迫方面具重要功能;另外,AP2/ERF类转录因子还参与了乙烯等介导的非生物信号传导。介绍AP2/ERF转录因子的结构分类特征、在不同植物中的数量分布,并阐述其在植物发育、次生代谢、生物与非生物胁迫和信号传导等方面的研究进展,以期为培育出兼具高产、抗病的实用型转基因作物提供理论依据。     

大叶杨干旱和水淹胁迫的转录组分析

大叶杨（Populus lasiocarpa）是中国特有的杨属物种,干旱和水淹是影响大叶杨生长和分布范围的两个关键因子。AP2/ERF转录因子家族在植物响应非生物胁迫中发挥重要作用。本研究采用转录组测序、生物信息学分析手段并结合分子实验验证初步鉴定了参与大叶杨干旱和水淹胁迫响应的关键基因。研究结果显示：（1）在大叶杨中分别鉴定到3,986/385个响应干旱/水淹胁迫的差异表达基因,其中包括237个同时响应干旱和水淹胁迫的差异表达基因。（2）在大叶杨中共鉴定到205个AP2/ERF家族成员,系统发育分析表明其在大叶杨中主要分为5个亚家族,并显著富集于差异表达基因中。（3）筛选部分胁迫前后差异表达的PlAP2/ERF基因进行qRT-PCR实验,经证实这些基因在大叶杨受到干旱/水淹胁迫时均可被诱导表达。综上,大叶杨在水淹胁迫下的差异表达基因数量明显少于干旱胁迫,AP2/ERF基因家族的部分基因参与到大叶杨干旱/水淹胁迫的应激表达过程。     

苹果DREB转录因子家族全基因组鉴定与分析

DREB转录因子属于AP2/ERF转录因子家族,能够与DRE/CRT顺式作用元件特异性结合,调控与逆境应答基因的表达,因而在植物应对低温、干旱、高盐等逆境胁迫中发挥重要作用。该研究利用苹果全基因组数据,通过生物信息学手段鉴定苹果DREB转录因子家族成员,并分析DREB转录因子家族保守域特点与功能及表达情况。结果表明:从苹果全基因组中共鉴定出60个DREB转录因子家族成员,与拟南芥和水稻相比基本一致,通过引入拟南芥DREB基因进行系统发生分析,进一步可以将其细分为6个亚组;结构域和保守元件分析表明,DREB基因家族含有一个AP2保守结构域;染色体定位表明,苹果DREB基因分布于11条染色体上,部分基因存在串联复制现象;基因结构分析显示,该亚家族基因不含内含子。利用同源拟南芥RNA-Seq数据分析结果表明,DREB转录因子家族对低温、ABA调节等非生物胁迫具有调控作用,同时在DREB亚家族中每个亚组响应不同的非生物胁迫;通过分析DREB基因在不同组织中的表达情况,结果显示DREB基因在植物根部中的表达量最强,其次是叶。     

杨树ERF11转录因子基因应答渗透胁迫表达分析

从小黑杨(Populus simonii×P.nigra)叶片中克隆出732 bp的ERF11转录因子基因(Potri.011G057000.1)cDNA,其编码的蛋白含有243个氨基酸,属于不稳定的亲水蛋白,不存在信号肽,不具跨膜能力,含AP2/ERF家族保守结构域。系统进化树表明,小黑杨ERF11蛋白与胡杨、麻疯树、蓖麻、橡胶树、木薯遗传距离较近,说明其亲缘关系较近。亚细胞定位结果表明,ERF11蛋白定位于细胞核中。ERF11基因表达具有组织特异性,并受胁迫诱导表达。其在根中相对表达水平明显比在茎和叶中表达水平高,在盐和甘露醇胁迫下,基因均表现上调表达。表明ERF11转录因子基因可能与植物应答高盐和干旱胁迫相关。     

ERF转录因子家族在甜菜块根发育中的功能分析

本研究以AP2/ERF保守结构域序列为检索序列在甜菜(Beta vulgaris)基因组数据中搜索并鉴定Bv ERF转录因子家族成员。通过系统分析,发现甜菜基因组中共有68个Bv ERF家族基因,系统进化分析将它们分为10个组。Bv ERF家族中15个成员与甜菜块根转入快速生长阶段有关,其中Bv_ammr可能与甜菜丰产性能有关,Bv_ignp、Bv_sqfr和Bv_khde与苗龄113 d时丰产型品种生长速率较高有关,Bv_wnjc可能与块根生长速率负相关,Bv_fefj基因与高糖型品系具有更多的维管束环有关。研究结果为进一步明确Bv ERF转录因子调控甜菜块根发育的分子机制提供参考。     

油橄榄AP2/ERF转录因子鉴定及水胁迫表达分析

为探究油橄榄AP2/ERF基因家族对水胁迫的响应机制,该研究对受干旱及水淹胁迫的‘佛奥''和‘TYZ-1号''2个品种的根和叶进行转录组测序,并对油橄榄中AP2/ERF转录因子的蛋白理化性质、基因结构及系统进化进行分析,同时分析与水胁迫相关的AP2/ERF转录因子在2个油橄榄品种中的基因表达差异且进行RT-qPCR验证。结果表明:(1)在油橄榄中鉴定得到110个AP2/ERF基因家族成员,该110个蛋白质所含氨基酸大小为173～717 bp,均不存在信号肽,为非分泌蛋白。将油橄榄AP2/ERF与模式植物拟南芥AP2/ERF蛋白构建进化树发现,油橄榄AP2/ERF蛋白分为AP2、RAV、ERF和Solosist 4大类,其中ERF分为ERF和DREB两个亚类,ERF包含ERF B1～ERF B6 6个子亚类,DREB包含DREB A1～DREB A6 6个子亚类,这与模式植物拟南芥AP2/ERF的分类一致,每个亚家族同时包含了油橄榄和拟南芥AP2/ERF蛋白,说明拟南芥和油橄榄的AP2/ERF家族在进化水平上有一定的相似性。(2)油橄榄AP2/ERF同一亚家族蛋白具有相同的基因结构及保守元件,结合基因表达和进化树中已知水分调控功能的基因初步推测OeAP2-75、OeAP2-97、OeAP2-101、OeAP2-23、 OeAP2-13与油橄榄水分调节密切相关,OeAP2-13、OeAP2-28、OeAP2-104、 OeAP2-75、 OeAP2-80、OeAP2-50在2个品种中都有不同的表达量,推测可能是‘佛奥''和‘TYZ-1号''抗性不同的原因。(3)利用RT-qPCR技术检测油橄榄AP2/ERF基因在不同胁迫下的表达变化,OeAP2-101、OeAP2-28、OeAP2-42受水胁迫诱导显著上调,这与转录组分析结果一致。该研究结果为油橄榄AP2/ERF家族基因的抗逆性表达及基因功能研究奠定了基础,为油橄榄选育抗旱和耐涝砧木品种提供了方法和理论依据。     

ERF转录因子调控生物胁迫反应的研究进展

ERF家族是植物所特有的APETALA2/乙烯响应因子(APETALA2/ethylene-responsive factor,AP2/ERF)转录因子家族的一个主要亚家族,其成员结构特点是仅含有1个58-60个氨基酸组成的AP2/ERF结构域。有关该家族成员的大多数研究集中在与寒、旱等非生物胁迫方面,最近越来越多的研究表明ERF在植物抵御病虫侵害等生物胁迫方面也发挥着重要作用。ERF亚家族成员通过结合下游互作基因启动子区域的GCC box元件,从而激活或抑制这些病程相关基因的表达。同时ERFs参与水杨酸(salicylic acid,SA)、茉莉酸(jasmonic acid,JA)、乙烯(ethylene,ET)及过氧化氢(H2O2)等多种激素的信号途径,通过相互促进/拮抗高效协调体内不同激素抵御病原菌的入侵,提高植物的抗病、抗虫性。本文综述了ERF转录因子的结构功能特点、在不同植物抗生物胁迫中的调控方式,以及其通过协调不同激素信号途径相互作用来提高植物抗病虫的最新研究进展,并对其应用前景进行了展望。     

AP2/EREBP 转录因子在植物发育和胁迫应答中的作用

AP2/EREBP (APETALA2/ethylene-res ponsive element binding proteins) 是一个起源古老的转录因子超家族, 它含有1个或2个由约60-70个氨基酸残基组成的非常保守的DNA结合域 (DNA-binding domain), 即AP2/ERF结构域。根据AP2/ERF结构域的数目, AP2/EREBP转录因子可以分为2个亚族: EREBP亚族(具有1个AP2/ERF结构域)和AP2亚族(具有2个AP2/ERF结构域)。AP2亚族转录因子有调控花、胚珠和种子发育的功能, 而EREBP亚族转录因子(包括DREB类和ERF类) 的主要功能是调节植物对激素(乙烯和ABA等)、病原和胁迫(低温、干旱及高盐)等的应答反应。本文讨论了AP2/EREBP转录因子在植物发育和胁迫应答中的研究进展。     

枣AP2/EREBP转录因子的全基因组鉴定及生物信息学分析

AP2/EREBP类转录因子为植物所特有的一类转录因子,参与了植物的发育和胁迫途径。本研究利用枣的全基因组数据库,通过生物信息学手段鉴定枣的AP2/EREBP家族的全部成员,并分析AP2/EREBP转录因子家族保守域特点与功能及表达情况,为枣的遗传改良抗性育种提供理论依据。利用在线软件EBI、TAIR、Hmmer3.0、SMART、Pfam、MEME、ProtParam、SOPMA、SignalP4.1Sever和String网站,以及ClustalX2、BioEdit、MEGA6.0、MapInspect软件,对枣的AP2/EREBP转录因子家族的类型、结构域、系谱进化、序列元件、蛋白理化性质及二级结构、信号肽分析、基因定位、基因芯片表达和蛋白功能联系预测进行了分析。分析表明,枣全基因组中有127条AP2/EREBP转录因子,根据其结构域划分为4个亚家族:AP2、RAV、DREB及ERF,各家族成员数分别为:17、6、50及54。进行多序列比对和系谱进化分析,发现2个亚家族各分为6个亚族。生物信息学分析表明,AP2/EREBP家族蛋白均为亲水性蛋白、非分泌型蛋白,且该家族蛋白富含酸性氨基酸;基因定位表明,该家族基因在12条染色体上呈不均匀分布,有染色存在串联复制现象;对127个枣AP2/EREBP蛋白之间的功能联系网络进行了系统预测,分析预测了相关基因的功能和多个基因间的联系;基因芯片表达表明,AP2/EREBP转录因子在低温、高温、光照、遮阴、紫外线、UV敏感处理及外源ABA处理条件下均能被诱导表达,且在根、叶、花和果4个器官中均有显著性表达。枣AP2/EREBP基因家族结构高度保守,对外界胁迫有显著性表达,可能在枣的生长发育及其在枣的多种生理反应信号转导中发挥着重要作用,且各家族成员之间在植物生长发育过程和胁迫响应及激素信号转导途径具有交叉作用。     

龙眼ERF家族成员鉴定及其在体胚发生早期的表达

为了解龙眼ERF家族的基本性质以及在体胚发生早期的表达规律,该研究对龙眼基因组鉴定出108个ERF家族成员进行基本理化性质分析与进化树构建,并结合龙眼转录组及miRNA、lncRNA数据库对DlERF家族成员与lncRNA、miRNA之间的关系预测与表达模式分析。结果显示：(1)理化性质及进化树分析表明,AP2/ERF结构域在龙眼中相对最为保守,DlERF对龙眼的抗胁迫能力以及对病菌的防御能力可能起着重要作用。(2)RNA Seq中的表达量分析可知,95个ERF基因在转录组中检测到表达,在愈伤组织(EC)、不完全胚性紧实结构(ICpEC)与球形胚(GE)阶段中分别存在31、11与53个ERF基因高表达。(3)qRT PCR结果显示,在龙眼体胚发生早期显著表达5个ERF基因中,Dlo_008317.1ERF15 2、Dlo_022310.1ERF1 1与Dlo_009939.1ERF98在GE阶段表达量显著高于EC与ICpEC阶段,Dlo_009070.1ERF106 3与Dlo_022634.1ERF22 1则分别在EC与ICpEC阶段高表达。(4)DlERF家族成员与相关lncRNA、miRNA表达量分析显示,LTCONS_00013739对靶基因Dlo_008317.1ERF15 2为正调控关系;Dlo_miR413与Dlo_miR1510a共同靶向Dlo_009070.1ERF106 3,从EC到GE阶段均表现出显著负调控关系,而Dlo_miR413与Dlo_008317.1ERF15 2的表达量表现出正相关,推测相比于Dlo_008317.1ERF15 2,Dlo_miR413更倾向于调控Dlo_009070.1ERF106 3;调控Dlo_009939.1ERF98相关的Dlo_miR408与Dlo_miR774b,从表达趋势来看,Dlo_miR774b在龙眼体胚发生早期过程能够负调控Dlo_009939.1ERF98;Dlo_miR399c与Dlo_022310.1ERF1 1在EC到GE阶段可能存在负调控关系。(5)不同梯度浓度的乙烯处理使得DlERF基因表达显著下调。这些发现提供了有关龙眼体胚发生早期DlERF的重要见解,从而为将来对龙眼体胚发生早期过程中DlERF的功能分析奠定了基础。     

雷公藤AP2/ERF转录因子基因的克隆与分析

该研究以雷公藤发状根为材料,根据雷公藤根转录组数据设计引物,采用RT-PCR方法克隆得到2个雷公藤AP2/ERF转录因子,分别命名为TwAP2/ERF1基因(GenBank登录号:GAVZ01042389.1)和TwAP2/ERF2基因(GenBank登录号:GAVZ01016765.1)。TwAP2/ERF1基因含有一个525bp开放阅读框(ORF),编码186个氨基酸;TwAP2/ERF2基因的ORF为789bp,编码262个氨基酸;2个基因编码的蛋白质均为亲水性蛋白质。系统进化分析表明,TwAP2/ERF1与油桐(Vernicia fordii)AP2/ERF(APQ47444.1)和木油桐(Vernicia montana)AP2/ERF(APQ47365.1)相似性较高,TwAP2/ERF2与毛果杨(Populus trichocarpa)AP2/ERF(XP_002304640.1)和樱桃(Prunus pseudocerasus)AP2/ERF(ALD84477.1)相似性较高。雷公藤发状根经MeJA诱导后,TwAP2/ERF1基因的相对表达量明显提高,并于处理后9h达到最高值,为对照表达量的16.77倍;而MeJA处理对TwAP2/ERF2基因的表达表现出抑制作用,但于处理后48h相对表达量有所提高。研究表明,雷公藤TwAP2/ERF1转录因子响应MeJA早期诱导正调控,推测其可能参与调控雷公藤植物次生代谢产物的生物合成,该研究结果为阐明雷公藤次生代谢物质的生物合成调控与利用现代生物技术提高雷公藤植物细胞中次生代谢物质的含量奠定了基础。     

云杉腐烂病菌Cytospora piceae全基因组测序及比较基因组分析

【目的】壳囊孢属(Cytospora)真菌引起的林木腐烂病和枝枯病,是一类重要的、分布广泛的枝干病害,可引起重大经济损失和生态破坏。通过全基因组测序和比较基因组学分析,探究不同腐烂病菌的全基因组特征,分析其与寄主选择和致病性相关的基因或基因家族的差异性,将有助于进一步揭示腐烂病菌与寄主互作的分子机制,为腐烂病的有效防治提供基础资料。【方法】采用PacBio测序技术对云杉腐烂病菌Cytosporapiceae进行了全基因组测序和组装,并通过比较基因组学方法,从基因组水平探究引起腐烂病的4种腐烂病菌的基因组的差异,分析其共有的和特有的与致病相关的基因家族。【结果】C. piceae基因组大小为39.25 Mb,GC含量为51.79%。基于单拷贝直系同源基因构建的系统发育树显示,C. piceae与Cytospora chrysosperma的进化关系相近,Valsamali和Valsapyri则更相近。比较基因组学分析表明,4种腐烂病菌均具有重复诱导的点突变(RIP)活性,其中,C. piceae的RIP活性最强。与其他3种腐烂病菌相似,与木质素降解相关的AA3和AA7家族在C. piceae中显著扩张,但木质素降解关键酶AA5家族均缺失;C. piceae和C. chrysosperma基因组中果胶降解关键酶GH28和CE8家族基因的数量与V. mali和V. pyri相近。4种腐烂病菌都含有较多数量的MFS(major facilitator superfamily)超家族转运蛋白和较少的ABC(ATP-binding cassette transporter)超家族转运蛋白,但C. piceae含有更多DHA2、PDR和MDR类转运蛋白。4种腐烂病菌的分泌蛋白的GO分类分子功能主要集中在水解酶活性,其中V. mali含有最多数量的该类别基因;而生物学过程则集中在碳水化合物代谢过程、果胶分解过程和氧化还原过程。在次生代谢核心基因中,C.piceae的PKS基因明显少于V.mali和V.pyri;在C.piceae含有的4个特异性次生代谢基因中,3个为NRPS基因。【结论】4种腐烂病菌含有的碳水化合物活性酶的种类和数量相似,且都具有较强的果胶降解能力。不同腐烂病菌的膜转运蛋白中多药转运体的选择性扩增,以及次生代谢核心基因中NRPS类基因的特异性存在和缺失,表明它们作为重要的致病因子很可能在腐烂病菌寄主选择中发挥了重要的作用。     

麻竹miR172a靶基因DlAP2的克隆及其表达

为了解开花麻竹(Dendrocalamus latiflorus)的Dl AP2基因功能,采用RT-PCR和RACE技术克隆了mi R172a靶基因AP2同源序列c DNA全长,命名为Dl AP2。结果表明,Dl AP2基因c DNA全长为1729 bp,包含5′端非编码区81 bp、开放阅读框1464 bp、3′端非编码区160 bp和24个碱基的Poly A尾巴,在编码框靠近3′端130 bp处有1个高度匹配mi R172a的结合位点(CTGCAGCATCATCAGGATTCT)。Dl AP2编码487个氨基酸的蛋白,具有两个AP2结构域,属于AP2/ERF家族AP2亚家族的AP2组,与来自其它单子叶植物的AP2蛋白均有较高同源性。RLM-5′RACE分析表明,mi R172a主要在靶序列的第11~12个碱基之间剪切靶基因Dl AP2的mi RNA。q RT-PCR结果表明,麻竹花芽中Dl AP2基因的表达规律与mi R172a表达变化正好相反,证明mi R172a对Dl AP2基因的表达具有调控作用。     

旱麦草EtAP2基因的克隆及其对干旱胁迫的响应表达

以旱麦草(Eremopyrum triticeum)为实验材料,利用RT-PCR技术从旱麦草叶片中克隆了1个AP2/ERF家族基因,命名为EtAP2(GenBank登录号KX622583)。EtAP2基因含有1 128bp开放阅读框,编码375个氨基酸,相对分子质量40.87kD,等电点为5.36。多序列比对和进化树分析表明,该基因编码蛋白具有2个AP2保守结构域,与小麦AP2/ERF家族蛋白具有较近的亲缘关系。实时荧光定量PCR分析表明,15%PEG 6000模拟干旱胁迫可诱导EtAP2基因在根和叶中表达,且在根中对干旱胁迫的响应大于叶片。研究表明,EtAP2可能参与旱麦草对干旱逆境胁迫应答的调节。     

互花米草NHX2基因的克隆与功能鉴定

NHX2属于CPA1基因家族,编码Na~+/H~+逆向转运蛋白,控制液泡膜中活性K~+的摄取,同时调节气孔的关闭。该研究以耐盐植物互花米草为材料,采用PCR技术克隆NHX2基因,并将其转入拟南芥进行相关功能鉴定。结果显示:(1)成功克隆获得互花米草NHX2基因CDS序列(1 602 bp),命名为SaNHX2,该基因编码533个氨基酸,SaNHX2蛋白的分子量约为58.65 kD,定位于细胞核和细胞膜,表明SaNHX2基因可能发挥转录调控的功能。(2) qRT-PCR结果显示,在ABA、NaCl和干旱胁迫处理下,互花米草叶和根中SaNHX2基因的表达量均上调。(3)为进一步鉴定其功能,成功构建植物表达载体,将SaNHX2基因转入拟南芥;经RT-PCR检测结果显示,SaNHX2基因在转基因植株中过表达;高盐胁迫处理后,转SaNHX2基因拟南芥的主根长度、叶绿素总量和相关胁迫应答基因表达量均高于转空载拟南芥,表明转SaNHX2基因拟南芥的耐盐能力显著增强。研究表明,SaNHX2基因可能在盐胁迫调节机制中发挥调控作用,可作为改良农作物耐盐的重要候选基因。     

沙棘HrTCP转录因子家族鉴定及其干旱胁迫下的表达分析

TCP家族是植物特有的响应高盐、干旱等非生物胁迫的重要转录因子。该研究基于沙棘转录组数据,利用生物信息学与qRT-PCR对HrTCP转录因子家族进行鉴定,预测其家族成员的结构和功能,为解析TCP转录因子调控沙棘抵御干旱胁迫的作用机制奠定基础。结果表明:(1)获得了11个HrTCP转录因子成员,并命名为HrTCP2/4/7/8/11/13/15/17/18/19/20,编码氨基酸序列长度在218～590之间,蛋白质相对分子量为23.44～61.78 kD;亚细胞定位预测发现,除HrTCP13/17/18蛋白定位于细胞质,其余8个蛋白均定位于细胞核。(2)在干旱(15%PEG-6000)和高盐(200 mmol/L NaCl)胁迫后HrTCP4/7/19/20基因表达量呈不同程度上升趋势,其中HrTCP20表达量极显著高于对照,分别是对照的24倍与23倍。(3)外源激素脱落酸(0.1 mmol/L ABA)和茉莉酸甲酯(0.1 mmol/L MeJA)处理后,HrTCP7/19/20基因表达量也均呈上升趋势,其中,ABA诱导下HrTCP19基因表达量达到最高,是对照的16倍,而MeJA诱导下HrTCP20基因表达量上升最高,是对照的5倍。研究发现,HrTCPs转录因子家族成员可受干旱、高盐和激素诱导表达,进而调控沙棘对干旱胁迫的响应。     

蒺藜苜蓿LBD转录因子基因家族全基因组分析

LBD是植物中所特有的转录因子基因家族,在调控植物侧生组织发育、营养代谢以及响应逆境胁迫等方面具有重要作用。该研究利用生物信息学手段,从全基因组水平筛选和鉴定了蒺藜苜蓿LBD基因家族,并对基因结构、系统进化、进化压力、保守域、染色体定位以及基因表达模式等进行了分析。研究结果共鉴定出2类5亚类共计56个蒺藜苜蓿LBD家族基因,在8条染色体上均有分布,但分布不均匀。该家族成员外显子数目都不超过2个,结构简单,基因间在进化时存在负向选择作用。基因表达模式分析发现,该家族成员的表达具有一定的时空特异性,并受干旱和氮素调控。该研究结果对蒺藜苜蓿LBD基因功能研究及进化分析具有重要的意义。