An integrative single-cell and spatial transcriptomics atlas highlights candidate regulatory factors in the development of gerbera capitulum

Resumen Técnico: Atlas Integrativo de Transcriptómica Espacial y de Célula Única del Desarrollo del Capítulo de Gerbera hybrida

Planteamiento del Problema
La familia Asteraceae, caracterizada por el capítulo (cabezuela floral), representa uno de los linajes de plantas con flores más exitosos evolutivamente. Aunque la complejidad morfológica del capítulo de Gerbera hybrida está bien documentada, los mecanismos reguladores a nivel de transcriptoma que gobiernan su desarrollo, particularmente durante la diferenciación de las floretes, siguen siendo poco comprendidos. Aunque se han identificado varios factores de transcripción clave de la familia MADS-box, el panorama regulador es incompleto y la dinámica espacial de la expresión génica dentro de la compacta estructura de la inflorescencia no se ha resuelto. Los estudios existentes carecen de una visión integrativa que combine la especificidad del tipo celular con el contexto espacial, lo que limita la capacidad de diseccionar las redes de regulación génica (GRN) que controlan el patrón del capítulo y el desarrollo vascular.

Metodología
Los autores emplearon un enfoque integrativo combinando la secuenciación de ARN de célula única (scRNA-seq) y la transcriptómica espacial de alta resolución (Stereo-seq) para construir un atlas dinámico del capítulo de G. hybrida.

• Recolección de Muestras: Se seleccionaron tres etapas sucesivas de desarrollo (S1–S3) de capítulos del cultivar de diploide doble G. hybrida SH6 para capturar la fase de diferenciación de las floretes. Se eliminaron las brácteas involucrales para facilitar el aislamiento de protoplastos.
• Secuenciación de Célula Única: Utilizando la plataforma 10x Genomics 3' v3, los autores generaron datos de scRNA-seq para las tres etapas. Tras el control de calidad, se retuvieron aproximadamente 25,790 células de alta calidad. Los datos fueron procesados mediante Scanpy, incluyendo normalización, detección de genes altamente variables, PCA, mapeo de difusión, y agrupamiento de Leiden.
• Transcriptómica Espacial: Se aplicó Stereo-seq a 19 muestras de capítulos que cubrían el rango S1–S3. Tras descartar los chips de baja calidad, se analizaron los datos de 14 muestras. Para abordar la dispersión de los datos y los límites de resolución, se agruparon dominios de 50×50 píxeles para crear "pseudo-células".
• Integración de Datos: Se utilizó el algoritmo cell2location para mapear el atlas de referencia de célula única sobre los datos de la transcriptómica espacial, infiriendo la distribución espacial de los tipos celulares.
• Análisis Secundario:
  • Anotación de Tipos Celulares: Los grupos (clusters) se anotaron utilizando la expresión de genes marcadores y el enriquecimiento de la Ontología Génica (GO) de los genes diferencialmente expresados (DEG).
  • Inferencia de Trayectorias: Se utilizó el paquete moscot para modelar las transiciones de tipos celulares y las trayectorias de desarrollo entre etapas.
  • Análisis de MADS-box: Un análisis filogenético identificó 90 genes MADS-box en el genoma de G. hybrida. El análisis de correlación espacial (basado en convolución) y el agrupamiento jerárquico se usaron para clasificar sus patrones de expresión.
  • Inferencia de Redes Reguladoras: El pipeline pySCENIC infirió las GRN a partir de los datos de scRNA-seq. Para el gen específico GAGL12, se identificaron posibles genes diana mediante la intersección de las predicciones de la GRN, el escaneo de motivos CArG-box (en regiones 3kb aguas arriba) y los DEG específicos de floema.
  • Validación Experimental: Se realizaron ensayos de levadura de dos híbridos (Y2H) para validar las interacciones proteína-proteína entre GAGL12 y otras proteínas de dominio MADS.

Resultados Clave

1. Atlas Celular y Dinámica: El estudio identificó 33 grupos agregados en 10 tipos celulares principales, incluyendo células en fase G2/M y S, epidermis, meristemos de inflorescencia/meristemos florales tempranos (IM/early-FM), meristemos florales (FM), endodermis, floema, procámbium/xilema, mesofilo fotosintético y primordios de brácteas.

   • Transiciones de Desarrollo: El análisis reveló un cambio significativo en la composición celular: las células IM/early-FM (abundantes en S1) disminuyeron a medida que el capítulo maduraba, mientras que las células FM aumentaron, lo cual es consistente con la conversión de meristemos en órganos florales.
   • Subtipos de Floema: Las células del floema se resolvieron en tres subpoblaciones (floema_a, b, c). El análisis de velocidad de ARN y de transición celular sugirió trayectorias de desarrollo distintas, donde el floema_b y el floema_c muestran progresión hacia estados terminales, mientras que el floema_a permanece estable.

2. Patrones de Expresión Espacial de Genes MADS-box: Los autores clasificaron 90 genes MADS-box en cuatro patrones distintos de expresión espacial:

   • Patrón 1 (Capítulo Completo): Incluye genes GSQUA (1–9) ampliamente expresados y GRCD2/7, sugiriendo roles en el desarrollo general del capítulo y de las floretes.
   • Patrón 2 (Superior): Enriquecido en meristemos y órganos florales, contiene genes involucrados en la identidad de órganos (ej., GAGA1/2, GGLO1, GDEF2, GRCD1/4/5/8).
   • Patrón 3 (Medio): Asociado con la expansión de tejidos y el soporte vascular. Este patrón destacó a GAGL12 y GANR1L10 como unicamente enriquecidos en las células del floema.
   • Patrón 4 (Base): Contiene únicamente GSOC2 y GSOC4 (clado TM3), potencialmente involucrados en la fotosíntesis y la tolerancia al calor.

3. Descubrimiento de GAGL12 como Regulador Vascular:

   • Especificidad Espacial: Se encontró que GAGL12 está específicamente enriquecido en las células putativas del floema, con una distribución espacial que coincide con las ubicaciones de floema inferidas en el capítulo.
   • Interacciones Proteicas: Los ensayos Y2H confirmaron que la proteína GAGL12 interactúa con otras proteínas de dominio MADS, incluyendo GRCD2/6/7, GSQUA1, GFLC6 y una proteína de tipo MIKC* (Ger_023266).
   • Predicción de Dianas: El análisis integrativo identificó genes diana putativos de GAGL12, incluyendo PDCB3 (proteína de unión a calosa de plasmodesmos 3) y PDLP6, ambos involucrados en la función de los plasmodesmos y el tráfico célula a célula.

Significación y Reivindicaciones
El artículo afirma proporcionar el primer atlas transcriptómico integrativo y dinámico del capítulo de G. hybrida durante la diferenciación de las floretes. Al combinar scRNA-seq y Stereo-seq, el estudio resuelve la heterogeneidad celular y el contexto espacial de la expresión génica que la secuenciación de volumen (bulk) no puede lograr.

Los autores enfatizan que este recurso sienta las bases para futuros estudios funcionales del control regulatorio de esta compleja estructura reproductiva. Específicamente, la identificación de GAGL12 como un candidato regulador de la vasculatura del capítulo, respaldada por su expresión específica en el floema, su interacción con reguladores vasculares conocidos y sus dianas predichas relacionadas con los plasmodesmos, ofrece una nueva hipótesis para el control molecular del desarrollo vascular en Asteraceae. El estudio también destaca numerosos genes MADS-box previamente no caracterizados con patrones espaciales específicos, proporcionando una hoja de ruta para la futura caracterización funcional. Los autores mantienen un tono modesto, señalando que, aunque las anotaciones de los tipos celulares están respaldadas por marcadores y enriquecimiento de GO, estas siguen siendo "putativas" a la espera de validación experimental en gerbera, y que las funciones biológicas exactas de los genes recién identificados aún están por revelarse.