Single-cell and pseudobulk transcriptomic profiling reveals immune-enriched disease- associated signatures in human pancreatic islets from type 2 diabetes

Resumen Técnico: Perfilado Transcriptómico de Célula Única y Pseudobulk en Islotes Pancreáticos Humanos de Diabetes Tipo 2

Planteamiento del Problema
La diabetes tipo 2 (DT2) se caracteriza por una disfunción progresiva de las células de los islotes pancreáticos, sin embargo, los mecanismos moleculares que subyacen a esta disfunción siguen siendo complejos debido a la heterogeneidad celular del microambiente de los islotes. Los enfoques tradicionales de transcriptómica bulk suelen ocultar las señales de la enfermedad específicas de cada tipo celular, particularmente cuando las alteraciones son sutiles o se limitan a poblaciones menos abundantes, como las células inmunes o subtipos endocrinos específicos. Si bien la secuenciación de ARN de célula única (scRNA-seq) ofrece una mayor resolución, estudios previos han destacado los desafíos para interpretar los cambios asociados a la DT2 debido a la variabilidad en la composición de los donantes, la profundidad de la secuenciación y el diseño analítico. Este estudio aborda la pregunta de cómo la DT2 altera el paisaje celular y los programas de transcripción específicos de cada tipo celular en los islotes pancreáticos humanos a resolución de célula única, buscando determinar si las señales de la enfermedad son globales o específicas de un compartimento.

Metodología
El estudio realizó un análisis secundario de datos de scRNA-seq disponibles públicamente (acceso GEO GSE221156) que comprenden 37 muestras de islotes pancreáticos humanos (20 donantes con DT2 y 17 donantes no diabéticos). El flujo de trabajo analítico incluyó:

1. Preprocesamiento y Control de Calidad (QC): Las matrices de recuentos brutos se procesaron utilizando el paquete Seurat en R. Un filtrado estricto a nivel de célula retuvo células con 500–10,000 genes detectados, 1,000–50,000 UMIs y <10% de contenido mitocondrial. Este proceso redujo el conjunto de datos de ~260,000 a 115,104 células de alta calidad.
2. Integración y Agrupamiento (Clustering): Para armonizar la variación técnica entre muestras preservando la heterogeneidad biológica, los autores utilizaron la función IntegrateLayers de Seurat v5 con Análisis de Correlación Canónica (CCA). El agrupamiento basado en grafos identificó 35 clústeres globales y 21 clústeres refinados dentro de un subconjunto enriquecido en islotes.
3. Anotación: Los tipos celulares se anotaron utilizando genes marcadores canónicos (p. ej., INS para células beta, PRSS1 para células acinares, PTPRC para células inmunes) y se validaron mediante feature plots, violin plots y heatmaps.
4. Análisis de Expresión Diferencial (DE):
   • A nivel de célula: Se realizaron pruebas de Wilcoxon rank-sum dentro de los tipos celulares anotados para identificar DEGs entre células con DT2 y no diabéticas.
   • Pseudobulk: Para tener en cuenta la replicación biológica y evitar tratar las células del mismo donante como observaciones independientes, los recuentos brutos se agregaron a nivel de muestra. Se utilizó DESeq2 para modelar estos recuentos de pseudobulk con un modelo lineal generalizado binomial negativo, probando los efectos del estado de la enfermedad dentro de compartimentos específicos.
5. Enriquecimiento Funcional: Los genes marcadores significativos y los DEGs exploratorios se sometieron a recursos de bioinformática DAVID para el análisis de sobrerepresentación funcional.

Resultados Clave

• Paisaje Global: El análisis integrado reveló una sustancial heterogeneidad celular, pero no hubo una separación global clara entre las muestras con DT2 y las no diabéticas. El estado de la enfermedad no fue el conductor dominante de la varianza global; en su lugar, la identidad celular y la variabilidad inter-donante estructuraron los datos.
• Análisis de DE a Nivel de Célula: No se identificaron genes robustos asociados a la enfermedad en las principales poblaciones celulares (tipo beta, tipo acinar, tipo ductal) tras la corrección de pruebas múltiples.
• Análisis de Pseudobulk:
  • En los compartimentos tipo beta, tipo acinar y tipo ductal, ningún gen cumplió con el umbral de significancia estricto (FDR < 0.01 y |log2FC| > 2).
  • En el compartimento enriquecido en células inmunes, se detectó una señal exploratoria. El gen de quimiocina CXCL3 estuvo sobreexpresado en las muestras con DT2 con un log2 fold change de 2.521 y un p-valor ajustado de 0.03854 (FDR < 0.05), aunque no alcanzó el umbral más estricto de FDR < 0.01.
  • Otros genes (p. ej., CHI3L1, IL33, STK38L) mostraron significancia nominal en células inmunes pero no superaron las correcciones de FDR de alta confianza.
• Enriquecimiento Funcional: El análisis de DAVID de los genes marcadores significativos destacó firmas relacionadas con la matriz extracelular, proteínas secretadas, digestión/absorción de proteínas y actividad proteolítica. La lista exploratoria de genes del compartimento enriquecido en inmunidad también se enriqueció para los términos "Secreted" (Secretado) y "Extracellular region" (Región extracelular).

Significancia y Reivindicaciones
El artículo sostiene que las alteraciones transcripcionales en los islotes humanos con DT2 son sutiles, específicas de cada compartimento y se detectan mejor mediante un análisis de célula tipo y consciente del donante en lugar de comparaciones globales. Los autores argumentan que la ausencia de genes expresados diferencialmente amplios y con significancia de FDR en los principales tipos celulares sugiere que la DT2 no induce un cambio transcriptómico global y uniforme. En cambio, los cambios relacionados con la enfermedad parecen estar distribuidos en nichos celulares específicos.

El estudio respalda un modelo donde la señal de la enfermedad es enmascarada por la variabilidad inter-donante y la heterogeneidad celular, lo que requiere enfoques de pseudobulk para detectar efectos sutiles. La identificación de CXCL3 en el compartimento enriquecido en células inmunes se presenta como una señal candidata que refleja la comunicación inflamatoria localizada dentro del microambiente de los islotes, consistente con el papel conocido de la inflamación en la patofisiología de la DT2. Sin embargo, los autores mantienen un tono modesto, caracterizando estos hallazgos como "exploratorios" y señalando que el conjunto de datos carece del poder estadístico para calificar definitivamente a los genes individuales como significativos tras la corrección estricta de FDR. El estudio concluye que se requieren cohortes más grandes, conjuntos de datos independientes y validación ortogonal para confirmar si estas señales representan mecanismos reproducibles de la patología de los islotes en la DT2 humana.