- Open Access
Hantavirus reservoir Oligoryzomys longicaudatus spatial distribution sensitivity to climate change scenarios in Argentine Patagonia
International Journal of Health Geographicsvolume 8, Article number: 44 (2009)
Oligoryzomys longicaudatus (colilargo) is the rodent responsible for hantavirus pulmonary syndrome (HPS) in Argentine Patagonia. In past decades (1967–1998), trends of precipitation reduction and surface air temperature increase have been observed in western Patagonia. We explore how the potential distribution of the hantavirus reservoir would change under different climate change scenarios based on the observed trends.
Four scenarios of potential climate change were constructed using temperature and precipitation changes observed in Argentine Patagonia between 1967 and 1998: Scenario 1 assumed no change in precipitation but a temperature trend as observed; scenario 2 assumed no changes in temperature but a precipitation trend as observed; Scenario 3 included changes in both temperature and precipitation trends as observed; Scenario 4 assumed changes in both temperature and precipitation trends as observed but doubled. We used a validated spatial distribution model of O. longicaudatus as a function of temperature and precipitation. From the model probability of the rodent presence was calculated for each scenario.
If changes in precipitation follow previous trends, the probability of the colilargo presence would fall in the HPS transmission zone of northern Patagonia. If temperature and precipitation trends remain at current levels for 60 years or double in the future 30 years, the probability of the rodent presence and the associated total area of potential distribution would diminish throughout Patagonia; the areas of potential distribution for colilargos would shift eastwards. These results suggest that future changes in Patagonia climate may lower transmission risk through a reduction in the potential distribution of the rodent reservoir.
According to our model the rates of temperature and precipitation changes observed between 1967 and 1998 may produce significant changes in the rodent distribution in an equivalent period of time only in certain areas. Given that changes maintain for 60 years or double in 30 years, the hantavirus reservoir Oligoryzomys longicaudatus may contract its distribution in Argentine Patagonia extensively.
The colilargo (Oligoryzomys longicaudatus Bennet 1832) belongs to a genus of small-sized mice classified in the New World Tribe Oryzomyini (Muridae: Sigmodontinae). It is a widespread rodent primarily found in woods and shrub lands in Chile and southwestern Argentina . This is a sigmodontine rodent of great importance because of its role as a major reservoir for the Andes Sout lineage of hantavirus, which produces hantavirus pulmonary syndrome (HPS) in southern South America [2, 3].
HPS is a severe and frequently fatal cardiopulmonary disease  caused by a virus of the genus Hantavirus, family Bunyaviridae. The first cases of HPS were reported in the United States in 1993. The etiologic agent, Sin Nombre virus (SNV), was associated with the sigmodontine rodent reservoir, Peromyscus maniculatus. Presently, more than 15 variants of hantavirus have been identified in the occidental hemisphere and many countries (Argentina, Bolivia, Chile, Brazil, United States of America, Panama, Paraguay and Uruguay) were affected by HPS with a lethality ranging from 1.8 to 47.2% . In all cases, HPS cases were associated with sigmodontine reservoirs.
Epidemiological analysis and planning preventive measures for zoonotic and vector-borne diseases require a knowledge of the geographic distribution and ecological conditions relevant to the circulation of a pathogen . Study of these factors have been conducted for vector, reservoir and/or disease transmission which helps in the prediction of potential effects of climate change on diseases like malaria, dengue, schistosomiasis, leishmaniasis and chagas [7–10] and their vectors .
Except for a few works in Chile [12, 13], the knowledge of colilargo ecology is anecdotal or brief at best. The Patagonian region (Figure 1), has suffered in the last century many environmental challenges that affect the rodent's distribution. Sheep ranching doubled between 1900' and 1950' but has since declined . Fire, considered a natural disturbance that shapes the south Andean forests, has also been influenced by humans since the late eighteen century . Grazing and fire together with temperature and rainfall have a complex relation with multiple feedbacks that influence the establishment and growth of woods, shrubs and grasses [16, 17], which in turn affect the colilargo populations. Eruptions in rodent populations (including colilargos) have been associated with bamboo flowering and subsequent seed production as well as with increased rainfall after El Niño . Also population dynamics were associated with the Antarctic Oscillation Index and Southern Oscillation Index in Chile; El Niño events (negative SOI values) had a lagged positive effect on population growth . In the United States not only the abundance of the reservoir but also increased risk of SNV transmission were related to rise in precipitation associated with El Niño . Excepting this last study, most have focused mainly on temporal relationships between climate indexes and population time series. Changes in the spatial distribution, particularly over southern South America, have yet to be explored extensively.
Recently, the spatial distribution of colilargos has been modelled in Argentina as a function of climatic, hydrological and vegetation cover variables [20, 21]. The modelling framework used, known as predictive habitat distribution or niche modelling, has rapidly increased in many fields such as biogeography, evolution, ecology, epidemiology, conservation, and invasive-species management . These models associate species and even community occurrences with environmental variables (biotic or abiotic) to predict their potential geographic distribution . The idea is that known occurrences of species across landscapes can be related to raster geographic information system coverages summarizing environmental variation across those landscapes to develop a quantitative picture of the ecologic distribution of the species . Although climate is only one of several factors that influence the distribution of the colilargo, in the Argentine model, mean annual temperature and cumulative precipitation resulted in the best predictors of its spatial distribution . In warmer areas, the presence of colilargos seamed to be associated with lower than normal precipitation values (i.e., north-eastern Patagonia). In colder areas, presence was associated with higher precipitation (i.e., western Andes, close to the Chilean border). The model, being parametric in nature, allows the recalculation of the probability of colilargos presence with different parameters of temperature and precipitation.
The Patagonian region, as other regions of the world, has experienced important climate variations during the last century (Figure 2). In particular, significant positive temperature trends in association with negative precipitation trends have been documented along the southern Andes [25, 26]. In addition, recent studies estimate climate change projections over South America for the twenty-first century  based on information provided by global climate simulations performed under greenhouse gas increment scenarios for the fourth assessment report of the Intergovernmental Panel on Climate Change (IPCC) . There is a general consensus among global climate models that climate changes projected over the southern Andes for the last part of the twenty-first century will be mainly associated with temperature increases and reductions in precipitation during all seasons. Nevertheless, Vera et al. , using an ensemble of 21 global climate model simulations, showed that the quantification of such regional potential climate change in the southern part of South America exhibits a large degree of uncertainty that inhibits its use in impact assessment studies. Therefore, the aim of our work is to evaluate the changes in the probability of the colilargo presence that might take place under four different potential climate change scenarios. In particular, we test the hypothesis that the future distribution of the hantavirus reservoir in Argentine Patagonia will be affected if future climate changes of similar magnitude and sign than those already observed in the twentieth century occur. We also examine the sensitivity of the potential distribution of the rodent to the uncertainties associated with those climate change scenarios.
The probability of the colilargo presence in a given place (presence probability) was calculated with a validated spatial distribution model which is as a function of temperature and precipitation . To evaluate changes in colilargo distribution its presence probability under current climate was compared to the probability modelled considering four possible future climate change scenarios. The four climate change scenarios were constructed using potential temperature and precipitation changes based on the linear extrapolation into the future of the changes recorded for the period between 1967 and 1998 in the Argentine Patagonia.
The colilargo distribution model for the Argentinean Patagonia is described in detail elsewhere . Briefly, the distribution model consisted of a generalized linear model , using logistic multiple regression. The probability of a site having colilargos (pO) was modelled as an S-shaped curve when LP is a first-order polynomial:
The linear predictor (PL) was related to the response variable (rv) through the logit link ln(rv/(1-rv)), which linearizes the response and restricts its values to lie between 0 and 1 . The LP was expressed as a linear combination of the environmental explanatory variables (xi): LP = a + bx1 + ... + cxi + e, with a, b, c being the parameters, and e a binomial error term with variance equal to μ (1-μ). The continuous environmental variables tested were mean annual temperature, mean annual cumulative precipitation, mean annual cumulative evapotranspiration, distance to permanent rivers, to any kind of water body, to roads, annual frost-day frequency, elevation above sea level, percent cover by trees, by grass, and by bare soil. The two factor variables tested where the phytogeographic classification of Argentina  and a detailed phytogeographic classification of Patagonia . The selection procedure was a manual upward stepwise selection based on a t-test (parameter/standard error with null model degrees of freedom) of each variable with α = 0.01 for retention because of the large number of variables considered. The final model included temperature (tt) and precipitation (pp) as the best predictors and the parameters after bootstrap re-sampling were:
with txp meaning both variables interaction. The equation was applied to each cell of the study area map with the geographic information system Arcview 3.2 . After transforming the variable with eq (1) the probability of the colilargo being present (presence probability) was obtained. Values lower than 0.5 indicated absence and values above or equal to 0.5 indicated presence. The model explained 59% of the total variance in presence and the external validation resulted in an 86% correct classification of the colilargo records .
Climate change scenarios
Considering that the best current distribution model was a function of temperature and precipitation, the climate change scenarios were built considering these two variables. Furthermore, as it was mentioned before, there is currently a lack of reliable quantification of future climate change at regional levels . In that sense, considering that IPCC reports and publications like  estimate that climate changes projected in Patagonia for the twenty-first century will be of the same sign than those observed in the second portion of the twentieth century, climate change scenarios were built considering the climate trends observed in the region during the last portion of the twentieth century.
Monthly mean surface air temperature from the Climatic Research Unit dataset  along with a gridded dataset of monthly rainfall means derived from meteorological stations were used to create the model. The linear trends of both variables were computed over the period 1967–1998 using a least squares technique. That period is long enough to estimate trends besides the climate variations associated with the natural climate variability. Those trends were mainly characterized by a generalized warming over the whole region with peaks along the Andes region (Figure 2a), and a significant decrease of precipitation along the Río Negro and Neuquén mountains (Figure 2b).
As mentioned in the Introduction, current global climate models project for the twenty-first century under a greenhouse gas increment scenario trends in both temperature and precipitation of the same sign as those observed during the last part of the twentieth century. Nevertheless, those climate models do not provide a realistic quantification of such change, due to the uncertainties inherent in such projections . Therefore, climate change scenarios for the region under study were estimated assuming that climatic trends observed over the period 1967–1998 would continue into the future. Based on the assumption of continued climate change, four different climate change scenarios were constructed: Scenario 1 assumed no change in precipitation but a future temperature change equal to that observed between 1967 and 1998; Scenario 2 assumed no change in temperature but a future precipitation change equal to that observed between 1967 and 1998; Scenario 3 assumed changes in both variables equal to those observed between 1967 and 1998; Scenario 4 assumed changes in both variables but doubled the values observed between 1967 and 1998.
Equations (1) and (2) of the distribution model were calculated cell by cell based on each scenario over the whole of Patagonia (square cells of approximately 19 km side). This created four maps of colilargo presence probability. To obtain a clearer view of presence probability change, the original map (obtained from 1967–1998 climate conditions) was subtracted from each scenario map (scenario model pO minus original model pO). This substraction of probabilities rendered 4 maps showing the "change in presence probability," positive values indicated an increase in presence probability according to the scenario and negative values a decrease.
To test the hypothesis of no change in presence probability due to the postulated climate change in the whole Patagonia, 230 random sites were sampled from the original model and from the four scenarios of climate change. A pair wise comparison between the original model and each scenario was used to assess significant differences in presence probability (Wilcoxon Rank Sum test). To determine if these changes could be spatially heterogeneous, ten homogeneous sub-areas were arbitrarily selected to evaluate local changes (Figure 1). Changes were analysed with a test for differences of proportions . The sub-areas were chosen in order to contain hantavirus transmission zones (1 & 6), zones without transmission but with the rodent widely present (3, 5, 7, 9), zones without transmission but with the rodent sparsely present (4, 8) and zones without transmission or rodent presence (2 & 10). The data about transmission corresponds to records between 1989 and 2002. After that date only the hospital of treatment was available, but potential place of infection and place of residence of the patient were not. The location of 59 cases out of 76 could be geolocated according to the potential place of infection. Sites separated by less than 5 km were clustered together resulting in 12 separate points. Finally, total areas of potential presence were compared between the original model and the significant scenarios. It was assumed that pixels with presence probability equal or higher than 0.5 indicated potential presence of the rodent, and those below 0.5 indicated potential absence. A contingency table was built with the number of pixels potentially present and potentially absent in the original model versus the scenario model. The differences were assessed with a χ2 test.
The colilargo presence probabilities were different among the four climate change scenarios, as shown in Figure 3. The presence probability based on Scenario 1 decreased in western Chubut and Santa Cruz provinces. Scenarios 2 and 3 also showed a decrease in those same areas as well as in western Neuquén. A rise in probability was observed in Scenarios 2 and 3 in a small patch in south-western Santa Cruz. On the other hand, Scenario 4 showed the greatest changes in presence probability. These changes occurred in the Andes range and in a fringe crossing the northern Patagonia from northwest to centre east. Rodent presence probability increased to 0.75 in patches aligned north-south from Neuquén to southern Santa Cruz. The similarities observed between scenarios 2 and 3 suggest that changes in precipitation may have a stronger effect on presence probability than temperature alone (Scenario 1).
According to the randomly sampled points only scenario 4 presented significant changes in presence probability (Z = -10.773, p < 0.0000). Thus the hypothesis of no change in presence probability is rejected for Scenario 4, which assumes a doubling of values in the trend observed between 1967 and 1998. The analysis for the sub-areas showed only significant reductions in presence probability (table 1). Scenario 4 showed most of the changes, with 7 out of 10 areas exhibiting significant probability decreases. Under this scenario the transmission area in Río Negro experienced a decrease in its colilargo presence probability by 87% (actual mean 0.997; Scenario 4 mean 0.129). Curiously, the only subset area that showed significant changes in another scenario was the transmission zone in Neuquén, where a fall of nearly 50% in colilargo presence probability was observed for scenarios 2 and 3 (actual mean 0.812; Scenario 2 mean 0.381 and Scenario 3 mean 0.412).
According to Scenario 4, the total surface of potential presence and potential absence of the rodent changed significantly (χ2 = 247.02, df = 1, p < 0.0001). 418,000 km2 would change from potential presence to absence and 82,000 km2 from potential absence to presence. This would imply a net reduction in potential presence surface of 336,000 km2. On the other hand, 120,000 km2 of potential presence and 1,106,000 of potential absence will not change their state (72% of the surface in total).
In this work, future potential distribution changes of the geographical distribution of a hantavirus reservoir, colilargo, were built assuming different climate change scenarios. These scenarios were estimated by extrapolating the climate trends observed in the Argentine Patagonia over the last part of the twentieth century. We used scenarios with spatially heterogeneous changes instead of a constant temperature or precipitation change for the whole study area. This proved to be important because a single value for a whole area might hide important local clusters of high presence probability. For example, using the mean precipitation change for Patagonia (20 mm) would have concealed rises of 380 mm and falls of nearly 900 mm in certain areas. These local changes were important, because they coincided with known areas of high transmission.
The starting point of this work was based on historical records of the species recorded during a long time span; we used the best information available about the rodent and the HPS cases. The temperature and precipitation data used therein did not exactly match the records compilation time. The model also assumed a static distribution of the species (i.e., in equilibrium with the environment), an assumption which is usually permitted for modelling purposes in spite of the large spatial extent of the study . It would be desirable to increase the temporal resolution of the rodent monitoring [12, 13, 37, 38] and to incorporate more variables in a study of transmission risk, including the relation between climate, fire, vegetation growth and seedling, rodent population dynamics  and hantavirus cases directly . But this is a difficult task to achieve in widespread areas and into the future, as it involves sequential relations (i.e., rain, flowering, seed release, population explosions) combined with spatially separated events (fires, cattle ranching, mast year) and many interactions. For example, in South America bamboo blooming events happen regularly, they occur every 12–14 years in a region but only repeat every 60 years in a certain spot . They not only influence rodent outbreaks but also affect the establishment of woods through favouring of fires and herbivory of the rodents on the seedlings . Also in southern Chile wood exploitation and plantations have produced landscape fragmentation which favours bamboo growth . Until spatio-temporal models together with extensive sampling are put together to include these variables and interactions we believe a summarizing model like the one presented herein is a step forward.
Colilargos are known to inhabit subantarctic forests, where woodlands are abundant, but they also occupy the steppe along scrublands adjacent to streams and roads .
Precipitation has a strong influence on the distribution of these woodlands, which might explain the observed subtle dominance of precipitation over temperature in influencing the distribution of the rodent. In Scenario 4 the higher presence probability areas for colilargos are predicted to shift eastwards, while the areas in the north eastern Patagonian steppe would become potential absence. This result can be explained by the combined changes of temperature and precipitation. In the western Andes temperature would rise and precipitation fall, probably diminishing the suitability of the environment for the rodent. However, to the east of these areas, the effects would result in the opposite (less augmentation of temperature and a rise in precipitation).
The relation between the colilargo distribution and HPS transmission is not straightforward. Colilargo populations experience dramatic outbreaks associated to bamboo blooming and mast seeding in southern Chile and Argentina; thus bamboo flowering may be used as a signal for forecasting colilargo outbreaks to predict transmission risk . But although in Chile these outbreaks were associated with an increased risk of HPS , in Argentina no cases were detected during an outbreak of colilargos in 1997 . In this same study there was no apparent correlation of rodent antibody prevalence with population density, or of rodent population density or antibody prevalence with numbers of human cases. Regarding human exposition risk, in Europe milder weather has been mentioned as a possible cause of more outdoor activities (recreational and occupational) and thus higher human exposure to infectious rodent excretions . Provided this temporal relation also held spatially, one could expect higher exposition in warmer areas.
In northern Argentina rural occupation was associated with higher hantavirus antibody prevalence, specially agricultural activities . But in the four corners U.S.A. most exposures occurred in or around the home . If we only consider studies about the Andes Sout lineage, in Chile the higher number of rodents positive for the virus were captured also in the domiciles and their surroundings . In southern Argentina seropositive colilargos were found in peridomestic rural environments  and a higher risk was anticipated in rural areas due to the lower trap success in the periurban areas . Also the subantarctic forest was the only ecosystem that harboured seropositive rodents in sylvan areas [37, 47]. Therefore the probability of humans contracting HPS may be greater in forest habitats, which is actually where most HPS cases have occurred. The distribution of reservoir species indicate the maximum potential extent of HPS  and in Argentine Patagonia HPS transmission zones matched the areas with the highest probability of colilargo presence . These findings suggest that presence probability of colilargos may indicate an approximate risk of transmission.
Our results suggest that if climate trends in the twenty-first century continue at the same rates as observed during the last portion of the twentieth century, significant changes in colilargo presence probability would be expected in 30 years from now in certain areas, and more extensively in 60 years (28% of the study area). The condition found with the Scenario 4 may be seen as either a steady maintenance of trends for the future 60 years or their doubling for the next 30 years. These results provide an estimate of the sensibility of the model, and suggest possible time periods that may be needed to observe changes in the rodent distribution under the assumptions of this work.
Also, the decline in presence probability of colilargos in the areas with HPS transmission and the net reduction in total area of potential presence suggest that future changes in Patagonia climate may lower transmission risk. In IPCC , the possibility of the incidence of some diseases decreasing due to climate change is mentioned only marginally. Most of the attention is centred on diseases that would increase their risk of exposure. Here we present a situation, in which, future potential climate change might reduce the extent of the distribution of a hantavirus reservoir and maybe the risk of its transmission.
Provided that temperature and precipitation follow the trends observed between 1963 and 1998 in Argentine Patagonia, the presence probability of O. longicaudatus would drop in some areas. Given that changes maintain for 60 years or double in 30 years, the hantavirus reservoir would contract its distribution, and although favourable areas might shift position eastwards, the net total area of potential presence would also fall.
Palma RE, Rivera-Milla E, Salazar-Bravo J, Torres-Pérez F, Pardiñas UFJ, Marquet PA, Spotorno AE, Meynard ASP, Yates TL: Phylogeography of Oligoryzomys longicaudatus (Rodentia: Sigmodontinae) in temperate South America. J Mammal. 2005, 86: 191-200. 10.1644/1545-1542(2005)086<0191:POOLRS>2.0.CO;2.
López N, Padula PJ, Rossi C, Lázaro ME, Franze-Fernandez MT: Genetic identification of a new hantavirus causing severe Pulmonay Syndrome in Argentina. Virology. 1996, 220: 223-226. 10.1006/viro.1996.0305.
Toro J, Vega JD, Khan AS, Mills JN, Padula P, Terry W, Yadón Z, Valderrama R, Ellis BA, Pavletic C: An Outbreak of Hantavirus Pulmonary Syndrome, Chile, 1997. Emerg Infect Dis. 1998, 4: 687-694.
Yahnke CJ, Meserve PL, Ksiazek TG, Mills JN: Patterns of infection with Laguna Negra virus in wild populations of Calomys laucha in the central Paraguayan chaco. Am J Trop Med Hyg. 2001, 65 (6): 768-776.
Puerta H, Cantillo C, Mills J, Hjelle B, Salazar-Bravo J, Mattar S: Hantavirus del nuevo mundo. Ecología y epidemiología de un virus emergente en Latinoamérica. Medicina Buenos Aire. 2006, 66: 343-356.http://www.scielo.org.ar/scielo.php?script=sci_abstract&pid=S0025-76802006000400013&lng=es&nrm=iso&tlng=en
Kosoy MY, Slonova RA, Mills JN, Mandel E, Childs JE: Community structure and prevalence of hantavirus infection in rodents: a geographic division of the enzootic area in far eastern Russia. J Vector Ecol. 1997, 22: 52-63.
Githeko A, Lindsay S, Confalonieri U, Patz J: Climate change and vector-borne diseases: a regional analysis. B World Health Organ. 2000, 78: 1136-1147.
Hales S, Wet Nd, Maindonald J, Woodward A: Potential effect of population and climate change on global distribution of dengue fever: an empirical model. Lancet. 2002, 360: 830-834. 10.1016/S0140-6736(02)09964-6.
Hopp M, Foley J: Global-scale relationships between climate and the dengue fever vector, Aedes aegypti. Climatic Change. 2001, 48: 441-463. 10.1023/A:1010717502442.
Martens W, Jetten T, Focks D: Sensitivity of malaria, schistosomiasis and dengue to global warming. Climatic Change. 1997, 35: 145-156. 10.1023/A:1005365413932.
Peterson AT, Shaw J: Lutzomyia vectors for cutaneous leishmaniasis in Southern Brazil: ecological niche models, predicted geographic distributions, and climate change effects. Int J Parasitol. 2003, 33: 919-931. 10.1016/S0020-7519(03)00094-8.
Murua R, Gonzalez LA, Lima M: Population dynamics of rice rats (a Hantavirus reservoir) in southern Chile: feedback structure and non-linear effects of climatic oscillations. Oikos. 2003, 102: 137-145. 10.1034/j.1600-0706.2003.12226.x.
Murua R, Briones M: Abundance of the sigmodontine mouse Oligoryzomys longicaudatus and patterns of tree seeding in the Chilean temperate forest. Mamm Biol. 2005
Aagesen D: Crisis and conservation at the end of the world: sheep ranching in Argentine Patagonia. Environ conserv. 2000, 27: 208-215. 10.1017/S0376892900000229.
González ME, Veblen TT, Sibold JS: Fire history of Araucaria-Nothofagus forests in Villarrica National Park, Chile. J Biogeography. 2005, 32: 1187-1202. 10.1111/j.1365-2699.2005.01262.x.
Kitzberger T, Veblen TT: Fire-induced changes in northern Patagonian landscapes. Landscape Ecol. 1999, 14: 1-15. 10.1023/A:1008069712826.
Tercero-Bucardo N, Kitzberger T, Veblen TT, Raffaele E: A field experiment on climatic and herbivore impacts on post-fire tree regeneration in north-western Patagonia. J Ecol. 2007, 95: 771-779. 10.1111/j.1365-2745.2007.01249.x.
Jaksic F, Lima M: Myths and facts on ratadas: Bamboo blooms, rainfall peaks and rodent outbreaks in South America. Austral Ecol. 2003, 28: 237-251. 10.1046/j.1442-9993.2003.01271.x.
Engelthaler DM, Mosley DG, Cheek JE, Levy CE, Komatsu KK, Ettestad P, Davis T, Tanda DT, Miller L, Frampton JW: Climatic and Environmental Patterns Associated with Hantavirus Pulmonary Syndrome, Four Corners Region, United States. Emerg Infect Dis. 1999, 5: 87-94.
Carbajo AE, Pardiñas UFJ: Spatial distribution model of the Hantavirus reservoir, the long-tailed colilargo (Oligoryzomys longicaudatus), in Argentina. J Mammal. 2007, 88: 1555-1568. 10.1644/06-MAMM-A-183R1.1.
Porcasi X, Calderón G, Lamfri M, Scavuzzo M, Sabattini MS, Polop JJ: Predictive distribution maps of rodent reservoir species of zoonoses in Southern America. Mastozool Neotrop. 2005, 12: 199-216.http://www.sarem.org.ar/archivos/revistas/12_2/porcasi.pdf
Anderson RP, Lew D, Peterson AT: Evaluating predictive models of species' distributions: criteria for selecting optimal models. Ecol Model. 2003, 162: 211-232. 10.1016/S0304-3800(02)00349-6.
Guisan A, Zimmermann NE: Predictive habitat distribution models in ecology. Ecol Model. 2000, 135: 147-186. 10.1016/S0304-3800(00)00354-9.
Peterson AT: Ecologic Niche Modeling and Spatial Patterns of Disease Transmission. Emerg Infect Dis. 2006, 12: 1822-1826.
Haylock MR, Peterson T, Sousa JRAd, Alves LM, Ambrizzi T, Baez J, Barbosa JI, Barros VR, Berlato MA, Bidegain M: Trends in total and extreme South American rainfall in1960–2000 and links with sea surface temperature. J Climate. 2006, 19: 1490-1512. 10.1175/JCLI3695.1.
Rusticucci M, Penalba O: Interdecadal changes in the precipitation seasonal cycle over Southern South America and their relationship with surface temperature. Cli Res. 2000, 16: 1-15. 10.3354/cr016001.
Vera C, Silvestri G, Liebmann B, González P: Climate change scenarios for seasonal precipitation in South America from IPCC-AR4 models. Geophys Res Lett. 2006, 33: L13707-10.1029/2006GL025759.
I.P.C.C.: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. 2007, Cambridge: Cambridge University Press
Vera CS, Gonzalez P, Silvestri G: About uncertainties in WCRP/CMIP3 climate simulations over South America. Proceedings of the 9th International Conference on Southern Hemisphere Meteorology and Oceanography; February 9–13; Melbourne, Australia. 2009, 9-13. http://www.bom.gov.au/events/9icshmo/manuscripts/TH1100_Vera.pdf
McCullagh P, Nelder JA: Generalized Linear Models. 1989, London: Chapman & Hall
Crawley MJ: GLIM for ecologists. 1993, Oxford: Blakwell Scientific Publications, 1
Cabrera AL, Willink A: Biogeografía de América Latina. 1973, Washington D.C.: OEA
León RJC, Bran D, Collantes M, Paruelo JM, Soriano A: Grandes unidades de vegetación de la Patagonia extra andina. Ecol Austral. 1998, 8: 125-144.
E.S.R.I: Arcview. 2000, Redlands, California, U.S.A.: Environmental Systems Research Institute, 3.2
Jones PD, New M, Parker DE, Martin S, Rigor IG: Surface air temperature and its variations over the last 150 years. Rev Geophys. 1999, 37: 173-199. 10.1029/1999RG900002.
Zar J: Biostatistical Analysis. 1999, New Jersey: Prentice Hall
Cantoni G, Padula PJ, Calderon G, Mills J, Herrero E, Sandoval P, Martinez V, Pini N, Larrieu E: Seasonal variation in prevalence of antibody to hantaviruses in rodents from southern Argentina. Trop Med Int Health. 2001, 6: 811-816. 10.1046/j.1365-3156.2001.00788.x.
Murua R, Gonzalez LA, Lima M: Second-order feedback and climatic effects determine the dynamics of a small rodent population in a temperate forest of South America. Popul Ecol. 2003, 45: 19-24. 10.1007/s10144-003-0135-y.
Clement J, Vercauteren J, Verstraeten WW, Ducoffre G, Barrios JM, Vandamme AM, Maes P, Van Ranst M: Relating increasing hantavirus incidences to the changing climate: the mast connection. Int J Health Geogr. 2009, 8: 1-10.1186/1476-072X-8-1.
González ME, Veblen TT, Donoso C, Valeria L: Tree regeneration responses in a lowland Nothofagus-dominated forest after bamboo dieback in South-Central Chile. Plant Ecol. 2002, 161: 59-73. 10.1023/A:1020378822847.
Murua R, Padula P: Ecología y evolución de hantavirus en el Cono Sur de América. Arch med vet. 2004, 36: 1-20.http://www.scielo.cl/scielo.php?script=sci_abstract&pid=S0301-732X2004000100001&lng=es&nrm=iso&tlng=en
Murua R, Gonzalez LA: Microhabitat selection in two Chilean cricetid rodents. Oecologia. 1982, 52: 12-15. 10.1007/BF00349005.
Piechotowski I, Brockmann SO, Schwarz C, Winter CH, Ranft U, Pfaff G: Emergence of hantavirus in South Germany: rodents, climate and human infections. Parasitol Res. 2008, 103 (Suppl): S131-S137. 10.1007/s00436-008-1055-8.
Pini N, Levis S, Calderón G, Ramirez J, Bravo D, Lozano E, Ripoll C, Jeor SS, Ksiazek TG, Barquez RM, Enria D: Hantavirus Infection in Humans and Rodents, Northwestern Argentina. Emerg Infect Dis. 2003, 9:
Hjelle B, Glass GE: Outbreak of Hantavirus Infection in the Four Corners Region of the United States in the Wake of the 1997–1998 El Niño-Southern Oscillation. J Infect Dis. 2000, 181: 1569-1573. 10.1086/315467.
Calderon G, Pini N, Bolpe J, Levis S, Mills J, Segura E, Guthmann N, Cantoni G, Becker J, Fonollat A: Hantavirus reservoir hosts associated with peridomestic habitats in Argentina. Emerg Infect Dis. 1999, 5: 792-797.
Piudo L, Monteverde M, Gonzalez Capria S, Padula PJ, Carmanchahi P: Distribution and abundance of sigmodontine rodents in relation to hantavirus in Neuquén, Argentina. J Vector Ecol. 2005, 30: 119-125.
Mills JN, Childs JE: Ecologic studies of rodent reservoires: their relevance for human health. Emerg infect dis. 1998, 4: 529-537.
This work was partially funded by PICT 2004 no. 20790 from the Agencia Nacional de Promoción Científica y Tecnológica. The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under Grant Agreement N° 212492 (CLARIS-LPB Project).
The authors declare that they have no competing interests.
AEC and CV designed the study. CV and PLG prepared and processed climatic data and maps. AEC made the modeling. All authors interpreted results, read and approved the final manuscript.