O Eurasiático Basal como miragem: um argumento steelman em favor da “ancestralidade fantasma” como artefato de seleção

TL;DR

• “Eurasiático Basal” (BE) é um componente de modelo introduzido para satisfazer restrições de correlação de frequências alélicas em grafos de DNA antigo, fortemente associado a um compartilhamento reduzido de alelos neandertais entre antigos habitantes do Oriente Próximo (Lazaridis et al. 2016).
• Grafos de mistura e estatísticas-f pressupõem que a maior parte da covariância em todo o genoma seja gerada por deriva neutra em um grafo; uma seleção forte e geograficamente estruturada pode acrescentar covariância sistemática que o modelo “explica” inventando ramos-fantasma (Maier et al. 2023; Schrider et al. 2016).
• A seleção purificadora contra a introgressão neandertal é real, em grande medida poligênica e mediada por seleção ligada; a depuração diferencial pode alterar as estatísticas de compartilhamento neandertal sem exigir uma linhagem profundamente divergente e livre de neandertais (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).
• Hipótese steelman: seleção maciça no Pleistoceno tardio/Holoceno em favor da cognição nas linhagens do Oriente Próximo/da Eurásia ocidental (poligênica, enriquecida em anotações e estruturada), somada à seleção ligada e à depuração de material arcaico, gera o padrão residual semelhante ao BE, e o ajuste do grafo converte esse viés em um fantasma “Eurasiático Basal”.
• Essa hipótese é testável: o BE deveria diminuir em partições com alta recombinação e intergênicas; o BE deveria apresentar intensidade variável ao longo do tempo em genomas datados do Oriente Próximo; e o sinal semelhante ao BE deveria acompanhar mais as anotações funcionais do que blocos de ancestralidade coerentes em termos de haplótipos.

“O espaço de grafos de mistura possíveis é enorme, e muitos grafos qualitativamente diferentes podem ajustar-se às mesmas estatísticas-f.”
— Maier et al., eLife (2023)

O que o Eurasiático Basal “é” no argumento original#

Na formulação canônica, o BE é uma linhagem não amostrada que se separou cedo das demais populações não africanas e posteriormente contribuiu substancialmente para populações antigas do Oriente Próximo; seu papel empírico central é explicar por que os eurasiáticos ocidentais exibem uma afinidade neandertal reduzida em relação aos asiáticos orientais e por que certos grupos antigos do Oriente Próximo se situam “fora” de árvores eurasiáticas simples (Lazaridis et al. 2016).

Dois detalhes são importantes para uma crítica steelman:

1. O BE não é observado diretamente. Ele é inferido porque sua inclusão faz com que um grafo de mistura se ajuste melhor aos padrões de correlações entre frequências alélicas (estatísticas-f).
2. O BE está vinculado a uma estatística de compartilhamento neandertal. Lazaridis et al. destacam explicitamente uma correlação negativa entre a ancestralidade BE estimada e uma estatística destinada a refletir a ancestralidade neandertal (por exemplo, uma forma de (f_4) tal como (f_4(\text{Test},\text{Mbuti};\text{Altai Neanderthal},\text{Denisovan}))) (Lazaridis et al. 2016).

Mesmo em sua discussão suplementar, Lazaridis et al. observam alternativas plausíveis para a maneira como a ancestralidade “semelhante à basal” poderia ter entrado no Oriente Próximo (por exemplo, sua associação com fluxo gênico proveniente da África) ou para uma presença mais antiga em contextos do Levante/da Arábia; isto é, a interpretação não é determinada de modo unívoco nem mesmo dentro do quadro favorável ao BE (Lazaridis et al. 2016, Supplementary Information).

O steelman começa aqui: se o BE é uma variável de compressão para a covariância residual, então qualquer processo não modelado que crie uma covariância residual semelhante pode fazer-se passar por BE.

O ponto de pressão da modelagem: estatísticas-f e grafos de mistura não podem representar a seleção

O que as estatísticas-f “pressupõem” (silenciosamente)#

Para as populações (A,B,C,D), uma estatística-f4 básica é (esquematicamente)

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

Se as mudanças nas frequências alélicas são impulsionadas principalmente por deriva neutra e mistura em um grafo, (f_4) possui invariâncias poderosas: sob a topologia correta (sem fluxo gênico através de determinados ramos), ela tende a ~0; os desvios codificam arestas de mistura. Esse é o fundamento sobre o qual se apoia o raciocínio em estilo qpGraph/qpAdm.

Mas uma premissa oculta crucial é: a covariância em todo o genoma que está sendo ajustada é composta majoritariamente por deriva + mistura.

Agora, introduza a verificação metodológica da realidade: muitos grafos de mistura distintos podem ajustar-se bem às mesmas estatísticas-f; a classe de modelos é subdeterminada, e “um grafo que se ajusta” não é “o grafo que de fato ocorreu” (Maier et al. 2023). Essa não identificabilidade não é uma nota de rodapé — é um convite para que termos de viés (como a seleção) sejam “explicados” como ancestralidade.

A seleção é um termo de viés que se parece com deriva#

A seleção altera as frequências alélicas. Se a seleção é:

• poligênica (muitos loci, efeitos pequenos),
• geograficamente estruturada (direções e magnitudes diferentes em populações diferentes),
• e ligada (arrastando variação neutra próxima por meio da paisagem de recombinação),

então a seleção gera diferenças correlacionadas nas frequências alélicas entre populações. Isto é, ela produz exatamente o tipo de estrutura de covariância que as estatísticas-f foram projetadas para interpretar como deriva e mistura em um grafo.

Isso não é especulativo: sabe-se que a seleção ligada confunde a inferência demográfica, chegando até a provocar escolha incorreta de modelos entre histórias demográficas quando a seleção não é modelada (Schrider et al. 2016). Revisões mais abrangentes enfatizam que seleção e demografia podem produzir padrões semelhantes, tornando difícil a inferência conjunta quando a seleção não é modelada explicitamente (Johri 2022).

Afirmação steelman: o BE é a maneira pela qual um grafo absorve um termo de covariância induzido por seleção não modelada.

A hipótese herética (explicitamente): “a seleção pela cognição criou uma distorção com a forma do BE”

O mundo que postulamos#

Suponha que, a partir de ~50.000 anos atrás — intensificando-se acentuadamente nos últimos ~15.000 — tenha havido seleção direta sobre um fenótipo distribuído de “consciência”: proficiência linguística, teoria da mente, controle executivo, manipulação simbólica, eficiência da aprendizagem cultural etc.

Suponha, além disso, que:

1. Ela seja poligênica. (Arquitetura fenotípica de alta dimensionalidade; muitos loci de pequeno efeito.)
2. Ela seja socialmente amplificada e espacialmente heterogênea. O gradiente de seleção depende da complexidade social, da densidade, das instituições, das hierarquias de prestígio e da estrutura de nichos — características que plausivelmente se intensificaram cedo e repetidamente no Sudoeste Asiático.
3. Ela interaja com a carga de introgressão. A seleção pela cognição incide fortemente sobre redes regulatórias e vias do desenvolvimento cerebral; segmentos arcaicos (neandertais) carregam uma carga elevada de alelos fracamente deletérios em razão do pequeno tamanho efetivo da população neandertal, e a seleção depura muitos segmentos arcaicos, especialmente nas proximidades de sítios funcionais.

A premissa (3) não é absurda; ela já é central nos modelos correntes de depleção da ancestralidade neandertal. A seleção purificadora contra a introgressão neandertal é mais bem explicada como atuando sobre muitos alelos fracamente deletérios, com a seleção ligada moldando a paisagem (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).

Também é relevante que as regiões neandertais sobreviventes apresentem, de modo geral, depleção das contribuições para a herdabilidade de muitos traços complexos (sugerindo seleção/restrição sistemática), com exceções como os traços de pele/cabelo (McArthur et al. 2021).

Como isso contrafaz o BE sem qualquer população-fantasma#

Queremos recriar dois “motivos” empíricos relacionados que o BE explica:

• Motivo 1: Os antigos habitantes do Oriente Próximo parecem “afastados” dos demais eurasiáticos no espaço de compartilhamento alélico, de uma maneira que o ajuste do grafo captura mediante a inserção de uma separação basal e de uma aresta de mistura.
• Motivo 2: A sua “ancestralidade BE” inferida apresenta correlação negativa com estatísticas de afinidade neandertal.

Eis o mecanismo steelman que gera ambos os motivos.

Mecanismo I: a seleção acrescenta um termo de covariância estruturada que os grafos interpretam como ancestralidade#

Decomponha o vetor de mudança na frequência alélica da população (X) da seguinte maneira:

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

Os grafos de mistura ajustam apenas os dois primeiros termos. Se (\Delta p^{\text{selection}}_X) for não desprezível e estruturado (por exemplo, maior nas linhagens do Oriente Próximo), o grafo ajustado “inventará” deriva adicional e/ou uma aresta de mistura fantasma para absorver a covariância residual.

A não identificabilidade facilita isso: há muitos grafos capazes de absorver uma dada forma residual preservando o ajuste às estatísticas-f (Maier et al. 2023).

Mecanismo II: a seleção pela cognição acelera a depuração de segmentos neandertais, criando estatísticas neandertais “semelhantes às basais”#

Nesse mundo, as linhagens da Eurásia ocidental (especialmente as do Oriente Próximo) sofrem uma seleção mais forte sobre um amplo conjunto de loci relevantes para a cognição e o funcionamento social. Como os segmentos neandertais são enriquecidos em alelos fracamente deletérios e se encontram reduzidos nas proximidades de regiões funcionais nos humanos atuais, a seleção tende a remover a ancestralidade neandertal próxima a genes e à arquitetura regulatória (Juric et al. 2016; Harris & Nielsen 2016; Sankararaman et al. 2014).

Assim, as linhagens do Oriente Próximo podem acabar apresentando:

• menor afinidade neandertal em todo o genoma,
• uma distribuição diferente de tratos arcaicos entre regiões funcionais e neutras,
• e um compartilhamento alélico alterado com outros eurasiáticos, que o ajuste do grafo explica como uma “diluição pelo BE”.

Isso gera a correlação negativa característica de BE com as estatísticas neandertais sem exigir uma linhagem-fantasma sem neandertais.

Ressalva (o steelman a reconhece): alguns trabalhos argumentam que seleção + demografia, por si sós, não podem explicar todas as diferenças na ancestralidade neandertal entre asiáticos orientais e europeus (Kim & Lohmueller 2015). A resposta herética é: “tudo bem — então a seleção não é o único fator; ela apenas precisa ser suficientemente forte e estruturada para contrafazer a forma residual que BE absorve.”

Mecanismo III: “Basal” também pode ser produzido por dinâmicas de retorno à África / dispersão inicial, e a seleção atribui o crédito de modo equivocado#

Há evidências de que os sinais neandertais aparentes na África podem ser explicados por migração de retorno e/ou fluxo gênico para os neandertais oriundo de humanos modernos iniciais provenientes da África (Chen et al. 2020). Isso importa porque componentes “relacionados à África” semelhantes a BE e padrões de compartilhamento com neandertais podem ser gerados por múltiplos processos interativos.

Posição steelman: em um mundo dominado pela seleção, um grafo pode atribuir peso explicativo a uma aresta “Eurasiático Basal” porque ela é o parâmetro mais simples disponível, mesmo que a verdadeira história geradora seja “(algum) fluxo gênico africano + depuração variável no tempo de ancestralidade neandertal + seleção poligênica estruturada”.

Por que a história “seleção sobre cognição → BE fantasma” não é apenas um jogo de palavras#

O conteúdo sério é que ela prevê distorções específicas e falseáveis no sinal de BE.

Previsão 1: BE deve ser mais forte em partições sensíveis à seleção#

Se BE é um componente de ancestralidade, ele deve aparecer amplamente pelo genoma (modulo o ruído da deriva). Se é um artefato de seleção, deve ser mais forte em regiões onde a seleção ligada é mais intensa:

• baixa recombinação;
• proximidade de genes;
• regiões conservadas / densas em elementos regulatórios.

Isso é uma extensão direta do fato conhecido de que a seleção ligada altera os padrões de diversidade neutra e pode induzir inferências demográficas equivocadas quando não é modelada (Schrider et al. 2016).

Previsão 2: BE deve mudar ao longo do tempo no Oriente Próximo de uma maneira semelhante à seleção#

Um componente de ancestralidade introduzido por um evento de mistura antigo deveria comportar-se (aproximadamente) como uma proporção estável, derivando lentamente, a menos que ocorra mistura posterior.

Um artefato de seleção deveria comportar-se mais como:

• uma ascensão gradual ao longo do Holoceno;
• possivelmente com pulsos alinhados a grandes transições sociais/ecológicas;
• e com efeitos mais fortes em populações que experimentam estratificação cultural intensificada e rápida diferenciação de nichos.

Isso é testável porque o DNA antigo do Oriente Próximo é estratificado temporalmente exatamente no intervalo em que a hipótese afirma que a seleção se intensifica (Lazaridis et al. 2016).

Previsão 3: o “caráter BE” deve acompanhar mais o enriquecimento de anotações funcionais do que blocos de ancestralidade coerentes em termos de haplótipos#

Uma população-fantasma implica tratos de ancestralidade com padrões coerentes de desequilíbrio de ligação (LD) e de haplótipos (tratos mais antigos são mais curtos, mas ainda assim mantêm o caráter de tratos).

Uma distorção decorrente de seleção poligênica implica:

• muitas mudanças pequenas e dispersas nas frequências;
• coerência haplotípica fraca em relação a uma aresta genuína de ancestralidade;
• enriquecimento em anotações relevantes para traços (embora testar anotações de “cognição” seja, por si só, traiçoeiro).

É nesse ponto que se deve ser extremamente cauteloso: os sinais de seleção poligênica são notoriamente sensíveis à estratificação e ao viés de GWAS (demonstrado de modo contundente para a altura) (Sohail et al. 2019). Em steelman, essa cautela é uma característica: a hipótese prevê um sinal enviesado em favor de determinadas anotações, mas também admite que testes ingênuos baseados em GWAS produzirão falsos positivos.

Previsão 4: instabilidade do espaço de modelos — BE aparece porque é um absorvedor conveniente#

Se BE é real, sua necessidade deve ser robusta em diferentes espaços de grafos razoáveis. Se BE é um absorvedor de resíduos induzidos pela seleção, então:

• muitos grafos diferentes se ajustarão de maneira semelhante;
• a posição de BE e sua proporção estimada poderão variar conforme as escolhas de modelagem;
• e BE deverá ser particularmente sensível ao ascertainment de SNPs e à filtragem.

Isso está de acordo com o alerta geral de que muitos grafos de mistura qualitativamente diferentes podem ajustar-se às mesmas estatísticas f, de modo que a interpretabilidade exige verificações de robustez, não um único “melhor grafo” (Maier et al. 2023).

Uma tabela comparativa concreta: ancestralidade-fantasma versus miragem de seleção#

Como testar isso de fato (sem trapacear)#

A seguir, apresenta-se um “desenho experimental” em steelman, destinado a refutar a hipótese caso ela esteja errada.

Teste A: estatísticas f particionadas e ajuste de grafos#

Calcular os padrões de f4 relevantes para BE (incluindo as estatísticas de afinidade neandertal usadas por Lazaridis) em:

1. janelas de alta recombinação, distantes de genes;
2. janelas de baixa recombinação / densas em regiões gênicas;
3. controles pareados de GC e mapeabilidade.

Se BE é uma miragem de seleção, o tamanho do efeito deve ser dependente da partição.

(O motivo pelo qual esperamos dependência da partição fundamenta-se no confundimento conhecido da inferência demográfica pela seleção ligada: Schrider et al. 2016.)

Teste B: inclinação da série temporal antiga#

Dentro de linhagens do Levante/Irã/Anatólia, abrangendo o Epipaleolítico → Neolítico → Idade do Bronze, estimar se os resíduos semelhantes a BE aumentam de modo monotônico ou pontuado, de uma maneira que não seja parcimoniosamente explicada apenas por mistura.

Isso utiliza o mesmo tipo de conjunto de dados que motivou BE, mas formula uma pergunta diferente: “o sinal comporta-se como ancestralidade ou como seleção variável no tempo?” (Lazaridis et al. 2016).

Teste C: paisagem de tratos arcaicos versus “proporção de BE”#

Se BE é real, a redução neandertal deve seguir principalmente as proporções de ancestralidade (diluição). Se BE é uma miragem de seleção, prevemos:

• redução mais forte nas proximidades de regiões funcionais condicionalmente à proporção arcaica geral;
• e divergência mais forte nas distribuições de tratos arcaicos em linhagens do Oriente Próximo/da Eurásia Ocidental.

A paisagem basal da redução da ancestralidade neandertal ao redor de genes está bem estabelecida (Sankararaman et al. 2014), e a formulação em termos de seleção purificadora é explícita tanto em Juric quanto em Harris/Nielsen (Juric et al. 2016; Harris & Nielsen 2016).

Teste D: triangulação entre métodos concorrentes (verificação de sanidade do qpAdm, baixas probabilidades a priori pré-estudo)#

A hipótese prevê fragilidade: se BE está absorvendo viés, então pipelines que examinam agressivamente muitos modelos de mistura podem apresentar falsos positivos inflados, a menos que sejam cuidadosamente controlados.

Avaliações recentes da triagem com qpAdm em paisagens complexas enfatizam preocupações relativas ao desempenho e à falsa descoberta quando o espaço de modelos é amplo e as probabilidades a priori pré-estudo são baixas (Flegontova et al. 2025; Williams et al. 2024).

Interpretação steelman: BE pode ser “descoberto” em parte porque o ecossistema analítico recompensa componentes de grafos que melhoram o ajuste, mesmo quando existem muitos grafos com ajuste equivalente e a seleção é um processo não modelado.

O que este steelman não afirma#

Para manter a heresia honesta:

• Ela não afirma que a seleção pode conjurar, do nada, qualquer estrutura populacional de tempos profundos.
• Ela não afirma que BE é falso; afirma que a população-fantasma específica de BE pode ser um absorvedor conveniente do modelo.
• Ela não depende de argumentos ingênuos de “seleção poligênica para o QI”; esses argumentos são empiricamente frágeis e facilmente confundidos pela estratificação (Sohail et al. 2019).

Ela afirma algo mais restrito e mais incisivo:

Em um mundo com seleção maciça, estruturada e poligênica sobre a cognição, combinada à depuração variável no tempo da carga arcaica, um modelo de grafo de mistura baseado apenas na deriva é estruturalmente forçado a inventar ancestralidade-fantasma para explicar uma covariância sistemática que não pode representar.

Perguntas frequentes#

P 1. Se BE é falso, por que aparece de modo “consistente”?
R. Porque o mesmo viés não modelado (seleção + seleção ligada + depuração de ancestralidade arcaica) pode ser consistente entre as análises, e os grafos de mistura são não identificáveis — muitos grafos se ajustam, e uma população-fantasma basal é um absorvedor particularmente eficiente de resíduos estruturados (Maier et al. 2023).

P 2. A correlação neandertal não é a prova definitiva de uma população-fantasma sem neandertais?

R. Não de maneira exclusiva: a seleção purificadora contra alelos neandertais introgressados é disseminada e, em grande medida, poligênica, de modo que regimes de seleção diferenciais podem alterar estatísticas de afinidade neandertal sem exigir diluição por uma linhagem distinta (Juric et al. 2016; Harris & Nielsen 2016).

P 3. O que falsificaria rapidamente a perspectiva da miragem da seleção?
R. Se o sinal semelhante ao BE permanecer forte em partições intergênicas com alta recombinação e se manifestar como ancestralidade coerente em haplótipos, consistente com um evento discreto de mistura, então a hipótese da miragem da seleção colapsa; a seleção ligada não deveria dominar essas regiões quase neutras (Schrider et al. 2016).

P 4. Como o retorno à África/a dispersão inicial afeta esse debate?
R. Isso acrescenta trajetórias demográficas adicionais que podem alterar a afinidade neandertal e os sinais “relacionados à África”; se elas estiverem presentes, um grafo cego à seleção poderá atribuir equivocadamente seus efeitos a um fantasma basal, em vez de a uma demografia mais complexa (Chen et al. 2020; Lazaridis et al. 2016, Supplement).

Notas de rodapé#

Fontes#

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