Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 16 February 2021
Sec. Veterinary Epidemiology and Economics

Modeling Intervention Scenarios During Potential Foot-and-Mouth Disease Outbreaks Within U.S. Beef Feedlots

  • 1Department of Diagnostic Medicine and Pathobiology, College of Veterinary Medicine, Kansas State University, Manhattan, KS, United States
  • 2Center for Outcomes Research and Epidemiology, College of Veterinary Medicine, Kansas State University, Manhattan, KS, United States

Foot-and-mouth disease (FMD) is a highly contagious disease of livestock and has severely affected livestock industries during the past two decades in previously FMD-free countries. The disease was eliminated in North America in 1953 but remains a threat for re-introduction. Approximately 44% of the on-feed beef cattle in the U.S. are concentrated in feedlots <32,000 heads, but little information is available on dynamics of FMD in large feedlots. Therefore, there is a need to explore possible management and intervention strategies that might be implemented during potential FMD outbreaks on feedlots. We used a within home-pen stochastic susceptible-latent-infectious-recovered (SLIR) FMD dynamics model nested in a meta-population model of home-pens in a feedlot. The combinatory model was previously developed to simulate foot-and-mouth disease virus (FMDv) transmission within U.S. beef feedlots. We evaluated three intervention strategies initiated on the day of FMD detection: stopping movements of cattle between home-pens and hospital-pen(s) (NH), barrier depopulation combined with NH (NH-BD), and targeted depopulation of at-risk home-pens combined with NH (NH-TD). Depopulation rates investigated ranged from 500 to 4,000 cattle per day. We evaluated the projected effectiveness of interventions by comparing them with the no-intervention FMD dynamics in the feedlot. We modeled a small-size (4,000 cattle), medium-size (12,000 cattle), and large-size (24,000 cattle) feedlots. Implementation of NH delayed the outbreak progression, but it did not prevent infection of the entire feedlot. Implementation of NH-BD resulted in depopulation of 50% of cattle in small- and medium-size feedlots, and 25% in large-size feedlots, but the intervention prevented infection of the entire feedlot in 40% of simulated outbreaks in medium-size feedlots, and in 8% in large-size feedlots. Implementation of NH-TD resulted in depopulation of up to 50% of cattle in small-size feedlots, 75% in medium-size feedlots, and 25% in large-size feedlots, but rarely prevented infection of the entire feedlot. Number of hospital-pens in the feedlot was shown to weakly impact the success of NH-TD. Overall, the results suggest that stopping cattle movements between the home-pens and hospital-pens, without or with barrier or targeted cattle depopulation, would not be highly effective to interrupt FMDv transmission within a feedlot.

Introduction

Foot-and-mouth disease (FMD) is a highly contagious disease that has affected several non-endemic countries in the past 20 years such as the United Kingdom, Japan, Uruguay, Argentina, the Netherlands, and France (17). In North America, the last outbreaks occurred in 1929 in the United States, 1952 in Canada, and 1946 in Mexico; FMD was eventually eliminated in North America in 1953 (8, 9). Nonetheless, there remains a threat of FMD re-introduction to the U.S. with animals infected or animal products contaminated with the foot-and-mouth disease virus (FMDv). The survival and infectivity of FMDv in fomites and excretions, and aspects of the spread of FMDv has been previously reviewed (1013). An FMD outbreak in the U.S. would cause catastrophic economic consequences for the livestock and associated industries, as has been previously suggested by modeling studies (1416).

At any given time the U.S. has ~13 million cattle on feed distributed in over 30,000 feedlots across the 48 states (17). Approximately 1% of those feedlots have a one-time capacity equal or greater than 32,000 head; however, they contain ~44% of the cattle on-feed population in the country (17). The main strategies used to control FMD during the course of an outbreak in non-endemic countries have been: livestock movement bans, depopulation of infected and susceptible animals in affected and at-risk areas, sanitary/biosecurity measures, surveillance zones, and emergency vaccination (25, 1821). To our knowledge, there are no studies that evaluate the effectiveness of on-farm (within-feedlot) intervention strategies during FMD outbreaks in large concentrated livestock operations such as beef cattle feedlots. The large concentration of cattle in this type of operation might represent a challenge to the success of any of the control strategies mentioned above. For this reason, there is a need to investigate possible management and intervention strategies that might be implemented during potential FMD outbreaks in large feedlots.

Some authors have used models to simulate potential between-farm FMDv transmission within the contiguous U.S., and to evaluate different intervention strategies such as livestock movement bans, depopulation, and vaccination (2226). These models do not focus on within-farm transmission dynamics, projecting which is necessary to assess on-farm management and control strategies in large compartmentalized feedlots. Others have used alternative methods to evaluate the potential use of control strategies at the farm level. For instance, the feasibility of cattle depopulation within a large feedlot has been described by McReynolds and Sanderson (27). They used a Delphi survey and facilitated expert discussion to investigate the complete depopulation of a large feedlot considering effectiveness, animal and human welfare, public perception, and availability of needed supplies. The study concluded that humane and timely depopulation of a large feedlot would be difficult. Given the difficulty of depopulation of large U.S. beef feedlots, alternative FMD control strategies may be required. In this study, we evaluated the projected impacts on FMD outbreak progression in a feedlot of several alternative intervention strategies. We used a mathematical model previously developed to simulate the FMDv transmission and clinical manifestation dynamics within U.S. beef feedlots. We investigate the projected effectiveness of intervention strategies such as within-feedlot cattle movement ban, and a barrier or targeted depopulation of subpopulations of the cattle after FMD detection in the feedlot.

Materials and Methods

Within-Farm FMD Dynamics Model and Feedlot Layouts

We used a previously published and fully described model by Cabezas et al. (28). The FMD infection and clinical disease dynamics within each home-pen and between home-pens are represented in Equations 1–7 and the model parameters are described in Table 1. Briefly, we modeled beef finishing cattle in five different feedlot layouts (see Table 2, and Supplementary Figures 1–5). The housing system most commonly represented in U.S. beef feedlots is cattle housed within pens (home-pens) with metal fences that allow cattle in contiguous home-pens to have direct nose to nose contact through fence-lines. These feedlots are open air with no roof and dirt floors. Cattle are monitored daily be pen riders who enter the pen and systematically ride through looking for sick cattle. Cattle from different home-pens can have directed contact at the hospital pen if they are pulled from their home-pens to receive treatment. In addition, direct contact between cattle from different home-pens might occur in the drovers alleys while cattle are pulled from home-pens to the hospital pen. Contiguous home-pens may also share a water trough. In each home-pen, we implemented a modified stochastic SLIR model with three I compartments of subclinical infectious one, subclinical infectious two, and clinical infectious cattle. The home-pens were nested in a home-pen meta-population in the feedlot. Two levels of FMDv transmission were modeled: within a home-pen and between home-pens. Within home-pen transmission was modeled via direct contact assuming homogeneous cattle mixing inside the home-pen. Between home-pen transmission was modeled via direct contact of cattle in hospital-pen(s), fence-line direct contact of cattle in contiguous home-pens, by pen-riders moving between home-pens, waterborne through drinking water troughs shared by home-pens, and airborne. More details about the model formulation and between home-pen transmission routes modeled can be found in Cabezas et al. (28). Transmission via cattle direct contact in the alleys of the feedlot, and via contaminated feed were not modeled. Simulations started with a proportion of FMD-latent cattle in an index home-pen located centrally within the feedlot.

TABLE 1
www.frontiersin.org

Table 1. Definitions and values of parameters used in modeling potential foot-and-mouth disease transmission, infection, and clinical manifestation dynamics on U.S. beef cattle feedlots.

TABLE 2
www.frontiersin.org

Table 2. Description of the feedlot size and layout, and the projected duration of the outbreak for the baseline no intervention scenario and NH [(hospital movement restrictions to stop mixing of cattle from different home-pens in the hospital-pen(s) beginning].

The feedlot layouts were: FS1—small-size feedlot with 4,000 cattle distributed in 20 home-pens, and operating with one hospital-pen; FM1—medium-size feedlot with 12,000 cattle distributed in 60 home-pens, and operating with one hospital-pen; FM2—medium-size feedlot with 12,000 cattle distributed in 60 home-pens, and operating with two hospital-pens (30 home-pens per hospital-pen); FL1—large-size feedlot with 24,000 cattle distributed in 120 home-pens, and operating with two hospital-pens (60 home-pens per hospital-pen); and FL2—large size feedlot with 24,000 cattle distributed in 120 home-pens, and operating four hospital-pen (30 home-pens per hospital-pen). In all feedlots modeled, every home-pen had 200 cattle. For those feedlots that operated with more than one hospital-pen, the hospital-pen received cattle from the section of home-pens that was in closest spatial proximity to the hospital-pen. Rows of pens were separated by a feed alley for delivering feed to each pen on one size and a drovers ally for driving cattle to the hospital on the other. No pens shared a feed trough. Pens across a feed or drovers alley could not have nose to nose contact. The detailed model formulation is described by (28). The feedlot layout diagrams are included in Supplementary Figures 1–5.

Susceptible:

dSdt=βwpS(I1+I2+I3)φSBin(φ(t1)S(t1), p_inf ​​_​​ hpl(t1)){Sβbp(I1+I2+I3)j; j present               0; otherwise}{Sβbp(I1+I2+I3)h; h present                0; otherwise}{Bin(S, 0.5);  j present, shares water-trough with i, and FMDv load in 1 L of the water ID50 per oral                                                                                          0; otherwise} {Bin(S, 0.5);  h present, shares water-trough with i, and FMDv load in 1 L of the water ID50 per oral                                                                                          0; otherwise}    (1){Bin[(FMDv_floorj×σID50 per oral), 0.5]; j present and (FMDv_floorj×σID50 per oral)                                                0; otherwise} {Bin(S, p_airi); k=1nI30               0; otherwise}μS

Latent:

dLdt=βwpS(I1+I2+I3)φL+Bin(φ(t1)S(t1), p_inf ​​_​​ hpl(t1))+{Sβbp(I1+I2+I3)j; j present               0; otherwise}+{Sβbp(I1+I2+I3)h; h present                0; otherwise}+{Bin(S, 0.5);  j present, shares water-trough with i, and FMDv load in 1 L of the water ID50 per oral                                                                                          0; otherwise}+ {Bin(S, 0.5);  h present, shares water-trough with i, and FMDv load in 1 L of the water ID50 per oral                                                                                          0; otherwise}+    (2){Bin[(FMDv_floorj×σID50 per oral), 0.5]; j present and (FMDv_floorj×σID50 per oral)                                                0; otherwise}+ {Bin(S, p_airi); k=1nI30               0; otherwise}δLμL

Subclinical infectious 1:

dI1dt=δL-θI1-φI1+φ(t-1)I1(t-1)-μI1    (3)

Subclinical infectious 2:

dI2dt=θI1-εI2-φI2+φ(t-1)I2(t-1)-μI2    (4)

Clinical infectious:

dI3dt=εI2-γI3-(φ+ς)I3+(φ(t-1)+ς)I3(t-1)-(μ+ψ)I3    (5)

Clinical non-infectious

dCdt=γI3-τC-(φ+ς)C+(φ(t-1)+ς)C(t-1)-(μ+ψ)C    (6)

Recovered

dRdt=τC-φR+φ(t-1)R(t-1)μR    (7)

Intervention Scenarios

The intervention scenarios were applied on the day of FMD detection in the feedlot. The detection was based on observational surveillance of clinical signs by pen-riders which are experienced personnel in feedlots to detect diseased cattle. The detection was assumed to occur when the proportion of FMD clinical cattle in the index home-pen reached a 3% prevalence threshold. The intervention scenarios were applied by modifying the parameter values from the baseline no-intervention scenario models (see Table 1).

Three intervention scenarios were investigated. In the no hospital (NH) scenario, cattle movements from home-pens to the hospital-pens were stopped to prevent mixing of cattle from different home-pens in the hospital-pen(s) beginning on the day of FMD detection. In the no hospital, barrier depopulation scenario (NH-BD), NH was combined with a barrier depopulation. On the day of FMD detection, cattle in the row of home-pens containing the index home-pen and in the home-pens in the adjacent rows were depopulated. No further depopulation was done following completion of the initial barrier depopulation. In the no hospital, targeted depopulation scenario (NH-TD), NH was combined with a trace-back targeted depopulation. On the day of FMD detection, cattle in home-pens that had contact with the hospital-pen(s) within 7 days prior to FMD detection were traced-back, and those home-pens were depopulated. No further depopulation was done following completion of the initial traceback based depopulation. We assumed a baseline depopulation rate of 1,000 cattle per day (or five home-pens each with 200 cattle), but also evaluated the impact of other depopulation rates. For FS1, we also modeled 2,000 cattle per day depopulation rate; for FM1 and FM2, we also modeled 500 and 2,000 cattle per day depopulation rates; and for FL1 and FL2, we also modeled 2,000 and 4,000 cattle per day depopulation rates. See Supplementary Figures 1–5 for a schematic representation of the intervention strategies in the feedlots modeled.

NH Intervention Scenario

The daily pulling rate of cattle from home-pens to the hospital-pen (φ) due to endemic infectious and non-infectious diseases was set to 0 to stop all movement to and mixing in the hospital pen starting on the day of FMD detection.

NH-BD Intervention Scenario

The barrier depopulation started on the day of FMD detection. The cattle mortality rate was set to 100% in a home-pen on the day of depopulation. The home-pens were depopulated in the inside-out order: the index home-pen, then home-pens in the row where the index home-pen is located, and then home-pens in rows adjacent to the index home-pen row. The number of home-pens depopulated each day was constrained to reflect the assumed maximum daily depopulated rate (head/day) for the feedlot.

NH-TD Intervention Scenario

The traceback based target depopulation started on the day of FMD detection. The cattle mortality rate was set to 100% in a home-pen on the day of depopulation. Only home-pens were depopulated that had cattle coming back from the hospital-pens within 7 days of FMD detection. We prioritized depopulation of home-pens based on the spatial location in the feedlots in relation to the index home-pen. Home-pens in the row containing the index home-pen were depopulated first, then home-pens in adjacent rows and so on. The number of home-pens depopulated each day was constrained to reflect the assumed maximum daily depopulated rate (head/day) for the feedlot.

Outbreak Metrics

We evaluated the following metrics in each feedlot modeled:

NH Intervention Scenario

(1) The projected duration of the outbreak for NH compared to the baseline no-intervention scenario defined as the time in days since the introduction of FMD latent cattle until the prevalence of infectious individuals within the feedlot is equal to 0. (2) The projected time to infection of all home-pens for NH compared to the baseline no-intervention scenario. A home-pen was considered infected when at least one animal become FMD latent during the outbreak. (3) The effectiveness of the intervention implemented in preventing FMD spread. We defined effectiveness of the intervention strategies as the percentage of simulations in which FMD transmission was interrupted and no further susceptible cattle in were infected after implementation of the intervention strategy in each feedlot and intervention modeled.

NH-BD Intervention Scenario

The effectiveness of the intervention implemented in preventing FMD spread. We defined effectiveness of the intervention strategies as the percentage of simulations in which FMD transmission was interrupted for each feedlot and intervention modeled. Success was defined as an iteration where no further pens were infected following the barrier depopulation.

NH-TD Intervention Scenario

The effectiveness of the intervention implemented in preventing FMD spread. We defined effectiveness of the intervention strategies as the percentage of simulations in which FMD transmission was interrupted for each feedlot and intervention modeled. Success was defined as an iteration where at least 1 pen remained uninfected at the end of the outbreak.

The baseline model verification and validation processes were previously conducted by Cabezas et al. (28). We finally compared the results of the intervention strategies to those of the baseline no-intervention scenario (results for the baseline no intervention scenario are reported in Cabezas et al. (28).

Model and Its Output Statistical Analysis Implementation

The model was implemented in Vensim® PLE Plus Version 6.4a (Ventana Systems Inc., Harvard, MA, USA). The output figures were done in R using the ggplot package and the schematics of the depopulation interventions in Microsoft Office Power Point® 365 ProPlus (Microsoft, Redmond, WA, USA). The statistical analysis of the model outputs was done in STATA® 13 (StataCorp LP, College Station, TX, USA). We provide median and percentiles of the results to best represent some non-normal outcome distributions.

Sensitivity Analysis

The target parameters included in the sensitivity analysis were the FMD latent, infectious, and subclinical periods, and the beta transmission parameter within the home-pens (see Table 1). We simulated the model for each feedlot size and layout, for a scenario in which FMD latent cattle are introduced in an index home-pen located centrally within the feedlot, and for each intervention scenario. The value of each target parameter was sampled for each of 2,000 Monte Carlo simulations. For each of the other parameters in the model, a single value was used for each of the 2,000 simulations (see Table 1).

For NH, we investigated the effect of the target parameters and the day of FMD detection on the projected duration of the outbreak. Using the outputs of the 2,000 model simulations for the feedlot size and layout, we used the Spearman rank correlation coefficient to estimate the correlation between each of the target parameters and the day of FMD detection, and between the number of home-pens depopulated and the projected outbreak duration. For NH-BD and NH-TD, we investigated the effect of the target parameters along with the number of hospital-pens in the feedlot, the day of FMD detection, and the number of pens depopulated (only in NH-TD) on the effectiveness of the interventions. Using the outputs of the 2,000 model simulations for the feedlot size and layout, we used the Spearman rank correlation coefficient to estimate the association of the target parameters, the day of FMD detection, and the number of home-pens depopulated with the projected outbreak duration. We also estimated the proportion of simulations with uninfected home-pens at the end of the outbreak. For NH-BD, a fixed number of home-pens was depopulated for each feedlot modeled (10 in FS1, and 30 for all others); then, if the remaining home-pens were not infected during the simulations (10 inf FS1, 30 for FM1 and FM2, and 90 for FL1 and FL2), we considered the intervention successful. For NH-TD, the number of home-pens depopulated was variable depending on how many home-pens had contact with the hospital-pen in the 7 days before FMD detection. So, we considered the intervention successful if at least one home-pen remained uninfected at the end of the outbreak. Descriptive statistics for the projected number of depopulated and uninfected home-pens during the simulations were summarized. Finally, we estimated the effect of modeling different depopulation rates on the outcome described above. For FS1, FM1, and FM2 we modeled depopulation rates of 500 and 2,000 cattle per day while we modeled depopulation rates of 2,000 and 4,000 cattle per day for FL1 and FL2.

Results

Outbreak Progression and Duration

The projected outbreak duration was only compared between NH and the baseline no intervention scenario in which depopulation was not implemented. The largest variation in the projected duration of the outbreak when NH was implemented was seen in FM1. Implementation of this strategy was found to significantly increase the projected duration of the outbreak when compared to the baseline no intervention scenario (82 days compared to 58 days for the baseline no intervention scenario). For FM2 and FL1 the projected median duration of the outbreak for the baseline no intervention scenario was 73 days, and implementation of NH increased it by 8 (FM2) to 11 days (FL1). For FS1 and FL2, there were no changes in the median projected duration of the outbreak when NH was implemented compared to the baseline no intervention scenario (49 days for FS1 and 84 days for FL2). See Table 2 for more detailed results.

Time to Infection of All Home-Pens

The projected time to infection was only compared between NH and the baseline no intervention scenario. The projected time to infect all home-pens since FMDv introduction when comparing NH to the baseline no intervention scenario was longest for FM1. All home-pens took a projected median of 22 days to become infected for the no intervention scenario compared to a projected median of 54 days when NH was implemented. The second largest difference was for FL1 in which the projected median time for all home-pens to become infected for the no intervention scenario was 37 days compared to 53 days for NH. For the rest of feedlots modeled the projected median time to infect all home-pens was 40 days for the baseline scenario compared to 54 days for NH (FM2), 46 days for the baseline scenario compared to 54 for NH (FL2), and 15 days for the baseline scenario compared to 18 days for NH (FS1).

Effectiveness of Different Depopulation Strategies Depending on the Maximum Depopulate Rate for the Feedlot

Implementation of NH was unsuccessful in preventing FMD infection in all feedlots modeled. NH delayed the projected time to infect the entire population but eventually all cattle were infected in all feedlots modeled (see Figures 1, 2 and Table 2 for more detailed information).

FIGURE 1
www.frontiersin.org

Figure 1. Line plots comparing the projected foot-and-mouth outbreak curves for clinical cattle count (Y-axis) during each day since FMDv introduction (X-axis) between the baseline no intervention scenario (28) and the hospital movement restriction (NH) scenario for each feedlot of size and layout. FS1—small-size feedlot with one hospital-pen, FM1—medium-size feedlot with one hospital-pen, FM2—medium-size feedlot with two hospital-pens, FL1—large-size feedlot with two hospital-pens; and FL2—large-size feedlot with four hospital-pens. The black dashed lines represent the 50th percentile of the distribution of the projected count of clinical cattle across 2,000 simulated outbreaks for the baseline no intervention scenario for each feedlot size and layout (28), red solid lines represent the 50th percentile of the distribution of the projected count of clinical cattle across 2,000 simulated outbreaks for NH—hospital movement restrictions to stop mixing of cattle from different home-pens in the hospital-pen(s) after the day of FMD detection.

FIGURE 2
www.frontiersin.org

Figure 2. Line plots comparing the projected cumulative number of infected home-pens (Y-axis) during each day since FMDv introduction (X-axis) between the baseline no intervention scenario (28), and the hospital movement restriction (NH) scenario for each feedlot of size and layout. FS1—small-size feedlot with one hospital-pen, FM1—medium-size feedlot with one hospital-pen, FM2—medium-size feedlot with two hospital-pens, FL1—large-size feedlot with two hospital-pens; and FL2—large-size feedlot with four hospital-pens. The black solid lines represent the 50th percentile of the distribution of the projected cumulative number of infected home-pens across 2,000 simulated outbreaks for the baseline no intervention scenario for each feedlot size and layout (28), red solid lines represent the 50th percentile of the distribution of the projected cumulative number of infected home-pens across 2,000 simulated outbreaks for NH—hospital movement restrictions to stop mixing of cattle from different home-pens in the hospital-pen(s) after the day of FMD detection.

Barrier depopulation (NH-BD) was the intervention with the highest probability of success in interrupting FMD infection resulting in no additional infected home-pens after its implementation (Table 3). For feedlots FM1 and FM2 (30 home-pens depopulated) continued transmission of FMDv infection was interrupted (30 home-pens uninfected following barrier depopulation) in 16% of simulations when depopulation was set at 500 cattle per day and in 38-40% of simulations if 1,000 or 2,000 cattle per day were depopulated. For FL1 and FL2 (30 pens depopulated), continued transmission of FMDv infection was interrupted (90 home-pens uninfected following barrier depopulation) in 41–42% of simulations when 4,000 cattle were depopulated per day; 34–38% of simulations when depopulation was set at 2,000 cattle per day but decreased to only 7–8% when depopulation was set at 1,000 cattle per day. NH-BD (10 home-pens depopulated) was never successful in interrupting FMD transmission in FS1 (no remaining uninfected home-pens) for either depopulation rate modeled (1,000 or 2,000 cattle per day) (see Table 3).

TABLE 3
www.frontiersin.org

Table 3. Projected effectiveness of NH-BD (NH combined with barrier depopulation) and NH-TD (NH combined with targeted depopulation) for the feedlot size and layout and the different depopulation rates modeled, 2000 simulations.

Targeted depopulation (NH-TD) was less effective compared to NH-BD (Table 3). In no case did targeted depopulation result in no additional infected home-pens following implementation. For this reason, success for NH-TD was defined as any uninfected home-pen at the end of the outbreak. For FM1, if depopulation was set at 500 cattle per day, the percentage of successful simulations was 82%; however, 43 home-pens (72%) were depopulated and a median number of 2 uninfected home-pens remained after the intervention. The percentage of successful simulations increased to over 90% and a median of 15 (25%) and 18 (30%) home-pens remain uninfected when 1,000 or 2,000 cattle were depopulated per day. For FM2, FL1, and FL2, the intervention was poorly effective with rare uninfected home-pens after the intervention with a median number of 25, 33, and 27 home-pens depopulated regardless the number of cattle depopulated per day. NH-TD was never successful in FS1 with a median number of 11 home-pens depopulated in the intervention (see Table 3).

Sensitivity Analysis

For NH, the sensitivity analysis showed that the duration of the latent period had a strong positive correlation (Spearman rank correlation coefficient ≥0.77 for all feedlots modeled) with the projected duration of the outbreak. The infectious period was moderately positively correlated (Spearman rank correlation coefficient was ~0.28 for FM1, FM2, FL1, and FL2 and 0.47 for FS1) with the projected duration of the outbreak. The subclinical period had a negative weak correlation with the projected duration of the outbreak in all feedlots modeled. The beta transmission parameter did not have a significant association (p > 0.05) with the projected duration of the outbreak for any of the feedlots modeled (see Table 4).

TABLE 4
www.frontiersin.org

Table 4. Target parameters investigated for associations with the projected duration of the outbreak after implementation of NH (hospital movement restrictions to stop mixing of cattle from different home-pens in the hospital-pen(s) beginning the day after FMD).

We summarized only the results of FM1, FM2, FL1, and FL2 sensitivity analysis for the depopulation scenarios because the interventions were never successful in interrupting FMDv transmission in FS1. For both NH-BD and NH-TD, simulations showed that the duration of FMD stages (latent, infectious, and subclinical period) were weakly correlated with having uninfected home-pens after the interventions for all of the feedlots and for all of the depopulation rates modeled. The number of hospital-pens did not show a significant correlation (p > 0.05) with having uninfected home-pens for any depopulation rate after NH-BD was implemented in medium-size feedlots. In contrast, it was found to be significantly (p < 0.05) and moderately correlated with having uninfected home-pens after NH-TD was implemented—The higher the depopulation rates the stronger the correlation. In large-size feedlots, for both NH-BD and NH-TD, the number of hospital pens were weakly correlated with having uninfected home-pens (Spearman rank correlation coefficient at most 0.05 for any of the feedlots and for any of the depopulation rates modeled).

The day of FMD detection was moderately negatively correlated with having uninfected home-pens after NH-BD for all feedlots modeled at higher depopulation rates (1,000 and 2,000 cattle per day for medium-size feedlots, and 2,000 and 4,000 cattle per day for large-size feedlots). In contrast, the day of FMD detection was weakly and positively correlated (ranged from 0.5 to 0.24) with having uninfected home-pens after NH-TD for all feedlots and all depopulation rates (see Table 5).

TABLE 5
www.frontiersin.org

Table 5. Target parameters investigated for associations with having uninfected home-pens after implementation of NH-BD (NH was combined with barrier depopulation) and NH-TD (NH was combined with targeted depopulation) for the feedlot size and layouts modeled.

Discussion

To our knowledge, our model is the first to describe the application of on-farm intervention strategies in the face of a potential FMD outbreak in U.S. beef feedlots. Based on our knowledge of the feedlot production system and previous experience of FMD epidemics in non-endemic countries, we evaluated the impact of movement restriction within the feedlot during the outbreak, and two partial depopulations strategies combined with movement restrictions on outbreak progression.

The interventions modeled had no effect on the projected duration of the outbreak and eventual infection of the entire feedlot (NH) or the number of remaining uninfected home-pens (NH-BD and NH-TD) in small sized FS1 feedlots. While feedlots with one-time head capacity of 4,000 or less represent ~27% of the cattle on-feed population, they represent up to 97% of the feedlots in the country (17). This model suggests that partial or targeted depopulation may have little effect on disease in these feedlots. NH was found to significantly decrease the projected outbreak progression in feedlots that operated with more home-pens per hospital-pen such as FM1 and FL1 (60 home-pens per hospital-pen). In feedlots that operated with fewer home-pens per hospital-pen such as FM2 and FL2 (30 home-pens per hospital-pen), the projected outbreak progression was delayed by NH but not as much as for FM1 and FL1. However, it is important to emphasize that, though delayed, the entire population in these feedlots was still infected. NH may be a useful strategy in medium- and large-size feedlots to delay infection progression while preparing logistics for other intervention strategies such as vaccination, however this was not assessed by the current model. In reality, complete stoppage of mixing of cattle from different home-pens in a hospital system might not be feasible. Any attempt to do so would likely require treatment within home-pens or use of portable hospital facilities that could move between home-pens. The impact of increased entry into home-pens or the use of portable hospital-pens on transmission within the feedlot was not assessed in this model. Also, complementary interventions to improve the success of movement restrictions such as well-defined biocontainment practices should be developed in advance as suggested by Brandt et al. (35). Temporary movement restriction during the period of targeted depopulation is also not assessed in this model but could be implemented.

Implementation of NH-BD (NH combined with barrier depopulation) under our assumptions was partially effective on medium- (50% of home-pens uninfected after the intervention) and large-size feedlots (25% of home-pens uninfected after the intervention) when higher depopulation rates were implemented (1,000 or more cattle per day for medium-size feedlots, and 2,000 or more cattle per day for large-size feedlots). However, this was following a depopulation of 50% of home-pens. We used an inside-out strategy of depopulation which means that depopulation was conducted starting with the index home-pen and then home-pens surrounding the index home-pen. An outside-in strategy could be explored to assess if there is any advantage.

Implementation of NH-TD (NH combined with targeted depopulation) was poorly effective in all feedlots modeled. In FM1, the intervention was partially successful to prevent infection in up to 30% of home-pens when the depopulation rate was 1,000 or 2,000 cattle per day, although 65–70% of home-pens in the feedlot had to be depopulated. For FM2, ~35–50% of home-pens were depopulated, and for FL1 and FL2 between 20 and 30% were depopulated but only a few uninfected home-pens were present at the end of the outbreak in successful simulations. We highlight that for all strategies we modeled an optimistic day of FMD detection of 3% clinical animals in the index home-pen by observational surveillance of pen-riders. While pen-riders are experienced personnel in detecting diseased animals (36, 37), clinical signs of FMD are very similar to other diseases which may confuse detection (38). Initial clinical detection will be followed by laboratory confirmation which can take up to several days depending on the logistics to collect and ship samples, and conduct the required tests to confirm FMDv suspicion as discussed by Sutmoller et al. (39) in their description of FMD outbreaks in the early 2000s. However, another study conducted by Walz et al. (40) modeled a detection threshold of 5% prevalence of clinical animals in beef herds with 5,000–50,000 one-time head capacity and conducted a sensitivity analysis testing 2.5 and 10% detection thresholds and found that the time-to-detection was not sensitive to those changes. Other methods for early detection such as the use of a surveillance test are not currently available but could be explored in the future. Other authors have discussed the potential use of real-time polymerase chain reaction (Rt-PCR) to test the saliva of animal in ropes in pens as a surveillance method to detect FMDv during the pre-clinical stage (41). Our model suggests that even with optimistic early detection of FMD within the feedlot, modeled methods are not sufficient to reliably stop an outbreak. Since interventions following our early detection were not sufficient, we did not model later detection times.

Depopulation strategies in previous outbreaks in Ireland, the Netherlands, and the UK attempted depopulation of affected cattle and/or susceptible cattle within 2 days (19, 42, 43). This might be feasible for these countries where the average herd size is <100 cattle (44); however, the experience in the UK demonstrated that the implementation of this policy was difficult to achieve (3, 45, 46). A requirement to complete depopulation on a large feedlot in 2 days is likely unrealistic even for the partial depopulations modeled by the current model. McReynolds and Sanderson (27) conducted a survey to investigate the feasibility of depopulation in large feedlots during a health emergency event such as an FMD outbreak and concluded that the methods explored were not viable to ensure a rapid, safe and humane depopulation.

We used a base depopulation rate of 1,000 cattle per day but also modeled higher depopulation rates than 1,000 cattle per day. We found that higher depopulation rates did not result in substantial differences in projected modeled outcomes compared to the base depopulation rate. Implementing higher depopulation rates than 1,000 cattle per day could be difficult to achieve, depending on available facilities. Hence, larger daily culling capacities modeled here may be optimistic. Moreover, facilities available in feedlots to dispose of depopulated carcasses might play a large role to limit rapid depopulation, however this limitation was not assessed in the current model. The cost of implementation should be explored in the future to make a more informed decision about the feasibility of implementation of partial or targeted depopulation strategies.

Another important factor to consider in our model is that for NH-TD we used a 100% accurate traceback prior to implementation of the depopulation strategies. Even good record-keeping in the feedlot likely will not achieve this level of trace-back accuracy. Since this level of accuracy was not effective, we did not explore less accurate traceback. There is only one report, to our knowledge, that addresses intervention strategies in a large feedlot (>14,000 cattle capacity) in South Africa (47). The authors reported that they adopted vaccination instead of depopulation due to the difficulties in maintaining bio-security measures during the depopulation of large number of animals. Other authors have suggested the potential implementation of selective depopulation which requires the culling of affected cattle based on the presence of clinical signs (27); however, the proportion of cattle to develop clinical FMD in a totally naïve population should be expected to be very high. An expert survey of FMD related parameters and clinical manifestation suggested that ~65–80% of cattle in U.S. beef feedlots might develop clinical FMD if infected by a high or low strain virulent strain, respectively (48).

In the sensitivity analysis, we found that the NH models were sensitive to changes in the latent and infectious periods. This suggests that introduction of high or low virulence strains could have an effect in the projected outbreak progression within the feedlots modeled. However, for NH-BD and NH-TD models, the duration of latent, infectious, and subclinical stages were less influential and other parameters had a larger effect on the projected outputs presented. For NH-BD, the projected day of FMD detection was the most influential parameters on probability of success in all feedlots modeled. This is not surprising since increased time to detection results in more time for spread between home-pens prior to implementation of interventions. However, for NH-TD, longer times for FMD detection were found to be associated with a higher likelihood to having uninfected home-pens. This can be explained by the fact that more home-pens have had contact with the hospital-pen in simulations with longer times to FMD detection, and therefore having uninfected home-pens after NH-TD might have been confounded by the number of home-pens depopulated by the intervention. For future models, different values for sensitivity of the observational surveillance, and delay in implementation of intervention due to FMD laboratory confirmation should be explored although we found that modeling an optimistic FMD detection threshold was already too late for the interventions to be highly successful. Finally, this model does not assess the potential effectiveness of vaccination and/or its combination with the intervention strategies modeled as a control option, and future work should evaluate the feasibility to implement vaccination.

Conclusions

We believe we have captured the important structure and management aspects of U.S. feedlot systems and best estimates of FMD transmission parameters. Even with some optimistic assumptions, the three intervention strategies modeled were not highly effective in controlling the outbreak or required depopulation of a large proportion of cattle. Still, the results of our model should be interpreted with caution. Little data is available to inform the biological behavior of FMD in an immunologically naïve cattle population in confined production systems. Refinement of the methods used to model shedding and transmission along with better quality of data is needed to produce more robust models. The strategies should also be measured with a financial component to evaluate the cost-effectiveness of their implementation. Restriction of cattle movements from home-pens to hospital-pen proved to considerably prolong the outbreak in larger feedlots. Strategies combining vaccination with such movement restriction or targeted depopulation should be investigated. Finally, the exploration of different intervention strategies is challenging in beef feedlots in the U.S. because there are few other countries in the world with a similar production system and the immunologically naïve cattle population, so there is substantial uncertainty in how severe an FMD outbreak will be if the virus is introduced into the country.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author's Note

This work was adapted from a doctoral dissertation chapter available at: https://krex.k-state.edu/dspace/handle/2097/40307.

Author Contributions

MS and VV conceived and designed the study. AC implemented the models and performed the sensitivity analyses. All authors contributed to the development, implementation, and analysis of the models and the output interpretation. All authors wrote the manuscript, read, and approved the final version for publication.

Funding

This material is based upon work supported by the Kansas Bioscience Authority under the federal Matching Program Grant Agreement.

Disclaimer

The views and conclusions contained in this manuscript are those of the authors and should not be interpreted as necessarily representing the official policies, either expressed or implied, of the funding agency. This article was prepared while VV was employed at Kansas State University. The opinions expressed in this article are the author's own and do not reflect the view of the National Institutes of Health, the Department of Health and Human Services, or the United States government.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.559785/full#supplementary-material

References

1. Anderson I. Foot and Mouth Disease 2001: Lessons to be Learned Inquiry Report. London: London: The Stationary Office (2002).

Google Scholar

2. Chmitelin I, Moutou F. Foot and mouth disease: lessons to be learned from the experience of France. Rev Sci Tech Off Int Epiz. (2002) 21:731–7. doi: 10.20506/rst.21.3.1370

CrossRef Full Text | Google Scholar

3. Davies G. The foot and mouth disease (FMD) epidemic in the United Kingdom 2001. Comp Immunol Microbiol Infect Die. (2002) 25:331–43. doi: 10.1016/S0147-9571(02)00030-9

CrossRef Full Text | Google Scholar

4. Bouma A, Elbers ARW, Dekker A, De Koejier A, Bartels C, Vellema P, et al. The foot-and-mouth disease epidemic in the Netherlands in 2001. Prev Vet Med. (2003) 57:155–66. doi: 10.1016/S0167-5877(02)00217-9

CrossRef Full Text | Google Scholar

5. Nishiura H, Omori R. An epidemiological analysis of the foot-and-mouth disease epidemic in Miyazaki, Japan, 2010. Transbound Emerg Dis. (2010) 57:396–403. doi: 10.1111/j.1865-1682.2010.01162.x

CrossRef Full Text | Google Scholar

6. Ellis-Iversen J, Smith RP, Gibbens JC, Sharpe CE, Dominguez M, Cook AJ. Risk factors for transmission of foot-and-mouth disease during an outbreak in southern England in 2007. Vet Rec. (2011) 168:128. doi: 10.1136/vr.c6364

CrossRef Full Text | Google Scholar

7. Clavijo A, Sanchez-Vazquez MJ, Buzanovsky LP, Martini M, Pompei JC, Cosivi O. Current status and future prospects to achieve foot-and-mouth disease eradication in South America. Transbound Emerg Dis. (2017) 64:31–6. doi: 10.1111/tbed.12345

CrossRef Full Text | Google Scholar

8. Mohler JR. The 1929 outbreak of foot-and-mouth disease in Southern California. In: U.S.D.o. Agriculture. Washington, DC: US. Dept. of Agriculture (1930).

Google Scholar

9. Graves JH. Foot and mouth disease: a constant threat to US livestock. JAVMA. (1979) 174:174–6.

Google Scholar

10. Bachrach HL. Foot-and-mouth disease. Annu Rev Microbiol. (1968) 22:201–44. doi: 10.1146/annurev.mi.22.100168.001221

CrossRef Full Text | Google Scholar

11. Sellers R. Quantitative aspects of the spread of foot and mouth disease. Vet Bull. (1971) 41:431–9.

Google Scholar

12. Bartley LM, Donnelly CA, Anderson RM. Review of foot-and-mouth disease virus survival in animal excretions and fomites. Vet Rec. (2002) 151:667–9. doi: 10.1136/vr.151.22.667

CrossRef Full Text | Google Scholar

13. Alexandersen S, Zhang Z, Donaldson AI, Garland AJM. The pathogenesis and diagnosis of foot-and-mouth disease. J Comp Pathol. (2003) 129:1–36. doi: 10.1016/S0021-9975(03)00041-0

CrossRef Full Text | Google Scholar

14. Paarlberg PL, Lee JG, Seitzinger AH. Potential revenue impact of an outbreak of foot-and-mouth disease in the United States. JAVMA. (2002) 220:988–92. doi: 10.2460/javma.2002.220.988

CrossRef Full Text | Google Scholar

15. Pendell DL, Marsh TL, Coble KH, Lusk JL, Szmania SC. Economic assessment of FMDv releases from the national bio and agro defense facility. PLoS ONE. (2015) 10:e0129134. doi: 10.1371/journal.pone.0129134

CrossRef Full Text | Google Scholar

16. Schroeder TC, Pendell DL, Sanderson MW, McReynolds S. Economic impact of alternative fmd emergency vaccination strategies in the midwestern United States. J. Agric Appl Econ. (2015) 47:47–76. doi: 10.1017/aae.2014.5

CrossRef Full Text | Google Scholar

17. USDA-NASS. (2017). Overview of U.S. Livestock, Poultry, and Aquaculture Production in 2016. Fort Collins, CO: USDA-APHIS.

Google Scholar

18. Kao RR. The role of mathematical modelling in the control of the 2001 FMD epidemic in the UK. Trends Microbiol. (2002) 10:279–86. doi: 10.1016/S0966-842X(02)02371-5

CrossRef Full Text | Google Scholar

19. Pluimers FH, Akkerman AM, Van der Wal P, Dekker A, Bianchi A. Lessons from the foot and mouth disease outbreak in the Netherlands in 2001. Rev Sci Tech Off Int Epiz. (2002) 21:711–21. doi: 10.20506/rst.21.3.1371

CrossRef Full Text | Google Scholar

20. Bouma A, Dekker A, De Jong MC. No foot-and-mouth disease virus transmission between individually housed calves. Vet Microbiol. (2004) 98:29–36. doi: 10.1016/j.vetmic.2003.10.016

CrossRef Full Text | Google Scholar

21. Park JH, Lee KN, Ko YJ, Kim SM, Lee HS, Shin YK, et al. Control of foot-and-mouth disease during 2010-2011 epidemic, South Korea. Emerg Infect Dis. (2013) 19:655–9. doi: 10.3201/eid1904.121320

CrossRef Full Text | Google Scholar

22. Bates TW, Thurmond MC, Carpenter TE. Description of an epidemic simulation model for use in evaluating strategies to control an outbreak of foot-and-mouth disease. AJVR. (2003) 64:195–204. doi: 10.2460/ajvr.2003.64.195

CrossRef Full Text | Google Scholar

23. Schoenbaum MA, Disney TW. Modeling alternative mitigation strategies for a hypothetical outbreak of foot-and-mouth disease in the United States. Prev Vet Med. (2003) 58:25–52. doi: 10.1016/S0167-5877(03)00004-7

CrossRef Full Text | Google Scholar

24. Ward MP, Highfield LD, Vongseng P, Graeme Garner M. Simulation of foot-and-mouth disease spread within an integrated livestock system in Texas, USA. Prev Vet Med. (2009) 88:286–97. doi: 10.1016/j.prevetmed.2008.12.006

CrossRef Full Text | Google Scholar

25. Tildesley MJ, Smith G, Keeling MJ. Modeling the spread and control of foot-and-mouth disease in pennsylvania following its discovery and options for control. Prev Vet Med. (2012) 104:224–39. doi: 10.1016/j.prevetmed.2011.11.007

CrossRef Full Text | Google Scholar

26. McReynolds SW, Sanderson MW, Reeves A, Hill AE. Modeling the impact of vaccination control strategies on a foot and mouth disease outbreak in the central United States. Prev Vet Med. (2014) 117:487–504. doi: 10.1016/j.prevetmed.2014.10.005

CrossRef Full Text | Google Scholar

27. McReynolds SW, Sanderson MW. Feasibility of depopulation of a large feedlot during a foot-and-mouth disease outbreak. JAVMA. (2014) 244:291–8. doi: 10.2460/javma.244.3.291

CrossRef Full Text | Google Scholar

28. Cabezas AH, Sanderson MW, Volkova VV. A meta-population model of potential foot-and-mouth disease transmission, clinical manifestation, detection within U.S. beef feedlots. Front Vet Sci. (2020) 7:527558. doi: 10.3389/fvets.2020.527558

CrossRef Full Text | Google Scholar

29. Chis Ster I, Dodd PJ, Ferguson NM. Within-farm transmission dynamics of foot and mouth disease as revealed by the 2001 epidemic in great Britain. Epidemics. (2012) 4:158–69. doi: 10.1016/j.epidem.2012.07.002

CrossRef Full Text | Google Scholar

30. Mardones F, Perez A, Sanchez J, Alkhamis M, Carpenter T. Parameterization of the duration of infection stages of serotype O foot-and-mouth disease virus: an analytical review and meta-analysis with application to simulation models. Vet Res. (2010) 41:45. doi: 10.1051/vetres/2010017

CrossRef Full Text | Google Scholar

31. DEFRA. (2005). Foot and Mouth Disease Ageing of Lesions. London: DEFRA.

Google Scholar

32. USDA. “Part IV: Health Management on U.S. Feedlots with a Capacity of 1,000 or More Head. (USDA–APHIS–VS–CEAH–NAHMS). Fort Collins, CO: USDA-APHIS.

Google Scholar

33. Schijven J, Rijs GBJ, De Roda Husman AM. Quantitative risk assessment of FMD virus transmission viawater. Risk Anal. (2005) 25:9. doi: 10.1111/j.0272-4332.2005.00563.x

CrossRef Full Text | Google Scholar

34. Boender GJ, van Roermund HJ, de Jong MC, Hagenaars TJ. Transmission risks and control of foot-and-mouth disease in The Netherlands: spatial patterns. Epidemics. (2010) 2:36–47. doi: 10.1016/j.epidem.2010.03.001

CrossRef Full Text | Google Scholar

35. Brandt AW, Sanderson MW, DeGroot BD, Thomson DU, Hollis LC. Biocontainment, biosecurity, and security practices in beef feedyards. JAVMA. (2008) 232:262–9. doi: 10.2460/javma.232.2.262

CrossRef Full Text | Google Scholar

36. Smith RA, Stokka GL, Radostits OM, FGriffin DD. Health and Production Management in Beeflots. Philadelphia, PA: W.B. Saunders (1994).

Google Scholar

37. Center for Food Security Public Health USDA United States Department of Agriculture and Service AA.a.P.H.I. FAD PReP Beef Feedlot Industry Manual. Ames, IA: Iowa State University of Science and Technology, USDA-APHIS (2011).

Google Scholar

38. Coetzer JAW, Tustin RC. Infectious Diseases of Livestock. Cape Town: Oxford University Press (2004).

Google Scholar

39. Sutmoller P, Barteling SS, Casas Olascoaga R, Sumption K. Control and eradication of foot-and-mouth disease. Virus Res. (2003) 91:101–44. doi: 10.1016/S0168-1702(02)00262-9

CrossRef Full Text | Google Scholar

40. Walz E, Middleton J, Sampedro F, VanderWaal K, Malladi S, Goldsmith T. Modeling the transmission of foot and mouth disease to inform transportation of infected carcasses to a disposal site during an outbreak event. Front Vet Sci. 6:501. doi: 10.3389/fvets.2019.00501

CrossRef Full Text | Google Scholar

41. Nelson N, Paton DJ, Gubbins S, Colenutt C, Brown E, Hodgson S, et al. Predicting the ability of preclinical diagnosis to improve control of farm-to-farm foot-and-mouth disease transmission in cattle. J Clin Microbiol. (2017) 55:1671–81. doi: 10.1128/JCM.00179-17

CrossRef Full Text | Google Scholar

42. Costelloe JA, Gaynor MC, Gaynor S, McAteer MJ, O'Reilly PJ. Control of foot and mouth disease: lessons from the experience of Ireland. Rev Sci Tech Off Int Epiz. (2002) 21:739–50. doi: 10.20506/rst.21.3.1369

CrossRef Full Text | Google Scholar

43. Scudamore JM, Harris DM. Control of foot and mouth disease: lessons from the experience of the outbreak in Great Britain in 2001. Rev Sci Tech Off Int Epiz. (2002) 21:699–710. doi: 10.20506/rst.21.3.1351

CrossRef Full Text | Google Scholar

44. AHBD. UK Yearbook 2016 Cattle. Kenilworth, UK: AHBD (2016).

Google Scholar

45. Ferguson NM, Donnelly CA, Anderson RM. The foot-and-mouth epidemic in great Britain: pattern of spread and impact of interventions. Science. (2001) 292:1155–60. doi: 10.1126/science.1061020

CrossRef Full Text | Google Scholar

46. Gibbens JC, Sharpe CE, Wilesmith JW, Mansley LM, Michalopoulou E, Ryan JBM, et al. Descriptive epidemiology of the 2001 foot-and-mouth disease epidemic in great Britain: the first five months. Vet Rec. (2001) 149:729–43. doi: 10.1136/vr.149.24.729

CrossRef Full Text | Google Scholar

47. Bruckner VK, Vosloo W, Du Plessis BJA, Kloeck PELG, Connoway L, Ekron MD, et al. Foot and mouth disease: the experience of South Africa. Rev Sci Tech Off Int Epiz. (2002) 21:751–64. doi: 10.20506/rst.21.3.1368

CrossRef Full Text | Google Scholar

48. Cabezas AH, Sanderson MW, Jaberi-Douraki M, Volkova VV. Clinical and infection dynamics of foot-and-mouth disease in beef feedlot cattle: an expert survey. Prev Vet Med. (2018) 158:160–8. doi: 10.1016/j.prevetmed.2018.08.007

CrossRef Full Text | Google Scholar

Keywords: foot-and-mouth disease, meta-population model, beef feedlot, beef cattle, intervention strategies, cattle depopulation

Citation: Cabezas AH, Sanderson MW and Volkova VV (2021) Modeling Intervention Scenarios During Potential Foot-and-Mouth Disease Outbreaks Within U.S. Beef Feedlots. Front. Vet. Sci. 8:559785. doi: 10.3389/fvets.2021.559785

Received: 07 May 2020; Accepted: 25 January 2021;
Published: 16 February 2021.

Edited by:

Francisco Ruiz-Fons, Consejo Superior de Investigaciones Científicas (CSIC), Spain

Reviewed by:

Anette Ella Boklund, University of Copenhagen, Denmark
Elisha Are, Stellenbosch University, South Africa

Copyright © 2021 Cabezas, Sanderson and Volkova. 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: Victoriya V. Volkova, vv88@vet.k-state.edu; Michael W. Sanderson, sandersn@vet.k-state.edu

Present address: Victoriya V. Volkova, National Institutes of Health, Bethesda, MD, United States

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.