MicroRNAs support a turtle + lizard clade

Despite much interest in amniote systematics, the origin of turtles remains elusive. Traditional morphological phylogenetic analyses place turtles outside Diapsida-amniotes whose ancestor had two fenestrae in the temporal region of the skull (among the living forms the tuatara, lizards, birds and crocodilians)-and allied with some unfenestrate-skulled (anapsid) taxa. Nonetheless, some morphological analyses place turtles within Diapsida, allied with Lepidosauria (tuatara and lizards). Most molecular studies agree that turtles are diapsids, but rather than allying them with lepidosaurs, instead place turtles near or within Archosauria (crocodilians and birds). Thus, three basic phylogenetic positions for turtles with respect to extant Diapsida are currently debated: (i) sister to Diapsida, (ii) sister to Lepidosauria, or (iii) sister to, or within, Archosauria. Interestingly, although these three alternatives are consistent with a single unrooted four-taxon tree for extant reptiles, they differ with respect to the position of the root. Here, we apply a novel molecular dataset, the presence versus absence of specific microRNAs, to the problem of the phylogenetic position of turtles and the root of the reptilian tree, and find that this dataset unambiguously supports a turtle + lepidosaur group. We find that turtles and lizards share four unique miRNA gene families that are not found in any other organisms' genome or small RNA library, and no miRNAs are found in all diapsids but not turtles, or in turtles and archosaurs but not in lizards. The concordance between our result and some morphological analyses suggests that there have been numerous morphological convergences and reversals in reptile phylogeny, including the loss of temporal fenestrae.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. MHC Class I genes in the Tuatara (Sphenodon spp.): evolution of the MHC in an ancient reptilian order

The major histocompatibility complex (MHC) is an extremely dynamic region of the genome, characterized by high polymorphism and frequent gene duplications and rearrangements. This has resulted in considerable differences in MHC organization and evolution among vertebrate lineages, particularly between birds and mammals. As nonavian reptiles are ancestral to both mammals and birds, they occupy an important phylogenetic position for understanding these differences. However, little is known about reptile MHC genes. To address this, we have characterized MHC class I sequences from the tuatara (Sphenodon spp.), the last survivor of an ancient order of reptiles, Sphenodontia. We isolated two different class I cDNA sequences, which share 93% sequence similarity with each other but are highly divergent from other vertebrate MHC genes. Southern blotting and polymerase chain reaction amplification of class I sequences from seven adult tuatara plus a family group indicate that these sequences represent at least two to three loci. Preliminary analysis of variation among individuals from an island population of tuatara indicates that these loci are highly polymorphic. Maximum likelihood analysis of reptile MHC class I sequences indicates that gene duplication has occurred within reptilian orders. However, the evolutionary relationships among sequences from different reptilian orders cannot be resolved, reflecting the antiquity of the major reptile lineages.     

AMNIOTE PHYLOGENY AND THE IMPORTANCE OF FOSSILS

Abstract— Several prominent cladists have questioned the importance of fossils in phylogenctic inference, and it is becoming increasingly popular to simply fit extinct forms, if they are considered at all, to a cladogram of Recent taxa. Gardiner's (1982) and Løvtrup's (1985) study of amniote phylogeny exemplifies this differential treatment, and we focused on that group of organisms to test the proposition that fossils cannot overturn a theory of relationships based only on the Recent biota. Our parsimony analysis of amniote phylogeny, special knowledge contributed by fossils being scrupulously avoided, led to the following best fitting classification, which is similar to the novel hypothesis Gardiner published: (lepidosaurs (turtles (mammals (birds, crocodiles)))). However, adding fossils resulted in a markedly different most parsimonious cladogram of the extant taxa: (mammals (turtles (lepidosaurs (birds, crocodiles)))). That classification is like the traditional hypothesis, and it provides a better fit to the stratigraphic record. To isolate the extinct taxa responsible for the latter classification, the data were successively partitioned with each phylogenetic analysis, and we concluded that: (1) the ingroup, not the outgroup, fossils were important; (2) synapsid, not reptile, fossils were pivotal; (3) certain synapsid fossils, not the earliest or latest, were responsible. The critical nature of the synapsid fossils seemed to lie in the particular combination of primitive and derived character slates they exhibited. Classifying those fossils, along with mammals, as the sister group to the lineage consisting of birds and crocodiles resulted in a relatively poor fit to data; one involving a 2—4 fold increase in evolutionary reversals! Thus, the importance of the critical fossils, collectively or individually, seems to reside in their relative primitive-ness, and the simplest explanation for their more conservative nature is that they have had less time to evolve. While fossils may be important in phylogenetic inference only under certain conditions, there is no compelling reason to prejudge their contribution. We urge systematists to evaluate fairly all of the available evidence.     

Microsatellite DNA loci identify individuals and provide no evidence for multiple paternity in wild tuatara (Sphenodon: Reptilia)

Six new microsatellite DNA loci are isolated from a genomic library of the sphenodontid reptile tuatara (Sphenodon) and presented here as a tool for identifying individuals for future paternity and kinship studies. These loci, in combination with four previously published loci, are sufficient to discriminate between clutch-mate siblings from Stephens Island and Brothers Island populations. These populations represent high and low levels of genetic diversity in tuatara populations respectively. An estimate of minimum number of fathers of each clutch found no evidence for multiple paternity in any clutch. These newly isolated loci complete the development of an array of genetic tools for use in tuatara to enhance ongoing conservation and management of wild, translocated and captive populations.     

Using Genes as Characters and a Parsimony Analysis to Explore the Phylogenetic Position of Turtles

The phylogenetic position of turtles within the vertebrate tree of life remains controversial. Conflicting conclusions from different studies are likely a consequence of systematic error in the tree construction process, rather than random error from small amounts of data. Using genomic data, we evaluate the phylogenetic position of turtles with both conventional concatenated data analysis and a “genes as characters” approach. Two datasets were constructed, one with seven species (human, opossum, zebra finch, chicken, green anole, Chinese pond turtle, and western clawed frog) and 4584 orthologous genes, and the second with four additional species (soft-shelled turtle, Nile crocodile, royal python, and tuatara) but only 1638 genes. Our concatenated data analysis strongly supported turtle as the sister-group to archosaurs (the archosaur hypothesis), similar to several recent genomic data based studies using similar methods. When using genes as characters and gene trees as character-state trees with equal weighting for each gene, however, our parsimony analysis suggested that turtles are possibly sister-group to diapsids, archosaurs, or lepidosaurs. None of these resolutions were strongly supported by bootstraps. Furthermore, our incongruence analysis clearly demonstrated that there is a large amount of inconsistency among genes and most of the conflict relates to the placement of turtles. We conclude that the uncertain placement of turtles is a reflection of the true state of nature. Concatenated data analysis of large and heterogeneous datasets likely suffers from systematic error and over-estimates of confidence as a consequence of a large number of characters. Using genes as characters offers an alternative for phylogenomic analysis. It has potential to reduce systematic error, such as data heterogeneity and long-branch attraction, and it can also avoid problems associated with computation time and model selection. Finally, treating genes as characters provides a convenient method for examining gene and genome evolution.     

Endocrinology of oviposition in the tuatara (Sphenodon punctatus)—II. Plasma arginine vasotocin concentrations during natural nesting

• 1.1. Plasma concentrations of the neurohypophysial arginine vasotocin (AVT) were measured during three stages of natural nesting [nest digging (N = 16), oviposition (N = 6), nest guarding (N = 6)] in the rare New Zealand reptile, the tuatara (Sphenodon punctatus).
• 2.2. Nest digging females (X ± 1 SE: 134.9 ± 15.0 pg/ml) exhibited elevated but significantly lower concentrations of plasma AVT compared to those observed during oviposition (216.0 ± 28.1 pg/ml). Nest guarding females (54.6 ± 29.8 pg/μl) had plasma AVT concentrations significantly lower than other stages.
• 3.3. These data are similar to those reporting elevated plasma AVT concentrations during oviposition in sea turtles but tuatara exhibit significantly lower plasma AVT values (10-fold difference). Plasma AVT concentrations reported for tuatara and birds are comparable.

     

Molecular clocks in reptiles: life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between body size and substitution rate for many molecular data sets. Both the generality and the cause of the negative body size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size, age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles, lizards, snakes, crocodile, and tuatara. Although this study was limited by the number of comparisons for which both sequence and life-history data were available, the results suggest that a negative body size trend in rate of molecular evolution may be a general feature of reptile molecular evolution, consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.     

Rapid molecular evolution in a living fossil

The tuatara of New Zealand is a unique reptile that coexisted with dinosaurs and has changed little morphologically from its Cretaceous relatives. Tuatara have very slow metabolic and growth rates, long generation times and slow rates of reproduction. This suggests that the species is likely to exhibit a very slow rate of molecular evolution. Our analysis of ancient and modern tuatara DNA shows that, surprisingly, tuatara have the highest rate of molecular change recorded in vertebrates. Our work also suggests that rates of neutral molecular and phenotypic evolution are decoupled.     

Phylogenetic position of turtles among amniotes: evidence from mitochondrial and nuclear genes

Maximum likelihood analysis, accounting for site-heterogeneity in evolutionary rate with the Gamma-distribution model, was carried out with amino acid sequences of 12 mitochondrial proteins and nucleotide sequences of mitochondrial 12S and 16S rRNAs from three turtles, one squamate, one crocodile, and eight birds. The analysis strongly suggests that turtles are closely related to archosaurs (birds+crocodilians), and it supports both Tree-2: (((birds, crocodilians), turtles), squamates) and Tree-3: ((birds, (crocodilians, turtles)), squamates). A more traditional Tree-1: (((birds, crocodilians), squamates), turtles) and a tree in which turtles are basal to other amniotes were rejected with high statistical significance. Tree-3 has recently been proposed by Hedges and Poling [Science 283 (1999) 998-1001] based mainly on nuclear genes. Therefore, we re-analyzed their data using the maximum likelihood method, and evaluated the total evidence of the analyses of mitochondrial and nuclear data sets. Tree-1 was again rejected strongly. The most likely hypothesis was Tree-3, though Tree-2 remained a plausible candidate.     

Molecular evidence for a clade of turtles.

Although turtles have been generally grouped with the most primitive reptile species, the origin and phylogenetic relationships of turtles have remained unresolved to date. To confirm the phylogenetic position of turtles in amniotes, we have cloned and determined the cDNA sequences encoding for skink lactate dehydrogenase (LDH)-A and LDH-B, snake LDH-A, and African clawed frog LDH-A; four alpha-enolase cDNA sequences from turtle, alligator, skink, and snake were also cloned and determined. All of these eight cDNA sequences, as well as the previously published LDH-A, LDH-B, and alpha-enolase of mammals, birds, reptiles, and African clawed frog, were analyzed by the phylogenetic tree reconstruction methods of neighbor-joining, maximum parsimony, and maximum likelihood. In the phylogenetic analyses, the turtle was found to be closely related to the alligator. Also, we found that the turtle had diverged after the divergence of squamates and birds. This departs from previous hypotheses of turtle evolution and further suggests that turtles are the latest of divergent reptiles, having been derived from an ancestor of crocodilian lineage within the last 200 million years.     

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.     

Deleterious Mutations of a Claw Keratin in Multiple Taxa of Reptiles

We have recently shown that homologs of mammalian hair keratins are expressed in the claws of the green anole lizard, Anolis carolinensis. To test whether reptilian hair keratin homologs are functionally associated with claws, we investigated the conservation of the prototypical reptilian hair keratin homolog, hard acidic keratin 1 (HA1), in representative species from all main clades of reptiles. A complete cDNA of HA1 was cloned from the claw-forming epidermis of the lacertid lizard Podarcis sicula, and partial HA1 gene sequences could be amplified from genomic DNA of tuatara, lizards, gekkos, turtles, and crocodiles. In contrast, the HA1 gene of the limbless slow worm, Anguis fragilis, and of two species of turtles contained at least one deleterious mutation. Moreover, an HA1 gene was undetectable in the softshell turtle, snakes, and birds. Mapping the presence and absence of HA1 onto the phylogenetic tree of sauropsids suggested that the HA1 gene has been lost independently in several lineages of reptiles. The species distribution of HA1 is compatible with the hypothesis of a primary function of HA1 in claws but also shows that the formation of reptilian claws does not strictly depend on this keratin.     

Sources and timing of calcium mobilization during embryonic development of the corn snake, Pantherophis guttatus

Embryos of oviparous Reptilia (=turtles, lepidosaurs, crocodilians and birds) extract calcium for growth and development from reserves in the yolk and eggshell. Yolk provides most of the calcium to embryos of lizards and snakes. In contrast, the eggshell supplies most of the calcium for embryonic development of turtles, crocodilians and birds. The yolk sac and chorioallantoic membrane of birds recover and transport calcium from the yolk and eggshell and homologous membranes of squamates (lizards and snakes) probably transport calcium from these two sources as well. We studied calcium mobilization by embryos of the snake Pantherophis guttatus during the interval of greatest embryonic growth and found that the pattern of calcium transfer was similar to other snakes. Calcium recovery from the yolk is relatively low until the penultimate embryonic stage. Calcium removal from the eggshell begins during the same embryonic stage and total eggshell calcium drops in each of the final 2 weeks prior to hatching. The eggshell supplies 28% of the calcium of hatchlings. The timing of calcium transport from the yolk and eggshell is coincident with the timing of growth of the yolk sac and chorioallantoic membrane and expression of the calcium binding protein, calbindin-D28K, in these tissues as reported in previous studies. In the context of earlier work, our findings suggest that the timing and mechanism of calcium transport from the yolk sac of P. guttatus is similar to birds, but that both the timing and mechanism of calcium transport by the chorioallantoic membrane differs. Based on the coincident timing of eggshell calcium loss and embryonic calcium accumulation, we also conclude that recovery of eggshell calcium in P. guttatus is regulated by the embryo.     

Transitional fossils and the origin of turtles

The origin of turtles is one of the most contentious issues in systematics with three currently viable hypotheses: turtles as the extant sister to (i) the crocodile–bird clade, (ii) the lizard–tuatara clade, or (iii) Diapsida (a clade composed of (i) and (ii)). We reanalysed a recent dataset that allied turtles with the lizard–tuatara clade and found that the inclusion of the stem turtle Proganochelys quenstedti and the ‘parareptile’ Eunotosaurus africanus results in a single overriding morphological signal, with turtles outside Diapsida. This result reflects the importance of transitional fossils when long branches separate crown clades, and highlights unexplored issues such as the role of topological congruence when using fossils to calibrate molecular clocks.     

Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink: statistical evidence for archosaurian affinity of turtles.

Turtles have highly specialized morphological characteristics, and their phylogenetic position has been under intensive debate. Previous molecular studies have not established a consistent and statistically well supported conclusion on this issue. In order to address this, complete mitochondrial DNA sequences were determined for the green turtle and the blue-tailed mole skink. These genomes possess an organization of genes which is typical of most other vertebrates, such as placental mammals, a frog, and bony fishes, but distinct from organizations of alligators and snakes. Molecular evolutionary rates of mitochondrial protein sequences appear to vary considerably among major reptilian lineages, with relatively rapid rates for snake and crocodilian lineages but slow rates for turtle and lizard lineages. In spite of this rate heterogeneity, phylogenetic analyses using amino acid sequences of 12 mitochondrial proteins reliably established the Archosauria (birds and crocodilians) and Lepidosauria (lizards and snakes) clades postulated from previous morphological studies. The phylogenetic analyses further suggested that turtles are a sister group of the archosaurs, and this untraditional relationship was provided with strong statistical evidence by both the bootstrap and the Kishino-Hasegawa tests. This is the first statistically significant molecular phylogeny on the placement of turtles relative to the archosaurs and lepidosaurs. It is therefore likely that turtles originated from a Permian-Triassic archosauromorph ancestor with two pairs of temporal fenestrae behind the skull orbit that were subsequently lost. The traditional classification of turtles in the Anapsida may thus need to be reconsidered.     

Timing of organogenesis support basal position of turtles in the amniote tree of life

Background  

The phylogenetic position of turtles is the most disputed aspect in the reconstruction of the land vertebrate tree of life. This controversy has arisen after many different kinds and revisions of investigations of molecular and morphological data. Three main hypotheses of living sister-groups of turtles have resulted from them: all reptiles, crocodiles + birds or squamates + tuatara. Although embryology has played a major role in morphological studies of vertebrate phylogeny, data on developmental timing have never been examined to explore and test the alternative phylogenetic hypotheses. We conducted a comprehensive study of published and new embryological data comprising 15 turtle and eight tetrapod species belonging to other taxa, integrating for the first time data on the side-necked turtle clade.     

Reduced MHC variation in a threatened tuatara species

The relationship between neutral and adaptive genetic diversity is important to understand in assessing the implications of a population bottleneck. Fitness-related genes, such as those of the major histocompatibility complex (MHC), may be influenced by selection, and so retain diversity even when it is lost at neutral markers. We measured MHC class I variation in an archaic reptile species Sphenodon guntheri [North Brother Island (NBI) tuatara], which naturally occurs on one 4 ha island in Cook Strait, New Zealand, and has low levels of microsatellite diversity. MHC variation in S. guntheri was compared with microsatellite DNA variation, and with MHC variation in a large population of Sphenodon punctatus (Cook Strait tuatara) on Stephens Island. The NBI population shows significantly decreased levels of genetic diversity compared with the Stephens Island population. Only three different MHC sequences and three genotypes were found on NBI, compared with 15 sequences and 21 genotypes in a similar sample size from Stephens Island. Two sequences appear to be unique to the NBI population. The assortment of sequence variants into genotypes suggests strong gametic disequilibrium between two MHC class I loci in S. guntheri , and only two haplotypes that were present in Hardy–Weinberg proportions were identified. MHC diversity in NBI tuatara appears to be largely influenced by genetic drift, consistent with a recent population bottleneck. This may compromise the ability of this population to respond to novel disease threats.