The current range of the Asian elephant is fragmented and restricted to southern Asia. Its historical range was far wider and extended from Anatolia and the Levant to Central China. The fossil record from these peripheral populations is scant and we know little of their relationship to modern Asian elephants. To gain a first insight to the genetic affinity of an E. maximus population that once inhabited Turkey we sequenced ca. 570 bp mtDNA from four individuals dating to ~3500 cal. BP. We show that these elephants carried a rare haplotype previously only observed in one modern elephant from Thailand. These results clarify the taxonomic identity of specimens with indeterminate morphologies and show that this ancient population groups within extant genetic variation. By placing the age of the common ancestor of this haplotype in the interval 3.7–58.7 kya (mean = 23.5 kya) we show that range-wide connectivity occurred at some time or times since the start of MIS 3, ~57 kya, probably reflecting range and population expansion during a favourable climatic episode. The genetic data do not distinguish natural versus anthropogenic origin of the Near Eastern Bronze Age population, but together with archaeological and paleoclimatic data they allow the possibility of a natural westward expansion around that time.
The Asian elephant (Elephas maximus, L. 1758) is an endangered species whose current distribution extends from the Indian subcontinent to Borneo (Figure 1A). Its range was far wider during the Pleistocene and early Holocene, extending from Anatolia and the Levant in the west, along a narrow land corridor by the southern Asian coastline, to the Indian subcontinent and as far eastward as southern and central China (Santiapillai and Jackson 1990, Sukumar 2011) (Figure 1A). Since the fossil record is scant, often comprising single fragmented teeth or long bones that can be confused with other elephantid species, their affinities and distributional change are poorly understood (Turvey et al. 2013). Another possibly complicating factor is human translocation of living elephants from their core distribution (e.g. India) as part of emerging trade networks and long-range war campaigns (Collon 1977, Albenda 2008).
|Sample ID||NHM ID||Species (morphology)||Element||DNA|
|46GGB05||EL001||cf. Elephas maximus||m3||No|
|46GGB20||EL005||cf. Elephas maximus||M2||No|
|46GGB01||EL006||cf. Elephas maximus||m3||Yes|
|46GGB07||EL009||cf. Elephas maximus||m1||No|
To gain insight to the E. maximus population that once inhabited the Near East we sequenced mitochondrial DNA (mtDNA) extracted from four ancient teeth excavated at Kahramanmaraş (Gavur Lake Swamp) in south-central Anatolia (Figure 1A, Table 1). The elephant remains have been directly 14C dated to approximately 3500 BP (ca. 1500 BCE) and represent one of the largest assemblages of Holocene E. maximus in the region (Albayrak and Lister 2012). The dates place the remains within the Late Bronze Age, when the region around Kahramanmaraş was likely an area of influence of both the Mittani and Hittite Kingdoms (Alaura 2016). The collection is primarily made up of loose teeth but also includes partial skulls and post-cranial elements (Figure 2). Although most specimens show strong morphological affinity to modern E. maximus, a number of molars could not be identified with certainty as E. maximus compared to the Pleistocene species Palaeoloxodon antiquus (straight-tusked elephant).
Here we demonstrate the phylogenetic position of the elephant teeth (determinate and indeterminate) in relation to the known mitochondrial phylogeography of modern Asian elephant populations (Vidya, Sukumar and Melnick 2009). Further, by estimating mtDNA coalescence times we discuss the position of the Turkish remains in the broader historical context of E. maximus.
We extracted DNA from nine elephantid teeth from Kahramanmaraş and sequenced 579 bp mitochondrial DNA, including 574 bp of the fragment analysed in Vidya, Sukumar and Melnick (2009), comprising the C-terminal of cyt-b, t-RNAThr, t-RNAPro and the hypervariable left domain of the control region. The sampling, DNA extraction and all pre-PCR work were carried out in a dedicated ancient DNA laboratory at the Natural History Museum, London, following strict protocols and recommendations (Gilbert et al. 2005). No work on other elephantid DNA had taken place previously in the ancient DNA laboratory.
Teeth chosen for DNA extraction were decontaminated first by applying a 1% dilute solution of bleach (sodium hypochlorite) on and around the surface area targeted for dentine excision, followed by UV irradiation at 254 nm for 20 min at a distance of <10 cm. The outer surface (~1 mm) was removed by abrasion and tooth powder was obtained using a Dremel drill at the lowest possible rpm setting. We obtained a total of 30–60 mg powder per specimen.
DNA was extracted and purified following a modified version of a silica-binding protocol (Yang et al. 1998, Ottoni et al. 2013). Tooth powder was suspended and digested overnight at constant rotation at 50°C in 2 mL of extraction buffer consisting of 0.425 M EDTA (pH8), 1 mM Tris–HCl (pH8), 0.05% w/v SDS, and 0.33 mg/ml Proteinase K. The digested solution was concentrated to ≤100 μL using Amicon Ultra-4 centrifugal filter units with a molecular weight cut-off of 30 kDa (Amicon® Ultra, Millipore). We isolated and purified the DNA on silica columns (QIAquick® PCR purification kit, Qiagen) following the manufacturer’s protocol and eluted in 100 μL EB buffer. One negative extraction control was processed alongside the nine tooth samples.
A series of novel PCR primers were designed in Primer3Plus (Untergasser et al. 2012) to amplify short (101–119 bp) and mostly overlapping mtDNA fragments (Table 2). Singleplex PCRs were carried out in 25 μL reactions using 1.0 U Smart-Taq Hot DNA Polymerase and 1X Smart-Taq Buffer with (NH4)2SO4 (Naxo), 2.5 mM MgCl2, 200 μM of each dNTP, 0.5 mg/mL final concentration of RSA, and 0.2 μM of each primer. The cycling conditions were 95°C for 15 min followed by 55 cycles of 94°C for 20 sec, 55–57°C for 20 sec, 72°C for 20 sec, and a final extension for 5 min at 72°C. We obtained two PCR amplified products/fragment/specimen (amplified at different times) and included one negative PCR control for every eight PCRs.
|Name||Sequence 5’–3’||Length bp||Fragment||Nucleotide positions in reference to NC_005129.2 (with primers/without primers)|
|P1F||CATGAATTGGCAGYCAACCA||20||1||15,154–15,265 (103 bp)/15,174–15,243 (63 bp)|
|P2F||CTTCTCCATTATTCTAGCTTTCC||23||2||15,221–15,336 (116 bp)/15,244–15,314 (71 bp)|
|P3F||CATCAAGTAACCCCTATAGTATAAGAC||27||3||15,275–15,391 (117 bp)/15,302–15,369 (68 bp)|
|P4F||AAGGGTATTCAGGGAAGAGG||20||4||15,343–15,449 (107 bp)/15,363–15,428 (66 bp)|
|P5F||CCTCGCTATCAATACCCAAAA||21||5||15,371–15,472 (102 bp)/15,392–15,452 (61 bp)|
|P6F||TAAATGCTCGTCCCCATACA||20||6||15,452–15,552 (101 bp)/15,472–15,530 (59 bp)|
|P7F||CCATACYATGTATAATCGTGCATCA||25||7||15,512–15,629 (118 bp)/15,537–15,622 (86 bp)|
|P8F||TCAATGTGTYRAGTCATATTYBTG||24||8||15,575–15,685 (111 bp)/15,599–15,664 (66 bp)|
|P9F||TCATGGATATTRTTYRCCTACGA||23||9||15,623–15,741 (119 bp)/15,646–15,721 (76 bp)|
|P10F||AAGCTCTTGATCGTRCATAGC||21||10||15,673–15,776 (104 bp)/15,694–15,756 (63 bp)|
PCR products were purified using the QIAquick® PCR purification kit (Qiagen), pooled into equimolar ratios following quantification on a Qubit® (Invitrogen), and constructed into barcoded sequencing libraries using the TruSeq Nano DNA sample prep kit (v2, Illumina) following manufacturer’s recommendations (but omitting initial shearing) at the sequencing facility at the Natural History Museum, London. Sequencing was performed in-house on the Illumina MiSeq platform the using MiSeq Reagent Kits v2 at ~10% of a full sequencing run. Libraries were de-multiplexed on the MiSeq platform and exported in FastQ format. FastQ files were trimmed in Geneious V7 (Kearse et al. 2012) first by removing regions with more than a 5% chance of an error per base, and secondly by retaining regions with no more than 2 bases with a quality of ≤30. Sequences were aligned against a reference E. maximus sequence (Genbank accession AY245813), and consensus sequences were called using majority consensus and highest cumulative quality scores. Bases with a quality of ≤30 were ignored and called as Ns. Base coverage ranged from ca. 1000–30,000.
We downloaded reference data from two Loxodonta cyclotis sequences (JQ438507 and KJ557424), two Loxodonta africana sequences (JQ438674 and JQ438767) (Ishida et al. 2013, Finch et al. 2014), as well as two Mammuthus primigenius sequences (GU984769 and KC427898) (Nyström et al. 2010, Palkopoulou et al. 2013) and aligned them with the Kahramanmaraş sequences using MAFFT v7.017 (Katoh et al. 2002). We then aligned and pruned the DNA sequences to previously published Elephas maximus data (Fernando et al. 2000, Fernando et al. 2003, Vidya, Sukumar and Melnick 2009). Of the 579 bp amplified from the Kahramanmaraş samples we used 570 bp for the phylogenetic analysis described below. However, the Kahramanmaraş sequence data includes 32 bp missing data resulting from the difficulty in designing PCR primers to amplify short and fully overlapping amplicons.
We first assessed the phylogenetic position of the ancient Turkish specimens using MrBayes 3.2.2 (Huelsenbeck and Ronquist 2001) and a non-partitioned data set of unique haplotypes (n = 38 including Loxodonta and Mammuthus; n = 32 Elephas). We ran three heated chains (temperature = 0.4) over 10,000,000 generations with a subsampling frequency of 10,000 and a burn-in length of 1,000,000 generations. The consensus tree was visualized in Geneious 7 (Kearse et al. 2012). The nucleotide substitution model for the phylogenetic analysis of the non-partitioned dataset in MrBayes was HKY+I+G, which was the best-fitting model for all partitions but one in the Partition Finder analysis described below.
We then estimated the age of the common ancestor of the node defining the β1-subclade in which we observe the Kahramanmaraş haplotype using BEAST 2.1.3 (Bouckaert et al. 2014). Nucleotides were first grouped into five different partitions: the 1st, 2nd, and 3rd codon positions respectively for cyt-b; the tRNAThr and tRNAPro; and the control region. Nucleotide substitution models were estimated for each partition using the Bayesian Information Criterion in Partition Finder (Lanfear et al. 2012). The best-fit model for 1st and 3rd codon positions of the cyt-b gene, the tRNAThr/tRNAPro, and the control region was HKY+I+G, and the K80 model for the 2nd codon position in cyt-b. The BEAST analysis included all DNA sequences compiled for this study: Loxodonta (n = 4), Mammuthus (n = 2), Kahramanmaraş (n = 4), and Elephas maximus (n = 534).
We used tip dating on ancient DNA sequences as follows:
To obtain new MCRA estimates we used the following node age calibrations based on mitogenome divergence times estimated previously using a ‘narrow range fossil calibration’ (Brandt et al. 2012):
We then placed uniform priors on the following clades and estimated their MRCA:
We used a strict clock model under a constant population size prior. The clock rate prior was set to default (uniform distribution ranging from 0 – infinity). Three runs of 30,000,000 generations were completed, reaching effective sample sizes (ESS) of >100 (though the vast majority were >200) for each parameter. Trees and trace files were sampled/logged every 3000 generations. The trace files were also visually inspected using Tracer V1.6 (Rambaut, Suchard and Drummond 2013) to ensure proper mixing and that the Markov chain had reached stationary distribution.
We successfully amplified and sequenced DNA from four specimens, yielding a success rate of ca. 44% (4/9), consistent with previously published data on ca. 120 bp mtDNA fragments sequenced from ancient pigs from Near East dating to approximately 3500 cal. BP (Ottoni et al. 2013). We furthermore observed damaged bases (C > T or G > A transitions) across all sequenced DNA fragments, consistent with cytosine deamination to uracil which is a common type of damage in ancient DNA molecules (Gilbert et al. 2003). The negative extraction and PCR controls did not produce PCR amplified products. These observations support the authenticity of our data.
The mtDNA phylogeny of extant E. maximus consists of two highly divergent clades (α and β). The four ancient elephants from Kahramanmaraş all carried an mtDNA haplotype that groups in the β1 clade of the major β clade of E. maximus (Figure 1B). A comparison with previously published data on NCBI GenBank shows that the Kahramanmaraş sequences are identical to only a single modern elephant from Thailand (Unpublished GenBank accession JQ287727). While this haplotype appears to be rare among extant populations, the β1 sub-clade is present across their range (except Peninsular Malaysia, Borneo and Sumatra), being particularly abundant in South and Central India – and varies in frequency from 3.2–100% (Vidya, Sukumar and Melnick 2009). The lack of a clear phylogeographic structure means that it would be problematic attempting to ‘link’ the Kahramanmaraş elephants to any specific modern population. We therefore conclude only that they group within present-day variation.
The α and β clades are estimated to share a most recent common ancestor (MRCA) some 1.35 [0.9–1.84] Ma (Brandt et al. 2012). The internal MRCAs for these clades have been estimated to 0.86 Ma (0.52–1.21 Ma) and 1.58 Ma (1.28–1.86 Ma) respectively (Vidya, Sukumar and Melnick 2009). Our analysis in BEAST provides somewhat younger dates and places the MRCA of the α-clade at 0.49 Ma (0.18–0.87 Ma) and the β-clade at 0.99 Ma (0.62–1.41 Ma). However, the MRCA for all E. maximus as well as the MRCA of the α and β clades are still older than the earliest morphologically distinct E. maximus fossils, which date to ca. 0.2–0.1 Ma (Maglio 1973, Lister et al. 2013), implying that the present-day mtDNA gene pool consists of lineages that originated in earlier, ancestral Elephas species (Vidya, Sukumar and Melnick 2009). Morphology as well as geographic and temporal distribution suggests that E. hysudricus is a likely ancestor of E. maximus (Maglio 1973, Lister et al. 2013). However, the deep divergences of the α and β clades might in theory reflect genetic contribution from more than one species, such as E. hysudrindicus or even P. namadicus, in addition to E. hysudricus (Lister et al. 2013). This possibility is congruent with a recent paper (Meyer et al. 2017) that analysed mitogenomes and partial autosomal genomes of P. antiquus and suggested possible historical gene flow between that species and E. maximus. Similarly, historical gene flow between African savanna and forest elephants (L. africana and L. cyclotis) explains the observation that some savanna elephants carry the deeply divergent ‘F-clade’ mitochondrial haplotypes and not ‘S-clade’ haplotypes, despite the former being characteristic of forest elephants while the latter is private to savanna elephants (Roca, Georgiadis and O’Brien 2004).
By estimating that the haplotype carried by the Kahramanmaraş elephants and a modern Thai elephant shared a common ancestor that lived 3.7–58.7 kyr ago (95% HPD; mean 23.5 kyr) we can conclude that it very likely belonged to a population of E. maximus. This implies that range-wide population connectivity, such as gene flow or migration, has taken place at some time or times since the start of MIS 3 some 57 kya, likely reflecting range and population expansion that led to the establishment of the Bronze Age southwest Asian population.
Deraniyagala (1950) defined a subspecies, E. maximus asurus, for the Near Eastern elephants, but their supposed characteristics were based on doubtful interpretation of Bronze Age illustrations, aside from the suggestion, based on rather few bone measurements, that the animals were of unusually large size (Deraniyagala 1955, Pfälzner 2013). While the majority of molars from Kahramanmaraş are of typical E. maximus morphology, a few are more ambiguous (Albayrak and Lister 2012) (Table 1). The sequencing results allow us to confirm the taxonomic identity of one tooth (46GGB01) identified as likely E. maximus but of unusual occlusal morphology, and two (46GGB03 and 46GGB09) whose specific identity was indeterminate between E. maximus and the P. antiquus. All are confirmed as belonging to the lineage of extant E. maximus. In light of the new genetic data and the morphological homogeneity of most of the assemblage (Albayrak and Lister 2012), it is parsimonious to conclude that the Kahramanmaraş population represents a population of E. maximus that harboured greater variation in dental morphology than extant populations (Albayrak and Lister 2012). Our mtDNA evidence provides no indication of genetic divergence from East Asian populations, but further modern and ancient genetic and morphological data – preferably including nuclear DNA sequences – are required to resolve this and other questions concerning the complex evolutionary history of E. maximus (Meyer et al. 2017).
The sharing of a haplotype between the Turkish fossils and a modern individual from the Far East supports the hypothesis that the Near-Eastern elephants were the western-most expansion of the core population of India and Southeast Asia. It leaves open the possibility that the Turkish population was established only shortly before its Bronze Age date. However, the estimated age of the common ancestor of the Kahramanmaraş haplotype, extending to 58.7 kyr at its 95% lower bound, means that the Near Eastern populations could alternatively have been established, or at least separated from other populations in southern Asia, as long ago as MIS 3, assuming that the Turkish population and their direct ancestors were isolated without gene flow. The split time between the population ancestral to the Kahramanmaraş elephants and the eastern core populations could in theory have been even longer in the presence of gene flow. The genetic data do not distinguish between these possibilities, nor therefore between natural range expansion and human introduction of elephants into the Near East. The latter has been posited mainly due to the lack of well-dated skeletal remains of elephants before 3500 BP, or ivory (which could have been humanly transported) before 4000 BP (Çakirlar and Ikram 2016). The use of captive elephants likely did not begin in India until ca. 4600 BP (Sukumar 2011), and deliberate long-distance movement of elephants has been thought impracticable in the Bronze Age (Pfälzner 2013). Although there are historical/pictorial sources referring to the movement of live elephants in the Bronze and Iron Ages in southwest Asia (Collon 1977), the westward movement of elephants from India for war or other purposes is first clearly recorded around ca. 2600 BP (Sukumar 2011). Of considerable significance is that the Kahramanmaraş assemblage, although not systematically excavated, appears very likely to represent a wild-living population. It comprises various elements of the skeleton of multiple individuals, some of them probably associated, with no cut marks, artefacts, or other sign of human activity (Alaura 2016, Yar et al. 2016) in stark contrast to Near-Eastern archaeological sites with occasional elephant remains (Çakirlar and Ikram 2016). This suggests a naturally established population, or conceivably a wild-living population founded by humanly-introduced animals.
Combining these data suggests that the most likely scenario is the westward expansion of the elephants’ range, presumably during favourable climatic episode(s), i.e. with sufficient warmth and moisture to provide the required vegetation and drinking water to support the animals along the route. These requirements would have been necessary whether the population arrived by natural expansion or human agency. Rowe et al. (2012) show cyclicity of wet and dry episodes in Turkey extending back into the Late Pleistocene, as do Rajagopalan et al. (1997) for the Indian subcontinent. Corresponding cyclical variation in the range of mammal species, including elephants, is likely. More recent studies show that the climate of the Near and Middle East underwent a shift from wetter conditions in the early Holocene to drier conditions from around 6500 BP, but that this was interrupted by an interval of moister conditions between around 5000–3500 BP (Finné et al. 2011). For example, Migowski et al. (2006) indicate a wet phase in the Levant ca. 5600–3500 BP, based on Dead Sea lake level changes. A variety of data from Lake Van in southeastern Turkey and Soreq Cave in Israel suggests that a moister climate prevailed in the region during the Chalcolithic and Early Bronze Age, although the region became significantly more arid from around 2000 cal. BC (ca. 4000 BP) (Bar-Matthews et al. 1999, Wick, Lemcke and Sturm 2003). Together, these data suggests that wetter climatic periods suitable for the immigration of elephants did prevail in the region at intervals during the Holocene, although the general climatic trend is one of increasing aridity.
Cyclical expansion and contraction of range is also supported by earlier fossil records of Elephas in the Near East, where the earliest dated finds are an Early Pleistocene molar tooth of Elephas sp. from Evron Quarry in Israel (Tchernov et al. 1994), and five Middle Pleistocene teeth of Elephas cf. hysudricus from Ma’ayan Baruch in Israel and ‘Ain Soda in Jordan (Lister et al. 2013).
The consensus DNA sequences are publically available on GenBank (accession numbers MF314177-MF314180).
Raman Sukumar and Liora Kolska Horwitz are thanked for very helpful discussion.
The study was supported by the Natural History Museum (SIF fund).
The authors have no competing interests to declare.
Albayrak, E and Lister, A M 2012 Dental remains of fossil elephants from Turkey. Quaternary International, 276: 198–211. DOI: https://doi.org/10.1016/j.quaint.2011.05.042
Bar-Matthews, M, Ayalon, A, Kaufman, A and Wasserburg, G J 1999 The Eastern Mediterranean paleoclimate as a reflection of regional events: Soreq cave, Israel. Earth and Planetary Science Letters, 166(1): 85–95. DOI: https://doi.org/10.1016/S0012-821X(98)00275-1
Bouckaert, R, Heled, J, Kühnert, D, Vaughan, T, Wu, C-H, Xie, D, Suchard, M A, Rambaut, A and Drummond, A J 2014 BEAST 2: A Software Platform for Bayesian Evolutionary Analysis. PLoS Comput Biol, 10(4): e1003537. DOI: https://doi.org/10.1371/journal.pcbi.1003537
Brandt, A L, Ishida, Y, Georgiadis, N J and Roca, A L 2012 Forest elephant mitochondrial genomes reveal that elephantid diversification in Africa tracked climate transitions. Molecular Ecology, 21(5): 1175–1189. DOI: https://doi.org/10.1111/j.1365-294X.2012.05461.x
Çakirlar, C and Ikram, S 2016 ‘When elephants battle, the grass suffers’. Power, ivory and the Syrian elephant. Levant, 48: 167–183. DOI: https://doi.org/10.1080/00758914.2016.1198068
Choudhury, A, Lahiri Choudhury, D K, Desai, A, Duckworth, J W, Easa, P S, Johnsingh, A J T, Fernando, P, Hedges, S, Gunawardena, M, Kurt, F, Karanth, U, Lister, A, Menon, V, Riddle, H, Rübel, A and Wikramanayake, E 2008 Elephas maximus, Asian Elephant, IUCN.
Collon, D 1977 Ivory. Iraq, 39(2): 219–222. DOI: https://doi.org/10.2307/4200069
Fernando, P, Pfrender, M, Encalada, S and Lande, R 2000 Mitochondrial DNA variation, phylogeography and population structure of the Asian elephant. Heredity, 84: 362–372. DOI: https://doi.org/10.1046/j.1365-2540.2000.00674.x
Fernando, P, Vidya, T N C, Payne, J, Stuewe, M, Davison, G, Alfred, R J, Andau, P, Bosi, E, Kilbourn, A and Melnick, D J 2003 DNA Analysis Indicates That Asian Elephants Are Native to Borneo and Are Therefore a High Priority for Conservation. PLOS Biology, 1(1): e6. DOI: https://doi.org/10.1371/journal.pbio.0000006
Finch, T M, Zhao, N, Korkin, D, Frederick, K H and Eggert, L S 2014 Evidence of Positive Selection in Mitochondrial Complexes I and V of the African Elephant. PLOS ONE, 9(4): e92587. DOI: https://doi.org/10.1371/journal.pone.0092587
Finné, M, Holmgren, K, Sundqvist, H S, Weiberg, E and Lindblom, M 2011 Climate in the eastern Mediterranean, and adjacent regions, during the past 6000 years – A review. Journal of Archaeological Science, 38(12): 3153–3173. DOI: https://doi.org/10.1016/j.jas.2011.05.007
Gilbert, M T P, Bandelt, H J, Hofreiter, M and Barnes, I 2005 Assessing ancient DNA studies. Trends in Ecology & Evolution, 20(10): 541–544. DOI: https://doi.org/10.1016/j.tree.2005.07.005
Gilbert, M T P, Hansen, A J, Willerslev, E, Rudbeck, L, Barnes, I, Lynnerup, N and Cooper, A 2003 Characterization of genetic miscoding lesions caused by postmortem damage. Am. J. Hum. Genet, 72: 48–61. DOI: https://doi.org/10.1086/345379
Huelsenbeck, J and Ronquist, F 2001 MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics, 17(8): 754–755. DOI: https://doi.org/10.1093/bioinformatics/17.8.754
Ishida, Y, Georgiadis, N J, Hondo, T and Roca, A L 2013 Triangulating the provenance of African elephants using mitochondrial DNA. Evolutionary Applications, 6(2): 253–265. DOI: https://doi.org/10.1111/j.1752-4571.2012.00286.x
Katoh, K K, Misawa, K, Kuma, K and Miyata, T 2002 MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research, 30: 3059–3066. DOI: https://doi.org/10.1093/nar/gkf436
Kearse, M, Moir, R, Wilson, A, Stones-Havas, S, Cheung, M, Sturrock, S, Buxton, S, Cooper, A, Markowitz, S, Duran, C, Thierer, T, Ashton, B, Mentjies, P and Drummond, A 2012 Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics, 28(12): 1647–1649. DOI: https://doi.org/10.1093/bioinformatics/bts199
Lanfear, R, Calcott, B, Ho, S Y W and Guindon, S 2012 PartitionFinder: Combined Selection of Partitioning Schemes and Substitution Models for Phylogenetic Analyses. Molecular Biology and Evolution, 29(6): 1695–1701. DOI: https://doi.org/10.1093/molbev/mss020
Lister, A M, Dirks, W, Assaf, A, Chazan, M, Goldberg, P, Applbaum, Y H, Greenbaum, N and Horwitz, L K 2013 New fossil remains of Elephas from the southern Levant: Implications for the evolutionary history of the Asian elephant. Palaeogeography Palaeoclimatology Palaeoecology, 386: 119–130. DOI: https://doi.org/10.1016/j.palaeo.2013.05.013
Maglio, V J 1973 Origin and evolution of the Elephantidae. Transactions Am phil Soc, 63(3): 1–149. DOI: https://doi.org/10.2307/1006229
Meyer, M, Palkopoulou, E, Baleka, S, Stiller, M, Penkman, K E H, Alt, K W, Ishida, Y, Mania, D, Mallick, S, Meijer, T, Meller, H, Nagel, S, Nickel, B, Ostritz, S, Rohland, N, Schauer, K, Schüler, T, Roca, A L, Reich, D, Shapiro, B and Hofreiter, M 2017 Palaeogenomes of Eurasian straight-tusked elephants challenge the current view of elephant evolution. eLife, 6: e25413. DOI: https://doi.org/10.7554/eLife.25413
Migowski, C, Stein, M, Prasad, S, Negendank, J F W and Agnon, A 2006 Holocene climate variability and cultural evolution in the Near East from the Dead Sea sedimentary record. Quaternary Research, 66(3): 421–431. DOI: https://doi.org/10.1016/j.yqres.2006.06.010
Nyström, V, Dalén, L, Vartanyan, S, Lidén, K, Ryman, N and Angerbjörn, A 2010 Temporal genetic change in the last remaining population of woolly mammoth. Proceedings of the Royal Society B: Biological Sciences. DOI: https://doi.org/10.1098/rspb.2010.0301
Olivier, R 1978 Distribution and Status of the Asian Elephant. Oryx, 14(4): 379–424. DOI: https://doi.org/10.1017/S003060530001601X
Ottoni, C, Girdland Flink, L, Evin, A, Geörg, C, De Cupere, B, Van Neer, W, Bartosiewicz, L, Linderholm, A, Barnett, R, Peters, J, Decorte, R, Waelkens, M, Vanderheyden, N, Ricaut, F-X, Çakırlar, C, Çevik, Ö, Hoelzel, A R, Mashkour, M, Mohaseb Karimlu, A F, Sheikhi Seno, S, Daujat, J, Brock, F, Pinhasi, R, Hongo, H, Perez-Enciso, M, Rasmussen, M, Frantz, L, Megens, H-J, Crooijmans, R, Groenen, M, Arbuckle, B, Benecke, N, Strand Vidarsdottir, U, Burger, J, Cucchi, T, Dobney, K and Larson, G 2013 Pig Domestication and Human-Mediated Dispersal in Western Eurasia Revealed through Ancient DNA and Geometric Morphometrics. Molecular Biology and Evolution, 30(4): 824–832. DOI: https://doi.org/10.1093/molbev/mss261
Palkopoulou, E, Dalén, L, Lister, A M, Vartanyan, S, Sablin, M, Sher, A, Edmark, V N, Brandström, M D, Germonpré, M, Barnes, I and Thomas, J A 2013 Holarctic genetic structure and range dynamics in the woolly mammoth. Proceedings of the Royal Society B: Biological Sciences, 280(1770). DOI: https://doi.org/10.1098/rspb.2013.1910
Pfälzner, P 2013 The Elephant Hunters of Bronze Age Syria. Cultures in Contact: From Mesopotamia to the Mediterranean in the Second Millennium, Aruz, B C J, Graff, S B and Rakic, Y (eds.), 112–131. New York, Metropolitan Museum of Art (Distributed by Yale University Press).
Rajagopalan, G, Sukumar, R, Ramesh, R, Pant, R K and Rajagopalan, G 1997 Late Quaternary vegetational and climatic changes from tropical peats in southern India – an extended record up to 40,000 years BP. Current Science, 73: 60–63.
Rambaut, A, Suchard, M and Drummond, A J 2013 Tracer v1.4: MCMC trace analyses tool. Available from: http://tree.bio.ed.ac.uk/software/tracer/.
Roca, A L, Georgiadis, N and O’Brien, S J 2004 Cytonuclear genomic dissociation in African elephant species. Nature Genetics, 37: 96. DOI: https://doi.org/10.1038/ng1485
Rowe, P J, Mason, J E, Andrews, J E, Marca, A D, Thomas, L, van Calsteren, P, Jex, C N, Vonhof, H B and Al-Omari, S 2012 Speleothem isotopic evidence of winter rainfall variability in northeast Turkey between 77 and 6 ka. Quaternary Science Reviews, 45: 60–72. DOI: https://doi.org/10.1016/j.quascirev.2012.04.013
Tchernov, E, Horwitz, L K, Ronen, A and Lister, A 1994 The faunal remains from Evron Quarry in relation to other Lower Paleolithic hominid sites in the southern Levant. Quaternary Research, 42: 328–339. DOI: https://doi.org/10.1006/qres.1994.1083
Turvey, S T, Tong, H, Stuart, A J and Lister, A M 2013 Holocene survival of late Pleistocene megafauna in China: a critical review of the evidence. Quaternary Science Reviews, 76: 156–166. DOI: https://doi.org/10.1016/j.quascirev.2013.06.030
Untergasser, A, Cutcutache, I, Koressaar, T, Ye, J, Faircloth, B C, Remm, M and Rozen, S G 2012 Primer3-new capabilities and interfaces. Nucleic Acids Research, 40(15). DOI: https://doi.org/10.1093/nar/gks596
Vidya, T N C, Sukumar, R and Melnick, D J 2009 Range-wide mtDNA phylogeography yields insights into the origins of Asian elephants. Proceedings of the Royal Society B-Biological Sciences, 276(1658): 893–902. DOI: https://doi.org/10.1098/rspb.2008.1494
Wick, L, Lemcke, G and Sturm, M 2003 Evidence of Lateglacial and Holocene climatic change and human impact in eastern Anatolia: high-resolution pollen, charcoal, isotopic and geochemical records from the laminated sediments of Lake Van, Turkey. The Holocene, 13(5): 665–675. DOI: https://doi.org/10.1191/0959683603hl653rp
Yang, D Y, Eng, B, Waye, J S, Dudar, J C and Saunders, S R 1998 Technical note: Improved DNA extraction from ancient bones using silica-based spin columns. American Journal of Physical Anthropology, 105(4): 539–543.