Assessment of climate change impact on the malaria vector Anopheles hyrcanus, West Nile disease, and incidence of melanoma in the Vojvodina Province (Serbia) using data from a regional climate model
Authors:
Dragutin T. Mihailović aff001; Dušan Petrić aff002; Tamaš Petrović aff003; Ivana Hrnjaković-Cvjetković aff004; Vladimir Djurdjevic aff006; Emilija Nikolić-Đorić aff007; Ilija Arsenić aff001; Mina Petrić aff008; Gordan Mimić aff011; Aleksandra Ignjatović-Ćupina aff002
Authors place of work:
Department of Field and Vegetable Crops, Faculty of Agriculture, University of Novi Sad, Novi Sad, Serbia
aff001; Department of Plant and Environment Protection, Faculty of Agriculture, University of Novi Sad, Novi Sad, Serbia
aff002; Department for virology, Scientific Veterinary Institute “Novi Sad”, Novi Sad, Serbia
aff003; Institute of Public Health of Vojvodina, Novi Sad, Serbia
aff004; Faculty of Medicine, University of Novi Sad, Novi Sad, Serbia
aff005; Institute of Meteorology, Faculty of Physics, University of Belgrade, Belgrade, Serbia
aff006; Department of Agricultural Economics, Faculty of Agriculture, University of Novi Sad, Novi Sad, Serbia
aff007; Avia-GIS NV, Zoersel, Belgium
aff008; Department of Physics, Faculty of Sciences, University of Novi Sad, Novi Sad, Serbia
aff009; Department of Physics and Astronomy, Faculty of Sciences, University of Gent, Gent, Belgium
aff010; BioSense Institute, University of Novi Sad, Novi Sad, Serbia
aff011
Published in the journal:
PLoS ONE 15(1)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0227679
Summary
Motivated by the One Health paradigm, we found the expected changes in temperature and UV radiation (UVR) to be a common trigger for enhancing the risk that viruses, vectors, and diseases pose to human and animal health. We compared data from the mosquito field collections and medical studies with regional climate model projections to examine the impact of climate change on the spreading of one malaria vector, the circulation of West Nile virus (WNV), and the incidence of melanoma. We analysed data obtained from ten selected years of standardised mosquito vector sampling with 219 unique location-year combinations, and 10 years of melanoma incidence. Trends in the observed data were compared to the climatic variables obtained by the coupled regional Eta Belgrade University and Princeton Ocean Model for the period 1961–2015 using the A1B scenario, and the expected changes up to 2030 were presented. Spreading and relative abundance of Anopheles hyrcanus was positively correlated with the trend of the mean annual temperature. We anticipated a nearly twofold increase in the number of invaded sites up to 2030. The frequency of WNV detections in Culex pipiens was significantly correlated to overwintering temperature averages and seasonal relative humidity at the sampling sites. Regression model projects a twofold increase in the incidence of WNV positive Cx. pipiens for a rise of 0.5°C in overwintering TOctober–April temperatures. The projected increase of 56% in the number of days with Tmax ≥ 30°C (Hot Days—HD) and UVR doses (up to 1.2%) corresponds to an increasing trend in melanoma incidence. Simulations of the Pannonian countries climate anticipate warmer and drier conditions with possible dominance of temperature and number of HD over other ecological factors. These signal the importance of monitoring the changes to the preparedness of mitigating the risk of vector-borne diseases and melanoma.
Keywords:
Europe – Horses – Climate change – West Nile virus – Climate modeling – Melanomas – Mosquitoes – Serbia
Introduction
Climate change is referred to as “the biggest global health threat of the 21st century” [1].The analysis of outputs from all general circulation models (GCM) suggests that the countries of the Pannonian Plain, including Serbia, are facing significant impacts of climate change, affecting all aspects of human life [2]. The authors of the manuscript (meteorology, entomology, veterinary medicine, and public health experts), have been working together since 2003, promoting the idea of multisectoral collaboration before the One Health Concept was officially inaugurated in the USA in 2007 [3], and endorsed by the European Union (EU) [4] as well as prominent organizations such as the World Health Organization (WHO), Food and Agriculture Organization (FAO), and the World Organization for Animal Health (OIE) in 2018 [5].
In this paper, the authors collected and analysed observed data over 31 years and related a subset to outputs from a Regional Climate Model (RCM). Vector-borne diseases and melanoma are significant climate-driven threats for which risk sources can be clearly defined [6]. Moreover, both present progressively growing environmental threats to the animal as well as human health in the countries of the Pannonian Plane [7,8,9].
The biology and distribution of mosquito vectors and their capacity to transmit mosquito-borne diseases are dependent on many factors such as global trade and travel, urbanisation, habitat destruction, pesticide application, host density, and climate. Anopheles hyrcanus and Culex pipiens are mosquito species that are vectors of malaria and West Nile virus (WNV) disease, respectively, the two vector-borne diseases distributed worldwide [10,11]. Malaria was eradicated from Serbia and other Balkan states during the last century. However, the spreading of its vectors (Anopheles mosquitoes) and the re-emergence of the disease in Greece [12] pose a threat to the South East and Central Europe once again. In 2018, Serbia was the second European country (after Italy) most affected by WNV disease (415 reported cases with 35 fatal outcomes). In Europe, the total number of reported human autochthonous WNV infections in 2018 (n = 2,083) exceeded, by far, the total number from the previous seven years (n = 1,832). During the same transmission season, outbreaks of West Nile fever among equids increased by 30% compared to the number of outbreaks in 2017. In total, 285 outbreaks among equids were reported by the EU Member States in 2018 [13]. Current evidence suggests that inter-annual and inter-decadal climate variability have a direct influence on the epidemiology of vector-borne diseases, with temperature and relative humidity as the principal abiotic factors influencing the life-cycles of the mosquito vector, the pathogen, the host and the interactions between them [14,15].
Melanoma is a malignant disease that has experienced a significant increase in incidence during the last few decades all over the world [16]. Climate change impact on melanoma should be considered as a synergy of changes in UV radiation (UVR) due to stratospheric ozone depletion and the long-term increase of air temperature leading to more prolonged exposure of individuals to UVR doses and consequently to a higher risk of melanoma [17]. Melanoma mortality in the Vojvodina Province (northern Serbia) (VPS) within the period 1985–2004 shows an evident increase, placing it amongst the most vulnerable regions in the world. Thus, Jovanović et al. [7] estimated and made the list of mortality rates from malignant melanoma for males (age-standardised rate/100,000) in Europe (39 countries) for the year 2000, using European Network of Cancer Registries (ENCR) data. This list shows that the VPS is among the top eleven states (six of them have parts in the Pannonian Plane) listed as the most endangered.
In this study, devoted to revealing the potential impact of climate change on animal and human health, we compared a considerable amount of previously unpublished ecological data obtained from the field and clinical surveys with climate change projections for the VPS, which is representative of the Central European low-altitude areas with a human-dominated landscape (Fig 1). We examined the effects of temperature on the spread and relative abundance of the malaria vector An. hyrcanus and the “microclimate” differentiation between sites with a specific frequency of WNV occurrence in Cx. pipiens. We also evaluated the impact of climate change on melanoma incidence as a synergy of changes in UVR doses and the long-term increase in the number of hot days (HD), with daily maximum temperature ≥ 30°C using the Eta Belgrade University and Princeton Ocean Model (EBU-POM) regional model data.
Materials and methods
For the assessment of the climate change and the impact of UVR doses, we used the climatic variables obtained by the coupled regional EBU-POM model for the historical period 1961–2000 and the period 2001–2030 according to the A1B scenario defined in Special Report on Emissions Scenarios (SRES) [18], from now on SRES-A1B. SRES scenarios, which defined future global greenhouse gases emissions, were extensively used in the Intergovernmental Panel on Climate Change (IPCC) Third, and Fourth Assessment Reports. The main storyline behind the A1B scenario is rapid economic growth, followed by a significant increase in greenhouse gases concentrations in the future. In the Fifth Assessment Report (AR5), the Representative Concentration Pathways (RCPs) are introduced, which are possible future concentration pathways without any storyline behind them [19]. Comparing SRES-A1B and RCPs in terms of the greenhouse gases concentrations, at the end of this century SRES-A1B is the closest to RCP6.0, but for the time horizon used in this study, up to 2030, the difference between any SRES or RCPs are relatively small.
Study area and climate
The VPS is situated in the northern part of Serbia and the southern part of the Pannonian lowland (18°51′–21°33′E, 44°37′–46°11′N and 75–641 m.a.s.l. (with the Fruška Gora Mountain in the south) as it is seen in Fig 1A and 1B). This region is the essential food production area in Serbia with a total surface area of 21,500 km2 and a population of about 2 million. This region has a continental climate, with elements of a sub-humid and warm climate (Cfwbx” according to Köppen classification [20]).
Models and formula used
The global and regional climate model
For climate simulations in this study, we used results of the EBU-POM model runs for the SRES-A1B scenario integrated over the period 2001–2030 [21]. The EBU-POM is a two-way, coupled RCM. The atmospheric part is the Eta/National Centres for Environmental Prediction (NCEP) limited area model (resolution 0.25° × 0.25° on 32 vertical levels; centred at 41.5° N, 15° E, with boundaries at ±19.9° W–E and ±13.0° S–N), while the oceanic part is the POM (resolution 0.20° × 0.20° on 21 vertical levels). The driving global circulation model (GCM) was the ECHAM5 model [22] coupled with the Max Planck Institute Ocean Model (MPI-OM) [23]. More details about model integrations and performed bias correction for VPS can be found in the paper by Mihailović et al. [2]. The POM model was set over the Mediterranean Sea without the Black Sea; for other open seas, the sea surface temperature from the GCM was used as a bottom boundary condition.
Environmental sampling
Mosquito vectors
We used standardised protocols to measure mosquito presence/absence, density, and infestation by WNV. Data were extracted from dry ice-baited trap samples, collected over 31 years at 166 different sites (745 sampled locations, S1 Table) in the VPS, to infer on the trends of local vector status and virus circulation in mosquitoes. In all years, mosquitoes were sampled from May to September, with different spatial intensity and time-frequency governed by the scale and scope of different research projects. For comparison with climate variables, we extracted data obtained in 10 years (1985–6, 2004–5 and 2010–5) for which a standardised surveillance protocol was in place. These periods have the highest number of particular location-year combinations (S1–S3 Tables).
Samples were collected by two different types of dry-ice baited suction traps. During 1985 and 1986 [24,25] by the miniature CDC light trap (CDC) and for 2004 and 2015 by the NS2 trap (our design of dry ice-baited suction trap without light). Both traps were operating without a light source (incandescent light proved not to be attractive/repellent for most mosquito species inhabiting the VPS [25]). The CDC trap has 3–5 times stronger suction power (operated by a 9V battery) than NS2 (operated by 3 x 1.2V batteries), meaning that the increase in density of species observed after 1986 could not be attributed to the change of the type of trap. Traps were operated from the afternoon until the morning of the next day (one trap night), with different periodicity. The specific location of the trap at each site was chosen by experienced entomologists to stabilise variation of the collected data.
We used three parameters to indicate An. hyrcanus spread and population growth in the period 1985–2015: i) the ratio of positive to total mosquito samplings per year; ii) the number of sites invaded (positive places where it was looked for, but was not found in the preceding sampling period, and the number of sites where was observed in both periods, i.e. established); and iii) the average number of specimens sampled in one trap during single sampling period from the afternoon of the starting day to the morning of the next day (Fig 2A and 2C). Here, we used data from 1,073 mosquito samples (1985–6, 2004–5 and 2014–5), obtained at 54 location over six years (142 unique location-year combinations) (S2 Table).
For Cx. pipiens, the period starting with the first detection of WNV in mosquitoes in Serbia, in 2010 [26], to 2015 was considered. To investigate the impact of microclimate on the complex interaction between Cx. pipiens and WNV, we used the following climatic parameters from the EBU-POM model outputs (covering the period 2006–2015) for 11 sites (GPS coordinates–S3 Table) in the VPS with different histories of WNV circulation: (i) mean annual temperature (Ta); (ii) overwintering temperature (TOA) for the period October–April; (iii) seasonal temperature (TMS) for the period May–September; and (iv) relative humidity (RMS) for the period May–September. For these sites, we examined the correlation between the frequency of WNV detections in Cx. pipiens at each site (from 2010 to 2015) and the corresponding period averages of climate time series for the same site. For detection of WNV, specimens were sampled, anaesthetised by dry ice, identified to species level [27] on dry ice-cooled paper, pooled according to date, location, sex, and species, transported on dry ice to the laboratory, and stored at -70°C before virus detection. Pool size did not exceed 50 mosquito specimens per pool. Mosquito pools were tested for WNV RNA presence by TaqMan-based one-step reverse transcription real-time PCR (RT-qPCR) that amplified both lineage 1 and 2 strains. Viral RNA was extracted using the commercial ISOLATE II RNA Mini Kit (Bioline, The Netherlands) according to the manufacturer's instruction. One-step RT-qPCR was conducted using the commercial kit RNA UltraSense™ One-Step qRT-PCR System (Life Technologies Corporation) with the primers and probe that targeted the nucleocapsid protein C gene regions of WNV, as described by Petrović et al. [8]. We analysed the yearly occurrence of the WNV positive Cx. pipiens mosquitoes sampled by dry ice-baited traps in the years 2010–2015 across 66 unique location-year combinations (S3 Table). Only traps positioned precisely at the same spot over the entire six-year period are considered for analysis. Numbers allocated to different places (Fig 3) indicate the number of years in the period 2010–2015 in which WNV was detected in sampled Cx. pipiens mosquitoes; e.g. 5 indicates that WNV positive Cx. pipiens were detected in five out of the six years in the samples collected from the same spot.
Melanoma incidence
In the analysis we have used two indicators: (i) melanoma incidence rate that is a measure of the number of new cases ("incidence") per unit of time ("rate") and (ii) cumulative incidence (“incidence proportion” that measures the number of new cases per person in the population over a defined period of time–often called risk or proportion). Melanoma incidence rate (per 100,000 people) for ten years 1995–2004 was based on the data obtained from the paper by Jovanović et al. [7]. From these data, we calculated the cumulative incidence [28]. We have used the model simulation to study the expected impact of climate change on UVR exposure of human skin for nine sites in VPS [PA (Palić), SO (Sombor), KI (Kikinda), NS (Novi Sad), BC (Bečej), ZR (Zrenjanin), SM (Sremska Mitrovica), BK (Bantaski Karlovac), and BG (Beograd)]. For calculating the daily doses of UVR (UVRD) in the study area sites, we have used the following empirical formula UVRD = 0.002507×Gd−5.985 (kJ/m2) derived by Malinović-Milićević et al. [29], where Gd is the daily sum of global solar radiation. First, we calculated UVRD from global radiation model outputs using the empirical formula for the nine sites for the period April–September, and then we found the relative change R(UVRD) of those doses as R(UVRD) = (UVRD-UVRDk)/ UVRDk where UVRDk is the dose for 1961–1990 reference period, while the UVRD is calculated for the period 2001–2030.
Statistics
We considered the papers by Mihailović et al. [2,30] in which Kolmogorov complexity measures [Kolmogorov complexity (KC), Kolmogorov complexity spectrum (KC spectrum), and the highest value of the KC spectrum (KCM)], and sample entropy (SE) [31] were used to quantify the regularity and complexity of air temperature and precipitation time series, obtained by the EBU-POM model, representing both deterministic chaos and stochastic processes. We considered the complexity of the EBU-POM model using the observed and modelled time series of temperature and precipitation. We computed the KC spectrum, KC, KCM and SE values for temperature and precipitation. The calculations were performed for the entire time interval 1961–1990: (i) on a daily basis with a size of n = 10,958 samples for temperature and (ii) on a monthly basis with a size n = 360 for the precipitation. The simulated time series of temperature and precipitation were obtained by the EBU-POM model for the given period. The observed time series of temperature and precipitations for two stations: Sombor (SO) (88 m.a.s.l.) and Novi Sad (NS) (84 m.a.s.l.) in the considered area, were taken from daily meteorological reports of the Republic Hydrometeorological Service of Serbia. For both sites, the modelled complexity is lower than the observed one, but with the reliability which is in the interval values allowed by the information measures (KC, KCM, and SE) [32,33,34]. The term model reliability we have used in the following context. The Lyapunov exponent (LLE) relates to the predictability of measured time series, which includes deterministic chaos as an inherent component. Model predictability is here understood as the degree to which a correct prediction of a system's state can be made either qualitatively or quantitatively. In a stochastic analysis, a random process is considered predictable if it is possible to infer the next state from previous observations. In many models, however, randomness is a phenomenon which “spoils” predictability [35]. Deterministic chaos does not mechanically denote total predictability but means that at least it improves the prognostic power. In contrast, stochastic trajectories cannot be projected into the future. If LLE>1 then time series is not chaotic, but is rather stochastic, and predictions cannot be based on chaos theory. However, if 0<LLE<1 it indicates the existence of chaos in time series. In that case, one can compute the approximate time [often called Lyapunov time (LT)] limit for which accurate prediction for a chaotic system is a function of LLE. It designates a period when a specific process (physical, mechanical, hydrological, quantum, or even biological) moves beyond the bounds of precise (or probabilistic) predictability and enters a chaotic mode. According to Frison and Abarbanel [36] that time can be calculated as LLE Δtlyap = 1/LLE. If LLE→0, implying that Δtlyap→∞, then long-term accurate predictions are possible. However, many climate time series are highly complex. Therefore, Δtlyap can be corrected for randomness in the following way. Similar to Δtlyap we can introduce a randomness time Δtrand = 1/KC (in time units, second, hour or day). Henceforth, we shall denote this quantity Kolmogorov time (KT), as it quantifies the period beyond which randomness significantly influences predictability. Then, the Lyapunov time corrected for randomness is defined as [0,Δtlyap]∩[0,Δtrand]. It can be stated that the KT designates the size of the time window within time series where complexity remains nearly unchanged. These findings mean that the models with a KC (and KCM) complexity lower than the measured time series complexity cannot always reconstruct some of the structures contained in the observed data. However, it does not mean that outputs from EBU-POM model do not correctly simulate climate elements since both sites’ values indicate the absence of stochastic influences, providing reliable projections of the climate elements [37]. Therefore, the EBU-POM model can be considered as a model having high performances and reliability in projection of temperature and precipitation, two of the most important elements used in different climate research.
For An. hyrcanus, the temperature trend was evaluated by the Mann-Kendall test using the R statistical package [38]. Field observed values on species distribution and density for the period 1985–2015 and forecasts of the numbers of sites invaded and specimens sampled for the period 2016–2030 based on linear trends were obtained by the Eviews 9.5 software [39]. For Cx. pipiens, the relationship between yearly frequency of WNV detection in mosquitoes, air temperature, and relative humidity (derived from the climate model) was estimated using Spearman’s Rank-Order Correlation and a Poisson regression model (Statistica 13 [40]). For melanoma, the linear regression model was used for modelling the cumulative incidence of melanoma versus the difference of the cumulative UVR doses for hot and warm days (Fig 4D). Analysis of residual distribution was done by Shapiro-Wilk test (Statistica 13 [40])
Results
Mosquito vectors
Fig 2A shows an evident linear trend of the mean annual temperature Ta for the period 1985–2030 (r = 0.467; p = 0.001; τ = 0.328) calculated from the EBU-POM regional model outputs for 29 representative sites in the VPS. All parameters that were chosen for the evaluation of the spread and population increase of An. hyrcanus were positively, but to a different extent, correlated to the time argument (periods in which sampling was performed since the beginning of monitoring in 1985) indicating a monotonic trend. The increase of parameters follows the trend of Ta (Fig 2A). The strongest correlation was found for the increase in the ratio of positive samplings (r = 0.986; p = 0.000307; τ = 0.828), followed by the number of mosquitoes per trap night (r = 0.919; p = 0.009639;τ = 0.733), and the number of sites invaded (r = 0.889; p = 0.01766; τ = 0.6). By 2030 we anticipate a further increase in numbers of invaded sites and adult females sampled, by 1.71 and 1.27 fold, respectively (Fig 2B and 2C).
Spearman rank-order correlation between the frequency of WNV detections in Cx. pipiens at 11 sites and the corresponding mean values of climate time series was the highest for TOA (r = 0.755; p = 0.00008), then for Ta (r = 0.616; p = 0.00294), RMS (r = 0.499; p = 0.02119), and TMS (r = 0.477; p = 0.02856). Fig 3A depicts the Poisson regression model for the dependence of a number of detections per site (frequency - λ) on TOA, which is statistically significant (p = 0.01393). The output of the model (lnλ = −7.923+1.533×TOA) indicates that for an increase of 0.5°C in TOA (presuming that all other factors needed for the circulation of WNV are kept constant), a twofold increase in the incidence of WNV positive Cx. pipiens could be projected. Fig 3B depicts that most of the sites with the high frequency of WNV occurrence (≥2) were distributed along the northwest-southeast axis of the VPS.
Melanoma incidence and UVR doses
Fig 4B shows the positive relative change of UVRD, remarkably covering an eastern, southern, western, and partly central area of VPS. Specifically, the projected increase is twofold going from the west and northwest (0.60%) towards the east and southeast where it reaches values of about 1.20%. The EBU-POM model (for nine sites) shows a significant expected increase of 56% in the number of HD days in the VPS (Fig 4A), compared to the period 1961–1990. Additionally, we observed a decrease of 1.1% in the number of days with maximum air temperature Tmax ≥ 25°C (Warm Days—WD). This prolongs the exposure of outdoor working adults to UVR and thus leads to an increase in melanoma risk. This risk becomes even more significant because of the increase in cumulative values of UVR doses (Fig 4C). Fig 4D depicts the cumulative incidence of melanoma for the period 1985–2004 with an increasing monotonic trend (r = 0.9712 and p = 0.000003).
Discussion
Here we presented an intriguing comparison of the impact of climate change on complex systems including mosquito vectors, pathogens, and melanoma, which are all indicators of the risk imposed on human health. Our objectives were to use historical, previously unpublished sets of entomological and published clinical data and examine the importance of temperature in contributing to the spreading of the malaria vector An. hyrcanus; to differentiate between sites with a specific frequency of WNV occurrence in Cx. pipiens, and to assess the impact of increasing UVR and HD on melanoma incidence using the EBU-POM regional model data. A similar approach was recently used in observing the dramatic decline in total flying insect biomass in protected areas in Germany [41].
Despite globalisation trends, researchers have become "closed" in their ever-smaller communication circles which are not limited by state or national borders but by the professional language and way of thinking. Thus, by the end of the 20th century, the scientific community has been faced with problems in communication within its confines. One of the principal reasons why vector-borne diseases (VBD) are so difficult to predict, is the complex interaction of multiple factors (vector, host, pathogen, environment including short-term weather patterns and long-term climate change) in space and time [42,43]. Only groups from multiple sectors that communicate and work together on specific aspects of VBD systems will be able to answer the most exciting and pressing problems in the field [42]. Authors of this paper started collaboration in 2003 comparing the climates of the foci of WNV circulation in USA (California Central Valley) and Europe (Bucharest area) with VPS. As compared climates showed quite similar patterns, colleagues from public health and veterinary joined the initial group of meteorologists and medical entomologists. With the idea to better draw upon the resources and insights of the various sectors we designed and implemented research and programmes to achieve better outcomes in the control of zoonoses (diseases that can spread between animals and humans, e.g., WNV disease). This led us to the following achievements: (i) the first detection of WNV in horses in Serbia in 2009 [44]; (ii) the first detection of WNV in mosquitoes in Serbia in 2010 [26]; (iii) the first detection of WNV in wild birds in Serbia in 2012 [45]; (iv) development and implementation of the national programme of WNV surveillance in mosquito, bird and horse population [8], combined with human surveillance in VPS from 2014; (v) increased visibility to ECDC, EFSA and WHO; (vi) the first detection of imported dengue human case in Serbia in 2016 [46]; and (vii) development and implementation of “One Health” programme in VPS from 2018. We are quite sure that much less would have been achieved without multidisciplinary communication and collaboration initiated in 2003, and this paper would not have been compiled.
The temperature trend over the period of observations used in this study and for the future time horizon following A1B scenario obtained with the EBU-POM regional climate model is within the multi-model ensemble (MME) spread of regional climate models with the similar configuration used in the ENSEMBLES project [47]. For the region of interest, temperature change for the period 2001–2030 with respect to the period 1961–1990, in the EBU-POM integration is 0.75°C and in the case of ENSEMBLES MME spread range is 0.5–1.5°C [48]. Furthermore, over period 2001–2030, uncertainty in temperature projections, related to the scenarios is just small fraction of total uncertainty [49], due to the fact that there is no significant difference in greenhouse gases concentrations for different SRES scenarios. Consequently, it can be considered that for the period 2001–2030, in case of other scenarios, the estimated range will be similar. Following this findings, even that selection of single model and single scenario does not allow estimate of the full uncertainty in the future, presented results that relay on temperature change, can be seen as plausible future realisation within total uncertainty range, and realistic response of the complex system for further increase in temperature as a forcing driver. In future studies, it will be beneficial to introduce results of other regional climate models, such as multi-model ensemble from the CORDEX framework, as a need for better understanding of different uncertainties in the results.
Mosquito vectors
Until the end of the 20th century, northern Serbia was considered the northern limit for the distribution of An. hyrcanus in Europe. The first detection in Serbia dates from 1979 [50] from the north part of VPS. We found it in the central part of the Province in 1985 [25] and since then have been noticing its continued spread. Several records north from Vojvodina, in Slovakia in 2004 [51], the Czech Republic in 2005 [52], and Austria in 2012 [53] confirm our observation. Due to its exophilic and exophagic behaviour, An. hyrcanus has never been considered as the primary vector of malaria in Europe. Its spread to higher latitudes, combined with the changes in human behaviour (increased outdoor leisure activities, the mobility of humans, number of seasonal workers in the field, number of migrants in Europe), might increase its vector capacity. The similar northern spread of population distribution range that was registered for An. maculipennis s.s. in Russia [54], and Culiseta longiareolata in southern (in 2012 [55]) and northern (in 2013 [56]) Austria might well represent the tendency described with our model.
The latest illustration of similar changes is the finding of Uranotaenia unguiculata, a thermophilic mosquito species frequently occurring in the Mediterranean basin, in northern Germany, some 300km north of the previous northern limit [57].
During the period 2001–2030 in which the spread and population growth of An. hyrcanus is expected, the intensity of UVR is likely to increase in the VPS (Fig 4A). Positive trends which are present in our observations might indicate that the findings on the negative influence of UVR and blue-light radiation (this radiation has a wavelength between approximately 380nm and 500nm; it has a very short wavelength, and so produces a higher amount of energy) on adult mosquitoes under laboratory conditions [58,59] could not be simply translated to the field. This experimental evidence does not mean unavoidably that the blue light radiation has a significant influence on adult mosquitoes in field conditions since they can actively escape over-exposure to radiation.
The WNV transmission cycle involves mosquito vectors and birds, but equines and humans are also susceptible to infection [60,61]. Although WNV infections have been described in a wide variety of vertebrates, birds are the main natural reservoir. Hundreds of wild and domestic avian species have been described as susceptible to WNV infection, but many of these showed only subclinical infection [62]. In horses, WNV infection is also frequently clinically unapparent, but around 10% of cases develop neurological disorders with up to 50% mortality rates [61,63]. An increasing number of severe outbreaks in horses have been reported in Europe in the past decade, including a large one that took place in northeast Italy in 2008 involving 251 stables with 794 cases and five deaths [63].
From the first detection of WNV in 8 out of 81 dead wild birds found in Serbia [45], each year WNV nucleic acid was detected in dead or captured wild birds during summertime [8]. Serological testing of horses sampled during 2009–2010 showed that 12% of 349 horses from the northern part of Serbia had neutralising WNV antibodies [44]. After that, each year horse WNV cases were detected by the positive serological response (IgG and IgM antibody seroconverted horses) [8] or as a clinical manifestation of West Nile neuroinvasive disease [64].
A positive association between WNV disease and temperature was already reported in Europe [15,65] where climate and landscape were critical predictors of WNV disease outbreaks [66]. Our focus was not on the number of human WNV cases, but the suitability of sites/microhabitats with different air temperatures for WNV circulation in mosquitoes, which may well correspond to a higher risk of transmission. We found that sites with higher TOA and Ta were characterized with the higher frequency of WNV presence in mosquitoes. Clustering of cases with an incidence higher than one in six years coincided with an area of a significant grouping of mosquito, bird, horse, and human cases in 2014 and 2015 [9] (Fig 5). This is in concurrence with Tran et al. [67] and Marcantonio et al. [66], who found that average summer temperatures are positively correlated with WNV human incidence. It seems that temperature in semi-urban areas is an essential environmental factor influencing WNV circulation (landscape suitability for reservoir host and mosquito vector, host availability, and precipitation/water availability are somewhat similar in investigated semi-urban areas of VPS), as it affects both mosquito vector abundance [9] and virus replication. Prediction of a two-fold increase in virus incidence for each 0.5°C increase in TOA indicates but does not necessarily means, that the number of human cases could increase too. Therefore, our findings support the statement that climate change is likely to intensify the re-emergence of WNV in Europe [68].
Melanoma incidence and UVR
According to World Health Organization (WHO) (1992) and many other authors [69,70], exposure to UVR radiation is considered to be a major etiological factor for all three forms of melanoma (i) basal cell carcinoma (BCC), (ii) squamous cell carcinoma (SCC), and (iii) malignant melanoma (MM). We found the correlation between MM and climate changes impact on UVR and also the number of HD. We see the impact as a modification of ambient UVR through influences on other variables such as clouds and aerosols. However, that impact might be more pronounced through the impact of changes in outdoor ambient temperature which will influence people’s behaviour and increase the time they spend outdoors, i.e. exposure to both higher UVR and higher temperatures [17]. Experiments with animals clearly show that increased temperatures enhance UVR-induced melanoma compared to the room temperature. In an intriguing study, van der Leun [71] postulated that long-term elevation of temperature by 2°C, as a consequence of climate change, would increase the carcinogenic effects of UVR by 10%. Our results for the UVR in the VPS are generally similar to the ones obtained by Malinović-Milićević et al. [72] and Malinović-Milićević and Radovanović [73], who reported the following changes: (i) the reduction of yearly averages for the total ozone of 3.44% and 3.21%, and (ii) increase in erythemal UVR dose of 6.9% and 9.7% for the periods 1990–1999 and 2000–2009, respectively.
According to Jovanović et al. [7], the incidence rate of MM cancer in VPS for the period 1985–2004 is higher than in central Serbia and is comparable with the majority of the central European countries with the highest melanoma incidence rate in the world [74]. However, most studies do not deal more quantitatively with the relationship between UVR doses and exposure during HD days, and as it has been stated above, the cumulative exposure to sunlight is probably the most critical risk factor for MM and SCC cancers, while BCC is more associated with intensive short-term exposure [75]. Thus, the increasing trend in the number of melanoma incidence in the VPS for the period 1985–2004 (Fig 4D) can be ascribed to (i) the increase in the number of HD days for about 55% and (ii) the increase in cumulative values of UVR doses for the period 1985–2030. In a cohort study, Wu et al. [76] considered the impact of long-term UV radiation flux on skin cancer risk. Comparing with participants in the lowest quintile of cumulative UV radiation flux in adulthood, they found that participants in the highest quintile had multivariable-adjusted risks (cumulative incidence). According to Vandenbroucke and Pearce [28], some studies where cumulative incidence is used can over-represent the trends.
From a statistical point of view, the linear regression model for modelling the cumulative incidence of melanoma versus the difference of the cumulative UVR doses for hot and warm days (Fig 4D) is acceptable. Parameters are statistically highly significant (r = 0.9712 and p = 0.000003) while analysis of residual distribution shows a good agreement with the normal distribution (Shapiro-Wilk test, W = 0.9608, p = 0.7952).
We hope that our results will indicate the importance of long-term monitoring/surveillance programs for providing crucial data to evidence the ongoing biological alteration triggered by climate change. Nonetheless, it is difficult to say how broadly our data represent the trends elsewhere. We believe that the specificity of the observations offers a unique window into the state of some of the planet's pressing threats to human health. Also, in the case where humans are exposed to UVR, due to the nature of their work (the VPS is an exclusively agricultural area), it is necessary to (i) establish a broader network for UVR measurements and warning centres and (ii) increase the awareness of the melanoma as a result of increased amount of UVR.
Supporting information
S1 Table [xlsx]
Overview of dry-ice trap samples sizes.
S2 Table [xlsx]
Number of the total trap nights, positive trap nights and specimens sampled at 54 selected sites in the Vojvodina Province, Serbia during the years 1985–86, 2004–5 and 2014–15.
S3 Table [docx]
Frequency of sampling of WNV infected mosquitoes (1–5 times) in the Vojvodina Province, Serbia, during the period 2010–2016.
Zdroje
1. Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, et al. Managing the health effects of climate change: lancet and University College London Institute for Global Health Commission. The Lancet. 2009;373: 1693–1733.
2. Mihailović DT, Lalić B, Drešković N, Mimić G, Djurdjević V, Jančić M. Climate change effects on crop yields in Serbia and related shifts of Köppen climate zones under the SRES-A1B and SRES-A2. Int J Climatol. 2015;35: 3320–3334. doi: 10.1002/joc.4209
3. King LJ, Anderson LR, Blackmore CG, Blackwell MJ, Lautner EA, Marcus LC, et al. Executive summary of the AVMA One Health Initiative Task Force report. J Am Vet Med Assoc. 2008;233: 259–261. doi: 10.2460/javma.233.2.259 18627228
4. Health—European External Action Service [Internet]. [cited 29 Mar 2019]. Available: https://eeas.europa.eu/headquarters/headquarters-homepage/2348/health_ku
5. One World, One Health: OIE—World Organisation for Animal Health [Internet]. [cited 27 Mar 2019]. Available: http://www.oie.int/en/for-the-media/editorials/detail/article/one-world-one-health/
6. Kjellstrom T, McMichael AJ. Climate change threats to population health and well-being: the imperative of protective solutions that will last. Glob Health Action. 2013;6. doi: 10.3402/gha.v6i0.20816 23561024
7. Jovanović M, M Bogdanović G, Mijatović Jovanović V, Jeremić P, Ac Nikolić E. Analysis of cutaneous melanoma in the province of Vojvodina. J BUON Off J Balk Union Oncol. 2009;14: 441–6.
8. Petrović T, Šekler M, Petrić D, Lazić S, Debeljak Z, Vidanović D, et al. Methodology and results of integrated WNV surveillance programmes in Serbia. PLOS ONE. 2018;13: e0195439. doi: 10.1371/journal.pone.0195439 29624622
9. Petrić D, Petrović T, Hrnjaković Cvjetković I, Zgomba M, Milošević V, Lazić G, et al. West Nile virus “circulation” in Vojvodina, Serbia: Mosquito, bird, horse and human surveillance. Mol Cell Probes. 2017;31: 28–36. doi: 10.1016/j.mcp.2016.10.011 27777104
10. Zeller H, Lenglet A, Bortel WV. West Nile virus: the need to strengthen preparedness in Europe. Eurosurveillance. 2010;15: 19647. doi: 10.2807/ese.15.34.19647-en 20807490
11. World Health Organization, World Malaria Report 2018. Geneva: The Organization; 2018.
12. Vakali A, Patsoula E, Spanakos G, Danis K, Vassalou E, Tegos N, et al. Malaria in Greece, 1975 to 2010. Euro Surveill Bull Eur Sur Mal Transm Eur Commun Dis Bull. 2012;17.
13. Epidemiological update: West Nile virus transmission season in Europe, 2018 [Internet]. [cited 09Sep 2019]. Available from: https://ecdc.europa.eu/en/news-events/epidemiological-update-west-nile-virus-transmission-season-europe-2018
14. Wilson N, Lush D, Baker MG. Meteorological and climate change themes at the 2010 International Conference on Emerging Infectious Diseases. Euro Surveill Bull Eur Sur Mal Transm Eur Commun Dis Bull. 2010;15.
15. Calzolari M, Pautasso A, Montarsi F, Albieri A, Bellini R, Bonilauri P, et al. West Nile virus surveillance in 2013 via mosquito screening in northern Italy and the influence of weather on virus circulation. PLoS One. 2015;10: e0140915. doi: 10.1371/journal.pone.0140915 26488475
16. Apalla Z, Lallas A, Sotiriou E, Lazaridou E, Ioannides D. Epidemiological trends in skin cancer. Dermatol Pract Concept. 2017;7: 1.
17. Kerr JB, Seckmeyer G. Surface ultraviolet radiation: past and future Scientific Assessment of Ozone Depletion 2002. Geneva; 2003. Report No.: 47.
18. Summary for Policy makers. In: Nakićenović N, Swart R, editors. Special Report on Emissions Scenarios. Cambridge: Cambridge University Press—Published for the Intergovernmental Panel on Climate Change. 2000. pp. 1–21.
19. Moss R, Babiker M, Brinkman S, Calvo E, Carter T, Edmonds J et al. Towards new scenarios for analysis of emissions, climate change, impacts, and response strategies. IPCC Expert Meeting Report on New Scenarios. Intergovernmental Panel on Climate Change, Geneva. 2008. pp. 1–132.
20. Kottek M, Grieser J, Beck C, Rudolf B, Rubel F. World Map of the Köppen-Geiger climate classification updated. Meteorologische Zeitschrift. 2006;15(3):259–263.
21. Đurđević V, Rajković B. Development of the EBU-POM coupled regional climate model and results from climate change experiments. In: Mihailovic DT, Lalic B, editors. Advances in environmental modeling and measurements. New York, USA: Nova Science Publishers Inc; 2012. pp. 23–32.
22. Roeckner E, Bäuml G, Bonaventura L. The atmospheric general circulation model ECHAM5. Part I: Model description Report. Max Planck Institute for Meteorology; 2003 p. 127. Report No.: 349.
23. Jungclaus JH, Keenlyside N, Botzet M, Haak H, Luo J-J, Latif M, et al. Ocean Circulation and Tropical Variability in the Coupled Model ECHAM5/MPI-OM. J Clim. 2006;19: 3952–3972. doi: 10.1175/JCLI3827.1
24. Srdić Ž, Zgomba M, Petrić D. Les moustiques (Dip.Culicidae) et la demoustication en Vojvodina, Yugoslavia. IV Congrès sur la protection de la santé humaine et des cultures en milieu tropical. Marseille, France; 1986. pp. 489–495.
25. Petrić D. Seasonal and daily activity of adult mosquitoes (Diptera, Culicidae) in Vojvodina Province. Ph.D. Thesis, University of Novi Sad. 1989.
26. Petrić D, Hrnjaković-Cvjetković I, Radovanov J, Cvjetković D, Jerant-Patić V, Milošević V, et al. West Nile virus surveillance in humans and mosquitoes and detection of cell fusing agent virus in Vojvodina Province (Serbia). HealthMED 2012;6(2):462–68.
27. Becker N, Petrić D, Zgomba M, Boase C, Madon M, Dahl C, et al. Mosquitoes and their control. 2nd ed. Berlin Heidelberg: Springer-Verlag; 2010.
28. Vandenbroucke JP, Pearce N. Incidence rates in dynamic populations. Int J Epidemiol. 2012;41(5):1472–1479. doi: 10.1093/ije/dys142 23045207
29. Malinovic-Milicevic S, Mihailovic DT, Lalic B, Dreskovic N. Thermal environment and UV-B radiation indices in the Vojvodina region, Serbia. Clim Res. 2013;57: 111–121. doi: 10.3354/cr01163
30. Mihailović DT, Mimić G, Arsenić I. Climate Predictions: The Chaos and Complexity in Climate Models. In: Advances in Meteorology [Internet]. 2014 [cited 27 Mar 2019]. doi: 10.1155/2014/878249
31. Richman JS, Moorman JR. Physiological time-series analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol. 2000;278: H2039–2049. doi: 10.1152/ajpheart.2000.278.6.H2039 10843903
32. Krzic A, Tosic I, Djurdjevic V, Veljovic K, Rajkovic B. Changes in some indices over Serbia according to the SRES A1B and A2 scenarios. Climate Research. 2011;49: 73–86.
33. Cavicchia L, Scoccimarro E, Gualdi S, Marson P, Ahrens B, Berthou S, et al. Mediterranean extreme precipitation: a multi-model assessment. Clim. Dyn. 2016. doi: 10.1007/s00382-016-3245-x
34. Dell’ Aquila A, Mariotti A, Bastin S, Calmanti S, Cavicchia L, Deque M, et al. Evaluation of simulated decadal variations over the Euro-Mediterranean region from ENSEMBLES to Med-CORDEX. Clim. Dyn. 2016. doi: 10.1007/s00382-016-3143-2
35. Mihailović DT, Nikolić-Đorić E, Malinović-Milićević S, Singh VP, Mihailović A, Stošić T et al. The Choice of an Appropriate Information Dissimilarity Measure for Hierarchical Clustering of River Streamflow Time Series, Based on Calculated Lyapunov Exponent and Kolmogorov Measures. Entropy. 2019; 21: 215.
36. Frison TW, Abarbanel HDI. Ocean gravity waves: A nonlinear analysis of observations. Geophys Res. 1997; 102(C1): 1051–1059. doi: 10.1029/96JC02993
37. Mihailovic DT, Mimic G, Arsenic I. Climate predictions: the chaos and complexity in climate models, Advances in Meteorology. 2014; doi: 10.1155/2014/878249
38. R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna: R Foundation for Statistical Computing; 2016. Available: http://www.R-project.org/
39. HS Global Inc. Eviews [Internet]. CA, USA: HS Global Inc.; 2016. Available: http://www.eviews.com
40. TIBCO Software Inc. Statistica [Internet]. CA, USA: TIBCO Software Inc.; Available: https://docs.tibco.com/products/tibco-statistica-13-3-0
41. Hallmann CA, Sorg M, Jongejans E, Siepel H, Hofland N, Schwan H, et al. More than 75 percent decline over 27 years in total flying insect biomass in protected areas. PloS One. 2017;12: e0185809. doi: 10.1371/journal.pone.0185809 29045418
42. Moore CG. Interdisciplinary research in the ecology of vector-borne diseases: Opportunities and needs. Journal of Vector Ecology 2008;33(2):218–224. doi: 10.3376/1081-1710-33.2.218 19263839
43. Zimmerman B. Engaging with Complexity: Thrive! A Plan for a Healthier Nova Scotia. 2014; [e-print] Available from: https://thrive.novascotia.ca/sites/default/files/Thrive-Summit-2014-Brenda-Zimmerman-En.pdf.
44. Lupulović D, Martin-Acebes MA, Lazić S, Alonso-Padilla J, Blazquez AB, Escribano-Romero E, et al. First serological evidence of West Nile virus activity in horses in Serbia. Vector Borne Zoonotic Dis. 2011;11(9):1303–5.40. doi: 10.1089/vbz.2010.0249 21438694
45. Petrović T, Blazquez AB, Lupulović D, Lazić G, Escribano-Romero E, Fabijan D, et al. Monitoring West Nile virus (WNV) infection in wild birds in Serbia during 2012: First isolation and characterisation of WNV strains from Serbia. Eurosurveillance. 2013;18(44):1–8.
46. Petrović V, Turkulov V, Ilić S, Milošević V, Petrović M, Petrić D, et al. First report of imported case of dengue fever in Republic of Serbia. Vol. 14, Travel Medicine and Infectious Disease. 2016. p. 60–1. doi: 10.1016/j.tmaid.2016.01.001 26822849
47. van der Linden P and Mitchell JFB, editors. ENSEMBLES Climate Change and its Impacts. Summary of research and results from the ENSEMBLES project. Exeter: Met Office Hadley Centre;. 2009.
48. Ministry of Environmental Protection of the Republic of Serbia. Second National Communication of the Republic of Serbia under the United Nations Framework Convention on Climate Change [Internet]. Belgrade: The Ministry; 2017 Aug. 162p [cited 2019 Sep 10]. Available from: http://www.klimatskepromene.rs/wp-content/uploads/2017/09/SNC_eng.pdf
49. Hawkins E, Sutton R. The Potential to Narrow Uncertainty in Regional Climate Predictions. Bull Amer Meteor. 2009; 90(8): 1095–1108. https://doi.org/10.1175/2009bams2607.1
50. Adamović Ž. Distribution and abundance of anopheline mosquitoes (Diptera, Culicidae) in Vojvodina, Serbia. Acta Entomol Jugosl. 1979;15: 82.
51. Halgoš J, Benková I. First record of Anopheles hyrcanus (Diptera: Culicidae) from Slovakia. Biológia. 2004;59: 68.
52. Votýpka J, Šeblová V, Rádrová J. Spread of the West Nile virus vector Culex modestus and the potential malaria vector Anopheles hyrcanus in central Europe. J Vector Ecol. 2008;33: 269–277. doi: 10.3376/1081-1710-33.2.269 19263846
53. Lebl K, Brugger K, Rubel F. Predicting Culex pipiens/restuans population dynamics by interval lagged weather data. Parasit Vectors. 2013;6: 129. doi: 10.1186/1756-3305-6-129 23634763
54. Novikov YM, Vaulin OV. Expansion of Anopheles maculipennis s.s. (Diptera: Culicidae) to northeastern Europe and northwestern Asia: Causes and Consequences. Parasit Vectors. 2014;7: 389. doi: 10.1186/1756-3305-7-389 25148893
55. Seidel B, Nowotny N, Duh D, Indra A, Hufnagl P, Allerberger F. First records of the thermophilic mosquito Culiseta longiareolata (Macquart, 1838) in Austria, 2012, and in Slovenia, 2013. J Eur Mosq Control Assoc. 2013;31: 17–20.
56. Zittra C, Waringer J, Werblow A, Melaun C, Fuehrer H-P. Reconfirmation of Culiseta (Allotheobaldia) longiareolata (Macquart 1838) (Diptera: Culicidae) in Austria. The first sequence-confirmed findings in northeastern Austria. Acta ZooBot Austria. 2014;150: 17–24.
57. Tippelt L, Walther D, Kampen H. The thermophilic mosquito species Uranotaenia unguiculata Edwards, 1913 (Diptera: Culicidae) moves north in Germany. Parasitol Res. 2017;116: 3437–3440. doi: 10.1007/s00436-017-5652-2 29103095
58. Riordan DF. Effects of ultraviolet radiation on adults of Aedes aegypti (L.) (Diptera: Culicidae). Mosq News. 1969;29. Available: https://www.cabdirect.org/cabdirect/abstract/19711000053
59. Hori M, Shibuya K, Sato M, Saito Y. Lethal effects of short-wavelength visible light on insects. Sci Rep. 2014;4: 7383. doi: 10.1038/srep07383 25488603
60. Kramer L, Li J, Shi PY. West Nile virus. Lancet Neurol. 2007;6:171–181. doi: 10.1016/S1474-4422(07)70030-3 17239804
61. Blitvich BJ. Transmission dynamics and changing epidemiology of West Nile virus. Anim Health Res Rev. 2008;9:71–86. doi: 10.1017/S1466252307001430 18348742
62. Komar N. West Nile virus: epidemiology and ecology in North America. Adv Vir Res. 2003;6:185–234.
63. Calistri P, Giovannini A, Hubalek Z, Ionescu A, Monaco F, Savini G et al. Epidemiology of West Nile in Europe and in the Mediterranean basin. Open Virol J. 2010;4(1):29–37.
64. Medić S, Lazić S, Petrović T, Petrić D, Samojlović M, Lazić G et al. Evidence of the first clinical case of equine neuroinvasive West Nile Disease in Serbia, 2018. Acta veterinaria 2019;69(1):123–130.
65. Paz S, Malkinson D, Green MS, Tsioni G, Papa A, Danis K, et al. Permissive Summer Temperatures of the 2010 European West Nile Fever Upsurge. PLOS ONE. 2013;8: e56398. doi: 10.1371/journal.pone.0056398 23431374
66. Marcantonio M, Rizzoli A, Metz M, Rosà R, Marini G, Chadwick E, et al. Identifying the Environmental Conditions Favouring West Nile Virus Outbreaks in Europe. Gourbiere S, editor. PLOS ONE. 2015;10: e0121158. doi: 10.1371/journal.pone.0121158 25803814
67. Tran A, Sudre B, Paz S, Rossi M, Desbrosse A, Chevalier V, et al. Environmental predictors of West Nile fever risk in Europe. Int J Health Geogr. 2014;13: 26. doi: 10.1186/1476-072X-13-26 24986363
68. Paz S. Climate change impacts on West Nile virus transmission in a global context. Philos Trans R Soc B Biol Sci. 2015;370: 20130561–20130561. doi: 10.1098/rstb.2013.0561 25688020
69. Narayanan DL, Saladi RN, Fox JL. Ultraviolet radiation and skin cancer. Int J Dermatol. 2010;49: 978–986. doi: 10.1111/j.1365-4632.2010.04474.x 20883261
70. Lacy K, Alwan W. Skin cancer. Medicine (Baltimore). 2013;41: 402–405. doi: 10.1016/j.mpmed.2013.04.008
71. Leun JC van der. Uv-Carcinogenesis. Photochem Photobiol. 1984;39: 861–868. doi: 10.1111/j.1751-1097.1984.tb08872.x 6379697
72. Malinovic-Milicevic S, Mihailovic DT, Radovanovic MM. Reconstruction of the erythemal UV radiation data in Novi Sad (Serbia) using the NEOPLANTA parametric model. Theor Appl Climatol. 2015;121: 131–138. doi: 10.1007/s00704-014-1223-y
73. Malinović-Milićević S and Radovanović M 2016 UV radiation and heat wives in Vojvodina (Book 87) (Belgrade: Serbian Academy of Sciences and Arts, Institute of Geography „Jovan Cvijić”) p ISBN.
74. de Vries E, Tyczynski J, Parkin MD. Cutaneous malignant melanoma in Europe. Eur Netw Cancer Regist Int Agency Res Cancer. 2003;ENCR Cancer Fact Sheets no.4 1–4.
75. Demers AA, Nugent Z, Mihalcioiu C, Wiseman MC, Kliewer EV. Trends of nonmelanoma skin cancer from 1960 through 2000 in a Canadian population. J Am Acad Dermatol. 2005;53: 320–328. doi: 10.1016/j.jaad.2005.03.043 16021129
76. Wu S, Han J, Laden F, Qureshi AA. Long-term ultraviolet flux, other potential risk factors, and skin cancer risk: a cohort study. Cancer epidemiology, biomarkers and prevention: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2014;23(6):1080–1089.
Článek vyšel v časopise
PLOS One
2020 Číslo 1
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Severity of misophonia symptoms is associated with worse cognitive control when exposed to misophonia trigger sounds
- Chemical analysis of snus products from the United States and northern Europe
- Calcium dobesilate reduces VEGF signaling by interfering with heparan sulfate binding site and protects from vascular complications in diabetic mice
- Effect of Lactobacillus acidophilus D2/CSL (CECT 4529) supplementation in drinking water on chicken crop and caeca microbiome