Respiratory diseases are recognized health problems in camelids1 and cause for referral to specialized veterinary hospitals. The basis of high-quality care is a firm understanding of the underlying mechanisms of disease. However, a paucity of information exists about respiratory function in camelids in comparison to other domestic animal species.
Respiratory function tests, commonly used in companion animal species, may improve the sensitivity and accuracy of diagnosis of respiratory diseases in camelids. For example, RL and Cdyn are used to describe the level of airway obstruction in horses with recurrent airway obstruction,2 and head-out plethysmography has similar application in conscious dogs.3 In horses with exercise intolerance and cough, airway responsiveness can be determined with a variety of techniques including measurements of RL and Cdyn,2 FOT,4-6 impulse oscillometry,7 or flowmetric plethysmography.8 Respiratory inductance plethysmography has further been used to describe breathing patterns in horses8 and more recently in a llama with dyspnea of unknown origin.9 In the latter report,9 abdominal (diaphragmatic) paradox was described, which corresponded precisely with neurologic deficits localized to the diaphragm causing paralysis.
The purpose of the study reported here was to examine the wider usefulness of respiratory function testing in llamas under conscious sedated and nonsedated conditions. We compared 2 methods of measuring pulmonary mechanics (ie, RL and Cdyn) and evaluated the effects of sedation on these parameters. In addition, we measured the effects of sedation on RIP assessment of breathing pattern in llamas. We hypothesized that the methods to derive RL and Cdyn would be correlated and that sedation would alter respiratory mechanical function and breathing pattern significantly.
Materials and Methods
Animals—The Institutional Animal Care and Use Committee at the Cummings School of Veterinary Medicine at Tufts University approved all procedures. Written client consent was obtained for all study llamas prior to enrollment in the study. Ten healthy client-owned adult llamas (3 male and 7 female) between 2 to 9 years (mean ± SD, 5 ± 2.3 years) of age with a mean body condition score of 5 out of 9 (body weight, 100 to 204 kg) were included in this study. None of the llamas had been previously exposed to pulmonary function testing or a laboratory environment. Following hospital admission in the morning, the llamas were acclimated for 4 to 6 hours indoors. All lung function testing was performed in a temperature-controlled environment (mean, 20°C and 71% humidity).
Study design—A complete history, physical examination, lateral radiographic views of the thorax, and a CBC were obtained prior to study enrollment to eliminate llamas with any evidence of preexisting respiratory disease. The following external body measurements were obtained by use of a flexible measuring tape with the llama in a standing position: height measured at withers, body length from point of shoulder to ischial tuberosity, and chest circumference over the eleventh intercostal space. All procedures were performed with the llama standing while lightly restrained in a commercially available chute.
The following tests were performed to measure pulmonary function: determination of FRC by the rebreathing technique that uses helium dilution as a test gas,10 assessment of breathing pattern (thoracic vs abdominal contribution to ventilation) by use of RIP in conjunction with pneumotachography,8,9 measurement of RL and Cdyn11 during tidal breathing, and the use of monofrequency FOT.12 All tests were performed in unsedated llamas (n = 10) and repeated following sedation in a subset (7) of these llamas for which owner consent was obtained for sedation. For sedation, llamas received xylazine hydrochloride (0.2 mg/kg, IM). All sedated llamas were tested standing while their heads were manually supported to maintain the mandible horizontal to the ground throughout the study period.
FRC determination via helium dilution—Functional residual capacity was measured by use of methods previously described for dogs10 and horses.13 All llamas were fitted with a low–dead space (50 mL), clear plastic face mask that was sealed around the bony portion of the muzzle approximately 10 cm behind the external nares, with a latex shroud. The face mask was connected to a low–dead space 3-way angled tap (120°) stopcocka and nondiffusible gas collection bag containing breathable test gas (20 mL/kg; 0.3% carbon monoxide, 10% helium, 20% O2, and 69.7% nitrogen). The initial and final concentrations of helium and the final concentration of CO2 following a rebreathing period of 45 to 50 seconds were determined by use of specific analyzers.b,c The dilution of helium (assumed to be a nonexchangeable gas) gave a measure of FRC according to the following formula10:


where Hei represents the initial concentration of helium, Hef represents the final concentration of helium, and DSins represents the instrument dead space (ie, bag volume in liters).10 Three consecutive measurements were obtained for each llama and averaged.
RIP and pneumotachography—The concurrent use of RIP and pneumotachography permits noninvasive assessment of breathing pattern and ventilation. The use of RIP has been described in a llama9 and has been used previously in humans,14,15 horses,8,16,17 and sheep.18 In summary, 2 elastic bandsd (4 cm wide, adult size), each containing a single conducting wire embedded in a sinusoidal fashion, were temporarily placed around the thorax (overlying the 11th intercostal space) and abdomen (behind the last rib) of each llama. To eliminate interference of coat, the fiber was separated to allow placement of bands in close proximity to the skin. Normal respiratory movements cause the bands to stretch (inspiration) and relax (expiration), thus altering the inductance properties of coiled wire, which runs through the center of these bands. Inductance decreases during inspiration and increases during expiration. An oscillating current produces a voltage across these bands that is linearly proportional to stretch and therefore proportional to the change in cross-sectional area of the thorax and abdomen. The voltage signals (thoracic from the thoracic band and abdominal from the abdominal band) were collected at 30 Hz, amplified, and digitized for processing by use of commercial data acquisition software.e It is assumed that the thoracic and abdominal contributions to ventilation are separate. Thus these 2 signals can be analyzed independently to evaluate breathing pattern by examining the relative thoracic and abdominal contributions to ventilation.
The analysis of phase angle of thoracic and abdominal signals has been previously described.17 Briefly, the phase angle was measured from the x-y (ie, abdominal vs thoracic) plots by use of the following equation19:


where m represents the width on the x-axis of the x-y plot at midway on y-axis and s represents the maximum width of the plot along the x-axis.19
Additionally, thoracic and abdominal volume signals were summed to obtain a third volume signal. Flow was also derived from RIP signals. The sum of thoracic and abdominal signals was differentiated offline by use of commercial softwaref to obtain the sum flow. The beginning of expiration and beginning of inspiration were defined by the upward- and downward-directed zero crossings of sum flow signals, respectively (Figure 1). Nasal flow was measured concurrently during tidal breathing by use of a heated pneumotachographg attached to the proximal port of the face mask. The pneumotachograph-transducerh was calibrated by use of a precision volume syringe (3-L syringe).a All measurements derived during RIP were based on the assumption that the external flow signal (ie, sum flow) equals the nasal flow signal. Because thoracic, abdominal, and sum signals are collected in an uncalibrated fashion, the gain setting for sum flow was chosen to closely match the volume signal derived from the pneumotachograph through integration (ie, volume = °flow). External flow (ie, sum flow) and nasal flow signals were subtracted in the same time domain to obtain differences in PIF and PEF. Similarly, area differences between these signals during the first 25% of exhaled volume and first 25% of inhaled volume were computed on a breath-by-breath basis. The peak flow and area differences represent the severity of gas compression (during expiration) and expansion (ie, rarefaction [during inspiration]).
RL and Cdyn measurements—Values of RL and Cdyn were derived via the previously described method of covariance ratios11 and from measurements of flow, VT, and transpulmonary pressure as previously described20 during spontaneous breathing in 10 llamas. An esophageal balloon (10 cm in length and 3.8 cm in diameter) sealed over the distal end of a polypropylene catheter (4 mm internal diameter, 5 mm outer diameter, and 90 cm length) was passed through the nares to the midthoracic level of the esophagus and inflated with 2 mL of air. Adequate catheter placement was verified by the identification of dPPlmax. The proximal end of the esophageal balloon catheter was exited through a hole within the proximal aspect of the face mask and connected to a differential pressure transducer.h Transpulmonary pressures were estimated by measuring the differential between face mask and esophageal pressures during spontaneous breathing. The nosepiece of the face mask was attached to a heated pneumotachograph.g Flow was thus detected during tidal breathing and electronically integrated to determine VT on a breath-by-breath basis. The analog signals derived from the pressure transducer and pneumotachograph were collected at 30 Hz, amplified, and digitized by use of pulmonary mechanics analysis software.i Measurements of flow, VT, expired minute volume, frequency, PIF, PEF, inspiratory time, expiratory time, and dPPlmax were recorded and displayed continuously. Flow, VT, and transpulmonary pressure were used to calculate RL by use of the isovolume method and Cdyn at 2 points of zero flow as described by Amdur and Mead.20
Monofrequency FOT—Monofrequency FOT (at 1 to 7 Hz) was used to determine respiratory system impedance, RRSosc, and respiratory system reactance, as well as fres as previously described.4 The RRSosc is comprised of the additive effects of chest wall resistance and RL.
Measurements obtained via FOT were conducted with the already described face mask and esophageal balloon catheter, but a different size pneumotacho-graph.j The nosepiece of the face mask was affixed to the pneumotachograph, with a more rostral T-piece containing 4 ports as follows: a small port for externally derived oscillating airflow (waveform sinusoidal, 1 to 7 Hz frequencies); a large side port for spontaneous breathing through a resistive element (2 cm H2O/L/s); a large port for the pneumotachograph; and, a lateral port (4 mm internal diameter) for measurement of airway opening pressure by use of a differential pressure transducer.
A sinusoidal airflow of desired frequency was generated by use of a 3-port proportional pneumatic valvek connected to a compressed air source (75 psi of pressure). Face mask pressure relative to atmospheric pressure (for measurement of respiratory system impedance) or relative to esophageal (pleural) pressure (for measurement of lung impedance) were determined by use of a differential pressure transducer.h Airflow to the face mask was measured by use of a pneumotachographj and differential pressure transducer.h Amplified pressure and flow signals were digitized at 25.6 Hz for 10 seconds. Signals were bandpass filtered (0.4-Hz-wide passband centered at the measurement frequency and 80-dB stopband attenuation) and divided into consecutive 5-second data segments with 50% overlap, from which resistance, reactance, and impedance were calculated by use of a personal computer, purpose-built controller, and digital processing system.l A coherence value was calculated to provide an indication of signal-to-noise ratio and linearity of the system. Only values with a coherence ≥ 0.9 were used for analysis.
The ZRS was calculated on the basis of the magnitude and phase shifts between pressure and flow sinusoidal waves at a given time and frequency as follows21:


where φ represents the phase shift between flow and pressure waves, Pao is pressure at the airway opening, V'ao is flow at the airway opening, ω represents angular frequency in radians, and j represents √−1. The phase angle was used via computer analysis software to calculate the components of respiratory system or pulmonary impedance, resistance, and reactance. Respiratory system resistance and reactance (XRS) were computed as follows21:


The forced oscillation system was calibrated each testing day by use of the wave tube principle.7 An open-ended polyvinyl chloride pipe (6.2 m long, 52.3 mm internal diameter) with a known impedance (Zref) was compared with the measured impedance by use of the forced oscillatory system. A correction factor (k) was derived from this comparison (k = Zref/ZRS) and used to correct the measured impedance (corrected impedance = kZRS) of llamas.
The fres was defined as the impulse frequency at which impedance measurements were entirely comprised of resistance (ie, reactance was equal to zero). By use of the components of reactance, inertance (I) and compliance (C), fres can be calculated as follows4:


Statistical analysis—Descriptive data analyses of respiratory mechanical parameters (FRC, VT, inspiratory time, expiratory time, expired minute volume, PIF, PEF, RL Cdyn, dPPlmax, RRSosc, fres, and RIP-derived variables) were reported as mean ± SD. The correlation between FRC and body dimensions (weight, height, length, and circumference) as well as between measurements obtained via FOT and esophageal balloon mechanics was evaluated by use of the Spearman rank correlation coefficient (R). The Bland-Altman test of agreement22 was applied to test agreement between measurements of Cdyn and RL obtained by use of FOT and esophageal balloon mechanics. These results were reported for Cdyn and RL as mean difference between the 2 techniques with the 95% confidence interval. Outlier tests were detected via extreme values analysis. The effect of sedation was evaluated by use of a paired samples t test, with an accepted significance level of P < 0.05. All analyses were performed by use of commercial software.m
Results
Physical characteristics of llamas and hematologic results—Measurements of weight (150 ± 40.3 kg), height at the withers (113.9 ± 12.1 cm), body length from point of the shoulder to ischial tuberosity (108.7 ± 12.6 cm), and chest circumference (135.5 ± 18.6 cm) were within the reference range for healthy adult llamas.23 Similarly, hematologic analysis revealed a mean total WBC count of 12.76 × 103 cells/μL (range, 9.5 × 103 cells/μL to 19.4 × 103 cells/μL) and mean Hct of 26% (range, 24% to 34%), values which were within reference range.24 Thoracic radiography did not reveal any changes consistent with respiratory disease in the study llamas.
Pulmonary mechanics in nonsedated llamas— Statistical analysis revealed a normal distribution for all data for nonsedated and sedated llamas. Baseline data for ventilatory parameters before sedation in 10 llamas were determined (Table 1). Mean FRC in nonsedated llamas was 5.60 ± 1.24 L, which translated to 39.18 ± 11.45 mL/kg for this group. A positive correlation was observed between FRC and body length in our study llamas (P = 0.016; ρ = 0.733). However, FRC values were not significantly correlated with body weight, height, or chest circumference. Breathing patterns were assessed by RIP (Figure 1). Peak differences between RIP and pneumotachographic flows were 2.48 ± 1.04 L/s during inspiration (ie, differences in PIF), and 2.61 ± 1.02 L/s during expiration (ie, differences in PEF). When flow was integrated over the first 25% of the inspired or expired breath, the derived RIP and pneumotachographic volume differences were 0.05 ± 0.04 L and 0.04 ± 0.14 L, respectively.
Mean ± SD baseline lung mechanic values of nonsedated adult llamas (n = 10).
Variable | Nonsedated llamas |
---|---|
f (breath s/min) | 42.38 ± 12.12 |
VT(L/s) | 1.03 ± 0.3 |
PIF(L/s) | 2.14 ± 0.65 |
PEF(L/s) | 2.01 ± 0.72 |
PEF:PIF | 0.95 ± 0.24 |
Te(s) | 0.89 ± 0.36 |
Ti(s) | 0.67 ± 0.15 |
Te:Ti | 1.31 ±0.4 |
VE (L/min) | 40.84 ± 12.95 |
f = Respiratory frequency. Ti = Inspiratory time. Te = Expiratory time. VE = Minute ventilation.

The RIP and pneumotachographic signals. A—Diagram depicting combined matching nasal flow (V'n; blue line) and sum flow (V'sum; red line) signals derived from elastic band measurements (sum flow = abdominal + thoracic components) in a healthy llama. The beginning of expiration (E) and beginning of inspiration (I) were defined by the upward- and downward-directed zero crossings (over horizontal double arrowheads) of sum flow signals, respectively. B—Thoracic (RIB) versus abdominal (ABD) contribution to breathing in a healthy llama. Both signals are completely in phase, and the llama has a much larger abdominal relative to thoracic contribution to ventilation.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676

The RIP and pneumotachographic signals. A—Diagram depicting combined matching nasal flow (V'n; blue line) and sum flow (V'sum; red line) signals derived from elastic band measurements (sum flow = abdominal + thoracic components) in a healthy llama. The beginning of expiration (E) and beginning of inspiration (I) were defined by the upward- and downward-directed zero crossings (over horizontal double arrowheads) of sum flow signals, respectively. B—Thoracic (RIB) versus abdominal (ABD) contribution to breathing in a healthy llama. Both signals are completely in phase, and the llama has a much larger abdominal relative to thoracic contribution to ventilation.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676
The RIP and pneumotachographic signals. A—Diagram depicting combined matching nasal flow (V'n; blue line) and sum flow (V'sum; red line) signals derived from elastic band measurements (sum flow = abdominal + thoracic components) in a healthy llama. The beginning of expiration (E) and beginning of inspiration (I) were defined by the upward- and downward-directed zero crossings (over horizontal double arrowheads) of sum flow signals, respectively. B—Thoracic (RIB) versus abdominal (ABD) contribution to breathing in a healthy llama. Both signals are completely in phase, and the llama has a much larger abdominal relative to thoracic contribution to ventilation.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676
Measured parameters of pulmonary mechanics (RL, Cdyn, dPPlmax, RLosc, and RRSosc) were determined (Table 2). Values of RL (conventional RL and RLosc) were significantly correlated (P = 0.025; ρ = 0.733) at 1 Hz (the input frequency closest to the spontaneous breathing frequency). The Bland-Altman test of agreement between RL and RLosc revealed a mean difference of 0.83 ± 0.69 cm H2O/L/s with a 95% confidence interval of 0.3 to 1.35 cm H2O/L/s. The Cdyn obtained via conventional pulmonary mechanics was compared with an index of Cdyn calculated from XLosc values at 1 Hz (XLosc values at 1 Hz, 0.54 ± 0.39; n = 9) as previously described.13 These compliance values (ie, Cdyn and Cdynosc) were positively correlated (P = 0.05; ρ = 0.667). However, mean compliance calculated from XLosc (Cdynosc; 0.52 ± 0.43 L/cm H2O) was significantly (P = 0.05) lower than mean Cdyn measured via conventional pulmonary mechanics. The Bland-Altman test of agreement between the calculated compliance and measured Cdyn revealed a mean difference of 0.43 ± 0.22 L/cm H2O with a 95% confidence interval of 0.26 to 0.60 L/cm H2O.
Mean ± SD respiratory function variables obtained via conventional pulmonary mechanics and FOT in nonsedated adult llamas.
Method | Variable | Oscillation frequency (Hz) | No. of Ilamas | Non-sedated of Ilamas |
---|---|---|---|---|
Conventional* | dPPImax (cm H2O) | NA | 10 | 6.55 ± 2.06 |
RL (cm H2O/L/s) | NA | 10 | 1.42 ± 0.54 | |
Cdyn (L/cm H20) | NA | 10 | 0.83 ± 0.40 | |
FOT† | RLosc(cm H2O/L/s) | 1 | 9 | 2.24 ± 0.90 |
2 | 10 | 2.08 ± 0.72 | ||
3 | 10 | 2.13 ± 0.61 | ||
5 | 10 | 2.33 ± 0.68 | ||
7 | 10 | 2.56 ± 0.39 | ||
FOT† | RRSosc (cm H2O/L/s) | 1 | 8 | 2.40 ± 0.93 |
1 | 8 | 2.31 ± 0.71 | ||
2 | 10 | 2.20 ± 0.62 | ||
3 | 10 | 2.71 ± 0.66 | ||
5 | 10 | 2.50 ± 0.47 |
Use of esophageal balloon and pneumotachography.
Only values with a coherence ≥ 0.9 were used for analysis via FOT.
NA= Not applicable.
Pulmonary mechanics in sedated llamas—The mean FRC in sedated llamas (n = 5) was 4.17 ± 0.60 L with a calculated FRC/kg of 30.35 ± 10.56 mL/kg. This represents a significant (P = 0.02) mean decrease in FRC of 0.83 ± 0.23 L (17%) following sedation with xylazine. Differences in PEF (n = 6) increased by 76% following sedation, with a mean difference in PEF of nonsedated llamas of 2.61 ± 1.02 L/s versus differences in PEF of sedated llamas of 4.99 ± 2.8 L/s. Differences in PIF values (n = 6) increased by 34% following sedation, with a mean differences in PIF of nonsedated llamas of 2.47 ± 1.04 L/s versus differences in PIF of sedated llamas of 3.74 ± 1.69 L/s.
Effects of sedation with xylazine on ventilatory parameters and tidal breathing mechanics were determined (Tables 3 and 4). Although measurements obtained via FOT and conventional pulmonary mechanics were attempted in 7 sedated llamas for which owner consent was available, results were excluded for 1 llama that did not remain standing following sedation with xylazine and 1 llama that developed epistaxis resulting from passage of the nasoesophageal balloon during conventional pulmonary mechanics. Data were also excluded if acceptable test coherence (ie, signal-to-noise ratio) was not met during the performance of FOT. This may be related to signal interference associated with the llama's normal respiratory frequency, glottic closure, cardiogenic oscillations, or abnormal signal transduction.
Mean ± SD lung mechanic values of adult llamas (n = 7) before (nonsedated) and after (sedated) IM administration of xylazine.
Variable | Llamas | ||
---|---|---|---|
Nonsedated | Sedated | P value | |
f(breaths/min) | 41.32 ± 12.9 | 17.67 ± 6.05 | 0.004 |
VT(L/s) | 1.06 ± 0.19 | 1.45 ± 0.2 | 0.002 |
PIF (L/s) | 2.24 ± 0.56 | 1.32 ± 0.48 | < 0.001 |
PEF(L/s) | 2.17 ± 0.84 | 1.27 ± 0.5 | 0.019 |
PEF:PIF | 0.95 ±0.27 | 0.97 ± 0.27 | 0.917 |
Te(s) | 0.96 ±0.45 | 2.08 ± 0.48 | 0.004 |
Ti(s) | 0.66 ± 0.14 | 1.57 ± 0.38 | 0.001 |
Te:Ti | 1.39 ± 0.49 | 1.35 ± 0.13 | 0.854 |
VE (L/min) | 43.95 ± 14.42 | 24.62 ± 9.31 |
See Table 1 for key.
Mean ± SD respiratory function variables obtained via conventional pulmonary mechanics and FOT in adult llamas before (nonsedated) and after (sedated) IM administration of xylazine.
Method | Variable | Oscillation frequency (Hz) | Llamas | P value | |||
---|---|---|---|---|---|---|---|
No. | Nonsedated | Sedated | Change (%) | ||||
Conventional* | dPPImax (cm H2O) | NA | 6 | 6.34 ± 1.24 | 8.80 ± 2.79 | 46 | 0.030 |
RL (cm H2O/L/s) | NA | 6 | 1.29 ± 0.5 | 3.18 ± 1.65 | 147 | 0.015 | |
Cdyn (L/cm H2O) | NA | 5 | 0.96 ± 0.44 | 0.62 ± 0.15 | -33 | 0.119 | |
FOT† | RLosc(cm H2O/L/s) | 1 | 6 | 2.04 ± 0.97 | 5.39 ± 2.45 | 175 | 0.021 |
2 | 7 | 1.92 ± 0.69 | 4.34 ± 1.91 | 130 | 0.005 | ||
3 | 7 | 1.95 ± 0.54 | 3.99 ± 2.18 | 176 | 0.034 | ||
5 | 7 | 2.08 ± 0.51 | 3.76 ± 1.49 | 80 | 0.009 | ||
7 | 7 | 2.57 ± 0.43 | 3.78 ± 1.1 | 53 | 0.023 | ||
FOT† | RRSosc (cm H2O/L/s) | 1 | 5 | 2.14 ± 0.84 | 5.14 ± 2.91 | 181 | 0.112 |
2 | 7 | 2.14 ± 0.61 | 5.39 ± 2.86 | 161 | 0.022 | ||
3 | 7 | 2.10 ± 0.64 | 4.41 ± 1.85 | 118 | 0.015 | ||
5 | 7 | 2.55 ± 0.58 | 4.06 + 1.28 | 64 | 0.019 | ||
7 | 7 | 2.51 ± 0.54 | 4.03 ± 1.02 | 55 | 0.019 |
See Table 2 for key.
A significant increase was found in dPPlmax (P = 0.03) and RL (P = 0.015) in sedated llamas. Although differences in Cdyn between sedated and nonsedated llamas were not significant, a decreased Cdyn was observed in sedated llamas. Sedation similarly induced a significant increase in the mean RLosc over all frequencies. It was also observed that sedation caused frequency dependence of RLosc whereby lower frequency resistances increased to a greater extent than higher frequency resistances. Mean fres in nonsedated llamas was 3.37 ± 0.68 Hz (n = 10), whereas fres could only be calculated for 3 sedated llamas. For these llamas, mean fres was 5.18 ± 2.37 Hz. For the other llamas (n = 4), respiratory system reactance values did not approach zero at frequencies between 1 to 7 Hz, and thus the influences of reactance on impedance could not be ruled out (Figure 2).

Box plots of XLosc versus oscillatory frequency in xylazine-sedated (white boxes) and nonsedated (gray boxes) llamas. Box represents the 25th to 75th percentile, horizontal line represents the median, and whiskers represent the range.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676

Box plots of XLosc versus oscillatory frequency in xylazine-sedated (white boxes) and nonsedated (gray boxes) llamas. Box represents the 25th to 75th percentile, horizontal line represents the median, and whiskers represent the range.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676
Box plots of XLosc versus oscillatory frequency in xylazine-sedated (white boxes) and nonsedated (gray boxes) llamas. Box represents the 25th to 75th percentile, horizontal line represents the median, and whiskers represent the range.
Citation: American Journal of Veterinary Research 68, 6; 10.2460/ajvr.68.6.676
Discussion
To our knowledge, this study presents for the first time the application of respiratory mechanical function tests in adult llamas. Some of the respiratory function tests discussed here have been used in other domestic species, results of which warrants comparison where possible. The ease with which respiratory function tests were performed in conscious llamas may provide an opportunity to use these tests as ancillary diagnostic tools. However, the significant alteration caused by sedation with xylazine may prohibit the use of A2-adrenoreceptor agonists for this purpose.
The mean FRC per kilogram of body weight in llamas from our study (40 mL/kg; 5.60 ± 1.24 L) falls within the range of values obtained in horses13,25,26 (44 mL/kg), cattle25 (39 mL/kg), dogs3,10,27 (34 to 98 mL/ kg), and sheep28-30 (25 to 52 mL/kg). Numerous factors influence FRC including abdominal fill, thoracic wall versus lung compliance, thorax size, active component of expiration, posture, and position.10,28 In this study, a positive correlation (P = 0.016) between body length and FRC was observed but not between FRC and body weight or circumference. Functional residual capacity has previously been positively correlated to body weight in horses.13 Earlier studies evaluating FRC in dogs10 and sheep29 found positive correlations between FRC and body weight, circumference, and length. Variation in body condition may explain some of the failure of allometric scaling on the basis of body weight and shape within species. The most extreme example of how body weight fails to predict FRC is in dogs, where body shape varies tremendously with breeds.3,10,27
Changes in FRC may serve as an index of pulmonary mechanical dysfunction. For example, loss of lung volume caused by pneumonia or space occupying masses in the chest, or decreased lung compliance secondary to pulmonary fibrosis, may all reduce FRC. In our study we demonstrate that FRC is obtained easily in conscious llamas.
Llamas have a greater abdominal (ie, diaphragmatic and abdominal muscle) than rib cage (ie, intercostal muscle) contribution to ventilation at rest, as evidenced by their larger abdominal signal, compared with thoracic signal. These findings are in contrast to those for humans and healthy horses. In humans, the contribution of rib cage motion predominates over abdominal movement during quiet breathing.31 Healthy horses have a relatively equal thoracic to abdominal contribution to ventilation at rest.16 This discrepancy may be explained by the relatively smaller chest size and possibly greater passive motion of the llama rib cage during respiration, in comparison with the adult horse.16 Abdominal and thoracic signals were synchronous during inspiration and expiration in all our study llamas, as expected for healthy llamas. Previous publications in veterinary medicine in which RIP was used include the measurement of breathing pattern during exercise,32 recurrent airway obstruction (heaves) and histamine challenge,8 proximal and distal airway obstructions evoked in foals,17 and diaphragmatic paralysis in a young llama.9 Results of our study expand our knowledge about the normal breathing pattern in llamas, which shows a remarkable dominance by the abdominal compartment. This may suggest that quiet breathing is dominated by the diaphragm, but does not imply that this relationship will persist during periods of increased ventilatory demand. Future studies are warranted to examine the llama breathing pattern under conditions of resistive load and hyperpnea.
In our study, mean VT in llamas was 7 mL/kg (1.03 L). This is less than the mean VT measured in horses (12 mL/kg)25,33,34 and dogs (12 mL/kg),3 but falls within the range observed in domestic ruminants such as cattle (8 mL/kg)25,34 and sheep (6 to 7 mL/kg).28 Variation in VT across species is related to absolute differences in lung volume as well as other anatomic features such as size and shape of the abdomen and relative abdominal mass25 and metabolic demand. For example, although similar in body size to the horse, cattle have a smaller lung volume and VT relative to body size.25 Although total lung volume was not measured in our study, llamas share certain anatomic similarities to ruminants because of their large abdominal compartment relative to thoracic compartment. Therefore, it could be speculated that VT relative to body size in llamas would be more similar to ruminants than monogastrics.
Measurements of RL and Cdyn in llamas are comparable to values observed in other species once corrections for differences in body size, lung volume, and breathing patterns have been made. Values of RL in llamas from our study are greater than those found previously in horses (1.16 ± 0.65 cm H2O/L/s) and less than values recorded in sheep (2.19 cm H2O/L/s). The RL is expected to decrease with increasing body size as a result of an overall increase in airway diameter,35 and comparisons of specific RL (RL and FRC) among llamas, horses, and sheep are similar.
In contrast, mean Cdyn in llamas is less than mean values recorded in horses (1.76 ± 0.55 L/cm H2O),5 and calculated specific Cdyn (Cdyn and FRC) is also lower than in horses.25 As Cdyn is related to lung volume, it should scale directly to body size. However, it is also determined by chest wall compliance. Results of a previous study comparing cattle and horses of similar body size also found a significantly lower specific Cdyn in cattle, compared with horses.25 In that study, relative to horses, cattle have decreased chest wall compliance (1.73 L/cm H2O vs 1.47 L/cm H2O) and smaller lung volumes. The authors attributed these differences to adaptations in response to the large ruminal contents in the abdomens of cattle. Other adaptations in cattle to decrease work of breathing in the face of decreased Cdyn include a smaller VT and increased respiratory frequency, compared with horses. Llamas are not true ruminants; however, the gastrointestinal anatomy and bulk within the abdomen of llamas is similar to ruminants, and measurements of VT and respiratory frequency in llamas are more similar to sheep (7 mL/kg29) and cattle (8 mL/kg25), as opposed to horses (12 mL/kg13). Although speculative, it is possible that the decreased relative Cdyn in llamas, compared with horses, is secondary to differences in chest wall compliance and anatomy similar to what is found in cattle and small ruminants.
A significant correlation (P = 0.025; ρ = 0.73) between RL measured by noninvasive FOT (ie, RLosc) and the conventional, isovolume method was observed, as expected. However, the mean RLosc at 1 Hz (2.0 cm H2O/L/s) was higher than the llamas' conventional RL (1.4 cm H2O/L/s), measured during spontaneous breathing frequencies (mean, 0.7 Hz). Although these results complement a report comparing RRSosc to RL in horses,4 this observation stands in general contrast to expected higher resistance values obtained via conventional techniques because RL tends to increase at lower frequencies. Measurements of Cdynosc and conventional Cdyn were also significantly correlated (P = 0.05; ρ = 0.67). However, values obtained via FOT were lower, compared with conventional measurements, with a mean difference of 0.43 ± 0.22 L/cmH2O (95% confidence interval, 0.26 to 0.6 L/cmH2O). A similar close correlation but disparity in absolute agreement between forced oscillation and conventional techniques has been previously reported in horses4 and calves36 and attributed to biological factors, including frequency dependence of compliance. Values of Cdynosc in the study reported here were obtained at 1 Hz, which is higher than the spontaneous breathing frequencies of our llamas, and might partially account for the differences between Cdynosc and conventional Cdyn measurements. Additionally, a decrease in Cdyn and increase in RL may reflect presence of airway obstruction and consequent unequal time constants in the lung. Although all llamas tolerated the pulmonary function procedures well, they were not accustomed to FOT and may have experienced glottic narrowing in response to airflow oscillation, similar to reports in humans and horses.4
Llamas had a higher fres (3.37 ± 0.68 Hz) relative to values reported in the horse (2.46 ± 0.02 Hz)6 and lower mean fres, compared with values recorded in calves37 (5 to 6 Hz), Beagles38 (6.17 ± 0.5 Hz), and humans39 (6 to 11 Hz). This is expected as a result of differences in body size where fres increases, with decreases in body size similar to increases in respiratory resistance. Results of numerous studies have previously revealed that fres increases in response to bronchoconstrictive respiratory diseases such as human chronic obstructive pulmonary disease and asthma40,41 and equine heaves5,8,12 and is related to higher airway resistance.
Sedation with xylazine significantly altered ventilatory parameters in llamas. Respiratory frequency and minute ventilation decreased, in conjunction with increased VT, a reduction in peak flows (ie, PIF and PEF), and prolongation of inspiratory and expiratory times, as found in horses.42,43 Measurements obtained via conventional pulmonary mechanics further revealed significant increases in RL (P = 0.013) and dPPlmax (P = 0.02), which was corroborated via FOT. The Cdyn decreased in sedated llamas, although this change was not significant. The effect of xylazine on respiratory mechanical function has been previously evaluated in studies on horses,42-45 sheep,46,47 and calves.37 Potential adverse effects of sedation with xylazine include increased upper airway resistance, decreased respiratory rate and minute volume, increased VT, bronchodilation or bronchoconstriction, and increased work of breathing.21,37,42-47 In several earlier studies,42,43 horses had increased upper airway resistance, primarily caused by nasal edema developing secondary to vascular congestion in the nares when the head was allowed to drop during sedation.42,43 Many of the adverse effects of xylazine have since been attenuated by controlling head and neck position and administering the lowest possible dose for sedation.21 Because of the anticipated potential for nasal edema in llamas, similar precautions were performed in our study.
Sheep appear to have a response to xylazine comparable to that of llamas, including evidence of bronchoconstriction, increased dPPlmax, decreased Cdyn, and increased RL.46,47 Similarly, in xylazine-sedated calves,37 analysis via FOT revealed a higher fres and RRSosc, compared with controls, with RRSosc demonstrating negative frequency dependence (ie, RRSosc was highest at lower frequencies) as a result of increased peripheral airway resistance. These reported results in ruminant species are qualitatively the same as our findings in llamas, where airway narrowing (increased resistance) was evidenced by significantly increased RLosc, negative frequency dependence of RRSosc, and increased fres in sedated llamas for which calculation was possible (n = 3 for fres). The 4 llamas for which fres could not be calculated under sedation had higher RL and RRSosc values, compared with the other llamas. Results of a previous study48 evaluating FOT in swine and calves revealed similar difficulties in calculating fres that were attributed to poor model fit secondary to shunting from increased face mask dead space and interference by soft tissue structures of the nasopharynx.Thus it is possible that in several of our sedated llamas, the degree of increased RRSosc impeded measurement of accurate impedance values (and thus fres). The magnitude of the observed airway closure in healthy sedated llamas strongly suggests that sedation with xylazine is contraindicated during respiratory function testing in llamas. On the basis of the excellent cooperation of untrained nonsedated llamas in our study, sedation should not be necessary, assuming proper chute restraint is available. Individual case reports have also reported good test compliance in sick llamas,9 although further experience needs to be gained under clinical conditions.
The respiratory mechanical function of llamas in our study was comparable to other animal species, and observed differences may be attributed to variations in body size or conformation. However, further evaluation of anatomic adaptations to a dry, high-altitude environment (ie, increased minute ventilation, narrowing of the nasal passage, and long neck) based on the native origin of camelids should be considered. In addition, a concise understanding of the gross and histologic lung anatomy of llamas and how this may influence the type and distribution of lung disease is important. Clear anatomic differences in lung anatomy exist between species. Although the gross lung anatomy of llamas is most similar to horses,49 a description of the histologic anatomy is currently not available. Grossly, equine lungs have simple lobation, compared with the well-developed lobe divisions found in cattle, sheep, pigs, and dogs. Furthermore, cattle, pigs, and sheep have well-developed lung lobulation, which is incomplete in horses and nonexistent in dogs.50 Because cattle and pigs lack collateral ventilation, lung disease, atelectasis, and airway obstruction generally localize to lobules. In horses with incomplete lobulation and dogs with absent lobulation, collateral ventilation is more developed and lung disease more diffusely distributed.50
In conclusion, noninvasive pulmonary function tests (ie, FOT and RIP) were well tolerated in untrained nonsedated llamas and validated through comparison to conventional (gold standard) techniques. Noninvasive pulmonary function testing has extensive clinical applications in the diagnosis and treatment of respiratory disease in llamas, although testing should not be performed with the use of sedation with xylazine to avoid bronchoconstriction.
ABBREVIATIONS
RL | Pulmonary resistance |
Cdyn | Dynamic compliance |
FOT | Forced oscillation technique |
RIP | Respiratory inductance plethysmography |
FRC | Functional residual capacity |
PIF | Peak flow during inspiration |
PEF | Peak flow during expiration |
VT | Tidal volume |
dPPlmax | Maximum change in transpulmonary pressure |
RRSosc | Respiratory system resistance measured via forced oscillation technique |
fres | Resonant frequency |
ZRS | Respiratory system impedance |
RLosc | RL measured via forced oscillation technique |
XLosc | Pulmonary reactance measured via forced oscillation technique |
Cdynosc | Dynamic compliance measured via forced oscillation technique |
Hans Rudolph, Kansas City, Mo.
He analyzer, PK Morgan Ltd, Chatham, Kent, England.
CO2 analyzer modelCD-3A, Amtek, Pittsburgh, Pa.
Respibands, Ambulatory Monitoring Systems, Ardsley, NY.
XA Biosystem, Buxco Electronics, Sharon, Conn.
Acknowledge, version 3.7, Biopac, Goleta, Calif.
Fleisch No. 3, OEM Medical, Lenoir, NC.
DP45-18, Validyne Engineering, Northbridge, Calif.
Buxco Biosystems XA, Buxco Electronics Inc, Sharon, Conn.
Fleisch No. 4, OEM Medical, Lenoir, NC.
Three-port proportional valve No. 602 00001, Joucomatic, Rueil, France.
On the Nose, Scientific Solutions, Eden Mills, ON, Canada.
SPSS, version 12, SPSS Inc, Chicago, Ill.
References
- 2↑
Derksen FJ, Robinson NE, Armstrong PJ, et al. Airway reactivity in ponies with recurrent airway obstruction (heaves). J Appl Physiol 1985;58:598–604.
- 3↑
Bedenice D, Rozanski E, Bach J, et al. Canine awake head-out plethysmography (HOP): characterization of external resistive loading and spontaneous laryngeal paralysis. Respir Physiol Neurobiol 2005;151:61–73.
- 4↑
Young SS, Tesarowski D. Respiratory mechanics of horses measured by conventional and forced oscillation techniques. J Appl Physiol 1994;76:2467–2472.
- 5↑
Mazan MR, Hoffman AM, Manjerovic N. Comparison of forced oscillation with the conventional method for histamine bronchoprovocation testing in horses. Am J Vet Res 1999;60:174–180.
- 6↑
Hoffman AM, Mazan MR, Ellenberg BS. Association between bronchoalveolar lavage cytologic features and airway reactivity in horses with a history of exercise intolerance. Am J Vet Res 1998;59:176–181.
- 7↑
van Erck E, Votion DM, Kirschovink N, et al. Use of impulse oscillometry system for testing pulmonary function during methacholine bronchoprovocation in horses. Am J Vet Res 2003;64:1414–1420.
- 8↑
Hoffman AM, Riedelberger K, Kupcinskas K, et al. Flowmetric comparison of respiratory inductance plethysmography and pneumotachography in horses. J Appl Physiol 2001;91:2767–2775.
- 9↑
Bedenice D, Mazan MR, Kuehn H, et al. Diaphragmatic paralysis due to phrenic nerve degeneration in a llama. J Vet Intern Med 2002;16:603–606.
- 10↑
Amis TC, Jones HA. Measurement of functional residual capacity and pulmonary carbon monoxide uptake in conscious Greyhounds. Am J Vet Res 1984;45:1447–1450.
- 11↑
Roy R, Powers SR, Kimball WR. Estimation of respiratory parameters by the method of covariance ratios. Comput Biomed Res 1974;7:21–39.
- 12↑
Young SS, Tesarowski D, Viel L. Frequency dependence of forced oscillatory respiratory mechanics in horses with heaves. J Appl Physiol 1997;82:983–987.
- 13↑
Hoffman A, Swanson L, Bruns S, et al. Effects of tension of the girth strap on respiratory system mechanics in horses at rest and during hyperpnea induced by administration of lobeline hydrochloride. Am J Vet Res 2005;66:1167–1174.
- 14
Cohen KP, Panescu D, Booske JH, et al. Design of an inductive plethysmograph for ventilation measurement. Physiol Meas 1994;15:217–229.
- 15
Cohn MA, Rao ASV, Broudy M, et al. The respiratory plethysmograph: a new noninvasive monitor of respiration. Bull Eur Physiopathol Respir 1982;18:643–658.
- 16↑
Korterba A, Kosch PC, Beech J, et al. Breathing strategy of the adult horse (Equus caballus) at rest. J Appl Physiol 1988;64:337–346.
- 17↑
Miller C, Hoffman A, Hunter J. Thoracoabdominal asynchrony failed to grade airway obstruction in foals. J Appl Physiol 2001;88:2081–2087.
- 18↑
Abraham WM, Watson H, Schneider A, et al. Noninvasive ventilatory monitoring by respiratory inductive plethysmography in conscious sheep. J Appl Physiol 1981;51:1658–1661.
- 19↑
Allen JL, Greenspan JS, Deoras KS. Interaction between chest wall motion and lung mechanics in normal infants and infants with BPD. Pediatr Pulmonol 1991;11:37–43.
- 20↑
Amdur MO, Mead J. Mechanics of respiration in unanesthetized guinea pigs. Am J Physiol 1958;192:364–368.
- 21↑
Hoffman AM. Clinical application of pulmonary function testing in horses. International Veterinary Information Service Web site. In:Lekeux P, ed.Equine respiratory diseases. Ithaca, New York: International Veterinary Information Service, 2002. Available at: www.ivis.org. Accessed Jan 24, 2007.
- 22↑
Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986;1:307–310.
- 23↑
Fowler ME. Clinical diagnosis: examination and procedures. In:Fowler ME, ed.Medicine and surgery of South American camelids. 2nd ed.Ames, Iowa: Iowa State University Press, 1998;69–88.
- 24↑
Fowler ME. Hemic and lymphatic systems. In:Fowler ME, ed.Medicine and surgery of South American camelids. 2nd ed.Ames, Iowa: Iowa State University Press, 1998;364–366.
- 25↑
Gallivan GJ, McDonell WN, Forrest JB. Comparative pulmonary mechanics in the horse and the cow. Res Vet Sci 1989;46:322–330.
- 26
Aguilera-Tejero E, Pascoe J, Smith BL, et al. Evaluation of a technique for detection of pulmonary hemorrhage in horses, using carbon monoxide uptake. Am J Vet Res 1994;55:1032–1036.
- 27
Dublin S, Wescott R. Functional residual capacity of normal unanesthetized Beagle dogs. Am J Vet Res 1969;30:2027–2030.
- 28↑
Bedenice D, Bar-Yishay E, Ingenito EP, et al. Head-out constant body plethysmography for measurement of specific airway resistance in conscious, sedated sheep. Am J Vet Res 2004;65:1259–1264.
- 29↑
Mundie TG, Dodd KT, Lagutchik M. Relationship of functional residual capacity to various body measurements in normal sheep. Lab Anim Sci 1992;42:589–592.
- 30
Wanner A, Reinhart ME. Respiratory mechanics in conscious sheep: response to methacholine. J Appl Physiol 1978;44:479–482.
- 31↑
Verschakelen JA, Demedts MG. Normal thoracoabdominal motions. Influence of sex, age, posture, and breath size. Am J Respir Crit Care Med 1995;151:399–405.
- 32↑
Marlin DJ, Schrotert RC, Cashman PM, et al. Movements of thoracic and abdominal compartments during ventilation at rest and during exercise. Equine Vet J Suppl 2002;34:384–390.
- 33
Mazan MR, Deveney EF, DeWitt S, et al. Energetic cost of breathing, body composition, and pulmonary function in horses with recurrent airway obstruction. J Appl Physiol 2004;97:91–97.
- 34
Gallivan GJ, McDonell WN, Forrest JB. Comparative ventilation and gas exchange in the horse and the cow. Res Vet Sci 1989;46:331–336.
- 35↑
Milsom WK. Comparative aspects of vertebrate pulmonary mechanics. In:Wood SC, ed.Comparative pulmonary physiology—current concepts. New York: Marcel Dekker, 1989;587–620.
- 36↑
Reinhold P, Close R, Lekeux P. Validation of the monofrequency forced oscillation technique for pulmonary function investigation in calves: an in vitro and in vivo study. Res Vet Sci 1992;53:165–171.
- 37↑
Gustin P, Dhem AR, Lekeux P, et al. Regulation of bronchomotor tone in conscious calves. J Vet Pharmacol Therap 1989;12:58–64.
- 38↑
Clercx C, Gustin P, Landser FJ, et al. Measurement of total respiratory impedance in dogs by the forced oscillation technique. Vet Res Commun 1993;17:227–239.
- 39↑
DiMango AM, Lopes AJ, Jansen JM. Changes in respiratory mechanics with increasing degrees of airway obstruction in COPD: detection by forced oscillometry technique. Respir Med 2006;100:399–410.
- 40
Clement J, Landser FJ, Van deWoestijne KP. Total respiratory resistance and reactance in patients with respiratory complaints with and without airways obstruction. Chest 1983;83:215–220.
- 41
Van Noord JA, Smeets J, Clement J, et al. Assessment of reversibility of airflow obstruction. Am J Respir Crit Care Med 1994;150:551–554.
- 42
Lavoie JP, Pascoe JR, Kurpershock CJ. Effects of xylazine on ventilation in horses. Am J Vet Res 1992;53:916–920.
- 43
Lavoie JP, Pascoe JR, Kurpershock CJ. Effect of head and neck position on respiratory mechanics in horses sedated with xylazine. Am J Vet Res 1992;53:1652–1657.
- 44
Broadstone RV, Gray P, Robinson NE, et al. Effects of xylazine on airway function in ponies with recurrent airway obstruction. Am J Vet Res 1992;53:1813–1817.
- 45
Lavoie JP, Phan ST, Blais D. Effects of a combination of detomidine and butorphanol on respiratory function in horses with or without chronic obstructive pulmonary disease. Am J Vet Res 1996;57:705–709.
- 46
Hsu WH, Hanson CE, Hembrough FB, et al. Effects of idazoxan, tolazoline, and yohimbine on xylazine induced respiratory changes and central nervous system depression in ewes. Am J Vet Res 1989;50:1570–1573.
- 47
Nolan A, Livingston A, Waterman A. The effects of alpha2 adrenoceptor agonists on airway pressure in anesthetized sheep. J Vet Pharmacol Ther 1986;9:157–163.
- 48↑
Klein C, Smith H-J, Reinhold P. Respiratory mechanics in conscious swine: effects of face mask, head position and bronchoconstriction evaluated by impulse oscillometery. Res Vet Sci 2003;75:71–81.
- 49↑
Fowler ME. Respiratory system. In:Fowler ME, ed.Medicine and surgery of South American camelids. 2nd ed. Ames, Iowa: Iowa State University Press, 1998;295–304.