Exploitation of data from breeding programs supports rapid implementation of genomic selection for key agronomic traits in perennial ryegrass

Genetic architecture of inter-specific and -generic grass hybrids by network analysis on multi-omics data

BACKGROUND

Understanding the mechanisms underlining forage production and its biomass nutritive quality at the omics level is crucial for boosting the output of high-quality dry matter per unit of land. Despite the advent of multiple omics integration for the study of biological systems in major crops, investigations on forage species are still scarce.

RESULTS

Our results identified substantial changes in gene co-expression and metabolite-metabolite network topologies as a result of genetic perturbation by hybridizing L. perenne with another species within the genus (L. multiflorum) relative to across genera (F. pratensis). However, conserved hub genes and hub metabolomic features were detected between pedigree classes, some of which were highly heritable and displayed one or more significant edges with agronomic traits in a weighted omics-phenotype network. In spite of tagging relevant biological molecules as, for example, the light-induced rice 1 (LIR1), hub features were not necessarily better explanatory variables for omics-assisted prediction than features stochastically sampled and all available regressors.

CONCLUSIONS

The utilization of computational techniques for the reconstruction of co-expression networks facilitates the identification of key omic features that serve as central nodes and demonstrate correlation with the manifestation of observed traits. Our results also indicate a robust association between early multi-omic traits measured in a greenhouse setting and phenotypic traits evaluated under field conditions.

Developing an integrated genomic selection approach beyond biomass for varietal protection and nutritive traits in perennial ryegrass (Lolium perenne L.)

Breeding target traits can be broadened to include nutritive value and plant breeder's rights traits in perennial ryegrass by using in-field regression-based spectroscopy phenotyping and genomic selection. Perennial ryegrass breeding has focused on biomass yield, but expansion into a broader set of traits is needed to benefit livestock industries whilst also providing support for intellectual property protection of cultivars. Numerous breeding objectives can be targeted simultaneously with the development of sensor-based phenomics and genomic selection (GS). Of particular interest are nutritive value (NV), which has been difficult and expensive to measure using traditional phenotyping methods, resulting in limited genetic improvement to date, and traits required to obtain varietal protection, known as plant breeder's rights (PBR) traits. In order to assess phenotyping requirements for NV improvement and potential for genetic improvement, in-field reflectance-based spectroscopy was assessed and GS evaluated in a single population for three key NV traits, captured across four timepoints. Using three prediction approaches, the possibility of targeting PBR traits using GS was evaluated for five traits recorded across three years of a breeding program. Prediction accuracy was generally low to moderate for NV traits and moderate to high for PBR traits, with heritability highly correlated with GS accuracy. NV did not show significant or consistent correlation between timepoints highlighting the need to incorporate seasonal NV into selection indexes and the value of being able to regularly monitor NV across seasons. This study has demonstrated the ability to implement GS for both NV and PBR traits in perennial ryegrass, facilitating the expansion of ryegrass breeding targets to agronomically relevant traits while ensuring necessary varietal protection is achieved.

Within and combined season prediction models for perennial ryegrass biomass yield using ground- and air-based sensor data

Across-season biomass assessment is crucial in the cultivar selection process to accurately evaluate the yield performance of lines under different growing conditions. However, it has been difficult to have an accurate, reliable, and repeated fresh biomass (FM) estimation of large populations of plants in the field without destructive harvesting, which incurs significant labor and operation costs. Sensor-based phenotyping platforms have advanced in the data collection of structural and vegetative information of plants, but the developed prediction models are still limited by low correlations at different growth stages and seasons. In this study, our objective was to develop and validate the global prediction models for across-season harvested fresh biomass (FM) yield based on the ground- and air-based sensor data including ground-based LiDAR, ground-based ultrasonic, and air-based multispectral camera to extract LiDAR plant volume (LV), LiDAR point density (LV_Den), height, and Normalized Difference Vegetative Index (NDVI). The study was conducted in a row-plot field trial with 480 rows (3 rows in a plot per cultivar) throughout the whole 2020 growing season up to the reproductive stage. We evaluated the performance of each plant parameter, their relationship, and the best subset prediction models using statistical stepwise selection at the row and plot levels through the seasonal and combined seasonal datasets. The best performing model: F M ~ L V ∗ L V _ D e n ∗ N D V I had a determination of coefficient R ² of at least 0.9 in vegetative stages and 0.8 in the reproductive stage. Similar results can be achieved in a simpler model with just two LiDAR variables- F M ~ L V ∗ L V _ D e n . In addition, LV and LV_Den showed a robust correlation with FM on their own over seasons and growth stages, while NDVI only performed well in some seasons. The simpler model based on only LiDAR data can be widely applied over season without the need of additional sensor data and may thus make the in-field across-season biomass assessment more feasible and practical for fast and cost-effective development of higher biomass yield cultivars.

Above-Ground Biomass Estimation in Oats Using UAV Remote Sensing and Machine Learning

Current strategies for phenotyping above-ground biomass in field breeding nurseries demand significant investment in both time and labor. Unmanned aerial vehicles (UAV) can be used to derive vegetation indices (VIs) with high throughput and could provide an efficient way to predict forage yield with high accuracy. The main objective of the study is to investigate the potential of UAV-based multispectral data and machine learning approaches in the estimation of oat biomass. UAV equipped with a multispectral sensor was flown over three experimental oat fields in Volga, South Shore, and Beresford, South Dakota, USA, throughout the pre- and post-heading growth phases of oats in 2019. A variety of vegetation indices (VIs) derived from UAV-based multispectral imagery were employed to build oat biomass estimation models using four machine-learning algorithms: partial least squares (PLS), support vector machine (SVM), Artificial neural network (ANN), and random forest (RF). The results showed that several VIs derived from the UAV collected images were significantly positively correlated with dry biomass for Volga and Beresford (r = 0.2-0.65), however, in South Shore, VIs were either not significantly or weakly correlated with biomass. For Beresford, approximately 70% of the variance was explained by PLS, RF, and SVM validation models using data collected during the post-heading phase. Likewise for Volga, validation models had lower coefficient of determination (R2 = 0.20-0.25) and higher error (RMSE = 700-800 kg/ha) than training models (R2 = 0.50-0.60; RMSE = 500-690 kg/ha). In South Shore, validation models were only able to explain approx. 15-20% of the variation in biomass, which is possibly due to the insignificant correlation values between VIs and biomass. Overall, this study indicates that airborne remote sensing with machine learning has potential for above-ground biomass estimation in oat breeding nurseries. The main limitation was inconsistent accuracy in model prediction across locations. Multiple-year spectral data, along with the inclusion of textural features like crop surface model (CSM) derived height and volumetric indicators, should be considered in future studies while estimating biophysical parameters like biomass.

Genomic Prediction of Complex Traits in Forage Plants Species: Perennial Grasses Case

The majority of forage grass species are obligate outbreeders. Their breeding classically consists of an initial selection on spaced plants for highly heritable traits such as disease resistances and heading date, followed by familial selection on swards for forage yield and quality traits. The high level of diversity and heterozygosity, and associated decay of linkage disequilibrium (LD) over very short genomic distances, has hampered the implementation of genomic selection (GS) in these species. However, next generation sequencing technologies in combination with the development of genomic resources have recently facilitated implementation of GS in forage grass species such as perennial ryegrass (Lolium perenne L.), switchgrass (Panicum virgatum L.), and timothy (Phleum pratense L.). Experimental work and simulations have shown that GS can increase significantly the genetic gain per unit of time for traits with different levels of heritability. The main reasons are (1) the possibility to select single plants based on their genomic estimated breeding values (GEBV) for traits measured at sward level, (2) a reduction in the duration of selection cycles, and less importantly (3) an increase in the selection intensity associated with an increase in the genetic variance used for selection. Nevertheless, several factors should be taken into account for the successful implementation of GS in forage grasses. For example, it has been shown that the level of relatedness between the training and the selection population is particularly critical when working with highly structured meta-populations consisting of several genetic groups. A sufficient number of markers should be used to estimate properly the kinship between individuals and to reflect the variability of major QTLs. It is also important that the prediction models are trained for relevant environments when dealing with traits with high genotype × environment interaction (G × E). Finally, in these outbreeding species, measures to reduce inbreeding should be used to counterbalance the high selection intensity that can be achieved in GS.

Genomic prediction in family bulks using different traits and cross-validations in pine

Genomic prediction integrates statistical, genomic, and computational tools to improve the estimation of breeding values and increase genetic gain. Due to the broad diversity in mating systems, breeding schemes, propagation methods, and unit of selection, no universal genomic prediction approach can be applied in all crops. In a genome-wide family prediction (GWFP) approach, the family is the basic unit of selection. We tested GWFP in two loblolly pine (Pinus taeda L.) datasets: a breeding population composed of 63 full-sib families (5–20 individuals per family), and a simulated population with the same pedigree structure. In both populations, phenotypic and genomic data was pooled at the family level in silico. Marker effects were estimated to compute genomic estimated breeding values (GEBV) at the individual and family (GWFP) levels. Less than six individuals per family produced inaccurate estimates of family phenotypic performance and allele frequency. Tested across different scenarios, GWFP predictive ability was higher than those for GEBV in both populations. Validation sets composed of families with similar phenotypic mean and variance as the training population yielded predictions consistently higher and more accurate than other validation sets. Results revealed potential for applying GWFP in breeding programs whose selection unit are family, and for systems where family can serve as training sets. The GWFP approach is well suited for crops that are routinely genotyped and phenotyped at the plot-level, but it can be extended to other breeding programs. Higher predictive ability obtained with GWFP would motivate the application of genomic prediction in these situations.

Calibration and validation of predicted genomic breeding values in an advanced cycle maize population

Model training on data from all selection cycles yielded the highest prediction accuracy by attenuating specific effects of individual cycles. Expected reliability was a robust predictor of accuracies obtained with different calibration sets. The transition from phenotypic to genome-based selection requires a profound understanding of factors that determine genomic prediction accuracy. We analysed experimental data from a commercial maize breeding programme to investigate if genomic measures can assist in identifying optimal calibration sets for model training. The data set consisted of six contiguous selection cycles comprising testcrosses of 5968 doubled haploid lines genotyped with a minimum of 12,000 SNP markers. We evaluated genomic prediction accuracies in two independent prediction sets in combination with calibration sets differing in sample size and genomic measures (effective sample size, average maximum kinship, expected reliability, number of common polymorphic SNPs and linkage phase similarity). Our results indicate that across selection cycles prediction accuracies were as high as 0.57 for grain dry matter yield and 0.76 for grain dry matter content. Including data from all selection cycles in model training yielded the best results because interactions between calibration and prediction sets as well as the effects of different testers and specific years were attenuated. Among genomic measures, the expected reliability of genomic breeding values was the best predictor of empirical accuracies obtained with different calibration sets. For grain yield, a large difference between expected and empirical reliability was observed in one prediction set. We propose to use this difference as guidance for determining the weight phenotypic data of a given selection cycle should receive in model retraining and for selection when both genomic breeding values and phenotypes are available.

Why and How to Switch to Genomic Selection: Lessons From Plant and Animal Breeding Experience

The present study is a transversal analysis of the interest in genomic selection for plant and animal species. It focuses on the arguments that may convince breeders to switch to genomic selection. The arguments are classified into three different “bricks.” The first brick considers the addition of genotyping to improve the accuracy of the prediction of breeding values. The second consists of saving costs and/or shortening the breeding cycle by replacing all or a portion of the phenotyping effort with genotyping. The third concerns population management to improve the choice of parents to either optimize crossbreeding or maintain genetic diversity. We analyse the relevance of these different bricks for a wide range of animal and plant species and sought to explain the differences between species according to their biological specificities and the organization of breeding programs.

A comparison of shared patterns of differential gene expression and gene ontologies in response to water-stress in roots and leaves of four diverse genotypes of Lolium and Festuca spp. temperate pasture grasses

Ryegrasses (Lolium spp.) and fescues (Festuca spp.) are closely related and widely cultivated perennial forage grasses. As such, resilience in the face of abiotic stresses is an important component of their traits. We have compared patterns of differentially expressed genes (DEGs) in roots and leaves of two perennial ryegrass genotypes and a single genotype of each of a festulolium (predominantly Italian ryegrass) and meadow fescue with the onset of water stress, focussing on overall patterns of DEGs and gene ontology terms (GOs) shared by all four genotypes. Plants were established in a growing medium of vermiculite watered with nutrient solution. Leaf and root material were sampled at 35% (saturation) and, as the medium dried, at 15%, 5% and 1% estimated water contents (EWCs) and RNA extracted. Differential gene expression was evaluated comparing the EWC sampling points from RNAseq data using a combination of analysis methods. For all genotypes, the greatest numbers of DEGs were identified in the 35/1 and 5/1 comparisons in both leaves and roots. In total, 566 leaf and 643 root DEGs were common to all 4 genotypes, though a third of these leaf DEGs were not regulated in the same up/down direction in all 4 genotypes. For roots, the equivalent figure was 1% of the DEGs. GO terms shared by all four genotypes were often enriched by both up- and down-regulated DEGs in the leaf, whereas generally, only by either up- or down-regulated DEGs in the root. Overall, up-regulated leaf DEGs tended to be more genotype-specific than down-regulated leaf DEGs or root DEGs and were also associated with fewer GOs. On average, only 5–15% of the DEGs enriching common GO terms were shared by all 4 genotypes, suggesting considerable variation in DEGs between related genotypes in enacting similar biological processes.

Genomic prediction and genomic heritability of grain yield and its related traits in a safflower genebank collection

Safflower, a minor oilseed crop, is gaining increased attention for food and industrial uses. Safflower genebank collections are an important genetic resource for crop enhancement and future breeding programs. In this study, we investigated the population structure of a safflower collection sourced from the Australian Grain Genebank and assessed the potential of genomic prediction (GP) to evaluate grain yield and related traits using single and multi-site models. Prediction accuracies (PA) of genomic best linear unbiased prediction (GBLUP) from single site models ranged from 0.21 to 0.86 for all traits examined and were consistent with estimated genomic heritability (h² ), which varied from low to moderate across traits. We generally observed a low level of genome × environment interactions (g × E). Multi-site g × E GBLUP models only improved PA for accessions with at least some phenotypes in the training set. We observed that relaxing quality filtering parameters for genotype-by-sequencing (GBS), such as missing genotype call rate, did not affect PA but upwardly biased h² estimation. Our results indicate that GP is feasible in safflower evaluation and is potentially a cost-effective tool to facilitate fast introgression of desired safflower trait variation from genebank germplasm into breeding lines.

The Impact of Alkaloid-Producing Epichloë Endophyte on Forage Ryegrass Breeding: A New Zealand Perspective

For 30 years, forage ryegrass breeding has known that the germplasm may contain a maternally inherited symbiotic Epichloë endophyte. These endophytes produce a suite of secondary alkaloid compounds, dependent upon strain. Many produce ergot and other alkaloids, which are associated with both insect deterrence and livestock health issues. The levels of alkaloids and other endophyte characteristics are influenced by strain, host germplasm, and environmental conditions. Some strains in the right host germplasm can confer an advantage over biotic and abiotic stressors, thus acting as a maternally inherited desirable ‘trait’. Through seed production, these mutualistic endophytes do not transmit into 100% of the crop seed and are less vigorous than the grass seed itself. This causes stability and longevity issues for seed production and storage should the ‘trait’ be desired in the germplasm. This makes understanding the precise nature of the relationship vitally important to the plant breeder. These Epichloë endophytes cannot be ‘bred’ in the conventional sense, as they are asexual. Instead, the breeder may modulate endophyte characteristics through selection of host germplasm, a sort of breeding by proxy. This article explores, from a forage seed company perspective, the issues that endophyte characteristics and breeding them by proxy have on ryegrass breeding, and outlines the methods used to assess the ‘trait’, and the application of these through the breeding, production, and deployment processes. Finally, this article investigates opportunities for enhancing the utilisation of alkaloid-producing endophytes within pastures, with a focus on balancing alkaloid levels to further enhance pest deterrence and improving livestock outcomes.

Natural variation and genomic prediction of growth, physiological traits, and nitrogen-use efficiency in perennial ryegrass under low-nitrogen stress

Genomic prediction of nitrogen-use efficiency (NUE) has not previously been studied in perennial grass species exposed to low-N stress. Here, we conducted a genomic prediction of physiological traits and NUE in 184 global accessions of perennial ryegrass (Lolium perenne) in response to a normal (7.5 mM) and low (0.75 mM) supply of N. After 21 d of treatment under greenhouse conditions, significant variations in plant height increment (ΔHT), leaf fresh weight (LFW), leaf dry weight (LDW), chlorophyll index (Chl), chlorophyll fluorescence, leaf N and carbon (C) contents, C/N ratio, and NUE were observed in accessions , but to a greater extent under low-N stress. Six genomic prediction models were applied to the data, namely the Bayesian method Bayes C, Bayesian LASSO, Bayesian Ridge Regression, Ridge Regression-Best Linear Unbiased Prediction, Reproducing Kernel Hilbert Spaces, and randomForest. These models produced similar prediction accuracy of traits within the normal or low-N treatments, but the accuracy differed between the two treatments. ΔHT, LFW, LDW, and C were predicted slightly better under normal N with a mean Pearson r-value of 0.26, compared with r=0.22 under low N, while the prediction accuracies for Chl, N, C/N, and NUE were significantly improved under low-N stress with a mean r=0.45, compared with r=0.26 under normal N. The population panel contained three population structures, which generally had no effect on prediction accuracy. The moderate prediction accuracies obtained for N, C, and NUE under low-N stress are promising, and suggest a feasible means by which germplasm might be initially assessed for further detailed studies in breeding programs.

Advancing crested wheatgrass [Agropyron cristatum (L.) Gaertn.] breeding through genotyping-by-sequencing and genomic selection

Crested wheatgrass [Agropyron cristatum (L.) Gaertn.] provides high quality, highly palatable forage for early season grazing. Genetic improvement of crested wheatgrass has been challenged by its complex genome, outcrossing nature, long breeding cycle, and lack of informative molecular markers. Genomic selection (GS) has potential for improving traits of perennial forage species, and genotyping-by-sequencing (GBS) has enabled the development of genome-wide markers in non-model polyploid plants. An attempt was made to explore the utility of GBS and GS in crested wheatgrass breeding. Sequencing and phenotyping 325 genotypes representing 10 diverse breeding lines were performed. Bioinformatics analysis identified 827, 3,616, 14,090 and 46,136 single nucleotide polymorphism markers at 20%, 30%, 40% and 50% missing marker levels, respectively. Four GS models (BayesA, BayesB, BayesCπ, and rrBLUP) were examined for the accuracy of predicting nine agro-morphological and three nutritive value traits. Moderate accuracy (0.20 to 0.32) was obtained for the prediction of heading days, leaf width, plant height, clump diameter, tillers per plant and early spring vigor for genotypes evaluated at Saskatoon, Canada. Similar accuracy (0.29 to 0.35) was obtained for predicting fall regrowth and plant height for genotypes evaluated at Swift Current, Canada. The Bayesian models displayed similar or higher accuracy than rrBLUP. These findings show the feasibility of GS application for a non-model species to advance plant breeding.

High-Throughput Genome-Wide Genotyping To Optimize the Use of Natural Genetic Resources in the Grassland Species Perennial Ryegrass (Lolium perenne L.)

The natural genetic diversity of agricultural species is an essential genetic resource for breeding programs aiming to improve their ecosystem and production services. A large natural ecotype diversity is usually available for most grassland species. This could be used to recombine natural climatic adaptations and agronomic value to create improved populations of grassland species adapted to future regional climates. However describing natural genetic resources can be long and costly. Molecular markers may provide useful information to help this task. This opportunity was investigated for Lolium perenne L., using a set of 385 accessions from the natural diversity of this species collected right across Europe and provided by genebanks of several countries. For each of these populations, genotyping provided the allele frequencies of 189,781 SNP markers. GWAS were implemented for over 30 agronomic and/or putatively adaptive traits recorded in three climatically contrasted locations (France, Belgium, Germany). Significant associations were detected for hundreds of markers despite a strong confounding effect of the genetic background; most of them pertained to phenology traits. It is likely that genetic variability in these traits has had an important contribution to environmental adaptation and ecotype differentiation. Genomic prediction models calibrated using natural diversity were found to be highly effective to describe natural populations for almost all traits as well as commercial synthetic populations for some important traits such as disease resistance, spring growth or phenological traits. These results will certainly be valuable information to help the use of natural genetic resources of other species.

The value of early-stage phenotyping for wheat breeding in the age of genomic selection

Genomic selection using data from an on-going breeding program can improve gain from selection, relative to phenotypic selection, by significantly increasing the number of lines that can be evaluated. The early stages of phenotyping involve few observations and can be quite inaccurate. Genomic selection (GS) could improve selection accuracy and alter resource allocation. Our objectives were (1) to compare the prediction accuracy of GS and phenotyping in stage-1 and stage-2 field evaluations and (2) to assess the value of stage-1 phenotyping for advancing lines to stage-2 testing. We built training populations from 1769 wheat breeding lines that were genotyped and phenotyped for yield, test weight, Fusarium head blight resistance, heading date, and height. The lines were in cohorts, and analyses were done by cohort. Phenotypes or GS estimated breeding values were used to determine the trait value of stage-1 lines, and these values were correlated with their phenotypes from stage-2 trials. This was repeated for stage-2 to stage-3 trials. The prediction accuracy of GS and phenotypes was similar to each other regardless of the amount (0, 50, 100%) of stage-1 data incorporated in the GS model. Ranking of stage-1 lines by GS predictions that used no stage-1 phenotypic data had marginally lower correspondence to stage-2 phenotypic rankings than rankings of stage-1 lines based on phenotypes. Stage-1 lines ranked high by GS had slightly inferior phenotypes in stage-2 trials than lines ranked high by phenotypes. Cost analysis indicated that replacing stage-1 phenotyping with GS would allow nearly three times more stage-1 candidates to be assessed and provide 0.84-2.23 times greater gain from selection. We conclude that GS can complement or replace phenotyping in early stages of phenotyping.

Genomic Predictive Ability for Foliar Nutritive Traits in Perennial Ryegrass

Forage nutritive value impacts animal nutrition, which underpins livestock productivity, reproduction and health. Genetic improvement for nutritive traits in perennial ryegrass has been limited, as they are typically expensive and time-consuming to measure through conventional methods. Genomic selection is appropriate for such complex and expensive traits, enabling cost-effective prediction of breeding values using genome-wide markers. The aims of the present study were to assess the potential of genomic selection for a range of nutritive traits in a multi-population training set, and to quantify contributions of family, location and family-by-location variance components to trait variation and heritability for nutritive traits. The training set consisted of a total of 517 half-sibling (half-sib) families, from five advanced breeding populations, evaluated in two distinct New Zealand grazing environments. Autumn-harvested samples were analyzed for 18 nutritive traits and maternal parents of the half-sib families were genotyped using genotyping-by-sequencing. Significant (P < 0.05) family variance was detected for all nutritive traits and genomic heritability (h2g ) was moderate to high (0.20 to 0.74). Family-by-location interactions were significant and particularly large for water soluble carbohydrate (WSC), crude fat, phosphorus (P) and crude protein. GBLUP, KGD-GBLUP and BayesCπ genomic prediction models displayed similar predictive ability, estimated by 10-fold cross validation, for all nutritive traits with values ranging from r = 0.16 to 0.45 using phenotypes from across two locations. High predictive ability was observed for the mineral traits sulfur (0.44), sodium (0.45) and magnesium (0.45) and the lowest values were observed for P (0.16), digestibility (0.22) and high molecular weight WSC (0.23). Predictive ability estimates for most nutritive traits were retained when marker number was reduced from one million to as few as 50,000. The moderate to high predictive abilities observed suggests implementation of genomic selection is feasible for most of the nutritive traits examined.

The role of forage management in addressing challenges facing Australasian dairy farming

Forage management underpins the viability of pastoral dairy systems. This review investigated recent developments in forage research and their potential to enable pastoral dairy systems to meet the challenges that will be faced over the next 10 years. Grazing management, complementary forages, pasture diversity, fertiliser use, chemical restriction, irrigation management and pasture breeding are considered. None of these areas of research are looking to increase production directly through increased inputs, but, rather, they aim to lift maximum potential production, defend against production decline or improve the efficiency of the resource base and inputs. Technology approaches consistently focus on improving efficiency, while genetic improvement or the use of complementary forages and species diversity aim to lift production. These approaches do not require additional labour to implement, but many will require an increase in skill level. Only a few areas will help address animal welfare (e.g. the use of selected complementary forages and novel endophytes) and only complementary forages will help address increased competition from non-dairy alternatives, by positively influencing the properties of milk. Overall, the diversity of activity and potential effects will provide managers of pastoral dairy systems with the best tools to respond to the production and environmental challenges they face over the next 10 years.

Multi-Trait Genomic Prediction Improves Predictive Ability for Dry Matter Yield and Water-Soluble Carbohydrates in Perennial Ryegrass

In perennial ryegrass (Lolium perenne L), annual and seasonal dry matter yield (DMY) and nutritive quality of herbage are high-priority traits targeted for improvement through selective breeding. Genomic prediction (GP) has proven to be a valuable tool for improving complex traits and may be further enhanced through the use of multi-trait (MT) prediction models. In this study, we evaluated the relative performance of MT prediction models to improve predictive ability for DMY and key nutritive quality traits, using two different training populations (TP1, n = 463 and TP2, n = 517) phenotyped at multiple locations. MT models outperformed single-trait (ST) models by 24% to 59% for DMY and 67% to 105% for nutritive quality traits, such as low, high, and total WSC, when a correlated secondary trait was included in both the training and test set (MT-CV2) or in the test set alone (MT-CV3) (trait-assisted genomic selection). However, when a secondary trait was included in training set and not the test set (MT-CV1), the predictive ability was not statistically significant (p > 0.05) compared to the ST model. We evaluated the impact of training set size when using a MT-CV2 model. Using a highly correlated trait (r_g = 0.88) as the secondary trait in the MT-CV2 model, there was no loss in predictive ability for DMY even when the training set was reduced to 50% of its original size. In contrast, using a weakly correlated secondary trait (r_g = 0.56) in the MT-CV2 model, predictive ability began to decline when the training set size was reduced by only 11% from its original size. Using a ST model, genomic predictive ability in a population unrelated to the training set was poor (r_p = -0.06). However, when using an MT-CV2 model, the predictive ability was positive and high (r_p = 0.76) for the same population. Our results demonstrate the first assessment of MT models in forage species and illustrate the prospects of using MT genomic selection in forages, and other outcrossing plant species, to accelerate genetic gains for complex agronomical traits, such as DMY and nutritive quality characteristics.

Development and Validation of a Model to Combine NDVI and Plant Height for High-Throughput Phenotyping of Herbage Yield in a Perennial Ryegrass Breeding Program

Sensor-based phenotyping technologies may offer a non-destructive, high-throughput and efficient assessment of herbage yield (HY) to replace current inefficient phenotyping methods. This paper assesses the feasibility of combining normalised difference vegetative index (NDVI) from multispectral imaging and ultrasonic sonar estimates of plant height to estimate HY of single plants in a large perennial ryegrass breeding program. For sensor calibration, fresh HY (FHY) and dry HY (DHY) were acquired destructively, and plant height was measured at four dates each in 2017 and 2018 from a selected subset of 480 plants. Global multiple linear regression models based on K-fold and random split cross-validation methods were used to evaluate the relationship between observed vs. predicted HY. The coefficient of determination (R2) = 0.67–0.68 and a root mean square error (RMSE) between 5.43–7.60 g was obtained for the validation of predicted vs. observed DHY. The mean absolute error (MAE) and mean percentage error (MPE) ranged between 3.59–5.44 g and 22–28%, respectively. For the FHY, R2 values ranged from 0.63 to 0.70, with an RMSE between 23.50 and 33 g, MAE between 15.11 and 24.34 g and MPE between ~22% and 31%. Combining NDVI and plant height is a robust method to enable high-throughput phenotyping of herbage yield in perennial ryegrass breeding programs.

Field Spectroscopy to Determine Nutritive Value Parameters of Individual Ryegrass Plants

The nutritive value (NV) of perennial ryegrass is an important driver of productivity for grazing stock; therefore, improving NV parameters would be beneficial to meat and dairy producers. NV is not actively targeted by most breeding programs due to NV measurement being prohibitively slow and expensive. Nondestructive spectroscopy has the potential to reduce the time and cost required to screen for NV parameters to make targeted breeding of NV practical. The application of a field spectrometer was trialed to gather canopy spectra of individual ryegrass plants to develop predictive models for eight NV parameters for breeding programs. The targeted NV parameters included acid detergent fibre, ash, crude protein, dry matter, in vivo dry matter digestibility, in vivo organic matter digestibility, neutral detergent fibre, and water-soluble carbohydrates. The models were developed with partial least square regression. Model predicted ranking of plants had R2 between (0.87 and 0.39) and lab rankings of highest preforming plants. The highest ranked plants, which are generally the selection target for breeding programs, were accurately identified with the canopy-based model at a speed, cost and accuracy that is promising for NV breeding programs.

Simultaneous selection for grain yield and protein content in genomics-assisted wheat breeding

KEY MESSAGE

Large genetic improvement can be achieved by simultaneous genomic selection for grain yield and protein content when combining different breeding strategies in the form of selection indices. Genomic selection has been implemented in many national and international breeding programmes in recent years. Numerous studies have shown the potential of this new breeding tool; few have, however, taken the simultaneous selection for multiple traits into account that is though common practice in breeding programmes. The simultaneous improvement in grain yield and protein content is thereby a major challenge in wheat breeding due to a severe negative trade-off. Accordingly, the potential and limits of multi-trait selection for this particular trait complex utilizing the vast phenotypic and genomic data collected in an applied wheat breeding programme were investigated in this study. Two breeding strategies based on various genomic-selection indices were compared, which (1) aimed to select high-protein genotypes with acceptable yield potential and (2) develop high-yielding varieties, while maintaining protein content. The prediction accuracy of preliminary yield trials could be strongly improved when combining phenotypic and genomic information in a genomics-assisted selection approach, which surpassed both genomics-based and classical phenotypic selection methods both for single trait predictions and in genomic index selection across years. The employed genomic selection indices mitigated furthermore the negative trade-off between grain yield and protein content leading to a substantial selection response for protein yield, i.e. total seed nitrogen content, which suggested that it is feasible to develop varieties that combine a superior yield potential with comparably high protein content, thus utilizing available nitrogen resources more efficiently.

Validation of Genotyping by Sequencing Using Transcriptomics for Diversity and Application of Genomic Selection in Tetraploid Potato

Potato is an important food crop due to its increasing consumption, and as a result, there is demand for varieties with improved production. However, the current status of breeding for improved varieties is a long process which relies heavily on phenotypic evaluation and dated molecular techniques and has little emphasis on modern genotyping approaches. Evaluation and selection before a cultivar is commercialized typically takes 10-15 years. Molecular markers have been developed for disease and pest resistance, resulting in initial marker-assisted selection in breeding. This study has evaluated and implemented a high-throughput transcriptome sequencing method for dense marker discovery in potato for the application of genomic selection. An Australian relevant collection of commercial cultivars was selected, and identification and distribution of high quality SNPs were examined using standard bioinformatic pipelines and a custom approach for the prediction of allelic dosage. As a result, a large number of SNP markers were identified and filtered to generate a high-quality subset that was then combined with historic phenotypic data to assess the approach for genomic selection. Genomic selection potential was predicted for highly heritable traits and the approach demonstrated advantages over the previously used technologies in terms of markers identified as well as costs incurred. The high-quality SNP list also provided acceptable genome coverage which demonstrates its applicability for much larger future studies. This SNP list was also annotated to provide an indication of function and will serve as a resource for the community in future studies. Genome wide marker tools will provide significant benefits for potato breeding efforts and the application of genomic selection will greatly enhance genetic progress.

Boosting Genetic Gain in Allogamous Crops via Speed Breeding and Genomic Selection

Breeding schemes that utilize modern breeding methods like genomic selection (GS) and speed breeding (SB) have the potential to accelerate genetic gain for different crops. We investigated through stochastic computer simulation the advantages and disadvantages of adopting both GS and SB (SpeedGS) into commercial breeding programs for allogamous crops. In addition, we studied the effect of omitting one or two selection stages from the conventional phenotypic scheme on GS accuracy, genetic gain, and inbreeding. As an example, we simulated GS and SB for five traits (heading date, forage yield, seed yield, persistency, and quality) with different genetic architectures and heritabilities (0.7, 0.3, 0.4, 0.1, and 0.3; respectively) for a tall fescue breeding program. We developed a new method to simulate correlated traits with complex architectures of which effects can be sampled from multiple distributions, e.g. to simulate the presence of both minor and major genes. The phenotypic selection scheme required 11 years, while the proposed SpeedGS schemes required four to nine years per cycle. Generally, SpeedGS schemes resulted in higher genetic gain per year for all traits especially for traits with low heritability such as persistency. Our results showed that running more SB rounds resulted in higher genetic gain per cycle when compared to phenotypic or GS only schemes and this increase was more pronounced per year when cycle time was shortened by omitting cycle stages. While GS accuracy declined with additional SB rounds, the decline was less in round three than in round two, and it stabilized after the fourth SB round. However, more SB rounds resulted in higher inbreeding rate, which could limit long-term genetic gain. The inbreeding rate was reduced by approximately 30% when generating the initial population for each cycle through random crosses instead of generating half-sib families. Our study demonstrated a large potential for additional genetic gain from combining GS and SB. Nevertheless, methods to mitigate inbreeding should be considered for optimal utilization of these highly accelerated breeding programs.