摘要
土地盐渍化对大豆的产量和品质产生了极大的负面影响,培育耐盐大豆品种是改善和提高盐胁迫下大豆产量和品质的有效途径之一。ERF转录因子对植物响应逆境胁迫十分重要,但在大豆中相关研究报道较少。基于已报道的盐胁迫处理下的RNA-seq数据、549份大豆重测序数据及耐盐指数数据,以及Soybean Expression Atlas数据库中大豆组织表达数据,在大豆中鉴定能够响应盐胁迫的ERF基因。同时,在549份大豆重测序数据中鉴定响应盐胁迫ERF基因的优异等位变异,并分析其驯化与人工选择规律。其中,鉴定出ERF15
盐渍化是土地退化的主要形式之一,目前已成为全球性环境问题,我国盐渍化土地面积占全国内总耕地面积的十分之一,其中华北、西北以及沿海地区较为严
在高盐的胁迫下,植物会发生多种生理变化,如细胞失水、离子毒性、离子失衡、氧化应激、气孔阻塞等,最终引起植物局部组织损伤甚至产生早衰现
乙烯响应因子(ERF, ethyleneresponsive factor)是植物特有的转录因子家族AP2/ERF (APETALA2/ethylene responsive factor)中的一个亚家族,其内包含一个AP2结构
Dong
从Soybean Expression Atlas数据库(https://venanciogroup.uenf.br/cgi-bin/gmax_atlas/index.cgi)中下载11个上调表达的ERF基因(Glyma.01G224100(ERF237)、Glyma.02G066200(ERF218)、Glyma.03G159800(ERF154)、Glyma.06G290000(ERF158)、Glyma.10G016500(ERF196)、Glyma.11G019000(ERF245)、Glyma.12G117000(ERF166)、Glyma.13G236500(ERF170)、Glyma.13G236600(ERF171)、Glyma.13G274100(ERF172)、Glyma.16G147500(ERF232))在大豆发育的7个关键部位(子叶、花、叶、荚、根、种子和芽)的表达数据,使用TBtools 软件绘制ERF家族基因在大豆不同组织器官中的表达量热图。
从 Phytozome网站数据(https://phytozome-next.jgi.doe.gov/)中下载大豆全基因组基本信息gff3文件,利用TBtools软件中的Gene Location Visualize From GTF/GFF插件,对8个在根中富集的ERF基因进行染色体可视化定位。
利用大豆中8个ERF基因的氨基酸序列,在Phytozome数据库中blast多个豆科植物的同源氨基酸序列,包括花生(arahy.Tifrunner.gnm1.ann1.QGFJ76)、鹰嘴豆(Ca_01387,Ca_20544,Ca_01387)、兵豆(Lcu.2RBY.3g058850,Lcu.2RBY.2g004730, Lcu.2RBY.4g046650,Lcu.2RBY.3g058850)、小扁豆(Ler.1DRT.3g064870,Ler.1DRT.2g005700,Ler. 1DRT.3g064870)、蒺藜苜蓿(Medtr2g014340,Medtr 4g078710)、棉豆(Pl06G0000303600.v1,Pl07G0000251100.v1,Pl11G0000268300.v1,Pl01G0000301900.v1)、菜豆(Phvul.007G217800,Phvul.011G125200,Phvul.001G157600)以及豇豆(Vigun07g122000,Vigun06g198800,Vigun11g100500,Vigun06g199000),并导入到MEGA-X软件,生成系统发育树。
前期已对实验室搜集的549份大豆材料(119 份野生大豆、204 份地方品种和 226 份栽培品种)的耐盐指数进行分析,获得两次独立重复的实验数
在549份大豆材料中鉴定ERF158、ERF166和ERF170的不同等位变异,并统计等位变异分别在野生大豆、地方品种和栽培品种的数量,并对不同等位变异分别在野生大豆、地方品种和栽培品种所占比例进行分析。从549份大豆材料全基因组数据中获得候选基因ERF158(Chr06:47878147.. 47879505)、ERF166(Chr12:11748442..11749494)、ERF170(Chr13:34675712..34677522)附近200 kb 的DNA序列,使用vcftools软件,每隔 2000 bp计算一次π值,选择目的基因附近的π值用于核苷酸多样性分析,π值越大,代表核苷酸多样性越丰富。
为了鉴定受盐胁迫诱导表达的关键基因,对栽培大豆品种W82的RNA-seq数

图 1 响应盐胁迫ERF基因的鉴定
Fig. 1 Identification of ERF genes response to salt stress in soybean
A:W82-Water、W82-NaCl和ERF家族基因之间的重叠基因(数字代表不同处理之间产生差异基因的数目);B:40个差异表达的ERF基因热图;C:11个上调基因的组织特异性表达,Glyma.01G224100:ERF237,Glyma.02G066200:ERF218,Glyma.03G159800:ERF154,Glyma.06G290000:ERF158,Glyma.10G016500:ERF196,Glyma.11G019000:ERF245,Glyma.12G117000:ERF166,Glyma.13G236500:ERF170,Glyma.13G236600:ERF171,Glyma.13G274100:ERF172,Glyma.16G147500:ERF232
A: Overlap between W82-Water, W82-NaCl and ERF genes (Numbers represent the number of differently genes that soybean seedings under NaCl and water); B: Heat map of differentially expressed 40 ERF genes under NaCl treatment; C: Tissue specific expression of 11 up-regulated genes,Glyma.01G224100:ERF237,Glyma.02G066200:ERF218,Glyma.03G159800:ERF154,Glyma.06G290000:ERF158,Glyma.10G016500:ERF196,Glyma.11G019000:ERF245,Glyma.12G117000:ERF166,Glyma.13G236500:ERF170,Glyma.13G236600:ERF171,Glyma.13G274100:ERF172,Glyma.16G147500:ERF232
由于已报道的ERF基因大多为响应逆境胁迫的正调控因
对8个在根中表达量较高的ERF基因,利用基因的位置信息,进行基因的染色体位置可视化,结果显示8个ERF基因分布在6条染色体上,其中ERF218在2号染色体上,ERF158在6号染色体上,ERF196在10号染色体上,ERF166在12号染色体上,ERF170、ERF171、ERF172 3个基因在13号染色体上,ERF232在16号染色体上(

图 2 ERF基因的染色体定位和系统发育树
Fig. 2 Chromosome location and phylogenetic tree of ERF genes
A:8个上调基因在染色体上的定位;B:大豆与其他常见豆科植物ERF蛋白的系统发育树,蓝色代表大豆ERF基因
A: Location of the 8 up-regulated genes on chromosomes; B: Phylogenetic tree of ERF protein of soybean and other common leguminous plants, the blue font represents ERF genes in soybean
构建了大豆ERF蛋白与其他豆科ERF蛋白的系统发育树,结果显示,大豆ERF158 和ERF166与棉豆(Phaseolus lunatus L.)、菜豆(Phaseolus vulgaris L.)和豇豆(Vigna unguiculata (Linn.) Walp.)的ERF170具有较高的同源性;大豆ERF170和ERF172与棉豆、兵豆、小扁豆ERF158蛋白以及豇豆、蒺藜苜蓿(Medicago truncatula)、花生(Arachis hypogaea L.)ERF166蛋白具有较高的同源性;大豆ERF171与豇豆、棉豆、菜豆的ERF171蛋白以及鹰嘴豆ERF172具有较高的同源性;大豆ERF196与豇豆、菜豆、棉豆ERF196蛋白聚为一支,同源性较高;大豆ERF218、ERF232与鹰嘴豆(Cicer arietinum)、兵豆(Lens culinaris Medic.)、小扁豆(Lens ervoides)的ERF218、ERF232同源性较高(
为了鉴定8个ERF基因的优异等位变异,在549份大豆重测序材料中,对8个ERF基因进行单倍型分析,以参考基因组W82型作为H1型,结果显示ERF218、ERF196、ERF171、ERF172和ERF232不存在等位变异;ERF158有2个等位变异,ERF15

图 3 ERF基因等位自然变异与耐盐指数的关联分析
Fig. 3 Association of natural variation in ERF gene alleles with salt tolerance
不同小写字母表示在0.05水平显著差异
Different lowercase above the histogram indicate significant difference at 0.05 level
通过单倍型与耐盐指数的关联结果发现,含有ERF15
为了确定ERF158、 ERF166和ERF170 3个基因的优异等位变异在驯化的过程中是否经历了人工选择,对ERF158、ERF166和ERF170优异等位变异在大豆品种中的分布进行分析,结果显示,在94.6%的野生大豆中均为ERF15

图 4 ERF158、ERF166和ERF170不同等位变异的品种分布和核苷酸多态性分析
Fig. 4 Variety distribution and domestication trajectories of different allelic variants of ERF158, ERF166 and ERF170
A:ERF158,ERF166和ERF170不同等位变异分别在野生大豆、地方品种和栽培品种中的分布比例。n代表存在ERF15
A: Proportions of ERF158, ERF166 and ERF170 alleles in wild soybean, local varieties and cultivated varieties, respectively. n represents the number of soybean varieties with ERF15
为了进一步确认ERF158、ERF166和ERF170的优异等位变异在大豆驯化的过程中是否丢失或受到人工选择,对549份大豆重测序数据中,ERF158、ERF166和ERF170附近基因组区域的核苷酸多样性进行分析。结果显示,携带ERF15
ERF转录因子是参与植物应对多种生物和非生物胁迫的重要转录因
在早期的作物驯化过程中,为了获得适合种植在田间的作物,农民针对个别性状进行了强烈的选择,导致栽培大豆中保留下来的基因越来越少,使现代育种中产生了遗传瓶颈效
伴随着驯化基因的强烈选择,大量控制作物优良性状的基因却发生了丢失现象,重新找回丢失的基因可以拓宽和丰富作物的遗传多样性。玉米中控制叶夹角的UPA2基因起源于玉米祖先大刍草,但在驯化的过程中发生了丢失现象,通过连续回交技术,将大刍草UPA2等位基因导入到现代玉米中,使玉米的叶夹角变小,能够高密度种植,提高了玉米的总产
由于自然气候以及人类活动等因素的影响,土地盐碱化日益严重,目前已成为全球面临的重大环境问
在盐诱导下的转录组数据中鉴定到11个上调ERF基因,其中8个ERF基因在根中表达量较高。ERF158、ERF166和ERF170存在优异等位变异,分别为ERF15
参考文献
范王涛. 土壤盐碱化危害及改良方法研究. 农业与技术, 2020, 40(23): 114-116 [百度学术]
Fan W T. Research on the harm of soil salinization and its improvement methods. Agriculture and Technology, 2020, 40(23): 114-116 [百度学术]
Dong L D, Cheng Q, Fang C, Kong L P, Yang H, Hou Z H, Li Y L, Nan H Y, Zhang Y H, Chen Q S, Zhang C B, Kou K, Su T, Wang L S, Li S C, Li H Y, Lin X Y, Tang Y, Zhao X H, Lu S J, Liu B H, Kong F J. Parallel selection of distinct Tof5 alleles drove the adaptation of cultivated and wild soybean to high latitude. Molecular Plant, 2021, 15(2): 308-321 [百度学术]
Cheng Q, Gan Z R, Wang Y P, Lu S J, Hou Z H, Li H Y, Xiang H T, Liu B H, Kong F J, Dong L D. The soybean gene J contributes to salt stress tolerance by up-regulating salt-responsive genes. Frontiers in Plant Science, 2020, 11: 272 [百度学术]
Bahmani K, Noori S A S, Darband A I, Akbari A. Molecular mechanisms of plant salinity tolerance: A review. Australian Journal of Crop Science, 2015, 9(4): 321-336 [百度学术]
Parker M B, Gascho G J, Gaines T P. Chloride toxicity of soybeans grown on Atlantic Coast Flatwoods Soils 1. Agronomy Journal, 1983, 75(3): 439-443 [百度学术]
Wang F, Chen H W, Li Q T, Wei W, Li W, Zhang W K, Ma B, Bi Y D, Lai Y C, Liu X L, Man W Q, Zhang J S, Chen S Y. GmWRKY27 interacts with GmMYB174 to reduce expression of GmNAC29 for stress tolerance in soybean plants. The Plant Journal, 2015, 83(2): 224-236 [百度学术]
Masaru O T, Hideaki S. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. The Plant Cell, 1995, 7(2): 173-182 [百度学术]
Pan Y, Seymour G B, Lu C, Hu Z L, Chen X Q, Chen G P. An ethylene response factor (ERF5) promoting adaptation to drought and salt tolerance in tomato. Plant Cell Reports, 2012, 31(2): 349-360 [百度学术]
Zhai Y, Wang Y, Li Y J, Lei T T, Yan F, Su L T, Li X W, Zhao Y, Sun X, Li J W, Wang Q Y. Isolation and molecular characterization of GmERF7, a soybean ethylene-response factor that increases salt stress tolerance in tobacco. Gene, 2013, 513(1): 174-183 [百度学术]
Wang H T, Ni D Q, Shen J C, Deng S S, Xuan H D, Wang C C, Xu J Y, Zhou L, Guo N, Zhao J M, Xing H. Genome-wide identification of the AP2/ERF gene family and functional analysis of GmAP2/ERF144 for drought tolerance in soybean. Frontiers in Plant Science, 2022, 13: 848766 [百度学术]
Dong L D, Hou Z H, Li H Y, Li Z B, Fang C, Kong L P, Li Y L, Du H, Li T, Wang L S, He M L, Zhao X H, Cheng Q, Kong F J, Liu B H. Agronomical selection on loss-of-function of GIGANTEA simultaneously facilitates soybean salt tolerance and early maturity. Journal of Integrative Plant Biology, 2022, 64(10): 1866-1882 [百度学术]
Mccormack M L, Dickie I A, Eissenstat D M, Fahey T J, Fernandez C W, Guo D L, Helmisaari H S, Hobbie E A, Iversen C M, Jackson R B, Jaana L K, Norby R J, Phillips R P, Pregitzer K S, Pritchard S G, Rewald B, Zadworny M. Redefining fine roots improves understanding of below-ground contributions to terrestrial biosphere processes. The New Phytologist, 2015, 207(3): 505-518 [百度学术]
Wu Y, Li X, Zhang J N, Zhao H Q, Tan S L, Xu W H, Pan J Q, Yang F, Pi E X. ERF subfamily transcription factors and their function in plant responses to abiotic stresses. Frontiers in Plant Science, 2022, 13: 1042084 [百度学术]
Shen L, Zhao E P, Liu R E, Yang X. Transcriptome analysis of eggplant under salt stress: AP2/ERF transcription factor SmERF1 acts as a positive regulator of salt stress. Plants, 2022, 11(17): 2205 [百度学术]
John F D, Brandon S G, Bruce D S. The molecular genetics of crop domestication. Cell, 2006, 127(7): 1309-1321 [百度学术]
Wang M, Li W Z, Fang C, Xu F, Liu Y C, Wang Z, Yang R, Zhang M, Liu S L, Lu S J, Lin T, Tang J Y, Wang Y Q, Wang H R, Lin H, Zhu B G, Chen M S, Kong F J, Liu B H, Zeng D L, Jackson S A, Chu C C, Tian Z X. Parallel selection on a dormancy gene during domestication of crops from multiple families. Nature Genetics, 2018, 50(10): 1435-1441 [百度学术]
Tian J G, Wang C L, Xia J L, Wu L S, Xu G H, Wu W H, Li D, Qin W C, Han X, Chen Q Y, Jin W W, Tian F. Teosinte ligule allele narrows plant architecture and enhances high-density maize yields. Science, 2019, 365(6454): 658-664 [百度学术]
李保国.新时代下盐碱地改良与利用的科学之路.中国农业综合开发, 2022(1): 8-9 [百度学术]
Li B G. The scientific road to the improvement and utilization of saline alkali land in the new era. Agricultural Comprehensive Development in China, 2022(1): 8-9 [百度学术]
高继平,林鸿宣.水稻耐盐机理研究的重要进展——耐盐数量性状基因SKC1的研究.生命科学, 2005(6): 563-565 [百度学术]
Gao J P, Lin H X. Important progress in the study of salt tolerance mechanism in rice: Study on the SKC1 gene for quantitative salt tolerance traits. Chinese Bulletin of Life Sciences, 2005(6): 563-565 [百度学术]
Liaqat A, Alfatih A, Jan S U, Sun L Q, Zhao P X, Xiang C B. Transcription elongation factor AtSPT4-2 positively modulates salt tolerance in Arabidopsis thaliana. BMC Plant Biology, 2023, 23(1): 49 [百度学术]
陈筱冉,张铭笑,严建萍,刘燕蓉,张万军.紫花苜蓿MsCCD4基因克隆及耐盐功能鉴定.草地学报, 2022, 30(10): 2572-2580 [百度学术]
Chen X R, Zhang M X, Yan J P, Liu Y R, Zhang W J. Cloning and identification of salt tolerance function of Alfalfa MsCCD4. Acta Agrestia Sinica, 2022, 30(10): 2572-2580 [百度学术]
Hou Z H, Li Y L, Cheng Y H, Li W W, Li T, Du H, Kong F J, Dong L D, Zheng D F, Feng N J, Liu B H, Cheng Q. Genome-wide analysis of DREB genes identifies a novel salt tolerance gene in wild soybean (Glycine soja). Frontiers in Plant Science, 2022, 13: 821647 [百度学术]
Guan R X, Yu L L, Li X X, Chang R Z, Gilliham M, Qiu L J. Selection of the salt tolerance gene GmSALT3 during six decades of soybean breeding in China. Frontiers in Plant Science, 2021, 12: 794241 [百度学术]