- Department of Clinical Sciences, Cummings School of Veterinary Medicine at Tufts University, North Grafton, MA, United States
Airway hyperresponsiveness (AHR) is linked to airway inflammation and is considered a key manifestation of mild/moderate equine asthma (EA). The study purpose was to determine whether two modalities of non-invasive lung function testing (FOM—forced oscillatory mechanics vs. FP—flowmetric plethysmography) establish the same clinical diagnosis of AHR in horses, using histamine bronchoprovocation. Nineteen horses (3–25 years, 335–650 kg) with clinical signs suggestive of mild/moderate equine asthma were enrolled. FOM and FP testing was performed in each horse on two consecutive days, using a randomized cross-over design. AHR was defined by the histamine dose needed to double FOM baseline resistance, or to achieve a 35% increase in FP delta flow. Bronchoalveolar lavage fluid (BALF) was subsequently collected and stained with modified Wright's and toluidine blue stains. Binary statistical tests (related samples T-test, Mann-Whitney U, Chi-square analyses) were performed to compare study groups, with P < 0.05 considered significant. Abnormal BALF cytology confirmed EA in 14/19 (73.7%) horses. Both FOM and FP revealed AHR in 7/14 (50%) of these EA horses. An additional 4/19 (21.1%) horses showed AHR based on FP but not FOM, including two horses with normal BALF cytology. A diagnosis of AHR was more often associated with FP than FOM (P = 0.013), although the prevalence of AHR was significantly higher in EA vs. non-EA horses, regardless of testing methodology. The phase angle between thoracic and abdominal components of breathing did not differ between test groups. In conclusion, FP diagnosed AHR more frequently than did FOM, including horses with no other diagnostic evidence of EA. Without further evaluation, these two testing modalities of AHR cannot be used interchangeably.
Introduction
Airway hyperresponsiveness (AHR) is linked to airway inflammation and is considered a key manifestation of mild/moderate equine asthma or inflammatory airway disease (hereafter termed EA) in horses with chronic cough or exercise intolerance (1) as well as in humans with the similar disease, asthma. Despite this well-established link, the exact mechanism is still not fully understood (2). The diagnosis of equine asthma (1) requires, in addition to appropriate history, a documentation of lower airway inflammation based on bronchoalveolar lavage fluid (BALF) cytology or abnormal pulmonary function testing (PFT) demonstrating AHR. Due to lack of availability of the latter, BALF cytology is more commonly used in clinical practice (1).
AHR is defined as an exaggerated bronchoconstrictive response to specific or non-specific inhaled stimuli that is present only in patients with asthma and not in normal individuals (3), and, in the horse, is diagnosed most commonly by identifying increased sensitivity and responsiveness to the bronchoconstrictor agent histamine. Various methods of pulmonary function testing to document this bronchoconstrictive response are available in horses; however, invasive methods (such as esophageal balloon-pneumotachography) are more difficult to use in client-owned animals. Two non-invasive testing methods that can be used to assess AHR in client-owned horses are forced oscillatory mechanics (FOM) and flowmetric plethysmography (FP). FOM superimposes oscillations of compressed room air on spontaneous breaths in awake horses to measure total respiratory system resistance (Rrs) which reflects the degree of airway obstruction. FP combines external sensors placed on the body surface (inductance bands) measuring thoracic and abdominal excursions (plethysmographic flow), and a pneumotachograph measuring airflow at the nares. Airway obstruction is determined by subtracting peak plethysmographic flow from peak nasal flow during expiration—a variable defined as delta flow (Δflow). This variable is an established measure of airway obstruction and, similar to respiratory system resistance, increases with bronchoconstriction (4).
When compared to the conventional method (esophageal balloon-pneumotachography), a 100% increase in Rrs (PC100Rrs; PC—provocative concentration) for FOM and a 35% increase in Δflow for FP were previously found to correlate with a 65% decrease in dynamic compliance (PC65Cdyn) (5, 6). Changes in dynamic compliance (Cdyn) (esophageal balloon-pneumotachography) correlate closely with inhaled histamine concentrations (7). As such, AHR is conventionally assessed by the histamine concentration needed to reach PC65Cdyn, which lies outside the range of variability and occurs on the steepest portion of the histamine dose-response curve (8). Normal horses are not considered to reach PC65Cdyn until concentrations of >8 mg/mL histamine have been inhaled (7). Likewise, normal horses experience a 100% increase in Rrs (PC100Rrs measured by FOM) (5) or a 35% increase in Δflow (PC35Δflow measured by FP) (4) at a nebulized histamine dose >8 mg/mL (FOM) and >6 mg/mL (FP) (9), respectively. In contrast, horses with EA may respond at histamine concentrations below these cut points, consistent with AHR (a key component of EA). However, whilst several studies have used these definitions of AHR (10, 11), other authors have used different test end-points, including PC75Rrs measured by FOM and PC50Δflow measured by FP (12).
Although both FOM and FP have been separately validated against conventional esophageal balloon-pneumotachography, FOM and FP have not been directly compared, compounding difficulties in evaluating studies using either methodology to assess AHR. While each system has its advantages and limitations, FP has the benefit for clinical practice of being non-invasive and easily portable. For example, the repeatability of FP over time was determined in a large population of university horses in the field setting (4), which would have been impractical with either FOM or the conventional methodology. Our current study therefore aims to validate FP as a method of portable lung function testing by demonstrating that comparable results may be achieved using either of the available non-invasive testing modalities (FOM vs. FP) to measure AHR. We hypothesized that when compared directly, FOM and FP would identify AHR at comparable doses of nebulized histamine in horses with clinical signs of equine asthma, thus validating the use of either method to diagnose AHR in clinical cases.
Materials and Methods
Twenty-two adult horses (200–650 kg, 2–25 years old) with clinical signs suggestive of EA (poor performance, exercise intolerance, cough, or nasal discharge) were recruited from a referral population for enrolment in a prospective crossover assessment of AHR. Horses with clinical signs of infectious disease (as determined by anamnesis, pyrexia, or physical examination) or uncontrolled severe equine asthma [recurrent airway obstruction (RAO) with signs of respiratory abnormalities at rest] were excluded from the analysis. The experimental protocol was approved by the Institutional Animal Care and Use Committee (IACUC) of the Cummings School of Veterinary Medicine at Tufts University and completed following written client consent. Detailed anamnesis was obtained using a standardized questionnaire, followed by a complete physical and respiratory examination. Further diagnostic procedures were performed if indicated (thoracic radiographs, lameness evaluation, blood sampling, echocardiography) to rule out other causes of poor performance or respiratory disease.
Horses were randomly assigned to undergo histamine bronchoprovocation using both FOM and FP with a 20–30 h wash-out period between tests. The horses were not accustomed to the mask or bands for either testing method prior to the study. All horses remained hospitalized throughout the study period and were fed according to their individual management protocols. After completion of lung function testing (day 2), a bronchoalveolar lavage (BAL) was performed, as specified below.
Forced Oscillatory Mechanics (FOM)
Mono-sinusoidal, multifrequency (1–3 Hz) FOM was used to measure total respiratory system resistance (Rrs) as previously described (13, 14). In brief, horses were sedated with 0.5 mg/kg xylazine1 intravenously (IV), administered 6 mL phenylephrine2 to each nasal passage via a small catheter, and fitted with a latex-sealed low dead space facemask, via which sinusoidal flow (generated using compressed air released through a proportional pneumatic valve) was superimposed over the horses' spontaneous breathing. A pneumotachograph3 was used to measure flow at the mask opening, and the difference between mask and atmospheric pressures recorded using a pressure transducer4. For each frequency, coherence was calculated, with values >0.9 being accepted for analysis. After baseline measurements were obtained, changes in Rrs were used to monitor effects of histamine aerosol challenge (see below) (14). During all measurements, the head rested on a stand in a neutral position.
Respiratory Inductive Plethysmography (FP)
FP was performed using a commercial system5 as previously described (6). Briefly, each horse was sedated with xylazine1 (0.5 mg/kg IV), administered 6 mL phenylephrine2 to each nasal passage via a small catheter, and fitted with a low-dead space face mask, pneumotachograph, and thoracic (11th intercostal space) and abdominal (directly behind the 18th rib) inductance bands. Airflow and thoracic/abdominal volume changes were simultaneously measured during spontaneous breathing, and mathematically converted to plethysmographic flow signals. Airway obstruction (delta flow, Δflow), thoraco-abdominal synchronicity (phase angle or Φ theta) and ventilatory parameters were recorded. Delta flow is the value recorded by the proprietary software by subtracting the flow signal generated by the thoracic and abdominal volume change at the body surface (measured by inductance bands) from the pneumotachograph flow at peak expiration (4, 6).
Thoraco-abdominal synchronicity, or phase angle (theta), is calculated by plotting the abdominal and thoracic inductance bands on x-y coordinates and using Lissajous figures as described previously (15, 16). Normal synchronous breathing has a phase angle close to 0°, which increases with increasing airway obstruction and resulting asynchrony, eventually reaching 180° which indicates complete paradoxical motion of the rib cage and abdomen. Respiratory rate was calculated breath by breath and recorded by the software, as were tidal volume, peak inspiratory and expiratory flows, minute ventilation, and inspiratory and expiratory times. During all measurements, the head rested on a stand in a neutral position.
Histamine Bronchoprovocation (for Both FOM and FP)
Histamine bronchoprovocation was performed as previously described (14). Initial baseline lung function parameters were recorded via FOM or FP (Rrs and Δflow, respectively). Nebulization of 2 mL of 0.9% saline (negative control or saline baseline) over 2 min using a low dead-space face mask with a portable air compressor6 and nebulizer7 was performed and repeated for subsequent increasing concentrations of histamine diphosphate8 in saline solution (2, 4, 8, and 16 mg/mL). Lung function measurements were repeated following each nebulized dose until either a 100% increase in Rrs from saline baseline (FOM PC100Rrs) or >50% increase in Δflow (FP PC35Δflow) was achieved, or the horse displayed a clinical reaction to histamine administration such as a notably increased respiratory rate or effort, or repeated coughing.
Following FOM testing, a dose-response curve was generated for Rrs at 1 Hz frequency to determine the histamine dose required to reach a doubling in saline baseline resistance (PC100Rrs) for each horse. Similarly, FP data were further evaluated using a breath-by-breath analysis, to discard signals that were obviously erroneous (e.g., first breath of recording was abnormally short, tidal volume non-physiologically small or large). The delta flow for each histamine concentration was calculated by the software by discarding the two highest values and then averaging the next five highest values. Respiratory rate and theta changes were generated manually in the same way. This enabled measurement of the peak response to each histamine dose, whilst still excluding outliers. The PC35Δflow (35% increase in delta flow from saline) was then recalculated from a manually-derived dose-response curve. Dose-response curves for respiratory rate and theta were similarly generated.
Using the dose-response curves generated, AHR was defined when the PC100Rrs was reached at ≤8 mg/mL for FOM, or the PC35Δflow was reached at ≤6 mg/mL for FP.
Bronchoalveolar Lavage and Evaluation of Bronchoalveolar Lavage Fluid
After all PFTs were completed, horses were sedated with detomidine9 (0.01 mg/kg) and butorphanol10 (0.01 mg/kg) IV and bronchoalveolar lavage performed using a commercial cuffed BAL tube11 or endoscope12 as previously described (17). A total of 120 mL of 0.3% lidocaine13 was used to anesthetize the upper respiratory tract and bronchi during passage of the tube or endoscope. Once wedged in a bronchus, two aliquots of 250 mL saline14 were instilled and then drawn back by suction at 10-cm H2O pressure15. The two samples were pooled and 10 mL added to a red top glass tube which was refrigerated and submitted within 4 h of collection for cytology. Slides were prepared after cytocentrifugation and stained with Wright-Giemsa and toluidine blue stains. Cells (n = 500) were classified by one of the authors (MRM) as the percentage of total cells that were macrophages, lymphocytes, neutrophils, eosinophils, and mast cells (400x magnification). Horses were classified as having an abnormal BALF if there were >5% neutrophils, >2% mast cells, or >1% eosinophils, and these horses were diagnosed with EA (1).
Statistical Analysis
Sample size calculations estimated that 16 horses would provide sufficient statistical power (alpha = 0.05 and power >85%) to reject the null hypothesis of test equivalence (FP vs. FOM) if the alternative is true. All continuous data were presented descriptively as mean ± SD or median ± range, and compared between test groups using related samples T-test or Mann-Whitney U-test, based on the normality of data distribution (Kolmogorov-Smirnov analysis). Chi squared (χ2) analyses were used to compare a diagnosis of AHR between lung function testing modalities. The coefficient of variation was calculated for initial baseline measurements (over 3 min) for each horse for the delta flow, theta, and respiratory rate, and the means compared using Mann-Whitney U-test. Statistical analyses were performed using statistical software16, with P < 0.05 considered significant.
Results
A total of 19 horses met the inclusion criteria and completed the study. Three additional horses were initially enrolled, but excluded due to evidence of either infectious respiratory disease or RAO at the time of examination. These were five mares and 14 geldings, with a mean age of 13 ± 5.5 years, weighing 507 ± 80 kg, in good to obese body condition (median BCS of 6, range 4–9; not recorded in one horse). Exercise level varied, with one racehorse, nine horses in a moderate level of exercise (eventers, showjumpers, etc.) and the remaining nine horses in light work (pleasure riding, etc.). Baseline characteristics are presented in Table 1. All horses were clinically normal at rest.
Table 1. Type and duration of clinical signs in horses presented for the evaluation of Equine Asthma (EA).
Equine asthma was diagnosed in 14 horses based on combination of history of exercise intolerance or poor performance along with abnormal BALF cytology, whereas five horses were included in the study with a diagnosis of upper airway dysfunction, obesity (leading to exercise intolerance), exercise induced pulmonary hemorrhage (EIPH), or unknown disease. Of those horses with a diagnosis of equine asthma based on clinical history/examination and BALF cytology, 9/14 (64%) had >5% neutrophils, 9/14 (64%) had >2% mast cells, and 4/14 (21%) had >1% eosinophils. An increase in more than one cell type was seen in 7/14 (50%) study animals, where a high percentage of eosinophils in BALF was always associated with a concurrent elevation in either neutrophils or mast cells.
Baseline respiratory system resistance (FOM) and delta flow (FP) are presented in Table 2. There was no significant difference in resting pulmonary function parameters between horses diagnosed with EA and those with normal BALF cytology. However, FOM testing established a diagnosis of airway hyperresponsiveness (100% increase in Rrs from baseline (PC100Rrs) or a clinical reaction at ≤8 mg/mL histamine) in seven horses, all of which had abnormal BALF cytology. Employing a cut point of PC75Rrs instead of PC100Rrs did not alter the diagnosis of AHR in any horse. The mean dose of histamine required to reach PC100Rrs was lower in horses with BALF cytology compatible with equine asthma compared to those with normal BALF cytology but the difference was not significant (Table 2). FP diagnosed AHR based on PC35Δflow ≤6 mg/mL histamine in 11 horses, of which 9/11 (82%) showed an abnormal BALF cytology, and 2/11 (18%) had a normal BALF cytology. As with FOM, the mean histamine dose required to reach PC35Δflow was lower in horses with BALF compatible with EA compared to those with normal BALF cytology, but the difference was not significant (Table 2). Overall, FP was more likely to diagnose AHR in horses with both normal and abnormal BALF compared to FOM (p = 0.013). The individual cellular profile of BALF was not statistically associated with a diagnosis of AHR based on FOM or FP (Table 3).
Table 2. Pulmonary function test results (median and range) for saline baseline and histamine provocation using FOM and FP.
Table 3. AHR as measured by FOM (PC100Rrs) and FP (PC35Δflow) compared to BALF cytology subtypes [median (range)].
The FP-measured parameter theta (phase angle) was not useful in distinguishing between horses with EA and those without, either during baseline measurements (Table 4) or using histamine provocation. Although theta did increase following histamine administration (saline baseline theta 18.88 ± 8.04; theta after final concentration of histamine 28 ± 14.53; p = 0.006), a significant association with percentage increase in theta and EA diagnosis was not present (using theta PC50 as a cut off for diagnosis of EA p = 0.500). In contrast, initial baseline respiratory rate was significantly different between EA and non-EA horses (p = 0.001) with a corresponding decrease in inspiratory and expiratory times in horses with EA (Table 4). Tidal volume was also lower in horses diagnosed with EA (p = 0.039). Respiratory rate also increased with histamine challenge as expected but the difference between groups did not reach statistical significance.
Table 4. FP derived parameters during initial baseline measurements, showing mean and SD for EA and non-EA horses and the p-value comparing the groups.
Both lung function methods were well-tolerated by all horses. Horses received 0.49 mg/kg (0.33–0.87 mg/kg) xylazine for FOM, and 0.50 mg/kg (0.33–0.99 mg/kg) for FP [median (range)]. Fourteen horses received the same total dose of sedation to perform both testing modalities. The remaining five horses required additional sedation to complete the tests, with one horse requiring a higher total dose of sedation to perform FOM, and the other four horses required a higher total dose of sedation for FP.
Discussion
This study demonstrates that FOM and FP cannot be used interchangeably to assess AHR in horses as test results are discordant, with AHR diagnosed in 11 horses using FP and in only seven of these horses using FOM. Moreover, FP was more likely to diagnose AHR even in horses without a clinical diagnosis of EA. Although both testing methods have previously been compared to conventional esophageal balloon-pneumotachography (5, 6), the latter testing method was not included in this study of client owned horses due to its degree of invasiveness and lack of use outside of research purposes (18).
The diagnosis of EA in this study was based on a clinical history of poor performance, exercise intolerance, or chronic cough along with abnormal BALF cytology rather than the documentation of AHR. While reliance on clinical history and BALF cytology for the diagnosis of EA is standard practice (1), it has been shown that a subset of horses with clinical signs suggestive of equine asthma but normal BALF cytology, exhibit pulmonary dysfunction based on histamine bronchoprovocation testing (9, 19). Moreover, the documented lack of improvement in BALF-based inflammation in multiple studies suggests that response to treatment should be assessed by PFT and not by re-evaluation of BALF cytology (20–24). These observations support the importance of pulmonary function testing with bronchoprovocation in the management of horses with equine asthma.
Airway hyperresponsiveness was a frequent finding in our study in horses with a complaint of cough or poor performance, with 11/19 [58%] affected, comparable to previous reports [24/45 (53%)] (9). However, a higher percentage of horses showed AHR in the absence of BALF abnormalities in the latter report [11/45 (24%)], compared to our current study [2/19 (11%)]. Various studies have associated increased AHR with increased mast cell concentrations in BALF (14, 25, 26), increased eosinophils (27) or both (28); while others saw no association between cell type and AHR (9, 10). Using open plethysmography in horses without any clinical history or evidence of EA, Cullimore et al. (29) documented AHR in 52% and BALF cytology compatible with EA in 42% of horses. However, AHR and BALF cytology were not correlated, and indeed, 53% of horses with normal BALF cytology demonstrated AHR. The results of our current study do not indicate a relationship between abnormal BALF cytology and AHR, irrespective of the cell type evaluated or testing method used.
The expectation that airway hyperresponsiveness must necessarily be conjoined with airway inflammation has been challenged in humans for some time (30, 31). In fact, recent studies in horses (9, 29) have supported that AHR may not be associated with airway inflammation. Certainly, airway hyperresponsiveness has been observed in non-asthmatic, atopic humans as well as in horses (32). This further underscores the importance of identifying accurate PFT methods for the diagnosis of EA, where both inflammatory and functional disturbances may independently contribute to the disease pathogenesis.
The reason for the difference between the two testing methods (FOM and FP) is not immediately obvious, and the human literature provides ample evidence that varying methods of lung function testing are not necessarily comparable. In people, the forced oscillatory technique was shown to be both more sensitive and more specific than spirometry [forced expiratory volume in one second (FEV1)] in assessing bronchial hyperresponsiveness (33). Likewise, body plethysmography, which is similar but not the same as flowmetric plethysmography (34), demonstrated evidence of AHR in a high number of clinically diagnosed asthmatic individuals who did not reach the diagnostic threshold required by spirometry. A close relationship between the dose-response curves obtained by plethysmography and a forced oscillatory technique was identified in human asthmatics in one study (35). Another study reported higher test sensitivity of body plethysmography compared to impulse oscillometry when investigating both methacholine and allergen challenges in people (36).
It is possible that the definition of AHR within each test requires adjustment of cutpoints (e.g., decreasing the sensitivity of FP by increasing the cutpoint to a 50% increase rather than a 35% increase in delta flow, PC35Δflow vs. PC50Δflow), or by similarly increasing the sensitivity of FOM. Looking first at FOM, the results of the current study would not have changed, and the same diagnosis would be made in each case regardless of whether an end point of PC75Rrs or PC100Rrs at ≤6 or ≤8 mg/mL was used. With respect to FP, a cut off of PC35Δflow for FP at ≤6 mg/mL histamine was used as this was found to correlate with conventional testing (PC65Cdyn) (6), resulting in 11 horses being diagnosed with AHR. However, had we used a cutpoint of ≤8 mg/mL histamine, three horses with EA based on history/clinical signs and BALF cytology but without AHR at PC35Δflow ≤6 mg/mL histamine, would have been placed in the AHR category. Alternatively, if PC50Δflow had been used (12) at ≤6 mg/mL histamine then one horse with AHR at PC35Δflow ≤6 mg/mL histamine would have been re-classified as normal. Altering the cut point for FP would not have changed the results for the two horses with AHR but normal BALF cytology. Regardless of alternative cutpoints, the lack of correlation between FOM and FP would not have changed.
In the four horses where the diagnosis of AHR differed between FOM and FP, only two were diagnosed with EA based on clinical assessment and BALF cytology. None of the four horses reacted with FOM, all having a PC100Rrs of >8 mg/mL histamine. Thus, FOM appears to be more specific in the diagnosis of AHR, whereas FP appears to be more sensitive.
Of the two horses not clinically diagnosed with EA, one was extremely over-conditioned (BCS 9/9) and solely presented with a clinical complaint of poor performance. In this case, obesity was determined to be the most likely cause. The second horse may have experienced post-viral AHR without cellular evidence of airway inflammation. The results of pulmonary function testing with FOM vs. PF were discordant in these two horses, with FOM demonstrating a normal baseline respiratory resistance and normal response to histamine (PC100Rrs = 16 and >8 mg/mL), and FP demonstrating marked hyperresponsiveness of the airways (PC35Δflow = 2.43 and 2.89 mg/mL, respectively). It might be expected that the obese horse would have abnormal pulmonary function, as obesity alters lung and chest properties in humans. The latter is based on both mechanical factors increasing airway closure and gas trapping, and inflammatory factors (such as leptin) increasing airway responsiveness (37–41). This suggests that an obese horse may be more likely to display AHR in the absence of equine asthma, but does not explain the disparity between FP and FOM for the obese horse in this study. Similarly, diverse methods such as body plethysmography and oscillometry yield mixed results in obese children, with an inconsistent finding of increased resistance (42–45). If obesity in this horse led to increased gas trapping, then this would increase the delta flow due to increased gas compression. The principle of plethysmography relies on the intrathoracic organs being incompressible (46), yet an increase in adipose tissues may alter this. An inability of the diaphragm to move caudally due to abdominal fat may also lead to increasing signal from the abdominal contribution to plethysmography flow as the horse attempts to breathe deeper in response to histamine. Altered body shape could also lead to increased likelihood of the bands displacing from their location during calibration. The latter is a significant concern in humans, leading to the sole use of a thoracic band in obese patients (47, 48). Similar to people (47), further evaluation of FP in obese horses is required in order ascertain its reliability.
The second horse without a clinical diagnosis of EA was the oldest horse of the study (21 years). It presented for poor performance and a short duration cough, showing an increased number of lymphocytes on BALF with no other abnormalities. Whilst increased lymphocytes may be an age related change, the degree of increase (72%) is still abnormal (10). The FP test showed a sudden increase in delta flow and theta at 4 mg/mL, but the respiratory rate and tidal volumes were unchanged (6–5 breaths per minute, 8.26–8.47 L respectively). The horse spontaneously returned to normal performance following evaluation with no therapy, and was thought to have suffered from a respiratory virus with prolonged recovery. While exaggerated airway response to non-specific stimuli is well-documented in humans (49), and IAD (mild/moderate equine asthma) has been previously associated with recent exposure to respiratory viruses in horses (12), the reason for the lack of concordance between the two diagnostic tests remains unclear.
The two horses with a clinical diagnosis of equine asthma but with a negative bronchoprovocation test using FOM had only subtle changes on BALF, with one having 6% neutrophils and the other having 5% mast cells. Examining the horse's dose response curve following FP may help explain the test disparity in the first horse (6% neutrophils), which showed a sudden increase in delta flow, theta and respiratory rate after nebulization with 8 mg/mL histamine, at the very start of the 2 min recording period. It is possible that this initial and transient peak was missed when testing with FOM due to the inherently greater time taken to initiate testing immediately following histamine nebulization. This horse, moreover, was one which would have been classified as normal on histamine bronchoprovocation if a cutpoint of PC50Δflow at ≤6 mg/mL of histamine had been used. All other horses with AHR in this study showed a consistent increase in delta flow throughout the recording period. This suggests that an advantage of FP over FOM is the ability to monitor the horse's response to histamine over a more prolonged time period, whereas a short transient peak reaction may be missed when testing with FOM. The final horse (with 5% mast cells) had a gradual increase in delta flow which was consistent throughout each recording period, but showed no change in theta or respiratory rate (respiratory rate 7 breaths/min, tidal volume minimally increased from 4.36 to 5.01 L). A definite reason for the difference in results between testing modalities of AHR could not be identified for this horse.
FOM and FP are testing methods that evaluate different aspects of pulmonary function, which may contribute to the differences in diagnosis of AHR between the two tests. FOM measures changes in oscillatory flow and pressure and, at lower frequencies (i.e., 1–3 Hz), computes respiratory system resistance at the level of the lower airways (18). In contrast, FP computes the difference between plethysmographic and pneumotachograph flows (Δflow), representing gas compression (expiration) and expansion (inspiration), as a measure of airway obstruction (18). FOM, as a non-invasive test of pulmonary function, has been previously validated against more invasive testing techniques (esophageal balloon-pneumotachography), and its use with histamine bronchoprovocation is well-described (5, 13, 25). Although the latter technique is routinely used in our hospital to diagnose EA in horses, FOM testing has specific limitations. Achieving good quality data and a high signal-to-noise ratio at 1 Hz oscillatory frequency (a wavelength that bests reflects lower airway function in large animals) can be challenging, due to interference from the horse's breathing (13). Respiratory interference increases at lower Hz frequencies that approach the horse's spontaneous breathing rate (18). The latter is quantified as coherence based on the measured signal-to-noise ratio, where Rrs with a coherence of <0.9 is rejected and eliminated from the analysis (14). This makes the methodology less useful for horses with an abnormal respiratory rate or effort at rest (e.g., RAO or upper respiratory abnormalities), as turbulent airflow and breathing rates approaching the 1 Hz frequency can result in data rejection and incomplete testing. Due to measuring Rrs of the entire respiratory system, including upper respiratory tract and the chest wall, abnormalities of the upper airway (such as epiglottic entrapment or laryngeal paresis) may lead to an increased Rrs due to increased turbulence in the absence of lower airway disease. Upper airway endoscopy should be used to identify or rule out such conditions. Inclusion of the chest wall in the measurement of Rrs is considered a limitation as it decreases the sensitivity (13), however it contributes <10% of the total value and remains constant (18), and as such is likely to have minimal effect on bronchoprovocation. Whilst the requirement for obtaining an appropriate coherence value of ≥0.9 can limit the use of FOM in some patients, it also enables elimination of abnormal and inaccurate breaths (such as sniffing or coughing) from being recorded, thus improving the reliability of the data. FOM is also limited in that it is not portable due to the requirement of pressurized gas, and is a system that is only available in a few centers, as it requires expertise for its use.
Despite its limitations, FOM remains a reliable method of lung function testing in horses that are clinically normal at rest, providing diagnostically useful information in both baseline measurements and following histamine bronchoprovocation. In this study, FOM diagnosed AHR only in horses with abnormal BALF cytology, suggesting reliability of positive results, and providing a clinically useful and non-invasive method of assessing response to treatment in these horses.
Whilst FOM is a useful and reliable testing method, FP has the benefit of portability and ease of testing under field conditions. FP also offers fully automated software requiring minimal training for its use: the latter is also a limitation. Real-time evaluation of breath-by-breath data during testing is not possible, although this is required due to the large number of erroneously recorded breaths [which would be eliminated from analysis using FOM due to poor coherence and from conventional testing by use of “filters” or inclusion criteria (18)]. These breaths often have higher delta flow values and contribute to the software produced PC value. The software also automatically stops the test once PC35Δflow is reached, which could result in erroneously completing the test at an inappropriately low concentration of histamine once final analysis of the data is completed. These concerns could be easily eliminated by alterations in the software. FP is based on the concept of whole body plethysmography, and replacing a box with inductive bands can lead to creation of artifacts. Before commencing testing the bands are calibrated to the pneumotachograph based on inspiratory volume, which may cause an overcorrection resulting in pre-existing abnormalities being less well-defined (6). However, further progression in response to histamine is still detected. This also means that the test may be more reliable for defining lower airway disease in horses with a structural upper airway abnormality, and across horses of different sizes as differences in tidal volume are eliminated, due to the horse acting as its own control. One of the limitations of the bands is movement or slipping following calibration, which is more likely to occur with the abdominal band due to body shape. Caudal movement of this band to a slightly narrower region would result in it being less stretched and therefore having greater movement during breathing and falsely increasing the delta flow. Alternatively, the band could move to an area with increased movement, leading to the same result, or to an area with less movement, falsely decreasing the delta flow. This artifact is reduced by careful placement and monitoring of the bands during testing, with re-calibration if movement of the band is suspected. Whilst the bands are calibrated to the individual horse accounting for differences in tidal volume, and tachypnea alone does not increase delta flow, it has been shown that hyperpnea increases delta flow due to accentuation of intrapulmonary gas compression (6). Progressive hyperpnea could occur during testing, because of mask placement (18), use of xylazine (50), or as a response to histamine nebulization. Thus, the occurrence of hyperpnea during testing could lead to a falsely low PC35Δflow with diagnosis of airway hyperresponsiveness. Monitoring of tidal volume during testing would thus be useful. Whilst FP has its limitations when used as a fully automated system, it provides useful information if the data is critically evaluated. Experience with its use eliminates some of the inaccuracies and artifacts that can occur.
One of the significant benefits of FP is the information provided alongside delta flow. The slightly increased minute ventilation, respiratory rate, and peak inspiratory flow, alongside a slight decrease in expiratory time observed in horses with mild/moderate equine asthma in this study is consistent with findings in horses with severe equine asthma (RAO) (15). Although the subtle differences between horses with equine asthma and normal horses may not allow for diagnosis based on this alone, it provides valuable additional information for both clinical evaluation and research.
It has been suggested that theta (phase angle) can be used to measure more subtle changes in thoracic and abdominal contributions to breathing pattern for the assessment of AHR. However, a previous study in foals (16) found no correlation between phase angle (theta) and other measures of lower airway constriction. In contrast, phase angle was significantly increased in horses with heaves (15) compared to normal controls. The latter study also identified that histamine bronchoprovocation led to variable alterations in phase angle in normal horses. However, in the present study the high variability within patients during baseline recording and following histamine administration, suggests that these values may not be accurate using the algorithm currently employed by the software. This is likely because the more the respiratory pattern deviates from sinusoidal, the less accurate Lissajous figure analysis is in estimating phase difference in thoracoabdominal synchrony. Examination of the respiratory pattern in all horses in this study showed a significant deviation from sinusoidal motion which would therefore render inaccurate the embedded Lissajous figure analysis (51). It is therefore not surprising that results of our current data and previous reports suggest that phase angle (theta) cannot be used to reliably assess AHR in response to histamine bronchoprovocation in horses with clinical signs suggestive of EA.
Limitations of this study include small sample size, potential confounding due to sedation, and potential variability of histamine dose. Sedation with xylazine, required for both forms of PFT used in this study, causes bronchodilation (52) and does affect lung mechanics in horses. Xylazine, however, does not block the bronchoconstrictive effects of histamine. This effect was mitigated as much as possible by standardizing with the same starting dose for each testing method (14). Also, following sedation each horse's head position was maintained in a neutral position (53), which may actually be easier to consistently achieve in the sedated horse (14). Using FOM, Rrs at 1 Hz frequency is stable for up to 45 min following xylazine administration at a similar dose used in this study (54). Thus, the effect on testing with FOM was likely minimal. However, other changes in respiratory parameters such as a decrease in respiratory rate and minute ventilation with an increase in tidal volume (50) may have had a greater effect on FP. Whilst a comparison of FOM in sedated and non-sedated horses has been performed (55), all studies using FP have used sedation in order to perform testing. The direct effect of xylazine on delta flow has therefore not been evaluated.
Histamine was nebulized over 2 min, as previous described for both FOM and FP (4–7, 9, 10, 14–16, 26). However, this may lead to variability of dose for individuals, or to different histamine doses for the same individual on consecutive days, since variations of tidal volume or respiratory rate may change the total amount of histamine inhaled (7). Although this is a limitation, it is likely that each horse responded similarly to nebulization on consecutive days, thus comparison of testing in each horse is still reliable.
Conclusion
Our study indicates that histamine bronchoprovocation using FP establishes a diagnosis of AHR more frequently than does testing with FOM, but that FOM is more likely to accord with a cytological diagnosis of equine asthma. As such, the two pulmonary function testing methods cannot be used interchangeably without further adjustments in current test algorithms or assessment of the effect of xylazine and hyperpnea on delta flow. Further evaluation of the use of FP in specific subgroups (such as obesity) may also be indicated.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by Clinical Sciences Research Committee (CSRC) at Cummings School of Veterinary Medicine. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
CD, DB, and MM were involved in the experimental design, patient recruitment, data collection, analyses, and manuscript writing. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank the Companion Animal Health Fund (CAHF) for financial support of this study (short-term training grant).
Footnotes
1. ^Xylazine 100 mg/mL, AnaSed LA, VetOne, Boise, ID.
2. ^Phenylephrine HCl 1%, Nose Drops, Cardinal Health, Dublin, OH.
3. ^Fleisch, No 4, OEM Medical, Lenoir, NC.
4. ^DP45-28, Validyne Engineering, Northridge, CA.
5. ^Open Pleth 4.6.7.0, Ambulatory Monitoring Inc.
6. ^ProNeb Turbo, Model 38B0201, Pari, Midlothian, VA.
7. ^LC Plus, Pari.
8. ^Histamine bisphosphate monohydrate, Sigma-Aldrich, St. Louis, MO.
9. ^Detomidine 10 mg/mL, Dormosedan, Zoetis Inc., Kalamazoo, MI.
10. ^Butorphanol 10 mg/mL, Torbugesic, Zoetis Inc., Kalamazoo, MI.
11. ^Bivona Medical Technologies, Gary, IL.
12. ^3 m length 12.9 outer diameter Endoscope, Fujinon, Wayne, NJ.
13. ^Lidocaine HCl 2%, Hospira Inc., Lake Forest, IL.
14. ^Normal saline solution, VEDCO Inc., St. Joseph, MO.
15. ^EasyVac, Precision Medical Inc., Northampton, PA.
16. ^SAS, SPSS.
References
1. Couetil LL, Cardwell JM, Gerber V, Lavoie JP, Leguillette R, Richard EA. Inflammatory airway disease of horses–revised consensus statement. J Vet Int Med. (2016) 30:503–15. doi: 10.1111/jvim.13824
2. Audrit KJ, Delventhal L, Aydin O, Nassenstein C. The nervous system of airways and its remodeling in inflammatory lung diseases. Cell Tissue Res. (2017) 367:571–90. doi: 10.1007/s00441-016-2559-7
3. Nensa F, Kotschy-Lang N, Smith H-J, Marek W, Merget R, editors. Assessment of Airway Hyperresponsiveness: Comparison of Spirometry and Body Plethysmography. Dordrecht: Springer Netherlands (2013).
4. Nolen-Walston RD, Kuehn H, Boston RC, Mazan MR, Wilkins PA, Bruns S, et al. Reproducibility of airway responsiveness in horses using flowmetric plethysmography and histamine bronchoprovocation. J Vet Int Med. (2009) 23:631–5. doi: 10.1111/j.1939-1676.2009.0307.x
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–80.
6. Hoffman A, Kuehn H, Riedelberger K, Kupcinskas R, Miskovic MB. Flowmetric comparison of respiratory inductance plethysmography and pneumotachography in horses. J Appl Physiol. (2001) 91:2767–75. doi: 10.1152/jappl.2001.91.6.2767
7. Klein HJ, Deegen E. Histamine inhalation provocation test: method to identify nonspecific airway reactivity in equids. Am J Vet Res. (1986) 47:1796–800.
8. Derksen FJ, Scott D, Robinson NE, Slocombe RF, Armstrong PJ. Intravenous histamine administration in ponies with recurrent airway obstruction (heaves). Am J Vet Res. (1985) 46:774–7.
9. Wichtel M, Gomez D, Burton S, Wichtel J, Hoffman A. Relationships between equine airway reactivity measured by flowmetric plethysmography and specific indicators of airway inflammation in horses with suspected inflammatory airway disease. Equine Vet J. (2016) 48:466–71. doi: 10.1111/evj.12482
10. Pacheco AP, Paradis MR, Hoffman AM, Hermida P, Sanchez A, Nadeau JA, et al. Age effects on blood gas, spirometry, airway reactivity, and bronchoalveolar lavage fluid cytology in clinically healthy horses. J Vet Int Med. (2014) 28:603–8. doi: 10.1111/jvim.12318
11. Hoffman A, Mazan MR, editors. How to test for airway hyper-reactivity: toward a better diagnosis of inflammatory airway disease. In: Annual Convention of the AAEP. Seattle, WA: IVIS (2005).
12. Houtsma A, Bedenice D, Pusterla N, Pugliese B, Mapes S, Hoffman AM, et al. Association between inflammatory airway disease of horses and exposure to respiratory viruses: a case control study. Multidiscip Respir Med. (2015) 10:33. doi: 10.1186/s40248-015-0030-3
13. Young SS, Tesarowski D. Respiratory mechanics of horses measured by conventional and forced oscillation techniques. J Appl Physiol. (1994) 76:2467–72. doi: 10.1152/jappl.1994.76.6.2467
14. Hoffman AM, Mazan MR, Ellenberg S. Association between bronchoalveolar lavage cytologic features and airway reactivity in horses with a history of exercise intolerance. Am J Vet Res. (1998) 59:176–81.
15. Hoffman AM, Oura TJ, Riedelberger KJ, Mazan MR. Plethysmographic comparison of breathing pattern in heaves (recurrent airway obstruction) versus experimental bronchoconstriction or hyperpnea in horses. J Vet Int Med. (2007) 21:184–92. doi: 10.1111/j.1939-1676.2007.tb02945.x
16. Miller C, Hoffman AM, Hunter J. Thoracoabdominal asynchrony failed to grade airway obstructions in foals. J Appl Physiol. (2000) 88:2081–7. doi: 10.1152/jappl.2000.88.6.2081
17. Mazan M. Respiratory secretions. In: Pusterla N, Higgins J, editors. Interpretation of Equine Laboratory Diagnostics. Hoboken, NJ: John Wiley & Sons (2017). p. 363.
18. Hoffman AM. Clinical Application of Pulmonary Function Testing in Horses. Ithaca, NY: IVIS (2002).
19. Couetil LL, Hoffman AM, Hodgson J, Buechner-Maxwell V, Viel L, Wood JL, et al. Inflammatory airway disease of horses. J Vet Int Med. (2007) 21:356–61. doi: 10.1111/j.1939-1676.2007.tb02975.x
20. Depecker M, Richard EA, Pitel PH, Fortier G, Leleu C, Courouce-Malblanc A. Bronchoalveolar lavage fluid in Standardbred racehorses: influence of unilateral/bilateral profiles and cut-off values on lower airway disease diagnosis. Vet J. (2014) 199:150–6. doi: 10.1016/j.tvjl.2013.10.013
21. Leguillette R, Tohver T, Bond SL, Nicol JA, McDonald KJ. Effect of dexamethasone and fluticasone on airway hyperresponsiveness in horses with inflammatory airway disease. J Vet Int Med. (2017) 31:1193–201. doi: 10.1111/jvim.14740
22. Lavoie JP, Leguillette R, Pasloske K, Charette L, Sawyer N, Guay D, et al. Comparison of effects of dexamethasone and the leukotriene D4 receptor antagonist L-708,738 on lung function and airway cytologic findings in horses with recurrent airway obstruction. Am J Vet Res. (2002) 63:579–85. doi: 10.2460/ajvr.2002.63.579
23. Gerber V, Schott Ii HC, Robinson NE. Owner assessment in judging the efficacy of airway disease treatment. Equine Vet J. (2011) 43:153–8. doi: 10.1111/j.2042-3306.2010.00156.x
24. Couetil LL, Art T, de Moffarts B, Becker M, Melotte D, Jaspar F, et al. Effect of beclomethasone dipropionate and dexamethasone isonicotinate on lung function, bronchoalveolar lavage fluid cytology, and transcription factor expression in airways of horses with recurrent airway obstruction. J Vet Int Med. (2006) 20:399–406. doi: 10.1111/j.1939-1676.2006.tb02875.x
25. Bedenice D, Mazan MR, Hoffman AM. Association between cough and cytology of bronchoalveolar lavage fluid and pulmonary function in horses diagnosed with inflammatory airway disease. J Vet Int Med. (2008) 22:1022–8. doi: 10.1111/j.1939-1676.2008.0109.x
26. Mazan MR, Vin R, Hoffman AM. Radiographic scoring lacks predictive value in inflammatory airway disease. Equine Vet J. (2005) 37:541–5. doi: 10.2746/042516405775314899
27. Hare JE, Viel L. Pulmonary eosinophilia associated with increased airway responsiveness in young racing horses. J Vet Int Med. (1998) 12:163–70. doi: 10.1111/j.1939-1676.1998.tb02112.x
28. Richard EA, Fortier GD, Denoix JM, Art T, Lekeux PM, Van Erck E. Influence of subclinical inflammatory airway disease on equine respiratory function evaluated by impulse oscillometry. Equine Vet J. (2009) 41:384–9. doi: 10.2746/042516409X366121
29. Cullimore AM, Secombe CJ, Lester GD, Robertson ID. Bronchoalveolar lavage fluid cytology and airway hyper-reactivity in clinically normal horses. Austral Vet J. (2018) 96:291–6. doi: 10.1111/avj.12721
30. Brusasco V, Crimi E, Pellegrino R. Airway hyperresponsiveness in asthma: not just a matter of airway inflammation. Thorax. (1998) 53:992–8. doi: 10.1136/thx.53.11.992
31. Baroffio M, Barisione G, Crimi E, Brusasco V. Noninflammatory mechanisms of airway hyper-responsiveness in bronchial asthma: an overview. Ther Adv Respir Dis. (2009) 3:163–74. doi: 10.1177/1753465809343595
32. Djukanovic R, Lai CK, Wilson JW, Britten KM, Wilson SJ, Roche WR, et al. Bronchial mucosal manifestations of atopy: a comparison of markers of inflammation between atopic asthmatics, atopic nonasthmatics and healthy controls. Eur Respir J. (1992) 5:538–44.
33. Schmekel B, Smith HJ. The diagnostic capacity of forced oscillation and forced expiration techniques in identifying asthma by isocapnic hyperpnoea of cold air. Eur Respir J. (1997) 10:2243–9. doi: 10.1183/09031936.97.10102243
34. Khalid I, Morris ZQ, Khalid TJ, Nisar A, DiGiovine B. Using spirometry to rule out restriction in patients with concomitant low forced vital capacity and obstructive pattern. Open Respir Med J. (2011) 5:44. doi: 10.2174/1874306401105010044
35. Chinet T, Pelle G, Macquin-Mavier I, Lorino H, Harf A. Comparison of the dose-response curves obtained by forced oscillation and plethysmography during carbachol inhalation. Eur Respir J. (1988) 1:600–5.
36. Naji N, Keung E, Kane J, Watson RM, Killian KJ, Gauvreau GM. Comparison of changes in lung function measured by plethymography and IOS after bronchoprovocation. Respir Med. (2013) 107:503–10. doi: 10.1016/j.rmed.2012.12.022
37. Peters U, Suratt BT, Bates JHT, Dixon AE. Beyond BMI: obesity and lung disease. Chest. (2018) 153:702–9. doi: 10.1016/j.chest.2017.07.010
38. Jones RL, Nzekwu MM. The effects of body mass index on lung volumes. Chest. (2006) 130:827–33. doi: 10.1378/chest.130.3.827
39. Zerah F, Harf A, Perlemuter L, Lorino H, Lorino AM, Atlan G. Effects of obesity on respiratory resistance. Chest. (1993) 103:1470–6. doi: 10.1378/chest.103.5.1470
40. Salome CM, King GG, Berend N. Physiology of obesity and effects on lung function. J Appl Physiol. (2010) 108:206–11. doi: 10.1152/japplphysiol.00694.2009
41. Chapman DG, Berend N, King GG, Salome CM. Increased airway closure is a determinant of airway hyperresponsiveness. Eur Respir J. (2008) 32:1563–9. doi: 10.1183/09031936.00114007
42. Parraguez Arevalo A, Rojas Navarro F, Ruz Cespedes M, Medina Gonzalez P, Escobar Cabello M, Munoz Cofre R. The impact of obesity on specific airway resistance and conductance among schoolchildren. Arch Argentinos Pediatria. (2018) 116:e227–33. doi: 10.5546/aap.2018.eng.e227
43. Kalhoff H, Breidenbach R, Smith HJ, Marek W. Impulse oscillometry in preschool children and association with body mass index. Respirology. (2011) 16:174–9. doi: 10.1111/j.1440-1843.2010.01906.x
44. Lauhkonen E, Koponen P, Nuolivirta K, Paassilta M, Toikka J, Saari A, et al. Obesity and bronchial obstruction in impulse oscillometry at age 5-7 years in a prospective post-bronchiolitis cohort. Pediatric Pulmonol. (2015) 50:908–14. doi: 10.1002/ppul.23085
45. Assumpcao MS, Ribeiro JD, Wamosy RMG, Figueiredo F, Parazzi PLF, Schivinski CIS. Impulse oscillometry and obesity in children. J Pediatr. (2018) 94:419–24. doi: 10.1016/j.jped.2017.06.007
46. Criee CP, Sorichter S, Smith HJ, Kardos P, Merget R, Heise D, et al. Body plethysmography–its principles and clinical use. Respir Med. (2011) 105:959–71. doi: 10.1016/j.rmed.2011.02.006
47. Hollier CA, Harmer AR, Maxwell LJ, Menadue C, Willson GN, Black DA, et al. Validation of respiratory inductive plethysmography (LifeShirt) in obesity hypoventilation syndrome. Respir Physiol Neurobiol. (2014) 194:15–22. doi: 10.1016/j.resp.2014.01.014
48. Retory Y, David P, de Picciotto C, Niedzialkowski P, Bonay M, Petitjean M. Validity of thoracic respiratory inductive plethysmography in high body mass index subjects. Respir Physiol Neurobiol. (2017) 242:52–8. doi: 10.1016/j.resp.2017.03.009
49. Lo D, Kennedy JL, Kurten RC, Panettieri RA Jr, Koziol-White CJ. Modulation of airway hyperresponsiveness by rhinovirus exposure. Respir Res. (2018) 19:208. doi: 10.1186/s12931-018-0914-9
50. Lavoie JP, Pascoe JR, Kurpershoek CJ. Effects of xylazine on ventilation in horses. Am J Vet Res. (1992) 53:916–20.
51. Chen YC, Hsiao TC. Instantaneous phase difference analysis between thoracic and abdominal movement signals based on complementary ensemble empirical mode decomposition. Biomed Eng Online. (2016) 15:112. doi: 10.1186/s12938-016-0233-7
52. Broadstone RV, Gray PR, Robinson NE, Derksen FJ. Effects of xylazine on airway function in ponies with recurrent airway obstruction. Am J Vet Res. (1992) 53:1813–7.
53. Lavoie JP, Pascoe JR, Kurpershoek CJ. Effect of head and neck position on respiratory mechanics in horses sedated with xylazine. Am J Vet Res. (1992) 53:1652–7.
Keywords: bronco-alveolar lavage, lung function testing, airway hyper responsiveness, horses, inflammatory airway disease, equine asthma
Citation: Dixon CE, Bedenice D and Mazan MR (2021) Comparison of Flowmetric Plethysmography and Forced Oscillatory Mechanics to Measure Airway Hyperresponsiveness in Horses. Front. Vet. Sci. 7:511023. doi: 10.3389/fvets.2020.511023
Received: 08 November 2019; Accepted: 18 September 2020;
Published: 22 February 2021.
Edited by:
Ali Mobasheri, University of Oulu, FinlandReviewed by:
Angela Briganti, University of Pisa, ItalyMarilena Bazzano, University of Camerino, Italy
Copyright © 2021 Dixon, Bedenice and Mazan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniela Bedenice, ZGFuaWVsYS5iZWRlbmljZSYjeDAwMDQwO3R1ZnRzLmVkdQ==