Genetic diversity and population structure of superior shea trees (Vitellaria paradoxa subsp. paradoxa) using SNP markers for the establishment of a core collection in Côte d’Ivoire

Abstract Background The shea tree is a well-known carbon sink in Africa that requires a sustainable conservation of its gene pool. However, the genetic structure of its population is not well studied, especially in Côte d’Ivoire. In this study, 333 superior shea tree genotypes conserved in situ in Côte d’Ivoire were collected and genotyped with the aim of investigating its genetic diversity and population structure to facilitate suitable conservation and support future breeding efforts to adapt to climate change effects. Results A total of 7,559 filtered high-quality single nucleotide polymorphisms (SNPs) were identified using the genotyping by sequencing technology. The gene diversity (HE) ranged between 0.1 to 0.5 with an average of 0.26, while the polymorphism information content (PIC) value ranged between 0.1 to 0.5 with an average of 0.24, indicating a moderate genetic diversity among the studied genotypes. The population structure model classified the 333 genotypes into three genetic groups (GP1, GP2, and GP3). GP1 contained shea trees that mainly originated from the Poro, Tchologo, and Hambol districts, while GP2 and GP3 contained shea trees collected from the Bagoué district. Analysis of molecular variance (AMOVA) identified 55% variance within populations and 45% variance within individuals, indicating a very low genetic differentiation (or very high gene exchange) between these three groups (F ST  = 0.004, gene flow Nm = 59.02). Morphologically, GP1 displayed spreading tree growth habit, oval nut shape, higher mean nut weight (10.62 g), wide leaf (limb width = 4.63 cm), and small trunk size (trunk circumference = 133.4 cm). Meanwhile, GP2 and GP3 showed similar morphological characteristics: erect and spreading tree growth habit, ovoid nut shape, lower mean nut weight (GP2: 8.89 g; GP3: 8.36 g), thin leaf (limb width = 4.45 cm), and large trunk size ( GP2: 160.5 cm, GP3: 149.1 cm). A core set of 100 superior shea trees, representing 30% of the original population size and including individuals from all four study districts, was proposed using the “maximum length sub-tree function” in DARwin v. 6.0.21. Conclusion These findings provide new knowledge of the genetic diversity and population structure of Ivorian shea tree genetic resources for the design of effective collection and conservation strategies for the efficient use of inbreeding.