Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and Bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records. © 2013 American Dairy Science Association.

Influence of the physiological state on infestation by Rhipicephalus microplus in dairy cows

The purpose of this study was to identify the influence of physiological state, season, breed, number of lactations, dairy productivity, and stocking rate on infestation by Rhipicephalus microplus in dairy cows. Two hundred cows were selected through proportional representative sampling at the Seropédica Experimental Station (Pesagro-Rio), Brazil. Fully or partially engorged R. microplus females measuring between 4.5 and 8.0. mm were counted on the right side of each animal. The tick infestation prevalence data were analyzed in relation to the risk factors, using multiple logistic regression. Associations between prevalence and its possible influencing determinants were measured using odds ratios. The average tick count of cows during the peripartum and lactation periods was significantly higher (p<0.05; OR=4.82) than the count in dry cows. Taurine animals showed significantly higher infestation (p<0.05; OR=3.28) than pure zebuine animals. Among cross-bred animals, the infestation was higher (p<0.05) in F1 animals (1/2 taurine × 1/2 zebuine) than in Girolando animals (5/8 zebuine × 3/8 taurine). Primiparous cows (p<0.05) had significantly heavier infestations than multiparous cows. Dairy production showed a positive correlation with tick infestation (p<0.03; OR=2.94), such that the most productive animals were the most parasitized ones. Animals kept at high densities were significantly more heavily infested (p<0.05, OR=6.32) than animals kept at low density. First-lactation and high-productivity taurine animals were more vulnerable to R. microplus, thus comprising a high-risk group in dairy herds. © 2012 Elsevier GmbH.

Short communication: Acute but transient increase in serum insulin reduces messenger RNA expression of hepatic enzymes associated with progesterone catabolism in dairy cows

The objective of this experiment was to evaluate the effects of glucose infusion on serum concentrations of glucose, insulin, and progesterone (P4), as well as mRNA expression of hepatic CYP2C19 and CYP3A4 in nonlactating, ovariectomized cows in adequate nutritional status. Eight Gir × Holstein cows were maintained on a low-quality Brachiaria brizantha pasture with reduced forage availability, but they individually received, on average, 3. kg/cow daily (as fed) of a corn-based concentrate from d -28 to 0 of the experiment. All cows had an intravaginal P4-releasing device inserted on d -14, which remained in cows until the end of the experiment (d 1). On d 0, cows were randomly assigned to receive, in a crossover design containing 2 periods of 24. h each (d 0 and 1), (1) an intravenous glucose infusion (GLUC; 0.5. g of glucose/kg of BW, over a 3-h period) or (2) an intravenous saline infusion (SAL; 0.9%, over a 3-h period). Cows were fasted for 12. h before infusions, and they remained fasted during infusion and sample collections. Blood samples were collected at 0, 3, and 6. h relative to the beginning of infusions. Liver biopsies were performed concurrently with blood collections at 0 and 3. h. After the last blood collection of period 1, cows received concentrate and returned to pasture. Cows gained BW (16.5 ± 3.6. kg) and BCS (0.08 ± 0.06) from d -28 to 0. Cows receiving GLUC had greater serum glucose and insulin concentrations at 3. h compared with SAL cohorts. No treatment effects were detected for serum P4 concentrations, although mRNA expression of CYP2C19 and CYP3A4 after the infusion period was reduced for cows in the GLUC treatment compared with their cohorts in the SAL treatment. In conclusion, hepatic CYP3A4 and CYP2C19 mRNA expression can be promptly modulated by glucose infusion followed by acute increases in circulating insulin, which provides novel insight into the physiological mechanisms associating nutrition and reproductive function in dairy cows. © 2013 American Dairy Science Association.

Effects of vaccination against reproductive diseases on reproductive performance of lactating dairy cows submitted to AI

Four experiments evaluated the effects of vaccination against bovine herpesvirus-1 (BoHV-1), bovine viral diarrhea virus (BVDV), and Leptospira spp. on reproductive performance of lactating dairy cows without (experiments 1, 2, and 3) or with previous vaccination against these diseases (experiment 4). Cows were assigned to a fixed-time AI protocol (FTAI; d -11 to 0) in all experiments, as well as AI 12. h upon estrus detection in experiment 3. Pregnancy status was determined with transrectal ultrasonography on d 30 and 71 (d 60 for experiment 3) after AI. Pregnancy loss was considered in cows pregnant on d 30 but non-pregnant on the subsequent evaluation. In experiment 1, 853 cows received (VAC) or not (CON) vaccination against BoHV-1, BVDV, and Leptospira spp. at the beginning of the FTAI (d -11) and 30. d after AI. Pregnancy loss was reduced (P=0.03) in VAC cows compared with CON. In experiment 2, 287 cows received VAC or CON 30. d prior to (d -41) and at the beginning (d -11) of the FTAI. Pregnancy rates on d 30 and 71 were greater (P≤0.03) in VAC cows compared with CON. In experiment 3, 1680 cows with more than 28. d in milk were randomly assigned to receive VAC or CON with doses administered 14. d apart, and inseminated within 15-135. d after the second dose. Pregnancy rates on d 30 and 60 were greater (P≤0.02) in VAC cows compared with CON. In experiment 4, 820 cows received (REVAC) or not (CON) revaccination against BoHV-1, BVDV, and Leptospira spp. at the beginning of the FTAI protocol (d -11). Pregnancy rates and loss were similar (P≥0.54) between treatments. Hence, vaccinating naïve cows against BoHV-1, BVDV, and Leptospira spp. improved reproductive efficiency in dairy production systems, particularly when both doses were administered prior to AI. © 2013 Elsevier B.V.

Analysis of the risk factors related to the immune humoral anti-Anaplasma marginale in dairy cattle

Anaplasma marginale is the most prevalent pathogen of cattle transmitted by ticks in the world. This study aimed to evaluate the risk factors for anaplasmosis in dairy cattle. Fifty dairy cattle from the herd of Empresa de Pesquisa Agropecuaria do Estado do Rio de Janeiro were selected by proportional stratified sampling. The risk factors evaluated were: physiological state, race pattern, number of lactations, milk production, infestation by Rhipicephalus microplus and animal density. Antibody activity against A. marginale was determined using the indirect enzyme-linked immunosorbent assay. The percentual values of seroprevalence for A. marginale were submitted to X2 test, and the level of minimum significance, to keep a factor in the model of logistic regression, was fixated in 5%. It was observed that pregnancy and lactation influenced significantly (p<0.05) in the seropositivity of the animals. Bos indicus animals had 5.21 times more chances of being seropositive than B. taurus animals. Primiparous female had 88% more chances of being seropositive than pluriparous female. Animals with high milk production were 63% more positive than low production animals. When infested by ticks the animals had 39% more chance of being seropositive to A. marginale. Bos indicus animals presented 5.21 times more chance of being seropositive than B. taurus animals. Primiparous females presented 88% more chance of being seropositive than the pluriparous ones. High milk production animals were 63% more positive than the low production ones. When infested by ticks the animals had 39% more chance of being seropositive to A. marginale. High density grazing provided for the animals 3.2 times more chances of being seropositive than low density ones. The herd was classified as erratic to A. marginale, even being placed in a steady enzootic area.

Hormonal treatment of dairy goats affected by hydrometra associated or not with ovarian follicular cyst

Hydrometra is considered a very important pathological condition, because it represents one of the main causes of temporary infertility in dairy goats. The objective was (i) to evaluate a protocol for the treatment of hydrometra associated (n=2) or not (n=17) with follicular ovarian cyst in 19 dairy goats and (ii) to assess its reproductive efficiency after treatment. For this purpose, 10. mg of dinoprost (PGF) divided in two equal doses were administered to all animals intravulvosubmucosally on Days 0 and 10. In addition, 500. IU hCG were administered on Day 7. Ultrasound exams were performed in all females from Days 0 to 3, 7 and 10 to 13 of treatment, in order to evaluate uterus drainage after each treatment. Goats were monitored for estrus after both treatments and mated after the second dose of PGF. Blood samples were collected from 11 goats to determine plasma progesterone concentrations before, during and after treatment. Of the 19 goats treated, 16 lost weight after the first dose, probably due to uterine discharge. Complete drainage of uterine fluid was observed in 11/19 (57.9%) and 17/19 (89.5%) after the first and second doses, respectively. Afterwards, we diagnosed 2 more goats with follicular cysts, for a total of 21.1% (4/19) of animals exhibiting hydrometra and ovarian cyst concomitantly. In one doe the diagnosis was on Day 2 and in the other on Day 11. All does showed progesterone concentrations superior to 1. ng/mL at Day 0, with an average of 10.6 ± 1.4. ng/mL. Out of the 10 goats mated, only two became pregnant after treatment, corresponding to 10.5% of the total (2/19). Although prostaglandin was effective to drain the uterine fluid and led to the onset of estrus, it did not improve the pregnancy rate. The use of hCG in female goats was not effective in luteinizing the cysts. It can be concluded that hydrometra alone or associated with ovarian follicular cyst may adversely affect goat reproductive performance. © 2012 Elsevier B.V.

Timing of prostaglandin F2α treatment in an estrogen-based protocol for timed artificial insemination or timed embryo transfer in lactating dairy cows

Objectives were to investigate progesterone concentrations and fertility comparing 2 different intervals from PGF2α treatment and induced ovulation in an estrogen-based ovulation synchronization protocol for timed artificial insemination (TAI) or timed embryo transfer (TET) in lactating dairy cows. A total of 1,058 lactating Holstein cows [primiparous (n=371) and multiparous (n=687)], yielding 34.1±0.33 kg of milk/d at various days in milk were randomly assigned to receive treatment with PGF2α on either d 7 or 8 of the following protocol: d 0: 2mg of estradiol benzoate + controlled internal drug release device; d 8: controlled internal drug release device removal + 1.0mg of estradiol cypionate; d 10: TAI or d 17: TET. Only cows with a corpus luteum at d 17 received an embryo and all cows received GnRH at TET. Pregnancy diagnoses were performed by detection (transrectal ultrasonography) of an embryo on d 28 or a fetus on d 60. Fertility [pregnancy per artificial insemination (P/AI) or pregnancy per embryo transfer (P/ET)] was affected by breeding technique (AI vs. ET) and time of PGF2α treatment (d 7 vs. 8) at the 28-d pregnancy diagnosis for TAI [32.9% (238) vs. 20.6% (168)] and TET cows [47% (243) vs. 40.7% (244)] and at the 60-d pregnancy diagnosis for TAI [30% (238) vs. 19.2% (168)] and TET cows [37.9% (243) vs. 33.5% (244)]. The progesterone (P4) concentration at d 10 altered fertility in TAI cows, with higher P/AI in cows with P4 concentration <0.1 ng/mL compared with cows with P4 concentration ≥0.1 ng/mL, and in ET cows, with higher P/ET in cows with P4 concentration <0.22 ng/mL compared with cows with P4 concentration ≥0.22 ng/mL. Prostaglandin F2α treatment at d 7 increased the percentage of cows with P4 <0.1 ng/mL on d 10 [39.4 (85) vs. 23.2 (54)]. Reducing the period between PGF2α and TAI from 72 to 48h in dairy cows resulted in a clear reduction in fertility in cows bred by TAI and a subtle negative effect in cows that received TET. The earlier PGF2α treatment benefits are most likely mediated through gamete transport, fertilization, or early embryo development and a more subtle effect of earlier PGF2α treatment that may be mediated through changes in the uterine or hormonal environment that manifests itself after ET on d 7. © 2013 American Dairy Science Association.

Dairy goat health management and milk production on organic and conventional system in Brazil

The aimed of this article is to measure risk factors on health and milk production on organic and conventional dairy goats in Brazil. Two experimental groups (organic and conventional) were evaluated simultaneously. The study design was completely randomized. The organic herd consisted of 25 goats and 15 kids. In the conventional production system, a dairy herd comprising 40 goats and 20 kids participated in the study. Data on milk production and health management were available from January 2007 to December 2009. The abortion rate in the conventional system was 5% (2/40) whereas in organic system no abortion was diagnosed (0/25). The mortality rate at weaning in the conventional system was 5% (2/40) and in the organic system was 8% (2/25). Milk production was lower (2.20 kg/day) in organic than conventional system (2.66 kg/day). Goats and kids in organic farm had a higher FEC (386±104 and 900±204, respectively) (p<0.05) than those in conventional farm (245±132 and 634±212, respectively). In addition, Saanen kids had higher FEC (p<0.001) than goats. Treatment with antiparasitic drugs was higher in conventional system (50%) than organic system (1.3%).

Detecting Loci under Recent Positive Selection in Dairy and Beef Cattle by Combining Different Genome-Wide Scan Methods

As the methodologies available for the detection of positive selection from genomic data vary in terms of assumptions and execution, weak correlations are expected among them. However, if there is any given signal that is consistently supported across different methodologies, it is strong evidence that the locus has been under past selection. In this paper, a straightforward frequentist approach based on the Stouffer Method to combine P-values across different tests for evidence of recent positive selection in common variations, as well as strategies for extracting biological information from the detected signals, were described and applied to high density single nucleotide polymorphism (SNP) data generated from dairy and beef cattle (taurine and indicine). The ancestral Bovinae allele state of over 440,000 SNP is also reported. Using this combination of methods, highly significant (P<3.17×10-7) population-specific sweeps pointing out to candidate genes and pathways that may be involved in beef and dairy production were identified. The most significant signal was found in the Cornichon homolog 3 gene (CNIH3) in Brown Swiss (P = 3.82×10-12), and may be involved in the regulation of pre-ovulatory luteinizing hormone surge. Other putative pathways under selection are the glucolysis/gluconeogenesis, transcription machinery and chemokine/cytokine activity in Angus; calpain-calpastatin system and ribosome biogenesis in Brown Swiss; and gangliosides deposition in milk fat globules in Gyr. The composite method, combined with the strategies applied to retrieve functional information, may be a useful tool for surveying genome-wide selective sweeps and providing insights in to the source of selection.

Effects of bovine somatotropin injection on serum concentrations of progesterone in non-lactating dairy cows

The objective of this experiment was to evaluate the effects of bovine somatotropin administration on serum concentrations of glucose, insulin, NEFA, IGF-I, and progesterone (P4) in ovariectomized non-lactating dairy cows receiving exogenous P4, as a model to estimate treatment effects on hepatic P4 degradation. Ten non-lactating, non-pregnant, and ovariectomized Gir×Holstein cows were assigned to the experiment (d -14 to 27). On d 0, cows were ranked by BW and BCS, and randomly assigned to one of two treatments: (1) bovine somatotropin (BST; n=5) or (2) saline control (control; n=5). Cows assigned to the BST treatment were administered s.c. injections containing 500. mg of sometribove zinc on d 0, 9, and 18 of the experiment, whereas control cows concurrently received a 10-mL s.c. injection of 0.9% saline. On d -2, cows were inserted with an intravaginal releasing device containing 1.9. g of P4, which remained in the cows until the end the experiment (d 27). Cow BW and BCS were assessed on d -14, 0, and 27. Blood samples were collected daily from d 0 to d 27, at 0 (immediately before), 1, and 2. h relative to concentrate feeding for determination of serum glucose, insulin, NEFA, P4, and IGF-I concentrations. Concentrations of glucose, NEFA, and insulin obtained prior to feeding (0. h) were used to determine pre-prandial revised quantitative insulin sensitivity check index (RQUICKI). No treatment effects were detected for BW (P=0.72) and BCS change (P=0.79) during the experiment. Beginning on d 2 of the experiment, BST cows had greater (P≤0.01) serum IGF-I concentrations compared with control cohorts (treatment×day interaction; P<0.01). Cows receiving BST had greater (P≤0.05) insulin concentrations compared with control cohorts from d 8 to d 11, d 16 and 17, as well as from d 19 to d 21 of the experiment (treatment×day interaction; P<0.01). Cows receiving BST had greater (P≤0.01) mean glucose and NEFA concentrations, as well as reduced (P<0.01) mean RQUICKI during the experiment compared with control cohorts. No treatment effects, however, were detected (P=0.73) for serum P4 concentrations. In conclusion, results from this experiment indicate that hepatic P4 catabolism is not directly regulated by circulating IGF-I, whereas BST administration decreases insulin sensitivity in non-lactating dairy cows in adequate nutritional status. © 2013 Elsevier B.V.

Effects of BST and high energy diet on gene expression in mammary parenchyma of dairy heifers

The objective of this study was to determine the effects of dietary energy and recombinant bovine somatotropin (bST) injection to identify genes that might control mammogenesis. Total RNA was extracted from the parenchymal tissue of 32 heifers randomly assigned to one of four treatments: two diets (a standard diet and a high energy, high protein diet), each with or without bST. To perform microarray experiments, RNA samples were pooled (2 animals/pool) before reverse transcription and labeling with Cy3 or Cy5. A 4-node loop design was used to examine the differential gene expression among treatments using a bovine-specific cDNA micro array (National Bovine Functional Genomics Consortium Library, NBFGC) containing 18,263 unique expressed sequence tags (EST). Significance levels of differential gene expression among treatments were assessed using a mixed model approach. Injection of bST altered the expression of 12 % of the genes on NBFGC slide related to tissue development, whereas 6% were altered by diet. Administration of bST increases the expression of genes positively related to cell proliferation and mammary parenchyma to a greater extent than a high energy diet. © 2013 Sociedade Brasileira de Zootecnia.

Effect of temperature on the lipolytic and proteolytic activity of Bacillus cereus isolated from dairy products

Bacillus cereus is a bacterium with deteriorating potential for dairy products, by being a psychrotrophic organism producer of lipases and proteases. This study evaluated the psychrotrophic behavior, lipolytic and proteolytic activity at 30°C, 10°C and 7°C of 86 strains of B. cereus lato sensu isolated from dairy products, marketed in Southern Brazil. It was also evaluated the optimal temperature for protease production. No strain grew at 7°C; but at 10°C, 84.9% of strains have grown. Only one strain had lipolytic activity at 30°C, and none at 7°C. At 10°C, 16.3% of strains produced lipases. All the strains presented proteolytic activity at 30°C; and at 10°C, 72.1% had this activity, and at 7°C, only 4.6%, an amount significantly lower (p < 0.05). The temperature of 20°C promoted the highest proteolytic activity, and at 10°C, the lowest activity. B. cereus can produce lipases and proteases at room and marginal chilling temperatures, causing technological defects in dairy products stored under these conditions. © 2008 IFRJ.

Effects of a single nucleotide polymorphism in the leptin gene on the productive traits of dairy buffaloes (Bubalus bubalis)

The gene responsible for coding the leptin hormone has been associated with productive and reproductive traits in cattle. In dairy cattle, different polymorphisms found in the leptin gene have been associated with several traits of economic interest, such as energy balance, milk yield and composition, live weight, fertility and dry matter consumption. The aim of this study was to detect genetic variability in the leptin gene of buffaloes and to test possible associations with milk yield, fat and protein percentages, age at first calving and first calving interval. Three genotypes (AA, AG and GG) were identified by polymerase chain reaction-restriction fragment length polymorphism, which presented genotypic frequencies of 0.30, 0.54 and 0.16, respectively. The allele frequencies were 0.57 for the A allele and 0.43 for the G allele. No significant effects were found in the present study, but there is an indicative that leptin gene affects lipid metabolism. © 2013 Springer Science+Business Media Dordrecht.

Nutrients balances and milk fatty acid profile of mid lactation dairy cows supplemented with unsaturated fatty acid

The objective was to evaluate the effect of unsaturated fatty acid sources supplementation on nutrients balances and milk fatty acid profile of mid lactation dairy cows. Twelve Brazilian Holstein cows in the mid lactation (mean of 128 days) and (580 ± 20kg of weight; mean ± SD) with milk yield of 25kg/d were assigned randomly into three 4 × 4 Latin square, fed the following diets: control (C); refined soybean oil; (SO); whole soybean raw (WS) and; calcium salts of unsaturated fatty acids (CSFA). Milk yield was 26.6; 26.4; 24.1 and 25.7 to the diets CO, SO, WS and CSFA respectively. Cows fed the WS treatment produced less milk (1.95kg/d of milk), fat and lactose than did cows fed the SO and CSFA. Cows fed the CSFA treatment showed less blood, urine (g/d) concentrations of N more energetic efficiency and intake of energy than did cows fed the SO treatment. Cows fed the unsaturated fatty acids sources showed more C18:2 cis-9, trans-11 CLA and trans-C18:1 FA concentration in milk than did cows fed the CO treatment. Diets with whole soybeans and soybeans oil provide more efficient digestive processes, and increase milk composition of unsaturated fatty acids.

Timed artificial insemination programs during the summer in lactating dairy cows: Comparison of the 5-d Cosynch protocol with an estrogen/progesterone-based protocol

The objective of this study was to compare a GnRH-based to an estrogen/progesterone (E2/P4)-based protocol for estrous cycle synchronization and fixed timed artificial insemination (TAI), both designed for synchronization of ovulation and to reduce the period from follicular emergence until ovulation in cows with a synchronized follicular wave. A total of 1,190 lactating Holstein cows (primiparous: n = 685 and multiparous: n = 505) yielding 26.5 ± 0.30 kg of milk/d at 177 ± 5.02 d in milk were randomly assigned to receive one of the following programs: 5-d Cosynch protocol [d -8: controlled internal drug release (CIDR) + GnRH; d -3: CIDR removal + PGF2α; d -2: PGF2α; d 0: TAI + GnRH] or E2/P4 protocol (d -10: CIDR + estradiol benzoate; d -3: PGF2α; d -2: CIDR removal + estradiol cypionate; d 0: TAI). Rectal temperature and circulating progesterone (P4) were measured on d -3, -2, 0 (TAI), and 7. The estrous cycle was considered to be synchronized when P4 was ≥1.0 ng/mL on d 7 in cows that had luteolysis (P4 ≤0.4 ng/mL on d 0). To evaluate the effects of heat stress, cows were classified by number of heat stress events: 0, 1, and 2-or-more measurements of elevated body temperature (≥39.1°C). Pregnancy success (pregnancy per artificial insemination, P/AI) was determined at d 32 and 60 after TAI. The cows in the 5-d Cosynch protocol had increased circulating P4 at the time of PGF2α injection (2.66 ± 0.13 vs. 1.66 ± 0.13 ng/mL). The cows in the E2/P4 protocol were more likely to be detected in estrus (62.8 vs. 43.4%) compared with the cows in the 5-d Cosynch protocol, and expression of estrus improved P/AI in both treatments. The cows in the 5-d Cosynch protocol had greater percentage of synchronized estrous cycle (78.2%), compared with cows in the E2/P4 protocol (70.7%). On d 60, the E2/P4 protocol tended to improve P/AI (20.7 vs. 16.7%) and reduced pregnancy loss from 32 to 60 d (11.0 vs. 19.6%), compared with the 5-d Cosynch protocol. In cows with their estrous cycle synchronized, the E2/P4 protocol had greater P/AI (25.6 vs. 17.7%) on d 60 and lower pregnancy loss from 32 to 60 d (6.7 vs. 21.7%) compared with cows in the 5-d Cosynch protocol. Follicle diameter affected pregnancy loss from 32 to 60 d only in the cows in the 5-d Cosynch protocol, with smaller follicles resulting in greater pregnancy loss. Pregnancy per AI at d 60 was different between protocols in the cows with 2 or more measurements of heat stress (5-d Cosynch = 12.2% vs. E2/P4 = 22.8%), but not in the cows without or with 1 heat stress measurement. In conclusion, the 5-d Cosynch protocol apparently produced better estrous cycle synchronization than the E2/P4 protocol but did not improve P/AI. The potential explanation for these results is that increased E2 concentrations during the periovulatory period can improve pregnancy success and pregnancy maintenance, and this effect appears to be greatest in heat-stressed cows when circulating E2 may be reduced. © 2013 American Dairy Science Association.