黄瓜DIR家族基因的全基因组鉴定及其表达分析

【目的】基于黄瓜基因组信息和转录组测序大数据,利用生物信息学手段,对黄瓜中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个亚族,...     

谷子抗锈病反应相关MYB转录因子的鉴定与表达

【目的】谷子锈病是影响其产量和品质的重要因素之一。鉴定谷子抗锈病相关的MYB转录因子,为谷子抗锈病机理研究奠定基础。【方法】 利用real-time PCR技术,检测9个SiMYBs转录因子基因在（1）谷子孕穗期根、茎、叶和穗4个部位的表达情况,（2）在谷子抗（resistance,R）、感（susceptibility,S）锈病反应120 h内的表达丰度差异,（3）在植株外接水杨酸（salicylic acid,SA）和茉莉酸甲酯（methyl jasmonate,MeJA）后24 h内的表达变化;比较SiMYBs基因在谷子R、S、SA与MeJA 4种反应中的表达模式;选择抗病相关SiMYBs基因进行转录激活活性检测和亚细胞定位。【结果】 SiMYB100在叶部表达量最高,其余8个基因均在根部表达量最高,SiMYB074最显著。5个基因的表达与抗病相关：SiMYB074和SiMYB202受锈菌侵染诱导表达,抗病反应早期的表达量明显高于感病反应;SiMYB041和SiMYB177在抗病反应早期下调表达,后期上升至接种前,而感病反应中保持低水平表达;SiMYB100在接种后的48 h内,抗病、感病反应表达模式相反。在响应外源激素SA和MeJA反应中,SiMYBs基因的表达量均发生不同程度的变化。4个基因（SiMYB074、SiMYB100、SiMYB174和SiMYB202）在R、SA与MeJA反应中的表达模式一致,不同于S反应。5个抗病相关SiMYBs基因具有转录激活活性,其编码蛋白质定位于细胞核中。【结论】 鉴定到SiMYB041、SiMYB074、SiMYB100、SiMYB177和SiMYB202 5个基因的表达与谷子抗锈病相关;SiMYB074与SiMYB100在谷子生长发育和抗病过程中都发挥一定的功能;SiMYB074、SiMYB100、SiMYB174和SiMYB202 4个基因可能通过SA和JA信号途径参与谷子的早期抗病反应。     

芥菜全基因组中Hsf基因家族鉴定与分析

【目的】通过鉴定与分析芥菜基因组中Hsf转录因子家族成员,为芥菜Hsf基因功能研究与遗传改良提高抗逆性提供理论依据。【方法】利用生物信息学方法分析芥菜Hsf热激转录因子家族成员的功能结构、保守基序、启动子顺式作用元件、系统进化、共线性及采用RNA-seq验证Hsf低温胁迫的基因表达。【结果】芥菜基因组中鉴定出71个Hsf基因家族成员,分布于18条染色体上,归类为3个亚家族,蛋白均含有DBD和HR-A/B结构域。BjuHsf启动子区域包含与逆境胁迫、激素、生长发育等相关顺式作用元件。系统进化与共线性分析表明,芥菜Hsf家族成员与大白菜Hsf家族成员具有更近的亲缘关系。在低温胁迫下,BjuHsf基因表达分析表明8个BjuHsf基因显著上调表达。【结论】这些显著差异表达BjuHsf基因与芥菜抗寒性极其相关,可作为芥菜耐寒性遗传改良的候选基因。     

苹果NLP（Nin-Like Protein）转录因子基因家族全基因组鉴定及表达模式分析

【目的】探究苹果NLP转录因子全基因组特征与表达模式,以便更深入地了解其结构特点与作用机制。【方法】基于本地BLAST数据库和Pfam数据库两大数据库,运用blastp和hmmsearch两种查询策略,对苹果全基因组范围内的NLP转录因子家族成员进行鉴定。经过严格筛选与确认后,对搜索结果进一步展开分析,主要分为3部分,包括NLP蛋白分析、NLP分析与苹果NLP表达分析,其中ProtParam、Clustal Omega、MEGA7、MEME 5.0.2、SOPMA、Phyre ²、WoLF PSORT、STRING等程序或软件用于蛋白质分析,基因分析使用MG2C、GSDS 2.0、PlantCARE、psRNATarget等在线工具,对于苹果NLP表达情况利用qRT-PCR进行定量检测。【结果】从苹果全部蛋白数据库中共筛选出6个NLP成员,系统进化分析可将它们分为I、II、III三类;蛋白二级结构以无规则卷曲为主,其次是α-螺旋,占比最小的是β-转角;亚细胞定位预测均定位于细胞核中,符合转录因子的特征;染色体定位显示5个基因（MDP0000584547除外）定位于4条染色体上;启动子分析发现大量与激素和逆境响应相关的顺式作用元件,暗示它们可能参与到激素和逆境信号的调控过程中,此外还识别到一个响应氮的GCN4作用元件,进一步说明该类转录因子与氮素有密不可分的关系;通过定量检测揭示苹果NLP家族在茎、叶等组织高表达的特征模式,并且表达分析结果也证实苹果NLP响应氮饥饿、干旱胁迫等。【结论】通过对苹果全基因组的分析,共识别6个NLP转录因子,对6个基因的结构与蛋白保守域分析,发现它们之间具有极高的相似性、保守性,同时又存在差异。类比拟南芥NLP蛋白的关联网络,推测与拟南芥NLP7同源性最高的MDP0000132856可能也具有复杂的功能。     

葡萄SAP家族的鉴定与表达分析

【目的】通过生物信息学分析明确SAP（Stress Associated Protein）家族在葡萄基因组中的数量、结构,利用qRT-PCR技术分析SAP家族的生物节律和在激素、非生物胁迫处理下的表达模式,为进一步探究SAP家族在葡萄中的作用奠定基础。【方法】在前期从山葡萄‘双优’（Vitis amurensis cv. Shuangyou）和欧洲葡萄‘红地球’（Vitis vinifera cv. Red Globe）0℃寒胁迫下转录组数据库中筛选出表达量明显上调的基因SAP15基础上,利用NCBI和葡萄基因组数据库中的BLAST功能,根据SAP家族的保守结构域,鉴定葡萄基因组中的SAP。利用生物信息学软件DNAMAN5.0、MEME、GSDS2.0、ExPASy和MEGA6等对VvSAPs序列、基因结构、蛋白质结构、理化性质、染色体定位和进化树等进行分析。采用qRT-PCR方法检测VvSAP家族的生物节律和在激素、逆境处理下的表达特征。【结果】从葡萄（Vitis vinifera）全基因组中鉴定得到15个SAP家族成员,均含有AN1保守结构域,大多含有A20保守结构域。按照其保守结构域和在染色体上的位置,依次命名为VvSAP1—VvSAP15,可分为Class Ⅰ、Class Ⅱ和Class Ⅲ三类。理化性质分析表明,VvSAP家族编码氨基酸数目介于109—293,理论等电点分布在7.99—9.68。染色体定位分析发现,该基因家族的15个基因分布于葡萄的9条染色体上,其中第8条染色体上分布最多,有3个VvSAP。亚细胞定位分析表明,VvSAP家族主要在细胞核、叶绿体和细胞骨架中进行表达。VvSAP家族的二级结构主要以无规则卷曲和α-螺旋为主。基因结构分析表明,VvSAP1—VvSAP12的DNA序列中均无内含子,VvSAP13—VvSAP15的DNA序列中有1个内含子,基因结构高度保守。qRT-PCR分析结果表明,在VvSAP家族成员中,VvSAP1和VvSAP9在无处理和各处理下（激素和非生物胁迫）均表达极低或不表达,初步鉴定为假基因。除VvSAP10在400 mmol·L ^-1 NaCl处理下呈下调表达,其余基因均在50 μmol·L ^-1 ABA、100 μmol·L ^-1SA和400 mmol·L ^-1NaCl处理下呈上调表达,其中,VvSAP10—VvSAP14显著响应50 μmol·L ^-1 ABA胁迫,处理24 h后相对表达量分别上调为0 h的37.19、36.63、21.69、58.34和267.35倍;VvSAP8和VvSAP11在NaCl处理4 h后相对表达量最高,分别为0 h的13.16和12.42倍;在4℃低温胁迫下,相对表达量上调最高的是VvSAP15,处理后8 h相对表达量为0 h的35.90倍。 【结论】从葡萄基因组中共鉴定出15个SAP家族成员,初步鉴定得出VvSAP1和VvSAP9为假基因,15个基因分别分布于9条染色体上,并且结构进化高度保守。所有成员均相响应逆境,且均有昼夜节律变化,但该基因家族成员在不同逆境胁迫下的响应程度和在不同时间点的表达模式存在一定的差异。     

茄子 TCP 转录因子的鉴定及胁迫处理下的表达分析

【目的】鉴定茄子 TCP 基因并对茄子 TCP 转录因子家族成员分析。确定其在胁迫、激素处理下的表达情况,为深入研究茄子 TCP 基因的功能奠定基础。【方法】从茄子基因组鉴定 SmTCP 基因家族成员,利用生物信息学方法对茄子 TCP 转录因子家族成员的理化性质、分布情况进行分析。以高温（42 ℃）、低温（4 ℃）、盐胁迫（200 μmol/L NaCl 溶液）、重金属（100 μmol/L CdCl2 溶液）、青枯菌、水杨酸及脱落酸处理茄子幼苗,利用 qRT-PCR 技术分析茄子 TCP 基因在 7 种处理及不同器官的表达情况。【结果】共鉴定到 32 个茄子 TCP 成员,并进行生物信息学分析。qRT-PCR 分析结果显示,茄子 SmTCP 基因受生物及环境胁迫响应上调,激素处理中,ABA 处理后茄子 TCP 整体出现上调,SA 处理后茄子 SmTCP 整体出现下调,茄子 TCP 在叶跟果中的的含量较高。【结论】茄子 TCP 基因受胁迫响应表达,SmTCP 成员响应胁迫表达情况及各器官中的表达量不同,Class Ⅰ类TCP 大多成员在热胁迫下存在明显上调响应。推测茄子 TCP 基因可能参与抵御逆境胁迫。     

紫花苜蓿MsWRKY42的分离、鉴定及其对非生物胁迫的响应

【目的】WRKY是转录因子大家族,在植物生长发育及逆境胁迫应答中发挥着核心调控作用。通过分析转录因子MsWRKY42在紫花苜蓿抗逆境过程中的作用,为进一步研究WRKY转录因子在紫花苜蓿抗逆分子调控中的作用奠定基础。【方法】通过同源比对方法从紫花苜蓿转录组基础数据库中获得MsWRKY42序列。使用MEGA-X对MsWRKY42蛋白序列及拟南芥序列进行多序列比对,构建系统发育树。通过PlantCARE预测分析MsWRKY42启动子的顺式作用元件。采用实时荧光定量PCR（qRT-PCR）方法分析MsWRKY42在紫花苜蓿不同组织中的表达量及其在NaCl（0.3 mol·L^-1）、PEG（15%）、4℃、40℃、低磷以及ABA（0.1×10^-3mol·L^-1）处理后的表达量变化。构建pCAMBIA1300-WRKY-GFP融合表达载体,通过农杆菌转化本氏烟草（Nicotiana benthamiana）,确定MsWRKY42蛋白的亚细胞定位。利用酵母单杂交系统检测MsWRKY42和启动子区域顺式作用元件W-box的体外结合活性。【结果】该基因包含1个1 692 bp的开放阅读框,编码563个氨基酸。多序列比对及系统进化树分析表明,该蛋白属于Ⅱb类WRKY家族成员,含有1个WRKY保守结构域和1个C₂H₂锌指基序,与拟南芥WRKY家族中的AtWRKY42相似度最高,故将其命名为MsWRKY42。在紫花苜蓿MsWRKY42启动子区域鉴定到多个顺式作用元件,主要包括胁迫响应、激素响应、生长发育等不同生命活动相关的调控作用元件。实时荧光定量PCR结果表明,MsWRKY42在紫花苜蓿各组织中均有不同程度的表达,其中根、叶中的表达量最高,茎、花、荚果中次之,芽中最低。经NaCl、PEG、低温、高温、低磷和ABA处理后的MsWRKY42表达量均有不同程度的上调。亚细胞定位结果显示,MsWRKY42蛋白定位在细胞核上。酵母单杂交系统检测结果显示,MsWRKY42能够与W-box顺式作用元件特异性结合。【结论】MsWRKY42为典型的WRKY转录因子,蛋白质定位在细胞核,能够与W-box顺式作用元件特异性结合;该基因在紫花苜蓿不同组织部位中均有表达,在根和叶中的表达量最高,且表达受NaCl、PEG、低温、高温胁迫和ABA激素的正向诱导。在紫花苜蓿中,该基因可能参与多种非生物胁迫反应。     

陆地棉BGLU基因家族成员的全基因组鉴定与表达分析

植物BGLU基因在生物和非生物胁迫防御中起着重要作用，利用生物信息学手段对陆地棉（Gossypium hirsutum） BGLU家族成员进行全基因组鉴定，对序列特征、保守域结构、染色体定位、启动子元件等家族特征进行了分析，并进一步结合转录组数据和qRT-PCR分析了BGLU基因家族的病菌胁迫响应模式。结果表明，共鉴定出53个陆地棉BGLU基因，根据系统发育进化关系分为6个亚族，部分基因在病菌胁迫下特异性表达。荧光定量和转录组测序结果表明，GhBGLU5、GhBGLU14、GhBGLU18、GhBGLU26、GhBGLU34在抗病棉花中的表达水平显著高于感病品种，推测这些基因正向调控棉花黄萎病抗性。上述结果为进一步分析棉花BGLU家族基因的功能以及分子机制提供一定的理论基础。     

低温胁迫下外源激素对番茄SlMYB102转录因子表达量的影响

探讨在低温胁迫下,6种外源激素对番茄转录因子SlMYB102基因表达量的影响。对长至四叶一心的番茄幼苗进行外源激素喷施和低温处理,激素处理采用水杨酸(SA)、生长素(IAA)、赤霉素(GA3)、6-苄氨基嘌呤(6-BA)、脱落酸(ABA)、茉莉酸甲酯(MeJA)进行叶面喷施,以喷施去离子水作为对照。分别在4℃低温处理后的0、1、3、6、9、12、24 h取样,利用实时荧光定量PCR法检测番茄转录因子SlMYB102基因的表达情况。结果表明,在ABA、GA3、IAA、SA、6-BA、MeJA的诱导下,SlMYB102基因对6种激素均有响应,并呈上调表达模式。SlMYB102转录因子作为MYB家族中R2R3-MYB亚家族的一员,可能参与这些激素对于逆境胁迫的调控网络。     

小麦IQM基因家族鉴定及非生物胁迫下表达分析

IQM基因是钙调素结合蛋白家族中的重要分支,在植物生长发育和应激反应中发挥重要作用。本研究利用生物信息学方法在小麦全基因组中鉴定出23个IQM基因家族成员,对其染色体位置、理化性质、系统进化关系、基因结构、蛋白保守结构域、启动子顺式作用元件和基因表达特性等进行了系统分析。结果表明,TaIQM基因家族成员随机分布在小麦18条染色体上,亚细胞定位结果显示所有基因均位于细胞核中,系统进化分析将其分为3个亚类,同一亚类间基因结构、蛋白保守结构域相似度较高,推测其功能相似;顺式作用元件分析表明,TaIQM基因家族成员启动子中含有多种与逆境胁迫及生长发育相关顺式作用元件;转录组数据分析显示,TaIQM基因家族成员在不同时期的根、茎、叶、穗和籽粒中表达量存在显著性差异;qRT-PCR分析显示,在小麦苗期地上部和地下部,TaIQM基因响应干旱、低温、高温、NaCl、ABA等多种胁迫,显示上调或下调表达。初步推断这些基因可能通过Ca²⁺信号通路参与非生物胁迫调控,研究结果为全面解析TaIQM基因结构与生物学功能、非生物胁迫响应分子机制提供依据。     

陆地棉转录因子基因GhMYB108的克隆及其在抗旱中的作用

【目的】MYB基因家族作为植物中最大的转录因子家族之一,在抵御逆境胁迫中发挥着重要的作用。克隆陆地棉MYB转录因子基因GhMYB108,并进行表达分析,验证其在干旱胁迫响应中的作用,为进一步研究GhMYB108调控陆地棉耐旱的分子机制奠定基础。【方法】根据干旱转录组数据分析,确定GhMYB108为干旱响应基因;运用聚合酶链式反应（PCR）从陆地棉根系cDNA中扩增目的基因;对GhMYB108进行基因结构特征、预测基因序列信息以及系统进化关系等生物信息学分析;利用Plant Care网站对获得的基因启动子序列进行分析;在不同逆境胁迫条件下,对GhMYB108的表达特性进行qRT-PCR分析;通过亚细胞定位确定GhMYB108蛋白在细胞中的位置;利用酵母试验验证其转录活性;使用病毒诱导的基因沉默技术（virus induced gene silencing,VIGS）沉默GhMYB108,并用qRT-PCR检测基因沉默效率。观察沉默株系在干旱处理前后的表型变化,并统计存活率,采用试剂盒测定相关生理生化指标;通过对棉花叶片喷施ABA与氟啶酮试验来分析GhMYB108与ABA的关系。【结果】从陆地棉中克隆了GhMYB108（Gh_A10G1563）,其全长879 bp,编码292个氨基酸,其蛋白质相对分子量为33.288 kD,等电点为6.037,多重序列比对和保守结构域分析,发现GhMYB108含有2个高度保守的MYB结合结构域,属于典型的R2R3型MYB转录因子。不同物种亲缘关系分析发现,GhMYB108与AtMYB108、AtMYB78和AtMYB2的同源性较高,属于同一亚族,且已有研究发现AtMYB108、AtMYB78和AtMYB2与干旱或ABA信号通路相关。GhMYB108定位于细胞核,且具有转录激活活性。在干旱和对照植株中,GhMYB108均在根中表达量最高,茎中表达量最低,并且受自然干旱、18% PEG 6000模拟干旱、盐胁迫和低温等非生物胁迫诱导表达。GhMYB108沉默之后,在自然干旱条件下,沉默植株出现临界表型,与对照相比,其萎蔫更严重,且存活率降低,一些生理生化指标也发生显著变化,如叶片失水率加快,丙二醛含量升高,叶片相对含水量和脯氨酸含量减少,过氧化氢酶（CAT）和过氧化物酶（POD）活性降低,且通过DAB与NBT染色发现植物体积累了更多过氧化氢（H₂O₂）和超氧阴离子（O₂^-）。通过对棉花叶片喷施激素ABA或氟啶酮发现GhMYB108可受ABA信号的正调控。【结论】GhMYB108正调控棉花抗旱性,且受ABA信号的正调控。     

玉米钙依赖蛋白激酶全基因组鉴定及抗旱表达分析

【目的】分析研究玉米钙依赖蛋白激酶全基因组鉴定及抗旱表达。【方法】使用NCBI Blast、DNAMAN和MotifScan等程序分析玉米钙依赖蛋白激酶CDPK的蛋白结构域和系统发育。采用Real-time RT-PCR方法分析CDPK基因在干旱胁迫处理玉米幼苗中的表达,评价CDPK基因对玉米干旱胁迫反应能力。【结果】鉴定出了在玉米中含有39个CDPK基因。将CDPK基因分为4个组。每个群体中的成员有共同的蛋白质基序和外显子及内含子结构,共同的进化起源。39个玉米CDPK基因根据它们的系统发育关系分组,并锚定在特定的玉米染色体上。多数CDPK基因在玉米干旱胁迫中的表达水平发生改变,CDPK基因参与对干旱胁迫的响应。【结论】从玉米基因组数据库中鉴定出了39个CDPK基因,大多数玉米CDPK基因在不同组织和不同发育阶段表现出不同的表达水平,CDPK基因在玉米发育中发挥不同的作用。多数CDPK基因在干旱胁迫下的表达量均发生变化。     

基于盐胁迫转录组信息的蚕豆F-box基因家族分析

【目的】通过生物信息学方法分析蚕豆F-box基因家族成员的分布、结构及进化,研究家族成员在不同处理时间条件下的表达模式及对盐胁迫的响应,为该类基因生物学功能和盐胁迫机制的研究提供参考。【方法】基于蚕豆盐胁迫转录组测序（RNA-seq）数据,利用NR、Swiss-prot和PFAM 3个数据库和NCBI网站,对蚕豆F-box基因进行筛选注释;利用Web Logo 3、Prot Comp 9.0、MEGA-X和MEME等软件进行保守结构域、亚细胞定位、系统进化树和Motif等生物信息学分析。基于盐胁迫转录组数据分析蚕豆（yz17134耐盐和yz17078不耐盐）F-box基因家族在盐胁迫下的差异表达模式,并采用实时荧光定量PCR技术（qRT-PCR）检测部分家族成员在16和24 h的表达情况。【结果】基于盐胁迫转录组测序数据,注释得到161个蚕豆F-box基因,均含有F-box保守结构域。根据C端结构域的不同,将其分成11个亚族：FBX、FBXFBA、FBXLRR、FBXPP2、FBXKelch、FBXTUB、FBXFBD、FBXDUF、FBXACTIN、FBXWD40和FBO。保守结构域分析表明,F-box保守基序中包含1个极度保守的色氨酸残基。比较分析蚕豆F-box家族和拟南芥F-box家族共同构建的进化树,发现同一C端结构域的基因大多聚集在一起。亚细胞定位预测结果显示,124个F-box基因定位于细胞外,37个定位于细胞核中。基因结构分析表明,蚕豆F-box家族基因的DNA序列中均无内含子,且均由UTR区和CDS区组成。基于盐胁迫转录组数据的F-box差异表达模式分析表明,蚕豆F-box基因在2个不同处理时间点上的表达各不相同,在盐处理16 h的表达较为明显。qRT-PCR分析结果表明,在F-box家族成员中,共存在5个差异基因。其中Vf056266.1、Vf062764.1和Vf024236.1在盐处理16 h的表达量均上调,Vf060904.1和Vf045761.1在盐处理16 h的表达量均下调。【结论】蚕豆F-box基因家族注释得到161个蚕豆F-box基因,分为11个亚族。其中5个重要的F-box基因在不同盐处理时间的表达量存在差异。     

玉米钙依赖蛋白激酶38(ZmCDPK38)的生物信息学及表达特性分析

[目的]研究ZmCDPK38与其它物种中同源基因之间的遗传进化关系,通过与亲缘关系较近且功能已知的同源基因比较,对ZmCDPK38的功能进行理论预判.确定ZmCDPK38基因在不同组织和不同胁迫条件下的表达特性.[方法]生物信息学分析使用在线分析工具及MEGAX软件;基因表达分析采用Real-time PCR.[结果]...     

干旱胁迫下两个玉米自交系雄穗花器官分化期的转录组分析

【目的】 运用现代分子生物学及高通量测序技术,在全基因组水平鉴定耐旱相关的基因,研究玉米雄穗花器官分化期对干旱胁迫的响应机理, 为抗旱育种将提供新的理论基础。【方法】 以干旱胁迫15 d和正常浇水的耐旱自交系“PHBA6”和干旱敏感自交系“吉63” 花器官分化期的雄穗进行Illumina HiSeq 2000高通量测序分析,将测序得到的Unigene与玉米参考基因组进行比对,利用 RPKM 来衡量样本间的基因表达丰度, 以|log₂foldchange|≥1,P<0.05,FDR≤0.001筛选出样本间的差异表达基因, 并通过 Gene Onto 1ogy (GO ) 数据库、 KEGG pathway 数据库对筛选出的差异表达基因的功能和参与的调控路径进行分析。【结果】 与正常水分处理相比,耐旱自交“PHBA6”和干旱敏感自交系“吉63”在干旱胁迫下分别有1 394个和1 300个显著差异表达基因,其中共有差异表达基因831个,包括上下调表达趋势是一致的822个基因和9个具有相反的表达模式基因。除这些基因外,各自分别还有563个和469个基因型特异的干旱响应基因,这些特异的干旱响应基因的G0 富集分析的条目并不完全相同,除共同富集在刺激响应,细胞壁改变外,耐旱自交“PHBA6”还富集在激素合成调控,机氮化合物代谢、蛋白结合上,特异的干旱响应基因的KEGG 显著性富集分析在2个自交系中也不完全相同,耐旱自交“PHBA6”还有部分差异基因富集在生物素代谢通路上,另外,筛选出耐旱自交“PHBA6”中在干旱胁迫下差异表达的3个脱水素基因和3个特有差异表达的bZIP、MYB和ERF转录因子。【结论】 胁迫相关基因的特异性响应可能是材料间耐旱性存在差异的主要原因,耐旱自交“PHBA6”中的3个脱水素基因和3个特有差异表达的bZIP、MYB和ERF转录因子基因在玉米耐旱育种中可能有重要的利用价值。     

红橘响应褐斑病菌侵染的转录组学分析

【目的】 明确红橘（Citrus reticulata Blanco, cv. Hongjv）接种褐斑病致病菌——链格孢菌橘致病型（Alternaria alternata tangerine pathotype）后基因种类和表达量在转录水平的变化规律,确定红橘响应该致病型侵染的关键基因。【方法】 采用链格孢菌橘致病型接种红橘离体叶片,28 h后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,以甜橙基因组为参考,以|log₂ fold change|≥1,q-value≤0.01为阈值选取感病和健康红橘叶片转录组的差异表达基因,应用GO数据库对差异表达基因（differentially expressed gene,DEG）进行功能分类,KEGG分析代谢途径,MapMan软件分析生物胁迫信号通路相关基因的表达变化。采用qRT-PCR方法对测序结果进行验证。【结果】 红橘接种链格孢菌橘致病型28 h后产生大量与胁迫相关的差异表达基因,获得上调差异基因5 173个,下调差异基因6 555个。GO功能分类显示差异基因主要与蛋白结合、膜、氧化还原过程等相关。通过KEGG富集和MapMan软件分析发现,红橘在受链格孢菌橘致病型胁迫的过程中基础代谢被严重破坏。乙烯、水杨酸和生长素等寄主防御相关的植物激素信号转导途径多个基因表达出现差异,其中乙烯起主导作用,乙烯受体ETR 3个成员被不同程度激活,下游激酶和乙烯响应因子均上调,生长素大部分关键信号基因、绝大部分生长素响应因子ARF和水杨酸合成途径的基因均下调表达。同时,黄酮醇、花青素、萜类化合物和生物碱合成相关基因受该菌诱导显著变化,萜类合成中大部分基因下调,而黄酮类合成相关上调基因数量和表达趋势均强于表达下调基因,有抗虫和抑菌作用的硫代葡萄糖苷基因呈上调趋势。进一步研究发现,大量参与抗逆过程的转录因子如WRKY、bZIP、ERF、MYB、NAC被诱导激活,其中大部分WRKY和bZIP转录因子受该菌正向调控,超过50%的ERF家族基因表达上调;在转录因子调控下,PTI及ETI响应基因如受体激酶、R蛋白、NBS抗病蛋白等大量表达,多个PR家族抗菌蛋白基因上调表达,22个抗氧化保护酶系统POD成员基因受到活性氧信号激发大量表达。以上结果表明链格孢菌橘致病型侵染对寄主内部生理状态产生显著影响。选取了19个与植物抗病相关基因进行qRT-PCR分析,其基因表达趋势与测序数据一致。【结论】 获得了红橘响应链格孢菌橘致病型侵染的差异表达基因及显著上调表达基因,其主要富集于代谢过程、应激反应及转录调控等条目中,这些基因的相互协同调控是红橘对该致病型产生防御反应的重要机制。     

独行菜抗冷相关基因的克隆及表达分析

【目的】独行菜在种子萌发和幼苗生长阶段可耐受一定程度的低温胁迫,研究独行菜的低温耐受机制。【方法】设计兼并引物,从独行菜冷诱导叶片的cDNA中克隆获得LaCBF3和LaCOR15a基因核心片段。利用半定量RT-PCR法分析两个基因在不同胁迫处理下的表达情况。【结果】获得442 bp LaCBF3及405 bp LaCOR15a核心片段,经比对与拟南芥CBF3、COR15a基因核心片段相似性分别为84.38%、81.8%。半定量RT-PCR结果显示在独行菜幼苗中,LaCOR15a基因不仅响应低温胁迫,还能迅速响应高盐、干旱及ABA处理,不同条件下表达各有差异。LaCBF3仅对低温胁迫响应迅速。【结论】LaCBF3、LaCOR15a在独行菜幼苗响应低温胁迫的分子机制中具有重要的调控作用,且LaCOR15a基因也应答干旱、高盐等非生物胁迫。