欢迎访问作物学报，今天是
作物学报

作物学报 ›› 2022, Vol. 48 ›› Issue (11): 2765-2773.doi: 10.3724/SP.J.1006.2022.14143

• • 上一篇    下一篇

ChIP-seq分析花生中AhGLK1与AhHDA1调控的下游靶基因网络

刘星1(), 苏良辰1, 李丽梅2, 李玲2,*()   

  1. 1遵义医科大学珠海校区生物工程系, 广东珠海 519041
    2华南师范大学生命科学学院 / 广东省植物发育生物工程重点实验室, 广东广州510631
  • 收稿日期:2021-08-10 接受日期:2021-11-29 出版日期:2022-11-12 网络出版日期:2021-12-17
  • 通讯作者: 李玲
  • 作者简介:第一作者联系方式: E-mail: liuxingas@126.com
  • 基金资助:
    本研究由贵州省科学技术厅基金项目(ZK2021-127);国家自然科学基金项目(31671600)

Downstream target gene network regulated by AhGLK1 and AhHDA1 using ChIP-seq in peanut

LIU Xing1(), SU Liang-Chen1, LI Li-Mei2, LI Ling2,*()   

  1. 1Department of Bioengineening, Zhuhai Campus, Zunyi Medical University, Zhuhai 519041, Guangdong, China
    2Guangdong Provincial Key Laboratory of Biotechnology for Plant Development / School of Life Sciences, South China Normal University, Guangzhou 510631, Guangdong, China
  • Received:2021-08-10 Accepted:2021-11-29 Published:2022-11-12 Published online:2021-12-17
  • Contact: LI Ling
  • Supported by:
    The Department of Science and Technology of Guizhou Province(ZK2021-127);The National Natural Science Foundation of China(31671600)

RichHTML

31

PDF

498

摘要：

花生(Arachis hypogaea L.)是重要的经济油料作物, 其生长发育、产量与品质受干旱影响。为深入了解花生的抗旱机理, 本研究通过ChIP-seq对组蛋白去乙酰化酶AhHDA1和转录因子AhGLK1富集的DNA序列进行分析, 揭示两者调控的下游靶基因网络。通过比对分析, GLK-IP获得6571万clean beads, HDA-IP获得6390万clean beads, Input获得7006万clean beads, 唯一比对率分别为74.97%、76.81%和76.75%。GLK-IP获得714个peak, HDA-IP获得543个peak。Peak在基因的外显子、内含子、上游、下游和基因间等功能元件均有分布。GO富集结果显示, AhGLK1-IP和AhHDA1-IP的peak相关基因在分子功能中的富集分别为35.1%和32.8%, 在生物学过程中的富集分别为39.3%和44.2%, 在细胞组分中的富集分别为25.5%和22.8%。KEGG信号通路富集结果显示, AhGLK1-IP相关基因显著富集在“代谢途径(metabolic pathways)”、“抗生素生物合成(biosynthesis of antibiotics)”、“二羧酸代谢(glyoxylate and dicarboxylate metabolism)”、“不同环境中微生物代谢(microbial metabolism in diverse environments)”、“碳代谢(carbon metabolism)”、“次生代谢生物合成(biosynthesis of secondary metabolites)”和“氨基酸生物合成(biosynthesis of amino acids)。而AhHDA1-IP相关基因在“N聚糖生物合成(N-glycan biosynthesis)”、“精氨酸和脯氨酸代谢(arginine and proline metabolism)”和“苯丙氨酸代谢(phenylalanine metabolism)”通路显著富集。AhGLK1-IP和AhHDA1-IP共同富集的peak有4个, 在AhGLK1-IP和AhHDA1-IP特异富集的基序(motif)中存在共同的保守序列AGAA/T。研究结果为深入认识AhGLK1和AhHDA1基因的功能和了解花生响应干旱胁迫和旱后恢复生长中的调控机制具有参考价值。

关键词: 花生, ChIP-seq, AhGLK1, AhHDA1, 靶基因

Abstract:

Peanut (Arachis hypogaea L.) is an important economic and edible oil crop whose growth, yield and quality are seriously affected by drought stress. In order to further investigate the mechanism of drought tolerance in peanut, chromatin immunoprecipitation sequence (ChIP-seq) was conducted to analyze immunoaffinity separation of specific binding DNA fragments with AhHDA1 (Arachis hypogaea L. histone deacetylase 1) and transcription factor AhGLK1 (Arachis hypogaea L. Golden2-like 1) for sequence identification, revealing the downstream target gene networks regulated by AhHDA1 and AhGLK1. Through alignment and analysis, 65.71, 63.90, and 70.06 million clean beads were obtained for GLKHDA and input IP from the sequencing run, respectively. The unique alignment of reads to the reference genome sequence were 74.97%, 76.81%, and 76.75%, respectively. GLK-IP obtained 714 peaks and HDA-IP obtained 543 peaks. Peaks were distributed into exons, introns, upstream, downstream, intergenic and other functional elements. GO enrichment revealed that 35.1% and 32.8% of the peak-related genes of AhGLK1-IP and AhHDA1-IP were enriched in “molecular function”. 39.3% and 44.2% were enriched in “biological process” and 25.5% and 22.8% were enriched in “cellular components”, respectively. KEGG pathway enrichment showed that AhGLK1-IP related genes were significantly enriched in metabolic, biosynthesis of antibiotics, glyoxylate and dicarboxylate metabolism, microbial metabolism in diverse environments, carbon metabolism, biosynthesis of secondary metabolites and biosynthesis of amino acids pathways. Whereas, AhHDA1-IP related genes were significantly enriched in N-glycan biosynthesis, arginine, proline and phenylalanine metabolism pathways. Moreover, there were four common peaks enriched both in AhGLK1-IP and AhHDA1-IP. A common conserved motif sequences (AGAA/T) was presented in the motif sequences specifically enriched by AhGLK1-IP and AhHDA1-IP. These results of the study have some reference value for further understanding the functions of AhGLK1 and AhHDA1 genes and the regulatory mechanism of peanuts in response to drought stress and post-drought recovery growth.

Key words: Arachis hypogaea L., chromatin immunoprecipitation sequence, AhGLK1, AhHDA1, target gene

表1

ChIP-seq数据产出统计"

样品名
Sample name		 序列总数
No. of total sequences		 碱基总数
No. of total bases		 GC
(%)		 Q20
(%)		 Q30
(%)
GLK-IP 65,712,234 9,856,835,100 41.50 97.60 93.96
HDA-IP 63,903,912 9,85,586,800 40.73 97.61 93.96
Input 70,062,804 10,509,420,600 40.49 97.43 93.58

表2

数据比对结果统计"

样品名
Sample		 总序列数
No. of total
sequences		 比对序列数
No. of aligned
sequences		 比对率
Aligned rate
(%)		 唯一比对序列数
No. of unique aligned
sequences		 唯一比对率
Unique aligned rate
(%)
GLK-IP 65,133,710 60,248,717 92.50 48,828,873 74.97
HDA-IP 63,362,432 59,643,659 94.13 48,666,448 76.81
Input 69,408,066 68,607,870 98.85 53,268,396 76.75

表3

Peak信息统计"

样品名
Sample name		 Peak数
Peak number		 Peak总长度
Peak total
length		 Peak平均长度
Peak average length		 Peak总序列深度
Peak total sequence depth		 Peak平均序列度
Peak average sequence degree		 基因组比例
Genome rate
(%)
GLK-IP 714 416,446 583 11,224 15 0.02
HDA-IP 543 298,118 549 7764 14 0.01

图1

AhGLK1和AhHDA1富集基因的peak分布"

图2

AhGLK1和AhHDA1富集基因GO分类分析"

图3

AhGLK1和AhHDA1富集基因KEGG pathway分析 A: AhGLK1-IP 相关基因; B: AhHDA1-IP相关基因。"

表4

AhGLK1和AhHDA1-IP相关基因在不同KEGG通路中的分布"

样品
Sample name		 KEGG通路
KEGG pathway		 基因个数
Gene number		 基因ID
Gene ID
AhGLK-IP 代谢途径
Metabolic pathways		 20 Arahy.14YEDZ, Arahy.46RTPJ, Arahy.46W7S4, Arahy.6Q1KS5, Arahy.8T1UF6, Arahy.C0Q7J8, Arahy.CQFZ3U, Arahy.D1ZPVC, Arahy.IC3HGI, Arahy.LT225G, Arahy.MY1TBE, Arahy.NDL7Z2, Arahy.S45BII, Arahy.S8DJY1, Arahy.SKV00U, Arahy.SP7CPF, Arahy.U4VZKL, Arahy.VD5UEC, Arahy.YS6PIM, Arahy.ZZNJ3G
抗生素生物合成
Biosynthesis of antibiotics		 8 Arahy.46W7S4, Arahy.8T1UF6, Arahy.C0Q7J8, Arahy.CQFZ3U, Arahy.NDL7Z2, Arahy.S45BII, Arahy.SP7CPF, Arahy.U4VZKL
二羧酸代谢
Glyoxylate and dicarboxylate metabolism		 3 Arahy.NDL7Z2, Arahy.S45BII, Arahy.VD5UEC
不同环境中微生物代谢
Microbial metabolism in diverse environments		 6 Arahy.CQFZ3U, Arahy.NDL7Z2, Arahy.S45BII, Arahy.SP7CPF, Arahy.U4VZKL, Arahy.VD5UEC
碳代谢
Carbon metabolism		 5 Arahy.CQFZ3U, Arahy.NDL7Z2, Arahy.S45BII, Arahy.U4VZKL, Arahy.VD5UEC
次生代谢生物合成
Biosynthesis of secondary metabolites		 12 Arahy.46W7S4, Arahy.6Q1KS5, Arahy.8T1UF6, Arahy.C0Q7J8, Arahy.CQFZ3U, Arahy.IC3HGI, Arahy.NDL7Z2, Arahy.S45BII, Arahy.SP7CPF, Arahy.U4VZKL, Arahy.YS6PIM, Arahy.ZZNJ3G
AhHDA-IP N聚糖生物合成
N-Glycan biosynthesis		 3 Arahy.6B4W14,Arahy.7LT28F, Arahy.YXCZ83
精氨酸和脯氨酸代谢
Arginine and proline metabolism		 3 Arahy.EJ6Q9X, Arahy.GRS0A1, Arahy.UIG78U
苯丙氨酸代谢
Phenylalanine metabolism		 2 Arahy.21PEZK, Arahy.SGZ2CH

图4

AhGLK1和AhHDA1共同富集peak分布"

图5

AhGLK1和AhHDA1富集基序序列"

[1] Luo H, Guo J, Ren X, Chen W, Huang L, Zhou X, Chen Y, Liu N, Xiong F, Lei Y, Liao B, Jiang H. Chromosomes A07 and A05 associated with stable and major QTLs for pod weight and size in cultivated peanut (Arachis hypogaea L.). Theor Appl Genet, 2018, 131: 267-282.
doi: 10.1007/s00122-017-3000-7
[2] 邓斌, 李玲, 李晓云, 苏良辰, 崔威韬, 冯锦欣. AhHDA1异源表达影响拟南芥植株干旱性. 华南师范大学学报(自然科学版), 2016, 48(5): 52-57.
Deng B, Li L, Li X Y, Su L C, Cui W T, Feng J X. Heterologous expression of AhHDA1 affects drought resistance of Arabidopsis plants. J South Chin Nor (Nat Sci Edn), 2016, 48(5): 52-57. (in Chinese with English abstract)
[3] Su L, Liu S, Liu X, Zhang B, Li M, Zeng L, Li L. Transcriptome profiling reveals histone deacetylase 1 gene overexpression improves flavonoid, isoflavonoid, and phenylpropanoid metabolism in Arachis hypogaea hairy roots. PeerJ, 2021, 9: e10976.
doi: 10.7717/peerj.10976
[4] 苏良辰. AhHDA1在花生响应干旱和ABA信号中的分子机制研究. 华南师范大学博士学位论文, 广东广州, 2017.
Su L C. Study on the Molecular Mechanism of AhHDA1 Responding to Drought and ABA Signal in Arachis hypogaea L. PhD Dissertation of South China Normal University, Guangzhou, Guangdong, China, 2017. (in Chinese with English abstract)
[5] 李媚娟, 苏良辰, 刘帅, 李晓云, 李玲. 花生AhHDA1互作蛋白AhGLK的筛选及特性分析. 作物学报, 2017, 43: 218-225.
doi: 10.3724/SP.J.1006.2017.00218
Li M J, Su L C, Liu S, Li X Y, Li L. Screening of AhHDA1 interacting-protein AhGLK and characterization in peanut (Arachis hypogaea L.). Acta Agron Sin, 2017, 43: 218-225. (in Chinese with English abstract)
doi: 10.3724/SP.J.1006.2017.00218
[6] Liu X, Li L, Li M, Su L, Lian S, Zhang B, Li X, Ge K, Li L. AhGLK1 affects chlorophyll biosynthesis and photosynthesis in peanut leaves during recovery from drought. Sci Rep, 2018, 8: 2250.
doi: 10.1038/s41598-018-20542-7
[7] Liu X, Li L, Zhang B, Zeng L, Li L. AhHDA1-mediated AhGLK1 promoted chlorophyll synthesis and photosynthesis regulates recovery growth of peanut leaves after water stress. Plant Sci, 2020, 294: 110461.
doi: 10.1016/j.plantsci.2020.110461
[8] 蔡大伟. 基于ChIP-seq全基因组识别毛竹笋尖与鞭笋尖组蛋白修饰位点. 福建农林大学硕士学位论文, 福建福州, 2017.
Cai D W. Genome-wide Mapping of Histone Modifications of Shoot Tips and Rhizome Shoot Tips in Moso bamboo by ChIP-seq. MS Thesis of FuJian Agriculture and Forestry University, Fuzhou, Fujian, China, 2017. (in Chinese with English abstract)
[9] Lismer A, Lambrot R, Lafleur C, Dumeaux V, Kimmins S. ChIP-seq protocol for sperm cells and embryos to assess environmental impacts and epigenetic inheritance. STAR Protoc, 2021, 2: 100602.
doi: 10.1016/j.xpro.2021.100602
[10] Liu S, Li M, Su L, Ge K, Li L, Li X, Liu X, Li L. Negative feedback regulation of ABA biosynthesis in peanut (Arachis hypogaea): a transcription factor complex inhibits AhNCED1 expression during water stress. Sci Rep, 2016, 6: 37943.
doi: 10.1038/srep37943
[11] Trapnell C, Williams B A, Pertea G, Mortazavi A, Kwan G, van Baren M J, Salzberg S L, Wold B J, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol, 2010, 28: 511-515.
doi: 10.1038/nbt.1621 pmid: 20436464
[12] Emes M J, Tetlow I J, Bowsher C G. Integration of metabolism within non-photosynthetic plastids, and with the cytosol. Regul Primary Metabolic Pathways Plants, 1999, 42: 117-136.
[13] 苏良辰, 钟钰婷, 李玲. AhHDA1对花生毛状根干旱的生理调节作用. 植物生理学报, 2017, 53: 1893-1900.
Su L C, Zhong Y T, Li L. Physiological regulation effect of AhHDA1 on peanut hairy roots to drought stress. PlantPhysiol J, 2017, 53: 1893-1900. (in Chinese with English abstract)
[14] Waters M T, Wang P, Korkaric M, Capper R G, Saunders N J, Langdale J A. GLK transcription factors coordinate expression of the photosynthetic apparatus in Arabidopsis. Plant Cell, 2009, 21: 1109-1128.
doi: 10.1105/tpc.108.065250 pmid: 19376934
[15] Savitch L V, Subramaniam R, Allard G C, Singh J. The GLK1 ‘regulon’ encodes disease defense related proteins and confers resistance to Fusarium graminearum in Arabidopsis. BiochemBiophys Res Commun, 2007, 359: 234-238.
doi: 10.1016/j.bbrc.2007.05.084
[16] Zubo Y O, Blakley I C, Franco-Zorrilla J M, Yamburenko M V, Solano R, Kieber J J, Loraine A E, Schaller G E. Coordination of chloroplast development through the action of the GNC and GLK transcription factor families. Plant Physiol, 2018, 178: 130-147.
doi: 10.1104/pp.18.00414 pmid: 30002259
[1] 郑玉珍, 齐飞艳, 孙子淇, 刘华, 秦利, 石磊, 王娟, 汪蒙蒙, 韩锁义, 徐静, 苗利娟, 黄冰艳, 董文召, 郑峥, 张新友. 花生籽仁总超长链脂肪酸和7种脂肪酸组分的QTL定位[J]. 作物学报, 2026, 52(6): 1646-1657.
[2] 陆一楚, 李振影, 买春海, 赵晓蕊, 王利祥. 夜间复合体关键基因AhLUX1正向调控花生结瘤的功能研究[J]. 作物学报, 2026, 52(6): 1658-1668.
[3] 谷春苗, 王润风, 黄璐, 刘浩, 鲁清, 李海芬, 李少雄, 何双呈, 洪彦彬, 陈小平, 谭斌, 余倩霞. 花生WOX基因家族的全基因组分析及不定芽再生候选基因的鉴定[J]. 作物学报, 2026, 52(5): 1326-1340.
[4] 于天一, 王春晓, 肖丽, 钟召迪, 王宣仓, 赵勇, 路亚, 吴月, 吴正锋. 不同结瘤特性花生品种氮素累积、产量及品质特性对氮肥用量的响应[J]. 作物学报, 2026, 52(3): 881-894.
[5] 张胜忠, 李国卫, 戈立江, 王菲菲, 胡晓辉, 苗华荣, 李燕, 钟文, 陈静. 花生机械脱壳损伤相关农艺指标筛选与QTL定位[J]. 作物学报, 2026, 52(2): 644-652.
[6] 王菲菲, 张胜忠, 杨贵华, 苗华荣, 胡晓辉, 张则林, 刘莎莎, 乔利仙, 单世华, 陈静. 331份花生种质苗期耐盐性综合评价和强耐盐种质鉴选[J]. 作物学报, 2026, 52(1): 279-294.
[7] 金欣欣, 宋亚辉, 苏俏, 杨永庆, 王瑾. 高产高油高油酸花生品种的生长发育及干物质生产特征[J]. 作物学报, 2026, 52(1): 191-201.
[8] 迟晓元, 刘庆, 张君, 赵旭红, 李美, 于天一, 潘丽娟, 许静, 姜骁, 殷祥贞, 马俊卿, 陈娜. 不同花生品种(系)耐盐碱性田间鉴定及各性状指标相关性研究[J]. 作物学报, 2026, 52(1): 85-98.
[9] 孙辰硕, 张月, 田泽锴, 晏立英, 康彦平, 陈玉宁, 王欣, 淮东欣, 王前前, 姜慧芳, 罗怀勇, 黄莉, 廖伯寿, 王志慧, 雷永. 花生种质果柄强度的遗传分化与主要影响因子分析[J]. 作物学报, 2026, 52(1): 118-130.
[10] 万书波, 张佳蕾, 高华鑫, 王才斌. 中国花生高产栽培研究进展与展望[J]. 作物学报, 2025, 51(7): 1703-1711.
[11] 郭腾达, 崔梦杰, 陈琳杰, 韩锁义, 郭敬坤, 吴晨迪, 付留洋, 黄冰艳, 董文召, 张新友. 花生磷脂酰肌醇转运蛋白基因AhSFH的克隆及其响应黄曲霉菌侵染的表达特征分析[J]. 作物学报, 2025, 51(6): 1489-1500.
[12] 李文佳, 廖泳俊, 黄璐, 鲁清, 李少雄, 陈小平, 金晶炜, 王润风. 花生开花时间的全基因组关联分析及候选基因筛选[J]. 作物学报, 2025, 51(5): 1400-1408.
[13] 林伟津, 郭泽佳, 刘浩, 李海芬, 王润风, 黄璐, 余倩霞, 陈小平, 洪彦彬, 李少雄, 鲁清. 花生荚果产量相关性状QTL定位与候选基因分析[J]. 作物学报, 2025, 51(4): 969-981.
[14] 方应浩, 周波, 陈茹梅, 杨文竹, 秦慧民. 基于RNA-seq和PER-seq联合分析探究ZmHDZ6表达调控网络[J]. 作物学报, 2025, 51(4): 958-968.
[15] 迟晓元, 毕竞男, 赵健鑫, 陈娜, 潘丽娟, 姜骁, 殷祥贞, 赵旭红, 马俊卿, 许静. 花生荚果力学特性评鉴及早熟种质筛选[J]. 作物学报, 2025, 51(4): 943-957.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
京ICP备15008789号-2