Haplotype-resolved ancestry profiling reveals mosaic architecture of the lemon genome

Resumen Técnico: El perfilado de la ascendencia con resolución de haplotipos revela la arquitectura de mosaico del genoma del limón

Planteamiento del Problema
El limón (Citrus × limon) es un aloploide altamente heterocigoto derivado de tres especies ancestrales: mandarina (C. reticulata), pomelo (C. maxima) y cidra (C. medica). Aunque el mosaicismo estructural de los genomas de los cítricos resultante de la evolución reticulada está bien documentado, las consecuencias funcionales de este origen multiparental a nivel de haplotipo permanecen inexploradas. Específicamente, no está claro cómo la división de los segmentos genómicos ancestrales se traduce en una divergencia funcional, como la expresión génica diferencial o las presiones de selección, dentro de los dos distintos haplotipos del genoma del limón. Los estudios existentes suelen depender de ensamblajes genómicos globales (bulk), lo que oscurece las contribuciones específicas de los linajes parentales individuales al fenotipo del híbrido.

Metodología
Los autores realizaron un reanálisis sistemático de ensamblajes de genomas de telómero a telómero (T2T) y transcriptomas públicamente disponibles sin generar nuevos datos de secuenciación. El estudio utilizó ensamblajes de alta calidad, sin brechas o a escala cromosómica, para el limón (dos haplotipos, Cli_hapA y Cli_hapB), la mandarina, el pomelo y la cidra.

El marco analítico consistió en cuatro pasos integrados:

1. Anotación Estandarizada: Se aplicó un flujo de trabajo unificado a los cinco haplotipos (los dos haplotipos del limón más los tres progenitores) para asegurar la comparabilidad cruzada. Esto implicó el enmascaramiento de repeticiones (EDTA) y la predicción de genes mediante una combinación de métodos basados en homología (GeMoMa), basados en transcriptomas (Hisat2/Iso-Seq) y de novo (Braker2), fusionados mediante EVidenceModeler.
2. Asignación de la Ascendencia Filogenómica: Se empleó un enfoque de mejor golpe recíproco (best-reciprocal-hit) utilizando DIAMOND (identidad ≥95%, E-value ≤1e-5) para asignar cada gen de los haplotipos del limón a uno de los tres progenitores (mandarina, pomelo o cidra) o clasificarlos como "Incierto" (UN). Se utilizaron análisis de ventana deslizante (ventanas físicas de 1 Mb y ventanas de 500 genes) para visualizar los bloques de ascendencia y sus transiciones.
3. Análisis de Divergencia Evolutiva: Las tasas de sustitución sinónima (Ks) y no sinónima (Ka) se calcularon para pares ortólogos entre los haplotipos del limón y los progenitores utilizando el modelo YN00. Los genes con Ka/Ks > 1 fueron identificados como candidatos para la selección positiva.
4. Perfilado de la Expresión de Alelo Específico (ASE): Los datos de RNA-seq de tejidos de hojas y brotes de limón se alinearon al genoma resuelto por haplotipos concatenado. Se extrajeron los conteos de lecturas que se mapean de forma única a cada alelo, y el sesgo de expresión se determinó mediante una prueba exacta binomial de dos colas (FDR < 0.05). Los niveles de expresión (TPM) se compararon entre pares de genes colineales de diferentes orígenes ancestrales.

Resultados Clave

• Redefinición de Haplotipos y Arquitectura de Mosaico: El estudio confirmó que el genoma del limón es un mosaico estructural. El Cli_hapA (renombrado como el haplotipo "tipo naranja amarga") deriva predominantemente de la mandarina (38.3%) y el pomelo (23.5%), mientras que el Cli_hapB (el haplotipo "tipo cidra") deriva principalmente de la cidra (47.4%). Las visualizaciones de ventana deslizante revelaron grandes bloques de ascendencia pura con transiciones abruptas, mostrando particularmente que el Cli_hapA está dominado por señales de mandarina/pomelo y el Cli_hapB por señales de cidra. Aproximadamente el 37–39% de los genes en ambos haplotipos no pudieron ser asignados con confianza a ningún progenitor ("Incierto").
• Divergencia de Secuencia y Selección: Las distribuciones de Ks corroboraron las asignaciones de ascendencia, mostrando una divergencia significativamente menor entre el Cli_hapA y la mandarina/pomelo en comparación con la cidra, y viceversa para el Cli_hapB. Se identificaron aproximadamente entre 120 y 130 genes bajo selección positiva (Ka/Ks > 1) por comparación, lo que sugiere una evolución adaptativa en linajes ancestrales específicos.
• Mosaicismo Funcional mediante Sesgo de Expresión: El estudio identificó 7,515 pares de genes significativamente sesgados (de un total de 17,799 pares expresados). Crucialmente, se observó una asimetría funcional: los genes derivados de la mandarina y el pomelo en el Cli_hapA se expresan en niveles significativamente más altos que sus ortólogos derivados de la cidra en el Cli_hapB. Esto indica que los componentes ancestrales difieren no solo en su presencia, sino también en su actividad transcripcional.
• Genes Candidatos para el Fitomejoramiento: Al intersectar los genes con sesgo de ascendencia con los genes bajo selección positiva, los autores identificaron una lista corta de 55 genes candidatos. Estos incluyen transportadores (MATE, ABC), chaperonas (BiP) y proteínas reguladoras (PPR, F-box) involucradas en el desarrollo del fruto, la tolerancia al estrés y la fertilidad.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera demostración con resolución de haplotipos de que la evolución reticulada de los cítricos produce no solo un mosaicismo estructural, sino también un mosaicismo funcional a nivel de subgenoma. El estudio establece que el haplotipo "tipo naranja amarga" porta alelos ancestrales (mandarina/pomelo) transcripcionalmente activos que son sistemáticamente más expresados que los alelos derivados de la cidra en el otro haplotipo.

Los autores posicionan su marco integrado —que combina genomas T2T, asignación de ascendencia filogenómica, análisis de selección y perfilado de ASE— como un paradigma transferible para estudiar híbridos multiparentales. Afirman que este enfoque ofrece un "inventario listo para marcadores" de alelos con sesgo de ascendencia, permitiendo a los fitomejoradores priorizar haplotipos específicos para rasgos como la resistencia a enfermedades, el sabor y la restauración de la fertilidad sin requerir nuevos datos de secuenciación. El trabajo cierra la breancia entre la genómica descriptiva y las estrategias de mejoramiento conscientes de los haplotipos.