Early lactation production, health, and welfare characteristics of cows selected for extended lactation

Association between days post-conception and lactation persistency in dairy cattle

Determining the optimal insemination moment for individual cows is complex, particularly when considering the effects of pregnancy on milk production. The effect of pregnancy on the absolute milk yield has already been reported in several studies. Currently, there is limited quantitative knowledge about the association between days post-conception (DPC) and lactation persistency, based on a lactation curve model, and, specifically, how persistency changes during pregnancy and relates to the days in milk at conception (DIMc). Understanding this association might provide valuable insights to determine the optimal insemination moment. This study, therefore, aimed to investigate the association between DPC and lactation persistency, with an additional focus on the influence of DIMc. Available milk production data from 2005 to 2022 were available for 23,908 cows from 87 herds located throughout the Netherlands and Belgium. Persistency was measured by a lactation curve characteristic decay, representing the time taken to halve milk production after peak yield. Decay was calculated for 8 DPC (0, 30, 60, 90, 120, 150, 180, and 210 d after DIMc) and served as the dependent variable. Independent variables included DPC, DIMc (≤60, 61-90, 91-120, 121-150, 151-180, 181-210, >210), parity group, DPC × parity group, DPC × DIMc, and variables from 30 d before DIMc as covariates. The results showed an increase in decay, which is to say, a decrease in persistency, during pregnancy for both parity groups, albeit in different ways. Specifically, from DPC 150 to DPC 210, multiparous cows showed a greater decline in persistency compared with primiparous cows. Furthermore, a later DIMc (cows conceiving later) was associated with higher persistency. Except for the early DIMc groups (DIMc <90), DIMc does not affect the change in persistency by gestation. The findings from this study contribute to a better understanding of how DPC and DIMc during lactation influence lactation persistency, enabling more informed decision-making by farmers who wish to take persistency into account in their reproduction management.

Customized voluntary waiting period before first insemination in primiparous dairy cows. Effect on milk production, fertility and health

Customized voluntary waiting period (VWP) before first insemination was tested in 18 commercial dairy herds in Sweden, to assess milk production, fertility and health in primiparous cows expected to be suited for extended VWP. Cow selection for extended VWP was based on 3 criteria in early lactation: 1) the 10% of cows with highest genomic persistency index, 2) cows with a difficult calving or disease during the first month of lactation, and 3) cows with higher yield during d 4-33 after calving than the herd average for primiparous cows. Cows meeting at least one of these criteria were randomly assigned to either the ExtExt treatment (extended VWP of at least 175 d, n = 174; calving interval (CInt) = 16.3 mo) or the ExtConv treatment (conventional VWP of maximum 100 d n = 173; CInt = 12.4 mo). Cows not meeting any of the criteria were assigned to the ConvConv treatment (conventional VWP, n = 183; CInt = 12.0 mo). There were no differences in milk yield per day in the CInt between treatments, although 305-d and whole-lactation (WL) milk yields were higher for ExtExt cows (10,371 and 13,803 kg) than ExtConv cows (9,812 and 10,257 kg). Milk yield at the last test milking before dry-off was lower in ExtExt compared with ExtConv cows (24.9 vs 28.3), however the results showed no difference in dry period length between the treatments. Regarding reproductive performance, the ExtExt cows had higher first service conception rate (FSCR; 60% vs. 45%) and lower number of inseminations per conception (NINS; 1.67 vs. 2.19), compared with the ExtConv cows. As expected, ConvConv cows had the lowest milk yield in 305-d, in WL, and per day in the CInt, however, FSCR and NINS did not differ between ConvConv cows and cows in the other 2 VWP treatments. Disease incidence was higher for cows in the ExtConv compared with the ConvConv treatment, but there was no difference between ExtExt and the 2 other VWP treatments. Further, no difference in proportion of cows with good udder health or culling rate was detected between any of the treatments, though due to low prevalence the study lacked power to draw major conclusions on these results. Thus prolonging VWP for suitable primiparous cows can produce benefits such as improved fertility in the form of higher FSCR and lower NINS, as well as lower dry-off yield, without compromising milk yield or prolonging dry period length.

Fitting mathematical functions to extended lactation curves and forecasting late-lactation milk yields of dairy cows

A 305-d lactation followed by a 60-d dry period has traditionally been considered economically optimal, yet dairy cows in modern intensive dairy systems are frequently dried off while still producing significant quantities of milk. Managing cows for an extended lactation has reported production, welfare, and economic benefits, but not all cows are suitable for an extended lactation. Implementation of an extended lactation strategy on-farm could benefit from use of a decision support system, based on a mathematical lactation model, that can identify suitable cows during early lactation that have a high likelihood of producing above a target milk yield (MY) at 305 d in milk (DIM). Therefore, our objectives were (1) to compare the suitability of 3 commonly used lactation models for modeling extended lactations (Dijkstra, Wood, and Wilmink) in primiparous and multiparous cows under a variety of lactation lengths, and (2) to determine the amount of early-lactation daily MY data needed to accurately forecast MY at d 305 by using the most suitable model and determine whether this is sufficient for identifying cows suitable for an extended lactation before the end of a typical voluntary waiting period (50-90 d). Daily MY data from 467 individual Holstein-Friesian lactations (DIM >305 d; 379 ± 65-d lactation length [mean ± SD]) were fitted by the 3 lactation models using a nonlinear regression procedure. The parameter estimates of these models, lactation characteristics (peak yield, time to peak yield, and persistency), and goodness-of-fit were compared between parity and different lactation lengths. The models had similar performance, and differences between parity groups were consistent with previous literature. Then, data from only the first i DIM for each individual lactation, where i was incremented by 30 d from 30 to 150 DIM and by 50 d from 150 to 300 DIM, were fitted by each model to forecast MY at d 305. The Dijkstra model was selected for further analysis, as it had superior goodness-of-fit statistics for i= 30 and 60. The data set was fit twice by the Dijkstra model, with parameter bounds either unconstrained or constrained. The quality of predictions of MY at d 305 improved with increasing data availability for both models and assisting the model fitting procedure with more biologically relevant constraints on parameters improved the predictions, but neither was reliable enough for practical use on-farm due to the high uncertainty of forecasted predictions. Using 90 d of data, the constrained model correctly classified 66% of lactations as being above or below a target MY at d 305 of 25 kg/d, with a probability threshold of 0.95. The proportion of correct classifications became smaller at lower targets of MY at d 305 and became greater when using more lactation days. Overall, further work is required to develop a model that can forecast late-lactation MY with sufficient accuracy for practical use. We envisage that a hybridized machine learning and mechanistic model that incorporates additional historical and genetic information with early-lactation MY could produce meaningful lactation curve forecasts.

Prediction of persistency for day 305 of lactation at the moment of the insemination decision

When deciding on the voluntary waiting period of an individual cow, it might be useful to have insight into the persistency for the remainder of that lactation at the moment of the insemination decision, especially for farmers who consider persistency in their reproduction management. Currently, breeding values for persistency are calculated for dairy cows but, to our knowledge, prediction models to accurately predict persistency at different moments of insemination are lacking. This study aimed to predict lactation persistency for DIM 305 at different insemination moments (DIM 50, 75, 100, and 125). Available cow and herd level data from 2005 to 2022 were collected for a total of 20,508 cows from 85 herds located in the Netherlands and Belgium. Lactation curve characteristics were estimated for every daily record using the data up to and including that day. Persistency was defined as the number of days it takes for the milk production to decrease by half during the declining stage of lactation, and calculated from the estimated lactation curve characteristic 'decay'. Four linear regression models for each of the selected insemination moment were built separately to predict decay at DIM 305 (decay-305). Independent variables included the lactation curve characteristics at the selected insemination moment, daily milk yield, age, calving season, parity group and other herd variables. The average decay-305 of primiparous cows was lower than that of multiparous cows (1.55 *10-3 vs. 2.41*10-3, equivalent to a persistency of 447 vs. 288 days, respectively). Results showed that our models had limitations in accurately predicting persistency, although predictions improved slightly at later insemination moments, with R2 values ranging between 0.27 and 0.41. It can thus be concluded that, based only on cow and herd milk production information, accurate prediction of persistency for DIM 305 is not feasible.

Graduate Student Literature Review: The challenge of drying-off high-yielding dairy cows

The cessation of lactation (i.e., dry-off) in dairy cattle is an area of research that has received much focus in recent years. The dry period is necessary to optimize tissue remodeling of the mammary gland, but represents a stressful event, incorporating several changes in daily routine, diet, and metabolism. Moreover, the high milk yields achieved by modern cows in late gestation exacerbate the need for relevant manipulations in the days around dry-off, as excessive accumulation of milk might jeopardize the success of the dry period, with potential negative effects on future lactation. Production levels over 15 kg/d are an additional risk factor for udder health, delay mammary involution, and worsen metabolic stress and inflammatory responses. Furthermore, the pressure to reduce antibiotic usage in farm animals has resulted in increased attention on the dry period, given that historically most dairy cattle were provided prophylactic intramammary antibiotic treatment at dry-off as a means to reduce the risk of intramammary infections in the subsequent lactation. Several strategies have been proposed over the years to cope with these challenges, aiming to gradually reduce milk yield before dry-off, promoting at the same time the start of mammary involution. Among them, the most common are based on feed or nutrient restriction, a decrease in milking frequency, or administration of prolactin inhibitors. These practices have different capacities to reduce milk yield through different mechanisms and entail several implications for udder health, animal welfare, behavior, endocrine status, metabolism, and inflammatory conditions. The present review aims to provide a comprehensive overview of the dry-off phase in high-yielding cows and of the impact of high milk production at dry-off, and to describe possible strategies that might be implemented by farmers and veterinarians to optimize this critical phase in an integrated way.

Effect of voluntary waiting period on metabolism of dairy cows during different phases of the lactation

An extended calving interval (CInt) by extending the voluntary waiting period (VWP) could be associated with altered metabolism in dairy cows. The aim of this study was first to evaluate the effects of VWP on metabolism and body condition during the first 305 d after the first calving in the experiment (calving 1), around the end of the VWP, and during pregnancy (280 d before calving 2). Second, the effects of the VWP on metabolism were determined from 2 wk before until 6 wk after calving 2. Third, individual cow characteristics were used to predict milk production and body condition of cows after different VWP. Holstein-Friesian cows (N = 154, 41 primiparous [PP], 113 multiparous [MP]) were blocked for parity, milk production, and lactation persistency, randomly assigned to a VWP of 50, 125, or 200 d (VWP50, VWP125, or VWP200) and followed from calving 1 until 6 wk after calving 2. In the first 6 wk after calving 1 and from 2 wk before until 6 wk after calving 2, weekly plasma samples were analyzed for nonesterified fatty acids (NEFA), β-hydroxybutyrate, glucose, insulin, and insulin-like growth factor 1 (IGF-1). From wk 7 after calving 1 until 2 wk before calving 2, insulin and IGF-1 were analyzed every 2 wk. Fat- and protein-corrected milk (FPCM) and body weight (BW) gain were measured weekly. Cows were divided in two parity classes based on calving 1 (PP and MP) and remained in these classes after calving 2. During pregnancy, MP cows in VWP200 had greater plasma insulin and IGF-1 concentration and lower FPCM compared with MP cows in VWP125 (insulin: 18.5 vs. 13.9 µU/mL, CI 13.0-19.7, P < 0.01; IGF-1: 198.5 vs. 175.3 ng/mL ± 5.3, P = 0.04; FPCM: 22.6 vs. 30.0 kg/d ± 0.8, P < 0.01) or VWP50 (insulin: 15.8 µU/mL, P < 0.01; IGF-1: 178.2 ng/mL, P < 0.01; FPCM: 26.6 kg/d, P < 0.01) and had a greater daily BW gain compared with cows in VWP50 (3.6 vs. 2.5 kg/d ± 0.2; P < 0.01). After calving 2, MP cows in VWP200 had greater plasma NEFA concentration (0.41 mmol/liter) compared with MP cows in VWP125 (0.30 mmol/liter, P = 0.04) or VWP50 (0.26 mmol/liter, P < 0.01). For PP cows, the VWP did not affect FPCM or body condition during the first lactation in the experiment, or metabolism after calving 2. Independent of the VWP, higher milk production and lower body condition before insemination were associated with higher milk production and lower body condition at the end of the lactation. Variation in these characteristics among cows could call for an individual approach for an extended VWP.

Effects of the Physiological Status and Diet on Blood Metabolic Parameters in Amiata Dairy Donkeys

Simple Summary

This study investigated the physiological responses of donkeys feeding two different diets during maintenance, pregnancy and lactation. To investigate how the metabolic state influences dairy production, we also assessed the relationships between the blood metabolic profile and milk quality. We found that pregnancy and the first months of lactation led to lower feed intake and increases in blood non-esterified fatty acids (NEFAs), which was linked to the difficulty that the jennies have in meeting the higher energy needs. The intake of the concentrate in the diet was associated with the increase in blood glucose, both during maintenance and lactation. Higher NEFA were related to lower de novo synthesized milk fatty acids, whereas higher plasma B-HBA were related to higher somatic cell count. This paper contributes to the knowledge of physiological responses of jennies during milk production phases and provides information for donkey milk producers.

Abstract

Body weight changes and blood metabolic parameters in jennies feeding two different diets and in three physiological statuses were investigated (maintenance vs. pregnancy; maintenance vs. lactation). The relationships between blood metabolic profile and milk quality were also evaluated. Fourteen jennies were allocated to two groups (1: pregnant/lactating; 2: non-pregnant, non–lactating). Pregnant jennies and maintenance jennies (during the first 10-week measurement period) fed a diet consisted of ad libitum grass hay (diet 1); lactating jennies and maintenance jennies (during the last 10-week measurement period) fed ad libitum grass hay plus 2 kg/head/day of concentrate (diet 2). Blood sampling was performed on the jennies of both groups; individual milk samples were also collected during the first 70 days in milk. Higher blood NEFA (p < 0.05) were found in pregnant compared to maintenance jennies (diet 1) (68 vs. 37 μmol/L). Lactating jennies showed higher (p < 0.01) average blood NEFA (268 vs. 26 μmol/L) and glucose (66 vs. 55 mg/dL) compared to the maintenance (diet 2). Blood glucose was positively correlated to milk fat (p < 0.05), while negative significant correlations between de novo milk fatty acids and NEFAs were observed. Positive correlations between plasma B-HBA and somatic cell count (p < 0.01) were also found.

Genetic analysis of lifetime productivity traits in goats

As part of the development of a breeding programme for dairy goats to support sustainable production in organic farming, the overall aim of the present study was to identify traits that can be used as selection criteria for lifetime productivity. The breeding goal is high lifetime productivity with a good milk composition and a good level of robustness in the animals, especially within grazing systems. The lifetime productivity traits analysed in this study were the length of the animal's productive life (LPL), the lifetime efficiency (LEF), and the animal's milk yield efficiency with respect the total number of lactating days (EDM); the average fat and protein content over the animal's lifetime, the fat-to-protein ratio (FPR), and the urea content (UC) were also included as indirect health traits and potential indicators of robustness in dairy goats. The traits' influencing factors, phenotypic and genetic correlations, and heritability were examined. Furthermore, factors influencing milk yield in the first 120 d of lactation during the animal's first lactation were analysed. The aim of investigating milk yield during the first lactation was to consider a connection between early performance recoding in the life of an animal and LPL, LEF, and EDM. In total, lactation and pedigree data from 9192 dairy goats of the common German Fawn (GF) and German White (GW) dairy breeds were used. Prerequisites were that the investigated birth cohorts had to have definitively completed their lifetime production, and a high proportion of goats had to have completed extended lactation. The data analysis showed that breed did not influence milk yield. The age at first kidding, the average number of kids born during the animal's lifetime, and the lactation length did influence the milk yield. This applies to the milk yield during the first 120 d of the first lactation as well as over the lifetime of an animal. Considering the influencing factors, the results showed that LPL was genetically and positively correlated with LEF and EDM ( r g = 0.65 ± 0.06 and 0.29 ± 0.07 respectively). The heritability of LPL, LEF, and EDM was 0.22 ± 0.02 , 0.29 ± 0.03 , and 0.44 ± 0.03 respectively. Regarding the lifetime milk composition, the heritability of protein and fat content, FPR, and UC was 0.63 ± 0.02 , 0.52 ± 0.02 , 0.32 ± 0.03 , and 0.47 ± 0.04 respectively. The heritability regarding the milk yield during the first 120 d of the first lactation was 0.34 ± 0.03 . We found that the milk yield during the first 120 d of the first lactation showed a genetic correlation with LPL, LEF, and EDM of 0.30 ± 0.08 , 0.82 ± 0.04 , and 0.89 ± 0.03 respectively. In summary, LPL, LEF, and EDM are suitable traits to indicate lifetime productivity in dairy goats. An additional indicator for lifetime productivity could be the milk yield during the first 120 d of the first lactation. Moreover, FPR and UC appear to be promising indicator traits for the health and robustness of dairy goats. The present study showed the importance of considering extended lactation in selective breeding programmes as well as the importance of modelling extended lactation in the breeding value estimation.

Assessment of Performance and Some Welfare Indicators of Cows in Vietnamese Smallholder Dairy Farms

Smallholder dairy farms (SDFs) are distributed widely across lowland and highland regions in Vietnam, but data on the productivity and welfare status of these cows remains limited. This cross-sectional study was conducted to describe and compare the productivity and welfare status of SDF cows across contrasting regions. It was conducted in autumn 2017 on 32 SDFs randomly selected from four typical but contrasting dairy regions (eight SDFs per region); a south lowland, a south highland, a north lowland, and a north highland region. Each farm was visited over a 24-h period (an afternoon followed by a morning milking and adjacent husbandry activities) to collect data of individual lactating cows (n = 345) and dry cows (n = 123), which included: milk yield and concentrations, body weight (BW), body condition score (BCS, 5-point scale, 5 = very fat), inseminations per conception, and level of heat stress experienced (panting score, 4.5-point scale, 0 = no stress). The high level of heat stress (96% of lactating cows were moderate to highly heat-stressed in the afternoon), low energy corrected milk yield (15.7 kg/cow/d), low percentage of lactating cows (37.3% herd), low BW (498 and 521 kg in lactating and dry cows, respectively), and low BCS of lactating cows (2.8) were the most important productivity and welfare concerns determined and these were most serious in the south lowland. By contrast, cows in the north lowland, a relatively hot but new dairying region, performed similarly to those in the south highland; a region historically considered to be one of the most suitable for dairy cows in Vietnam due to its cool environment. This indicates the potential to mitigate heat stress through new husbandry strategies. Cows in the north highland had the highest BW (535 and 569 kg in lactating and dry cows, respectively) and the highest energy corrected milk yield (19.2 kg/cow/d). Cows in all regions were heat-stressed during the daytime, although less so in the highlands compared to the lowlands. Opportunities for research into improving the productivity and welfare of Vietnamese SDF cows are discussed.

Fertility and milk production on commercial dairy farms with customized lactation lengths

Drying-off, calving, and start of lactation are critical transition events for a dairy cow. As a consequence, most animal health issues occur during these periods. By extending the voluntary waiting period for first insemination after calving, calving interval (CInt) can be extended, with possible positive effects for fertility and health. Some cows might be better suited for an extended CInt than others, due to differences in milk yield level, lactation persistency, or health status, which would justify a customized CInt based on individual cow characteristics. This study aims to investigate 13 farms with customized CInt, with respect to calving to first service interval (CFSI), accomplished CInt, services per conception (SC), conception rate at first artificial insemination (CR1AI), peak yield, lactation persistency, 305-d yield, and effective lactation yield. In total, 4,858 complete lactations of Holstein Friesian cows between 2014 and 2019 from the 13 farms were grouped by parity (1 or 2+) and CFSI (CFSI class; CFSI-1 < 84; 84 ≤ CFSI-2 < 140; 140 ≤ CFSI-3 < 196; 196 ≤ CFSI-4 < 252, CFSI-5 ≥ 252 d) or CInt (CInt class; CInt-1 < 364; 364 ≤ CInt-2 < 420; 420 ≤ CInt-3 < 476; 476 ≤ CInt-4 < 532, CInt-5 ≥ 532 d). Cow inseminations, available for 11 out of 13 farms (3,597 complete lactations), were grouped by parity (1 and 2+) and CFSI class or CInt class. The fertility and milk production characteristics were analyzed with generalized and general linear mixed models. The CFSI class was not associated with SC, but extended CInt class was associated with increased SC (CInt-1-5; 1.11-3.70 SC). More than 50% of cows in the CFSI class <84 d ended up in longer than expected CInt (>364 d), showing that these cows were not able to conceive for the desired CInt. More than 50% of cows in CInt classes 3 and higher (CInt ≥ 420 d) had an earlier first insemination before successful insemination (CFSI class 1; <196 d), showing that these extended CInt classes consisted of both cows with an extended waiting period for first insemination and cows that failed to conceive at earlier insemination(s). On most farms, lactation persistency was greatest in CInt class 1 (<364 d), probably related to the low peak yield in this class. When this shortest CInt class was excluded, persistency increased with extended CInt classes on most farms. Although at the majority of farms 305-d yield was greater in CInt ≥ 532 d, effective lactation yield at most farms was greatest in CInt from 364 to 531 d, especially for multiparous cows. Based on the results of this study, a CInt between 364 and 531 days seems most optimal for milk production, when high-yielding cows were selected.

Effects of micronutrient supplementation on performance and epigenetic status in dairy cows

The postpartum period is crucial in dairy cows and is marked by major physiological and metabolic changes that affect milk production, immune response and fertility. Nutrition remains the most important lever for limiting the negative energy balance and its consequences on general health status in highly selected dairy cows. In order to analyze the effect of a commercial micronutrient on intrinsic parameters, performances and the epigenome of dairy cows, 2 groups of 12 Holstein cows were used: 1 fed a standard diet (mainly composed of corn silage, soybean meal and non-mineral supplement) and the other 1 fed the same diet supplemented with the commercial micronutrient (µ-nutrient supplementation) for 4 weeks before calving and 8 weeks thereafter. Milk production and composition, BW, body condition score (BCS), DM intake (DMI) and health (calving score, metritis and mastitis) were recorded over the study period. Milk samples were collected on D15 and D60 post-calving for analyses of casein, Na+ and K+ contents and metalloprotease activity. Milk leukocytes and milk mammary epithelial cells (mMECs) were purified and counted. The viability of mMECs was assessed, together with their activity, through an analysis of gene expression. At the same time points, peripheral blood mononuclear cells (PBMCs) were purified and counted. Using genomic DNA extracted from PBMCs, mMECs and milk leukocytes, we assessed global DNA methylation (Me-CCGG) to evaluate the epigenetic imprinting associated with the µ-nutrient-supplemented diet. The µ-nutrient supplementation increased BCS and BW without modifying DMI or milk yield and composition. It also improved calving condition, reducing the time interval between calving and first service. Each easily collectable cell type displayed a specific pattern of Me-CCGG with only subtle changes associated with lactation stages in PBMCs. In conclusion, the response to the µ-nutrient supplementation improved the body condition without alteration of global epigenetic status in dairy cows.

Review: extended lactation in dairy cattle

This paper reviews the effects of extended lactation (EXT) as a strategy in dairy cattle on milk production and persistency, reproduction, milk quality, lifetime performance of the cow and finally the economic effects on herd and farm levels as well as the impact on emission of greenhouse gas at product level. Primiparous cows are able to produce equal or more milk per feeding day during EXT compared with a standard 305-d lactation, whereas results for multiparous cows are inconsistent. Cows managed for EXT can achieve a higher lifetime production while delivering milk with unchanged or improved quality properties. Delaying insemination enhances mounting behaviour and allows insemination after the cow's energy balance has become positive. However, in most cases EXT has no effect or a non-significant positive effect on reproduction. The EXT strategy sets off a cascade of effects at herd and farm level. Thus, the EXT strategy leads to fewer calvings and thereby expected fewer diseases, fewer replacement heifers and fewer dry days per cow per year. The optimal lifetime scenario for milk production was modelled to be an EXT of 16 months for first parity cows followed by an EXT of 10 months for later lactations. Modelling studies of herd dynamics indicate a positive effect of EXT on lifetime efficiency (milk per dry matter intake), mainly originating from benefits of EXT on daily milk yield in primiparous cows and the reduced number of replacement heifers. Consequently, EXT also leads to reduced total meat production at herd level. For the farmer, EXT can give the same economic return as a traditional lactation period. At farm level, EXT can contribute to a reduction in the environmental impact of dairy production, mainly as a consequence of the reduced production of beef. A wider dissemination of the EXT concept will be supported by methods to predict which cows may be most suitable for EXT, and clarification of how milking frequency and feeding strategy through the lactation can be organised to support milk yield and an appropriate body condition at the next calving.

Tracing recombinant bovine somatotropin ab(use) through transcriptomics: the potential of bovine somatic cells in a multi-dose longitudinal study

In the European Union, the use of recombinant bovine somatotropin (rbST) in dairy cattle is forbidden. Monitoring rbST (ab)use by its direct detection in animal matrices still remains a challenging task. New monitoring methods based on indirect detection of the substance are necessary. A new transcriptomic system based on the use of high-throughput real-time PCR in combination with somatic cells was developed to control rbST administration in dairy animals. A total of nine cows, separated into control and rbST-treated groups, were included in the study. A subcutaneous injection containing 500 mg of rbST was administered to the treated group every 14 days, up to a total of 12 doses. Milk somatic cells (MSCs) were sampled from each animal at different time points throughout 8 months of study. It was possible to obtain the transcriptomic profile of 18 genes in MSCs of rbST-treated and control groups, and using univariate and multivariate statistical analysis control and treated animals were discriminated. The transcription of CCND1, IGF-1R, TNF and IL-1β genes resulted strongly influenced by rbST treatment. The combination of MSCs, transcriptomic tools and statistical analysis has allowed the selection of four genes as potential biomarkers that could be used in a transcriptomic panel for monitoring rbST administration in cows.