甘蓝型油菜P5CR同源基因的克隆表达及其多态性分析

马秀奇，张晓娟 ，孙晓敏，陈 乔，张 羽，何开罡，权 莹，宋建民

（1.陕西理工大学生物科学与工程学院，陕西 汉中 723001；2.汉中市农业技术推广与培训中心，陕西 汉中723001；3.汉中职业技术学院，陕西 汉中 723001；4.陕西省资源生物重点实验室，陕西 汉中 723000）

0 引言

【研究意义】油菜作为重要的油料作物，是我国乃至世界范围内食用油供给的主要来源。油菜除了油用，还可作为蔬菜、饲料、观赏作物及蜜源，具有重要的经济价值。我国油菜种植面积及分布居油料作物之首[1]。因此，油菜安全生产对国民经济意义重大。由真菌-核盘菌引起的菌核病给油菜生产带来严重威胁，引起油菜产量及品质大幅度降低[2,3]。我国长江中下游地区气候湿润，加之近些年油菜种植密度增加，加重了该病害的发生[4]。克隆油菜抗菌核病基因，研究油菜抗菌核病分子机理，对通过分子设计育种途径选育性状优良的抗病油菜新品种意义重大[5]。【前人研究进展】吡咯林-5-羧酸还原酶（Pyrroline-5-carboxylate reductase，P5CR ）是真核生物中重要的一种管家蛋白，催化脯氨酸生物合成的最后一步反应，负责在细胞质中将脯氨酸、精氨酸、鸟氨酸和谷氨酸中间代谢产物吡咯啉-5-羧酸（Δ1-pyrroline-5-carboxylate，P5C）还原成为脯氨酸，在氨基酸代谢的调节中发挥重要作用。实验室前期通过大豆菌核病全基因组关联分析得出P5CR基因与菌核病抗性强关联[6]，推测该基因可能与油菜菌核病抗性有关。Senthil-Kuma等研究表明鸟氨酸氨基转移酶（δOAT）和脯氨酸脱氢酶（ProDH1和ProDH2）通过调控P5C代谢诱导的超敏反应进而参与非寄主抗性[7]。周婉莹等对甘蓝型油菜A10、C09上的P5CR同源基因进行克隆及多态性分析，结果得出，A10染色体中的P5CR同源基因上的多态性位点与菌核病抗性没有显著关联，C09染色体上P5CR同源基因在感病材料中存在一段A10染色体上P5CR同源基因片段的插入，导致前者ORF改变，进而可能造成该基因抗病性丧失。推测C09染色体中的P5CR同源基因与油菜对菌核病的抗性相关[8]。【本研究切入点】目前，关于P5CR基因与油菜抗菌核病关系研究报道较少。甘蓝型油菜A03、C03染色体上的P5CR同源基因与油菜菌核病抗性关系及抗性相关多态性位点尚不清楚。【拟解决的关键问题】本研究以对菌核病具有较好抗性的中油821、沪16油菜以及感病油菜川秦6R及南12R为研究材料，设计特异引物，对甘蓝型油菜C03、A03染色体上的P5CR同源基因进行PCR扩增、克隆测序及表达研究，分析P5CR基因中的SNP位点及其与油菜菌核病抗性的相关性，为油菜菌核病抗性基因克隆、油菜抗菌核病机制揭示及通过分子标记辅助育种选育抗病油菜新品种提供理论依据。

1 材料与方法

1.1 材料

本研究选用高抗、高感菌核病油菜为试验材料（表1），抗病材料记为R，感病材料记为S。油菜材料于2021年种植于汉中勉县黄沙镇试验田（106°7＇91.26″E，33°13＇77.72″N）。大肠杆菌 DH5α（TaKaRa，中国大连）由本实验室保存。核盘菌（Sclerotinia sclerotiorum）由汉中市农业技术推广与培训中心惠赠。TransStart fast Pfu fly DNA polymerse（全式金公司，北京）、植物总RNA提取试剂盒（庄萌生物，北京）、EvoM-MLV反转录试剂盒II（艾科瑞生物，湖南）、CTAB、琼脂糖、异丙醇、无水乙醇、Goldview核酸染料购自西安沃尔森公司。

表1 油菜样品的编号、品种名称和表现型Table 1 Codes, variety names, and phenotypes of B.napus samples

1.2 油菜基因组DNA提取及检测

采集供试油菜材料幼嫩叶片。采取改良CTAB法[9]提取油菜基因组DNA。通过1%的琼脂糖凝胶检测基因组DNA完整性。

1.3 引物设计、PCR扩增及克隆测序

从拟南芥基因组数据库（https://www.arabidopsis.org/）中下载拟南芥P5CR基因（AT5G14800）蛋白序列，通过NCBI 数据库（https://www.ncbi.nlm.nih.gov/gene/）中的Blast软件，在甘蓝型油菜数据库中进行同源基因比对（比对指标为序列覆盖度达90%以上，同源性大于80%，e-value小于1e-5）。针对甘蓝型油菜C03及A03染色体上P5CR同源基因设计特异引物（表2）。引物由西安沃尔森科技有限公司合成。

表2 甘蓝型油菜C03及A03染色体上P5CR同源基因引物序列信息Table 2 Primer sequence of P5CR on chromosomes C03 and A03 of B.napus

PCR反应体系共25 μL，其中TransStart fast Pfu fly DNA polymerse 0.5 μL， 5×TransStrat Fast Pfu Fly buffer 5 μL， 2.5 mmol·L-1dNTP 2 μL， 50 ng·μL-1DNA 模板 0.5 μL，ddH2O 17 μL。PCR 预变性采取 94℃ 4 min，接着94 ℃变性1 min，然后58 ℃退火30 s，最后72 ℃延伸2 min，循环35次，72 ℃延伸7 min后4 ℃保存。PCR产物经1%琼脂糖凝胶电泳分离后胶回收目标产物，并与PMD19-T载体（TaKaRa）连接，转化大肠杆菌DH5α，将转化体系涂布于LB固体培养基上（含IPTG 200 mg·L-1，XGal 20 mg·L-1） 37 ℃ 倒 置 培 养 过 夜 ， 进 行 菌 落PCR鉴定，将阳性克隆送测序。

1.4 序列分析

高抗、高感菌核病油菜材料P5CR同源基因测序结果通过DNASTAR Lasergene中的Megalign软件进行序列比对，分析抗、感材料之间的多态性位点。同时，从NCBI网站中甘蓝型油菜数据库中下载甘蓝型油菜C03及A03染色体上P5CR基因的参考序列，与测序结果进行多重序列比对。利用DNAMAN软件将P5CR基因翻译成蛋白序列，并对抗、感油菜材料中P5CR基因编码蛋白序列进行比对分析。

1.5 核盘菌培养、接种及样品采集

将核盘菌菌核置于马铃薯葡萄糖琼脂（PDA）培养基上22 ℃暗培养4 d。接着，将边缘的菌丝转移至新的PDA培养基上生长。2 d后，菌丝体直径约为6 mm，用于接种试验。于油菜初花期，采集抗菌核病油菜材料中油821及感病材料南12R相同部位相同大小的叶片，通过叶片离体接种法对其进行核盘菌接种[10]。分别采集接种前（即接种0 h）及接种后6 h、12 h、24 h及48 h抗、感菌核病油菜叶片，迅速放入液氮中速冻后保存于-80 ℃冰箱，每个样品取3个重复。

1.6 RNA提取及反转录

总RNA提取采用植物总RNA试剂盒提取（提取步骤参照试剂盒说明书），通过琼脂糖凝胶电泳检测RNA完整性，采用Nanodrop核酸蛋白检测仪测定RNA浓度。采取Evo M-MLV反转录试剂盒II进行反转录。用ddH2O将合成的cDNA稀释10倍，冷冻保存备用。

1.7 qPCR检测P5CR基因表达

根据A03及C03上P5CR同源基因序列分别设计特异引物，通过定量PCR分析A03及C03染色体上P5CR基因在抗、感菌核病油菜材料接种前及接种后6 h、12 h、24 h、48 h的表达量。本研究中的内参基因为Ubc21。P5CR基因及内参引物序列信息见表2。

2 结果与分析

2.1 甘蓝型油菜A03及C03染色体上P5CR同源基因的扩增

将拟南芥P5CR基因在甘蓝型油菜基因组数据库中比对，共得到4个P5CR同源基因，分别位于A10、C09、C03及A03染色体上（基因序列号分别为106 372 121、106 416 421、106 417 090、106 431 715）。课题组前期已对A10、C09染色体上P5CR同源基因进行了克隆[8]，本研究对A03及C03染色体上P5CR基因进行克隆及序列分析。

在抗、感菌核病油菜材料中，利用引物P5CRA03F1/P5CRA03R1对A03染色体上的P5CR同源基因进行扩增，得到一条2000 bp左右的特异性条带，与预期条带大小（1977 bp）相符（图1A）。接着，对PCR产物进行胶回收、克隆及测序。

利用引物P5CRC03F2/P5CRC03R2对抗、感菌核病油菜材料C03染色体上P5CR同源基因进行扩增，扩增产物为大小介于1 200～2 000 bp的特异性条带（图1B），与本研究预期扩增产物大小（1 889 bp）相符。对PCR产物进行胶回收、克隆测序。

2.2 抗、感菌核病油菜P5CR基因多态性分析

抗、感菌核病油菜材料A03染色体上P5CR基因扩增产物经克隆测序、序列分析得出，该基因全长1 526 bp，有6个内含子，7个外显子。该基因在抗、感菌核病油菜材料中共存在15个多态性位点，分别位于A03染色体上P5CR同源基因的21、34、70、130、153、234、235、580、631、736、766、803、841、1 546及1 574位。其中，234位在感病油菜材料中为C、在抗病油菜材料中为A，235位在感病油菜中为C，与参考基因序列一致，在抗病油菜中突变为T。推测这2个多态性位点与油菜菌核病抗性相关（图2）。

抗、感菌核病油菜材料C03染色体上P5CR基因扩增产物经克隆测序、序列分析得出，该基因全长1 457 bp，有6个内含子，7个外显子。抗、感材料中该基因共存在7个多态性位点，分别位于该基因的 564、579、636、671、673、676、679 位。其中3个多态性位点与材料的抗、感性状相关（图3），即抗病材料中该基因564、636及679位由C突变为T，推测这三个位点与油菜抗病性相关。

A03染色体上P5CR基因编码蛋白序列分析得出，该基因编码蛋白质含276个氨基酸。在抗、感菌核病油菜中存在多态性的位点234及235均位于该基因第一个内含子中，不引起编码氨基酸序列改变。抗、感菌核病油菜材料C03染色体上P5CR基因编码蛋白序列分析得出，该基因编码蛋白质含144个氨基酸。分析抗、感材料中可能与油菜菌核病抗性相关的位点得出，抗性材料的P5CR基因564位SNP为同义突变，不引起编码氨基酸改变，其余两个SNP位点（636及679位）位于内含子中，不编码蛋白质。

2.3 qPCR分析抗、感菌核病油菜中P5CR基因表达

为进一步分析P5CR基因是否参与油菜对菌核病的抗性，通过qPCR分析A03及C03染色体上P5CR基因在抗、感菌核病油菜材料中接种核盘菌前及接种后不同时间的表达。结果显示，A03上的P5CR基因在接种前后的抗、感菌核病油菜材料中表达模式不同。抗病油菜材料中，A03上的P5CR基因在接种前后变化明显，在接种后24 h该基因的表达量最高，达到未接种时的2倍多。感病材料中，A03上的P5CR基因在接种前后表达量差异不大（图4A）。C03染色体上的P5CR基因在接种前后抗、感菌核病油菜材料中同样存在不同的表达模式，其在抗病材料中接种后24 h表达量最高，为未接种时的3倍多。感病材料中，C03染色体上的P5CR基因在接种后6 h、12 h表达量相对于未接种时略有所提高（图4B）。综上，A03及C03染色体上的P5CR基因在抗病材料接种后24 h表达量显著升高，说明该基因可能参与了油菜对菌核病的抗性。感病材料中，C03染色体上的P5CR基因虽然在接种后6 h、12 h表达量相对于未接种时略有所提高，但整体变化不大，导致其对菌核病未能表现出较好的抗性。

3 讨论

油菜是世界范围内广泛种植的油料作物。油菜菌核病是世界范围内尤其是温带地区普遍发生的一种严重的油菜病害。我国的东南沿海地区及长江流域，该病的田间发病率较高，导致油菜产量降低10%～30%，甚至高达80%，给油菜生产造成很大损失[11]。近年来，由于油菜机械化高密度种植，该病有不断加重的趋势。克隆油菜抗菌核病基因，开发抗病基因分子标记，通过分子标记辅助育种方式选育抗病油菜新品种是防治油菜菌核病最根本有效的途径。

P5CR是脯氨酸生物合成途径的关键酶。脯氨酸作为一种抗渗透胁迫调节剂，在植物渗透调节中发挥重要作用。研究表明，脯氨酸的积累是植物对各种胁迫（包括盐和渗透胁迫）的普遍反应，脯氨酸累积量与植物的胁迫耐受能力呈正相关[12,13]。首个P5CR基因在大豆中被克隆到，该基因全长1.2 kb，编码蛋白为单体蛋白[14]。之后相继在拟南芥、水稻、黑麦草、木薯等植物中分离到P5CR基因[15-18]。付莉莉等（2016）克隆了木薯中P5CR基因，对该基因在干旱胁迫下的表达分析表明，P5CR基因表达受胁迫调控。P5CR基因还响应低温、盐和ABA等胁迫条件[19]。在低温胁迫条件下，AtP5CR的表达量显著提高[20]。拟南芥中，P5CR基因表达受盐胁迫调控[21]。Cao等（2015）克隆了黑麦草（Lolium perenneL.）P5CR基因，其全长1 047 bp，受盐、PEG、低温和ABA处理诱导[15]。

目前，关于P5CR基因与植物抗病性的研究报道较少。甘蓝型油菜基因组中共有4个P5CR同源基因，课题组前期克隆了C09、A10上P5CR同源基因。本研究对C03、A03染色体上的P5CR同源基因进行克隆测序及多态性分析，研究得出甘蓝型油菜C03染色体上P5CR同源基因全长1 457 bp，抗、感油菜材料中，C03染色体上P5CR同源基因共存在7个多态性位点，其中3个多态性位点（564、636及679位）推测与油菜抗病性相关。A03染色体上的P5CR同源基因全长1 526 bp，其在抗、感菌核病油菜材料中共存在15个多态性位点，其中，234、235位2个多态性位点推测与油菜菌核病抗性相关。进一步，P5CR基因编码蛋白质预测分析显示，抗性材料中C03染色体上P5CR基因564位SNP为同义突变，不引起编码氨基酸改变，其余两个SNP位点（636及679位）位于内含子中。A03上与菌核病抗性相关的两个位点234、235位均位于内含子中。通过qPCR分析A03、C03染色体上P5CR在抗、感油菜接种前及接种后不同时间的表达量。结果显示，在抗病材料中，A03及C03染色体上的P5CR基因在接种后24 h表达量均明显提高，说明A03、C03染色体上的P5CR基因均可能与油菜菌核病抗性相关。结合基因测序及定量PCR结果，本研究得出P5CR基因可能是通过基因表达水平的变化参与油菜抗菌核病反应。下一步，研究将扩大样本量，筛选出更多的高抗、高感菌核病油菜材料，对P5CR基因进行克隆及测序，从而进一步验证抗、感菌核病油菜材料中P5CR基因上的多态性位点与油菜菌核病抗性相关性。同时，研究将通过转基因、RNA干扰等实验对P5CR基因进行功能验证。

本研究为揭示油菜抗菌核病分子机理、油菜抗菌核病分子标记开发及油菜抗病新品种选育奠定基础。