Implications of climate change on the distribution of the tick vector Ixodes scapularis and risk for Lyme disease in the Texas-Mexico transboundary region
Parasites & Vectors volume 7, Article number: 199 (2014)
Disease risk maps are important tools that help ascertain the likelihood of exposure to specific infectious agents. Understanding how climate change may affect the suitability of habitats for ticks will improve the accuracy of risk maps of tick-borne pathogen transmission in humans and domestic animal populations. Lyme disease (LD) is the most prevalent arthropod borne disease in the US and Europe. The bacterium Borrelia burgdorferi causes LD and it is transmitted to humans and other mammalian hosts through the bite of infected Ixodes ticks. LD risk maps in the transboundary region between the U.S. and Mexico are lacking. Moreover, none of the published studies that evaluated the effect of climate change in the spatial and temporal distribution of I. scapularis have focused on this region.
The area of study included Texas and a portion of northeast Mexico. This area is referred herein as the Texas-Mexico transboundary region. Tick samples were obtained from various vertebrate hosts in the region under study. Ticks identified as I. scapularis were processed to obtain DNA and to determine if they were infected with B. burgdorferi using PCR. A maximum entropy approach (MAXENT) was used to forecast the present and future (2050) distribution of B. burgdorferi-infected I. scapularis in the Texas-Mexico transboundary region by correlating geographic data with climatic variables.
Of the 1235 tick samples collected, 109 were identified as I. scapularis. Infection with B. burgdorferi was detected in 45% of the I. scapularis ticks collected. The model presented here indicates a wide distribution for I. scapularis, with higher probability of occurrence along the Gulf of Mexico coast. Results of the modeling approach applied predict that habitat suitable for the distribution of I. scapularis in the Texas-Mexico transboundary region will remain relatively stable until 2050.
The Texas-Mexico transboundary region appears to be part of a continuum in the pathogenic landscape of LD. Forecasting based on climate trends provides a tool to adapt strategies in the near future to mitigate the impact of LD related to its distribution and risk for transmission to human populations in the Mexico-US transboundary region.
Globally, Lyme disease (LD), also known as Lyme borreliosis is considered an emerging transboundary zoonotic tick borne disease[1–7]. Foci of LD exist in the U.S, Europe, and Asia. In the US LD is the most prevalent arthropod borne infection with 30,831 cases reported to the Centers for Disease Control and Prevention (CDC) in 2012 and it is estimated that over 300,000 people become infected every year. The increase in LD cases during the last decade has prompted its classification as an emerging infectious disease. Several hard tick species in the genus Ixodes are recognized as common vectors of the spirochete Borrelia burgdorferi, the causative agent of LD. Ixodes scapularis and I. pacificus are the known competent vectors in the US while I. persulcatus and I. ricinus are the documented vectors in Eurasia[3, 10–13]. The majority of confirmed LD cases in the US are reported from the Northeast and upper-Midwest regions. Consequently, almost all LD studies done to date are centered in the Northeast and Midwest regions of the US where this disease is more prevalent, with no major studies done in Southern US (SUS). Moreover, reports of LD in Mexico continue to increase[14–16]. In this geographic region I. scapularis is also considered as the main vector of B. burgdorferi, although the involvement of other hard tick species as competent vectors is suspected.
Similar to other arthropod borne diseases, LD is a complex system subject to shifts in ecological processes that influence vector biology and the epidemiology of B. burgdorferi infection in reservoir hosts and humans[18–22]. Molecular plasticity is a feature that allows the adaptation of B. burgdorferi to invertebrate and vertebrate hosts where it grows[23, 24]. The spirochete is maintained in the environment by different vertebrate hosts with varying degrees of competence. In the forests of eastern North America the white-footed mouse, Peromyscus leucopus, is its main reservoir[25, 26]. On the other hand, white-tailed deer (Odocoileus virginianus, WTD) are the primary reproductive host for I. scapularis in the US, but they are reservoir-incompetent for B. burgdorferi[19, 27]. A number of studies have suggested that vertebrate species diversity affects the risk of contracting human LD[28–32]. In addition, climate change has been predicted to affect the geographic distribution of tick vectors[33–39] and therefore the risk of LD infection. Species distribution models (SDMs) based on machine-learning algorithms and Geographic Information Systems (GIS) platforms have been used to predict areas of potential distribution of insects that serve as vectors of zoonotic infectious agents such as Trypanosoma cruzi, a parasite that causes Chagas disease[40, 41]. These analyses typically show that climatic factors significantly influence the potential geographic distributions of vectors and reservoirs. Additionally, temperature may have a strong influence on vector species behavior[42, 43]. For example, in the case of triatomine insects, temperatures higher than 30°C and low humidity increase triatomine feeding rates helping them to avoid dehydration. In addition, these insect vectors may develop shorter life cycles and higher population densities in domestic life cycles when indoor temperature increases. High temperatures can also speed up the development of T. cruzi in vectors. In regards to arthropod tick vectors, it has been observed that the spread of the tick vector I. scapularis in North American and Canada is followed by an increased invasion of the bacterial pathogen B. burgdorferi and therefore an increase in the risk of disease transmission. In addition, recent studies in Canada have shown that the distribution of I. scapularis has been significantly affected by climate, and in particular by temperature[46, 47]. Populations of this tick vector have been established and are spreading north. In this report, authors forecast, that with a continuous increase of temperature, the basic reproductive number (R 0 ) of I. scapularis in this region is predicted to increase, and therefore the risk for Lyme disease.
The limited number of ecological studies in other regions like the Southern US and Gulf Coast has prevented a comparative analysis of ecological factors driving the pathogenic landscape of LD in the US. In some parts of the world the ecology and epidemiology of LD remain to be fully understood. For instance, a LD-like syndrome has been described in Brazil, and ticks other than I. scapularis have been associated with the transmission of B. burgdorferi in Mexico. Thus, LD is considered a transboundary zoonotic disease in that it can reach epidemic proportions in different regions of the globe regardless of political borders. There is unequally distributed knowledge about the ecology of this disease among the regions in which it occurs, even within the US.
Knowledge gaps in our understanding of the epidemiology of LD in the Southern US exist. A recent survey could not detect B. burgdorferi from I. scapularis ticks collected from deer in Tennessee. During the late 1980’s and early 1990’s B. burgdorferi was identified in Texas using classical culturing techniques and immunofluorescence of tick midguts utilizing monoclonal anti-OspA antibody with an incidence as low as 1%[52–55]. With more advanced molecular techniques the incidence of B. burgdorferi in ticks collected from human subjects over a 4-year period was also near 1%. However, this last study was based on a passive surveillance set up, in which ticks coming from humans were collected, and no active sampling of ticks in the state of Texas was carried out.
Human risk for infection with B. burgdorferi over the continental US (east of the 100th meridian) has been predicted using the density of I. scapularis infected nymphs (DIN)[56, 57]. Under this scenario southern US states were portrayed as a low risk region given the non-appearance of host seeking I. scapularis nymphs at sampled sites[56, 57]. Worth noting is that southern states had fewer sampling sites (e.g., 9 in Texas, zero in Louisiana, 4 in Oklahoma) while other areas of the country (the northeast and midwest) were better represented (e.g. 19 in Illinois, 20 in Wisconsin, 13 in New York)[56, 57]. In striking contrast to the conclusion of this purported null risk of acquiring LD in southern states, a steady number of LD cases have been reported in these low risk areas every year (http://www.cdc.gov/lyme/stats/index.html). Other relevant epidemiological aspects that require attention include human movement because in some LD cases the disease can be acquired in a region different from the one where they are reported and differences in tick phenology and behavior at different geographic regions may exist. For instance, the questing behavior of I. scapularis nymphs show variation among northern and southern sites, with southern nymphs rarely climbing to the top of vegetation likely reducing the ability to detect these nymphs at southern sites (J. Tsao, personal communication). These limitations might explain why the models utilized could not explain the variation in the distribution of the disease observed in low prevalence areas. Consequently, the drivers for variation in the distribution of LD cases observed in low prevalence areas remain to be identified.
In Mexico, a national serosurvey of human serum samples reported a B. burgdorferi sero-prevalence of 1.1%. The Mexican states of Tamaulipas, Nuevo León, and Coahuila in the Texas-Mexico border region presented the highest sero-prevalence (6.4%) compared with the rest of the country. Also, Ixodes ticks infected with B. burgdorferi sensu stricto occur in the same states, and the infection has been recently documented in mouse reservoirs as well as in WTD. Distribution models of potential tick vectors in Mexico point to a wide distribution range that overlaps not only northeastern Mexican states along the border with the US, but also extends to central Mexico[59, 61]. These studies, together with confirmed clinical cases of LD, acquired in parks near Mexico City[14, 62] demonstrate the existence of a enzootic cycle responsible for LD in Mexico.
The large transboundary region shared between the US and Mexico present several risk factors for zoonotic arthropod borne diseases that can be exacerbated due to climate change[6, 63, 64]. This situation emphasizes the need to evaluate the risk for LD from a regional standpoint. In addition, the understanding of vector borne disease ecology has improved significantly in recent years due to the advancement in molecular biology, Geographic Information Systems (GIS) and SDM[65, 66].
In this study a distribution model was constructed for I. scapularis covering the eastern sector of the transboundary region comprising the states of Texas in the US, and Tamaulipas, Nuevo León, and Coahuila in Mexico. Future distribution models for the tick vector were forecasted under climate change scenarios within the defined study area. Additionally, data on the prevalence of B. burgdorferi infection in I. scapularis sampled at sites in Texas and the neighboring northeastern Mexican states are presented. This international research collaboration affords the opportunity to foster a binational approach for the adoption of control and prevention strategies to manage the risk, and reduce the burden of LD in human populations living along the US-Mexico border.
The Texas-Mexico transboundary region considered here includes sectors on both sides of the Rio Grande river, east of the Chihuahuan desert, which share a common biogeographic history and represents distinct ecozones. This region encompasses the Tamaulipan scrub biome. Based on the updated Köppen–Geiger climate classification[67, 68], this biome climate includes four categories. The majority of the biome is classified as an arid-steppe (BSh), which changes to a temperate classification (Cfa) in Texas, whereas in Mexico it includes small areas with climate generally classified as arid desert (BWh), or temperate (Cwa). The region consists of four ecoregions including Gulf Prairies and Marshes, South Texas Plains (also known as Rio Grande Plains), Edwards Plateau (northeast of the river), and the Tamaulipan brushlands (southwest of the river). Elevation ranges from sea-level to around 1000 feet near the most inland locations. Annual rainfall may be 30-50 inches along the Gulf Coast, but declines to 14-16 inches in the inland reaches of the South Texas Plains-Edwards Plateau. Along the Gulf Coast summers are long, hot and humid, i.e. subtropical, and winters are very mild (>300 frost free days), while progressively more inland locations have hotter, drier summers, and winters subject to more freezing temperatures. The region is subject to the fluctuations of tropical storms from the Caribbean and Gulf of Mexico, as well as periodic and sometimes extended periods of drought. Within the region considered in our study we obtained and analyzed a total of 109 I. scapularis ticks from 16 unique geographical locations (Figure 1).
Tick collection and classification
In Texas, ticks were collected during 2011 and 2012 from three Texas Parks and Wildlife Management Areas (WMAs) covering five counties (Anderson, Cameron, Hidalgo, Mason and Willacy), animal shelters (Brazos County), veterinary clinics from 13 counties (Bee, Collin, Edwards, Fort Bent, Harris, Kerr, Maverick, Montague, Montgomery, Robertson, Tarrant, Travis, and Uvalde), and the Big Thicket National Park in Tyler County. Counties were selected based on the presence of I. scapularis ticks previously reported in the literature[13, 19, 54–56, 73]. The main sources of ticks were hunter-harvested WTD (Odocoileus virginianus) or gemsbok (Oryx gazella) brought to check stations at each WMA mentioned above. Questing ticks as well as tick samples from humans and dogs were also acquired from WMAs in Texas. Ticks collected from human subjects and their dogs at the hunter-check stations were voluntarily provided and no personal identifiers were recorded. Only location, date and the fact that the ticks were feeding on human subjects were recorded. All human subjects were properly informed about the scientific nature of the project.
For each tick specimen we gathered associated information (host species, time and date, GPS coordinates, ZIP code, County) to allow for cross-referencing our findings with the reported human LD cases in Texas. Results from those analyses will be presented somewhere else (manuscript in preparation). All ticks collected were placed in vials containing 70% ethanol and transported to the Texas A&M University College of Veterinary Medicine and Biomedical Sciences for their codification, identification using dichotomous keys[74–77], and subsequently processed for DNA extraction and B. burgdorferi detection.
In Mexico, I. scapularis ticks were collected in the northeastern states of Tamaulipas, Nuevo Leon, and Coahuila. Tick specimens were obtained from four different regions where vegetation was sampled for questing ticks and from collected wild mammals (eastern cottontail, the hispid cotton rat, the painted spiny pocket mouse, prairie dogs, ringtail, black-tailed jackrabbit and jaguar) and domestic animals (bovids, horses, and dogs)[17, 78]. Ticks were preserved in 70% ethanol and archived. Two institutions were involved in the collections, the Entomology Laboratory of Escuela Nacional de Ciencias Biologicas, IPN, Mexico, which collected in forested and rural areas during 1990-1998, and the Emerging Infectious Disease Laboratory, UIMEIP, Centro Medico Nacional SXXI, IMSS, which collected ticks in forested areas during 1999-2007[17, 78]. An additional record of I. scapularis was obtained from the report of this species occurring on a female jaguar (Panthera onca) from the Reserva de la Biosfera, El Cielo in Tamaulipas Mexico. All wild mammal collections were done under the approved permits by Dirección General de Vida Silvestre in Mexico.
Borrelia burgdorferi identification
DNA was extracted from individual I. scapularis ticks collected during the sampling period described above using the High Pure PCR template preparation kit for genomic DNA (Roche, Indianapolis, IN) following manufacturer’s recommendations with modifications to adapt the protocol to ticks. Briefly, individual ticks were homogenized utilizing the bead mill BeadRuptor 24 (Omni International, Inc., Kennesaw, GA) in 200 μl of Phosphate Buffer Saline (PBS) plus 200 μl of lysis buffer and 800 μg of Proteinase K in 2 ml screw cap tubes for 5 minutes at a 5.65 m/s intensity (equivalent to 210 × g). During this step 1.4 mm ceramic beads were used to homogenize unengorged ticks while 2.8 mm ceramic beads were used with engorged ticks. After homogenizing, each sample was centrifuged at 10,000 × g to pellet debris. Supernatants were collected and 200 μl of binding buffer was added to each sample and the mixture incubated at 70°C for 10 minutes followed by the addition of 100 μl of isopropanol. The mixture was run through filter columns at 8,000 × g for 1 minute. DNA bound to the filter was washed and eluted following manufacturer recommendations. Eluted DNA concentration and purity were measured using a NanoDrop (Thermo Scientific, Inc., Willinton, DE). DNA was used for the detection of B. burgdorferi sensu lato by PCR using the flagella gene (flaB) as previously described by Jaulhac et al.. Amplicons were separated on a 1.5% agarose gel and positive bands were sent for sequencing (Eton Biosciences, San Diego, CA). Sequences were analyzed through BLAST® using MacVector 12.6 software (MacVector Inc., Cary, NC) to facilitate species and strain identification. A second nested PCR for the intergenic region (IGR) 16SrRNA-23SrRNA using the primers and conditions described by Bunikis et al. was done to all flaB positive samples with B. burgdorferi identification. IGR amplicons were also submitted for sequencing and results analyzed using the MacVector 12.6 software to facilitate strain identification. Water served as negative control, and genomic DNA isolated from B. burgdorferi B31 A3 strain culture as positive control in all PCR steps done in this study.
Modeling methods and evaluation
To construct the potential distribution for I. scapularis at the study area we compiled a database with unique locations (e.g., different geographic coordinates) of presence records that cover most of the known range of this species. Locations were obtained from the scientific literature and from our own fieldwork (Texas and Northern Mexico). Our modelling methods required identification of precise locations where this tick species was found but a prevalent issue regarding tick records is that most lack an adequate level of geographical precision or reliable taxonomic identification. For example, the most extensive collection of reliable records for I. scapularis, covering most of its range, is limited to county level precision and thus it is too rough for our purposes. Instead, we accessed data from several published studies that included standardized collecting techniques to sample 304 sites throughout 37 states from the continental US east of the 100th parallel. This study consisted mostly of state parks or natural preserves[56, 57, 81, 82]. I. scapularis was found at 94 of these sites. This was the basis for our locations. For each park or natural area where I. scapularis was recorded we identified a geographical centroid that was verified to occur within natural vegetation (via Google Earth) and the corresponding geographic coordinates were used as a unique location. Since most of the parks are small (<10,000 Ha), the level of resolution is substantially better than relying on county level records.
Present and projected future potential distributions for the target species were computed using these presence records for the species (longitude/latitude) and using climatic parameters as exploratory variables, using a maximum entropy algorithm (MaxEnt)[66, 83–85]. MaxEnt predicts a probability distribution ranging from 0 (least suitable) to 1 (most suitable) of habitat suitability over a given study area[84, 86]. We used MaxEnt version 3.3.3 k (http://www.cs.princeton.edu/~schapire/maxent/) with the default modeling parameters (convergence threshold = 105, maximum iterations = 500, regularization value β = auto) following conditions standardized elsewhere. We used the bioclimatic variables available athttp://www.worldclim.org as predictors (Table 1).
These bioclimatic variables have a resolution of approximately 1 × 1 km2 grid cells. To avoid the use of highly correlated variables, we ran a correlation analysis and eliminated one of the variables in each pair with a Pearson correlation value > 0.8 (Additional file1: Table S1). Sixteen variables were used to construct the model (Table 1). To avoid the use of artifactual distribution limits for our target species (e.g., political boundaries), we used as base map to construct our models a biogeographic regions map (Figure 2) and then projected the models to Texas and Northeastern Mexico (Tamaulipas, Nuevo Leon, Coahuila). Biogeographic regions are areas of animal and plant distribution having similar or shared characteristics, thus this map is considered essential to construct a more realistic model of potential distribution for I. scapularis. Model results were processed and visualized using ArcGIS 10. The Area Under the Curve (AUC) of Receiver Operating Characteristics plots (ROC) was calculated to evaluate the models. ROC is a threshold–independent measure that evaluates the sensitivity (probability that the model produces a positive result in a positive locality) versus the specificity (probability that the model produces a negative result in a negative locality) of a model when presented with new data. A ROC plot is obtained by plotting all sensitivity values on the y–axis against their equivalent (1–specificity) values for all available decision thresholds on the x–axis. The theoretically perfect result is AUC = 1, whereas a test performing no better than random yields AUC = 0.5. The AUC was calculated in MaxEnt.
Projection under future scenarios
Future distribution models were developed for the year 2050 using two climate scenarios from the Intergovernmental Panel on Climate Change (IPCC): A2A and B2A. These scenarios assume a heterogeneous world, but with B2A being more environmentally focused than A2A. B2 describes a world in which the emphasis is on local solutions to economic, social, and environmental sustainability. It is a prospective view of the world with continuously increasing global population, at a rate lower than A2, and intermediate levels of economic development. Different Global Circulation Models (GCM) exist, and since no single one can be considered superior, we selected three different GCM to compare results for potential distributions: the Canadian Centre for Climate Modeling and Analysis (CCCMA), the Commonwealth Scientific and Industrial Research Organization (CSIRO) and the Hadley Centre for Climate Change (HADCM3) models.
We estimated a “stable” suitable habitat considering present and future maps (both IPCC scenarios and the three GCM). First, we reclassified the potential distribution models for the present and future to obtain binary maps (presence = 1, absence = 0) using the “minimum presence threshold”. This is the most conservative threshold that includes suitable habitat predictions for all the locations where we have collected I. scapularis ticks. The binary maps were then projected onto UTM with a spatial resolution of 1 km2. The “stable” habitat was estimated as the area that remained constant under each scenario and GCM.
During this study, a total of 1235 tick specimens were collected in the Texas-México transboundary region. One hundred and nine specimens (8.83%) were identified as I. scapularis, which were collected questing or from WTD, dogs, gemsbok, and a cat in Texas, and eastern cottontail (Silvilagus floridanus), painted spiny pocket mouse (Lyomis pictus) and jaguar (Pantera Onca) in Mexico, as shown in Table 2. Forty five percent of the I. scapularis were infected with B. burgdorferi by PCR and sequencing of the flagella and IGR amplicons.
Presence of I. scapularis and B. burgdorferi infection in the Texas-Mexico transboundary region
A total of 574 tick specimens were collected from twenty counties across the state of Texas. Of those, 74 (12.89%) I. scapularis ticks were collected in eight of the counties sampled (Table 2). Fifty percent of I. scapularis ticks showed infection with B. burgdorferi sensu stricto (Table 2). This accounts for 6.445% of the total number of ticks tested in Texas in the current study. Even though the majority of the B. burgdorferi infected I. scapularis ticks were found on the most Eastern counties in Texas (Anderson and Tyler Counties), they were also detected in Central (Mason County), and South Texas (Cameron County, located on the border with Mexico). The central and south Texas areas had been previously considered as unsuitable for I. scapularis[56, 65, 89]. Only one location in Texas (Hidalgo, Table 2) showed no positive I. scapularis ticks for B. burgdorferi. When analyzing the intergenic region 16S-23S rRNA of B. burgdorferi detected in the ticks collected in Texas, the sequences tend to group in two main clusters (Figure 3). These results suggest that strains similar to the positive control used in this study, B. burgdorferi B31 strain, were rarely detected in the area of study.
Of the 661 ticks collected from 11 localities in México, 203 were identified as Ixodes sp and only 35 as I. scapularis. Twelve I. scapularis were infected with B. burgdorferi and were collected in six of the 7 localities (Table 2). The infected Ixodes ticks were predominantly from four rural localities in San Josesito district, state of Nuevo Leon, and only one locality in Tampico district, state of Tamaulipas. A single I. scapularis tick was collected in Coahuila, which was not infected with B. burgdorferi.
Present and future suitable habitat for I. scapularis in the Texas-Mexico region
One hundred and ten unique locations, at which I. scapularis is already present, were used to develop models of present and future suitable habitats for I. scapularis (Figure 2). From those, 94 came from the literature and 16 from our fieldwork: 9 were from Texas, and 7 were from NE Mexico states along the Texas-Mexico border (Table 2). The model of potential distribution for I. scapularis had an AUC =0.831. Comparison of present (Figure 4) and future (Figure 5A-F) suitable habitat shows potential shifts in the distribution of suitable habitat for I. scapularis due to climate change. As shown in Figure 4, the area of I. scapularis distribution is expected to expand in Northeastern Texas as well as along the coastal areas, whereas in Mexico the model predicts maintenance of the distribution area in the three northeastern states. Interestingly, there is a continuous line along the gulf coast that will remain a highly suitable habitat for this tick vector in all of the forecasting maps generated (Figure 5A-F).
The climatic variables that contributed the most to the model are presented in Table 3. Briefly, isothermality accounts for 20% together with precipitation of the wettest quarter (18%), the maximum temperature of the warmest month (14.6%), and the precipitation observed in the wettest month (11.5%). These parameters account for 64.2% of the average model. The correlation between each one of these four relevant environmental variables and the probability for the presence of I. scapularis confirm the close tie of this species to precipitation and temperature. Considering these environmental parameters, the suitable area for the survival of I. scapularis that will be maintained stable regardless of the forecasting model used will cover an area of 569,910 Km2 along Eastern Texas all the way to Northern Mexico and along the Gulf Coast (Figure 6).
As an emerging transboundary zoonotic tick-borne disease, LD affects thousands of humans and domestic animals around the world[5, 7, 90]. Transboundary zoonotic diseases can maintain a dynamic focus and have pathogens circulating in a geographic region encompassing multiple geopolitical boundaries. The sector of the US-Mexico border shared by south Texas and the states of Tamaulipas, Nuevo León, and Coahuila, is inhabited by millions of people and includes a biogeographic region linking the Neotropic and Nearctic ecozones where foci of vector-borne diseases like LD exist[91–94]. The distribution map presented here documents the transboundary aspect for the risk of contracting LD in the region encompassing Texas and northeast Mexico by taking into account the presence of the tick vector I. scapularis infected with B. burgdorferi.
There are a number of studies evaluating the distribution of I. scapularis, the competent vector for the transmission of LD[28, 47, 51, 56, 57, 82, 95–98], and how its distribution affects the risk for human disease. The analysis in most of these studies is focused in the Northeastern and Midwest US as well as in Canada. A climate based model classified each of 3,628 0.5° grid cells, covering the US, as places where established, suitable and unsuitable conditions to support I. scapularis exist, thus providing an estimation of its geographic distribution in the US. This model portrays most of East Texas as suitable habitats or where established populations can exist. Our current species distribution model of I. scapularis overlaps greatly over East Texas with the climate-based model but also expands the range of habitat identified as suitable for I. scapularis towards Central and South Texas through the Gulf Coast. Additionally, our model expands into northern Mexico as a continuous geographic landmark to maintain a stable I. scapularis population at this border region. Visualization of the predicted suitable area depicted in Figure 5 is further evidence that the Texas-Mexico transboundary region is part of a continuum in the pathogenic landscape of LD serving as a corridor linking B. burgdorferi with tick vectors and mammalian reservoirs in the Neotropic and Nearctic ecozones. It remains to be determined if migratory birds are the only animals dispersing infected I. scapularis in the region, which could establish LD foci. More detailed molecular studies of Borrelia spp. identified in I. scapularis and other ticks will enhance our understating of the diversity among species in the B. burgdorferi sensu lato complex and the tick vectors involved in the epidemiology of Lyme borreliosis in the Texas-Mexico transboundary region[5, 61, 100].
No specific distribution model exists for I. scapularis in Mexico but a distribution model for the genus Ixodes, generated with similar methodologies to our model, also predicts a wide distribution covering northeastern Mexico. Although our model also overlaps greatly with this genus level model, especially over the mountainous regions of Tamaulipas, Nuevo Leon, and Coahuila, it expands the distribution of Ixodes over the coastal plain of Tamaulipas thus connecting to the range in Texas along the Gulf Coast area. The level of resolution of our model is substantially greater than in previous climate based models[89, 101] because of the smaller grid cell used (1 Km × 1 Km cells in our model versus roughly 55 × 42 Km cells), thus allowing a more precise prediction of sites where I. scapularis currently occurs. To our knowledge, this report is the first attempt to model how the range of I. scapularis in the Texas-Mexico transboundary region will be affected under specific climatic change scenarios. Only one location in Texas showed no positive I. scapularis ticks for B. burgdorferi. We attribute this to the low number of ticks (only 7) collected with only one I. scapularis. Geographical areas in which a species occur are determined by several ecological, evolutionary, and historical factors. These factors can be related to (1) abiotic aspects (e.g., elevation, climate); 2) biotic aspects (i.e., species interactions like competition and parasitism); 3) historical characteristics (e.g., barriers, speciation process); 4) dispersal capabilities; 5) accessible regions for dispersal; 6) evolutionary capability of a species to adapt to new conditions; and 7) anthropogenic influences (e.g., changes in land use, translocation of organisms;[102, 103]). The most frequently used factors to estimate species distributions are climatic variables[104, 105]. Climate can limit the distribution of animals directly by affecting growth or survival (e.g., lower and upper lethal temperatures), and indirectly via interacting species (e.g., food sources, pathogens, competitors, or predators). While testing the model, two locations in the region under study, had low values of suitable habitat predicted by MaxEnt for I. scapularis. These two localities are in West Texas (Mason County) and in the Mexican state of Coahuila. These are two of the most western locations, sampled, where I. scapularis was detected in vegetation and feeding on exotic ungulates inhabiting gated properties. Consequently, these locations are considered of great biological relevance, because they could represent the western boundary of the geographic region suitable for I. scapularis. This observation highlights the need for ecological studies to understand the phenology and host preference in locations marginally suitable for the blacklegged tick.
Although our species distribution model is congruent with the spatial distribution of LD cases in humans for Texas (see maps inhttp://www.cdc.gov/lyme/stats/index.html), a direct correlation cannot be established because additional vector density and B. burgdorferi prevalence data need to be gathered, and human movement has epidemiological consequences since infection can be acquired in a region different from the one where it is reported. The density of infected I. scapularis nymphs (DIN) has been suggested as a strong predictor for LD risk of infection[56, 57, 82, 89, 101]. For Texas this risk of infection has been deemed low because no I. scapularis nymphs were detected using standardized sampling techniques over nine localities in the state[56, 82, 106]. This contrasts with the continuous reports of human LD cases and the detection of canine LD in the state[107–109]. As an initial approach to understand where infected I. scapularis ticks occurred in Texas and Mexico, we performed a combination of passive and active searches of different developmental stages of this competent LD vector. In our search we did find B. burgdorferi infected adult I. scapularis ticks in areas were human LD has been reported and where canine LD has been observed. The prevalence of I. scapularis infected with B. burgdorferi detected was higher than the 1% recorded in the early 1990’s using classical culture and microscopy techniques for the detection of B. burgdorferi[52, 54, 55], and than that observed more recently through a passive surveillance study and molecular techniques. The adult ticks found in this preliminary study were either questing on vegetation or feeding on wildlife (WTD and gemsbok) and companion animals (mainly dogs). Therefore, our results suggest that infected I. scapularis nymphs are present in Texas, but sampling methods standardized for the conditions in Northeast and Midwest US sites might not be optimal to sample a representative I. scapularis nymph population questing in the Texas-Mexico transboundary region. Variation in environmental conditions influencing questing behavior may significantly impact the type of host ticks encounter and may lead to differential host use within a particular study area[110–113]. Therefore, further studies testing different sampling procedures, including different times of the day and seasons are currently being developed by our team. Our goal is to determine the phenology of I. scapularis in the Southern US and Mexico and the questing behavior of the different developmental stages to determine how that will affect the risk for LD in humans and companion animals.
In forested regions of the northeastern US a strong determinant of nymph density, and percentage of Borrelia infected nymphs, is the composition of the mammalian assemblage present at a given area[25, 26]. At locales where the main competent reservoir, P. leucopus, is the dominant component of the assemblage, a higher prevalence of Borrelia infected nymphs is present, and that prevalence diminishes as other small mammal species become more common[25, 26]. P. leucopus has a wide distribution range from Canada to south Mexico and although it is the most common small mammal in mixed forests (deciduous and coniferous) at eastern states it can be more restricted in the western and southern parts of its range. This rodent uses diverse habitats and the plant communities it inhabits in the northeastern US are strikingly different from the more xeric habitats in the south, to the point that it is rarely found in forests but is more common in highly disturbed habitats in Mexico. Given these differences of habitat use between northern and southern localities, the relative abundance of P. leucopus may vary between different regions. Assessing host range for I. scapularis nymphs and the reservoir competence for B. burgdorferi of other small mammals is needed to understand the pathogenic landscape of LD in the Texas-Mexico transboundary region.
The tick collection effort reported here was initiated in Texas during a historical drought in the Southern US (http://droughtmonitor.unl.edu), with record high temperatures and record low precipitations affecting 100% of the Texas with “abnormally dry” to “severely dry” conditions that continued in some parts of the state until 2013. Even under these conditions we were able to collect I. scapularis ticks. This observation suggests that there were areas where microhabitat conditions allowed blacklegged ticks to overcome high temperatures and desiccation. Under this scenario continuous sampling of our area of study will be essential to determine where and how the immature stages of I. scapularis survive. In addition, continuous sampling of more localities will help validate and refine the generated maps that forecast the current and future distribution of this tick vector in the Texas-México transboundary region. However, public safety issues in northern Mexico have prevented further fieldwork. Current efforts are focused on additional sampling of Texas locations utilizing different sampling techniques to enhance detection of I. scapularis nymphs. Ongoing samplings of small mammal assemblages in Texas will help us understand which are the vertebrate hosts for the immature stages of I. scapularis in Southern latitudes.
In our analyses isothermality, maximum temperature of the warmest month, precipitation of the wettest month, and the precipitation of the wettest quarter were the environmental factors that contributed the most to the model. The correlation between each one of these four relevant environmental variables and the probability for the presence of I. scapularis confirm the close tie of this species to precipitation and temperature[58, 62]. These results correspond with the fact that I. scapularis survival depends on the temperature and humidity of a particular region[18, 36, 115]. Moreover, average precipitation in a certain area can explain the fluctuation of LD cases reported to CDC. The total number of LD cases reported appears to follow the 2-year tick density pattern, in which a year with more rain and higher number of ticks are followed by a peak in LD cases two years later. Taken together, our results suggest that similar parameters to those observed in other parts of the country govern the distribution of I. scapularis in Texas. However, attention is required to adapt the timing and methodologies for tick sampling campaigns to optimize the collection of the different developmental stages of this tick species. Several combinations of climatic variables were explored to generate our current model, but no major changes in the distribution of infected I. scapularis were observed over the region of interest. More research is required to ascertain similarities and differences in the ecology of LD between the northeast and upper-midwest parts of the US and the Texas-Mexico transboundary region and their effect on the risk of exposure to B. burgdorferi among human and domestic animal populations. These efforts are expected to enhance the robustness of the I. scapularis distribution model presented here. In combination with genetic data, consideration of the entire black legged tick range in a refined version of our model would allow testing if genetic variation of I. scapularis is structured according to geography, to host use, or both. Geographic and host-based structure has been documented for I. ricinus and other Ixodes species[111, 112, 116–118].
The model presented here shows that, the distribution range of I. scapularis, the competent vector for the transmission of LD can potentially be widely occurring in Eastern Texas and Northern parts of Mexico along the Gulf Coast as a continuous region. Extreme climate change conditions predict the maintenance of suitable areas for I. scapularis similar to those currently observed and potentially expanding towards central Texas. The ecology of LD remains to be fully understood in the Texas-Mexico transboundary region. Under these conditions, it will be necessary to continue the study of the distribution of this tick vector, at the western and southern edges of its known range. Consequently, research on the range of hosts the larvae and nymphs of I. scapularis can parasitize is required to enhance our understanding of the risk for LD transmission in the Southern US and México. The enzootic cycle of B. burgdorferi in this region remains to be fully understood. A binational and collaborative approach to such scientific efforts is recommended to develop and adapt control and prevention strategies that will help manage the risk, and reduce the burden of LD in human populations living along the US-Mexico border.
Parola P, Raoult D: Ticks and tickborne bacterial diseases in humans: an emerging infectious threat. Clin Infect Dis. 2001, 32 (6): 897-928. 10.1086/319347.
Stanek G, Reiter M: The expanding Lyme Borrelia complex–clinical significance of genomic species?. Clin Microbiol Infect. 2011, 17 (4): 487-493. 10.1111/j.1469-0691.2011.03492.x.
Stanek G, Wormser GP, Gray J, Strle F: Lyme borreliosis. Lancet. 2012, 379 (9814): 461-473. 10.1016/S0140-6736(11)60103-7.
Margos G, Vollmer SA, Ogden NH, Fish D: Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato. Infect Gen Evol. 2011, 11 (7): 1545-1563. 10.1016/j.meegid.2011.07.022.
Rudenko N, Golovchenko M, Grubhoffer L, Oliver JH: Updates on Borrelia burgdorferi sensu lato complex with respect to public health. Ticks Tick Borne Dis. 2011, 2 (3): 123-128. 10.1016/j.ttbdis.2011.04.002.
Siembieda JL, Kock RA, McCracken TA, Newman SH: The role of wildlife in transboundary animal diseases. Anim Health Res Rev. 2011, 12 (1): 95-111. 10.1017/S1466252311000041.
Dubrey SW, Bhatia A, Woodham S, Rakowicz W: Lyme disease in the United Kingdom. Postgrad Med J. 2014, 90 (1059): 33-42. 10.1136/postgradmedj-2012-131522.
CDC: Notice to readers: final 2012 reports of nationally notifiable infectious diseases. MMWR Morb Mort Wky Rep. 2013, 62 (33): 669-682.
Kuehn BM: CDC estimates 300,000 US cases of Lyme disease annually. JAMA. 2013, 310 (11): 1110-10.1001/jama.2013.278331.
Mannelli A, Bertolotti L, Gern L, Gray J: Ecology of Borrelia burgdorferi sensu lato in Europe: transmission dynamics in multi-host systems, influence of molecular processes and effects of climate change. FEMS Microbiol Rev. 2012, 36 (4): 837-861. 10.1111/j.1574-6976.2011.00312.x.
Radolf JD, Caimano MJ, Stevenson B, Hu LT: Of ticks, mice and men: understanding the dual-host lifestyle of Lyme disease spirochaetes. Nat Rev Microbiol. 2012, 10 (2): 87-99.
Jaulhac B, Heller R, Limbach FX, Hansmann Y, Lipsker D, Monteil H, Sibilia J, Piemont Y: Direct molecular typing of Borrelia burgdorferi sensu lato species in synovial samples from patients with lyme arthritis. J Clin Microbiol. 2000, 38 (5): 1895-1900.
Williamson PC, Billingsley PM, Teltow GJ, Seals JP, Turnbough MA, Atkinson SF: Borrelia, Ehrlichia, and Rickettsia spp. in ticks removed from persons, Texas, USA. Emerg Infect Dis. 2010, 16 (3): 441-446. 10.3201/eid1603.091333.
Gordillo-Perez G, Torres J, Solorzano-Santos F, de Martino S, Lipsker D, Velazquez E, Ramon G, Onofre M, Jaulhac B: Borrelia burgdorferi infection and cutaneous Lyme disease, Mexico. Emerg Infect Dis. 2007, 13 (10): 1556-1558. 10.3201/eid1310.060630.
Gordillo G, Torres J, Solorzano F, Cedillo-Rivera R, Tapia-Conyer R, Munoz O: Serologic evidences suggesting the presence of Borrelia burgdorferi infection in Mexico. Arch Med Res. 1999, 30 (1): 64-68. 10.1016/S0188-0128(98)00015-3.
Gordillo-Perez G, Torres J, Solorzano-Santos F, Garduno-Bautista V, Tapia-Conyer R, Munoz O: Seroepidemiologic study of Lyme’s borreliosis in Mexico City and the northeast of the Mexican Republic. Salud Publica Mex. 2003, 45 (5): 351-355. 10.1590/S0036-36342003000500004.
Gordillo-Perez G, Vargas M, Solorzano-Santos F, Rivera A, Polaco OJ, Alvarado L, Munoz O, Torres J: Demonstration of Borrelia burgdorferi sensu stricto infection in ticks from the northeast of Mexico. Clin Microbiol Infect. 2009, 15 (5): 496-498. 10.1111/j.1469-0691.2009.02776.x.
Ostfeld RS: Lyme disease: the ecology of a complex system. 2011, New York: Oxford University Press
Hoen AG, Margos G, Bent SJ, Diuk-Wasser MA, Barbour A, Kurtenbach K, Fish D: Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events. Proc Natl Acad Sci U S A. 2009, 106 (35): 15013-15018. 10.1073/pnas.0903810106.
Kurtenbach K, Hanincova K, Tsao JI, Margos G, Fish D, Ogden NH: Fundamental processes in the evolutionary ecology of Lyme borreliosis. Nat Rev Microbiol. 2006, 4 (9): 660-669. 10.1038/nrmicro1475.
Kurtenbach K, De Michelis S, Sewell HS, Etti S, Schafer SM, Holmes E, Hails R, Collares-Pereira M, Santos-Reis M, Hanincova K, Labuda M, Bormane A, Donaghy M: The key roles of selection and migration in the ecology of Lyme borreliosis. Int J Med Microbiol. 2002, 291 (Suppl 33): 152-154.
Levi T, Kilpatrick AM, Mangel M, Wilmers CC: Deer, predators, and the emergence of Lyme disease. Proc Natl Acad Sci U S A. 2012, 109 (27): 10942-10947. 10.1073/pnas.1204536109.
Bunikis J, Garpmo U, Tsao J, Berglund J, Fish D, Barbour AG: Sequence typing reveals extensive strain diversity of the Lyme borreliosis agents Borrelia burgdorferi in North America and Borrelia afzelii in Europe. Microbiol. 2004, 150 (Pt 6): 1741-1755.
Tsao JI: Reviewing molecular adaptations of Lyme borreliosis spirochetes in the context of reproductive fitness in natural transmission cycles. Vet Res. 2009, 40 (2): 36-10.1051/vetres/2009019.
LoGiudice K, Duerr ST, Newhouse MJ, Schmidt KA, Killilea ME, Ostfeld RS: Impact of host community composition on Lyme disease risk. Ecology. 2008, 89 (10): 2841-2849. 10.1890/07-1047.1.
LoGiudice K, Ostfeld RS, Schmidt KA, Keesing F: The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc Natl Acad Sci U S A. 2003, 100 (2): 567-571. 10.1073/pnas.0233733100.
Rosen ME, Hamer SA, Gerhardt RR, Jones CJ, Muller LI, Scott MC, Hickling GJ: Borrelia burgdorferi not detected in widespread Ixodes scapularis (Acari: Ixodidae) collected from white-tailed deer in Tennessee. J Med Entomol. 2012, 49 (6): 1473-1480. 10.1603/ME11255.
Ogden NH, Tsao JI: Biodiversity and Lyme disease: dilution or amplification?. Epidemics. 2009, 1 (3): 196-206. 10.1016/j.epidem.2009.06.002.
Kurtenbach K, Sewell HS, Ogden NH, Randolph SE, Nuttall PA: Serum complement sensitivity as a key factor in Lyme disease ecology. Infect Immun. 1998, 66 (3): 1248-1251.
Brisson D, Brinkley C, Humphrey PT, Kemps BD, Ostfeld RS: It takes a community to raise the prevalence of a zoonotic pathogen. Interdiscip Perspect Infect Dis. 2011, doi:10.1155/2011/741406
Brisson D, Dykhuizen DE, Ostfeld RS: Conspicuous impacts of inconspicuous hosts on the Lyme disease epidemic. Proc Biol Sci. 2008, 275 (1631): 227-235. 10.1098/rspb.2007.1208.
Brisson D, Vandermause MF, Meece JK, Reed KD, Dykhuizen DE: Evolution of northeastern and midwestern Borrelia burgdorferi, United States. Emerg Infect Dis. 2010, 16 (6): 911-917. 10.3201/eid1606.090329.
Ogden NH, Bigras-Poulin M, Hanincová K, Maarouf A, O’Callaghan CJ, Kurtenbach K: Projected effects of climate change on tick phenology and fitness of pathogens transmitted by the North American tick Ixodes scapularis. J Theor Biol. 2008, 254 (3): 621-632. 10.1016/j.jtbi.2008.06.020.
Ogden NH, Maarouf A, Barker IK, Bigras-Poulin M, Lindsay LR, Morshed MG, O’Callaghan CJ, Ramay F, Waltner-Toews D, Charron DF: Climate change and the potential for range expansion of the Lyme disease vector Ixodes scapularis in Canada. Int J Parasitol. 2006, 36 (1): 63-70. 10.1016/j.ijpara.2005.08.016.
Ogden NH, Bigras-Poulin M, O’Callaghan CJ, Barker IK, Lindsay LR, Maarouf A, Smoyer-Tomic KE, Waltner-Toews D, Charron D: A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. Int J Parasitol. 2005, 35 (4): 375-389. 10.1016/j.ijpara.2004.12.013.
Estrada-Peña A, Sánchez N, Estrada-Sánchez A: An assessment of the distribution and spread of the tick Hyalomma marginatum in the Western palearctic under different climate scenarios. Vector Borne Zoonotic Dis. 2012, 12 (9): 758-768. 10.1089/vbz.2011.0771.
Gray JS, Dautel H, Estrada-Pena A, Kahl O, Lindgren E: Effects of climate change on ticks and tick-borne diseases in Europe. Interdiscip Perspect Infect Dis. 2009, 2009: 593232-
Estrada-Pena A, Horak IG, Petney T: Climate changes and suitability for the ticks Amblyomma hebraeum and Amblyomma variegatum (Ixodidae) in Zimbabwe (1974-1999). Vet Parasitol. 2008, 151 (2-4): 256-267. 10.1016/j.vetpar.2007.11.014.
Margos G, Piesman J, Lane RS, Ogden NH, Sing A, Straubinger RK, Fingerle V: Borrelia kurtenbachii sp. nov.: a widely distributed member of the Borrelia burgdorferi sensu lato species complex in North America. Int J Syst Evol Microbiol. 2014, 64 (Pt1): 128-130.
Garza M, Feria ATP, Casillas EA, Sanchez-Cordero V, Rivaldi C-L, Sarkar S: Projected future distributions of vectors of Trypanosoma cruzi in North America under climate change scenarios. PLoS Negl Trop Dis. 2014, In Press
Sarkar S, Strutz SE, Frank DM, Rivaldi CL, Sissel B, Sanchez-Cordero V: Chagas disease risk in Texas. PLoS Negl Trop Dis. 2010, 4 (10): doi:10.1371/journal.pntd.0000836
Carcavallo RU, Curto De Casas SI: Some health impacts of global warming in South America: vector-borne diseases. J Epidemiol. 1996, 6 (4): S153-S157.
IPCC: Summary for Policymakers. Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Edited by: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL. 2007, Cambridge, United Kingdom and New York, NY, USA,http://www.ipcc.ch/publications_and_data/ar4/wg1/en/contents.html,
Asin S, Catala S: Development of Trypanosoma cruzi in Triatoma infestans: influence of temperature and blood consumption. J Parasitol. 1995, 81 (1): 1-7.
Wu X, Duvvuri VR, Lou Y, Ogden NH, Pelcat Y, Wu J: Developing a temperature-driven map of the basic reproductive number of the emerging tick vector of Lyme disease Ixodes scapularis in Canada. J Theor Biol. 2013, 319: 50-61.
Ogden NH, Radojevic M, Wu X, Duvvuri VR, Leighton PA, Wu J: Estimated effects of projected climate change on the basic reproductive number of the Lyme disease vector. Environ Health Perspect. 2014, (in Press) doi:10.1289/ehp.1307799
Hamer SA, Tsao JI, Walker ED, Hickling GJ: Invasion of the lyme disease vector Ixodes scapularis: implications for Borrelia burgdorferi endemicity. EcoHealth. 2010, 7 (1): 47-63. 10.1007/s10393-010-0287-0.
Patz JA, Graczyk TK, Geller N, Vittor AY: Effects of environmental change on emerging parasitic diseases. Int J Parasitol. 2000, 30 (12-13): 1395-1405. 10.1016/S0020-7519(00)00141-7.
Mantovani E, Costa IP, Gauditano G, Bonoldi VL, Higuchi ML, Yoshinari NH: Description of Lyme disease-like syndrome in Brazil: is it a new tick borne disease or Lyme disease variation?. Braz J Med Biol Res. 2007, 40 (4): 443-456. 10.1590/S0100-879X2006005000082.
Aviña A, CTiwari C, Williamson P, Oppong J, Atkinson S: Geospatial Analysis of Environmental Health. Edited by: Maantay JA, McLalterty S. 2011, Springer, 357-371. A Spatially Explicit Environmental Health Surveillance Framework for tick-borne diseases, Chapter 18, Volume 4.
Hamer SA, Hickling GJ, Sidge JL, Rosen ME, Walker ED, Tsao JI: Diverse Borrelia burgdorferi strains in a bird-tick cryptic cycle. Appl Environ Microb. 2011, 77 (6): 1999-2007. 10.1128/AEM.02479-10.
Rawlings JA: Lyme disease in Texas. Zentralblatt fur Bakteriologie, Mikrobiologie, und Hygiene Series A, Medical microbiology, infectious diseases, virology, parasitology. 1987, 263 (3): 483-487. 10.1016/S0176-6724(87)80115-3.
Rawlings JA, Fournier PV, Teltow GJ: Isolation of Borrelia spirochetes from patients in Texas. J Clin Microbiol. 1987, 25 (7): 1148-1150.
Rawlings JA, Teltow GJ: Prevalence of Borrelia (Spirochaetaceae) spirochetes in Texas ticks. J Med Entomol. 1994, 31 (2): 297-301.
Teltow GJ, Fournier PV, Rawlings JA: Isolation of Borrelia burgdorferi from arthropods collected in Texas. Am J Trop Med Hyg. 1991, 44 (5): 469-474.
Diuk-Wasser MA, Hoen AG, Cislo P, Brinkerhoff R, Hamer SA, Rowland M, Cortinas R, Vourc’h G, Melton F, Hickling GJ, Tsao JI, Bunikis J, Barbour AG, Kitron U, Piesman J, Fish D: Human risk of infection with Borrelia burgdorferi, the Lyme disease agent, in eastern United States. Am J Trop Med Hyg. 2012, 86 (2): 320-327. 10.4269/ajtmh.2012.11-0395.
Pepin KM, Eisen RJ, Mead PS, Piesman J, Fish D, Hoen AG, Barbour AG, Hamer S, Diuk-Wasser MA: Geographic variation in the relationship between human Lyme disease incidence and density of infected host-seeking Ixodes scapularis nymphs in the Eastern United States. Am J Trop Med Hyg. 2012, 86 (6): 1062-1071. 10.4269/ajtmh.2012.11-0630.
Gordillo-Pérez G, Torres J, Solórzano-Santos F, Garduño-Bautista V, Tapia-Conyer R, Muñóz O: Estudio seroepidemiológico de Borreliosis de Lyme en la Ciudad de México y el Noreste de la República Mexicana. Salud Pública Mexicana. 2003, 45: 351-355.
Gordillo-Pérez G, Vargas-Sandoval M, Sosa-Gutierrez C, Minero-Gonzalez E, Schoeder-Lima E, Parra-Montiel G, Guzmán-Cornejo C, Sánchez-Cordero V, Torres J: Prevalencia de Infección de Borrelia burgdorferi y Ehrlichia spp en garrapatas y roedores provenientes de tres parques nacionales del Centro de la República Mexicana. Acarologia Latinoamerica. 2012, 1: 291-295.
Martinez A, Salinas A, Martinez F, Cantu A, Miller DK: Serosurvey for selected disease agents in white-tailed deer from Mexico. J Wildl Dis. 1999, 35 (4): 799-803. 10.7589/0090-3558-35.4.799.
Illoldi-Rangel P, Rivaldi CL, Sissel B, Trout Fryxell R, Gordillo-Perez G, Rodriguez-Moreno A, Williamson P, Montiel-Parra G, Sanchez-Cordero V, Sarkar S: Species distribution models and ecological suitability analysis for potential tick vectors of Lyme disease in Mexico. J Trop Med. 2012, 2012: 959101-
Gordillo-Pérez MG, Solórzano-Santos F: Enfermedad de Lyme: experiencia en niños mexicanos. Boletin Médico del Hospital Infantil de México. 2010, 67: 164-176.
Bosch J, Munoz MJ, Martinez M, de la Torre A, Estrada-Pena A: Vector-borne pathogen spread through ticks on migratory birds: a probabilistic spatial risk model for South-Western europe. Transbound Emerg Dis. 2013, 60 (5): 403-415. 10.1111/j.1865-1682.2012.01358.x.
Cartín-Rojas A: International Trade from Economic and Policy Perspective. Edited by: Bobem V. 2012, InTech, 143-166. doi:10.5772/48151, Transboundary Animal Diseases and International Trade, Chapter 7.
Piesman J, Eisen L: Prevention of tick-borne diseases. Annu Rev Entomol. 2008, 53: 323-343. 10.1146/annurev.ento.53.103106.093429.
Porretta D, Mastrantonio V, Amendolia S, Gaiarsa S, Epis S, Genchi C, Bandi C, Otranto D, Urbanelli S: Effects of global changes on the climatic niche of the tick Ixodes ricinus inferred by species distribution modelling. Parasit Vectors. 2013, 6: 271-10.1186/1756-3305-6-271.
Rubel F, Kottek M: Observed and projected climate shifts 1901-2100 depicted by world maps of the Köppen-Geiger climate classification. Meteorol Z. 2010, 19 (2): 135-141. 10.1127/0941-2948/2010/0430.
KOTTEK M, GRIESER J, BECK C, RUDOLF B, RUBEL F: World Map of the Köppen-Geiger climate classification updated. Meteorol Z. 2006, 15 (3): 259-263. 10.1127/0941-2948/2006/0130.
Correll DS, Johnston MC: Manual of the vascular plants of Texas. 1970, Renner: Texas Research Foundation
Archer S, Scifres CJ, Bassham CR, Maggio R: Autogenic succession in a subtropical savanna: conversion of grassland to thornwood land. Ecol Monogr. 1988, 58 (2): 111-127. 10.2307/1942463.
de Leon AA P, Teel PD, Auclair AN, Messenger MT, Guerrero FD, Schuster G, Miller RJ: Integrated strategy for sustainable cattle fever tick eradication in USA is required to mitigate the impact of global change. Front Physiol. 2012, 3: 195-
CONABIO: Actividades económicas en México por municipio, 2010’, escala: 1:250000. Comisión Nacional para el Conocimiento y Uso de la Biodiversidad. Edited by: Datos estadísticos del 2010 dINdEyGI. 2012, México D.F
Billings AN, Rawlings JA, Walker DH: Tick-borne diseases in Texas: a 10-year retrospective examination of cases. Texas Med. 1998, 94 (12): 66-76.
Keirans JE, Clifford CM: The genus Ixodes in the United States: a scanning electron microscope study and key to the adults. J Med Entomol Suppl. 1978, 2: 1-149.
Keirans JE, Durden LA: Illustrated key to nymphs of the tick genus Amblyomma (Acari: Ixodidae) found in the United States. J Med Entomol. 1998, 35 (4): 489-495.
Keirans JE, Litwak TR: Pictorial key to the adults of hard ticks, family Ixodidae (Ixodida: Ixodoidea), east of the Mississippi River. J Med Entomol. 1989, 26 (5): 435-448.
United States. Animal and Plant Health Inspection Service. Veterinary Services, Strickland RK, Diamant G: Ticks of veterinary importance. 1976, Washington: Animal and Plant Health Inspection Service, U.S. Dept. of Agriculture: for sale by the Supt. of Docs., U.S. Govt. Print. Off
Vargas M, Gordillo-Pérez G, Solórzano S, Rivera A, Polaco O, Muñoz O, Torres J: Evidences of Borrelia burgdorferi in ticks of the Northeast of Mexico. Entomologia Mexicana. 2007, 6: 830-835.
Almazan C, Castro-Arellano I, Camacho-Puga E: Black-Legged ticks (Ixodes scapularis) on the jaguar (Panthera onca). Southwest Nat. 2013, 58 (1): 126-127. 10.1894/0038-4909-58.1.126.
Dennis DT, Nekomoto TS, Victor JC, Paul WS, Piesman J: Reported distribution of Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae) in the United States. J Med Entomol. 1998, 35 (5): 629-638.
Diuk-Wasser MA, Brown HE, Andreadis TG, Fish D: Modeling the spatial distribution of mosquito vectors for West Nile virus in Connecticut. USA Vector Borne Zoonotic Dis. 2006, 6 (3): 283-295. 10.1089/vbz.2006.6.283.
Diuk-Wasser MA, Vourc’h G, Cislo P, Hoen AG, Melton F, Hamer SA, Rowland M, Cortinas R, Hickling GJ, Tsao JI, Barbour AG, Kilton U, Piesman J, Fish D: Field and climate-based model for predicting the density of host-seeking nymphal Ixodes scapularis, an important vector of tick-borne disease agents in the eastern United States. Global Ecol Biogeogr. 2010, 19 (4): 504-514.
Phillips SJ: Transferability, sample selection bias and background data in presence-only modelling: a response to Peterson et al. (2007). Ecography. 2008, 31 (2): 272-278. 10.1111/j.0906-7590.2008.5378.x.
Phillips SJ, Dudik M: Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008, 31 (2): 161-175. 10.1111/j.0906-7590.2008.5203.x.
Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ: A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011, 17 (1): 43-57. 10.1111/j.1472-4642.2010.00725.x.
Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190 (3-4): 231-259. 10.1016/j.ecolmodel.2005.03.026.
Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Powell GVN, Underwood EC, D’Amico JA, Itoua I, Strand HE, Morrison JC, Loucks CJ, Allnutt TF, Ricketts TH, Kura Y, Lamoreux JF, Wettengel WW, Hedao P, Kassem KR: Terrestrial ecoregions of the worlds: a new map of life on Earth. Bioscience. 2001, 51 (11): 933-938. 10.1641/0006-3568(2001)051[0933:TEOTWA]2.0.CO;2.
Fielding HA, Bell JF: A review of methods for the assessment of predictions errors in conservation presence/absence models. Environ Conserv. 1998, 24: 38-49.
Brownstein JS, Holford TR, Fish D: A climate-based model predicts the spatial distribution of the Lyme disease vector Ixodes scapularis in the United States. Environ Health Perspect. 2003, 111 (9): 1152-1157. 10.1289/ehp.6052.
Dantas-Torres F, Chomel BB, Otranto D: Ticks and tick-borne diseases: a one health perspective. Trends Parasitol. 2012, 28 (10): 437-446. 10.1016/j.pt.2012.07.003.
Hotez PJ, Gurwith M: Europe’s neglected infections of poverty. Int J Infect Dis. 2011, 15 (9): e611-e619. 10.1016/j.ijid.2011.05.006.
Hotez PJ: Fighting neglected tropical diseases in the southern United States. BMJ. 2012, 345: e6112-10.1136/bmj.e6112.
Andrus J, Bottazzi ME, Chow J, Goraleski KA, Fisher-Hoch SP, Lambuth JK, Lee BY, Margolis HS, McCormick JB, Melby P, Murray KO, Rico-Hesser R, Valenzuela JG, Hotez PJ: Ears of the Armadillo: global health research and neglected diseases in Texas. PLoS Negl Trop Dis. 2013, 7 (6): e2021-10.1371/journal.pntd.0002021.
Escalante T, Rodríguez-Tapia G, Szumik C, Morrone JJ, Rivas M: Delimitation of the Nearctic region according to mammalian distributional patterns. J Mammal. 2010, 91 (6): 1381-1388. 10.1644/10-MAMM-A-136.1.
Killilea ME, Swei A, Lane RS, Briggs CJ, Ostfeld RS: Spatial dynamics of lyme disease: a review. EcoHealth. 2008, 5 (2): 167-195. 10.1007/s10393-008-0171-3.
Gatewood AG, Liebman KA, Vourc’h G, Bunikis J, Hamer SA, Cortinas R, Melton F, Cislo P, Kitron U, Tsao J, Barbour AG, Fish D, Diuk-Wasser MA: Climate and tick seasonality are predictors of Borrelia burgdorferi genotype distribution. Appl Environ Microbiol. 2009, 75 (8): 2476-2483. 10.1128/AEM.02633-08.
Rydzewski J, Mateus-Pinilla N, Warner RE, Hamer S, Weng HY: Ixodes scapularis and Borrelia burgdorferi among diverse habitats within a natural area in east-central Illinois. Vector Borne Zoonotic Dis. 2011, 11 (10): 1351-1358. 10.1089/vbz.2010.0160.
Hamer SA, Hickling GJ, Sidge JL, Walker ED, Tsao JI: Synchronous phenology of juvenile Ixodes scapularis, vertebrate host relationships, and associated patterns of Borrelia burgdorferi ribotypes in the midwestern United States. Ticks Tick Borne Dis. 2012, 3 (2): 65-74. 10.1016/j.ttbdis.2011.11.004.
Reed KD, Meece JK, Henkel JS, Shukla SK: Birds, migration and emerging zoonoses: west nile virus, lyme disease, influenza A and enteropathogens. Clinl Med Res. 2003, 1 (1): 5-12. 10.3121/cmr.1.1.5.
Clark KL, Leydet B, Hartman S: Lyme Borreliosis in human patients in Florida and Georgia, USA. Int J Med Sci. 2013, 10 (7): 915-931. 10.7150/ijms.6273.
Brownstein JS, Holford TR, Fish D: Effect of climate change on Lyme disease risk in North America. EcoHealth. 2005, 2 (1): 38-46. 10.1007/s10393-004-0139-x.
Pulliam H: On the relationship between niche and distribution. Ecol Lett. 2000, 3: 349-361. 10.1046/j.1461-0248.2000.00143.x.
Soberón J, Peterson AT: Interpretation of models of fundamental ecological niches and species’ distributional areas. Biodivers Informatics. 2005, 2: 1-10.
Brown JH, Stevens GC, Kaufman DM: The geographic range: size, shape, boundaries, and internal structure. Annu Rev Ecol Syst. 1996, 27: 597-623. 10.1146/annurev.ecolsys.27.1.597.
Austin MP: Spatial prediction of species distribution: an interface between ecological theory and statistical modelling. Ecol Model. 2002, 157: 101-118. 10.1016/S0304-3800(02)00205-3.
Diuk-Wasser MA, Gatewood AG, Cortinas MR, Yaremych-Hamer S, Tsao J, Kitron U, Hickling G, Brownstein JS, Walker E, Piesman J, Fish D: Spatiotemporal patterns of host-seeking Ixodes scapularis nymphs (Acari: Ixodidae) in the United States. J Med Entomol. 2006, 43 (2): 166-176. 10.1603/0022-2585(2006)043[0166:SPOHIS]2.0.CO;2.
Salinas-Melendez JA, Avalos-Ramirez R, Riojas-Valdez VM, Martinez-Munoz A: Serological survey of canine borreliosis. Rev Latinoam Microbiol. 1999, 41 (1): 1-3.
Bowman D, Little SE, Lorentzen L, Shields J, Sullivan MP, Carlin EP: Prevalence and geographic distribution of Dirofilaria immitis, Borrelia burgdorferi, Ehrlichia canis, and Anaplasma phagocytophilum in dogs in the United States: results of a national clinic-based serologic survey. Vet Parasitol. 2009, 160 (1–2): 138-148.
Salinas-Melendez JA, Zarate-Ramos JJ, Avalos-Ramirez R, Hernandez-Escareno JJ, Guzman-Acosta G, Riojas-Valdes VM, Segura-Correa JC: Prevalence of Antibodies Against Borrelia burgdorferi in Dogs from Monterrey, Mexico. J Anim Vet Adv. 2011, 10 (20): 2720-2723.
Healy JA, Cross TF, Healy A: The alpha-Gpdh polymorphism in the tick Ixodes ricinus: similar diurnal trends in genotypic composition in Irish and Swedish population samples. Exp Appl Acarol. 2004, 32 (1-2): 111-118.
Radulovic Z, Milutinovic M, Tomanovic S, Mihaljica D, Cakic S, Stamenkovic-Radak M, Andelkovic M: Seasonal and spatial occurrence of glycerol-3-phosphate dehydrogenase variability in Ixodes ricinus (Acari: Ixodidae) populations. J Med Entomol. 2012, 49 (3): 497-503. 10.1603/ME11101.
Radulovic Z, Milutinovic M, Andelkovic M, Vujcic Z, Tomanovic S, Bozic N, Marinkovic D: Allozyme polymorphism of Mdh and alpha-Gpdh in Ixodes ricinus populations: comparison of borreliae-infected and uninfected ticks. Exp Appl Acarol. 2006, 40 (2): 113-121. 10.1007/s10493-006-9028-z.
Schulze TL, Jordan RA, Hung RW, Schulze CJ: Effectiveness of the 4-Poster passive topical treatment device in the control of Ixodes scapularis and Amblyomma americanum (Acari: Ixodidae) in New Jersey. Vector Borne Zoonotic Dis. 2009, 9 (4): 389-400. 10.1089/vbz.2008.0160.
Wilson DE, Ruff S: The Smithsonian Book of North American mammals. 1999, Washington: Smithsonian Institution Press
Anderson JF, Magnarelli LA: Biology of ticks. Infect Dis Clin North Am. 2008, 22 (2): 195-215. 10.1016/j.idc.2007.12.006. v
Kempf F, Boulinier T, De Meeus T, Arnathau C, McCoy KD: Recent evolution of host-associated divergence in the seabird tick Ixodes uriae. Mol Ecol. 2009, 18 (21): 4450-4462. 10.1111/j.1365-294X.2009.04356.x.
Kempf F, De Meeus T, Vaumourin E, Noel V, Taragel’ova V, Plantard O, Heylen DJ, Eraud C, Chevillon C, McCoy KD: Host races in Ixodes ricinus, the European vector of Lyme borreliosis. Infect Gen Evol. 2011, 11 (8): 2043-2048. 10.1016/j.meegid.2011.09.016.
Dietrich M, Beati L, Elguero E, Boulinier T, McCoy KD: Body size and shape evolution in host races of the tick Ixodes uriae. Biol J Linn Soc. 2013, 108: 323-334. 10.1111/j.1095-8312.2012.02021.x.
Authors want to thank the Department of Texas Parks and Wildlife for their help providing tick samples from a number of wildlife management areas. Also authors would like to thank Veterinary clinics across the state of Texas, and the Brazos Valley animal shelter (currently Aggieland Humane Society) and Drs. Thomas Craig and Karen Snowden for their help collecting ticks for this study. USDA is an equal opportunity provider and employer. Support for this study came from the Department of Veterinary Pathobiology, Texas A&M University, the AgriLife-TVMDL seed grant to the project entitled “Improving diagnostic methods for Lyme disease, and epidemiology of human and animal infections in TX” and the bi-national cooperative grant TAMU-CONACYT-052 ‘Typing virulent isolates of Lyme disease agents in central Mexico’ to MDEG. CONACYT–SALUD 1999-30694-M and 2008-87686. ICA received a Research Enhancement Program Grant from Texas State University. JT is a recipient of an exclusivity scholarship from Fundacion IMSS, Mexico.
The authors declare that they have no competing interests.
TPF participated in the design of the study, carried out the organization of data bases, developed the distribution models presented in the paper and participated in the organization and drafting of this manuscript. ICA and APL participated in the study design, data base analysis, and in the organization and drafting of this manuscript. GGP participated in the design of the study, coordinated the sample collection in Mexico, provided the Mexican database, and participated in the organization and drafting of this manuscript. ALC carried out the data analysis and participated in the generation of distribution models presented. MVS participated in the design of the fieldwork carried out in Mexico, identified tick species and provided the distribution data. AG participated in the design of the field work carried out in Texas and collected samples, analyzed DNA by PCR, analyzed sequencing sequences and provided data for Figure 1 of this paper. JT, coordinated the analysis of Mexican samples and participated in the organization and drafting of this manuscript, RFM participated in the organization and drafting of this manuscript. MDEG participated in the design of the study, coordinated and participated in the sample collection and analysis of the Texas samples, provided the Texan database, coordinated the organization and drafting of this manuscript as well as the acquisition of funds to carry out the study. All authors read and approved the final manuscript.
Teresa P Feria-Arroyo, Ivan Castro-Arellano, Guadalupe Gordillo-Perez contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Feria-Arroyo, T.P., Castro-Arellano, I., Gordillo-Perez, G. et al. Implications of climate change on the distribution of the tick vector Ixodes scapularis and risk for Lyme disease in the Texas-Mexico transboundary region. Parasites Vectors 7, 199 (2014). https://doi.org/10.1186/1756-3305-7-199