A prospective cohort study examining the association of claw anatomy and sole temperature with the development of claw horn disruption lesions in dairy cattle.

Foot characteristics have been linked to the development of sole lesions (sole hemorrhage and sole ulcers) and white line lesions, also known as claw horn disruption lesions (CHDL). The objective of this study was to examine the association of claw anatomy and sole temperature with the development of CHDL. A cohort of 2,352 cows was prospectively enrolled from 4 UK farms and assessed at 3 time points: before calving (T1-precalving), immediately after calving (T2-calving), and in early lactation. At each time point body condition score was recorded, a thermography image of each foot was taken for sole temperature measurement, the presence of CHDL was assessed by veterinary surgeons, and an ultrasound image was taken to retrospectively measure the digital cushion and sole horn thickness. Additionally, at the postcalving time point, foot angle and heel depth were recorded. Four multivariable logistic regression models were fit to separately examine the relationship of precalving and postcalving explanatory variables with the development of either white line lesions or sole lesions. Explanatory variables tested included digital cushion thickness, sole horn thickness, sole temperature, foot angle, and heel depth. Farm, parity, body condition score, and presence of lesion at the time of measurement were also included in the models. A thicker digital cushion shortly after calving was associated with decreased odds of cows developing sole lesions during early lactation (odds ratio [OR]: 0.74, 95% confidence interval [CI]: 0.65-0.84). No association was found between digital cushion thickness and development of white line lesions. Sole temperature after calving was associated with increased odds of the development of sole lesions (OR: 1.03, 95% CI: 1.02-1.05), and sole temperature before and after calving was associated with the development of white line lesions (T1-precalving OR: 1.04, 95% CI: 1.01-1.07; T2-calving OR: 0.96, 95% CI: 0.93-0.99). Neither foot angle nor heel depth was associated with the development of either lesion type. However, an increased sole horn thickness after calving reduced the odds of cows developing sole lesions during early lactation (OR: 0.88, 95% CI: 0.83-0.93), highlighting the importance of maintaining adequate sole horn when foot trimming. Before calving, animals with a lesion at the time of measurement and a thicker sole were more likely to develop a sole lesion (OR: 1.23, 95% CI: 1.09-1.40), compared with those without a sole lesion. The results presented here suggest that white line and sole lesions may have differing etiopathogenesis. Results also confirm the association between the thickness of the digital cushion and the development of sole lesions, highlight the association between sole horn thickness and sole lesions, and challenge the potential importance of foot angle and heel depth in the development of CHDL.

A prospective cohort study of factors associated with the digital cushion thickness in dairy cattle.

Key factors such as stage of lactation, parity and body fat reserves have been associated with the digital cushion thickness, however, there are discrepancies between the results of previously published studies. The objective of this study was to examine the association of stage of lactation, body fat reserves, parity, and lesion incidence with the digital cushion thickness (DCT) in a large cohort of intensively monitored cows. Across 4 UK farms, 2,352 cows were prospectively enrolled and assessed at 4 time points; before calving (T1-Precalving), immediately post-calving (T2-Calving), in early lactation (T3-Early) and late lactation (T4-Late). At each time point body condition score was recorded, the presence of sole lesions (sole ulcers and sole hemorrhage) and white line lesions were assessed by veterinarians, and an ultrasound image was taken to retrospectively measure the back-fat thickness in the pelvic (BFT) region and the digital cushion on the hind left lateral claw. Mixed effects multivariable linear regression models, with the cow as a random effect were fit to examine the association between explanatory variables and the DCT. Explanatory variables tested were farm, parity, stage of lactation, body condition score, BFT, height, the presence of a lesion at the time of measurement, the chronicity of a lesion during early lactation, predicted maximum daily milk yield and the rate of milk production rise in early lactation. Stage of lactation and farm were both associated with the DCT, however an interaction was present and this DCT pattern of change was farm dependent. Two distinct patterns emerged; one indicated the nadir to occur shortly after calving, the other indicated the nadir to occur during early lactation. Neither back fat thickness nor BCS were significantly associated with the DCT. Heifers displayed thinner digital cushions compared with multiparous cows, however, this effect was dependent on the stage of lactation, with heifers having a thinner digital cushion up until late lactation, by which time the DCT was commensurate with multiparous animals. Sole lesions and white line lesions at the time of measurement were associated with the DCT (sole lesion; Estimate: -0.07mm, 95% CI: -0.14-0.00, P = 0.039, white line lesion; Estimate: 0.28mm, 95% CI: 0.15-0.42, P < 0.001).

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach.

BACKGROUND: The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds. RESULTS: The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne). CONCLUSIONS: Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

Single-step genome-wide association analyses of claw horn lesions in Holstein cattle using linear and threshold models.

BACKGROUND: Lameness in dairy cattle is primarily caused by foot lesions including the claw horn lesions (CHL) of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). This study investigated the genetic architecture of the three CHL based on detailed animal phenotypes of CHL susceptibility and severity. Estimation of genetic parameters and breeding values, single-step genome-wide association analyses, and functional enrichment analyses were performed. RESULTS: The studied traits were under genetic control with a low to moderate heritability. Heritability estimates of SH and SU susceptibility on the liability scale were 0.29 and 0.35, respectively. Heritability of SH and SU severity were 0.12 and 0.07, respectively. Heritability of WL was relatively lower, indicating stronger environmental influence on the presence and development of WL than the other two CHL. Genetic correlations between SH and SU were high (0.98 for lesion susceptibility and 0.59 for lesion severity), whereas genetic correlations of SH and SU with WL also tended to be positive. Candidate quantitative trait loci (QTL) were identified for all CHL, including some on Bos taurus chromosome (BTA) 3 and 18 with potential pleiotropic effects associated with multiple foot lesion traits. A genomic window of 0.65 Mb on BTA3 explained 0.41, 0.50, 0.38, and 0.49% of the genetic variance for SH susceptibility, SH severity, WL susceptibility, and WL severity, respectively. Another window on BTA18 explained 0.66, 0.41, and 0.70% of the genetic variance for SH susceptibility, SU susceptibility, and SU severity, respectively. The candidate genomic regions associated with CHL harbour annotated genes that are linked to immune system function and inflammation responses, lipid metabolism, calcium ion activities, and neuronal excitability. CONCLUSIONS: The studied CHL are complex traits with a polygenic mode of inheritance. Most traits exhibited genetic variation suggesting that animal resistance to CHL can be improved with breeding. The CHL traits were positively correlated, which will facilitate genetic improvement for resistance to CHL as a whole. Candidate genomic regions associated with lesion susceptibility and severity of SH, SU, and WL provide insights into a global profile of the genetic background underlying CHL and inform genetic improvement programmes aiming at enhancing foot health in dairy cattle.

Genetic parameters of sole lesion recovery in Holstein cows.

Sole hemorrhage and sole ulcers, referred to as sole lesions, are important causes of lameness in dairy cattle. The objective of this study was to estimate the genetic parameters of a novel trait reflecting how well cows recovered from sole lesions and the genetic correlation of this trait with overall susceptibility to sole lesions. A cohort of Holstein dairy cows was prospectively enrolled on 4 farms and assessed at 4 timepoints: before calving, immediately after calving, in early lactation, and in late lactation. At each timepoint, sole lesions were recorded at the claw level by veterinary surgeons and used to define 2 binary traits: (1) susceptibility to sole lesions-whether animals were affected with sole lesions at least once during the study or were unaffected at every assessment, and (2) sole lesion recovery-whether sole lesions healed between early and late lactation. Animals were genotyped and pedigree details extracted from the national database. Analyses were conducted with BLUPF90 software in a single-step framework; genetic parameters were estimated from animal threshold models using Gibbs sampling. The genetic correlation between both traits was approximated as the correlation between genomic estimated breeding values, adjusting for their reliabilities. A total of 2,025 animals were used to estimate the genetic parameters of sole lesion susceptibility; 44% of animals recorded a sole lesion at least once during the study period. The heritability of sole lesion susceptibility, on the liability scale, was 0.25 (95% highest density interval = 0.16-0.34). A total of 498 animals were used to estimate the genetic parameters of sole lesion recovery; 71% of animals had recovered between the early and late lactation assessments. The heritability of sole lesion recovery, on the liability scale, was 0.27 (95% highest density interval = 0.02-0.52). The approximate genetic correlation between each trait was -0.11 (95% confidence interval = -0.20 to -0.02). Our results indicate that recovery from sole lesions is heritable. If this finding is corroborated in further studies, it may be possible to use selective breeding to reduce the frequency of chronically lame cows. As sole lesion recovery appears to be weakly genetically related to sole lesion susceptibility, successful genetic improvement of sole lesion recovery would benefit from selection on this trait directly.

Serum 1H nuclear magnetic resonance-based metabolomics of sole lesion development in Holstein cows.

Sole hemorrhage and sole ulcers, referred to as sole lesions, are important causes of lameness in dairy cattle. We aimed to compare the serum metabolome of dairy cows that developed sole lesions in early lactation with that of cows that remained unaffected. We prospectively enrolled a cohort of 1,169 Holstein dairy cows from a single dairy herd and assessed animals at 4 time points: before calving, immediately after calving, early lactation, and late lactation. Sole lesions were recorded by veterinary surgeons at each time point, and serum samples were collected at the first 3 time points. Cases were defined by the presence of sole lesions in early lactation and further subdivided by whether sole lesions had been previously recorded; unaffected controls were randomly selected to match cases. Serum samples from a case-control subset of 228 animals were analyzed with proton nuclear magnetic resonance spectroscopy. Spectral signals, corresponding to 34 provisionally annotated metabolites and 51 unlabeled metabolites, were analyzed in subsets relating to time point, parity cohort, and sole lesion outcome. We used 3 analytic methods (partial least squares discriminant analysis, least absolute shrinkage and selection operator regression, and random forest) to determine the predictive capacity of the serum metabolome and identify informative metabolites. We applied bootstrapped selection stability, triangulation, and permutation to support the inference of variable selection. The average balanced accuracy of class prediction ranged from 50 to 62% depending on the subset. Across all 17 subsets, 20 variables had a high probability of being informative; those with the strongest evidence of being associated with sole lesions corresponded to phenylalanine and 4 unlabeled metabolites. We conclude that the serum metabolome, as characterized by proton nuclear magnetic resonance spectroscopy, does not appear able to predict sole lesion presence or future development of lesions. A small number of metabolites may be associated with sole lesions although, given the poor prediction accuracies, these metabolites are likely to explain only a small proportion of the differences between affected and unaffected animals. Future metabolomic studies may reveal underlying metabolic mechanisms of sole lesion etiopathogenesis in dairy cows; however, the experimental design and analysis need to effectively control for interanimal and extraneous sources of spectral variation.

Including genotypic information in genetic evaluations increases the accuracy of sheep breeding values.

The impact of inclusion of genome-wide genotypes into breeding value predictions for UK Texel sheep is addressed in this article. The main aim was to investigate the level of change in the accuracy values for EBVs when information from animal genotypes is incorporated into the genetic evaluations. New genetic parameters for a range of lamb growth, carcass composition and health traits are described and applied in the estimation of conventional breeding values (EBVs) for almost 822,000 animals as well as genomic breeding values (gEBVs) after adding 10,143 genotypes. Principal component analyses showed that there are no major distinct groups; hence, the population is mainly homogenous and genetically well-linked. Results suggested that the highest change in accuracy was observed for the animals that are not phenotyped but have good links to the reference population. This was seen especially for the lowly heritable health traits thereby proving that the use of genotypes in breeding values estimation may accelerate the genetic gain by producing more accurate values especially for young, un-phenotyped animals.

Genetic parameters and genome-wide association study of digital cushion thickness in Holstein cows.

The digital cushion is linked to the development of claw horn lesions (CHL) in dairy cattle. The objectives of this study were to (1) estimate genetic parameters for digital cushion thickness (DCT), (2) estimate the genetic correlation between DCT and CHL, and (3) identify candidate genes associated with DCT. A cohort of 2,352 Holstein dairy cows were prospectively enrolled on 4 farms and assessed at 4 time points: before calving, immediately after calving, in early lactation, and in late lactation. At each time point, CHL was recorded by veterinary surgeons, and ultrasonographic images of the digital cushion were stored and retrospectively measured at 2 anatomical locations. Animals were genotyped and pedigree details extracted from the national database. Genetic parameters were estimated following a single-step approach implemented in AIREMLF90. Four traits were analyzed: the 2 DCT measurements, sole lesions (sole hemorrhage and sole ulcers), and white line lesions. All traits were analyzed with univariate linear mixed models; bivariate models were fit to estimate the genetic correlation between traits within and between time points. Single-marker and window-based genome-wide association analyses of DCT traits were conducted at each time point; candidate genes were mapped near (<0.2 Mb) or within the genomic markers or windows with the largest effects. Heritability estimates of DCT ranged from 0.14 to 0.44 depending on the location of DCT measurement and assessment time point. The genetic correlation between DCT and sole lesions was generally negative, notably between DCT immediately after calving and sole lesions in early or late lactation, and between DCT in early or late lactation and sole lesion severity in early or late lactation. Digital cushion thickness was not genetically correlated with white line lesions. A polygenic background to DCT was found; genes associated with inflammation, fat metabolism, and bone development were mapped near or within the top markers and windows. The moderate heritability of DCT provides an opportunity to use selective breeding to change DCT in a population. The negative genetic correlation between DCT and sole lesions at different stages of production lends support to current hypotheses of sole lesion pathogenesis. Highlighted candidate genes provide information regarding the complex genetic background of DCT in Holstein cows, but further studies are needed to explore and corroborate these findings.

Across-country genetic evaluation of meat sheep from Ireland and the United Kingdom.

Genetic evaluations in sheep have proven to be an effective way of increasing farm profitability. Much research has previously been conducted on producing within-country genetic evaluations; however, to date, no across-country sheep genetic evaluations have been produced between Ireland and the UK. The objective of the present study was to examine the feasibility of an across-country genetic evaluation of live body weight and carcass composition traits for Texel sheep raised in Ireland and the UK. The benefit of genetic selection based on across-country genetic evaluations, in comparison with within-country genetic evaluations, was also quantified. Animal traits included early-life and postweaning live body weights, and muscle and fat depth ultrasound measurements. Irish and UK data were combined, common animals with progeny with records in both countries were identified and a series of bivariate analyses were performed separately for each trait to produce across-country genetic evaluations. Fixed effects included contemporary group, age at first lambing of the dam, parity of the dam (Ireland), dam age at lamb's birth (UK), a gender by age of the lamb interaction, a birth type by rearing type of the lamb interaction and country of birth of the lamb. Random effects included the animal additive genetic, dam maternal, litter common environment and residual effect. The model for postweaning weight, muscle depth and fat depth included only the animal additive genetic and litter common environmental random effects. Genetic correlations between the two countries ranged from 0.82 to 0.88 for the various traits. Across-country breeding values were estimated for all animals and response to selection was predicted using the top 10 and top 20 sires in both within- and across-country analyses for the two countries. Overall, results showed that rates of genetic gain could potentially increase from between 2.59% and 19.63% from selection based on across-country genetic evaluations compared to within-country evaluations alone. Across-country evaluations are feasible and would be of significant benefit to both the Irish and UK sheep industries. In order to realize these potential gains though, there would need to be a switch in emphasis by sheep breeders towards using objective traits as their primary selection criteria.

A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep.

BACKGROUND: The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change. RESULTS: We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h2 = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response. CONCLUSIONS: Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest.

Breeding strategies for improving smallholder dairy cattle productivity in Sub-Saharan Africa.

Breeding strategies for smallholder dairy farming systems in Sub-Saharan Africa (SSA) were simulated and evaluated considering cow traits identified as priorities by farmers in different agro-ecological zones. These traits were related to cow milk yield, fertility, temperament, feed intake and disease resistance. The first breeding strategy was based on continuous importation of genetically superior exotic dairy sire semen to SSA and crossing with local females leading to a gradual upgrade of the indigenous population. The second strategy assumed that semen from elite exotic bulls would be imported to SSA and used on indigenous cows to produce F1 animals. Thereafter, elite animals would be selected from within the F1 and each subsequent generation to establish a new synthetic breed. The third strategy was to improve the indigenous population by genetically selecting the best sires available domestically. Results showed positive genetic progress for all breeding goal traits. After 15 generations of selection, the genetic response of the importation strategy exceeded the corresponding genetic response of the synthetic breed strategy by 20%-60%. The former also exceeded the genetic response of the indigenous breed improvement strategy by 43%-75%. Potentially there is an opportunity for breeders to choose an appropriate breeding strategy that fits a specific need of smallholder dairy farmers.

Breeding strategies for animal resilience to weather variation in meat sheep.

BACKGROUND: The alteration in weather patterns expected due to climate change will affect farm animal performance, probably resulting in lower quantity and quality of available products. A potential mitigation strategy would be to breed selected animals for enhanced resilience to climate change. In this context, resilience would reflect stable animal performance in spite of weather variation. The objectives of this study were to (i) derive and characterise novel animal resilience phenotypes, (ii) investigate their genetic profiles and (iii) assess the impact of integrating them in breeding strategies for genetic improvement in meat sheep. RESULTS: Random regression models were used to jointly analyse live body weight measured in different time points throughout the growth phases of 4469 Scottish Blackface sheep and weather variables during the same period to derive novel resilience phenotypes. The genetic analysis of these phenotypes revealed significant genetic variance and heritability, and an antagonistic genetic correlation with some animal performance traits. Simulated breeding strategies demonstrated that a relative emphasis of 10% on resilience compared to other traits would enhance performance stability against weather volatility without compromising animal growth. CONCLUSIONS: Novel resilience traits exhibited sufficient genetic variation to be amenable to genetic improvement with selective breeding and are recommended to be included in future breeding goals.

ACAA2 and FASN polymorphisms affect the fatty acid profile of Chios sheep milk.

The objective of the research reported in this research communication was the identification and association of single nucleotide polymorphisms (SNP) in the ovine DGAT1, FASN, SCD1 and ACAA2 genes with milk fat percentage and fatty acid (FA) content. Three consecutive monthly milk samplings were obtained from a total of 429 purebred Chios ewes during mid-lactation. Genotypic data were jointly analyzed with 1184 fat content and 37 718 FA percentage records using mixed models. The 3' untranslated region (UTR) of the DGAT1 gene and the 5' and 3'UTRs of the SCD1 gene appeared to be monomorphic. The FASN g.14777C>T SNP on exon 31 was associated with C13:0 and the ACAA2 g.2982T>C SNP on the 3'UTR was associated with C9:0, C11:0, C12:1 cis-9, C13:0 and the ω6/ω3 index, while fat percentage was not affected by the identified SNPs. The results could be useful for breeding programs aiming to improve the quality and nutritional value of ovine milk.

Genetic analysis of novel phenotypes for farm animal resilience to weather variability.

BACKGROUND: Climate change is expected to have a negative impact on food availability. While most efforts have been directed to reducing greenhouse gas emissions, complementary strategies are necessary to control the detrimental effects of climate change on farm animal performance. The objective of this study was to develop novel animal resilience phenotypes using reaction norm slopes, and examine their genetic and genomic parameters. A closely monitored dairy goat population was used for this purpose. RESULTS: Individual animals differed in their response to changing atmospheric temperature and a temperature-humidity index. Significant genetic variance and heritability estimates were derived for these animal resilience phenotypes. Furthermore, some resilience traits had a significant unfavourable genetic correlation with animal performance. Genome-wide association analyses identified several candidate genes related to animal resilience to environment change. CONCLUSIONS: Heritable variation exists among dairy goats in their production response to fluctuating weather variables. Results may inform future breeding programmes aimed to ensure efficient animal performance under changing climatic conditions.

Immune-associated traits measured in milk of Holstein-Friesian cows as proxies for blood serum measurements.

Previous work has highlighted that immune-associated (IA) traits measurable in blood are associated with health, productivity, and reproduction in dairy cows. The aim of the present study was to determine relationships between IA traits measured in blood serum and those simultaneously measured in milk as well as their association with disease phenotypes. All animals were Holstein-Friesian cows from the Langhill research herd (n = 546) housed at the SRUC Dairy Research Centre in Scotland. Milk and serum samples were collected on 20 separate occasions between July 2010 and March 2015 and analyzed by ELISA for haptoglobin (Hp), tumor necrosis factor-α (TNF-α), and natural antibodies binding keyhole limpet hemocyanin (NAbKLH) and lipopolysaccharide (NAbLPS). Data were analyzed using mixed linear models that included pedigree information. Analyses revealed positive phenotypic correlations between milk and serum NAb (0.59 ≤ r ≤ 0.77), Hp (r = 0.37), and TNF-α (r = 0.12). Milk and serum NAb were also found to have a strong genetic correlation (0.81 ≤ r ≤ 0.94) and were genetically correlated with cow lameness (0.66 and 0.79 for milk NAbKLH and serum NAbLPS, respectively). Clinical mastitis was found to be phenotypically correlated with both milk and serum Hp (0.09 ≤ r ≤ 0.23). Serum Hp was also strongly genetically correlated with other cellular IA traits such as percent NKp46+ (a natural killer cell marker; 0.35) and percent peripheral blood mononuclear cells (PBMC; -0.90). Similarly, genetic correlations were found to exist between serum TNF-α and percent NKp46+ (0.22), percent PBMC (0.41), and percent lymphocytes (0.47). Excluding serum Hp, all milk and serum IA traits were repeatable, ranging from 0.11 (milk Hp) to 0.43 (serum NAbLPS). Between-animal variation was highest in milk and serum NAb (0.34-0.43) and significant estimates of heritability were also observed in milk and serum NAb (0.17-0.37). Our findings show that certain IA traits, such as NAbKLH and NAbLPS, found in milk and serum are strongly correlated and highlight the potential of using routinely collected milk samples as a less invasive and cost-effective source of informative data for predictive modeling of animal IA traits.

Genomic regions underlying susceptibility to bovine tuberculosis in Holstein-Friesian cattle.

BACKGROUND: The significant social and economic loss as a result of bovine tuberculosis (bTB) presents a continuous challenge to cattle industries in the UK and worldwide. However, host genetic variation in cattle susceptibility to bTB provides an opportunity to select for resistant animals and further understand the genetic mechanisms underlying disease dynamics. METHODS: The present study identified genomic regions associated with susceptibility to bTB using genome-wide association (GWA), regional heritability mapping (RHM) and chromosome association approaches. Phenotypes comprised de-regressed estimated breeding values of 804 Holstein-Friesian sires and pertained to three bTB indicator traits: i) positive reactors to the skin test with positive post-mortem examination results (phenotype 1); ii) positive reactors to the skin test regardless of post-mortem examination results (phenotype 2) and iii) as in (ii) plus non-reactors and inconclusive reactors to the skin tests with positive post-mortem examination results (phenotype 3). Genotypes based on the 50 K SNP DNA array were available and a total of 34,874 SNPs remained per animal after quality control. RESULTS: The estimated polygenic heritability for susceptibility to bTB was 0.26, 0.37 and 0.34 for phenotypes 1, 2 and 3, respectively. GWA analysis identified a putative SNP on Bos taurus autosomes (BTA) 2 associated with phenotype 1, and another on BTA 23 associated with phenotype 2. Genomic regions encompassing these SNPs were found to harbour potentially relevant annotated genes. RHM confirmed the effect of these genomic regions and identified new regions on BTA 18 for phenotype 1 and BTA 3 for phenotypes 2 and 3. Heritabilities of the genomic regions ranged between 0.05 and 0.08 across the three phenotypes. Chromosome association analysis indicated a major role of BTA 23 on susceptibility to bTB. CONCLUSION: Genomic regions and candidate genes identified in the present study provide an opportunity to further understand pathways critical to cattle susceptibility to bTB and enhance genetic improvement programmes aiming at controlling and eradicating the disease.

Estimating genetic and phenotypic parameters of cellular immune-associated traits in dairy cows.

Data collected from an experimental Holstein-Friesian research herd were used to determine genetic and phenotypic parameters of innate and adaptive cellular immune-associated traits. Relationships between immune-associated traits and production, health, and fertility traits were also investigated. Repeated blood leukocyte records were analyzed in 546 cows for 9 cellular immune-associated traits, including percent T cell subsets, B cells, NK cells, and granulocytes. Variance components were estimated by univariate analysis. Heritability estimates were obtained for all 9 traits, the highest of which were observed in the T cell subsets percent CD4+, percent CD8+, CD4+:CD8+ ratio, and percent NKp46+ cells (0.46, 0.41, 0.43 and 0.42, respectively), with between-individual variation accounting for 59 to 81% of total phenotypic variance. Associations between immune-associated traits and production, health, and fertility traits were investigated with bivariate analyses. Strong genetic correlations were observed between percent NKp46+ and stillbirth rate (0.61), and lameness episodes and percent CD8+ (-0.51). Regarding production traits, the strongest relationships were between CD4+:CD8+ ratio and weight phenotypes (-0.52 for live weight; -0.51 for empty body weight). Associations between feed conversion traits and immune-associated traits were also observed. Our results provide evidence that cellular immune-associated traits are heritable and repeatable, and the noticeable variation between animals would permit selection for altered trait values, particularly in the case of the T cell subsets. The associations we observed between immune-associated, health, fertility, and production traits suggest that genetic selection for cellular immune-associated traits could provide a useful tool in improving animal health, fitness, and fertility.

A practical approach to detect ancestral haplotypes in livestock populations.

BACKGROUND: The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle. RESULTS: The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content. CONCLUSION: In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.

Genome-wide association studies of immune, disease and production traits in indigenous chicken ecotypes.

BACKGROUND: The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek's disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses. RESULTS: The large majority of birds was found to have antibody titres for all pathogens and were infected with both parasites, suggesting almost universal exposure. We derived significant moderate to high heritabilities for IBDV, MDV and PM antibody titres, cestodes infestation, body weight and BCS. We identified single nucleotide polymorphisms (SNPs) with genome-wide significance for each trait. Based on these associations, we identified for each trait, pathways, networks and functional gene clusters that include plausible candidate genes. Selective sweep analyses revealed a locus on chromosome 18 associated with viral antibody titres and resistance to Eimeria parasitism that is within a positive selection signal. We found no significant genetic correlations between production, immune and disease traits, implying that selection for altered antibody response and/or disease resistance will not affect production. CONCLUSIONS: We confirmed the presence of genetic variability and identified SNPs significantly associated with immune, disease and production traits in indigenous village chickens. Results underpin the feasibility of concomitant genetic improvement for enhanced antibody response, resistance to parasitism and productivity within and across indigenous chicken ecotypes.

Association of lameness with milk yield and lactation curves in Chios dairy ewes.

The objective of the study was twofold: (i) to quantify the differences in daily milk yield (DMY) and total milk yield (TMY) between lame and non-lame dairy ewes and (ii) to determine the shape of lactation curves around the lameness incident. The overall study was a prospective study of lameness for the surveyed sheep population, with a nested study including the selection of matching controls for each lame ewe separately. Two intensively reared flocks of purebred Chios ewes and a total of 283 ewes were used. Data, including gait assessment and DMY records, were collected on a weekly basis during on-farm visits across the milking period. A general linear model was developed for the calculation of lactation curves of lame and non-lame ewes, whereas one-way ANOVA was used for the comparisons between lame ewes and their controls. Lameness incidence was 12.4 and 16.8% on Farms A and B, respectively. Average DMY in lame ewes was significantly lower (213.8 g, P < 0.001) compared with the rest of the flock, where DMY averaged 1.340 g. The highest DMY reduction in lame ewes was observed during the week 16 of the milking period (P < 0.001), whereas the reduction of DMY, for lame ewes, remained significant at P < 0.001 level from week 8 to week 28 of milking. Comparisons between lame and controls revealed that at the week of lameness diagnosis a significant DMY reduction (P ≤ 0.001) was observed in lame ewes (about 32.5%), which was maximised 1 week later (35.8%, P ≤ 0.001) and continued for several weeks after recovery, resulting in 19.3% lower TMY for lame ewes for the first 210 d of the milking period (P < 0.01). Moreover, at flock level, TMY for non-lame and lame ewes, as calculated by the general linear model, was 318.9 and 268.0 kg, respectively. The results of this study demonstrate evidence of significant financial losses in dairy sheep due to lameness which, however, need to be accurately estimated in further, more detailed, analyses.