解 盼 刘 蔚 康 郁 华 玮 钱论文,2,3 官春云,2,3,* 何 昕,2,3,*

1 湖南农业大学南方粮油作物协同创新中心, 湖南长沙 410128; 2湖南农业大学油料作物研究所, 湖南长沙 410128; 3国家油料作物改良中心湖南分中心, 湖南长沙 410128; 4 湖南农业大学风景园林与艺术设计学院, 湖南长沙 410128

两熟制、三熟制是我国南方地区推广的主要种植模式[1], 为适应南方地区“稻–稻–油”轮作生产模式的变化趋势, 需要种植早熟、特早熟油菜品种。但大量的研究和生产实践表明, 油菜早熟性与抗逆性和丰产性之间始终存在着较大的矛盾, 早熟与低产、抗性差之间存在着连锁关系[2]。早熟品种生育进程快, 现蕾、抽薹、开花都较早, 抗冻性差并且花期耐春寒能力较弱, 一旦遭遇冰冻或持续的低温阴雨天气, 产量损失严重[3]。如2008年1月, 我国遭受了严重的低温冰雪灾害, 对冬油菜生产带来了严重的影响, 受灾面积约占全国冬油菜面积77.8%, 产量损失10.88%。在受灾区随机抽样调查发现, 受冻害影响较严重的品种主要为未经当地审定的早熟品种,而经过正规审定的油菜品种的受灾程度较轻[4]。此外, 2010年春季“倒春寒”引起的持续低温阴雨天气给我国长江中下游油菜生产带来严重影响, 平均减产10%~20%, 部分地区减产高达30%[5]。因此, 需要我们加强对耐寒、抗冻性强的早熟油菜品种选育研究。

低温作为世界最为严重的自然灾害之一, 影响植物的生长发育以及地理分布, 严重时甚至会造成植物死亡。目前, 低温信号途径研究比较清楚的是ICE1(inducer ofCBFexpression 1)–CBF(C-repeat(CRT)-binding factors)–COR(cold responsive)信号途径。研究表明,ICE1通过特异性结合CBF3启动子上的MYC结合位点(CANNTG), 正调控CBF3基因的表达[6]。CBF转录因子中的AP2 domain可直接与冷调节COR基因启动子中的CRT/DRE顺式元件结合, 进而诱导COR基因的表达。在拟南芥中,CBF1/2/3串联排列在第4条染色体上, 受冷诱导表达[7-10]。而CBF4不受冷诱导, 但在拟南芥中过表达CBF4可增强植物的抗冻性和耐旱性[11]。在拟南芥植株上过表达CBF1/2/3均明显增强其抗冻性[8-9,12]。与野生型相比,敲除CBF1和CBF3基因拟南芥植株抗冻性降低约60%[13]。CBF1/2/3三突变体在非冷驯化时没有或具有轻微冻敏感表型, 而在冷驯化后表现出强烈的冻敏感表型[14-16]。通过对CBF1/2/3三突变体进行RNA-seq分析发现,CBFs突变影响了全转录组10%~20%的COR基因表达[15]。大量研究证明, CBF基因在植物抗寒方面是一个保守的信号路径。抗冷相关基因(CBF4、COR6、COR15和COR25)在抗冻性强的甘蓝型油菜1801中的表达量显著高于在抗冻性弱的C20中的[17]。Maryam等[18]在葡萄中发现4个受冷胁迫诱导表达的CBF基因, 这些CBF在冷处理后几分钟开始表达, 其中CBF2持续表达几个小时,CBF4持续表达到第10天。

研究人员开始利用CBF基因改良作物的抗寒、抗冻性。在冷胁迫下, 超表达拟南芥CBF3基因可以增加木薯、烟草、水稻、油菜和地被石竹中脯氨酸的积累, 减少丙二醛的产生以及电解质外渗, 从而增强它们的抗寒性[19-22]。此外,OsCBF基因转入水稻后, 脯氨酸及可溶性糖的含量明显提高, 使得水稻的抗寒性、抗旱性和抗盐性都有所提高[23]。但是有研究发现, 超表达甘蓝型油菜BnCBF5和BnCBF17可以提高植株抗冻性, 但会严重抑制植株的生长发育,并推迟开花时间[24]。

为充分挖掘CBF基因, 利用CBF基因解决甘蓝型油菜早熟品种与抗冻性之间的矛盾, 本研究对拟南芥、甘蓝、白菜以及甘蓝型油菜的CBF家族进行系统地比较分析, 鉴定甘蓝型油菜、白菜和甘蓝CBF基因, 并对这些CBF成员的基因结构、分子特性、基因染色体分布、系统进化树、组织表达、不同逆境、不同激素以及不同低温胁迫处理下的表达模式进行比较分析, 筛选候选抗冻基因, 以期为进一步了解甘蓝型油菜CBF家族基因的生物学功能及其对低温尤其是冻害响应奠定基础, 同时为CBF家族基因在其他物种中的生物信息学研究提供参考。

1 材料与方法

1.1 试验材料

将甘蓝型油菜半冬性品种中双11号(Zhongshuang 11, ZS11)的种子在滤纸上发芽, 移栽到装有土壤或蛭石的花盆中, 在生长室内生长6周。生长室中条件设置如下: 22℃ 14 h光照/20℃10 h黑暗, 相对湿度60%。在8:00使用100 µmol L-1脱落酸(abscisic acid, ABA)、100 μmol L-1茉莉酸甲酯(methyl jasmonate, MeJA)、1 mmol L-1水杨酸(salicylic acid, SA)和10 μg mL-1乙烯利(ethephon, ETH)溶液喷洒叶片进行激素处理, 并在处理后的3、6和12 h时间点采集叶片样品。高温胁迫和低温胁迫, 在8:00—20:00期间, 生长室分别设置为40℃和4℃, 相对湿度60%。在胁迫处理的3、h和12 h时间点采集叶片样品。盐胁迫和PEG胁迫处理, 分别在8:00—20:00期间用200 mmol L-1NaCl或20%聚乙二醇6000(PEG-6000)溶液灌溉幼苗, 并在处理期间的3、6和12 h时间点采集根系样本。对照幼苗用去离子水处理, 每个处理重复3次, 每个处理用3~4株幼苗。采集的叶和根样本立即在液氮中冷冻, 然后在-80℃下保存以提取RNA[24]。

在光照培养箱中培育2个早熟的半冬性油菜品种(HX17和HX58)的幼苗, 生长室条件设置如下: 20℃14 h光照(6:00—20:00)/16℃10 h黑暗(20:00—6:00),相对湿度60%, 培养4周; 冷驯化冷害与冻害处理:在4℃条件下处理14 d后, 分别放入4℃和-4℃处理12 h; 非冷驯化冷害与冻害处理: 20℃ 14 h光照(6:00—20:00)/16℃10 h黑暗(20:00—6:00)处理6周后, 分别放入4℃和-4℃处理12 h; 其中, 4℃冷害/-4℃冻害处理12 h的时段为20:00—8:00 (黑暗10 h/光照2 h), 在处理后的上午8:00收集从顶部往下数的第3片叶子并立即保存在-80℃, 直到提取RNA[25]。

1.2 白菜、甘蓝及甘蓝型油菜CBF基因的筛选鉴定

从Ensembl Plants (http://plants.ensembl.org/index.html)数据库中获取拟南芥、白菜(Brassica rapa)、甘蓝(Brassica oleracea)、甘蓝型油菜(Brassica napus)的CBF基因全长DNA序列、CDS序列、蛋白的氨基酸序列及相关信息。首先在Ensembl Plants中下载拟南芥CBF基因的核酸序列和蛋白的氨基酸序列, 利用BLASTP比对甘蓝型油菜、白菜和甘蓝基因组, 以E<10–5为阈值分别筛选出对应的CBF候选同源基因,并利用BLASTP比对NCBI非冗余数据库鉴定候选基因是否为CBF同源基因。

1.3 生物信息学分析

利用DNAMan6.0软件对拟南芥、白菜、甘蓝和甘蓝型油菜CBF蛋白序列比对分析, 并使用MEGA7中的Neighbour-Joining (NJ)系统发育方法构建系统发育树[26]。利用MEME (http://MEME.nbcr.net/ MEME/cgi-bin/MEME.cgi)对拟南芥、白菜、甘蓝和甘蓝型油菜CBF进行保守结构域分析。使用Compute-pI/Mw工具(http://web.expasy.org/ComputepI/)和ProtComp 9.0 (http://linux1.softberry.com/)预测分子量(molecular weight, MW)、等电点(Isoelectric point prediction, pI)和亚细胞定位。使用TBtools软件(TBtools_v.6695)展示基因结构、具体的染色体分布情况。参考以往的研究报道[27], 在100 kb基因组范围内排列有2个及以上基因, 我们将其定义为串联重复基因。

1.4 表达热图分析

参考He等[24-25]以及本课题组尚未公开的数据进行热图分析, 采用TBtools软件(TBtools.v.6695)作图。

1.5 BnaCBF基因表达情况检测

参照植物总RNA提取试剂盒说明书提取植物总RNA, 用反转录试剂盒合成cDNA。根据甘蓝型油菜BnaCBF参考基因序列, 通过Primer premier 6.0软件设计qRT-PCR特异性引物, 引物序列见表1。采用Bio-Rad CFX96-M实时荧光定量PCR仪(Bio-Rad, 美国)检测目标基因的表达。参照荧光定量试剂盒说明书配置20 μL的反应体系, 含2×Trans start, Tip Green qPCR Super Mix 10 μL、cDNA模板2 μL、正反向引物(0.2 μmol L-1)各0.5 μL、去离子水7 μL。qRT-PCR反应程序为95℃ 5 min; 95℃ 10 s, 60℃ 45 s, 40个循环。选用BnaActin作为内参基因, 采用2-ΔΔCt方法分析基因表达变化情况。

表1 本研究所用qRT-PCR引物Table 1 Primers for qRT-PCR used in this study

2 结果与分析

2.1 甘蓝型油菜、白菜和甘蓝CBF基因鉴定

利 用4个AtCBF (At4g25490、At4g25470、At4g25480和AT5G51990)核苷酸和氨基酸序列分别在甘蓝型油菜、白菜和甘蓝的基因组数据库中进行同源序列搜索, 并结合SMART数据库、InterPro和NCBI数据库中的保守域(conservative domain, CD)搜索分析, 最终获得11个甘蓝型油菜CBF基因(BnaCBF)、6个白菜CBF基因(BraCBF)、5个甘蓝CBF基因(BoCBF)。其中11个甘蓝型油菜CBF基因编码区长度在576~663之间, 编码189~220个氨基酸, 预测等电点(Isoelectric point prediction, pI)在4.58~6.14之间, 为酸性蛋白, 分子量在20.57~24.56 kD之间。亚细胞定位预测表明, 这11个甘蓝型油菜CBF基因在细胞核中, 具体信息详见表2。

2.2 拟南芥、白菜、甘蓝及甘蓝型油菜CBF基因的序列比对分析

将鉴定得到的CBF家族蛋白的保守结构域通过DNAMAN软件比对(图1)发现, 蛋白序列保守性较高, 甘蓝型油菜、白菜、甘蓝及拟南芥的26个CBF序列比对一致性为61.08%, 这些基因都具有相同的保守序列AP2 domain (图1红色框出部分)。与其他成员相比,BnaCnng49280D、Bra010463和BnaA08g 30930D3个基因在AP2domain的前半部分有27个氨基酸序列缺失。

2.3 拟南芥、白菜、甘蓝及甘蓝型油菜的CBF基因进化树分析

为揭示甘蓝型油菜及其二倍体祖先种之间CBF家族基因的进化关系, 本研究利用MEGA7.0软件对拟南芥、白菜、甘蓝和甘蓝型油菜的CBF蛋白序列进行系统发育分析, 构建系统进化树(图2-A)。26个CBF蛋白被分为2个亚组, 与前人研究结果一致。亚组I中含有包括拟南芥AtCBF1/2/3在内的21个基因,其中, 5个白菜CBF基因、4个甘蓝CBF基因以及9个甘蓝型油菜CBF基因, 它们与拟南芥AtCBF1/2/3基因比对一致性分别为79.77%、68.38%和72.98%; 亚组II中含有包括拟南芥AtCBF4在内的5个基因, 其中2个来自甘蓝型油菜, 白菜和甘蓝各1个, 它们与拟南芥CBF4比对一致性分别为91.07%、80.36%和79.91%。进一步分析发现, 除白菜基因Bra010461外,其他白菜和甘蓝基因在甘蓝型油菜中均存在直系同源基因; 11个甘蓝型油菜CBF基因中, 除BnaA08g 30910D外, 其他10个基因都能在白菜和甘蓝中找到同源基因。

2.4 基因结构分析和MEME分析

为进一步揭示甘蓝型油菜、白菜和甘蓝及拟南芥CBF基因家族的进化关系, 本研究对其基因结构(图2-B)和motif进行了预测分析(图2-C)。CBF基因具有相对保守的基因结构和蛋白功能结构域,CBFs内含子的数目在0~2之间, 亚组I中, 有7个CBF基因有内含子, 其中白菜中2个CBF基因、甘蓝型油菜1个CBF基因有1个内含子和2个外显子, 甘蓝型油菜中4个CBF基因有2个内含子和3个外显子, 其他14个CBF基因只有1个外显子没有内含子, 与拟南芥AtCBF1/2/3基因结构一致, 进化树关系也更近; 亚组II中5个CBF基因没有内含子, 基因结构一致, 具有保守性。进一步分析发现, 亚组I中白菜中Bra010460、Bra010463含有1个内含子, 甘蓝型油菜中与这2个基因进化同源关系BnaA08g30950D和BnaA08g30930D含有2个内含子, 因此推断这4个基因在进化过程中可能发生了变异; 另一方面, 甘蓝型油菜中BnaUnng01150D有1个内含子,BnaCnng49280D有2个内含子, 而它们对应的同源基因Bo8g050800和Bo8g050740没有内含子, 推测这2个甘蓝型油菜CBF基因进化过程可能发生了变异;同时我们在进化树分析中推断BnaA08g30910D可能是BnaA08g30950D的拷贝基因, 两者基因结构一致。

?

Motif预测分析结果表明, Motif预测分析与进化树比对分析结果吻合, 可明显分为2个亚组, 2组保守基序组成高度相似。图2-C展示了保守性较高的6条基序在26条CBF蛋白序列中的分布情况, 可见各蛋白高度保守, 但不同亚组之间蛋白序列包含的保守基序的种类和数目存在一定的差异, 2个亚组内蛋白保守基序的组成模式更为相似。26条CBF蛋白序列都含有Motif 6、Motif 1、Motif 4和Motif 3, 由此推断这些区域可能是CBF家族蛋白共性的、保守的功能基序。进一步分析发现, 亚组I中的22个成员保守基序组成模式较为一致, 除了亚组I中BnaUnng01150D、BnaA08g30930D不含有Motif5, 其余成员保守基序组成模式与拟南芥CBF1/2/3蛋白保守基序完全一致, 均为Motif6-Motif2-Motif1-Motif4-Motif5-Motif3。亚组II中, 所有成员保守基序组成模式完全一致, 均为Motif6-Motif2-Motif1-Motif4-Motif3; 植物中的CBF家族蛋白含有1个高度保守的AP2domain, AP2domain包括Motif2-Motif1-Motif4保守基序。与亚组I相比, 亚组II少了1个保守基序Motif5,由此推测Motif5为亚组I的特征Motif。

2.5 染色体定位分析

ATCBF1/2/3基因在Chr4染色体上串联,AtCBF4在Chr5染色体上(图3); 6个BrCBF基因分别在A03、A08和A10染色体上分布, 其中Bra010460、Bra010461和Bra010463在A08号染色体上串联(图3); 6个BoCBF基因被定位在C03、C07、C08、C09染色体上, 其中Bo8g050800和Bo8g050740在C08染色体上串联(图3); 根据Region Comparison(http://plants.ensembl.org/)分析, 我们对每个CBF成员之间的共线性进行比较分析发现,BnaA08g 30950D和BnaA03g13620D与AtCBF1同 源,BnaUnng01150D、BnaA08g30910D、BnaCnng492 80D、BnaC03g71900D、BnaAnng-34260D和BnaC07g 39680D与AtCBF2同源,BnaA08g30930D与AtCBF3同源,BnaA10g07630D和BnaC09g28190D与AtCBF4同源,这些BnaCBFs基因被定位到19条染色体中的9条上(图3), 其中BnaA08g30910D、BnaA08g30930D和BnaA08g30950D在染色体A08 random上串联重复。

2.6 甘蓝型油菜BnaCBF基因在不同组织中的表达模式分析

利用本课题组前期完成的中双11材料的转录组数据, 分析了甘蓝型油菜CBF基因家族在各组织器官中表达情况, 以及在不同非生物胁迫和激素处理下的表达情况。在正常情况下, 亚组I中9个BnaCBF成 员 可 分 为Ia类 和Ib类, 其 中Ia类5个CBF成 员(BnaA08g30950D、BnaUnng01150D、BnaA08g309 10D、BnaA08g30930D和BnaCnng49280D)在各个组织器官中表达量均较低或者不表达; Ib类4个成员(BnaA03g13620D、BnaC03g71900D、BnaAnng34260D和BnaC07g39680D)分别在茎、叶片和雌蕊中表达,其中BnaAnng34260D在雌蕊中表达量较高; 亚组II中 的2个CBF4成 员(BnaA10g07630D和BnaC 09g28190D)在根、叶片、花萼和雌蕊中均有表达, 并且在根部表达量明显高于其他组织器官(图4)。

2.7 甘蓝型油菜BnaCBF基因在不同非生物逆境和不同激素处理表达模式分析

所有BnaCBF在逆境处理下均存在不同程度的表达。在低温下, Ia类5个基因(BnaA08g30950D、BnaUnng01150D、BnaA08g30910D、BnaA08g30930D和BnaCnng49280D)对冷胁迫迅速响应表达, 并在3 h后表达量显著下降, 意味着这5个基因是早期冷诱导表达基因。Ib类4个基因(BnaA03g13620D、BnaC03g71900D、BnaAnng34260D和BnaC07g39680D)对冷胁迫响应迅速且持续时间长, 冷处理12 h仍有较高表达量; 而亚组II两个CBF4基因(BnaA 10g07630D和BnaC09g28190D)在冷胁迫条件下表现为早期高表达, 而后降低, 6 h后恢复正常。BnaA10g07630D和BnaC09g28190D两个基因在根中表达量明显高于叶片, 并且参与盐胁迫和渗透胁迫响应, 而其他9个基因在根中表达变化不大(图5)。值得注意的是, Ib类4个CBF基因的表达表现出较明显的生物钟节律变化, 它们在11:00 (Leaf_3h)表达量较低, 14:00表达量较之前升高, 20:00表达量几乎为零, 表现出明显的昼夜节律; 亚组II中的2个CBF4基因(BnaA10g07630D和BnaC09g28190D)也呈现出明显的昼夜节律, 不同于Ib类CBF基因, 其在11:00和14:00几乎不表达, 但在20:00有较高的表达量。其他ABA、MeJA和ET等激素处理后表达量与对照无明显变化。

2.8 不同低温处理下甘蓝型油菜BnaCBF基因表达模式分析

为进一步分析这些CBF成员对不同低温的响应,本研究利用本课题组前期完成的冷驯化和非冷驯化材料中不同低温处理的HX-17和HX-58材料的转录组数据, 对甘蓝型油菜BnaCBF基因家族的表达情况进一步分析(图6)。亚组I中有3个CBF基因在冷胁迫下表达较弱甚至不表达, 冻害胁迫下表达量升高;其他6个成员在4℃冷害处理和-4℃冻害胁迫下均有不同程度响应, 且-4℃冻害胁迫下表达量更高。亚组II 2个CBF4成员仅在冻害胁迫条件响应表达, 意味着甘蓝型油菜中CBF基因家族对冻害响应更强烈。分析发现, 冷驯化和非冷驯化的材料对低温胁迫和冻害胁迫的响应表达的影响无明显差异。

2.9 BnaCBF基因相对表达量结果分析

为验证上图中的RNA-seq数据, 本研究从油菜CBF基因家族亚组Ia、Ib和II中挑选4个具有代表性的基因进行qRT-PCR验证, 分别是亚组Ia类基因BnaCnng49280D, Ib类基因BnaAnng34260D以及亚组II的BnaA10g07630D和BnaC09g28190D。结果显示, 表达情况与图5和图6中完全吻合, 说明表达谱数据可靠(图7和图8)。

3 讨论

甘蓝型油菜是我国最主要的栽培品种, 在湖南、湖北、江西、浙江、云南等多地大面积种植[4]。南方地区推广的稻–稻–油轮作生产模式需要种植早熟、特早熟油菜品种, 但早熟型油菜推广面临抗冻性差、产量受影响等问题[2]。CBF基因家族在植物抗冻性研究中具有重要作用[28], 因此, 系统的比较分析甘蓝型油菜CBF基因以及在油菜二倍体祖先中的进化关系, 对解析CBF基因调控油菜低温及冻害胁迫反应的机制具有重要的意义。

CBF/DREB1依赖的冷信号通路在过去的20年中取得了长足的进展。从最早发现的一种对干旱、寒冷和高盐胁迫有响应的顺式作用元件C-反应/脱水响应元件(CRT/DRE)开始[29], 研究者陆续从拟南芥、油菜、小麦(Triticum aestivum)、黑麦(Secale cereale)、番茄(Lycopersicon esculentum)、水稻(Oryza sativa)及玉米(Zea mays)等植物中分离克隆出许多CBF的ESTs序列, 发现其都具有冷诱导特性[30-32],过表达拟南芥及其他CBF基因均可以增强植物的抗冻性[12,19]。大量研究证实, 在植物中这些CBF是参与冷胁迫信号传导, 并在不同物种间具有保守性。

本研究发现,AtCBF在甘蓝型油菜、白菜和甘蓝基因组都存在直系同源基因, 且这些CBF成员序列高度保守, 根据基因结构和Motif等结果分析, 推测BnaA08g30910D可能是BnaA08g30950D的拷贝复制。甘蓝型油菜CBF家族蛋白的结构域保守性较高,所有成员都有且只有1个AP2 domain, 在此结构域的上下游, 各有一段保守的氨基酸序列PKKP/PKK PAGR (RAGRxxKFx ETRHP)和DSAWR[33], 该保守序列在CBF发挥其转录因子功能中扮演着非常重要的角色。

基因中内含子和外显子的组成对于其进化分析研究具有非常重要的作用[34], 本研究对CBF基因家族各成员的基因结构分析推断含有内含子的CBF基因在进化过程中可能发生了变异。Motif预测分析发现, 所有的CBF基因中Motif1、Motif3、Motif4和Motif6这4个Motif是CBF基因家族的特征性Motif。保守基序分析发现,BnaCnng49280D、BnaA08g 30930D和Bra010463三个基因缺失的27个氨基酸序列位于AP2 domain上游, 结构域分析发现该缺失部分在Motif2区段, 由此推测可能会影响它们结合靶基因启动子中顺式元件CRT/DRE的能力, 从而影响它们激活下游靶基因(如COR)的表达。染色体定位分析表明, 11个甘蓝型油菜CBF基因在染色体上的分布是不均匀的, 这可能是由于在全基因组加倍过程中,一些同源基因可能发生丢失或扩张。

CBF参与多种逆境响应和植株生长发育等过程[35-36]。本研究发现, 甘蓝型油菜CBF基因的表达强烈受低温诱导, 且具有明显的组织特异性, 由此推断甘蓝型油菜CBF基因可能在功能上产生细化。A亚基因组BnaA03g13620D、BnaAnng34260DCBF的表达量明显高于对应的C亚基因组BnaC03g 71900D、BnaC07g39680D的表达量, 说明A、C亚基因组的基因对表达量存在明显偏移不对称。

Haake等[11]研究表明,CBF1/2/3(CBF1/DREB1b、CBF2/DREB1c、CBF3/DREB1a)属于冷诱导CBF基因, 还有一个CBF基因(CBF4)不受冷诱导表达; 但过表达CBF4会增强植物的抗冻性和耐旱性的。本研究发现, 11个BnaCBF中, 与AtCBF1/2/3同源的9个基因对冷胁迫迅速响应, 并且冻害胁迫比冷胁迫响应更为强烈;AtCBF4同源的2个BnaCBF基因BnaA10g07630D和BnaC09g 28190D轻微冷胁迫响应, 但强烈响应冻害胁迫。由于前人主要以非典型的越冬植物如水稻、拟南芥、番茄等作为研究材料, 相关冻害研究开展的较少, 但由于甘蓝型油菜本身是越冬性作物, 生育期通常会经历一段较长的低温冻害时期, 因此BnaCBFs基因在低温时的诱导表达没有冻害胁迫下强烈。另外, 本研究中发现, 冷驯化的冻害和非冷驯化冻害下表达差异不大, 说明BnaCBFs基因驯化和非冷驯化条件下均起作用。

4 结论

本研究鉴定出甘蓝型油菜中11个CBF基因家族成员, 并将BnaCBFs分为亚组I (CBF1/2/3)和亚组II (CBF4) 2个亚组。所有甘蓝型油菜CBF基因受低温诱导表达, 其中亚组Ib的4个基因BnaA03g 13620D、BnaC03g71900D、BnaAnng34260D、BnaC07g 39680D对冷胁迫响应迅速且持续时间长, 亚组II中2个基因BnaA10g07630D和BnaC09g28190D对冷胁迫响应表达较弱, 但它们在根中表达量明显高于叶片, 并且参与盐胁迫和渗透胁迫响应; BnaCBF对冻害响应更强烈。本研究将为进一步了解甘蓝型油菜CBF家族基因的生物学功能、低温胁迫尤其是冻害胁迫研究及耐寒、抗冻品种培育奠定了基础。