Skip to main content

Molecular decay of enamel matrix protein genes in turtles and other edentulous amniotes

Abstract

Background

Secondary edentulism (toothlessness) has evolved on multiple occasions in amniotes including several mammalian lineages (pangolins, anteaters, baleen whales), birds, and turtles. All edentulous amniote clades have evolved from ancestors with enamel-capped teeth. Previous studies have documented the molecular decay of tooth-specific genes in edentulous mammals, all of which lost their teeth in the Cenozoic, and birds, which lost their teeth in the Cretaceous. By contrast with mammals and birds, tooth loss in turtles occurred in the Jurassic (201.6-145.5 Ma), providing an extended time window for tooth gene degradation in this clade. The release of the painted turtle and Chinese softshell turtle genomes provides an opportunity to recover the decayed remains of tooth-specific genes in Testudines.

Results

We queried available genomes of Testudines (Chrysemys picta [painted turtle], Pelodiscus sinensis [Chinese softshell turtle]), Aves (Anas platyrhynchos [duck], Gallus gallus [chicken], Meleagris gallopavo [turkey], Melopsittacus undulatus [budgerigar], Taeniopygia guttata [zebra finch]), and enamelless mammals (Orycteropus afer [aardvark], Choloepus hoffmanni [Hoffmann’s two-toed sloth], Dasypus novemcinctus [nine-banded armadillo]) for remnants of three enamel matrix protein (EMP) genes with putative enamel-specific functions. Remnants of the AMBN and ENAM genes were recovered in Chrysemys and retain their original synteny. Remnants of AMEL were recovered in both testudines, although there are no shared frameshifts. We also show that there are inactivated copies of AMBN, AMEL and ENAM in representatives of divergent avian lineages including Galloanserae, Passeriformes, and Psittaciformes, and that there are shared frameshift mutations in all three genes that predate the basal split in Neognathae. Among enamelless mammals, all three EMP genes exhibit inactivating mutations in Orycteropus and Choloepus.

Conclusions

Our results highlight the power of combining fossil and genomic evidence to decipher macroevolutionary transitions and characterize the functional range of different loci involved in tooth development. The fossil record and phylogenetics combine to predict the occurrence of molecular fossils of tooth-specific genes in the genomes of edentulous amniotes, and in every case these molecular fossils have been discovered. The widespread occurrence of EMP pseudogenes in turtles, birds, and edentulous/enamelless mammals also provides compelling evidence that in amniotes, the only unique, non-redundant function of these genes is in enamel formation.

Background

Gnathostomes are unique among vertebrates in utilizing jaws and teeth for food acquisition and processing. Divergent gnathostome lineages exhibit a wide array of dental modifications in association with diverse dietary specializations. Multiple iterations of secondary tooth loss have also evolved in gnathostomes, perhaps most famously in birds, turtles, and several mammalian lineages including baleen whales, pangolins, and anteaters [13]. There are also mammalian species with enamelless teeth, i.e., sloths, armadillos, aardvarks, pygmy and dwarf sperm whales [2].

Tooth development is an intricate process that encompasses a complex series of epithelial-mesenchymal interactions involving growth factors, transcription factors, and signal receptors that affect tooth shape, tooth number, and cusp number [49]. Several hundred genes are associated with tooth development [6], including members of the hedgehog (Hh), fibroblast growth factor (Fgf) and bone morphogenic protein (Bmp) families that mediate epithelio-mesenchymal signaling interactions [10]. Most or all of the genes that are involved in early tooth development are thought to be pleiotropic and have additional functions outside of tooth formation. However, there are also genes with putative tooth-specific functions that are expressed later in development by dentin-forming odontoblasts and/or enamel-forming ameloblasts [1117]. These genes affect the physical properties of teeth including enamel thickness and structure. The durability of teeth, in combination with tooth-specific genes that impact the physical structure of dentin and enamel, make teeth a model system for studying the coevolution of morphological change in the fossil record and molecular change in the genome.

Previous studies have shown that the genes encoding three enamel matrix proteins (EMPs), enamelin, amelogenin, and ameloblastin, have become pseudogenized in one or more edentulous/enamelless mammals by frameshift mutations and/or stop codons [13] (Table 1). Further, the molecular decay of tooth-specific genes in placental mammal lineages has been shown to parallel the morphological degeneration of enamel in the fossil record [13]. The retention of inactivated EMP genes in mammalian lineages is not surprising given that all edentulous and enamelless taxa appear to have originated in the Cenozoic [2].

Table 1 Pseudogene remnants of enamel matrix protein (EMP) genes in edentulous and enamelless amniotes

By contrast, tooth loss in both birds and turtles occurred in the Mesozoic, providing a longer time period for the molecular decay of tooth-specific genes in these lineages. In birds, the presence of a horny beak and gizzard have presumably compensated for edentulism in food acquisition and processing [22]. Edentulism in the ancestry of modern birds occurred in the Cretaceous, between 125 and 65.5 Ma, and also evolved independently in several lineages of extinct birds that are on the stem to Neornithes (crown group birds) [22]. Even though tooth loss in Neornithes occurred in the Cretaceous, the chicken (Gallus gallus) genome retains pseudogenized copies of both AMEL and ENAM[18, 20, 21].

Edentulism in Testudines (turtles) occurred even earlier, at least as far back as the Late Jurassic [23, 24]. The oldest turtle is Odontochelys semitestacea from the Triassic of China (~220 Ma) [25]. Marginal and palatal teeth are both present in this taxon [25]. Slightly younger (~210 Ma) is Proganochelys quenstedti from the Upper Triassic [23, 26]. Marginal teeth are absent in Proganochelys, but palatal teeth are present. Phylogenetic analyses suggest that Proganochelys, like Odontochelys, is a stem testudine [23]. The oldest testudines with unambiguous crown-group affinities, including tooth loss on the palatines, vomer, and pterygoids, are from the Late Jurassic [23].

To date, remnants of pseudogenized copies of tooth-specific genes in Testudines have not been reported. Girondot and Sire [27] were unable to amplify fragments of amelogenin in turtles with degenerate PCR primers, and it remains unclear if vestiges of tooth-specific genes are retained in testudine genomes given the antiquity of tooth loss in this clade. On the other hand, phylogenomic data suggest that Testudines have slower rates of nuclear gene evolution than birds and mammals (fig. 1 in [28]). The recent release of two testudine genomes, Chrysemys picta (painted turtle) and Pelodiscus sinensis (Chinese softshell turtle), provides an opportunity to screen for remnants of tooth-specific genes in turtles. These taxa, both of which belong to Cryptodira (hidden neck turtles), index the basal cladogenic event among crown cryptodires and leave only Pleurodira (sideneck turtles) unrepresented among the three oldest crown-testudine lineages. Cryptodira and Pleurodira have traditionally been regarded as reciprocally monophyletic based on morphology [2931] and large molecular data sets [32, 33]. An exception is Barley et al.’s [32] coalescence analysis, which recovered a basal split between softshells plus the pig-nosed turtle (Carettochelys insculpta) and Pleurodira plus other Cryptodira. However, Barley et al.’s [32] concatenation analysis recovered Cryptodira and Pleurodira. Moreover, the most complete testudine phylogeny with a variety of outgroups [33] supports the monophyly of both Cryptodira and Pleurodira.

Here, we report the results of querying the genomes of Chrysemys and Pelodiscus for inactivated remnants of three EMP genes, ENAM, AMEL, and AMBN, all of which are hypothesized to have tooth-specific functions. We also report the results of querying genome sequences of five birds (Gallus gallus [chicken], Meleagris gallopavo [turkey], Taeniopygia guttata [zebra finch], Anas platyrhynchos [duck], Melopsittacus undulatus [budgerigar]) and three enamelless mammals (Orycteropus afer [aardvark], Choloepus hoffmanni [Hoffmann’s two-toed sloth], Dasypus novemcinctus [nine-banded armadillo]) for remnants of these EMP genes. The occurrence of pseudogenized remnants of EMP genes in highly divergent amniotes, including turtles, birds, and edentulous/enamelless mammals, would provide robust evidence for the hypothesis that the only essential, non-redundant function of AMEL, AMBN, and ENAM is in enamel formation. Moreover, the search for shared frameshift mutations in EMP genes provides an opportunity to date the timing of enamel loss in these clades.

Methods

The genomes of Chrysemys picta (painted turtle; PreEnsembl) and Pelodiscus sinensis (Chinese softshell turtle; Ensembl 68) were queried with BLASTN using crocodylian mRNA sequences for three enamel matrix proteins genes: Paleosuchus palpebrosus (Cuvier’s dwarf caiman) AMEL (AF095568), Caiman crocodilus (spectacled caiman) AMBN (AY043290), and Crocodylus niloticus (Nile crocodile) ENAM (GU344683). Intron-exon boundaries of mRNA sequences were determined by blasting crocodylian mRNAs against the recently released genome sequence of Alligator mississippiensis (American alligator). Crocodylian sequences were chosen as probes based on recent evidence that provides strong support for a sister group relationship between Testudines and Archosauria (crocodylians, birds) [27, 34]. AMBN has 11 exons in Alligator, all of which include protein-coding regions. AMEL includes six exons in Alligator, five of which (1–5) contain protein-coding regions. Finally, ENAM includes nine exons in Alligator, eight of which (2–9) include protein-coding regions. Remnants of EMP genes that were recovered from Chrysemys were subsequently used to query the Pelodiscus genome, and remnants of EMP genes from Pelodiscus were used to query the Chrysemys genome. Five bird genomes (Anas platyrhynchos [duck], Gallus gallus [chicken], Meleagris gallopavo [turkey], Melopsittacus undulatus [budgerigar], Taeniopygia guttata [zebra finch]) were queried with crocodylian mRNAs and virtual pseudogene mRNAs for AMEL (EU340348) and ENAM (GU198360) from Gallus, as well as with positive hits from the initial round of screening for the bird genomes. The NCBI genomes of three enamelless mammals (Orycteropus afer [aardvark], Choloepus hoffmanni [Hoffmann’s two-toed sloth], Dasypus novemcinctus [nine-banded armadillo]) were queried with gene sequences from mammals that retain enamel-capped teeth (i.e., Bos taurus, Homo sapiens, Echinops telfairi, Elephantulus edwardii, Trichechus manatus). Sequences were aligned with Se-Al [35].

Selection analyses on branches were performed with PAML 4.5 [36] to estimate dN/dS values for an assemblage of 20 complete or nearly complete mammalian AMBN sequences, all of which were derived from Ensembl 68 (Ailuropoda melanoleuca, Bos taurus, Callithrix jacchus, Equus caballus, Oryctolagus cuniculus) or NCBI (Canis lupus familiaris [XM_539304], Choloepus hoffmanni [ABVD01066982, TI 1336370099, ABVD01066984, TI 1338793682, TI 1368206053, ABVD01066985, ABVD01066986], Dasypus novemcinctus [AAGV03237582], Echinops telfairi [AAIY02097297, AAIY02097298], Elephantulus edwardii [AMGZ01205754, AMGZ01205755], Homo sapiens [NM_016519], Loxodonta africana [AAGU03054364, AAGU03054363], Orycteropus afer [ALYB01124783], Pan paniscus [XM_001160377], Pan troglodytes [XM_001160377], Pongo abelii [XM_002814835], Procavia capensis [ABRQ01419216, ABRQ01419217, ABRQ01419218, ABRQ01419219, ABRQ01419220], Sus scrofa [NM_214037]). Sequences were aligned with Se-Al [35]. Branch analyses were performed with four branch categories: Orycteropus, Choloepus, Dasypus, and all other branches. Orycteropus, Choloepus and Dasypus were recognized as separate branches in the PAML runs based on the hypothesized independent loss of enamel in these lineages [2]. PAML analyses were conducted with codon frequency models 1, 2, and 3. We used a composite species tree based on Meredith et al. [37] and Springer et al. [38]. Approximate divergence dates and fossil ages in Figure 1 are taken from the following sources: Amniota and Testudines to Archosauria [39]; Pelodiscus to Chrysemys[40]; Odontochelys[25]; Neognathae [41]; Galloanserae and Gallus to Meleagris[42]; Taeniopygia to Melopsittacus[43]; Ichthyornis[44]; Xenarthra [45]; Xenarthra to Orycteropus[37]; Eomaia[46].

Results

ENAM Gene

Blast searches with Crocodylus niloticus (Nile crocodile) ENAM (GU344683) identified almost the entirety of exon 9 in Chrysemys picta (Additional file 1). Remnants of this exon in C. picta were identified on scaffold JH584398, which has not yet been mapped onto a chromosome. The ENAM sequence for exon 9 in Chrysemys includes numerous frameshift indels (Figure 1, Table 2, Additional files 1, 2), one of which is a SINE insertion that shares homology with members of the Cry family of SINEs in Cryptodira [47] (Additional files 1, 2). An additional SINE insertion occurs 3’ to the stop codon. Remnants of exon 9 in Chrysemys comprise ~3840 bp after excluding the SINE insertions. Chrysemys exon 9 retains a “TAA” stop codon in the same location as Crocodylus. By contrast with the results of blast searches against the Chrysemys genome, blast searches with Crocodylus ENAM resulted in no significant hits to the Pelodiscus sinensis genome. Likewise, blast searches with the recovered ENAM pseudogene segment from Chrysemys did not find significant matches in the Pelodiscus genome.

Figure 1
figure 1

Timetree of edentulous and enamelless amniote taxa for which genomic sequences are available (Ensembl, PreEnsembl, NCBI) with a mapping of inferred inactivating mutations (frameshifts, premature stop codons, splice site mutations) in three EMP genes to particular branches in the tree. Branches with inactivating mutations are marked by vertical bars (pink = ENAM [E], blue = AMBN [B], green = AMEL [A]). Extinct stem taxa with teeth (Odontochelys, Ichthyornis, Eomaia) are also shown. Mutations were mapped to branches by Fitch parsimony with delayed transformation optimization. Teeth with enamel, teeth without enamel, and edentulism are denoted with black, gray, and white-filled circles, respectively, at terminal and internal nodes. In combination with fossil evidence, dN/dS ratios and frameshift mutations in ENAM suggest that enamel was lost independently in the common ancestor of Pilosa (sloths and anteaters) and in multiple armadillo lineages, including Dasypus[2]. Remnants of AMBN and ENAM were not recovered from the Pelodiscus genome. References for divergence dates and fossil ages are provided in Methods.

Table 2 Examples of inactivating mutations in enamel matrix protein genes among edentulous and enamelless amniotes with genome sequences

Blast searches with the virtual pseudogene mRNA of Gallus ENAM exon 9 recovered homologous fragments of this exon in Meleagris gallopavo (turkey), Taeniopygia guttata (zebra finch), Anas platyrhynchos (duck), and Melopsittacus undulatus (budgerigar). The 3’ region of M. gallopavo ENAM was recovered on chromosome 1, whereas the 5’ region was recovered on the Z chromosome of this species. All of the exon 9 sequences were characterized by frameshift mutations (Figure 1), including a 406 bp deletion near the 3’ end of the coding sequence of exon 9 that is shared by all five birds (Additional file 1).

ENAM sequences for Orycteropus afer, Choloepus hoffmanni and Dasypus novemcinctus were reported by Meredith et al. [2] and are not duplicated here. Orycteropus ENAM includes three frameshift mutations [2]. Among xenarthrans, Choloepus ENAM includes numerous frameshift mutations whereas Dasypus ENAM includes a single frameshift mutation that is located near the carboxy-terminal end of the coding sequence in exon 10 (= exon 9 of some non-mammalian vertebrates) [2].

AMBN Gene

Blast searches with Caiman crocodilus (spectacled caiman) AMBN (AY043290) recovered portions of exons 1, 2, 3, 6, 8, 9, 10, and 11, in this order, on Chrysemys scaffold JH584398 (Additional files 3 and 4). The AMBN exons are located immediately upstream of ENAM exon 9 on scaffold JH584398. Frameshift mutations occur in exons 6, 10, and 11 (Figure 1, Table 2, Additional files 3, 4). AMBN sequences were not found in Pelodiscus.

Among birds, remnants of AMBN were discovered in Anas (scaffold 247: exons 2, 3, 5, 6, 7, 8), Gallus (Z chromosome: exons 2, 3, 5, 6), Taeniopygia (Z chromosome; exons 6, 7, 8), and Melopsittacus (scaffold JH556633: exons 7, 8). Frameshift mutations occur in Anas (exons 2, 6, 8), Gallus (exon 3) and Taeniopygia (exon 6) (Additional file 3). A ten-bp frameshift deletion is shared by Anas, Gallus, and Taeniopygia, which together index the common ancestry of Neognathae, although this region was not recovered in Meleagris and Melopsittacus. There were no frameshifts in Melopsittacus exons 7 and 8, but a stop codon occurs in the latter exon.

Among mammals, we recovered exons 1–5, 7, and 10–13 in Orycteropus afer (Additional file 5). Exon 5 contains a stop codon and there is a 13-bp frameshift deletion in exon 13 (Table 2). Exon 6 has been deleted along with portions of introns 5 and 6 (Table 2). Exons 8 and 9 are 39-bp duplications of exon 7 that occur in primates [48] and their absence in Orycteropus is not unexpected. With the exception of Elephantulus edwardii, which has two duplications of exon 7 as in Homo, other afrotherians (Loxodonta africana, Trichechus manatus, Procavia capensis, Chrysochloris asiatica, Echinops telfairi) are similar to Orycteropus in possessing only exon 7. We also recovered sequences for exons 1–7 and 10–13 in Dasypus AMBN (Additional file 5). The recovered protein-coding sequence is intact, although the presumed start codon is ten codons downstream of the start codon in the human sequence and occurs in exon 2 rather than exon 1, as is also the case for Loxodonta africana (African elephant) (Additional file 6). Complete coding sequences for nine exons (1–5, 7, 10–12) and partial sequences for exons 6 and 13 were recovered from the Choloepus genome (Additional file 5). The putative start codon in Choloepus is located in exon 2 as in Dasypus. There are no frameshift mutations or stop codons in the available Choloepus sequence, but intron 2 exhibits an inactivating mutation at the acceptor splice site (AG to AT; see TI# 1338556193 in Trace Archives) (Table 2). By contrast, all of the donor and accepter splice sites exhibit canonical GT and AG motifs, respectively, in 16 mammals with enamel-capped teeth (taxon names provided in Methods).

Selection analyses based on an alignment of 20 mammalian AMBN sequences (Additional file 6) with codon frequency (CF) models 1, 2, and 3 all suggest elevated dN/dS ratios on the Orycteropus (CF1 = 0.71, CF2 = 0.77, CF3 = 0.89), Choloepus (CF1 = 0.86, CF2 = 0.90, CF3 = 1.04) and Dasypus (CF1 = 0.64, CF2 = 0.71, CF3 = 0.74) branches relative to the background rate (CF1 = 0.48, CF2 = 0.51, CF3 = 0.56) in mammalian taxa with enamel-covered teeth.

AMEL Gene

Blast searches with Paleosuchus palpebrosus (Cuvier’s dwarf caiman) AMEL (AF095568) recovered most of exons 1, 3, 4, and 5 in Pelodiscus, with ancestral synteny preserved, on scaffold JH208023 (Additional files 7 and 8). Exons 1 and 5 were identified on scaffold JH584884 in Chrysemys (Additional files 7 and 8). Chrysemys and Pelodiscus both retain the same start codon as in crocodylians. Frameshift mutations are present in sequences of both testudine species (exon 1 in Chrysemys, exons 1, 3, and 4 in Pelodiscus), although there are no shared frameshifts in these two exons (Figure 1, Table 2, Additional files 7 and 8). In addition, Pelodiscus and Chrysemys share exon-flanking sequences that are 3’ to exon 1 and 5’ to exon 5. The sequence that flanks exon 1 is > 1000 bp and includes regions that are homologous to intron 1 and possibly exon 2 of Alligator mississippiensis (Additional file 6). Intron 1 retains the canonical “GT” splice donor site in both testudines (Additional file 9). The homology of exon 2 is less certain, but nevertheless of potential importance as there are multiple shared frameshift mutations in Chrysemys and Pelodiscus (Additional file 9). The sequence that flanks exon 5 is > 600 bp long and is homologous to a portion of intron 4 of A. mississippiensis with the canonical “AG” splice acceptor site in Chrysemys (“AA” in Pelodiscus) (Additional file 10).

In birds, blast searches with Gallus AMEL pseudogene mRNA (EU340348) recovered homologous DNA sequences in Anas (exons 1–4 with original synteny on scaffold 415), Meleagris (exons 1–4 with original synteny on chromosome 1), Melopsittacus (exons 1–4 with original synteny on AGAI01063440), and Taeniopygia (exons 2–4 with original synteny on chromosome 1) (Additional file 7). All of the birds except for Anas share a common start codon with crocodylians. Frameshift mutations occur in all avian taxa, including a one-bp deletion in exon 4 that is shared by all five birds (Figure 1).

Among mammals, there are two stop codons in exon 6 of Orycteropus AMEL (Figure 1, Table 2, Additional file 11). Choloepus AMEL shows replacement of the methionine initiation codon in exon 2 by a threonine codon, an 11 bp frameshift deletion in exon 2, and a premature stop codon followed by a frameshift deletion in exon 6 (Figure 1, Table 2, Additional file 11). All of the inactivating mutations in Choloepus are corroborated by chromatograms in NCBI’s Trace Archives. The sequence for Dasypus is intact through exon 6 (excepting exon 4, see below), but a deletion incorporates the 3’ end of intron 6, including the acceptor splice site, and the first 14 bp of exon 7, including the stop codon. Orycteropus, Choloepus, and Dasypus also lack a functional copy of exon 4, either because of stop codons (Orycteropus, Choloepus) or a donor splice site mutation (GT to AT) in intron 4 that is adjacent to this exon. However, a functional copy of exon 4 is variably present in Mammalia [48] and is missing from six other afrotherians (Chrysochloris, Echinops, Elephantulus, Loxodonta, Procavia, Trichechus) with genome sequences because of stop codons, frameshift mutations, and/or splice site mutations (data not shown).

Discussion

EMP Genes in Turtles

Crown Testudines comprise 319 extant species [49], all of which are toothless and instead have a keratinized beak that in combination with strong jaw muscles allows testudines to tear food and capture prey [21]. The oldest testudine fossil is the Late Triassic Odontochelys semitestacea from China [25]. Odontochelys lacked a beak and instead retained both marginal and palatal teeth [25]. Proganochelys quenstedti is also known from the Late Triassic, and is more derived than Odontochelys. Proganochelys lacked marginal teeth, but palatal teeth were present. The morphology of the mandibles suggests that Proganochelys possessed a keratinized beak [26]. It remains unclear if palatal teeth, which occur on the vomer, palatine, and pterygoid bones, were lost in the common ancestor of all living turtles or independently in cryptodires and pleurodires owing to alternate phylogenetic hypotheses wherein taxa with palatal teeth such as Kayentachelys apix are positioned as stem testudines or as stem cryptodires [23, 24, 2931, 50, 51]. In either case palatal teeth were lost no later than the Late Jurassic in the ancestry of crown testudines.

The retention of pseudogenized copies of three EMP genes (AMEL, AMBN, ENAM) in Chrysemys is perhaps surprising given the antiquity of tooth loss in turtles and the extended time window for large-scale deletions and/or rearrangements to erase or scramble the genomic instructions for enamel production in this edentulous taxon. However, remnants of AMBN and ENAM exon 9 in Chrysemys occur in juxtaposition to each other on the same contig. Kawasaki’s [52] reconstructions of SCPP gene order in both stem tetrapods and stem amniotes suggest that AMBN is immediately upstream of ENAM. Thus, it appears that remnants of this original gene order are present in the painted turtle genome. The protein-coding regions of AMBN and ENAM in Chrysemys have been battered by inactivating mutations, including SINE insertions in exon 9 of ENAM, but nevertheless retain unambiguous signatures of their heritage. The only EMP gene that was recovered in Pelodiscus is AMEL. It remains unclear whether AMBN and ENAM were completely deleted from the Pelodiscus genome during evolutionary history, are difficult to recognize because of numerous mutations, or are missing from the current assembly of the Pelodiscus genome due to incomplete sequencing assembly. Pelodiscus AMEL is relatively complete and includes exons 1, 3, 4, and 5 in their original syntenic order. By contrast, we only discovered exons 1 and 5 of Chrysemys AMEL. The two exons represent the only overlapping EMP fragments that were discovered in both Pelodiscus and Chrysemys. Frameshift mutations in AMEL occur in both species, but the only shared frameshift mutations occur in putative exon 2, which has equivocal homology with the same exon in crocodylians. Nevertheless, if these frameshift mutations are genuine then there is molecular evidence for the inactivation of enamel production in the common ancestry of Cryptodira.

The genetic and developmental basis of edentulism in turtles involves both tooth-specific and pleiotropic genes. Results presented here demonstrate that three tooth-specific EMP genes were pseudogenized in the ancestry of modern cryptodires. Previously, Tokita et al. [53] examined patterns of gene expression in the developing lower jaw of Pelodiscus sinensis and reported impairment of Shh signaling in the oral epithelium along with early-stage arrest of odontoblast development by abrogation of Msx2 expression in dental mesenchyme. It remains unclear if changes in gene expression or EMP pseudogenization occurred first in the evolutionary history of turtles.

EMP Genes in Birds

AMEL and ENAM pseudogenes have previously been reported in Gallus[20], which belongs to Galliformes. We provide extended evidence for inactivated EMP genes in birds and show that there are remnants of AMEL and ENAM in four other birds with complete genome sequences: Meleagris (Galliformes), Anas (Anseriformes), Melopsittacus (Psittaciformes), and Taeniopygia (Passeriformes) (Table 1). In addition, we recovered molecular evidence for an inactivated copy of AMBN in Anas, Gallus, Taeniopygia, and Melopsittacus. It is noteworthy that there are putative frameshift mutations in AMBN, AMEL, and ENAM that are shared by representatives of Galliformes (Gallus, Meleagris), Anseriformes (Anas), and Neoaves (Taeniopygia, Melopsittacus), which together index the deepest split in Neognathae [54]. It will be important to determine if these inactivating mutations are shared with Palaeognathae (tinamous and ratites), in which case molecular evidence would be consistent with the loss of enamel in the common ancestor of Neornithes rather than independently in Neognathae and Palaeognathae.

EMP Genes in Mammals

It has previously been shown that ENAM is a pseudogene in a wide array of edentulous and enamelless mammals [1, 2] (Table 1). Further, AMEL and AMBN sequences have been reported for baleen whales and in both cases are inactivated in multiple species [1, 3]. However, the functional versus pseudogene status of AMEL and AMBN remains to be investigated in most edentulous and enamelless taxa.

Here, we provide evidence of inactivating mutations in AMEL in Orycteropus afer (aardvark) and Choloepus hoffmanni (Hoffmann’s two-toed sloth), both of which have enamelless teeth comprised of dentin. Multiple frameshift mutations in Choloepus AMEL are confirmed by Trace Archives chromatograms (Additional file 8). The two stop codons in Orycteropus are based on Illumini Hi-Seq sequencing technology with 44X genome coverage.

Previous studies suggest that AMBN may play a role in dentin formation/regeneration [55, 56]. However, this hypothesis is contradicted by our finding that there are inactivating mutations in Orycteropus and Choloepus AMBN. Moreover, dN/dS ratios indicate relaxed purifying selection in these taxa. Similarly, McGowen [19] reported a premature stop codon in exon 13 (JF504758) of Monodon monoceros (narwhal). Narwhals have vestigial teeth, and in the case of males a single enlarged tusk. Both the vestigial teeth and tusk are composed of dentin and cementum, but no enamel [5759]. The occurrence of a premature stop codon in Monodon AMBN, in conjunction with inactivating mutations in two other lineages with enamelless teeth (Orycteropus, Choloepus), is consistent with the hypothesis that the only essential, non-redundant role of AMBN is in enamel formation. Along these lines, ameloblastin knockout mice do not form an enamel layer, and lack dentin defects [13, 15].

By contrast with Choloepus, where enamel loss occurred more than 50 million years ago, enamel degeneration in Dasypus occurred more recently [2]. Dasypus retains a vestigial enamel or enamel-like substance that covers the dentin of some teeth and is quickly worn off. Moreover, there are no frameshifts in ENAM that are shared by Dasypus and other armadillo genera [2]. The single frameshift in Dasypus ENAM precludes a functional, full-length enamelin protein, but the location of the frameshift is close to the 3’ end of the protein-coding sequence in exon 10, and shorter protein products are possible [2]. Our finding that AMBN is intact in Dasypus, and that only the last amino acid is truncated from the coding region of AMEL, is consistent with the presence of vestigial ‘enamel’ in this species [2]. We also note that a partial AMBN sequence (JF701624) for Kogia breviceps (pygmy sperm whale) has an intact protein-coding region even though K. breviceps has enamelless teeth and frameshift mutations in the ENAM gene. It will be important in future studies to sequence the remaining exons of Kogia AMBN to determine if the missing exonic regions contain frameshift mutations and/or other inactivating mutations.

Conclusions

Edentulism has evolved independently in multiple lineages of living amniotes including turtles, birds, echidnas, baleen whales, anteaters, and pangolins. There are also mammals with enamelless teeth including pygmy sperm whale, narwhal, sloths, armadillos, and aardvarks. In every case these edentulous or enamelless forms have descended from ancestors with enamel-capped teeth. Thus, amniote diversity provides a natural laboratory for testing hypotheses of tooth-specific gene function [3]. Moreover, this laboratory includes multiple, replicated experiments. AMBN, AMEL, and ENAM have all been postulated to have tooth-specific or even enamel-specific gene functions [1117], although pleiotropic functions have been suggested for AMBN[9, 48, 60, 61] and AMEL[62, 63]. The widespread occurrence of EMP pseudogenes in turtles, birds, and several mammalian lineages (Table 1) provides compelling evidence that the only unique, non-redundant function of these genes in amniotes is in enamel formation: functional copies of these genes have not been retained by natural selection when enamel production was abrogated independently in distantly related lineages. By contrast, representative sequences from amniotes with enamel-capped teeth retain intact coding sequences for AMBN[1, 48]; also see Additional file 11, AMEL[3, 48, 64], and ENAM[1, 2, 18, 37, 65].

The evolution of tooth loss in multiple amniote lineages also provides a model system for integrating the fossil record, phylogenetics, and genomics. This system allows for reciprocal hypothesis testing and provides a multifaceted, synthetic view on macroevolutionary transitions in testudines, birds, and edentulous/enamelless mammals [1, 2]. The fossil record and phylogenetics combine to predict the occurrence of molecular fossils of tooth-specific genes in the genomes of edentulous and enamelless amniotes, and molecular fossils have been discovered in all lineages that have been investigated. We suggest that these striking patterns may prove especially useful for educating the public on the convergence of evidence that clearly documents the evolutionary process over deep time.

References

  1. Deméré TA, McGowen MR, Berta A, Gatesy J: Morphological and molecular evidence for a stepwise evolutionary transition from teeth to baleen in mysticete whales. Syst Biol. 2008, 57: 15-37. 10.1080/10635150701884632.

    Article  PubMed  Google Scholar 

  2. Meredith RW, Gatesy J, Murphy WJ, Ryder OA, Springer MS: Molecular decay of the tooth gene enamelin (ENAM) mirrors the loss of enamel in the fossil record of placental mammals. PLoS Genet. 2009, 5: e1000634-10.1371/journal.pgen.1000634.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Meredith RW, Gatesy J, Cheng J, Springer MS: Pseudogenization of the tooth gene enamelysin (MMP20) in the common ancestor of extant baleen whales. Proc R Soc B. 2011, 278: 993-1002. 10.1098/rspb.2010.1280.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Scarel-Caminaga RM, Pasetto S, da Silva ER, Peres RCR: Genes and tooth development: reviewing the structure and function of some key players. Braz J Oral Sci. 2003, 2: 339-347.

    Google Scholar 

  5. Chai Y, Maxson RE: Recent advances in craniofacial morphogenesis. Dev Dynamics. 2006, 235: 2353-2375. 10.1002/dvdy.20833.

    Article  Google Scholar 

  6. Thesleff I: The genetic basis of tooth development and dental defects. Am J Med Gen. 2006, 140A: 2530-2535. 10.1002/ajmg.a.31360.

    Article  Google Scholar 

  7. Kapadia H, Mues G, D’Sousa R: Genes affecting tooth morphogenesis. Orthod Craniofacial Res. 2007, 10: 105-113. 10.1111/j.1601-6343.2007.00395.x.

    Article  CAS  Google Scholar 

  8. Bei M: Molecular genetics of tooth development. Curr Opin Gen Dev. 2009, 19: 504-510. 10.1016/j.gde.2009.09.002.

    Article  CAS  Google Scholar 

  9. Catón J, Tucker AS: Current knowledge of tooth development: patterning and mineralization of the murine dentition. J Anat. 2009, 214: 502-515. 10.1111/j.1469-7580.2008.01014.x.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Fleischmannova J, Matalova E, Tucker AS, Sharpe PT: Mouse models of tooth abnormalities. Eur J Oral Sci. 2008, 116: 1-10. 10.1111/j.1600-0722.2007.00504.x.

    Article  PubMed  Google Scholar 

  11. Snead ML, Lau EC, Zeichner-David M, Fincham AG, Woo SLC, Slavkin HC: DNA sequence for cloned cDNA for murine amelogenin reveal the amino acid sequence for enamel-specific protein. Biochem Biophys Res Comm. 1985, 129: 812-818. 10.1016/0006-291X(85)91964-3.

    Article  CAS  PubMed  Google Scholar 

  12. Rajpar MH, Harley K, Laing C, Davies RM, Dixon MJ: Mutation of the gene encoding the enamel-specific protein, enamelin, causes autosomal-dominant amelogenesis imperfecta. Hum Mol Gen. 2001, 10: 1673-1677. 10.1093/hmg/10.16.1673.

    Article  CAS  PubMed  Google Scholar 

  13. Fukumoto S, Kiba T, Hall B, Iehara N, Nakamura T, Longenecker G, Krebsbach PH, Nanci A, Kulkarni AB, Yamada Y: Ameloblastin is a cell adhesion molecule required for maintaining the differentiation state of ameloblasts. J Cell Biol. 2004, 167: 973-983. 10.1083/jcb.200409077.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Fukumoto S, Arakaki M, Iwamoto T, Yamada A, Miyamoto R, Naruse M, Nakamura T: Epithelial cell lines in the field of dental research: review. Interface of Oral Health Science 2011. Edited by: Sasaki K, Suzuki O, Takahashi N. 2012, New York: Springer, 327-333.

    Chapter  Google Scholar 

  15. Hu JC-C, Hu Y, Smith CE, McKee MD, Wright JT, Yamakoshi Y, Papagerakis P, Hunter GK, Feng JQ, Yamakoshi F, Simmer JP: Enamel defects and ameloblast-specific expression in Enam knock-out/lacZ knock-in mice. J Biol Chem. 2008, 283: 10858-10871. 10.1074/jbc.M710565200.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Hatakeyama J, Fukumoto S, Nakamura T, Haruyama N, Suzuki S, Hatakeyama Y, Shum L, Gibson CW, Yamada Y, Kulkarni AB: Synergistic roles of amelogenin and ameloblastin. J Dent Res. 2009, 88: 318-322. 10.1177/0022034509334749.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Chan H-C, Estrella NMRP, Milkovich RN, Kim J-W, Simmer JP, Hu JC-C: Target gene analyses of 39 amelogenesis imperfecta kindreds. Eur J Oral Sci. 2011, 119 (Suppl 1): 311-323.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Al-Hashimi N, Lafont A-G, Delgado S, Kawasaki K, Sire J-Y: The enamelin genes in lizard, crocodile, and frog and the pseudogene in the chicken provide new insights on enamelin evolution in tetrapods. Mol Biol Evol. 2010, 27: 2078-2094. 10.1093/molbev/msq098.

    Article  CAS  PubMed  Google Scholar 

  19. McGowen MR: Toward the resolution of an explosive radiation – a multilocus phylogeny of oceanic dolphins (Delphinidae). Mol Phylogenet Evol. 2011, 60: 345-357. 10.1016/j.ympev.2011.05.003.

    Article  PubMed  Google Scholar 

  20. Sire J-Y, Delgado SC, Girondot M: Hen’s teeth with enamel cap: from dream to impossibility. BMC Evol Biol. 2008, 8: 246-10.1186/1471-2148-8-246.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Davit-Béal T, Tucker AS, Sire J-Y: Loss of teeth and enamel in tetrapods: fossil record, genetic data and morphological adaptations. J Anat. 2009, 214: 477-501. 10.1111/j.1469-7580.2009.01060.x.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Louchart A, Viriot L: From snout to beak: the loss of teeth in birds. Trends Ecol Evol. 2011, 26: 663-673. 10.1016/j.tree.2011.09.004.

    Article  PubMed  Google Scholar 

  23. Joyce WG: Phylogenetic relationships of Mesozoic turtles. Bull Peabody Mus Nat Hist. 2007, 48: 3-102. 10.3374/0079-032X(2007)48[3:PROMT]2.0.CO;2.

    Article  Google Scholar 

  24. Anquetin J, Barrett PM, Jones MEH, Moore-Fay S, Evans SE: A new stem turtle from the Middle Jurassic of Scotland: new insights into the evolution and palaeoecology of basal turtles. Proc Roy Soc B. 2009, 276: 879-886. 10.1098/rspb.2008.1429.

    Article  Google Scholar 

  25. Li C, Wu X-C, Rieppel O, Wang L-T, Zhao L-J: An ancestral turtle from the Late Triassic of southwestern China. Nature. 2008, 456: 497-501. 10.1038/nature07533.

    Article  CAS  PubMed  Google Scholar 

  26. Gaffney ES: The comparative osteology of the Triassic turtle Proganochelys. Bull Amer Mus Nat Hist. 1990, 194: 1-263.

    Google Scholar 

  27. Girondot M, Sire J-Y: Evolution of the amelogenin gene in toothed and toothless vertebrates. Eur J Oral Sci. 1998, 106 (Suppl 1): 501-508.

    Article  CAS  PubMed  Google Scholar 

  28. Chiari Y, Cahais V, Galtier N, Delsuc F: Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria). BMC Biol. 2012, 10: 65-10.1186/1741-7007-10-65.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Gaffney ES, Hutchison JH, Jenkins FA, Meeker LJ: Modern turtle origins: the oldest cryptodire. Science. 1987, 237: 289-291. 10.1126/science.237.4812.289.

    Article  CAS  PubMed  Google Scholar 

  30. Gaffney ES, Meylan PA, Wyss A: A computer assisted analysis of the relationships of the higher categories of turtles. Cladistics. 1991, 7: 313-335. 10.1111/j.1096-0031.1991.tb00041.x.

    Article  Google Scholar 

  31. Gaffney ES, Jenkins FA: The cranial morphology of Kayentachelys, an Early Jurassic cryptodire, and the early history of turtles. Acta Zool. 2010, 91: 335-368.

    Google Scholar 

  32. Barley AJ, Spinks PQ, Thomson RC, Shaffer HB: Fourteen nuclear genes provide phylogenetic resolution for difficult nodes in the turtle tree of life. Mol Phylogenet Evol. 2010, 55: 1189-1194. 10.1016/j.ympev.2009.11.005.

    Article  PubMed  Google Scholar 

  33. Guillon J-M, Guéry L, Hulin V, Girondot M: A large phylogeny of turtles (Testudines) using molecular data. Contrib Zool. 2012, 81: 147-158.

    Google Scholar 

  34. Crawford NG, Faircloth BC, McCormack JE, Brumfield RT, Winker K, Glenn TC: More than 1000 ultraconserved elements provide evidence that turtles are the sister group of archosaurs. Biol Lett. 2012, 8: 783-786. 10.1098/rsbl.2012.0331.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Rambaut A: Se-Al: Sequence Alignment Editor, 2.0a11. 1996, http://evolve.zoo.ox.ac.uk,

    Google Scholar 

  36. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.

    Article  CAS  PubMed  Google Scholar 

  37. Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simão TLL, Stadler T, Rabosky DL, Honeycutt RL, Flynn JJ, Ingram CM, Steiner S, Williams TL, Robinson T, Burk-Herrick A, Westerman M, Ayoub NA, Springer MS, Murphy WJ: Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammalian diversification. Science. 2011, 334: 521-524. 10.1126/science.1211028.

    Article  CAS  PubMed  Google Scholar 

  38. Springer MS, Meredith RW, Gatesy J, Emerling CA, Park J, Rabosky DL, Stadler T, Steiner C, Ryder OA, Janečka JE, Fisher CA, Murphy WJ: Macroevolutionary dynamics and historical biogeography of primate diversification inferred from a supermatrix. PLoS ONE. 2012, 7: e45921-10.1371/journal.pone.0045921.

    Article  Google Scholar 

  39. Shedlock AM, Edwards SC: Amniotes (Amniota). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford Univ Press, 375-379.

    Google Scholar 

  40. Shaffer HB: Turtles (Testudines). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford Univ Press, 398-401.

    Google Scholar 

  41. Pereira SL, Baker AJ: A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock. Mol Biol Evol. 2006, 23: 1731-1740. 10.1093/molbev/msl038.

    Article  CAS  PubMed  Google Scholar 

  42. Pereira SL, Baker AJ: A molecular timescale for galliform birds accounting for uncertainty in time estimates and heterogeneity of rates of DNA substitutions across lineages and sites. Mol Phylogenet Evol. 2006, 38: 499-509. 10.1016/j.ympev.2005.07.007.

    Article  CAS  PubMed  Google Scholar 

  43. van Tuinen M: Advanced birds (Neoaves). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford Univ Press, 419-422.

    Google Scholar 

  44. Clarke JA: Morphology, phylogenetic taxonomy, and systematics ofIchthyornisandApatornis(Avialae: Ornithurae). Am Mus Novitat. 2006, 286: 1-179.

    Google Scholar 

  45. Delsuc F, Douzery EJP: Armadillos, anteaters, and sloths (Xenarthra). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford Univ Press, 475-478.

    Google Scholar 

  46. Ji Q, Luo Z-X, Yuan C-X, Wible JR, Zhang J-P, Georgi JA: The earliest known eutherian mammal. Nature. 2002, 416: 816-822. 10.1038/416816a.

    Article  CAS  PubMed  Google Scholar 

  47. Sasaki T, Takahashi K, Nikaido M, Miura S, Yasukawa Y, Okada N: First application of the SINE (short interspersed repetitive element) method to infer phylogenetic relationships in reptiles: an example from the turtle superfamily Testudinoidea. Mol Biol Evol. 2004, 21: 705-715. 10.1093/molbev/msh069.

    Article  CAS  PubMed  Google Scholar 

  48. Sire J-Y, Davit-Béal T, Delgado S, Gu X: The origin and evolution of enamel mineralization genes. Cells Tissue Organs. 2007, 186: 25-48. 10.1159/000102679.

    Article  Google Scholar 

  49. Turtle Taxonomy Working Group: An annotated list of modern turtle terminal taxa with comments on areas of taxonomic instability and recent change. Chelonian Res Monogr. 2007, 4: 173-199.

    Google Scholar 

  50. Sterli J, Joyce WG: The cranial anatomy of the lower Jurassic turtleKayentachelys apix. Acta Palaentol Polon. 2007, 52: 675-694.

    Google Scholar 

  51. Sterli J: A new, nearly complete stem turtle from the Jurassic of South America with implications for turtle evolution. Biol Lett. 2008, 4: 286-289. 10.1098/rsbl.2008.0022.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Kawasaki K: The SCPP gene repertoire in bony vertebrates and graded differences in mineralized tissues. Dev Gene Evol. 2009, 219: 147-157. 10.1007/s00427-009-0276-x.

    Article  CAS  Google Scholar 

  53. Tokita M, Chaeychrmsri W, Siruntawineti J: Developmental basis of toothlessness in turtles: insight into convergent evolution of vertebrate morphology. Evolution. 2012, 10.1111/j.1558-5646.2012.01752.x.

    Google Scholar 

  54. Hackett SJ, Kimball RT, Reddy S, Bowie RCK, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han K-L, Harshman J, Huddleston CJ, Marks BD, Miglia KJ, Moore WS, Sheldon FH, Steadman DW, Witt CC, Yuri T: A phylogenomic study of birds reveals their evolutionary history. Science. 2008, 320: 1763-1768. 10.1126/science.1157704.

    Article  CAS  PubMed  Google Scholar 

  55. Nakamura Y, Slaby I, Spahr A, Pezeshki G, Matsumoto K, Lyngstadaas SP: Ameloblastin fusion protein enhances pulpal healing and dentin formation in porcine teeth. Calcif Tissue Int. 2006, 2006 (78): 278-284.

    Article  Google Scholar 

  56. Tamburstuen MV, Reseland JE, Spahr A, Brookes SJ, Kvalheim G, Slaby I, Snead ML, Lyngstadaas SP: Ameloblastin expression and putative autoregulation in mesenchymal cells suggest a role in early bone formation and repair. Bone. 2011, 48: 406-413. 10.1016/j.bone.2010.09.007.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Turner P: Some observations on the dentition of the narwhal (Monodon monoceros). J Anat Physiol. 1872, 7 (Pt 1): 75-79.

    CAS  PubMed Central  PubMed  Google Scholar 

  58. Ishiyama M: Enamel structure in odontocete whales. Scanning Microsc. 1987, 1: 1071-1079.

    CAS  PubMed  Google Scholar 

  59. Nweeia MT, Eichmiller FC, Hauschka PV, Tyler E, Mead JG, Potter CW, Angnatsiak DP, Richard PR, Orr JR, Black SR: Vestigial tooth anatomy and tusk nomenclature forMonodon monoceros. Anat Rec. 2012, 295: 1006-1016. 10.1002/ar.22449.

    Article  Google Scholar 

  60. Zeichner-David M, Chen L-S, Hsu Z, Reyna J, Catón J, Bringas P: Amelogenin and ameloblastin show growth-factor like activity in periodontal ligament cells. Eur J Oral Sci. 2006, 114 (Suppl 1): 244-253.

    Article  CAS  PubMed  Google Scholar 

  61. Vymětal J, Slabý I, Spahr A, Vondrášek J, Lyngstadaas SP: Bioinformatic analysis and molecular modeling of human ameloblastin suggest a two-domain intrinsically unstructured calcium-binding protein. Eur J Oral Sci. 2008, 116: 124-134. 10.1111/j.1600-0722.2008.00526.x.

    Article  PubMed  Google Scholar 

  62. Deutsch D, Haze-Filderman A, Blumenfeld A, Dafni L, Leiser Y, Shay B, Gruenbaum-Cohen Y, Rosenfeld E, Fermon E, Zimmermann B, Haegewald S, Bernimoulin J-P, Taylor AL: Amelogenin, a major structural protein in mineralizing enamel, is also expressed in soft tissues: brain and cells of the hematopoietic system. Eur J Oral Sci. 2006, 114 (Suppl 1): 183-189.

    Article  CAS  PubMed  Google Scholar 

  63. Haze A, Taylor AL, Blumenfeld A, Rosenfeld E, Leiser Y, Dafni L, Shay B, Gruenbaum-Cohen Y, Fermon E, Haegewald S, Bernimoulin J-P, Deutsch D: Amelogenin expression in long bone and cartilage cells and in bone marrow progenitor cells. Anat Rec. 2007, 290: 455-460. 10.1002/ar.20520.

    Article  CAS  Google Scholar 

  64. Delgado S, Ishiyama Y, Sire J-Y: Validation of amelogenesis imperfecta inferred from amelogenin evolution. J Dent Res. 2007, 86: 326-330. 10.1177/154405910708600405.

    Article  CAS  PubMed  Google Scholar 

  65. Al-Hashimi N, Sire J-Y, Delgado S: Evolutionary analysis of mammalian enamelin, the largest enamel protein, supports a crucial role for the 32-kDa peptide and reveals selective adaptation in rodents and primates. J Mol Evol. 2009, 69: 635-656. 10.1007/s00239-009-9302-x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by NSF (EF0629860 to MSS and JG; DEB0743724 to JG). We thank four anonymous referees for helpful comments on an earlier draft of this manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mark S Springer.

Additional information

Competing interests

The authors declare that they have no competing financial interests.

Authors’ contributions

JG, RWM, and MSS conceived and designed the study. RWM and MSS performed the analyses. MSS wrote the manuscript. JG and RWM provided comments on the manuscript. JG and MSS prepared figures. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2012_2228_MOESM1_ESM.nex

Additional file 1: Nexus alignment of ENAM exon 9 sequences for two crocodylians, one testudine, and five birds. The stop codon in crocodylians and Chrysemys picta occurs at positions 3788–3790. A 406 bp deletion that is shared by all five birds occurs near the 3’ end of the coding region at positions 3353–3758. Additional frameshift mutations are shared by four of five birds including representatives of Anseriformes, Galliformes, and Passeriformes (2634–2637). Cry SINE elements [47] in Chrysemys occur at positions 994–1138 and 4782–4994. Several long strings of Ns in Chrysemys picta were abridged to shorter strings of ten Ns. (NEX 91 KB)

12862_2012_2228_MOESM2_ESM.pdf

Additional file 2: Schematic alignment showing regions of ENAM exon 9 that were identified in Chrysemys picta . Green bars below sequence numbers show regions of sequence similarity between C. picta and Alligator mississippiensis. Red rectangles = frameshift mutations in C. picta; light blue rectangles = indels within the coding sequence that are in multiples of three base pairs; CRY = Cry SINE insertion; red arrow = position of stop codon in A. mississippiensis. (PDF 221 KB)

12862_2012_2228_MOESM3_ESM.nex

Additional file 3: Nexus alignment of AMBN exon sequences for two crocodylians, one testudine, and one bird. The start and stop codons in crocodylians occur at positions 127–129 (exon 2) and 1374–1376 (exon 11), respectively. (NEX 111 KB)

12862_2012_2228_MOESM4_ESM.pdf

Additional file 4: Schematic alignment showing regions of AMBN exons 1–11 that were identified in Chrysemys picta . Green bars below sequence numbers show regions of sequence similarity between C. picta and Alligator mississippiensis. Dark blue rectangles = numbered exons in A. mississippiensis; red rectangles = frameshift mutations in C. picta; light blue rectangles = indels within the coding sequence that are in multiples of three bp; dark blue arrow = position of start codon in A. mississippiensis; red arrow = position of stop codon in A. mississippiensis. (PDF 123 KB)

12862_2012_2228_MOESM5_ESM.nex

Additional file 5: Nexus alignment of AMBN exon sequences for Homo sapiens , Orycteropus afer and six other afrotherians, Dasypus novemcinctus , and Choloepus hoffmanni . The start and stop codons in Homo occur at positions 121–123 (exon 1) and 1468–1470 (exon 13), respectively. The presumed start codon in Dasypus occurs at positions 151–153. Dasypus shares the same stop codon as Homo at positions 1468–1470, whereas the presumed stop codon in Choloepus is at positions 1477–1479. (NEX 77 KB)

12862_2012_2228_MOESM6_ESM.nex

Additional file 6: Nexus alignment of protein-coding AMBN sequences for 20 mammals that was used in selection analyses with PAML 4.5 [36].(NEX 43 KB)

12862_2012_2228_MOESM7_ESM.nex

Additional file 7: Nexus alignment of AMEL exon sequences in nexus format for two crocodylians, two testudines, and five birds. The start and stop codons in crocodylians occur at positions 69–71 (exon 1) and 731–733 (exon 5), respectively. A putative one-bp frameshift deletion that is shared by all five birds species occurs at position 242 in exon 4. (NEX 49 KB)

12862_2012_2228_MOESM8_ESM.pdf

Additional file 8: Schematic alignment showing regions of AMEL exons 1–5 that were identified in Chrysemys picta and Pelodiscus sinensis . Bars below sequence numbers show regions of sequence similarity between Alligator mississippiensis and both testudines (green) or A. mississippiensis and one testudine (chartreuse). Dark blue rectangles = numbered exons in A. mississippiensis; red rectangles = frameshift mutations in C. picta; light blue rectangles = indels within the coding sequence that are in multiples of three bp; dark blue arrow = position of start codon in A. mississippiensis; red arrow = position of stop codon in A. mississippiensis. (PDF 63 KB)

12862_2012_2228_MOESM9_ESM.nex

Additional file 9: Nexus alignment of Alligator mississippiensis AMEL exon 1 (positions 1–114), intron 1 (positions 115–813), and exon 2 (838–894) with exon 1 and flanking sequences from Chrysemys picta and Pelodiscus sinensis .(NEX 17 KB)

12862_2012_2228_MOESM10_ESM.nex

Additional file 10: Nexus alignment of AMEL intron 4 (positions 1–662) and exon 5 (positions 663–844) sequences for Chrysemys picta and Pelodiscus sinensis .(NEX 11 KB)

12862_2012_2228_MOESM11_ESM.nex

Additional file 11: Nexus alignment of AMEL sequences (exons 2–7) for Homo sapiens , Orycteropus afer and six other afrotherians ( Echinops telfairi, Chrysochloris asiatica, Elephantulus edwardii, Trichechus manatus, Loxodonta africana, Procavia capensis ), Dasypus novemcinctus , and Choloepus hoffmanni . The start and stop codons in Homo occur at positions 13–15 (exon 2) and 668–670 (exon 7), respectively. The mutated stop codon and frameshift mutation in exon 2 of Choloepus are clear in gnl|ti|1325289135 and gnl|ti|1311883086. The stop codon in exon 4 of Choloepus is confirmed by traces gnl|ti|1325336699 and gnl|ti|1312095734. The frameshift mutation in exon 6 of Choloepus is evident in trace gnl|ti|1363579107. The deletion that includes the 3’ region of intron 6 and first 14 bp of exon 7 in Dasypus is evident in gn|ti| 1886397515. (NEX 21 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Meredith, R.W., Gatesy, J. & Springer, M.S. Molecular decay of enamel matrix protein genes in turtles and other edentulous amniotes. BMC Evol Biol 13, 20 (2013). https://doi.org/10.1186/1471-2148-13-20

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-13-20

Keywords