Advertisement
Research| Volume 101, ISSUE 3, P2423-2432, March 2018

Herd-level prevalence of the ultrasonographic lung lesions associated with bovine respiratory disease and related environmental risk factors

  • S. Buczinski
    Correspondence
    Corresponding author
    Affiliations
    Département des Sciences Cliniques, Faculté de Médecine Vétérinaire, Université de Montréal, Saint-Hyacinthe J2S 2M2, Québec, Canada
    Search for articles by this author
  • Author Footnotes
    2 Current address: Vankleek Hill Veterinary Services, Vankleek Hill, Ontario, Canada.
    M.E. Borris
    Footnotes
    2 Current address: Vankleek Hill Veterinary Services, Vankleek Hill, Ontario, Canada.
    Affiliations
    Département des Sciences Cliniques, Faculté de Médecine Vétérinaire, Université de Montréal, Saint-Hyacinthe J2S 2M2, Québec, Canada
    Search for articles by this author
  • J. Dubuc
    Affiliations
    Département des Sciences Cliniques, Faculté de Médecine Vétérinaire, Université de Montréal, Saint-Hyacinthe J2S 2M2, Québec, Canada
    Search for articles by this author
  • Author Footnotes
    2 Current address: Vankleek Hill Veterinary Services, Vankleek Hill, Ontario, Canada.
Open ArchivePublished:December 28, 2017DOI:https://doi.org/10.3168/jds.2017-13459

      ABSTRACT

      Bovine respiratory disease complex (BRD) is a major calf disease during the preweaning period. Thoracic ultrasound (TUS) has been recently described as a reliable tool for assessing BRD-associated lung lesions. The objectives of this study were to define the herd-level prevalence of lung consolidation assessed by TUS (CONSTUS). A total of 39 Québec dairy herds were randomly chosen to participate in this cross-sectional study. Between 6 and 12 preweaned calves were examined for signs of CONSTUS defined by any site with ≥3 cm consolidated lung tissue during 1 visit in summer and 1 visit in winter. Herd questionnaire focused on calf health and housing data [airborne bacteria (aerobic, coliform, yeasts and mold counts), air drafts, temperature, hygrometry, and ammonia levels] were also collected during these visits looking for potential association with CONSTUS prevalence. The median herd-level of CONSTUS prevalence (interquartile range) were 8% (0–22%) in summer and 15% (0–35%) in winter. Multivariable analyses showed that season was associated with an increased CONSTUS prevalence [10.0 vs. 19.3%; odds ratio (OR) = 2.16], as well as group housing during preweaned period (9.6 vs. 20.0%; OR = 2.37) and perceived BRD problem by the farmer (10.6 vs. 18.3%; OR = 1.89). Despite tremendous changes in calves' environment between winter and summer, none of the housing variables were associated with the CONSTUS prevalence in this study. Based on the observed CONSTUS prevalence ranges in the present study (herds in the lower prevalence quartile), an achievable goal of none of 12 calves with consolidation ≥3 cm can potentially be used to define a low-BRD risk herd.

      Key words

      INTRODUCTION

      Bovine respiratory disease complex (BRD) is a multifactorial disease resulting from interaction between the host, the microbiological agents, and environmental factors. In dairy calves, enzootic pneumonia is classically used to depict the typical pattern of occurrence of several cases that regularly occur (
      • Van Donkersgoed J.
      • Ribble C.S.
      • Boyer L.G.
      • Townsend H.G.
      Epidemiological study of enzootic pneumonia in dairy calves in Saskatchewan.
      ). Besides the specific effect related to morbidity and mortality costs, BRD also has a negative effect on growth and long-term productivity (
      • Guterbock W.M.
      The impact of BRD: the current dairy experience.
      ). Predisposing factors of BRD include inadequate transfer of passive immunity (
      • Van Donkersgoed J.
      • Ribble C.S.
      • Boyer L.G.
      • Townsend H.G.
      Epidemiological study of enzootic pneumonia in dairy calves in Saskatchewan.
      ;
      • Furman-Fratczak K.
      • Rzasa A.
      • Stefaniak T.
      The influence of colostral immunoglobulin concentration in heifer calves' serum on their health and growth.
      ;
      • Windeyer M.C.
      • Leslie K.E.
      • Godden S.M.
      • Hodgins D.C.
      • Lissemore K.D.
      • LeBlanc S.J.
      Factors associated with morbidity, mortality, and growth of dairy heifer calves up to 3 months of age.
      ), inadequate vaccination protocols (
      • Gorden P.J.
      • Plummer P.
      Control, management, and prevention of bovine respiratory disease in dairy calves and cows.
      ), as well as diet restriction (e.g., ingestion of low milk amounts) which can depress immune function (
      • Khan M.A.
      • Weary D.M.
      • von Keyserlingk M.A.
      Invited review: Effects of milk ration on solid feed intake, weaning, and performance in dairy heifers.
      ). The environmental factors commonly reported to be associated with greater risk of BRD include group housing (
      • Svensson C.
      • Liberg P.
      The effect of group size on health and growth rate of Swedish dairy calves housed in pens with automatic milk-feeders.
      ), poor bedding quality (
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ), large group size (
      • Svensson C.
      • Liberg P.
      The effect of group size on health and growth rate of Swedish dairy calves housed in pens with automatic milk-feeders.
      ), direct contact with older animals (
      • Lundborg G.K.
      • Svensson E.C.
      • Oltenacu P.A.
      Herd-level risk factors for infectious diseases in Swedish dairy calves aged 0–90 days.
      ), and poor microbiological air quality (
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ).
      The clinical signs of BRD can be insidious (
      • McGuirk S.M.
      • Peek S.F.
      Timely diagnosis of dairy calf respiratory disease using a standardized scoring system.
      ); although helpful, clinical signs are not perfectly accurate to diagnose BRD (
      • Ollivett T.L.
      • Caswell J.L.
      • Nydam D.V.
      • Duffield T.
      • Leslie K.E.
      • Hewson J.
      • Kelton D.
      Thoracic ultrasonography and bronchoalveolar lavage fluid analysis in Holstein calves with subclinical lung lesions.
      ). The clinical presentation of the disease may be difficult to diagnose correctly despite the existence of several screening strategies based on clinical scoring (
      • McGuirk S.M.
      Disease management of dairy calves and heifers. T.
      ;
      • Aly S.S.
      • Love W.J.
      • Williams D.R.
      • Lehenbauer T.W.
      • Van Eenennaam A.
      • Drake C.
      • Kass P.H.
      • Farver T.B.
      Agreement between bovine respiratory disease scoring systems for pre-weaned dairy calves.
      ). Thus, thoracic ultrasound (TUS) has been recently mentioned as an interesting diagnostic tool that could be implemented on farms for a rapid scanning of dairy calves' lungs (
      • Ollivett T.L.
      • Buczinski S.
      On-farm use of ultrasonography for bovine respiratory disease.
      ). Observation of pleural and lung anomalies can be used as a proxy for diagnosis of BRD with a higher sensitivity and specificity than clinical scoring (
      • Buczinski S.
      • Ollivett T.L.
      • Dendukuri N.
      Bayesian estimation of the accuracy of the calf respiratory scoring chart and ultrasonography for the diagnosis of bovine respiratory disease in pre-weaned dairy calves.
      ) and a good correlation with necropsy findings (
      • Ollivett T.L.
      • Caswell J.L.
      • Nydam D.V.
      • Duffield T.
      • Leslie K.E.
      • Hewson J.
      • Kelton D.
      Thoracic ultrasonography and bronchoalveolar lavage fluid analysis in Holstein calves with subclinical lung lesions.
      ). Moreover, extensive lung lesions assessed by TUS at weaning are associated with an increased risk of being prematurely culled before the first calving (
      • Adams E.A.
      • Buczinski S.
      Short communication: Ultrasonographic assessment of lung consolidation postweaning and survival to the first lactation in dairy heifers.
      ;
      • Teixeira A.G.
      • McArt J.A.
      • Bicalho R.C.
      Thoracic ultrasound assessment of lung consolidation at weaning in Holstein dairy heifers: Reproductive performance and survival.
      ). In a randomized controlled trial assessing the efficacy of intranasal vaccination in dairy calves, animals that were found with at least 1 lesion ≥3 cm of lung consolidation depth had 0.1 lbs/d (45.4 g/d) less of daily weight gain than nonconsolidated calves over a time period of 56 d (
      • Ollivett T.
      Understanding the diagnosis and risk factors for respiratory disease in dairy calves.
      ). It would therefore be of interest to gather more data on ultrasonographic lesions in preweaning dairy calves as a way to assess calf health and productivity.
      Unfortunately, data quantifying calf- and herd-level prevalence of lung consolidation assessed by TUS (CONSTUS) are scant. In a previous study performed in 13 Québec dairy farms (
      • Buczinski S.
      • Forté G.
      • Francoz D.
      • Bélanger A.M.
      Comparison of thoracic auscultation, clinical score, and ultrasonography as indicators of bovine respiratory disease in preweaned dairy calves.
      ), 56 of 106 calves (53%) had ultrasonographic evidence of lung consolidation (at least 1 cm of consolidated lung). At the herd level, the prevalence of consolidation in this previous study varied from 0 to 78% (median 50%; interquartile range 50–60%). However, this study was performed on a small number of herds obtained by convenience sampling and focused at enrolling herds with enzootic BRD problems; therefore, the findings cannot be extrapolated to the general population of dairy farms. For these reasons, it is difficult to know the prevalence of CONSTUS in preweaning dairy calves at the herd level. Knowing the latter information as well as potential herd-level risk factors associated with such prevalence could be valuable to assess respiratory health in preweaning calves of dairy farms.
      Poor air quality and ventilation have classically been described as a risk factor for BRD in dairy calves at the herd level (
      • Pritchard D.G.
      • Carpenter C.A.
      • Morzaria S.P.
      • Harkness J.W.
      • Richards M.S.
      • Brewer J.I.
      Effect of air filtration on respiratory disease in intensively housed veal calves.
      ;
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ). The assessment of air quality can be performed in multiple ways (
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ). Air drafts as well as ammonia or airborne microbiological flora have been associated with an increased risk of BRD at the herd level (
      • Lundborg G.K.
      • Svensson E.C.
      • Oltenacu P.A.
      Herd-level risk factors for infectious diseases in Swedish dairy calves aged 0–90 days.
      ;
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ). In a previous study using an air-sampler for aerobic and coliform culture, the authors reported that aerobic counts, but not coliform counts, were associated with BRD prevalence in a naturally ventilated barn during the winter (
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ). Mold and fungi are also commonly found in livestock buildings (
      • Blom J.Y.
      • Madsen E.
      • Krogh H.
      • Wolstrup J.
      Numbers of airborne bacteria and fungi in calf houses.
      ) and have been associated with allergic respiratory diseases in horses or human (
      • Beeler-Marfisi J.
      • Clark M.E.
      • Wen X.
      • Sears W.
      • Huber L.
      • Ackerley C.
      • Viel L.
      • Bienzle D.
      Experimental induction of recurrent airway obstruction with inhaled fungal spores, lipopolysaccharide, and silica microspheres in horses.
      ;
      • Pomés A.
      • Chapman M.D.
      • Wunschmann S.
      Indoor Allergens and Allergic Respiratory Disease.
      ). Recently, on-farm culture systems have been developed and adapted for a practical use to assess etiologic agents of cows' mastitis (
      • Cameron M.
      • Keefe G.P.
      • Roy J.P.
      • Dohoo I.R.
      • MacDonald K.A.
      • McKenna S.L.
      Evaluation of a 3M Petrifilm on-farm culture system for the detection of intramammary infection at the end of lactation.
      ). The Petrifilm culture system has also been validated for the assessment of airborne bacteria, mold, and yeasts and was recently described as a relatively inexpensive alternative to assess microbiological quality of farm environment (
      • Buczinski S.
      • Borris M.E.
      • Dubuc J.
      Assessment of dairy calves' microbiological environment using Petrifilm bacteriology plates.
      ).
      Our hypotheses were that the prevalence of CONSTUS would be higher in winter than in summer and that it would be associated with the same herd risk factors as those associated with clinical BRD. The main objective of our study was to report the herd-level prevalence of CONSTUS in dairy calves from a large number of herds. A secondary objective was to describe environmental or herd factors associated with the proportion of consolidated calves on farms.

      MATERIALS AND METHODS

      A cross-sectional study was performed at the bovine ambulatory clinic of the Faculté de médecine vétérinaire of the Université de Montréal (St-Hyacinthe, QC, Canada) during 2 time periods, July to August 2015 (summer) and December 2015 to January 2016 (winter). The herd was the unit of interest in this study. Herd selection was randomly performed within all the herds using veterinary services from the bovine ambulatory clinic (n = 161; source population). For this purpose, a list of all the herds from the clinic was built. For practical considerations, all herds with limited size (<50 lactating cows) were excluded (not having enough calves to examine within a 2-mo period). The random function of Microsoft Excel (Microsoft Corp., Redmond, WA) was used to select herd from that list. Selected dairy producers were then asked to confirm their participation through an individual phone call.

      Sample Size Calculation

      The herd-level sample size estimation (i.e., herd number needed) for the main objective (prevalence) was estimated to be 32 herds and was based on the following assumptions: prevalence expected = 10%, α = 5% (for a 95% CI), and precision of estimate = 7.5%. This 10% prevalence was based on our best guess and a study in male dairy calves (n = 209) with low BRD risk systematically screened with ultrasound at 2 to 3 wk of age, where prevalence of lung lesion was 12% (
      • Berman J.
      • Francoz D.
      • Dubuc J.
      • Buczinski S.
      A randomised clinical trial of a metaphylactic treatment with tildipirosin for bovine respiratory disease in veal calves.
      ). The herd-level sample size estimation for the secondary objective (risk factors) was estimated to be 40 herds and was based on finding a difference of 35% in prevalence (absence of risk factor = 5%; presence of risk factor = 40%), and accounting for α and β errors of 5 and 20%, respectively. Therefore, recruitment of 40 herds was targeted. We were able to recruit a total of 39 herds for the study, which represented a 24% participation rate from the source population. In each participating herd, 6 to 12 preweaning calves were targeted to be sampled during a farm visit. This within-herd sample size calculation was based the following assumptions: prevalence expected = 5%, α = 25% (for a 75% CI), and precision of estimate = 10%.

      Assessment of Herds and Calf Rearing Characteristics

      During each study period, every herd was visited once by a veterinarian and a research technician for data collection. Specific information concerning calf management during the preweaning period was collected from the producer during the farm visit. The average numbers of milking cows during the sampling periods were obtained, as well as vaccination protocols used for the adults cows and for the preweaning calves. The total number of preweaning calves at the day of the visit was also recorded. If more than 12 preweaning calves were present at the time of the visit, the sampling among available calves was performed taking into account the weighted proportion of calves housed in similar pen or housing conditions to be representative of herd structure at this time. Specific information collected concerning the preweaning period were the nature of the liquid diet [milk replacer (MR) vs. waste milk], the maximal daily volume received by a calf during this period, as well as the criteria for weaning (age, weight, or both). The feeding technique (individual feeding such as bottle, bucket with or without nipple vs. automatic milk feeder) and the presence of free access to water for all calves were noted. The main bedding component and the nesting score (1 to 3 scale) in up to 5 calves were recorded as previously described (
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ). The calves that were lying down during the visit were scored and an average score was obtained by dividing the sum of the scores by the number of calves assessed. The type of contact with other animals was also quantified. Contact with other animals was considered as absent or minor if calves were raised individually and had a contact limited to their penmates. The contact was considered significant if animals were group-raised or had direct contact with other animals. The age of the animals that were in contact with preweaning calves was also taken into account. The farms were classified as having preweaning calves with no direct contact with other animals or preweaning calves with direct contact with animals of their age or direct or indirect contact with other animals (postweaning calves, heifers, or cows). The owners were also asked at the day of the visit if they had any concern on current pneumonia problems in preweaning calves.

      Assessment of Microbiological Air Quality

      The assessment of airborne bacteria, yeasts, and mold counts was performed using on farm culturing systems (Petrifilms, 3M Canada, London, ON, Canada). Three different types of plates were used. The total aerobic count plate (TAC), coliform count (CC), and yeasts and mold plates (YM) were used to assess the microbiological air quality in the barn, as previously reported (
      • Buczinski S.
      • Borris M.E.
      • Dubuc J.
      Assessment of dairy calves' microbiological environment using Petrifilm bacteriology plates.
      ). The plates were used according to the manufacturer specification for air sampling, incubation, and reading. The plates were rehydrated prior farm visit with 1 mL of sterile water solution and kept refrigerated up to 7 d. As the study objective was to focus on calves' environment, Petrifilm samples were only taken at the calf housing site. The site of air quality assessment needed to be representative of typical calf environment on that particular farm. In the case of external hutches, only 1 random hutch was sampled. In the case of calves raised attached in the front or back of the cows, because we found no difference between calf environment and feeding alley, only 1 sample was taken at the calf site. In the case of individual pens, 1 occupied pen was randomly sampled. When calves were raised in groups with automatic milk feeders, sampling was performed in the pen approximatively 2 m far from the milk line. The samples at the calf site were always taken in occupied pens, boxes, or hutches. To avoid any movement of the calf that could contaminate Petrifilm plates, the calves were attached and moved away from the individual pen or hutch during the sampling period. One operator was present in the pen of group-raised calves to limit risk of contamination due to calf movements. The samples were always taken at 20 to 30 cm from the ground to be closely associated with what a calf breathes when being down. The system used was a roof plastic base (Oatey, Cleveland, OH), where the plates were kept in place using metal binders avoiding any contact with growth medium. The Petrifilm manufacturer recommends air contact up to 15 min to assess airborne bacteria in factories. However, based on a previous study, we used a 10-min sampling period, which was considered as a good compromise contact time (
      • Buczinski S.
      • Borris M.E.
      • Dubuc J.
      Assessment of dairy calves' microbiological environment using Petrifilm bacteriology plates.
      ). Contact time assessment was assessed using chronometer function included in a smartphone (iPhone 5, Apple, Cupertino, CA).
      The Petrifilm plates were incubated (Boeker digital incubator, Boeker Scientific, Feasterville, PA) and read according to the manufacturer recommendations (https://www.3mcanada.ca/3m/en_ca/food-safety-testing-ca/products?Ntt=petrifilm&LC=en_CA&co=cc&gsaAction=scBR&type=cc). Briefly, the aerobic count was read 48 h after incubation and the coliform plate 24 h after incubation at 35°C. The yeast and mold plate were read 3 to 5 d after incubation at room temperature (20–25°C). In case the plate count exceeds 1,000 (i.e., too numerous to be counted), the value of 1,001 was given.

      Assessment of Air Draft, Temperature, Hygrometry, and Ammonia

      The sites that were used for microbiological sampling were also used for air draft assessment with a portative anemometer (Velocicalc 9565, TSI Inc., Shoreview, MN) assessing wind speed (0–50 m/s; accuracy ±0.015 m/s), relative hygrometry (0–95% relative humidity, accuracy ±3%), and temperature (−10 to 60°C, accuracy ±0.3°C). The sampling period was 2 min in duration, with the probe oriented to the ventilation system at the 20 to 30 cm from the ground. The recordings were stored in a laptop and minimal, maximal, and mean values were kept for further analysis. The ammonia level was also recorded at the same sites using a hand-held ammonia analyzer (GasAlertNH3, Honeywell Analytics, Lincolnshire, IL) with a range detection of 0 to 50 ppm. An increased ammonia level was considered if ammonia concentration was ≥5 ppm to take into account recent recommendations (between 3.5 and 7 ppm) in Canadian dairy farms (
      • Kaufman J.
      • Linington M.
      • Osborne V.
      • Wagner-Riddle C.
      • Wright T.
      Short communication: Field study of air ammonia concentrations in Ontario dairy calf housing microenvironments.
      ).

      Calf Data Measurements

      During the farm visits, each calf was examined and a systematic thoracic ultrasonography was performed for screening 3 different lung areas: the right and left lung area caudal to the heart, and the right cranial part of the lung (
      • Ollivett T.L.
      • Buczinski S.
      On-farm use of ultrasonography for bovine respiratory disease.
      ). Alcohol was sprayed on the intercostal space and the area of interest was screened using a linear probe (Draminski iScan, Olsztyn, Poland), focusing on any area of consolidation in particular. The maximal depth of consolidation (cm) was noted for each calf at each site, as described elsewhere, using the 1-cm2 grid on the screen (
      • Buczinski S.
      • Forté G.
      • Francoz D.
      • Bélanger A.M.
      Comparison of thoracic auscultation, clinical score, and ultrasonography as indicators of bovine respiratory disease in preweaned dairy calves.
      ). A calf was considered to have ultrasonographic evidence of consolidation if consolidation depth was ≥3 cm, as that depth has been previously associated with decreased preweaning ADG (
      • Ollivett T.
      Understanding the diagnosis and risk factors for respiratory disease in dairy calves.
      ).

      Statistical Analyses

      The statistical analyses were performed using commercial software such as SAS (version 9.4, SAS Institute Inc., Cary, NC) and MedCalc (version 13.1.0, Ostend, Belgium). Descriptive statistics concerning herd characteristics and prevalence of ultrasonographic consolidation were calculated using a depth of 3 cm as cutoff (<3 vs. ≥3 cm).

      Variable Selection

      The correlation between the different collected variables was performed using Spearman correlation (CORR procedure of SAS). Multicollinearity was considered for a Spearman correlation ≥0.7, and only 1 variable was kept for the analysis. For continuous data (bacterial and yeasts and mold counts), normality was assessed using Shapiro Wilk test in UNIVARIATE procedure of SAS. As the data were right-skewed, a natural logarithm transformation was made to improve data normality for TAC, CC, and YM counts before using them for univariable analysis. The association between the proportion of calves with lung lesions and calf age in weeks (from wk 1 to ≥9) was done using multiple comparisons test for proportion with Tukey correction for multiple comparisons (%COMPPROP macro of SAS). The median age of sampled calves in every herd was added as a potential covariable that could affect the prevalence of CONSTUS.

      Herd-Level Models

      The season was considered as a fixed effect and forced in the models. The within-herd prevalence of consolidated calves was assessed using a mixed logistic regression model using the GLIMMIX procedure of SAS (event/trial syntax, logit link). The farm was considered as the cluster of replication during both seasons [1]. The general linear mixed model was described as (
      • Stroup W.W.
      Generalized Linear Mixed Models: Modern Concepts, Methods and Applications.
      )
      Losgitπ=βX+Fb,
      [1]


      where π is the prevalence of calves with consolidation (calves consolidated/calves screened) in the ith farm; β is the vector of model parameters; X is a matrix of farm covariates that also include season as a fixed effect; the term b (0; σF2) is the random effect for farm vector F.
      Univariable analyses were performed (all the different tested variables are summarized in Table 1, Table 2). Categorical variables were offered specifying the CLASS statement and continuous transformed variables (log-transformed colony-forming unit plates and median age of calves per herd) were assessed using the same univariable mixed model taking into account farm clustering (GLIMMIX procedure). The variables with P < 0.2 using an F-test were offered to the multivariable model using a backward elimination strategy until all remaining variable had P-values <0.05. A 2-way interaction term (season × Petrifilm microbiological count) was also tested, as it was anticipated that season could potentially interact with microbiologic air quality. A question was asked to the farmers about their perception of having a BRD problem at the time of the farm visits. This information was used as a potential variable that could be positively associated with the proportion of consolidated calves. We also built a model without including this covariate.
      Table 1Descriptive statistics of environment characteristic in farms in winter and summer and univariable association with the prevalence of ultrasonographic consolidation (≥3 cm depth) in preweaned calves
      ParameterSummerWinterP-value (between season)
      P-values for comparing winter and summer environmental characteristics using a nonparametric Wilcoxon rank sum test for paired sample (winter vs. summer).
      P-value (univar)
      P-value resulting for univariable analysis. The counts were log-transformed and their association was tested in relation with prevalence of ultrasonographic lung consolidation (≥3 cm depth) using the event/trial syntax of the GLIMMIX procedure of SAS (version 9.4, SAS Institute Inc., Cary, NC).
      MeanMedianRangeMeanMedianRange
      Aerobic count (cfu/plate)43836590–10012501984–720<0.010.80
      Coliform count (cfu/plate)741–35210–18<0.010.11
      Yeasts count (cfu/plate)23120–30033171–3100.190.37
      Mold count (cfu/plate)60601–1501360–84<0.010.03
      Ammonia (ppm)0.400–51.3700–70.02NI
      NI (not included in the univariable analyses): these variables were not included in univariable analysis because of low variability observed (ammonia) or high correlation with season (temperature, hygrometry) ≥0.7.
      Mean air speed (m/s)0.210.100.03–1.570.090.080.01–0.30<0.01NI
      Maximal air speed (m/s)0.430.300.13–1.780.200.190.05–0.38<0.010.08
      Mean barn temperature (°C)22.822.818.2–28.513.613.57.8–18.0<0.01NI
      Maximal barn temperature (°C)23.223.218.5–29.513.9148–18.2<0.01NI
      Mean relative hygrometry (%)59.659.347.9–71.626.926.518.7–40.1<0.01NI
      Maximal relative hygrometry (%)61.260.650.1–74.027.326.719.0–40.5<0.01NI
      1 P-values for comparing winter and summer environmental characteristics using a nonparametric Wilcoxon rank sum test for paired sample (winter vs. summer).
      2 P-value resulting for univariable analysis. The counts were log-transformed and their association was tested in relation with prevalence of ultrasonographic lung consolidation (≥3 cm depth) using the event/trial syntax of the GLIMMIX procedure of SAS (version 9.4, SAS Institute Inc., Cary, NC).
      3 NI (not included in the univariable analyses): these variables were not included in univariable analysis because of low variability observed (ammonia) or high correlation with season (temperature, hygrometry) ≥0.7.
      Table 2Univariable association of various herd-level characteristics and the proportion of ultrasonographic consolidation (defined as ≥3 cm depth) in preweaned calves
      ItemProportion (%) of calves with consolidation per herd (LSM ± SEM
      LSM of fixed effect using the GLIMMIX procedure of SAS (SAS version 9.4, SAS Institute Inc., Cary, NC).
      )
      P-value
      P-values obtained with univariable analysis and type III tests of fixed effect considering herd random effect using the GLIMMIX procedure of SAS.
      Season
       Summer8.9 ± 2.1<0.01
       Winter17.4 ± 3.5
      Feeding
       No automated milk feeder10.3 ± 2.80.73
       Automated milk feeder12.4 ± 5.5
       Maximal milk allowance <8 L11.51 ± 4.90.86
       Maximal milk allowance ≥8 L10.5 ± 2.9
       Milk replacer10.9 ± 2.80.93
       Whole milk ± milk replacer10.3 ± 5.5
       All calves have water access12.7 ± 4.60.57
       No water access or limited water access9.9 ± 2.8
      Housing and management
       Individual housing until weaning6.7 ± 2.1<0.01
       Group housing before weaning19.5 ± 4.8
       No direct or indirect contact with older animals9.2 ± 3.00.21
       Direct or indirect contact with heifers20.8 ± 7.7
       Direct or indirect contact with adults8.9 ± 3.5
       Vaccination of preweaned calves for BRD
      Bovine respiratory disease complex.
      16.2 ± 4.30.25
       No vaccination of preweaned calves for BRD10.8 ± 2.7
       Air draft <0.3 m/s during the visit14.5 ± 3.00.08
       Air draft ≥0.3 m/s during the visit8.3 ± 2.8
       No perceived BRD problem during preweaning11.2 ± 2.20.05
       Perceived BRD problem during preweaning19.4 ± 4.3
      1 LSM of fixed effect using the GLIMMIX procedure of SAS (SAS version 9.4, SAS Institute Inc., Cary, NC).
      2 P-values obtained with univariable analysis and type III tests of fixed effect considering herd random effect using the GLIMMIX procedure of SAS.
      3 Bovine respiratory disease complex.
      Fit statistics for conditional distribution of the generalized linear mixed model was considered adequate if no sign of overdispersion was observed, so long as the conditional distribution of Pearson chi-squared or degrees of freedom was close to 1 (
      • Stroup W.W.
      Generalized Linear Mixed Models: Modern Concepts, Methods and Applications.
      ). The least squares means were calculated for all significant variables in the final models and tested with Tukey-Kramer test.

      RESULTS

      A total of 39 farms participated in the study, representing a 24.2% participation rate from the source population of herds (n = 161). These farms had an average of 104 milking cows (median = 80 cows; range = 55–280 cows). Thirty-eight herds (97.4%) used calf age as the sole criterion for weaning. Twenty herds (51.3%) used group housing during the preweaning period and 8 (20.5%) used an automated milk feeder. The maximal volume of milk or milk replacer that was allowed during preweaning varied greatly from a farm to another (median = 8 L, range = 4–16 L). The majority of farms (n = 32; 82.1%) used milk replacer, whereas some used raw milk or wasted milk for calf feeding (n = 7; 17.9%). Water access was available for all calves in 25 farms (64%). In 14 farms (36%), all calves had no direct contact with heifers or cows, whereas in 8 farms (21%) calves had a contact with heifers; in the 11 remaining farms a contact with adult was also possible. The percentage of farms where the owner perceived BRD as a problem were 23 and 26% for summer and winter, respectively. Eleven of 39 farms (28.2%) used vaccination for BRD pathogens in preweaning calves, whereas 38 herds used cow vaccination against respiratory viruses (including bovine respiratory syncytial virus, bovine herpesvirus type 1, bovine viral diarrhea virus, and parainfluenza 3 virus).
      A total of 608 calves were enrolled during the entire study (2 sampling periods for each herd): 319 calves (median = 8 per herd) during summer and 289 calves (median = 7 per herd) during winter samplings. Only 9 herds (23%) had 12 or more preweaning calves available at the time of the visit (12 to 27 available preweaning calves). In the remaining herds, all present preweaning calves were therefore included. The 50th percentile of the median age of sampled calves per herd was 38 d in summer [interquartile range (IQR) = 28–50 d] and 44 d in winter (IQR = 33–50 d). The median proportion of ultrasonographic consolidation was 8% (IQR = 0–22%) in summer and 15% (IQR = 0–35%) in winter. The seasonal distribution is presented in Figure 1.
      Figure thumbnail gr1
      Figure 1Seasonal distribution of the prevalence of preweaned calves with lung consolidation at the herd level. The distribution of proportion of consolidated calves (consolidation defined as at least 1 site of consolidation ≥3 cm) in 39 dairy herds during summer (black) and winter (white) season is represented. The prevalence of consolidation was higher in winter than in summer (P < 0.01; see details in ).
      The Petrifilm plate counts [median (ranges) in cfu/plate] were 365 (90–1,001) for TAC, 4 (1–35) for CC; and 12 (0–300) for yeasts and 60 (1–150) for molds in summer. During the winter, the TAC were 197.5 (4–720), CC 1 (0–18), and yeasts 16.5 (1–310) and mold 6 (0–84) and are summarized in Table 1. We found a seasonal difference for TAC (P < 0.01), CC (P < 0.01), and molds (P < 0.01) but not for yeast counts (P = 0.19).
      The univariable statistical analysis results are summarized in Table 2, Table 3. Despite the fact that the probability of consolidation was significantly lower for calves in their first week of life when compared with calves in their fifth or more week (Figure 2), the age of calves in the herd was not associated with the odds of consolidation when controlling for herd effect (P = 0.37). Two multivariable models were built including the farmer's perception of whether the herd had a problem with BRD or not (Table 3). In both models, the season and group housing during preweaning were associated with an increased odd of consolidation.
      Table 3Effect of herd-level risk factors on the prevalence of ultrasonographic consolidation (≥3 cm depth) in preweaned calves
      ItemLevelβ
      Regression coefficient parameter.
      SE βOR
      Odds ratio.
      95% CI (OR)LSM ± SEMP-value
      Model with BRD
      Bovine respiratory disease complex.
      Intercept2.180.34<0.01
       SeasonSummerReferent10.0 ± 2.3<0.01
      Winter0.770.242.161.32–3.5319.3 ± 3.7
       Perceived BRD problemNoReferent10.6 ± 2.20.05
      Yes0.640.311.891.00–3.5518.3 ± 4.5
       Group housing during preweaningNoReferent9.6 ± 2.80.04
      Yes0.860.392.371.06–5.2920.0 ± 4.2
      Model excluding BRD perceptionIntercept−2.830.36<0.01
       SeasonSummer08.6 ± 2.1<0.01
      Winter0.760.242.151.32–3.5016.9 ± 3.3
       Group housing during preweaningNo08.0 ± 2.40.03
      Yes0.940.412.551.10–5.9018.1 ± 4.0
      1 Regression coefficient parameter.
      2 Odds ratio.
      3 Bovine respiratory disease complex.
      Figure thumbnail gr2
      Figure 2Distribution of the proportion of calves with lung consolidation based on age stratification. The number of calves per week of age as well as the proportion (±SD) of calves with lung lesion (defined with a lesion of consolidation ≥3 cm) are presented. The significant differences (P < 0.05) between weeks adjusted for multiple comparisons are presented using brackets.

      DISCUSSION

      Thoracic ultrasound has recently gained in interest due to its potential for helping in BRD diagnosis (
      • Ollivett T.L.
      • Buczinski S.
      On-farm use of ultrasonography for bovine respiratory disease.
      ). Our study found that the prevalence of lung lesions using thoracic ultrasound is relatively common in Québec dairy herds and is increased in herds with a perceived problem of BRD, herds that are group housed (vs. individually housed), and during winter season (vs. summer). These herd-level risk factors have been previously associated with increased odds of respiratory disease or increased respiratory sounds when assessed by a veterinarian in a Swedish study (
      • Lundborg G.K.
      • Svensson E.C.
      • Oltenacu P.A.
      Herd-level risk factors for infectious diseases in Swedish dairy calves aged 0–90 days.
      ). Unfortunately, and because of the wide variation of calf-raising systems in these relatively small farms, it was not possible to consistently assess the stocking density as a potential confounding factor for group-housed calves. Previous studies have shown that the size of the group rather than stocking density was associated with increased odds of respiratory disease at the calf level (
      • Svensson C.
      • Lundborg K.
      • Emanuelson U.
      • Olsson S.O.
      Morbidity in Swedish dairy calves from birth to 90 days of age and individual calf-level risk factors for infectious diseases.
      ).
      Not surprisingly, our results showed that having a perceived BRD problem in preweaning calves was associated with an increased herd prevalence of ultrasonographic pulmonary consolidation. Producers' perception of BRD problems has not been extensively studied; however, perception of infectious disease may be affected by various factors (
      • Frössling J.
      • Nöremark M.
      Differing perceptions—Swedish farmers' views of infectious disease control.
      ), including the accuracy of case detection. The diagnosis of enzootic pneumonia is a challenge in the absence of a reliable and accurate test. Application of clinical scoring charts, although appealing, may be difficult to implement in small farms on a regular basis (
      • McGuirk S.M.
      Disease management of dairy calves and heifers. T.
      ;
      • Love W.J.
      • Lehenbauer T.W.
      • Kass P.H.
      • Van Eenennaam A.L.
      • Aly S.S.
      Development of a novel clinical scoring system for on-farm diagnosis of bovine respiratory disease in pre-weaned dairy calves.
      ) due to the limited number of working people. Moreover, the sensitivity of producer detection or clinical scoring systems to detect BRD is generally lower than the specificity (
      • Sivula N.
      • Ames T.
      • Marsh W.
      • Werdin R.
      Descriptive epidemiology of morbidity and mortality in Minnesota dairy heifer calves.
      ;
      • Love W.J.
      • Lehenbauer T.W.
      • Van Eenennaam A.L.
      • Drake C.M.
      • Kass P.H.
      • Farver T.B.
      • Aly S.S.
      Sensitivity and specificity of on-farm scoring systems and nasal culture to detect bovine respiratory disease complex in preweaned dairy calves.
      ), which can explain why perception of a problem would be more related to the presence of a true problem (higher positive predictive value) than the opposite (i.e., not having a perceived BRD problem would not rule out a true BRD problem).
      One of the limitations of thoracic ultrasound is that it is an indicator of cumulative incidence of BRD episodes for a particular calf. Indeed, when a calf faced a BRD episode, depending on the etiology (viral, bacterial, or mixed), lung lesions may appear as soon as a few hours, such as after an experimental challenge with Mannheimia haemolytica (
      • Ollivett T.
      • Hewson J.
      • Schubotz R.
      • Caswell J.
      Ultrasonographic progression of lung consolidation after experimental infection with Mannheimia haemolytica in Holstein calves.
      ). However, little is known on the persistence of lung consolidation after BRD healing of an initial lung injury. For this reason, the proportion of calves with ultrasonographic lung lesions represent a mix of acute and chronic lung lesions. No specific ultrasonographic findings have been reported to differentiate between acute and chronic lesions at this time.
      At the calf level, it has been previously shown that extensive lung lesions at weaning have been associated with increased hazards of being culled before the first calving (
      • Adams E.A.
      • Buczinski S.
      Short communication: Ultrasonographic assessment of lung consolidation postweaning and survival to the first lactation in dairy heifers.
      ;
      • Teixeira A.G.
      • McArt J.A.
      • Bicalho R.C.
      Thoracic ultrasound assessment of lung consolidation at weaning in Holstein dairy heifers: Reproductive performance and survival.
      ) and decreased reproduction efficiency after first calving (
      • Teixeira A.G.
      • McArt J.A.
      • Bicalho R.C.
      Thoracic ultrasound assessment of lung consolidation at weaning in Holstein dairy heifers: Reproductive performance and survival.
      ). For these reasons, it is also of interest to accurately monitor lung lesions as a potential indicator of replacement heifer productivity. However, because decisions to change specific heifer management practices are made at the farm level, it was the objective of our study to be able to use thoracic ultrasound as a farm-level tool rather than as an individual test. An acceptable herd-level prevalence of consolidation threshold still needs to be established, but, based on the present study results, the first quartile of herds (i.e., 0% of calves with consolidation ≥3 cm) would be a practical recommendation to define herds without pneumonia problems. In contrast, the last quartiles of herds (with consolidation proportion >22% in summer and 35% in winter) could be a tentative reasonable guess for defining herds with abnormal proportions of consolidation. One of the limitations of our study was inherent to Québec dairy herd characteristics; as the herds are relatively small, it may be challenging to obtain a high number of calves to obtain a precise estimate of consolidation proportion in a limited period. From our study, a rule of thumb would be to have at least 6 calves with negative ultrasonographic examination to rule out a problem in winter [with a false-negative risk of 7.5% (1–0.35), with 0.35 defining the lower bound of the last quartile of winter prevalence, therefore assuming that every calf has the same probability of being consolidated] or 9 negative calves in summer [false-negative risk of 10% (1–0.22), with 0.22 defining the lower bound of the last quartile of summer prevalence]. Interestingly, despite the fact that very young calves (1st and 2nd week of life) had lower probability of presenting lung lesions, the median age of sampled calves in the herd was not associated with the herd prevalence of consolidation when accounting for herd clustering; this suggests that herd-specific risk factors were more important than calf age distribution. Even if the median age of calves was not associated with the prevalence of consolidation in the present study, the probability of consolidation was low in calves less than 2 wk old, indicating that choosing calves ≥3 wk old would increase the chances of finding lung consolidation. Care should therefore be taken extrapolating our findings, where the majority of calves are <2 wk old at time of ultrasound examination.
      It is also important to keep in mind that finding calves with lung lesions does not mean that lung consolidation should be used as a tool for treatment purpose. Studies in horse farms with enzootic Rhodococcus equi pulmonary infections have shown that relying only on lung lesions to treat foals increased antimicrobial use with no beneficial effect on foal health (
      • Venner M.
      • Astheimer K.
      • Lammer M.
      • Giguere S.
      Efficacy of mass antimicrobial treatment of foals with subclinical pulmonary abscesses associated with Rhodococcus equi.
      ). Such information is currently missing in cattle.
      We noted a seasonal effect on the prevalence of ultrasonographic lung consolidation, but we did not find any association between the calves' specific environment parameters and the prevalence of consolidation. Mold and coliform counts, as well as maximal air speed, were the only environmental variables that were selected during the univariable analysis, but they were partially confounded with season effect. As the environment assessment was done only once, and the observed consolidation could be related to a BRD episode that occurred days or weeks before the examination, the punctual measurement of environmental parameters (air speed, aerobic, coliform, YM counts, ammonia) may not reflect accurately the environmental conditions that the calves had before the day of data collection. Moreover, the limited numbers of herds may also be a limitation of the present study. Interestingly, ammonia levels were consistently low throughout our study and even lower than previously reported in winter (
      • Seedorf J.
      • Hartung J.
      Survey of ammonia concentrations in livestock buildings.
      ;
      • Lago A.
      • McGuirk S.M.
      • Bennett T.B.
      • Cook N.B.
      • Nordlund K.V.
      Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
      ) or in summer (
      • Kaufman J.
      • Linington M.
      • Osborne V.
      • Wagner-Riddle C.
      • Wright T.
      Short communication: Field study of air ammonia concentrations in Ontario dairy calf housing microenvironments.
      ). The inflammatory effect of high levels of ammonia on the respiratory tract has been mentioned in calves (
      • Dewes H.F.
      • Goodall G.
      Some preliminary observations on the possible relationship between ammonia production from soiled bedding in calf rearing sheds and calf illness.
      ), but in a more recent Swedish study at the herd-level ammonia levels <6 ppm were counterintuitively associated with an increased odd of respiratory disease (odds ratio = 2.46; 95% CI = 1.16–4.93;
      • Lundborg G.K.
      • Svensson E.C.
      • Oltenacu P.A.
      Herd-level risk factors for infectious diseases in Swedish dairy calves aged 0–90 days.
      ). The association between this air quality parameter and BRD deserves future research.
      In conclusion, our study showed that lung consolidation lesions detected by thoracic ultrasonography were present (median herd prevalence) in 8 and 15% of preweaning calves at the herd level in summer and winter, respectively. The prevalence of lung lesions was higher during winter and in farms that used group-housing during preweaning. Even if no specific herd-level prevalence threshold can be definitely proven from the present study, thoracic ultrasonography in preweaning dairy calves has potential to be used as a monitoring tool for respiratory health assessment in calves. From the present study, having no evidence of lung consolidation (≥3 cm lung consolidation) in up to 12 preweaning calves should be a goal for defining low-BRD risk farms. As BRD risks and herd dynamics or management are subject of change, assessment of ultrasonographic lung lesions should be done periodically at the herd level.

      Acknowledgments

      The authors thank the Zoetis clinical research fund (bovine ambulatory clinic) of the Université de Montréal (St-Hyacinthe, QC, Canada) for partially financing the study as well as 3M-Diagnostics (London, ON, Canada) for providing some of the Petrifilm used in that study for environment assessment. Parts of the material used in this study (anemometer, ammoniac detector) were obtained from a grant from the Fondation du Centre Hospitalier Universitaire Vétérinaire of the Université de Montréal (St-Hyacinthe, Québec, Canada). The authors also want to thank Clotilde Cros and Anne-Laure Lescarret (École Nationale Vétérinaire de Toulouse, France) for their valuable help in data collection during the summer part of the study and Jean-Philippe Pelletier (Faculté de Médecine Vétérinaire, St-Hyacinthe, Québec, Canada) for his valuable help during the winter sampling period.

      REFERENCES

        • Adams E.A.
        • Buczinski S.
        Short communication: Ultrasonographic assessment of lung consolidation postweaning and survival to the first lactation in dairy heifers.
        J. Dairy Sci. 2016; 99 (26686717): 1465-1470
        • Aly S.S.
        • Love W.J.
        • Williams D.R.
        • Lehenbauer T.W.
        • Van Eenennaam A.
        • Drake C.
        • Kass P.H.
        • Farver T.B.
        Agreement between bovine respiratory disease scoring systems for pre-weaned dairy calves.
        Anim. Health Res. Rev. 2014; 15 (25424381): 148-150
        • Beeler-Marfisi J.
        • Clark M.E.
        • Wen X.
        • Sears W.
        • Huber L.
        • Ackerley C.
        • Viel L.
        • Bienzle D.
        Experimental induction of recurrent airway obstruction with inhaled fungal spores, lipopolysaccharide, and silica microspheres in horses.
        Am. J. Vet. Res. 2010; 71 (20513185): 682-689
        • Berman J.
        • Francoz D.
        • Dubuc J.
        • Buczinski S.
        A randomised clinical trial of a metaphylactic treatment with tildipirosin for bovine respiratory disease in veal calves.
        BMC Vet. Res. 2017; 13 (28615070): 176
        • Blom J.Y.
        • Madsen E.
        • Krogh H.
        • Wolstrup J.
        Numbers of airborne bacteria and fungi in calf houses.
        Nord. Vet. Med. 1984; 36 (6493973): 215-220
        • Buczinski S.
        • Borris M.E.
        • Dubuc J.
        Assessment of dairy calves' microbiological environment using Petrifilm bacteriology plates.
        Bov. Pract. 2017; 51: 54-58
        • Buczinski S.
        • Forté G.
        • Francoz D.
        • Bélanger A.M.
        Comparison of thoracic auscultation, clinical score, and ultrasonography as indicators of bovine respiratory disease in preweaned dairy calves.
        J. Vet. Intern. Med. 2014; 28 (24236441): 234-242
        • Buczinski S.
        • Ollivett T.L.
        • Dendukuri N.
        Bayesian estimation of the accuracy of the calf respiratory scoring chart and ultrasonography for the diagnosis of bovine respiratory disease in pre-weaned dairy calves.
        Prev. Vet. Med. 2015; 119 (25794838): 227-231
        • Cameron M.
        • Keefe G.P.
        • Roy J.P.
        • Dohoo I.R.
        • MacDonald K.A.
        • McKenna S.L.
        Evaluation of a 3M Petrifilm on-farm culture system for the detection of intramammary infection at the end of lactation.
        Prev. Vet. Med. 2013; 111 (23597619): 1-9
        • Dewes H.F.
        • Goodall G.
        Some preliminary observations on the possible relationship between ammonia production from soiled bedding in calf rearing sheds and calf illness.
        N. Z. Vet. J. 1995; 43 (16031801): 37-41
        • Frössling J.
        • Nöremark M.
        Differing perceptions—Swedish farmers' views of infectious disease control.
        Vet. Med. Sci. 2016; 2 (29067181): 54-68
        • Furman-Fratczak K.
        • Rzasa A.
        • Stefaniak T.
        The influence of colostral immunoglobulin concentration in heifer calves' serum on their health and growth.
        J. Dairy Sci. 2011; 94 (22032377): 5536-5543
        • Gorden P.J.
        • Plummer P.
        Control, management, and prevention of bovine respiratory disease in dairy calves and cows.
        Vet. Clin. North Am. Food Anim. Pract. 2010; 26 (20619182): 243-259
        • Guterbock W.M.
        The impact of BRD: the current dairy experience.
        Anim. Health Res. Rev. 2014; 15 (25351327): 130-134
        • Kaufman J.
        • Linington M.
        • Osborne V.
        • Wagner-Riddle C.
        • Wright T.
        Short communication: Field study of air ammonia concentrations in Ontario dairy calf housing microenvironments.
        Can. J. Anim. Sci. 2015; 95: 539-542
        • Khan M.A.
        • Weary D.M.
        • von Keyserlingk M.A.
        Invited review: Effects of milk ration on solid feed intake, weaning, and performance in dairy heifers.
        J. Dairy Sci. 2011; 94 (21338773): 1071-1081
        • Lago A.
        • McGuirk S.M.
        • Bennett T.B.
        • Cook N.B.
        • Nordlund K.V.
        Calf respiratory disease and pen microenvironments in naturally ventilated calf barns in winter.
        J. Dairy Sci. 2006; 89 (16960078): 4014-4025
        • Love W.J.
        • Lehenbauer T.W.
        • Kass P.H.
        • Van Eenennaam A.L.
        • Aly S.S.
        Development of a novel clinical scoring system for on-farm diagnosis of bovine respiratory disease in pre-weaned dairy calves.
        PeerJ. 2014; 2 (24482759): e238
        • Love W.J.
        • Lehenbauer T.W.
        • Van Eenennaam A.L.
        • Drake C.M.
        • Kass P.H.
        • Farver T.B.
        • Aly S.S.
        Sensitivity and specificity of on-farm scoring systems and nasal culture to detect bovine respiratory disease complex in preweaned dairy calves.
        J. Vet. Diagn. Invest. 2016; 28 (26796957): 119-128
        • Lundborg G.K.
        • Svensson E.C.
        • Oltenacu P.A.
        Herd-level risk factors for infectious diseases in Swedish dairy calves aged 0–90 days.
        Prev. Vet. Med. 2005; 68 (15820112): 123-143
        • McGuirk S.M.
        Disease management of dairy calves and heifers. T.
        Vet. Clin. North Am. Food Anim. Pract. 2008; 24 (18299036): 139-153
        • McGuirk S.M.
        • Peek S.F.
        Timely diagnosis of dairy calf respiratory disease using a standardized scoring system.
        Anim. Health Res. Rev. 2014; 15 (25410122): 145-147
        • Ollivett T.
        Understanding the diagnosis and risk factors for respiratory disease in dairy calves.
        (PhD Thesis) Population Medicine Department, University of Guelph, Guelph, Canada2014
        • Ollivett T.
        • Hewson J.
        • Schubotz R.
        • Caswell J.
        Ultrasonographic progression of lung consolidation after experimental infection with Mannheimia haemolytica in Holstein calves.
        J. Vet. Intern. Med. 2013; 27: 673
        • Ollivett T.L.
        • Buczinski S.
        On-farm use of ultrasonography for bovine respiratory disease.
        Vet. Clin. North Am. Food Anim. Pract. 2016; 32 (26922110): 19-35
        • Ollivett T.L.
        • Caswell J.L.
        • Nydam D.V.
        • Duffield T.
        • Leslie K.E.
        • Hewson J.
        • Kelton D.
        Thoracic ultrasonography and bronchoalveolar lavage fluid analysis in Holstein calves with subclinical lung lesions.
        J. Vet. Intern. Med. 2015; 29 (26332345): 1728-1734
        • Pomés A.
        • Chapman M.D.
        • Wunschmann S.
        Indoor Allergens and Allergic Respiratory Disease.
        Curr. Allergy Asthma Rep. 2016; 16 (27184001): 43
        • Pritchard D.G.
        • Carpenter C.A.
        • Morzaria S.P.
        • Harkness J.W.
        • Richards M.S.
        • Brewer J.I.
        Effect of air filtration on respiratory disease in intensively housed veal calves.
        Vet. Rec. 1981; 109 (7303449): 5-9
        • Seedorf J.
        • Hartung J.
        Survey of ammonia concentrations in livestock buildings.
        J. Agric. Sci. 1999; 133: 433-437
        • Sivula N.
        • Ames T.
        • Marsh W.
        • Werdin R.
        Descriptive epidemiology of morbidity and mortality in Minnesota dairy heifer calves.
        Prev. Vet. Med. 1996; 27: 155-171
        • Stroup W.W.
        Generalized Linear Mixed Models: Modern Concepts, Methods and Applications.
        CRC Press, Boca Raton, FL2012
        • Svensson C.
        • Liberg P.
        The effect of group size on health and growth rate of Swedish dairy calves housed in pens with automatic milk-feeders.
        Prev. Vet. Med. 2006; 73 (16191449): 43-53
        • Svensson C.
        • Lundborg K.
        • Emanuelson U.
        • Olsson S.O.
        Morbidity in Swedish dairy calves from birth to 90 days of age and individual calf-level risk factors for infectious diseases.
        Prev. Vet. Med. 2003; 58 (12706057): 179-197
        • Teixeira A.G.
        • McArt J.A.
        • Bicalho R.C.
        Thoracic ultrasound assessment of lung consolidation at weaning in Holstein dairy heifers: Reproductive performance and survival.
        J. Dairy Sci. 2017; 100 (28215891): 2985-2991
        • Van Donkersgoed J.
        • Ribble C.S.
        • Boyer L.G.
        • Townsend H.G.
        Epidemiological study of enzootic pneumonia in dairy calves in Saskatchewan.
        Can. J. Vet. Res. 1993; 57 (8269363): 247-254
        • Venner M.
        • Astheimer K.
        • Lammer M.
        • Giguere S.
        Efficacy of mass antimicrobial treatment of foals with subclinical pulmonary abscesses associated with Rhodococcus equi.
        J. Vet. Intern. Med. 2013; 27 (23278131): 171-176
        • Windeyer M.C.
        • Leslie K.E.
        • Godden S.M.
        • Hodgins D.C.
        • Lissemore K.D.
        • LeBlanc S.J.
        Factors associated with morbidity, mortality, and growth of dairy heifer calves up to 3 months of age.
        Prev. Vet. Med. 2014; 113 (24269039): 231-240