Genetics of animal health and disease in cattle

Optimizing Cost-Effective gene expression phenotyping approaches in cattle using 3' mRNA sequencing

BACKGROUND

Genetic and genomic selection programs require large numbers of phenotypes observed for animals in shared environments. Direct measurements of phenotypes like meat quality, methane emission, and disease susceptibility are difficult and expensive to measure at scale but are critically important to livestock production. Our work leans on our understanding of the "Central Dogma" of molecular genetics to leverage molecular intermediates as cheaply-measured proxies of organism-level phenotypes. The rapidly declining cost of next-generation sequencing presents opportunities for population-level molecular phenotyping. While the cost of whole transcriptome sequencing has declined recently, its required sequencing depth still makes it an expensive choice for wide-scale molecular phenotyping. We aim to optimize 3' mRNA sequencing (3' mRNA-Seq) approaches for collecting cost-effective proxy molecular phenotypes for cattle from easy-to-collect tissue samples (i.e., whole blood). We used matched 3' mRNA-Seq samples for 15 Holstein male calves in a heat stress trail to identify the (1) best library preparation kit (Takara SMART-Seq v4 3' DE and Lexogen QuantSeq) and (2) optimal sequencing depth (0.5 to 20 million reads/sample) to capture gene expression phenotypes most cost-effectively.

RESULTS

Takara SMART-Seq v4 3' DE outperformed Lexogen QuantSeq libraries across all metrics: number of quality reads, expressed genes, informative genes, differentially expressed genes, and 3' biased intragenic variants. Serial downsampling analyses identified that as few as 8.0 million reads per sample could effectively capture most of the between-sample variation in gene expression. However, progressively more reads did provide marginal increases in recall across metrics. These 3' mRNA-Seq reads can also capture animal genotypes that could be used as the basis for downstream imputation. The 10 million read downsampled groups called an average of 109,700 SNPs and 11,367 INDELs, many of which segregate at moderate minor allele frequencies in the population.

CONCLUSION

This work demonstrates that 3' mRNA-Seq with Takara SMART-Seq v4 3' DE can provide an incredibly cost-effective (< 25 USD/sample) approach to quantifying molecular phenotypes (gene expression) while discovering sufficient variation for use in genotype imputation. Ongoing work is evaluating the accuracy of imputation and the ability of much larger datasets to predict individual animal phenotypes.

Single nucleotide polymorphisms, gene expression and evaluation of immunological, antioxidant, and pathological parameters associated with bacterial pneumonia in Barki sheep

BACKGROUND

In sheep, pneumonia is a major concern because of its high morbidity, mortality, and economic impact. It results from various infectious agents, including bacteria, viruses, and environmental stressors, that weaken the immune system.

OBJECTIVE

The objective of this study was to monitor nucleotide sequence variations, gene expression, and serum biomarkers of inflammation and oxidative stress in sheep with pneumonia. Additionally, this study aimed to identify various bacterial strains and virulent gene combinations in pneumonic sheep, as confirmed by PCR.

METHODOLOGY

The enrolled animals were categorized as follows: 50 apparently healthy ewes, considered the control group, and 150 infected ewes with pneumonia. The infected ewes included 100 sporadic cases from the Center for Sustainable Development of Matrouh Resources, Desert Research Center, Matrouh, Egypt, and 50 ewes from the slaughterhouse, all exhibiting respiratory symptoms such as coughing, serous to mucopurulent nasal discharge, fever, and abnormal lung sounds. Blood samples were collected to assess various biochemical parameters, detect SNPs, and analyse the expression of specific immunological and antioxidant-related genes. Nasopharyngeal and lung swabs were taken from the affected ewes for bacteriological analysis, and lung samples were collected for histological examination.

RESULTS

Phenotypic characterization and identification revealed the presence of Klebsiella pneumoniae, Pasteurella multocida, Mannheimia haemolytica, Pseudomonas spp., Mycoplasma, Streptococcus, and Escherichia coli, with frequencies of 40%, 28.6%, 34%, 18%, 44%, 29.3%, and 20%, respectively. Additionally, virulence genes for Klebsiella pneumoniae, iutA and fimH, were detected at rates of 39% and 68%, respectively, whereas the toxA gene for Pseudomonas spp. was present in 59.2% of the cases. Nucleotide sequence variations in immunity- and antioxidant-related genes were observed between healthy and pneumonic ewes. The genes encoding IL-1α, IL1B, IL6, TNF-α, LFA-1, CR2, IL17, IL13, DEFB123, SCART1, ICAM1, NOS, and HMOX1 were significantly upregulated in pneumonia-affected ewes compared with resistant ewes. Conversely, the genes encoding IL10, SOD1, CAT, GPX1, and NQO1 were downregulated. Further analysis of the serum profile revealed a significant (P < 0.05) increase in IL-1α, IL-1β, IL-6, TNF-α, NO and MDA along with a significant (P < 0.05) decrease in the serum levels of C3, C4, CAT, GPx, GR and IL-10 in diseased ewes compared with healthy ewes. Histopathological examination revealed that the infected sheep exhibited broncho-interstitial pneumonia and purulent to fibrino-purulent bronchopneumonia.

CONCLUSIONS

This study revealed the significant presence of various pathogens and virulence factors in infected sheep, along with distinct immunological and antioxidant gene expression patterns. The altered serum profile and gene regulation in pneumonia-affected ewes underscore the complex immune response and potential biomarkers for disease susceptibility and resistance.

Insights into garlic (Allium Sativum)'s nutrigenomics-associated fly-repellent potency in cattle

Despite effective control of flies using synthetic pesticides, fly resistance and environmental contamination have led to the inadequacy of this strategy. The use of integrated pest management approaches has since been advocated in contemporary research to sustainably control fly populations. Recent studies have found garlic (Allium Sativum) and its derivative bioactive compounds to possess insect-repellent attributes among other key health and production enhancing properties. This highlights the potential of garlic as a botanical pesticide to control flies in cattle. Moreover, the ability of cattle to naturally repel flies is influenced by animal genetic predisposition. The dietary garlic supplementation and gene interaction in disease resistance could also be an influential factor in repelling flies in cattle. Transcriptomics has emerged as a valuable tool in animal breeding and genetics which allows identification of trait-associated genes and understanding of complex interactions between dietary nutrients and animal genome expression. This paper explores the nutrigenomic effects of garlic supplementation on cattle and its contribution towards fly repellence efficacy in cattle. It was concluded that garlic supplementation in cattle diets could offer a sustainable approach to managing fly infestations in cattle farming. These findings underscore the importance of further research to validate these assertions and optimise the use of garlic to control flies in cattle under different production systems.

Genetic advancements and future directions in ruminant livestock breeding: from reference genomes to multiomics innovations

Ruminant livestock provide a rich source of products, such as meat, milk, and wool, and play a critical role in global food security and nutrition. Over the past few decades, genomic studies of ruminant livestock have provided valuable insights into their domestication and the genetic basis of economically important traits, facilitating the breeding of elite varieties. In this review, we summarize the main advancements for domestic ruminants in reference genome assemblies, population genomics, and the identification of functional genes or variants for phenotypic traits. These traits include meat and carcass quality, reproduction, milk production, feed efficiency, wool and cashmere yield, horn development, tail type, coat color, environmental adaptation, and disease resistance. Functional genomic research is entering a new era with the advancements of graphical pangenomics and telomere-to-telomere (T2T) gap-free genome assembly. These advancements promise to improve our understanding of domestication and the molecular mechanisms underlying economically important traits in ruminant livestock. Finally, we provide new perspectives and future directions for genomic research on ruminant genomes. We suggest how ever-increasing multiomics datasets will facilitate future studies and molecular breeding in livestock, including the potential to uncover novel genetic mechanisms underlying phenotypic traits, to enable more accurate genomic prediction models, and to accelerate genetic improvement programs.

Analysis of Potential Genes, Oxidative, Metabolic, and Hormonal Markers Associated with Postpartum Disorder Susceptibility in Barki Sheep (Ovis aries)

This study purpose was to determine the gene expression as well as serum profile of acute phase proteins (APPs) and hormonal indicators linked to Barki sheep's susceptibility to postpartum issues. Three equal-sized groups (each with fifty ewes) were created from the blood of 150 adult Barki ewes: the control group (CG), the inflammatory postpartum disorders group (IPG), and the non-inflammatory postpartum disorders group (NIPG). The expression levels of the oxidative stress (PGC-1α, SIRT1, GCLC, GCLM, and EPAS1) and metabolic (FBXL12, KPNA7, and LRRK1) genes were significantly higher in postpartum disorders sheep than in resistant ones. Ewes with inflammatory postpartum illnesses showed significantly higher levels of the examined markers than did the non-inflammatory and control groups. The serum profile analysis also revealed that the levels of Fb, Cp, Hp, SAA, cortisol, TIBC, UIBC, and ferritin were significantly higher in the IPG than in the NIPG and CG. Serum insulin, iron, transferrin, and Tf Sat.% levels, however, were all markedly lower. On the basis of the variance in the genes being studied and the modulation in the serum indicators being studied, it should be possible to monitor the health status in postpartum problems of sheep.

Mendelian randomisation to uncover causal associations between conformation, metabolism, and production as potential exposure to reproduction in German Holstein dairy cattle

BACKGROUND

Reproduction is vital to welfare, health, and economics in animal husbandry and breeding. Health and reproduction are increasingly being considered because of the observed genetic correlations between reproduction, health, conformation, and performance traits in dairy cattle. Understanding the detailed genetic architecture underlying these traits would represent a major step in comprehending their interplay. Identifying known, putative or novel associations in genomics could improve animal health, welfare, and performance while allowing further adjustments in animal breeding.

RESULTS

We conducted genome-wide association studies for 25 different traits belonging to four different complexes, namely reproduction (n = 13), conformation (n = 6), production (n = 3), and metabolism (n = 3), using a cohort of over 235,000 dairy cows. As a result, we identified genome-wide significant signals for all the studied traits. The obtained summary statistics collected served as the input for a Mendelian randomisation approach (GSMR) to infer causal associations between putative exposure and reproduction traits. The study considered conformation, production, and metabolism as exposure and reproduction as outcome. A range of 139 to 252 genome-wide significant SNPs per combination were identified as instrumental variables (IVs). Out of 156 trait combinations, 135 demonstrated statistically significant effects, thereby enabling the identification of the responsible IVs. Combinations of traits related to metabolism (38 out of 39), conformation (68 out of 78), or production (29 out of 39) were found to have significant effects on reproduction. These relationships were partially non-linear. Moreover, a separate variance component estimation supported these findings, strongly correlating with the GSMR results and offering suggestions for improvement. Downstream analyses of selected representative traits per complex resulted in identifying and investigating potential physiological mechanisms. Notably, we identified both trait-specific SNPs and genes that appeared to influence specific traits per complex, as well as more general SNPs that were common between exposure and outcome traits.

CONCLUSIONS

Our study confirms the known genetic associations between reproduction traits and the three complexes tested. It provides new insights into causality, indicating a non-linear relationship between conformation and reproduction. In addition, the downstream analyses have identified several clustered genes that may mediate this association.

Application of Genomic Selection in Beef Cattle Disease Prevention

Genomic applications in beef cattle disease prevention have gained traction in recent years, offering new strategies for improving herd health and reducing economic losses in the livestock industry. Advances in genomics, including identification of genetic markers linked to disease resistance, provide powerful tools for early detection, selection, and management of cattle resistant to infectious diseases. By incorporating genomic technologies such as whole-genome sequencing, genotyping, and transcriptomics, researchers can identify specific genetic variants associated with resistance to pathogens like bovine respiratory disease and Johne's disease. These genomic insights allow for more accurate breeding programs aimed at enhancing disease resistance and overall herd resilience. Genomic selection, in particular, enables identification of individuals with superior genetic traits for immune function, reducing the need for antibiotic treatments and improving animal welfare. Moreover, precision medicine, powered by genomic data, supports development of tailored health management strategies, including targeted vaccination plans and antimicrobial stewardship. Incorporation of genomic tools in beef cattle management also offers the potential for early disease detection, facilitating proactive interventions that reduce the spread of infections. Despite challenges like cost, data interpretation and integration into current management systems, the potential advantages of genomic applications in disease prevention are substantial. As these technologies advance, they are anticipated to have crucial roles in improving sustainability (by enhancing herd performance), profitability (by improving overall herd longevity), and biosecurity (by decreasing the likelihood of disease outbreaks) of beef cattle production systems worldwide.

Nucleotide sequence variants, gene expression and serum profile of immune and antioxidant markers associated with brucellosis resistance/susceptibility in Shami goat

Brucellosis is a highly contagious zoonotic bacterial disease. It has considerable negative consequences on the animal production industry worldwide. The objective of this study was to investigate the genetic and molecular variations in Shami goat susceptible to Brucella infection. Blood samples were collected from fifty mature Shami goats (30 Brucella-infected does and 20 non-infection). DNA was extracted and selected parts the immunity; solute carrier family 11 member 1 (SLC11A1), toll-like receptor 1 (TLR1), toll-like receptor 9 (TLR9), SP110 nuclear body protein (SP110), the adenosine A3 receptor (ADORA3), caspase activating recruitment domain 15 (CARD15) and interferon regulatory factor 3 (IRF3), antioxidant glutathione peroxidase 1 (GPX1), nitric oxide synthase (NOS), NAD(P)H dehydrogenase [quinone] 1 (NQO1) and transcription factor NF-E2-related factor 2 (Nrf2) and erythritol related transketolase (TKT), ribose 5-phosphate isomerase (RPIA) and Adenosine monophosphate deaminase (AMPD) genes were sequenced. Likewise, the levels of gene expressions were investigated. The results identified polymorphic variants between healthy and infected does. Levels of gene expression of SLC11A1, TLR1, TLR9, SP110, ADORA3, CARD15, IRF3, HMOX1, TKT, RPIA and AMPD were significantly (P < 0.05) up regulated in the infected compared to the non-infected ones. On the other hand, GPX1, NOS, NQO1 and Nrf2 genes were significantly (P < 0.05) downregulated in the infected compared to the non-infected does. The results of serum profile indicated that there is a significant (P < 0.05) increase in the activities of AST, ALT, GGT, LDH, ALP as well as serum level of globulin, triglycerides, cholesterol, MDA, NO, IL-1β, TNF-α, IgM, IgG, haptoglobin and amyloid A. On the other hand, there were significant reductions in the glucose, total protein albumin, urea, calcium, inorganic phosphorus, sodium, copper, zinc, iron, TAC, GSH, SOD, GPx, IL-10 and fibrinogen in the infected compared to the non-infected does. Our results provide valuable information about the serum profile variations and putative genetic markers for Brucella infection in goats. This could be utilized in controlling goat brucellosis through selective breeding of natural resistant animals.

Genetic Variation of Hypothalamic-Pituitary-Adrenal Axis Activity in Farm Animals and Beyond

BACKGROUND

Many experimental data in several species clearly demonstrate the important genetic contribution to variations in HPA axis activity. The influence of corticosteroid hormones on adaptive processes and on production traits such as growth rate, feed efficiency, carcass composition, and meat quality is a strong impetus to the search for the molecular bases of these differences for efficient genetic selection.

SUMMARY

Three main sources of genetic variability have been documented so far in farm animal species, the adrenal cortex sensitivity to ACTH-regulating corticosteroid hormone production, the bioavailability of corticosteroid hormones and especially corticosteroid-binding globulin capacity, and glucocorticoid receptor function. The effect of single mutations may be dependent on the genetic background, and genetic variation of cortisol levels may have different functional consequences depending on the molecular mechanisms responsible for this change.

KEY MESSAGES

Understanding the genetic basis of HPA axis activity allows the development of genomic tools and breeding technologies aimed at improving adaptive capacity and stress tolerance in farm animals and their use as valuable models for the genetic study of the HPA axis and the correlation with adaptation, metabolism, and other functions regulated by adrenal hormones, and associated pathologies (obesity, cardiovascular, etc.). The next step will be to explore HPA axis variability from a system genetics perspective including the multiple sources of variation and their interactions. This multifactorial approach is a prerequisite to the use of the HPA axis phenotypes in the genetic selection for more productive and robust animals, with a high level of production of quality products.

Nucleotide sequence variants, gene expression and serum profile of immune and antioxidant markers associated with bacterial diarrhea susceptibility in Barki lambs

BACKGROUND

Despite the fact that diarrhea is more accurately described as a clinical symptom than a disease. Diarrhea is one of the most important issues in ovine medicine, particularly in lambs, and because of high morbidity and mortality rate, sluggish growth performance, and veterinary costs, it is believed to be a major source of economic loss. Salmonella and enterotoxigenic Escherichia coli are the most common and commercially significant agents responsible for diarrhea.

OBJECTIVE

The objective of this study was to monitor the nucleotide sequence variations, gene expression, serum inflammatory and oxidative stress biomarkers in diarrheic lambs. Another aim was to identify different pathotypes and virulence genes of Salmonella and E. coli causing diarrhea.

METHODOLOGY

Blood samples were taken from 50 Barki who were diarrheal and 50 who appeared to be healthy, and then divided in 3 portions, with EDTA added to the first part for CBC, DNA and RNA extraction. The second sample received 5000 I.U. of heparin calcium, and a clean plain tube was used for the third component. The second and third sections were centrifuged to extract serum and plasma until the biochemical and immunological analysis was completed. Fecal samples were collected for bacteriological examination, and the bacteria were identified by PCR analysis. PCR-DNA sequencing was conducted for immune (SELL, JAK2, SLC11A1, IL10, FEZF1, NCF4, LITAF, SBD2, NFKB, TNF-α, IL1B, IL6, LGALS, and CATH1), antioxidant (SOD1, CAT, GPX1, GST, Nrf2, Keap1, HMOX1, and NQO1), and GIT health (CALB1, GT, and MUC2) genes in healthy and diarrheic lambs.

RESULTS

Virulent genetic markers of pathogenic characteristics of E. coli (astA, Vt2e (Stx2e), CFA/I, groES and luxS) and Salmonella (invA, SopB, bcfC and avrA) were detected in all diarrheic lambs. PCR-DNA sequencing of immune, antioxidant and intestinal health genes found eleven single nucleotide polymorphisms (SNPs) linked to either diarrhea resistance or susceptibility in Barki lambs. Transcript levels of immune, antioxidant, and GIT health (CALB1, GT, and MUC2) genes varied between healthy and diarrheic lambs. Nucleotide sequence variation of the genes under inquiry between reference sequences in GenBank and those of the animals under investigation verified all identified SNPs. Significant (P = 0.001) erythrocytosis, neutrophilic leukocytosis, with lymphocytopenia were observed in diarrheic lambs. Significant (P = 0.001) increases in serum IL-1α, IL-1β, IL-6, TNF-α (90.5 ± 1.7, 101.8 ± 1.7, 72.3 ± 6.6, 71.26 ± 4.89 Pg/ml, respectively), serum Fb, Cp, Hp, SAA (230.7 ± 12.4 mg/dl, 6.5 ± 0.07 mg/dl, 2.5 ± 0.09 g/dl, 7.4 ± 0.4 mg/L, respectively), free radicals (MDA, NO), cortisol (6.91 ± 0.18 μg/dl) and growth hormone, with significant (P = 0.001) decreases in serum IL-10 (81.71 ± 1.05 Pg/ml), antioxidants (CAT, GPx), insulin, triiodothyronine (T3) and thyroxine (T4) in diarrheic lambs.

CONCLUSIONS

The study's findings provided credence to the theory that marker-assisted selection (MAS) could be used to predict and prevent diarrhea in Barki sheep by selecting lambs based on SNPs in genes linked to inflammation, antioxidants, and intestinal health. In order to establish an efficient management protocol and determine the most susceptible risk period for disease occurrence, gene expression profiles of the genes under investigation, pro-inflammatory cytokines and acute phase proteins may also be utilized as proxy biomarkers for lamb enteritis.

Analysis of potential genes, immunological and antioxidant profiles associated with trypanosomiasis susceptibility in dromedary camels

Trypanosomiasis is associated with tissue damage and may trigger an immunological response. These tissue lesions are linked to metabolic issues and oxidative stress. The current study aimed to investigate the immunological, antioxidant, and metabolic changes that may be connected to camel trypanosomiasis. Blood samples were collected from 54 camels and allocated into two groups: The control group (35 camels) and the infected group (19 camels). The genes TLR2, TLR5, IL-17, MARCHF3, RASGRP1, EPS15L1, PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, DIO3, keap1, and OXSR1 were significantly up-regulated in trypanosomiasis camels. However, down-regulation was observed for the genes Nrf2, PRDX6, and NDUFS5. PCR-DNA sequencing was used to identify nucleotide sequence polymorphisms in the immune (TLR2, TLR5, IL-17, MARCHF3, RASGRP1, and EPS15L1), metabolic (PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, and DIO3), and antioxidant (Nrf2, Keap1, PRDX6, NDUFS5, and OXSR1) genes between healthy and trypanosomiasis-affected camels. Exploring the serum profile also showed a significant (P ˂ 0.05) increase in Hp, SAA, Cp, IL-1β, IL-6, IL 10, TNF-α, and MDA, with significant (P ˂ 0.05) reduction in the serum levels of CAT, SOD, GSH, T3, and T4 in diseased camels compared with healthy ones. Our findings confirm the significance of nucleotide variations, gene expression patterns, and the biochemical profile of the investigated markers as indicators for the susceptibility of trypanosomiasis in dromedary camels and may be utilized to create management strategies.

The Transcript Levels and the Serum Profile of Biomarkers Associated with Clinical Endometritis Susceptibility in Buffalo Cows

Determining the gene expression and serum profile of the indicators linked to clinical endometritis susceptibility in Egyptian buffalo cows was the aim of this investigation. The buffalo cows that were enrolled were divided into two groups: forty infected buffalo cows with clinical endometritis and forty seemingly healthy buffalo cows that served as the control group. For the purposes of gene expression and biochemical analysis, ten milliliters of blood was obtained via jugular venipuncture from each buffalo cow. TLR4, IL-8, IL-17, NFKB, SLCA11A1, NCF4, Keap1, HMOX1, OXSR1, ST1P1, and SERP1 were manifestly expressed at much higher levels in the buffaloes with endometritis. On the other hand, the genes that encode SOD, CAT, NDUFS6, Nrf2, and PRDX2 were down-regulated. There was a significant (p < 0.05) elevation of the serum levels of non-esterified fatty acids (NEFAs), beta hydroxy butyric acid (BHBA), triglycerides (TGs), globulin, creatinine, and cortisol, along with a reduction in the serum levels of glucose, cholesterol, total protein albumin, urea, estrogen (E2), progesterone (P4), follicle-stimulating hormone (FSH), luteinizing hormone (LH), thyroxine (T4), prostaglandin F2 α (PGF2α), calcium, iron, and selenium, in the endometritis group in comparison with the control. However, no significant change was observed in the values of phosphorus, magnesium, copper, or zinc in either group. Within the selective breeding of naturally resistant animals, the variation in the genes under study and the changes in the serum profiles of the indicators under investigation may serve as a reference guide for reducing endometritis in Egyptian buffalo cows.

Doppler ultrasonographic scan, gene expression and serum profile of immune, APPs and antioxidant markers in Egyptian buffalo-cows with clinical endometritis

The objective of this study was to elaborate Doppler ultrasonographic scan, genetic resistance and serum profile of markers associated with endometritis susceptibility in Egyptian buffalo-cows. The enrolled animals were designed as; twenty five apparently healthy buffalo-cows considered as a control group and twenty five infected buffalo with endometritis. There were significant (p < 0.05) increased of cervical diameter, endometrium thickness, uterine horn diameter, TAMEAN, TAMAX and blood flow through middle uterine artery with significant decrease of PI and RI values in endometritis buffalo-cows. Gene expression levels were considerably higher in endometritis-affected buffaloes than in resistant ones for the genes A2M, ADAMTS20, KCNT2, MAP3K4, MAPK14, FKBP5, FCAMR, TLR2, IRAK3, CCl2, EPHA4, and iNOS. The RXFP1, NDUFS5, TGF-β, SOD3, CAT, and GPX genes were expressed at substantially lower levels in endometritis-affected buffaloes. The PCR-DNA sequence verdicts of healthy and affected buffaloes revealed differences in the SNPs in the amplified DNA bases related to endometritis for the investigated genes. However, MAP3K4 elicited a monomorphic pattern. There was a significant decrease of red blood cells (RBCs) count, Hb and packed cell volume (PCV) with neutrophilia, lymphocytosis and monocytosis in endometritis group compared with healthy ones. The serum levels of Hp, SAA, Cp, IL-6, IL-10, TNF-α, NO and MDA were significantly (P˂0.05) increased, along with reduction of CAT, GPx, SOD and TAC in buffalo-cows with endometritis compared to healthy ones. The variability of Doppler ultrasonographic scan and studied genes alongside alterations in the serum profile of investigated markers could be a reference guide for limiting buffalo endometritis through selective breeding of natural resistant animals.

Veterinary systems biology for bridging the phenotype-genotype gap via computational modeling for disease epidemiology and animal welfare

Veterinary systems biology is an innovative approach that integrates biological data at the molecular and cellular levels, allowing for a more extensive understanding of the interactions and functions of complex biological systems in livestock and veterinary science. It has tremendous potential to integrate multi-omics data with the support of vetinformatics resources for bridging the phenotype-genotype gap via computational modeling. To understand the dynamic behaviors of complex systems, computational models are frequently used. It facilitates a comprehensive understanding of how a host system defends itself against a pathogen attack or operates when the pathogen compromises the host's immune system. In this context, various approaches, such as systems immunology, network pharmacology, vaccinology and immunoinformatics, can be employed to effectively investigate vaccines and drugs. By utilizing this approach, we can ensure the health of livestock. This is beneficial not only for animal welfare but also for human health and environmental well-being. Therefore, the current review offers a detailed summary of systems biology advancements utilized in veterinary sciences, demonstrating the potential of the holistic approach in disease epidemiology, animal welfare and productivity.

New insights into the genetic predisposition of brucellosis and its effect on the gut and vaginal microbiota in goats

Goats contribute significantly to the global food security and industry. They constitute a main supplier of meat and milk for large proportions of people in Egypt and worldwide. Brucellosis is a zoonotic infectious disease that causes a significant economic loss in animal production. A case-control genome-wide association analysis (GWAS) was conducted using the infectious status of the animal as a phenotype. The does that showed abortion during the last third period of pregnancy and which were positive to both rose bengal plate and serum tube agglutination tests, were considered as cases. Otherwise, they were considered as controls. All animals were genotyped using the Illumina 65KSNP BeadChip. Additionally, the diversity and composition of vaginal and fecal microbiota in cases and controls were investigated using PCR-amplicone sequencing of the V4 region of 16S rDNA. After applying quality control criteria, 35,818 markers and 66 does were available for the GWAS test. The GWAS revealed a significantly associated SNP (P = 5.01 × 10-7) located on Caprine chromosome 15 at 29 megabases. Four other markers surpassed the proposed threshold (P = 2.5 × 10-5). Additionally, fourteen genomic regions accounted for more than 0.1% of the variance explained by all genome windows. Corresponding markers were located within or in close vicinity to several candidate genes, such as ARRB1, RELT, ATG16L2, IGSF21, UBR4, ULK1, DCN, MAPB1, NAIP, CD26, IFIH1, NDFIP2, DOK4, MAF, IL2RB, USP18, ARID5A, ZAP70, CNTN5, PIK3AP1, DNTT, BLNK, and NHLRC3. These genes play important roles in the regulation of immune responses to the infections through several biological pathways. Similar vaginal bacterial community was observed in both cases and controls while the fecal bacterial composition and diversity differed between the groups (P < 0.05). Faeces from the control does showed a higher relative abundance of the phylum Bacteroidota compared to cases (P < 0.05), while the latter showed more Firmicutes, Spirochaetota, Planctomycetota, and Proteobacteria. On the genus level, the control does exhibited higher abundances of Rikenellaceae RC9 gut group and Christensenellaceae R-7 group (P < 0.05), while the infected does revealed higher Bacteroides, Alistipes, and Prevotellaceae UCG-003 (P < 0.05). This information increases our understanding of the genetics of the susceptibility to Brucella in goats and may be useful in breeding programs and selection schemes that aim at controlling the disease in livestock.

Invited Review: Increasing Milk Yield and Negative Energy Balance: A Gordian Knot for Dairy Cows?

The continued increase in milk production during the last century has not been accompanied by an adequate dry matter intake (DMI) by cows, which therefore experience a negative energy balance (NEB). NEB is low and of minor importance at low milk yield (MY), such as for the nutrition of one calf, and under these circumstances is considered "natural". MY and low DMI around parturition are correlated and are the reason for the genetic correlation between increasing MY and increasing NEB up to 2000 MJ or more for 2-3 months postpartum in high-genetic-merit dairy cows. The extension and duration of NEB in high-producing cows cannot be judged as "natural" and are compensated by the mobilization of nutrients, particularly of fat. The released non-esterified fatty acids (NEFAs) overwhelm the metabolic capacity of the cow and lead to the ectopic deposition of NEFAs as triglycerides (TGs) in the liver. The subsequent lipidosis and the concomitant hampered liver functions cause subclinical and clinical ketosis, both of which are associated with "production diseases", including oxidative and endoplasmatic stress, inflammation and immunosuppression. These metabolic alterations are regulated by homeorhesis, with the priority of the physiological function of milk production. The prioritization of one function, namely, milk yield, possibly results in restrictions in other physiological (health) functions under conditions of limited resources (NEB). The hormonal framework for this metabolic environment is the high concentration of growth hormone (GH), the low concentration of insulin in connection with GH-dependent insulin resistance and the low concentration of IGF-1, the so-called GH-IGF-1 axis. The fine tuning of the GH-IGF-1 axis is uncoupled because the expression of the growth hormone receptor (GHR-1A) in the liver is reduced with increasing MY. The uncoupled GH-IGF-1 axis is a serious impairment for the GH-dependent stimulation of gluconeogenesis in the liver with continued increased lipolysis in fat tissue. It facilitates the pathogenesis of lipidosis with ketosis and, secondarily, "production diseases". Unfortunately, MY is still increasing at inadequate DMI with increasing NEB and elevated NEFA and beta-hydroxybutyric acid concentrations under conditions of low glucose, thereby adding health risks. The high incidences of diseases and of early culling and mortality in dairy cows are well documented and cause severe economic problems with a waste of resources and a challenge to the environment. Moreover, the growing public concerns about such production conditions in agriculture can no longer be ignored.

Assessment of Genetic and Health Management of Tunisian Holstein Dairy Herds with a Focus on Longevity

In Tunisia, the recognition of the possibility of including longevity and disease resistance in dairy cattle selection objectives has been hypothesized as a useful strategy by both researchers and producers. However, in this paper, the state of the art, with a focus on health and longevity, is reviewed. Along the same lines, the heritability for the milk traits, fertility traits, and longevity of Tunisian Holstein dairy cows complies with the literature. Therefore, the influence of genetics on some diseases of the dairy cow was investigated. In addition, a decreasing efficiency in cow fertility has been observed over the last few years. The results showed that the risk of culling increased with common diseases. When analyzed with the Weibull model, functional lifespan was strongly influenced by milk yield; therefore, the risk increased with a reduced milk yield. In her first three lactations, the relative risk of selection increased gradually with lactation. Thus, the risk of thinning is highest at the beginning and end of the first feeding and the end of her second feeding. In conclusion, the risk of culling was reduced in parity. The factors that influence the life of the herd, such as health, husbandry, environmental conditions, and management, are often ignored when evaluating longevity.

Individual Genomic Loci and mRNA Levels of Immune Biomarkers Associated with Pneumonia Susceptibility in Baladi Goats

The effectiveness of breeding for inherent disease resistance in animals could be considerably increased by identifying the genes and mutations that cause diversity in disease resistance. One hundred and twenty adult female Baladi goats (sixty pneumonic and sixty apparently healthy) were used in this study. DNA and RNA were extracted from blood samples collected from the jugular vein of each goat. SLC11A1, CD-14, CCL2, TLR1, TLR7, TLR8, TLR9, β defensin, SP110, SPP1, BP1, A2M, ADORA3, CARD15, IRF3, and SCART1 SNPs that have been previously found to be associated with pneumonia resistance/susceptibility were identified via PCR-DNA sequencing. The pneumonic and healthy goats differed significantly, according to a Chi-square analysis of the discovered SNPs. The mRNA levels of the studied immune markers were noticeably greater in the pneumonic goats than in the healthy ones. The findings could support the significance of the use of immune gene expression profiles and nucleotide variations as biomarkers for the susceptibility/resistance to pneumonia and provide a practical management technique for Baladi goats. These results also suggest a potential strategy for lowering pneumonia in goats by employing genetic markers linked to an animal’s ability to fend off infection in selective breeding.

Progress and opportunities through use of genomics in animal production

The rearing of farmed animals is a vital component of global food production systems, but its impact on the environment, human health, animal welfare, and biodiversity is being increasingly challenged. Developments in genetic and genomic technologies have had a key role in improving the productivity of farmed animals for decades. Advances in genome sequencing, annotation, and editing offer a means not only to continue that trend, but also, when combined with advanced data collection, analytics, cloud computing, appropriate infrastructure, and regulation, to take precision livestock farming (PLF) and conservation to an advanced level. Such an approach could generate substantial additional benefits in terms of reducing use of resources, health treatments, and environmental impact, while also improving animal health and welfare.

Genetic Resistance of Bovines to Theileriosis

Diseases caused by ticks have a high impact on the health, welfare, and productivity of livestock species. They are also an important cause of economic losses in farms worldwide. An example of such diseases is theileriosis, which can be controlled by drugs or vaccines, although these are not fully efficient. Therefore, there is a need to develop alternative and more sustainable and efficient complementary strategies. These may involve the identification and selection of animals more resistant to the disease. Several previous studies have identified significant differences in resistance between different breeds, with resistant breeds typically identified as those native to the region where they are being studied, and susceptible as those from exotic breeds. These studies have indicated that resistance traits are intrinsically related to the modulation of the immune response to infection. This review aims to systematize the general knowledge about theileriosis, emphasize resistance to this disease as a sustainable control strategy, and identify which traits of resistance to the disease are already known in cattle.

Keeping profitable cows in the herd: A lifetime cost-benefit assessment to support culling decisions

Increasing the productive lifespan of dairy cows is important to achieve a sustainable dairy industry, but making strategic culling decisions based on cow profitability is challenging for farmers. The objective of this study was to carry out a lifetime cost-benefit analysis based on production and health records and to explore different culling decisions among farmers. The cost-benefit analysis was conducted for 22 747 dairy cows across 114 herds in Quebec, Canada for which feed costs and the occurrence of diseases were reported. Costs and revenues related to productive lifespan were compared among cohorts of cows that left their respective herd at the end of their last completed lactation or stayed for a complete additional lactation. Hierarchical clustering analysis was carried out based on costs and revenues to explore different culling decisions among farmers. Our results showed that the knowledge of lifetime cumulative costs and revenues was of great importance to identify low-profitable cows at an earlier lactation, while only focusing on current lactation costs and revenues can lead to an erroneous assessment of profitability. While culling decisions were mostly based on current lactation costs and revenues and disregarded the occurrence of costly events on previous lactations, there was variation among farmers as we identified three different culling decision clusters. Monitoring cumulative costs and revenues would help farmers to identify low-profitable cows at an earlier lactation and make the decision to increase herd productive lifespan and farm profitability by keeping the most profitable cows.

Repeatability of Health and Welfare Traits and Correlation with Performance Traits in Dairy Goats Reared under Low-Input Farming Systems

The objectives of the study were to estimate the repeatability of health and welfare traits and investigate their association with performance in three breeds of dairy goats reared under low-input farming systems in Greece. A total of 1210 goats of Eghoria (n = 418), Skopelos (n = 429), and Damascus (n = 363) breeds were assessed. Udder health, parasitic resistance, welfare, milk yield and quality, and body condition score were recorded monthly for two milking periods. Udder health records included somatic cell count (SCC) and total viable count (TVC). Based on combinations of SCC and TVC and thresholds set at >106 cells/mL and >2 × 104 cfu/mL, respectively, additional udder health phenotypes were defined. Parasitism included myiasis, tick infestation, gastrointestinal nematode (GIN) and cestode faecal egg count (FEC), and lungworm faecal larval count (FLC). Infection with each of the endoparasites was defined based on FEC/FLC. Welfare assessment parameters included the presence of ear and horn injuries, ocular and nasal discharge, body and udder abscesses, injury and lesions on the skin of different regions, diarrhoea, hernias, overgrown hooves, arthritis, lameness, and udder asymmetry. Trait repeatability and animal correlations were estimated. Significant (p < 0.05) repeatability was reported for all udder health and most welfare traits in all breeds, GIN and cestode FEC, and GIN and lungworm infection in Eghoria, and myiasis in Skopelos. Correlations of health and most of welfare traits with performance were non-significant or favourable. Overall, results demonstrate potential to improve health and welfare of the studied breeds without compromising performance.

Whole-Genome comparative analysis reveals genetic mechanisms of disease resistance and heat tolerance of tropical Bos indicus cattle breeds

The Bos indicus cattle breeds have been naturally selected over thousands of years for disease resistance and thermo-tolerance. However, a genetic mechanism of these specific inherited characteristics needs to be discovered. Hence, in this study, the whole-genome comparative analysis of Bos indicus cattle breeds of Kangayam, Tharparkar, Sahiwal, Red Sindhi, and Hariana of the Indian subcontinent was conducted. The genetic variants identification analysis revealed a total of 15,58,51,012 SNPs and 1,00,62,805 InDels in the mapped reads across all Bos indicus cattle breeds. The functional annotation of 17,252 genes that comprised both, SNPs and InDels, of high functional impact on proteins, has been carried out. The functional annotation results revealed the pathways that were involved in the innate immune response including toll-like receptors, a retinoic acid-inducible gene I like receptors, NOD-like receptors, Jak-STAT signaling pathways, and the non-synonymous variants in the candidate immune genes. Further, we also identified several pathways involved in heat shock response, hair and skin properties, oxidative stress response, osmotic stress response, thermal sweating, feed intake, metabolism, and the non-synonymous variants in the candidate thermo-tolerant genes. These pathways and genes were directly or indirectly contributing to the disease resistance and thermo-tolerance adaptations of Bos indicus cattle breeds.

Genome-wise engineering of ruminant nutrition- nutrigenomics: applications, challenges, and future perspectives – a review

Use of genomic information in ruminant production systems can help relieve concerns related to food security and sustainability of production. Nutritional genomics (i.e., Nutrigenomics) is a field of research that is interested in all types of reciprocal interactions between nutrients and genomes of organisms, i.e., variable patterns of gene expression and effect of genetic variations on the nutritional environment. Devising a revolutionizing analytical approach to traditional ruminant nutrition research, the relatively novel area of ruminant nutrigenomics has several studies concerning different aspects of animal production systems. This paper aims to review the current nutrigenomics research in the frame of how nutrition of ruminants can be modified accounting for individual genetic backgrounds and gene/diet relationships behind productivity, quality, efficiency, disease resistance, fertility, and GHG emissions. Furthermore, current challenges facing ruminant nutrigenomics are evaluated and future directions for the novel area are strongly argued by this review.

Single nucleotide polymorphisms, gene expression and serum profile of immune and antioxidant markers associated with postpartum disorders susceptibility in Barki sheep

The objective of this study was to explore the immunological and antioxidant alterations associated with ovine postpartum disorders. Blood samples were collected from 90 adult Barki ewes and allocated into three equal-sized groups (30 ewes each): control group (CG), inflammatory postpartum disorders group (IPG) and non-inflammatory postpartum disorders group (NIPG). PCR-DNA sequencing approach was carried out for TLR4 (256-bp) and SOD (456-bp) genes, and nucleotide sequence variations were noticed to be associated with postpartum disorders resistance/susceptibility. Gene expression profile was also evaluated and levels of IL5, IL6, IL1-ß, TNF alpha, TLR4 and Tollip were significantly up-regulated in ewes affected with postpartum disorders than resistant ones, while SOD and CAT genes pattern elicited an opposite trend. Exploring serum profile also showed a significant increase of IL-1α, IL-1β, IL-6, TNF-α, MDA and NO in IPG compared to their correspond values in NIPG and CG. However, serum levels of IL-10, CAT, GSH and GPx were significantly decreased. This study highlights that SNPs in TLR4 and SOD genes could be genetic markers for postpartum disorders resistance/susceptibility in Barki ewes. Gene expression alongside serum profiles of antioxidant markers could also be used to follow-up the immune status of ewes to build up an effective management protocol.

Including milk production, conformation, and functional traits in multivariate models for genetic evaluation of lameness

Lameness is a serious health and welfare issue that can negatively affect the economic performance of cows, especially on pasture-based dairy farms. However, most genetic predictions (GP) of lameness have low accuracy because lameness data are often incomplete as data are collected voluntarily by farmers in countries such as Australia. The objective of this study was to find routinely measured traits that are correlated with lameness and use them in multivariate evaluation models to improve the accuracy of GP for lameness. We used health events and treatments associated with lameness recorded by Australian farmers from 2002 to early 2019. The lameness incidence rates in Holstein and Jersey cows were 3.3% and 4.6%, respectively. We analyzed the records of 36 other traits (milk production, conformation, fertility, and survival traits) to estimate genetic correlations with lameness. The estimated heritability ± standard error (and repeatability ± standard error) for lameness in both Holstein and Jersey breeds were very low: 0.007 ± 0.002 (and 0.029 ± 0.002) and 0.005 ± 0.003 (and 0.027 ± 0.006), respectively, in univariate sire models. For the GP models, we tested including measurements of overall type to prediction models for Holsteins, stature and body length for Jersey, and milk yield and fertility traits for both breeds. The average accuracy of GP, calculated from prediction error variances, were 0.38 and 0.24 for Holstein and Jersey sires, respectively, when estimated using univariate sire models and both increased to 0.43 using multivariate sire models. In conclusion, we found that the accuracy of GP for lameness could be improved by including genetically correlated traits in a multivariate model. However, to further improve the accuracy of predictions of lameness, precise identification and recording incidences of hoof or leg disorder, or large-scale recording of locomotion and claw scores by trained personnel should be considered.

Genotype by heat stress interactions for production and functional traits in dairy cows from an across-generation perspective

The aim of this study was to analyze time-lagged heat stress (HS) effects during late gestation on genetic co(variance) components in dairy cattle across generations for production, female fertility, and health traits. The data set for production and female fertility traits considered 162,492 Holstein Friesian cows from calving years 2003 to 2012, kept in medium-sized family farms. The health data set included 69,986 cows from calving years 2008 to 2016, kept in participating large-scale co-operator herds. Production traits were milk yield (MKG), fat percentage (fat%), and somatic cell score (SCS) from the first official test-day in first lactation. Female fertility traits were the nonreturn rate after 56 d (NRR56) in heifers and the interval from calving to first insemination (ICFI) in first-parity cows. Health traits included clinical mastitis (MAST), digital dermatitis (DD), and endometritis (EM) in the early lactation period in first-parity cows. Meteorological data included temperature and humidity from public weather stations in closest herd distance. The HS indicator was the temperature-humidity index (THI) during dams' late gestation, also defined as in utero HS. For the genetic analyses of production, female fertility, and health traits in the offspring generation, a sire-maternal grandsire random regression model with Legendre polynomials of order 3 for the production and of order 2 for the fertility and health traits on prenatal THI, was applied. All statistical models additionally considered a random maternal effect. THI from late gestation (i.e., prenatal climate conditions), influenced genetic parameter estimates in the offspring generation. For MKG, heritabilities and additive genetic variances decreased in a wave-like pattern with increasing THI. Especially for THI >58, the decrease was very obvious with a minimal heritability of 0.08. For fat% and SCS, heritabilities increased slightly subjected to prenatal HS conditions at THI >67. The ICFI heritabilities differed marginally across THI [heritability (h²) = 0.02-0.04]. For NRR56, MAST, and DD, curves for heritabilities and genetic variances were U-shaped, with largest estimates at the extreme ends of the THI scale. For EM, heritability increased from THI 25 (h² = 0.13) to THI 71 (h² = 0.39). The trait-specific alterations of genetic parameters along the THI gradient indicate pronounced genetic differentiation due to intrauterine HS for NRR56, MAST, DD, and EM, but decreasing genetic variation for MKG and ICFI. Genetic correlations smaller than 0.80 for NRR56, MAST, DD, and EM between THI 65 with corresponding traits at remaining THI indicated genotype by environment interactions. The lowest genetic correlations were identified when considering the most distant THI. For MKG, fat%, SCS, and ICFI, genetic correlations throughout were larger than 0.80, disproving concerns for any genotype by environment interactions. Variations in genetic (co)variance components across prenatal THI may be due to epigenetic modifications in the offspring genome, triggered by in utero HS. Epigenetic modifications have a persistent effect on phenotypic responses, even for traits recorded late in life. However, it is imperative to infer the underlying epigenetic mechanisms in ongoing molecular experiments.

Genetic and nongenetic factors associated with lactation length in seasonal-calving, pasture-based dairy cows

Lactation yield estimates standardized to common lactation lengths of 270-d or 305-d equivalents are commonly used in management decision support tools and dairy cow genetic evaluations. The use of such measurements to quantify the (genetic) merit of individual cows fails to penalize cows that do not reach the standardized lactation length, or indeed reward cows that lactate for more than the standardized lactation length. The objective of the present study was to quantify the genetic and nongenetic factors associated with lactation length in seasonal-calving, pasture-based dairy cows. A total of 616,350 lactation length records from 285,598 Irish cows were used. Linear mixed models were used to quantify the associations between lactation length and calving month, parity, age at calving, previous dry period length, calving difficulty score, heterosis, recombination loss, breed, and herd size, as well as to estimate the genetic and residual variance components of lactation length. The median lactation length in the edited data set was 288 d, with 27% of cows achieving lactations of at least 305 d. Relative to cows calving in January, the lactations of cow calving in February, March, or April was, on average, 4.2, 12.7, and 21.9 d shorter, respectively. The lactation length of a first parity cow was, on average, 7.8, 8.6, and 8.4 d shorter than that of second, third, and fourth parity cows, respectively. Norwegian Red and Montbéliarde cows had, on average, a 4.7- and 1.6-d shorter lactation than Holstein-Friesian cows, respectively. The heritability estimate, coefficient of genetic variation, and repeatability estimate of lactation length were 0.02, 1.2%, and 0.04, respectively. Based on the genetic standard deviation for lactation length estimated in the present study (3.3 d), cows ranked in the top 20% for genetic merit for lactation length would be expected to have lactations 9.2 d longer than cows in the bottom 20%, demonstrating exploitable genetic variability. Given the vast array of genetic and nongenetic factors associated with lactation length, an approach which combines improved management practices and selective breeding may be an efficient and effective strategy to lengthen lactations.

Genome wide association study of passive immunity and disease traits in beef-suckler and dairy calves on Irish farms

Calves with lower concentrations of immunoglobulin G (IgG) in their blood, have a greater risk of developing diseases. There is a lack of knowledge on genetic markers known to be associated with immunological variability or disease resistance. Therefore, the objective of this study was to identify SNP markers associated with passive immunity measures (serum IgG, serum protein, albumin, globulin and total protein concentrations, total solids Brix percentage, zinc sulphate turbidity units) and disease (pneumonia, diarrhoea, crude illness) traits in Irish commercial beef-suckler and dairy calves through genome wide association studies (GWAS). Genotyping was performed on DNA samples from beef-suckler (n = 698) and dairy (n = 1178) calves, using the IDBv3 chip. Heritability of passive immunity associated traits (range 0.02-0.22) and the disease traits (range 0.03-0.20) were low-to-moderate. Twenty-five and fifteen SNPs approached genome wide significance (P < 5 × 10^-5) for the passive immunity and the disease traits, respectively. One SNP "ARS-BFGL-BAC-27914" reached Bonferroni genome wide significance (P < 1.15 × 10^-6) for an association with serum IgG concentration in beef calves. Further work will evaluate these SNPs in larger cattle populations and assess their contribution to genomic selection breeding strategies, aimed towards producing more disease resistant livestock.

Genomic Study of Babesia bovis Infection Level and Its Association With Tick Count in Hereford and Braford Cattle

Bovine babesiosis is a tick-borne disease caused by intraerythrocytic protozoa and leads to substantial economic losses for the livestock industry throughout the world. Babesia bovis is considered the most pathogenic species, which causes bovine babesiosis in Brazil. Genomic data could be used to evaluate the viability of improving resistance against B. bovis infection level (IB) through genomic selection, and, for that, knowledge of genetic parameters is needed. Furthermore, genome-wide association studies (GWAS) could be conducted to provide a better understanding of the genetic basis of the host response to B. bovis infection. No previous work in quantitative genetics of B. bovis infection was found. Thus, the objective of this study was to estimate the genetic correlation between IB and tick count (TC), evaluate predictive ability and applicability of genomic selection, and perform GWAS in Hereford and Braford cattle. The single-step genomic best linear unbiased prediction method was used, which allows the estimation of both breeding values and marker effects. Standard phenotyping was conducted for both traits. IB quantifications from the blood of 1,858 animals were carried using quantitative PCR assays. For TC, one to three subsequent tick counts were performed by manually counting adult female ticks on one side of each animal's body that was naturally exposed to ticks. Animals were genotyped using the Illumina BovineSNP50 panel. The posterior mean of IB heritability, estimated by the Bayesian animal model in a bivariate analysis, was low (0.10), and the estimations of genetic correlation between IB and TC were also low (0.15). The cross-validation genomic prediction accuracy for IB ranged from 0.18 to 0.35 and from 0.29 to 0.32 using k-means and random clustering, respectively, suggesting that genomic predictions could be used as a tool to improve genetics for IB, especially if a larger training population is developed. The top 10 single nucleotide polymorphisms from the GWAS explained 5.04% of total genetic variance for IB, which were located on chromosomes 1, 2, 5, 6, 12, 17, 18, 16, 24, and 26. Some candidate genes participate in immunity system pathways indicating that those genes are involved in resistance to B. bovis in cattle. Although the genetic correlation between IB and TC was weak, some candidate genes for IB were also reported in tick infestation studies, and they were also involved in biological resistance processes. This study contributes to improving genetic knowledge regarding infection by B. bovis in cattle.

How mechanistic modelling supports decision making for the control of enzootic infectious diseases

Controlling enzootic diseases, which generate a large cumulative burden and are often unregulated, is needed for sustainable farming, competitive agri-food chains, and veterinary public health. We discuss the benefits and challenges of mechanistic epidemiological modelling for livestock enzootics, with particular emphasis on the need for interdisciplinary approaches. We focus on issues arising when modelling pathogen spread at various scales (from farm to the region) to better assess disease control and propose targeted options. We discuss in particular the inclusion of farmers' strategic decision-making, the integration of within-host scale to refine intervention targeting, and the need to ground models on data.

European perspectives on efforts to reduce antimicrobial usage in food animal production

New regulations on veterinary medicines and medicated feed will substantially influence antimicrobial prescribing and usage throughout Europe into the future. These regulations have been informed by a very large body of work, including the substantial progress towards reduced antimicrobial usage in food animal production in a number of member states of the European Union (EU). This paper seeks to summarise European perspectives on efforts to reduce antimicrobial usage in food animal production. Work within the EU is informed by the global action plan of the World Health Organization, which includes a strategic objective to optimise the use of antimicrobial medicines in human and animal health. There is ongoing measurement of trends in antimicrobial usage and resistance throughout the EU, and detailed information on strategies to reduce the need to use antimicrobials in food animal production. Substantial scientific progress has been made on the measurement of antimicrobial usage, including at herd-level, and on the objective assessment of farm biosecurity. In a number of EU member states, monitoring systems for usage are well-established, allowing benchmarking for veterinarians and farms, and monitoring of national and industry-level trends. Several countries have introduced restrictions on antimicrobial prescribing and usage, including strategies to limit conflicts of interest around antimicrobial prescribing and usage. Further, a broad range of measures are being used across member states to reduce the need for antimicrobial usage in food animal production, focusing both at farm level and nationally. Veterinarians play a central role in the reduction of antimicrobial usage in farm animals. Ireland's National Action Plan on Antimicrobial Resistance 2017-20 (iNAP) provides an overview of Ireland's commitment to the development and implementation of a holistic, cross-sectoral 'One Health' approach to the problem of antimicrobial resistance. The new regulations offer an important springboard for further progress, in order to preserve the efficacy of existing antimicrobials, which are a critical international resource.

Human Relationships with Domestic and Other Animals: One Health, One Welfare, One Biology

Excessive human population growth, uncontrolled use of natural resources, including deforestation, mining, wasteful systems, biodiversity reduction by agriculture, and damaging climate change affect the existence of all animals, including humans. This discussion is now urgent and people are rethinking their links with the animals we use for clothing, food, work, companionship, entertainment, and research. The concepts of one health, one welfare, and one biology are discussed as a background to driving global change. Nothing should be exploited without considering the ethics of the action and the consequences. This review concerns domesticated animals, including those used for human consumption of meat, eggs, and milk; horses kept for work; and dogs kept for company. Animal welfare includes health, emotional state, and comfort while moving and resting, and is affected by possibilities to show behavior and relationships with others of the same species or with humans. We show some examples of the relations between humans and domesticated animals in the environmental context, including zoonotic diseases, and consider the consequences and the new paradigms resulting from current awareness.

Economics of reducing antibiotic usage for clinical mastitis and metritis through genomic selection

Antibiotics use (ABU) in animal agriculture has been implicated in the emergence of antibiotic resistance, a global public health threat. Economically justifiable antibiotic reduction strategies can motivate farmers to reduce ABU for clinical mastitis (CM) and metritis, the most common reasons for ABU on dairy farms. Our objective was to quantify the reduction in incidence of CM, metritis, and ABU, and the herd performance of a representative US herd that uses genomic selection for Lifetime Net Merit 2018 (NM$) selection index, compared with genetic selection based only on the mastitis (MAST) or metritis resistance (METR) trait or a health trait subindex (HTH$). A stochastic dynamic simulation model of a 1,000-cow herd with multi-trait genetics for 19 correlated traits included in the NM$ affected the performance of animals. The incidence of CM and metritis for each animal was affected by the genetic and environmental components of the MAST or METR, along with a standard phenotypic function that calculated the daily underlying herd probability to contract CM or metritis. Selection decisions were made based on genomic estimated breeding values of the traits of interest. A strategy named AI-NM$ based decisions on the NM$ trait so that the correlated genetic trends in MAST and METR are improved. Three other strategies named AI-MAST, AI-METR, and AI-HTH$ maximized respectively MAST, METR, and HTH$ genetic merit, but with a tradeoff in NM$ genetic merit. The cumulative true breeding values (TBV) of NM$ for 15 yr showed a difference of $4,947 per cow between the AI-NM$ (best strategy for NM$) and AI-METR (worst strategy for NM$). However, the 15-yr cumulative TBV of MAST was 26.50 percentage points (PP) higher in AI-MAST, and 18.5 PP higher for METR in AI-METR, compared with AI-NM$. As a result, the 15-yr cumulative phenotypic CM and metritis incidence was respectively 94.03 PP and 58 PP lower in AI-MAST and AI-METR compared with AI-NM$. Therefore the corresponding 15-yr cumulative ABU decreased by 42% in AI-MAST and by 53% in AI-METR. We found that AI-MAST had the lowest CM incidence across the 15 yr, whereas AI-METR had the lowest incidence of metritis and the smallest total ABU for 15 yr. To achieve the lowest incidence of CM, metritis, and ABU strategies AI-MAST, AI-METR, and AI-HTH$ had to incur 15-yr discounted cumulative losses per cow of $1,486, $1,434, and $1,130, respectively, compared with AI-NM$. Hence, AI-NM$ had the best financial performance, despite having slightly higher incidence of CM, metritis, and ABU.

Routine Herd Health Data as Cow-Based Risk Factors Associated with Lameness in Pasture-Based, Spring Calving Irish Dairy Cows

Simple Summary

Dairy cow lameness is considered one of the most important animal welfare and economic concerns for the dairy industry. Cow-based risk factors for lameness are not well described, especially in comparison to herd-level risk factors such as housing environment and roadway condition. This study investigates the use of routinely gathered herd health data as cow-based risk factors for lameness in dairy cows. A total of 1715 cows in 10 pasture-based Irish dairy herds were monitored for lameness during the spring and summer of 2013 and 2014 as part of the University College Dublin herd health programme. Herd health monitoring data analysed to identify potential risk factors for lameness consisted of milk production data, genetic merit information, calving details, peri-parturient disease records and body condition scores. This analysis showed a significant effect of increasing parity, lower body condition score at calving and excessive body condition loss after calving on the risk of cows being diagnosed as lame during the lactation. In conclusion, routinely gathered herd health monitoring data can be used to identify cows at increased risk of being lame and to implement lameness control strategies.

Abstract

Herd-level risk factors related to the cow’s environment have been associated with lameness. Uncomfortable stall surface and inadequate depth of bedding as well as abrasive alley way surface are contributing factors to increased levels of lameness. Access to pasture has been found as having a beneficial effect on cows’ locomotion. However, dairy cattle managed under grazing conditions are exposed to a different set of risk factors for lameness, mainly associated with cow tracks. Cow-based risk factors for lameness are not as clearly defined as the herd level risk factors. The objective of the present study was to use routine herd health monitoring data to identify cow-based risk factors for lameness and utilise this information to indicate cows at risk of developing lameness in the first 150 days of lactation. Lameness data were gathered from 10 pasture-based dairy herds. A total of 1715 cows were monitored, of which 1675 cows were available for analysis. Associations between lameness status and potential cow-level risk factors were determined using multivariable logistic regression. Parity 3 and 4 + cows showed odd ratios (OR’s) for lameness of 3.92 and 8.60 respectively (95% confidence interval (CI) 2.46–6.24; 5.68–13.0). Maximum loss of Body condition score (BCS) after calving exhibits OR’s for lameness of 1.49 (95% CI 1.08–2.04) if cows lost 0.5 in BCS after calving and 2.26 (95% CI 1.30–3.95) for cows losing more than 0.5 BCS. Animals calving in BCS 3.25 and ≥ 3.5 had correlating OR’s of 0.54 (95% CI 0.34–0.87) and 0.33 (95% CI 0.16–0.65) for being lame compared to cows calving with BCS ≤ 2.75. Data gathered as part of herd health monitoring can be used in conjunction with lameness records to identify shortcomings in lameness management. Findings and recommendations on lameness management can be formulated from readily available information on cow-based risk factors for lameness.

Genetic parameters for variability in the birth of persistently infected cattle following likely in utero exposure to bovine viral diarrhea virus

Genetic selection is an inexpensive and complementary strategy to traditional methods of improving animal health and welfare. Nonetheless, endeavors to incorporate animal health and welfare traits in international breeding programs have been hampered by the availability of informative phenotypes. The recent eradication program for bovine viral diarrhea (BVD) in the Republic of Ireland has provided an opportunity to quantify the potential benefits that genetic selection could offer BVD eradication programs elsewhere, as well as inform possible eradication programs for other diseases in the Republic of Ireland. Using a dataset of 188,085 Irish calves, the estimated direct and maternal heritability estimates for the birth of persistently infected calves following likely in utero exposure to BVD virus ranged from not different from zero (linear model) to 0.29 (SE = 0.075; threshold model) and from essentially zero (linear model) to 0.04 (SE = 0.033; threshold model), respectively. The corresponding genetic SD for the direct and maternal effect of the binary trait (0, 1) ranged from 0.005 (linear model) to 0.56 (threshold model) units and ranged from 0.00008 (linear model) to 0.20 (threshold model) units, respectively. The coefficient of direct genetic variation based on the linear model was 2.56% indicating considerable genetic variation could be exploited. Based on results from the linear model in the present study, there is the potential to reduce the incidence of persistent infection in cattle by on average 0.11 percentage units per year which is cumulative and permanent. Therefore, genetic selection can contribute to reducing the incidence of persistent infection in cattle. Moreover, where populations are free from persistent infection, inclusion of the estimated genetic merit for BVD in national breeding indexes could contribute to a preservation of a BVD-free status. Results from the present study can be used to inform breeding programs of the potential genetic gains achievable. Moreover, the approaches used in the present study can be applied to other diseases when data become available.

A review of paratuberculosis in dairy herds - Part 1: Epidemiology

Bovine paratuberculosis is a chronic infectious disease of cattle caused by Mycobacterium avium subspecies paratuberculosis (MAP). This is the first in a two-part review of the epidemiology and control of paratuberculosis in dairy herds. Paratuberculosis was originally described in 1895 and is now considered endemic among farmed cattle worldwide. MAP has been isolated from a wide range of non-ruminant wildlife as well as humans and non-human primates. In dairy herds, MAP is assumed to be introduced predominantly through the purchase of infected stock with additional factors modulating the risk of persistence or fade-out once an infected animal is introduced. Faecal shedding may vary widely between individuals and recent modelling work has shed some light on the role of super-shedding animals in the transmission of MAP within herds. Recent experimental work has revisited many of the assumptions around age susceptibility, faecal shedding in calves and calf-to-calf transmission. Further efforts to elucidate the relative contributions of different transmission routes to the dissemination of infection in endemic herds will aid in the prioritisation of efforts for control on farm.

Impact of Genetic Selection for Increased Cattle Resistance to Bovine Tuberculosis on Disease Transmission Dynamics

Bovine tuberculosis (bTB) poses a challenge to animal health and welfare worldwide. Presence of genetic variation in host resistance to Mycobacterium bovis infection makes the trait amenable to improvement with genetic selection. Genetic evaluations for resistance to infection in dairy cattle are currently available in the United Kingdom (UK), enabling genetic selection of more resistant animals. However, the extent to which genetic selection could contribute to bTB eradication is unknown. The objective of this study was to quantify the impact of genetic selection for bTB resistance on cattle-to-cattle disease transmission dynamics and prevalence by developing a stochastic genetic epidemiological model. The model was used to implement genetic selection in a simulated cattle population. The model considered various levels of selection intensity over 20 generations assuming genetic heterogeneity in host resistance to infection. Our model attempted to represent the dairy cattle population structure and current bTB control strategies in the UK, and was informed by genetic and epidemiological parameters inferred from data collected from UK bTB infected dairy herds. The risk of a bTB breakdown was modeled as the percentage of herds where initially infected cows (index cases) generated secondary cases by infecting herd-mates. The model predicted that this risk would be reduced by half after 4, 6, 9, and 15 generations for selection intensities corresponding to genetic selection of the 10, 25, 50, and 70% most resistant sires, respectively. In herds undergoing bTB breakdowns, genetic selection reduced the severity of breakdowns over generations by reducing both the percentage of secondary cases and the duration over which new secondary cases were detected. Selection of the 10, 25, 50, and 70% most resistant sires reduced the percentage of secondary cases to <1% in 4, 5, 7, and 11 generations, respectively. Similarly, the proportion of long breakdowns (breakdowns in which secondary cases were detected for more than 365 days) was reduced by half in 2, 2, 3, and 4 generations, respectively. Collectively, results suggest that genetic selection could be a viable tool that can complement existing management and surveillance methods to control and ultimately eradicate bTB.

Genome-wide pleiotropy and shared biological pathways for resistance to bovine pathogens

Host genetic architecture is a major factor in resistance to pathogens and parasites. The collection and analysis of sufficient data on both disease resistance and host genetics has, however, been a major obstacle to dissection the genetics of resistance to single or multiple pathogens. A severe challenge in the estimation of heritabilities and genetic correlations from pedigree-based studies has been the confounding effects of the common environment shared among relatives which are difficult to model in pedigree analyses, especially for health traits with low incidence rates. To circumvent this problem we used genome-wide single-nucleotide polymorphism data and implemented the Genomic-Restricted Maximum Likelihood (G-REML) method to estimate the heritabilities and genetic correlations for resistance to 23 different infectious pathogens in calves and cows in populations undergoing natural pathogen challenge. Furthermore, we conducted gene-based analysis and generalized gene-set analysis to understand the biological background of resistance to infectious diseases. The results showed relatively higher heritabilities of resistance in calves than in cows and significant pleiotropy (both positive and negative) among some calf and cow resistance traits. We also found significant pleiotropy between resistance and performance in both calves and cows. Finally, we confirmed the role of the B-lymphocyte pathway as one of the most important biological pathways associated with resistance to all pathogens. These results both illustrate the potential power of these approaches to illuminate the genetics of pathogen resistance in cattle and provide foundational information for future genomic selection aimed at improving the overall production fitness of cattle.

Genetic correlations between endo-parasite phenotypes and economically important traits in dairy and beef cattle

Parasitic diseases have economic consequences in cattle production systems. Although breeding for parasite resistance can complement current control practices to reduce the prevalence globally, there is little knowledge of the implications of such a strategy on other performance traits. Records on individual animal antibody responses to Fasciola hepatica, Ostertagia ostertagi, and Neospora caninum were available from cows in 68 dairy herds (study herds); national abattoir data on F. hepatica-damaged livers were also available from dairy and beef cattle. After data edits, 9,271 dairy cows remained in the study herd dataset, whereas 19,542 dairy cows and 68,048 young dairy and beef animals had a record for the presence or absence of F. hepatica-damaged liver in the national dataset. Milk, reproductive, and carcass phenotypes were also available for a proportion of these animals as well as their contemporaries. Linear mixed models were used to estimate variance components of antibody responses to the three parasites; covariance components were estimated between the parasite phenotypes and economically important traits. Heritability of antibody responses to the different parasites, when treated as a continuous trait, ranged from 0.07 (O. ostertagi) to 0.13 (F. hepatica), whereas the coefficient of genetic variation ranged from 4% (O. ostertagi) to 20% (F. hepatica). The antibody response to N. caninum was genetically correlated with the antibody response to both F. hepatica (-0.29) and O. ostertagi (-0.67); a moderately positive genetic correlation existed between the antibody response to F. hepatica and O. ostertagi (0.66). Genetic correlations between the parasite phenotypes and the milk production traits were all close to zero (-0.14 to 0.10), as were the genetic correlations between F. hepatica-damaged livers and the carcass traits of carcass weight, conformation, and fat score evaluated in cows and young animals (0.00 to 0.16). The genetic correlation between F. hepatica-damaged livers in cows and milk somatic cell score was 0.32 (SE = 0.20). Antibody responses to F. hepatica and O. ostertagi had favorable genetic correlations with fertility traits, but conversely, antibody response to N. caninum and F. hepatica-damaged livers were unfavorably genetically correlated with fertility. This study provides the necessary information to undertake national multitrait genetic evaluations for parasite phenotypes.

Evaluation of metrics for benchmarking antimicrobial use in the UK dairy industry

The issue of antimicrobial resistance is of global concern across human and animal health. In 2016, the UK government committed to new targets for reducing antimicrobial use (AMU) in livestock. Although a number of metrics for quantifying AMU are defined in the literature, all give slightly different interpretations. This paper evaluates a selection of metrics for AMU in the dairy industry: total mg, total mg/kg, daily dose and daily course metrics. Although the focus is on their application to the dairy industry, the metrics and issues discussed are relevant across livestock sectors. In order to be used widely, a metric should be understandable and relevant to the veterinarians and farmers who are prescribing and using antimicrobials. This means that clear methods, assumptions (and possible biases), standardised values and exceptions should be published for all metrics. Particularly relevant are assumptions around the number and weight of cattle at risk of treatment and definitions of dose rates and course lengths; incorrect assumptions can mean metrics over-represent or under-represent AMU. The authors recommend that the UK dairy industry work towards the UK-specific metrics using the UK-specific medicine dose and course regimens as well as cattle weights in order to monitor trends nationally.