Genomic impact of the second plague pandemic on three human populations

Resumen Técnico: Impacto Genómico de la Segunda Pandemia de Peste en Tres Poblaciones Humanas

Planteamiento del Problema
La segunda pandemia de peste (principios del siglo XIV–principios del siglo XIX), causada por Yersinia pestis, provocó un colapso demográfico masivo en toda Eurasia, con tasas de mortalidad estimadas entre el 40 % y el 60 % en muchas regiones europeas. Mientras que los registros históricos documentan cambios socioeconómicos profundos, las consecuencias genéticas de esta pandemia siguen siendo objeto de debate. Estudios previos han propuesto que la pandemia impulsó una fuerte selección positiva en loci específicos relacionados con la inmunidad (p. ej., ERAP2, CCR5, genes TLR) y alteró los patrones de migración. Sin embargo, las investigaciones previas de ADN antiguo (aDNA) se han visto limitadas por tamaños de muestra pequeños, un enfoque en un solo sitio o la dependencia de la captura dirigida de regiones genómicas específicas en lugar de la secuenciación del genoma completo. Este estudio aborda la necesidad de un conjunto de datos más grande, con mayor cobertura y geográficamente diverso para probar rigurosamente la selección impulsada por la peste y cuantificar los cambios demográficos en la composición de la ascendencia antes y después de la llegada de la pandemia.

Metodología
Los autores analizaron un conjunto de datos de 529 genomas completos antiguos de alta calidad (profundidad media de 8.8x) de tres sitios urbanos del norte de Europa: Lund (Suecia, N=337), Trondheim (Noruega, N=135) y Vilnius (Lituania, N=57). Los individuos fueron estratificados temporalmente en grupos "pre-Peste Negra" (pre-BD) y "post-Peste Negra" (post-BD) basándose en el contexto arqueológico y la datación por radiocarbono.

El flujo de trabajo analítico incluyó:

1. Control de Calidad e Imputación: Tras filtrar por baja cobertura y contaminación, los genotipos fueron imputados mediante GLIMPSE2 utilizando un panel de referencia de 50,129 individuos contemporáneos (deCODE 50k), lo que resultó en aproximadamente 99.6 millones de loci.
2. Escaneos de Selección: Se realizaron estudios de asociación de genoma completo (GWAS) utilizando modelos lineales mixtos (GEMMA) para detectar diferencias en la frecuencia alélica entre los grupos pre-BD y post-BD. Las covariables incluyeron el sexo y los primeros siete componentes principales (PCs) de un PCA de Europa Occidental para controlar la estructura poblacional. Un metaanálisis combinó los resultados de Lund y Trondheim para aumentar la potencia estadística.
3. Verificación de Replicación: El estudio probó específicamente las frecuencias alélicas en loci reportados previamente como bajo selección (p. ej., ERAP2, TICAM2, NFATC1, CTLA4, RIPK2, CCR5-delta32 y variantes del clúster TLR).
4. Análisis de Ascendencia y Demografía:
   • PCA y Dispersión: Los individuos antiguos fueron proyectados en un PCA de 9,052 personas actuales de Europa Occidental. Se calcularon la dispersión multivariante (varianza generalizada) y las distancias de Mahalanobis ($D_M$) para evaluar la diversidad de la ascendencia y la proximidad a los proxies regionales modernos.
   • Identidad por Descendencia (IBD): Se calculó el intercambio de IBD (>6 cM para Escandinavia, >7 cM para Lituania) entre individuos antiguos y grandes paneles de referencia actuales (OmniExpress para Escandinavia; UK Biobank para el este de Europa/Báltico) para identificar valores atípicos de ascendencia y orígenes locales frente a no locales.
   • Clasificación de Valores Atípicos: Se utilizó un marco de dos pasos que combina el PCA basado en $D_M$ y el intercambio de IBD para clasificar a los individuos como "valores atípicos" (ascendencia no local) y subdividir además a los locales en "locales" y "no locales" dentro de Escandinavia.

Resultados Clave

1. Ausencia de una Fuerte Selección Positiva: El estudio no encontró evidencia de una fuerte selección positiva impulsada por la exposición a Y. pestis que fuera estadísticamente significativa a nivel genómico.
   • Ninguna variante cumplió con los umbrales de significancia genómica informados funcionalmente en el metaanálisis de Lund y Trondheim.
   • Los loci candidatos reportados previamente (incluido ERAP2) no mostraron cambios significativos en la frecuencia alélica entre los grupos pre-BD y post-BD.
   • El análisis de potencia indicó que el estudio tenía suficiente potencia (>80 %) para detectar cambios de magnitud en la frecuencia reportados en estudios previos (p. ej., para ERAP2), lo que sugiere que si tal selección fuerte ocurrió, habría sido detectada. Los autores señalan que cambios sutiles en la frecuencia (<0.12–0.26) o la selección en variantes raras requerirían tamaños de muestra aún mayores para ser detectados.
2. Reducción en la Diversidad de la Ascendencia: En contraste con la falta de señales de selección, el estudio identificó un declive marcado en la diversidad de la ascendencia en los tres sitios tras la llegada de la Peste Negra.
   • Dispersión del PCA: Los individuos pre-BD exhibieron una mayor dispersión en el espacio del PCA en comparación con los individuos post-BD (prueba de Box's M: $P < 2.2 \times 10^{-16}$ para Lund).
   • Declive de Valores Atípicos: La proporción de valores atípicos de ascendencia (individuos con ascendencia no local) disminuyó significativamente de pre-BD a post-BD. En Lund, los valores atípicos disminuyeron del 52 % al 9 % ($P = 1.37 \times 10^{-13}$); en Vilnius, del 47 % al 5 % ($P = 3.71 \times 10^{-4}$).
   • Cambios Específicos de la Ascendencia: El Lund pre-BD mostraba una alta diversidad con contribuciones significativas de las Islas Británicas, Irlanda, Alemania y la región de Benelux. El post-BD, la población se volvió más homogénea, agrupándose más cerca de las poblaciones actuales de Escandinavia. En Trondheim, los valores atípicos pre-BD incluían un número notable de individuos con ascendencia islandesa, la cual declinó tras la Peste Negra.
   • Local vs. No Local: Mientras que la inmigración de larga distancia (valores atípicos) disminuyó, la proporción de escandinavos no locales (migración dentro de la región) se mantuvo estadísticamente sin cambios entre periodos.

Significancia y Afirmaciones
El artículo afirma que, si bien la segunda pandemia de peste fue una catástrofe demográfica, no dejó una firma detectable de fuerte selección positiva en los genomas humanos de estas poblaciones del norte de Europa, desafiando la hipótesis de que la pandemia actuó como un motor primario de evolución adaptativa rápida en los genes inmunitarios.

En su lugar, el estudio destaca que la pandemia coincidió con una transición sociohistórica más amplia hacia sociedades más insulares. Los autores atribuyen la reducción observada en la migración de larga distancia y la diversidad de la ascendencia a una confluencia de factores:

• El fin de la Era Vikinga.
• El inicio de la Pequeña Edad de Hielo.
• El proceso de cristianización, que inicialmente atrajo a diversos misioneros y emisarios (contribuyendo a la diversidad pre-BD) pero que pudo haber estabilizado posteriormente a las poblaciones locales.
• La propia pandemia, que probablemente interrumpió las redes comerciales, eclesiásticas y políticas que facilitaban la movilidad a larga distancia.

Los autores concluyen que el impacto genético de la pandemia en estas regiones fue primordialmente demográfico (reduciendo la migración y la diversidad) en lugar de selectivo (alterando las frecuencias de variantes protectoras). Subrayan que detectar señales de selección sutiles requeriría tamaños de muestra sustancialmente mayores o estrategias de muestreo que se centren en individuos que murieron durante la fase aguda de la pandemia frente a los supervivientes.