Ponente
Descripción
La fragmentación del bosque de niebla andino reduce el flujo génico e incrementa la deriva y la endogamia en poblaciones aisladas. La palma de cera del Quindío (Ceroxylon quindiuense), árbol nacional de Colombia categorizada Vulnerable por la UICN, ha perdido más del $50\%$ de su población en tres generaciones (Minambiente, Plan de conservación, manejo y uso sostenible de la palma de cera del Quindío, Universidad Nacional de Colombia, 2015). Su ciclo reproductivo de hasta 80 años implica que los efectos genéticos de la fragmentación se manifiestan con décadas de retraso. Estudios previos documentaron alta estructuración genética en la especie (González-Rivillas, N. et al., BioRxiv, 443960, 2018; Chacón-Vargas, K. et al., Biodivers. Conserv. 29, 283–302, 2019), pero no identificaron rutas críticas de flujo génico. La teoría de grafos permite complementar esta aproximación modelando poblaciones como nodos y su conectividad como aristas dirigidas y ponderadas (Dyer, R.J. & Nason, J.D., Mol. Ecol. 13, 1713–1727, 2004).
Se caracterizó la diversidad y estructura genética de $212$ individuos de ocho poblaciones distribuidas en la Cordillera Central (Quindío, Tolima, Caldas y Valle del Cauca), genotipadas con diez marcadores microsatélites. Se empleó análisis multivariado (DAPC) e inferencia de ancestría (snmf) para identificar grupos genéticos y asignar individuos a unidades poblacionales. A partir de los $F_{ST}$ pareados se construyó un grafo dirigido ponderado, con direccionalidad inferida desde $H_e$ y aristas filtradas bajo el umbral $1 − F_{ST} ≥ 0,85$ (equivalente a $F_{ST} ≤ 0,15$), validado por la distribución bimodal de los valores pareados y sustentado en los rangos de Wright. Sobre el grafo filtrado se calcularon PageRank, betweenness ponderado, fuerza de nodos y robustez por eliminación secuencial usando la eficiencia global:
\begin{equation}
E(G) = \frac{1}{n(n-1)} \sum_{i \neq j} \frac{1}{d_{ij}}
\end{equation}
donde $d_{ij}$ es la distancia del camino más corto entre los nodos $i$ y $j$. La relación distancia-diferenciación se evaluó con el test de Mantel.
Se identificaron $337$ alelos con $H_e=0,768, Ho=0,384$ y $F_{IS}=0,500$. La diferenciación fue significativa ($F_{ST}=0,171; p < 0,001$). El DAPC y snmf $(K=4)$ identificaron dos grupos: G1: Anaime(AN), Cocora1(CO) flanco occidental , Pijao(PI), San Félix(SF), Tenerife(TE) y G2: (Cocora2(CX) flanco oriental, La Línea(LI), Toche(TO)), CO-CX se encuentran separados por el río Quindío. El grafo filtrado reveló un núcleo de conectividad activa en G1 y tres nodos completamente aislados en G2. CO, AN y TE actuaron como emisores netos; PI y SF como receptores. Anaime (caída en eficiencia global: $28,52\%$) y Pijao ($24,39\%$) resultaron los nodos estructuralmente más críticos. El test de Mantel no detectó aislamiento por distancia ($r=−0,306; p=0,735$), siendo la paradoja CO-CX el caso más ilustrativo: poblaciones geográficamente contiguas, pero genéticamente desconectadas.
La integración de genética para la conservación y teoría de grafos reveló una arquitectura modular con conectividad restringida entre grupos. AN y PI son los nodos críticos del sistema, cuya perdida comprometería la cohesión genética del núcleo conectado. Se prioriza la restauración del corredor AN-PI-CO y la protección in situ de LI, TO y CX como unidades evolutivas de conservación.