Admixture analysis using genotyping-by-sequencing reveals genetic relatedness and parental lineage distribution in highbush blueberry genotypes and cross derivatives

Blueberries (Vaccinium section Cyanococcus) are perennial shrubs widely cultivated for their edible fruits. In this study, we performed admixture and genetic relatedness analysis of northern highbush (NHB, primarily V. corymbosum) and southern highbush (SHB, V. corymbosum introgressed with V. darrowii, V. virgatum, or V. tenellum) blueberry genotypes, and progenies of the BNJ16-5 cross (V. corymbosum × V. darrowii). Using genotyping-by-sequencing (GBS), we generated more than 334 million reads (75 bp). The GBS reads were aligned to the V. corymbosum cv. Draper v1.0 reference genome sequence, and ~2.8 million reads were successfully mapped. From the alignments, we identified 2,244,039 single-nucleotide polymorphisms, which were used for principal component, haplotype, and admixture analysis. Principal component analysis revealed three main groups: (1) NHB cultivars, (2) SHB cultivars, and (3) BNJ16-5 progenies. The overall fixation index (F_ST) and nucleotide diversity for NHB and SHB cultivars indicated wide genetic differentiation, and haplotype analysis revealed that SHB cultivars are more genetically diverse than NHB cultivars. The admixture analysis identified a mixture of various lineages of parental genomic introgression. This study demonstrated the effectiveness of GBS-derived single-nucleotide polymorphism markers in genetic and admixture analyses to reveal genetic relatedness and to examine parental lineages in blueberry, which may be useful for future breeding plans.