GWAS and genetic and phenotypic correlations of plasma metabolites with complete blood count traits in healthy young pigs reveal implications for pig immune response

Introduction: In this study estimated genetic and phenotypic correlations between fifteen complete blood count (CBC) traits and thirty-three heritable plasma metabolites in young healthy nursery pigs. In addition, it provided an opportunity to identify candidate genes associated with variation in metabolite concentration and their potential association with immune response, disease resilience, and production traits.

Methods: The blood samples were collected from healthy young pigs and Nuclear Magnetic Resonance (NMR) was used to quantify plasma metabolites. CBC was determined using the ADVIA® 2120i Hematology System. Genetic correlations of metabolite with CBC traits and single step genome-wide association study (ssGWAS) were estimated using the BLUPF90 programs.

Results: Results showed low phenotypic correlation estimates between plasma metabolites and CBC traits. The highest phenotypic correlation was observed between lactic acid and plasma basophil concentration (0.36 ± 0.04; p < 0.05). Several significant genetic correlations were found between metabolites and CBC traits. The plasma concentration of proline was genetically positively correlated with hemoglobin concentration (0.94 ± 0.03; p < 0.05) and L-tyrosine was negatively correlated with mean corpuscular hemoglobin (MCH; −0.92 ± 0.74; p < 0.05). The genomic regions identified in this study only explained a small percentage of the genetic variance of metabolites levels that were genetically correlated with CBC, resilience, and production traits.

Discussion: The results of this systems approach suggest that several plasma metabolite phenotypes are phenotypically and genetically correlated with CBC traits, suggesting that they may be potential genetic indicators of immune response following disease challenge. Genomic analysis revealed genes and pathways that might interact to modulate CBC, resilience, and production traits.

Accuracy of genomic prediction of shell quality in a White Leghorn line

Several genomic methods were applied for predicting shell quality traits recorded at 4 different hen ages in a White Leghorn line. The accuracies of genomic prediction of single-step GBLUP and single-trait Bayes B were compared with predictions of breeding values based on pedigree-BLUP under single-trait or multitrait models. Breaking strength (BS) and dynamic stiffness (Kdyn) measurements were collected on 18,524 birds from 3 consecutive generations, of which 4,164 animals also had genotypes from an Affymetrix 50K panel containing 49,591 SNPs after quality control edits. All traits had low to moderate heritability, ranging from 0.17 for BS to 0.34 for Kdyn. The highest accuracies of prediction were obtained for the multitrait single-step model. The use of marker information resulted in higher prediction accuracies than pedigree-based models for almost all traits. A genome-wide association study based on a Bayes B model was conducted to detect regions explaining the largest proportion of genetic variance. Across all 8 shell quality traits analyzed, 7 regions each explaining over 2% of genetic variance and 54 regions each explaining over 1% of genetic variance were identified. The windows explaining a large proportion of genetic variance overlapped with several potential candidate genes with biological functions linked to shell formation. A multitrait repeatability model using a single-step method is recommended for genomic evaluation of shell quality in layer chickens.

Investigating the genetic determination of clutch traits in laying hens

Clutch traits were proposed as a more detailed description of egg-laying patterns than simple total egg production. In this study, egg production of 23,809 Rhode Island Red (RIR) and 22,210 White Leghorn (WL) hens was described in terms of number of clutches, average and maximum clutch size, age at first egg, total saleable egg production, and percentage of egg defects. Genetic parameters were estimated using a six-trait animal model. Of the phenotyped birds, 1433 RIR hens and 1515 WL hens were genotyped with line specific 50K Affymetrix Axiom single nucleotide polymorphism chips to perform genome-wide association analyses. Moderate heritabilities were estimated for clutch traits of 0.20 to 0.42 in the RIR line and 0.29 to 0.41 in the WL line. Average and maximum clutch size was positively genetically correlated with total saleable egg number in both lines. Genome-wide association analysis identified seven regions that were associated with egg production in the RIR line and 12 regions in the WL line. The regions identified were line and trait specific, except for one region on chromosome 6 from 28 to 29 Mb that influenced number of clutches and maximum and average clutch size in WL hens. Regions associated with egg production identified here overlapped with 260 genes, with some strong positional candidates based on gene ontology including WASH1, which is involved in oocyte maturation, NPVF, involved in regulation of follicle-stimulating hormone secretion, and FOXO3, involved in oocyte maturation and ovulation from the ovarian follicle. Confirmation of the role of these genes in regulation of egg production pattern will require further studies.

Effect of lower-energy, higher-fiber diets on pigs divergently selected for residual feed intake when fed higher-energy, lower-fiber diets

Residual feed intake (RFI) is the difference between observed and predicted feed intake of an animal, based on growth and maintenance requirements. In Yorkshire pigs, divergent selection for increased (Low RFI) and decreased (High RFI) RFI was carried out over 10 generations (G) while feeding a corn- and soybean-meal-based, higher-energy, lower-fiber (HELF) diet. In G8 to G10, representing 4 replicates, barrows and gilts (n = 649) of the RFI lines were fed the HELF diet and a diet incorporating coproducts that were lower in energy and higher in dietary fiber (LEHF). The diets differed in ME, 3.32 vs. 2.87 Mcal/kg, and in neutral detergent fiber (NDF), 9.4% vs. 25.9%, respectively. The impact of the LEHF diet on 1) performance and growth, 2) diet digestibility, 3) genetic parameter estimates, and 4) responses to selection for RFI, when fed the HELF, was assessed. In general, the LEHF diet reduced the performance of both lines. When fed the HELF diet, the Low RFI pigs had lower (P < 0.05) ADFI (-12%), energy intake (-12%), ADG (-6%), and backfat depth (-12%); similar (P > 0.05) loin muscle area (LMA; +5%); and greater (P < 0.05) feed efficiency (i.e., 8% higher G:F and 7% lower RFI) than the High RFI line. These patterns of line differences were still present under the LEHF diet but differences for ADFI (-11%), energy intake (-10%), G:F (+2%), and RFI (-6%) were reduced compared to the HELF diet. Apparent total tract digestibility (ATTD) of the HELF and LEHF diets was assessed using 116 barrows and gilts from G8. When fed the HELF diet, ATTD of DM, GE, N, and NDF were similar between lines (P ≥ 0.27), but when fed the LEHF diet, the Low RFI pigs had greater digestibility (7%, 7%, 10%, and 32%) than the High RFI line (P ≤ 0.04). To measure responses to selection for RFI and estimate genetic parameters, data from all 10 generations were used (HELF; n = 2,310; LEHF, n = 317). Heritability estimates of performance traits ranged from 0.19 to 0.63, and genetic correlations of traits between diets were high and positive, ranging from 0.87 (RFI) to 0.99 (LMA). By G10, RFI in the Low RFI line was 3.86 and 1.50 genetic SD lower than in the High RFI line when fed the HELF and LEHF diets, respectively. Taken together, the results of this study demonstrate that responses to selection for RFI when fed a HELF diet are not fully realized when pigs are fed an extremely LEHF diet. Thus, feeding diets that differ from those used for selection may not maximize genetic potential for feed efficiency.

Effect of a major quantitative trait locus for porcine reproductive and respiratory syndrome (PRRS) resistance on response to coinfection with PRRS virus and porcine circovirus type 2b (PCV2b) in commercial pigs, with or without prior vaccination for PRRS

A major QTL for host response to porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) infection was identified in a previous study. Single nucleotide polymorphism WUR10000125 (WUR), which is in complete linkage disequilibrium with the putative causative mutation, can be used as a tag SNP for the QTL. However, the effect of WUR following PRRS vaccination and/or coinfection with other pathogens is not known. Therefore, objectives of this study were to estimate the effect of WUR on host response following PRRS vaccination and coinfection of PRRSV with porcine circovirus type 2b (PCV2b), to estimate genetic parameters for host response to vaccination and coinfection, and to estimate the effect of previously identified candidate SNP under PRRSV-only or PCV2b-only infection on host response to coinfection. Data from 2 trials, comprising a total of 396 commercial crossbred nursery pigs from a single genetic source, were used for all analyses. Pigs were preselected based on WUR genotype: approximately half AA and half AB, where B is the favorable and dominant allele. At weaning, pigs were shipped to Kansas State University, where half of the pigs were vaccinated with a PRRS modified live virus vaccine. Four weeks later, all pigs were coinfected with field strains of PRRSV and PCV2b and followed for 42 d. Body weight and serum viremia measurements were collected following vaccination and coinfection to calculate ADG and viral load (VL), respectively. Average heritability estimates for PRRS VL, PCV2b VL, and ADG were 0.29, 0.09, and 0.40, respectively. After vaccination, AB pigs had lower vaccination VL ( = 0.03) and faster gain ( = 0.004) than AA pigs, as expected. After coinfection, AB pigs had lower PRRSV VL ( < 0.001) but did not significantly differ from AA pigs in growth rate ( = 0.86). For PCV2b VL, suggestive evidence of an interaction between vaccination and WUR genotype ( = 0.11) was detected, where AB pigs had significantly lower PCV2b VL when vaccinated ( = 0.007) but not when they were not vaccinated ( = 0.87). In addition to WUR, several PRRS-associated SNP and a PCV2b-associated SNP had significant effects on host response to coinfection. In conclusion, marker-assisted selection based on WUR genotype alone, or along with other candidate SNP for PRRSV and PCV2b infection, is a promising strategy to select for improved host response to not just PRRS but also coinfection of PRRSV with PCV2b and perhaps other pathogens.

A major gene for host response to porcine reproductive and respiratory syndrome is not unfavorably associated with overall performance under nonchallenging conditions in commercial pig lines

A QTL for host response to porcine reproductive and respiratory syndrome (PRRS) was identified in a previous study. The SNP WUR10000125 (WUR) is used as a tag SNP for this QTL. The favorable (B) allele at this SNP is in low frequency in commercial populations, possibly because this allele is unfavorably associated with an important trait under nonchallenging conditions and, therefore, may have been selected against. Therefore, objectives of this study were to estimate the effect of WUR on traits under selection in commercial lines under nonchallenging conditions and to estimate the effect of WUR genotype of parents on performance of crossbred progeny in a commercial-like environment. Data were collected on 4 purebred lines: a Landrace dam line (D1), a Large White dam line (D2), a synthetic sire line (S1), and a Pietrain sire line (S2). Traits analyzed included total number born, number stillborn, farrowing survival, lactation survival, litter mortality, daily feed intake (DFI), backfat, average test daily gain (TGR), average lifetime daily gain (LGR), and Topigs Norsvin selection index (TSI) value, indicative of overall economic value. Deregressed EBV were calculated for each trait (except TSI) and analyzed within line. In the S1 line, AB and BB pigs had significantly lower TGR (P = 0.002) and LGR (P = 0.001) than AA pigs but also lower DFI (P = 0.004). Conversely, AB and BB pigs had significantly higher DFI (P < 0.001) and AB pigs had significantly higher TGR (P = 0.03) than AA pigs in the S2 line. The effect of WUR on TSI was not significant for any line (P ≥ 0.15). Analyses of phenotypic records collected on crossbred progeny of S1 sires and D1 × D2 F₁ females showed no significant effect of parent WUR genotype on DFI, backfat, TGR, or LGR (P ≥ 0.07). In conclusion, the effect of WUR was nonsignificant for most traits but the magnitude and direction of the effect differed by trait and by line. The favorable allele for host response to PRRS was associated with greater DFI and a tendency for greater TGR in the S2 line, but the opposite direction of effect was detected for the S1 line. Regardless of the effect on individual traits, no significant effect of WUR on TSI was detected for any line. Therefore, selecting for the B allele is expected to result in progeny with increased resistance to PRRS without compromising overall economic value under normal, nonchallenging conditions.

Genomewide association analysis of sow lactation performance traits in lines of Yorkshire pigs divergently selected for residual feed intake during grow-finish phase

Lactation is an economically and biologically important phase in the life cycle of sows. Short generation intervals in nucleus herds and low heritability of traits associated with lactation along with challenges associated with collecting accurate lactation performance phenotypes emphasize the importance of using genomic tools to examine the underlying genetics of these traits. We report the first genomewide association study (GWAS) on traits associated with lactation and efficiency in 2 lines of Yorkshire pigs that were divergently selected for residual feed intake during grow-finish phase. A total of 862 farrowing records from 2 parities were analyzed using a Bayesian whole genome variable selection model (Bayes B) to locate 1-Mb regions that were most strongly associated with each trait. The GWAS was conducted separately for parity 1 and 2 records. Marker-based heritabilities ranged from 0.03 to 0.39 for parity 1 traits and from 0.06 to 0.40 for parity 2 traits. For all traits studied, around 90% of genetic variance came from a large number of genomic regions with small effects, whereas genomic regions with large effects were found to be different for the same trait measured in parity 1 and 2. The highest percentage of genetic variance explained by a 1-Mb window for each trait ranged from 0.4% for feed intake during lactation to 4.2% for back fat measured at farrowing in parity 1 sows and from 0.2% for lactation feed intake to 5.4% for protein mass loss during lactation in parity 2 sows. A total of thirteen 1-Mb nonoverlapping windows were found to explain more than 1.5% of genetic variance for either a single trait or across multiple traits. These 1-Mb windows were on chromosomes 2, 3, 6, 7, 8, 11, 14, 15, 17, and 18. The major positional candidate genes within 1 Mb upstream and downstream of these windows were , (SSC2), (SSC6) (SSC7), (SSC8), (SSC11), (SSC14), (SSC17). Further validation studies on larger populations are required to validate these findings and to improve our understanding of the biology and complex genetic architecture of traits associated with sow lactation performance.

Estimation of genetic parameters for traits associated with reproduction, lactation, and efficiency in sows

Increased milk production due to high litter size, coupled with low feed intake, results in excessive mobilization of sow body reserves during lactation, which can have detrimental effects on future reproductive performance. A possibility to prevent this is to improve sow lactation performance genetically, along with other traits of interest. The aim of this study was to estimate breed-specific genetic parameters (by parity, between parities, and across parities) for traits associated with lactation and reproduction in Yorkshire and Landrace sows. Performance data were available for 2,107 sows with 1 to 3 parities (3,424 farrowings total). Sow back fat, loin depth and BW at farrowing, sow feed intake (SFI), and body weight loss (BWL) during lactation showed moderate heritabilities (0.21 to 0.37) in both breeds, whereas back fat loss (BFL), loin depth loss (LDL), and litter weight gain (LWG) showed low heritabilities (0.12 to 0.18). Among the efficiency traits, sow lactation efficiency showed extremely low heritability (near zero) in Yorkshire sows but a slightly higher (0.05) estimate in Landrace sows, whereas sow residual feed intake (SRFI) and energy balance traits showed moderate heritabilities in both breeds. Genetic correlations indicated that SFI during lactation had strong negative genetic correlations with body resource mobilization traits (BWL, BFL, and LDL; -0.35 to -0.70), and tissue mobilization traits in turn had strong positive genetic correlations with LWG (+0.24 to +0.54; < 0.05). However, SFI did not have a significant genetic correlation with LWG. These genetic correlations suggest that SFI during lactation is predominantly used for reducing sow body tissue losses, rather than for milk production. Estimates of genetic correlations for the same trait measured in parities 1 and 2 ranged from 0.64 to 0.98, which suggests that first and later parities should be treated as genetically different for some traits. Genetic correlations estimated between traits in parities 1 and 2 indicated that BWF and BWL measured in parity 1 can be used as indicator traits for SFI and SRFI measured in parities 1 and 2. In conclusion, traits associated with lactation in sows have a sizable genetic component and show potential for genetic improvement.

Repeatability vs. multiple-trait models to evaluate shell dynamic stiffness for layer chickens

Shell quality is one of the most important traits for improvement in layer chickens. Proper consideration of repeated records can increase the accuracy of estimated breeding values and thus genetic improvement of shell quality. The objective of this study was to compare different models for genetic evaluation of the collected data. For this study, 81,646 dynamic stiffness records on 21,321 brown egg layers and 93,748 records on 24,678 white egg layers from 4 generations were analyzed. Across generations, data were collected at 2 to 4 ages (at approximately 26, 42, 65, and 86 wk), with repeated records at each age. Seven models were compared, including 5 repeatability models with increasing complexity, a random regression model, and a multitrait model. The models were compared using Akaike Information Criteria with significance testing of nested models with a Log Likelihood Ratio test. Estimates of heritability were 0.31-0.36 for the brown line and 0.23-0.26 for the white line, but repeatability was higher for the model with age-specific permanent environment effects (0.59 for both lines) than for the model with an overall permanent environmental effects (0.47 for the brown and 0.41 for the white line). The model that allowed for permanent environmental effect within age and heterogeneous residual variance between ages resulted in improved fit compared to the traditional model that fits single permanent environment and residual effects, but was inferior in fit and predictive ability to the full multiple-trait model. The random regression model had better fit to the data than repeatability models but slightly worse than the multiple-trait model. For traits with repeated records at different ages, repeatability within and across ages as well as genetic correlations should be considered while choosing the number of records collected per individual as well as the model for genetic evaluation.

Genomewide association of piglet responses to infection with one of two porcine reproductive and respiratory syndrome virus isolates

Porcine reproductive and respiratory syndrome (PRRS) is a devastating disease in the swine industry. Identification of host genetic factors that enable selection for improved performance during PRRS virus (PRRSV) infection would reduce the impact of this disease on animal welfare and production efficiency. We conducted genomewide association study (GWAS) analyses of data from 13 trials of approximately 200 commercial crossbred nursery-age piglets that were experimentally infected with 1 of 2 type 2 isolates of PRRSV (NVSL 97-7985 [NVSL] and KS2006-72109 [KS06]). Phenotypes analyzed were viral load (VL) in blood during the first 21 d after infection (dpi) and weight gain (WG) from 0 to 42 dpi. We accounted for the previously identified QTL in the region on SSC4 in our models to increase power to identify additional regions. Many regions identified by single-SNP analyses were not identified using Bayes-B, but both analyses identified the same regions on SSC3 and SSC5 to be associated with VL in the KS06 trials and on SSC6 in the NVSL trials ( < 5 × 10); for WG, regions on SSC5 and SSC17 were associated in the NVSL trials ( < 3 × 10). No regions were identified with either method for WG in the KS06 trials. Except for the region on SSC4, which was associated with VL for both isolates (but only with WG for NVSL), identified regions did not overlap between the 2 PRRSV isolate data sets, despite high estimates of the genetic correlation between isolates for traits based on these data. We also identified genomic regions whose associations with VL or WG interacted with either PRRSV isolate or with genotype at the SSC4 QTL. Gene ontology (GO) annotation terms for genes located near moderately associated SNP ( < 0.003) were enriched for multiple immunologically (VL) and metabolism- (WG) related GO terms. The biological relevance of these regions suggests that, although it may increase the number of false positives, the use of single-SNP analyses and a relaxed threshold also increased the identification of true positives. In conclusion, although only the SSC4 QTL was associated with response to both PRRSV isolates, genes near associated SNP were enriched for the same GO terms across PRRSV isolates, suggesting that host responses to these 2 isolates are affected by the actions of many genes that function together in similar biological processes.

Effect of selection for residual feed intake during the grow/finish phase of production on sow reproductive performance and lactation efficiency

As feed costs continue to rise and efficiency during finishing is emphasized, the impact of selecting for more efficient grow/finish pigs on reproductive performance and feed efficiency of sows must be evaluated. Therefore, the objectives of this study were to evaluate correlated responses for sow reproductive performance and lactation feed efficiency to selection for residual feed intake (RFI) during the grow/finish phase of production (RFI) in 2 selection lines of pigs developed at Iowa State University (Ames, IA) and to estimate heritabilities of these traits. One line was selected over 7 generations for decreased RFIG/F (low RFI [LRFI] line) and the other line was randomly selected for 5 generations and then selected for increased RFIG/F (high RFI [HRFI] line). After 7 generations of selection, LRFI sows had 1.0 more piglets farrowed ( = 0.11) compared with HRFI sows, 1.3 more pigs born alive ( < 0.05), similar farrowing survival, 0.4 fewer mummies ( < 0.01), and more piglets weaned, both by litter (1.6 more; < 0.01) and by sow (1.1 more; < 0.01). Low RFI sows consumed 25 kg less feed and lost 9.8 kg more BW, 7.0 kg more fat mass, and 3.1 mm more backfat than HRFI sows ( < 0.001) during lactation. Although LRFI sows had a greater negative energy balance (-19.8 vs. -8.0 MJ ME/d; < 0.001), they had better RFI during lactation (-28.6 vs. 8.2 kg; < 0.0001), and the trend was for LRFI sows to have better lactation efficiency (61.3 vs. 57.8%; = 0.47) than HRFI sows. Heritabilities for sow weights, sow body composition, sow maintenance requirements (estimated from BW), and piglet birth weight were high ( > 0.4, SE < 0.07). Traits pertaining to piglet growth during lactation and mobilization of body tissue of the sow were moderately heritable (0.2 < < 0.4, SE < 0.07). In conclusion, selection for decreased RFI has favorably affected piglet performance and lactation efficiency but has unfavorably affected sow body condition loss and energy balance during lactation. These results indicate that pigs selected for increased efficiency during grow-finish are better able to direct resources where needed during other life history phases, that is, reproduction and lactation.

Assessing peripheral blood cell profile of Yorkshire pigs divergently selected for residual feed intake

The cost of feed is a serious issue in the pork industry, contributing about 65 to 75% of the total production cost. To prevent economic losses and decreased productivity of the herd, it is important to select for animals that eat less for the same lean gain, or more efficient animals. Residual feed intake (RFI) is the difference between observed feed intake and expected feed intake based on estimated maintenance and production requirements. Selection for decreased RFI, or more efficient animals, is a potential solution to higher feed costs in pig production. However, animals that are highly selected for decreased RFI may have reduced energy input to the immune system and fail to withstand diseases and stressors after infection that negatively impact profitability. The objective of this study was to evaluate differences in circulating blood cell profiles at a young age between 2 lines of Yorkshire pigs that were divergently selected for RFI as well as the heritability of these traits, to investigate effects of selection for RFI on immune system parameters, and to identify potential biomarkers for feed efficiency. Previous work has shown that the 2 lines had diverged for IGF-1 in serum in young pigs and, therefore, this stage was investigated for other potential physiological differences. Blood samples were drawn for a complete blood count (CBC) analysis from 517 gilts and barrows, ages 35 to 42 d, across the 2 lines. In general, the low-RFI line had lower numbers of specific types of white blood cells but higher hemoglobin concentration and red blood cell volume compared to the high-RFI line. No significant correlations were found between CBC traits and RFI across and within the lines (0.05 < < 0.1). Of the 15 CBC traits that were measured, 3 were highly heritable (0.56 < < 0.62), 9 were moderately heritable (0.12 < < 0.47), and 3 were lowly heritable ( < 0.12), suggesting a substantial genetic component for CBC traits and that selection for CBC traits could be effective. Our results also show that selection for RFI has significantly impacted the number of circulating blood cells. In this experiment, we studied only healthy animals that were not under known pathogen challenge; therefore, our results cannot be directly applied to a disease challenge situation. Future work will be to challenge the animals and determine the effect of challenge on CBC levels.

Composition and quality characteristics of carcasses from pigs divergently selected for residual feed intake on high- or low-energy diets

The objective was to determine the extent to which feeding low-energy, high-fiber (LEHF) and high-energy, low-fiber (HELF) diets impacts meat quality and carcass composition of pigs divergently selected for residual feed intake (RFI). Two experiments were conducted in the divergently selected Iowa State University RFI lines: Exp. 1 evaluated carcasses of generation (G) 8 pigs fed on commercial feeders; Exp. 2 evaluated composition, pork quality, sensory, and postmortem proteolysis of pigs fed on electronic single-space feeders in G 8 and 9. Pigs (N = 177) in Exp. 1 were randomly assigned a pen (mixed sex and line; N = 8). Groups (n = 3) of pigs were slaughtered at a mean BW of 121.5 kg. Pigs in Exp. 2 (G8: n = 158; G9: n = 157) were randomly assigned to 1 of 6 pens of each diet per G. Pigs from G8 were slaughtered at a mean BW of 122.5 kg and G9 at a mean of 128.4 kg. Data were analyzed using the mixed procedure of SAS. Fixed effects were line, diet, sex, and all appropriate interactions. Random effects were group, pen, litter, and sire and covariate of off-test BW. For Exp. 2, G was added as a fixed effect and sensory day was added as a random effect when applicable. In Exp. 1, carcasses from low RFI (LRFI) pigs were leaner and had less fat depth (P < 0.01). Carcasses from pigs fed the LEHF diet had a lighter HCW and greater estimated percent lean than pigs fed HELF diet (P < 0.01). In Exp. 2, LRFI pigs on the HELF diet had the greatest loin depth (P < 0.01). Chops from HRFI pigs had greater drip loss, color scores, lean tissue a*, and percent lipid and lesser percent moisture than LRFI ( P< 0.05). Chops from pigs on the LEHF diet had lesser muscle L* values and greater percent moisture than chops from pigs fed the HELF diet (P < 0.05). Chops from LRFI pigs were juicer than those from HRFI pigs (P < 0.05). Protein extracted at d 2 postmortem from LRFI pigs on the LEHF diet had a greater 38 kDa desmin degradation product than protein from LRFI pigs fed the HELF diet (P < 0.05). Day 5 postmortem extracted protein from HRFI pigs had greater 38 kDa desmin degradation product than LRFI (P = 0.05). Pigs fed LEHF (P < 0.01) had adipose with a greater iodine value than adipose from HELF fed pigs. Pork sensory quality from pigs differentially selected for residual feed intake was not influenced by energy content of the diet the pigs were fed.

The effect of selection for residual feed intake on scale-activity and scale-exit scores in Yorkshire gilts

The objectives of the present study were to use scale-activity and scale-exit scores to (1) determine the extent to which selection for improved residual feed intake (RFI) correlates with individual gilt temperament, (2) determine whether the timing of assessment affects the scores, and (3) determine whether gilts habituate to the process of weighing. To achieve these objectives, 192 grow–finish purebred Yorkshire gilts (39.7 ± 5.7 kg; 98 ± 6 days old) were utilised. Ninety-six gilts were from the 5th generation of a line selected for low residual feed intake (LRFI), and 96 gilts were from the 5th generation of a randomly selected control line (CRFI). Gilt activity and exit behaviour were scored on the weigh scale. The activity score was taken at two time points, namely, t = 0 (immediately on the back gate closing on the weight scale) and t = 15 (15 s later). Activity was assessed using a 5-point scoring system (1 = calm, minimal movement; 5 = continuous rapid movement and an escape attempt), and exit score was assessed using a 3-point scoring system (1 = no encouragement needed, full exit; 3 = encouragement needed). Gilts were weighed once every 2 weeks (each weighing considered a testing round) for a maximum of eight scores per gilt. Statistical differences were identified in all rounds for activity between genetic lines, except Rounds 2, 4 and 5. For both lines, gilt-activity scores decreased over rounds. The LRFI line began with a lower activity score (2.31 vs 2.65 ± 0.13), but did not experience as great of a drop in average score over rounds as did the CRFI gilts (0.62 vs 1.21). By the end of the trial, the CRFI gilts scored lower than the LRFI gilts. Activity was greater at t = 15 than at t = 0. Although timing affected the score, the pattern was similar, and so consistency will be more important than timing procedure. For exit scores, in Rounds 5 and 7, the LRFI line scored lower than the CRFI line, and there was a trend for the LRFI line to score lower in Round 6. Across all gilts, the mean exit score increased slowly throughout the trial. In conclusion, selection for lower RFI in purebred Yorkshires has an effect on activity score but this relationship is complex and warrants further research. The observation that gilts become habituated to the process of weighing supports the conclusion that the process is not aversive.

Investigation of the efficacy of albumin removal procedures on porcine serum proteome profile

Improving the ability to predict livestock performance using biomarkers will provide a benefit for livestock genetic evaluation and improvement. The most practical biological sample to screen for development of biomarkers is serum due to the ease of collection. However, protein profiles in serum are complex and dynamic. Strategies are needed to manage variation in serum proteins used for biomarker identification. Albumin is the most abundant protein in serum, comprising over 50% of the overall protein content, and has historically been depleted from serum before biomarker identification. The objective of this study was to investigate the use of gel-based proteomic techniques to evaluate the need for porcine albumin depletion in biomarker identification. Albumin is known to bind many proteins in the blood, thus potential biomarkers could be removed during albumin depletion. Using two-dimensional difference in gel electrophoresis (2D-DIGE), we show whole serum can be used for biomarker discovery. The data obtained show that albumin removal methods are effective for porcine sera. Over 85% of the protein spots resolved on at least half of the gels were changed in abundance between whole and albumin depleted sera. Of the 204 protein spots significantly altered in abundance, 59 were changed over 400%. However, albumin removal also altered the serum proteome in an unpredictable manner; in the depleted sera, 86 protein spots were increased in abundance and 118 were decreased. Furthermore, the abundance of 59.4% of the protein spots in the albumin depleted samples had a larger standard error than whole sera. However, the resolution of albumin in 2D-DIGE analysis of whole sera permitted the detection and quantification of substantial numbers of proteins. Thus, it is proposed that whole serum can be used in a gel-based proteomics system for the identification of porcine biomarkers.

Analysis of blood leukocytes in a naturally occurring immunodeficiency of pigs shows the defect is localized to B and T cells

Severe combined immunodeficiency (SCID) is the result of a set of inherited genetic defects which render components of the immune response nonfunctional. In Arabian horses, Jack Russell terriers, and mice, the disorder is a consequence of the absence of T and B lymphocytes, while natural killer (NK) cell and other leukocyte populations remain intact. Preliminary analysis of a naturally acquired form of inherited SCID in a line of pigs showed several defects in the architecture and composition of secondary lymphoid organs. In this study, a quantitative assessment of lymphocyte populations in affected and normal littermates showed depleted T or B lymphocyte populations in affected pigs; however, NK cells and neutrophils were present in numbers comparable to unaffected littermates. The results indicate that the immune defect in pigs shares the same features as other SCID-affected species.

Genetic analysis of reproductive traits and antibody response in a PRRS outbreak herd

Porcine reproductive and respiratory syndrome (PRRS) is the most economically significant disease impacting pig production in North America, Europe, and Asia, causing reproductive losses such as increased rates of stillbirth and mummified piglets. The objective of this study was to explore the genetic basis of host response to the PRRS virus (PRRSV) in a commercial multiplier sow herd before and after a PRRS outbreak, using antibody response and reproductive traits. Reproductive data comprising number born alive (NBA), number alive at 24 h (NA24), number stillborn (NSB), number born mummified (NBM), proportion born dead (PBD), number born dead (NBD), number weaned (NW), and number of mortalities through weaning (MW) of 5,227 litters from 1,967 purebred Landrace sows were used along with a pedigree comprising 2,995 pigs. The PRRS outbreak date was estimated from rolling averages of farrowing traits and was used to split the data into a pre-PRRS phase and a PRRS phase. All 641 sows in the herd during the outbreak were blood sampled 46 d after the estimated outbreak date and were tested for anti-PRRSV IgG using ELISA (sample-to-positive [S/P] ratio). Genetic parameters of traits were estimated separately for the pre-PRRS and PRRS phase data sets. Sows were genotyped using the PorcineSNP60 BeadChip, and genome-wide association studies (GWAS) were performed using method Bayes B. Heritability estimates for reproductive traits ranged from 0.01 (NBM) to 0.12 (NSB) and from 0.01 (MW) to 0.12 (NBD) for the pre-PRRS and PRRS phases, respectively. S/P ratio had heritability (0.45) and strong genetic correlations with most traits, ranging from -0.72 (NBM) to 0.73 (NBA). In the pre-PRRS phase, regions associated with NSB and PBD explained 1.6% and 3% of the genetic variance, respectively. In the PRRS phase, regions associated with NBD, NSB, and S/P ratio explained 0.8%, 11%, and 50.6% of the genetic variance, respectively. For S/P ratio, 2 regions on SSC 7 (SSC7) separated by 100 Mb explained 40% of the genetic variation, including a region encompassing the major histocompatibility complex, which explained 25% of the genetic variance. These results indicate a significant genomic component associated with PRRSV antibody response and NSB in this data set. Also, the high heritability and genetic correlation estimates for S/P ratio during the PRRS phase suggest that S/P ratio could be used as an indicator of the impact of PRRS on reproductive traits.

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.

Accuracy of genomic prediction using an evenly spaced, low-density single nucleotide polymorphism panel in broiler chickens

One approach for cost-effective implementation of genomic selection is to genotype training individuals with a high-density (HD) panel and selection candidates with an evenly spaced, low-density (ELD) panel. The purpose of this study was to evaluate the extent to which the ELD approach reduces the accuracy of genomic estimated breeding values (GEBV) in a broiler line, in which 1,091 breeders from 3 generations were used for training and 160 progeny of the third generation for validation. All birds were genotyped with an Illumina Infinium platform HD panel that included 20,541 segregating markers. Two subsets of HD markers, with 377 (ELD-1) or 766 (ELD-2) markers, were selected as ELD panels. The ELD-1 panel was genotyped using KBiosciences KASPar SNP genotyping chemistry, whereas the ELD-2 panel was simulated by adding markers from the HD panel to the ELD-1 panel. The training data set was used for 2 traits: BW at 35 d on both sexes and hen house production (HHP) between wk 28 and 54. Methods Bayes-A, -B, -C and genomic best linear unbiased prediction were used to estimate HD-marker effects. Two scenarios were used: (1) the 160 progeny were ELD-genotyped, and (2) the 160 progeny and their dams (117 birds) were ELD-genotyped. The missing HD genotypes in ELD-genotyped birds were imputed by a Gibbs sampler, capitalizing on linkage within families. In scenario (1), the correlation of GEBV for BW (HHP) of the 160 progeny based on observed HD versus imputed genotypes was greater than 0.94 (0.98) with the ELD-1 panel and greater than 0.97 (0.99) with the ELD-2 panel. In scenario (2), the correlation of GEBV for BW (HHP) was greater than 0.92 (0.96) with the ELD-1 panel and greater than 0.95 (0.98) with the ELD-2 panel. Hence, in a pedigreed population, genomic selection can be implemented by genotyping selection candidates with about 400 ELD markers with less than 6% loss in accuracy. This leads to substantial savings in genotyping costs, with little sacrifice in accuracy.

Validation and further characterization of a major quantitative trait locus associated with host response to experimental infection with porcine reproductive and respiratory syndrome virus

Infectious diseases are costly to the swine industry; porcine reproductive and respiratory syndrome (PRRS) is the most devastating. In earlier work, a quantitative trait locus associated with resistance/susceptibility to PRRS virus was identified on Sus scrofa chromosome 4 using approximately 560 experimentally infected animals from a commercial cross. The favorable genotype was associated with decreased virus load and increased weight gain (WG). The objective here was to validate and further characterize the association of the chromosome 4 region with PRRS resistance using data from two unrelated commercial crossbred populations. The validation populations consisted of two trials each of approximately 200 pigs sourced from different breeding companies that were infected with PRRS virus and followed for 42 days post-infection. Across all five trials, heritability estimates were 0.39 and 0.34 for viral load (VL; area under the curve of log-transformed viremia from 0 to 21 days post-infection) and WG to 42 days post-infection respectively. Effect estimates of SNP WUR10000125 in the chromosome 4 region were in the same directions and of similar magnitudes in the two new trials as had been observed in the first three trials. Across all five trials, the 1-Mb region on chromosome 4 explained 15 percent of genetic variance for VL and 11 percent for WG. The effect of the favorable minor allele at SNP WUR10000125 was dominant. Ordered genotypes for SNP WUR10000125 showed that the effect was present irrespective of whether the favorable allele was paternally or maternally inherited. These results demonstrate that selection for host response to PRRS virus infection could reduce the economic impact of PRRS.

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.

Evidence of decreased muscle protein turnover in gilts selected for low residual feed intake

The objective of this study was to evaluate the contribution of muscle protein turnover (synthesis and degradation) to the biological basis for genetic differences in finisher pigs selected for residual feed intake (RFI). Residual feed intake is defined as the difference between expected feed intake (based on the achieved rate of BW gain and backfat depth of individual pigs) and the observed feed intake of the individual pig. We hypothesized that protein turnover would be reduced in pigs selected for low RFI. Twelve gilts from a line selected for 7 generations for low RFI and 12 from a contemporary line selected for 2 generations for high RFI were paired by age and BW and fed a standard corn-soybean diet for 6 wk. Pigs were euthanized, muscle and liver samples were collected, and insulin signaling, protein synthesis, and protein degradation proteins were analyzed for expression and activities. Muscle from low RFI pigs tended to have less μ- and m-calpain activities (P = 0.10 and 0.09, respectively) and had significantly greater calpastatin activity and a decreased μ-calpain:calpastatin activity ratio (P < 0.05). Muscle from low RFI pigs had less 20S proteasome activity compared with their high RFI counterparts (P < 0.05). No differences in insulin signaling intermediates and translation initiation signaling proteins [mammalian target of rapamycin (mTOR) pathway] were observed (P > 0.05). Postmortem proteolysis was determined in the LM from the eighth generation of the low RFI pigs versus their high RFI counterparts (n = 9 per line). Autolysis of μ-calpain was decreased in the low RFI pigs and less troponin-T degradation product was observed at 3 d postmortem (P < 0.05), indicating slowed postmortem proteolysis during aging in the low RFI pigs. These data provide significant evidence that less protein degradation occurs in pigs selected for reduced RFI, and this may account for a significant portion of the increased efficiency observed in these animals.