Single-cell RNA sequencing of organ-specific endothelial landscapes and intercellular networks in endocrine organs

Resumen Técnico: Secuenciación de ARN de célula única del paisaje endotelial específico de órganos y redes intercelulares en órganos endocrinos

Planteamiento del Problema
Las glándulas endocrinas dependen de una vasculatura especializada para orquestar el tráfico hormonal y la homeostasis tisular. Si bien las células endoteliales (CE) en estos órganos exhiben características estructurales distintivas —como fenestraciones con diafragmas que facilitan la permeabilidad selectiva—, el panorama transcriptómico subyacente a la heterogeneidad de las CE en los órganos endocrinos sigue siendo insuficientemente comprendido. Los atlas previos de secuenciación de ARN de célula única (scRNA-seq) a menudo han excluido a los órganos endocrinos debido a desafíos técnicos, dejando ambigua la heterogeneidad molecular y funcional de las CE dentro de estos tejidos. Además, la relación entre las firmas moleculares y la heterogeneidad estructural (fenestrada continua, no fenestrada continua y sinusoidal) no ha sido mapeada de manera exhaustiva en un contexto endocrino.

Metodología
Los autores generaron un atlas exhaustivo de scRNA-seq de células endoteliales de seis órganos: tres órganos endocrinos (tiroides, glándula adrenal, glándula pituitaria), una referencia exocrina/estructural (riñón) y dos referencias estructurales (hígado y corteza cerebral).

• Generación y Procesamiento de Datos: Tras el filtrado de calidad y la corrección de lotes (batch correction), se analizaron 52,970 células. Se realizó un agrupamiento (clustering) para identificar las principales líneas celulares, seguido del aislamiento y re-agrupamiento de las poblaciones endoteliales.
• Clasificación Estructural: Las células endoteliales se categorizaron en subtipos arteriales (aCEs), venosos (vCEs) y capilares (cCEs). Se excluyeron las no-cCEs (aCEs, vCEs y cCEs enriquecidas en Hba) para centrarse en la heterogeneidad de las cCEs organotípicas.
• Análisis Computacional:
  • Análisis de Redes Funcionales: Se utilizó el Análisis de Sobre-Representación (ORA) para construir redes de términos funcionales, visualizando las relaciones entre vías biológicas a través de los subtipos estructurales.
  • NMF de Consenso (cNMF): Se utilizó para identificar programas transcripcionales y puntuaciones de uso de módulos dentro de las cCEs.
  • Inferencia de Regulones: Se empleó pySCENIC para reconstruir redes de factores de transcripción y sus dianas, y para cuantificar la actividad de los regulones.
  • Interacción Célula-Célula (CCI): Se utilizó el marco de trabajo LIANA (CellPhoneDB) para inferir interacciones ligando-receptor entre las células endoteliales y otras líneas celulares principales (epiteliales, inmunes, fibroblastos, células murales).
  • Análisis Comparativo: Se calcularon los índices de similitud de Jaccard para evaluar los solapamientos de conjuntos de genes específicos de cada órgano. El Análisis de Variación de Conjuntos de Genes (GSVA) cuantificó la actividad de las vías inflamatorias.
  • Validación: Los hallazgos se validaron contra conjuntos de datos de scRNA-seq humano y datos transcriptómicos de volumen (bulk) de The Cancer Genome Atlas (TCGA) para el análisis de supervivencia en el carcinoma papilar de tiroides (PTC).

Resultos Clave

1. Heterogeneidad Estructural y Molecular: El estudio confirmó que las cCEs pueden estratificarse en tres subtipos estructurales con firmas moleculares distintas:
   • CEs Fenestradas Continuas: Enriquecidas en órganos endocrinos y riñón; caracterizadas por la expresión de Emcn, Plvap y Ehd4, reguladas por Meox1, Ebf1, Sox4 e Irf7. Muestran un enriquecimiento en la señalización de VEGF/VEGFR2.
   • CEs Sinusoidales: Presentes en el hígado; caracterizadas por Mrc1 y Stab2, reguladas por Maf, Meis1 y Nr2f2.
   • CEs Continuas No Fenestradas: Presentes en el cerebro; caracterizadas por Cldn5, Slco1a4 y Spock2, reguladas por Foxf2 y Foxc1.
2. Poblaciones Pan-cEC: Se identificaron dos poblaciones de cCE universales presentes en todos los órganos:
   • Common EC1: Enriquecida en genes de proteínas ribosomales y metabolismo de aminoácidos.
   • Common EC2: Enriquecida en genes mitocondriales y fosforilación oxidativa.
     Ambas poblaciones convergen en el regulador transcripcional TEAD4 (eje Hippo–YAP/TAZ), lo que sugiere una arquitectura regulatoria metabólica y de citoesqueleto central.
3. Especialización Organotípica: El agrupamiento jerárquico reveló que los órganos endocrinos (tiroides, pituitaria, adrenal) y el riñón comparten firmas génicas específicas, mientras que el hígado y el cerebro forman grupos distintos.
   • Adrenal/Riñón: Asociados con la remodelación vascular (Npr3, Pbx1).
   • Pituitaria/Tiroides: Asociados con la respuesta al estrés (c-Fos, Igfbp7).
   • Tiroides: Únicamente enriquecida para genes de identidad linfática (Foxp2, Prdm1) y actividad antioxidante (Rbp7).
4. Interfaz Endotelio-Inmune: Un hallazgo principal es la regulación diferencial de la comunicación endotelio-inmune.
   • Órganos Inmunoactivos: La glándula adrenal, el riñón, el hígado y la pituitaria exhiben una fuerte interacción endotelio-inmune, caracterizada por regulones proinflamatorios (p. ej., Rarb) e interacciones que involucran vías de TNF/IFN y quimiocinas de tráfico de leucocitos.
   • Órganos Inmunoprivilegiados: El cerebro y la tiroides muestran una interfaz endotelio-inmune notablemente restringida. Las cCEs de la tiroides activan reguladores antiinflamatorios (Klf4, Prdm1) y muestran señales de interacción débiles con las células T, reflejando el fenotipo de inmunoprivilegio del cerebro.
5. Remodelación por Envejecimiento y Cáncer: El análisis de datos humanos de la tiroides reveló que el estado de inmunoprivilegio se ve comprometido durante el envejecimiento y el cáncer.
   • Envejecimiento: El tejido tiroideo normal de individuos de edad avanzada muestra una regulación al alza de las firmas de IFN-α/γ y un aumento en las interacciones endotelio-inmunes en comparación con el tejido joven.
   • Cáncer (PTC): El carcinoma papilar de tiroides presenta una pérdida de la interfaz vascular protectora, con firmas enriquecidas de inflamación activa (TNF-α/NF-κB, EMT, angiogénesis) e interacciones fuertes entre cCE y células T/mieloides. El par ligando-receptor MIF-ACKR3 fue identificado como específico de las interacciones en PTC y se asoció significativamente con la supervivencia de los pacientes.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona el primer atlas de célula única de las células endoteliales enfocado específicamente en órganos endocrinos. El trabajo define la base molecular de la especialización vascular endocrina, demostrando que la heterogeneidad endotelial no es meramente una adaptación estructural, sino una necesidad funcional moldeada por interacciones recíprocas con las células parenquimatosas.

Una contribución clave es la identificación de la tiroides como un órgano inmunoprivilegiado a nivel de célula única, un hallazgo que previamente solo se había hipotetizado en reportes históricos. El estudio postula que la ruptura de esta interfaz endotelio-inmune protectora durante el envejecimiento y la tumorigénesis contribuye a la patogénesis tiroidea. Al delinear los nichos vasculares organotípicos y las redes intercelulares, el artículo ofrece una nueva perspectiva sobre el papel de la plasticidad endotelial en el mantenimiento de la homeostasis orgánica y su desregulación en la enfermedad. Los hallazgos están respaldados por la validación transespecie en conjuntos de datos humanos, lo que sugiere relevancia traslacional para la comprensión de las contribuciones vasculares a los trastornos endocrinos y el cáncer.