中国农业科学 ›› 2024, Vol. 57 ›› Issue (6): 1034-1048.doi: 10.3864/j.issn.0578-1752.2024.06.002

• 作物遗传育种·种质资源·分子遗传学 • 上一篇    下一篇

干旱条件下玉米转座子插入关联的表观调控分析

高晨曦(), 郝陆洋, 胡悦, 李永祥, 张登峰, 李春辉, 宋燕春, 石云素, 王天宇, 黎裕(), 刘旭洋()   

  1. 中国农业科学院作物科学研究所/作物基因资源与育种全国重点实验室,北京 100081
  • 收稿日期:2023-08-30 接受日期:2023-10-23 出版日期:2024-03-25 发布日期:2024-03-25
  • 通信作者:
    黎裕,E-mail:
    刘旭洋,E-mail:
  • 联系方式: 高晨曦,E-mail:18745115576@163.com。
  • 基金资助:
    国家重点研发计划(2020YFE0202300)

Analysis of Transposable Element Associated Epigenetic Regulation under Drought in Maize

GAO ChenXi(), HAO LuYang, HU Yue, LI YongXiang, ZHANG DengFeng, LI ChunHui, SONG YanChun, SHI YunSu, WANG TianYu, LI Yu(), LIU XuYang()   

  1. Institute of Crop Sciences, Chinese Academy of Agricultural Sciences/State Key Laboratory of Crop Gene Resources and Breeding, Beijing 100081
  • Received:2023-08-30 Accepted:2023-10-23 Published:2024-03-25 Online:2024-03-25

摘要:

【目的】 干旱是全球范围影响玉米生产的最主要胁迫因素之一。解析抗旱性的遗传基础与分子机制为玉米的抗旱改良提供依据。【方法】 利用代表性玉米自交系,以叶片相对含水量和散粉-吐丝间隔为指标开展田间抗旱性精准鉴定。筛选2个抗旱性极端差异的自交系,开展基因组重测序分析和转座子插入鉴定;利用全基因组重亚硫酸盐测序(WGBS)方法分析不同水分处理下叶片和根系组织的DNA甲基化水平;同时利用转录组测序方法对相同样品的基因表达进行分析;通过比较分析获得2个材料间的转座子插入缺失变异、差异甲基化区域和差异表达基因,并综合分析这三者间的相关关系。针对前期克隆的玉米抗旱基因ZCN7,分析该基因区域转座子插入缺失变异介导的DNA甲基化和基因表达变化情况。【结果】 在田间干旱处理下,自交系H082183的叶片相对含水量和散粉-吐丝间隔均与正常处理没有显著差异,而旅28在所有试验材料中表现最低的叶片相对含水量和最大的散粉-吐丝间隔。利用H082183和旅28这两个抗旱性极端差异的玉米自交系开展基因组重测序和转座子插入分析,分别检测到333 754和333 296个转座子插入,其中,有89 954个转座子插入在2个自交系间具有多态性。基因组DNA甲基化分析表明,转座子、内含子和启动子区域较外显子和非编码区呈现较高的CG和CHG甲基化水平,经差异甲基化分析,在2个自交系间共检测到41 352个差异甲基化区域,其中60%的差异甲基化区域位于转座子插入缺失变异的上下游5 kb范围内。基因表达水平与基因的CG和CHG甲基化水平负相关,在2个自交系干旱下的叶片和根系中分别鉴定到4 196和3 500个差异表达基因,其中19.5%和19.7%与差异甲基化区域关联。通过对抗旱相关基因ZCN7的研究,发现该基因34 kb区间内的3个LTR类转座子插入,造成自交系旅28在干旱和正常处理下的CG和CHG甲基化显著高于抗旱自交系H082183,并且H082183中ZCN7表达量也显著高于旅28。【结论】 揭示了转座子介导的表观遗传调控在玉米响应干旱胁迫中的重要作用,进一步扩展了转座子变异和DNA甲基化调控抗旱基因ZCN7表达的分子机制。

关键词: 玉米(Zea mays L.), 抗旱性, 转座子, DNA甲基化, 转录组

Abstract:

【Objective】 Drought is one of the serious abiotic stresses influencing maize production worldwide. Understanding the molecular mechanisms underlying drought tolerance is of great importance in maize improvement. 【Method】 In the present study, representative maize inbred lines were selected for field drought experiment and the drought tolerance was estimated based on leaf relative water content and anthesis-silking interval. Two inbred lines with contrasting drought tolerance was used for genome resequencing and transposable element insertions were identified. The DNA methylation level of leaf and root tissues under different water treatments of the two lines were measured using Whole Genome Bisulfite Sequencing (WGBS). And the gene expression profiles of these samples were detected by RNA sequencing. The integrative atlas of transposable element insertion/deletion variants, differentially methylated regions and differentially expressed genes in the two lines were constructed. In addition, the transposable element insertion/deletion variant mediated epigenetic regulation of ZCN7, which has been conferred the drought tolerance function in our previous study, was analyzed. 【Result】 The field experiment showed inbred line H082183 showed highest drought tolerance, in which the leaf relative water content and anthesis-silking interval had no significant difference between drought and well-watered treatments. While the Lü28 displayed lowest leaf relative water content and largest anthesis-silking interval under drought. Thus, these two lines were selected for further analysis. A total of 333 754 and 333 296 transposable element insertions were identified in the genome of H082183 and Lü28, respectively. And 89 954 transposable element insertions were polymorphism between two lines. The transposable element insertions, introns and promoters showed higher CG and CHG methylation level than exons and untranslated regions. Furthermore, 41 352 differentially methylated regions were identified between H082183 and Lü28. And 60% of the differentially methylated regions were located in the transposable element insertion\deletion variants and 5 kb flanking regions. The gene expression level showed negatively correlated with CG and CHG methylation. Differentially expression analysis between H082183 and Lü28 obtained 4 196 and 3 500 differentially expressed genes in leaf and root under drought, respectively. The 19.5% and 19.7% of these differentially expressed genes were located in differentially methylated regions. Three LTR transposable element insertions were identified in the 34 kb region of ZCN7 in Lü28 but absent in the genome H082183. And the DNA methylation levels of CG and CHG in this genomic region were significantly higher in Lü28 than H082183 under both drought and well-watered environments, which conferred higher ZCN7 expression in the drought tolerant line H082183. 【Conclusion】 Our results highlight the important role of interplay of transposable element insertions, DNA methylation and gene expression under drought. And gene expression regulation mechanism of ZCN7 relied on the transposable element insertion/deletion variants mediated DNA methylation was proposed.

Key words: maize (Zea mays L.), drought tolerance, transposable element, DNA methylation, transcriptome