The power of relative rates tests depends on the data

Diversification dynamics of a common deep-sea octocoral family linked to the Paleocene-Eocene thermal maximum

The deep-sea has experienced dramatic changes in physical and chemical variables in the geological past. However, little is known about how deep-sea species richness responded to such changes over time and space. Here, we studied the diversification dynamics of one of the most diverse octocorallian families inhabiting deep sea benthonic environments worldwide and sustaining highly diverse ecosystems, Primnoidae. A newly dated species-level phylogeny was constructed to infer their ancestral geographic locations and dispersal rates initially. Then, we tested whether their global and regional (the Southern Ocean) diversification dynamics were mediated by dispersal rate and abiotic factors as changes in ocean geochemistry. Finally, we tested whether primnoids showed changes in speciation and extinction at discrete time points. Our results suggested primnoids likely originated in the southwestern Pacific Ocean during the Lower Cretaceous ∼112 Ma, with further dispersal after the physical separation of continental landmasses along the late Mesozoic and Cenozoic. Only the speciation rate of the Southern Ocean primnoids showed a significant correlation to ocean chemistry. Moreover, the Paleocene-Eocene thermal maximum marked a significant increase in the diversification of primnoids at global and regional scales. Our results provide new perspectives on the macroevolutionary and biogeographic patterns of an ecologically important benthic organism typically found in deep-sea environments.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

A Not-So-Long Introduction to Computational Molecular Evolution

In this chapter, we give a not-so-long and self-contained introduction to computational molecular evolution. In particular, we present the emergence of the use of likelihood-based methods, review the standard DNA substitution models, and introduce how model choice operates. We also present recent developments in inferring absolute divergence times and rates on a phylogeny, before showing how state-of-the-art models take inspiration from diffusion theory to link population genetics, which traditionally focuses at a taxonomic level below that of the species, and molecular evolution. Although this is not a cookbook chapter, we try and point to popular programs and implementations along the way.

Accounting for Calibration Uncertainty: Bayesian Molecular Dating as a "Doubly Intractable" Problem

This study introduces a new Bayesian technique for molecular dating that explicitly accommodates for uncertainty in the phylogenetic position of calibrated nodes derived from the analysis of fossil data. The proposed approach thus defines an adequate framework for incorporating expert knowledge and/or prior information about the way fossils were collected in the inference of node ages. Although it belongs to the class of "node-dating" approaches, this method shares interesting properties with "tip-dating" techniques. Yet, it alleviates some of the computational and modeling difficulties that hamper tip-dating approaches. The influence of fossil data on the probabilistic distribution of trees is the crux of the matter considered here. More specifically, among all the phylogenies that a tree model (e.g., the birth-death process) generates, only a fraction of them "agree" with the fossil data. Bayesian inference under the new model requires taking this fraction into account. However, evaluating this quantity is difficult in practice. A generic solution to this issue is presented here. The proposed approach relies on a recent statistical technique, the so-called exchange algorithm, dedicated to drawing samples from "doubly intractable" distributions. A small example illustrates the problem of interest and the impact of uncertainty in the placement of calibration constraints in the phylogeny given fossil data. An analysis of land plant sequences and multiple fossils further highlights the pertinence of the proposed approach.

Population Size and the Rate of Language Evolution: A Test Across Indo-European, Austronesian, and Bantu Languages

What role does speaker population size play in shaping rates of language evolution? There has been little consensus on the expected relationship between rates and patterns of language change and speaker population size, with some predicting faster rates of change in smaller populations, and others expecting greater change in larger populations. The growth of comparative databases has allowed population size effects to be investigated across a wide range of language groups, with mixed results. One recent study of a group of Polynesian languages revealed greater rates of word gain in larger populations and greater rates of word loss in smaller populations. However, that test was restricted to 20 closely related languages from small Oceanic islands. Here, we test if this pattern is a general feature of language evolution across a larger and more diverse sample of languages from both continental and island populations. We analyzed comparative language data for 153 pairs of closely-related sister languages from three of the world's largest language families: Austronesian, Indo-European, and Niger-Congo. We find some evidence that rates of word loss are significantly greater in smaller languages for the Indo-European comparisons, but we find no significant patterns in the other two language families. These results suggest either that the influence of population size on rates and patterns of language evolution is not universal, or that it is sufficiently weak that it may be overwhelmed by other influences in some cases. Further investigation, for a greater number of language comparisons and a wider range of language features, may determine which of these explanations holds true.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

The Complete Chloroplast Genome of Guadua angustifolia and Comparative Analyses of Neotropical-Paleotropical Bamboos

To elucidate chloroplast genome evolution within neotropical-paleotropical bamboos, we fully characterized the chloroplast genome of the woody bamboo Guadua angustifolia. This genome is 135,331 bp long and comprises of an 82,839-bp large single-copy (LSC) region, a 12,898-bp small single-copy (SSC) region, and a pair of 19,797-bp inverted repeats (IRs). Comparative analyses revealed marked conservation of gene content and sequence evolutionary rates between neotropical and paleotropical woody bamboos. The neotropical herbaceous bamboo Cryptochloa strictiflora differs from woody bamboos in IR/SSC boundaries in that it exhibits slightly contracted IRs and a faster substitution rate. The G. angustifolia chloroplast genome is similar in size to that of neotropical herbaceous bamboos but is ~3 kb smaller than that of paleotropical woody bamboos. Dissimilarities in genome size are correlated with differences in the lengths of intergenic spacers, which are caused by large-fragment insertion and deletion. Phylogenomic analyses of 62 taxa yielded a tree topology identical to that found in preceding studies. Divergence time estimation suggested that most bamboo genera diverged after the Miocene and that speciation events of extant species occurred during or after the Pliocene.

Plasmodium falciparum field isolates from areas of repeated emergence of drug resistant malaria show no evidence of hypermutator phenotype

Multiple transcontinental waves of drug resistance in Plasmodium falciparum have originated in Southeast Asia before spreading westward, first into the rest of Asia and then to sub-Saharan Africa. In vitro studies have suggested that hypermutator P. falciparum parasites may exist in Southeast Asia and that an increased rate of acquisition of new mutations in these parasites may explain the repeated emergence of drug resistance in Southeast Asia. This study is the first to test the hypermutator hypothesis using field isolates. Using genome-wide SNP data from human P. falciparum infections in Southeast Asia and West Africa and a test for relative rate differences we found no evidence of increased relative substitution rates in P. falciparum isolates from Southeast Asia. Instead, we found significantly increased substitution rates in Mali and Bangladesh populations relative to those in populations from Southeast Asia. Additionally we found no association between increased relative substitution rates and parasite clearance following treatment with artemisinin derivatives.

Flying rocks and flying clocks: disparity in fossil and molecular dates for birds

Major disparities are recognized between molecular divergence dates and fossil ages for critical nodes in the Tree of Life, but broad patterns and underlying drivers remain elusive. We harvested 458 molecular age estimates for the stem and crown divergences of 67 avian clades to explore empirical patterns between these alternate sources of temporal information. These divergence estimates were, on average, over twice the age of the oldest fossil in these clades. Mitochondrial studies yielded older ages than nuclear studies for the vast majority of clades. Unexpectedly, disparity between molecular estimates and the fossil record was higher for divergences within major clades (crown divergences) than divergences between major clades (stem divergences). Comparisons of dates from studies classed by analytical methods revealed few significant differences. Because true divergence ages can never be known with certainty, our study does not answer the question of whether fossil gaps or molecular dating error account for a greater proportion of observed disparity. However, empirical patterns observed here suggest systemic overestimates for shallow nodes in existing molecular divergence dates for birds. We discuss underlying biases that may drive these patterns.

Rate of language evolution is affected by population size

The effect of population size on patterns and rates of language evolution is controversial. Do languages with larger speaker populations change faster due to a greater capacity for innovation, or do smaller populations change faster due to more efficient diffusion of innovations? Do smaller populations suffer greater loss of language elements through founder effects or drift, or do languages with more speakers lose features due to a process of simplification? Revealing the influence of population size on the tempo and mode of language evolution not only will clarify underlying mechanisms of language change but also has practical implications for the way that language data are used to reconstruct the history of human cultures. Here, we provide, to our knowledge, the first empirical, statistically robust test of the influence of population size on rates of language evolution, controlling for the evolutionary history of the populations and formally comparing the fit of different models of language evolution. We compare rates of gain and loss of cognate words for basic vocabulary in Polynesian languages, an ideal test case with a well-defined history. We demonstrate that larger populations have higher rates of gain of new words whereas smaller populations have higher rates of word loss. These results show that demographic factors can influence rates of language evolution and that rates of gain and loss are affected differently. These findings are strikingly consistent with general predictions of evolutionary models.

Signals of historical interlocus gene conversion in human segmental duplications

Standard methods of DNA sequence analysis assume that sequences evolve independently, yet this assumption may not be appropriate for segmental duplications that exchange variants via interlocus gene conversion (IGC). Here, we use high quality multiple sequence alignments from well-annotated segmental duplications to systematically identify IGC signals in the human reference genome. Our analysis combines two complementary methods: (i) a paralog quartet method that uses DNA sequence simulations to identify a statistical excess of sites consistent with inter-paralog exchange, and (ii) the alignment-based method implemented in the GENECONV program. One-quarter (25.4%) of the paralog families in our analysis harbor clear IGC signals by the quartet approach. Using GENECONV, we identify 1477 gene conversion tracks that cumulatively span 1.54 Mb of the genome. Our analyses confirm the previously reported high rates of IGC in subtelomeric regions and Y-chromosome palindromes, and identify multiple novel IGC hotspots, including the pregnancy specific glycoproteins and the neuroblastoma breakpoint gene families. Although the duplication history of a paralog family is described by a single tree, we show that IGC has introduced incredible site-to-site variation in the evolutionary relationships among paralogs in the human genome. Our findings indicate that IGC has left significant footprints in patterns of sequence diversity across segmental duplications in the human genome, out-pacing the contributions of single base mutation by orders of magnitude. Collectively, the IGC signals we report comprise a catalog that will provide a critical reference for interpreting observed patterns of DNA sequence variation across duplicated genomic regions, including targets of recent adaptive evolution in humans.

Low diversity in the mitogenome of sperm whales revealed by next-generation sequencing

Large population sizes and global distributions generally associate with high mitochondrial DNA control region (CR) diversity. The sperm whale (Physeter macrocephalus) is an exception, showing low CR diversity relative to other cetaceans; however, diversity levels throughout the remainder of the sperm whale mitogenome are unknown. We sequenced 20 mitogenomes from 17 sperm whales representative of worldwide diversity using Next Generation Sequencing (NGS) technologies (Illumina GAIIx, Roche 454 GS Junior). Resequencing of three individuals with both NGS platforms and partial Sanger sequencing showed low discrepancy rates (454-Illumina: 0.0071%; Sanger-Illumina: 0.0034%; and Sanger-454: 0.0023%) confirming suitability of both NGS platforms for investigating low mitogenomic diversity. Using the 17 sperm whale mitogenomes in a phylogenetic reconstruction with 41 other species, including 11 new dolphin mitogenomes, we tested two hypotheses for the low CR diversity. First, the hypothesis that CR-specific constraints have reduced diversity solely in the CR was rejected as diversity was low throughout the mitogenome, not just in the CR (overall diversity π = 0.096%; protein-coding 3rd codon = 0.22%; CR = 0.35%), and CR phylogenetic signal was congruent with protein-coding regions. Second, the hypothesis that slow substitution rates reduced diversity throughout the sperm whale mitogenome was rejected as sperm whales had significantly higher rates of CR evolution and no evidence of slow coding region evolution relative to other cetaceans. The estimated time to most recent common ancestor for sperm whale mitogenomes was 72,800 to 137,400 years ago (95% highest probability density interval), consistent with previous hypotheses of a bottleneck or selective sweep as likely causes of low mitogenome diversity.

Epitope of titin A-band-specific monoclonal antibody Tit1 5 H1.1 is highly conserved in several Fn3 domains of the titin molecule. Centriole staining in human, mouse and zebrafish cells

Background

Previously we have reported on the development of a new mouse anti-titin monoclonal antibody, named MAb Titl 5 H1.1, using the synthetic peptide N-AVNKYGIGEPLESDSVVAK-C which corresponds to an amino acid sequence in the A-region of the titin molecule as immunogen. In the human skeletal muscles, MAb Titl 5 H1.1 reacts specifically with titin in the A-band of the sarcomere and in different non-muscle cell types with nucleus and cytoplasm, including centrioles. In this report we have studied the evolutionary aspects of the binding of MAb Tit1 5 H1.1 with its target antigen (titin).

Results

We have specified the epitope area of MAb Tit1 5 H1.1 by subpeptide mapping to the hexapeptide N-AVNKYG-C. According to protein databases this amino acid sequence is located in the COOH-terminus of several different Fn3 domains of the A-region of titin molecule in many organisms, such as human being, mouse, rabbit, zebrafish (Danio rerio), and even in sea squirt (Ciona intestinalis). Our immunohisto- and cytochemical studies with MAb Tit1 5 H1.1 in human, mouse and zebrafish tissues and cell cultures showed a striated staining pattern in muscle cells and also staining of centrioles, cytoplasm and nuclei in non-muscle cells.

Conclusions

The data confirm that titin can play, in addition to the known roles in striated muscle cells also an important role in non-muscle cells as a centriole associated protein. This phenomenon is highly conserved in the evolution and is related to Fn3 domains of the titin molecule. Using titin A-band-specific monoclonal antibody MAb Tit1 5 H1.1 it was possible to locate titin in the sarcomeres of skeletal muscle cells and in the centrioles, cytoplasm and nuclei of non-muscle cells in phylogenetically so distant organisms as Homo sapiens, Mus musculus and zebrafish (Danio rerio).

Cryptic species in putative ancient asexual darwinulids (Crustacea, Ostracoda)

BACKGROUND

Fully asexually reproducing taxa lack outcrossing. Hence, the classic Biological Species Concept cannot be applied.

METHODOLOGY/PRINCIPAL FINDINGS

We used DNA sequences from the mitochondrial COI gene and the nuclear ITS2 region to check species boundaries according to the evolutionary genetic (EG) species concept in five morphospecies in the putative ancient asexual ostracod genera, Penthesilenula and Darwinula, from different continents. We applied two methods for detecting cryptic species, namely the K/θ method and the General Mixed Yule Coalescent model (GMYC). We could confirm the existence of species in all five darwinulid morphospecies and additional cryptic diversity in three morphospecies, namely in Penthesilenula brasiliensis, Darwinula stevensoni and in P. aotearoa. The number of cryptic species within one morphospecies varied between seven (P. brasiliensis), five to six (D. stevensoni) and two (P. aotearoa), respectively, depending on the method used. Cryptic species mainly followed continental distributions. We also found evidence for coexistence at the local scale for Brazilian cryptic species of P. brasiliensis and P. aotearoa. Our ITS2 data confirmed that species exist in darwinulids but detected far less EG species, namely two to three cryptic species in P. brasiliensis and no cryptic species at all in the other darwinulid morphospecies.

CONCLUSIONS/SIGNIFICANCE

Our results clearly demonstrate that both species and cryptic diversity can be recognized in putative ancient asexual ostracods using the EG species concept, and that COI data are more suitable than ITS2 for this purpose. The discovery of up to eight cryptic species within a single morphospecies will significantly increase estimates of biodiversity in this asexual ostracod group. Which factors, other than long-term geographic isolation, are important for speciation processes in these ancient asexuals remains to be investigated.

Fast molecular evolution associated with high active metabolic rates in poison frogs

Molecular evolution is simultaneously paced by mutation rate, genetic drift, and natural selection. Life history traits also affect the speed of accumulation of nucleotide changes. For instance, small body size, rapid generation time, production of reactive oxygen species (ROS), and high resting metabolic rate (RMR) are suggested to be associated with faster rates of molecular evolution. However, phylogenetic correlation analyses failed to support a relationship between RMR and molecular evolution in ectotherms. In addition, RMR might underestimate the metabolic budget (e.g., digestion, reproduction, or escaping predation). An alternative is to test other metabolic rates, such as active metabolic rate (AMR), and their association with molecular evolution. Here, I present comparative analyses of the associations between life history traits (i.e., AMR, RMR, body mass, and fecundity) with rates of molecular evolution of and mitochondrial loci from a large ectotherm clade, the poison frogs (Dendrobatidae). My results support a strong positive association between mass-specific AMR and rates of molecular evolution for both mitochondrial and nuclear loci. In addition, I found weaker and genome-specific covariates such as body mass and fecundity for mitochondrial and nuclear loci, respectively. No direct association was found between mass-specific RMR and rates of molecular evolution. Thus, I provide a mechanistic hypothesis of the link between AMRs and the rate of molecular evolution based on an increase in ROS within germ line cells during periodic bouts of hypoxia/hyperoxia related to aerobic exercise. Finally, I propose a multifactorial model that includes AMR as a predictor of the rate of molecular evolution in ectothermic lineages.

The essentials of computational molecular evolution

In this chapter, we give a brief yet self-contained introduction to computational molecular evolution. In particular, we present the emergence of the use of likelihood-based methods, review the standard DNA substitution models, and introduce how model choice operates. We also present recent developments in inferring absolute dates and rates on a phylogeny and show how state-of-the-art models take inspiration from diffusion theory to link population genetics, which traditionally focuses at a taxonomic level under that of species, and molecular evolution.

Gone with the plate: the opening of the Western Mediterranean basin drove the diversification of ground-dweller spiders

BACKGROUND

The major islands of the Western Mediterranean--Corsica, Sardinia, and the Balearic Islands--are continental terrenes that drifted towards their present day location following a retreat from their original position on the eastern Iberian Peninsula about 30 million years ago. Several studies have taken advantage of this well-dated geological scenario to calibrate molecular rates in species for which distributions seemed to match this tectonic event. Nevertheless, the use of external calibration points has revealed that most of the present-day fauna on these islands post-dated the opening of the western Mediterranean basin. In this study, we use sequence information of the cox1, nad1, 16S, L1, and 12S mitochondrial genes and the 18S, 28S, and h3 nuclear genes, along with relaxed clock models and a combination of biogeographic and fossil external calibration points, to test alternative historical scenarios of the evolutionary history of the ground-dweller spider genus Parachtes (Dysderidae), which is endemic to the region.

RESULTS

We analyse 49 specimens representing populations of most Parachtes species and close relatives. Our results reveal that both the sequence of species formation in Parachtes and the estimated divergence times match the geochronological sequence of separation of the main islands, suggesting that the diversification of the group was driven by Tertiary plate tectonics. In addition, the confirmation that Parachtes diversification matches well-dated geological events provides a model framework to infer substitution rates of molecular markers. Divergence rates estimates ranged from 3.5% My(-1) (nad1) to 0.12% My(-1) (28S), and the average divergence rate for the mitochondrial genes was 2.25% My(-1), very close to the "standard" arthropod mitochondrial rate (2.3% My(-1)).

CONCLUSIONS

Our study provides the first unequivocal evidence of terrestrial endemic fauna of the major western Mediterranean islands, whose origin can be traced back to the Oligocene separation of these islands from the continent. Moreover, our study provides useful information on the divergence rate estimates of the most commonly used genes for phylogenetic inference in non-model arthropods.

The uneven phylogeny and biogeography of Erodium (Geraniaceae): radiations in the Mediterranean and recent recurrent intercontinental colonization

BACKGROUND AND AIMS

The genus Erodium is a common feature of Mediterranean-type climates throughout the world, but the Mediterranean Basin has significantly higher diversity than other areas. The aim here is to reveal the biogeographical history of the genus and the causes behind the evolution of the uneven distribution.

METHODS

Seventy-eight new nrITS sequences were incorporated with existing plastid data to explore the phylogenetic relationships and biogeography of Erodium using several reconstruction methods. Divergence times for major clades were calculated and contrasted with other previously published information. Furthermore, topological and temporal diversification rate shift analyses were employed using these data.

KEY RESULTS

Phylogenetic relationships among species are widely congruent with previous plastid reconstructions, which refute the classical taxonomical classification. Biogeographical reconstructions point to Asia as the ancestral area of Erodium, arising approx. 18 MYA. Four incidences of intercontinental dispersal from the Mediterranean Basin to similar climates are demonstrated. Increases in diversification were present in two independent Erodium lineages concurrently. Two bursts of diversification (3 MYA and 0·69 MYA) were detected only in the Mediterranean flora.

CONCLUSIONS

Two lineages diverged early in the evolution of the genus Erodium: (1) subgenus Erodium plus subgenus Barbata subsection Absinthioidea and (2) the remainder of subgenus Barbata. Dispersal across major water bodies, although uncommon, has had a major influence on the distribution of this genus and is likely to have played as significant role as in other, more easily dispersed, genera. Establishment of Mediterranean climates has facilitated the spread of the genus and been crucial in its diversification. Two, independent, rapid radiations in response to the onset of drought and glacial climate change indicate putative adaptive radiations in the genus.

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.

Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria

Background

Because bacteria do not have a robust fossil record, attempts to infer the timing of events in their evolutionary history requires comparisons of molecular sequences. This use of molecular clocks is based on the assumptions that substitution rates for homologous genes or sites are fairly constant through time and across taxa. Violation of these conditions can lead to erroneous inferences and result in estimates that are off by orders of magnitude. In this study, we examine the consistency of substitution rates among a set of conserved genes in diverse bacterial lineages, and address the questions regarding the validity of molecular dating.

Results

By examining the evolution of 16S rRNA gene in obligate endosymbionts, which can be calibrated by the fossil record of their hosts, we found that the rates are consistent within a clade but varied widely across different bacterial lineages. Genome-wide estimates of nonsynonymous and synonymous substitutions suggest that these two measures are highly variable in their rates across bacterial taxa. Genetic drift plays a fundamental role in determining the accumulation of substitutions in 16S rRNA genes and at nonsynonymous sites. Moreover, divergence estimates based on a set of universally conserved protein-coding genes also exhibit low correspondence to those based on 16S rRNA genes.

Conclusion

Our results document a wide range of substitution rates across genes and bacterial taxa. This high level of variation cautions against the assumption of a universal molecular clock for inferring divergence times in bacteria. However, by applying relative-rate tests to homologous genes, it is possible to derive reliable local clocks that can be used to calibrate bacterial evolution.

Reviewers

This article was reviewed by Adam Eyre-Walker, Simonetta Gribaldo and Tal Pupko (nominated by Dan Graur).

Big and slow: phylogenetic estimates of molecular evolution in baleen whales (suborder mysticeti)

Baleen whales are the largest animals that have ever lived. To develop an improved estimation of substitution rate for nuclear and mitochondrial DNA for this taxon, we implemented a relaxed-clock phylogenetic approach using three fossil calibration dates: the divergence between odontocetes and mysticetes approximately 34 million years ago (Ma), between the balaenids and balaenopterids approximately 28 Ma, and the time to most recent common ancestor within the Balaenopteridae approximately 12 Ma. We examined seven mitochondrial genomes, a large number of mitochondrial control region sequences (219 haplotypes for 465 bp) and nine nuclear introns representing five species of whales, within which multiple species-specific alleles were sequenced to account for within-species diversity (1-15 for each locus). The total data set represents >1.65 Mbp of mitogenome and nuclear genomic sequence. The estimated substitution rate for the humpback whale control region (3.9%/million years, My) was higher than previous estimates for baleen whales but slow relative to other mammal species with similar generation times (e.g., human-chimp mean rate > 20%/My). The mitogenomic third codon position rate was also slow relative to other mammals (mean estimate 1%/My compared with a mammalian average of 9.8%/My for the cytochrome b gene). The mean nuclear genomic substitution rate (0.05%/My) was substantially slower than average synonymous estimates for other mammals (0.21-0.37%/My across a range of studies). The nuclear and mitogenome rate estimates for baleen whales were thus roughly consistent with an 8- to 10-fold slowing due to a combination of large body size and long generation times. Surprisingly, despite the large data set of nuclear intron sequences, there was only weak and conflicting support for alternate hypotheses about the phylogeny of balaenopterid whales, suggesting that interspecies introgressions or a rapid radiation has obscured species relationships in the nuclear genome.

Reconsidering the generation time hypothesis based on nuclear ribosomal ITS sequence comparisons in annual and perennial angiosperms

Background

Differences in plant annual/perennial habit are hypothesized to cause a generation time effect on divergence rates. Previous studies that compared rates of divergence for internal transcribed spacer (ITS1 and ITS2) sequences of nuclear ribosomal DNA (nrDNA) in angiosperms have reached contradictory conclusions about whether differences in generation times (or other life history features) are associated with divergence rate heterogeneity. We compared annual/perennial ITS divergence rates using published sequence data, employing sampling criteria to control for possible artifacts that might obscure any actual rate variation caused by annual/perennial differences.

Results

Relative rate tests employing ITS sequences from 16 phylogenetically-independent annual/perennial species pairs rejected rate homogeneity in only a few comparisons, with annuals more frequently exhibiting faster substitution rates. Treating branch length differences categorically (annual faster or perennial faster regardless of magnitude) with a sign test often indicated an excess of annuals with faster substitution rates. Annuals showed an approximately 1.6-fold rate acceleration in nucleotide substitution models for ITS. Relative rates of three nuclear loci and two chloroplast regions for the annual Arabidopsis thaliana compared with two closely related Arabidopsis perennials indicated that divergence was faster for the annual. In contrast, A. thaliana ITS divergence rates were sometimes faster and sometimes slower than the perennial. In simulations, divergence rate differences of at least 3.5-fold were required to reject rate constancy in > 80 % of replicates using a nucleotide substitution model observed for the combination of ITS1 and ITS2. Simulations also showed that categorical treatment of branch length differences detected rate heterogeneity > 80% of the time with a 1.5-fold or greater rate difference.

Conclusion

Although rate homogeneity was not rejected in many comparisons, in cases of significant rate heterogeneity annuals frequently exhibited faster substitution rates. Our results suggest that annual taxa may exhibit a less than 2-fold rate acceleration at ITS. Since the rate difference is small and ITS lacks statistical power to reject rate homogeneity, further studies with greater power will be required to adequately test the hypothesis that annual and perennial plants have heterogeneous substitution rates. Arabidopsis sequence data suggest that relative rate tests based on multiple loci may be able to distinguish a weak acceleration in annual plants. The failure to detect rate heterogeneity with ITS in past studies may be largely a product of low statistical power.

Parallel rate heterogeneity in chloroplast and mitochondrial genomes of Brazil nut trees (Lecythidaceae) is consistent with lineage effects

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication-dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.

Multigene phylogeny of the Mustelidae: resolving relationships, tempo and biogeographic history of a mammalian adaptive radiation

BACKGROUND

Adaptive radiation, the evolution of ecological and phenotypic diversity from a common ancestor, is a central concept in evolutionary biology and characterizes the evolutionary histories of many groups of organisms. One such group is the Mustelidae, the most species-rich family within the mammalian order Carnivora, encompassing 59 species classified into 22 genera. Extant mustelids display extensive ecomorphological diversity, with different lineages having evolved into an array of adaptive zones, from fossorial badgers to semi-aquatic otters. Mustelids are also widely distributed, with multiple genera found on different continents. As with other groups that have undergone adaptive radiation, resolving the phylogenetic history of mustelids presents a number of challenges because ecomorphological convergence may potentially confound morphologically based phylogenetic inferences, and because adaptive radiations often include one or more periods of rapid cladogenesis that require a large amount of data to resolve.

RESULTS

We constructed a nearly complete generic-level phylogeny of the Mustelidae using a data matrix comprising 22 gene segments (approximately 12,000 base pairs) analyzed with maximum parsimony, maximum likelihood and Bayesian inference methods. We show that mustelids are consistently resolved with high nodal support into four major clades and three monotypic lineages. Using Bayesian dating techniques, we provide evidence that mustelids underwent two bursts of diversification that coincide with major paleoenvironmental and biotic changes that occurred during the Neogene and correspond with similar bursts of cladogenesis in other vertebrate groups. Biogeographical analyses indicate that most of the extant diversity of mustelids originated in Eurasia and mustelids have colonized Africa, North America and South America on multiple occasions.

CONCLUSION

Combined with information from the fossil record, our phylogenetic and dating analyses suggest that mustelid diversification may have been spurred by a combination of faunal turnover events and diversification at lower trophic levels, ultimately caused by climatically driven environmental changes. Our biogeographic analyses show Eurasia as the center of origin of mustelid diversity and that mustelids in Africa, North America and South America have been assembled over time largely via dispersal, which has important implications for understanding the ecology of mustelid communities.

Metabolic rate does not calibrate the molecular clock

Rates of molecular evolution vary widely among lineages, but the causes of this variation remain poorly understood. It has been suggested that mass-specific metabolic rate may be one of the key factors determining the rate of molecular evolution, and that it can be used to derive "corrected" molecular clocks. However, previous studies have been hampered by a paucity of mass-specific metabolic rate data and have been largely limited to vertebrate taxa. Using mass-specific metabolic rate measurements and DNA sequence data for >300 metazoan species for 12 different genes, we find no evidence that mass-specific metabolic rate drives substitution rates. The mechanistic basis of the metabolic rate hypothesis is discussed in light of these findings.

Dates from the molecular clock: how wrong can we be?

Large discrepancies have been found in dates of evolutionary events obtained using the molecular clock. Twofold differences have been reported between the dates estimated from molecular data and those from the fossil record; furthermore, different molecular methods can give dates that differ 20-fold. New software attempts to incorporate appropriate allowances for this uncertainty into the calculation of the accuracy of date estimates. Here, we propose that these innovations represent welcome progress towards obtaining reliable dates from the molecular clock, but warn that they are currently unproven, given that the causes and pattern of the discrepancies are the subject of ongoing research. This research implies that many previous studies, even some of those using recently developed methods, might have placed too much confidence in their date estimates, and their conclusions might need to be revised.

Evolution of rattlesnakes (Viperidae; Crotalus) in the warm deserts of western North America shaped by Neogene vicariance and Quaternary climate change

During Pleistocene, the Laurentide ice sheet rearranged and diversified biotic distributions in eastern North America, yet had minimal physical impact in western North America where lineage diversification is instead hypothesized to result from climatic changes. If Pleistocene climatic fluctuations impacted desert species, the latter would reflect patterns of restricted gene flow concomitant with indications of demographic bottlenecks. Accordingly, molecular evidence for refugia should be present within these distributions and for subsequent range expansions as conditions improved. We sought answers to these questions by evaluating mitochondrial DNA (mtDNA) sequences from four species of rattlesnakes [Crotalus mitchellii (speckled rattlesnake), Crotalus cerastes (sidewinder), Crotalus tigris (tiger rattlesnake), Crotalus ruber (red diamond rattlesnake)] with distributions restricted to desert regions of southwestern North America. We inferred relationships using parsimony and maximum likelihood, tested intraspecific clades for population expansions, applied an isolation-with-migration model to determine bi-directional migration rates (m) among regions, and inferred divergence times for species and clades by applying a semiparametric penalized likelihood approach to our molecular data. Evidence for significant range expansion was present in two of eight regions in two species (Crotalus mitchellii pyrrhus, C. tigris region north). Two species (C. cerastes, C. mitchellii) showed a distribution concomitant with northward displacement of Baja California from mainland México, followed by vicariant separation into subclades. Effects of Pleistocene climate fluctuations were found in the distributions of all four species. Three regional diversification patterns were identified: (i) shallow genetic diversity that resulted from Pleistocene climatic events (C. tigris, C. ruber); (ii) deep Pleistocene divisions indicating allopatric segregation of subclades within refugia (C. mitchellii, C. cerastes); and (iii) lineage diversifications that extended to Pliocene or Late Miocene (C. mitchellii, C. cerastes). Clade-diversifying and clade-constraining effects impacted the four species of rattlesnakes unequally. We found relatively high levels of molecular diversification in the two most broadly distributed species (C. mitchellii, C. cerastes), and lower levels of genetic diversification in the two species (C. tigris, C. ruber) whose ranges are relatively more restricted. Furthermore, in several cases, the distributions of subspecies were not congruent with our molecular information. We suggest regional conservation perspectives for southwestern deserts cannot rely upon subspecies as biodiversity surrogates, but must instead employ a molecular and deep historical perspective as a primary mechanism to frame biodiversity reserves within this region.

[Molecular dating in the genomic era]

The comparison of DNA and protein sequences of extant species might be informative for reconstructing the chronology of evolutionary events on Earth. A phylogenetic tree inferred from molecular data directly depicts the evolutionary affinities of species and indirectly allows estimating the age of their origin and diversification. Molecular dating is achieved by assuming the molecular clock hypothesis, i.e., that the rate of change of nucleotide and amino acid sequences is on average constant over geological time. If paleontological calibrations are available, then absolute divergence times of species can be estimated. However, three major difficulties potentially hamper molecular dating : (1) a limited sample of genes and organisms, (2) a limited number of fossil references, and (3) pervasive variations of molecular evolutionary rates among genomes and species. To circumvent these problems, different solutions have been recently proposed. Larger data sets are built with more genes and more species sampled through the mining of an increasing number of genomes. Moreover, independent key fossils are identified to calibrate molecular clocks, and the uncertainty on their age is integrated in subsequent analyses. Finally, models of molecular rate variations are constructed, and incorporated in the so-called relaxed molecular clock approaches. As an illustration of these improvements, we mention that the debated age of the animal (bilaterian metazoans) diversification may have occurred between 642-761 million years ago (Mya), roughly 100 Ma before the Cambrian explosion. Among mammals, the initial diversification of major placental groups may have taken place around 100 Mya, well before the Cretaceous/Tertiary boundary marking the extinction of dinosaurs.

Intercontinental biogeography of subfamily Orontioideae (Symplocarpus, Lysichiton, and Orontium) of Araceae in eastern Asia and North America

Symplocarpus, Lysichiton, and Orontium (Orontioideae) are three of the few north temperate genera of the primarily tropical Araceae. Symplocarpus is disjunctly distributed in eastern Asia (3 spp.) and eastern North America (1 sp.); Lysichiton has an intercontinental discontinuous distribution in eastern Asia (1 sp.) and northwestern North America (1 sp.); and the monotypic Orontium is restricted to eastern North America. Phylogenetic analysis of the trnL-F and ndhF sequences supports (1) the monophyly of both Symplocarpus and Lysichiton, (2) the sister-group relationship of Symplocarpus and Lysichiton, and (3) the clade of Orontium, Symplocarpus, and Lysichiton. Although Symplocarpus shows a much wider disjunction than Lysichiton, the estimated divergence time of the former [4.49+/-1.69 or 6.88+/-4.18 million years ago (mya)] was similar to that of the latter (4.02+/-1.60 or 7.18+/-4.33 mya) based on the penalized likelihood and the Bayesian dating methods, respectively. Eastern Asia was suggested to be the ancestral area of the Symplocarpus-Lysichiton clade based on the dispersal-vicariance analysis. Our biogeographic results support independent migrations of Symplocarpus and Lysichiton across the Bering land bridge in the late Tertiary (Pliocene/late Miocene). Fossil evidence suggests Orontioideae dated back to the late Cretaceous in the temperate Northern Hemisphere (72 mya). The relative rate test shows similar substitution rates of the trnL-F sequences between the proto and the true aroids, although the latter has substantially higher species diversity. The proto Araceae perhaps suffered from a higher rate of extinction in the temperate zone associated with periods of climatic cooling in the Tertiary.

There is no universal molecular clock for invertebrates, but rate variation does not scale with body size

The existence of a universal molecular clock has been called into question by observations that substitution rates vary widely between lineages. However, increasing empirical evidence for the systematic effects of different life history traits on the rate of molecular evolution has raised hopes that rate variation may be predictable, potentially allowing the "correction" of the molecular clock. One such example is the body size trend observed in vertebrates; smaller species tend to have faster rates of molecular evolution. This effect has led to the proposal of general predictive models correcting for rate heterogeneity and has also been invoked to explain discrepancies between molecular and paleontological dates for explosive radiations in the fossil record. Yet, there have been no tests of an effect in any nonvertebrate taxa. In this study, we have tested the generality of the body size effect by surveying a wide range of invertebrate metazoan lineages. DNA sequences and body size data were collected from the literature for 330 species across five phyla. Phylogenetic comparative methods were used to investigate a relationship between average body size and substitution rate at both interspecies and interfamily comparison levels. We demonstrate significant rate variation in all phyla and most genes examined, implying a strict molecular clock cannot be assumed for the Metazoa. Furthermore, we find no evidence of any influence of body size on invertebrate substitution rates. We conclude that the vertebrate body size effect is a special case, which cannot be simply extrapolated to the rest of the animal kingdom.

Phylogeography of a subterranean amphipod reveals cryptic diversity and dynamic evolution in extreme environments

Extreme conditions in subsurface are suspected to be responsible for morphological convergences, and so to bias biodiversity assessment. Subterranean organisms are also considered as having poor dispersal abilities that in turn generate a large number of endemic species when habitat is fragmented. Here we test these general hypotheses using the subterranean amphipod Niphargus virei. All our phylogenetic analyses (Bayesian, maximum likelihood and distance), based on two independent genes (28S and COI), revealed the same tripartite structure. N. virei populations from Benelux, Jura region and the rest of France appeared as independent evolutionary units. Molecular rates estimated via global or Bayesian relaxed clock suggest that this split is at least 13 million years old and accredit the cryptic diversity hypothesis. Moreover, the geographical distribution of these lineages showed some evidence of recent dispersal through apparent vicariant barrier. In consequence, we argue that future analyses of evolution and biogeography in subsurface, or more generally in extreme environments, should consider dispersal ability as an evolving trait and morphology as a potentially biased marker.

Asymmetrical Evolution of Cytochrome bd Subunits

Functionally linked genes generally evolve at similar rates and the knowledge of this particular feature of genomic evolution has been used as the basis for the phylogenetic profiling method. We illustrate here an exception to this rule in the evolution of the cytochrome bd complex. This is a two-component oxidase complex, with the subunits I and II known to be widely present in bacteria. The subunits within the cytochrome bd complex are under the same evolutionary pressure and most likely behave in the same evolutionary manner. However, the sequence similarity of genes encoding subunit II varies considerably across species. Genes encoding subunit II evolve 1.2 times faster on most of the branches of their phylogeny than subunit I genes. Furthermore, the genes encoding subunit II in Oceanobacillus iheyensis, Bacillus halodurans, and Staphylococcus species do not have detectable homologues within E. coli due to their large divergence. Together, the two subunits of cytochrome bd reveal an interesting example of an asymmetric pattern of evolutionary change.

Relaxed molecular clocks for dating historical plant dispersal events

Age estimation from molecular sequences has emerged as a powerful tool for inferring when a plant lineage arrived in a particular area. Knowing the tenure of lineages within a region is key to understanding the evolution of traits, the evolution of biotic interactions, and the evolution of floras. New analytical methods model change in substitution rates along individual branches of a phylogenetic tree by combining molecular data with time constraints, usually from fossils. These "relaxed clock" approaches can be applied to several gene regions that need not all have the same substitution rates, and they can also incorporate multiple simultaneous fossil calibrations. Since 1995, at least 100 plant biogeographic studies have used molecular-clock dating, and about a fifth has used relaxed clocks. Many of these report evidence of long-distance dispersal. Meta-analyses of studies from the same geographic region can identify directional biases because of prevailing wind or water currents and the relative position and size of landmasses.

An Empirical Assessment of Long-Branch Attraction Artefacts in Deep Eukaryotic Phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.

Molecular clocks: four decades of evolution

During the past four decades, the molecular-clock hypothesis has provided an invaluable tool for building evolutionary timescales, and has served as a null model for testing evolutionary and mutation rates in different species. Molecular clocks have also influenced the development of theories of molecular evolution. As DNA-sequencing technologies have progressed, the use of molecular clocks has increased, with a profound effect on our understanding of the temporal diversification of species and genomes.

Dissection of phylogenetic relationships among 19 rapidly growing Mycobacterium species by 16S rRNA, hsp65, sodA, recA and rpoB gene sequencing

The current classification of non-pigmented and late-pigmenting rapidly growing mycobacteria (RGM) capable of producing disease in humans and animals consists primarily of three groups, the Mycobacterium fortuitum group, the Mycobacterium chelonae-abscessus group and the Mycobacterium smegmatis group. Since 1995, eight emerging species have been tentatively assigned to these groups on the basis of their phenotypic characters and 16S rRNA gene sequence, resulting in confusing taxonomy. In order to assess further taxonomic relationships among RGM, complete sequences of the 16S rRNA gene (1483-1489 bp), rpoB (3486-3495 bp) and recA (1041-1056 bp) and partial sequences of hsp65 (420 bp) and sodA (441 bp) were determined in 19 species of RGM. Phylogenetic trees based upon each gene sequence, those based on the combined dataset of the five gene sequences and one based on the combined dataset of the rpoB and recA gene sequences were then compared using the neighbour-joining, maximum-parsimony and maximum-likelihood methods after using the incongruence length difference test. Combined datasets of the five gene sequences comprising nearly 7000 bp and of the rpoB+recA gene sequences comprising nearly 4600 bp distinguished six phylogenetic groups, the M. chelonae-abscessus group, the Mycobacterium mucogenicum group, the M. fortuitum group, the Mycobacterium mageritense group, the Mycobacterium wolinskyi group and the M. smegmatis group, respectively comprising four, three, eight, one, one and two species. The two protein-encoding genes rpoB and recA improved meaningfully the bootstrap values at the nodes of the different groups. The species M. mucogenicum, M. mageritense and M. wolinskyi formed new groups separated from the M. chelonae-abscessus, M. fortuitum and M. smegmatis groups, respectively. The M. mucogenicum group was well delineated, in contrast to the M. mageritense and M. wolinskyi groups. For phylogenetic organizations derived from the hsp65 and sodA gene sequences, the bootstrap values at the nodes of a few clusters were <70 %. In contrast, phylogenetic organizations obtained from the 16S rRNA, rpoB and recA genes were globally similar to that inferred from combined datasets, indicating that the rpoB and recA genes appeared to be useful tools in addition to the 16S rRNA gene for the investigation of evolutionary relationships among RGM species. Moreover, rpoB gene sequence analysis yielded bootstrap values higher than those observed with recA and 16S rRNA genes. Also, molecular signatures in the rpoB and 16S rRNA genes of the M. mucogenicum group showed that it was a sister group of the M. chelonae-abscessus group. In this group, M. mucogenicum ATCC 49650(T) was clearly distinguished from M. mucogenicum ATCC 49649 with regard to analysis of the five gene sequences. This was in agreement with phenotypic and biochemical characteristics and suggested that these strains are representatives of two closely related, albeit distinct species.

Initial Diversification of Living Amphibians Predated the Breakup of Pangaea

The origin and divergence of the three living orders of amphibians (Anura, Caudata, Gymnophiona) and their main lineages are one of the most hotly debated topics in vertebrate evolution. Here, we present a robust molecular phylogeny based on the nuclear RAG1 gene as well as results from a variety of alternative independent molecular clock calibrations. Our analyses suggest that the origin and early divergence of the three living amphibian orders dates back to the Palaeozoic or early Mesozoic, before the breakup of Pangaea, and soon after the divergence from lobe-finned fishes. The resulting new biogeographic scenario, age estimate, and the inferred rapid divergence of the three lissamphibian orders may account for the lack of fossils that represent plausible ancestors or immediate sister taxa of all three orders and the heretofore paradoxical distribution of some amphibian fossil taxa. Furthermore, the ancient and rapid radiation of the three lissamphibian orders likely explains why branch lengths connecting their early nodes are particularly short, thus rendering phylogenetic inference of implicated relationships especially difficult.

Genome-scale evidence of the nematode-arthropod clade

The most extensive phylogenetic analysis carried out to date, including 11 complete genomes, is shown to support the Ecdysozoa hypothesis in the open-ended debate of the Coelomata-Ecdysozoa evolutionary problem.

Background

The issue of whether coelomates form a single clade, the Coelomata, or whether all animals that moult an exoskeleton (such as the coelomate arthropods and the pseudocoelomate nematodes) form a distinct clade, the Ecdysozoa, is the most puzzling issue in animal systematics and a major open-ended subject in evolutionary biology. Previous single-gene and genome-scale analyses designed to resolve the issue have produced contradictory results. Here we present the first genome-scale phylogenetic evidence that strongly supports the Ecdysozoa hypothesis.

Results

Through the most extensive phylogenetic analysis carried out to date, the complete genomes of 11 eukaryotic species have been analyzed in order to find homologous sequences derived from 18 human chromosomes. Phylogenetic analysis of datasets showing an increased adjustment to equal evolutionary rates between nematode and arthropod sequences produced a gradual change from support for Coelomata to support for Ecdysozoa. Transition between topologies occurred when fast-evolving sequences of Caenorhabditis elegans were removed. When chordate, nematode and arthropod sequences were constrained to fit equal evolutionary rates, the Ecdysozoa topology was statistically accepted whereas Coelomata was rejected.

Conclusions

The reliability of a monophyletic group clustering arthropods and nematodes was unequivocally accepted in datasets where traces of the long-branch attraction effect were removed. This is the first phylogenomic evidence to strongly support the 'moulting clade' hypothesis.

Sociality and the Rate of Molecular Evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.