Disease Transmission in Cattle and Horses: Spatial Analysis of Landscape Features Associated with Vesicular Stomatitis

 

Susan A. Maroney and Jerome E. Freier

U.S. Department of Agriculture, Animal and Plant Health Inspection Service

Centers for Epidemiology and Animal Health

 

Abstract

 

Three major epidemics of vesicular stomatitis have occurred recently in the southwestern United States. Infection with vesicular stomatitis virus causes lesions to form around the head and hooves, making it difficult for infected animals to feed and results in significant weight losses. This disease can have important economic impacts on both dairy and beef cattle industries in the form of reduced productivity and the prevention of animal exports through international trade restrictions. GIS and spatial analysis methods were applied to individual case sites to establish patterns of clustering and to identify landscape features possibly correlated with disease transmission. Vegetation, terrain, hydrologic, and transportation factors were analyzed. Several spatial features were identified as common to infected animal locations and this information is now being used to predict the possible direction and spread of future outbreaks of this disease.

 

Introduction

 

Background

Vesicular stomatitis virus (VSV)  (Family Rhabdovirus, Genus Vesiculovirus) occurs in North and South America. This virus is endemic in Mexico and Central America and occurs sporadically in the U.S. There are two serotypes of VSV: New Jersey and Indiana. Immunity to one serotype does not protect against infection with the other. VSV primarily affects equines and bovines, and manifests itself in fever and blister-like lesions on the lips, tongue, palate, teats, and/or hooves of the affected animal. The result is lethargy, difficulty eating or suckling young, and weight loss. There is no specific treatment, but the disease is self-limiting and resolves itself in about two weeks (Fraser 1991). Animals are quarantined while clinical symptoms are present. The disease can be transmitted by arthropod vectors or by mechanical methods, including direct contact.

There are two reasons for concern about this disease: its similarity to Foot-and-Mouth Disease and its economic impacts. The symptoms of VSV are very similar to those of Foot-and-Mouth Disease, a highly contagious devastating bovine disease eradicated from the U.S. in 1929. In all susceptible animals, ELISA testing is used to confirm that the animal was infected with VSV. VSV causes decreased productivity in beef and dairy cattle and weight loss in cattle and horses. The presence of VSV in a herd triggers export trade restrictions preventing U.S producers from selling animals abroad. Domestically, noninfected states restrict the movement of animals from the infected states (USDA 1996).

Three epidemics of Vesicular Stomatitis Virus (VSV) have occurred in the U.S. in the past five years, in 1995, 1997, and 1998. The focus of this study is the 1995 outbreak. In 1995, the states involved were New Mexico, Colorado, Wyoming, Utah, Texas, and Arizona.  There were 329 confirmed cases in horses, cattle, three sheep, one pig, and one llama. The outbreak consisted entirely of New Jersey strain virus and lasted from May to November 1995. Cases were reported to a USDA emergency operations center. An investigator was sent to collect tissue samples from each suspected case. Geographic coordinates for each case were determined. Cases of the disease have not been reported in the U.S. in the past eighteen months. It is important to note that in all three outbreaks the only animals tested were those that exhibited clinical signs of disease. Unaffected animals on a premises at the time of the veterinarian's visit were not tested. Also, reporting of VSV cases to veterinary health authorities is voluntary. As a result, the true prevalence of this disease in the affected area is unknown.

Current Study

This study focused on the 1995 VSV outbreak consisting of 329 cases occurring in New Mexico, Colorado, Wyoming, Utah, Texas, and Arizona (Figure 1).

Figure 1. Extent of the 1995 VSV Outbreak

 

 

The mission of the USDA Veterinary Services is to protect and promote animal agriculture in the U.S.  Integral to this mission is tracking of livestock animal diseases. The intent of this study is to create a landscape epidemiology model, which relates disease occurrence to environmental factors. In such a model, spatial features like terrain, vegetation, and hydrology that are characteristic of the disease areas, are defined and used to identify areas at risk for new outbreaks.

Little is known about the method of transmission for VSV.  Insects as a vector of transmission is suspected, perhaps black flies (Diptera: Simuliidae) or phlebotomine sand flies (Diptera: Pyschodidae). Distributions and densities of vectors are difficult to obtain. These vectors usually live on a wide range of host animals including horses, prairie dogs, and wood rats. Mechanical means of transmission include transportation of infected animals and by fomites. The disease spreads rapidly with animal-to-animal contact.

With the cessation of new cases coinciding with the onset of cold weather, it appears that outbreaks are seasonal. This is indicative of a vector-borne disease. The first cases in a year are typically seen in late May and the last ones diagnosed in early November. These months coincide with the end and then the onset of cold weather in the Southwest.

 

Data Collection

Case and Cluster Definitions

Cases were identified from reports to veterinarians by owners of sick animals. The premises were located by the coordinates of the front gate. Case locations were documented at the time of the veterinarian's visit with a GPS receiver, or with a digital atlas after the fact. Coordinates collected with a GPS receiver were not differentially corrected. Polygon fitting was used to identify clusters from the complete data set of case points. Three case clusters were identified using a minimum polygon method (Figure 2). The minimum polygon algorithm of polygon fitting is a parametric statistical method for “connecting the dots.”  It draws line segments to connect the outermost dots, after defining outliers and setting a maximum line segment length.   The three case clusters identified (called the Rio Grande Cluster, the Grand Junction Cluster, and the Cortez Cluster) are the basis for the remainder of the data gathering. The Rio Grande cluster contains 173 cases, the Grand Junction 104 cases, and the Cortez cluster 21 cases. 

Figure 2. Clusters of cases in the 1995 VSV outbreak

Vegetation

Normalized Difference Vegetation Index (NDVI) Biweekly Composite images derived from Advanced Very High Resolution Radiometer (AVHRR) data were acquired from the USGS for March-December 1995. Brightness values, indicative of vegetation biomass, are calculated every two weeks from the visible band and the near infrared band to create the NDVIs.  The brightness value for each case location in each cluster is extracted from each NDVI for the period March - December 1995 and exported to a database table. The resolution of NDVIs is 1.1 kilometers. Although small scale for environmental analysis, NDVIs can give an idea of biomass changes in the area around each case over the months in question.

 

Hydrology

For each cluster a hydrology data layer was imported from Bureau of the Census TIGER/Line files for the counties containing VSV cases. ArcView's Geoprocessing Wizard was used to calculate distance from each case in the cluster to the nearest water feature.  Many of the TIGER water features are identified as ditches or canals, as differentiated from rivers and streams. This distinction is important in locating potential vector habitats; difference in velocity of water flow and vegetation on the banks create suitable habitat for different vectors.

Terrain

Terrain features include elevation, slope, and aspect of each case in the three clusters. USGS digital elevation model (DEM) data were acquired for each cluster.  They were mosaicked using the Mosaic GRID Themes extension for ArcView found on the Esri web site.  These 10-meter resolution data layers provide point elevation data for each case in a cluster. The elevation value of the cell under each case was exported to a table. In addition, slope and aspect for the cell under each case location were calculated and exported to a table in preparation for spatial analysis.

Spatial Statistical Analysis

The data layers described above were all compiled into several database tables for use as input to statistical analysis. Preliminary statistics have been completed to date; detailed analysis and comparison with other outbreaks are continuing. Tools used include ArcView with Spatial Analyst, S+ for ArcView, Lotus Approach, and EpiInfo (a database and statistics program for epidemiology).

 

Temporal Analysis

Figure 3 reveals a temporal aspect to the cases. In this figure, the cases in all three clusters are graphed against quarantine date. Note how the early cases (June and July 1995) appear at lower latitudes (more southerly) and the later cases (September and October 1995) occur further north. This northward movement of the virus as the warm season progresses holds true for the 1997 and 1998 outbreaks.

Figure 3. Comparison of latitude and quarantine date

Vegetation

 

The incubation period for VSV is approximately two to eight days, or more (Fraser, 1991).  The NDVI value for each two-week time period that included a quarantine date for a case was extracted from the NDVI data for the Southwest for March-December 1995.  To cover the incubation period, the value for the two-week NDVI period immediately before the one containing in the quarantine date was also extracted.   A histogram of these two NDVI periods for the Rio Grande cluster is shown in Figure 4. Note that the NDVI values have been scaled; instead of –1 to 0 to +1, they have been scaled to 0 to 100 to 200.  The values from 100 to 200 represent the vegetated pixels.

 


Figure 4. Comparison of NDVI values for Rio Grande cluster, for NDVI period prior to quarantine date (left) and during quarantine date (right)




The histograms show a relatively normal distribution of NDVI values.  The cases within ± one standard deviation (68%) of the mean have a vegetative value in a fairly narrow band, 120-140 for the NDVI period prior to the quarantine date (left graph) and values 123-139 for the NDVI period covering the quarantine date (right graph).  Together, these two graphs describe the vegetation change during the incubation period for the VSV cases in the Rio Grande cluster.  They may represent a spectral signature for VSV risk.  Controls would make it possible to confirm whether these NDVI signatures apply only for diseased animals, or for all animals in the area.

 

Hydrology

Since water bodies provide habitats for possible vectors of VSV, the distance to water features is an important variable. This study considers moving or flowing water bodies only. The mean distance to water for the cases in the Rio Grande cluster is 200 meters. For the Cortez cluster the mean distance is 271 meters, and for the Grand Junction cluster it is 289 meters.  Without controls, it is difficult to discern if these distances are significant for diseased animals or is typical for all livestock in the area.  Figure 4 shows the distribution of distances for the cases in the Rio Grande cluster.

Figure 5. Distance in meters to nearest water feature for Rio Grande cluster

Moving water bodies can be subdivided into two types of features: streams or rivers and canals or ditches. Each type may support different vegetation and, thus, different vectors. Table 1 shows the types of water features for all cases in each cluster.

 

 

Feature type

Rio Grande (RG) cases

Percentage of RG cases closest to streams vs. canals

Cortez (CR) cases

Percentage of CR cases closest to streams vs. canals

Grand Junction (GJ) cases

Percentage of GJ cases closest to streams vs. canals

Stream or river

88 

50.9%

17 

80.9%

63 

60.6%

Canal or ditch

85 

49.1%

4 

19.1%

41 

39.4%

Totals  

173 

  

21 

  

104 

  

Table 1. Types of water features for cases in each cluster.

 

In the Cortez and Grand Junction areas, cases are more likely to occur close to a stream or river. The type of water feature did not appear to have an impact on case location in the Rio Grande area.

 

Terrain

Elevation was derived from DEM data.  Slope and aspect were then calculated from the elevation data using functions in ArcView's Spatial Analyst module.  Average elevation for the Rio Grande cluster is 1671 meters, for the Cortez cluster 1936 meters, and for the Grand Junction cluster 1545 meters. The average slope for the Rio Grande cluster is 3.1° with a range of 0° to 61.6°. For the Cortez cluster, the average slope is 4.7° with a range of 0° to 21.6°. For the Grand Junction cluster is 2.5° with a range of 0° to 30.6°.  For all three clusters the average aspect is between east and south-southeast. Again, it is difficult to tell is these attributes are specific to the diseased animals or typical to all animals in the area.

Species Distribution

Of the 173 cases in the Rio Grande cluster, 116 (66.5%) were equines, and the remainder were bovines and one camelid.  For the Cortez cluster, the composition of the 21 cases is 13 equines (61.9%), and the remainder bovines. For the 104 cases in the Grand Junction cluster, 96 (92.3%) were equines and the remainder bovines.

Livestock Density

The National Agriculture Statistics Service's 1992 Census of Agriculture provides data about livestock animal population on a county-by-county basis.  In 1995 the VSV cases did not occur in counties with the highest horse or cattle populations. It remains a possibility that horses and cattle are densely populated in local areas of a county where cases occurred.

Results and Conclusions

A preliminary landscape model of VSV-susceptible areas in the Southwest U.S. consists of several factors. Cases occur near moving water (within several hundred meters).  In the Southwest U.S., livestock are usually found in irrigated pastures.  The NDVI vegetation signature would include a peak in the 120-140 range shortly before the cases would appear.  The cases first appear with the start of the warm season and disappear with the advent of cold weather. The cases initially occur in the southern part of the affected region and progress northward, following the seasonal green-up of vegetation. It is not yet clear whether terrain features like elevation hold definitive clues to VSV outbreaks.

With regard to data collection, several conclusions can be drawn. The lack of controls makes it difficult to distinguish the environmental factors that indicate virus susceptibility from those that are associated with all livestock in the area.  Sufficient hardware and software and personnel trained in spatial data collection are needed to gather coordinate data as part of field visits during or shortly after an outbreak. There is also a need for standard protocols for spatial data collection.

Future Work

Future work will include additional data gathering and additional data analysis. This effort may include the acquisition of National GAP Analysis Program (A Geographic Approach to Planning for Biological Diversity) data on species distribution for potential hosts for the three clusters, and research into home ranges of possible vectors. Also, color infrared photography from NASA's High Altitude Photography Program will provide additional data on vegetation in the three clusters.  The Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) data set from the Environmental Protection Agency's Office of Water will provide stream flow rates for specific segments of the water features in each cluster for different times of year.

Existing and derived data will be used in statistical spatial analysis, possibly including regression analysis and surface modeling.  Analysis will look more closely at the cases within each cluster, to discern possible subclustering.  The 1998 outbreak will be used to refine and test ongoing results and conclusions. Coordinates were not collected for the 1997 outbreak, making those data less useful for spatial analysis.  The intended outcome is to provide a predictive capability for VSV based on environmental features to aid in disease prevention and management.

 

References

Fraser, Clarence M. et.al. editors 1991.  Merck Veterinary Manual, 7th ed. Rahway, N.J.: Merck & Co., Inc., 372pp.

USDA. 1996. Information about Vesicular Stomatitis for the Beef Producer  [On-line] Available: http://www.aphis.usda.gov/oa/pubs/fsvsbeef.html

 

Author Information

 

Susan Maroney

GIS Specialist

Veterinary Services, Centers for Epidemiology and Animal Health

Animal and Plant Health Inspection Service, U.S. Department of Agriculture

555 S. Howes Street

Fort Collins, CO 80526

voice (970) 490-7894

fax (970) 490-7999

susan.a.maroney@usda.gov

 

Jerome Freier

GIS Research Analyst

Veterinary Services, Centers for Epidemiology and Animal Health

Animal and Plant Health Inspection Service, U.S. Department of Agriculture

555 S. Howes Street

Fort Collins, CO 80526

voice (970) 490-7974

fax (970) 490-7999

jerome.e.freier@usda.gov