利用2b-RAD测序结合HRM分析技术开发与梨矮生性状相关的DNA分子标记

doi:10.3864/j.issn.0578-1752.2017.15.014

摘要/Abstract

摘要： 【目的】果树的矮生性状是一种重要的农艺性状，对果树的集约化栽培具有重要意义。来自西洋梨实生变异品种‘Le Nain Vert’的矮生性状受控于一个单显性基因PcDw，目前关于该基因的序列信息等还不清楚。本研究的目的是开发与其紧密连锁的DNA分子标记，为鉴定该基因提供依据。【方法】根据分离群体分组分析的原理，以‘矮生梨’×‘茌梨’和‘2-3’×‘绿宝石’2个F₁杂交分离群体为试材，应用IIB型限制性内切酶的RAD技术（Restriction association site DNA，2b-RAD）对2对矮生型/普通型对比基因池进行基因组测序分析，在对比基因池间筛选出单核苷酸多态性（Single nucleotide polymorphism，SNP）位点，从中筛查出位于PcDw，以确定其与PcDw位点的连锁关系。【结果】对4个样本（即4个对比基因池）的2b-RAD标签测序文库的测序结果共产生67 186 260条reads，平均每个样本测序reads数为16 796 565。将原始reads进行质量过滤后的统计结果表明，每个样本获得平均unique标签数目为86 810，平均测序深度为77定位染色体上的SNPs，用高分辨率熔解曲线分析技术（High-resolution melting analysis，HRM）在群体上进一步检测验证×，该测序深度能够达到准确分型的标准。SOAP软件定位结果表明，4个测序文库中含有酶切位点的高质量reads占测序原始reads的70%以上，表明测序质量较好。在来自2个不同群体的矮生型基因池与普通型基因池间进行对比分析，初步筛选出SNP位点1 317个，其中有8个位于PcDw的定位染色体scaffold00074上。用HRM技术对这8个SNP标记在群体上的进一步检测结果表明，在‘矮生梨’×‘茌梨’群体上有2个SNP标记、在‘2-3’×‘绿宝石’群体上有4个SNP标记表现出与PcDw位点共分离的特性，根据其扩增子的熔解曲线形状差异，可有效区分矮生型和普通型表型。在来自‘矮生梨’×‘茌梨’群体的215个杂种后代和来自‘2-3’×‘绿宝石’群体的168个杂种后代中，未发现有标记与性状的重组类型。【结论】2b-RAD测序技术与HRM分析技术相结合，进行果树重要农艺性状分子标记的开发，是一种行之有效的方法。基于这一策略，本研究鉴定获得了4个与西洋梨矮生性状单显性基因PcDw连锁的SNP标记。

关键词: 梨, 矮生性状, PcDw, SNP标记, 2b-RAD, HRM

Abstract: 【Objective】As an important agronomic trait, the dwarf character of tree is significant to fruit intensive culture. The dwarf traits of pear originated from variety ‘Le Nain Vert’, a chance seedling of Pyrus communis, is determined by a dominant gene PcDw, which sequence information is still unknown. Development of DNA molecular markers tightly linked to the PcDw locus could provide important information for identifying of this gene.【Method】Using 2 different F₁ populations obtained from ‘Aishengli’×‘Chili’ and ‘2-3’×‘Lvbaoshi’, respectively, as plant materials, 2 pairs of bulks for dwarf type/standard type were sequenced through IIB restriction association site DNA (2b-RAD) technology, and then single nucleotide polymorphism (SNP) markers were detected between the dwarf and standard bulks based on the principle of bulked segregant analysis. After that, those SNPs located on the PcDw mapped chromosome were selected and tested on the whole population to confirm their linkage relationship to the PcDw gene by high-resolution melting (HRM) analysis.【Result】A total of 67 186 260 reads were produced from the 2b-RAD sequencing of the 4 samples (2 pairs of bulks). That means the average reads per sample was 16 796 565. Statistical analysis of the filtered raw reads showed that the average unique tags per sample were 86 810, and the average sequencing depth was 77× which was enough for accurate genotyping. For the 4 libraries, mapping through SOAP software indicated that the high quality reads containing the restriction site accounted for more than 70% of the raw reads, which means the high quality of sequencing. Totally 1 317 polymorphic SNP markers were screened between the dwarf bulks and standard bulks derived from the 2 different populations, and 8 of them were mapped on the PcDw-localized chromosome scaffold00074. Further detection of the 8 SNPs through HRM analysis showed that 2 and 4 SNP markers co-segregated with the PcDw locus were identified from population ‘Aishengli’×‘Chili’ and population ‘2-3’×‘Lvbaoshi’, respectively. According to the different shapes of the melting curves of amplicons, the dwarf/normal phenotype could be distinguished effectively. The recombinants for each marker and the target trait were not found in both of the two populations, which contained 215 and 168 progenies, respectively. 【Conclusion】The combination of 2b-RAD sequencing and HRM analysis technology is an efficient method for exploiting molecular markers of important agronomic traits in fruit trees. In this study, a total of 4 SNP markers co-segregated with the PcDw locus were identified based on this strategy.

Key words: pear, dwarf traits; PcDw, SNP marker, 2b-RAD, HRM

肖玉雄，王彩虹，田义轲，杨绍兰，李鼎立，张海月. 利用2b-RAD测序结合HRM分析技术开发与梨矮生性状相关的DNA分子标记[J]. 中国农业科学, 2017, 50(15): 3006-3012.

XIAO YuXiong, WANG CaiHong, TIAN YiKe, YANG ShaoLan, LI DingLi, ZHANG HaiYue. Development of DNA Molecular Markers for the Dwarf Trait in Pear Through the Method of 2b-RAD Sequencing and HRM Analysis[J]. Scientia Agricultura Sinica, 2017, 50(15): 3006-3012.

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: https://www.chinaagrisci.com/CN/10.3864/j.issn.0578-1752.2017.15.014

https://www.chinaagrisci.com/CN/Y2017/V50/I15/3006

参考文献

[1] Chen B Y, Wang C H, Tian Y K, Chu Q, Hu C. Anatomical characteristics of young stems and mature leaves of dwarf pear. Scientia Horticulturae, 2015, 186: 172-179.

[2] Rivalta L, Dradi M, Rosati C. Thirty years of pear breeding activity at ISF Forli, Italy. Acta Horticulturae, 2002, 596(1): 233-238.

[3] Wang C H, Tian Y K, Buck E J, Gardiner S E, Dai H Y, Jia Y. Genetic mapping of PcDw determining pear dwarf trait. Journal of the American Society for Horticultural Science, 2011, 136(1): 48-53.

[4] 李炜, 田义轲, 王彩虹, 白牡丹, 侯董亮. 通过HRM技术筛查与梨矮生性状决定位点PcDw紧密连锁的SNP标记. 园艺学报, 2015, 42(2): 214-220.

Li W, Tian Y K, Wang C H, Bai M D, Hou D L. Screening of SNP markers tightly linked to PcDw locus determining pear dwarf trait using HRM technology. Acta Horticulturae Sinica, 2015, 42(2): 214-220. (in Chinese)

[5] Davey J W, Blaxter M L. RADSeq: Next-generation population genetics. Brief Funct Genomics, 2011, 9(5): 416-423.

[6] Miller M R, Dunham J P, Amores A, Cresko W A, Johnson E A. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Research, 2007, 17(2): 240-248.

[7] Baird N A, Etter P D, Atwood T S, Currey M C, Shiver A L, Lewis Z A, Selker E U, Cresko W A, Johnson E A. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE, 2008, 3(10): e3376.

[8] Peterson B K, Weber J N, Kay E H, Fisher H S, Hoekstra H E. Double digest RADseq: An inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE, 2012, 7(5): e37135.

[9] Elshire R J, Glaubitz J C, Sun Q, Poland J A, Kawamoto K, Buckler E S, Mitchell S E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE, 2011, 6(5): e19379.

[10] Wang S, Meyer E, Mckay J K, Matz M V. 2b-RAD: A simple and flexible method for genome-wide genotyping. Nature Methods, 2012, 9(8): 808-810.

[11] Dou J, Zhao X, Fu X, Jiao W, Wang N, Zhang L, Hu X, Wang S, Bao Z. Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions. Biology Direct, 2012, 7(1): 1-9.

[12] Chutimanitsakun Y, Nipper R W, Cuesta-Marcos A, Cistué L, Corey A, Filichkina T, Johnson E A, Hayes P M. Construction and application for QTL analysis of a restriction site associated DNA (RAD) linkage map in barley. BMC Genomics, 2011, 12(1): 4.

[13] Pfender W F, Saha M C, Johnson E A, Slabaugh M B. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theoretical and Applied Genetics, 2011, 122(8): 1467-1480.

[14] Chagné D, Crowhurst R N, Pindo M, Thrimawithana A, Deng C, Ireland H, Fiers M, Dzierzon H, Cestaro A, Fontana P, Bianco L, Lu A, Storey R, Knäbell M, Saeedl M, Montanari1 S, Kim Y K, Nicolini D, Larger S, Stefani E, Allan A C, Bowen J, Harvey I, Johnston J, Malnoy M, Troggio M, Perchepied L, Sawyer G, Wiedowl C, Won K, Viola R, Hellens R P, Brewer L, Bus V G , Schaffer R J, Gardiner S E, Velasco R. The draft genome sequence of european pear (Pyrus communis L. ‘Bartlett’). PLoS ONE, 2014, 9: e92644.

[15] 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 USA, 1991, 88: 9828-9832.

[16] Li R, Yu C, Li Y, Lam T W, Yiu S M, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics, 2009, 25 (15): 1966-1967.

[17] Baumgartner I O, Kellerhals M, Costa F, Dondini L, Pagliarani G, Gregori R, Tartarini S, Leumann L, Laurens F, Patocchi A. Development of SNP-based assays for disease resistance and fruit quality traits in apple (Malus × domestica Borkh.) and validation in breeding pilot studies. Tree Genetics & Genomes, 2016, 12(3): 1-21.

[18] Zeballos J L, Abidi W, Giménez R, Monforte A J, Moreno M Á, Gogorcena Y. Mapping QTLs associated with fruit quality traits in peach [Prunus persica (L.) Batsch] using SNP maps. Tree Genetics & Genomes, 2016, 12(3): 1-17.

[19] Wu J, Li L T, Li M, Khan M A, Li X G, Chen H, Yin H, Zhang S L. High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers. Journal of Experimental Botany, 2014, 65(20): 5771-5781.

[20] SUN R, CHANG Y, YANG F, WANG Y, LI H, ZHAO Y, CHEN D, WU T, ZHANG X, HAN Z. A dense SNP genetic map constructed using restriction site-associated DNA sequencing enables detection of QTLs controlling apple fruit quality. BMC Genomics, 2015, 16(1): 747.

[21] 李全林, 王义发, 韩晴, 袁政, 沈雪芳. 糯玉米‘沪五彩花糯1号’品系特异性2b-RAD分子标记的开发及应用. 植物生理学报, 2016, 52(5): 669-677.

Li Q L, Wang Y F, Han Q, Yuan Z, Shen X F. Development and application of molecular markers for event-speci?c identi?cation of waxy corn ‘Huwucaihuanuo 1’ based on 2b-RAD technique. Plant Physiology Journal, 2016, 52(5): 669-677. (in Chinese)

[22] Jiao W, Fu X, Dou J, Li H, Su H, Mao J, Yu Q, Zhang L, Hu X, Huang X, Wang Y, Wang S, Bao Z. High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc. DNA Research, 2013, 21(1): 85-101.

[23] Seetharam A S, Stuart G W. Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments. PeerJ, 2013, 1: e226.

[24] Fletcher R S, Mullen J L, Yoder S, Bauerle W L, Reuning G, Sen S, Meyer E, Juenger T E, McKay J K. Development of a next-generation NIL library in Arabidopsis thaliana for dissecting complex traits. BMC Genomics, 2013, 14(1): 655-655.

[25] Guo Y, Yuan H, Fang D, Song L, Liu Y, Liu Y, Wu L, Yu J, Li Z, Xu X. An improved 2b-RAD approach (2b-RAD) offering genotyping tested by a rice (Oryza sativa L.) F2 population. BMC Genomics, 2014, 15 (1): 956.

[26] Li C, Li Y, Bradbury P J, Wu X, Shi Y, Song Y, Zhang D, Rodgers-Melnick E, Buckler E S, Zhang Z. Construction of high-quality recombination maps with low-coverage genomic sequencing for joint linkage analysis in maize. BMC Biology, 2015, 13: 78.

[27] Pecoraro C, Babbucci M, Villamor A, Franch R, Papetti C, Leroy B, Ortega-Garcia S, Muir J, Rooker J, Arocha F, Murua H, Zudaire I, Chassot E, Bodin N, Tinti F, Bargelloni L, Cariani A. Methodological assessment of 2b-RAD genotyping technique for population structure inferences in yellowfin tuna (Thunnus albacares). Marine Genomics, 2015, 25: 43-48.

[28] Pauletto M, Carraro L, Babbucci M, Lucchini R, Bargelloni L, Cardazzo B. Extending RAD tag analysis to microbial ecology: A comparison between multilocus sequence typing and 2b-RAD to investigate listeriamonocytogenes genetic structure. Molecular Ecology Resources, 2016, 16(3): 823-835.

[29] Lu Z, Niu L, Chagné D, Cui G, Pan L, Foster T, Zhang R P, Zeng W F, Wang Z Q. Fine mapping of the temperature-sensitive semi-dwarf (Tssd) locus regulating the internode length in peach (Prunus persica). Molecular Breeding, 2016, 36(2): 1-11.

[30] Wang C H, Li W, Tian Y K, Hou D L, Bai M D. Development of molecular markers for genetic and physical mapping of the PcDw locus in pear (Pyrus communis L.). Journal of Horticultural Science & Biotechnology, 2016, 91(3): 299-307.

[31] 赵俊生, 杨晓燕, 曾祥有, 钟声, 方静, 罗剑斌, 曾运友, 向旭. 利用SNP 分子标记分析化橘红种质资源. 分子植物育种, 2016, 14(5): 1203-1211.

Zhao J S, Yang X Y, Zeng X Y, Zhong S, Fang J, Luo J B, Zeng Y Y, Xiang X. Analysis on germplasm resources of Exocarpium citri grandis using SNP molecular markers. Molecular Plant Breeding, 2016, 14(5): 1203-1211.