Acumulativamente, el extenso número de muestras reunidas en este documento, combinado con la disponibilidad de filogenias paternas muy resueltas que se remontan a muchas generaciones, han permitido que ciertas inferencias sobre el haplogrupo levita Asquenazí R1a linaje fundador se formulen con un alto grado de confianza. Un total de 71 individuos que declararon una herencia paterna ashkenazi Levita o ashkenazi no Levita se determinó al inicio del estudio que pertenecían al haplogrupo R1a-M582, y en base a sus perfiles STR, que representaban un amplio rango de la variación encontrada dentro de ese haplogrupo. Sorprendentemente, todos los cromosomas 71 Y podrían reasignarse a un solo clado en expansión anidado dentro de R1a-M582 y etiquetado en este documento R1a-Y2619, siguiendo una de las seis variantes encontradas en ese nivel (Fig. 2 y Figura complementaria S1). Las únicas muestras notificadas anteriormente que compartían este clado eran cinco muestras de herencia judía asquenazí para las que no se informó el estatus de levita26. En consecuencia, todos los individuos R1a-Y2619, ya sean judíos o no judíos, asquenazíes o no Asquenazíes, Levitas o no Levitas, son los descendientes masculinos directos de la línea paterna de un ancestro masculino común que vivió ~1,743 ybp (Tabla 1). Como los machos contemporáneos de todas las ramas de R1a-Y2619 muestreadas hasta ahora llevan uno de los muchos apellidos de los levitas, se puede argumentar fuertemente que este antepasado masculino auto-afiliado como levita y puede haber llevado el apellido patronímico Levita. La magnitud del evento fundacional se puede estimar calculando el número de individuos contemporáneos que se espera porten un cromosoma Y R1a-Y2619. Dos grupos de muestras independientes de judíos asquenazíes en los que se desconocía el estatus de levita han estimado de manera similar el porcentaje del haplogrupo paterno R1a-Y2619, R1a-M17/M198, en la población asquenazí en 9.6% 20 y 11,5%, respectivamente27. El documento anterior informó que el haplogrupo R1a-M582 representaba el 7,9% de la población asquenazí20. Aquí, mostramos que todas las muestras ashkenazi pertenecientes al haplogrupo R1a-M582 pueden reclasificarse como R1a-Y2619. Sobre la base de un tamaño de población ashkenazi de ~4,000,00013 hombres, de los cuales aproximadamente el 7,9% son R1a-Y2619, habría ~300,000 hombres ashkenazíes descendiendo en su línea masculina directa de un solo ancestro relativamente reciente, con muchos de esos hombres auto-afiliados como levitas.
Dos factores nos permitieron lograr edades de coalescencia e intervalos de confianza más precisos que los calculados previamente para la línea de levita asquenazí R1a18,20. En primer lugar, el gran número de cromosomas Y completos disponibles (Tabla Suplementaria S1), y en segundo lugar, los datos disponibles del pedigrí de Horowitz que permitieron un punto de calibración interno (Fig. 3). Los registros genealógicos bien documentados del pedigrí de Horowitz demuestran cuán robustas pueden ser las secuencias completas del cromosoma Y en la reconstrucción de los ancestros paternos. El punto de calibración se basó en tres muestras de Horowitz con genealogías bien documentadas que fusionaron 402 ybp. Dos individuos adicionales, auto-afiliados con la dinastía Horowitz, pero sin genealogías claras agrupadas estrechamente con estas tres muestras. Ninguna de las muestras de los clados hermanos inmediatos se identificó como parte de la dinastía Horowitz. Sorprendentemente, la edad molecular del nodo que comprende las cinco muestras de Horowitz se calculó en ~690 ybp, mientras que el tiempo de coalescencia genealógica sugiere que es de 546 ybp. Esto sugiere que el nodo de Horowitz R1a-YP268 (g.23133909G > A) podría representar la secuencia reconstruida real del fundador nombrado de este pedigrí (Fig. 3b). Sin embargo, es importante tener en cuenta que las edades de coalescencia estimadas son aproximadas, ya que algunas variaciones metodológicas pueden conducir a resultados algo diferentes. Por ejemplo, la aplicación de nuestra metodología de coalescencia a los datos obtenidos solo de la plataforma Illumina, habría producido una edad de coalescencia de ~1,398 ybp para el clado de Levita Ashkenazi R1a-Y2619 y ~645 años ybp para el clado R1a-YP268. Estas diferencias podrían atribuirse a la capacidad de utilizar tramos superpuestos más largos del cromosoma Y al comparar muestras secuenciadas por la misma maquinaria y parámetros. No calculamos un tiempo de coalescencia basado en las variaciones de STR, ya que hemos demostrado previamente que no proporcionan más información para los vestidos recién fusionados6.
El origen propuesto en Oriente Medio del linaje levita asquenazí basado en lo que antes era un número relativamente limitado de muestras notificadas20, ahora se puede considerar firmemente validado. Si bien las frecuencias más altas del haplogrupo R1a se encuentran en Europa Oriental18,28, nuestros datos revelaron una rica variación del haplogrupo R1a fuera de Europa, que está filogenéticamente separado de las ramas típicamente europeas de R1a. Evidentemente, el R1a-Y2619 está bien anidado dentro de una plétora de clados hermanos filogenéticamente cercanos de Oriente Medio, muestreados en azeríes iraníes, un Kerman, un Yazidi y un hombre de Iberia. Esto proporciona la evidencia necesaria para su origen (Fig. 2 y Figura complementaria S1). Sin embargo, la ruta de migración exacta del R1a-Y2619 a Europa sigue siendo difícil de alcanzar. La mayoría de los registros históricos sugieren dos rutas principales de migración judía a Europa (Fig. 1)29,30. Se considera que la judería asquenazí fue fundada como resultado de la migración judía a través de Italia al Valle del Rin, y luego a Polonia. Se considera que los judíos sefardíes (españoles) emigraron junto con la expansión islámica gradual al norte de África y luego a España14. Debido a que el R1a-Y2619 es de origen de Oriente Medio, es posible que su introducción a Europa haya sido por una o ambas rutas. Naturalmente, el fuerte evento fundacional de R1a-Y2619 entre los judíos asquenazíes, junto con la presencia de todas las ramas conocidas de R1a-Y2619 en los judíos asquenazíes, tienta a inferir que su ruta de migración desde el Levante estaba directamente relacionada con los fundadores asquenazíes. Sin embargo, estos hechos podrían reflejar simplemente una expansión dentro de los judíos asquenazíes en lugar de una prueba de la primera llegada de los R1a-Y2619 con o dentro de la población asquenazí. Más confuso es el hecho de que los registros genealógicos de la dinastía rabínica Horowitz, que ahora se muestra que lleva el cromosoma Y R1a-Y2619, sugieren su presencia en la Península Ibérica en el siglo XV y probablemente antes (Fig. 3)17. De hecho, las migraciones judías repetidas que podrían haber llevado cromosomas R1a-Y2619 Y a Cataluña se documentan desde el siglo IV y durante la expansión musulmana a Iberia31. Además, debido a que Cataluña volvió a ser territorio cristiano desde el año 800 d. C., los antepasados proto-Horowitz R1a-Y2619 también podrían representar la migración de judíos asquenazíes a Iberia. En consecuencia, la presencia de R1a-Y2619 en España en el siglo XV no pudo establecer pruebas de la primera llegada del linaje R1a-Y2619 a la Península Ibérica, ya que esto podría reflejar simplemente movimientos repetidos y no registrados de judíos de ida y vuelta entre Europa Oriental y Occidental y la Península Ibérica. Evidencias previas de ADNmt y marcadores autosómicos ya han sugerido un probable flujo genético entre judíos asquenazíes y sefardíes dentro de Europa32, 33, 34. En conjunto, nuestros resultados tienden a favorecer una única ruta de entrada a Europa como parte de la migración y expansión ashkenazi en Europa. Debido a que el tiempo de coalescencia de todos los Levitas contemporáneos R1a-Y2619 es de ~1,743 ybp, bien dentro del tiempo de la Diáspora del exilio romano, y cada una de las ramas de R1a-Y2619 se encuentra en judíos asquenazíes, nuestros resultados son inconsistentes con un escenario de rápida expansión en el Levante seguido de una propagación a través de múltiples rutas a Europa. Los resultados de las muestras de levita no asquenazí R1a-Y2619 también sugieren una ruta de expansión única. Si el linaje levita R1a-Y2619 hubiera sido tan prevalente en el Levante, se hubiera esperado que los Levitas no asquenazíes se separaran de las muestras de Levitas asquenazíes antes de la Diáspora y formaran un grupo diferente. Sin embargo, los resultados muestran que los Levitas no asquenazíes caen dentro de múltiples sub-ramas coalescentes relativamente recientes de R1a-Y2619 (Fig. 2 y Figura complementaria S1). Este patrón es más compatible con el flujo continuo de genes de la población asquenazí a la población no asquenazí durante la Diáspora, en lugar de múltiples rutas de entrada para el haplogrupo R1a-Y2619 a Europa desde el Levante.
El apoyo adicional para una única ruta de migración y la posterior expansión dentro de los judíos asquenazíes, emana de los patrones observados de expansión de haplogrupos asquenazíes adicionales. En consecuencia, es importante distinguir entre un evento fundacional específico de Levitas asquenazíes para el haplogrupo R1a-Y2619 que podría haber sido el resultado de un estatus socioeconómico favorable o de otro tipo, y una expansión general de los linajes fundacionales del cromosoma Y asquenazíes, incluidos los Levitas asquenazíes R1a-Y2619. Para este propósito, hemos estudiado un puñado de otros haplogrupos descritos anteriormente como prevalentes entre los judíos Asquenazí19,24,35. Nuestros resultados muestran que la expansión de R1a-Y2619 entre los judíos asquenazíes no es específica ni exclusiva de este haplogrupo. Por ejemplo, la fusión del haplogrupo G-M377 (Figura Suplementaria S4) y Q-M242 (Figura Suplementaria S5), que se sabe que representan cada uno el 5% de la variación paterna ashkenazi 24, se fusionan a ~1.223 y ~1.672 años, respectivamente (Tabla 1). La estrategia del cromosoma Y adoptada aquí nos permitió resolver los linajes del haplogrupo E en sus diminutas sub-ramas (Figura Suplementaria S3). Las edades de coalescencia del haplogrupo E-Z838, E – PF3780, E – B923 y E-B933, conocidas por representar acumulativamente el 20% de la variación paterna ashkenazí24, se estimaron en ~1.200–1.600 años. Este patrón de múltiples eventos fundacionales, no observado entre los judíos Españoles36, proporciona una base adicional de que el ancestro levita asquenazí R1a-Y2619 entró en Europa a través de la ruta asquenazí en lugar de a través de la expansión judía a la Península Ibérica. Otros patrones también son claramente visibles en la ascendencia paterna judía asquenazí. El Haplogrupo T-M70, frecuente en Oriente Medio, también está presente entre los judíos Asquenazíes (Figura complementaria S6). Las muestras genotipadas mostraron fracturas profundamente arraigadas que probablemente apuntan a la preservación de una diversidad antigua de este haplogrupo en el Levante, que se remonta al Pleistoceno. Mientras tanto, el patrón observado para el haplogrupo R1b-M269 (Figura Suplementaria S7), prevalente en Europa occidental, sugiere principalmente una introgresión repetida puntuada de los cromosomas Y europeos a la comunidad ashkenazi, y es compatible con informes anteriores36.
Además, comparamos el linaje fundador más frecuente encontrado entre los machos Cohen asquenazíes, anidados dentro del haplogrupo J1a-P58, con el linaje levita asquenazí R1a-Y2619. Evidentemente, los miembros de la casta levita R1a-Y2619 y la casta cohen J1a-P58 no comparten un ancestro masculino común dentro del marco de tiempo de la narración bíblica. Como en el caso de los Levitas Asquenazíes, el linaje de Cohen Asquenazí J1a-P58 formó un cúmulo apretado anidado dentro de un conjunto de muestras de Oriente Medio que confirma su origen, y fue compartido en nuestro estudio solo con Cohens no asquenazíes (Figura Suplementaria S2). Sin embargo, a diferencia del patrón obtenido para los Levitas asquenazíes R1a-Y2619, el grupo se une ~2,570 yp, lo que apunta al inicio de su expansión en el período previo a la Diáspora. Secuencias extensas de cromosomas Y completos a partir de un mayor número de muestras de Cohen de más comunidades judías en el contexto de otras poblaciones levantinas, incluida la determinación de variantes específicas de familias y clanes, serían muy informativas para abordar el hallazgo de que los individuos de Cohen ashkenazi y no Ashkenazi comparten una distribución superpuesta de linajes. En particular, el estudio del linaje Cohen más dominante anidado con el haplogrupo prevalente J1-P58, junto con aportes históricos expertos, podría otorgar una visión crítica a la comprensión de los hebreos en el Viejo Mundo. Además, los antiguos estudios de ADN del Levante pueden ofrecer información directa. De hecho, un estudio reciente reveló la presencia de cromosomas Y J1a-P58 y J2-M12, frecuentes entre los judíos contemporáneos, en dos muestras de cananeos que datan de 3,700 ybp37.
Es difícil simular el número máximo de hombres que habrían necesitado estar presentes en la fundación de la Judería asquenazí para producir el fuerte patrón de eventos fundacionales observados actualmente en el grupo de cromosomas Y asquenazí. Faltan índices fundamentales como las frecuencias de los respectivos haplogrupos en las poblaciones ancestrales de deme y el grado de introgresión que ha ocurrido a lo largo de las generaciones. Sin embargo, las frecuencias contemporáneas de estos linajes fundadores del cromosoma Y, en vista de sus frecuencias mucho más bajas en poblaciones no judías de Oriente Medio y Europa, sugieren que tales linajes deben haber estado presentes en el inicio de la Judería asquenazí, cuando el escaso número de hombres comprendía la población fundadora, y permitiendo que la deriva desempeñara un papel en el establecimiento de sus altas frecuencias actuales. De hecho, la reconstrucción de la reciente historia judía asquenazí a partir de genomas enteros sugirió un cuello de botella de solo 350 individuales26, 38. Es importante notar que si bien este cuello de botella no coincide necesariamente con el tamaño efectivo de la población masculina fundacional y los eventos para los judíos asquenazíes, sí nos dice que el ancestro levita asquenazí R1a-Y2619 fue probablemente uno de los hombres fundadores sobre los que se aplicó el cuello de botella. En consecuencia, los datos de secuencias completas de ADN mitocondrial apoyan la misma noción de un número limitado de linajes maternos fundadores principales33.
En conjunto, la magnitud de los datos presentados aquí facilita una resolución casi completa de la historia genómica de los Levitas dentro de los judíos asquenazíes. Se puede argumentar fuertemente que los levitas asquenazíes contemporáneos R1a-Y2619 descienden de un solo ancestro levita que llegó a Europa desde el Levante. La expansión de su linaje masculino directo comenzó en un período de tiempo compatible con el patrón de expansión observado para varios padres fundadores adicionales de la judería asquenazí. Por lo tanto, además de proporcionar información sobre una genealogía masculina única, estos hallazgos son un ejemplo importante y altamente resuelto de la historia demográfica más general de los linajes judíos ashkenazi masculinos y femeninos, en los que un número relativamente pequeño de fundadores hacen contribuciones desproporcionadamente grandes a la variación genómica de la población ashkenazi contemporánea. Se puede argumentar además que el hecho de que los hombres levitas no asquenazíes R1a-Y2619 estén bien distribuidos dentro de la filogenia de los Levitas asquenazíes, en lugar de agruparse por separado, da más credibilidad al escenario de que el hombre R1a-Y2619 entró en Europa con los fundadores asquenazíes. La pregunta más enigmática – el momento y la ubicación en que el fundador del Levita ashkenazi R1a-Y2619 obtuvo el estatus de Levita – sigue sin resolverse. Esta pregunta podría estar más allá del alcance de los estudios genéticos que utilizan la variación genética contemporánea, dada la ausencia hasta la fecha de cualquier hombre probado cuyas líneas se ramificaron entre el tiempo del ancestro masculino directo compartido de todos los hombres R1a-M582, ~3,143 años, y el ancestro masculino directo compartido de los hombres levitas ashkenazi R1a-Y2619 ~1,743 años. Los conocimientos históricos o arqueológicos futuros podrían proporcionar los medios para investigar más a fondo este problema.