Previamente, informamos que los niveles séricos del metabolito intestinal del triptófano, ácido indol-3-propiónico (IPA), son más bajos en pacientes con fibrosis hepática. En este estudio, investigamos el transcriptoma y el metiloma del ADN en hígados obesos en relación con los niveles séricos de IPA, así como su papel en la inducción de la inactivación fenotípica de las células estrelladas hepáticas (HSC) in vitro.
El estudio incluyó a 116 pacientes obesos sin diabetes mellitus tipo 2 (DM2) (edad: 46,8 ± 9,3 años; IMC: 42,7 ± 5,0 kg/m²) sometidos a cirugía bariátrica en el Centro de Cirugía Bariátrica de Kuopio (KOBS). Se midieron los niveles circulantes de IPA mediante cromatografía líquida-espectrometría de masas (LC-MS), se realizó un análisis del transcriptoma hepático mediante secuenciación de ARN total y se realizó un análisis de metilación del ADN con el chip Infinium HumanMethylation450 BeadChip. Se utilizaron células estrelladas hepáticas humanas (LX-2) para los experimentos in vitro.
Los niveles séricos de IPA se correlacionaron con la expresión de genes implicados en las vías apoptósicas, mitofágicas y de longevidad en el hígado. El gen AKT serina/treonina quinasa 1 (AKT1) fue el gen interactuante más abundante y dominante en los perfiles de transcripción hepática y metilación del ADN. El tratamiento con IPA indujo apoptosis, disminuyó la respiración mitocondrial y alteró la morfología celular y la dinámica mitocondrial al modular la expresión de genes que regulan la fibrosis, la apoptosis y la supervivencia de las células LX-2.
En conjunto, estos datos respaldan que el IPA tiene efectos terapéuticos potenciales y puede inducir la apoptosis y cambiar el fenotipo de las HSC hacia un estado inactivo, ampliando así la posibilidad de inhibir la fibrosis hepática al interferir con la activación de las HSC y el metabolismo mitocondrial.
La prevalencia de la obesidad y el síndrome metabólico se ha asociado con una incidencia creciente de la enfermedad del hígado graso metabólicamente asociada (MASLD); la enfermedad afecta del 25% al ​​30% de la población general [1]. La principal consecuencia de la etiología de la MASLD es la fibrosis hepática, un proceso dinámico caracterizado por la acumulación continua de matriz extracelular fibrosa (ECM) [2]. Las principales células implicadas en la fibrosis hepática son las células estrelladas hepáticas (HSC), que presentan cuatro fenotipos conocidos: quiescente, activado, inactivado y senescente [3, 4]. Las HSC pueden activarse y transdiferenciarse de una forma quiescente a células proliferativas similares a fibroblastos con altas demandas de energía, con mayor expresión de actina de músculo liso α (α-SMA) y colágeno tipo I (Col-I) [5, 6]. Durante la reversión de la fibrosis hepática, las HSC activadas se eliminan por apoptosis o inactivación. Estos procesos incluyen la regulación negativa de los genes fibróticos y la modulación de los genes prosupervivencia (como las vías de señalización NF-κB y PI3K/Akt) [7, 8], así como cambios en la dinámica y la función mitocondrial [9].
Se ha descubierto que los niveles séricos del metabolito del triptófano, ácido indol-3-propiónico (IPA), producido en el intestino, disminuyen en enfermedades metabólicas humanas, incluyendo la MASLD [10–13]. El IPA se asocia con la ingesta de fibra dietética, es conocido por sus efectos antioxidantes y antiinflamatorios, y atenúa el fenotipo de esteatohepatitis no alcohólica (NASH) inducido por la dieta in vivo e in vitro [11–14]. Parte de la evidencia proviene de nuestro estudio previo, que demostró que los niveles séricos de IPA fueron menores en pacientes con fibrosis hepática que en pacientes obesos sin fibrosis hepática en el Estudio de Cirugía Bariátrica de Kuopio (KOBS). Además, demostramos que el tratamiento con IPA podría reducir la expresión de genes que son marcadores clásicos de la adhesión celular, la migración celular y la activación de células madre hematopoyéticas en un modelo de células estrelladas hepáticas humanas (LX-2) y es un posible metabolito hepatoprotector [15]. Sin embargo, aún no está claro cómo el IPA induce la regresión de la fibrosis hepática activando la apoptosis de las HSC y la bioenergética mitocondrial.
En este estudio, demostramos que el IPA sérico se asocia con la expresión de genes enriquecidos en las vías de apoptosis, mitofagia y longevidad en el hígado de personas obesas pero sin diabetes tipo 2 (KOBS). Además, descubrimos que el IPA puede inducir la depuración y degradación de células madre hematopoyéticas (CMH) activadas mediante la vía de inactivación. Estos resultados revelan un nuevo papel para el IPA, convirtiéndolo en una posible diana terapéutica para promover la regresión de la fibrosis hepática.
Un estudio previo en la cohorte KOBS mostró que los pacientes con fibrosis hepática tenían niveles circulantes de IPA más bajos en comparación con los pacientes sin fibrosis hepática [15]. Para excluir el posible efecto de confusión de la diabetes tipo 2, reclutamos a 116 pacientes obesos sin diabetes tipo 2 (edad media ± DE: 46,8 ± 9,3 años; IMC: 42,7 ± 5,0 kg/m²) (Tabla 1) del estudio KOBS en curso como población de estudio [16]. Todos los participantes dieron su consentimiento informado por escrito y el protocolo del estudio fue aprobado por el Comité de Ética del Hospital del Condado de North Savo de acuerdo con la Declaración de Helsinki (54/2005, 104/2008 y 27/2010).
Las muestras de biopsia hepática se obtuvieron durante la cirugía bariátrica y fueron evaluadas histológicamente por patólogos experimentados según los criterios descritos previamente [17, 18]. Los criterios de evaluación se resumen en la Tabla Suplementaria S1 y se describieron previamente [19].
Se analizaron muestras de suero en ayunas mediante cromatografía líquida-espectrometría de masas (LC-MS) para análisis metabolómico (n = 116). Las muestras se analizaron con un sistema UHPLC-qTOF-MS (LC 1290, qTOF-MS 6540, Agilent Technologies, Waldbronn, Karlsruhe, Alemania), como se describió previamente19. La identificación del alcohol isopropílico (AIP) se basó en el tiempo de retención y la comparación del espectro MS/MS con estándares puros. La intensidad de la señal del AIP (área del pico) se consideró en todos los análisis posteriores [20].
La secuenciación completa del ARN del hígado se realizó con Illumina HiSeq 2500 y los datos se preprocesaron como se describió previamente [19, 21, 22]. Realizamos un análisis de expresión diferencial dirigido de las transcripciones que afectan la función/biogénesis mitocondrial utilizando 1957 genes seleccionados de la base de datos MitoMiner 4.0 [23]. El análisis de la metilación del ADN hepático se realizó con Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA, EE. UU.) utilizando la misma metodología que se describió previamente [24, 25].
Las células estrelladas hepáticas humanas (LX-2) fueron amablemente proporcionadas por el Prof. Stefano Romeo y se cultivaron y mantuvieron en medio DMEM/F12 (Biowest, L0093-500, 1% Pen/Strep; Lonza, DE17-602E, 2% FBS; Gibco, 10270-106). Para seleccionar la dosis de trabajo de IPA, las células LX-2 se trataron con diferentes concentraciones de IPA (10 μM, 100 μM y 1 mM; Sigma, 220027) en medio DMEM/F12 durante 24 h. Además, para investigar la capacidad del IPA para inactivar las HSC, las células LX-2 se cotrataron con 5 ng/ml de TGF-β1 (R&D systems, 240-B-002/CF) y 1 mM de IPA en medio sin suero durante 24 h. Para los controles del vehículo correspondientes, se utilizó 4 nM HCL que contenía 0,1 % de BSA para el tratamiento con TGF-β1 y se utilizó 0,05 % DMSO para el tratamiento con IPA, y ambos se utilizaron juntos para el tratamiento combinado.
La apoptosis se evaluó utilizando el kit de detección de apoptosis con anexina V FITC con 7-AAD (Biolegend, San Diego, CA, EE. UU., n.° de cat. 640922) según las instrucciones del fabricante. En resumen, se cultivaron LX-2 (1 × 10⁻¹ células/pocillo) durante la noche en placas de 12 pocillos y luego se trataron con múltiples dosis de IPA o IPA y TGF-β1. Al día siguiente, se recogieron las células flotantes y adherentes, se tripsinizaron, se lavaron con PBS, se resuspendieron en tampón de unión a anexina V y se incubaron con FITC-anexina V y 7-AAD durante 15 min.
Las mitocondrias de células vivas se tiñeron para determinar su actividad oxidativa con Mitotracker™ Red CMXRos (MTR) (Thermo Fisher Scientific, Carlsbad, CA). Para los ensayos de MTR, las células LX-2 se incubaron a densidades iguales con IPA y TGF-β1. Después de 24 h, las células vivas se tripsinizaron, se lavaron con PBS y se incubaron con 100 μM de MTR en medio sin suero a 37 °C durante 20 min, como se describió previamente [26]. Para el análisis de la morfología de las células vivas, se analizaron el tamaño celular y la complejidad citoplasmática mediante los parámetros de dispersión frontal (FSC) y dispersión lateral (SSC), respectivamente.
Todos los datos (30.000 eventos) se adquirieron utilizando NovoCyte Quanteon (Agilent) y se analizaron utilizando el software NovoExpress® 1.4.1 o FlowJo V.10.
La tasa de consumo de oxígeno (OCR) y la tasa de acidificación extracelular (ECAR) se midieron en tiempo real utilizando un analizador de flujo extracelular Seahorse (Agilent Technologies, Santa Clara, CA) equipado con un sistema de estrés celular Seahorse XF Mito, según las instrucciones del fabricante. Brevemente, se sembraron 2 × 10⁻¹ células LX-2 por pocillo en placas de cultivo celular XF96. Tras la incubación durante la noche, las células se trataron con isopropanol (IPA) y TGF-β1 (Métodos complementarios 1). El análisis de datos se realizó con el software Seahorse XF Wave, que incluye el generador de informes de prueba de fenotipo de energía celular Seahorse XF. A partir de esto, se calculó el índice de salud bioenergética (BHI) [27].
El ARN total se transcribió a ADNc. Para métodos específicos, véase la referencia [15]. Se utilizaron los niveles de ARNm de la proteína ácida ribosomal 60S P0 (RPLP0) y la ciclofilina A1 (PPIA) humanas como controles génicos constitutivos. Se utilizó el sistema de PCR en tiempo real QuantStudio 6 pro (Thermo Fisher, Landsmeer, Países Bajos) con el kit TaqMan™ Fast Advanced Master Mix (Applied Biosystems) o el kit Sensifast SYBR Lo-ROX (Bioline, BIO 94050). Se calculó el pliegue de expresión génica relativa utilizando los parámetros comparativos de ciclado del valor Ct (ΔΔCt) y el método ∆∆Ct. Los detalles de los cebadores se proporcionan en las Tablas Suplementarias S2 y S3.
El ADN nuclear (ADNnc) y el ADN mitocondrial (ADNmt) se extrajeron utilizando el kit de sangre y tejido DNeasy (Qiagen), como se describió previamente [28]. La cantidad relativa de ADNmt se calculó calculando la relación entre cada región diana de ADNmt y la media geométrica de las tres regiones de ADN nuclear (ADNmt/ADNnc), como se detalla en los Métodos Suplementarios 2. Los detalles de los cebadores para el ADNmt y el ADNnc se proporcionan en la Tabla Suplementaria S4.
Las células vivas se tiñeron con Mitotracker™ Red CMXRos (MTR) (Thermo Fisher Scientific, Carlsbad, CA) para visualizar las redes mitocondriales intercelulares e intracelulares. Las células LX-2 (1 × 10⁻¹ células/pocillo) se cultivaron en portaobjetos de vidrio en las placas de cultivo con fondo de vidrio correspondientes (Ibidi GmbH, Martinsried, Alemania). Después de 24 h, las células LX-2 vivas se incubaron con 100 μM de MTR durante 20 min a 37 °C y los núcleos celulares se tiñeron con DAPI (1 μg/ml, Sigma-Aldrich) como se describió previamente [29]. Las redes mitocondriales se visualizaron con un microscopio invertido Zeiss Axio Observer (Carl Zeiss Microimaging GmbH, Jena, Alemania) equipado con un módulo confocal Zeiss LSM 800 a 37 °C en una atmósfera humidificada con un 5 % de CO2 y un objetivo de 63×NA 1.3. Se adquirieron diez imágenes de la serie Z para cada tipo de muestra. Cada serie Z contiene 30 secciones, cada una con un espesor de 9,86 μm. Para cada muestra, se adquirieron imágenes de diez campos de visión diferentes con el software ZEN 2009 (Carl Zeiss Microimaging GmbH, Jena, Alemania), y el análisis de la morfología mitocondrial se realizó con el software ImageJ (v1.54d) [30, 31], de acuerdo con los parámetros detallados en el Método Suplementario 3.
Las células se fijaron con glutaraldehído al 2% en tampón fosfato 0,1 M, seguidas de una fijación con una solución de tetróxido de osmio al 1% (Sigma Aldrich, MO, EE. UU.), se deshidrataron gradualmente con acetona (Merck, Darmstadt, Alemania) y finalmente se incluyeron en resina epoxi. Se prepararon secciones ultrafinas y se tiñeron con acetato de uranilo al 1% (Merck, Darmstadt, Alemania) y citrato de plomo al 1% (Merck, Darmstadt, Alemania). Se obtuvieron imágenes ultraestructurales con un microscopio electrónico de transmisión JEM 2100F EXII (JEOL Ltd, Tokio, Japón) a un voltaje de aceleración de 80 kV.
La morfología de las células LX-2 tratadas con IPA durante 24 h se analizó mediante microscopía de contraste de fases con un aumento de 50x utilizando un microscopio óptico invertido Zeiss (Zeiss Axio Vert.A1 y AxioCam MRm, Jena, Alemania).
Los datos clínicos se expresaron como media ± desviación estándar o mediana (rango intercuartil: RIQ). Se utilizó un análisis de varianza unidireccional (variables continuas) o la prueba χ² (variables categóricas) para comparar las diferencias entre los tres grupos de estudio. Se utilizó la tasa de falsos positivos (TFP) para corregir las pruebas múltiples, y los genes con una TFP < 0,05 se consideraron estadísticamente significativos. Se utilizó el análisis de correlación de Spearman para correlacionar la metilación del ADN CpG con la intensidad de la señal del IPA, con valores de p nominales (p < 0,05).
Se realizó un análisis de vías mediante una herramienta web de análisis de conjuntos de genes (WebGestalt) para 268 transcripciones (p nominal < 0,01), 119 transcripciones asociadas a mitocondrias (p nominal < 0,05) y 4350 CpG de 3093 transcripciones hepáticas asociadas con los niveles séricos circulantes de IPA. Se utilizó la herramienta gratuita Venny DB (versión 2.1.0) para identificar genes superpuestos, y StringDB (versión 11.5) para visualizar las interacciones proteína-proteína.
Para el experimento LX-2, se analizó la normalidad de las muestras mediante la prueba de D'Agostino-Pearson. Los datos se obtuvieron de al menos tres réplicas biológicas y se sometieron a un ANOVA de una vía con la prueba post hoc de Bonferroni. Un valor p inferior a 0,05 se consideró estadísticamente significativo. Los datos se presentan como media ± DE, y el número de experimentos se indica en cada figura. Todos los análisis y gráficos se realizaron con el software estadístico GraphPad Prism 8 para Windows (GraphPad Software Inc., versión 8.4.3, San Diego, EE. UU.).
En primer lugar, investigamos la asociación de los niveles séricos de IPA con los transcritos hepáticos, corporales y mitocondriales. En el perfil de transcritos totales, el gen con mayor asociación con los niveles séricos de IPA fue MAPKAPK3 (FDR = 0,0077; proteína quinasa activada por mitógenos-proteína quinasa 3 activada); en el perfil de transcritos relacionado con las mitocondrias, el gen con mayor asociación fue AKT1 (FDR = 0,7621; AKT serina/treonina quinasa 1) (Archivo adicional 1 y Archivo adicional 2).
A continuación, analizamos las transcripciones globales (n = 268; p < 0,01) y las transcripciones asociadas a las mitocondrias (n = 119; p < 0,05), identificando finalmente la apoptosis como la vía canónica más significativa (p = 0,0089). Para las transcripciones mitocondriales asociadas con los niveles séricos de IPA, nos centramos en la apoptosis (FDR = 0,00001), la mitofagia (FDR = 0,00029) y las vías de señalización del TNF (FDR = 0,000006) (Figura 1A, Tabla 2 y Figuras suplementarias 1A-B).
Análisis superpuesto de transcripciones globales, asociadas a mitocondrias y metilación de ADN en hígado humano en asociación con niveles séricos de IPA. A representa 268 transcripciones globales, 119 transcripciones asociadas a mitocondrias y transcripciones de ADN metilado que se asignan a 3092 sitios CpG asociados con niveles séricos de IPA (valores p < 0,01 para transcripciones globales y ADN metilado, y valores p < 0,05 para transcripciones mitocondriales). Las principales transcripciones superpuestas se muestran en el medio (AKT1 y YKT6). B El mapa de interacción de los 13 genes con la puntuación de interacción más alta (0,900) con otros genes se construyó a partir de los 56 genes superpuestos (región de línea negra) que se asociaron significativamente con niveles séricos de IPA utilizando la herramienta en línea StringDB. Verde: Genes asignados al componente celular de Gene Ontology (GO): mitocondrias (GO:0005739). AKT1 es la proteína con la puntuación más alta (0,900) en interacciones con otras proteínas, según los datos (basados ​​en minería de texto, experimentos, bases de datos y coexpresión). Los nodos de la red representan proteínas y los bordes representan conexiones entre proteínas.
Dado que los metabolitos de la microbiota intestinal pueden regular la composición epigenética mediante la metilación del ADN [32], investigamos si los niveles séricos de IPA se asociaban con la metilación del ADN hepático. Descubrimos que los dos principales sitios de metilación asociados con los niveles séricos de IPA se encontraban cerca de la región 3 rica en prolina-serina (C19orf55) y del miembro 6 (pequeño) de la familia de proteínas de choque térmico B (HSPB6) (Archivo adicional 3). La metilación del ADN de 4350 CpG (p < 0,01) se correlacionó con los niveles séricos de IPA y se enriqueció en vías reguladoras de la longevidad (p = 0,006) (Figura 1A, Tabla 2 y Figura Suplementaria 1C).
Para comprender los mecanismos biológicos subyacentes a las asociaciones entre los niveles séricos de IPA, las transcripciones globales, las transcripciones asociadas a las mitocondrias y la metilación del ADN en el hígado humano, realizamos un análisis de superposición de los genes identificados en el análisis de la vía anterior (Figura 1A). Los resultados del análisis de enriquecimiento de la vía de los 56 genes superpuestos (dentro de la línea negra en la Figura 1A) mostraron que la vía de la apoptosis (p = 0,00029) destacó dos genes comunes a los tres análisis: AKT1 y YKT6 (homólogo de YKT6 v-SNARE), como se muestra en el diagrama de Venn (Figura Suplementaria 2 y Figura 1A). Curiosamente, encontramos que AKT1 (cg19831386) y YKT6 (cg24161647) se correlacionaron positivamente con los niveles séricos de IPA (Archivo adicional 3). Para identificar posibles interacciones proteicas entre productos génicos, seleccionamos 13 genes con la mayor puntuación de región común (0,900) entre 56 genes superpuestos como entrada y construimos un mapa de interacción. Según el nivel de confianza (confianza marginal), el gen AKT1 con la mayor puntuación (0,900) se ubicó en la cima (Figura 1B).
Con base en el análisis de la vía, encontramos que la apoptosis era la vía principal, por lo que investigamos si el tratamiento con IPA afectaría la apoptosis de las HSC in vitro. Previamente demostramos que diferentes dosis de IPA (10 μM, 100 μM y 1 mM) no eran tóxicas para las células LX-2 [15]. Este estudio mostró que el tratamiento con IPA a 10 μM y 100 μM aumentó el número de células viables y necróticas. Sin embargo, en comparación con el grupo de control, la viabilidad celular disminuyó a una concentración de IPA de 1 mM, mientras que la tasa de necrosis celular se mantuvo sin cambios (Figura 2A, B). A continuación, para encontrar la concentración óptima para inducir la apoptosis en las células LX-2, probamos IPA a 10 μM, 100 μM y 1 mM durante 24 h (Figura 2A-E y Figura Suplementaria 3A-B). Curiosamente, IPA 10 μM y 100 μM disminuyeron la tasa de apoptosis (%), sin embargo, IPA 1 mM aumentó la apoptosis tardía y la tasa de apoptosis (%) en comparación con el control y, por lo tanto, fue elegido para experimentos adicionales (Figuras 2A–D).
El IPA induce la apoptosis de las células LX-2. Se empleó el método de doble tinción con anexina V y 7-AAD para cuantificar la tasa apoptótica y la morfología celular mediante citometría de flujo. Las células BA se incubaron con IPA 10 μM, 100 μM y 1 mM durante 24 h o con F–H TGF-β1 (5 ng/ml) y IPA 1 mM en medio sin suero durante 24 h. A: células vivas (anexina V -/ 7AAD-); B: células necróticas (anexina V -/ 7AAD+); C, F: tempranas (anexina V +/ 7AAD-); D, G: tardías (anexina V+/7AAD.+); E, H: porcentaje del total de células apoptóticas tempranas y tardías en la tasa apoptótica (%). Los datos se expresan como media ± DE, n = 3 experimentos independientes. Las comparaciones estadísticas se realizaron mediante ANOVA unidireccional con prueba post hoc de Bonferroni. *p < 0,05; ****p < 0,0001
Como se ha demostrado previamente, 5 ng/ml de TGF-β1 pueden inducir la activación de las células madre hematopoyéticas (HSC) al aumentar la expresión de genes marcadores clásicos [15]. Las células LX-2 se trataron con 5 ng/ml de TGF-β1 y 1 mM de IPA en combinación (Fig. 2E-H). El tratamiento con TGF-β1 no modificó la tasa de apoptosis; sin embargo, el tratamiento conjunto con IPA aumentó la apoptosis tardía y la tasa de apoptosis (%) en comparación con el tratamiento con TGF-β1 (Fig. 2E-H). Estos resultados indican que 1 mM de IPA puede promover la apoptosis en células LX-2 independientemente de la inducción de TGF-β1.
Investigamos más a fondo el efecto del IPA sobre la respiración mitocondrial en células LX-2. Los resultados mostraron que 1 mM de IPA redujo los parámetros de la tasa de consumo de oxígeno (OCR): respiración no mitocondrial, respiración basal y máxima, fuga de protones y producción de ATP, en comparación con el grupo control (Figuras 3A, B), mientras que el índice de salud bioenergética (BHI) no se modificó.
El IPA reduce la respiración mitocondrial en células LX-2. La curva de respiración mitocondrial (OCR) se presenta como parámetros de respiración mitocondrial (respiración no mitocondrial, respiración basal, respiración máxima, fuga de protones, generación de ATP, SRC y BHI). Las células A y B se incubaron con 10 μM, 100 μM y 1 mM de IPA durante 24 h, respectivamente. Las células C y D se incubaron con TGF-β1 (5 ng/ml) y 1 mM de IPA en medio sin suero durante 24 h, respectivamente. Todas las mediciones se normalizaron al contenido de ADN utilizando el kit CyQuant. BHI: índice de salud bioenergética; SRC: capacidad de reserva respiratoria; OCR: tasa de consumo de oxígeno. Los datos se presentan como media ± desviación estándar (DE), n = 5 experimentos independientes. Las comparaciones estadísticas se realizaron mediante ANOVA de una vía y la prueba post hoc de Bonferroni. *p < 0,05; **p < 0,01; y ***p < 0,001
Para comprender mejor el efecto del IPA en el perfil bioenergético de las células LX-2 activadas con TGF-β1, analizamos la fosforilación oxidativa mitocondrial mediante OCR (Fig. 3C, D). Los resultados mostraron que el tratamiento con TGF-β1 pudo reducir la respiración máxima, la capacidad de reserva respiratoria (SRC) y el BHI en comparación con el grupo control (Fig. 3C, D). Además, el tratamiento combinado disminuyó la respiración basal, la fuga de protones y la producción de ATP, pero la SRC y el BHI fueron significativamente mayores que en el grupo tratado con TGF-β1 (Fig. 3C, D).
También realizamos la "Prueba de fenotipo de energía celular" proporcionada por el software Seahorse (Fig. suplementaria 4A–D). Como se muestra en la Fig. suplementaria 3B, los potenciales metabólicos de OCR y ECAR disminuyeron después del tratamiento con TGF-β1; sin embargo, no se observaron diferencias en los grupos de tratamiento combinado y con IPA en comparación con el grupo control. Además, los niveles basales y de estrés de OCR disminuyeron después del tratamiento combinado y con IPA en comparación con el grupo control (Fig. suplementaria 4C). Curiosamente, se observó un patrón similar con la terapia combinada, donde no se observaron cambios en los niveles basales y de estrés de ECAR en comparación con el tratamiento con TGF-β1 (Fig. suplementaria 4C). En las HSC, la reducción de la fosforilación oxidativa mitocondrial y la capacidad del tratamiento combinado para restaurar SCR y BHI después de la exposición al tratamiento con TGF-β1 no alteraron el potencial metabólico (OCR y ECAR). En conjunto, estos resultados indican que el IPA puede reducir la bioenergética en las HSC, lo que sugiere que el IPA puede inducir un perfil energético más bajo que desplaza el fenotipo de las HSC hacia la inactivación (Figura complementaria 4D).
El efecto del IPA en la dinámica mitocondrial se investigó mediante la cuantificación tridimensional de la morfología mitocondrial y las conexiones de red, así como la tinción de MTR (Figura 4 y Figura suplementaria 5). La Figura 4 muestra que, en comparación con el grupo de control, el tratamiento con TGF-β1 disminuyó el área superficial media, el número de ramificaciones, la longitud total de las ramificaciones y el número de uniones de las ramificaciones (Figuras 4A y B) y cambió la proporción de mitocondrias de morfología esférica a intermedia (Figura 4C). Solo el tratamiento con IPA disminuyó el volumen mitocondrial medio y cambió la proporción de mitocondrias de morfología esférica a intermedia en comparación con el grupo de control (Figura 4A). Por el contrario, la esfericidad, la longitud media de las ramificaciones y la actividad mitocondrial evaluada mediante MTR dependiente del potencial de membrana mitocondrial (Figuras 4A y E) se mantuvieron sin cambios y estos parámetros no difirieron entre los grupos. En conjunto, estos resultados sugieren que el tratamiento con TGF-β1 e IPA parece modular la forma y el tamaño mitocondrial, así como la complejidad de la red en las células LX-2 vivas.
El IPA altera la dinámica mitocondrial y la abundancia de ADN mitocondrial en células LX-2. A. Imágenes confocales representativas de células LX-2 vivas incubadas con TGF-β1 (5 ng/ml) y 1 mM de IPA durante 24 h en medio sin suero que muestran redes mitocondriales teñidas con Mitotracker™ Red CMXRos y núcleos teñidos de azul con DAPI. Todos los datos contenían al menos 15 imágenes por grupo. Adquirimos 10 imágenes de pila Z para cada tipo de muestra. Cada secuencia del eje Z contenía 30 cortes, cada uno con un grosor de 9,86 μm. Barra de escala: 10 μm. B. Objetos representativos (solo mitocondrias) identificados aplicando umbral adaptativo a la imagen. Se realizó un análisis cuantitativo y una comparación de las conexiones de la red morfológica mitocondrial para todas las células de cada grupo. C. Frecuencia de las proporciones de forma mitocondrial. Los valores cercanos a 0 indican formas esféricas y los valores cercanos a 1 indican formas filamentosas. D El contenido de ADN mitocondrial (ADNmt) se determinó como se describe en Materiales y Métodos. E El análisis de CMXRos con Mitotracker™ Red se realizó mediante citometría de flujo (30.000 eventos) como se describe en Materiales y Métodos. Los datos se presentan como media ± DE, n = 3 experimentos independientes. Las comparaciones estadísticas se realizaron mediante ANOVA de una vía y la prueba post hoc de Bonferroni. *p < 0,05; **p < 0,01; ***p < 0,001; ****p < 0,0001
A continuación, analizamos el contenido de ADNmt en las células LX-2 como indicador del número de mitocondrias. En comparación con el grupo control, el contenido de ADNmt aumentó en el grupo tratado con TGF-β1 (Figura 4D). En comparación con el grupo tratado con TGF-β1, el contenido de ADNmt disminuyó en el grupo de tratamiento combinado (Figura 4D), lo que sugiere que el IPA podría reducir el contenido de ADNmt y posiblemente el número de mitocondrias, así como la respiración mitocondrial (Figura 3C). Además, el IPA pareció reducir el contenido de ADNmt en el tratamiento combinado, pero no afectó la actividad mitocondrial mediada por MTR (Figuras 4A-C).
Investigamos la asociación del IPA con los niveles de ARNm de genes asociados con la fibrosis, la apoptosis, la supervivencia y la dinámica mitocondrial en células LX-2 (Figuras 5A-D). En comparación con el grupo control, el grupo tratado con TGF-β1 mostró una mayor expresión de genes como la cadena α2 del colágeno tipo I (COL1A2), la actina α de músculo liso (αSMA), la metaloproteinasa de matriz 2 (MMP2), el inhibidor tisular de la metaloproteinasa 1 (TIMP1) y el gen similar a la dinamina 1 (DRP1), lo que indica un aumento de la fibrosis y la activación. Además, en comparación con el grupo control, el tratamiento con TGF-β1 redujo los niveles de ARNm del receptor nuclear de pregnano X (PXR), la caspasa 8 (CASP8), MAPKAPK3, el inhibidor de la α de células B, el potenciador del péptido ligero del gen del factor nuclear κ (NFκB1A) y el inhibidor de la subunidad β de la cinasa del factor nuclear κB (IKBKB) (Figura 5A–D). En comparación con el tratamiento con TGF-β1, el tratamiento combinado con TGF-β1 e IPA redujo la expresión de COL1A2 y MMP2, pero aumentó los niveles de ARNm de PXR, TIMP1, linfoma de células B-2 (BCL-2), CASP8, NFκB1A, NFκB1-β e IKBKB. El tratamiento con IPA disminuyó significativamente la expresión de MMP2, proteína X asociada a Bcl-2 (BAX), AKT1, proteína de atrofia óptica 1 (OPA1) y proteína de fusión mitocondrial 2 (MFN2), mientras que la expresión de CASP8, NFκB1A, NFκB1B e IKBKB aumentó en comparación con el grupo control. Sin embargo, no se encontraron diferencias en la expresión de caspasa-3 (CASP3), factor activador de peptidasa apoptótica 1 (APAF1), proteína de fusión mitocondrial 1 (MFN1) y proteína de fisión 1 (FIS1). En conjunto, estos resultados sugieren que el tratamiento con IPA modula la expresión de genes asociados con la fibrosis, la apoptosis, la supervivencia y la dinámica mitocondrial. Nuestros datos sugieren que el tratamiento con IPA reduce la fibrosis en células LX-2; al mismo tiempo, estimula la supervivencia al cambiar el fenotipo hacia la inactivación.
El IPA modula la expresión de genes responsables de fibroblastos, apoptosis, viabilidad y dinámica mitocondrial en células LX-2. Los histogramas muestran la expresión de ARNm en relación con el control endógeno (RPLP0 o PPIA) tras la inducción de células LX-2 con TGF-β1 e IPA en medio sin suero durante 24 h. A indica fibroblastos, B células apoptóticas, C células supervivientes y D expresión génica responsable de la dinámica mitocondrial. Los datos se presentan como media ± desviación estándar (DE), n = 3 experimentos independientes. Las comparaciones estadísticas se realizaron mediante ANOVA de una vía y la prueba post hoc de Bonferroni. *p < 0,05; **p < 0,01; ***p < 0,001; ****p < 0,0001
Luego, los cambios en el tamaño celular (FSC-H) y la complejidad citoplasmática (SSC-H) se evaluaron mediante citometría de flujo (Figura 6A,B), y los cambios en la morfología celular después del tratamiento con IPA se evaluaron mediante microscopía electrónica de transmisión (TEM) y microscopía de contraste de fase (Figura 6A-B suplementaria). Como se esperaba, las células en el grupo tratado con TGF-β1 aumentaron de tamaño en comparación con el grupo control (Figura 6A,B), mostrando la expansión clásica del retículo endoplasmático rugoso (RE*) y los fagolisosomas (P), lo que indica la activación de las células madre hematopoyéticas (HSC) (Figura 6A suplementaria). Sin embargo, en comparación con el grupo tratado con TGF-β1, el tamaño celular, la complejidad citoplasmática (Fig. 6A,B) y el contenido de RE* disminuyeron en el grupo de tratamiento combinado de TGF-β1 e IPA (Fig. 6A suplementaria). Además, el tratamiento con IPA redujo el tamaño celular, la complejidad citoplasmática (Figs. 6A, B) y el contenido de P y ER* (Fig. 6A suplementaria) en comparación con el grupo control. Asimismo, el contenido de células apoptóticas aumentó tras 24 h de tratamiento con IPA en comparación con el grupo control (flechas blancas, Fig. 6B suplementaria). En conjunto, estos resultados sugieren que 1 mM de IPA puede estimular la apoptosis de las células madre hematopoyéticas (HSC) y revertir los cambios en los parámetros morfológicos celulares inducidos por TGF-β1, regulando así el tamaño y la complejidad celular, lo cual podría estar asociado con la inactivación de las HSC.
El IPA altera el tamaño celular y la complejidad citoplasmática en células LX-2. Imágenes representativas del análisis por citometría de flujo. El análisis empleó una estrategia de selección específica para células LX-2: SSC-A/FSC-A para definir la población celular, FSC-H/FSC-A para identificar dobletes y SSC-H/FSC-H para el análisis del tamaño celular y la complejidad. Las células se incubaron con TGF-β1 (5 ng/ml) y 1 mM de IPA en medio sin suero durante 24 h. Las células LX-2 se distribuyeron en los cuadrantes inferior izquierdo (SSC-H-/FSC-H-), superior izquierdo (SSC-H+/FSC-H-), inferior derecho (SSC-H-/FSC-H+) y superior derecho (SSC-H+/FSC-H+) para el análisis del tamaño celular y la complejidad citoplasmática. B. La morfología celular se analizó mediante citometría de flujo utilizando FSC-H (dispersión frontal, tamaño celular) y SSC-H (dispersión lateral, complejidad citoplasmática) (30.000 eventos). Los datos se presentan como media ± DE, n = 3 experimentos independientes. Las comparaciones estadísticas se realizaron mediante ANOVA de una vía y la prueba post hoc de Bonferroni. *p < 0,05; **p < 0,01; ***p < 0,001 y ****p < 0,0001
Los metabolitos intestinales como el IPA se han convertido en un tema de investigación de gran interés, lo que sugiere que podrían descubrirse nuevas dianas en la microbiota intestinal. Por lo tanto, resulta interesante que el IPA, un metabolito que hemos vinculado a la fibrosis hepática en humanos [15], haya demostrado ser un posible compuesto antifibrótico en modelos animales [13, 14]. En este estudio, demostramos por primera vez una asociación entre el IPA sérico y la transcriptómica hepática global y la metilación del ADN en personas obesas sin diabetes tipo 2 (DT2), destacando la apoptosis, la mitofagia y la longevidad, así como un posible gen candidato, AKT1, que regula la homeostasis hepática. Otra novedad de nuestro estudio es que demostramos la interacción del tratamiento con IPA con la apoptosis, la morfología celular, la bioenergética mitocondrial y la dinámica en células LX-2, lo que indica un espectro energético más bajo que desplaza el fenotipo de las HSC hacia la inactivación, convirtiendo al IPA en un candidato potencial para mejorar la fibrosis hepática.
Descubrimos que la apoptosis, la mitofagia y la longevidad fueron las vías canónicas más importantes enriquecidas en genes hepáticos asociados con el IPA sérico circulante. La interrupción del sistema de control de calidad mitocondrial (MQC) puede conducir a disfunción mitocondrial, mitofagia y apoptosis, promoviendo así la aparición de MASLD[33, 34]. Por lo tanto, podemos especular que el IPA puede estar involucrado en el mantenimiento de la dinámica celular y la integridad mitocondrial a través de la apoptosis, la mitofagia y la longevidad en el hígado. Nuestros datos mostraron que dos genes fueron comunes en los tres ensayos: YKT6 y AKT1. Vale la pena señalar que YKT6 es una proteína SNARE involucrada en el proceso de fusión de la membrana celular. Desempeña un papel en la autofagia y la mitofagia al formar un complejo de iniciación con STX17 y SNAP29 en el autofagosoma, promoviendo así la fusión de autofagosomas y lisosomas[35]. Además, la pérdida de la función de YKT6 resulta en mitofagia deteriorada[36], mientras que la regulación positiva de YKT6 se asocia con la progresión del carcinoma hepatocelular (CHC), mostrando una mayor supervivencia celular[37]. Por otro lado, AKT1 es el gen interactuante más importante y desempeña un papel importante en enfermedades hepáticas, incluyendo la vía de señalización PI3K/AKT, el ciclo celular, la migración celular, la proliferación, la adhesión focal, la función mitocondrial y la secreción de colágeno[38–40]. La vía de señalización PI3K/AKT activada puede activar las células madre hematopoyéticas (HSC), que son las células responsables de la producción de la matriz extracelular (ECM), y su desregulación puede contribuir a la aparición y progresión de la fibrosis hepática[40]. Además, AKT es uno de los factores clave de supervivencia celular que inhibe la apoptosis celular dependiente de p53, y la activación de AKT puede estar asociada con la inhibición de la apoptosis de las células hepáticas[41, 42]. Los resultados obtenidos sugieren que el IPA podría estar involucrado en la apoptosis asociada a las mitocondrias hepáticas, al afectar la decisión de los hepatocitos entre entrar en apoptosis o sobrevivir. Estos efectos podrían estar regulados por los genes candidatos AKT y/o YKT6, cruciales para la homeostasis hepática.
Nuestros resultados mostraron que 1 mM de IPA indujo apoptosis y disminuyó la respiración mitocondrial en células LX-2 independientemente del tratamiento con TGF-β1. Cabe destacar que la apoptosis es una vía principal para la resolución de la fibrosis y la activación de las células madre hematopoyéticas (HSC), y también es un evento clave en la respuesta fisiológica reversible de la fibrosis hepática [4, 43]. Además, la restauración de BHI en células LX-2 después del tratamiento combinado proporcionó nuevos conocimientos sobre el papel potencial del IPA en la regulación de la bioenergética mitocondrial. En condiciones de reposo e inactivas, las células hematopoyéticas normalmente utilizan la fosforilación oxidativa mitocondrial para producir ATP y tienen baja actividad metabólica. Por otro lado, la activación de las HSC mejora la respiración y la biosíntesis mitocondrial para compensar las demandas energéticas de entrar en el estado glucolítico [44]. El hecho de que el IPA no afectara el potencial metabólico y el ECAR sugiere que la vía glucolítica tiene menos prioridad. De manera similar, otro estudio mostró que 1 mM de IPA fue capaz de modular la actividad de la cadena respiratoria mitocondrial en cardiomiocitos, la línea celular de hepatocitos humanos (Huh7) y las células endoteliales de la vena umbilical humana (HUVEC); Sin embargo, no se encontró ningún efecto del IPA en la glucólisis en cardiomiocitos, lo que sugiere que el IPA puede afectar la bioenergética de otros tipos de células [45]. Por lo tanto, especulamos que 1 mM de IPA puede actuar como un desacoplador químico leve, ya que puede reducir significativamente la expresión génica fibrogénica, la morfología celular y la bioenergética mitocondrial sin cambiar la cantidad de ADNmt [46]. Los desacopladores mitocondriales pueden inhibir la fibrosis inducida por cultivo y la activación de HSC [47] y reducir la producción de ATP mitocondrial regulada o inducida por ciertas proteínas como las proteínas desacopladoras (UCP) o la translocasa de nucleótido de adenina (ANT). Dependiendo del tipo de célula, este fenómeno puede proteger a las células de la apoptosis y/o promover la apoptosis [46]. Sin embargo, se necesitan más estudios para dilucidar el papel del IPA como desacoplador mitocondrial en la inactivación de células madre hematopoyéticas.
Luego investigamos si los cambios en la respiración mitocondrial se reflejan en la morfología mitocondrial en células LX-2 vivas. Curiosamente, el tratamiento con TGF-β1 altera la proporción mitocondrial de esférica a intermedia, con una disminución de la ramificación mitocondrial y un aumento de la expresión de DRP1, un factor clave en la fisión mitocondrial [48]. Además, la fragmentación mitocondrial se asocia con la complejidad general de la red, y la transición de la fusión a la fisión es crítica para la activación de las células madre hematopoyéticas (HSC), mientras que la inhibición de la fisión mitocondrial conduce a la apoptosis de las HSC [49]. Por lo tanto, nuestros resultados indican que el tratamiento con TGF-β1 puede inducir una disminución en la complejidad de la red mitocondrial con una menor ramificación, que es más común en la fisión mitocondrial asociada con las células madre hematopoyéticas (HSC) activadas. Además, nuestros datos mostraron que el IPA podría cambiar la proporción de mitocondrias de forma esférica a intermedia, reduciendo así la expresión de OPA1 y MFN2. Estudios han demostrado que la regulación negativa de OPA1 podría causar una disminución del potencial de membrana mitocondrial y desencadenar la apoptosis celular[50]. Se sabe que MFN2 media la fusión mitocondrial y la apoptosis[51]. Los resultados obtenidos sugieren que la inducción de células LX-2 por TGF-β1 y/o IPA parece modular la forma y el tamaño mitocondriales, así como el estado de activación y la complejidad de la red.
Nuestros resultados indican que el tratamiento combinado de TGFβ-1 e IPA podría reducir el ADNmt y los parámetros morfológicos celulares al regular la expresión del ARNm de genes relacionados con la fibrosis, la apoptosis y la supervivencia en células que evaden la apoptosis. De hecho, el IPA disminuyó la expresión del ARNm de AKT1 y de genes importantes de fibrosis como COL1A2 y MMP2, pero aumentó la expresión de CASP8, que se asocia con la apoptosis. Nuestros resultados mostraron que, tras el tratamiento con IPA, la expresión de BAX disminuyó y la expresión del ARNm de las subunidades de la familia TIMP1, BCL-2 y NF-κB aumentó, lo que sugiere que el IPA podría estimular señales de supervivencia en células madre hematopoyéticas (HSC) que evaden la apoptosis. Estas moléculas podrían actuar como señales pro-supervivencia en células madre hematopoyéticas activadas, lo que podría estar asociado con una mayor expresión de proteínas antiapoptóticas (como Bcl-2), una menor expresión de BAX proapoptótica y una compleja interacción entre TIMP y NF-κB [5, 7]. El IPA ejerce sus efectos a través de PXR, y encontramos que el tratamiento combinado con TGF-β1 e IPA aumentó los niveles de expresión de ARNm de PXR, lo que indica la supresión de la activación de las HSC. Se sabe que la señalización de PXR activado inhibe la activación de las HSC tanto in vivo como in vitro [52, 53]. Nuestros resultados indican que el IPA puede participar en la eliminación de las HSC activadas al promover la apoptosis, reducir la fibrosis y el metabolismo mitocondrial, y mejorar las señales de supervivencia, que son procesos típicos que convierten el fenotipo de las HSC activadas en uno inactivado. Otra posible explicación para el mecanismo y la función potenciales del IPA en la apoptosis es que elimina las mitocondrias disfuncionales principalmente a través de la mitofagia (vía intrínseca) y la vía de señalización extrínseca del TNF (Tabla 1), que está directamente relacionada con la vía de señalización de supervivencia de NF-κB (Figura Suplementaria 7). Curiosamente, los genes enriquecidos relacionados con el IPA son capaces de inducir señales proapoptóticas y prosupervivencia en la vía apoptótica [54], lo que sugiere que el IPA podría inducir la vía apoptótica o la supervivencia al interactuar con estos genes. Sin embargo, aún no se conoce con certeza cómo el IPA induce la apoptosis o la supervivencia durante la activación de las CMH ni sus mecanismos.
El IPA es un metabolito microbiano que se forma a partir del triptófano de la dieta a través de la microbiota intestinal. Diversos estudios han demostrado que posee propiedades antiinflamatorias, antioxidantes y reguladoras epigenéticas en el entorno intestinal.[55] Diversos estudios han demostrado que el IPA puede modular la función de la barrera intestinal y reducir el estrés oxidativo, lo que podría contribuir a sus efectos fisiológicos locales.[56] De hecho, el IPA se transporta a los órganos diana a través de la circulación, y dado que comparte una estructura de metabolito principal similar con los derivados del triptófano, la serotonina y el indol, el IPA ejerce acciones metabólicas que resultan en destinos metabólicos competitivos.[52] El IPA puede competir con los metabolitos derivados del triptófano por los sitios de unión en enzimas o receptores, lo que podría alterar las vías metabólicas normales. Esto resalta la necesidad de realizar más estudios sobre su farmacocinética y farmacodinamia para comprender mejor su ventana terapéutica.[57] Queda por ver si esto también puede ocurrir en las células madre hematopoyéticas (CMH).
Reconocemos que nuestro estudio tiene algunas limitaciones. Para examinar específicamente las asociaciones relacionadas con el IPA, excluimos a los pacientes con diabetes mellitus tipo 2 (DM2). Reconocemos que esto limita la amplia aplicabilidad de nuestros hallazgos a pacientes con diabetes mellitus tipo 2 y enfermedad hepática avanzada. Aunque la concentración fisiológica de IPA en suero humano es de 1–10 μM [11, 20], se eligió una concentración de IPA de 1 mM con base en la concentración no tóxica más alta [15] y la tasa más alta de apoptosis, sin diferencia en el porcentaje de la población de células necróticas. Aunque se utilizaron niveles suprafisiológicos de IPA en este estudio, actualmente no hay consenso con respecto a la dosis efectiva de IPA [52]. Aunque nuestros resultados son significativos, el destino metabólico más amplio del IPA sigue siendo un área activa de investigación. Además, nuestros hallazgos sobre la asociación entre los niveles séricos de IPA y la metilación del ADN de las transcripciones hepáticas se obtuvieron no solo de células madre hematopoyéticas (HSC), sino también de tejidos hepáticos. Optamos por utilizar células LX-2 humanas basándonos en nuestros hallazgos previos del análisis del transcriptoma, que indican que el IPA está asociado con la activación de las células madre hematopoyéticas (CMH) [15], y que las CMH son las principales células implicadas en la progresión de la fibrosis hepática. El hígado está compuesto por múltiples tipos de células, por lo que se deben considerar otros modelos celulares, como el sistema de cocultivo de hepatocitos, CMH y células inmunitarias, combinado con la activación de caspasas y la fragmentación del ADN, así como el mecanismo de acción, incluyendo el nivel de proteínas, para estudiar el papel del IPA y su interacción con otros tipos de células hepáticas.