Crops ›› 2023, Vol. 39 ›› Issue (6): 69-78.doi: 10.16035/j.issn.1001-7283.2023.06.010

Previous Articles     Next Articles

Preliminary Identification of Candidate Genes with Resistance to Powdery Mildew in Oil Flax Based on BSA-Seq

Ye Chunlei1(), Wang Wei1(), Chen Jun1, Chen Chen1, Luo Junjie1, Wang Yi2, Zhang Jianping3()   

  1. 1Bio-Technology Institute, Gansu Academy of Agricultural Sciences, Lanzhou 730070, Gansu, China
    2College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
    3Institute of Crops, Gansu Academy of Agricultural Sciences, Lanzhou 730070, Gansu, China
  • Received:2022-05-20 Revised:2022-09-08 Online:2023-12-15 Published:2023-12-15

RichHTML

4

PDF (PC)

91

Abstract:

Powdery mildew is one of the main diseases of oil flax. However, due to the lack of disease resistant germplasm resources and the lack of relevant gene mining and utilization, breeding of oil flax with powdery mildew resistance has been seriously restricted. In view of this, oil flax germplasm H98-1 with high resistance to powdery mildew and another oil flax germplasm D367 with high susceptibility to powdery mildew obtained by EMS mutation were used as parents, and used mixed pool grouping analysis (BSA) to perform whole genome resequencing on the parents and the resistant and susceptible F2 mixed pool. Candidate regions for disease resistance genes were identified through association analysis. BLAST software was applied to deeply annotate the coding genes within the candidate interval and identify candidate genes for resistance to powdery mildew. The results showed that the candidate regions related to the resistance to powdery mildew trait of H98-1 were mainly located in the intervals of 12 436 438-12 438 246, 12 475 667-12 476 462, 12 478 211-12 479 109 on chromosome 2, and 14 120 916-18 541 859 on chromosome 4 of oil flax, with a total length of 4.43Mb and containing 844 genes which are mainly involved in secondary metabolite biosynthesis, transport and catabolism, protein transport, signal transduction and other biological processes. The 11 genes may be related to the resistance of oil flax to powdery mildew which screened out by gene annotation, metabolic pathway analysis and literature analysis.

Key words: Oil flax, Powdery mildew, BSA, Candidate genes

Fig.1

Veen diagram of numbers of SNP and InDel among the tested samples (a) The Number of SNPs; (b) The Number of InDels"

Table 1

Annotation of variation sites of the tested samples"

变异类型
Variation type		 变异位点信息
Information of variation site		 H98-1 D367 抗病池
Disease-resistance pool		 感病池
Disease-susceptible pool
SNP 编码序列 63 580 51 467 65 069 61 267
基因上游 213 290 165 255 218 843 207 330
基因下游 165 650 130 082 168 902 159 594
基因区间 160 364 129 773 163 657 165 461
内含子 87 418 70 699 90 075 85 010
5'非翻译区 1528 1123 1469 1387
5'非翻译区 3715 2692 3616 3388
剪切供体突变 108 91 110 104
剪切受体突变 111 101 130 125
剪切位点区域突变 2932 2396 2959 2819
同义编码突变 31 290 25 525 31 942 30 032
非同义编码突变 31 283 25 102 32 085 30 258
起始密码子获得 50 32 43 43
起始密码子丢失 50 40 53 50
终止密码子获得 543 436 562 526
终止密码子丢失 334 295 344 325
InDel 编码序列 3479 2961 3675 3520
基因上游 41 329 34 072 44 359 42 494
基因下游 32 005 26 545 34 201 32 639
基因区间 21 841 18 153 23 327 23 461
内含子 20 219 17 296 21 277 20 449
5'非翻译区 521 404 545 513
5'非翻译区 1085 842 1082 1054
剪切供体突变 78 74 89 86
剪切受体突变 74 50 69 66
剪切位点区域突变 691 624 733 713
起始密码子丢失 23 17 21 20
终止密码子获得 42 40 41 40
终止密码子丢失 39 34 38 41
密码子插入 176 154 181 187
密码子删除 356 317 391 367

Table 2

Candidate region of powdery mildew resistance genes of flax germplasm H98-1"

染色体
Chromosome		 关联区域起始位置
Starting position of associated region		 关联区域终止位置
Ending position of associated region		 大小
Size (Mb)		 基因数量
Gene number
chr2 12433995 12433995 0.00 0
chr2 12436438 12438246 0.00 1
chr2 12439649 12440010 0.00 0
chr2 12442385 12442531 0.00 0
chr2 12475667 12476462 0.00 1
chr2 12478211 12479109 0.00 1
chr2 13270712 13279515 0.01 0
chr2 13388502 13388502 0.00 0
chr4 14120916 18541859 4.42 841
合计Total － － 4.43 844

Table 3

Gene function annotation of SNPs and InDels in candidate region"

功能注释数据库
Annotated-databases		 候选区域基因数
Number of candidate region genes		 非同义突变基因数
Number of nonsynonymous mutant genes		 移码突变基因数
Number of frameshift mutation genes
COG 292 116 36
GO 673 261 67
KEGG 586 221 58
NR 804 303 76
NT 815 306 77
SwissProt 627 239 63
trEMBL 805 303 76
合计Total 815 306 77

Fig.2

COG annotation classification map of genes in candidate region"

Table 4

Genes in candidate region annotated by GO database"

基因功能分类
Gene function classfication		 基因功能
Gene function		 基因数
Number of genes
细胞成分
Cellular component		 细胞外区域 12
细胞 86
细胞膜 75
细胞连接 3
膜封闭腔 6
大分子复合体 15
细胞器 64
胞外区部分 4
细胞器部分 20
细胞膜部分 61
细胞部分 86
突触 1
共质体 1
分子功能
Molecular function		 核酸结合转
录因子活性		 7
催化活性 127
受体活性 2
结构分子活性 3
转运子活性 8
结合 142
电子载体活性 2
抗氧化物活性 2
酶调节子活性 3
分子转导活性 2
生物过程
Biological process		 复制 8
免疫系统过程 1
代谢过程 141
细胞过程 120
复制过程 10
信号 7
多细胞生物过程 13
发育过程 13
生长 2
单组织过程 83
应激反应 44
定位 19
多组织过程 9
生物调控 52
细胞成分组织
或生物发生		 19

Fig.3

Distribution of gene pathways in candidate regions"

Table 5

Candidate genes information in candidate region"

序号
Number		 候选基因ID
Candidate gene ID		 染色体位置
Chromosome
position		 起始位点
Starting
position		 终止位点
Ending
position		 基因长度
Gene
length (bp)		 基因功能
Gene function
1 L.us.o.g.scaffold187.45 Chr4 16526088 16541454 15366 编码抗病蛋白L6
2 L.us.o.g.scaffold187.47 Chr4 16552817 16557562 4745 为拟南芥HOS1的同源基因，编码E3泛素蛋白连接酶
3 L.us.o.g.scaffold187.55 Chr4 16589848 16591005 1157 为拟南芥At5g62660的同源基因，推测编码F-box蛋白
4 L.us.o.g.scaffold187.89 Chr4 16705005 16708737 3732 编码抗病蛋白L6
5 L.us.o.g.scaffold187.98 Chr4 16754356 16759010 4654 编码抗病蛋白L6
6 L.us.o.g.scaffold187.99 Chr4 16762357 16763315 958 编码抗病蛋白L6
7 L.us.o.g.scaffold187.100 Chr4 16763361 16766387 3026 编码抗病蛋白L6
8 L.us.o.g.scaffold187.121 Chr4 16871587 16875166 3579 编码抗病蛋白L6
9 L.us.o.g.scaffold214.52 Chr4 17370672 17372100 1428 为拟南芥At3g58860的同源基因，推测编码F-box蛋白
10 L.us.o.g.scaffold318.36 Chr4 17211936 17222836 10900 编码抗病蛋白L6
11 L.us.o.g.scaffold322.22 Chr4 16970278 16972939 2661 编码抗病蛋白L6
[1] 国家胡麻产业技术体系. 中国现代农业产业可持续发展战略研究——胡麻分册. 北京: 中国农业出版社, 2016.
[2] Bekhit E D A, Shavandi A, Jodjaja T, et al. Flaxseed: Composition, detoxification, utilization, and opportunities. Biocatalysis and Agricultural Biotechnology, 2018, 13:129-152.
doi: 10.1016/j.bcab.2017.11.017
[3] Cullis C A. Genetics and genomics of Linum. Cham: Springer Nature Switzerland AG, 2019, 23:215-225.
[4] 王炜, 叶春雷, 陈琛, 等. 亚麻白粉病研究进展. 中国油料作物学报, 2019, 41(3):478-484.
[5] Kaushal P K, Srivas S R. Inheritance of resistance to powdery mildew (Oidium lini Skoric) in linseed (Linum usitatissimum L.). Current Science, 1974, 43(11):353-354.
[6] Sran R S, Paul S, Kumar A, et al. Genetics of resistance to rust and powdery mildew in linseed (Linum usitatissimum L.). Indian Phytopathology, 2021, 74:633-637.
doi: 10.1007/s42360-021-00349-9
[7] Rashid K, Duguid S. Inheritance of resistance to powdery mildew in flax. Canadian Journal of Plant Pathology, 2005, 27(3):404- 409.
doi: 10.1080/07060660509507239
[8] 杨学, 赵云, 关凤芝, 等. 亚麻品系980 1-1对白粉病的抗性遗传分析. 植物病理学报, 2008, 38(6):656-658.
[9] 张倩. 亚麻抗白粉病基因的定位. 哈尔滨: 黑龙江大学, 2015.
[10] 罗俊杰, 叶春雷, 欧巧明, 等. 抗白粉病胡麻种质资源田间鉴定与筛选. 植物保护, 2019, 45(5):259-262.
[11] Michelmore R W, Paran I, Kesseli R V. Identification of markers linked to disease-resistance genes by bulked segregant analysis: A rapid method to detect markers in specific genomic regions by using segregating populations. Proceedings of the National Academy of Sciences of the United States of America, 1991, 88 (2):9828-9832.
[12] Marks R A, Hotaling S, Frandsen P B, et al. Representation and participation across 20 years of plant genome sequencing. Nature Plants, 2021, 7(12):1571-1578.
doi: 10.1038/s41477-021-01031-8 pmid: 34845350
[13] Zou C, Wang P X, Xu Y B. Bulked sample analysis in genetics, genomics and crop improvement. Plant Biotechnology Journal, 2016, 14(10):1941-1955.
doi: 10.1111/pbi.12559 pmid: 26990124
[14] Zhang K J, Wang X, Zhu W W, et al. Complete resistance to powdery mildew and partial resistance to downy mildew in a Cucumis hystrix introgression line of cucumber were controlled by a co-localized locus. Theoretical and Applied Genetics, 2018, 131(10):2229-2243.
doi: 10.1007/s00122-018-3150-2
[15] Li C L, Ling F L, Su G H, et al. Location and mapping of the NCLB resistance genes in maize by bulked segregant analysis (BSA) using whole genome re-sequencing. Molecular Breeding, 2020, 40(10):92.
doi: 10.1007/s11032-020-01171-3
[16] 赵改会, 李书宇, 詹杰鹏, 等. 甘蓝型油菜角果数突变体基因的定位及候选基因分析. 作物学报, 2022, 48(1):27-39.
doi: 10.3724/SP.J.1006.2022.04281
[17] Zhang J P, Qi Y N, Wang L M, et al. Genomic comparison and population diversity analysis provide insights into the domestication and improvement of flax. iScience, 2020, 23(4):100967.
doi: 10.1016/j.isci.2020.100967
[18] McKenna A, Hanna M, Banks E, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research, 2010, 20(9):1297- 1303.
doi: 10.1101/gr.107524.110 pmid: 20644199
[19] Cingolani P, Platts A, Wang L L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain W1118; iso-2; iso-3. Fly, 2012, 6(2):80-92.
[20] Fekih R, Takagi H, Tamiru M, et al. MutMap+: genetic mapping and mutant identification without crossing in rice. PLoS ONE, 2013, 8(7):e68529.
doi: 10.1371/journal.pone.0068529
[21] Hill J T, Demarest B L, Bisgrove B W, et al. MMAPPR: mutation mapping analysis pipeline for pooled RNA-seq. Genome Research, 2013, 23(4):687-697.
doi: 10.1101/gr.146936.112 pmid: 23299975
[22] 张之昊, 王俊, 刘章雄, 等. 基于BSA-Seq技术挖掘大豆中黄622的多小叶基因. 作物学报, 2020, 46(12):1839-1849.
doi: 10.3724/SP.J.1006.2020.04075
[23] 甘露, 马含月, 高京草, 等. 瓜类蔬菜白粉病抗性诱导及抗性遗传研究进展. 中国瓜菜, 2021, 34(3):1-6.
[24] 周军, 徐如宏, 谢鑫, 等. 小麦抗白粉病及分子标记研究进展. 湖北农业科学, 2020, 59(6):10-15.
[25] 向贵生, 张真建, 王其刚, 等. 月季白粉病及其抗性研究进展. 江苏农业科学, 2017, 45(10):9-15.
[26] 张辉, 贾霄云, 高凤云, 等. 胡麻. 北京: 中国农业科学技术出版社, 2021.
[27] Li W, Essuman K, Anderson R G, et al. TIR domains of plant immune receptors are NAD+-cleaving enzymes that promote cell death. Science, 2019, 365(6455):799-803.
doi: 10.1126/science.aax1771
[28] Bernoux M, Ve T, Williams S, et al. Structural and functional analysis of a plant resistance protein TIR domain reveals interfaces for self-association, signaling, and autoregulation. Cell Host & Microbe, 2011, 9(3):200-211.
[29] Jiang B, Li M, Cheng Y, et al. Genetic mapping of powdery mildew resistance genes in soybean by high-throughput genome- wide sequencing. Theoretical and Applied Genetics, 2019, 132 (6):1833-1845.
doi: 10.1007/s00122-019-03319-y
[30] 胡玉慈, 姚立萍, 张童, 等. 黄毛草莓编码NB-ARC结构域的FnCN基因和启动子克隆及FnCN基因表达分析. 植物资源与环境学报, 2020, 29(4):1-10.
[31] Moffett P, Farnham G, Peart J, et al. Interaction between domains of a plant NBS-LRR protein in disease resistance-related cell death. The EMBO Journal, 2002, 21(17):4511-4519.
doi: 10.1093/emboj/cdf453
[32] 文志丰. 中国野生葡萄编码NB-ARC结构的抗白粉病基因克隆及功能分析. 杨凌: 西北农林科技大学, 2016.
[33] Hulbert S H, Webb C A, Smith S M, et al. Resistance gene complexes: evolution and utilization. Annual Review of Phytopathology, 2003, 39:285-312.
doi: 10.1146/phyto.2001.39.issue-1
[34] 金彦龙, 李艳军, 张新宇, 等. 利用SLAF-Seq结合BSA方法分子标记‘小白冬麦’抗白粉病基因mlxbd. 西北农业学报, 2019, 28(6):914-921.
[35] Zhu Y, Li Y, Fei F, et al. E3 ubiquitin ligase gene CMPG1-V from Haynaldia villosa L. contributes to powdery mildew resistance in common wheat (Triticum aestivum L.). The Plant Journal, 2015, 84(1):154-168.
doi: 10.1111/tpj.2015.84.issue-1
[36] Wang G, Yin H, Qiao X, et al. F-box genes: Genome-wide expansion, evolution and their contribution to pollen growth in pear (Pyrus bretschneideri). Plant Science, 2016, 253:164-175.
doi: 10.1016/j.plantsci.2016.09.009
[37] Jiang H Y, Wang C C, Ping L, et al. Pattern of LRR nucleotide variation in plant resistance genes. Plant Science, 2007, 173(2):253-261.
doi: 10.1016/j.plantsci.2007.05.010
[38] Paquis S, Mazeyrat-Gourbeyre F, Fernandez O, et al. Characterization of a F-box gene up-regulated by phytohormones and upon biotic and abiotic stresses in grapevine. Molecular Biology Reports, 2011, 38(5):3327-3337.
doi: 10.1007/s11033-010-0438-y pmid: 21104020
[1] Hou Jingjing, Jin Fang, Zhao Li, Wang Bin. Comprehensive Evaluation of Agronomic and Quality Traits of 16 New Oil Flax Lines [J]. Crops, 2022, 38(5): 42-48.
[2] Xiang Chao, Sun Suli, Zhu Zhendong, Zong Xuxiao, Yang Tao, Liu Rong, Yang Mei, Xian Dongfeng, Yang Xiuyan. Resistance and Molecular Identification to Powdery Mildew of Pea Germplasms in Sichuan [J]. Crops, 2021, 37(3): 51-56.
[3] Feng Xuejin,Guo Xiujuan,Yang Jianchun,Wang Liqin. Effects of Spraying Selenium Fertilizer on Selenium Content, Yield and Quality of Flax Seed [J]. Crops, 2019, 35(3): 155-157.
[4] Shuyan Wang,Bing Han,Simin Zhou,Jun Xu. Correlation Analysis between the Expression of Stearoyl Acyl Carrier Protein Dehydrogenase Gene and the Soluble Sugar and Fat Content of Oil Flax [J]. Crops, 2016, 32(4): 56-61.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!