Figures
Abstract
Background
Intensification of world trade is responsible for an increase in the number of alien species introductions. Human-mediated dispersal promotes not only introductions but also expansion of the species distribution via long-distance dispersal. Thus, understanding the role of anthropogenic pathways in the spread of invading species has become one of the most important challenges nowadays.
Methodology/Principal Findings
We analysed the invasion pattern of the pinewood nematode in China based on invasion data from 1982 to 2005 and monitoring data on 7 locations over 15 years. Short distance spread mediated by long-horned beetles was estimated at 7.5 km per year. Infested sites located further away represented more than 90% of observations and the mean long distance spread was estimated at 111–339 km. Railways, river ports, and lakes had significant effects on the spread pattern. Human population density levels explained 87% of the variation in the invasion probability (P<0.05). Since 2001, the number of new records of the nematode was multiplied by a factor of 5 and the spread distance by a factor of 2. We combined a diffusion model to describe the short distance spread with a stochastic, individual based model to describe the long distance jumps. This combined model generated an error of only 13% when used to predict the presence of the nematode. Under two climate scenarios (stable climate or moderate warming), projections of the invasion probability suggest that this pest could expand its distribution 40–55% by 2025.
Conclusions/Significance
This study provides evidence that human-induced dispersal plays a fundamental role in the spread of the pinewood nematode, and appropriate control measures should be taken to stop or slow its expansion. This model can be applied to Europe, where the nematode had been introduced later, and is currently expanding its distribution. Similar models could also be derived for other species that could be accidentally transported by humans.
Citation: Robinet C, Roques A, Pan H, Fang G, Ye J, Zhang Y, et al. (2009) Role of Human-Mediated Dispersal in the Spread of the Pinewood Nematode in China. PLoS ONE 4(2): e4646. https://doi.org/10.1371/journal.pone.0004646
Editor: Sean Rands, University of Bristol, United Kingdom
Received: November 4, 2008; Accepted: January 29, 2009; Published: February 27, 2009
Copyright: © 2009 Robinet et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: We greatly acknowledge support for this work from the National Natural Science Foundation of China (30730074 and 30525009), a French grant from the “Partenariat Hubert Curien of the programme de recherches avancees franco-chinois” (phc pra AHH610548C- INVACHINE ), and the EU projects ALARM GOCE-CT-2003-506675 (FP6 Integrated Project, Assessing large-scale environmental risks with tested methods; www.alarmproject.net [69]) and PRATIQUE KBBE-2007-212459 (FP7 Project, Enhancements of pest risk analysis techniques). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Intensification of world trade is a major factor explaining the increase in the number of exotic species introduced everywhere in the world [1]. Invertebrates, mostly insects and nematodes, represent the greatest proportion (∼90%) of the organisms intercepted every year in Europe [2]–[3]. Insects are also the dominant intercepted organisms entering the United-States, representing 73–84% of the total [4]. Invasions depend on biological and climatic factors, but another important factor that has been mostly neglected until now is the economic factor [5]. In China, the rapid increase in the number of introduced species since the 1970s is most likely due to economic growth first, and then, to a lesser extent, climatic factors [5]. While only a small fraction of these biological invaders are able to survive and establish a new population [6], some of them cause a dramatic impact in the newly invaded environment, resulting in great ecological damage and economic loss [7]–[8]. Understanding invasion processes for assessing the risk of biological invasions and improving the management of invasive species has thus become a major challenge [9]–[10].
In this paper, we will deal with the invasion processes of an insect-vectored nematode in China. This nematode constitutes a model organism because: (1) it is spreading very rapidly throughout China, largely affecting tree survival, (2) Chinese territory is large enough to study long distance dispersal, and (3) China is a hotspot of international trade since it entered the World Trade Organization on December 11, 2001.
The pinewood nematode, Bursaphelenchus xylophilus (Steiner & Buhrer) Nickle, is native to North America [11]. Since the beginning of the 20th century, it has been accidentally introduced into several Asian countries (Japan- 1905, China- 1982, Taiwan- 1985, Korea- 1988). Then, it reached Europe where it was discovered in Portugal in 1999 [12]–[13]. Despite an intensive containment program in the Setúbal Peninsula where it was restricted since its arrival, the pinewood nematode is at present expanding to other parts of Portugal [14]–[16]).
In China, the nematode was discovered for the first time in Nanjing in 1982. Another infested area was found at more than 1000 km away from the previous one, near Hong-Kong, in 1988, supposedly resulting from a separate introduction. Genetic diversity and phylogenetic analysis suggest that, with the exception of populations in Mingguang, nematode infestations probably originated from Nanjing [17]. The pinewood nematode is vectored by long-horned beetles, mainly the Japanese pine sawyer, Monochamus alternatus Hope (Coleoptera: Cerambycidae). Almost 70 tree species have been listed as susceptible to the pinewood nematode but Pinus spp. constitute the most susceptible host plants [18]–[23]. Some of these pine species, especially the most affected one, Masson pine (P. massoniana), are widely distributed throughout China [24].
Ecophysiological models have already been developed to describe the host-vector interaction and the spread of the nematode among pine stands in Japan [25]–[27]. These models accurately describe the transmission of the disease based on empirical data. Long distance dispersal has been considered in terms of the proportion of beetles with the greatest dispersal capabilities (average long distance dispersal: 1.8 km [26]). However, there was no research describing the effects of human-mediated dispersal at larger scales until now.
Two measures can be used to calculate a confidence interval of the dispersal rate: (1) the distance from the nearest neighbour, and (2) the distance from the introduction point (see [28]–[29] for the dispersal rate of Argentine ants). When individuals (from plant or animal species) disperse randomly, reaction-diffusion models are generally used to describe the natural expansion of the populations (e.g. [30]–[33]). Their mathematical properties are well-known [32], [34]–[35]: they can generate a travelling wave with a constant asymptotic speed , with D the diffusion coefficient and ε the growth rate [35]–[36].
Human activities, such as logging or trade which requires wooden packaging material, increase the risk for accidental transportation of infested materials and may be responsible for the rapid spread of the nematode. Therefore any connection between infested and non-infested areas (e.g., via highways, railways, rivers, or an electric power network with wooden poles) probably increases the risk of invasion (GF Fang, pers. com.). Transportation hubs have been reported to play an important role in occasional long-distance, human mediated dispersal [37]. Thus, we also investigated the role of anthropogenic pathways and we developed a model combining short- and long-distance spread to gain a greater understanding of the dispersal mechanism and predict any future range expansion.
Human population density could explain a large part of long-distance dispersal for some insect species, as has been shown for the horse chestnut leafminer, Cameraria ohridella Deschka & Dimic, in Europe [38]–[39], and the emerald ash borer, Agrilus planipennis Fairmaire, in North America [40]. Because the movements of people, cars and trucks also increase the probability of moving infected beetles or infested wood, we carefully investigated the effect of the human population density as an indicator of human-mediated dispersal risk.
Many studies on the climatic tolerance of the carrier beetle have been conducted in Japan and China. Two main thresholds associated with the beetle survival and distribution have been determined: (1) the mean air temperature in July should be above 21.3°C, and (2) the mean temperature in January should be above −10°C [41]–[42]. Moreover, the pinewilt disease has never occurred in North America or Japan when the mean air temperature of the warmest month is lower than 20°C [11]. Territorial expansion of alien species often results not only from human-mediated dispersal but also from climate change [43]. The distribution of many plant and animal species is already affected by climate change (e.g. [44]–[47]). A temperature rise of 2°C could considerably increase the risk of the pine wilt disease development in north western Spain [48]. Consequently, we tested the potential effects of a temperature increase on the spread of this nematode in China.
Methods
Data description
Pinewood nematode.
Data was collected by the General Forest Pest Control Station (Shenyang, China), the administration in charge of managing invasive species in China. Two data sets were analyzed: (1) a history of invasion between 1982 and 2005: a list of 171 locations invaded by the pinewood nematode and the year of the initial observation for each location (Fig. 1A), and (2) a monitoring data set summarizing the number of dead trees attributed to the pinewood nematode over 15 years at 7 locations (Table 1). Random samples were collected with Baermann funnels and the samples were examined microscopically. We attributed the death of a tree to the pinewood nematode according to the following typical symptoms: browning needles, many sawyer oviposition scars, declining of resin excretion, and blue xylem stains.
The pinewood nematode invasion in China between 1982 and 2005 (Panel A). The first pinewood nematode observation in China (in Nanjing) is represented by a white star. Spatial distribution of 10 susceptible tree species (Panel B). Black dots represent locations already infested. Potential anthropogenic pathways (Panel C), and spatial distribution of human-population density in 2000 (Panel D).
Host trees.
We considered the 10 native susceptible species with the largest distributions in China (Fig. 1B [24], [49]–[50]): Pinus armandii Franch, Pinus bungeana Zucc. ex Endl., Pinus densiflora Sieb. & Zucc., Pinus kesiya Royle ex Gordon var. langbianensis (A. Chev.) Gaussen ex N.-S. Bui, Pinus koraiensis Sieb. & Zucc., Pinus massoniana Lamb., Pinus sylvestris L. var. mongolica Litv., Pinus tabuliformis Carr., Pinus taiwanensis Hayata, and Pinus yunnanensis Franch. We used GIS maps to represent the ranges of these susceptible tree species. The maps were provided by the Institute of Geographic Sciences and Natural Resources Research (Chinese Academy of Sciences, 11A Datun Road, Chaoyang District, Beijing 100101, People's Republic of China).
Potential anthropogenic pathways and human-population density.
To investigate the role of anthropogenic pathways, our study utilized the GIS maps of rivers, river ports, lakes, railways (Fig. 1C), and the spatial distribution of human-population density in 2000 (Fig. 1D). These data were also provided by the Institute of Geographic Sciences and Natural Resources Research (Beijing).
Data analysis and models
Climate suitability.
The mean temperature in January (TJan) and the mean temperature in July (TJul) were estimated from a multiple regression using latitude, longitude and elevation as predictive variables, and the data of 279 weather stations in China, over 1951–1980, were used to fit the regression model [51] (data were provided by the Department of Geography, Mainz University, Germany, http://www.webgis-china.de/). Then, these temperatures were interpolated on a grid cell (720×600 cells, with a cell size Δx = Δy = 0.04°, also used hereafter for simulations of the model) covering our study area (20°00–45°00 N, 95°00–125°00 E). For this purpose, we used a digital elevation model derived from the USGS/NASA SRTM data [52]. Based on the developmental thresholds previously defined for the carrier beetle, we determined the areas where the climate was suitable (TJul≥21.3°C and TJan≥−10°C) and we assumed that, once arrived at a given location, the pinewood nematode could establish only if the climate was suitable, otherwise the invasion would fail. To explore the effects of global warming, we also tested the effects of a temperature increase (see the model description for further details).
Short-distance dispersal model.
We developed a linear regression between the distance spread and years using the invasion dataset from the area of the first introduction point in China (up to 100 km from Nanjing). To assess short-distance dispersal ability, we considered the spread from this single point because, in other locations, migrants could come from many sources and it was not possible to distinguish short and long distance dispersal. Then, we applied a reaction-diffusion model to describe this short-distance dispersal:(1)where N is the population density of the nematode being carried by the beetles, t the time, x and y the geographic coordinates, D the diffusion coefficient and ε the growth rate.
To estimate the parameter ε, we calculated the growth rate for each year, at the 7 locations using the monitoring data set. For this, we assumed that the ratio of dead trees was a good indicator of the pinewood nematode growth rate [53]. Thus, we calculated εi = ln( Ti / Ti–1), where Ti is the number of dead trees in year i, at a given location. Based on this value and the observed rate of population expansion, we determined the value of the diffusion coefficient D using the formula .
Long-distance dispersal model.
Since the source of the new populations was uncertain, we identified the points which could potentially be invaded via short-distance spread. They were located at a distance that could be traversed in one generation (i.e. one year) from previously infested areas. That distance was estimated from the observed expansion rate at one of the original infestation areas near Nanjing. We thus quantified the relative importance of short-distance and long-distance dispersal. Then, to analyze long-distance dispersal in-depth, we used two methods to estimate the distance spread: (1) distance from the nearest neighbour, and (2) distance from the introduction point. Since there were probably two introduction points (in Nanjing area in 1982, and in Hong-Kong area in 1988), we calculated the distance from both introduction points (if the area was infested after 1988) and selected the shortest distance. For the locations infested before 1988, we considered only the distance from the first introduction point in Nanjing area.
Then, we investigated the effects of potential anthropogenic pathways such as rivers, river ports, lakes, railways and the human population density on the long-distance dispersal pattern. Since it was difficult to determine the appropriate spatial scale to detect the effects of these pathways, we considered three neighbourhoods around each invasion point: N120, a grid cell composed of 121×121 cells (60 cells on the left, on the right, above and below the invasion point, representing ca. 245.103 km2), N60, a grid cell composed of 61×61 cells (30 cells in each direction, representing ca. 62.103 km2), and N30, a grid cell composed of 16×16 cells (15 cells in each direction, representing ca. 16.103 km2). We calculated the proportion of cells occupied by each pathway (river, railway and lake), the number of ports and the mean human population density in these neighbourhoods. We made these calculations using the invasion dataset and also a second dataset (of the same size, n = 156) selected randomly over China in order to test whether the spread patterns were significantly different from what we would expect at random. We logarithmically transformed these numbers and proportions, and used the Jarque Bera test [54] to check whether the data were normally distributed. Since most of them were not normally distributed, we used a non parametric test (Wilcoxon test) to compare both datasets (observations and randomly generated dataset).
We also calculated the correlation between spread distance and the human population density in the infested area, and we calculated the invasion probability as a function of the human population density. Since the human population is heterogeneously distributed over China, we calculated a corrected invasion probability taking into account the frequency of human densities as follows: (1) we counted the number of infested areas in which the human population density was in the interval [2500 n, 2500 n+2500[, for n = 0, …, 15; (2) we counted the number of cells (from the grid cell used for interpolating temperature data) in which the human population in China was in the same interval; (3) we divided the first number by the second, and expressed it finally as a proportion. Thus, the probability of the invasion of a given area was adjusted with regard to the frequency distribution of human population density in China.
To define the long-distance dispersal kernel, the previous probability (associated with the human population density) was multiplied by the probability to disperse at a given distance. Since long-distance dispersal was noticeably different before and after 2001, we applied a Gaussian kernel with different parameters. These parameters were estimated using the least-squares method.
Combination of short and long distance dispersal models.
To explicitly model the pinewood nematode expansion, short- and long-distance spread models were combined. Based on the infestation pattern around Nanjing, we determined the number of years, ND, an area should be infested before it can provide long-distance dispersers and we included this delay in the model.
Invasion of the Hong-Kong area was first observed in 1988 with 6 infested locations, therefore the initial introduction probably took place much earlier. In 1988, infested sites were concentrated in a 27.3 km diameter region. Based on both the observed short-distance dispersal rate of the beetles (estimated around Nanjing) and the infested area around Hong-Kong reported in 1988, we estimated the year of first introduction around Hong-Kong.
A diffusion model (Eq. 1) was applied to describe short distance dispersal, with parameters D and ε estimated in a previous section. The population density was rounded to 0 in areas where climate was not suitable. The cell was considered infested when the population density reached 1.
Then, we randomly chose the number of long distance dispersers using a normal distribution with the previously determined mean and standard deviation. Parameters of this normal law were different before and after 2001. In reality, even if the main pattern is characterized by a sudden increase of invasions since 2001, the number of long distance dispersers would probably increase with the number of infested sites. Though, to avoid longer computation time, we assumed that the number of dispersers was independent of the number of infested sites.
For each disperser, we randomly chose a cell, infested for at least 6 years, from which it could emigrate, and we randomly selected the cell in which the disperser can settle following the long-distance dispersal kernel (Gaussian kernel multiplied by the probability to invade a location given the human population density, and then normalized). If the climate was not suitable in this cell, we assumed that this invasion failed. Distribution of host trees was not considered in the model because locations of infested sites were not correlated with host tree distribution.
To test the effects of an increase in temperature on the nematode spread, we made two types of simulations: (1) temperatures remain constant ( = mean temperatures over 1951–1980), and (2) temperatures increase linearly. According to climate model projections, global surface temperature will probably rise by 1.1°C–6.4°C between 1980–1999 and 2090–2099 [55]. Consequently, we tested a linear warming of 0.03°C per year, which is close to the mean predicted warming. For simplification, no stochasticity was introduced in this climate scenario.
We made 300 replicate simulations from 1982 to 2005, and we calculated the invasion probability in each cell as the number of times the model predicted the infestation of the cell, divided by 300. The spatial distribution of the invasion probability was compared to the pinewood nematode observations in 2005 to evaluate the fitting success of the model. An independent dataset is generally required to assess the model performance, but due to the complex dispersal processes, it was not possible to select a limited subset of this data and use the model to produce meaningful results. To complete this study, a projection for 2025 was calculated, with and without a temperature increase, using the infested locations (observations) up to 2005 as the initial population in the model. We made 300 replicate simulations for each model (with and without climate warming).
All these analyses and simulations were performed using R language [56]. The coastline extractor from the National Geophysical Data Center, NOAA Satellite and Information Service (http://www.ngdc.noaa/ngdc.htlm), provided a convenient file of China's border for the use in R. Maps were generated using ArcView 9.2, ESRI.
Results
Analysis of invasion data
The number of newly infested areas per year was significantly greater after 2001 (t test: t = 5.7031, df = 4.638, P = 0.003) (Fig. 1A). The mean number of newly infested areas was nearly five times the number of pre- 2001 infestations: the mean±SD was 3.89±3.18 before 2001, and 19.40±5.86 after. Of the 171 infested areas, only 93 (54%) were located within the distribution of susceptible trees documented here.
Climate data
Temperatures in both January and July were weakly correlated (R2<0.30). The following regressions explained 96% and 89% of the variation in the mean temperature in January (TJan) and the mean temperature in July (TJul), respectively:(2)
(3)where Lon is longitude (decimal degrees), Lat is latitude (decimal degrees) and Elev is elevation (m). All P-values were highly significant (F3,273 = 2100, P<0.001 for Eq.2, and F3,273 = 723.3, P<0.001 for Eq. 3).
Suitable climate area was consistent with the pinewood nematode observations (Fig. 2A). A large part in the south-eastern China is favourable for the nematode establishment. With a 3°C warming, the suitable area could expand by 40%, especially in the northern and north-western regions (Fig. 2A), but the Tibetan Plateau and the surrounding mountainous area, and also the northern part of Yunnan and the western part of Sichuan would remain quite unfavourable.
Potentially favourable areas according to normal temperatures over 1951–1980 (dark red), and 3°C of temperature increase (dark+light red is the potential expansion area) (Panel A). White dots represent locations already infested. Invasion probability predicted by the dispersal model: in 2005 (Panel B), in 2025 under the assumption of a stable climate (Panel C), and in 2025 under the assumption of a constant warming (+0.03°C/yr) (Panel D). White dots represent locations infested until 2005, and in Panel A, white dots with a black point inside represent infested locations where the predicted invasion probability in 2005 is zero.
Short-distance dispersal model
The analysis of the distance spread from the first introduction point with years revealed different patterns before and after 1987 (Fig. 3A). Indeed, the relationship between distance and years was linear during the first few years (Pearson's correlation = 0.98; t test: t = 15.2899, df = 8, P<0.001), while there was no such relationship later (Pearson's correlation = 0.10; t test: t = 0.4405, df = 19, P = 0.665). Between 1982 and 1987, the pinewood nematode invaded surrounding areas at a constant speed: c = 7.5 km/year. Furthermore, the range of maximum growth rates was ε = 0.51–2.17. When we considered the maximum value, the net reproductive rate was exp(ε) = 8.76. Using the formula with ε = 2.17 and c = 7.5 km per year, we found D = 6.480 km2/year, and we used these estimates in the short-distance dispersal model.
Distance spread (up to 100 km) from the first introduction point around Nanjing (Panel A). Corrected invasion probability as a function of human population density (Panel B). The grey line is the regression line. Grey dots were discarded from this analysis. Long-distance dispersal kernels (Panel C): observed probabilities over 1982–2005 (dash line), observed probabilities before and after 2001 (grey and black dots, respectively), and estimated probabilities before and after 2001 (grey and black lines, respectively).
Long-distance dispersal model
After eliminating local spread (spread within 7.5 km), 156 sites of the 171 appeared to result from long-distance dispersal. Thus, short-distance dispersal represented only 8.8% of the pinewood nematode dispersal. The mean number of long distance dispersals per year was close to the total number of dispersals, the mean±SD was 3.41±3.27 before 2001, and 18.2±5.89 after. These values were used as parameters of a normal law to estimate the number of long distance dispersers each year.
Estimates of long-distance dispersal were significantly higher with the second method (distance from an introduction point: mean = 338.62 km, SE = 23.85 km, median = 271.23 km, n = 150) in comparison with the first method (distance from the nearest neighbour: mean = 111.06 km, SE = 11.68 km, median = 53.83 km, n = 150) (paired t-test, log transformed, t = 18.73, d.f. = 149, P<0.001). Records around Hong-Kong in 1988 (n = 6) were excluded from this analysis because they were considered as a secondary introduction from abroad, independent from populations originating from Nanjing.
Using both methods, distance spread before 2001 was significantly lower than after this date (t-test, log-transformed, t = −3.80, d.f. = 143.18, P<0.001 for method 1, and t = −7.84, d.f. = 104.16, P<0.001 for method 2; the mean±SE were 73.39±10.38 km before 2001 and 135.48±12.09 km after, for method 1; and 171.17±13.46 km before 2001 and 447.18±24.96 km after, for method 2). The mean distance spread approximately doubled.
Effects of rivers were not significant whereas railways, ports, lakes and human population densities had a significant effect on the spread pattern (Wilcoxon test, P<0.001). Infested areas were more likely to be located nearby these pathways than points chosen at random (Table 2). The median human population density in areas where the pinewood nematode was present was three times larger than the median human population density in areas selected at random (ca. 3000 inhabitants per km2 vs 1000; Table 2).
No direct relationship was observed between the spread distance and the human population density. Spread distance was not correlated directly with the human population density, neither before 2001 (t-test, log-transformed, t = −1.22, d.f. = 57, P = 0.227 for method 1, and t = 0.03, d.f. = 57, P = 0.98 for method 2) nor after 2001 (t-test, log-transformed, t = 0.37, d.f. = 89, P = 0.71 for method 1, and t = 0.76, d.f. = 89, P = 0.45 for method 2). However, the probability that the pinewood nematode would invade a given area increased with the human population density when we introduced a correction related to the frequency of human population and discarded outliers (three points with a null invasion probability at more than 30000 inhabitants per km2) (Fig. 3B). Human population density (H) explained 87% of the variation in invasion probability (IP):(4)The P-value was highly significant (F1,11 = 76.83, P<0.001).
The following regression explained 99% of the variation in the dispersal probability (γLDD) before 2001 (with m = 0, σ = 252.62) and 88% after 2001 (with m = 377.22, σ = 240.45) (Fig. 3C):(5)All P-values were highly significant (F1,15 = 1032.8, P<0.001 before 2001, and F2,14 = 42.854, P<0.001 after 2001).
Combination of short and long distance dispersal models
The number of years an area should be infested before it can provide long-distance dispersers was estimated to ND = 6 years. Based on the invasion pattern near Hong-Kong (6 infested locations discovered in 1988, concentrated in a 27.3 km diameter region), and the results from natural dispersal (7.5 km per year), we estimated that the introduction took place ca. 4 years before, in 1984, and we selected one of the 6 locations to represent the second introduction point in China . As a result, in the model, we considered two introduction points from abroad: Nanjing in 1982 (32°05′ N, 118°48′ E) and around Hong-Kong in 1984 (23°06′ N, 113°26′ E).
The mean invasion probability predicted for locations already infested was 0.59. Only 22 invaded sites (12.9%) were located in a cell where the model predicted an invasion probability in 2005 of 0 (Fig. 2B). They were mainly located in the central part of China and on the coast line. Discarding these points, the mean invasion probability in infested locations was 0.68.
The model predicted an expansion of the potential invasion area (defined by P>0) by 47% in 2025 under a stable climate (Fig. 2C), and by 55% under climate warming (Fig. 2D). The invasion probability dramatically increases by 2025, and reached nearly 1 in areas where the invasion probability was predicted to be >0 in 2005. The most susceptible area (defined by P>0.95) could expand by 642% in 2025 with a stable climate, and by 662% with climate warming.
Discussion
The model presented here assesses the role of anthropogenic pathways in the spread of the pinewood nematode. The range expansion from the two source populations, Nanjing and Hong-Kong, was clearly identified, and according to the predicted probability of invasion, the pinewood nematode could become established in all of south-eastern China. The estimated spread rate of the beetle (7.5 km/year) was much higher than that previously observed in Japan (around 2 km/ year [57]–[58], and 4.2 km/year [26]). The net reproductive rate for the carrier beetles (exp(ε) = 8.76) was roughly of the same order than previously reported ( = 10.7 [27]).
The method used to determine an annual spread rate did not give accurate long distance estimates but, instead provided a reliable confidence interval (the minimum estimates resulting from method 1 and the maximum estimates resulting method 2). Our analysis revealed a clear difference in the invasion patterns before and after 2001. This breaking point coincided with the entry of China into the World Trade Organization. Although there is no evidence of cause-effect relationships, international trade has probably enhanced internal trade and wood exchange within China. In addition, survey techniques for monitoring pest species have been improved considerably in recent years (e.g. the Global Positioning System), and innovative methods have been developed to reliably detect new pinewood nematode infestations [59]. Improved methods of detection and monitoring may explain why many more infested areas have been discovered since 2000.
Errors in the model projections can be divided into three types: (1) errors in the accuracy of geo-referencing infested areas as well as the precise location of coast lines. These errors should diminish with improvements in geo-referencing methods and technology. (2) Errors resulting from undetected nematode introductions whether they originate as a result of maritime trade or from incursions from neighbouring countries like Myanmar. This problem points out the need for monitoring pest species on a global scale and the exchange of data among countries. (3) Errors in central China where the choice of dispersal kernels is compromised by anthropogenic activity in high population areas. For instance, we should investigate whether the construction of the gigantic Three Gorges Dam (30°49′N; 111°00′E) on the Yangtze River (which flows through Nanjing) has increased the risk of invasion in this central part of China due to intensified transportation of materials.
Climate warming would have little effect on the potential distribution of the pinewood nematode and is negligible when compared to the effects of human-mediated dispersal. In the model, we assume that temperatures will increase linearly and at the same rate throughout China, but important local variability in temperature increase is likely to occur among bioclimatic regions of China [60] (Beijing Climate Center, http://bcc.cma.gov.cn/). This warming was applied to both of the climate variables considered in our model, mean temperature in January (TJan) and mean temperature in July (TJul) even though winter and summer temperatures may not increase at the same rate [61]. Here we focussed on the effects of an increase in the mean temperature because temperature thresholds had previously been defined [41]–[42]. Other components of climate change will probably affect the pinewood nematode and the carrier beetle. For example, changes in mean precipitation, frequent heat-waves, droughts, floods and storms could potentially play an important role in the population abundance and distribution of many species [62]. Various climatic variables and their combinations should be considered to rigorously assess the effects of climate change on insect populations. Even if all the assumptions we have made are not realistic, they may help to simplify the problem and more detailed models could be developed in the future to more accurately determine the effects of climate change.
Rivers have no apparent effect on nematode spread. The Chinese riparian system is very large and uniformly distributed (Fig. 1C), so it might be difficult to determine its effects using our method. Apart from rivers, other pathways could be involved in the long distance dispersal. The number of human inhabitants is probably a good indicator of the risk for accidental transportation because there is also a good correlation between the human population density and the exchange of wood and wood products [63].
Distribution of susceptible trees cannot, by itself, explain the invasion pattern of the pinewood nematode in China. In this study, we considered the natural distribution of 10 susceptible Pinus species which are probably the most widespread host trees. The observed discrepancy between host distribution and infested sites may be the result of nematode colonization of non-native pine plantations and human-made plantations. Their distribution was not included in our analysis. A reliable documentation of this distribution should improve the accuracy of our model predictions. Inclusion of host-tree density data would lead to much better predictions than a simple presence/absence indicator because an invasion may fail due to Allee effects when the density of susceptible trees is too low [25]. Moreover, with better estimates of host-tree density, we may be able to combine our dispersal model with the disease transmission model developed by Yoshimura et al. [25]. The resulting model may provide us with more reliable predictions. Because phytosanitary certificates are required for the transfer of wood between provinces, it would be interesting to include these records in our model as an additional variable in the data analysis.
This model was developed to document and quantify the spread of the pinewood nematode spread in China, but it may also have of great potential for predicting the invasion patterns in other areas. In July 1985, the European Plant Protection Organization placed the pinewood nematode on the A1 list of quarantined pests [64]. The risk of invasion of European countries is very high [65] and the disease could become a major threat in European Scots and maritime pine forests [13]. The nematode was discovered for the first time in Portugal in 1999, and containment measures were immediately applied in the Setúbal Peninsula where it was successfully isolated for a few years. Genetic analysis revealed that the nematode found in this region was probably introduced from eastern Asia [66]. Despite governmental actions, the pinewood nematode was recently detected several hundred kilometres from the original infestation source, probably due to human-mediated transportation [16]. A national eradication programme, called PROLUNP, was recently launched for controlling this pest [14], [15]. Consequently, our model could be useful in assessing the invasion risk in Portugal and, at a larger scale, in Europe, and in identifying the most probable invasion pathways so as to improve the management strategies.
Modelling long-distance dispersal is quite difficult due to strong stochasticities, but such predictive models are essential for improving risk assessment and eradication strategies. Gypsy moth, Lymantria dispar (L.), is an invasive species which is currently expanding its distribution in North America mainly due to human transportation. Its spread rate was considerably reduced due to the creation of a barrier zone [67]. Eradication strategies and quarantine have also been applied in North America to contain the spread of the emerald ash borer, Agrilus planipennis [40]. Our model suggests that the pinewood nematode will probably extend its distribution further into the central and northern provinces of China. Because human mediated dispersal is a fundamental factor in this expansion, authorities should intensively control exchange of wood and wood products between infested and non-infested areas. Various management strategies are available: quarantine to avoid new introductions, eradication to avoid establishment, and a barrier zone to slow the spread [68]. Control measures should be able to break the relationships among the host-tree, the pinewood nematode and the carrier beetle [13], for instance by eliminating or reducing vector beetle populations. Complete eradication of the beetles in a barrier zone is not necessary to control the nematode range expansion because of Allee effects [26]. If the population density is reduced under a certain threshold, then the range expansion speed should decrease and isolated colonies could eventually go extinct naturally.
Acknowledgments
Elevation data used to interpolate temperature over the study area were kindly provided by the International Centre for Tropical Agriculture (CIAT, http://www.ciat.cgiar.org/). We thank Dr. Bin Shao (Institute of Geographic Sciences and Natural Resources Research, Beijing, China) for providing GIS maps. We thank Dr. Andrew M. Liebhold (USDA Forest Service, WV, USA) for providing insightful comments on an earlier version of the manuscript, and Olivier Denux (INRA, France) and Regis Young (USDA Forest Service, WV, USA) for helpful technical assistance.
Author Contributions
Conceived and designed the experiments: AR JHS. Analyzed the data: CR. Wrote the paper: CR. Supervised the research: JHS AR. Developed the model: CR. Did invasion data collection (pinewood nematode): HP GF JY. Gave assistance for data collection (abiotic factors) and translation: YZ.
References
- 1. Westphal MI, Browne M, MacKinnon K, Noble I (2008) The link between international trade and the global distribution of invasive alien species. Biol Invasions 10: 391–398. DOI: 10.1007/s10530-007-9138-5.
- 2. Roques A, Auger-Rozenberg M-A (2006) Tentative analysis of the interceptions of nonindigenous organisms in Europe during 1995–2004. EPPO Bull 36: 490–496. DOI: 10.1111/j.1365-2338.2006.01049.x.
- 3.
Roques A, Rabitsch W, Rasplus JY, Lopez-Vaamonde C, Nentwig W, et al. (2009) Alien terrestrial invertebrates of Europe. In: Nentwig W, Hulme P, Pysek P, Vila M, editors. Handbook of alien species in Europe. Springer Netherlands. In press. DOI: 10.1007/978-1-4020-8280-1.
- 4. McCullough DG, Work TT, Cavey JF, Liebhold AM, Marshall D (2006) Interceptions of nonindigenous plant pests at US ports of entry and border crossings over a 17-year period. Biol Invasions 8: 611–630. DOI: 10.1007/s10530-005-1798-4.
- 5. Lin W, Zhou G, Cheng X, Xu R (2007) Fast economic development accelerates biological invasions in China. PLoS ONE 2(11): e1208. DOI:10.1371/journal.pone.0001208.
- 6. Williamson M, Fitter A (1996) The varying success of invaders. Ecology 77: 1661–1666. DOI: 10.2307/2265769.
- 7.
Wan FH, Guo JY, Wang DH (2002) Alien invasive species in China: current status, research development, management strategies and risk assessment frame. In: Wang DH, Jeffrey AM, editors. International workshop on biodiversity and management for alien species. Beijing: China Environmental Science Press. pp. 77–102.
- 8. Pimentel D, Zuniga R, Morrison D (2005) Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecol Econ 52: 273–288. DOI:10.1016/j.ecolecon.2004.10.002.
- 9. Andersen MC, Adams H, Hope B, Powell M (2004) Risk analysis for invasive species: general framework and research needs. Risk Anal 24: 893–900. DOI: 10.1111/j.0272-4332.2004.00487.x.
- 10. Hayes KR, Barry SC (2008) Are there any consistent predictors of invasion success? Biol Invasions 10: 483–506. DOI: 10.1007/s10530-007-9146-5.
- 11. Rutherford TA, Mamiya Y, Webster JM (1990) Nematode-induced pine wilt disease: factors influencing its occurrence and distribution. For Sci 36: 145–155.
- 12. Mota MM, Braasch H, Bravo MA, Penas AC, Burgermeister W, et al. (1999) First report of Bursaphelenchus xylophilus in Portugal and in Europe. Nematology 1: 727–734. DOI: 10.1163/156854199508757.
- 13. Suzuki K (2002) Pine wilt disease – a threat to pine forest in Europe. Dendrobiology 48: 71–74.
- 14.
Rodrigues J (2006) Eradication program for the pinewood nematode in Portugal. Pine wilt disease: a worldwide threat to forest ecosystems, International symposium. 10–14 July, Lisbon. Available: http://www.nemalab.uevora.pt. Accessed 2009 Jan 5.
- 15.
Vieira PR, Mota MM (2008) Pine wilt disease: a worldwide threat to forest ecosystems. Springer Netherlands.
- 16.
Rego A, Alves RP, dos Santos Carvalho F, de Campos JV, Simeõn AM, et al. (2008) Protocolo de colaboração (report of the Portuguese Ministry of Agriculture). Available: http://www.afn.min-agricultura.pt/portal/proteccao-florestal/prolunp/resource/ficheiros/noticias/Protocolo-PROLUNP-Centro-Final.pdf. Accessed 2009 Jan 19.
- 17. Sun J, Yang S-Y, Cui C-L, Zhang C-X, Lin M-S, et al. (2008) Possible transmission routes of Bursaphelenchus xylophilus in China based on molecular data. J Nanjing Agr Univ (in Chinese with English summary) 31: 55–60.
- 18. Gao TY (1987) The fitness analysis of the pinewood nematode in China. Plant Quarantine (in Chinese) 1: 215–219.
- 19. Wang YY (1988) A brief view of the pine wilt disease in China. For Pest Dis (in Chinese) 1: 45–46.
- 20. Yang BJ, Wang QL (1989) Distribution of the pinewood nematode in China and susceptibility of some Chinese and exotic pines to the nematode. Can J For Res 19: 1527–1530. DOI:10.1139/x89-232.
- 21.
Song YS (2001) Integrated management for the pine wilt disease in China. Shenyang: Center for Forest Pest Protection of the Forestry Department of Administration, Shengyang, China (in Chinese).
- 22. Wang MX, Cheng LC, Song YS (2001) Risk analysis of impact of pine wood nematode (Bursaphelenchus xylophilus) on Hunan forestry and ecological environment. For Pest Dis (in Chinese, English summary) 2: 42–45.
- 23. Wang F, Yu SF, Feng SM, Zhao YX, Xu ZH (2002) Risk assessment for pinewood nematode introduced to Yunnan. J Yunnan Agricult Univ (in Chinese, English summary) 17: 421–422.
- 24.
Fu LG, Li N, Mills RR (1999) Pinaceae in flora of China. Tome 4. In: Wu ZY, Raven PH, editors. Cycadaceae through Fagaceae. Beijing: Science Press & St Louis: Missouri Botanical Garden Press.
- 25. Yoshimura A, Kawasaki K, Takasu F, Togashi K, Futai K, et al. (1999) Modeling the spread of pine wilt disease caused by nematodes with pine sawyers as vector. Ecology 80: 1691–1702. DOI: 10.1890/0012-9658(1999)080[1691:MTSOPW]2.0.CO;2.
- 26. Takasu F, Yamamoto N, Kawasaki K, Togashi K, Kishi Y, et al. (2000) Modeling the expansion of an introduced tree disease. Biol Invasions 2: 141–150. DOI: 10.1023/A:1010048725497.
- 27. Togashi K, Shigesada N (2006) Spread of the pinewood nematode vectored by the Japanese pine sawyer: modelling and analytical approaches. Popul Ecol 48: 271–283. DOI: 10.1007/s10144-006-0011-7.
- 28. Suarez AV, Holway DA, Case TJ (2001) Patterns of spread in biological invasions dominated by long-distance jump dispersal: insights from Argentine ants. Proc Natl Acad Sci USA 98: 1095–1100. DOI: 10.1073/pnas.031573098.
- 29. Ward DF, Harris RJ, Stanley MC (2005) Human-mediated range expansion of Argentine ants Linepithema humile (Hymenoptera: Formicidae) in New Zealand. Sociobiology 45: 1–7.
- 30. Skellam JG (1951) Random dispersal in theoretical populations. Biometrika 38: 196–218. DOI: 10.1093/biomet/38.1-2.196.
- 31. Kareiva PM (1983) Local movement in herbivorous insects: applying a passive diffusion model to mark-recapture field experiments. Oecologia 57: 322–327. DOI: 10.1007/BF00377175.
- 32.
Shigesada N, Kawasaki K (1997) Biological invasions: Theory and Practice. Oxford: Oxford University Press.
- 33.
Murray JD (2002) Mathematical biology I: an introduction. Third edition. Interdisciplinary Applied Mathematics, 17. New York: Springer-Verlag.
- 34. Fisher RA (1937) The wave of advance of advantageous genes. Ann Eugenics 7: 355–369.
- 35. Kolmogorov AN, Petrovski IG, Piskonov NS (1937) Study of the diffusion equation with growth of the quantity of matter and its application to a biological problem. Bull Univ Moscow, Ser Int, Sec A 1: 1–25.
- 36. Aronson DG, Weinberger HF (1978) Multidimensional diffusion arising in population genetics. Adv Math 30: 33–58. DOI: 10.1016/0001-8708(78)90130-5.
- 37. Floerl O, Inglis GJ, Dey K, Smith A (2008) The importance of transport hubs in stepping-stone invasions. J Appl Ecol. DOI: 10.1111/j.1365-2664.2008.01540.x.
- 38. Gilbert M, Grégoire J-C, Freise J, Heitland W (2004) Long-distance dispersal and human population density allow the prediction of invasive patterns in the horse chestnut leafminer Cameraria ohridella. J Anim Ecol 73: 459–468. DOI 10.1111/j.0021-8790.2004.00820.x.
- 39. Gilbert M, Guichard S, Freise J, Grégoire J-C, Heitland W, Straw N, Tilbuty C, Augustin S (2005) Forecasting Cameraria ohridella invasion dynamics in recently invaded countries: from validation to prediction. J Appl Ecol 42: 805–813. DOI: 10.1111/j.1365-2664.2005.01074.x.
- 40. Muirhead JR, Leung B, van Overdijk C, Kelly DW, Nandakumar K, et al. (2006) Modelling local and long-distance dispersal of invasive emerald ash borer Agrilus planipennis (Coleoptera) in North America. Divers Distrib 12: 71–79. DOI: 10.1111/j.1366-9516.2006.00218.x.
- 41. Kobayashi F, Yamane A, Ikeda T (1984) The Japanese pine sawyer beetle as the vector of pine wilt disease. Annu Rev Entomol 29: 115–135. DOI: 10.1146/annurev.en.29.010184.000555.
- 42. Ma R-Y, Hao S-G, Kong W-N, Sun J-H, Kang L (2006) Cold hardiness as a factor for assessing the potential distribution of the Japanese pine sawyer Monochamus alternatus (Coleoptera: Cerambycidae). Ann For Sci 63: 449–456. DOI: 10.1051/forest:2006025.
- 43. Dukes JS, Mooney HA (1999) Does global change increase the success of biological invaders? Trends Ecol Evol 14: 135–139. DOI: 10.1016/S0169-5347(98)01554-7.
- 44. Parmesan C, Ryrholm N, Stefanescu C, Hill JK, Thomas CD, et al. (1999) Poleward shifts in geographical ranges of butterfly species associated with regional warming. Nature 399: 579–583. DOI:10.1038/21181.
- 45. Thuiller W, Lavorel S, Araújo MB, Sykes MT, Prentice IC (2005) Climate change threats to plant diversity in Europe. Proc Natl Acad Sci USA 102: 8245–8250. DOI: 10.1073/pnas.0409902102.
- 46. Araújo MB, Thuiller W, Parson RG (2006) Climate warming and the decline of amphibians and reptiles in Europe. J Biogeogr 33: 1712–1728. DOI: 10.1111/j.1365-2699.2006.01482.x.
- 47. Hickling R, Roy DB, Hill JK, Fox R, Thomas CD (2006) The distributions of a wide range of taxonomic groups are expanding polewards. Global Change Biol 12: 450–455. DOI: 10.1111/j.1365-2486.2006.01116.x.
- 48.
Pérez G, Díez JJ, Ibeas F, Pajares JA (2008) Assessing pine wilt disease risk under a climate change scenario in Northwestern Spain. In: Bravo F, LeMay V, Jandl R, von Gadow K, editors. Managing forest ecosystems: the challenge of climate change. Dordrecht: Springer. pp. 269–282. DOI: 10.1007/978-1-4020-8343-3_15.
- 49.
Wang Y-A (1987) Research on regionalization of plant species (2), research on regionalization of Pinus massoniana. Forestry Regionalization Office of State Forestry Administration, P.R. China. Technical report.
- 50.
Ying T-S, Chen M-L, Chang H-C (2004) Atlas of the Gymnosperms of China. Beijing: China Science and Technology Press. pp. 70–97.
- 51.
Domroes M, Peng G (1988) The climate of China. Berlin: Springer.
- 52.
Jarvis A, Reuter HI, Nelson A, Guevara E (2006) Hole-filled seamless SRTM data V3, International Centre for Tropical Agriculture (CIAT). Available: http://srtm.csi.cgiar.org. Accessed 2009 Jan 5.
- 53.
Yang B-J, Zhu K-G (1995) Epidemics and management of pine wilt disease in China (in Chinese, English summary). Beijing: China Forestry Publishing House.
- 54. Jarque CM, Bera AK (1980) Efficient tests for normality, homoscedasticity and serial independence of regression residuals. Econ Lett 6: 255–259. DOI: 10.1016/0165-1765(80)90024-5.
- 55.
Meehl GA, Stocker TF, Collins WD, Friedlingstein P, Gaye AT, et al. (2007) Global climate projections. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, et al., editors. Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press.
- 56.
R Development Core Team (2007) R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. ISBN 3-900051-07-0, URL: http://www.R-project.org.
- 57.
Ido N, Kobayashi K (1977) Dispersal of Monochamus alternatus. Studies on the control of pine wilt disease. Tokyo: Secretariat of Agriculture, Forestry and Fisheries Research Council, Ministry of Agriculture, Forestry and Fisheries. pp. 87–88.
- 58. Fujioka H (1993) A report on the habitat of Monochamus alternatus Hope in Akita prefecture. Bull Akita For Tech Cent 2: 40–56.
- 59. Zhao L, Wei W, Liu X, Kang L, Sun J (2007) A novel rapid sampling method for pinewood nematode, Bursaphelenchus xylophilus (Nematoda: Parasitaphelenchidae). Can J For Res 37: 1867–1872. DOI: 10.1139/X07-049.
- 60.
Christensen JH, Hewitson B, Busuioc A, Chen A, Gao X, et al. (2007) Regional climate projections. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, et al., editors. Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press.
- 61. Hu ZZ, Yang S, Wu R (2003) Long-term climate variations in China and global warming signals. J Geophys Res – Atmos 108: 4614. DOI: 10.1029/2003JD003651.
- 62. Easterling DR, Meehl GA, Parmesan C, Changnon SA, Karl TR, et al. (2000) Climate extremes: observations, modeling, and impacts. Science 289: 2068–2074. DOI 10.1126/science.289.5487.2068.
- 63. Weber E, Li B (2008) Plant invasions in China: what is to be expected in the wake of economic development. Bioscience 58: 437–444. DOI: 10.1641/B580511.
- 64. Smith IM (1985) Pests and disease problems in European forests. FAO Plant Prot Bull 33: 159–164.
- 65. Evans HF, McNamara DG, Braasch H, Chadoeuf J, Magnusson C (1996) Pest risk analysis (PRA) for the territories of the European Union (as PRA area) on Bursaphelenchus xylophilus and its vectors in the genus Monochamus. EPPO Bull 26: 199–249. DOI: 10.1111/j.1365-2338.1996.tb00594.x.
- 66. Vieira P, Burgermeister W, Mota M, Metge K, Silva G (2007) Lack of genetic variation of Bursaphelenchus xylophilus in Portugal revealed by RAPD-PCR analyses. J Nematol 39: 118–126.
- 67. Sharov AA, Liebhold AM (1998) Model of slowing the spread of gypsy moth (Lepidoptera: Lymantriidae) with a barrier zone. Ecol Appl 8: 1170–1179. DOI: 10.1890/1051-0761(1998)008[1170:MOSTSO]2.0.CO;2.
- 68. Liebhold AM, Tobin PC (2008) Population ecology of insect invasions and their management. Annu Rev Entomol 53: 387–408. DOI: 10.1146/annurev.ento.52.110405.091401.
- 69. Settele J, Hammen V, Hulme P, Karlson U, Klotz S, et al. (2005) ALARM – Assessing LArge-scale environmental Risks for biodiversity with tested Methods. GAIA 14: 69–72.