缢蛏(Sinonovacula constricta)是我国传统的四大海水养殖贝类之一，养殖过程中易受病原菌感染而造成大量死亡和经济损失。Toll样受体(Toll-like receptor， TLR)是目前研究最多的模式识别受体，在无脊椎动物的先天性免疫中发挥着重要作用。研究缢蛏TLR分子特征及其免疫功能对于预防和控制病原菌感染至关重要。本研究从缢蛏基因组序列中鉴定出42个TLR (ScTLR)基因，其中，33个编码典型的TLR蛋白，其余9个为TLR样蛋白。根据蛋白结构域特点，将典型的TLR蛋白分为2大类4个亚类，一类为多半胱氨酸簇TLR (mccTLR)，包括P型TLR (1个)和sPP型TLR (7个)；另一类为单半胱氨酸簇TLR (sccTLR)，包括sP型TLR (16个)和Ls型TLR (9个)。与其他9种软体动物TLR基因的类型和数目比较结果显示，缢蛏基因组中未发现其他软体动物中存在的V型和Twin-TIR型TLR，起源古老的mccTLR类群在软体动物中没有大规模的扩张，而起源较近的sccTLR类群却发生了大规模的扩张。荧光定量PCR分析显示，6种ScTLR在检测的7种组织中普遍表达，且在血细胞、鳃和肝胰腺组织中高表达。最后，副溶血弧菌(Vibrio parahemolyticus)胁迫12 h和48 h的缢蛏鳃和肝胰腺转录组数据分析显示，副溶血弧菌感染前后9个TLR基因在鳃或肝胰腺中的表达量发生显著变化，6个基因(ctg118.25、ctg118.26、ctg356.25、ctg774.6、ctg681.6和ctg1513.5)在感染12 h或48 h后的鳃组织中显著差异表达，前5个基因表达量上调，仅ctg1513.5表达量下调；3个基因(ctg467.9、ctg2496.3和ctg903.17)在感染12 h或48 h后的肝胰腺中显著差异表达，其中，ctg467.9和ctg2496.3表达量下调，ctg903.17表达量上调。综上所述，本研究结果将为深入探讨不同TLR基因在缢蛏天然免疫中发挥的作用提供研究基础。

腺苷三磷酸结合盒转运蛋白(ATP binding cassette transporter， ABC transporter)家族是最古老的膜蛋白家族之一，广泛存在于原核生物和真核生物体内，利用ATP水解释放的能量对氨基酸、脂质、抗生素等物质进行跨膜运输，参与生物体内多种生理活动。目前，对软体动物ABC转运蛋白家族的鉴定，仅在虾夷扇贝(Patinopecten yessoensis)、栉孔扇贝(Chlamys farreri)和长牡蛎(Crassostrea gigas) 3种双壳类中有系统研究，对缢蛏(Sinonovacula constricta) ABC转运蛋白家族的鉴定和表达模式分析尚未见报道。利用缢蛏基因组和转录组数据，运用本地及NCBI在线BLAST程序、FGENESH+、SMART、ExPASy、MEGA X、Mev 4.90和MapInspect等生物信息学分析工具，在全基因组水平系统地鉴定出52个缢蛏ABC转运蛋白，并对转运蛋白基因外显子数目、染色体定位等信息进行分析；通过系统进化分析，将缢蛏ABC转运蛋白家族分为8个亚家族，即ABCA～ABCH；通过比较分析不同物种ABC转运蛋白家族，推测基因串联复制事件对软体动物ABC转运蛋白家族成员数目的增加有影响。ABC转运蛋白基因在缢蛏不同发育时期和组织的表达分析显示，ABCC和ABCG亚家族多个基因在缢蛏稚贝时期表达量达到最高值，ABC转运蛋白基因在鳃、肝胰腺与性腺中表达种类较多。软体动物作为第二大动物类群，目前仅有3个物种的ABC转运蛋白得到了系统分类，缢蛏ABC转运蛋白基因家族的系统鉴定与表达模式分析为研究软体动物ABC转运蛋白基因的进化和缢蛏ABC转运蛋白功能提供了基础资料。

[目的]对木本模式植物毛果杨PLD基因家族的进化中的选择压力、启动子中顺式作用元件、组织表达特性以及盐胁迫下表达模式进行分析,为挖掘PtrPLD在非生物胁迫中作用提供参考。[方法]利用拟南芥PLD基因家族蛋白序列比对得到毛果杨基因组同源基因,再经过保守结构域鉴定后确定PtrPLD基因;利用软件ClustalW和MEGA对PtrPLD和AtPLD基因的氨基酸序列进行比对和系统进化分析;利用MEME、PlantmPLoc、 ExPasy等软件工具分析PtrPLD基因及编码蛋白的特征;利用Tbtools软件分析同源基因的Ka/Ks值;利用Plantcare在线工具分析PtrPLD启动子中顺式作用元件;利用Phytozome转录组数据库以及qRT-PCR分析PtrPLD组织表达特性;利用qRT-PCR分析各组织中PtrPLD对盐胁迫响应情况。[结果]PtrPLD家族的16个基因可分为C2-PLD和PX/PH-PLD 2个亚家族,分别包含13个和3个基因;有7对旁系同源基因且它们之间的Ka/Ks均远小于1;PtrPLD家族基因启动子区含有大量非生物胁迫和激素响应元件,其中,PtrPLDδ4启动子共含有9种、20个元件;PtrPLD家族编码的蛋白均含有Motif 1～4,且同一进化分支上的基因编码蛋白序列高度保守。PtrPLD基因家族表达特性分析表明,PtrPLD家族基因在根、茎和叶中具有表达特异性,并且多数成员主要在根部表达;NaCl胁迫下,PtrPLD家族基因在根、茎和叶中表达量在0～72 h内均表现为先上升后下降再上升的变化趋势。[结论]PtrPLD家族基因在毛果杨响应盐胁迫过程中有着重要作用,本项研究对于今后PtrPLD家族基因生物学功能的鉴定与非生物逆境胁迫响应基因资源的挖掘具有推动作用。

【目的】基于黄瓜基因组信息和转录组测序大数据，利用生物信息学手段，对黄瓜中DIR基因家族进行鉴定，并分析其在不同组织器官和胁迫响应过程中的表达模式，为后续深入研究黄瓜DIR基因的生物学功能奠定重要基础。【方法】基于已报道的DIR基因HMM模型文件，利用HMMER软件包的hmmsearch程序从黄瓜蛋白数据库中筛选出可能的DIR基因ID，并利用在线工具Pfam和SMART进行验证，最终确定黄瓜DIR家族基因。利用ExPASy、TBtools、GSDS、MEME、MEGA、MCScanX和Circos等工具分析黄瓜DIR基因家族成员的理化特征、染色体定位、基因结构、系统进化树和共线性。基于黄瓜在不同组织和胁迫响应下的转录组测序大数据，利用黄瓜V3版本基因组信息进行转录组分析，检索黄瓜DIR基因在不同转录组测序分析中的表达情况，利用TBtools软件绘制表达热图，分析黄瓜DIR基因在不同组织和胁迫响应过程中的表达模式。【结果】在黄瓜中鉴定到23个DIR家族基因，分布于7条染色体上，编码氨基酸个数在78—684，分子量8.70—73.82 kD；系统进化分析将黄瓜DIR基因家族划分为3个亚族，每个亚族中的基因结构和motif基本一致；共线性分析发现黄瓜中有12个DIR基因与拟南芥中的19个DIR基因存在27种线性关系，黄瓜中有12个DIR基因与水稻中的11个DIR基因存在19种线性关系，而黄瓜中另外8个DIR基因比较保守，既不与拟南芥中的DIR基因存在共线性，也不与水稻中的DIR基因存在共线性；组织特异性表达分析发现有些黄瓜DIR基因在根、茎、花、果实、叶片等所有组织器官中的表达量均较低或不表达，有些黄瓜DIR基因在所有组织器官中的表达量均较高，而有些黄瓜DIR基因只在特定组织中表达，但在其他组织中不表达或低表达，表明不同的黄瓜DIR基因具有组织特异性表达模式；黄瓜DIR家族基因在胁迫响应过程中的表达模式分析发现CsaV3＿4G023490在黄瓜非生物胁迫和生物胁迫响应过程中均发生上调表达，表明该基因在黄瓜生长发育过程中发挥着重要作用。【结论】在黄瓜中共鉴定到23个DIR家族基因，分为3个亚族，每个亚族内的基因成员保守性高，不同亚族间的基因结构和蛋白保守结构域有所不同。黄瓜DIR基因在不同组织器官和胁迫响应下的表达模式具有差异性，协同调控了黄瓜的生长发育。

为明确芥蓝TCP家族相关基因在叶片发育中的功能，基于甘蓝基因组数据，分析鉴定出芥蓝TCP基因家族有40个成员，参考拟南芥同源TCP基因注释及拟南芥数据库基因的相对表达量，初步选取BoTCP21和BoTCP25基因，采用qRT-PCR分析。结果表明：BoTCP21和BoTCP25均在叶片中大量表达，但二者在真叶不同部位的表达模式存在差异。为了进一明确TCP基因家族中参与叶片发育的基因，采用生物信息学方法并结合qRT-PCR对BoTCP家族进行了全面分析。结果显示：40个BoTCP都属于非分泌亲水性蛋白，分别定位在8条染色体上；系统进化树构建和亲缘关系分析将芥蓝中的BoTCP成员归为3类，其中PCF亚家族有19个成员，CYC亚家族8个成员，CIN亚家族13个成员。qRT-PCR结果表明：BoTCP21和BoTCP25在真叶和薹叶中具有较高的表达量，BoTCP14在开花结荚的植株根部表达量较高，而BoTCP16则在抽薹植株的根部表达量最高，说明BoTCP家族成员在组织部位表达存在时空特异性且广泛参与了植株的形态建成和器官发育。

【目的】从全基因组水平上鉴定玉米大斑病菌(Setosphaeria turcica)Homeobox转录因子家族及其分布,阐明该家族的序列及进化特征,分析该家族基因在病菌不同生长发育时期的表达规律。【方法】利用生物信息学手段搜索玉米大斑病菌全基因组数据库,鉴定Homeobox转录因子家族;采用MEGA 5.0软件进行系统进化树分析;利用在线工具GSDS(gene structure display server)(http://gsds1.cbi.pku.edu.cn/index.php)绘制基因结构图;

【目的】对毛果杨Populus trichocarpa ZHD (PtrZHD)家族进行生物信息学以及干旱胁迫下表达特性分析，为研究PtrZHD在干旱胁迫中的功能提供参考。【方法】利用生物信息学方法从全基因组水平鉴定出毛果杨ZHD家族全部成员，并对其进化、理化性质、基因结构、保守基序、启动子顺式作用元件和表达特性进行分析。【结果】毛果杨ZHD家族包括21个成员，可分为7个亚家族；有8对同源基因，且非同义替换率(K_a)/同义替换率(K_s)值远小于1。该家族成员理化性质存在差异，但其结构较为保守，均含有Motif 1；启动子区含有数量不等的激素和非生物胁迫响应元件，不同基因之间响应元件的种类存在差异。在毛果杨PtrZHDs中，分别有1、7和13个基因在根、茎和叶组织中具有偏好性表达特征；PtrZHD家族成员对干旱胁迫的响应具有组织和时间表达特异性，在根、茎和叶部组织中各成员的表达量不同，但随着干旱胁迫时间的增加均呈先上升后下降的趋势。【结论】PtrZHD家族基因对干旱胁迫有不同程度的响应，可调控毛果杨对干旱胁迫的应答。图6表2参27

[目的 ]本研究通过系统检测沙棘UGT基因家族成员,探讨其结构特征及潜在功能,为解析沙棘类黄酮糖苷生物合成机制及其积累模式奠定基础。[方法 ]利用BLASTP和hmmsearch搜索沙棘基因组,并通过Pfam、CDD和SMART数据库验证保守结构域。使用Prot-Param、MUSCLE、MAGA7.0、MEME、MCScanX等工具分析蛋白理化性质、系统发育、蛋白基序和基因结构及基因复制事件。利用转录组数据和实时荧光定量PCR分析沙棘UGT基因在果实发育过程中的表达模式。[结果 ]本研究从沙棘基因组中鉴定出89个沙棘UGT基因,全部包含UGT基因家族保守结构域PSPG box。通过分析发现:沙棘UGT蛋白长度为266～533氨基酸,平均分子量50.00 KDa,平均等电点5.89。根据系统发育关系,可将89个沙棘UGT分为16个组,其中,A组包含最多的沙棘UGT基因家族成员。除7号染色体外,UGT基因共分布在11条沙棘染色体上。研究发现,串联重复是导致沙棘UGT基因家族扩张的主要复制事件。转录组分析和实时荧光定量PCR表明,大部分基因具有广泛的果实发育阶段特异性表达。[结论 ]本研究提供了沙棘UGT基因家族的完整信息,为进一步研究沙棘基因家族成员的生物学功能提供了数据参考和理论依据。

【目的】探明瓠瓜KNOX基因家族特征，在瓠瓜全基因组水平上对KNOX基因家族成员进行鉴定和分析，研究其组织表达模式。【方法】利用生物信息学方法，对瓠瓜KNOX家族基因成员进行鉴定、编码蛋白理化性质分析及家族成员结构域和保守基序分析，利用MEGA 5.05软件构建系统进化树，采用实时荧光定量PCR技术，分析KNOX基因在瓠瓜不同组织(根、茎、叶、花、果实)中的表达情况。【结果】鉴定到12个瓠瓜KNOX家族基因，这些基因分布在6条染色体上，且均定位于细胞核。大多数瓠瓜KNOX家族蛋白含有同源异型盒蛋白特有的4种结构域KNOX1、KNOX2、ELK和Homeobox＿KN。系统进化分析将瓠瓜KNOX基因分为KNOX ClassⅠ和KNOX ClassⅡ2大类群，2个类群又可进一步分为5个分支，其中ClassⅠ-C为瓠瓜特有的进化分支。KNOX基因启动子区有多个激素响应元件和植物生长相关元件。组织表达分析发现，瓠瓜KNOX基因具有不同的表达模式，KNOX ClassⅠ类基因表达具有组织特异性，在根、茎和果实中表达量较高，在其余组织中表达量较低或不表达，其中Lsi04G014450在根和茎中表达量最高，Lsi11G006380在果实中表达量最高；KNOX ClassⅡ类基因在所有组织中均有较高表达。【结论】瓠瓜KNOX基因家族有12个成员，分布于6条染色体上，在进化关系上可分为5组，具有瓠瓜特有的进化分支，在不同组织中的表达具有特异性。

为了探究大麻二酚酸合成酶基因(CBDAS)家族在大麻中的功能,利用生物信息学方法对大麻二酚酸合成酶基因(CBDAS)进行了全基因组鉴定,并系统分析其理化特性、进化发育、基因结构、启动子顺式结合元件及表达模式。结果表明,大麻CBDAS基因家族包含5个成员,仅分布在2,7号染色体上;理化性质分析显示,大麻CBDAS基因家族编码的氨基酸数目为541～545 aa,分子质量为介于61.49～62.41 ku,等电点为6.90～9.00;系统进化分析显示,来自植物、动物和微生物的33个CBDAS基因可分为3个家族,多数CsCBDAS基因家族成员属于同一亚家族;启动子区顺式作用元件预测表明,CsCBDAS基因家族成员均富含光响应元件;组织特异表达分析显示,CsCBDAS1和CsCBDAS2在雌蕊中表达最高,CsCBDAS4、CsCBDAS5在根中表达量最高;大麻CBDAS基因在高大麻二酚(CBD)和低CBD含量品种的表达模式不同,CsCBDAS1表达量在高大麻二酚(CBD)品种中高于低CBD含量品种;外源重金属Cd处理CsCBDAS4和CsCBDAS5基因表达量显著上调表达;CsCBDAS1、CsCBDAS2及CsCBDAS5表达量在黑暗和光照处理间表现出显著差异,光照处理下,CBDAS1基因表达量显著低于黑暗处理,而CsCBDAS2及CsCBDAS5的表达量高于黑暗处理。这些基因可能参与大麻的生长发育与大麻素的合成,该研究为后续大麻CBDAS基因家族的功能研究奠定了基础。

为了进一步理解己糖转运蛋白基因家族的结构特征以及其在植物逆境信号转导过程中的作用,基于小桐子基因组数据,鉴定到15个小桐子己糖转运蛋白基因,并对其基因结构、进化关系、密码子偏性及低温表达特性进行了分析。结果表明,15个小桐子HXT基因聚类为5个亚组,包含25个外显子,且在基因上游调控区包含ABA、乙烯、干旱及低温等逆境响应元件。编码蛋白质都富含疏水性氨基酸,具有典型的12个α-螺旋跨膜区。另外,部分小桐子HXT基因存在交叉共线性关系及基因倍增现象,且密码子的第3位碱基偏好使用U或A。荧光定量分析显示,Jc

【目的】探明瓠瓜KNOX基因家族特征,在瓠瓜全基因组水平上对KNOX基因家族成员进行鉴定和分析,研究其组织表达模式。【方法】利用生物信息学方法,对瓠瓜KNOX家族基因成员进行鉴定、编码蛋白理化性质分析及家族成员结构域和保守基序分析,利用MEGA 5.05软件构建系统进化树,采用实时荧光定量PCR技术,分析KNOX基因在瓠瓜不同组织(根、茎、叶、花、果实)中的表达情况。【结果】鉴定到12个瓠瓜KNOX家族基因,这些基因分布在6条染色体上,且均定位于细胞核。大多数瓠瓜KNOX家族蛋白含有同源异型盒蛋白特有的4种结构域KNOX1、KNOX2、ELK和Homeobox＿KN。系统进化分析将瓠瓜KNOX基因分为KNOX ClassⅠ和KNOX ClassⅡ2大类群,2个类群又可进一步分为5个分支,其中ClassⅠ-C为瓠瓜特有的进化分支。KNOX基因启动子区有多个激素响应元件和植物生长相关元件。组织表达分析发现,瓠瓜KNOX基因具有不同的表达模式,KNOX ClassⅠ类基因表达具有组织特异性,在根、茎和果实中表达量较高,在其余组织中表达量较低或不表达,其中Lsi04G014450在根和茎中表达量最高,Lsi11G006380在果实中表达量最高;KNOX ClassⅡ类基因在所有组织中均有较高表达。【结论】瓠瓜KNOX基因家族有12个成员,分布于6条染色体上,在进化关系上可分为5组,具有瓠瓜特有的进化分支,在不同组织中的表达具有特异性。

[目的 ]本研究通过系统检测沙棘UGT基因家族成员,探讨其结构特征及潜在功能,为解析沙棘类黄酮糖苷生物合成机制及其积累模式奠定基础。[方法 ]利用BLASTP和hmmsearch搜索沙棘基因组,并通过Pfam、CDD和SMART数据库验证保守结构域。使用Prot-Param、MUSCLE、MAGA7.0、MEME、MCScanX等工具分析蛋白理化性质、系统发育、蛋白基序和基因结构及基因复制事件。利用转录组数据和实时荧光定量PCR分析沙棘UGT基因在果实发育过程中的表达模式。[结果 ]本研究从沙棘基因组中鉴定出89个沙棘UGT基因,全部包含UGT基因家族保守结构域PSPG box。通过分析发现:沙棘UGT蛋白长度为266～533氨基酸,平均分子量50.00 KDa,平均等电点5.89。根据系统发育关系,可将89个沙棘UGT分为16个组,其中,A组包含最多的沙棘UGT基因家族成员。除7号染色体外,UGT基因共分布在11条沙棘染色体上。研究发现,串联重复是导致沙棘UGT基因家族扩张的主要复制事件。转录组分析和实时荧光定量PCR表明,大部分基因具有广泛的果实发育阶段特异性表达。[结论 ]本研究提供了沙棘UGT基因家族的完整信息,为进一步研究沙棘基因家族成员的生物学功能提供了数据参考和理论依据。

[目的]在全基因组水平对小麦(Triticum aestivum L.)OSCA家族成员(TaOSCAs)进行鉴定,以明确TaOSCA在小麦中的潜在生物学功能。[方法]利用BLASTP结合保守结构域筛选的方法鉴定小麦TaOSCA蛋白。用生物信息学方法分析TaOSCA基因家族成员系统发育关系,基因结构及跨膜结构域;利用公开的转录组数据分析TaOSCA在小麦不同发育阶段不同器官/组织的表达谱;以小麦品种中国春为材料,进行PEG-6000模拟的渗透胁迫处理,用实时荧光定量PCR检测叶片TaOSCAs在渗透胁迫下表达模式。[结果]从小麦中鉴定到35个TaOSCA蛋白;TaOSCAs分为TaOSCAⅠ、TaOSCAⅡ、TaOSCAⅢ和TaOSCAⅣ4个亚家族,分别有15,16,3和1个成员。基因结构分析结果表明,除TaOSCA4.1-4B外,其余TaOSCA均包含多个外显子。跨膜结构域分析结果表明,除TaOSCA2.3-4B有5个跨膜结构域(TMs)外,其余TaOSCA成员均包含8～11个TMs。TaOSCAs在小麦不同发育时期各器官中的表达结果表明,TaOSCA1.1-3B/3D、TaOSCA1.2-1A/1B/1D、TaOSCA1.3-1A、TaOSCA1.6-1A/1D和TaOSCA3.1-2B/2D在全生育期高水平表达,TaOSCA1.4-2A/2B/2D在不同发育阶段的根组织中均高表达,TaOSCA1.5-1B/1D在生殖生长期穗中高表达,TaOSCA1.6-1B、TaOSCA2.6-3B/3D在生殖生长期籽粒中高表达,表明这些基因可能在根、穗及籽粒发育过程中发挥重要功能。渗透胁迫处理下,TaOSCA2.4和TaOSCA4.1表达下调,TaOSCA2.5的表达无显著变化,TaOSCA1.1、TaOSCA1.2等11个TaOSCA基因表达上调,说明这11个基因可能在小麦苗期响应渗透胁迫中发挥作用。[结论]TaOSCA基因可能在小麦响应渗透胁迫的过程中发挥重要作用。

脱落酸(ABA)是一种内源性激素,在调节植物的生长发育、非生物胁迫和生物胁迫反应方面发挥着重要作用。Pyrabactin resistance 1-like(PYR/PYL)蛋白在ABA信号传导途径中起核心调控作用。然而,在蒺藜苜蓿(Medicagotruncatula)中并没有关于这个家族的细节。本研究根据拟南芥(Arabidopsis thaliana)中的14个PYL基因,利用BLAST方法鉴定蒺藜苜蓿基因组中的PYL家族成员。结果表明,鉴定得到14个蒺藜苜蓿PYL基因家族成员,且所有成员均含有PYR-PYL-RCAR-like功能域。预测启动子中的顺式元件发现所有蒺藜苜蓿PYL基因家族成员均含有与胁迫以及激素相关的元件。基因本体论(GO)富集分析与KEGG通路分析显示其主要参与到脱落酸信号通路及对外界刺激的反应。基因表达模式分析显示,其在逆境胁迫下表现出特定的表达模式。本研究为进一步研究蒺藜苜蓿PYL基因奠定了基础。