El euroasiático basal como espejismo: un argumento steelman a favor de que la «ascendencia fantasma» es un artefacto de la selección

TL;DR

• «Basal Eurasian» (BE) es un componente del modelo introducido para satisfacer restricciones de correlación de frecuencias alélicas en grafos de ADN antiguo, estrechamente vinculado con una menor compartición de alelos neandertales en los primeros habitantes del Cercano Oriente (Lazaridis et al. 2016).
• Los grafos de mezcla y los estadísticos f suponen que la mayor parte de la covarianza genómica es generada por deriva neutral en un grafo; una selección intensa y geográficamente estructurada puede añadir covarianza sistemática que el modelo «explica» inventando ramas fantasma (Maier et al. 2023; Schrider et al. 2016).
• La selección purificadora contra la introgresión neandertal es real, en gran medida poligénica y mediada por selección ligada; la depuración diferencial puede modificar las estadísticas de compartición neandertal sin requerir un linaje profundamente divergente y carente de ascendencia neandertal (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).
• Hipótesis steelman: una selección masiva durante el Pleistoceno tardío y el Holoceno a favor de la cognición en linajes del Cercano Oriente y Eurasia occidental (poligénica, enriquecida en anotaciones y estructurada), junto con la selección ligada y la depuración de ascendencia arcaica, genera el patrón residual similar al de BE, y el ajuste del grafo convierte ese sesgo en un fantasma «Basal Eurasian».
• Esta hipótesis es comprobable: BE debería reducirse en particiones con alta recombinación e intergénicas; BE debería mostrar una intensidad variable a lo largo del tiempo en genomas del Cercano Oriente fechados; y la señal similar a BE debería seguir las anotaciones funcionales más que segmentos de ascendencia coherentes en términos de haplotipos.

«El espacio de los grafos de mezcla posibles es enorme, y muchos grafos cualitativamente diferentes pueden ajustarse a los mismos estadísticos f.»
— Maier et al., eLife (2023)

Lo que «es» Basal Eurasian en el argumento original#

En la formulación canónica, BE es un linaje no muestreado que se separó tempranamente de otros no africanos y que posteriormente contribuyó sustancialmente a las poblaciones antiguas del Cercano Oriente; su función empírica clave es explicar por qué los euroasiáticos occidentales muestran una menor afinidad neandertal en relación con los asiáticos orientales y por qué ciertos grupos antiguos del Cercano Oriente se sitúan «fuera» de los árboles euroasiáticos simples (Lazaridis et al. 2016).

Dos detalles son importantes para una crítica steelman:

1. BE no se observa directamente. Se infiere porque añadirlo permite que un grafo de mezcla se ajuste mejor a los patrones de correlaciones entre frecuencias alélicas (estadísticas f).
2. BE está ligado a una estadística de compartición neandertal. Lazaridis et al. destacan explícitamente una correlación negativa entre la ascendencia BE estimada y una estadística destinada a reflejar la ascendencia neandertal (por ejemplo, una forma de (f_4) como (f_4(\text{Test},\text{Mbuti};\text{Altai Neanderthal},\text{Denisovan}))) (Lazaridis et al. 2016).

Incluso en su discusión suplementaria, Lazaridis et al. señalan alternativas plausibles sobre cómo podría haber ingresado al Cercano Oriente una ascendencia «de tipo basal» (por ejemplo, mediante su asociación con un flujo génico procedente de África) o sobre su presencia más temprana en contextos del Levante y Arabia; es decir, la interpretación no queda determinada de manera unívoca ni siquiera dentro del marco favorable a BE (Lazaridis et al. 2016, Supplementary Information).

La formulación steelman comienza aquí: si BE es una variable de compresión para la covarianza residual, entonces cualquier proceso no modelado que cree una covarianza residual similar puede hacerse pasar por BE.

El punto de presión del modelado: los estadísticos f y los grafos de mezcla no pueden representar la selección

Lo que los estadísticos f «suponen» (silenciosamente)#

Para las poblaciones (A,B,C,D), un estadístico f4 básico es (esquemáticamente)

[ f_4(A,B;C,D)=\mathbb{E}[(p_A-p_B)(p_C-p_D)]. ]

Si los cambios en las frecuencias alélicas son impulsados principalmente por la deriva neutral y la mezcla en un grafo, (f_4) presenta invariancias poderosas: bajo la topología adecuada (sin flujo génico a través de ciertas ramas), tiende a ~0; las desviaciones codifican aristas de mezcla. Esta es la base sobre la que se apoya el razonamiento al estilo de qpGraph/qpAdm.

Pero una premisa oculta crucial es la siguiente: la covarianza genómica que se ajusta proviene principalmente de la deriva y la mezcla.

Ahora incorporemos la comprobación metodológica de la realidad: muchos grafos de mezcla distintos pueden ajustarse bien a los mismos estadísticos f; la clase de modelos está subdeterminada, y «un grafo que se ajusta» no es «el grafo que ocurrió» (Maier et al. 2023). Esta no identificabilidad no es una nota al pie: es una invitación a que los términos de sesgo (como la selección) sean «explicados» como ascendencia.

La selección es un término de sesgo que parece deriva#

La selección modifica las frecuencias alélicas. Si la selección es:

• poligénica (muchos loci, efectos pequeños),
• geográficamente estructurada (direcciones y magnitudes diferentes en distintas poblaciones),
• y ligada (arrastrando variación neutral cercana mediante el paisaje de recombinación),

entonces la selección genera diferencias correlacionadas en las frecuencias alélicas entre poblaciones. Es decir, produce exactamente el tipo de estructura de covarianza que los estadísticos f están diseñados para interpretar como deriva y mezcla en un grafo.

Esto no es especulativo: se sabe que la selección ligada confunde la inferencia demográfica e incluso puede provocar una selección incorrecta entre modelos de historias demográficas cuando la selección no se modela (Schrider et al. 2016). Revisiones más amplias subrayan que la selección y la demografía pueden producir patrones similares, lo que dificulta la inferencia conjunta a menos que la selección se modele explícitamente (Johri 2022).

Afirmación steelman: BE es la manera en que un grafo absorbe un término de covarianza inducido por la selección que no fue modelado.

La hipótesis herética (explícitamente): «la selección a favor de la cognición creó una distorsión con forma de BE»

El mundo que postulamos#

Supongamos que, desde hace ~50,000 años —con una intensificación marcada en los últimos ~15,000—, existe selección directa sobre un fenotipo distribuido de «conciencia»: competencia lingüística, teoría de la mente, control ejecutivo, manipulación simbólica, eficiencia del aprendizaje cultural, etc.

Supongamos además que:

1. Es poligénica. (Arquitectura de rasgo de alta dimensionalidad; muchos loci con efectos pequeños.)
2. Está amplificada socialmente y es espacialmente heterogénea. El gradiente de selección depende de la complejidad social, la densidad, las instituciones, las jerarquías de prestigio y la estructura del nicho; características que plausiblemente se intensificaron de manera temprana y repetida en el suroeste de Asia.
3. Interactúa con la carga de introgresión. La selección cognitiva ejerce una fuerte presión sobre redes reguladoras y vías del desarrollo cerebral; los segmentos arcaicos (neandertales) portan una carga elevada de alelos débilmente deletéreos debido al pequeño tamaño efectivo de las poblaciones neandertales, y la selección depura muchos segmentos arcaicos, especialmente cerca de sitios funcionales.

La premisa (3) no es descabellada; ya ocupa un lugar central en los modelos convencionales del agotamiento de la ascendencia neandertal. La selección purificadora contra la introgresión neandertal se explica mejor como una acción sobre muchos alelos débilmente deletéreos, con la selección ligada dando forma al paisaje (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).

También es pertinente que las regiones neandertales supervivientes presentan, en términos generales, un agotamiento de las contribuciones a la heredabilidad de muchos rasgos complejos (lo que sugiere una selección o restricción sistemática), con excepciones como los rasgos de piel y cabello (McArthur et al. 2021).

Cómo esto imita a BE sin ninguna población fantasma#

Queremos recrear dos «motivos» empíricos vinculados que BE explica:

• Motivo 1: Los primeros habitantes del Cercano Oriente parecen «alejarse» de los demás euroasiáticos en el espacio de compartición de alelos, de una manera que el ajuste del grafo captura insertando una separación basal y una arista de mezcla.
• Motivo 2: Su «ascendencia BE» inferida se correlaciona negativamente con las estadísticas de afinidad neandertal.

Este es el mecanismo steelman que genera ambos motivos.

Mecanismo I: la selección añade un término de covarianza estructurada que los grafos interpretan como ascendencia#

Sea el vector de cambio en la frecuencia alélica de la población (X) descompuesto como:

[ \Delta p_X = \Delta p^{\text{drift}}_X + \Delta p^{\text{admixture}}_X + \Delta p^{\text{selection}}_X. ]

Los grafos de mezcla ajustan únicamente los dos primeros términos. Si (\Delta p^{\text{selection}}_X) no es insignificante y está estructurado (por ejemplo, es mayor en los linajes del Cercano Oriente), el grafo ajustado «inventará» deriva adicional o una arista de mezcla fantasma para absorber la covarianza residual.

La no identificabilidad facilita esto: hay muchos grafos capaces de absorber una forma residual determinada y, al mismo tiempo, conservar un buen ajuste a los estadísticos f (Maier et al. 2023).

Mecanismo II: la selección cognitiva acelera la depuración de segmentos neandertales, creando estadísticas neandertales «de tipo basal»#

En este mundo, los linajes de Eurasia occidental —especialmente los del Cercano Oriente— experimentan una selección más intensa sobre un conjunto amplio de loci relevantes para la cognición y el funcionamiento social. Puesto que los segmentos neandertales están enriquecidos en alelos débilmente deletéreos y se encuentran agotados alrededor de regiones funcionales en los humanos actuales, la selección tiende a eliminar la ascendencia neandertal cercana a genes y a la arquitectura reguladora (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).

Así, los linajes del Cercano Oriente pueden terminar con:

• una menor afinidad neandertal en todo el genoma,
• una distribución diferente de los tractos arcaicos entre regiones funcionales y neutrales,
• y una compartición de alelos alterada con otros euroasiáticos, que el ajuste del grafo explica como una «dilución por BE».

Esto genera la correlación negativa característica de BE con las estadísticas neandertales sin requerir un linaje fantasma sin neandertales.

Advertencia (steelman lo reconoce): algunos trabajos sostienen que la selección y la demografía por sí solas no pueden explicar todas las diferencias en la ascendencia neandertal entre los asiáticos orientales y los europeos (Kim & Lohmueller 2015). La respuesta herética es: «de acuerdo; entonces la selección no es el único factor; solo necesita ser lo bastante intensa y estructurada para imitar la forma residual que BE absorbe».

Mecanismo III: lo «basal» también puede producirse mediante dinámicas de retorno a África / dispersión temprana, y la selección asigna el mérito de forma errónea#

Hay evidencia de que las señales neandertales aparentes en África pueden explicarse por migración de retorno y/o flujo génico hacia los neandertales proveniente de los primeros humanos modernos que salieron de África (Chen et al. 2020). Esto importa porque los componentes «relacionados con África» semejantes a BE y los patrones de compartición con neandertales pueden generarse mediante múltiples procesos interactuantes.

Posición steelman: en un mundo dominado por la selección, un grafo puede asignar peso explicativo a una arista «Basal Eurasian» porque es el ajuste más sencillo disponible, aun cuando la historia generativa verdadera sea «(cierto) flujo génico africano + depuración de neandertales variable en el tiempo + selección poligénica estructurada».

Por qué la historia «selección sobre la cognición → BE fantasma» no es solo un juego de palabras#

El contenido serio es que predice distorsiones específicas y falsables en la señal de BE.

Predicción 1: BE debería ser más intensa en particiones sensibles a la selección#

Si BE es un componente de ascendencia, debería aparecer ampliamente en todo el genoma (salvo el ruido de la deriva). Si es un artefacto de selección, debería ser más intensa en las regiones donde la selección ligada es más fuerte:

• baja recombinación,
• proximidad a genes,
• regiones conservadas / densas en elementos reguladores.

Esto es una extensión directa del hecho conocido de que la selección ligada modifica los patrones de diversidad neutra y puede inducir a error en la inferencia demográfica cuando no se modela (Schrider et al. 2016).

Predicción 2: BE debería cambiar a través del tiempo en el Cercano Oriente de una manera semejante a la selección#

Un componente de ascendencia introducido por un acontecimiento antiguo de mezcla debería comportarse (aproximadamente) como una proporción estable, con una deriva lenta salvo que ocurra una mezcla posterior.

Un artefacto de selección debería comportarse más bien como:

• una rampa a lo largo del Holoceno,
• posiblemente con pulsos alineados con grandes transiciones sociales/ecológicas,
• y con efectos más intensos en poblaciones que experimentan una estratificación cultural intensificada y una diferenciación rápida de nichos.

Esto es comprobable porque el ADN antiguo del Cercano Oriente está estratificado temporalmente precisamente en el intervalo en el que la hipótesis sostiene que la selección se intensifica (Lazaridis et al. 2016).

Predicción 3: el carácter «BE» debería correlacionarse más con el enriquecimiento en anotaciones funcionales que con fragmentos de ascendencia coherentes en cuanto a haplotipos#

Una población fantasma implica tractos de ascendencia con patrones coherentes de desequilibrio de ligamiento (LD) y haplotipos (los tractos más antiguos son más cortos, pero siguen teniendo forma de tractos).

Una distorsión por selección poligénica implica:

• muchos cambios pequeños y dispersos en las frecuencias,
• una coherencia haplotípica débil en relación con una arista genuina de ascendencia,
• enriquecimiento en anotaciones relevantes para rasgos (aunque poner a prueba las anotaciones de «cognición» es, en sí mismo, peligroso).

Aquí es donde se debe proceder con extrema cautela: las señales de selección poligénica son notoriamente sensibles a la estratificación y al sesgo de los estudios de asociación del genoma completo (GWAS) (esto se demostró de manera contundente para la estatura) (Sohail et al. 2019). En forma steelman, esta cautela es una característica: la hipótesis predice una señal sesgada hacia determinadas anotaciones, pero también admite que las pruebas ingenuas basadas en GWAS producirán falsos positivos.

Predicción 4: inestabilidad del espacio de modelos: BE aparece porque es un absorbente conveniente#

Si BE es real, su necesidad debería ser robusta en distintos espacios de grafos razonables. Si BE es un absorbente de los residuales inducidos por la selección, entonces:

• muchos grafos diferentes tendrán ajustes similares,
• la posición de BE y su proporción estimada pueden variar según las decisiones de modelado,
• y BE debería ser particularmente sensible al ascertainment de SNP y al filtrado.

Esto concuerda con la advertencia general de que muchos grafos de mezcla cualitativamente distintos pueden ajustarse a los mismos estadísticos f, por lo que la interpretabilidad requiere comprobaciones de robustez, no un único «mejor grafo» (Maier et al. 2023).

Una tabla comparativa concreta: fantasma de ascendencia frente a espejismo de selección#

Cómo poner esto realmente a prueba (sin hacer trampa)#

A continuación se presenta un «diseño experimental» en forma steelman, cuyo objetivo es refutar la hipótesis si es errónea.

Prueba A: estadísticos f particionados y ajuste de grafos#

Calcular los patrones de f4 relevantes para BE (incluidas las estadísticas de afinidad neandertal que utiliza Lazaridis) en:

1. ventanas de alta recombinación y alejadas de los genes,
2. ventanas de baja recombinación / con alta densidad génica,
3. controles emparejados de GC y de mapeabilidad.

Si BE es un espejismo de selección, el tamaño del efecto debería ser dependiente de la partición.

(La razón por la que esperamos una dependencia de la partición se fundamenta en la confusión conocida que la selección ligada introduce en la inferencia demográfica: Schrider et al. 2016.)

Prueba B: pendiente de la serie temporal antigua#

Dentro de linajes del Levante/Irán/Anatolia, a lo largo del Epipaleolítico → Neolítico → Edad del Bronce, estimar si los residuales semejantes a BE aumentan de manera monotónica o puntuada, de un modo que no pueda explicarse parsimoniosamente solo por la mezcla.

Esto utiliza el mismo tipo de conjunto de datos que motivó BE, pero plantea una pregunta diferente: «¿la señal se comporta como ascendencia o como una selección variable en el tiempo?» (Lazaridis et al. 2016).

Prueba C: paisaje de tractos arcaicos frente a la «proporción de BE»#

Si BE es real, la disminución de neandertales debería seguir principalmente las proporciones de ascendencia (dilución). Si BE es un espejismo de selección, predecimos:

• una disminución más intensa cerca de regiones funcionales condicionada a la proporción arcaica general,
• y una divergencia más intensa en las distribuciones de tractos arcaicos en los linajes del Cercano Oriente y de Eurasia occidental.

El paisaje de referencia de la disminución de ascendencia neandertal alrededor de los genes está bien establecido (Sankararaman et al. 2014), y el marco de la selección purificadora se explicita tanto en Juric como en Harris/Nielsen (Juric et al. 2016; Harris & Nielsen 2016).

Prueba D: triangulación entre métodos (comprobación de sensatez de qpAdm, bajas probabilidades previas al estudio)#

La hipótesis predice fragilidad: si BE absorbe sesgo, los flujos de trabajo que examinan de manera agresiva muchos modelos de mezcla pueden mostrar tasas infladas de falsos positivos, a menos que se controlen cuidadosamente.

Evaluaciones recientes del cribado de qpAdm en paisajes complejos hacen hincapié en las preocupaciones relativas al desempeño y al descubrimiento falso cuando el espacio de modelos es grande y las probabilidades previas al estudio son bajas (Flegontova et al. 2025; Williams et al. 2024).

Interpretación steelman: BE podría ser «descubierta» en parte porque el ecosistema analítico recompensa los componentes de los grafos que mejoran el ajuste, aun cuando existan muchos grafos con ajustes equivalentes y la selección sea un proceso no modelado.

Lo que este steelman no afirma#

Para mantener la herejía dentro de límites honestos:

• No afirma que la selección pueda conjurar de la nada una estructura poblacional arbitraria de tiempos profundos.
• No afirma que BE sea falsa; afirma que el fantasma específico de BE puede ser un absorbente conveniente del modelo.
• No se basa en argumentos ingenuos de «selección poligénica para el CI»; estos son empíricamente frágiles y fácilmente confundibles por la estratificación (Sohail et al. 2019).

Afirma algo más estrecho y preciso:

En un mundo con una selección masiva, estructurada y poligénica sobre la cognición, junto con una depuración variable en el tiempo de la carga arcaica, un modelo de grafo de mezcla basado únicamente en la deriva se ve estructuralmente obligado a inventar ascendencia fantasma para dar cuenta de una covariación sistemática que no puede representar.

Preguntas frecuentes#

P 1. Si BE es falsa, ¿por qué aparece de manera «consistente»?
R. Porque el mismo sesgo no modelado (selección + selección ligada + depuración de la carga arcaica) puede ser consistente entre análisis, y los grafos de mezcla no son identificables: muchos grafos pueden ajustarse, y un fantasma basal es un absorbente particularmente eficiente de residuales estructurados (Maier et al. 2023).

P 2. ¿Acaso la correlación neandertal no es la prueba definitiva de un fantasma sin neandertales?

A. No de manera unívoca: la selección purificadora contra los alelos neandertales introgresados está ampliamente distribuida y es en gran medida poligénica, por lo que distintos regímenes de selección pueden alterar las estadísticas de afinidad neandertal sin requerir dilución por parte de un linaje separado (Juric et al. 2016; Harris & Nielsen 2016).

P 3. ¿Qué refutaría rápidamente la perspectiva del espejismo de la selección?
R. Si la señal similar a BE permanece fuerte en particiones intergénicas con alta recombinación y se manifiesta como una ascendencia coherente a nivel de haplotipos, compatible con un evento discreto de mezcla, entonces la hipótesis del espejismo de la selección colapsa; la selección ligada no debería dominar esas regiones casi neutrales (Schrider et al. 2016).

P 4. ¿Cómo afecta el retorno a África/la dispersión temprana a este debate?
R. Añade trayectorias demográficas adicionales que pueden modificar la afinidad neandertal y las señales «relacionadas con África»; si estas están presentes, un grafo ciego a la selección podría atribuir erróneamente sus efectos a un linaje fantasma basal, en lugar de atribuirlos a una demografía más compleja (Chen et al. 2020; Lazaridis et al. 2016, Suplemento).

Notas al pie#

Fuentes#

1. Lazaridis, I., et al. “Genomic insights into the origin of farming in the ancient Near East.” Nature 536 (2016). doi:10.1038/nature19310
2. Lazaridis, I., et al. Supplementary Information for “Genomic insights into the origin of farming in the ancient Near East.” (2016).
3. Maier, R., et al. “On the limits of fitting complex models of population history to f-statistics.” eLife (2023). doi:10.7554/eLife.85492
4. Schrider, D. R., Shanku, A. G., & Kern, A. D. “Effects of linked selective sweeps on demographic inference and model selection.” Genetics (2016). doi:10.1534/genetics.116.190223
5. Johri, P. “On the prospect of achieving accurate joint estimation of population history and the distribution of fitness effects.” Genome Biology and Evolution (2022). doi:10.1093/gbe/evac088
6. Juric, I., Aeschbacher, S., & Coop, G. “The Strength of Selection against Neanderthal Introgression.” PLOS Genetics (2016). doi:10.1371/journal.pgen.1006340
7. Harris, K., & Nielsen, R. “The Genetic Cost of Neanderthal Introgression.” Genetics (2016). doi:10.1534/genetics.116.186890
8. Sankararaman, S., et al. “The genomic landscape of Neanderthal ancestry in present-day humans.” Nature (2014). doi:10.1038/nature12961
9. McArthur, E., et al. “Quantifying the contribution of Neanderthal introgression to the heritability of complex traits.” Nature Communications (2021). doi:10.1038/s41467-021-24582-y
10. Chen, L., et al. “Identifying and Interpreting Apparent Neanderthal Ancestry in African Individuals.” Cell (2020). doi:10.1016/j.cell.2020.01.012
11. Kim, B. Y., & Lohmueller, K. E. “Selection and reduced population size cannot explain higher amounts of Neandertal ancestry in East Asian than in European human populations.” American Journal of Human Genetics (2015). doi:10.1016/j.ajhg.2014.12.029
12. Sohail, M., et al. “Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies.” eLife (2019). doi:10.7554/eLife.39702
13. Flegontova, O., et al. “Performance of qpAdm-based screens for genetic admixture under complex demography and low pre-study odds.” Genetics (2025). doi:10.1093/genetics/iyaf047
14. Williams, M. P., et al. “Testing times: disentangling admixture histories in recent and complex demographic scenarios.” Genetics (2024). doi:10.1093/genetics/iyae110