In order to study the effect of mutation rate heterogeneity on patterns of DNA polymorphism, we simulated samples of DNA sequences with gamma-distributed nucleotide substitution rates in stationary and expanding populations. We find that recent population expansions and mutation rate heterogeneity have similar effects on several polymorphism indicators, like the shape and the mean of the observed pairwise difference distribution, or the number of segregating sites. The inferred size of population expansion thus appears overestimated if nucleotides have dissimilar substitution rates. Interestingly, population expansion and uneven mutation rates have contrasting effects on Tajima's D statistic when acting separately, and the consequence on the associated test of selective neutrality is investigated. The patterns of polymorphism of several human populations analyzed for the mitochondrial control region are examined, mainly showing the difficulty in quantifying the respective contribution of past demographic history and uneven mutation rates from a single sampled evolutionary process. However, substitution rates appear more heterogeneous in the second hypervariable segment of the control region than in the first segment.

The current paradigm holds that cyanobacteria, which evolved oxygenic photosynthesis more than 2 billion years ago, are still the major light harvesters driving primary productivity in open oceans. Here we show that tiny unicellular eukaryotes belonging to the photosynthetic lineage of the Haptophyta are dramatically diverse and ecologically dominant in the planktonic photic realm. The use of Haptophyta-specific primers and PCR conditions adapted for GC-rich genomes circumvented biases inherent in classical genetic approaches to exploring environmental eukaryotic biodiversity and led to the discovery of hundreds of unique haptophyte taxa in 5 clone libraries from subpolar and subtropical oceanic waters. Phylogenetic analyses suggest that this diversity emerged in Paleozoic oceans, thrived and diversified in the permanently oxygenated Mesozoic Panthalassa, and currently comprises thousands of ribotypic species, belonging primarily to low-abundance and ancient lineages of the “rare biosphere.” This extreme biodiversity coincides with the pervasive presence in the photic zone of the world ocean of 19′-hexanoyloxyfucoxanthin (19-Hex), an accessory photosynthetic pigment found exclusively in chloroplasts of haptophyte origin. Our new estimates of depth-integrated relative abundance of 19-Hex indicate that haptophytes dominate the chlorophyll a -normalized phytoplankton standing stock in modern oceans. Their ecologic and evolutionary success, arguably based on mixotrophy, may have significantly impacted the oceanic carbon pump. These results add to the growing evidence that the evolution of complex microbial eukaryotic cells is a critical force in the functioning of the biosphere.

The molecular clock, i.e., constancy of the rate of evolution over time, is commonly assumed in estimating divergence dates. However, this assumption is often violated and has drastic effects on date estimation. Recently, a number of attempts have been made to relax the clock assumption. One approach is to use maximum likelihood, which assigns rates to branches and allows the estimation of both rates and times. An alternative is the Bayes approach, which models the change of the rate over time. A number of models of rate change have been proposed. We have extended and evaluated models of rate evolution, i.e., the lognormal and its recent variant, along with the gamma, the exponential, and the Ornstein–Uhlenbeck processes. These models were first applied to a small hominoid data set, where an empirical Bayes approach was used to estimate the hyperparameters that measure the amount of rate variation. Estimation of divergence times was sensitive to these hyperparameters, especially when the assumed model is close to the clock assumption. The rate and date estimates varied little from model to model, although the posterior Bayes factor indicated the Ornstein-Uhlenbeck process outperformed the other models. To demonstrate the importance of allowing for rate change across lineages, this general approach was used to analyze a larger data set consisting of the 18S ribosomal RNA gene of 39 metazoan species. We obtained date estimates consistent with paleontological records, the deepest split within the group being about 560 million years ago. Estimates of the rates were in accordance with the Cambrian explosion hypothesis and suggested some more recent lineage-specific bursts of evolution.

Multicellular animals, or Metazoa, appear in the fossil records between 575 and 509 million years ago (MYA). At odds with paleontological evidence, molecular estimates of basal metazoan divergences have been consistently older than 700 MYA. However, those date estimates were based on the molecular clock hypothesis, which is almost always violated. To relax this hypothesis, we have implemented a Bayesian approach to describe the change of evolutionary rate over time. Analysis of 22 genes from the nuclear and the mitochondrial genomes under the molecular clock assumption produced old date estimates, similar to those from previous studies. However, by allowing rates to vary in time and by taking small species-sampling fractions into account, we obtained much younger estimates, broadly consistent with the fossil records. In particular, the date of protostome–deuterostome divergence was on average 582 ± 112 MYA. These results were found to be robust to specification of the model of rate change. The clock assumption thus had a dramatic effect on date estimation. However, our results appeared sensitive to the prior model of cladogenesis, although the oldest estimates (791 ± 246 MYA) were obtained under a suboptimal model. Bayes posterior estimates of evolutionary rates indicated at least one major burst of molecular evolution at the end of the Precambrian when protostomes and deuterostomes diverged. We stress the importance of assumptions about rates on date estimation and suggest that the large discrepancies between the molecular and fossil dates of metazoan divergences might partly be due to biases in molecular date estimation.

The use of genomic data and the rise of phylogenomics have radically changed our view of the eukaryotic tree of life at a high taxonomic level by identifying 4–6 “supergroups.” Yet, our understanding of the evolution of key innovations within each of these supergroups is limited because of poor species sampling relative to the massive diversity encompassed by each supergroup. Here we apply a multigene approach that incorporates a wide taxonomic diversity to infer the time line of the emergence of strategic evolutionary transitions in the haptophytes, a group of ecologically and biogeochemically significant marine protists that belong to the Chromalveolata supergroup. Four genes (SSU, LSU, tufA, and rbcL) were extensively analyzed under several Bayesian models to assess the robustness of the phylogeny, particularly with respect to 1) data partitioning; 2) the origin of the genes (host vs. endosymbiont); 3) across-site rate variation; and 4) across-lineage rate variation. We show with a relaxed clock analysis that the origin of haptophytes dates back to 824 million years ago (Ma) (95% highest probability density 1,031–637 Ma). Our dating results show that the ability to calcify evolved earlier than previously thought, between 329 and 291 Ma, in the Carboniferous period and that the transition from mixotrophy to autotrophy occurred during the same time period. Although these two transitions precede a habitat change of major diversities from coastal/neritic waters to the pelagic realm (291–243 Ma, around the Permian/Triassic boundary event), the emergence of calcification, full autotrophy, and oceanic lifestyle seem mutually independent.