Variation in protein metabolism biomarkers during the transition period and associations with health, colostrum quality, reproduction, and milk production traits in Holstein cows

The aims of this study were to assess (1) the variation of protein metabolism biomarkers and factors affecting them during the transition period, (2) the association of each biomarker with skeletal muscle reserves and their changes, and (3) the association of these biomarkers with postpartum health, colostrum quality, reproduction, and milk production. For this purpose, 238 multiparous Holstein cows from 6 herds were used in a prospective cohort study. Plasma concentrations of 3-methylhistidine (3-MH) and 1-methylhistidine (1-MH) and serum concentrations of total protein (TP), albumin (ALB), urea nitrogen (BUN), and creatinine (SCR) were determined for each cow at -21, -7, 7, 21, and 28 d relative to calving. Clinical diseases were recorded during the first 28 d postcalving, and presence of subclinical ketosis (scKET) was investigated at 7 and 21 d. Colostrum quality was estimated by Brix refractometry. Reproduction data by 150 d in milk (DIM) and milk production records were also available. Linear mixed models including the fixed effects of time point, herd, parity, body condition score (-21 d), duration of dry period and postparturient diseases were fitted to assess the variation in each biomarker's concentration. The association between the biomarkers' concentration during the prepartum period with the odds for each postparturient disease and for a combined trait (CD_1-28), defined as the presence of at least one clinical condition during the first 28 d after calving, were assessed with separate binary logistic models for time points -21 d and -7 d. The relationship of each biomarker's concentration with longissimus dorsi thickness (LDT) and the changes in LDT (ΔLDT) was assessed with pairwise correlations. Separate general linear models were used to assess the association of each biomarker with colostrum Brix values and milk production traits. Finally, the associated hazard for first artificial insemination (AI) and for pregnancy by 150 DIM (PREG_150DIM) was assessed with Cox proportional hazard models, whereas odds for pregnancy to the first AI (PREG_1stAI) were assessed with binary logistic models. The level of 3-MH was affected mainly by herd, time points, and their interaction. Higher 3-MH was associated with increased odds for metritis and CD_1-28, increased hazard for PREG_150 DIM and with increased milk production. 1-Methylhistidine was affected mainly by herd, scKET and occurrence of displaced abomasum. Higher 1-MH was associated with better colostrum quality, increased odds for scKET, increased hazard for first AI by 150 DIM and with decreased milk production. Both 3-MH and 1-MH were weakly to moderately negatively correlated with LDT and moderately to strongly negatively correlated with ΔLDT at the corresponding time periods. Additionally, higher TP was associated with increased odds for metritis and CD_1-28 and increased milk production, while higher ALB was associated with increased odds for scKET and increased milk production. Moreover, higher BUN was associated with decreased odds for scKET, increased odds for PREG_1stAI and increased milk production. Higher SCR was associated with decreased odds for retained fetal membranes, metritis, and CD_1-28. Periparturient protein metabolism is significantly associated with postpartum health, colostrum quality, reproduction, and milk production; mechanisms involved require further investigation.

Association between a genetic index for digital dermatitis resistance and the presence of digital dermatitis, heel horn erosion and interdigital hyperplasia in Holstein cows

Digital dermatitis (DD) is a polybacterial disease endemic to most UK dairy farms. It poses a major financial and welfare threat and is characterized by high incidence and recurrence rates. We aimed to investigate the association between the UK estimated breeding value for resistance to digital dermatitis, the Digital Dermatitis Index (DDI) and the frequency of DD, heel horn erosion (HHE), and interdigital hyperplasia (IH) in a population of Holstein dairy cows. We enrolled and genotyped 2,352 cows from 4 farms in a prospective cohort study. Foot lesion records were recorded by veterinary surgeons for each animal at 4 time points during a production cycle, starting at approximately 2 mo before calving and ending in late lactation. Importantly, these records were not used in the calculation of the DDI. Lesion records were matched to the animal's own DDI (n = 2,101) and their sire's DDI (n = 1,812). Digital Dermatitis Index values in our study population ranged from -1.41 to +1.2 and were transformed to represent distance from the mean expressed in standard deviations. The relationship between the DDI and the presence of DD was investigated using a logistic regression model, with farm, parity, and a farm-parity interaction fitted as covariates. A multivariable logistic regression model was fitted to evaluate the relationship between HHE and DDI with farm fitted as a covariate. Finally, a univariable logistic regression model with DDI as explanatory variable was used to investigate the relationship between IH and DDI. The odds ratio of an animal being affected by DD was 0.69 for one standard deviation (SD) increase in the animal's DDI (95% confidence interval (CI) = 0.63-0.76). The odds of HHE and IH were 0.69 (95%CI = 0.62-0.76) and 0.58 (95%CI = 0.49-0.68) respectively for one SD increase in DDI. The adjusted probability of DD was 32% (95% CI = 27-36%) for cows with mean DDI value of 0 while it was 24% (95% CI = 20-29%) in cows with a DDI value of +1. Sire DDI breeding values were standardized in the same way and then binned into terciles creating an ordinal variable representing bulls of high, medium, and low genetic merit for DD resistance. The daughters of low genetic merit bulls were at 2.05 (95% CI = 1.60-2.64), 1.96 (95% CI = 1.53-2.50), and 2.85 (95% CI = 1.64-5.16) times greater odds of being affected by DD, HHE, and IH respectively compared with the daughters of high genetic merit bulls. The results of this study highlight the potential of digital dermatitis genetic indexes to aid herd management of DD, and suggest that breeding for resistance to DD, alongside environmental and management control practices, could reduce the prevalence of the disease.

Genetic profile of adaptive immune traits and relationships with parasite resistance and productivity in Scottish Blackface sheep

Gastrointestinal (GI) parasites cause significant production losses in grazing ruminants which can be mitigated by breeding animals resistant to disease. Lymphocyte cytokine production and parasite-specific Immunoglobulin A (IgA) are adaptive immune traits associated with immunity to GI parasites. To explore the utility of these traits for selective breeding purposes, this study estimated the genetic parameters of the immune traits in sheep and assessed their relationship with disease and productivity traits. Whole blood stimulation assays were performed on 1 040 Scottish Blackface lambs at two months of age in 2016-2017. Blood was stimulated with either pokeweed mitogen (PWM), a non-specific activator of lymphocytes, and Teladorsagia circumcincta (T-ci) larval antigen to activate parasite-specific T lymphocytes. The type of adaptive immune response was determined by quantifying production of cytokines interferon-gamma (IFN-γ), interleukin (IL)-4, and IL-10, which relate to T-helper type (Th) 1, Th2 and regulatory T cell responses, respectively. Serum T-ci specific IgA was also quantified. Heritabilities were estimated for each immune trait by univariate analyses. Genetic and phenotypic correlations were estimated between different immune traits, and between immune traits vs. disease and productivity traits that were recorded at three months of age. Disease phenotypes were expressed as faecal egg counts (FEC) of nematode parasites (Strongyles and Nematodirus), faecal oocyst counts (FOC) of coccidian parasites, and faecal soiling score; production was measured as lamb live weight. Significant genetic variation was observed in all immune response traits. Heritabilities of cytokine production varied from low (0.14 ± 0.06) to very high (0.77 ± 0.09) and were always significantly greater than zero (P < 0.05). IgA heritability was found to be moderate (0.41 ± 0.09). Negative associations previously identified between IFN-γ production and FOC, and IL-4 production and strongyle FEC, were not evident in this study, potentially due to the time-lag between immune and parasitology measures. Instead, a positive genetic correlation was found between FOC and PWM-induced IFN-γ production, while a negative genetic correlation was found between FOC and T-ci induced IL-10. Live weight was negatively genetically correlated with IFN-γ responses. Overall, IFN-γ and IL-4 responses were positively correlated, providing little evidence of cross-regulation of Th1 and Th2 immunity within individual sheep. Furthermore, T-ci specific IgA was highly positively correlated with PWM-induced IL-10, indicating a possible role for this cytokine in IgA production. Our results suggest that while genetic selection for adaptive immune response traits is possible and may be beneficial for parasite control, selection of high IFN-γ responsiveness may negatively affect productivity.

Genome-wide association studies of parasite resistance, productivity and immunology traits in Scottish Blackface sheep

Gastrointestinal parasitism represents a global problem for grazing ruminants, which can be addressed sustainably by breeding animals to be more resistant against infection by parasites. The aim of this study was to assess the genetic architecture underlying traits associated with gastrointestinal parasite resistance, immunological profile and production in meat sheep, and identify and characterise candidate genes affecting these traits. Data on gastrointestinal parasite infection (faecal egg counts for Strongyles (FECS) and Nematodirus (FECN) and faecal oocyst counts for Coccidia, along with faecal soiling scores (DAG), characterised by the accumulation of faeces around the perineum) and production (live weight (LWT)) were gathered from a flock Scottish Blackface lambs at three and four months of age. Data on the immune profile were also collected from a subset of these lambs at two and five months of age. Immune traits included the production of Interferon-γ (IFN-γ), Interleukin (IL)-4 and IL-10 following stimulation of whole blood with pokeweed mitogen (PWM) or antigen from the gastric parasite Teladorsagia circumcincta (T-ci), and serum levels of T. circumcincta-specific immunoglobulin A (IgA). Animals were genotyped with genome-wide DNA arrays, and a total of 1 766 animals and 45 827 Single Nucleotide Polymorphisms (SNPs) were retained following quality control and imputation. Genome-wide association studies were performed for 24 traits. The effects of individual markers with significant effects were estimated, and the genotypic effect solutions were used to estimate additive and dominance effects, and the proportion of additive genetic variance attributed to each SNP locus. A total of 15 SNPs were associated at least at a suggestive level with FECS, FECN, DAG, IgA, PWM-induced IFN-γ and IL-4, and T-ci-induced IL-10. This study uncovered 52 genes closely related to immune function in proximity to these SNPs. A number of genes encoding C-type lectins and killer cell lectin-like family members were close to a SNP associated with FECN, while several genes encoding IL-1 cytokine family members were found to be associated with IgA. Potential candidate genes belonging to or in close proximity with the Major Histocompatibility Complex (MHC) were revealed, including Homeostatic Iron Regulator and butyrophilin coding genes associated with IFN-γ(PWM), and IL-17 coding genes associated with IgA. Due to the importance of the MHC in the control of immune responses, these genes may play an important role in resistance to parasitic infections. Our results reveal a largely complex and polygenic genetic profile of the studied traits in this Scottish Blackface sheep population.

Associations among post-partum rumen fill and motility, subclinical ketosis and fertility in Holstein dairy cows

This prospective observational study aimed to investigate the association of rumen fill and motility in post-partum Holstein cows with their future reproductive performance and subclinical ketosis (SCK). The study population consisted of two independent data sets: the first (DS1) included 237 cows from 6 herds and the second one (DS2) 709 cows from 9 herds. Rumen Fill Score (RFS) was transformed into a 3 level-trait, representing very low, low and adequate dry matter intake, respectively. A binary Rumen Contraction Score (RCS) was defined as: 0: <2 contractions/2 min, impaired rumen motility and 1: ≥2 contractions/2 min, normal rumen motility. A combined binary trait based on RFS and RCS (RFCS) was also established, representing unsatisfactory and satisfactory rumen function. Three SCK traits were defined, based on 3 different thresholds, SCK_I: BHB≥1,000 mmol/L, SCK_II: BHB≥1,100 mmol/L and SCK_III: BHB≥1,200 mmol/L. Scores were assessed and blood samples collected on day 7 (DS1) or day 8 (DS2), postpartum. Kaplan-Meier survival analysis, multivariable Cox proportional hazards models and Generalized Linear Mixed Models were performed to evaluate the association of rumen and SCK traits with reproduction. Herd, parity, calving season and several postparturient diseases were also included as potential explanatory variables. Mean days from calving to pregnancy after the 1st artificial insemination (AI) and from calving to pregnancy (all AIs) were shorter for levels of rumen traits representing adequate DMI and normal rumen motility; in most cases these differences were statistically significant in both datasets. Cows with adequate DMI and normal rumen motility (only in DS2) had greater hazard (hazard ratio [HR] = 1.84 and 1.61, for RFS and RFCS, respectively) and odds (odds ratio [OR] = 2.49 and 1.98, for RFS and RFCS, respectively) for pregnancy at 1st AI. Assessment of the association of examined rumen traits with hazard and odds for pregnancy at all AIs yielded statistically significant results in both datasets. For RFS, RCS and RFCS, HRs ranged from 1.57 to 3.31 and ORs from 1.95 to 4.83. No statistically significant associations with hazard and odds for pregnancy at 1st or all AIs were detected, for any of the 3 SCK traits, in either dataset. Overall, the combined RFCS trait constantly identified more than twice the number of cows with future reproductive problems than a positive SCK blood test.

Identification of temporal patterns of environmental heat stress of Holstein dairy heifers raised in Mediterranean climate during their in-utero and post-natal life periods and modelling their effects on age at first calving

A retrospective study was conducted to evaluate temporal patterns of environmental heat stress during the in-uterus period of development (IUP) and the 3-month post-natal (PN) period of dairy heifers, and to estimate their association with the age at first calving (AFC). Data from 30 dairy herds in Northern Greece including 9098 heifers were extracted from National Cattle Database. Data (2005-2019) regarding 230,100 farm-specific ambient daily temperature and relative humidity records, were obtained from ERA5-Land. Average monthly Temperature-Humidity-Index values (THI; low≤68, and high>68) were calculated and matched for each heifer to their IUP and PN. Subsequently, Cluster Analysis was used with monthly THIs as predictors to allocate heifers to THI clusters. The association of clusters with AFC was assessed with Generalized Linear Mixed Model analysis, an extended form of multiple linear regression. Finally, 8 Heat Stress Clusters (HSC; namely HSC-1 to HSC-8) were identified. Compared to HSC-8 (8th-9th IUP months and 1st PN month) heifers of HSC-5 (4th-7th IUP months) and HSC-6 (6th-8th IUP months) calved 13.8 and 17.8 days later, respectively (P < 0.01-0.001). Moreover, when AFC was treated as a binary variable, heifers of HSC-5 and HSC-6 had 1.15 and 1.34 (P < 0.01-0.001) higher risk of calving for the first time later than 787 days compared to HSC-8, respectively.

Genome-wide association study of health and production traits in meat sheep

Genotypes are currently widely used in animal breeding programmes to enhance the speed of genetic progress. With sufficient data, a Genome-Wide Association Study (GWAS) can be performed to identify informative markers. The aim of this study was to investigate the genetic background of health (footrot and mastitis) and production (birth weight, weaning weight, scan weight, and fat and muscle depth) traits using the available phenotypic and Single Nucleotide Polymorphism (SNP) data collected on the UK Texel sheep population. Initially, 10 193 genotypes were subject to quality control, leaving 9 505 genotypes for further analysis. Selected genotypes, recorded on four different Illumina chip types from low density (15 k SNPs) to high density (606 006 SNPs), were imputed to a subset of 45 686 markers from 50 k array, distributed on 27 chromosomes. Phenotypes collected on 32 farms across the UK for footrot and mastitis and extracted from the UK National database (iTexel) for the production traits were used along with pre-estimated variance components to obtain de-regressed breeding values and used to perform GWAS. Results showed three SNPs being significant on the genome-wise level ('OAR8_62240378.1' on chromosome 8 for birth weight, 's14444.1' on chromosome 19 for weaning weight and 's65197.1' on chromosome 23 for scan weight). Fourteen subsequent SNPs were found to be significant at the chromosome-wise level. These SNPs are located within or close to previously reported QTLs impacting on animal health (such as faecal egg count or somatic cell count) and production (such as body or carcass weight and fat amount). These results indicate that the studied traits are highly polygenic with complex genetic architecture.

Novel genomic markers and genes related to reproduction in prolific Chios dairy sheep: a genome-wide association study

Genetic architecture of sheep reproduction is increasingly gaining scientific interest due to the major impact on sheep production systems. In the present study, we conducted pedigree-based analyses and genome-wide association studies using the Illumina Ovine SNP50K BeadChip to explore the genetic mechanisms underlying the reproduction of the highly prolific Chios dairy sheep. First lambing age, total prolificacy and maternal lamb survival were selected as representative reproductive traits and estimated as significantly heritable (h² = 0.07-0.21) with no evident genetic antagonism among traits. We identified novel genome-wide and suggestive significant single-nucleotide polymorphisms (SNPs) on chromosomes 2 and 12 associated with age at first lambing. The new variants detected on chromosome 2 span a region of 357.79 kb with high pairwise linkage disequilibrium estimates (r² = 0.8-0.9). Functional annotation analysis revealed candidate genes, such as the collagen-type genes and the Myostatin gene, that participate in osteogenesis, myogenesis, skeletal and muscle mass development resembling the functionality of major genes affecting the ovulation rate and prolificacy. Additional functional enrichment analysis associated the collagen-type genes with multiple uterine-related disfunctions, such as cervical insufficiency, uterine prolapse and abnormalities of the uterine cervix. Several genes (e.g., KAZN, PRDM2, PDPN, LRRC28) localised close to the SNP marker on chromosome 12 were grouped in annotation enrichment clusters majorly involved in developmental and biosynthetic pathways, apoptosis, and nucleic acid-templated transcription. Our findings may further contribute to unravel the genomic regions that are important for sheep reproduction and could be incorporated into future selective breeding programmes.

The bovine foot skin microbiota is associated with host genotype and the development of infectious digital dermatitis lesions

BACKGROUND

Bovine Digital Dermatitis (BDD) is a prevalent infectious disease, causing painful foot skin lesions and lameness in cattle. We describe herein the bovine foot skin microbiota and its associations with BDD using 16S rRNA gene amplicon and shotgun metagenomic sequencing on samples from 259 dairy cows from three UK dairy farms.

RESULTS

We show evidence of dysbiosis, and differences in taxonomy and functional profiles in the bovine foot skin microbiome of clinically healthy animals that subsequently develop BDD lesions, compared to those that do not. Our results suggest that taxonomical and functional differences together with alterations in ecological interactions between bacteria in the normal foot skin microbiome may predispose an animal to develop BDD lesions. Using genome-wide association and regional heritability mapping approaches, we provide first evidence for interactions between host genotype and certain members of the foot skin microbiota. We show the existence of significant genetic variation in the relative abundance of Treponema spp. and Peptoclostridium spp. and identify regions in the bovine genome that explain a significant proportion of this variation.

CONCLUSIONS

Collectively this work shows early changes in taxonomic and functional profiles of the bovine foot-skin microbiota in clinically healthy animals which are associated with subsequent development of BDD and could be relevant to prevention of disease. The description of host genetic control of members of the foot skin microbiota, combined with the association of the latter with BDD development offer new insights into a complex relationship that can be exploited in selective breeding programmes. Video Abstract.

Skeletal muscle and adipose tissue reserves and mobilisation in transition Holstein cows: Part 1 Biological variation and affecting factors

Nutrient deficit during the periparturient period leads to mobilisation of body energy and protein reserves. Research regarding fat reserves and mobilisation is extensive, while, on the contrary, investigation of muscle mobilisation during the periparturient period is limited. The aim of this cohort study was to simultaneously investigate the biological variation of skeletal muscle and subcutaneous fat reserves together with their mobilisation in transition Holstein cows of different herds, using ultrasonography, and to assess potential affecting factors. For this purpose, ultrasound measurements of longissimus dorsi muscle thickness (LDT) and backfat thickness (BFT) from 238 multiparous cows of six dairy farms were obtained at six time points across the transition period (from 21 days pre- to 28 days postpartum). Concentrations of serum creatinine and non-esterified fatty acids were determined in order to confirm the loss of muscle mass and adipose tissue, respectively. Cases of clinical postparturient diseases and subclinical ketosis (scKET) during the first 28 days postcalving were recorded. Cows mobilised on average 32.8% and 37.3% of LDT and BFT reserves, respectively. Large between-cow variation was observed for both the onset and the degree of mobilisation. Time point, initial body condition score and parity were the most important predictors of LDT variation. Cows diagnosed with metritis (MET) had lower LDT postpartum and mobilised more muscle depth compared to cows not diagnosed with MET. Initial BCS, time point, initial BW (estimated by heart girth measurement) and parity were the most important predictors of BFT variation. Cows diagnosed with MET mobilised more backfat between -7d and 7d compared to cows not diagnosed with MET. Cows with scKET mobilised more backfat between 7- and 21 days postpartum compared to healthy ones. Variation of subcutaneous fat and skeletal muscle reserves during the transition period was large and affected by herd and several cow-level factors.

Skeletal muscle and adipose tissue reserves and mobilisation in transition Holstein cows: Part 2 association with postpartum health, reproductive performance and milk production

The aim of this study was, for the first time, to simultaneously assess the association of skeletal muscle and subcutaneous fat reserves and their mobilisation, measured by ultrasonography, with the incidence of specific postparturient health, reproduction, and milk production traits. For this purpose, ultrasound measurements of longissimus dorsi thickness (LDT) and backfat thickness (BFT) from 238 multiparous cows from 6 dairy farms were obtained at 6 time points during the transition period (from 21 days pre- to 28 days postpartum). In each case, LDT and BFT measurements at each time point and LDT and BFT mobilisation variables at each study period were assessed simultaneously. Cases of specific clinical postparturient diseases and subclinical ketosis were recorded. An additional disease trait was used, defined as the presence or absence of at least one clinical condition after calving (CD_1-28). The associated disease odds with LDT/BFT variables were assessed with binary logistic regression models. The associated hazard for 1st artificial insemination (AI) and for pregnancy by 150 days-in-milk (PREG_150DIM) was assessed with Cox proportional hazard models. Moreover, binary logistic models were used to assess the associated odds for pregnancy to 1^stAI (PREG_1^stAI). Finally, association with 30d, 100d and 305d milk yield was assessed with linear regression models. Increased muscle depth during transition was negatively associated with odds for metritis and CD_1-28, while associations with odds for subclinical ketosis were inconclusive. Moreover, increased LDT reserves were associated with greater hazard for 1st AI by 150 days-in-milk, but results were inconclusive regarding odds for PREG_1^stAI. Increased LDT mobilisation was associated with increased odds for metritis. Increased BFT reserves were positively associated with odds for metritis, CD_1-28 and subclinical ketosis and with decreased hazard for PREG_150DIM. Increased BFT mobilisation was associated with increased odds for subclinical ketosis and with decreased odds for PREG_1^stAI and decreased hazard for PREG_150DIM. Cows with moderate BFT reserves performed better. Finally, increased BFT mobilisation during -21d to -7d from parturition was associated with less milk by 30d and 100d. On the contrary, increased BFT mobilisation during -7d to 7d was associated with more milk by 305d. Metabolism of muscle and fat tissue during transition period was differently associated with different postparturient health, reproduction and milk production traits. In general, greater muscle mass and moderate fat reserves with limited muscle and fat mobilisation were associated with better performance.

Discovery of hidden pedigree errors combining genomic information with the genomic relationship matrix in Texel sheep

Genomic variants such as Single Nucleotide Polymorphisms and animal pedigree are now used widely in routine genetic evaluations of livestock in many countries. The use of genomic information not only can be used to enhance the accuracy of prediction but also to verify pedigrees for animals that are extensively managed using natural mating and enabling multiple-sire mating groups to be used. By so doing, the rate of genetic gain is enhanced, and any bias associated with incorrect pedigrees is removed. This study used a set of 8 764 sheep genotypes to verify the pedigree based on both the conventional opposing homozygote method as well as a novel method when combined with the inclusion of the genomic relationship matrix (GRM). The genomic relationship coefficients between verified pairs of animals showed on average a relationship of 0.50 with parent, 0.25 with grandparent, 0.13 with great grandparent, 0.50 with full-sibling and 0.27 with half-sibling. Minimum obtained values from these verified pairs were then used as thresholds to determine the pedigree for unverified pairs of animals, to detect potential errors in the pedigree. Using a case study from a population partially genotyped UK sheep, the results from this study illustrate a powerful way to resolve parentage inconsistencies, when combining the conventional 'opposing homozygote' method using genomic information together with GRM for pedigree checking. In this way, previously undetected pedigree errors can be resolved.

A study on the use of thermal imaging as a diagnostic tool for the detection of digital dermatitis in dairy cattle

Our aims were to (1) determine how interdigital skin temperature (IST), measured using infrared thermography, was associated with different stages of digital dermatitis (DD) lesions and (2) develop and validate models that can use IST measurements to identify cows with an active DD lesion. Between March 2019 and March 2020, infrared thermographic images of hind feet were taken from 2,334 Holstein cows across 4 farms. We recorded the maximum temperature reading from infrared thermographic images of the interdigital skin between the heel bulbs on the hind feet. Pregnant animals were enrolled approximately 1 to 2 mo precalving, reassessed 1 wk after calving, and again at approximately 50 to 100 d postpartum. At these time points, IST and the clinical stage of DD (M-stage scoring system: M1-M4.1) were recorded in addition to other data such as the ambient environmental temperature, height, body condition score, parity, and the presence of other foot lesions. A mixed effect linear regression model with IST as the dependent variable was fitted. Interdigital skin temperature was associated with DD lesions; compared to healthy feet, IST was highest in feet with M2 lesions, followed by M1 and M4.1 lesions. Subsequently, the capacity of IST measurements to detect the presence or absence of an active DD lesion (M1, M2, or M4.1) was explored by fitting logistic regression models, which were tested using 10-fold validation. A mixed effect logistic regression model with the presence of active DD as the dependent variable was fitted first. The average area under the curve for this model was 0.80 when its ability to detect presence of active DD was tested on 10% of the data that were not used for the model's training; an average sensitivity of 0.77 and an average specificity of 0.67 was achieved. This model was then restricted so that only explanatory variables that could be practically recorded in a nonresearch, external setting were included. Validation of this model demonstrated an average area under the curve of 0.78, a sensitivity of 0.88, and a specificity of 0.66 for 1 of the time points (precalving). Lower sensitivity and specificity were achieved for the other 2 time points. Our study adds further evidence to the relationship between DD and foot skin temperature using a large data set with multiple measurements per animal. Additionally, we highlight the potential for infrared thermography to be used for routine on-farm diagnosis of active DD lesions.

Subclinical hypocalcaemia follows specific time-related and severity patterns in post-partum Holstein cows

Hypocalcaemia is a common metabolic disorder of post-parturient dairy cows; enhancement of our knowledge regarding Ca dynamics would improve cow health and welfare. This study investigated the presence of time- and severity-related subclinical hypocalcaemia (SCHCa) patterns in Holstein cows after calving and their association with diseases during the first week of lactation. A total of 938 cows from nine herds were blood sampled on days 1, 2, 4 and 8 post-partum, and serum Ca concentration was measured with atomic absorption. Subclinical ketosis (serum β-hydroxybutyrate≥1.2 mmol/L) and cases of retained foetal membranes, metritis (MET), mastitis, ketosis and displaced abomasum (DA) were recorded. Using receiver operating characteristic (ROC) analysis, the SCHCa cut-off was identified at serum Ca ≤ 2.09 mmol/L. Hierarchical and two-step cluster analysis classified cows to seven clusters based on test day cow Ca records and overall SCHCa status. Two clusters (NORM and HIGH) included 318 normocalcaemic cows and five clusters (SCH-1 to -5) those that were at least once subclinically hypocalcaemic (n = 620). A second ROC analysis was performed in order to distinguish mild from severe cases of SCHCa in these 620 cows; this cut-off was identified at 1.93 mmol/L. The associated risk of disease with Ca clustership was assessed with generalized linear mixed models. Overall incidence of SCHCa was 66.1%. Clusters SCH-1 and SCH-2 included cows with short-term SCHCa of day 1 and day 2, respectively, while SCH-3 included cows with severe and prolonged SCHCa extending to day 4 and beyond. Clusters SCH-4 and SCH-5 included cows with delayed SCHCa, which culminated on days 4 and 8, respectively. Compared to NORM cows in HIGH had lower risk of MET and no cases of DA. Cows in SCH-3 had higher risk of being diagnosed with retained foetal membranes, DA or any disease during the study period. Cows in SCH-5 had higher risk of being diagnosed with ketosis, subclinical ketosis or any disease. In conclusion, there are multiple normocalcaemic and hypocalcaemic patterns that are differently associated with disease risk.

Joint genetic analysis for dairy cattle performance across countries in sub-Saharan Africa

This study assessed the feasibility of across-country genetic evaluation of dairy cattle in sub-Saharan Africa where data on livestock production are scarce. Genetic parameters were estimated for the 305-day milk yield in the first lactation and across five lactations, for age at first calving and for interval between first and second calving. Estimated breeding values of individual animals for these traits were calculated. There were records from 2 333, 25 208, and 5 929 Holstein cows in Kenya, South Africa, and Zimbabwe, and 898 and 65134 Jersey cows from Kenya and South Africa. Genetic gain from sire selection within and across countries. was predicted Genetic links between countries were determined from sires with daughters that had records in two or more countries, and from common ancestral sires across seven generations on both the maternal and paternal sides of the pedigree. Each country was treated as a trait in the across-country evaluation. The results showed that genetic variance and heritability were not always estimable within country, but were significantly different from zero in the across-country evaluation. In all three countries, there was greater genetic gain in all traits from an across-country genetic evaluation owing to greater accuracy of selection compared with within country. Kenya stood to benefit most from an across-country evaluation, followed by Zimbabwe, then South Africa. An across-country breeding programme using joint genetic evaluation would be feasible, provided that there were genetic links across countries, and would provide a platform for accelerated genetic progress through selection and germplasm exchange between sub-Saharan African countries.Keywords: across-country genetic evaluation, genetic connectedness, genetic progress

Genetic analyses of live weight and carcass composition traits in purebred Texel, Suffolk and Charollais lambs

Lamb live weight is one of the key drivers of profitability on sheep farms. Previous studies in Ireland have estimated genetic parameters for live weight and carcass composition traits using a multi-breed population rather than on an individual breed basis. The objective of the present study was to undertake genetic analyses of three lamb live weight and two carcass composition traits pertaining to purebred Texel, Suffolk and Charollais lambs born in the Republic of Ireland between 2010 and 2017, inclusive. Traits (with lamb age range in parenthesis) considered in the analyses were: pre-weaning weight (20 to 65 days), weaning weight (66 to 120 days), post-weaning weight (121 to 180 days), muscle depth (121 to 180 days) and fat depth (121 to 180 days). After data edits, 137 402 records from 50 372 lambs across 416 flocks were analysed. Variance components were derived using animal linear mixed models separately for each breed. Fixed effects included for all traits were contemporary group, age at first lambing of the dam, parity of the dam, a gender by age of the lamb interaction and a birth type by rearing type of the lamb interaction. Random effects investigated in the pre-weaning and weaning weight analyses included animal direct additive genetic, dam maternal genetic, litter common environment, dam permanent environment and residual variances. The model of analysis for post-weaning, muscle and fat depth included an animal direct additive genetic and litter common environment effect only. Significant direct additive genetic variation existed in all cases. Direct heritability for pre-weaning weight ranged from 0.14 to 0.30 across the three breeds. Weaning weight had a direct heritability ranging from 0.17 to 0.27 and post-weaning weight had a direct heritability ranging from 0.15 to 0.27. Muscle and fat depth heritability estimates ranged from 0.21 to 0.31 and 0.15 to 0.20, respectively. Positive direct correlations were evident for all traits. Results revealed ample genetic variation among animals for the studied traits and significant differences between breeds to suggest that genetic evaluations could be conducted on a per-breed basis.

Analysis of genome-wide DNA arrays reveals the genomic population structure and diversity in autochthonous Greek goat breeds

Goats play an important role in the livestock sector in Greece. The national herd consists mainly of two indigenous breeds, the Eghoria and Skopelos. Here, we report the population structure and genomic profiles of these two native goat breeds using Illumina’s Goat SNP50 BeadChip. Moreover, we present a panel of candidate markers acquired using different genetic models for breed discrimination. Quality control on the initial dataset resulted in 48,841 SNPs kept for downstream analysis. Principal component and admixture analyses were applied to assess population structure. The rate of inbreeding within breed was evaluated based on the distribution of runs of homozygosity in the genome and respective coefficients, the genomic relationship matrix, the patterns of linkage disequilibrium, and the historic effective population size. Results showed that both breeds exhibit high levels of genetic diversity. Level of inbreeding between the two breeds estimated by the Wright’s fixation index F_ST was low (Fst = 0.04362), indicating the existence of a weak genetic differentiation between them. In addition, grouping of farms according to their geographical locations was observed. This study presents for the first time a genome-based analysis on the genetic structure of the two indigenous Greek goat breeds and identifies markers that can be potentially exploited in future selective breeding programs for traceability purposes, targeted genetic improvement schemes and conservation strategies.

Genetic parameters of colostrum traits in Holstein dairy cows

The main objective of this study was to assess the genetic background of colostrum yield and quality traits after calving in Holstein dairy cows. The secondary objective was to investigate genetic and phenotypic correlations among laboratory-based and on-farm-measured colostrum traits. The study was conducted in 10 commercial dairy herds located in northern Greece. A total of 1,074 healthy Holstein cows with detailed pedigree information were examined from February 2015 to September 2016. All cows were clinically examined on the day of calving and scored for body condition. All 4 quarters were machine-milked, and a representative and composite colostrum sample was collected and examined. Colostrum total solids (TS) content was determined on-farm using a digital Brix refractometer. Colostrum fat, protein, and lactose contents were determined using an infrared milk analyzer, and energy content was calculated using National Research Council (2001) equations. Dry period length (for cows of parity ≥2), milk yield of previous 305-d lactation (for cows of parity ≥2), age at calving, parity number, season of calving, time interval between calving and first colostrum milking, and milk yield were recorded. Each trait (colostrum yield and quality traits) was analyzed with a univariate mixed model, including fixed effects of previously mentioned factors and the random animal additive genetic effect. All available pedigrees were included in the analysis, bringing the total animal number to 5,662. Estimates of (co)variance components were used to calculate heritability for each trait. Correlations among colostrum traits were estimated with bivariate analysis using the same model. Mean percentage (±SD) colostrum TS, fat, protein, and lactose contents were 25.8 ± 4.7, 6.4 ± 3.3, 17.8 ± 4.0, and 2.2 ± 0.7%, respectively; mean energy content was 1.35 ± 0.3 Mcal/kg and mean colostrum yield was 6.18 ± 3.77 kg. Heritability estimates for the above colostrum traits were 0.27, 0.21, 0.19, 0.15, 0.22, and 0.04, respectively. Several significant genetic and phenotypic correlations were derived. The genetic correlation of TS content measured on-farm with colostrum protein was practically unity, whereas the correlation with energy content was moderate (0.61). Fat content had no genetic correlation with TS content; their phenotypic correlation was positive and low. Colostrum yield was not correlated genetically with any of the other traits. In conclusion, colostrum quality traits are heritable and can be amended with genetic selection.

The genomic architecture of mastitis resistance in dairy sheep

Background

Mastitis is the most prevalent disease in dairy sheep with major economic, hygienic and welfare implications. The disease persists in all dairy sheep production systems despite the implementation of improved management practises. Selective breeding for enhanced mastitis resistance may provide the means to further control the disease. In the present study, we investigated the genetic architecture of four mastitis traits in dairy sheep. Individual animal records for clinical mastitis occurrence and three mastitis indicator traits (milk somatic cell count, total viable bacterial count in milk and the California mastitis test) were collected monthly throughout lactation for 609 ewes of the Greek Chios breed. All animals were genotyped with a custom-made 960-single nucleotide polymorphism (SNP) DNA array based on markers located in quantitative trait loci (QTL) regions for mastitis resistance previously detected in three other distinct dairy sheep populations.

Results

Heritable variation and strong positive genetic correlations were estimated for clinical mastitis occurrence and the three mastitis indicator traits. SNP markers significantly associated with these mastitis traits were confirmed on chromosomes 2, 3, 5, 16 and 19. We identified pathways, molecular interaction networks and functional gene clusters for mastitis resistance. Candidate genes within the detected regions were identified based upon analysis of an ovine transcriptional atlas and transcriptome data derived from milk somatic cells. Relevant candidate genes implicated in innate immunity included SOCS2, CTLA4, C6, C7, C9, PTGER4, DAB2, CARD6, OSMR, PLXNC1, IDH1, ICOS, FYB, and LYFR.

Conclusions

The results confirmed the presence of animal genetic variability in mastitis resistance and identified genomic regions associated with specific mastitis traits in the Chios sheep. The conserved genetic architecture of mastitis resistance between distinct dairy sheep breeds suggests that across-breed selection programmes would be feasible.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3982-1) contains supplementary material, which is available to authorized users.

Variants in the 3' untranslated region of the ovine acetyl-coenzyme A acyltransferase 2 gene are associated with dairy traits and exhibit differential allelic expression

The acetyl-CoA acyltransferase 2 (ACAA2) gene encodes an enzyme of the thiolase family that is involved in mitochondrial fatty acid elongation and degradation by catalyzing the last step of the respective β-oxidation pathway. The increased energy needs for gluconeogenesis and triglyceride synthesis during lactation are met primarily by increased fatty acid oxidation. Therefore, the ACAA2 enzyme plays an important role in the supply of energy and carbon substrates for lactation and may thus affect milk production traits. This study investigated the association of the ACAA2 gene with important sheep traits and the putative functional involvement of this gene in dairy traits. A single nucleotide substitution, a T to C transition located in the 3' untranslated region of the ACAA2 gene, was used in mixed model association analysis with milk yield, milk protein yield and percentage, milk fat yield and percentage, and litter size at birth. The single nucleotide polymorphism was significantly associated with total lactation production and milk protein percentage, with respective additive effects of 6.81 ± 2.95 kg and -0.05 ± 0.02%. Additionally, a significant dominance effect of 0.46 ± 0.21 kg was detected for milk fat yield. Homozygous TT and heterozygous CT animals exhibited higher milk yield compared with homozygous CC animals, whereas the latter exhibited increased milk protein percentage. Expression analysis from age-, lactation-, and parity-matched female sheep showed that mRNA expression of the ACAA2 gene from TT animals was 2.8- and 11.8-fold higher in liver and mammary gland, respectively. In addition, by developing an allelic expression imbalance assay, it was estimated that the T allele was expressed at an average of 18% more compared with the C allele in the udder of randomly selected ewes. We demonstrated for the first time that the variants in the 3' untranslated region of the ovine ACAA2 gene are differentially expressed in homozygous ewes of each allele and exhibit allelic expression imbalance within heterozygotes in a tissue-specific manner, supporting the existence of cis-regulatory DNA variation in the ovine ACAA2 gene. This is the first study reporting differential allelic imbalance expression of a candidate gene associated with milk production traits in dairy sheep.

Genetic analysis of carcass traits in beef cattle using random regression models

Livestock mature at different rates depending, in part, on their genetic merit; therefore, the optimal age at slaughter for progeny of certain sires may differ. The objective of the present study was to examine sire-level genetic profiles for carcass weight, carcass conformation, and carcass fat in cattle of multiple beef and dairy breeds, including crossbreeds. Slaughter records from 126,214 heifers and 124,641 steers aged between 360 and 1,200 d and from 86,089 young bulls aged between 360 and 720 d were used in the analysis; animals were from 15,127 sires. Variance components for each trait across age at slaughter were generated using sire random regression models that included quadratic polynomials for fixed and random effects; heterogeneous residual variances were assumed across ages. Heritability estimates across genders ranged from 0.08 (±0.02) to 0.34 (±0.02) for carcass weight, from 0.24 (±0.02) to 0.42 (±0.01) for conformation, and from 0.16 (±0.03) to 0.40 (±0.02) for fat score. Genetic correlations within each trait across ages weakened as the interval between ages compared lengthened but were all >0.64, suggesting a similar genetic background for each trait across different ages. Eigenvalues and eigenfunctions of the additive genetic covariance matrix revealed genetic variability among animals in their growth profiles for carcass traits, although most of the genetic variability was associated with the height of the growth profile. At the same age, a positive genetic correlation (0.60 to 0.78; SE ranged from 0.01 to 0.04) existed between carcass weight and conformation, whereas negative genetic correlations existed between fatness and both conformation (-0.46 to 0.08; SE ranged from 0.02 to 0.09) and carcass weight (-0.48 to -0.16; SE ranged from 0.02 to 0.14) at the same age. The estimated genetic parameters in the present study indicate genetic variability in the growth trajectory in cattle, which can be exploited through breeding programs and used in decision support tools.

Short communication: Identification of variation in the ovine prolactin gene of Chios sheep with a cost-effective sequence-based typing assay

The present study identified single nucleotide polymorphisms (SNP) in the coding and untranslated regions of the ovine prolactin gene of Chios sheep. By developing a cost-effective direct sequence-based typing assay, around 600 bp of reliable sequencing data and clear identification of heterozygous positions was achieved. Five SNP were found, located in exons 2 (KC764410:g.567G>A, g.625C>T, g.683C>A) and 3 (KC764410:g.2015C>A, g.2101G>A), whereas the remaining exons were monomorphic. The identified SNP were synonymous, with the exception of the g.567G>A SNP, which results in an Arg to His amino acid change. As the sequencing cost of the sequence-based typing assay was 20 orders of magnitude lower compared with a standard Sanger method, the assay was also used as a genotyping tool. The identified polymorphism was genotyped for 247 ewes and was subsequently used in mixed model association analyses of milk yield, milk fat content, and litter size at birth. The association analysis revealed a significant dominance effect of 0.17 ± 0.07 of the g.2015C>A SNP on milk fat percentage, whereas a dominance effect of -21.33 ± 10.51 of the same SNP on total lactation milk yield was also estimated. The g.2015C>A SNP explained 2.47 and 3.68% of the total phenotypic variance of milk yield and milk fat percentage, respectively, whereas the corresponding values for the animal variance were 7.14 and 11.75%. A suggestive association of the nonsynonymous g.567G>A SNP with litter size at birth was also detected.

Genetic evaluation for bovine tuberculosis resistance in dairy cattle

Genetic evaluations for resistance to bovine tuberculosis (bTB) were calculated based on British national data including individual animal tuberculin skin test results, postmortem examination (presence of bTB lesions and bacteriological culture for Mycobacterium bovis), animal movement and location information, production history, and pedigree records. Holstein cows with identified sires in herds with bTB breakdowns (new herd incidents) occurring between the years 2000 and 2014 were considered. In the first instance, cows with a positive reaction to the skin test and a positive postmortem examination were defined as infected. Values of 0 and 1 were assigned to healthy and infected animal records, respectively. Data were analyzed with mixed models. Linear and logit function heritability estimates were 0.092 and 0.172, respectively. In subsequent analyses, breakdowns were split into 2-mo intervals to better model time of exposure and infection in the contemporary group. Intervals with at least one infected individual were retained and multiple intervals within the same breakdown were included. Healthy animal records were assigned values of 0, and infected records a value of 1 in the interval of infection and values reflecting a diminishing probability of infection in the preceding intervals. Heritability and repeatability estimates were 0.115 and 0.699, respectively. Reliabilities and across time stability of the genetic evaluation were improved with the interval model. Subsequently, 2 more definitions of "infected" were analyzed with the interval model: (1) all positive skin test reactors regardless of postmortem examination, and (2) all positive skin test reactors plus nonreactors with positive postmortem examination. Estimated heritability was 0.085 and 0.089, respectively; corresponding repeatability estimates were 0.701 and 0.697. Genetic evaluation reliabilities and across time stability did not change. Correlations of genetic evaluations for bTB with other traits in the current breeding goal were mostly not different from zero. Correlation with the UK Profitable Lifetime Index was moderate, significant, and favorable. Results demonstrated the feasibility of a national genetic evaluation for bTB resistance. Selection for enhanced resistance will have a positive effect on profitability and no antagonistic effects on current breeding goal traits. Official genetic evaluations are now based on the interval model and the last bTB trait definition.

Genetic parameters of subclinical macromineral disorders and major clinical diseases in postparturient Holstein cows

The main objective of this study was to assess the genetic parameters of subclinical disorders associated with subclinical hypocalcemia, hypophosphatemia, subclinical hypomagnesemia, hypokalemia, and hyperphosphatemia, as well as major clinical diseases after calving in Holstein cows. The secondary objective was to estimate the associated genetic and phenotypic correlations among these subclinical and clinical conditions after calving in Holstein cows. The study was conducted in 9dairy herds located in Northern Greece. None of the herds used any kind of preventive measures for milk fever (MF). A total of 1,021 Holstein cows with pedigree information were examined from November 2010 until November 2012. The distribution across parities was 466 (parity 1), 242 (parity 2), 165 (parity 3), and 148 (parity 4 and above) cows. All cows were subjected to a detailed clinical examination and blood was sampled on d 1, 2, 4, and 8 after calving. Serum concentrations of Ca, P, Mg, and K were measured in all samples, whereas β-hydroxybutyrate (BHB) was measured only for d 8. The final data set included 4,064 clinical and 16,848 biochemical records (4,020 Ca, 4,019 P, 4,020Mg, 3,792K, and 997 BHB). Data of 1,988 observations of body condition score at d 1 and 8 were also available. All health traits were analyzed with a univariate random regression model. The genetic analysis for macromineral-related disorders included 986 cows with no obvious signs of MF (35 cows with MF were excluded). Analysis for other health traits included all 1,021 cows. A similar single record model was used for the analysis of BHB. Genetic correlations among traits were estimated with a series of bivariate analyses. Statistically significant daily heritabilities of subclinical hypocalcemia (0.13-0.25), hypophosphatemia (0.18-0.33), subclinical hypomagnesemia (0.11-0.38), and hyperphosphatemia (0.14-0.22) were low to moderate, whereas that of hypokalemia was low (0.08-0.10). The heritability of body condition score was 0.20±0.10. Statistically significant daily heritabilities of clinical diseases were those of MF (0.07-0.11), left displaced abomasum (0.19-0.31), and mastitis (0.15-0.41). Results suggest that these health disorders are heritable traits and could be minimized with proper genetic selection. Statistically significant phenotypic correlations were estimated for the first time between macromineral concentrations and almost all transition cow metabolic and infectious health disorders.

The genomic architecture of resistance to Campylobacter jejuni intestinal colonisation in chickens

BACKGROUND

Campylobacter is the leading cause of foodborne diarrhoeal illness in humans and is mostly acquired from consumption or handling of contaminated poultry meat. In the absence of effective licensed vaccines and inhibitors, selection for chickens with increased resistance to Campylobacter could potentially reduce its subsequent entry into the food chain. Campylobacter intestinal colonisation levels are influenced by the host genetics of the chicken. In the present study, two chicken populations were used to investigate the genetic architecture of avian resistance to colonisation: (i) a back-cross of two White Leghorn derived inbred lines [(61 x N) x N] known to differ in resistance to Campylobacter colonisation and (ii) a 9(th) generation advanced intercross (61 x N) line.

RESULTS

The level of colonisation with Campylobacter jejuni following experimental infection was found to be a quantitative trait. A back-cross experiment using 1,243 fully informative single nucleotide polymorphism (SNP) markers revealed quantitative trait loci (QTL) on chromosomes 7, 11 and 14. In the advanced intercross line study, the location of the QTL on chromosome 14 was confirmed and refined and two new QTLs were identified located on chromosomes 4 and 16. Pathway and re-sequencing data analysis of the genes located in the QTL candidate regions identified potential pathways, networks and candidate resistance genes. Finally, gene expression analyses were performed for some of the candidate resistance genes to support the results.

CONCLUSION

Campylobacter resistance in chickens is a complex trait, possibly involving the Major Histocompatibility Complex, innate and adaptive immune responses, cadherins and other factors. Two of the QTLs for Campylobacter resistance are co-located with Salmonella resistance loci, indicating that it may be possible to breed simultaneously for enhanced resistance to both zoonoses.

Comparison of inference methods of genetic parameters with an application to body weight in broilers

Abstract. REML (restricted maximum likelihood) has become the standard method of variance component estimation in animal breeding. Inference in Bayesian animal models is typically based upon Markov chain Monte Carlo (MCMC) methods, which are generally flexible but time-consuming. Recently, a new Bayesian computational method, integrated nested Laplace approximation (INLA), has been introduced for making fast non-sampling-based Bayesian inference for hierarchical latent Gaussian models. This paper is concerned with the comparison of estimates provided by three representative programs (ASReml, WinBUGS and the R package AnimalINLA) of the corresponding methods (REML, MCMC and INLA), with a view to their applicability for the typical animal breeder. Gaussian and binary as well as simulated data were used to assess the relative efficiency of the methods. Analysis of 2319 records of body weight at 35 days of age from a broiler line suggested a purely additive animal model, in which the heritability estimates ranged from 0.31 to 0.34 for the Gaussian trait and from 0.19 to 0.36 for the binary trait, depending on the estimation method. Although in need of further development, AnimalINLA seems a fast program for Bayesian modeling, particularly suitable for the inference of Gaussian traits, while WinBUGS appeared to successfully accommodate a complicated structure between the random effects. However, ASReml remains the best practical choice for the serious animal breeder.

Genetic parameters for predicted methane production and laser methane detector measurements

Enteric ruminant methane is the most important greenhouse gas emitted from the pastoral agricultural systems. Genetic improvement of livestock provides a cumulative and permanent impact on performance, and using high-density SNP panels can increase the speed of improvement for most traits. In this study, a data set of 1,726 dairy cows, collected since 1990, was used to calculate a predicted methane emission (PME) trait from feed and energy intake and requirements based on milk yield, live weight, feed intake, and condition score data. Repeated measurements from laser methane detector (LMD) data were also available from 57 cows. The estimated heritabilities for PME, milk yield, DMI, live weight, condition score, and LMD data were 0.13, 0.25, 0.11, 0.92, 0.38, and 0.05, respectively. There was a high genetic correlation between DMI and PME. No SNP reached the Bonferroni significance threshold for the PME traits. One SNP was within the 3 best SNP for PME at wk 10, 20, 30, and 40. Genomic prediction accuracies between dependent variable and molecular breeding value ranged between 0.26 and 0.30. These results are encouraging; however, more work is required before a PME trait can be implemented in a breeding program.

Short communication: Genetic characterization of digital cushion thickness

Dairy cow lameness is a serious animal welfare issue. It is also a significant cause of economic losses, reducing reproductive efficiency and milk production and increasing culling rates. The digital cushion is a complex structure composed mostly of adipose tissue located underneath the distal phalanx and has recently been phenotypically associated with incidence of claw horn disruption lesions (CHDL); namely, sole ulcers and white line disease. The objective of this study was to characterize digital cushion thickness genetically and to investigate its association with body condition score (BCS), locomotion score (LOCO), CHDL, and milk production. Data were collected from 1 large closely monitored commercial dairy farm located in upstate New York; 923 dairy cows were used. Before trimming, the following data were collected by a member of the research team: BCS, cow height measurement, and LOCO. Presence or not of CHDL (sole ulcer or white line disease, or both) was recorded at trimming. Immediately after the cows were hoof trimmed, they underwent digital sonographic B-mode examination for the measurement of digital cushion thickness. Factors such as parity number, stage of lactation, calving date, mature-equivalent 305-d milk yield (ME305MY), and pedigree information were obtained from the farm's dairy management software (DairyCOMP 305; Valley Agricultural Software, Tulare, CA). Univariate animal models were used to obtain variance component estimations for each studied trait (CHDL, BCS, digital cushion thickness average, LOCO, height, and ME305MY) and a 6-variate analysis was conducted to estimate the genetic, residual, and phenotypic correlations between the studied traits. The heritability estimate of DCTA was 0.33±0.09, whereas a statistically significant genetic correlation was estimated between DCTA and CHDL (-0.60±0.29). Of the other genetic correlations, significant estimates were derived for BCS with LOCO (-0.49±0.19) and ME305MY (-0.48±0.20). Digital cushion thickness is moderately heritable and genetically strongly correlated with CHDL.

Model comparison and estimation of genetic parameters for body weight in commercial broilers

Maniatis, G., Demiris, N., Kranis, A., Banos, G. and Kominakis, A. 2013. Model comparison and estimation of genetic parameters for body weight in commercial broilers. Can. J. Anim. Sci. 93: 67–77. The availability of powerful computing and advances in algorithmic efficiency allow for the consideration of increasingly complex models. Consequently, the development and application of appropriate statistical procedures for model evaluation is becoming increasingly important. This paper is concerned with the application of an alternative model determination criterion (conditional Akaike Information Criterion, cAIC) in a large dataset comprising 203 323 body weights of broilers, pertaining to 7 (BW7) and 35 (BW35) days of age. Seven univariate and seven bivariate models were applied. Direct genetic, maternal genetic and maternal environmental (c2) effects were estimated via REML. The model evaluation criteria included conditional Akaike Information Criterion (cAIC), Bayesian Information Criterion (BIC) and the standard Akaike Information Criterion (henceforth marginal; mAIC). According to cAIC the best-fitting model included direct genetic, maternal genetic and c2 effects. Maternal heritabilities were low (0.10 and 0.03) compared to the direct heritabilities (0.17 and 0.21), while c2 was 0.05 and 0.04 for BW7 and BW35, respectively. BIC and mAIC favoured a model that additionally included a direct-maternal genetic covariance, resulting in highly negative direct-maternal genetic correlations (−0.47 and −0.64 for BW7 and BW35, respectively) and higher direct heritabilities (0.25 and 0.28 for BW7 and BW35, respectively). Results suggest that cAIC can select different animal models than mAIC and BIC with different biological properties.

Description and typology of intensive Chios dairy sheep farms in Greece

The aim was to assess the intensified dairy sheep farming systems of the Chios breed in Greece, establishing a typology that may properly describe and characterize them. The study included the total of the 66 farms of the Chios sheep breeders' cooperative Macedonia. Data were collected using a structured direct questionnaire for in-depth interviews, including questions properly selected to obtain a general description of farm characteristics and overall management practices. A multivariate statistical analysis was used on the data to obtain the most appropriate typology. Initially, principal component analysis was used to produce uncorrelated variables (principal components), which would be used for the consecutive cluster analysis. The number of clusters was decided using hierarchical cluster analysis, whereas, the farms were allocated in 4 clusters using k-means cluster analysis. The identified clusters were described and afterward compared using one-way ANOVA or a chi-squared test. The main differences were evident on land availability and use, facility and equipment availability and type, expansion rates, and application of preventive flock health programs. In general, cluster 1 included newly established, intensive, well-equipped, specialized farms and cluster 2 included well-established farms with balanced sheep and feed/crop production. In cluster 3 were assigned small flock farms focusing more on arable crops than on sheep farming with a tendency to evolve toward cluster 2, whereas cluster 4 included farms representing a rather conservative form of Chios sheep breeding with low/intermediate inputs and choosing not to focus on feed/crop production. In the studied set of farms, 4 different farmer attitudes were evident: 1) farming disrupts sheep breeding; feed should be purchased and economies of scale will decrease costs (mainly cluster 1), 2) only exercise/pasture land is necessary; at least part of the feed (pasture) must be home-grown to decrease costs (clusters 1 and 4), 3) providing pasture to sheep is essential; on-farm feed production decreases costs (mainly cluster 3), and 4) large-scale farming (feed production and cash crops) does not disrupt sheep breeding; all feed must be produced on-farm to decrease costs (mainly cluster 3). Conducting a profitability analysis among different clusters, exploring and discovering the most beneficial levels of intensified management and capital investment should now be considered.

Validation of mid-infrared spectrometry in milk for predicting body energy status in Holstein-Friesian cows

Cow energy balance is known to be associated with cow health and fertility; therefore, routine access to data on energy balance can be useful in both management and breeding decisions to improve cow performance. The objective of this study was to determine if individual cow milk mid-infrared spectra (MIR) could be useful to predict cow energy balance across contrasting production systems. Direct energy balance was calculated as the differential between energy intake and energy output in milk and maintenance (maintenance was predicted using body weight). Body energy content was calculated from (change in) body weight and body condition score. Following editing, 2,992 morning, 2,742 midday, and 2,989 evening milk MIR records from 564 lactations on 337 Scottish cows, managed in a confinement system on 1 of 2 diets, were available. An additional 844 morning and 820 evening milk spectral records from 338 lactations on 244 Irish cows offered a predominantly grazed grass diet were also available. Equations were developed to predict body energy status using the milk spectral data and milk yield as predictor variables. Several different approaches were used to test the robustness of the equations calibrated in one data set and validated in another. The analyses clearly showed that the variation in the validation data set must be represented in the calibration data set. The accuracy (i.e., square root of the coefficient of multiple determinations) of predicting, from MIR, direct energy balance, body energy content, and energy intake was 0.47 to 0.69, 0.51 to 0.56, and 0.76 to 0.80, respectively. This highlights the ability of milk MIR to predict body energy balance, energy content, and energy intake with reasonable accuracy. Very high accuracy, however, was not expected, given the likely random errors in the calculation of these energy status traits using field data.

The use of mid-infrared spectrometry to predict body energy status of Holstein cows

Energy balance, especially in early lactation, is known to be associated with subsequent health and fertility in dairy cows. However, its inclusion in routine management decisions or breeding programs is hindered by the lack of quick, easy, and inexpensive measures of energy balance. The objective of this study was to evaluate the potential of mid-infrared (MIR) analysis of milk, routinely available from all milk samples taken as part of large-scale milk recording and milk payment operations, to predict body energy status and related traits in lactating dairy cows. The body energy status traits investigated included energy balance and body energy content. The related traits of body condition score and energy intake were also considered. Measurements on these traits along with milk MIR spectral data were available on 17 different test days from 268 cows (418 lactations) and were used to develop the prediction equations using partial least squares regression. Predictions were externally validated on different independent subsets of the data and the results averaged. The average accuracy of predicting body energy status from MIR spectral data was as high as 75% when energy balance was measured across lactation. These predictions of body energy status were considerably more accurate than predictions obtained from the sometimes proposed fat-to-protein ratio in milk. It is not known whether the prediction generated from MIR data are a better reflection of the true (unknown) energy status than the actual energy status measures used in this study. However, results indicate that the approach described may be a viable method of predicting individual cow energy status for a large scale of application.

Characterization of the PRNP gene locus in Chios dairy sheep and its association with milk production and reproduction traits

The objective of this study was to examine the prion protein gene locus (PRNP) in Chios sheep. PRNP is linked with scrapie resistance in small ruminants. Here, its impact on milk production (test-day and total lactation yield) and reproduction (age at first lambing, conception rate at first service, and prolificacy) was assessed. Genotyping at codons 136, 154 and 171 (classical scrapie) and 141 (atypical scrapie) was performed using DNA from milk somatic cells and PCR-RFLP analysis. A total of 1013 Chios ewes raised in 23 flocks were used. This constituted a random sample of the national breeding population. A total of 15 genotypes and 6 alleles linked to codons 136, 154 and 171 were detected. All animals were homozygous for the leucine allele at codon 141. Linear mixed models were used to assess the impact of PRNP genotypes and alleles on milk production and reproduction traits. The TRQ allele, whose association with such traits was assessed for the first time, had an adverse effect on age at first lambing. All other PRNP alleles, including ARR, which is associated with increased resistance to classical scrapie, had no significant effect on the traits studied. No significant associations of the PRNP genotypes with production and reproduction traits were observed. It was concluded that selection for scrapie-resistant sheep is not expected to affect the ongoing breeding programme that aims to enhance the milk yield and reproduction of the Chios breed.

Effect of polymorphisms at the STAT5A and FGF2 gene loci on reproduction, milk yield and lameness of Holstein cows

The objective of this study was to investigate the association of polymorphisms at the STAT5A and FGF2 gene loci with fertility, lactation milk yield and lameness in dairy cattle. Five hundred and eighteen primiparous Holstein cows were included in the study. Several reproductive traits were considered including conception rate (0/1) at first insemination, conception rate (0/1) in a 305-day lactation, number of inseminations per conception, interval (days) from calving to conception for cows that conceived in the first 305 days of lactation and age at first calving. Milk yield and lameness incidence were also recorded. Genotyping was performed using PCR-RFLP. The effect of allele substitution at each gene locus on reproductive traits, milk yield and lameness was assessed with single-trait mixed linear models. No significant associations were found between reproduction traits and any of the studied polymorphisms, apart from age at first calving, for which STAT5A polymorphism had a suggestive effect (P = 0.077). In addition, no significant effect of any polymorphism on lameness was found. Replacement of the C by G allele at the STAT5A locus was associated with a significant (P<0.05) increase in lactation milk yield, suggesting that this locus could be considered in gene assisted selection for the genetic improvement of milk production.

Impact of single nucleotide polymorphisms in leptin, leptin receptor, growth hormone receptor, and diacylglycerol acyltransferase (DGAT1) gene loci on milk production, feed, and body energy traits of UK dairy cows

The impact of 9 single nucleotide polymorphisms (SNP) in the leptin (LEP), leptin receptor (LEPR), growth hormone receptor (GHR), and diacylglycerol acyltransferase (DGAT1) gene loci on daily milk production, feed intake, and feed conversion, and weekly measures of live weight, BCS, and body energy traits was evaluated using genetic and phenotypic data on 571 Holstein cows raised at the Langhill Dairy Cattle Research Center in Scotland. Six SNP were typed on the LEP gene and 1 on each of the other 3 loci. Of the 6 LEP SNP, 3 were in very high linkage disequilibrium, meaning there is little gain in typing all of them in the future. Seven LEP haplotypes were identified by parsimony-based analyses. Random-regression allele-substitution models were used to assess the impact of each SNP allele or haplotype on the traits of interest. Diacylglycerol acyltransferase had a significant effect on milk yield, whereas GHR significantly affected feed intake, feed conversion, and body energy traits. There was also evidence of dominance in allelic effects on milk yield and BCS. The LEP haplotype CCGTTT (corresponding to leptin SNP C207T, C528T, A1457G, C963T, A252T, and C305T, respectively) significantly affected milk yield and feed and dry matter intake. Animals carrying this haplotype produced 3.13 kg more milk daily and consumed 4.64 kg more feed. Furthermore, they tended to preserve more energy than average. Such results may be used to facilitate genetic selection in animal breeding programs.