BEAST: Bayesian evolutionary analysis by sampling trees

Rabies virus in raccoons, Ohio, 2004

In 2004, the raccoon rabies virus variant emerged in Ohio beyond an area where oral rabies vaccine had been distributed to prevent westward spread of this variant. Our genetic investigation indicates that this outbreak may have begun several years before 2004 and may have originated within the vaccination zone.

Langerhans cells in porcine skin

Langerhans cells (LCs) are resident dendritic cells (DCs) of skin and mucosal epithelium. The standard for identifying skin DCs as LCs is expression of langerin (CD207), a surface protein that mediates Birbeck granule (BG) formation upon internalization. Reports of BGs in porcine skin DC are contradictory, due to lack of langerin detection. Here, we present the sequence of porcine langerin/CD207, showing that the predicted porcine protein shares 75%/86% amino acid identity/similarity with human. Langerin mRNA was detected in porcine skin DCs by PCR and langerin protein was detected in both isolated skin DCs and skin sections by immunostaining. Approximately, 50-70% of skin DCs expressed langerin, demonstrating that the majority of porcine skin DCs are LCs. The full length sequence combined with the identification of antibodies reactive with porcine langerin, facilitates the study of LCs in swine, and advances the use of swine for studying skin diseases and infectious disease processes involving skin.

Global distribution of novel rhinovirus genotype

Global surveillance for a novel rhinovirus genotype indicated its association with community outbreaks and pediatric respiratory disease in Africa, Asia, Australia, Europe, and North America. Molecular dating indicates that these viruses have been circulating for at least 250 years.

Mitochondrial genomes reveal an explosive radiation of extinct and extant bears near the Miocene-Pliocene boundary

BACKGROUND

Despite being one of the most studied families within the Carnivora, the phylogenetic relationships among the members of the bear family (Ursidae) have long remained unclear. Widely divergent topologies have been suggested based on various data sets and methods.

RESULTS

We present a fully resolved phylogeny for ursids based on ten complete mitochondrial genome sequences from all eight living and two recently extinct bear species, the European cave bear (Ursus spelaeus) and the American giant short-faced bear (Arctodus simus). The mitogenomic data yield a well-resolved topology for ursids, with the sloth bear at the basal position within the genus Ursus. The sun bear is the sister taxon to both the American and Asian black bears, and this clade is the sister clade of cave bear, brown bear and polar bear confirming a recent study on bear mitochondrial genomes.

CONCLUSION

Sequences from extinct bears represent the third and fourth Pleistocene species for which complete mitochondrial genomes have been sequenced. Moreover, the cave bear specimen demonstrates that mitogenomic studies can be applied to Pleistocene fossils that have not been preserved in permafrost, and therefore have a broad application within ancient DNA research. Molecular dating of the mtDNA divergence times suggests a rapid radiation of bears in both the Old and New Worlds around 5 million years ago, at the Miocene-Pliocene boundary. This coincides with major global changes, such as the Messinian crisis and the first opening of the Bering Strait, and suggests a global influence of such events on species radiations.

The effect of ancient DNA damage on inferences of demographic histories

The field of ancient DNA (aDNA) is casting new light on many evolutionary questions. However, problems associated with the postmortem instability of DNA may complicate the interpretation of aDNA data. For example, in population genetic studies, the inclusion of damaged DNA may inflate estimates of diversity. In this paper, we examine the effect of DNA damage on population genetic estimates of ancestral population size. We simulate data using standard coalescent simulations that include postmortem damage and show that estimates of effective population sizes are inflated around, or right after, the sampling time of the ancestral DNA sequences. This bias leads to estimates of increasing, and then decreasing, population sizes, as observed in several recently published studies. We reanalyze a recently published data set of DNA sequences from the Bison (Bison bison/Bison priscus) and show that the signal for a change in effective population size in this data set vanishes once the effects of putative damage are removed. Our results suggest that population genetic analyses of aDNA sequences, which do not accurately account for damage, should be interpreted with great caution.

Radiation of southern African daisies: biogeographic inferences for subtribe Arctotidinae (Asteraceae, Arctotideae)

The majority of the approximately 80-90 species in subtribe Arctotidinae occur in southern Africa with the centre of diversity in the winter-rainfall region. Three species are restricted to afromontane eastern Africa and three species are endemic to Australia. To investigate biogeographic and phylogenetic relationships within Arctotidinae, sequence data from four cpDNA regions (psbA-trnH, trnT-trnL and trnL-trnF spacers and trnL intron) and the ITS nrDNA region for 59 Arctotidinae species were analyzed with parsimony and Bayesian-inference approaches. Eight well-supported major lineages were resolved. The earliest-diverging extant lineages are afromontane or inhabit mesic habitats, whereas almost all sampled taxa from the winter-rainfall and semi-arid areas have diverged more recently. Molecular dating estimated that the major clades diverged during the Miocene and Pliocene, which is coincident with the trend of increasing rainfall seasonality, aridification and vegetation changes in southwestern Africa. Trans-oceanic dispersal to Australia was estimated to have occurred during the Pliocene.

Molecular ecology and natural history of simian foamy virus infection in wild-living chimpanzees

Identifying microbial pathogens with zoonotic potential in wild-living primates can be important to human health, as evidenced by human immunodeficiency viruses types 1 and 2 (HIV-1 and HIV-2) and Ebola virus. Simian foamy viruses (SFVs) are ancient retroviruses that infect Old and New World monkeys and apes. Although not known to cause disease, these viruses are of public health interest because they have the potential to infect humans and thus provide a more general indication of zoonotic exposure risks. Surprisingly, no information exists concerning the prevalence, geographic distribution, and genetic diversity of SFVs in wild-living monkeys and apes. Here, we report the first comprehensive survey of SFVcpz infection in free-ranging chimpanzees (Pan troglodytes) using newly developed, fecal-based assays. Chimpanzee fecal samples (n = 724) were collected at 25 field sites throughout equatorial Africa and tested for SFVcpz-specific antibodies (n = 706) or viral nucleic acids (n = 392). SFVcpz infection was documented at all field sites, with prevalence rates ranging from 44% to 100%. In two habituated communities, adult chimpanzees had significantly higher SFVcpz infection rates than infants and juveniles, indicating predominantly horizontal rather than vertical transmission routes. Some chimpanzees were co-infected with simian immunodeficiency virus (SIVcpz); however, there was no evidence that SFVcpz and SIVcpz were epidemiologically linked. SFVcpz nucleic acids were recovered from 177 fecal samples, all of which contained SFVcpz RNA and not DNA. Phylogenetic analysis of partial gag (616 bp), pol-RT (717 bp), and pol-IN (425 bp) sequences identified a diverse group of viruses, which could be subdivided into four distinct SFVcpz lineages according to their chimpanzee subspecies of origin. Within these lineages, there was evidence of frequent superinfection and viral recombination. One chimpanzee was infected by a foamy virus from a Cercopithecus monkey species, indicating cross-species transmission of SFVs in the wild. These data indicate that SFVcpz (i) is widely distributed among all chimpanzee subspecies; (ii) is shed in fecal samples as viral RNA; (iii) is transmitted predominantly by horizontal routes; (iv) is prone to superinfection and recombination; (v) has co-evolved with its natural host; and (vi) represents a sensitive marker of population structure that may be useful for chimpanzee taxonomy and conservation strategies.

Temporal and spatial dynamics of human immunodeficiency virus type 1 circulating recombinant forms 08_BC and 07_BC in Asia

Human immunodeficiency virus type 1 (HIV-1) CRF08_BC and CRF07_BC are two major recombinants descended from subtypes B' and C. Despite their massive epidemic impact in China, their migration patterns and divergence times remain unknown. Phylogenetic and population genetic analyses were performed on 228 HIV-1 sequences representing CRF08_BC, CRF07_BC, and subtype C strains from different locations across China, India, and Myanmar. Genome-specific rates of evolution and divergence times were estimated using a Bayesian Markov chain Monte Carlo framework under various evolutionary models. CRF08_BC originated in 1990.3 (95% credible region [CR], 1988.6 to 1991.9) in Yunnan province before spreading to Guangxi (south) and Liaoning (northeast) around 1995. Inside Guangxi region, the eastward expansion of CRF08_BC continued from Baise city (west) to Binyang (central) between 1997 and 1998 and later spread into Pingxiang around 1999 in the south, mainly through injecting drug users. Additionally, CRF07_BC diverged from its common ancestor in 1993.3 (95% CR, 1991.2 to 1995.2) before crossing the border into southern Taiwan in late 1990s. Phylogenetic analysis indicates that both CRF08_BC and CRF07_BC can trace their origins to Yunnan. The parental Indian subtype C lineage likely entered China around 1981.2 (95% CR, 1976.7 to 1985.9). Using a multiple unlinked locus model, we also showed that the dates of divergence calculated in this study may not be significantly affected by intrasubtype recombination among different lineages. This is the first phylodynamic study depicting the spatiotemporal dynamics of HIV/AIDS in East Asia.

The prehistory of potyviruses: their initial radiation was during the dawn of agriculture

BACKGROUND

Potyviruses are found world wide, are spread by probing aphids and cause considerable crop damage. Potyvirus is one of the two largest plant virus genera and contains about 15% of all named plant virus species. When and why did the potyviruses become so numerous? Here we answer the first question and discuss the other.

METHODS AND FINDINGS

We have inferred the phylogenies of the partial coat protein gene sequences of about 50 potyviruses, and studied in detail the phylogenies of some using various methods and evolutionary models. Their phylogenies have been calibrated using historical isolation and outbreak events: the plum pox virus epidemic which swept through Europe in the 20th century, incursions of potyviruses into Australia after agriculture was established by European colonists, the likely transport of cowpea aphid-borne mosaic virus in cowpea seed from Africa to the Americas with the 16th century slave trade and the similar transport of papaya ringspot virus from India to the Americas.

CONCLUSIONS/SIGNIFICANCE

Our studies indicate that the partial coat protein genes of potyviruses have an evolutionary rate of about 1.15x10(-4) nucleotide substitutions/site/year, and the initial radiation of the potyviruses occurred only about 6,600 years ago, and hence coincided with the dawn of agriculture. We discuss the ways in which agriculture may have triggered the prehistoric emergence of potyviruses and fostered their speciation.

Rapid evolutionary dynamics of zucchini yellow mosaic virus

Zucchini yellow mosaic virus (ZYMV) is an economically important virus of cucurbit crops. However, little is known about the rate at which this virus has evolved within members of the family Cucurbitaceae, or the timescale of its epidemiological history. Herein, we present the first analysis of the evolutionary dynamics of ZYMV. Using a Bayesian coalescent approach we show that the coat protein of ZYMV has evolved at a mean rate of 5.0 x 10(-4) nucleotide substitutions per site, per year. Notably, this rate is equivalent to those observed in animal RNA viruses. Using the same approach we show that the lineages of ZYMV sampled here have an ancestry that dates back no more than 800 years, suggesting that human activities have played a central role in the dispersal of ZYMV. Finally, an analysis of phylogeographical structure provides strong evidence for the in situ evolution of ZYMV within individual countries.

Diversification of myco-heterotrophic angiosperms: evidence from Burmanniaceae

BACKGROUND

Myco-heterotrophy evolved independently several times during angiosperm evolution. Although many species of myco-heterotrophic plants are highly endemic and long-distance dispersal seems unlikely, some genera are widely dispersed and have pantropical distributions, often with large disjunctions. Traditionally this has been interpreted as evidence for an old age of these taxa. However, due to their scarcity and highly reduced plastid genomes our understanding about the evolutionary histories of the angiosperm myco-heterotrophic groups is poor.

RESULTS

We provide a hypothesis for the diversification of the myco-heterotrophic family Burmanniaceae. Phylogenetic inference, combined with biogeographical analyses, molecular divergence time estimates, and diversification analyses suggest that Burmanniaceae originated in West Gondwana and started to diversify during the Late Cretaceous. Diversification and migration of the species-rich pantropical genera Burmannia and Gymnosiphon display congruent patterns. Diversification began during the Eocene, when global temperatures peaked and tropical forests occurred at low latitudes. Simultaneous migration from the New to the Old World in Burmannia and Gymnosiphon occurred via boreotropical migration routes. Subsequent Oligocene cooling and breakup of boreotropical flora ended New-Old World migration and caused a gradual decrease in diversification rate in Burmanniaceae.

CONCLUSION

Our results indicate that extant diversity and pantropical distribution of myco-heterotrophic Burmanniaceae is the result of diversification and boreotropical migration during the Eocene when tropical rain forest expanded dramatically.

Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.

Intraspecific patterns of mitochondrial variation in natural population fragments of a localized desert dung beetle species, Pachysoma gariepinum (Coleoptera: Scarabaeidae)

Phylogenetic, population, and coalescent methods were used to examine the genetic structuring of Pachysoma gariepinum, a flightless dung beetle species endemic to the arid west coast of southern Africa that exhibits interrupted south to north morphological clinal variation along a distributional gradient. Mitochondrial cytochrome oxidase I sequence data of 67 individuals from 5 localities revealed the presence of 3 geographically distinct evolutionary lineages (with an overall nucleotide divergence of 5.7% and a per-locality divergence of 1.9-3.8%) which display significant levels of genetic structuring. The separation of the lineages was estimated to have occurred between 2.2 and 5.7 million years ago--which is the late Miocene, early Plio-Pleistocene era--possibly in response to the ebb and flow of the Orange and Holgat River systems as well as the interactions between the moving and stable sand dune systems. Moreover the species' current range appears to have been influenced by the formation of advective fog resulting in a constant source of water in an area with low precipitation thereby allowing for the beetles to radiate to areas that were previously inhospitable. Fu's F-statistics and population parameters based on recent mutations indicated that little to no recent population growth has occurred. This together with changing anthropogenic factors and the recovery of 3 geographically discrete management units, points to a need for census data in order to monitor and conserve the genetic diversity of this species.

Molecular phylogeny of reed beetles (Col., Chrysomelidae, Donaciinae): the signature of ecological specialization and geographical isolation

The Donaciinae consist of approximately 165 species predominantly occurring in the northern hemisphere. We analysed mitochondrial and nuclear DNA (COI, EF-1alpha) of 46 species to investigate their phylogeny and to discuss general topics in the context of insect herbivory (generalists versus specialists, ecological speciation). Phylogenetic reconstructions from various methodical approaches yielded very similar results. Clades corresponding to the traditional tribes/genera were recovered. Within the genus Donacia, species groups with characteristic host plant preference were identified. Estimated divergence times are discussed on the background of geological events. The origin of the Donaciinae is dated to 75-100 million years before present, after which they quickly diversified into the main groups. An initial split of those groups occurred in the Palaeocene. In the Eocene and Oligocene, major lineages specialized on certain host plants, where they radiated in the Miocene. This radiation was enforced by geographic isolation brought about by the final separation of America and Europe, after which there arose continental lineages within three larger species groups. In their evolution based on ecological specialization with a recently superimposed geographic isolation, the Donaciinae follow a pattern of specialists arising from generalists. Host plant shifts show that such a specialization is not necessarily an 'evolutionary dead-end'.

Density-dependent cladogenesis in birds

A characteristic signature of adaptive radiation is a slowing of the rate of speciation toward the present. On the basis of molecular phylogenies, studies of single clades have frequently found evidence for a slowdown in diversification rate and have interpreted this as evidence for density dependent speciation. However, we demonstrated via simulation that large clades are expected to show stronger slowdowns than small clades, even if the probability of speciation and extinction remains constant through time. This is a consequence of exponential growth: clades, which, by chance, diversify at above the average rate early in their history, will tend to be large. They will also tend to regress back to the average diversification rate later on, and therefore show a slowdown. We conducted a meta-analysis of the distribution of speciation events through time, focusing on sequence-based phylogenies for 45 clades of birds. Thirteen of the 23 clades (57%) that include more than 20 species show significant slowdowns. The high frequency of slowdowns observed in large clades is even more extreme than expected under a purely stochastic constant-rate model, but is consistent with the adaptive radiation model. Taken together, our data strongly support a model of density-dependent speciation in birds, whereby speciation slows as ecological opportunities and geographical space place limits on clade growth.

It is probable that the number of species that a given region can support is limited; however, it is unclear whether the limit is approached sufficiently in nature such that the rate at which new species form slows down. Using the pattern of phylogenetic branching, a technique that estimates evolutionary relationships based on molecular data, we demonstrate that in large clades of birds, there is a decrease in the per-lineage probability of speciation as the number of species in the clade increase. We also show that this pattern can arise even if speciation and extinction occur randomly through time. This is because large clades are likely, by chance, to have rapidly speciated early in their history, and will relax back to the average speciation rate later on. We account for this effect, and we still find evidence that, as a clade grows to large size, the per-lineage probability of speciation declines. These results strongly suggest that speciation rates are slowed as environments fill up with competitors.

Molecular evidence provides strong evidence that speciation rates slow down through time.

Rapid sequence change and geographical spread of human parvovirus B19: comparison of B19 virus evolution in acute and persistent infections

Parvovirus B19 is a common human pathogen maintained by horizontal transmission between acutely infected individuals. However, B19 virus can also be detected in tissues throughout the life of the host, although little is understood about the nature of such persistence. In the current study, we created large VP1/2 sequence data sets of plasma- and tissue (autopsy)-derived variants of B19 virus with known sample dates to compare the rates of sequence change in exogenous virus populations with those in persistently infected individuals. By using linear regression and likelihood-based methods (such as the BEAST program), we found that plasma-derived B19 virus showed a substitution rate of 4 x 10(-4) and an unconstrained (synonymous)-substitution rate of 18 x 10(-4) per site per year, several times higher than previously estimated and within the range of values for mammalian RNA viruses. The underlying high mutation frequency implied by these substitution rates may enable rapid adaptive changes that are more commonly ascribed to RNA virus populations. These revised estimates predict that the last common ancestor for currently circulating genotype 1 variants of B19 virus existed around 1956 to 1959, fitting well with previous analyses of the B19 virus "bioportfolio" that support a complete cessation of genotype 2 infections and their replacement by genotype 1 infections in the 1960s. In contrast, the evolution of B19 virus amplified from tissue samples was best modeled by using estimated dates of primary infection rather than sample dates, consistent with slow or absent sequence change during persistence. Determining what epidemiological or biological factors led to such a complete and geographically extensive population replacement over this short period is central to further understanding the nature of parvovirus evolution.

The genomic and epidemiological dynamics of human influenza A virus

The evolutionary interaction between influenza A virus and the human immune system, manifest as 'antigenic drift' of the viral haemagglutinin, is one of the best described patterns in molecular evolution. However, little is known about the genome-scale evolutionary dynamics of this pathogen. Similarly, how genomic processes relate to global influenza epidemiology, in which the A/H3N2 and A/H1N1 subtypes co-circulate, is poorly understood. Here through an analysis of 1,302 complete viral genomes sampled from temperate populations in both hemispheres, we show that the genomic evolution of influenza A virus is characterized by a complex interplay between frequent reassortment and periodic selective sweeps. The A/H3N2 and A/H1N1 subtypes exhibit different evolutionary dynamics, with diverse lineages circulating in A/H1N1, indicative of weaker antigenic drift. These results suggest a sink-source model of viral ecology in which new lineages are seeded from a persistent influenza reservoir, which we hypothesize to be located in the tropics, to sink populations in temperate regions.

Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics

Kingman's coalescent process opens the door for estimation of population genetics model parameters from molecular sequences. One paramount parameter of interest is the effective population size. Temporal variation of this quantity characterizes the demographic history of a population. Because researchers are rarely able to choose a priori a deterministic model describing effective population size dynamics for data at hand, nonparametric curve-fitting methods based on multiple change-point (MCP) models have been developed. We propose an alternative to change-point modeling that exploits Gaussian Markov random fields to achieve temporal smoothing of the effective population size in a Bayesian framework. The main advantage of our approach is that, in contrast to MCP models, the explicit temporal smoothing does not require strong prior decisions. To approximate the posterior distribution of the population dynamics, we use efficient, fast mixing Markov chain Monte Carlo algorithms designed for highly structured Gaussian models. In a simulation study, we demonstrate that the proposed temporal smoothing method, named Bayesian skyride, successfully recovers "true" population size trajectories in all simulation scenarios and competes well with the MCP approaches without evoking strong prior assumptions. We apply our Bayesian skyride method to 2 real data sets. We analyze sequences of hepatitis C virus contemporaneously sampled in Egypt, reproducing all key known aspects of the viral population dynamics. Next, we estimate the demographic histories of human influenza A hemagglutinin sequences, serially sampled throughout 3 flu seasons.

Rice yellow mottle virus, an RNA plant virus, evolves as rapidly as most RNA animal viruses

The rate of evolution of an RNA plant virus has never been estimated using temporally spaced sequence data, by contrast to the information available on an increasing range of animal viruses. Accordingly, the evolution rate of Rice yellow mottle virus (RYMV) was calculated from sequences of the coat protein gene of isolates collected from rice over a 40-year period in different parts of Africa. The evolution rate of RYMV was estimated by pairwise distance linear regression on five phylogeographically defined groups comprising a total of 135 isolates. It was further assessed from 253 isolates collected all over Africa by Bayesian coalescent methods under strict and relaxed molecular clock models and under constant size and skyline population genetic models. Consistent estimates of the evolution rate between 4 x 10(-4) and 8 x 10(-4) nucleotides (nt)/site/year were obtained whatever method and model were applied. The synonymous evolution rate was between 8 x 10(-4) and 11 x 10(-4) nt/site/year. The overall and synonymous evolution rates of RYMV were within the range of the rates of 50 RNA animal viruses, below the average but above the distribution median. Experimentally, in host change studies, substitutions accumulated at an even higher rate. The results show that an RNA plant virus such as RYMV evolves as rapidly as most RNA animal viruses. Knowledge of the molecular clock of plant viruses provides methods for testing a wide range of biological hypotheses.

A Bayesian evaluation of human mitochondrial substitution rates

Accurate estimates of mitochondrial substitution rates are central to molecular studies of human evolution, but meaningful comparisons of published studies are problematic because of the wide range of methodologies and data sets employed. These differences are nowhere more pronounced than among rates estimated from phylogenies, genealogies, and pedigrees. By using a data set comprising mitochondrial genomes from 177 humans, we estimate substitution rates for various data partitions by using Bayesian phylogenetic analysis with a relaxed molecular clock. We compare the effect of multiple internal calibrations with the customary human-chimpanzee split. The analyses reveal wide variation among estimated substitution rates and divergence times made with different partitions and calibrations, with evidence of substitutional saturation, natural selection, and significant rate heterogeneity among lineages and among sites. Collectively, the results support dates for migration out of Africa and the common mitochondrial ancestor of humans that are considerably more recent than most previous estimates. Our results also demonstrate that human mitochondrial genomes exhibit a number of molecular evolutionary complexities that necessitate the use of sophisticated analytical models for genetic analyses.

Major evolutionary transitions in ant agriculture

Agriculture is a specialized form of symbiosis that is known to have evolved in only four animal groups: humans, bark beetles, termites, and ants. Here, we reconstruct the major evolutionary transitions that produced the five distinct agricultural systems of the fungus-growing ants, the most well studied of the nonhuman agriculturalists. We do so with reference to the first fossil-calibrated, multiple-gene, molecular phylogeny that incorporates the full range of taxonomic diversity within the fungus-growing ant tribe Attini. Our analyses indicate that the original form of ant agriculture, the cultivation of a diverse subset of fungal species in the tribe Leucocoprineae, evolved approximately 50 million years ago in the Neotropics, coincident with the early Eocene climatic optimum. During the past 30 million years, three known ant agricultural systems, each involving a phylogenetically distinct set of derived fungal cultivars, have separately arisen from the original agricultural system. One of these derived systems subsequently gave rise to the fifth known system of agriculture, in which a single fungal species is cultivated by leaf-cutter ants. Leaf-cutter ants evolved remarkably recently ( approximately 8-12 million years ago) to become the dominant herbivores of the New World tropics. Our analyses identify relict, extant attine ant species that occupy phylogenetic positions that are transitional between the agricultural systems. Intensive study of those species holds particular promise for clarifying the sequential accretion of ecological and behavioral characters that produced each of the major ant agricultural systems.

How to summarize estimates of ancestral divergence times

The use of molecular sequence data has increased interest in trying to date evolutionary events, with researchers wanting both an estimate of the divergence time and a confidence interval for that estimate. However, two methodological issues have recently been raised with respect to precision of the estimates: (i) the time of the ancestral event is over-estimated; and (ii) the confidence interval is asymmetrical. I argue that if the estimates of divergence time are considered to be samples from a lognormal probability distribution, then this would explain both of these problems. This implies that divergence times should be presented using geometric means rather than arithmetic means, both for estimates and for their confidence intervals. I present analyses based on both computer simulations and empirical data to show that this approach is effective for both single-gene and multiple-gene data sets. Treating divergence time as a lognormal variable thus provides a simple unifying framework for dealing with many of the problems associated with the estimation of divergence (and possibly coalescence) times. Use of this approach (based on geometric means) can, unfortunately, lead to very different biological conclusions compared to the currently used calculation methods (based on arithmetic means).

Tempo and mode of diversification of lake Tanganyika cichlid fishes

BACKGROUND

Understanding the causes of disparities in species diversity across taxonomic groups and regions is a fundamental aim in evolutionary biology. Addressing these questions is difficult because of the need for densely sampled phylogenies and suitable empirical systems.

METHODOLOGY/PRINCIPAL FINDINGS

Here we investigate the cichlid fish radiation of Lake Tanganyika and show that per lineage diversification rates have been more than six times slower than in the species flocks of Lakes Victoria and Malawi. The result holds even at peak periods of diversification in Lake Tanganyika, ruling out the age of the lake as an explanation for slow average rates, and is robust to uncertainties over the calibration of cichlid radiations in geological time. Moreover, Lake Tanganyika lineages, irrespective of different biological characteristics (e.g. sexually dichromatic versus sexually monochromatic clades), have diversified at similar rates, falling within typical estimates across a range of plant and animal clades. For example, the mostly sexually dichromatic haplochromines, which have speciated explosively in Lakes Victoria and Malawi, have displayed modest rates in Lake Tanganyika (where they are called Tropheini).

CONCLUSION/SIGNIFICANCE

Our results show that either the Lake Tanganyika environment is less conducive for cichlid speciation or the remarkable diversifying abilities of the haplochromines were inhibited by the prior occupancy of older radiations. Although the results indicate a dominant role for the environment in shaping cichlid diversification, differences in the timing of diversification among the Tanganyikan tribes indicate that biological differences were still important for the dynamics of species build-up in the lake. While we cannot resolve the timing of the radiation relative to the origin of the lake, because of the lack of robust geological date calibrations for cichlids, our results are consistent with a scenario that the different clades reflect independent adaptive radiations into different broad niches in the lake.

Bayesian coalescent analysis reveals a high rate of molecular evolution in GB virus C

GB virus C/hepatitis G (GBV-C) is an RNA virus of the family Flaviviridae. Despite replicating with an RNA-dependent RNA polymerase, some previous estimates of rates of evolutionary change in GBV-C suggest that it fixes mutations at the anomalously low rate of approximately 10(-7) nucleotide substitution per site, per year. However, these estimates were largely based on the assumption that GBV-C and its close relative GBV-A (New World monkey GB viruses) codiverged with their primate hosts over millions of years. Herein, we estimated the substitution rate of GBV-C using the largest set of dated GBV-C isolates compiled to date and a Bayesian coalescent approach that utilizes the year of sampling and so is independent of the assumption of codivergence. This revealed a rate of evolutionary change approximately four orders of magnitude higher than that estimated previously, in the range of 10(-2) to 10(-3) sub/site/year, and hence in line with those previously determined for RNA viruses in general and the Flaviviridae in particular. In addition, we tested the assumption of host-virus codivergence in GBV-A by performing a reconciliation analysis of host and virus phylogenies. Strikingly, we found no statistical evidence for host-virus codivergence in GBV-A, indicating that substitution rates in the GB viruses should not be estimated from host divergence times.

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.

Post-Messinian evolutionary relationships across the Sicilian channel: mitochondrial and nuclear markers link a new green toad from Sicily to African relatives

BACKGROUND

Little attention has been paid to the consequences of the last landbridge between Africa and Sicily on Mediterranean biogeography. Previous paleontological and scarce molecular data suggest possible faunal exchange later than the well-documented landbridge in the Messinian (5.3 My); however, a possible African origin of recent terrestrial Sicilian fauna has not been thoroughly tested with molecular methods. To gain insight into the phylogeography of the region, we examine two mitochondrial and two nuclear markers (one is a newly adapted intron marker) in green toads (Bufo viridis subgroup) across that sea barrier, the Strait of Sicily.

RESULTS

Extensive sampling throughout the western Mediterranean and North Africa revealed a deep sister relationship between Sicilian (Bufo siculus n.sp.) and African green toads (B. boulengeri) on the mitochondrial and nuclear level. Divergence times estimated under a Bayesian-coalescence framework (mtDNA control region and 16S rRNA) range from the Middle Pliocene (3.6 My) to Pleistocene (0.16 My) with an average (1.83 to 2.0 My) around the Pliocene/Pleistocene boundary, suggesting possible land connections younger than the Messinian (5.3 My). We describe green toads from Sicily and some surrounding islands as a new endemic species (Bufo siculus). Bufo balearicus occurs on some western Mediterranean islands (Corsica, Sardinia, Mallorca, and Menorca) and the Apennine Peninsula, and is well differentiated on the mitochondrial and nuclear level from B. siculus as well as from B. viridis (Laurenti), whose haplotype group reaches northeastern Italy, north of the Po River. Detection of Calabrian B. balearicus haplotypes in northeastern Sicily suggests recent invasion. Our data agree with paleogeographic and fossil data, which suggest long Plio-Pleistocene isolation of Sicily and episodic Pleistocene faunal exchange across the Strait of Messina. It remains unknown whether both species (B. balearicus, B. siculus) occur in sympatry in northern Sicily.

CONCLUSION

Our findings on green toads give the first combined mitochondrial and nuclear sequence evidence for a phylogeographic connection across the Strait of Sicily in terrestrial vertebrates. These relationships may have implications for comparative phylogeographic research on other terrestrial animals co-occurring in North Africa and Sicily.

Grafting the molecular phylogenetic tree with morphological branches to reconstruct the evolutionary history of the genus Zaprionus (Diptera: Drosophilidae)

A molecular phylogeny for the drosophilid genus Zaprionus was inferred using a mitochondrial (CO-II) and a nuclear (Amyrel) gene using 22 available species. The combined molecular tree does not support the current classification, dubbed phylogenetic, based entirely upon a morphocline of forefemoral ornamentation. For species for which DNA was not available, phylogenetic positioning was only assigned using morphological characters. In order to avoid conflict between DNA and morphology in the combined analyses (supermatrix method), we developed a new method in which few morphological characters were sampled according to an a priori homoplasy assessment on the consensus molecular tree. At each internal node of the tree, a number of synapomorphies was determined, and species with no molecular sequences were grafted thereon. Analogously to tree vocabulary, we called our method 'morphological grafting'. New species groups and complexes were then defined in the light of our findings. Further, divergence times were estimated under a relaxed molecular clock, and historical biogeography was reconstructed under a maximum likelihood model. Zaprionus appears to be of recent origin in the Oriental region during the Late Miocene ( approximately 10 MYA), and colonization of Africa started shortly after ( approximately 7 MYA) via the maritime route of the Indian Ocean Islands. Most of the morphological and ecological diversification took place, later, in Western Africa during the Quaternary cyclic climatic changes. Furthermore, some species became recent invaders, with one, Zaprionus indianus, has successfully invaded South and North America during the last decade.

Shrimps down under: evolutionary relationships of subterranean crustaceans from Western Australia (Decapoda: Atyidae: Stygiocaris)

BACKGROUND

We investigated the large and small scale evolutionary relationships of the endemic Western Australian subterranean shrimp genus Stygiocaris (Atyidae) using nuclear and mitochondrial genes. Stygiocaris is part of the unique cave biota of the coastal, anchialine, limestones of the Cape Range and Barrow Island, most of whose nearest evolutionary relations are found in coastal caves of the distant North Atlantic. The dominance of atyids in tropical waters and their food resources suggest they are pivotal in understanding these groundwater ecosystems.

METHODOLOGY/PRINCIPLE FINDINGS

Our nuclear and mitochondrial analyses all recovered the Mexican cave genus Typhlatya as the sister taxon of Stygiocaris, rather than any of the numerous surface and cave atyids from Australia or the Indo-Pacific region. The two described Stygiocaris species were recovered as monophyletic, and a third, cryptic, species was discovered at a single site, which has very different physiochemical properties from the sites hosting the two described species.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that Stygiocaris and Typhlatya may descend from a common ancestor that lived in the coastal marine habitat of the ancient Tethys Sea, and were subsequently separated by plate tectonic movements. This vicariant process is commonly thought to explain the many disjunct anchialine faunas, but has rarely been demonstrated using phylogenetic techniques. The Cape Range's geological dynamism, which is probably responsible for the speciation of the various Stygiocaris species, has also led to geographic population structure within species. In particular, Stygiocaris lancifera is split into northern and southern groups, which correspond to population splits within other sympatric subterranean taxa.

Pandemic human viruses cause decline of endangered great apes

Commercial hunting and habitat loss are major drivers of the rapid decline of great apes [1]. Ecotourism and research have been widely promoted as a means of providing alternative value for apes and their habitats [2]. However, close contact between humans and habituated apes during ape tourism and research has raised concerns that disease transmission risks might outweigh benefits [3-7]. To date only bacterial and parasitic infections of typically low virulence have been shown to move from humans to wild apes [8, 9]. Here, we present the first direct evidence of virus transmission from humans to wild apes. Tissue samples from habituated chimpanzees that died during three respiratory-disease outbreaks at our research site, Côte d'Ivoire, contained two common human paramyxoviruses. Viral strains sampled from chimpanzees were closely related to strains circulating in contemporaneous, worldwide human epidemics. Twenty-four years of mortality data from observed chimpanzees reveal that such respiratory outbreaks could have a long history. In contrast, survey data show that research presence has had a strong positive effect in suppressing poaching around the research site. These observations illustrate the challenge of maximizing the benefit of research and tourism to great apes while minimizing the negative side effects.

Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages

BACKGROUND

Determining an absolute timescale for avian evolutionary history has proven contentious. The two sources of information available, paleontological data and inference from extant molecular genetic sequences (colloquially, 'rocks' and 'clocks'), have appeared irreconcilable; the fossil record supports a Cenozoic origin for most modern lineages, whereas molecular genetic estimates suggest that these same lineages originated deep within the Cretaceous and survived the K-Pg (Cretaceous-Paleogene; formerly Cretaceous-Tertiary or K-T) mass-extinction event. These two sources of data therefore appear to support fundamentally different models of avian evolution. The paradox has been speculated to reflect deficiencies in the fossil record, unrecognized biases in the treatment of genetic data or both. Here we attempt to explore uncertainty and limit bias entering into molecular divergence time estimates through: (i) improved taxon (n = 135) and character (n = 4594 bp mtDNA) sampling; (ii) inclusion of multiple cladistically tested internal fossil calibration points (n = 18); (iii) correction for lineage-specific rate heterogeneity using a variety of methods (n = 5); (iv) accommodation of uncertainty in tree topology; and (v) testing for possible effects of episodic evolution.

RESULTS

The various 'relaxed clock' methods all indicate that the major (basal) lineages of modern birds originated deep within the Cretaceous, although temporal intraordinal diversification patterns differ across methods. We find that topological uncertainty had a systematic but minor influence on date estimates for the origins of major clades, and Bayesian analyses assuming fixed topologies deliver similar results to analyses with unconstrained topologies. We also find that, contrary to expectation, rates of substitution are not autocorrelated across the tree in an ancestor-descendent fashion. Finally, we find no signature of episodic molecular evolution related to either speciation events or the K-Pg boundary that could systematically mislead inferences from genetic data.

CONCLUSION

The 'rock-clock' gap has been interpreted by some to be a result of the vagaries of molecular genetic divergence time estimates. However, despite measures to explore different forms of uncertainty in several key parameters, we fail to reconcile molecular genetic divergence time estimates with dates taken from the fossil record; instead, we find strong support for an ancient origin of modern bird lineages, with many extant orders and families arising in the mid-Cretaceous, consistent with previous molecular estimates. Although there is ample room for improvement on both sides of the 'rock-clock' divide (e.g. accounting for 'ghost' lineages in the fossil record and developing more realistic models of rate evolution for molecular genetic sequences), the consistent and conspicuous disagreement between these two sources of data more likely reflects a genuine difference between estimated ages of (i) stem-group origins and (ii) crown-group morphological diversifications, respectively. Further progress on this problem will benefit from greater communication between paleontologists and molecular phylogeneticists in accounting for error in avian lineage age estimates.

Climate, phylogeny and the ecological distribution of C4 grasses

'C4 photosynthesis' refers to a suite of traits that increase photosynthesis in high light and high temperature environments. Most C4 plants are grasses, which dominate tropical and subtropical grasslands and savannas but are conspicuously absent from cold growing season climates. Physiological attributes of C4 photosynthesis have been invoked to explain C4 grass biogeography; however, the pathway evolved exclusively in grass lineages of tropical origin, suggesting that the prevalence of C4 grasses in warm climates could be due to other traits inherited from their non-C4 ancestors. Here we investigate the relative influences of phylogeny and photosynthetic pathway in determining the ecological distributions of C4 grasses in Hawaii. We find that the restriction of C4 grasses to warmer areas is due largely to their evolutionary history as members of a warm-climate grass clade, but that the pathway does appear to confer a competitive advantage to grasses in more arid environments.

First DNA sequences from Asian cave bear fossils reveal deep divergences and complex phylogeographic patterns

Until recently, cave bears were believed to have only inhabited Europe. However, recent morphological evidence suggests that cave bears' geographic range extended as far east as Transbaikalia, Eastern Siberia. These Asian cave bears were morphologically distinct from European cave bears. However, how they related to European lineages remains unclear, stressing the need to assess the phylogenetic and phylogeographic relationship between Asian cave bears and their European relatives. In this work, we address this issue using a 227 base-pair fragment of the mitochondrial control region obtained from nine fossil bone samples from eight sites from the Urals, Caucasus, Altai Mountains, Ukraine and Yana River region in Eastern Siberia. Results of the phylogenetic analyses indicate that (i) the cave bear from the Yana River is most closely related to cave bears from the Caucasus region; (ii) the Caucasus/Yana group of bears is genetically very distinct from both European cave bears and brown bears, suggesting that these bears could represent an independent species; and (iii) the Western European cave bear lineage reached at least temporarily to the Altai Mountains, 7000 km east of their known centre of distribution. These results suggest that the diversity of cave bears was greater than previously believed, and that they could survive in a much wider range of ecological conditions than previously assumed. They also agree with recent studies on other extinct and extant species, such as wolves, hyenas and steppe bison, which have also revealed higher genetic and ecological diversity in Pleistocene populations than previously known.

A phylogeny and timescale for the living genera of kangaroos and kin (Macropodiformes:Marsupialia) based on nuclear DNA sequences

Kangaroos and kin (Macropodiformes) are the most conspicuous elements of the Australasian marsupial fauna. The approximately 70 living species can be divided into three families: (1) Hypsiprymnodontidae (the musky rat kangaroo); (2) Potoroidae (potoroos and bettongs); and (3) Macropodidae (larger kangaroos, wallabies, banded hare wallaby and pademelons). Here we examine macropodiform relationships using protein-coding portions of the ApoB, BRCA1, IRBP, Rag1 and vWF genes via maximum parsimony, maximum likelihood and Bayesian methods. We estimate times of divergence using two different relaxed molecular clock methods to present a timescale for macropodiform evolution and reconstruct ancestral states for grades of dental organisation. We find robust support for a basal split between Hypsiprymnodontidae and the other macropodiforms, potoroid monophyly and macropodid monophyly, with Lagostrophus as the sister-taxon to all other macropodids. Our divergence estimates suggest that kangaroos diverged from Phalangeroidea in the early Eocene, that crown-group Macropodiformes originated in the late Eocene or early Oligocene and that the potoroid–macropodid split occurred in the late Oligocene or early Miocene followed by rapid cladogenesis within these families 5 to 15 million years ago. These divergence estimates coincide with major geological and ecological changes in Australia. Ancestral state reconstructions for grades of dental organisation suggest that the grazer grade evolved independently on two different occasions within Macropodidae.

Genetic variation and phylogeography of free-living mouse species (genus Mus) in the Balkans and the Middle East

This work presents a study of the distribution and pattern of variation throughout the ranges of three free-living mouse species of the genus Mus-M. macedonicus, M. spicilegus, and a M. cypriacus - based on sequencing of two segments of the mitochondrial DNA (mtDNA) control region. The study shows a similar level of variability in the three species and suggests their recent population expansion. The highest proportion of variation is found within populations indicating low genetic structuring. Phylogenetic analysis confirms the significant divergence of a mitochondrial lineage of M. macedonicus from Israel, recently described as a new subspecies, M. macedonicus spretoides. Conversely, no genetic hiatus is revealed between European and Asian populations of M. macedonicus macedonicus. Although phylogenetic relationships among M. spicilegus populations could not be unravelled precisely, the results suggest a recent westward expansion of the species. The mtDNA divergence between M. macedonicus and M. spicilegus is 7.3%, suggesting their split between c. 700,000 and 1 million years ago. These dates correspond with a coalescent estimate about 720,000 years ago. On the other hand, M. cypriacus appeared almost twice as divergent from the former species (4.5%) as from the latter (8.8%) suggesting a divergence of c. 430,000-610,000 years ago (coalescent approximately 490,000 years ago) and 830,000-1.2 million years ago (coalescent approximately 780,000 years ago), respectively. Approximate times of population expansion have also been estimated for all taxa and groups of populations. Existence of several glacial refuges and various colonization scenarios are discussed; since all estimated divergence times fall within interglacial periods it seems that climatic oscillations did not play a crucial role in the evolution of the three species.

Inference of demographic history from genealogical trees using reversible jump Markov chain Monte Carlo

BACKGROUND

Coalescent theory is a general framework to model genetic variation in a population. Specifically, it allows inference about population parameters from sampled DNA sequences. However, most currently employed variants of coalescent theory only consider very simple demographic scenarios of population size changes, such as exponential growth.

RESULTS

Here we develop a coalescent approach that allows Bayesian non-parametric estimation of the demographic history using genealogies reconstructed from sampled DNA sequences. In this framework inference and model selection is done using reversible jump Markov chain Monte Carlo (MCMC). This method is computationally efficient and overcomes the limitations of related non-parametric approaches such as the skyline plot. We validate the approach using simulated data. Subsequently, we reanalyze HIV-1 sequence data from Central Africa and Hepatitis C virus (HCV) data from Egypt.

CONCLUSIONS

The new method provides a Bayesian procedure for non-parametric estimation of the demographic history. By construction it additionally provides confidence limits and may be used jointly with other MCMC-based coalescent approaches.