Document Type : Original Article
1 Department of Biology, College of Education, University of Garmian, Kalar, Kurdistan Region, Iraq
2 College of Medicine, University of Garmian, Kalar, Kurdistan Region, Iraq
Toll-like receptors (TLRs), encoded by innate immune genes, functions to detect microbial ligands. Studies found that polymorphisms in TLR genes among populations are associated with diseases. Among the classes of TLRs, TLR-4 is important for exploring bacterial lipopolysaccharide. TLR-4 Thre399Ile is one of the most important non-synonymous variant which varies among different background populations. Kurds, an ethnic population descended mainly from indigenous inhabitants of Zagros mountain survived for millennia, with nearly 40 million populations, are distributed in Iraq, Turkey, Iran and Syria. There are no studies on genetic variations of TLR4 in this population, particularly in Iraqi Kurdistan. The aim of this study is to find the percentage frequency of TLR4 Thr399Ile variant in Kurdish populations in Iraq, particularly in Sulaymaniyah province. The percentage of single nucleotide polymorphisms (SNPs) of the TLR4 gene was determined by polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP). The results showed that the heterozygous variant of TLR4 Thr399Ile is 7.1% (n=85) in Sulaymaniyah populations. However, no homozygous mutant variant was found; this suggests that it is either absent or seldom among Kurds in the region. This study emphasises the perceptiveness of the TLR SNP in Sulaymaniyah populations that recommends future study for linking this genetic variation with both infectious and immunological diseases, in addition to historical, anthropological and archaeological studies.
Innate immune genes, for example those encode receptors that sense microbes, are evolutionarily ancient genes which are conserved and varied possibly due to genetic drift or selective pressure for fighting infections such as bacteria (like tuberculosis), fungi, viruses, or parasites like malaria. These pathogen recognition receptors, known as toll like receptors (TLRs), are a class of 10 evolutionarily transmembrane receptors, named TLR1-10. Among the TLRs, TLR4 plays an important role in innate immune response against Gram negative bacteria. The in vitro functions and structures of TLR4 has been previously studied that is found on the innate immune cells recognising bacterial lipopolysaccharide (LPS) which can trigger inflammatory cytokines through NF-kB pathways . Thus, studying this gene is needed for finding correlations between healthy and diseased populations.
Human genome project found 2.1 million single nucleotide polymorphisms (SNPs) with 1% of them related to protein functions and revealed that genetic variations are often common among the genes responsible for immunity . Non- synonymous SNPs have been identified in TLR4 gene that modify amino acid sequences in the gene’s leucine-rich repeat domain which recognizes microbial ligands. Among these, a SNP (rs4986791), in which an amino acid threonine (Thr399) is substituted by Isoleucine (Ile399), is known as TLR4 Thr399Ile variant (1196C/T). Research has been conducted on TLR SNPs in most ethnic background populations in the world. Nonetheless, investigation of TLRs is lacking in Iraqi Kurdish populations.
Despite of associations of TLR4 Thr399Ile variant with bacterial infections, this variant has also been linked with prognostic chemotherapy and cancer genetic marker , cancer genetic risk factor , and a damaging impact on functions and structures of the protein . Additionally, this SNP seems to be an important variant for studying ethnic populations since its frequency varies in different populations in different continents including Asia, Africa and Europe. For instance, it is absent or very rare in Chinese  or Korean populations , and infrequent in Africans, however, it is significantly prevalent among Caucasian, European and Middle eastern populations .
Kurds, the most ancient indigenous people, lived mainly in Zagros mountains and often northern Mesopotamia for thousands years . Approximately, 40 million Kurdish populations are currently distributed in Turkey, Iran, Iraq and Syria, with 100 thousand Kurds have also been refuged all over the world, including European countries. Kurdish people, from prehistoric times until Ottoman period, appear to be homogeneous people as they have been endogamous and tribal nomads inside Kurdistan lands . Kurdistan has also an ancient history of first village around 12,000 years ago . However, research on this area was previously less focused due to wars and internal conflicts. Within last decade, the region has become the most prominent location for archaeological discoveries in the middle east . Kurdistan, as a middle eastern region between Arabia and Eurasia, might be considered as an important site for ancient human migrations because the migrations outside of Africa spread via Arabic Peninsula to Eurasia  and possibly through Kurdistan where ancient caves and archaic villages have been found [8,11]. Therefore, the aim of this study is to explore the percentage frequency of TLR4 Thr399Ile variant in Kurdish populations of Iraqi Kurdistan, particularly in Sulaymaniyah province. The hypothesis is whether this polymorphism is either common or rare or absent in Kurdish population.
2. Materials and methods
2.1. Study area and population
EDTA preserved venous blood were randomly collected from 85 staffs and students in the University of Garmian. The collected population samples represent populations of Sulaymaniyah province, Kurdistan Regional Government (KRG) of Iraq. The sample distributions among Sulaymaniyah province are shown in the Supplementary Table 1. Only Kurdish people, whose family originally descended from Kurdish ancestry, were included in this study. Others with different backgrounds, such as Arabic, were excluded. An informed consent form was filled by each participant and the research was approved by an ethical committee of Department of Biology, College of Education, University of Garmian (Ethical committee No.2: 07/01/2020). The study was conducted in accordance with the Declaration of Helsinki. The collected blood samples were kept in a refrigerator until total genomic DNA were extracted. This study was conducted in April 2019 and its geographical area is shown in Fig. 1.
Figure 1: Map of the study area (Sulaymaniyah Province). The map was designed by Landsatlook viewer (USGS Products, Data available from the U.S. Geological Survey). The black line indicates the geographical area of Sulaymaniyah province
2.2. Genomic DNA extraction:
Total DNA extraction kit (GENET BIO CO., Daejeon, KR) was used by mixing each blood sample (200 µl) with 20 µl proteinase K in a 1.5 ml tube and incubated at 56°C for 10 minutes. Following manufacturer’s instructions, purified genomic DNA samples were eluted with 200 µl elution buffer. The extracted genomic DNA samples were kept at -30°C.
2.3. Genotyping of TLR4 Thre399Ile SNP:
Polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) were used to amplify the single nucleotide polymorphism of TLR4 Thr399Ile variant gene and then digested with Hinf1 restriction enzymes.
Forward primer, TLR4 Thr399Ile F: GGTTGCTGTTCTCAAAGTGATTTTGGGAGAA; and reverse primer, TLR4 Thr399Ile R: ACCTGAAGACTGGAGAGTGAGTTAAATGCT, (Macrogen Co., Seoul, KR) were used for the SNP rs4986791 . The details of the TLR4 Thre399Ile variant gene is described in Fig.2.
Figure. 2. Homo sapiens TLR4 mRNA for toll-like receptor 4, complete cds: GenBank: AB445638.1. When the primer is designed, the highlighted c (red) is altered to g (red) in the F primer to fit with the restriction site Hinf1 (GANTC)= gaatc. N= any nucleotide. This SNP is known as rs4986791 or 1196 C/T which means the nucleotide cytosine 1196 (C1196) is mutated to threonine 1996 (T1196). The wild type codon (aCc) is Threonine (Thre) and its product size= 406 bp. The SNP codon (aTc) is Isoleucine (Ile) and its product size= 406-29= 377 bp. Thus, the amino acid Thr399 is changed to Ile399. Therefore, the SNP is also called TLR4 Thre399Ieu
PCR reactions were prepared by mixing 10 µl Add start mastermix 2 x (Add Bio Inc, Korea) with 1 µl of each primer (5 pM) and 10 µl genomic DNA. The reactions were run in a thermal cycler (Master cycler nexus, Eppendorf AG, Hamburg, Germany) by initial denaturation at 95°C for 5 minutes; followed by 40 cycles of 95°C for 30 seconds, annealing temperature at 62 for 45 seconds and 72°C for 30 seconds; and finally extended at 72°C for 5 minutes. Ten (10) µl PCR products were digested by 1 µl of Hinf1 restriction endonucleases (New England Biolabs, Ipswich, MA, USA) in 4 µl Hinf1 buffer 4x and 10 µl nuclease free water (Promega) incubated at 37°C for 2 hours. The digested products were mixed with gel loading dye, Purple (6X) (New England Biolabs, Ipswich, MA, USA) and explored on 2% agarose gel, stained with 10 µl Prime safe dye (GENET BIO CO., Daejeon, KR), run with 1x TBE buffer at 100 v, 80 mA for 2 hours. The wild type product size (406 bp) is cleaved with Hinf1 restriction enzyme to form 377 bp and 29 bp in case of the mutant variant.
To test the Hardy-Weinberg equilibrium, Gene Calculators online program (Chi-square, P>0.05) was used (https://www.genecalculators.net/pq-chwe-check.html).
3. Results and discussion
Current study investigated TLR4 Thr399Ile variant in Kurdish populations in Sulaymaniyah province using PCR-RFLP. The results of this study showed that 6 out of 85 (7.1 %) of the samples have carried the heterozygous mutant alleles of the TLR4 Thr399Ile variant. An example was shown in Supplementary Figure 1. No homozygous variant was identified and the genotype frequency was not deviated from Hardy-Weinberg equilibrium (Chi-square P=0.735). This variant is infrequent among other populations in middle eastern countries. For example, 5.8 % of TLR4 Thr399Ile heterozygous mutant was determined in healthy Iranian populations and homozygous mutants were only found in chronic cutaneous leishmaniasis cases . Likewise, no homozygous mutants of TLR4 Thr399Ile were found in Turkish populations ,but the heterozygous variants were 9% . A study also suggested weak selective pressure of the TLR4 Thr399Ile variant in fighting infections since the alleles are heterogeneous among Iranian ethnic populations including Iranian Kurds . It has been shown that TLR polymorphisms are associated with both reduced responses to bacterial infections and susceptibility to Gram negative bacterial infections due to LPS signaling impairment . Several studies have focused on laboratory mice, or different ethnic populations including Caucasians and Koreans , Iranians [12,14], sub-Saharan Africa, Europe and East Asia . For instance, the TLR4 Thr399Ile SNPs (rs4986791) is only 0.79% in Africans, 8.51% in Europeans, and 0% in East Asians and this SNP is considered as ‘possibly damaging’ to the function and structure of the protein . Interestingly, in Caucasian populations, TLR4 Thr399Ile is 17.7% heterozygous variant with no any homozygous mutants were found . Up to my best knowledge by the time of writing this manuscript (08/11/2019), no articles were found in Pubmed- NCBI data base when searching for key words as ‘TLR4 SNP Iraq’. Although, when searched for ‘TLR4 polymorphism Iraq’ only one article was published that was about TLR4 Asp299Gly SNP (rs4986790) and the samples were taken from Jordanian populations but one of its authors is from Iraq . Similarly, no data were found for TLR4 SNP in Syria as a country of having Kurdish populations.
In spite of its limitations, this is the first time, TLR4 gene polymorphisms are investigated in Kurdistan Region of Iraq. However, when searching in google for ‘TLR SNPs Iraq, in other provinces of Iraq, where the majority of populations are almost Arabs, several articles have been published in several journals which were not cited in this manuscript due to their poor quality. The limitation of the current study is that the other SNP, TLR4 Asp299Gly variant was not investigated. It is worth to mention that this SNP is often co-segregated with the SNP investigated in this study. The co-segregated TLR4 Asp299Gly/Thr399Ile mutants play an important role in inflammatory response to bacterial LPS challenge similar to the wildtype variants known as ‘functional neutrality’ . Both segregated mutations may have been originated from Africa approximately 60,000 years ago, spreading out through Arabian Peninsula and middle eastern route to Europe as a result of genetic drift . Further researches are required to investigate both SNPs in Kurdish ethnic populations in all Kurdish inhabiting provinces and countries including Iraq, Iran, Turkey and Syria. Particularly, comparing the mutants among the main Kurdish tribes will give an idea about the selective pressure and homogeneity of the TLR4 SNP alleles in Kurdish ethnic populations.
Genetic variations, together with historical, demographic, linguistic, anthropological and archaeological findings, can be exploited as a population marker to study human evolution . Discovery of both Neanderthal and Denisovan like haplotypes of TLR1, 6, and 10 in modern human indicates gene flows among archaic human species . Investigations of gene polymorphisms in Kurdish people might be of archaeological and anthropological interests since the area was a homeland of ancient civilizations and archaic human species. For instance, in Shanidar, (an ancient cave in Zagros mountain belonging to Iraqi Kurdistan), Neanderthal bones were found and dated back to about 40,000-70, 000 years ago . This is apparently the period when both homo sapiens and Homo Neanderthals encountered and interbred [19,20] Research about ancient DNA has not been recorded yet in human species in the Kurdistan caves. It might be fascinating for both archaeologists and evolutionary geneticists to seek polymorphisms in innate immune genes between the two human species.
In conclusions, Iraqi Kurdish populations, lived in Sulaymaniyah province, bear only the heterozygous variant of TLR4 Thr399Ile (7.1%). This suggested the rarity of the homozygous mutant in the province. Further study with large size samples is recommended to confirm this, taking samples from other provinces in Kurdistan region. This also suggests that TLR4 Thr399Ile SNP in Kurdish people is neither very common as in Caucasian nor rare as in East Asian nor uncommon as in African people. This result also suggests that the Kurds has TLR4 Thr399Ile SNP frequency similar to that of Europeans or Middle Eastern nations as in Iran; it supports both Indo-Iranian and Endo-European language families. It seems to be exciting to compare the TLR4 SNPs between modern and ancient human species found in the caves of Kurdistan. Finally, this investigation unlocks a main gate for further research on genetic associations of this innate immune gene with infectious, immunological or non-infectious diseases, and it is also anthropologically and genetically important for studying Kurdish population.
Dr. Sirwan M. Amin, Mr. Mariwan AbdulRahman, Ms. Sakar Nariman, Ms. Soma Mahmood, Mr, Chya Mustafa Othman, and Ms. Dekan Kamaran helped me during sample collections. I acknowledge all staffs and students who donated the blood samples.
Supplementary information related to this article can be found at: (Supplementary ) Niranji, Sh. S. Passer 1 (2020) 32- 36
- S. Akira, S. Uematsu, O. Takeuchi, Pathogen recognition and innate immunity, Cell. 124 (2006) 783–801. https://doi.org/10.1016/j.cell.2006.02.015.
- J.C. Venter, M.D. Adams, E.W. Myers, P.W. Li, R.J. Mural, G.G. Sutton, H.O. Smith, M. Yandell, C.A. Evans, R.A. Holt, J.D. Gocayne, P. Amanatides, R.M. Ballew, D.H. Huson, J.R. Wortman, Q. Zhang, C.D. Kodira, X.H. Zheng, L. Chen, M. Skupski, G. Subramanian, P.D. Thomas, J. Zhang, G.L. Gabor Miklos, C. Nelson, S. Broder, A.G. Clark, J. Nadeau, V.A. McKusick, N. Zinder, A.J. Levine, R.J. Roberts, M. Simon, C. Slayman, M. Hunkapiller, R. Bolanos, A. Delcher, I. Dew, D. Fasulo, M. Flanigan, L. Florea, A. Halpern, S. Hannenhalli, S. Kravitz, S. Levy, C. Mobarry, K. Reinert, K. Remington, J. Abu-Threideh, E. Beasley, K. Biddick, V. Bonazzi, R. Brandon, M. Cargill, I. Chandramouliswaran, R. Charlab, K. Chaturvedi, Z. Deng, V. Di Francesco, P. Dunn, K. Eilbeck, C. Evangelista, A.E. Gabrielian, W. Gan, W. Ge, F. Gong, Z. Gu, P. Guan, T.J. Heiman, M.E. Higgins, R.-R. Ji, Z. Ke, K.A. Ketchum, Z. Lai, Y. Lei, Z. Li, J. Li, Y. Liang, X. Lin, F. Lu, G. V Merkulov, N. Milshina, H.M. Moore, A.K. Naik, V.A. Narayan, B. Neelam, D. Nusskern, D.B. Rusch, S. Salzberg, W. Shao, B. Shue, J. Sun, Z.Y. Wang, A. Wang, X. Wang, J. Wang, M.-H. Wei, R. Wides, C. Xiao, C. Yan, A. Yao, J. Ye, M. Zhan, W. Zhang, H. Zhang, Q. Zhao, L. Zheng, F. Zhong, W. Zhong, S.C. Zhu, S. Zhao, D. Gilbert, S. Baumhueter, G. Spier, C. Carter, A. Cravchik, T. Woodage, F. Ali, H. An, A. Awe, D. Baldwin, H. Baden, M. Barnstead, I. Barrow, K. Beeson, D. Busam, A. Carver, A. Center, M.L. Cheng, L. Curry, S. Danaher, L. Davenport, R. Desilets, S. Dietz, K. Dodson, L. Doup, S. Ferriera, N. Garg, A. Gluecksmann, B. Hart, J. Haynes, C. Haynes, C. Heiner, S. Hladun, D. Hostin, J. Houck, T. Howland, C. Ibegwam, J. Johnson, F. Kalush, L. Kline, S. Koduru, A. Love, F. Mann, D. May, S. McCawley, T. McIntosh, I. McMullen, M. Moy, L. Moy, B. Murphy, K. Nelson, C. Pfannkoch, E. Pratts, V. Puri, H. Qureshi, M. Reardon, R. Rodriguez, Y.-H. Rogers, D. Romblad, B. Ruhfel, R. Scott, C. Sitter, M. Smallwood, E. Stewart, R. Strong, E. Suh, R. Thomas, N.N. Tint, S. Tse, C. Vech, G. Wang, J. Wetter, S. Williams, M. Williams, S. Windsor, E. Winn-Deen, K. Wolfe, J. Zaveri, K. Zaveri, J.F. Abril, R. Guigó, M.J. Campbell, K. V Sjolander, B. Karlak, A. Kejariwal, H. Mi, B. Lazareva, T. Hatton, A. Narechania, K. Diemer, A. Muruganujan, N. Guo, S. Sato, V. Bafna, S. Istrail, R. Lippert, R. Schwartz, B. Walenz, S. Yooseph, D. Allen, A. Basu, J. Baxendale, L. Blick, M. Caminha, J. Carnes-Stine, P. Caulk, Y.-H. Chiang, M. Coyne, C. Dahlke, A.D. Mays, M. Dombroski, M. Donnelly, D. Ely, S. Esparham, C. Fosler, H. Gire, S. Glanowski, K. Glasser, A. Glodek, M. Gorokhov, K. Graham, B. Gropman, M. Harris, J. Heil, S. Henderson, J. Hoover, D. Jennings, C. Jordan, J. Jordan, J. Kasha, L. Kagan, C. Kraft, A. Levitsky, M. Lewis, X. Liu, J. Lopez, D. Ma, W. Majoros, J. McDaniel, S. Murphy, M. Newman, T. Nguyen, N. Nguyen, M. Nodell, S. Pan, J. Peck, M. Peterson, W. Rowe, R. Sanders, J. Scott, M. Simpson, T. Smith, A. Sprague, T. Stockwell, R. Turner, E. Venter, M. Wang, M. Wen, D. Wu, M. Wu, A. Xia, A. Zandieh, X. Zhu, The Sequence of the Human Genome, Science (80-. ). 291 (2001) 1304–1351. https://doi.org/10.1126/science.1058040.
- C. Bergmann, H.S. Bachmann, A. Bankfalvi, R. Lotfi, C. Pütter, C.A. Wild, P.J. Schuler, J. Greve, T.K. Hoffmann, S. Lang, A. Scherag, G.F. Lehnerdt, Toll-like receptor 4 single-nucleotide polymorphisms Asp299Gly and Thr399Ile in head and neck squamous cell carcinomas, J. Transl. Med. 9 (2011) 1–9. https://doi.org/10.1186/1479-5876-9-139.
- L. Ding, Q. Jiang, G. Li, J. Shen, J. Du, X. Lu, X. Xiong, Comprehensive assessment of association between TLR4 gene polymorphisms and cancer risk: a systematic meta-analysis, Oncotarget. 8 (2017) 100593–100602. https://doi.org/doi: 10.18632/oncotarget.21543.
- L.B. Barreiro, M. Ben-Ali, H. Quach, G. Laval, E. Patin, J.K. Pickrell, C. Bouchier, M. Tichit, O. Neyrolles, B. Gicquel, J.R. Kidd, K.K. Kidd, A. Alcaïs, J. Ragimbeau, S. Pellegrini, L. Abel, J.L. Casanova, L. Quintana-Murci, Evolutionary dynamics of human toll-like receptors and their different contributions to host defense, PLoS Genet. 5 (2009). https://doi.org/10.1371/journal.pgen.1000562.
- Y.S. Kim, Y.J. Hwang, S.Y. Kim, S.M. Yang, K.Y. Lee, I.B. Park, Rarity of TLR4 Asp299Gly and Thr399Ile polymorphisms in the Korean population, Yonsei Med. J. 49 (2008) 58–62. https://doi.org/10.3349/ymj.2008.49.1.58.
- M. Izady, THE KURDS: A Concise Handbook, Taylor & Francis, New York, 1992.
- R. Matthews, W. Matthews, A. Richardson, K.R. Raheem, S. Walsh, K.R. Aziz, R. Bendrey, J. Whitlam, C. M., A. Bogaard, I. Iversen, D. Mudd, S. Elliott, The early Neolithic of Iraqi Kurdistan: current research at Bestansur, Shahrizor Plain, CNRS, 2019.
- K. Kopanias, J. MacGinnis, J. Ur, Archaeological Projects in the Kurdistan Region in Iraq, (2015). https://dash.harvard.edu/handle/1/14022526.
- S.J. Armitage, S.A. Jasim, A.E. Marks, A.G. Parker, V.I. Usik, H.P. Uerpmann, The southern route “out of Africa”: Evidence for an early expansion of modern humans into Arabia, Science (80-. ). 331 (2011) 453–456. https://doi.org/10.1126/science.1199113.
- T.D. Stewart, The Neanderthal Skeletal Remains from Shanidar Cave, Iraq: A Summary of Findings to Date, Proc. Am. Philos. Soc. 121 (1977) 121–165. http://www.jstor.org/stable/986524.
- S. Ajdary, M.M. Ghamilouie, M.H. Alimohammadian, F. Riazi-Rad, S.R. Pakzad, Toll-like receptor 4 polymorphisms predispose to cutaneous leishmaniasis, Microbes Infect. 13 (2011) 226–231. https://doi.org/10.1016/j.micinf.2010.10.018.
- M. Guven, Z. Ismailoglu, B. Batar, S. Unal, I. Onaran, B. Karadag, Z. Ongen, The effect of genetic polymorphisms of TLR2 and TLR4 in Turkish patients with coronary artery disease, Gene. 558 (2015) 99–102. https://doi.org/10.1016/j.gene.2014.12.047.
- M. Ioana, B. Ferwerda, S. Farjadian, L. Ioana, A. Ghaderi, M. Oosting, L.A.B. Joosten, J.W.M. Van Der Meer, G. Romeo, D. Luiselli, D. Dediu, M.G. Netea, High variability of TLR4 gene in different ethnic groups in Iran, Innate Immun. 18 (2012) 492–502. https://doi.org/10.1177/1753425911423043.
- J. Karasneh, M. Bani-Hani, A. Alkhateeb, A. Hassan, F. Alzoubi, M. Thornhill, TLR2, TLR4 and CD86 gene polymorphisms in recurrent aphthous stomatitis, J. Oral Pathol. Med. 44 (2015) 857–863. https://doi.org/10.1111/jop.12298.
- B. Ferwerda, M.B.B. Mccall, S. Alonso, E.J. Giamarellos-bourboulis, M. Mouktaroudi, N. Izagirre, D. Syafruddin, G. Kibiki, T. Cristea, A. Hijmans, L. Hamann, S. Israel, Evolutionary Pressure During Migration of Modern Humans, Pnas. 104 (2007) 16645–16650.
- L.L. Cavalli-Sforza, M.W. Feldman, The application of molecular genetic approaches to the study of human evolution, Nat. Genet. 33 (2003) 266–275. https://doi.org/10.1038/ng1113.
- M. Dannemann, A.M. Andrés, J. Kelso, Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors, Am. J. Hum. Genet. 98 (2016) 22–33. https://doi.org/10.1016/j.ajhg.2015.11.015.
- S. Sankararaman, N. Patterson, H. Li, S. Pääbo, D. Reich, The Date of Interbreeding between Neandertals and Modern Humans, PLoS Genet. 8 (2012). https://doi.org/10.1371/journal.pgen.1002947.
- M. Kuhlwilm, I. Gronau, M.J. Hubisz, C. de Filippo, J. Prado-Martinez, M. Kircher, Q. Fu, H.A. Burbano, C. Lalueza-Fox, M. de la Rasilla, A. Rosas, P. Rudan, D. Brajkovic, Ž. Kucan, I. Gušic, T. Marques-Bonet, A.M. Andrés, B. Viola, S. Pääbo, M. Meyer, A. Siepel, S. Castellano, Ancient gene flow from early modern humans into Eastern Neanderthals, Nature. 530 (2016) 429. https://doi.org/10.1038/nature16544.