Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia

Were molecular data available for extinct taxa, questions regarding the origins of many groups could be settled in short order. As this is not the case, various strategies have been proposed to combine paleontological and neontological data sets. The use of fossil dates as node age calibrations for divergence time estimation from molecular phylogenies is commonplace. In addition, simulations suggest that the addition of morphological data from extinct taxa may improve phylogenetic estimation when combined with molecular data for extant species, and some studies have merged morphological and molecular data to estimate combined evidence phylogenies containing both extinct and extant taxa. However, few, if any, studies have attempted to estimate divergence times using phylogenies containing both fossil and living taxa sampled for both molecular and morphological data. Here, I infer both the phylogeny and the time of origin for Lissamphibia and a number of stem tetrapods using Bayesian methods based on a data set containing morphological data for extinct taxa, molecular data for extant taxa, and molecular and morphological data for a subset of extant taxa. The results suggest that Lissamphibia is monophyletic, nested within Lepospondyli, and originated in the late Carboniferous at the earliest. This research illustrates potential pitfalls for the use of fossils as post hoc age constraints on internal nodes and highlights the importance of explicit phylogenetic analysis of extinct taxa. These results suggest that the application of fossils as minima or maxima on molecular phylogenies should be supplemented or supplanted by combined evidence analyses whenever possible.     

Phylogeny of the owlet-nightjars (Aves: Aegothelidae) based on mitochondrial DNA sequence

The avian family Aegothelidae (Owlet-nightjars) comprises nine extant species and one extinct species, all of which are currently classified in a single genus, Aegotheles. Owlet-nightjars are secretive nocturnal birds of the South Pacific. They are relatively poorly studied and some species are known from only a few specimens. Furthermore, their confusing morphological variation has made it difficult to cluster existing specimens unambiguously into hierarchical taxonomic units. Here we sample all extant owlet-nightjar species and all but three currently recognized subspecies. We use DNA extracted primarily from museum specimens to obtain mitochondrial gene sequences and construct a molecular phylogeny. Our phylogeny suggests that most species are reciprocally monophyletic, however A. albertisi appears paraphyletic. Our data also suggest splitting A. bennettii into two species and splitting A. insignis and A. tatei as suggested in another recent paper.     

Congruent Avian Phylogenies Inferred from Mitochondrial and Nuclear DNA Sequences

Recent molecular studies addressing the phylogenetic relationships of avian orders have had conflicting results. While studies using nuclear DNA sequences tend to support traditional taxonomic views, also supported by morphological data [(paleognaths (galloanseres (all other birds)))], with songbirds forming a clade within Neoaves (all other birds), analyses with complete mtDNA genomes have resulted in topologies that place songbirds as one of the earliest-diverging avian lineages. Considering that over half of the extant bird species are songbirds, these different results have very different implications for our understanding of avian evolution. We analyzed data sets comprising nearly 4 kb of mitochondrial DNA (mtDNA) (complete 12S, ND1, ND2, and cytochrome b) plus 600 bp of the nuclear gene c-mos for 15 birds that were chosen to represent all major avian clades and to minimize potential long-branch attraction problems; we used a partition-specific maximum likelihood approach. Our results show congruence with respect to the ingroup among phylogenies obtained with mtDNA and the nuclear gene c-mos, separately or combined. The data sets support a traditional avian taxonomy, with paleognaths (ratites and tinamous) occupying a basal position and with songbirds more derived and forming a monophyletic group. We also show that, for mtDNA studies, turtles may be a better outgroup for birds than crocodilians because of their slower rate of sequence evolution.     

A phylogenetic supertree of oscine passerine birds (Aves: Passeri)

Oscine passerine birds make up almost half of all avian diversity. Relationships within the group, and its classification, have long been controversial. Over the last 10 years numerous molecular phylogenies have been published. We compiled source phylogenies from 99 published studies to construct an oscine supertree. We aimed to illustrate weak and strong parts of the phylogeny and set targets for future phylogenetic work and therefore preferred a heuristic approach where we judged the adequacy of taxon sampling and molecular method of each source tree instead of using matrices and automated tree-building programs. We present an estimate of the phylogenetic relationships of 1723 extant and one extinct species of oscine passerine birds (Aves: Passeri) — more than 37% of the total. We included 34/35 (97%) families, 38/39 (97%) subfamilies and 40/43 (93%) tribes. Overall resolution is 83% of a fully bifurcating tree. The basal lineages are all distributed in the Australo-Papuan region, but several more distal lineages dispersed out of this region and radiated in other parts of the world. However, taxa of the Australian region suffer from larger evolutionary gaps and the deep branches of the Sylvioidea and nine South American primaried oscines are still poorly resolved.     

When rediscovery is not enough: taxonomic uncertainty hinders conservation of a critically endangered bird

In 2003, birds similar to the extinct New Zealand storm-petrel Oceanites maorianus were observed in Hauraki Gulf NZ, raising the possibility of rediscovery after 150 years. O. maorianus has and continues to be surrounded by taxonomic uncertainty, being variously described as a distinct genus, a distinct species, or merely a plumage variant. This uncertainty has hindered conservation planning and funding for the species. Here we examine the taxonomic identity of the rediscovered birds and museum specimens using phylogenetic analysis of mitochondrial (1143 bp cytochrome b) and nuclear (890 bp β-fibrinogen) genes. Using cytochrome b sequence amplified from the 150+ year old specimens, we found that the extant and museum O. maorianus were the same taxon (0.01 genetic distance), with both differing from all other storm-petrel taxa. Using both genes, we examined the phylogenetic affinities of O. maorianus to the Oceanitinae and Hydrobatinae storm-petrels finding that O. maorianus was more closely aligned to Fregetta (0.08-0.09) than Oceanites (0.11-0.12), thereby confirming its status as a distinct taxon, not a plumage variant of O. oceanicus. Our analysis verifies that the previously presumed extinct New Zealand storm-petrel has been rediscovered and can now be assigned a conservation priority commensurate with its critically endangered status.     

Artificial neural networks can learn to estimate extinction rates from molecular phylogenies

Molecular phylogenies typically consist of only extant species, yet they allow inference of past rates of extinction, because recently originated species are less likely to be extinct than ancient species. Despite the simple structure of the assumed underlying speciation-extinction process, parametric functions to estimate extinction rates from phylogenies turned out to be complex and often difficult to derive. Moreover, these parametric functions are specific to a particular process (e.g. complete species level phylogeny with constant birth and death rates) and a particular type of data (e.g. times between bifurcations). Here, it is shown that artificial neural networks can substitute for parametric estimation functions once they have been sufficiently trained on simulated data. This technique can in principle be used for different processes and data types, and because it circumvents the time-consuming and difficult task of deriving parametric estimation functions, it may greatly extend the possibilities to make macro-evolutionary inferences from molecular phylogenies. This novel approach is explained, applied to estimate speciation and extinction rates from a molecular phylogeny of the reef fish genus Naso (Acanturidae), and its performance is compared to that of maximum likelihood estimation.     

No substitute for real data: A cautionary note on the use of phylogenies from birth–death polytomy resolvers for downstream comparative analyses

The statistical estimation of phylogenies is always associated with uncertainty, and accommodating this uncertainty is an important component of modern phylogenetic comparative analysis. The birth–death polytomy resolver is a method of accounting for phylogenetic uncertainty that places missing (unsampled) taxa onto phylogenetic trees, using taxonomic information alone. Recent studies of birds and mammals have used this approach to generate pseudoposterior distributions of phylogenetic trees that are complete at the species level, even in the absence of genetic data for many species. Many researchers have used these distributions of phylogenies for downstream evolutionary analyses that involve inferences on phenotypic evolution, geography, and community assembly. I demonstrate that the use of phylogenies constructed in this fashion is inappropriate for many questions involving traits. Because species are placed on trees at random with respect to trait values, the birth–death polytomy resolver breaks down natural patterns of trait phylogenetic structure. Inferences based on these trees are predictably and often drastically biased in a direction that depends on the underlying (true) pattern of phylogenetic structure in traits. I illustrate the severity of the phenomenon for both continuous and discrete traits using examples from a global bird phylogeny.     

V.PhyloMaker: an R package that can generate very large phylogenies for vascular plants

We present V.PhyloMaker, a freely available package for R designed to generate phylogenies for vascular plants. The mega‐tree implemented in V.PhyloMaker (i.e. GBOTB.extended.tre), which was derived from two recently published mega‐trees and includes 74 533 species and all families of extant vascular plants, is the largest dated phylogeny for vascular plants. V.PhyloMaker can generate phylogenies for very large species lists (the largest species list that we tested included 314 686 species). V.PhyloMaker generates phylogenies at a fast speed, much faster than other phylogeny‐generating packages. Our tests of V.PhyloMaker show that generating a phylogeny for 60 000 species requires less than six hours. V.PhyloMaker includes an approach to attach genera or species to their close relatives in a phylogeny. We provide a simple example in this paper to show how to use V.PhyloMaker to generate phylogenies.     

The demography of extinction in eastern North American birds

Species are being lost at an unprecedented rate during the Anthropocene. Progress has been made in clarifying how species traits influence their propensity to go extinct, but the role historical demography plays in species loss or persistence is unclear. In eastern North America, five charismatic landbirds went extinct last century, and the causes of their extinctions have been heavily debated. Although these extinctions are most often attributed to post-colonial human activity, other factors such as declining ancestral populations prior to European colonization could have made these species particularly susceptible. We used population genomic data from these extinct birds and compared them with those from four codistributed extant species. We found extinct species harboured lower genetic diversity and effective population sizes than extant species, but both extinct and non-extinct birds had similar demographic histories of population expansion. These demographic patterns are consistent with population size changes associated with glacial–interglacial cycles. The lack of support for overall population declines during the Pleistocene corroborates the view that, although species that went extinct may have been vulnerable due to low diversity or small population size, their disappearance was driven by human activities in the Anthropocene.     

Sampling-through-time in birth-death trees

I consider the constant rate birth-death process with incomplete sampling. I calculate the density of a given tree with sampled extant and extinct individuals.This density is essential for analyzing datasets which are sampled through time. Such datasets are common in virus epidemiology as viruses in infected individuals are sampled through time. Further, such datasets appear in phylogenetics when extant and extinct species data is available.I show how the derived tree density can be used (i) as a tree prior in a Bayesian method to reconstruct the evolutionary past of the sequence data on a calender-timescale, (ii) to infer the birth- and death-rates for a reconstructed evolutionary tree, and (iii) for simulating trees with a given number of sampled extant and extinct individuals which is essential for testing evolutionary hypotheses for the considered datasets.     

Evolution of fossil and living spider flies based on morphological and molecular data (Diptera,Acroceridae)

The phylogeny of spider flies is presented based on an analysis of DNA sequence data combined with morphological characters for both living and fossil species. We sampled 40 extant and extinct genera across all major lineages of Acroceridae, which were compared with outgroup taxa from various lower brachyceran families. In all, 81 morphological characters of 60 extant and 10 extinct ingroup species were combined with 7.1 kb of DNA sequences of two nuclear (CAD and 28S rDNA) and two mitochondrial genes (COI and 16S rDNA). Results strongly support the monophyly of Acroceridae, with major clades contained within classified here in five extant subfamilies (Acrocerinae, Cyrtinae stat. rev. , Ogcodinae stat. rev. , Panopinae and Philopotinae) and one extinct subfamily, Archocyrtinae. The evolution of important spider fly traits is discussed, including genitalia and wing venation. The status of the enigmatic Psilodera Gray and Pterodontia Gray as members of the Panopinae is confirmed based on both molecular and morphological data.     

Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds

Using standard phylogenetic methods, it can be hard to resolve the order in which speciation events took place when new lineages evolved in the distant past and within a short time frame. As an example, phylogenies of galliform birds (including well-known species such as chicken, turkey, and quail) usually show low bootstrap support values at short internal branches, reflecting the rapid diversification of these birds in the Eocene. However, given the key role of chicken and related poultry species in agricultural, evolutionary, general biological and disease studies, it is important to know their internal relationships. Recently, insertion patterns of transposable elements such as long and short interspersed nuclear element markers have proved powerful in revealing branching orders of difficult phylogenies. Here we decipher the order of speciation events in a group of 27 galliform species based on insertion events of chicken repeat 1 (CR1) transposable elements. Forty-four CR1 marker loci were identified from the draft sequence of the chicken genome, and from turkey BAC clone sequence, and the presence or absence of markers across species was investigated via electrophoretic size separation of amplification products and subsequent confirmation by DNA sequencing. Thirty markers proved possible to type with electrophoresis of which 20 were phylogenetically informative. The distribution of these repeat elements supported a single homoplasy-free cladogram, which confirmed that megapodes, cracids, New World quail, and guinea fowl form outgroups to Phasianidae and that quails, pheasants, and partridges are each polyphyletic groups. Importantly, we show that chicken is an outgroup to turkey and quail, an observation which does not have significant support from previous DNA sequence- and DNA-DNA hybridization-based trees and has important implications for evolutionary studies based on sequence or karyotype data from galliforms. We discuss the potential and limitations of using a genome-based retrotransposon approach in resolving problematic phylogenies among birds.     

Back-casting sociality in extinct species: new perspectives using mass death assemblages and sex ratios

Despite 150 years of interest in the ecology of dinosaurs, mammoths, proto-hominids and other extinct vertebrates, a general framework to recreate patterns of sociality has been elusive. Based on our recent discovery of a contemporary heterospecific mass death assemblage in the Gobi Desert (Mongolia), we fit predictions about gender-specific associations and group living in extant ungulates to extinct ones. We relied on comparative data on sex-ratio variation and body-size dimorphism, basing analyses on 38 additional mass mortality sites from Asia, Africa, Europe and North America that span 50 million years. Both extant and extinct species died in aggregations with biased adult sex ratios, but the skew (from 1:1) was greater for extinct dimorphic taxa, suggesting that sociality in these extinct species can be predicted from spatial and demographic traits of extant ones. However, extinct rhinos, horses and zebras were inconsistent with predictions about adult sex ratios, which underscores the inherent difficulty in backcasting historic patterns to some monomorphic taxa. These findings shed light not only on the sociality of extinct species but provide a sound, although limited, footing for interpretation of modern death assemblages within the context of the emerging science of taphonomy and palaeobehaviour.     

Reconstructing Dryopteris "semicristata" (Dryopteridaceae): Molecular profiles of tetraploids verify their undiscovered diploid ancestor

? Premise of the study: Discovering missing ancestors is essential to understanding the evolutionary history of biodiversity on Earth. Evidence from extinct species can provide links for reconstructing intricate patterns of reticulate relationships among extant descendents. When fossils are unavailable and other evidence yields competing hypotheses to explain species ancestry, data from proteins and DNA can help resolve conflicts and generate novel perspectives. The identity of a parent shared by two tetraploid species in the cosmopolitan fern genus Dryopteris has remained elusive for more than 50 years. Based on available data, four hypotheses were developed previously, each providing a different resolution to this uncertainty. ? Methods: New molecular evidence from studies of isozymes and restriction site analysis of chloroplast DNA tested the competing hypotheses about the diploid ancestors of these two extant Dryopteris polyploids. ? Key results: The results falsify two of the hypotheses, resolve the uncertainty in the third, and support the fourth. ? Conclusions: Our data validate the prior existence of Dryopteris "semicristata," which was proposed 38 years ago as a diploid progenitor of the allotetraploids D. cristata and D. carthusiana but has never been collected. After developing a phylogeny using the new molecular data, we describe a plausible morphology for D. "semicristata" by extrapolating likely character states from related extant species.     

Mammalian phylogeny: comparison of morphological and molecular results

In an attempt to resolve the "bushy" part at the root of the eutherian tree, 182 nondental morphological characters from 100 species (79 extant and 21 extinct; 98 mammalian and 2 nonmammalian) were analyzed using two maximum-parsimony tree-building algorithms. Parallel analyses of 2,258 pairwise immunodiffusion comparisons with chicken antisera on 101 mammalian species and of amino acid sequence data of alpha and beta hemoglobins and other published protein sequences were also carried out. The morphological and molecular phylogenies agree in depicting the infraclass Eutheria as consisting of five major clades (thus resolving part of the "bush"). Rates of evolution were also found to be similar in the two types of phylogenies.     

Chromosomal phylogeny and evolution of gibbons (Hylobatidae)

Although human and gibbons are classified in the same primate superfamily (Hominoidae), their karyotypes differ by extensive chromosome reshuffling. To date, there is still limited understanding of the events that shaped extant gibbon karyotypes. Further, the phylogeny and evolution of the twelve or more extant gibbon species (lesser apes, Hylobatidae) is poorly understood, and conflicting phylogenies have been published. We present a comprehensive analysis of gibbon chromosome rearrangements and a phylogenetic reconstruction of the four recognized subgenera based on molecular cytogenetics data. We have used two different approaches to interpret our data: (1) a cladistic reconstruction based on the identification of ancestral versus derived chromosome forms observed in extant gibbon species; (2) an approach in which adjacent homologous segments that have been changed by translocations and intra-chromosomal rearrangements are treated as discrete characters in a parsimony analysis (PAUP). The orangutan serves as an "outgroup", since it has a karyotype that is supposed to be most similar to the ancestral form of all humans and apes. Both approaches place the subgenus Bunopithecus as the most basal group of the Hylobatidae, followed by Hylobates, with Symphalangus and Nomascus as the last to diverge. Since most chromosome rearrangements observed in gibbons are either ancestral to all four subgenera or specific for individual species and only a few common derived rearrangements at subsequent branching points have been recorded, all extant gibbons may have diverged within relatively short evolutionary time. In general, chromosomal rearrangements produce changes that should be considered as unique landmarks at the divergence nodes. Thus, molecular cytogenetics could be an important tool to elucidate phylogenies in other species in which speciation may have occurred over very short evolutionary time with not enough genetic (DNA sequence) and other biological divergence to be picked up.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Influence of Evolutionary Allometry on Rates of Morphological Evolution and Disparity in strictly Subterranean Moles (Talpinae,Talpidae, Lipotyphla,Mammalia)

The adaptation to a particular function could directly influence the morphological evolution of an anatomical structure as well as its rates. The humeral morphology of moles (subfamily Talpinae) is highly modified in response to intense burrowing and fully fossorial lifestyle. However, little is known of the evolutionary pathways that marked its diversification in the two highly fossorial moles tribes Talpini and Scalopini. We used two-dimensional landmark-based geometric morphometrics and comparative methods to understand which factors influenced the rates and patterns of the morphological evolution of the humerus in 53 extant and extinct species of the Talpini (22 extant plus 12 extinct) and Scalopini (six extant plus 13 extinct) tribes, for a total of 623 humeri. We first built a synthetic phylogeny of extinct and extant taxa of the subfamily Talpinae based on all the available information from known phylogenies, molecular data, and age ranges of fossil records. We tested for evolutionary allometry by means of multivariate regression of shape on size variables. Evolutionary allometric trajectories exhibited convergence of humeral shape between the two tribes, even when controlling for phylogeny, though a significant differences in the evolutionary rates was found between the two tribes. Talpini, unlike Scalopini, seem to have reached a robust fossorial morphology early during their evolution, and their shape disparity did not change, if it did not decrease, through time. Furthermore, the basal Geotrypus spp. clearly set apart from the other highly fossorial moles, exhibiting a significant acceleration of evolutionary shifts toward higher degree of fossorial adaptation. Our observations support the hypothesis that the evolution of allometry may reflect a biological demand (in this case functional) that constrains the rates of evolution of anatomical structures.     

The phylogeny of Cetartiodactyla: the importance of dense taxon sampling, missing data, and the remarkable promise of cytochrome b to provide reliable species-level phylogenies

We perform Bayesian phylogenetic analyses on cytochrome b sequences from 264 of the 290 extant cetartiodactyl mammals (whales plus even-toed ungulates) and two recently extinct species, the 'Mouse Goat' and the 'Irish Elk'. Previous primary analyses have included only a small portion of the species diversity within Cetartiodactyla, while a complete supertree analysis lacks resolution and branch lengths limiting its utility for comparative studies. The benefits of using a single-gene approach include rapid phylogenetic estimates for a large number of species. However, single-gene phylogenies often differ dramatically from studies involving multiple datasets suggesting that they often are unreliable. However, based on recovery of benchmark clades-clades supported in prior studies based on multiple independent datasets-and recovery of undisputed traditional taxonomic groups, Cytb performs extraordinarily well in resolving cetartiodactyl phylogeny when taxon sampling is dense. Missing data, however, (taxa with partial sequences) can compromise phylogenetic accuracy, suggesting a tradeoff between the benefits of adding taxa and introducing question marks. In the full data, a few species with a short sequences appear misplaced, however, sequence length alone seems a poor predictor of this phenomenon as other taxa with equally short sequences were not conspicuously misplaced. Although we recommend awaiting a better supported phylogeny based on more character data to reconsider classification and taxonomy within Cetartiodactyla, the new phylogenetic hypotheses provided here represent the currently best available tool for comparative species-level studies within this group. Cytb has been sequenced for a large percentage of mammals and appears to be a reliable phylogenetic marker as long as taxon sampling is dense. Therefore, an opportunity exists now to reconstruct detailed phylogenies of most of the major mammalian clades to rapidly provide much needed tools for species-level comparative studies.     

Is a large-scale DNA-based inventory of ancient life possible?

A complete DNA-based inventory of the Earth's present biota using large-scale high-throughput DNA sequencing of signature region(s) (DNA barcoding) is an ambitious proposal rivaling the Human Genome Project. We examine whether this approach will also enable us to assess the past diversity of the earth's biota. To test this, we sequenced the 5' terminus of the mitochondrial cytochrome c oxidase I (COI) gene of individuals belonging to a group of extinct ratite birds, the moa of New Zealand. Moa comprised a large number of taxa that radiated in isolation on this oceanic landmass. Using a phylogenetic approach based on a large data set including protein coding and 12S DNA sequences as well as morphology, we now have precise information about the number of moa species that once existed. We show that each of the moa species detected using this extensive data set has a unique COI barcode(s) and that they all show low levels of within-species COI variation. Consequently, we conclude that COI sequences accurately identify the species discovered using the larger data set. Hence, more generally, this study suggests that DNA barcoding might also help us detect other extinct animal species and that a large-scale inventory of ancient life is possible.     

A refined model of body mass and population density in flightless birds reconciles extreme bimodal population estimates for extinct moa

Flightless birds were once the largest and heaviest terrestrial fauna on many archipelagos around the world. Robust approaches for estimating their population parameters are essential for understanding prehistoric insular ecosystems and extinction processes. Body mass and population density are negatively related for extant flightless bird species, providing a method for quantifying densities and population sizes of extinct flightless species. Here we assemble an updated global data set of body mass and population densities for extant flightless birds and estimate the relationship between these variables. We use generalised least squares models that account for phylogenetic relatedness and incorporate the effects of limiting factors (e.g. habitat suitability) on population density. We demonstrate the applicability of this allometric relationship to extinct species by estimating densities for each of the nine species of moa (Dinornithiformes) and generating a combined spatially explicit map of total moa density across New Zealand. To compare our density estimates with those previously published, we summed individual species' abundances to generate a mean national density of 2.02–9.66 birds km⁻² for low- and high-density scenarios, respectively. Our results reconcile the extreme bimodality of previous estimates (< 2 birds km⁻² and > 10 birds km⁻²) and are comparable to contemporary densities of large herbivorous wild mammals introduced into New Zealand about 150 yr ago. The revised moa density has little effect on the harvest rates required to bring about extinction within 150–200 yr, indicating that rapid extinction was an inevitable response to human hunting, irrespective of the initial population of moa.