Resumen Técnico: Entropía de la Herencia en Árboles de Linaje Celular

1. Planteamiento del Problema

Cuantificar el impacto de la transmisión hereditaria dentro de los árboles de linaje celular sigue siendo un desafío fundamental en biología. Si bien el perfilado de célula única permite la caracterización de los estados celulares, combinar esto con el rastreo de linaje clásico para seguir las ramificaciones hereditarias es difícil. Un desafío específico surge con la inactividad (células que dejan de dividirse, incluyendo la senescencia, diferenciación o quiescencia). A diferencia de la proliferación activa, donde las mutaciones hereditarias amplifican el fenotipo en la progenie, la inactividad es autolimitada: una vez que una célula se vuelve inactiva, su rama se termina, borrando la información descendente. Esto dificulta determinar si la emergencia de la inactividad es impulsada por eventos estocásticos aleatorios o por factores hereditarios (epigenéticos) que regulan la salida del ciclo celular. Además, distinguir entre mutaciones genéticas y la herencia epigenética es crucial, particularmente en colonias de células madre donde la heterogeneidad inter-colonia afecta el potencial terapéutico.

2. Metodología

Los autores introducen un método independiente del modelo basado en la entropía de Shannon aplicada a la topología de los árboles de linaje, centráéndose específicamente en la distribución de las células inactivas.

A. Adquisición de Datos

El estudio utiliza el rastreo de linaje de célula única de alta definición de 32 colonias clonales de Células Estromales de Médula Ósea (BMSC) humanas. Los experimentos se realizaron mediante microscopía de contraste de fase por lapso de tiempo hasta que todas las células replicantes alcanzaron la séptima generación ($k_{max}=7$). Las células se clasificaron como:

• Activas: Mitosis observada.
• Inactivas (G0/Muertas): Sin división durante 84 horas o muerte celular.
• Hojas: Células en $k=7$ rastreadas solo hasta su nacimiento.

B. Desequilibrio de Inactividad

Para cuantificar el impacto topológico de la inactividad, los autores definen el desequilibrio de inactividad ($I_m$) para cada mitosis $m$ en el árbol:
$$I_m = |M_{m}^{izquierda} - M_{m}^{derecha}|$$
donde $M$ representa el número de hojas faltantes en $k_{max}$ en las ramas izquierda y derecha. Un gran desequilibrio en un nodo específico sugiere una divergencia en la probabilidad de inactividad entre las dos ramas hijas, lo que podría indicar un evento de mutación aguas arriba.

C. Entropía de la Herencia

La métrica central es la entropía de la herencia ($S$), calculada a partir de los pesos normalizados ($w_m$) de los desequilibrios en todas las mitosis del linaje:
$$w_m = \frac{I_m}{\sum_n I_n}, \quad S = -\sum_m w_m \log w_m$$

• Entropía Alta: Los desequilibrios están dispersos de manera homogénea en el árbol, lo que sugiere fluctuaciones aleatorias y no hereditarias.
• Entropía Baja: Los desequilibrios se concentran en unos pocos nodos, indicando un evento hereditario específico (mutación) que alteró la probabilidad de inactividad para una rama descendente específica.

D. Conjuntos Nulos y Validación

Para determinar la significancia estadística, la entropía del árbol biológico se compara contra tres conjuntos nulos no hereditarios:

1. Galton-Watson (GW): Probabilidad constante de inactividad ($q$) para todas las células.
2. GW No Homogéneo (NHGW): Probabilidad de inactividad dependiente de la generación ($q_k$).
3. Remezclado Aleatorio (RR): El método preferido. Esto implica desordenar las posiciones de las células dentro de cada generación preservando el número total de células inactivas por generación. Esto corta todos los vínculos hereditarios madre-hija manteniendo las restricciones estructurales del árbol.

El Valor-P se define como la fracción de árboles del conjunto nulo con una entropía $S \leq S_{biológica}$. Un valor-P bajo indica una fuerte evidencia de herencia.

E. Ajuste de Modelo (SKH)

Los autores proponen un modelo Estocástico de Kirkwood-Holliday (SKH) para describir los datos. En este modelo:

• Las células silvestres mutan a un estado "mutado" con una probabilidad $p$.
• Las células mutadas (y sus descendientes) tienen una probabilidad $q$ de volverse inactivas.
• Crucialmente, el modelo incorpora un retraso de mutación ($l$), el número de generaciones entre el evento de mutación y la expresión fenotípica de la inactividad.

3. Contribuciones Clave

1. Entropía de la Herencia: Una métrica novedosa y agnóstica al modelo para detectar la estructura hereditaria en árboles de linaje basándose únicamente en la topología de las ramas truncadas (inactivas).
2. Validación de Modelos Nulos: Un estudio de benchmarking riguroso que demuestra que mientras los modelos GW y NHGW fallan en distinguir árboles hereditarios de árboles no hereditarios no homogéneos (produciendo falsos positivos), el método de Remezclado Aleatorio (RR) es robusto, sensible y agnóstico al modelo.
3. Cuantificación del Retraso de Mutación: La capacidad de localizar la mitosis específica donde cambió la probabilidad de inactividad y medir el retraso de mutación (el retraso entre el cambio epigenético y la expresión fenotípica).
4. Modelo Estocástico de Kirkwood-Holliday: Un modelo estocástico generalizado que captura con éxito las características hereditarias de los linajes de BMSC, incluyendo el retraso entre la mutación y la inactividad.

4. Resultados

• Estructura Hereditaria: En el conjunto de datos experimentales de 32 colonias de BMSC, 21 colonias (75%) mostraron un carácter hereditario significativo ($P < 0.05$) utilizando la prueba RR. Las colonias restantes tenían muy pocas células inactivas para producir un resultado estadísticamente significativo, lo que no significa necesariamente una falta de herencia.
• Análisis de Entropía: Los árboles biológicos exhibieron consistentemente una entropía significativamente menor que sus contrapartes remezcladas aleatoriamente, probando que la distribución de las células inactivas no es aleatoria sino que sigue un patrón hereditario.
• Retraso de Mutación: El retraso de mutación promedio fue de aproximadamente 2 generaciones. Esto confirma que el cambio epigenético que conduce a la inactividad ocurre aguas arriba de la expresión fenotípica real.
• Ajuste del Modelo: El modelo SKH, ajustado a los datos, arrojó parámetros $p \approx 0.1$ (probabilidad de mutación) y $q \approx 0.46$ (probabilidad de inactividad para células mutadas). Los árboles sintéticos generados por SKH con estos parámetros se asemejaron estrechamente a la estructura topológica de los linajes experimentales de BMSC.
• Espacial vs. Topológico: El análisis confirmó que el agrupamiento de células inactivas no es un artefacto de la proximidad física (correlación espacial) en la placa de cultivo, sino una característica topológica del linaje.

5. Significado y Reivindicaciones

El artículo afirma que la heterogeneidad observada en las colonias de BMSC humanas es impulsada por factores hereditarios epigenéticos que regulan la salida del ciclo celular, en lugar de mecanismos puramente estocásticos o genéticos.

• Robustez: La prueba basada en la entropía proporciona una forma simple, robusta y agnóstica al modelo para evaluar el origen hereditario del crecimiento diferencial en cualquier árbol de linaje.
• Relevancia Clínica: Los hallazgos sugieren que la heterogeneidad inter-colonia observada en los trasplantes de BMSC (un obstáculo importante para las terapias de regeneración ósea) tiene un origen hereditario epigenético. Comprender y controlar estos mecanismos hereditarios podría mejorar la consistencia de las terapias con células madre.
• Avance Metodológico: Al establecer un "retraso de mutación", el método permite a los investigadores rastrear hacia atrás el origen de los cambios fenotípicos incluso cuando el fenotipo (inactividad) es auto-borrable.
• Direcciones Futuras: Los autores señalan que, aunque el método actual se centra en la inactividad, la correlación entre la topología del linaje y los tiempos de división sugiere que la duración del ciclo celular también puede estar regulada por mecanismos hereditarios epigenéticos, aunque esto requiere más estudio.

Los autores mantienen la modestia, reconociendo que el modelo SKH es una simplificación y que modelos más complejos podrían capturar mejor la varianza en los retrasos de mutación. Enfatizan que el método detecta cambios en la probabilidad; si un linaje no tiene variación en la expresión de inactividad, la prueba no detectará la herencia, incluso si existe un mecanismo hereditario.