Effects of applying compound fertilizer on yield, nutrient absorption and oligosaccharide accumulation of Morinda officinalis
-
摘要:目的
研究复合肥施用量对巴戟天产量、养分吸收和主要有效成分寡糖累积量的影响，为巴戟天的施肥管理提供理论依据。
方法以巴戟天为研究对象，通过盆栽试验，设置每株施用复合肥0、2、4、6、8和10 g (CK、T1、T2、T3、T4和T5)共6个处理，测定巴戟天在不同施肥处理下的产量、养分含量和寡糖含量，分析各指标间的相关性，并运用主成分分析方法进行综合评价。
结果施用复合肥可以显著提高巴戟天产量和N、P、K含量，并有利于提高寡糖的累积量。其中，单株施肥量为6 g处理的巴戟天产量和寡糖总累积量最高，与CK相比，分别提升114.47%和75.36%。施肥在一定程度上降低了巴戟天的寡糖含量，与CK相比，降幅最高达35.65%。巴戟天对复合肥中养分的利用效率表现为N > K > P，其中T3处理N、P、K累积量均高于其他处理。影响巴戟天寡糖累积量的最关键因素是产量，巴戟天产量与全株N/P/K含量、4种寡糖累积量呈显著或极显著正相关，但耐斯糖、1F−果呋喃糖基耐斯糖含量与产量、全株N/P/K含量呈显著或极显著负相关。主成分分析的综合得分由高到低为T3 > T5 > T4 > T2 > T1 > CK，T3处理最佳。
结论建议巴戟天的种植以N肥为主，适量补充K和P肥，适当施肥可显著提高巴戟天药材产量和寡糖累积量。
Abstract:ObjectiveTo study the effects of compound fertilizer application amount on yield, nutrient absorption and accumulation of main effective component oligosaccharides of Morinda officinalis, and provide a theoretical basis for the fertilization management of M. officinalis.
MethodUsing M. officinalis as the research object, through the pot experiment, a total of six treatments were set up to apply 0, 2, 4, 6, 8 and 10 g compound fertilizer (CK, T1, T2, T3, T4 and T5) per plant. The yield, nutrient content and oligosaccharide content of M. officinalis under different fertilization treatments were determined. Then we analyzed the correlation between the indicators, and comprehensively evaluated them by principal component analysis.
ResultThe application of compound fertilizer significantly improved the yield of M. officinalis, increased the N, P and K contents, and contributed to the increase in the accumulation of oligosaccharides. The highest yield and total oligosaccharide accumulation were achieved in the treatment of 6 g compound fertilizer per plant (T3), compared with no fertilizer (CK), with increases of 114.47% and 75.36% respectively. Fertilizer application reduced the oligosaccharide content of M. officinalis to some extent, with a maximum decrease of 35.65% compared to CK. The utilization efficiency of nutrient in compound fertilizer by M. officinalis showed N > K > P, and the accumulations of N, P and K in T3 treatment were higher than those of other treatments. The yield was the most critical factor affecting the accumulation of oligosaccharides in M. officinalis. The yield of M. officinalis was significantly or highly significantly positively correlated with the whole plant N, P, K contents and the accumulation of four oligosaccharides, but the contents of nystose and 1F-fructofuranosyl nystose were significantly or highly significantly negatively correlated with the yield and the whole plant N, P, K contents. The comprehensive score of principal component analysis from high to low was T3>T5>T4>T2>T1>CK, with T3 treatment being the best.
ConclusionIt is suggested that the cultivation of M. officinalis should focus on N fertilizer, and supplement K and P fertilizer appropriately. Proper fertilization can significantly increase the yield of M. officinalis and the accumulation of oligosaccharides.
-
Keywords:
- Morinda officinalis /
- Compound fertilizer /
- Yield /
- Oligosaccharide /
- Nutrient absorption
-
大豆Glycine max是重要的粮油作物，在世界范围内广泛种植，为食品和畜牧业提供了重要的蛋白和油脂来源，也是生物柴油生产的主要原料[1-2]。我国南方水热条件充足，但土壤大多呈酸性，酸性土壤中磷的有效性低是限制大豆生产的重要因子[3-4]。研究发现改变根系形态构型是大豆适应低磷胁迫的重要机制之一[5-6]。低磷条件下，大豆根冠比、根长度、根毛密度以及侧根数量均显著提高。根系形态构型的这些变化，增加了植物根系与土壤接触的面积，促进了植株对土壤磷的吸收和利用[7]。另外，赵静等[8]研究表明，浅根型大豆品种的生物量、磷含量和产量均明显高于深根型大豆品种，说明浅根型大豆对耕层土壤磷的吸收效率更高。因此，根系的向重性反应对于植株利用土壤中异质分布的磷至关重要[9-12]。
生长素是植物根系形态建成中的重要调控因子[13]。前期研究显示，生长素的不对称分布在植物的根系生长和根系向重性反应中均起到了重要的调控作用[13-14]。植物的PIN(PIN-FORMED)基因家族是介导植物根系生长素不对称分布的重要因子[15]。模式植物拟南芥Arabidopsis thaliana的基因组中存在8个AtPIN同源基因[16]。其中，AtPIN2介导了生长素从根尖侧皮层细胞到根伸长区皮层细胞的极性运输[17]。研究显示，在拟南芥根系向重性反应的过程中，根尖细胞的AtPIN2蛋白发生不对称分布，使生长素在根系向地面的伸长区积累，抑制了向地侧伸长区细胞的伸长，最终导致根系生长弯曲，产生向重性反应[14, 18]。因此，AtPIN2突变后，突变体根部下侧的生长素分布受到抑制，根系的向地性功能明显缺失[19]。水稻、番茄、玉米、毛白杨、葡萄和桉树等植物的PIN2同源基因也有类似的保守结构或功能[20-23]。例如，水稻的同源基因OsPIN2可以回补拟南芥突变体Atpin2的表型缺陷，且OsPIN2的回补表达可以恢复突变体Atpin2在向重情况下根尖生长素的极性分布，说明OsPIN2通过影响根尖生长素的极性转运调控了根系生长夹角，在水稻根系的重力反应中起重要作用[20, 24-25]。
虽然关于PIN2基因在其他植物中的功能已取得较大的研究进展，但大豆中PIN2的同源基因是否也具有类似的功能，且是否参与了低磷胁迫对大豆根系形态建成的调控尚未清楚。因此，本研究首先对大豆GmPIN基因家族进行生物信息学分析；然后，选取与拟南芥AtPIN2同源性较高的GmPIN2b作为研究对象，利用实时荧光定量PCR测定该基因响应低磷胁迫的表达模式；进一步在拟南芥Atpin2突变体中回补表达该基因，并对转基因植株的生长和根系向重性进行分析。
1. 材料与方法
1.1 材料
本试验采用的大豆材料为‘粤春03-3’，供试拟南芥为哥伦比亚野生型Col-0和Atpin2突变体CS8058。拟南芥Atpin2突变体由浙江大学生命科学学院毛传澡教授课题组惠赠。
1.2 试验方法
1.2.1 生物信息学分析
以AtPIN2氨基酸序列为参考，在植物基因组数据库Phytozome ( http://www.phytozome.org/index.php)中通过BLAST比对获得拟南芥、大豆和水稻Oryza sativaL.等植物的PIN同源基因的氨基酸序列，再利用软件MEGA5.05的邻接法(Neighbor-joining method)生成进化树。在SOYBASE数据库( https://www.soybase.org/soyseq/)中获得GmPIN基因的组织表达数据，在已发表的大豆根系转录组测序结果[26]中获得不同磷处理条件下GmPIN基因的表达水平的数据，通过TBtools软件将这些数据制作成表达模式热图。
1.2.2 大豆低磷胁迫处理
挑选种皮无破损、大小均一的大豆种子，利用10%(φ)次氯酸钠溶液表皮消毒后，在沙子中萌发。待幼苗子叶张开，将幼苗转移至Hoagland营养液中进行不同磷水平的处理。其中，高磷处理的KH2PO4浓度为250 μmol/L，低磷处理的KH2PO4浓度为25 μmol/L。处理3、6、12和15 d后，分析植株鲜质量、总根长和根冠比，其中根冠比的计算方法为地下部鲜质量与地上部鲜质量的比值；同时收取根系样品用液氮速冻后用于基因表达水平的测定。每个处理设4次重复。
1.2.3 大豆GmPIN2b基因的克隆
从Phytozome数据库中获得大豆GmPIN2b(Glyma.17G057300)的CDS序列，根据该序列设计正向引物(ATGATCACTGGTAAGGATATTTATG)和反向引物(TTAAACTCCAAGCAGCACGTA)。以大豆‘粤春03-3’根系cDNA为模板，扩增GmPIN2b基因，将所得的PCR产物连接T载体，测序验证，选取阳性克隆保存备用。
1.2.4 回补表达GmPIN2b的拟南芥Atpin2突变体材料的构建
以拟南芥根系基因组DNA为模板，设计正向引物(TATGACCATGATTACGAATTCGCCCCTAGACCAAAATAATATGG)和反向引物(CATAAATATCCTTACCAGTGATCATTTTGATTTACTTTTTCCGGC)扩增拟南芥AtPIN2基因，翻译起始密码子上游3 229 bp的启动子片段。以大豆根系cDNA为模板，设计正向引物(GCTACAAACCCTGCATTAACCATGATCACTGGTAAGGATATTTATG)和反向引物(CCGGGTACCGAGCTCGAATTCTTAAACTCCAAGCAGCACGTA)扩增GmPIN2b基因编码框全长。将扩增产物进行凝胶电泳，利用琼脂糖凝胶DNA回收试剂盒(美基生物，上海)进行纯化回收。用SalI对目标载体pCAMBIA1300进行酶切，获得线性化的载体。用Exnase® Multis(诺唯赞生物，中国)将上述PCR纯化产物与线性化pCAMBIA1300进行多片段连接，转化大肠埃希菌，经过测序无误后抽取质粒转化农杆菌GV3101。
将GV3101平板划线活化，挑取单克隆接种于添加卡那霉素和利福平添加剂量均为10 μg/mL的10 mL YEP培养液中，放置培养箱于180 r/min、28 ℃条件下培养过夜。之后取3 mL菌液转入100 mL YEP培养液中扩大培养，然后5000 r/min离心10 min，收集菌体，用等体积的1/4MS培养液重悬菌体并加入0.02%(φ)SilwetL-77配成拟南芥转化液。利用花絮浸染法转化拟南芥突变体Atpin2，获得pAtPIN2::GmPIN2b回补表达的拟南芥转基因株系，T3代纯合后用于进一步试验分析。拟南芥AtPIN2基因的检测引物为At-F(正向引物)：CCTTGCTTGGTCCCTTGTCT，At-R(反向引物)：ATCGCAAACCCTGCTACTGA；pAtPIN2::GmPIN2b回补表达的拟南芥转基因株系的检测引物为At-F1(正向引物)：GTGCATGAACGACTACAATGG，Gm-R1(反向引物)：GAAAGAAGTGGAACTGCGAAC；转基因拟南芥中GmPIN2基因表达水平的检测引物为Gm-F2(正向引物)：GATTTCTATGCCATGTTCGCAA，Gm-R2(反向引物)：CTCAATCAATTCGTGAACAGCT。
1.2.5 GmPIN2b回补拟南芥Atpin2突变体表型和根系向重性的分析
将野生型、Atpin2突变体以及pAtPIN2::GmPIN2b转基因株系的种子表皮灭菌后，在1/2MS固体培养基中萌发。待幼苗根长至1.5 cm左右，将生长一致的幼苗转移到新的1/2MS固体培养基中进行不同磷水平处理。其中高磷处理含KH2PO4 1 250 μmol/L，低磷处理含KH2PO4 25 μmol/L；培养箱温度为24 ℃；光周期为16 h/8 h。磷处理8 d后，分析植株鲜质量、主根长和侧根数。另外，将长势一致的幼苗转移到新的1/2MS固体培养基中进行上述磷处理，然后将培养皿逆时针旋转90°放置，处理30 h后扫描植株根系，并用ImageJ(National Institutes of Health，美国)测定根系弯曲角度。每个处理设4次重复，每次重复8棵植株。
1.2.6 植物RNA的提取及反转录
取约0.1~0.2 g植物鲜样用液氮研磨，根据Trizol(Omega，美国)的说明书在样品中加入预冷的1 mL Trizol震荡混匀；再加入0.2 mL氯仿，震荡混匀后，室温静置5 min。离心后吸取0.5 mL上清液，往上清液中加入等体积的异丙醇，冰上静置后离心，所得的沉淀用75%(φ)乙醇溶液洗涤2次，沉淀风干后加入无RNA酶的灭菌水溶解RNA。用分光光度计(IMPLEN，美国)测定RNA浓度，然后参照反转录试剂盒(Promega，美国)说明书，合成cDNA第1链，反应结束后，−20 ℃保存备用。
1.2.7 实时荧光定量PCR
设计GmPIN2b的定量PCR引物为F(正向引物)：AATTGCATTGCCCATAACCATACTC，R(反向引物)：ATTTCCCTTCATCCACCCCGC；大豆看家基因GmEF1a的定量PCR引物为F(正向引物)：TGCAAAGGAGGCTGCTAACT，R(反向引物)：CAGCATCACCGTTCTTCAAA和拟南芥AtEF1a的定量PCR引物为F(正向引物)：GTCGATTCTGGAAAGTCGACC，R(反向引物)：AATGTCAATGGTGATACCACGC。将所得第1链稀释20倍做为定量PCR反应模板；适量原液做梯度稀释为标准曲线模板。定量PCR反应体系：10 μL 2× Go Taq®qPCR master mix(Promega，美国)；正向引物和反向引物各0.4 μL；2 μL模板；0.2 μL CXR Reference；补H2O至 20 μL。反应程序：95 ℃预变性10 min；95 ℃ 15 s，58 ℃ 15 s，72 ℃ 1 min，40个循环。反应结束后用实时荧光定量PCR仪(Thermo Fisher，美国)携带的Real-Time Analysis Software 6.0计算每个样品的表达量。相对表达量为目的基因的表达量与看家基因表达量的比值。
1.2.8 数据统计与分析
本文所有数据均使用Microsoft Excel 2016(Microsoft Company，美国)和IBM SPSS Statistics V21.0(IBM SPSS，美国)进行处理及统计分析。利用TBtools软件分析基因的表达模式和制作基因表达模式热图。
2. 结果与分析
2.1 大豆GmPIN的系统进化分析及组织表达模式分析
利用MEAG 5.05对大豆、水稻和拟南芥的PIN蛋白进行系统进化树的构建。进化树分析结果显示，拟南芥、水稻和大豆的PIN家族在进化树上可分为7个亚族(图1)。其中，大豆GmPIN2a和GmPIN2b与拟南芥AtPIN2位于同一亚族，二者相似性超过90%。对SOYBASE数据库中23个大豆GmPIN基因的组织表达水平构建表达模式热图，结果(图2)显示，GmPIN1e、GmPIN3b、GmPIN1a和GmPIN1c等基因在各组织部位中均有较高的表达水平；而GmPIN2b、GmPIN2a与GmPIN5b等基因在根部的表达水平较高。
![]()
图 1 GmPIN系统进化分析Figure 1. Phylogenetic analysis of GmPINAt: Arabidopsis thaliana, Gm: Glycine max, Os: Oryza sativa![]()
图 2 GmPIN基因组织表达模式1：新叶，2：花，3：1 cm豆荚，4：花后10 d的豆荚壳，5：花后14 d的豆荚壳，6：花后10 d的种子，7：花后14 d的种子，8：花后21 d的种子，9：花后25 d的种子，10：花后28 d的种子，11：花后35 d的种子，12：花后42 d的种子，13：根，14：根瘤Figure 2. Tissue expression pattern of GmPIN genes1: Young leaf, 2: Flower, 3: Pod of 1 cm length, 4: Pod shell of 10 d after flowering, 5: Pod shell of 14 d after flowering, 6: Seed of 10 d after flowering, 7: Seed of 14 d after flowering, 8: Seed of 21 d after flowering, 9: Seed of 25 d after flowering, 10: Seed of of 28 d after flowering, 11: Seed of 35 d after flowering, 12: Seed of 42 d after flowering, 13: Root, 14: Root nodule2.2 低磷胁迫对GmPIN家族成员在大豆根系中表达模式的影响
利用已报道的大豆根系转录组数据[26]，对低磷调控的19个GmPIN基因的表达模式进行了分析。结果(图3)表明，GmPIN9d、GmPIN8d、GmPIN8a、GmPIN2a与GmPIN2b这5个基因均受低磷胁迫上调表达，其中，低磷条件下GmPIN2b的表达量是高磷对照的1.2倍；相反，GmPIN3d、GmPIN9b、GmPIN5a与GmPIN6b这4个基因受低磷胁迫下调表达，其中，GmPIN3d的表达水平在低磷条件下仅为高磷对照的50%；除此之外，其他GmPIN基因成员的表达水平不受低磷胁迫的影响。
![]()
图 3 低磷胁迫下大豆根系中GmPIN基因表达模式Figure 3. Expression pattern ofGmPIN genes in soybean roots in response to low phosphorus stress2.3 不同低磷胁迫时间对大豆生长和根系中GmPIN2b表达水平的影响
根据SOYBASE数据库中GmPIN同源基因的组织表达模式，选取在根组织(包括根和根瘤)中均有表达的GmPIN2b作为研究对象，进一步分析不同磷处理时间对该基因表达水平的影响。植物表型分析结果(图4)显示，低磷处理12和15 d后，与高磷对照相比，大豆单株鲜质量分别降低了11%和33%(P<0.01)(图4A)；单株总根长分别为高磷对照的1.18和1.32倍(P<0.05)(图4B)；根冠比极显著增加，分别为高磷对照的1.43和1.99倍(P<0.01)(图4C)。定量PCR结果显示，低磷处理6、12和15 d后，大豆根系中GmPIN2b的表达水平均显著上调，分别是高磷对照的4.19、1.77和1.93倍(P<0.05)(图4D)。
![]()
图 4 不同低磷处理时间对大豆生长和GmPIN2b表达水平的影响“*”和“**”表示相同胁迫时间高磷对照和低磷处理在P<0.05和P<0.01水平差异显著(t检验)Figure 4. Soybean growth and GmPIN2b expression level in response to different periods of low phosphorus treatment“*” and “**” indicate significant differences at P<0.05 and P<0.01 levels between high phosphorus control and low phosphorus treatment at the same stress time (t test)2.4 转基因拟南芥检测
将拟南芥AtPIN2启动子与大豆GmPIN2b的编码区融合，构建了pAtPIN2::GmPIN2b表达载体，转化拟南芥突变体Atpin2，获得了回补表达GmPIN2b的转基因材料。对转基因材料进行了转基因鉴定和GmPIN2b表达水平的分析，结果显示，转基因材料扩增出pAtPIN2启动子融合GmPIN2b基因CDS的片段，表明GmPIN2b回补表达转基因材料为阳性植株(图5)。GmPIN2b在阳性株系RA1和RA2的根系和叶片中均有较高表达水平，与野生型相比呈显著差异(P<0.001)(图6)。
![]()
图 5 GmPIN2b阳性转基因植株的检测M：Marker，WT：拟南芥Col-0野生型，Atpin2：拟南芥突变体，RA1和RA2：pAtPIN2::GmPIN2b回补表达株系Figure 5. Detection of GmPIN2b positive transgenic plantsM: Marker, WT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b![]()
图 6 拟南芥回补表达株系根系和叶片中GmPIN2b的表达水平WT：拟南芥Col-0野生型，Atpin2：拟南芥突变体，RA1和RA2：pAtPIN2::GmPIN2b回补表达株系；“***”表示回补表达株系与野生型在P<0.001水平差异显著(t检验)Figure 6. Expression levels of GmPIN2b in the shoots and roots of Arabidopsis complement expression linesWT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b; “***” indicates significant differences at P<0.001 level between complement expression lines and wild-type (t test)2.5 回补表达GmPIN2b对转基因拟南芥生长的影响
野生型、pAtPIN2::GmPIN2b转基因株系和Atpin2突变体等材料的种子在1/2MS培养基上萌发5 d后，进行不同磷水平的处理，8 d后对其单株鲜质量、主根长度和侧根数等进行分析。与高磷对照相比，低磷胁迫下植株较小，主根生长均受到明显抑制(图7)。结果显示，低磷条件下，无论是野生型还是转基因材料，单株鲜质量和主根长均比高磷对照显著降低(P<0.001)(表1)。无论在高磷还是低磷条件下，AtPIN2基因的突变均显著降低了拟南芥单株的鲜质量和主根长，而回补表达GmPIN2能够部分恢复Atpin2的鲜质量和主根长突变表型(图8A、8B)。其中，回补表达株系RA1和RA2的植株鲜质量在高磷和低磷条件下分别比突变体高27%、28%和18%、16%，其主根长度分别比突变体提高33%、39%和43%、18%。另外，低磷显著减少了拟南芥植株的一级侧根数(P<0.001)(表1)。但在高磷条件下，植株一级侧根数在不同材料之间无显著差别；在低磷条件下，Atpin2突变体的一级侧根数显著低于野生型，而回补表达GmPIN2b后，转基因株系RA1和RA2的一级侧根数均比突变体显著增加，分别为突变体的1.24和1.20倍(图8C)。
![]()
图 7 回补表达GmPIN2b对转基因拟南芥生长的影响WT：拟南芥Col-0野生型，Atpin2：拟南芥突变体，RA1和RA2：pAtPIN2::GmPIN2b回补表达株系；标尺 = 1 cmFigure 7. Effect of GmPIN2bcomplement expression on growth of transgenic ArabidopsisWT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b; Bar = 1 cm表 1 不同磷处理和转基因回补表达对拟南芥生理指标和根系向重性影响的双因素方差分析Table 1. Two-way ANOVA of the effects of different phosphorus treatments and transgenic complement expression on Arabidopsis physiological indicators and root gravity指标
IndicatorF1) 磷处理
Phosphorus treatment转基因回补表达
Transgenic
complement expression磷处理×转基因回补表达
Phosphorus treatment ×
transgenic complement expression单株鲜质量 Fresh weight per plant 210.442*** 35.232*** 23.287*** 单株主根长 Main root length per plant 454.853*** 41.777*** 17.297*** 单株一级侧根数
Primary lateral root number per plant44.212*** 5.271* 1.383 弯曲角度 Angle of curvature 4.841* 24.205*** 1.549 1)“*”：P<0.05，“**”：P<0.01，“***”：P<0.001(Duncan’s法)
1)“*”: P<0.05, “**”:P<0.01, “***”:P<0.001(Duncan’s method)![]()
图 8 回补表达GmPIN2b对转基因拟南芥生长指标的影响WT：拟南芥Col-0野生型，Atpin2：拟南芥突变体，RA1和RA2：pAtPIN2::GmPIN2b回补表达株系；各小图中相同磷处理柱子上的不同小写字母表示株系间差异显著(P < 0.05, Duncan’s法)Figure 8. Effects of GmPIN2bcomplement expression on growth indexes of transgenic ArabidopsisWT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b; Different lowercase letters on the columns of the same phosphorus treatment in each figure indicate significant differences among different lines (P < 0.05, Duncan’s method)2.6 回补表达GmPIN2b对转基因拟南芥根系向重性的影响
将萌发了5 d、生长一致的幼苗转移到新的1/2MS固体培养基中进行不同磷水平处理，将培养皿逆时针旋转90°放置，处理30 h后观察拟南芥主根的弯曲角度。从拟南芥的表型可见，高磷条件下野生型根系的弯曲程度较大，突变体Atpin2较小；而低磷条件下，所有植株的根系生长均明显受阻，但野生型根系的弯曲程度较其他根系明显(图9)。拟南芥根系向重性受磷处理显著影响(P<0.05)，在转基因回补表达株系间具有极显著差异(P<0.001)(表1)。无论在高磷还是低磷条件下，RA1、RA2和Atpin2的主根弯曲角度均显著小于野生型(图10)。高磷条件下，RA1和RA2主根弯曲的角度显著高于Atpin2，分别是Atpin2的2.33和2.10倍；而低磷条件下，回补表达株系的主根弯曲角度稍微高于Atpin2，差异不显著(图10)。
![]()
图 9 回补表达GmPIN2b对转基因拟南芥根系表型的影响WT：拟南芥Col-0野生型，Atpin2：拟南芥突变体，RA1和RA2：pAtPIN2::GmPIN2b回补表达株系；标尺=1.5 cmFigure 9. Effect of GmPIN2bcomplement expression on root phenotype of transgenic ArabidopsisWT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b; Bar=1.5 cm![]()
图 10 回补表达GmPIN2b对转基因拟南芥主根弯曲角度的影响WT：拟南芥 Col-0 野生型，Atpin2：拟南芥突变体，RA1 和 RA2：pAtPIN2::GmPIN2b 回补表达株系；各小图中相同磷处理柱子上的不同小写字母表示株系间差异显著 (P < 0.05, Duncan’s法)Figure 10. Effect of GmPIN2b complement expression on curvature angle of transgenic ArabidopsisWT: Arabidopsis Col-0 wild-type, Atpin2: Arabidopsis mutant, RA1 and RA2: Complement expression lines of pAtPIN2::GmPIN2b; Different lowercase letters on the columns of the same phosphorus treatment indicate significant differences among different lines (P < 0.05, Duncan’s method)3. 讨论与结论
生长素作为一种重要的生长调节激素参与植物根系的生长发育[13-14]。植物生长素输出载体PIN基因家族在介导植物生长素不对称分布、调控植株形态建成的过程中起着关键的作用。例如，拟南芥的基因组中存在8个PIN同源基因，这些基因参与拟南芥植株发育的重要调控过程，包括调控胚和花粉的发育以及维管组织的分化[27-28]，介导生长素在茎尖分生组织的环流和从地上部到根系的极性运输[27, 29-30]，调控细胞内生长素的动态平衡[28, 31]，参与根的向重性生长等[19, 32]。本研究对大豆的PIN同源基因进行分析，发现大豆基因组中存在23个GmPIN同源基因，其中大豆GmPIN2b与拟南芥AtPIN2位于同一亚族，具有较高的相似性。AtPIN2参与了根系向地生长的调控[13-14]，但是大豆GmPIN调控根系生长和向重性的功能未见报道。因此，本研究选取与AtPIN2较为相似的GmPIN2b作为研究对象，获得了回补表达GmPIN2b的拟南芥Atpin2突变体的转基因材料。该转基因材料能够部分恢复Atpin2突变体根系向重性的表型。该结果与前人的研究结果相似，例如，水稻OsPIN2可以回补突变体Atpin2主根的向重性功能缺陷，且OsPIN2的回补表达可以恢复Atpin2在向重情况下根尖生长素的极性分布[20]。GmPIN2b能够回补Atpin2突变体的根系向重性这一结果进一步说明了植物PIN2在调控植物根系向重性反应的过程中具有一定的功能保守性。
早期研究表明，植物根系形态构型的改变是植物适应低磷胁迫的重要机制[33]。本研究显示，低磷处理12 d后的大豆总根长高于高磷对照，这与前人的研究结果相同[12, 33]，说明低磷胁迫条件下，大豆通过促进根系生长、提高植株根冠比，从而提高植株对外界环境磷的吸收和利用。另外，低磷胁迫显著提高了GmPIN2b在大豆根系的表达水平，回补表达GmPIN2b能够部分回补拟南芥Atpin2突变体的向重性。这些结果暗示了GmPIN2b介导的生长素不对称分布可能参与了低磷胁迫调控大豆根系生长的过程；因此，本研究观察了回补表达GmPIN2b的拟南芥Atpin2突变体在不同磷水平条件下根系的生长情况。研究发现无论高磷还是低磷条件下，与Atpin2突变体相比，回补表达GmPIN2b均明显增加了转基因植株的鲜质量和主根长；尤其是低磷条件下回补表达GmPIN2b显著提高了突变体的一级侧根数。最近有研究显示，油菜素甾醇促进PIN2蛋白从根尖到伸长区的积累，导致侧根数增加[19, 34]，说明PIN2蛋白也参与了侧根发育的调控；因此，低磷条件下，大豆可能通过提高GmPIN2b的表达水平，改变主根尖生长素的向基运输以及侧根原基中生长素的浓度，从而参与低磷条件下侧根的发育。
综上所述，低磷胁迫提高了大豆根系中GmPIN2b的表达，回补表达该基因能够部分回补Atpin2突变体的生物量、侧根数和根系向重性；因此，GmPIN2b在根系形态建成响应低磷胁迫的过程中可能起着重要的调控作用。本研究结果为进一步揭示大豆根系适应低磷胁迫的分子机制提供了重要的理论依据。
-
![]()
图 1 复合肥不同施用量对巴戟天不同部位N、P、K含量的影响
各小图中，相同部位柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 1. Effects of different application amounts of compound fertilizer on contents of N, P and K in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns of the same part indicate significant differences among different treatments (P<0.05, Duncan’s method)
![]()
图 2 复合肥不同施用量对巴戟天各部位N、P、K单株累积量的影响
各小图中，不同小写字母表示处理间地上部、根部养分累积量差异显著，不同大写字母表示处理间全株养分累积量差异显著(P<0.05，Duncan’s法)
Figure 2. Effects of different application amounts of compound fertilizer on the accumulation of N, P and K per plant in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns indicate significant differences in aboveground and root nutrient accumulation among treatments, different capital letters indicate significant differences in whole plant nutrient accumulation among treatments (P < 0.05, Duncan’s method)
![]()
图 3 复合肥不同施用量对巴戟天4种寡糖含量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 3. Effects of different application amounts of compound fertilizer on contents of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 4 复合肥不同施用量对巴戟天4种寡糖单株累积量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 4. Effects of different application amounts of compound fertilizer on accumulation of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 5 巴戟天不同指标的相关性分析
“*”和“**”分别表示在P<0.05和P<0.01水平显著相关(Pearson法)
Figure 5. Correlation analyses of different indicators for Morinda officinalis
“*” and “**” indicate significant correlations at P < 0.05 and P < 0.01 levels, respectively (Pearson method)
![]()
图 6 复合肥不同施用量的主成分分析
PC1：第一主成分，PC2：第二主成分
Figure 6. Principal component analyses of different application amounts of compound fertilizer
PC1: The first principal component, PC2: The second principal component
表 1 复合肥不同施用量对巴戟天生长指标的影响1)
Table 1 Effects of different application amounts of compound fertilizer on growth indexes of Morinda officinalis
处理
Treatment地径/mm
Ground diameter总生物量/g
Total biomass根部产量/g
Root yieldCK 5.19±0.35a 6.12±0.31d 0.76±0.04c T1 5.25±0.14a 8.00±0.14c 0.89±0.05c T2 6.32±0.54a 9.65±0.61bc 1.11±0.08b T3 6.07±0.31a 13.96±0.78a 1.63±0.04a T4 6.22±0.44a 10.88±0.24b 1.22±0.08b T5 5.72±0.25a 10.40±0.82b 1.26±0.05b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among different treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 2 复合肥不同施用量的肥料利用效率1)
Table 2 Fertilizer utilization efficiency of different application amounts of compound fertilizer
% 处理 Treatment N P K T1 36.97±1.37a 2.69±0.31ab 9.16±0.53b T2 27.14±4.41b 2.31±0.07b 8.04±2.05b T3 41.35±2.37a 3.23±0.20a 16.62±2.63a T4 20.97±1.29bc 2.11±0.17b 7.27±1.00b T5 18.38±1.20c 1.18±0.10c 6.11±0.56b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 3 主成分分析的载荷系数
Table 3 Loading coefficients for principal component analysis
评价指标
Evaluation
index第一主成分
The first
principal
component第二主成分
The second
principal
component总生物量 Total biomass 0.952 −0.083 产量 Yield 0.959 −0.154 全株N含量
Whole plant N content0.843 0.353 全株P含量
Whole plant P content0.708 0.475 全株K含量
Whole plant K content0.961 −0.141 蔗糖积累量
Sucrose accumulation0.946 −0.189 1−蔗果三糖积累量
1-Kestose accumulation0.935 −0.232 耐斯糖积累量
Nystose accumulation0.755 −0.372 1F−果呋喃糖基耐斯糖累积量
1F-Fructofuranosyl nystose accumulation0.618 0.655 特征值 Eigenvalue 6.682 1.057 方差贡献率/%
Contribution of variance74.246 11.748 累计贡献率/%
Cumulative contribution rate74.246 85.994
下载: 导出CSV
表 4 复合肥不同施用量的主成分分析综合得分
Table 4 Comprehensive score of principal component analysis of different application amounts of compound fertilizer
处理
Treatment主成分分值
Principal component value综合得分(F)
Comprehensive
score排序
RankingF1 F2 CK −3.92 −0.49 −2.96 6 T1 −2.28 0.27 −1.66 5 T2 −0.01 −0.86 −0.11 4 T3 3.48 −0.52 2.53 1 T4 1.07 1.17 0.93 3 T5 1.65 0.42 1.27 2
下载: 导出CSV
-
[1] 国家药典委员会. 中华人民共和国药典: 1部[M]. 北京: 中国医药科技出版社, 2020: 83. [2] YANG Z, HU J, ZHAO M. Isolation and quantitative determination of inulin-type oligosaccharides in roots of Morinda officinalis[J]. Carbohydrate Polymers, 2011, 83(4): 1997-2004. doi: 10.1016/j.carbpol.2010.11.006
[3] CHEN D, LI N, LIN L, et al. Confocal mirco-Raman spectroscopic analysis of the antioxidant protection mechanism of the oligosaccharides extracted from Morinda officinalis on human sperm DNA[J]. Journal of Ethnopharmacology, 2014, 153(1): 119-124. doi: 10.1016/j.jep.2014.01.021
[4] YANG X, CHEN D, CHEN T, et al. Oligosaccharides from Morinda officinalis slow the progress of aging mice by regulating the key microbiota-metabolite pairs[J]. Evidence-Based Complementary and Alternative Medicine, 2019, 2019: 9306834. doi: 10.1155/2019/9306834.
[5] ZHU J, PENG Q, XU Y, et al. Morinda officinalis oligosaccharides ameliorate depressive-like behaviors in poststroke rats through upregulating GLUT3 to improve synaptic activity[J]. FASEB Journal, 2020, 34(10): 13376-13395. doi: 10.1096/fj.201902546RR
[6] 葛阳, 康传志, 万修福, 等. 生产中氮肥施用及其对中药材产量和质量的影响[J]. 中国中药杂志, 2021, 46(8): 1883-1892. [7] 刘哲, 钟海蓉, 威则日沙, 等. 氮、磷、钾配施对华重楼根茎的产量和有效成分含量的影响[J]. 中草药, 2019, 50(24): 6103-6113. [8] 邹庆军, 汪涛, 郭巧生, 等. N、P、K调控对野菊药用部位性状及活性成分含量的影响[J]. 中草药, 2019, 50(2): 535-540. doi: 10.7501/j.issn.0253-2670.2019.02.036 [9] 苏文华, 张光飞, 王泽明, 等. 氮、磷和钾肥对灯盏花生长和有效成分积累的影响[J]. 中草药, 2009, 40(12): 1963-1966. [10] 苏现明, 王洪庆, 陈若芸, 等. 巴戟天属植物化学成分及药理活性研究进展[J]. 中药材, 2017, 40(4): 986-991. [11] LIANG J, LIANG J, HAO H, et al. The extracts of Morinda officinalis and its hairy roots attenuate dextran sodium sulfate-induced chronic ulcerative colitis in mice by regulating inflammation and lymphocyte apoptosis[J]. Frontiers in Immunology, 2017, 8: 905. doi: 10.3389/fimmu.2017.00905.
[12] 代文豪. 发酵法炮制巴戟天工艺优化及其活性成分变化研究[D]. 广州: 华南农业大学, 2018. [13] 李倩. 巴戟天炮制过程中化学成分变化规律的研究[D]. 广州: 广州中医药大学, 2015. [14] 朱露林, 沈燚, 易晓薇, 等. 不同产地及种植年限巴戟天中6种寡糖类成分的定量分析[J]. 中草药, 2021, 52(14): 4348-4352. [15] 冼丽铧, 吴道铭, 隆曼迪, 等. 土壤养分、种植年限和种植方式对巴戟天寡糖含量的影响[J]. 华南农业大学学报, 2021, 42(3): 75-85. [16] 邵玲, 李美映, 黎学荣, 等. 肇庆巴戟天两种种植模式的植物学性状及药效成分研究[J]. 中药材, 2019, 42(11): 2480-2485. [17] 隆曼迪. 种植环境及栽培管理措施对巴戟天生长和品质的影响[D]. 广州: 华南农业大学, 2019. [18] 吴礼树. 土壤肥料学[M]. 北京: 中国农业出版社, 2011. [19] 邵镪钎, 李丹, 蒋攀, 等. 氮、磷、钾配施对川明参产量和品质的影响[J]. 中草药, 2018, 49(16): 3926-3932. doi: 10.7501/j.issn.0253-2670.2018.16.029 [20] 梁琴, 陈兴福, 李瑶, 等. 化肥与有机肥配施对川芎产量的影响[J]. 中药材, 2015, 38(10): 2015-2020. [21] 陈晓玉, 贺超, 闫滨滨, 等. N、P、K配施对丹参产量和质量的影响[J]. 中草药, 2019, 50(3): 722-730. [22] 范巧佳, 刘灵, 郑顺林, 等. 春季施氮时期和数量对川芎产量、生物碱和阿魏酸的影响[J]. 四川农业大学学报, 2013, 31(2): 136-139. [23] 贾袭伟, 晋小军, 李慧, 等. 磷酸二氢钾施肥量对款冬生长及产量、质量的影响[J]. 中成药, 2022, 44(1): 314-317. [24] DIMKPA C O, FUGICE J, SINGH U, et al. Development of fertilizers for enhanced nitrogen use efficiency: Trends and perspectives[J]. Science of the Total Environment, 2020, 731: 139113. doi: 10.1016/j.scitotenv.2020.139113.
[25] 洪丕征, 刘世荣, 王晖, 等. 正常和环割条件下不同形态氮素添加对红椎幼苗光合特性的影响[J]. 生态学报, 2016, 36(14): 4485-4495. [26] 潘超美, 詹若挺, 丁平, 等. 不同有机肥对巴戟天生长及矿质养分吸收的影响[J]. 中药材, 2002, 25(10): 699-701. doi: 10.3321/j.issn:1001-4454.2002.10.002 [27] 漆小雪, 韦霄, 陈宗游, 等. 黄花蒿干物质的积累及青蒿素与N、P、K量的动态变化研究[J]. 中草药, 2011, 42(12): 2541-2544. [28] 董娟娥, 张康健, 梁宗锁. 植物次生代谢与调控[M]. 杨凌: 西北农林科技大学出版社, 2009. [29] 冯嘉仪, 谢姗宴, 吴道铭, 等. 氮磷钾配施对银杏果实和外种皮产量及品质的影响[J]. 生态学杂志, 2021, 40(6): 1650-1659. [30] 胡佳栋, 毛歌, 张志伟, 等. 不同施肥处理对党参产量和次生代谢物含量的影响研究[J]. 中国中药杂志, 2017, 42(15): 2946-2953. [31] 卢挺, 杨全, 唐晓敏, 等. 氮磷钾配比施肥对广金钱草产量及质量的影响[J]. 广西植物, 2014, 34(3): 426-430. [32] 吴向维. 基于巴戟天寡糖的喘可治注射剂质量评价研究[D]. 广州: 广州中医药大学, 2014. [33] 邓秋林, 杨正明, 陈雨, 等. 氮磷钾配施对瓦布贝母产量及总生物碱质量分数的影响[J]. 西北农业学报, 2019, 28(7): 1138-1146. [34] 鲁泽刚, 朱永全, 卢迎春, 等. 氮、磷、钾施用对灯盏花产量和主要提取物收获量的影响[J]. 植物科学学报, 2019, 37(1): 55-62. -
期刊类型引用(5)
1. 程功，焦洁洁，蒋建潮，薛逸琳，郑炳松，闫道良. 植物对氮磷亏缺的生长及生理、分子响应研究进展. 湖北林业科技. 2024(05): 65-72 . 
百度学术
2. 李雅雪，盘耀亮，彭光粉，田江，陆星，梁翠月. GmNTLs调控大豆根系响应低磷胁迫的功能研究. 华南农业大学学报. 2023(02): 221-229 . 本站查看
3. 薛迎斌，宋佳，李枭艺，李小豪，陈经烨，伍萍珍，朱胜男，刘颖. 大豆GmMADS4基因克隆、亚细胞定位及功能分析. 华南农业大学学报. 2023(03): 420-429 . 本站查看
4. 吕铭滔，龚海光，黄永芳，龚勇军. 磷对铝胁迫油茶芽苗初生根保护酶的影响. 经济林研究. 2023(03): 271-277 . 百度学术
5. 单玉静，王召军，程泽华，余婧，雷波，闫筱筱，张洪映，崔红. NtSPX4基因敲除对烟草苗期磷吸收和生长发育的影响研究. 中国烟草学报. 2022(01): 78-85 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 174
- HTML全文浏览量: 27
- PDF下载量: 17
- 被引次数: 7
