Jesus Arango

Breeding value prediction for production traits in layer chickens using pedigree or genomic relationships in a reduced animal model

BackgroundGenomic selection involves breeding value estimation of selection candidates based on high-density SNP genotypes. To quantify the potential benefit of genomic selection, accuracies of estimated breeding values (EBV) obtained with different methods using pedigree or high-density SNP genotypes were evaluated and compared in a commercial layer chicken breeding line.MethodsThe following traits were analyzed: egg production, egg weight, egg color, shell strength, age at sexual maturity, body weight, albumen height, and yolk weight. Predictions appropriate for early or late selection were compared. A total of 2,708 birds were genotyped for 23,356 segregating SNP, including 1,563 females with records. Phenotypes on relatives without genotypes were incorporated in the analysis (in total 13,049 production records).The data were analyzed with a Reduced Animal Model using a relationship matrix based on pedigree data or on marker genotypes and with a Bayesian method using model averaging. Using a validation set that consisted of individuals from the generation following training, these methods were compared by correlating EBV with phenotypes corrected for fixed effects, selecting the top 30 individuals based on EBV and evaluating their mean phenotype, and by regressing phenotypes on EBV.ResultsUsing high-density SNP genotypes increased accuracies of EBV up to two-fold for selection at an early age and by up to 88% for selection at a later age. Accuracy increases at an early age can be mostly attributed to improved estimates of parental EBV for shell quality and egg production, while for other egg quality traits it is mostly due to improved estimates of Mendelian sampling effects. A relatively small number of markers was sufficient to explain most of the genetic variation for egg weight and body weight.

Genome-wide association analysis and genetic architecture of egg weight and egg uniformity in layer chickens.

The pioneering work by Professor Soller et al., among others, on the use of genetic markers to analyze quantitative traits has provided opportunities to discover their genetic architecture in livestock by identifying quantitative trait loci (QTL). The recent availability of high-density single nucleotide polymorphism (SNP) panels has advanced such studies by capitalizing on population-wide linkage disequilibrium at positions across the genome. In this study, genomic prediction model Bayes-B was used to identify genomic regions associated with the mean and standard deviation of egg weight at three ages in a commercial brown egg layer line. A total of 24,425 segregating SNPs were evaluated simultaneously using over 2900 genotyped individuals or families. The corresponding phenotypic records were represented as individual measurements or family means from full-sib progeny. A novel approach using the posterior distribution of window variances from the Monte Carlo Markov Chain samples was used to describe genetic architecture and to make statistical inferences about regions with the largest effects. A QTL region on chromosome 4 was found to explain a large proportion of the genetic variance for the mean (30%) and standard deviation (up to 16%) of the weight of eggs laid at specific ages. Additional regions with smaller effects on chromosomes 2, 5, 6, 8, 20, 23, 28 and Z showed suggestive associations with mean egg weight and a region on chromosome 13 with the standard deviation of egg weight at 26-28 weeks of age. The genetic architecture of the analyzed traits was characterized by a limited number of genes or genomic regions with large effects and many regions with small polygenic effects. The region on chromosome 4 can be used to improve both the mean and standard deviation of egg weight by marker-assisted selection.

Persistence of accuracy of genomic estimated breeding values over generations in layer chickens

BackgroundThe predictive ability of genomic estimated breeding values (GEBV) originates both from associations between high-density markers and QTL (Quantitative Trait Loci) and from pedigree information. Thus, GEBV are expected to provide more persistent accuracy over successive generations than breeding values estimated using pedigree-based methods. The objective of this study was to evaluate the accuracy of GEBV in a closed population of layer chickens and to quantify their persistence over five successive generations using marker or pedigree information.MethodsThe training data consisted of 16 traits and 777 genotyped animals from two generations of a brown-egg layer breeding line, 295 of which had individual phenotype records, while others had phenotypes on 2,738 non-genotyped relatives, or similar data accumulated over up to five generations. Validation data included phenotyped and genotyped birds from five subsequent generations (on average 306 birds/generation). Birds were genotyped for 23,356 segregating SNP. Animal models using genomic or pedigree relationship matrices and Bayesian model averaging methods were used for training analyses. Accuracy was evaluated as the correlation between EBV and phenotype in validation divided by the square root of trait heritability.ResultsPedigree relationships in outbred populations are reduced by 50% at each meiosis, therefore accuracy is expected to decrease by the square root of 0.5 every generation, as observed for pedigree-based EBV (Estimated Breeding Values). In contrast the GEBV accuracy was more persistent, although the drop in accuracy was substantial in the first generation. Traits that were considered to be influenced by fewer QTL and to have a higher heritability maintained a higher GEBV accuracy over generations. In conclusion, GEBV capture information beyond pedigree relationships, but retraining every generation is recommended for genomic selection in closed breeding populations.

Extent and Consistency of Linkage Disequilibrium and Identification of DNA Markers for Production and Egg Quality Traits in Commercial Layer Chicken Populations

BackgroundThe genome sequence and a high-density SNP map are now available for the chicken and can be used to identify genetic markers for use in marker-assisted selection (MAS). Effective MAS requires high linkage disequilibrium (LD) between markers and quantitative trait loci (QTL), and sustained marker-QTL LD over generations. This study used data from a 3,000 SNP panel to assess the level and consistency of LD between single nucleotide polymorphisms (SNPs) over consecutive years in two egg-layer chicken lines, and analyzed one line by two methods (SNP-wise association and genome-wise Bayesian analysis) to identify markers associated with egg-quality and egg-production phenotypes.ResultsThe LD between markers pairs was high at short distances (r2 > 0.2 at < 2 Mb) and remained high after one generation (correlations of 0.80 to 0.92 at < 5 Mb) in both lines. Single- and 3-SNP regression analyses using a mixed model with SNP as fixed effect resulted in 159 and 76 significant tests (P < 0.01), respectively, across 12 traits. A Bayesian analysis called BayesB, that fits all SNPs simultaneously as random effects and uses model averaging procedures, identified 33 SNPs that were included in the model >20% of the time (φ > 0.2) and an additional ten 3-SNP windows that had a sum of φ greater than 0.35. Generally, SNPs included in the Bayesian model also had a small P-value in the 1-SNP analyses.ConclusionHigh LD correlations between markers at short distances across two generations indicate that such markers will retain high LD with linked QTL and be effective for MAS. The different association analysis methods used provided consistent results. Multiple single SNPs and 3-SNP windows were significantly associated with egg-related traits, providing genomic positions of QTL that can be useful for both MAS and to identify causal mutations.

Genome-wide association study for egg production and quality in layer chickens

Discovery of genes with large effects on economically important traits has for many years been of interest to breeders. The development of SNP panels which cover the whole genome with high density and, more importantly, that can be genotyped on large numbers of individuals at relatively low cost, has opened new opportunities for genome-wide association studies (GWAS). The objective of this study was to find genomic regions associated with egg production and quality traits in layers using analysis methods developed for the purpose of whole genome prediction. Genotypes on over 4500 birds and phenotypes on over 13,000 hens from eight generations of a brown egg layer line were used. Birds were genotyped with a custom 42K Illumina SNP chip. Recorded traits included two egg production and 11 egg quality traits (puncture score, albumen height, yolk weight and shell colour) at early and late stages of production, as well as body weight and age at first egg. Egg weight was previously analysed by Wolc et al. (2012). The Bayesian whole genome prediction model--BayesB (Meuwissen et al. 2001) was used to locate 1 Mb regions that were most strongly associated with each trait. The posterior probability of a 1 Mb window contributing to genetic variation was used as the criterion for suggesting the presence of a quantitative trait locus (QTL) in that window. Depending upon the trait, from 1 to 7 significant (posterior probability >0.9) 1 Mb regions were found. The largest QTL, a region explaining 32% of genetic variance, was found on chr4 at 78 Mb for body weight but had pleiotropic effects on other traits. For the other traits, the largest effects were much smaller, explaining <7% of genetic variance, with regions on chromosomes 2, 12 and 17 explaining above 5% of genetic variance for albumen height, shell colour and egg production, respectively. In total, 45 of 1043 1 Mb windows were estimated to have a non-zero effect with posterior probability > 0.9 for one or more traits.

Genetic parameters of egg defects and egg quality in layer chickens

Genetic parameters were estimated for egg defects, egg production, and egg quality traits. Eggs from 11,738 purebred brown-egg laying hens were classified as salable or as having one of the following defects: bloody, broken, calcium deposit, dirty, double yolk, misshapen, pee-wee, shell-less, and soft shelled. Egg quality included albumen height, egg weight, yolk weight, and puncture score. Body weight, age at sexual maturity, and egg production were also recorded. Heritability estimates of liability to defects using a threshold animal model were less than 0.1 for bloody and dirty; between 0.1 and 0.2 for pee-wee, broken, misshapen, soft shelled, and shell-less; and above 0.2 for calcium deposit and double yolk. Quality and production traits were more heritable, with estimates ranging from 0.29 (puncture score) to 0.74 (egg weight). High-producing hens had a lower frequency of egg defects. High egg weight and BW were associated with an increased frequency of double yolks, and to a lesser extent, with more shell quality defects. Estimates of genetic correlations among defect traits that were related to shell quality were positive and moderate to strong (0.24-0.73), suggesting that these could be grouped into one category or selection could be based on the trait with the highest heritability or that is easiest to measure. Selection against defective eggs would be more efficient by including egg defect traits in the selection criterion, along with egg production rate of salable eggs and egg quality traits.

Variation in the ovocalyxin-32 gene in commercial egg-laying chickens and its relationship with egg production and egg quality traits.

Avian eggshell quality is an important trait for commercial egg production, as the eggshell is the primary packaging material and antimicrobial barrier for the internal food resource. Strong eggshells are essential to ensure that eggs can reach their final destination without damage. Ovocalyxin-32 (OCX32) is a matrix protein found within the outer layers of the eggshell and in the cuticle. Numerous reports in the literature have identified association between variants in the gene encoding this protein, OCX32, and various eggshell quality traits. Thus, OCX32 is a candidate gene for selection for eggshell traits in commercial poultry populations. Sequencing of exons 2-6 of the OCX32 gene in eight elite brown and white eggshell commercial egg-laying lines revealed 28 SNPs and one SNP/indel. Eighteen of these SNPs were predicted to alter the amino acid sequence of the protein. Clusters of SNPs in complete linkage disequilibrium were found in both exons 2 and 6. A total of 19 different versions or protein-sequence haplotypes of the OCX32 protein were inferred, revealing considerable variation within commercial lines. Genotypes for 13 of the SNPs were determined for 330-1819 individuals per line. Trait association studies revealed a significant effect of OCX32 on shell color in white egg lines and line-specific significant effects on albumen height, early egg weight, puncture score, and yolk weight. Three of the lines showed a significant change in OCX32 frequency over time, indicating selection pressure for certain variants of this gene during the breeding program.

Pedigree and genomic analyses of feed consumption and residual feed intake in laying hens

Efficiency of production is increasingly important with the current escalation of feed costs and demands to minimize the environmental footprint. The objectives of this study were 1) to estimate heritabilities for daily feed consumption and residual feed intake and their genetic correlations with production and egg-quality traits; 2) to evaluate accuracies of estimated breeding values from pedigree- and marker-based prediction models; and 3) to localize genomic regions associated with feed efficiency in a brown egg layer line. Individual feed intake data collected over 2-wk trial periods were available for approximately 6,000 birds from 8 generations. Genetic parameters were estimated with a multitrait animal model; methods BayesB and BayesCπ were used to estimate marker effects and find genomic regions associated with feed efficiency. Using pedigree information, feed efficiency was found to be moderately heritable (h(2) = 0.46 for daily feed consumption and 0.47 for residual feed intake). Hens that consumed more feed and had greater residual feed intake (lower efficiency) had a genetic tendency to lay slightly more eggs with greater yolk weights and albumen heights. Regions on chromosomes 1, 2, 4, 7, 13, and Z were found to be associated with feed intake and efficiency. The accuracy from genomic prediction was higher and more persistent (better maintained across generations) than that from pedigree-based prediction. These results indicate that genomic selection can be used to improve feed efficiency in layers.

Analysis of egg production in layer chickens using a random regression model with genomic relationships

Random regression models allow for analysis of longitudinal data, which together with the use of genomic information are expected to increase accuracy of selection, when compared with analyzing average or total production with pedigree information. The objective of this study was to estimate variance components for egg production over time in a commercial brown egg layer population using genomic relationship information. A random regression reduced animal model with a marker-based relationship matrix was used to estimate genomic breeding values of 3,908 genotyped animals from 6 generations. The first 5 generations were used for training, and predictions were validated in generation 6. Daily egg production up to 46 wk in lay was accumulated into 85,462 biweekly (every 2 wk) records for training, of which 17,570 were recorded on genotyped hens and the remaining on their nongenotyped progeny. The effect of adding additional egg production data of 2,167 nongenotyped sibs of selection candidates [16,037 biweekly (every 2 wk) records] to the training data was also investigated. The model included a 5th order Legendre polynomial nested within hatch-week as fixed effects and random terms for coefficients of quadratic polynomials for genetic and permanent environmental components. Residual variance was assumed heterogeneous among 2-wk periods. Models using pedigree and genomic relationships were compared. Estimates of residual variance were very similar under both models, but the model with genomic relationships resulted in a larger estimate of genetic variance. Heritability estimates increased with age up to mid production and decreased afterward, resulting in an average heritability of 0.20 and 0.33 for pedigree and genomic models. Prediction of total egg number was more accurate with the genomic than with the pedigree-based random regression model (correlation in validation 0.26 vs. 0.16). The genomic model outperformed the pedigree model in most of the 2-wk periods. Thus, results of this study show that random regression reduced animal models can be used in breeding programs using genomic information and can result in substantial improvements in the accuracy of selection for trajectory traits.

QTL mapping of egg albumen quality in egg layers.

BackgroundA fresh, good quality egg has a firm and gelatinous albumen that anchors the yolk and restricts growth of microbiological pathogens. As the egg ages, the gel-like structure collapses, resulting in thin and runny albumen. Occasionally thin albumen is found in a fresh egg, giving the impression of a low quality product. A mapping population consisting of 1599 F2 hens from a cross between White Rock and Rhode Island Red lines was set up, to identify loci controlling albumen quality. The phenotype for albumen quality was evaluated by albumen height and in Haugh units (HU) measured on three consecutive eggs from each F2 hen at the age of 40 weeks. For the fine-mapping analysis, albumen height and HU were used simultaneously to eliminate contribution of the egg size to the phenotype.ResultsLinkage analysis in a small population of seven half-sib families (668 F2) with 162 microsatellite markers spread across 27 chromosomes revealed two genome-wide significant regions with additive effects for HU on chromosomes 7 and Z. In addition, two putative genome-wide quantitative trait loci (QTL) regions were identified on chromosomes 4 and 26. The QTL effects ranged from 2 to 4% of the phenotypic variance. The genome-wide significant QTL regions on chromosomes 7 and Z were selected for fine-mapping in the full set composed of 16 half-sib families. In addition, their existence was confirmed by an association analysis in an independent commercial Hy-Line pure line.ConclusionsWe identified four chicken genomic regions that affect albumen quality. Our results also suggest that genes that affect albumen quality act both directly and indirectly through several different mechanisms. For instance, the QTL regions on both fine-mapped chromosomes 7 and Z overlapped with a previously reported QTL for eggshell quality, indicating that eggshell membranes may play a role in albumen quality.