Macroevolutionary patterns in Rhynchocephalia: is the tuatara (Sphenodon punctatus) a living fossil?

The tuatara, Sphenodon punctatus, known from 32 small islands around New Zealand, has often been noted as a classic ‘living fossil’ because of its apparently close resemblance to its Mesozoic forebears and because of a long, low‐diversity history. This designation has been disputed because of the wide diversity of Mesozoic forms and because of derived adaptations in living Sphenodon. We provide a testable definition for ‘living fossils’ based on a slow rate of lineage evolution and a morphology close to the centroid of clade morphospace. We show that through their history since the Triassic, rhynchocephalians had heterogeneous rates of morphological evolution and occupied wide morphospaces during the Triassic and Jurassic, and these then declined in the Cretaceous. In particular, we demonstrate that the extant tuatara underwent unusually slow lineage evolution, and is morphologically conservative, being located near the centre of the morphospace for all Rhynchocephalia.     

Temperature‐dependent mutational robustness can explain faster molecular evolution at warm temperatures,affecting speciation rate and global patterns of species diversity

Distribution of species across the Earth shows strong latitudinal and altitudinal gradients with the number of species decreasing with declining temperatures. While these patterns have been recognized for well over a century, the mechanisms generating and maintaining them have remained elusive. Here, we propose a mechanistic explanation for temperature‐dependent rates of molecular evolution that can influence speciation rates and global biodiversity gradients. Our hypothesis is based on the effects of temperature and temperature‐adaptation on stability of proteins and other catalytic biomolecules. First, due to the nature of physical forces between biomolecules and water, stability of biomolecules is maximal around + 20°C and decreases as temperature either decreases or increases. Second, organisms that have adapted to cold temperatures have evolved especially flexible (but unstable) proteins to facilitate catalytic reactions in cold, where molecular movements slow down. Both these effects should result in mutations being on average more detrimental at cold temperatures (i.e. lower mutational robustness in cold). At high temperatures, destabilizing water–biomolecule interactions, and the need to maintain structures that withstand heat denaturation, should decrease mutational robustness similarly. Decreased mutational robustness at extreme temperatures will slow down molecular evolution, as a larger fraction of new mutations will be removed by selection. Lower mutational robustness may also select for reduced mutation rates, further slowing down the rate of molecular evolution. As speciation requires the evolution of epistatic incompatibilities that prevent gene flow among incipient species, slow rate of molecular evolution at extreme temperatures will directly slow down the rate at which new species arise. The proposed mechanism can thus explain why molecular evolution is faster at warm temperatures, contributing to higher speciation rate and elevated species richness in environments characterized by stable and warm temperatures.     

Erratic rates of molecular evolution and incongruence of fossil and molecular divergence time estimates in Ostracoda (Crustacea)

Dating evolutionary origins of taxa is essential for understanding rates and timing of evolutionary events, often inciting intense debate when molecular estimates differ from first fossil appearances. For numerous reasons, ostracods present a challenging case study of rates of evolution and congruence of fossil and molecular divergence time estimates. On the one hand, ostracods have one of the densest fossil records of any metazoan group. However, taxonomy of fossil ostracods is controversial, owing at least in part to homoplasy of carapaces, the most commonly fossilized part. In addition, rates of evolution are variable in ostracods. Here, we report evidence of extreme variation in the rate of molecular evolution in different ostracod groups. This rate is significantly elevated in Halocyprid ostracods, a widespread planktonic group, consistent with previous observations that planktonic groups show elevated rates of molecular evolution. At the same time, the rate of molecular evolution is slow in the lineage leading to Manawa staceyi, a relict species that we estimate diverged approximately 500 million years ago from its closest known living relative. We also report multiple cases of significant incongruence between fossil and molecular estimates of divergence times in Ostracoda. Although relaxed clock methods improve the congruence of fossil and molecular divergence estimates over strict clock models, incongruence is present regardless of method. We hypothesize that this observed incongruence is driven largely by problems with taxonomy of fossil Ostracoda. Our results illustrate the difficulty in consistently estimating lineage divergence times, even in the presence of a voluminous fossil record.     

Molecular clocks in reptiles: life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between body size and substitution rate for many molecular data sets. Both the generality and the cause of the negative body size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size, age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles, lizards, snakes, crocodile, and tuatara. Although this study was limited by the number of comparisons for which both sequence and life-history data were available, the results suggest that a negative body size trend in rate of molecular evolution may be a general feature of reptile molecular evolution, consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.     

Evolutionary rate variation in Old World monkeys

We analysed over 8 million base pairs of bacterial artificial chromosome-based sequence alignments of four Old World monkeys and the human genome. Our findings are as follows. (i) Genomic divergences among several Old World monkeys mirror those between well-studied hominoids. (ii) The X-chromosome evolves slower than autosomes, in accord with ‘male-driven evolution’. However, the degree of male mutation bias is lower in Old World monkeys than in hominoids. (iii) Evolutionary rates vary significantly between lineages. The baboon branch shows a particularly slow molecular evolution. Thus, lineage-specific evolutionary rate variation is a common theme of primate genome evolution. (iv) In contrast to the overall pattern, mutations originating from DNA methylation exhibit little variation between lineages. Our study illustrates the potential of primates as a model system to investigate genome evolution, in particular to elucidate molecular mechanisms of substitution rate variation.     

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.     

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a ??molecular clock??. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50?years, evolutionary biologists have obtained an in-depth understanding of this clock??s nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to ??relaxed?? molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively ??slow?? sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage??s sequence neutral substitution rate, enabling species to leapfrog between phenotypic ??islands??. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.     

Slow molecular evolution in 18S rDNA, rbcL and nad5 genes of mosses compared with higher plants

The evolutionary potential of bryophytes (mosses, liverworts and hornworts) has been debated for decades. Fossil record and biogeographical distribution patterns suggest very slow morphological evolution and the retainment of several ancient traits since the split with vascular plants some 450 million years ago. Many have argued that bryophytes may evolve as rapidly as higher plants on the molecular level, but this hypothesis has not been tested so far. Here, it is shown that mosses have experienced significantly lower rates of molecular evolution than higher plants within 18S rDNA (nuclear), rbcL (chloroplast) and nad5 (mitochondrial) genes. Mosses are on an average evolving 2-3 times slower than ferns, gymnosperms and angiosperms; and also green algae seem to be evolving faster than nonvascular plants. These results support the observation of a general correlation between morphological and molecular evolutionary rates in plants and also show that mosses are 'evolutionary sphinxes' regarding both morphological and molecular evolutionary potential.     

Understanding the overdispersed molecular clock

Rates of molecular evolution at some protein-encoding loci are more irregular than expected under a simple neutral model of molecular evolution. This pattern of excessive irregularity in protein substitutions is often called the "overdispersed molecular clock" and is characterized by an index of dispersion, R(T) > 1. Assuming infinite sites, no recombination model of the gene R(T) is given for a general stationary model of molecular evolution. R(T) is shown to be affected by only three things: fluctuations that occur on a very slow time scale, advantageous or deleterious mutations, and interactions between mutations. In the absence of interactions, advantageous mutations are shown to lower R(T); deleterious mutations are shown to raise it. Previously described models for the overdispersed molecular clock are analyzed in terms of this work as are a few very simple new models. A model of deleterious mutations is shown to be sufficient to explain the observed values of R(T). Our current best estimates of R(T) suggest that either most mutations are deleterious or some key population parameter changes on a very slow time scale. No other interpretations seem plausible. Finally, a comment is made on how R(T) might be used to distinguish selective sweeps from background selection.     

Identification and partial characterization of α2-macroglobulin from the tuatara (Sphenodon punctatus)

α₂-Macroglobulin (α₂-M), a large molecular mass proteinase-binding protein, was identified in plasma from tuatara (Sphenodon), a rare reptile endemic to New Zealand. In this genus, α₂-M constitutes 11–13% of total plasma protein (∼2.2–3.9 mg/ml). Analysis of blood samples collected at approximately monthly intervals from individual tuatara indicated that the plasma level of α₂-M remains fairly constant. The subunits of tuatara α₂-M have an apparent molecular mass of ∼160 kDa as determined by SDS-polyacrylamide gel electrophoresis and the intact protein is an oligomer that contains inter-chain disulfide bonds. N-terminal sequence analyses of tuatara α₂-M revealed a distinct similarity to α-macroglobulins of other vertebrates and that at least two types of α₂-M subunits are present in plasma of tuatara.     

Evolutionary rates and species diversity in flowering plants

Genetic change is a necessary component of speciation, but the relationship between rates of speciation and molecular evolution remains unclear. We use recent phylogenetic data to demonstrate a positive relationship between species numbers and the rate of neutral molecular evolution in flowering plants (in both plastid and nuclear genes). Rates of protein and morphological evolution also correlate with the neutral substitution rate, but not with species numbers. Our findings reveal a link between the rate of neutral molecular change within populations and the evolution of species diversity.     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

Body size effects and rates of cytochrome b evolution in tube-nosed seabirds

Variation in rates of molecular evolution now appears to be widespread. The demonstration that body size is correlated with rates of molecular evolution suggests that physiological and ecological factors may be involved in molecular rate variation, but large-scale comparative studies are still lacking. Here, we use complete cytochrome b sequences from 85 species of tube-nosed seabirds (order Procellariiformes) and 5 outgroup species of penguins (order Sphenisciformes) to test for an association between body mass and rates of molecular evolution within the former avian order. Cladistic analysis of the 90 sequences estimates a phylogeny largely consistent with the traditional taxonomy of the Procellariiformes. The Diomedeidae, Procellariidae, and Pelecanoididae are monophyletic, while the Hydrobatidae are basal and paraphyletic. However, the two subfamilies within the Hydrobatidae (Hydrobatinae and Oceanitinae) are monophyletic. A likelihood ratio test detects significant deviation from clocklike evolution in our data. Using a sign test for an association between body mass and branch length in the seabird phylogeny, we find that larger taxa tend to have shorter terminal branch lengths than smaller taxa. This observation suggests that rates of mitochondrial DNA evolution are slower for larger taxa. Rate calibrations based on the fossil record reveal concordant body size effects. We interpret these results as evidence for a metabolic rate effect, as the species in this order exhibit large differences in metabolic rates, which are known to be highly correlated with body mass in this group. Our results support previous findings of body size effects and show that this effect can be significant even within a single avian order. This suggests that even lineage-specific molecular clocks may not be tenable if calibrations involve taxa with different metabolic rates.      

Why coelacanths are not ‘living fossils’

A series of recent studies on extant coelacanths has emphasised the slow rate of molecular and morphological evolution in these species. These studies were based on the assumption that a coelacanth is a ‘living fossil’ that has shown little morphological change since the Devonian, and they proposed a causal link between low molecular evolutionary rate and morphological stasis. Here, we have examined the available molecular and morphological data and show that: (i) low intra‐specific molecular diversity does not imply low mutation rate, (ii) studies not showing low substitution rates in coelacanth are often neglected, (iii) the morphological stability of coelacanths is not supported by paleontological evidence. We recall that intra‐species levels of molecular diversity, inter‐species genome divergence rates and morphological divergence rates are under different constraints and they are not necessarily correlated. Finally, we emphasise that concepts such as ‘living fossil’, ‘basal lineage’, or ‘primitive extant species’ do not make sense from a tree‐thinking perspective. Editor's suggested further reading in BioEssays Tree thinking for all biology: the problem with reading phylogenies as ladders of progress Abstract     

Nuclear mitochondrial pseudogenes as molecular outgroups for phylogenetically isolated taxa: a case study in Sphenodon

'Living fossil' taxa, by definition, have no close relatives, and therefore no outgroup to provide a root to phylogenetic trees. We identify and use a molecular outgroup in the sole extant lineage of sphenodontid reptiles, which separated from other reptiles 230 million years ago. We isolated and sequenced a partial nuclear copy of the mitochondrial cytochrome b gene. We confirm the copy is indeed not mitochondrial, is older than all extant mitochondrial copies in Sphenodon (tuatara), and is therefore useful as a molecular outgroup. Under phylogenetic analysis, the nuclear copy places the root of the tuatara mitochondrial gene tree between the northern and the southern (Cook Strait) groups of islands of New Zealand that are the last refugia for Sphenodon. This analysis supports a previous mid-point rooted mitochondrial gene tree. The mitochondrial DNA tree conflicts with allozyme analyses which place a Cook Strait population equidistant to all northern and other Cook Strait populations. This population on North Brother Island is the only natural population of extant S. guntheri; thus, we suggest that the current species designations of tuatara require further investigation.     

Seasonal and spatial dynamics of ectoparasite infestation of a threatened reptile, the tuatara (Sphenodon punctatus)

Abstract The conservation of threatened vertebrate species and their threatened parasites requires an understanding of the factors influencing their distribution and dynamics. This is particularly important for species maintained in conservation reserves at high densities, where increased contact among hosts could lead to increased rates of parasitism. The tuatara (Sphenodon punctatus) (Reptilia: Sphenodontia) is a threatened reptile that persists at high densities in forests (～ 2700 tuatara/ha) and lower densities in pastures and shrubland (< 200 tuatara/ha) on Stephens Island, New Zealand. We investigated the lifecycles and seasonal dynamics of infestation of two ectoparasites (the tuatara tick, Amblyomma sphenodonti, and trombiculid mites, Neotrombicula sp.) in a mark‐recapture study in three forest study plots from November 2004 to March 2007, and compared infestation levels among habitat types in March 2006. Tick loads were lowest over summer and peaked from late autumn (May) until early spring (September). Mating and engorgement of female ticks was highest over spring, and larval tick loads subsequently increased in early autumn (March). Nymphal tick loads increased in September, and adult tick loads increased in May. Our findings suggest the tuatara tick has a 2‐ or 3‐year lifecycle. Mite loads were highest over summer and autumn, and peaked in March. Prevalences (proportion of hosts infected) and densities (estimated number of parasites per hectare) of ticks were similar among habitats, but tick loads (parasites per host) were higher in pastures than in forests and shrub. The prevalence and density of mites was higher in forests than in pasture or shrub, but mite loads were similar among habitats. We suggest that a higher density of tuatara in forests may reduce the ectoparasite loads of individuals through a dilution effect. Understanding host–parasite dynamics will help in the conservation management of both the host and its parasites.     

Reconciling extreme branch length differences: decoupling time and rate through the evolutionary history of filmy ferns

The rate of molecular evolution is not constant across the Tree of Life. Characterizing rate discrepancies and evaluating the relative roles of time and rate along branches through the past are both critical to a full understanding of evolutionary history. In this study, we explore the interactions of time and rate in filmy ferns (Hymenophyllaceae), a lineage with extreme branch length differences between the two major clades. We test for the presence of significant rate discrepancies within and between these clades, and we separate time and rate across the filmy fern phylogeny to simultaneously yield an evolutionary time scale of filmy fern diversification and reconstructions of ancestral rates of molecular evolution. Our results indicate that the branch length disparity observed between the major lineages of filmy ferns is indeed due to a significant difference in molecular evolutionary rate. The estimation of divergence times reveals that the timing of crown group diversification was not concurrent for the two lineages, and the reconstruction of ancestral rates of molecular evolution points to a substantial rate deceleration in one of the clades. Further analysis suggests that this may be due to a genome-wide deceleration in the rate of nucleotide substitution.     

INCREASED RATES OF MOLECULAR EVOLUTION IN AN EQUATORIAL PLANT CLADE: AN EFFECT OF ENVIRONMENT OR PHYLOGENETIC NONINDEPENDENCE?

A recent study of environmental effects on rates of molecular evolution in the plant subgenus Mearnsia shows that species occurring in more equatorial latitudes have higher rates of substitution in rDNA sequences as compared to their more southerly congeners (Wright et al. 2003). However, we believe that the statistical approach employed by Wright et al. (2003) insufficiently accounts for the phylogenetic nonindependence of the species examined, given that all six equatorial species of Mearnsia form a clade. To distinguish between the effect of latitude and that of phylogenetic nonindependence, we have employed a variety of comparative approaches that use independent contrasts to test for an effect of environment across this entire subgenus. We find very little evidence for an effect of latitude on rate of molecular evolution using these approaches and believe that the shared evolutionary history of the clade is a plausible explanation of the apparent rate difference between equatorial and subequatorial Mearnsia species.     

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric‐based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.