Next Article in Journal
Pollination Parameter Optimization and Field Verification of UAV-Based Pollination of ‘Kuerle Xiangli’
Next Article in Special Issue
Association Mapping for Quantitative Trait Loci Controlling Superoxide Dismutase, Flavonoids, Anthocyanins, Carotenoids, γ-Oryzanol and Antioxidant Activity in Rice
Previous Article in Journal
A New Approach on Making European Agriculture More Efficient under Uncertainty Conditions
Previous Article in Special Issue
Alternative Splicing (AS) Dynamics in Dwarf Soybean Derived from Cross of Glycine max and Glycine soja
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Pedigree-Based BLUP and Phenotypic Mass Selection for Developing Elite Inbred Lines, Based on Field and Simulated Data

by
José Marcelo Soriano Viana
*,
Kaio Olimpio das Graças Dias
and
Jean Paulo Aparecido da Silva
Department of General Biology, Federal University of Viçosa, Viçosa 36570-900, MG, Brazil
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(10), 2560; https://doi.org/10.3390/agronomy12102560
Submission received: 6 September 2022 / Revised: 29 September 2022 / Accepted: 12 October 2022 / Published: 19 October 2022
(This article belongs to the Special Issue Genetics, Genomics and Breeding of Cereals and Grain Legumes)

Abstract

:
Pedigree-based best linear unbiased prediction (pBLUP) is an effective method for genetic evaluation. The objective of this study was to assess the efficacy of pBLUP to develop superior inbred lines, using field and simulated data. This study was based on the data of 30,986 plants from six S0 generations and eight inbred progeny trials (S1–S4), performed over eight years. We measured the grain yield and the expansion volume (EV). We also simulated the EV and grain yield of 5000 plants (S0–S4) in two in silico populations, assuming 300 and 400 genes distributed across 10 chromosomes of 200 cM, respectively, and 10% of pleiotropic genes. We selected selfed plants based on the predicted additive value for EV by fitting the individual model in inbred population. The use of pBLUP provided total genetic gains in EV in the range of 1 to 45%, inversely proportional to the level of improvement of the population, and indirect changes in grain yield in the range of −17 to 3%. Only from the analysis of the in silico populations, assuming selection based on the true additive value and genetic gain computed from the genotypic values, there was evidence that pBLUP is superior to phenotypic mass selection.

1. Introduction

The best linear unbiased prediction (BLUP) is a statistical approach that allows the prediction of a non-observable random variable, as the additive genetic value [1]. Its features make BLUP suitable for genetic evaluation under a wide variety of situations [2]. Advances in computing algorithms have led animal breeders to accept BLUP as a standard procedure for genetic evaluation in animal breeding, based on frequentist and Bayesian approaches [3]. The initial and the most significant application of BLUP in plant breeding is the prediction of non-assessed single crosses [4]. It is impressive that the initial proposition was based on relationship information from molecular markers. In our opinion, the most significant aspect of BLUP in animal and plant breeding, as well as in human genetics, is the prediction of the additive or genotypic values of individuals that cannot be measured for a sex-limited trait or that are too young to be measured for a trait, or of a subset of individuals that were not measured for a trait because quantification is too expansive. However, this remarkable feature depends on a key basic aspect of BLUP: relationship information [5].
The most important models for genetic assessment in animal breeding, such as the animal, sire, and gametic models, among others, have application in plant breeding, but denominated individual, half-sib progeny, and full-sib progeny models, respectively [6,7]. The individual model is used for the prediction of individual additive value using individual phenotype. The half-sib progeny model aims to predict the common parent additive value using the average phenotypic value of their progeny. The full-sib progeny model is used for the prediction of the parents’ additive values using the average phenotypic value of their progeny. Since 2001, when Meuwissen et al. [8] proposed the use of a high number of single nucleotide polymorphisms (SNPs) and a limited number of phenotypic records to predict the genetic merits of a sample of selection candidates (genomic selection), BLUP allowed animal and plant breeders to reach a higher level of efficacy, permitting fast genetic gains under a decreasing cost, keeping comparable accuracies relative to previous more expensive procedures, as progeny test [9,10]. Note, however, that the only innovation in genomic BLUP (GBLUP) was the use of a genomic relationship matrix (observed) instead of a pedigree-based relationship matrix (expected) [11]. However, this was a very significant innovation. The use of SNPs for measuring relationship allows discrimination among equivalent relatives and prediction of non-relatives based on identity by state and linkage disequilibrium between SNPs and genes [12].
In recent years, few studies, especially on plant breeding, have published results from genetic assessment based exclusively on pedigree-based BLUP. In these few papers, pedigree-based BLUP was regularly compared with a least squares approach. In most papers using GBLUP, especially in animal breeding, this approach was commonly compared with pedigree-based BLUP. In a survey based on the results from 50 recently published papers (2018–2022), in approximately 50% of the studies the authors emphasized that GBLUP, especially the single-step procedure, provided superior prediction accuracy [13]. Velazco et al. [14] observed in sorghum that a combined pedigree-genomic relationship matrix provided the largest increases in predictive ability and the largest reductions in prediction bias. Comparing pedigree-based and genomic BLUP, Suontama et al. [15] concluded that a significant increment in the breeding value prediction accuracy and genetic gains was possible by implementing genomic selection for wood properties in E. nitens. On the contrary, in the study of Kainer et al. [16], pedigree-based BLUP performed best across eight traits of E. polybractea. Using a simulated dairy cattle population, Seno et al. [17] did not observe significant differences between pedigree-based and genomic BLUP.
Thus, pedigree-based BLUP is an effective method for genetic evaluation in plant and animal breeding. To our knowledge, most studies using pedigree-based BLUP in animal and plant populations focused on the estimation of genetic parameters (components of genetic variances, heritability, and genetic/genotypic correlation) and assessment of the selection process efficacy, from the estimation of prediction accuracy and genetic gains. In plant breeding, these investigations generally involved a single cycle of selection. In animal breeding, these studies involved a huge amount of data from several generations of a single population. By using field and simulated datasets over four generations of selfing, our objective was to assess the efficacy of pedigree-based BLUP to develop superior inbred lines.

2. Materials and Methods

2.1. Populations and Trials

This study was based on data from six S0 generations and eight inbred progeny trials, performed over the growing seasons 2013/2014 (S1 progeny test) to 2020/2021 (S4 progeny test) in an experimental area at Federal University of Viçosa, Minas Gerais, Brazil. Due to a low stand in the flowering time or a reduced number of selfed plants, the S1 and S4 progeny trials were replicated twice (2013/2014 and 2014/2015, and 2019/2020 and 2020/2021). Note that the progeny were replicated, not the plants. Distinct groups of S2 and S3 progeny were assessed over two years (2015/2016 and 2016/2017, and 2017/2018 and 2018/2019, respectively). The S1 progeny were obtained over the growing seasons 2011/2012 and 2012/2013 from a tropical synthetic and five temperate populations. The synthetic was generated from the recombination of 30 inbred lines from the tropical populations Viçosa and Beija-Flor. The other populations (UFV-1 to 5) were obtained from recombining the single crosses P622, P625, P802, AP2501, and AP4502, developed by the Agricultural Alumni Seed, Indiana, USA. All trials were designed as completely randomized with 40–70 replicated inbred lines (twice). Each plot corresponded to a 5 m row with an ideal stand of 25 plants, spaced by 0.2 m. The space between plots was 0.9 m. In each plot, at least one plant was selfed, if available. The number of progeny assessed per population ranged from 32 to 840. The number of selfed plants per progeny ranged between 1 and 47 (four, on average). For each selfed plant (S0 to S4), we measured the grain weight of the selfed ear (g) and the expansion volume (EV; mL/g), a measure of popcorn quality. To assess EV, we used a 27 L microwave oven (900 W) and samples of 10 g. The processed dataset includes information from 30,986 plants with 145 and 1689 missing data for grain weight and EV, respectively.
Regarding the selection process, over the years we did not process the data and the criterion for selecting a selfed plant was an EV higher than or equal to 40 mL/g. Because we accomplished mass selection based on EV, it was possible to compare pedigree-based BLUP and phenotypic mass selection. To carry out a more accurate assessment on the relative efficacy of pedigree-based BLUP, compared with phenotypic mass selection, we included two other populations, generated using simulation. These populations will be referred to as in silico. We used the software REALbreeding (available upon request), which has been developed in Xojo (https://www.xojo.com/, accessed on 5 September 2022). REALbreeding has been used to provide simulated data in investigations in the areas of genomic selection, GWAS, QTL mapping, linkage disequilibrium, population structure, heterotic grouping/genetic diversity, and quantitative genetics. The current version allows the inclusion of digenic epistasis, genotype x environment interaction, and multiple traits, including pleiotropy. Details on the theoretical background and the stages for simulating individual genotypes for markers and genes and phenotypes can be found in a series of papers published by the main author (the software’s developer) since 2013, including Viana and Garcia [18]. REALbreeding can also be used for research on human genetics, animal genetics and breeding, population genetics, and evolution.
We simulated EV (mL/g) and grain yield (g/0.18 m2) assuming 300 and 400 genes distributed across 10 chromosomes of 200 cM, respectively, 10% of pleiotropic genes (30 genes), and no epistasis. The minimum and maximum genotypic values for EV and grain yield were 60 and 15, and 90 and 20, respectively. The minimum and maximum phenotypic values were 70 and 5, and 100 and 10. The non-inbred populations were defined assuming 0.4 and 0.6 as the average frequency of the favorable genes for EV in the in silico tropical (low EV) and in silico temperate (high EV). Concerning grain yield, for both in silico populations, we assumed 0.3 as the average frequency of the favorable genes. We assumed bidirectional dominance (average degree of dominance: 0.01) for EV and positive dominance (average degree of dominance: 0.57) for grain yield. The sample size from S0 to S4 was 1000 plants, assuming 0 to 10 selfed plants per progeny, with an average of four. We maximized the linkage disequilibrium (LD) in the populations and the number of positive (in silico tropical) and negative (in silico temperate) LD values (for D, r2, and D’). This led to significant positive and negative genotypic correlations between the traits. The heritabilities at the individual level for EV and grain yield were 30 and 20%, respectively.

2.2. Efficacy of Pedigree-Based BLUP and Statistical Analysis

Assume a single biallelic gene (A/a) and a non-inbred population in Hardy–Weinberg equilibrium. The genotype probabilities are p2 AA: 2pq Aa: q2 aa, where p and q are the allele frequencies. The population mean is M = m + (p − q)a + 2pqd. The parameters m, a, and d are the mean of the genotypic values of the homozygotes, the deviation between the genotypic value of the homozygote of greater expression and m, and the dominance deviation, respectively [19]. After n − 1 generations of selfing under selection, the genotype probabilities in generation n can be defined as p2 + pqF + ΔAA AA: 2pq(1 − F) + ΔAa Aa: q2 + pqF + Δaa aa, where Δ is the change in the genotype probability due to selection and F is the inbreeding coefficient. Note that ΔAA + ΔAa + Δaa = 0. Assuming no selection (ΔAA = ΔAa = Δaa = 0), the inbred population mean is MF = M − 2Fpqd, where −2Fpqd is the change in the non-inbred population mean due to inbreeding. If the selection changed the allelic frequencies, the mean of the generation n is MFs = MF + ΔAA(m + a) + ΔAa(m + d) + Δaa(m − a) = MF + [(ΔAA − Δaa)a + ΔAad] = MF + Δg, where Δg = MFs − MF is the true genetic gain. Using phenotypic values, we have a predicted genetic gain of Dg = (avPFs − avPF)·h2, where avP stand for the average phenotypic value and h2 is the heritability. This is a standard function for the estimated genetic gain: selection differential times heritability.
We adopted a selection process based on the predicted additive value by fitting the individual model in inbred population [20]. The percentages of selected plants by phenotypic mass selection were approximately 60, 50, 60, and 90, from S1 to S4. However, in this study, we assumed 50, 45 and 40% of selected S1, S2, and S3 plants for both selection processes (Figure 1). Because there was selection based on EV, we computed the increment in the genetic gain due to pedigree-based BLUP, relative to the gain attributable to phenotypic mass selection, given by iDg = (avPFsb − avPFsp)·h2, where the first mean refers to the individuals selected by pedigree-based BLUP and the second mean refers to the individuals selected by phenotypic mass selection. Because REALbreeding provides the parametric population means and genetic variances and covariances, and the true additive and dominance genetic values per generation, in the in silico populations the relative efficacy of pedigree-based BLUP selection was based on the parametric genotypic values (iΔg = avGsb − avGsp), where avGsb and avGsp are the average genotypic values of the selected plants under pedigree-based BLUP and phenotypic mass selection, respectively.
The covariance between relatives in inbred populations, assuming linkage equilibrium and no epistasis, was derived by Cockerham [21]. Assuming biallelic genes, the covariance between relatives with self-fertilization depends on the additive ( σ A 2 ) and dominance ( σ D 2 ) variances, the covariance between the a and d deviations (D1), and the variance of the dominance deviations (D2). Ignoring the components D1 and D2, the diagonal and off-diagonal elements of the additive and dominance relationship matrices were computed from the covariances c g g g = ( 1 + F g ) σ A 2 + ( 1 F g 2 ) σ D 2 and c t g g = ( 1 + F t ) σ A 2 + [ ( 1 + F t ) ( 1 F g ) ( 1 F g ) / ( 1 F t ) ] σ D 2 , respectively, where t, g, and g’ are the generations. We used t as the first common ancestor (S0 plant). Thus, for example, c 011 is the covariance between two S1 plants derived from the same S0 plant and c 012 is the covariance between a S1 and a S2 plant derived from the same S0 plant. Although the elements of the additive and dominance relationship matrices are easy to compute, since they are a function of the inbreeding coefficient, the process is not available in known R packages such as AGHmatrix, ggroups, pedigree, and pedigreemm. We used a REALbreeding tool (Relationship) to compute the additive and dominance relationship matrices.
We performed the statistical analyses fitting the linear mixed model y   =   μ 1 +   X β   +   Z u   + ϵ using ASReml-R v.4 [22], where y is the vector of phenotypes, μ is the general mean, β is the matrix of fixed effects of years and generations, u is the vector of additive values within generations, with u   N ( 0 ,   Σ g ) , and ϵ is the error vector, with ϵ   N ( 0 ,   I σ ϵ 2 ) . Σ g is a diagonal matrix with an additive variance for each generation. 1 is a vector of ones and X and Z are incidence matrices. We estimated the components of variance using the Residual Maximum Likelihood (REML) [23].

3. Results

The in silico tropical population shows similar means for EV and grain yield with Synthetic UFV. The in silico temperate population shows an EV comparable to the temperate populations. However, both in silico populations show much lower phenotypic variability relative to all other populations (Table 1). In the in silico populations, because of bidirectional dominance, the inbreeding had no significant effect on the EV mean of the inbred generations. Due to positive dominance, the grain yield decreased by 16% with inbreeding, on average. The inbreeding increased the genotypic variability for both traits. Note also that the inbreeding did not significantly affect the genotypic correlation between the traits. The values were significantly negative and positive for the temperate and tropical in silico populations, respectively. These results are consistent with no selection.
In regard to the popcorn populations, there is some evidence that the phenotypic mass selection improved EV, with a significant EV increase in Synthetic UFV and a substantial decrease in the phenotypic variability over the selfing generations (Table 1). This is also supported by the considerable decrease in the estimates of the additive variances for EV over the selfing generations. However, there is, in general, no evidence of significant indirect changes in grain yield. The results for grain yield are, in general, those expected with inbreeding. However, changes occurred in the estimated genetic correlations over the selfed generations. Because only BLUP provided the estimates of the additive variances for EV and grain yield, the accuracies for phenotypic mass selection and pedigree-based BLUP are very similar (Table 2). Thus, the experimental evidence is that both processes are equally efficient. Based on the same two accuracies, this is also true for the in silico populations. Comparing the accuracies in the tropical and temperate in silico populations, the general conclusion is that both selection processes have similar efficiencies.
In complete agreement with the accuracy analysis, the estimated increments in the genetic gain for EV with pedigree-based BLUP selection, relative to phenotypic mass selection, indicate a similar efficacy for the popcorn populations. The total increments were close to zero, in the range of −0.14 to 0.10 mL/g. However, slightly higher efficacy occurred in the in silico populations. The total increments were 0.30 and 0.50 mL/g (Table 2). The true total increments were 1.26 and 1.70 mL/g for the tropical and temperate populations, respectively. These values are associated with correlations between the predicted and parametric additive values of approximately 0.7, 0.7–0.8, and 0.8 in the generations S1, S2, and S3, respectively. However, selecting based on the parametric additive values instead of the predicted additive values, the true total increments were 2.54 and 3.53 mL/g, respectively.
In regard to the indirect increments for grain yield, we observed small values in the popcorn populations, in the range of −0.36 to 0.44 g/0.18 m2 (Table 2). In the in silico populations, the true indirect increments were 0.47 and −0.29 g/0.18 m2 for the tropical and temperate populations, respectively. The true changes due to inbreeding were 1.59 and −2.69 g/0.18 m2, respectively. To confirm a higher efficacy of pedigree-based BLUP, relative to phenotypic mass selection, observed from a single simulation, we also computed the true increments in other nine simulations of the in silico temperate population, assuming selection based on the predicted and true additive values (Table 3). The average direct total increments, 0.34 and 1.74 mL/g for selection based on the predicted and true additive values, confirm that pedigree-based BLUP is superior to phenotypic mass selection. The average total indirect increments in grain yield are also negative with small magnitude, indicating that pedigree-based BLUP provides slightly different indirect changes in grain yield, relative to phenotypic mass selection. It is important to emphasize that the average true direct increments with pedigree-based BLUP selection in the S1, S2, and S3 generations are proportional to the correlation between the predicted and parametric additive values. In S3, the accuracy was 26% lower than the maximum value and the true direct increment was 54% lower relative to the maximum value.
Only in Synthetic UFV we observed a significant estimated genetic gain for EV, with pedigree-based BLUP selection (approximately 14 mL/g) (Figure 2). The genetic gains in the other populations ranged from approximately 0.5 to 4.6 mL/g. This cannot be solely due to the reduced S4 size (76 plants), since the number of S4 plants in the in silico temperate population was only 53. Remember that the proportion of selected plants was kept constant in each population. Note that, as theoretically expected, higher gains were observed in the tropical populations (less improved). As expected, the selection process decreased the variability in the populations in similar proportions, between −58 and −92%. The selection for EV in Synthetic UFV also provided the highest indirect change in grain yield, of approximately −18 g/0.18 m2. This is surprising since the genotypic correlations ranged between −0.1 and 0.2, approximately. In the in silico populations, with significant positive and negative genotypic correlations, the indirect changes were also negative. The change in the variability for grain yield was not consistent, varying from −67 to 28%.

4. Discussion

In a vast amount of studies on animal and plant breeding, distinct statistical approaches were compared, such as least squares (no relationship information), pedigree-based BLUP, and GBLUP, especially single-step GBLUP in recent years. Based on the results, it cannot be definitely stated that pedigree-based BLUP is superior to least squares and that GBLUP is superior to pedigree-based BLUP, in any situation [25,26,27,28]. However, taking into account the BLUP principles, the theoretical advantages of including relationship information, and the huge amount of favorable field results from genetic assessment in plant and animal breeding [5], there is no reason to not use BLUP for genetic evaluation and for the prediction of complex traits. Using a simulated population, Seno, Guidolin, Aspilcueta-Borquis, do Nascimento, da Silva, de Oliveira and Munari [17] investigated the efficacy of phenotypic selection, pedigree-based BLUP, and GBLUP over 25 years of selection. They observed equivalence for the two BLUP approaches concerning the rates of gain and average inbreeding coefficient. However, both methods yielded superior genetic gains relative to phenotypic selection.
When pedigree-based BLUP was applied to the recurrent genetic assessment of inbred lines, it was not possible to fit the additive-dominance model. This occurred not because of a problem with the additive and dominance relationship matrices (both well-conditioned) but due to singularities in the average information matrix. Theoretically, fitting the complete model would be advantageous since in inbred populations additive and dominance genetic values are correlated [19]. However, as shown from the analysis of the simulated data, even ignoring dominance, the correlation between the predicted and true additive values ranged from 0.5 (S0) to 0.8 (S4), proportional to the amount of relationship information. Velazco, Malosetti, Hunt, Mace, Jordan and van Eeuwijk [14] observed that the impact of including genealogy information to improve predictions was stronger for the lower heritability traits, such as grain yield and stay green. The analysis of the simulated data also showed comparable estimated and true values for the genotypic correlation between EV and grain yield.
For the popcorn and the in silico populations, we generally observed equivalence between the prediction accuracies based on estimated heritability and on prediction error variance [24]. However, only in the third selection cycle did we observe a positive correlation between the prediction accuracy and the genetic gain with pedigree-based BLUP (0.48). This does not mean that estimating accuracy is useless, since it is an appropriate measure for assessing and comparing selection processes and statistical approaches. Viana et al. [29] observed that the efficacy for identifying the superior 9% untested single crosses achieved 0.71 under a prediction accuracy of 0.92. In the study of Viana et al. [30], the efficacy of identification of the superior 6–9% of the untested testcrosses reached 0.76 when the prediction accuracy was 0.94. Note that, because recurrent selection provides theoretically lower genetic gains over selection cycles, only with a decrease in the prediction accuracy due to a decline in the genetic variability, under similar environmental conditions, will the correlation be positive. In the in silico populations, the prediction accuracies showed comparable or higher values for the correlation between the predicted and parametric additive values over the selection cycles, associated with decreasing true genetic gains. The genetic gains due to pedigree-based BLUP in S1 to S3 were 1.8, 1.4, and 1.3 mL/g for the tropical population and 1.5, 1.3, and 0.9 mL/g for the temperate population. Using simulation, Jibrila et al. [31] observed a slight decrease in the genetic gains in subsequent generations by increasing the intensity of preselection. They assumed genetic gain as the difference between the average true breeding values of the selected individuals in two subsequent generations.
Probably due to sampling, the analysis of the increments in the genetic gain with pedigree-based BLUP using the predicted additive value, relative to phenotypic mass selection, showed that the procedures were equivalent. However, the analysis of the increments in the non-replicated tropical and the replicated temperate in silico populations, assuming selection based on the predicted and true additive values, clearly showed that pedigree-based BLUP was superior to mass phenotypic selection. The breeders should not be disappointed with the magnitude of the genetic gains with pedigree-based BLUP because they are proportional to the percentage of selected plants. Note that even assuming 50, 45, and 40% of selected plants, there was a significant decrease in the genetic variability for EV in the populations. Using pedigree-based BLUP, Suontama, Klapste, Telfer, Graham, Stovold, Low, McKinley and Dungey [15] observed genetic gains for several growth, form, and wood quality traits across two E. nitens orchards. For most traits, the genetic gains were comparable (at least 80%) to those obtained with GBLUP.
Concerning the indirect changes in grain yield, the analysis of the alterations in the in silico population means due to inbreeding under no selection and inbreeding under indirect selection shows that inbreeding and indirect selection induced comparable changes. Because negative correlation in the in silico temperate population, the decrease in grain yield due to inbreeding and indirect selection was 1.9 times greater than the decrease due to inbreeding. On the contrary, because of the positive correlation in the in silico tropical population, the decrease in grain yield due to inbreeding was 2.1 times greater than the decrease due to inbreeding and indirect selection. El-Attrouny et al. [32] investigated the effects of pedigree-based selection for Japanese quail bodyweight at four weeks on bodyweight and bodyweight gain across four generations and observed significant indirect genetic gains. Because of moderate heritabilities, the strong positive genotypic correlation between spawn weight and number of eggs in the spawn and individual egg size traits, and the lack of correlation between the number of eggs in the spawn with egg size traits and female body weight, D’Ambrosio et al. [33] concluded that selection for growth will not induce indirect improvements in female reproduction traits. In the study of Cobo et al. [34] with rainbow trout, the genetic progress curves of the index MAT and litter weight (LW) indicated that the selection based on MAT gave a positive correlated response on LW. This occurred due to a high correlation between the traits (0.85).
Concluding, the use of pedigree-based BLUP in eight populations, over four generations of selfing, provided total genetic gains in EV in the range of 1 to 45%, inversely proportional to the level of improvement of the reference population, and indirect changes in grain yield in the range of −17 to 3%, values that correspond to approximately half of the values assuming inbreeding and indirect selection. The analysis of the in silico populations, assuming selection based on the true additive value and genetic gain computed from the genotypic values, proved that pedigree-based BLUP is superior to the phenotypic mass selection. The average direct total increments were 0.34 and 1.74 mL/g for selection based on the predicted and true additive values, respectively. The equivalence between both selection processes observed in the popcorn populations are due to sampling, as observed in the replications of the in silico temperate population. Based on our results, then, we strongly recommend the use of pedigree-based BLUP for genetic evaluation of inbred progeny.

Author Contributions

Conceptualization, J.M.S.V.; methodology, J.M.S.V. and K.O.d.G.D.; software, J.M.S.V.; formal analysis, J.P.A.d.S.; writing—original draft preparation, J.M.S.V.; writing—review and editing, K.O.d.G.D. and J.P.A.d.S.; funding acquisition, J.M.S.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The dataset is available at https://doi.org/10.6084/m9.figshare.19928330 (accessed on 5 September 2022).

Acknowledgments

We thank the National Council for Scientific and Technological Development (CNPq), the Brazilian Federal Agency for Support and Evaluation of Graduate Education (Capes), and the Foundation for Research Support of Minas Gerais State (Fapemig) for financial support.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Henderson, C.R. Sire evaluation and genetic trends. J. Anim. Sci. 1973, 1973, 10–41. [Google Scholar] [CrossRef]
  2. Henderson, C.R. General flexibility of linear model techniques for sire evaluation. J. Dairy Sci. 1974, 57, 963–972. [Google Scholar] [CrossRef]
  3. Blasco, A. The Bayesian controversy in animal breeding. J. Anim. Sci. 2001, 79, 2023–2046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bernardo, R. Prediction of Maize Single-Cross Performance Using RFLPs and Information from Related Hybrids. Crop Sci. 1994, 34, 20–25. [Google Scholar] [CrossRef]
  5. Gianola, D.; Rosa, G.J. One hundred years of statistical developments in animal breeding. Annu. Rev. Anim. Biosci. 2015, 3, 19–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Viana, J.M.S.; Faria, V.R.; Fonseca e Silva, F.; Vilela de Resende, M.D. Best Linear Unbiased Prediction and Family Selection in Crop Species. Crop Sci. 2011, 51, 2371–2381. [Google Scholar] [CrossRef]
  7. Viana, J.M.S.; de Almeida, I.F.; Vilela de Resende, M.D.; Faria, V.R.; Fonseca e Silva, F. BLUP for genetic evaluation of plants in non-inbred families of annual crops. Euphytica 2010, 174, 31–39. [Google Scholar] [CrossRef]
  8. Meuwissen, T.H.E.; Hayes, B.J.; Goddard, M.E. Prediction of Total Genetic Value Using Genome-Wide Dense Marker Maps. Genetics 2001, 157, 1819–1829. [Google Scholar] [CrossRef]
  9. Van Eenennaam, A.L.; Weigel, K.A.; Young, A.E.; Cleveland, M.A.; Dekkers, J.C.M. Applied Animal Genomics: Results from the Field. Annu. Rev. Anim. Biosci. 2014, 2, 105–139. [Google Scholar] [CrossRef] [Green Version]
  10. Meuwissen, T.; Hayes, B.; Goddard, M. Accelerating Improvement of Livestock with Genomic Selection. Annu. Rev. Anim. Biosci. 2013, 1, 221–237. [Google Scholar] [CrossRef]
  11. Gianola, D.; Cecchinato, A.; Naya, H.; Schon, C.C. Prediction of Complex Traits: Robust Alternatives to Best Linear Unbiased Prediction. Front. Genet. 2018, 9, 195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Vela-Avitua, S.; Meuwissen, T.H.E.; Luan, T.; Odegard, J. Accuracy of genomic selection for a sib-evaluated trait using identity-by-state and identity-by-descent relationships. Genet. Sel. Evol. 2015, 47, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Scholtens, M.; Lopez-Villalobos, N.; Lehnert, K.; Snell, R.; Garrick, D.; Blair, H.T. Advantage of including Genomic Information to Predict Breeding Values for Lactation Yields of Milk, Fat, and Protein or Somatic Cell Score in a New Zealand Dairy Goat Herd. Animals 2021, 11, 24. [Google Scholar] [CrossRef] [PubMed]
  14. Velazco, J.G.; Malosetti, M.; Hunt, C.H.; Mace, E.S.; Jordan, D.R.; van Eeuwijk, F.A. Combining pedigree and genomic information to improve prediction quality: An example in sorghum. Theor. Appl. Genet. 2019, 132, 2055–2067. [Google Scholar] [CrossRef] [Green Version]
  15. Suontama, M.; Klapste, J.; Telfer, E.; Graham, N.; Stovold, T.; Low, C.; McKinley, R.; Dungey, H. Efficiency of genomic prediction across two Eucalyptus nitens seed orchards with different selection histories. Heredity 2019, 122, 370–379. [Google Scholar] [CrossRef] [Green Version]
  16. Kainer, D.; Stone, E.A.; Padovan, A.; Foley, W.J.; Kulheim, C. Accuracy of Genomic Prediction for Foliar Terpene Traits in Eucalyptus polybractea. G3-Genes Genomes Genet. 2018, 8, 2573–2583. [Google Scholar] [CrossRef] [Green Version]
  17. Seno, L.D.; Guidolin, D.G.F.; Aspilcueta-Borquis, R.R.; do Nascimento, G.B.; da Silva, T.B.R.; de Oliveira, H.N.; Munari, D.P. Genomic selection in dairy cattle simulated populations. J. Dairy Res. 2018, 85, 125–132. [Google Scholar] [CrossRef] [Green Version]
  18. Viana, J.M.S.; Garcia, A.A.F. Significance of linkage disequilibrium and epistasis on genetic variances in noninbred and inbred populations. BMC Genom. 2022, 23, 286. [Google Scholar] [CrossRef]
  19. Falconer, D.S.; Mackay, T.F.C. Introduction to Quantitative Genetics, 4th ed.; Longman: London, UK, 1996. [Google Scholar]
  20. Viana, J.M.S.; Faria, V.R.; Fonseca e Silva, F.; Vilela de Resende, M.D. Combined selection of progeny in crop breeding using best linear unbiased prediction. Can. J. Plant Sci. 2012, 92, 553–562. [Google Scholar] [CrossRef] [Green Version]
  21. Cockerham, C.C. Covariances of relatives from self-fertilization. Crop Sci. 1983, 23, 1177–1180. [Google Scholar] [CrossRef]
  22. Butler, D.G.; Cullis, B.R.; Gilmour, A.R.; Gogel, B.G.; Thompson, R. ASReml-R Reference Manual Version 4; VSN International Ltd.: Hemel Hempstead, UK, 2017. [Google Scholar]
  23. Patterson, H.D.; Thompson, R. Recovery of inter-block information when block sizes are unequal. Biometrika 1971, 58, 545–554. [Google Scholar] [CrossRef]
  24. Mrode, R.A. Linear Models for the Prediction of Animal Breeding Values, 2nd ed.; CABI Publishing: Wallingford, UK, 2005. [Google Scholar]
  25. Mehrban, H.; Naserkheil, M.; Lee, D.; Ibanez-Escriche, N. Multi-Trait Single-Step GBLUP Improves Accuracy of Genomic Prediction for Carcass Traits Using Yearling Weight and Ultrasound Traits in Hanwoo. Front. Genet. 2021, 12, 692356. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, J.; Wang, J.; Li, Q.H.; Wang, Q.; Wen, J.; Zhao, G.P. Comparison of the Efficiency of BLUP and GBLUP in Genomic Prediction of Immune Traits in Chickens. Animals 2020, 10, 419. [Google Scholar] [CrossRef] [PubMed]
  27. Aguilar, I.; Fernandez, E.N.; Blasco, A.; Ravagnolo, O.; Legarra, A. Effects of ignoring inbreeding in model-based accuracy for BLUP and SSGBLUP. J. Anim. Breed. Genet. 2020, 137, 356–364. [Google Scholar] [CrossRef]
  28. Dunne, F.L.; Kelleher, M.M.; Walsh, S.W.; Berry, D.P. Characterization of best linear unbiased estimates generated from national genetic evaluations of reproductive performance, survival, and milk yield in dairy cows. J. Dairy Sci. 2018, 101, 7625–7637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Viana, J.M.S.; Pereira, H.D.; Mundim, G.B.; Piepho, H.P.; Silva, F.F.E. Efficiency of genomic prediction of non-assessed single crosses. Heredity 2018, 120, 283–295. [Google Scholar] [CrossRef] [Green Version]
  30. Viana, J.M.S.; Pereira, H.D.; Piepho, H.P.; Silva, F.F.E. Efficiency of Genomic Prediction of Nonassessed Testcrosses. Crop Sci. 2019, 59, 2020–2027. [Google Scholar] [CrossRef]
  31. Jibrila, I.; ten Napel, J.; Vandenplas, J.; Veerkamp, R.F.; Calus, M.P.L. Investigating the impact of preselection on subsequent single-step genomic BLUP evaluation of preselected animals. Genet. Sel. Evol. 2020, 52, 42. [Google Scholar] [CrossRef]
  32. El-Attrouny, M.M.; Manaa, E.A.; Ramadan, S.I. Genetic evaluation and selection correlated response of growth traits in Japanese quail. S. Afr. J. Anim. Sci. 2020, 50, 325–333. [Google Scholar] [CrossRef]
  33. D’Ambrosio, J.; Morvezen, R.; Brard-Fudulea, S.; Bestin, A.; Perez, A.A.; Guemen, D.; Poncet, C.; Haffray, P.; Dupont-Nivet, M.; Phocas, F. Genetic architecture and genomic selection of female reproduction traits in rainbow trout. BMC Genom. 2020, 21, 558. [Google Scholar] [CrossRef]
  34. Cobo, E.; Raoul, J.; Bodin, L. Genetic parameters of litter weight, an alternative criterion to prolificacy and pre-weaning weight for selection of French meat sheep. Livest. Sci. 2021, 250, 104596. [Google Scholar] [CrossRef]
Figure 1. Sample size and number of selected plants in S1, S2, and S3 generations, in each popcorn and in silico populations, under pedigree-based BLUP and phenotypic mass selection.
Figure 1. Sample size and number of selected plants in S1, S2, and S3 generations, in each popcorn and in silico populations, under pedigree-based BLUP and phenotypic mass selection.
Agronomy 12 02560 g001
Figure 2. Changes in the populations means (a,b) and variances (c,d) over generations (S0 = 1) with pedigree-based BLUP, using the predicted or true additive (for the in silico populations) values in the generations S1, S2, and S3, relative to EV ((a,c); mL/g) and grain yield ((b,d); g/0.18 m2).
Figure 2. Changes in the populations means (a,b) and variances (c,d) over generations (S0 = 1) with pedigree-based BLUP, using the predicted or true additive (for the in silico populations) values in the generations S1, S2, and S3, relative to EV ((a,c); mL/g) and grain yield ((b,d); g/0.18 m2).
Agronomy 12 02560 g002
Table 1. Number of progeny, sample size, mean (M/MF), phenotypic ( σ P 2 ), additive ( σ A 2 ), dominance ( σ D 2 ), and error ( σ E 2 ) variances, covariance between additive and dominance values ( σ A , D ), and genotypic correlation ( σ g ), per generation, in the popcorn (estimates) and in silico (estimates in the 2nd row) populations, relative to EV (mL/g) and grain yield (g/0.18 m2).
Table 1. Number of progeny, sample size, mean (M/MF), phenotypic ( σ P 2 ), additive ( σ A 2 ), dominance ( σ D 2 ), and error ( σ E 2 ) variances, covariance between additive and dominance values ( σ A , D ), and genotypic correlation ( σ g ), per generation, in the popcorn (estimates) and in silico (estimates in the 2nd row) populations, relative to EV (mL/g) and grain yield (g/0.18 m2).
Pop.Gen.Prog.SizeEVGrain Yield σ g
M/MF σ P 2 σ A 2 σ D 2 σ A , D σ E 2 M/MF σ P 2 σ A 2 σ D 2 σ A , D σ E 2
Synthetic S0-41730.860.144.25--11.3452.9309.5155.13--148.16−0.05
S13216541.079.397.42--11.3440.6228.250.93--148.16−0.21
S25929039.351.043.91--11.3432.3209.836.10--148.160.57
S37536039.023.810.82--11.3430.4201.339.29--148.16−0.11
S47636743.311.80.50--11.3436.9168.211.86--148.160.11
UFV-1S0-37840.876.232.11--17.9042.8195.252.92--141.26−0.25
S121580443.973.069.59--17.9039.3186.030.61--141.26−0.40
S2317113040.754.037.82--17.9035.9225.746.58--141.260.02
S3276113439.625.43.85--17.9029.8166.615.68--141.26−0.58
S4427198143.117.30.26--17.9036.2177.017.62--141.26−1.00
UFV-2S0-39142.074.534.96--20.2137.9215.844.68--142.84−0.17
S116849045.267.047.00--20.2137.6194.029.96--142.84−0.13
S216961741.248.224.72--20.2136.6219.630.72--142.84−0.16
S318370639.224.92.30--20.2129.5177.117.08--142.84−0.20
S4315146842.720.30.70--20.2136.7173.817.23--142.840.03
UFV-3S0-53042.048.426.06--16.8237.1145.90.00--146.75-
S131094644.759.846.83--16.8241.0214.842.22--146.75−0.23
S2357133040.455.936.47--16.8235.4225.728.09--146.750.13
S3337142139.327.77.18--16.8229.9167.114.45--146.75−0.44
S4502243342.917.10.80--16.8236.4183.519.15--146.75−0.02
UFV-4S0-25240.756.824.13--20.8437.5205.171.60--133.33−0.10
S110136943.572.253.12--20.8441.5206.347.47--133.33−0.03
S211651740.661.638.44--20.8435.5218.132.94--133.330.12
S314460339.030.35.32--20.8430.6170.017.74--133.33−0.29
S4218110742.620.20.50--20.8436.2154.512.18--133.33−0.41
UFV-5S0-75342.452.835.73--12.1140.6168.439.58--126.92−0.53
S1546203743.967.578.68--12.1136.6132.69.63--126.92−1.00
S2600191340.751.741.12--12.1135.1231.350.00--126.92−0.26
S3533216339.829.813.21--12.1128.8151.615.92--126.92−0.63
S4840391443.212.60.59--12.1137.0191.346.52--126.92−1.00
In silicoS0 100042.914.04.040.170.009.8448.226.54.400.300.0018.82−0.65
temperate 42.813.41.11--12.3248.423.20.09--23.16−1.00
S1247100042.518.45.380.16−0.0112.8744.930.26.010.54−0.2823.93−0.61
42.516.33.13--12.3244.516.743.44--23.16−0.74
S2250100042.320.05.910.12−0.0213.9643.231.86.690.44−0.4325.09−0.65
42.116.64.05--12.3243.316.323.60--23.16−0.60
S3220100042.220.66.130.09−0.0314.4042.432.27.000.35−0.5025.39−0.65
42.016.53.36--12.3242.117.375.16--23.16−0.59
S4243100042.120.96.230.08−0.0314.5842.032.47.140.30−0.5325.49−0.52
41.616.43.57--12.3242.116.884.00--23.16−0.36
In silicoS0-100033.914.84.310.120.0010.3349.232.26.080.360.0025.760.78
tropical 34.015.32.77--12.4949.032.92.04--30.871.00
S1243100033.819.35.710.11−0.0113.5145.541.68.190.64−0.2833.060.75
34.120.45.34--12.4945.941.06.28--30.871.00
S2227100033.721.06.250.08−0.0214.6843.643.99.060.49−0.4334.760.80
33.921.04.87--12.4943.944.98.29--30.870.75
S3237100033.921.76.480.06−0.0215.1542.744.59.430.37−0.5035.200.81
33.922.55.15--12.4942.643.36.74--30.870.67
S4217100033.622.06.580.05−0.0215.3642.244.79.600.30−0.5335.340.82
34.023.75.45--12.4942.942.46.11--30.870.62
Table 2. Accuracies 1 of phenotypic mass selection and pedigree-based BLUP based on EV, percent of selected plants (%S), and increments of direct (iDgd; mL/g) and indirect (iDgi; g/plant) genetic gains per generation with pedigree-based BLUP, relative to phenotypic mass selection, in the popcorn and in silico populations.
Table 2. Accuracies 1 of phenotypic mass selection and pedigree-based BLUP based on EV, percent of selected plants (%S), and increments of direct (iDgd; mL/g) and indirect (iDgi; g/plant) genetic gains per generation with pedigree-based BLUP, relative to phenotypic mass selection, in the popcorn and in silico populations.
Pop.Gen.Ac1Ac2Ac3Ac4Ac5%SiDgdiDgi
SyntheticS0-0.89--0.880--
S1-0.95--0.92500.000.00
S2-0.89--0.8545−0.580.16
S3-0.70--0.55400.45−0.17
S4-0.20--- 2---
UFV-1S0-0.80--- 20--
S1-0.89--0.62500.13−0.08
S2-0.82--- 245−0.20−0.38
S3-0.42--- 2400.060.10
S4-0.12--- 2---
UFV-2S0-0.80--0.750--
S1-0.84--0.8550−0.010.03
S2-0.74--0.6945−0.04−0.05
S3-0.32--- 2400.020.46
S4-0.18--- 2---
UFV-3S0-0.78--0.760--
S1-0.86--0.86500.00−0.01
S2-0.83--0.83450.05−0.01
S3-0.55--- 2400.000.18
S4-0.21--- 2---
UFV-4S0-0.73--0.680--
S1-0.85--0.8650−0.25−0.10
S2-0.80--0.80450.06−0.07
S3-0.45--- 2400.05−0.04
S4-0.15--- 2---
UFV-5S0-0.86--0.840--
S1-0.93--0.9450−0.010.00
S2-0.88--0.8745−0.040.00
S3-0.72--0.49400.15−0.01
S4-0.21--- 2---
In silicoS00.540.290.520.520.280--
temperateS10.540.450.520.670.42500.07−0.04
0.46 3−0.31 3
1.16 4−0.68 4
S20.540.500.490.720.45450.14−0.04
0.41 3−0.16 3
0.32 4−0.27 4
S30.540.460.510.780.34400.29−0.20
0.83 3−0.96 3
2.05 4−1.74 4
S40.550.470.480.780.34---
In silicoS00.540.430.560.560.420--
tropicalS10.540.550.530.710.58500.020.02
0.15 30.23 3
0.75 40.48 4
S20.550.530.540.790.49450.03−0.01
0.15 3−0.06 3
0.54 40.30 4
S30.550.540.560.820.48400.260.31
0.96 30.86 3
1.25 40.81 4
S40.550.550.570.830.48---
1 Ac1: root square of the parametric narrow sense heritability; Ac2: root square of the estimated narrow sense heritability; Ac3: correlation between the parametric additive value and the phenotypic value; Ac4: correlation between the predicted and parametric additive values; Ac5 = 1 ( P E V / σ ^ A 2 ) , where PEV is the prediction error variance [24]. 2 P E V / σ ^ A 2 > 1 . 3 Increment due to selection based on predicted additive value and gains computed using the genotypic value. 4 Increment due to selection based on true additive value and gains computed using the genotypic value.
Table 3. Average additive and error variances per generation for EV, percent of selected plants (%S), accuracy 1, and increments in the direct (mL/g) and indirect genetic gains (g/0.18 m2) with pedigree-based BLUP based on EV, relative to phenotypic mass selection, in the in silico temperate population (minimum and maximum values between brackets).
Table 3. Average additive and error variances per generation for EV, percent of selected plants (%S), accuracy 1, and increments in the direct (mL/g) and indirect genetic gains (g/0.18 m2) with pedigree-based BLUP based on EV, relative to phenotypic mass selection, in the in silico temperate population (minimum and maximum values between brackets).
Gen. σ A 2 σ E 2 Ac%SiDgdiDgi
S01.99 [0.83, 3.42]11.80 [10.71, 12.94]0.54 [0.52, 0.58]---
S14.29 [3.06, 5.37]11.80 [10.71, 12.94]0.71 [0.63, 0.75]50−0.48 [−0.93, 0.59] 20.25 [−0.37, 0.74]
0.08 [−0.28, 1.16] 3−0.03 [−0.68, 0.35]
S24.36 [3.16, 5.35]11.80 [10.71, 12.94]0.75 [0.71, 0.81]450.32 [−0.04, 1.03] 2−0.23 [−0.77, 0.11]
0.33 [0.06, 0.66] 3−0.28 [−0.71, 0.04]
S34.57 [3.15, 6.76]11.80 [10.71, 12.94]0.78 [0.72, 0.82]400.50 [−0.40, 0.93] 2−0.33 [−1.12, 0.59]
1.33 [0.80, 2.21] 3−0.76 [−1.74, 0.18]
S44.32 [3.03, 5.87]11.80 [10.71, 12.94]0.80 [0.73, 0.83]---
1 Ac: correlation between the predicted and parametric additive values. 2 Increment due to selection based on predicted additive value and gains computed using the genotypic value. 3 Increment due to selection based on true additive value and gains computed using the genotypic value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Viana, J.M.S.; Dias, K.O.d.G.; Silva, J.P.A.d. Comparative Analysis of Pedigree-Based BLUP and Phenotypic Mass Selection for Developing Elite Inbred Lines, Based on Field and Simulated Data. Agronomy 2022, 12, 2560. https://doi.org/10.3390/agronomy12102560

AMA Style

Viana JMS, Dias KOdG, Silva JPAd. Comparative Analysis of Pedigree-Based BLUP and Phenotypic Mass Selection for Developing Elite Inbred Lines, Based on Field and Simulated Data. Agronomy. 2022; 12(10):2560. https://doi.org/10.3390/agronomy12102560

Chicago/Turabian Style

Viana, José Marcelo Soriano, Kaio Olimpio das Graças Dias, and Jean Paulo Aparecido da Silva. 2022. "Comparative Analysis of Pedigree-Based BLUP and Phenotypic Mass Selection for Developing Elite Inbred Lines, Based on Field and Simulated Data" Agronomy 12, no. 10: 2560. https://doi.org/10.3390/agronomy12102560

APA Style

Viana, J. M. S., Dias, K. O. d. G., & Silva, J. P. A. d. (2022). Comparative Analysis of Pedigree-Based BLUP and Phenotypic Mass Selection for Developing Elite Inbred Lines, Based on Field and Simulated Data. Agronomy, 12(10), 2560. https://doi.org/10.3390/agronomy12102560

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop