基于转录组的洋葱bHLH转录因子家族鉴定与生物信息学分析

bHLH转录因子家族是真核生物中重要的转录因子家族，在植物生长发育、次生代谢和逆境应答中发挥重要作用。本研究利用转录组数据对洋葱基因家族成员进行筛选和鉴定，并对家族成员进行理化性质、保守基序、亚细胞定位、系统发育和蛋白互作生物信息学分析，同时对bHLH基因家族在洋葱幼苗中的表达情况进行分析。结果表明，从洋葱转录组共鉴定到36个bHLH基因，预测编码蛋白质含有的氨基酸数量为127～610个，且均为亲水性蛋白，系统发育分析表明，洋葱bHLH基因家族分为12个亚族，同一亚族的大多数基因具有相似的保守基序，大多数bHLH基因在洋葱幼苗中表达，不同基因的表达量水平差异较大。本研究结果为进一步研究洋葱bHLH转录因子家族的生物学功能奠定了基础。     

玉米脱水素基因家族的鉴定与分析

脱水素(dehydrin,DHN)属于LEA蛋白第二家族成员,是一种植物中广泛存在的亲水性蛋白,在干旱、低温和高盐等非生物胁迫的环境中起到重要作用。通过生物信息学方法对玉米全基因组DHN家族成员进行了鉴定,并进一步对其系统发育、基因结构、染色体定位和基因复制以及表达模式进行了系统分析。结果显示,在玉米DHN基因家族中共有5个家族成员,多物种系统发育进化树分析、基因结构以及基序分析都表明DHN家族成员在进化上具有高度的保守性。基因复制和物种间微共线性分析表明,在5个玉米DHN基因中存在着1对片段复制基因(ZmDHN1-ZmDHN2),玉米、高粱和水稻3个物种间存在2对直系同源基因(ZmDHN2-Sb04g032250.1,ZmDHN2-Os02g44870.1)。通过转录组表达数据分析表明,玉米DHN家族基因在不同发育时期具有不同的组织表达模式;同时,诱导表达模式分析表明ZmDHN基因的表达受到盐和干旱胁迫的显著诱导。该研究结果将为进一步鉴定玉米DHN家族重要的基因成员并对其开展功能分析奠定基础。     

栀子AP2基因家族成员鉴定及其在盐胁迫下的表达模式

栀子耐干旱、耐瘠薄,适应性强。AP2/ERF(APETALA2/乙烯反应元件结合因子)基因家族成员在调控植物胁迫中起着关键作用,分析AP2转录因子家族成员对探究栀子耐盐机制具有重要意义。基于全基因组数据,鉴定栀子AP2基因家族成员并进行理化性质和系统发育分析,明确栀子在不同浓度盐胁迫下的表达模式。结果表明：在栀子中鉴定出113个AP2基因家族成员,所编码蛋白平均分子质量为26 323.08 u,等电点为4.51～10.2,脂肪指数为30.47～71.02,蛋白不稳定且亲水;蛋白二级结构以无规则卷曲(54.82%)为主,大部分定位在细胞核中,AP2和DREB亚家族成员定位在细胞质中的概率较大;AP2基因家族成员分属AP2(11个)、DREB(38个)、RAV(1个)和ERF(63个)亚家族;AP2基因家族成员进化的主要动力是纯化选择;大部分AP2基因家族成员的基因复制类型属于分散复制或片段重复;AP2基因家族成员启动区域中含有较多G-box、ABRE和CAT-box顺式作用元件;ERF亚家族成员Ⅹ类基因GjERF19、GjERF21和GjERF46积极响应盐胁迫,这3个基因将为栀子盐胁迫...     

玉米β-甘露聚糖酶基因家族的生物信息学分析

[目的]鉴定出玉米基因组中的β-甘露聚糖酶基因(MAN),并分析其进化机制和表达模式。[方法]利用拟南芥MAN蛋白序列在玉米基因组数据库以及NCBI非冗余蛋白质数据库中比对搜索,获取玉米MAN家族全部成员;通过检索植物基因组复制数据库查找玉米MAN基因的片段复制事件;利用GSDS 2.0服务器分析基因外显子/内含子结构;采用MEGA7软件极大似然法构建系统发育树;利用芯片数据分析基因在不同组织器官中的表达模式。[结果]玉米基因组中含有6个MAN基因,分布于5条染色体上。没有玉米MAN基因起源于片段复制或串联重复事件。系统发育分析将植物MAN蛋白分为3个亚族(I～Ⅲ),说明其序列已发生分化。玉米MAN基因结构较为保守,蛋白质均含有甘露聚糖酶的核心基序。玉米MAN基因具有不同的表达模式,部分基因在发育叶片中表达量较高,部分基因在减数分裂期雄穗中表达量较高,而其余基因在胚芽鞘、初生根、授粉前果穗、花丝和发育种子中表达量较高。[结论]该研究结果为进一步研究玉米MAN基因功能奠定了基础。     

西瓜OSCA基因家族全基因组鉴定及胁迫响应分析

【目的】对西瓜OSCA基因家族进行成员鉴定、系统发育分析以及胁迫响应分析,为揭示西瓜OSCA基因生物学功能和抗逆分子机制提供参考依据。【方法】运用生物信息学的方法从西瓜全基因组数据库中鉴定出西瓜OSCA基因家族成员,根据染色体定位信息进行命名。对基因家族成员染色体位置、基因结构、共线性、启动子、蛋白结构和系统发育等进行全面分析。利用转录组数据和qRT-PCR分析西瓜OSCA基因在胁迫条件下的表达情况。【结果】西瓜OSCA基因家族有10个成员,不均等地分布在6条染色体上。西瓜OSCA蛋白分子量变化范围为24.59~91.97kD,氨基酸数量为214~809aa,多数定位于质膜上。共线性分析结果表明,西瓜OSCA基因在进化过程中发生了片段复制事件。西瓜OSCA基因家族分为三大类群,同一类群成员的外显子一内含子的组成模式、蛋白保守基序排列、蛋白二级结构数量比例和蛋白三级结构模型均相似。西瓜、水稻、番茄和拟南芥的OSCA基因被分为6个亚族,每个亚族均有西瓜OSCA基因分布。西瓜OSCA基因启动子区域含有厌氧诱导、冷胁迫应答、光反应和干旱胁迫应答等多种非生物胁迫响应元件。在干旱、低温和盐胁迫下,西瓜OSCA基因家族各成员表达量均有不同程度变化,其中ClaOSCA3和ClaOSCA5在3种不同胁迫条件下表达量均有显著差异（P<0.05）。【结论】西瓜OSCA基因在进化过程中具有一定保守性,与拟南芥、水稻和番茄OSCA基因存在较近亲缘关系。西瓜OSCA基因家族内部存在功能分化,ClaOSCA3和ClaOSCA5可能是胁迫响应机制中的重要抗逆基因。     

大麦TIFY基因家族成员鉴定及表达分析

【目的】对大麦TIFY基因家族成员进行鉴定及表达分析,为进一步探究TIFY基因家族在大麦生长发育与胁迫响应中的作用机理打下基础。【方法】基于TIFY家族蛋白的保守域特征,利用HMMER从大麦中鉴定TIFY基因家族成员,利用采用生物信息学软件对其理化性质、保守基序、特征结构域、顺式作用元件、基因结构、系统进化及表达模式进行预测分析。【结果】从大麦中鉴定出15个HvTIFYs基因（HvTIFY1～HvTIFY15）,分布于5条染色体上,且大多数基因在染色体上成簇分布。15个HvTIFYs蛋白均具有TIFY家族蛋白的特征结构域（TIFY）,根据所含保守结构域的不同,可分为ZML（4个）和JAZ亚族（11个）,且亲水性蛋白（14个）和偏碱性蛋白（11个）居多,但均定位于细胞核;二级结构相似度较高,均由α-螺旋、β-转角和无规则卷曲组成,除HvTIFY7蛋白外,其余蛋白二级结构所占比排序:无规则卷曲>α-螺旋>β-转角。HvTIFYs基因结构存在明显差异,其中,JAZ亚族11个基因的内含子数为0~6; ZML亚族4个基因的内含子数为6～7个,系统发育进化树上相邻分支的基因具有较相似的基因结构。HvTIFYs基因启动子区域富含光、激素和胁迫等顺式作用元件,种类及分布均呈多样性。5个物种的79条TIFY蛋白分为4个组,恰好与TIFY家族的4个亚族对应,其中,ZML、TIFY和JAZ亚族包含单、双子叶植物的TIFY蛋白,而PPD亚族仅含有双子叶植物的TIFY蛋白。15个HvTIFYs基因在不同组织器官中的表达量存在明显差异,其中HvTIFY1、HvTIFY2和HvTIFY8基因在8个组织中的表达量均较高,HvTIFY10和HvTIFY15基因表达量中等,HvTIFY6基因表达量较低; HvTIFY11基因不表达。15个基因在根的不同组织中对盐胁迫的敏感程度不同。【结论】从大麦中鉴定出的15个HvTIFYs基因存在一定的功能分化,具有明显的组织和时空特异性,推测其在大麦逆境响应和激素调节中具有重要调控作用。     

水稻CIPK基因家族的鉴定及OsCIPK5受稻瘟病菌诱导的qRT-PCR分析

【目的】植物通过启动一系列信号传导过程来应对外部环境,这些过程通常涉及多种蛋白激酶,包括钙调神经磷酸酶B样蛋白互作激酶(calcineurin B-like protein-interacting protein kinases, CIPKs)。为更加清晰全面地了解水稻CIPK基因家族,本研究根据最新的基因组测序数据对水稻基因组中的CIPKs进行了鉴定。【方法】通过探讨拟南芥和水稻中CIPKs蛋白家族的结构特点,结合生物信息学和qRT-PCR技术系统分析水稻中CIPKs家族蛋白的结构。结合转录组数据,比较了粳稻云引受稻瘟病菌诱导后的表达情况。【结果】根据最新的水稻基因组数据,鉴定出31个水稻OsCIPK基因。系统发育树分析结果表明,31个OsCIPK基因可分为5个亚家族,这些亚家族具有不同的外显子-内含子和UTR的结构特点。从广谱抗稻瘟病品种粳稻云引受稻瘟病菌诱导的基因表达谱的趋势聚类中筛选出了OsCIPK5基因并对其进行了表达分析,结果表明,云引中OsCIPK5基因受稻瘟病菌的诱导表达。【结论】内含子缺失和片段重复在水稻OsCIPK基因家族的扩展中起到重要作用,同时OsCIPK5受到稻瘟病菌的诱导表达。     

玉米COBRA家族成员全基因组鉴定与表达模式分析

COBRA基因编码糖基磷脂酰肌醇锚定的植物特异性蛋白，在初生和次生细胞壁的纤维素生物合成中发挥重要作用。通过生物信息学方法，从玉米基因组中发现了9个COBRA家族基因，并对其基因结构、系统发育和表达模式等进行了分析。结果表明，9个COBRA家族基因都含有CCVS保守结构域，并且均定位于细胞膜上。系统进化分析结果显示，该家族可以分为2个亚族，每个亚族内的基因具有相似的基因结构和理化性质。基因表达分析结果表明，所有COBRA家族成员响应多种非生物胁迫，且当植株受到紫外照射处理时，ZmCOBL1～ZmCOBL3和ZmCOBL7～ZmCOBL9这6个基因均表现出上调应答。     

海南冬季诱导火龙果开花的补光条件

以红肉火龙果(Hylocereus undulatus Britt)‘金都1号’为材料,通过不同的光谱及光照时间处理,研究不同的补光条件对诱导海南冬季火龙果开花的影响。结果表明:1)组合黄光和红光能够诱导火龙果成花,蓝光不能诱导火龙果成花。2)补光时间的长短和光谱的交互作用对火龙果成花有显著影响。3)不同月份,诱导火龙果成花所需的补光时间不同。采用15W组合黄光,在2月初至2月下旬,需补光5 h诱导成花;在2月下旬至3月上旬,补光4 h可诱导成花;在3月上旬到3月底,补光3 h就可诱导火龙果成花。     

陆地棉TCP家族基因鉴定及组织表达分析

【目的】鉴定、分析陆地棉TCP家族基因,为进一步研究TCP基因在棉花植株内的生物学功能奠定基础。【方法】基于2019年最新组装的陆地棉TM-1参考基因组,通过Pfam网站获取TCP家族HMM模型文件,使用HMMER网站鉴定陆地棉TCP基因,CottonFGD网站获取陆地棉TCP基因组蛋白序列。利用MapInspect进行染色体定位、MEGA 7.0进行多序列比对聚类分析、MEME网站motif预测、TBtools软件鉴定基因结构以及陆地棉TCP基因组织特异性表达分析。【结果】共鉴定出63个陆地棉TCP基因,其中39个Class I亚族,24个ClassⅡ亚族;ClassⅡ亚族中包括17个CIN亚类,7个CYC/TB1亚类。染色体定位发现该63个TCP基因在陆地棉22条染色体上均有分布,其中A组染色体上分布有33个基因、D组上分布有30个基因。所有TCP蛋白均包含TCP结构域。TCP基因外显子、内含子结构及长度在同一亚家族内具有相似性。TCP家族基因中有12个基因在陆地棉纤维中优势表达,32个基因在花器官中优势表达,16个基因在营养器官：根、茎和叶中优势表达。【结论】获得了与陆地棉生长发...     

西洋参bHLH转录因子家族生物信息学分析

bHLH是真核生物中一类编码转录调控因子的基因家族,在植物生长发育中具有重要的调控作用。通过对西洋参根、叶、花蕾样品进行转录组测序,筛选其中的bHLH家族基因,利用生物信息学方法对筛选到的基因进行表达模式、蛋白结构、细胞定位及系统进化分析。结果显示:从西洋参转录组中共鉴定得到62个西洋参bHLH家族基因,这些bHLH在根、叶、花蕾中具有4种表达模式。西洋参bHLH亚细胞定位预测主要在细胞核内,所有表达产物均包含HLH结构域。进化分析显示,西洋参中62个bHLH基因可分为18个亚类,其中第Ⅺ亚族中包含最多的bHLH成员,共包含11个bHLH蛋白;第Ⅲf、Ⅳa、Ⅲ(a+c)、Ⅰb(2)、Orphans、Ⅳb亚族中的成员数量均为1个,其他亚族为2～9个。     

闽楠bZIP基因家族鉴定和脱落酸处理下的表达分析

【目的】对闽楠Phoebe bournei bZIP (PbbZIP)转录因子家族成员鉴定,分析其对脱落酸(ABA)信号的响应水平。【方法】基于生物信息学方法,对PbbZIP基因家族进行了全基因组鉴定,分析其蛋白理化特性、基因结构、进化关系、启动子顺式作用元件和ABA处理下的表达分析。【结果】从闽楠12条染色体共鉴定出63个PbbZIPs基因,分为12亚族,不同亚族基因结构和基序差异显著,但同亚族高度保守。PbbZIP基因多数定位在细胞核,其编码蛋白长度为110～835个氨基酸,等电点为4.48～11.95,疏水性为-1.19～-0.19。分布于12条染色体的27对PbbZIPs基因存在共线性关系,是PbbZIP基因家族扩张的主要模式。PbbZIP基因上游启动子区域存在多种与非生物胁迫相关的作用元件,其中ABA、水杨酸、茉莉酸甲酯的响应元件较多。基因表达分析结果表明：2 mmol·L^-1ABA处理闽楠1～72 h,17个PbbZIPs基因在叶和根中被ABA信号不同程度地诱导表达,普遍上调,且根基因表达水平普遍低于叶片。【结论】63个PbbZIPs基因序列高度保守,...     

亚洲棉Argonaute蛋白家族的生物信息学预测与分析

【目的】在全基因组水平上鉴定亚洲棉(Gossypium arboretum L.)AGO蛋白家族成员,分析其进化关系,为亚洲棉遗传改良提供基因资源。【方法】用AGO蛋白保守域在亚洲棉基因组数据库中进行同源比对,并对AGO蛋白家族的进化关系、基因结构、结构域、染色体定位、三级结构及表达模式等进行生物信息学分析。【结果】从亚洲棉蛋白数据库中共鉴定出11个AGO蛋白家族成员,其编码区长度为2 949~9 807bp,定位在6条染色体上,其中第9号染色体上分布有4个基因;11个AGO家族成员分为5个亚族,亚族内各成员间具有很高的保守性;AGO家族成员均具有典型的α-螺旋跨膜结构,属于碱性蛋白;从表达模式看,AGO家族成员在纤维发育不同时期和不同组织中的表达量存在较大差异,其中5个基因在棉纤维次生壁后期高表达,在其他时期表达水平较低,存在组织特异性。【结论】成功鉴定了亚洲棉AGO蛋白家族,明确了其基本结构的进化关系及基因表达特点。     

割手密GST基因家族的分析及其对阿特拉津的响应

从全基因组的角度分析了割手密谷胱苷肽-S-转移酶(GST)基因家族成员的特征,共鉴定出129个GST。通过对割手密GST蛋白的系统发育分析,发现这些GST可以划分为4个亚家族,且每个亚家族内的蛋白排列紧凑,说明同亚族内部基因进化关系均比较亲近。此外,对在阿特拉津处理下GST基因家族中10个F亚族基因的表达水平进行检测,结果表明,根中的SpGSTF1、SpGSTF9、SpGSTF14、SpGSTF16、SpGSTF18、SpGSTF31在阿特拉津胁迫下表达显著上调,推测这6个基因可能在响应阿特拉津胁迫的过程中起重要作用。研究结果有助于揭示割手密GST家族的进化和功能演替,并对进一步了解割手密GST对阿特拉津的解毒响应机制提供理论参考,同时也为培养抗除草剂甘蔗品种提供基因资源。     

拟南芥HMGs家族成员全基因组分析

高迁移率族蛋白是一种非组蛋白染色体蛋白质.文中通过生物信息学的方法鉴定到17个HMGs家族基因,利用在线工具分析了HMGs家族基因的组成以及染色体定位信息,利用Clustal-W以及MEGA4.0软件构建了拟南芥HMGs家族的系统发育树,利用通过查询基因芯片结果分析了该家族基因的不同组织部位的表达情况,并推测了部分基因的功能,为进一步研究HMGs家族基因的功能提供了参考.     

红肉火龙果TST基因家族的生物信息学分析

为分离参与红肉火龙果果实可溶性糖积累的功能基因,基于转录组测序数据克隆液泡糖转运蛋白（Tonoplast sugar transporter,TST）基因家族,进行序列同源性、跨膜结构、亚细胞定位和系统发育等生物信息学分析。结果表明：红肉火龙果含有4个TST基因,开放读码框（open reading frame,ORF）长度为2 127～2 235 bp,编码长度为708～744 个氨基酸的序列。预测等电点为4.84～5.36,相对分子量为75.7～80.4 kDa;红肉火龙果TST家族氨基酸序列与甜菜和藜麦TST家族具有74.9%～83.4%的相似性;红肉火龙果TST成员均含11～12个跨膜结构域,定位于质体、叶绿体、液泡或内质网等细胞器。系统进化分析表明HpTST1属于TST1类,HpTST2.1和HpTST2.2属于TST2类,HpTST3属于TST3类。     

红肉火龙果TST基因家族的生物信息学分析

为分离参与红肉火龙果果实可溶性糖积累的功能基因,基于转录组测序数据克隆液泡糖转运蛋白(Tonoplast sugar transporter,TST)基因家族,进行序列同源性、跨膜结构、亚细胞定位和系统发育等生物信息学分析。结果表明:红肉火龙果含有4个TST基因,开放读码框(open reading frame,ORF)长度为2 127～2 235bp,编码长度为708～744个氨基酸的序列。预测等电点为4.84～5.36,相对分子量为75.7～80.4kDa;红肉火龙果TST家族氨基酸序列与甜菜和藜麦TST家族具有74.9%～83.4%的相似性;红肉火龙果TST成员均含11～12个跨膜结构域,定位于质体、叶绿体、液泡或内质网等细胞器。系统进化分析表明HpTST1属于TST1类,HpTST2.1和HpTST2.2属于TST2类,HpTST3属于TST3类。     

黄瓜DnaJ基因家族鉴定及对高温胁迫的表达响应

[目的]本文旨在研究黄瓜热激蛋白DnaJ(CsDnaJ)的家族基因鉴定,并分析其在高温胁迫环境下的响应模式。[方法]利用生物信息学软件及黄瓜的全基因组测序数据库、公共数据库对CsDnaJ进行基因家族成员、染色体位置、蛋白特性、结构特征、进化及复制模式和组织表达特性分析,同时通过实时定量PCR技术分析其在高温胁迫环境下的响应模式。[结果]黄瓜CsDnaJ基因家族共存在81个成员基因,其中Chr 3含DnaJ基因最多,为16个,Chr 4、Chr 7最少,均为9个;CsDnaJ基因所编码的蛋白序列长度为113～2 551 aa,且至少含有1个完整的J结构域核心区序列;亚细胞位置预测显示大部分蛋白定位于细胞核及叶绿体中。基因结构显示CsDnaJ基因均含有1～22个外显子。通过进化分析发现CsDnaJ蛋白可明显区分为9个亚簇,每个亚簇均存在拟南芥同源蛋白。共线性分析发现,CsDnaJ基因家族共存在17个共线基因,组成10对复制基因对,其中存在5对片段重复事件,3对串联重复事件,2对随机重复事件。组织表达特异性结果显示CsDnaJ基因在叶、茎、子房中的表达量略高于其在雄花、雌花、根及卷须中的表达量。荧光定量PCR结果显示,高温环境(昼/夜温度为42℃/35℃)可以诱导CsDnaJ基因的高表达,从而参与植物响应高温环境胁迫进程。[结论]黄瓜CsDnaJ基因家族核心结构域高度保守,在黄瓜响应高温胁迫中发挥重要作用。     

干旱胁迫下大苞萱草bHLH转录因子家族鉴定与分析

【目的】为探究大苞萱草bHLH基因家族成员特性，基于对干旱胁迫下大苞萱草叶和根转录组测序结果，鉴定并分析大苞萱草bHLH基因家族成员。【方法】利用生物信息学方法对大苞萱草bHLH转录因子家族基因进行系统发育、理化性质、二级结构、保守结构域及基因表达等分析。【结果】共鉴定出55个大苞萱草bHLH家族基因并分为11个亚族；bHLH蛋白理化性质差异较大，总平均亲水性为负值，均为亲水性蛋白；亚细胞定位预测大苞萱草bHLH蛋白主要分布在细胞核中；基因表达分析表明，叶中有26个上调基因，26个下调基因，根中有32个上调基因，22个下调基因，差异基因HmbHLH50的表达量变化显著，可能与大苞萱草的抗旱能力相关。【结论】本研究为挖掘bHLH转录因子家族基因的功能奠定基础，也为深入解析大苞萱草的抗旱机制提供理论依据。     

芝麻AQP家族的全基因组序列鉴定及其特征分析

【目的】从芝麻全基因组中分离鉴定水通道蛋白AQP（aquaporin）家族基因,并进行系统进化关系、连锁群定位、基因结构、跨膜结构域和亚细胞定位预测、保守性氨基酸残基以及青枯雷尔氏菌诱导表达分析,为芝麻AQP的功能研究与利用奠定基础。【方法】通过生物信息学手段,结合芝麻基因组注释信息,鉴定芝麻AQP家族成员序列信息,并用InterPro逐一进行验证。利用ClustalW2对芝麻、拟南芥和水稻的AQP以及马铃薯的XIPs进行多序列比对,用MEGA6.0构建进化树。通过MapInspect和Gene Structure Display Server 2.0进行连锁群定位和基因结构分析。采用ProtParam、WoLF PSORT和TMHMM Server v.2.0在线工具预测芝麻水通道蛋白的分子质量和等电点、亚细胞定位及跨膜结构域。通过芝麻、拟南芥和水稻AQP及马铃薯XIPs的蛋白多序列比对结果推测NPA基序、ar/R滤器及P1-P5的氨基酸残基。利用前期研究获得的转录组测序结果进行青枯雷尔氏菌诱导表达分析,并通过qRT-PCR技术对差异表达较为明显的12个芝麻AQP进行验证。【结果】系统分析鉴定了36个芝麻AQP家族基因,根据多序列比对及系统进化分析将其分为5个亚家族：13个质膜内在蛋白（PIP）、12个液胞膜内在蛋白（TIP）、8个类NOD26膜内在蛋白（NIP）、2个膜内在小分子碱性蛋白（SIP）和1个未知内在蛋白（XIP）,其中有34个基因定位在12个连锁群上。同一亚家族成员在基因结构、蛋白序列、亚细胞定位预测及保守性氨基酸残基等方面都较为相似。青枯雷尔氏菌诱导表达分析显示,NIPs、SIPs和XIPs表达无明显变化,部分PIPs和TIPs能够响应青枯菌诱导。SiPIP1;2、SiPIP1;3、SiPIP2;3和SiPIP2;4受青枯菌诱导后表达上调,SiPIP1;3和SiPIP2;3为持续上调,而SiPIP1;2和SiPIP2;4的表达先下调后上调;与之相反,表达下调较为显著的有SiPIP1;4、SiPIP2;1、SiPIP2;6、SiTIP1;1、SiTIP1;3、SiTIP2;1及SiTIP2;2。上述差异表达基因的qRT-PCR验证结果与转录组测序结果一致。【结论】通过全基因组分析,在芝麻中鉴定出36个AQP家族基因,分为5个亚家族,分布于12个连锁群上,大部分基因具有1-4个内含子（除SiNIP1;2有7个内含子外）。根据ar/R滤器和P1-P5氨基酸残基组成类型,预测不同亚家族AQP可能识别的底物。青枯菌诱导表达分析表明,部分PIPs和TIPs成员的表达发生了显著变化。