Genetic contributions and their optimization

Evolutionary computing to assemble standing genetic diversity and achieve long-term genetic gain

Loss of genetic diversity in elite crop breeding pools can severely limit long-term genetic gains and limit ability to make gains in new traits, like heat tolerance, that are becoming important as the climate changes. Here, we investigate and propose potential breeding program applications of optimal haplotype stacking (OHS), a selection method that retains useful diversity in the population. OHS selects sets of candidates containing, between them, haplotype segments with very high segment breeding values for the target trait. We compared the performance of OHS, a similar method called optimal population value (OPV), truncation selection on genomic estimated breeding values (GEBVs), and optimal contribution selection (OCS) in stochastic simulations of recurrent selection on founder wheat genotypes. After 100 generations of intercrossing and selection, OCS and truncation selection had exhausted the genetic diversity, while considerable diversity remained in the OHS population. Gain under OHS in these simulations ultimately exceeded that from truncation selection or OCS. OHS achieved faster gains when the population size was small, with many progeny per cross. A promising hybrid strategy, involving a single cycle of OHS in the first generation followed by recurrent truncation selection, substantially improved long-term gain compared with truncation selection and performed similarly to OCS. The results of this study provide initial insights into where OHS could be incorporated into breeding programs.

Optimizing selection based on BLUPs or BLUEs in multiple sets of genotypes differing in their population parameters

KEY MESSAGE

Selection response in truncation selection across multiple sets of candidates hinges on their post-selection proportions, which can deviate grossly from their initial proportions. For BLUPs, using a uniform threshold for all candidates maximizes the selection response, irrespective of differences in population parameters. Plant breeding programs typically involve multiple families from either the same or different populations, varying in means, genetic variances and prediction accuracy of BLUPs or BLUEs for true genetic values (TGVs) of candidates. We extend the classical breeder's equation for truncation selection from single to multiple sets of genotypes, indicating that the expected overall selection response ( Δ G Tot ) for TGVs depends on the selection response within individual sets and their post-selection proportions. For BLUEs, we show that maximizingΔ G Tot requires thresholds optimally tailored for each set, contingent on their population parameters. For BLUPs, we prove thatΔ G Tot is maximized by applying a uniform threshold across all candidates from all sets. We provide explicit formulas for the origin of the selected candidates from different sets and show that their proportions before and after selection can differ substantially, especially for sets with inferior properties and low proportion. We discuss implications of these results for (a) optimum allocation of resources to training and prediction sets and (b) the need to counteract narrowing the genetic variation under genomic selection. For genomic selection of hybrids based on BLUPs of GCA of their parent lines, selecting distinct proportions in the two parent populations can be advantageous, if these differ substantially in the variance and/or prediction accuracy of GCA. Our study sheds light on the complex interplay of selection thresholds and population parameters for the selection response in plant breeding programs, offering insights into the effective resource management and prudent application of genomic selection for improved crop development.

Efficient genomics-based 'end-to-end' selective tree breeding framework

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an "end-to-end" selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals' breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits' correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method's superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Assessing the potential of genetic resource introduction into elite germplasm: a collaborative multiparental population for flint maize

KEY MESSAGE

Implementing a collaborative pre-breeding multi-parental population efficiently identifies promising donor x elite pairs to enrich the flint maize elite germplasm. Genetic diversity is crucial for maintaining genetic gains and ensuring breeding programs' long-term success. In a closed breeding program, selection inevitably leads to a loss of genetic diversity. While managing diversity can delay this loss, introducing external sources of diversity is necessary to bring back favorable genetic variation. Genetic resources exhibit greater diversity than elite materials, but their lower performance levels hinder their use. This is the case for European flint maize, for which elite germplasm has incorporated only a limited portion of the diversity available in landraces. To enrich the diversity of this elite genetic pool, we established an original cooperative maize bridging population that involves crosses between private elite materials and diversity donors to create improved genotypes that will facilitate the incorporation of original favorable variations. Twenty donor × elite BC1S2 families were created and phenotyped for hybrid value for yield related traits. Crosses showed contrasted means and variances and therefore contrasted potential in terms of selection as measured by their usefulness criterion (UC). Average expected mean performance gain over the initial elite material was 5%. The most promising donor for each elite line was identified. Results also suggest that one more generation, i.e., 3 in total, of crossing to the elite is required to fully exploit the potential of a donor. Altogether, our results support the usefulness of incorporating genetic resources into elite flint maize. They call for further effort to create fixed diversity donors and identify those most suitable for each elite program.

Enhancing climate change resilience in agricultural crops

Climate change threatens global food and nutritional security through negative effects on crop growth and agricultural productivity. Many countries have adopted ambitious climate change mitigation and adaptation targets that will exacerbate the problem, as they require significant changes in current agri-food systems. In this review, we provide a roadmap for improved crop production that encompasses the effective transfer of current knowledge into plant breeding and crop management strategies that will underpin sustainable agriculture intensification and climate resilience. We identify the main problem areas and highlight outstanding questions and potential solutions that can be applied to mitigate the impacts of climate change on crop growth and productivity. Although translation of scientific advances into crop production lags far behind current scientific knowledge and technology, we consider that a holistic approach, combining disciplines in collaborative efforts, can drive better connections between research, policy, and the needs of society.

Comparison of genomic-enabled cross selection criteria for the improvement of inbred line breeding populations

A crucial step in inbred plant breeding is the choice of mating design to derive high-performing inbred varieties while also maintaining a competitive breeding population to secure sufficient genetic gain in future generations. In practice, the mating design usually relies on crosses involving the best parental inbred lines to ensure high mean progeny performance. This excludes crosses involving lower performing but more complementary parents in terms of favorable alleles. We predicted the ability of crosses to produce putative outstanding progenies (high mean and high variance progeny distribution) using genomic prediction models. This study compared the benefits and drawbacks of 7 genomic cross selection criteria (CSC) in terms of genetic gain for 1 trait and genetic diversity in the next generation. Six CSC were already published, and we propose an improved CSC that can estimate the proportion of progeny above a threshold defined for the whole mating plan. We simulated mating designs optimized using different CSC. The 835 elite parents came from a real breeding program and were evaluated between 2000 and 2016. We applied constraints on parental contributions and genetic similarities between selected parents according to usual breeder practices. Our results showed that CSC based on progeny variance estimation increased the genetic value of superior progenies by up to 5% in the next generation compared to CSC based on the progeny mean estimation (i.e. parental genetic values) alone. It also increased the genetic gain (up to 4%) and/or maintained more genetic diversity at QTLs (up to 4% more genic variance when the marker effects were perfectly estimated).

A simulation study comparing advanced marker-assisted selection with genomic selection in tree breeding programs

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future.

Genomic prediction in hybrid breeding: II. Reciprocal recurrent genomic selection with full-sib and half-sib families

KEY MESSAGE

Genomic prediction of GCA effects based on model training with full-sib rather than half-sib families yields higher short- and long-term selection gain in reciprocal recurrent genomic selection for hybrid breeding, if SCA effects are important. Reciprocal recurrent genomic selection (RRGS) is a powerful tool for ensuring sustainable selection progress in hybrid breeding. For training the statistical model, one can use half-sib (HS) or full-sib (FS) families produced by inter-population crosses of candidates from the two parent populations. Our objective was to compare HS-RRGS and FS-RRGS for the cumulative selection gain ([Formula: see text]), the genetic, GCA and SCA variances ([Formula: see text],[Formula: see text], [Formula: see text]) of the hybrid population, and prediction accuracy ([Formula: see text]) for GCA effects across cycles. Using SNP data from maize and wheat, we simulated RRGS programs over 10 cycles, each consisting of four sub-cycles with genomic selection of [Formula: see text] out of 950 candidates in each parent population. Scenarios differed for heritability [Formula: see text] and the proportion [Formula: see text] of traits, training set (TS) size ([Formula: see text]), and maize vs. wheat. Curves of [Formula: see text] over selection cycles showed no crossing of both methods. If [Formula: see text] was high, [Formula: see text] was generally higher for FS-RRGS than HS-RRGS due to higher [Formula: see text]. In contrast, HS-RRGS was superior or on par with FS-RRGS, if [Formula: see text] or [Formula: see text] and [Formula: see text] were low. [Formula: see text] showed a steeper increase and higher selection limit for scenarios with low [Formula: see text], high [Formula: see text] and large [Formula: see text]. [Formula: see text] and even more so [Formula: see text] decreased rapidly over cycles for both methods due to the high selection intensity and the role of the Bulmer effect for reducing [Formula: see text]. Since the TS for FS-RRGS can additionally be used for hybrid prediction, we recommend this method for achieving simultaneously the two major goals in hybrid breeding: population improvement and cultivar development.

Impact of kinship matrices on genetic gain and inbreeding with optimum contribution selection in a genomic dairy cattle breeding program

BACKGROUND

Genomic selection has increased genetic gain in dairy cattle, but in some cases it has resulted in higher inbreeding rates. Therefore, there is need for research on efficient management of inbreeding in genomically-selected dairy cattle populations, especially for local breeds with a small population size. Optimum contribution selection (OCS) minimizes the increase in average kinship while it maximizes genetic gain. However, there is no consensus on how to construct the kinship matrix used for OCS and whether it should be based on pedigree or genomic information. VanRaden's method 1 (VR1) is a genomic relationship matrix in which centered genotype scores are scaled with the sum of 2p(1-p) where p is the reference allele frequency at each locus, and VanRaden's method 2 (VR2) scales each locus with 2p(1-p), thereby giving greater weight to loci with a low minor allele frequency. We compared the effects of nine kinship matrices on genetic gain, kinship, inbreeding, genetic diversity, and minor allele frequency when applying OCS in a simulated small dairy cattle population. We used VR1 and VR2, each using base animals, all genotyped animals, and the current generation of animals to compute reference allele frequencies. We also set the reference allele frequencies to 0.5 for VR1 and the pedigree-based relationship matrix. We constrained OCS to select a fixed number of sires per generation for all scenarios. Efficiency of the different matrices were compared by calculating the rate of genetic gain for a given rate of increase in average kinship.

RESULTS

We found that: (i) genomic relationships were more efficient than pedigree-based relationships at managing inbreeding, (ii) reference allele frequencies computed from base animals were more efficient compared to reference allele frequencies computed from recent animals, and (iii) VR1 was slightly more efficient than VR2, but the difference was not statistically significant.

CONCLUSIONS

Using genomic relationships for OCS realizes more genetic gain for a given amount of kinship and inbreeding than using pedigree relationships when the number of sires is fixed. For a small genomic dairy cattle breeding program, we recommend that the implementation of OCS uses VR1 with reference allele frequencies estimated either from base animals or old genotyped animals.

Clonal diploid and autopolyploid breeding strategies to harness heterosis: insights from stochastic simulation

KEY MESSAGE

Reciprocal recurrent selection sometimes increases genetic gain per unit cost in clonal diploids with heterosis due to dominance, but it typically does not benefit autopolyploids. Breeding can change the dominance as well as additive genetic value of populations, thus utilizing heterosis. A common hybrid breeding strategy is reciprocal recurrent selection (RRS), in which parents of hybrids are typically recycled within pools based on general combining ability. However, the relative performances of RRS and other breeding strategies have not been thoroughly compared. RRS can have relatively increased costs and longer cycle lengths, but these are sometimes outweighed by its ability to harness heterosis due to dominance. Here, we used stochastic simulation to compare genetic gain per unit cost of RRS, terminal crossing, recurrent selection on breeding value, and recurrent selection on cross performance considering different amounts of population heterosis due to dominance, relative cycle lengths, time horizons, estimation methods, selection intensities, and ploidy levels. In diploids with phenotypic selection at high intensity, whether RRS was the optimal breeding strategy depended on the initial population heterosis. However, in diploids with rapid-cycling genomic selection at high intensity, RRS was the optimal breeding strategy after 50 years over almost all amounts of initial population heterosis under the study assumptions. Diploid RRS required more population heterosis to outperform other strategies as its relative cycle length increased and as selection intensity and time horizon decreased. The optimal strategy depended on selection intensity, a proxy for inbreeding rate. Use of diploid fully inbred parents vs. outbred parents with RRS typically did not affect genetic gain. In autopolyploids, RRS typically did not outperform one-pool strategies regardless of the initial population heterosis.

A method for partitioning trends in genetic mean and variance to understand breeding practices

BACKGROUND

In breeding programmes, the observed genetic change is a sum of the contributions of different selection paths represented by groups of individuals. Quantifying these sources of genetic change is essential for identifying the key breeding actions and optimizing breeding programmes. However, it is difficult to disentangle the contribution of individual paths due to the inherent complexity of breeding programmes. Here we extend the previously developed method for partitioning genetic mean by paths of selection to work both with the mean and variance of breeding values.

METHODS

First, we extended the partitioning method to quantify the contribution of different paths to genetic variance assuming that the breeding values are known. Second, we combined the partitioning method with the Markov Chain Monte Carlo approach to draw samples from the posterior distribution of breeding values and use these samples for computing the point and interval estimates of partitions for the genetic mean and variance. We implemented the method in the R package AlphaPart. We demonstrated the method with a simulated cattle breeding programme.

RESULTS

We show how to quantify the contribution of different groups of individuals to genetic mean and variance and that the contributions of different selection paths to genetic variance are not necessarily independent. Finally, we observed that the partitioning method under the pedigree-based model has some limitations, which suggests the need for a genomic extension.

CONCLUSIONS

We presented a partitioning method to quantify sources of change in genetic mean and variance in breeding programmes. The method can help breeders and researchers understand the dynamics in genetic mean and variance in a breeding programme. The developed method for partitioning genetic mean and variance is a powerful method for understanding how different selection paths interact within a breeding programme and how they can be optimised.

Assessment of long-term trends in genetic mean and variance after the introduction of genomic selection in layers: a simulation study

Nucleus-based breeding programs are characterized by intense selection that results in high genetic gain, which inevitably means reduction of genetic variation in the breeding population. Therefore, genetic variation in such breeding systems is typically managed systematically, for example, by avoiding mating the closest relatives to limit progeny inbreeding. However, intense selection requires maximum effort to make such breeding programs sustainable in the long-term. The objective of this study was to use simulation to evaluate the long-term impact of genomic selection on genetic mean and variance in an intense layer chicken breeding program. We developed a large-scale stochastic simulation of an intense layer chicken breeding program to compare conventional truncation selection to genomic truncation selection optimized with either minimization of progeny inbreeding or full-scale optimal contribution selection. We compared the programs in terms of genetic mean, genic variance, conversion efficiency, rate of inbreeding, effective population size, and accuracy of selection. Our results confirmed that genomic truncation selection has immediate benefits compared to conventional truncation selection in all specified metrics. A simple minimization of progeny inbreeding after genomic truncation selection did not provide any significant improvements. Optimal contribution selection was successful in having better conversion efficiency and effective population size compared to genomic truncation selection, but it must be fine-tuned for balance between loss of genetic variance and genetic gain. In our simulation, we measured this balance using trigonometric penalty degrees between truncation selection and a balanced solution and concluded that the best results were between 45° and 65°. This balance is specific to the breeding program and depends on how much immediate genetic gain a breeding program may risk vs. save for the future. Furthermore, our results show that the persistence of accuracy is better with optimal contribution selection compared to truncation selection. In general, our results show that optimal contribution selection can ensure long-term success in intensive breeding programs using genomic selection.

Genomics in animal breeding from the perspectives of matrices and molecules

BACKGROUND

This paper describes genomics from two perspectives that are in use in animal breeding and genetics: a statistical perspective concentrating on models for estimating breeding values, and a sequence perspective concentrating on the function of DNA molecules.

MAIN BODY

This paper reviews the development of genomics in animal breeding and speculates on its future from these two perspectives. From the statistical perspective, genomic data are large sets of markers of ancestry; animal breeding makes use of them while remaining agnostic about their function. From the sequence perspective, genomic data are a source of causative variants; what animal breeding needs is to identify and make use of them.

CONCLUSION

The statistical perspective, in the form of genomic selection, is the more applicable in contemporary breeding. Animal genomics researchers using from the sequence perspective are still working towards this the isolation of causative variants, equipped with new technologies but continuing a decades-long line of research.

Genomic prediction and selection response for grain yield in safflower

In plant breeding programs, multiple traits are recorded in each trial, and the traits are often correlated. Correlated traits can be incorporated into genomic selection models, especially for traits with low heritability, to improve prediction accuracy. In this study, we investigated the genetic correlation between important agronomic traits in safflower. We observed the moderate genetic correlations between grain yield (GY) and plant height (PH, 0.272–0.531), and low correlations between grain yield and days to flowering (DF, −0.157–0.201). A 4%–20% prediction accuracy improvement for grain yield was achieved when plant height was included in both training and validation sets with multivariate models. We further explored the selection responses for grain yield by selecting the top 20% of lines based on different selection indices. Selection responses for grain yield varied across sites. Simultaneous selection for grain yield and seed oil content (OL) showed positive gains across all sites with equal weights for both grain yield and oil content. Combining g×E interaction into genomic selection (GS) led to more balanced selection responses across sites. In conclusion, genomic selection is a valuable breeding tool for breeding high grain yield, oil content, and highly adaptable safflower varieties.

Genomic selection strategies for clonally propagated crops

For genomic selection in clonally propagated crops with diploid (-like) meiotic behavior to be effective, crossing parents should be selected based on genomic predicted cross-performance unless dominance is negligible. For genomic selection (GS) in clonal breeding programs to be effective, parents should be selected based on genomic predicted cross-performance unless dominance is negligible. Genomic prediction of cross-performance enables efficient exploitation of the additive and dominance value simultaneously. Here, we compared different GS strategies for clonally propagated crops with diploid (-like) meiotic behavior, using strawberry as an example. We used stochastic simulation to evaluate six combinations of three breeding programs and two parent selection methods. The three breeding programs included (1) a breeding program that introduced GS in the first clonal stage, and (2) two variations of a two-part breeding program with one and three crossing cycles per year, respectively. The two parent selection methods were (1) parent selection based on genomic estimated breeding values (GEBVs) and (2) parent selection based on genomic predicted cross-performance (GPCP). Selection of parents based on GPCP produced faster genetic gain than selection of parents based on GEBVs because it reduced inbreeding when the dominance degree increased. The two-part breeding programs with one and three crossing cycles per year using GPCP always produced the most genetic gain unless dominance was negligible. We conclude that (1) in clonal breeding programs with GS, parents should be selected based on GPCP, and (2) a two-part breeding program with parent selection based on GPCP to rapidly drive population improvement has great potential to improve breeding clonally propagated crops.

Developing Methods for Maintaining Genetic Diversity in Novel Aquaculture Species: The Case of Seriola lalandi

Developing sound breeding programs for aquaculture species may be challenging when matings cannot be controlled due to communal spawning. We developed a genotyping-by-sequencing marker panel of 300 SNPs for parentage testing and sex determination by using data from an in-house reference genome as well as a 90 K SNP genotyping array based on different populations of yellowtail kingfish (Seriola lalandi). The minimum and maximum distance between adjacent marker pairs were 0.7 Mb and 13 Mb, respectively, with an average marker spacing of 2 Mb. Weak evidence of the linkage disequilibrium between adjacent marker pairs was found. The results showed high panel performance for parental assignment, with probability exclusion values equaling 1. The rate of false positives when using cross-population data was null. A skewed distribution of genetic contributions by dominant females was observed, thus increasing the risk of higher rates of inbreeding in subsequent captive generations when no parentage data are used. All these results are discussed in the context of breeding program design, using this marker panel to increase the sustainability of this aquaculture resource.

Accuracy of Selection in Early Generations of Field Pea Breeding Increases by Exploiting the Information Contained in Correlated Traits

Accuracy of predicted breeding values (PBV) for low heritability traits may be increased in early generations by exploiting the information available in correlated traits. We compared the accuracy of PBV for 10 correlated traits with low to medium narrow-sense heritability (h2) in a genetically diverse field pea (Pisum sativum L.) population after univariate or multivariate linear mixed model (MLMM) analysis with pedigree information. In the contra-season, we crossed and selfed S1 parent plants, and in the main season we evaluated spaced plants of S0 cross progeny and S2+ (S2 or higher) self progeny of parent plants for the 10 traits. Stem strength traits included stem buckling (SB) (h2 = 0.05), compressed stem thickness (CST) (h2 = 0.12), internode length (IL) (h2 = 0.61) and angle of the main stem above horizontal at first flower (EAngle) (h2 = 0.46). Significant genetic correlations of the additive effects occurred between SB and CST (0.61), IL and EAngle (-0.90) and IL and CST (-0.36). The average accuracy of PBVs in S0 progeny increased from 0.799 to 0.841 and in S2+ progeny increased from 0.835 to 0.875 in univariate vs MLMM, respectively. An optimized mating design was constructed with optimal contribution selection based on an index of PBV for the 10 traits, and predicted genetic gain in the next cycle ranged from 1.4% (SB), 5.0% (CST), 10.5% (EAngle) and -10.5% (IL), with low achieved parental coancestry of 0.12. MLMM improved the potential genetic gain in annual cycles of early generation selection in field pea by increasing the accuracy of PBV.

Optimal Contribution Selection Improves the Rate of Genetic Gain in Grain Yield and Yield Stability in Spring Canola in Australia and Canada

Crop breeding must achieve higher rates of genetic gain in grain yield (GY) and yield stability to meet future food demands in a changing climate. Optimal contributions selection (OCS) based on an index of key economic traits should increase the rate of genetic gain while minimising population inbreeding. Here we apply OCS in a global spring oilseed rape (canola) breeding program during three cycles of S_0,1 family selection in 2016, 2018, and 2020, with several field trials per cycle in Australia and Canada. Economic weights in the index promoted high GY, seed oil, protein in meal, and Phoma stem canker (blackleg) disease resistance while maintaining plant height, flowering time, oleic acid, and seed size and decreasing glucosinolate content. After factor analytic modelling of the genotype-by-environment interaction for the additive effects, the linear rate of genetic gain in GY across cycles was 0.059 or 0.087 t ha^-1 y^-1 (2.9% or 4.3% y^-1) based on genotype scores for the first factor (f₁) expressed in trait units or average predicted breeding values across environments, respectively. Both GY and yield stability, defined as the root-mean-square deviation from the regression line associated with f₁, were predicted to improve in the next cycle with a low achieved mean parental coancestry (0.087). These methods achieved rapid genetic gain in GY and other traits and are predicted to improve yield stability across global spring canola environments.

Fine-tuning genomic and pedigree inbreeding rates in equine population with a deep and reliable stud book: the case of the Pura Raza Española horse

BACKGROUND

Estimating inbreeding, which is omnipresent and inevitable in livestock populations, is a primary goal for management and animal breeding especially for those interested in mitigating the negative consequences of inbreeding. Inbreeding coefficients have been historically estimated by using pedigree information; however, over the last decade, genome-base inbreeding coefficients have come to the forefront in this field. The Pura Raza Española (PRE) horse is an autochthonous Spanish horse breed which has been recognised since 1912. The total PRE population (344,718 horses) was used to estimate Classical (F), Ballou's ancestral, Kalinowski's ancestral, Kalinowski's new and the ancestral history coefficient values. In addition, genotypic data from a selected population of 805 PRE individuals was used to determine the individual inbreeding coefficient using SNP-by-SNP-based techniques (methods of moments -F_HOM-, the diagonal elements of the genomic -F_G-, and hybrid matrixes -F_H-) and ROH measures (F_RZ). The analyse of both pedigree and genomic based inbreeding coefficients in a large and robust population such as the PRE horse, with proven parenteral information for the last 40 years and a high degree of completeness (over 90% for the last 70 years) will allow us to understand PRE genetic variability better and the correlations between the estimations will give the data greater reliability.

RESULTS

The mean values of the pedigree-based inbreeding coefficients ranged from 0.01 (F for the last 3 generations -F3-) to 0.44 (ancestral history coefficient) and the mean values of genomic-based inbreeding coefficients varied from 0.05 (F_RZ for three generations, F_H and F_HOM) to 0.11 (F_RZ for nine generations). Significant correlations were also found between pedigree and genomic inbreeding values, which ranged between 0.58 (F3 with F_HOM) and 0.79 (F with F_RZ). In addition, the correlations between F_RZ estimated for the last 20 generations and the pedigree-based inbreeding highlight the fact that fewer generations of genomic data are required when comparing total inbreeding values, and the opposite when ancient values are calculated.

CONCLUSIONS

Ultimately, our results show that it is still useful to work with a deep and reliable pedigree in pedigree-based genetic studies with very large effective population sizes. Obtaining a satisfactory parameter will always be desirable, but the approximation obtained with a robust pedigree will allow us to work more efficiently and economically than with massive genotyping.

New cycle, same old mistakes? Overlapping vs. discrete generations in long-term recurrent selection

Background

Recurrent selection is a foundational breeding method for quantitative trait improvement. It typically features rapid breeding cycles that can lead to high rates of genetic gain. Usually, generations are discrete in recurrent selection, which means that breeding candidates are evaluated and considered for selection for only one cycle. Alternately, generations can overlap, with breeding candidates considered for selection as parents for multiple cycles. With recurrent genomic selection but not phenotypic selection, candidates can be re-evaluated by using genomic estimated breeding values without additional phenotyping of the candidates themselves. Therefore, it may be that candidates with true high breeding values discarded in one cycle due to underestimation of breeding value could be identified and selected in subsequent cycles. The consequences of allowing generations to overlap in recurrent selection are unknown. We assessed whether maintaining overlapping and discrete generations led to differences in genetic gain for phenotypic, genomic truncation, and genomic optimum contribution recurrent selection by stochastic simulation.

Results

With phenotypic selection, overlapping generations led to decreased genetic gain compared to discrete generations due to increased selection error bias. Selected individuals, which were in the upper tail of the distribution of phenotypic values, tended to also have high absolute error relative to their true breeding value compared to the overall population. Without repeated phenotyping, these individuals erroneously believed to have high value were repeatedly selected across cycles, leading to decreased genetic gain. With genomic truncation selection, overlapping and discrete generations performed similarly as updating breeding values precluded repeatedly selecting individuals with inaccurately high estimates of breeding values in subsequent cycles. Overlapping generations did not outperform discrete generations in the presence of a positive genetic trend with genomic truncation selection, as individuals from previous breeding cycles typically had truly lower breeding values than candidates from the current generation. With genomic optimum contribution selection, overlapping and discrete generations performed similarly, but overlapping generations slightly outperformed discrete generations in the long term if the targeted inbreeding rate was extremely low.

Conclusion

Maintaining discrete generations in recurrent phenotypic selection leads to increased genetic gain, especially at low heritabilities, by preventing selection error bias. With genomic truncation selection and genomic optimum contribution selection, genetic gain does not differ between discrete and overlapping generations assuming non-genetic effects are not present. Overlapping generations may increase genetic gain in the long term with very low targeted rates of inbreeding in genomic optimum contribution selection.

Supplementary information

The online version contains supplementary material available at 10.1186/s12864-022-08929-3.

Empirical comparison of genomic and phenotypic selection for resistance to Fusarium ear rot and fumonisin contamination in maize

GS and PS performed similarly in improving resistance to FER and FUM content. With cheaper and faster genotyping methods, GS has the potential to be more efficient than PS. Fusarium verticillioides is a common maize (Zea mays L.) pathogen that causes Fusarium ear rot (FER) and produces the mycotoxin fumonisin (FUM). This study empirically compared phenotypic selection (PS) and genomic selection (GS) for improving FER and FUM resistance. Three intermating generations of recurrent GS were conducted in the same time frame and from a common base population as two generations of recurrent PS. Lines sampled from each PS and GS cycle were evaluated in three North Carolina environments in 2020. We observed similar cumulative responses to GS and PS, representing decreases of about 50% of mean FER and FUM compared to the base population. The first cycle of GS was more effective than later cycles. PS and GS both achieved about 70% of predicted total gain from selection for FER, but only about 26% of predicted gains for FUM, suggesting that heritability for FUM was overestimated. We observed a 20% decrease in genetic marker variation from PS and 30% decrease from GS. Our greatest challenge was our inability to quickly obtain dense and consistent set of marker genotypes across generations of GS. Practical implementation of GS in individual small-scale breeding programs will require cheaper and faster genotyping methods, and such technological advances will present opportunities to significantly optimize selection and mating schemes for future GS efforts beyond what we were able to achieve in this study.

Assessing Population Structure and Genetic Diversity in U.S. Suffolk Sheep to Define a Framework for Genomic Selection

Long-term sustainability of breeds depends on having sufficient genetic diversity for adaptability to change, whether driven by climatic conditions or by priorities in breeding programs. Genetic diversity in Suffolk sheep in the U.S. was evaluated in four ways: 1) using genetic relationships from pedigree data [(n=64,310 animals recorded in the U.S. National Sheep Improvement Program (NSIP)]; 2) using molecular data (n=304 Suffolk genotyped with the OvineHD BeadChip); 3) comparing Australian (n=109) and Irish (n=55) Suffolk sheep to those in the U.S. using molecular data; and 4) assessing genetic relationships (connectedness) among active Suffolk flocks (n=18) in NSIP. By characterizing genetic diversity, a goal was to define the structure of a reference population for use for genomic selection strategies in this breed. Pedigree-based mean inbreeding level for the most recent year of available data was 5.5%. Ten animals defined 22.8% of the current gene pool. The effective population size (Ne) ranged from 27.5 to 244.2 based on pedigree and was 79.5 based on molecular data. Expected (HE) and observed (HO) heterozygosity were 0.317 and 0.306, respectively. Model-based population structure included 7 subpopulations. From Principal Component Analysis, countries separated into distinct populations. Within the U.S. population, flocks formed genetically disconnected clusters. A decline in genetic diversity over time was observed from both pedigree and genomic-based derived measures with evidence of population substructure as measured by FST. Using these measures of genetic diversity, a framework for establishing a genomic reference population in U.S. Suffolk sheep engaged in NSIP was proposed.

The long-term effects of genomic selection: 1. Response to selection, additive genetic variance, and genetic architecture

Background

Genomic selection has revolutionized genetic improvement in animals and plants, but little is known about its long-term effects. Here, we investigated the long-term effects of genomic selection on response to selection, genetic variance, and the genetic architecture of traits using stochastic simulations. We defined the genetic architecture as the set of causal loci underlying each trait, their allele frequencies, and their statistical additive effects. We simulated a livestock population under 50 generations of phenotypic, pedigree, or genomic selection for a single trait, controlled by either only additive, additive and dominance, or additive, dominance, and epistatic effects. The simulated epistasis was based on yeast data.

Results

Short-term response was always greatest with genomic selection, while response after 50 generations was greater with phenotypic selection than with genomic selection when epistasis was present, and was always greater than with pedigree selection. This was mainly because loss of genetic variance and of segregating loci was much greater with genomic and pedigree selection than with phenotypic selection. Compared to pedigree selection, selection response was always greater with genomic selection. Pedigree and genomic selection lost a similar amount of genetic variance after 50 generations of selection, but genomic selection maintained more segregating loci, which on average had lower minor allele frequencies than with pedigree selection. Based on this result, genomic selection is expected to better maintain genetic gain after 50 generations than pedigree selection. The amount of change in the genetic architecture of traits was considerable across generations and was similar for genomic and pedigree selection, but slightly less for phenotypic selection. Presence of epistasis resulted in smaller changes in allele frequencies and less fixation of causal loci, but resulted in substantial changes in statistical additive effects across generations.

Conclusions

Our results show that genomic selection outperforms pedigree selection in terms of long-term genetic gain, but results in a similar reduction of genetic variance. The genetic architecture of traits changed considerably across generations, especially under selection and when non-additive effects were present. In conclusion, non-additive effects had a substantial impact on the accuracy of selection and long-term response to selection, especially when selection was accurate.

Temporal and genomic analysis of additive genetic variance in breeding programmes

Genetic variance is a central parameter in quantitative genetics and breeding. Assessing changes in genetic variance over time as well as the genome is therefore of high interest. Here, we extend a previously proposed framework for temporal analysis of genetic variance using the pedigree-based model, to a new framework for temporal and genomic analysis of genetic variance using marker-based models. To this end, we describe the theory of partitioning genetic variance into genic variance and within-chromosome and between-chromosome linkage-disequilibrium, and how to estimate these variance components from a marker-based model fitted to observed phenotype and marker data. The new framework involves three steps: (i) fitting a marker-based model to data, (ii) sampling realisations of marker effects from the fitted model and for each sample calculating realisations of genetic values and (iii) calculating the variance of sampled genetic values by time and genome partitions. Analysing time partitions indicates breeding programme sustainability, while analysing genome partitions indicates contributions from chromosomes and chromosome pairs and linkage-disequilibrium. We demonstrate the framework with a simulated breeding programme involving a complex trait. Results show good concordance between simulated and estimated variances, provided that the fitted model is capturing genetic complexity of a trait. We observe a reduction of genetic variance due to selection and drift changing allele frequencies, and due to selection inducing negative linkage-disequilibrium.

Long-term impact of conventional and optimal contribution conservation methods on genetic diversity and genetic gain in local pig breeds

There are rich and vast genetic resources of indigenous pig breeds in the world. Currently, great attention is paid to either crossbreeding or conservation of these indigenous pig breeds, and insufficient attention is paid to the combination of conservation and breeding along with their long-term effects on genetic diversity. Therefore, the objective of this study is to compare the long-term effects of using conventional conservation and optimal contribution selection methods on genetic diversity and genetic gain. A total of 11 different methods including conventional conservation and optimal contribution selection methods were investigated using stochastic simulations. The long-term effects of using these methods were evaluated in terms of genetic diversity metrices such as expected heterozygosity (He) and the rate of genetic gain. The results indicated that the rates of true inbreeding in these conventional conservation methods were maintained at around 0.01. The optimal contribution selection methods based either on the pedigree (POCS) or genome (GOCS) information showed more genetic gain than conventional methods, and POCS achieved the largest genetic gain. Furthermore, the effect of using GOCS methods on most of the genetic diversity metrics was slightly better than the conventional conservation methods when the rate of true inbreeding was the same, but this also required more sires used in OCS methods. According to the rate of true inbreeding, there was no significant difference among these conventional methods. In conclusion, there is no significant difference in different ways of selecting sows on inbreeding when we use different conventional conservation methods. Compared with conventional methods, POCS method could achieve the most genetic gain. However, GOCS methods can not only achieve higher genetic gain, but also maintain a relatively high level of genetic diversity. Therefore, GOCS is a better choice if we want to combine conservation and breeding in actual production in the conservation farms.

Genomic mating in outbred species: predicting cross usefulness with additive and total genetic covariance matrices

Diverse crops are both outbred and clonally propagated. Breeders typically use truncation selection of parents and invest significant time, land, and money evaluating the progeny of crosses to find exceptional genotypes. We developed and tested genomic mate selection criteria suitable for organisms of arbitrary homozygosity level where the full-sibling progeny are of direct interest as future parents and/or cultivars. We extended cross variance and covariance variance prediction to include dominance effects and predicted the multivariate selection index genetic variance of crosses based on haplotypes of proposed parents, marker effects, and recombination frequencies. We combined the predicted mean and variance into usefulness criteria for parent and variety development. We present an empirical study of cassava (Manihot esculenta), a staple tropical root crop. We assessed the potential to predict the multivariate genetic distribution (means, variances, and trait covariances) of 462 cassava families in terms of additive and total value using cross-validation. Most variance (89%) and covariance (70%) prediction accuracy estimates were greater than zero. The usefulness of crosses was accurately predicted with good correspondence between the predicted and the actual mean performance of family members breeders selected for advancement as new parents and candidate varieties. We also used a directional dominance model to quantify significant inbreeding depression for most traits. We predicted 47,083 possible crosses of 306 parents and contrasted them to those previously tested to show how mate selection can reveal the new potential within the germplasm. We enable breeders to consider the potential of crosses to produce future parents (progeny with top breeding values) and varieties (progeny with top own performance).

Genomic selection of parents and crosses beyond the native gene pool of a breeding program

Genomic selection has become a valuable tool for selecting cultivar candidates in many plant breeding programs. Genomic selection of elite parents and crossing combinations with germplasm developed outside a breeding program has, however, hardly been explored until now. The aim of this study was to assess the potential of this method for commonly ranking and selecting elite germplasm developed within and beyond a given breeding program. A winter wheat (Triticum aestivum L.) population consisting of 611 in-house and 87 externally developed lines was used to compare training population compositions and statistical models for genomically predicting baking quality in this framework. Augmenting training populations with lines from other breeding programs had a larger influence on the prediction ability than adding in-house generated lines when aiming to commonly rank both germplasm sets. Exploiting preexisting information of secondary correlated traits resulted likewise in more accurate predictions both in empirical analyses and simulations. Genotyping germplasm developed beyond a given breeding program is moreover a convenient way to clarify its relationships with a breeder's own germplasm because pedigree information is oftentimes not available for this purpose. Genomic predictions can thus support a more informed diversity management, especially when integrating simply to phenotype correlated traits to partly circumvent resource reallocations for a costly phenotyping of germplasm from other programs.

Multivariate genomic analysis and optimal contributions selection predicts high genetic gains in cooking time, iron, zinc, and grain yield in common beans in East Africa

Common bean (Phaseolus vulgaris L.) is important in African diets for protein, iron (Fe), and zinc (Zn), but traditional cultivars have long cooking time (CKT), which increases the time, energy, and health costs of cooking. Genomic selection was used to predict genomic estimated breeding values (GEBV) for grain yield (GY), CKT, Fe, and Zn in an African bean panel of 358 genotypes in a two-stage analysis. In Stage 1, best linear unbiased estimates (BLUE) for each trait were obtained from 898 genotypes across 33 field trials in East Africa. In Stage 2, BLUE in a training population of 141 genotypes were used in a multivariate genomic analysis with genome-wide single nucleotide polymorphism data from the African bean panel. Moderate to high genomic heritability was found for GY (0.45 ± 0.10), CKT (0.50 ± 0.15), Fe (0.57 ± 0.12), and Zn (0.61 ± 0.13). There were significant favorable genetic correlations between Fe and Zn (0.91 ± 0.06), GY and Fe (0.66 ± 0.17), GY and Zn (0.44 ± 0.19), CKT and Fe (-0.57 ± 0.21), and CKT and Zn (-0.67 ± 0.20). Optimal contributions selection (OCS), based on economic index of weighted GEBV for each trait, was used to design crossing within four market groups relevant to East Africa. Progeny were predicted by OCS to increase in mean GY by 12.4%, decrease in mean CKT by 9.3%, and increase in mean Fe and Zn content by 6.9 and 4.6%, respectively, with low achieved coancestry of 0.032. Genomic selection with OCS will accelerate breeding of high-yielding, biofortified, and rapid cooking African common bean cultivars.

Strategies to Assure Optimal Trade-Offs Among Competing Objectives for the Genetic Improvement of Soybean

Plant breeding is a decision-making discipline based on understanding project objectives. Genetic improvement projects can have two competing objectives: maximize the rate of genetic improvement and minimize the loss of useful genetic variance. For commercial plant breeders, competition in the marketplace forces greater emphasis on maximizing immediate genetic improvements. In contrast, public plant breeders have an opportunity, perhaps an obligation, to place greater emphasis on minimizing the loss of useful genetic variance while realizing genetic improvements. Considerable research indicates that short-term genetic gains from genomic selection are much greater than phenotypic selection, while phenotypic selection provides better long-term genetic gains because it retains useful genetic diversity during the early cycles of selection. With limited resources, must a soybean breeder choose between the two extreme responses provided by genomic selection or phenotypic selection? Or is it possible to develop novel breeding strategies that will provide a desirable compromise between the competing objectives? To address these questions, we decomposed breeding strategies into decisions about selection methods, mating designs, and whether the breeding population should be organized as family islands. For breeding populations organized into islands, decisions about possible migration rules among family islands were included. From among 60 possible strategies, genetic improvement is maximized for the first five to 10 cycles using genomic selection and a hub network mating design, where the hub parents with the largest selection metric make large parental contributions. It also requires that the breeding populations be organized as fully connected family islands, where every island is connected to every other island, and migration rules allow the exchange of two lines among islands every other cycle of selection. If the objectives are to maximize both short-term and long-term gains, then the best compromise strategy is similar except that the mating design could be hub network, chain rule, or a multi-objective optimization method-based mating design. Weighted genomic selection applied to centralized populations also resulted in the realization of the greatest proportion of the genetic potential of the founders but required more cycles than the best compromise strategy.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

How genomic selection has increased rates of genetic gain and inbreeding in the Australian national herd, genomic information nucleus, and bulls

Genomic selection has been commonly used for selection for over a decade. In this time, the rate of genetic gain has more than doubled in some countries, while inbreeding per year has also increased. Inbreeding can result in a loss of genetic diversity, decreased long-term response to selection, reduced animal performance and ultimately, decreased farm profitability. We quantified and compared changes in genetic gain and diversity resulting from genomic selection in Australian Holstein and Jersey cattle populations. To increase the accuracy of genomic selection, Australia has had a female genomic reference population since 2013, specifically designed to be representative of commercial populations and thus including both Holstein and Jersey cows. Herds that kept excellent health and fertility data were invited to join this population and most their animals were genotyped. In both breeds, the rate of genetic gain and inbreeding was greatest in bulls, and then the female genomic reference population, and finally the wider national herd. When comparing pre- and postgenomic selection, the rates of genetic gain for the national economic index has increased by ~160% in Holstein females and ~100% in Jersey females. This has been accompanied by doubling of the rates of inbreeding in female populations, and the rate of inbreeding has increased several fold in Holstein bulls since the widespread use of genomic selection. Where cow genotype data were available to perform a more accurate genomic analysis, greater rates of pedigree and genomic inbreeding were observed, indicating actual inbreeding levels could be underestimated in the national population due to gaps in pedigrees. Based on current rates of genetic gain, the female reference population is progressing ahead of the national herd and could be used to infer and track the future inbreeding and genetic trends of the national herds.

Estimation of Population-Specific Genetic Parameters Important for Long-Term Optimum Contribution Selection-Case Study on a Dairy Istrian Sheep Breed

The Istrian sheep breed has been subjected to selection for dairy traits for more than two decades. However, a detailed study of some important population-specific parameters such as effective population size (Ne) and connectedness between flocks has never been carried out. The aim of the study was to examine the above parameters in dairy Istrian sheep subjected to a national selection program. The Ne was estimated as the mean rate of increase in coancestry, and connectedness was determined using four different statistics. The Ne was estimated at 73 animals with pedigree constraints imposed on 4 equivalent generations and 3 full generations. Analysis of ΔNe ("sliding window approach") revealed a negative ΔNe indicating a progressive loss of genetic variability (ΔNeNEG≥4 = -6.6, p < 0.01; ΔNeNFG≥3 = -4.9, p > 0.05). The overall connectedness (r¯ ~ 0.0001) was below the acceptable level for unbiased ranking of the animals belonging to different flocks (ri,j  = 0.05). OCS appears to be the best option for the long-term survival (self-sufficiency) of the breed, but genetic links between flocks need to be strengthened to allow unbiased ranking of the animals based on the estimated breeding values.

Gaining Acceptance of Novel Plant Breeding Technologies

Ensuring the sustainability of agriculture under climate change has led to a surge in alternative strategies for crop improvement. Advances in integrated crop breeding, social acceptance, and farm-level adoption are crucial to address future challenges to food security. Societal acceptance can be slow when consumers do not see the need for innovation or immediate benefits. We consider how best to address the issue of social licence and harmonised governance for novel gene technologies in plant breeding. In addition, we highlight optimised breeding strategies that will enable long-term genetic gains to be achieved. Promoted by harmonised global policy change, innovative plant breeding can realise high and sustainable productivity together with enhanced nutritional traits.

Long-term comparison between index selection and optimal independent culling in plant breeding programs with genomic prediction

In the context of genomic selection, we evaluated and compared breeding programs using either index selection or independent culling for recurrent selection of parents. We simulated a clonally propagated crop breeding program for 20 cycles using either independent culling or an economic index with two unfavourably correlated traits under selection. Cycle time from crossing to selection of parents was kept the same for both strategies. Both methods led to increasingly unfavourable genetic correlations between traits and, compared to independent culling, index selection led to larger changes in the genetic correlation between the two traits. When linkage disequilibrium was not considered, the two methods had similar losses of genetic diversity. Two independent culling approaches were evaluated, one using optimal culling levels and one using the same selection intensity for both traits. Optimal culling levels outperformed the same selection intensity even when traits had the same economic importance. Therefore, accurately estimating optimal culling levels is essential for maximizing gains when independent culling is performed. Once optimal culling levels are achieved, independent culling and index selection lead to comparable genetic gains.

Changes in Allele Frequencies When Different Genomic Coancestry Matrices Are Used for Maintaining Genetic Diversity

A main objective in conservation programs is to maintain genetic variability. This can be achieved using the Optimal Contributions (OC) method that optimizes the contributions of candidates to the next generation by minimizing the global coancestry. However, it has been argued that maintaining allele frequencies is also important. Different genomic coancestry matrices can be used on OC and the choice of the matrix will have an impact not only on the genetic variability maintained, but also on the change in allele frequencies. The objective of this study was to evaluate, through stochastic simulations, the genetic variability maintained and the trajectory of allele frequencies when using two different genomic coancestry matrices in OC to minimize the loss of diversity: (i) the matrix based on deviations of the observed number of alleles shared between two individuals from the expected numbers under Hardy–Weinberg equilibrium (θ_LH); and (ii) the matrix based on VanRaden’s genomic relationship matrix (θ_VR). The results indicate that the use of θ_LH resulted in a higher genetic variability than the use of θ_VR. However, the use of θ_VR maintained allele frequencies closer to those in the base population than the use of θ_LH.

Heterosis and Hybrid Crop Breeding: A Multidisciplinary Review

Although hybrid crop varieties are among the most popular agricultural innovations, the rationale for hybrid crop breeding is sometimes misunderstood. Hybrid breeding is slower and more resource-intensive than inbred breeding, but it allows systematic improvement of a population by recurrent selection and exploitation of heterosis simultaneously. Inbred parental lines can identically reproduce both themselves and their F1 progeny indefinitely, whereas outbred lines cannot, so uniform outbred lines must be bred indirectly through their inbred parents to harness heterosis. Heterosis is an expected consequence of whole-genome non-additive effects at the population level over evolutionary time. Understanding heterosis from the perspective of molecular genetic mechanisms alone may be elusive, because heterosis is likely an emergent property of populations. Hybrid breeding is a process of recurrent population improvement to maximize hybrid performance. Hybrid breeding is not maximization of heterosis per se, nor testing random combinations of individuals to find an exceptional hybrid, nor using heterosis in place of population improvement. Though there are methods to harness heterosis other than hybrid breeding, such as use of open-pollinated varieties or clonal propagation, they are not currently suitable for all crops or production environments. The use of genomic selection can decrease cycle time and costs in hybrid breeding, particularly by rapidly establishing heterotic pools, reducing testcrossing, and limiting the loss of genetic variance. Open questions in optimal use of genomic selection in hybrid crop breeding programs remain, such as how to choose founders of heterotic pools, the importance of dominance effects in genomic prediction, the necessary frequency of updating the training set with phenotypic information, and how to maintain genetic variance and prevent fixation of deleterious alleles.

Application of multi-trait Bayesian decision theory for parental genomic selection

In all breeding programs, the decision about which individuals to select and intermate to form the next selection cycle is crucial. The improvement of genetic stocks requires considering multiple traits simultaneously, given that economic value and net genetic merits depend on many traits; therefore, with the advance of computational and statistical tools and genomic selection (GS), researchers are focusing on multi-trait selection. Selection of the best individuals is difficult, especially in traits that are antagonistically correlated, where improvement in one trait might imply a reduction in other(s). There are approaches that facilitate multi-trait selection, and recently a Bayesian decision theory (BDT) has been proposed. Parental selection using BDT has the potential to be effective in multi-trait selection given that it summarizes all relevant quantitative genetic concepts such as heritability, response to selection and the structure of dependence between traits (correlation). In this study, we applied BDT to provide a treatment for the complexity of multi-trait parental selection using three multivariate loss functions (LF), Kullback–Leibler (KL), Energy Score, and Multivariate Asymmetric Loss (MALF), to select the best-performing parents for the next breeding cycle in two extensive real wheat data sets. Results show that the high ranking lines in genomic estimated breeding value (GEBV) for certain traits did not always have low values for the posterior expected loss (PEL). For both data sets, the KL LF gave similar importance to all traits including grain yield. In contrast, the Energy Score and MALF gave a better performance in three of four traits that were different than grain yield. The BDT approach should help breeders to decide based not only on the GEBV per se of the parent to be selected, but also on the level of uncertainty according to the Bayesian paradigm.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Keep Garfagnina alive. An integrated study on patterns of homozygosity, genomic inbreeding, admixture and breed traceability of the Italian Garfagnina goat breed

The objective of this study was to investigate the genetic diversity of the Garfagnina (GRF) goat, a breed that currently risks extinction. For this purpose, 48 goats were genotyped with the Illumina CaprineSNP50 BeadChip and analyzed together with 214 goats belonging to 9 other Italian breeds (~25 goats/breed), whose genotypes were available from the AdaptMap project [Argentata (ARG), Bionda dell'Adamello (BIO), Ciociara Grigia (CCG), Di Teramo (DIT), Garganica (GAR), Girgentana (GGT), Orobica (ORO), Valdostana (VAL) and Valpassiria (VSS)]. Comparative analyses were conducted on i) runs of homozygosity (ROH), ii) admixture ancestries and iii) the accuracy of breed traceability via discriminant analysis on principal components (DAPC) based on cross-validation. ROH analyses was used to assess the genetic diversity of GRF, while admixture and DAPC to evaluate its relationship to the other breeds. For GRF, common ROH (more than 45% in GRF samples) was detected on CHR 12 at, roughly 50.25-50.94Mbp (ARS1 assembly), which spans the CENPJ (centromere protein) and IL17D (interleukin 17D) genes. The same area of common ROH was also present in DIT, while a broader region (~49.25-51.94Mbp) was shared among the ARG, CCG, and GGT. Admixture analysis revealed a small region of common ancestry from GRF shared by BIO, VSS, ARG and CCG breeds. The DAPC model yielded 100% assignment success for GRF. Overall, our results support the identification of GRF as a distinct native Italian goat breed. This work can contribute to planning conservation programmes to save GRF from extinction and will improve the understanding of the socio-agro-economic factors related with the farming of GRF.

Selection and drift reduce genetic variation for milk yield in Manech Tête Rousse dairy sheep

Decreases in genetic variance over generations reduce future genetic gain. We studied the evolution of genetic variance in the dairy sheep breed Manech Tête Rousse, which has been selected for increasingly complex objectives, including, in this order, milk yield, milk contents, scrapie resistance, and somatic cell score. We estimated base population genetic variance and genetic variance by sex and per year of birth from 1981 to 2014. The data consisted of 1,842,295 milk yield records (from 1978 to 2017) and a pedigree including 530,572 females (96% of them with records) and 3,798 artificial insemination males. As a measure of drift, we computed average relationships for each cohort from which we derived expected reduction of variance due to increased relationships. The difference between observed and expected reductions in genetic variances is the reduction in genetic variance due to selection. Average relationships increased steadily but slowly in both sexes. For females, genetic variance reduced with time until a plateau was reached at around 90% of the initial genetic variance. The reduction due to relationships (roughly 3% cumulated in 30 yr) was smaller than that due to selection (roughly 10% across the last years). A smaller loss due to selection was seen in recent years, possibly due to a change in selection objectives. These results agree well with theoretical expectations. The pattern of the evolution of genetic variance in males was similar to that for females but with a stronger reduction because of strong selection of AI males at birth. We conclude that the reductions in genetic variation due to selection and drift agree with expectations, and none of the reductions are very strong in this population because of control of inbreeding and smooth changes in selection objectives over time.

A chickpea genetic variation map based on the sequencing of 3,366 genomes

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources¹. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level². Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

Optimal Contribution Selection in Highly Fecund Species With Overlapping Generations

Optimal contributions approaches to parental selection in closed breeding populations aim to maximize genetic gains, while restraining long-term inbreeding. The adoption of optimal contribution selection (OCS) in highly fecund outcrossing species presents a number of challenges not applicable to species of low fecundity (e.g., livestock) for which they were developed. This is particularly true if overlapping-generations or rolling-front breeding strategies are applied, in which case the number of individuals per family in juvenile (i.e., sexually immature) age groups is not necessarily known but is likely to be large. In these circumstances, conventional OCS procedures must be modified or a large number of dummy individuals defined, making computations onerous. Here, an approach to OCS is presented that involves the use of “between-family relationship matrices” instead of “between-individual relationship matrices.” The method is applicable to breeding programs involving highly fecund outcrossing species with overlapping generations, including circumstances where the number of juveniles per family is unknown but large.

RG, Bauck S. Genomic mating as sustainable breeding for Chinese indigenous Ningxiang pigs

An important economic reason for the loss of local breeds is that they tend to be less productive, and hence having less market value than commercial breeds. Nevertheless, local breeds often have irreplaceable values, genetically and sociologically. In the breeding programs with local breeds, it is crucial to balance the selection for genetic gain and the maintaining of genetic diversity. These two objectives are often conflicting, and finding the optimal point of the trade-off has been a challenge for breeders. Genomic selection (GS) provides a revolutionary tool for the genetic improvement of farm animals. At the same time, it can increase inbreeding and produce a more rapid depletion of genetic variability of the selected traits in future generations. Optimum-contribution selection (OCS) represents an approach to maximize genetic gain while constraining inbreeding within a targeted range. In the present study, 515 Ningxiang pigs were genotyped with the Illumina Porcine SNP60 array or the GeneSeek Genomic Profiler Porcine 50K array. The Ningxiang pigs were found to be highly inbred at the genomic level. Average locus-wise inbreeding coefficients were 0.41 and 0.37 for the two SNP arrays used, whereas genomic inbreeding coefficients based on runs of homozygosity were 0.24 and 0.25, respectively. Simulated phenotypic data were used to assess the utility of genomic OCS (GOCS) in comparison with GS without inbreeding control. GOCS was conducted under two scenarios, selecting sires only (GOCS_S) or selecting sires and dams (GOCS_SD), while kinships were constrained on selected parents. The genetic gain for average daily body weight gain (ADG) per generation was between 18.99 and 20.55 g with GOCS_S, and between 23.20 and 28.92 with GOCS_SD, and it varied from 25.38 to 48.38 g under GS without controlling inbreeding. While the rate of genetic gain per generation obtained using GS was substantially larger than that obtained by the two scenarios of genomic OCS in the beginning generations of selection, the difference in the genetic gain of ADG between GS and GOCS reduced quickly in latter generations. At generation ten, the difference in the realized rates of genetic gain between GS and GOCS_SD diminished and ended up with even a slightly higher genetic gain with GOCS_SD, due to the rapid loss of genetic variance with GS and fixation of causative genes. The rate of inbreeding was mostly maintained below 5% per generation with genomic OCS, whereas it increased to between 10.5% and 15.3% per generation with GS. Therefore, genomic OCS appears to be a sustainable strategy for the genetic improvement of local breeds such as Ningxiang pigs, but keeping mind that a variety of GOCS methods exist and the optimal forms remain to be exploited further.

Management of Genetic Diversity in the Era of Genomics

Management of genetic diversity aims to (i) maintain heterozygosity, which ameliorates inbreeding depression and loss of genetic variation at loci that may become of importance in the future; and (ii) avoid genetic drift, which prevents deleterious recessives (e.g., rare disease alleles) from drifting to high frequency, and prevents random drift of (functional) traits. In the genomics era, genomics data allow for many alternative measures of inbreeding and genomic relationships. Genomic relationships/inbreeding can be classified into (i) homozygosity/heterozygosity based (e.g., molecular kinship matrix); (ii) genetic drift-based, i.e., changes of allele frequencies; or (iii) IBD-based, i.e., SNPs are used in linkage analyses to identify IBD segments. Here, alternative measures of inbreeding/relationship were used to manage genetic diversity in genomic optimal contribution (GOC) selection schemes. Contrary to classic inbreeding theory, it was found that drift and homozygosity-based inbreeding could differ substantially in GOC schemes unless diversity management was based upon IBD. When using a homozygosity-based measure of relationship, the inbreeding management resulted in allele frequency changes toward 0.5 giving a low rate of increase in homozygosity for the panel used for management, but not for unmanaged neutral loci, at the expense of a high genetic drift. When genomic relationship matrices were based on drift, following VanRaden and as in GCTA, drift was low at the expense of a high rate of increase in homozygosity. The use of IBD-based relationship matrices for inbreeding management limited both drift and the homozygosity-based rate of inbreeding to their target values. Genetic improvement per percent of inbreeding was highest when GOC used IBD-based relationships irrespective of the inbreeding measure used. Genomic relationships based on runs of homozygosity resulted in very high initial improvement per percent of inbreeding, but also in substantial discrepancies between drift and homozygosity-based rates of inbreeding, and resulted in a drift that exceeded its target value. The discrepancy between drift and homozygosity-based rates of inbreeding was caused by a covariance between initial allele frequency and the subsequent change in frequency, which becomes stronger when using data from whole genome sequence.

Selection criteria for high-yielding and early-flowering bread wheat hybrids under heat stress

Hybrid performance during wheat breeding can be improved by analyzing genetic distance (GD) among wheat genotypes and determining its correlation with heterosis. This study evaluated the GD between 16 wheat genotypes by using 60 simple sequence repeat (SSR) markers to classify them according to their relationships and select those with greater genetic diversity, evaluate the correlation of the SSR marker distance with heterotic performance and specific combining ability (SCA) for heat stress tolerance, and identify traits that most influence grain yield (GY). Eight parental genotypes with greater genetic diversity and their 28 F1 hybrids generated using diallel crossing were evaluated for 12 measured traits in two seasons. The GD varied from 0.235 to 0.911 across the 16 genotypes. Cluster analysis based on the GD estimated using SSRs classified the genotypes into three major groups and six sub-groups, almost consistent with the results of principal coordinate analysis. The combined data indicated that five hybrids showed 20% greater yield than mid-parent or better-parent. Two hybrids (P2 × P4) and (P2 × P5), which showed the highest performance of days to heading (DH), grain filling duration (GFD), and GY, and had large genetic diversity among themselves (0.883 and 0.911, respectively), were deemed as promising heat-tolerant hybrids. They showed the best mid-parent heterosis and better-parent heterosis (BPH) for DH (-11.57 and -7.65%; -13.39 and -8.36%, respectively), GFD (12.74 and 12.17%; 12.09 and 10.59%, respectively), and GY (36.04 and 20.04%; 44.06 and 37.73%, respectively). Correlation between GD and each of BPH and SCA effects based on SSR markers was significantly positive for GFD, hundred kernel weight, number of kernels per spike, harvest index, GY, and grain filling rate and was significantly negative for DH. These correlations indicate that the performance of wheat hybrids with high GY and earliness could be predicted by determining the GD of the parents by using SSR markers. Multivariate analysis (stepwise regression and path coefficient) suggested that GFD, hundred kernel weight, days to maturity, and number of kernels per spike had the highest influence on GY.

Genetic Structure Analysis of the Pura Raza Español Horse Population through Partial Inbreeding Coefficient Estimation

The aim of this work was to analyze genetic parameters such as the inbreeding coefficient (F), relatedness coefficient (AR) and partial inbreeding coefficient (Fij) of the whole PRE population, and the ancestors which account for 50% of the total genetic variability of the current population, from genealogical information. The average F of the whole PRE population (328,706 animals) has decreased from 8.45% to 7.51% in the least 20 years. The Fij was estimated for the whole PRE population, resulting in a database of 58,772,533 records containing one record for each Fij that each animal receives from a certain common ancestor (CA). A total of 10,244 CAs contributed to the Fij with an average of 5370 descendants, with each descendant having an average of 170 CAs. Over the generations, the number of CAs has increased, while the proportion of Fij by each one has decreased. In addition, the contributions of the more influential ancestors have changed. The increased census, the limited use of artificial insemination and our increased knowledge about inbreeding depression and the animals' breeding values allow breeders to select horses more for their functionality and conformation than for their pedigree reputation, which is the cause of all these changes.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

Origin Specific Genomic Selection: A Simple Process To Optimize the Favorable Contribution of Parents to Progeny

Modern crop breeding is in constant demand for new genetic diversity as part of the arms race with genetic gain. The elite gene pool has limited genetic variation and breeders are trying to introduce novelty from unadapted germplasm, landraces and wild relatives. For polygenic traits, currently available approaches to introgression are not ideal, as there is a demonstrable bias against exotic alleles during selection. Here, we propose a partitioned form of genomic selection, called Origin Specific Genomic Selection (OSGS), where we identify and target selection on favorable exotic alleles. Briefly, within a population derived from a bi-parental cross, we isolate alleles originating from the elite and exotic parents, which then allows us to separate out the predicted marker effects based on the allele origins. We validated the usefulness of OSGS using two nested association mapping (NAM) datasets: barley NAM (elite-exotic) and maize NAM (elite-elite), as well as by computer simulation. Our results suggest that OSGS works well in its goal to increase the contribution of favorable exotic alleles in bi-parental crosses, and it is possible to extend the approach to broader multi-parental populations.