FieldSimR: an R package for simulating plot data in multi-environment field trials.

This paper presents a general framework for simulating plot data in multi-environment field trials with one or more traits. The framework is embedded within the R package FieldSimR, whose core function generates plot errors that capture global field trend, local plot variation, and extraneous variation at a user-defined ratio. FieldSimR's capacity to simulate realistic plot data makes it a flexible and powerful tool for a wide range of improvement processes in plant breeding, such as the optimisation of experimental designs and statistical analyses of multi-environment field trials. FieldSimR provides crucial functionality that is currently missing in other software for simulating plant breeding programmes and is available on CRAN. The paper includes an example simulation of field trials that evaluate 100 maize hybrids for two traits in three environments. To demonstrate FieldSimR's value as an optimisation tool, the simulated data set is then used to compare several popular spatial models for their ability to accurately predict the hybrids' genetic values and reliably estimate the variance parameters of interest. FieldSimR has broader applications to simulating data in other agricultural trials, such as glasshouse experiments.

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.

Genomic selection strategies for clonally propagated crops.

KEY MESSAGE: 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.

Genetic architecture and major genes for backfat thickness in pig lines of diverse genetic backgrounds.

BACKGROUND: Backfat thickness is an important carcass composition trait for pork production and is commonly included in swine breeding programmes. In this paper, we report the results of a large genome-wide association study for backfat thickness using data from eight lines of diverse genetic backgrounds. METHODS: Data comprised 275,590 pigs from eight lines with diverse genetic backgrounds (breeds included Large White, Landrace, Pietrain, Hampshire, Duroc, and synthetic lines) genotyped and imputed for 71,324 single-nucleotide polymorphisms (SNPs). For each line, we estimated SNP associations using a univariate linear mixed model that accounted for genomic relationships. SNPs with significant associations were identified using a threshold of p < 10-6 and used to define genomic regions of interest. The proportion of genetic variance explained by a genomic region was estimated using a ridge regression model. RESULTS: We found significant associations with backfat thickness for 264 SNPs across 27 genomic regions. Six genomic regions were detected in three or more lines. The average estimate of the SNP-based heritability was 0.48, with estimates by line ranging from 0.30 to 0.58. The genomic regions jointly explained from 3.2 to 19.5% of the additive genetic variance of backfat thickness within a line. Individual genomic regions explained up to 8.0% of the additive genetic variance of backfat thickness within a line. Some of these 27 genomic regions also explained up to 1.6% of the additive genetic variance in lines for which the genomic region was not statistically significant. We identified 64 candidate genes with annotated functions that can be related to fat metabolism, including well-studied genes such as MC4R, IGF2, and LEPR, and more novel candidate genes such as DHCR7, FGF23, MEDAG, DGKI, and PTN. CONCLUSIONS: Our results confirm the polygenic architecture of backfat thickness and the role of genes involved in energy homeostasis, adipogenesis, fatty acid metabolism, and insulin signalling pathways for fat deposition in pigs. The results also suggest that several less well-understood metabolic pathways contribute to backfat development, such as those of phosphate, calcium, and vitamin D homeostasis.

Regional association analysis coupled with transcriptome analyses reveal candidate genes affecting seed oil accumulation in Brassica napus.

KEY MESSAGE: Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions in combination with co-expression analysis reveal candidate genes affecting oil accumulation in Brassica napus. One of the breeding goals in rapeseed production is to enhance the seed oil content to cater to the increased demand for vegetable oils due to a growing global population. To investigate the genetic basis of variation in seed oil content, we used 60 K Brassica Infinium SNP array along with phenotype data of 203 Chinese semi-winter rapeseed accessions to perform a genome-wide analysis of haplotype blocks associated with the oil content. Nine haplotype regions harbouring lipid synthesis/transport-, carbohydrate metabolism- and photosynthesis-related genes were identified as significantly associated with the oil content and were mapped to chromosomes A02, A04, A05, A07, C03, C04, C05, C08 and C09, respectively. Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions combined with transcriptome datasets from 13 accessions was further performed on these nine haplotype regions. This revealed natural variation in the BnTGD3-A02 and BnSSE1-A05 gene regions correlated with the phenotypic variation of the oil content within the A02 and A04 chromosome haplotype regions, respectively. Moreover, co-expression network analysis revealed that BnTGD3-A02 and BnSSE1-A05 were directly linked with fatty acid beta-oxidation-related gene BnKAT2-C04, thus forming a molecular network involved in the potential regulation of seed oil accumulation. The results of this study could be used to combine favourable haplotype alleles for further improvement of the seed oil content in rapeseed.

Multi-omics-based prediction of hybrid performance in canola.

KEY MESSAGE: Complementing or replacing genetic markers with transcriptomic data and use of reproducing kernel Hilbert space regression based on Gaussian kernels increases hybrid prediction accuracies for complex agronomic traits in canola. In plant breeding, hybrids gained particular importance due to heterosis, the superior performance of offspring compared to their inbred parents. Since the development of new top performing hybrids requires labour-intensive and costly breeding programmes, including testing of large numbers of experimental hybrids, the prediction of hybrid performance is of utmost interest to plant breeders. In this study, we tested the effectiveness of hybrid prediction models in spring-type oilseed rape (Brassica napus L./canola) employing different omics profiles, individually and in combination. To this end, a population of 950 F1 hybrids was evaluated for seed yield and six other agronomically relevant traits in commercial field trials at several locations throughout Europe. A subset of these hybrids was also evaluated in a climatized glasshouse regarding early biomass production. For each of the 477 parental rapeseed lines, 13,201 single nucleotide polymorphisms (SNPs), 154 primary metabolites, and 19,479 transcripts were determined and used as predictive variables. Both, SNP markers and transcripts, effectively predict hybrid performance using (genomic) best linear unbiased prediction models (gBLUP). Compared to models using pure genetic markers, models incorporating transcriptome data resulted in significantly higher prediction accuracies for five out of seven agronomic traits, indicating that transcripts carry important information beyond genomic data. Notably, reproducing kernel Hilbert space regression based on Gaussian kernels significantly exceeded the predictive abilities of gBLUP models for six of the seven agronomic traits, demonstrating its potential for implementation in future canola breeding programmes.

Modeling Illustrates That Genomic Selection Provides New Opportunities for Intercrop Breeding.

Intercrop breeding programs using genomic selection can produce faster genetic gain than intercrop breeding programs using phenotypic selection. Intercropping is an agricultural practice in which two or more component crops are grown together. It can lead to enhanced soil structure and fertility, improved weed suppression, and better control of pests and diseases. Especially in subsistence agriculture, intercropping has great potential to optimize farming and increase profitability. However, breeding for intercrop varieties is complex as it requires simultaneous improvement of two or more component crops that combine well in the field. We hypothesize that genomic selection can significantly simplify and accelerate the process of breeding crops for intercropping. Therefore, we used stochastic simulation to compare four different intercrop breeding programs implementing genomic selection and an intercrop breeding program entirely based on phenotypic selection. We assumed three different levels of genetic correlation between monocrop grain yield and intercrop grain yield to investigate how the different breeding strategies are impacted by this factor. We found that all four simulated breeding programs using genomic selection produced significantly more intercrop genetic gain than the phenotypic selection program regardless of the genetic correlation with monocrop yield. We suggest a genomic selection strategy which combines monocrop and intercrop trait information to predict general intercropping ability to increase selection accuracy in the early stages of a breeding program and to minimize the generation interval.

GWAS and co-expression network combination uncovers multigenes with close linkage effects on the oleic acid content accumulation in Brassica napus.

BACKGROUND: Strong artificial and natural selection causes the formation of highly conserved haplotypes that harbor agronomically important genes. GWAS combination with haplotype analysis has evolved as an effective method to dissect the genetic architecture of complex traits in crop species. RESULTS: We used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with oleic acid (C18:1) in rapeseed. Six haplotype regions were identified as significantly associated with oleic acid (C18:1) that mapped to chromosomes A02, A07, A08, C01, C02, and C03. Additionally, whole-genome sequencing of 50 rapeseed accessions revealed three genes (BnmtACP2-A02, BnABCI13-A02 and BnECI1-A02) in the A02 chromosome haplotype region and two genes (BnFAD8-C02 and BnSDP1-C02) in the C02 chromosome haplotype region that were closely linked to oleic acid content phenotypic variation. Moreover, the co-expression network analysis uncovered candidate genes from these two different haplotype regions with potential regulatory interrelationships with oleic acid content accumulation. CONCLUSIONS: Our results suggest that several candidate genes are closely linked, which provides us with an opportunity to develop functional haplotype markers for the improvement of the oleic acid content in rapeseed.

How Population Structure Impacts Genomic Selection Accuracy in Cross-Validation: Implications for Practical Breeding.

Over the last two decades, the application of genomic selection has been extensively studied in various crop species, and it has become a common practice to report prediction accuracies using cross validation. However, genomic prediction accuracies obtained from random cross validation can be strongly inflated due to population or family structure, a characteristic shared by many breeding populations. An understanding of the effect of population and family structure on prediction accuracy is essential for the successful application of genomic selection in plant breeding programs. The objective of this study was to make this effect and its implications for practical breeding programs comprehensible for breeders and scientists with a limited background in quantitative genetics and genomic selection theory. We, therefore, compared genomic prediction accuracies obtained from different random cross validation approaches and within-family prediction in three different prediction scenarios. We used a highly structured population of 940 Brassica napus hybrids coming from 46 testcross families and two subpopulations. Our demonstrations show how genomic prediction accuracies obtained from among-family predictions in random cross validation and within-family predictions capture different measures of prediction accuracy. While among-family prediction accuracy measures prediction accuracy of both the parent average component and the Mendelian sampling term, within-family prediction only measures how accurately the Mendelian sampling term can be predicted. With this paper we aim to foster a critical approach to different measures of genomic prediction accuracy and a careful analysis of values observed in genomic selection experiments and reported in literature.

Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola.

A major challenge of plant biology is to unravel the genetic basis of complex traits. We took advantage of recent technical advances in high-throughput phenotyping in conjunction with genome-wide association studies to elucidate genotype-phenotype relationships at high temporal resolution. A diverse Brassica napus population from a commercial breeding programme was analysed by automated non-invasive phenotyping. Time-resolved data for early growth-related traits, including estimated biovolume, projected leaf area, early plant height and colour uniformity, were established and complemented by fresh and dry weight biomass. Genome-wide SNP array data provided the framework for genome-wide association analyses. Using time point data and relative growth rates, multiple robust main effect marker-trait associations for biomass and related traits were detected. Candidate genes involved in meristem development, cell wall modification and transcriptional regulation were detected. Our results demonstrate that early plant growth is a highly complex trait governed by several medium and many small effect loci, most of which act only during short phases. These observations highlight the importance of taking the temporal patterns of QTL/allele actions into account and emphasize the need for detailed time-resolved analyses to effectively unravel the complex and stage-specific contributions of genes affecting growth processes that operate at different developmental phases.

Genetic insights into underground responses to Fusarium graminearum infection in wheat.

The ongoing global intensification of wheat production will likely be accompanied by a rising pressure of Fusarium diseases. While utmost attention was given to Fusarium head blight (FHB) belowground plant infections of the pathogen have largely been ignored. The current knowledge about the impact of soil borne Fusarium infection on plant performance and the underlying genetic mechanisms for resistance remain very limited. Here, we present the first large-scale investigation of Fusarium root rot (FRR) resistance using a diverse panel of 215 international wheat lines. We obtained data for a total of 21 resistance-related traits, including large-scale Real-time PCR experiments to quantify fungal spread. Association mapping and subsequent haplotype analyses discovered a number of highly conserved genomic regions associated with resistance, and revealed a significant effect of allele stacking on the stembase discoloration. Resistance alleles were accumulated in European winter wheat germplasm, implying indirect prior selection for improved FRR resistance in elite breeding programs. Our results give first insights into the genetic basis of FRR resistance in wheat and demonstrate how molecular parameters can successfully be explored in genomic prediction. Ongoing work will help to further improve our understanding of the complex interactions of genetic factors influencing FRR resistance.

Effective Genomic Selection in a Narrow-Genepool Crop with Low-Density Markers: Asian Rapeseed as an Example.

Genomic selection (GS) has revolutionized breeding for quantitative traits in plants, offering potential to optimize resource allocation in breeding programs and increase genetic gain per unit of time. Modern high-density single nucleotide polymorphism (SNP) arrays comprising up to several hundred thousand markers provide a user-friendly technology to characterize the genetic constitution of whole populations and for implementing GS in breeding programs. However, GS does not build upon detailed genotype profiling facilitated by maximum marker density. With extensive genome-wide linkage disequilibrium (LD) being a common characteristic of breeding pools, fewer representative markers from available high-density genotyping platforms could be sufficient to capture the association between a genomic region and a phenotypic trait. To examine the effects of reduced marker density on genomic prediction accuracy, we collected data on three traits across 2 yr in a panel of 203 homozygous Chinese semiwinter rapeseed ( L.) inbred lines, broadly encompassing allelic variability in the Asian genepool. We investigated two approaches to selecting subsets of markers: a trait-dependent strategy based on genome-wide association study (GWAS) significance thresholds and a trait-independent method to detect representative tag SNPs. Prediction accuracies were evaluated using cross-validation with ridge-regression best linear unbiased predictions (rrBLUP). With semiwinter rapeseed as a model species, we demonstrate that low-density marker sets comprising a few hundred to a few thousand markers enable high prediction accuracies in breeding populations with strong LD comparable to those achieved with high-density arrays. Our results are valuable for facilitating routine application of cost-efficient GS in breeding programs.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture.

KEY MESSAGE: Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F1-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Exploring and Harnessing Haplotype Diversity to Improve Yield Stability in Crops.

In order to meet future food, feed, fiber, and bioenergy demands, global yields of all major crops need to be increased significantly. At the same time, the increasing frequency of extreme weather events such as heat and drought necessitates improvements in the environmental resilience of modern crop cultivars. Achieving sustainably increase yields implies rapid improvement of quantitative traits with a very complex genetic architecture and strong environmental interaction. Latest advances in genome analysis technologies today provide molecular information at an ultrahigh resolution, revolutionizing crop genomic research, and paving the way for advanced quantitative genetic approaches. These include highly detailed assessment of population structure and genotypic diversity, facilitating the identification of selective sweeps and signatures of directional selection, dissection of genetic variants that underlie important agronomic traits, and genomic selection (GS) strategies that not only consider major-effect genes. Single-nucleotide polymorphism (SNP) markers today represent the genotyping system of choice for crop genetic studies because they occur abundantly in plant genomes and are easy to detect. SNPs are typically biallelic, however, hence their information content compared to multiallelic markers is low, limiting the resolution at which SNP-trait relationships can be delineated. An efficient way to overcome this limitation is to construct haplotypes based on linkage disequilibrium, one of the most important features influencing genetic analyses of crop genomes. Here, we give an overview of the latest advances in genomics-based haplotype analyses in crops, highlighting their importance in the context of polyploidy and genome evolution, linkage drag, and co-selection. We provide examples of how haplotype analyses can complement well-established quantitative genetics frameworks, such as quantitative trait analysis and GS, ultimately providing an effective tool to equip modern crops with environment-tailored characteristics.