全 文 :麦类作物学报 2014,34(7):929-935
Journal of Triticeae Crops doi:10.7606/j.issn.1009-1041.2014.07.10
网络出版时间:2014-7-7
网络出版地址:http://www.cnki.net/kcms/doi/10.7606/j.issn.1009-1041.2014.07.02.html
粗山羊草和二穗短柄草SOD基因成员的鉴定与比较分析
收稿日期:2014-03-04 修回日期:2014-04-03
基金项目:浙江省本科院校中青年学科带头人学术攀登项目(pd2013420);台州市科技计划项目(121KY16);台州学院生态学浙江
省重点学科开放课题(EKD2013-03)
第一作者E-mail:jiangming1973@139.com
蒋 明,管 铭,李金枝,金建峰
(台州学院生命科学学院,浙江省植物进化生态学与保护重点实验室,浙江椒江318000)
摘 要:为给超氧化物歧化酶(SOD)基因克隆和功能研究奠定基础,以粗山羊草和二穗短柄草的基因组
序列为材料,利用生物信息学手段对SOD基因家族成员进行鉴定和分析。结果表明,两种植物总共鉴定出
12个SOD基因,其中二穗短柄草7个,粗山羊草5个;它们的编码区全长327~2 145bp,编码108~714个氨
基酸,保守基序数量1~4个,等电点4.84~9.25;Cu/Zn-SOD、Fe-SOD与 Mn-SOD的数量分别为8、3和
1个,在粗山羊草基因组中没有鉴定到 Mn-SOD基因;3种不同类型的SOD在进化树上处于不同分支。除
Bradi2g30580.2和Bradi1g50550.1 外,其他二穗短柄草的SOD基因在粗山羊草中均有同源基因。
关键词:粗山羊草;二穗短柄草;超氧化物歧化酶;基因家族;鉴定
中图分类号:S512.9;S330 文献标识码:A 文章编号:1009-1041(2014)07-0929-07
Identification and Comparative Analysis of SOD Gene Members
fromAegilops tauschi and Brachypodium distachyon
JIANG Ming,GUAN Ming,LI Jinzhi,JIN Jianfeng
(Colege of Life Science,Taizhou University,Zhejiang Provincial Key Laboratory of Plant
Evolutionary Ecology and Conservation,Jiaojiang,Zhejiang 318000,China)
Abstract:Superoxide dismutase(SOD)enzymes act as important antioxidants thus protect cels from
being injured by reactive oxygen species,and they play a vital role in stress defense.To provide a ba-
sis for gene cloning and function analysis,genome sequences of two Gramineae plants,Aegilops taus-
chii and Brachypodium distachyon,were chosen to identify SOD gene family members applying bioin-
formatics methods.Results showed that there were 12SOD genes in total,and 7of them were from
B.distachyon while the other 5fromA.tauschii.The complete coding sequences were 327~2 145bp
in length encoding 108~714amino acids with conserved motif number of 1~4;Their isoelectric point
values ranged from 4.84to 9.25.The amount of Cu/Zn-SOD,Fe-SOD and Mn-SOD were 8,3and 1,
respectively,and no Mn-SOD was indentified fromA.tauschii;Three categories of SOD were located
at three different phylogenetic clades.Homologous SOD genes of B.distachyon could be found in
A.tauschii except for Bradi2g30580.2 and Bradi1g50550.1.
Key words:Aegilops tauschii;Brachypodium distachyon;Superoxide dismutase;Gene family;Iden-
tification
植物在正常生长发育和逆境条件下均会产生
活性氧(Reactive O2species,ROS),ROS包括超
氧阴离子、羟自由基、过氧化氢和过氧化物自由基
等,它们攻击脂类、碳水化合物、核酸和蛋白质,造
成细胞膜脂过氧化、碱基突变、DNA链断裂和蛋
白质损伤等不利后果[1]。超氧化物歧化酶(Su-
peroxide dismutase,SOD,EC 1.15.1.1)广泛存
在于原核和真核生物中,是一类重要的抗氧化剂,
它们通过歧化作用将过氧化物还原成氧和过氧化
氢,从而保护细胞免受活性氧的损伤,在逆境反应
中起着重要作用[2]。根据金属离子辅基的不同,
SOD在高等植物中分为三类,分别是 Cu/Zn-
SOD、Fe-SOD和Mn-SOD,Cu/Zn-SOD存在于细
胞质、叶绿体、过氧化物体和质外体,Fe-SOD分
布在叶绿体和过氧化物体中,而 Mn-SOD定位于
线粒体和过氧化物体[3]。近年来,已从水稻
(Oryza sativa)、大麦(Hordeum vulgare)、小麦
(Triticum aestivum)、獐茅(Aeluropus littora-
lis)和玉米(Zea mays)等禾本科植物中克隆到了
SOD基因[4-8]。
粗山羊草(Aegilops tauschii,2n=2x=14,
DD)为禾本科(Gramineae)小麦族(Triticeae)山
羊草属植物,是普通小麦(T.aestivum,2n=6x
=42,AABBDD)D染色体组的供体[9]。粗山羊
草遗传变异丰富、生长迅速、抗逆性强,是小麦品
种改良的重要基因资源[10]。研究人员在粗山羊
草形态特征、遗传改良、种质资源多样性和基因克
隆等方面做了大量研究[11-14]。最近,粗山羊草的
全基因组测序已经完成,为研究D基因组的多样
性、小麦品种改良和基因功能研究奠定了基
础[15]。二穗短柄草(Brachypodium distachyon,
2n=2x=10)为禾本科短柄草属植物,具有植株
小、生育期短和易于遗传转化等优点,是麦类作物
比较基因组和功能基因组研究的重要模式种,已
于2010年完成基因组测序[16-17]。本研究拟以粗
山羊草和二穗短柄草的基因组序列为材料,利用
生物信息学手段对SOD进行全基因组鉴定和比
较分析,为基因克隆与功能研究奠定基础。
1 材料与方法
1.1 材 料
粗山羊草基因和编码蛋白序列下载自ftp://
climb.genomics.cn,二穗短柄草基因和蛋白质序
列下载自ftp://brachypodium.org/brachypodi-
um.org;SOD的隐马尔可夫模型(Hidden Mark-
ov Model,HMM)文件 PF00080、PF02777 和
PF00081下载自Pfam蛋白家族数据库网站ht-
tp://pfam.sanger.ac.uk。
1.2 方 法
利用 HMMER 3.0的hmmsearch程序[18],
分别以PF00080、PF02777和PF00081为种子文
件在粗山羊草与二穗短柄草的蛋白质数据库中获
得候选序列;将候选序列逐条用SMART在线工
具(http://smart.embl.de)[19]分析,去除没有
SOD结构域的蛋白质序列,同时获得各自的编码
区;分子量和等电点的预测用在线工具Compute
pI/Mw (http://web.expasy.org/compute _
pi/)[20]。
Clustal X 1.81软件[21]用来对齐所有的SOD
序列;用 Mega 3.1软件[22]构建进化树,方法为邻
接法,自举检测次数为 1 000。利用 MEME
SUITE 在 线 工 具 (http://meme.sdsc.edu/
meme/cgi-bin/meme.cgi)预测SOD蛋白的保守
模体[23]。
2 结果与分析
2.1 SOD基因家族成员鉴定结果
在 Windows环境下,利用hmmsearch程序
在本地数据库中搜索 SOD 候选序列,并用
SMART在线工具剔除没有SOD结构域的序列,
共得到12条目标序列,其中二穗短柄草7条,粗
山羊草5条(表1)。NCBI比对结果表明,Bra-
di1g69680.1、Bradi3g43070.3、Bradi1g18340.3、Bra-
di5g18900.3、AEGTA29877、AEGTA00133、AEG-
TA30722和AEGTA28552为Cu/Zn-SOD,它们各有
1个SOD结构域;Bradi1g51140.1、Bradi1g50550.1、
Bradi2g30580.2 和AEGTA12403 为Fe/Mn-SOD,
它们均有2个SOD结构域。
2.2 SOD基因的特点
二穗短柄草SOD基因分布于第1、2、3和5
号染色体,其中1号染色体最多,有4个,其余染
色体各1个;粗山羊草SOD基因分布于5个不同
的scaffold上。基因编码区(Coding DNA Se-
quence,CDS)长度327~2 145bp,编码108~
714个氨基酸,其中,粗山羊草AEGTA12403 的
CDS最长,AEGTA00133 次之,为 2 100bp,
AEGTA28552最短。在线工具Compute pI/Mw
预测分子量和等电点的结果表明,12个SOD基
因编码蛋白的分子量为10 870.36~80 773.35
D,其中 AEGTA12403 的分子量最大,AEG-
TA00133 次 之 (为 75 368.19 D),而 AEG-
TA28552 的最小;等电点范围为4.84~9.25,
·039· 麦 类 作 物 学 报 第34卷
Bradi2g30580.2最大,Bradi1g50550.1 次之(为 9.19),AEGTA30722 最小。
表1 粗山羊草和二穗短柄草SOD基因家族成员的特征
Table 1 Characteristics of SOD gene family members fromAegilops tauschi and Brachypodium distachyon
基因名称
Gene name
位置
Locus
编码区
长度
Coding
sequence
length/bp
编码蛋
白长度
Protein
sequence
length/aa
编码蛋白
分子量
Molecular
weight of
decuced
protein/D
SOD结构
域位置
Location of
SOD
domain
等电点
Isoelectric
point
Bradi1g69680.1 chr01_68164994~68169194 495 164 16 527.39 9~158 6.38
Bradi3g43070.3 chr03_44386396~44388692 621 206 20 906.43 54~199 5.79
Bradi1g18340.3 chr01_14737275~14733706 330 109 11 319.48 4~102 4.91
Bradi5g18900.3 chr05_22011902~22017609 978 325 34 769.44 179~304 5.38
Bradi1g51140.1 chr01_49562344~49564604 771 256 29 230.21 43~128、133~239 8.65
Bradi1g50550.1 chr01_49166110~49163572 1 197 398 43 916.14 122~207、212~333 9.19
Bradi2g30580.2 chr02_30235286~30237900 552 183 19 737.72 29~110、117~181 9.25
AEGTA29877 scaffold72866_4558~8520 522 173 17 966.99 10~157 6.52
AEGTA00133 scaffold121720_7996~17205 2 100 699 75 368.19 550~695 5.67
AEGTA30722 scaffold72512_17307~21165 732 243 25 774.20 97~223 4.84
AEGTA28552 scaffold15106_27387~28276 327 108 10 870.36 29~108 6.81
AEGTA12403 scaffold58517_30346~36018 2 145 714 80 773.35 505~585、591~697 5.91
2.3 SOD基序的特征
序列间的相似区域可能具有特定的功能,利
用MEME SUITE在线工具可从一组序列中发现
保守区域并生成基序(motif),用图示的方式呈
现,字母的高度表示该氨基酸在模体中出现的相
对次数。搜索结果表明,粗山羊草和二穗短柄草
共有3类基序,它们的正则表达式分别为[HW]G
[FK]H[IQ]H[EA][FY][GV]DT[TL]NG
[CA]、TG[PKS][HV][FY]NP[AH][GDN]
[KEY]SH[GR][AK]P和[FY]L[TD][QY]E
[DN]D[GR]P[TA][TY]VS[GA]R[IV][ST]
GL[AV][PS](图1)。除AEGTA28552、AEG-
TA12403 和Bradi2g30580.2 外,其余SOD的基
序数均为3个,其中AEGTA12403 的基序数为
4,AEGTA28552 的基序为2,而Bradi2g30580.2
只有1个基序(图4)。
2.4 系统发育分析
用Clustal X 1.81软件对比二穗短柄草和粗
山羊草的Cu/Zn-SOD序列与Fe/Mn-SOD(图2
和图3),结果表明,8条Cu/Zn-SOD蛋白质序列
的氨基酸组成差异较大,但个别位点十分保守,如
+550位(缬氨酸)、+566位(甘氨酸)、+587位
(甘氨酸)和+588位(亮氨酸)等;4条Fe/Mn-
SOD序列中,Bradi1g51140.1、Bradi1g50550.1
和AEGTA12403 为 Fe-SOD,C端序列相对保
守,同源性高,而Bradi2g30580.2 为 Mn-SOD,与
另3条序列的差异较大,但这4条多肽链的部分
区域十分保守,如+513~+516均为 ALEP残
基,+579~ +580 为 NH,+609~ +612 为
DFGS,+629~+635为GSGWVWI。
用 Mega 3.1软件构建了进化树,建树方法
为邻接法(图4)。结果表明,3个Fe-SOD(Bra-
di1g51140.1、Bradi1g50550.1 和AEGTA12403 )
聚为一组,8个Cu/Zn-SOD(Bradi1g69680.1、Bra-
di3g43070.3、Bradi1g18340.3、Bradi5g18900.3、
AEGTA-29877、AEGTA00133、AEGTA30722和
AEGTA-28552 )聚为一组,而唯一的 Mn-SOD
蛋白Bradi2g30580.2 单独处于一个分支,再与3
个Fe-SOD聚为一组。8个Cu/Zn-SOD中,Bra-
di1g69680.1 与AEGTA29877、Bradi3g43070.3
与AEGTA28552、Bradi1g18340.3 与AEGTA-
00133、Bradi5g18900.3 与AEGTA30722 两两聚
于同一小分支,它们互为同源序列;在 Fe-SOD
中,二穗短柄草的Bradi1g51140.1 与粗山羊草的
AEGTA12403 聚 于 一 小 分 支,再 与 Bra-
di1g50550.1 聚为一组。
·139·第7期 蒋 明等:粗山羊草和二穗短柄草SOD基因成员的鉴定与比较分析
X轴上的数字表示次数,字母的高度表示在模体中出现的相对次数;Y轴上的数字表示氨基酸残基的位置
The number on x axis indicates the number,and the height of a letter indicates its relative frequency at the given position(x-axis)in
the motif;the number on y axis indicates the position of amino acid
图1 粗山羊草和二穗短柄草的保守模体
Fig.1 Conserved motifs of Aegilops tauschi and Brachypodium distachyon
黑色与灰色底纹分别表示氨基酸具100%和大于60%相似性;数字表示氨基酸残基所处位置;连字符代表缺失;序列下面的大写
字母表示序列中氨基酸残基完全一致,小写字母表示至少6个氨基酸一致。图2同此
Black and gray shading indicate the similarity among amino acid residues of 100%and more than 60%,respectively;numbers refer
to amino acid positions;hyphens represent gaps;uppercase letters refer to amino acids identical in al amino acid residues,and lowercase
ones indicate identical in at least 6residues.The same are as in Fig.3
图2 粗山羊草和二穗短柄草Cu/Zn-SOD的序列比对
Fig.2 Sequence alignment of Cu/Zn-SOD fromAegilops tauschi and Brachypodium distachyon
·239· 麦 类 作 物 学 报 第34卷
图3 粗山羊草和二穗短柄草Fe/Mn-SOD的序列比对
Fig.3 Sequence alignment of Fe/Mn-SOD fromAegilops tauschi and Brachypodium distachyon
图4 粗山羊草和二穗短柄草SOD的系统发育分析
Fig.4 Phylogenetic tree of SOD genes fromAegilops
tauschi and Brachypodium distachyon
3 讨 论
分子克隆、序列表达标签(Expressed Se-
quence Tags,EST)和基因组测序等为SOD的全
基因组鉴定奠定了基础。近年来,研究人员对拟
南芥(Arabidopis thaliana)[24]、毛果杨(Populus
trichocarpa)[25]和巴尔干苦苣苔(Haberlea rho-
dopensis)[26]等植物的SOD基因家族进行了鉴定
分析。Kliebenstein等[24]从拟南芥中鉴定出7个
成员,包括3个Cu/Zn-SOD、3个Fe-SOD及1个
Mn-SOD;Molina-Rueda等[25]以毛果杨基因组序
列为材料,共鉴定出12个成员,其中7个为Cu/
Zn-SOD、3个为Fe-SOD,另2个为 Mn-SOD;利
用RT-PCR和 RACE方法,Apostolova等[26]从
复苏植物巴尔干苣苔中分离到7个SOD基因,包
括4个Cu/Zn-SOD、1个Fe-SOD、1个 Mn-SOD
及1个假基因。
在基因家族的进化过程中,经历诸如染色体
倍增、基因复制、碱基插入和缺失等事件,家族成
员数量、基因大小和序列组成等可能产生差
异[27-29]。本研究以两种禾本科植物的基因组序列
·339·第7期 蒋 明等:粗山羊草和二穗短柄草SOD基因成员的鉴定与比较分析
为材料,共得到12个SOD基因,它们的编码区全
长差异较大,为327~2 145bp。从粗山羊草中鉴
定出5个SOD基因,其中4个为Cu/Zn-SOD,1
个为Fe-SOD;从二穗短柄草鉴定出7个SOD基
因,其中4个为Cu/Zn-SOD,2个为Fe-SOD,另
一个为 Mn-SOD。从水稻品种Cpslo17和獐茅中
克隆得到的OsCu/ZnSOD与AlSOD均编码152
个氨基酸,序列长度与同为Cu/Zn-SOD的二穗
短柄草Bradi1g69680.1 接近[4,7]。大麦的 Hv-
SOD和小麦的SOD1.1均编码201个氨基酸,长
度与Bradi3g43070.3 相差5个氨基酸残基,它们
均为Cu/Zn-SOD[5-6]。玉米SOD4和SOD4A序
列十分相似,仅19和29位氨基酸不同,蛋白质序
列的长度均为153,它们是Cu/Zn类型的SOD,
与水稻OsCu/ZnSOD和獐茅AlSOD的长度仅差
一个残基[8],但其与粗山羊草中鉴定得到的AE-
GTA00133 和AEGTA12403 在序列长度上差异
较大,后两者的氨基酸长度分别为699和714。
与毛果杨相比,粗山羊草和二穗短柄草的SOD基
因数量相对较少,但与拟南芥和巴尔干苦苣苔相
仿。在 Cu/Zn-SOD 中,粗 山 羊 草 的 AEG-
TA29877、AEGTA28552、AEGTA00133 和AE-
GTA3072 在二穗短柄草中均有对应的同源序列,
在进化树上两两组合,它们的序列长度、组成和结
构相似;但在Fe-SOD中,Bradi1g50550.1 在粗
山羊草中没有同源序列;同样,粗山羊草中缺少二
穗短柄草Bradi2g30580.2 基因的同源序列,这可
能是进化过程中基因丢失所致[30]。
本研究以粗山羊草和二穗短柄草的基因组序
列为材料,共鉴定出12个SOD基因。SOD的全
基因组鉴定与比较分析,为粗山羊草和二穗短柄
草SOD基因的分子克隆、表达分析和功能鉴定奠
定了基础。
参考文献:
[1]Du X M (杜秀敏),Yin W X(殷文璇),Zhang H (张 慧),et
al.The researching progress of superoxide dismutase [J].
Journal of Chinese Biotechnology(中国生物工程杂志),
2003,23(1):48-50(in Chinese).
[2]Dionisi O,Galeotti T,Terranova T,et al.Superoxide radicals
and hydrogen peroxide formation in mitochondria from normal
and neoplastic tissues[J].Biochimica et Biophysica Acta,
1975,403(2):292-300.
[3]Corpas F J,Fernández-Ocana A,Carreras A,et al.The ex-
pression of different superoxide dismutase forms is cel-type
dependent in olive(Olea europaea L.)leaves[J].Plant and
Cell Physiology,2006,47(7):984-994.
[4]Wang F(王 峰),Wang H B(王宏斌),Wang J F(王金发).
Sequence and expression analysis of cytoplasmic copper/zinc
superoxide dismutase gene in rice[J].Journal of Tropical
and Subtropical Botany(热带亚热带植物学报),2007,15(2):
101-106(in Chinese).
[5]Abu-Romman S M,Shatnawi M.Isolation and expression a-
nalysis of chloroplastic copper/zinc superoxide dismutase gene
in barley[J].South African Journal of Botany,2011,77(2):
328-334.
[6]Wu G,Wilen R W,Robertson A J,et al.Isolation,chromo-
somal localization,and differential expression of mitochondrial
manganese superoxide dismutase and chloroplastic copper/zinc
superoxide dismutase genes in wheat[J].Plant Physiology,
1999,120(2):513-520.
[7]Modarresi M,Nematzadeh G A,Moradian F,et al.Identifica-
tion and cloning of the Cu/Zn superoxide dismutase gene from
halophyte plant Aeluropus littoralis[J].Genetika,2012,48
(1):130-134.
[8]Cannon R E,Scandalios J G.Two cDNAs encode two nearly i-
dentical Cu/Zn superoxide dismutase proteins in maize[J].
Molecular Genetics and Genomics,1989,219(1-2):1-8.
[9]ter Steege M W,den Ouden F M,Lambers H,et al.Genetic
and physiological architecture of early vigor in Aegilops taus-
chii,the D-genome donor of hexaploid wheat.A quantitative
trait loci analysis[J].Plant Physiology,2005,139(2):1078-
1094.
[10]Malik R,Smith C M,Brown-Guedira G L,et al.Assessment
of Aegilops tauschii for resistance to biotypes of wheat curl
mite(Acari:Eriophyidae)[J].Journal of Economic Ento-
mology,2003,96(4):1329-1333.
[11]Naghavi M R,Amirian R.Morphological characterization of
accessions of Aegilops tauschii[J].International Journal of
Agriculture and Biology,2005,7(3):392-394.
[12]Atiq-ur-Rahman R,Ashraf M,Rizwan S,et al.Evaluation of
synthetic hexaploid wheats(Triticum turgidum L.x Ae-
gilops tauschii L.)and their durum parents for stripe rust
(Puccinia striiformis Westend.f.sp.tritici Erikson)Resist-
ance[J].Revista Mexicana de Fitopatología,2007,25(2):
152-160.
[13]Su Y R(苏亚蕊),Zhang D L(张大乐),Xu S M (徐守明),
et al.Genetic diversity and differentiation in different Ae-
gilops tauschii populations revealed by SSR[J].ACTA Eco-
logica Sinica(生态学报),2010,30(4):969-975(in Chinese
with English abstract).
[14]Li W,Liu A J,Sheng Y Z,et al.Molecular Characterization
of Waxy Gene in Aegilops tauschii [J].Asian Journal of
Plant Sciences,2013,12:34-39.
[15]Jia J,Zhao S,Kong X,et al.Aegilops tauschii draft genome
sequence reveals a gene repertoire for wheat adaptation[J].
Nature,2013,496(7443):91-95.
[16]Chen J Y(陈军营),Li X N(李香妞),Zhao Y D(赵一丹),
·439· 麦 类 作 物 学 报 第34卷
et al.A new model plant in Gramineae-Brachypodium dis-
tachyon L.[J].Plant Physiology Communications(植物生
理学通讯),2008,44(4):781-784(in Chinese with English
abstract).
[17]International Brachypodium Initiative,Genome sequencing
and analysis of the model grass Brachypodium distachyon
[J].Nature,2010,463(7282):763-788.
[18]Finn R D,Clements J,Eddy S R.HMMER web server:inter-
active sequence similarity searching[J].Nucleic Acids Re-
search,2011,39(Web Server issue):W29-W37.
[19]Letunic I,Doerks T,Bork P.SMART 7:recent updates to
the protein domain annotation resource[J].Nucleic Acids
Research,2012,40(Database issue):D302-D305.
[20]Artimo P,Jonnalagedda M,Arnold K,et al.ExPASy:SIB
bioinformatics resource portal[J].Nucleic Acids Research,
2012,40(W1):W597-W603.
[21]Thompson J D,Gibson T J,Plewniak F,et al.The CLUST-
AL_X windows interface:flexible strategies for multiple se-
quence alignment aided by quality analysis tools[J].Nucleic
Acids Research,1997,25(24):4876-4882.
[22]Kumar S,Tamura K,Nei M.MEGA3:Integrated software
for molecular evolutionary genetics analysis and sequence a-
lignment[J].Briefings in Bioinformatics,2004,5(2):150-
163.
[23]Alscher R G,Erturk N,Heath L S.Role of superoxide dis-
mutases(SODs)in controling oxidative stress in plants[J].
Journal of Experimental Botany,2002,53(372):1331-
1341.
[24]Kliebenstein D J,Monde R A,Last R L.Superoxide dis-
mutase in Arabidopsis:an eclectic enzyme family with dispa-
rate regulation and protein localization[J].Plant Physiolo-
gy,1998,118(2):637-650.
[25]Molina-Rueda J J,Tsai C J,Kirby E G.The Populus super-
oxide dismutase gene family and its responses to drought
stress in transgenic poplar overexpressing apine cytosolic
glutamine synthetase(GS1a)[J].PLoS One,2013,8(2):
e56421.
[26]Apostolova E,Rashkova M,Anachkov N,et al.Molecular
cloning and characterization of cDNAs of the superoxide dis-
mutase gene family in the resurrection plant Haberlea rho-
dopensis[J].Plant Physiology and Biochemistry,2012,55:
85-92.
[27]Hahn M W,De Bie T,Stajich J E,et al.Estimating the tempo
and mode of gene family evolution from comparative genomic
data[J].Genome Research,2005,15(8):1153-1160.
[28]Guo Y L.Gene family evolution in green plants with empha-
sis on the origination and evolution of Arabidopsis thaliana
genes[J].Plant Journal,2013,73(6):941-951.
[29]Demuth J P,Hahn M W.The life and death of gene families
[J].Bioessays,2009,31(1):29-39.
·539·第7期 蒋 明等:粗山羊草和二穗短柄草SOD基因成员的鉴定与比较分析