Radial Interaction Tomography: Recognizing Non-Transitive Evolutionary Games from One Range-Expansion Image

Resumen Técnico: Tomografía de Interacción Radial

Formulación del Problema
El artículo aborda un problema de visión por computadora inverso: la recuperación de la historia de las interacciones evolutivas de pares a partir de una única imagen de un punto final de una expansión microbiana. Mientras que trabajos previos utilizaron la propiedad de "frente congelado" de las colonias en expansión (donde el radio se correlaciona con el tiempo de deposición) para estimar la aptitud relativa constante o modelar el vagabundeo estocástico, estos enfoques asumen una jerarquía de aptitud escalar transitiva. Este artículo investiga si una única imagen de punto final puede revelar interacciones no transitivas (por ejemplo, dominancia cíclica como piedra-papel-tijeras o antagonismo mediado por bacteriocinas) que no pueden explicarse mediante un ordenamiento de aptitud escalar. El desafío central es distinguir las señales geométricas de la competencia cíclica de los artefactos del ruido de imagen, el vagabundeo de los límites o las variaciones de la aptitud escalar.

Metodología
Los autores proponen la Tomografía de Interacción Radial, un flujo de trabajo determinista que trata la imagen de la colonia como un objeto de reconocimiento de patrones en lugar de depender de incrustaciones (embeddings) aprendidas o redes neuronales.

1. Observable Geométrico: El método extrae un campo de flujo de frontera indexado por radio, $q_{ij}(\rho)$, definido como la derivada del ángulo desenrollado de un sector de la frontera con respecto a una coordenada de deposición $\rho$ (típicamente $\log(r/r_0)$).
2. Descomposición de Hodge: El vector de flujo de borde observado se descompone en un componente compatible con escalares (gradientes de un potencial de nodo, que representa la aptitud transitiva) y un componente residual (flujos no potenciales, que representan interacciones cíclicas). Esta descomposición está ponderada por una matriz de precisión para manejar el ruido.
3. Análisis de Identificabilidad: El artículo establece límites estrictos sobre lo que es recuperable:
   • No identificabilidad del Reloj Radial: El tiempo físico no puede recuperarse sin una calibración independiente; solo la coordenada de deposición es identificable.
   • No identificabilidad del Mecanismo: La imagen revela el flujo de frontera por pares, pero no velocidades de frente aisladas ni mecanismos microscópicos sin datos adicionales.
   • Diseño de Contacto Completo: Un teorema demuestra que para observar cada par en un juego antisimétrico de $n$ genotipos en una frontera no restringida, se requiere un mínimo de $\binom{n}{2}$ (para $n$ impar) o $\binom{n}{2} + n/2$ (para $n$ par) fronteras. Para $n=4$, se requieren al menos ocho fronteras, no seis.
4. Prueba Estadística: El flujo de trabajo emplea una prueba de razón de verosimilitud generalizada Gaussiana exacta (GLRT) para detectar residuos cíclicos, junto con correcciones de Bonferroni para el escaneo de intervalos a lo largo de la trayectoria radial.

Contribuciones Clave

• Formulación de la Tomografía de Interacción Radial: Definición del problema inverso de recuperar patrones transitivos y cíclicos indexados por radio a partir de una única imagen.
• Separación del Observable y la Interpretación: Desacoplamiento explícito del campo de flujo de borde geométrico de las interpretaciones biológicas (como la velocidad de frente aislada), evitando inferencias inválidas.
• Teorema de Diseño de Contacto: Establecimiento del número mínimo de fronteras de sectores requeridos para resolver completamente las interacciones por pares en un diseño circular.
• Benchmarking Determinista: Creación de una suite de pruebas que incluye imágenes de punto final analíticas, recorridos de ida y vuelta a nivel de píxel, pruebas de estrés de nulidad escalar y trazado de imágenes públicas, todo ello sin redes neuronales.
• Control de Diseño Activo: Extensión del resultado inverso a una capa de control donde los residuos cíclicos recuperados sirven como señales de forzamiento para el control de reacción-difusión adimensional y la optimización del diseño de cepas.

Resultados Experimentales
El artículo valida el método a través de varios benchmarks deterministas y estocásticos:

• Recuperación a Nivel de Píxel: En un "recorrido de ida y vuelta de píxeles" (renderizando una imagen borrosa/con ruido y recuperando flujos), el método logró un RMSE de flujo de borde de $\approx 2.86 \times 10^{-3}$ y un RMSE de norma cíclica de $\approx 2.95 \times 10^{-3}$.
• Pruebas de Estrés de Nulidad Escalar: Cuando la verdad de terreno es puramente transitiva (escalar), el método identifica correctamente falsos positivos solo cuando los factores de molestia (ruido, extinción, distorsión) son severos. Bajo alta molestia, el sistema está diseñado para abstenerse en lugar de inferir ecología no transitiva (la tasa de abstención aumenta apropiadamente, por ejemplo, una tasa de detección falsa de 0.564 bajo la peor molestia).
• Validación Mecanicista: Utilizando un simulador de frente congelado no basado en aprendizaje con cuatro genotipos, el método recuperó con éxito un anillo estable y detectó un residuo cíclico con una norma máxima de 0.358, coincidiendo con la interacción cíclica programada.
• Respuesta a la Dosis: Un barrido de las escalas de antagonismo cíclico programadas mostró una fuerte correlación ($r=0.776$) entre la escala programada y el residuo cíclico recuperado.
• Control Activo: Los residuos recuperados se utilizaron para impulsar un barrido de control de reacción-difusión. Un diseño refinado por estado adjunto redujo la densidad de la población no deseada en un 99.35% mientras retenía el 161.39% de la densidad beneficiosa.
• Robustez: Una prueba de robustez de Monte Carlo con 1,000,000 de trayectorias confirmó que los diseños seleccionados mantenían al menos un 40% de reducción de estado adverso y un 85% de retención de la cepa beneficiosa bajo perturbaciones de parámetros.
• Auditoría de Imagen Pública: El flujo de trabajo se aplicó a una figura pública de Weinstein et al., pasando como un "control de sanidad de visión por computadora", pero señalando explícitamente que las anotaciones de publicación impiden su uso como verdad de terreno cuantitativa bruta.

Significancia y Reivindicaciones
El artículo afirma que una única imagen de punto final contiene información geométrica suficiente para reconstruir historias de interacción por pares y detectar firmas no transitivas, siempre que el diseño experimental cumpla con criterios específicos de completitud de contacto. La contribución principal es la reivindicación de reconocimiento visual falsable: la capacidad de separar los observables geométricos de las interpretaciones biológicas posteriores.

Los autores enfatizan que, si bien el método recupera con éxito historias de interacción de "frente congelado" e impulsa simulaciones de diseño activo, no pretende reemplazar los experimentos físicos. La capa de "diseño-control activo" se presenta como un objetivo de simulación adimensional, y el artículo establece explícitamente que la traducción biológica requiere calibración independiente, tensores de fenotipo medidos y réplicas físicas. El trabajo sirve como un marco riguroso para interpretar imágenes de expansión de rango, definiendo los límites de lo que se puede inferir de la geometría por sí sola y proporcionando una herramienta determinista para la generación de hipótesis en la ecología evolutiva.