• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

静原鸡AK1基因序列分析及真核表达载体的构建

张娟, 虎红红, 邓占钊, 母童, 顾亚玲, 辛国省

张娟, 虎红红, 邓占钊, 等. 静原鸡AK1基因序列分析及真核表达载体的构建[J]. 华南农业大学学报, 2021, 42(2): 17-25. DOI: 10.7671/j.issn.1001-411X.202006037
引用本文: 张娟, 虎红红, 邓占钊, 等. 静原鸡AK1基因序列分析及真核表达载体的构建[J]. 华南农业大学学报, 2021, 42(2): 17-25. DOI: 10.7671/j.issn.1001-411X.202006037
ZHANG Juan, HU Honghong, DENG Zhanzhao, et al. Sequence analysis of AK1 gene from Jingyuan chicken and construction of its eukaryotic expression vector[J]. Journal of South China Agricultural University, 2021, 42(2): 17-25. DOI: 10.7671/j.issn.1001-411X.202006037
Citation: ZHANG Juan, HU Honghong, DENG Zhanzhao, et al. Sequence analysis of AK1 gene from Jingyuan chicken and construction of its eukaryotic expression vector[J]. Journal of South China Agricultural University, 2021, 42(2): 17-25. DOI: 10.7671/j.issn.1001-411X.202006037

静原鸡AK1基因序列分析及真核表达载体的构建

基金项目: 国家自然科学基金(31860621)
详细信息
    作者简介:

    张娟(1982—),女,博士,E-mail: zhangjkathy@126.com

    通讯作者:

    辛国省(1979—),男,副教授,博士,E-mail: gsxin@nxu.edu.cn

  • 中图分类号: S831

Sequence analysis of AK1 gene from Jingyuan chicken and construction of its eukaryotic expression vector

  • 摘要:
    目的 

    探究宁夏地方品种静原鸡前期转录组测序筛选出的AK1基因编码蛋白的结构与功能,构建其真核表达载体。

    方法 

    根据GenBank上已公布的原鸡AK1基因序列,针对其CDS区设计特异性引物,通过克隆测序对AK1基因CDS区SNPs进行快速筛查,构建AK1基因的真核表达载体,并进行编码区生物信息学功能分析。

    结果 

    AK1基因编码区全长585 bp,编码194个氨基酸;AK1蛋白无跨膜结构域,存在2个CpG岛,18个磷酸化位点,10个抗原表位,为稳定的水溶性蛋白,空间结构以α–螺旋和无规则卷曲为主。亚细胞定位结果显示,AK1蛋白主要位于细胞质内;GO富集分析发现,AK1基因同样富集在细胞质中,与亚细胞定位结果一致。基因共表达分析发现,在与AK1互作的基因中,AK1AMPD1PKM2基因存在共表达,共表达系数分别为0.116和0.063。成功构建出AK1-pEGFP-N1载体。

    结论 

    该研究结果为后续AK1蛋白功能以及AK1作为肌苷酸相关基因的深入研究提供了科学依据。

    Abstract:
    Objective 

    To explore the structure and function of protein encoded by AK1 gene screened by transcriptome sequencing of Ningxia local breed Jingyuan chicken, and construct its eukaryotic expression vector.

    Method 

    According to the original chicken AK1 gene sequence published on GenBank, specific primers were designed for its CDS region, and the CDS region SNPs of AK1 gene were screened quickly by cloning and sequencing. The eukaryotic expression vector of AK1 gene was constructed. The bioinformatics function of the coding region was analyzed.

    Result 

    The full length of AK1 gene coding region was 585 bp, which encoded 194 amino acids. AK1 protein had no transmembrane domain, and it had two CpG islands, 18 phosphorylation sites and 10 antigenic epitopes. AK1 protein was a stable water-soluble protein, and the spatial structure was mainly α-helix and irregular crimp. The results of subcellular localization showed that AK1 protein was mainly located in the cytoplasm, and GO enrichment analysis showed that AK1 gene was also enriched in the cytoplasm, which was consistent with the result of subcellular localization. Gene co-expression analysis of genes interacting with AK1 showed that AMPD1 and PKM2 coexpressed with AK1 gene, and the co-expression coefficients were 0.116 and 0.063, respectively. The AK1-pEGFP-N1 vector was successfully constructed.

    Conclusion 

    The results of this study provide a scientific basis for further study of the function of AK1 protein and AK1 as an inosinic acid-related gene.

  • 近年来,在我国经济迅速发展的过程中,矿山开采和工业废物等带来的土壤和环境重金属污染日益严重,对农业生产和人们的身体健康构成严重威胁。其中,重金属Cd是毒性很强的环境污染物,通过空气沉降、灌溉和施肥等途径进入土壤,引起土壤污染。土壤中高含量Cd对植物的毒性主要表现为植物褪绿、坏疽、落叶、生长量减少、根伸长受阻等[1]。Cd污染农业环境后会不同程度地影响水稻生长及籽粒Cd含量,由此对以稻米为食的人们的生活与健康造成严重危害,如20世纪60年代日本稻田Cd污染导致的骨痛病。水稻在各个生长阶段都会受到土壤中Cd的胁迫,正常的生长发育受到严重影响。水稻种子萌发过程中,Cd对种子发芽率影响较小,但会显著抑制胚根及胚芽的生长,对胚根的抑制作用明显大于胚芽[2-4]。水稻幼苗期遭受Cd胁迫后,根系生长受到抑制,苗高和根长增长减慢[5]。水稻籽粒Cd积累在不同品种间存在显著差异,研究发现了控制籽粒Cd积累的主效基因[6-7]。籽粒低Cd积累是水稻耐Cd研究的重点,从水稻育种和生产需要来看,各个生育时期均耐Cd的材料是最理想的材料。因此,前期耐Cd性与籽粒低Cd积累的相关性研究有助于早期筛选籽粒低Cd积累的材料。程旺大等[8]利用苗期耐Cd性和籽粒Cd含量差异明显的4个晚粳稻品种研究了Cd胁迫对水稻生长和营养代谢的影响,发现Cd胁迫下各品种的产量、每株穗数、每穗粒数、结实率、粒质量及地上部干物质质量等显著降低,苗期耐Cd性强的品种受到的抑制作用小于耐Cd性弱的品种,耐性相近的品种中,籽粒Cd富集程度较高的品种也严重受害。该研究预示可以探索用苗期耐Cd性评价不同品种籽粒Cd含量的差异。

    水稻耐Cd性是多基因控制的数量性状,在耐Cd种质资源筛选、QTL定位、耐Cd基因克隆与功能等研究方面已取得一定进展[6-7, 9-21],上述研究有助于解析水稻耐Cd性的遗传控制基因及作用机理,为育种改良水稻耐Cd性提供重要的指导依据。野生稻是栽培稻的祖先,在长期自然进化过程中适应了多种自然环境,积累了较强的抗逆性,是天然的基因库[22]。但是目前野生稻耐Cd性研究的报道还较少。广东省植物分子育种重点实验室前期构建了以展颖野生稻Oryza glumaepatula为供体的染色体单片段代换系(Single-segment substitution lines, SSSLs)文库[23],本研究对部分SSSLs在芽期进行了不同浓度CdCl2溶液的胁迫处理,对发芽第7天的种子根长及芽长进行分析,以相对根长与相对芽长为指标,筛选出多份芽期耐Cd性强的SSSLs,并根据代换片段的信息,进行初步的QTL鉴定分析。本研究有助于扩大水稻耐Cd资源的范围,获得新的优异基因,为水稻耐Cd育种提供参考。

    栽培稻品种‘华粳籼74’、‘新黄占’、‘中花11’、‘台中65’、‘粤晶丝苗2号’和‘象牙香占’用于确定芽期耐Cd试验中Cd溶液处理浓度与培养条件及Cd耐性分析指标。以‘华粳籼74’为受体、展颖野生稻IRGC104387为供体构建的SSSLs材料共37份。

    挑选颖壳光亮、大小均匀、无斑点、饱满无损伤的水稻种子,用蒸馏水洗涤1遍后用70%(φ)乙醇溶液浸泡20 s,倒去乙醇溶液,用蒸馏水洗去残留的乙醇溶液,重复3次,用培养皿法进行发芽试验。栽培稻品种试验中,CdCl2浓度设为0(对照,只加蒸馏水)、30、50、100和500 μmol·L−1。根据栽培稻品种的浓度梯度的试验结果,将野生稻SSSLs CdCl2处理浓度设为0(对照,只加蒸馏水)、50、100 μmol·L−1。每个培养皿放置20~30粒种子,1个培养皿为1次重复,每份材料每个浓度重复3次。浸种时,在各培养皿中加入相同体积蒸馏水或CdCl2溶液以淹没种子,放在(25±1) ℃的光照培养箱中培养,每日观察,待露白后,将露白的种子转移到垫有2层无菌滤纸的培养皿中。每天补充处理溶液,保持滤纸湿润。每日光照/黑暗时间设为16/8 h,光照强度为5 000 lx。

    以胚芽达到种子长度一半作为发芽的标准。于发芽后第7天统计发芽率,并用直尺测量已发芽种子的根长与芽长。测量根长时,测量最长根长。每份材料每个浓度测30棵幼苗。统计分析不同材料在不同CdCl2浓度处理下的根长、芽长,并计算根芽比、相对根长及相对芽长,以反映其耐Cd特性。

    $$ {\text{根芽比}}={\text{根长}}/{\text{芽长}}{\text{,}} $$
    $$\begin{aligned} {\text{相对根长}}= & {\text{不同}}{{\rm{CdCl}}_2}{\text{浓度处理的根长}}/{\text{对应材料}}\\ & {\text{在对照条件下的根长}}{\text{,}} \end{aligned} $$
    $$ \begin{aligned} {\text{相对芽长}}= & {\text{不同}}{{\rm{CdCl}}_2}{\text{浓度处理的芽长}}/{\text{对应材料}}\\ & {\text{在对照条件下的芽长}}{\text{。}} \end{aligned} $$

    P=0.01为差异显著标准,将供试材料的性状值进行单因素方差分析及Duncan’s多重比较,筛选出性状值与‘华粳籼74’差异显著的SSSLs,以这些SSSLs进行耐Cd性QTL鉴定。在本研究中,分别对衡量耐Cd性的相对根长与相对芽长进行QTL鉴定。对有重叠片段的SSSLs通过代换作图[24]的方法确定QTL的最小区间。QTL的命名参考McCouch等[25]的命名原则。参考Eshed等[26]的方法计算QTL的加性效应及表型贡献率:

    $$ {\text{加性效应}}= ({\rm{SSSL}}{\text{表型值}}-{\text{‘华粳籼}}74{\text{’表型值}})/2{\text{,}} $$
    $$ \begin{aligned} {\text{加性效应表型贡献率}}= & {\text{加性效应}}/{\text{‘华粳籼}}74{\text{’}}\\ & {\text{表型值}}\times 100{\text{%}}{\text{。}} \end{aligned} $$

    数据用Excel软件处理,方差分析用SPSS 21.0软件进行单因素方差分析(One-way ANOVA),多重比较采用Duncan’ s字母标记法,QTL在染色体上的分布图用Mapchart2.2软件绘制。

    对‘华粳籼74’、‘新黄占’、‘中花11’、‘台中65’、‘粤晶丝苗2号’、‘象牙香占’6个栽培稻品种的种子进行了5个CdCl2浓度的处理。分析发芽后第7 天发芽率、根长、芽长、根芽比4个性状指标。其中在500 μmol·L−1浓度下,水稻种子几乎不发芽(图1)。在0、30、50和100 μmol·L−1 4个浓度下,6个栽培稻品种的种子发芽率为91%~95%,虽然在P=0.05水平,50和100 μmol·L−1浓度下的发芽率与0 μmol·L−1下的发芽率有显著差异,但30、50和100 μmol·L−1 3个浓度下发芽率没有显著差异(图2)。6个栽培稻品种的根芽伸长生长明显受到CdCl2的抑制。随着CdCl2浓度升高,根长和芽长均显著变短;在不同CdCl2浓度处理下,各品种种子根长较芽长更短,表明种子根的生长受到更大程度的抑制。在对照条件下,6个品种的根芽比接近1,CdCl2处理下的根芽比均小于0.5,比对照显著变小(图2)。总体来看,CdCl2胁迫处理条件下,供试水稻品种种子根、芽的伸长生长均明显受到抑制,相同处理浓度下CdCl2对根的抑制程度大于对芽的抑制程度。

    图  1  部分栽培稻品种种子萌发后第7天的根芽生长状况
    Figure  1.  Root and shoot growth of some cultivated rice varieties on the 7th day after germination
    图  2  栽培稻品种在CdCl2各处理浓度下的发芽结果
    相同性状不同柱子上的小写字母表示不同处理浓度间差异显著(P<0.05,Duncan’ s法)
    Figure  2.  Germination results of cultivated rice varieties at different CdCl2 concentrations
    Different lowercase letters on the columns of the same trait indicate significant differences (P<0.05,Duncan’ s test)

    此外,不同品种的根芽生长对Cd胁迫的反应表现出差异(图3)。首先,正常发芽条件下,发芽第7天供试品种种子根长与芽长有显著差异。‘台中65’和‘新黄占’的根长最长,‘中花11’的根长最短;‘象牙香占’的芽长最长,‘华粳籼74’的芽长最短。其次,在Cd胁迫条件下,在不同CdCl2处理浓度下6个品种种子的根长均比对照短,且品种间差异显著,随着处理浓度升高,根长越来越短(图3A);不同品种间及不同CdCl2处理浓度间芽长的差异较小(图3B)。这说明根长能更好地反映栽培稻品种间耐Cd性差异。就相对长度而言,各处理浓度下各品种的相对芽长均大于相对根长(图4)。

    图  3  栽培稻品种在不同CdCl2处理浓度下的根长与芽长
    相同CdCl2浓度不同柱子上的小写字母表示不同品种间差异显著(P<0.05,Duncan’ s法)
    Figure  3.  Root and shoot length of cultivated rice varieties at different CdCl2concentrations
    Different lowercase letters on the columns of the same CdCl2 concentration indicate significant differences (P<0.05,Duncan’ s test)
    图  4  栽培稻品种在不同CdCl2处理浓度下的相对根长与相对芽长
    Figure  4.  Relative root and shoot length of cultivorted rice varieties at different CdCl2 concentrations

    以上分析结果表明,6个供试品种芽期耐Cd性存在显著差异,籼稻品种(‘华粳籼74’、‘新黄占’、‘粤晶丝苗2号’、‘象牙香占’)较粳稻品种(‘中花11’、‘台中65’)更耐Cd胁迫。就耐Cd性的分析而言,本研究CdCl2处理浓度为50、100 μmol·L−1时,相对根长和相对芽长能够较好反映水稻品种间的耐Cd性差异。因此,在本研究SSSLs耐Cd处理筛选及QTL鉴定时,CdCl2处理浓度采用50和100 μmol·L−1

    对37个SSSLs材料和受体亲本‘华粳籼74’做了多个批次的发芽试验。P=0.01水平的方差分析结果表明,在CdCl2胁迫下,供试材料种子发芽第7天的根长、相对根长、芽长及相对芽长等存在显著差异(数据尚未发表)。以根长与芽长为性状指标时,在P=0.01水平,CdCl2浓度为50和100 μmol·L−1时,供试SSSLs的平均根长分布范围分别为2.48~4.93 cm及1.77~3.03 cm,分别有14个和4个SSSLs的根长显著长于‘华粳籼74’(图5A);供试SSSLs的平均芽长分布范围分别为2.76~3.88 cm及2.58~3.52 cm,50 μmol·L−1处理浓度下有2个SSSLs的芽长显著长于‘华粳籼74’(图5B)。由于本研究中各基因型材料在对照条件下发芽第7天的根长与芽长存在显著差异,因此以各材料的相对根长与相对芽长为指标进行耐Cd SSSLs筛选与QTL鉴定。

    图  5  不同SSSLs在不同CdCl2处理浓度下根长和芽长的差异
    “*”表示供试材料的相关性状值与‘华粳籼74’在P=0.01水平差异显著(单因素方差分析)
    Figure  5.  Differences of root and shoot length in SSSLs at different CdCl2 concentrations
    “*” represents significant difference of relative character between experimental material and ‘Huajingxian 74’ at P=0.01 level (One-way ANOVA)

    P=0.01水平,以相对根长为性状指标时,CdCl2浓度为50 μmol·L−1时筛选出7个相对根长显著长于‘华粳籼74’的SSSLs,分别为A1、A12、A124、A14、A18、A36、A9;CdCl2浓度为100 μmol·L−1时筛选出6个相对根长显著长于‘华粳籼74’的SSSLs,分别为A16、A14、A124、A18、A36、A9(图6A)。以相对芽长为性状指标时,CdCl2浓度为50 μmol·L−1时筛选出5个相对芽长显著长于‘华粳籼74’的SSSLs,分别为A9、A14、A18、A19、A36;CdCl2浓度为100 μmol·L−1时筛选出4个相对芽长显著长于‘华粳籼74’的SSSLs,分别为A9、A14、A18、A19(图6B)。

    图  6  不同SSSLs在不同CdCl2处理浓度下相对根长和相对芽长的差异
    “*”表示供试材料的相关性状值与‘华粳籼74’在P=0.01水平差异显著(单因素方差分析)
    Figure  6.  Differences of relative root and shoot length in SSSLs at different CdCl2 concentrations
    “*” represents significant difference of relative character between experimental material and ‘Huajingxian 74’ at P=0.01 level (One-way ANOVA)

    P=0.01水平、2个Cd处理浓度下,A9、A14、A18、A36及A124相对根长均显著长于‘华粳籼74’,表明其耐Cd性强且稳定。根据SSSLs QTL鉴定的原理,认为以上这些SSSLs的代换片段上携带有来自展颖野生稻的耐Cd QTLs(表1图7)。此外,根据各SSSLs代换片段的染色体分布及重叠关系,A9、A12、A14及A36的代换片段均在1号染色体上,片段重叠区间是RM490-RM572,因此,这4个SSSLs携带相同的1个相对根长QTL qRRL1-1。A18、A16、A1及A124分别在第2、3、3及6号染色体上携带QTLs qRRL2-1qRRL3-1qRRL3-2qRRL6-1。其中qRRL1-1qRRL2-1qRRL6-1在2个CdCl2胁迫浓度均能鉴定到,而qRRL3-1只在100 μmol·L−1浓度出现,qRRL3-2只在50 μmol·L−1浓度出现,说明qRRL1-1qRRL2-1qRRL6-1是相对稳定的QTLs。从各QTL的加性效应来看,在50 μmol·L−1 CdCl2浓度条件下,各QTL的加性效应分布范围为0.19~0.60,表型贡献率为31.18%~100.59%。qRRL1-1的平均加性效应为0.33,表型贡献率为55.92%,为遗传效应最大的QTL。在100 μmol·L−1 CdCl2浓度条件下,各QTL的加性效应及表型贡献率均小于50 μmol·L−1 CdCl2,其中qRRL1-1的平均加性效应为0.25,表型贡献率为50.84%,依然是遗传效应最大的QTL。

    表  1  相对根长及相对芽长QTLs及其遗传效应1)
    Table  1.  Identified QTLs and genetic effects for relative root and shoot length
    性状
    Character
    SSSLs QTL 染色体
    Chromosome
    位置
    Position
    加性效应
    Additive effect
    表型贡献率/%
    Phenotypic contribution rate
    50 μmol·L−1
    CdCl2
    100 μmol·L−1
    CdCl2
    50 μmol·L−1
    CdCl2
    100 μmol·L−1
    CdCl2
    相对根长
    Relative root length
    A36 qRRL1-1 1 PSM012-RM572 0.35 0.25 59.08 50.91
    A14 qRRL1-1 1 RM490-RM572 0.20 0.09 32.86 18.99
    A9 qRRL1-1 1 RM490-RM572 0.60 0.41 100.59 82.63
    A12 qRRL1-1 1 RM490-RM449 0.19 31.18
    A18 qRRL2-1 2 PSM375-end 0.33 0.21 54.71 42.53
    A16 qRRL3-1 3 end-RM489 0.10 20.71
    A1 qRRL3-2 3 RM3280-PSM379 0.22 36.64
    A124 qRRL6-1 6 PSM388-RM454 0.27 0.12 45.80 24.14
    相对芽长
    Relative shoot length
    A19 qRSL1-1 1 end-PSM024 0.08 0.08 8.40 7.98
    A36 qRSL1-2 1 PSM012-RM572 0.07 7.43
    A14 qRSL1-2 1 RM490-RM572 0.14 0.10 15.02 10.27
    A9 qRSL1-2 1 RM490-RM572 0.18 0.14 18.95 14.70
    A18 qRSL2-1 2 PSM375-end 0.12 0.14 12.76 14.29
     1) “—”表示未检测到QTL
     1) “—” represents no QTL detected
    下载: 导出CSV 
    | 显示表格
    图  7  耐Cd QTLs染色体分布
    Figure  7.  Chromosome distribution of identified QTLs for Cd tolerance

    P=0.01水平、CdCl2浓度为50和100 μmol·L−1时,分别筛选出5个(A9、A14、A18、A19、A36)和4个(A9、A14、A18、A19)SSSLs的相对芽长显著长于‘华粳籼74’,本研究认为这些SSSLs携带相对芽长的QTLs,具体见表1图7。A19在1号染色体短臂端end-PSM024区间携带QTL qRSL1-1;A36、A14及A9的代换片段在1号染色体上有重叠区间RM490-RM572,因此这3个SSSLs在重叠区间有1个共同的QTL,为qRSL1-2;A18在2号染色体长臂端PSM375-end区间携带QTL qRSL2-1。对各QTL遗传效应分析表明,3个QTLs加性效应均较小,分布范围为0.08~0.18,表型贡献率为7.43%~18.95%。其中A9携带的qRSL1-2在2个胁迫处理条件下均具有最大的加性效应和表型贡献率,这与其在相对根长中的QTL表现一致。

    土壤Cd污染严重影响水稻的正常生长,而大米又是我国居民Cd摄入的主要途径,约占Cd摄入总量的56%[27],Cd的摄入对人体健康破坏特别大[28],所以水稻耐Cd研究以及培育高耐Cd水稻品种十分关键。目前关于水稻种子萌发及根芽生长过程对Cd胁迫响应的研究报道不多,孟桂元等[3]研究表明Cd胁迫对水稻种子的发芽率、发芽指数影响不显著,对种子活力指数及根芽生长具有显著影响;Cd胁迫对根的抑制作用明显大于对芽的抑制。Cd对水稻生长的影响存在浓度效应和时间效应,对根系的抑制程度大于对苗的抑制程度[29]。本研究以相对根长和相对芽长为指标分析不同基因型栽培稻品种芽期耐Cd性,结果表明,6个供试品种的相对根长在2个CdCl2处理浓度下耐Cd性顺序相同,即:‘台中65’<‘中花11’<‘新黄占’<‘象牙香占’<‘华粳籼74’<‘粤晶丝苗2号’;供试品种相对芽长在2个Cd处理浓度下耐Cd性的变化存在差别。陈毓瑾等[30]对30个常规稻品种在发芽期进行2 mg·L−1 Cd胁迫处理,以发芽率、发芽指数、幼苗根长和幼苗芽长等指标进行耐Cd性研究及低Cd积累种质的筛选,结果表明Cd胁迫对水稻种子的根长具有明显的抑制效应,对种子发芽率、发芽指数、幼苗芽长等指标无显著影响,水稻根长可以作为评价水稻种质耐Cd性的参考指标,水稻品种间根长指数有显著差异,‘IR24’是耐Cd性最强的品种。同时,孙亚莉等[4]以50个水稻品种(系)为材料,研究了不同Cd浓度胁迫对种子发芽率、发芽指数、活力指数、胚芽长、胚根长、胚芽鲜质量与胚芽干质量等性状的影响,表明0.5 mmol·L−1 Cd胁迫浓度对各水稻品种种子的发芽率和发芽指数影响不显著,不同Cd胁迫浓度对水稻种子胚根的抑制作用明显大于胚芽,50个供试品种(系)分为Cd胁迫敏感型、中间型和耐受型等。在水稻芽期耐Cd性指标分析及选择上,本研究结果与陈毓瑾等[30]及孙亚莉等[4]的结果较为一致,说明以相对根长为指标筛选耐Cd种质是可靠的。

    近年来,水稻耐Cd QTL的鉴定与定位研究取得一定进展,主要是苗期耐Cd QTL定位[15-18]及籽粒低Cd积累含量QTL/基因的定位与克隆[6-7, 14, 19-21]等研究,但未见芽期耐Cd QTL定位的报道。在水稻幼苗期耐Cd QTL定位研究中,林辉锋等[15]以‘Lemont’和‘Dular’杂交建立的RIL群体,对水培条件下水稻幼苗耐Cd胁迫的QTL进行了初步定位,利用根长和株高作为抗性指标分别检测到QTLs qRAP-1aqSAP-1b,均位于1号染色体RM243-RM578区间(在染色体上的物理区间:7 971 722~8 411 145 bp),经与本研究1号染色体定位QTL qRRL1-1qRSL1-2的标记区间RM490-RM572 (6 677 153~9 866 595 bp)进行比对,发现RM243-RM578包含在标记区间RM490-RM572内,因此,本研究鉴定的相对根长QTL qRRL1-1与林辉锋等[15]定位的qRAP-1a很可能是同一个QTL,相对芽长QTL qRSL1-2与林辉锋等[15]定位的qSAP-1b很可能是同一个QTL。该区间与李炜星等[18]利用籼稻品种‘昌恢121’与粳稻品种‘越光’构建的RIL群体在幼苗期以根长作为抗性指标所定位到的QTL的位置(1号染色体RM8111处)非常接近,因此,认为在1号染色体RM490-RM572附近应该存在控制水稻耐Cd性的QTL。此外,陈志德等[17]利用‘韭菜青’בIR26’构建的RIL群体在幼苗耐Cd胁迫QTL定位中,在2号染色体长臂端标记RM6312-RM3860 (物理区间:35 390 710~35 431 667 bp) 之间定位了一个控制相对苗高的QTL qRSH-2,该区间包含在本研究鉴定的相对芽长QTL qRSL2-1的区间PSM375 -end(物理区间:33 175 495 bp~35 680 037 bp)内,表明水稻2号染色体该区间存在与水稻地上部耐Cd相关的QTL。综合来看,本研究鉴定到的其余5个QTLs,即qRRL2-1qRRL3-1qRRL3-2qRRL6-1qRSL1-1,所在染色体区间与前人已定位的水稻苗期或芽期耐Cd QTL不存在重叠区域,表明这5个QTLs可能为新的QTLs,且各QTL来自展颖野生稻的等位基因产生了增效作用,增强了SSSL的耐Cd性。水稻不同生育期耐Cd性及其籽粒Cd积累特性极其复杂,目前研究者们对其遗传及生理机制的研究尚在进行。对本文鉴定的QTLs进一步进行精细定位与克隆将有助于揭示展颖野生稻SSSLs芽期耐Cd性的机理;此外,对携带芽期耐Cd QTL的SSSLs进一步开展苗期耐Cd性以及籽粒Cd积累的研究也有助于探索QTL的多效性,为开展水稻耐Cd性的分子育种奠定基础。

  • 图  1   静原鸡不同组织总RNA提取

    1、2:胸肌;3、4:腿肌;M:Trans2K Plus DNA Marker

    Figure  1.   Total RNA extraction from different tissues of Jingyuan chicken

    1, 2: Chest muscle; 3, 4: Leg muscle; M: Trans2K Plus DNA Marker

    图  2   静原鸡AK1基因PCR扩增电泳图

    M: DL2000 DNA Marker;1~5:扩增片段

    Figure  2.   Electrophoresis of PCR amplification of AK1 gene in Jingyuan chicken

    M: DL2000 DNA Marker; 1–5: Amplified fragments

    图  3   静原鸡AK1基因加酶切位点PCR扩增电泳图

    M:DL2000 DNA Marker;1~5:扩增片段

    Figure  3.   Electrophoresis of PCR amplification of AK1 gene with restriction site in Jingyuan chicken

    M: DL2000 DNA Marker; 1–5: Amplified fragments

    图  4   静原鸡AK1-pEGFP-N1载体酶切电泳图

    M:DL2000 DNA Marker;1~7:扩增片段

    Figure  4.   Restrictive enzyme digestion and electrophoresis of AK1-pEGFP-N1 vector in Jingyuan chicken

    M: DL2000 DNA Marker; 1–7: Amplified fragments

    图  5   AK1基因编码区全长序列与对应的氨基酸序列

    Figure  5.   The full length coding sequence of AK1 gene and corresponding amino acid sequence

    图  6   AK1基因CDS区核苷酸组成

    Figure  6.   Nucleotide composition in CDS region of AK1 gene

    图  7   静原鸡AK1蛋白的疏水性/亲水性分析

    Figure  7.   Hydrophobic/hydrophilic analysis for AK1 protein in Jingyuan chicken

    图  8   静原鸡AK1蛋白磷酸化位点

    Figure  8.   The phosphorylation sites of AK1 protein in Jingyuan chicken

    图  9   AK1蛋白B抗原表位预测

    Figure  9.   Prediction of B epitope of AK1 protein

    图  10   静原鸡AK1蛋白质二级结构

    Figure  10.   The secondary structure of AK1 protein in Jingyuan chicken

    图  11   静原鸡AK1蛋白质预测三级结构

    Figure  11.   The predicted tertiary structure of AK1 protein in Jingyuan chicken

    表  1   静原鸡AK1基因CDS区引物序列

    Table  1   Primer sequence for CDS region in AK1 gene of Jingyuan chicken

    基因
    Gene
    登录号
    Accession number
    引物序列(5′→3′)1)
    Primer sequence
    产物长度/bp
    Product length
    退火温度/℃
    Annealing temperature
    AK1 M37901 F:TCCTCCACCCAGACAGCA
    R:ACGGGAAAGAGCCAAACA
    741 57.6
    AK13 M37901 F:ccg GAATTCTCCTCCACCCAGACAGCA
    R:cgc GGATCCACGGGAAAGAGCCAAACA
    759 62.7
     1) 带有下划线的序列为酶切位点,AK13上、下游引物分别插入EcoRⅠ和BamHⅠ酶切位点
     1) Underlined sequences are restriction sites, and forward and reverse primers of AK13 contained EcoR Ⅰ and BamHⅠ restriction sites respectively
    下载: 导出CSV

    表  2   AK1基因编码蛋白的氨基酸组成

    Table  2   The composition of amino acid coded by AK1 gene

    氨基酸
    Amino Acid
    数量
    Number
    频率/%
    Frequency
    Ala 10 5.2
    Arg 10 5.2
    Asn 2 1.0
    Asp 11 5.7
    Cys 2 1.0
    Gln 6 3.1
    Glu 18 9.3
    Pro 6 3.1
    Thr 13 6.7
    Val 15 7.7
    Gly 19 9.8
    His 4 2.1
    Ile 9 4.6
    Leu 20 10.3
    Lys 22 11.3
    Met 5 2.6
    Phe 5 2.6
    Ser 9 4.6
    Tyr 8 4.1
    下载: 导出CSV

    表  3   AK1基因的GO注释分析

    Table  3   GO annotation analysis of AK1 gene

    项目
    Item
    GO条目
    GO term
    描述
    Description
    基因数1)
    Gene number
    错误发现率
    False discovery rate
    生物学过程
    Biological process
    GO:0009167 嘌呤核糖核苷单磷酸代谢过程
    Purine ribonucleoside monophosphate metabolic process
    4(41) 2.78×10−6
    GO:0009165 核苷酸生物合成的过程
    Nucleotide biosynthetic process
    4(40) 2.78×10−6
    GO:0009150 嘌呤核糖核苷酸的代谢过程
    Purine ribonucleotide metabolic process
    4(49) 2.78×10−6
    GO:0009168 嘌呤核糖核苷单磷酸生物合成过程
    Purine ribonucleoside monophosphate biosynthetic process
    3(25) 7.56×10−6
    GO:0009152 嘌呤核糖核苷酸的生物合成过程
    Purine ribonucleotide biosynthetic process
    3(29) 8.98×10−6
    GO:0017144 药物代谢过程
    Drug metabolic process
    3(68) 7.47×10−5
    GO:0006188 肌苷酸生物合成过程
    IMP biosynthetic process
    2(7) 8.04×10−5
    GO:0046040 肌苷酸代谢过程
    IMP metabolic process
    2(9) 1.20×10−4
    GO:0006165 二磷酸核苷磷酸化
    Nucleoside diphosphate phosphorylation
    2(14) 2.40×10−4
    GO:0009142 核苷三磷酸生物合成过程
    Nucleoside triphosphate biosynthetic process
    2(21) 4.20×10−4
    GO:0046034 ATP代谢过程
    ATP metabolic process
    2(31) 8.40×10−4
    分子功能
    Molecular function
    GO:0003824 催化活性
    Catalytic activity
    5(642) 2.00×10−3
    GO:0016301 激酶活性
    Kinase activity
    2(98) 7.50×10−3
    GO:0005524 ATP结合
    ATP binding
    3(202) 7.50×10−3
    GO:0016787 水解酶活性
    Hydrolase activity
    2(271) 3.13×10−2
    细胞组分
    Cellular component
    GO:0005829 细胞溶质
    Cytosol
    3(264) 1.43×10−2
     1) 括号内的数值表示富集到这个条目的所有基因数,括号外的数值表示差异显著的基因数
     1) The value in brackets represents the number of all genes enriched to this GO term, and the value outside brakets represents the number of significantly different genes
    下载: 导出CSV
  • [1] 王建保. 朝那鸡列入区级畜禽遗传保护品种[J]. 固原日报, 2009(1): 9-25.
    [2] 张娟, 母童, 赵平, 等. 静原鸡ELOVL5基因遗传多样性研究[J]. 浙江农业学报, 2019, 31(2): 200-206. doi: 10.3969/j.issn.1004-1524.2019.02.04
    [3] 母童, 张娟, 赵平, 等. 静原鸡ELOVL2ELOVL5基因表达的组织特异性研究[J]. 浙江农业学报, 2017, 29(8): 1290-1296. doi: 10.3969/j.issn.1004-1524.2017.08.09
    [4]

    JONES M E, LIPMANN F. Aceto-CoA-kinase: ATP + CoA + Ac ⇄Ac~S CoA + 5-AMP + PP[J]. Methods in Enzymology, 1955, 1: 585-591. doi: 10.1016/0076-6879(55)01101-4

    [5]

    ZELEZNIKAR R J, HEYMAN R A, GRAEFF R M, et al. Evidence for compartmentalized adenylate kinase catalysis serving a high energy phosphoryl transfer function in rat skeletal muscle[J]. Journal of Biological Chemistry, 1990, 265: 300-311.

    [6]

    VAN ROMPAY A R, JOHANSSON M, KARLSSON A. Phosphorylation of nucleosides and nucleoside analogs by mammalian nucleoside monophosphate kinases[J]. Pharmacolocy & Therapeutics, 2000, 87(2/3): 189-198.

    [7]

    TANABE T, YAMADA M, NOMA T, et al. Tissue- specific and developmentally regulated expression of the genes encoding adenylate kinase isozymes[J]. International Journal of Biochemistry & Cell Biology, 1993, 113(2): 200-207.

    [8]

    SCHULZ G E, ELZINGA M, MARX F, et al. Three-dimensional structure of adenylkinase[J]. Nature, 1974, 250(5462): 120-123. doi: 10.1038/250120a0

    [9]

    SACHSENHEIMER W, SCHULZ G E. Two conformations of crystalline adenylatekinase[J]. Journal of Molecular Biology, 1977, 114(1): 23-36. doi: 10.1016/0022-2836(77)90280-7

    [10]

    HANCOCK C R, JANSSEN E, TERJUNG R L. Skeletal muscle contractile performance and ADP accumulation in adenylate kinase-deficient mice[J]. American Journal of Physiology-Cell Physiology, 2005, 288(6): C1287-C1297. doi: 10.1152/ajpcell.00567.2004

    [11]

    CARRASCO A J, DZEJA P P, ALEKSEEV A E, et al. Adenylate kinase phosphotransfer communicates cellular energetic signals to ATP-sensitive potassium channels[J]. Proceedings of the National Academy of Sciences of the United States of America, 2001, 98(13): 7623-7628. doi: 10.1073/pnas.121038198

    [12] 易继财. 生物类专业生物信息学课程教学探索: 华南农业大学生物类专业生物信息学课程的教改实践与思考[J]. 安徽农业科学, 2018, 46(26): 231-233. doi: 10.3969/j.issn.0517-6611.2018.26.071
    [13] 吴风瑞, 陈德宇, 姜双林, 等. PBL教学法在硕士研究生生物信息学教学中的应用[J]. 阜阳师范学院学报(自然科学版), 2017, 34(2): 109-111.
    [14]

    PANAYIOTOU C, SOLAROLI N, KARLSSON A. The many isoforms of human adenylate kinases[J]. International Journal of Biochemistry & Cell Biology, 2014, 49: 75-83.

    [15]

    JACOBASCH G. Biochemical and genetic basis of red cell enzyme deficiencies[J]. Best Practice & Research Clinical Haematology, 2000, 13(1): 1-20.

    [16]

    ABRUSCI P, CHIARELLI L R, GALIZZI A, et al. Erythrocyte adenylate kinase deficiency: Characterization of recombinant mutant forms and relationship with nonspherocytic hemolytic anemia[J]. Experimental Hematology, 2007, 35(8): 1182-1189. doi: 10.1016/j.exphem.2007.05.004

    [17]

    PUCAR D, BAST P, GUMINA R J, et al. Adenylate kinase AK1 knockout heart: Energetics and functional performance under ischemia-reperfusion[J]. American Journal of Physiology Heart and Circulatory Physiology, 2002, 283(2): H776-H782. doi: 10.1152/ajpheart.00116.2002

    [18]

    SUMINAMI Y, KISHI F, TORIGOE T, et al. Structure and complete nucleotide sequence of the gene encoding chicken cytosolic adenylate kinase[J]. Journal of Biochemistry, 1988, 103(4): 611-617. doi: 10.1093/oxfordjournals.jbchem.a122315

    [19]

    PAI E F, SACHSENHEIMER W, SCHIRMER R H, et al. Substrate positions and induced-fit in crystalline adenylate kinase[J]. Journal of Molecular Biology, 1977, 114(1): 37-45. doi: 10.1016/0022-2836(77)90281-9

    [20]

    FRY D C, KUBY S A, MILDVAN A S. ATP-binding site of adenylate kinase: Mechanistic implications of its homology with ras-encoded p21, F1-ATPase, and other nucleotide-binding proteins[J]. Proceedings of the National Academy of Sciences of the United States of America, 1986, 83(4): 907-911. doi: 10.1073/pnas.83.4.907

    [21]

    GURUPRASAD K, REDDY B V, PANDIT M W. Correlation between stability of a protein and its dipeptide composition: A novel approach for predicting in vivo stability of a protein from its primary sequence[J]. Protein Engineering, 1990, 4(2): 155-161. doi: 10.1093/protein/4.2.155

    [22]

    ZOR T, MAYR B M, DYSON H J, et al. Roles of phosphorylation and helix propensity in the binding of the KIX domain of CREB-binding protein by constitutive (c-myb) and inducible (CREB) activators[J]. Journal of Biological Chemistry, 2002, 277(44): 42241-42248. doi: 10.1074/jbc.M207361200

    [23]

    KYTE J, DOOLITTLE R F. A simple method for displaying the hydropathic character of a protein[J]. Journal of Molecular Biology, 1982, 157(1): 105-132. doi: 10.1016/0022-2836(82)90515-0

    [24] 郑海军, 朱荣, 葛春蕾, 等. 人白细胞介素−29的生物信息学分析[J]. 中国生物制品学杂志, 2013, 26(2): 209-212.
    [25] 赵骁, 张峰波, 王红英, 等. 生物信息分析细粒棘球绦虫EgA31蛋白T细胞及B细胞的优势抗原表位[J]. 中国组织工程研究, 2019, 23(7): 1078-1083. doi: 10.3969/j.issn.2095-4344.1001
    [26]

    CHOO H J, KIM B W, KWON O B, et al. Secretion of adenylate kinase 1 is required for extracellular ATP synthesis in C2C12 myotubes[J]. Experimental and Molecular Medicine, 2008, 40(2): 220-228. doi: 10.3858/emm.2008.40.2.220

  • 期刊类型引用(1)

    1. 赵薇,曹国伟,周子航,王卫振,辛国省,蔡正云,顾亚玲,张娟. 鸡肉品质相关调控基因研究进展. 农业生物技术学报. 2022(02): 370-378 . 百度学术

    其他类型引用(1)

图(11)  /  表(3)
计量
  • 文章访问数:  942
  • HTML全文浏览量:  12
  • PDF下载量:  1116
  • 被引次数: 2
出版历程
  • 收稿日期:  2020-06-17
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2021-03-09

目录

/

返回文章
返回