Estimation of primate speciation dates using local molecular clocks

Incorporating gene-specific variation when inferring and evaluating optimal evolutionary tree topologies from multilocus sequence data

Because of the increase of genomic data, multiple genes are often available for the inference of phylogenetic relationships. The simple approach for combining multiple genes from the same taxon is to concatenate the sequences and then ignore the fact that different positions in the concatenated sequence came from different genes. Here, we discuss two criteria for inferring the optimal tree topology from data sets with multiple genes. These criteria are designed for multigene data sets where gene-specific evolutionary features are too important to ignore. One criterion is conventional and is obtained by taking the sum of log-likelihoods over all genes. The other criterion is obtained by dividing the log-likelihood for a gene by its sequence length and then taking the arithmetic mean over genes of these ratios. A similar strategy could be adopted with parsimony scores. The optimal tree is then declared to be the one for which the sum or the arithmetic mean is maximized. These criteria are justified within a two-stage hierarchical framework. The first level of the hierarchy represents gene-specific evolutionary features, and the second represents site-specific features for given genes. For testing significance of the optimal topology, we suggest a two-stage bootstrap procedure that involves resampling genes and then resampling alignment columns within resampled genes. An advantage of this procedure over concatenation is that it can effectively account for gene-specific evolutionary features. We discuss the applicability of the two-stage bootstrap idea to the Kishino-Hasegawa test and the Shimodaira-Hasegawa test.

Molecular evolution of recombination hotspots and highly recombining pseudoautosomal regions in hominoids

We examined the effects of recombination on the molecular evolution of noncoding regions in pseudoautosomal regions (PARs) and recombination hotspots in hominoids. The PAR-linked regions analyzed had on average longer branch lengths than those of the recombination hotspots. Moreover, contrary to previous observations, we found no correlation between recombination rate and silent site divergence in our data set and little change in the GC content during recent hominoid evolution. This suggests that the current rate of recombination is not a good indicator of the past rates of recombination for these highly recombining regions. Furthermore, human recombination hotspots show increased AT to GC substitutions in the human lineage, while no such pattern is detected for PAR-linked regions. Together, these observations suggest that recombination hotspots in hominoids are transient in the evolutionary time-scale. Interestingly, the 16p13.3 recombination hotspot locus violates a local molecular clock, though the locus appears to be noncoding and should evolve neutrally. We hypothesize that sudden changes in recombination rate have caused the changes in substitution rate at this locus.

The Effect of Paralogous Lineages on the Application of Reconciliation Analysis by Cophylogeny Mapping

Paralogy defines similarity caused by duplication rather than common descent and is well known in the case of paralogous gene copies within a single genome. The term is here extended to paralogous lineages of associates within a single host. The phylogenies of four genera within the Herpesviridae were reconciled with host phylogenies using cophylogenetic mapping. The observed correspondence for each pair of phylogenies was evaluated through randomization of the viral phylogeny and demonstrated to be greater than expected by chance. A simulation study was then carried out to assess the influence of paralogous lineages on the efficacy of reconciliation analysis. Combining viral taxa from different genera that infected common hosts introduced incongruence into the cophylogenies and reduced both the minimum and maximum observed number of codivergence events relative to the initial analysis of orthologous clades. However, at an average sample size this did not alter the fundamental significance of observed correspondence. With smaller sample sizes, the number of orthologous taxa selected at random from the pool of taxa was reduced. False-negative results then increased in proportion from 0.02 to 0.33. These results demonstrated that reconciliation analysis is robust under conditions of paralogy at "normal" sample sizes but is adversely affected by a combination of paralogy and low sample size. Consideration of phylogenies for Papillomavirus, Atadenovirus, and Mastadenovirus suggest that paralogous lineages may be a widespread phenomenon among DNA viruses and that duplication irrespective of host speciation is an important cause of viral diversification.

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.

ANGIOSPERM DIVERGENCE TIMES: THE EFFECT OF GENES, CODON POSITIONS, AND TIME CONSTRAINTS

An understanding of the evolution of modern terrestrial ecosystems requires an understanding of the dynamics associated with angiosperm evolution, including the timing of their origin and diversification into their extraordinary present-day diversity. Molecular estimates of angiosperm age have varied widely, and many substantially predate the Early Cretaceous fossil appearance of the group. In this study, the effect of different genes, codon positions, and chronological constraints on node ages are examined on divergence time estimates across seed plants, with a special focus on angiosperms. Penalized likelihood was used to estimate divergence times on a phylogenetic hypothesis for seed plants derived from Bayesian analysis, with branch lengths estimated with maximum likelihood. The plastid genes atpB, psaA, psbB, and rbcL were used individually and in combination, using first and second, third, and the three codon positions, including and excluding age constraints on 20 nodes derived from a critical examination of the land-plant fossil record. The optimal level of rate smoothing according to each unconstrained and constrained dataset was obtained with penalized likelihood. Tests for a molecular clock revealed significantly unclocklike rates in all datasets. Addition of fossil constraints resulted in even greater departures from constancy. Consistently with significant deviations from a clock, estimated optimal smoothing values were low, but a strict correlation between rate heterogeneity and optimal smoothing value was not found. Age estimates for nodes across the phylogeny varied, sometimes substantially, with gene and codon position. Nevertheless, estimates based on the four concatenated genes are very similar to the mean of the four individual gene estimates. For any given node, unconstrained age estimates are more variable than constrained estimates and are frequently younger than well-substantiated fossil members of the clade. Constrained estimates of ages of clades are older than unconstrained estimates and oldest fossil representatives, sometimes substantially so. Angiosperm age estimates decreased as rate smoothing increased. Whereas the range of unconstrained angiosperm age estimates spans the fossil age of the clade, the range of constrained estimates is narrower (and older) than the earliest angiosperm fossils. Results unambiguously indicate the relevance of constraints in reducing the variability of ages derived from different partitions of the data and diminishing the effect of the smoothing parameter. Constrained optimizations of divergence times and substitution rates across the phylogeny suggest appreciably different evolutionary dynamics for angiosperms and for gymnosperms. Whereas the gymnosperm crown group originated shortly after the origin of seed plants, a long time elapsed before the origin of crown group angiosperms. Although absolute age estimates of angiosperms and angiosperm clades are older than their earliest fossils, the estimated pace of phylogenetic diversification largely agrees with the rapid appearance of angiosperm lineages in stratigraphic sequences.

Using models of nucleotide evolution to build phylogenetic trees

Molecular phylogenetics and its applications are popular and useful tools for making comparative investigations in genetics; however, estimating phylogenetic trees is not always straightforward. Some phylogenetic estimators use an explicit model of nucleotide evolution to estimate evolutionary parameters such as branch lengths and tree topology. There are many models to choose from, and use of the optimal model for a particular data set is important to avoid a loss of power and accuracy in phylogenetic estimations. Here, we review some molecular evolutionary forces and the parameters included in some common models of evolution used to interpret resulting patterns of molecular variation. We present some statistical methods of selecting a particular model of nucleotide evolution, and provide an empirical example of model selection. Statistical model selection strikes a balance between the bias introduced by some models and the increased variance of parameter estimates that results from using other models.

Assessing the quality of molecular divergence time estimates by fossil calibrations and fossil-based model selection

Estimates of species divergence times using DNA sequence data are playing an increasingly important role in studies of evolution, ecology and biogeography. Most work has centred on obtaining appropriate kinds of data and developing optimal estimation procedures, whereas somewhat less attention has focused on the calibration of divergences using fossils. Case studies with multiple fossil calibration points provide important opportunities to examine the divergence time estimation problem in new ways. We discuss two cross-validation procedures that address different aspects of inference in divergence time estimation. 'Fossil cross-validation' is a procedure used to identify the impact of different individual calibrations on overall estimation. This can identify fossils that have an exceptionally large error effect and may warrant further scrutiny. 'Fossil-based model cross-validation' is an entirely different procedure that uses fossils to identify the optimal model of molecular evolution in the context of rate smoothing or other inference methods. Both procedures were applied to two recent studies: an analysis of monocot angiosperms with eight fossil calibrations and an analysis of placental mammals with nine fossil calibrations. In each case, fossil calibrations could be ranked from most to least influential, and in one of the two studies, the fossils provided decisive evidence about the optimal molecular evolutionary model.

Evolution of the species-rich Cape flora

The Cape Floristic Region ('fynbos biome') has very high levels of plant species diversity and endemism. Much of this diversity is concentrated in a relatively small number of clades centered in the region (Cape clades), and these form a vegetation called 'fynbos'. The general explanation for the origin of this diversity is that much of it evolved in the Pliocene and Late Miocene in response to progressive aridification. We present a phylogenetic analysis of an almost complete species sample of the largest clade of Restionaceae, the third largest Cape clade. This indicates that the radiation of the Restionaceae started between 20 and 42 Myr ago, and since then there were no, or at most gradual, changes in the speciation rate in this clade. For seven other clades, the estimated starting dates for their radiation ranges from 7 to 20 Myr ago. Combining the radiation patterns for these clades shows that ca. 15% of the modern species evolved during the Pleistocene, and almost 40% since the beginning of the Pliocene. We suggest that these clades might have radiated in response to the fynbos vegetation increasing its extent in the Cape as a result of climatic change.

Genomic data support the hominoid slowdown and an Early Oligocene estimate for the hominoid-cercopithecoid divergence

Several lines of indirect evidence suggest that hominoids (apes and humans) and cercopithecoids (Old World monkeys) diverged around 23-25 Mya. Importantly, although this range of dates has been used as both an initial assumption and as a confirmation of results in many molecular-clock analyses, it has not been critically assessed on its own merits. In this article we test the robusticity of the 23- to 25-Mya estimate with approximately 150,000 base pairs of orthologous DNA sequence data from two cercopithecoids and two hominoids by using quartet analysis. This method is an improvement over other estimates of the hominoid-cercopithecoid divergence because it incorporates two calibration points, one each within cercopithecoids and hominoids, and tests for a statistically appropriate model of molecular evolution. Most comparisons reject rate constancy in favor of a model incorporating two rates of evolution, supporting the "hominoid slowdown" hypothesis. By using this model of molecular evolution, the hominoid-cercopithecoid divergence is estimated to range from 29.2 to 34.5 Mya, significantly older than most previous analyses. Hominoid-cercopithecoid divergence dates of 23-25 Mya fall outside of the confidence intervals estimated, suggesting that as much as one-third of ape evolution has not been paleontologically sampled. Identifying stem cercopithecoids or hominoids from this period will be difficult because derived features that define crown catarrhines need not be present in early members of these lineages. More sites that sample primate habitats from the Oligocene of Africa are needed to better understand early ape and Old World monkey evolution.

Evidence for multiple lateral transfers of the circadian clock cluster in filamentous heterocystic cyanobacteria Nostocaceae

Cyanobacteria are the first prokaryotes reported to show circadian rhythmicity, which is regulated by a cluster of three genes: kaiA, kaiB, and kaiC. Phylogenetic analysis of the kaiBC cluster in filamentous cyanobacteria of the family Nostocaceae including Nodularia spumigena and Nostoc linckia from Arubotaim Cave, Mt. Sedom, Israel, indicated that this cluster has experienced multiple lateral transfers. The transfers have occurred in different periods of the species' evolution. The data obtained suggest that lateral transfers of the circadian clock cluster in filamentous cyanobacteria have been common and might have adaptive significance.

Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context

The lemurs of Madagascar are a unique radiation of primates that show an extraordinary diversity of lifestyles, morphologies and behaviours. However, very little is known about the relative antiquity of lemuriform clades due to the lack of terrestrial fossils for the Tertiary of Madagascar. Here, we employ a Bayesian method to estimate divergence dates within the lemuriform radiation using several unlinked gene loci and multiple fossil calibrations outside the lemuriform clade. Two mitochondrial genes (cytochrome oxidase II and cytochrome b), two nuclear introns (transthyretin intron 1 and von Willebrand factor gene intron 11) and one nuclear exon (interphotoreceptor retinoid binding protein, exon 1) are used in separate and combined analyses. The genes differ in taxon sampling and evolutionary characteristics but produce congruent date estimates. Credibility intervals narrow considerably in combined analyses relative to separate analyses due to the increased amount of data. We also test the relative effects of multiple vs. single calibration points, finding that, when only single calibration points are employed, divergence dates are systematically underestimated. For the mitochondrial DNA data set, we investigate the effects of sampling density within the mouse lemur radiation (genus Microcebus). When only two representative species are included, estimated dates throughout the phylogeny are more recent than with the complete-species sample, with basal nodes less affected than recent nodes. The difference appears to be due to the manner in which priors on node ages are constructed in the two analyses. In nearly all analyses, the age of the lemuriform clade is estimated to be approximately 62-65 Ma, with initial radiation of mouse lemurs and true lemurs (genus Eulemur) occurring approximately 8-12 Ma. The antiquity of the mouse lemur radiation is surprising given the near uniform morphology among species. Moreover, the observation that mouse lemurs and true lemurs are of similar ages suggests discrepancies in rates of morphological, behavioural and physiological evolution in the two clades, particularly with regard to characteristics of sexual signalling. These differences appear to correlate with the nocturnal vs. diurnal lifestyles, respectively, of these two primate groups.

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.

Molecular evidence on plant divergence times

Estimation of divergence times from sequence data has become increasingly feasible in recent years. Conflicts between fossil evidence and molecular dates have sparked the development of new methods for inferring divergence times, further encouraging these efforts. In this paper, available methods for estimating divergence times are reviewed, especially those geared toward handling the widespread variation in rates of molecular evolution observed among lineages. The assumptions, strengths, and weaknesses of local clock, Bayesian, and rate smoothing methods are described. The rapidly growing literature applying these methods to key divergence times in plant evolutionary history is also reviewed. These include the crown group ages of green plants, land plants, seed plants, angiosperms, and major subclades of angiosperms. Finally, attempts to infer divergence times are described in the context of two very different temporal settings: recent adaptive radiations and much more ancient biogeographic patterns.

Origin and evolution of the light-dependent protochlorophyllide oxidoreductase (LPOR) genes

Light-dependent NADPH-protochlorophyllide oxidoreductase (LPOR) is a nuclear-encoded chloroplast protein in green algae and higher plants which catalyzes the light-dependent reduction of protochlorophyllide to chlorophyllide. Light-dependent chlorophyll biosynthesis occurs in all oxygenic photosynthetic organisms. With the exception of angiosperms, this pathway coexists with a separate light-independent chlorophyll biosynthetic pathway, which is catalyzed by light-independent protochlorophyllide reductase (DPOR) in the dark. In contrast, the light-dependent function of chlorophyll biosynthesis is absent from anoxygenic photosynthetic bacteria. Consequently, the question is whether cyanobacteria are the ancestors of all organisms that conduct light-dependent chlorophyll biosynthesis. If so, how did photosynthetic eukaryotes acquire the homologous genes of LPOR in their nuclear genomes? The large number of complete genome sequences now available allow us to detect the evolutionary history of LPOR genes by conducting a genome-wide sequence comparison and phylogenetic analysis. Here, we show the results of a detailed phylogenetic analysis of LPOR and other functionally related enzymes in the short chain dehydrogenase/reductase (SDR) family. We propose that the LPOR gene originated in the cyanobacterial genome before the divergence of eukaryotic photosynthetic organisms. We postulated that the photosynthetic eukaryotes obtained their LPOR homologues through endosymbiotic gene transfer.

Phylogenomics of Eukaryotes: Impact of Missing Data on Large Alignments

Resolving the relationships between Metazoa and other eukaryotic groups as well as between metazoan phyla is central to the understanding of the origin and evolution of animals. The current view is based on limited data sets, either a single gene with many species (e.g., ribosomal RNA) or many genes but with only a few species. Because a reliable phylogenetic inference simultaneously requires numerous genes and numerous species, we assembled a very large data set containing 129 orthologous proteins ( approximately 30,000 aligned amino acid positions) for 36 eukaryotic species. Included in the alignments are data from the choanoflagellate Monosiga ovata, obtained through the sequencing of about 1,000 cDNAs. We provide conclusive support for choanoflagellates as the closest relative of animals and for fungi as the second closest. The monophyly of Plantae and chromalveolates was recovered but without strong statistical support. Within animals, in contrast to the monophyly of Coelomata observed in several recent large-scale analyses, we recovered a paraphyletic Coelamata, with nematodes and platyhelminths nested within. To include a diverse sample of organisms, data from EST projects were used for several species, resulting in a large amount of missing data in our alignment (about 25%). By using different approaches, we verify that the inferred phylogeny is not sensitive to these missing data. Therefore, this large data set provides a reliable phylogenetic framework for studying eukaryotic and animal evolution and will be easily extendable when large amounts of sequence information become available from a broader taxonomic range.

Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Reconstructing the chronology of mammalian evolution is a debated issue between molecule- and fossil-based inferences. A methodological limitation of molecules is the evolutionary rate variation among lineages, precluding the application of the global molecular clock. We considered 2422 first and second codon positions of the combined ADRA2B, IRBP, and vWF nuclear genes for a well-documented set of placentals including an extensive sampling of rodents. Using seven independent calibration points and a maximum-likelihood framework, we evaluated whether molecular and paleontological estimates of mammalian divergence dates may be reconciled by the local molecular clocks approach, allowing local constancy of substitution rates with variations at larger phylogenetic scales. To handle the difficulty of choosing among all possible rate assignments for various lineages, local molecular clocks were based on the results of branch-length and two-cluster tests. Extensive lineage-specific variation of evolutionary rates was detected, even among rodents. Cross-calibrations indicated some incompatibilities between divergence dates based on different paleontological references. To decrease the impact of a single calibration point, estimates derived from independent calibrations displaying only slight reciprocal incompatibility were averaged. The divergence dates inferred for the split between mice and rats (approximately 13-19 Myr) was younger than previously published molecular estimates. The most recent common ancestors of rodents, primates and rodents, boreoeutherians, and placentals were estimated to be, respectively, approximately 60, 70, 75, and 78 Myr old. Global clocks, local clocks, and quartet dating analyses suggested a Late Cretaceous origin of the crown placental clades followed by a Tertiary radiation of some placental orders like rodents.

Unraveling the evolutionary radiation of the thoracican barnacles using molecular and morphological evidence: a comparison of several divergence time estimation approaches

The Thoracica includes the ordinary barnacles found along the sea shore and is the most diverse and well-studied superorder of Cirripedia. However, although the literature abounds with scenarios explaining the evolution of these barnacles, very few studies have attempted to test these hypotheses in a phylogenetic context. The few attempts at phylogenetic analyses have suffered from a lack of phylogenetic signal and small numbers of taxa. We collected DNA sequences from the nuclear 18S, 28S, and histone H3 genes and the mitochondrial 12S and 16S genes (4,871 bp total) and data for 37 adult and 53 larval morphological characters from 43 taxa representing all the extant thoracican suborders (except the monospecific Brachylepadomorpha). Four Rhizocephala (highly modified parasitic barnacles) taxa and a Rhizocephala + Acrothoracica (burrowing barnacles) hypothetical ancestor were used as the outgroup for the molecular and morphological analyses, respectively. We analyzed these data separately and combined using maximum likelihood (ML) under "hill-climbing" and genetic algorithm heuristic searches, maximum parsimony procedures, and Bayesian inference coupled with Markov chain Monte Carlo techniques under mixed and homogeneous models of nucleotide substitution. The resulting phylogenetic trees answered key questions in barnacle evolution. The four-plated Iblomorpha were shown as the most primitive thoracican, and the plateless Heteralepadomorpha were placed as the sister group of the Lepadomorpha. These relationships suggest for the first time in an invertebrate that exoskeleton biomineralization may have evolved from phosphatic to calcitic. Sessilia (nonpedunculate) barnacles were depicted as monophyletic and appear to have evolved from a stalked (pedunculate) multiplated (5+) scalpelloidlike ancestor rather than a five-plated lepadomorphan ancestor. The Balanomorpha (symmetric sessile barnacles) appear to have the following relationship: (Chthamaloidea(Coronuloidea(Tetraclitoidea, Balanoidea))). Thoracican divergence times were estimated under ML-based local clock, Bayesian, and penalized likelihood approaches using an 18S data set and three calibration points: Heteralepadomorpha = 530 million years ago (MYA), Scalpellomorpha = 340 MYA, and Verrucomorpha = 120 MYA. Estimated dates varied considerably within and between approaches depending on the calibration point. Highly parameterized local clock models that assume independent rates (r > or = 15) for confamilial or congeneric species generated the most congruent estimates among calibrations and agreed more closely with the barnacle fossil record. Reasonable estimates were also obtained under the Bayesian procedure of Kishino et al. (2001, Mol. Biol. Evol. 18:352-361) but using multiple calibrations. Most of the dates estimated under the Bayesian procedure of Aris-Brosou and Yang (2002, Syst. Biol. 51:703-714) and the penalized likelihood method using single and/or multiple calibrations were inconsistent among calibrations and did not fit the fossil record.

Evolution of the fish rhabdovirus viral haemorrhagic septicaemia virus

Viral haemorrhagic septicaemia (VHS) caused by the rhabdovirus VHSV is economically the most important viral disease in European rainbow trout farming. Until 1989, this virus was mainly isolated from freshwater salmonids but in the last decade, it has also been isolated from an increasing number of free-living marine fish species. To study the genetic evolution of VHSV, the entire G gene from 74 isolates was analysed. VHSV from wild marine species caught in the Baltic Sea, Skagerrak, Kattegat, North Sea, and English Channel and European freshwater isolates, appeared to share a recent common ancestor. Based on the estimated nucleotide substitution rate, the ancestor of the European fresh water isolates was dated some 50 years ago. This finding fits with the initial reports in the 1950s on clinical observations of VHS in Danish freshwater rainbow trout farms. The study also indicates that European marine VHSV and the North American marine line separated approx. 500 years ago. The codon substitution rate among the freshwater VHSV isolates was found to be 2.5 times faster than among marine isolates. The data support the hypothesis of the marine environment being the original reservoir of VHSV and that the change in host range (to include rainbow trout) may have occurred several times. Virus from the marine environment will therefore continue to represent a threat to the trout aquaculture industry.

African Haplogroup L mtDNA sequences show violations of clock-like evolution

A set of 96 complete mtDNA sequences that belong to the three major African haplogroups (L1, L2, and L3) was analyzed to determine if mtDNA has evolved as a molecular clock. Likelihood ratio tests (LRTs) were carried out with each of the haplogroups and with combined haplogroup sequence sets. Evolution has not been clock-like, neither for the coding region nor for the control region, in combined sets of African haplogroup L mtDNA sequences. In tests of individual haplogroups, L2 mtDNAs showed violations of a molecular clock under all conditions and in both the control and coding regions. In contrast, haplogroup L1 and L3 sequences, both for the coding and control regions, show clock-like evolution. In clock tests of individual L2 subclades, the L2a sequences showed a marked violation of clock-like evolution within the coding region. In addition, the L2a and L2c branch lengths of both the coding and control regions were shorter relative to those of the L2b and L2d sequences, a result that indicates lower levels of sequence divergence. Reduced median network analyses of the L2a sequences indicated the occurrence of marked homoplasy at multiple sites in the control region. After exclusion of the L2a and L2c sequences, African mtDNA coding region evolution has not significantly departed from a molecular clock, despite the results of neutrality tests that indicate the mitochondrial coding region has evolved under nonneutral conditions. In contrast, control region evolution is clock-like only at the haplogroup level, and it thus appears to have evolved essentially independently from the coding region. The results of the clock tests, the network analyses, and the branch length comparisons all caution against the use of simple mtDNA clocks.

CALIBRATION AGE AND QUARTET DIVERGENCE DATE ESTIMATION

The date of a single divergence point--between living alligators and crocodiles--was estimated with quartet dating using calibrations of widely divergent ages. For five mitochondrial sequence datasets, there is a clear relationship between calibration age and quartet estimate--quartets based on two relatively recent calibrations support younger divergence estimates than do quartets based on two older calibrations. Some of the estimates supported by young quartets are impossibly young and exclude the first appearance of the group in the fossil record as too old. The older estimates--those based on two relatively old calibrations--may be overestimates, and those based on one old and one recent calibration support divergence estimates very close to fossil data. This suggests that quartet dating methods may be most effective when calibrations are applied from different parts of a clade's history.

Molecular Phylogenetic Dating of Asterid Flowering Plants Shows Early Cretaceous Diversification

We present a phylogenetic dating of asterids, based on a 111-taxon tree representing all major groups and orders and 83 of the 102 families of asterids, with an underlying data set comprising six chloroplast DNA markers totaling 9914 positions. Phylogenetic dating was done with semiparametric rate smoothing by penalized likelihood. Confidence intervals were calculated by bootstrapping. Six reference fossils were used for calibration. To explore the effects of various sources of error, we repeated the analyses with alternative dating methods (nonparametric rate smoothing and the Langley-Fitch clock-based method), alternative tree topologies, reduced taxon sampling (22 of the 111 taxa deleted), partitioning the data into three genes and three noncoding regions, and calibrating with single reference fossils. The analyses with alternative topologies, reduced taxon sampling, and coding versus noncoding sequences all yielded small or in some cases no deviations. The choice of method influenced the age estimates of a few nodes considerably. Calibration with reference fossils is a critical issue, and use of single reference fossils yielded different results depending on the fossil. The bootstrap confidence intervals were generally small. Our results show that asterids and their major subgroups euasterids, campanulids, and lamiids diversified during the Early Cretaceous. Cornales, Ericales, and Aquifoliales also have crown node ages from the Early Cretaceous. Dipsacales and Solanales are from the Mid-Cretaceous, the other orders of core campanulids and core lamiids from the Late Cretaceous. The considerable diversity exhibited by asterids almost from their first appearance in the fossil record also supports an origin and first phase of diversification in the Early Cretaceous.

Influence of Tertiary paleoenvironmental changes on the diversification of South American mammals: a relaxed molecular clock study within xenarthrans

BACKGROUND

Comparative genomic data among organisms allow the reconstruction of their phylogenies and evolutionary time scales. Molecular timings have been recently used to suggest that environmental global change have shaped the evolutionary history of diverse terrestrial organisms. Living xenarthrans (armadillos, anteaters and sloths) constitute an ideal model for studying the influence of past environmental changes on species diversification. Indeed, extant xenarthran species are relicts from an evolutionary radiation enhanced by their isolation in South America during the Tertiary era, a period for which major climate variations and tectonic events are relatively well documented.

RESULTS

We applied a Bayesian approach to three nuclear genes in order to relax the molecular clock assumption while accounting for differences in evolutionary dynamics among genes and incorporating paleontological uncertainties. We obtained a molecular time scale for the evolution of extant xenarthrans and other placental mammals. Divergence time estimates provide substantial evidence for contemporaneous diversification events among independent xenarthran lineages. This correlated pattern of diversification might possibly relate to major environmental changes that occurred in South America during the Cenozoic.

CONCLUSIONS

The observed synchronicity between planetary and biological events suggests that global change played a crucial role in shaping the evolutionary history of extant xenarthrans. Our findings open ways to test this hypothesis further in other South American mammalian endemics like hystricognath rodents, platyrrhine primates, and didelphid marsupials.

Molecular cloning of the estrogen and progesterone receptors of the American alligator

Steroid hormones perform many essential roles in vertebrates during embryonic development, reproduction, growth, water balance, and responses to stress. The estrogens are essential for normal reproductive activity in female and male vertebrates and appear to have direct actions during sex determination in some vertebrates. To begin to understand the molecular mechanisms of estrogen action in alligators, we have isolated cDNAs encoding the estrogen receptors (ER) from the ovary. Degenerate PCR primers specific to ER were designed and used to amplify alligator ovary RNA. Two different DNA fragments (ERalpha and ERbeta) were obtained and the full-length alligator ERalpha cDNA was obtained using 5' and 3' RACE. The inferred amino acid sequence of alligator ERalpha (aERalpha) was very similar to the chicken ERalpha (91% identity), although phylogenetic analyses suggested profound differences in the rate of sequence evolution for vertebrate ER sequences. We also isolated partial DNA fragments encoding ERbeta and the progesterone receptor (PR) of the alligator, both of which show strong sequence similarities to avian ERbeta and PR. We examined the expression levels of these three steroid receptors (ERalpha, ERbeta, and PR) in the ovary of juvenile alligators and observed detectable levels of all three receptors. Quantitative RT-PCR showed that gonadal ERalpha transcript levels in juvenile alligators decreased after E2 treatment whereas ERbeta and PR transcripts were not changed. These results provide tools that will allow future studies examining the regulation and ontogenic expression of steroid receptors in alligators and expand our knowledge of vertebrate steroid receptor evolution.

Site stripping based on likelihood ratio reduction is a useful tool to evaluate the impact of non-clock-like behavior on viral phylogenetic reconstructions

The site stripping for clock detection procedure was implemented in the recently developed maximum likelihood framework for estimating evolutionary rates and divergence times in measurably evolving populations. The method was used to investigate the effect of rate variability on estimating divergence times in non-clock-like trees for human immunodeficiency viruses and hepatitis C viruses. We validate our approach by comparing dated coalescent nodes in molecular phylogenies with known dates of transmission. Our method was able to rapidly recover clock-like behavior and to indicate the presence and direction of a bias when estimates of divergence times using the unstripped data were flawed.

A decade of progress in plant molecular phylogenetics

Over the past decade, botanists have produced several thousand phylogenetic analyses based on molecular data, with particular emphasis on sequencing rbcL, the plastid gene encoding the large subunit of Rubisco (ribulose bisphosphate carboxylase). Because phylogenetic trees retrieved from the three plant genomes (plastid, nuclear and mitochondrial) have been highly congruent, the "Angiosperm Phylogeny Group" has used these DNA-based phylogenetic trees to reclassify all families of flowering plants. However, in addition to taxonomy, these major phylogenetic efforts have also helped to define strategies to reconstruct the "tree of life", and have revealed the size of the ancestral plant genome, uncovered potential candidates for the ancestral flower, identified molecular living fossils, and linked the rate of neutral substitutions with species diversity. With an increased interest in DNA sequencing programmes in non-model organisms, the next decade will hopefully see these phylogenetic findings integrated into new genetic syntheses, from genomes to taxa.

Molecular systematics of dormice (Rodentia: Gliridae) and the radiation of Graphiurus in Africa

The phylogenetic relationships among the Gliridae (order Rodentia) were assessed using 3430 nucleotides derived from three nuclear fragments (beta-spectrin non-erythrocytic 1, thyrotropin and lecithin cholesterol acyl transferase) and one mitochondrial gene (12S rRNA). We included 14 glirid species, representative of seven genera of the three recognized subfamilies (Graphiurinae, Glirinae and Leithiinae) in our analysis. The molecular data identified three evolutionary lineages that broadly correspond to the three extant subfamilies. However, the data suggest that the genus Muscardinus, previously regarded as falling within the Glirinae, should be included in the Leithiinae. Molecular dating using local molecular clocks and partitioned datasets allowed an estimate of the timing of cladogenesis within the glirids. Graphiurus probably diverged early in the group's evolution (40-50 Myr ago) and the three subfamilies diverged contemporaneously, probably in Europe. The radiation within Graphiurus is more recent, with the colonization of Africa by this lineage estimated at ca. 8-10 Myr ago.

Molecular estimation of eulipotyphlan divergence times and the evolution of "Insectivora"

"Insectivores" are one of the key groups in understanding mammalian origins. For years, systematics of "Lipotyphla" taxa remained extremely unstable and challenged. Today, with the application of molecular techniques, "Lipotyphla" appears to be a paraphyletic assemblage that encompasses hedgehogs, shrews, and moles (i.e., Eulipotyphla-a member of Laurasiatheria), and golden moles and tenrecs (i.e., Afrosoricida-a member of Afrotheria). Based on nuclear genes and on this well-established phylogenetic framework, we estimated Bayesian relaxed molecular clock divergence times among major lineages of "Lipotyphla." Crown placental mammals are shown to diversify 102+/-6 million years ago (Mya; mean+/-one standard-deviation), followed by Boreoeutheria (94+/-6 Mya), Laurasiatheria (85+/-5 Mya), and Eulipotyphla (73+/-5), with moles separating from hedgehogs+shrews just at the K/T boundary (65+/-5 Mya). During the Early and Middle Eocene, all extant eulipotyphlan subfamilies originated: Uropsilinae (52+/-5 Mya), and Desmaninae, Talpinae, Erinaceinae, Hylomyinae, Soricinae, and Crocidurinae (38-42+/-5 Mya). Afrosoricida separated from Macroscelidae 69+/-5 Mya, golden moles from tenrecs 63+/-5 Mya, and the diversification within tenrecs occurred 43+/-5 Mya. Divergence times are shown to be in reasonably good agreement with the fossil record of eulipotyphlans, but not with the one of afrosoricid "insectivores." Eulipotyphlans diversification might have been sculpted by variations in paleoclimates of the cenozoic era.

A neutral explanation for the correlation of diversity with recombination rates in humans

One of the most striking findings to emerge from the study of genomic patterns of variation is that regions with lower recombination rates tend to have lower levels of intraspecific diversity but not of interspecies divergence. This uncoupling of variation within and between species has been widely interpreted as evidence that natural selection shapes patterns of genetic variability genomewide. We revisited the relationship between diversity, divergence, and recombination in humans, using data from closely related species and better estimates of recombination rates than previously available. We show that regions that experience less recombination have reduced divergence to chimpanzee and to baboon, as well as lower levels of diversity. This observation suggests that mutation and recombination are associated processes in humans, so that the positive correlation between diversity and recombination may have a purely neutral explanation. Consistent with this hypothesis, diversity levels no longer increase significantly with recombination rates after correction for divergence to chimpanzee.

Origin and evolution of circadian clock genes in prokaryotes

Regulation of physiological functions with approximate daily periodicity, or circadian rhythms, is a characteristic feature of eukaryotes. Until recently, cyanobacteria were the only prokaryotes reported to possess circadian rhythmicity. It is controlled by a cluster of three genes: kaiA, kaiB, and kaiC. Using sequence data of approximately 70 complete prokaryotic genomes from the various public depositories, we show here that the kai genes and their homologs have quite a different evolutionary history and occur in Archaea and Proteobacteria as well. Among the three genes, kaiC is evolutionarily the oldest, and kaiA is the youngest and likely evolved only in cyanobacteria. Our data suggest that the prokaryotic circadian pacemakers have evolved in parallel with the geological history of the earth, and that natural selection, multiple lateral transfers, and gene duplications and losses have been the major factors shaping their evolution.

Higher-level systematics of rodents and divergence time estimates based on two congruent nuclear genes

Phylogenetic analysis of over 4600 aligned nucleotide sequences from two nuclear genes, growth hormone receptor and BRCA1, provided congruent phylogenies depicting relationships among the major lineages of rodents. Separate and combined analyses resulted in five major conclusions: (1) strong support for a monophyletic Myodonta (containing the superfamilies Muroidea + Dipodoidea), with subfamily Gerbillinae being more closely related to Murinae than is Sigmodontinae; (2) a sister-group relationship between the family Castoridae and the superfamily Geomyoidea; (3) monophyly of Ctenohystrica (containing the suborders Sciuravida and Hystricognatha); (4) a near polytomy among Myodonta (suborder Myomorpha), Pedetes (family Pedetidae, suborder Anomaluromorpha), Castoridae (suborder Sciuromorpha) + Geomyoidea (suborder Myomorpha), and Ctenohystrica; and (5) basal position of a monophyletic group containing Graphiurus (family Gliridae, suborder Myomorpha) + two members of the Sciuromorpha (Sciuridae + Aplodontidae). Divergence dates among rodents and primates were also estimated using the combined data. Applying a global molecular clock and a primate calibration point, divergence dates among rodents exceeded fossil-based dates but were generally compatible with other molecule-based dates estimated under similar conditions. However, when a relaxed molecular clock was applied, estimated divergence dates were highly compatible with the fossil record.

Mitochondrial phylogeography of the Woodmouse (Apodemus sylvaticus) in the Western Palearctic region

We sequenced 965 base pairs of the mitochondrial DNA cytochrome b from 102 woodmice (Apodemus sylvaticus) collected from 40 European localities. The aims of the study were to answer the following questions. (i) Did the Mediterranean peninsulas play a role as refuge for woodmice? (ii) Is genetic variability of A. sylvaticus higher in the Mediterranean region compared with northern Europe? (iii) Are the patterns of the postglacial colonization of Europe by woodmice similar to those presently recognized for other European species? The results provide a clear picture of the impact of the Quaternary glaciations on the genetic and geographical structure of the woodmouse. Our analyses indicate a higher genetic variability of woodmice in the Mediterranean peninsulas compared to northern Europe, suggesting a role of the former as refuge regions for this small mammal. An original pattern of postglacial colonization is proposed where the Iberian and southern France refuge populations colonized almost all European regions. The Sicilian population appears to be very differentiated and highly variable. This emphasizes the importance of this island as a 'hot spot' for the intraspecific genetic diversity of the woodmouse. Finally, woodmice in North Africa originated from southwestern Europe, most probably as a result of a recent anthropogenic introduction.

Time scale of eutherian evolution estimated without assuming a constant rate of molecular evolution

Controversies over the molecular clock hypothesis were reviewed. Since it is evident that the molecular clock does not hold in an exact sense, accounting for evolution of the rate of molecular evolution is a prerequisite when estimating divergence times with molecular sequences. Recently proposed statistical methods that account for this rate variation are overviewed and one of these procedures is applied to the mitochondrial protein sequences and to the nuclear gene sequences from many mammalian species in order to estimate the time scale of eutherian evolution. This Bayesian method not only takes account of the variation of molecular evolutionary rate among lineages and among genes, but it also incorporates fossil evidence via constraints on node times. With denser taxonomic sampling and a more realistic model of molecular evolution, this Bayesian approach is expected to increase the accuracy of divergence time estimates.

Molecular systematics and biogeography of the Neotropical monkey genus, Alouatta

We take advantage of the broad distribution of howler monkeys from Mexico to Argentina to provide a historical biogeographical analysis on a regional scale that encompasses the entire Neotropics. The phylogenetic relationships among 9 of the 10 recognized Alouatta species were inferred using three mitochondrial and two nuclear genes. The nuclear gene regions provided no phylogenetic resolution among howler monkey species, and were characterized by very low levels of sequence divergence between Alouatta and the Ateles outgroup. The mtDNA genes, on the other hand, produced a well-resolved phylogeny, which indicated that the earliest split among howler monkeys separated cis- and trans-Andean clades. Eight monophyletic mtDNA haplotype clades were identified, representing six named species in South America, including Alouatta seniculus, Alouatta sara, Alouatta macconelli, Alouatta caraya, Alouatta belzebul, and Alouatta guariba, and two in Mesoamerica, Alouatta pigra and Alouatta palliata. Molecular clock-based estimates of branching times indicated that contemporary howler monkey species originated in the late Miocene and Pliocene, not the Pleistocene. The causes of Alouatta diversification were more difficult to pin down, although we posit that the initial cis-, trans-Andean split in the genus was caused by the late Miocene completion of the northern Andes. Riverine barriers to dispersal and putative forest refuges can neither be discounted nor distinguished as causes of speciation in many cases, and one, the other or both have likely played a role in the diversification of South American howler monkeys. Finally, we estimated the separation of Mesoamerican A. pigra and A. palliata at 3Ma, which corresponds to the completion date of the Panama Isthmus promoting a role for this earth history event in the speciation of Central American howler monkeys.

Likelihood and Bayes estimation of ancestral population sizes in hominoids using data from multiple loci

Polymorphisms in an ancestral population can cause conflicts between gene trees and the species tree. Such conflicts can be used to estimate ancestral population sizes when data from multiple loci are available. In this article I extend previous work for estimating ancestral population sizes to analyze sequence data from three species under a finite-site nucleotide substitution model. Both maximum-likelihood (ML) and Bayes methods are implemented for joint estimation of the two speciation dates and the two population size parameters. Both methods account for uncertainties in the gene tree due to few informative sites at each locus and make an efficient use of information in the data. The Bayes algorithm using Markov chain Monte Carlo (MCMC) enjoys a computational advantage over ML and also provides a framework for incorporating prior information about the parameters. The methods are applied to a data set of 53 nuclear noncoding contigs from human, chimpanzee, and gorilla published by Chen and Li. Estimates of the effective population size for the common ancestor of humans and chimpanzees by both ML and Bayes methods are approximately 12,000-21,000, comparable to estimates for modern humans, and do not support the notion of a dramatic size reduction in early human populations. Estimates published previously from the same data are several times larger and appear to be biased due to methodological deficiency. The divergence between humans and chimpanzees is dated at approximately 5.2 million years ago and the gorilla divergence 1.1-1.7 million years earlier. The analysis suggests that typical data sets contain useful information about the ancestral population sizes and that it is advantageous to analyze data of several species simultaneously.

A novel use of equilibrium frequencies in models of sequence evolution

Current mathematical models of amino acid sequence evolution are often applied in variants that match their expected amino acid frequencies to those observed in a data set under analysis. This has been achieved by setting the instantaneous rate of replacement of a residue i by another residue j proportional to the observed frequency of the resulting residue j. We describe a more general method that maintains the match between expected and observed frequencies but permits replacement rates to be proportional to the frequencies of both the replaced and resulting residues, raised to powers other than 1. Analysis of a database of amino acid alignments shows that the description of the evolutionary process in a majority (approximately 70% of 182 alignments) is significantly improved by use of the new method, and a variety of analyses indicate that parameter estimation with the new method is well-behaved. Improved evolutionary models increase our understanding of the process of molecular evolution and are often expected to lead to improved phylogenetic inferences, and so it seems justified to consider our new variants of existing standard models when performing evolutionary analyses of amino acid sequences. Similar methods can be used with nucleotide substitution models, but we have not found these to give corresponding significant improvements to our ability to describe the processes of nucleotide sequence evolution.