Unless otherwise stated, statistical results presented here are from multivariable analyses.
3.1. Descriptive data
Assessment of the tutorial feedback forms confirmed the veterinarians had understood the material presented and the protocol details, but revealed wide differences of opinion e.g. over what constituted a maintenance fluid dose rate. There were less than 10 questions for the contact person regarding logistical issues and one question for the investigator regarding concurrent disease and case inclusion.
Based on a randomly chosen subset for data entry accuracy prior to analyses, an error rate of 0.5% was detected when 4190 data points were checked.
A total of 506 cases were originally assessed from 42 practices with an average of 12 cases per practice completed successfully. One practice had very few valid cases (Vic) while others well exceeded their case/time quota. Not all data points were available for all cases due to many reasons e.g. questions omitted, conflicting and obvious errors in the interpretation of terms such as “fluid use” and “type of examination” (e.g. outside a cage versus on a table), data submitted past the closure date, final discharge data missing etc. Missing data were both minor (e.g. dates incorrect) and major (e.g. no therapy data, different presumptions of fluid use – TAS dilution verses hydration).
Complete data sets for 447 cases were available. A total of 59 other cases (11.7%) were initially scored as incomplete, but at various points in the progression of analyses over time, different numbers of data points were validly available from these cases, and so varying n values appear for different analyses. The data were estimated to contain over 100,000 data points and were archived under each of the 42 practices.
With regard to all completed data, a total of 416 of 447 dogs survived. Mortality rate was 6.9% (n=31), including three dogs that were euthanased due to progressive disease. Time to recovery data (n=360) showed early improvement in the majority of cases, but further deterioration after hospitalisation (“lag phase”) was observed in 9% (37) cases before signs of improvement were seen (
Figure 1).
3.2. Analyses
3.2.1. Animal and tick factors
The only variable (animal or tick) that showed any association with mortality in these dogs was age. The median age was 5.4 years (mean 4.5 years) with a range of 0.09 to 17.5 years. The median age of dogs that died was 7 years with 4.5 for survivors. Dogs over 7 years had a greater risk of death (p=0.03), even after adjusting for all concurrent compensated diseases and therapies, including cardiopulmonary disease. There was an 11% increase in the probability of death from years 6 to 7 and a 37% increase for years 6 to 9.
No significant relationship existed between mortality rate and gender (entire or desexed; p=0.89), breed (p=0.86), or body weight (p=0.98), although a trend was seen with heavier dogs taking longer to be discharged at the acceptable recovery level for that practice.
Brachycephalic (non-cross) dogs (n=58) were presented to 33 clinics. A range of brachycephalic breeds (n=12) were involved with the most common being Cavalier King Charles Spaniels and Staffordshire Bull Terriers, with Shi Tzu the next most common. Cases per clinic ranged from 1 to 6 with the average per clinic of 1.7, but the most common being 1 case (in 15 clinics). There was no difference (p=0.86) in deterioration (5.17% vs 5.74%), time to improvement (p=0.86) or time in hospital (p=0.43) for these brachycephalic breeds when compared to all other breeds in this general practice population.
Table 1 outlines the tick factors that showed no association with mortality rate. Tick location was consistent with previous reports [8, 14], with 81% of ticks found from the front legs forward. Owners detected 74% of ticks. Eighty percent of cases involved only one tick, 10% of cases involved two ticks and the remaining 10% of cases involved three or more ticks. Of 455 ticks, 6.6% were dead at presentation. Many (34.6%, 156) measured 4mm dorsal width, 22% measured 3mm and 22% measured 5mm, but mean tick size did not vary between geographical areas. There were no associations between VAS scores and the number or size of ticks present.
3.1.2. Predictors of outcome – NMJ and VAS
Of over 6000 NMJ and VAS scores recorded, univariable analysis showed that VAS respiratory scores were most predictive of a poor prognosis, with the upper quartile (D) having a 36.4% cumulative mortality and an OR of 30.3, relative to the first quartile (
Table 2). The VAS paralysis score had similar predictive results to the NMJ score at all score levels e.g., VAS paralysis D had an OR of 15.5 and a mortality rate of 30.3% whereas the NMJ 4 score had an OR of 14 and a mortality rate of 32.0% (
Table 2). The VAS toxicity assessment had highly significant (p<0.0001) OR scores for the D quartile (19.9), with D group mortality rate at 39.1% (p<0.001).
There were minimal differences in mortality rates between VAS scores A and B and NMJ 1 and 2, but probabilities increased substantially to VAS C and NMJ 3 and then exponentially to VAS D and NMJ 4 i.e. A=B, <C, <<D. This trend applied to all four averaged assessment scores for both mortality rates and OR values (
Table 2).
The most reliable of the four classifications for predicting mortality was VAS respiratory (p=0.005) and the NMJ score (p=0.003). The OR values for both these assessments increased from the lower to the upper quartile; for VAS respiratory they were 1, 4.1, 7.2, 19.3 (p=0.005) and for the NMJ score they were 1, 0.2, 0.8, 4.26 (p=0.003). The conjoined OR for the most severe categories of both assessments (i.e. upper quartile D VAS resp and class 4 NMJ) was 16.6. Thus, a D4 classification had approximately a 16-fold risk of death.
Analysis of the VAS-toxicity assessed at admission by veterinarians produced OR values (A to D) of 1, 15.6, 15.1 and 79.3 with associated mortalities of 0.7, 10.3, 10.0 and 37.8, respectively.
Higher clinical scores (C+D>A; 3+4>1) of the three VAS assessments (Toxic p<0.0001, OR >3.1; Respiratory p<0.004, OR >2.0; Paralysis p<0.0001, OR 3.9) and NMJ scores (p<0.0001, OR 3) were also all associated with longer recovery times.
Figure 2 shows the distribution of disease severity as measured by the four assessment methods. There is some variability between the assessment methods for the three lower quartiles, but the upper quartile was more consistently scored by all four methods. VAS respiratory could clearly differentiate the four categories of severity, reinforcing the effectiveness (as above) of this scoring system to best predict mortality.
3.1.3. Predictors of outcome –specific clinical signs.
Dogs that died had an average of 3.6 ±2.4 clinical signs (with each additional sign increasing the risk of mortality by 23%), whereas those that lived had an average of 2.5±2.0 signs. However, when each of the pre-defined clinical signs (e.g. breathing pattern, polypnoea, etc) were assessed against NMJ scores (
Figure 3), no progressive prognostic pattern could be detected that would aid case predictability. When each clinical sign was assessed in relation to time (i.e., day 1 and day 2), no consistent pattern (
Figure 4A,B) could be seen that would verify the use of the respiratory score previously used in general practice and in other tick studies [
21].
The presence of inspiratory dyspnoea (p=0.002, OR 5.3, CI 1.9-13.0) and crackles (P=0.003, OR 5.7, CI 1.9-15.0) on day 1 were the most significant prognostic indicators of mortality. Mortality was associated with 21.4% of cases that presented with crackles compared to 4.6% of cases that did not so present, and with 19.4% of cases presenting with inspiratory dyspnea compared to 4.4% of cases that did not. Dogs with either of these specific signs had a five-fold risk of death.
On day 2, the presence of expiratory dyspnoea (15% dogs; p=0.0005, OR 4.2, CI 1.5–11.6) and an expiratory wheeze (p=0.0001, OR 9.1, CI 2.9–28.1) were the most significant prognostic indicators of mortality, but the presence of inspiratory or expiratory stridor was not associated with mortality (p=0.19 and p=0.07, respectively).
Retching was associated with both a higher probability of mortality (14% dogs; p=0.02, OR 3.1, CI 1.2-7.3) and longer hospitalisation times (p=0.02, OR 0.38, CI 0.2 – 0.9), while grunting (p=0.05, OR 1.7, CI 0.7-8.7), and gagging (p=0.001, OR.1.92 , CI 0.7-6.6) were associated with delayed recovery.
3.1.4. Predictors of outcome –facial expressions of distress
The results of the analyses of subjective observations of defined facial expressions were profound; anxiety (p=0.003, OR 5; mortality of 18.9% compared to non-anxiety mortality of 4.5%), glazed eyes (p=0.015, OR 5.5; mortality of 22.2% compared to 5.0% non-glazed) and fatigue (p<0.001, OR 10.4; mortality of 31.8% compared to 4.3% non-fatigued) were all highly predictive of death and the OR approximately doubled (P=0.001) for cases with greater than one facial sign (one sign P<0.001, OR 4.9). Increased risk of death was associated with additional facial signs. These facial signs are presumably caused by various aspects of disease expression, and when present, increased the risk of mortality by five to 10-fold and were a reliable indicator of a poor prognosis.
3.1.5. Geographic variation in clinical disease
Data from the 42 sites were widely grouped into four areas – north, south, coastal and inland and analysed to determine whether clinical signs varied between locations. Clinical signs varied (p=0.0068) between areas with respect to the mean VAS toxicity for each area. The frequency of defined respiratory signs also varied between areas. However, overall disease severity (p=0.4) and mortality (p=0.55) did not vary between these areas, and neither did time to improvement (p>0.05).
3.1.6. Pharmaceutical product use and supportive treatments
The only product that was positively associated with mortality was TAS (p<0.001). Mortality rate was not affected by the brand (n=4) of TAS used (p=0.4), the dose / volume administered (p=0.9), or the time of TAS administration (p=0.5 to 0.9). Mean TAS usage in dogs that died versus those that survived was 14.1mL (+/- 7.6) and 14.3 mL (+/- 8.4) (p=0.90), and 1.1 (+/-0.5) and 1.0 (+/-0.5) mL/kg (p=0.73) respectively for TAS dosage. TAS was used as the only therapy in 34 mild cases (VAS resp A-B; NMJ 1-2) and all survived. In a cohort of 407 dogs, TAS use plus one other drug was associated with a mortality of 6.1% (n=25). Three mild cases (VAS resp A, NMJ 1) had no TAS and all survived.
In clinical cases up to NMJ (3) or VAS (C), antibiotic use (but not one specific class) was negatively correlated with mortality (p=0.018, OR 3.08, CI 1.23 – 7.42), as more dogs receiving antibiotics died compared with those that did not receive antibiotics. However, in the most severe classification of VAS Resp, the use of antibiotics was positively correlated to outcome (p=0.018, OR 3.1, CI 1.2 – 9.4). For all four assessment methods, the more severe the case the more likely antibiotics were used (NMJ 4>1 and VAS D>A; 2.7 to 10 times more likely; p=0.02 to <0.001).
Diuretic use (in 13.1% of cases) was also negatively correlated to mortality (p=0.02; OR 3.29; CI 1.2 – 8.02). Of the 51 cases so treated, 44 recovered and seven died (13.7%). In the remaining cases (with a similar drug use profile) that were not treated with diuretics (n=39), only 5% died. However, as with antibiotics, diuretics were more likely to be used in the most severe cases (e.g. VAS Respiratory Score D, 5.4 times more than A; p=0.0014 and NMJ 4, 4.7 times more than NMJ 1; p=0.0071).
Many different pharmaceutical products were used on the 506 dogs. In general, more classes of products were used as severity increased. Risk of death was associated with multiple product use, increasing by 23% with each additional drug class (p=0.003, OR 1.23, CI= 1.04 – 1.43). Although other drug types used (acepromazine, atropine, steroids, antihistamines, antiemetics, and S8 anti-anxiety and sedation drugs) showed no significant associations with mortality (p values from 0.1 to 0.8), the use of antihistamines (p<0.0001), antiemetics (0.015) and diuretics (p=0.001) were all associated with longer hospitalisation times. Analysis of individual cases showed that the use of antiemetics did not necessarily follow clinical indications: antiemetics were used in only 40% of dogs where vomit was observed and were used in 6% of dogs where vomit was not observed.
Similarly, associations with delayed recovery were found for supportive care such as clipping (0.02, OR 1.7, CI 1.1-2.7) eye care (p=0.014, OR 1.7, CI 0.36-1.55) and bladder care (p=0.001, OR 3.7, CI 0.34-7.8) and frequency of examination (p=0.0001, OR 0.55, CI 0.19-1.48).
Fluid use data were hard to assess because the protocol format did not clearly enough define the term, leading to a poor distinction between fluids used to treat dehydration, fluid used for maintenance, for TAS dilution, or used for multiple purposes. Of the 493 dogs assessed, 12% were considered to be dehydrated based on varying combinations of skin turgor, packed cell volume, and total plasma protein, but the degree of dehydration was not recorded. The rate of intravenous fluid administration was noted in 82 cases and varied widely: 16% were administered between one third and half maintenance; 43% were administered between 1.5 and 3 times maintenance and the remainder were delivered at maintenance rates. Fluid therapy to correct dehydration (n=59) significantly improved recovery time (98% dogs; p=0.005; OR 7, CI 1.5–126).
Manual ventilation was used in only 5 of 506 cases, which were too few to analyse. None of the clinics involved had access to intensive care facilities.
Of the cases that were hospitalised, 50% had not received tick preventative product prior to hospitalisation. Acaricide use after hospital admission (n=253) resulted in a 2.2 times (OR) advantage with respect to hospitalisation time (p=0.027; CI 1.0 – 3.9).