2024年6月13日发(作者:)
doi:10.15889/j.issn.1002-1302.2023.11.007
郑 玲,闫晓曼.甜瓜GRAS家族全基因组的鉴定与表达分析[J].江苏农业科学,2023,51(11):53-59.
甜瓜GRAS家族全基因组的鉴定与表达分析
郑 玲,闫晓曼
(洛阳师范学院生命科学学院,河南洛阳471934)
摘要:为了研究甜瓜(CucumismeloL.)GRAS基因家族的结构特征、进化模式和生长发育过程中的表达模式,对
甜瓜GRAS全基因组的基本信息、染色体位置、共线性、系统发育树、基因结构、表达模式、顺式元件进行了全面分析。
结果表明,从甜瓜中鉴定出37个GRAS基因;这些基因在甜瓜的11条染色体上呈不均匀分布,并在6、7、8号染色体
上出现串联重复;甜瓜GRAS家族可分为HAM、LS、SCR、SCL3、DELLA、LISCL、SHR、PAT18个亚族;基因结构分析显
示,含有1~3个内含子的基因有8个,没有内含子的基因有29个;启动子分析,鉴定到生长素、脱落酸、赤霉素、低温、
干旱、光信号响应作用元件以及参与分生组织和胚乳表达的元件;表达模式分析结果表明,HAM、DELLA亚族的基因
几乎在所有组织中都有表达,LS亚族在根和茎中表达含量较高,PAT1亚族在甜瓜的整个生长发育过程中发挥着重要
作用。得出结论,甜瓜全基因组共获得GRAS基因37个,共可分为8个亚族,不同的亚族具有不同的基因结构、蛋白
质结构和表达模式;不同亚族在甜瓜生长发育的不同阶段具有差异性表达,它们共同参与调控甜瓜的发育过程。
关键词:甜瓜;GRAS家族;生物信息学;表达模式
中图分类号:S652.01 文献标志码:A 文章编号:1002-1302(2023)11-0053-07
1]
GRAS转录因子可以调控基因的转录
[
,其名
[2]
AI、RGA、SCR。GRAS蛋白由400~770称来自G
3-4]
个氨基酸残基组成
[
,C端区域可以分为5个基
R-E、W-G、W-W这3个保守氨基酸组成SAW基
6]
。GRAS转录因子的氨基末端序列差异较大,决序
[
7]
定了GRAS蛋白在植物体中可以发挥不同的作用
[
。
序:LHRI、LHRII、VHIID、PFYRE、SAW,其中VHIID
是GRAS蛋白的核心结构,含有非常保守的P-N-
5]
H-D-Q-L单位,以L-R-I-T-G结尾
[
。富
近年来的研究表明,拟南芥有34个GRAS基
8]9]10]
因
[
,水稻有57个
[
、番茄有53个
[
、杨树有106
11]12]13]
个
[
、大白菜有48个
[
、葡萄有43个
[
、无油樟
14]15]
有34个
[
、苹果有127个
[
。GRAS转录因子在
亮氨酸基序LHRI、LHRII分列于VHIID的两侧,
PFYRE基序由3对蛋白序列特征P、FY、RE组成,
信号转导、根径向伸长、腋生枝分生组织、植物的应
13,16]
激反应过程中发挥着关键作用
[
。
收稿日期:2022-08-22
基金项目:洛阳师范学院2020年度校级骨干教师培养计划;河南省
2021年度省级精品在线开放课程“植物学”;2021年度河南省省派
科技特派员计划。
作者简介:郑 玲(1981—),女,河南南阳人,硕士,副教授,主要从事
植物学研究。E-mail:linda2266@qq.com。
州:华南农业大学,2016:47-48.
[11]潘 健,温海帆,何欢乐,等.黄瓜ERF基因家族鉴定及其在雌
J].中国农业科学,2020,53(1):133-花芽分化中的表达分析[
147.
[12]SunHY,DaiHX,WangXY,etal.Physiologicalandproteomic
analysisofselenium-mediatedtolerancetoCdstressincucumber
(CucumissativusL.)[J].EcotoxicologyandEnvironmentalSafety,
2016,133:114-126.
[13]SunHY,WangXY,WangYN,etal.Alleviationofcadmium
toxicityincucumber(Cucumissativus)seedlingsbytheapplication
CucumismeloL.)又称香瓜,在我国广泛甜瓜(
种植,具有药用价值和食用价值。甜瓜基因组数据
在2012年已公布,但对甜瓜GRAS(CmGRAS)家族
17]
。本研究对甜瓜全基因组的研究还少之又少
[
GRAS转录因子家族在基本信息、染色体定位、共线
櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄
ofselenium[J].SpanishJournalofAgriculturalResearch,2016,14
(4):e1105.
[14]赵普莹.ERFs调控黄瓜果实贮藏冷害及其与膜脂代谢关系研
D].广州:华南农业大学,2018:1-3.究[
[15]蒙林平.油菜素内酯诱导黄瓜幼苗疫病抗性的转录组学研究
[D].广州:华南农业大学,2018:1-2.
[16]李巧丽.外源油菜素内酯对黄瓜耐盐性的调控[D].兰州:西
北师范大学,2020:1-2.
[17]王凯旋.ERFD].扬州:扬
Ⅶ
调控黄瓜耐涝性的分子机制解析[
州大学,2021:64-66.
Copyright©博看网. All Rights Reserved.
性、基因结构、表达模式、顺式元件和系统进化树方
面进行分析,旨在为研究甜瓜GRAS家族基因的功
能提供重要参考。
1 材料与方法
1.1 甜瓜GRAS家族的鉴定
本试验的研究对象是甜瓜GRAS家族。于
2021年3月到6月在洛阳师范学院生命科学学院
生物信息学实验室(214)电脑上进行软件操作和信
息处理。在PlantTFDB数据库(http://planttfdb.gao-
lab.org/)和葫芦科数据库(http://cucurbitgenomics.
org/)分别获得甜瓜GRAS家族的信息,将两者比对
nterPro并删除重复序列,确定候选基因。利用I
(http://www.ebi.ac.uk/interpro/)在线验证候选基
因,最终确定甜瓜GRAS家族成员。使用ExPASy
(http//www.expasy.org/)网站获取甜瓜GRAS蛋白
的等电点(PI)、长度、分子量等理化性质信息。利
用Plant-PLo进行亚细胞定位。
1.2 甜瓜GRAS家族基因在染色体上的分布
在MelonetDB网站(https://melonet-db.dna.
affrc.g-o.jp/ap/top)获取甜瓜家族基因的染色体
信息,制成模板,用MapInspect软件制图。
1.3 甜瓜GRAS家族系统树的构建
利用MEGA6.0软件里提供的ClustalW工具
RAS氨基酸序列进行多序列比对,对拟南芥、甜瓜G
然后采用邻接法构建系统进化树。
1.4 甜瓜GRAS转录因子家族的共线性分析
在甜瓜数据库中下载甜瓜GRAS家族的蛋白质
序列,在TBtools上做比对,再和甜瓜的GRAS基因
文件进行比对,通过TBtools软件进行共线性分析。
1.5 甜瓜GRAS家族的基因结构和保守基序分析
从甜瓜数据库获得甜瓜GRAS家族的CDS序
列和GENE序列,通过GSDS2.0网站(http://gsds.
gao-lab.org/index.php)制图。利用MEME网站
(http://meme-suite.org/tools/meme)对CmGRAS
的保守基序进行分析。
1.6 甜瓜GRAS家族的启动子分析
从甜瓜数据库(CuGenDB)中下载GRAS家族
lantCARE网站(http://基因的启动子序列,通过P
bioinformatics.psb.ugent.be/webtool-s/plantcare/
html/)对其进行预测,得到制作甜瓜GRAS家族启
动子的相关信息,手动操作去除重复值。在GSDS
20网站进行检索,制图。
1.7 甜瓜GRAS家族基因表达模式分析
在MelonetDB网站根据甜瓜的基因序列号检
索出每个基因在甜瓜不同器官发育过程中的表达
数据,选取甜瓜基因在根、茎、叶、雄花、雌花、果实、
利用TBtools软件绘种子这7种器官中的表达数据,
制基因表达热图。
2 结果与分析
2.1 甜瓜GRAS家族成员的鉴定
利用InterPro鉴定出37个GRAS候选基因,在
ExPASy网站分析得出等电点、长度、蛋白质的分子
量、亲水性平均值、脂溶指数、不稳定系数等数据,
再利用Plant-PLo网站对甜瓜GRAS转录因子进
行亚细胞定位。CmGRAS蛋白长度为104
(CmGRAS14)~846(CmGRAS33)aa;蛋白质分子量
为11950.08(CmGRAS14)~92420.94
(CmGRAS33)ku;等电点范围为4.78(CmGRAS7)~
8.73(CmGRAS14),其中2个大于7,呈碱性,35个
,呈酸性;37个CmGRAS的氨基酸序列均为小于7
不稳定蛋白,不稳定系数大于40;仅CmGRAS14的
平均亲水系数为正值,为疏水蛋白,其余均为亲水
7个GRAS基因全蛋白;亚细胞定位显示,甜瓜的3
部在线粒体上。
2.2 甜瓜GRAS家族基因的染色体定位与共线性
分析
CmGRAS基因在11条染色体上分布不均,其中
8号染色体上的基因数量最多,有6个;而10号、12
号染色体上的基因数量最少,只有1个(图1)。在
同一条染色体上的物理位置小于或等于100kb的
18]
同源基因就被认为是串联重复的基因
[
。通过分
析染色体定位图和共线性分析图(图1、图2),发现
在6、7、8号染色体上存在串联重复的基因,总共发
8对共线基因(CmGRAS25/CmGRAS31;现了2
CmGRAS31/CmGRAS17;CmGRAS10/CmGRAS30等),5
个串联重复基因(CmGRAS4、CmGRAS5、CmGRAS20、
CmGRAS21、CmGRAS22)。
2.3 甜瓜GRAS家族系统进化树的分析
8]
参照拟南芥GRAS家族的亚族分类标准
[
,可
将甜瓜的GRAS基因分为8个亚族(图3)。甜瓜
GRAS基因在亚家族中分布不均匀。HAM亚族是
最大的亚族,包括7个CmGRAS成员;其次是SHR、
PAT1亚族,有6个成员;最少的是SCL3亚族,只有
2个成员。甜瓜GRAS家族存在2对同源旁系基因,
Copyright©博看网. All Rights Reserved.
AtGRAS6/CmGRAS16等;SHR亚族中的2对,
AtGRAS27/CmGRAS5、AtGRAS17/CmGRAS32;LISCL
亚族中的1对,AtGRAS14/CmGRAS9;DELLA亚族中
的1对,AtGRAS9/CmGRAS30;SCL3亚族中的1对,
AtGRAS5/CmGRAS17;SCR亚族中的1对,
AtGRAS21/CmGRAS33;LS亚族中的1对,AtGRAS7/
CmGRAS14;HAM亚族中的3对,例如AtGRAS8/
CmGRAS15。
2.4 甜瓜GRAS家族基因结构和保守基序分析
CmGRAS基因结构分析结果(图4)表明,甜瓜
GRAS基因的最大内含子数为3个,29个GRAS基
5个基因(CmGRAS6、CmGRAS8、因没有内含子,
CmGRAS19、CmGRAS24、CmGRAS29)只包含1个内
含子,CmGRAS14包含2个内含子,CmGRAS27、
CmGRAS33包含3个内含子。保守基序分析结果
分别是CmGRAS6/CmGRAS22、CmGRAS21/)表明,37个CmGRAS都有Motif5、Motif6,体现(图4
了GRAS转录因子的高度保守。CmGRAS14中所含
保守基序最少,只有2个(Motif5、Motif6)。除了
CmGRAS20。还确定了15对姐妹基因,分别是PAT1
亚族中的5对,例如AtGRAS30/CmGRAS36、
Copyright©博看网. All Rights Reserved.
CmGRAS14,其他蛋白序列都有Motif1、Motif3、
Motif4、Motif5、Motif6、Motif7。
2.5 甜瓜GRAS家族基因表达模式分析
甜瓜GRAS家族基因在不同组织中都有表达
(图5)。CmGRAS1、CmGRAS28、CmGRAS26、
CmGRAS31、CmGRAS36在根中表达量较高,
CmGRAS1、CmGRAS30、CmGRAS21、CmGRAS13在茎
中表达量较高,CmGRAS1、CmGRAS30在叶中表达量
Copyright©博看网. All Rights Reserved.
较高,CmGRAS2在花中表达量较高,CmGRAS25、
CmGRAS23、CmGRAS3在果实中表达量较高,
CmGRAS3在种子中表达量较高。CmGRAS25、
CmGRAS30、CmGRAS1、CmGRAS34都位于DELLA亚
族,在根、茎、叶、花、果实、种子中都有明显表达,推
ELLA亚族调控了甜瓜多个器官的发育。PAT1测D
亚族在花、果实、种子中表达量较高,推测它们可能
参与了果实和种子的形成。LISCL亚族中的
CmGRAS26和HAM亚族中的CmGRAS21、
CmGRAS13、CmGRAS28在各个组织中都有表达,推
测HAM亚族和SCL亚族广泛参与了甜瓜的生长发
育过程。SCL3亚族中的CmGRAS17在茎中有少量
表达,推测其可能参与了茎的生长发育过程。
2.6 甜瓜GRAS家族基因顺式元件分析
甜瓜GRAS基因中的顺式元件可分为4组(图
6):(1)与光反应有关的,包括ACE(参与光响应的
I-box(光响应元件的一部分);(2)与顺式元件)、
植物抗逆性有关的,包括TC-rich(与植物应激反
应有关的元件)、LTR(与低温胁迫有关的响应元
件)、ARE(应对缺氧环境做出反应的作用元件)、MBS
(应对干旱极端环境做出响应的元件);(3)与植物激
BRE(脱落酸传递和表达相关素反应有关的,包括A
的作用元件)、TATC-box(赤霉素信号传导元件)、
AUXRR-core(生长素传递和表达的响应元件);(4)
ircadion(与昼夜规律与时空基因表达有关的,包括c
有关的作用元件)、CAT-box(参与分生组织表达的
元件)、GCN4-motif(与胚乳表达作用相关的元件)。
CmGRAS17、CmGRAS14、CmGRAS7、CmGRAS21、
CmGRAS37、CmGRAS35、CmGRAS4、CmGRAS29、
CmGRAS2、CmGRAS31、CmGRAS36这11个CmGRAS
基因中同时出现了ABRE,ABRE分布在所有的亚
AUXRR-core存在于CmGRAS36、CmGRAS15、族中;
CmGRAS33、CmGRAS8中。11个CmGRAS基因都出
现了MBS,该元件除了不在SCR亚族中分布,其余
各亚族均有分布;TATC-box只在CmGRAS4、
CmGRAS5、CmGRAS31中有分布;ARE分布在每个亚
S亚族的CmGRAS6中分布最多;LTR分族中,在L
布在DELLA、HAM、SCR、LISCL、SHR这5个亚族
中,其中在CmGRAS30中分布最多;CAT-box只在
CmGRAS6、CmGRAS18、CmGRAS13、CmGRAS24、
Copyright©博看网. All Rights Reserved.
CmGRAS36中分布。
3 讨论
GRAS家族是一种重要的植物特异性转录调节
剂,在植物的生长和发育过程中有非常重要的作
8,19]
用
[
。甜瓜GRAS基因中,78%的基因没有内含
族的顺式作用元件中也发现参与了小麦光反应和
干旱胁迫反应的元件,表明GRAS基因在植物光反
23]
。此外,甘薯应及抗胁迫反应中发挥着关键作用
[
中大多数GRAS基因(85.7%)含有与高温压力有
关的响应元件,大约25%的基因具有LTR,推测它
们可能对低温有较强的反应,27%的基因含有参与
24]
盐胁迫反应的GT1-motif元件
[
。
子,22%的基因含有1~3个内含子,可能是因为古
老基因的水平已经下降或现存的基因中不存在内
20]
含子基因的复制现象
[
,而水稻、扁豆的基因中无
21]
内含子基因的比例却很少,没有超过20%
[
。此
HR亚族和SCR亚族调控根部辐射拟南芥的S
25]
形态的建成
[
,在甜瓜的SHR亚族中CmGRAS5是
AtGRAS27的姐妹基因,SCR亚族中CmGRAS33/
AtGRAS21是一对姐妹基因。基因表达分析结果表
CmGRAS5、CmGRAS33在甜瓜的根部高度表达,明,
推测其可能参与甜瓜根尖分生组织的生长,在维持
根的生长中发挥关键作用。拟南芥的HAM亚族作
26]
用于植株顶部分生组织的发育
[
。甜瓜的HAM亚
外大多数的CmGRAS成员具有相似的外显子-内
mGRAS基因的结构高度保守。甜含子结构,表明C
瓜37个GRAS基因的等电点范围在4.78
(CmGRAS7)~8.73(CmGRAS14)之间,其中2个大
于7,呈碱性,35个小于7,呈酸性,与大麻GRAS基
22]
因等电点的性质
[
基本一致。这表明,不同种类的
族中具有3对姐妹基因CmGRAS15/AtGRAS8、
CmGRAS10/AtGRAS24、CmGRAS28/AtGRAS2。甜瓜
GRAS家族基因的表达分析结果显示,CmGRAS10、
CmGRAS7、CmGRAS21、CmGRAS20这4个基因在茎
中表达量较高,表明甜瓜的HAM亚族可能与拟南
芥的HAM亚族具有类似的功能,也作用于甜瓜顶
GRAS蛋白可能存在相似的理化性质。
RAS顺式作用元件分析结果表明,ABRE甜瓜G
和TC-rich在每个亚族中都有分布,AUXRR-core
只在PAT1、SCR、HAM、DELLA亚族中分布,TATC-
box只在SHR、PAT1亚族中分布。在小麦GRAS家
Copyright©博看网. All Rights Reserved.
部分生组织的生长发育。
[J].BMCPlantBiology,2015,15:209.
拟南芥PAT1亚族在光信号的转导中起调控作
27]
。对甜瓜GRAS基因家族进行启动子分析发用
[
[13]LuJX,WangT,XuZD,etal.Genome-wideanalysisofthe
GRASgenefamilyinPrunusmume[J].MolecularGeneticsand
,2015,290(1):303-317.Genomics
[14]CenciA,RouardM.EvolutionaryanalysesofGRAStranscription
factorsinangiosperms[J].FrontiersinPlantScience,2017,8:273.
[15]ZhangHL,CaoYP,ShangC,etal.Genome-wide
characterizationofGRASfamilygenesinMedicagotruncatula
revealstheirevolutionarydynamicsandfunctionaldiversification
[J].PLoSOne,2017,12(9):e0185439.
[16]ZengX,LingH,ChenXM,etal.Genome-wideidentification,
phylogenyandfunctionanalysisofGRASgenefamilyinDendrobium
Orchidaceae)[J].Gene,2019,705:5-15.catenatum(
[17]Garcia-MasJ,BenjakA,SanseverinoW,etal.Thegenomeof
melon(CucumismeloL.)[J].ProceedingsoftheNational
,2012,109AcademyofSciencesoftheUnitedStatesofAmerica
(29):11872-11877.
[18]李 晨,赵雪惠,王庆杰,等.桃GRAS家族全基因组鉴定与响
应UV-B的表达模式分析[J].中国农业科学,2019,52(24):
4567-4581.
[19]张 斌,陈丽娟,李其华,等.栽培大豆GRAS转录因子家族基
因鉴定及其盐胁迫下表达模式分析[J].江苏农业学报,2021,
37(2):296-309.
[20]SunXL,XueB,JonesWT,etal.Afunctionallyrequired
unfoldomefromtheplantKingdom:intrinsicallydisorderedN-
terminaldomainsofGRASproteinsareinvolvedinmolecular
[J].PlantMolecularBiology,recognitionduringplantdevelopment
2011,77(3):205-223.
[21]ZouM,GuoBC,HeSP.Therolesandevolutionarypatternsof
intronlessgenesindeuterostomes[J].ComparativeandFunctional
Genomics,2011,2011:680673.
[22]JainM,KhuranaP,TyagiAK,etal.Genome-wideanalysisof
].Functional&intronlessgenesinriceandArabidopsis[J
IntegrativeGenomics,2008,8(1):69-78.
[23]李亚飞,阳文龙,顾晶晶,等.小麦GRAS基因家族的全基因组
鉴定与分析[J].麦类作物学报,2019,39(5):549-559.
[24]LiuMY,HuangL,MaZT,etal.Genome-wideidentification,
expressionanalysisandfunctionalstudyoftheGRASgenefamilyin
Tartarybuckwheat(Fagopyrumtataricum)[J].BMCPlant
,2019,19(1):342.Biology
[25]Laura,LaurenzioD.TheSCARECROWgeneregulatesan
asymmetriccelldivisionthatisessentialforgeneratingtheradial
Arabidopsisroot[J].Cell,1996,86(3):423-433.organizationofthe
[26]LiXY,QianQ,FuZM,etal.Controloftilleringinrice[J].
,2003,422(6932):618-621.Nature
[27]BolleC,KonczC,ChuaNH.PAT1,anewmemberoftheGRAS
,isinvolvedinphytochromeAsignaltransduction[J].Genesfamily
&Development,2000,14(10):1269-1278.
[28]PengJ,CarolP,RichardsDE,etal.TheArabidopsisGAIgene
definesasignalingpathwaythatnegativelyregulatesgibberellin
responses[J].Genes&Development,1997,11(23):3194-3205.
现,TATC-box光响应元件也在PAT1亚族有分布,
推测甜瓜的PAT1亚族也作用于光信号的传导。拟
ELLA亚族成员参与了赤霉素信号的传南芥的D
28]
导
[
,甜瓜GRAS家族基因的表达分析结果显示,
甜瓜的DELLA亚族基因各组织中都有表达,推测其
通过赤霉素的信号传导途径调控甜瓜的生长。
参考文献:
[1]SunXL,JonesWT,RikkerinkEHA.GRASproteins:theversatile
rolesofintrinsicallydisorderedproteinsinplantsignalling[J].The
BiochemicalJournal,2012,442(1):1-12.
[2]PyshLD,Wysocka-DillerJW,CamilleriC,etal.TheGRASgene
familyinArabidopsis:sequencecharacterizationandbasicexpression
CARECROW-LIKEgenes[J].ThePlantJournal:analysisoftheS
forCellandMolecularBiology,1999,18(1):111-119.
[3]BolleC.TheroleofGRASproteinsinplantsignaltransductionand
[J].Planta,2004,218(5):683-692.development
[4]ChenY,ZhuPP,WuSY,etal.Identificationandexpression
analysisofGRAStranscriptionfactorsinthewildrelativeofsweet
Ipomoeatrifida[J].BMCGenomics,2019,20(1):911.potato
[5]TianCG,WanP,SunSH,etal.Genome-wideanalysisofthe
rabidopsis[J].PlantMolecularGRASgenefamilyinriceandA
,2004,54(4):519-532.Biology
[6]LeeMH,KimB,SongSK,etal.Large-scaleanalysisoftheGRAS
genefamilyinArabidopsisthaliana[J].PlantMolecularBiology,
2008,67(6):659-670.
[7]GrebT,ClarenzO,SchferE,etal.Molecularanalysisofthe
LATERALSUPPRESSORgeneinArabidopsisrevealsaconserved
controlmechanismforaxillarymeristemformation[J].Genes&
,2003,17(9):1175-1187.Development
[8]LiuXY,WidmerA.Genome-widecomparativeanalysisofthe
opulus,Arabidopsisandrice[J].PlantGRASgenefamilyinP
MolecularBiologyReporter,2014,32(6):1129-1145.
[9]MaHS,LiangD,ShuaiP,etal.Thesalt-anddrought-inducible
poplarGRASproteinSCL7conferssaltanddroughttolerancein
Arabidopsisthaliana[J].JournalofExperimentalBotany,2010,61
(14):4011-4019.
[10]AbarcaD,PizarroA,HernándezI,etal.TheGRASgenefamilyin
:transcriptexpressionpatternsassociatedwiththematuration-pine
relateddeclineofcompetencetoformadventitiousroots[J].BMC
PlantBiology,2014,14:354.
[11]SongXM,LiuTK,DuanWK,etal.Genome-wideanalysisofthe
rassicarapassp.GRASgenefamilyinChinesecabbage(B
pekinensis)[J].Genomics,2014,103(1):135-146.
[12]HuangW,XianZQ,KangX,etal.Genome-wideidentification,
phylogenyandexpressionanalysisofGRASgenefamilyintomato
Copyright©博看网. All Rights Reserved.


发布评论