INTRODUCCIÓN
En la actualidad, la osteoporosis se define como "una enfermedad sistémica del esqueleto que se caracteriza por una masa ósea disminuida y una alteración estructural del tejido óseo que determina una disminución en la resistencia ósea derivando en un incremento importante de la fragilidad y susceptibilidad a las fracturas"1. Los factores de riesgo clásicos asociados al desarrollo de la osteoporosis son la edad, antecedentes de fractura previa o historial familiar de osteoporosis, y deficiencia estrogénica prolongada en el tiempo2.
Por otro lado, uno de los factores de riesgo de mortalidad más importante a nivel poblacional es la contaminación del aire. Las consecuencias derivadas de la contaminación del aire suponen un elevado impacto económico y social. En 2015, los gastos derivados de la morbimortalidad asociada a contaminación alcanzaron los 21 mil millones de dólares a nivel mundial, y se estima un número de muertes prematuras entre 6 y 9 millones de personas en 2060 debidas a la contaminación de aire exterior3. Se ha demostrado que la contaminación atmosférica tiene un impacto directo en la salud provocando diversos efectos adversos4. La relación entre la contaminación del aire y las enfermedades pulmonares ambientales se ha reportado en numerosos estudios4-6. Así mismo, existe una fuerte evidencia científica que relaciona la pobre calidad del aire en distintos entornos laborales con el desarrollo de distintas enfermedades respiratorias7,8.
En los años 80 se publicaron los primeros estudios que ponían de manifiesto una asociación entre calidad del aire y la calidad ósea observándose una incidencia de fracturas significativamente mayor en personas pertenecientes a áreas urbanas frente a aquellas residentes en áreas rurales9-11. Desde entonces, los estudios en esta línea han aumentado de forma notable a lo largo de los últimos años determinando que la exposición prolongada a la contaminación del aire está asociada a una disminución de la calidad ósea12-14, siendo considerado como un factor de riesgo modificable de osteoporosis y de fracturas15. Sin embargo, una revisión sistemática publicada recientemente, apunta a la existencia de asociaciones inconsistentes entre la contaminación atmosférica y el riesgo de osteoporosis, que podrían explicarse por la heterogeneidad en las características de los sujetos que participan en los estudios publicados en este sentido16.
A pesar de que la evidencia científica sugiere la existencia de una relación entre enfermedades pulmonares y osteoporosis, actualmente no se presta especial atención a la calidad del aire como un potencial problema de salud ósea. Considerando el incremento en la contaminación del aire durante los últimos años, es importante profundizar en este tema con el objetivo de conocer las vías cruzadas entre los efectos nocivos de la contaminación del aire y la salud ósea.
En este contexto, la identificación de marcadores genéticos robustos que estén presentes tanto en enfermedades pulmonares ambientales y ocupacionales relacionadas con la contaminación y la calidad del aire, como en enfermedades relacionadas con la calidad ósea, permitirá el desarrollo de estrategias preventivas y terapéuticas enfocadas a la población de mayor riesgo. Para este tipo de estudios, la bioinformática juega un importante papel, permitiendo identificar mediante el uso de distintas herramientas y algoritmos, moléculas que pueden actuar como biomarcadores comunes entre distintas patologías. En este sentido, el propósito de este trabajo fue construir y analizar una red de interacciones proteína-proteína relacionando genes implicados en distintas enfermedades pulmonares ambientales y ocupacionales relacionadas con la contaminación y la calidad del aire, con genes asociados a desarrollo de osteoporosis, para identificar biomarcadores comunes en ambas enfermedades.
MATERIAL Y MÉTODOS
Recopilación de datos
En primer lugar, se realizó una búsqueda bibliográfica para obtener una selección de términos unificados para determinar las enfermedades pulmonares y ocupacionales. Las bases de datos de terminología médica unificada consultadas fueron Unified Medical Language System (UMLS, https://www.nlm.nih.gov/research/umls/index.html) y Medical Subject Headings (MeSH, https://www.nlm.nih.gov/mesh/meshhome.html). Posteriormente, cada una de las enfermedades seleccionadas, fueron usadas para determinar los genes asociados a enfermedad (GDA) mediante la utilización de la plataforma DisGeNET (http://www.disgenet.org/home/), la cual integra información sobre relación entre genes y enfermedades a partir de varias fuentes de datos públicos y de literatura sobre expresión génica, biomarcadores, asociación entre variantes y enfermedades, polimorfismos de un solo nucleótido y asociaciones de fenotipos clínicos con las enfermedades correspondientes17.
Selección de genes semilla
Una vez obtenidos tanto los genes asociados a las diferentes enfermedades pulmonares (de origen ambiental u ocupacional) y a la osteoporosis, se realizó un filtrado inicial seleccionando aquellos genes comunes a osteoporosis y con al menos una de las enfermedades pulmonares enmarcadas dentro de este estudio. En una segunda fase, se seleccionaron aquellos genes vinculados con la osteoporosis y con dos enfermedades pulmonares para generar la red de interacción de proteínas (PPI).
Construcción de redes PPI
Para construir la red PPI de los genes asociados a las enfermedades objeto de estudio, se utilizó la aplicación STRING (versión 1.6.0, released: 8 Sep 2020, http://apps.cytoscape.org/apps/stringapp)18 dentro de la plataforma Cytoscape (versión 3.8.2, https://cytoscape.org/)19. El valor de punto de corte de confianza se estableció en 0,7, considerándose significativas únicamente las interacciones con un valor mayor o igual a 0,7. Las proteínas por debajo de este punto de corte y aquellas sin relación con otras proteínas, fueron descartadas.
Estudio de centralidad e identificación de genes
La obtención de los genes centrales, es decir, aquellos nodos importantes con un elevado número de interacciones hacia otros nodos, fueron seleccionados utilizando el complemento CytoHubba20 en Cytoscape. Se utilizaron cinco métodos de cálculo: Degree, Betweennes, Maximal Clique Centrality (MCC), Bottleneck y Closeness21. Los genes comunes derivados de estos seis algoritmos representan genes candidatos clave con importantes funciones reguladoras biológicas.
RESULTADOS
Recopilación de datos
Tras la revisión de términos médicos unificados, fueron seleccionadas las enfermedades incluidas en la tabla 1. Estas se utilizaron para comparar los genes asociados a cada una de ellas con la osteoporosis. No fueron incluidas enfermedades en las que no se hubieran identificado genes asociados a la enfermedad.
Selección de genes semilla
Tras la aplicación de los criterios establecidos para la inclusión de genes en la red PPI, se obtuvieron 157 genes vinculados con osteoporosis y al menos con dos enfermedades pulmonares. La lista completa puede visualizarse en el material suplementario.
Construcción de redes PPI
La red PPI generada incluyó 152 proteínas (nodos) y 1.378 relaciones (ejes). De este conjunto, 12 proteínas fueron descartadas por no estar vinculadas con otras proteínas (figura 1).
Estudio de centralidad e identificación de genes
Se determinaron 10 genes centrales mediante la aplicación CytoHubba a partir de los métodos de centralidad detallados en la tabla 2. A continuación se seleccionaron los 5 genes que aparecían con más frecuencia en los diferentes algoritmos de centralidad (tabla 3), los cuales compartían su presencia en todos los algoritmos utilizados, evidenciando su importancia dentro de la red de interacción construida.
DISCUSIÓN
Cualquier informe anual sobre la calidad del aire revela conclusiones nada halagüeñas. La Organización Mundial de la Salud advierte constantemente sobre el incumplimiento de los estándares recomendados para lograr una calidad de aire saludable. Una pobre calidad del aire tiene una gran repercusión sanitaria llegando a provocar 4,2 millones de muertes prematuras al año en todo el mundo, asociadas principalmente a la exposición a pequeñas partículas con un diámetro menor o igual a 2,5 micrones22. Los principales efectos sobre la salud que tiene la exposición a un entorno con aire de mala calidad o contaminado son conocidos y han sido demostrados científicamente durante muchos años. Entre ellos cabe destacar el cáncer de pulmón, enfermedades respiratorias y cardiovasculares, ictus, etc. Sin embargo, existen otros efectos menos conocidos o con escasa evidencia científica asociada a la falta de estudios en esta área.
Debido al envejecimiento de la población, la osteoporosis es considerada como una de las enfermedades de mayor prevalencia en la población a nivel mundial. Actualmente, la asociación entre una mala calidad del aire y la osteoporosis no está bien definida en la literatura científica, pero comienza a tener una gran importancia debido a la repercusión sanitaria que supone el elevado porcentaje de población que presenta problemas osteoporóticos o baja calidad ósea.
Los resultados de este estudio muestran un conjunto de genes asociados a la presencia tanto de enfermedades pulmonares (derivadas de contaminación ambiental o asociadas a determinadas actividades ocupacionales) como de osteoporosis. Las enfermedades pulmonares son causadas principalmente por una exposición repetida en el tiempo de irritantes químicos, alérgenos o toxinas, pudiendo provocar efectos duraderos en el individuo.
La búsqueda bibliográfica de enfermedades pulmonares de origen ambiental u ocupacional, ofreció un conjunto de enfermedades que, sin estar a priori relacionadas con la osteoporosis, podrían compartir uno o varios genes vinculados con procesos osteoporóticos. Así, encontramos la antracosis, una afección pulmonar producida por una exposición e inhalación de polvo con alta concentración en carbón; la asbestosis, en la que el amianto es la principal partícula inhalada; la bisinosis, producida por la inhalación del polvo de algodón entre otras partículas de origen vegetal; la enfermedad pulmonar obstructiva crónica (EPOC), consecuencia de una exposición a largo plazo a irritantes como el aire contaminado, vapores químicos o polvo, entre otros; el mesotelioma, una forma de cáncer específico del mesotelio pulmonar que suele estar relacionado principalmente con la exposición ocupacional a asbestos; la silicosis, provocada por la inhalación de polvo con partículas de sílice; las enfermedades pulmonares intersticiales, que engloban un conjunto de enfermedades pulmonares la mayoría de las cuales están causadas por cicatrizaciones progresivas del tejido pulmonar debidas a exposiciones prolongadas a sustancias peligrosas (como determinadas sustancias orgánicas, maderas, metales, agentes infecciosos como los virus, fármacos, etc.); el asma ocupacional, producido por la inflamación de las vías respiratorias pulmonares a causa de inhalación de sustancias producidas en el lugar de trabajo como polvo de madera, hongos, y/o sustancias químicas, entre otras y neumonitis por hipersensibilidad, provocada por una exposición a una gran cantidad de partículas orgánicas como hongos o bacterias, lo que provoca una importante respuesta inflamatoria23.
Este conjunto de enfermedades derivadas de una mala calidad del aire se ha relacionado con la osteoporosis, reflejando una alta relación a nivel genético entre los ambos tipos de enfermedades. Los resultados del presente estudio son consistentes con estudios previos sobre asociación de enfermedades respiratorias con osteoporosis12-14. Además, se han obtenido genes especialmente relevantes que podrían considerarse como potenciales marcadores comunes entre enfermedades pulmonares y osteoporóticas.
En la biología de sistemas, el descubrimiento de las principales proteínas y sus correspondientes vías moleculares en enfermedades complejas está en auge gracias a los análisis de redes PPI y análisis de enriquecimiento. El enfoque de biología de sistemas para investigar la biología asociada a enfermedades está revolucionando la comprensión de vías celulares y redes de genes que subyacen a la aparición de enfermedades, facilitando de esta manera el descubrimiento y desarrollo de nuevos fármacos y terapias gracias a la identificación de marcadores en la progresión de la enfermedad.
En redes biomoleculares, donde los genes y/o las proteínas son nodos y las interacciones moleculares son los bordes o aristas que interconectan los nodos, la importancia de un nodo se puede medir por el efecto producido en los cambios de función de la red de interacción tras la eliminación de ese nodo. Estos nodos esenciales son los llamados hubs. En este estudio se han determinado 5 hubs que podrían perturbar de manera consistente la red de interacción y eso podría ocasionar efectos importantes en la aptitud del organismo. Los hubs que han sido seleccionados y que relacionan enfermedades pulmonares de origen ambiental u ocupacional con osteoporosis son AKT1, ALB, IL6, TP53 y VEGFA, tabla 4. Sin embargo, estudios futuros podrían sugerir nuevos hubs que presentan un score más bajo en la clasificación realizada en este estudio.
El papel que desempeñan estos hubs en los procesos osteoporóticos y en las enfermedades pulmonares es bastante conocido y corrobora la importancia de los genes seleccionados en este trabajo para relacionar las enfermedades objeto de estudio.
En el caso del gen AKT1, que codifica la proteína denominada quinasa AKT1, se encuentra en diversos tipos de células desempeñando funciones importantes en muchas vías de señalización; en el caso de la osteoporosis, AKT1 actúa como regulador negativo de la diferenciación de los osteoblastos y como mediador positivo de la osteoclastogénesis. Se considera por tanto, un potencial objetivo como diana terapéutica para mejorar la diferenciación de osteoblastos y formación de masa ósea limitando al mismo tiempo el desarrollo de los osteoclastos y la resorción ósea24,25. Esta proteína también tiene implicación a nivel pulmonar, participando en la aparición y progresión de fibrosis pulmonar, al promover la diferenciación de los miofibroblastos y la deposición de proteínas en el espacio extracelular y por otra parte, a través de la regulación de TSP1. Debido a esto, AKT1 ha sido propuesta como biomarcador de la enfermedad pulmonar intersticial26, de manera que el uso de inhibidores de AKT1 como triciribine, se plantea como potencial estrategia terapéutica para el tratamiento de estas enfermedades. Estos datos corroboran su papel dual en ambas patologías.
Por otra parte, se identificó el gen ALB que codifica la albúmina. Los bajos niveles de esta proteína se han asociado a un mayor riesgo de fracturas debido a una baja densidad mineral ósea (DMO), existiendo una correlación positiva entre la albúmina sérica y la DMO27,28. En relación a enfermedades pulmonares, se ha mostrado una asociación entre albúmina y enfermedad pulmonar intersticial, observándose niveles significativamente inferiores de albúmina en pacientes afectados con enfermedad pulmonar intersticial29.
La interleucina 6 codificada por el gen IL6, se considera un factor clave de la osteoporosis postmenopáusica debido a su capacidad para activar los osteoclastos e inducir la reabsorción ósea30. Además, ha sido identificada como objetivo prometedor para el tratamiento osteoporótico debido a su papel fundamental como inhibidor de la osteogénesis y como predictor de la pérdida de masa ósea31,32. En relación a enfermedades pulmonares, podemos destacar estudios que correlacionan la sobreexpresión de IL6 con EPOC33,34, aunque existen otros que apuntan en sentido contrario35. Con respecto a silicosis, existen varios estudios que muestran niveles aumentados de esta citoquina desde estados tempranos de la enfermedad siendo considerado como un potencial biomarcador para el diagnóstico precoz y el tratamiento de los pacientes con silicosis36. Así mismo, existen estudios que relacionan el aumento de niveles de IL6 con el desarrollo de antracosis. En este sentido, varios estudios muestran una asociación de ciertos polimorfismos de IL6 con el desarrollo de antracosis y/o antrosilicosis así como con efectos genotóxicos asociados37,38.
Con respecto a la proteína tumoral p53, codificada por el gen TP53, existen varios estudios que la relacionan tanto con osteoporosis como con distintas enfermedades pulmonares. En este contexto, varios estudios muestran que los niveles séricos aumentados de p53 se asocian con una disminución de masa ósea, y la supresión de p53 revierte parcialmente la disminución de la DMO in vitro e in vivo39. En el caso de enfermedades pulmonares, juega un papel protagonista en el desarrollo de EPOC, observándose niveles incrementados de esta proteína en pacientes afectados40. Así mismo, se ha observado una asociación entre p53 y mesotelioma41. En este sentido, se ha reportado que aproximadamente el 15% de pacientes afectados con mesotelioma presenta algún tipo de mutación en el gen TP53. Del mismo modo, se ha observado que aproximadamente en el 25% de mesoteliomas humanos existen co-deleciones de los genes TP53 junto con PTEN y/o CDKN2A/p14ARF, asociándose las pérdidas cooperativas de estos genes con el desarrollo de una proporción significativa de estas neoplasias agresivas42. Estos hallazgos pueden dar lugar al establecimiento de estrategias terapéuticas adecuadas para el tratamiento conjunto de ambos tipos de enfermedades.
Por último, el gen VEGFA que codifica el factor de crecimiento endotelial vascular A se sabe que tiene un papel importante en la biología del hueso participando en la osificación endocondral. Además, VEGFA se expresa a niveles elevados en los precursores de los osteoblastos y puede estimular la diferenciación osteogénica en varios tipos de células43. Por otra parte, se ha observado que niveles disminuidos de VEGFA se asocian con osteoporosis mediante un mecanismo intracelular mediado posiblemente por la regulación del factor de transcripción RUNX2, habiéndose llegado a considerar como objetivo terapéutico de la enfermedad44. Este gen también ha mostrado una actividad relevante en enfermedades pulmonares, siendo un factor de crecimiento autocrino en el mesotelioma y un potente mitógeno para las células mesoteliales45. VEGFA está implicado en la angiogénesis y estimula la neovascularización de los tumores. En el mesotelioma maligno se han detectado niveles elevados de VEGFA y de su receptor mediante inmunohistoquímica46,47 y se han correlacionado con la densidad de microvasos, el aumento de la necrosis tumoral y una peor supervivencia48. En este sentido, fármacos dirigidos al bloqueo de VEGFA, como el bevacizumab, han demostrado recientemente su eficacia en el tratamiento del mesotelioma49.
A pesar de la cantidad de información del papel de estos genes en determinadas enfermedades tanto pulmonares como osteoporóticas, hay pocos estudios que relacionen ambos tipos de enfermedades siendo necesario, por tanto, profundizar en estos aspectos con objeto de identificar posibles dianas terapéuticas comunes. La asociación entre enfermedades pulmonares (desencadenadas por contaminación atmosférica o por una mala calidad del aire en entornos laborales) y calidad ósea podría abrir la puerta al diseño de estrategias conductuales orientadas a la modificación del estilo de vida que implique una mejoría en la salud ósea. Por otra parte, es importante considerar a la población con algún tipo de enfermedad pulmonar, como población de riesgo de fragilidad ósea. Por todo ello, el conocimiento que aporta este tipo de estudios revela la importancia de estudiar la DMO en personas con enfermedades pulmonares para establecer medidas terapéuticas y preventivas tempranas con el objeto de reducir el riesgo de fracturas en esta población vulnerable.