Persistence of the coral Porites harrisoni in the Persian/Arabian Gulf is underpinned by strong selection on few genomic loci

Resumen Técnico: Persistencia de Porites harrisoni en el Golfo Pérsico/Arábico

Planteamiento del Problema
Los corales constructores de arrecifes se están acercando a sus límites térmicos superiores, con proyecciones que sugieren pérdidas catastróficas bajo un calentamiento oceánico acelerado. El Golfo Pérsico/Arábico (GPA) representa un entorno extremo donde las temperaturas del agua de mar en verano superan los 34 °C, sobrepasando los límites térmicos de la mayoría de los arrecifes tropicales por más de 2 °C. A pesar de estas duras condiciones, el coral pétreo Porites harrisoni persiste en el GPA. Sin embargo, la base genómica de esta tolerancia térmica, la estructura genética de la población que conecta a las poblaciones del GPA con entornos vecinos más benignos (como el Golfo de Omán, GO) y la especificidad de las simbiosis huésped-alga en este contexto, siguen siendo poco comprendidas. Comprender estos mecanismos es crítico para evaluar la viabilidad de estrategias de intervención activa, como el flujo génico asistido, para mejorar la resiliencia de los corales en otros lugares.

Metodología
El estudio empleó un enfoque de genómica de poblaciones para investigar P. harrisoni en tres sitios: dos en el extremo térmicamente extremo del sur del GPA (Arrecife Al Saada y Arrecife Saadiyat) y uno en el ambientalmente benigno Golfo de Omán (Al Aqah).

• Muestreo y Secuenciación: Se muestrearon un total de 119 colonias de coral en 2022. Se realizó la secuenciación de genoma completo (WGS) mediante shotgun en una plataforma Illumina NovaSeq6000, obteniendo una cobertura media de 45X.
• Identificación de Variantes y Filtrado: Las lecturas se procesaron siguiendo las mejores prácticas de GATK. Se aplicaron dos canales de filtrado de SNPs distintos: uno para el análisis de la estructura poblacional (eliminando variantes de baja frecuencia y clones) y uno más permisivo para la detección de barridos selectivos (reteniendo variantes raras para capturar señales de fijación).
• Estructura Poblacional: Los análisis incluyeron Análisis de Componentes Principales (PCA), ADMIXTURE para la estimación de la ascendencia, cálculos de $F_{ST}$ por pares y Stairway Plot v2 para la reconstrucción de la historia demográfica.
• Asociación Huésped-Simbionte: Los genotipos del huésped se correlacionaron con los ensamblajes de algas simbiontes (perfiles de tipo ITS2) utilizando pruebas de Mantel para evaluar las asociaciones genotipo-genotipo.
• Firmas de Selección: Las regiones candidatas bajo selección positiva se identificaron integrando las estadísticas de $F_{ST}$ de Weir & Cockerham (diferenciación genética) y la homozigosidad extendida de haplotipos entre poblaciones (XP-EHH). Los SNPs atípicos fueron anotados mediante SnpEff, y se evaluaron la diversidad de nucleótidos ($\pi$) y el decaimiento del desequilibrio de ligamiento (LD) para caracterizar la naturaleza de los barridos selectivos.

Resultos Clave

• Estructura Poblacional: La diferenciación genética general entre los tres sitios fue baja ($F_{ST}$ en un rango de 0.021 a 0.041), lo que indica flujo génico a través del Estrecho de Ormuz. Sin embargo, el análisis de ADMIXTURE ($K=2$) reveló que las poblaciones del GPA comparten un linaje ancestral distinto al del GO, siendo que el sitio más al sur del GPA (SA) mostró una mezcla mínima. La reconstrucción demográfica sugiere una reducción significativa en el tamaño efectivo de la población ($N_e$) que coincide con la formación del GPA durante el Holoceno (~6,000 años atrás), indicativo de un efecto fundador o un cuello de botella.
• Especificidad Huésped-Simbionte: Aunque la similitud genética del huésped no dictó estrictamente la identidad del simbionte, se encontró una asociación moderada pero significativa entre las distancias genéticas del huésped y de Cladocopium. Las colonias del GPA estaban predominantemente asociadas con el termotolerante Cladocopium thermophilum, mientras que las colonias del GO albergaron Durusdinium trenchii y Symbiodinium microadriaticum.
• Firmas Genómicas de Selección: El estudio identificó dos regiones genómicas específicas bajo fuerte selección positiva en los corales del GPA que no se observaron en los del GO:
  1. Contig 29: Exhibió los valores de $F_{ST}$ más altos (hasta 0.93) y haplotipos extendidos. Esta región contiene genes relacionados con la remodelación de la matriz extracelular (ECM) (p. ej., COL4A1, COL2A1), señalización de TGF-$\beta$ (ACVR2B, BMPR2) y regulación del ciclo celular (CDK2, ORC4).
  2. Contig 62: Mostró haplotipos excepcionalmente largos (XP-EHH normalizado hasta 10.47) y una alta diferenciación ($F_{ST}$ hasta 0.69). Esta región está enriquecida con genes involucrados en la respuesta inmune, respuesta al estrés (p. ej., KDELC1), apoptosis (PDCD6) y procesamiento de ARN.
• Naturaleza de la Selección: Los gráficos de frecuencia alélica mostraron una casi fijación de alelos específicos en los corales del GPA, pero no en los del GO. La falta de una reducción pronunciada en la diversidad de nucleótidos ($\pi$) junto con haplotipos extendidos sugiere un barrido selectivo "suave", donde la selección actuó sobre la variación genética existente en lugar de sobre una única mutación novedosa.

Significancia y Reivindicaciones
El artículo sostiene que la persistencia de P. harrisoni en el extremo del GPA está sustentada por una fuerte selección direccional que actúa sobre un pequeño conjunto de loci genómicos. La identificación de "SNPs no negociables" —alelos llevados a una casi fijación en las poblaciones del GPA— destaca el papel esencial de estas variantes genéticas específicas para sobrevivir al estrés térmico y de salinidad extremos.

El estudio afirma proporcionar el primer paso para identificar objetivos genómicos para futuras estrategias de intervención activa, como el flujo génico asistido mediante la cría selectiva, para mejorar la resiliencia de los corales. Sin embargo, los autores enfatizan una restricción crítica: la adaptación de los corales del GPA está estrechamente vinculada a las condiciones ambientales específicas del Golfo (incluyendo la alta salinidad). Por consiguiente, la transferencia de estos genotipos a entornos más fríos y de menor salinidad (como el GO) puede no ser viable a largo plazo. Los hallazgos subrayan que, si bien la diversidad genética regional es un recurso valioso, el estrecho régimen selectivo de los entornos extremos puede limitar el potencial adaptativo futuro, haciendo que las estrategias de mitigación climática global sean más urgentes que las intervenciones locales por sí solas.