Exploring the diversity and genetic structure of the U.S. National Cultivated Strawberry Collection

Abstract

The cultivated strawberry (Fragaria × ananassa) arose through a hybridization of two wild American octoploid strawberry species in a French garden in the 1750s. Since then, breeders have developed improved cultivars adapted to different growing regions. Diverse germplasm is crucial to meet the challenges strawberry breeders will continue to address. The USDA-ARS National Clonal Germplasm Repository (NCGR) in Corvallis, Oregon maintains the U.S. strawberry collection. Recent developments in high-throughput genotyping for strawberry can provide new insights about the diversity and structure of the collection, germplasm management, and future breeding strategies. Genotyping was conducted on 539 F. × ananassa accessions using either the iStraw35 or FanaSNP 50 K Axiom array. Data for markers shared by the two arrays were curated for call quality, missing data, and minor allele frequency resulting in 4033 markers for structure assessment, diversity analysis, pedigree confirmation, core collection development, and the identification of haplotypes associated with desirable traits. The F. × ananassa collection was equally diverse across the different geographic regions represented. K-means clustering, sNMF, and UPGMA hierarchal clustering revealed seven to nine sub-populations associated with different geographic breeding centers. Two 100 accession core collections were created. Pedigree linkages within the collection were confirmed. Finally, accessions containing disease resistance-associated haplotypes for FaRCa1, FaRCg1, FaRMp1, and FaRPc2 were identified. These new core collections will allow breeders and researchers to more efficiently utilize the F. × ananassa collection. The core collections and other accessions of interest can be requested for research from the USDA-ARS NCGR via the Germplasm Resources Information Network (https://www.ars-grin.gov/).

References

[1] Darrow GM, ed. The Strawberry: History, Breeding and Physiology . Rinehart and Winston: Holt; 1966.

[2] Food and Agriculture Organization of the United Nations. 2018 Crops Report. 2018. http://www.fao.org/faostat/en/#data/QC.

[3] USDA-NASS. Strawberry production data. 2020. https://www.nass.usda.gov/.

[4] McCouch S, Baute GJ, Bradeen J et al. Feeding the future. Nature. 2013; 499: 23-4.

[5] Byrne PF, Volk GM, Gardner C et al. Sustaining the future of plant breeding: the critical role of the USDA-ARS national plant germplasm system. Crop Sci. 2018; 58: 451-68.

[6] De Beukelaer H, Davenport GF, Fack V . Core hunter 3: flexible core subset selection. BMC Bioinformatics. 2018; 19: 203.

[7] Odong TL, Jansen J, van Eeuwijk FA et al. Quality of core collections for effective utilisation of genetic resources review, discussion and interpretation. Theor Appl Genet. 2013; 126: 289-305.

[8] Langridge P, Waugh R . Harnessing the potential of germplasm collections. Nat Genet. 2019; 51: 200-1.

[9] Brown AHD. The use of plant genetic resources. In: Brown AHD, Frankel OH, Marshall DR, Williams JT, eds. The Case for Core Collections.New York, NY, U.S.A: Cambridge Univ. Press, 1989.

[10] Sharma S, Upadhyaya HD, Varshney RK et al. Pre-breeding for diversification of primary gene pool and genetic enhancement of grain legumes. Front Plant Sci. 2013; 4: 309.

[11] Edger PP, VanBuren R, Colle M et al. Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity . Gigascience. 2018; 7: gix124.

[12] Edger PP, Poorten TJ, VanBuren R et al. Origin and evolution of the octoploid strawberry genome. Nat Genet. 2019; 51: 541-7.

[13] Bassil NV, Davis TM, Zhang H et al. Development and preliminary evaluation of a 90 K axiom^® SNP array for the Allo-octoploid cultivated strawberry Fragaria × ananassa . BMC Genomics. 2015; 16: 155.

[14] Hardigan MA, Feldmann MJ, Lorant A et al. Genome Synteny has been conserved among the octoploid progenitors of cultivated strawberry over millions of years of evolution. Font Plant Sci. 2020; 10: 1789.

[15] Verma S, Davis TM, Zhang H et al. Development and evaluation of the axiom^® IStraw35 384HT array for the Allo-octoploid cultivated strawberry Fragaria × ananassa . Acta Hortic. 2016; 1156: 75-82.

[16] Chambers AH, Pillet J, Plotto A et al. Identification of a strawberry flavor gene candidate using an integrated genetic-genomic-analytical chemistry approach. BMC Genomics. 2014; 15: 217.

[17] Iezzoni AF, McFerson J, Luby J et al. RosBREED: bridging the chasm between discovery and application to enable DNA-informed breeding in rosaceous crops. Hortic Res. 2020; 7: 177.

[18] Perrotte J, Gaston A, Potier A et al. Narrowing down the single homoeologous FAPFRU locus controlling flowering in cultivated octoploid strawberry using a selective mapping strategy. Plant Biotechnol J. 2016; 14: 2176-89.

[19] Sánchez-Sevilla JF, Cruz-Rus E, Valpuesta V et al. Deciphering gamma-decalactone biosynthesis in strawberry fruit using a combination of genetic mapping, RNA-Seq and eQTL analyses. BMC Genomics. 2014; 15: 218.

[20] Verma S, Zurn JD, Salinas N et al. Clarifying sub-genomic positions of QTLs for flowering habit and fruit quality in U.S. strawberry (Fragaria × ananassa) breeding populations using pedigree-based QTL analysis . Hortic Res. 2017; 4: 17062.

[21] Zorrilla-Fontanesi Y, Rambla JL, Cabeza A et al. Genetic analysis of strawberry fruit aroma and identification of O-methyltransferase FaOMT as the locus controlling natural variation in mesifurane content. Plant Physiol. 2012; 159: 851-70.

[22] Anciro A, Mangandi J, Verma S et al. FaRCg1: a quantitative trait locus conferring resistance to Colletotrichum crown rot caused by Colletotrichum gloeosporioides in octoploid strawberry. Theor Appl Genet. 2018; 131: 2167-77.

[23] Mangandi J, Verma S, Osorio L et al. Pedigree-based analysis in a multiparental population of Octoploid strawberry reveals QTL alleles conferring resistance to Phytophthora cactorum. G3 (Bethesda). 2017; 7: 1707-19.

[24] Nelson JR, Verma S, Bassil NV et al. Discovery of three loci increasing resistance to charcoal rot caused by Macrophomina phaseolina in octoploid strawberry. G3 (Bethesda). 2021; 11: jkab037.

[25] Pincot DDA, Poorten TJ, Hardigan MA et al. Genome-wide association mapping uncovers Fw1, a dominant gene conferring resistance to Fusarium wilt in strawberry. G3 (Bethesda). 2018; 8: 1817-28.

[26] Salinas N, Verma S, Peres N et al. FaRCa1: a major subgenome-specific locus conferring resistance to Colletotrichum acutatum in strawberry. Theor Appl Genet. 2019; 132: 1109-20.

[27] Hardigan MA, Lorant A, Pincot DDA et al. Unraveling the complex hybrid ancestry and domestication history of cultivated strawberry. Mol Biol Evol. 2021; 38: 2285-305.

[28] Pincot DDA, Ledda M, Feldmann MJ et al. Social network analysis of the genealogy of strawberry: retracing the wild roots of heirloom and modern cultivars. G3 (Bethesda). 2021; 11: jkab015.

[29] Hardigan MA, Poorten TJ, Acharya CB et al. Domestication of temperate and coastal hybrids with distinct ancestral gene selection in octoploid strawberry. Plant Genome. 2018; 11: 180049.

[30] Gilmore BS, Bassil NV, Hummer KE . DNA extraction protocols from dormant buds for twelve woody plant genera. J Am Pomol Soc. 2011; 65: 201-7.

[31] R Core Team . R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2020. https://www.R-project.org/.

[32] Jombart T . Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008; 24: 1403-5.

[33] Jombart T, Ahmed I . Adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011; 27: 3070-1.

[34] Frichot E, François O . LEA: an R package for landscape and ecological association studies. Methods Ecol Evol. 2015; 6: 925-9.

[35] Frichot E, Mathieu F, Trouillon T et al. Fast and efficient estimation of individual ancestry coefficients. Genetics. 2014; 196: 973-83.

[36] Jakobsson M, Rosenberg NA . CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007; 23: 1801-6.

[37] Pritchard JK, Stephens M, Donnelly P . Inference of population structure using multilocus genotype data. Genetics. 2000; 155: 945-59.

[38] Earl DA, vonHoldt BM . STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012; 4: 359-61.

[39] Evanno G, Regnaut S, Goudet J . Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005; 14: 2611-20.

[40] Ramasamy RK, Ramasamy S, Bindroo BB et al. STRUCTURE PLOT: a program for drawing elegant STRUCTURE bar plots in user friendly interface. Springer Plus. 2014; 3: 431.

[41] Prevosti A, Ocaña J, Alonso G . Distances between populations of Drosophila subobscura, based on chromosome arrangement frequencies. Theor Appl Genet. 1975; 45: 231-41.

[42] Kamvar ZN, Brooks JC, Grünwald NJ . Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front Genet. 2015; 6: 208.

[43] Kamvar ZN, Tabima JF, Grünwald NJ . Poppr: an R package for genetic analysis of populations with clonal partially clonal, and/or sexual reproduction. PeerJ. 2014; 2: e281.

[44] Simpson EH . Measurement of diversity. Nature. 1949; 163: 688.

[45] Nei M . Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978; 89: 583-90.

[46] Goudet J . HIERFSTAT, a package for R to compute and test hierarchical F-statistics. Mol Ecol Notes. 2005; 5: 184-6.

[47] Odong TL, van Heerwaarden J, Jansen J et al. Statistical techniques for defining reference sets of accessions and microsatellite markers. Crop Sci. 2011; 51: 2401-11.

[48] Oh Y, Zurn JD, Bassil N et al. The strawberry DNA testing handbook. Hort Science. 2019; 54: 2267-70.

[49] van Dijk T, Pagliarani G, Pikunova A et al. Genomic rearrangements and signatures of breeding in the Allo-octoploid strawberry as revealed through an allele dose based SSR linkage map. BMC Plant Biol. 2014; 14: 55.

[50] Jones OR, Wang J . COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010; 10: 551-5.

[51] Browning SR, Browning BL . Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007; 81: 1084-97.

[52] Browning BL, Zhou Y, Browning SR . A one-penny imputed genome from next-generation reference panels. AJHG. 2018; 103: 338-48.

[53] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[54] Gil-Ariza DJ, Amaya I, López-Aranda JM et al. Impact of plant breeding on the genetic diversity of cultivated strawberry as revealed by expressed sequence tag-derived simple sequence repeat markers. J Amer Soc Hort Sci. 2009; 134: 337-47.

[55] Chu J, Zhao Y, Beier S et al. Suitability of single-nucleotide polymorphism arrays versus genotyping-by-sequencing for genebank genomics in wheat. Front Plant Sci. 2020; 11: 42.

[56] Geibel J, Reimer C, Weigend S et al. How array design creates SNP ascertainment bias. PLoS One. 2021; 16: e0245178.

[57] Zurn JD, Nyberg A, Montanari S et al. A new SSR fingerprinting set and its comparison to existing SSR- and SNP-based genotyping platforms to manage Pyrus germplasm resources. Tree Genet Genomes. 2020; 16: 72.

[58] Kalinowski ST . The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure. Heredity. 2011; 106: 625-32.

[59] Puechmaille SJ . The program STRUCTURE does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem. Mol Ecol Resour. 2016; 16: 608-27.

[60] Wang J . The computer program STRUCTURE for assigning individuals to populations: easy to use but easier to misuse. Mol Ecol Resour. 2017; 17: 981-90.

[61] Schwartz MK, McKelvey KS . Why sampling scheme matters: the effect of sampling scheme on landscape genetic results. Conservation genetics (Print). 2009; 10: 441-52.

[62] Hancock JF, Finn CE, Luby JJ et al. Reconstruction of the strawberry, Fragaria ×ananassa, using genotypes of F. virginiana and F. chiloensis. Science. 2010; 45: 1006-13.

[63] Hancock JF, Hokanson SC, Finn CE, Hummer KE . Introducing a supercore collection of wild octoploid strawberries. Acta Hortic. 2002; 567: 77-9.

[64] Zurn J, Ivors KL, Cole GS et al. Assessing cultivated strawberries and the Fragaria supercore for resistance to soilborne pathogens. J Am Pom Soc. 2020; 74: 18-23.

[65] Bassil N, Bidani A, Nyberg A et al. Microsatellite markers confirm identity of blueberry (Vaccinium spp.) plants in the USDA-ARS National Clonal Germplasm Repository collection . Genet Resour Crop Evol. 2020; 67: 393-409.

[66] Zurn JD, Carter KA, Yin MH et al. Validating blackberry seedling pedigrees and developing an improved multiplexed microsatellite fingerprinting set. J Amer Soc Hort Sci. 2018; 143: 381-90.

[67] Osorio LF, Gezan SA, Verma S et al. Independent validation of genomic prediction in strawberry over multiple cycles. Front Genet. 2021; 11: 596258.

[68] Gaut BS, Long AD . The lowdown on linkage disequilibrium. Plant Cell. 2003; 15: 1502-6.

[69] Nordborg M, Tavaré S . Linkage disequilibrium: what history has to tell us. Trends Genet. 2002; 18: 83-90.

[70] Yu G, Smith DK, Zhu H et al. GGTREE: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017; 8: 28-36.

[71] Shaw PD, Graham M, Kennedy J et al. Helium: visualization of large scale plant pedigrees. BMC Bioinformatics. 2014; 15: 259.