Sister Group Relationship of Turtles to the Bird-Crocodilian Clade Revealed by Nuclear DNA–Coded Proteins

Functional characterization of follicle-stimulating hormone and luteinizing hormone receptors in cloudy catshark, Scyliorhinus torazame

The follicle-stimulating hormone receptor (FSHR) and luteinizing hormone receptor (LHR) in cloudy catshark were cloned, and recombinant FSHR and LHR were expressed for characterization. Ventral lobe extract (VLE) from the pituitary contains homologous FSH and LH, and it stimulated the cAMP signaling of FSHR and LHR dose-dependently. Two transcript variants of LHR (LHR-L with exon 10 and LHR-S without) were identified, and LHR-S was the dominant form with higher basal cAMP activity without VLE stimulation. Among various developmental stages of follicles, FSHR expression was mainly associated with the pre-vitellogenic and early white follicles. When follicles were recruited into vitellogenesis, the expression of FSHR decreased while of LHR was upregulated reciprocally, suggesting that LHR may also be responsible for the control of vitellogenesis in chondrichthyans. The expression of LHR-L was upregulated among maturing follicles before ovulation, indicating LHR-L could have a specific role in receiving the LH surge signal for final maturation. Plasma LH-like activity was transiently increased prior to the progesterone (P4)-surge and testosterone-drop at the beginning of P4-phase, supporting a pituitary control of follicle-maturation via LH signaling in chondrichthyans. The expression of follicular LHR was downregulated during the P4-phase when LH-like activity was high, indicating that the LH-dependent downregulation of LHR is conserved in chondrichthyans as it is in other vertebrate lineages. (213 words).

Evolution of ancient satellite DNAs in extant alligators and caimans (Crocodylia, Reptilia)

BACKGROUND

Crocodilians are one of the oldest extant vertebrate lineages, exhibiting a combination of evolutionary success and morphological resilience that has persisted throughout the history of life on Earth. This ability to endure over such a long geological time span is of great evolutionary importance. Here, we have utilized the combination of genomic and chromosomal data to identify and compare the full catalogs of satellite DNA families (satDNAs, i.e., the satellitomes) of 5 out of the 8 extant Alligatoridae species. As crocodilian genomes reveal ancestral patterns of evolution, by employing this multispecies data collection, we can investigate and assess how satDNA families evolve over time.

RESULTS

Alligators and caimans displayed a small number of satDNA families, ranging from 3 to 13 satDNAs in A. sinensis and C. latirostris, respectively. Together with little variation both within and between species it highlighted long-term conservation of satDNA elements throughout evolution. Furthermore, we traced the origin of the ancestral forms of all satDNAs belonging to the common ancestor of Caimaninae and Alligatorinae. Fluorescence in situ experiments showed distinct hybridization patterns for identical orthologous satDNAs, indicating their dynamic genomic placement.

CONCLUSIONS

Alligators and caimans possess one of the smallest satDNA libraries ever reported, comprising only four sets of satDNAs that are shared by all species. Besides, our findings indicated limited intraspecific variation in satellite DNA, suggesting that the majority of new satellite sequences likely evolved from pre-existing ones.

Cross-species chromosome painting and repetitive DNA mapping illuminate the karyotype evolution in true crocodiles (Crocodylidae)

Crocodilians have maintained very similar karyotype structures and diploid chromosome numbers for around 100 million years, with only minor variations in collinearity. Why this karyotype structure has largely stayed unaltered for so long is unclear. In this study, we analyzed the karyotypes of six species belonging to the genera Crocodylus and Osteolaemus (Crocodylidae, true crocodiles), among which the Congolian endemic O. osborni was included and investigated. We utilized various techniques (differential staining, fluorescence in situ hybridization with repetitive DNA and rDNA probes, whole chromosome painting, and comparative genomic hybridization) to better understand how crocodile chromosomes evolved. We studied representatives of three of the four main diploid chromosome numbers found in crocodiles (2n = 30/32/38). Our data provided new information about the species studied, including the identification of four major chromosomal rearrangements that occurred during the karyotype diversification process in crocodiles. These changes led to the current diploid chromosome numbers of 2n = 30 (fusion) and 2n = 38 (fissions), derived from the ancestral state of 2n = 32. The conserved cytogenetic tendency in crocodilians, where extant species keep near-ancestral state, contrasts with the more dynamic karyotype evolution seen in other major reptile groups.

The Evolution of Sociality and the Polyvagal Theory

The polyvagal theory (PT), offered by Porges (2021), proposes that the autonomic nervous system (ANS) was repurposed in mammals, via a "second vagal nerve", to suppress defensive strategies and support the expression of sociality. Three critical assumptions of this theory are that (1) the transition of the ANS was associated with the evolution of 'social' mammals from 'asocial' reptiles; (2) the transition enabled mammals, unlike their reptilian ancestors, to derive a biological benefit from social interactions; and (3) the transition forces a less parsimonious explanation (convergence) for the evolution of social behavior in birds and mammals, since birds evolved from a reptilian lineage. Two recently published reviews, however, provided compelling evidence that the social-asocial dichotomy is overly simplistic, neglects the diversity of vertebrate social systems, impedes our understanding of the evolution of social behavior, and perpetuates the erroneous belief that one group-non-avian reptiles-is incapable of complex social behavior. In the worst case, if PT depends upon a transition from 'asocial reptiles' to 'social mammals,' then the ability of PT to explain the evolution of the mammalian ANS is highly questionable. A great number of social behaviors occur in both reptiles and mammals. In the best case, PT has misused the terms 'social' and 'asocial'. Even here, however, the theory would still need to identify a particular suite of behaviors found in mammals and not reptiles that could be associated with, or explain, the transition of the ANS, and then replace the 'asocial' and 'social' labels with more specific descriptors.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Creating morphological diversity in reptilian temporal skull region: A review of potential developmental mechanisms

Reptilian skull morphology is highly diverse and broadly categorized into three categories based on the number and position of the temporal fenestrations: anapsid, synapsid, and diapsid. According to recent phylogenetic analysis, temporal fenestrations evolved twice independently in amniotes, once in Synapsida and once in Diapsida. Although functional aspects underlying the evolution of tetrapod temporal fenestrations have been well investigated, few studies have investigated the developmental mechanisms responsible for differences in the pattern of temporal skull region. To determine what these mechanisms might be, we first examined how the five temporal bones develop by comparing embryonic cranial osteogenesis between representative extant reptilian species. The pattern of temporal skull region may depend on differences in temporal bone growth rate and growth direction during ontogeny. Next, we compared the histogenesis patterns and the expression of two key osteogenic genes, Runx2 and Msx2, in the temporal region of the representative reptilian embryos. Our comparative analyses suggest that the embryonic histological condition of the domain where temporal fenestrations would form predicts temporal skull morphology in adults and regulatory modifications of Runx2 and Msx2 expression in osteogenic mesenchymal precursor cells are likely involved in generating morphological diversity in the temporal skull region of reptiles.

Revisiting the Karyotypes of Alligators and Caimans (Crocodylia, Alligatoridae) after a Half-Century Delay: Bridging the Gap in the Chromosomal Evolution of Reptiles

Although crocodilians have attracted enormous attention in other research fields, from the cytogenetic point of view, this group remains understudied. Here, we analyzed the karyotypes of eight species formally described from the Alligatoridae family using differential staining, fluorescence in situ hybridization with rDNA and repetitive motifs as a probe, whole chromosome painting (WCP), and comparative genome hybridization. All Caimaninae species have a diploid chromosome number (2n) 42 and karyotypes dominated by acrocentric chromosomes, in contrast to both species of Alligatorinae, which have 2n = 32 and karyotypes that are predominantly metacentric, suggesting fusion/fission rearrangements. Our WCP results supported this scenario by revealing the homeology of the largest metacentric pair present in both Alligator spp. with two smaller pairs of acrocentrics in Caimaninae species. The clusters of 18S rDNA were found on one chromosome pair in all species, except for Paleosuchus spp., which possessed three chromosome pairs bearing these sites. Similarly, comparative genomic hybridization demonstrated an advanced stage of sequence divergence among the caiman genomes, with Paleosuchus standing out as the most divergent. Thus, although Alligatoridae exhibited rather low species diversity and some level of karyotype stasis, their genomic content indicates that they are not as conserved as previously thought. These new data deepen the discussion of cytotaxonomy in this family.

wQFM: Highly Accurate Genome-scale Species Tree Estimation from Weighted Quartets

MOTIVATION

Species tree estimation from genes sampled from throughout the whole genome is complicated due to the gene tree-species tree discordance. Incomplete lineage sorting (ILS) is one of the most frequent causes for this discordance, where alleles can coexist in populations for periods that may span several speciation events. Quartet-based summary methods for estimating species trees from a collection of gene trees are becoming popular due to their high accuracy and statistical guarantee under ILS. Generating quartets with appropriate weights, where weights correspond to the relative importance of quartets, and subsequently amalgamating the weighted quartets to infer a single coherent species tree can allow for a statistically consistent way of estimating species trees. However, handling weighted quartets is challenging.

RESULTS

We propose wQFM, a highly accurate method for species tree estimation from multi-locus data, by extending the quartet FM (QFM) algorithm to a weighted setting. wQFM was assessed on a collection of simulated and real biological datasets, including the avian phylogenomic dataset which is one of the largest phylogenomic datasets to date. We compared wQFM with wQMC, which is the best alternate method for weighted quartet amalgamation, and with ASTRAL, which is one of the most accurate and widely used coalescent-based species tree estimation methods. Our results suggest that wQFM matches or improves upon the accuracy of wQMC and ASTRAL.

AVAILABILITY

wQFM is available in open source form at https://github.com/Mahim1997/wQFM-2020.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phylogeny and evolutionary history of the amniote egg

We review morphological features of the amniote egg and embryos in a comparative phylogenetic framework, including all major clades of extant vertebrates. We discuss 40 characters that are relevant for an analysis of the evolutionary history of the vertebrate egg. Special attention is given to the morphology of the cellular yolk sac, the eggshell, and extraembryonic membranes. Many features that are typically assigned to amniotes, such as a large yolk sac, delayed egg deposition, and terrestrial reproduction have evolved independently and convergently in numerous clades of vertebrates. We use phylogenetic character mapping and ancestral character state reconstruction as tools to recognize sequence, order, and patterns of morphological evolution and deduce a hypothesis of the evolutionary history of the amniote egg. Besides amnion and chorioallantois, amniotes ancestrally possess copulatory organs (secondarily reduced in most birds), internal fertilization, and delayed deposition of eggs that contain an embryo in the primitive streak or early somite stage. Except for the amnion, chorioallantois, and amniote type of eggshell, these features evolved convergently in almost all major clades of aquatic vertebrates possibly in response to selective factors such as egg predation, hostile environmental conditions for egg development, or to adjust hatching of young to favorable season. A functionally important feature of the amnion membrane is its myogenic contractility that moves the (early) embryo and prevents adhering of the growing embryo to extraembryonic materials. This function of the amnion membrane and the liquid-filled amnion cavity may have evolved under the requirements of delayed deposition of eggs that contain developing embryos. The chorioallantois is a temporary embryonic exchange organ that supports embryonic development. A possible evolutionary scenario is that the amniote egg presents an exaptation that paved the evolutionary pathway for reproduction on land. As shown by numerous examples from anamniotes, reproduction on land has occurred multiple times among vertebrates-the amniote egg presenting one "solution" that enabled the conquest of land for reproduction.

Origin and Evolution of the Turtle Body Plan

The origin of turtles and their uniquely shelled body plan is one of the longest standing problems in vertebrate biology. The unfulfilled need for a hypothesis that both explains the derived nature of turtle anatomy and resolves their unclear phylogenetic position among reptiles largely reflects the absence of a transitional fossil record. Recent discoveries have dramatically improved this situation, providing an integrated, time-calibrated model of the morphological, developmental, and ecological transformations responsible for the modern turtle body plan. This evolutionary trajectory was initiated in the Permian (>260 million years ago) when a turtle ancestor with a diapsid skull evolved a novel mechanism for lung ventilation. This key innovation permitted the torso to become apomorphically stiff, most likely as an adaption for digging and a fossorial ecology. The construction of the modern turtle body plan then proceeded over the next 100 million years following a largely stepwise model of osteological innovation.

Functional morphology, diversity, and evolution of yolk processing specializations in embryonic reptiles and birds

Evolution of the terrestrial, amniotic egg of vertebrates required new mechanisms by which yolk material could be processed for embryonic use. Recent studies on each of the major extant reptile groups have revealed elaborate morphological specializations for yolk processing, features that differ dramatically from those of birds. In the avian pattern, liquid yolk is housed in a yolk sac whose endodermal lining absorbs and digests yolk material and sends resultant nutrients into the blood circulation. In snakes, lizards, turtles, and crocodilians, as documented herein, the yolk sac becomes invaded by endodermal cells that proliferate and phagocytose yolk material. Blood vessels then invade, and the endodermal cells become arranged around them, forming elongated "spaghetti-like" strands that fill the yolk sac cavity. This pattern provides an effective means by which yolk material is cellularized, digested, and transported by vitelline vessels to the developing embryo. Phylogenetically, the (non-avian) "reptilian" pattern was ancestral for sauropsids and was modified or replaced in ancestors to birds. This review postulates that evolution of the "avian" pattern involved increased reliance on extracellular digestion of yolk, allowing embryonic development to occur more rapidly than in typical reptiles. Comparative studies of yolk processing that draw on morphological, biochemical, molecular approaches are needed to explain how and why the "reptilian" pattern was replaced in birds or their archosaurian ancestors.

Distribution of GFAP in Squamata: Extended Immunonegative Areas, Astrocytes, High Diversity, and Their Bearing on Evolution

Squamata is one of the richest and most diverse extant groups. The present study investigates the glial fibrillary acidic protein (GFAP)-immunopositive elements of five lizard and three snake species; each represents a different family. The study continues our former studies on bird, turtle, and caiman brains. Although several studies have been published on lizards, they usually only investigated one species. Almost no data are available on snakes. The animals were transcardially perfused. Immunoperoxidase reactions were performed with a mouse monoclonal anti-GFAP (Novocastra). The original radial ependymoglia is enmeshed by secondary, non-radial processes almost beyond recognition in several brain areas like in other reptiles. Astrocytes occur but only as complementary elements like in caiman but unlike in turtles, where astrocytes are absent. In most species, extended areas are free of GFAP—a meaningful difference from other reptiles. The predominance of astrocytes and the presence of areas free of GFAP immunopositivity are characteristic of birds and mammals; therefore, they must be apomorphic features of Squamata, which appeared independently from the evolution of avian glia. However, these features show a high diversity; in some lizards, they are even absent. There was no principal difference between the glial structures of snakes and lizards. In conclusion, the glial structure of Squamata seems to be the most apomorphic one among reptiles. The high diversity suggests that its evolution is still intense. The comparison of identical brain areas with different GFAP contents in different species may promote understanding the role of GFAP in neuronal networks. Our findings are in accordance with the supposal based on our previous studies that the GFAP-free areas expand during evolution.

Polarized Sonic Hedgehog Protein Localization and a Shift in the Expression of Region-Specific Molecules Is Associated With the Secondary Palate Development in the Veiled Chameleon

Secondary palate development is characterized by the formation of two palatal shelves on the maxillary prominences, which fuse in the midline in mammalian embryos. However, in reptilian species, such as turtles, crocodilians, and lizards, the palatal shelves of the secondary palate develop to a variable extent and morphology. While in most Squamates, the palate is widely open, crocodilians develop a fully closed secondary palate. Here, we analyzed developmental processes that underlie secondary palate formation in chameleons, where large palatal shelves extend horizontally toward the midline. The growth of the palatal shelves continued during post-hatching stages and closure of the secondary palate can be observed in several adult animals. The massive proliferation of a multilayered oral epithelium and mesenchymal cells in the dorsal part of the palatal shelves underlined the initiation of their horizontal outgrowth, and was decreased later in development. The polarized cellular localization of primary cilia and Sonic hedgehog protein was associated with horizontal growth of the palatal shelves. Moreover, the development of large palatal shelves, supported by the pterygoid and palatine bones, was coupled with the shift in Meox2, Msx1, and Pax9 gene expression along the rostro-caudal axis. In conclusion, our results revealed distinctive developmental processes that contribute to the expansion and closure of the secondary palate in chameleons and highlighted divergences in palate formation across amniote species.

STELAR: a statistically consistent coalescent-based species tree estimation method by maximizing triplet consistency

Background

Species tree estimation is frequently based on phylogenomic approaches that use multiple genes from throughout the genome. However, estimating a species tree from a collection of gene trees can be complicated due to the presence of gene tree incongruence resulting from incomplete lineage sorting (ILS), which is modelled by the multi-species coalescent process. Maximum likelihood and Bayesian MCMC methods can potentially result in accurate trees, but they do not scale well to large datasets.

Results

We present STELAR (Species Tree Estimation by maximizing tripLet AgReement), a new fast and highly accurate statistically consistent coalescent-based method for estimating species trees from a collection of gene trees. We formalized the constrained triplet consensus (CTC) problem and showed that the solution to the CTC problem is a statistically consistent estimate of the species tree under the multi-species coalescent (MSC) model. STELAR is an efficient dynamic programming based solution to the CTC problem which is highly accurate and scalable. We evaluated the accuracy of STELAR in comparison with SuperTriplets, which is an alternate fast and highly accurate triplet-based supertree method, and with MP-EST and ASTRAL – two of the most popular and accurate coalescent-based methods. Experimental results suggest that STELAR matches the accuracy of ASTRAL and improves on MP-EST and SuperTriplets.

Conclusions

Theoretical and empirical results (on both simulated and real biological datasets) suggest that STELAR is a valuable technique for species tree estimation from gene tree distributions.

How do embryonic turtles process yolk? Evidence from the Snapping Turtle, Chelydra serpentina (Chelydridae)

Compared with amphibians, oviparous reptiles and birds lay large eggs that contain abundant yolk. Because the yolk is extracellular, it must be taken up by cells of the yolk sac and metabolized so that products of yolk digestion can be transported to the embryo to fuel development. In birds, the yolk is processed by cells that line the inside of the yolk sac. A very different developmental pattern recently has been demonstrated in lizards and snakes in which the yolk sac cavity is converted to a compact mass of blood vessels lined by endodermal cells. In this study, we used scanning electron microscopy to determine which of these developmental patterns (if either) occurs in a representative chelonian, the North American Snapping Turtle (Chelydra serpentina (Linnaeus, 1758)). Our observations reveal that yolk-filled endodermal cells progressively fill the yolk sac cavity. These cells become organized around anastomosing blood vessels, forming elongated strands that are morphologically well suited for yolk digestion and vascular transport of nutrients. This developmental pattern shares features with that of squamates, but it differs markedly from that of birds. These observations indicate that mechanisms of yolk processing in lizards and snakes have an evolutionary history that pre-dates the squamate clade.

Asymmetric paralog evolution between the "cryptic" gene Bmp16 and its well-studied sister genes Bmp2 and Bmp4

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Quantitative histological models suggest endothermy in plesiosaurs

Background

Plesiosaurs are marine reptiles that arose in the Late Triassic and survived to the Late Cretaceous. They have a unique and uniform bauplan and are known for their very long neck and hydrofoil-like flippers. Plesiosaurs are among the most successful vertebrate clades in Earth’s history. Based on bone mass decrease and cosmopolitan distribution, both of which affect lifestyle, indications of parental care, and oxygen isotope analyses, evidence for endothermy in plesiosaurs has accumulated. Recent bone histological investigations also provide evidence of fast growth and elevated metabolic rates. However, quantitative estimations of metabolic rates and bone growth rates in plesiosaurs have not been attempted before.

Methods

Phylogenetic eigenvector maps is a method for estimating trait values from a predictor variable while taking into account phylogenetic relationships. As predictor variable, this study employs vascular density, measured in bone histological sections of fossil eosauropterygians and extant comparative taxa. We quantified vascular density as primary osteon density, thus, the proportion of vascular area (including lamellar infillings of primary osteons) to total bone area. Our response variables are bone growth rate (expressed as local bone apposition rate) and resting metabolic rate (RMR).

Results

Our models reveal bone growth rates and RMRs for plesiosaurs that are in the range of birds, suggesting that plesiosaurs were endotherm. Even for basal eosauropterygians we estimate values in the range of mammals or higher.

Discussion

Our models are influenced by the availability of comparative data, which are lacking for large marine amniotes, potentially skewing our results. However, our statistically robust inference of fast growth and fast metabolism is in accordance with other evidence for plesiosaurian endothermy. Endothermy may explain the success of plesiosaurs consisting in their survival of the end-Triassic extinction event and their global radiation and dispersal.

Pangenomic Definition of Prokaryotic Species and the Phylogenetic Structure of Prochlorococcus spp

The pangenome is the collection of all groups of orthologous genes (OGGs) from a set of genomes. We apply the pangenome analysis to propose a definition of prokaryotic species based on identification of lineage-specific gene sets. While being similar to the classical biological definition based on allele flow, it does not rely on DNA similarity levels and does not require analysis of homologous recombination. Hence this definition is relatively objective and independent of arbitrary thresholds. A systematic analysis of 110 accepted species with the largest numbers of sequenced strains yields results largely consistent with the existing nomenclature. However, it has revealed that abundant marine cyanobacteria Prochlorococcus marinus should be divided into two species. As a control we have confirmed the paraphyletic origin of Yersinia pseudotuberculosis (with embedded, monophyletic Y. pestis) and Burkholderia pseudomallei (with B. mallei). We also demonstrate that by our definition and in accordance with recent studies Escherichia coli and Shigella spp. are one species.

Gene Expression Profiling in the Injured Spinal Cord of Trachemys scripta elegans: An Amniote with Self-Repair Capabilities

Slider turtles are the only known amniotes with self-repair mechanisms of the spinal cord that lead to substantial functional recovery. Their strategic phylogenetic position makes them a relevant model to investigate the peculiar genetic programs that allow anatomical reconnection in some vertebrate groups but are absent in others. Here, we analyze the gene expression profile of the response to spinal cord injury (SCI) in the turtle Trachemys scripta elegans. We found that this response comprises more than 1000 genes affecting diverse functions: reaction to ischemic insult, extracellular matrix re-organization, cell proliferation and death, immune response, and inflammation. Genes related to synapses and cholesterol biosynthesis are down-regulated. The analysis of the evolutionary distribution of these genes shows that almost all are present in most vertebrates. Additionally, we failed to find genes that were exclusive of regenerating taxa. The comparison of expression patterns among species shows that the response to SCI in the turtle is more similar to that of mice and non-regenerative Xenopus than to Xenopus during its regenerative stage. This observation, along with the lack of conserved “regeneration genes” and the current accepted phylogenetic placement of turtles (sister group of crocodilians and birds), indicates that the ability of spinal cord self-repair of turtles does not represent the retention of an ancestral vertebrate character. Instead, our results suggest that turtles developed this capability from a non-regenerative ancestor (i.e., a lineage specific innovation) that was achieved by re-organizing gene expression patterns on an essentially non-regenerative genetic background. Among the genes activated by SCI exclusively in turtles, those related to anoxia tolerance, extracellular matrix remodeling, and axonal regrowth are good candidates to underlie functional recovery.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

Development of the turtle plastron, the order-defining skeletal structure

The dorsal and ventral aspects of the turtle shell, the carapace and the plastron, are developmentally different entities. The carapace contains axial endochondral skeletal elements and exoskeletal dermal bones. The exoskeletal plastron is found in all extant and extinct species of crown turtles found to date and is synaptomorphic of the order Testudines. However, paleontological reconstructed transition forms lack a fully developed carapace and show a progression of bony elements ancestral to the plastron. To understand the evolutionary development of the plastron, it is essential to know how it has formed. Here we studied the molecular development and patterning of plastron bones in a cryptodire turtle Trachemys scripta We show that plastron development begins at developmental stage 15 when osteochondrogenic mesenchyme forms condensates for each plastron bone at the lateral edges of the ventral mesenchyme. These condensations commit to an osteogenic identity and suppress chondrogenesis. Their development overlaps with that of sternal cartilage development in chicks and mice. Thus, we suggest that in turtles, the sternal morphogenesis is prevented in the ventral mesenchyme by the concomitant induction of osteogenesis and the suppression of chondrogenesis. The osteogenic subroutines later direct the growth and patterning of plastron bones in an autonomous manner. The initiation of plastron bone development coincides with that of carapacial ridge formation, suggesting that the development of dorsal and ventral shells are coordinated from the start and that adopting an osteogenesis-inducing and chondrogenesis-suppressing cell fate in the ventral mesenchyme has permitted turtles to develop their order-specific ventral morphology.

Palaeohistological Evidence for Ancestral High Metabolic Rate in Archosaurs

Metabolic heat production in archosaurs has played an important role in their evolutionary radiation during the Mesozoic, and their ancestral metabolic condition has long been a matter of debate in systematics and palaeontology. The study of fossil bone histology provides crucial information on bone growth rate, which has been used to indirectly investigate the evolution of thermometabolism in archosaurs. However, no quantitative estimation of metabolic rate has ever been performed on fossils using bone histological features. Moreover, to date, no inference model has included phylogenetic information in the form of predictive variables. Here we performed statistical predictive modeling using the new method of phylogenetic eigenvector maps on a set of bone histological features for a sample of extant and extinct vertebrates, to estimate metabolic rates of fossil archosauromorphs. This modeling procedure serves as a case study for eigenvector-based predictive modeling in a phylogenetic context, as well as an investigation of the poorly known evolutionary patterns of metabolic rate in archosaurs. Our results show that Mesozoic theropod dinosaurs exhibit metabolic rates very close to those found in modern birds, that archosaurs share a higher ancestral metabolic rate than that of extant ectotherms, and that this derived high metabolic rate was acquired at a much more inclusive level of the phylogenetic tree, among non-archosaurian archosauromorphs. These results also highlight the difficulties of assigning a given heat production strategy (i.e., endothermy, ectothermy) to an estimated metabolic rate value, and confirm findings of previous studies that the definition of the endotherm/ectotherm dichotomy may be ambiguous.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Molecular characterization of insulin from squamate reptiles reveals sequence diversity and possible adaptive evolution

The Squamata are the most adaptive and prosperous group among ectothermic amniotes, reptiles, due to their species-richness and geographically wide habitat. Although the molecular mechanisms underlying their prosperity remain largely unknown, unique features have been reported from hormones that regulate energy metabolism. Insulin, a central anabolic hormone, is one such hormone, as its roles and effectiveness in regulation of blood glucose levels remain to be examined in squamates. In the present study, cDNAs coding for insulin were isolated from multiple species that represent various groups of squamates. The deduced amino acid sequences showed a high degree of divergence, with four lineages showing obviously higher number of amino acid substitutions than most of vertebrates, from teleosts to mammals. Among 18 sites presented to comprise the two receptor binding surfaces (one with 12 sites and the other with 6 sites), substitutions were observed in 13 sites. Among them was the substitution of HisB10, which results in the loss of the ability to hexamerize. Furthermore, three of these substitutions were reported to increase mitogenicity in human analogues. These substitutions were also reported from insulin of hystricomorph rodents and agnathan fishes, whose mitogenic potency have been shown to be increased. The estimated value of the non-synonymous-to-synonymous substitution ratio (ω) for the Squamata clade was larger than those of the other reptiles and aves. Even higher values were estimated for several lineages among squamates. These results, together with the regulatory mechanisms of digestion and nutrient assimilation in squamates, suggested a possible adaptive process through the molecular evolution of squamate INS. Further studies on the roles of insulin, in relation to the physiological and ecological traits of squamate species, will provide an insight into the molecular mechanisms that have led to the adaptivity and prosperity of squamates.

Recent insights into the morphological diversity in the amniote primary and secondary palates

The assembly of the upper jaw is a pivotal moment in the embryonic development of amniotes. The upper jaw forms from the fusion of the maxillary, medial nasal, and lateral nasal prominences, resulting in an intact upper lip/beak and nasal cavities; together called the primary palate. This process of fusion requires a balance of proper facial prominence shape and positioning to avoid craniofacial clefting, whilst still accommodating the vast phenotypic diversity of adult amniotes. As such, variation in craniofacial ontogeny is not tolerated beyond certain bounds. For clarity, we discuss primary palatogenesis of amniotes into in two categories, according to whether the nasal and oral cavities remain connected throughout ontogeny or not. The transient separation of these cavities occurs in mammals and crocodilians, while remaining connected in birds, turtles and squamates. In the latter group, the craniofacial prominences fuse around a persistent choanal groove that connects the nasal and oral cavities. Subsequently, all lineages except for turtles, develop a secondary palate that ultimately completely or partially separates oral and nasal cavities. Here, we review the shared, early developmental events and highlight the points at which development diverges in both primary and secondary palate formation.

Comparative Genomics Identifies Epidermal Proteins Associated with the Evolution of the Turtle Shell

The evolution of reptiles, birds, and mammals was associated with the origin of unique integumentary structures. Studies on lizards, chicken, and humans have suggested that the evolution of major structural proteins of the outermost, cornified layers of the epidermis was driven by the diversification of a gene cluster called Epidermal Differentiation Complex (EDC). Turtles have evolved unique defense mechanisms that depend on mechanically resilient modifications of the epidermis. To investigate whether the evolution of the integument in these reptiles was associated with specific adaptations of the sequences and expression patterns of EDC-related genes, we utilized newly available genome sequences to determine the epidermal differentiation gene complement of turtles. The EDC of the western painted turtle (Chrysemys picta bellii) comprises more than 100 genes, including at least 48 genes that encode proteins referred to as beta-keratins or corneous beta-proteins. Several EDC proteins have evolved cysteine/proline contents beyond 50% of total amino acid residues. Comparative genomics suggests that distinct subfamilies of EDC genes have been expanded and partly translocated to loci outside of the EDC in turtles. Gene expression analysis in the European pond turtle (Emys orbicularis) showed that EDC genes are differentially expressed in the skin of the various body sites and that a subset of beta-keratin genes within the EDC as well as those located outside of the EDC are expressed predominantly in the shell. Our findings give strong support to the hypothesis that the evolutionary innovation of the turtle shell involved specific molecular adaptations of epidermal differentiation.

Sound localization in the alligator

In early tetrapods, it is assumed that the tympana were acoustically coupled through the pharynx and therefore inherently directional, acting as pressure difference receivers. The later closure of the middle ear cavity in turtles, archosaurs, and mammals is a derived condition, and would have changed the ear by decoupling the tympana. Isolation of the middle ears would then have led to selection for structural and neural strategies to compute sound source localization in both archosaurs and mammalian ancestors. In the archosaurs (birds and crocodilians) the presence of air spaces in the skull provided connections between the ears that have been exploited to improve directional hearing, while neural circuits mediating sound localization are well developed. In this review, we will focus primarily on directional hearing in crocodilians, where vocalization and sound localization are thought to be ecologically important, and indicate important issues still awaiting resolution.

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.

Disk covering methods improve phylogenomic analyses

Motivation

With the rapid growth rate of newly sequenced genomes, species tree inference from multiple genes has become a basic bioinformatics task in comparative and evolutionary biology. However, accurate species tree estimation is difficult in the presence of gene tree discordance, which is often due to incomplete lineage sorting (ILS), modelled by the multi-species coalescent. Several highly accurate coalescent-based species tree estimation methods have been developed over the last decade, including MP-EST. However, the running time for MP-EST increases rapidly as the number of species grows.

Results

We present divide-and-conquer techniques that improve the scalability of MP-EST so that it can run efficiently on large datasets. Surprisingly, this technique also improves the accuracy of species trees estimated by MP-EST, as our study shows on a collection of simulated and biological datasets.

Evaluating Summary Methods for Multilocus Species Tree Estimation in the Presence of Incomplete Lineage Sorting

Species tree estimation is complicated by processes, such as gene duplication and loss and incomplete lineage sorting (ILS), that cause discordance between gene trees and the species tree. Furthermore, while concatenation, a traditional approach to tree estimation, has excellent performance under many conditions, the expectation is that the best accuracy will be obtained through the use of species tree estimation methods that are specifically designed to address gene tree discordance. In this article, we report on a study to evaluate MP-EST-one of the most popular species tree estimation methods designed to address ILS-as well as concatenation under maximum likelihood, the greedy consensus, and two supertree methods (Matrix Representation with Parsimony and Matrix Representation with Likelihood). Our study shows that several factors impact the absolute and relative accuracy of methods, including the number of gene trees, the accuracy of the estimated gene trees, and the amount of ILS. Concatenation can be more accurate than the best summary methods in some cases (mostly when the gene trees have poor phylogenetic signal or when the level of ILS is low), but summary methods are generally more accurate than concatenation when there are an adequate number of sufficiently accurate gene trees. Our study suggests that coalescent-based species tree methods may be key to estimating highly accurate species trees from multiple loci.

Toward consilience in reptile phylogeny: miRNAs support an archosaur, not lepidosaur, affinity for turtles

Understanding the phylogenetic position of crown turtles (Testudines) among amniotes has been a source of particular contention. Recent morphological analyses suggest that turtles are sister to all other reptiles, whereas the vast majority of gene sequence analyses support turtles as being inside Diapsida, and usually as sister to crown Archosauria (birds and crocodilians). Previously, a study using microRNAs (miRNAs) placed turtles inside diapsids, but as sister to lepidosaurs (lizards and Sphenodon) rather than archosaurs. Here, we test this hypothesis with an expanded miRNA presence/absence dataset, and employ more rigorous criteria for miRNA annotation. Significantly, we find no support for a turtle + lepidosaur sister-relationship; instead, we recover strong support for turtles sharing a more recent common ancestor with archosaurs. We further test this result by analyzing a super-alignment of precursor miRNA sequences for every miRNA inferred to have been present in the most recent common ancestor of tetrapods. This analysis yields a topology that is fully congruent with our presence/absence analysis; our results are therefore in accordance with most gene sequence studies, providing strong, consilient molecular evidence from diverse independent datasets regarding the phylogenetic position of turtles.

Evolution of the vertebrate Pax4/6 class of genes with focus on its novel member, the Pax10 gene

The members of the paired box (Pax) family regulate key developmental pathways in many metazoans as tissue-specific transcription factors. Vertebrate genomes typically possess nine Pax genes (Pax1-9), which are derived from four proto-Pax genes in the vertebrate ancestor that were later expanded through the so-called two-round (2R) whole-genome duplication. A recent study proposed that pax6a genes of a subset of teleost fishes (namely, acanthopterygians) are remnants of a paralog generated in the 2R genome duplication, to be renamed pax6.3, and reported one more group of vertebrate Pax genes (Pax6.2), most closely related to the Pax4/6 class. We propose to designate this new member Pax10 instead and reconstruct the evolutionary history of the Pax4/6/10 class with solid phylogenetic evidence. Our synteny analysis showed that Pax4, -6, and -10 originated in the 2R genome duplications early in vertebrate evolution. The phylogenetic analyses of relationships between teleost pax6a and other Pax4, -6, and -10 genes, however, do not support the proposed hypothesis of an ancient origin of the acanthopterygian pax6a genes in the 2R genome duplication. Instead, we confirmed the traditional scenario that the acanthopterygian pax6a is derived from the more recent teleost-specific genome duplication. Notably, Pax6 is present in all vertebrates surveyed to date, whereas Pax4 and -10 were lost multiple times in independent vertebrate lineages, likely because of their restricted expression patterns: Among Pax6-positive domains, Pax10 has retained expression in the adult retina alone, which we documented through in situ hybridization and quantitative reverse transcription polymerase chain reaction experiments on zebrafish, Xenopus, and anole lizard.

Antimicrobial peptides in reptiles

Reptiles are among the oldest known amniotes and are highly diverse in their morphology and ecological niches. These animals have an evolutionarily ancient innate-immune system that is of great interest to scientists trying to identify new and useful antimicrobial peptides. Significant work in the last decade in the fields of biochemistry, proteomics and genomics has begun to reveal the complexity of reptilian antimicrobial peptides. Here, the current knowledge about antimicrobial peptides in reptiles is reviewed, with specific examples in each of the four orders: Testudines (turtles and tortosises), Sphenodontia (tuataras), Squamata (snakes and lizards), and Crocodilia (crocodilans). Examples are presented of the major classes of antimicrobial peptides expressed by reptiles including defensins, cathelicidins, liver-expressed peptides (hepcidin and LEAP-2), lysozyme, crotamine, and others. Some of these peptides have been identified and tested for their antibacterial or antiviral activity; others are only predicted as possible genes from genomic sequencing. Bioinformatic analysis of the reptile genomes is presented, revealing many predicted candidate antimicrobial peptides genes across this diverse class. The study of how these ancient creatures use antimicrobial peptides within their innate immune systems may reveal new understandings of our mammalian innate immune system and may also provide new and powerful antimicrobial peptides as scaffolds for potential therapeutic development.

The evolutionary origin of the turtle shell and its dependence on the axial arrest of the embryonic rib cage

Turtles are characterized by their possession of a shell with dorsal and ventral moieties: the carapace and the plastron, respectively. In this review, we try to provide answers to the question of the evolutionary origin of the carapace, by revising morphological, developmental, and paleontological comparative analyses. The turtle carapace is formed through modification of the thoracic ribs and vertebrae, which undergo extensive ossification to form a solid bony structure. Except for peripheral dermal elements, there are no signs of exoskeletal components ontogenetically added to the costal and neural bones, and thus the carapace is predominantly of endoskeletal nature. Due to the axial arrest of turtle rib growth, the axial part of the embryo expands laterally and the shoulder girdle becomes encapsulated in the rib cage, together with the inward folding of the lateral body wall in the late phase of embryogenesis. Along the line of this folding develops a ridge called the carapacial ridge (CR), a turtle-specific embryonic structure. The CR functions in the marginal growth of the carapacial primordium, in which Wnt signaling pathway might play a crucial role. Both paleontological and genomic evidence suggest that the axial arrest is the first step toward acquisition of the turtle body plan, which is estimated to have taken place after the divergence of a clade including turtles from archosaurs. The developmental relationship between the CR and the axial arrest remains a central issue to be solved in future.

Unidirectional photoreceptor-to-Müller glia coupling and unique K+ channel expression in Caiman retina

Background

Müller cells, the principal glial cells of the vertebrate retina, are fundamental for the maintenance and function of neuronal cells. In most vertebrates, including humans, Müller cells abundantly express Kir4.1 inwardly rectifying potassium channels responsible for hyperpolarized membrane potential and for various vital functions such as potassium buffering and glutamate clearance; inter-species differences in Kir4.1 expression were, however, observed. Localization and function of potassium channels in Müller cells from the retina of crocodiles remain, hitherto, unknown.

Methods

We studied retinae of the Spectacled caiman (Caiman crocodilus fuscus), endowed with both diurnal and nocturnal vision, by (i) immunohistochemistry, (ii) whole-cell voltage-clamp, and (iii) fluorescent dye tracing to investigate K⁺ channel distribution and glia-to-neuron communications.

Results

Immunohistochemistry revealed that caiman Müller cells, similarly to other vertebrates, express vimentin, GFAP, S100β, and glutamine synthetase. In contrast, Kir4.1 channel protein was not found in Müller cells but was localized in photoreceptor cells. Instead, 2P-domain TASK-1 channels were expressed in Müller cells. Electrophysiological properties of enzymatically dissociated Müller cells without photoreceptors and isolated Müller cells with adhering photoreceptors were significantly different. This suggests ion coupling between Müller cells and photoreceptors in the caiman retina. Sulforhodamine-B injected into cones permeated to adhering Müller cells thus revealing a uni-directional dye coupling.

Conclusion

Our data indicate that caiman Müller glial cells are unique among vertebrates studied so far by predominantly expressing TASK-1 rather than Kir4.1 K⁺ channels and by bi-directional ion and uni-directional dye coupling to photoreceptor cells. This coupling may play an important role in specific glia-neuron signaling pathways and in a new type of K⁺ buffering.

A common developmental plan for neocortical gene-expressing neurons in the pallium of the domestic chicken Gallus gallus domesticus and the Chinese softshell turtle Pelodiscus sinensis

The six-layered neocortex is a unique characteristic of mammals and likely provides the neural basis of their sophisticated cognitive abilities. Although all mammalian species share the layered structure of the neocortex, the sauropsids exhibit an entirely different cytoarchitecture of the corresponding pallial region. Our previous gene expression study revealed that the chicken pallium possesses neural subtypes that express orthologs of layer-specific genes of the mammalian neocortex. To understand the evolutionary steps leading toward animal group-specific neuronal arrangements in the pallium in the course of amniote diversification, we examined expression patterns of the same orthologs and a few additional genes in the pallial development of the Chinese softshell turtle Pelodiscus sinensis, and compared these patterns to those of the chicken. Our analyses highlighted similarities in neuronal arrangements between the two species; the mammalian layer 5 marker orthologs are expressed in the medial domain and the layer 2/3 marker orthologs are expressed in the lateral domain in the pallia of both species. We hypothesize that the mediolateral arrangement of the neocortical layer-specific gene-expressing neurons originated in their common ancestor and is conserved among all sauropsid groups, whereas the neuronal arrangement within the pallium could have highly diversified independently in the mammalian lineage.

The endoskeletal origin of the turtle carapace

The turtle body plan, with its solid shell, deviates radically from those of other tetrapods. The dorsal part of the turtle shell, or the carapace, consists mainly of costal and neural bony plates, which are continuous with the underlying thoracic ribs and vertebrae, respectively. Because of their superficial position, the evolutionary origins of these costo-neural elements have long remained elusive. Here we show, through comparative morphological and embryological analyses, that the major part of the carapace is derived purely from endoskeletal ribs. We examine turtle embryos and find that the costal and neural plates develop not within the dermis, but within deeper connective tissue where the rib and intercostal muscle anlagen develop. We also examine the fossils of an outgroup of turtles to confirm that the structure equivalent to the turtle carapace developed independently of the true osteoderm. Our results highlight the hitherto unravelled evolutionary course of the turtle shell.

The evolutionary origins of the costal and neural bony plates of the turtle shell have long remained elusive. Here the authors show, through comparative morphological and embryological analyses, that the most of the carapace is derived from endoskeletal ribs.

Divergent palate morphology in turtles and birds correlates with differences in proliferation and BMP2 expression during embryonic development

During embryonic development, amniotes typically form outgrowths from the medial sides of the maxillary prominences called palatal shelves or palatine processes. In mammals the shelves fuse in the midline and form a bony hard palate that completely separates the nasal and oral cavities. In birds and lizards, palatine processes develop but remain unfused, leaving a natural cleft. Adult turtles do not possess palatine processes and unlike other amniotes, the internal nares open into the oral cavity. Here we investigate craniofacial ontogeny in the turtle, Emydura subglobosa to determine whether vestigial palatine processes develop and subsequently regress, or whether development fails entirely. We found that the primary palate in turtles develops similarly to other amniotes, but secondary palate ontogeny diverges. Using histology, cellular dynamics and in situ hybridization we found no evidence of palatine process development at any time during ontogeny of the face in the turtle. Furthermore, detailed comparisons with chicken embryos (the model organism most closely related to turtles from a molecular phylogeny perspective), we identified differences in proliferation and gene expression patterns that correlate with the differences in palate morphology. We propose that, in turtles, palatine process outgrowth is never initiated due to a lack of mesenchymal bone morphogenetic protein 2 (BMP2) expression in the maxillary mesenchyme, which in turn fails to induce the relatively higher cellular proliferation required for medial tissue outgrowth. It is likely that these differences between turtles and birds arose after the divergence of the lineage leading to modern turtles.

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.

New insights into non-avian dinosaur reproduction and their evolutionary and ecological implications: linking fossil evidence to allometries of extant close relatives

It has been hypothesized that a high reproductive output contributes to the unique gigantism in large dinosaur taxa. In order to infer more information on dinosaur reproduction, we established allometries between body mass and different reproductive traits (egg mass, clutch mass, annual clutch mass) for extant phylogenetic brackets (birds, crocodiles and tortoises) of extinct non-avian dinosaurs. Allometries were applied to nine non-avian dinosaur taxa (theropods, hadrosaurs, and sauropodomorphs) for which fossil estimates on relevant traits are currently available. We found that the reproductive traits of most dinosaurs conformed to similar-sized or scaled-up extant reptiles or birds. The reproductive traits of theropods, which are considered more bird-like, were indeed consistent with birds, while the traits of sauropodomorphs conformed better to reptiles. Reproductive traits of hadrosaurs corresponded to both reptiles and birds. Excluding Massospondyluscarinatus, all dinosaurs studied had an intermediary egg to body mass relationship to reptiles and birds. In contrast, dinosaur clutch masses fitted with either the masses predicted from allometries of birds (theropods) or to the masses of reptiles (all other taxa). Theropods studied had probably one clutch per year. For sauropodomorphs and hadrosaurs, more than one clutch per year was predicted. Contrary to current hypotheses, large dinosaurs did not have exceptionally high annual egg numbers (AEN). Independent of the extant model, the estimated dinosaur AEN did not exceed 850 eggs (75,000 kg sauropod) for any of the taxa studied. This estimated maximum is probably an overestimation due to unrealistic assumptions. According to most AEN estimations, the dinosaurs studied laid less than 200 eggs per year. Only some AEN estimates obtained for medium to large sized sauropods were higher (200-400 eggs). Our results provide new (testable) hypotheses, especially for reproductive traits that are insufficiently documented or lacking from the fossil record. This contributes to the understanding of their evolution.

Skeletal gene expression in the temporal region of the reptilian embryos: implications for the evolution of reptilian skull morphology

Reptiles have achieved highly diverse morphological and physiological traits that allow them to exploit various ecological niches and resources. Morphology of the temporal region of the reptilian skull is highly diverse and historically it has been treated as an important character for classifying reptiles and has helped us understand the ecology and physiology of each species. However, the developmental mechanism that generates diversity of reptilian skull morphology is poorly understood. We reveal a potential developmental basis that generates morphological diversity in the temporal region of the reptilian skull by performing a comparative analysis of gene expression in the embryos of reptile species with different skull morphology. By investigating genes known to regulate early osteoblast development, we find dorsoventrally broadened unique expression of the early osteoblast marker, Runx2, in the temporal region of the head of turtle embryos that do not form temporal fenestrae. We also observe that Msx2 is also uniquely expressed in the mesenchymal cells distributed at the temporal region of the head of turtle embryos. Furthermore, through comparison of gene expression pattern in the embryos of turtle, crocodile, and snake species, we find a possible correlation between the spatial patterns of Runx2 and Msx2 expression in cranial mesenchymal cells and skull morphology of each reptilian lineage. Regulatory modifications of Runx2 and Msx2 expression in osteogenic mesenchymal precursor cells are likely involved in generating morphological diversity in the temporal region of the reptilian skull.

Middle ear cavity morphology is consistent with an aquatic origin for testudines

The position of testudines in vertebrate phylogeny is being re-evaluated. At present, testudine morphological and molecular data conflict when reconstructing phylogenetic relationships. Complicating matters, the ecological niche of stem testudines is ambiguous. To understand how turtles have evolved to hear in different environments, we examined middle ear morphology and scaling in most extant families, as well as some extinct species, using 3-dimensional reconstructions from micro magnetic resonance (MR) and submillimeter computed tomography (CT) scans. All families of testudines exhibited a similar shape of the bony structure of the middle ear cavity, with the tympanic disk located on the rostrolateral edge of the cavity. Sea Turtles have additional soft tissue that fills the middle ear cavity to varying degrees. When the middle ear cavity is modeled as an air-filled sphere of the same volume resonating in an underwater sound field, the calculated resonances for the volumes of the middle ear cavities largely fell within testudine hearing ranges. Although there were some differences in morphology, there were no statistically significant differences in the scaling of the volume of the bony middle ear cavity with head size among groups when categorized by phylogeny and ecology. Because the cavity is predicted to resonate underwater within the testudine hearing range, the data support the hypothesis of an aquatic origin for testudines, and function of the middle ear cavity in underwater sound detection.