Search

Background: Genomic selection has the potential to accelerate genetic gain in perennial ryegrass breeding, provided complex traits such as forage yield can be predicted with sufficient accuracy. Methods: In this study, we compared modelling approaches and feature selection strategies to evaluate the accuracy of genomic prediction models for seasonal forage yield production. Results: Overall, model selection had limited impact on predictive ability when using the full data set. For a baseline genomic best linear unbiased prediction model, the highest mean predictive accuracy was obtained for spring grazing (0.78), summer grazing (0.62) and second cut silage (0.56). In terms of feature selection strategies, using uncorrelated single‐nucleotide polymorphisms (SNPs) had no impact on predictive ability, allowing for a potential decrease of the data set dimensions. With a genome‐wide association study, we found a significant SNP marker for spring grazing, located in the genic region annotated as coding for an enzyme responsible for fucosylation of xyloglucans—major components of the plant cell wall. We also presented an approach to increase interpretability of genomic prediction models with the use of Gene Ontology enrichment analysis. Conclusions: Approaches for feature selection will be relevant in development of low‐cost genotyping platforms in support of routine and cost‐effective implementation of genomic selection. [ABSTRACT FROM AUTHOR]

To address the lack of a truly portable, universal reference mapping population for perennial ryegrass, we have been developing a recombinant inbred line (RIL) mapping population of perennial ryegrass derived via single seed descent from a well-characterized F2 mapping population based on genetically distinct inbred parents in which the natural self-incompatibility (SI) system of perennial ryegrass has been overcome. We examined whether it is possible to create a genotyping by sequencing (GBS) based genetic linkage map in a small population of the F6 generation of this population. We used 41 F6 genotypes for GBS with PstI/MspI-based libraries. We successfully developed a genetic linkage map comprising 6074 SNP markers, placing a further 22080 presence and absence variation (PAV) markers on the map. We examined the resulting genetic map for general and RIL specific features. Overall segregation distortion levels were similar to those experienced in the F2 generation, but segregation distortion was reduced on linkage group 6 and increased on linkage group 7. Residual heterozygosity in the F6 generation was observed at a level of 5.4%. There was a high proportion of chromosomes (30%) exhibiting the intact haplotype of the original inbred parents of the F1 genotype from which the population is derived, pointing to a tendency for chromosomes to assort without recombining. This could affect the applicability of these lines and might make them more suitable for situations where repressed recombination is an advantage. Inter- and intra-chromosomal linkage disequilibrium (LD) analysis suggested that the map order was robust. We conclude that this RIL population, and subsequent F7 and F8 generations will be useful for genetic analysis and phenotyping of agronomic and biological important traits in perennial ryegrass. [ABSTRACT FROM AUTHOR]

Background: Heading and aftermath heading are important traits in perennial ryegrass because they impact forage quality. So far, genome-wide association analyses in this major forage species have only identified a small number of genetic variants associated with heading date that overall explained little of the variation. Some possible reasons include rare alleles with large phenotypic affects, allelic heterogeneity, or insufficient marker density. We established a genome-wide association panel with multiple genotypes from multiple full-sib families. This ensured alleles were present at the frequency needed to have sufficient statistical power to identify associations. We genotyped the panel via partial genome sequencing and performed genome-wide association analyses with multi-year phenotype data collected for heading date, and aftermath heading. Results: Genome wide association using a mixed linear model failed to identify any variants significantly associated with heading date or aftermath heading. Our failure to identify associations for these traits is likely due to the extremely low linkage disequilibrium we observed in this population. However, using single marker analysis within each full-sib family we could identify markers and genomic regions associated with heading and aftermath heading. Using the ryegrass genome we identified putative orthologs of key heading genes, some of which were located in regions of marker-trait associations. Conclusion: Given the very low levels of LD, genome wide association studies in perennial ryegrass populations are going to require very high SNP densities. Single marker analysis within full-sibs enabled us to identify significant marker-trait associations. One of these markers anchored proximal to a putative ortholog of TFL1, homologues of which have been shown to play a key role in continuous heading of some members of the rose family, Rosaceae. [ABSTRACT FROM AUTHOR]

• Background and Aims High density genetic linkage maps that are extensively anchored to assembled genome sequences of the organism in question are extremely useful in gene discovery. To facilitate this process in perennial ryegrass (Lolium perenne L.), a high density single nucleotide polymorphism (SNP)- and presence/absence variant (PAV)-based genetic linkage map has been developed in an F2 mapping population that has been used as a reference population in numerous studies. To provide a reference sequence to which to align genotyping by sequencing (GBS) reads, a shotgun assembly of one of the grandparents of the population, a tenth-generation inbred line, was created using Illumina-based sequencing. • Methods The assembly was based on paired-end Illumina reads, scaffolded by mate pair and long jumping distance reads in the range of 3-40 kb, with >200-fold initial genome coverage. A total of 169 individuals from an F2 mapping population were used to construct PstI-based GBS libraries tagged with unique 4-9 nucleotide barcodes, resulting in 284 million reads, with approx. 1.6 million reads per individual. A bioinformatics pipeline was employed to identify both SNPs and PAVs. A core genetic map was generated using high confidence SNPs, to which lower confidence SNPs and PAVs were subsequently fitted in a straightforward binning approach. • Key Results The assembly comprises 424 750 scaffolds, covering 1.11 Gbp of the 2.5 Gbp perennial ryegrass genome, with a scaffold N50 of 25 212 bp and a contig N50 of 3790 bp. It is available for download, and access to a genome browser has been provided. Comparison of the assembly with available transcript and gene model data sets for perennial ryegrass indicates that approx. 570 Mbp of the gene-rich portion of the genome has been captured. An ultra-high density genetic linkage map with 3092 SNPs and 7260 PAVs was developed, anchoring just over 200Mb of the reference assembly. • Conclusions The combined genetic map and assembly, combined with another recently released genome assembly, represent a significant resource for the perennial ryegrass genetics community. [ABSTRACT FROM AUTHOR]

Background: Genomic selection (GS) can accelerate genetic gains in breeding programmes by reducing the time it takes to complete a cycle of selection. Puccinia coronata f. sp lolli (crown rust) is one of the most widespread diseases of perennial ryegrass and can lead to reductions in yield, persistency and nutritional value. Here, we used a large perennial ryegrass population to assess the accuracy of using genome wide markers to predict crown rust resistance and to investigate the factors affecting predictive ability. Results: Using these data, predictive ability for crown rust resistance in the complete population reached a maximum of 0.52. Much of the predictive ability resulted from the ability of markers to capture genetic relationships among families within the training set, and reducing the marker density had little impact on predictive ability. Using permutation based variable importance measure and genome wide association studies (GWAS) to identify and rank markers enabled the identification of a small subset of SNPs that could achieve predictive abilities close to those achieved using the complete marker set. Conclusion: Using a GWAS to identify and rank markers enabled a small panel of markers to be identified that could achieve higher predictive ability than the same number of randomly selected markers, and predictive abilities close to those achieved with the entire marker set. This was particularly evident in a sub-population characterised by having on-average higher genome-wide linkage disequilibirum (LD). Higher predictive abilities with selected markers over random markers suggests they are in LD with QTL. Accuracy due to genetic relationships will decay rapidly over generations whereas accuracy due to LD will persist, which is advantageous for practical breeding applications. [ABSTRACT FROM AUTHOR]