Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity.

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Symposium review: Exploiting homozygosity in the era of genomics-Selection, inbreeding, and mating programs.

The advent of genomic selection paved the way for an unprecedented acceleration in genetic progress. The increased ability to select superior individuals has been coupled with a drastic reduction in the generation interval for most dairy populations, representing both an opportunity and a challenge. Homozygosity is now rapidly accumulating in dairy populations. Currently, inbreeding depression is managed mostly by culling at the farm level and by controlling the overall accumulation of homozygosity at the population level. A better understanding of how homozygosity and recessive load are related will guarantee continued genetic improvement while curtailing the accumulation of harmful recessives and maintaining enough genetic variability to ensure the possibility of selection in the face of changing environmental conditions. In this review, we present a snapshot of the current dairy selection structure as it relates to response to selection and accumulation of homozygosity, briefly outline the main approaches currently used to manage inbreeding and overall variability, and present some approaches that can be used in the short term to control accumulation of harmful recessives while maintaining sustained selection pressure.

Analysis of milk leukocyte differential measures for use in management practices for decreased mastitis incidence.

The aim of this study was to assess the usefulness of measures derived from milk leukocyte differential (MLD) in practices that improve fresh cow mastitis monitoring and decrease mastitis incidence. Quarter milk samples were collected from Holstein and Jersey cows on d 4 and 11 postcalving. Samples were analyzed using MLD, whereby cell counts and quarter infection diagnosis were obtained. Measures derived from MLD included cell scores (total leukocyte, neutrophil, macrophage, and lymphocyte scores), cell proportions (neutrophil, macrophage, and lymphocyte percentages), cell thresholds (total leukocyte, neutrophil, macrophage, and lymphocyte thresholds), and MLD diagnosis at different threshold settings (A, B, and C). Microbiological culturing of milk samples was used to determine infection status to compare the MLD diagnosis and serve as an indicator of infection. Measures derived from the microbiological analysis included occurrence of major pathogens, minor pathogens, and infection. Data analysis was based on a linear mixed model, which was used on all measures for the estimation of the fixed effects of breed, lactation number, day of sample collection, time of sampling, and quarter location, and the random effects of animal and week of sampling. All the fixed effects studied were significant for one or more of the analyzed measures. The results of this study showed that MLD-derived measures justify further study on their use for management practices for mastitis screening and prevention in early lactation.

A comparison of accuracy validation methods for genomic and pedigree-based predictions of swine litter size traits using Large White and simulated data.

The objective of this study was to compare and determine the optimal validation method when comparing accuracy from single-step GBLUP (ssGBLUP) to traditional pedigree-based BLUP. Field data included six litter size traits. Simulated data included ten replicates designed to mimic the field data in order to determine the method that was closest to the true accuracy. Data were split into training and validation sets. The methods used were as follows: (i) theoretical accuracy derived from the prediction error variance (PEV) of the direct inverse (iLHS), (ii) approximated accuracies from the accf90(GS) program in the BLUPF90 family of programs (Approx), (iii) correlation between predictions and the single-step GEBVs from the full data set (GEBVFull ), (iv) correlation between predictions and the corrected phenotypes of females from the full data set (Yc ), (v) correlation from method iv divided by the square root of the heritability (Ych ) and (vi) correlation between sire predictions and the average of their daughters' corrected phenotypes (Ycs ). Accuracies from iLHS increased from 0.27 to 0.37 (37%) in the Large White. Approximation accuracies were very consistent and close in absolute value (0.41 to 0.43). Both iLHS and Approx were much less variable than the corrected phenotype methods (ranging from 0.04 to 0.27). On average, simulated data showed an increase in accuracy from 0.34 to 0.44 (29%) using ssGBLUP. Both iLHS and Ych approximated the increase well, 0.30 to 0.46 and 0.36 to 0.45, respectively. GEBVFull performed poorly in both data sets and is not recommended. Results suggest that for within-breed selection, theoretical accuracy using PEV was consistent and accurate. When direct inversion is infeasible to get the PEV, correlating predictions to the corrected phenotypes divided by the square root of heritability is adequate given a large enough validation data set.

A heuristic method to identify runs of homozygosity associated with reduced performance in livestock.

Although, for the most part, genome-wide metrics are currently used in managing livestock inbreeding, genomic data offer, in principle, the ability to identify functional inbreeding. Here, we present a heuristic method to identify haplotypes contained within a run of homozygosity (ROH) associated with reduced performance. Results are presented for simulated and swine data. The algorithm comprises 3 steps. Step 1 scans the genome based on marker windows of decreasing size and identifies ROH genotypes associated with an unfavorable phenotype. Within this stage, multiple aggregation steps reduce the haplotype to the smallest possible length. In step 2, the resulting regions are formally tested for significance with the use of a linear mixed model. Lastly, step 3 removes nested windows. The effect of the unfavorable haplotypes identified and their associated haplotype probabilities for a progeny of a given mating pair or an individual can be used to generate an inbreeding load matrix (ILM). Diagonals of ILM characterize the functional individual inbreeding load (IIL). We estimated the accuracy of predicting the phenotype based on IIL. We further compared the significance of the regression coefficient for IIL on phenotypes with genome-wide inbreeding metrics. We tested the algorithm using simulated scenarios (12 scenarios), combining different levels of linkage disequilibrium (LD) and number of loci impacting a quantitative trait. Additionally, we investigated 9 traits from 2 maternal purebred swine lines. In simulated data, as the LD in the population increased, the algorithm identified a greater proportion of the true unfavorable ROH effects. For example, the proportion of highly unfavorable true ROH effects identified rose from 32 to 41% for the low- to the high-LD scenario. In both simulated and real data, the haplotypes identified were contained within a much larger ROH (9.12-12.1 Mb). The IIL prediction accuracy was greater than 0 across all scenarios for simulated data (mean of 0.49 [95% confidence interval 0.47-0.52] for the high-LD scenario) and for nearly all swine traits (mean of 0.17 [SD 0.10]). On average, across simulated and swine data sets, the IIL regression coefficient was more closely related to progeny performance than any genome-wide inbreeding metric. A heuristic method was developed that identified ROH genotypes with reduced performance and characterized the combined effects of ROH genotypes within and across individuals.

The relationship between different measures of feed efficiency and feeding behavior traits in Duroc pigs.

Utilization of feed in livestock species consists of a wide range of biological processes, and therefore, its efficiency can be expressed in various ways, including direct measurement, such as daily feed intake, as well as indicator measures, such as feeding behavior. Measuring feed efficiency is important to the swine industry, and its accuracy can be enhanced by using automated feeding systems, which record feed intake and associated feeding behavior of individual animals. Each automated feeder space is often shared among several pigs and therefore raises concerns about social interactions among pen mates with regard to feeding behavior. The study herein used a data set of 14,901 Duroc boars with individual records on feed intake, feeding behavior, and other off-test traits. These traits were modeled with and without the random spatial effect of Pen_Room, a concatenation of room and pen, or random social interaction among pen mates. The nonheritable spatial effect of common Pen-Room was observed for traits directly measuring feed intake and accounted for up to 13% of the total phenotypic variance in the average daily feeding rate. The social interaction effect explained larger proportions of phenotypic variation in all the traits studied, with the highest being 59% for ADFI in the group of feeding behaviors, 73% for residual feed intake (RFI; RFI4 and RFI6) in the feed efficiency traits, and 69% for intramuscular fat percentage in the off-test traits. After accounting for the social interaction effect, residual BW gain and RFI and BW gain (RIG) were found to have the heritability of 0.38 and 0.18, respectively, and had strong genetic correlations with growth and off-test traits. Feeding behavior traits were found to be moderately heritable, ranging from 0.14 (ADFI) to 0.52 (average daily occupation time), and some of them were strongly correlated with feed efficiency measures; for example, there was a genetic correlation of 0.88 between ADFI and RFI6. Our work suggested that accounting for the social common pen effect was important for estimating genetic parameters of traits recorded by the automated feeding system. Residual BW gain and RIG appeared to be two robust measures of feed efficiency. Feeding behavior measures are worth further investigation as indicators of feed efficiency.

Geno-Diver: A combined coalescence and forward-in-time simulator for populations undergoing selection for complex traits.

Geno-Diver is a combined coalescence and forward-in-time simulator designed to simulate complex traits with a quantitative and/or fitness component and implement multiple selection and mating strategies utilizing pedigree or genomic information. The simulation is carried out in two steps. The first step generates whole-genome sequence data for founder individuals. A variety of trait architectures can be generated for quantitative and fitness traits along with their covariance. The second step generates new individuals forward-in-time based on a variety of selection and mating scenarios. Genetic values are predicted for individuals utilizing pedigree or genomic information. Relationship matrices and their associated inverses are generated using computationally efficient routines. We benchmarked Geno-Diver with a previous simulation program and described how to simulate a traditional quantitative trait along with a quantitative and fitness trait. A user manual with examples, source code in C++11 and executable versions of Geno-Diver for Linux are freely available at https://github.com/jeremyhoward/Geno-Diver.

Genotype by environment (climate) interaction improves genomic prediction for production traits in US Holstein cattle.

Genotype by environment interaction (G × E) in dairy cattle productive traits has been shown to exist, but current genetic evaluation methods do not take this component into account. As several environmental descriptors (e.g., climate, farming system) are known to vary within the United States, not accounting for the G × E could lead to reranking of bulls and loss in genetic gain. Using test-day records on milk yield, somatic cell score, fat, and protein percentage from all over the United States, we computed within herd-year-season daughter yield deviations for 1,087 Holstein bulls and regressed them on genetic and environmental information to estimate variance components and to assess prediction accuracy. Genomic information was obtained from a 50k SNP marker panel. Environmental effect inputs included herd (160 levels), geographical region (7 levels), geographical location (2 variables), climate information (7 variables), and management conditions of the herds (16 total variables divided in 4 subgroups). For each set of environmental descriptors, environmental, genomic, and G × E components were sequentially fitted. Variance components estimates confirmed the presence of G × E on milk yield, with its effect being larger than main genetic effect and the environmental effect for some models. Conversely, G × E was moderate for somatic cell score and small for milk composition. Genotype by environment interaction, when included, partially eroded the genomic effect (as compared with the models where G × E was not included), suggesting that the genomic variance could at least in part be attributed to G × E not appropriately accounted for. Model predictive ability was assessed using 3 cross-validation schemes (new bulls, incomplete progeny test, and new environmental conditions), and performance was compared with a reference model including only the main genomic effect. In each scenario, at least 1 of the models including G × E was able to perform better than the reference model, although it was not possible to find the overall best-performing model that included the same set of environmental descriptors. In general, the methodology used is promising in accounting for G × E in genomic predictions, but challenges exist in identifying a unique set of covariates capable of describing the entire variety of environments.

The use of multiple imputation for the accurate measurements of individual feed intake by electronic feeders.

Obtaining accurate individual feed intake records is the key first step in achieving genetic progress toward more efficient nutrient utilization in pigs. Feed intake records collected by electronic feeding systems contain errors (erroneous and abnormal values exceeding certain cutoff criteria), which are due to feeder malfunction or animal-feeder interaction. In this study, we examined the use of a novel data-editing strategy involving multiple imputation to minimize the impact of errors and missing values on the quality of feed intake data collected by an electronic feeding system. Accuracy of feed intake data adjustment obtained from the conventional linear mixed model (LMM) approach was compared with 2 alternative implementations of multiple imputation by chained equation, denoted as MI (multiple imputation) and MICE (multiple imputation by chained equation). The 3 methods were compared under 3 scenarios, where 5, 10, and 20% feed intake error rates were simulated. Each of the scenarios was replicated 5 times. Accuracy of the alternative error adjustment was measured as the correlation between the true daily feed intake (DFI; daily feed intake in the testing period) or true ADFI (the mean DFI across testing period) and the adjusted DFI or adjusted ADFI. In the editing process, error cutoff criteria are used to define if a feed intake visit contains errors. To investigate the possibility that the error cutoff criteria may affect any of the 3 methods, the simulation was repeated with 2 alternative error cutoff values. Multiple imputation methods outperformed the LMM approach in all scenarios with mean accuracies of 96.7, 93.5, and 90.2% obtained with MI and 96.8, 94.4, and 90.1% obtained with MICE compared with 91.0, 82.6, and 68.7% using LMM for DFI. Similar results were obtained for ADFI. Furthermore, multiple imputation methods consistently performed better than LMM regardless of the cutoff criteria applied to define errors. In conclusion, multiple imputation is proposed as a more accurate and flexible method for error adjustments in feed intake data collected by electronic feeders.

Benchmarking dairy herd health status using routinely recorded herd summary data.

Genetic improvement of dairy cattle health through the use of producer-recorded data has been determined to be feasible. Low estimated heritabilities indicate that genetic progress will be slow. Variation observed in lowly heritable traits can largely be attributed to nongenetic factors, such as the environment. More rapid improvement of dairy cattle health may be attainable if herd health programs incorporate environmental and managerial aspects. More than 1,100 herd characteristics are regularly recorded on farm test-days. We combined these data with producer-recorded health event data, and parametric and nonparametric models were used to benchmark herd and cow health status. Health events were grouped into 3 categories for analyses: mastitis, reproductive, and metabolic. Both herd incidence and individual incidence were used as dependent variables. Models implemented included stepwise logistic regression, support vector machines, and random forests. At both the herd and individual levels, random forest models attained the highest accuracy for predicting health status in all health event categories when evaluated with 10-fold cross-validation. Accuracy (SD) ranged from 0.61 (0.04) to 0.63 (0.04) when using random forest models at the herd level. Accuracy of prediction (SD) at the individual cow level ranged from 0.87 (0.06) to 0.93 (0.001) with random forest models. Highly significant variables and key words from logistic regression and random forest models were also investigated. All models identified several of the same key factors for each health event category, including movement out of the herd, size of the herd, and weather-related variables. We concluded that benchmarking health status using routinely collected herd data is feasible. Nonparametric models were better suited to handle this complex data with numerous variables. These data mining techniques were able to perform prediction of health status and could add evidence to personal experience in herd management.

Variance component estimates for alternative litter size traits in swine.

Litter size at d 5 (LS5) has been shown to be an effective trait to increase total number born (TNB) while simultaneously decreasing preweaning mortality. The objective of this study was to determine the optimal litter size day for selection (i.e., other than d 5). Traits included TNB, number born alive (NBA), litter size at d 2, 5, 10, 30 (LS2, LS5, LS10, LS30, respectively), litter size at weaning (LSW), number weaned (NW), piglet mortality at d 30 (MortD30), and average piglet birth weight (BirthWt). Litter size traits were assigned to biological litters and treated as a trait of the sow. In contrast, NW was the number of piglets weaned by the nurse dam. Bivariate animal models included farm, year-season, and parity as fixed effects. Number born alive was fit as a covariate for BirthWt. Random effects included additive genetics and the permanent environment of the sow. Variance components were plotted for TNB, NBA, and LS2 to LS30 using univariate animal models to determine how variances changed over time. Additive genetic variance was minimized at d 7 in Large White and at d 14 in Landrace pigs. Total phenotypic variance for litter size traits decreased over the first 10 d and then stabilized. Heritability estimates increased between TNB and LS30. Genetic correlations between TNB, NBA, and LS2 to LS29 with LS30 plateaued within the first 10 d. A genetic correlation with LS30 of 0.95 was reached at d 4 for Large White and at d 8 for Landrace pigs. Heritability estimates ranged from 0.07 to 0.13 for litter size traits and MortD30. Birth weight had an h of 0.24 and 0.26 for Large White and Landrace pigs, respectively. Genetic correlations among LS30, LSW, and NW ranged from 0.97 to 1.00. In the Large White breed, genetic correlations between MortD30 with TNB and LS30 were 0.23 and -0.64, respectively. These correlations were 0.10 and -0.61 in the Landrace breed. A high genetic correlation of 0.98 and 0.97 was observed between LS10 and NW for Large White and Landrace breeds, respectively. This would indicate that NW could possibly be used as an effective maternal trait, given a low level of cross-fostering, to avoid back calculating litter size traits from piglet records. Litter size at d 10 would be a compromise between gain in litter size at weaning and minimizing the potentially negative effects of the nurse dam and direct additive genetics of the piglets, as they are expected to increase throughout lactation.

Short communication: Genomic selection for hoof lesions in first-parity US Holsteins.

Hoof lesions contributing to lameness are crucial economic factors that hinder the profitability of dairy enterprises. Producer-recorded hoof lesions data of US Holsteins were categorized into infectious (abscess, digital and interdigital dermatitis, heel erosion, and foot rot) and noninfectious (korn, corkscrew, sole and toe ulcer, sole hemorrhage, white line separation, fissures, thin soles, and upper leg lesions) categories of hoof lesions. Pedigree- and genomic-based univariate analyses were conducted to estimate the variance components and heritability of infectious and noninfectious hoof lesions. A threshold sire model was used with fixed effects of year-seasons and random effects of herd and sire. For genomic-based analysis, a single-step procedure was conducted, incorporating H matrix to estimate genomic variance components and heritability for hoof lesions. The pedigree-based analysis produced heritability estimates of 0.11 (±0.05) for infectious hoof lesions and 0.08 (±0.05) for noninfectious hoof lesions. The single-step genomic analysis produced heritability estimates of 0.14 (±0.06) for infectious hoof lesions and 0.12 (±0.08) for noninfectious hoof lesions. Approximated genetic correlations between hoof lesion traits and hoof type traits along with productive life and net merit were all low and ranged between -0.25 and 0.14. Sire reliabilities increased, on average, by 0.24 and 0.18 for infectious and noninfectious hoof lesions, respectively, with incorporation of genomic data.

Inferring causal relationships between reproductive and metabolic health disorders and production traits in first-lactation US Holsteins using recursive models.

Health disorders in dairy cows have a substantial effect on the profitability of a dairy enterprise because of loss in milk sales, culling of unhealthy cows, and replacement costs. Complex relationships exist between health disorders and production traits. Understanding the causal structures among these traits may help us disentangle these complex relationships. The principal objective of this study was to use producer-recorded data to explore phenotypic and genetic relationships among reproductive and metabolic health disorders and production traits in first-lactation US Holsteins. A total of 77,004 first-lactation daughters' records of 2,183 sires were analyzed using recursive models. Health data contained information on reproductive health disorders [retained placenta (RP); metritis (METR)] and metabolic health disorders [ketosis (KETO); displaced abomasum (DA)]. Production traits included mean milk yield (MY) from early lactation (mean MY from 6 to 60 d in milk and from 61 to 120 d in milk), peak milk yield (PMY), day in milk of peak milk yield (PeakD), and lactation persistency (LP). Three different sets of traits were analyzed in which recursive effects from each health disorder on culling, recursive effects of one health disorder on another health disorder and on MY, and recursive effects of each health disorder on production traits, including PeakD, PMY, and LP, were assumed. Different recursive Gaussian-threshold and threshold models were implemented in a Bayesian framework. Estimates of the structural coefficients obtained between health disorders and culling were positive; on the liability scale, the structural coefficients ranged from 0.929 to 1.590, confirming that the presence of a health disorder increased culling. Positive recursive effects of RP to METR (0.117) and of KETO to DA (0.122) were estimated, whereas recursive effects from health disorders to production traits were negligible in all cases. Heritability estimates of health disorders ranged from 0.023 to 0.114, in accordance with previous studies. Similarly, genetic correlations obtained between health disorders were moderate. The results obtained suggest that reproductive and metabolic health disorder and culling due to metabolic and reproductive diseases have strong causal relationships. Based on these results, we concluded that a health disorder (either reproductive or metabolic) occurring in early lactation has a moderate causal effect on the reproductive or metabolic health disorder occurring in later lactation. In addition, direct, indirect, and overall effects of reproductive and metabolic health disorders on milk yields for cows that avoid culling are weak.

Progression of Gene Expression Changes following a Mechanical Injury to Articular Cartilage as a Model of Early Stage Osteoarthritis.

An impact injury model of early stage osteoarthritis (OA) progression was developed using a mechanical insult to an articular cartilage surface to evaluate differential gene expression changes over time and treatment. Porcine patellae with intact cartilage surfaces were randomized to one of three treatments: nonimpacted control, axial impaction (2000 N), or a shear impaction (500 N axial, with tangential displacement to induce shear forces). After impact, the patellae were returned to culture for 0, 3, 7, or 14 days. At the appropriate time point, RNA was extracted from full-thickness cartilage slices at the impact site. Quantitative real-time PCR was used to evaluate differential gene expression for 18 OA related genes from four categories: cartilage matrix, degradative enzymes and inhibitors, inflammatory response and signaling, and cell apoptosis. The shear impacted specimens were compared to the axial impacted specimens and showed that shear specimens more highly expressed type I collagen (Col1a1) at the early time points. In addition, there was generally elevated expression of degradative enzymes, inflammatory response genes, and apoptosis markers at the early time points. These changes suggest that the more physiologically relevant shear loading may initially be more damaging to the cartilage and induces more repair efforts after loading.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: II. Genomewide association.

Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.

The effect of breed and sex on sulfamethazine, enrofloxacin, fenbendazole and flunixin meglumine pharmacokinetic parameters in swine.

Drug use in livestock has received increased attention due to welfare concerns and food safety. Characterizing heterogeneity in the way swine populations respond to drugs could allow for group-specific dose or drug recommendations. Our objective was to determine whether drug clearance differs across genetic backgrounds and sex for sulfamethazine, enrofloxacin, fenbendazole and flunixin meglumine. Two sires from each of four breeds were mated to a common sow population. The nursery pigs generated (n = 114) were utilized in a random crossover design. Drugs were administered intravenously and blood collected a minimum of 10 times over 48 h. A non-compartmental analysis of drug and metabolite plasma concentration vs. time profiles was performed. Within-drug and metabolite analysis of pharmacokinetic parameters included fixed effects of drug administration date, sex and breed of sire. Breed differences existed for flunixin meglumine (P-value<0.05; Cl, Vdss ) and oxfendazole (P-value<0.05, AUC0→∞ ). Sex differences existed for oxfendazole (P-value < 0.05; Tmax ) and sulfamethazine (P-value < 0.05, Cl). Differences in drug clearance were seen, and future work will determine the degree of additive genetic variation utilizing a larger population.

Genomic selection for producer-recorded health event data in US dairy cattle.

Emphasizing increased profit through increased dairy cow production has revealed a negative relationship of production with fitness and health traits. Decreased cow health can affect herd profitability through increased rates of involuntary culling and decreased or lost milk sales. The development of genomic selection methodologies, with accompanying substantial gains in reliability for low-heritability traits, may dramatically improve the feasibility of genetic improvement of dairy cow health. Producer-recorded health information may provide a wealth of information for improvement of dairy cow health, thus improving profitability. The principal objective of this study was to use health data collected from on-farm computer systems in the United States to estimate variance components and heritability for health traits commonly experienced by dairy cows. A single-step analysis was conducted to estimate genomic variance components and heritabilities for health events, including cystic ovaries, displaced abomasum, ketosis, lameness, mastitis, metritis, and retained placenta. A blended H matrix was constructed for a threshold model with fixed effects of parity and year-season and random effects of herd-year and sire. The single-step genomic analysis produced heritability estimates that ranged from 0.02 (standard deviation = 0.005) for lameness to 0.36 (standard deviation = 0.08) for retained placenta. Significant genetic correlations were found between lameness and cystic ovaries, displaced abomasum and ketosis, displaced abomasum and metritis, and retained placenta and metritis. Sire reliabilities increased, on average, approximately 30% with the incorporation of genomic data. From the results of these analyses, it was concluded that genetic selection for health traits using producer-recorded data are feasible in the United States, and that the inclusion of genomic data substantially improves reliabilities for these traits.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: I. Genetic parameter estimation and accuracy of genomic prediction.

The efficiency of producing salable products in the pork industry is largely determined by costs associated with feed and by the amount and quality of lean meat produced. The objectives of this paper were 1) to explore heritability and genetic correlations for growth, feed efficiency, and real-time ultrasound traits using both pedigree and marker information and 2) to assess accuracy of genomic prediction for those traits using Bayes A prediction models in a Duroc terminal sire population. Body weight at birth (BW at birth) and weaning (BW at weaning) and real-time ultrasound traits, including back fat thickness (BF), muscle depth (MD), and intramuscular fat content (IMF), were collected on the basis of farm protocol. Individual feed intake and serial BW records of 1,563 boars obtained from feed intake recording equipment (FIRE; Osborne Industries Inc., Osborne, KS) were edited to obtain growth, feed intake, and feed efficiency traits, including ADG, ADFI, feed conversion ratio (FCR), and residual feed intake (RFI). Correspondingly, 1,047 boars were genotyped using the Illumina PorcineSNP60 BeadChip. The remaining 516 boars, as an independent sample, were genotyped with a low-density GGP-Porcine BeadChip and imputed to 60K. Magnitudes of heritability from pedigree analysis were moderate for growth, feed intake, and ultrasound traits (ranging from 0.44 ± 0.11 for ADG to 0.58 ± 0.09 for BF); heritability estimates were 0.32 ± 0.09 for FCR but only 0.10 ± 0.05 for RFI. Comparatively, heritability estimates using marker information by Bayes A models were about half of those from pedigree analysis, suggesting "missing heritability." Moderate positive genetic correlations between growth and feed intake (0.32 ± 0.05) and back fat (0.22 ± 0.04), as well as negative genetic correlations between growth and feed efficiency traits (-0.21 ± 0.08, -0.05 ± 0.07), indicate selection solely on growth traits may lead to an undesirable increase in feed intake, back fat, and reduced feed efficiency. Genetic correlations among growth, feed intake, and FCR assessed by a multiple-trait Bayes A model resulted in increased genetic correlation between ADG and ADFI, a negative correlation between ADFI and FCR, and a positive correlation between ADG and FCR. Accuracies of genomic prediction for the traits investigated, ranging from 9.4% for RFI to 36.5% for BF, were reported that might provide new insight into pig breeding and future selection programs using genomic information.

Short communication: Genetic analysis of dairy bull fertility from field data of Brown Swiss cattle.

The aim of this study was to estimate heritability and repeatability of dairy bull fertility in Italian Brown Swiss cattle. Bull fertility indicators were calving per service and nonreturn rate at 56 d after service. Data included 124,206 inseminations carried out by 86 technicians on 28,873 heifers and cows in 1,400 herds. Services were recorded from 1999 to 2008 and were performed with semen from 306 AI Brown Swiss bulls. Data were analyzed with a threshold animal model, which included the fixed effects of parity by class of days in milk of the inseminated cow (age at insemination for heifers), year-season of insemination, and status of the service bull at the time of insemination (progeny testing or proven), and the random effects of herd, technician, additive genetic, and permanent environment of inseminated heifer/cow and service bull, and residual. Also, genetic covariance between heifer/cow and service bull effects was considered in the model. Heritability and repeatability were 0.0079 and 0.0100 for nonreturn rate at 56 d after service, and 0.0153 and 0.0202 for calving per service, respectively. The low estimates obtained in the present study indicate that selection for male fertility using field data is hardly pursuable.