猪中性粒细胞磷酸二酯酶基因表达及活性研究

 【目的】通过对猪中性粒细胞磷酸二酯酶（PDE）基因表达、活性及特异性抑制剂实验检测，研究其存在的主要PDE亚型及在体外的活性，为探讨PDE在中性粒细胞的分布、活性变化及对cAMP、cGMP调节与中性粒细胞功能关系提供依据。【方法】采用反转录聚合酶链反应（RT-PCR）检测PDE基因在猪中性粒细胞的表达，以高效液相色谱（HPLC）检测环核苷酸在PDE反应前后的含量变化，计算PDE活性。【结果】在检测的18个PDE亚型中，PDE1B、2A、3A、3B、4A、4B、4C、4D、5A、7A、7B、8A、8B、9A、11A 共15个PDE mRNA在猪中性粒细胞表达，特异性抑制剂实验表明PDE4、PDE5分别为水解cAMP、cGMP的主要PDE亚型；cAMP/cGMP-PDE在10～40 μl范围内，与其活性存在良好的线性关系（r＝0.9929/0.9992）。【结论】猪中性粒细胞至少存在15个PDE亚型，其中PDE4、PDE5为主要存在的2个亚型，PDE样品量在一定范围内与其活性存在良好的线性关系，可作为筛选PDE型中性粒细胞功能调节剂的方法。     

猪附红细胞体SYBR GreenⅠ荧光定量PCR检测方法的建立

为建一种快速、高效、灵敏的猪附红细胞体(M.suis)检测方法。本研究根据GenBank中的猪附红细胞体α-烯醇化酶(Eno)基因(序列号:AB265823.1),设计并合成了1对引物P5和P6,通过常规PCR方法对Eno基因中的一段序列进行扩增,并将其纯化后的PCR产物克隆入pMD-18T载体中,从而构建重组质粒。以阳性重组质粒为模板,建立了猪附红细胞体SYBR GreenⅠ荧光定量PCR检测方法,并对其灵敏性、特异性进行检测。在1.44×10~6~1.44×10~2拷贝/μL范围内具有较好的线性关系,相关系数R~2=0.965,检测到的浓度下限为1.44拷贝/μL。特异性试验中仅对猪附红细胞体的检测结果呈阳性,而对猪肺炎支原体、弓形虫、猪链球菌、猪呼吸与繁殖综合征病毒的检测结果均呈阴性。重复性试验结果表明,其组内和组间变异系数均小于2.0%。用该检测方法对49份猪附红细胞体疑似血液样品进行检测,所建立的荧光定量PCR方法检测的阳性率为40.8%(20/49),高出常规PCR方法检测的阳性率1.33倍。SYBR GreenⅠ荧光定量PCR检测方法的成功建立,为猪附红细胞体感染的早期诊断及定量分析提供了可靠依据。     

胡黄连与牛蒡子有效成分对猪磷酸二酯酶7A2活性的影响

从牛蒡子、胡黄连中筛选对猪磷酸二酯酶（PDE） 7A2具有抑制作用的中药单体成分,为PDE亚型抑制剂研究积累资料.使用高效液相色谱法（HPLC）测定绿原酸、牛蒡子苷、β-谷甾醇和胡黄连苷-Ⅱ对体外原核表达的猪PDE7A2的抑制作用.结果表明：在反应质量浓度为1,5,10,50,100 μmol/L时,绿原酸对猪PDE7A2的活性抑制率分别为12.55％,46.29％,56.17％,68.52％,47.12％;牛蒡子苷抑制率分别为21.60％,41.35％,70.16％,88.27％,80.07％;β-谷甾醇抑制率分别为40.53％,54.52％,71.81％,96.50％,94.86％;胡黄连苷-Ⅱ抑制率分别为69.34％,94.03％,96.50％,103.09％,101.44％.绿原酸、牛蒡子苷、β-谷甾醇和胡黄连苷-Ⅱ对体外表达的猪PDE7A2均有抑制作用,其中胡黄连苷-Ⅱ对猪PDE7A2的抑制作用最强.     

广西巴马小型猪PPARγ-2基因多态性与2型糖尿病易感性的相关性

[目的]探讨PPARγ-2基因多态性与广西巴马小型猪2型糖尿病易感性的关系。[方法]利用高糖高脂日粮饲喂24头广西巴马小型猪,利用PCR方法扩增PPARγ-2基因外显子2部分序列,测定其空腹血糖和胰岛素(INS)含量并进行糖耐量试验。[结果]克隆的PPARγ-2基因外显子2部分序列存在1处SNP位点(19813A/G),有AG(7头)和AA(17头)2种基因型。AG型广西巴马小型猪的空腹胰岛素含量和糖耐量试验中60 min血糖值显著高于AA型(P0.05),而AG型广西巴马小型猪2型糖尿病发生率(71.4%)明显高于AA型(5.9%)。[结论]广西巴马小型猪PPARγ-2基因外显子2的19813A/G多态性与2型糖尿病易感性有关。     

奶牛GPR109A基因SYBR GreenI实时荧光定量PCR检测方法的建立

【目的】建立奶牛G蛋白偶联受体109A基因(GPR109A)的实时荧光定量PCR检测方法,为开展奶牛产后脂类代谢紊乱及酮病发生机理研究奠定基础。【方法】根据Gen Bank中奶牛GPR109A基因和β-actin基因(内参基因)的m RNA保守序列设计合成两对引物,以奶牛肝脏组织c DNA为模板,PCR扩增目的基因后与载体连接构建重组质粒,以SYBR Green I实时荧光定量PCR进行扩增,绘制标准曲线,并对其特异性、敏感性、重复性等进行检验。【结果】GPR109A基因和β-actin基因的实时荧光定量PCR标准曲线方程分别为y=-3.309x+10.92(R2=0.9985)和y=-3.289x+9.794(R2=0.9997),扩增效率分别为101.0%和100.0%。特异性试验结果显示,GPR109A基因和β-actin基因的溶解曲线峰单一,无引物二聚体和非特异性扩增产物生成;敏感性试验结果显示,线性扩增范围广,可检测到1.8×102个拷贝的GPR109A基因;重复性试验结果显示,各扩增反应的变异系数(CV)小于或等于1.7%。【结论】基于GPR109A基因建立的SYBR GreenⅠ实时荧光定量PCR具有操作简便、敏感性高、特异性强和重复性好等特点,可用于牛源性GPR109A基因表达量的快速检测和定量分析。     

FAM213B基因启动子在猪子宫内膜细胞的转录调控

【目的】分析猪FAM213B基因启动子转录活性区域,并检测转录因子NFκB对启动子活性及猪FAM213B基因表达影响,为深入了解FAM213B基因的转录调控机制影响前列腺素合成、母猪妊娠等繁殖活动奠定基础。【方法】采集卵泡期子宫,分离子宫内膜上皮组织,经胶原酶方法获得猪原代子宫内膜细胞,用于猪FAM213B基因启动子活性检测。参照笔者所在课题组前期研究获得的猪FAM213B基因m RNA全序列(Gen Bank登录号:KX444503)以及5￠调控区序列(Gen Bank登录号:100134955),通过PCR扩增猪FAM213B基因较长启动子序列,并进行测序鉴定;在此基础上,PCR扩增带有Mlu I和Xho I限制性酶切位点的FAM213B基因启动子7个5'端缺失片段,并通过双荧光素酶报告基因载体系统构建猪FAM213B基因启动子7个不同的5'端缺失载体。将7个载体经无内毒素处理后与p RL-TK质粒用阳离子脂质体法一起转染猪子宫内膜细胞,用双荧光素酶报告基因系统进行Luciferase活性检测,比较各启动子片段转录活性。生物信息学分析FAM213B基因启动子潜在转录因子结合位点,通过Ch IP试验验证FAM213B基因启动子与潜在转录因子NFκB的相互作用。构建NFκB1和Rel A的超表达载体,化学合成NFκB1和Rel A的干扰si RNA片段,分别转染猪子宫内膜细胞,通过双荧光素酶报告基因系统进行Luciferase活性检测和荧光定量PCR分别检测超表达和干扰表达NFκB1和Rel A对猪FAM213B基因启动子活性和m RNA表达影响。【结果】PCR和测序获得猪FAM213B基因启动子长度为2 261bp(-2178/+83)。生物信息学预测结果表明猪FAM213B基因启动子区存在潜在的CREB、CCAAT增强子结合蛋白、E-box因子的结合位点,炎性因子NFκB潜在结合位点分别位于-1143/-1132和-664/-655区间。启动子报告基因活性检测结果表明:重组载体P2(-1352/+30)荧光活性最高,极显著高于P1(-1 760/+30)(P0.01),在-1 760/-1 352区域存在负调控元件;而P2显著高于P3(-919/+30)(P0.05),表明在-1352/-919区域存在正调控元件;P3(-919/+30)活性极显著高于P4(-604/+80)(P0.01),表明在-919/-604区域存在正调控元件;P4、P5、P6和P7活性差异不显著,-1352/-919区域为核心启动子区,-1352/-604对于维持该启动子较高转录活性起着重要的作用。Ch IP结果表明启动子区-1143/-1132存在NFκB1结合位点,-664/-655存在Rel A结合位点。超表达载体pc DNA3.1-NFκB1与P2(-1352/+30)启动子片段重组载体共转染猪子宫内膜细胞后,P2启动子活性极显著高于对照组(P0.01),pc DNA3.1-NFκB1载体单独转染猪子宫内膜细胞后FAM213B基因m RNA表达量显著高于对照组(P0.05);超表达载体pc DNA3.1-Rel A与P3(-919/+30)启动子片段重组体共转染猪子宫内膜细胞,P3启动子活性极显著低于对照组(P0.05),pc DNA3.1-Rel A载体单独转染猪子宫内膜细胞后,FAM213B基因的mRNA表达量显著低于对照组。转染NFκB1和Rel A的si RNA片段干扰片段,FAM213B基因启动子活性和m RNA表达量则表现与超表达相反的结果。【结论】获得了猪FAM213B基因启动子核心启动子序列为-1352/-919区域,NFκB是FAM213B基因启动子的转录因子,NFκB成员NFκB1和Rel A在猪子宫内膜细胞中对FAM213B基因的表达起调控作用。     

猪非典型性瘟病毒与猪流行性腹泻病毒双重PCR方法的建立与应用

为了快速准确地检测出猪非典型性瘟病毒(APPV)与猪流行性腹泻病毒(PEDV),根据Gen Bank中发布的PEDV S2基因和APPV NS3基因的序列,分别选取PEDV S2基因和APPV NS3基因的保守序列,设计、合成了1对特异性引物。通过对体系进行优化,分别扩增出190 bp和500 bp的特异性片段。建立的双重RT-PCR方法对猪繁殖与呼吸综合征病毒(PRRSV)、猪瘟病毒(CSFV)、猪轮状病毒(RV)的检测均为阴性。对PEDV和APPV的最低检出量分别为1μl 2.65×10~5和1μl 3.56×10~4,对四川遂宁、眉山等地采集到的腹泻、肌肉震颤的病料进行检测,PEDV、APPV阳性病料检出率分别为41. 18%和11. 76%,经分别与特异性检测APPV和PEDV的单一RT-PCR法检测结果进行比较分析,符合率均为100%。与单一RT-PCR检测方法具有相同特异性和重复性,可用于临床检测。     

猪圆环病毒2型荧光定量PCR检测试剂盒的研制及应用

为了开发针对猪圆环病毒2型(PCV2)的诊断技术及其试剂盒,对PCV2的ORF2基因设计相应的特异性引物和探针,采用荧光定量PCR方法检测PCV2。试验表明:该方法循环阈值(Ct)与标准DNA模板在10~1—10~7拷贝/μL浓度范围内具有极好的线性关系,相关系数为0.9999;重复性好,批内试验中Ct值的变异系数在0.59%—1.05%,批间试验的Ct值变异系数为1.9%—4.2%;研制的检测试剂盒与猪瘟病毒(HCV)、猪呼吸与繁殖障碍综合征病毒(PRRSV)、猪伪狂犬病病毒(PRV)不发生交叉反应,检测极限和定量极限分别为10和100拷贝数。用此方法对1 000个临床样品进行检测,结果710个被检测为阳性。结果显示,所研制的PCV2诊断试剂盒灵敏度高,特异性强,稳定性好,为PCV2的检测和流行病学的调查提供了一种简便快速的诊断方法。     

猪胚胎附植期Eph A4的组织表达及其多态性对产仔数的影响

为研究Eph A4对猪繁殖性能的影响,在猪妊娠第13天(附植前期)、第18天(中期)和第24天(后期),用PCR-SSCP方法检测了7个品种的950头母猪群中Eph A4基因的多态性,检测了Eph A4基因在大白猪子宫内膜等组织中的mRNA和蛋白表达量,并与猪的产仔数进行了关联分析.结果表明:Eph A4基因在子宫内膜附着点的mRNA表达量要高于任何其他组织的,且妊娠第13、18和24天的表达量差异不显著,但妊娠猪极显著低于空怀猪的(P＜0.01);Eph A4的蛋白表达量在附植前、中和后期的表达量两两间差异均显著(P＜0.05),且妊娠第18天最高,妊娠第24天最低.Eph A4基因第3外显子上检测到了2个突变位点EphA4_1和EphA4_2.EphA4_1是一个A→G的突变,使谷氨酸变成了赖氨酸,G为优势等位基因,且GG基因型关联的产仔数显著高于AG或GG基因型的(P＜0.05);EphA4_2是一个T→C的突变,C为优势等位基因,且CC基因型关联的产仔数显著高于TC或TT基因型的(P＜0.05).上述结果表明,Eph A4的表达对猪的胚胎附植发挥着一定的调控作用,EphA4_1位点GG基因型和EphA4_2位点CC基因型可以作为产仔数性状的标记基因型.     

犬中枢神经不同部位磷酸二酯酶基因表达和活性测定

为探明磷酸二酯酶(phosphodiesterase,PDE)在犬中枢神经活动中的作用,采用反转录聚合酶链反应检测18个PDE亚型在犬中枢神经不同部位中(大脑、小脑、延髓和脊髓各段)的表达分布,以高效液相色谱法检测环核苷酸在酶反应前后的含量变化,计算PDE活性。结果显示:在所测18个PDE亚型中,犬中枢神经的8个测定部位均有PDE1A,PDE1B,PDE1C,PDE2A,PDE3A,PDE3B,PDE4A,PDE4B,PDE4C,PDE4D,PDE5A,PDE7A,PDE7B,PDE8A,PDE8B,PDE9A,PDE11A表达,PDE10A均呈阴性表达,不同部位之间表达量差异较大,其中PDE4在各部位均呈优势表达;8个测定部位cAMP-PDE活性均高于cGMP-PDE活性,且不同部位PDE活性存在差异,其中脑组织cAMP-PDE活性和cGMP-PDE活性均高于脊髓各段。结果提示PDE可能在由cAMP和cGMP介导的中枢神经经信号传导中发挥着重要作用。     

玉米大斑病菌环腺苷酸磷酸二酯酶基因克隆及表达分析

【目的】获得玉米大斑病菌(Setosphaeria turcica)环腺苷酸磷酸二酯酶(c AMP phosphodiesterase,PDE)基因,并分析其在病菌侵染结构发育过程及侵染寄主早期阶段中的表达,为深入探索c AMP信号途径调控病菌致病性的作用机制打下基础。【方法】根据酿酒酵母(Saccharomyces cerevisiae)、白色念珠菌(Candida albicans)、灰霉病菌(Botrytis cinerea)、稻瘟病菌(Magnaporthe oryzae)、绿僵菌(Metarhizium anisopliae)和黑曲霉(Aspergillus niger)6种真菌PDE基因的保守序列,利用简并引物PCR对基因保守片段进行扩增,结合RACE和Genome Walking技术获得基因全长及其侧翼序列;利用MEGA 5.0软件对PDE蛋白预测编码产物进行多重序列比对,并采用邻近法构建系统发育树;利用GSDS分析基因结构,Prot Param分析理化性质,SOMPA预测二级结构,SMART数据库在线分析保守结构域;收集人造疏水介质诱导下侵染结构发育不同时期以及侵染感病寄主叶片不同时间的玉米大斑病菌材料,利用q RT-PCR技术研究StPDE在病菌侵染结构形成过程中不同阶段的转录水平。【结果】玉米大斑病菌基因组中存在1个高亲和力磷酸二酯酶基因(StH-PDE)和1个低亲和力磷酸二酯酶基因(StL-PDE)。其中,StH-PDE全长3 208 bp,包含5个内含子和6个外显子,ORF为2 898 bp。StL-PDE全长5054 bp,包含4个内含子和5个外显子,ORF为3 090 bp。StH-PDE和StL-PDE分别含有保守的Ⅰ型和Ⅱ型c AMP磷酸二酯酶催化结构域。不同的植物病原真菌中PDE同源基因分别呈现出高度的相似性,其中,StH-PDE与Magnaporthe grisea、Cordyceps militaris、Metarhizium acridum等病原真菌的高亲和力磷酸二酯酶基因聚于同一进化支,StL-PDE与Ascochyta rabiri、Scedosporium apiospermum、Fusarium oxysporum、Metarhizium album等病原真菌的低亲和力磷酸二酯酶基因聚于相同进化支。人造疏水介质诱导下,与菌丝相比,StH-PDE和StL-PDE在分生孢子中的表达水平均显著上调,其中StH-PDE和StL-PDE分别上调了约2倍和52倍。孢子萌发初期两个基因表达水平显著下调,随着萌发的进行,表达水平缓慢回升,至附着胞形成阶段,基因表达出现了第2次高峰,随后表达水平再次下调。在整个过程中,StH-PDE的最高表达水平则出现在附着胞形成时期,达到了菌丝体中的近7倍、分生孢子中的近2倍,而StL-PDE的表达水平始终低于其在分生孢子中的表达水平。StH-PDE和StL-PDE在病菌侵染寄主过程中的表达水平变化趋势与其在人工疏水介质诱导下的表现基本一致。在孢子萌发早期表达显著下调,随着时间延长,缓慢回升,至接种后18和24 h StH-PDE表达水平上调,超过萌发初期。【结论】在病菌侵染结构发育过程中,StH-PDE和StL-PDE均呈现下调-上调-下调的表达水平变化趋势。StH-PDE在附着胞时期表达水平最高,StL-PDE在分生孢子中的相对表达水平最高。     

Luteolin prevents f MLP-induced neutrophils adhesion via suppression of LFA-1 and phosphodiesterase 4 activity

Luteolin is an active ingredient found early from Folium perillae and Flos lonicerae, and has a specific inhibition on phosphodiesterase 4(PDE4) activity in vitro. Researches show luteolin has pharmacological effects of anti-inflammation, anti-anaphylaxis, antitumor, antioxidant, protection of nervous system and so on, and has mainly been used for the treatment of respiratory inflammatory diseases, cancer and cardiovascular disease in clinic. PDE4, specific to hydrolyze cyclic AMP(c AMP), is considered to be a new anti-inflammatory target due to the decisive role on c AMP signal in inflammatory cells such as neutrophils. In order to explore the anti-inflammatory mechanism, we further studied the effects of luteolin on the activity and expression of PDE4, the expression of lymphocyte function-associated antigen-1(LFA-1) and macrophage-1(MAC-1) in neutrophils, and the adhesion of neutrophils and endothelial cells. The results showed that luteolin had a dose-dependent inhibition on both bare PDE4 activity and PDE4 in cultured neutrophils, and had an obviously promotive effect on gene expressions of PDE4 A, 4B and 4D in later period. Luteolin had a significant inhibitory effect on neutrophils adhesion and LFA-1 expression in early stage, and had no obvious effect on MAC-1 expression. Therefore, luteolin can inhibit LFA-1 expression of neutrophils, then inhibit the adhesion of neutrophils and endothelial cells, and the mechanism is at least related with the inhibition of PDE4 activity.     

Influence of gallic acid on porcine neutrophils phosphodiesterase 4, IL-6, TNF-α and rat arthritis model

Our previous studies showed that the anti-inflammatory effects of Paeonia lactiflora roots extract may be mediated, at least in part, through its gallic acid content, and this effect may be regulated in part by an inhibition on cAMP-phosphodiesterase (PDE). To explore the anti-inflammatory effect and mechanism, the influence of gallic acid on neutrophils PDE4 activity and expression, TNF-α and IL-6 content and rat arthritis model were further studied. PDE4 activity and gene express were calculated respectively by substrate cAMP change examined with HPLC and real-time RT-PCR. The concentration of IL-6 and TNF-α in supernatant were assayed by ELISA method. Model of rat arthritis was caused by complete Freund’s adjuvant. Results showed that gallic acid had a dose-dependent restraint on PDE4 activity of neutrophils in vitro, promoted significantly PDE4A expression (P<0.01), and had no influence on the expressions of PDE4B and 4D. However, PDE4C expression was not detected. Gallic acid could promote IL-6 release (P<0.05), and inhibit TNF-α release of neutrophils (P<0.05). The experiment in vivo showed that gallic acid had obvious restraint on local inflammation of animal model (P<0.05). Therefore, the anti-inflammatory effect of gallic acid may be mediated in part through an inhibition on PDE4 activity and further an increase of IL-6 and a decrease of TNF-α of neutrophils, and this effect seemed to have no relationship with PDE4 expression.     

ToMV外壳蛋白互作IP-L蛋白的亚细胞定位及表达分析

【目的】前期研究显示To MV外壳蛋白(coat protein,CP)可以与烟草蛋白L(CP-interacting protein-L,IP-L)相互作用,本研究旨在明确To MV CP与IP-L的共定位情况、IP-L的组织表达和在To MV侵染条件下烟草IP-L及其编码蛋白的表达变化情况,为进一步明确IP-L的功能提供依据。【方法】通过双酶切法从笔者实验室构建保存的p GBKT7-CP和p GADT7-IP-L重组载体上切下To MV CP和IP-L目的片段,构建融合蛋白植物表达载体p PZP-IP-L-N-EGFP和p PZP-CP-N-Ds Red及原核表达载体p EGX-IP-L。通过热激法将融合表达的植物表达载体转化至农杆菌EHA105,在烟草表皮细胞瞬时表达IP-L-N-EGFP和CP-N-Ds Red,共聚焦荧光显微镜下观察IP-L和To MV CP共定位情况。用实时荧光定量RT-PCR(q RT-PCR)分析IP-L在本氏烟各组织部位的表达量。原核表达载体p EGX-IP-L,转化原核表达菌BL21,优化GST-IP-L可溶性表达条件后大量表达该蛋白,用GST亲和层析法纯化获得可溶性GST-IP-L蛋白,免疫家兔制备多克隆抗体。ELISA测定抗体效价,Western印迹法明确抗体的特异性后,利用该抗体分析To MV侵染条件下番茄IP-L蛋白表达情况,用q RT-PCR检测IP-L表达变化与蛋白水平是否一致。【结果】To MV CP在本氏烟叶片表皮细胞的细胞质和细胞质膜以及叶绿体均有分布,而IP-L只在本氏烟叶片表皮细胞质膜表达,二者共定位在细胞质膜。IP-L在本氏烟叶片中特异高表达,显著高于茎、根和花中的相对表达量。在温度为30℃,IPTG诱导浓度为0.3 mmol·L-1条件下表达出大小为42.8 k D的可溶性GST-IP-L融合蛋白。纯化获得约4.2 mg可溶性蛋白,免疫家兔制备了效价为1/6400的多克隆抗体。Western印迹结果表明,该抗体可以与IP-L的原核产物特异性结合。在To MV侵染本氏烟1、3、7 d后,Western印迹分析表明IP-L在叶片内表达量随接种时间呈明显上调趋势。q RT-PCR检测结果与Western印迹结果一致,显示IP-L在To MV侵染本氏烟第7天后的叶片内表达量是健康对照的3倍多,差异达到显著水平。【结论】To MV CP和IP-L共定位于本氏烟表皮细胞质膜,IP-L在本氏烟叶片内特异高表达。制备的IP-L多克隆抗体具有良好的免疫反应活性,可以用于IP-L表达量检测;To MV侵染烟草可诱导IP-L及其编码蛋白表达量显著上调。     

Identification of suitable reference genes in leaves and roots of rapeseed(Brassica napus L.)under different nutrient deficiencies

Nutrient deficiency stresses often occur simultaneously in soil. Thus,it's necessary to investigate the mechanisms underlying plant responses to multiple stresses through identification of some key stress-responsive genes. Quantitative real-time PCR(qRT-PCR)is essential for detecting the expression of the interested genes,of which the selection of suitable reference genes is a crucial step before qRT-PCR. To date,reliable reference genes to normalize qRT-PCR data under different nutrient deficiencies have not been reported in plants. In this study,expression of ten candidate reference genes was detected in leaves and roots of rapeseed(Brassica napus L.)after implementing different nutrient deficiencies for 14 days. These candidate genes,included two traditionally used reference genes and eight genes selected from an RNASeqdataset. Two software packages(Ge Norm,Norm Finder)were employed to evaluate candidate gene stability. Results showed that VHA-E1 was the highest-ranked gene in leaves of nutrient-deficient rapeseed,while VHA-G1 and UBC21 were most stable in nutrient-deficient roots. When rapeseed leaves and roots were combined,UBC21,HTB1,VHA-G1 and ACT7 were most stable among all samples. To evaluate the stabilities of the highest-ranked genes,the relative expression of two target genes,Bn Trx1;1 and Bn Pht1;3 were further determined. The results showed that the relative expression of Bn Trx1;1 depended on reference gene selection,suggesting that it's necessary to evaluate the stability of reference gene prior to qRT-PCR. This study provides suitable reference genes for gene expression analysis of rapeseed responses to different nutrient deficiencies,which is essential for elucidation of mechanisms underlying rapeseed responses to multiple nutrient deficiency stresses.