基于BSA-Seq技术的甜瓜抗霜霉病InDel标记开发

doi:10.6048/j.issn.1001-4330.2021.12.014

摘要/Abstract

摘要： 目的获得与甜瓜抗霜霉病基因连锁的分子标记或功能标记,用于甜瓜分子标记辅助选择育种。为进一步克隆甜瓜抗霜霉病主效基因奠定基础。方法以野生高抗霜霉病资源DM-4和感病自交系DM-2为亲本,构建F₂代分离群体,利用BSA-seq对甜瓜抗霜霉基因进行QTL定位,开发InDel分子标记,并在双亲及其F₂群体单株进行筛选验证。结果甜瓜抗霜霉病主效QTL位于第9号染色体。在候选区间开发了37个InDel分子标记,有9个PCR产物大小在抗、感双亲表现出差异,用F₂单株验证,结合田间表型鉴定结果,引物InDel 15和InDel 20准确率分别为95 %和98 %,引物InDel 15和InDel 20在抗病亲本中扩增产物大小分别为251和349 bp,在感病亲本中分别为231和324 bp。结论引物InDel 15和InDel 20可以作为分子标记进行抗霜霉病辅助选育,加快了甜瓜抗霜霉病育种速度。

关键词: 甜瓜; 霜霉病; InDel分子标记; BSA

Abstract:

【Objective】 o obtain molecular markers of melon downy mildew resistance genes, and improve the accuracy and scientific city of melon downy mildew resistance traits selection. 【Methods】 The F₂ generation segregation population was constructed by hybridization between wild high-resistant downy mildew resource DM-4 and susceptible inbred line DM-2BSA-seq was used to locate and regulate the candidate interval of downy mildew resistance gene in melon InDel molecular markers. parents and their F₂ individuals were used for screening and verification.【Results】 Genes regulating melon downy mildew resistance are located on 9th chromosome. 37 InDel molecular markers were developed in the candidate region, and 9 of them showed differences in the size of PCR products between resistant and susceptible parents. Further verification was carried out with F₂ plants. Combined with field phenotypic identification results, the accuracy rates of primer InDel 15 and InDel 20 were 95 % and 98 %, respectively. The amplified products of primer InDel 15 and InDel 20 in resistant parents were 251 bp and 349 bp, respectively, and 231 bp and 324 bp, respectively.【Conclusion】 Primers InDel 15 and InDel 20 can be used as molecular markers for downy mildew resistance breeding, which can accelerate the breeding speed of melon downy mildew resistance and lay the foundation for further cloning the main genes of melon downy mildew resistance.

Key words: Cucumis melo L.; downy mildew; In Del marker; BSA

中图分类号:

S652

凌悦铭, 李寐华, 杨永, 杨文莉, 范蓉, 伊鸿平, 张红, 张学军. 基于BSA-Seq技术的甜瓜抗霜霉病InDel标记开发[J]. 新疆农业科学, 2021, 58(12): 2265-2273.

LING Yueming, LI Meihua, YANG Yong, YANG Wenli, FAN Rong, YI Hongping, ZHANG Hong, ZHANG Xuejun. Development of InDel Marker for Melon Resistance to Downy Mildew Based on BSA Resequencing[J]. Xinjiang Agricultural Sciences, 2021, 58(12): 2265-2273.

图/表 8

参考文献 26

[1]	王志丹. 中国甜瓜产业经济发展研究[D]. 中国农业科学院, 2014.
	WANG Zhidan. Study on the Development of Muskmelon Industry Economy in China[D]. Chinese Beijing: Academy of Agricultural Sciences, 2014.
[2]	Wallace E C, D'Arcangelo K N, Quesada-Ocampo L M. Population Analyses Reveal Two Host-Adapted Clades of Pseudoperonospora cubensis, the Causal Agent of Cucurbit Downy Mildew, on Commercial and Wild Cucurbits[J]. Phytopathology, 2020,110(9):1-10.
[3]	Ulrich Gisi, Helge Sierotzki. Fungicide modes of action and resistance in downy mildews[J]. European Journal of Plant Pathology, 2008,122(1):157-167.
[4]	Yigal Cohen, Kyle M,. Van den Langenberg, Todd C. Wehner, et al. Resurgence of Pseudoperonospora cubensis: The Causal Agent of Cucurbit Downy Mildew[J]. Phytopathology, 2015,105(7):998-1012.
[5]	贺玉花, 徐永阳, 徐志红, 等. 甜瓜霜霉病抗性研究进展[J]. 果树学报, 2014,31(2):324-334.
	HE Yuhua, XU Yongyang, XU Zhihong, et al. Research progress of resistance to downy mildew in melon[J]. Journal of Fruit Science, 2014,31(2):324-334.
[6]	Takagi H, Tamiru M, Abe A, et al. MutMap accelerates breeding of a salt-tolerant rice cultivar[J]. Nat Biotechnol, 2015,33:445-449.
[7]	Takagi H, Abe A, Yoshide K, et al. QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations[J]. The Plant Journal, 2013,74:174-183.
[8]	Lu H F, Liu T, Jo?l K, et al. QTL-seq identifies an early flowering QTL located near flowering Locus T in cucumber[J]. Theor Appl Genet, 2014,127:1491-1499.
[9]	Illa-Berenguer E, Houten J V, Huang Z J, et al. Rapid and reliable identification of tomato fruit weight and locule number loci by QTL-seq[J]. Theor Appl Genet, 2015,128:1329-1342.
[10]	Zhong C, Sun S, Li Y, et al. Next-generation sequencing to identify candidate genes and develop diagnostic markers for a novel Phytophthora resistance gene, RpsHC18, in soybean[J]. Theor Appl Genet, 2018,131:525-538.
[11]	张学军, 郭丽霞, 马新力, 等. 甜瓜种质资源抗霜霉病和蔓枯病苗期筛选与评价[J]. 中国瓜菜, 2014,27(S1):38-40.
	ZHANG Xuejun, GUO Lixia, MA Xinli, et al. Screening and evaluation of melon germplasm resources against downy mildew and vine blight at seedling stage[J]. China Cucurbits and Vegetables, 2014,27(S1):38-40.
[12]	Criswell A D, Wehner T C, Klosinska U, et al. Use of sporulation and other leaf and vine traits for evaluation of resistance to downy mildew in cucumber. In: Pitrat M (ed)[C]. Proceedings of the IXth EUCARPIA meeting on genetics and breeding of Cucurbitaceae, INRA, Avignon, France, 2008, 433-440.
[13]	Épinat C,. and Pitrat M. Inheritance of resistance to downy mildew (Pseudoperonospora cubensis) in muskmelon (Cucumis melo). I. Analysis of a 8 x 8 diallel table[J]. Agronomie, 1994,14(4):239-248.
[14]	Épinat C,. and Pitrat M. Inheritance of resistance to downy mildew (Pseudoperonospora cubensis) in muskmelon (Cucumis melo). II. Generation means analysis of 5 genitors[J]. Agronomie, 1994,14(4):249-257.
[15]	Heng Li, and Richard Durbin. Fast and accurate short read alignment with Burrows-Wheeler transform[J]. Bioinformatics, 2009,25(14):1754-1760.
[16]	Lebeda A, Stepankova J, Krskova M, et al. Resistance in Cucumis melo germplasm to Pseudoperonospora cubensis pathotypes. In: Advances in Downy Mildew Research[C]. Proceedings of the 2nd International Downy Mildews Symposium, 2nd International Downy Mildews Symposium, Olomouc, Czech Republic, 2007: 157-167.
[17]	Ivanoff S S, Resistance of cantaloupes to downy mildew and the melon aphid[J]. Journal of Heredity, 1944,35(2):35-39.
[18]	李靖, 利容千, 袁文静. 黄瓜感染霜霉病菌叶片中一些酶活性的变化[J]. 植物病理学报, 1991,(4):39-45.
	LI Jing, LI Rongqian, YUAN Wenjing. ON the Change of Enzyme Activites of Cucumber Leaf Infected By Pseudoperonospora Cubensis(Berk. Et Ctrt) Rosws[J]. Acta Phytopathologica Sinica, 1991,(4):39-45.
[19]	杨柳燕, 徐永阳, 徐志红, 等. 甜瓜霜霉病抗性遗传及SRAP分子标记[J]. 江苏农业学报, 2012,28(5):1200-1202.
	YANG Liuyan, XU Yongyang, XU Zhihong, et al. Inheritance of downy mildew resistance in melon and SRAP marker linked to resistant genes[J]. Jiangsu Journal of Agricultural Sciences, 2012,28(5):1200-1202.
[20]	贺玉花, 徐永阳, 徐志红, 等. 甜瓜霜霉病抗性基因的SSR标记[J]. 江苏农业科学, 2014,42(7):54-55.
	HE Yuhua, XU Yongyang, XU Zhihong, et al. SSR markers of downy mildew resistance genes in melon[J]. Jiangsu Agricultural Sciences, 2014,42(7):54-55.
[21]	张学军, 宁雪飞, 杨永, 等. 甜瓜PI390452霜霉病抗性基因的QTL定位[J]. 新疆农业科学, 2016,53(12):2157-2165.
	ZHANG Xuejun, NING Xuefei, YANG Yong, et al. QTL Mapping of Resistance Genes to Downy Mildew in Cucumis melo ssp. melo PI390452[J]. Xinjiang Agricultural Sciences, 2016,53(12):2157-2165.
[22]	张学军, 杨俊涛, 李寐华, 等. 甜瓜PI414723抗霜霉病基因SSR分子标记[J]. 华北农学报, 2016,31(S1):86-91.
	ZHANG Xuejun, YANG Juntao, LI Meihua, et al. SSR Markers Linked of the Resistance Gene for Downy Mildew in Cucumis melo ssp. melo PI414723[J]. Acta Agriculturae Boreali-Sinica, 2016,31(S1):86-91.
[23]	徐婷婷, 汪巧玲, 邹淑琼, 等. 基于高通量测序的大麦InDel标记开发及应用[J]. 作物学报, 2020,46(9):1340-1355.
	XU Tingting, WANG Qiaoling, ZOU Shuqiong, et al. Development and application of InDel markers based on high throughput sequencing in barley[J]. Acta Agronomica Sinica, 2020,46(9):1340-1355.
[24]	Islam M R, Hossain M R, Jesse M I, et al. Development of Molecular Marker Linked with Bacterial Fruit Blotch Resistance in Melon (Cucumis melo L.)[J].Genes, 2020, 11(2):220-.
[25]	Hassan M Z, Robin A, Rahim M A, et al. Gummy Stem Blight Resistance in Melon: Inheritance Pattern and Development of Molecular Markers[J]. International Journal of Molecular Sciences, 2018. 19(10):2914.
[26]	Sathishkumar N, Hoy-Taek K, Kumar T S, et al. Whole Genome Re-Sequencing and Characterization of Powdery Mildew Disease-Associated Allelic Variation in Melon[J]. Plos One, 2016,11(6):e0157524.

样本名 Sample	原始数据量 Raw Base (bp)	有效数据量 Clean Base (bp)	碱基错误率 Error Rate(%)	Q20(%)	Q30(%)	GC含量 GC Content(%)
DM-4	5 230 484 100	5 205 287 100	0.01	97.74	94.94	37.39
DM-2	5 954 570 100	5 926 054 500	0.01	97.83	95.10	37.23
DM-R	9 302 519 100	9 266 289 000	0.01	97.93	95.27	37.17
DM-S	9 931 454 700	9 878 084 400	0.01	97.81	95.08	37.12

样本名 Sample	原始数据量 Raw Base (bp)	有效数据量 Clean Base (bp)	碱基错误率 Error Rate(%)	Q20(%)	Q30(%)	GC含量 GC Content(%)
DM-4	5 230 484 100	5 205 287 100	0.01	97.74	94.94	37.39
DM-2	5 954 570 100	5 926 054 500	0.01	97.83	95.10	37.23
DM-R	9 302 519 100	9 266 289 000	0.01	97.93	95.27	37.17
DM-S	9 931 454 700	9 878 084 400	0.01	97.81	95.08	37.12

样本名 Sample	比对序列数 Mapped reads	有效序列数 Total reads	比对率 Mapping rate(%)	平均测序深度 Average depth(X)	1X覆盖度 Coverage at least 1X(%)
DM-2	36 714 906	39 507 030	92.93	13.11	96.38
DM-4	32 024 597	34 701 914	92.28	11.57	96.37
DM-S	61 510 228	65 853 896	93.40	21.29	97.75
DM-R	57 254 273	61 775 260	92.68	19.85	97.70

样本名 Sample	比对序列数 Mapped reads	有效序列数 Total reads	比对率 Mapping rate(%)	平均测序深度 Average depth(X)	1X覆盖度 Coverage at least 1X(%)
DM-2	36 714 906	39 507 030	92.93	13.11	96.38
DM-4	32 024 597	34 701 914	92.28	11.57	96.37
DM-S	61 510 228	65 853 896	93.40	21.29	97.75
DM-R	57 254 273	61 775 260	92.68	19.85	97.70

名称 Name	正向引物序列5’- 3’ Forward primers sequence 5’ - 3’	反向引物5’- 3’ Reverse primers sequence 5’ - 3’	是否具有多态性 Polymorphism or not
InDel 1	AGGTAAAGCCCCAGAGAGGA	CCTCCCGTGAAGTTTGGAGT	否
InDel 2	AGTGAATGTGGTCATTGTGCT	CCTCCCGTGAAGTTTGGAGT	是
InDel 3	ACACAGCACAACCCAGACAA	CTACGACCAACCCCCACATG	是
InDel 4	GCACTCTCCAGCCTTGAACT	GAGCAAAGGTACCGCTGTCT	是
InDel 5	TTTTCATTCGGATTCCTTGTTGT	AAGCTATGAGACCTACGCGC	否
InDel 6	TCAAAGTGACCCCATCCAGTG	TGTTTTGAGTTCAACGGGGGA	否
InDel 7	GTGACCCACGATAGGTACACC	TCAAGTTTGGGAGATTGGATCA	否
InDel 8	AGCATTCTCCCAGTTTCCCC	TGGGAACAAGAAAGTCAGGCA	否
InDel 9	ACACAATGAACCAAGATATTTGCCA	TGTTTGTAAGCCACTCGACT	是
InDel 10	ACGACACTAAATGACAGTCCGT	CCAAAATGGAAACCCATGCCC	否
InDel 11	GCAAAGTTGACCCTCTGTTCG	CAGCGTTCTAGGGGTTGGTT	否
InDel 12	TGTGAGAATGACAGCCCACA	CACATCTCTCAAGAGGCGGA	否
InDel 13	CCTCGCTCTATTTGCATTCACA	TTTATACGCGTACTAGTTCTCTCA	否
InDel 14	AAAGTCAGGATGGCCGAGTG	TCTCTTTCATCAGCGGCGTT	是
InDel 15	TCACAGCTCACAATCACAGGT	AGGGGAAGCTGGAGAAGACA	是
InDel 16	ACCTGCCGTCTATCGTGTTT	TGAGGTCGGTCCAAACAAGT	否
InDel 17	GCTGCCGTTTCGTTTTCAGT	ACGACCAAGAAACCCAGGAA	否
InDel 18	CCGTTGATTGCGCTAACACA	TGCCGCACGAGTATTCACTT	是
InDel 19	AGGAGTTTGGGTATGACACGAC	TTCTCCCTTTACCCGTTGCC	是
InDel 20	TGGGGTTGCTTGTGGATAGT	AGGCGAACACCTTTTGATTCA	是
InDel 21	GGAAGAGGAAGAGGCAGAGG	ACGAATCATTTCCTCCTCCGA	否
InDel 22	TGCTCTTGAGGCTAGGGCTA	TAGAGAGGTCCGAGCCCATC	否
InDel 23	TCATTGGCTAGAACTTCGACCA	CGAAAGTGATGGCCATGGGA	否
InDel 24	AGTGGTTGGTGGTGCATAGG	GCAACACTTTGATCCCCACA	否
InDel 25	ACCTCCAAGAAATGAATTCAGGA	TGAACGACGACGACCAATCT	否
InDel 26	CTTGCAAGGGCTAATGGTGC	ACGAATCAAAGGTGTAGCCA	否
InDel 27	ACAGTTTAGCTTCCATCCCTACA	TCGTTTAGGTCATTGTCCTCCA	否