Spatial aspects of MRSA epidemiology:A case study using stochastic simulation, kernel estimation and SaTScan

The identification of disease clusters in space or space-time is of vital importance for public health policy and action. In the case of methicillin-resistant Staphylococcus aureus (MRSA), it is particularly important to distinguish between community and health care-associated infections, and to identify reservoirs of infection. 832 cases of MRSA in the West Midlands (UK) were tested for clustering and evidence of community transmission, after being geo-located to the centroids of UK unit postcodes (postal areas roughly equivalent to Zip+4 zip code areas). An age-stratified analysis was also carried out at the coarser spatial resolution of UK Census Output Areas. Stochastic simulation and kernel density estimation were combined to identify significant local clusters of MRSA (p<0.025), which were supported by SaTScan spatial and spatio-temporal scan. In order to investigate local sampling effort, a spatial 'random labelling' approach was used, with MRSA as cases and MSSA (methicillin-sensitive S. aureus) as controls. Heavy sampling in general was a response to MRSA outbreaks, which in turn appeared to be associated with medical care environments. The significance of clusters identified by kernel estimation was independently supported by information on the locations and client groups of nursing homes, and by preliminary molecular typing of isolates. In the absence of occupational/ lifestyle data on patients, the assumption was made that an individual's location and consequent risk is adequately represented by their residential postcode. The problems of this assumption are discussed, with recommendations for future data collection.

Spatial and temporal clusters of Barmah Forest virus disease in Queensland, Australia

Objective To identify the spatial and temporal clusters of Barmah Forest virus (BFV) disease in Queensland in Australia, using geographical information systems (GIS) and spatial scan statistic (SaTScan). Methods We obtained BFV disease cases, population and statistical local areas boundary data from Queensland Health and Australian Bureau of Statistics respectively during 1992-2008 for Queensland. A retrospective Poisson-based analysis using SaTScan software and method was conducted in order to identify both purely spatial and space-time BFV disease high-rate clusters. A spatial cluster size of a proportion of the population and a 200km circle radius and varying time windows from 1 month to 12 months were chosen (for the space-time analysis). Results The spatial scan statistic detected a most likely significant purely spatial cluster (including 23 SLAs) and a most likely significant space-time cluster (including 24 SLAs) in approximately the same location. Significant secondary clusters were also identified from both the analyses in several locations. Conclusions This study provides evidence of the existence of statistically significant BFV disease clusters in Queensland, Australia. The study also demonstrated the relevance and applicability of SaTScan in analysing on-going surveillance data to identify clusters to facilitate the development of effective BFV disease prevention and control strategies in Queensland, Australia.

Space-time clusters of dengue fever in Bangladesh

Objective:  To examine the space-time clustering of dengue fever (DF) transmission in Bangladesh using geographical information system and spatial scan statistics (SaTScan). Methods:  We obtained data on monthly suspected DF cases and deaths by district in Bangladesh for the period of 2000–2009 from Directorate General of Health Services. Population and district boundary data of each district were collected from national census managed by Bangladesh Bureau of Statistics. To identify the space-time clusters of DF transmission a discrete Poisson model was performed using SaTScan software. Results:  Space-time distribution of DF transmission was clustered during three periods 2000–2002, 2003–2005 and 2006–2009. Dhaka was the most likely cluster for DF in all three periods. Several other districts were significant secondary clusters. However, the geographical range of DF transmission appears to have declined in Bangladesh over the last decade. Conclusion:  There were significant space-time clusters of DF in Bangladesh over the last decade. Our results would prompt future studies to explore how social and ecological factors may affect DF transmission and would also be useful for improving DF control and prevention programs in Bangladesh.

Outbreak detection algorithms for seasonal disease data : a case study using Ross River virus disease

Background Detection of outbreaks is an important part of disease surveillance. Although many algorithms have been designed for detecting outbreaks, few have been specifically assessed against diseases that have distinct seasonal incidence patterns, such as those caused by vector-borne pathogens. Methods We applied five previously reported outbreak detection algorithms to Ross River virus (RRV) disease data (1991-2007) for the four local government areas (LGAs) of Brisbane, Emerald, Redland and Townsville in Queensland, Australia. The methods used were the Early Aberration Reporting System (EARS) C1, C2 and C3 methods, negative binomial cusum (NBC), historical limits method (HLM), Poisson outbreak detection (POD) method and the purely temporal SaTScan analysis. Seasonally-adjusted variants of the NBC and SaTScan methods were developed. Some of the algorithms were applied using a range of parameter values, resulting in 17 variants of the five algorithms. Results The 9,188 RRV disease notifications that occurred in the four selected regions over the study period showed marked seasonality, which adversely affected the performance of some of the outbreak detection algorithms. Most of the methods examined were able to detect the same major events. The exception was the seasonally-adjusted NBC methods that detected an excess of short signals. The NBC, POD and temporal SaTScan algorithms were the only methods that consistently had high true positive rates and low false positive and false negative rates across the four study areas. The timeliness of outbreak signals generated by each method was also compared but there was no consistency across outbreaks and LGAs. Conclusions This study has highlighted several issues associated with applying outbreak detection algorithms to seasonal disease data. In lieu of a true gold standard, a quantitative comparison is difficult and caution should be taken when interpreting the true positives, false positives, sensitivity and specificity.

Hot spot detection and spatio-temporal dynamics of dengue in Queensland, Australia

Dengue has been a major public health concern in Australia since it re-emerged in Queensland in 1992-1993. This study explored spatio-temporal distribution and clustering of locally-acquired dengue cases in Queensland State, Australia and identified target areas for effective interventions. A computerised locally-acquired dengue case dataset was collected from Queensland Health for Queensland from 1993 to 2012. Descriptive spatial and temporal analyses were conducted using geographic information system tools and geostatistical techniques. Dengue hot spots were detected using SatScan method. Descriptive spatial analysis showed that a total of 2,398 locally-acquired dengue cases were recorded in central and northern regions of tropical Queensland. A seasonal pattern was observed with most of the cases occurring in autumn. Spatial and temporal variation of dengue cases was observed in the geographic areas affected by dengue over time. Tropical areas are potential high-risk areas for mosquito-borne diseases such as dengue. This study demonstrated that the locally-acquired dengue cases have exhibited a spatial and temporal variation over the past twenty years in tropical Queensland, Australia. There is a clear evidence for the existence of statistically significant clusters of dengue and these clusters varied over time. These findings enabled us to detect and target dengue clusters suggesting that the use of geospatial information can assist the health authority in planning dengue control activities and it would allow for better design and implementation of dengue management programs.

First Description of a Cluster of Acute/Subacute Paracoccidioidomycosis Cases and Its Association with a Climatic Anomaly

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Avaliação de fatores de risco de cães sororreagentes à leptospira spp. e sua distribuição espacial, em área territorial urbana

We aimed to evaluate the risk factor of serum reactive dogs to antileptospire agglutinin and their spatial distribution in an urban area. We collected 1,000 blood samples from dogs at 20 immunization centers, homogeneously distributed in the urban area (32 km²) of Botucatu, Sao Paulo, Brazil. Diagnosis was made by microscopic serum agglutination using 24 serovars of Leptospira spp. Statistical analysis was performed by Goodman's test with 5% of significance. Spatial clusters were tested by spatial scan statistic using SaTScan. We found that 17.9% dogs reacted to Leptospira spp. (p<0.0001). Dogs with street access (22.14%) were more reactive (p<0.05) than those without (14.83%). The scan test indicated only one significant cluster influenced by age. Incorporating each one of the co-variables gender, breed, and street access to age, we found that street access was the most important. Street access and age were the most important risk factors in the large number of reactive animals in the cluster, with the castellonis serovar being the most reactive in urban dogs. The identification of a cluster with more reactive dogs than expected allows local preventive measures to be taken.

Temporal, spatial and spatiotemporal analysis of the occurrence of visceral leishmaniasis in humans in the City of Birigui, State of São Paulo, from 1999 to 2012

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Suicide rates and income in Sao Paulo and Brazil: a temporal and spatial epidemiologic analysis from 1996 to 2008

Background: In a classical study, Durkheim noted a direct relation between suicide rates and wealth in the XIX century France. Since that time, several studies have verified this relationship. It is known that suicide rates are associated with income, although the direction of this association varies worldwide. Brazil presents a heterogeneous distribution of income and suicide across its territory; however, evaluation for an association between these variables has shown mixed results. We aimed to evaluate the relationship between suicide rates and income in Brazil, State of Sao Paulo (SP), and City of SP, considering geographical area and temporal trends. Methods: Data were extracted from the National and State official statistics departments. Three socioeconomic areas were considered according to income, from the wealthiest (area 1) to the poorest (area 3). We also considered three regions: country-wide (27 Brazilian States and 558 Brazilian micro-regions), state-wide (645 counties of SP State), and city-wide (96 districts of SP city). Relative risks (RR) were calculated among areas 1, 2, and 3 for all regions, in a cross-sectional approach. Then, we used Joinpoint analysis to explore the temporal trends of suicide rates and SaTScan to investigate geographical clusters of high/low suicide rates across the territory. Results: Suicide rates in Brazil, the State of SP, and the city of SP were 6.2, 6.6, and 5.4 per 100,000, respectively. Taking suicide rates of the poorest area (3) as reference, the RR for the wealthiest area was 1.64, 0.88, and 1.65 for Brazil, State of SP, and city of SP, respectively (p for trend <0.05 for all analyses). Spatial cluster of high suicide rates were identified at Brazilian southern (RR = 2.37), state of SP western (RR = 1.32), and city of SP central (RR = 1.65) regions. A direct association between income and suicide were found for Brazil (OR = 2.59) and the city of SP (OR = 1.07), and an inverse association for the state of SP (OR = 0.49). Conclusions: Temporospatial analyses revealed higher suicide rates in wealthier areas in Brazil and the city of SP and in poorer areas in the State of SP. We further discuss the role of socioeconomic characteristics for explaining these discrepancies and the importance of our findings in public health policies. Similar studies in other Brazilian States and developing countries are warranted.

Suicide rates and income in São Paulo and Brazil: a temporal and spatial epidemiologic analysis from 1996 to 2008

Abstract Background In a classical study, Durkheim noted a direct relation between suicide rates and wealth in the XIX century France. Since that time, several studies have verified this relationship. It is known that suicide rates are associated with income, although the direction of this association varies worldwide. Brazil presents a heterogeneous distribution of income and suicide across its territory; however, evaluation for an association between these variables has shown mixed results. We aimed to evaluate the relationship between suicide rates and income in Brazil, State of São Paulo (SP), and City of SP, considering geographical area and temporal trends. Methods Data were extracted from the National and State official statistics departments. Three socioeconomic areas were considered according to income, from the wealthiest (area 1) to the poorest (area 3). We also considered three regions: country-wide (27 Brazilian States and 558 Brazilian micro-regions), state-wide (645 counties of SP State), and city-wide (96 districts of SP city). Relative risks (RR) were calculated among areas 1, 2, and 3 for all regions, in a cross-sectional approach. Then, we used Joinpoint analysis to explore the temporal trends of suicide rates and SaTScan to investigate geographical clusters of high/low suicide rates across the territory. Results Suicide rates in Brazil, the State of SP, and the city of SP were 6.2, 6.6, and 5.4 per 100,000, respectively. Taking suicide rates of the poorest area (3) as reference, the RR for the wealthiest area was 1.64, 0.88, and 1.65 for Brazil, State of SP, and city of SP, respectively (p for trend <0.05 for all analyses). Spatial cluster of high suicide rates were identified at Brazilian southern (RR = 2.37), state of SP western (RR = 1.32), and city of SP central (RR = 1.65) regions. A direct association between income and suicide were found for Brazil (OR = 2.59) and the city of SP (OR = 1.07), and an inverse association for the state of SP (OR = 0.49). Conclusions Temporospatial analyses revealed higher suicide rates in wealthier areas in Brazil and the city of SP and in poorer areas in the State of SP. We further discuss the role of socioeconomic characteristics for explaining these discrepancies and the importance of our findings in public health policies. Similar studies in other Brazilian States and developing countries are warranted.

Toxics Release Inventory facilities and childhood cancer: Geographic information systems based approach

Purpose. To examine the association between living in proximity to Toxics Release Inventory (TRI) facilities and the incidence of childhood cancer in the State of Texas. ^ Design. This is a secondary data analysis utilizing the publicly available Toxics release inventory (TRI), maintained by the U.S. Environmental protection agency that lists the facilities that release any of the 650 TRI chemicals. Total childhood cancer cases and childhood cancer rate (age 0-14 years) by county, for the years 1995-2003 were used from the Texas cancer registry, available at the Texas department of State Health Services website. Setting: This study was limited to the children population of the State of Texas. ^ Method. Analysis was done using Stata version 9 and SPSS version 15.0. Satscan was used for geographical spatial clustering of childhood cancer cases based on county centroids using the Poisson clustering algorithm which adjusts for population density. Pictorial maps were created using MapInfo professional version 8.0. ^ Results. One hundred and twenty five counties had no TRI facilities in their region, while 129 facilities had at least one TRI facility. An increasing trend for number of facilities and total disposal was observed except for the highest category based on cancer rate quartiles. Linear regression analysis using log transformation for number of facilities and total disposal in predicting cancer rates was computed, however both these variables were not found to be significant predictors. Seven significant geographical spatial clusters of counties for high childhood cancer rates (p<0.05) were indicated. Binomial logistic regression by categorizing the cancer rate in to two groups (<=150 and >150) indicated an odds ratio of 1.58 (CI 1.127, 2.222) for the natural log of number of facilities. ^ Conclusion. We have used a unique methodology by combining GIS and spatial clustering techniques with existing statistical approaches in examining the association between living in proximity to TRI facilities and the incidence of childhood cancer in the State of Texas. Although a concrete association was not indicated, further studies are required examining specific TRI chemicals. Use of this information can enable the researchers and public to identify potential concerns, gain a better understanding of potential risks, and work with industry and government to reduce toxic chemical use, disposal or other releases and the risks associated with them. TRI data, in conjunction with other information, can be used as a starting point in evaluating exposures and risks. ^

Catchment area of oral cleft cases for two major hospitals in Harris County

Background. Orofacial clefts are among the most common birth defects and considered to be of complex etiology with both genetic and environmental factors.^ Objectives. The purpose of this study was to describe maternal and infant characteristics, examine the catchment area, and determine if there are any geospatial patterns among infants with an orofacial cleft delivered at two major hospitals in Harris County, The Woman's Hospital of Texas and Memorial Hermann Hospital-Texas Medical Center, from January 1, 2003 through December 31, 2007.^ Methods. Data were obtained from two major hospitals in Harris County and included all babies delivered in the period from 2003 through 2007 with an orofacial cleft. Residential addresses were mapped using MapInfo GIS software and the cluster analysis performed with SaTScan software.^ Results. Ninety-nine cases were identified spanning nine counties. 26% of cases resided within a 5-mile radius of the Texas Medical Center. Birth rates ranged from 1.4 to 16.5 per 10,000 total births. A cluster was identified in southwest Harris County, however, it was not significant (p=0.066).^ Conclusion. This study encourages further focus on linking cleft cases to environmental factors in order to determine potential risks. ^

A community-based assessment of the spatiotemporal distribution of influenza in a U.S.-Mexico border county during the spring 2009 novel H1N1 swine-origin influenza A pandemic

This study retrospectively evaluated the spatial and temporal disease patterns associated with influenza-like illness (ILI), positive rapid influenza antigen detection tests (RIDT), and confirmed H1N1 S-OIV cases reported to the Cameron County Department of Health and Human Services between April 26 and May 13, 2009 using the space-time permutation scan statistic software SaTScan in conjunction with geographical information system (GIS) software ArcGIS 9.3. The rate and age-adjusted relative risk of each influenza measure was calculated and a cluster analysis was conducted to determine the geographic regions with statistically higher incidence of disease. A Poisson distribution model was developed to identify the effect that socioeconomic status, population density, and certain population attributes of a census block-group had on that area's frequency of S-OIV confirmed cases over the entire outbreak. Predominant among the spatiotemporal analyses of ILI, RIDT and S-OIV cases in Cameron County is the consistent pattern of a high concentration of cases along the southern border with Mexico. These findings in conjunction with the slight northward space-time shifts of ILI and RIDT cluster centers highlight the southern border as the primary site for public health interventions. Finally, the community-based multiple regression model revealed that three factors—percentage of the population under age 15, average household size, and the number of high school graduates over age 25—were significantly associated with laboratory-confirmed S-OIV in the Lower Rio Grande Valley. Together, these findings underscore the need for community-based surveillance, improve our understanding of the distribution of the burden of influenza within the community, and have implications for vaccination and community outreach initiatives.^

Ross River virus disease clusters and spatial relationship with mosquito biting exposure in Redland Shire, southern Queensland, Australia

The spatial heterogeneity in the risk of Ross River virus (family Togaviridae, genus Alphavirus, RRV) disease, the most common mosquito-borne disease in Australia, was examined in Redland Shire in southern Queensland, Australia. Disease cases, complaints from residents of intense mosquito biting exposure, and human population data were mapped using a geographic information system. Surface maps of RRV disease age-sex standardized morbidity ratios and mosquito biting complaint morbidity ratios were created. To determine whether there was significant spatial variation in disease and complaint patterns, a spatial scan analysis method was used to test whether the number of cases and complaints was distributed according to underlying population at risk. Several noncontiguous areas in proximity to productive saline water habitats of Aedes vigilax (Skuse), a recognized vector of RRV, had higher than expected numbers of RRV disease cases and complaints. Disease rates in human populations in areas which had high numbers of adult Ae. vigilax in carbon dioxide- and octenol-baited light traps were up to 2.9 times those in areas that rarely had high numbers of mosquitoes. It was estimated that targeted control of adult Ae. vigilax in these high-risk areas could potentially reduce the RRV disease incidence by an average of 13.6%. Spatial correlation was found between RRV disease risk and complaints from residents of mosquito biting. Based on historical patterns of RRV transmission throughout Redland Shire and estimated future human population growth in areas with higher than average RRV disease incidence, it was estimated that RRV incidence rates will increase by 8% between 2001 and 2021. The use of arbitrary administrative areas that ranged in size from 4.6 to 318.3 km2, has the potential to mask any small scale heterogeneity in disease patterns. With the availability of georeferenced data sets and high-resolution imagery, it is becoming more feasible to undertake spatial analyses at relatively small scales.