Background: In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation. Results: We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed $$f_{4}$$ f 4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds. Conclusions: Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Cultivated peanut (Arachis hypogaea L.) is a widely grown oilseed crop worldwide; however, the events leading to its origin and diversification are not fully understood. Here by combining chloroplast and whole-genome sequence data from a large germplasm collection, we show that the two subspecies of A. hypogaea (hypogaea and fastigiata) likely arose from distinct allopolyploidization and domestication events. Peanut genetic clusters were then differentiated in relation to dissemination routes and breeding efforts. A combination of linkage mapping and genome-wide association studies allowed us to characterize genes and genomic regions related to main peanut morpho-agronomic traits, namely flowering pattern, inner tegument color, growth habit, pod/seed weight and oil content. Together, our findings shed light on the evolutionary history and phenotypic diversification of peanuts and might be of broad interest to plant breeders.

Almond [Prunus dulcisMiller (D.A. Webb)] is the main tree nut species worldwide. Here, genotyping-by-sequencing (GBS) was applied to 149 almond cultivars from the ex situ collections of the Italian Council for Agricultural Research (CREA) and the Spanish National Research Council (CSIC), leading to the detection of 93,119 single-nucleotide polymorphisms (SNPs). The study of population structure outlined four distinct genetic groups and highlighted diversification between the Mediterranean and Californian gene pools. Data on SNP diversity and runs of homozygosity (ROHs) allowed the definition of kinship, inbreeding, and linkage disequilibrium (LD) decay in almond cultivated germplasm. Four-year phenotypic observations, gathered on 98 cultivars of the CREA collection, were used to perform a genome-wide association study (GWAS) and, for the first time in a crop species, homozygosity mapping (HM), resulting in the identification of genomic associations with nut, shell, and seed weight. Both GWAS and HM suggested that loci controlling nut and seed weight are mostly independent. Overall, this study provides insights on the almond cultivation history and delivers information of major interest for almond genetics and breeding. In a broader perspective, our results encourage the use of ROHs in crop science to estimate inbreeding, choose parental combinations minimizing the risk of inbreeding depression, and identify genomic footprints of selection for specific traits.

Drought stress is one of the key plant stresses reducing grain yield in cereal crops worldwide. Although it is not a breeding target in Northern Europe, the changing climate and the drought of 2018 have increased its significance in the region. A key challenge, therefore, is to identify novel germplasm with higher drought tolerance, a task that will require continuous characterization of a large number of genotypes. The aim of this work was to assess if phenotyping systems with low-cost consumer-grade digital cameras can be used to characterize germplasm for drought tolerance. To achieve this goal, we built a proximal phenotyping cart mounted with digital cameras and evaluated it by characterizing 142 winter wheat genotypes for drought tolerance under field conditions. The same genotypes were additionally characterized for seedling stage traits by imaging under controlled growth conditions. The analysis revealed that under field conditions, plant biomass, relative growth rates, and Normalized Difference Vegetation Index (NDVI) from different growth stages estimated by imaging were significantly correlated to drought tolerance. Under controlled growth conditions, root count at the seedling stage evaluated by imaging was significantly correlated to adult plant drought tolerance observed in the field. Random forest models were trained by integrating measurements from field and controlled conditions and revealed that plant biomass and relative growth rates at key plant growth stages are important predictors of drought tolerance. Thus, based on the results, it can be concluded that the consumer-grade cameras can be key components of affordable automated phenotyping systems to accelerate pre-breeding for drought tolerance.

Genotyping by sequencing (GBS) was used to analyze relationships among cowpea and asparagus bean landraces from southern Italy and to assess the utility of this technology to study taxonomy in a wider panel, including V. unguiculata cultigroups, subspecies, and other Vigna species. The analysis of SNPs derived from GBS highlighted that, among the cowpea landraces, the African samples were separated from the other material, while, for the Italian landraces, a certain clustering depending on seed color/pattern was observed in the dendrogram. When examining the V. unguiculata species complex, a clear separation between the two groups of wild subspecies, i.e., the allogamous wild perennials and the perennial out/inbreds, could be observed, the former representing the more ancestral wild progenitors of V. unguiculata. The species V. vexillata appeared more closely related to V. unguiculata than to the other Vigna species analyzed.