Search

Introduction: The establishment of Aedes albopictus in new areas in Europe has changed the risk of local dengue transmission represented by imported human cases. The risk of transmission is determined by the distribution of travelers arriving from dengue-endemic areas and the distribution of Ae. albopictus as potential vectors of dengue in Spain., Methods: Environmental, entomological, epidemiological, demographic, tourism and travel data were analyzed to produce a series of maps to represent: the distribution of Ae. albopictus across municipalities; the risk of expansion of Ae. albopictus based on a species distribution model; the calculated index of travelers from dengue-endemic areas (IDVZE) per province; the percentage contribution of each municipality to the total number of cases in Spain. The maps were then added using map algebra, to profile the spatial risk of autochthonous dengue in Spain at a municipal level from 2016 to 2018., Results: Ae. albopictus was detected in 983 municipalities. The calculated IDVZE varied from 0.23 to 10.38, with the highest IDVZE observed in Madrid. The overall risk of autochthonous cases oscillated between 0.234 and 115, with the very high risk and high risk areas detected in the Mediterranean region, mainly in the Levantine coast and some parts of the Balearic Islands. Most of the interior of the peninsula was characterized as low risk., Conclusion: Prevention and control measures to mitigate the risk of autochthonous dengue should be prioritized for municipalities in the high risk areas integrating early detection of imported dengue cases and vector control., (Copyright © 2023 The Author(s). Published by Elsevier España, S.L.U. All rights reserved.)

Malaria remains one of the most important infectious diseases globally due to its high incidence and mortality rates. The influx of infected cases from endemic to non-endemic malaria regions like Europe has resulted in a public health concern over sporadic local outbreaks. This is facilitated by the continued presence of competent Anopheles vectors in non-endemic countries.We modelled the potential distribution of the main malaria vector across Spain using the ensemble of eight modelling techniques based on environmental parameters and the Anopheles maculipennis s.l. presence/absence data collected from 2000 to 2020. We then combined this map with the number of imported malaria cases in each municipality to detect the geographic hot spots with a higher risk of local malaria transmission.The malaria vector occurred preferentially in irrigated lands characterized by warm climate conditions and moderate annual precipitation. Some areas surrounding irrigated lands in northern Spain (e.g. Zaragoza, Logroño), mainland areas (e.g. Madrid, Toledo) and in the South (e.g. Huelva), presented a significant likelihood of A. maculipennis s.l. occurrence, with a large overlap with the presence of imported cases of malaria.While the risk of malaria re-emergence in Spain is low, it is not evenly distributed throughout the country. The four recorded local cases of mosquito-borne transmission occurred in areas with a high overlap of imported cases and mosquito presence. Integrating mosquito distribution with human incidence cases provides an effective tool for the quantification of large-scale geographic variation in transmission risk and pinpointing priority areas for targeted surveillance and prevention.

BackgroundQ fever is a bacterial zoonosis caused by Coxiella burnetii. Spain has the highest number of notified human cases in Europe. Small ruminants are a key reservoir for the pathogen, transmission from animals to humans is usually airborne.AimWe aimed at exploring temporal and spatial epidemiological patterns of sporadic and outbreak cases of Q fever in four Spanish regions with the highest number of notified cases.MethodsWe extracted data on Q fever cases in the Canary Islands, Basque Country, La Rioja and Navarre between 2016 and 2022 from the Spanish National Epidemiological Surveillance Network. We calculated standardised incidence ratios (SIR), spatial relative risks (sRR) and posterior probabilities (PP) utilising Besag-York-Mollié models.ResultsThere were 1,059 notifications, with a predominance of males aged 30-60 years. In Basque Country, La Rioja and Navarre area, 11 outbreaks were reported, while no in the Canary Islands. A seasonal increase in incidence rates was observed between March and June. In the Canary Islands, elevated sRR was seen in La Palma, Gran Canaria, Lanzarote and Fuerteventura. In Basque Country, La Rioja and Navarre area, the highest sRR was identified in the south of Biscay province.ConclusionGoats were the main source for humans in outbreaks reported in the literature. Seasonal increase may be related to the parturition season of small ruminants and specific environmental conditions. Local variations in sRR within these regions likely result from diverse environmental factors. Future One Health-oriented studies are essential to deepen our understanding of Q fever epidemiology.

Background: Understanding the patterns of disease importation through international travel is paramount for effective public health interventions and global disease surveillance. While global airline network data have been used to assist in outbreak prevention and effective preparedness, accurately estimating how these imported cases disseminate locally in receiving countries remains a challenge., Objective: This study aimed to describe and understand the regional distribution of imported cases of dengue and malaria upon arrival in Spain via air travel., Methods: We have proposed a method to describe the regional distribution of imported cases of dengue and malaria based on the computation of the "travelers' index" from readily available socioeconomic data. We combined indicators representing the main drivers for international travel, including tourism, economy, and visits to friends and relatives, to measure the relative appeal of each region in the importing country for travelers. We validated the resulting estimates by comparing them with the reported cases of malaria and dengue in Spain from 2015 to 2019. We also assessed which motivation provided more accurate estimates for imported cases of both diseases., Results: The estimates provided by the best fitted model showed high correlation with notified cases of malaria (0.94) and dengue (0.87), with economic motivation being the most relevant for imported cases of malaria and visits to friends and relatives being the most relevant for imported cases of dengue., Conclusions: Factual descriptions of the local movement of international travelers may substantially enhance the design of cost-effective prevention policies and control strategies, and essentially contribute to decision-support systems. Our approach contributes in this direction by providing a reliable estimate of the number of imported cases of nonendemic diseases, which could be generalized to other applications. Realistic risk assessments will be obtained by combining this regional predictor with the observed local distribution of vectors., (©David García-García, Beatriz Fernández-Martínez, Frederic Bartumeus, Diana Gómez-Barroso. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 27.05.2024.)

Tuberculosis (TB) remains a global leading cause of death, necessitating an investigation into its unequal distribution. Sun exposure, linked to vitamin D (VD) synthesis, has been proposed as a protective factor. This study aimed to analyse TB rates in Spain over time and space and explore their relationship with sunlight exposure. An ecological study examined the associations between rainfall, sunshine hours, and TB incidence in Spain. Data from the National Epidemiological Surveillance Network (RENAVE in Spanish) and the Spanish Meteorological Agency (AEMET in Spanish) from 2012 to 2020 were utilized. Correlation and spatial regression analyses were conducted. Between 2012 and 2020, 43,419 non-imported TB cases were reported. A geographic pattern (north-south) and distinct seasonality (spring peaks and autumn troughs) were observed. Sunshine hours and rainfall displayed a strong negative correlation. Spatial regression and seasonal models identified a negative correlation between TB incidence and sunshine hours, with a four-month lag. A clear spatiotemporal association between TB incidence and sunshine hours emerged in Spain from 2012 to 2020. VD levels likely mediate this relationship, being influenced by sunlight exposure and TB development. Further research is warranted to elucidate the causal pathway and inform public health strategies for improved TB control.

Introduction: The main preventive measure against invasive meningococcal disease is vaccination. The aim of our study was to evaluate the acceptability of the meningococcal B (MenB) vaccine and socioeconomic inequalities in the access to the vaccine in the Community of Madrid in the period prior to its introduction in the immunization schedule., Materials and Methods: We conducted an observational and ecological descriptive study in the cohort of children born between 2016 and 2019 using population-based electronic records. We calculated the vaccination coverage and analysed factors associated with vaccination status, determined the spatial distribution of vaccination coverage and the deprivation index (DI) and assessed the association between them by means of spatial regression., Results: We observed an increasing trend in primary vaccination coverage, from 44% in the cohort born in 2016 to 68% in the 2019 cohort. We found a statistically significant association between vaccination status and the DI (OR of primary vaccination in areas with DI5 compared to areas with DP1, 0.38; 95% confidence interval, 0.39-0.50; P<.001). The spatial analysis showed an inverse correlation between the DI and vaccination coverage., Conclusions: The rise in the coverages of the MenB vaccine shows acceptance by the population. The association between socioeconomic level and vaccination coverage confirms the existence of health inequality and underlines the importance including this vaccine in the immunization schedule., (Copyright © 2023. Published by Elsevier España, S.L.U.)

Background: Social restrictions and vaccination seem to have shaped the pandemic development in Europe, but the influence of geographical position is still debated. This study aims to verify whether the pandemic spread through Europe following a particular direction, during the period between the start of the pandemic and November 2021. The existence of a spatial gradient for epidemic intensity is also hypothesized., Methods: Daily COVID-19 epidemiological data were extracted from Our World in Data COVID-19 database, which also included vaccination and non-pharmacological interventions data. Latitude and longitude of each country's centroid were used as geographic variables. Epidemic periods were delimited from epidemic surge data. Multivariable linear and Cox's regression models were performed for each epidemic period to test if geographical variables influenced surge dates. Generalized additive models (GAM) were used to test the spatial gradient hypothesis with three epidemic intensity measures., Results: Linear models suggest a possible west-east shift in the first epidemic period and features a significant association of NPIs with epidemic surge delay. Neither latitude nor longitude had significant associations with epidemic surge timing in both second and third periods. Latitude displays strong negative associations with all epidemic intensity measures in GAM models. Vaccination was also negatively associated with intensity., Conclusions: A longitudinal spread of the pandemic in Europe seems plausible, particularly concerning the first wave. However, a recurrent trend was not observed. Southern Europe countries may have experienced increased transmissibility and incidence, despite climatic conditions apparently unfavourable to the virus., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Martínez-Portillo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In Spain, the largest human West Nile virus (WNV) outbreak among humans was reported in 2020, constituting the second most important outbreak in Europe that season. Extremadura (southwestern Spain) was one of the affected areas, reporting six human cases. The first autochthonous human case in Spain was reported in Extremadura in 2004, and no other human cases were reported until 2020. In this work, we describe the first WNV human outbreak registered in Extremadura, focusing on the most important clinical aspects, diagnostic results, and control actions which followed. In 2020, from September to October, human WNV infections were diagnosed using a combination of molecular and serological methods (an in-house specific qRT-PCR and a commercial ELISA for anti-WNV IgM and IgG antibodies) and by analysing serum, urine, and/or cerebrospinal fluid samples. Serological positive serum samples were further tested using commercial kits against related flaviviruses Usutu and Tick-borne encephalitis in order to analyse serological reactivity and to confirm the results by neutralisation assays. In total, six cases of WNV infection (five with neuroinvasive disease and one with fever) were identified. Clinical presentation and laboratory findings are described. No viral RNA was detected in any of the analysed samples, but serological cross-reactivity was detected against the other tested flaviviruses. Molecular and serological methods for WNV detection in various samples as well as differential diagnosis are recommended. The largest number of human cases of WNV infection ever registered in Extremadura, Spain, occurred in 2020 in areas where circulation of WNV and other flaviviruses has been previously reported in humans and animals. Therefore, it is necessary to enhance surveillance not only for the early detection and implementation of response measures for WNV but also for other emerging flaviviruses that could be endemic in this area., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Macias, Martín, Pérez-Olmeda, Fernández-Martínez, Gómez-Barroso, Fernández, Ramos, Herrero, Rodríguez, Delgado, Sánchez-Seco, Galán, Corbacho, Jimenez, Montero-Peña, Valle and Vázquez.)

We analysed the transmission of the human mpox virus in Spain by estimating the effective reproduction number of the disease from official surveillance data. Our computations show that this decreased steadily after an initial burst phase, dropping below 1 on July 12, and thus the outbreak was expected to reduce in the following weeks. Differences in trends were found across geographical regions of the country and across MSM and heterosexual populations.

Human mobility drives the geographical diffusion of infectious diseases at different scales, but few studies focus on mobility itself. Using publicly available data from Spain, we define a Mobility Matrix that captures constant flows between provinces by using a distance-like measure of effective distance to build a network model with the 52 provinces and 135 relevant edges. Madrid, Valladolid and Araba/Álaba are the most relevant nodes in terms of degree and strength. The shortest routes (most likely path between two points) between all provinces are calculated. A total of 7 mobility communities were found with a modularity of 63%, and a relationship was established with a cumulative incidence of COVID-19 in 14 days (CI14) during the study period. In conclusion, mobility patterns in Spain are governed by a small number of high-flow connections that remain constant in time and seem unaffected by seasonality or restrictions. Most of the travels happen within communities that do not completely represent political borders, and a wave-like spreading pattern with occasional long-distance jumps (small-world properties) can be identified. This information can be incorporated into preparedness and response plans targeting locations that are at risk of contagion preventively, underscoring the importance of coordination between administrations when addressing health emergencies.

Objectives: To study the spatio-temporal distribution of cases of invasive pneumococcal disease (IPD) due to serotypes resistant to erythromycin and its relationship with community consumption of macrolides and childhood vaccination coverage., Methods: We selected IPD cases in adults over 59 years old, residents in the Community of Madrid (MC), notified in the period 2007-2016. The variables studied were obtained from the Vaccination Information Systems and the Pharmaceutical Service. The cut-off point (minimum inhibitory erythromycin concentration > 0.5 mg/L) of the EUCAST classification was used to define erythromycin resistant serotypes. We used JointPoint to estimate the incidence trends by erythromycin resistant serotypes included in the 13-valent vaccine (STPCV13) and not included in it (STnoPCV13). The association of these incidences with the community consumption of macrolides and vaccination coverage was made using Poisson models. Statistical scanning was used for the detection of temporal-spaces clusters of cases., Results: 1936 cases were identified, of which 427 erythromycin resistant serotypes were identified. The incidence of all cases due to resistant serotypes was decreasing (AAPC: -5,40%). During the period studied, the incidence of cases due to erythromycin resistant STPCV13 was decreasing with an annual percentage change (APC): -13.8 and was inversely associated with childhood vaccination coverage (IRR 0.641), while that of cases due to erythromycin resistant STnoPCV13 was ascending (APC): 4.5; and was not associated with coverage. 1 cluster was detected by STnoPCV13 and none by STPCV13 after the date of inclusion of the 13-valent in the childhood vaccination calendar., Conclusions: The decrease in IPD due to resistant STPCV13 was associated with an increase in childhood vaccination coverage. The presence of clusters due to STnoPCV13 after the date of inclusion of the 13-valent vaccine in the childhood vaccination calendar indicates serotypes replacement. The increase in cases of resistant STnoPCV13 could be related to the replacement of vaccine serotypes in nasopharyngeal colonization, facilitated by the consumption of macrolides still at high levels in MC., (Copyright © 2021 Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica. Published by Elsevier España, S.L.U. All rights reserved.)

Background: Several environmental factors seem to be involved in childhood leukaemia incidence. Traffic exposure could increase the risk while urban green spaces (UGS) exposure could reduce it. However, there is no evidence how these two factors interact on this infant pathology., Objectives: to evaluate how residential proximity to UGS could be an environmental protective factor against traffic exposure on childhood leukaemia incidence., Methods: A population-based case control study was conducted across thirty Spanish regions during the period 2000-2018. It included 2526 incident cases and 15,156, individually matched by sex, year-of-birth, and place-of-residence. Using the geographical coordinates of the participants' home residences, a 500 m proxy for exposure to UGS was built. Annual average daily traffic (AADT) was estimated for all types of roads 100 m near the children's residence. Odds ratios (ORs) and 95% confidence intervals (95% CIs), UGS, traffic exposure, and their possible interactions were calculated for overall childhood leukaemia, and the acute lymphoblastic (ALL) and acute myeloblastic leukaemia (AML) subtypes, with adjustment for socio-demographic covariates., Results: We found an increment of childhood leukaemia incidence related to traffic exposure, for every 100 AADT increase the incidence raised 1.1% (95% CI: 0.58-1.61%). UGS exposure showed an incidence reduction for the highest exposure level, Q5 (OR = 0.63; 95% CI = 0.54-0.72). Regression models with both traffic exposure and UGS exposure variables showed similar results but the interaction was not significant., Conclusions: Despite their opposite effects on childhood leukaemia incidence individually, our results do not suggest a possible interaction between both exposures. This is the first study about the interaction of these two environmental factors; consequently, it is necessary to continue taking into account more individualized data and other possible environmental risk factors involved.

Introduction/objective: West Nile virus (WNV) is one of the most widely distributed flaviviruses worldwide. It is considered an endemic and emerging pathogen in different areas of the Europe and Mediterranean countries (MR). Mosquitoes of the genus Culex spp. are the main vectors, and birds its main vertebrate hosts. It can occasionally infect mammals, including humans. Different environmental factors can influence its distribution and transmission through its effects on vector or host populations. Our objective was to determine environmental factors associated with changes in vector distribution and WNV transmission in Europe and MR., Material & Methods: Systematic peer review of articles published between 2000 and 2020. We selected studies on WNV, and its vectors carried out in Europe and MR. The search included terms referring to climatic and environmental factors., Results: We included 65 studies, of which 21 (32%) were conducted in Italy. Culex spp. was studied in 26 papers (40%), humans in 19 papers (29%) and host animals (mainly horses) in 16 papers (25%), whereas bird reservoirs were addressed in 5 studies (8%). A significant positive relationship was observed between changes in temperature and precipitation patterns and the epidemiology of WNV, although contrasting results were found among studies. Other factors positively related to WNV dynamics were the normalized difference vegetation index (NDVI] and expansion of anthropized habitats., Conclusion: The epidemiology of WNV seems to be related to climatic factors that are changing globally due to ongoing climate change. Unfortunately, the complete zoonotic cycle was not analyzed in most papers, making it difficult to determine the independent impact of environment on the different components of the transmission cycle. Given the current expansion and endemicity of WNV in the area, it is important to adopt holistic approaches to understand WNV epidemiology and to improve WNV surveillance and control., Competing Interests: None., (© 2022 Published by Elsevier B.V.)

The aim of our study was to describe the results of the epidemiological surveillance of hepatitis A infections in Spain in the context of the 2016/2017 European outbreak, particularly of hepatitis A outbreaks reported in the MSM population, incorporating the results of a spatio-temporal analysis of cases. Hepatitis A cases and outbreaks reported in 2016-2017 to the National Epidemiological Surveillance Network were reviewed: outbreaks in which some of the cases belonged to the MSM group were described, and clusters of hepatitis A cases in men and women were analysed using a space-time scan statistic. Twenty-six outbreaks were identified, with a median size of two cases per outbreak, with most of the outbreak-related cases belonging to the 15-44 years-old group. Nearly 85% occurred in a household setting, and in all outbreaks, the mode of transmission was direct person-to-person contact. Regarding space-time analysis, twenty statistically significant clusters were identified in the male population and eight in the female population; clusters in men presented a higher number of observed cases and affected municipalities, as well as a higher percentage of municipalities classified as large urban areas. The elevated number of cases detected in clusters of men indicates that the number of MSM-related outbreaks may be higher than reported, showing that spatio-temporal analysis is a complementary, useful tool which may improve the detection of outbreaks in settings where epidemiological investigation may be more challenging.

Background: To analyse differences in confirmed cases, hospitalisations and deaths due to COVID-19 related to census section socioeconomic variables.  METHODS: Ecological study in the 12 largest municipalities in Andalusia (Spain) during the first three epidemic waves of the COVID-19 (02/26/20-03/31/21), covering 2,246 census sections (unit of analysis) and 3,027,000 inhabitants. Incidence was calculated, standardised by age and sex, for infection, hospitalisation and deaths based on average gross income per household (AGI) for the census tracts in each urban area. Association studied using a Poisson Bayesian regression model with random effects for spatial smoothing., Results: There were 140,743 cases of COVID-19, of which 12,585 were hospitalised and 2,255 died. 95.2% of cases were attributed to the second and third waves, which were jointly analysed. We observed a protective effect of income for infection in 3/12 cities. Almeria had the largest protective effect (smoothed relative risk (SRR) = 0.84 (0.75-0.94 CI 95%). This relationship reappeared with greater magnitude in 10/12 cities for hospitalisation, lowest risk in Algeciras SRR = 0.41 (0.29-0.56). The pattern was repeated for deaths in all urban areas and reached statistical significance in 8 cities. Lowest risk in Dos Hermanas SRR = 0.35 (0.15-0.81)., Conclusions: Income inequalities by geographical area were found in the incidence of COVID-19. The strengths of the association increased when analysing the severe outcomes of hospitalisations and, above all, deaths., (© 2022. The Author(s).)

A number of studies suggest that meteorological conditions are related to the risk of Legionnaires' disease (LD) but the findings are not consistent. A systematic review was conducted to investigate the association of weather with sporadic LD and highlight the key meteorological conditions related to this outcome. PubMed, EMBASE, The Cochrane Library and OpenGrey were searched on 26-27 March 2020 without date, language or location restrictions. Key words included "legionellosis", "legionnaires' disease", combined with "meteorological conditions", "weather", "temperature", "humidity", "rain", "ultraviolet rays", "wind speed", etc. Studies were excluded if they did not examine the exposure of interest, the outcome of interest and their association or if they only reported LD outbreak cases. The study was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines and it was registered in PROSPERO (#CRD42020168869). There were 811 articles, of which 17 were included in the review. The studies investigated different meteorological variables and most of them examined the combined effect of several variables. The most commonly examined factors were precipitation and temperature, followed by relative humidity. The studies suggested that increased precipitation, temperature and relative humidity were positively associated with the incidence of LD. There was limited evidence that higher wind speed, pressure, visibility, UV radiation and longer sunshine duration were inversely linked with the occurrence of LD. A period of increased but not very high temperatures, followed by a period of increased precipitation, favour the occurrence of LD. Increased awareness of the association of temperature and precipitation and LD occurrence among clinicians and public health professionals can improve differential diagnosis for cases of sporadic community-acquired pneumonia and at the same time contribute to improving LD surveillance., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Background: There is scarce information focused on the effect of weather conditions and air pollution on specific acute viral respiratory infections, such as rhinovirus (RV), with a wide clinical spectrum of severity., Objective: The aim of this study was to analyze the association between episodes of severe respiratory tract infection by RV and air pollutant concentrations (NO x and SO 2 ) in the reference area of a pediatric university hospital., Methods: An analysis of temporal series of daily values of NO x and SO 2 , weather variables, circulating pollen and mold spores, and daily number of admissions in the pediatric intensive care unit (PICU) with severe respiratory RV infection (RVi) in children between 6 months and 18 years was performed. Lagged variables for 0-5 days were considered. The study spanned from 2010 to 2018. Patients with comorbidities were excluded., Results: One hundred and fifty patients were admitted to the PICU. Median age was 19 months old (interquartile range [IQR]: 11-47). No relationship between RV-PICU admissions and temperature, relative humidity, cumulative rainfall, or wind speed was found. Several logistic regression models with one pollutant and two pollutants were constructed but the best model was that which included average daily NO x concentrations. Average daily NO x concentrations were related with the presence of PICU admissions 3 days later (odds ratio per IQR-unit increase: 1.64, 95% confidence interval: 1.20-2.25))., Conclusions: This study has shown a positive correlation between NO x concentrations at Lag 3 and children's PICU admissions with severe RV respiratory infection. Air pollutant data should be taken into consideration when we try to understand the severity of RVis., (© 2022 The Authors. Pediatric Pulmonology published by Wiley Periodicals LLC.)

BackgroundAfter a national lockdown during the first wave of the COVID-19 pandemic in Spain, regional governments implemented different non-pharmaceutical interventions (NPIs) during the second wave.AimTo analyse which implemented NPIs significantly impacted effective reproduction number (R t ) in seven Spanish provinces during 30 August 2020-31 January 2021.MethodsWe coded each NPI and levels of stringency with a 'severity index' (SI) and computed a global SI (mean of SIs per six included interventions). We performed a Bayesian change point analysis on the R t curve of each province to identify possible associations with global SI variations. We fitted and compared several generalised additive models using multimodel inference, to quantify the statistical effect on R t of the global SI (stringency) and the individual SIs (separate effect of NPIs).ResultsThe global SI had a significant lowering effect on the R t (mean: 0.16 ± 0.05 units for full stringency). Mandatory closing times for non-essential businesses, limited gatherings, and restricted outdoors seating capacities (negative) as well as curfews (positive) were the only NPIs with a significant effect. Regional mobility restrictions and limited indoors seating capacity showed no effect. Our results were consistent with a 1- to 3-week-delayed R t as a response variable.ConclusionWhile response measures implemented during the second COVID-19 wave contributed substantially to a decreased reproduction number, the effectiveness of measures varied considerably. Our findings should be considered for future interventions, as social and economic consequences could be minimised by considering only measures proven effective.

Background: Blue spaces have been a key part of human evolution, providing resources and helping economies develop. To date, no studies have been carried out to explore how they may be linked to paediatric oncological diseases., Objectives: To explore the possible relationship of residential proximity to natural and urban blue spaces on childhood leukaemia., Methods: A population-based case-control study was conducted in four regions of Spain across the period 2000-2018. A total of 936 incident cases and 5616 controls were included, individually matched by sex, year of birth and place of residence. An exposure proxy with four distances (250 m, 500 m, 750 m, and 1 km) to blue spaces was built using the geographical coordinates of the participants' home residences. Odds ratios (ORs) and 95% confidence intervals (95%CIs) for blue-space exposure were calculated for overall childhood leukaemia, and the acute lymphoblastic (ALL) and acute myeloblastic leukaemia (AML) subtypes, with adjustment for socio-demographic and environmental covariates., Results: A decrease in overall childhood leukaemia and ALL-subtype incidence was found as we came nearer to children's places of residence, showing, for the study as a whole, a reduced incidence at 250 m (odds ratio (OR) = 0.77; 95%CI = 0.60-0.97), 500 m (OR = 0.78; 95%CI = 0.65-0.93), 750 m (OR = 0.80; 95%CI = 0.69-0.93), and 1000 m (OR = 0.84; 95%CI = 0.72-0.97). AML model results showed an increasing incidence at closest to subjects' homes (OR at 250m = 1.06; 95%CI=0.63-1.71)., Conclusions: Our results suggest a possible association between lower childhood leukaemia incidence and blue-space proximity. This study is a first approach to blue spaces' possible effects on childhood leukaemia incidence; consequently, it is necessary to continue studying these spaces-while taking into account more individualised data and other possible environmental risk factors.

When responding to infectious disease outbreaks, rapid and accurate estimation of the epidemic trajectory is critical. However, two common data collection problems affect the reliability of the epidemiological data in real time: missing information on the time of first symptoms, and retrospective revision of historical information, including right censoring. Here, we propose an approach to construct epidemic curves in near real time that addresses these two challenges by 1) imputation of dates of symptom onset for reported cases using a dynamically-estimated "backward" reporting delay conditional distribution, and 2) adjustment for right censoring using the NobBS software package to nowcast cases by date of symptom onset. This process allows us to obtain an approximation of the time-varying reproduction number (Rt) in real time. We apply this approach to characterize the early SARS-CoV-2 outbreak in two Spanish regions between March and April 2020. We evaluate how these real-time estimates compare with more complete epidemiological data that became available later. We explore the impact of the different assumptions on the estimates, and compare our estimates with those obtained from commonly used surveillance approaches. Our framework can help improve accuracy, quantify uncertainty, and evaluate frequently unstated assumptions when recovering the epidemic curves from limited data obtained from public health systems in other locations., Competing Interests: I have read the journal’s policy and some of the co-authors of this manuscript have the following competing interests: ML discloses honoraria/consulting from Merck, Affinivax, Sanofi-Pasteur, Bristol Myers-Squibb, and Antigen Discovery; research funding (institutional) from Pfizer, and an unpaid scientific advice to Janssen, Astra-Zeneca, One Day Sooner, and Covaxx (United Biomedical). MS discloses having received institutional research support from Johnson and Johnson.The rest of co-authors declare no competing interest.

Introduction: Scabies is a neglected disease stablished worldwide with a fairy well determined incidence. In high-income countries, it often causes outbreaks affecting the residents and staff of institutions and long-term facilities, usually hard to detect and control due to the difficult diagnosis and notification delay. This study aim at characterizing the affected population, geographical distribution, and evolution of scabies in Spain from 1997-2019 as well as to describe the main environments of transmission using different data sources., Methods: We carried out a nationwide retrospective study using four databases, which record data from different perspectives: hospital admissions, patients attended at primary healthcare services, outbreaks, and occupational diseases. We described the main characteristics from each database and calculated annual incidences in order to evaluate temporal and geographical patterns. We also analyzed outbreaks and occupational settings to characterize the main transmission foci and applied Joinpoint regression models to detect trend changes., Results: The elderly was the most frequent collective among the hospital admitted patients and notified cases in outbreaks, while children and young adults were the most affected according to primary care databases. The majority of the outbreaks occurred in homes and nursing homes; however, the facilities with more cases per outbreak were military barracks, healthcare settings and nursing homes. Most occupational cases occurred also in healthcare and social services settings, being healthcare workers the most common affected professional group. We detected a decreasing trend in scabies admissions from 1997 to 2014 (annual percentage change -APC- = -11.2%) and an increasing trend from 2014 to 2017 (APC = 23.6%). Wide geographical differences were observed depending on the database explored., Discussion: An increasing trend in scabies admissions was observed in Spain since 2014, probably due to cutbacks in social services and healthcare in addition to worsen of living conditions as a result of the 2008 economic crisis, among other reasons. The main transmission foci were healthcare and social settings. Measures including enhancing epidemic studies and national registries, reinforcing clinical diagnosis and early detection of cases, hygiene improvements and training of the staff and wide implementation of scabies treatment (considering mass drug administration in institutions outbreaks) should be considered to reduce the impact of scabies among most vulnerable groups in Spain., Competing Interests: The authors have declared that no competing interests exist.

Objectives: To study the spatio-temporal distribution of cases of invasive pneumococcal disease (IPD) due to serotypes resistant to erythromycin and its relationship with community consumption of macrolides and childhood vaccination coverage., Methods: We selected IPD cases in adults over 59 years old, residents in the Community of Madrid (MC), notified in the period 2007-2016. The variables studied were obtained from the Vaccination Information Systems and the Pharmaceutical Service. The cut-off point (minimum inhibitory erythromycin concentration > 0.5 mg/L) of the EUCAST classification was used to define erythromycin resistant serotypes. We used JointPoint to estimate the incidence trends by erythromycin resistant serotypes included in the 13-valent vaccine (STPCV13) and not included in it (STnoPCV13). The association of these incidences with the community consumption of macrolides and vaccination coverage was made using Poisson models. Statistical scanning was used for the detection of temporal-spaces clusters of cases., Results: 1936 cases were identified, of which 427 erythromycin resistant serotypes were identified. The incidence of all cases due to resistant serotypes was decreasing (AAPC: -5,40%). During the period studied, the incidence of cases due to erythromycin resistant STPCV13 was decreasing with an annual percentage change (APC): -13.8 and was inversely associated with childhood vaccination coverage (IRR 0.641), while that of cases due to erythromycin resistant STnoPCV13 was ascending (APC): 4.5; and was not associated with coverage. 1 cluster was detected by STnoPCV13 and none by STPCV13 after the date of inclusion of the 13-valent in the childhood vaccination calendar., Conclusions: The decrease in IPD due to resistant STPCV13 was associated with an increase in childhood vaccination coverage. The presence of clusters due to STnoPCV13 after the date of inclusion of the 13-valent vaccine in the childhood vaccination calendar indicates serotypes replacement. The increase in cases of resistant STnoPCV13 could be related to the replacement of vaccine serotypes in nasopharyngeal colonization, facilitated by the consumption of macrolides still at high levels in MC., (Copyright © 2021 Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica. Publicado por Elsevier España, S.L.U. All rights reserved.)

Background: On June 21st de-escalation measures and state-of-alarm ended in Spain after the COVID-19 first wave. New surveillance and control strategy was set up to detect emerging outbreaks., Aim: To detect and describe the evolution of COVID-19 clusters and cases during the 2020 summer in Spain., Methods: A near-real time surveillance system to detect active clusters of COVID-19 was developed based on Kulldorf's prospective space-time scan statistic (STSS) to detect daily emerging active clusters., Results: Analyses were performed daily during the summer 2020 (June 21st - August 31st) in Spain, showing an increase of active clusters and municipalities affected. Spread happened in the study period from a few, low-cases, regional-located clusters in June to a nationwide distribution of bigger clusters encompassing a higher average number of municipalities and total cases by end-August., Conclusion: STSS-based surveillance of COVID-19 can be of utility in a low-incidence scenario to help tackle emerging outbreaks that could potentially drive a widespread transmission. If that happens, spatial trends and disease distribution can be followed with this method. Finally, cluster aggregation in space and time, as observed in our results, could suggest the occurrence of community transmission.

Mosquito borne diseases are a group of infections that affect humans. Emerging or reemerging diseases are those that (re)occur in regions, groups or hosts that were previously free from these diseases: dengue virus; chikungunya virus; Zika virus; West Nile fever and malaria. In Europe, these infections are mostly imported; however, due to the presence of competent mosquitoes and the number of trips both to and from endemic areas, these pathogens are potentially emergent or re-emergent. Present and future climatic conditions, as well as meteorological, environmental and demographic aspects are risk factors for the distribution of different vectors and/or diseases. This review aimed to identify and analyze the existing literature on the transmission of mosquito borne diseases and those factors potentially affecting their transmission risk of them in six southern European countries with similar environmental conditions: Croatia, France, Greece, Italy, Portugal and Spain. In addition, we would identify those factors potentially affecting the (re)introduction or spread of mosquito vectors. This task has been undertaken with a focus on the environmental and climatic factors, including the effects of climate change. We undertook a systematic review of the vectors, diseases and their associations with climactic and environmental factors in European countries of the Mediterranean region. We followed the PRISMA guidelines and used explicit and systematic methods to identify, select and critically evaluate the studies which were relevant to the topic. We identified 1302 articles in the first search of the databases. Of those, 160 were selected for full-text review. The final data set included 61 articles published between 2000 and 2017.39.3% of the papers were related with dengue, chikungunya and Zika virus or their vectors. Temperature, precipitation and population density were key factors among others. 32.8% studied West Nile virus and its vectors, being temperature, precipitation and NDVI the most frequently used variables. Malaria have been studied in 23% of the articles, with temperature, precipitation and presence of water indexes as the most used variables. The number of publications focused on mosquito borne diseases is increasing in recent years, reflecting the increased interest in that diseases in southern European countries. Climatic and environmental variables are key factors on mosquitoes' distribution and to show the risk of emergence and/or spread of emergent diseases and to study the spatial changes in that distributions., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

During the first months of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) epidemic in 2020, Spain implemented an initial lockdown period on 15 March followed by a strengthened lockdown period on 30 March when only essential workers continued to commute to work. However, little is known about the epidemic dynamics in different age groups during these periods.We used the daily number of coronavirus 2019 cases (by date of symptom onset) reported to the National Epidemiological Surveillance Network among individuals aged 15-19 years through 65-69 years. For each age group g, we computed the proportion PrE(g) of individuals in age group g among all reported cases aged 15-69 years during the pre-lockdown period (1-10 March 2020) and the corresponding proportion PrL(g) during two lockdown periods (initial: 25 March-3 April; strengthened: 8-17 April 2020). For each lockdown period, we computed the proportion ratios PR(g) = PrL(g)/PrE(g). For each pair of age groups g1, g2, PR(g1)>PR(g2) implies a relative increase in the incidence of detected SARS-CoV-2 infection in the age group g1 compared with g2 for the lockdown period vs. the pre-lockdown period.For the initial lockdown period, the highest PR values were in age groups 50-54 years (PR = 1.21; 95% CI: 1.12,1.30) and 55-59 years (PR = 1.19; 1.11,1.27). For the second lockdown period, the highest PR values were in age groups 15-19 years (PR = 1.26; 0.95,1.68) and 50-54 years (PR = 1.20; 1.09,1.31).Our results suggest that different outbreak control measures led to different changes in the relative incidence by age group. During the initial lockdown period, when non-essential work was allowed, individuals aged 40-64 years, particularly those aged 50-59 years, had a higher relative incidence compared with the pre-lockdown period. Younger adults/older adolescents had an increased relative incidence during the later, strengthened lockdown. The role of different age groups during the epidemic should be considered when implementing future mitigation efforts.

Background: The aim of this study is to explore whether transmission of M. leprae has ceased in Spain, based upon the patterns and trends of notified cases., Methodology: Data on new cases reported to the National Leprosy Registry between the years 2003-2018 were extracted. In absence of detailed travel history, cases were considered "autochthonous" or "imported" based on whether they were born within or outside of Spain. These data were analyzed by age, sex, clinical type, country of origin, and location of residence at time of notification., Principal Findings: Data were available on 61 autochthonous and 199 imported cases since 2003. There were clear declines in incidence in both groups, and more imported than autochthonous cases every year since 2006. Autochthonous cases were more frequently multibacillary and had older age at diagnosis compared to imported cases. All the autochthonous cases had been born before 1985 and were more than 25 years old at diagnosis. Male-to-female ratio increased with time for autochthonous cases (except for the last time period). The imported cases originated from 25 countries, half of them from Brasil and Paraguay. Autochthonous cases were mainly distributed in the traditionally endemic regions, especially Andalucía and the eastern Mediterranean coast., Conclusions: Autochthonous and imported cases have different epidemiologic patterns in Spain. There was a clear decline in incidence rates of autochthonous disease, and patterns consistent with those reported from other regions where transmission has ceased. Autochthonous transmission of M. leprae is likely to have now effectively stopped in Spain., Competing Interests: The authors have declared that no competing interests exist.

Changes in environmental conditions, whether related or not to human activities, are continuously modifying the geographic distribution of vectors, which in turn affects the dynamics and distribution of vector-borne infectious diseases. Determining the main ecological drivers of vector distribution and how predicted changes in these drivers may alter their future distributions is therefore of major importance. However, the drivers of vector populations are largely specific to each vector species and region. Here, we identify the most important human-activity-related and bioclimatic predictors affecting the current distribution and habitat suitability of the mosquito Culex pipiens and potential future changes in its distribution in Spain. We determined the niche of occurrence (NOO) of the species, which considers only those areas lying within the range of suitable environmental conditions using presence data. Although almost ubiquitous, the distribution of Cx. pipiens is mostly explained by elevation and the degree of urbanization but also, to a lesser extent, by mean temperatures during the wettest season and temperature seasonality. The combination of these predictors highlights the existence of a heterogeneous pattern of habitat suitability, with most suitable areas located in the southern and northeastern coastal areas of Spain, and unsuitable areas located at higher altitude and in colder regions. Future climatic predictions indicate a net decrease in distribution of up to 29.55%, probably due to warming and greater temperature oscillations. Despite these predicted changes in vector distribution, their effects on the incidence of infectious diseases are, however, difficult to forecast since different processes such as local adaptation to temperature, vector-pathogen interactions, and human-derived changes in landscape may play important roles in shaping the future dynamics of pathogen transmission., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Background: Mediterranean spotted fever (MSF) is a zoonotic disease caused by Rickettsia conorii and transmitted by Rhipicephalus sanguineus sensu lato. The aim of this study is to understand the epidemiology and trends regarding the disease in Spain, based on notifications to the Spanish National Epidemiology Surveillance Network (RENAVE) and the National Hospital Discharge Database (CMBD) between 2005 and 2015., Methods: We carried out a retrospective cross-sectional study of the cases and the outbreaks reported to the RENAVE and of those found in the CMBD between January 1 st , 2005 and December 31 st , 2015. We studied the characteristics of the cases and analyzed their spatio-temporal distribution., Results: 1603 cases notified to the RENAVE and 1789 cases registered in the CMBD were analyzed. The most affected group were men between 45 and 64. There were 8 MSF outbreaks during the study period. RENAVE registered lower rates until 2012, when it was decided that MSF in Spain would become a notifiable disease. Across the temporal series we saw that there was seasonality with an increase in cases in summer, and an overall upward trend according to the RENAVE data and descending according to the CMBD. The geographic distribution was heterogeneous throughout the territory, with maximum rates in La Rioja at 1.87 cases and 2.01 cases per 100,000 inhabitants according to the RENAVE and the CMBD, respectively., Conclusions: It is of great importance to continue monitoring the disease since it appears to be endemic throughout Spain. There is a need for a common strategy on monitoring and reporting, which would facilitate a more accurate picture on the MSF epidemiological scenario. Entomological information will be of added value., Competing Interests: Declaration of competing Interest None., (Copyright © 2019 Elsevier GmbH. All rights reserved.)

Background: Malaria was eliminated in Spain in 1964. Since then, more than 10,000 cases of malaria have been reported, mostly in travellers and migrants, making it the most frequently imported disease into this country. In order to improve knowledge on imported malaria cases characteristics, the two main malaria data sources were assessed: the national surveillance system and the hospital discharge database (CMBD)., Methods: Observational study using prospectively gathered surveillance data and CMBD records between 2002 and 2015. The average number of hospitalizations per year was calculated to assess temporal patterns. Socio-demographic, clinical and travel background information were analysed. Bivariate and multivariable statistical methods were employed to evaluate hospitalization risk, fatal outcome, continent of infection and chemoprophylaxis failure and their association with different factors., Results: A total of 9513 malaria hospital discharges and 7421 reported malaria cases were identified. The number of reported cases was below the number of hospitalizations during the whole study period, with a steady increase trend in both databases since 2008. Males aged 25-44 were the most represented in both data sources. Most frequent related co-diagnoses were anaemia (20.2%) and thrombocytopaenia (15.4%). The risks of fatal outcome increased with age and were associated with the parasite species (Plasmodium falciparum). The main place of infection was Africa (88.9%), particularly Equatorial Guinea (33.2%). Most reported cases were visiting friends and relatives (VFRs) and immigrants (70.2%). A significant increased likelihood of hospitalization was observed for children under 10 years (aOR:2.7; 95% CI 1.9-3.9), those infected by Plasmodium vivax (4.3; 95% CI 2.1-8.7) and travellers VFRs (1.4; 95% CI 1.1-1.7). Only 4% of cases reported a correct regime of chemoprophylaxis. Being male, over 15 years, VFRs, migrant and born in an endemic country were associated to increased risk of failure in preventive chemotherapy., Conclusions: The joint analysis of two data sources allowed for better characterization of imported malaria profile in Spain. Despite the availability of highly effective preventive measures, the preventable burden from malaria is high in Spain. Pre-travel advice and appropriately delivered preventive messages needs to be improved, particularly in migrants and VFRs.

Introduction and Objective: Zika virus disease is a challenge for public health due to its rapid spread and potential foetal complications. Although it is imported in Spain, there is a risk of autochthonous transmission due to Aedes albopictus presence. Zika disease and congenital cases have been under surveillance since 2016. The objective of this study is to explore the epidemiology of disease and pregnancies result., Material and Methods: A descriptive study was carried out into cases reported to the National Surveillance Network (RENAVE) during the 30/11/2015 to 31/12/2017 period. The case definition and the survey are included in the RENAVE protocol. The variables were: date; notifying region (Autonomous Community (AC)); pregnancy and its evolution; case classification; mode of transmission; country or region of infection; socio-demographical, clinical and microbiological data. A descriptive analysis of the cases and their distribution according to the other variables was carried out., Results: A total of 512 cases were reported by 17 ACs. 507 were non-congenital, of which 327 (64.5%) were women (52.5% of childbearing age). 403 cases (79.5%) corresponded to 2016 and 193 (38.1%) resided in regions with A. albopictus presence between May and October. 96.1% of imported cases were infected in America (51.7% while visiting relatives). Three cases (3.9%) of congenital Zika virus infection were detected among 77 pregnant women., Conclusions: The evolution of reported cases was in accordance with that of the epidemic in America. The largest group of travellers was young women who travelled to Latin America on family visits. Pregnancy monitoring resulted in the identification of Zika related foetal complications., (Copyright © 2019 Elsevier España, S.L.U. All rights reserved.)

Background: Environmental exposures are related to the risk of some types of cancer, and children are the most vulnerable group of people. This study seeks to present the methodological approaches used in the papers of our group about risk of childhood cancers in the vicinity of pollution sources (industrial and urban sites). A population-based case-control study of incident childhood cancers in Spain and their relationship with residential proximity to industrial and urban areas was designed. Two methodological approaches using mixed multiple unconditional logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs) were developed: (a) "near vs. far" analysis, where possible excess risks of cancers in children living near ("near") versus those living far ("far") from industrial and urban areas were assessed; and (b) "risk gradient" analysis, where the risk gradient in the vicinity of industries was assessed. For each one of the two approaches, three strategies of analysis were implemented: "joint", "stratified", and "individualized" analysis. Incident cases were obtained from the Spanish Registry of Childhood Cancer (between 1996 and 2011)., Results: Applying this methodology, associations between proximity (≤ 2 km) to specific industrial and urban zones and risk (OR; 95% CI) of leukemias (1.31; 1.04-1.65 for industrial areas, and 1.28; 1.00-1.53 for urban areas), neuroblastoma (2.12; 1.18-3.83 for both industrial and urban areas), and renal (2.02; 1.16-3.52 for industrial areas) and bone (4.02; 1.73-9.34 for urban areas) tumors have been suggested., Conclusions: The two methodological approaches were used as a very useful and flexible tool to analyze the excess risk of childhood cancers in the vicinity of industrial and urban areas, which can be extrapolated and generalized to other cancers and chronic diseases, and adapted to other types of pollution sources.

Background: Current evidence suggests that childhood leukaemia can be associated with residential traffic exposure; nevertheless, more results are needed to support this conclusion., Objectives: To ascertain the possible effects of residential proximity to road traffic on childhood leukaemia, taking into account traffic density, road proximity and the type of leukaemia (acute lymphoid leukaemia or acute myeloid leukaemia)., Methods: We conducted a population-based case-control study of childhood leukaemia in Spain, covering the period 1990-2011. It included 1061 incidence cases gathered from the Spanish National Childhood Cancer Registry and those Autonomous Regions with 100% coverage, and 6447 controls, individually matched by year of birth, sex and autonomous region of residence. Distances were computed from the respective participant's residential locations to the different types of roads and four different buffers. Using logistic regression, odds ratios (ORs) and 95% confidence intervals (95%CIs), were calculated for four different categories of distance to roads., Results: Cases of childhood leukaemia had more than three-fold increased odds of living at <50 m of the busiest motorways compared to controls (OR = 2.90; 95%CI = 1.30-6.49). The estimates for acute lymphoid leukaemia (ALL) were slightly higher (OR = 2.95; 95%CI = 1.22-7.14), while estimates for cases with the same address at birth and at diagnosis were lower (OR = 2.40; 95%CI = 0.70-8.30)., Conclusions: Our study agrees with the literature and furnishes some evidence that living near a busy motorway could be a risk factor for childhood leukaemia., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Introduction: Onchocerciasis is a chronic neglected tropical disease caused by the filarial nematode Onchocerca volvulus, which is endemic in Equatorial Guinea. The aim was to estimate the current spatial distribution of onchocerciasis, and its related factors, in Bioko Island after several years of mass drug administration and vector control activities, by using GIS technics., Methodology: The survey was carried out within the framework of a wider research project entitled "Strengthening the National Programme for Control of Onchocerciasis and other Filariasis in Equatorial Guinea". A structured questionnaire was designed to cover basic socio-demographic information and risk factors for onchocerciasis and the coordinates of household. the hydrographic network data to calculate the positive onchocerciasis rate was used. Poisson generalized linear model was used to explore the association between onchocerciasis and the following covariates: distance to the river, preventive practices, water source and household´s main source of income. Two different cluster analysis methods were used: Getis-Ord Gi statistic and SaTScan™ purely spatial statistic estimator., Results: The risk of onchocerciasis was higher for those who drank water from external sources (RR 25.3) than for those who drank home tap water (RR 8.0). The clusters with z-score higher were located at the east of the island. For 5 km and 1 km distances, one significant cluster in the east was detected (RR 5.91 and RR 7.15)., Conclusion: No environmental factors related with onchocerciasis were found, including proximity to rivers. This could be partially explained by the fact that the vector was eliminated in 2005., Competing Interests: No Conflict of Interest is declared, (Copyright (c) 2018 Diana Gómez-Barroso, Laura Moya, Zaida Herrador, Belén García, Justino Nguema, Policarpo Ncogo, Pilar Aparicio, Agustín Benito.)

Background: There is growing concern of how climate change could affect public health, due to the increase number of extreme climate events. Hence, the study of the role that climate events play on the distribution of waterborne diseases, like Hepatitis A, could be key for developing new prevention approaches., Objective: To investigate the association between climate factors and Hepatitis A in Spain between 2010 and 2014., Methods: Weekly Hepatitis A cases between 2010 and 2014 were obtained from the Spanish Epidemiology Surveillance Network. Climate variables (weekly cumulative rainfall, rainy days, storm days and snow days) were obtained from National Climatic Data Center (NOAA satellite and information Service of USA). Each municipality was assigned to the nearest weather station (N=73). A Mixed-Effects Poisson regression was performed to estimate Incidence Rate Ratios (IRR), including a time lag of 2, 3 and 4weeks (most probable incubation period for Hepatitis A)., Results: Rainfall higher than 90th percentile (extreme precipitation) was associated with increased number of Hepatitis A cases 2weeks (IRR=1.24 CI 95%=1.09-1.40) and 4weeks after the event (IRR=1.15 CI 95%=1.01-1.30). An extra rainy day increased the risk of Hepatitis A two weeks after (IRR=1.03 CI 95%=1.01-1.05). We found higher risk of Hepatitis A two weeks after each extra storm day (IRR=1.06 CI 95%=1.00-1.12), and lower risk with 3 and 4weeks' lag (IRR=0.93 CI 95%=0.88-0.99 for lag3; IRR=0.94 CI 95%=0.88-0.99 for lag 4)., Conclusions: There is an increased risk of Hepatitis A 2weeks after water-related climate events. Including meteorological information in surveillance systems might improve to develop early prevention strategies for waterborne diseases., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Background: The transmission of malaria is intense in the majority of the countries of sub-Saharan Africa, particularly in those that are located along the Equatorial strip. The present study aimed to describe the current distribution of malaria prevalence among children and its environment-related factors as well as to detect malaria spatial clusters in the district of Bata, in Equatorial Guinea., Methods: From June to August 2013 a representative cross-sectional survey using a multistage, stratified, cluster-selected sample was carried out of children in urban and rural areas of Bata District. All children were tested for malaria using rapid diagnostic tests (RDTs). Results were linked to each household by global position system data. Two cluster analysis methods were used: hot spot analysis using the Getis-Ord Gi statistic, and the SaTScan™ spatial statistic estimates, based on the assumption of a Poisson distribution to detect spatial clusters. In addition, univariate associations and Poisson regression model were used to explore the association between malaria prevalence at household level with different environmental factors., Results: A total of 1416 children aged 2 months to 15 years living in 417 households were included in this study. Malaria prevalence by RDTs was 47.53%, being highest in the age group 6-15 years (63.24%, p < 0.001). Those children living in rural areas were there malaria risk is greater (65.81%) (p < 0.001). Malaria prevalence was higher in those houses located <1 km from a river and <3 km to a forest (IRR: 1.31; 95% CI 1.13-1.51 and IRR: 1.44; 95% CI 1.25-1.66, respectively). Poisson regression analysis also showed a decrease in malaria prevalence with altitude (IRR: 0.73; 95% CI 0.62-0.86). A significant cluster inland of the district, in rural areas has been found., Conclusions: This study reveals a high prevalence of RDT-based malaria among children in Bata district. Those households situated in inland rural areas, near to a river, a green area and/or at low altitude were a risk factor for malaria. Spatial tools can help policy makers to promote new recommendations for malaria control.

Background: Central nervous system tumors (CNS) are the most frequent solid tumor in children. Causes of CNS tumors are mainly unknown and only 5% of the cases can be explained by genetic predisposition. We studied the effects of environmental exposure on the incidence of CNS tumors in children by subtype, according to exposure to industrial and/or urban environment, exposure to crops and according to socio-economic status of the child., Methods: We carried out a population-based case-control study of CNS tumors in Spain, covering 714 incident cases collected from the Spanish Registry of Childhood Tumors (period 1996-2011) and 4284 controls, individually matched by year of birth, sex, and autonomous region of residence. We built a covariate to approximate the exposure to industrial and/or urban environment and a covariate for the exposure to crops (GCI) using the coordinates of the home addresses of the children. We used the 2001 Census to obtain information about socio-economic status (SES). We fitted logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals (95%CIs)., Results: The results for all CNS tumors showed an excess risk (OR = 1.37; 95%CI = 1.09-1.73) for SES, i.e., children living in the least deprived areas had 37% more risk of CNS tumor than children living in the most deprived areas. For GCI, an increase of 10% in crop surface in the 1-km buffer around the residence implied an increase of 22% in the OR (OR = 1.22; 95%CI = 1.15-1.29). Children living in the intersection of industrial and urban areas could have a greater risk of CNS tumors than children who live outside these areas (OR = 1.20; 95%CI = 0.82-1.77). Living in urban areas (OR = 0.90; 95%CI = 0.65-1.24) or industrial areas (OR = 0.96; 95%CI = 0.81-1.77) did not seem to increase the risk for all CNS tumors together. By subtype, Astrocytomas, Intracranial and intraspinal embryonal tumors, and other gliomas showed similar results., Conclusion: Our results suggest that higher socioeconomic status and exposure to crops could increase the risk of CNS tumors in children.

Few epidemiologic studies have explored risk factors for bone tumors in children, and the role of environmental factors needs to be analyzed. Our objective was to ascertain the association between residential proximity to industrial plants and urban areas and risk of bone tumors in children, taking into account industrial groups and toxic pollutants released. A population-based case-control study of childhood bone cancer in Spain was carried out, covering 114 incident cases obtained from the Spanish Registry of Childhood Tumors (between 1996 and 2011), and 684 controls individually matched by sex, year of birth, and autonomous region of residence. Distances from the subject's residences to the 1271 industries and the 30 urban areas (towns) with ≥75,000 inhabitants located in the study area were computed. Unconditional logistic regression models were fitted to estimate odds ratios (ORs) and 95% confidence intervals (95%CIs) for categories of distance (from 1km to 3km) to industrial and urban areas, with adjustment for matching variables and sociodemographic indicators. Excess risk (OR; 95%CI) of bone tumors in children was detected for children close to industrial facilities as a whole (2.33; 1.17-4.63 at 3km) - particularly surface treatment of metals (OR=2.50; 95%CI=1.13-5.56 at 2km), production and processing of metals (OR=3.30; 95%CI=1.41-7.77 at 2.5km), urban waste-water treatment plants (OR=4.41; 95%CI=1.62-11.98 at 2km), hazardous waste (OR=4.63; 95%CI=1.37-15.61 at 2km), disposal or recycling of animal waste (OR=4.73; 95%CI=1.40-15.97 at 2km), cement and lime (OR=3.89; 95%CI=1.19-12.77 at 2.5km), and combustion installations (OR=3.85; 95%CI=1.39-10.66 at 3km)-, and urban areas (4.43; 1.80-10.92). These findings support the need for more detailed exposure assessment of certain toxics released by these facilities., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Background: Few epidemiologic studies have explored risk factors for rare tumors in children, and the role of environmental factors needs to be assessed., Objectives: To ascertain the effect of residential proximity to both industrial and urban areas on childhood cancer risk, taking industrial groups into account., Methods: We conducted a population-based case-control study of five childhood cancers in Spain (retinoblastoma, hepatic tumors, soft tissue sarcomas, germ cell tumors, and other epithelial neoplasms/melanomas), including 557 incident cases from the Spanish Registry of Childhood Tumors (period 1996-2011), and 3342 controls individually matched by year of birth, sex, and region of residence. Distances were computed from the residences to the 1271 industries and the 30 urban areas with ≥75,000 inhabitants located in the study area. Using logistic regression, odds ratios (ORs) and 95% confidence intervals (95%CIs) for categories of distance to industrial and urban pollution sources were calculated, with adjustment for matching variables and socioeconomic confounders., Results: Children living near industrial and urban areas as a whole showed no excess risk for any of the tumors analyzed. However, isolated statistical associations (OR; 95%CI) were found between retinoblastoma and proximity to industries involved in glass and mineral fibers (2.49; 1.01-6.12 at 3km) and organic chemical industries (2.54; 1.10-5.90 at 2km). Moreover, soft tissue sarcomas registered the lower risks in the environs of industries as a whole (0.59; 0.38-0.93 at 4km)., Conclusions: We have found isolated statistical associations between retinoblastoma and proximity to industries involved in glass and mineral fibers and organic chemical industries., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Background: Neuroblastoma is the most common extracranial solid tumor in children but its etiology is not clearly understood. While a small fraction of cases might be attributable to genetic factors, the role of environmental pollution factors needs to be assessed., Objectives: To ascertain the effect of residential proximity to both industrial and urban areas on neuroblastoma risk, taking into account industrial groups and toxic substances released., Methods: We conducted a population-based case-control study of neuroblastoma in Spain, including 398 incident cases gathered from the Spanish Registry of Childhood Tumors (period 1996-2011), and 2388 controls individually matched by year of birth, sex, and region of residence. Distances were computed from the respective subject's residences to the 1271 industries and the 30 urban areas with ≥75,000 inhabitants located in the study area. Using logistic regression, odds ratios (ORs) and 95% confidence intervals (95%CIs) for categories of distance (from 1km to 5km) to industrial and urban pollution sources were calculated, with adjustment for matching variables and socioeconomic confounders., Results: Excess risk (OR; 95%CI) of neuroblastoma was detected for the intersection between industrial and urban areas: (2.52; 1.20-5.30) for industrial distance of 1km, and (1.99; 1.17-3.37) for industrial distance of 2km. By industrial groups, excess risks were observed near 'Production of metals' (OR=2.05; 95%CI=1.16-3.64 at 1.5km), 'Surface treatment of metals' (OR=1.89; 95%CI=1.10-3.28 at 1km), 'Mines' (OR=5.82; 95%CI=1.04-32.43 at 1.5km), 'Explosives/pyrotechnics' (OR=4.04; 95%CI=1.31-12.42 at 4km), and 'Urban waste-water treatment plants' (OR=2.14; 95%CI=1.08-4.27 at 1.5km)., Conclusions: These findings support the need for more detailed exposure assessment of certain substances released by these industries., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Background: Childhood cancer is the main cause of disease-related death in children in Spain. Although little is known about the etiology, environmental factors are potential explanations for a fraction of the cases. Previous studies have shown pesticides to be associated with childhood cancer. The difficulty of collecting personal environmental exposure data is an important limitation; this lack of information about pesticides motivates the development of new methods to subrogate this exposure. We developed a crop exposure index based on geographic information to study the relationship between exposure to different types of crops and risk of childhood tumors., Methods: We conducted a population-based case-control study of childhood cancer covering 3350 cases and 20,365 controls in two Spanish regions. We used CORINE Land Cover to obtain data about agricultural land use. We created a 1 km buffer around every child and calculated the percentage of crop surface within the buffer (Global Crop Index) for total crops and for individual types of crops. We fitted mixed multiple unconditional logistic regression models by diagnostic group., Results: We found excess of risk among children living in the proximity of crops. For total crops our results showed excesses of risk for almost all diagnostic groups and increasing risk with increasing crop index value. Analyses by region and individual type of crop also showed excess of risk., Conclusion: The results suggest that living in the proximity of cultivated land could be a risk factor for several types of cancer in children.

Background: Few risk factors for childhood renal tumors are well established. While a small fraction of cases might be attributable to susceptibility genes and congenital anomalies, the role of environmental factors needs to be assessed., Objectives: To explore the possible association between residential proximity to environmental pollution sources (industrial and urban areas, and agricultural crops) and childhood renal cancer, taking into account industrial groups and toxic substances released., Methods: We conducted a population-based case-control study of childhood renal cancer in Spain, including 213 incident cases gathered from the Spanish Registry of Childhood Tumors (period 1996-2011), and 1278 controls individually matched by year of birth, sex, and region of residence. Distances were computed from the respective subject's residences to the 1271 industries, the 30 urban areas with ≥75,000 inhabitants, and the agricultural crops located in the study area. Using logistic regression, odds ratios (ORs) and 95% confidence intervals (95%CIs) for categories of distance to pollution sources were calculated, with adjustment for matching variables and socioeconomic confounders., Results: Excess risk (OR; 95%CI) of childhood renal tumors was observed for children living near (≤2.5km) industrial installations as a whole (1.97; 1.13-3.42) - particularly glass and mineral fibers (2.69; 1.19-6.08), galvanization (2.66; 1.14-6.22), hazardous waste (2.59; 1.25-5.37), ceramic (2.35; 1.06-5.21), surface treatment of metals (2.25; 1.24-4.08), organic chemical industry (2.22; 1.15-4.26), food and beverage sector (2.19; 1.18-4.07), urban and waste-water treatment plants (2.14; 1.07-4.30), and production and processing of metals (1.98; 1.03-3.82) -, and in the proximity of agricultural crops (3.16; 1.54-8.89 for children with percentage of crop surface ≥24.35% in a 1-km buffer around their residences)., Conclusions: Our study provides some epidemiological evidence that living near certain industrial areas and agricultural crops may be a risk factor for childhood renal cancer., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Objective: To estimate the areas of greatest density of road traffic accidents with fatalities at 24 hours per km(2)/year in Spain from 2008 to 2011, using a geographic information system., Methods: Accidents were geocodified using the road and kilometer points where they occurred. The average nearest neighbor was calculated to detect possible clusters and to obtain the bandwidth for kernel density estimation., Results: A total of 4775 accidents were analyzed, of which 73.3% occurred on conventional roads. The estimated average distance between accidents was 1,242 meters, and the average expected distance was 10,738 meters. The nearest neighbor index was 0.11, indicating that there were aggregations of accidents in space. A map showing the kernel density was obtained with a resolution of 1 km(2), which identified the areas of highest density., Conclusions: This methodology allowed a better approximation to locating accident risks by taking into account kilometer points. The map shows areas where there was a greater density of accidents. This could be an advantage in decision-making by the relevant authorities., (Copyright © 2014 SESPAS. Published by Elsevier Espana. All rights reserved.)

Background: Few risk factors for the childhood leukemia are well established. While a small fraction of cases of childhood leukemia might be partially attributable to some diseases or ionizing radiation exposure, the role of industrial and urban pollution also needs to be assessed., Objectives: To ascertain the possible effect of residential proximity to both industrial and urban areas on childhood leukemia, taking into account industrial groups and toxic substances released., Methods: We conducted a population-based case-control study of childhood leukemia in Spain, covering 638 incident cases gathered from the Spanish Registry of Childhood Tumors and for those Autonomous Regions with 100% coverage (period 1990-2011), and 13,188 controls, individually matched by year of birth, sex, and autonomous region of residence. Distances were computed from the respective subject's residences to the 1068 industries and the 157 urban areas with ≥10,000 inhabitants, located in the study area. Using logistic regression, odds ratios (ORs) and 95% confidence intervals (95%CIs) for categories of distance to industrial and urban pollution sources were calculated, with adjustment for matching variables., Results: Excess risk of childhood leukemia was observed for children living near (≤2.5 km) industries (OR=1.31; 95%CI=1.03-1.67) - particularly glass and mineral fibers (OR=2.42; 95%CI=1.49-3.92), surface treatment using organic solvents (OR=1.87; 95%CI=1.24-2.83), galvanization (OR=1.86; 95%CI=1.07-3.21), production and processing of metals (OR=1.69; 95%CI=1.22-2.34), and surface treatment of metals (OR=1.62; 95%CI=1.22-2.15) - , and urban areas (OR=1.36; 95%CI=1.02-1.80)., Conclusions: Our study furnishes some evidence that living in the proximity of industrial and urban sites may be a risk factor for childhood leukemia., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Objective: An excess of mortality was detected in Spain in February and March 2012 by the Spanish daily mortality surveillance system and the «European monitoring of excess mortality for public health action» program. The objective of this article was to determine whether this excess could be attributed to influenza in this period., Methods: Excess mortality from all causes from 2006 to 2012 were studied using time series in the Spanish daily mortality surveillance system, and Poisson regression in the European mortality surveillance system, as well as the FluMOMO model, which estimates the mortality attributable to influenza. Excess mortality due to influenza and pneumonia attributable to influenza were studied by a modification of the Serfling model. To detect the periods of excess, we compared observed and expected mortality., Results: In February and March 2012, both the Spanish daily mortality surveillance system and the European mortality surveillance system detected a mortality excess of 8,110 and 10,872 deaths (mortality ratio (MR): 1.22 (95% CI:1.21-1.23) and 1.32 (95% CI: 1.29-1.31), respectively). In the 2011-12 season, the FluMOMO model identified the maximum percentage (97%) of deaths attributable to influenza in people older than 64 years with respect to the mortality total associated with influenza (13,822 deaths). The rate of excess mortality due to influenza and pneumonia and respiratory causes in people older than 64 years, obtained by the Serfling model, also reached a peak in the 2011-2012 season: 18.07 and 77.20, deaths per 100,000 inhabitants, respectively., Conclusion: A significant increase in mortality in elderly people in Spain was detected by the Spanish daily mortality surveillance system and by the European mortality surveillance system in the winter of 2012, coinciding with a late influenza season, with a predominance of the A(H3N2) virus, and a cold wave in Spain. This study suggests that influenza could have been one of the main factors contributing to the mortality excess observed in the winter of 2012 in Spain., (Copyright © 2014 SESPAS. Published by Elsevier Espana. All rights reserved.)

Background: Childhood cancer was the leading cause of death among children aged 1-14 years for 2012 in Spain. Leukemia has the highest incidence, followed by tumors of the central nervous system (CNS) and lymphomas (Hodgkin lymphoma, HL, and Non-Hodgkin's lymphoma, NHL). Spatial distribution of childhood cancer cases has been under concern with the aim of identifying potential risk factors., Objective: The two objectives are to study overall spatial clustering and cluster detection of cases of the three main childhood cancer causes, looking to increase etiological knowledge., Methods: We ran a case-control study. The cases were children aged 0 to 14 diagnosed with leukemia, lymphomas (HL and NHL) or CNS neoplasm in five Spanish regions for the period 1996-2011. As a control group, we used a sample from the Birth Registry matching every case by year of birth, autonomous region of residence and sex with six controls. We geocoded and validated the address of the cases and controls. For our two objectives we used two different methodologies. For the first, for overall spatial clustering detection, we used the differences of K functions from the spatial point patterns perspective proposed by Diggle and Chetwynd and the second, for cluster detection, we used the spatial scan statistic proposed by Kulldorff with a level for statistical significance of 0.05., Results: We had 1062 cases of leukemia, 714 cases of CNS, 92 of HL and 246 of NHL. Accordingly we had 6 times the number of controls, 6372 controls for leukemia, 4284 controls for CNS, 552 controls for HL and 1476 controls for NHL. We found variations in the estimated empirical D(s) for the different regions and cancers, including some overall spatial clustering for specific regions and distances. We did not find statistically significant clusters., Conclusions: The variations in the estimated empirical D(s) for the different regions and cancers could be partially explained by the differences in the spatial distribution of the population; however, according to the literature, we cannot discard environmental hazards or infections agents in the etiology of these cancers.

Background: Cardiovascular disease is the leading cause of death worldwide, but its spatial distribution is not homogeneous. The objective of this study is to analyze the spatial pattern of mortality from these diseases for men and women, in the populated urban area (AUP) of the municipality of Madrid, and to identify spatial aggregations., Methods: An ecological study was carried out by census tract, for men and women in 2010. Standardized Mortality Ratio (SMR), Relative Risk Smoothing (RRS) and Posterior Probability (PP) were calculated to consider the spatial pattern of the disease. To identify spatial clusters the Moran index (Moran I) and the Local Index of Spatial Autocorrelation (LISA) were used. The results were mapped., Results: SMR higher than 1.1 was observed mainly in central areas among men and in peripheral areas among women. The PP that RRS was higher than 1 surpassed 0.8 in the center and in the periphery, in both men and women. Moran's I was 0.04 for men and 0.03 for women (p <0.05 in both cases)., Conclusions: Sex differences were observed in the spatial distribution of mortality cases. RME RRS and PP maps showed a heterogeneous pattern in men, whereas in women a clearer pattern was detected, with a relatively higher risk in peripheral areas of the AUP. The LISA method showed similar patterns to those previously observed.

Cluster detection has become an important part of the agenda of epidemiologists and public health authorities, the identification of high- and low-risk areas is fundamental in the definition of public health strategies and in the suggestion of potential risks factors. Currently, there are different cluster detection techniques available, the most popular being those using windows to scan the areas within the studied region. However, when these areas are heterogeneous in populations' sizes, scan window methods can lead to inaccurate conclusions. In order to perform cluster detection over heterogeneously populated areas, we developed a method not based on scanning windows but instead on standard mortality ratios (SMR) using irregular spatial aggregation (ISA). Its extension, i.e. irregular spatial aggregation with covariates (ISAC), includes covariates with residuals from Poisson regression. We compared the performance of the method with the flexible shaped spatial scan statistic (FlexScan) using mortality data for stomach and bladder cancer for 8,098 Spanish towns. The results show a collection of clusters for stomach and bladder cancer similar to that detected by ISA and FlexScan. However, in general, clusters detected by FlexScan were bigger and include towns with SMR, which were not statistically significant. For bladder cancer, clusters detected by ISAC differed from those detected by ISA and FlexScan in shape and location. The ISA and ISAC methods could be an alternative to the traditional scan window methods for cluster detection over aggregated data when the areas under study are heterogeneous in terms of population. The simplicity and flexibility of the methods make them more attractive than methods based on more complicated algorithms.

Using mortality data from National Institute of Statistics in Spain, we analyzed trends of infectious disease mortality rates in Spain during 1980-2011 to provide information on surveillance and control of infectious diseases. During the study period, 628,673 infectious disease-related deaths occurred, the annual change in the mortality rate was -1.6%, and the average infectious disease mortality rate was 48.5 deaths/100,000 population. Although the beginning of HIV/AIDS epidemic led to an increased mortality rate, a decreased rate was observed by the end of the twentieth century. By codes from the International Classification of Diseases, 9th revision, the most frequent underlying cause of death was pneumonia. Emergence and reemergence of infectious diseases continue to be public health problems despite reduced mortality rates produced by various interventions. Therefore, surveillance and control systems should be reinforced with a goal of providing reliable data for useful decision making.

A retrospective, space-time study of whooping cough cases reported to the Public Health Agency of Barcelona, Spain between the years 2000 and 2011 is presented. It is based on 633 individual whooping cough cases and the 2006 population census from the Spanish National Statistics Institute, stratified by age and sex at the census tract level. Cluster identification was attempted using space-time scan statistic assuming a Poisson distribution and restricting temporal extent to 7 days and spatial distance to 500 m. Statistical calculations were performed with Stata 11 and SatScan and mapping was performed with ArcGis 10.0. Only clusters showing statistical significance (P <0.05) were mapped. The most likely cluster identified included five census tracts located in three neighbourhoods in central Barcelona during the week from 17 to 23 August 2011. This cluster included five cases compared with the expected level of 0.0021 (relative risk = 2436, P <0.001). In addition, 11 secondary significant space-time clusters were detected with secondary clusters occurring at different times and localizations. Spatial statistics is felt to be useful by complementing epidemiological surveillance systems through visualizing excess in the number of cases in space and time and thus increase the possibility of identifying outbreaks not reported by the surveillance system.