www.nature.com/scientificreports
SCIENTIFIC
OPEN
Received: 21 March 2017 Accepted: 26 September 2017 Published online: 11 October 2017
Trade routes and plague transmission in pre-industrial
Europe
Ricci P. H. Yue1, Harry F. Lee1'2 & ConnorY. H. Wu3
Numerous historical works have mentioned that trade routes were to blame for the spread of plague in European history, yet this relationship has never been tested by quantitative evidence. Here, we resolve the hypothetical role of trade routes through statistical analysis on the geo-referenced major trade routes in the early modern period and the 6,656 geo-referenced plague outbreak records in AD1347-1760. Ordinary Least Square (OLS) estimation results showthat majortrade routes played a dominant role in spreading plague in pre-industrial Europe. Furthermore, the negative correlation between plague outbreaks and their distance from major trade ports indicates the absence of a permanent plague focus in the inland areas of Europe. Majortrade routes decided the major plague outbreak hotspots, while navigable rivers determined the geographic pattern of sporadic plague cases. A case study in Germany indicates that plague penetrated further into Europe through the local trade route network. Based on our findings, we propose the mechanism of plague transmission in historical Europe, which is imperative in demonstrating how pandemics were spread in recent human history.
Plague is initiated by the flea-borne bacterium Yersinia pestis, which circulates mainly on rodents and other mammal hosts through the rodents associated fleas1'2. Normally, the bloodsucking fleas acquire Yersinia pestis from an infected rodent. The bacterium will quickly multiply and cluster, leading to the blockage of the alimentary canal in the fleas' guts3. When the infected flea jumps onto another mammal, preferably rodents, it will transmit the bacteria to the new host by regurgitating the clotted blood from the blockage of the alimentary canal4. If an infected flea attempts to feed on a human, it will transmit Yersinia pestis to the human and lead to human plague in the form of bubonic plague or pulmonary plague. Traditional thought suggested that the clustering of Yersinia pestis rarely happened on human fleas5. Yet, recent evidence revealed that not only rodent fleas (Xenopsylla cheopis as a classic example) are to blame for plague transmission, human fleas (Pulex irritans) and cat fleas (Ctenocephalides felis) are also likely to play a role in disseminating Yersinia pestis6'7. Furthermore, laboratory results illustrated the experimental possibility of oral route transmission of plague8 and epidemiological records suggested plague infection through consumption of contaminated meat9'10.
Despite improvements in sanitation and medical advancements in the course of human history, plague remains a major threat to human beings11'12. Over the past few decades, thousands of cases of human plague have occurred around the world, particularly in Africa13. Some researchers suggested that plague may reign over our planet again when global climate change makes some places on earth become wetter and hotter1416. Considering the widespread wildlife reservoirs of plague foci, together with the quick spread, rapid clinical course, inherent communicability, and high mortality rate of plague, the risk of plague outbreak should never be underestimated, although the number of human plague cases is relatively low compared to other infectious diseases at present17. Flashing back in history, plague caused three great pandemics, in which 200 million people perished18. Given the nature of plague and its notorious history, the international community should be more prepared for the re-emergence of plague. Ironically, remarkably little research has been done to elucidate how plague spread through metastatic spatial domains and the mechanism behind its distribution.
The idea that infectious diseases are spread by transportation routes has been supported by proof from various studies on infectious diseases1922. In an earlier study, the authors also hypothesized the possible role of trade routes in plague distribution through examining the linkage between navigable rivers and plague outbreaks23. Using the
department of Geography, The University of Hong Kong, Hong Kong, China, international Center for China Development Studies, The University of Hong Kong, Hong Kong, China, department of Population Health Sciences, Virginia-Maryland College of Veterinary Medicine, VirginiaTech, Virginia, USA. Correspondence and requests for materials should be addressed to R.P.H.Y. (email: ricciyue@connect.hku.hk) or H.F.L. (email: harrylee@hku.hk)
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
1
www.nature.com/scientificreports/
geo-referenced Old World trade route database prepared by Ciolek24 and the historical human plague outbreak distribution in Europe by Biintgen et al.25, this study examined the extent to which major trade routes shaped the dispersal of plague in Europe between AD 1347 and 1760. The study timespan was determined by the availability of data before the Industrial Revolution in Europe (see SI Text for details). During the late medieval period in Europe, long-distance maritime trade routes were developed as a result of the commercial revival of Italy26. Combined with the already active overland trade routes in Western Europe27, they formed a comprehensive network that linked European cities together. Commodity goods, agricultural products, and luxury items were circulated amongst major port cities and the continental hinterland beyond until the Industrial Revolution introduced a new era of trading patterns. It was also a period when the agricultural sectors of Europe had a surplus of labor, and peasants were pushed to the city, generating rapid urban growth28. Human movement became more frequent than in previous eras. Under the influence of these conditions, infectious disease gained access and penetrated every corner of Europe. Would it be possible that trade routes was one of the determining factor in shaping the pattern of plague outbreak in historical Europe? We sought to answer this question in the present study.
Other scholars have long identified the influence of trade routes on plague transmission in historical Europe20'23,29,30. Yet, no scientific consensus could be reached concerning the coherence of plague outbreaks and major trade route patterns. Here, we based on our statistical results to prove that there would be more plague outbreaks when a city is closer to the major trade routes. Moreover, the evidence found in our investigation of the plague/trade-port relationship did not indicate any sign of permanent plague reservoirs in the inland areas of Europe in history. Plague was imported from trade ports or it originated from somewhere linked to the maritime trading system. Sporadic plague outbreak was used as an indicator to show the different roles of trade routes and navigable rivers in plague transmission. A specific case study of Germany indicated that a local trade route was more important in distributing plague outbreak. Altogether, the results allowed us to propose a potential plague transmission mechanism.
Results
We combined records of 6,656 plague outbreak cases in historical Europe and North Africa and the trade route database that geo-referenced the major overland and maritime trade routes during the early modern period (Fig. 1). To examine whether trade routes were related to the plague outbreak patterns during our study period, we started by checking whether plague hotspots were also key trade nodes. According to a recent study by Schmid et al.30, there was never any permanent plague reservoir in pre-industrial Europe. This implied that a plague outbreak at any given place in our study area was transmitted from a nearby outbreak. If human movements and the circulation of goods provided an ideal channel for the spread of plague, a city's proximity to key trade nodes would determine its likelihood to become a plague hotspot.
The top 20 cities with the highest year count of plague outbreaks between AD 1347 and 1760 are listed in Table 1. Thirteen of these cities were key trade nodes that linked trade routes together. These cities were geographically dispersed and spread in seven countries. As stated by Vogler et al.31, maritime ports were common plague outbreak centers when "plague ships" introduced infected rodent hosts and flea vectors to the cities. In addition, even though half of the trade nodes as documented in our dataset are port cities, only six port cities (in six countries) are listed in Table 1. The same result was also obtained with a sensitivity test (Table S2). Briefly, plague hotspots were mostly trade nodes. Yet, there seemed to be no evidence to support that those hotspots were necessarily port cities. We dissimilated this pattern by searching for the relationship between plague outbreak and its distance to trade route and trade port.
If plague was spread along trade routes, the closer one got to a trade route, the easier it would be to become infected. Therefore, we calculated the distance between each plague outbreak and its closest trade route and estimated its relationship with the accumulated count of each plague outbreak point. Several control variables were included in our Ordinary Least Square (OLS) regressions to investigate the robustness of the relationship. Our results are summarized in Table 2. Model 1 was the base model for estimating the relationship between logged distance to trade route and plague count. The association was negative, which supported our hypothesis. In Model 2, the regression was run with both the time fixed-effects and the region fixed-effects to control the influence of any observable or unobservable predictors over the dependent variable. Hence, bias on omitted variables was greatly reduced. The relationship between logged distance to trade route and plague count remained negatively significant. In Model 3, the regression was run with additional geographical controls. As climate is one of the most dominant factors in determining the prevalence of plague13,15,16,32,33, we controlled the effect of temperature on plague reoccurrence by incorporating the elevation and the latitude of plague location into the regression. In addition, the longitude variable also helped to capture the difference in plague distribution between the eastern part and the western part of Europe34. The association continued to remain negatively significant. Further control on climatic zone did not affect the negative correlation between plague count and its distance from trade route (Table SI). As there might be difference between North Africa and Europe, as well as between coastal city and inland city, in terms of their plague recurrence, we included a North Africa dummy and a coastal city dummy in Model 4. Our results showed that no significant difference is seen between North Africa and Europe. In parallel, although coastal cities were more prone to plague outbreak, it did not distort the negative correlation between plague count and its distance to trade route. In Model 5, we included vegetation cover and normalized population density, which symbolizes the degree of urban development over different plague outbreak points, as control variables in the regression. In Model 6, we included per capita Gross Domestic Product (GDP), Consumer Price Index (CPI), and normal laborer wages as control variables in the regression. The specifications in Models 5 and 6 were set according to the practices of traditional public health studies. The rapid and unplanned urbanization process would promote the spread of epidemics35, while landscape contexts were associated with plague occurrence36. On the other hand, accumulated wealth and improvements in the living standard in early modern Europe might dampen plague reoccurrence. Nevertheless, the inclusion of these control variables did not alter the significance and robustness of our estimated relationship.
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
2
www.nature.com/scientificreports/
Legend
Port city with plague outbreaks Frequency of plague outbreaks
0 40+
0 31-40
• 21-30
• 5-20
• 1-4
• Port city without plague outbreaks —— Trade routes Plague outbreaks FREQUENCY
Q 40+ 0 31-40
• 21-30 © 5-20
• 1 -4
Figure 1. Spatial distribution of plague outbreak in Europe, and Northern Africa, AD1347-1760. Plague outbreaks are related to the patterns of trade routes, both overland and maritime, and also major trade ports in pre-industrial Europe. Cities with recorded plague outbreaks are marked with red dots, with the size of dots referring to the number of plague outbreak during the study period (See legends). The blue lines indicate the major trade route in early modern Europe. The black dots identify the locations of major trade ports with plague outbreak over the study period. Major trade ports without plague outbreak over the study period are labeled in grey dots. Trade routes and trade ports at countries with no plague record are omitted. From our results, more plague outbreaks happened in the periphery of trade routes and trade ports. The map is generated in ArcGIS version 10.1 (www.esri.com/software/arcgis).
The OLS estimates revealed that distance to trade route has a high explanatory power to the distribution of plague outbreak. Despite its large sample size (n = 6,656), the R2 ranged from 0.41 to 0.44 in Models 2 to 6. The negativity of the association implies that cities closer to trade routes were more vulnerable to plague reoccurrence. On the other hand, being further away from trade routes was a good way to escape from plague.
We performed several sensitivity tests to further check the robustness of our OLS regression results (Table S3). It was shown that the relationship between plague outbreak and trade route was highly significant and remained negative in different temporal domains of our study period (AD 1347-1449, AD1450-1549, AD1550-1649, and AD1650-1760). In addition, the relationship was robust in our different specification of spatial domains, in which the cases in Russia and Africa were excluded, or only the cases in continental Europe and the six major plague outbreak countries were included (see SI for more details). The above results implied that the pattern of plague outbreak was determined by the trade route patterns for the entire study area, which was consistent over the study period. Hence, the relationship should be independent of cultural, demographic, economic factors or the possible spatial bias of plague database as suggested by Alfani37. Otherwise, the relationship over different regions and time periods should differ.
The pattern of plague outbreak in historical Europe was related to the distribution of major trade routes at that time. However, was plague circulating within villages and cities and being retained as a reservoir, as suggested by Ell38? Or was plague repeatedly introduced from major trade ports to inland areas, as suggested by Schmid et al?0? Here we used the distance to major trade ports as an indicator to determine whether there was a permanent plague focus in historical Europe (Fig. 1). A few plausible scenarios should be considered. If there was a permanent plague focus in historical Europe, plague would circulate between cities and villages, and not necessarily be transmitted through trade ports. Therefore, we should not detect any significant relationship between plague outbreak and distance to trade port. However, there might also be a scenario in which the permanent plague focus was so strong that it kept on exporting plague from inland to other parts of Europe. In this case, the relationship being examined in this section might become positively significant. At the same time, it might be possible that trade ports were always the entrance point or starting point for plague transmission, which moved inland until that wave of plague outbreak was no longer able to sustain itself and faded. In such a case, we should detect a negatively significant correlation between distance to trade port and plague outbreak. The OLS estimation results for
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
3
www.nature.com/scientificreports/
Name of City	Count of Year with Plague Outbreak	Key Trade Node	Major Port City
London	132	Yes	Yes
Alger	90	Yes	Yes
Paris	89	Yes	No
Toulouse	88	Yes	No
Bourg-en-Bresse	77	No	No
Amiens	75	No	No
Nantes	73	No	No
Bordeaux	68	Yes	Yes
Venezia	63	Yes	Yes
Strasbourg	59	Yes	No
Rouen	57	No	No
Limoges	57	Yes	No
Basel	56	Yes	No
Angers	55	No	No
Dijon	53	Yes	No
Troyes	52	No	No
Bremen	46	No	No
Barcelona	46	Yes	Yes
Gdansk	44	Yes	Yes
Milan	44	Yes	No
Table 1. List of top 20 plague hotspot cities and their role in old world trade route in Europe, AD 1347-1760. Key trade node refers to the major transportation node connecting medieval Europe as indicated by Evans and Brooke41 and Spufford42. Major port city indicates whether it also functioned in maritime trade route.
the above scenarios were reported in Table S4. It was clear that distance to major trade ports was negatively correlated to plague outbreak (p = 0.000, F = 114.38). The negative relationship highlights that more plague outbreaks were recorded when distance to trade port decreased, suggesting that there should not be any permanent plague reservoir in the inland part of Europe. We also conducted a detailed robustness check on how the pattern of plague outbreak was controlled by the distance to trade ports over different specifications in spatial and temporal domains (Table S5). Overall, the results showed a constantly highly significant and negative relationship between distance to trade port and plague. The relationship was robust and did not vary in any timespan or region that we singled out from the database. This further confirmed the estimation that plague was repeatedly introduced to the inland of Europe through maritime trade ports in our study period. Plague was introduced to Europe mainly by maritime trade routes and entered overland trade routes until that wave of plague outbreak faded away or strengthened again by a new introduction of plague supply. Or it was possible that any of the major ports could be a strong permanent plague reservoir. As Davis39 has suggested, in the "upper half" of Europe, the climate was a disadvantage to black rat reproduction. The population of black and brown rats could never be sustained as an epidemic focus without a continuous supply of new rats by ships40.
In the previous analysis, we verified that plagues were repeatedly introduced to Europe through major trade ports, following the trade route to various trade nodes as well as to the cities around the trade routes. This, however, gave rise to another question regarding the extent of the impact of trade route to plague outbreak. Is trade route responsible for shaping all the plague outbreak patterns in historical Europe? To answer this question, we narrowed down our investigation to the sporadic plague outbreak points that had four or fewer outbreaks during our study period (Fig. 2). As shown in Fig. 2, sporadic plague outbreaks (frequency < 4) are scattered around on the map. These sporadic cases did not follow the pattern as described in our previous analysis. The associated statistical result also indicated the relationship between trade routes and plague outbreak to be insignificant (Table S6).
However, there must be certain ways for animal hosts to cause these sporadic outbreaks scattered around the European continent. It might be attributable to some less active transportation routes (e.g., navigable rivers) that connected these sporadic cases with the plague hotspots. We calculated the correlation between distance to the closest navigable river and plague outbreak for these sporadic cases (Table S7). We found that although navigable rivers were not as capable as major trade routes in influencing the total plague transmission (Table S8), they did account for the pattern of sporadic plague distribution in historical Europe (p < 0.005; F = 6.13). There were more outbreaks of plague in those cities located closer to the navigable rivers. The result was also verified by robustness checks in various geographical specifications. It confirmed that local river channels, instead of major trade routes, were more significant in determining the distribution of sporadic plague outbreak cases in early modern Europe.
The major trade route database by Evans and Brooke41 and Spufford42 might have neglected the trade route-plague relationship at the local level, as they primarily address the major trade routes and trade ports at the continental scale. Here, we employed a very fine-gained trade route database in Germany, which consists of both major and local trade routes (Fig. 3), to see whether our findings about the trade route-plague relationship could be supported. Our results showed that the parts of Germany that were closer to trade routes, measured
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
4
www.nature.com/scientificreports/
	Model 1	Model 2	Model 3	Model 4	Model 5	Model 6
Log(Distance from city centre)	-6.086034*** (3 = -0.4665029)	-6.181281*** (3 = -0.4738037)	-6.175487*** (3=-0.4733596)	-6.289981*** (3=-0.4821357)	-6.175596*** (3=-0.473368)	-6.174889*** (3=-0.4733137)
Latitude			-0.1978809 (3=-0.0892771)	-0.0283914 (3=-0.070624)	-0.2012694 (3=-0.0500661)	-0.203408 (3 = 0.0505981)
Longitude			1.349016*** 03 = 0.2379364)	1.212566*** (3=0.2138698)	1.349467*** (3 = 0.2380161)	1.368032*** (3 = 0.2412906)
Elevation			-0.0175064*** (3=-0.1049527)	-0.023406*** (3=-0.1403218)	-0.0174867*** (3=-0.1048349)	-0.0174736*** (3=-0.104756)
Coast indicator				-9.778023*** (3=-0.1281232)		
North Africa Indicator				6.565075 (3=0.0484845)		
Vegetation Cover					0.2457425 (3 = 0.0913583)	0.3816018 (3 = 0.1418658)
Normalized Population Density					0023217 (3 = 0.0699493)	0.0025352 (3 = 0.0763817)
Per Capita GDP						-0.0101291 (3=-0.0857906)
CPI						2.320846 (3 = 0.0386653)
Normal Wage						0.7168337 (3=-0.346733)
Time Fixed effect	No	Yes	Yes	Yes	Yes	Yes
Regional Fixed effect	No	Yes	Yes	Yes	Yes	Yes
Number obs.	6656	6656	6656	6656	6656	6656
R2	0.2176	0.4144	0.4339	0.4434	0.4340	0.4345
Table 2. OLS estimates of relationship between plague outbreak and trade route in Europe, AD 1347-1760. Notes. The dependent variable of checks is the total number of plague reoccurrence. ***p < 0.005; **p < 0.01; *p<0.05.
Figure 2. Sporadic plague outbreak (N < 5) did not follow the pattern of major trade routes. The red spots represent the locations of sporadic outbreak of plague (856 locations). The blue lines indicate the major trade route within our study period. The black dots identify the major trade ports with plague outbreak. The trade ports with no reported plague outbreak within our study period are labeled in grey dots. The map is generated in ArcGIS version 10.1 (www.esri.com/software/arcgis).
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
5
www.nature.com/scientificreports/
Legend
Figure 3. Distribution and frequency of plague outbreak in relation to the local Holy Roman Empire trade route in Germany, AD1347-1760. It can be seen that locations with more plague recurrence (as referred by the size of red dots) are closer to the local Holy Roman Empire trade route (blue lines). The strength of recurrence fades in according to the distance away from these trade routes as suggested by the statistical analysis (Table 3). The map is generated in ArcGIS version 10.1 (www.esri.com/software/arcgis).
Model	Log (distance)	Time fixed effect	Number of obs.	F	R2
Chief trade route database by Evans and Brooke41 and Spufford42					
1	-1.843461*** (-0.2278522)	Yes	728	2.19	0.4595
Local trade route dataset of Germany in Holy Roman Empire time by Davies et al47					
2	-4.930627*** (-0.1991022)	Yes	728	6.63	0.7199
Navigable	rivers				
3	-3.534741*** (0.5771989)	Yes	728	1.97	0.4326
Table 3. OLS estimates of relationship between plague outbreak and chief/local trade route and navigable river in Germany, AD 1347-1760. Notes. The dependent variable of model is the total number of plague reoccurrence. ***p < 0.005; **p < 0.01; *p < 0.05.
by logged distance, tended to be the areas suffering repeated plague occurrence (Model 2 in Table 3). Also, we repeated our analysis by using the major trade route data (in Germany only) complied by Evans and Brooke41 and Spufford42 (Model 1 in Table 3) for cross-validation (Model 2 in Table 3). The same conclusion was reached. In addition, we tested for the role of navigable rivers in plague transmission in Germany (Model 3 in Table 3), and the navigable rivers were shown to be less influential than trade routes (Models 1 and 2) in plague transmission, which is consistent with our results for the whole of Europe (Table 2 and Table S8). The above findings indicated that local trade routes, which are linked with and nested within the major trade routes networks, also contributed to plague distribution significantly. Yet, the absence of fine-grained local trade route databases constrained our investigation on the topic in Germany only.
Discussions
Prior to the Industrial Revolution, long distance human movement was mainly confined to the trading of goods along certain overland trade routes, navigable rivers, and maritime trade routes43. Plague was spread primarily by its rodent host or occasionally vectorized by other animal hosts. However, these hosts did not move across Europe themselves, but were transported by humans along the major trade routes. By combing data from historical trade routes and plague records, we found that the geographic pattern of plague was determined by major trade routes in early modern Europe.
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
6
www.nature.com/scientificreports/
Figure 4. Possible plague spreading pattern from port to inland in Europe, AD1347-1760. Plagues were carried from other permanent plague focus to major trade ports in Europe. The contagion will go further to the hinterland by major trade routes or navigable river connecting the major trade ports. By transporting through the major trade route, the contagion will eventually focus at major trade node, resulting in the formation of plague hotspot in historical Europe. The pathway from major trade route to plague hotspot would also pass through local trade route. Certain amount of contagion would enter nearby navigable rivers from major trade route or major trade node. The navigable rivers would further carry the contagion inland and create sporadic cases all over the European continent.
There were five key findings in our statistical results. First, places closer to trade routes were more prone to plague outbreak and thus, plague reoccurrence. Second, plague was repeatedly introduced to several key trade ports and spread further inland in early modern Europe. Third, we found no sign of a permanent plague reservoir according to the distribution of plague outbreak. Forth, localized river navigation systems instead of trade routes accounted for the geographic distribution of sporadic plague outbreaks. Fifth, the case study of Germany suggested that local trade routes could be a significant explanatory factor to plague transmission. Based on these findings, we proposed a hypothetical corollary for plague transmission in historical Europe (Fig. 4).
Concurring with the work of Schmid et al.30, we found no evidence for a permanent plague focus in the inland of Europe. Thus, our corollary started with the assumption that plague must have been imported from a permanent plague focus to/at major trade ports. The disease probably rode on its animal hosts and approached key maritime trade ports or inland trade ports by ship. These urbanized ports and commercial hubs provided perfect conditions, such as grain warehouses, high humidity, and dense human settlements, for the establishment of epidemic epi-centers. The contagion would then pass on to major trade nodes through major trade routes. Those major trade nodes which linked up multiple trade routes would have a higher probability to become plague hotspots, as they were often connected with infected ports, or they were the infected ports themselves. Also, contagion would be likely to stop at these trade nodes, transforming them as additional epi-centers. As described by Benedictow44 in his work about plague transmission on the Silk Road, infected caravans would stop in caravanserai, transmitting plague to other caravans and distributing the disease across trade routes. The same phenomenon would probably apply to historical Europe. These commercial hubs would be plague hotspots and the disease was further disseminated through major trade routes such as navigable rivers, overland routes, or routes across coastal ports. A piece of supporting evidence to this plausible plague transmission mechanism was that in historical Europe, cities were small and there were no suburbs linking cities together43. Probably, only major trade routes would provide enough contagion density to sustain plague transmission. For other routes that were poorly built for transportation, infected people would barely survive the trip29. Plague would follow these major trade routes and spread to the periphery. As such, cities closer to these trade routes would have a bigger chance of infection.
It might be argued that maritime trade routes were not the sole way for plague to conquer the European continent. Plague might be repeatedly introduced to Europe through the overland trade routes such as the Silk Road. In such case, a plague pandemic would transfer from inland to port and then it was exported to other places of Europe. However, this plague transmission pathway could not be supported by our statistical results. Our robustness check showed that plague was either imported to Europe by maritime trade routes or a specific trade port was developed as a plague reservoir.
According to our results, major trade routes did not account for sporadic plague cases. In fact, the imperative role of navigable rivers in connecting cities after the medieval era has been mentioned by Edwards and Hindle45 and Jones46 and is also highlighted by our statistical analyses in spreading plague. Inland waterways provided further penetration for the contagion from various major trade routes to the hinterland. Our results proved that the further away from navigable rivers, the number of plague outbreak dropped in a statistically significant manner.
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
7
www.nature.com/scientificreports/
Certain contagions might have left the major trade routes and spread to other settlements through inland waterways. Although these sporadic cases seemed to be randomly distributed, they were indeed anchored with trade routes.
The above explanation might not totally account for the pattern of sporadic plague outbreak in historical Europe as indicated by the low F and R2 values in the statistical analysis (F = 6.13, R2 = 0.0748, Table S7). Certain sporadic cases in our database would be attributable to other factors such as war or undocumented localized trade/communication routes. Also, sporadic events which cannot be addressed by our research methods might actually cluster in certain temporal and spatial scales. However, when we looked at long distance travel or major carriers of plague in a long-temporal and large-spatial perspective, major trade route largely accounts for the distribution of plague outbreaks.
Ideally, we would prefer to test our model with local trade route at high spatial resolution in Europe against plague distribution. Unfortunately, these data are unavailable at the moment. Given this shortcoming, we were only able to issue a case study of Germany as evidence of the role of local trade routes in plague transmission. The result further validated the possible linkage between trade route and plague recurrence, indicating that trade routes in higher resolution and local context were connected closer to plague outbreak. One possible explanation for this was that plague entered these local trade routes through the major trade route network, although it remained unclear whether this explanation could be applied in other countries.
To sum up, this study illustrated a plausible pathway for plague transmission in Europe in AD1347-1760. It had important implications in explaining how the plague outbreak pattern was shaped, and how the plague hot-spots were generated, by major trade routes. The sign of permanent plague reservoirs in the inland of historical Europe could not be substantiated with strong evidence. Yet, it remains possible, according to the statistical result, that a plague reservoir once existed at the trade ports. The correlation between navigable rivers and sporadic plague outbreak cases might supplement the current explanation of plague distribution, which might provide new insight in examining the geographic patterns of plague. We did not exclude other plausible explanations for plague distribution, as there might be other factors favoring or hindering plague outbreaks in different temporal and spatial domains. Future works should focus on whether our proposed mechanism behaves differently in other spatio-temporal settings. Also, assessing whether a clustering effect in time or space exists might help explain sporadic cases, which is useful in forecasting future plague outbreaks. Most importantly, the case study result of Germany augments the high resolution trade route dataset in explaining more of the relationship between trade route and infectious disease in history. Our findings may contribute to plague prevention and mitigation, especially in the Third World where living conditions and transportation means are similar to those in pre-industrial Europe. Both major transportation routes and navigable waterways are high risk areas to be closely monitored.
Methods
Data. Our geo-referenced trade routes data were retrieved from the Old World Trade Routes Project built by Ciolek24. We combined the datasets of Evans and Brooke41 and Spufford42 to reconstruct the major trade route network from late medieval Europe to early modern time. The local trade route in Germany is originated from the study of Davies et al.47. Our geo-referenced plague outbreak data came from the database of Biintgen et al.25. Their work was originated from the literature review done by Biraben48. The geo-referenced database also provided coordinates for latitude and longitude of each plague outbreak case. Elevation data came from measurement through ArcGIS. Vegetation cover data of historical Europe come from the study of Kaplan etal.49. Population density was calculated from dividing the historical demographic figure from McEvedy and Jones50 by current regional area51. It was further normalized according to the data of historical urban area49. We retrieved the per capita Gross Domestic Product (GDP) data for our study area from the database of Bolt and Zanden52 and Maddison53. Consumer Price Index (CPI) and normal laborer's wages were obtained from the study of Allen54. North Africa indicator included plague outbreaks that happened in North Africa and Turkey. Coastal indicators referred to plague outbreaks that happened at a point less than 5 km to the current coastline. Data on distance of river were acquired by measurement in ArcGIS. The definition of navigable river was based on the study of McGrail55 and Eckoldt56 and therefore we only included rivers wider than 5 m and those rivers that have a connection with other cities. We only measure rivers within a 10 km radius of plague outbreak points. For those cases where rivers were too small or too far away from plague outbreak points, 10,000 m was manually set as the distance. Please refer to SI Text for the details of our datasets.
Ordinary Least Square (OLS) estimation. We hypothesize the pattern of plague outbreak to be determined by trade routes. To validate the robustness of the above relationship, we include various control variables in our regression models. Our base regression model is:
Pc = cc + P(logDi) + X'fi + C'i4> + £,-, (1)
where Pc is the total number of plague count at an individual point, cc is the intercept of the equation, f] is the coefficient of the association, logDi stands for the logged distance to the closest major trade route from each individual plague outbreak point, Xt is the vector of controlled variables in the model, C, is the country fixed-effects estimator that controls the differences among regions, ei is the error term. Time fixed-effects estimator is only added to those models with time variants. The relationship between trade route and plague outbreak remains significant even if the time-fixed effects estimator is excluded. The above method is also applied in the remaining parts of our statistical analyses.
Robustness checks. We perform various robustness checks to validate the significance of our results in different temporal and spatial domains. We built different regression models to explore the impact of trade routes on
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
8
www.nature.com/scientificreports/
plague outbreak patterns in different regional settings. This helps to avoid the potential regional discrepancies that are not controlled in our models. Furthermore, factors such as technological improvement are hard to quantify, but they are related to the performance of trade routes over time. Therefore, we slice our data into different time sections to see whether trade routes were still significant in shaping plague outbreak pattern over time. The results of these robustness checks are presented in the Appendix.
References
1. Perry, R. D. & Fetherston, J. D. Yersinia pestis-etiologic agent of plague. Clinical microbiology reviews 10,35-66 (1997).
2. Gage, K. L. & Kosoy, M. Y. Natural history of plague: perspectives from more than a century of research. Annu Rev Entomol 50, 505-528, https://doi.org/10.1146/annurev.ento.50.071803.130337 (2005).
3. Burroughs, A. L. Sylvatic plague studies. The vector efficiency of nine species of fleas compared with Xenopsylla cheopis. Journal of Hygiene 45,371-396 (1947).
4. Hinnebusch, B. J., Perry, R. D. & Schwan, T. G. Role of the Yersinia pestis hemin storage (hms) locus in the transmission of plague by fleas. Science 273, 367 (1996).
5. Hinnebusch, B. J. The evolution of flea-borne transmission in Yersinia pestis. Current issues in molecular biology 7,197-212 (2005).
6. Laudisoit, A. et al. Plague and the human flea, Tanzania. Emerging infectious diseases 13,687 (2007).
7. Eisen, R. J. et al. Early-phase transmission of Yersinia pestis by cat fleas (Ctenocephalides felis) and their potential role as vectors in a plague-endemic region of Uganda. The American journal of tropical medicine and hygiene 78, 949-956 (2008).
8. Butler, T, Fu, Y, Furman, L., Almeida, C. & Almeida, A. Experimental Yersinia pestis infection in rodents after intragastric inoculation and ingestion of bacteria. Infection and immunity 36,1160-1167 (1982).
9. Saeed, A. A. B., Al-Hamdan, N. A. & Fontaine, R. E. Plague from eating raw camel liver. Emerging infectious diseases 11,1456 (2005).
10. Malek, M. A., Bitam, I. & Drancourt, M. Plague in Arab Maghreb, 1940-2015: A Review. Frontiers in public health 4 (2016).
11. Zietz, B. P. & Dunkelberg, H. The history of the plague and the research on the causative agent Yersinia pestis. International journal of hygiene and environmental health 207,165-178 (2004).
12. WHO. Interregional Meeting on Prevention and Control of Plague: Antananarivo, Madagascar, 1-11 April 2006. (World Health Organization, 2008).
13. Kreppel, K. S. et al. A non-stationary relationship between global climate phenomena and human plague incidence in Madagascar. PLoS Negl Trop Dis 8, e3155 (2014).
14. Enscore, R. E. et al. Modeling relationships between climate and the frequency of human plague cases in the south-western United States, 1960-1997'. American Journal of Tropical Medicine and Hygiene 66,186-196 (2002).
15. Kausrud, K. L. et al. Climatically driven synchrony of gerbil populations allows large-scale plague outbreaks. Proc Biol Sci 274, 1963-1969, https://doi.org/10.1098/rspb.2007.0568 (2007).
16. Stenseth, N. C. et al. Plague dynamics are driven by climate variation. Proc Natl Acad Sci USA 103, 13110-13115, https://doi. org/10.1073/pnas.0602447103 (2006).
17. Stenseth, N. C. et al. Plague: past, present, and future. PLoS Med 5, e3 (2008).
18. Sherman, I. In the Power o/Plagues (ed I. Sherman) Ch. 4, 66-87 (ASM Press, 2006).
19. Messina, J. P. etal. Global spread of dengue virus types: mapping the 70 year history. Trends in microbiology 22,138-146 (2014).
20. Xu, L. et al. Wet climate and transportation routes accelerate spread of human plague. Proceedings of the Royal Society B: Biological Sciences 281,20133159 (2014).
21. Tatem, A. J., Rogers, D. J. & Hay, S. I. Estimating the malaria risk of African mosquito movement by air travel. Malaria journal 5, 57 (2006).
22. Morse, S. S. Factors in the emergence of infectious diseases. (Springer, 2001).
23. Yue, R. P. H., Lee, H. F. & Wu, C. Y. H. Navigable rivers facilitated the spread and recurrence of plague in pre-industrialEurope. Scientific Reports, doi:srep34867 (2016).
24. Ciolek, T. M. In Old World Trade Routes(OWTRAD) Projects (ed OWTRAD Dromographic Digital Data Archives(ODDDA)) (Canberra: http://www.ciolek.com, 1999).
25. Buntgen, U.,Ginzler, C, Esper, J., Tegel, W. & McMichael, A. J. Digitizing historical plague. Clinical infectious diseases 55,1586-1588 (2012).
26. Tracy, J. D. The Rise of Merchant Empires: LongDistance Trade in the Early Modern World 1350-1750. Vol. 1 (Cambridge University Press, 1993).
27. Munro, J. H. The Low Countries' export trade in textiles with the Mediterranean basin, 1200-1600: a cost-benefit analysis of comparative advantages in overland and maritime trade routes. (1999).
28. Hilton, R. Bond men made free: Medieval peasant movements and the English rising of 1381. (Routledge, 2003).
29. Benedictow, O. J. The Black Death, 1346-1353: the complete history. (Boydell & Brewer, 2004).
30. Schmid, B. V. et al. Climate-driven introduction of the Black Death and successive plague re introductions intoEurope. Proceedings of the National Academy of Sciences 111, 3020-3025 (2015).
31. Vogler, A. J. etal. A decade of plague in Mahajanga, Madagascar: insights into the global maritime spread of pandemic plague. MBio 4, e00623-00612 (2013).
32. Zhang, Z. et al. Relationship between increase rate of human plague in China and global climate index as revealed by cross-spectral and cross-wavelet analyses. Integr Zool 2,144-153, https://doi.Org/10.llll/j.1749-4877.2007.00061.x (2007).
33. Snail, T, O'Hara, R. B., Ray, C. & Collinge, S. K. Climate-driven spatial dynamics of plague among prairie dog colonies. Am Nat 171, 238-248, https://doi.org/10.1086/525051 (2008).
34. Nunn, N. The Long-Term Effects of Africa's Slave Trades. The Quarterly Journal of Economics, 139-176 (2008).
35. Moore, M., Gould, P. & Keary, B. S. Global urbanization and impact on health. International journal of hygiene and environmental health 206,269-278 (2003).
36. Collinge, S. K. et al. Landscape structure and plague occurrence in black-tailed prairie dogs on grasslands of the western USA. Landscape ecology 20, 941-955 (2005).
37. Alfani, G. Plague in seventeenth-century Europe and the decline of Italy: an epidemiological hypothesis European Review of Economic History, 408-430 (2013).
38. Ell, S. R. Immunity as a factor in the epidemiology of medieval plague. Review of Infectious Diseases 6,866-879 (1984).
39. Davis, D. E. The scarcity of rats and the Black Death: an ecological history. The Journal of Interdisciplinary History 16, 455-470 (1986).
40. Hufthammer, A. K. & Walloe, L. Rats cannot have been intermediate hosts for Yersinia pestis during medieval plague epidemics in Northern Europe. Journal of Archaeological Science 40,1752-1759 (2013).
41. Evans, J. & Brooke, C. N. L. The flowering of the Middle Ages. 360 (Thames & Hudson, 1967).
42. Spufford, P. Power and profit: the merchant in medieval Europe. (Thames & Hudson London:, 2002).
43. Rodrigue, J.-P, Comtois, C. &Slack,B. The geography of transport systems. (Routledge, 2013).
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
9
www.nature.com/scientificreports/
44. Benedictow, O. J. Yersinia pestis, the Bacterium of Plague, Arose in East Asia. Did it Spread Westwards via the Silk Roads, the Chinese Maritime Expeditions of Zheng He or over the Vast Eurasian Populations of Sylvatic (Wild) Rodents? Journal of Asian History 47,1-31 (2013).
45. Edwards, J. F. & Hindle, B. P. The transportation system of medieval England and Wales. Journal of Historical Geography 17,123-134 (1991).
46. Jones, E. T. River navigation in medieval England. Journal of Historical Geography 26,60-75 (2000).
47. Davies, N & Moorhouse, R. Microcosm: A portrait of a Central European City. (Pimlico Random House, 2003).
48. Biraben, J.-N. Les hommes etlapeste en France et dans lespays europeens et mediterraneens. Vol. 2 (Mouton Paris, 1975).
49. Kaplan, J. O., Krumhardt, K. M. & Zimmermann, N. The prehistoric and preindustrial deforestation of Europe. Quaternary Science Reviews 28,3016-3034 (2009).
50. McEvedy, C. & Jones, R. Atlas of world population history. (Penguin Books Ltd, Harmondsworth, Middlesex, England., 1978).
51. DK-Publishing, I. DK world atlas. 338 (DK Publishing, 1997).
52. Bolt, J. & Zanden, J. L. The Maddison Project: collaborative research on historical national accounts. The Economic History Review 67,627-651 (2014).
53. Maddison, A. The world economy volume 1: A millennial perspective volume 2: Historical statistics. (Academic Foundation, 2007).
54. Allen, R. Allen-Unger Database: European Commodity Prices 1260-1914. On-line: http://www2.history.ubc.ca/unger/htmfiles/new grain.htm (2007).
55. McGrail, S. Early Ships and Seafaring. European Water Transport. Barnsley: Pen & Sword Archaeology (2014).
56. Eckoldt, M. Navigation on small rivers in Central Europe in Roman and Medieval times. Jnternational Journal of Nautical Archaeology 13, 3-10 (1984).
Acknowledgements
This research was supported by the Hui Oi-Chow Trust Fund (201502172003 and 201602172006) and Research Grants Council of The Government of the Hong Kong Special Administrative Region of the People's Republic of China (HKU745113H and 17610715).
Author Contributions
R.P.H.Y. and H.F.L. designed research; R.P.H.Y. and H.F.L. performed research; R.P.H.Y, H.F.L., and C.Y.H.W. analyzed data; and R.P.H.Y. and H.F.L. wrote the paper.
Additional Information
Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-13481-2. Competing Interests: The authors declare that they have no competing interests.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
l/JJv 0 I Open Access This article is licensed under a Creative Commons Attribution 4.0 International ^■s^M License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.Org/licenses/by/4.0/.
© The Author(s) 2017
SCIENTIFIC REPORTS | 7: 12973 | DOI:10.1038/s41598-017-13481-2
10