Kumulativně umožnil rozsáhlý počet zde shromážděných vzorků v kombinaci s dostupností velmi dobře rozlišených otcovských fylogenezí sahajících mnoho generací zpět formulovat určité závěry o zakladatelské linii aškenázských Levitů haploskupiny R1a s vysokou mírou jistoty. Na počátku studie bylo zjištěno, že celkem 71 jedinců hlásících se k aškenázskému levitskému nebo aškenázskému nelevitskému otcovskému původu patří do haploskupiny R1a-M582 a na základě jejich STR profilů reprezentuje široký rozsah variability zjištěné v rámci této haploskupiny. Je pozoruhodné, že všech 71 chromozomů Y bylo možné přeřadit do jediného rozšiřujícího se kladu vnořeného do haploskupiny R1a-M582 a označeného zde jako R1a-Y2619 podle jedné ze šesti variant zjištěných na této úrovni (obr. 2 a doplňkový obrázek S1). Jedinými dříve popsanými vzorky, které sdílely tento klad, bylo pět vzorků aškenázského židovského původu, u nichž nebyl uveden status levita26. V souladu s tím jsou všichni jedinci R1a-Y2619, ať už se sami hlásí k židovství, nebo ne, ať už jsou aškenázští, nebo ne, ať už jsou levité, nebo ne, přímými mužskými potomky po otcovské linii jednoho společného mužského předka, který žil ~1 743 ybp (tabulka 1). Vzhledem k tomu, že současní muži ze všech dosud odebraných větví R1a-Y2619 nesou jedno z mnoha levitských příjmení, lze s velkou pravděpodobností tvrdit, že tento mužský předek se sám označoval za levitu a mohl nést patronymické příjmení Levita. Rozsah zakladatelské události lze odhadnout výpočtem počtu současných jedinců, u nichž se očekává, že ponesou Y chromozom R1a-Y2619. Dva nezávislé soubory vzorků aškenázských Židů, u nichž nebyl znám status levity, podobně odhadly procento otcovské haploskupiny R1a-Y2619, R1a-M17/M198, v aškenázské populaci na 9,6 %20 a 11,5 %27. První z těchto prací uvádí, že haploskupina R1a-M582 představuje 7,9 % celkové aškenázské populace20. Zde ukazujeme, že všechny aškenázské vzorky patřící do haploskupiny R1a-M582 lze překlasifikovat na R1a-Y2619. Na základě velikosti aškenázské populace ~4 000 00013 mužů, z nichž asi 7,9 % je R1a-Y2619, by bylo ~300 000 aškenázských mužů pocházejících v přímé mužské linii z jediného relativně nedávného předka, přičemž mnozí z těchto mužů se sami označují za levity.
Dva faktory nám umožnily dosáhnout přesnějšího koalescenčního stáří a intervalů spolehlivosti, než jaké byly dříve vypočteny pro aškenázskou levitskou linii R1a18,20 . Zaprvé velký počet celých dostupných chromozomů Y (doplňková tabulka S1) a zadruhé dostupná data z rodokmenu Horowitz, která umožnila vnitřní kalibrační bod (obr. 3). Dobře zdokumentované genealogické záznamy Horowitzova rodokmenu ukazují, jak robustní mohou být sekvence celých chromozomů Y při rekonstrukci otcovských předků. Kalibrační bod byl založen na třech Horowitzových vzorcích s dobře zdokumentovanými genealogiemi s koalescencí 402 ybp. Dva další jedinci, kteří se sami hlásí k dynastii Horowitzů, ale nemají jasné genealogie, se těsně spojili s těmito třemi vzorky. Žádný ze vzorků v bezprostředních sesterských klanech se sám neidentifikoval jako součást dynastie Horowitzů. Pozoruhodné je, že molekulární stáří uzlu zahrnujícího všech pět vzorků Horowitzů bylo vypočteno na ~690 ybp, zatímco genealogická doba koalescence naznačuje, že je 546 ybp. To naznačuje, že uzel Horowitz R1a-YP268 (g.23133909G > A) by mohl představovat skutečnou rekonstruovanou sekvenci celého Y jmenovaného zakladatele tohoto rodokmenu (obr. 3b). Je však důležité poznamenat, že odhadované koalescenční stáří je přibližné, protože několik metodických odchylek může vést k poněkud odlišným výsledkům. Například při použití naší metodiky koalescence pouze na data získaná z platformy Illumina by bylo získáno stáří koalescence ~1398 let ybp pro aškenázský levitský rod R1a-Y2619 a ~645 let ybp pro rod R1a-YP268. Tyto rozdíly lze přičíst možnosti použití delších překrývajících se úseků chromozomu Y při porovnávání vzorků sekvenovaných stejnými stroji a parametry. Nepočítali jsme dobu koalescence na základě STR variací, protože jsme již dříve ukázali, že neposkytují žádné další informace pro nedávno koalescenční klad6.
Navržený blízkovýchodní původ aškenázské linie levitů založený na dříve relativně omezeném počtu hlášených vzorků20 lze nyní považovat za pevně potvrzený. Zatímco nejvyšší frekvence haploskupiny R1a se vyskytují ve východní Evropě18,28, naše data odhalila bohatou variabilitu haploskupiny R1a mimo Evropu, která je fylogeneticky oddělená od typicky evropských větví R1a. Je zřejmé, že R1a-Y2619 je dobře vnořena do plejády fylogeneticky blízkých sesterských skupin z Blízkého východu, jejichž vzorky byly odebrány u íránských Azerů, Kermánců, Jezídů a jednoho muže z Pyrenejského poloostrova. To poskytuje potřebné důkazy o jejím původu (obr. 2 a doplňkový obrázek S1). Přesná migrační cesta R1a-Y2619 do Evropy však zůstává neobjasněna. Většina historických záznamů naznačuje dvě hlavní cesty židovské migrace do Evropy (obr. 1)29,30. Aškenázští Židé jsou považováni za výsledek židovské migrace přes Itálii do údolí Rýna a poté do Polska. Sefardští (španělští) Židé jsou považováni za migranty, kteří migrovali spolu s postupnou islámskou expanzí do severní Afriky a poté do Španělska14. Vzhledem k tomu, že R1a-Y2619 je blízkovýchodního původu, je možné, že k jeho zavlečení do Evropy došlo některou z těchto cest nebo oběma. Výrazný výskyt R1a-Y2619 mezi aškenázskými Židy spolu s přítomností všech známých větví R1a-Y2619 u aškenázských Židů samozřejmě svádí k závěru, že jeho migrační cesta z Levanty přímo souvisela se zakladateli Aškenázů. Tyto skutečnosti však mohou odrážet spíše jen expanzi v rámci aškenázských Židů než důkaz prvního příchodu R1a-Y2619 s aškenázskou populací nebo do ní. Ještě více matoucí je skutečnost, že genealogické záznamy rabínské dynastie Horowitzů, kteří jsou nyní prokazatelně nositeli chromozomu Y R1a-Y2619, naznačují jejich přítomnost na Pyrenejském poloostrově v 15. století a pravděpodobně i dříve (obr. 3)17 . Ve skutečnosti jsou opakované migrace Židů, kteří mohli být nositeli chromozomu Y R1a-Y2619, do Katalánska doloženy od 4. století a během muslimské expanze na Pyrenejský poloostrov31. Navíc vzhledem k tomu, že Katalánsko bylo od roku 800 n. l. opět křesťanským územím, mohli proto-horowitzští předci R1a-Y2619 představovat také migraci aškenázských Židů na Pyrenejský poloostrov. Přítomnost R1a-Y2619 ve Španělsku v 15. století by tedy nemohla být důkazem pro první příchod linie R1a-Y2619 na Pyrenejský poloostrov, protože by mohla jednoduše odrážet opakované a nezaznamenané přesuny Židů tam a zpět mezi východní a západní Evropou a Pyrenejským poloostrovem. Již předchozí důkazy z mtDNA a autozomálních markerů naznačily pravděpodobný tok genů mezi aškenázskými a sefardskými Židy v rámci Evropy32,33,34. Celkově se naše výsledky přiklánějí k jediné cestě vstupu do Evropy v rámci aškenázské migrace a expanze v Evropě. Protože doba koalescence všech současných levitů R1a-Y2619 je ~1 743 ybp, tedy v době římské exilové diaspory, a každá z větví R1a-Y2619 se vyskytuje u aškenázských Židů, naše výsledky neodpovídají scénáři rychlé expanze v Levantě a následného šíření více cestami do Evropy. Výsledky z ne-škenázských vzorků R1a-Y2619 levitů také naznačují jedinou cestu expanze. Pokud by linie R1a-Y2619 levitů byla v Levantě již stejně rozšířená, očekávalo by se, že se ne-aschenázští levité oddělí od aškenázských vzorků levitů před dobou diaspory a vytvoří jiný shluk. Výsledky však ukazují, že neaškenázští levité spadají do několika relativně nedávno koalescentních podvětví R1a-Y2619 (obr. 2 a doplňkový obrázek S1). Tento vzorec je spíše slučitelný s kontinuálním tokem genů z aškenázské populace do ne-aschenázské populace během diaspory než s více cestami vstupu haploskupiny R1a-Y2619 do Evropy z Levanty.
Další podpora pro jedinou migrační cestu a následnou expanzi v rámci aškenázských Židů, vychází z pozorovaných vzorců expanze dalších aškenázských haploskupin. V souladu s tím je důležité rozlišovat mezi specifickou aškenázskou zakladatelskou událostí pro haploskupinu R1a-Y2619, která mohla být výsledkem příznivého socioekonomického nebo jiného postavení, a obecnou expanzí aškenázských zakladatelských linií chromozomu Y, včetně aškenázských levitů R1a-Y2619. Za tímto účelem jsme studovali několik dalších haploskupin, které byly dříve popsány jako převládající mezi aškenázskými Židy19,24,35 . Naše výsledky ukazují, že rozšíření R1a-Y2619 mezi aškenázskými Židy není specifické nebo jedinečné pro tuto haploskupinu. Například koalescence haploskupin G-M377 (Doplňkový obrázek S4) a Q-M242 (Doplňkový obrázek S5), o nichž je známo, že každá z nich představuje 5 % aškenázské otcovské variace24 , koalescence nastává v čase ~1223 ybp, respektive ~1672 ybp (tabulka 1). Zde přijatá strategie chromozomu Y nám umožnila rozdělit linie haploskupiny E na její drobné podskupiny (doplňkový obrázek S3). Stáří koalescence haploskupin E-Z838, E- PF3780, E- B923 a E-B933, o nichž je známo, že kumulativně představují 20 % aškenázské otcovské variability24 , bylo odhadnuto na ~1 200-1 600 ybp. Tento vzorec vícenásobných zakladatelských událostí, který nebyl pozorován u španělských Židů36 , poskytuje další podporu pro to, že aškenázský předek R1a-Y2619 Levita se do Evropy dostal aškenázskou cestou, nikoliv prostřednictvím židovské expanze na Pyrenejský poloostrov. U aškenázských židovských otcovských předků jsou jasně patrné i další vzorce. Haploskupina T-M70, rozšířená na Blízkém východě, je přítomna i mezi aškenázskými Židy (doplňkový obrázek S6). Genotypované vzorky vykazovaly hluboce zakořeněné štěpení, které pravděpodobně poukazuje na zachování starobylé diverzity této haploskupiny v Levantě, sahající až do pleistocénu. Mezitím vzor pozorovaný u haploskupiny R1b-M269 (Doplňkový obrázek S7), převažující v západní Evropě, naznačuje především opakovanou bodovou introgresi evropských chromozomů Y do aškenázské komunity a je v souladu s předchozími zprávami36.
Dále jsme porovnávali nejčastější zakladatelskou linii nalezenou mezi aškenázskými muži Cohenů, vnořenou do haploskupiny J1a-P58, s aškenázskou levitskou linií R1a-Y2619. Je zřejmé, že příslušníci levitské kasty R1a-Y2619 a cohenovské kasty J1a-P58 nemají společného mužského předka v časovém rámci biblického vyprávění. Stejně jako v případě aškenázských levitů tvořila aškenázská linie Cohenů J1a-P58 těsný shluk vnořený do blízkovýchodní sady vzorků potvrzující její původ a v naší studii byla společná pouze s neaškenázskými Coheny (doplňkový obrázek S2). Na rozdíl od vzorce získaného pro aškenázské levity R1a-Y2619 však shluk koaleskuje ~2570 ybp, což ukazuje na počátek jeho expanze v období před diasporou. Rozsáhlé úplné sekvence chromozomu Y z většího počtu vzorků Cohenů z více židovských komunit na pozadí jiných levantských populací, včetně zjištění rodinných a klanových specifických variant, by byly velmi informativní při řešení zjištění, že aškenázští a neaškenázští jedinci Cohenů sdílejí překrývající se distribuci linií. Zejména studium nejdominantnější cohenovské linie vnořené do převládající haploskupiny J1-P58 by spolu s odbornými historickými podklady mohlo poskytnout zásadní vhled do chápání Hebrejců ve Starém světě. Kromě toho mohou přímé informace nabídnout studie starověké DNA v Levantě. Nedávná studie totiž odhalila přítomnost chromozomů Y J1a-P58 i J2-M12, častých mezi současnými Židy, ve dvou kanaánských vzorcích datovaných do období 3700 let před naším letopočtem37.
Je obtížné simulovat maximální počet mužů, který by musel být přítomen při založení aškenázských Židů, aby vznikl silný vzorec zakladatelských událostí, který je v současné době pozorován v aškenázském fondu chromozomů Y. V případě, že by se jednalo o chromozomy J1a-P58 a J2-M12, bylo by to obtížné. Chybí základní ukazatele, jako jsou frekvence příslušných haploskupin v populacích předků dem a rozsah introgrese, ke které došlo v průběhu generací. Současné frekvence těchto zakladatelských linií chromozomu Y však vzhledem k jejich mnohem nižším frekvencím v nežidovských populacích Blízkého východu a Evropy naznačují, že tyto linie musely být přítomny již při vzniku aškenázských Židů, kdy zakladatelskou populaci tvořil mizivý počet mužů, a umožňují, aby při vzniku jejich současných vysokých frekvencí hrál roli drift. Rekonstrukce nedávné historie aškenázských Židů z celých genomů skutečně naznačila úzké hrdlo pouhých 350 jedinců26,38. Je důležité poznamenat, že ačkoli se toto úzké hrdlo nemusí nutně shodovat se zakládající efektivní velikostí mužské populace a událostmi u aškenázských Židů, říká nám, že aškenázský levitský předek R1a-Y2619 byl pravděpodobně mezi zakládajícími muži, na které se úzké hrdlo vztahovalo. Stejně tak údaje z kompletních sekvencí mitochondriální DNA podporují stejnou představu o omezeném počtu hlavních zakládajících mateřských linií33.
Souhrnně vzato, rozsah zde prezentovaných údajů usnadňuje téměř úplné vyřešení genomické pohádky o levitech v rámci aškenázských Židů. Lze důrazně tvrdit, že současní aškenázští levité R1a-Y2619 pocházejí z jediného levitského předka, který přišel do Evropy z Levanty. Expanze jeho přímé mužské linie začala v časovém rámci kompatibilním s modelem expanze pozorovaným u několika dalších zakladatelů aškenázských Židů. Kromě toho, že tato zjištění umožňují nahlédnout do rodokmenu jednoho muže, jsou tedy důležitým a velmi vyřešeným příkladem obecnější demografické historie aškenázských mužských a ženských židovských linií, v níž relativně malý počet zakladatelů neúměrně výrazně přispívá k současné genomické variabilitě aškenázské populace. Lze dále tvrdit, že skutečnost, že neaškenázští levité R1a-Y2619 jsou dobře rozloženi v rámci aškenázské levitské fylogeneze, a nikoliv seskupeni odděleně, dodává větší důvěryhodnost scénáři, že muži R1a-Y2619 se do Evropy dostali spolu s aškenázskými zakladateli. Nejzáhadnější otázka – kdy a kde získal zakladatel aškenázského levitského rodokmenu R1a-Y2619 status levity – zůstává nevyřešena. Tato otázka může být mimo rámec genetických studií využívajících současnou genetickou variabilitu, vzhledem k tomu, že doposud nebyli testováni žádní muži, jejichž linie se rozvětvily mezi dobou společného přímého mužského předka všech mužů R1a-M582, ~3 143 ybp, a společným přímým mužským předkem aškenázských levitů R1a-Y2619 ~1 743 ybp. Budoucí historické nebo archeologické poznatky by mohly poskytnout prostředky k dalšímu zkoumání této otázky.
.