An epidemiological and gis-based analysis of mortality in selected rare diseases
- Autores: Germán Sánchez Díaz
- Directores de la Tesis: Francisco Javier Escobar Martínez (dir. tes.), Verónica Alonso Ferreira (codir. tes.)
- Lectura: En la Universidad de Alcalá ( España ) en 2019
- Idioma: español
- Tribunal Calificador de la Tesis: Francisco Javier Salas Rey (presid.), Diana Gómez Barroso (secret.), Luisa Nicia Borrell Mejía (voc.)
- Programa de doctorado: Programa de Doctorado en Tecnologías de la Información Geográfica por la Universidad de Alcalá
- Materias:
- Enlaces
- Tesis en acceso abierto en: TESEO
- Resumen
- español
Introducción: Las Enfermedades Raras (ER) son aquellas que no superan los 5 casos por cada 10000 habitantes en la Unión Europea (UE). Se estima que podrían estar afectadas entre 27 y 36 millones de personas por ER en la UE (entre un 6 y un 8% de la población), por lo que se ha convertido en un problema de Salud Pública. Las características propias de este conjunto de más 6000 enfermedades (heterogeneidad, bajo número de casos, dispersión de casos, dificultades de codificación, entre otras) hacen más difícil su investigación epidemiológica. En el caso de la mortalidad, de la que existen datos a nivel nacional de base poblacional y estandarizada en su metodología de registro, su seguimiento es de gran utilidad ya que su reducción es un objetivo perseguido por las políticas de Salud Pública. Los Sistemas de Información Geográfica, gracias a su capacidad para la integración, análisis y representación de datos geoespaciales, se muestran como una valiosa herramienta para estudiar su distribución y evolución espacio-temporal. El conocimiento de patrones geográficos de la mortalidad y su relación con otras variables espaciales contribuyen a describir de forma más precisa su epidemiología. Además, permite aportar importantes elementos de base a la hora de establecer relaciones causales entre factores ambientales y las enfermedades analizadas.
Objetivos: El principal objetivo de esta tesis doctoral es profundizar en el conocimiento de la mortalidad debida a ER en España desde una perspectiva geográfica y epidemiológica. Los objetivos específicos se detallan a continuación:
Objetivo 1: Identificar las problemáticas actuales en la elección de la unidad geográfica para el trabajo con ER y proporcionar recomendaciones para escoger la unidad más apropiada buscando el mejor resultado epidemiológico y cartográfico.
Objetivo 2: Evaluar las tendencias temporales de la mortalidad debida a ER a lo largo de tres décadas en España. Para lograr este objetivo se seleccionarán las siguientes ER: Enfermedad de Huntington (EH) y Granulomatosis con poliangeitis (GPA).
Objetivo 3: Identificar los patrones geográficos de la mortalidad debida a ER para observar su variación entre las diferentes unidades geográficas en España. Este objetivo será desarrollado con las siguientes enfermedades: EH, GPA y enfermedades de la neuronas motoras (ENM).
Objetivo 4: Describir las características de los metales pesados emitidos a cuencas hidrográficas, así como explorar las posibles asociaciones entre su localización y la mortalidad atribuida a ENM con el fin de ayudar a entender su etiología.
Métodos: Esta tesis doctoral se circunscribe al territorio español y a un periodo máximo comprendido entre los años 1984 y 2016. Se utilizó como base del estudio la estadística de defunciones proporcionada de forma anual por el Instituto Nacional de Estadística. Se seleccionaron aquellos fallecimientos atribuidos a cada una de las tres ER, teniendo en cuenta la causa básica de defunción. Los datos de emisión de metales pesados a las cuencas hidrográficas, desde complejos industriales así como otras instalaciones, fueron recogidos del European Pollutant Release and Transfer Register (E-PRTR). Los ficheros de las unidades geográficas e hidrografía en formato shapefile fueron descargados del Instituto Geográfico Nacional. En primer lugar se analizó la idoneidad de las unidades de agregación disponibles (municipios, comarcas y provincias) mediante la comparación epidemiológica y cartográfica de diferentes indicadores como la Razón de Mortalidad Estandarizada (RME), RME suavizada, así como indicadores locales de asociación espacial (LISA). En cuanto a indicadores epidemiológicos de mortalidad, se analizó la evolución temporal de la EH y GPA mediante las tasas de mortalidad ajustada por edad (TAE) así como TAE suavizada. Las tendencias fueron evaluadas mediante una regresión joinpoint. Respecto a la distribución geográfica de la mortalidad para las tres enfermedades seleccionadas, se calcularon las RMEs y sus respectivos intervalos de confianza al 95%. La RME se define como el ratio entre los fallecimientos observados y los esperados en una unidad geográfica. También se calcularon las RMEs suavizadas y su probabilidad posterior. El proceso de suavizado tiene en cuenta los valores observados y esperados de las unidades geográficas vecinas y ayuda a encontrar patrones espaciales en el riesgo de fallecimiento. Por último, se codificaron las instalaciones que habían emitido metales pesados a ríos. Para estudiar su influencia en los fallecimientos atribuidos a ENM, se consideraron municipios expuestos aquellos que estaban incluidos 20 kilómetros aguas abajo desde el punto de emisión, y se comparó la mortalidad con los no expuestos mediante un modelo de regresión de Poisson. Todos los resultados epidemiológicos de carácter espacial fueron representados cartográficamente.
Resultados: El Problema de la Unidad Espacial Modificable (PUEM) se mostró relevante al analizar diferentes unidades de agregación para el estudio de las ER, ya que se han encontrado diferencias significativas al comparar indicadores epidemiológicos entre las tres unidades de análisis. El número de casos por enfermedad es determinante para elegir el nivel de agregación, siendo el nivel comarcal en España el que mejor muestra la variabilidad de la mortalidad manteniendo un buen detalle espacial. Respecto a la evolución de las TAE de mortalidad para las ER analizadas, se han observado cambios significativos en las dos enfermedades analizadas (EH y GPA). La TAE para la EH se incrementó anualmente un 3,44% (p<0,001) durante todo el periodo, mientras que en GPA se detectó un incremento anual del 20,6% desde 1984 hasta 1992 y, posteriormente un descenso del 1,91% (p<0,05). Tanto para GPA como EH, la edad media de fallecimiento aumentó un 0,78% y un 0,59% anual (p<0,05) durante sus respectivos periodos de estudio. Respecto a la mortalidad espacial, las tres ER estudiadas mostraron variabilidad en el territorio español. En EH se encontraron RME y RME suavizadas más altas en algunas comarcas del suroeste peninsular. Para GPA se encontró un mayor riesgo de fallecimiento en dos comarcas del norte y dos del sur. En el caso de ENM, se encontró un mayor riesgo de defunción en municipios localizados al norte del país en comparación con municipios situados en la mitad sur. En el análisis de regresión de Poisson de la variación de mortalidad por ENM y emisión de metales pesados a cuencas hidrográficas, se encontró que la mortalidad por este grupo de enfermedades fue un 18,4% mayor en los municipios expuestos que en los no expuestos, tanto de forma global así como en los análisis individuales por tipo de metal.
Conclusiones: Las ER suscitan cada vez mayor interés debido al alto impacto que tienen sobre los pacientes afectados, familiares y sociedad en general. Por ello, es necesario mejorar la información disponible mediante la investigación a través de perspectivas multidisciplinares. Esta tesis doctoral ha enriquecido el conocimiento de la variabilidad temporal y espacial de la mortalidad atribuida a ER añadiendo un enfoque geográfico. Se han analizado algunas de las dificultades que presenta el estudio epidemiológico de estas enfermedades poco prevalentes debido al bajo número de casos y la necesidad de agregación espacial que provoca como consecuencia el PUEM. También se ha aportado información relevante sobre la variabilidad temporal y geográfica de la mortalidad para las ER seleccionadas, además de explorar posibles asociaciones con factores geográficos. Los resultados de este estudio ofrecen una información muy valiosa para la planificación sanitaria, ya que apoyan y justifican acciones de prevención y de asistencia en aquellas zonas donde se han detectado más desigualdades en el riesgo de defunción debido a estas enfermedades. Además, al relacionarlo con factores ambientales sospechosos se ofrecen pistas para futuros estudios que puedan ahondar en el descubrimiento de las causas que provocan las ER.
- English
Introduction: Rare diseases (RD) are defined as those affecting no more than 5 in 10,000 people within the European Union (EU). It is estimated that between 27 and 36 million people might be affected by RD in the EU (which amounts to between 6 and 8% of its total population), with them consequently having become a Public Health issue. The epidemiological research of the over 6,000 accounted for RD proves harder than that of more frequent diseases, due to their specific characteristics (e.g. heterogeneity, low number of cases, case dispersion or difficulties in codification, to name a few). When it comes to mortality, the availability of national population-based data and standardized registries are very useful when researching RD, since the reduction of mortality is a public health policy goal. Geographic Information Systems allowing for the integration, analysis and representation of geospatial data have proved to be a valuable tool with which to study the mortality distribution and its spatio-temporal evolution. Knowing the geographical patterns of mortality and their relationships with other spatial variables helps more accurately describe its epidemiology, as well as providing important baseline elements when searching for causal relationships between environmental factors and the diseases object of analysis.
Objective: The general aim of this doctoral thesis is to deepen the knowledge on mortality attributed to RDs in Spain from a geographical and epidemiological perspective. The specific objectives are the following:
Objective 1: To identify difficulties related to the optimal geographic units of representation when working with RDs and to provide sound advice for the election of the most appropriate unit, aiming for a better epidemiological and cartographic result.
Objective 2: To evaluate the Spanish population’s time trends in mortality due to RDs over three decades. This target will be achieved by selecting Huntington’s disease (HD) and Granulomatosis with polyangiitis (GPA).
Objective 3: To identify geographic patterns in mortality attributed to RDs, in order to trace variations in mortality among spatial units throughout Spain. This objective will be tackled through the analysis of GPA, HD and Motor neuron diseases (MND) as the RDs of reference.
Objective 4: To describe the characteristics of heavy metal emissions to river basins, and to explore the possible correlation that exists between their location and the mortality attributed to MND, in order to help understand the etiology of this group of diseases.
Methodology: This doctoral thesis has the Spanish territory as its study framework. It covers the years 1984 to 2016, at maximum, and uses the annual death registry provided by the National Statistics Institute as the basis for analysis. The decease-related data relevant to each of the three analyzed RDs were selected taking into account the underlying cause of death, whilst data relative to the heavy metal emissions to river basins were collected from the European Pollutant Release and Transfer Register (E-PRTR). Shapefiles containing information on geographic units and hydrography were downloaded from the National Mapping Agency. As a starting point, the suitability of the available aggregation units (municipalities, districts and provinces) was analyzed, by means of the epidemiological and geospatial comparison of different indicators, such as the Standardized Mortality Ratio (SMR), smoothed SMR and local indicators of spatial association (LISA). Secondly, regarding epidemiological mortality indicators, time trends in HD and GPA were analyzed by calculating age-adjusted mortality rates, and subsequently smoothed. Time trends were assessed by a joinpoint regression. As to the geographical distribution of mortality for the three selected diseases, the SMRs and their 95% confidence intervals were calculated. SMR is defined as the existing ratio between observed and expected deaths within a geographical unit. Smoothed SMR and its posterior probability were also calculated. The smoothing process takes into account the observed and expected values of the neighboring geographical units and helps identify spatial patterns of death risk. Lastly, the sites that release heavy metal emissions into the rivers were geocoded. In order to study their influence on MND deaths, we took exposed municipalities to be those downstream from the point of emission, within 20 km of it. The mortality of the exposed municipalities was compared to that of the non-exposed ones through a Poisson regression model. The obtained epidemiological results of a spatial nature were represented cartographically.
Results: The modifiable areal unit problem (MAUP) proved relevant when analyzing different aggregation units for the study of RDs, since significant differences were found when comparing epidemiological indicators among the three existing units. The number of cases occurring due to disease was a determining factor when choosing the most appropriate level of aggregation. The district level (comarca) has proved the most appropriate to display the variability of mortality in Spain, while still displaying relevant spatial detail. Regarding the evolution of mortality rates for the analyzed RDs, significant time trends were identified in both HD and GPA within the periods of study. The mortality rate from HD increased annually by 3.44% (p <0.001) throughout, while that for GPA displayed an annual increase of 20.6% between 1984 and 1992, with a 1.91% annual increase thereafter. (p <0.05). The average ages of death from GPA and EH rose by 0.78% and 0.59% (p <0.05), respectively, within the studied timeframe. As for the spatial analysis of RD-related mortality, the three reference diseases showed geographical variability within the Spanish territory. Regarding HD, higher SMRs and smoothed SMRs were found in some districts to the Southwest of the peninsula. GPA posed an increased death risk for the inhabitants of two specific districts in the North of Spain, and for a further two in the South, while in the case of MND, a comparison between municipalities showed those located in the North to be at a higher risk of death from the disease than those to the South.
The Poisson regression analysis conducted on the variation of mortality attributed to MND in relation to the emission of heavy metals to river basins showed that MND mortality rates were 18.4% higher in the exposed municipalities than it was in non-exposed ones. The figure remained relevant both globally and for the analyses carried out individually for each type of metal.
Conclusions: RDs are of increasing interest due to the high impact they have on affected patients, their families and society in general. Therefore, it is necessary to improve the information available on them through research, carried out over multidisciplinary perspectives. This doctoral thesis has enriched the body of knowledge on spatial variability in mortality from RDs, by adding a geographical approach to its analysis. Some of the difficulties presented by the epidemiological study of these low prevalence diseases have been analyzed. The low number of cases and the need for spatial aggregation have the MAUP as a consequence. Relevant information regarding the temporal and geographic variability in mortality for selected RDs has been provided, as well as having explored the possible associations with geographic factors. The results of this study offer very valuable information for health planning, since they support and justify the need for preventive and assistive action in those areas where there have been more inequalities in the risk of death due to these diseases. In addition, by establishing a correlation between mortality from RDs and suspect environmental factors, we offer clues for future studies that may delve into the discovery of RD causes.
- español
