Johne's disease (JD) is an infectious enteric disease in ruminants, causing substantial economic loss annually worldwide. This work aimed to estimate JD's genetic parameters and the phenotypic and genetic trends by incorporating recent data. It also explores the feasibility of a national genetic evaluation for JD susceptibility in Holstein cattle in the United States (US). The data were extracted from a Johne's disease data repository, maintained at the Council on Dairy Cattle Breeding (CDCB), and initially supplied by 2 dairy records processing centers. The data comprised 365,980 Holstein cows from 1,048 herds participating in a voluntary control program for JD. Two protocol kits, Idexx Paratuberculosis Screening Ab Test (IDX) and Parachek® 2 (PCK), were used to analyze milk samples with the Enzyme-Linked Immunosorbent Assay (ELISA) technique. Test results from the first 5 parities were considered. An animal was considered infected if it had at least one positive outcome. The overall average of JD incidence was 4.72% in these US Holstein cattle. Genotypes of 78,964 SNP markers were used for 25,000 animals randomly selected from the phenotyped population. Variance components and genetic parameters were estimated based on 3 models, namely, a pedigree-only threshold model (THR), a single-step threshold model (ssTHR), and a single-step linear model (ssLR). The posterior heritability estimates of JD susceptibility were low to moderate: 0.11 to 0.16 based on the 2 threshold models and 0.05 to 0.09 based on the linear model. The average reliability of estimated breeding values of JD susceptibility using single-step analysis for animals with or without phenotypes varied from 0.18 (THR) to 0.22 (ssLR) for IDX and from 0.14 (THR) to 0.18 (ssTHR and ssLR) for PCK. Despite no prior direct genetic selection against JD, the estimated genetic trends of JD susceptibility were negative and highly significant (P-value <0.01). The correlations of bulls' predicted transmitting abilities with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a non-overlapping genetic selection process with traits in current genetic evaluations. Our results suggest the feasibility of reducing the JD incidence rate by incorporating it into the national genetic evaluation programs.
The Resilient Dairy Genome Project (RDGP) is an international large-scale applied research project that aims to generate genomic tools to breed more resilient dairy cows. In this context, improving feed efficiency and reducing greenhouse gases from dairy is a high priority. The inclusion of traits related to feed efficiency (e.g., dry matter intake [DMI]) or greenhouse gases (e.g., methane emissions [CH4]) relies on available genotypes as well as high quality phenotypes. Currently, 7 countries (i.e., Australia, Canada, Denmark, Germany, Spain, Switzerland, and United States) contribute with genotypes and phenotypes including DMI and CH4. However, combining data are challenging due to differences in recording protocols, measurement technology, genotyping, and animal management across sources. In this study, we provide an overview of how the RDGP partners address these issues to advance international collaboration to generate genomic tools for resilient dairy. Specifically, we describe the current state of the RDGP database, data collection protocols in each country, and the strategies used for managing the shared data. As of February 2022, the database contains 1,289,593 DMI records from 12,687 cows and 17,403 CH4 records from 3,093 cows and continues to grow as countries upload new data over the coming years. No strong genomic differentiation between the populations was identified in this study, which may be beneficial for eventual across-country genomic predictions. Moreover, our results reinforce the need to account for the heterogeneity in the DMI and CH4 phenotypes in genomic analysis.
Greenhouse Gases , Female , Animals , Cattle , Genomics , Genotype , Australia , Methane
Various methods have been proposed to estimate daily yield from partial yields, primarily to deal with unequal milking intervals. This paper offers an exhaustive review of daily milk yields, the foundation of lactation records. Seminal advancements in the late 20th century concentrated on two main adjustment metrics: additive additive correction factors (ACF) and multiplicative correction factors (MCF). An ACF model provides additive adjustments to two times AM or PM milk yield, which then becomes the estimated daily yields, whereas an MCF is a ratio of daily yield to the yield from a single milking. Recent studies highlight the potential of alternative approaches, such as exponential regression and other nonlinear models. Biologically, milk secretion rates are not linear throughout the entire milking interval, influenced by the internal mammary gland pressure. Consequently, nonlinear models are appealing for estimating daily milk yields as well. MCFs and ACFs are typically determined for discrete milking interval classes. Nonetheless, large discrete intervals can introduce systematic biases. A universal solution for deriving continuous correction factors has been proposed, ensuring reduced bias and enhanced daily milk yield estimation accuracy. When leveraging test-day milk yields for genetic evaluations in dairy cattle, two predominant statistical models are employed: lactation and test-day yield models. A lactation model capitalizes on the high heritability of total lactation yields, aligning closely with dairy producers' needs because the total amount of milk production in a lactation directly determines farm revenue. However, a lactation yield model without harnessing all test-day records may ignore vital data about the shapes of lactation curves needed for informed breeding decisions. In contrast, a test-day model emphasizes individual test-day data, accommodating various intervals and recording plans and allowing the estimation of environmental effects on specific test days. In the United States, the patenting of test-day models in 1993 used to restrict the use of test-day models to regional and unofficial evaluations by the patent holders. Estimated test-day milk yields have been used as if they were accurate depictions of actual milk yields, neglecting possible estimation errors. Its potential consequences on subsequent genetic evaluations have not been sufficiently addressed. Moving forward, there are still numerous questions and challenges in this domain.
This study compared 3 correlational (best prediction, linear regression, and feed-forward neural networks) and 2 causal models (recursive structural equation model and recurrent neural networks) for estimating lactation milk yields. The correlational models assumed associations between test-day milk yields (health conditions), while the casual models postulated unidirectional recursive effects between these test-day variables. Wood lactation curves were used to simulate the data and served as a benchmark model. Individual Wood lactation curves provided an excellent parametric interpretation of lactation dynamics, with their prediction accuracies depending on the coverage of the lactation curve dynamics. Best prediction outperformed other models in the absence of mastitis but was suboptimal when mastitis was present and unaccounted for. Recurrent neural networks yielded the highest accuracy when mastitis was present. Although causal models facilitated the inference about the causality underlying lactation, precisely capturing the causal relationships was challenging because the underlying biology was complex. Misspecification of recursive effects in the recursive structural equation model resulted in a loss of accuracy. Hence, modeling causal relationships does not necessarily guarantee improved accuracies. In practice, a parsimonious model is preferred, balancing model complexity and accuracy. In addition to the choice of statistical models, the proper accounting for factors and covariates affecting milk yields is equally crucial.
In the United States, lactation milk yields are not measured directly but are calculated from the test-day milk yields. Still, test-day milk yields are estimated from partial yields obtained from single milkings. Various methods have been proposed to estimate test-day milk yields, primarily to deal with unequal milking intervals dating back to the 1970s and 1980s. The Wiggans model is a de facto method for estimating test-day milk yields in the United States, which was initially proposed for cows milked 3 times daily, assuming a linear relationship between a proportional test-day milk yield and milking interval. However, the linearity assumption did not hold precisely in Holstein cows milked twice daily because of prolonged and uneven milking intervals. The present study reviewed and evaluated the nonlinear models that extended the Wiggans model for estimating daily or test-day milk yields. These nonlinear models, except step functions, demonstrated smaller errors and greater accuracies for estimated test-day milk yields compared with the conventional methods. The nonlinear models offered additional benefits. For example, the locally weighted regression model (e.g., locally estimated scatterplot smoothing) could utilize data information in scalable neighborhoods and weigh observations according to their distance in milking interval time. General additive models provide a flexible, unified framework to model nonlinear predictor variables additively. Another drawback of the conventional methods is a loss of accuracy caused by discretizing milking interval time into large bins while deriving multiplicative correction factors for estimating test-day milk yields. To overcome this problem, we proposed a general approach that allows milk yield correction factors to be derived for every possible milking interval time, resulting in more accurately estimated test-day milk yields. This approach can be applied to any model, including nonparametric models.
Dairying , Milk , Female , Cattle , Animals , Time Factors , Dairying/methods , Lactation , Nonlinear Dynamics
Cows are typically milked 2 or more times on a test-day, but not all these milkings are sampled and weighed. The initial approach estimated a test-day yield with doubled morning (AM) or evening (PM) yield in the AM-PM milking plans, assuming equal AM and PM milking intervals. However, AM and PM milking intervals can vary, and milk secretion rates may be different between day and night. Statistical methods have been proposed to estimate daily yields in dairy cows, focusing on various yield correction factors in 2 broad categories: additive correction factors (ACF) and multiplicative correction factors (MCF). The ACF are evaluated by the average differences between AM and PM milk yield for various milking interval classes, coupled with other categorical variables. We show that an ACF model is equivalent to a regression model of daily yield on categorical regressor variables, and a continuous variable for AM or PM yield with a fixed regression coefficient of 2.0. Similarly, a linear regression model can be implemented as an ACF model with the regression coefficient for AM or PM yield estimated from the data. The linear regression models improved the accuracy of the estimates compared with the ACF models. The MCF are ratios of daily yield to yield from single milkings, but their statistical interpretations vary. Overall, MCF were more accurate for estimating daily milk yield than ACF. The MCF have biological and statistical challenges. Systematic biases occurred when ACF or MCF were computed on discretized milking interval classes, leading to accuracy loss. An exponential regression model was proposed as an alternative model for estimating daily milk yields, which improved the accuracy. Characterization of ACF and MCF showed how they improved the accuracy compared with doubling AM or PM yield as the daily milk yield. All the methods performed similarly with equal AM and PM milkings. The methods were explicitly described to estimate daily milk yield in AM and PM milking plans. Still, the principles generally apply to cows milked more than 2 times a day and apply similarly to the estimation of daily fat and protein yields with some necessary modifications.
Cost-effective milking plans have been adapted to supplement the standard supervised twice-daily monthly testing scheme since the 1960s. Various methods have been proposed to estimate daily milk yields (DMY), focusing on yield correction factors. The present study evaluated the performance of existing statistical methods, including a recently proposed exponential regression model, for estimating DMY using 10-fold cross-validation in Holstein and Jersey cows. The initial approach doubled the morning (AM) or evening (PM) yield as estimated DMY in AM-PM plans, assuming equal 12-h AM and PM milking intervals. However, in reality, AM milking intervals tended to be longer than PM milking intervals. Additive correction factors (ACF) provided additive adjustments beyond twice AM or PM yields. Hence, an ACF model equivalently assumed a fixed regression coefficient or a multiplier of "2.0" for AM or PM yields. Similarly, a linear regression model was viewed as an ACF model, yet it estimated the regression coefficient for a single milk yield from the data. Multiplicative correction factors (MCF) represented daily to partial milk yield ratios. Hence, multiplying a yield from single milking by an appropriate MCF gave a DMY estimate. The exponential regression model was analogous to an exponential growth function with the yield from single milking as the initial state and the rate of change tuned by a linear function of milking interval. In the present study, all the methods had high precision in the estimates, but they differed considerably in biases. Overall, the MCF and linear regression models had smaller squared biases and greater accuracies for estimating DMY than the ACF models. The exponential regression model had the greatest accuracies and smallest squared biases. Model parameters were compared. Discretized milking interval categories led to a loss of accuracy of the estimates. Characterization of ACF and MCF revealed their similarities and dissimilarities and biases aroused by unequal milking intervals. The present study focused on estimating DMY in AM-PM milking plans. Yet, the methods and relevant principles are generally applicable to cows milked more than two times a day.
Genetic selection has been an effective strategy to improve calving traits including stillbirth in dairy cattle. The primary objectives of the present study were to characterize stillbirth data and determine the feasibility of implementing routine genetic evaluations of stillbirth in five non-Holstein dairy breeds, namely Ayrshire, Guernsey, Milking Shorthorn, Brown Swiss, and Jersey. An updated sire-maternal grandsire threshold model was used to estimate genetic parameters and genetic values for stillbirth. Stillbirth data with the birth years of dams from 1995 to 2018 were extracted from the United States national calving ease database maintained by the Council on Dairy Cattle Breeding. The extracted stillbirth records varied drastically among the five dairy breeds. There were approximately 486K stillbirth records for Jersey and more than 80K stillbirth records for Brown Swiss. The direct and maternal heritability estimates of stillbirth were 6.0% (4.5-7.6%) and 4.7% (3.3-6.1%) in Jersey and 6.8% (3.2-10.5%) and 1.1% (0.6-2.9%) in Brown Swiss. The estimated genetic correlations between direct and maternal genetic effects for stillbirth were -0.15 (-0.38 to -0.08) in Jersey and -0.35 (-0.47 to -0.12) in Brown Swiss. The estimated genetic parameters for stillbirth in these two breeds were within close ranges of previous studies. The reliabilities of predicted transmitting abilities in Jersey and Brown Swiss increased substantially, thanks to the substantial increase in available stillbirth data in the past 10 years. The stillbirth records for Ayrshire, Guernsey, and Milking Shorthorn, which ranged approximately between 3K and 12K, are insufficient to implement reliable routine genetic evaluations of stillbirth in these three dairy breeds. Estimated genetic (co)variances and genetic values deviated considerably from the reported ranges of previous studies, and the reliabilities of predicted transmitting abilities were low in these three breeds. In conclusion, routine genetic evaluations of stillbirth are feasible in Brown Swiss and Jersey. However, reliable genetic evaluations of stillbirth in Ayrshire, Guernsey, and Milking Shorthorn require further data collection on stillbirth.
There has been increasing interest in residual feed intake (RFI) as a measure of net feed efficiency in dairy cattle. Residual feed intake phenotypes are obtained as residuals from linear regression encompassing relevant factors (i.e., energy sinks) to account for body tissue mobilization. By rearranging the single-trait linear regression, we showed a causal RFI interpretation underlying the linear regression for RFI. It postulates recursive effects in energy allocation from energy sinks on dry matter intake, but the feedback or simultaneous effects are nonexistent. A Bayesian recursive structural equation model was proposed for directly predicting RFI and energy sinks and estimating relevant genetic parameters simultaneously. A simplified Markov chain Monte Carlo algorithm was described. The recursive model is asymptotically equivalent to one-step linear regression for RFI, yet extends the analytical capacity to multiple-trait analysis.
Dry-cured meats products are consumed in various regions of the world and, consumers are increasingly demanding better quality and safety of these products. Some fungal species can produce mycotoxins in drycured meats, such as aflatoxins and ochratoxins, which, when ingested, can produce carcinogenic and mutagenic effects in humans. Contamination of these products can occur at different points of the production chain, from the field (animal contaminated with feed) to the production or storage of the final product. Although the presence of mycotoxins in drycured meats has been reported in several regions of the world, the presence of these contaminants are not legislated in most countries. Therefore, it is important to put in place methods to identify and reduce the contamination of dry-cured meats, minimizing the consumption and deleterious effects caused by mycotoxins. This review aimed to describe mycotoxin-producing fungi, mycotoxins, the relationship between human food and animal feed; legislation; incidence, identification and control methods for mycotoxins in dry-cured meats intended for human consumption.
Food Contamination/analysis , Meat Products/analysis , Mycotoxins/analysis , Animal Feed/analysis , Animals , Humans
Two experiments were carried out to assess the effects of electric and magnetic fields (EMF) on blood thyroxine (T4) in dairy cattle. In experiment 1, 16 lactating pregnant Holstein cows were exposed to 10 kV/m, 30 microTesla (microT) EMF. The animals were divided into two groups of eight animals each. Each group was exposed to EMF according to one of two treatment sequences of three periods of 28 days each. Sequence 1 was EMF OFF-ON-OFF and sequence 2 was EMF ON-OFF-ON. During the last day of each treatment period, blood samples were collected every 4 h for 24 h to estimate T4 plasma concentrations. In experiment 2, 16 nonlactating, nonpregnant, multiparous Holsteins were exposed to 10 kV/m, 30 microT EMF. The animals were divided into two groups of eight animals each. Each group was exposed to EMF according to one of the two treatment sequences described above, except that each period amounted to the number of days corresponding to one estrous cycle. During treatment, blood samples were collected every other day for T4 analysis. In both experiments, the light cycle emulated a short photoperiod (8 h light/16 h dark). During the ON periods, the animals were exposed to EMF for 16 h, 8 h of the light period plus the first 8 h of during the dark period. In experiment 1, exposed animals did not have any change in T4 plasma concentrations due to treatment (P = .0968), but, the time of sample collection revealed a significant difference (P = .0012). In experiment 2, the effect of period (P = .0009) and the treatment by days interaction (P = .0003) were statistically significant. We conclude that a worst case scenario exposure of dairy cattle to 10 kV/m, 30 microT EMF influences, in a moderate fashion, the blood levels of thyroxine.
Cattle/blood , Electricity , Electromagnetic Fields , Thyroxine/blood , Animals , Dairying , Female , Pregnancy
An exposure chamber was designed to study the effects of electric and magnetic fields (EMF) on oestrous cycles, hormonal profile during gestation, pineal function, quantity and quality of milk production, feed intake, and central nervous system of dairy cattle. The chamber was 15 x 10 x 3 m; and the control system was fully computerized so that the field intensities can be varied and monitored continuously, on site or remotely. During exposure to EMF, milk production, feed consumption, and health were monitored closely and blood and cerebral spinal fluid were continuously sampled. The chamber characteristics allow use of a wide range of exposure such as electric fields (0-30 kV/m) and magnetic fields (0-100 microT) at frequencies ranging from 45 to 3000 Hz.
Cattle/physiology , Electricity , Environment, Controlled , Lactation/physiology , Lactation/radiation effects , Milk/radiation effects , Radiation Monitoring/instrumentation , Animals , Equipment Design , Equipment Failure Analysis , Female , Radiation Dosage , Radiation Monitoring/methods , Relative Biological Effectiveness
