Hostname: page-component-7b9c58cd5d-v2ckm Total loading time: 0 Render date: 2025-03-15T22:30:14.857Z Has data issue: false hasContentIssue false

Higher sika deer density is associated with higher local abundance of Haemaphysalis longicornis nymphs and adults but not larvae in central Japan

Published online by Cambridge University Press:  23 May 2013

Hideharu Tsukada*
Affiliation:
NARO Institute of Livestock and Grassland Science, Shiono, Miyota, Nagano 389-0201, Japan
Yoshio Nakamura
Affiliation:
NARO Institute of Animal Health, Kannondai, Tsukuba, Ibaraki 305-0856, Japan
Tsugihiko Kamio
Affiliation:
Kyushu Research Station, NARO Institute of Animal Health, Chuzan, Kagoshima 891-0105, Japan
Hisashi Inokuma
Affiliation:
Department of Clinical Veterinary Medicine, Obihiro University of Agriculture and Veterinary Medicine, Inada-cho, Obihiro, Hokkaido 080-8555, Japan
Yasuko Hanafusa
Affiliation:
NARO Institute of Animal Health, Kannondai, Tsukuba, Ibaraki 305-0856, Japan
Naoko Matsuda
Affiliation:
Tochigi Prefectural Citizen’s Park Management Office, Nagai, Yaita, Tochigi 329-2514, Japan
Tetsuya Maruyama
Affiliation:
Department of Environment and Forestry, Tochigi Prefectural Office, Tochigi 320-8501, Japan
Takahiro Ohba
Affiliation:
Shizuoka Prefecture Forestry and Forest Products Research Institute, Hamamatsu 434-0016, Japan
Koji Nagata
Affiliation:
Kanagawa Prefecture Natural Environment Conservation Center, Nanasawa, Atsugi-shi, Kanagawa 243-0121, Japan
*
*Author for correspondence Phone: +81 (0)267 32 0757 Fax: +81 (0)267 32 2318 E-mail: htsuka@affrc.go.jp
Rights & Permissions [Opens in a new window]

Abstract

Haemaphysalis longicornis (Acari: Ixodidae) is one of the most common and important arthropod disease vectors in Japan, carrying Japanese spotted fever and bovine theileriosis. The recent expansion of sika deer (Cervus nippon, Artiodactyla: Cervidae) populations, the most common wild host of H. longicornis, has also caused concern about increasing the risk of vector-borne diseases in Japan. We used generalized linear mixed model analysis to determine the relative contribution of deer density and other biological and abiotic factors on the abundance of H. longicornis ticks questing at each developmental stage. A total of 6223 H. longicornis adults, nymphs, and larvae were collected from 70 sites in three regions of central Japan. The abundance of questing adult and nymphal ticks was associated with deer density and other biotic and abiotic factors. However, the abundance of questing larvae showed no association with deer density but did show an association with other biotic and abiotic factors. These findings show that a high density of deer along with other biotic and abiotic factors is associated with increased risk of vector-borne diseases through amplified local abundance of questing nymphal and adult H. longicornis. Further, questing larvae abundance is likely regulated by environmental conditions and is likely correlated with survival potential or the distribution of other host species.

Type
Research Paper
Copyright
Copyright © Cambridge University Press 2013 

Introduction

Haemaphysalis longicornis (Acari: Ixodidae) is one of the most common and important arthropod disease vectors for human and animal health in Japan. This species is known to be a carrier of the causative agents of Japanese spotted fever (Mahara, Reference Mahara1997) and bovine theileriosis (Ishihara, Reference Ishihara1968; Kamio et al., Reference Kamio, Fujisaki and Minami1989). This species has a wide distribution in Far East Asia and has also been introduced into the Australian and Oceania regions (Hoogstraal et al., Reference Hoogstraal, Roberts, Kohls and Tipton1968). Its distribution extends over all of Japan (Yamaguchi et al., Reference Yamaguchi, Tipton, Keegan and Toshioka1971), and it has a variety of animal hosts, ranging from small to large mammals and birds (Yamaguchi et al., Reference Yamaguchi, Tipton, Keegan and Toshioka1971).

Grazing cattle are the major hosts of H. longicornis in grasslands, including on pastures in Japan (Kitaoka et al., Reference Kitaoka, Morii and Fujisaki1975; Yamane et al., Reference Yamane, Nishiguchi, Kobayashi and Zeniya2006), and sika deer (Cervus nippon, Artiodactyla: Cervidae) in forested areas is considered to be one of the most important wild host species for this tick species at all developmental stages (Mori et al., Reference Mori, Tsunoda and Fujimagari1995; Inokuma et al., Reference Inokuma, Fujimoto, Hosoi, Tanaka, Fujisaki, Okuda and Onishi2002; Yamauchi et al., Reference Yamauchi, Tabara, Kanamori, Kawabata, Arai, Katayuama, Fujita, Yano, Takada and Itagaki2009). Sika deer is the only native species of the family Cervidae in Japan (Ohdachi et al., Reference Ohdachi, Ishibashi, Iwasa and Saitoh2009) and has recently expanded in range and numbers (Takatsuki, Reference Takatsuki2009), causing serious problems such as agricultural damage (Yamane et al., Reference Yamane, Nakamura and Tsukada2009; Tsukada, Reference Tsukada2011), destruction of natural vegetation (Takatsuki, Reference Takatsuki2009), and deterioration of ecosystem function (Takatsuki & Maeyama, Reference Takatsuki and Maeyama2007). Likewise, the expansion of areas with tick-borne endemic diseases and an increase in disease occurrences related to the drastic changes in distribution and density of the host sika deer are of concern from the perspective of animal health (Inokuma et al., Reference Inokuma, Fujimoto, Hosoi, Tanaka, Fujisaki, Okuda and Onishi2002, Reference Inokuma, Seino, Suzuki, Kaji, Takahashi, Igota and Inoue2008). Although sika deer has not been confirmed to be a reservoir of the causative agent of bovine theileriosis (Inokuma et al., Reference Inokuma, Tsuji, Kim, Fujimoto, Nagata, Hosoi, Arai, Ishihara and Okuda2004) and Japanese spotted fever (Inokuma et al., Reference Inokuma, Seino, Suzuki, Kaji, Takahashi, Igota and Inoue2008), sika deer can be an amplifier host of the tick carrying these pathogens (Yamauchi et al., Reference Yamauchi, Tabara, Kanamori, Kawabata, Arai, Katayuama, Fujita, Yano, Takada and Itagaki2009; Inokuma & Yokoyama, Reference Inokuma and Yokoyama2010).

While some studies on tick fauna have shown that there is a higher abundance of adult and nymphal H. longicornis in areas with higher densities of sika deer (Yamauchi et al., Reference Yamauchi, Tabara, Kanamori, Kawabata, Arai, Katayuama, Fujita, Yano, Takada and Itagaki2009), another has shown that the abundance of larval H. longicornis is also negatively affected by the decrease in plant coverage, plant height, and litter depth caused by the presence of sika deer in high density (Tsunoda, Reference Tsunoda2007b ). These previous studies suggest that sika deer density can have different effects on the abundances of questing H. longicornis in different stages. In addition, the abundance of H. longicornis is influenced by various abiotic factors, including temperature (Fujisaki et al., Reference Fujisaki, Kitaoka and Morii1975; Heath, Reference Heath1981; Sutherst & Bourne, Reference Sutherst and Bourne1991; Fujimoto, Reference Fujimoto1993), rainfall (Yamane et al., Reference Yamane, Nishiguchi, Kobayashi and Zeniya2006), humidity (Heath, Reference Heath1981; Fujimoto, Reference Fujimoto1988; Tsunoda, Reference Tsunoda2008), and photoperiod (Fujimoto, Reference Fujimoto2000, Reference Fujimoto2008).

In order to develop effective control measures against bovine theileriosis among grazing cattle on livestock farms in Japan, it is important to determine the relative contribution of sika deer density and other biological and abiotic factors, including climatic and geographical parameters, on H. longicornis individuals questing at each stage. Here we use generalized linear mixed model (GLMM) analysis to evaluate these effects at three sites in the central region of Honshu.

Materials and methods

Study area

The study was conducted in the three regions in central Japan (fig. 1): Tochigi Prefecture and Izu Peninsula in Shizuoka Prefecture in 2008, and in Kanagawa Prefecture in 2009. These regions were selected from areas in which population density surveys of sika deer were already being conducted by each respective prefectural government. Within each sampling area, 27 (mean size, 1.41 km2; range, 0.78–2.09 km2), 21 (size data are not available; line transects totaling 300 m were surveyed in each area) and 22 (mean size, 1.65 km2; range, 1.20–2.09 km2) sampling sites, respectively, were used for analysis (fig. 1). These sampling sites are all located in forested areas of various types (Supplementary data) and did not include any pasture land used for grazing cattle.

Fig. 1. Locations of tick sampling sites in three regions of central Japan: Tochigi Prefecture, Izu Peninsula, and Kanagawa Prefecture. Sampling sites in each region are represented by the same symbol, and thin lines indicate the prefecture borders.

Deer population density

We used data from prefectural surveys of sika deer in Tochigi Prefecture, Izu Peninsula, and Kanagawa Prefecture.

Tochigi Prefecture

In Tochigi Prefecture, deer population surveys had been conducted annually at 30 sites by the block count method (Maruyama & Furubayashi, Reference Maruyama and Furubayashi1983) from August to October since 2003. We selected 27 sites for our tick sampling survey, excluding three sites that were inaccessible for tick sampling. We used averaged density of sika deer at each sampling site over the period of 2003–2007 for GLMM analysis because we assumed that there was annual fluctuation in the population as well as sampling error.

Izu Peninsula

On Izu Peninsula in Shizuoka Prefecture, sampling surveys to estimate the density of sika deer had been performed annually at about 10–20 sites out of 72 sites since 2001 using the fecal pellet count method and analysis by the program FUNRYU (Iwamoto et al., Reference Iwamoto, Sakata, Nakazono, Kaoka, Ikeda, Nishishita, Tokida and Doi2000) along line transects (300 m site−1). We selected the 21 sites where deer fecal pellet count surveys had been performed twice from December 2006 to January 2007 and February to March 2007 for our tick sampling survey, excluding one sampling site that was difficult to access (fig. 1). We used sika deer densities averaged over two samplings for each site for later GLMM analysis, and these data were calculated using the program FUNRYU (Iwamoto et al., Reference Iwamoto, Sakata, Nakazono, Kaoka, Ikeda, Nishishita, Tokida and Doi2000).

Kanagawa Prefecture

In Kanagawa Prefecture, deer population surveys had been conducted annually at 11–24 sites out of 52 sites by the block count method (Maruyama & Furubayashi, Reference Maruyama and Furubayashi1983) during the winter (December–January) from 2003 to 2008 (Wildlife Management Office, 2009). We selected 22 out of 27 sites where the deer population survey had been conducted in 2008 for our tick sampling survey, excluding five sites that were inaccessible for tick sampling (fig. 1). We used averaged densities of sika deer at each sampling site in the period of 2004–2008 for our later analysis because we assumed that that there was annual fluctuation in the population as well as sampling error.

Ticks

The locations of the deer sampling sites in the three regions were entered into in a handheld GPS (GPSMAP 60CSx, Garmin, Ltd, Olathe, KS) in advance of tick sampling, and the exact location of each sampling survey was recorded. Questing ticks were collected by flagging with a white cotton flannel cloth (0.7×0.7 m) on the vegetation for 30 min at each tick sampling survey site in spring (May–June), summer (July–August), and autumn (September–October) in order to sample for seasonal activity and tick stages (Fujimoto, Reference Fujimoto2000, Reference Fujimoto2004; Tsunoda, Reference Tsunoda2007a ). Tick-sampling surveys were conducted daily (weather permitting) for about a week in each season, and sampling was conducted in the daytime but with variation in the timing among the sites due to logistical limitations. The ticks collected during each sampling survey were placed in vials with 70% ethanol for later examination in the laboratory. Ticks were examined for morphological characteristics, identified by stage, sex, and species according to references (Fujita & Takada, Reference Fujita and Takada2007; Yamaguchi et al., Reference Yamaguchi, Tipton, Keegan and Toshioka1971), and then enumerated.

Climatic and geographical data

A total of six climatic and geographical variables were analyzed as potential abiotic factors for the local abundance of ticks: (1) altitude, (2) mean temperature, (3) rainfall, (4) duration of sunshine without shading, (5) global solar radiation without shading, and (6) depth of snow cover. Data sets for these parameters were obtained from the Mesh Climatic Data 2000 (Japan Meteorological Agency, 2002). The data sets include various climatic variables in a 1-km mesh for the period 1971–2000 that were used as Japanese climate normals until May 2011 by the Japan Meteorological Agency.

Model development

We adopted a model selection approach to assess the factors influencing the abundance of ticks at each stage (Burnham & Anderson, Reference Burnham and Anderson2002; Johnson & Omland, Reference Johnson and Omland2004). We made several sets of working hypotheses, including that local H. longicornis tick abundance at each stage (response variables: abundance of larval, nymphal, and adult H. longicornis) was influenced by three biological factors: (1) density of sika deer, (2) abundance of H. longicornis in the previous stage, and (3) abundance of Haemaphysalis megaspinosa ticks of the same stage, and six climatic-geographical factors: (1) altitude, (2) mean temperature, (3) rainfall, (4) duration of sunshine without shading, (5) global solar radiation without shading, and (6) depth of snow cover. To investigate these working hypotheses, we applied a GLMM with a negative binomial error structure to the data, which were log transformed with the statistical software R (ver. 2.13.2, R Development Core Team, 2011) using the glmmADMB package (Skaug et al., Reference Skaug, Fournier, Nielsen, Magnusson and Bolker2012). In these models, the three regions (Tochigi Prefecture, Izu Peninsula, and Kanagawa Prefecture) were regarded as a random effect, along with the density of sika deer, which had random variable intercepts and slopes. Since the methods of estimating the density of sika deer differed among the regions, we assumed that the association between tick abundance and deer density also varied among the regions. As a preliminary analysis, we set the abundance of ticks in each stage as the response variable in a generalized linear model (GLM) and assumed three error distributions: normal, Poisson, and negative binominal. We confirmed that the GLM with a negative binominal error distribution showed the lowest Akaike information criterion (AIC) values in all stages.

We used AIC adjusted for small sample size (AICc; Burnham & Anderson, Reference Burnham and Anderson2002) to identify the subsets of parameters for the GLMM (table 1). AICc was calculated for all possible additive models starting with the full model, and models were ranked by AICc values. From the differences in AICc values (ΔAICc), AIC weights (ω i ) were calculated. For ΔAICc<2, which indicates approximately equal parsimony of models, all variables considered in the full model were ranked according to importance (predictor weights, ω+(j); Burnham & Anderson, Reference Burnham and Anderson2002). Unless a single model had a ω i >0.9, all other models were considered when interpreting the data (Burnham & Anderson, Reference Burnham and Anderson2002). To estimate the relative importance of predictor variables, x j (relative variable importance, RVI), the sum of the Akaike weights across the best models (ΔAICc<2) in the sets where variable j occurs was calculated. Using the RVI, the importance of all the variables was ranked (Burnham & Anderson, Reference Burnham and Anderson2002).

Table 1. Species and number of ixodid ticks collected by the flagging method in three regions of central Japan.

Results

Tick species and stage distribution

In total, 6144, 9980, and 12,533 ticks were collected in Tochigi Prefecture, Izu Peninsula, and Kanagawa Prefecture, respectively. The ticks (all Acari: Ixodidae) were morphologically classified into 4 genera and 12 species (table 1). H. megaspinosa (63.8%) was the most abundant species, followed by H. longicornis (21.9%), Haemaphysalis flava (11.3%), Haemaphysalis kitaokai (1.4%), Ixodes ovatus (0.6%), and Ixodes persulcatus (0.4%). The three major species comprised 96.9% of all ticks. Other species (all Acari: Ixodidae), such as Amblyomma testudinarium, Haemaphysalis japonica, Ixodes nipponensis, Ixodes monospinosus, Dermacentor taiwanensis, and Ixodes turdus, were collected in fewer numbers. The three major species were found in all three regions, but some species, such as I. persulcatus, I. monospinosus, and I. turdus, were site-specific.

The number of H. longicornis ticks in each stage was plotted against the deer density (fig. 2). This simple nonparametric correlation analysis showed that the abundance of nymphs and adults was positively correlated with deer density (nymph, Kendall's τ = 0.28, P < 0.001; adult, Kendall's τ = 0.19, P < 0.05), but the abundance of larvae was not (Kendall's τ = 0.05, P>0.5). The seasonal abundance of H. longicornis varied with stage, with larvae, nymphs, and adults being collected mainly in autumn, spring, and summer, respectively (table 2).

Fig. 2. Relationship between H. longicornis tick abundance by stage (larval, nymphal, and adult) and sika deer density at each tick-sampling site.

Table 2. Total number of H. longicornis larvae, nymphs, and adults collected by season

GLMM analysis of factors influencing the abundance of questing ticks

A total of 9, 12, and 3 best-ranked models (that is, ΔAICc < 2) were identified for larval, nymphal, and adult stages, respectively (table 3). No single model emerged as the top-ranking model (that is, ω i  > 0.9) among all stages. Therefore, the differences among these best-ranked models are considered not to be substantial (Burnham & Anderson, Reference Burnham and Anderson2002). The ω i of the null model including only the intercept were almost all 0 for all tick stages, and the differences of AICc between the least AICc of the best model and the null model were 7.2, 13.6, and 19.2 for larval, nymphal, and adult stages, respectively. The differences of AICc among these stages indicate that the best-ranked models are better than the null model for all tick stages.

Table 3. AICc ranking of full models using biological and abiotic factors to estimate dependence of total number of H. longicornis larvae, nymphs, and adults in GLMMs. Null models including only the intercept are also shown for each tick stage.

K, Number of parameters; ΔAICc, difference in AICc between best and the actual model; ω i , Akaike's weight; Evidence ratios, ratio of Akaike's weight of the best and the actual model. Bold letters indicate models with ΔAICc<2. LarvHl: tick abundance of larval H. longicornis; NympHl: tick abundance of nymphal H. longicornis; AdulHl: tick abundance of adult H. longicornis; Deer Dens: density of sika deer; LarvHm: tick abundance of larval H. megaspinosa; NympHm: tick abundance of nymphal H. megaspinosa; Alt: altitude; Temp: mean temperature; Rain: rainfall; DayHour: duration of sunshine without shading; Solar: global solar radiation without shading; SnowDep: depth of snow cover.

The estimated coefficients for the best-ranked models for each tick stage and the RVI of each predictor variable are shown in table 4. Deer density was positively associated with nymphal and adult but not with larval stages. For larvae, snow depth, altitude, and the larval abundance of H. megaspinosa had negative associations and solar radiation had a positive association with tick abundance, all showing relatively high RVIs. Mean temperature and duration of sunshine without shading had negative associations, while rainfall and adult abundance had positive associations with larval abundance, all showing relatively low RVIs. For nymphs, larval abundance, duration of sunshine without shading, and rainfall were positively associated and snow depth was negatively associated with tick abundance, all showing high RVIs. Nymphal abundance of H. megaspinosa and altitude had negative associations while solar radiation without shading had a positive association with the nymphal abundance, all showing low RVIs. In adults, nymphal abundance and altitude were positively associated and snow depth was negatively associated with tick abundance, all showing high RVIs.

Table 4. A parameter estimate of each factor relative to the best model (table 3) in the case of equally parsimonious models. Estimates of the relative importance of predictor variables (RVI) were calculated by summing the Akaike weights across the set of all models in which the variable occurs. Model structure IDs correspond to those in table 3.

See table 3 for model structure's numbers; RVI: relative variable importance; Int: Intercept; LarvHl: tick abundance of larval H. longicornis; NympHl: tick abundance of nymphal H. longicornis; AdulHl: tick abundance of adult H. longicornis; Deer Dens: density of sika deer; LarvHm: tick abundance of larval H. megaspinosa; NympHm: tick abundance of nymphal H. megaspinosa; Alt: altitude; Temp: mean temperature; Rain: rainfall; DayHour: duration of sunshine without shading; Solar: global solar radiation without shading; SnowDep: depth of snow cover

Discussion

The density of sika deer was confirmed to be positively associated with the local abundance of questing nymphal and adult H. longicornis through statistical analysis using model comparison with GLMM. These results are well anticipated. Since the range of movement of H. longicornis is limited to only several meters (Mori & Tsunoda, Reference Mori and Tsunoda1996), its distribution depends almost entirely on host movement. Ticks are virtually sessile organisms that depend greatly on host movement to act as a kind of vehicle (Ruiz-Fons & Gilbert, Reference Ruiz-Fons and Gilbert2010). The local abundance of tick species, including various species of ticks that target deer, was previously found to be affected by host density (Rand et al., Reference Rand, Lubelczyk, Lavigne, Elias, Holman, Lacombe and Smith2003; Stafford et al., Reference Stafford, Denicola and Kilpatrick2003; Elias et al., Reference Elias, Lubelczyk, Rand, Lacombe, Holman and Smith2006; Tagliapietra et al., Reference Tagliapietra, Rosà, Arnoldi, Cagnacci, Capelli, Montarsi, Hauffe and Rizzoli2011). However, questing larval abundance was not associated with deer density and was also not strongly associated with questing adult tick abundance. This result is unexpected at first glance; however, it indicates that questing larval abundance is much more strongly regulated by local environmental conditions, which are likely correlated with survival potential rather than by the major supplier factors (i.e., host and adult abundance). As the host range of H. longicornis is known to be relatively large (Hoogstraal et al., Reference Hoogstraal, Roberts, Kohls and Tipton1968), other host species (i.e., other than sika deer) may serve as suppliers of adult ticks that in turn lay eggs. The small sample size of adult ticks collected in this study might also affect the weak association with adult tick abundance. Larval abundance was shown to be negatively associated with the larval abundance of H. megaspinosa, which is another major parasite of wild sika deer (Inokuma et al., Reference Inokuma, Fujimoto, Hosoi, Tanaka, Fujisaki, Okuda and Onishi2002; Tsunoda & Tatsuzawa, Reference Tsunoda and Tatsuzawa2004) and quests high on plants by climbing up on a stalk of grass or on a bush and waiting for a host to pass as does H. longicornis. Tsunoda (Reference Tsunoda2007c ) revealed that the height in the grass above which H. longicornis larvae and nymphs quest is significantly decreased by the presence of H. megaspinosa. This previous interspecific interaction found between these two species was corroborated by our current results.

Larval abundance was strongly associated with several abiotic factors, including snow depth, solar radiation, altitude, and temperature. The positive solar radiation and negative temperature associations were unique to the larval stage, as low temperature could reduce larval activity and egg survival (Fujimoto, Reference Fujimoto1993), and it could prolong the duration of egg development and hatching (Fujisaki et al., Reference Fujisaki, Kitaoka and Morii1975; Sutherst & Bourne, Reference Sutherst and Bourne1991). Solar radiation, however, had not been shown to affect larval H. longicornis abundance in previous studies. There are few reports of the positive effect of solar radiation on larval feeding and adult questing activity of American dog tick, Dermacentor variabilis (Acari: Ixodidae) (Atwood & Sonenshine, Reference Atwood and Sonenshine1967), and there is a report of a positive effect of solar radiation with low temperature but a negative effect with high temperature on Ixodes ricinus (Acari: Ixodidae)–man interactions in Demark (Jensen & Jespersen, Reference Jensen and Jespersen2005). Although we could not identify the exact reason for the positive association between solar radiation and questing larval tick abundance in our study, solar radiation may indirectly affect larval abundance by suppressing harmful fungal activity through the damaging effects of ultraviolet B (UVB) radiation (Inglis et al., Reference Inglis, Goettel, Butt, Strasser, Butt, Jackson and Magan2001). Some species of eumycete and deuteromycete fungi are known to attack and kill ixoid ticks (Benjamin et al., Reference Benjamin, Zhioua and Ostfeld2002; Ostfeld et al., Reference Ostfeld, Price, Hornbostel, Benjamin and Keesing2006). In fact, the entomopathogenic fungi Beauveria bassiana (Sordariomycetes: Cordycipitaceae) and Metarhizium anisopliae (Sordariomycetes: Clavicipitaceae) were found to cause 100% mortality in Rhipicephalus appendiculatus and Amblyomma variegatum (both Acari: Ixodidae) larvae in laboratory tests (Kaaya & Hassan, Reference Kaaya and Hassan2000). Strong UVB radiation often severely reduces the persistence of fungal conidia in the field (Fernandes et al., Reference Fernandes, Bittencourt and Roberts2012). For example, the entomopathogenic fungi B. bassiana can be sterilized under exposure to direct sunlight for less than 3 h (Yamashiro, Reference Yamashiro2000). Although approximately 90% of solar UVB radiation can be absorbed by ozone, water vapor, oxygen, and carbon dioxide through the atmosphere (World Health Organization, 2002), global solar radiation is highly correlated with UVB radiation on the ground in Japan (Motokuni et al., 1999; Shimizu et al., Reference Shimizu, Yoshizawa, Egashira, Hidaka and Yamamoto2006).

Except for the association with adult abundance and deer density, H. longicornis nymphal abundance was positively associated with the duration of sunshine and rainfall. Fujimoto (Reference Fujimoto1995, Reference Fujimoto2000, Reference Fujimoto2008) experimentally clarified that host attachment activity of unfed nymphs was significantly suppressed by exposure to the short-day photoperiod and suggested that nymphs show low activity in autumn around Saitama Prefecture, central Japan. Such behavioral diapause or quiescence in unfed nymphs was also related to the positive association with the duration of sunshine on the questing nymphal abundance at our study sites. H. longicornis larvae were greatly affected by relative humidity, as they cannot molt at relative humidity <74% (Fujimoto, Reference Fujimoto1988). Many other tick species are vulnerable to desiccation in off-host conditions (Needham & Teel, Reference Needham and Teel1991; Chilton & Bull, Reference Chilton and Bull1993; Sutherst & Bourne, Reference Sutherst and Bourne2006), and their abundance and activities were also affected by relative humidity (Oorebeek & Kleindorfer, Reference Oorebeek and Kleindorfer2008; Ruiz-Fons & Gilbert, Reference Ruiz-Fons and Gilbert2010; Tagliapietra et al., Reference Tagliapietra, Rosà, Arnoldi, Cagnacci, Capelli, Montarsi, Hauffe and Rizzoli2011). Such sensitivity to moisture explains the positive effect of rainfall on questing nymphal abundance observed here.

The abundance of questing H. longicornis adults was positively correlated with altitude, which is contrary to the negative effect observed for larvae. As temperature decreases with higher altitude, the negative effect of altitude on larvae could be attributed to the negative effect of temperature. In other tick species, the negative effect of altitude on abundance could be explained by delayed activation of already molted individuals (Jouda et al., Reference Jouda, Perret and Gern2004). Low temperature has also been reported to result in delayed development of eggs and larvae of H. longicornis (Fujisaki et al., Reference Fujisaki, Kitaoka and Morii1975; Sutherst & Bourne, Reference Sutherst and Bourne1991). In contrast, temperature did not show a clear effect on questing adult abundance. Another possible explanation is an indirect effect of host movement on questing adult distribution. The distribution of adult ticks partly depends on the movement of hosts infested with engorged nymphs. The unfed nymphs are usually found questing high on plants from spring to early summer (Tsunoda, Reference Tsunoda2007a ) and then infest a primary host, sika deer. In central Japan from early spring to early summer, some populations of sika deer move to areas that are higher than 2000 m above sea level (asl) to forage in alpine meadows and then move down to lower altitudes for overwintering (Izumiyama & Mochizuki, Reference Izumiyama and Mochizuki2008). Similar seasonal movements of sika deer have also been observed in the northern part of Japan (Takatsuki et al., Reference Takatsuki, Suzuki and Higashi2000). Such seasonal host movement may indirectly affect adult abundance through the positive effect of altitude. In a spatially stratified large-scale tick study, significant negative effects of altitude were observed on host-seeking Ixodes scapularis (Acari: Ixodidae) abundance while controlling for temperature-derived variables (Diuk-Wasser et al., Reference Diuk-Wasser, Vourc'h, Cislo, Hoen, Melton, Hamer, Rowland, Cortinas, Hickling, Tsao, Barbour, Kitron, Piesman and Fish2010). The abundance of this tick was also reported to be related to deer abundance in previous studies (Daniels et al., Reference Daniels, Fish and Schwartz1993; Duffy et al., Reference Duffy, Campbell, Clark, Dimotta and Gurney1994; Stafford et al., Reference Stafford, Denicola and Kilpatrick2003).

Snow depth showed a negative association with abundance at all tick stages. Snow cover appears to act as insulation and thus protects against low temperature and desiccation, giving a positive effect on the abundances of other tick species (McEnroe, Reference McEnroe1984; Drew & Samuel, Reference Drew and Samuel1986). However, in this study, no tick sampling was conducted during times of snow cover, and the opposite effect was observed. A possible explanation is that sika deer are limited to less snowy areas due to the heavy load on the hooves (Takatsuki, Reference Takatsuki1992; Igota et al., Reference Igota, Sakuragi, Uno, McCullough, Kaji and Takatsuki2008; Yabe & Takatsuki, Reference Yabe, Takatsuki, McCullough, Kaji and Takatsuki2009). Although larval abundance is not directly affected by deer density, snow depth might have an indirect effect through its effects on deer or other host species distributions.

H. longicornis is a known vector of Japanese spotted fever and bovine theileriosis. In the case of Japanese spotted fever, the causative agent, Rickettsia japonica (Alphaproteobacteria: Rickettsiaceae), infects all stages of H. longicornis. However, in a review of 101 cases of human infestation with this tick reviewed by Okino et al. (Reference Okino, Ushirogawa, Matoba and Hatsushika2008), most cases were caused by adults (n = 83), while larvae and nymphs were responsible for only 10 and 2 cases, respectively. In the case of bovine theileriosis, only adults and nymphs can transmit the causal agent, Theileria sergenti (Aconoidasida: Theileriidae), because transovarial transmission of T. sergenti does not occur in this tick species.

Our GLMM analysis showed that a high density of sika deer along with specific biological and abiotic factors could increase the risk of vector-borne Japanese spotted fever and bovine theileriosis through increased local abundance of questing H. longicornis nymphs and adults, but not larvae. Therefore, the positive relationship between the abundance of questing adult and nymphal H. longicornis and the density of sika deer revealed in our study suggests that the recent expansion of sika deer in Japan is directly correlated with the expansion of Japanese spotted fever and bovine theileriosis. Sika deer can serve as an amplifier host of the tick carrying these diseases.

The supplementary material for this article can be found at http://www.journals.cambridge.org/BER

Acknowledgements

Dr Hiromi Fujita and Dr Hideo Ootake kindly provided advice on tick species identification. This study was supported by a Grant-in-Aid for Research and Development Projects for Application and Promoting New Policy of Agriculture, Forestry and Fisheries (grant number 1910, 2007–2009) by the Ministry of Agriculture Forestry and Fisheries, Japan.

References

Atwood, E.L. & Sonenshine, D.E. (1967) Activity of American dog tick Dermacentor variabilis (Acarina – Ixodidae) in relation to solar energy changes. Annals of the Entomological Society of America 60, 354362.CrossRefGoogle ScholarPubMed
Benjamin, M.A., Zhioua, E. & Ostfeld, R.S. (2002) Laboratory and field evaluation of the entomopathogenic fungus Metarhizium anisopliae (Deuteromycetes) for controlling questing adult Ixodes scapularis (Acari: Ixodidae). Journal of Medical Entomology 39, 723728.Google Scholar
Burnham, K.P. & Anderson, D.R. (2002) Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. New York, Springer-Verlag.Google Scholar
Chilton, N.B. & Bull, C.M. (1993) A comparison of the off-host survival times of larvae and nymphs of two species of reptile ticks. International Journal for Parasitology 23, 693696.Google Scholar
Daniels, T.J., Fish, D. & Schwartz, I. (1993) Reduced abundance of Ixodes scapularis (Acari: Ixodidae) and Lyme disease risk by deer exclusion. Journal of Medical Entomology 30, 10431049.Google Scholar
Diuk-Wasser, M.A., Vourc'h, G., Cislo, P., Hoen, A.G., Melton, F., Hamer, S.A., Rowland, M., Cortinas, R., Hickling, G.J., Tsao, J.I., Barbour, A.G., Kitron, U., Piesman, J. & Fish, D. (2010) 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 Ecology and Biogeography 19, 504514.Google Scholar
Drew, M.L. & Samuel, W. (1986) Reproduction of the winter tick, Dermacentor albipictus, under field conditions in Alberta, Canada. Canadian Journal of Zoology 64, 714721.Google Scholar
Duffy, D.C., Campbell, S.R., Clark, D., Dimotta, C. & Gurney, S. (1994) Ixodes scapularis (Acari: Ixodidae) deer tick mesoscale populations in natural areas: effects of deer, area, and location. Journal of Medical Entomology 31, 152158.Google Scholar
Elias, S.P., Lubelczyk, C.B., Rand, P.W., Lacombe, E.H., Holman, M.S. & Smith, R.P. (2006) Deer browse resistant exotic-invasive understory: an indicator of elevated human risk of exposure to Ixodes scapularis (Acari: Ixodidae) in southern coastal Maine woodlands. Journal of Medical Entomology 43, 11421152.Google Scholar
Fernandes, E.K.K., Bittencourt, V.R.E.P. & Roberts, D.W. (2012) Perspectives on the potential of entomopathogenic fungi in biological control of ticks. Experimental Parasitology 130, 300305.Google Scholar
Fujimoto, K. (1988) Ecological studies on ixodid ticks : 5. The effect of humidity on the oviposition and development of Haemaphysalis longicornis and H. flava (Acarina : Ixodidae). Japanese Journal of Sanitary Zoology 39, 2733 (in Japanese with English summary).Google Scholar
Fujimoto, K. (1993) Tolerances of Haemaphysalis longicornis, H. flava, and Ixodes ovatus (Acari : Ixodidae) to high and low temperatures in relation to distribution patterns of these ticks in the Chichibu Mountains. Japanese Journal of Environmental Entomology and Zoology 5, 158165.Google Scholar
Fujimoto, K. (1995) Effect of photoperiod on the host attachment and development in the tick, Haemaphysalis longicornis . Japanese Journal of Sanitary Zoology 46, 345348.Google Scholar
Fujimoto, K. (2000) The host-feeding activity of Haemaphysalis longicornis nymphs (Acari : Ixodidae) in autumn observed under experimental conditions. Japanese Journal of Sanitary Zoology 51, 105107.Google Scholar
Fujimoto, K. (2004) The host-feeding activity of Haemaphysalis longicornis larvae (Acari: lxodidae) in autumn observed under experimental conditions, with special reference to behavioral diapause. Bulletin of Saitama Medical School Junior College 15, 2931.Google Scholar
Fujimoto, K. (2008) An experimental analysis on the autumn activity of the new generation of Haemaphysalis longicornis nymphs (Acari: Ixodidae) : the autumn activities of the nymphs that developed outdoors in August and September. Japanese Journal of Sanitary Zoology 59, 1517.Google Scholar
Fujisaki, K., Kitaoka, S. & Morii, T. (1975) Effect of different combinations of temperature and humidity on the oviposition of Haemaphysalis longicornis and Boophilus microplus . Bulletin of the National Institute of Animal Health 70, 2734.Google Scholar
Fujita, H. & Takada, N. (2007) Identification of immature ticks in Japan. pp. 5368 in Organizing Committee of SADI (Ed.) Acari and Emerging Reemerging Infectious Disease. Tokyo, Zenkoku Noson Kyoiku Kyokai Publishing. (in Japanese).Google Scholar
Heath, A.C. (1981) The temperature and humidity preferences of Haemaphysalis longicornis, Ixodes holocyclus and Rhipicephalus sanguineus (Ixodidae): studies on engorged larvae. International Journal for Parasitology 11, 169175.Google Scholar
Hoogstraal, H., Roberts, F.H., Kohls, G.M. & Tipton, V.J. (1968) Review of Haemaphysalis (Kaiseriana) longicornis Neumann (resurrected) of Australia, New Zealand, New Caledonia, Fiji, Japan, Korea, and Northeastern China and USSR, and its parthenogenetic and bisexual populations (Ixodoidea, Ixodidae). Journal of Parasitology 54, 11971213.Google Scholar
Igota, H., Sakuragi, M. & Uno, H. (2008) Seasonal migration of sika deer on Hokkaido Island, Japan. pp. 251272 in McCullough, D., Kaji, K. & Takatsuki, S. (Eds) Sika Deer: Biology and Management of Native and Introduced Populations. Tokyo, Springer.Google Scholar
Inglis, G.D., Goettel, M.S., Butt, T. & Strasser, H. (2001) Use of hyphomycetous fungi for managing insect pests. pp. 2369 in Butt, T.M., Jackson, C. & Magan, N. (Eds) Fungi as Biocontrol Agents: Progress Problems and Potential. Oxford, CAB International.Google Scholar
Inokuma, H., Fujimoto, T., Hosoi, E., Tanaka, S., Fujisaki, K., Okuda, M. & Onishi, T. (2002) Tick infestation of sika deer (Cervus nippon) in the western part of Yamaguchi Prefecture, Japan. Journal of Veterinary Medical Science 64, 615617.Google Scholar
Inokuma, H., Tsuji, T., Kim, S.J., Fujimoto, T., Nagata, M., Hosoi, E., Arai, S., Ishihara, C., Okuda, M. (2004) Phylogenetic analysis of Theileria sp. from sika deer, Cervus nippon, in Japan. Veterinary Parasitology 120, 339345.Google Scholar
Inokuma, H., Seino, N., Suzuki, M., Kaji, K., Takahashi, H., Igota, H. & Inoue, S. (2008) Detection of Rickettsia helvetica DNA from peripheral blood of Sika deer (Cervus nippon yesoensis) in Japan. Journal of Wildlife Diseases 44, 164167.Google Scholar
Inokuma, Y. & Yokoyama, N. (2010) Study on factors to decrease risk for Theileria orientalis infection of cattle in pastureland in Eastern Hokkaido. Japanese Journal of Clinical Bovine Parasitology 1, 1317 (in Japansese with English Summary).Google Scholar
Ishihara, T. (1968) Bovine piroplasmosis in Japan. Japanese Agricultural Research Quarterly 3, 2354.Google Scholar
Iwamoto, T., Sakata, T., Nakazono, T., Kaoka, H., Ikeda, K., Nishishita, Y., Tokida, K. & Doi, T. (2000) Improvement of the pellet count method for the estimation of sika deer density. Mammalian Science 40, 117 (in Japanese with English summary).Google Scholar
Izumiyama, S. & Mochizuki, T. (2008) Seasonal range use of Sika deer which inhabits the sub-alpine zone in the Southern Japan Alps. Bulletin Shinshu University Alpine Field Center 30, 2532 (in Japanese with English summary).Google Scholar
Japan Meteorological Agency (2002) Mesh Climatic Data 2000. Tokyo, Japan Meteorological Business Support Center.Google Scholar
Jensen, P. & Jespersen, J. (2005) Five decades of tick-man interaction in Denmark - an analysis. Experimental and Applied Acarology 35, 131146.CrossRefGoogle ScholarPubMed
Johnson, J.B. & Omland, K.S. (2004) Model selection in ecology and evolution. Trends in Ecology and Evolution 19, 101108.Google Scholar
Jouda, F., Perret, J.L. & Gern, L. (2004) Ixodes ricinus density, and distribution and prevalence of Borrelia burgdorferi sensu lato infection along an altitudinal gradient. Journal of Medical Entomology 41, 162169.Google Scholar
Kaaya, G.P. & Hassan, S. (2000) Entomogenous fungi as promising biopesticides for tick control. Experimental and Applied Acarology 24, 913926.CrossRefGoogle Scholar
Kamio, T., Fujisaki, K. & Minami, T. (1989) Correlation between the infection rate of the vector tick, Haemaphysalis longicornis and the parasitaemia of cattle infected with Theileria sergenti . Annals of Tropical Medicine and Parasitology 83, 7783.Google Scholar
Kitaoka, S., Morii, T. & Fujisaki, K. (1975) Consideration on the grazing cattle-tick relationship in Japan with special reference to large wild mammals and the deer-tick ecosystem in the Tanzawa Mountains. Bulletin of the National Institute of Animal Health 70, 3542 (in Japanese).Google Scholar
Mahara, F. (1997) Japanese spotted fever: report of 31 cases and review of the literature. Emerging infectious diseases 3, 105111.CrossRefGoogle ScholarPubMed
Maruyama, N. & Furubayashi, K. (1983) Preliminary examination of block count method for estimating numbers of sika deer in Fukadake. Journal of the Mammalogical Society of Japan 9, 274278.Google Scholar
McEnroe, W. (1984) The effect of snowcover on an American dog tick, Dermacentor variabilis (Say.)(Acari, Ixodidae) population under a harsh winter environment. Zeitschrift für Angewandte Entomologie 97, 481484.Google Scholar
Mori, K. & Tsunoda, T. (1996) Migration of Larvae in the Tick, Haemaphysalis longicornis NEUMANN (Acarina, Ixodidae). Japanese Journal of Environmental Entomology and Zoology 7, 211213 (in Japanese with English summary).Google Scholar
Mori, K., Tsunoda, T., & Fujimagari, M. (1995) Ixodid ticks on sika deer Cervus nippon Temminck in Chiba prefecture. Medical Entomology and Zoology 46, 313316 (in Japanese with English summary).Google Scholar
Motokuni, T., Tarumi, H., Miyamura, H., Kameide, N. & Sunaga, T. (1999) Analysis of DUV values of ground arrival ultra violet radiation based on weather classification: part2 seasonal change of daily accumulations of UV-B/DUV and relation to solar radiation. pp. 677678. in Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan. D-1, Environmental Engineering I. Hiroshima, Architectural Institute of Japan (in Japanese).Google Scholar
Needham, G.R. & Teel, P.D. (1991) Off-host physiological ecology of ixodid ticks. Annual Review of Entomology 36, 659681.Google Scholar
Ohdachi, S.D., Ishibashi, Y., Iwasa, M.A. & Saitoh, T. (2009) The Wild Mammals of Japan. Kyoto, Shoukadoh Book Sellers.Google Scholar
Okino, T., Ushirogawa, H., Matoba, K. & Hatsushika, R. (2008) Bibliographical studies on human cases of hard tick (Acarina: Ixodidae) bites in Japan (2) Cases of Haemaphysalis longicornis and H. fulava infestation. Kawasaki Medical Journal 34, 185201 (in Japanese with English summary).Google Scholar
Oorebeek, M. & Kleindorfer, S. (2008) Climate or host availability: what determines the seasonal abundance of ticks? Parasitology Research 103, 871875.Google Scholar
Ostfeld, R.S., Price, A., Hornbostel, V.L., Benjamin, M.A. & Keesing, F. (2006) Controlling ticks and tick-borne zoonoses with biological and chemical agents. BioScience 56, 383394.Google Scholar
Rand, P.W., Lubelczyk, C., Lavigne, G.R., Elias, S., Holman, M.S., Lacombe, E.H. & Smith, R.P. (2003) Deer density and the abundance of Ixodes scapularis (Acari: Ixodidae). Journal of Medical Entomology 40, 179184.Google Scholar
R Development Core Team (2011) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, Available online at http://www.R-project.org/.Google Scholar
Ruiz-Fons, F. & Gilbert, L. (2010) The role of deer as vehicles to move ticks, Ixodes ricinus, between contrasting habitats. International Journal for Parasitology 40, 10131020.Google Scholar
Shimizu, G., Yoshizawa, K., Egashira, K., Hidaka, T. & Yamamoto, T. (2006) Measurements of solar UV radiation for 1996–2009. in Pref. Yamanashi. Annual Report of the Yamanashi Institute for Public Health 50, 6166 (in Japanese).Google Scholar
Skaug, H., Fournier, D., Nielsen, A., Magnusson, A. & Bolker, B. (2012) Generalized Linear Mixed Models using AD Model Builder. R package version 0.7.2.12. Available online at http://glmmadmb.r-forge.r-project.org.Google Scholar
Stafford, K.C., Denicola, A.J. & Kilpatrick, H.J. (2003) Reduced abundance of Ixodes scapularis (Acari : Ixodidae) and the tick parasitoid Ixodiphagus hookeri (Hymenoptera : Encyrtidae) with reduction of white-tailed deer. Journal of Medical Entomology 40, 642652.Google Scholar
Sutherst, R. & Bourne, A. (2006) The effect of desiccation and low temperature on the viability of eggs and emerging larvae of the tick, Rhipicephalus (Boophilus) microplus (Canestrini) (Ixodidae). International Journal for Parasitology 36, 193200.Google Scholar
Sutherst, R.W. & Bourne, A.S. (1991) Development, survival, fecundity and behaviour of Haemaphysalis (Kaiseriana) longicornis (Ixodidae) at two locations in southeast Queensland. International Journal for Parasitology 21, 661672.Google Scholar
Tagliapietra, V., Rosà, R., Arnoldi, D., Cagnacci, F., Capelli, G., Montarsi, F., Hauffe, H. & Rizzoli, A. (2011) Saturation deficit and deer density affect questing activity and local abundance of Ixodes ricinus (Acari, Ixodidae) in Italy. Veterinary Parasitology 183, 114124.Google Scholar
Takatsuki, S. (1992) Foot morphology and distribution of Sika deer in relation to snow depth in Japan. Ecological Research 7, 1923.Google Scholar
Takatsuki, S. (2009) Effects of sika deer on vegetation in Japan: a review. Biological Conservation 142, 19221929.Google Scholar
Takatsuki, S. & Maeyama, E. (2007) What is happening in the sika deer habiat in western Tokyo: degradation of Okutama Mountains. Biological Science 59, 5056 (in Japanese).Google Scholar
Takatsuki, S., Suzuki, K. & Higashi, H. (2000) Seasonal elevational movements of sika deer on Mt. Goyo, northern Japan. Mammal Study 2, 107114.Google Scholar
Tsukada, H. (2011) Effect of grazing on habitat and behavior of wild mammals. Journal of Integrated Field Science 8, 4954.Google Scholar
Tsunoda, T. (2007 a) Seasonal change in habitat use by Haemaphysalis longicornis (Acari: Ixodidae): plant, litter, and soil. Journal of Medical Entomology 44, 575579.Google Scholar
Tsunoda, T. (2007 b) Ticks of the Tanzawa Mountains. pp. 360362 in The Research Group of the Tanzawa Mountains (Ed.) Results of the Sceinetific Research on the Tanzawa Mountains. Sagamihara, Hiraoka Environmental Science Laboratory (in Japanese).Google Scholar
Tsunoda, T. (2007 c) Interspecific and intraspecific associations of two species of hard ticks, Haemaphysalis longicornis and Haemaphysalis megaspinosa, in relation to questing site. Journal of Parasitology 93, 531541.Google Scholar
Tsunoda, T. (2008) Influence of aggregation and relative humidity on longevity of unfed bush tick, Haemaphysalis longicornis Neumann (Acari: Ixodidae). Journal of Parasitology 94, 990992.Google Scholar
Tsunoda, T. & Tatsuzawa, S. (2004) Questing height of nymphs of the bush tick, Haemaphysalis longicornis, and its closely related species, H. mageshimaensis: correlation with body size of the host. Parasitology 128, 503509.Google Scholar
Wildlife Management Office (2009) A Report on the Present Population Status of Sika Deer in Kanagawa Prefecture in 2008. Kanagawa, Wildlife Management Office (in Japanese).Google Scholar
World Health Organization (2002) Global Solar UV Index: A Practical Guide. Geneva, World Health Organization.Google Scholar
Yabe, T. & Takatsuki, S. (2009) Migratory and sedentary behavior patterns of sika deer in Honshu and Kyushu, Japan. pp. 273283 in McCullough, D., Kaji, K. & Takatsuki, S. (Eds) Sika Deer: Biology and Management of Native and Introduced Populations. Tokyo, Springer.Google Scholar
Yamaguchi, N., Tipton, V., Keegan, H. & Toshioka, S. (1971) Ticks in Japan, Korea and Ryukyuislands. Brigham Young University Science Bulletin Biological Series 15, 1226.Google Scholar
Yamane, I., Nishiguchi, A., Kobayashi, S. & Zeniya, Y. (2006) Cross-sectional survey of ixodid tick species on grazing cattle in Japan. Experimental and Applied Acarology 38, 6774.Google Scholar
Yamane, I., Nakamura, Y. & Tsukada, H. (2009) A Report on National Survey for Cattle Grazing Fields. Tsukuba, National Institute of Animal Health (in Japanese).Google Scholar
Yamashiro, T. (2000) The propaties of an entomopathogenic fungus: the effects of ultraviolet and agricultural chemicals on sterillization of Beauveria bassiana . Tochigi-ken Sangyō Sentā kenkyū yōhō [Bulletin of Tochigi-ken Sericultural Center Resarch] 42, 5356 (in Japanese).Google Scholar
Yamauchi, T., Tabara, K., Kanamori, H., Kawabata, H., Arai, S., Katayuama, T., Fujita, H., Yano, Y., Takada, N. & Itagaki, A. (2009) Tick fauna associated with sika deer density in the Shimane peninsula, Honshu, Japan. Medical Entomology and Zoology 60, 297304.Google Scholar
Figure 0

Fig. 1. Locations of tick sampling sites in three regions of central Japan: Tochigi Prefecture, Izu Peninsula, and Kanagawa Prefecture. Sampling sites in each region are represented by the same symbol, and thin lines indicate the prefecture borders.

Figure 1

Table 1. Species and number of ixodid ticks collected by the flagging method in three regions of central Japan.

Figure 2

Fig. 2. Relationship between H. longicornis tick abundance by stage (larval, nymphal, and adult) and sika deer density at each tick-sampling site.

Figure 3

Table 2. Total number of H. longicornis larvae, nymphs, and adults collected by season

Figure 4

Table 3. AICc ranking of full models using biological and abiotic factors to estimate dependence of total number of H. longicornis larvae, nymphs, and adults in GLMMs. Null models including only the intercept are also shown for each tick stage.

Figure 5

Table 4. A parameter estimate of each factor relative to the best model (table 3) in the case of equally parsimonious models. Estimates of the relative importance of predictor variables (RVI) were calculated by summing the Akaike weights across the set of all models in which the variable occurs. Model structure IDs correspond to those in table 3.