花生FAD基因家族的全基因组鉴定与表达模式分析

脂肪酸去饱和酶(FAD)是合成不饱和脂肪酸的关键酶,在植物维持各项生命活动、应对生物与非生物胁迫方面具有重要作用。本研究对花生基因组中的FAD基因家族进行鉴定与生物信息学分析,并利用花生转录组数据对花生FAD(Ah FAD)基因的表达模式及可变剪接形式进行分析。结果显示,在花生基因组中共有31个Ah FAD基因,分为6类,Ah SAD基因亚家族有12条序列,AhFAD2有6条序列,Ah FAD3有4条序列,Ah FAD4/5有3条序列,Ah FAD6有2条序列,Ah FAD7/8有4条序列。聚类分析表明,Ah FAD有两大分支:游离Ah SAD分支和膜结合Ah FAD分支;膜结合Ah FAD分支中Ah FAD4/5起源最早,ω-3 Ah FAD由ω-6Ah FAD进化而来;单双子叶聚类于不同分支。表达模式分析表明,Ah SAD、AhFAD2和Ah FAD3在种子中表达水平较高,Ah FAD4/5、Ah FAD6、Ah FAD7/8在叶中表达量最高。可变剪接分析表明,Ah FAD亚家族之间可变剪接形式差异明显,其中有12个Ah FAD基因发生了可变剪接。TTS和TSS是发生最多的可变剪接事件,其次是IR,最少的是AE和ES。     

花生DGAT基因家族的生物信息学分析

含油量是花生最重要的品质性状之一。花生油的主要成分是三酰甘油(TAG),二酰甘油酰基转移酶(DGAT)是催化合成TAG的关键酶类,其活性高低与含油量呈正相关。本研究利用生物信息学方法对花生DGAT(Ah DGAT)基因家族的进化、基因结构特点、表达模式以及可变剪接等进行了分析。进化分析结果显示,Ah DGAT1、Ah DGAT2和Ah DGAT3亲缘关系较近,与Ah WSD/DGAT的亲缘关系稍远。Ah DGAT1具有17~18个外显子,Ah DGAT2具有7~9个外显子,Ah DGAT3具有2~3个外显子,而Ah WSD/DGAT具有3~6个外显子。Ah DGAT1的跨膜区数目为7~11个,Ah DGAT2和Ah WSD/DGAT为0~2个,Ah DGAT3没有跨膜结构。另外,4个亚家族均具有自己特有的保守结构域和表达模式。Ah DGAT1和Ah DGAT2在种子中表达量较高,而Ah DGAT3和Ah WSD/DGAT在根和叶中的表达量较高。Ah DGAT1的可变剪接形式最多,其次为Ah WSD/DGAT。花生的4个Ah DGAT亚家族具有不同的特点,表明它们可能在进化中具有不同的祖先。     

基于RNA-seq数据分析玉米抗虫响应基因的可变剪接事件

采用Illumina Hiseq~(TM)分别对黏虫取食的玉米叶片及对照组材料进行转录组测序,鉴定和分析响应黏虫取食基因的可变剪接事件。结果表明,对照组中鉴定出15 701个基因对应79 363个可变剪接事件,黏虫取食组中鉴定到11 791个基因的39 385个可变剪接事件。2组玉米基因组在不同的可变剪接类型中,均是以第1个外显子可变剪切、最后1个外显子可变剪切、单内含子滞留和可变5′或3′端剪切4种类型为主。对2组发生可变剪接事件的基因进行比对发现,黏虫处理组新增加了1 121个可变剪接基因,差异表达分析鉴定出177个表达显著差异的可变剪接基因。GO功能富集分析表明,发生可变剪接的差异基因主要富集于氧化还原酶活性、转移酶活性、DNA结合、ATP结合、代谢过程、以DNA为模板的转录及调控相关的功能,表明这些可变剪接差异表达基因参与了玉米的抗虫响应过程。     

草鱼CC趋化因子基因CCL24及其可变剪接

CC趋化因子是一类能够促进动物体内炎症部位的各种白细胞的补充、激活和黏附的趋化性细胞因子家族,是鱼类天然免疫系统的重要组成部分。分析了从草鱼肠道cDNA文库中筛选到的CC趋化因子基因CCL24,并克隆了阅读框区域内的基因组序列。序列分析表明,CCL24基因由4个外显子和3个内含子组成;外显子拼接的序列与CCL24cDNA序列完全一致,编码95个氨基酸,有2个相邻的半胱氨酸(CC),为典型的CC趋化因子亚家族成员。此外,还发现草鱼CCL24基因存在可变剪接现象,第1内含子没有被剪切掉的非正常转录本翻译后可能产生没有CC趋化因子活性的24个氨基酸的短肽,而且这种转录本在精巢、头肾、皮肤、肝胰脏、肠道、肌肉、肾脏等组织均检测到;而正常剪接的CCL24转录本,仅在肾脏和肠道中检测到,且其表达量要明显低于非正常剪接的转录本。     

草鱼CC趋化因子基因CCL24及其可变剪接

CC趋化因子是一类能够促进动物体内炎症部位的各种白细胞的补充、激活和黏附的趋化性细胞因子家族,是鱼类天然免疫系统的重要组成部分。分析了从草鱼肠道cDNA文库中筛选到的CC趋化因子基因CCL24,并克隆了阅读框区域内的基因组序列。序列分析表明,CCL24基因由4个外显子和3个内含子组成;外显子拼接的序列与CCL24cDNA序列完全一致,编码95个氨基酸,有2个相邻的半胱氨酸(CC),为典型的CC趋化因子亚家族成员。此外,还发现草鱼CCL24基因存在可变剪接现象,第1内含子没有被剪切掉的非正常转录本翻译后可能产生没有CC趋化因子活性的24个氨基酸的短肽,而且这种转录本在精巢、头肾、皮肤、肝胰脏、肠道、肌肉、肾脏等组织均检测到;而正常剪接的CCL24转录本,仅在肾脏和肠道中检测到,且其表达量要明显低于非正常剪接的转录本。     

辣椒ARF基因家族的鉴定与表达分析（英文）

为进一步研究辣椒ARF (Auxin Response Factor)基因家族在其生长发育及逆境胁迫中的作用,利用生物信息学方法从辣椒基因组中鉴定出20个ARF基因,这些基因不均匀地分布在辣椒12条染色体上。进化关系显示辣椒ARF基因可分为Class I和II两大类,进一步可细分为Ia、Ib、Ic、Id、Ie、IIa和IIb等7个亚类。基因结构图表明,该基因家族由1～15个外显子构成,且各基因在不同发育时期不同组织中的表达量不同,具有一定的特异性。qRT-PCR结果表明,高盐胁迫可以明显激活或抑制ARF基因家族的表达。     

花生KNOX基因家族的全基因组鉴定及组织表达分析

为研究花生KNOX基因家族在生长发育以及非生物胁迫响应中的功能，本研究采用生物信息学手段在栽培种花生基因组水平上对KNOX基因家族成员进行鉴定，并对其理化性质、基因结构、蛋白质保守结构域、系统进化、启动子顺式作用元件以及组织与逆境下的表达模式进行分析。结果显示：花生基因组中共有26个KNOX基因家族成员，分布在17条染色体上，绝大多数编码酸性蛋白，均为亲水性蛋白，并全部定位于细胞核。系统发育分析可将花生KNOX家族分为ClassⅠ和ClassⅡ两个亚家族，并进一步分为4个亚组。启动子分析发现一些与生长发育、激素和胁迫响应相关的顺式作用调控元件。全生育期组织表达分析结果显示，ClassⅠ类基因表达组织特异性明显，主要在花、营养茎尖、生殖茎尖、果针等组织中高表达；ClassⅡ类基因时空表达更广泛。在7种非生物胁迫处理下，部分基因表达量变化较大，其中5个基因响应特定胁迫，RT-qPCR验证了2个特异性响应PEG胁迫的基因。本研究为进一步探究花生KNOX基因在生长发育、逆境响应中的功能奠定了基础。     

观赏海棠'火焰'果皮中可变剪接基因分析

[目的]研究观赏海棠'火焰'果皮中可变剪接基因.[方法]利用已有苹果属观赏海棠品种'火焰'果皮转录组测序数据,分析可变剪接类型和数量,并对发生可变剪接的基因进行KEGG分析,研究观赏海棠可变剪接基因与其发育过程的关系.[结果]在观赏海棠果皮发育过程中,主要发生的可变剪接是第一个外显子可变剪接和最后一个外显子可变剪接两种类型.可变剪接在果皮发育前期发生较多,同时KEGG通路分析表明可变剪接主要发生在嘌呤代谢、氧化磷酸化、糖酵解/糖异生等通路过程中.[结论]可变剪接在观赏海棠果皮的整个发育过程中普遍存在.观赏海棠果皮发育过程中共有的可变剪接基因主要与磷酸化和糖代谢有关.     

山羊5个组织的全长转录组构建与分析

【目的】构建山羊（Capra hircus）全长转录组,旨在改进参考基因组的功能注释。【方法】利用PacBio长读段测序技术对南江黄羊肝、脾、肾周脂肪、背最长肌和睾丸5个组织的混合样本进行高通量测序;对原始数据进行校正、聚类去冗余,最后用KEGG数据库对预测转录本序列进行功能注释。【结果】质控后获得28 432条高质量一致性全长序列,其中99.49%的序列被成功比对到山羊基因组;得到源自9 879个基因的19 115个转录本。进一步分析得到表达基因总数的41.39%发生了可变剪接;其中起始外显子可变剪接占比最高（36.47%）,外显子互斥占比最低（1.78%）。另外,144条序列未被比对到参考基因组上,这些序列的平均GC含量高达56.50%;利用KEGG数据库对其进行功能注释,揭示28个转录本与免疫系统相关。【结论】获得的测序数据全面涵盖混合样本中表达的转录本,证明了山羊基因组转录的复杂性。研究结果总体上提高了每个基因平均表达的转录本数量,证明可变剪接的广泛存在;未成功比对的转录本序列GC含量较高且许多转录本来自免疫基因的表达。     

简州大耳羊早期胎儿到新出生阶段背最长肌中可变剪接的鉴定与分析

可变剪接广泛存在于动植物的各个组织中,是形成动植物体内蛋白质多样性和复杂性的重要机制之一。以简州大耳羊为研究对象,采集母羊妊娠期第45天、60天、105天胎儿以及出生后第3天羔羊的背最长肌组织,进行转录组测序。测序数据经比对、组装后,共鉴定出来自于4 562个(16.77%)已注释基因的137 308个可变剪接事件。其中,所占比例最高的类型为外显子跳跃(30.64%),而3′端可变剪接位点、内含子保留、5′端可变剪接位点的比例分别为19.47%、18.07%和10.92%。对RBMs、HNRNPs和SRSFs 3类可变剪接调控因子中58个基因表达量分析表明,19个基因的表达量呈现显著差异(P0.05);其中15个基因在妊娠期第45天和第60天的表达量显著高于其他2个阶段。GO功能注释和富集分析结果表明,4个阶段均发生可变剪接的基因显著富集在细胞表达调控和物质合成等过程(P0.05)。值得注意的是,在妊娠期第60天,可变剪接事件、可变剪接基因的数量和可变剪接因子表达量均显著高于其他阶段,并且特异性基因显著富集在有丝分裂和大分子合成相关过程。上述结果表明,妊娠期第60天胎儿的成肌过程非常活跃,暗示可变剪接在山羊背最长肌发育过程中具有重要作用,为深入探索可变剪接在山羊肌肉发育中的作用提供理论依据。     

陆地棉扩展蛋白基因的鉴定与特征分析

目的 扩展蛋白（Expansin）是细胞壁的重要组成部分,在植物的生长发育及逆境胁迫应答等方面均发挥着重要作用。基于全基因组水平系统鉴定陆地棉Expansin基因家族,并通过生物信息学及表达模式分析,为揭示扩展蛋白基因在棉花生长发育中的功能及后续利用奠定基础。方法 利用BLAST和HMMER在陆地棉基因组中搜索并鉴定扩展蛋白基因家族成员;利用ClustalW、MEGA、MCScanX、Prot Param、MEME、SignalP、Euk-mPLoc、FancyGene和DnaSP等软件对其基因序列和蛋白序列进行生物信息学分析。通过RNA-seq数据分析扩展蛋白基因的表达模式和部分同源基因间表达差异,利用qRT-PCR验证部分扩展蛋白基因的表达谱。结果 陆地棉基因组中含有46个EXPA基因、8个EXPB基因、6个EXLA基因和12个EXLB基因,合计72个Expansin基因;四倍体陆地棉中扩展蛋白成员的数量几乎是二倍体棉种（亚洲棉与雷蒙德氏棉）的2倍。除GhA02和GhD06 2个染色体外,其余各染色体上均分布有数目不等的扩展蛋白基因（2—4个）,具有部分同源关系的染色体GhA08和GhD08分别有5个和8个扩展蛋白基因。系统发育树显示,各亚家族成员聚集成群,并且大部分的末端分支均由来源于3个物种的4个（亚）基因组的4个基因组成,如EXPA亚家族的Cotton_A_28454/Gh_A03G0885/Gh_D02G1269/Gorai. 005G142200等等,4个基因之间具有同线性关系。亚细胞定位发现陆地棉所有的扩展蛋白均位于细胞外。基因结构分析显示,扩展蛋白基因由3—5个外显子组成,外显子-内含子结构在进化上高度保守且与氨基酸序列的多样性一致,且在外显子上存在密码子偏好性。RNA-seq数据显示,不同基因在不同时空条件下存在特异性表达,如GhEXPA19A和GhEXPA19D相比其他基因在纤维10 DPA和20 DPA中的表达量很高;在不同的组织（如子叶、新叶、老叶、苞叶）中,GhEXPA24D具有较高的表达量。部分同源基因之间具有不同的表达模式,显示它们之间功能的异化与互补。qRT-PCR结果与RNA-seq数据基本吻合,如GhEXLA3A和GhEXLA3D在纤维发育的伸长阶段高量表达。GhEXPA19D和GhEXLA2D在3DPA的胚珠中表达活跃。结论 陆地棉基因组中含72个扩展蛋白基因,其在DNA水平和氨基酸水平具有一致的结构多样性和进化保守性,在转录水平具有各异的表达模式,显示出家族内成员间功能上的异化与互补。     

荆芥HD-Zip基因家族的全基因组鉴定及分析

  目的  鉴定荆芥Schizonepeta tenuifolia的HD-Zip基因家族，利用生物信息学方法分析其在全基因组中的分布和相关特征以及在不同时期中的表达规律，为该家族基因的进一步研究奠定基础。  方法  根据已经表征的HD-Zip基因，筛选荆芥基因组内的HD-Zip基因序列，利用MEME、PlantCARE、NCBI、MEGA X、MCScanX、Circos等在线网站及软件对蛋白序列进行基本理化性质分析、进化树构建、染色体定位、基因结构分析、共线性基因分析等。  结果  在荆芥全基因组中共鉴定到42条HD-Zip基因序列，它们可被分为4个亚家族，分别含有16、7、5、14个基因，亚家族之间的基因长度、结构及保守基序差异显著，但在亚家族内部保守，荆芥基因组与拟南芥Arabidopsis thaliana基因组共线性分析发现有37对基因，可能具有相似的生物学功能。荆芥的4个亚家族基因的顺式元件中均高频出现了光响应、脱落酸响应、MeJA响应等元件，在不同生长时期的叶片及部位的转录组数据中具有不同的表达趋势，Ⅰ亚家族主要在幼叶中表达，Ⅱ和Ⅲ亚家族主要在根中在表达，Ⅳ亚家族主要在叶中表达。  结论  在荆芥基因组中共获得42条HD-Zip基因序列，被分为4个亚家族(HD-ZipⅠ~Ⅳ)，亚家族内部高度保守，亚家族之间差异显著，其基因结构、保守结构域及表达模式不同。亚家族Ⅰ和Ⅱ，亚家族Ⅲ和Ⅳ亲缘关系更近，HD-Zip基因具有组织表达差异性，协同调控了荆芥的生长发育和次生代谢。图9表1参25     

花生热激蛋白AhHSP70与热激因子AhHSF基因的克隆及表达分析

热激蛋白(heat shock proteins,HSPs)和热激转录因子(heat shock factors,HSFs)在植物热胁迫信号转导和耐热性的产生过程中发挥了重要作用。本研究从花生转录组文库中筛选到HSP70、HSF的c DNA片段,通过序列比对在花生全基因组序列中获得这两个基因的基因组序列,根据序列信息设计引物,以花生叶片c DNA为模板扩增全长ORF并进行生物信息学分析。结果显示,Ah HSP70基因的ORF全长为1 962 bp,编码653个氨基酸,分子质量为71.45 k D,理论等电点p I为4.93;Ah HSF基因的ORF全长为1 212 bp,编码403个氨基酸,分子质量为46.03 k D,理论等电点p I为4.85。Ah HSP70与Ah HSF均不具有信号肽,为可溶性蛋白,二级结构中有大量无规则卷曲。利用这两个基因的氨基酸序列分别与来源于其他物种HSP70、HSF的氨基酸序列进行同源比对,并构建进化树进行亲缘关系分析,结果表明,Ah HSP70与大豆Gm HSP70亲缘关系较近,与番茄Sl HSP70亲缘关系比较远;Ah HSF与菜豆Pa HSF亲缘关系较近,而与蒺藜苜蓿Mt HSF的亲缘关系较远。利用实时荧光定量PCR对Ah HSP70和Ah HSF在热胁迫情况下的表达进行分析,结果表明这两个基因在42℃高温条件下表达量显著升高,Ah HSP70在高温胁迫3 h后表达明显升高,热处理24 h和48 h后,表达量为对照的50倍和135倍;转录因子Ah HSF在高温胁迫3 h后表达明显升高,高温处理6 h后达到最高,随后下降。本研究初步验证了花生热激蛋白和热激因子基因在花生响应高温胁迫中的作用。     

哈萨克羊和藏绵羊尾部脂肪组织差异表达基因和可变剪接的鉴定与分析

为鉴定哈萨克羊和藏绵羊尾部脂肪组织的差异表达基因和可变剪接,分别提取6只(3公,3母)无亲缘关系的2岁成年哈萨克羊与藏绵羊尾部脂肪组织总RNA,构建2个高通量测序文库,利用TopHat2软件将短序列定位到绵羊参考基因组,采用cuffdiff软件分析基因表达丰度和差异表达基因。结果表明,共有19个基因在2个品种绵羊尾脂中的表达丰度均高于1 000FPKM,包括ADIPOQ、SPARC、VIM等。其中,2个品种尾部脂肪组织中表达量最高的转录本均定位于线粒体基因组5 331~14 154bp,该区域主要包含COX1、COX2、ATP8等基因,其表达丰度在哈萨克羊和藏绵羊中分别高达29 724.20和28 818.80FPKM。而FABP4则是表达丰度最高的核基因,在哈萨克羊和藏绵羊中的表达量分别为3 194.03和2 667.09。差异表达基因分析共鉴定出627个差异表达的基因,与藏绵羊尾部脂肪组织的表达量相比,哈萨克羊的尾部脂肪组织中共有221个上调基因和406个下调基因。其中,HBB基因表达量的变化倍数最大,下调约98倍。差异表达基因的富集分析结果表明,哈萨克羊尾脂中上调表达的基因主要参与脂肪酸代谢过程。另外,在2个品种尾脂中鉴定出11 870个可变剪接事件,其中383个为差异剪接事件。对相关差异表达基因和可变剪接的进一步研究将有助于揭示绵羊肥脂尾形成的分子机制。     

辣椒紫色酸性磷酸酶基因家族的鉴定与表达特征

紫色酸性磷酸酶(purple acid phosphatases,PAPs)属于金属磷酸酯酶家族,参与了植物大量的生理反应,特别是磷素的吸收与利用。本研究利用已经发表的辣椒基因组数据,采用同源比对方法对辣椒PAP基因家族进行鉴定与表达分析。通过鉴定发现辣椒基因组中至少含有25个PAP基因成员,氨基酸序列长度在264~639 aa,具有1~12个内含子,不均匀地分布在12条染色体上,仅发现1对基因以串联重复的形式存在。系统发育关系表明,辣椒、拟南芥和水稻PAP基因家族成员可以分为3个家族,8个亚家族,每个家族成员的数目是变异的,而且存在6对直系同源基因和20对旁系同源基因,这表明PAP基因家族的存在早于辣椒、拟南芥和水稻的分化。基于RNA-seq数据库和q RT-PCR方法,对辣椒PAP基因的组织表达及果实发育的不同时期进行分析,结果表明PAP家族基因具有组织表达特异性和发育时期表达特异性。     

甜瓜GRAS家族全基因组的鉴定与表达分析

为了研究甜瓜(Cucumis melo L.)GRAS基因家族的结构特征、进化模式和生长发育过程中的表达模式，对甜瓜GRAS全基因组的基本信息、染色体位置、共线性、系统发育树、基因结构、表达模式、顺式元件进行了全面分析。结果表明，从甜瓜中鉴定出37个GRAS基因；这些基因在甜瓜的11条染色体上呈不均匀分布，并在6、7、8号染色体上出现串联重复；甜瓜GRAS家族可分为HAM、LS、SCR、SCL3、DELLA、LISCL、SHR、PAT1 8个亚族；基因结构分析显示，含有1～3个内含子的基因有8个，没有内含子的基因有29个；启动子分析，鉴定到生长素、脱落酸、赤霉素、低温、干旱、光信号响应作用元件以及参与分生组织和胚乳表达的元件；表达模式分析结果表明，HAM、DELLA亚族的基因几乎在所有组织中都有表达，LS亚族在根和茎中表达含量较高，PAT1亚族在甜瓜的整个生长发育过程中发挥着重要作用。得出结论，甜瓜全基因组共获得GRAS基因37个，共可分为8个亚族，不同的亚族具有不同的基因结构、蛋白质结构和表达模式；不同亚族在甜瓜生长发育的不同阶段具有差异性表达，它们共同参与调控甜瓜的发育过程。     

马铃薯ZF-HD基因家族的鉴定与分析

为了进一步研究ZF-HD基因家族在马铃薯中作用,本研究利用各类基因家族分析软件对马铃薯ZF-HD家族成员进行全基因组鉴定与分析。结果表明,马铃薯ZF-HD家族共包含20个成员,大部分基因都只含有2个内含子。19个马铃薯ZF-HD成员分布于6条染色体上,有三对重复基因,其中两对为重复串联基因,另一对为片段重复基因。根据进化关系和保守结构域分布可以分为ZHD和MIF和UN三个亚家族,同一亚家族成员之间motifs特征相似。组织差异表达表明,几乎所有的ZF-HD在马铃薯中都在叶片中表达。在非生物胁迫中,在干旱(利用甘露醇模拟)和高盐胁迫下马铃薯ZF-HD基因均高表达。本研究进一步揭示了ZF-HD在马铃薯中的保守性和参与马铃薯生长与抗逆的潜在功能,为后续研究奠定基础。     

香蕉钾转运体HAK/KUP/KT家族鉴定及其在果实发育和低钾胁迫下的表达分析

钾是香蕉生长发育中需求量最大的营养元素,对果实产量、品质和贮藏性都有重要的影响。从香蕉基因组中鉴定出24个HAK/KUP/KT(MaKUP)基因,并对其理化性质、染色体定位、基因结构、保守功能域、进化关系、组织表达和低钾胁迫下的表达特征进行分析。结果表明,MaKUP编码的肽链平均有772个氨基酸,分子量为63.54～93.82 ku,等电点为5.44～9.30,蛋白不稳定指数为32.69～44.29,脂肪族指数为105.34～116.90,总平均亲水性为0.352～0.549,含有8～11个外显子,可分为4个亚族;表达分析发现15个MaKUP基因家族成员在香蕉不同发育时期的果实中高表达,6个MaKUP基因家族成员在低钾胁迫中上调表达,表明MaKUP在香蕉果实发育和钾吸收中具有重要的作用。     

杨树磷酸肌醇特异性磷脂酶C基因家族鉴定与分析

为了对已测序杨树Populus trichocarpa Torr.Gray中磷酸肌醇特异性磷脂酶C(phosphoinositidespecific phospholipase C,PI-PLC)家族基因进行系统分析,利用杨树基因组数据库,通过生物信息学手段,鉴定杨树PI-PLC家族基因的基因结构、染色体定位和编码蛋白,通过蛋白序列比对进行进化和分类分析。结果表明,杨树中含有7个PI-PLC家族基因,含有8~10个外显子,分布于杨树的4条染色体上。MEME和Pfam保守结构域分析显示,杨树PI-PLC蛋白均含有4个保守的EF、X、Y和C2结构域。蛋白质进化树分析结果表明PI-PLC可分为2个亚家族。     

鹰嘴豆NF-Y基因家族的全基因组序列鉴定及生物信息学分析

为了揭示鹰嘴豆NF-Y基因家族的功能,从鹰嘴豆全基因组中鉴定NF-Y基因,并利用生物信息学分析软件对其进行系统发育进化、染色体定位、基因结构、保守结构域以及microRNA的靶位点预测分析。结果表明,共鉴定出41个鹰嘴豆NF-Y基因家族成员,在进化上分为3个亚家族:CarNF-YA(8个)、CarNF-YB(22个)和CarNF-YC(11个);这些基因分布于7条(Ca1—Ca7)染色体上,另外还有3条序列未定位到染色体上。鹰嘴豆NF-Y家族基因序列大部分都含有1～5个内含子和1～6个外显子,少数序列如CarNF-YB9和CarNF-YC6等仅含有外显子结构;编码的蛋白质中共发现3个保守基序(Motif1、Motif2和Motif3),其中Motif1、Motif2在CarNF-YB和CarNF-YC亚家族中均有分布,而Motif3仅分布在CarNF-YA亚家族中。鹰嘴豆和拟南芥NF-Y蛋白家族成员系统进化树表明,Car/At NF-YA可以分成6类,分别由6、1、3、3、1、4种蛋白质组成; Car/At NF-YB可以分成4类,分别由21、8、2、4种蛋白质组成; Car/At NF-YC亚基可以分成3类,分别由17、4、4个蛋白质组成。对鹰嘴豆NF-Y基因家族成员的miRNA169靶位点预测发现,41个成员中仅有CarNF-YA亚家族中的6个成员受到miRNA169调控。