Advertisement
Research Article| Volume 96, ISSUE 3, P1701-1714, March 2013

Feeding canola meal to dairy cows: A meta-analysis on lactational responses

Open ArchivePublished:January 07, 2013DOI:https://doi.org/10.3168/jds.2012-5740

      Abstract

      The objective of this meta-analysis was to determine the effects of the substitution of a protein source by canola meal (CM) on lactational responses (CM minus control) in dairy cows. The study included 49 comparisons of isonitrogenous (±1.0% dietary CP content) dietary treatments published since 1975 (27 experiments). The CM intake ranged from 1.0 to 4.0 kg/d (standard deviation = 0.65). Regressions were forced through the origin, weighted by sample size, and controlled for changes in dry matter intake and in dietary concentrations of CP and ether extract. Milk yield and milk protein yield responded positively to the substitution of a protein source by CM. The response in milk protein yield was affected by the type of protein source that was substituted; the positive response was half less when CM was substituted for soybean meal compared with substitution of CM for other protein sources. The latter effect was in part related to a positive response on milk protein percentage when CM replaced protein sources other than soybean meal. The response in efficiency of N utilization (milk N yield/N intake) was positive to the substitution of a protein source by CM. Negative changes in supply of metabolizable protein (MP) estimated from the 2001 National Research Council model were associated with positive responses in milk protein yield with CM substitution, a finding contrary to the expected positive relationship between supply of MP and milk protein yield. In conclusion, a protein supplement can be substituted by CM with positive effects on milk and milk protein production. These data also indicate an underestimation of MP supply associated with CM inclusion in dairy rations using the National Research Council model.

      Key words

      Introduction

      Cruciferous seeds and plants contain glucosinolates, which can yield breakdown products potentially goitrogenic, hepatotoxic, or pungent (
      • Fenwick G.R.
      The assessment of a new protein source—Rapeseed.
      ). Goitrogenic compounds reduce the availability of iodine to the animal and diminish the synthesis of thyroxine, which is involved in the hormonal mechanisms of milk production (
      • Grenet N.
      • Journet M.
      Rapeseed oil meal in animal feeding. III. Influence of the processing method and the proportion of oil meal in concentrate fed to dairy cows on feed intake, milk production and composition.
      ). A lower grain intake was also reported with diets containing rapeseed meal high in glucosinolate, tannin, or phytate content (
      • Papas A.
      • Ingalls J.R.
      • Cansfield P.
      Effects of Tower and 1821 rapeseed meals and Tower gums on milk yield, milk composition and blood parameters of lactating dairy cows.
      ;
      • Thomke S.
      Review of rapeseed meal in animal nutrition: Ruminant animals.
      ). Therefore, this type of rapeseed meal was not considered a suitable feedstuff for livestock (
      • Abbadi A.
      • Leckband G.
      Rapeseed breeding for oil content, quality, and sustainability.
      ).
      In 1969, a low-glucosinolate trait was identified in the Polish rapeseed spring variety Bronowski and that discovery started an international backcrossing breeding program to introduce the new trait into high-yielding erucic acid-free rapeseed varieties (
      • Abbadi A.
      • Leckband G.
      Rapeseed breeding for oil content, quality, and sustainability.
      ). A few years later, Canadian plant breeders were able to license the first rapeseed cultivar with low levels of both glucosinolates and erucic acid (Tower;
      • Stefansson B.R.
      • Kondra Z.P.
      Tower summer rape.
      ). The term canola was coined to specifically identify rapeseed varieties containing less than 2% erucic acid in the oil portion and less than 30 µmol of glucosinolates/g in the meal portion (i.e., double-zero or double-low rapeseed;
      • Bell J.M.
      Nutrients and toxicants in rapeseed meal: A review.
      ;
      • Newkirk R.
      ).
      Since the mid 1980s, Canadian canola producers have consistently grown more canola in a greater area than ever before (
      • Harker K.N.
      • O’Donovan J.T.
      • Turkington T.K.
      • Blackshaw R.E.
      • Lupwayi N.Z.
      • Smith E.G.
      • Klein-Gebbinck H.
      • Dosdall L.M.
      • Hall L.M.
      • Willenborg C.J.
      • Kutcher H.R.
      • Malhi S.S.
      • Vera C.L.
      • Gan Y.
      • Lafond G.P.
      • May W.E.
      • Grant C.A.
      • McLaren D.L.
      High-yield no-till canola production on the Canadian prairies.
      ). Consequently, the use of its coproduct, canola meal (CM), also increased tremendously as a ruminant feed in Canada (
      • Christensen D.A.
      • McKinnon P.J.
      Canola meal for beef and dairy cattle.
      ) as well as in other countries (

      Emanuelson, M. 1989. Rapeseed Products of Double Low Cultivars to Dairy Cows: Effects of Long-Term Feeding and Studies on Rumen Metabolism. Volume 189 of Rapport. Swedish University of Agricultural Sciences, Department Animal Nutrition Management, Uppsala, Sweden.

      ;
      • Tuori M.
      Rapeseed meal as a supplementary protein for dairy cows on grass silage-based diet, with the emphasis on the Nordic AAT-PBV feed protein evaluation system.
      ;
      • Huhtanen P.
      Supply of nutrients and productive responses in dairy cows given diets based on restrictively fermented silage.
      ;

      Moss, A. 2002. The effects of long-term feeding of extracted rapeseed meal and whole rapeseed on the physical and financial performance, health and welfare of high yielding dairy cows. HGCA Project Report OS59. Home Grown Cereals Authority, London, UK.

      ;
      • Tan S.H.
      • Mailer R.J.
      • Blanchard C.L.
      • Agboola S.O.
      Canola proteins for human consumption: Extraction, profile, and functional properties.
      ). In Western Canada and parts of the United States, CM is the principal source of protein for dairy cow diets because it is readily available and provides a high-quality protein (
      • Hickling D.
      Maximized utilization of canola co-products in livestock industry.
      ;
      • Mulrooney C.N.
      • Schingoethe D.J.
      • Kalscheur K.F.
      • Hippen A.R.
      Canola meal replacing distillers grains with solubles for lactating dairy cows.
      ).
      Over the years, many manuscripts have been published comparing lactational performances of dairy cattle obtained with CM and other protein supplements, especially soybean meal (SBM). In an extensive review,
      • Hill R.
      Rapeseed meal in the diets of ruminants.
      concluded that rapeseed meal low in glucosinolates could be used as freely as SBM and that milk production and composition were equally satisfactory from either protein supplement. More recently,
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      provided evidence that milk production was as good or better and that milk protein secretion was higher with CM compared with SBM in grass silage-based diets. Positive responses in milk protein yield might be related to an improved efficiency of microbial synthesis and reduced ruminal CP degradability, based on the findings of a meta-analysis on omasal studies (
      • Broderick G.A.
      • Huhtanen P.
      • Ahvenjärvi S.
      • Reynal S.M.
      • Shingfield K.J.
      Quantifying ruminal nitrogen metabolism using the omasal sampling technique in cattle—A meta-analysis.
      ). As suggested by
      NRC
      , the estimated supply of MP, either from RUP or microbial protein, is the main factor that determines milk protein secretion in lactating dairy cows.
      Therefore, the first objective of this meta-analysis was to evaluate lactational performances when CM substituted another protein source in dairy rations. The second objective was to determine if lactational responses were affected by experimental conditions or factors such as the type of forage in the diet or the protein source that was substituted. A third objective was to evaluate if changes in milk protein yield were in line with changes in estimated supply of MP (
      NRC
      ).

      Materials and Methods

      Database

      A comprehensive literature review was conducted in October 2011 and updated in February 2012 to track experiments in which CM (or a low-glucosinolate rapeseed meal) was replaced by another protein supplement in dairy cows. A large body of literature was compiled after a thorough search of computerized-stored databases [e.g., Scopus (http://www.scopus.com/home.url), PubMed (http://www.ncbi.nlm.nih.gov/pubmed), and Canola Council database (http://www.canolacouncil.org/research/)]. References in articles and review manuscripts were scrutinized in the search for new articles. The literature search was not limited to studies published in North American journals or to publications in the English language.
      Four a priori criteria were used to select eligible studies: (1) N intake reported, (2) isonitrogenous diets [i.e. difference in dietary CP concentration not exceeding 1.0% (DM basis) between the control and the CM diets], (3) no CM in the control diet and a proportion of CM of 0.05 (DM basis) or more in the treated diet, and (4) treated diets with a high proportion of CM in the protein supplement portion of the diet. For the latter criterion, a protein supplement was defined as a concentrate feed with a CP concentration of 20% or more (DM basis) and a high proportion of CM was defined as 0.85 or more (weight basis) or 0.75 or more (N basis). The breakpoint was set at 0.85 to ensure a minimal bias of responses by the complementary protein source in the treated diets.
      Among the eligible studies, some were deemed outliers and excluded from the database for the reasons outlined below. The heat-damaged and dark-colored CM used in the study of
      • Fisher L.J.
      • Walsh D.S.
      Substitution of rapeseed meal for soybean meal as a source of protein for lactating cows.
      was excluded because it depressed the apparent DM digestibility and the lactational performances. Initial milk production was low (i.e., 13.5 kg/cow per day) in the study of
      • Khorasani G.R.
      • Robinson P.H.
      • Kennelly J.J.
      Evaluation of solvent and expeller linseed meals as protein sources for dairy cattle.
      and diets had a low CP concentration (i.e., 12.5%) in the study of
      • Munger A.
      Use of rapeseed oilmeal concentrates in dairy cow feeding.
      compared with the other studies in the database; therefore, these 2 studies were also excluded from the eligible database. The study of
      • Mancuso W.A.
      • Santini F.J.
      • Gagliostro G.A.
      • Rearte D.H.
      Supplementacíon con harina de girasol o de colza “00” a vacas lecheras alimentadas con dietas en base a silaje de maíz.
      had a large gap between its dietary proportion of CM (i.e., 30%) and the upper proportion of CM in the database (i.e., 17.2%; Table 1); therefore, that study was not included in the database. Finally, the urea-based diet (1.9% DM) in
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      reduced DMI and lactational performances compared with other diets based on true protein sources; that urea-based diet was also excluded from the database.
      Table 1Descriptive statistics of intakes of protein source and canola meal (CM); dietary proportion of CM, CP, and ether extract (EE) in control (Ctrl) and CM diets; and the difference between diets (Δ; CM minus Ctrl)
      Δ = CM diet mean minus Ctrl diet mean. A positive value represents an increase in response to CM and a negative value a decrease.
      Variable
      Protein source was defined as a concentrated feed containing at least 20% CP (DM basis).
      Dietn
      Overall, the database included 88 different treatment means that generated 49 pairs of Ctrl and CM.
      MeanSDMinimumMaximum
      Protein source (kg/cow per day)Ctrl442.00.930.75.0
      CM442.30.651.04.0
      Δ49+0.430.599−1.801.10
      CM intake (kg/cow per day)Ctrl440.00.000.00.0
      CM442.30.651.04.0
      Δ49+2.400.7131.044.01
      CM proportion (100 g/kg, DM basis)Ctrl440.00.000.00.0
      CM441.170.2670.511.72
      Δ49+1.170.2670.511.72
      CP (%, DM basis)Ctrl4416.61.5614.020.3
      CM4416.41.5713.620.1
      Δ490.000.371−0.751.00
      EE (%, DM basis)Ctrl443.700.9721.965.72
      CM443.751.0192.096.11
      Δ49+0.190.364−0.801.07
      1 Δ = CM diet mean minus Ctrl diet mean. A positive value represents an increase in response to CM and a negative value a decrease.
      2 Protein source was defined as a concentrated feed containing at least 20% CP (DM basis).
      3 Overall, the database included 88 different treatment means that generated 49 pairs of Ctrl and CM.
      With the exception of

      Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

      , authors of all eligible manuscripts were contacted to confirm the use of CM before inclusion of the experiment in the database if the glucosinolate concentration of the rapeseed meal was not specified in the article. The final database contained 27 experiments (88 diets) reported in 22 publications using these exclusion criteria:
      • Ingalls J.R.
      • Sharma H.R.
      Feeding of Bronowski, Span and commercial rapeseed meals with or without addition of molasses or flavor in the rations of lactating cows.
      ,
      • Laarveld B.
      • Christensen D.A.
      Rapeseed meal in complete feeds for dairy cows.
      ,
      • Sharma H.R.
      • Ingalls J.R.
      • McKirdy J.A.
      Effects of feeding a high level of Tower rapeseed meal in dairy rations on feed intake and milk production.
      ,
      • Papas A.
      • Ingalls J.R.
      • Cansfield P.
      Effects of Tower and 1821 rapeseed meals and Tower gums on milk yield, milk composition and blood parameters of lactating dairy cows.
      ,
      • DePeters E.J.
      • Bath D.L.
      Canola meal versus cottonseed meal as a protein supplement in dairy rations.
      ,
      • Robinson P.H.
      • Kennelly J.J.
      Influence of intake of rumen undegradable protein on milk production of late lactation Holstein cows.
      ,
      • Huhtanen P.
      • Khalili H.
      • Näsi M.
      A comparison of untreated and formaldehyde-treated barley distiller's solubles and rapeseed meal as protein supplements in dairy cows given grass silage ad libitum.
      ,
      • Macleod G.K.
      Canola meal as a protein supplement in corn based dairy rations.
      ,
      • McClean C.
      • Laarveld B.
      Effect of somatotropin and protein supplement on thyroid function of dairy cattle.
      ,
      • Tuori M.
      Rapeseed meal as a supplementary protein for dairy cows on grass silage-based diet, with the emphasis on the Nordic AAT-PBV feed protein evaluation system.
      ,
      • Emanuelson M.
      • Ahlin K.-Å.
      • Wiktorsson H.
      Long-term feeding of rapeseed meal and full-fat rapeseed of double low cultivars to dairy cows.
      ,
      • Huhtanen P.
      • Jaakkola S.
      • Saarisalo E.
      The effects of concentrate energy source on the milk production of dairy cows given a grass silage-based diet.
      ,
      • Piepenbrink M.S.
      • Schingoethe D.J.
      • Brouk M.J.
      • Stegeman G.A.
      Systems to evaluate the protein quality of diets fed to lactating cows.
      ,
      • Dewhurst R.J.
      • Aston K.
      • Fisher W.J.
      • Evans R.T.
      • Dhanoa M.S.
      • McAllan A.B.
      Comparison of energy and protein sources offered at low levels in grass-silage-based diets for dairy cows.
      ,
      • Khalili H.
      • Sairanen A.
      • Hissa K.
      • Huhtanen P.
      Effects of type and treatment of grain and protein source on dairy cow performance.
      ,
      • Shingfield K.J.
      • Jaakkola S.
      • Huhtanen P.
      Effects of level of nitrogen fertilizer application and various nitrogenous supplements on milk production and nitrogen utilization of dairy cows given grass silage-based diets.
      ,
      • Vanhatalo A.
      • Pahkala E.
      • Salo-Väänänen P.
      • Korhonen H.
      • Piironen V.
      • Huhtanen P.
      Rapeseed and soybean as protein supplements of dairy cows fed grass silage based diets.
      ,

      Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

      ,
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      ,
      • Mulrooney C.N.
      • Schingoethe D.J.
      • Kalscheur K.F.
      • Hippen A.R.
      Canola meal replacing distillers grains with solubles for lactating dairy cows.
      ,
      • Christen K.A.
      • Schingoethe D.J.
      • Kalscheur K.F.
      • Hippen A.R.
      • Karges K.K.
      • Gibson M.L.
      Response of lactating dairy cows to high protein distillers grains or 3 other protein supplements.
      , and
      • Oba M.
      • Penner G.B.
      • Whyte T.D.
      • Wierenga K.
      Effects of feeding triticale dried distillers grains plus solubles as a nitrogen source on productivity of lactating dairy cows.
      .

      Rations and Nitrogen Fractions of Canola Meal

      Dietary characteristics of rations were estimated using the
      NRC
      model with feed ingredients as similar as possible to those reported in the publications. The chemical composition of each feed ingredient was used when reported [e.g., CP, NDF, ADF, ether extract (EE), ash, and N fractions]; otherwise, table values of
      NRC
      were used.
      Nitrogen fractions of CM were required to estimate MP supply. The oil extraction process affects the N fractions and the rumen CP degradability of CM. The processing method of CM, either mechanical or solvent extraction (
      • Khorasani G.R.
      • Robinson P.H.
      • Kennelly J.J.
      Evaluation of solvent and expeller linseed meals as protein sources for dairy cattle.
      ), was seldom reported in the publications. When missing, the processing method of CM was assumed to be mechanical extraction with no heating, unless the oil content was reported and was less than 4%, in which case CM was presumed to be solvent extracted (
      • Lusas E.E.
      • Riaz M.N.
      Fats and oils in feedstuffs and pet foods.
      ). Nitrogen fractions of CM or rapeseed meal from in situ methodology and published in
      NRC
      were used to estimate rumen CP degradability when not reported in the manuscript. Data reported in Table 15-1 of
      NRC
      for the unheated mechanically extracted CM and data reported in Table 15-2b of
      NRC
      for the heated mechanically extracted CM and the solvent-extracted rapeseed meal (unheated or heated) were used. It should be pointed out that data reported in Table 15-2b of
      NRC
      are for less commonly used feedstuffs and, therefore, means are less reliable because they are derived from a fewer number of observations.

      Meta-Analyses

      Coding of Data

      This meta-analysis was conducted on responses (Δ: the value from the CM diet minus the value from the control diet). A publication could contain more than 1 feeding trial or experiment and each experiment could contrast more than 1 CM diet to 1 or more control diets. For example,
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      compared CM with solvent-extracted SBM and cottonseed meal (the urea diet was excluded based on previous criteria); therefore, 2 isonitrogenous pairs of diets were constituted by coupling SBM and cottonseed meal with CM. As previously mentioned, a pair of diets was deemed isonitrogenous as long as the difference in dietary CP content did not exceed 1.0% (DM basis). From the 88 diets in the final database, 49 pairs of control and treated diets were generated.

      Responses

      Responses included DMI in kg/cow per day; yields of milk, 4% FCM (

      Gaines, W. L. 1928. The energy basis of measuring milk yield in dairy cows. Illinois Agric. Expt. Sta., Bull. No. 308:401–438. University of Illinois Agricultural Experiment Station, Urbana.

      ), and ECM (
      • Sjaunja L.O.
      • Baevre L.
      • Junkkarinen L.
      • Pedersen J.
      • Setälä J.
      A Nordic proposal for an energy corrected milk (ECM) formula.
      ) in kilograms per cow per day; percentages of milk protein, fat, and lactose; yields of milk protein, fat, and lactose in grams per cow per day; feed efficiency (ECM/DMI, kg/kg); and apparent N efficiency (N in milk/N intake, g/kg). Yields of milk protein, fat, and lactose, as well as feed and apparent N efficiencies were calculated from milk production, milk percentages, and intake data. Production data in
      • Khalili H.
      • Sairanen A.
      • Hissa K.
      • Huhtanen P.
      Effects of type and treatment of grain and protein source on dairy cow performance.
      were kindly provided by P. Huhtanen (Swedish University of Agricultural Sciences, Umeå, Sweden, personal communication).

      Independent Covariates

      The explanatory covariate of interest was the dietary proportion of CM and was forced in each model. The dietary proportion of CM was expressed in 100 g/kg (DM basis); therefore, the coefficients of CM presented in the current study represent the response expected for each increment of 1 unit of 10% inclusion of CM in the ration (e.g., 2 kg of CM per 20 kg of DMI). Other units of CM were explored to report the results [e.g., intake of CM (kg/d, DM basis), intake of N from CM (g/d), and proportion of CM in the concentrate]. Given that the CP content of CM was reported only in half the publications, N intake from CM could not be used as the unit of CM. Furthermore, expressing the unit of CM as the proportion of CM in the diet (DM basis) appeared to be the best choice because DMI and the amount of CM ingested were reported in all experiments and the ratio took into account the differences in DMI that existed among trials.
      Besides CM, the independent covariates also included differences in DMI (ΔDMI) and in dietary CP concentration (ΔCP) because both covariates were reported in all experiments and could act as potential confounders. Differences in dietary EE concentration (ΔEE) were also considered as a covariate. These were reported in 60% of experiments, whereas ΔEE were estimated using the
      NRC
      model in the other 40%. To determine if ΔEE should be added or not to the models of the whole database, ΔEE was first tested in the subset database of experiments that reported ΔEE (i.e., 60% of experiments). If ΔEE was a significant covariate (i.e., P ≤ 0.10), then it was added to the models of the whole database; if not significant, then the models in the whole database only included the following independent covariates: CM, ΔDMI, and ΔCP. By controlling for ΔDMI, ΔCP, and ΔEE, the responses to CM substitution were assumed to be at equal DMI and equal dietary CP and EE concentrations.

      Factors

      Each response was analyzed in relation to experimental conditions or factors that could explain some within-experiment variability (e.g., type of protein source substituted, main type of forage in the diet, and origin of the study). The type of protein source substituted was encoded into 2 classes: SBM as a single source (n = 16 studies or 23 experiments) or other protein sources (e.g., corn gluten meal, cottonseed meal, distillers product, or a mixture of more than 1 protein source including SBM; n = 13 studies or 26 experiments). The main type of forage was categorized into 2 classes: grass and (or) legume (n = 16 studies or 30 experiments) or other forages (e.g., corn or barley silage, alone or mixed with a similar proportion of grass or legume forage; n = 11 studies or 19 experiments). The origin of publication was encoded into 2 classes: North America (n = 17 studies or 27 experiments) or Europe (n = 10 studies or 22 experiments). Other factors were analyzed but are not reported in this paper because they had minor influence on the responses (e.g., stage of lactation, CP and concentrate levels in the control diet, initial milk production level, and initial milk protein and fat concentrations).

      Statistical Model

      Regressions between the responses and the amount of inclusion of CM were forced through the origin (i.e., no intercept in the model) because responses were assumed to be zero when CM was zero (
      • Glasser F.
      • Ferlay A.
      • Chilliard Y.
      Oilseed lipid supplements and fatty acid composition of cow milk: A meta-analysis.
      ). Indeed, coefficients of CM are similar whether obtained from a regression of response values or from a within-experiment analysis on absolute values (Figure 1A and B). Only linear relationships could be determined in the within-experiment analysis because only 1 study had more than 2 data points (i.e.,
      • Robinson P.H.
      • Kennelly J.J.
      Influence of intake of rumen undegradable protein on milk production of late lactation Holstein cows.
      ).
      Figure thumbnail gr1
      Figure 1Positive within-experiment relationship between milk yield (kg/cow per day) and the dietary proportion of canola meal (CM; 100 g/kg, DM basis) based on absolute values (Milk) or based on responses (Δ Milk; CM diet minus control diet, where each point represents a pair of joined data from panel A). P-value in superscript: panel A: Milk = +0.62≤0.001 (±0.083) × CM + 28.8≤0.001 (±0.07); n = 98; panel B: Δ Milk = +0.62≤0.001 (±0.084) × CM; n = 49.
      All regressions were weighted by the number of cows in the CM diet. The covariates ΔDMI, ΔCP, or ΔEE were excluded from the model if not significant (P > 0.10) or if the variance inflation factor indicated multicollinearity. Furthermore, significant covariates other than CM were deleted from the model as long as they did not modify meaningfully the coefficient of CM (
      • Kleinbaum D.G.
      • Kupper L.L.
      • Muller K.E.
      • Nizam A.
      Confounding and interaction in regression.
      ). Outliers were identified by evaluating their influence on Cook’s distance (
      • Cook R.D.
      Influential observations in linear regression.
      ) and were removed if they changed the relationship. Graphical examinations were used at each stage of the meta-analysis to detect extreme data.
      Each experimental condition or factor was tested separately. A categorical factor was deemed significant if the P-value of the interaction between the factor and CM was ≤0.10; upon significance, the database was divided by the classes of the factor and the response was retested in each subset database as described previously. The independence of distribution of CM diets across the classes of 2 factors was tested with the Pearson chi-square test. Residual values (i.e., differences between model-predicted and measured value of the studied parameter) were evaluated graphically to assess the extent to which each model accounted for the variation in the observed data. Because the experimental conditions were not randomly distributed and were specific to each study, variables were considered as fixed factors; therefore, the range of inference for this meta-analysis is limited to the domain of the specific experiments in the data set (
      • Sauvant D.
      • Schmidely P.
      • Daudin J.J.
      • St-Pierre N.R.
      Meta-analyses of experimental data: Application in animal nutrition.
      ). All analyses were conducted with Minitab (version 16.1.1.0, 2010; Minitab Inc., State College, PA). For comparison purpose, mixed models of SAS (

      Littell, R. C., G. A. Milliken, W. W. Stroup, and R. D. Wolfinger. 1996. SAS® system for mixed models. SAS Institute, Cary, NC.

      ) with study as a random effect were conducted and root mean squared errors (RMSE) were adjusted for random study effects as described by
      • St-Pierre N.R.
      Invited review: Integrating quantitative findings from multiple studies using mixed model methodology.
      .
      The project was submitted to and approved by the Institutional Committee for Animal Care of the Sherbrooke Research Centre (Quebec, Canada), although no animals were required for this study.

      Results

      Meta-Design

      The composition of CM was reported in half the publications and was similar to that of the mechanically extracted CM in
      NRC
      : CP = 39.4 ± 3.65% (range: 32.3 to 44.1%; n = 14); NDF = 27.7 ± 5.85% (range: 19.9 to 37.1%; n = 9); ADF = 18.2 ± 2.39% (range: 12.9 to 21.0%; n = 13); and EE = 5.0 ± 3.40% (range: 1.1 to 11.7%; n = 12).
      Descriptive statistics of diet composition and of lactational performances of cows fed control diets (no CM) and CM diets, and the resulting responses to the substitution of a protein source by CM are listed in Tables 1 and 2, respectively. The number of control and CM diets is lower than the number of differences or responses because control and CM diets included only unrepeated diets (refer to coding of data). The within-experiment relationships in all experiments and among significant factors are presented in Table 3.
      Table 2Descriptive statistics of lactational performance in control (Ctrl) and canola meal (CM) diets, and the resulting responses to substitution of a protein source by CM (Δ; CM minus Ctrl)
      Δ is the change in response of each dependent variable to the substitution of a protein source by CM. A positive value represents an increase in response to CM and a negative value a decrease.
      Variable
      Dry matter intake, MY = milk yield, 4% FCM (Gaines, 1928), and ECM (Sjaunja et al., 1990) in kilograms per cow per day; protein, fat, and lactose = percentages of milk protein, fat, and lactose; MPY, MFY, and MLY = yields of milk protein, fat, and lactose; ECM/DMI = efficiency of milk production; milk N/N intake = apparent N efficiency.
      Dietn
      Overall, the database included 88 individual treatment means that generated 49 pairs of Ctrl and CM (for milk lactose % and yield: 62 treatment means and 34 pairs).
      MeanSDMinimumMaximum
      DMI (kg/cow per day)Ctrl4420.32.8415.226.0
      CM4420.12.6214.725.7
      Δ49+0.290.508−1.11.6
      MY (kg/cow per day)Ctrl4429.36.2018.440.5
      CM4428.96.3619.141.4
      Δ49+0.730.774−1.02.3
      4% FCM (kg/cow per day)Ctrl4428.65.3319.736.9
      CM4428.05.8019.338.5
      Δ49+0.520.848−1.32.6
      ECM (kg/cow per day)Ctrl4428.55.3419.536.9
      CM4427.95.7019.537.6
      Δ49+0.520.798−1.22.3
      Protein (%)Ctrl443.190.1882.843.58
      CM443.230.1852.833.51
      Δ49+0.0160.0659−0.110.15
      MPY (g/cow per day)Ctrl44927175.75931,260
      CM44926178.86391,282
      Δ49+28.728.03−2897
      Fat (%)Ctrl443.900.5502.845.08
      CM443.830.5262.424.89
      Δ49−0.0480.154−0.420.31
      MFY (g/cow per day)Ctrl441,126214.37601,422
      CM441,097238.07431,466
      Δ49+15.146.32−78135
      Lactose (%)Ctrl334.850.1634.425.05
      CM294.840.1594.385.00
      Δ34−0.0220.0478−0.130.07
      MLY (g/cow per day)Ctrl331,500283.69882,017
      CM291,508305.89322,016
      Δ34+30.646.81−57119
      ECM/DMI (kg/kg)Ctrl441.400.1891.031.72
      CM441.390.2061.041.81
      Δ49+0.0070.0398−0.0800.098
      Milk N/N intake (g/kg)Ctrl4427233.9209347
      CM4427733.3214353
      Δ49+4.59.91−2423
      1 Δ is the change in response of each dependent variable to the substitution of a protein source by CM. A positive value represents an increase in response to CM and a negative value a decrease.
      2 Dry matter intake, MY = milk yield, 4% FCM (

      Gaines, W. L. 1928. The energy basis of measuring milk yield in dairy cows. Illinois Agric. Expt. Sta., Bull. No. 308:401–438. University of Illinois Agricultural Experiment Station, Urbana.

      ), and ECM (
      • Sjaunja L.O.
      • Baevre L.
      • Junkkarinen L.
      • Pedersen J.
      • Setälä J.
      A Nordic proposal for an energy corrected milk (ECM) formula.
      ) in kilograms per cow per day; protein, fat, and lactose = percentages of milk protein, fat, and lactose; MPY, MFY, and MLY = yields of milk protein, fat, and lactose; ECM/DMI = efficiency of milk production; milk N/N intake = apparent N efficiency.
      3 Overall, the database included 88 individual treatment means that generated 49 pairs of Ctrl and CM (for milk lactose % and yield: 62 treatment means and 34 pairs).
      Table 3Lactational responses to substitution of a protein supplement by canola meal (CM) in the full database (All) and in databases separated by the classes of significant experimental conditions or factors
      Response to substitution of a protein supplement by CM was tested within experiment by forcing the regression through the origin (no intercept); regressions were weighted by sample size and controlled for differences in DMI (ΔDMI) and in dietary concentrations of CP (ΔCP) and ether extract (ΔEE). Only factors where the P-value of interaction with CM was ≤0.10 and models where the P-value of CM was ≤0.10 are presented.
      Lactational

      response
      Δ = CM diet mean minus control diet mean; DMI, MY = milk yield, 4% FCM (Gaines, 1928), and ECM (Sjaunja et al., 1990) in kilograms per cow per day; protein, fat, and lactose = percentages of milk protein, fat, and lactose; MPY, MFY, and MLY = yields of milk protein, fat, and lactose in grams per cow per day; ECM/DMI = efficiency of milk production (kg/kg); milk N/N intake = apparent N efficiency (g/kg).
      Database
      Database used: full database (All) or subset database restricted to the experiments for each class of factor. Factors = type of protein: soybean meal (SBM) as a single source or other protein sources (e.g., corn gluten meal, cottonseed meal, and distillers product); type of forages: grass and (or) legume or other forages (e.g., corn or barley silage, alone or mixed with a grass or legume forage); origin of publication: North America or Europe.
      Term
      CM = dietary proportion of CM (100g/kg, DM basis); the coefficient of CM represents the response associated with an increment of 10% CM in the diet, DM basis;
      R
      Δ = CM diet mean minus control diet mean; DMI, MY = milk yield, 4% FCM (Gaines, 1928), and ECM (Sjaunja et al., 1990) in kilograms per cow per day; protein, fat, and lactose = percentages of milk protein, fat, and lactose; MPY, MFY, and MLY = yields of milk protein, fat, and lactose in grams per cow per day; ECM/DMI = efficiency of milk production (kg/kg); milk N/N intake = apparent N efficiency (g/kg).
      P-value

      model
      n
      ΔDMIAll+0.24
      P≤0.001,
      (±0.054) × CM
      0.29<0.00149
      ΔMYAll+0.62
      P≤0.001,
      (±0.084) × CM
      0.53<0.00149
      Δ4% FCMAll+0.50
      P≤0.001,
      (±0.089) × CM
      0.39<0.00149
      SBMNS
      NS = equation is not reported when CM is not significant (P>0.10).
      23
      Other sources+0.59
      P≤0.001,
      (±0.121) × CM
      0.49<0.00126
      ΔECMAll+0.49
      P≤0.001,
      (±0.084) × CM
      0.42<0.00149
      SBMNS23
      Other sources+0.57
      P≤0.001,
      (±0.117) × CM
      0.48<0.00126
      ΔProteinAll+0.020
      P≤0.05 (SE in parenthesis).
      (±0.0081) × CM
      0.110.01849
      SBMNS23
      Other sources+0.047
      P≤0.001,
      (±0.0111) × CM
      0.42<0.00126
      Grass/legume+0.029
      P≤0.01,
      (±0.0090) × CM + 0.13
      P≤0.001,
      (±0.034) × ΔCP
      0.43<0.00130
      Other foragesNS18
      North AmericaNS26
      Europe+0.032
      P≤0.05 (SE in parenthesis).
      (±0.0131) × CM
      0.220.02522
      ΔMPYAll+26.2
      P≤0.001,
      (±3.13) × CM
      0.59<0.00149
      SBM+16.8
      P≤0.001,
      (±3.40) × CM
      0.54<0.00122
      Other sources+32.2
      P≤0.001,
      (±4.33) × CM
      0.69<0.00126
      ΔFatAllNS48
      ΔMFYAll+16.4
      P≤0.01,
      (±4.95) × CM
      0.190.00249
      SBMNS23
      Other sources+24.2
      P≤0.01,
      (±6.53) × CM
      0.350.00126
      ΔLactoseAll−0.022
      P≤0.01,
      (±0.0063) × CM
      0.270.00134
      SBMNS12
      Other sources−0.022
      P≤0.01,
      (±0.0077) × CM
      0.300.01120
      Grass/legumeNS21
      Other forages−0.038
      P≤0.01,
      (±0.0112) × CM
      0.490.00513
      ΔMLYAllNS34
      Grass/legume+25.7
      P≤0.01,
      (±9.35) × CM + 49
      P≤0.05 (SE in parenthesis).
      (±20.0) × ΔDMI
      0.71<0.00121
      Other foragesNS13
      North AmericaNS15
      Europe+44.6
      P≤0.001,
      (±7.45) × CM
      0.67<0.00119
      ΔECM/DMIAll+0.020
      P≤0.001,
      (±0.0048) × CM − 0.05
      P≤0.001,
      (±0.011) × ΔDMI
      0.34<0.00149
      SBMNS23
      Other sources+0.026
      P≤0.01,
      (±0.0071) × CM − 0.06
      P≤0.01,
      (±0.016) × ΔDMI
      0.420.00126
      ΔMilk N/N intakeAll+7.0
      P≤0.001,
      (±1.01) × CM − 11
      P≤0.001,
      (±2.3) × ΔDMI − 13
      P≤0.001,
      (±3.0) × ΔCP
      0.60<0.00149
      SBMNS23
      Other sources+5.8
      P≤0.001,
      (±1.42) × CM − 12
      P≤0.01,
      (±3.1) × ΔDMI + 13
      P≤0.05 (SE in parenthesis).
      (±4.2) × ΔEE
      0.66<0.00126
      1 Response to substitution of a protein supplement by CM was tested within experiment by forcing the regression through the origin (no intercept); regressions were weighted by sample size and controlled for differences in DMI (ΔDMI) and in dietary concentrations of CP (ΔCP) and ether extract (ΔEE). Only factors where the P-value of interaction with CM was ≤0.10 and models where the P-value of CM was ≤0.10 are presented.
      2 Δ = CM diet mean minus control diet mean; DMI, MY = milk yield, 4% FCM (

      Gaines, W. L. 1928. The energy basis of measuring milk yield in dairy cows. Illinois Agric. Expt. Sta., Bull. No. 308:401–438. University of Illinois Agricultural Experiment Station, Urbana.

      ), and ECM (
      • Sjaunja L.O.
      • Baevre L.
      • Junkkarinen L.
      • Pedersen J.
      • Setälä J.
      A Nordic proposal for an energy corrected milk (ECM) formula.
      ) in kilograms per cow per day; protein, fat, and lactose = percentages of milk protein, fat, and lactose; MPY, MFY, and MLY = yields of milk protein, fat, and lactose in grams per cow per day; ECM/DMI = efficiency of milk production (kg/kg); milk N/N intake = apparent N efficiency (g/kg).
      3 Database used: full database (All) or subset database restricted to the experiments for each class of factor. Factors = type of protein: soybean meal (SBM) as a single source or other protein sources (e.g., corn gluten meal, cottonseed meal, and distillers product); type of forages: grass and (or) legume or other forages (e.g., corn or barley silage, alone or mixed with a grass or legume forage); origin of publication: North America or Europe.
      4 CM = dietary proportion of CM (100 g/kg, DM basis); the coefficient of CM represents the response associated with an increment of 10% CM in the diet, DM basis;
      *** P ≤ 0.001,
      ** P ≤ 0.01,
      * P ≤ 0.05 (SE in parenthesis).
      5 NS = equation is not reported when CM is not significant (P > 0.10).
      As mentioned previously, mixed regression model analyses were also performed for comparison purposes (Supplemental Table 1, available online at http://www.journalofdairyscience.org/). Responses were similar in fixed and mixed analyses when the same independent covariates were included in the models. Differences between fixed and mixed models were detected only on milk fat responses because different covariates were retained in the final models. Milk fat concentration did not respond to CM substitution in the fixed model because ΔEE (P = 0.073) was also retained in the final model; on the other hand, milk fat concentration responded negatively to CM substitution in the mixed model, as ΔEE was not significant. Milk fat yield responded positively to CM substitution in the fixed model (Table 3) but not in the mixed model, as ΔDMI (P < 0.001) and ΔCP (P = 0.050) were also retained in the final model. Based on the results observed between both approaches and on the fact that variables were coined as fixed factors in our meta-analysis, only results from the fixed regression models are reported and discussed in the article.
      Overall, the daily intake of CM did not exceed 4.0 kg/cow and the CM proportion in the treated diets averaged 11.7% of total DMI or 1.17 (SD 0.267) when the unit was expressed as 100 g/kg (Table 1). The dietary proportion of CM did not exceed 17.2% in the database. As a prerequisite for inclusion of studies in the database, the proportion of CM in the protein supplement portion of CM needed to be above 0.85 in all diets; in fact, it was 1.00 in all diets except 5.

      Analysis of Lactational Responses

      Hereafter, “response to the substitution of a protein source by CM” will be referred more conveniently as “response to substitution.” As mentioned previously, the coefficients of CM in Table 3 are the responses associated with each increment of 100 g of CM/kg of diet (DM basis) or 10% CM in the diet at equal DMI and equal dietary CP and EE concentrations. In the text below, responses are also discussed for a cow fed the maximum proportion of CM in the database (i.e., 17.2%).

      DMI

      Responses in DMI were mostly positive (35 out of 49). Responses in DMI to substitution increased (P < 0.001) on average by 0.24 kg/cow per day per 10% of CM inclusion in the ration; therefore, DMI would increase by 0.41 kg/cow per day for a cow fed a diet with 17.2% of CM.

      Yields of Milk, 4% FCM, and ECM, and Milk Efficiency

      Most milk yield responses were positive (42 out of 49; Figure 1, panel B). As a result, milk yield responded positively (P < 0.001) to the substitution of a protein supplement by CM (Table 3). In line with the response in milk yield, 4% FCM, ECM, and milk efficiency (ΔECM/DMI) also responded positively (P < 0.001) to substitution (Table 3). For a cow fed CM at 17.2% in the diet, yields of milk, 4% FCM, and ECM would increase by 1.07, 0.85, and 0.84 kg/cow per day, respectively, and milk efficiency would increase by 0.035 kg/kg.

      Milk Composition

      Overall, the response in milk protein percentage was positive (P < 0.02) and would increase by 0.034% for a cow fed CM at 17.2% in the diet (Table 3). Milk fat percentage did not respond to CM substitution but tended to be associated negatively with changes in EE concentrations [P-value in superscript: ΔEE: −0.02NS (±0.0020) × CM − 0.110.07 (±0.058) × ΔEE (R2 = 0.18; n = 48)]. This effect could result from alterations in rumen biohydrogenation of dietary PUFA (
      • Bauman D.E.
      • Griinari J.M.
      Nutritional regulation of milk fat synthesis.
      ). Conversely to the response in milk protein percentage, the response in milk lactose percentage was negative (P = 0.001) and would decrease by 0.038% for a cow fed CM at 17.2% in the diet.

      Yields of Milk Components and Apparent N Efficiency

      Responses in yields of milk protein and fat to substitution were positive (P ≤ 0.002; Table 3). For a cow fed CM at 17.2% in the diet, milk protein and fat yields would increase by 45 and 28 g/cow per day, respectively. The response in milk protein yield is better due to the combined positive effect on milk yield and milk protein percentage. The response in apparent N efficiency (Δmilk N/N intake) to substitution was positive (P < 0.001) mainly due to the positive effect on milk protein secretion (Table 3). The apparent N efficiency would increase by 12 g of milk N/kg of N intake for a cow fed CM at 17.2% in the diet.

      Analysis of Factors

      Type of Protein Source

      Several lactation parameters were affected by the type of protein source that was substituted by CM. The response in milk yield to substitution was not affected by the type of protein source; however, responses in 4% FCM and ECM to substitution were positive only in experiments in which CM replaced protein sources other than SBM (Pinteraction = 0.092 and 0.061 for 4% FCM and ECM, respectively; Table 3). The response in milk protein percentage to substitution was positive only in experiments using protein sources other than SBM (Pinteraction < 0.001; Figure 2); accordingly, the response in milk protein yield to substitution was positive but 2 times higher in experiments using protein sources other than SBM compared with those using SBM (Pinteraction = 0.010; Table 3). In line with responses in 4% FCM and ECM, the response in milk fat yield to substitution was positive only in experiments using protein sources other than SBM compared with those using SBM (Pinteraction = 0.027; Table 3). The response in milk lactose percentage to substitution was negative only in experiments using protein sources other than SBM (Pinteraction = 0.035; Table 3); consequently, the response in milk lactose yield was not affected by the type of protein source.
      Figure thumbnail gr2
      Figure 2Effect of type of protein substituted by canola meal (CM) on responses in milk protein percentage [Δ Protein (%); CM diet minus control diet] (Pinteraction < 0.001): soybean meal [• and solid line; P-value in superscript: −0.010.11 (±0.008) × CM; R2 = 0.11; n = 23] versus other protein sources [○ and dotted line; +0.05≤0.001 (±0.011) × CM; R2 = 0.42; n = 26].
      The response in milk efficiency and in apparent N efficiency to substitution were positive only in experiments using protein sources other than SBM (Pinteraction = 0.10 and 0.040, respectively; Table 3), in accordance with earlier findings on the response in ECM and in milk protein yield to substitution.

      Type of Forage

      The main type of forage in the diet affected the responses in milk protein and lactose percentages to substitution. Milk protein percentage responded positively to substitution only in experiments based on grass and (or) legume forages (Pinteraction = 0.008; Table 3 and Figure 3). The response in milk lactose percentage to substitution was negative only in studies based on corn or barley silage, alone or mixed with a grass or legume forage (Pinteraction = 0.005; Table 3); as a result, the response in milk lactose yield to substitution was positive only in studies based on grass and (or) legume forages (Pinteraction = 0.073; Table 3).
      Figure thumbnail gr3
      Figure 3Effect of type of forages on responses in milk protein percentage [Δ Protein (%); canola meal (CM) diet minus control diet] (Pinteraction = 0.008): grass/legume [• and solid line; P-value in superscript: +0.030.004 (±0.009) × CM + 0.130.001 (±0.034) × ΔCP; R2 = 0.43; n = 30] versus other forages [ and dotted line; −0.01NS (±0.009) × CM + 0.080.005 (±0.026) × ΔCP; R2 = 0.40; n = 18].

      Origin of the Publication

      One can argue that the origin of a publication is not a valid factor to consider in a meta-analysis. However, this factor was taken into consideration in the current meta-analysis because major differences were noted between North American and European studies for the type of protein supplement that was substituted and main forage fed to the cows as outlined below. For example, SBM was used mostly in North America (16 out of 23 studies) whereas other protein sources were used in Europe (14 out of 21 studies; Pearson chi-square, P = 0.016). The experiments based on grass and (or) legume forages were conducted mostly in Europe (21 out of 29 studies) compared with those based on other types of forages (e.g., corn or barley silage), which were conducted entirely in North America (all 15 studies; Pearson chi-square, P < 0.001). Therefore, the origin of study can be seen as a mean to compare diets based on corn or barley silage plus SBM (i.e., typical North American diets in the database) with diets based on grass and (or) legume forages plus corn gluten meal, cottonseed meal, distillers product, or protein mixes (i.e., typical European diets in the database).
      The response in milk protein percentage was positive only in European studies (Pinteraction = 0.004), in accordance with earlier findings on type of protein source and on type of forage because positive responses were reported for grass and (or) legume forages and for supplementation of protein sources other than SBM (i.e., typical European diets in the database). The origin of publication was not associated with the response in milk lactose percentage to substitution; however, the response in milk lactose yield to substitution was positive only in European studies (Pinteraction < 0.001).
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      recently published a meta-analysis comparing CM and SBM in dairy cows fed grass silage-based diets. Seven experiments with grass-based diets supplemented with either SBM or CM were extracted from the whole database. Milk yield and milk protein yield responded positively to substitution [P-value in superscript: +1.160.005 (±0.263) × CM (R2 = 0.76; n = 7) and +35.20.009 (±9.19) × CM (R2 = 0.71; n = 7), respectively], in line with the results of
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      . These responses are greater than those observed in the whole database, which suggest that CM feeding is a good protein supplement for grass-based diets. Indeed, His is the first-limiting AA for milk protein synthesis in cows fed grass silage-based diets supplemented with cereal concentrates (
      • Vanhatalo A.
      • Huhtanen P.
      • Toivonen V.
      • Varvikko T.
      Response of dairy cows fed grass silage diets to abomasal infusions of histidine alone or in combinations with methionine and lysine.
      ) and His concentration is high in CM (
      NRC
      ).

      Estimated MP Supply

      Overall, estimated total MP supplies (g/cow per day) were +2,037 (SD 385.6; range 1,454 to 1,889; n = 44) in control diets and +1,893 (SD 322.3; range 1,395 to 2,778; n = 44) in CM diets (data not shown). Unexpectedly, given the positive responses in milk yield and milk protein yield to CM substitution, changes in supplies of total MP and MP from RUP (g/cow per day;
      NRC
      ) responded negatively to CM substitution [P-value in superscript: −71.7≤0.001 (±12.11) × CM (R2 = 0.42, n = 49) and −86.7≤0.001 (±12.20) × CM (R2 = 0.51, n = 49), respectively]. Changes in estimated supply of MP from bacteria did not respond to CM substitution. When looking at the relationship between changes in estimated total MP supply and responses in milk protein yield (Figure 4), 60% of data points are mismatched (i.e., negative changes in MP supply are associated with positive changes in milk protein yield). Furthermore, this observation holds up even when the database is restricted to the experiments in which the CP concentration of CM is reported (i.e., closed circles in Figure 4), the latter providing a more accurate estimate of diet composition and MP.
      Figure thumbnail gr4
      Figure 4Relationship between responses in milk protein yield [Δ MPY; canola meal (CM) diet minus control diet] and changes in predicted MP supply from the
      NRC
      model (• = CP concentration of CM reported in the publication; = CP content of CM not reported in the publication).

      Discussion

      The current meta-analysis differs from that of
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      on the types of studies used, inclusion criteria and other methodological aspects; as a result, only 5 publications are common to both databases.
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      aimed at avoiding the absence of effect sometimes observed in experiments comparing CM and SBM supplemented at a single rate of inclusion, especially when requirements for production were exceeded in both diets. Therefore, a criteria for inclusion in their database was that CM or SBM was used at different rates in the diet (i.e., 2 or more dietary CP concentrations), but did not require comparison of both CM and SBM in the same experiment. Their meta-analysis was limited to grass silage-based diets because few experiments were found in which CM was included at different rates in corn- or alfalfa silage-based diets. A within-experiment analysis was conducted separately for CM and SBM on the effect of increasing dietary CP content on the performance of dairy cows. The results for CM and SBM were then compared using a t-test. Regardless of the major methodological differences used in both meta-analyses, it is important to note that similar results were obtained whether the meta-analysis was conducted on varying dietary CP concentrations or on isonitrogenous diets (i.e., current meta-analysis).
      Experiments using a Latin square or a crossover design were not a priori excluded from the database because only 5 studies used a randomized block design. The effect of design (Latin square vs. randomized block) was tested for each response and found not significant. Latin square or crossover designs represent a challenge in nutritional meta-analyses because of potential carryover effects from one period to another. Such carryover effects result in smaller variance of responses and they reduce the potential of finding significant differences between estimates. An example of such carryover effects was provided by
      • Cant J.P.
      • Fredeen A.H.
      • MacIntyre T.
      • Gunn J.
      • Crowe N.
      Effect of fish oil and monensin on milk composition in dairy cows.
      who used a 4 × 4 Latin square design with 21-d periods to determine the effect of redfish oil and monensin sodium on milk composition. The authors detected docosahexaenoic acid in the milk of 2 cows treated with monensin sodium at a level similar to that of cows fed redfish oil in the previous period and concluded that carryover effects could not be excluded in their experiment.
      • Cant J.P.
      • Fredeen A.H.
      • MacIntyre T.
      • Gunn J.
      • Crowe N.
      Effect of fish oil and monensin on milk composition in dairy cows.
      postulated that docosahexaenoic acid was stored preferentially in adipose tissues and subsequently released over a time frame of a few weeks. Although carryover effects cannot be completely excluded from protein substitution studies, they are probably of much shorter duration. For example, responses in milk protein concentration to jugular infusion of essential AA (EAA) mixtures reach a plateau quickly (e.g., within 3 d of treatment;
      • Metcalf J.A.
      • Crompton L.A.
      • Wray-Cahen D.
      • Lomax M.A.
      • Sutton J.D.
      • Beever D.E.
      • MacRae J.C.
      • Bequette B.J.
      • Backwell F.R.C.
      • Lobley G.E.
      Responses in milk constituents to intravascular administration of two mixtures of amino acids to dairy cows.
      ) and are restored to initial values soon after treatment cessation (from our own observations).
      The current meta-analysis was conducted on responses rather than on absolute values for 3 reasons. First, as mentioned previously, the inclusion of a quadratic term was not feasible in the within-experiment analysis because only 1 study had more than 2 data points. Second, both methods were tested in parallel and yielded similar within-experiment results given that the regression on responses was forced through the origin. Third, as can be seen in Figure 1A and B, graphs based on responses allow a better visualization of the relationships between 2 variables, especially when all experiments share the same absolute zero value.
      It should be pointed out that coefficient of determination values reported in Table 3 might appear low. However, this was anticipated given the time span covered by the literature search, the different feeding regimens, and the experimental conditions varying among trials. As expected, higher coefficient of determination values were obtained when the response was tested in subset databases of significant factors or with the inclusion of 2 or more covariates.
      Response in DMI to substitution was positive for the range of dietary CM inclusion in the database [i.e., up to 17% (DM basis)]. Experiments with CM concentrations greater than 20% can be found in the literature but they were not eligible in this meta-analysis because either N intake was missing (e.g.,
      • Andersen P.E.
      Double low rapeseed meal in diets to dairy cows.
      ) or CM was not compared with another protein source [e.g.,
      • Rae R.C.
      • Ingalls J.R.
      • McKirdy J.A.
      Response of dairy cows to formaldehyde-treated canola meal during early lactation.
      and
      • Ha J.K.
      • Kennelly J.J.
      Effect of protein on nutrient digestion and milk production by Holstein cows.
      ].
      • Emanuelson M.
      • Ahlin K.-Å.
      • Wiktorsson H.
      Long-term feeding of rapeseed meal and full-fat rapeseed of double low cultivars to dairy cows.
      cited earlier European work in which DMI was decreased with dietary proportions of CM greater than 25%, irrespective of the glucosinolate content in the meal. More recently,
      • Mancuso W.A.
      • Santini F.J.
      • Gagliostro G.A.
      • Rearte D.H.
      Supplementacíon con harina de girasol o de colza “00” a vacas lecheras alimentadas con dietas en base a silaje de maíz.
      fed cows with CM or sunflower meal at 30% of the diet and reported no difference in DMI and milk production. Given that our findings are valid only for the range of CM inclusion in the database, more research is warranted to determine the responses associated with an inclusion rate of CM greater than 17%.
      Even though small differences in DMI might be difficult to highlight in a single study due to lack of statistical power, significant increases were nonetheless reported when cottonseed meal (
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      ) and SBM (
      • Vanhatalo A.
      • Pahkala E.
      • Salo-Väänänen P.
      • Korhonen H.
      • Piironen V.
      • Huhtanen P.
      Rapeseed and soybean as protein supplements of dairy cows fed grass silage based diets.
      ;

      Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

      ) were replaced by CM. The positive effect of CP supplementation on DMI has been attributed to the increased diet digestibility (
      • Oldham J.D.
      Protein-energy interrelationships in dairy cows.
      ). However,
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      reported that CM feeding increased DMI more than SBM in grass silage-based diets but with no difference in OM digestibility between CM and SBM.
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      postulated that DMI might increase in response to a greater energy demand due to the improved lactational performance (“pull effect”) likely induced by a greater supply of AA or a more balanced profile of AA when CM was fed in place of SBM.
      Overall, yields of milk, 4% FCM, and ECM responded positively to substitution of a protein source by CM as most responses were positive. The magnitude of these responses was within ±3 kg/cow per day. Trends or significant increases were reported in the following studies from our database:
      • Laarveld B.
      • Christensen D.A.
      Rapeseed meal in complete feeds for dairy cows.
      ,
      • Papas A.
      • Ingalls J.R.
      • Cansfield P.
      Effects of Tower and 1821 rapeseed meals and Tower gums on milk yield, milk composition and blood parameters of lactating dairy cows.
      ,
      • Dewhurst R.J.
      • Aston K.
      • Fisher W.J.
      • Evans R.T.
      • Dhanoa M.S.
      • McAllan A.B.
      Comparison of energy and protein sources offered at low levels in grass-silage-based diets for dairy cows.
      ,
      • Shingfield K.J.
      • Vanhatalo A.
      • Huhtanen P.
      Comparison of heat-treated rapeseed expeller and solvent-extracted soya-bean meal as protein supplements for dairy cows given grass silage-based diets.
      ,
      • Vanhatalo A.
      • Pahkala E.
      • Salo-Väänänen P.
      • Korhonen H.
      • Piironen V.
      • Huhtanen P.
      Rapeseed and soybean as protein supplements of dairy cows fed grass silage based diets.
      , and

      Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

      .
      Overall, the substitution of a protein source by CM had a positive effect on the responses in milk protein percentage and in milk protein yield. Response in milk protein percentage to substitution was positive despite a potential dilution effect caused by the increase in milk production. Although significant, the regression model accounted poorly (R2 = 0.11) for the variance, suggesting that the response in milk protein percentage to substitution was variable across experiments. Indeed, a better fit to the data was obtained in subset databases separated by the type of protein source or type of forage. For example, the regression model had a better fit to the data when the database was restricted to the experiments in which protein sources other than SBM were substituted by CM (R2 = 0.42). Trends or significant increases in milk protein percentage were reported when CM replaced cottonseed meal (
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      ), dried distillers grains with solubles (
      • Christen K.A.
      • Schingoethe D.J.
      • Kalscheur K.F.
      • Hippen A.R.
      • Karges K.K.
      • Gibson M.L.
      Response of lactating dairy cows to high protein distillers grains or 3 other protein supplements.
      ), untreated barley distillers solubles (
      • Huhtanen P.
      • Khalili H.
      • Näsi M.
      A comparison of untreated and formaldehyde-treated barley distiller's solubles and rapeseed meal as protein supplements in dairy cows given grass silage ad libitum.
      ), or a mixture of SBM, cottonseed meal, and coconut cake (
      • Emanuelson M.
      • Ahlin K.-Å.
      • Wiktorsson H.
      Long-term feeding of rapeseed meal and full-fat rapeseed of double low cultivars to dairy cows.
      ).
      Unlike the response in milk protein percentage, a much better fit to the data was obtained for the response in milk protein yield to substitution across the whole data set (R2 = 0.59). As expected, the response was positive and significant in many publications:
      • Huhtanen P.
      • Khalili H.
      • Näsi M.
      A comparison of untreated and formaldehyde-treated barley distiller's solubles and rapeseed meal as protein supplements in dairy cows given grass silage ad libitum.
      ,
      • Emanuelson M.
      • Ahlin K.-Å.
      • Wiktorsson H.
      Long-term feeding of rapeseed meal and full-fat rapeseed of double low cultivars to dairy cows.
      ,
      • Shingfield K.J.
      • Vanhatalo A.
      • Huhtanen P.
      Comparison of heat-treated rapeseed expeller and solvent-extracted soya-bean meal as protein supplements for dairy cows given grass silage-based diets.
      ,
      • Vanhatalo A.
      • Pahkala E.
      • Salo-Väänänen P.
      • Korhonen H.
      • Piironen V.
      • Huhtanen P.
      Rapeseed and soybean as protein supplements of dairy cows fed grass silage based diets.
      ,

      Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

      ,
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      , and
      • Oba M.
      • Penner G.B.
      • Whyte T.D.
      • Wierenga K.
      Effects of feeding triticale dried distillers grains plus solubles as a nitrogen source on productivity of lactating dairy cows.
      .
      The positive responses in milk production and in milk protein percentage and yield suggest a more efficient N metabolism, which might be related to a greater supply of EAA or an improvement in the duodenal profile of supplied EAA (e.g., His, Lys, and Met). In a comprehensive meta-analysis of experiments studying substitution of various protein sources (
      • Ipharraguerre I.R.
      • Clark J.H.
      Impacts of the source and amount of crude protein on the intestinal supply of nitrogen fractions and performance of dairy cows.
      ), milk protein yield was depressed by the replacement of solvent-extracted SBM by corn byproducts. The authors suggested that the effect could be attributed to the poor quality of the supplied MP due to the low Lys content of corn byproducts. Supporting the EAA profile hypothesis, the response in milk protein percentage was positive in the current meta-analysis only when CM substituted protein sources other than SBM (i.e., corn byproducts that would offer a profile of EAA not as good as SBM or CM).
      Improving the quantity of potentially limiting AA or the profile of EAA in duodenal flow can be achieved by stimulating rumen microbial synthesis or by feeding protein sources that would deliver a digestible RUP fraction with an EAA profile closer to that of milk.
      • Clark J.H.
      Lactational responses to postruminal administration of proteins and amino acids.
      and
      • Choung J.-J.
      • Chamberlain D.G.
      The effects of abomasal infusions of casein or soya-bean-protein isolate on the milk production of dairy cows in mid-lactation.
      demonstrated that the duodenal EAA profile should be similar to that of milk protein to increase milk protein percentage. Microbial protein has a favorable AA profile that closely matches AA requirements for milk and milk protein synthesis (
      NRC
      ). Canola meal can stimulate rumen microbial synthesis because it provides an adequate supply of RDP, which is important because N forms other than ammonia are needed for maximum microbial growth (
      • Griswold K.E.
      • Hoover W.H.
      • Miller T.K.
      • Thayne W.V.
      Effect of form of nitrogen on growth of ruminal microbes in continuous culture.
      ). For instance, preformed free AA and peptides can be readily incorporated into the microbial mass (
      • Bach A.
      • Calsamiglia S.
      • Stern M.D.
      Nitrogen metabolism in the rumen.
      ). The EAA profile of the digestible RUP fraction is also important as it is complementary to the supply of microbial protein. Even though canola meal supplies less duodenal RUP compared with most protein sources, it contains relatively high concentrations of Met, Lys, and His (often quoted as being limiting EAA in dairy rations) compared with other common protein sources (e.g., corn gluten meal or distillers products;
      NRC
      ).
      • Chibisa G.E.
      • Christensen D.A.
      • Mutsvangwa T.
      Effects of replacing canola meal as the major protein source with wheat dried distillers grains with solubles on ruminal function, microbial protein synthesis, omasal flow, and milk production in cows.
      found that omasal flow of bacterial NAN, as a percentage of total NAN, increased when CM replaced increasing amounts of wheat dried distillers grains with solubles, which contains much less RDP.
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      compared CM with cottonseed meal (i.e., a poor RDP source) and reported higher MUN and urinary urea N excretion with CM compared with cottonseed meal, in accordance with the high RDP content of the CM diet estimated using the
      NRC
      model or measured in vivo. In contrast with this apparent loss in N efficiency, milk protein percentage and yield were greater with CM compared with cottonseed meal, thus reflecting an overall improved N metabolism. In a companion paper,
      • Brito A.F.
      • Broderick G.A.
      • Reynal S.M.
      Effects of different protein supplements on omasal nutrient flow and microbial protein synthesis in lactating dairy cows.
      reported that omasal NAN flows (% of N intake) were similar for CM and cottonseed meal but that CM supplied less RUP (% of DMI) and more microbial NAN (% of total NAN) compared with cottonseed meal. These results suggest that CM improved milk protein percentage and yield due to better quality of EAA in omasal NAN flow and that the AA profile of RUP from cottonseed meal might not be complementary to that of microbial protein. Indeed, CM has a greater Lys content compared with cottonseed meal (
      NRC
      ) and, although omasal Lys flows were similar for CM and cottonseed meal (
      • Brito A.F.
      • Broderick G.A.
      • Reynal S.M.
      Effects of different protein supplements on omasal nutrient flow and microbial protein synthesis in lactating dairy cows.
      ), metabolizable Lys could differ between diets depending on the intestinal digestibility coefficient applied for Lys.
      Several responses to substitution were clearly better when CM replaced protein sources other than SBM (as opposed to CM vs. SBM): greater response in 4% FCM and ECM, greater response in milk protein percentage and yield, as well as greater milk and apparent N efficiencies. The response in milk protein yield was 2 times greater when CM replaced protein sources other than SBM: this suggests that the EAA profile of SBM might be more balanced or complementary to that of microbial protein compared with the EAA profile of other protein sources (e.g., distillers grain, corn gluten meal, and cottonseed meal). Indeed, it must be acknowledged that substitution of CM by other protein sources may have altered not just the supply of AA but also the provision of carbohydrates, FA, minerals, and vitamins, all of them potentially contributing to the effects observed. However, the scarcity of details about diet composition precluded a thorough examination of the contribution of these factors on the responses to CM substitution.
      In a further attempt to determine why CM had an overall positive effect on milk production, estimated MP supply was calculated with the
      NRC
      model, as diets might be isonitrogenous but have different MP supply. A good relationship was expected between changes in MP supply and changes in milk protein secretion, given that MP supply is a key factor that determines milk protein secretion in lactating dairy cows (
      NRC
      ;
      • Broderick G.A.
      • Huhtanen P.
      • Ahvenjärvi S.
      • Reynal S.M.
      • Shingfield K.J.
      Quantifying ruminal nitrogen metabolism using the omasal sampling technique in cattle—A meta-analysis.
      ). However, MP supply appeared underestimated in CM diets, as most positive changes in milk protein yield were coupled with negative changes in estimated MP supply.
      • Huhtanen P.
      • Hetta M.
      • Swensson C.
      Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
      also had concerns about accuracy of most feed evaluation systems in estimating MP supply and reported that these systems usually underestimate MP concentration of CM relative to SBM. A low estimation of MP supply by CM can result from low estimations of either microbial synthesis or RUP, and also from low intestinal digestibility of RUP (% of RUP).
      • Brito A.F.
      • Broderick G.A.
      • Reynal S.M.
      Effects of different protein supplements on omasal nutrient flow and microbial protein synthesis in lactating dairy cows.
      reported a lower estimation of microbial protein synthesis by 14% using
      NRC
      model compared with in vivo estimations using 15N in omasal samples. The authors speculated that the discrepancy might be due to low purine recovery, thus biasing estimations of total purines, which are the primary microbial marker used for the microbial protein measurements in many studies included in the
      NRC
      database. Furthermore,
      • Brito A.F.
      • Broderick G.A.
      • Reynal S.M.
      Effects of different protein supplements on omasal nutrient flow and microbial protein synthesis in lactating dairy cows.
      pointed out that variations in prediction of microbial synthesis might occur because many factors affecting microbial synthesis were not accounted for in the
      NRC
      model (e.g., form of N available in the rumen, source and amount of carbohydrate, amount and source of fat, ruminal pH, and dilution rate). This would negatively affect estimated MP supply of high-RDP protein sources such as CM, especially when compared with less-degradable protein sources (e.g., corn gluten meal, cottonseed meal, and distillers products). The intestinal digestibility of RUP of the expeller CM used by
      NRC
      is relatively low (75% of RUP) compared with that of other protein sources (e.g., 93% for mechanically or solvent-extracted SBM and 92% for cottonseed meal). For example, in
      • Brito A.F.
      • Broderick G.A.
      Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
      , changes in estimated MP supply were approximately −140 and −270 g/cow per day (calculated values not reported in the article) when CM replaced solvent SBM and cottonseed meal, respectively. These negative MP supplies were associated with positive responses in milk protein yield. Increasing the intestinal digestibility of RUP of CM from 75 to 85% of RUP would add approximately 60 g to the supply of MP from CM. The intestinal digestibility of RUP of CM should be investigated further, given that it could account for a nonnegligible part of the discrepancy observed between changes in estimated MP supply and responses in milk protein yield.

      Conclusions

      The present meta-analysis was based on experiments using isonitrogenous diets (±1% dietary CP between the control and the treatment diets) in which CM substituted a protein source in dairy rations. Dry matter intake responded positively to substitution; therefore, all regressions were controlled for differences in DMI and differences in dietary CP and EE concentrations between the control and the CM diet. Yields of milk, 4% FCM, and ECM, as well as yields of milk protein and fat responded positively to substitution. Furthermore, milk protein percentage also responded positively to substitution. Among the experimental conditions or factors that significantly influenced the responses to substitution, it is noteworthy to mention the type of forage (e.g., grass or legume forages vs. corn or barley silage) and the type of protein substituted (SBM or other protein sources). Milk protein percentage responded positively when cows were fed grass or legume forages but the response was not significant when cows were fed other forages. On the other hand, responses in yields of 4% FCM and ECM, milk protein percentage and yield, milk efficiency, and apparent N efficiency were all positive or greater when CM replaced a protein source based on corn gluten meal, cottonseed meal, distillers product, or a protein mix that included or not some SBM. Milk protein yield and apparent N efficiency responded positively to substitution, suggesting an increment in duodenal MP supplied or a better profile of duodenal EAA. However, most positive responses in milk protein yield were coupled with negative changes in estimated MP supply, reflecting a systematic underestimation of MP supply associated with CM inclusion in dairy rations in the current database using the
      NRC
      model.

      Acknowledgments

      Sincere thanks to the Canola Council of Canada (Winnipeg, MB, Canada) and Agriculture and Agri-Food Canada (Sherbrooke, QC, Canada) for funding this research.

      Appendix A. Supplementary data

      References

        • Abbadi A.
        • Leckband G.
        Rapeseed breeding for oil content, quality, and sustainability.
        Eur. J. Lipid Sci. Technol. 2011; 113: 1198-1206
        • Andersen P.E.
        Double low rapeseed meal in diets to dairy cows.
        in: Sørensen H. Advances in the Production and Utilization of Cruciferous Crops. Martinus Nijhoff Publishers, The Hague, the Netherlands1985: 218-221
        • Bach A.
        • Calsamiglia S.
        • Stern M.D.
        Nitrogen metabolism in the rumen.
        J. Dairy Sci. 2005; 88: E9-E21
        • Bauman D.E.
        • Griinari J.M.
        Nutritional regulation of milk fat synthesis.
        Annu. Rev. Nutr. 2003; 23: 203-227
        • Bell J.M.
        Nutrients and toxicants in rapeseed meal: A review.
        J. Anim. Sci. 1984; 58: 996-1010
        • Brito A.F.
        • Broderick G.A.
        Effects of different protein supplements on milk production and nutrient utilization in lactating dairy cows.
        J. Dairy Sci. 2007; 90: 1816-1827
        • Brito A.F.
        • Broderick G.A.
        • Reynal S.M.
        Effects of different protein supplements on omasal nutrient flow and microbial protein synthesis in lactating dairy cows.
        J. Dairy Sci. 2007; 90: 1828-1841
        • Broderick G.A.
        • Huhtanen P.
        • Ahvenjärvi S.
        • Reynal S.M.
        • Shingfield K.J.
        Quantifying ruminal nitrogen metabolism using the omasal sampling technique in cattle—A meta-analysis.
        J. Dairy Sci. 2010; 93: 3216-3230
        • Cant J.P.
        • Fredeen A.H.
        • MacIntyre T.
        • Gunn J.
        • Crowe N.
        Effect of fish oil and monensin on milk composition in dairy cows.
        Can. J. Anim. Sci. 1997; 77: 125-131
        • Chibisa G.E.
        • Christensen D.A.
        • Mutsvangwa T.
        Effects of replacing canola meal as the major protein source with wheat dried distillers grains with solubles on ruminal function, microbial protein synthesis, omasal flow, and milk production in cows.
        J. Dairy Sci. 2012; 95: 824-841
        • Choung J.-J.
        • Chamberlain D.G.
        The effects of abomasal infusions of casein or soya-bean-protein isolate on the milk production of dairy cows in mid-lactation.
        Br. J. Nutr. 1993; 69: 103-115
        • Christen K.A.
        • Schingoethe D.J.
        • Kalscheur K.F.
        • Hippen A.R.
        • Karges K.K.
        • Gibson M.L.
        Response of lactating dairy cows to high protein distillers grains or 3 other protein supplements.
        J. Dairy Sci. 2010; 93: 2095-2104
        • Christensen D.A.
        • McKinnon P.J.
        Canola meal for beef and dairy cattle.
        in: Clandinin D.R. Canola Meal for Livestock and Poultry. Canola Council of Canada, Winnipeg, MB, Canada1989: 19-22
        • Clark J.H.
        Lactational responses to postruminal administration of proteins and amino acids.
        J. Dairy Sci. 1975; 58: 1178-1197
        • Cook R.D.
        Influential observations in linear regression.
        J. Am. Stat. Assoc. 1979; 74: 169-174
        • DePeters E.J.
        • Bath D.L.
        Canola meal versus cottonseed meal as a protein supplement in dairy rations.
        J. Dairy Sci. 1986; 69: 148-154
        • Dewhurst R.J.
        • Aston K.
        • Fisher W.J.
        • Evans R.T.
        • Dhanoa M.S.
        • McAllan A.B.
        Comparison of energy and protein sources offered at low levels in grass-silage-based diets for dairy cows.
        Anim. Sci. 1999; 68: 789-799
      1. Emanuelson, M. 1989. Rapeseed Products of Double Low Cultivars to Dairy Cows: Effects of Long-Term Feeding and Studies on Rumen Metabolism. Volume 189 of Rapport. Swedish University of Agricultural Sciences, Department Animal Nutrition Management, Uppsala, Sweden.

        • Emanuelson M.
        • Ahlin K.-Å.
        • Wiktorsson H.
        Long-term feeding of rapeseed meal and full-fat rapeseed of double low cultivars to dairy cows.
        Livest. Prod. Sci. 1993; 33: 199-214
        • Fenwick G.R.
        The assessment of a new protein source—Rapeseed.
        Proc. Nutr. Soc. 1982; 41: 277-288
        • Fisher L.J.
        • Walsh D.S.
        Substitution of rapeseed meal for soybean meal as a source of protein for lactating cows.
        Can. J. Anim. Sci. 1976; 56: 233-242
      2. Gaines, W. L. 1928. The energy basis of measuring milk yield in dairy cows. Illinois Agric. Expt. Sta., Bull. No. 308:401–438. University of Illinois Agricultural Experiment Station, Urbana.

        • Glasser F.
        • Ferlay A.
        • Chilliard Y.
        Oilseed lipid supplements and fatty acid composition of cow milk: A meta-analysis.
        J. Dairy Sci. 2008; 91: 4687-4703
        • Grenet N.
        • Journet M.
        Rapeseed oil meal in animal feeding. III. Influence of the processing method and the proportion of oil meal in concentrate fed to dairy cows on feed intake, milk production and composition.
        Ann. Zootech. 1971; 20: 437-449
        • Griswold K.E.
        • Hoover W.H.
        • Miller T.K.
        • Thayne W.V.
        Effect of form of nitrogen on growth of ruminal microbes in continuous culture.
        J. Anim. Sci. 1996; 74: 483-491
        • Ha J.K.
        • Kennelly J.J.
        Effect of protein on nutrient digestion and milk production by Holstein cows.
        J. Dairy Sci. 1984; 67: 2302-2307
        • Harker K.N.
        • O’Donovan J.T.
        • Turkington T.K.
        • Blackshaw R.E.
        • Lupwayi N.Z.
        • Smith E.G.
        • Klein-Gebbinck H.
        • Dosdall L.M.
        • Hall L.M.
        • Willenborg C.J.
        • Kutcher H.R.
        • Malhi S.S.
        • Vera C.L.
        • Gan Y.
        • Lafond G.P.
        • May W.E.
        • Grant C.A.
        • McLaren D.L.
        High-yield no-till canola production on the Canadian prairies.
        Can. J. Plant Sci. 2012; 92: 221-233
        • Hickling D.
        Maximized utilization of canola co-products in livestock industry.
        Proc. 29th Western Nutrition Conference, Edmonton, AB, Canada. University of Alberta, Edmonton, AB, Canada2008 (Pages 3–14)
        • Hill R.
        Rapeseed meal in the diets of ruminants.
        Nutr. Abstr. Rev. 1991; 61 (Series B): 139-155
        • Huhtanen P.
        Supply of nutrients and productive responses in dairy cows given diets based on restrictively fermented silage.
        Agric. Food Sci. 1998; 7: 219-250
        • Huhtanen P.
        • Hetta M.
        • Swensson C.
        Evaluation of canola meal as a protein supplement for dairy cows: A review and a meta-analysis.
        Can. J. Anim. Sci. 2011; 91: 529-543
        • Huhtanen P.
        • Jaakkola S.
        • Saarisalo E.
        The effects of concentrate energy source on the milk production of dairy cows given a grass silage-based diet.
        Anim. Sci. 1995; 60: 31-40
        • Huhtanen P.
        • Khalili H.
        • Näsi M.
        A comparison of untreated and formaldehyde-treated barley distiller's solubles and rapeseed meal as protein supplements in dairy cows given grass silage ad libitum.
        J. Agric. Sci. Finland. 1991; 63: 455-463
        • Ingalls J.R.
        • Sharma H.R.
        Feeding of Bronowski, Span and commercial rapeseed meals with or without addition of molasses or flavor in the rations of lactating cows.
        Can. J. Anim. Sci. 1975; 55: 721-729
        • Ipharraguerre I.R.
        • Clark J.H.
        Impacts of the source and amount of crude protein on the intestinal supply of nitrogen fractions and performance of dairy cows.
        J. Dairy Sci. 2005; 88: E22-E37
        • Khalili H.
        • Sairanen A.
        • Hissa K.
        • Huhtanen P.
        Effects of type and treatment of grain and protein source on dairy cow performance.
        Anim. Sci. 2001; 72: 573-584
        • Khorasani G.R.
        • Robinson P.H.
        • Kennelly J.J.
        Evaluation of solvent and expeller linseed meals as protein sources for dairy cattle.
        Can. J. Anim. Sci. 1994; 74: 479-485
        • Kleinbaum D.G.
        • Kupper L.L.
        • Muller K.E.
        • Nizam A.
        Confounding and interaction in regression.
        in: Kleinbaum G. Kupper L.L. Nizam A. Muller K.E. Applied Regression Analysis and Other Multivariable Methods. Duxbury Press, Pacific Grove, CA1998: 186-211
        • Laarveld B.
        • Christensen D.A.
        Rapeseed meal in complete feeds for dairy cows.
        J. Dairy Sci. 1976; 59: 1929-1935
      3. Littell, R. C., G. A. Milliken, W. W. Stroup, and R. D. Wolfinger. 1996. SAS® system for mixed models. SAS Institute, Cary, NC.

        • Lusas E.E.
        • Riaz M.N.
        Fats and oils in feedstuffs and pet foods.
        in: Shahidi F. Bailey’s Industrial Oil & Fat Product. 6th. John Wiley & Sons, New York, NY2005: 317-390
        • Macleod G.K.
        Canola meal as a protein supplement in corn based dairy rations.
        Research on Canola Seed, Oil, and Meal. Ninth Project Report. Canola Council of Canada, Winnipeg, MB, Canada1991 (Pages 255–264)
        • Mancuso W.A.
        • Santini F.J.
        • Gagliostro G.A.
        • Rearte D.H.
        Supplementacíon con harina de girasol o de colza “00” a vacas lecheras alimentadas con dietas en base a silaje de maíz.
        Rev. Arg. Prod. Anim. 2002; 22: 169-184
        • McClean C.
        • Laarveld B.
        Effect of somatotropin and protein supplement on thyroid function of dairy cattle.
        Can. J. Anim. Sci. 1991; 71: 1053-1061
        • Metcalf J.A.
        • Crompton L.A.
        • Wray-Cahen D.
        • Lomax M.A.
        • Sutton J.D.
        • Beever D.E.
        • MacRae J.C.
        • Bequette B.J.
        • Backwell F.R.C.
        • Lobley G.E.
        Responses in milk constituents to intravascular administration of two mixtures of amino acids to dairy cows.
        J. Dairy Sci. 1996; 79: 1425-1429
      4. Moss, A. 2002. The effects of long-term feeding of extracted rapeseed meal and whole rapeseed on the physical and financial performance, health and welfare of high yielding dairy cows. HGCA Project Report OS59. Home Grown Cereals Authority, London, UK.

        • Mulrooney C.N.
        • Schingoethe D.J.
        • Kalscheur K.F.
        • Hippen A.R.
        Canola meal replacing distillers grains with solubles for lactating dairy cows.
        J. Dairy Sci. 2009; 92: 5669-5676
        • Munger A.
        Use of rapeseed oilmeal concentrates in dairy cow feeding.
        Revue Suisse d’Agriculture. 1996; 28: 133-136
        • NRC
        Nutrient Requirements of Dairy Cattle. 7th rev. National Academy of Sciences, Washington, DC2001
        • Newkirk R.
        Canola Meal: Feed Industry Guide. 4th. Canadian International Grains Institute, Winnipeg, MB, Canada2009
        • Oba M.
        • Penner G.B.
        • Whyte T.D.
        • Wierenga K.
        Effects of feeding triticale dried distillers grains plus solubles as a nitrogen source on productivity of lactating dairy cows.
        J. Dairy Sci. 2010; 93: 2044-2052
        • Oldham J.D.
        Protein-energy interrelationships in dairy cows.
        J. Dairy Sci. 1984; 67: 1090-1114
        • Papas A.
        • Ingalls J.R.
        • Cansfield P.
        Effects of Tower and 1821 rapeseed meals and Tower gums on milk yield, milk composition and blood parameters of lactating dairy cows.
        Can. J. Anim. Sci. 1978; 58: 671-679
        • Piepenbrink M.S.
        • Schingoethe D.J.
        • Brouk M.J.
        • Stegeman G.A.
        Systems to evaluate the protein quality of diets fed to lactating cows.
        J. Dairy Sci. 1998; 81: 1046-1061
        • Rae R.C.
        • Ingalls J.R.
        • McKirdy J.A.
        Response of dairy cows to formaldehyde-treated canola meal during early lactation.
        Can. J. Anim. Sci. 1983; 63: 905-915
      5. Rinne, M., K. Kuoppala, S. Ahvenjärvi, and A. Vanhatalo. 2006. Rypsi soijaa parempi lypsylehmien valkuaistäydennys myös apilapitoista säilörehua syötettäessä. In Maataloustieteen päivät 2006. Suomen Maataloustieteellisen Tiedote 21:1–7. A. Hopponen, ed. Accessed Oct. 5, 2011. http://www.smts.fi/esit06/1002.pdf

        • Robinson P.H.
        • Kennelly J.J.
        Influence of intake of rumen undegradable protein on milk production of late lactation Holstein cows.
        J. Dairy Sci. 1988; 71: 2135-2142
        • Sauvant D.
        • Schmidely P.
        • Daudin J.J.
        • St-Pierre N.R.
        Meta-analyses of experimental data: Application in animal nutrition.
        Animal. 2008; 2: 1203-1214
        • Sharma H.R.
        • Ingalls J.R.
        • McKirdy J.A.
        Effects of feeding a high level of Tower rapeseed meal in dairy rations on feed intake and milk production.
        Can. J. Anim. Sci. 1977; 57: 653-662
        • Shingfield K.J.
        • Jaakkola S.
        • Huhtanen P.
        Effects of level of nitrogen fertilizer application and various nitrogenous supplements on milk production and nitrogen utilization of dairy cows given grass silage-based diets.
        Anim. Sci. 2001; 73: 541-554
        • Shingfield K.J.
        • Vanhatalo A.
        • Huhtanen P.
        Comparison of heat-treated rapeseed expeller and solvent-extracted soya-bean meal as protein supplements for dairy cows given grass silage-based diets.
        Anim. Sci. 2003; 77: 305-317
        • Sjaunja L.O.
        • Baevre L.
        • Junkkarinen L.
        • Pedersen J.
        • Setälä J.
        A Nordic proposal for an energy corrected milk (ECM) formula.
        in: Gaillon P. Chabert Y. Proceedings of the 27th Session of the International Committee of Recording and Productivity of Milk Animal. EEAP Publication No. 50. Pudoc, Wageningen, the Netherlands1990: 156-157
        • St-Pierre N.R.
        Invited review: Integrating quantitative findings from multiple studies using mixed model methodology.
        J. Dairy Sci. 2001; 84: 741-755
        • Stefansson B.R.
        • Kondra Z.P.
        Tower summer rape.
        Can. J. Plant Sci. 1975; 55: 343-344
        • Tan S.H.
        • Mailer R.J.
        • Blanchard C.L.
        • Agboola S.O.
        Canola proteins for human consumption: Extraction, profile, and functional properties.
        J. Food Sci. 2011; 76: R16-R28
        • Thomke S.
        Review of rapeseed meal in animal nutrition: Ruminant animals.
        J. Am. Oil Chem. Soc. 1981; 58: 805-810
        • Tuori M.
        Rapeseed meal as a supplementary protein for dairy cows on grass silage-based diet, with the emphasis on the Nordic AAT-PBV feed protein evaluation system.
        Agric. Sci. Finland. 1992; 1: 375-439
        • Vanhatalo A.
        • Huhtanen P.
        • Toivonen V.
        • Varvikko T.
        Response of dairy cows fed grass silage diets to abomasal infusions of histidine alone or in combinations with methionine and lysine.
        J. Dairy Sci. 1999; 82: 2674-2685
        • Vanhatalo A.
        • Pahkala E.
        • Salo-Väänänen P.
        • Korhonen H.
        • Piironen V.
        • Huhtanen P.
        Rapeseed and soybean as protein supplements of dairy cows fed grass silage based diets.
        in: Sørensen H. Sørensen J.C. Sørensen S. Bellostas Muguerza N. Bjergegaard C. Proc. 11th Intl. Rapeseed Congress: Enhanced Value of Cruciferous Oilseed Crops by Optimal Production and Use of the High Quality Seed Components. The Royal Veterinary and Agricultural University, Copenhagen, Denmark2003: 1238-1240