Genomic selection in plant breeding: Key factors shaping two decades of progress.

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

Novel SNP markers for flowering and seed quality traits in faba bean (Vicia faba L.): characterization and GWAS of a diversity panel.

Faba bean (Vicia faba L.) is a legume crop grown in diverse climates worldwide. It has a high potential for increased cultivation to meet the need for more plant-based proteins in human diets, a prerequisite for a more sustainable food production system. Characterization of diversity panels of crops can identify variation in and genetic markers for target traits of interest for plant breeding. In this work, we collected a diversity panel of 220 accessions of faba bean from around the world consisting of gene bank material and commercially available cultivars. The aims of this study were to quantify the phenotypic diversity in target traits to analyze the impact of breeding on these traits, and to identify genetic markers associated with traits through a genome-wide association study (GWAS). Characterization under field conditions at Nordic latitude across two years revealed a large genotypic variation and high broad-sense heritability for eleven agronomic and seed quality traits. Pairwise correlations showed that seed yield was positively correlated to plant height, number of seeds per plant, and days to maturity. Further, susceptibility to bean weevil damage was significantly higher for early flowering accessions and accessions with larger seeds. In this study, no yield penalty was found for higher seed protein content, but protein content was negatively correlated to starch content. Our results showed that while breeding advances in faba bean germplasm have resulted in increased yields and number of seeds per plant, they have also led to a selection pressure towards delayed onset of flowering and maturity. DArTseq genotyping identified 6,606 single nucleotide polymorphisms (SNPs) by alignment to the faba bean reference genome. These SNPs were used in a GWAS, revealing 51 novel SNP markers significantly associated with ten of the assessed traits. Three markers for days to flowering were found in predicted genes encoding proteins for which homologs in other plant species regulate flowering. Altogether, this work enriches the growing pool of phenotypic and genotypic data on faba bean as a valuable resource for developing efficient breeding strategies to expand crop cultivation.

Haplotype-tagged SNPs improve genomic prediction accuracy for Fusarium head blight resistance and yield-related traits in wheat.

KEY MESSAGE: Linkage disequilibrium (LD)-based haplotyping with subsequent SNP tagging improved the genomic prediction accuracy up to 0.07 and 0.092 for Fusarium head blight resistance and spike width, respectively, across six different models. Genomic prediction is a powerful tool to enhance genetic gain in plant breeding. However, the method is accompanied by various complications leading to low prediction accuracy. One of the major challenges arises from the complex dimensionality of marker data. To overcome this issue, we applied two pre-selection methods for SNP markers viz. LD-based haplotype-tagging and GWAS-based trait-linked marker identification. Six different models were tested with preselected SNPs to predict the genomic estimated breeding values (GEBVs) of four traits measured in 419 winter wheat genotypes. Ten different sets of haplotype-tagged SNPs were selected by adjusting the level of LD thresholds. In addition, various sets of trait-linked SNPs were identified with different scenarios from the training-test combined and only from the training populations. The BRR and RR-BLUP models developed from haplotype-tagged SNPs had a higher prediction accuracy for FHB and SPW by 0.07 and 0.092, respectively, compared to the corresponding models developed without marker pre-selection. The highest prediction accuracy for SPW and FHB was achieved with tagged SNPs pruned at weak LD thresholds (r2 < 0.5), while stringent LD was required for spike length (SPL) and flag leaf area (FLA). Trait-linked SNPs identified only from training populations failed to improve the prediction accuracy of the four studied traits. Pre-selection of SNPs via LD-based haplotype-tagging could play a vital role in optimizing genomic selection and reducing genotyping costs. Furthermore, the method could pave the way for developing low-cost genotyping methods through customized genotyping platforms targeting key SNP markers tagged to essential haplotype blocks.

The mosaic oat genome gives insights into a uniquely healthy cereal crop.

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

Splicing variation at a FLOWERING LOCUS C homeolog is associated with flowering time variation in the tetraploid Capsella bursa-pastoris.

The long-term fates of duplicate genes are well studied both empirically and theoretically, but how the short-term evolution of duplicate genes contributes to phenotypic variation is less well known. Here, we have studied the genetic basis of flowering time variation in the disomic tetraploid Capsella bursa-pastoris. We sequenced four duplicate candidate genes for flowering time and 10 background loci in samples from western Eurasia and China. Using a mixed-model approach that accounts for population structure, we found that polymorphisms at one homeolog of two candidate genes, FLOWERING LOCUS C (FLC) and CRYPTOCHROME1 (CRY1), were associated with natural flowering time variation. No potentially causative polymorphisms were found in the coding region of CRY1; however, at FLC two splice site polymorphisms were associated with early flowering. Accessions harboring nonconsensus splice sites expressed an alternatively spliced transcript or did not express this FLC homeolog. Our results are consistent with the function of FLC as a major repressor of flowering in Arabidopsis thaliana and imply that nonfunctionalization of duplicate genes could provide an important source of phenotypic variation.

New DArT markers for oat provide enhanced map coverage and global germplasm characterization.

BACKGROUND: Genomic discovery in oat and its application to oat improvement have been hindered by a lack of genetic markers common to different genetic maps, and by the difficulty of conducting whole-genome analysis using high-throughput markers. This study was intended to develop, characterize, and apply a large set of oat genetic markers based on Diversity Array Technology (DArT). RESULTS: Approximately 19,000 genomic clones were isolated from complexity-reduced genomic representations of pooled DNA samples from 60 oat varieties of global origin. These were screened on three discovery arrays, with more than 2000 polymorphic markers being identified for use in this study, and approximately 2700 potentially polymorphic markers being identified for use in future studies. DNA sequence was obtained for 2573 clones and assembled into a non-redundant set of 1770 contigs and singletons. Of these, 705 showed highly significant (Expectation < 10E-10) BLAST similarity to gene sequences in public databases. Based on marker scores in 80 recombinant inbred lines, 1010 new DArT markers were used to saturate and improve the 'Kanota' x 'Ogle' genetic map. DArT markers provided map coverage approximately equivalent to existing markers. After binning markers from similar clones, as well as those with 99% scoring similarity, a set of 1295 non-redundant markers was used to analyze genetic diversity in 182 accessions of cultivated oat of worldwide origin. Results of this analysis confirmed that major clusters of oat diversity are related to spring vs. winter type, and to the presence of major breeding programs within geographical regions. Secondary clusters revealed groups that were often related to known pedigree structure. CONCLUSION: These markers will provide a solid basis for future efforts in genomic discovery, comparative mapping, and the generation of an oat consensus map. They will also provide new opportunities for directed breeding of superior oat varieties, and guidance in the maintenance of oat genetic diversity.

The evolution and diversification of S-locus haplotypes in the Brassicaceae family.

Self-incompatibility (SI) in the Brassicaceae plant family is controlled by the SRK and SCR genes situated at the S locus. A large number of S haplotypes have been identified, mainly in cultivated species of the Brassica and Raphanus genera, but recently also in wild Arabidopsis species. Here, we used DNA sequences from the SRK and SCR genes of the wild Brassica species Brassica cretica, together with publicly available sequence data from other Brassicaceae species, to investigate the evolutionary relationships among S haplotypes in the Brassicaceae family. The results reveal that wild and cultivated Brassica species have similar levels of SRK diversity, indicating that domestication has had but a minor effect on S-locus diversity in Brassica. Our results also show that a common set of S haplotypes was present in the ancestor of the Brassica and Arabidopsis genera, that only a small number of haplotypes survived in the Brassica lineage after its separation from Arabidopsis, and that diversification within the two Brassica dominance classes occurred after the split between the two lineages. We also found indications that recombination may have occurred between the kinase domain of SRK and the SCR gene in Brassica.

Molecular population genetics of the SRK and SCR self-incompatibility genes in the wild plant species Brassica cretica (Brassicaceae).

Self-incompatibility (SI) in plants is a classic example of a trait evolving under strong frequency-dependent selection. As a consequence, population genetic theory predicts that the S locus, which controls SI, should maintain numerous alleles, display a high level of nucleotide diversity, and, in structured populations, show a lower level of among-population differentiation compared to neutral loci. Population-level investigations of DNA sequence variation at the S locus have recently been carried out in the genus Arabidopsis, largely confirming results from theoretical models of S-locus evolutionary dynamics, but no comparable studies have been done in wild Brassica species. In this study, we sequenced parts of the S-locus genes SRK and SCR, two tightly linked genes that are directly involved in the determination of SI specificity in samples from four natural populations of the wild species Brassica cretica. The amount and distribution of nucleotide diversity, as well as the frequency spectrum of putative functional haplotypes, observed at the S locus in B. cretica fit very well with expectations from theoretical models, providing strong evidence for frequency-dependent selection acting on the S locus in a wild Brassica species.

Unequal segregation of SRK alleles at the S locus in Brassica cretica.

In the Brassicaceae plant family, which includes the Arabidopsis and Brassica genera, self-incompatibility (SI) is controlled by genes at the S locus. Using experimental crosses, we studied the pattern of inheritance of S-locus alleles in the wild species Brassica cretica. Four full-sib families were established and unequal segregation of alleles at the SRK SI gene was found in one family. The segregation distortion acted in favour of a recessive (class II) allele and was best explained by some form of gametic-level selection. Our findings are discussed in the light of theoretical predictions of differential accumulation of deleterious mutations among S-locus alleles.

Polyploid speciation did not confer instant reproductive isolation in Capsella (Brassicaceae).

Polyploid formation is a major mode of sympatric speciation in flowering plants. Unlike other speciation processes, polyploidization is often assumed to confer instant reproductive isolation. Shared polymorphism across ploidy levels has therefore often been attributed to multiple polyploid origins, whereas the alternative hypothesis of introgressive hybridization has rarely been rigorously tested. Here, we sequence 12 nuclear loci representing 6 genes duplicated by polyploidy in 92 accessions of the tetraploid Capsella bursa-pastoris together with the corresponding loci in 21 accessions of its close diploid relative Capsella rubella. In C. bursa-pastoris accessions from western Eurasia, where the 2 species occur in partial sympatry, we find higher levels of nucleotide diversity than in accessions from eastern Eurasia, where C. rubella does not grow. Furthermore, haplotypes are shared across ploidy levels at 4 loci in western but not in eastern Eurasia. We test whether haplotype sharing is due to retention of ancestral polymorphism or due to hybridization and introgression using a coalescent-based isolation-with-migration model. In western but not in eastern Eurasia, there is evidence for unidirectional gene flow from C. rubella to C. bursa-pastoris. An independent estimate of the timing of dispersal of C. bursa-pastoris to eastern Eurasia indicates that it probably predated introgression. Our results show that polyploid speciation need not result in immediate and complete reproductive isolation, that postpolyploidization hybridization and introgression can contribute significantly to genetic variation in a newly formed polyploid, and that divergence population genetic analysis constitutes a powerful way of testing hypotheses on polyploid speciation.

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae).

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.

Bayesian inference of evolutionary history from chloroplast microsatellites in the cosmopolitan weed Capsella bursa-pastoris (Brassicaceae).

Besides showing an extraordinary degree of phenotypic variability, Capsella bursa-pastoris (Brassicaceae) is also one of the world's most common plant species and a serious weed in many countries. We have employed a coalescent-based Bayesian analysis of chloroplast microsatellite data to infer demographic and evolutionary parameters of this species. Two different demographic models applied to data from seven chloroplast microsatellite loci among 59 accessions show that the effective population size of C. bursa-pastoris is very small indicating a rapid expansion of the species, a result that is in accordance with fossil and historical data. Against this background, analysis of flowering time variation among accessions suggests that ecotypic differentiation in flowering time has occurred recently in the species' history. Finally, our results also indicate that mononucleotide repeat loci in the chloroplast genome can deteriorate in relatively short periods of evolutionary time.

Coalescence times and the Meselson effect in asexual eukaryotes.

In asexual eukaryotes, the two allelic gene copies at a locus are expected to become highly divergent as a result of the independent accumulation of mutations in the absence of segregation. If sexual reproduction was abandoned millions of generations ago, intra-individual allelic divergences can be significantly larger than in species that reproduce sexually. Owing to the disputed existence of truly ancient asexual species, this so-called 'Meselson effect' has been put forward as a means of confirming the complete loss of sexual reproduction. Very few attempts have, however, been made at quantifying the effect of sexual reproduction on the degree of divergence between gene copies in an asexual population. Here, I describe how asexual reproduction can be regarded as a special case of population subdivision. Using a slightly modified version of the standard two-deme structured coalescent, I derive the expected coalescence time for a pair of gene copies in an asexual population and show that the Meselson effect is compatible with low rates of sexual reproduction.