The evolution of multiple colour mechanisms is correlated with diversification in sunbirds (Nectariniidae)

How and why certain groups become speciose is a key question in evolutionary biology. Novel traits that enable diversification by opening new ecological niches are likely important mechanisms. However, ornamental traits can also promote diversification by opening up novel sensory niches and thereby creating novel inter-specific interactions. More specifically, ornamental colours may enable more precise and/or easier species recognition, and may act as key innovations by increasing the number of species-specific patterns and promoting diversification. While the influence of colouration on diversification is well-studied, the influence of the mechanisms that produce those colours (e.g. pigmentary, nanostructural) is less so, even though the ontogeny and evolution of these mechanisms differ. We estimated a new phylogenetic tree for 121 sunbird species and combined colour data of 106 species with a range of phylogenetic tools to test the hypothesis that the evolution of novel colour mechanisms increases diversification in sunbirds, one of the most colourful bird clades. Results suggest that (1) the evolution of novel colour mechanisms expands the visual sensory niche, increasing the number of achievable colours. (2) Structural colouration diverges more readily across the body than pigment-based colouration, enabling an increase in colour complexity. (3) Novel colour mechanisms might minimize trade-offs between natural and sexual selection such that colour can function both as camouflage and conspicuous signal. (4) Despite structural colours being more colourful and mobile, only melanin-based colouration is positively correlated with net diversification. Together, these findings explain why colour distances increase with increasing number of sympatric species, even though packing of colour space would predict otherwise.

Ancestral reconstruction of the MotA stator subunit reveals that conserved residues far from the pore are required to drive flagellar motility

The bacterial flagellar motor (BFM) is a rotary nanomachine powered by the translocation of ions across the inner membrane through the stator complex. The stator complex consists of two membrane proteins: MotA and MotB (in H⁺-powered motors), or PomA and PomB (in Na⁺-powered motors). In this study, we used ancestral sequence reconstruction (ASR) to probe which residues of MotA correlate with function and may have been conserved to preserve motor function. We reconstructed 10 ancestral sequences of MotA and found four of them were motile in combination with contemporary Escherichia coli MotB and in combination with our previously published functional ancestral MotBs. Sequence comparison between wild-type (WT) E. coli MotA and MotA-ASRs revealed 30 critical residues across multiple domains of MotA that were conserved among all motile stator units. These conserved residues included pore-facing, cytoplasm-facing, and MotA-MotA intermolecular facing sites. Overall, this work demonstrates the role of ASR in assessing conserved variable residues in a subunit of a molecular complex.

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNA^Ile, ND1, tRNA^Leu2, 16S, tRNA^Val, 12S, tRNA^Phe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae)

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4-10x, Veronica: 12-18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.

Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales

Evolutionary models account for either population- or species-level processes but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth–death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology.[Caninae; fossilized birth–death; molecular clock; multispecies coalescent; phylogenetics; species trees.]

Flagellar export apparatus and ATP synthetase: Homology evidenced by synteny predating the Last Universal Common Ancestor

We report evidence further supporting homology between proteins in the F₁ F_O -ATP synthetase and the bacterial flagellar motor (BFM). BFM proteins FliH, FliI, and FliJ have been hypothesized to be homologous to F_O -b + F₁ -δ, F₁ -α/β, and F₁ -γ, with similar structure and interactions. We conduct a further test by constructing a gene order dataset, examining the order of fliH, fliI, and fliJ genes across the phylogenetic breadth of flagellar and nonflagellar type 3 secretion systems, and comparing this to published surveys of gene order in the F₁ F_O -ATP synthetase, its N-ATPase relatives, and the bacterial/archaeal V- and A-type ATPases. Strikingly, the fliHIJ gene order was deeply conserved, with the few exceptions appearing derived, and exactly matching the widely conserved F-ATPase gene order atpFHAG, coding for subunits b-δ-α-γ. The V/A-type ATPases have a similar conserved gene order. Our results confirm homology between these systems, and suggest a rare case of synteny conserved over billions of years, predating the Last Universal Common Ancestor (LUCA).

Trait-dependent dispersal in rails (Aves: Rallidae): Historical biogeography of a cosmopolitan bird clade

The ability of lineages to disperse over evolutionary timescales may be influenced by the gain or loss of traits after adaptation to new ecological conditions. For example, rails (Aves: Rallidae) have many cases of flightless insular endemic species that presumably evolved after flying ancestors dispersed over large ocean barriers and became isolated. Nonetheless, the details of how flying and its loss have influenced the clade's historical biogeography are unknown, as is the importance of other predictors of dispersal such as the geographic distance between regions. Here, we used a dated phylogeny of 158 species of rails to compare trait-dependent and trait-independent biogeography models in BioGeoBEARS. We evaluated a probabilistic historical biogeographical model that allows geographic range and flight to co-evolve and influence dispersal ability on a phylogeny. The best-fitting dispersal model was a trait-dependent dispersal (DEC + j + x + t₂₁ + m₁) that accrued 85.2% of the corrected Akaike Information Criterion (AICc) model weight. The distance-dependence parameter, x was estimated at -0.54, ranging from -0.49 to -0.65 across models, suggesting that a doubling of dispersal distance results in an approximately 31% decrease in dispersal rate (2^-0.54 = 0.69). The estimated rate of loss of flight (t₂₁) was similar across all models (~0.029 loss events per lineage per million years). The multiplier on dispersal rate when a lineage is non-flying, m₁, is estimated to be 0.38 under this model. Surprisingly, the estimate of m₁ was not 0.0, probably because the loss of flight is so common in the rails that entire clades of flightless species are found in the data, forcing the model to attribute some dispersal to flightless lineages. These results indicate that long-distance dispersal over macroevolutionary timespans can be modelled, rather than simply attributed to chance, allowing support for different hypotheses to be quantified and limitations to be identified. Overall, by combining new analytical methods with a comprehensive phylogeny, we use a quantitative framework to show how traits influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Ancestral Sequence Reconstructions of MotB Are Proton-Motile and Require MotA for Motility

The bacterial flagellar motor (BFM) is a nanomachine that rotates the flagellum to propel many known bacteria. The BFM is powered by ion transit across the cell membrane through the stator complex, a membrane protein. Different bacteria use various ions to run their BFM, but the majority of BFMs are powered by either proton (H+) or sodium (Na+) ions. The transmembrane (TM) domain of the B-subunit of the stator complex is crucial for ion selectivity, as it forms the ion channel in complex with TM3 and TM4 of the A-subunit. In this study, we reconstructed and engineered thirteen ancestral sequences of the stator B-subunit to evaluate the functional properties and ionic power source of the stator proteins at reconstruction nodes to evaluate the potential of ancestral sequence reconstruction (ASR) methods for stator engineering and to test specific motifs previously hypothesized to be involved in ion-selectivity. We found that all thirteen of our reconstructed ancient B-subunit proteins could assemble into functional stator complexes in combination with the contemporary Escherichia coli MotA-subunit to restore motility in stator deleted E. coli strains. The flagellar rotation of the thirteen ancestral MotBs was found to be Na+ independent which suggested that the F30/Y30 residue was not significantly correlated with sodium/proton phenotype, in contrast to what we had reported previously. Additionally, four among the thirteen reconstructed B-subunits were compatible with the A-subunit of Aquifex aeolicus and able to function in a sodium-independent manner. Overall, this work demonstrates the use of ancestral reconstruction to generate novel stators and quantify which residues are correlated with which ionic power source.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Statistical Comparison of Trait-Dependent Biogeographical Models Indicates That Podocarpaceae Dispersal Is Influenced by Both Seed Cone Traits and Geographical Distance

The ability of lineages to disperse long distances over evolutionary timescales may be influenced by the gain or loss of traits adapted to enhance local, ecological dispersal. For example, some species in the southern conifer family Podocarpaceae have fleshy cones that encourage bird dispersal, but it is unknown how this trait has influenced the clade’s historical biogeography, or its importance compared with other predictors of dispersal such as the geographic distance between regions. We answer these questions quantitatively by using a dated phylogeny of 197 species of southern conifers (Podocarpaceae and their sister family Araucariaceae) to statistically compare standard, trait-independent biogeography models with new BioGeoBEARS models where an evolving trait can influence dispersal probability, and trait history, biogeographical history, and model parameters are jointly inferred. We validate the method with simulation-inference experiments. Comparing all models, those that include trait-dependent dispersal accrue 87.5% of the corrected Akaike Information Criterion (AICc) model weight. Averaged across all models, lineages with nonfleshy cones had a dispersal probability multiplier of 0.49 compared with lineages with fleshy cones. Distance is included as a predictor of dispersal in all credible models (100% model weight). However, models with changing geography earned only 22.0% of the model weight, and models submerging New Caledonia/New Zealand earned only 0.01%. The importance of traits and distance suggests that long-distance dispersal over macroevolutionary timespans should not be thought of as a highly unpredictable chance event. Instead, long-distance dispersal can be modeled, allowing statistical model comparison to quantify support for different hypotheses.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae)

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Taxon cycle predictions supported by model-based inference in Indo-Pacific trap-jaw ants (Hymenoptera: Formicidae: Odontomachus)

Nonequilibrium dynamics and non-neutral processes, such as trait-dependent dispersal, are often missing from quantitative island biogeography models despite their potential explanatory value. One of the most influential nonequilibrium models is the taxon cycle, but it has been difficult to test its validity as a general biogeographical framework. Here, we test predictions of the taxon cycle model using six expected phylogenetic patterns and a time-calibrated phylogeny of Indo-Pacific Odontomachus (Hymenoptera: Formicidae: Ponerinae), one of the ant genera that E.O. Wilson used when first proposing the hypothesis. We used model-based inference and a newly developed trait-dependent dispersal model to jointly estimate ancestral biogeography, ecology (habitat preferences for forest interiors, vs. "marginal" habitats, such as savannahs, shorelines, disturbed areas) and the linkage between ecology and dispersal rates. We found strong evidence that habitat shifts from forest interior to open and disturbed habitats increased macroevolutionary dispersal rate. In addition, lineages occupying open and disturbed habitats can give rise to both island endemics re-occupying only forest interiors and taxa that re-expand geographical ranges. The phylogenetic predictions outlined in this study can be used in future work to evaluate the relative weights of neutral (e.g., geographical distance and area) and non-neutral (e.g., trait-dependent dispersal) processes in historical biogeography and community ecology.

Diversification across biomes in a continental lizard radiation

Ecological opportunity is a powerful driver of evolutionary diversification, and predicts rapid lineage and phenotypic diversification following colonization of competitor-free habitats. Alternatively, topographic or environmental heterogeneity could be key to generating and sustaining diversity. We explore these hypotheses in a widespread lineage of Australian lizards: the Gehyra variegata group. This clade occurs across two biomes: the Australian monsoonal tropics (AMT), where it overlaps a separate, larger bodied clade of Gehyra and is largely restricted to rocks; and in the larger Australian arid zone (AAZ) where it has no congeners and occupies trees and rocks. New phylogenomic data and coalescent analyses of AAZ taxa resolve lineages and their relationships and reveal high diversity in the western AAZ (Pilbara region). The AMT and AAZ radiations represent separate radiations with no difference in speciation rates. Most taxa occur on rocks, with small geographic ranges relative to widespread generalist taxa across the vast central AAZ. Rock-dwelling and generalist taxa differ morphologically, but only the lineage-poor central AAZ taxa have accelerated evolution. This accords with increasing evidence that lineage and morphological diversity are poorly correlated, and suggests environmental heterogeneity and refugial dynamics have been more important than ecological release in elevating lineage diversity.

Including autapomorphies is important for paleontological tip-dating with clocklike data, but not with non-clock data

Tip-dating, where fossils are included as dated terminal taxa in Bayesian dating inference, is an increasingly popular method. Data for these studies often come from morphological character matrices originally developed for non-dated, and usually parsimony, analyses. In parsimony, only shared derived characters (synapomorphies) provide grouping information, so many character matrices have an ascertainment bias: they omit autapomorphies (unique derived character states), which are considered uninformative. There has been no study of the effect of this ascertainment bias in tip-dating, but autapomorphies can be informative in model-based inference. We expected that excluding autapomorphies would shorten the morphological branchlengths of terminal branches, and thus bias downwards the time branchlengths inferred in tip-dating. We tested for this effect using a matrix for Carboniferous-Permian eureptiles where all autapomorphies had been deliberately coded. Surprisingly, date estimates are virtually unchanged when autapomorphies are excluded, although we find large changes in morphological rate estimates and small effects on topological and dating confidence. We hypothesized that the puzzling lack of effect on dating was caused by the non-clock nature of the eureptile data. We confirm this explanation by simulating strict clock and non-clock datasets, showing that autapomorphy exclusion biases dating only for the clocklike case. A theoretical solution to ascertainment bias is computing the ascertainment bias correction (Mk_parsinf), but we explore this correction in detail, and show that it is computationally impractical for typical datasets with many character states and taxa. Therefore we recommend that palaeontologists collect autapomorphies whenever possible when assembling character matrices.

Inferring node dates from tip dates in fossil Canidae: the importance of tree priors

Tip-dating methods are becoming popular alternatives to traditional node calibration approaches for building time-scaled phylogenetic trees, but questions remain about their application to empirical datasets. We compared the performance of the most popular methods against a dated tree of fossil Canidae derived from previously published monographs. Using a canid morphology dataset, we performed tip-dating using BEAST v. 2.1.3 and MrBayes v. 3.2.5. We find that for key nodes (Canis, approx. 3.2 Ma, Caninae approx. 11.7 Ma) a non-mechanistic model using a uniform tree prior produces estimates that are unrealistically old (27.5, 38.9 Ma). Mechanistic models (incorporating lineage birth, death and sampling rates) estimate ages that are closely in line with prior research. We provide a discussion of these two families of models (mechanistic versus non-mechanistic) and their applicability to fossil datasets.

Recent origin and rapid speciation of Neotropical orchids in the world's richest plant biodiversity hotspot

The Andean mountains of South America are the most species-rich biodiversity hotspot worldwide with c. 15% of the world's plant species, in only 1% of the world's land surface. Orchids are a key element of the Andean flora, and one of the most prominent components of the Neotropical epiphyte diversity, yet very little is known about their origin and diversification. We address this knowledge gap by inferring the biogeographical history and diversification dynamics of the two largest Neotropical orchid groups (Cymbidieae and Pleurothallidinae), using two unparalleled, densely sampled orchid phylogenies (including more than 400 newly generated DNA sequences), comparative phylogenetic methods, geological and biological datasets. We find that the majority of Andean orchid lineages only originated in the last 20-15 million yr. Andean lineages are derived from lowland Amazonian ancestors, with additional contributions from Central America and the Antilles. Species diversification is correlated with Andean orogeny, and multiple migrations and recolonizations across the Andes indicate that mountains do not constrain orchid dispersal over long timescales. Our study sheds new light on the timing and geography of a major Neotropical diversification, and suggests that mountain uplift promotes species diversification across all elevational zones.

Empirical and Bayesian approaches to fossil-only divergence times: A study across three reptile clades

Estimating divergence times on phylogenies is critical in paleontological and neontological studies. Chronostratigraphically-constrained fossils are the only direct evidence of absolute timing of species divergence. Strict temporal calibration of fossil-only phylogenies provides minimum divergence estimates, and various methods have been proposed to estimate divergences beyond these minimum values. We explore the utility of simultaneous estimation of tree topology and divergence times using BEAST tip-dating on datasets consisting only of fossils by using relaxed morphological clocks and birth-death tree priors that include serial sampling (BDSS) at a constant rate through time. We compare BEAST results to those from the traditional maximum parsimony (MP) and undated Bayesian inference (BI) methods. Three overlapping datasets were used that span 250 million years of archosauromorph evolution leading to crocodylians. The first dataset focuses on early Sauria (31 taxa, 240 chars.), the second on early Archosauria (76 taxa, 400 chars.) and the third on Crocodyliformes (101 taxa, 340 chars.). For each dataset three time-calibrated trees (timetrees) were calculated: a minimum-age timetree with node ages based on earliest occurrences in the fossil record; a 'smoothed' timetree using a range of time added to the root that is then averaged over zero-length internodes; and a tip-dated timetree. Comparisons within datasets show that the smoothed and tip-dated timetrees provide similar estimates. Only near the root node do BEAST estimates fall outside the smoothed timetree range. The BEAST model is not able to overcome limited sampling to correctly estimate divergences considerably older than sampled fossil occurrence dates. Conversely, the smoothed timetrees consistently provide node-ages far older than the strict dates or BEAST estimates for morphologically conservative sister-taxa when they sit on long ghost lineages. In this latter case, the relaxed-clock model appears to be correctly moderating the node-age estimate based on the limited morphological divergence. Topologies are generally similar across analyses, but BEAST trees for crocodyliforms differ when clades are deeply nested but contain very old taxa. It appears that the constant-rate sampling assumption of the BDSS tree prior influences topology inference by disfavoring long, unsampled branches.

Crown group Oxyphotobacteria postdate the rise of oxygen

The rise of oxygen ca. 2.3 billion years ago (Ga) is the most distinct environmental transition in Earth history. This event was enabled by the evolution of oxygenic photosynthesis in the ancestors of Cyanobacteria. However, long-standing questions concern the evolutionary timing of this metabolism, with conflicting answers spanning more than one billion years. Recently, knowledge of the Cyanobacteria phylum has expanded with the discovery of non-photosynthetic members, including a closely related sister group termed Melainabacteria, with the known oxygenic phototrophs restricted to a clade recently designated Oxyphotobacteria. By integrating genomic data from the Melainabacteria, cross-calibrated Bayesian relaxed molecular clock analyses show that crown group Oxyphotobacteria evolved ca. 2.0 billion years ago (Ga), well after the rise of atmospheric dioxygen. We further estimate the divergence between Oxyphotobacteria and Melainabacteria ca. 2.5-2.6 Ga, which-if oxygenic photosynthesis is an evolutionary synapomorphy of the Oxyphotobacteria-marks an upper limit for the origin of oxygenic photosynthesis. Together, these results are consistent with the hypothesis that oxygenic photosynthesis evolved relatively close in time to the rise of oxygen.

Topology, divergence dates, and macroevolutionary inferences vary between different tip-dating approaches applied to fossil theropods (Dinosauria)

Dated phylogenies of fossil taxa allow palaeobiologists to estimate the timing of major divergences and placement of extinct lineages, and to test macroevolutionary hypotheses. Recently developed Bayesian 'tip-dating' methods simultaneously infer and date the branching relationships among fossil taxa, and infer putative ancestral relationships. Using a previously published dataset for extinct theropod dinosaurs, we contrast the dated relationships inferred by several tip-dating approaches and evaluate potential downstream effects on phylogenetic comparative methods. We also compare tip-dating analyses to maximum-parsimony trees time-scaled via alternative a posteriori approaches including via the probabilistic cal3 method. Among tip-dating analyses, we find opposing but strongly supported relationships, despite similarity in inferred ancestors. Overall, tip-dating methods infer divergence dates often millions (or tens of millions) of years older than the earliest stratigraphic appearance of that clade. Model-comparison analyses of the pattern of body-size evolution found that the support for evolutionary mode can vary across and between tree samples from cal3 and tip-dating approaches. These differences suggest that model and software choice in dating analyses can have a substantial impact on the dated phylogenies obtained and broader evolutionary inferences.

Bayesian analysis of a morphological supermatrix sheds light on controversial fossil hominin relationships

The phylogenetic relationships of several hominin species remain controversial. Two methodological issues contribute to the uncertainty-use of partial, inconsistent datasets and reliance on phylogenetic methods that are ill-suited to testing competing hypotheses. Here, we report a study designed to overcome these issues. We first compiled a supermatrix of craniodental characters for all widely accepted hominin species. We then took advantage of recently developed Bayesian methods for building trees of serially sampled tips to test among hypotheses that have been put forward in three of the most important current debates in hominin phylogenetics--the relationship between Australopithecus sediba and Homo, the taxonomic status of the Dmanisi hominins, and the place of the so-called hobbit fossils from Flores, Indonesia, in the hominin tree. Based on our results, several published hypotheses can be statistically rejected. For example, the data do not support the claim that Dmanisi hominins and all other early Homo specimens represent a single species, nor that the hobbit fossils are the remains of small-bodied modern humans, one of whom had Down syndrome. More broadly, our study provides a new baseline dataset for future work on hominin phylogeny and illustrates the promise of Bayesian approaches for understanding hominin phylogenetic relationships.

Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial-interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia - fossil records, species distribution models and phylogeographic surveys - in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Bayesian analysis of congruence of core genes in Prochlorococcus and Synechococcus and implications on horizontal gene transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. "Universal proteins" have been used to infer the organismal phylogeny, but have been criticized as being only the "tree of one percent." Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees - a "treecloud" - for each, as well as for a concatenated dataset based on putative "universal proteins." We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the "tree of one percent" is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.

Genome duplication and multiple evolutionary origins of complex migratory behavior in Salmonidae

Multiple rounds of whole genome duplication have repeatedly marked the evolution of vertebrates, and correlate strongly with morphological innovation. However, less is known about the behavioral, physiological and ecological consequences of genome duplication, and whether these events coincide with major transitions in vertebrate complexity. The complex behavior of anadromy - where adult fishes migrate up rivers from the sea to their natal site to spawn - is well known in salmonid fishes. Some hypotheses suggest that migratory behavior evolved as a consequence of an ancestral genome duplication event, which permitted salinity tolerance and osmoregulatory plasticity. Here we test whether anadromy evolved multiple times within salmonids, and whether genome duplication coincided with the evolution of anadromy. We present a method that uses ancestral character simulation data to plot the frequency of character transitions over a time calibrated phylogenetic tree to provide estimates of the absolute timing of character state transitions. Furthermore, we incorporate extinct and extant taxa to improve on previous estimates of divergence times. We present the first phylogenetic evidence indicating that anadromy evolved at least twice from freshwater salmonid ancestors. Results suggest that genome duplication did not coincide in time with changes in migratory behavior, but preceded a transition to anadromy by 55-50 million years. Our study represents the first attempt to estimate the absolute timing of a complex behavioral trait in relation to a genome duplication event.

Bayesian analysis of biogeography when the number of areas is large

Historical biogeography is increasingly studied from an explicitly statistical perspective, using stochastic models to describe the evolution of species range as a continuous-time Markov process of dispersal between and extinction within a set of discrete geographic areas. The main constraint of these methods is the computational limit on the number of areas that can be specified. We propose a Bayesian approach for inferring biogeographic history that extends the application of biogeographic models to the analysis of more realistic problems that involve a large number of areas. Our solution is based on a "data-augmentation" approach, in which we first populate the tree with a history of biogeographic events that is consistent with the observed species ranges at the tips of the tree. We then calculate the likelihood of a given history by adopting a mechanistic interpretation of the instantaneous-rate matrix, which specifies both the exponential waiting times between biogeographic events and the relative probabilities of each biogeographic change. We develop this approach in a Bayesian framework, marginalizing over all possible biogeographic histories using Markov chain Monte Carlo (MCMC). Besides dramatically increasing the number of areas that can be accommodated in a biogeographic analysis, our method allows the parameters of a given biogeographic model to be estimated and different biogeographic models to be objectively compared. Our approach is implemented in the program, BayArea.

Whole mitochondrial genome sequencing of domestic horses reveals incorporation of extensive wild horse diversity during domestication

BACKGROUND

DNA target enrichment by micro-array capture combined with high throughput sequencing technologies provides the possibility to obtain large amounts of sequence data (e.g. whole mitochondrial DNA genomes) from multiple individuals at relatively low costs. Previously, whole mitochondrial genome data for domestic horses (Equus caballus) were limited to only a few specimens and only short parts of the mtDNA genome (especially the hypervariable region) were investigated for larger sample sets.

RESULTS

In this study we investigated whole mitochondrial genomes of 59 domestic horses from 44 breeds and a single Przewalski horse (Equus przewalski) using a recently described multiplex micro-array capture approach. We found 473 variable positions within the domestic horses, 292 of which are parsimony-informative, providing a well resolved phylogenetic tree. Our divergence time estimate suggests that the mitochondrial genomes of modern horse breeds shared a common ancestor around 93,000 years ago and no later than 38,000 years ago. A Bayesian skyline plot (BSP) reveals a significant population expansion beginning 6,000-8,000 years ago with an ongoing exponential growth until the present, similar to other domestic animal species. Our data further suggest that a large sample of wild horse diversity was incorporated into the domestic population; specifically, at least 46 of the mtDNA lineages observed in domestic horses (73%) already existed before the beginning of domestication about 5,000 years ago.

CONCLUSIONS

Our study provides a window into the maternal origins of extant domestic horses and confirms that modern domestic breeds present a wide sample of the mtDNA diversity found in ancestral, now extinct, wild horse populations. The data obtained allow us to detect a population expansion event coinciding with the beginning of domestication and to estimate both the minimum number of female horses incorporated into the domestic gene pool and the time depth of the domestic horse mtDNA gene pool.

Biological design in science classrooms

Although evolutionary biology is replete with explanations for complex biological structures, scientists concerned about evolution education have been forced to confront "intelligent design" (ID), which rejects a natural origin for biological complexity. The content of ID is a subset of the claims made by the older "creation science" movement. Both creationist views contend that highly complex biological adaptations and even organisms categorically cannot result from natural causes but require a supernatural creative agent. Historically, ID arose from efforts to produce a form of creationism that would be less vulnerable to legal challenges and that would not overtly rely upon biblical literalism. Scientists do not use ID to explain nature, but because it has support from outside the scientific community, ID is nonetheless contributing substantially to a long-standing assault on the integrity of science education.

From The Origin of Species to the origin of bacterial flagella

In the recent Dover trial, and elsewhere, the 'Intelligent Design' movement has championed the bacterial flagellum as an irreducibly complex system that, it is claimed, could not have evolved through natural selection. Here we explore the arguments in favour of viewing bacterial flagella as evolved, rather than designed, entities. We dismiss the need for any great conceptual leaps in creating a model of flagellar evolution and speculate as to how an experimental programme focused on this topic might look.