The phylogenetic trunk: maximal inclusion of taxa with missing data in an analysis of the lepospondyli (Vertebrata, Tetrapoda)

Returning to the roots: resolution, reproducibility, and robusticity in the phylogenetic inference of Dissorophidae (Amphibia: Temnospondyli)

The phylogenetic relationships of most Paleozoic tetrapod clades remain poorly resolved, which is variably attributed to a lack of study, the limitations of inference from phenotypic data, and constant revision of best practices. While refinement of phylogenetic methods continues to be important, any phylogenetic analysis is inherently constrained by the underlying dataset that it analyzes. Therefore, it becomes equally important to assess the accuracy of these datasets, especially when a select few are repeatedly propagated. While repeat analyses of these datasets may appear to constitute a working consensus, they are not in fact independent, and it becomes especially important to evaluate the accuracy of these datasets in order to assess whether a seeming consensus is robust. Here I address the phylogeny of the Dissorophidae, a speciose clade of Paleozoic temnospondyls. This group is an ideal case study among temnospondyls for exploring phylogenetic methods and datasets because it has been extensively studied (eight phylogenetic studies to date) but with most (six studies) using a single matrix that has been propagated with very little modification. In spite of the conserved nature of the matrix, dissorophid studies have produced anything but a conserved topology. Therefore, I analyzed an independently designed matrix, which recovered less resolution and some disparate nodes compared to previous studies. In order to reconcile these differences, I carefully examined previous matrices and analyses. While some differences are a matter of personal preference (e.g., analytical software), others relate to discrepancies with respect to what are currently considered as best practices. The most concerning discovery was the identification of pervasive dubious scorings that extend back to the origins of the widely propagated matrix. These include scores for skeletal features that are entirely unknown in a given taxon (e.g., postcrania in Cacops woehri) and characters for which there appear to be unstated working assumptions to scoring that are incompatible with the character definitions (e.g., scoring of taxa with incomplete skulls for characters based on skull length). Correction of these scores and other pervasive errors recovered a distinctly less resolved topology than previous studies, more in agreement with my own matrix. This suggests that previous analyses may have been compromised, and that the only real consensus of dissorophid phylogeny is the lack of one.

Revised Description of the Early Permian Recumbirostran “Microsaur” Nannaroter mckinziei Based on New Fossil Material and Computed Tomographic Data

The early Permian Richards Spur locality of Oklahoma has produced abundant material of numerous terrestrial fossil tetrapods, including various “microsaurs,” several of which are considered to belong to the clade Recumbirostra. We present a new partial skull of the recumbirostran “microsaur” Nannaroter mckinziei; through computed tomography (CT) analysis of both this new specimen and the holotype, we provide an updated description of the taxon. This new description provides novel information regarding several regions that could not be examined previously due to either being absent in the holotype or difficult to access. This includes missing and obscured aspects of the skull roof, braincase, lower jaw, and the palatal region. Furthermore, the new information obtained from this description was used to update phylogenetic character codings of Nannaroter, and a revised phylogenetic analysis was conducted. The results of this updated analysis are congruent with those of other recent phylogenetic analyses of recumbirostran “microsaurs.” This new information adds to the ever-growing body of early tetrapod CT data, which has been, and will continue to be, important in revealing details regarding early tetrapod anatomy, interrelationships, paleoecology, and evolution.

Size matters: the effects of ontogenetic disparity on the phylogeny of Trematopidae (Amphibia: Temnospondyli)

Trematopids are a clade of terrestrial Permo-Carboniferous temnospondyl amphibians. The intrarelationships of this clade are poorly known. This is largely attributable to a substantial disparity in size between type specimens, which range from the small-bodied lectotype of Mattauschia laticeps (< 4 cm skull length) to the large-bodied holotype of Acheloma cumminsi (> 15 cm skull length). Inferred correlation of size disparity with ontogenetic disparity has led previous workers either to omit taxa in phylogenetic analyses or to forgo an analysis altogether. Here, I take a specimen-level approach and multiple subsampling permutations to explore the phylogeny of the Trematopidae as a case study for assessing the effects of ontogenetic disparity on phylogenetic reconstruction in temnospondyls. The various analyses provide evidence that ontogenetic disparity confounds the phylogenetic inference of trematopids but without a directional bias. Tree topologies of most permutations are poorly resolved and weakly supported, reflecting character conflict that results from the inability of the analyses to differentiate retained plesiomorphies from juvenile features. These findings urge caution in the interpretation of phylogenetic analyses for which ontogenetic disparity exists, but is unaccounted for, and provide a strong impetus for more directed exploration of the interplay of ontogeny and phylogeny across Temnospondyli.

Phylogeny of Paleozoic limbed vertebrates reassessed through revision and expansion of the largest published relevant data matrix

The largest published phylogenetic analysis of early limbed vertebrates (Ruta M, Coates MI. 2007. Journal of Systematic Palaeontology 5:69-122) recovered, for example, Seymouriamorpha, Diadectomorpha and (in some trees) Caudata as paraphyletic and found the "temnospondyl hypothesis" on the origin of Lissamphibia (TH) to be more parsimonious than the "lepospondyl hypothesis" (LH)-though only, as we show, by one step. We report 4,200 misscored cells, over half of them due to typographic and similar accidental errors. Further, some characters were duplicated; some had only one described state; for one, most taxa were scored after presumed relatives. Even potentially continuous characters were unordered, the effects of ontogeny were not sufficiently taken into account, and data published after 2001 were mostly excluded. After these issues are improved-we document and justify all changes to the matrix-but no characters are added, we find (Analysis R1) much longer trees with, for example, monophyletic Caudata, Diadectomorpha and (in some trees) Seymouriamorpha; Ichthyostega either crownward or rootward of Acanthostega; and Anthracosauria either crownward or rootward of Temnospondyli. The LH is nine steps shorter than the TH (R2; constrained) and 12 steps shorter than the "polyphyly hypothesis" (PH-R3; constrained). Brachydectes (Lysorophia) is not found next to Lissamphibia; instead, a large clade that includes the adelogyrinids, urocordylid "nectrideans" and aïstopods occupies that position. As expected from the taxon/character ratio, most bootstrap values are low. Adding 56 terminal taxa to the original 102 increases the resolution (and decreases most bootstrap values). The added taxa range in completeness from complete articulated skeletons to an incomplete lower jaw. Even though the lissamphibian-like temnospondyls Gerobatrachus, Micropholis and Tungussogyrinus and the extremely peramorphic salamander Chelotriton are added, the difference between LH (R4; unconstrained) and TH (R5) rises to 10 steps, that between LH and PH (R6) to 15; the TH also requires several more regains of lost bones than the LH. Casineria, in which we tentatively identify a postbranchial lamina, emerges rather far from amniote origins in a gephyrostegid-chroniosuchian grade. Bayesian inference (Analysis EB, settings as in R4) mostly agrees with R4. High posterior probabilities are found for Lissamphibia (1.00) and the LH (0.92); however, many branches remain weakly supported, and most are short, as expected from the small character sample. We discuss phylogeny, approaches to coding, methods of phylogenetics (Bayesian inference vs. equally weighted vs. reweighted parsimony), some character complexes (e.g. preaxial/postaxial polarity in limb development), and prospects for further improvement of this matrix. Even in its revised state, the matrix cannot provide a robust assessment of the phylogeny of early limbed vertebrates. Sufficient improvement will be laborious-but not difficult.

Diversity change during the rise of tetrapods and the impact of the 'Carboniferous rainforest collapse'

The Carboniferous and early Permian were critical intervals in the diversification of early four-limbed vertebrates (tetrapods), yet the major patterns of diversity and biogeography during this time remain unresolved. Previous estimates suggest that global tetrapod diversity rose continuously across this interval and that habitat fragmentation following the 'Carboniferous rainforest collapse' (CRC) drove increased endemism among communities. However, previous work failed to adequately account for spatial and temporal biases in sampling. Here, we reassess early tetrapod diversity and biogeography with a new global species-level dataset using sampling standardization and network biogeography methods. Our results support a tight relationship between observed richness and sampling, particularly during the Carboniferous. We found that subsampled species richness initially increased into the late Carboniferous, then decreased substantially across the Carboniferous/Permian boundary before slowly recovering in the early Permian. Our analysis of biogeography does not support the hypothesis that the CRC drove endemism; instead, we found evidence for increased cosmopolitanism in the early Permian. While a changing environment may have played a role in reducing diversity in the earliest Permian, our results suggest that the CRC was followed by increased global connectivity between communities, possibly reflecting both reduced barriers to dispersal and the diversification of amniotes.

A primitive actinopterygian braincase from the Tournaisian of Nova Scotia

The vertebrate fossil record of the earliest Carboniferous is notoriously poorly sampled, obscuring a critical interval in vertebrate evolution and diversity. Recent studies of diversity across the Devonian-Carboniferous boundary have proposed a vertebrate mass extinction at the end-Devonian, and recent phylogenies suggest that the origin of the actinopterygian crown may have occurred in the earliest Carboniferous, as part of a broader recovery fauna. However, the data necessary to test this are limited. Here, we describe a partial actinopterygian skull, including diagnostic elements of the posterior braincase, from the Tournaisian Horton Bluff Formation of Blue Beach, Nova Scotia. The braincase surprisingly shows a confluence of characters common in Devonian taxa but absent in Mississippian forms, such as an open spiracular groove; lateral dorsal aortae that pass through open broadly separated, parallel grooves in the ventral otoccipital region, posterior to the articulation of the first infrapharyngobranchial and an intertemporal-supratemporal complex. Phylogenetic analysis places it deep within the actinopterygian stem, among Devonian moythomasiids and mimiids, suggesting more phylogenetically inclusive survivorship of stem group actinopterygians across the end-Devonian mass extinction. With a high lineage survivorship in tetrapods and lungfish across the Devonian-Carboniferous boundary and high vertebrate diversity at Tournaisian localities, this hints at a more gradual turnover between Devonian and Carboniferous vertebrate faunas.

Cranial Morphology of the Carboniferous-Permian Tetrapod Brachydectes newberryi (Lepospondyli, Lysorophia): New Data from µCT

Lysorophians are a group of early tetrapods with extremely elongate trunks, reduced limbs, and highly reduced skulls. Since the first discovery of this group, general similarities in outward appearance between lysorophians and some modern lissamphibian orders (specifically Urodela and Gymnophiona) have been recognized, and sometimes been the basis for hypotheses of lissamphibian origins. We studied the morphology of the skull, with particular emphasis on the neurocranium, of a partial growth series of the lysorophian Brachydectes newberryi using x-ray micro-computed tomography (μCT). Our study reveals similarities between the braincase of Brachydectes and brachystelechid recumbirostrans, corroborating prior work suggesting a close relationship between these taxa. We also describe the morphology of the epipterygoid, stapes, and quadrate in this taxon for the first time. Contra the proposals of some workers, we find no evidence of expected lissamphibian synapomorphies in the skull morphology in Brachydectes newberryi, and instead recognize a number of derived amniote characteristics within the braincase and suspensorium. Morphology previously considered indicative of taxonomic diversity within Lysorophia may reflect ontogenetic rather than taxonomic variation. The highly divergent morphology of lysorophians represents a refinement of morphological and functional trends within recumbirostrans, and is analogous to morphology observed in many modern fossorial reptiles.

Phylogenetic Stability, Tree Shape, and Character Compatibility: A Case Study Using Early Tetrapods

Phylogenetic tree shape varies as the evolutionary processes affecting a clade change over time. In this study, we examined an empirical phylogeny of fossil tetrapods during several time intervals, and studied how temporal constraints manifested in patterns of tree imbalance and character change. The results indicate that the impact of temporal constraints on tree shape is minimal and highlights the stability through time of the reference tetrapod phylogeny. Unexpected values of imbalance for Mississippian and Pennsylvanian time slices strongly support the hypothesis that the Carboniferous was a period of explosive tetrapod radiation. Several significant diversification shifts take place in the Mississippian and underpin increased terrestrialization among the earliest limbed vertebrates. Character incompatibility is relatively high at the beginning of tetrapod history, but quickly decreases to a relatively stable lower level, relative to a null distribution based on constant rates of character change. This implies that basal tetrapods had high, but declining, rates of homoplasy early in their evolutionary history, although the origin of Lissamphibia is an exception to this trend. The time slice approach is a powerful method of phylogenetic analysis and a useful tool for assessing the impact of combining extinct and extant taxa in phylogenetic analyses of large and speciose clades.

Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum

A genome-wide detection of phylogenetic signals by next generation sequencing (NGS) has recently emerged as a promising genomic approach for phylogenetic analysis of non-model organisms. Here we explored the use of a multiplexed shotgun sequencing method to assess the phylogenetic relationships of 18 Linum samples representing 16 species within four botanical sections of the flax genus Linum. The whole genome DNAs of 18 Linum samples were fragmented, tagged, and sequenced using an Illumina MiSeq. Acquired sequencing reads per sample were further separated into chloroplast, mitochondrial and nuclear sequence reads. SNP calls upon genome-specific sequence data sets revealed 6143 chloroplast, 2673 mitochondrial, and 19,562 nuclear SNPs. Phylogenetic analyses based on three-genome SNP data sets with and without missing observations showed congruent three-genome phylogenetic signals for four botanical sections of the Linum genus. Specifically, two major lineages showing a separation of Linum-Dasylinum sections and Linastrum-Syllinum sections were confirmed. The Linum section displayed three major branches representing two major evolutionary stages leading to cultivated flax. Cultivated flax and its immediate progenitor were formed as its own branch, genetically more closely related to L. decumbens and L. grandiflorum with chromosome count of eight, and distantly apart from six other species with chromosome count of nine. Five species of the Linastrum and Syllinum sections were genetically more distant from cultivated flax, but they appeared to be more closely related to each other, even with variable chromosome counts. These findings not only provide the first evidence of congruent three-genome phylogenetic pathways within the Linum genus, but also demonstrate the utility of the multiplexed shotgun sequencing in acquisition of three-genome phylogenetic signals of non-model organisms.

Vertebral Development in Paleozoic and Mesozoic Tetrapods Revealed by Paleohistological Data

Basal tetrapods display a wide spectrum of vertebral centrum morphologies that can be used to distinguish different tetrapod groups. The vertebral types range from multipartite centra in stem-tetrapods, temnospondyls, and seymouriamorphs up to monospondylous centra in lepospondyls and have been drawn upon for reconstructing major evolutionary trends in tetrapods that are now considered textbook knowledge. Two modes of vertebral formation have been postulated: the multipartite vertebrae formed first as cartilaginous elements with subsequent ossification. The monospondylous centrum, in contrast, was formed by direct ossification without a cartilaginous precursor. This study describes centrum morphogenesis in basal tetrapods for the first time, based on bone histology. Our results show that the intercentra of the investigated stem-tetrapods consist of a small band of periosteal bone and a dense network of endochondral bone. In stereospondyl temnospondyls, high amounts of calcified cartilage are preserved in the endochondral trabeculae. Notably, the periosteal region is thickened and highly vascularized in the plagiosaurid stereospondyls. Among “microsaur” lepospondyls, the thickened periosteal region is composed of compact bone and the notochordal canal is surrounded by large cell lacunae. In nectridean lepospondyls, the periosteal region has a spongy structure with large intertrabecular spaces, whereas the endochondral region has a highly cancellous structure. Our observations indicate that regardless of whether multipartite or monospondylous, the centra of basal tetrapods display first endochondral and subsequently periosteal ossification. A high interspecific variability is observed in growth rate, organization, and initiation of periosteal ossification. Moreover, vertebral development and structure reflect different lifestyles. The bottom-dwelling Plagiosauridae increase their skeletal mass by hyperplasy of the periosteal region. In nectrideans, the skeletal mass decreases, as the microstructure is spongy and lightly built. Additionally, we observed that vertebral structure is influenced by miniaturization in some groups. The phylogenetic information that can be drawn from vertebral development, however, is limited.

Cranial Morphology of the Brachystelechid 'Microsaur' Quasicaecilia texana Carroll Provides New Insights into the Diversity and Evolution of Braincase Morphology in Recumbirostran 'Microsaurs'

Recumbirostran ‘microsaurs,’ a group of early tetrapods from the Late Carboniferous and Early Permian, are the earliest known example of adaptation to head-first burrowing in the tetrapod fossil record. However, understanding of the diversity of fossorial adaptation within the Recumbirostra has been hindered by poor anatomical knowledge of the more divergent forms within the group. Here we report the results of μCT study of Quasicaecilia texana, a poorly-known recumbirostran with a unique, broad, shovel-like snout. The organization of the skull roof and braincase of Quasicaecilia is found to be more in line with that of other recumbirostrans than previously described, despite differences in overall shape. The braincase is found to be broadly comparable to Carrolla craddocki, with a large presphenoid that encompasses much of the interorbital septum and the columella ethmoidalis, and a single compound ossification encompassing the sphenoid, otic, and occipital regions. The recumbirostran braincase conserves general structure and topology of braincase regions and cranial nerve foramina, but it is highly variable in the number of ossifications and their extent, likely associated with the reliance on braincase ossifications to resist compression during sediment compaction and mechanical manipulation by epaxial and hypaxial musculature. Expansion of the deep ventral neck musculature in Quasicaecilia, autapomorphic among recumbirostrans, may reflect unique biomechanical function, and underscores the importance of future attention to the role of the cervical musculature in contextualizing the origin and evolution of fossoriality in recumbirostrans.

Skeletal Morphogenesis of Microbrachis and Hyloplesion (Tetrapoda: Lepospondyli), and Implications for the Developmental Patterns of Extinct, Early Tetrapods

The ontogeny of extant amphibians often is used as a model for that of extinct early tetrapods, despite evidence for a spectrum of developmental modes in temnospondyls and a paucity of ontogenetic data for lepospondyls. I describe the skeletal morphogenesis of the extinct lepospondyls Microbrachis pelikani and Hyloplesion longicostatum using the largest samples examined for either taxon. Nearly all known specimens were re-examined, allowing for substantial anatomical revisions that affect the scoring of characters commonly used in phylogenetic analyses of early tetrapods. The palate of H. longicostatum is re-interpreted and suggested to be more similar to that of M. pelikani, especially in the nature of the contact between the pterygoids. Both taxa possess lateral lines, and M. pelikani additionally exhibits branchial plates. However, early and rapid ossification of the postcranial skeleton, including a well-developed pubis and ossified epipodials, suggests that neither taxon metamorphosed nor were they neotenic in the sense of branchiosaurids and salamanders. Morphogenetic patterns in the foot suggest that digit 5 was developmentally delayed and the final digit to ossify in M. pelikani and H. longicostatum. Overall patterns of postcranial ossification may indicate postaxial dominance in limb and digit formation, but also more developmental variation in early tetrapods than has been appreciated. The phylogenetic position and developmental patterns of M. pelikani and H. longicostatum are congruent with the hypothesis that early tetrapods lacked metamorphosis ancestrally and that stem-amniotes exhibited derived features of development, such as rapid and complete ossification of the skeleton, potentially prior to the evolution of the amniotic egg.

Micro-CT Study of Rhynchonkos stovalli (Lepospondyli, Recumbirostra), with Description of Two New Genera

The Early Permian recumbirostran lepospondyl Rhynchonkos stovalli has been identified as a possible close relative of caecilians due to general similarities in skull shape as well as similar robustness of the braincase, a hypothesis that implies the polyphyly of extant lissamphibians. In order to better assess this phylogenetic hypothesis, we studied the morphology of the holotype and three specimens previously attributed to R. stovalli. With the use of micro-computed x-ray tomography (μCT) we are able to completely describe the external and internal cranial morphology of these specimens, dramatically revising our knowledge of R. stovalli and recognizing two new taxa, Aletrimyti gaskillae gen et sp. n. and Dvellacanus carrolli gen et sp. n. The braincases of R. stovalli, A. gaskillae, and D. carrolli are described in detail, demonstrating detailed braincase morphology and new information on the recumbirostran supraoccipital bone. All three taxa show fossorial adaptations in the braincase, sutural articulations of skull roof bones, and in the lower jaw, but variation in cranial morphology between these three taxa may reflect different modes of head-first burrowing behaviors and capabilities. We revisit the homology of the supraoccipital, median anterior bone, and temporal bone of recumbirostrans, and discuss implications of alternate interpretations of the homology of these elements. Finally, we evaluate the characteristics previously used to unite Rhynchonkos stovalli with caecilians in light of these new data. These proposed similarities are more ambiguous than previous descriptions suggest, and result from the composite nature of previous descriptions, ambiguities in external morphology, and functional convergence between recumbirostrans and caecilians for head-first burrowing.

A new species of the basal "kangaroo" Balbaroo and a re-evaluation of stem macropodiform interrelationships

Exceptionally well-preserved skulls and postcranial elements of a new species of the plesiomorphic stem macropodiform Balbaroo have been recovered from middle Miocene freshwater limestone deposits in the Riversleigh World Heritage Area of northwestern Queensland, Australia. This constitutes the richest intraspecific sample for any currently known basal "kangaroo", and, along with additional material referred to Balbaroo fangaroo, provides new insights into structural variability within the most prolific archaic macropodiform clade--Balbaridae. Qualitative and metric evaluations of taxonomic boundaries demonstrate that the previously distinct species Nambaroo bullockensis is a junior synonym of B. camfieldensis. Furthermore, coupled Maximum Parsimony and Bayesian phylogenetic analyses reveal that our new Balbaroo remains represent the most derived member of the Balbaroo lineage, and are closely related to the middle Miocene B. camfieldensis, which like most named balbarid species is identifiable only from isolated jaws. The postcranial elements of Balbaroo concur with earlier finds of the stratigraphically oldest balbarid skeleton, Nambaroo gillespieae, and suggest that quadrupedal progression was a primary gait mode as opposed to bipedal saltation. All Balbaroo spp. have low-crowned bilophodont molars, which are typical for browsing herbivores inhabiting the densely forested environments envisaged for middle Miocene northeastern Australia.

New phylogenomic and comparative analyses provide corroborating evidence that Myxozoa is Cnidaria

Myxozoa, a diverse group of morphologically simplified endoparasites, are well known fish parasites causing substantial economic losses in aquaculture. Despite active research, the phylogenetic position of Myxozoa remains ambiguous. After obtaining the genome and transcriptome data of the myxozoan Thelohanellus kitauei, we examined the phylogenetic position of Myxozoa from three different perspectives. First, phylogenomic analyses with the newly sequenced genomic data strongly supported the monophyly of Myxozoa and that Myxozoa is sister to Medusozoa within Cnidaria. Second, we detected two homologs to cnidarian-specific minicollagens in the T. kitauei genome with molecular characteristics similar to cnidarian-specific minicollagens, suggesting that the minicollagen homologs in T. kitauei may have functions similar to those in Cnidaria and that Myxozoa is Cnidaria. Additionally, phylogenetic analyses revealed that the minicollagens in myxozoans and medusozoans have a common ancestor. Third, we detected 11 of the 19 proto-mesodermalgenes in the T. kitauei genome, which were also present in the cnidarian Hydra magnipapillata, indicating Myxozoa is within Cnidaria. Thus, our results robustly support Myxozoa as a derived cnidarian taxon with an affinity to Medusozoa, helping to understand the diversity of the morphology, development and life cycle of Cnidaria and its evolution.

Recovering missing data: estimating position and size of caudal vertebrae in Staurikosaurus pricei Colbert, 1970

Missing data is a common problem in paleontology. It makes it difficult to reconstruct extinct taxa accurately and restrains the inclusion of some taxa on comparative and biomechanical studies. Particularly, estimating the position of vertebrae on incomplete series is often non-empirical and does not allow precise estimation of missing parts. In this work we present a method for calculating the position of preserved middle sequences of caudal vertebrae in the saurischian dinosaur Staurikosaurus pricei, based on the length and height of preserved anterior and posterior caudal vertebral centra. Regression equations were used to estimate these dimensions for middle vertebrae and, consequently, to assess the position of the preserved middle sequences. It also allowed estimating these dimensions for non-preserved vertebrae. Results indicate that the preserved caudal vertebrae of Staurikosaurus may correspond to positions 1-3, 5, 7, 14-19/15-20, 24-25/25-26, and 29-47, and that at least 25 vertebrae had transverse processes. Total length of the tail was estimated in 134 cm and total body length was 220-225 cm.

Osmotic tolerance and habitat of early stegocephalians: indirect evidence from parsimony, taphonomy, palaeobiogeography, physiology and morphology

There are probably many reasons for the widespread belief that temnospondyls and other early stegocephalians were largely restricted to freshwater, but three of the contributing factors will be discussed below. First, temnospondyls have been called amphibians (and thought to be more closely related to extant amphibians than to amniotes). Some authors may have simply concluded that, like extant amphibians, temnospondyls could not live in oceans and seas. Second, under some phylogenies, temnospondyls are more closely related to anurans (and possibly urodeles) than to gymnophionans and could be expected, for parsimony reasons, to share the intolerance of all extant amphibians to saltwater. Similarly, ‘lepospondyls’ are often thought to be more closely related to gymnophionans than to anurans, and could also be expected to lack saltwater tolerance. Third, extant lungfishes live exclusively in freshwater, and early sarcopterygians have long been thought to share this habitat. These interpretations probably explain the widespread belief that early amphibians and early stem-tetrapods were largely restricted to freshwater. However, these three interpretations have been refuted or questioned by recent investigations. A review of the evidence suggests that several (perhaps most) early stegocephalians tolerated saltwater, even although they also lived in freshwater.

Evolution of ossification sequences in salamanders and urodele origins assessed through event-pairing and new methods

Ossification sequences of the skull in extant Urodela and in Permo-Carboniferous Branchiosauridae have already been used to study the origin of lissamphibians. But most of these studies did not consider some recent methods developed to analyze the developmental sequences within a phylogenetic framework. Here, we analyze the ossification sequences of 24 cranial bones of 23 extant species of salamanders using the event-pairing method. This reveals new developmental synapomorphies for several extant salamander taxa and ancestral sequences for Urodela under four alternative reference phylogenies. An analysis with the 12 bones for which ossification sequence data are available in urodeles and in the branchiosaurid Apateon is also performed in order to compare the ancestral condition of the crown-group of Urodela to the sequence of Apateon. This reveals far more incompatibilities than previously suggested. The similarities observed between some extant salamanders and branchiosaurids may result from extensive homoplasy, as the extreme variation observed in extant Urodela suggests, or be plesiomorphic, as the conservation of some ossification patterns observed in other remotely related vertebrates like actinopterygians suggests. We propose a new, simpler method based on squared-change optimization to estimate the relative timing of ossification of various bones of hypothetical ancestors, and use independent-contrasts analysis to estimate the confidence intervals around these times. Our results show that the uncertainty of the ancestral ossification sequence of Urodela is much greater than event-pairing suggests. The developmental data do not allow to conclude that branchiosaurids are closely related to salamanders and their limited taxonomic distribution in Paleozoic taxa precludes testing hypotheses about lissamphibian origins. This is true regardless of the analytical method used (event-pairing or our new method based on squared-change parsimony). Simulations show that the new analytical method is generally more powerful to detect evolutionary shifts in developmental timing, and has lower Type I error rate than event-pairing. It also makes fewer errors in ancestral character value or state assignment than event-pairing.

Paleontology, genomics, and combined-data phylogenetics: can molecular data improve phylogeny estimation for fossil taxa?

The genomics revolution offers great promise for resolving the phylogeny of living taxa, but does it offer any benefits for reconstructing relationships among extinct (fossil) taxa? Superficially, the answer would seem to be "no," given that molecular data cannot be obtained for most fossil taxa. However, because fossil taxa often interdigitate among living taxa on the Tree of Life, molecular data may indirectly enhance phylogenetic accuracy for fossil taxa in the context of a combined analysis of morphological and molecular data for living and fossil taxa. Here, I use simulations to assess accuracy for fossil taxa in a mixed analysis of living and fossil taxa, before and after addition of molecular data to the living taxa. The results show conditions where the accuracy for fossil taxa is greatly increased by adding molecular data, sometimes by as much as 100%. In other cases, the increase is negligible, such as when fossil taxa greatly outnumber living taxa in the analysis. However, there were few cases where accuracy was significantly decreased by the addition of the molecular data, suggesting that this practice may range from highly beneficial to mostly harmless. Overall, the results suggest that improvements in molecular phylogenetics can potentially benefit phylogeny reconstruction for fossil taxa.

The armoured dissorophid Cacops from the Early Permian of Oklahoma and the exploitation of the terrestrial realm by amphibians

Cacops, one of the most distinctive Paleozoic amphibians, is part of a clade of dissorophoid temnospondyls that diversified in the equatorial region of Pangea during the Late Carboniferous and Early Permian, persisting into the Late Permian in Central Russia and China. Dissorophids were a successful group of fully terrestrial, often spectacularly armoured predators, the only amphibians apparently able to coexist with amniotes when the latter started to dominate terrestrial ecosystems. In this paper, we describe excellent new skulls from the Early Permian of Oklahoma attributed to Cacops, Cacops morrisi sp. nov. and provide for the first time detailed information about this iconic dissorophid. These specimens show anatomical and ontogenetic features that will impact on future studies on the evolution of terrestriality in tetrapods. For example, the large, posteriorly closed tympanic embayment has fine striations on an otherwise smooth surface, documenting the oldest known clear evidence for the presence of a tympanic membrane in the fossil record, a structure that is used for hearing airborne sound in extant tetrapods. The skull of C. morrisi also has several features associated with predatory behaviour, indicating that this dissorophid may have been one of the top terrestrial predators of its time.

A hierarchical model for incomplete alignments in phylogenetic inference

MOTIVATION

Full-length DNA and protein sequences that span the entire length of a gene are ideally used for multiple sequence alignments (MSAs) and the subsequent inference of their relationships. Frequently, however, MSAs contain a substantial amount of missing data. For example, expressed sequence tags (ESTs), which are partial sequences of expressed genes, are the predominant source of sequence data for many organisms. The patterns of missing data typical for EST-derived alignments greatly compromise the accuracy of estimated phylogenies.

RESULTS

We present a statistical method for inferring phylogenetic trees from EST-based incomplete MSA data. We propose a class of hierarchical models for modeling pairwise distances between the sequences, and develop a fully Bayesian approach for estimation of the model parameters. Once the distance matrix is estimated, the phylogenetic tree may be constructed by applying neighbor-joining (or any other algorithm of choice). We also show that maximizing the marginal likelihood from the Bayesian approach yields similar results to a profile likelihood estimation. The proposed methods are illustrated using simulated protein families, for which the true phylogeny is known, and one real protein family.

AVAILABILITY

R code for fitting these models are available from: http://people.bu.edu/gupta/software.htm.

Ever Since Owen: Changing Perspectives on the Early Evolution of Tetrapods

The traditional notion of a gap between fishes and amphibians has been closed by a wealth of fish-like fossil tetrapods, many discovered since the mid 1980s. This review summarizes these discoveries and explores their significance relative to changing ideas about early tetrapod phylogeny, biogeography, and ecology. Research emphasis can now shift to broader-based questions, including the whole of the early tetrapod radiation, from the divergence from other lobed-finned fishes to the origins of modern amphibians and amniotes. The fish-to-tetrapod morphological transition occurred within the Upper Devonian; the divergence of modern tetrapod groups is an Early Carboniferous event. Modern tetrapods emerged in the aftermath of one of the five major extinction episodes in the fossil record, but the earlier Devonian tetrapod radiation is not well understood. Tetrapod limbs, paired fins, and comparative developmental data are reviewed; again, research emphasis needs to change to explore the origins of tetrapod diversity.

Phylogenomic analyses support the monophyly of Taphrinomycotina, including Schizosaccharomyces fission yeasts

Several morphologically dissimilar ascomycete fungi including Schizosaccharomyces, Taphrina, Saitoella, Pneumocystis, and Neolecta have been grouped into the taxon Taphrinomycotina (Archiascomycota or Archiascomycotina), originally based on rRNA phylogeny. These analyses lack statistically significant support for the monophyly of this grouping, and although confirmed by more recent multigene analyses, this topology is contradicted by mitochondrial phylogenies. To resolve this inconsistency, we have assembled phylogenomic mitochondrial and nuclear data sets from four distantly related taphrinomycotina taxa: Schizosaccharomyces pombe, Pneumocystis carinii, Saitoella complicata, and Taphrina deformans. Our phylogenomic analyses based on nuclear data (113 proteins) conclusively support the monophyly of Taphrinomycotina, diverging as a sister group to Saccharomycotina + Pezizomycotina. However, despite the improved taxon sampling, Taphrinomycotina continue to be paraphyletic with the mitochondrial data set (13 proteins): Schizosaccharomyces species associate with budding yeasts (Saccharomycotina) and the other Taphrinomycotina group as a sister group to Saccharomycotina + Pezizomycotina. Yet, as Schizosaccharomyces and Saccharomycotina species are fast evolving, the mitochondrial phylogeny may be influenced by a long-branch attraction (LBA) artifact. After removal of fast-evolving sequence positions from the mitochondrial data set, we recover the monophyly of Taphrinomycotina. Our combined results suggest that Taphrinomycotina is a legitimate taxon, that this group of species diverges as a sister group to Saccharomycotina + Pezizomycotina, and that phylogenetic positioning of yeasts and fission yeasts with mitochondrial data is plagued by a strong LBA artifact.

Using ESTs for phylogenomics: can one accurately infer a phylogenetic tree from a gappy alignment?

BACKGROUND

While full genome sequences are still only available for a handful of taxa, large collections of partial gene sequences are available for many more. The alignment of partial gene sequences results in a multiple sequence alignment containing large gaps that are arranged in a staggered pattern. The consequences of this pattern of missing data on the accuracy of phylogenetic analysis are not well understood. We conducted a simulation study to determine the accuracy of phylogenetic trees obtained from gappy alignments using three commonly used phylogenetic reconstruction methods (Neighbor Joining, Maximum Parsimony, and Maximum Likelihood) and studied ways to improve the accuracy of trees obtained from such datasets.

RESULTS

We found that the pattern of gappiness in multiple sequence alignments derived from partial gene sequences substantially compromised phylogenetic accuracy even in the absence of alignment error. The decline in accuracy was beyond what would be expected based on the amount of missing data. The decline was particularly dramatic for Neighbor Joining and Maximum Parsimony, where the majority of gappy alignments contained 25% to 40% incorrect quartets. To improve the accuracy of the trees obtained from a gappy multiple sequence alignment, we examined two approaches. In the first approach, alignment masking, potentially problematic columns and input sequences are excluded from from the dataset. Even in the absence of alignment error, masking improved phylogenetic accuracy up to 100-fold. However, masking retained, on average, only 83% of the input sequences. In the second approach, alignment subdivision, the missing data is statistically modelled in order to retain as many sequences as possible in the phylogenetic analysis. Subdivision resulted in more modest improvements to alignment accuracy, but succeeded in including almost all of the input sequences.

CONCLUSION

These results demonstrate that partial gene sequences and gappy multiple sequence alignments can pose a major problem for phylogenetic analysis. The concern will be greatest for high-throughput phylogenomic analyses, in which Neighbor Joining is often the preferred method due to its computational efficiency. Both approaches can be used to increase the accuracy of phylogenetic inference from a gappy alignment. The choice between the two approaches will depend upon how robust the application is to the loss of sequences from the input set, with alignment masking generally giving a much greater improvement in accuracy but at the cost of discarding a larger number of the input sequences.

Fossils Impact as Hard as Living Taxa in Parsimony Analyses of Morphology

Systematists disagree whether data from fossils should be included in parsimony analyses. In a handful of well-documented cases, the addition of fossil data radically overturns a hypothesis of relationships based on extant taxa alone. Fossils can break up long branches and preserve character combinations closer in time to deep splitting events. However, fossils usually require more interpretation than extant taxa, introducing greater potential for spurious codings. Moreover, because fossils often have more "missing" codings, they are frequently accused of increasing numbers of MPTs, frustrating resolution and reducing support. Despite the controversy, remarkably little is known about the effects of fossils more generally. Here we provide the first systematic study, investigating empirically the behavior of fossil and extant taxa in 45 published morphological data sets. First-order jackknifing is used to determine the effects that each terminal has on inferred relationships, on the number of MPTs, and on CI' and RI as measures of homoplasy. Bootstrap leaf stabilities provide a proxy for the contribution of individual taxa to the branch support in the rest of the tree. There is no significant difference in the impact of fossil versus extant taxa on relationships, numbers of MPTs, and CI' or RI. However, adding individual fossil taxa is more likely to reduce the total branch support of the tree than adding extant taxa. This must be weighed against the superior taxon sampling afforded by including judiciously coded fossils, providing data from otherwise unsampled regions of the tree. We therefore recommend that investigators should include fossils, in the absence of compelling and case specific reasons for their exclusion.

Fossils, Molecules, Divergence Times, and the Origin of Lissamphibians

A review of the paleontological literature shows that the early dates of appearance of Lissamphibia recently inferred from molecular data do not favor an origin of extant amphibians from temnospondyls, contrary to recent claims. A supertree is assembled using new Mesquite modules that allow extinct taxa to be incorporated into a time-calibrated phylogeny with a user-defined geological time scale. The supertree incorporates 223 extinct species of lissamphibians and has a highly significant stratigraphic fit. Some divergences can even be dated with sufficient precision to serve as calibration points in molecular divergence date analyses. Fourteen combinations of minimal branch length settings and 10 random resolutions for each polytomy give much more recent minimal origination times of lissamphibian taxa than recent studies based on a phylogenetic analyses of molecular sequences. Attempts to replicate recent molecular date estimates show that these estimates depend strongly on the choice of calibration points, on the dating method, and on the chosen model of evolution; for instance, the estimate for the date of the origin of Lissamphibia can lie between 351 and 266 Mya. This range of values is generally compatible with our time-calibrated supertree and indicates that there is no unbridgeable gap between dates obtained using the fossil record and those using molecular evidence, contrary to previous suggestions.

Phylogeny of Extant and Fossil Juglandaceae Inferred from the Integration of Molecular and Morphological Data Sets

It is widely acknowledged that integrating fossils into data sets of extant taxa is imperative for proper placement of fossils, resolution of relationships, and a better understanding of character evolution. The importance of this process has been further magnified because of the crucial role of fossils in dating divergence times. Outstanding issues remain, including appropriate methods to place fossils in phylogenetic trees, the importance of molecules versus morphology in these analyses, as well as the impact of potentially large amounts of missing data for fossil taxa. In this study we used the angiosperm clade Juglandaceae as a model for investigating methods of integrating fossils into a phylogenetic framework of extant taxa. The clade has a rich fossil record relative to low extant diversity, as well as a robust molecular phylogeny and morphological database for extant taxa. After combining fossil organ genera into composite and terminal taxa, our objectives were to (1) compare multiple methods for the integration of the fossils and extant taxa (including total evidence, molecular scaffolds, and molecular matrix representation with parsimony [MRP]); (2) explore the impact of missing data (incomplete taxa and characters) and the evidence for placing fossils on the topology; (3) simulate the phylogenetic effect of missing data by creating "artificial fossils"; and (4) place fossils and compare the impact of single and multiple fossil constraints in estimating the age of clades. Despite large and variable amounts of missing data, each of the methods provided reasonable placement of both fossils and simulated "artificial fossils" in the phylogeny previously inferred only from extant taxa. Our results clearly show that the amount of missing data in any given taxon is not by itself an operational guideline for excluding fossils from analysis. Three fossil taxa (Cruciptera simsonii, Paleoplatycarya wingii, and Platycarya americana) were placed within crown clades containing living taxa for which relationships previously had been suggested based on morphology, whereas Polyptera manningii, a mosaic taxon with equivocal affinities, was placed firmly as sister to two modern crown clades. The position of Paleooreomunnea stoneana was ambiguous with total evidence but conclusive with DNA scaffolds and MRP. There was less disturbance of relationships among extant taxa using a total evidence approach, and the DNA scaffold approach did not provide improved resolution or internal support for clades compared to total evidence, whereas weighted MRP retained comparable levels of support but lost crown clade resolution. Multiple internal minimum age constraints generally provided reasonable age estimates, but the use of single constraints provided by extinct genera tended to underestimate clade ages.

SCaFoS: a tool for selection, concatenation and fusion of sequences for phylogenomics

BACKGROUND

Phylogenetic analyses based on datasets rich in both genes and species (phylogenomics) are becoming a standard approach to resolve evolutionary questions. However, several difficulties are associated with the assembly of large datasets, such as multiple copies of a gene per species (paralogous or xenologous genes), lack of some genes for a given species, or partial sequences. The use of undetected paralogous or xenologous genes in phylogenetic inference can lead to inaccurate results, and the use of partial sequences to a lack of resolution. A tool that selects sequences, species, and genes, while dealing with these issues, is needed in a phylogenomics context.

RESULTS

Here, we present SCaFoS, a tool that quickly assembles phylogenomic datasets containing maximal phylogenetic information while adjusting the amount of missing data in the selection of species, sequences and genes. Starting from individual sequence alignments, and using monophyletic groups defined by the user, SCaFoS creates chimeras with partial sequences, or selects, among multiple sequences, the orthologous and/or slowest evolving sequences. Once sequences representing each predefined monophyletic group have been selected, SCaFos retains genes according to the user's allowed level of missing data and generates files for super-matrix and super-tree analyses in several formats compatible with standard phylogenetic inference software. Because no clear-cut criteria exist for the sequence selection, a semi-automatic mode is available to accommodate user's expertise.

CONCLUSION

SCaFos is able to deal with datasets of hundreds of species and genes, both at the amino acid or nucleotide level. It has a graphical interface and can be integrated in an automatic workflow. Moreover, SCaFoS is the first tool that integrates user's knowledge to select orthologous sequences, creates chimerical sequences to reduce missing data and selects genes according to their level of missing data. Finally, applying SCaFoS to different datasets, we show that the judicious selection of genes, species and sequences reduces tree reconstruction artefacts, especially if the dataset includes fast evolving species.

Challenges of phylogenetic analyses of aDNA sequences

One of the crucial steps of authentication of aDNA sequences is phylogenetic consistency. Amplified sequences should fit into the phylogenetic framework of their supposed origin. An inherent property of aDNA sequences however, is their short sequence length. Additionally, genes for aDNA studies are often chosen by their preservation potential rather than by phylogenetically informative content. This poses potential challenges regarding their analyses, and might result in an inaccurate reflection of the supposed phylogenetic history of the sequence or organism under study. In this paper some fundamental problems of phylogenetic analysis and interpretation of aDNA datasets are discussed. Suggestions for character sampling and treatment of missing data are made. The publication is the result of a talk from the 1st PAMINSA Meeting in Rio de Janeiro, July 2005.

Can incomplete taxa rescue phylogenetic analyses from long-branch attraction?

Taxon sampling may be critically important for phylogenetic accuracy because adding taxa can help to subdivide misleading long branches. Although the idea that added taxa can break up long branches was exemplified by a study of "incomplete" fossil taxa, the issue of taxon completeness (i.e., proportion of missing data) has been largely ignored in most subsequent discussions of taxon sampling and long-branch attraction. In this article, I use simulations to test the ability of incomplete taxa to subdivide long branches and improve phylogenetic accuracy in situations of potential long-branch attraction. The results show that for most methods and conditions examined, adding taxa that are only 50% complete may provide similar benefits to adding the same number of complete taxa (suggesting that the advantages of increased taxon sampling may be obtained with less data than previously considered). For parsimony, taxa that are less complete (5% to 25% complete) may often have limited ability to rescue analyses from long-branch attraction. In contrast, highly incomplete taxa can be surprisingly beneficial when using model-based methods. The results also suggest the importance of model-based methods in phylogenetic analyses that combine molecular and fossil data.

Revision of Nicrophorus in part: new species and inferred phylogeny of the nepalensis-group based on evidence from morphology and mitochondrial DNA (Coleoptera:Silphidae:Nicrophorinae)

Carrion beetles of the genus Nicrophorus Fabricius, 1775 (Silphidae) are well known for their biparental brood care and monopolisation of small vertebrate carcasses in subterranean crypts. Although the taxonomy of New World species has received modern attention, the fauna of Asia, primarily of the nepalensis-group of species, has not. Herein we revise this species-group and describe as new the following seven species: Nicrophorus charon Sikes & Madge (Sulawesi), Nicrophorus herscheli Sikes & Madge (Sumatra), Nicrophorus insignis Sikes & Madge (Flores Island), Nicrophorus melissae Sikes & Madge (Nepal, Bhutan), Nicrophorus reticulatus Sikes & Madge (Guadalcanal), Nicrophorus schawalleri Sikes & Madge (Gansu, Shaanxi, Sichuan Province) and Nicrophorus trumboi Sikes & Madge (Nepal, Bhutan). We obtained a preliminary phylogeny using morphology and mtDNA (COII). This was inferred using maximum likelihood and Bayesian methods with the Mkv and GTR+I+G models (parsimony was rejected by the Akaike information criterion for being excessively parameter-rich). The phylogenetic signal in the morphological dataset was not strong and results were confounded by a ‘long-branch’ species, N. reticulatus. The signal was stronger in the combined dataset and the COII-only dataset. The molecular phylogeny supported the new status of species N. trumboi and N. melissae. Support was found for a mainland origin of the group with subsequent radiations into the Malay Archipelago.

The evolution of body size, Cope's rule and the origin of amniotes

The evolution of body size in tetrapods is assessed using a database that includes 107 early stegocephalian species ranging in time from the Frasnian (Upper Devonian) to the Tatarian (Upper Permian). All analyses use methods that incorporate phylogenetic information (topology and branch lengths). In all tests, the impact of alternative topologies and branch lengths are assessed. Previous reports that raised doubts about the accuracy of squared-change parsimony assessment of ancestral character value appear to have used datasets in which there was no phylogenetic signal. Hence, squared-change parsimony may be more reliable than suggested in recent studies, at least when a phylogenetic signal is present in the datasets of interest. Analysis using random taxon reshuffling on three reference phylogenies shows that cranial and presacral length include a strong phylogenetic signal. Character optimization of body size in stegocephalians using squared-change parsimony on a time-calibrated phylogeny incorporating branch length information is used to test a previously published scenario on the origin of amniotes and of the amniotic egg that implies that the ancestors of amniotes were small (no more than 10 cm in snout-vent length), and that their size increased subsequent to the appearance of the amniotic egg. The optimization suggests that first amniotes were somewhat larger than previously hypothesized; the estimated snout-vent length is about 24 cm, and the lower end of the 95% confidence interval of the phylogeny that yields the smallest inferred size suggests that no ancestor of amniotes measured less than 12 cm in snout-vent length. Character optimization, permutational multiple linear regressions, and independent contrast analyses show that Cope's rule of phyletic size increase applies to early reptiliomorphs but that it does not apply to early stegocephalians globally.

Phylogenomics of Eukaryotes: Impact of Missing Data on Large Alignments

Resolving the relationships between Metazoa and other eukaryotic groups as well as between metazoan phyla is central to the understanding of the origin and evolution of animals. The current view is based on limited data sets, either a single gene with many species (e.g., ribosomal RNA) or many genes but with only a few species. Because a reliable phylogenetic inference simultaneously requires numerous genes and numerous species, we assembled a very large data set containing 129 orthologous proteins ( approximately 30,000 aligned amino acid positions) for 36 eukaryotic species. Included in the alignments are data from the choanoflagellate Monosiga ovata, obtained through the sequencing of about 1,000 cDNAs. We provide conclusive support for choanoflagellates as the closest relative of animals and for fungi as the second closest. The monophyly of Plantae and chromalveolates was recovered but without strong statistical support. Within animals, in contrast to the monophyly of Coelomata observed in several recent large-scale analyses, we recovered a paraphyletic Coelamata, with nematodes and platyhelminths nested within. To include a diverse sample of organisms, data from EST projects were used for several species, resulting in a large amount of missing data in our alignment (about 25%). By using different approaches, we verify that the inferred phylogeny is not sensitive to these missing data. Therefore, this large data set provides a reliable phylogenetic framework for studying eukaryotic and animal evolution and will be easily extendable when large amounts of sequence information become available from a broader taxonomic range.