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

Dinosaurs: Comparative Cytogenomics of Their Reptile Cousins and Avian Descendants

Reptiles known as dinosaurs pervade scientific and popular culture, while interest in their genomics has increased since the 1990s. Birds (part of the crown group Reptilia) are living theropod dinosaurs. Chromosome-level genome assemblies cannot be made from long-extinct biological material, but dinosaur genome organization can be inferred through comparative genomics of related extant species. Most reptiles apart from crocodilians have both macro- and microchromosomes; comparative genomics involving molecular cytogenetics and bioinformatics has established chromosomal relationships between many species. The capacity of dinosaurs to survive multiple extinction events is now well established, and birds now have more species in comparison with any other terrestrial vertebrate. This may be due, in part, to their karyotypic features, including a distinctive karyotype of around n = 40 (~10 macro and 30 microchromosomes). Similarity in genome organization in distantly related species suggests that the common avian ancestor had a similar karyotype to e.g., the chicken/emu/zebra finch. The close karyotypic similarity to the soft-shelled turtle (n = 33) suggests that this basic pattern was mostly established before the Testudine-Archosaur divergence, ~255 MYA. That is, dinosaurs most likely had similar karyotypes and their extensive phenotypic variation may have been mediated by increased random chromosome segregation and genetic recombination, which is inherently higher in karyotypes with more and smaller chromosomes.

The evolution of dermal shield vascularization in Testudinata and Pseudosuchia: phylogenetic constraints versus ecophysiological adaptations

Studies on living turtles have demonstrated that shells are involved in the resistance to hypoxia during apnea via bone acidosis buffering; a process which is complemented with cutaneous respiration, transpharyngeal and cloacal gas exchanges in the soft-shell turtles. Bone acidosis buffering during apnea has also been identified in crocodylian osteoderms, which are also known to employ heat transfer when basking. Although diverse, many of these functions rely on one common trait: the vascularization of the dermal shield. Here, we test whether the above ecophysiological functions played an adaptive role in the evolutionary transitions between land and aquatic environments in both Pseudosuchia and Testudinata. To do so, we measured the bone porosity as a proxy for vascular density in a set of dermal plates before performing phylogenetic comparative analyses. For both lineages, the dermal plate porosity obviously varies depending on the animal lifestyle, but these variations prove to be highly driven by phylogenetic relationships. We argue that the complexity of multi-functional roles of the post-cranial dermal skeleton in both Pseudosuchia and Testudinata probably is the reason for a lack of obvious physiological signal, and we discuss the role of the dermal shield vascularization in the evolution of these groups. This article is part of the theme issue 'Vertebrate palaeophysiology'.

The Complete Mitochondrial Genome of Platysternon megacephalum peguense and Molecular Phylogenetic Analysis

Platysternon megacephalum is the only living representative species of Platysternidae and only three subspecies remain: P. m. megalorcephalum, P. m. shiui, and P. m. peguense. However, previous reports implied that P. m. peguense has distinct morphological and molecular features. The characterization of the mitogenome has been accepted as an efficient means of phylogenetic and evolutionary analysis. Hence, this study first determined the complete mitogenome of P. m. peguense with the aim to identify the structure and variability of the P. m. peguense mitogenome through comparative analysis. Furthermore, the phylogenetic relationship of the three subspecies was tested. Based on different tRNA gene loss and degeneration of these three subspecies, their rearrangement pathways have been inferred. Phylogenetic analysis showed that P. m. peguense is a sister group to (P. m. megalorcephalum and P. m. shiui). Furthermore, the divergence time estimation of these three subspecies coincided with the uplift of the Tibetan Plateau. This study shows that the genetic distances between P. m. peguense and the other two subspecies are comparable to interspecific genetic distances, for example within Mauremys. In general, this study provides new and meaningful insights into the evolution of the three Platysternidae subspecies.

Information Criteria for Comparing Partition Schemes

When inferring phylogenies, one important decision is whether and how nucleotide substitution parameters should be shared across different subsets or partitions of the data. One sort of partitioning error occurs when heterogeneous subsets are mistakenly lumped together and treated as if they share parameter values. The opposite kind of error is mistakenly treating homogeneous subsets as if they result from distinct sets of parameters. Lumping and splitting errors are not equally bad. Lumping errors can yield parameter estimates that do not accurately reflect any of the subsets that were combined whereas splitting errors yield estimates that did not benefit from sharing information across partitions. Phylogenetic partitioning decisions are often made by applying information criteria such as the Akaike information criterion (AIC). As with other information criteria, the AIC evaluates a model or partition scheme by combining the maximum log-likelihood value with a penalty that depends on the number of parameters being estimated. For the purpose of selecting an optimal partitioning scheme, we derive an adjustment to the AIC that we refer to as the AIC$^{(p)}$ and that is motivated by the idea that splitting errors are less serious than lumping errors. We also introduce a similar adjustment to the Bayesian information criterion (BIC) that we refer to as the BIC$^{(p)}$. Via simulation and empirical data analysis, we contrast AIC and BIC behavior to our suggested adjustments. We discuss these results and also emphasize why we expect the probability of lumping errors with the AIC$^{(p)}$ and the BIC$^{(p)}$ to be relatively robust to model parameterization.

Seeing red to being red: conserved genetic mechanism for red cone oil droplets and co-option for red coloration in birds and turtles

Avian ketocarotenoid pigments occur in both the red retinal oil droplets that contribute to colour vision and bright red coloration used in signalling. Turtles are the only other tetrapods with red retinal oil droplets, and some also display red carotenoid-based coloration. Recently, the CYP2J19 gene was strongly implicated in ketocarotenoid synthesis in birds. Here, we investigate CYP2J19 evolution in relation to colour vision and red coloration in reptiles using genomic and expression data. We show that turtles, but not crocodiles or lepidosaurs, possess a CYP2J19 orthologue, which arose via gene duplication before turtles and archosaurs split, and which is strongly and specifically expressed in the ketocarotenoid-containing retina and red integument. We infer that CYP2J19 initially functioned in colour vision in archelosaurs and conclude that red ketocarotenoid-based coloration evolved independently in birds and turtles via gene regulatory changes of CYP2J19 Our results suggest that red oil droplets contributed to colour vision in dinosaurs and pterosaurs.

The origin of turtles: a paleontological perspective

The origin of turtles and their unusual body plan has fascinated scientists for the last two centuries. Over the course of the last decades, a broad sample of molecular analyses have favored a sister group relationship of turtles with archosaurs, but recent studies reveal that this signal may be the result of systematic biases affecting molecular approaches, in particular sampling, non-randomly distributed rate heterogeneity among taxa, and the use of concatenated data sets. Morphological studies, by contrast, disfavor archosaurian relationships for turtles, but the proposed alternative topologies are poorly supported as well. The recently revived paleontological hypothesis that the Middle Permian Eunotosaurus africanus is an intermediate stem turtle is now robustly supported by numerous characters that were previously thought to be unique to turtles and that are now shown to have originated over the course of tens of millions of years unrelated to the origin of the turtle shell. Although E. africanus does not solve the placement of turtles within Amniota, it successfully extends the stem lineage of turtles to the Permian and helps resolve some questions associated with the origin of turtles, in particular the non-composite origin of the shell, the slow origin of the shell, and the terrestrial setting for the origin of turtles.

Complete mitochondrial genome of Cuora trifasciata (Chinese three-striped box turtle), and a comparative analysis with other box turtles

Cuora trifasciata has become one of the most critically endangered species in the world. The complete mitochondrial genome of C. trifasciata (Chinese three-striped box turtle) was determined in this study. Its mitochondrial genome is a 16,575-bp-long circular molecule that consists of 37 genes that are typically found in other vertebrates. And the basic characteristics of the C. trifasciata mitochondrial genome were also determined. Moreover, a comparison of C. trifasciata with Cuora cyclornata, Cuora pani and Cuora aurocapitata indicated that the four mitogenomics differed in length, codons, overlaps, 13 protein-coding genes (PCGs), ND3, rRNA genes, control region, and other aspects. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 12 protein-coding genes of the genus Cuora indicated the phylogenetic position of C. trifasciata within Cuora. The phylogenetic analysis also showed that C. trifasciata from Vietnam and China formed separate monophyletic clades with different Cuora species. The results of nucleotide base compositions, protein-coding genes and phylogenetic analysis showed that C. trifasciata from these two countries may represent different Cuora species.

The hooked element in the pes of turtles (Testudines): a global approach to exploring primary and secondary homology

The hooked element in the pes of turtles was historically identified by most palaeontologists and embryologists as a modified fifth metatarsal, and often used as evidence to unite turtles with other reptiles with a hooked element. Some recent embryological studies, however, revealed that this element might represent an enlarged fifth distal tarsal. We herein provide extensive new myological and developmental observations on the hooked element of turtles, and re-evaluate its primary and secondary homology using all available lines of evidence. Digital count and timing of development are uninformative. However, extensive myological, embryological and topological data are consistent with the hypothesis that the hooked element of turtles represents a fusion of the fifth distal tarsal with the fifth metatarsal, but that the fifth distal tarsal dominates the hooked element in pleurodiran turtles, whereas the fifth metatarsal dominates the hooked element of cryptodiran turtles. The term 'ansulate bone' is proposed to refer to hooked elements that result from the fusion of these two bones. The available phylogenetic and fossil data are currently insufficient to clarify the secondary homology of hooked elements within Reptilia.

Production of chimeras between the Chinese soft-shelled turtle and Peking duck through transfer of early blastoderm cells

Chimeras are useful models for studies of developmental biology and cell differentiation. Intraspecies and interspecies germline chimeras have been produced in previous studies, but the feasibility of producing chimeras between animals of two different classes remains unclear. To address this issue, we attempted to produce chimeras between the Chinese soft-shelled turtle and the Peking duck by transferring stage X blastoderm cells to recipient embryos. We then examined the survival and development of the PKH26-labeled donor cells in the heterologous embryos. At early embryonic stages, both turtle and duck donor cells that were labeled with PKH26 were readily observed in the brain, neural tube, heart and gonads of the respective recipient embryos. Movement of turtle donor-derived cells was observed in the duck host embryos after 48 h of incubation. Although none of the hatchlings presented a chimeric phenotype, duck donor-derived cells were detected in a variety of organs in the hatchling turtles, particularly in the gonads. Moreover, in the hatched turtles, mRNA expression of tissue-specific duck genes MEF2a and MEF2c was detected in many tissues, including the muscle, heart, small and large intestines, stomach and kidney. Similarly, SPAG6 mRNA was detected in a subset of turtle tissues, including the gonad and the small and large intestines. These results suggest that duck donor-derived cells can survive and differentiate in recipient turtles; however, no turtle-derived cells were detected in the hatched ducks. Our findings indicate that chimeras can be produced between animals of two different classes.

Intra-genomic GC heterogeneity in sauropsids: evolutionary insights from cDNA mapping and GC(3) profiling in snake

BACKGROUND

Extant sauropsids (reptiles and birds) are divided into two major lineages, the lineage of Testudines (turtles) and Archosauria (crocodilians and birds) and the lineage of Lepidosauria (tuatara, lizards, worm lizards and snakes). Karyotypes of these sauropsidan groups generally consist of macrochromosomes and microchromosomes. In chicken, microchromosomes exhibit a higher GC-content than macrochromosomes. To examine the pattern of intra-genomic GC heterogeneity in lepidosaurian genomes, we constructed a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 183 cDNA clones by fluorescence in situ hybridization, and examined the correlation between the GC-content of exonic third codon positions (GC3) of the genes and the size of chromosomes on which the genes were localized.

RESULTS

Although GC3 distribution of snake genes was relatively homogeneous compared with those of the other amniotes, microchromosomal genes showed significantly higher GC3 than macrochromosomal genes as in chicken. Our snake cytogenetic map also identified several conserved segments between the snake macrochromosomes and the chicken microchromosomes. Cross-species comparisons revealed that GC3 of most snake orthologs in such macrochromosomal segments were GC-poor (GC3 < 50%) whereas those of chicken orthologs in microchromosomes were relatively GC-rich (GC3 ≥ 50%).

CONCLUSION

Our results suggest that the chromosome size-dependent GC heterogeneity had already occurred before the lepidosaur-archosaur split, 275 million years ago. This character was probably present in the common ancestor of lepidosaurs and but lost in the lineage leading to Anolis during the diversification of lepidosaurs. We also identified several genes whose GC-content might have been influenced by the size of the chromosomes on which they were harbored over the course of sauropsid evolution.

The head and neck anatomy of sea turtles (Cryptodira: Chelonioidea) and skull shape in Testudines

BACKGROUND

Sea turtles (Chelonoidea) are a charismatic group of marine reptiles that occupy a range of important ecological roles. However, the diversity and evolution of their feeding anatomy remain incompletely known.

METHODOLOGY/PRINCIPAL FINDINGS

Using computed tomography and classical comparative anatomy we describe the cranial anatomy in two sea turtles, the loggerhead (Caretta caretta) and Kemp's ridley (Lepidochelys kempii), for a better understanding of sea turtle functional anatomy and morphological variation. In both taxa the temporal region of the skull is enclosed by bone and the jaw joint structure and muscle arrangement indicate that palinal jaw movement is possible. The tongue is relatively small, and the hyoid apparatus is not as conspicuous as in some freshwater aquatic turtles. We find several similarities between the muscles of C. caretta and L. kempii, but comparison with other turtles suggests only one of these characters may be derived: connection of the m. adductor mandibulae internus into the Pars intramandibularis via the Zwischensehne. The large fleshy origin of the m. adductor mandibulae externus Pars superficialis from the jugal seems to be a characteristic feature of sea turtles.

CONCLUSIONS/SIGNIFICANCE

In C. caretta and L. kempii the ability to suction feed does not seem to be as well developed as that found in some freshwater aquatic turtles. Instead both have skulls suited to forceful biting. This is consistent with the observation that both taxa tend to feed on relatively slow moving but sometimes armoured prey. The broad fleshy origin of the m. adductor mandibulae externus Pars superficialis may be linked to thecheek region being almost fully enclosed in bone but the relationship is complex.

Efficiency of nuclear and mitochondrial markers recovering and supporting known amniote groups

We have analysed the efficiency of all mitochondrial protein coding genes and six nuclear markers (Adora3, Adrb2, Bdnf, Irbp, Rag2 and Vwf) in reconstructing and statistically supporting known amniote groups (murines, rodents, primates, eutherians, metatherians, therians). The efficiencies of maximum likelihood, Bayesian inference, maximum parsimony, neighbor-joining and UPGMA were also evaluated, by assessing the number of correct and incorrect recovered groupings. In addition, we have compared support values using the conservative bootstrap test and the Bayesian posterior probabilities. First, no correlation was observed between gene size and marker efficiency in recovering or supporting correct nodes. As expected, tree-building methods performed similarly, even UPGMA that, in some cases, outperformed other most extensively used methods. Bayesian posterior probabilities tend to show much higher support values than the conservative bootstrap test, for correct and incorrect nodes. Our results also suggest that nuclear markers do not necessarily show a better performance than mitochondrial genes. The so-called dependency among mitochondrial markers was not observed comparing genome performances. Finally, the amniote groups with lowest recovery rates were therians and rodents, despite the morphological support for their monophyletic status. We suggest that, regardless of the tree-building method, a few carefully selected genes are able to unfold a detailed and robust scenario of phylogenetic hypotheses, particularly if taxon sampling is increased.

Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria)

BACKGROUND

The morphological peculiarities of turtles have, for a long time, impeded their accurate placement in the phylogeny of amniotes. Molecular data used to address this major evolutionary question have so far been limited to a handful of markers and/or taxa. These studies have supported conflicting topologies, positioning turtles as either the sister group to all other reptiles, to lepidosaurs (tuatara, lizards and snakes), to archosaurs (birds and crocodiles), or to crocodilians. Genome-scale data have been shown to be useful in resolving other debated phylogenies, but no such adequate dataset is yet available for amniotes.

RESULTS

In this study, we used next-generation sequencing to obtain seven new transcriptomes from the blood, liver, or jaws of four turtles, a caiman, a lizard, and a lungfish. We used a phylogenomic dataset based on 248 nuclear genes (187,026 nucleotide sites) for 16 vertebrate taxa to resolve the origins of turtles. Maximum likelihood and Bayesian concatenation analyses and species tree approaches performed under the most realistic models of the nucleotide and amino acid substitution processes unambiguously support turtles as a sister group to birds and crocodiles. The use of more simplistic models of nucleotide substitution for both concatenation and species tree reconstruction methods leads to the artefactual grouping of turtles and crocodiles, most likely because of substitution saturation at third codon positions. Relaxed molecular clock methods estimate the divergence between turtles and archosaurs around 255 million years ago. The most recent common ancestor of living turtles, corresponding to the split between Pleurodira and Cryptodira, is estimated to have occurred around 157 million years ago, in the Upper Jurassic period. This is a more recent estimate than previously reported, and questions the interpretation of controversial Lower Jurassic fossils as being part of the extant turtles radiation.

CONCLUSIONS

These results provide a phylogenetic framework and timescale with which to interpret the evolution of the peculiar morphological, developmental, and molecular features of turtles within the amniotes.

The development of three long universal nuclear protein-coding locus markers and their application to osteichthyan phylogenetics with nested PCR

BACKGROUND

Universal nuclear protein-coding locus (NPCL) markers that are applicable across diverse taxa and show good phylogenetic discrimination have broad applications in molecular phylogenetic studies. For example, RAG1, a representative NPCL marker, has been successfully used to make phylogenetic inferences within all major osteichthyan groups. However, such markers with broad working range and high phylogenetic performance are still scarce. It is necessary to develop more universal NPCL markers comparable to RAG1 for osteichthyan phylogenetics.

METHODOLOGY/PRINCIPAL FINDINGS

We developed three long universal NPCL markers (>1.6 kb each) based on single-copy nuclear genes (KIAA1239, SACS and TTN) that possess large exons and exhibit the appropriate evolutionary rates. We then compared their phylogenetic utilities with that of the reference marker RAG1 in 47 jawed vertebrate species. In comparison with RAG1, each of the three long universal markers yielded similar topologies and branch supports, all in congruence with the currently accepted osteichthyan phylogeny. To compare their phylogenetic performance visually, we also estimated the phylogenetic informativeness (PI) profile for each of the four long universal NPCL markers. The PI curves indicated that SACS performed best over the whole timescale, while RAG1, KIAA1239 and TTN exhibited similar phylogenetic performances. In addition, we compared the success of nested PCR and standard PCR when amplifying NPCL marker fragments. The amplification success rate and efficiency of the nested PCR were overwhelmingly higher than those of standard PCR.

CONCLUSIONS/SIGNIFICANCE

Our work clearly demonstrates the superiority of nested PCR over the conventional PCR in phylogenetic studies and develops three long universal NPCL markers (KIAA1239, SACS and TTN) with the nested PCR strategy. The three markers exhibit high phylogenetic utilities in osteichthyan phylogenetics and can be widely used as pilot genes for phylogenetic questions of osteichthyans at different taxonomic levels.

More than 1000 ultraconserved elements provide evidence that turtles are the sister group of archosaurs

We present the first genomic-scale analysis addressing the phylogenetic position of turtles, using over 1000 loci from representatives of all major reptile lineages including tuatara. Previously, studies of morphological traits positioned turtles either at the base of the reptile tree or with lizards, snakes and tuatara (lepidosaurs), whereas molecular analyses typically allied turtles with crocodiles and birds (archosaurs). A recent analysis of shared microRNA families found that turtles are more closely related to lepidosaurs. To test this hypothesis with data from many single-copy nuclear loci dispersed throughout the genome, we used sequence capture, high-throughput sequencing and published genomes to obtain sequences from 1145 ultraconserved elements (UCEs) and their variable flanking DNA. The resulting phylogeny provides overwhelming support for the hypothesis that turtles evolved from a common ancestor of birds and crocodilians, rejecting the hypothesized relationship between turtles and lepidosaurs.

Reptilian-transcriptome v1.0, a glimpse in the brain transcriptome of five divergent Sauropsida lineages and the phylogenetic position of turtles

Background

Reptiles are largely under-represented in comparative genomics despite the fact that they are substantially more diverse in many respects than mammals. Given the high divergence of reptiles from classical model species, next-generation sequencing of their transcriptomes is an approach of choice for gene identification and annotation.

Results

Here, we use 454 technology to sequence the brain transcriptome of four divergent reptilian and one reference avian species: the Nile crocodile, the corn snake, the bearded dragon, the red-eared turtle, and the chicken. Using an in-house pipeline for recursive similarity searches of >3,000,000 reads against multiple databases from 7 reference vertebrates, we compile a reptilian comparative transcriptomics dataset, with homology assignment for 20,000 to 31,000 transcripts per species and a cumulated non-redundant sequence length of 248.6 Mbases. Our approach identifies the majority (87%) of chicken brain transcripts and about 50% of de novo assembled reptilian transcripts. In addition to 57,502 microsatellite loci, we identify thousands of SNP and indel polymorphisms for population genetic and linkage analyses. We also build very large multiple alignments for Sauropsida and mammals (two million residues per species) and perform extensive phylogenetic analyses suggesting that turtles are not basal living reptiles but are rather associated with Archosaurians, hence, potentially answering a long-standing question in the phylogeny of Amniotes.

Conclusions

The reptilian transcriptome (freely available at http://www.reptilian-transcriptomes.org) should prove a useful new resource as reptiles are becoming important new models for comparative genomics, ecology, and evolutionary developmental genetics.

Cloning of estrogen receptor alpha and aromatase cDNAs and gene expression in turtles (Chrysemys picta and Pseudemys scripta) exposed to different environments

Reproductive changes have been observed in painted turtles (Chrysemys picta) from a site with known contamination located on Cape Cod, MA, USA. We hypothesize that these changes are caused by exposure to endocrine-disrupting compounds and that genes playing a significant role in reproduction are affected. cDNA sequences were determined for estrogen receptor alpha and aromatase in the painted turtle. These genes were measured in our study animals using quantitative PCR. Adult turtles were trapped from our study site (Moody Pond, MP) or a reference site (Washburn Pond) and exposed to laboratory environments containing soil from either site. The red-eared slider (Pseudemys scripta), a pond turtle closely related to the painted turtle, was used to assess neonatal exposure to soil and water from the sites. Our results show an increase in hepatic estrogen receptor, which suggests exposure to estrogenic contaminants. Female turtles from MP appear to have a long-term effect on hepatic ER. Other findings were apparent age-dependent differences in expression of aromatase and ER in the brains of neonate and year-old juvenile turtles. Phylogenetic analyses of the cDNA sequences further support the hypothesis that turtles are in a sister clade to birds and crocodilians.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.

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.

Comparative anatomy, homologies and evolution of the pectoral and forelimb musculature of tetrapods with special attention to extant limbed amphibians and reptiles

The main aim of the present work is to synthesize the information obtained from our dissections of the pectoral and forelimb muscles of representative members of the major extant taxa of limbed amphibians and reptiles and from our review of the literature, in order to provide an account of the comparative anatomy, homologies and evolution of these muscles in the Tetrapoda. The pectoral and forelimb musculature of all these major taxa conform to a general pattern that seems to have been acquired very early in the evolutionary history of tetrapods. Although some muscles are missing in certain taxa, and a clear departure from this general pattern is obviously present in derived groups such as birds, the same overall configuration is easily distinguishable in these taxa. Among the most notable anatomical differences between the groups, one that seems to have relevant evolutionary and functional implications, concerns the distal insertion points of the forearm musculature. In tetrapods, the muscles of the radial and ulnar complexes of the forearm are pleisomorphically mainly inserted onto the radius/ulna or onto the more proximal carpal bones, but in mammals some of these muscles insert more distally onto bones such as the metacarpals. Interestingly, a similar trend towards a more distal insertion of these muscles is also found in some non-mammalian tetrapod taxa, such as some anurans (e.g. Phyllomedusa). This may be correlated with the acquisition of more subtle digital movement abilities in these latter taxa.

Analyzing the relationship between sequence divergence and nodal support using Bayesian phylogenetic analyses

Determining the appropriate gene for phylogeny reconstruction can be a difficult process. Rapidly evolving genes tend to resolve recent relationships, but suffer from alignment issues and increased homoplasy among distantly related species. Conversely, slowly evolving genes generally perform best for deeper relationships, but lack sufficient variation to resolve recent relationships. We determine the relationship between sequence divergence and Bayesian phylogenetic reconstruction ability using both natural and simulated datasets. The natural data are based on 28 well-supported relationships within the subphylum Vertebrata. Sequences of 12 genes were acquired and Bayesian analyses were used to determine phylogenetic support for correct relationships. Simulated datasets were designed to determine whether an optimal range of sequence divergence exists across extreme phylogenetic conditions. Across all genes we found that an optimal range of divergence for resolving the correct relationships does exist, although this level of divergence expectedly depends on the distance metric. Simulated datasets show that an optimal range of sequence divergence exists across diverse topologies and models of evolution. We determine that a simple to measure property of genetic sequences (genetic distance) is related to phylogenic reconstruction ability in Bayesian analyses. This information should be useful for selecting the most informative gene to resolve any relationships, especially those that are difficult to resolve, as well as minimizing both cost and confounding information during project design.

Genome evolution in Reptilia: in silico chicken mapping of 12,000 BAC-end sequences from two reptiles and a basal bird

BACKGROUND

With the publication of the draft chicken genome and the recent production of several BAC clone libraries from non-avian reptiles and birds, it is now possible to undertake more detailed comparative genomic studies in Reptilia. Of interest in particular are the genomic events that transformed the large, repeat-rich genomes of mammals and non-avian reptiles into the minimalist chicken genome. We have used paired BAC end sequences (BESs) from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta) and emu (Dromaius novaehollandiae) to investigate patterns of sequence divergence, gene and retroelement content, and microsynteny between these species and chicken.

RESULTS

From a total of 11,967 curated BESs, we successfully mapped 725, 773 and 2597 sequences in alligator, turtle, and emu, respectively, to sites in the draft chicken genome using a stringent BLAST protocol. Most commonly, sequences mapped to a single site in the chicken genome. Of 1675, 1828 and 2936 paired BESs obtained for alligator, turtle, and emu, respectively, a total of 34 (alligator, 2%), 24 (turtle, 1.3%) and 479 (emu, 16.3%) pairs were found to map with high confidence and in the correct orientation and with BAC-sized intermarker distances to single chicken chromosomes, including 25 such paired hits in emu mapping to the chicken Z chromosome. By determining the insert sizes of a subset of BAC clones from these three species, we also found a significant correlation between the intermarker distance in alligator and turtle and in chicken, with slopes as expected on the basis of the ratio of the genome sizes.

CONCLUSION

Our results suggest that a large number of small-scale chromosomal rearrangements and deletions in the lineage leading to chicken have drastically reduced the number of detected syntenies observed between the chicken and alligator, turtle, and emu genomes and imply that small deletions occurring widely throughout the genomes of reptilian and avian ancestors led to the ~50% reduction in genome size observed in birds compared to reptiles. We have also mapped and identified likely gene regions in hundreds of new BAC clones from these species.

Species identification using the cytochrome b gene of commercial turtle shells

Turtle shells and their gelled products are familiar in some countries as foods, tonics and medicines. These shells may come from endangered and protected species, requiring the identification of the species present to enforce national and international legislation. We report on the design of five combinations of primer pairs for the identification of turtle shells and shell fragments used as ornaments, food products and medicines. The types of samples used are those encountered frequently and will typically contain highly degraded DNA. The success rate for species identification using the test described is dependent upon the choice of primer sets used and the length of the expected amplification product. Gelled products were simulated by the process of decoction for up to 12 h, after which all the turtle species could be identified from the liquid samples. This study establishes a method for the identification of commercial turtle shells and illustrates a simulated case using gelled products.

Trends in the evolution of reptilian chromosomes

Reptiles are a karyologically heterogeneous group, where some orders and suborders exhibit characteristics similar to those of anamniotes and others share similarities with homeotherms. The class also shows different evolutionary trends, for instance in genome and chromosome size and composition. The turtle DNA base composition is similar to that of mammals, whereas that of lizards and snakes is more similar to that of anamniotes. The major karyological differences between turtles and squamates are the size and composition of the genome and the rate at which chromosomes change. Turtles have larger and more variable genome sizes, and a greater amount of middle repetitive DNA that differs even among related species. In lizards and snakes size of the genome are smaller, single-copy DNA is constant within each suborder, and differences in repetitive DNA involve fractions that become increasingly heterogeneous with widening phylogenetic distance. With regard to variation in karyotype morphology, turtles and crocodiles show low variability in chromosome number, morphology, and G-banding pattern. Greater variability is found among squamates, which have a similar degree of karyotypic change-as do some mammals, such as carnivores and bats-and in which there are also differences among congeneric species. An interesting relationship has been highlighted in the entire class Reptilia between rates of change in chromosomes, number of living species, and rate of extinction. However, different situations obtain in turtles and crocodiles on the one hand, and squamates on the other. In the former, the rate of change in chromosomes is lower and the various evolutionary steps do not seem to have entailed marked chromosomal variation, whereas squamates have a higher rate of change in chromosomes clearly related to the number of living species, and chromosomal variation seems to have played an important role in the evolution of several taxa. The different evolutionary trends in chromosomes observed between turtles and crocodiles on the one hand and squamates on the other might depend on their different patterns of G-banding.

Calibration Choice, Rate Smoothing, and the Pattern of Tetrapod Diversification According to the Long Nuclear Gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.

Ribosomal RNA genes and deuterostome phylogeny revisited: More cyclostomes, elasmobranchs, reptiles, and a brittle star

This is an expanded study of the relationships among the deuterostome animals based on combined, nearly complete 28S and 18S rRNA genes (>3925 nt.). It adds sequences from 20 more taxa to the approximately 45 sequences used in past studies. Seven of the new taxa were sequenced here (brittle star Ophiomyxa, lizard Anolis, turtle Chrysemys, sixgill shark Hexanchus, electric ray Narcine, Southern Hemisphere lamprey Geotria, and Atlantic hagfish Myxine for 28S), and the other 13 were from GenBank and the literature (from a chicken, dog, rat, human, three lungfishes, and several ray-finned fishes, or Actinopterygii). As before, our alignments were based on secondary structure but did not account for base pairing in the stems of rRNA. The new findings, derived from likelihood-based tree-reconstruction methods and by testing hypotheses with parametric bootstrapping, include: (1) brittle star joins with sea star in the echinoderm clade, Asterozoa; (2) with two hagfishes and two lampreys now available, the cyclostome (jawless) fishes remain monophyletic; (3) Hexanchiform sharks are monophyletic, as Hexanchus groups with the frilled shark, Chlamydoselachus; (4) turtle is the sister taxon of all other amniotes; (5) bird is closer to the lizard than to the mammals; (6) the bichir Polypterus is in a monophyletic Actinopterygii; (7) Zebrafish Danio is the sister taxon of the other two teleosts we examined (trout and perch); (8) the South American and African lungfishes group together to the exclusion of the Australian lungfish. Other findings either upheld those of the previous rRNA-based studies (e.g., echinoderms and hemichordates group as Ambulacraria; orbitostylic sharks; batoids are not derived from any living lineage of sharks) or were obvious (monophyly of mammals, gnathostomes, vertebrates, echinoderms, etc.). Despite all these findings, the rRNA data still fail to resolve the relations among the major groups of deuterostomes (tunicates, Ambulacraria, cephalochordates and vertebrates) and of gnathostomes (chondrichthyans, lungfishes, coelacanth, actinopterygians, amphibians, and amniotes), partly because tunicates and lungfishes are rogue taxa that disrupt the tree. Nonetheless, parametric bootstrapping showed our RNA-gene data are only consistent with these dominant hypotheses: (1) deuterostomes consist of Ambulacraria plus Chordata, with Chordata consisting of tunicates and 'vertebrates plus cephalochordates'; and (2) lungfishes are the closest living relatives of tetrapods.

The complete mitochondrial genome of the Korean soft-shelled turtle Pelodiscus sinensis (Testudines, Trionychidae)

We isolated Korean soft-shelled turtle, Pelodiscus sinensis, mitochondrial DNA by long-polymerase chain reaction (long-PCR) with conserved primers and sequenced this mitochondrial genome (mitogenome) with primer walking using flanking sequences. The P. sinensis mitochondrial DNA has 17,042 bp and its structural organization is conserved compared to those of other reptiles and mammals. To unveil the phylogenetic relationship of the turtles, we used the NJ, MP, and ML analysis methods after inferring those sequences from the mitochondrial 16S rRNA gene. We also compared two P. sinensis variants from Korea and China using the mitochondrial genome. In this study, we report the basic characteristics of the P. sinensis mitochondrial genome, including structural organization and base composition of the rRNAs, tRNAs and protein-coding genes, as well as characteristics of tRNAs. These features are applicable for the study of phylogenetic relationships in turtles.

Phylogenomics of nonavian reptiles and the structure of the ancestral amniote genome

We report results of a megabase-scale phylogenomic analysis of the Reptilia, the sister group of mammals. Large-scale end-sequence scanning of genomic clones of a turtle, alligator, and lizard reveals diverse, mammal-like landscapes of retroelements and simple sequence repeats (SSRs) not found in the chicken. Several global genomic traits, including distinctive phylogenetic lineages of CR1-like long interspersed elements (LINEs) and a paucity of A-T rich SSRs, characterize turtles and archosaur genomes, whereas higher frequencies of tandem repeats and a lower global GC content reveal mammal-like features in Anolis. Nonavian reptile genomes also possess a high frequency of diverse and novel 50-bp unit tandem duplications not found in chicken or mammals. The frequency distributions of approximately 65,000 8-mer oligonucleotides suggest that rates of DNA-word frequency change are an order of magnitude slower in reptiles than in mammals. These results suggest a diverse array of interspersed and SSRs in the common ancestor of amniotes and a genomic conservatism and gradual loss of retroelements in reptiles that culminated in the minimalist chicken genome. The sequences reported in this paper have been deposited in the GenBank database (accession nos. CZ 250707-CZ 257443 and DX 390731-DX 389174).

Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion

In recent years, avian systematics has been characterized by a diminished reliance on morphological cladistics of modern taxa, intensive palaeornithogical research stimulated by new discoveries and an inundation by analyses based on DNA sequences. Unfortunately, in contrast to significant insights into basal origins, the broad picture of neornithine phylogeny remains largely unresolved. Morphological studies have emphasized characters of use in palaeontological contexts. Molecular studies, following disillusionment with the pioneering, but non-cladistic, work of Sibley and Ahlquist, have differed markedly from each other and from morphological works in both methods and findings. Consequently, at the turn of the millennium, points of robust agreement among schools concerning higher-order neornithine phylogeny have been limited to the two basalmost and several mid-level, primary groups. This paper describes a phylogenetic (cladistic) analysis of 150 taxa of Neornithes, including exemplars from all non-passeriform families, and subordinal representatives of Passeriformes. Thirty-five outgroup taxa encompassing Crocodylia, predominately theropod Dinosauria, and selected Mesozoic birds were used to root the trees. Based on study of specimens and the literature, 2954 morphological characters were defined; these characters have been described in a companion work, approximately one-third of which were multistate (i.e. comprised at least three states), and states within more than one-half of these multistate characters were ordered for analysis. Complete heuristic searches using 10 000 random-addition replicates recovered a total solution set of 97 well-resolved, most-parsimonious trees (MPTs). The set of MPTs was confirmed by an expanded heuristic search based on 10 000 random-addition replicates and a full ratchet-augmented exploration to ascertain global optima. A strict consensus tree of MPTs included only six trichotomies, i.e. nodes differing topologically among MPTs. Bootstrapping (based on 10 000 replicates) percentages and ratchet-minimized support (Bremer) indices indicated most nodes to be robust. Several fossil Neornithes (e.g. Dinornithiformes, Aepyornithiformes) were placed within the ingroup a posteriori either through unconstrained, heursitic searches based on the complete matrix augmented by these taxa separately or using backbone-constraints. Analysis confirmed the topology among outgroup Theropoda and achieved robust resolution at virtually all levels of the Neornithes. Findings included monophyly of the palaeognathous birds, comprising the sister taxa Tinamiformes and ratites, respectively, and the Anseriformes and Galliformes as monophyletic sister-groups, together forming the sister-group to other Neornithes exclusive of the Palaeognathae (Neoaves). Noteworthy inferences include: (i) the sister-group to remaining Neoaves comprises a diversity of marine and wading birds; (ii) Podicipedidae are the sister-group of Gaviidae, and not closely related to the Phoenicopteridae, as recently suggested; (iii) the traditional Pelecaniformes, including the shoebill (Balaeniceps rex) as sister-taxon to other members, are monophyletic; (iv) traditional Ciconiiformes are monophyletic; (v) Strigiformes and Falconiformes are sister-groups; (vi) Cathartidae is the sister-group of the remaining Falconiformes; (vii) Ralliformes (Rallidae and Heliornithidae) are the sister-group to the monophyletic Charadriiformes, with the traditionally composed Gruiformes and Turniciformes (Turnicidae and Mesitornithidae) sequentially paraphyletic to the entire foregoing clade; (viii) Opisthocomus hoazin is the sister-taxon to the Cuculiformes (including the Musophagidae); (ix) traditional Caprimulgiformes are monophyletic and the sister-group of the Apodiformes; (x) Trogoniformes are the sister-group of Coliiformes; (xi) Coraciiformes, Piciformes and Passeriformes are mutually monophyletic and closely related; and (xii) the Galbulae are retained within the Piciformes. Unresolved portions of the Neornithes (nodes having more than one most-parsimonious solution) comprised three parts of the tree: (a) several interfamilial nodes within the Charadriiformes; (b) a trichotomy comprising the (i) Psittaciformes, (ii) Columbiformes and (iii) Trogonomorphae (Trogoniformes, Coliiformes) + Passerimorphae (Coraciiformes, Piciformes, Passeriformes); and (c) a trichotomy comprising the Coraciiformes, Piciformes and Passeriformes. The remaining polytomies were among outgroups, although several of the highest-order nodes were only marginally supported; however, the majority of nodes were resolved and met or surpassed conventional standards of support. Quantitative comparisons with alternative hypotheses, examination of highly supportive and diagnostic characters for higher taxa, correspondences with prior studies, complementarity and philosophical differences with palaeontological phylogenetics, promises and challenges of palaeogeography and calibration of evolutionary rates of birds, and classes of promising evidence and future directions of study are reviewed. Homology, as applied to avian examples of apparent homologues, is considered in terms of recent theory, and a revised annotated classification of higher-order taxa of Neornithes and other closely related Theropoda is proposed. (c) 2007 The Linnean Society of London, Zoological Journal of the Linnean Society, 2007, 149, 1-95.

Homologies of thelongissimus,iliocostalis, and hypaxial muscles in the anterior presacral region of extant diapsida

Homologies of muscles of the m. longissimus and m. iliocostalis groups in the dorsal and cervical regions, as well as those of the subvertebral muscles and mm. intercostales externi that continue from the dorsal into the cervical regions, in extant Diapsida are proposed based on detailed dissections and published accounts of lepidosaurs, crocodylians, and birds. The morphology of tendons and innervation patterns suggest that the avian "m. iliocostalis" in the dorsal region include the homologs of both m. longissimus and m. iliocostalis in non-avian diapsids. The conserved nature of the morphology of tendons in palaeognath birds also revealed that the avian mm. intertransversarii in the cervical region consist of muscles of the both m. longissimus and m. iliocostalis groups despite having been treated as a single series of muscles, and thus are not homologous with muscles of the same name in Lepidosauria or Crocodylia. The avian mm. inclusi that lie medial to mm. intertransversarii are homologous with mm. intercostales externi in Lepidosauria and mm. intercostales externi and m. scalenus combined in Crocodylia. Innervation patterns suggest that a muscle ("m. iliocostalis capitis") connecting the atlas rib and occiput in Crocodylia includes contributions from the subvertebral layer and m. cucullaris complex, and possibly m. iliocostalis as well. The present findings may serve as a basis for revising the currently used avian nomenclature so that it will reflect homologies of muscles with their non-avian counterparts.

Phylogenomic Investigation of CR1 LINE Diversity in Reptiles

It is unlikely that taxonomically diverse phylogenetic studies will be completed rapidly in the near future for nonmodel organisms on a whole-genome basis. However, one approach to advancing the field of "phylogenomics" is to estimate the structure of poorly known genomes by mining libraries of clones from suites of taxa, rather than from single species. The present analysis adopts this approach by taking advantage of megabase-scale end-sequence scanning of reptilian genomic clones to characterize diversity of CR1-like LINEs, the dominant family of transposable elements (TEs) in the sister group of mammals. As such, it helps close an important gap in the literature on the molecular systematics and evolution of retroelements in nonavian reptiles. Results from aligning more than 14 Mb of sequence from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta), Bahamian green anole (Anolis smaragdinus), Tuatara (Sphenodon punctatus), Emu (Dromaius novaehollandiae), and Zebra Finch (Taeniopygia guttata) against a comprehensive library approximately 3000 TE-encoding peptides reflect an increasing abundance of LINE and non-long-terminal-repeat (non-LTR) retrotransposon repeat types with the age of common ancestry among exemplar reptilian clades. The hypothesis that repeat diversity is correlated with basal metabolic rate was tested using comparative methods and a significant nonlinear relationship was indicated. This analysis suggests that the age of divergence between an exemplary clade and its sister group as well as metabolic correlates should be considered in addition to genome size in explaining patterns of retroelement diversity. The first phylogenetic analysis of the largely unexplored chicken repeat 1 (CR1) 3' reverse transcriptase (RT) conserved domains 8 and 9 in nonavian reptiles reveals a pattern of multiple lineages with variable branch lengths, suggesting presence of both old and young elements and the existence of several distinct well-supported clades not apparent from previous characterization of CR1 subfamily structure in birds and the turtle. This mode of CR1 evolution contrasts with historical patterns of LINE 1 diversification in mammals and hints toward the existence of a rich but still largely unexplored diversity of nonavian retroelements of importance to advancing both comparative vertebrate genomics and amniote systematics.

Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes

All snake species exhibit genetic sex determination with the ZZ/ZW type of sex chromosomes. To investigate the origin and evolution of snake sex chromosomes, we constructed, by FISH, a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 109 cDNA clones. Eleven of the 109 clones were localized to the Z chromosome. All human and chicken homologues of the snake Z-linked genes were located on autosomes, suggesting that the sex chromosomes of snakes, mammals, and birds were all derived from different autosomal pairs of the common ancestor. We mapped the 11 Z-linked genes of E. quadrivirgata to chromosomes of two other species, the Burmese python (Python molurus bivittatus) and the habu (Trimeresurus flavoviridis), to investigate the process of W chromosome differentiation. All and 3 of the 11 clones were localized to both the Z and W chromosomes in P. molurus and E. quadrivirgata, respectively, whereas no cDNA clones were mapped to the W chromosome in T. flavoviridis. Comparative mapping revealed that the sex chromosomes are only slightly differentiated in P. molurus, whereas they are fully differentiated in T. flavoviridis, and E. quadrivirgata is at a transitional stage of sex-chromosome differentiation. The differentiation of sex chromosomes was probably initiated from the distal region on the short arm of the protosex chromosome of the common ancestor, and then deletion and heterochromatization progressed on the sex-specific chromosome from the phylogenetically primitive boids to the more advanced viperids.

Species identification of Kachuga tecta using the cytochrome b gene

A DNA technique has been established for the identification to species level of tortoises. The test on the shell of the animal was used to identify samples from the species Kachuga tecta. A total of 100 tortoise shell specimens collected from the National Council of Agriculture (COA), Taiwan, were used in this study. Primer pairs were designed to amplify partial DNA fragments of cytochrome b within the mitochondrial genome. The DNA data showed that among the 100 samples, there were four distinct haplotype DNA sequences, within which there were a total of 90 variable sites. Between haplotypes I and II, there was only 1 nucleotide difference at position 228. Between haplotypes I and III, 65 nucleotide differences were observed; haplotypes I and IV, 62 nucleotide differences; and haplotypes III and IV, 56 nucleotide differences were observed. There were 66 and 63 nucleotide differences between haplotypes II and III and haplotypes II and IV respectively. All four haplotypes were compared with the DNA sequences held at the GenBank and EMBL databases. The most similar species were K. tecta (haplotype I and II), Morenia ocellata (haplotype III) and Geoclemys hamiltonii (haplotype IV), and their respective mtDNA similarities were 99.5%, 99.3%, 89.9% and 99.5%. However, as haplotype III was only 89.9% homologous with M. ocellata, it would seem that this haplotype shows only a limited relationship with a similar species registered currently in these databases. The method established by this study is an additional method for the identification of samples protected under Convention International Trade in Endangered Species (CITES) and will improve the work for the preservation of the endangered species.

cDNA-based gene mapping and GC3 profiling in the soft-shelled turtle suggest a chromosomal size-dependent GC bias shared by sauropsids

Mammalian and avian genomes comprise several classes of chromosomal segments that vary dramatically in GC-content. Especially in chicken, microchromosomes exhibit a higher GC-content and a higher gene density than macrochromosomes. To understand the evolutionary history of the intra-genome GC heterogeneity in amniotes, it is necessary to examine the equivalence of this GC heterogeneity at the nucleotide level between these animals including reptiles, from which birds diverged. We isolated cDNAs for 39 protein-coding genes from the Chinese soft-shelled turtle, Pelodiscus sinensis, and performed chromosome mapping of 31 genes. The GC-content of exonic third positions (GC3) of P. sinensis genes showed a heterogeneous distribution, and exhibited a significant positive correlation with that of chicken and human orthologs, indicating that the last common ancestor of extant amniotes had already established a GC-compartmentalized genomic structure. Furthermore, chromosome mapping in P. sinensis revealed that microchromosomes tend to contain more GC-rich genes than GC-poor genes, as in chicken. These results illustrate two modes of genome evolution in amniotes: mammals elaborated the genomic configuration in which GC-rich and GC-poor regions coexist in individual chromosomes, whereas sauropsids (reptiles and birds) refined the chromosomal size-dependent GC compartmentalization in which GC-rich genomic fractions tend to be confined to microchromosomes.

Topics in herpesvirus genomics and evolution

Herpesviruses comprise an abundant, widely distributed group of large DNA viruses of humans and other vertebrates, and overall are among the most extensively studied large DNA viruses. Many herpesvirus genome sequences have been determined, and interpreted in terms of gene contents to give detailed views of both ubiquitous and lineage-specific functions. Availability of gene sequences has also enabled evaluations of evolutionary relationships. For herpesviruses of mammals, a robust phylogenetic tree has been constructed, which shows many features characteristic of synchronous development of virus and host lineages over large evolutionary timespans. It has also emerged that three distinct groupings of herpesviruses exist: the first containing viruses with mammals, birds and reptiles as natural hosts; the second containing viruses of amphibians and fish; and the third consisting of a single invertebrate herpesvirus. Within each of the first two groups, the genomes show clear evidence of descent from a common ancestor, but relationships between the three groups are extremely remote. Detailed analyses of capsid structures provide the best evidence for a common origin of the three groups. At a finer level, the structure of the capsid shell protein further suggests an element of common origin between herpesviruses and tailed DNA bacteriophages.

Comparative ontogeny and phylogeny of the upper jaw skeleton in amniotes

The morphology, position, and presence of the upper jaw bones vary greatly across amniote taxa. In this review, we compare the development and anatomy of upper jaw bones from the three living amniote groups: reptiles, birds, and mammals. The study of reptiles is particularly important as comparatively little is known about the embryogenesis of the jaw in this group. Our review covers the ontogeny and phylogeny of membranous bones in the face. The aim is to identify conserved embryonic processes that may exist among the three major amniote groups. Finally, we discuss how temporal and spatial regulation of preosseous condensations and ossification centers can lead to variation in the morphology of amniote upper jaw bones.

Highly conserved linkage homology between birds and turtles: bird and turtle chromosomes are precise counterparts of each other

The karyotypes of birds, turtles and snakes are characterized by two distinct chromosomal components, macrochromosomes and microchromosomes. This close karyological relationship between birds and reptiles has long been a topic of speculation among cytogeneticists and evolutionary biologists; however, there is scarcely any evidence for orthology at the molecular level. To define the conserved chromosome synteny among humans, chickens and reptiles and the process of genome evolution in the amniotes, we constructed comparative cytogenetic maps of the Chinese soft-shelled turtle (Pelodiscus sinensis) and the Japanese four-striped rat snake (Elaphe quadrivirgata) using cDNA clones of reptile functional genes. Homology between the turtle and chicken chromosomes is highly conserved, with the six largest chromosomes being almost equivalent to each other. On the other hand, homology to chicken chromosomes is lower in the snake than in the turtle. Turtle chromosome 6q and snake chromosome 2p represent conserved synteny with the chicken Z chromosome. These results suggest that the avian and turtle genomes have been well conserved during the evolution of the Arcosauria. The avian and snake sex Z chromosomes were derived from different autosomes in a common ancestor, indicating that the causative genes of sex determination may be different between birds and snakes.

Molecular and cytogenetic characterization of site-specific repetitive DNA sequences in the Chinese soft-shelled turtle (Pelodiscus sinensis, Trionychidae)

A novel family of repetitive DNA sequences that are components of constitutive heterochromatin were cloned from BglI-digested genomic DNA of the Chinese soft-shelled turtle (Pelodiscus sinensis, Trionychidae), and characterized by filter hybridization and chromosome in-situ hybridization. The BglI-family of repetitive sequences were classified into four types by their genome organization and chromosomal distribution as follows: the repeated sequences located on (1) two pairs of microchromosomes, (2) four pairs of microchromosomes,(3) about half the number of microchromosomes and (4) the interstitial region of the short arm of chromosome 2. The presence of microchromosome-specific repetitive sequences has also been reported in the Struthioniformes and Galliformes, suggesting that turtle chromosomes retain some similarity to the chromosome structure as well as the karyotypes of avian species.

Integration of Morphological Data Sets for Phylogenetic Analysis of Amniota: The Importance of Integumentary Characters and Increased Taxonomic Sampling

Several mutually exclusive hypotheses have been advanced to explain the phylogenetic position of turtles among amniotes. Traditional morphology-based analyses place turtles among extinct anapsids (reptiles with a solid skull roof), whereas more recent studies of both morphological and molecular data support an origin of turtles from within Diapsida (reptiles with a doubly fenestrated skull roof). Evaluation of these conflicting hypotheses has been hampered by nonoverlapping taxonomic samples and the exclusion of significant taxa from published analyses. Furthermore, although data from soft tissues and anatomical systems such as the integument may be particularly relevant to this problem, they are often excluded from large-scale analyses of morphological systematics. Here, conflicting hypotheses of turtle relationships are tested by (1) combining published data into a supermatrix of morphological characters to address issues of character conflict and missing data; (2) increasing taxonomic sampling by more than doubling the number of operational taxonomic units to test internal relationships within suprageneric ingroup taxa; and (3) increasing character sampling by approximately 25% by adding new data on the osteology and histology of the integument, an anatomical system that has been historically underrepresented in morphological systematics. The morphological data set assembled here represents the largest yet compiled for Amniota. Reevaluation of character data from prior studies of amniote phylogeny favors the hypothesis that turtles indeed have diapsid affinities. Addition of new ingroup taxa alone leads to a decrease in overall phylogenetic resolution, indicating that existing characters used for amniote phylogeny are insufficient to explain the evolution of more highly nested taxa. Incorporation of new data from the soft and osseous components of the integument, however, helps resolve relationships among both basal and highly nested amniote taxa. Analysis of a data set compiled from published sources and data original to this study supports monophyly of Amniota, Synapsida, Reptilia, Parareptilia, Eureptilia, Eosuchia, Diapsida, Neodiapsida, Sauria, Lepidosauria, and Archosauriformes, as well as several more highly nested divisions within the latter two clades. Turtles are here resolved as the sister taxon to a monophyletic Lepidosauria (squamates + Sphenodon), a novel phylogenetic position that nevertheless is consistent with recent molecular and morphological studies that have hypothesized diapsid affinities for this clade.

Molecular cloning of P450 aromatase from the leopard gecko and its expression in the ovary

In this study, we identified the cDNA of P450 aromatase in the leopard gecko, a lizard with temperature-dependent sex determination. The cDNA encodes a putative protein of 505 amino acids. The deduced amino acid sequence of leopard gecko aromatase cDNA showed 80% identity with that of turtles, 70% with humans and 77% with chickens. This is the first report of the identification of P450 aromatase cDNA in squamata species. It has been reported that this gene is expressed in different layers of cells in the ovary of mammalian species and avian species. Thus, we also investigated cells expressing the mRNA of this gene in the ovary of the leopard gecko by RT-PCR and in situ hybridization. The mRNA expression of leopard gecko P450 aromatase was localized in both the thecal and granulosa cell layers in the ovary. The expression in thecal and granulosa cell layers was examined in the largest follicle, second largest follicle and third largest follicle by RT-PCR. A higher level of mRNA expression was observed in the granulosa cell layer of the second largest follicle than in other cell layers. This result may reflect the characteristics of follicles in species with automonochronic ovulation.