A next-generation marker genotyping platform (AmpSeq) in heterozygous crops: a case study for marker-assisted selection in grapevine

Abstract

Marker-assisted selection (MAS) is often employed in crop breeding programs to accelerate and enhance cultivar development, via selection during the juvenile phase and parental selection prior to crossing. Next-generation sequencing and its derivative technologies have been used for genome-wide molecular marker discovery. To bridge the gap between marker development and MAS implementation, this study developed a novel practical strategy with a semi-automated pipeline that incorporates trait-associated single nucleotide polymorphism marker discovery, low-cost genotyping through amplicon sequencing (AmpSeq) and decision making. The results document the development of a MAS package derived from genotyping-by-sequencing using three traits (flower sex, disease resistance and acylated anthocyanins) in grapevine breeding. The vast majority of sequence reads (⩾99%) were from the targeted regions. Across 380 individuals and up to 31 amplicons sequenced in each lane of MiSeq data, most amplicons (83 to 87%) had <10% missing data, and read depth had a median of 220–244×. Several strengths of the AmpSeq platform that make this approach of broad interest in diverse crop species include accuracy, flexibility, speed, high-throughput, low-cost and easily automated analysis.

References

[1] Myles S. Improving fruit and wine: what does genomics have to offer? Trends Genet 2013; 29: 190-196.

[2] Ru S, Main D, Evans K, Peace C . Current applications, challenges, and perspectives of marker-assisted seedling selection in Rosaceae tree fruit breeding. Tree Genet Genomes 2015; 11: 1-12.

[3] van Nocker S, Gardiner SE . Breeding better cultivars, faster: applications of new technologies for the rapid deployment of superior horticultural tree crops. Hortic Res 2014; 1: 14022.

[4] Bus VM, Esmenjaud D, Buck E, Laurens F . Application of genetic markers in rosaceous crops. In: Folta K, Gardiner S (eds). Genetics and Genomics of Rosaceae vol. 6. Springer International Publishing: New York, NY, USA. 2009 pp 563-599.

[5] Karaagac E, Vargas A, de Andrés M, Carreño I, Ibáñez J, Carreño J et al. Marker-assisted selection for seedlessness in table grape breeding. Tree Genet Genomes 2012; 8: 1003-1015.

[6] Eibach R, Zyprian E, Welter L, Topfer R . The use of molecular markers for pyramiding resistance genes in grapevine breeding. Vitis 2007; 46: 120-124.

[7] Zini E, Raffeiner M, Di Gaspero G, Eibach R, Grando MS, Letschka T . Applying a defined set of molecular markers to improve selections of resistant grapevine accessions. Acta Hortic 2015; 1082: 73-78.

[8] Collard BCY, Mackill DJ . Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Philos Trans R Soc Lond B Biol Sci 2008; 363: 557-572.

[9] Varshney RK, Thudi M, Nayak SN, Gaur PM, Kashiwagi J, Krishnamurthy L et al. Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.) . Theor Appl Genet 2014; 127: 445-462.

[10] Reisch B, Owens C, Cousins P . Grape. In: Badenes ML, Byrne DH (eds). Fruit Breeding vol. 8. Springer International Publishing: New York, NY, USA. 2012, 225-262.

[11] Shirasawa K, Hirakawa H . DNA marker applications to molecular genetics and genomics in tomato. Breeding Sci 2013; 63: 21-30.

[12] Boudehri K, Bendahmane A, Cardinet G, Troadec C, Moing A, Dirlewanger E . Phenotypic and fine genetic characterization of the D locus controlling fruit acidity in peach. BMC Plant Biol 2009; 9: 59.

[13] Alpert KB, Tanksley SD . High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato. Proc Natl Acad Sci USA 1996; 93: 15503-15507.

[14] De Franceschi P, Stegmeir T, Cabrera A, van der Knaap E, Rosyara UR, Sebolt AM et al. Cell number regulator genes in Prunus provide candidate genes for the control of fruit size in sweet and sour cherry. Mol Breed 2013; 32: 311-326.

[15] Spindel J, Wright M, Chen C, Cobb J, Gage J, Harrington S et al. Bridging the genotyping gap: using genotyping by sequencing (GBS) to add high-density SNP markers and new value to traditional bi-parental mapping and breeding populations. Theor Appl Genet 2013; 126: 2699-2716.

[16] Korte A, Farlow A . The advantages and limitations of trait analysis with GWAS: a review. Plant Methods 2013; 9: 29.

[17] Huang XH, Han B . Natural variations and genome-wide association studies in crop plants. Annu Rev Plant Biol 2014; 65: 531-551.

[18] Fodor A, Segura V, Denis M, Neuenschwander S, Fournier-Level A, Chatelet P et al. Genome-wide prediction methods in highly diverse and heterozygous species: Proof-of-concept through simulation in grapevine. PLoS One 2014; 9: e110436.

[19] Khan MA, Korban SS . Association mapping in forest trees and fruit crops. J Exp Bot 2012; 63: 4045-4060.

[20] Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK et al. Genetic structure and domestication history of the grape. Proc Natl Acad Sci USA 2011; 108: 3530-3535.

[21] Myles S, Chia JM, Hurwitz B, Simon C, Zhong GY, Buckler E et al. Rapid genomic characterization of the genus Vitis . PLoS One 2010; 5: e8219.

[22] Cattonaro F, Testolin R, Scalabrin S, Morgante M, Gaspero GD . Genetic diversity in the grapevine germplasm. In: Tuberosa R, Graner A, Frison E (eds). Genomics of Plant Genetic Resources. Springer International Publishing: New York, NY, USA. 2014 pp 683-704.

[23] Arnold C, Rossetto M, McNally J, Henry RJ . The application of SSRs characterized for grape (Vitis vinifera) to conservation studies in Vitaceae . Am J Bot 2002; 89: 22-28.

[24] Dirlewanger E, Denoyes-Rothan B, Yamamoto T, Chagné D . Genomics tools across Rosaceae species. In: Folta K, Gardiner S (eds). Genetics and Genomics of Rosaceae vol. 6. Springer International Publishing: New York, NY, USA. 2009 pp 539-561.

[25] Mnejja M, Garcia-Mas J, Audergon JM, Arus P . Prunus microsatellite marker transferability across rosaceous crops. Tree Genet Genomes 2010; 6: 689-700.

[26] Deschamps S, Llaca V, May GD . Genotyping-by-sequencing in plants. Biology 2012; 1: 460-483.

[27] Fechter I, Hausmann L, Daum M, Sorensen TR, Viehover P, Weisshaar B et al. Candidate genes within a 143 kb region of the flower sex locus in Vitis . Mol Genet Genomics 2012; 287: 247-259.

[28] Picq S, Santoni S, Lacombe T, Latreille M, Weber A, Ardisson M et al. A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines . BMC Plant Biol 2014; 14: 229.

[29] Grattapaglia D, Silva-Junior OB, Kirst M, de Lima BM, Faria DA, Pappas GJ Jr . High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species. BMC Plant Biol 2011; 11: 65.

[30] Yang HA, Li CD, Lam HM, Clements J, Yan GJ, Zhao SC . Sequencing consolidates molecular markers with plant breeding practice. Theor Appl Genet 2015; 128: 779-795.

[31] Moragues M, Comadran J, Waugh R, Milne I, Flavell AJ, Russell JR . Effects of ascertainment bias and marker number on estimations of barley diversity from high-throughput SNP genotype data. Theor Appl Genet 2010; 120: 1525-1534.

[32] Thomson MJ . High-throughput SNP genotyping to accelerate crop improvement. Plant Breed Biotech 2014; 2: 195-212.

[33] Miller AJ, Matasci N, Schwaninger H, Aradhya MK, Prins B, Zhong GY et al. Vitis phylogenomics: hybridization intensities from a SNP array outperform genotype calls. PLoS One 2013; 8: e78680.

[34] Myles S, Mahanil S, Harriman J, Gardner KM, Franklin JL, Reisch BI et al. Genetic mapping in grapevine using SNP microarray intensity values. Mol Breed 2015; 35: 88-100.

[35] Luby JJ, Shaw DV . Does marker-assisted selection make dollars and sense in a fruit breeding program? Hort Sci 2001; 36: 872-879.

[36] Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES et al. A robust, simple Genotyping-by-Sequencing (GBS) approach for high diversity species. PLoS One 2011; 6: e19379.

[37] He JF, Zhao XQ, Laroche A, Lu ZX, Liu HK, Li ZQ . Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding. Front Plant Sci 2014; 5: 484.

[38] Beissinger TM, Hirsch CN, Sekhon RS, Foerster JM, Johnson JM, Muttoni G et al. Marker density and read depth for genotyping populations using Genotyping-by-Sequencing. Genetics 2013; 193: 1073-1081.

[39] Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q et al. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS One 2014; 9: e90346.

[40] Swarts K, Li HH, Navarro JAR, An D, Romay MC, Hearne S et al. Novel methods to optimize genotypic imputation for low-coverage, next-generation sequence data in crop plants. Plant Genome 2014; 7(3): 1-12.

[41] Bianco L, Cestaro A, Sargent DJ, Banchi E, Derdak S, Di Guardo M et al. Development and validation of a 20 K Single Nucleotide Polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh) . PLoS One 2014; 9: e110377.

[42] Hyma KE, Barba P, Wang M, Londo JP, Acharya CB, Mitchell SE et al. Heterozygous mapping strategy (HetMappS) for high resolution Genotyping-By-Sequencing markers: a case study in grapevine. PLoS One 2015; 10: e0134880.

[43] Dalbò MA, Ye GN, Weeden NF, Wilcox WF, Reisch BI . Marker-assisted selection for powdery mildew resistance in grapes. J Am Soc Hortic Sci 2001; 126: 83-89.

[44] Reisch BI, Mahanil S, Consolie N, Luce RS, Wallace PG, Cadle-Davidson L . Examination of marker-assisted selection for powdery mildew resistance. Acta Hortic 2014; 1046: 151-155.

[45] Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES . TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 2007; 23: 2633-2635.

[46] Koressaar T, Remm M . Enhancements and modifications of primer design program Primer3. Bioinformatics 2007; 23: 1289-1291.

[47] Cadle-Davidson L, Chicoine DR, Consolie NH . Variation within and among Vitis spp. for foliar resistance to the powdery mildew pathogen Erysiphe necator . Plant Dis 2011; 95: 202-211.

[48] Feechan A, Kocsis M, Riaz S, Zhang W, Gadoury DM, Walker MA et al. Strategies for RUN1 deployment using RUN2 and REN2 to manage grapevine powdery mildew informed by studies of race specificity. Phytopathology 2015; 105: 1104-1113.

[49] Manns DC, Mansfield AK . A core-shell column approach to a comprehensive high-performance liquid chromatography phenolic analysis of Vitis vinifera L. and interspecific hybrid grape juices, wines, and other matrices following either solid phase extraction or direct injection . J Chromatogr A 2012; 1251: 111-121.

[50] Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 2007; 449: 463-467.

[51] Quinlan AR, Hall IM . BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 2010; 26: 841-842.

[52] Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA et al. The variant call format and VCFtools. Bioinformatics 2011; 27: 2156-2158.

[53] R core team : A language and environment for statistical computing ver v3.1.2. R Foundation for Statistical Computing: Vienna, Austria. 2014.

[54] Hanley JA . Receiver operating characteristic (ROC) methodology: the state of the art. Crit Rev Diagn Imaging 1989; 29: 307-335.

[55] Grzybowski M, Younger JG . Statistical methodology: III. Receiver operating characteristic (ROC) curves. Acad Emerg Med 1997; 4: 818-826.

[56] Pepe MS . Receiver operating characteristic methodology. J Am Stat Assoc 2000; 95: 308-311.

[57] Lawal B. Probability and probability distributions. In: Applied Statistical Methods in Agriculture, Health and Life Sciences. Springer International Publishing: New York, NY, USA. 2014 pp 59-113.

[58] Utro F, Haiminen N, Livingstone D, Cornejo OE, Royaert S, Schnell RJ et al. iXora: exact haplotype inferencing and trait association. BMC Genet 2013; 14: 48-63.

[59] Delaneau O, Zagury JF, Marchini J . Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods 2013; 10: 5-6.

[60] Knapp SJ . Marker-assisted selection as a strategy for increasing the probability of selecting superior genotypes. Crop Sci 1998; 38: 1164-1174.

[61] Bernardo R. Molecular markers and selection for complex traits in plants: Learning from the last 20 years. Crop Sci 2008; 48: 1649-1664.

[62] Guimarães EP . Food and Agriculture Organization of the United Nations. In: Marker-Assisted Selection: Current Status and Future Perspectives in Crops, Livestock, Forestry and Fish. Food and Agriculture Organization of the United Nations: Rome, Italy. 2007.

[63] Hospital F, Moreau L, Lacoudre F, Charcosset A, Gallais A . More on the efficiency of marker-assisted selection. Theor Appl Genet 1997; 95: 1181-1189.

[64] Iftekharuddaula KM, Salam MA, Newaz MA, Ahmed HU, Collard BCY, Septiningsih EM et al. Comparison of phenotypic versus marker-assisted background selection for the SUB1 QTL during backcrossing in rice. Breeding Sci 2012; 62: 216-222.

[65] Xu Y. Molecular Plant Breeding. CABI Publishing: Wallingford, CT, USA. 2010.

[66] Castro P, Rubio J, Madrid E, Fernandez-Romero MD, Millan T, Gil J . Efficiency of marker-assisted selection for ascochyta blight in chickpea. J Agr Sci 2015; 153: 56-67.

[67] Hackett CA, Broadfoot LB . Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps. Heredity 2003; 90: 33-38.

[68] Zhang L, Wang S, Li H, Deng Q, Zheng A, Li S et al. Effects of missing marker and segregation distortion on QTL mapping in F₂ populations . Theor Appl Genet 2010; 121: 1071-1082.

[69] Remington DL, Thornsberry JM, Matsuoka Y, Wilson LM, Whitt SR, Doeblay J et al. Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc Natl Acad Sci USA 2001; 98: 11479-11484.

[70] Di Gaspero G, Cipriani G, Adam-Blondon AF, Testolin R . Linkage maps of grapevine displaying the chromosomal locations of 420 microsatellite markers and 82 markers for R-gene candidates. Theor Appl Genet 2007; 114: 1249-1263.

[71] Di Gaspero G, Cipriani G, Marrazzo MT, Andreetta D, Castro MJP, Peterlunger E et al. Isolation of (AC)n-microsatellites in Vitis vinifera L. and analysis of genetic background in grapevines under Marker-assisted selection . Mol Breed 2005; 15: 11-20.

[72] Merdinoglu D, Butterlin G, Bevilacqua L, Chiquet V, Adam-Blondon AF, Decroocq S . Development and characterization of a large set of microsatellite markers in grapevine (Vitis vinifera L.) suitable for multiplex PCR . Mol Breed 2005; 15: 349-366.

[73] Calus MPL, Meuwissen THE, Windig JJ, Knol EF, Schrooten C, Vereijken ALJ et al. Effects of the number of markers per haplotype and clustering of haplotypes on the accuracy of QTL mapping and prediction of genomic breeding values. Genet Sel Evol 2009; 41: 11.

[74] Hayes BJ, Chamberlain AJ, McPartlan H, Macleod I, Sethuraman L, Goddard ME . Accuracy of marker-assisted selection with single markers and marker haplotypes in cattle. Genet Res 2007; 89: 215-220.

[75] Cervantes-Martinez C, Brown JS . A haplotype-based method for QTL mapping of F₁ populations in outbred plant species . Crop Sci 2004; 44: 1572-1583.

[76] Feng XH, Keim D, Wanjugi H, Coulibaly I, Fu Y, Schwarz J et al. Development of molecular markers for genetic male sterility in Gossypium hirsutum . Mol Breed 2015; 35: 141.

[77] Fernandez L, Le Cunff L, Tello J, Lacombe T, Boursiquot JM, Fournier-Level A et al. Haplotype diversity of VvTFL1A gene and association with cluster traits in grapevine (V. vinifera) . BMC Plant Biol 2014; 14: 209.

[78] Rafalski A. Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol 2002; 5: 94-100.

[79] Jones ES, Sullivan H, Bhattramakki D, Smith JSC . A comparison of simple sequence repeat and single nucleotide polymorphism marker technologies for the genotypic analysis of maize (Zea mays L.) . Theor Appl Genet 2007; 115: 361-371.

[80] Lambert SG, Asenstorfer RE, Williamson NM, Iland PG, Jones GP . Copigmentation between malvidin-3-glucoside and some wine constituents and its importance to colour expression in red wine. Food Chem 2011; 125: 106-115.

[81] Chagne D, Crowhurst RN, Troggio M, Davey MW, Gilmore B, Lawley C et al. Genome-wide SNP detection, validation, and development of an 8 K SNP array for apple. PLoS One 2012; 7: e31745.

[82] Vezzulli S, Micheletti D, Riaz S, Pindo M, Viola R, This P et al. A SNP transferability survey within the genus Vitis . BMC Plant Biol 2008; 8: 128.

[83] Gasperini D, Greenland A, Hedden P, Dreos R, Harwood W, Griffiths S . Genetic and physiological analysis of Rht8 in bread wheat: an alternative source of semi-dwarfism with a reduced sensitivity to brassinosteroids. J Exp Bot 2012; 63: 4419-4436.

[84] Mahanil S, Ramming D, Cadle-Davidson M, Owens C, Garris A, Myles S et al. Development of marker sets useful in the early selection of Ren4 powdery mildew resistance and seedlessness for table and raisin grape breeding. Theor Appl Genet 2012; 124: 23-33.

[85] Fennell A. Genomics and functional genomics of winter low temperature tolerance in temperate fruit crops. Crit Rev Plant Sci 2014; 33: 125-140.

[86] Fennell A, Mathiason K, Luby J . Genetic segregation for indicators of photoperiod control of dormancy induction in Vitis species . Acta Hortic 2005; 689: 533-539.

[87] Poland JA, Rife TW . Genotyping-by-Sequencing for plant breeding and genetics. Plant Genome 2012; 5: 92-102.

[88] Mathiason K, He D, Grimplet J, Venkateswari J, Galbraith DW, Or E et al. Transcript profiling in Vitis riparia during chilling requirement fulfillment reveals coordination of gene expression patterns with optimized bud break . Funct Integr Genomics 2009; 9: 81-96.

[89] Sreekantan L, Mathiason K, Grimplet J, Schlauch K, Dickerson JA, Fennell AY . Differential floral development and gene expression in grapevines during long and short photoperiods suggests a role for floral genes in dormancy transitioning. Plant Mol Biol 2010; 73: 191-205.

[90] Garris A, Clark L, Owens C, McKay S, Luby J, Mathiason K et al. Mapping of photoperiod-induced growth cessation in the wild grape Vitis riparia . J Am Soc Hortic Sci 2009; 134: 261-272.