Digital Cushion Fatty Acid Composition and Lipid Metabolism Gene Network Expression in Holstein Dairy Cows Fed a High-Energy Diet

Associations of cow and farm characteristics with cow-level lameness using data from an extensive cross-sectional study across 3 structurally different dairy regions in Germany

The aim of the present study was to evaluate the associations between milk recording data, body condition score (BCS), housing factors, management factors, and lameness in freestall-housed dairy cows in 3 structurally different regions in Germany. These regions substantially vary regarding herd size, breeds, access to pasture, farm management (family run or company owned), and percentage of organic farms. The data used was collected in a large cross-sectional study from 2016 to 2019. A total of 58,144 cows from 651 farms in 3 regions of Germany (North, East, and South) was scored for locomotion and body condition. Additionally, data on milk yield, milk composition, breed, age, as well as information on housing and management were retrieved. One mixed-logistic regression model was fitted per region to evaluate the association of the data with the target variable "lame" and to allow for a comprehensive reflection across different kinds of farming types. In all regions, undercondition (BCS lower than recommended for the lactation stage; North: odds ratio [OR] 2.15, CI 1.96-2.34; East: OR 2.66, CI 2.45-2.88; South: OR 2.45, CI 2.01-2.98) and mid-lactation stage (102-204 d in milk; North: OR 1.15, CI 1.05-1.27; East: OR 1.24, CI 1.17-1.32; South: OR 1.38, CI 1.18-1.62) were associated with higher odds for lameness, whereas overcondition (BCS higher than recommended for the lactation stage; North: OR 0.51, CI 0.44-0.60; East: OR 0.51, CI 0.48-0.54; South: OR 0.65, CI 0.54-0.77) and parity of 1 or 2 was associated with lower odds (parity 1 = North: OR 0.32, CI 0.29-0.35; East: OR 0.19, CI 0.18-0.20; South: OR 0.28, CI 0.24-0.33; parity 2 = North: OR 0.51, CI 0.47-0.46; East: OR 0.41, CI 0.39-0.44; South: OR 0.49, CI 0.42-0.57), irrespective of the regional production characteristics. Low energy-corrected milk yield was associated with higher odds for lameness in South and North (North: OR 1.16, CI 1.05-1.27; South: OR 1.43, CI 1.22-1.69). Further factors such as pasture access for cows (North: OR 0.64, CI 0.50-0.82; and South: OR 0.65, CI 0.47-0.88), milk protein content (high milk protein content = North: OR 1.34, CI 1.18-1.52; East: OR 1.17, CI 1.08-1.28; low milk protein content = North: OR 0.79, CI 0.71-0.88; East: OR 0.84, CI 0.79-0.90), and breed (lower odds for "other" [other breeds than German Simmental and German Holstein] in East [OR 0.47, CI 0.42-0.53] and lower odds both for German Holstein and "other" in South [German Holstein: OR 0.62, CI 0.43-0.90; other: OR 0.46, CI 0.34 - 0.62]) were associated with lameness in 2 regions, respectively. The risk of ketosis (higher odds in North: OR 1.11, CI 1.01-1.22) and somatic cell count (higher odds in East: increased (>39.9 cells × 1,000/mL): OR 1.10; CI 1.03-1.17; high (>198.5 cells × 1,000/mL): OR 1.08; CI 1.01-1.06) altered the odds for lameness in 1 region, respectively. Cows from organic farms had lower odds for lameness in all 3 regions (North: OR 0.18, CI 0.11-0.32; East: OR 0.39, CI 0.28-0.56; South: OR 0.45, CI 0.29-0.68). As the dairy production systems differed substantially between the different regions, the results of this study can be viewed as representative for a wide variety of loose-housed dairy systems in Europe and North America. The consistent association between low BCS and lameness in all regions aligns with the previous literature. Our study also suggests that risk factors for lameness can differ between geographically regions, potentially due to differences in which dairy production system is predominantly used and that region-specific characteristics should be taken into account in comparable future projects.

Antagonistic effects of nano-selenium on broilers hepatic injury induced by Cr(VI) poisoning in AMPK pathway

Cr (chromium, with common valence states of III and VI) is one of the common broiler feed additives. Liver injury and metabolic disorders could be caused by Cr_(VI) (hexavalent chromium) poisoning in broilers. Oxidative damage and metabolic disorders of organisms caused by heavy metals could be antagonized by nano-Se (nano-selenium). Nano-Se was chosen to study the antagonism of Cr_(VI) poisoning in broilers. AMPK (Adenosine 5^,-monophosphate-activated protein kinase) is known as a "cell energy regulator" and plays a key regulatory role in carbohydrate and lipid metabolism. AMPK pathway and ACACA/CPT1A two genes were selected to study the prevention and treatment of nano-Se on Cr_(VI) poisoning in broilers and its molecular mechanism. For this purpose, 180 1-day-old AA (Arbor Acres) broilers were selected and randomly divided into 6 groups (n = 30) for further testing. After feeding as planned for 35 days, the livers of such broilers were taken for further examination including histopathological examination, differential gene expression analysis, and further validation on both mRNA and protein levels using related techniques like RT-qPCR, western blot, and immunohistochemistry (IHC). The histopathological examination suggested that the liver cells of the Cr_(VI) poisoning group were more severely injured than the nano-Se addition group. RT-qPCR results showed that the relative expression of ACACA gene in the Cr_(VI) poisoning group was significantly increased (P < 0.05), while the CPT1A gene's expression was significantly decreased (P < 0.01). Those results were reversed in the nano-Se addition group. Western blot results were consistent with RT-qPCR and both suggested antagonism of nano-Se on Cr_(VI). Through morphological and histopathological observation, as well as the measurement of the mRNA and protein expression levels of ACACA and CPT1A genes in AMPK pathway, it was confirmed that nano-Se has certain preventive and protective effects on Cr_(VI) poisoning in broiler chickens. Furthermore, the adverse effects of Cr_(VI) on carbohydrate and lipid metabolism in broilers can be antagonized by nano-Se through AMPK pathway. A new method and experimental basis were provided to the future study of Cr_(VI) poisoning in broilers.

Exploring physiological and genetic variation of digital cushion thickness in Holstein and Jersey cows and bulls

The objectives of our study were to expand phenotypic characterization of digital cushion thickness (DCT) to the Jersey breed and include mature bulls and to identify breed-specific quantitative trait loci (QTL) for DCT within the Jersey or Holstein breeds and common QTL across breeds and sexes to better understand biological regulation and guide future use of marker-assisted selection. In a cohort of 698 cows and 85 bulls (Holstein and Jersey) from 8 farms in New York State, data were collected on DCT, body condition score (BCS), withers height, and sacral height. All animals underwent sonographic examination of the digital cushion evaluated at the sole ulcer site for the right front and hind feet. Linear mixed models were conducted on DCT separately for cows and bulls with fixed effects of time point, breed, age, digit, and BCS group. The models included random effects to control for the random subset of animals per farm, repeated measures, and multiple measurements from each animal. The phenotypic results indicated that DCT varied by time point, breed, age, digit, and BCS group for cows and by breed, age, digit, and BCS group for bulls. For the genotypic study, 616 cow DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA), whereas 76 bull DNA samples were genotyped on different platforms ranging from 5K to 150K. Multiple genome-wide association studies were conducted to highlight pertinent phenotyping parameters and genetic markers for genomic selection. Data were separated into 8 data sets based on different combinations of breed and sex. Each data set was assessed for quality of markers and samples before conducting genome-wide association studies for DCT, testing the inheritance models and genetic variation of digit, foot, and average thickness. Ten markers passed the Bonferroni correction threshold and 9 passed false discovery rate from 10 genome-wide association studies using a combination of the covariates breed, sex, genotyping batch plate, age, BCS, withers height, and sacral height. Of the 43 candidate genes, 8 novel biologically plausible genes were identified on Bos taurus autosomes 3, 4, 7, and 9: SFRS18 and LRRFIP1 function in fat deposition, whereas AHR, BZW2, EFNA5, USP45, and VAV3 effect bone growth, and SOSTDC1 is related to epidermal keratinocyte function. The genetic markers associated with DCT in this study were explored for variation between cows and bulls within and across breeds for their potential use in marker-assisted selection.

Phenotypic characterization and genome-wide association studies of digital cushion thickness in Holstein cows

The bovine digital cushion is a compression pad between the distal phalanx and sole and has been associated with claw horn disruption lesions. Digital cushion thickness (DCT) is estimated to be moderately heritable. Therefore, the objectives of our study were to examine influences of management and environment on DCT and to identify genetic markers and candidate genes associated with DCT. In a cohort of 502 Holsteins from 5 farms in New York State, DCT and body condition score (BCS) were collected twice, at <137 d prepartum and from 86 to 127 d in milk, corresponding to periods when the digital cushion is thickest and thinnest, respectively, as determined by previous research. Cows underwent sonographic examination of the digital cushion evaluated at the typical sole ulcer site for the right front and hind foot. Linear mixed models were conducted on DCT with the fixed effects of time point, digit, wither height, sacral height, BCS group, and multiple farm system variables separately and included random effects to control for the random subset of cows per farm, repeated measures, and multiple measurements from each cow. The phenotypic results indicated that DCT varied by sample time point, sacral height, parity, digit, BCS group, and wither height. For the genotypic study, 447 DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA). Quality assessment of markers and samples provided a final data set of 431 samples and 579,449 markers. Genome-wide association studies were conducted for DCT testing inheritance models and genetic variation of digit, foot, time point, and average thickness. One marker passed the Bonferroni correction threshold and 26 passed false discovery rate from 4 genome-wide association studies with covariates of sequencing batch plate, parity group, BCS, wither height, and sacral height. Ten candidate genes were identified, with 2 genes on Bos taurus autosomes 24 and 29 involved in biological functions related to the digital cushion: MC4R and DLG2 were related to fat deposition and bone growth, respectively. The genetic markers discovered in this study have the opportunity to be used in breeding programs using genomic selection to select against claw horn disruption lesions and lameness due to associations between the markers and DCT. Further studies on the biologically plausible candidate genes may identify causative genetic variants and how they relate to DCT through gene regulation, expression, structure, or copy number variation.

Factors associated with lameness prevalence in lactating cows housed in freestall and compost-bedded pack dairy farms in southern Brazil

The aim of this cross-sectional study was to investigate factors associated with lameness in dairy cows on intensive farms in southern Brazil. Farms (freestall: n = 38; compost-bedded pack: n = 12) having on average 274 (range: 41-901) lactating cows were visited once in 2016 (March to October). Potential risk factors for lameness at the cow, pen and herd levels were investigated through inspection of facilities, examination of cows and the use of data collected on routine management practices. All milking cows on each farm were assessed for gait score and BCS (n = 13,716). Associations between lameness, days in milk (DIM), BCS, parity, and test-day milk yield were investigated in 16 farms with available data (n = 5,301 cows). Mixed-effects linear and logistic regressions were used to model the data. Within-herd lameness prevalence was 41.1% (range: 13.8-64.5, SD = 11.3). First- and second-lactation cows after 120 DIM and older cows after 335 DIM were more likely to be observed lame than early lactation cows. Greater parity and low BCS (≤ 2.75) were associated with increased odds of lameness. Severely lame cows had lower milk yield (on average 1.3 kg/d) than non-lame cows. Freestall dairies using mattresses as stall base had greater within-pen (95% CI: 52-69%) and herd-level (38-57%) lameness prevalence than compost-bedded farms. Higher lameness prevalence was observed on farms having slippery feed bunk alley floors (32-49%) and shortened dry periods (< 60 days: 32-42%). First-lactation pens had lower predicted within-pen lameness prevalence (0-4%) and special-needs pens higher (52-73%) compared with the prevalence observed in compost-bedded fresh-cow pens (19-41%). Freestall pens using sawdust deep-bedding had greater (46-68%) within-pen lameness prevalence, but the prevalence in barns using sand was not different from compost-bedded farms. Given the high lameness prevalence observed in this study, there is a great opportunity for implementation of lameness prevention programs. Providing walking surfaces with high traction that facilitates mobility and using soft and deep-bedded material, such as compost and sand (and avoiding the use of mattresses) may reduce lameness prevalence in the types of dairy farms visited in this study.

Lipidomic Profiling Reveals Disruption of Lipid Metabolism in Valproic Acid-Induced Hepatotoxicity

Valproic acid (VPA) is one of the most widely prescribed antiepileptic drugs, as VPA-induced hepatotoxicity is one of the most severe adverse reaction that can lead to death. The objective of this study was to gain an understanding of dysregulated lipid metabolism in mechanism of hepatotoxicity. Nontargeted lipidomics analysis with liquid chromatography-quadrupole-time-of-flight mass spectrometry (LC-Q-TOF/MS) was performed to explore differential lipids from the patient serum and L02 cells. Lipidomics data interpretation was augmented by gene expression analyses for the key enzymes in lipid metabolism pathways. From patient serum lipidomics, pronouncedly changed lipid species between abnormal liver function (ALF) patients and normal liver function (NLF) patients were identified. Among these lipid species, LPCs, Cers, and SMs were markedly reduced in the ALF group and showed negative relationships with liver injury severity [alanine aminotransferase (ALT) levels], while significantly increased triacylglycerols (TAG) with higher summed carbon numbers demonstrated a positive relationship with ALT levels. Regarding lipidomics in hepatic L02 cells, TAG was markedly elevated after VPA exposure, especially in TAGs with more than 53 summed carbons. Besides, gene expression analysis revealed dysregulated lipid metabolism in VPA-treated L02 cells. Peroxime proliferators-activated receptor (PPARγ) pathway played an important role in VPA-induced lipid disruption through inducing long-chain fatty acid uptake and TAG synthesis, which was also regulated by Akt pathway. Our findings present that VPA-induced lipid metabolism disruption might lead to lipotoxicity in the liver. This approach is expected to be applicable for other drug-induced toxicity assessments.