Enhancing bovine abortion surveillance: A learning experience.

Abortions and perinatal mortalities (APM) substantially affect cattle industry efficiency. Various infectious and noninfectious factors have been associated with bovine APM worldwide. Infections are often considered pivotal due to their abortifacient potential, leading laboratories to primarily investigate relevant infectious agents for APM cases. Some infectious causes, such as Brucella abortus, have also a zoonotic impact, necessitating monitoring for both animal and human health. However, underreporting of bovine APM is a global issue, affecting early detection of infectious and zoonotic causes. Previous studies identified factors influencing case submission, but regional characteristics may affect results. In Belgium, farmers are obliged to report cases of APM within the context of a national brucellosis monitoring program. The inclusion criteria for this monitoring program cover abortions (gestation length of 42-260 d) and perinatal mortalities of (pre)mature calves following a gestation length of more than 260 d, which were stillborn or died within 48 h after birth. The objective of the present study was to describe the evolution in submission of APM cases within a mandatory abortion monitoring program in relation to subsidized initiatives in the northern part of Belgium. Based on the proportion of APM submissions versus the proportion of bovine reproductive females, an APM proportion (APMPR) was calculated, and factors at both animal and herd level that may influence this APMPR were explored by using linear models. This evaluation revealed that the APMPR increased with the introduction of an extensive analytical panel of abortifacient agents and a free on-farm sample collection from 0.44% to 0.94%. Additionally, an increase of the APMPR was associated with an outbreak of an emerging abortifacient pathogen (Schmallenberg virus; 1.23%), and the introduction of a mandatory eradication program for bovine viral diarrhea virus (BVDv; 1.20%). The APMPR was higher in beef compared with dairy cattle, and it was higher in winter compared with fall, spring, and summer. Smaller herds categorized in the first quartile had a higher APMPR compared with larger herds. Herds that submitted an APM in the previous year had a higher APMPR in the next year compared with herds without an APM submission. Finally, herds for which there was evidence of the presence of BVDv had a higher APMPR compared with herds without evidence of the presence of BVDv. In conclusion, the number of APM submissions increased after the introduction of a free on-farm sample collection and an extensive pathogen screening panel. Production type (beef), season (winter), smaller herd size, previous APM, and presence of BVDv seemed to have a positive effect on APMPR. However, even under mandatory circumstances, APM still seems to be underreported, since the APMPR was lower than the expected minimal rate of 2%. Therefore, further research is necessary to identify the drivers that convince farmers to submit APM cases to improve submission rates and ensure an efficient monitoring program for APM and eventually associated zoonotic pathogens.

Milk yield residuals and their link with the metabolic status of dairy cows in the transition period.

The transition period is one of the most challenging periods in the lactation cycle of high-yielding dairy cows. It is commonly known to be associated with diminished animal welfare and economic performance of dairy farms. The development of data-driven health monitoring tools based on on-farm available milk yield development has shown potential in identifying health-perturbing events. As proof of principle, we explored the association of these milk yield residuals with the metabolic status of cows during the transition period. Over 2 yr, 117 transition periods from 99 multiparous Holstein-Friesian cows were monitored intensively. Pre- and postpartum dry matter intake was measured and blood samples were taken at regular intervals to determine ß-hydroxybutyrate, nonesterified fatty acids (NEFA), insulin, glucose, fructosamine, and IGF1 concentrations. The expected milk yield in the current transition period was predicted with 2 previously developed models (nextMILK and SLMYP) using low-frequency test-day (TD) data and high-frequency milk meter (MM) data from the animal's previous lactation, respectively. The expected milk yield was subtracted from the actual production to calculate the milk yield residuals in the transition period (MRT) for both TD and MM data, yielding MRTTD and MRTMM. When the MRT is negative, the realized milk yield is lower than the predicted milk yield, in contrast, when positive, the realized milk yield exceeded the predicted milk yield. First, blood plasma analytes, dry matter intake, and MRT were compared between clinically diseased and nonclinically diseased transitions. MRTTD and MRTMM, postpartum dry matter intake and IGF1 were significantly lower for clinically diseased versus nonclinically diseased transitions, whereas ß-hydroxybutyrate and NEFA concentrations were significantly higher. Next, linear models were used to link the MRTTD and MRTMM of the nonclinically diseased cows with the dry matter intake measurements and blood plasma analytes. After variable selection, a final model was constructed for MRTTD and MRTMM, resulting in an adjusted R2 of 0.47 and 0.73, respectively. While both final models were not identical the retained variables were similar and yielded comparable importance and direction. In summary, the most informative variables in these linear models were the dry matter intake postpartum and the lactation number. Moreover, in both models, lower and thus also more negative MRT were linked with lower dry matter intake and increasing lactation number. In the case of an increasing dry matter intake, MRTTD was positively associated with NEFA concentrations. Furthermore, IGF1, glucose, and insulin explained a significant part of the MRT. Results of the present study suggest that milk yield residuals at the start of a new lactation are indicative of the health and metabolic status of transitioning dairy cows in support of the development of a health monitoring tool. Future field studies including a higher number of cows from multiple herds are needed to validate these findings.

An in-depth investigation of the microbiota and its virulence factors associated with severe udder cleft dermatitis lesions.

Udder cleft dermatitis (UCD) is a skin condition affecting the anterior parts of the udder in dairy cattle. In the present study, we aimed to shed light on the microbiota in severe UCD lesions versus healthy udder skin by putting forward a taxonomic and functional profile based on a virulence factor analysis. Through shotgun metagenomic sequencing, we found a high proportion of bacteria in addition to a low abundance of archaea. A distinct clustering of healthy udder skin versus UCD lesion samples was shown by applying principal component analysis and (sparse) partial least squares analysis on the metagenomic data. Proteobacteria, Bacillota, and Actinomycetota were among the most abundant phyla in healthy udder skin samples. In UCD samples, Bacteroidota was the most abundant phylum. At genus level, Bifidobacterium spp. had the highest relative abundance in healthy skin samples, whereas Porphyromonas spp. and Corynebacterium spp. had the highest relative abundance in UCD samples. In the differential abundance analysis, Porphyromonas spp. and Bacteroides spp. were significantly differentially abundant in UCD samples, whereas Bifidobacterium spp., Staphylococcus sp. AntiMn-1, and Staphylococcus equorum were more commonly found in healthy samples. Moreover, the abundance of several treponeme phylotypes was significantly higher in lesion samples. The streptococcal cysteine protease speB was among the most abundant virulence factors present in severe UCD lesions, while a plethora of virulence factors such as the antitoxin relB were downregulated, possibly contributing to creating the ideal wound climate for the dysbiotic community. Network analysis showed healthy lesion samples had a large network ofpositive, correlations between the abundances of beneficial species such as Aerococcus urinaeequi and Bifidobacterium angulatum, indicating that the healthy skin microbiome forms an active protective bacterial network, which is disrupted in case of UCD. In UCD samples, a smaller microbial network mainly consisting of positive correlations between the abundances of Bacteroides fragilis and anaerobic Bacteroidota was exposed. Moreover, a high correlation between the taxonomic data and virulence factors was revealed, concurrently with 2 separate networks of microbes and virulence factors. One network, matching with the taxonomic findings in the healthy udder skin samples, showcased a community of harmless or beneficial bacteria, such as Bifidobacterium spp. and Butyrivibrio proteoclasticus, associated with hcnB, hcnC, relB, glyoxalase, and cupin 2. The other network, corresponding to UCD samples, consisted of pathogenic or facultative pathogenic and mainly anaerobic bacteria such as Treponema spp., Mycoplasmopsis spp., and bovine gammaherpesvirus 4, that correlated with virulence factors SpvB, fhaB, and haemagglutination activity domain-associated factor. Our results point toward a dysbiotic community with a notable decrease in diversity and evenness, with a loss of normal skin inhabitants and innocuous or useful species making way for predominantly anaerobic, facultative pathogens. The shift in the abundance of virulence factors such as fhaB and SpvB could play a role in the manifestation of a local micro-environment favorable to the microbiome associated with udder skin lesions. Lastly, the presence of specific networks between microbial species, and between microbes and virulence factors was shown.

The effect of maternal supply of rumen-protected protein to Holstein Friesian cows during the dry period on the transfer of passive immunity and colostral microbial composition.

The objective of this study was to analyze if maternal supply of rumen-protected protein during the dry period can affect the IgG concentration and microbial composition of colostrum and the IgG absorption and fecal microbial composition in the calf. Seventy-four multiparous Holstein Friesian (HF) dairy cows were stratified per parity and randomly assigned to one of 2 different dry period diets, a diet with a low crude protein (CP) level (LP) and a diet with a high CP level (HP) by addition of rumen-undegraded protein (RUP; formaldehyde-treated soybean meal, Mervobest, Nuscience, Drongen, Belgium). Colostrum was collected within 1 h after calving and IgG concentration was quantified by radial immunodiffusion analysis. Forty-nine calves (23 female and 26 male) were enrolled in the trial with a 2 × 2 factorial design, with prenatal and postnatal treatment as the 2 independent variables. This led to 4 experimental groups: LPLP, LPHP, HPLP, and HPHP, in which the first 2 letters refer to the prenatal treatment (diet of the dam) and the last 2 refer to the postnatal treatment (diet of the colostrum-producing cow). Calves received 3× 2 L of colostrum within 2, 6, and 24 h after birth. Meconium and feces were collected solely from female calves (n = 18) by digital palpation of the rectum, immediately after birth and before colostrum administration and at d 3 of age. Microbial DNA was extracted from meconium (n = 9), feces (n = 15), and colostrum (n = 49). Amplicon sequencing of the bacterial V3-V4 region of the 16S rRNA gene was performed for characterization of the bacterial communities. Colostrum IgG concentration was higher in cows that were supplemented with RUP, especially in cows entering their second lactation (LSM ± SEM 61.3 ± 2.3 vs. 55.2 ± 2.8 g of IgG/L). Calves born out of LP cows that received colostrum from HP cows (LPHP) had a lower serum IgG level compared with HPHP and LPLP calves (LSM ± SEM 14.2 ± 1.3 vs. 18.8 ± 1.2 and 20.9 ± 1.3 g of IgG/L in HPHP and LPLP, respectively). The most abundant phyla in colostrum were Proteobacteria (48.2%), Firmicutes (24.8%), Bacteroidetes (9.5%), and Actinobacteria (5.0%). The most abundant phyla in calf meconium and feces were Firmicutes (42.5 and 47.5%), Proteobacteria (21.7% and 33.7%), Bacteroidetes (16.8% and 15.7%), and Actinobacteria (2.9% and 3.1%). There was no difference in the overall microbial communities between colostrum from HP and LP cows. However, 2 genera (both members of the family Lachnospiraceae) were more abundant in colostrum from HP cows compared with LP cows. The microbial composition of meconium, feces and colostrum differed from each other. Fecal samples were more similar to each other and are characterized by a lower intersample diversity compared with colostrum and meconium samples. To conclude, increasing the CP level by addition of RUP in the dry period diet affected the colostrum IgG concentration and the transfer of passive immunity, but did not change the overall microbial composition of colostrum nor of meconium and feces in the calf.

Exploring the microbial composition of Holstein Friesian and Belgian Blue colostrum in relation to the transfer of passive immunity.

For centuries, multicellular organisms have lived in symbiosis with microorganisms. The interaction with microorganisms has been shown to be very beneficial for humans and animals. During a natural birth, the initial inoculation with bacteria occurs when the neonate passes through the birth canal. Colostrum and milk intake are associated with the acquisition of a healthy gut flora. However, little is known about the microbial composition of bovine colostrum and the possible beneficial effects for the neonatal calf. In this prospective cohort study, the microbial composition of first-milking colostrum was analyzed in 62 Holstein Friesian (HF) and 46 Belgian Blue (BB) cows by performing amplicon sequencing of the bacterial V3-V4 region of the 16S rRNA gene. Calves received, 3 times, 2 L of their dam's colostrum within 24 h after birth. Associations between colostral microbial composition and its IgG concentration, as well as each calf's serum IgG levels, were analyzed. Colostrum samples were dominated by the phyla Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria. The 10 most abundant genera in the complete data set were Acinetobacter (16.2%), Pseudomonas (15.1%), a genus belonging to the Enterobacteriaceae family (4.9%), Lactococcus (4.0%), Chryseobacterium (3.9%), Staphylococcus (3.6%), Proteus (1.9%), Streptococcus (1.8%), Enterococcus (1.7%), and Enhydrobacter (1.5%). The remaining genera (other than these top 10) accounted for 36.5% of the counts, and another 8.7% were unidentified. Bacterial diversity differed significantly between HF and BB samples. Within each breed, several genera were found to be differentially abundant between colostrum of different quality. Moreover, in HF, the bacterial composition of colostrum leading to low serum IgG levels in the calf differed from that of colostrum leading to high serum IgG levels. Results of the present study indicate that the microbes present in colostrum are associated with transfer of passive immunity in neonatal calves.

Effect of overconditioning on the hepatic global gene expression pattern of dairy cows at the end of pregnancy.

Overconditioning is a risk factor for upregulated pre- and postpartum fat mobilization. Therefore, we hypothesized that overconditioning at the end of pregnancy leads to the accumulation of lipids in the liver and modifications of the hepatic gene expression pattern. The aim of this study was to evaluate the effect of normal- versus overconditioning on the hepatic transcriptomic profile of dairy cows at the end of pregnancy. Ten dry multiparous Holstein cows were killed 2 wk before expected calving. Body condition score (BCS) and backfat thickness (BFT) were evaluated, and blood samples for nonesterified fatty acids (NEFA) were taken before cows were killed. After cows were killed, liver biopsy samples were collected for further assessment of total lipids and RNA sequencing. Five cows were classified as normal-conditioned (median BCS = 3, range 2.75-3.5) and 5 as overconditioned (median BCS = 4, range 4-5). Regression models confirmed that normal-conditioned cows had lower BFT (1.29 ± 0.29 cm; least squares means ± standard error) and serum NEFA (0.16 ± 0.04 mmol/L) in comparison to overconditioned cows (3.14 ± 0.43 cm and 0.38 ± 0.07 mmol/L for BFT and NEFA, respectively). Total liver lipid percentage tended to be lower in normal- versus overconditioned cows (4.63 ± 0.40% and 6.06 ± 0.44%, respectively). In comparison to the mean liver lipid percentage of the normal- and overconditioned cows, 1 overconditioned cow had a relatively low (5.21%) and 1 normal-conditioned cow had a relatively high (6.07%) liver lipid percentage. Differentially expressed genes analysis (edgeR quasi-likelihood method) showed that normal-conditioned cows presented 11 upregulated and 12 downregulated genes in comparison to overconditioned cows. Linear discriminant analysis effects size revealed 133 differentially expressed genes between normal- versus overconditioned cows. Notably, the liver of normal-conditioned cows had upregulated genes associated with liver functionality (ALB, SELENOP, IGF1, and IGF2). On the other hand, overconditioned cows had upregulated genes associated with the acute-phase response (C3, HPX, and, LBP). High basal lipolysis in overconditioned cows at the end of pregnancy increased liver lipid content, and this may alter the hepatic gene expression pattern to a pro-inflammatory state.

Low microbial biomass within the reproductive tract of mid-lactation dairy cows: A study approach.

The microbiome from the reproductive tract is being investigated for its putative effect on fertility, embryo development, and health status of the human or animal host postpartum. Besides the presence of a vaginal microbiome, recent studies have claimed the existence and putative role of the uterine microbiome. Yet, the extremely low bacterial numbers and high eukaryotic/prokaryotic DNA ratio make this a highly challenging environment to study with next-generation sequencing (NGS) techniques. Here, we describe the methodological challenges that are typically encountered when performing an accurate analysis of low microbial biomass samples, illustrated by data of our own observational study. In terms of the research question, we compared the microbial composition throughout different parts of the reproductive tract of clinically healthy, mid-lactation Holstein-Friesian cows. Samples were collected from 5 dairy cows immediately after killing. Swabs were taken from the vagina, and from 4 pre-established locations of the uterine endometrium. In addition to the conventional DNA extraction blank controls, sterile swabs rubbed over disinfected disposable gloves and the disinfected surface of the uterus (tunica serosa) before incision were taken as sampling controls. The DNA extraction, DNA quantification, quantitative PCR of the 16S rRNA genes, and 16S rRNA gene sequencing were performed. In terms of NGS data analysis, we performed prevalence-based filtering of putative contaminant operational taxonomic units (OTU) using the decontam R package. Although the bacterial composition differed between the vagina and uterus, no differences in bacterial community structure (α and ß diversity) were found among the different locations in the uterus. At phylum level, uterine samples had a greater relative abundance of Proteobacteria, and a lesser relative abundance of Firmicutes than vaginal samples. The number of shared OTU between vagina and uterus was limited, suggesting the existence of bacterial transmission routes other than the transcervical one to the uterus. The mid-lactation bovine genital tract is a low microbial biomass environment, which makes it difficult to distinguish between its constitutive versus contaminant microbiome. The integration of key controls is therefore strictly necessary to decrease the effect of accidentally introduced contaminant sequences and improve the reliability of results in samples with low microbial biomass.

Invited review: MicroRNAs in bovine colostrum-Focus on their origin and potential health benefits for the calf.

Colostrum is the first milk produced by a cow after she gives birth. Compared with mature milk, it has a high concentration of immunoglobulin G. Calves are born without circulating antibodies, thus ingestion of colostrum is necessary to protect the calf against pathogens in the first challenging weeks of life. In addition to the life-saving supply of antibodies, colostrum contains minerals, vitamins, growth factors, and immune cells. Recently, microRNAs (miRNAs) were added to that list. MicroRNAs are short, non-coding RNA molecules that can regulate gene expression at the post-transcriptional level. They are thought to act as key regulators of diverse biological and developmental processes. Colostrum contains higher amounts of miRNAs than mature milk; immune- and development-related miRNAs are prominent. Their expression pattern in milk is likely to be influenced by maternal nutrition and environment. The fat content of the maternal diet appears to have a major effect on expression of miRNAs in milk and in the neonate. The immunological state of the mammary gland seems to affect miRNA expression as well. In cows diagnosed with subclinical mastitis, alterations in the expression of miRNAs in milk have been observed. It is believed that miRNAs in colostrum and milk are signaling molecules passed from mother to newborn. They are packaged in extracellular vesicles, which makes them resistant to the harsh conditions in the gastrointestinal tract. Therefore, they can reach the small intestine, where they are absorbed and transferred into the bloodstream. MicroRNAs are important for the development of the intestines. For example, miRNAs stimulate cell viability, proliferation, and stem cell activity of the intestinal epithelium. Furthermore, miRNAs seem to act as key players in the development of the complete immune system. They can, among other things, regulate B- and T-cell differentiation and affect interleukin production of macrophages. The abundance of miRNAs in colostrum and milk and the possibility for their absorption in the intestines of the neonate supports the hypothesis that these tiny molecules are important for the development of the newborn. The probable relation of diet to the expression of miRNAs by the mother creates a possible avenue to optimize expression of miRNAs and improve neonatal maturation.

A large-scale study on the effect of age at first calving, dam parity, and birth and calving month on first-lactation milk yield in Holstein Friesian dairy cattle.

Milk yield during first lactation is an important economical trait. Age at first calving (AFC) is considered an important predictor of subsequent milk yield. In addition, both season of birth, as well as season of calving, have been shown to influence milk production, with conflicting results. Finally, higher parity of the dam has been associated with a lower performance of the offspring. The aim of the present study was to assess the effect of the above-mentioned factors based on a large-scale study and to rank the most important determinants for first-lactation milk yield. Data on 3,810,678 Holstein Friesian heifers, born in Belgium and the Netherlands between 2000 and 2015, were provided by Cooperative CRV and CRV BV (Arnhem, the Netherlands) and consisted of birth dates, calving dates, and first-lactation productions. In addition, herd, sire, and dam information was provided. Linear regression models were built with herd-calving year and sire as random effects and 305-d energy-corrected milk (ECM) yield during first lactation as outcome variable. Birth month, calving month, parity of the dam, and AFC were included as fixed effects in the model and a dominance analysis was performed to rank the associated factors according to importance. Results revealed AFC to be the most important factor (R2 = 0.047), with an increase in ECM up to an age of 33 mo. Calving month was a more important predictor than birth month (R2 = 0.010 vs. R2 = 0.002, respectively), with the highest first-lactation production in heifers calving in October to December, and the lowest in heifers calving in June and July. Birth month had a limited effect on first-lactation milk yield (R2 = 0.002), potentially masked by rearing strategies during early life. Finally, parity of the dam ≥3 was associated with a reduced ECM of the offspring (R2 = 0.002). In conclusion, our results show AFC to be an important determinant of milk yield during first lactation. In addition, seasonal patterns in milk production are seen, which should be further explored to identify the underlying mechanism.

Genome-wide association for metabolic clusters in early-lactation Holstein dairy cows.

The aim of this study was to detect the genomic region or regions associated with metabolic clusters in early-lactation Holstein cows. This study was carried out in 2 experiments. In experiment I, which was carried out on 105 multiparous Holstein cows, animals were classified through k-means clustering on log-transformed and standardized concentrations of blood glucose, insulin-like growth factor I, free fatty acids, and ß-hydroxybutyrate at 14 and 35 d in milk (DIM), into metabolic clusters, either balanced (BAL) or other (OTR). Forty percent of the animals were categorized in the BAL group, and the remainder were categorized as OTR. The cows were genotyped for a total of 777,962 SNP. A genome-wide association study was performed, using a case-control approach through the GEMMA software, accounting for population structure. We found 8 SNP (BTA11, BTA23, and BTAX) associated with the predicted metabolic clusters. In experiment II, carried out on 4,267 second-parity Holstein cows, milk samples collected starting from the first week until 50 DIM were used to determine Fourier-transform mid-infrared (FT-MIR) spectra and subsequently to classify the animals into the same metabolic clusters (BAL vs. OTR). Twenty-eight percent of the animals were categorized in the BAL group, and the remainder were classified in the OTR category. Although daily milk yield was lower in BAL cows, we found no difference in daily fat- and protein-corrected milk yield in cows from the BAL metabolic cluster compared with those in the OTR metabolic cluster. In the next step, a single-step genomic BLUP was used to identify the genomic region(s) associated with the predicted metabolic clusters. The results revealed that prediction of metabolic clusters is a highly polygenic trait regulated by many small-sized effects. The region of 36,258 to 36,295 kb on BTA27 was the highly associated region for the predicted metabolic clusters, with the closest genes to this region (ANK1 and miR-486) being related to hematopoiesis, erythropoiesis, and mammary gland development. The heritability for metabolic clustering was 0.17 (SD 0.03), indicating that the use of FT-MIR spectra in milk to predict metabolic clusters in early-lactation across a large number of cows has satisfactory potential to be included in genetic selection programs for modern dairy cows.

Prediction of metabolic clusters in early-lactation dairy cows using models based on milk biomarkers.

The aim of this study was to describe metabolism of early-lactation dairy cows by clustering cows based on glucose, insulin-like growth factor I (IGF-I), free fatty acid, and ß-hydroxybutyrate (BHB) using the k-means method. Predictive models for metabolic clusters were created and validated using 3 sets of milk biomarkers (milk metabolites and enzymes, glycans on the immunogamma globulin fraction of milk, and Fourier-transform mid-infrared spectra of milk). Metabolic clusters are used to identify dairy cows with a balanced or imbalanced metabolic profile. Around 14 and 35 d in milk, serum or plasma concentrations of BHB, free fatty acids, glucose, and IGF-I were determined. Cows with a favorable metabolic profile were grouped together in what was referred to as the "balanced" group (n = 43) and were compared with cows in what was referred to as the "other balanced" group (n = 64). Cows with an unfavorable metabolic profile were grouped in what was referred to as the "imbalanced" group (n = 19) and compared with cows in what was referred to as the "other imbalanced" group (n = 88). Glucose and IGF-I were higher in balanced compared with other balanced cows. Free fatty acids and BHB were lower in balanced compared with other balanced cows. Glucose and IGF-I were lower in imbalanced compared with other imbalanced cows. Free fatty acids and BHB were higher in imbalanced cows. Metabolic clusters were related to production parameters. There was a trend for a higher daily increase in fat- and protein-corrected milk yield in balanced cows, whereas that of imbalanced cows was higher. Dry matter intake and the daily increase in dry matter intake were higher in balanced cows and lower in imbalanced cows. Energy balance was continuously higher in balanced cows and lower in imbalanced cows. Weekly or twice-weekly milk samples were taken and milk metabolites and enzymes (milk glucose, glucose-6-phosphate, BHB, lactate dehydrogenase, N-acetyl-ß-d-glucosaminidase, isocitrate), immunogamma globulin glycans (19 peaks), and Fourier-transform mid-infrared spectra (1,060 wavelengths reduced to 15 principal components) were determined. Milk biomarkers with or without additional cow information (days in milk, parity, milk yield features) were used to create predictive models for the metabolic clusters. Accuracy for prediction of balanced (80%) and imbalanced (88%) cows was highest using milk metabolites and enzymes combined with days in milk and parity. The results and models of the present study are part of the GplusE project and identify novel milk-based phenotypes that may be used as predictors for metabolic and performance traits in early-lactation dairy cows.

Association between metabolic diseases and the culling risk of high-yielding dairy cows in a transition management facility using survival and decision tree analysis.

The objective of this study was to assess the association between individual metabolic diseases (MD) and multiple MD (MD+) in the transition period (±3 wk relative to calving) and the culling risk in the first 120 d in milk (DIM) in Holstein-Friesian dairy cows. Health records from a transition management facility in Germany with 1,946 calvings were analyzed in a 1-yr cohort via survival analysis and a decision tree model. The recorded MD were milk fever (MF), retained placenta (RP), metritis (METR), ketosis (KET), displaced abomasum (DA), twinning (TWIN), and clinical mastitis (MAST). The overall culling within 120 DIM was 18%. The 120 DIM culling risk for healthy cows (64.8% of the total) was 13%, whereas it was 25% for MD (24.5%) and 33% for MD+ (10.7%) cows. The 120 DIM culling risk (%) for each MD and MD+, respectively, was 34.6 and 48 for MF and MF+, 15 and 31.5 for RP and RP+, 9.4 and 22.2 for METR and METR+, 30.7 and 37.3 for KET and KET+, 56.1 and 46.8 for DA and DA+, 30.3 and 34 for TWIN and TWIN+, and 36.6 and 27.8 for MAST and MAST+. Moreover, the incidence risk (%) for each MD and MD+, respectively, was 4 and 2.6 for MF and MF+, 1 and 2.8 for RP and RP+, 8.7 and 6 for METR and METR+, 4.5 and 6.1 for KET and KET+, 0.8 and 2.4 for DA and DA+, 1.7 and 2.7 for TWIN and TWIN+, and 3.6 and 1.8 for MAST and MAST+. Setting the healthy cows as the referent, the 120 DIM hazard ratios (HR) for culling were MD 2.1, MD+ 2.9, MF 3.3, MF+ 4.6, RP+ 2.7, METR+ 1.8, KET 2.6, KET+ 3.3, DA 5.5, DA+ 4.5, TWIN 2.8, TWIN+ 3.0, MAST 3.1, and MAST+ 2.3. According to both decision tree and random forest analyses, MF was the most significant disease influencing survival, followed by DA, MAST, METR, and TWIN. In conclusion, the presence of MD or MD+ during the transition period was associated with increased culling risk in the first 120 DIM. The culling hazard was greater when an MD was complicated with another MD. In this study performed in a well-managed large farm, uncomplicated cases of RP (HR = 1.2) and METR (HR = 0.7) did not have an influence on the 120 DIM culling risk. Interestingly, both decision tree and random forest analyses pointed to MF and DA as main culling reasons in the first 120 DIM in the present dairy herd.

Effects of overconditioning on pancreatic insulin secretory capacity, fat infiltration, and the number and size of islets in dairy cows at the end of the dry period.

In the present study, we tested the hypothesis that overconditioning in dairy cows at the end of the dry period leads to infiltration of fat and alterations of the insulin secretory capacity of the pancreas. Pregnant Holstein Friesian dairy cows were selected based on body condition score (BCS) at the start of the dry period. Body condition score varied between cows to have optimal conditioned (2.5 < BCS ≤3.5, n = 5) and overconditioned (3.5 < BCS ≤5, n = 5) cows. All animals underwent an intravenous glucose tolerance test (IVGTT) at an average of 260 d of gestation to measure the pancreatic insulin secretory capacity and assess peripheral insulin sensitivity regarding glucose metabolism. Eight days after the IVGTT, animals were slaughtered. The pancreas was dissected and weighed and tissue samples were taken for histological analysis. Results revealed that overconditioning in dairy cows led to fat infiltration in the pancreas and an increase in size of pancreatic islets expressed relative to the total area of pancreatic tissue. In addition, results revealed a positive correlation between serum fatty acid concentration and peak insulin concentration and area and number of pancreatic islets expressed relative to the total area of pancreatic tissue. The IVGTT revealed that overconditioned animals have a higher insulin secretory capacity of the pancreas, as demonstrated by higher peak insulin concentration, higher acute insulin response to glucose, and higher area under the curve (AUC) for insulin compared with optimal conditioned cows. A higher AUC for glucose during the first 60 min following administration of the glucose bolus in overconditioned cows indicates an insulin-resistant state regarding glucose metabolism. Our results suggest that the pancreas of overconditioned dairy cows at the end of gestation compensates for the concomitantly elevated level of peripheral insulin resistance by greater secretion of insulin.

Influence of adipocyte size and adipose depot on the number of adipose tissue macrophages and the expression of adipokines in dairy cows at the end of pregnancy.

The aim of this study was to determine the number of adipose tissue macrophages (ATM) and the mRNA expression of adipokines [adiponectin (ADIPOQ), leptin (LEP), interleukin 6 (IL6), tumor necrosis factor (TNF), and interleukin 10 (IL10)] in different adipose depots from cows with a variable body condition score (BCS) at the end of the dry period. We hypothesized that the number of ATM and the expression of these adipokines depend on adipocyte size and the anatomical location of the adipose depot. Subcutaneous, omental, mesenteric, perirenal, and intrapelvic adipose tissue samples were taken immediately after euthanasia of 10 Holstein Friesian dairy cows (upcoming parity 2 to 5, age 3.9 ± 1.4 yr; mean ± standard deviation) at the end of pregnancy (actual days of pregnancy at the moment of euthanasia: 269 ± 5 d). During the dry period, all animals received similar diets to meet but not exceed requirements. Five animals were considered to have a normal BCS (2.5-3.5) and 5 animals were considered to be over-conditioned (BCS = 3.75-5). Body weight of the animals at the moment of euthanasia was 717 ± 77 kg. Expression of the different genes was determined by reverse transcription quantitative real-time PCR. Adipocyte size was determined by measuring the area of 100 adipocytes on histological sections. Average adipocyte area was 10,475 ± 1,019, 8,500 ± 780, 10,383 ± 1,227, 11,466 ± 1,039, and 11,087 ± 1,632 µm2 for the subcutaneous, mesenteric, omental, intrapelvic, and perirenal adipose depot, respectively. Immunohistochemistry using anti-bovine CD172a antibodies was performed to determine the proportion of ATM (the number of CD172a-positive cells per 100 adipocytes, given as a percentage). Expression of LEP, IL6, and TNF was positively associated with adipocyte size, whereas no association could be detected between ADIPOQ and IL10 with the size of the adipocytes. The omental adipose depot was especially infiltrated with ATM (1.92 ± 0.55, 1.10 ± 0.33, and 8.28 ± 2.24% for the subcutaneous, mesenteric, and omental adipose depot, respectively). The proportion of ATM was positively associated with the size of the adipocytes in the omental and mesenteric adipose depot. Expression of ADIPOQ, LEP, IL6, TNF, and IL10 differed among depots, which suggests differences in inflammatory characteristics depending on the anatomical location of depots. In conclusion, the results of the present study confirm the adipose tissue as a potential source of inflammatory mediators and demonstrate ATM infiltration, especially in the omental adipose depot.

Cytological endometritis diagnosed at artificial insemination in repeat breeder dairy cows.

This study aimed to (i) assess the prevalence of cytological endometritis (CYTO) diagnosed at artificial insemination (AI); (ii) evaluate the effect of CYTO on the pregnancy outcome of the same AI sample; and (iii) determine the risk factors associated with CYTO diagnosed at AI in repeat breeder (RB) dairy cows. We analysed the productive and reproductive performances of 146 RB Holstein-Friesian cows. To obtain a CYTO sample at AI, we used the cytotape technique. Generalized mixed effect models were computed to find the risk factors associated with the pregnancy and CYTO outcome. Based on ≥1% PMN cut-off point, the CYTO prevalence at AI in RB cows was 25.3%. The overall pregnancy at AI was 44.2%. The conception rate in CYTO-positive (n = 37) RB cows was 29.7% versus 49.5% for CYTO-negative (n = 109) cows. A RB cow diagnosed CYTO positive at AI had 0.47 [odds ratio (OR)] odds to become pregnant in comparison with a CYTO-negative cow. Cows that produced more milk than their counterparts in this study had increased odds (OR = 1.01) to be CYTO positive at AI. A novel risk factor positively associated with CYTO diagnosed at AI in RB cows was the level of daily milk urea (OR = 1.11). To conclude, CYTO at the moment of AI had a significantly negative effect on the pregnancy outcome in RB dairy cows. However, as only one of fourth of RB cows is affected with CYTO at AI, it may not be considered a key element associated with the RB syndrome.

Cytological endometritis at artificial insemination in dairy cows: Prevalence and effect on pregnancy outcome.

The aims of our field study in dairy cows were (1) to consolidate cytotape (CT) as a valid technique to diagnose cytological endometritis (CYTO) during artificial insemination (AI); (2) to establish a cutoff point concerning the polymorphonuclear cells (PMN) proportion to diagnose CYTO at AI; (3) to assess the prevalence of CYTO at AI; and (4) to evaluate the effect of CYTO on the pregnancy outcome of that AI. The investigation was performed using 1,625 AI-CT samples harvested from 873 Holstein-Friesian cows from 18 dairy farms in the Flemish region of Belgium. The CT device consisted of adapting a 1.5-cm piece of paper tape on the top of a conventional AI catheter covered with a double guard sheet, allowing an endometrial cytology sample to be taken when performing an AI. A receiving operator characteristic curve was built to assess the threshold level above which the PMN proportion significantly affected the AI success. Multilevel generalized mixed-effect models were built to identify factors affecting the pregnancy outcome of the AI under investigation. Only 7 samples (0.4%) harvested in 5 cows were discarded because of low-quality parameters. The cutoff point for CYTO at AI was set at ≥1% PMN (sensitivity=33.8%, specificity=88.6%). Prevalence of CYTO at AI was 27.8%. The conception rate for CYTO-positive samples was 32.7%, whereas it was 47% for CYTO-negative samples. A CYTO-negative AI had 1.8 [odds ratio (OR)] more chances to become pregnant than a CYTO-positive one. Other factors identified as detrimental for the pregnancy outcome were body condition score ≤1.5 (OR=0.6), relative 305-d milk yield (OR=0.9), dystocia (OR=0.3), parity ≥2 (OR=0.7), and warm months of the year. In conclusion, CT is a consolidated technique to diagnose CYTO at AI, PMN 1% is the threshold level to diagnose CYTO at AI, around one-quarter of inseminated uteri suffer from CYTO, and affected uteri having a significantly lower chance to become pregnant from that insemination.

Novel approaches to assess the quality of fertility data stored in dairy herd management software.

Scientific journals and popular press magazines are littered with articles in which the authors use data from dairy herd management software. Almost none of such papers include data cleaning and data quality assessment in their study design despite this being a very critical step during data mining. This paper presents 2 novel data cleaning methods that permit identification of animals with good and bad data quality. The first method is a deterministic or rule-based data cleaning method. Reproduction and mutation or life-changing events such as birth and death were converted to a symbolic (alphabetical letter) representation and split into triplets (3-letter code). The triplets were manually labeled as physiologically correct, suspicious, or impossible. The deterministic data cleaning method was applied to assess the quality of data stored in dairy herd management from 26 farms enrolled in the herd health management program from the Faculty of Veterinary Medicine Ghent University, Belgium. In total, 150,443 triplets were created, 65.4% were labeled as correct, 17.4% as suspicious, and 17.2% as impossible. The second method, a probabilistic method, uses a machine learning algorithm (random forests) to predict the correctness of fertility and mutation events in an early stage of data cleaning. The prediction accuracy of the random forests algorithm was compared with a classical linear statistical method (penalized logistic regression), outperforming the latter substantially, with a superior receiver operating characteristic curve and a higher accuracy (89 vs. 72%). From those results, we conclude that the triplet method can be used to assess the quality of reproduction data stored in dairy herd management software and that a machine learning technique such as random forests is capable of predicting the correctness of fertility data.

Short communication: Limitations of glucose tolerance tests in the assessment of peripheral tissue insulin sensitivity during pregnancy and lactation in dairy heifers.

The aim of the present study was to point at the limitations of glucose tolerance tests (GTT) to assess peripheral tissue insulin sensitivity in dairy heifers in different physiological states (pregnancy and lactation). Intravenous GTT were performed in 5 nonpregnant, nonlactating heifers, 5 heifers at the end of pregnancy (12-7 d before calving), and 5 lactating primiparous cows (11-14 d after calving). Glucose and insulin concentrations were determined and area under the curve (AUC) and clearance rate of glucose and insulin were calculated. Additionally, data were analyzed using the minimal model to derive the insulin sensitivity parameter (Si). Basal glucose and insulin concentrations were greater in the nonpregnant, nonlactating heifers. Clearance rate of glucose and Si were lowest, whereas the AUC for glucose was greatest in the nonpregnant, nonlactating heifers. Insulin concentrations during the GTT were greater for the nonpregnant, nonlactating heifers. Results from the GTT in pregnant heifers and lactating primiparous cows are biased by the fact that a large part of the glucose disappearance during an intravenous GTT occurs independently of insulin by the pregnant uterus or the lactating mammary gland. As such, greater AUC of glucose, lower clearance rate of glucose, or lower Si derived from GTT performed in nonpregnant, nonlactating dairy heifers in the present study might indicate decreased peripheral tissue insulin sensitivity of the glucose metabolism or decreased insulin-independent glucose disappearance. Based on the results from a GTT, it is impossible to discriminate between both metabolic pathways. It can be concluded that parameters derived from GTT are not suited to compare peripheral tissue insulin sensitivity of the glucose metabolism between dairy heifers in different physiological states due to the large variation in insulin secretion and the substantial difference in insulin-independent glucose disposal associated with these physiological states.

Innovative look at dairy heifer rearing: Effect of prenatal and post-natal environment on later performance.

As heifer rearing is a costly investment, dairy farmers have been stimulated to maximize early growth of their calves, mainly by enhanced liquid feeding. However, the long-term effects of this "accelerated growth" are largely unknown. Studies recently performed at Ghent University indicate that in dairy cattle, certain maternal factors (such as young age and high milk yield) and environmental factors (such as high ambient temperatures) create a suboptimal environment for the developing foetus, altering the phenotype of the newborn calf. According to the "thrifty phenotype hypothesis," these metabolic alterations prepare the newborn for similar ("matching") conditions after birth, enhancing its survival during periods of limited feeding. Yet, when an abundance of nutrients is available in post-natal life (e.g., during periods of enhanced feeding), the "mismatch" between pre- and post-natal environment results in an early catch-up growth, with potential negative consequences. The aim of the article was to discuss this mismatch between pre- and post-natal environment in dairy calves. Previous studies, especially in human medicine, have shown catch-up growth to be associated with obesity, fertility problems, metabolic diseases and a reduced lifespan. Hence, we hypothesize that, by applying programs of accelerated growth, our current management system accentuates the mismatch between the pre- and post-natal environment in dairy calves. We can conclude that, although more research is necessary, the current findings point towards a more individual approach when rearing dairy heifers.

Environmental and maternal factors associated with gross placental morphology in dairy cattle.

This article reports on a study of gross placental morphology of 282 expelled placentas from 89 primi- and 193 multiparous Holstein dams immediately after calving and examines associations with environmental factors such as typical herd features and season of calving, and maternal factors such as age at calving, level of milk yield at conception and cumulative amount of milk produced during gestation. The highest correlation between calf measurements and placental characteristics was found between the weight of the calf and the total cotyledonary surface (r = .643; p < .001), confirming the high importance of the cotyledonary surface available for nutrient transfer to the developing foetus. Younger age in adolescent and smaller heart girth in multiparous dams were associated with a higher cotyledon number, suggesting placental compensation in dams with lower capacities in terms of dry matter intake. No significant association between milk yield during gestation in multiparous animals and gross placental morphometrics could be detected, indicating that factors such as the amount of milk produced during gestation affect placental development less than foetal weight close to term. Therefore, placental growth may be sustained at the expense of other tissues in an attempt to maintain pregnancy and minimize the adverse consequences for the foetus. This study offers evidence concerning factors affecting the placental surface size for nutrient transfer from dam to calf in dairy cattle based on gross morphometrics, but needs confirmation from studies in which this surface size is more profoundly assessed by measuring the branches present in the cotyledonary villi.