Serum hepatokines in dairy cows: periparturient variation and changes in energy-related metabolic disorders

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Hepatic Lipidosis in Ruminants

Hepatic lipidosis (ie, fatty liver) occurs primarily during the first weeks of lactation in dairy cows because of excessive lipolysis overwhelming the concomitant capacity for beta-oxidation and hepatic export of triglycerides. Besides economic losses due to reduced lactational and reproductive performance, close associations with concomitantly occurring infectious and metabolic health disorders, in particular ketosis, exist. Hepatic lipidosis is not only a consequence from the postpartal negative energy balance but also acts as a disease component for further health disorders.

Performance and Metabolic, Inflammatory, and Oxidative Stress-Related Parameters in Early Lactating Dairy Cows with High and Low Hepatic FGF21 Expression

Induction of FGF21 expression in the liver and a significant increase in plasma FGF21 concentration have been demonstrated in cows during early lactation, but knowledge about the function of FGF21 in dairy cows remains limited. In order to improve the understanding of the physiological role of FGF21 in dairy cows, the present study aimed to investigate differences in metabolic pathways between dairy cows with high and low hepatic expression of FGF21 at week 1 of lactation (n = 8/group) by liver transcriptomics, targeted plasma metabolomics, and analysis of inflammatory and oxidative stress-related parameters. Dry matter intake, energy balance, milk yield, and energy-corrected milk yield at days 8−14 postpartum did not differ between cows with high and low hepatic FGF21 expression. However, cows with high FGF21 expression showed an upregulation of genes involved in endoplasmic reticulum stress, inflammation, and nuclear factor E2-related factor 2 (Nrf2)-dependent cytoprotection compared to cows with low FGF21 expression at week 1 postpartum (p < 0.05). Concentrations of important antioxidants (tocopherols, β-carotene, and glutathione) in the liver and plasma, trolox equivalent antioxidant capacity in plasma, concentrations of oxidative stress-related compounds (thiobarbituric acid-reactive substances and protein carbonyls), and levels of most acute phase proteins at week 1 postpartum did not differ between cows with high or low FGF21 expression. Moreover, among a total of >200 metabolites assayed in the plasma, concentrations of only 7 metabolites were different between cows with high or low FGF21 expression (p < 0.05). Overall, the results showed that cows with high and low FGF21 hepatic expression had only moderate differences in metabolism, but FGF21 might be important in the adaptation of dairy cows to stress conditions during early lactation.

FGF21 Reduces Lipid Accumulation in Bovine Hepatocytes by Enhancing Lipid Oxidation and Reducing Lipogenesis via AMPK Signaling

During the periparturient period, dairy cows suffer drastic metabolic stress because of plasma increased non-esterified fatty acids (NEFAs) that stem from a negative energy balance. Fibroblast growth factor 21 (FGF21) is a hepatokine that activates the AMP-activated protein kinase (AMPK) signaling pathway to maintain intracellular energy balance and tissue integrity via the promotion of catabolism and the inhibition of anabolic regulation. FGF21 treatment caused a 50% reduction in triglyceride (TG) content in liver in dairy cows. However, it is not clear whether FGF21 regulates lipid metabolism in bovine liver. The purpose of this study was to evaluate the influence of FGF21 on lipid metabolism via AMPK signaling in bovine hepatocytes. The hepatocytes isolated from calves were treated with different concentrations of FGF21 or co-treated with AMPK inhibitor (BML-275). Herein, the study showed that FGF21 significantly reduced TG content in a dose–response manner and promoted very-low-density lipoprotein (VLDL) secretion via an up-regulation of the proteins (ApoB 100, ApoE and MTTP) involved in VLDL secretion. Otherwise, the genes associated with lipid transport (LDLR and CD36) and lipid oxidation (PPARGC1A, ACOX1 and CPT1A), were up-regulated following FGF21 treatment. Moreover, FGF21 treatment inhibited lipogenesis via SREBF1, ACACA, FASN and ACLY inhibition. After being co-treated with the AMPK inhibitor, FGF21-induced changes were reversed in some genes. In conclusion, these results indicate that FGF21 adaptively regulates energy metabolism for a negative impact on lipogenesis, strengthens lipid oxidation, and inhibited lipid transportation via AMPK signaling in bovine hepatocytes. The present data suggest the possibility that FGF21 has potential value in alleviating perinatal metabolic diseases in dairy cows, and specific research in vivo should be studied in more detail.

Metabogenomic analysis to functionally annotate the regulatory role of long non-coding RNAs in the liver of cows with different nutrient partitioning phenotype

Long non-coding RNAs (lncRNAs) hold gene regulatory potential, but require substantial further functional annotation in livestock. Applying two metabogenomic approaches by combining transcriptomic and metabolomic analyses, we aimed to identify lncRNAs with potential regulatory function for divergent nutrient partitioning of lactating crossbred cows and to establish metabogenomic interaction networks comprising metabolites, genes and lncRNAs. Through correlation analysis of lncRNA expression with transcriptomic and metabolomic data, we unraveled lncRNAs that have a putative regulatory role in energy and lipid metabolism, the urea and tricarboxylic acid cycles, and gluconeogenesis. Especially FGF21, which correlated with a plentitude of differentially expressed genes, differentially abundant metabolites, as well as lncRNAs, suggested itself as a key metabolic regulator. Notably, lncRNAs in close physical proximity to coding-genes as well as lncRNAs with natural antisense transcripts appear to perform a fine-tuning function in gene expression involved in metabolic pathways associated with different nutrient partitioning phenotypes.

The crucial role of lysine in the hepatic metabolism of growing Holstein dairy heifers as revealed by LC-MS-based untargeted metabolomics

The objective of this experiment was to evaluate the effect of supplementing rumen-protected Lys based on a Lys-deficient diet on liver metabolism in growing Holstein heifers. The experiment was conducted for 3 months with 36 Holstein heifers (initial body weight: 200 ± 9.0 kg; 7-month-old). Heifers were randomly assigned to 2 diets based on corn, soybean meal, alfalfa hay, and wheat bran: control, Lys-deficient diet (LD; 0.66% Lys in diet), and Lys-adequate diet (LA; 1.00% Lys in diet). The results showed no difference in growth performance between the 2 groups (P > 0.05). However, there was a clear trend of increasing feed conversion rate with Lys supplementation (0.05 < P < 0.01). The serum urea nitrogen concentration was significantly decreased, and the aspartate aminotransferase-to-alanine aminotransferase ratio was significantly decreased by Lys supplementation (P < 0.05). Moreover, growing heifers fed a Lys-adequate diet had lower levels of urine nitrogen excretion and higher levels of the biological value of nitrogen (P < 0.05). Metabolomic analysis revealed that 5 types of phosphatidylcholine and 3 types of ceramide were significantly increased and enriched in sphingolipid metabolism and glycerophospholipid metabolism (P < 0.05). His, Leu, and Asp levels were significantly decreased in the liver following Lys supplementation (P < 0.05). In conclusion, Lys supplementation may promote the synthesis of body tissue proteins, as evidenced by significantly decreased amino acids in the liver and urine N excretion, it also improves hepatic lipid metabolism by providing lipoprotein precursors.

NEFA Promotes Autophagosome Formation through Modulating PERK Signaling Pathway in Bovine Hepatocytes

During the perinatal period, the abnormally high plasma non-esterified fatty acids (NEFA) concentration caused by the negative energy balance (NEB) can impose a significant metabolic stress on the liver of dairy cows. Endoplasmic reticulum (ER) stress is an important adaptive response that can serve to maintain cell homeostasis in the event of stress. The protein kinase R-like endoplasmic reticulum kinase (PERK) pathway is the most rapidly activated cascade when ER stress occurs in cells and has an important impact on the regulation of hepatic lipid metabolism and autophagy modulation. However, it is unknown whether NEFA can affect autophagy through modulating the PERK pathway, under NEB conditions. In this study, we provide evidence that NEFA treatment markedly increased lipid accumulation, the phosphorylation level of PERK and eukaryotic initiation factor 2α (eIF2α), and the expression of glucose-regulated protein 78 (Grp78), activating transcription factor 4 (ATF4), and C/EBP homologous protein (CHOP). More importantly, NEFA treatment can cause a substantial increase in the protein levels of autophagy-related gene 7 (ATG7), Beclin-1 (BECN1), sequestosome-1 (p62), and microtubule-associated protein 1 light chain 3 (LC3)-II, and in the number of autophagosomes in primary bovine hepatocytes. The addition of GSK2656157 (PERK phosphorylation inhibitor) can significantly inhibit the effect of NEFA on autophagy and can further increase lipid accumulation. Overall, our results indicate that NEFA could promote autophagy via the PERK pathway in bovine hepatocytes. These findings provide novel evidence about the potential role of the PERK signaling pathway in maintaining bovine hepatocyte homeostasis.

Fibroblast growth factor 21 in dairy cows: current knowledge and potential relevance

Fibroblast growth factor 21 (FGF21) has been identified as an important regulator of carbohydrate and lipid metabolism, which plays an important role for metabolic regulation, particularly under conditions of energy deprivation or stress conditions. Dairy cows are subjected to a negative energy balance and various kinds of stress particularly during the periparturient phase and during early lactation. It has been shown that the plasma concentration of FGF21 in dairy cows is dramatically increased at parturition and remains high during the first weeks of lactation. This finding suggests that FGF21 might exert similar functions in dairy cows than in other species, such as mice or humans. However, the role of FGF21 in dairy cows has been less investigated so far. Following a brief summary of the previous findings about the function of FGF21 in humans and mice, the present review aims to present the current state of knowledge about the role of FGF21 in dairy cows. The first part of the review deals with the tissue localization of FGF21 and with conditions leading to an upregulation of FGF21 expression in the liver of dairy cows. In the second part, the influence of nutrition on FGF21 expression and the role of FGF21 for metabolic diseases in dairy cows is addressed. In the third part, findings of exogenous FGF21 application on metabolism in dairy cows are reported. Finally, the potential relevance of FGF21 in dairy cows is discussed. It is concluded that FGF21 might be of great importance for metabolic adaptation to negative energy balance and stress conditions in dairy cows. However, further studies are needed for a better understanding of the functions of FGF21 in dairy cows.

Oxidative status in dairy goats: periparturient variation and changes in subclinical hyperketonemia and hypocalcemia

BACKGROUND

A better comprehension of the redox status during the periparturient period may facilitate the development of management and nutritional solutions to prevent subclinical hyperketonemia (SCHK) and subclinical hypocalcemia (SCHC) in dairy goats. We aimed to evaluate the variation in the redox status of dairy goats with SCHK and SCHC during their periparturient periods. Guanzhong dairy goats (n = 30) were assigned to SCHK (n = 10), SCHC (n = 10), and healthy (HEAL, n = 10) groups based on their blood β-hydroxybutyrate (BHBA) and calcium (Ca) concentrations. Blood were withdrawn from goats every week from 3 weeks before the expected parturition date to 3 weeks post-kidding. On the same day, the body condition scores (BCS) were evaluated, and the milk yield was recorded for each goat. The metabolic profile parameters and the indicators of oxidative status were determined by using the standard biochemical techniques.

RESULTS

In comparison with the HEAL goats, SCHK and SCHC goats presented with a more dramatic decline of BCS post-kidding and a significant decrease in the milk yield at 2- and 3-weeks postpartum, ignoring the obvious increase at 1-week postpartum. The levels of non-esterified fatty acids (NEFA) peaked at parturition, exhibiting significantly higher levels from 1-week prepartum to the parturition day in the SCHK and SCHC groups. The malondialdehyde (MDA) concentration was increased in the SCHK goats from 1-week antepartum until 3-weeks postpartum, with its concentration being significantly higher in the SCHC goats at parturition. The hydrogen peroxide (H2O2) concentration was significantly lower in the SCHK and SCHC goats from 2-weeks antepartum to 1-week post-kidding. The total antioxidant capacity (T-AOC) and the superoxide dismutase (SOD) level were decreased at 1-week antepartum in the SCHK and SCHC goats, respectively. The glutathione peroxidase (GSH-Px) level was increased in the SCHK and SCHC goats during the early lactation period.

CONCLUSIONS

The SCHK and SCHC goats exerted more efforts to maintain their redox homeostasis and to ensure the production performance than the HEAL goats during their periparturient period, probably owing to more intense fat mobilization and lipid peroxidation in the former.

Exploring the liver fibrosis induced by deltamethrin exposure in quails and elucidating the protective mechanism of resveratrol

Deltamethrin (DLM) is widely used in agriculture and the prevention of human insect-borne diseases. However, the molecular mechanism of DLM induced liver injury remains unclear to date. This study investigated the potential molecular mechanism that DLM induced liver fibrosis in quails. Japanese quails received resveratrol (500 mg/kg) daily with or without DLM (45 mg/kg) exposure for 12 weeks. Histopathology, transmission electron microscopy, biochemical indexes, TUNEL, quantitative real-time PCR, and western blot analysis were performed. DLM exposure induced hepatic steatosis, oxidative stress, inflammation, and apoptosis. Most importantly, the Nrf2/TGF-β1/Smad3 signaling pathway played an important role on DLM-induced liver fibrosis in quails. Interestingly, the addition of resveratrol, an Nrf2 activator, alleviates oxidative stress and inflammation response by activating Nrf2, thereby inhibits the liver fibrosis induced by DLM in quails. Collectively, these findings demonstrate that chronic exposure to DLM induces oxidative stress via the Nrf2 expression inhibition and apoptosis, and then results in liver fibrosis in quails by the activation of NF-κB/TNF-α and TGF-β1/Smad3 signaling pathway.

Metabolomic and Proteomic Profiles Associated With Ketosis in Dairy Cows

Ketosis is a common metabolic disease in dairy cows during early lactation. However, information about the metabolomic and proteomic profiles associated with the incidence and progression of ketosis is still limited. In this study, an integrated metabolomics and proteomics approach was performed on blood serum sampled from cows diagnosed with clinical ketosis (case, ≥ 2.60 mmol/L plasma β-hydroxybutyrate; BHBA) and healthy controls (control, < 1.0 mmol/L BHBA). Samples were taken 2 weeks before parturition and 2 weeks after parturition from 19 animals (nine cases, 10 controls). All serum samples (n = 38) were subjected to Liquid Chromatography-Mass Spectrometry (LC-MS) based metabolomic analysis, and 20 samples underwent Data-Independent Acquisition (DIA) LC-MS based proteomic analysis. A total of 97 metabolites and 540 proteins were successfully identified, and multivariate analysis revealed significant differences in both metabolomic and proteomic profiles between cases and controls. We investigated clinical ketosis-associated metabolomic and proteomic changes using statistical analyses. Correlation analysis of statistically significant metabolites and proteins showed 78 strong correlations (correlation coefficient, R ≥ 0.7) between 38 metabolites and 25 proteins, which were then mapped to pathways using IMPaLA. Results showed that ketosis altered a wide range of metabolic pathways, such as metabolism, metabolism of proteins, gene expression and post-translational protein modification, vitamin metabolism, signaling, and disease related pathways. Findings presented here are relevant for identifying molecular targets for ketosis and biomarkers for ketosis detection during the transition period.

Angiopoietin-like protein 4 promotes very-low-density lipoprotein assembly and secretion in bovine hepatocytes in vitro

In dairy cows, fatty liver is one of the most common metabolic diseases that occurs during the periparturient period. Angiopoietin-like protein 4 (ANGPTL4) is a well-known downstream target of peroxisome proliferator-activated receptors (PPARs), which regulate the glucose and fatty acid metabolisms. The inhibition of lipoprotein lipase (LPL) activity interferes with the storage of triglycerides (TG) in adipocytes, which plays an essential role in lipid metabolism in rodents. However, it remains unclear whether ANGPTL4 is involved in the pathological process of fatty liver in dairy cows as a result of the regulation of the hepatocellular lipid transport system. This study intended to investigate the effect of ANGPTL4 on the very-low-density lipoprotein (VLDL) assembly and secretion in bovine hepatocytes. Bovine hepatocytes were isolated using a modified two-step perfusion and collagenase digestion process, and treated with different concentrations of ANGPTL4 (0, 4, 12, and 24 ng/ml) for 24 hr. The results showed that a high concentration of ANGPTL4 could significantly increase the extracellular concentration of VLDL while reducing the intracellular content of TG. Thus, it was confirmed that ANGPTL4 could promote the transport of TG in the form of VLDL by partially regulating the expression of related proteins in hepatocytes, thereby contributing to the partial adaptive regulation of lipid transport in dairy cows.

Mitochondrial dysfunction and endoplasmic reticulum stress in calf hepatocytes are associated with fatty acid-induced ORAI calcium release-activated calcium modulator 1 signaling

The store-operated Ca²⁺ entry (SOCE) moiety ORAI calcium release-activated calcium modulator 1 (ORAI1) located in the endoplasmic reticulum (ER) participates in key cellular functions such as protein folding, transport, and secretion, and lipid metabolism. We used an in vitro approach to test whether exogenous fatty acids alter ORAI1 signaling and to explore potential consequences on mitochondrial dysfunction and ER stress. First, hepatocytes isolated from 4 healthy female calves (1 d old, 40-50 kg) were challenged with a 1.2 mM mixture of oleic, linoleic, palmitic, stearic, and palmitoleic acids for 0.5, 1, 3, 6, 9, and 12 h to measure oxidative stress [intracellular reduced glutathione (GSH), superoxide dismutase (SOD), malondialdehyde (MDA), and hydrogen peroxide] and ER stress (protein abundance of PERK, IRE, ATF6, and GRP78). Concentrations of GSH and SOD decreased at 0.5 h, and MDA and hydrogen peroxide increased at 1 h; ER stress proteins increased at 6 h. To determine whether ER stress was caused by oxidative stress, primary calf hepatocytes were treated with the same 1.2 mM fatty acid mix or the reactive oxygen species (ROS) inhibitor N-acetylcysteine (NAC) for 6 h. We found that NAC prevented an increase in ER stress protein abundance. Next, the role of ORAI1 on ER stress was measured by transfecting hepatocytes with small interfering (si)ORAI1 or the ORAI1 inhibitor BTP2, followed by a challenge with 1.2 mM fatty acids for 3 h. Without inhibiting ORAI1, exogenous fatty acids upregulated ORAI1 mRNA and protein abundance, oxidative stress, ER stress proteins, and protein abundance of marker indicators of an opened mitochondrial permeability transition pore (mPTP). Inhibition with BPT2 or silencing via siORAI1 abrogated oxidative stress, including increased GSH concentration and SOD activity, decreased MDA, hydrogen peroxide, and ROS concentration; ER stress protein abundance was downregulated, and mitochondrial function was restored. Last, changes in markers of mPTP opening were evaluated by culturing hepatocytes for 6 h with the sarcoendoplasmic Ca²⁺ ATPase inhibitor thapsigargin or the calcium ionophore ionomycin. We detected an increase in VDAC1, CLPP, and CypD protein abundance, all of which indicated opening of the mPTP. Overall, data from these in vitro studies suggest that ORAI1 mediates ER stress induced by high concentrations of fatty acids, in part through alleviating mitochondrial dysfunction caused by oxidative stress.

Clinical Ketosis-Associated Alteration of Gene Expression in Holstein Cows

Ketosis is one of the most prevalent transition metabolic disorders in dairy cows, and has been intrinsically influenced by both genetic and nutritional factors. However, altered gene expression with respective to dairy cow ketosis has not been addressed yet, especially at the genome-wide level. In this study, we recruited nine Holsteins diagnosed with clinical ketosis and ten healthy controls, for which whole blood samples were collected at both prepartum and postpartum. Four groups of blood samples were defined: from cows with ketosis at prepartum (PCK, N = 9) and postpartum (CK, N = 9), respectively, and controls at prepartum (PHC, N = 10) and postpartum (HC, N = 10). RNA-Seq approach was used for investigating gene expression, by which a total of 27,233 genes were quantified with four billion high-quality reads. Subsequently, we revealed 75 and four differentially expressed genes (DEGs) between sick and control cows at postpartum and prepartum, respectively, which indicated that sick and control cows had similar gene expression patterns at prepartum. Meanwhile, there were 95 DEGs between postpartum and prepartum for sick cows, which showed depressed changes of gene expression during this transition period in comparison with healthy cows (428 DEGs). Functional analyses revealed the associated DEGs with ketosis were mainly involved in biological stress response, ion homeostasis, AA metabolism, energy signaling, and disease related pathways. Finally, we proposed that the expression level of STX1A would be potentially used as a new biomarker because it was the only gene that was highly expressed in sick cows at both prepartum and postpartum. These results could significantly help us to understand the underlying molecular mechanisms for incidence and progression of ketosis in dairy cows.

Ketotic cows display a different serum nonesterified fatty acid composition

The experiments reported in this research communication aimed to compare the serum nonesterified fatty acid (NEFA) composition in ketotic cows and healthy cows during the perinatal period. NEFAs play significant roles in etiology and pathology of ketosis. We hypothesized that ketotic cows will display a different serum NEFA composition compared to healthy controls, and fatty acid related indicators for ketosis prediction can be screened. Pre-partum healthy cows were recruited, and blood samples were collected on -7, 3, 7, 14 and 21 d postpartum. Cows were further divided into a healthy control group (C group, n = 6) and a ketosis group (K group, n = 6) if blood β-hydroxybutyric acid levels exceeded 1.2 mm during the experiment. NEFA composition was then analyzed by means of Gas Chromatography-Mass Spectrometer (GC-MS). Only C12 : 0% was significantly higher in C group than K group on 7 d pre-partum (P < 0.05), when the cows were not diagnosed with ketosis. Five fatty acids displayed statistical differences in composition between C and K group (P < 0.05), namely C12 : 0, C16 : 0, C17 : 0, C18 : 1n9 and C22 : 1n9. Saturates%, unsaturates%, mono-unsaturates% and saturates/unsaturates were also different between C and K group (P < 0.05). Of note, C18 : 1n9/C12 : 0 and C18 : 1n9/C22 : 1n9 in K group were significantly higher than those in controls on 7 d pre-partum (P < 0.05). It is suggested that the ratios show potential as indicators for prediction of ketosis.

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.