作物杂志,2023, 第4期: 31–37 doi: 10.16035/j.issn.1001-7283.2023.04.005

• 遗传育种·种质资源·生物技术 • 上一篇    下一篇

基于转录组测序揭示玉米抗倒伏相关基因和代谢通路

刘松涛1(), 田再民1, 刘子刚1, 高志佳2, 张静2, 贺东刚3, 黄智鸿1, 兰鑫1   

  1. 1河北北方学院/河北省农产品食品质量安全分析检测重点实验室,075000,河北张家口
    2河北巡天农业科技有限公司,075000,河北张家口
    3河北兆育种业集团有限公司,050000,河北石家庄
  • 收稿日期:2022-01-02 修回日期:2022-04-02 出版日期:2023-08-15 发布日期:2023-08-15
  • 通讯作者: 智鸿,主要从事作物高产栽培研究,E-mail:hbnuhzh@163.com
  • 作者简介:刘松涛,主要从事作物遗传育种研究,E-mail:1129600629@qq.com
  • 基金资助:
    河北省现代农业产业技术体系(玉米体系岗位专家)(HBCT2018020203);河北省科技支撑重点研发项目(18226334D);张家口市科学技术局项目(1911012C)

Transcriptomic Analysis to Reveal Lodging Resistance Genes and Metabolism Pathways in Maize (Zea mays L.)

Liu Songtao1(), Tian Zaimin1, Liu Zigang1, Gao Zhijia2, Zhang Jing2, He Donggang3, Huang Zhihong1, Lan Xin1   

  1. 1Hebei North University/Key Laboratory of Hebei Province Agricultural Products Food Quality and Safety Analysis and Testing, Zhangjiakou 075000, Hebei, China
    2Hebei Universe Agricultural Science and Technology Co., Ltd., Zhangjiakou 075000, Hebei, China
    3Hebei Zhaoyu Seed Industry Group Co., Ltd., Shijiazhuang 050000, Hebei, China
  • Received:2022-01-02 Revised:2022-04-02 Online:2023-08-15 Published:2023-08-15

摘要:

倒伏是影响玉米产量的重要因素之一。以抗倒性不同的3个玉米品种[京农科728:高抗倒性(H),金农738:中抗倒性(M),先玉335:低抗倒性(L)]为材料进行转录组学分析,挖掘玉米抗倒性相关基因。结果表明,3个比较组共鉴定到10 093个差异表达基因(differentially expressed genes,DEGs),其中先玉335与京农科728比较组鉴定到的DEGs最多,为7779个。GO功能富集分析表明,3个玉米品种抗倒性不同可能与富集到相同GO条目的DEGs数量不同有关。代谢通路富集分析表明,L-vs-H、M-vs-H分组显著富集到苯丙素生物合成、次生代谢产物的生物合成和类黄酮生物合成途径,而L-vs-M分组显著富集到光合作用-天线蛋白、植物-病原体相互作用途径。茎秆显微结构表明,先玉335单个维管束面积最小,茎秆表皮细胞厚度最薄,京农科728单个维管束面积最大,表皮细胞厚度最厚。结果进一步明确了玉米抗倒伏相关的基因与代谢通路,为定向克隆和抗倒伏新品种的分子设计育种奠定了基础。

关键词: 玉米, 转录组测序, 显微结构, 差异表达基因, 抗倒伏

Abstract:

Lodging is one of the important factors affecting maize production. In this study, three maize varieties with different lodging resistance [Jingnongke 728 with high lodging resistance (H), Jinnong 738 with medium lodging resistance (M) and Xianyu 335 with low lodging resistance (L)] were used for transcriptome analysis to explore the genes related to maize lodging resistance. The results showed that a total of 10093 differentially expressed genes (DEGs) were identified in the three comparison groups, and the most DEGs of 7779 were identified in Xianyu 335 and Jingnongke 728 groups. GO functional enrichment analysis showed that the different lodging resistance of the three maize varieties may be related to the different number of DEGs enriched into the same GO terms. Metabolic pathway enrichment analysis showed that L-vs-H, M-vs-H were significantly enriched in phenylpropanoid biosynthesis, secondary metabolite biosynthesis and flavonoid biosynthesis, while L-vs-M was significantly enriched in photosynthesis antenna protein and plant-pathogen interaction pathways. The stem microstructure showed that Xianyu 335 had the smallest single vascular bundle area and the thinnest stem epidermal cell thickness, while Jingnongke 728 had the largest single vascular bundle area and the thickest epidermal cell thickness. The results of this study further clarified the genes and metabolic pathways related to lodging resistance in maize, which laid a foundation for directional cloning, and also provided a foundation for molecular design breeding of new lodging resistant varieties.

Key words: Maize, RNA-seq, Microscopic structure, Differentially expressed genes, Lodging-resistant

表1

荧光定量引物序列

基因名称Gene name 前引物Forward primer 后引物Reverse primer 产物长度Product size (bp)
LOC100281532 GGAGGAGATGATGGGCAGC CTCGATCTTCACCAGGGGC 79
LOC100273579 AGAAGTCGCTGAGCCTGAAC TCTGCATCAGCGGGTAGTTG 141
LOC100272756 AGGAGGACAAGTCCGTGGAG TATCGATCTTGTCGAGGCCG 174
LOC541914 CCGTTAACCTGTCGAGGCTT GTCTCCAACCTTCCAGCTCC 123
LOC100272970 ATCAGCGTACATCGCGTCC AACTGAGGTGGGCTCTGTGT 94
LOC100194371 CGGTTCTGCTTCAAGACGAT CGGTAGTAGCTCCTTGGGTG 113
LOC100281042 TGGTGGAGGAGTACAGGAGG GTGCCGTTCATCATGCTGTC 172
LOC100282047 GTCGAACAGCGAGGAGTACC GTCGAACCAGAGGTGGAACC 163
LOC100283318 CTACTTCATCTCGGAGGGGC GAGTAGGTGTGGAAGTCGGC 177
LOC100284765 CCACCATCAGTAGCGGTCG GCCCGTGATGTTGCTGGA 168

表2

转录组测序数据统计

样品名称
Sample name
原始数据
Raw read
有效数据
Clean read
碱基所占百分比
Q30(%)
GC含量
GC content(%)
单一比对率
Unique mapped
多比对率
Multi mapped
JNK728-1 46743554 44778802 94.64 53.21 38131199(85.15%) 2833048(3.58%)
JNK728-2 42213792 40278042 93.90 53.24 34173234(84.84%) 2543399(3.56%)
JNK728-3 45140086 42175780 94.60 53.16 35820137(84.93%) 2740012(3.64%)
JN738-1 45628206 43858376 94.36 53.69 38079538(86.82%) 2519952(3.37%)
JN738-2 49113472 46927518 94.46 53.78 40671536(86.67%) 2782255(3.41%)
JN738-3 47537856 45732776 93.90 53.72 39628836(86.65%) 2690579(3.43%)
XY335-1 45920308 44134008 94.30 54.57 38518679(87.28%) 2222552(3.04%)
XY335-2 44670476 42529790 94.55 54.84 37079717(87.19%) 2186084(3.04%)
XY335-3 46723572 44596496 94.02 55.08 38796470(86.99%) 2319389(3.06%)

图1

转录组测序样本PCA分析

表3

差异表达基因统计

分组
Group
DEGs数量
DEGs number
下调
Down-regulated
上调
Up-regulated
L-vs-H 7779 4411 3368
L-vs-M 5373 2896 2477
M-vs-H 4905 2831 2074

图2

分组鉴定到的DEGs Venn图分析

图3

DEGs的GO富集分析 颜色标尺(红―白)表示富集程度的高低(P<0.05)

图4

DEGs代谢通路富集分析

图5

3个玉米品种茎秆的显微结构 (a)显微镜可视范围内茎秆维管束数目;(b)单个维管束面积;(c)茎秆表皮细胞厚度

图6

3个玉米品种茎秆的显微结构 不同字母表示差异显著(P<0.05)

图7

10个DEGs的RNA-seq测序结果的qRT-PCR验证

[1] 王聪, 马青美, 郭新梅, 等. 玉米茎秆抗倒伏相关性状的遗传效应分析. 玉米科学, 2021, 29(4):9-17.
[2] 丰光, 黄长玲, 邢锦丰. 玉米抗倒伏的研究进展. 作物杂志, 2008(4):12-14.
[3] 王夏青, 宋伟, 张如养, 等. 玉米茎秆抗倒伏遗传的研究进展. 中国农业科学, 2021, 54(11):2261-2272.
doi: 10.3864/j.issn.0578-1752.2021.11.002
[4] 马天军. 玉米倒伏种类及防治对策研究. 新农业, 2021(21):14.
[5] 李妍妍, 景希强, 丰光, 等. 玉米倒伏的主要相关因素研究进展. 辽宁农业科学, 2013(4):47-51.
[6] 徐天成, 周富亮, 黄联润, 等. 玉米抗倒伏性相关研究进展. 安徽农学通报, 2019, 25(11):29-32.
[7] Cloninger FR, Zuber MS. Methods for evaluating stalk quality in corn. Phytopathology, 1970, 60:295-300.
doi: 10.1094/Phyto-60-295
[8] 冯素伟, 姜小苓, 胡铁柱, 等. 不同小麦品种茎秆显微结构与抗倒强度关系研究. 中国农学通报, 2012, 28(36):57-62.
[9] LivakKJ, SchmittgenTD. Analysis of relative gene expression data using real-time quantitative PCR andthe 2-∆∆CT method. Methods, 2001, 25:402-408.
doi: 10.1006/meth.2001.1262 pmid: 11846609
[10] 殷静. 大豆茎倒伏突变体基因定位、转录组测序及2个木质素合成基因的克隆研究.南京:南京农业大学, 2013.
[11] 杨硕, 郑云霄, 黄亚群, 等. 不同玉米自交系茎秆上部维管束数目的差异分析. 植物遗传资源学报, 2021, 22(2):466-475.
[12] 张巫军, 段秀建, 姚雄, 等. 烯效唑对遮阴下重穗型水稻渝香203茎秆形态结构和抗倒伏性的影响. 四川农业大学学报, 2019, 37(6):755-761.
[13] 胡丹. 甜荞抗倒伏相关性状的遗传分析及木质素合成特征. 重庆:西南大学, 2016.
[14] 刘丹. 木质素相关基因表达调控对油菜木质素合成及抗菌核病和抗倒伏性的影响. 北京: 中国农业科学院, 2008.
[1] 袁刘正, 王会强, 王秋岭, 朱世蝶, 赵月强, 袁曼曼, 王会涛, 张运栋, 柳家友, 袁永强. 遮阴条件下玉米自交系配合力与遗传效应分析[J]. 作物杂志, 2023, (4): 104–109
[2] 郑飞, 陈静, 崔亚坤, 孔令杰, 孟庆长, 李杰, 刘瑞响, 张美景, 赵文明, 陈艳萍. 淮北不同生态区丰产稳产宜机收玉米新品种筛选[J]. 作物杂志, 2023, (4): 110–117
[3] 王丽萍, 白岚方, 王天昊, 王宵璇, 白云鹤, 王玉芬. 不同施氮水平对青贮玉米植株氮素积累和转运的影响[J]. 作物杂志, 2023, (4): 165–173
[4] 李雨新, 卢敏, 赵久然, 王荣焕, 徐田军, 吕天放, 蔡万涛, 张勇, 薛洪贺, 刘月娥. 京津唐夏玉米区生产现状调研分析[J]. 作物杂志, 2023, (4): 174–181
[5] 谢文锦, 李方明, 李宁, 杨海龙, 付俊, 张中伟, 高旭东, 丰光. 基于GGE双标图的北方地区鲜食糯玉米产量和品质性状及试点鉴别力分析[J]. 作物杂志, 2023, (4): 85–90
[6] 温胜慧, 杨俊伟, 王洋, 李公建, 翁建峰, 段灿星, 贾鑫, 王建军. 玉米抗真菌病害基因挖掘与分子育种利用研究进展[J]. 作物杂志, 2023, (3): 1–11
[7] 常青, 李立军, 渠佳慧, 张艳丽, 韩冬雨, 赵鑫瑶. 土默川平原麦后复种饲用玉米‖油菜模式的增产优势及氮素利用效率[J]. 作物杂志, 2023, (3): 167–174
[8] 郭书磊, 王瑛, 魏良明, 张新, 刘焱, 吴伟华, 卢道文, 雷晓兵, 王振华, 鲁晓民. 不同生态条件下玉米产量影响因素分析[J]. 作物杂志, 2023, (3): 205–214
[9] 高沐甜, 邱冠杰, 朱通通, 李瑞莲, 邓敏, 罗红兵, 黄成. 玉米―大刍草渗入系群体的剑叶遗传基础解析[J]. 作物杂志, 2023, (3): 51–57
[10] 李忠南, 王越人, 马艺文, 相洋, 邬生辉, 曲海涛, 李福林, 张淑琴, 李光发. 基于玉米DH系分离群体叶鞘、花丝、花药与穗轴颜色性状的遗传分析[J]. 作物杂志, 2023, (3): 75–79
[11] 张盼盼, 李川, 张美微, 赵霞, 黄璐, 刘京宝, 乔江方. 氮肥减施下添加硝化抑制剂对夏玉米植株及土壤氮素积累分配及产量的影响[J]. 作物杂志, 2023, (2): 145–150
[12] 崔淑娜, 王晔, 卢雨晴, 潘金豹, 张秋芝. 玉米穗三叶与产量的相关和通径分析[J]. 作物杂志, 2023, (2): 201–206
[13] 孟亚轩, 姚旭航, 周宝元, 刘颖慧, 袁进成, 马玮, 赵明. 青贮玉米混合青贮研究进展[J]. 作物杂志, 2023, (2): 24–29
[14] 于永涛, 张楠, 谢利华, 李光玉, 刘建华, 李武, 李高科, 胡建广. 消费者在甜玉米种质食味品质鉴评中的偏好性初探[J]. 作物杂志, 2023, (1): 14–19
[15] 张东霞, 秦安振. 冬小麦-夏玉米作物蒸散量及其水热关系研究[J]. 作物杂志, 2022, (6): 145–151
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!