Plant mating systems shape patterns of genetic diversity and impact the long-term success of populations. As such, they are relevant to the design of seed collections aiming to maximise genetic diversity (e.g., germplasm conservation, ecological restoration). However, for most species, little is known empirically about how variation in mating systems and genetic diversity is distributed. We investigated the relationship between genetic diversity and mating systems in two functionally similar, co-occurring species of Hakea (Proteaceae), and evaluated the extent to which genetic diversity was captured in seeds. We genotyped hundreds of seedlings and mother plants via DArTseq, and developed novel implementations of two approaches to inferring the mating system from SNP data. A striking contrast in patterns of genetic diversity between H. sericea and H. teretifolia was revealed, consistent with a contrast in their mating systems. While both species had mixed mating systems, H. sericea was found to be habitually selfing, while H. teretifolia more evenly employed both selfing and outcrossing. In both species, seed collection schemes maximised genetic diversity by increasing the number of maternal lines and sites sampled, but twice as many sites were needed for the selfing species to capture equivalent levels of genetic variation at a regional scale.

The mating system is a central parameter of plant biology because it shapes their ecological and evolutionary properties. Therefore, determining ecological variables that influence the mating system is important for a deeper understanding of the functioning of plant populations. Here, using old concepts and recent statistical developments, we propose a new statistical tool to make inferences about ecological determinants of outcrossing in natural plant populations. The method requires codominant genotypes of seeds collected from maternal plants within different locations. Using extensive computer simulations, we demonstrated that the method is robust to the issues expected for real-world data, including the Wahlund effect, inbreeding and genotyping errors such as allele dropout and allele misclassification. Furthermore, we showed that the estimates of ecological effects and outcrossing rates can be severely biased if genotyping errors and genetic differentiation are not treated explicitly. Application of the new method to the case study of a dioecious tree (Taxus baccata) allowed revealing that female trees that grow in lower local densities have a greater tendency towards mating with relatives. Moreover, we also demonstrated that biparental inbreeding is higher in populations that are characterized by a longer mean distance between trees and a smaller mean trunk perimeter. We found these results to agree with both the theoretical predictions and the history of English yew.