Branch length estimation and divergence dating: estimates of error in Bayesian and maximum likelihood frameworks

Molecular dating of the blood pigment hemocyanin provides new insight into the origin of animals

The Neoproterozoic included changes in oceanic redox conditions, the configuration of continents and climate, extreme ice ages (Sturtian and Marinoan), and the rise of complex life forms. A much-debated topic in geobiology concerns the influence of atmospheric oxygenation on Earth and the origin and diversification of animal lineages, with the most widely popularized hypotheses relying on causal links between oxygen levels and the rise of animals. The vast majority of extant animals use aerobic metabolism for growth and homeostasis; hence, the binding and transportation of oxygen represent a vital physiological task. Considering the blood pigment hemocyanin (Hc) is present in sponges and ctenophores, and likely to be present in the common ancestor of animals, we investigated the evolution and date of Hc emergence using bioinformatics approaches on both transcriptomic and genomic data. Bayesian molecular dating suggested that the ancestral animal Hc gene arose approximately 881 Ma during the Tonian Period (1000-720 Ma), prior to the extreme glaciation events of the Cryogenian Period (720-635 Ma). This result is corroborated by a recently discovered fossil of a putative sponge ~890 Ma and modern molecular dating for the origin of metazoans of ~1,000-650 Ma (but does contradict previous inferences regarding the origin of Hc ~700-600 Ma). Our data reveal that crown-group animals already possessed hemocyanin-like blood pigments, which may have enhanced the oxygen-carrying capacity of these animals in hypoxic environments at that time or acted in the transport of hormones, detoxification of heavy metals, and immunity pathways.

Complex co-evolutionary relationships between cicadas and their symbionts

Previous evidence suggests that cicadas lacking Hodgkinia may harbour the yeast-like fungal symbionts (YLS). Here, we reinforce an earlier conclusion that the pathogenic ancestor of YLS independently infected different cicada lineages instead of the common ancestor of Cicadidae. Five independent replacement events in the loss of Hodgkinia/acquisition of YLS and seven other replacement events of YLS (from an Ophiocordyceps fungus to another Ophiocordyceps fungus) are hypothesised to have occurred within the sampled cicada taxa. The divergence time of YLS lineages was later than that of corresponding cicada lineages. The rapid shift of diversification rates of YLS and related cicada-parasitizing Ophiocordyceps began at approximately 32.94 Ma, and the diversification rate reached the highest value at approximately 24.82 Ma, which corresponds to the cooling climate changes at the Eocene-Oligocene boundary and the Oligocene-Miocene transition respectively. Combined with related acquisition/replacement events of YLS occurred during the cooling-climate periods, we hypothesise that the cooling-climate changes impacted the interactions between cicadas and related Ophiocordyceps, which coupled with the unusual life cycle and the differentiation of cicadas may finally led to the diversification of YLS in Cicadidae. Our results contribute to a better understanding of the evolutionary transition of YLS from entomopathogenic fungi in insects.

Evolution of a key enzyme of aerobic metabolism reveals Proterozoic functional subunit duplication events and an ancient origin of animals

The biological toolkits for aerobic respiration were critical for the rise and diversification of early animals. Aerobic life forms generate ATP through the oxidation of organic molecules in a process known as Krebs' Cycle, where the enzyme isocitrate dehydrogenase (IDH) regulates the cycle's turnover rate. Evolutionary reconstructions and molecular dating of proteins related to oxidative metabolism, such as IDH, can therefore provide an estimate of when the diversification of major taxa occurred, and their coevolution with the oxidative state of oceans and atmosphere. To establish the evolutionary history and divergence time of NAD-dependent IDH, we examined transcriptomic data from 195 eukaryotes (mostly animals). We demonstrate that two duplication events occurred in the evolutionary history of NAD-IDH, one in the ancestor of eukaryotes approximately at 1967 Ma, and another at 1629 Ma, both in the Paleoproterozoic Era. Moreover, NAD-IDH regulatory subunits β and γ are exclusive to metazoans, arising in the Mesoproterozoic. Our results therefore support the concept of an ''earlier-than-Tonian'' diversification of eukaryotes and the pre-Cryogenian emergence of a metazoan IDH enzyme.

Advantages of an easy-to-use DNA extraction method for minimal-destructive analysis of collection specimens

Here we present and justify an approach for minimal-destructive DNA extraction from historic insect specimens for next generation sequencing applications. An increasing number of studies use insects from museum collections for biodiversity research. However, the availability of specimens for molecular analyses has been limited by the degraded nature of the DNA gained from century-old museum material and the consumptive nature of most DNA extraction procedures. The method described in this manuscript enabled us to successfully extract DNA from specimens as old as 241 years using a minimal-destructive approach. The direct comparison of the DNeasy extraction Kit and the Monarch® PCR & DNA Clean-up Kit showed a significant increase of 17.3-fold higher DNA yield extracted with the Monarch Oligo protocol on average. By using an extraction protocol originally designed for oligonucleotide clean-up, we were able to combine overcoming the restrictions by target fragment size and strand state, with minimising time consumption and labour-intensity. The type specimens used for the minimal-destructive DNA extraction exhibited no significant external change or post-extraction damage, while sufficient DNA was retrieved for analyses.

Maximize Resolution or Minimize Error? Using Genotyping-By-Sequencing to Investigate the Recent Diversification of Helianthemum (Cistaceae)

A robust phylogenetic framework, in terms of extensive geographical and taxonomic sampling, well-resolved species relationships and high certainty of tree topologies and branch length estimations, is critical in the study of macroevolutionary patterns. Whereas Sanger sequencing-based methods usually recover insufficient phylogenetic signal, especially in recently diversified lineages, reduced-representation sequencing methods tend to provide well-supported phylogenetic relationships, but usually entail remarkable bioinformatic challenges due to the inherent trade-off between the number of SNPs and the magnitude of associated error rates. The genus Helianthemum (Cistaceae) is a species-rich and taxonomically complex Palearctic group of plants that diversified mainly since the Upper Miocene. It is a challenging case study since previous attempts using Sanger sequencing were unable to resolve the intrageneric phylogenetic relationships. Aiming to obtain a robust phylogenetic reconstruction based on genotyping-by-sequencing (GBS), we established a rigorous methodological workflow in which we i) explored how variable settings during dataset assembly have an impact on error rates and on the degree of resolution under concatenation and coalescent approaches, ii) assessed the effect of two extreme parameter configurations (minimizing error rates vs. maximizing phylogenetic resolution) on tree topology and branch lengths, and iii) evaluated the effects of these two configurations on estimates of divergence times and diversification rates. Our analyses produced highly supported topologically congruent phylogenetic trees for both configurations. However, minimizing error rates did produce more reliable branch lengths, critically affecting the accuracy of downstream analyses (i.e. divergence times and diversification rates). In addition to recommending a revision of intrageneric systematics, our results enabled us to identify three highly diversified lineages in Helianthemum in contrasting geographical areas and ecological conditions, which started radiating in the Upper Miocene.

Near-Full-Length Genetic Characterization of a Novel HIV-1 Unique Recombinant with Similarities to A1, CRF01_AE, and CRFO2_AG Viruses in Yaoundé, Cameroon

Variations in the HIV genome influence HIV/AIDS epidemiology. We report here a novel HIV-1 unique recombinant form (URF) isolated from an HIV-infected female (NACMR092) in Cameroon, based on the analyses of near-full-length viral genome (partial gag, full-length pol, env, tat, rev, vif, vpr, vpu, and nef genes, and partial 3'-long terminal repeat). Phylogeny, recombination breakpoints, and recombination map analyses showed that NACMR092 was infected with a mosaic URF that had eight breakpoints (two in gag, one in pol, one in vpr, two in env, and two in the nef regions), nine subgenomic regions, and included fragments that had important similarities with HIV-1 subtypes A1, CRF02_AG, and CRF01_AE. This novel mosaic URF underscores complex recombination events occurring between HIV-1 subtypes circulating in Cameroon. Continued monitoring and detection of such recombinants and accurate classification of HIV genotype is important for tracking viral molecular epidemiology and antigenic diversity.

Undersampling Genomes has Biased Time and Rate Estimates Throughout the Tree of Life

Genomic data drive evolutionary research on the relationships and timescale of life but the genomes of most species remain poorly sampled. Phylogenetic trees can be reconstructed reliably using small data sets and the same has been assumed for the estimation of divergence time with molecular clocks. However, we show here that undersampling of molecular data results in a bias expressed as disproportionately shorter branch lengths and underestimated divergence times in the youngest nodes and branches, termed the small sample artifact. In turn, this leads to increasing speciation and diversification rates towards the present. Any evolutionary analyses derived from these biased branch lengths and speciation rates will be similarly biased. The widely used timetrees of the major species-rich studies of amphibians, birds, mammals, and squamate reptiles are all data-poor and show upswings in diversification rate, suggesting that their results were biased by undersampling. Our results show that greater sampling of genomes is needed for accurate time and rate estimation, which are basic data used in ecological and evolutionary research.

A Phylogenomic Framework and Divergence History of Cephalochordata Amphioxus

Amphioxus, or cephalochordates, are often used as the living invertebrate proxy of vertebrate ancestors and are widely used as evolutionary biology models of chordates. However, their phylogeny, divergence history, and speciation characteristics remain poorly understood, and phylogenomic studies to explore these problems lacking entirely from the literature. Here, we determined a new transcriptome of Branchiostoma japonicum. Combined with mass sequences of all other 18 species, a 19-way phylogeny was constructed via multiple methods (ML, BI, PhyloBayes, and ASTRAL), consistently supporting a phylogeny of [(B. belcheri + B. japonicum) + (B. lanceolatum + B. floridae) + Asymmetron lucayanum] in amphioxus. Congruent phylogenetic signals were found across mitochondrial genes, 12S RNA, and complete mitochondrial genomes according to previous reports, indicating that 12S RNA may have potential as a molecular marker for phylogenetic analysis in amphioxus. Molecular dating analysis indicated a radiation of the cephalochordates during the Cretaceous (∼104-61 million years ago), supporting an association between the diversification and speciation of cephalochordates with continental drift and associated changes in their respective habitats during this time. The identified functional enrichment analysis for species-specific domains indicated that their function mainly involves immune response, apoptosis, and lipid metabolism and utilization, signaling that pathogens and changes of energy requirements are an important driving force for amphioxus speciation. This study represents the first large-scale phylogenomic analysis of most major amphioxus genera based on phylogenomic data, providing a new perspective on both phylogeny and divergence speciation of cephalochordates.

Impact of the tree prior on estimating clock rates during epidemic outbreaks

Genetic sequencing data of pathogens allow one to quantify the evolutionary rate together with epidemiological dynamics using Bayesian phylodynamic methods. Such tools are particularly useful for obtaining a timely understanding of newly emerging epidemic outbreaks. During the West African Ebola virus disease epidemic, an unusually high evolutionary rate was initially estimated, promoting discussions regarding the potential danger of the strain quickly evolving into an even more dangerous virus. We show here that such high evolutionary rates are not necessarily real but can stem from methodological biases in the analyses. While most analyses of epidemic outbreak data are performed such that these biases may be present, we suggest a solution to overcome these biases in the future.

Bayesian phylogenetics aims at estimating phylogenetic trees together with evolutionary and population dynamic parameters based on genetic sequences. It has been noted that the clock rate, one of the evolutionary parameters, decreases with an increase in the sampling period of sequences. In particular, clock rates of epidemic outbreaks are often estimated to be higher compared with the long-term clock rate. Purifying selection has been suggested as a biological factor that contributes to this phenomenon, since it purges slightly deleterious mutations from a population over time. However, other factors such as methodological biases may also play a role and make a biological interpretation of results difficult. In this paper, we identify methodological biases originating from the choice of tree prior, that is, the model specifying epidemiological dynamics. With a simulation study we demonstrate that a misspecification of the tree prior can upwardly bias the inferred clock rate and that the interplay of the different models involved in the inference can be complex and nonintuitive. We also show that the choice of tree prior can influence the inference of clock rate on real-world Ebola virus (EBOV) datasets. While commonly used tree priors result in very high clock-rate estimates for sequences from the initial phase of the epidemic in Sierra Leone, tree priors allowing for population structure lead to estimates agreeing with the long-term rate for EBOV.

Molecular Phylogenies indicate a Paleo-Tibetan Origin of Himalayan Lazy Toads (Scutiger)

The Himalaya presents an outstanding geologically active orogen and biodiversity hotspot. However, our understanding of the historical biogeography of its fauna is far from comprehensive. Many taxa are commonly assumed to have originated from China-Indochina and dispersed westward along the Himalayan chain. Alternatively, the “Tibetan-origin hypothesis” suggests primary diversification of lineages in Paleo-Tibet, and secondary diversification along the slopes of the later uplifted Greater Himalaya. We test these hypotheses in high-mountain megophryid anurans (Scutiger). Extensive sampling from High Asia, and analyses of mitochondrial (2839 bp) and nuclear DNA (2208 bp), using Bayesian and Maximum likelihood phylogenetics, suggest that the Himalayan species form a distinct clade, possibly older than those from the eastern Himalaya-Tibet orogen. While immigration from China-Indochina cannot be excluded, our data may indicate that Himalayan Scutiger originated to the north of the Himalaya by colonization from Paleo-Tibet and then date back to the Oligocene. High intraspecific diversity of Scutiger implies limited migration across mountains and drainages along the Himalaya. While our study strengthens support for a “Tibetan-origin hypothesis”, current sampling (10/22 species; 1 revalidated: S. occidentalis) remains insufficient to draw final conclusions on Scutiger but urges comparative phylogeographers to test alternative, geologically supported hypotheses for a true future understanding of Himalayan biogeography.

On the young age of intraspecific herbaceous taxa

Dated phylogenies rarely include the divergence times of sister intraspecific taxa, and when they do little is said about this subject. We show that over 90% of the intraspecific plant taxa found in a literature search are estimated to be 5 million yr old or younger, with only 4% of taxa estimated to be over 10 million yr old or older. A Bayesian analysis of intraspecific taxon ages indicates that indeed these taxa are expected to be < 10 million yr old. This result for the young age of intraspecific taxa is consistent with the earlier observation that post-pollination reproductive barriers develop between 5 and 10 million yr after lineage splitting, thus leading to species formation. If lineages have not graduated to the species level of divergence by 10 million yr or so, they are likely to have gone extinct by that time as a result of narrow geographical distributions, narrow niche breadths, and relatively small numbers across populations.

Molecular clocks and the early evolution of metazoan nervous systems

The timing of early animal evolution remains poorly resolved, yet remains critical for understanding nervous system evolution. Methods for estimating divergence times from sequence data have improved considerably, providing a more refined understanding of key divergences. The best molecular estimates point to the origin of metazoans and bilaterians tens to hundreds of millions of years earlier than their first appearances in the fossil record. Both the molecular and fossil records are compatible, however, with the possibility of tiny, unskeletonized, low energy budget animals during the Proterozoic that had planktonic, benthic, or meiofaunal lifestyles. Such animals would likely have had relatively simple nervous systems equipped primarily to detect food, avoid inhospitable environments and locate mates. The appearance of the first macropredators during the Cambrian would have changed the selective landscape dramatically, likely driving the evolution of complex sense organs, sophisticated sensory processing systems, and diverse effector systems involved in capturing prey and avoiding predation.

Phylogenetic incongruence and the evolutionary origins of cardenolide-resistant forms of Na(+) ,K(+) -ATPase in Danaus butterflies

Many distantly related insect species are specialized feeders of cardenolide-containing host plants such as milkweed (Asclepias spp.). Previous studies have revealed frequent, parallel substitution of a functionally important amino acid substitution (N122H) in the alpha subunit of Na(+) ,K(+) -ATPase in a number of these species. This substitution facilitates the ability of these insects to feed on their toxic hosts and sequester cardenolides for their own use in defense. Among milkweed butterflies of the genus Danaus, the previously established phylogeny for this group suggests that N122H arose independently and fixed in two distinct lineages. We reevaluate this conclusion by examining Danaus phylogenetic relationships using >400 orthologous gene sequences assembled from transcriptome data. Our results indicate that the three Danaus species known to harbor the N122H substitution are more closely related than previously thought, consistent with a single, common origin for N122H. However, we also find evidence of both incomplete lineage sorting and post-speciation genetic exchange among these butterfly species, raising the possibility of collateral evolution of cardenolide-insensitivity in this species group.

Inflation of Molecular Clock Rates and Dates: Molecular Phylogenetics, Biogeography, and Diversification of a Global Cicada Radiation from Australasia (Hemiptera: Cicadidae: Cicadettini)

Dated phylogenetic trees are important for studying mechanisms of diversification, and molecular clocks are important tools for studies of organisms lacking good fossil records. However, studies have begun to identify problems in molecular clock dates caused by uncertainty of the modeled molecular substitution process. Here we explore Bayesian relaxed-clock molecular dating while studying the biogeography of ca. 200 species from the global cicada tribe Cicadettini. Because the available fossils are few and uninformative, we calibrate our trees in part with a cytochrome oxidase I (COI) clock prior encompassing a range of literature estimates for arthropods. We show that tribe-level analyses calibrated solely with the COI clock recover extremely old dates that conflict with published estimates for two well-studied New Zealand subclades within Cicadettini. Additional subclade analyses suggest that COI relaxed-clock rates and maximum-likelihood branch lengths become inflated relative to EF-1[Formula: see text] intron and exon rates and branch lengths as clade age increases. We present corrected estimates derived from: (i) an extrapolated EF-1[Formula: see text] exon clock derived from COI-calibrated analysis within the largest New Zealand subclade; (ii) post hoc scaling of the tribe-level chronogram using results from subclade analyses; and (iii) exploitation of a geological calibration point associated with New Caledonia. We caution that considerable uncertainty is generated due to dependence of substitution estimates on both the taxon sample and the choice of model, including gamma category number and the choice of empirical versus estimated base frequencies. Our results suggest that diversification of the tribe Cicadettini commenced in the early- to mid-Cenozoic and continued with the development of open, arid habitats in Australia and worldwide. We find that Cicadettini is a rare example of a global terrestrial animal group with an Australasian origin, with all non-Australasian genera belonging to two distal clades. Within Australia, we show that Cicadettini is more widely distributed than any other cicada tribe, diverse in temperate, arid and monsoonal habitats, and nearly absent from rainforests. We comment on the taxonomic implications of our findings for thirteen cicada genera.

Potential for bias and low precision in molecular divergence time estimation of the Canopy of Life: an example from aquatic bird families

Uncertainty in divergence time estimation is frequently studied from many angles but rarely from the perspective of phylogenetic node age. If appropriate molecular models and fossil priors are used, a multi-locus, partitioned analysis is expected to equally minimize error in accuracy and precision across all nodes of a given phylogeny. In contrast, if available models fail to completely account for rate heterogeneity, substitution saturation and incompleteness of the fossil record, uncertainty in divergence time estimation may increase with node age. While many studies have stressed this concern with regard to deep nodes in the Tree of Life, the inference that molecular divergence time estimation of shallow nodes is less sensitive to erroneous model choice has not been tested explicitly in a Bayesian framework. Because of available divergence time estimation methods that permit fossil priors across any phylogenetic node and the present increase in efficient, cheap collection of species-level genomic data, insight is needed into the performance of divergence time estimation of shallow (<10 MY) nodes. Here, we performed multiple sensitivity analyses in a multi-locus data set of aquatic birds with six fossil constraints. Comparison across divergence time analyses that varied taxon and locus sampling, number and position of fossil constraint and shape of prior distribution showed various insights. Deviation from node ages obtained from a reference analysis was generally highest for the shallowest nodes but determined more by temporal placement than number of fossil constraints. Calibration with only the shallowest nodes significantly underestimated the aquatic bird fossil record, indicating the presence of saturation. Although joint calibration with all six priors yielded ages most consistent with the fossil record, ages of shallow nodes were overestimated. This bias was found in both mtDNA and nDNA regions. Thus, divergence time estimation of shallow nodes may suffer from bias and low precision, even when appropriate fossil priors and best available substitution models are chosen. Much care must be taken to address the possible ramifications of substitution saturation across the entire Tree of Life.

Revisiting the phylogenetic relationships, biogeography, and taxonomy of spider monkeys (genus Ateles) in light of new molecular data

Spider monkeys (Ateles) are one of the most endangered groups of primates in the Neotropics. The genus is widely distributed from Mexico to the north of Bolivia and includes many morphologically distinct forms in terms of pelage color and patterning. The taxonomy, phylogenetic relationships, and biogeographic history of the genus have been subject to much debate, making scientific communication difficult and creating challenges for conservation actions. We extracted DNA from samples of all currently recognized species of spider monkeys collected from across the geographic range of the genus, sequenced ∼3.5 kilobases of coding sequence from the mitochondrial genome, and used this large dataset to (a) infer the phylogenetic relationships among the different forms of spider monkeys, (b) evaluate whether currently recognized species of spider monkeys form reciprocally monophyletic groups that are concordant with contemporary classifications, and (c) estimate divergence dates among the different lineages of Ateles. We found that all proposed species of spider monkeys for which we have samples from multiple localities indeed appear to form monophyletic groups. However, in contrast to previous studies, several of our analyses robustly inferred Ateles marginatus from northeast Brazil as the sister taxon to all other spider monkeys. A Bayesian dating analysis suggests that the most recent common ancestor of extant Ateles dates to ∼6.7 Ma, in the late Miocene, and most species-level splits within the genus took place in the late Pliocene, suggesting that the modern diversity in spider monkeys cannot be explained principally by isolation and divergence of populations in forest refugia during the Pleistocene. Based on our new phylogenetic inference and dating analysis, we propose a revised biogeographic scenario for the evolution of this genus.

The influence of ignoring secondary structure on divergence time estimates from ribosomal RNA genes

Genes coding for ribosomal RNA molecules (rDNA) are among the most popular markers in molecular phylogenetics and evolution. However, coevolution of sites that code for pairing regions (stems) in the RNA secondary structure can make it challenging to obtain accurate results from such loci. While the influence of ignoring secondary structure on multiple sequence alignment and tree topology has been investigated in numerous studies, its effect on molecular divergence time estimates is still poorly known. Here, I investigate this issue in Bayesian Markov Chain Monte Carlo (BMCMC) and penalized likelihood (PL) frameworks, using empirical datasets from dragonflies (Odonata: Anisoptera) and glass sponges (Porifera: Hexactinellida). My results indicate that highly biased inferences under substitution models that ignore secondary structure only occur if maximum-likelihood estimates of branch lengths are used as input to PL dating, whereas in a BMCMC framework and in PL dating based on Bayesian consensus branch lengths, the effect is far less severe. I conclude that accounting for coevolution of paired sites in molecular dating studies is not as important as previously suggested, as long as the estimates are based on Bayesian consensus branch lengths instead of ML point estimates. This finding is especially relevant for studies where computational limitations do not allow the use of secondary-structure specific substitution models, or where accurate consensus structures cannot be predicted. I also found that the magnitude and direction (over- vs. underestimating node ages) of bias in age estimates when secondary structure is ignored was not distributed randomly across the nodes of the phylogenies, a phenomenon that requires further investigation.

The timetable for allopolyploidy in flowering plants

BACKGROUND

Our understanding of the processes and dynamics of allopolyploid speciation, the long-term consequences of ploidal change, and the genetic and chromosomal changes in new emerged allopolyploids has substantially increased during the past few decades. Yet we remain uncertain about the time since lineage divergence when two taxa are capable of spawning such entities. Indeed, the matter has seemed intractable. Knowledge of the window of opportunity for allopolyploid production is very important because it provides temporal insight into a key evolutionary process, and a temporal reference against which other modes of speciation may be measured.

SCOPE

This Viewpoint paper reviews and integrates published information on the crossability of herbaceous species and the fertility of their hybrids in relation to species' divergence times. Despite limitations in methodology and sampling, the estimated times to hybrid sterility are somewhat congruent across disparate lineages. Whereas the waiting time for hybrid sterility is roughly 4-5 million years, the waiting time for cross-incompatibility is roughly 8-10 million years, sometimes considerably more. Strict allopolyploids may be formed in the intervening time window. The progenitors of several allopolyploids diverged between 4 and 6 million years before allopolyploid synthesis, as expected. This is the first study to propose a general temporal framework for strict allopolyploidy. This Viewpoint paper hopefully will stimulate interest in studying the tempo of speciation and the tempo of reproductive isolation in general.

Assessing parameter identifiability in phylogenetic models using data cloning

The success of model-based methods in phylogenetics has motivated much research aimed at generating new, biologically informative models. This new computer-intensive approach to phylogenetics demands validation studies and sound measures of performance. To date there has been little practical guidance available as to when and why the parameters in a particular model can be identified reliably. Here, we illustrate how Data Cloning (DC), a recently developed methodology to compute the maximum likelihood estimates along with their asymptotic variance, can be used to diagnose structural parameter nonidentifiability (NI) and distinguish it from other parameter estimability problems, including when parameters are structurally identifiable, but are not estimable in a given data set (INE), and when parameters are identifiable, and estimable, but only weakly so (WE). The application of the DC theorem uses well-known and widely used Bayesian computational techniques. With the DC approach, practitioners can use Bayesian phylogenetics software to diagnose nonidentifiability. Theoreticians and practitioners alike now have a powerful, yet simple tool to detect nonidentifiability while investigating complex modeling scenarios, where getting closed-form expressions in a probabilistic study is complicated. Furthermore, here we also show how DC can be used as a tool to examine and eliminate the influence of the priors, in particular if the process of prior elicitation is not straightforward. Finally, when applied to phylogenetic inference, DC can be used to study at least two important statistical questions: assessing identifiability of discrete parameters, like the tree topology, and developing efficient sampling methods for computationally expensive posterior densities.

Winding up the molecular clock in the genus Carabus (Coleoptera: Carabidae): assessment of methodological decisions on rate and node age estimation

BACKGROUND

Rates of molecular evolution are known to vary across taxa and among genes, and this requires rate calibration for each specific dataset based on external information. Calibration is sensitive to evolutionary model parameters, partitioning schemes and clock model. However, the way in which these and other analytical aspects affect both the rates and the resulting clade ages from calibrated phylogenies are not yet well understood. To investigate these aspects we have conducted calibration analyses for the genus Carabus (Coleoptera, Carabidae) on five mitochondrial and four nuclear DNA fragments with 7888 nt total length, testing different clock models and partitioning schemes to select the most suitable using Bayes Factors comparisons.

RESULTS

We used these data to investigate the effect of ambiguous character and outgroup inclusion on both the rates of molecular evolution and the TMRCA of Carabus. We found considerable variation in rates of molecular evolution depending on the fragment studied (ranging from 5.02% in cob to 0.26% divergence/My in LSU-A), but also on analytical conditions. Alternative choices of clock model, partitioning scheme, treatment of ambiguous characters, and outgroup inclusion resulted in rate increments ranging from 28% (HUWE1) to 1000% (LSU-B and ITS2) and increments in the TMRCA of Carabus ranging from 8.4% (cox1-A) to 540% (ITS2). Results support an origin of the genus Carabus during the Oligocene in the Eurasian continent followed by a Miocene differentiation that originated all main extant lineages.

CONCLUSIONS

The combination of several genes is proposed as the best strategy to minimise both the idiosyncratic behaviors of individual markers and the effect of analytical aspects in rate and age estimations. Our results highlight the importance of estimating rates of molecular evolution for each specific dataset, selecting for optimal clock and partitioning models as well as other methodological issues potentially affecting rate estimation.

High relatedness in a social amoeba: the role of kin-discriminatory segregation

A major challenge for social theory is to explain the importance of kin discrimination for the evolution of altruism. One way to assess the importance of kin discrimination is to test its effects on increasing relatedness within groups. The social amoeba Dictyostelium discoideum aggregates to form a fruiting body composed of dead stalk and live spores. Previous studies of a natural population showed that where D. discoideum occurs in the soil, multiple clones are often found in the same small soil samples. However, actual fruiting bodies usually contain only one clone. We here performed experiments to gauge the effect of kin-discriminatory segregation on increasing relatedness. We mixed co-occurring clones from this population using a relatedness level found in small soil samples. We found a lower proportion of uniclonal fruiting bodies and a lower level of relatedness compared with natural fruiting bodies. We found that the amount of relatedness increase attributable to kin-discriminatory segregation was small. These findings suggest a relatively minor influence of kin-discriminatory segregation on relatedness in D. discoideum. We discuss our results comparing with the results of previous studies, including those of wild clones and laboratory mutants. We ask why wild clones of D. discoideum exhibit a low degree of kin-discriminatory segregation, and what alternative factors might account for high relatedness in D. discoideum.

The shape of mammalian phylogeny: patterns, processes and scales

Mammalian phylogeny is far too asymmetric for all contemporaneous lineages to have had equal chances of diversifying. We consider this asymmetry or imbalance from four perspectives. First, we infer a minimal set of 'regime changes'-points at which net diversification rate has changed-identifying 15 significant radiations and 12 clades that may be 'downshifts'. We next show that mammalian phylogeny is similar in shape to a large set of published phylogenies of other vertebrate, arthropod and plant groups, suggesting that many clades may diversify under a largely shared set of 'rules'. Third, we simulate six simple macroevolutionary models, showing that those where speciation slows down as geographical or niche space is filled, produce more realistic phylogenies than do models involving key innovations. Lastly, an analysis of the spatial scaling of imbalance shows that the phylogeny of species within an assemblage, ecoregion or larger area always tends to be more unbalanced than expected from the phylogeny of species at the next more inclusive spatial scale. We conclude with a verbal model of mammalian macroevolution, which emphasizes the importance to diversification of accessing new regions of geographical or niche space.

Patterns of mammalian diversification in recent evolutionary times: global tendencies and methodological issues

Changes in diversification patterns estimated from phylogenetic trees are an important source of information about the dynamics of evolution. To study the diversification of mammals, we reconstructed phylogenetic trees of 29 families and fitted both constant-rate and variable-rate models of diversification. In addition, we investigated the effect of clock models and phylogenetic reconstruction problems on diversification analyses. We observed, first, that none of the families increased its diversification rate during the last few million years, including the Pleistocene. Furthermore, we detected a decrease in diversification that, after application of different tests, was significant only for a minority of families. However, when diversification variation was analysed in a combined tree of all families, a global decline in diversification became significant. Therefore, although distorted by some methodological artefacts, we found an underlying signal of gradually decreasing diversification that suggests that ecological factors may have shaped the recent diversification of mammals.

Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches

Evolutionary and biogeographic studies increasingly rely on calibrated molecular clocks to date key events. Although there has been significant recent progress in development of the techniques used for molecular dating, many issues remain. In particular, controversies abound over the appropriate use and placement of fossils for calibrating molecular clocks. Several methods have been proposed for evaluating candidate fossils; however, few studies have compared the results obtained by different approaches. Moreover, no previous study has incorporated the effects of nucleotide saturation from different data types in the evaluation of candidate fossils. In order to address these issues, we compared three approaches for evaluating fossil calibrations: the single-fossil cross-validation method of Near, Meylan, and Shaffer (2005. Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am. Nat. 165:137-146), the empirical fossil coverage method of Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am. Nat. 171:726-742), and the Bayesian multicalibration method of Sanders and Lee (2007. Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol. Lett. 3:275-279) and explicitly incorporate the effects of data type (nuclear vs. mitochondrial DNA) for identifying the most reliable or congruent fossil calibrations. We used advanced (Caenophidian) snakes as a case study; however, our results are applicable to any taxonomic group with multiple candidate fossils, provided appropriate taxon sampling and sufficient molecular sequence data are available. We found that data type strongly influenced which fossil calibrations were identified as outliers, regardless of which method was used. Despite the use of complex partitioned models of sequence evolution and multiple calibrations throughout the tree, saturation severely compressed basal branch lengths obtained from mitochondrial DNA compared with nuclear DNA. The effects of mitochondrial saturation were not ameliorated by analyzing a combined nuclear and mitochondrial data set. Although removing the third codon positions from the mitochondrial coding regions did not ameliorate saturation effects in the single-fossil cross-validations, it did in the Bayesian multicalibration analyses. Saturation significantly influenced the fossils that were selected as most reliable for all three methods evaluated. Our findings highlight the need to critically evaluate the fossils selected by data with different rates of nucleotide substitution and how data with different evolutionary rates affect the results of each method for evaluating fossils. Our empirical evaluation demonstrates that the advantages of using multiple independent fossil calibrations significantly outweigh any disadvantages.

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.

A multigenic perspective on phylogenetic relationships in the largest family of salamanders, the Plethodontidae

Despite several recent studies, the phylogeny of plethodontid salamanders is not yet fully resolved and the phylogenetic positions of several key genera, especially Aneides, Hemidactylium, Hydromantes and Karsenia, are contentious. Here we present a combined dataset of complete mitochondrial genomes and three nuclear loci for 20 species (16 genera) of plethodontids, representing all major clades in the family. The combined dataset without mitochondrial third codon positions provides a fully resolved, statistically well-supported tree. In this topology two major clades are recovered. A northern clade includes Aneides, Desmognathus, Ensatina, Hydromantes, Karsenia, Phaeognathus and Plethodon, with Plethodon being the sister taxon to the rest of the clade. Hydromantes and Karsenia are sister taxa, and Aneides is recovered as the sister taxon to Ensatina. Desmognathus+Phaeognathus form the sister taxon to Aneides+Ensatina. An eastern/southern clade comprises two subclades. One subclade, the spelerpines (Eurycea, Gyrinophilus, Pseudotriton, Stereochilus, Urspelerpes) is the sister taxon to a subclade comprising Hemidactylium, Batrachoseps and the tropical plethodontids (represented by Bolitoglossa, Nototriton and Thorius). In this topology Hemidactylium is well-supported as the sister taxon to Batrachoseps. Only when mitochondrial third codon positions are included using maximum likelihood analysis is Hemidactylium recovered as the sister taxon to Batrachoseps+tropical genera. Hypothesis testing of alternative topologies supports these conclusions. On the basis of these results we propose a conservative taxonomy for Plethodontidae.

Cleaning up the 'Bigmessidae': molecular phylogeny of scleractinian corals from Faviidae, Merulinidae, Pectiniidae and Trachyphylliidae

BACKGROUND

Molecular phylogenetic studies on scleractinian corals have shown that most taxa are not reflective of their evolutionary histories. Based principally on gross morphology, traditional taxonomy suffers from the lack of well-defined and homologous characters that can sufficiently describe scleractinian diversity. One of the most challenging clades recovered by recent analyses is 'Bigmessidae', an informal grouping that comprises four conventional coral families, Faviidae, Merulinidae, Pectiniidae and Trachyphylliidae, interspersed among one another with no apparent systematic pattern. There is an urgent need for taxonomic revisions in this clade, but it is vital to first establish phylogenetic relationships within the group. In this study, we reconstruct the evolutionary history of 'Bigmessidae' based on five DNA sequence markers gathered from 76 of the 132 currently recognized species collected from five reef regions in the central Indo-Pacific and the Atlantic.

RESULTS

We present a robust molecular phylogeny of 'Bigmessidae' based on the combined five-gene data, achieving a higher degree of resolution compared to previous analyses. Two Pacific species presumed to be in 'Bigmessidae' are more closely related to outgroup clades, suggesting that other unsampled taxa have unforeseen affinities. As expected, nested within 'Bigmessidae' are four conventional families as listed above, and relationships among them generally corroborate previous molecular analyses. Our more resolved phylogeny supports a close association of Hydnophora (Merulinidae) with Favites + Montastraea (Faviidae), rather than with the rest of Merulinidae, i.e., Merulina and Scapophyllia. Montastraea annularis, the only Atlantic 'Bigmessidae' is sister to Cyphastrea, a grouping that can be reconciled by their septothecal walls, a microstructural feature of the skeleton determined by recent morphological work. Characters at the subcorallite scale appear to be appropriate synapomorphies for other subclades, which cannot be explained using macromorphology. Indeed, wide geographic sampling here has revealed more instances of possible cryptic taxa confused by evolutionary convergence of gross coral morphology.

CONCLUSIONS

Numerous examples of cryptic taxa determined in this study support the assertion that diversity estimates of scleractinian corals are erroneous. Fortunately, the recovery of most 'Bigmessidae' genera with only minor degrees of paraphyly offers some hope for impending taxonomic amendments. Subclades are well defined and supported by subcorallite morphological features, providing a robust framework for further systematic work.

The evolutionary history of Beet necrotic yellow vein virus deduced from genetic variation, geographical origin and spread, and the breaking of host resistance

Beet necrotic yellow vein virus (BNYVV) is an economically important pathogen of sugar beet and has been found worldwide, probably as the result of recent worldwide spread. The BNYVV genome consists of four or five RNA components. Here, we report analysis of sequence variation in the RNA3-p25, RNA4-p31, RNA2-CP, and RNA5-p26 genes of 73 worldwide isolates. The RNA3-p25 gene encodes virulence and avirulence factors. These four sets of gene sequences each fell into two to four groups, of which the three groups of p25 formed eight subgroups with different geographical distributions. Each of these subgroup isolates (strains) could have arisen from four original BNYVV population and their mixed infections. The genetic diversity for BNYVV was relatively small. Selection pressure varied greatly depending on the BNYVV gene and geographical location. Isolates of the Italy strain, in which p25 was subject to the strongest positive selection, were able to overcome the Rz1-host resistance gene to differing degrees, whereas other geographically limited strains could not. Resistance-breaking variants were generated by p25 amino acid changes at positions 67 and 68. Our studies suggest that BNYVV originally evolved in East Asia and has recently become a pathogen of cultivated sugar beet followed by the emergence of new resistance-breaking variants.

Circadian input kinases and their homologs in cyanobacteria: evolutionary constraints versus architectural diversification

The circadian input kinase A (cikA) gene encodes a protein relaying environmental signal to the central circadian oscillator in cyanobacteria. The CikA protein has a variable architecture and usually consists of four tandemly arrayed domains: GAF, histidine kinase (HisKA), histidine kinase-like ATPase (HATPase_c), and a pseudo-receiver (REC). Among them, HisKA and HATPase_c are the least polymorphic, and REC is not present in heterocystic filamentous cyanobacteria. CikA contains several conserved motifs that are likely important for circadian function. There are at least three types of circadian systems, each of which possesses a different set of circadian genes. The originally described circadian system (kaiABC system) possesses both cikA and kaiA, while the others lack either only cikA (kaiABC (Delta)) or both (kaiBC). The results we obtained allowed us to approximate the time of the cikA origin to be about 2600-2200 MYA and the time of its loss in the species with the kaiABC (Delta) or kaiBC system between 1100 and 600 MYA. Circadian specialization of CikA, as opposed to its non-circadian homologs, is a result of several factors, including the unique conserved domain architecture and high evolutionary constraints of some domains and regions, which were previously identified as critical for the circadian function of the gene.

Variation in DNA substitution rates among lineages erroneously inferred from simulated clock-like data

BACKGROUND

The observation of variation in substitution rates among lineages has led to (1) a general rejection of the molecular clock model, and (2) the suggestion that a number of biological characteristics of organisms can cause rate variation. Accurate estimates of rate variation, and thus accurate inferences regarding the causes of rate variation, depend on accurate estimates of substitution rates. However, theory suggests that even when the substitution process is clock-like, variable numbers of substitutions can occur among lineages because the substitution process is stochastic. Furthermore, substitution rates along lineages can be misestimated, particularly when multiple substitutions occur at some sites. Although these potential causes of error in rate estimation are well understood in theory, such error has not been examined in detail; consequently, empirical studies that estimate rate variation among lineages have been unable to determine whether their results could be impacted by estimation error.

METHODOLOGY/PRINCIPAL FINDINGS

To evaluate the extent to which error in rate estimation could erroneously suggest rate variation among lineages, we examined rate variation estimated for datasets simulated under a molecular clock on trees with equal and variable branch lengths. Thus, any apparent rate variation in these datasets reflects error in rate estimation rather than true differences in the underlying substitution process. We observed substantial rate variation among lineages in our simulations; however, we did not observe rate variation when average substitution rates were compared between different clades.

CONCLUSIONS/SIGNIFICANCE

Our results confirm previous theoretical work suggesting that observations of among lineage rate variation in empirical data may be due to the stochastic substitution process and error in the estimation of substitution rates, rather than true differences in the underlying substitution process among lineages. However, conclusions regarding rate variation drawn from rates averaged across multiple branches are likely due to real, systematic variation in rates between groups.