Advertisement
Research| Volume 105, ISSUE 3, P2558-2571, March 2022

Download started.

Ok

Suitability of milk lactate dehydrogenase and serum albumin for pathogen-specific mastitis detection in automatic milking systems

  • M. Khatun
    Correspondence
    Corresponding author
    Affiliations
    School of Life and Environmental Sciences and Sydney Institute of Agriculture, The University of Sydney, Camden 2570, New South Wales, Australia

    Bangladesh Agricultural University, Mymensingh, Bangladesh, 2202

    Veterinary Physiology, University of Bern, Bremgartenstrasse 109a, 3012 Bern, Switzerland
    Search for articles by this author
  • P.C. Thomson
    Affiliations
    School of Life and Environmental Sciences and Sydney Institute of Agriculture, The University of Sydney, Camden 2570, New South Wales, Australia
    Search for articles by this author
  • S.C. García
    Affiliations
    School of Life and Environmental Sciences and Sydney Institute of Agriculture, The University of Sydney, Camden 2570, New South Wales, Australia
    Search for articles by this author
  • R.M. Bruckmaier
    Affiliations
    Veterinary Physiology, University of Bern, Bremgartenstrasse 109a, 3012 Bern, Switzerland
    Search for articles by this author
Open AccessPublished:January 05, 2022DOI:https://doi.org/10.3168/jds.2021-20475

      ABSTRACT

      In response to intramammary infection (IMI), blood-derived leukocytes are transferred into milk, which can be measured as an increase of somatic cell count (SCC). Additionally, pathogen-dependent IgG increases in milk following infection. The IgG transfer into milk is associated with the opening of the blood-milk barrier, which is much more pronounced during gram-negative than gram-positive IMI. Thus, milk IgG concentration may help to predict the pathogen type causing IMI. Likewise, lactate dehydrogenase (LDH) and serum albumin (SA) cross the blood-milk barrier with IgG if its integrity is reduced. Because exact IgG analysis is complicated and difficult to automate, LDH activity and SA concentration aid as markers to predict the IgG transfer into milk in automatic milking systems (AMS). This study was conducted to test the hypothesis that LDH and SA in milk correlate with the IgG transfer, and in combination with SCC these factors allow the differentiation between gram-positive and gram-negative IMI or even more precisely the infection-causing pathogen. Further, the expression of these parameters in foremilk before (BME) and after (AME) milk ejection was tested. In the AMS, quarter milk samples (n = 686) from 48 Holstein-Friesian cows were collected manually BME and AME, followed by an aseptic sample for bacteriological culture. Mixed models were used to (1) predict the concentration of IgG transmitted from blood into milk based on LDH and SA; (2) use principal component analysis to evaluate joint patterns of SCC (cells/mL), IgG (mg/mL), LDH (U/L), and SA (mg/mL) and use the principal component scores to compare gram-positive, gram-negative, and control IMI types and BME versus AME samples; and (3) predict gram-positive and gram-negative IMI by inclusion of combined SCC-LDH and SCC-SA as predictors in the model. Overall, the SA and LDH had similar ability to predict IgG transmission from blood into milk. Comparing the areas under the curve (AUC) of the receiver operator characteristic curves, the SCC-LDH versus SCC-SA had lower gram-positive (AUC = 0.984 vs. 0.986) but similar gram-negative (AUC = 0.995 vs. 0.998) IMI prediction ability. The SCC, IgG, LDH, and SA were greater in gram-negative than in gram-positive IMI (BME and AME) in early lactation. All measured factors had higher values in milk samples taken BME than AME. In conclusion, LDH and SA could be used as replacement markers to indicate the presence of IgG transfer from blood into milk; in combination with SCC, both SA and LDH are suitable for differentiating IMI type, and BME is better for mastitis detection in AMS.

      Key words

      INTRODUCTION

      Mastitis in dairy cows is an inflammation of the mammary gland that affects animal welfare and has a huge economic impact on the dairy industry worldwide (
      • Halasa T.
      • Huijps K.
      • Østerås O.
      • Hogeveen H.
      Economic effects of bovine mastitis and mastitis management: A review.
      ;
      • Huijps K.
      • Lam T.J.
      • Hogeveen H.
      Costs of mastitis: facts and perception.
      ). Early detection of mastitis including the pathogen causing the disease is the basis for a rapid and efficient treatment decision for the respective mastitis type to minimize the use of antibiotics, to increase animal health and welfare, and to improve the economic return (
      • Milner P.
      • Page K.L.
      • Hillerton J.E.
      The effects of early antibiotic treatment following diagnosis of mastitis detected by a change in the electrical conductivity of milk.
      ;
      • Lehmann M.
      • Wall S.K.
      • Wellnitz O.
      • Bruckmaier R.M.
      Changes in milk L-lactate, lactate dehydrogenase, serum albumin, and IgG during milk ejection and their association with somatic cell count.
      ).
      In the immunopathogenesis of IMI during mastitis, the causal pathogens (mostly gram-positive or gram-negative bacteria, or both) induce a response of the innate immune system, but depending on the pathogen type, also of the acquired or specific immune system in the host (
      • Pyörälä S.
      Indicators of inflammation in the diagnosis of mastitis.
      ;
      • Schmitz S.
      • Pfaffl M.W.
      • Meyer H.H.D.
      • Bruckmaier R.M.
      Short-term changes of mRNA expression of various inflammatory factors and milk proteins in mammary tissue during LPS-induced mastitis.
      ;
      • Wellnitz O.
      • Arnold E.T.
      • Bruckmaier R.M.
      Lipopolysaccharide and lipoteichoic acid induce different immune responses in the bovine mammary gland.
      ). The innate immune response is mainly based on the cellular components consisting of circulatory leukocytes, mainly PMN, which migrate from blood into the mammary tissue and milk (
      • Wellnitz O.
      • Bruckmaier R.M.
      The innate immune response of the bovine mammary gland to bacterial infection.
      ). This paracellular invasion of PMN is driven by chemoattractants mainly produced by macrophages in the milk as well as mammary epithelial cells that are in contact with bacteria and their pathogen-associated molecular patterns (
      • Burton J.L.
      • Erskine R.J.
      Immunity and mastitis. Some new ideas for an old disease.
      ;
      • Wagner S.A.
      • Jones D.E.
      • Apley M.D.
      Effect of endotoxic mastitis on epithelial cell numbers in the milk of dairy cows.
      ;
      • Bruckmaier R.M.
      • Wellnitz O.
      Triennial Lactation Symposium/BOLFA: Pathogen-specific immune response and changes in the blood–milk barrier of the bovine mammary gland.
      ). To support the innate immune response against the invading mastitis pathogens, there is an elicitation of acquired immune response and transfer of IgG from blood into milk (
      • Burton J.L.
      • Erskine R.J.
      Immunity and mastitis. Some new ideas for an old disease.
      ). It has been shown that the transferred IgG reflects the whole spectrum of antibodies that are present in blood (i.e., not only those directed against the invaded pathogen;
      • Burton J.L.
      • Erskine R.J.
      Immunity and mastitis. Some new ideas for an old disease.
      ;
      • Lehmann M.
      • Wellnitz O.
      • Bruckmaier R.M.
      Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
      ;
      • Wellnitz O.
      • Arnold E.T.
      • Lehmann M.
      • Bruckmaier R.M.
      Short communication: Differential immunoglobulin transfer during mastitis challenge by pathogen-specific components.
      ). It has been shown that the loss of blood-milk barrier integrity enhances the IgG transfer into milk, which is much stronger in response to pathogen-associated molecular patterns from gram-negative than gram-positive bacteria (
      • Wellnitz O.
      • Arnold E.T.
      • Lehmann M.
      • Bruckmaier R.M.
      Short communication: Differential immunoglobulin transfer during mastitis challenge by pathogen-specific components.
      ). If specific IgG against the respective pathogen are present because of a previous infection or vaccination, these are considered to opsonize the pathogen and thus facilitate the elimination of the pathogen by the PMN (
      • Burton J.L.
      • Erskine R.J.
      Immunity and mastitis. Some new ideas for an old disease.
      ). Because of the differential transfer of IgG depending on the invaded pathogen, the increase of IgG in milk is also considered as a potential marker for pathogen-specific IMI detection if combined with SCC. However, the exact determination of IgG concentration in milk is complex and is therefore not suitable for automated and fast cow-side analyses in both conventional and automatic milking systems (AMS). Earlier studies indicate that other blood constituents such as lactate dehydrogenase (LDH) and serum albumin (SA) are transferred from blood into milk concomitantly with IgG; these constituents can be more easily analyzed and may serve as markers for the IgG transfer (
      • Lehmann M.
      • Wellnitz O.
      • Bruckmaier R.M.
      Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
      ;
      • Wall S.K.
      • Gross J.J.
      • Kessler E.C.
      • Villez K.
      • Bruckmaier R.M.
      Blood-derived proteins in milk at start of lactation: Indicators of active or passive transfer.
      ). Hence, the suitability of LDH and SA to indicate IgG transfer via the blood-milk barrier needs to be studied.
      An additional issue is the use of the ideal milk sample that allows the most sensitive detection of mastitis based on SCC, IgG, and related potential markers. All current strategies of mastitis detection in AMS are based on measuring various parameters in foremilk samples obtained after milk ejection (AME) whereas foremilk before milk ejection (BME) cannot be separated. This is due to immediately starting the teat-cleaning process that induces oxytocin release and milk ejection before foremilk is removed (
      • Mačuhová J.
      • Tančin V.
      • Bruckmaier R.M.
      Oxytocin release, milk ejection and milk removal in a multi-box automatic milking system.
      ;
      • Dzidic A.
      • Macuhova J.
      • Bruckmaier R.M.
      Effects of cleaning duration and water temperature on oxytocin release and milk removal in an automatic milking system.
      ). In our previous study, we have found that milk ejection decreases the power of mastitis detection by LDH (
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      ). Hence, the influence of milk ejection on SCC, IgG, and SA and the combination of these parameters for pathogen-specific IMI detection have not been tested before. Therefore, the objectives of the present study were to investigate (1) if LDH and SA could be used as alternative markers for IgG in automated mastitis detection and the benefit of the combined use of LDH or SA with SCC to differentiate gram-positive and gram-negative IMI; (2) the differences in in vivo gram-positive and gram-negative IMI for the SCC, IgG, LDH, and SA responses; and (3) the influences of milk ejection on SCC, IgG, LDH, and SA to identify the suitability of foremilk either BME or AME for more efficient automatic mastitis detection in AMS.

      MATERIALS AND METHODS

      All procedures involving animals were approved by the Animal Ethics Committee of the University of Sydney. The experiment was conducted at the University of Sydney ‘Corstorphine' pasture-based dairy farm located in Camden, NSW, Australia, with 85 ha of effective grazing land. The pasture land is mostly covered with annual ryegrass (Lolium multiflorum) oversown on kikuyu (Pennisetum clandestinum). A partial mixed ration containing primarily brewer's grain, orange pulp, and pasture silage (lucerne hay, oaten hay) and ∼7 kg of DM/cow grain-based commercial pelleted concentrate (18% protein) were provided as a supplement. The calving system was year-round and an automatic rotary system with a 24-unit platform and 5 robotic arms (DeLaval Automatic Milking Rotary) was used for milking the cows.

      Milk Sample Collection from Cows

      The selection criteria to collect milk samples from the cows used in this experiment are described in
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      . In brief, 48 Holstein-Friesian cows were selected for milk sampling based on their having a relatively high electrical conductivity (EC; ≥7.5 mS/cm, at milking temperature, 38°C) in any of the 4 quarters. They had an average of 4.34 ± 1.99 (mean ± SD) lactations and an average of 161 ± 117 DIM. Any cow with clinical mastitis in single or multiple quarters were treated with antibiotics after sampling. About 110 mL of foremilk samples was collected manually from individual quarters in 3 steps: 50 mL of strict foremilk before udder cleaning, taken within 60 s from all quarters; 50 mL of foremilk samples after cleaning the teats with warm water containing Iodophor LF12 solution; followed by an aseptic sample for culture after dipping the teats in Iodophor LF12 solution and cleaned with a 70% alcohol-soaked gauze (
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      ). Several cows were sampled on several days as a result of the EC-threshold criterion imposed (7 cows × 2 times, 5 cows × 3 times, 4 cows × 5 times, and 1 cow × 6 times). One cow had only 3 functional quarters and this resulted in 686 samples from 48 cows for the laboratory analysis and 343 samples for the culture test.

      Laboratory Analysis

      Milk samples collected in the first 2 steps (BME and AME) were analyzed for SCC, total IgG, LDH activity, and SA measurements. The procedures of SCC and LDH activity measurements are described in
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      .

      SCC

      The SCC was quantified using the principles of laser-based flow cytometry following the manufacturer's protocol (Bentley 2000 Instruments). A Bentley 2000 Autoanalyzer was used to measure SCC (cells × 1,000/mL).

      Total IgG Detection

      The procedure to detect total IgG is described in
      • Hernández-Castellano L.E.
      • Wall S.K.
      • Stephan R.
      • Corti S.
      • Bruckmaier R.
      Milk somatic cell count, lactate dehydrogenase activity, and immunoglobulin G 2 concentration associated with mastitis caused by different pathogens: A field study.
      . In brief, total IgG concentration was measured using a commercial ELISA kit specific for bovine IgG (Bethyl Laboratories). Samples were blocked in 5% fish skin gelatin (Sigma-Aldrich) diluted in double-distilled water. Samples were diluted in wash buffer (50 mM Tris, 0.14 M NaCl, 0.05% Tween 20, adjusted to pH 8.0) to ensure that samples were within the range of the standard curve. The standard curve was adjusted to 400, 300, 150, 75, 37.5, 18.75, and 9.375 ng/mL. Absorbance measurements were read on the Synergy Mx plate reader (Bio Tec Instruments). The inter- and intraassay coefficients of variation were ∼3 and ∼8%, respectively. The minimum detectable concentration was 9.375 ng/mL and samples were performed in duplicate.

      LDH Activity

      The LDH was measured by a commercial kit LDH IFCC (Axon Lan AG) and an automated analyzer (Cobas Mira, Roche Diagnostics) with minimum detectable activity 5 (U/L).

      SA Measurement

      The concentration of SA in milk samples was analyzed by ELISA using a bovine-specific commercial kit (Bethyl Laboratories) according to the manufacturer's instructions. Milk samples were diluted in wash buffer (50 mM Tris, 0.14 M NaCl, 0.05% Tween 20, adjusted to pH 8.0) to ensure the samples were in range of the standards. Absorbance measurements were read on the Synergy Mx plate reader (BioTek Instruments). The standard curve was 400, 200, 100, 50, 25, 12.5, and 6.25 (ng/mL), and the limit of detection was 6.25 ng/mL. The inter- and intraassay coefficients of variation were 4.58 and 6.99%, respectively. All analyses were performed in duplicate.

      Bacteriological Culture

      The aseptic milk samples collected on the third step were used for bacteriological culture (
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      ). Briefly, gram-positive species were identified by their growth in sheep blood agar, Enterococcal agar, Rambach agar, bile esculin test, the Christie–Atkins–Munch-Petersen (CAMP) test, and the leukocyte alkaline phosphatase test. Gram-negative coliforms were identified by their growth on sheep blood agar, MacConkey's agar, Gram staining, and potassium hydroxide (slime) test. Samples that did not yield microbial growth following 48 h of incubation were classified as control (no growth). Isolation of 2 or more bacteria genera from the same sample was considered as nondiagnostic or mixed.

      Statistical Analysis

      Data were analyzed using ASReml-R version 4 (
      • Butler D.G.
      • Cullis B.R.
      • Gilmour A.R.
      • Gogel B.G.
      • Thompson R.
      ASReml-R Reference Manual Version 4.
      ) built under R version 3.5.2 (https://www.r-project.org), in the form of mixed models as described below. As the distributions of SCC, IgG, LDH, and SA were positively skewed, they were log (base e) transformed before analysis to stabilize the variance and achieve normality of the outcome variables, or to reduce the leverage of very large values when used as explanatory variables. In total, 662 quarter observations were included in the analysis but due to insufficient volume to measure SCC in 24 milk samples, those observations were not possible to include in SCC analysis.

      Prediction of IgG Transmission from Blood into Milk by LDH and SA

      To compare the IgG transmission from blood into milk the following linear mixed model (LMM) was used:
      ln(IgG) = β0 + β1ln(LDH) + β2ln(SA) + s[ln(LDH)] + uC + uCQ + ε,
      [1]


      where ln(IgG) = log (base e) transformed IgG; ln(LDH) = log (base e) transformed lactate dehydrogenase as fixed effects together with nonlinear spline term (
      • Verbyla A.P.
      • Cullis B.R.
      • Kenward M.G.
      • Welham S.J.
      The analysis of designed experiments and longitudinal data by using smoothing splines.
      ) s[ln(LDH)]; ln(SA) = log (base e) transformed serum albumin as fixed effect (nonlinear spline was not significant); β0 is the constant; β1 and β2 specify the overall linear effects of the variables; uC and uCQ are cow and quarter nested within cow as random effects; and ε is a random error.

      Comparison of LDH Activity and SA to Predict Gram-Positive and Gram-Negative IMI

      To compare the effectiveness of LDH activity and SA to predict gram-positive IMI versus control, and gram-negative IMI versus control, both of these variables after log-transformation were standardized or re-scaled; that is, x'=(x-x¯)/SDx to compare variables across different scales, using the scale function in R. The standardized LDH activity and SA were assessed individually, as well as separate models for LDH and SA to predict gram-positive versus gram-negative IMI (i.e., 4 models). The following logistic generalized linear mixed model (GLMM) was fitted to the IMI status data, incorporating a spline function to allow for a possible nonlinear response (on the log-odds scale) of the explanatory variable:
      ln[π/(1 − π)] = β0 + β1x′ + s(x′) + uC + uCQ,
      [2]


      where π is the probability of mastitis (gram-positive IMI vs. control, or gram-negative IMI vs. control) and x′ is the standardized ln(LDH) and ln(SA). Each is included as a fixed linear effect, together with a nonlinear spline term, s(x′), specified as a random effect in the model, and also included uC and uCQ as cow and quarter nested within cow as random effects. Rescaling allowed displaying 2 predictors and 2 response variables with the corresponding fitted values on the same plot, allowing comparison of their effects.

      Assessment of Combined SCC-LDH and SCC-SA to Predict Gram-Positive and Gram-Negative IMI by Receiver Operating Characteristic Curve

      The following logistic GLMM were fitted separately to the gram-positive and gram-negative IMI status data with model terms as described in the previous model:
      ln[π/(1-π)]=β0+β1ln(SCC)+β2ln(LDH)+uC+uCQ,


      and
      ln[π/(1-π)]=β0+β1ln(SCC)+β2ln(SA)+uC+uCQ,
      [3]


      where π is the probability that a particular quarter is gram-positive or gram-negative IMI versus control; ln(SCC), ln(LDH), and ln(SA) are predictors, as a fixed effect; and uC and uCQ are the random cow and quarter nested within cow as random effects. Using output from these fitted models, construction of receiver operating characteristic (ROC) curves was then performed using the pROC package in R (version 4.0.2) (
      • Robin X.
      • Turck N.
      • Hainard A.
      • Tiberti N.
      • Lisacek F.
      • Sanchez J.
      • Mueller M.
      pROC: An open-source package for R and S+ to analyze and compare ROC curves.
      ) as described by
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      . In brief, the ROC assessment graphically illustrates the diagnostic test as a plot of sensitivity (Se) versus the complement of specificity (Sp; 1 − Sp) for varying cut points (
      • Hanley J.A.
      • McNeil B.J.
      The meaning and use of the area under a receiver operating characteristic (ROC) curve.
      ). The generated area under the curve (AUC) value from the ROC curve is used to measure diagnostic test performance, classified as excellent (0.9 to 1), good (0.8 to 0.9), fair (0.7 to 0.8), poor (0.6 to 0.7), or fail (0.5 to 0.6,
      • Swets J.A.
      Measuring the accuracy of diagnostic systems.
      ). The AUC values of SCC-LDH versus SCC-SA for gram-positive and gram-negative IMI were compared by using the roc.test (DeLong's test for 2 correlated ROC curves) in the pROC package. Test performance was also evaluated using Youden's index (J = Se + Sp − 1), selecting a cutoff point at which the index is maximized (
      • Ruopp M.D.
      • Perkins N.J.
      • Whitcomb B.W.
      • Schisterman E.F.
      Youden Index and optimal cut-point estimated from observations affected by a lower limit of detection.
      ).

      Correlation

      To estimate the correlations between pairs of variables (on the logarithmic scale) adjusting for the fixed effects, the following bivariate linear mixed models were fitted to all possible 6 pairs of variables (SCC vs. IgG, SCC vs. LDH, SCC vs. SA, IgG vs. LDH, IgG vs. SA, and LDH vs. SA):
      (y1y2)=(μ1μ2)+(Culture1Culture2)+(Cow1Cow2)+(Cow·Quarter1Cow·Quarter2)+(ε1ε2),
      [4]


      where y1 and y2 are the pair of variables as listed above (log-transformed), μ1 and μ2 are the overall trait means, Culture1 and Culture2 are the fixed effects for each variable of the gram-positive and gram-negative bacteria (e.g., Aerococcus sp., Bacillus, CNS, coagulase-positive Staphylococcus, coliform, Corynebacterium sp., Enterococcus faecalis, environmental Streptococcus, Streptococcus agalactiae, Streptococcus dysgalactiae, Streptococcus uberis, and Trueperella pyogenes) for variables 1 and 2, with random effects Cow1, Cow2, Cow·Quarter1, Cow·Quarter2, and random errors ε1 and ε2, with correlations also between the pairs of random effects. This is specified through the following models for the variance-covariance structure, where it is assumed that
      (Cow1Cow2)~N((00),(σC12σC12σC122σC22)),


      (Cow·Quarter1Cow·Quarter2)~N((00),(σCQ12σCQ12σCQ122σCQ22)),and


      (ε1ε2)~N((00),(σε12σε12σε122σε22)).


      Using the REML estimates of variance and covariance components, the correlation between a pair of traits (y1, y2), adjusting for the fixed effects, is calculated as follows:
      r12=cov(y1,y2)[var(y1)×var(y2)12]=σ^C12+σ^CQ12+σ^ε12[(σ^C1+σ^CQ12+σ^ε12)×(σ^C2+σ^CQ22+σ^ε22)]12.


      In the 3 covariance matrices, the diagonal terms represent the variances (σ2) for each trait, whereas the off-diagonal terms represent the covariances (σ) between the traits. Model fitting was conducted using ASReml R version 4 (
      • Butler D.G.
      • Cullis B.R.
      • Gilmour A.R.
      • Gogel B.G.
      • Thompson R.
      ASReml-R Reference Manual Version 4.
      ). From these pairwise correlations, a correlation matrix of all 4 traits was constructed and tested for validity (positive definite form). (Note that it was not possible to fit a multivariate model with 4 traits simultaneously in ASReml-R).

      Difference in SCC, IgG, LDH, and SA Between Gram-Positive and Gram-Negative IMI

      The differences in gram-positive and gram-negative IMI were assessed using the following linear mixed model:
      Y = β0 + culture + DIM + culture × DIM + uC + uCQ + ε,
      [5]


      where Y is the response variable ln(SCC), ln(IgG), ln(LDH), or ln(SA); culture = gram-positive, gram-negative IMI, and control as fixed effect; DIM = early (0 to 100 d), mid (101 to 200 d), and late (>200 d) lactation period as a fixed effect, with culture × DIM being the interaction; uC and uCQ are the cow and quarter nested within cow as random effects, and ε is the random error. Predicted means were calculated with corresponding standard error.

      Principal Component Analysis

      To explore all 4 milk parameters [i.e., ln(SCC), ln(IgG), ln(LDH), and ln(SA)] simultaneously, a principal component analysis (PCA) was undertaken on these data (variables scaled to have a unit standard deviation), using the prcomp function in R. Loadings of the first 2 principal components were interpreted.

      Principal Component Scores Comparing Gram-Positive and Gram-Negative IMI

      The following linear mixed model was used to predict mastitis types based on principal component (PC) scores:
      Y = β0 + culture + (βDIM + γCulture.DIM)DIM + s(DIM) + uC + uCQ + ε,
      [6]


      where Y is the response variable, either PC1 scores or PC2 scores; DIM = days in milk as linear fixed effect together with spline term s(DIM) to allow for a nonlinear response over the lactation; βDIM specifies the overall linear effect of DIM; γCulture.DIM is used to specify the culture × DIM interaction; and other terms are as defined in previous models. Predicted means were calculated with corresponding standard error.

      Effect of Milk Ejection on SCC, IgG, LDH Activity, and SA

      The differences in 2 different milk samples collected BME versus AME were assessed using the following linear mixed model:
      Y = β0 + time + βDIMDIM + s(DIM) + uC + uCQ + ε,
      [7]


      where Y is the response variable [ln(SCC), ln(IgG), ln(LDH), ln(SA)], time = BME versus AME as a fixed effect, and other terms are as defined previously. Predicted means were calculated with corresponding standard error.

      RESULTS

      Bacteriological Culture and Electrical Conductivity

      The gold standard of this study was bacteriological culture. Out of 343 bacteriological culture tests, 169 quarters were culture positive (e.g., gram-positive, gram-negative, and mixed growth) and 174 quarters had no growth (control). There were 157 (45.77%) gram-positive quarter samples (n = 39 cows), 6 (1.75%) gram-negative quarter samples (n = 4 cows), and 174 (50.73%) control quarter samples (n = 41 cows). Each culture was used as a gold standard for individual quarter milk samples collected in 2 steps, resulting in a doubling of the number of the gram-positive, gram-negative, and control samples in the analysis. Among the gram-positive, gram-negative, and control quarters, 121, 3, and 21 quarters had >7.5 mS/cm EC and 35, 3, and 145 quarters had <7.5 mS/cm EC, respectively.

      Prediction of IgG Transmission from Blood into Milk by LDH and SA

      Overall, both SA (P < 0.001) and LDH (P = 0.001) were positively associated with IgG transmission from blood into milk based on combined BME and AME milk samples from the LMM (Eq. [1]; Table 1). The positive association of IgG with LDH was nonlinear (Figure 1), whereas the association was linear with SA (on log scale, Figure 2). Only SA had significant (P < 0.001) positive associations with IgG transmission based on separate BME and AME samples.
      Table 1Comparison of lactate dehydrogenase (LDH) and serum albumin (SA) to predict IgG transmission from blood into milk by linear mixed model
      Linear mixed models to predict IgG including LDH and SA, with random effect estimates for each cow and cow-quarter.
      Item
      BME = sampling times before ejection; AME = sampling times after ejection.
      Coefficient ± SEP-valueCow (σ2)
      σ2 = variance estimate ± SE.
      Quarter (σ2)
      σ2 = variance estimate ± SE.
      BME + AME (n = 662)
       LDH (log base e, U/L)0.05 ± 0.010.0010.18 ± 0.040.02 ± 0.00
       SA (log base e, mg/mL)0.46 ± 0.03<0.001
      BME (n = 331)
       LDH (log base e, U/L)0.02 ± 0.020.4020.19 ± 0.040.02 ± 0.00
       SA (log base e, mg/mL)0.50 ± 0.03<0.001
      AME (n = 331)
       LDH (log base e, U/L)0.09 ± 0.020.0050.17 ± 0.040.06 ± 0.04
       SA (log base e, mg/mL)0.42 ± 0.04<0.001
      1 Linear mixed models to predict IgG including LDH and SA, with random effect estimates for each cow and cow-quarter.
      2 BME = sampling times before ejection; AME = sampling times after ejection.
      3 σ2 = variance estimate ± SE.
      Figure thumbnail gr1
      Figure 1Prediction of IgG (mg/mL, estimated value ± SE) transmission by lactate dehydrogenase (LDH, U/L, plotted on a logarithmic x-axis scale). This mixed model included ln(LDH) and ln(serum albumin, SA) as fixed effects together with nonlinear term ln(LDH), with cow and quarter nested within cow as random effects (n = 662 quarters). Because of positively skewed distributions, IgG, LDH, and SA data were log-transformed.
      Figure thumbnail gr2
      Figure 2Prediction of IgG (mg/mL, estimated value ± SE) transmission by serum albumin (SA, mg/mL). This mixed model included ln(lactate dehydrogenase, LDH) and ln(SA) as fixed effects together with nonlinear term ln(LDH), with cow and quarter nested within cow as random effects (n = 662 quarters). Because of positively skewed distributions, IgG, LDH, and SA data were log-transformed.

      Comparison of LDH Activity and SA to Predict Gram-Positive and Gram-Negative IMI

      After rescaling (standardizing), ln(LDH) activity and ln(SA) showed positive associations with the probability of gram-positive IMI versus control and gram-negative IMI versus control (Figure 3) based on the GLMM in Eq. [2]. The standardized ln(SA) (P = 0.001) and standardized ln(LDH) (P = 0.008) had similar gram-negative IMI prediction ability. However, standardized ln(SA) (P = 0.008) had better gram-positive IMI prediction ability than standardized ln(LDH) (P = 0.058).
      Figure thumbnail gr3
      Figure 3Comparison of probabilities of mastitis (predicted value ± SE) for varying lactate dehydrogenase (LDH, solid lines) and serum albumin (SA, dashed lines), with separate predictions for gram-positive (n = 304), gram-negative IMI (n = 12) and control (n = 334) quarters. Each of the 4 logistic models included standardized LDH or standardized SA (expressed as number of SD away from the mean, after log-transformation) as fixed effects, with cow or quarter nested within cow as random effects, and with a model for gram-positive versus control, and gram-negative versus control. Because of positively skewed distributions, LDH and SA data were log-transformed.

      Assessment of Individual and Combined SCC-LDH and SCC-SA to Predict Gram-Positive and Gram-Negative IMI by ROC Curve

      The combined SCC-LDH and SCC-SA (GLMM in Eq. [3]) had significantly greater IMI prediction ability than individual ln(SCC), ln(LDH) and ln(SA) separately (Table 2). In the ROC evaluation, we obtained excellent (AUC >0.9) gram-positive versus control and gram-negative versus control IMI prediction ability for both SCC-LDH and SCC-SA. The SCC-LDH vs. SCC-SA had lower gram-positive (AUC = 0.984 vs. 0.986, P = 0.010) but similar gram-negative (AUC = 0.995 vs. 0.998, P = 0.134) IMI prediction ability, as specified by the AUC values, respectively.
      Table 2Analysis of receiver operating characteristic curves, sensitivity (Se), and specificity (Sp) at optimum cutoff value for prediction of gram-positive and gram-negative IMI based on model (logistic generalized linear mixed model
      Logistic generalized linear mixed models included 2 variables (SCC, and lactate dehydrogenase activity or serum albumin), with cow and cow-quarter random effects. In both gram-positive and gram-negative mastitis, the control was quarters without any bacterial growth (n = 329).
      )
      Predictor
      LDH = lactate dehydrogenase; SA = serum albumin.
      Culture
      Culture: gram-negative = IMI caused by gram-negative coliform bacteria (n = 10); gram-positive = IMI caused by gram-positive bacteria (n = 295).
      AUC
      AUC = area under the curve.
      (95% CI)
      SeSpP-value
      SCC-LDHGram-positive0.984 (0.978–0.991)0.9630.9180.0105
      Gram-negative0.995 (0.990–1.000)1.0000.9880.134
      Comparison of AUC values with combined SCC-SA.
      SCC-SAGram-positive0.986 (0.980–0.992)0.9630.924
      Gram-negative0.998 (0.994–1.000)1.0000.994
      SCCGram-positive0.649 (0.606–0.692)0.7120.568<0.001
      Comparison of AUC values with combined SCC-LDH.
      Gram-negative0.806 (0.655–0.957)0.7000.8270.013
      Comparison of AUC values with combined SCC-LDH.
      LDHGram-positive0.608 (0.564–0.652)0.6070.565<0.001
      Comparison of AUC values with combined SCC-LDH.
      Gram-negative0.883 (0.773–0.993)0.8000.8480.046
      Comparison of AUC values with combined SCC-LDH.
      SAGram-positive0.608 (0.564–0.652)0.4270.726<0.001
      Comparison of AUC values with combined SCC-SA.
      Gram-negative0.849 (0.708–0.990)0.6000.9880.037
      Comparison of AUC values with combined SCC-SA.
      1 Logistic generalized linear mixed models included 2 variables (SCC, and lactate dehydrogenase activity or serum albumin), with cow and cow-quarter random effects. In both gram-positive and gram-negative mastitis, the control was quarters without any bacterial growth (n = 329).
      2 LDH = lactate dehydrogenase; SA = serum albumin.
      3 Culture: gram-negative = IMI caused by gram-negative coliform bacteria (n = 10); gram-positive = IMI caused by gram-positive bacteria (n = 295).
      4 AUC = area under the curve.
      5 Comparison of AUC values with combined SCC-SA.
      6 Comparison of AUC values with combined SCC-LDH.

      Correlations

      In bivariate mixed model analysis, after adjusting for fixed effects (Eq. [4]), there were very strong positive correlations (all calculated on a log scale) between LDH activity with SCC (r = 0.80, P < 0.001) followed by LDH with SA (r = 0.70, P < 0.001; Table 3).
      Table 3Model-based correlations (± SE) between SCC, IgG, lactate dehydrogenase (LDH), and serum albumin (SA) estimated from bivariate linear mixed model (n = 642 quarters)
      Itemln(SCC)ln(IgG)ln(LDH)ln(SA)
      ln(SCC)1
      ln(IgG)0.39 ± 0.07
      P < 0.001 for all correlation coefficients shown.
      1
      ln(LDH)0.80 ± 0.02
      P < 0.001 for all correlation coefficients shown.
      0.53 ± 0.05
      P < 0.001 for all correlation coefficients shown.
      1
      ln(SA)0.47 ± 0.05
      P < 0.001 for all correlation coefficients shown.
      0.62 ± 0.06
      P < 0.001 for all correlation coefficients shown.
      0.70 ± 0.03
      P < 0.001 for all correlation coefficients shown.
      1
      *** P < 0.001 for all correlation coefficients shown.

      Difference in SCC, IgG, LDH, and SA Between Gram-Positive and Gram-Negative IMI

      Based on the LMM in Eq. [5], the mean SCC, IgG, LDH, and SA between gram-positive and gram-negative IMI in early (0–100 DIM), mid (101–200 DIM), and late (>200 DIM) lactation period are presented in Table 4. Overall., the differences in SCC, IgG, LDH, and SA in gram-positive and gram-negative IMI type were significant in the early- to mid-lactation period. There was a significant interaction between DIM and pathogen type for SCC (P < 0.001), IgG (P < 0.001), LDH activity (P < 0.001), and SA (P < 0.001). The ranges of DIM included in the analysis for different IMI type are presented in Table 5.
      Table 4Model-based mean (± SE) SCC, total IgG, lactate dehydrogenase (LDH), and serum albumin (SA) in milk samples collected before (BME) and after (AME) ejection times from quarters infected by gram-positive and gram-negative (coliform) bacteria by linear mixed model
      Linear mixed models with the random effects for cow and cow-quarter.
      Item
      Early = lactation stages within 0 to 100 DIM; mid = lactation period within 101 to 200 DIM; late = lactation period >200 DIM.
      Time
      Time: BME = sampling times before ejection; AME = sampling times after ejection.
      and culture
      Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; control = quarters without any bacterial growth. In the case of SCC, the total number of observations was 634 (gram-positive: 295, gram-negative: 10, control: 329) due to lack of SCC information in a few cows.
      BME + AMEBMEAME
      Gram-positive (n = 304)Gram-negative (n = 12)Control (n = 334)Gram-positive (n = 152)Gram-negative (n = 6)Control (n = 167)Gram-positive (n = 152)Gram-negative (n = 6)Control (n = 167)
      SCC (loge cells/mL)
       Early5.99 ± 0.29
      Means within a row with different superscripts differ significantly (P < 0.05).
      7.78 ± 0.66
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.55 ± 0.26
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.26 ± 0.32
      Means within a row with different superscripts differ significantly (P < 0.05).
      8.12 ± 1.4
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.56 ± 0.28
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.99 ± 0.32
      Means within a row with different superscripts differ significantly (P < 0.05).
      8.05 ± 0.82
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.30 ± 0.27
      Means within a row with different superscripts differ significantly (P < 0.05).
       Mid6.39 ± 0.30
      Means within a row with different superscripts differ significantly (P < 0.05).
      4.43 ± 0.76
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.52 ± 0.28
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.91 ± 0.35
      Means within a row with different superscripts differ significantly (P < 0.05).
      2.72 ± 1.22
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.61 ± 0.33
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.18 ± 0.34
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.60 ± 1.15
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.22 ± 0.33
      Means within a row with different superscripts differ significantly (P < 0.05).
       Late6.33 ± 0.27
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.10 ± 0.58
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.10 ± 0.25
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.70 ± 0.29
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.49 ± 0.91
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.46 ± 0.28
      Means within a row with different superscripts differ significantly (P < 0.05).
      6.12 ± 0.29
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.98 ± 0.85
      Means within a row with different superscripts differ significantly (P < 0.05).
      5.74 ± 0.28
      Means within a row with different superscripts differ significantly (P < 0.05).
      IgG (mg/mL)
       Early1.08 ± 0.13
      Means within a row with different superscripts differ significantly (P < 0.05).
      3.02 ± 0.55
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.92 ± 0.11
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.10 ± 0.15
      Means within a row with different superscripts differ significantly (P < 0.05).
      2.89 ± 0.67
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.00 ± 0.13
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.02 ± 0.13
      Means within a row with different superscripts differ significantly (P < 0.05).
      3.09 ± 0.71
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.85 ± 0.10
      Means within a row with different superscripts differ significantly (P < 0.05).
       Mid0.92 ± 0.11
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.52 ± 0.33
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.89 ± 0.10
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.94 ± 0.13
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.52 ± 0.45
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.90 ± 0.12
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.87 ± 0.11
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.66 ± 0.51
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.84 ± 0.11
      Means within a row with different superscripts differ significantly (P < 0.05).
       Late0.82 ± 0.09
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.72 ± 0.13
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.72 ± 0.08
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.90 ± 0.12
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.81 ± 0.19
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.80 ± 0.10
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.77 ± 0.10
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.63 ± 0.15
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.70 ± 0.08
      Means within a row with different superscripts differ significantly (P < 0.05).
      LDH activity (U/L)
       Early261 ± 50
      Means within a row with different superscripts differ significantly (P < 0.05).
      3,482 ± 1,525
      Means within a row with different superscripts differ significantly (P < 0.05).
      141 ± 25
      Means within a row with different superscripts differ significantly (P < 0.05).
      326 ± 75
      Means within a row with different superscripts differ significantly (P < 0.05).
      4,676 ± 3,068
      Means within a row with different superscripts differ significantly (P < 0.05).
      149 ± 30
      Means within a row with different superscripts differ significantly (P < 0.05).
      245 ± 48
      Means within a row with different superscripts differ significantly (P < 0.05).
      3,799 ± 2,107
      Means within a row with different superscripts differ significantly (P < 0.05).
      115 ± 20
      Means within a row with different superscripts differ significantly (P < 0.05).
       Mid208 ± 43
      Means within a row with different superscripts differ significantly (P < 0.05).
      622 ± 357
      Means within a row with different superscripts differ significantly (P < 0.05).
      177 ± 35
      Means within a row with different superscripts differ significantly (P < 0.05).
      261 ± 68
      Means within a row with different superscripts differ significantly (P < 0.05).
      761 ± 725
      Means within a row with different superscripts differ significantly (P < 0.05).
      213 ± 53
      Means within a row with different superscripts differ significantly (P < 0.05).
      170 ± 37
      Means within a row with different superscripts differ significantly (P < 0.05).
      871 ± 696
      Means within a row with different superscripts differ significantly (P < 0.05).
      141 ± 29
      Means within a row with different superscripts differ significantly (P < 0.05).
       Late177 ± 33
      Means within a row with different superscripts differ significantly (P < 0.05).
      189 ± 82
      Means within a row with different superscripts differ significantly (P < 0.05).
      171 ± 29
      Means within a row with different superscripts differ significantly (P < 0.05).
      212 ± 46
      Means within a row with different superscripts differ significantly (P < 0.05).
      218 ± 154
      Means within a row with different superscripts differ significantly (P < 0.05).
      221 ± 46
      Means within a row with different superscripts differ significantly (P < 0.05).
      153 ± 28
      Means within a row with different superscripts differ significantly (P < 0.05).
      178 ± 105
      Means within a row with different superscripts differ significantly (P < 0.05).
      137 ± 24
      Means within a row with different superscripts differ significantly (P < 0.05).
      SA (mg/mL)
       Early0.29 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.37 ± 0.30
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.19 ± 0.02
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.31 ± 0.04
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.37 ± 0.45
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.19 ± 0.02
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.28 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.46 ± 0.40
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.18 ± 0.02
      Means within a row with different superscripts differ significantly (P < 0.05).
       Mid0.32 ± 0.04
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.17 ± 0.32
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.27 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.30 ± 0.05
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.12 ± 0.50
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.27 ± 0.04
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.28 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      1.27 ± 0.48
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.24 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
       Late0.25 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.24 ± 0.05
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.23 ± 0.02
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.29 ± 0.04
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.28 ± 0.09
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.27 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.26 ± 0.03
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.23 ± 0.07
      Means within a row with different superscripts differ significantly (P < 0.05).
      0.23 ± 0.02
      Means within a row with different superscripts differ significantly (P < 0.05).
      a–c Means within a row with different superscripts differ significantly (P < 0.05).
      1 Linear mixed models with the random effects for cow and cow-quarter.
      2 Early = lactation stages within 0 to 100 DIM; mid = lactation period within 101 to 200 DIM; late = lactation period >200 DIM.
      3 Time: BME = sampling times before ejection; AME = sampling times after ejection.
      4 Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; control = quarters without any bacterial growth. In the case of SCC, the total number of observations was 634 (gram-positive: 295, gram-negative: 10, control: 329) due to lack of SCC information in a few cows.
      Table 5Descriptive statistics of DIM for gram-positive and gram-negative IMI compared with control
      Culture
      Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; control = quarters without any bacterial growth.
      DIM
      MinimumFirst quarterMedianMeanThird quarterMaximum
      Control256124144231398
      Gram-positive156179178273459
      Gram-negative19119126247260
      1 Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; control = quarters without any bacterial growth.
      Based on BME and AME together and AME samples, gram-negative IMI had significantly greater (P < 0.05) SCC than control in the early-lactation period. However, based on BME and AME together, and BME samples, gram-negative IMI had significantly lower SCC than control and gram-positive IMI in the mid-lactation period. The gram-negative IMI had significantly (P < 0.05) greater IgG than control and gram-positive IMI in the early- and mid-lactation period in all sample types (e.g., BME and AME together, only BME and only AME).
      Based on BME and AME together, and only AME samples, the LDH was significantly (P < 0.05) greater in gram-negative IMI than control and gram-positive IMI in the early- and mid-lactation period. Based on only BME, gram-negative IMI had significantly (P < 0.05) greater LDH than control and gram-positive IMI in the early-lactation period.
      Based on BME and AME together, the SA was significantly (P < 0.05) greater in gram-negative IMI than control in early- and mid-lactation period, also than gram-positive IMI in early-lactation period. Gram-negative IMI had significantly greater SA than control and gram-positive IMI in early- and mid-lactation period when comparing the only BME and only AME samples.

      PCA

      Exploring associations between all 4 parameters, based on the PCA, 90% of the variation in the 4-dimensional space was accounted for by the first 2 PC (PC1 and PC2; Table 6). In PC1, all 4 variables had approximately equal positive loadings (∼0.5) and is effectively an average of these 4 parameters. In PC2, there was a contrast between a higher impact positive SCC loading versus negative IgG loading and a lesser impact positive LDH activity versus negative SA loading. From the plot of PC1 versus PC2, there is evidence of clustering of results on types of IMI (Figure 4).
      Table 6Variation explained by principal component (PC) analysis and its component loadings (n = 634)
      ItemPC1PC2
      Variation explained0.740.15
      Variable
       SCC (loge, cells/mL)0.490.61
       IgG (loge, mg/mL)0.45−0.69
       Lactate dehydrogenase (loge, U/L)0.540.28
       Serum albumin (loge, mg/mL)0.51−0.27
      Figure thumbnail gr4
      Figure 4Plots of the first 2 principal component scores (PC1, PC2) for gram-positive (n = 295) and gram-negative (n = 10) IMI and control (n = 329) quarters. Each infection type is indicated by the different point shapes.

      Association Between PC Scores and Gram-Positive and Gram-Negative IMI

      Based on the analysis of PC1 scores using the LMM (Eq. [6]), gram-negative IMI had significantly (P < 0.05) greater scores than control and gram-positive IMI (Table 7). Given the loadings of PC1, this indicates that gram-negative IMI samples tend to have higher values of SCC, IgG, LDH, and SA, and hence are useful to predict gram-negative IMI. Based on PC2 scores, gram-negative IMI also had significantly (P < 0.05) lower scores than control and gram-positive IMI. Hence, samples with lower PC2 scores (i.e., greater IgG but lower SCC) were more likely to be associated with coliform IMI. On the other hand, higher PC2 scores (i.e., greater SCC but lower IgG) were more likely in gram-positive IMI samples.
      Table 7Model-based means (± SE) of principal component scores (PC1, PC2) from a linear mixed model
      Linear mixed models with the random effects for cow and cow-quarter.
      comparing IMI and control groups
      The principal component analysis was conducted on the 4 variables SCC, IgG, lactate dehydrogenase, and serum albumin responses.
      Culture
      Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; Control = quarters without any bacterial growth.
      PC1 scorePC2 score
      Control (n = 329)−0.01 ± 0.35
      Means within a column with different superscripts differ significantly (P < 0.05).
      0.01 ± 0.15
      Means within a column with different superscripts differ significantly (P < 0.05).
      Gram-positive (n = 295)0.29 ± 0.36
      Means within a column with different superscripts differ significantly (P < 0.05).
      −0.03 ± 0.15
      Means within a column with different superscripts differ significantly (P < 0.05).
      Gram-negative (n = 10)1.12 ± 0.46
      Means within a column with different superscripts differ significantly (P < 0.05).
      −0.65 ± 0.21
      Means within a column with different superscripts differ significantly (P < 0.05).
      a,b Means within a column with different superscripts differ significantly (P < 0.05).
      1 Linear mixed models with the random effects for cow and cow-quarter.
      2 The principal component analysis was conducted on the 4 variables SCC, IgG, lactate dehydrogenase, and serum albumin responses.
      3 Culture: gram-negative = IMI caused by gram-negative coliform bacteria; gram-positive = IMI caused by gram-positive bacteria; Control = quarters without any bacterial growth.

      Effect of Milk Ejection on SCC, IgG, LDH Activity, and SA

      Based in the LMM specified by Eq. [7], the milk ejection had significant effects on SCC, IgG, LDH activity, and SA in mastitic quarters caused my gram-positive and gram-negative bacteria as well as in control (Table 8). Overall, the mastitic quarters had greater SCC, IgG, LDH, and SA than the healthy quarters, particularly in BME samples. When comparing without separating mastitic and healthy quarters, the BME samples had significantly greater SCC, LDH, IgG, and SA than the AME samples (all P < 0.001). In comparison, considering only mastitic quarters, BME samples had significantly greater SCC (P < 0.001), IgG (P < 0.001), LDH activity (P < 0.001), and SA (P = 0.010) than AME samples. Similar comparison in control quarters also resulted in significant greater SCC (P < 0.001), IgG (P < 0.001), LDH activity (P < 0.001), and SA (P = 0.020) in BME than AME samples.
      Table 8Model-based means (± SE) of SCC, total IgG, lactate dehydrogenase (LDH), and serum albumin (SA) in milk samples before and after ejection from linear mixed models
      Linear mixed models to calculate predicted means of 4 outcome variables (SCC, IgG, lactate dehydrogenase activity, and serum albumin), with random effect estimates for each cow and cow-quarter.
      Culture
      Culture: quarters with IMI: IMI caused by gram-positive or gram-negative bacteria; control = quarters without any bacterial growth.
      Time
      BME = sampling times before ejection; AME = sampling times after ejection.
      P-value
      BMEAME
      Quarters with IMI and control (n = 662)
       SCC (loge cells/mL)6.31 ± 0.175.89 ± 0.17<0.001
       IgG (mg/mL)0.96 ± 0.080.84 ± 0.07<0.001
       LDH (U/L)226 ± 27159 ± 19<0.001
       SA (mg/mL)0.27 ± 0.020.24 ± 0.02<0.001
      Quarters with IMI (n = 328)
       SCC (loge cells/mL)6.77 ± 0.196.34 ± 0.19<0.001
       IgG (mg/mL)1.02 ± 0.110.91 ± 0.10<0.001
       LDH (U/L)313 ± 54211 ± 36<0.001
       SA (mg/mL)0.33 ± 0.030.29 ± 0.030.010
      Control (n = 334)
       SCC (loge cells/mL)6.04 ± 0.215.61 ± 0.21<0.001
       IgG (mg/mL)0.91 ± 0.080.80 ± 0.07<0.001
       LDH (U/L)178 ± 20130 ± 14<0.001
       SA (mg/mL)0.23 ± 0.020.21 ± 0.020.002
      1 Linear mixed models to calculate predicted means of 4 outcome variables (SCC, IgG, lactate dehydrogenase activity, and serum albumin), with random effect estimates for each cow and cow-quarter.
      2 Culture: quarters with IMI: IMI caused by gram-positive or gram-negative bacteria; control = quarters without any bacterial growth.
      3 BME = sampling times before ejection; AME = sampling times after ejection.

      DISCUSSION

      To our knowledge this is the first study that has investigated the suitability of SA concentration and LDH activity in milk to be used as potential markers of IgG transmission from blood into milk in vivo in naturally occurring mastitis, and based on that, to distinguish between pathogen types that cause mastitis. Both SA and, in particular, LDH can be analyzed with established laboratory methods, which would allow a potential automatization. In contrast, the exact measurement of IgG by ELISA appears to be more complicated, mostly requiring different steps of sample dilution, and an automatization for on-farm use is likely not possible (
      • Lehmann M.
      • Wellnitz O.
      • Bruckmaier R.M.
      Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
      ;
      • Wall S.K.
      • Wellnitz O.
      • Hernández-Castellano L.E.
      • Ahmadpour A.
      • Bruckmaier R.M.
      Supraphysiological oxytocin increases the transfer of immunoglobulins and other blood components to milk during lipopolysaccharide- and lipoteichoic acid–induced mastitis in dairy cows.
      ). Our study has found that both LDH and SA could be markers of IgG transfer via the blood-milk barrier, and both are suitable to differentiate gram-positive and gram-negative IMI in combination with the SCC results in milk samples. Further, foremilk BME is more informative for better mastitis detection in AMS than AME.
      In this study, prediction of IgG using LDH supports previous findings (in experimentally challenged mammary glands) that LDH could be a marker for IgG transfer via the leaky blood-milk barrier in vivo in naturally occurring mastitis (
      • Lehmann M.
      • Wellnitz O.
      • Bruckmaier R.M.
      Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
      ). The reason for investigating SA in this in vivo study is that SA appears earlier than LDH (2 vs. 5 h) in experimentally induced mastitis as reported in other studies (
      • Wall S.K.
      • Wellnitz O.
      • Hernández-Castellano L.E.
      • Ahmadpour A.
      • Bruckmaier R.M.
      Supraphysiological oxytocin increases the transfer of immunoglobulins and other blood components to milk during lipopolysaccharide- and lipoteichoic acid–induced mastitis in dairy cows.
      ;
      • Wellnitz O.
      • Zbinden C.
      • Huang X.
      • Bruckmaier R.M.
      Short communication: Differential loss of bovine mammary epithelial barrier integrity in response to lipopolysaccharide and lipoteichoic acid.
      ). In addition to the transfer via the blood-milk barrier, both SA and LDH may originate from mammary epithelial cell production that may cause an increase of SA (
      • Shamay A.
      • Homans R.
      • Fuerman Y.
      • Levin I.
      • Barash H.
      • Silanikove N.
      • Mabjeesh S.J.
      Expression of albumin in nonhepatic tissues and its synthesis by the bovine mammary gland.
      ), and mammary cell damage may contribute to the increased concentration of LDH in milk (
      • Wellnitz O.
      • Zbinden C.
      • Huang X.
      • Bruckmaier R.M.
      Short communication: Differential loss of bovine mammary epithelial barrier integrity in response to lipopolysaccharide and lipoteichoic acid.
      ;
      • Cheng J.
      • Zhang J.
      • Yang J.
      • Yi B.
      • Liu G.
      • Zhou M.
      • Kastelic J.P.
      • Han B.
      • Gao J.
      Klebsiella pneumoniae infection causes mitochondrial damage and dysfunction in bovine mammary epithelial cells.
      ). The potential contribution of the 2 sources of SA and LDH may support the possibility to distinguish between IMI-causing pathogen types. The better gram-positive IMI prediction ability of SA than LDH further specifies the utility of SA to differentiate the mastitis type. Our findings demonstrate that SA could be used as an alternative marker to LDH to indicate IgG transfer via the blood-milk barrier. Based on several studies, the combination of multiple sensor information is useful to detect and differentiate mastitis types more efficiently in AMS (
      • Hogeveen H.
      • Kamphuis C.
      • Steeneveld W.
      • Mollenhorst H.
      Sensors and clinical mastitis-the quest for the perfect alert.
      ;
      • Steeneveld W.
      • Vernooij J.C.M.
      • Hogeveen H.
      Effect of sensor systems for cow management on milk production, somatic cell count, and reproduction.
      ;
      • Hernández-Castellano L.E.
      • Wall S.K.
      • Stephan R.
      • Corti S.
      • Bruckmaier R.
      Milk somatic cell count, lactate dehydrogenase activity, and immunoglobulin G 2 concentration associated with mastitis caused by different pathogens: A field study.
      ). In this study, the greater gram-positive and gram-negative IMI prediction ability of combined SCC-LDH and SCC-SA compared with SCC, LDH, and SA alone, supports previous findings (
      • Hernández-Castellano L.E.
      • Wall S.K.
      • Stephan R.
      • Corti S.
      • Bruckmaier R.
      Milk somatic cell count, lactate dehydrogenase activity, and immunoglobulin G 2 concentration associated with mastitis caused by different pathogens: A field study.
      ;
      • Khatun M.
      • Bruckmaier R.M.
      • Thomson P.C.
      • House J.
      • García S.C.
      Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
      ). Based on similar gram-positive and gram-negative IMI prediction ability by SCC-LDH and SCC-SA, our recommendation is SCC-SA could be the alternate of SCC-LDH to differentiate mastitis type inline in AMS.
      Overall, the greater SCC, IgG, LDH, and SA in gram-negative coliform IMI compared with gram-positive IMI might be due to severe damage of the blood-milk barrier allowing abundant migration as found previously in experimentally induced mastitis (
      • Wellnitz O.
      • Arnold E.T.
      • Bruckmaier R.M.
      Lipopolysaccharide and lipoteichoic acid induce different immune responses in the bovine mammary gland.
      ;
      • Lehmann M.
      • Wellnitz O.
      • Bruckmaier R.M.
      Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
      ;
      • Wall S.K.
      • Wellnitz O.
      • Hernández-Castellano L.E.
      • Ahmadpour A.
      • Bruckmaier R.M.
      Supraphysiological oxytocin increases the transfer of immunoglobulins and other blood components to milk during lipopolysaccharide- and lipoteichoic acid–induced mastitis in dairy cows.
      ). The greater SCC, IgG, LDH, and SA in coliform IMI compared with gram-positive IMI was also supported by the PCA where significantly greater PC1 scores (with approximately equal positive loadings of SCC, IgG, LDH, and SA) were found. The differences in SCC, IgG, LDH, and SA responses for gram-positive and gram-negative IMI might come from differential activation of the immune system by them (
      • Wellnitz O.
      • Arnold E.T.
      • Bruckmaier R.M.
      Lipopolysaccharide and lipoteichoic acid induce different immune responses in the bovine mammary gland.
      ). The greatest correlation of LDH with SCC (r = 0.80) indicates the concurrent stronger innate immune response and migration of the leucocytes into milk with greater damage of the mammary epithelial cells for LDH production (
      • Wagner S.A.
      • Jones D.E.
      • Apley M.D.
      Effect of endotoxic mastitis on epithelial cell numbers in the milk of dairy cows.
      ;
      • Wellnitz O.
      • Zbinden C.
      • Huang X.
      • Bruckmaier R.M.
      Short communication: Differential loss of bovine mammary epithelial barrier integrity in response to lipopolysaccharide and lipoteichoic acid.
      ;
      • Hernández-Castellano L.E.
      • Wall S.K.
      • Stephan R.
      • Corti S.
      • Bruckmaier R.
      Milk somatic cell count, lactate dehydrogenase activity, and immunoglobulin G 2 concentration associated with mastitis caused by different pathogens: A field study.
      ). Additionally, the LDH control mean was 5.26× greater than the IMI groups, whereas the corresponding ratio for IgG was 1.63× and for SA was 2.48×; these might be due to local LDH release from affected mammary epithelial cells, as IgG and SA originate exclusively from blood (
      • Stelwagen K.
      • Prosser C.G.
      • Davis S.R.
      • Farr V.C.
      • Politis I.
      • White J.H.
      • Zavizion B.
      Effect of milking frequency and somatotropin on the activity of plasminogen activator, plasminogen, and plasmin in bovine milk.
      ;
      • Wellnitz O.
      • Zbinden C.
      • Huang X.
      • Bruckmaier R.M.
      Short communication: Differential loss of bovine mammary epithelial barrier integrity in response to lipopolysaccharide and lipoteichoic acid.
      ).
      In mastitic quarters, greater SCC, IgG, LDH activity, and SA in BME samples compared with AME samples could be due to a protection mechanism against invading pathogens in the Furstenberg's rosette area of the teat (
      • Nickerson S.C.
      • Pankey J.W.
      Cytologic observations of the bovine teat end.
      ). Such differences in control (culture-negative healthy quarters with average SCC <120,000 cells/mL) might also be due to different pathways of transfer of SCC, IgG, LDH, and SA from blood into alveolar secretory tissue compared with teat tissue as reported for SCC previously (
      • Sarikaya H.
      • Werner-Misof C.
      • Atzkern M.
      • Bruckmaier R.M.
      Distribution of leucocyte populations, and milk composition, in milk fractions of healthy quarters in dairy cows.
      ). Moreover, in mastitic and control quarters, the difference in SCC, IgG, LDH, and SA, particularly in the early- and mid-lactation period, might be due to the unequal distribution of the blood-milk barrier, which has also been shown for other parameters such as electrolytes, electrical conductivity, with generally more pronounced changes of milk composition in BME foremilk (
      • Bruckmaier R.M.
      • Weiss D.
      • Wiedemann M.
      • Schmitz S.
      • Wendl G.
      Changes of physicochemical indicators during mastitis and the effects of milk ejection on their sensitivity.
      ). In the healthy quarters, the greater SCC responses in BME compared with AME samples was in contrast to previous findings where lower SCC was reported in the BME samples taken within 40 s of udder stimulation (
      • Bruckmaier R.M.
      • Ontsouka C.E.
      • Blum J.W.
      Fractionized milk composition in dairy cows with subclinical mastitis.
      ). The reason for such variation could be that in this study the BME samples were taken within ∼60 s of udder touch and there were 8 quarters with prior dipping before BME sampling, leading to reduced concentration of SCC in AME samples (
      • Lehmann M.
      • Wall S.K.
      • Wellnitz O.
      • Bruckmaier R.M.
      Changes in milk L-lactate, lactate dehydrogenase, serum albumin, and IgG during milk ejection and their association with somatic cell count.
      ). In mid-lactation period, the lower SCC and LDH (BME samples) responses in gram-negative IMI were due to a lower number of coliform cases (2 out of 12), as the majority of the coliform IMI (6 out of 12) were found in the early-lactation period (within 1 to 57 DIM). It should be noted that the number of gram-negative coliform observations was not very large (n = 12) in this study, though there was sufficient statistical power to obtain statistical significance when comparing with the gram-positive IMI. Further studies should consider a herd approach to identify diversified and larger numbers of specific pathogens by a bulk-tank milk PCR test and investigation of the unequal distribution of the blood-milk barrier density in the entire mammary gland by measuring markers in fractionized milk samples collected in a noninvasive way.

      CONCLUSIONS

      Both LDH and SA could be used as replacement markers to indicate IgG transfer via the blood-milk barrier, and both are suitable to differentiate gram-positive and gram-negative IMI when combined with information generated from SCC. Further, SCC, IgG, LDH, and SA responses differ in vivo between gram-positive and gram-negative IMI. Foremilk BME is more sensitive to changes of blood-milk barrier integrity caused by IMI because milk ejection affects SCC, IgG, LDH, and SA in milk. Therefore, in AMS it would be advantageous for sensitive detection of IMI to perform the respective measurements in foremilk BME. Earlier differentiation of IMI type based on BME will be expected to minimize the antibiotic used for clinical mastitis.

      ACKNOWLEDGMENTS

      The authors thank Katrina Bosward from the University of Sydney for assistance in Coxiella burnetii antibody detection. We also thank Yolande Zbinden (laboratory staff of the University of Bern, Switzerland) for excellent technical assistance in the laboratory activities. The authors thank Ian Chapman and John Garrod (farm staff of the University of Sydney) and Kim McKean and Oliver Roberts (former farm managers of the University of Sydney) for their assistance in collecting the samples and data. This project was supported by FutureDairy (Camden, NSW, Australia), the Dairy Research Foundation of The University of Sydney (Australia) and partly by the University of Bern (Switzerland). The first author is a recipient of an Australian Endeavour PhD Scholarship. The authors have not stated any conflicts of interest.

      REFERENCES

        • Bruckmaier R.M.
        • Ontsouka C.E.
        • Blum J.W.
        Fractionized milk composition in dairy cows with subclinical mastitis.
        Vet. Med. (Praha). 2004; 49: 283-290
        • Bruckmaier R.M.
        • Weiss D.
        • Wiedemann M.
        • Schmitz S.
        • Wendl G.
        Changes of physicochemical indicators during mastitis and the effects of milk ejection on their sensitivity.
        J. Dairy Res. 2004; 71 (15354578): 316-321
        • Bruckmaier R.M.
        • Wellnitz O.
        Triennial Lactation Symposium/BOLFA: Pathogen-specific immune response and changes in the blood–milk barrier of the bovine mammary gland.
        J. Anim. Sci. 2017; 95 (29293747): 5720-5728
        • Burton J.L.
        • Erskine R.J.
        Immunity and mastitis. Some new ideas for an old disease.
        Vet. Clin. North Am. Food Anim. Pract. 2003; 19 (12682934): 1-45
        • Butler D.G.
        • Cullis B.R.
        • Gilmour A.R.
        • Gogel B.G.
        • Thompson R.
        ASReml-R Reference Manual Version 4.
        VSN International Ltd, 2017
        • Cheng J.
        • Zhang J.
        • Yang J.
        • Yi B.
        • Liu G.
        • Zhou M.
        • Kastelic J.P.
        • Han B.
        • Gao J.
        Klebsiella pneumoniae infection causes mitochondrial damage and dysfunction in bovine mammary epithelial cells.
        Vet. Res. 2021; 52 (33568211): 17
        • Dzidic A.
        • Macuhova J.
        • Bruckmaier R.M.
        Effects of cleaning duration and water temperature on oxytocin release and milk removal in an automatic milking system.
        J. Dairy Sci. 2004; 87 (15545378): 4163-4169
        • Halasa T.
        • Huijps K.
        • Østerås O.
        • Hogeveen H.
        Economic effects of bovine mastitis and mastitis management: A review.
        Vet. Q. 2007; 29 (17471788): 18-31
        • Hanley J.A.
        • McNeil B.J.
        The meaning and use of the area under a receiver operating characteristic (ROC) curve.
        Radiology. 1982; 143 (7063747): 29-36
        • Hernández-Castellano L.E.
        • Wall S.K.
        • Stephan R.
        • Corti S.
        • Bruckmaier R.
        Milk somatic cell count, lactate dehydrogenase activity, and immunoglobulin G 2 concentration associated with mastitis caused by different pathogens: A field study.
        Schweiz. Arch. Tierheilkd. 2017; 159 (28475483): 283-290
        • Hogeveen H.
        • Kamphuis C.
        • Steeneveld W.
        • Mollenhorst H.
        Sensors and clinical mastitis-the quest for the perfect alert.
        Sensors (Basel). 2010; 10 (22163637): 7991-8009
        • Huijps K.
        • Lam T.J.
        • Hogeveen H.
        Costs of mastitis: facts and perception.
        J. Dairy Res. 2008; 75 (18226298): 113-120
        • Khatun M.
        • Bruckmaier R.M.
        • Thomson P.C.
        • House J.
        • García S.C.
        Suitability of somatic cell count, electrical conductivity, and lactate dehydrogenase activity in foremilk before versus after alveolar milk ejection for mastitis detection.
        J. Dairy Sci. 2019; 102 (31351709): 9200-9212
        • Lehmann M.
        • Wall S.K.
        • Wellnitz O.
        • Bruckmaier R.M.
        Changes in milk L-lactate, lactate dehydrogenase, serum albumin, and IgG during milk ejection and their association with somatic cell count.
        J. Dairy Res. 2015; 82 (25467384): 129-134
        • Lehmann M.
        • Wellnitz O.
        • Bruckmaier R.M.
        Concomitant lipopolysaccharide-induced transfer of blood-derived components including immunoglobulins into milk.
        J. Dairy Sci. 2013; 96 (23219120): 889-896
        • Mačuhová J.
        • Tančin V.
        • Bruckmaier R.M.
        Oxytocin release, milk ejection and milk removal in a multi-box automatic milking system.
        Livest. Prod. Sci. 2003; 81: 139-147
        • Milner P.
        • Page K.L.
        • Hillerton J.E.
        The effects of early antibiotic treatment following diagnosis of mastitis detected by a change in the electrical conductivity of milk.
        J. Dairy Sci. 1997; 80 (9178126): 859-863
        • Nickerson S.C.
        • Pankey J.W.
        Cytologic observations of the bovine teat end.
        Am. J. Vet. Res. 1983; 44 (6625294): 1433-1441
        • Pyörälä S.
        Indicators of inflammation in the diagnosis of mastitis.
        Vet. Res. 2003; 34 (14556695): 565-578
        • Robin X.
        • Turck N.
        • Hainard A.
        • Tiberti N.
        • Lisacek F.
        • Sanchez J.
        • Mueller M.
        pROC: An open-source package for R and S+ to analyze and compare ROC curves.
        BMC Bioinformatics. 2011; 12 (21414208): 77
        • Ruopp M.D.
        • Perkins N.J.
        • Whitcomb B.W.
        • Schisterman E.F.
        Youden Index and optimal cut-point estimated from observations affected by a lower limit of detection.
        Biom. J. 2008; 50 (18435502): 419-430
        • Sarikaya H.
        • Werner-Misof C.
        • Atzkern M.
        • Bruckmaier R.M.
        Distribution of leucocyte populations, and milk composition, in milk fractions of healthy quarters in dairy cows.
        J. Dairy Res. 2005; 72 (16223466): 486-492
        • Schmitz S.
        • Pfaffl M.W.
        • Meyer H.H.D.
        • Bruckmaier R.M.
        Short-term changes of mRNA expression of various inflammatory factors and milk proteins in mammary tissue during LPS-induced mastitis.
        Domest. Anim. Endocrinol. 2004; 26 (14757184): 111-126
        • Shamay A.
        • Homans R.
        • Fuerman Y.
        • Levin I.
        • Barash H.
        • Silanikove N.
        • Mabjeesh S.J.
        Expression of albumin in nonhepatic tissues and its synthesis by the bovine mammary gland.
        J. Dairy Sci. 2005; 88 (15653522): 569-576
        • Steeneveld W.
        • Vernooij J.C.M.
        • Hogeveen H.
        Effect of sensor systems for cow management on milk production, somatic cell count, and reproduction.
        J. Dairy Sci. 2015; 98 (25841965): 3896-3905
        • Stelwagen K.
        • Prosser C.G.
        • Davis S.R.
        • Farr V.C.
        • Politis I.
        • White J.H.
        • Zavizion B.
        Effect of milking frequency and somatotropin on the activity of plasminogen activator, plasminogen, and plasmin in bovine milk.
        J. Dairy Sci. 1994; 77 (7699135): 3577-3583
        • Swets J.A.
        Measuring the accuracy of diagnostic systems.
        Science. 1988; 240 (3287615): 1285-1293
        • Verbyla A.P.
        • Cullis B.R.
        • Kenward M.G.
        • Welham S.J.
        The analysis of designed experiments and longitudinal data by using smoothing splines.
        Appl. Stat. 1999; 48: 269-311
        • Wagner S.A.
        • Jones D.E.
        • Apley M.D.
        Effect of endotoxic mastitis on epithelial cell numbers in the milk of dairy cows.
        Am. J. Vet. Res. 2009; 70 (19496671): 796-799
        • Wall S.K.
        • Gross J.J.
        • Kessler E.C.
        • Villez K.
        • Bruckmaier R.M.
        Blood-derived proteins in milk at start of lactation: Indicators of active or passive transfer.
        J. Dairy Sci. 2015; 98 (26298756): 7748-7756
        • Wall S.K.
        • Wellnitz O.
        • Hernández-Castellano L.E.
        • Ahmadpour A.
        • Bruckmaier R.M.
        Supraphysiological oxytocin increases the transfer of immunoglobulins and other blood components to milk during lipopolysaccharide- and lipoteichoic acid–induced mastitis in dairy cows.
        J. Dairy Sci. 2016; 99 (27592421): 9165-9173
        • Wellnitz O.
        • Arnold E.T.
        • Bruckmaier R.M.
        Lipopolysaccharide and lipoteichoic acid induce different immune responses in the bovine mammary gland.
        J. Dairy Sci. 2011; 94 (22032363): 5405-5412
        • Wellnitz O.
        • Arnold E.T.
        • Lehmann M.
        • Bruckmaier R.M.
        Short communication: Differential immunoglobulin transfer during mastitis challenge by pathogen-specific components.
        J. Dairy Sci. 2013; 96 (23295119): 1681-1684
        • Wellnitz O.
        • Bruckmaier R.M.
        The innate immune response of the bovine mammary gland to bacterial infection.
        Vet. J. 2012; 192 (22498784): 148-152
        • Wellnitz O.
        • Zbinden C.
        • Huang X.
        • Bruckmaier R.M.
        Short communication: Differential loss of bovine mammary epithelial barrier integrity in response to lipopolysaccharide and lipoteichoic acid.
        J. Dairy Sci. 2016; 99 (27060811): 4851-4856