Use of Host Feeding Behavior and Gut Microbiome Data in Estimating Variance Components and Predicting Growth and Body Composition Traits in Swine.

The purpose of this study was to investigate the use of feeding behavior in conjunction with gut microbiome sampled at two growth stages in predicting growth and body composition traits of finishing pigs. Six hundred and fifty-one purebred boars of three breeds: Duroc (DR), Landrace (LR), and Large White (LW), were studied. Feeding activities were recorded individually from 99 to 163 days of age. The 16S rRNA gene sequences were obtained from each pig at 123 ± 4 and 158 ± 4 days of age. When pigs reached market weight, body weight (BW), ultrasound backfat thickness (BF), ultrasound loin depth (LD), and ultrasound intramuscular fat (IMF) content were measured on live animals. Three models including feeding behavior (Model_FB), gut microbiota (Model_M), or both (Model_FB_M) as predictors, were investigated. Prediction accuracies were evaluated through cross-validation across genetic backgrounds using the leave-one-breed-out strategy and across rearing environments using the leave-one-room-out approach. The proportions of phenotypic variance of growth and body composition traits explained by feeding behavior ranged from 0.02 to 0.30, and from 0.20 to 0.52 when using gut microbiota composition. Overall prediction accuracy (averaged over traits and time points) of phenotypes was 0.24 and 0.33 for Model_FB, 0.27 and 0.19 for Model_M, and 0.40 and 0.35 for Model_FB_M for the across-breed and across-room scenarios, respectively. This study shows how feeding behavior and gut microbiota composition provide non-redundant information in predicting growth in swine.

Impact of including the cause of missing records on genetic evaluations for growth in commercial pigs.

It is of interest to evaluate crossbred pigs for hot carcass weight (HCW) and birth weight (BW); however, obtaining a HCW record is dependent on livability (LIV) and retained tag (RT). The purpose of this study is to analyze how HCW evaluations are affected when herd removal and missing identification are included in the model and examine if accounting for the reasons for missing traits improves the accuracy of predicting breeding values. Pedigree information was available for 1,965,077 purebred and crossbred animals. Records for 503,716 commercial three-way crossbred terminal animals from 2014 to 2019 were provided by Smithfield Premium Genetics. Two pedigree-based models were compared; model 1 (M1) was a threshold-linear model with all four traits (BW, HCW, RT, and LIV), and model 2 (M2) was a linear model including only BW and HCW. The fixed effects used in the model were contemporary group, sex, age at harvest (for HCW only), and dam parity. The random effects included direct additive genetic and random litter effects. Accuracy, dispersion, bias, and Pearson correlations were estimated using the linear regression method. The heritabilities were 0.11, 0.07, 0.02, and 0.04 for BW, HCW, RT, and LIV, respectively, with standard errors less than 0.01. No difference was observed in heritabilities or accuracies for BW and HCW between M1 and M2. Accuracies were 0.33, 0.37, 0.19, and 0.23 for BW, HCW, RT, and LIV, respectively. The genetic correlation between BW and RT was 0.34 ± 0.03, and between BW and LIV was 0.56 ± 0.03. Similarly, the genetic correlation between HCW and RT was 0.26 ± 0.04, and between HCW and LIV was 0.09 ± 0.05, respectively. The positive and moderate genetic correlations between BW and other traits imply a heavier BW resulted in a higher probability of surviving to harvest. Genetic correlations between HCW and other traits were lower due to the large quantity of missing records. Despite the heritable and correlated aspects of RT and LIV, results imply no major differences between M1 and M2; hence, it is unnecessary to include these traits in classical models for BW and HCW.

Genotype-by-environment interactions for reproduction, body composition, and growth traits in maternal-line pigs based on single-step genomic reaction norms.

BACKGROUND: There is an increasing need to account for genotype-by-environment (G × E) interactions in livestock breeding programs to improve productivity and animal welfare across environmental and management conditions. This is even more relevant for pigs because selection occurs in high-health nucleus farms, while commercial pigs are raised in more challenging environments. In this study, we used single-step homoscedastic and heteroscedastic genomic reaction norm models (RNM) to evaluate G × E interactions in Large White pigs, including 8686 genotyped animals, for reproduction (total number of piglets born, TNB; total number of piglets born alive, NBA; total number of piglets weaned, NW), growth (weaning weight, WW; off-test weight, OW), and body composition (ultrasound muscle depth, MD; ultrasound backfat thickness, BF) traits. Genetic parameter estimation and single-step genome-wide association studies (ssGWAS) were performed for each trait. RESULTS: The average performance of contemporary groups (CG) was estimated and used as environmental gradient in the reaction norm analyses. We found that the need to consider heterogeneous residual variance in RNM models was trait dependent. Based on estimates of variance components of the RNM slope and of genetic correlations across environmental gradients, G × E interactions clearly existed for TNB and NBA, existed for WW but were of smaller magnitude, and were not detected for NW, OW, MD, and BF. Based on estimates of the genetic variance explained by the markers in sliding genomic windows in ssGWAS, several genomic regions were associated with the RNM slope for TNB, NBA, and WW, indicating specific biological mechanisms underlying environmental sensitivity, and dozens of novel candidate genes were identified. Our results also provided strong evidence that the X chromosome contributed to the intercept and slope of RNM for litter size traits in pigs. CONCLUSIONS: We provide a comprehensive description of G × E interactions in Large White pigs for economically-relevant traits and identified important genomic regions and candidate genes associated with GxE interactions on several autosomes and the X chromosome. Implementation of these findings will contribute to more accurate genomic estimates of breeding values by considering G × E interactions, in order to genetically improve the environmental robustness of maternal-line pigs.

Determining the stability of accuracy of genomic estimated breeding values in future generations in commercial pig populations.

Genomic information has a limited dimensionality (number of independent chromosome segments [Me]) related to the effective population size. Under the additive model, the persistence of genomic accuracies over generations should be high when the nongenomic information (pedigree and phenotypes) is equivalent to Me animals with high accuracy. The objective of this study was to evaluate the decay in accuracy over time and to compare the magnitude of decay with varying quantities of data and with traits of low and moderate heritability. The dataset included 161,897 phenotypic records for a growth trait (GT) and 27,669 phenotypic records for a fitness trait (FT) related to prolificacy in a population with dimensionality around 5,000. The pedigree included 404,979 animals from 2008 to 2020, of which 55,118 were genotyped. Two single-trait models were used with all ancestral data and sliding subsets of 3-, 2-, and 1-generation intervals. Single-step genomic best linear unbiased prediction (ssGBLUP) was used to compute genomic estimated breeding values (GEBV). Estimated accuracies were calculated by the linear regression (LR) method. The validation population consisted of single generations succeeding the training population and continued forward for all generations available. The average accuracy for the first generation after training with all ancestral data was 0.69 and 0.46 for GT and FT, respectively. The average decay in accuracy from the first generation after training to generation 9 was -0.13 and -0.19 for GT and FT, respectively. The persistence of accuracy improves with more data. Old data have a limited impact on the predictions for young animals for a trait with a large amount of information but a bigger impact for a trait with less information.

Influence of host genetics in shaping the rumen bacterial community in beef cattle.

In light of recent host-microbial association studies, a consensus is evolving that species composition of the gastrointestinal microbiota is a polygenic trait governed by interactions between host genetic factors and the environment. Here, we investigated the effect of host genetic factors in shaping the bacterial species composition in the rumen by performing a genome-wide association study. Using a common set of 61,974 single-nucleotide polymorphisms found in cattle genomes (n = 586) and corresponding rumen bacterial community composition, we identified operational taxonomic units (OTUs), Families and Phyla with high heritability. The top associations (1-Mb windows) were located on 7 chromosomes. These regions were associated with the rumen microbiota in multiple ways; some (chromosome 19; position 3.0-4.0 Mb) are associated with closely related taxa (Prevotellaceae, Paraprevotellaceae, and RF16), some (chromosome 27; position 3.0-4.0 Mb) are associated with distantly related taxa (Prevotellaceae, Fibrobacteraceae, RF16, RFP12, S24-7, Lentisphaerae, and Tenericutes) and others (chromosome 23; position 0.0-1.0) associated with both related and unrelated taxa. The annotated genes associated with identified genomic regions suggest the associations observed are directed toward selective absorption of volatile fatty acids from the rumen to increase energy availability to the host. This study demonstrates that host genetics affects rumen bacterial community composition.

The impact of reducing the frequency of animals genotyped at higher density on imputation and prediction accuracies using ssGBLUP1.

The largest gains in accuracy in a genomic selection program come from genotyping young selection candidates who have not yet produced progeny and who might, or might not, have a phenotypic record recorded. To reduce genotyping costs and to allow for an increased amount of genomic data to be available in a population, young selection candidates may be genotyped with low-density (LD) panels and imputed to a higher density. However, to ensure that a reasonable imputation accuracy persists overtime, some parent animals originally genotyped at LD must be re-genotyped at a higher density. This study investigated the long-term impact of selectively re-genotyping parents with a medium-density (MD) SNP panel on the accuracy of imputation and on the genetic predictions using ssGBLUP in a simulated beef cattle population. Assuming a moderately heritable trait (0.25) and a population undergoing selection, the simulation generated sequence data for a founder population (100 male and 500 female individuals) and 9,000 neutral markers, considered as the MD panel. All selection candidates from generation 8 to 15 were genotyped with LD panels corresponding to a density of 0.5% (LD_0.5), 2% (LD_2), and 5% (LD_5) of the MD. Re-genotyping scenarios chose parents at random or based on EBV and ranged from 10% of male parents to re-genotyping all male and female parents with MD. Ranges in average imputation accuracy at generation 15 were 0.567 to 0.936, 0.795 to 0.985, and 0.931 to 0.995 for the LD_0.5, LD_2, and LD_5, respectively, and the average EBV accuracies ranged from 0.453 to 0.735, 0.631 to 0.784, and 0.748 to 0.807 for LD_0.5, LD_2, and LD_5, respectively. Re-genotyping parents based on their EBV resulted in higher imputation and EBV accuracies compared to selecting parents at random and these values increased with the size of LD panels. Differences between re-genotyping scenarios decreased when the density of the LD panel increased, suggesting fewer animals needed to be re-genotyped to achieve higher accuracies. In general, imputation and EBV accuracies were greater when more parents were re-genotyped, independent of the proportion of males and females. In practice, the relationship between the density of the LD panel used and the target panel must be considered to determine the number (proportion) of animals that would need to be re-genotyped to enable sufficient imputation accuracy.

The impact of clustering methods for cross-validation, choice of phenotypes, and genotyping strategies on the accuracy of genomic predictions.

For genomic predictors to be of use in genetic evaluation, their predicted accuracy must be a reliable indicator of their utility, and thus unbiased. The objective of this paper was to evaluate the accuracy of prediction of genomic breeding values (GBV) using different clustering strategies and response variables. Red Angus genotypes (n = 9,763) were imputed to a reference 50K panel. The influence of clustering method [k-means, k-medoids, principal component (PC) analysis on the numerator relationship matrix (A) and the identical-by-state genomic relationship matrix (G) as both data and covariance matrices, and random] and response variables [deregressed estimated breeding values (DEBV) and adjusted phenotypes] were evaluated for cross-validation. The GBV were estimated using a Bayes C model for all traits. Traits for DEBV included birth weight (BWT), marbling (MARB), rib-eye area (REA), and yearling weight (YWT). Adjusted phenotypes included BWT, YWT, and ultrasonically measured intramuscular fat percentage and REA. Prediction accuracies were estimated using the genetic correlation between GBV and associated response variable using a bivariate animal model. A simulation mimicking a cattle population, replicated 5 times, was conducted to quantify differences between true and estimated accuracies. The simulation used the same clustering methods and response variables, with the addition of 2 genotyping strategies (random and top 25% of individuals), and forward validation. The prediction accuracies were estimated similarly, and true accuracies were estimated as the correlation between the residuals of a bivariate model including true breeding value (TBV) and GBV. Using the adjusted Rand index, random clusters were clearly different from relationship-based clustering methods. In both real and simulated data, random clustering consistently led to the largest estimates of accuracy, while no method was consistently associated with more or less bias than other methods. In simulation, random genotyping led to higher estimated accuracies than selection of the top 25% of individuals. Interestingly, random genotyping seemed to overpredict true accuracy while selective genotyping tended to underpredict accuracy. When forward in time validation was used, DEBV led to less biased estimates of GBV accuracy. Results suggest the highest, least biased GBV accuracies are associated with random genotyping and DEBV.

Symposium review: The genomic architecture of inbreeding: How homozygosity affects health and performance.

Inbreeding depression is a growing concern in livestock because it can detrimentally affect animal fitness, health, and production levels. Genomic information can be used to more effectively capture variance in Mendelian sampling, thereby enabling more accurate estimation of inbreeding, but further progress is still required. The calculation of inbreeding for herd management purposes is largely still done using pedigree information only, although inbreeding coefficients calculated in this manner have been shown to be less accurate than genomic inbreeding measures. Continuous stretches of homozygous genotypes, so called runs of homozygosity, have been shown to provide a better estimate of autozygosity at the genomic level than conventional measures based on inbreeding coefficients calculated through conventional pedigree information or even genomic relationship matrices. For improved and targeted management of genomic inbreeding at the population level, the development of methods that incorporate genomic information in mate selection programs may provide a more precise tool for reducing the detrimental effects of inbreeding in dairy herds. Additionally, a better understanding of the genomic architecture of inbreeding and incorporating that knowledge into breeding programs could significantly refine current practices. Opportunities to maintain high levels of genetic progress in traits of interest while managing homozygosity and sustaining acceptable levels of heterozygosity in highly selected dairy populations exist and should be examined more closely for continued sustainability of both the dairy cattle population as well as the dairy industry. The inclusion of precise genomic measures of inbreeding, such as runs of homozygosity, inbreeding, and mating programs, may provide a path forward. In this symposium review article, we describe traditional measures of inbreeding and the recent developments made toward more precise measures of homozygosity using genomic information. The effects of homozygosity resulting from inbreeding on phenotypes, the identification and mapping of detrimental homozygosity haplotypes, management of inbreeding with genomic data, and areas in need of further research are discussed.

The impact of selective genotyping on the response to selection using single-step genomic best linear unbiased prediction.

Across the majority livestock species, routinely collected genomic and pedigree information has been incorporated into evaluations using single-step methods. As a result, strategies that reduce genotyping costs without reducing the response to selection are important as they could have substantial economic impacts on breeding programs. Therefore, the objective of the current study was to investigate the impact of selectively genotyping selection candidates on the selection response using simulation. Populations were simulated to mimic the genome and population structure of a swine and cattle population undergoing selection on an index comprised of the estimated breeding values (EBV) for 2 genetically correlated quantitative traits. Ten generations were generated and genotyping began generation 7. Two phenotyping scenarios were simulated that assumed the first trait was recorded early in life on all individuals and the second trait was recorded on all versus a random subset of the individuals. The EBV were generated from a bivariate animal model. Multiple genotyping scenarios were generated that ranged from not genotyping any selection candidates, a proportion of the selection candidates based on either their index value or chosen at random, and genotyping all selection candidates. An interim index value was utilized to decide who to genotype for the selective genotype strategy. The interim value assumed only the first trait was observed and the only genotypic information available was on animals in previous generations. Within each genotyping scenario 25 replicates were generated. Within each genotyping scenario the mean response per generation and the degree to which EBV were inflated/deflated was calculated. Across both species and phenotyping strategies, the plateau of diminishing returns was observed when 60% of the selection candidates with the largest index values were genotyped. When randomly genotyping selection candidates, either 80 or 100% of the selection candidates needed to be genotyped for there not to be a reduction in the index response. Across both populations, no differences in the degree that EBV were inflated/deflated for either trait 1 or 2 were observed between nongenotyped and genotyped animals. The current study has shown that animals can be selectively genotyped in order to optimize the response to selection as a function of the cost to conduct a breeding program using single-step genomic best linear unbiased prediction.

The impact of truncating data on the predictive ability for single-step genomic best linear unbiased prediction.

Simulated and swine industry data sets were utilized to assess the impact of removing older data on the predictive ability of selection candidate estimated breeding values (EBV) when using single-step genomic best linear unbiased prediction (ssGBLUP). Simulated data included thirty replicates designed to mimic the structure of swine data sets. For the simulated data, varying amounts of data were truncated based on the number of ancestral generations back from the selection candidates. The swine data sets consisted of phenotypic and genotypic records for three traits across two breeds on animals born from 2003 to 2017. Phenotypes and genotypes were iteratively removed 1 year at a time based on the year an animal was born. For the swine data sets, correlations between corrected phenotypes (Cp) and EBV were used to evaluate the predictive ability on young animals born in 2016-2017. In the simulated data set, keeping data two generations back or greater resulted in no statistical difference (p-value > 0.05) in the reduction in the true breeding value at generation 15 compared to utilizing all available data. Across swine data sets, removing phenotypes from animals born prior to 2011 resulted in a negligible or a slight numerical increase in the correlation between Cp and EBV. Truncating data is a method to alleviate computational issues without negatively impacting the predictive ability of selection candidate EBV.

Genetic Parameter Estimates for Metabolizing Two Common Pharmaceuticals in Swine.

In livestock, the regulation of drugs used to treat livestock has received increased attention and it is currently unknown how much of the phenotypic variation in drug metabolism is due to the genetics of an animal. Therefore, the objective of the study was to determine the amount of phenotypic variation in fenbendazole and flunixin meglumine drug metabolism due to genetics. The population consisted of crossbred female and castrated male nursery pigs (n = 198) that were sired by boars represented by four breeds. The animals were spread across nine batches. Drugs were administered intravenously and blood collected a minimum of 10 times over a 48 h period. Genetic parameters for the parent drug and metabolite concentration within each drug were estimated based on pharmacokinetics (PK) parameters or concentrations across time utilizing a random regression model. The PK parameters were estimated using a non-compartmental analysis. The PK model included fixed effects of sex and breed of sire along with random sire and batch effects. The random regression model utilized Legendre polynomials and included a fixed population concentration curve, sex, and breed of sire effects along with a random sire deviation from the population curve and batch effect. The sire effect included the intercept for all models except for the fenbendazole metabolite (i.e., intercept and slope). The mean heritability across PK parameters for the fenbendazole and flunixin meglumine parent drug (metabolite) was 0.15 (0.18) and 0.31 (0.40), respectively. For the parent drug (metabolite), the mean heritability across time was 0.27 (0.60) and 0.14 (0.44) for fenbendazole and flunixin meglumine, respectively. The errors surrounding the heritability estimates for the random regression model were smaller compared to estimates obtained from PK parameters. Across both the PK and plasma drug concentration across model, a moderate heritability was estimated. The model that utilized the plasma drug concentration across time resulted in estimates with a smaller standard error compared to models that utilized PK parameters. The current study found a low to moderate proportion of the phenotypic variation in metabolizing fenbendazole and flunixin meglumine that was explained by genetics in the current study.

Invited review: Inbreeding in the genomics era: Inbreeding, inbreeding depression, and management of genomic variability.

Traditionally, pedigree-based relationship coefficients have been used to manage the inbreeding and degree of inbreeding depression that exists within a population. The widespread incorporation of genomic information in dairy cattle genetic evaluations allows for the opportunity to develop and implement methods to manage populations at the genomic level. As a result, the realized proportion of the genome that 2 individuals share can be more accurately estimated instead of using pedigree information to estimate the expected proportion of shared alleles. Furthermore, genomic information allows genome-wide relationship or inbreeding estimates to be augmented to characterize relationships for specific regions of the genome. Region-specific stretches can be used to more effectively manage areas of low genetic diversity or areas that, when homozygous, result in reduced performance across economically important traits. The use of region-specific metrics should allow breeders to more precisely manage the trade-off between the genetic value of the progeny and undesirable side effects associated with inbreeding. Methods tailored toward more effectively identifying regions affected by inbreeding and their associated use to manage the genome at the herd level, however, still need to be developed. We have reviewed topics related to inbreeding, measures of relatedness, genetic diversity and methods to manage populations at the genomic level, and we discuss future challenges related to managing populations through implementing genomic methods at the herd and population levels.

Gene co-expression network analysis identifies porcine genes associated with variation in metabolizing fenbendazole and flunixin meglumine in the liver.

Identifying individual genetic variation in drug metabolism pathways is of importance not only in livestock, but also in humans in order to provide the ultimate goal of giving the right drug at the right dose at the right time. Our objective was to identify individual genes and gene networks involved in metabolizing fenbendazole (FBZ) and flunixin meglumine (FLU) in swine liver. The population consisted of female and castrated male pigs that were sired by boars represented by 4 breeds. Progeny were randomly placed into groups: no drug (UNT), FLU or FBZ administered. Liver transcriptome profiles from 60 animals with extreme (i.e. fast or slow drug metabolism) pharmacokinetic (PK) profiles were generated from RNA sequencing. Multiple cytochrome P450 (CYP1A1, CYP2A19 and CYP2C36) genes displayed different transcript levels across treated versus UNT. Weighted gene co-expression network analysis identified 5 and 3 modules of genes correlated with PK parameters and a portion of these were enriched for biological processes relevant to drug metabolism for FBZ and FLU, respectively. Genes within identified modules were shown to have a higher transcript level relationship (i.e. connectivity) in treated versus UNT animals. Investigation into the identified genes would allow for greater insight into FBZ and FLU metabolism.

Characterization and management of long runs of homozygosity in parental nucleus lines and their associated crossbred progeny.

BACKGROUND: In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome. METHODS: We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated. RESULTS: We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (-0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions. CONCLUSIONS: We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.

Investigation of regions impacting inbreeding depression and their association with the additive genetic effect for United States and Australia Jersey dairy cattle.

BACKGROUND: Variation in environment, management practices, nutrition or selection objectives has led to a variety of different choices being made in the use of genetic material between countries. Differences in genome-level homozygosity between countries may give rise to regions that result in inbreeding depression to differ. The objective of this study was to characterize regions that have an impact on a runs of homozygosity (ROH) metric and estimate their association with the additive genetic effect of milk (MY), fat (FY) and protein yield (PY) and calving interval (CI) using Australia (AU) and United States (US) Jersey cows. METHODS: Genotyped cows with phenotypes on MY, FY and PY (n = 6751 US; n = 3974 AU) and CI (n = 5816 US; n = 3905 AU) were used in a two-stage analysis. A ROH statistic (ROH4Mb), which counts the frequency of a SNP being in a ROH of at least 4 Mb was calculated across the genome. In the first stage, residuals were obtained from a model that accounted for the portion explained by the estimated breeding value. In the second stage, these residuals were regressed on ROH4Mb using a single marker regression model and a gradient boosted machine (GBM) algorithm. The relationship between the additive and ROH4Mb of a region was characterized based on the (co)variance of 500 kb estimated genomic breeding values derived from a Bayesian LASSO analysis. Phenotypes to determine ROH4Mb and additive effects were residuals from the two-stage approach and yield deviations, respectively. RESULTS: Associations between yield traits and ROH4Mb were found for regions on BTA13, BTA23 and BTA25 for the US population and BTA3, BTA7, BTA17 for the AU population. Only one association (BTA7) was found for CI and ROH4Mb for the US population. Multiple potential epistatic interactions were characterized based on the GBM analysis. Lastly, the covariance sign between ROH4Mb and additive SNP effect of a region was heterogeneous across the genome. CONCLUSION: We identified multiple genomic regions associated with ROH4Mb in US and AU Jersey females. The covariance of regions impacting ROH4Mb and the additive genetic effect were positive and negative, which provides evidence that the homozygosity effect is location dependent.

Differential Gene Expression across Breed and Sex in Commercial Pigs Administered Fenbendazole and Flunixin Meglumine.

Characterizing the variability in transcript levels across breeds and sex in swine for genes that play a role in drug metabolism may shed light on breed and sex differences in drug metabolism. The objective of the study is to determine if there is heterogeneity between swine breeds and sex in transcript levels for genes previously shown to play a role in drug metabolism for animals administered flunixin meglumine or fenbendazole. Crossbred nursery female and castrated male pigs (n = 169) spread across 5 groups were utilized. Sires (n = 15) of the pigs were purebred Duroc, Landrace, Yorkshire or Hampshire boars mated to a common sow population. Animals were randomly placed into the following treatments: no drug (control), flunixin meglumine, or fenbendazole. One hour after the second dosing, animals were sacrificed and liver samples collected. Quantitative Real-Time PCR was used to measure liver gene expression of the following genes: SULT1A1, ABCB1, CYP1A2, CYP2E1, CYP3A22 and CYP3A29. The control animals were used to investigate baseline transcript level differences across breed and sex. Post drug administration transcript differences across breed and sex were investigated by comparing animals administered the drug to the controls. Contrasts to determine fold change were constructed from a model that included fixed and random effects within each drug. Significant (P-value <0.007) basal transcript differences were found across breeds for SULT1A1, CYP3A29 and CYP3A22. Across drugs, significant (P-value <0.0038) transcript differences existed between animals given a drug and controls across breeds and sex for ABCB1, PS and CYP1A2. Significant (P <0.0038) transcript differences across breeds were found for CYP2E1 and SULT1A1 for flunixin meglumine and fenbendazole, respectively. The current analysis found transcript level differences across swine breeds and sex for multiple genes, which provides greater insight into the relationship between flunixin meglumine and fenbendazole and known drug metabolizing genes.

Genome-wide association study on legendre random regression coefficients for the growth and feed intake trajectory on Duroc Boars.

BACKGROUND: Feed intake and growth are economically important traits in swine production. Previous genome wide association studies (GWAS) have utilized average daily gain or daily feed intake to identify regions that impact growth and feed intake across time. The use of longitudinal models in GWAS studies, such as random regression, allows for SNPs having a heterogeneous effect across the trajectory to be characterized. The objective of this study is therefore to conduct a single step GWAS (ssGWAS) on the animal polynomial coefficients for feed intake and growth. RESULTS: Corrected daily feed intake (DFI Adj) and average daily weight measurements (DBW Avg) on 8981 (n=525,240 observations) and 5643 (n=283,607 observations) animals were utilized in a random regression model using Legendre polynomials (order=2) and a relationship matrix that included genotyped and un-genotyped animals. A ssGWAS was conducted on the animal polynomials coefficients (intercept, linear and quadratic) for animals with genotypes (DFIAdj: n=855; DBWAvg: n=590). Regions were characterized based on the variance of 10-SNP sliding windows GEBV (WGEBV). A bootstrap analysis (n=1000) was conducted to declare significance. Heritability estimates for the traits trajectory ranged from 0.34-0.52 to 0.07-0.23 for DBWAvg and DFIAdj, respectively. Genetic correlations across age classes were large and positive for both DBWAvg and DFIAdj, albeit age classes at the beginning had a small to moderate genetic correlation with age classes towards the end of the trajectory for both traits. The WGEBV variance explained by significant regions (P<0.001) for each polynomial coefficient ranged from 0.2-0.9 to 0.3-1.01% for DBWAvg and DFIAdj, respectively. The WGEBV variance explained by significant regions for the trajectory was 1.54 and 1.95% for DBWAvg and DFIAdj. Both traits identified candidate genes with functions related to metabolite and energy homeostasis, glucose and insulin signaling and behavior. CONCLUSIONS: We have identified regions of the genome that have an impact on the intercept, linear and quadratic terms for DBWAvg and DFIAdj. These results provide preliminary evidence that individual growth and feed intake trajectories are impacted by different regions of the genome at different times.

Characterizing homozygosity across United States, New Zealand and Australian Jersey cow and bull populations.

BACKGROUND: Dairy cattle breeding objectives are in general similar across countries, but environment and management conditions may vary, giving rise to slightly different selection pressures applied to a given trait. This potentially leads to different selection pressures to loci across the genome that, if large enough, may give rise to differential regions with high levels of homozygosity. The objective of this study was to characterize differences and similarities in the location and frequency of homozygosity related measures of Jersey dairy cows and bulls from the United States (US), Australia (AU) and New Zealand (NZ). RESULTS: The populations consisted of a subset of genotyped Jersey cows born in US (n = 1047) and AU (n = 886) and Jersey bulls progeny tested from the US (n = 736), AU (n = 306) and NZ (n = 768). Differences and similarities across populations were characterized using a principal component analysis (PCA) and a run of homozygosity (ROH) statistic (ROH45), which counts the frequency of a single nucleotide polymorphism (SNP) being in a ROH of at least 45 SNP. Regions that exhibited high frequencies of ROH45 and those that had significantly different ROH45 frequencies between populations were investigated for their association with milk yield traits. Within sex, the PCA revealed slight differentiation between the populations, with the greatest occurring between the US and NZ bulls. Regions with high levels of ROH45 for all populations were detected on BTA3 and BTA7 while several other regions differed in ROH45 frequency across populations, the largest number occurring for the US and NZ bull contrast. In addition, multiple regions with different ROH45 frequencies across populations were found to be associated with milk yield traits. CONCLUSION: Multiple regions exhibited differential ROH45 across AU, NZ and US cow and bull populations, an interpretation is that locations of the genome are undergoing differential directional selection. Two regions on BTA3 and BTA7 had high ROH45 frequencies across all populations and will be investigated further to determine the gene(s) undergoing directional selection.

Beef cattle body temperature during climatic stress: a genome-wide association study.

Cattle are reared in diverse environments and collecting phenotypic body temperature (BT) measurements to characterize BT variation across diverse environments is difficult and expensive. To better understand the genetic basis of BT regulation, a genome-wide association study was conducted utilizing crossbred steers and heifers totaling 239 animals of unknown pedigree and breed fraction. During predicted extreme heat and cold stress events, hourly tympanic and vaginal BT devices were placed in steers and heifers, respectively. Individuals were genotyped with the BovineSNP50K_v2 assay and data analyzed using Bayesian models for area under the curve (AUC), a measure of BT over time, using hourly BT observations summed across 5-days (AUC summer 5-day (AUCS5D) and AUC winter 5-day (AUCW5D)). Posterior heritability estimates were moderate to high and were estimated to be 0.68 and 0.21 for AUCS5D and AUCW5D, respectively. Moderately positive correlations between direct genomic values for AUCS5D and AUCW5D (0.40) were found, although a small percentage of the top 5% 1-Mb windows were in common. Different sets of genes were associated with BT during winter and summer, thus simultaneous selection for animals tolerant to both heat and cold appears possible.