The straight-tusked elephants Palaeoloxodon spp. were widespread across Eurasia during the Pleistocene. Phylogenetic reconstructions using morphological traits have grouped them with Asian elephants (Elephas maximus), and many paleontologists place Palaeoloxodon within Elephas. Here, we report the recovery of full mitochondrial genomes from four and partial nuclear genomes from two P. antiquus fossils. These fossils were collected at two sites in Germany, Neumark-Nord and Weimar-Ehringsdorf, and likely date to interglacial periods ~120 and ~244 thousand years ago, respectively. Unexpectedly, nuclear and mitochondrial DNA analyses suggest that P. antiquus was a close relative of extant African forest elephants (Loxodonta cyclotis). Species previously referred to Palaeoloxodon are thus most parsimoniously explained as having diverged from the lineage of Loxodonta, indicating that Loxodonta has not been constrained to Africa. Our results demonstrate that the current picture of elephant evolution is in need of substantial revision.https://doi.org/10.7554/eLife.25413.001
Understanding how extinct species are related to each other or to their living relatives is often a difficult task. Many extinct species have been identified only from incomplete fragments of some of their bones. However, even if complete skeletons have been found, determining the relationships between species can be tricky because researchers often have to rely solely on the shapes of the bones.
It is sometimes possible to retrieve DNA sequences from fossil bones. This is easier with younger fossils and those that have been recovered from cold environments. Ancient DNA sequences have been retrieved from only a few fossils older than 100,000 years, but such DNA sequences can be tremendously useful in determining how different species are related to each other.
Today there are three living elephant species: the African forest elephant, the African savanna elephant and the Asian elephant. However, there are many extinct elephant species. For example, the European straight-tusked elephant went extinct at least 30,000 years ago, although most of the fossils that have been discovered are at least 100,000 years old. Straight-tusked elephants are generally assumed to be closely related to the Asian elephant, but this conclusion had been based solely on reconstructing skeletons.
Meyer et al. have now obtained DNA sequences from fossils of four straight-tusked elephants ranging from around 120,000 to 240,000 years in age. These sequences were analysed to determine how straight-tusked elephants are related to the three living elephant species and the extinct mammoth, the DNA sequences for which can be found in public databases. The analyses revealed that straight-tusked elephants are in fact most closely related to the African forest elephant, not the Asian elephant as previously thought.
This result completely changes our picture of elephant evolution and suggests that it is extremely difficult to determine elephant relationships based on the shape of their skeleton alone. It also shows that the African elephant lineage was not restricted to the African continent (the place where all elephant lineages originated), but that it also left Africa.
Overall, the results presented by Meyer et al. confirm that DNA sequences are of critical importance for understanding the evolution of animals. Future research should include obtaining DNA sequences from additional extinct elephant species as well as careful re-evaluation of skeletal measurements for reconstructing elephant evolution.https://doi.org/10.7554/eLife.25413.002
In the late Miocene in Africa, the last of several major radiations within Proboscidea gave rise to the family Elephantidae, which comprises living elephants and their extinct relatives including mammoths (genus Mammuthus) and various dwarf elephant species from Mediterranean islands. The three living elephant species (the African savanna elephant, Loxodonta africana, the African forest elephant, L. cyclotis and the Asian elephant, Elephas maximus), represent the last remnants of this family and of the formerly much more widely distributed and species-rich order Proboscidea. Apart from mammoths, the elephant genus with the most abundant fossil record in Eurasia is Palaeoloxodon (straight-tusked elephants; Figure 1), which appears in Eurasia around 0.75 million years ago (Ma) (Lister, 2016). Based on morphological analyses, Palaeoloxodon is widely accepted as being more closely related to the extant Asian elephant than to mammoths or extant African elephants (Shoshani et al., 2007; Todd, 2010) and is often subsumed into the genus Elephas (Maglio, 1973; Sanders et al., 2010). Across its range from Western Europe to Japan, Palaeoloxodon probably comprised several species (Shoshani et al., 2007), and, based on morphological comparisons, all of them are considered to be derived from the African Palaeoloxodon (or Elephas) recki (Maglio, 1973; Saegusa and Gilbert, 2008), which was the predominant proboscidean lineage in Africa during the Pliocene and Pleistocene but went extinct around 100 thousand years ago (ka) (Owen-Smith, 2013). Straight-tusked elephants may have survived in mainland Eurasia until around 35 ka, although the youngest reliably dated remains are from the last interglacial, 115–130 ka (Stuart, 2005).
Recent technological progress has pushed back the temporal limit of ancient DNA research, enabling, for example, recovery of a low coverage genome of a ~700,000 year-old horse preserved in permafrost (Orlando et al., 2013). For more temperate regions, however, evidence of DNA preservation reaching far beyond the last glacial period is still limited to a single locality, Sima de los Huesos in Spain, where DNA has been recovered from ~430 ka old hominin and bear remains (Dabney et al., 2013; Meyer et al., 2016). While genetic analyses of the extinct interglacial fauna remain a challenging undertaking, recent advances in ancient DNA extraction (Dabney et al., 2013) and sequencing library construction (Meyer et al., 2012) have improved access to highly degraded DNA.
To better understand the evolutionary relationships between the extinct straight-tusked elephants and other elephant species, we attempted DNA extraction and sequencing from several P. antiquus fossils, four of which we investigated in depth. Three of these, which were all unambiguously assigned to P. antiquus based on their morphology, were from Neumark-Nord (NN) 1 in Germany, a fossil-rich site that has been proposed to date to MIS 5e (~120 ka) or MIS 7 (~244 ka) or both (Mania, 2010; Schüler, 2010; Penkman, 2010). This site has yielded one of the largest collections of P. antiquus remains known to date. The fourth fossil was recovered during recent active mining in the travertine deposits of Weimar-Ehringsdorf (WE), Germany, a quarry that has for more than a century yielded a rich collection of fossils representing a typical European interglacial fauna (Kahlke, 1975). Weimar-Ehringsdorf is best known for the discovery of Neanderthal remains in the early 20th century, and the assemblage is dated to MIS 7 (Mallick and Frank, 2002). The Paleoloxodon bone fragment from Weimar-Ehringsdorf is morphologically undiagnostic with respect to species. However, it was found in the Lower Travertine, which was dated to ~233 ka (Schüler, 2003) and where P. antiquus is the only elephantid found so far. We performed DNA extraction, library preparation, hybridization capture and high-throughput sequencing on all four fossils (Figure 2—source data 1) and obtained full mitochondrial genome sequences for all of them (Figure 2—figure supplement 1). All sequences show short fragment lengths (Figure 2—figure supplement 2) and signals of cytosine deamination compatible with the old age of the specimens (Figure 2—figure supplement 3).
We inferred a phylogeny using the four Paleoloxodon mitochondrial genomes and mitochondrial genomes from 16 M. primigenius, 2 E. maximus and 13 Loxodonta individuals. The latter were chosen for a diversity of haplotypes, including forest elephant derived (‘F-clade’) haplotypes as well as ‘S-clade’ haplotypes found only among savanna elephants (Debruyne, 2005). For calibration, we used an estimated divergence of the African elephant lineage from that of Asian elephants and mammoths of 6.6–8.6 Ma (Rohland et al., 2007). Surprisingly, P. antiquus did not cluster with E. maximus, as hypothesized from morphological analyses. Instead, it fell within the mito-genetic diversity of extant L. cyclotis, with very high statistical support (Figure 2). The four straight-tusked elephants did not cluster together within this mitochondrial clade, but formed two separate lineages that share a common ancestor with an extant L. cyclotis lineage 0.7–1.6 Ma (NN) and 1.5–3.0 Ma (WE) ago, respectively. However, mitochondrial DNA represents a single, maternally inherited locus and does not reflect the full evolutionary history of populations or species. Furthermore, the transfer of mitochondrial DNA between hybridizing species is not unusual when gene flow is strongly male-mediated (Petit and Excoffier, 2009; Li et al., 2016; Cahill et al., 2013), as is the case with elephants. For example, mitochondrial sequences of the F-clade have also been found in some L. africana individuals (Debruyne, 2005) despite the very substantial divergence of their nuclear genomes (Roca et al., 2005; Rohland et al., 2010), a pattern that has been attributed to mitochondrial gene flow from forest to savanna elephants (Roca et al., 2005).
We therefore performed shotgun sequencing of DNA libraries prepared from the two best-preserved NN individuals (a petrous bone and a molar) to recover nuclear DNA sequences. When mapped to the L. africana reference genome, 39% and 28% of the sequence reads generated from these specimens were identified as elephant, respectively (Figure 2—source data 1). A neighbor-joining phylogenetic tree based on ~770M (petrous) and ~210M (molar) base pairs of P. antiquus nuclear DNA placed P. antiquus and L. cyclotis as sister taxa to the exclusion of L. africana (Figure 2). A tree with identical topology was obtained using coding sequences only and a maximum likelihood approach (Figure 2—figure supplement 4). Despite the high sequence error rates associated with the low-coverage genomes generated from the P. antiquus specimens, all nodes in the nuclear trees show maximal bootstrap support. The mitochondrial and nuclear phylogenies thus support a sister group relationship between P. antiquus and L. cyclotis.
Despite their geographical proximity, the WE and NN specimens are found in different positions in the mitochondrial tree. Given that the three NN specimens show highly similar mitogenome sequences, we considered whether the sites date to different interglacials. Electron Spin Resonance (ESR) dating of tooth enamel has been applied to both sites, suggesting an age of ~117 ka (range 97–142 ka [Schüler, 2010]) for the NN1 layers from which our samples originate and of 233 ka (range 216–250 ka [Schüler, 2003]) for the WE specimen. In order to better estimate the age of the NN1 site, we used amino acid racemization of snail opercula (Penkman et al., 2011), including samples from the continental Eemian type-site of Amersfoort (Zagwijn, 1961) (Figure 2—source data 1), which is correlated with MIS 5e (Cleveringa et al., 2000). The NN1 opercula show similar (perhaps slightly lower) levels of biomolecular degradation compared to Amersfoort, suggesting an Eemian age for NN1 (Figure 2—figure supplement 5). Importantly, NN1 shows very similar levels of amino acid racemization in intra-crystalline protein decomposition as a second site at Neumark-Nord (NN2), indicating that both sites are of the same age. Considerable evidence supports an Eemian age for NN2, including palaeomagnetic data that shows a correlation with the MIS 5e Blake event and thermoluminescence dating of flint to ~126 ± 6 ka (Sier et al., 2011). These results therefore indicate an Eemian age also for NN1. Since the WE specimen likely dates to the previous interglacial, this suggests that the very different mitogenomes between WE and NN1 may reflect the contraction and re-expansion of the range of P. antiquus across glacial cycles.
Our results have implications both for understanding elephant evolutionary history and for the use of morphological data to decipher phylogenetic relationships among elephants. The strongly supported mitochondrial and nuclear DNA phylogenies clearly demonstrate that Palaeoloxodon antiquus is more closely related to Loxodonta than to Elephas (Figure 3), suggesting that Elephas antiquus should not be used synonymously for Paleoloxodon antiquus when referring to the taxon. The new phylogeny suggests a remarkable degree of evolutionary transformation, from an ancestor that possessed the features of the cranium and dentition of Loxodonta, shared by both L. africana and L. cyclotis (Maglio, 1973; Sanders et al., 2010), to a descendant that is highly similar to Elephas (sensu stricto, the lineage of the Asian elephant) in many morphological features. However, it should be noted that currently available genomic data from elephantids only allow for reconstructing the broad picture of elephant evolution. More complex evolutionary scenarios are conceivable, which might explain the presence of some Elephas-like traits in P. antiquus. These could for example involve gene flow, as has been shown for L. africana and L. cyclotis based on mitochondrial evidence (Roca et al., 2005). In addition, the very large effective population size of the forest elephants (Rohland et al., 2010) could have allowed the retention of ancestral traits by incomplete lineage sorting.
In summary, the molecular results presented here urge for a re-examination of morphology across the Elephantidae. This is especially important as the fossil record for elephants dates back several million years, well beyond the survival of ancient DNA. If, for example, P. recki, which was the most abundant Pleistocene elephant species in Africa, is indeed ancestral to P. antiquus and thus also represents a member of the Loxodonta lineage, the interpretation of the fossil record of elephantids in Africa is in strong need of revision. Furthermore, in contrast to the genera Mammuthus and Elephas, which also had their origin in Africa, the lineage of Loxodonta is generally assumed never to have left Africa. Although Osborn, 1942 placed Palaeoloxodon in the Loxodontinae on the basis of several cranial characters, later authors (Shoshani et al., 2007; Todd, 2010) have rejected this placement in favor of a placement in Elephantinae, restricting Loxodonta (and Loxodontinae) geographically to Africa. However, our data reveal that the Loxodonta lineage (as Paleoloxodon) also colonized the Eurasian continent. Last, the finding that L. africana is genetically more distant from L. cyclotis than is P. antiquus strongly supports previous evidence that urged recognition of L. cyclotis and L. africana as distinct species and underlines the importance of conservation efforts directed toward African forest elephants.
In January 2014, a fragment of an elephant long bone was discovered during work at the Ehringsdorf quarries. The specimen was removed from the lower travertine (~3m above the base and 2.5 m below the Pariser horizon), which has been dated by micro probe U/Th-series dating of primary travertine (Mallick and Frank, 2002) and ESR dating of tooth enamel (Schüler, 2003) to ~233 ka. A piece of the bone (inventory number 14/18–1) was transferred to the ancient DNA laboratory at the MPI-EVA in Leipzig. Ten extracts were prepared using between 36 and 53 mg of bone (totaling 425 mg) material following the method of Dabney et al. 2013 (Dabney et al., 2013). From these extracts, 30 libraries were prepared using single-stranded library preparation (Gansauge and Meyer, 2013) with input volumes of 4, 8 and 12 µl DNA extract (of 25 µl extract volume), respectively. The number of library molecules was determined by digital droplet PCR using Bio-Rad's (Hercules, CA) QX200 system with EvaGreen chemistry (QX200 ddPCR EvaGreen Supermix, Bio-Rad) and primers IS7 and IS8 (Meyer and Kircher, 2010) following the manufacturer’s instructions (Figure 2—source data 1). Libraries were then amplified using AccuPrime Pfx DNA polymerase (Thermo Fisher Scientific, Waltham, MA) (Dabney and Meyer, 2012) and labeled with two sample-specific indices (Kircher et al., 2012).
Ancient DNA work on the Neumark-Nord specimens was carried out in the ancient DNA laboratory at the University of Potsdam. Initially, eight specimens were screened for the presence of elephant DNA of which the two best preserved ones (individual 23, Landesmuseum Halle museum inventory number HK 2007:25.285,117; a molar fragment [NEU2A] and a fragmentary upper jaw [NEU8B]; inventory number HK 92:990) were selected for further analyses. In addition, in 2013, the petrous bone of individual 30 (NEPEC; inventory number HK 2007:25:280 = E15.1.96) was sampled and also used in the analysis. Fourteen DNA extracts were prepared from individual 30, and six from each of the two other specimens, using approximately 50 mg of bone powder in each extraction. DNA extraction and library preparation were performed as described above but including Archaeoglobus fulgidus uracil-DNA glycosylase in library preparation (Gansauge and Meyer, 2013), which removes the majority of uracils that are typically present in ancient DNA fragments. In addition, to maximize yields in library preparation, two extracts (25 µl each) from each specimen were combined and 40 µl were used as input for library preparation. Reaction volumes in steps 1–3 of the protocol were doubled to accommodate larger input volumes of extract. Optimal amplification cycle numbers were established using qPCR (PikoReal Real-Time PCR system, Thermo Fisher Scientific) with primers IS7 and IS8 (Gansauge and Meyer, 2013). Libraries were then amplified and labeled with one sample-specific index. After purification (MinElute PCR purification kit, Qiagen, Germany), the different libraries for each sample were pooled.
52-mer capture probes for the enrichment of mtDNA sequences from elephants were designed using the published mtDNA genome sequences of African forest elephant (NC_020759), Asian elephant (NC_005129), African savanna elephant (NC_000934) and the mastodon (NC_009574), with one probe starting at each position in these genomes. Probes containing simple repeats longer than 24 bp (repetition of the same 1–8 bp sequence motif) were removed. Single-stranded biotinylated DNA probes were generated as described elsewhere (Fu et al., 2013) and used for two successive rounds of hybridization capture following a bead-based protocol (Maricic et al., 2010). Enriched libraries were pooled and sequenced on one lane of a HiSeq2500 (Illumina, San Diego, CA) in paired-end mode (2 × 76 cycles plus 2 × 7 cycles index reads; Weimar-Ehringsdorf libraries) or on an Illumina NextSeq 500 (2 × 76 cycles plus 1 × 8 cycles index read; Neumark-Nord libraries).
Sequences were assigned to their source library requiring perfect matches to one of the expected indices or index pairs and overlap-merged to reconstruct full-length molecule sequences (Renaud et al., 2014). Due to the different properties of the data obtained from Weimar-Ehringsdorf and Neumark-Nord with regard to sequence length distribution and damage patterns (Figure 2—figure supplements 2 and 3), two different strategies were used for mapping and consensus calling. To minimize the loss of alignments due to the high frequencies of damaged-induced substitutions in the Weimar-Ehringsdorf data, mapping to the L. cyclotis mtDNA genome (JN673264) was performed as previously described for the Sima de los Huesos mtDNA assemblies (Dabney et al., 2013), using BWA and allowing up to five C to T substitutions but not more than three of other types. The sequences from Neumark-Nord were mapped with ‘ancient’ parameters as described elsewhere (Meyer et al., 2012). PCR duplicates were removed with bam-rmdup (Stenzel, 2014 Biohazard, available from https://bitbucket.org/ustenzel/biohazard) by calling a consensus from sequences with identical alignment start and end coordinates. Sequences shorter than 30 bp were discarded. An overview of the DNA extracts, libraries and sequences generated in this study is provided in Figure 2—source data 1.
When visually inspecting the Weimar-Ehringsdorf sequence alignments, we identified several regions in the mitochondrial genome where more than one sequence variant was present. Based on BLAST searches on a subset of these sequences, we found that they derived from present-day human or microbial contamination. We thus aligned all sequences to the identified contaminant genomes (GenBank accession nos. NC_012920, AF365635 and CP008889) and removed sequences that showed a greater similarity to one of the contaminants than to the African forest elephant mtDNA. No removal of contaminant sequences was necessary for the Neumark-Nord samples. To minimize the impact of damage-derived C to T substitutions on consensus calling, all T occurring in the first and last three positions of the Weimar-Ehringsdorf sequences were substituted by N. Next, a position-based tabular output was generated from the alignment files using the ‘mpileup’ function of SAMtools (Li et al., 2009). This file was used to call the consensus at positions with minimum sequence coverage of 3 if the sequences were in at least 67.0% agreement. At three positions in the mtDNA genome (positions 384, 8467, 8469) with low consensus support, we spotted obvious alignment errors in one or all specimens and determined the consensus base manually. Apart from a ~500 bp stretch of repetitive sequence in the D-loop, which cannot be reconstructed with short DNA fragments, only four positions remain undetermined in the Weimar-Ehringsdorf sequence and even fewer (between none and three) in the Neumark-Nord sequences.
We estimated mitochondrial phylogenies using the software BEAST (Drummond et al., 2012) v 1.8.2 and a data set including 31 complete mitochondrial genomes (GenBank accession nos.; L. cyclotis: JN673264, JN673263, KJ557424, KJ557423, KY616976, KY616979, KY616978; L. africana: WA4020, KR0014, KR0138, NC000934, DQ316069, AB443879; E. maximus: NC005129, DQ316068; M. primigenius: DQ316067, NC007596, EU155210, EU153449, EU153455, EU153456, EU153458, EU153445, EU153446, EU153447, EU153448, EU153452, EU153453, EU153454, JF912200; M. columbi: JF912199). For three of the L. cyclotis individuals (LO3505, LO3508 and DS1511) and three L. africana individuals (WA4020, KR0014, KR0138), only partial mitochondrial sequences were previously published. Full genome sequences were obtained using previously collected samples (Ishida et al., 2013) and the amplification and sequencing strategy detailed by Brandt et al. (2012), except that additional primers were used in sequencing (Figure 2—source data 1). The complete mitochondrial genome sequences were partitioned prior to analysis into four partitions, representing concatenated genes (with ND6 reversed), tRNAs, rRNAs, and the control region, and analyses were performed with and without the control region fragment. All BEAST analyses were performed assuming the flexible skygrid coalescent model (Gill et al., 2013) and the uncorrelated lognormal relaxed molecular clock (Drummond et al., 2006). We calibrated the molecular clock using the ages of ancient tips and a lognormal prior with a mean of 7.6 million years and standard deviation of 500,000 years for the divergence of the Loxodonta and Elephas/Mammuthus lineages (Rohland et al., 2007). Ages of the ancient samples were sampled from normal distributions derived from stratigraphic and previously estimated radiometric dates: Neumark-Nord: 142–92 ka (Schüler, 2010); Ehringsdorf 250–216 ka (Mallick and Frank, 2002). Separate evolutionary rates and models of nucleotide substitution, as estimated using jModelTest (Posada, 2008), were estimated for each partition in the alignment. We ran two MCMC chains for 60 million iterations each, with trees and model parameters sampled every 6000 iterations. Chain convergence and parameter sampling were examined by eye using Tracer v 1.6 (Rambaut A, Suchard MA, Xie D & Drummond AJ (2014) Tracer v1.6, available from http://beast.bio.ed.ac.uk/Tracer). The first 10% of samples were discarded from each run after which the two runs were combined. Trees were summarized and maximum clade credibility (MCC) trees identified using TreeAnnotator v 1.8.2, which is distributed as part of the BEAST package. MCC trees were edited and annotated using FigTree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/).
Libraries from the three Neumark-Nord samples and another sample (not included in this study) were pooled in equimolar concentrations and shotgun-sequenced on an Illumina NextSeq 500 (2 × 76 bp cycles) at Harvard Medical School. Following determination of endogenous content and complexity in each library, two of them (NEPEC from the petrous bone and NEU2A from the molar fragment) were chosen for additional sequencing and were pooled together with another sample (not included in this study) for one NextSeq 500 run.
Sequences were assigned to their source library according to their index allowing for one mismatch. Adapters were trimmed and paired-end sequences were merged with SeqPrep 1.1 (https://github.com/jstjohn/SeqPrep) using default parameters but with a modification in the source code to retain the best base quality scores in the merged region. Merged sequences shorter than 30 bp were discarded. Alignment to the African savanna elephant reference genome (loxAfr4; downloaded from ftp://ftp.broadinstitute.org/distribution/assemblies/mammals/elephant/loxAfr4/) was performed with BWA’s version 0.7.8 (Li and Durbin, 2009) using ‘ancient’ parameters and SAMtools’ v.0.1.19 ‘samse’ command (Li et al., 2009). A custom script was used to remove duplicates, which takes into account the alignment coordinates of both ends of the sorted sequences and their orientation.
From the first sequencing run, 47% and 35% of the sequences from the libraries from the petrous bone and the molar fragment (NEPEC and NEU2A, respectively) aligned to the reference genome while only 0.5% of the sequences from the third library (NEU8B) aligned to the reference genome. The high percentage of mapped sequences in the petrous sample is consistent with previous reports on the superior DNA preservation in this part of the skeleton (Gamba et al., 2014). Following the second sequencing run, the total endogenous content of the first two libraries was estimated to 39% and 28%, with an average sequence length of 39 bp and 38 bp, respectively (Figure 2—source data 1). The average depth of coverage was 0.65-fold for NEPEC and 0.14-fold for NEU2A. Both of them showed low frequencies of C to T substitutions at the 5’ and 3’ end, which are characteristic for Afu UDG-treated single-stranded libraries (Gansauge and Meyer, 2013), except for in CpG context, where deamination of 5-methylcytosine leaves thymine and not uracil (Figure 2—figure supplement 3).
We also processed sequencing data from an African forest elephant (SL0001) that was sequenced to high-coverage at the Broad Institute and re-processed sequencing data of an Asian elephant from (Lynch et al., 2015). We trimmed adapters with SeqPrep 1.1 using default parameters and aligned paired-end reads to loxAfr4 using BWA’s ‘aln’ algorithm and SAMtools’ ‘sampe’ command. Duplicate reads were removed with SAMtools’s ‘rmdup’. Moreover, we used the high-coverage genome of a woolly mammoth (Wrangel) from (Palkopoulou et al., 2015). The woolly mammoth alignments were re-processed for removal of duplicate reads with the custom script mentioned above.
To determine the phylogenetic relationships between the two P. antiquus specimens and other members of the Elephantidae family, we called pseudo-haploid consensus sequences for all autosomes of the two P. antiquus samples (~770 and~210 million sites, respectively). Sites with base quality below 30 and reads with mapping quality below 30 were filtered out. To exclude post-mortem damage-derived C to T substitutions, we trimmed 2 bp from the ends of all reads. We included regions of the loxAfr4 genome for which at least 90% of all possible 35-mers do not find a match at another position allowing for up to one mismatch, similar to the mappability filter described in (Prüfer et al., 2014). We used a majority-allele calling rule that required at least one read aligned at each position of the genome. Using the same approach, we called sequences for an Asian elephant (Uno [Lynch et al., 2015]), a woolly mammoth (Wrangel [Palkopoulou et al., 2015]) and an African forest elephant (SL0001; Broad Institute). We also used the reference sequence loxAfr4, as an African savanna elephant. We estimated the number of differences per base-pair for pairwise comparisons of all sequences and constructed a distance matrix, from which we built a Neighbor-joining (NJ) tree using PHYLIP version 3.696 (Felsenstein, 2005). To obtain support values for the nodes of the tree, we performed a bootstrap analysis (100 replicates) by splitting all autosomes in blocks of 5 Mb and randomly sampling blocks with replacement and built a majority-rule consensus tree.
We also extracted coding DNA sites (CDS) of protein-coding genes using the Ensembl 87 release for the loxAfr3 genome (downloaded from http://www.ensembl.org/) from each elephant genome sequence. CDS mapping to unknown chromosomes as well as CDS containing partial codons were excluded, resulting in a total of 86,212 CDS. Multiple sequence alignments were generated for each gene using MAFFT –ginsi (Katoh et al., 2002; Katoh and Standley, 2013) with 1000 iterations, which were concatenated into a single fasta file. Maximum likelihood phylogenetic analysis was performed with RAxML v8.2 (Stamatakis, 2014) with the GTRGAMMA model of nucleotide substitutions and 100 bootstrap trees. The resulting phylogeny is identical in its topology to that of the NJ tree with 100% bootstrap support (Figure 2—figure supplement 4).
Amino acid racemization (AAR) analyses were undertaken on the intra-crystalline protein from four individual Bithynia tentaculata opercula from the Eemian type-site, Amersfoort (Cleveringa et al., 2000): Amersfoort-1, upper depth 27.71, lower depth 28.50 (NEaar 2982–3, 3972 and 4681) and compared with previously published data from a single horizon at Neumark-Nord 1 (15.5.87/2, Schluffmudde, 25 cm under Anmoor = surface of the lower shore area; NEaar 5698–5703 [Penkman, 2010]) and several horizons from Neumark-Nord 2 (Sier et al., 2011). All samples were prepared using procedures of isolating the intra-crystalline protein by bleaching (Penkman et al., 2008). Two subsamples were then taken from each shell; one fraction was directly demineralized and the free amino acids analyzed (referred to as the 'free' amino acids, FAA, F), and the second was treated to release the peptide-bound amino acids, thus yielding the 'total' amino acid concentration, referred to as the ‘total hydrolysable amino acid fraction (THAA, H*). Samples were analyzed in duplicate by RP-HPLC. During preparative hydrolysis, both asparagine and glutamine undergo rapid irreversible deamination to aspartic acid and glutamic acid, respectively (Hill, 1965). It is therefore not possible to distinguish between the acidic amino acids and their derivatives and they are reported together as Asx and Glx, respectively. The D/L values of aspartic acid/asparagine, glutamic acid/glutamine, alanine and valine (D/L Asx, Glx, Ala, Val) are then assessed to provide an overall estimate of intra-crystalline protein decomposition (Penkman et al., 2011).
The Eemian stratotype locality at Amersfoort in the central Netherlands: a re-evaluation of old and new dataNetherlands Journal of Geosciences 79:197–216.https://doi.org/10.1017/S0016774600023659
A case study of apparent conflict between molecular phylogenies: the interrelationships of African elephantsCladistics : The International Journal of the Willi Hennig Society 21:31–50.https://doi.org/10.1111/j.1096-0031.2004.00044.x
Bayesian phylogenetics with BEAUti and the BEAST 1.7Molecular Biology and Evolution 29:1969–1973.https://doi.org/10.1093/molbev/mss075
PHYLIP (Phylogeny Inference Package), version Version 3.6Distributed by the Author. Department of Genome Sciences, University of Washington, Seattle.
Genome flux and stasis in a five millennium transect of European prehistoryNature Communications 5:5257.https://doi.org/10.1038/ncomms6257
Improving Bayesian population dynamics inference: a coalescent-based model for multiple lociMolecular Biology and Evolution 30:713–724.https://doi.org/10.1093/molbev/mss265
Triangulating the provenance of African elephants using mitochondrial DNAEvolutionary Applications 6:253–265.https://doi.org/10.1111/j.1752-4571.2012.00286.x
Das Pleistozän von Weimar-Ehringsdorf, Teil 2Abhandl. Zentr. Geol. Inst 23:1–596.
MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transformNucleic Acids Research 30:3059–3066.https://doi.org/10.1093/nar/gkf436
MAFFT multiple sequence alignment software version 7: improvements in performance and usabilityMolecular Biology and Evolution 30:772–780.https://doi.org/10.1093/molbev/mst010
Dating the arrival of straight-tusked elephant (Palaeoloxodon spp.) in EurasiaBulletin Du Musée d'anthropologie Préhistorique De Monaco. Supplément N° 6:123–128.
A new technique for precise uranium-series dating of travertine micro-samplesGeochimica Et Cosmochimica Acta 66:4261–4272.https://doi.org/10.1016/S0016-7037(02)00999-7
The positioning of the warm period of Neumark-Nord and its elephant-fauna in the context of the history of the earthIn: H Meller, editors. Elefantenreich – Eine Fossilwelt in Europa. Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt – Landesmuseum für Vorgeschichte Halle. pp. 65–69.
Illumina sequencing library preparation for highly multiplexed target capture and sequencingCold Spring Harbor Protocols 2010:pdb.prot5448.https://doi.org/10.1101/pdb.prot5448
Stegodontoidea, ElephantoideaProboscidea: A Monograph of the Discovery, Evolution, Migration and Extinction of the Mastodonts and Elephants of the World, Stegodontoidea, Elephantoidea, II, New York, The American Museum Press.
Neumark-Nord 1: Preliminary results of the amino acid analysisIn: H Meller, editors. Elefantenreich – Eine Fossilwelt in Europa. Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt – Landesmuseum für Vorgeschichte Halle. pp. 75–78.
leeHom: adaptor trimming and merging for Illumina sequencing readsNucleic Acids Research 42:e141.https://doi.org/10.1093/nar/gku699
Cytonuclear genomic dissociation in African elephant speciesNature Genetics 37:96–100.https://doi.org/10.1038/ng1485
Homo erectus in Africa, Pleistocene Evidence From the Middle AwashW Henry, W. H Gilbert, B Asfaw, editors. Univ. of California Press.
Cenozoic Mammals of AfricaL Werdelin, W. J Sanders, editors. Univ. of California Press.
ESR dating of a new Palaeolithic find layer of the travertine site of Weimar-Ehringsdorf (Central Germany)Terra Nostra 2:233–235.
Elefantenreich – Eine Fossilwelt in Europa71–74, ESR dating of tooth enamel samples from the archaeological find horizons of Neumark-Nord, Elefantenreich – Eine Fossilwelt in Europa, Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt – Landesmuseum für Vorgeschichte Halle.
Relationships within the Elephantinae using hyoid charactersQuaternary International 169-170:174–185.https://doi.org/10.1016/j.quaint.2007.02.003
New phylogenetic analysis of the family Elephantidae based on cranial-dental morphologyThe Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology 293:74–90.https://doi.org/10.1002/ar.21010
Vegetation, climate and radiocarbon datings in the Late Pleistocene of the Netherlands. Part I: Eemian and Early WeichselianMededelingen Van De Geologische Stichting, Nieuwe Serie 14:15–45.
Antonis RokasReviewing Editor; Vanderbilt University, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Palaeogenomes of Eurasian straight-tusked elephants challenge the current view of elephant evolution" for consideration by eLife. Your article has been favorably evaluated by Diethard Tautz (Senior Editor) and three reviewers, one of whom is a member of our Board of Reviewing Editors. The following individuals involved in review of your submission have agreed to reveal their identity: William Sanders (Reviewer #2); Thomas Gilbert (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The compact manuscript by Meyer and colleagues reports the sequencing of four full mitochondrial genomes and two partial nuclear genomes from Palaeoloxodon antiquus fossils. P. antiquus is a species that is thought to have been quite widespread throughout Eurasia in the Pleistocene. Phylogenetic analyses of these mitochondrial and nuclear data yield the surprising conclusion that the closest extant relatives of P. antiquus are not the Asian elephants (genus Elephas), as suggested by morphology, but the African elephants (genus Loxodonta). If the view offered by the genetic data holds, the current interpretation of the evolution of elephant morphology would need to be substantially revised. In addition, the fossils from which ancient DNA was extracted are quite old (~120-244 thousand years ago), so it is remarkable that such old samples yielded so much data.
1) It is not true that no one has suggested that relatives of Loxodonta made it out of Africa; Kretzoi did so emphatically in 1950, and others (e.g., van den Bergh; Markov and Saegusa) have accepted the possibility without endorsing the notion. These studies should be discussed and cited in the revised manuscript.
2) The authors do not treat the paleontological/morphological hypothesis about the relationships of P. antiquus/Loxodonta/Elephas as a competing hypothesis – they dismiss it in a single sentence, saying morphological features showing P. antiquus to be closely related to Elephas (such as E. recki) are convergences. And then go on to act as if elephant systematics were in a shamble and need serious revision – this simply isn't the case. They should present their view as a competing hypothesis and say if further work supports their hypothesis, then it has serious implications for elephant systematics and here is what those implications are. But, morphology isn't as labile as the authors suggest, and they do not provide compelling selective reasons for why a forest elephant such as P. antiquus would converge morphologically on a savanna grazer such as E. recki.
For example, one good test would be to study the abundant postcrania of Loxodonta, P. antiquus, and E. recki – there are clear differences in many aspects of the postcranium, especially in manus and pes elements, between Loxodonta and Elephas – see where Palaeoloxodon falls out. Additionally, it will be essential for the authors to discuss how additional species that were not sampled could be critical in resolving the proposed conundrum between the morphological and genetic results. Examples of such species include Palaeoloxodon namadicus, a straight-tusked species from Asia, Elephas planifrons and E. hysudricus (and other fossil species of Elephas, such as E. playcephalus) or any fossil species of Loxodonta, such as the abundant material from Laetoli of L. exoptata. Finally, it would be interesting to know if the authors tried but failed to sample any of the Elephas ekorensis – E. recki – E. iolensis lineage.
3) Execution of the phylogenetic analyses. The authors used two different methods on the mtDNA (Bayesian inference) and nuclear (neighbor-joining) data. It is not clear to me why that was the case. Bayesian inference is quite well established. In contrast, neighbor-joining is great for generating quick-and-dirty trees during the data exploration phase of experiments, but not very good in exploring tree space (the model of sequence evolution employed is also not described). I would recommend the authors also use BEAST to analyze the nuclear data (and at the same time estimate divergence times on that data set as well). If the use of BEAST is problematic due to computational constraints (if that's the case, the authors should include a brief description of those constraints in the manuscript), then my suggestion would be to use maximum likelihood to infer the nuclear phylogeny (e.g., use RAxML with the PROTGAMMAAUTOF option) and then r8s to calculate the divergence times.
4) Interpretation of the phylogenetic results, which are shown in Figure 2. The branch lengths of the African forest elephant – Palaeoloxodon clade are rather different between the mitochondrial and the nuclear phylogenetic trees. For example, the three NN samples are nearly identical in the mtDNA phylogeny but quite divergent in the nuclear one even though in general the nuclear genes appear to be somewhat slower evolving than the mitochondrial ones. In contrast, the key internode supporting the monophyly of the African forest elephant-Palaeoloxodon clade is quite long in the mtDNA phylogeny but very short in the nuclear phylogeny. Please discuss these two "discrepancies" and provide potential explanations.https://doi.org/10.7554/eLife.25413.028
- Matthias Meyer
- Mathias Stiller
- Sarah Nagel
- Beth Shapiro
- Yasuko Ishida
- Alfred L. Roca
- Kirsty E. H. Penkman
- Kirsty E. H. Penkman
- Michael Hofreiter
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank the Elephant Genome Sequencing Consortium, especially Elinor Karlsson and Jessica Alfoldi, for permission to use the African forest elephant genome data. We thank Chloé Piot, Ronny Barr and Susanne Sicker for help with preparing the figures, Adam Brandt and Tolulope Perrin-Stowe for sequence information and Adrian Lister for comments on the manuscript. Mitochondrial genome sequences were deposited in GenBank under accession nos. KY499555 - KY499558, KY616976, KY616978 and KY616979, nuclear sequence alignments are available through the European Nucleotide Archive (ENA) project no. PRJEB18563.
MM is supported by the Max Planck Society. BS was supported by the Gordon and Betty Moore Foundation. AR and YI were supported by the US Fish and Wildlife Service. The amino acid analyses were supported by Wellcome Trust and Leverhulme Trust funding to KP, with thanks to Sheila Taylor for technical support. D.R. is an investigator of the Howard Hughes Medical Institute. MH is supported by ERC consolidator grant 310763 GeneFlow.
- Antonis Rokas, Reviewing Editor, Vanderbilt University, United States
- Received: January 25, 2017
- Accepted: May 6, 2017
- Version of Record published: June 6, 2017 (version 1)
© 2017, Meyer et al
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.