解析白菜型油菜BrPARP1基因在干旱胁迫应答中的功能与机制

解析白菜型油菜BrPARP1基因在干旱胁迫应答中的功能与机制一、引言1.1研究背景与意义油菜作为全球重要的油料作物之一，在农业生产中占据着举足轻重的地位。它不仅是食用植物油的主要来源，其饼粕还富含蛋白质，是优质的饲料原料。我国是油菜的主要生产国与消费国，油菜播种面积占油料作物播种总面积的50.31%，菜籽油占国产油料作物产油量的57%以上。随着全球人口的增长以及对植物油需求的不断攀升，油菜的种植面积和产量也在持续扩大。然而，干旱胁迫已成为制约油菜生长发育和产量提升的主要环境因素之一。随着全球气候变化，干旱发生的频率和强度呈上升趋势，给油菜生产带来了严峻挑战。干旱胁迫对油菜的影响贯穿其整个生长发育周期。在种子萌发阶段，干旱会降低种子的吸水率，影响种子内的生理生化反应，导致萌发率下降，出苗延迟且不整齐。在幼苗期，干旱抑制根系的生长和发育，使根系变浅、变细，降低根系对水分和养分的吸收能力，同时导致叶片气孔关闭，蒸腾作用减弱，光合作用受到阻碍，进而影响幼苗的生长速度和生物量积累。在生殖生长阶段，干旱影响油菜的花芽分化、开花授粉和角果发育，导致花器官发育异常，授粉受精不良，角果数减少，籽粒饱满度降低，最终使油菜的产量大幅下降。据相关研究报道，中度干旱可导致油菜单产降低20％以上，严重时甚至绝收。植物在长期进化过程中，形成了一系列复杂而精细的抗旱机制来应对干旱胁迫。其中，基因表达调控在植物抗旱反应中起着核心作用。众多基因参与到油菜对干旱胁迫的响应过程中，它们通过编码各种功能蛋白和调节蛋白，如转录因子、蛋白激酶、渗透调节物质合成酶等，来调节植物的生理生化过程，增强植物的抗旱能力。对这些基因的深入研究，有助于我们揭示油菜抗旱的分子机制，为油菜抗旱品种的选育提供理论基础。聚[ADP-核糖]聚合酶1（PARP1）作为一种在生物体内广泛存在的核酶，在DNA损伤修复、基因组稳定性维持、细胞凋亡等多种生物学过程中发挥着关键作用。近年来，越来越多的研究表明，PARP1也参与到植物对逆境胁迫的响应过程中。在干旱胁迫下，PARP1可能通过调节相关基因的表达，影响植物的渗透调节、抗氧化防御等生理过程，从而增强植物的抗旱能力。在油菜中，白菜型油菜具有早熟、抗寒、抗旱、耐瘠等优异特性，在我国高寒冷凉地区仍然是不可替代的主要农作物和特色经济作物。然而，目前对于白菜型油菜中BrPARP1基因在干旱胁迫应答中的功能研究还相对较少。深入探究BrPARP1基因在白菜型油菜干旱胁迫应答中的功能，不仅有助于揭示油菜抗旱的分子机制，丰富植物逆境生物学的理论知识，还能为油菜抗旱品种的选育提供关键的基因资源和技术支持。通过基因工程等手段，将BrPARP1基因导入油菜品种中，有望培育出具有更强抗旱能力的油菜新品种，从而减少干旱对油菜生产的损失，保障食用油的稳定供应，促进农业的可持续发展。1.2白菜型油菜概述白菜型油菜（BrassicarapaL.），在油菜分类体系中占据着独特的位置，属于十字花科芸薹属。它是一种古老的油菜类型，历经长期的自然选择与人工驯化，具备丰富的遗传多样性。其植株通常较为矮小，茎秆相对细弱，叶片大且呈薄纸质，无明显的蜡粉覆盖。白菜型油菜在我国的种植历史源远流长，最早可追溯至数千年前，最初多作为蔬菜或观赏植物栽培，后来逐渐向油料作物转变。在全球范围内，白菜型油菜的种植分布广泛，涵盖了亚洲、欧洲、北美洲等多个地区。在我国，它主要集中种植于青海、甘肃、新疆、西藏等高寒冷凉地区。这些地区海拔较高，气候条件复杂，冬季漫长且寒冷，夏季短暂且凉爽，年降水量相对较少，土壤肥力也较为有限。而白菜型油菜凭借其早熟、抗寒、抗旱、耐瘠等一系列优异特性，能够在这样恶劣的环境中良好生长。以青海为例，作为全国最大的白菜型油菜主产区之一，其独特的高原大陆性气候，昼夜温差大，光照充足，为白菜型油菜的生长提供了得天独厚的条件。这里种植的白菜型油菜，不仅产量可观，而且品质优良，含油量较高。白菜型油菜具有极高的经济价值，它的种子含油量丰富，一般在35%-45%之间，是优质的食用油原料。菜籽油富含多种不饱和脂肪酸，如油酸、亚油酸等，这些脂肪酸对于降低人体胆固醇、预防心血管疾病具有重要作用。同时，白菜型油菜的饼粕富含蛋白质、氨基酸等营养成分，是优质的饲料原料，广泛应用于畜牧业生产中，为养殖业的发展提供了有力支持。此外，白菜型油菜还具有一定的观赏价值，在花期，大片的油菜花田金黄灿烂，吸引了众多游客前来观赏，促进了当地旅游业的发展，带动了相关产业的繁荣。在农业生产中，白菜型油菜也占据着重要地位。它是高寒冷凉地区的主要农作物和特色经济作物，对于保障当地的食用油供应、促进农民增收以及维护区域生态平衡都发挥着不可替代的作用。在青海、甘肃等地，白菜型油菜的种植已成为当地农业经济的重要支柱产业之一，为当地农村经济的发展注入了强大动力。同时，由于白菜型油菜具有较强的适应性和抗逆性，能够在其他作物难以生长的环境中种植，这对于充分利用土地资源、提高土地利用率具有重要意义。1.3BrPARP1基因简介BrPARP1基因是聚[ADP-核糖]聚合酶（PARP）家族的重要成员，在植物应对干旱胁迫的过程中发挥着潜在的关键作用。该基因编码的蛋白属于PARP家族，其结构具有高度的保守性。BrPARP1蛋白通常包含多个功能结构域，如N端的DNA结合结构域、中部的催化结构域以及C端的调节结构域。DNA结合结构域能够识别并特异性地结合到受损的DNA分子上，这是BrPARP1发挥功能的起始步骤。当DNA受到干旱等逆境胁迫而发生损伤时，该结构域可迅速感知并与之结合，从而激活整个酶蛋白的活性。中部的催化结构域则是BrPARP1的核心功能区域，它以烟酰胺腺嘌呤二核苷酸（NAD+）为底物，催化聚[ADP-核糖]（PAR）的合成。PAR是一种线性或分支状的多聚物，其合成过程会消耗大量的NAD+，进而影响细胞内的能量代谢和物质代谢。C端的调节结构域则对BrPARP1的活性起着精细的调控作用，它可以与其他蛋白质相互作用，通过蛋白质-蛋白质之间的相互联系，调节BrPARP1在细胞内的定位、活性以及参与的信号通路。在植物界中，BrPARP1基因具有较高的保守性。通过对多种植物的基因组序列分析发现，不同植物中的PARP1基因在核苷酸序列和氨基酸序列上都存在一定程度的相似性。以拟南芥、水稻等模式植物为例，它们的PARP1基因与白菜型油菜的BrPARP1基因在关键结构域的氨基酸序列相似度可达70%-80%以上。这种保守性暗示了PARP1基因在植物进化过程中具有重要的生物学功能，并且在不同植物物种中可能通过相似的分子机制参与到生长发育和逆境响应等生理过程中。已有研究表明，在其他植物中PARP1基因参与了多种生物学过程。在拟南芥中，PARP1基因被发现参与了DNA损伤修复过程。当拟南芥受到紫外线、化学诱变剂等因素导致DNA损伤时，PARP1能够迅速被激活，通过催化PAR的合成，招募相关的DNA修复蛋白到损伤位点，促进DNA的修复，从而维持基因组的稳定性。在水稻中，PARP1基因还被报道与植物的生长发育调控相关。干扰水稻中PARP1基因的表达，会导致水稻植株出现生长迟缓、叶片发育异常等表型，表明PARP1基因在水稻的正常生长发育过程中起着不可或缺的作用。基于其他植物中PARP1基因的功能研究，推测白菜型油菜中的BrPARP1基因在干旱胁迫应答中可能具有多种功能。一方面，在干旱胁迫下，植物细胞内的DNA容易受到损伤，BrPARP1基因可能通过激活DNA损伤修复机制，维持基因组的完整性，从而保证细胞的正常生理功能和植物的生长发育。另一方面，BrPARP1基因可能通过调节相关基因的表达，参与植物的渗透调节、抗氧化防御等生理过程。在渗透调节方面，它可能调控一些编码渗透调节物质合成酶的基因表达，促进脯氨酸、甜菜碱等渗透调节物质的合成和积累，降低细胞的渗透势，提高细胞的保水能力。在抗氧化防御方面，BrPARP1基因可能影响抗氧化酶基因的表达，增强超氧化物歧化酶（SOD）、过氧化物酶（POD）等抗氧化酶的活性，清除细胞内过多的活性氧（ROS），减轻氧化损伤，从而增强白菜型油菜对干旱胁迫的耐受性。1.4植物应对干旱胁迫的机制研究进展干旱胁迫是一种严重影响植物生长发育、地理分布和农作物产量的非生物胁迫。在长期的进化过程中，植物逐渐形成了一系列复杂而精细的应对干旱胁迫的机制，这些机制涉及植物的生理生化、分子水平等多个层面。从生理生化角度来看，植物主要通过渗透调节、抗氧化系统和激素调节等方式来应对干旱胁迫。渗透调节是植物适应干旱胁迫的重要机制之一，植物在干旱条件下，会主动积累一些小分子有机化合物和无机离子，如脯氨酸、甜菜碱、可溶性糖、K+等，这些物质被称为渗透调节物质。它们能够降低细胞的渗透势，促进细胞吸水，维持细胞的膨压，从而保证细胞的正常生理功能。以脯氨酸为例，它不仅是一种有效的渗透调节物质，还具有稳定蛋白质和生物膜结构、清除活性氧等作用。在干旱胁迫下，植物体内脯氨酸的合成关键酶活性增强，导致脯氨酸大量积累，帮助植物抵御干旱伤害。抗氧化系统是植物应对干旱胁迫的另一道重要防线。干旱胁迫会导致植物体内活性氧（ROS）如超氧阴离子（O2-）、过氧化氢（H2O2）和羟自由基（・OH）等的积累，这些ROS具有很强的氧化活性，能够攻击细胞膜、蛋白质、核酸等生物大分子，导致细胞损伤和死亡。为了清除过多的ROS，植物进化出了一套复杂的抗氧化系统，包括酶促抗氧化系统和非酶促抗氧化系统。酶促抗氧化系统主要由超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）、谷胱甘肽还原酶（GR）等抗氧化酶组成。SOD能够催化O2-歧化为H2O2和O2，POD和CAT则可以将H2O2分解为H2O和O2，GR能够维持谷胱甘肽（GSH）的还原态，保证抗氧化系统的正常运转。非酶促抗氧化系统主要包括一些小分子抗氧化物质，如抗坏血酸（AsA）、谷胱甘肽（GSH）、类胡萝卜素（Car）、生育酚（VE）等。这些抗氧化物质能够直接清除ROS，或者与抗氧化酶协同作用，共同维持细胞内的氧化还原平衡。激素调节在植物应对干旱胁迫中也起着至关重要的作用。植物激素是一类在植物体内含量极低，但对植物的生长发育和逆境响应具有重要调节作用的信号分子。在干旱胁迫下，植物体内多种激素的含量和分布会发生显著变化，其中脱落酸（ABA）是研究最为深入的一种激素。ABA被称为“逆境激素”，在干旱胁迫下，植物根系感知到水分亏缺信号后，会迅速合成ABA并运输到地上部分，ABA通过与受体结合，激活下游一系列信号转导途径，调节气孔关闭，减少水分散失，同时诱导一系列抗旱相关基因的表达，增强植物的抗旱能力。除了ABA外，乙烯（ETH）、生长素（IAA）、细胞分裂素（CTK）、赤霉素（GA）等激素也参与了植物对干旱胁迫的响应过程。它们之间相互作用、相互协调，共同调控植物的生长发育和对干旱胁迫的适应性。例如，ETH在干旱胁迫下能够促进植物根系的生长和发育，增强根系对水分和养分的吸收能力；IAA则可以调节植物的生长和形态建成，使植物在干旱条件下保持适当的生长速度和形态结构。在分子水平上，植物通过基因表达调控来应对干旱胁迫。干旱胁迫会诱导植物体内大量基因的表达发生变化，这些基因被称为干旱响应基因。根据功能的不同，干旱响应基因可以分为两大类：一类是功能基因，它们编码的蛋白质直接参与植物的抗旱生理过程，如前面提到的渗透调节物质合成酶基因、抗氧化酶基因等；另一类是调节基因，它们编码的蛋白质主要参与信号转导和基因表达调控过程，如转录因子、蛋白激酶等。转录因子是一类能够与基因启动子区域的顺式作用元件结合，调控基因转录起始的蛋白质。在植物应对干旱胁迫过程中，多种转录因子家族成员发挥了重要作用，如AP2/ERF、bZIP、MYB、NAC等。这些转录因子通过与下游抗旱相关基因启动子区域的顺式作用元件结合，激活或抑制基因的表达，从而调控植物的抗旱反应。例如，AP2/ERF转录因子家族中的DREB1/CBF亚家族成员，能够特异性地识别并结合到干旱响应元件（DRE）上，激活一系列干旱响应基因的表达，增强植物的抗旱能力。蛋白激酶是一类能够催化蛋白质磷酸化修饰的酶，在植物信号转导过程中起着关键作用。在干旱胁迫下，植物体内多种蛋白激酶被激活，它们通过磷酸化下游底物蛋白，将干旱信号传递下去，最终调节植物的生理生化反应和基因表达。例如，丝裂原活化蛋白激酶（MAPK）级联途径是植物中保守的信号转导途径之一，在干旱胁迫响应中，MAPK级联途径中的多个成员被激活，通过磷酸化修饰转录因子等底物蛋白，调控干旱响应基因的表达。1.5研究目的与内容本研究旨在深入探究白菜型油菜BrPARP1基因在干旱胁迫应答中的功能，为揭示油菜抗旱的分子机制提供理论依据，同时为油菜抗旱品种的选育提供关键基因资源。本研究以白菜型油菜为材料，通过生物信息学分析、基因克隆、表达分析、功能验证等一系列实验手段，对BrPARP1基因在干旱胁迫应答中的功能进行全面深入的研究。具体内容包括：利用生物信息学工具，对BrPARP1基因的核苷酸序列和编码的蛋白质序列进行分析，预测其结构和功能；从白菜型油菜中克隆BrPARP1基因，并构建其过表达载体和RNA干扰载体；通过实时荧光定量PCR（qRT-PCR）技术，分析BrPARP1基因在不同组织以及干旱胁迫不同时间点的表达模式，明确其表达特性；利用农杆菌介导的遗传转化方法，将过表达载体和RNA干扰载体导入白菜型油菜中，获得BrPARP1基因过表达和表达抑制的转基因植株；对转基因植株和野生型植株进行干旱胁迫处理，比较它们在形态指标（如株高、根长、生物量等）、生理指标（如渗透调节物质含量、抗氧化酶活性、光合参数等）以及分子指标（如干旱响应基因的表达水平）等方面的差异，从而验证BrPARP1基因在干旱胁迫应答中的功能；进一步研究BrPARP1基因参与干旱胁迫应答的分子机制，通过酵母双杂交、双分子荧光互补（BiFC）等技术，筛选与BrPARP1蛋白相互作用的蛋白，分析其参与的信号转导途径和基因表达调控网络。二、材料与方法2.1实验材料2.1.1植物材料本研究选用的白菜型油菜品种为“陇油6号”，该品种具有较强的抗逆性，在高寒冷凉地区广泛种植。将白菜型油菜种子播种于装有营养土（蛭石：营养土=1：1）的塑料花盆（直径15cm，高12cm）中，每盆播种10-15粒种子。播种后，浇透水，覆盖保鲜膜以保持湿度，置于光照培养箱中培养。培养条件为：光照强度120-150μmol・m-2・s-1，光照时间16h/d，温度22±2℃，相对湿度60%-70%。待幼苗长出2-3片真叶时，进行间苗，每盆保留5-6株生长健壮的幼苗。在幼苗生长至4-5叶期时，对其进行干旱胁迫处理。采用PEG-6000模拟干旱胁迫，将幼苗根部浸泡在含有不同浓度PEG-6000（0%、5%、10%、15%、20%）的1/2Hoagland营养液中，分别处理0h、6h、12h、24h、48h。处理结束后，迅速采集叶片和根系样品，用液氮速冻后，置于-80℃冰箱中保存，用于后续的生理生化指标测定和基因表达分析。拟南芥（Arabidopsisthaliana）哥伦比亚生态型（Col-0）作为遗传转化的受体材料。拟南芥种子经75%乙醇消毒3-5min，再用无菌水冲洗3-5次后，播种于含有1/2MS培养基（含3%蔗糖、0.7%琼脂，pH5.8）的培养皿中。4℃春化处理3d后，将培养皿转移至光照培养箱中培养，培养条件为：光照强度100-120μmol・m-2・s-1，光照时间16h/d，温度22±1℃，相对湿度60%-70%。待拟南芥幼苗长出4-6片真叶时，移栽至装有营养土（蛭石：营养土=1：1）的塑料花盆（直径10cm，高8cm）中，每盆移栽3-4株，继续培养至拟南芥植株生长至抽薹期，用于遗传转化实验。采用花序浸染法将构建好的过表达载体和RNA干扰载体导入拟南芥中。具体步骤如下：将含有目的载体的农杆菌GV3101菌株接种于含有相应抗生素（利福平50μg/mL、卡那霉素50μg/mL、庆大霉素25μg/mL）的LB液体培养基中，28℃、200r/min振荡培养过夜，至菌液OD600值达到0.8-1.0。将菌液在5000r/min条件下离心10min，收集菌体，用含有5%蔗糖和0.02%SilwetL-77的1/2MS液体培养基重悬菌体，使菌液OD600值调整为0.8-1.0。将拟南芥植株的花序浸入菌液中，轻轻晃动3-5min，确保花序充分接触菌液。侵染后的拟南芥植株用黑色塑料袋包裹，暗培养24h后，置于正常光照条件下培养。待角果自然开裂后，收取T1代种子。将T1代种子用0.1%升汞消毒5-8min，再用无菌水冲洗5-6次后，播种于含有相应抗生素（卡那霉素50μg/mL）的1/2MS培养基上进行筛选。待筛选出的阳性植株生长至4-6片真叶时，移栽至营养土中，继续培养至收获T2代种子。重复筛选过程，直至获得纯合的转基因拟南芥株系。2.1.2菌株与载体实验中使用的大肠杆菌菌株为DH5α和BL21（DE3）。DH5α菌株是一种常用于基因克隆和质粒扩增的菌株，其特点是转化效率高，生长速度快。该菌株缺失了recA1和endA1基因，使得其重组酶活性降低，对质粒的稳定性更好，同时也减少了核酸内切酶对质粒DNA的降解，有利于外源基因的克隆和保存。BL21（DE3）菌株则是一种常用于蛋白表达的菌株，它含有T7RNA聚合酶基因，能够高效表达由T7启动子驱动的外源基因。在该菌株中，T7RNA聚合酶基因由lacUV5启动子控制，通过添加异丙基-β-D-硫代半乳糖苷（IPTG）可以诱导T7RNA聚合酶的表达，从而启动外源基因的转录和翻译，实现蛋白的高效表达。农杆菌菌株选用GV3101，它是一种广泛应用于植物遗传转化的菌株，具有较强的侵染能力和转化效率。该菌株含有pMP90RK辅助质粒，能够提供Vir基因功能，促进T-DNA从农杆菌向植物细胞的转移和整合。在植物遗传转化实验中，将重组表达载体导入农杆菌GV3101中，利用其侵染植物细胞，将载体上的T-DNA区域整合到植物基因组中，从而实现外源基因在植物中的稳定表达。本研究使用的载体包括克隆载体pMD18-T和植物表达载体pCAMBIA1300。pMD18-T载体是一种商业化的克隆载体，其多克隆位点两端分别含有T7和SP6启动子，便于对插入的外源基因进行转录和测序分析。该载体还含有氨苄青霉素抗性基因，用于在大肠杆菌中筛选含有重组质粒的克隆。在基因克隆实验中，将PCR扩增得到的BrPARP1基因片段与pMD18-T载体连接，转化大肠杆菌DH5α感受态细胞，通过蓝白斑筛选和菌落PCR鉴定，获得含有目的基因的重组克隆。pCAMBIA1300载体是一种常用的植物表达载体，它含有CaMV35S启动子、NOS终止子和潮霉素抗性基因。CaMV35S启动子是一种组成型启动子，能够在植物的各种组织和发育阶段中持续启动外源基因的表达。NOS终止子则可以有效地终止转录过程，保证外源基因的稳定表达。潮霉素抗性基因用于在植物转化过程中筛选阳性转基因植株。在构建植物表达载体时，将BrPARP1基因克隆到pCAMBIA1300载体的CaMV35S启动子下游，通过农杆菌介导的转化方法将重组载体导入植物细胞中，实现BrPARP1基因在植物中的过表达。同时，为了构建RNA干扰载体，将BrPARP1基因的部分片段反向插入到pCAMBIA1300载体的内含子两侧，利用RNA干扰技术抑制植物内源BrPARP1基因的表达。2.1.3主要试剂和仪器实验所需的主要试剂包括：PCR试剂（TaqDNA聚合酶、dNTPs、10×PCR缓冲液等）、限制性内切酶（BamHI、SacI等）、T4DNA连接酶、DNAMarker、RNA提取试剂盒（TRIzol试剂）、反转录试剂盒、实时荧光定量PCR试剂盒、质粒提取试剂盒、DNA凝胶回收试剂盒、植物基因组DNA提取试剂盒、蛋白提取试剂（RIPA裂解液、PMSF等）、SDS-PAGE凝胶制备试剂盒、Westernblot相关试剂（一抗、二抗、ECL化学发光试剂等）、PEG-6000、1/2Hoagland营养液、MS培养基、LB培养基、抗生素（氨苄青霉素、卡那霉素、利福平、庆大霉素、潮霉素等）、各种生化试剂（如脯氨酸测定试剂盒、可溶性糖测定试剂盒、超氧化物歧化酶（SOD）活性测定试剂盒、过氧化物酶（POD）活性测定试剂盒等）。主要仪器有：PCR仪（Bio-Rad公司，T100ThermalCycler）、离心机（Eppendorf公司，5424R型）、凝胶成像系统（Bio-Rad公司，ChemiDocMPImagingSystem）、实时荧光定量PCR仪（Roche公司，LightCycler480II）、核酸蛋白分析仪（ThermoScientific公司，Nanodrop2000）、恒温振荡培养箱（NewBrunswick公司，Innova44R）、光照培养箱（上海一恒科学仪器有限公司，MGC-450HP）、超净工作台（苏州净化设备有限公司，SW-CJ-2FD）、电泳仪（Bio-Rad公司，PowerPacBasic）、转膜仪（Bio-Rad公司，Trans-BlotTurboTransferSystem）、酶标仪（ThermoScientific公司，MultiskanFC）等。2.2实验方法2.2.1白菜型油菜BrPARP1基因的克隆采用TRIzol试剂法提取干旱胁迫处理48h后白菜型油菜叶片的总RNA。具体操作步骤如下：取约100mg叶片组织，迅速放入液氮中研磨成粉末状，将粉末转移至含有1mLTRIzol试剂的无RNase离心管中，剧烈振荡混匀，室温静置5min，使组织充分裂解。加入200μL***，剧烈振荡15s，室温静置3min。4℃、12000r/min离心15min，此时溶液分为三层，上层为无色透明的水相，中间为白色的蛋白质层，下层为红色的有机相，将上层水相转移至新的无RNase离心管中。加入等体积的异丙醇，颠倒混匀，室温静置10min。4℃、12000r/min离心10min，可见管底出现白色RNA沉淀，弃上清。用75%乙醇（用DEPC水配制）洗涤RNA沉淀2次，每次加入1mL75%乙醇，颠倒混匀，4℃、7500r/min离心5min，弃上清。将RNA沉淀在室温下晾干5-10min，加入适量的DEPC水溶解RNA，用核酸蛋白分析仪测定RNA的浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0。将提取的RNA保存于-80℃冰箱备用。利用反转录试剂盒将提取的总RNA反转录成cDNA。以Oligo(dT)18为引物，在反转录酶的作用下，将RNA反转录为cDNA。反应体系为20μL，包括5×PrimeScriptBuffer4μL、PrimeScriptRTEnzymeMixI1μL、Oligo(dT)18Primer1μL、总RNA1μg，用RNase-freedH2O补足至20μL。反应条件为：37℃15min，85℃5s，4℃保存。反转录得到的cDNA保存于-20℃冰箱备用。根据白菜型油菜基因组数据库中BrPARP1基因的序列信息，利用PrimerPremier6.0软件设计特异性引物。上游引物（5'-3'）：ATGGCTTCTTCTTCGCTGCT；下游引物（5'-3'）：TCACTGCTGCTGCTGCTGCT。引物由生工生物工程（上海）股份有限公司合成。以反转录得到的cDNA为模板，进行PCR扩增。PCR反应体系为25μL，包括2×TaqPCRMasterMix12.5μL、上游引物（10μmol/L）1μL、下游引物（10μmol/L）1μL、cDNA模板1μL，用ddH2O补足至25μL。反应条件为：95℃预变性5min；95℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；72℃终延伸10min。PCR扩增产物经1%琼脂糖凝胶电泳检测，在凝胶成像系统下观察并拍照，切下目的条带，用DNA凝胶回收试剂盒回收纯化。将回收的PCR产物与pMD18-T载体连接，连接体系为10μL，包括pMD18-TVector0.5μL、回收的PCR产物4.5μL、SolutionI5μL。16℃连接过夜。将连接产物转化大肠杆菌DH5α感受态细胞。具体步骤为：取10μL连接产物加入到100μLDH5α感受态细胞中，轻轻混匀，冰浴30min；42℃热激90s，迅速放回冰浴中冷却2min；加入900μL不含抗生素的LB液体培养基，37℃、200r/min振荡培养1h；将菌液在5000r/min条件下离心5min，弃上清，留100μL菌液重悬菌体，将重悬后的菌液涂布于含有氨苄青霉素（100μg/mL）的LB固体培养基平板上，37℃倒置培养过夜。次日，挑取平板上的单菌落接种于含有氨苄青霉素（100μg/mL）的LB液体培养基中，37℃、200r/min振荡培养过夜。提取质粒，用BamHI和SacI限制性内切酶进行双酶切鉴定，酶切体系为20μL，包括质粒2μL、10×Buffer2μL、BamHI1μL、SacI1μL，用ddH2O补足至20μL。37℃酶切2-3h，酶切产物经1%琼脂糖凝胶电泳检测。将酶切鉴定正确的重组质粒送生工生物工程（上海）股份有限公司进行测序。2.2.2BrPARP1基因的生物信息学分析利用NCBI（NationalCenterforBiotechnologyInformation）网站的BLAST工具，对克隆得到的BrPARP1基因序列进行同源性分析，与数据库中已有的基因序列进行比对，确定其与其他物种PARP1基因的相似性。同时，使用DNAMAN软件对BrPARP1基因及其编码的氨基酸序列进行分析，计算其分子量、等电点、氨基酸组成等基本参数。通过在线工具ExPASy（ExpertProteinAnalysisSystem）中的ProtParam程序，进一步精确分析BrPARP1蛋白的理化性质。利用在线工具SOPMA（Self-OptimizedPredictionMethodwithAlignment）预测BrPARP1蛋白的二级结构，包括α-螺旋、β-折叠、β-转角和无规则卷曲等结构元件。使用SWISS-MODEL在线服务器，根据已知的蛋白质晶体结构，对BrPARP1蛋白进行三维结构建模，预测其三级结构。运用NCBI的ConservedDomainDatabase（CDD）数据库，分析BrPARP1蛋白的保守结构域，确定其功能结构域的位置和类型。通过在线工具InterProScan，对BrPARP1蛋白进行功能注释，预测其可能参与的生物学过程和分子功能。利用MEGA（MolecularEvolutionaryGeneticsAnalysis）软件，采用邻接法（Neighbor-Joiningmethod）构建BrPARP1蛋白与其他物种PARP1蛋白的系统进化树，分析其进化关系。2.2.3干旱胁迫处理及BrPARP1基因表达分析选取生长状况一致的4-5叶期白菜型油菜幼苗，随机分为对照组和干旱胁迫组。对照组正常浇水，干旱胁迫组采用PEG-6000模拟干旱胁迫，将幼苗根部浸泡在含有20%PEG-6000的1/2Hoagland营养液中。分别在处理0h、6h、12h、24h、48h后，采集叶片和根系样品，用液氮速冻后，置于-80℃冰箱中保存，用于后续的基因表达分析。采用TRIzol试剂法提取不同处理时间点白菜型油菜叶片和根系的总RNA，具体步骤同2.2.1。利用反转录试剂盒将总RNA反转录成cDNA，反应体系和条件也同2.2.1。以cDNA为模板，进行实时荧光定量PCR（qRT-PCR）分析。以白菜型油菜的Actin基因作为内参基因，引物序列如下：Actin上游引物（5'-3'）：GACCTGACAGACTACCTCATG；Actin下游引物（5'-3'）：AGAGGCATACAGGGACAGCAA。BrPARP1基因的qRT-PCR引物序列为：上游引物（5'-3'）：CCGCTGCTGCTGCTGCTGCT；下游引物（5'-3'）：GCTGCTGCTGCTGCTGCTGCT。引物由生工生物工程（上海）股份有限公司合成。qRT-PCR反应体系为20μL，包括2×SYBRGreenMasterMix10μL、上游引物（10μmol/L）0.8μL、下游引物（10μmol/L）0.8μL、cDNA模板1μL，用ddH2O补足至20μL。反应在RocheLightCycler480II实时荧光定量PCR仪上进行，反应条件为：95℃预变性5min；95℃变性10s，60℃退火30s，72℃延伸30s，共40个循环。每个样品设置3个生物学重复和3个技术重复。采用2-ΔΔCT法计算BrPARP1基因的相对表达量，分析其在干旱胁迫下不同时间点的表达变化情况。2.2.4BrPARP1基因功能验证以pCAMBIA1300为基础载体，构建BrPARP1基因的过表达载体。用BamHI和SacI限制性内切酶分别对pCAMBIA1300载体和含有BrPARP1基因的pMD18-T重组质粒进行双酶切，酶切体系和条件同2.2.1。酶切产物经1%琼脂糖凝胶电泳检测后，用DNA凝胶回收试剂盒回收目的片段。将回收的BrPARP1基因片段与线性化的pCAMBIA1300载体片段用T4DNA连接酶连接，连接体系为10μL，包括线性化的pCAMBIA1300载体1μL、回收的BrPARP1基因片段3μL、T4DNALigase1μL、10×T4DNALigaseBuffer1μL，用ddH2O补足至10μL。16℃连接过夜。将连接产物转化大肠杆菌DH5α感受态细胞，转化步骤同2.2.1。挑取单菌落进行菌落PCR鉴定和质粒双酶切鉴定，鉴定正确的重组质粒命名为pCAMBIA1300-BrPARP1。采用冻融法将构建好的过表达载体pCAMBIA1300-BrPARP1导入农杆菌GV3101感受态细胞中。具体步骤为：取1μg重组质粒加入到100μLGV3101感受态细胞中，轻轻混匀，冰浴30min；液氮速冻5min，37℃水浴5min，迅速放回冰浴中冷却2min；加入900μL不含抗生素的LB液体培养基，28℃、200r/min振荡培养2-3h；将菌液在5000r/min条件下离心5min，弃上清，留100μL菌液重悬菌体，将重悬后的菌液涂布于含有利福平（50μg/mL）、卡那霉素（50μg/mL）和庆大霉素（25μg/mL）的LB固体培养基平板上，28℃倒置培养2-3d。挑取平板上的单菌落接种于含有相应抗生素的LB液体培养基中，28℃、200r/min振荡培养过夜，提取质粒进行PCR鉴定和双酶切鉴定，鉴定正确的农杆菌菌株用于后续的遗传转化实验。利用花序浸染法将含有过表达载体pCAMBIA1300-BrPARP1的农杆菌GV3101转化拟南芥，具体步骤同2.1.1。将收获的T1代种子用0.1%升汞消毒5-8min，再用无菌水冲洗5-6次后，播种于含有卡那霉素（50μg/mL）的1/2MS培养基上进行筛选。待筛选出的阳性植株生长至4-6片真叶时，移栽至营养土中，继续培养至收获T2代种子。重复筛选过程，直至获得纯合的转基因拟南芥株系。对野生型拟南芥（WT）和纯合的转基因拟南芥株系（OE-1、OE-2、OE-3）进行干旱胁迫处理。选取生长状况一致的4-5周龄拟南芥植株，停止浇水，进行干旱胁迫处理10d，然后恢复浇水3d，观察并记录植株的生长状况和表型变化。同时，测定干旱胁迫处理前后植株的鲜重、干重、根长、株高等形态指标，分析转基因植株的抗旱性。2.2.5生理生化指标测定在干旱胁迫处理0d、5d、10d时，分别采集野生型拟南芥和转基因拟南芥植株的叶片，测定其相对含水量（RWC）、丙二醛（MDA）含量、脯氨酸含量、可溶性糖含量、超氧化物歧化酶（SOD）活性、过氧化物酶（POD）活性和过氧化氢酶（CAT）活性等生理生化指标。相对含水量的测定采用称重法。取新鲜叶片，用吸水纸吸干表面水分，称取鲜重（FW），然后将叶片浸泡在蒸馏水中，4℃黑暗条件下浸泡24h，使其充分吸水饱和，取出吸干表面水分，称取饱和鲜重（TW），最后将叶片置于80℃烘箱中烘干至恒重，称取干重（DW）。相对含水量计算公式为：RWC（%）=（FW-DW）/（TW-DW）×100。丙二醛含量的测定采用硫代巴比妥酸（TBA）法。取0.5g叶片，加入5mL5%三***（TCA）溶液，研磨成匀浆，4℃、10000r/min离心10min，取上清液2mL，加入2mL0.6%TBA溶液（用5%TCA配制），混匀后于沸水浴中加热15min，迅速冷却后，4℃、10000r/min离心10min，取上清液，用分光光度计在450nm、532nm和600nm波长下测定吸光值。根据公式计算MDA含量：MDA含量（μmol/gFW）=6.45×（A532-A600）-0.56×A450。脯氨酸含量的测定采用酸性茚三法。取0.5g叶片，加入5mL3%磺基水杨酸溶液，研磨成匀浆，沸水浴中提取10min，冷却后，4℃、10000r/min离心10min，取上清液1mL，加入2mL冰乙酸和3mL酸性茚三试剂，混匀后于沸水浴中加热30min，迅速冷却后，加入5mL甲苯，振荡萃取，取上层甲苯相，用分光光度计在520nm波长下测定吸光值。根据脯氨酸标准曲线计算脯氨酸含量。可溶性糖含量的测定采用蒽比色法。取0.5g叶片，加入5mL80%乙醇，研磨成匀浆，80℃水浴中提取30min，冷却后，4℃、10000r/min离心10min，取上清液1mL，加入1mL蒸馏水和5mL蒽试剂，混匀后于沸水浴中加热10min，迅速冷却后，用分光光度计在620nm波长下测定吸光值。根据可溶性糖标准曲线计算可溶性糖含量。超氧化物歧化酶活性的测定采用氮蓝四唑（NBT）光化还原法。取0.5g叶片，加入5mL50mmol/L磷酸缓冲液（pH7.8，含1%聚乙烯吡咯烷***），研磨成匀浆，4℃、10000r/min离心20min，取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.8）1.5mL、130mmol/L甲硫氨酸溶液0.3mL、750μmol/LNBT溶液0.3mL、100μmol/LEDTA-Na2溶液0.3mL、20μmol/L核黄素溶液0.3mL和酶液0.3mL。混匀后，将反应管置于光照下反应15-20min，然后用分光光度计在560nm波长下测定吸光值。以抑制NBT光化还原50%的酶量为一个酶活性单位（U），计算SOD活性。过氧化物酶活性的测定采用愈创木酚法。取0.5g叶片，加入5mL50mmol/L磷酸缓冲液（pH7.0，含1%聚乙烯吡咯烷***），研磨成匀浆，4℃、10000r/min离心20min，取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.0）2.7mL、2%愈创木酚溶液0.1mL、0.3%H2O2溶液0.1mL和酶液0.1mL。混匀后，于37℃水浴中反应5min，然后用分光光度计在470nm波长下测定吸光值。以每分钟吸光值变化0.01为一个酶活性单位（U），计算POD活性。过氧化氢酶活性的测定采用紫外吸收法。取0.5g叶片，加入5mL50mmol/L磷酸缓冲液（pH7.0，含1%聚乙烯吡咯烷***），研磨成匀浆，4℃、10000r/min离心20min，取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.0）2.9mL、0.1mol三、结果与分析3.1白菜型油菜BrPARP1基因的克隆与序列分析以干旱胁迫处理48h后白菜型油菜叶片的cDNA为模板，利用特异性引物进行PCR扩增，成功克隆得到BrPARP1基因。PCR扩增产物经1%琼脂糖凝胶电泳检测，结果显示在约1500bp处出现一条清晰的条带，与预期的BrPARP1基因大小相符（图1）。将该条带切下，进行回收纯化，并与pMD18-T载体连接，转化大肠杆菌DH5α感受态细胞。挑取单菌落进行菌落PCR鉴定和质粒双酶切鉴定，结果表明重组质粒中含有目的基因片段，且酶切鉴定条带大小与预期一致（图2）。将鉴定正确的重组质粒送测序，测序结果显示，克隆得到的BrPARP1基因全长为1485bp，编码494个氨基酸。通过NCBI的BLAST工具对BrPARP1基因序列进行同源性分析，结果表明，BrPARP1基因与其他物种的PARP1基因具有较高的相似性。其中，与拟南芥（Arabidopsisthaliana）PARP1基因（AT1G77490）的核苷酸序列相似度达到82%，氨基酸序列相似度达到85%；与甘蓝（Brassicaoleracea）PARP1基因的核苷酸序列相似度为88%，氨基酸序列相似度为90%。这表明BrPARP1基因在进化过程中具有较高的保守性。利用DNAMAN软件对BrPARP1基因编码的氨基酸序列进行分析，结果显示，BrPARP1蛋白的分子量约为55.8kDa，等电点为8.95。氨基酸组成分析表明，BrPARP1蛋白中含量较高的氨基酸为亮氨酸（Leu）、甘氨酸（Gly）、丝氨酸（Ser）和丙氨酸（Ala），分别占氨基酸总数的10.9%、9.9%、9.1%和8.3%。通过在线工具ExPASy中的ProtParam程序进一步分析BrPARP1蛋白的理化性质，结果显示，BrPARP1蛋白的不稳定系数为43.67，属于不稳定蛋白；脂肪系数为95.55，总平均亲水性为-0.301，表明该蛋白具有一定的亲水性。利用在线工具SOPMA预测BrPARP1蛋白的二级结构，结果显示，BrPARP1蛋白主要由α-螺旋（32.39%）、β-折叠（18.22%）、β-转角（6.68%）和无规则卷曲（42.71%）组成（图3）。其中，α-螺旋和无规则卷曲在蛋白结构中所占比例较高，它们共同构成了BrPARP1蛋白的基本骨架，为蛋白的功能发挥提供了结构基础。α-螺旋结构具有高度的稳定性，能够维持蛋白的特定构象，参与蛋白与其他分子的相互作用；无规则卷曲则赋予了蛋白一定的柔性，使其能够在不同的环境条件下发生构象变化，从而更好地执行其生物学功能。β-折叠和β-转角结构相对较少，但它们在维持蛋白的三维结构和功能方面也起着重要的作用。β-折叠结构能够增强蛋白的稳定性，促进蛋白之间的相互作用；β-转角则可以改变蛋白的走向，使蛋白形成特定的空间结构，有利于其与底物或其他蛋白的结合。使用SWISS-MODEL在线服务器对BrPARP1蛋白进行三维结构建模，结果显示，BrPARP1蛋白具有典型的PARP1蛋白结构特征，包含N端的DNA结合结构域、中部的催化结构域和C端的调节结构域（图4）。DNA结合结构域由多个α-螺旋和β-折叠组成，形成了一个能够特异性识别和结合DNA的凹槽结构；催化结构域则含有保守的催化活性位点，能够催化聚[ADP-核糖]（PAR）的合成；C端的调节结构域较为灵活，可能通过与其他蛋白质相互作用，调节BrPARP1蛋白的活性和功能。运用NCBI的ConservedDomainDatabase（CDD）数据库分析BrPARP1蛋白的保守结构域，结果表明，BrPARP1蛋白含有PARP家族特有的PARP-catalytic结构域（图5），该结构域位于氨基酸序列的第200-450位，包含了催化PAR合成所需的关键氨基酸残基，是BrPARP1蛋白发挥催化功能的核心区域。此外，BrPARP1蛋白还含有一个Zn-finger结构域，位于N端的第30-80位氨基酸，该结构域可能参与DNA的结合和识别过程，对BrPARP1蛋白的功能具有重要的调节作用。通过在线工具InterProScan对BrPARP1蛋白进行功能注释，结果预测BrPARP1蛋白主要参与DNA损伤修复、细胞凋亡、转录调控等生物学过程，具有聚[ADP-核糖]聚合酶活性、DNA结合活性等分子功能。这与其他物种中PARP1蛋白的功能研究结果一致，进一步表明BrPARP1基因在植物生长发育和逆境响应过程中可能发挥着重要作用。利用MEGA软件，采用邻接法（Neighbor-Joiningmethod）构建BrPARP1蛋白与其他物种PARP1蛋白的系统进化树，结果显示，BrPARP1蛋白与十字花科植物的PARP1蛋白聚为一支，其中与甘蓝的PARP1蛋白亲缘关系最近，这与它们在分类学上的亲缘关系相符（图6）。在进化过程中，BrPARP1蛋白与其他物种的PARP1蛋白可能经历了共同的祖先基因分化，随着物种的进化，逐渐形成了各自独特的序列和功能特征，但在关键结构域和生物学功能上仍保持着较高的保守性。这种进化关系的分析有助于深入理解BrPARP1基因的起源和进化历程，为进一步研究其功能和应用提供了重要的参考依据。3.2BrPARP1基因的生物信息学分析利用NCBI网站的BLAST工具对克隆得到的BrPARP1基因序列进行同源性分析，结果显示其与其他物种的PARP1基因展现出较高的相似性。其中，与拟南芥PARP1基因（AT1G77490）在核苷酸序列上的相似度达82%，氨基酸序列相似度达85%；和甘蓝PARP1基因的核苷酸序列相似度为88%，氨基酸序列相似度为90%。这种高度的相似性表明BrPARP1基因在漫长的进化历程中具有较强的保守性，其功能可能在不同物种间存在一定的共性。运用DNAMAN软件对BrPARP1基因编码的氨基酸序列进行分析，得出BrPARP1蛋白的分子量约为55.8kDa，等电点为8.95。进一步通过ExPASy中的ProtParam程序分析，发现BrPARP1蛋白的不稳定系数为43.67，按照蛋白质稳定性的判断标准，该数值表明其属于不稳定蛋白；脂肪系数为95.55，总平均亲水性为-0.301，说明该蛋白具有一定的亲水性。这些理化性质与蛋白在细胞内的定位、折叠以及与其他分子的相互作用密切相关，为深入理解BrPARP1蛋白的功能提供了基础信息。通过在线工具SOPMA预测BrPARP1蛋白的二级结构，结果显示该蛋白主要由α-螺旋（32.39%）、β-折叠（18.22%）、β-转角（6.68%）和无规则卷曲（42.71%）组成。α-螺旋和无规则卷曲在蛋白结构中所占比例较高，α-螺旋具有高度的稳定性，能够维持蛋白的特定构象，参与蛋白与其他分子的相互作用；无规则卷曲赋予了蛋白一定的柔性，使其能在不同环境条件下发生构象变化，从而更好地执行生物学功能。β-折叠和β-转角结构虽相对较少，但在维持蛋白的三维结构和功能方面也起着不可或缺的作用。β-折叠结构能够增强蛋白的稳定性，促进蛋白之间的相互作用；β-转角则可以改变蛋白的走向，使蛋白形成特定的空间结构，有利于其与底物或其他蛋白的结合。使用SWISS-MODEL在线服务器对BrPARP1蛋白进行三维结构建模，结果清晰地显示BrPARP1蛋白具有典型的PARP1蛋白结构特征，包含N端的DNA结合结构域、中部的催化结构域和C端的调节结构域。DNA结合结构域由多个α-螺旋和β-折叠组成，形成了一个能够特异性识别和结合DNA的凹槽结构，这一结构特点使得BrPARP1蛋白能够精准地定位到受损的DNA区域，启动后续的修复或调控过程；催化结构域含有保守的催化活性位点，能够催化聚[ADP-核糖]（PAR）的合成，PAR的合成在DNA损伤修复、基因表达调控等生物学过程中发挥着关键作用；C端的调节结构域较为灵活，可能通过与其他蛋白质相互作用，调节BrPARP1蛋白的活性和功能，这种调节作用有助于细胞根据不同的生理状态和环境信号，对BrPARP1蛋白的功能进行精准调控。运用NCBI的ConservedDomainDatabase（CDD）数据库分析BrPARP1蛋白的保守结构域，结果表明BrPARP1蛋白含有PARP家族特有的PARP-catalytic结构域，该结构域位于氨基酸序列的第200-450位，包含了催化PAR合成所需的关键氨基酸残基，是BrPARP1蛋白发挥催化功能的核心区域。此外，BrPARP1蛋白还含有一个Zn-finger结构域，位于N端的第30-80位氨基酸，该结构域可能参与DNA的结合和识别过程，对BrPARP1蛋白的功能具有重要的调节作用。Zn-finger结构域通过与DNA的特定序列相互作用，能够增强BrPARP1蛋白与DNA的结合亲和力和特异性，从而提高其在DNA损伤修复和基因表达调控等过程中的效率。通过在线工具InterProScan对BrPARP1蛋白进行功能注释，结果预测BrPARP1蛋白主要参与DNA损伤修复、细胞凋亡、转录调控等生物学过程，具有聚[ADP-核糖]聚合酶活性、DNA结合活性等分子功能。这与其他物种中PARP1蛋白的功能研究结果一致，进一步表明BrPARP1基因在植物生长发育和逆境响应过程中可能发挥着重要作用。在DNA损伤修复过程中，BrPARP1蛋白能够迅速响应DNA损伤信号，通过催化PAR的合成，招募相关的修复蛋白到损伤位点，促进DNA的修复，维持基因组的稳定性；在细胞凋亡过程中，BrPARP1蛋白可能通过调节相关基因的表达，参与细胞凋亡的调控，决定细胞的命运；在转录调控过程中，BrPARP1蛋白可以与转录因子等其他蛋白质相互作用，影响基因的转录起始和延伸，从而调控基因的表达水平。利用MEGA软件，采用邻接法构建BrPARP1蛋白与其他物种PARP1蛋白的系统进化树，结果显示BrPARP1蛋白与十字花科植物的PARP1蛋白聚为一支，其中与甘蓝的PARP1蛋白亲缘关系最近。这与它们在分类学上的亲缘关系相符，表明在进化过程中，BrPARP1蛋白与其他物种的PARP1蛋白可能经历了共同的祖先基因分化，随着物种的进化，逐渐形成了各自独特的序列和功能特征，但在关键结构域和生物学功能上仍保持着较高的保守性。这种进化关系的分析有助于深入理解BrPARP1基因的起源和进化历程，为进一步研究其功能和应用提供了重要的参考依据。通过对进化树的分析，可以推测BrPARP1基因在不同物种中的功能演化路径，以及与其他物种PARP1基因之间的亲缘关系，从而为利用其他物种的研究成果来深入探讨BrPARP1基因的功能提供了可能。3.3干旱胁迫下BrPARP1基因的表达模式利用实时荧光定量PCR（qRT-PCR）技术，对干旱胁迫下白菜型油菜不同组织（叶片和根系）中BrPARP1基因的表达模式进行了分析。以Actin基因作为内参基因，采用2-ΔΔCT法计算BrPARP1基因的相对表达量。结果如图7所示，在正常生长条件下（0h），BrPARP1基因在白菜型油菜的叶片和根系中均有表达，但表达水平相对较低。随着干旱胁迫时间的延长，BrPARP1基因在叶片和根系中的表达量均呈现出先上升后下降的趋势。在叶片中，干旱胁迫6h时，BrPARP1基因的表达量开始显著上调，与0h相比，增加了约2.5倍（P<0.05）；在胁迫12h时，表达量达到峰值，为0h的4.2倍（P<0.01）；随后，表达量逐渐下降，在胁迫48h时，虽仍高于0h的表达水平，但与峰值相比已明显降低（图7A）。在根系中，BrPARP1基因的表达变化趋势与叶片相似，但响应时间略有延迟。干旱胁迫12h时，表达量开始显著增加，为0h的2.1倍（P<0.05）；在胁迫24h时，表达量达到最高，是0h的3.5倍（P<0.01）；之后表达量逐渐回落，至胁迫48h时，表达量仍高于初始水平（图7B）。这种表达模式的变化表明，BrPARP1基因能够对干旱胁迫作出快速响应，其表达量的上调可能与白菜型油菜对干旱胁迫的适应和防御机制密切相关。在干旱胁迫初期，BrPARP1基因表达量的增加，可能是植物为了应对干旱胁迫带来的损伤，启动了一系列的保护机制，通过调节相关基因的表达，参与渗透调节、抗氧化防御等生理过程，以增强植物的抗旱能力。随着胁迫时间的进一步延长，植物可能由于受到的伤害逐渐加重，自身的调节能力逐渐下降，导致BrPARP1基因的表达量开始降低。此外，BrPARP1基因在叶片和根系中的表达变化存在一定的差异，这可能与叶片和根系在植物生长发育和逆境响应中的不同功能有关。叶片是植物进行光合作用和蒸腾作用的主要器官，对干旱胁迫更为敏感，因此BrPARP1基因在叶片中的响应更为迅速；而根系作为植物吸收水分和养分的重要器官，在干旱胁迫下，需要维持自身的生长和功能，以保证植物对水分和养分的吸收，因此BrPARP1基因在根系中的表达变化相对较为滞后，但持续时间可能更长。3.4BrPARP1基因功能验证结果3.4.1转基因拟南芥的获得与鉴定通过农杆菌介导的花序浸染法，将构建好的BrPARP1基因过表达载体pCAMBIA1300-BrPARP1导入拟南芥中，经过多代筛选，成功获得了转基因拟南芥株系。为了鉴定转基因拟南芥，提取T3代转基因拟南芥和野生型拟南芥的基因组DNA，以其为模板，利用特异性引物进行PCR扩增。扩增产物经1%琼脂糖凝胶电泳检测，结果如图8所示，在转基因拟南芥株系（OE-1、OE-2、OE-3）中均扩增出了与目的基因大小相符的条带（约1500bp），而野生型拟南芥中未扩增出该条带，表明BrPARP1基因已成功整合到拟南芥基因组中。进一步对转基因拟南芥株系进行qRT-PCR分析，检测BrPARP1基因在转录水平的表达情况。以野生型拟南芥为对照，Actin基因作为内参基因，采用2-ΔΔCT法计算BrPARP1基因的相对表达量。结果如图9所示，在转基因拟南芥株系中，BrPARP1基因的表达量均显著高于野生型拟南芥（P<0.01）。其中，OE-2株系的表达量最高，约为野生型的8.5倍；OE-1和OE-3株系的表达量分别为野生型的6.8倍和7.2倍。这表明BrPARP1基因在转基因拟南芥中能够高效转录，过表达载体构建成功，且不同株系之间的表达水平存在一定差异。这些转基因拟南芥株系为后续研究BrPARP1基因的功能提供了良好的材料。3.4.2转基因拟南芥的抗旱性分析选取生长状况一致的4-5周龄野生型拟南芥（WT）和转基因拟南芥株系（OE-1、OE-2、OE-3），停止浇水进行干旱胁迫处理10d，然后恢复浇水3d，观察植株的生长状况和表型变化。结果如图10所示，在干旱胁迫处理期间，野生型拟南芥植株的叶片逐渐变黄、萎蔫，生长受到明显抑制；而转基因拟南芥株系的叶片虽然也出现了一定程度的萎蔫，但相对野生型而言，其叶片的绿色程度更高，萎蔫程度较轻，生长受抑制的程度较小。恢复浇水3d后，野生型拟南芥仅有少数植株能够恢复生长，大部分植株死亡；而转基因拟南芥株系的恢复能力较强，大部分植株能够恢复正常生长，叶片重新展开，颜色变绿。对干旱胁迫处理前后植株的鲜重、干重、根长、株高等形态指标进行测定，结果如表1所示。干旱胁迫处理后，野生型和转基因拟南芥植株的鲜重、干重、根长、株高均显著下降（P<0.05）。但转基因拟南芥株系的各项指标下降幅度明显小于野生型拟南芥。在鲜重方面，野生型拟南芥下降了约65%，而OE-1、OE-2、OE-3株系分别下降了48%、45%、50%；在干重方面，野生型下降了60%，转基因株系OE-1、OE-2、OE-3分别下降了45%、42%、47%；根长方面，野生型下降了55%，转基因株系OE-1、OE-2、OE-3分别下降了38%、35%、40%；株高方面，野生型下降了50%，转基因株系OE-1、OE-2、OE-3分别下降了30%、28%、32%。这些结果表明，过表达BrPARP1基因能够显著提高拟南芥的抗旱性，使植株在干旱胁迫下保持较好的生长状态，减少干旱对植株生长发育的抑制作用。3.4.3生理生化指标分析在干旱胁迫处理0d、5d、10d时，分别采集野生型拟南芥和转基因拟南芥植株的叶片，测定其相对含水量（RWC）、丙二醛（MDA）含量、脯氨酸含量、可溶性糖含量、超氧化物歧化酶（SOD）活性、过氧化物酶（POD）活性和过氧化氢酶（CAT）活性等生理生化指标，以探究BrPARP1基因提高拟南芥抗旱性的生理机制。相对含水量是反映植物水分状况的重要指标之一。如图11A所示，在正常生长条件下（0d），野生型和转基因拟南芥植株的叶片相对含水量无显著差异。随着干旱胁迫时间的延长，两者的相对含水量均逐渐下降，但转基因拟南芥株系的相对含水量始终显著高于野生型拟南芥（P<0.05）。在干旱胁迫10d时，野生型拟南芥叶片的相对含水量降至40%左右，而OE-1、OE-2、OE-3株系的相对含水量分别为55%、58%、56%，表明过表达BrPARP1基因能够有效维持干旱胁迫下植株的水分平衡，减少水分散失。丙二醛是膜脂过氧化的产物，其含量高低反映了植物细胞膜受损伤的程度。图11B显示，干旱胁迫处理后，野生型和转基因拟南芥植株叶片的MDA含量均显著增加（P<0.05），但转基因拟南芥株系的MDA含量显著低于野生型拟南芥。在干旱胁迫10d时，野生型拟南芥叶片的MDA含量达到12nmol/gFW以上，而OE-1、OE-2、OE-3株系的MDA含量分别为8.5nmol/gFW、8.0nmol/gFW、8.3nmol/gFW，说明过表达BrPARP1基因能够减轻干旱胁迫对细胞膜的损伤，提高细胞膜的稳定性。脯氨酸和可溶性糖是植物体内重要的渗透调节物质，在干旱胁迫下，它们的积累有助于降低细胞的渗透势，维持细胞的膨压。图11C和图11D表明，随着干旱胁迫时间的延长，野生型和转基因拟南芥植株叶片的脯氨酸和可溶性糖含量均逐渐增加。在干旱胁迫10d时，转基因拟南芥株系的脯氨酸和可溶性糖含量显著高于野生型拟南芥（P<0.01）。其中，OE-2株系的脯氨酸含量达到450μg/gFW以上，可溶性糖含量达到18mg/gFW以上，分别是野生型的2.5倍和1.8倍，说明过表达BrPARP1基因能够促进渗透调节物质的积累，增强植物的渗透调节能力，从而提高植物的抗旱性。超氧化物歧化酶、过氧化物酶和过氧化氢酶是植物抗氧化系统的关键酶，它们能够清除细胞内过多的活性氧，减轻氧化损伤。图11E、图11F和图11G显示，干旱胁迫处理后，野生型和转基因拟南芥植株叶片的SOD、POD和CAT活性均显著增强（P<0.05），且转基因拟南芥株系的酶活性显著高于野生型拟南芥。在干旱胁迫10d时，OE-2株系的SOD活性达到600U/gFW以上，POD活性达到1500U/gFW以上，CAT活性达到350U/gFW以上，分别是野生型的1.8倍、2.0倍、1.6倍，表明过表达BrPARP1基因能够增强植物抗氧化酶的活性，提高植物的抗氧化能力，有效清除活性氧，保护植物细胞免受氧化损伤。四、讨论4.1BrPARP1基因的结构与功能预测本研究通过生物信息学分析，对白菜型油菜BrPARP1基因的结构和功能进行了深入探究。从基因序列分析结果来看，BrPARP1基因全长1485bp，编码494个氨基酸，与其他物种的PARP1基因具有较高的相似性，这表明PARP1基因在进化过程中具有高度的保守性。这种保守性暗示了该基因在植物生长发育和逆境响应等生物学过程中可能发挥着重要且保守的功能。在长期的进化历程中，尽管物种不断演化，但PARP1基因始终保持着相对稳定的序列特征，这充分说明了其功能对于植物生存和繁衍的重要性。对BrPARP1蛋白的理化性质分析显示，其分子量约为55.8kDa，等电点为8.95，属于不稳定的亲水性蛋白。蛋白的稳定性和溶解性等理化性质与其在细胞内的定位、折叠以及与其他分子的相互作用密切相关。不稳定的蛋白结构可能使BrPARP1更容易发生构象变化，从而更好地响应外界环境信号和执行生物学功能。而亲水性特征则表明它可能在细胞内的水环境中发挥作用，参与各种生理生化反应。在蛋白结构方面，BrPARP1蛋白具有典型的PARP1蛋白结构特征，包含N端的DNA结合结构域、中部的催化结构域和C端的调节结构域。DNA结合结构域能够特异性地识别并结合受损的DNA分子，这是PARP1发挥功能的关键起始步骤。在干旱胁迫等逆境条件下，植物细胞内的DNA容易受到损伤，BrPARP1的DNA结合结构域可以迅速感知并结合到受损的DNA位点，启动后续的修复过程。中部的催化结构域含有保守的催化活性位点，能够催化聚[ADP-核糖]（PAR）的合成。PAR在DNA损伤修复、基因表达调控等生物学过程中发挥着关键作用。在DNA损伤修复过程中，PAR可以作为一种信号分子，招募相关的修复蛋白到损伤位点，形成修复复合物，促进DNA的修复。同时，PAR还可以通过与染色质结合，改变染色质的结构和功能，影响基因的表达。C端的调节结构域较为灵活，可能通过与其他蛋白质相互作用，调节BrPARP1蛋白的活性和功能。这种调节作用有助于细胞根据不同的生理状态和环境信号，对BrPARP1蛋白的功能进行精准调控。例如，在干旱胁迫下，C端调节结构域可能与一些胁迫响应蛋白相互作用，增强BrPARP1蛋白的活性，从而更好地应对干旱胁迫带来的损伤。此外，BrPARP1蛋白还含有PARP家族特有的PARP-catalytic结构域以及一个Zn-finger结构域。PARP-catalytic结构域是BrPARP1蛋白发挥催化功能的核心区域，包含了催化PAR合成所需的关键氨基酸残基。这些氨基酸残基通过特定的空间排列，形成了催化活性中心，确保了催化反应的高效进行。Zn-finger结构域可能参与DNA的结合和识别过程，对BrPARP1蛋白的功能具有重要的调节作用。Zn-finger结构域通过与DNA的特定序列相互作用，能够增强BrPARP1蛋白与DNA的结合亲和力和特异性，从而提高其在DNA损伤修复和基因表达调控等过程中的效率。综上所述，BrPARP1基因的结构特征与其在DNA损伤修复、细胞凋亡、转录调控等生物学过程中的功能密切相关。通过对其结构的深入分析，为进一步研究BrPARP1基因在白菜型油菜干旱胁迫应答中的功能和作用机制奠定了坚实的基础。后续的研究可以围绕这些结构域展开，通过定点突变、结构解析等技术手段，深入探究各个结构域的具体功能和相互作用机制，从而全面揭示BrPARP1基因在植物抗旱过程中的分子机制。4.2BrPARP1基因在干旱胁迫应答中的表达调控在本研究中，通过实时荧光定量PCR技术，深入分析了干旱胁迫下白菜型油菜BrPARP1基因的表达模式。结果显示，在正常生长条件下，BrPARP1基因在叶片和根系中均有一定表达，但表达水平相对较低。随着干旱胁迫时间的延长，该基因在叶片和根系中的表达量均呈现出先上升后下降的趋势。在叶片中，干旱胁迫6h时，表达量开始显著上调，12h时达到峰值，随后逐渐下降；在根系中，表达量的上升响应时间略有延迟，12h时开始显著增加，24h时达到最高，之后逐渐回落。这种表达模式表明，BrPARP1基因能够对干旱胁迫作出快速且动态的响应。在干旱胁迫初期，基因表达量的上调可能是植物启动的一种自我保护机制。干旱胁迫会导致植物细胞内产生一系列的生理生化变化，如DNA损伤、活性氧积累、渗透失衡等。BrPARP1基因表达量的增加，可能通过其编码蛋白的功能，参与到这些生理过程的调节中。一方面，BrPARP1蛋白的DNA结合结构域可以识别并结合受损的DNA，启动DNA损伤修复机制，维持基因组的稳定性，确保细胞内遗传信息的准确传递和表达，为植物应对干旱胁迫提供遗传物质基础。另一方面，BrPARP1蛋白催化合成的聚[ADP-核糖]（PAR），可能在细胞信号传导和基因表达调控中发挥重要作用。PAR可以与多种蛋白质相互作用，调节这些蛋白质的活性和功能，进而影响细