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 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.

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.

The recent outbreak of the Olive Quick Decline Syndrome (OQDS), caused by Xylella fastidiosa subsp. pauca (Xf), is dramatically altering ecosystem services in the peninsula of Salento (Apulia Region, southeastern Italy). Here we report the accomplishment of several exploratory missions in the Salento area, resulting in the identification of thirty paucisymptomatic or asymptomatic plants in olive orchards severely affected by the OQDS. The genetic profiles of such putatively resistant plants (PRPs), assessed by a selection of ten simple sequence repeat (SSR) markers, were compared with those of 141 Mediterranean cultivars. Most (23) PRPs formed a genetic cluster (K1) with 22 Italian cultivars, including ‘Leccino’ and ‘FS17’, previously reported as resistant to Xf. The remaining PRPs displayed relatedness with genetically differentiated germplasm, including a cluster of Tunisian cultivars. Markedly lower colonization levels were observed in PRPs of the cluster K1 with respect to control plants. Field evaluation of four cultivars related to PRPs allowed the definition of partial resistance in the genotypes ‘Frantoio’ and ‘Nocellara Messinese’. Some of the PRPs identified in this study might be exploited in cultivation, or as parental clones of breeding programs. In addition, our results indicate the possibility to characterize resistance to Xf in cultivars genetically related to PRPs.

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.

Almond [Prunus dulcis Miller (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.