Functional characterization of phosphorus deficiency-responsive GmNTLs in soybean roots
-
摘要:目的
低磷和铝毒胁迫是酸性土壤中限制作物生产的重要因素。植物NTL转录因子参与调控多种环境胁迫(包括铝毒胁迫)的适应性机制，本文探究GmNTLs调控大豆Glycine max 根系响应低磷胁迫的功能。
方法通过RT-qPCR分析大豆15个GmNTLs基因在根系响应低磷胁迫的表达模式，进一步构建了GmNTL1/4/7/8/10/12共6个GmNTLs基因的拟南芥超量表达材料，探究GmNTL成员在拟南芥根系中响应低磷胁迫的功能。
结果系统进化树及组织表达模式分析结果表明，GmNTLs家族分3个亚族，各亚族成员在大豆中组织表达模式不同。RT-qPCR结果表明，低磷处理12 d显著提高了GmNTL1/4/7/8/10/12在大豆根系中的表达。在拟南芥中超量表达不同GmNTL基因对低磷的响应不同。高磷处理下，超量表达GmNTL4/10/12拟南芥的鲜质量显著增加；低磷处理时，超量表达GmNTL4显著提高拟南芥鲜质量，而超量表达GmNTL1/12拟南芥的鲜质量显著降低。同时，仅超量表达GmNTL12拟南芥的主根长显著缩短，而超量表达其他基因对拟南芥植株的主根长无明显影响。
结论GmNTLs参与大豆根系对低磷胁迫的响应，该结果可为培育磷高效型大豆品种提供数据支持。
Abstract:ObjectiveLow phosphorus (P) availability and aluminum (Al) toxicity constrain crop production in acid soils. NTL transcription factors play an important role in the mechanisms of plant response to various abiotic stresses, including Al toxicity. This study focused on analyzing the function of GmNTLs in soybean roots responding to P deficiency.
MethodExpression pattern were performed on 15 members of the GmNTL family in the soybean by RT-qPCR assays. We further investigated the function of some GmNTLs members in adaptation to low P by overexpressing the genes of GmNTL1/4/7/8/10/12 in Arabidopsis thaliana.
ResultThe phylogenetic analysis and the tissue expression analysis of each subfamily members revealed that 15 GmNTLs were divided into three subgroups and different GmNTL family members had different tissue expression pattern in soybean. The RT-qPCR results showed that the expression levels of GmNTL1/4/7/8/10/12 in soybean roots significantly increased after 12 days of low P treatment. Overexpression of different GmNTLs in Arabidopsis showed different responses to low P. The fresh weight of transgenic Arabidopsis overexpressing GmNTL4/10/12 significantly increased compared to control lines under high P treatment. Overexpression of GmNTL4 significantly improved fresh weight of transgenic Arabidopsis plant under low P deficiency; Whereas, the plant fresh weight of transgene lines overexpressing GmNTL1/12 significantly decreased. Overexpression of GmNTL12 reduced the primary root length of transgenic plants; Whereas, overexpression of other GmNTLs had no significant effect on primary root length.
ConclusionGmNTLs involve in the response of soybean roots to low P stress, and these results can provide a theoretical basis for cultivating soybean varieties with high P efficiency.
-
Keywords:
- Glycine max /
- Root /
- Low phosphorus stress /
- GmNTLs /
- NTL transcription factor
-
随着经济全球化的加深，我国进口货物和进口国家数量不断增长，同时，许多有害杂草伴随着国际流通进入我国[1]。研究表明，进口额每增长1个百分点，植物疫情在海关截获量约增长2个百分点[2]。近年来有多起南方三棘果Emex australis在海关截获的报告。例如，2012—2015年广东黄埔口岸在主要进境粮食中截获的有害生物中，南方三棘果截获高达90批次排第7位。此外在广东的肇庆和汕头口岸以及江苏口岸也都有截获报道[3-4]。
南方三棘果是蓼科亦模属的冬型一年生杂草，喜生长在湿润和半湿润热带、亚热带和温带气候区，被我国列入进口植物检疫对象。目前在我国内地尚未发现，仅在台湾有分布记录。南方三棘果单株植株可以产生1000粒以上的果实，种子在土壤中休眠多年，仍有繁殖能力，且具有繁殖速度快、传播迅速、危害性大和防治难度高的特点，一旦大面积入侵我国，难以根治，将会给农业生产造成巨大破坏[5-7]。
随着信息化技术的发展，目前用于预测物种潜在分布区的模型有CLIMEX、GARP、DIVA-GIS、BIOCLIM和MaxEnt等。相比于其他模型，MaxEnt生态位模型在分布数据较少的情况下，预测精确度和准确度高[8-9]。其通过将研究区域的环境因子转换为栅格数据，形成多个栅格单元，把物种已知分布区域的栅格单位作为样点，然后以分布点区域的环境变量数据对所有栅格单元进行数理统计，得到当前和未来一段时间内研究区域的每一个栅格单元分布的概率以及环境变量数值与分布概率之间的关系[10-12]。
由于各物种的研究和关注程度不同，物种在世界范围内详实的分布记录存在差异，因此物种的分布数据获得量存在差异[12]。南方三棘果属于数据获得量较少的物种，因此选择MaxEnt作为预测南方三棘果在我国潜在适生区的模型。全球气温呈现升高趋势，会影响植物的潜在分布区域变化[13-14]。气候系统模式能够很好地反映气候变化的趋势，对未来一段时间的气候变化进行准确预测[15-17]。因此本文采用国家气候中心4种气候情景下的未来气候数据对南方三棘果在我国的适生区进行分析预测。
本研究将MaxEnt模型与ArcGIS软件相结合，得到南方三棘果潜在分布区的可视化地图，以期能够为有关部门在潜在分布区的相关地区开展入侵风险评估、检验检疫和监测预警等提供理论参考，对防控南方三棘果传入、定殖和扩散，保障农业生产和生物多样性健康以及生态安全和经济贸易安全等具有重要意义[18]。
1. 材料与方法
1.1 技术路线
技术路线包括查找和处理数据、模型优化与构建以及制作可视化地图3部分，按照“物种分布点查找、处理→查找影响分布的环境因子以及SPSS分析→MaxEnt模型参数优化→构建MaxEnt模型→ArcGIS制作可视化图”的研究方案进行，数据分析贯穿于整个研究过程。
1.2 相关软件和程序
MaxEnt软件3.4.1版、ArcGIS软件10.4.1版、SPSS Statistics软件20.0版、R语言3.6.3版、ENMeval程序包以及Microsoft Excel 2016。
1.3 环境数据及处理
从世界气候数据网站( https://www.worldclim.org/)下载2.1版、分辨率为10 min(340 km2)的1970—2000年和国家气候中心气候系统模式(BCC-CSM2-MR)4个气候情景模式(SSP1-2.6、SSP 2-4.5、SSP 3-7.0和SSP5-8.5)2021—2040年世界范围内19个生物气候变量数据(Bio1~Bio19)以及高程海拔数据(Elev)。采用ArcGIS软件从高程海拔数据中提取出坡度(Slo)和坡向(Asp)[19-20]。一共得到温度、降水和海拔特征的环境变量数据22个(表1)。
表 1 环境变量因子Table 1. Environmental variable factors变量 Variable 描述 Description 变量 Variable 描述 Description Bio1 年平均温度/℃ Annual mean temperature Bio12 年均降水量/mm Annual mean precipitation Bio2 昼夜温差月均值/℃ Monthly mean diurnal temperature range Bio13 最湿月降水量/mm Precipitation of the wettest month Bio3 等温性/% (Bio2/Bio7) ×100 Isothermality Bio14 最干月降水量/mm Precipitation of the driest month Bio4 温度季节变化标准差/℃ Standard deviation of seasonal variation of temperature Bio15 降水量变异系数 Standard deviation of the precipitation Bio5 最暖月最高温/℃ Maximum temperature of the warmest month Bio16 最湿季度降水量/mm Precipitation of the wettest quarter Bio6 最冷月最低温/℃ Minimum temperature of the coldest month Bio17 最干季度降水量/mm Precipitation of the driest quarter Bio7 年均温度变化范围/℃ Temperature annual range Bio18 最暖季度降水量/mm Precipitation of the warmest quarter Bio8 最湿季度平均温度/℃ Mean temperature of the wettest quarter Bio19 最冷季度降水量/mm Precipitation of the coldest quarter Bio9 最干季度平均温度/℃ Mean temperature of the driest quarter Elev 高程海拔 Elevation Bio10 最暖季度平均温度/℃ Mean temperature of the warmest quarter Slo 坡度/° Slope Bio11 最冷季度平均温度/℃ Mean temperature of the coldest quarter Asp 坡向 Aspect 参考文献[21-22]方法，用ArcGIS将所有环境变量坐标系及基准面统一定义，并转出为ASCⅡ格式，以待MaxEnt模型预测分布；使用R语言将TIFF格式的环境变量因子统一转换为ASC格式，以待使用R语言ENMeval 程序包对MaxEnt模型进行优化。为了避免变量间的高度共线性导致的过拟合和减少运算量，对环境数据进行相关性分析。将物种分布点与环境变量输入MaxEnt模型运算10次，获得各个环境变量的贡献率，将环境数据和物种分布点数据在ArcGIS处理为“*．txt”格式文件，将文件利用Excel转化为“*．csv”格式文件，最后将分析数据导入SPSS进行双变量相关性分析，作双尾检验。对相关系数(R)>0.8的环境变量，选择贡献率较高的作为主要变量进行下一步分析，舍去贡献率较低项Bio5、Bio6、Bio7、Bio10、Bio12、Bio16和Bio17。舍去贡献率为0的环境变量Bio8。由于南方三棘果为冬型一年生草本植物，在最湿月和最暖季度是种子形态[23-24]，因此舍去Bio13和Bio18。最终得到12个环境变量因子(表1)，分别为Bio1、Bio2、Bio3、Bio4、Bio9、Bio11、Bio14、Bio15、Bio19、Elev、Slo和Asp。
1.4 南方三棘果分布数据及处理
在全球生物多样性信息交换平台( https://www. gbif.org/)等平台数据库检索物种分布相关记录，通过发表的文献、专著、文件、报告等查询相关分布信息及记录[25-27]。一共得到1473个分布数据点，去除重复、无效、缺失、错误分布数据后，得到1205个分布点数据。利用下载好的气候数据用ENMeval程序包对南方三棘果分布点进行过拟合分析，最终得到南方三棘果自然分布点647个。
1.5 地图数据
从标准地图服务系统( http://bzdt.ch.mnr.gov.cn/)下载GS(2022)4312号标准地图。
1.6 模型优化
使用R语言ENMeval程序包将南方三棘果分布点分为2部分，75%随机样本数据用于训练集，25%用于测试集，对MaxEnt软件中调控倍频(Regularization multiplier，RM)和特征组合(Feature combination，FC)参数进行优化，利用数据包设定RM从0.1到4，每次增加0.1，共40个调控倍频，同时设定29个特征组合，即H、L 、LH、LP、LPH、LPT、LPTH、LQ、LQH、LQP、LQPH、LQPT、LQPTH、LQT、LQTH、LT、P、PH、PT、PTH、Q、QH、QP、QPH、QPT、QT、QTH、T、TH，其中L为线性(Linear)，Q为二次型(Quadratic)，H为片段化(Hinge)，P为乘积型(Product)，T为阈值型(Threshold)，ENMeval 语言包将上述 1160 种参数组合进行检验，最终采用最小信息准则(The minimum information criterion AICc value)模型评估受试者工作特征曲线(Receiver operating curve，ROC)下的面积(Area under curve，AUC)差异检验模型的拟合度与复杂度，当AICc值最低(AICc=0)时模型参数组合最佳[28]，用于MaxEnt模型建模。
1.7 构建 MaxEnt 预测模型
参考前人所述方法[29-33]，将处理后的当前环境变量和物种分布点数据输入MaxEnt模型，模型参数设置随机选择25%的样本作为测试集，75%的样本数据为训练集，同时勾选制作响应曲线功能以获得分布概率与气候因子的关系，勾选刀切法检验环境变量的重要程度，RM与FC取优化后的参数，其他参数为默认设置。以AUC评价模型的准确性。AUC越大，表明模型的可靠性越高，0.9以上表示模型结果优秀可靠，为降低模型的随机性和模拟结果的不确定性，模型重复运行10次，以运行的均值作为模型运算结果。
1.8 绘制物种分布图
将MaxEnt模型生成的ASC格式物种分布数据导入ArcGIS软件，设置坐标系为地理坐标系GCS_WGS_1984，通过掩膜提取获得南方三棘果在我国的分布图，对预测结果利用自然间断点分级法进行重分类，得到南方三棘果的非适生区、低适生区、中适生区和高适生区。
2. 结果与分析
2.1 MaxEnt生态位模型优化结果
基于12个环境变量因子和南方三棘果的647个分布点数据，MaxEnt软件初始默认RM=1，FC为LQPH，AICc=142.939251，经ENMeval语言包优化后RM=1，FC为LQPT，AICc=0，表明其为最优参数组合。在优化参数下进行模型预测，采用ROC分析法进行模型准确度检验。经10次重复训练得到的AUC > 0.9( 图1)，模拟结果优秀可靠。
2.2 影响分布环境变量因子贡献率及适宜范围
使用12个环境因子建立模型，模型参数为ENMeval语言包优化结果，RM=1，FC为LQPT，采用刀切法评估各环境因子对南方三棘果地理分布的重要性。贡献率超过15%的包括最冷季度平均气温(Bio11，27.7%)、等温性(Bio3，24.9%)、温度季节变化标准差(Bio4，18.2%)，这3个因子的贡献率合计达70.8%。贡献率低于15%的包括年平均温度(Bio1，14.0%)、最干月降水量(Bio14，6.5%)、降水量变异系数(Bio15，6.5%)、最冷季度降水量(Bio19，2.7%)、高程海拔(Elev，1.4%)、坡度(Slo，0.6%)、最干季度平均温度(Bio9，0.3%)、昼夜温差月均值(Bio2，0.2%)、坡向(Asp，0.2%)。
刀切法检验表明，仅使用单一环境变量时，正则化训练增益和测试增益都表明最重要的环境因子为最冷季度平均温度(Bio11)(图2)。
![]()
图 2 环境变量重要性的刀切法检验结果Asp: 坡向；Slo: 坡度；Elev: 高程海拔；Bio9: 最干季度平均温度；Bio4: 温度季节变化标准差；Bio3: 等温性 ( Bio2/ Bio7) ×100；Bio2: 昼夜温差月均值；Bio19: 最冷季度降水量；Bio15: 降水量变异系数；Bio14: 最干月降水量；Bio11: 最冷季度平均温度；Bio1: 年平均温度Figure 2. Results of Jackknife for the importance of the environmental variablesAsp: Aspect; Slo: Slope; Elev: Elevation; Bio9: Mean temperature of the driest quarter; Bio4: Standard deviation of seasonal variation of temperature; Bio3: Isothermality; Bio2: Mean diurnal range; Bio19: Precipitation of the coldest quarter; Bio15: Precipitation seasonality; Bio14: Precipitation of the driest month; Bio11: Mean temperature of the coldest quarter; Bio1: Annual mean temperature一般认为存在概率≥0.5对应的环境因子范围是物种适宜分布的阈值。根据响应曲线确定影响南方三棘果分布环境因子的阈值：最冷季度平均温度9.35~12.76 ℃、等温性45.2%~53.72%以及温度季节变化标准差337.06~616.40 ℃(图3)。
2.3 南方三棘果在当前气候下的潜在适生分布区
将南方三棘果分布点和环境变量数据导入MaxEnt模型进行模拟，得到南方三棘果在当前(2022年)气候下的潜在分布区(图4)。总适生面积约69.13万km2，高适生区、中适生区、低适生区分别为11.64万、13.37万、44.02万km2。南方三棘果的适生范围主要集中在我国云南、广西、广东、福建以及西藏局部，其中，高适生区和中适生区主要位于云南、海南以及西藏局部，低适生区主要位于广西、广东、台湾、福建和四川。
![]()
图 4 南方三棘果在中国的适生区(2022年)Figure 4. Current suitable regions of Emex australis in China (2022)2.4 南方三棘果在未来气候模型下的潜在适生分布区
如表4所示，SSP1-2.6气候情景下，2040年南方三棘果总适生区面积约63.54万km2，较2022年分布面积减少约5.48万km2，其中，低适生区面积减少5.93万km2，中适生区面积增加 6.57万km2，高适生区面积减少6.12万km2；SSP2-4.5气候情景下，2040年南方三棘果总适生区面积约65.10万km2，较2022年分布面积减少约3.92万km2，其中，低适生区面积减少2.31万km2，中适生区面积增加2.31万km2，高适生区面积减少3.92万km2；SSP3-7.0气候情景下， 2040年南方三棘果总适生区面积约63.34万km2，较2022年分布面积减少约5.68万km2，其中，低适生区面积减少9.74万km2，中适生区面积增加6.51万km2，高适生区面积减少2.45万km2；SSP5-8.5情景下，2040年南方三棘果总适生区面积约62.90万km2，较2022年分布面积减少约6.12万km2，其中，低适生区面积减少5.42万km2，中适生区面积增加3.39万km2，高适生区面积减少4.09万km2。
表 4 不同时期气候情景下南方三棘果在中国潜在适生区的面积Table 4. Areas of the potential suitable region for Emex australis in China under different climate scenarios in different periods×104 km2 年份 Year 气候情景 Climate scenario 低适生区 Low suitable region 中适生区 Middle suitable region 高适生区 High suitable region 总适生区 Total suitable region 2022 44.02 13.36 11.63 69.02 2040 SSP1-2.6 38.09 19.93 5.51 63.54 2040 SSP2-4.5 41.71 15.67 7.71 65.10 2040 SSP3-7.0 34.27 19.88 9.18 63.34 2040 SSP5-8.5 38.59 16.76 7.54 62.90 如图5所示，对不同气候情景下南方三棘果适生区的空间格局进行对比分析，结果表明，2040年，未来4种情景模式下的适生区都呈现轻微缩小趋势，广东和广西的部分低适生区变为非适生区，云南的高适生区变为中适生区，海南的高、中、低适生区都呈减少趋势。从单个情景模式下看也有一些地区变得更加适应南方三棘果分布，SSP1-2.6气候情景模式下西藏出现局部高适生区，四川的局部低适生区转变为中适生区；SSP2-4.5气候情景模式下福建的一些非适生区变为低适生区；SSP3-7.0气候情景模式下四川部分低适生区转变为中适生区，局部地区变为高适生区，西藏出现局部高适生区；SSP5-8.5气候情景模式下，四川和西藏的一些非适生区变为低适生区，其中，西藏出现局部高适生区。这些变化区域为南方三棘果分布的敏感区。
![]()
图 5 4种气候情景模式下南方三棘果在中国的适生区Figure 5. The suitable growth reqion of Emex australis in China under four climate scenarios3. 结论与讨论
本模型的AUC测试值 > 0.9，模型结果优秀可靠。最冷季度平均气温(Bio11)、等温性(Bio3)和温度季节变化标准差(Bio4)是影响南方三棘果地理分布的主要环境因子。正则化训练增益和测试增益以及贡献率都表明最冷季度平均气温(Bio11)对南方三棘果分布影响最大，贡献率为27.7%。环境因子响应曲线表明，当最冷季度平均气温为9.35~12.76 ℃时，南方三棘果发生概率≥0.5。南方三棘果的适宜分布区总面积约为69.02万km 2，适生区域主要集中于云南、广东、广西和福建地区，其中高、中适生区主要集中于云南。
通过对2022—2040年国家气候中心气候系统模式(BCC-CSM2-MR)的4种未来气候情景下的南方三棘果潜在分布区研究发现，总体上适生区呈现轻微减少趋势。四川、西藏、广东和广西属于南方三棘果分布的敏感区，易随气候变化而导致适生区等级以及适生面积发生变化。广东、广西、海南和云南的一些适生区在未来一段时间会出现转变为非适生区或降低适生等级的情况。福建部分地区在SSP2-4.5情景模式下会更适宜南方三棘果生存，西藏和四川部分地区在SSP1-2.6、SSP3-7.0和SSP5-8.5情景下会更适宜南方三棘果生存。
其中，西藏的适生区尤为特殊，主要集中在藏南地区的错那县和墨脱县以及朗县最南部。这一部分地区由于地形原因导致气候复杂，高原气候和亚热带、热带季风气候同时存在，加上地处边疆，适生区西边与不丹接壤，南边与印度接壤，导致藏南外来物种入侵严重，值得相关部门进一步关注。中国科学院植物研究所的植物物种信息系统( http://www.iplant.cn/)显示三叶鬼针草Bidens pilosa、胜红蓟Ageratum conyzoides等原产于热带地区的外来杂草在墨脱县都有分布记录[34]。国家气候中心气候系统模式数据显示未来这一地区的温度季节变化标准差(Bio4)为390~550 ℃，等温性(Bio3)为39%~45%，最冷季度平均温度(Bio11)为0.9~19 ℃，这种气候适宜南方三棘果分布，因此，对这里需要加大监测力度，阻止南方三棘果进入藏南地区，并防止由藏南地区入侵我国西南和华南地区。
南方三棘果是具有严重危害的检疫杂草，对农业生产会造成巨大破坏[35-36]。云南是我国葡萄十大产地之一，南方三棘果在我国的高、中适生区主要集中于云南，南方三棘果一旦定殖于葡萄园，可能会有部分种子在机械化收割时与葡萄一起收割，影响葡萄的品相和食品安全性[37]。藏南是西藏的主要农区，藏南地区生产的小麦是西藏的主要粮食作物之一，研究表明南方三棘果茎的水提取物和植株残留物会抑制小麦种子萌发和早期幼苗的生长，降低小麦产量[38-39]。油菜是广西、广东和福建的主要油料作物，在农田里南方三棘果与油料作物存在种间竞争关系[40]。
外来入侵杂草传入与人员和货物的流动密切相关，作为国与国之间管控人员和货物流通的重要节点，海关和动植物检验检疫局实施的植物检疫是防止外来杂草传入的最重要方式。在适生区域内以及邻近适生区域的海关和动植物检验检疫局应当加大对南方三棘果的检疫力度，防止南方三棘果的传入[41]。尤其值得注意的是，广东、广西和福建的海关每年都有大量的货物进口，加上这些地区存在大面积的南方三棘果适生区，海关和动植物检验检疫局的植物疫情数据也表明在广东和福建的海关有多起南方三棘果的截获记录，因此，这些地区是南方三棘果传入的高风险地区。
有害杂草入侵的定殖和扩散阶段在适生区范围内进行，其中，外部气候环境的变化是影响其扩散速度的重要因素，试验表明最冷季度平均温度是影响南方三棘果分布的最重要环境因子，在适生区气候适宜其生存的重要年份要加大南方三棘果的监测，有效防止其定殖、扩散，以免影响农作物的质量和产量[42]。生物入侵杂草的分布与气候环境、地形因素、土壤类型以及植被类型等多种因素相关，在使用MaxEnt生态位模型预测入侵杂草潜在分布区时，目前尚未有统一的影响因子选择。朱煜建等[43]在预测疣果匙荠Bunias orientalis适生区时选择了气候环境、地形因素以及土壤类型的环境变量因子。郭燕青等[44]在预测假臭草Praxelis clematidea时仅选择了气候因素有关因子。本研究选用与环境和地形有关的变量因子，缺少土壤类型、物种类型等与南方三棘果分布有关的其他变量，这可能对适生区的预测带来一定的偏差。伴随着我国与世界各国日益频繁的贸易往来，生物入侵难以避免。近些年假高粱Sorghum halepense和豚草Ambrosia artemisiifolia等入侵生物对我国农业生产和生态安全造成了巨大破坏。采用MaxEnt生态位模型将信息学技术与生物的生态学特性结合，对各种入侵杂草的潜在分布区进行预测，对我国的口岸检疫和入侵杂草的防控具有指导意义。
-
![]()
图 1 大豆GmNTL蛋白系统进化树及组织表达模式分析
1：新叶；2：花；3：根；4：根瘤；5：花后14 d豆荚壳；6：花后14 d种子
Figure 1. Phylogenetic tree and tissue expression analysis of GmNTL proteins in soybean
1: Young leaf; 2: Flower; 3: Root; 4: Nodule; 5: Pod shell after flowering for 14 days; 6: Seed after flowering for 14 days
![]()
图 2 不同磷处理时间对大豆生长的影响
Figure 2. Effect of different phosphorus treatment time on soybean growth
![]()
图 3 不同磷处理时间对大豆生物量和根长的影响
“*”“**”“***”分别表示相同胁迫时间不同磷浓度处理在P < 0.05、0.01、0.001水平差异显著(t检验)
Figure 3. Effect of different phosphorus treatment time on soybean biomass and root length
“*” “**” “***” indicate significant differences between different phosphorus concentration treatments at the same stress time at P < 0.05, 0.01 and 0.001 levels, respectively (t test)
![]()
图 4 不同磷处理时间对大豆根系GmNTLs表达模式的影响
“*”“**”“***”分别表示相同胁迫时间不同磷浓度处理在P < 0.05、0.01、0.001水平差异显著(t检验)
Figure 4. Effect of different phosphorus treatment time on relative expression of GmNTLs in soybean roots
“*” “**” “***” indicate significant differences between different phosphorus concentration treatments at the same stress time at P < 0.05, 0.01 and 0.001 levels, respectively (t test)
![]()
图 5 不同磷处理浓度下大豆根系GmNTLs相对表达量
各小图中，柱子上方的不同小写字母表示不同磷处理浓度间差异显著(P < 0.05，Duncan’s法)
Figure 5. Relative expression of GmNTLs in soybean roots under different phosphorus concentrations
In each figure, different lowercase letters on bars indicate significant differences among different phosphorus concentrations (P < 0.05, Duncan’s method)
![]()
图 6 不同磷浓度下转基因拟南芥生长情况
WT：拟南芥Col-0野生型；OX1和OX2：超量表达株系；标尺=1 cm
Figure 6. Growth of transgenic Arabidopsis under different phosphorus concentrations
WT: Arabidopsis Col-0 wild-type; OX1 and OX2: Overexpressed lines; Bar = 1 cm
![]()
图 7 不同磷浓度下转基因拟南芥的生物量
WT：拟南芥Col-0野生型，OX1和OX2：超量表达株系；“*”“**”“***”分别表示相同磷浓度下超量表达株系和野生型在P < 0.05、0.01、0.001水平差异显著（t检验）
Figure 7. Biomass of transgenic Arabidopsis under different phosphorus concentrations
WT: Arabidopsis Col-0 wild-type, OX1 and OX2: Overexpressed lines; “*” “**” “***” indicate significant differences between overexpressed lines and wild-type under the same phosphorus concentration at P < 0.05, 0.01 and 0.001 levels, respectively (t test)
![]()
图 8 不同磷浓度下转基因拟南芥的主根长
WT：拟南芥Col-0野生型，OX1和OX2：超量表达株系；“*”“**”“***”分别表示相同磷浓度下超量表达株系和野生型在P < 0.05、0.01、0.001水平差异显著(t检验)
Figure 8. Primary root length of transgenic Arabidopsis under different phosphorus concentrations
WT: Arabidopsis Col-0 wild-type, OX1 and OX2: Overexpressed lines; “*” “**” “***” indicate significant differences between overexpressed lines and wild-type under the same phosphorus treatment concentration at P < 0.05, 0.01 and 0.001 levels, respectively (t test)
-
[1] OOKA H, SATOH K, DOI K, et al. Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana[J]. DNA Research, 2003, 10(6): 239-247. doi: 10.1093/dnares/10.6.239
[2] KIM S Y, KIM S G, KIM Y S, et al. Exploring membrane-associated NAC transcription factors in Arabidopsis: Implications for membrane biology in genome regulation[J]. Nucleic Acids Research, 2007, 35(1): 203-213. doi: 10.1093/nar/gkl1068
[3] 赵翠珠, 刘振华, 赵赫, 等. 植物NAC膜结合转录因子的研究进展[J]. 生命科学, 2012, 24(1): 74-80. [4] CHI Y H, MELENCION S, ALINAPON C V, et al. The membrane-tethered NAC transcription factor, AtNTL7, contributes to ER-stress resistance in Arabidopsis[J]. Biochemical and Biophysical Research Communications, 2017, 488(4): 641-647. doi: 10.1016/j.bbrc.2017.01.047
[5] ZHAO J, LIU J S, MENG F N, et al. ANAC005 is a membrane-associated transcription factor and regulates vascular development in Arabidopsis[J]. Journal of Integrative Plant Biology, 2016, 58(5): 442-451. doi: 10.1111/jipb.12379
[6] WANG D, YU Y, LIU Z, et al. Membrane-bound NAC transcription factors in maize and their contribution to the oxidative stress response[J]. Plant Science, 2016, 250: 30-39. doi: 10.1016/j.plantsci.2016.05.019
[7] KIM M J, PARK M, SEO P J, et al. Controlled nuclear import of the transcription factor NTL6 reveals a cytoplasmic role of SnRK2.8 in the drought-stress response[J]. Biochemical Journal, 2012, 448(3): 353-363. doi: 10.1042/BJ20120244
[8] KIM S, LEE A, YOON H, et al. A membrane-bound NAC transcription factor NTL8 regulates gibberellic acid-mediated salt signaling in Arabidopsis seed germination[J]. Plant Journal, 2008, 55: 77-88. doi: 10.1111/j.1365-313X.2008.03493.x
[9] SEO P J, KIM M J, PARK J, et al. Cold activation of a plasma membrane-tethered NAC transcription factor induces a pathogen resistance response in Arabidopsis[J]. Plant Journal, 2010, 61(4): 661-671. doi: 10.1111/j.1365-313X.2009.04091.x
[10] KOCHIAN L V, HOEKENGA O A, PINEROS M A. How do crop plants tolerate acid soils? Mechanisms of aluminum tolerance and phosphorous efficiency[J]. Annual Review of Plant Biology, 2004, 55: 459-493. doi: 10.1146/annurev.arplant.55.031903.141655
[11] 吴佩, 李浩, 早浩龙, 等. 植物对缺磷和铝毒协同进化应答的分子生理机制[J]. 生物技术通报, 2020, 36(7): 170-181. [12] LIAO H, WAN H, SHAFF J, et al. Phosphorus and aluminum interactions in soybean in relation to aluminum tolerance, exudation of specific organic acids from different regions of the intact root system[J]. Plant Physiology, 2006, 141(2): 674-684. doi: 10.1104/pp.105.076497
[13] LIANG C, PINEROS M A, TIAN J, et al. Low pH, aluminum, and phosphorus coordinately regulate malate exudation through GmALMT1 to improve soybean adaptation to acid soils[J]. Plant Physiology, 2013, 161(3): 1347-1361. doi: 10.1104/pp.112.208934
[14] CHEN W, TANG L, WANG J, et al. Research advances in the mutual mechanisms regulating response of plant roots to phosphate deficiency and aluminum toxicity[J]. International Journal of Molecular Sciences, 2022, 23(3): 1137. doi: 10.3390/ijms23031137.
[15] 于芮, 王斌, 王建忠, 等. 我国大豆市场发展现状及建议[J]. 合作经济与科技, 2021(19): 92-93. doi: 10.3969/j.issn.1672-190X.2021.19.037 [16] 刘国选, 陈康, 陆星, 等. 大豆GmPIN2b调控根系响应低磷胁迫的功能研究[J]. 华南农业大学学报, 2021, 42(4): 33-41. doi: 10.7671/j.issn.1001-411X.202010014 [17] LI S, WANG N, JI D, et al. Evolutionary and functional analysis of membrane-bound NAC transcription factor genes in soybean[J]. Plant Physiology, 2016, 172(3): 1804-1820. doi: 10.1104/pp.16.01132
[18] 向凤宁, 王楠. 大豆威廉姆斯82中NAC膜结合转录因子基因GmNTL1及其应用: CN106167801A[P]. 2018-08-28[2022-04-16]. [19] 向凤宁, 王楠. 大豆威廉姆斯82中NAC膜结合转录因子基因GmNTL7及其应用: CN106367423B[P]. 2017-02-01[2022-04-16]. [20] LIN Y, LIU G, XUE Y, et al. Functional characterization of aluminum (Al)-responsive membrane-bound NAC transcription factors in soybean roots[J]. International Journal of Molecular Sciences, 2021, 22(23): 12854. doi: 10.3390/ijms222312854.
[21] PARK J, KIM Y S, KIM S G, et al. Integration of auxin and salt signals by the NAC transcription factor NTM2 during seed germination in Arabidopsis[J]. Plant Physiology, 2011, 156(2): 537-549. doi: 10.1104/pp.111.177071
[22] YOON H K, KIM S G, KIM S Y, et al. Regulation of leaf senescence by NTL9-mediated osmotic stress signaling in Arabidopsis[J]. Molecules and Cells, 2008, 25(3): 438-445.
[23] DUAN M, ZHANG R, ZHU F, et al. A lipid-anchored NAC transcription factor is translocated into the nucleus and activates glyoxalase I expression during drought stress[J]. Plant Cell, 2017, 29(7): 1748-1772. doi: 10.1105/tpc.17.00044
[24] PERET B, CLEMENT M, NUSSAUME L, et al. Root developmental adaptation to phosphate starvation: Better safe than sorry[J]. Trends in Plant Science, 2011, 16(8): 442-450. doi: 10.1016/j.tplants.2011.05.006
[25] PÉREZ-TORRES C, LÓPEZ-BUCIO J, CRUZ-RAMÍREZ A, et al. Phosphate availability alters lateral root development in Arabidopsis by modulating auxin sensitivity via a mechanism involving the TIR1 auxin receptor[J]. Plant Cell, 2008, 20(12): 3258-3272.
[26] SHEN C, WANG S, ZHANG S, et al. OsARF16, a transcription factor, is required for auxin and phosphate starvation response in rice (Oryza sativa L. )[J]. Plant Cell and Environment, 2013, 36(3): 607-620. doi: 10.1111/pce.12001
[27] BALZERGUE C, DARTEVELLE T, GODON C, et al. Low phosphate activates STOP1-ALMT1 to rapidly inhibit root cell elongation[J]. Nature Communications, 2017, 8: 15300. doi: 10.1038/ncomms15300.
[28] HAM B K, CHEN J, YAN Y, et al. Insights into plant phosphate sensing and signaling[J]. Current Opinion in Biotechnology, 2018, 49: 1-9. doi: 10.1016/j.copbio.2017.07.005
[29] CHEN Y N, SLABAUGH E, BRANDIZZI F. Membrane-tethered transcription factors in Arabidopsis thaliana: Novel regulators in stress response and development[J]. Current Opinion in Plant Biology, 2008, 11(6): 695-701. doi: 10.1016/j.pbi.2008.10.005
[30] LEE S, LEE H J, HUH S U, et al. The Arabidopsis NAC transcription factor NTL4 participates in a positive feedback loop that induces programmed cell death under heat stress conditions[J]. Plant Science, 2014, 227: 76-83. doi: 10.1016/j.plantsci.2014.07.003
-
期刊类型引用(2)
1. 徐慧娟，邝智祥，左珂菁，吴承洋，郑泓昊，谢青梅，陈伟国. tva受体基因起始密码子突变对清远麻鸡感染K亚群禽白血病病毒的影响. 华南农业大学学报. 2025(03): 351-357 . 
本站查看
2. 郑文静，田亚东，李转见，李东华. 禽白血病抗病育种研究进展. 中国畜禽种业. 2024(01): 42-46 . 百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 556
- HTML全文浏览量: 29
- PDF下载量: 382
- 被引次数: 2
