Phylogenetic relationships and divergence times among dormice (Rodentia, Gliridae) based on three nuclear genes

We examined phylogenetic relationships among six species representing three subfamilies, Glirinae, Graphiurinae and Leithiinae with sequences from three nuclear protein-coding genes (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). Phylogenetic trees reconstructed from maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian-inference (BI) analyses showed the monophyly of Glirinae (Glis and Glirulus) and Leithiinae (Dryomys, Eliomys and Muscardinus) with strong support, although the branch length maintaining this relationship was very short, implying rapid diversification among the three subfamilies. Divergence time estimates were calculated from ML (local clock model) and Bayesian-dating method using a calibration point of 25 Myr (million years) ago for the divergence between Glis and Glirulus, and 55 Myr ago for the split between lineages of Gliridae and Sciuridae on the basis of fossil records. The results showed that each lineage of Graphiuros, Glis, Glirulus and Muscardinus dates from the Late Oligocene to the Early Miocene period, which is mostly in agreement with fossil records. Taking into account that warm climate harbouring a glirid-favoured forest dominated from Europe to Asia during this period, it is considered that this warm environment triggered the prosperity of the glirid species through the rapid diversification. Glirulus japonicas is suggested to be a relict of this ancient diversification during the warm period.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Determining divergence times with a protein clock: Update and reevaluation

A recent study of the divergence times of the major groups of organisms as gauged by amino acid sequence comparison has been expanded and the data have been reanalyzed with a distance measure that corrects for both constraints on amino acid interchange and variation in substitution rate at different sites. Beyond that, the availability of complete genome sequences for several eubacteria and an archaebacterium has had a great impact on the interpretation of certain aspects of the data. Thus, the majority of the archaebacterial sequences are not consistent with currently accepted views of the Tree of Life which cluster the archaebacteria with eukaryotes. Instead, they are either outliers or mixed in with eubacterial orthologs. The simplest resolution of the problem is to postulate that many of these sequences were carried into eukaryotes by early eubacterial endosymbionts about 2 billion years ago, only very shortly after or even coincident with the divergence of eukaryotes and archaebacteria. The strong resemblances of these same enzymes among the major eubacterial groups suggest that the cyanobacteria and Gram-positive and Gram-negative eubacteria also diverged at about this same time, whereas the much greater differences between archaebacterial and eubacterial sequences indicate these two groups may have diverged between 3 and 4 billion years ago.

Using rare mutations to estimate population divergence times: A maximum likelihood approach

In this paper we propose a method to estimate by maximum likelihood the divergence time between two populations, specifically designed for the analysis of nonrecurrent rare mutations. Given the rapidly growing amount of data, rare disease mutations affecting humans seem the most suitable candidates for this method. The estimator RD, and its conditional version RDc, were derived, assuming that the population dynamics of rare alleles can be described by using a birth–death process approximation and that each mutation arose before the split of a common ancestral population into the two diverging populations. The RD estimator seems more suitable for large sample sizes and few alleles, whose age can be approximated, whereas the RDc estimator appears preferable when this is not the case. When applied to three cystic fibrosis mutations, the estimator RD could not exclude a very recent time of divergence among three Mediterranean populations. On the other hand, the divergence time between these populations and the Danish population was estimated to be, on the average, 4,500 or 15,000 years, assuming or not a selective advantage for cystic fibrosis carriers, respectively. Confidence intervals are large, however, and can probably be reduced only by analyzing more alleles or loci.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Phylogenetics of Olea (Oleaceae) based on plastid and nuclear ribosomal DNA sequences: tertiary climatic shifts and lineage differentiation times.

BACKGROUND AND AIMS: The genus Olea (Oleaceae) includes approx. 40 taxa of evergreen shrubs and trees classified in three subgenera, Olea, Paniculatae and Tetrapilus, the first of which has two sections (Olea and Ligustroides). Olive trees (the O. europaea complex) have been the subject of intensive research, whereas little is known about the phylogenetic relationships among the other species. To clarify the biogeographical history of this group, a molecular analysis of Olea and related genera of Oleaceae is thus necessary. METHODS: A phylogeny was built of Olea and related genera based on sequences of the nuclear ribosomal internal transcribed spacer-1 and four plastid regions. Lineage divergence and the evolution of abaxial peltate scales, the latter character linked to drought adaptation, were dated using a Bayesian method. KEY RESULTS: Olea is polyphyletic, with O. ambrensis and subgenus Tetrapilus not sharing a most recent common ancestor with the main Olea clade. Partial incongruence between nuclear and plastid phylogenetic reconstructions suggests a reticulation process in the evolution of subgenus Olea. Estimates of divergence times for major groups of Olea during the Tertiary were obtained. CONCLUSIONS: This study indicates the necessity of revising current taxonomic boundaries in Olea. The results also suggest that main lines of evolution were promoted by major Tertiary climatic shifts: (1) the split between subgenera Olea and Paniculatae appears to have taken place at the Miocene-Oligocene boundary; (2) the separation of sections Ligustroides and Olea may have occurred during the Early Miocene following the Mi-1 glaciation; and (3) the diversification within these sections (and the origin of dense abaxial indumentum in section Olea) was concomitant with the aridification of Africa in the Late Miocene.

Cryptic vicariance in Gulf of California fishes parallels vicariant patterns found in Baja, California mammals and reptiles

Comparisons across multiple taxa can often clarify the histories of biogeographic regions. In particular, historic barriers to movement should affect multiple species and, thus, result in a pattern of concordant intraspecific genetic divisions among species. A striking example of such comparative phylogeography is the recent observation that populations of many small mammals and reptiles living on the Baja, California peninsula have a large genetic break between northern and southern peninsular populations. In the present study, I demonstrate that five species of near-shore fishes living on the Baja coastline of the Gulf of California share this genetic pattern. The simplest explanation for this concordant genetic division within both terrestrial and marine vertebrates is that the Baja peninsula was fragmented by a Plio-Pleistocene marine seaway and that this seaway posed a substantial barrier to movement for near-shore fishes. The genetic divisions within Gulf of California fishes also coincide with recognized biogeographic regions based on fish community composition and several environmental factors. It is likely that adaptation to regional environments and present-day oceanographic circulation limits gene exchange between biogeographic regions and helps maintain evidence of past vicariance.

Differential patterns of Male and Female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.

Allochronic taxa as an alternative model to explain circumantarctic disjunctions

Most biogeographical studies propose that southern temperate faunal disjunctions are either the result of vicariance of taxa originated in Gondwana or the result of transoceanic dispersal of taxa originated after the breakup of Gondwana. The aim of this paper is to show that this is a false dichotomy. Antarctica retained a mild climate until mid-Cenozoic and had lasting connections, notably with southern South America and Australia. Both taxa originally Gondwanan and taxa secondarily on Gondwanan areas were subjected to tectonic-induced vicariance, and there is no need to invoke ad hoc transoceanic dispersal, even for post-Gondwanan taxa. These different elements with circumantarctic distributions are here called `allochronic taxa` - taxa presently occupying the same area, but whose presence in that area does not belong to the same time period. This model allows accommodation of conflicting sources of evidence now available for many groups with circumantarctic distributions. The fact that the species from both layers are mixed up in the current biodiversity implies the need to use additional sources of evidence - such as biogeographical, palaeontological, geological and molecular - to discriminate which are the original Gondwanan and which are post-Gondwanan elements in austral landmasses.

Radiation of the Australian flora: What can comparisons of molecular phylogenies across multiple taxa tell us about the evolution of diversity in present-day communities?

The Australian fossil record shows that from ca. 25 Myr ago, the aseasonal-wet biome (rainforest and wet heath) gave way to the unique Australian sclerophyll biomes dominated by eucalypts, acacias and casuarinas. This transition coincided with tectonic isolation of Australia, leading to cooler, drier, more seasonal climates. From 3 Myr ago, aridification caused rapid opening of the central Australian and zone. Molecular phylogenies with dated nodes have provided new perspectives on how these events could have affected the evolution of the Australian flora. During the Mid-Cenozoic (25-10 Myr ago) period of climatic change, there were rapid radiations in sclerophyll taxa, such as Banksia, eucalypts, pea-flowered legumes and Allocasuarina. At the same time, taxa restricted to the aseasonal-wet biome (Nothofagus, Podocarpaceae and Araucariaceae) did not radiate or were depleted by extinction. During the Pliocene aridification, two Eremean biome taxa (Lepidium and Chenopodiaceae) radiated rapidly after dispersing into Australia from overseas. It is clear that the biomes have different histories. Lineages in the aseasonal-wet biome are species poor, with sister taxa that are species rich, either outside Australia or in the sclerophyll biomes. In conjunction with the fossil record, this indicates depletion of the Australian aseasonal-wet biome from the Mid-Cenozoic. In the sclerophyll biomes, there have been multiple exchanges between the southwest and southeast, rather than single large endemic radiations after a vicariance event. There is need for rigorous molecular phylogenetic studies so that additional questions can be addressed, such as how interactions between biomes may have driven the speciation process during radiations. New studies should include the hither-to neglected monsoonal tropics.

Not so ancient: the extant crown group of Nothofagus represents a post-Gondwanan radiation

This study uses a molecular-dating approach to test hypotheses about the biogeography of Nothofagus. The molecular modelling suggests that the present-day subgenera and species date from a radiation that most likely commenced between 55 and 40 Myr ago. This rules out the possibility of a reconciled all-vicariance hypothesis for the biogeography of extant Nothofagus. However, the molecular dates for divergences between Australasian and South American taxa are consistent with the rifting of Australia and South America from Antarctica. The molecular dates further suggest a dispersal of subgenera Lophozonia and Fuscospora between Australia and New Zealand after the onset of the Antarctic Circumpolar Current and west wind drift. It appears likely that the New Caledonian lineage of subgenus Brassospora diverged from the New Guinean lineage elsewhere, prior to colonizing New Caledonia. The molecular approach strongly supports fossil-based estimates that Nothofagus diverged from the rest of Fagales more than 84 Myr ago. However, the mid-Cenozoic estimate for the diversification of the four extant subgenera conflicts with the palynological interpretation because pollen fossils, attributed to all four extant subgenera, were widespread across the Weddellian province of Gondwana about 71 Myr ago. The discrepancy between the pollen and molecular dates exists even when confidence intervals from several sources of error are taken into account. In contrast, the molecular age estimates are consistent with macrofossil dates. The incongruence between pollen fossils and molecular dates could be resolved if the early pollen types represent extinct lineages, with similar types later evolving independently in the extant lineages.

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus ( Sigmodontinae), Dobrava virus ( Murinae), Puumala virus ( Arvicolinae), and Tula virus ( Arvicolinae). Our results reveal that hantaviruses exhibit shortterm substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

Evolutionary history of dog rabies in Brazil

Although dogs are considered to be the principal transmitter of rabies in Brazil, dog rabies had never been recorded in South America before European colonization. In order to investigate the evolutionary history of dog rabies virus (RABV) in Brazil, we performed a phylogenetic analysis of carnivore RABV isolates from around the world and estimated the divergence times for dog RABV in Brazil. Our estimate for the time of introduction of dog RABV into Brazil was the late-19th to early-20th century, which was later than the colonization period but corresponded to a period of increased immigration from Europe to Brazil. In addition, dog RABVs appeared to have spread to indigenous animals in Brazil during the latter half of the 20th century, when the development and urbanization of Brazil occurred. These results suggest that the movement of rabid dogs, along with human activities since the 19th century, promoted the introduction and expansion of dog RABV in Brazil.

Oligocene CO2 decline promoted C4 photosynthesis in grasses.

C4 photosynthesis is an adaptation derived from the more common C3 photosynthetic pathway that confers a higher productivity under warm temperature and low atmospheric CO2 concentration [1, 2]. C4 evolution has been seen as a consequence of past atmospheric CO2 decline, such as the abrupt CO2 fall 32-25 million years ago (Mya) [3-6]. This relationship has never been tested rigorously, mainly because of a lack of accurate estimates of divergence times for the different C4 lineages [3]. In this study, we inferred a large phylogenetic tree for the grass family and estimated, through Bayesian molecular dating, the ages of the 17 to 18 independent grass C4 lineages. The first transition from C3 to C4 photosynthesis occurred in the Chloridoideae subfamily, 32.0-25.0 Mya. The link between CO2 decrease and transition to C4 photosynthesis was tested by a novel maximum likelihood approach. We showed that the model incorporating the atmospheric CO2 levels was significantly better than the null model, supporting the importance of CO2 decline on C4 photosynthesis evolvability. This finding is relevant for understanding the origin of C4 photosynthesis in grasses, which is one of the most successful ecological and evolutionary innovations in plant history.

The Evolution of Trypanosomes Infecting Humans and Primates

Based on phylogenetic analysis of 18S rRNA sequences and clade taxon composition, this paper adopts a biogeographical approach to understanding the evolutionary relationships of the human and primate infective trypanosomes, Trypanosoma cruzi, T. brucei, T. rangeli and T. cyclops. Results indicate that these parasites have divergent origins and fundamentally different patterns of evolution. T. cruzi is placed in a clade with T. rangeli and trypanosomes specific to bats and a kangaroo. The predominantly South American and Australian origins of parasites within this clade suggest an ancient southern super-continent origin for ancestral T. cruzi, possibly in marsupials. T. brucei clusters exclusively with mammalian, salivarian trypanosomes of African origin, suggesting an evolutionary history confined to Africa, while T. cyclops, from an Asian primate appears to have evolved separately and is placed in a clade with T. (Megatrypanum) species. Relating clade taxon composition to palaeogeographic evidence, the divergence of T. brucei and T. cruzi can be dated to the mid-Cretaceous, around 100 million years before present, following the separation of Africa, South America and Euramerica. Such an estimate of divergence time is considerably more recent than those of most previous studies based on molecular clock methods. Perhaps significantly, Salivarian trypanosomes appear, from these data, to be evolving several times faster than Schizotrypanum species, a factor which may have contributed to previous anomalous estimates of divergence times.