基于BSA-seq技术挖掘芝麻株高相关候选基因

doi:10.7668/hbnxb.20192363

摘要/Abstract

摘要： 为了进一步挖掘芝麻株高相关基因，为适机收芝麻新品种选育提供理论指导，以冀航芝1号和DW607为亲本构建F₂群体，并构建以株高为目标性状的极端混池，利用BSA-seq技术，采用ED和Δ（SNP-index或InDel-index）2种方法挖掘株高相关的染色体区段，注释区段内的基因信息，利用GO和KEGG等数据库分析注释基因的功能。结果发现，亲本之间共获得298 634个SNP和76 360个InDel，混池之间共获得24 048个SNP和9 630个InDel；基于SNP标记的ED方法关联到5个染色体区段，ΔSNP-index方法关联到3个染色体区段，两者的交集有3个；基于InDel标记的ED方法关联到5个染色体区段，ΔInDel-index方法关联到8个染色体区段，两者的交集有8个；4个染色体区段同时被SNP和InDel标记关联到，共注释到330个基因，比对到的前20个KEGG通路主要包括植物激素信号转导和能量代谢；GO富集结果表明，有18个基因参与生长素响应，可能是参与株高调控的关键基因。

关键词: 芝麻, 株高, 混池测序, 候选基因

Abstract: In order to further explore sesame plant height related genes and provide theoretical guidance for the selection of machine harvesting sesame varieties, Jihangzhi 1 and DW607 were used as parents to construct an F₂ population. Bulked segregant pools were constructed with the plant height as the target trait using the F₂ population, and BSA-seq was carried out. Candidate chromosome regions were obtained using ED and ΔSNP/InDel-index methods. KEGG pathway and GO annotation were used for gene function analysis. A total of 298 634 SNPs and 76 360 InDels were found between the parents, and 24 048 SNPs and 9 360 InDels were found between the pools. Based on the SNP markers, the ED method obtained 5 chromosome regions, the ΔSNP-index method obtained 3 chromosome regions, and 3 regions were common between them. Based on the InDel markers, the ED method obtained 5 chromosome regions, the ΔInDels-index method obtained 8 chromosome regions, and 8 regions were common between them. Four regions were common between the SNP and the InDel merkers results, and 330 genes were annotated in these 4 regions. The first 20 KEGG pathways mainly involved plant hormone signal transduction and energy metabolism. By GO annotation analysis, 18 genes were involved in auxin response, which might be vital in plant height regulation.

Key words: Sesame, Plant height, BSA-seq, Candidate gene

中图分类号:

S565.03
Q78

崔彦芹, 郭元章, 侯少锋, 李思达, 关中波, 徐桂真. 基于BSA-seq技术挖掘芝麻株高相关候选基因[J]. 华北农学报, 2021, 36(S1): 8-15. doi: 10.7668/hbnxb.20192363.

CUI Yanqin, GUO Yuanzhang, HOU Shaofeng, LI Sida, GUAN Zhongbo, XU Guizhen. Candidate Gene Detection for Sesame Plant Height Based on BSA-seq[J]. ACTA AGRICULTURAE BOREALI-SINICA, 2021, 36(S1): 8-15. doi: 10.7668/hbnxb.20192363.

http://www.hbnxb.net/CN/Y2021/V36/IS1/8

参考文献

[1] Wei W L,Zhang Y X,Lü H,Li D H,Wang L H,Zhang X R. Association analysis for quality traits in a diverse panel of Chinese sesame(Sesamum indicum L.)germplasm[J]. Journal of Integrative Plant Biology,2013,55(8):745-758.doi:10.1111/jipb.12049.
[2] Johnson L A,Suleiman T M,Lusas E W. Sesame protein:A review and prospectus[J]. Journal of the American Oil Chemists' Society,1979,56(3):463-468.doi:10.1007/BF02671542.
[3] 丁霞. 芝麻株高相关性状遗传分析和QTL定位[D].北京:中国农业科学院,2013.
Ding X. Inheritance analysis and QTL mapping on plant height related traits in sesame(Sesamum indicum L.)[D].Beijing:Chinese Academy of Agricultural Sciences,2013.
[4] 王瑞元. 中国为全球芝麻产业的发展作出了重要贡献[J].中国油脂,2019,44(12):1-2.doi:10.3969/j.issn.1003-7969.2019.12.001.
Wang R Y. China has made important contributions to the development of global sesame industry[J]. China Oils and Fats,2019,44(12):1-2.
[5] 丁霞,王林海,张艳欣,黎冬华,高媛,危文亮,王蕾,张秀荣. 芝麻核心种质株高构成相关性状的遗传变异及关联定位[J].中国油料作物学报,2013,35(3):262-270.doi:10.7505/j.issn.1007-9084.2013.03.006.
Ding X,Wang L H,Zhang Y X,Li D H,Gao Y,Wei W L,Wang L,Zhang X R. Genetic variation and associated mapping for traits related to plant height constitutions in core collections of sesame(Sesamum indicum L.)[J]. Chinese Journal of Oil Crop Sciences,2013,35(3):262-270.
[6] 钟明志,魏淑红,彭正松,杨在君. 小麦Rht矮秆基因研究和应用综述[J].分子植物育种,2018,16(20):6670-6677.doi:10.13271/j.mpb.016.006670.
Zhong M Z,Wei S H,Peng Z S,Yang Z J. A review of the research and application of Rht dwarf genes in wheat[J]. Molecular Plant Breeding,2018,16(20):6670-6677.
[7] Liu F,Wang P D,Zhang X B,Li X F,Yan X H,Fu D H,Wu G. The genetic and molecular basis of crop height based on a rice model[J]. Planta,2018,247(1):1-26.doi:10.1007/s00425-017-2798-1.
[8] Sakamoto T,Matsuoka M. Generating high-yielding varieties by genetic manipulation of plant architecture[J]. Current Opinion in Biotechnology,2004,15(2):144-147.doi:10.1016/j.copbio.2004.02.003.
[9] 朱晓凤. 芝麻遗传图谱构建和株高相关性状QTL分析[D].北京:中国农业科学院,2015.doi:10.7666/d.Y2787588.
Zhu X F. Construction of genetic linkage map and QTL analysis for plant height associated traits in sesame(Sesame indicum L.)[D].Beijing:Chinese Academy of Agricultural Sciences,2015.
[10] Wei X,Liu K Y,Zhang Y X,Feng Q,Wang L H,Zhao Y,Li D H,Zhao Q,Zhu X D,Zhu X F,Li W J,Fan D L,Gao Y,Lu Y Q,Zhang X M,Tang X M,Zhou C C,Zhu C R,Liu L F,Zhong R C,Tian Q L,Wen Z R,Weng Q J,Han B,Huang X H,Zhang X R. Genetic discovery for oil production and quality in sesame[J]. Nature Communications,2016,6(1):1-10.doi:10.1038/ncomms9609.
[11] Zhang H Y,Miao H M,Li C,Wei L B,Duan Y H,Ma Q,Kong J J,Xu F F,Chang S X. Ultra-dense SNP genetic map construction and identification of SiDt gene controlling the determinate growth habit in Sesamum indicum L.[J]. Scientific Reports,2016,6:31556.doi:10.1038/srep31556.
[12] Wang L H,Xia Q J,Zhang Y X,Zhu X D,Zhu X F,Li D H,Ni X M,Gao Y,Xiang H T,Wei X,Yu J Y,Quan Z W,Zhang X R. Updated sesame genome assembly and fine mapping of plant height and seed coat color QTLs using a new high-density genetic map[J]. BMC Genomics,2016,17:31.doi:10.1186/s12864-015-2316-4.
[13] Miao H M,Li C,Duan Y H,Wei L B,Ju M,Zhang H Y. Identification of a Sidwf1 gene controlling short internode length trait in the sesame dwarf mutant dw607[J]. Theoretical and Applied Genetics,2020,133(1):73-86.doi:10.1007/s00122-019-03441-x.
[14] Wang L H,Yu S,Tong C B,Zhao Y Z,Liu Y,Song C,Zhang Y X,Zhang X D,Wang Y,Hua W,Li D H,Li D,LiF,Yu J Y,Xu C Y,Han X L,Huang S M,Tai S H,Wang J Y,Xu X,Li Y R,Liu S Y,Varshney R K,Wang J,Zhang X R. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis[J]. Genome Biology,2014,15(2):39.doi:10.1186/gb-2014-15-2-r39.
[15] Li H,Durbin R. Fast and accurate long read alignment with Burrows-Wheeler transform[J]. Bioinformatics,2010,26(5):1754-1760.doi:10.1093/bioinformatics/btp324.
[16] McKenna A,Hanna M,Banks E,Sivachenko A,Cibulskis K,Kernytsky A,Garimella K,Altshuler D,Gabriel S,Daly M,DePristo M A. The Genome Analysis Toolkit:A MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome Research,2010,20(9):1297-1303.doi:10.1101/gr.107524.110.
[17] Cingolani P,Platts A,Wang L L,Coon M,Nguyen T,Wang L,Land S J,Lu X Y,Ruden D M. A program for annotating and predicting the effects of single nucleotide polymorphisms,SnpEff[J]. Fly,2012,6(2):80-92.doi:10.4161/fly.19695.
[18] Hill J T,Demarest B L,Bisgrove B W,Gorsi B,Su Y C,Yost H J. MMAPPR:mutation mapping analysis pipeline for pooled RNA-seq[J]. Genome Research,2013,23(4):687-697.doi:10.1101/gr.146936.112.
[19] 张之昊,王俊,刘章雄,邱丽娟. 基于BSA-Seq技术挖掘大豆中黄622的多小叶基因[J].作物学报,2020,46(12):1839-1849.doi:10.3724/SP.J.1006.2020.04075.
Zhang Z H,Wang J,Liu Z X,Qiu L J. Mapping of an incomplete dominant gene controlling multifoliolate leaf by BSA-seq in soybean(Glycine max L.)[J]. Acta Agronomica Sinica,2020,46(12):1839-1849.
[20] Altschul S F,Madden T L,Schäffer A A,Zhang J H,Zhang Z,Miller W,Lipman D J. Gapped BLAST and PSI-BLAST:A new generation of protein database search programs[J]. Nucleic Acids Research,1997,25(17):3389-3402.doi:10.1093/nar/25.17.3389.
[21] Ashburner M,Ball C A,Blake J A,Botstein D,Butler H,Cherry J M,Davis A P,Dolinski K,Dwight S S,Eppig J T,Harris M A,Hill D P,Issel-Tarver L,Kasarskis A,Lewis S,Matese J C,Richardson J E,Ringwald M,Rubin G M,Sherlock G. Gene ontology:Tool for the unification of biology[J]. Nature Genetics,2000,25(1):25-29.doi:10.1038/75556.
[22] Kanehisa M,Goto S,Kawashima S,Okuno Y,Hattori M. The KEGG resource for deciphering the genome[J]. Nucleic Acids Research,2004,32(Database issue):D277-D280.doi:10.1093/nar/gkh063.
[23] Anilakumar K.R,Pal A,Khanum F,Bawa A S. Nutritional,medicinal and industrial uses of sesame(Sesamum indicum L.)seeds-an overview[J]. Agriculturae Conspectus Scientificus,2010.75(4):159-168.doi:10.1083/jcb.77.3.685.
[24] Hedden P. The genes of the green revolution[J]. Trends in Genetics,2003,19(1):5-9.doi:10.1016/S0168-9525(02)00009-4.
[25] Win K T,Vegas J,Zhang C Y,Song K,Lee S. QTL mapping for downy mildew resistance in cucumber via bulked segregant analysis using next-generation sequencing and conventional methods[J]. Theoretical and Applied Genetics,2017,130(1):199-211.doi:10.1007/s00122-016-2806-z.
[26] Chen W,Yao J B,Chu L,Yuan Z W,Li Y,Zhang Y S. Genetic mapping of the nulliplex-branch gene(gb_nb1)in cotton using next-generation sequencing[J]. Theoretical and Applied Genetics,2015,128(3):539-547.doi:10.1007/s00122-014-2452-2.
[27] Klein H,Xiao Y G,Conklin P A,Govindarajulu R,Kelly J A,Scanlon M J,Whipple C J,Bartlett M. Bulked-segregant analysis coupled to whole genome sequencing(BSA-seq)for rapid gene cloning in maize[J]. G3,2018,8(11):3583-3592.doi:10.1534/g3.118.200499.
[28] Xu D D,Sun D,Diao Y L,Liu M X,Gao J,Wu B,Yuan X M,Lu P,Zhang Z W,Zhang J,Guo G G. Fast mapping of a chlorophyll b synthesis-deficiency gene in barley(Hordeum vulgare L.)via bulked-segregant analysis with reduced-representation sequencing[J]. The Crop Journal,2019,7(1):58-64.doi:10.1016/j.cj.2018.07.002.
[29] Liu H Y,Zhou F,Zhou T,Yang Y X,Zhao Y Z. Fine mapping of a novel male-sterile mutant showing wrinkled-leaf in sesame by BSA-Seq technology[J]. Industrial Crops and Products,2020,156:112862.doi:10.1016/j.indcrop.2020.112862.
[30] 朱晓凤,黎冬华,王林海,张艳欣,高媛,魏鑫,张秀荣.矮秆与高秆芝麻株高建成中内源激素含量变化比较分析[J].中国油料作物学报,2015,37(1):83-89.doi:10.7505/j.issn.1007-9084.2015.01.013.
Zhu X F,Li D H,Wang L H,Zhang Y X,Gao Y,Wei X,Zhang X R. Phytohormone change during plant height development between dwarfs and high genotypes of sesame[J]. Chinese Journal of Oil Crop Sciences,2015,37(1):83-89.

选择文件类型/文献管理软件名称

选择包含的内容