Skip to main content

Genomic prediction of crown rust resistance in Lolium perenne



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.


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.


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.


Perennial ryegrass (Lolium perenne L.) is the predominant forage species grown in temperate regions of the world [1]. Puccinia coronata f. sp. lolli (crown rust) is one of the most widespread diseases of perennial ryegrass and can lead to a reduction in forage nutritive value, yield and persistency [24]. Poor quality, rust infected swards can impact animal performance and well-being [57]. Developing resistant cultivars is the most viable option for disease control and it has been shown that resistance to crown rust is conferred by both quantitative and qualitative inheritance [811]. As an obligate out-crossing species, perennial ryegrass germplasm has high variation for disease resistance that can be utilized to develop resistant cultivars [1113]. Phenotypic recurrent selection is typically used to develop cultivars with improved resistance and selection is often carried out on spaced plants [9, 11, 12, 14]. There is a high correlation between spaced plants and swards for disease resistance and indirect selection for disease resistance on spaced plants can improve resistance in sward conditions [15]. However, with the advancements in molecular marker development over the last decade, efforts to use marker assisted breeding strategies have been pursued. One such strategy involves identifying quantitative trait loci (QTL) in bi-parental mapping populations and using markers to efficiently backcross the QTL into elite breeding material [16]. Although QTLs explaining significant phenotypic variation for crown rust resistance were mapped onto linkage group (LG) 1-5 and 7 [1723], it is unclear if any of these QTLs were successfully introduced into breeding material. Genome wide association studies (GWAS) are another approach to identify markers linked to QTL. In this case breeding populations can directly be used to identify marker-trait associations, although identified markers tended to explain a small proportion of the total additive genetic variance, resulting in smaller genetic gains [2426].

Genomic selection (GS) was first proposed by Meuwissen et al. [27], as a method to capture complete additive genetic variance using genome wide markers. GS is a form of marker assisted breeding, which accounts for all marker effects across the entire genome to calculate genomic estimated breeding values (GEBVs), which are used to select individual plants for advancement [26]. Use of genome-wide markers will include small effect loci and is ideal for complex traits with low to moderate heritability. In GS, a training population is genotyped with genome wide markers and phenotyped for the trait under selection and models to predict breeding values from marker data are developed. Implementing GS for complex traits like yield and quality is a primary objective of many perennial ryegrass breeding programmes. In contrast to yield and quality traits, the cost (labour and time) of phenotyping for disease resistance is much lower. However, it is important that any GS approaches targeting yield and quality improvements also ensure adequate disease resistance is maintained, particularly where multiple rounds of marker based selections are performed between field evaluations. Opportunities for GS in perennial ryegrass were first reviewed by Hayes et al. [28], and the earliest empirical study was done by Fè et al. [29] for heading date, which confirmed the superiority of GS over marker assisted selection. Later Fè et al. [30], Grinberg et al. [31] and Byrne et al. [32] reported high predictive ability for important agronomical traits in perennial ryegrass. In particular, predictive ability for crown rust reached up to 0.58 [30] when genotypes and phenotypes were evaluated on F2 families. In this study, we evaluated predictive ability for crown rust resistance on individual plants in a large perennial ryegrass population, and assessed factors contributing to predictive ability, such as training population size and marker density. We also performed a GWAS to identify a small to moderately sized panel of markers with good predictive ability for crown rust resistance.


Plant material, phenotyping and genotyping

The training population consists of 30 diploid perennial ryegrass families that have been described previously [32, 33]. Each family consists of 60 genotypes making up a population of 1800 individuals. The complete population consists of ten cultivars, eight full-sib families, eight half-sib families and four ecotypes. Plants were established in a glasshouse and later transplanted to the field in 2013 at Oak Park, Carlow, Ireland (52° 5134.2N, 6° 5503.0W). Plants were grown in two replicates in a partially balanced incomplete block design. Each block consists of 60 test genotypes and 5 check genotypes and was surrounded by a 1.5 m sward consisting of a four way mix of crown rust susceptible perennial ryegrass cultivars. Crown rust was recorded in the years 2014 and 2015 as mean percentage disease score on each plant. Briefly, percentage disease score was obtained by combining scores of percentage of leaves with infection and average percentage of infection on diseased leaves. Scoring was done at multiple time points in September to November without any harvest cuts between scorings (Table 1). We are trying to develop genomic models to identify plants with good resistance to crown rust across the season, and we decided to use all time points for constructing a quantitative summary for crown rust resistance. To do this we calculated AUDPC for each spaced plant in both years. Using multiple time point data, we generated area under disease progress curve (AUDPC) as follows:

$$ A_{k} = \sum_{i=1}^{N_{i}-1} \frac{(y_{i}+y_{i+1})}{2} (t_{i+1} - t_{i}) $$
Table 1 Mean percentage disease score for crown rust resistance at different time points (TP) in Year1 (2014) and Year2 (2015)

where y i is the extent of infection (percentage disease score) at ith observation and t i is the time point at ith observation. N i is the number of individuals in the data set.

Variance components for crown rust were estimated using the R package lme4 (linear mixed-effects models using ’eigen’ and S4) [34]. Broad sense heritability was estimated as follows:

$$ H^{2} = \frac{\sigma_{g}^{2}}{\left(\sigma_{g}^{2}\right) + \left(\sigma_{g*yr}^{2}\right)/2 + \left(\sigma_{res}^{2}\right)/4} $$

where \(\sigma _{g}^{2}\) is the total genetic variance among individuals, \(\sigma _{g*yr}^{2}\) is the variance associated with genotype by year interaction and \(\sigma _{res}^{2}\) is residual variance. With genotype as random effect and year and checks as fixed effects, conditional modes (BLUPs) were calculated in lme4 and used as input for genomic prediction.

Genotyping was done using genotyping by sequencing (GBS) approach described by Elshire et al. [35] and data analysed as described in Byrne et al. [32]. Briefly, genomic DNA was extracted from leaf samples and GBS libraries were prepared using the restriction enzyme ApeKI, libraries were amplified and sequenced on an Illumina Hiseq2000. Panels of SNPs were identified in the complete population, as well as in all sub populations separately (half-sibs, full-sibs, ecotypes, cultivars). Individuals with very low sequencing coverage and/or largely missing phenotypic data were eliminated from the analysis giving a final population for analysis of 1582 individuals. Missing marker data was imputed using mean imputation.

Genomic prediction models

We used four statistical algorithms for genomic prediction, ridge regression best linear unbiased prediction (rrBLUP) [27], Bayes B [36] and Bayesian Lasso [37], and random forest [38]. rrBLUP is a mixed model approach, which was initially proposed for GS. We used an R package called rrBLUP [39] for fitting the mixed model as follows

$$ y = \mu + Xg + \epsilon $$

where μ is the overall mean, X is the marker matrix, g is the matrix of marker effects, ε is a vector of residual effects and y is a vector of conditional modes for crown rust. We also evaluated two Bayesian approaches, Bayes B [36] and Bayesian Lasso [37], which were both implemented using the R package BGLR [40] with the following parameters: number of iterations = 5000, burn-in = 500 and thinning = 5. Random forest is a machine-learning tool, in which series of regression trees were grown independently to the largest extent possible using subsets of bootstrap samples. At each split of the tree, a random subset of variables is selected to identify the best split. We implemented random forest using the R package randomForest [41], setting the number of variables at each split to 1/3 of the total variables, and using a terminal node size of five and minimum of 500 trees per forest. We also used random forest to rank variables using the variable importance measure, a permutation based measure in which variables are ranked based on the mean decrease in accuracy.

Cross validation scheme

We evaluated genomic prediction models using Monte-Carlo cross-validation by randomly assigning plants into training (70%) and test (30%) sets and the procedure was repeated 100 times and the resulting accuracies were then averaged. This approach to cross-validation has previously been used to evaluate genomic prediction models [42, 43]. Predictive ability and bias were assessed in the complete population and in each sub-population. Predictive ability (r p ) was determined as the Pearson’s correlation coefficient between observed phenotypic value and predicted phenotype. Bias was evaluated by regressing observed phenotypic value on predictions. We reduced training population size and marker density in order to identify the impact of training population size and marker number on predictive ability. To compare predictive ability for traits with contrasting genetic architecture we compared heading date, a highly heritable trait, with crown rust. Predictive ability for heading date has already been shown to be high (0.81) in this population [32]. We re-analyzed data for heading date according to methods described above and made a comparison with crown rust. To evaluate the impact of leaving related material out of the training set we also performed cross validation by leaving one family out. In this approach one complete family (up to 60 individuals) is left out of the training set and only used for testing. This was repeated so that each family in turn is used as a test set.

Genome wide association

A mixed linear model (MLM) was also used for association mapping, implemented in the R package rrBLUP [39]. Population structure and family relatedness was accounted for with a kinship matrix calculated by rrBLUP from the input genotypic data. We accounted for multiple testing using a Bonferroni correction and markers passing an α level 0.05 threshold were considered statistically significant.

Results and discussion

Phenotypic analysis for crown rust

The mean percentage disease score for crown rust infection in the population increased over time in both evaluation years as infection levels accumulated (Table 1). In both years, evaluations were carried out in the period from September to November during a time when disease pressure tends to be at its greatest [15, 44]. The highest mean percentage disease score was seen in late October 2015 and was more than double the highest mean percentage disease score from 2014 (Table 1). In addition to plant health and level of host resistance, crown rust infection is influenced by various environmental factors, such as temperature, relative humidity, and light [4547]. The latency period is reduced and spore production increased as temperature increases [45], and it has been shown that when temperatures exceed 25°C, the susceptibility of previously resistant cultivars can be increased [46]. It has already been shown that there is variability within pathogen populations, and different races can be found within and between locations. It is also possible that the composition of a pathogen population can change over short periods of time and plants that are resistant at one point in time will become susceptible as the pathogen population shifts or evolves.

AUDPC values ranged from 0 to 1371 and the Pearson correlation co-efficient between replicates within years was moderate (0.69 in 2014 and 0.59 in 2015). However, the Pearson correlation co-efficient between years was low (0.28), and there was a significant genotype by year interaction (F(1761)= 3.025, MSE = 60676, p = 0.0001). The broad sense heritability for crown rust infection was moderate (0.36), which is in line with previous estimates of heritability calculated in other populations [11, 48]. Overall there is a good phenotypic variation for crown rust infection among and within the 30 families/cultivars/ecotypes making up the entire population (Fig. 1). Plants were placed into one of four categories (sub-populations) based on mating type or origin, these were (i) full-sib families, (ii) half-sib families, (iii) cultivars, and (iv) ecotypes. In general the ecotypes were more susceptible to crown rust infection than cultivars or breeding material (Fig. 1), which presents a challenge for the incorporation of ecotypes into breeding programmes. The broad-sense heritability calculated in each sub-population varied between 0.17 in the cultivars to 0.44 in the full-sib families.

Fig. 1

Phenotypic variation for crown rust resistance in the complete population, grouped according to sub-population types: cultivars (CS), ecotypes (ES), full-sibs (FS) and half-sibs (HS). Broad sense heritability (H2) in complete population and sub-populations is highlighted over the figure

Crown rust infection is typically evaluated in breeding programmes by growing spaced plants or potted plants from a population and visually scoring the level of crown rust infection. A mean score is assigned to each family and used to aid selection of the top performing families from which to construct the synthetics cultivars. During construction of synthetics a spaced plant nursery may be established to evaluate heading date and crown rust resistance before selecting individual genotypes from which to construct synthetics (within family selection). In practice, this has a time cost of 2 to 3 years (establishment, evaluation, selection and recombining), and using molecular markers offers an opportunity to reduce this to one year in those selection cycles where GEBVs are predicted. This depends on our ability to accurately predict traits such as crown rust from genomic data.

Predicting crown rust resistance with genomic data

We evaluated four algorithms for prediction of crown rust infection from genomic data, rrBLUP, Bayes B, Bayesian Lasso, and random forest. The mean predictive ability after cross-validation within the complete population was 0.52 using rrBLUP, 0.52 using Bayesian Lasso, 0.51 using Bayes B, and 0.49 using random forest (Additional file 1: Figure S1). rrBLUP was computationally faster, and therefore results from all further analysis are only reported for rrBLUP. The predictive ability of 0.52 is in line with previous estimates reported in perennial ryegrass where predictions were based on mean genotypes and phenotypes of F2 families [30]. Predictive ability did not differ depending on whether the equations were developed using phenotypes from the last time point scored or the AUDPC values incorporating all time points. This indicates that a single scoring each year would have sufficed. However, the importance of evaluating crown rust in more than one year was emphasised by the low correlation between scores in 2014 and 2015.

When we calculate the predictive ability within each of the sub-populations (cultivars, half-sib families, full-sib families, and ecotypes), the highest predictive ability for crown rust was obtained using plants from full-sib families (0.54) and the lowest predictive ability for crown rust was obtained with the plants from the ecotypes (0.24) (Fig. 2). Generally, traits with higher heritability achieve higher predictive abilities [49, 50], and we see that here where crown rust measurements taken in the full-sib families had the highest broad-sense heritability and the highest predictive ability. In general, there was a good correlation between predictive ability and both phenotypic variance and heritability. This relationship between phenotypic variance and predictive ability has been observed previously [49, 51].

Fig. 2

Predictive ability in different population types. Complete population (CP), cultivars (CS), ecotypes (ES), full-sibs (FS) and half-sibs (HS) are listed on x-axis, predictive ability (left) and bias (right) on y-axis. Crown rust is in red and heading date in blue

We also evaluated the predictive ability using a leave-one-family-out cross validation scheme. The complete population is comprised of 30 families/cultivars/ecotypes, each with up to 60 individual genotypes. The predictive ability was assessed in the complete population by selectively leaving one family out of the training set and using it for testing. In addition to crown rust we also evaluated predictive ability for heading date phenotypes previously reported [32]. The predictive ability for both crown rust (r p =0.02, min =−0.36, max = 0.36) and heading date (r p =0.29, min =−0.14, max = 0.65) varied greatly depending on which family was left out, and having related material in the training set (shared parentage) greatly improved predictive ability.

Effect of training population size and marker density on predictive ability

As we reduced the number of individuals in the training population we saw a decrease in predictive ability and an increasingly upward bias in the variance of predictions for both crown rust resistance and heading date (Fig. 3). The drop in predictive ability was more pronounced as we reduced the training population size for crown rust resistance than it was for heading date. The predictive ability for crown rust resistance when using 90% of the population as a training set was 0.52 and the predictive ability was 0.38 when using just 10% of the population. Irrespective of the trait, as the training population size increased there was an increase in predictive ability which is consistent with similar correlations between training population size and predictive ability reported previously for perennial ryegrass [29, 30] and other crops [5154]. Useful linkage disequilibrium (LD) only extends over short distances in perennial ryegrass and it has been suggested that this is the result of a very large past effective population size [28]. This impacts both the size of the reference population and marker density required to achieve high accuracies when predicting traits from genomic data. The fact that we are able to achieve high predictive abilities with relatively small training populations is likely a result of strong genetic structure and differentiation in our diverse population and the use of the marker data to capture genetic relationships [55].

Fig. 3

Effect of training population size on predictive ability. Training population is varied from 90% (1423 individuals) to 10% (158 individuals) on x-axis and predictive ability (left), bias (right) on y-axis. Crown rust is in red and heading date in blue

The limited LD also affects the number of markers required to obtain high predictive accuracies, and given the extent of LD in the broader perennial ryegrass population, marker numbers in excess of one million have been suggested for achieving high accuracies [28]. When we reduced marker number in the complete population and the various sub-populations we observed very little impact on the predictive ability for either trait (Table 2). Reducing the marker set to 5% of the total available had virtually no impact on predictive ability in all cases. This would support our observation that much of the predictive ability can be derived from makers capturing familial relationships. When marker number dropped below 5% (10878) predictive ability for both traits in the complete population began to drop. However, even with 0.05% (109) of markers the mean predictive ability was 0.30 for crown rust resistance and 0.52 for heading date. Knowing the contribution of genetic relationships to predictive ability is important because it will change over generations. In contrast, predictive ability due to LD has greater persistence over generations and is therefore preferential [55]. Schemes for implementing genomic selection in perennial ryegrass that pursue a reduction in effective population size from the outset have been proposed. Such schemes would lead to an increase in the extent of LD and ensure that predictive ability due to LD can be captured using a reasonable number of markers and a reference population size that is feasible in breeding programmes.

Table 2 Predictive ability (r p ) and bias for crown rust (CR) and heading date (HD) by selecting random markers of 100 to 0.05%, in complete population (CP), cultivars (CS), full-sibs (FS) and half-sibs (HS)

Identifying SNPs associated with crown rust resistance

The cost of genotyping impacts the number of selection candidates that can be evaluated and therefore impacts the selection intensity. Different approaches to low density SNP genotyping for genomic selection have been proposed. These include variable selection methods to identify a small subset of markers in strong LD with the trait [56] or using a small random subset of markers to impute from low-to-high density [57]. Until a chromosome scale assembly of the perennial ryegrass genome becomes available the latter remains a challenge. We used both permutation based variable importance measures and GWAS analysis to identify a subset of markers capable of predicting crown rust resistance. Using permutation based variable importance measures we were able to rank markers by mean decrease in accuracy and select the top ranked markers for use in genomic prediction. In the case of GWAS we ranked SNPs based on significance and again selected the top ranked markers for use in genomic prediction. All variable importance measures and GWAS were identified and ranked in the training set and used to predict phenotypes in the test set via cross-validation. When we used the top 100 ranked markers from the permutation based variable importance measures, the mean predictive ability of 100 iterations was 0.42 (ranging from 0.36 to 0.48). When we used the top 100 ranked markers from the GWAS analysis, the mean predictive ability of 100 iterations was 0.36 (ranging from 0.25 to 0.44). In both cases the mean predictive ability with selected markers is higher than the predictive ability with random markers, which was 0.28 (ranging from 0.18 to 0.39). The lower predictive ability using GWAS marker selection is not surprising considering that we corrected for population structure using a kinship matrix, and we are more reliant on identifying markers in LD with the trait. As discussed above, the predictive ability of these markers is expected to be more persistent over subsequent generations. Using GWAS selected markers it is clear to see that they are superior to randomly selected markers up to the point, beyond which adding more markers does not improve predictive ability in either case (Fig. 4). The ability of a GWAS within each sub-population to identify and select a small set of SNPs with excellent predictive ability varied, and in some cases was little better than random SNP selection (Fig. 5). The GWAS on plants originating from IBERs bred cultivars identified a small set of twenty SNPs with 77% of the predictive ability achieved with 20,000 SNPs. The power of a GWAS to identify markers with high predictive ability was much greater within the population made up of IBERs plants than within cultivars, and full-sib families where twenty SNPs could only achieve 46 and 48% of the predictive ability with 20,000 SNPs, respectively. On average LD is higher within the sub-population with IBERs plants, which may explain the greater ability to identify markers associated with crown-rust resistance.

Fig. 4

Predictive ability of selected markers versus random markers in the complete population. Markers were selected based on the ranking from genome wide association studies and compared with random markers of similar size

Fig. 5

Comparing predictive ability of selected versus random markers. Markers were selected based on the ranking from genome wide association studies in cultivars, full-sibs and IBERS material and compared with random markers of similar size

In order to characterise the markers associated with crown rust resistance we repeated the GWAS analysis without division of genotypes into training and testing sets. We carried out GWAS using the complete population and found 29 markers significantly associated with crown rust resistance after correction for multiple testing (Additional file 2: Table S1). Using the perennial ryegrass genome [58] as a reference, we located all markers within 22 genomic scaffolds that contained 50 predicted genes. Using the Genome Zippper [58, 59], we anchored ten scaffolds onto LG2, 3, 4, 5 and 7 (Additional file 3: Table S2). Similarly, we did GWAS on IBERS material and found 24 markers associated with crown rust resistance (Additional file 2: Table S1). All markers were located within 16 genomic scaffolds containing 44 predicted genes. Out of 16 scaffolds we were able to place seven scaffolds onto LG3, 5 and 7 (Additional file 3: Table S2). We found five common scaffolds between the complete population and the IBERS and only two of these scaffolds were mapped, onto LG3. On LG3 five markers were anchored within 60.4-61.21 cM. Genes present on these scaffolds were coding for domains including Mon1, Aquaporin, DUF1635, Nucleoredoxin, Beta-glucan export ATP-binding/permease protein, BRASSINOSTEROID INSENSITIVE 1-associated receptor kinase 1, Alpha N-terminal protein methyltransferase 1. Gene function of these domains plays a key role in ATP-binding, membrane proteins, enzyme catalysis and pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) (Additional file 4: Table S3) [60].

Using small subsets of trait associated markers may be an effective strategy for within-family prediction of traits such as heading date, crown rust resistance and some quality traits. Predicting heading date from markers would enable plants to be matched in heading date to ensure sufficient cross-pollination when constructing synthetic cultivars [32]. Combining these with markers to predict crown rust resistance would also avoid the inclusion of plants with high levels of susceptibility, and furthermore prediction models can be based on multi-year evaluations. It is clear from the phenotypic data presented here that there is substantial within family variation for crown rust resistance. Opportunities already exist to genotype small to moderate sized marker panels in 1000s of samples at low cost [61]. Using these approaches small fragments (200-300 bp) are amplified and sequenced at hundreds of loci. These amplicons can be used as short haplotypes in marker aided selection strategies. An assay can be developed to target loci in linkage with QTL for heading date [32], crown rust resistance, and other traits such as quality parameters. The assay can also include a suite of loci with a good distribution throughout the genome and be deployed for among-and-within-full-sib-family selection (Additional file 5: Figure S2). Once high yielding families are identified in field trials, within family selection for crown rust resistance and forage quality can be performed at a high selection intensity with the molecular marker assay. Furthermore, plants can be selected to be synchronous in flowering time.


Our findings show that predicting crown rust resistance in perennial ryegrass can be achieved with high accuracy using AUDPC scores on spaced plants. However, there was no difference in predictive ability when equations were developed using phenotypes from the last time-point scored or the AUDPC values, meaning a single time point was adequate to evaluate the crown rust susceptibility of the spaced plants. Much of the predictive ability comes from markers capturing familial relationships, highlighted by the observation that there was no drop in predictive ability when going from the entire marker set down to only 5% (10,878) of the marker set. Accuracy due to genetic relationships will decay rapidly over generations whereas accuracy due to LD will persist. Using a GWAS we attempted to identify and rank markers in LD with QTL. This enabled a small panel of markers to be identified that had higher predictive ability than the same number of randomly selected markers, and had predictive abilities close to those achieved with the entire marker set.



Area under disease progress curve


Genomic estimated breeding value


Genomic selection


Genome wide association study


Linkage disequilibrium


Quantitative trait locus


  1. 1

    Humphreys M, Feuerstein U, Vandewalle M, Baert J. Ryegrasses. In: Fodder Crops and Amenity Grasses. New York: Springer: 2010. p. 211–60.

    Google Scholar 

  2. 2

    Potter LR. Effect of crown rust on regrowth, competitive ability and nutritional quality of perennial and Italian ryegrasses. Plant Pathol. 1987; 36(4):455–61.

    Article  Google Scholar 

  3. 3

    Valk P, Proveniers MCG, Pertijs JH, Lamers JTW, Dun CMP, Smeekens JCM. Late heading of perennial ryegrass caused by introducing an Arabidopsis homeobox gene. Plant Breed. 2004; 123(6):531–5.

    Article  Google Scholar 

  4. 4

    Price T. Ryegrass rust in victoria. Plant Protect Quart. 1987; 2:189.

    Google Scholar 

  5. 5

    Lancashire JA, Latch GCM. The effect of crown rust (Puccinia coronata Corda) on the yield and botanical composition of two ryegrass/white clover pastures. N Z J Agric Res. 1970; 13(2):279–86.

    Article  Google Scholar 

  6. 6

    Plummer RM, Hall RL, Watt TA. The influence of crown rust (Puccinia coronata) on tiller production and survival of perennial ryegrass (Lolium perenne) plants in simulated swards. Grass Forage Sci. 1990; 45(1):9–16.

    Article  Google Scholar 

  7. 7

    Smit HJ, Tas BM, Taweel HZ, Tamminga S, Elgersma A. Effects of perennial ryegrass (Lolium perenne L,) cultivars on herbage production, nutritional quality and herbage intake of grazing dairy cows. Grass Forage Sci. 2005; 60(3):297–309.

    Article  Google Scholar 

  8. 8

    Hayward MD. Genetic control of resistance to crown rust (Puccinia coronato Corda) in Lolium perenne L, and its implications in breeding. Theor Appl Genet. 1977; 51(2):49–53.

    Article  PubMed  CAS  Google Scholar 

  9. 9

    Kimbeng CA. Genetic basis of crown rust resistance in perennial ryegrass, breeding strategies, and genetic variation among pathogen populations: a review. Anim Prod Sci. 1999; 39(3):361–78.

    Article  Google Scholar 

  10. 10

    McVeigh K. Breeding for resistance to crown rust (Puccinia coronata Corda var lolii Brown) in turf-type perennial ryegrass (Lolium perenne). New Brunswick: Rutgers University; 1975. PhD thesis.

    Google Scholar 

  11. 11

    Reheul D, Ghesquiere A. Breeding perennial ryegrass with better crown rust resistance. Plant Breed. 1996; 115(6):465–9.

    Article  Google Scholar 

  12. 12

    Gibbs JG. Field resistance in Lolium sp to leaf rust (Puccinia coronata). Nature. 1966; 209(5021):420.

    Article  Google Scholar 

  13. 13

    McGrath S, Hodkinson TR, Charles TM, Zen D, Barth S. Variation in inflorescence characters and inflorescence development in ecotypes and cultivars of lolium perenne l. Grass Forage Sci. 2010; 65(4):398–409.

    Article  Google Scholar 

  14. 14

    Wilkins PW. Specialisation of crown rust on highly and moderately resistant plants of perennial ryegrass. Ann Appl Biol. 1978; 88(1):179–84.

    Article  Google Scholar 

  15. 15

    Easton HS, Cooper BM, Frasers TJ, Widdup KH. Crown rust on perennial ryegrass in field trials. In: Proceedings of the New Zealand Grassland Association, vol. 50.New Zealand: New Zealand Grassland Association: 1989. p. 253–4.

    Google Scholar 

  16. 16

    Collard BC, Mackill DJ. Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Phil Trans R Soc B Biol Sci. 2008; 363(1491):557–72.

    Article  CAS  Google Scholar 

  17. 17

    Dumsday JL, Smith KF, Forster JW, Jones ES. SSR-based genetic linkage analysis of resistance to crown rust (Puccinia coronata f sp lolii) in perennial ryegrass (Lolium perenne). Plant Pathol. 2003; 52(5):628–37.

    Article  Google Scholar 

  18. 18

    Dracatos PM, Cogan NOI, Dobrowolski MP, Sawbridge TI, Spangenberg GC, Smith KF, Forster JW. Discovery and genetic mapping of single nucleotide polymorphisms in candidate genes for pathogen defence response in perennial ryegrass (Lolium perenne L,). Theor Appl Genet. 2008; 117(2):203–19.

    Article  PubMed  CAS  Google Scholar 

  19. 19

    Dracatos PM, Cogan NO, Sawbridge TI, Gendall AR, Smith KF, Spangenberg GC, Forster JW. Molecular characterisation and genetic mapping of candidate genes for qualitative disease resistance in perennial ryegrass (Lolium perenne L,). BMC Plant Biol. 2009; 9(1):1.

    Article  CAS  Google Scholar 

  20. 20

    Muylle H, Baert J, Van Bockstaele E, Moerkerke B, Goetghebeur E, Roldàn-Ruiz I. Identification of molecular markers linked with crown rust (Puccinia coronata f sp lolii) resistance in perennial ryegrass (Lolium perenne) using AFLP markers and a bulked segregant approach. Euphytica. 2005; 143(1-2):135–44.

    Article  CAS  Google Scholar 

  21. 21

    Muylle H, Baert J, Van Bockstaele E, Pertijs J, Roldàn-Ruiz I. Four QTLs determine crown rust (Puccinia coronata f sp lolii) resistance in a perennial ryegrass (Lolium perenne) population. Heredity. 2005; 95(5):348–57.

    Article  PubMed  CAS  Google Scholar 

  22. 22

    Schejbel B, Jensen LB, Xing Y, Lübberstedt T. QTL analysis of crown rust resistance in perennial ryegrass under conditions of natural and artificial infection. Plant Breed. 2007; 126(4):347–52.

    Article  Google Scholar 

  23. 23

    Thorogood D, Paget MF, Humphreys MO, Turner LB, Armstead I, Roderick H. QTL analysis of crown rust resistance in perennial ryegrass-implications for breeding. Int Turfgrass Soc Res J. 2001; 9:218–23.

    Google Scholar 

  24. 24

    Hayes B, Goddard M. Genome-wide association and genomic selection in animal breeding this article is one of a selection of papers from the conference “exploiting genome-wide association in oilseed brassicas: a model for genetic improvement of major OECD crops for sustainable farming”. Genome. 2010; 53(11):876–83.

    Article  PubMed  CAS  Google Scholar 

  25. 25

    Heffner EL, Sorrells ME, Jannink J-L. Genomic selection for crop improvement. Crop Sci. 2009; 49(1):1–12.

    Article  CAS  Google Scholar 

  26. 26

    Jannink J-L, Lorenz AJ, Iwata H. Genomic selection in plant breeding: from theory to practice. Brief Funct Genom. 2010; 9(2):166–77.

    Article  CAS  Google Scholar 

  27. 27

    Meuwissen TH, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001; 157(4):1819–29.

    PubMed  PubMed Central  CAS  Google Scholar 

  28. 28

    Hayes BJ, Cogan NO, Pembleton LW, Goddard ME, Wang J, Spangenberg GC, Forster JW. Prospects for genomic selection in forage plant species. Plant Breed. 2013; 132(2):133–43.

    Article  Google Scholar 

  29. 29

    Fè D, Cericola F, Byrne S, Lenk I, Ashraf BH, Pedersen MG, Roulund N, Asp T, Janss L, Jensen CS, et al. Genomic dissection and prediction of heading date in perennial ryegrass. BMC Genomics. 2015; 16(1):1.

    Article  CAS  Google Scholar 

  30. 30

    Fè D, Ashraf BH, Pedersen MG, Janss L, Byrne S, Roulund N, Lenk I, Didion T, Asp T, Jensen CS, et al. Accuracy of genomic prediction in a commercial perennial ryegrass breeding program. Plant Genome. 2016;9(3).

  31. 31

    Grinberg NF, Lovatt A, Hegarty M, Lovatt A, Skøt KP, Kelly R, Blackmore T, Thorogood D, King RD, Armstead I, et al. Implementation of genomic prediction in Lolium perenne (L) breeding populations. Front Plant Sci. 2016;7(133).

  32. 32

    Byrne S, Conaghan P, Barth S, Arojju SK, Casler M, Thibauld M, Velmurugan J, Milbourne D. Using variable importance measures to identify a small set of SNPs to predict heading date in perennial ryegrass. Sci Rep. 2017;7(3566).

  33. 33

    Arojju SK, Barth S, Milbourne D, Conaghan P, Velmurugan J, Hodkinson TR, Byrne SL. Markers associated with heading and aftermath heading in perennial ryegrass full-sib families. BMC Plant Biol. 2016; 16(1):160.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34

    Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. arXiv preprint arXiv:14065823. 2014.

  35. 35

    Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE. 2011; 6(5):19379.

    Article  CAS  Google Scholar 

  36. 36

    Habier D, Fernando RL, Kizilkaya K, Garrick DJ. Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011; 12(1):186.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37

    Park T, Casella G. The bayesian lasso. J Am Stat Assoc. 2008; 103(482):681–6.

    Article  CAS  Google Scholar 

  38. 38

    Breiman L. Random forests. Mach Learn. 2001; 45(1):5–32.

    Article  Google Scholar 

  39. 39

    Endelman JB. Ridge regression and other kernels for genomic selection with R package rrBLUP. The Plant Genome. 2011; 4(3):250–5.

    Article  Google Scholar 

  40. 40

    de los Campos G, Pérez-Rodríguez P. Bayesian generalized linear regression. R package version 1.0. 4. 2014.

  41. 41

    Liaw A, Wiener M. Classification and regression by randomForest. R news. 2002; 2(3):18–22.

    Google Scholar 

  42. 42

    Annicchiarico P, Nazzicari N, Li X, Wei Y, Pecetti L, Brummer EC. Accuracy of genomic selection for alfalfa biomass yield in different reference populations. BMC Genom. 2015; 16(1):1020.

    Article  CAS  Google Scholar 

  43. 43

    Crossa J, Beyene Y, Kassa S, Pérez P, Hickey JM, Chen C, De Los Campos G, Burgueño J, Windhausen VS, Buckler E, et al. Genomic prediction in maize breeding populations with genotyping-by-sequencing. G3: Genes, Genomes, Genetics. 2013; 3(11):1903–26.

    Article  CAS  Google Scholar 

  44. 44

    Schubiger FX, Baert J, Cagas B, Cernoch V, Chosson JF, Czembor E, Eickmeyer F, Feuerstein U, Hartmann S, Jakesova H, et al. The eucarpia multi-site rust evaluation–results 2007. In: Sustainable Use of Genetic Diversity in Forage and Turf Breeding. New York: Springer: 2010. p. 331–40.

    Google Scholar 

  45. 45

    Roderick HW, Thomas BJ. Infection of ryegrass by three rust fungi (Puccinia coronata, P graminis and P loliina) and some effects of temperature on the establishment of the disease and sporulation. Plant Pathol. 1997; 46(5):751–61.

    Article  Google Scholar 

  46. 46

    Roderick HW, Thorogood D, Adomako B. Temperature-dependent resistance to crown rust infection in perennial ryegrass, Lolium perenne. Plant Breed. 2000; 119(1):93–5.

    Article  Google Scholar 

  47. 47

    Schein RD, Rotem J. Temperature and humidity effects on uredospore viability. Mycologia. 1965; 57(3):397–403.

    Article  Google Scholar 

  48. 48

    Fè D, Pedersen MG, Jensen CS, Jensen J. Genetic and environmental variation in a commercial breeding program of perennial ryegrass. Crop Sci. 2015; 55(2):631–40.

    Article  Google Scholar 

  49. 49

    Muranty H, Troggio M, Sadok IB, Al Rifaï M, Auwerkerken A, Banchi E, Velasco R, Stevanato P, Van De Weg WE, Di Guardo M, et al. Accuracy and responses of genomic selection on key traits in apple breeding. Hortic Res. 2015; 2:15060.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 50

    Wimmer V, Lehermeier C, Albrecht T, Auinger H-J, Wang Y, Schön CC. Genome-wide prediction of traits with different genetic architecture through efficient variable selection. Genetics. 2013; 195(2):573–87.

    Article  PubMed  PubMed Central  Google Scholar 

  51. 51

    Isidro J, Jannink J-L, Akdemir D, Poland J, Heslot N, Sorrells ME. Training set optimization under population structure in genomic selection. Theor Appl Genet. 2015; 128(1):145–58.

    Article  PubMed  Google Scholar 

  52. 52

    Lorenz AJ, Smith KP, Jannink JL. Potential and optimization of genomic selection for fusarium head blight resistance in six-row barley. Crop Sci. 2012; 52(4):1609–21.

    Article  Google Scholar 

  53. 53

    Riedelsheimer C, Endelman JB, Stange M, Sorrells ME, Jannink J-L, Melchinger AE. Genomic predictability of interconnected biparental maize populations. Genetics. 2013; 194(2):493–503.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54

    Tayeh N, Klein A, Le Paslier MC, Jacquin F, Houtin H, Rond C, Chabert-Martinello M, Magnin-Robert JB, Marget P, Aubert G, et al. Genomic prediction in pea: effect of marker density and training population size and composition on prediction accuracy. Front Plant Sci. 2015; 6:941.

    PubMed  PubMed Central  Google Scholar 

  55. 55

    Habier D, Fernando RL, Dekkers JCM. The impact of genetic relationship information on genome-assisted breeding values. Genetics. 2007; 177(4):2389–97.

    PubMed  PubMed Central  CAS  Google Scholar 

  56. 56

    Zhang Z, Ding X, Liu J, Zhang Q, de Koning D-J. Accuracy of genomic prediction using low-density marker panels. J Dairy Sci. 2011; 94(7):3642–50.

    Article  PubMed  CAS  Google Scholar 

  57. 57

    Habier D, Fernando RL, Dekkers JCM. Genomic selection using low-density marker panels. Genetics. 2009; 182(1):343–53.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. 58

    Byrne SL, Nagy I, Pfeifer M, Armstead I, Swain S, Studer B, Mayer K, Campbell JD, Czaban A, Hentrup S, et al. A synteny-based draft genome sequence of the forage grass lolium perenne. Plant J. 2015; 84(4):816–26.

    Article  PubMed  CAS  Google Scholar 

  59. 59

    Pfeifer M, Martis M, Asp T, Mayer KF, Lübberstedt T, Byrne S, Frei U, Studer B. The perennial ryegrass genomezipper: targeted use of genome resources for comparative grass genomics. Plant Physiol. 2013; 161(2):571–82.

    Article  PubMed  CAS  Google Scholar 

  60. 60

    Postel S, Küfner I, Beuter C, Mazzotta S, Schwedt A, Borlotti A, Halter T, Kemmerling B, Nürnberger T. The multifunctional leucine-rich repeat receptor kinase bak1 is implicated in arabidopsis development and immunity. Eur J Cell Biol. 2010; 89(2):169–74.

    Article  PubMed  CAS  Google Scholar 

  61. 61

    Campbell NR, Harmon SA, Narum SR. Genotyping-in-thousands by sequencing (GT-seq): A cost effective SNP genotyping method based on custom amplicon sequencing. Mol Ecol Resour. 2015; 15(4):855–67.

    Article  PubMed  CAS  Google Scholar 

  62. 62

    Proost S, Van Bel M, Vaneechoutte D, Van de Peer Y, Inzé D, Mueller-Roeber B, Vandepoele K. Plaza 30: an access point for plant comparative genomics. Nucleic Acids Res. 2015; 43(D1):974–81.

    Article  CAS  Google Scholar 

Download references


Authors would like to thank Jean-Baptiste Enjelvin, Olivia Aylesbury, Mary O’Sullivan, Padraig Dunne and Michael Murphy for help with phenotyping and Sean Murray and Helena Meally for genomic DNA extraction.


This work received funding from the Irish Department of Agriculture Food and the Marine DAFM (RSF 11/S/109) and Teagasc core funding. SKA is supported by a Teagasc PhD Walsh Fellowship. SLB has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 658031. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

All sequence data associated with this project has been deposited in the NCBI Short Read Archive (SRA) under the BioProject accession PRJNA352789. Phenotypic data is available on Figshare (, along with an R script template used for cross-validation (

Author information




DM, PC, SB, TRH, MC conceived and designed the study. SKA and SLB performed the data analysis and drafted the initial manuscript. SKA, SLB, DM, PC, TRH, MC, TM and SB contributed to interpretation of data and preparation of the final manuscript. All authors read and approved the final version.

Corresponding author

Correspondence to Stephen L. Byrne.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1

Figure S1. Predictive ability and bias for crown rust using various algorithms for genomic prediction. (PDF 101 kb)

Additional file 2

Table S1. List of markers associated with crown rust resistance based on genome wide association studies in complete population and IBERS material. (XLSX 14 kb)

Additional file 3

Table S2. List of genomic scaffolds where all the significant markers from genome wide association studies were located. Scaffolds were placed onto linkage group with the aid of Genome Zipper [59]. (XLSX 11 kb)

Additional file 4

Table S3. List of predicted proteins on the genomic scaffolds. Markers located on these scaffolds were associated with crown rust resistance. BLAST was done on the predicted protein sequences using PLAZA [62] to obtain the gene function. (XLSX 17 kb)

Additional file 5

Figure S2. Among-and-within-full-sib-family selection that incorporates an inexpensive genotyping assay to implement within-family selection using a high selection intensity. (PDF 227 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Arojju, S., Conaghan, P., Barth, S. et al. Genomic prediction of crown rust resistance in Lolium perenne. BMC Genet 19, 35 (2018).

Download citation


  • Genomic selection
  • Crown rust
  • Perennial ryegrass
  • Genetic relationship
  • GWAS