Modelling the ecological-functional diversification of marine Metazoa on geological time scales

The ecological traits and functional capabilities of marine animals have changed significantly since their origin in the late Precambrian. These changes can be analysed quantitatively using multi-dimensional parameter spaces in which the ecological lifestyles of species are represented by particular combinations of parameter values. Here, we present models that describe the filling of this multi-dimensional 'ecospace' by ecological lifestyles during metazoan diversification. These models reflect varying assumptions about the processes that drove ecological diversification; they contrast diffusive expansion with driven expansion and niche conservatism with niche partitioning. Some models highlight the importance of interactions among organisms (ecosystem engineering and predator-prey escalation) in promoting new lifestyles or eliminating existing ones. These models reflect processes that were not mutually exclusive; rigorous analyses will continue to reveal their applicability to episodes in metazoan history.     

The evolution of South American endemic canids: a history of rapid diversification and morphological parallelism

The origin of endemic South American canid fauna has been traditionally linked with the rise of the Isthmus of Panama, suggesting that diversification of the dog fauna on this continent occurred very rapidly. Nevertheless, despite its obvious biogeographic appeal, the tempo of Canid evolution in South America has never been studied thoroughly. This issue can be suitably tackled with the inference of a molecular timescale. In this study, using a relaxed molecular clock method, we estimated that the most recent common ancestor of South American canids lived around 4 Ma, whereas all other splits within the clade occurred after the rise of the Panamanian land bridge. We suggest that the early diversification of the ancestors of the two main lineages of South American canids may have occurred in North America, before the Great American Interchange. Moreover, a concatenated morphological and molecular analysis put some extinct canid species well within the South American radiation, and shows that the dental adaptations to hypercarnivory evolved only once in the South American clade.     

Molecular evidence for the diversification of extant lichens in the late cretaceous and tertiary

A molecular clock based on ITS sequence data from the lichen genera Biatora and Phyllopsora is calibrated with the help of paleoclimatic data and evidence of forest history. The clock indicates that diversification within Biatora started as early as in the Late Cretaceous and took place during periods of climatic cooling, when new types of forest evolved and spread in the Northern Hemisphere. Arctic-alpine species of the genus appear to be of considerable age, dating back to the Late Eocene-Oligocene climatic cooling. By using calibrated phylogenies of epiphytic lichens it may become possible to date many paleoenvironmental events, for which little fossil evidence exists.     

A general comparison of relaxed molecular clock models

Several models have been proposed to relax the molecular clock in order to estimate divergence times. However, it is unclear which model has the best fit to real data and should therefore be used to perform molecular dating. In particular, we do not know whether rate autocorrelation should be considered or which prior on divergence times should be used. In this work, we propose a general bench mark of alternative relaxed clock models. We have reimplemented most of the already existing models, including the popular lognormal model, as well as various prior choices for divergence times (birth-death, Dirichlet, uniform), in a common Bayesian statistical framework. We also propose a new autocorrelated model, called the "CIR" process, with well-defined stationary properties. We assess the relative fitness of these models and priors, when applied to 3 different protein data sets from eukaryotes, vertebrates, and mammals, by computing Bayes factors using a numerical method called thermodynamic integration. We find that the 2 autocorrelated models, CIR and lognormal, have a similar fit and clearly outperform uncorrelated models on all 3 data sets. In contrast, the optimal choice for the divergence time prior is more dependent on the data investigated. Altogether, our results provide useful guidelines for model choice in the field of molecular dating while opening the way to more extensive model comparisons.     

Stochastic models of molecular evolution and the estimation of phylogeny and rates of nucleotide substitution in the hominoid primates

We have estimated phylogenetic patterns and rates of nucleotide substitution in the hominoid primates using two different probabilistic models of molecular evolution as applied to three different data sets of nucleic acid sequences. The orang-utan was found to be the out-group of the other hominoids examined. Within the African apes and human clade the sister-group relationship of chimpanzee and human was found to be statistically the best, although the magnitude of the error estimates (a reflection of random statistical fluctuations) makes this conclusion tentative. The ψν-globin data sets were found to be statistically the most consistent and gave estimates of the times of divergence of chimpanzee and human from gorilla and of chimpanzee from human as 7·7 ± 1·5 Ma (Millions of years ago) and 7·4 ± 1·5 Ma respectively, although the speculative nature of these estimates is emphasized. In all cases the calibration point was the assumed divergence of the orang-utan from the remaining hominoids at 14·5 Ma. There was no statistically significant evidence of a slowdown in nucleotide substitution rate for the human lineage, or among the hominoids as a whole with respect to the Old and New World monkeys. We advocate the continued use and development of stochastic models of molecular evolution as a basis for phylogenetic estimation. On this basis one can choose between competing hypotheses of relationship in a statistical manner and can provide estimates of the errors involved in such estimations. The assumptions of all stochastic models are open to test and future refinement.     

Evolutionary origins of germline segregation in Metazoa: evidence for a germ stem cell lineage in the coral Orbicella faveolata (Cnidaria,Anthozoa)

The ability to segregate a committed germ stem cell (GSC) lineage distinct from somatic cell lineages is a characteristic of bilaterian Metazoans. However, the occurrence of GSC lineage specification in basally branching Metazoan phyla, such as Cnidaria, is uncertain. Without an independently segregated GSC lineage, germ cells and their precursors must be specified throughout adulthood from continuously dividing somatic stem cells, generating the risk of propagating somatic mutations within the individual and its gametes. To address the potential for existence of a GSC lineage in Anthozoa, the sister-group to all remaining Cnidaria, we identified moderate- to high-frequency somatic mutations and their potential for gametic transfer in the long-lived coral Orbicella faveolata (Anthozoa, Cnidaria) using a 2b-RAD sequencing approach. Our results demonstrate that somatic mutations can drift to high frequencies (up to 50%) and can also generate substantial intracolonial genetic diversity. However, these somatic mutations are not transferable to gametes, signifying the potential for an independently segregated GSC lineage in O. faveolata. In conjunction with previous research on germ cell development in other basally branching Metazoan species, our results suggest that the GSC system may be a Eumetazoan characteristic that evolved in association with the emergence of greater complexity in animal body plan organization and greater specificity of stem cell functions.     

Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes

With about 24,000 extant species, teleosts are the largest group of vertebrates. They constitute more than 99% of the ray-finned fishes (Actinopterygii) that diverged from the lobe-finned fish lineage (Sarcopterygii) about 450 MYA. Although the role of genome duplication in the evolution of vertebrates is now established, its role in structuring the teleost genomes has been controversial. At least two hypotheses have been proposed: a whole-genome duplication in an ancient ray-finned fish and independent gene duplications in different lineages. These hypotheses are, however, based on small data sets and lack adequate statistical and phylogenetic support. In this study, we have made a systematic comparison of the draft genome sequences of Fugu and humans to identify paralogous chromosomal regions ("paralogons") in the Fugu that arose in the ray-finned fish lineage ("fish-specific"). We identified duplicate genes in the Fugu by phylogenetic analyses of the Fugu, human, and invertebrate sequences. Our analyses provide evidence for 425 fish-specific duplicate genes in the Fugu and show that at least 6.6% of the genome is represented by fish-specific paralogons. We estimated the ages of Fugu duplicate genes and paralogons using the molecular clock. Remarkably, the ages of duplicate genes and paralogons are clustered, with a peak around 350 MYA. These data strongly suggest a whole-genome duplication event early during the evolution of ray-finned fishes, probably before the origin of teleosts.     

Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds

We implement a Bayesian Markov chain Monte Carlo algorithm for estimating species divergence times that uses heterogeneous data from multiple gene loci and accommodates multiple fossil calibration nodes. A birth-death process with species sampling is used to specify a prior for divergence times, which allows easy assessment of the effects of that prior on posterior time estimates. We propose a new approach for specifying calibration points on the phylogeny, which allows the use of arbitrary and flexible statistical distributions to describe uncertainties in fossil dates. In particular, we use soft bounds, so that the probability that the true divergence time is outside the bounds is small but nonzero. A strict molecular clock is assumed in the current implementation, although this assumption may be relaxed. We apply our new algorithm to two data sets concerning divergences of several primate species, to examine the effects of the substitution model and of the prior for divergence times on Bayesian time estimation. We also conduct computer simulation to examine the differences between soft and hard bounds. We demonstrate that divergence time estimation is intrinsically hampered by uncertainties in fossil calibrations, and the error in Bayesian time estimates will not go to zero with increased amounts of sequence data. Our analyses of both real and simulated data demonstrate potentially large differences between divergence time estimates obtained using soft versus hard bounds and a general superiority of soft bounds. Our main findings are as follows. (1) When the fossils are consistent with each other and with the molecular data, and the posterior time estimates are well within the prior bounds, soft and hard bounds produce similar results. (2) When the fossils are in conflict with each other or with the molecules, soft and hard bounds behave very differently; soft bounds allow sequence data to correct poor calibrations, while poor hard bounds are impossible to overcome by any amount of data. (3) Soft bounds eliminate the need for "safe" but unrealistically high upper bounds, which may bias posterior time estimates. (4) Soft bounds allow more reliable assessment of estimation errors, while hard bounds generate misleadingly high precisions when fossils and molecules are in conflict.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

Biogeography and evolution of the Galapagos: integration of the biological and geological evidence

Biogeographic tracks are mapped for Galapagos endemics representing 25 plant and animal taxa and including organisms with good and poor means of dispersal. These patterns confirm standard biogeographic tracks linking Galapagos with Central America, western North and South America, the Caribbean, Asia and Australasia. Discovery of the Galapagos Gore in the 1970s corroborates the biogeographic prediction for a major tectonic centre associated with the Galapagos. The biogeographic model developed by Croizat in 1958 of Galapagos colonization involving an ancestral biota inhabiting eastern Pacific geosynclinal forelands is congruent with plate tectonic models supporting a Pacific island arc origin for western American terranes. American relatives of Galapagos endemics may have originated within an eastern Pacific paleogeography rather than representing centres of origin for dispersal to the Galapagos. Galapagos colonization by an eastern Pacific biota between late Cretaceous and mid-Tertiary has significant implications for understanding the tempo and mode for both the origins of island biota and general models of evolutionary differentiation. Popular assertions that overwater dispersal represents the only viable origin for the entire Galapagos biota is no longer biogeographically or geologically tenable.     

Phylogeography of the piranha genera Serrasalmus and Pygocentrus: implications for the diversification of the Neotropical ichthyofauna

The phylogenetic relationships within the piranhas were assessed using mitochondrial sequences with the aim of testing several hypotheses proposed to explain the origin of Neotropical diversity (palaeogeography, hydrogeology and museum hypotheses). Sequences of the ribosomal 16S gene (510 bp) and control region (980 bp) were obtained from 15 localities throughout the main South American rivers for 21 of the 28 extant piranha species. The results indicate that the genus Serrasalmus is monophyletic and comprises three major clades. The phylogeographical analyses of these clades allowed the identification of five vicariant events, extensive dispersal and four lineage duplications suggesting the occurrence of sympatric speciation. Biogeographical patterns are consistent with the prediction made by the museum hypothesis that lineages from the Precambrian shields are older than those from the lowlands of the Amazon. The vicariant events inferred here match the distribution of the palaeoarches and several postdispersal speciation events are identified, thereby matching the predictions of the palaeogeography and hydrogeology hypotheses, respectively. Molecular clock calibration of the control region sequences indicates that the main lineages differentiated from their most recent common ancestor at 9 million years ago in the proto Amazon-Orinoco and the present rate of diversification is the highest reported to date for large carnivorous Characiformes. The present results emphasize that an interaction among geology, sea-level changes, and hydrography created opportunities for cladogenesis in the piranhas at different temporal and geographical scales.     

On the rate of DNA sequence evolution inDrosophila

Summary Analysis of the rate of nucleotide substitution at silent sites inDrosophila genes reveals three main points. First, the silent rate varies (by a factor of two) among nuclear genes; it is inversely related to the degree of codon usage bias, and so selection among synonymous codons appears to constrain the rate of silent substitution in some genes. Second, mitochondrial genes may have evolved only as fast as nuclear genes with weak codon usage bias (and two times faster than nuclear genes with high codon usage bias); this is quite different from the situation in mammals where mitochondrial genes evolve approximately 5–10 times faster than nuclear genes. Third, the absolute rate of substitution at silent sites in nuclear genes inDrosophila is about three times hihger than the average silent rate in mammals.     

Accuracy of rate estimation using relaxed-clock models with a critical focus on the early metazoan radiation

In recent years, a number of phylogenetic methods have been developed for estimating molecular rates and divergence dates under models that relax the molecular clock constraint by allowing rate change throughout the tree. These methods are being used with increasing frequency, but there have been few studies into their accuracy. We tested the accuracy of several relaxed-clock methods (penalized likelihood and Bayesian inference using various models of rate change) using nucleotide sequences simulated on a nine-taxon tree. When the sequences evolved with a constant rate, the methods were able to infer rates accurately, but estimates were more precise when a molecular clock was assumed. When the sequences evolved under a model of auto-correlated rate change, rates were accurately estimated using penalized likelihood and by Bayesian inference using lognormal and exponential models of rate change, while other models did not perform as well. When the sequences evolved under a model of uncorrelated rate change, only Bayesian inference using an exponential rate model performed well. Collectively, the results provide a strong recommendation for using the exponential model of rate change if a conservative approach to divergence time estimation is required. A case study is presented in which we use a simulation-based approach to examine the hypothesis of elevated rates in the Cambrian period, and it is found that these high rate estimates might be an artifact of the rate estimation method. If this bias is present, then the ages of metazoan divergences would be systematically underestimated. The results of this study have implications for studies of molecular rates and divergence dates.     

The evolution of ecomorphological traits within the Abrothrichini (Rodentia: Sigmodontinae): a Bayesian phylogenetics approach

The generally accepted hypothesis regarding the origin of fossorial mammals proposes adaptive convergence from open environments towards the use of subterranean environments. We evaluated this hypothesis for South American mole-mice using conventional and Bayesian frameworks, with independent evidence. By using a molecular approach based on Cytochrome b and IRBP sequences, we evaluated phylogenetic relationships, time of origin, the ancestral trait of fossoriality, and ancestral distributions of species belonging to the Andean Clade (Rodentia: Sigmodontinae). Our results indicate that the Andean Clade is highly sustained; with one clade grouping all fossorial forms and another grouping all cursorial species. We hypothesized that fossoriality originated in the Miocene/Pliocene transition, in the Temperate Forests of southern South America. We conclude that the origin of fossorial ecomorphological traits did not necessarily occur under a general model of open environments, the origin of these traits depends on the ecological-historical relationship of the taxon with the environment.     

Quantifying the impact of protein tertiary structure on molecular evolution

To investigate the evolutionary impact of protein structure, the experimentally determined tertiary structure and the protein-coding DNA sequence were collected for each of 1,195 genes. These genes were studied via a model of sequence change that explicitly incorporates effects on evolutionary rates due to protein tertiary structure. In the model, these effects act via the solvent accessibility environments and pairwise amino acid interactions that are induced by tertiary structure. To compare the hypotheses that structure does and does not have a strong influence on evolution, Bayes factors were estimated for each of the 1,195 sequences. Most of the Bayes factors strongly support the hypothesis that protein structure affects protein evolution. Furthermore, both solvent accessibility and pairwise interactions among amino acids are inferred to have important roles in protein evolution. Our results also indicate that the strength of the relationship between tertiary structure and evolution has a weak but real correlation to the annotation information in the Gene Ontology database. Although their influences on rates of evolution vary among protein families, we find that the mean impacts of solvent accessibility and pairwise interactions are about the same.     

Analysis of the overdispersed clock in the short-term evolution of hepatitis C virus: Using the E1/E2 gene sequences to infer infection dates in a single source outbreak

The assumption of a molecular clock for dating events from sequence information is often frustrated by the presence of heterogeneity among evolutionary rates due, among other factors, to positively selected sites. In this work, our goal is to explore methods to estimate infection dates from sequence analysis. One such method, based on site stripping for clock detection, was proposed to unravel the clocklike molecular evolution in sequences showing high variability of evolutionary rates and in the presence of positive selection. Other alternatives imply accommodating heterogeneity in evolutionary rates at various levels, without eliminating any information from the data. Here we present the analysis of a data set of hepatitis C virus (HCV) sequences from 24 patients infected by a single individual with known dates of infection. We first used a simple criterion of relative substitution rate for site removal prior to a regression analysis. Time was regressed on maximum likelihood pairwise evolutionary distances between the sequences sampled from the source individual and infected patients. We show that it is indeed the fastest evolving sites that disturb the molecular clock and that these sites correspond to positively selected codons. The high computational efficiency of the regression analysis allowed us to compare the site-stripping scheme with random removal of sites. We demonstrate that removing the fast-evolving sites significantly increases the accuracy of estimation of infection times based on a single substitution rate. However, the time-of-infection estimations improved substantially when a more sophisticated and computationally demanding Bayesian method was used. This method was used with the same data set but keeping all the sequence positions in the analysis. Consequently, despite the distortion introduced by positive selection on evolutionary rates, it is possible to obtain quite accurate estimates of infection dates, a result of especial relevance for molecular epidemiology studies.     

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.     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.