Estimation of primate speciation dates using local molecular clocks

The sequence, structure and evolutionary features of HOTAIR in mammals

Background

An increasing number of long noncoding RNAs (lncRNAs) have been identified recently. Different from all the others that function in cis to regulate local gene expression, the newly identified HOTAIR is located between HoxC11 and HoxC12 in the human genome and regulates HoxD expression in multiple tissues. Like the well-characterised lncRNA Xist, HOTAIR binds to polycomb proteins to methylate histones at multiple HoxD loci, but unlike Xist, many details of its structure and function, as well as the trans regulation, remain unclear. Moreover, HOTAIR is involved in the aberrant regulation of gene expression in cancer.

Results

To identify conserved domains in HOTAIR and study the phylogenetic distribution of this lncRNA, we searched the genomes of 10 mammalian and 3 non-mammalian vertebrates for matches to its 6 exons and the two conserved domains within the 1800 bp exon6 using Infernal. There was just one high-scoring hit for each mammal, but many low-scoring hits were found in both mammals and non-mammalian vertebrates. These hits and their flanking genes in four placental mammals and platypus were examined to determine whether HOTAIR contained elements shared by other lncRNAs. Several of the hits were within unknown transcripts or ncRNAs, many were within introns of, or antisense to, protein-coding genes, and conservation of the flanking genes was observed only between human and chimpanzee. Phylogenetic analysis revealed discrete evolutionary dynamics for orthologous sequences of HOTAIR exons. Exon1 at the 5' end and a domain in exon6 near the 3' end, which contain domains that bind to multiple proteins, have evolved faster in primates than in other mammals. Structures were predicted for exon1, two domains of exon6 and the full HOTAIR sequence. The sequence and structure of two fragments, in exon1 and the domain B of exon6 respectively, were identified to robustly occur in predicted structures of exon1, domain B of exon6 and the full HOTAIR in mammals.

Conclusions

HOTAIR exists in mammals, has poorly conserved sequences and considerably conserved structures, and has evolved faster than nearby HoxC genes. Exons of HOTAIR show distinct evolutionary features, and a 239 bp domain in the 1804 bp exon6 is especially conserved. These features, together with the absence of some exons and sequences in mouse, rat and kangaroo, suggest ab initio generation of HOTAIR in marsupials. Structure prediction identifies two fragments in the 5' end exon1 and the 3' end domain B of exon6, with sequence and structure invariably occurring in various predicted structures of exon1, the domain B of exon6 and the full HOTAIR.

A distance-based least-square method for dating speciation events

Distance-based phylogenetic methods are widely used in biomedical research. However, there has been little development of rigorous statistical methods and software for dating speciation and gene duplication events by using evolutionary distances. Here we present a simple, fast and accurate dating method based on the least-squares (LS) method that has already been widely used in molecular phylogenetic reconstruction. Dating methods with a global clock or two different local clocks are presented. Single or multiple fossil calibration points can be used, and multiple data sets can be integrated in a combined analysis. Variation of the estimated divergence time is estimated by resampling methods such as bootstrapping or jackknifing. Application of the method to dating the divergence time among seven ape species or among 35 mammalian species including major mammalian orders shows that the estimated divergence time with the LS criterion is nearly identical to those obtained by the likelihood method or Bayesian inference.

Genomic and evolutionary features of the SPI-1 type III secretion system that is present in Xanthomonas albilineans but is not essential for xylem colonization and symptom development of sugarcane leaf scald

Xanthomonas albilineans is the causal agent of sugarcane leaf scald. Interestingly, this bacterium, which is not known to be insect or animal associated, possesses a type III secretion system (T3SS) belonging to the injectisome family Salmonella pathogenicity island 1 (SPI-1). The T3SS SPI-1 of X. albilineans shares only low similarity with other available T3SS SPI-1 sequences. Screening of a collection of 128 plant-pathogenic bacteria revealed that this T3SS SPI-1 is present in only two species of Xanthomonas: X. albilineans and X. axonopodis pv. phaseoli. Inoculation of sugarcane with knockout mutants showed that this system is not required by X. albilineans to spread within xylem vessels and to cause disease symptoms. This result was confirmed by the absence of this T3SS SPI-1 in an X. albilineans strain isolated from diseased sugarcane. To investigate the importance of the T3SS SPI-1 during the life cycle of X. albilineans, we analyzed T3SS SPI-1 sequences from 11 strains spanning the genetic diversity of this species. No nonsense mutations or frameshifting indels were observed in any of these strains, suggesting that the T3SS SPI-1 system is maintained within the species X. albilineans. Evolutionary features of T3SS SPI-1 based on phylogenetic, recombination, and selection analyses are discussed in the context of the possible functional importance of T3SS SPI-1 in the ecology of X. albilineans.

Site-specific time heterogeneity of the substitution process and its impact on phylogenetic inference

Background

Model violations constitute the major limitation in inferring accurate phylogenies. Characterizing properties of the data that are not being correctly handled by current models is therefore of prime importance. One of the properties of protein evolution is the variation of the relative rate of substitutions across sites and over time, the latter is the phenomenon called heterotachy. Its effect on phylogenetic inference has recently obtained considerable attention, which led to the development of new models of sequence evolution. However, thus far focus has been on the quantitative heterogeneity of the evolutionary process, thereby overlooking more qualitative variations.

Results

We studied the importance of variation of the site-specific amino-acid substitution process over time and its possible impact on phylogenetic inference. We used the CAT model to define an infinite mixture of substitution processes characterized by equilibrium frequencies over the twenty amino acids, a useful proxy for qualitatively estimating the evolutionary process. Using two large datasets, we show that qualitative changes in site-specific substitution properties over time occurred significantly. To test whether this unaccounted qualitative variation can lead to an erroneous phylogenetic tree, we analyzed a concatenation of mitochondrial proteins in which Cnidaria and Porifera were erroneously grouped. The progressive removal of the sites with the most heterogeneous CAT profiles across clades led to the recovery of the monophyly of Eumetazoa (Cnidaria+Bilateria), suggesting that this heterogeneity can negatively influence phylogenetic inference.

Conclusion

The time-heterogeneity of the amino-acid replacement process is therefore an important evolutionary aspect that should be incorporated in future models of sequence change.

Against expectation: a short sequence with high signal elucidates cone snail phylogeny

A short (259 nucleotide) conserved intronic sequence (CIS) is surprisingly informative for delineating deep phylogenetic relationships in cone snails. Conus species previously have been assigned to clades based on the evidence from mitochondrial 12S and 16S rRNA gene sequences (1129 bp). Despite their length, these genes lack the phylogenetic information necessary to resolve the relationships among the clades. Here we show that the relationships can be inferred from just 46 sites in the very short CIS sequence (a portion of "intron 9" of the γ-glutamyl carboxylase gene). This is counterintuitive because in short sequences sampling error (noise) often drowns out phylogenetic signal. The intron 9 CIS is rich in synapomorphies that define the divergence patterns among eight clades of worm- and fish-hunting Conus, and it contains almost no homoplasy. Parsimony, maximum likelihood and Bayesian analyses of the combined sequences (mt rRNA+CIS) confirm most of the relationships among 23 Conus sequences. This phylogeny implies that fish-hunting behavior evolved at least twice during the history of Conus-once among New World species and independently in the Indo-Pacific clades.

Post-glacial expansion into the Paleozoic Plateau: evidence of an Ozarkian refugium for the Ozark minnow Notropis nubilus (Teleostei: Cypriniformes)

Genetic variation was examined within the Ozark minnow Notropis nubilus using complete mtDNA cytochrome b gene sequences from 160 individuals representing 30 localities to test hypotheses on the origin of the distribution. Phylogenetic analyses revealed three strongly supported clades of haplotypes consistent with geographic distributions: a clade from the Western Ozarks, a clade from the Southern Ozarks and a clade from the Northern Ozarks and upper Mississippi River basin. The estimated mean ages of these clades indicated that they diverged during pre-Illinoian glacial cycles extending from the late Pliocene into the early Pleistocene. Results of demographic analyses based on coalescent approaches supported the hypothesis that the Paleozoic Plateau was not a refugium for N. nubilus during periodic glacial advances. There is evidence of a genetic signature of northern expansion into the Paleozoic Plateau from a Southern Ozarkian refugium. Populations expanded out of drainages in the Northern Ozarks into the Paleozoic Plateau during the late Pleistocene. Subsequently, the two regions were isolated due to the recent extirpation of intervening populations caused by the loss of suitable habitat.

A comparison of global, gene-specific, and relaxed clock methods in a comparative genomics framework: dating the polyploid history of soybean (Glycine max)

It is widely recognized that many genes and lineages do not adhere to a molecular clock, yet molecular clocks are commonly used to date divergences in comparative genomic studies. We test the application of a molecular clock across genes and lineages in a phylogenetic framework utilizing 12 genes linked in a 1-Mb region on chromosome 13 of soybean (Glycine max); homoeologous copies of these genes formed by polyploidy in Glycine; and orthologous copies in G. tomentella, Phaseolus vulgaris, and Medicago truncatula. We compare divergence dates estimated by two methods each in three frameworks: a global molecular clock with a single rate across genes and lineages using full and approximate likelihood methods based on synonymous substitutions, a gene-specific clock assuming rate constancy over lineages but allowing a different rate for each gene, and a relaxed molecular clock where rates may vary across genes and lineages estimated under penalized likelihood and Bayesian inference. We use the cumulative variance across genes as a means of quantifying precision. Our results suggest that divergence dating methods produce results that are correlated, but that older nodes are more variable and more difficult to estimate with precision and accuracy. We also find that models incorporating less rate heterogeneity estimate older dates of divergence than more complex models, as node age increases. A mixed model nested analysis of variance testing the effects of framework, method, and gene found that framework had a significant effect on the divergence date estimates but that most variation among dates is due to variation among genes, suggesting a need to further characterize and understand the evolutionary phenomena underlying rate variation within genomes, among genes, and across lineages.

Bayesian random local clocks, or one rate to rule them all

Background

Relaxed molecular clock models allow divergence time dating and "relaxed phylogenetic" inference, in which a time tree is estimated in the face of unequal rates across lineages. We present a new method for relaxing the assumption of a strict molecular clock using Markov chain Monte Carlo to implement Bayesian modeling averaging over random local molecular clocks. The new method approaches the problem of rate variation among lineages by proposing a series of local molecular clocks, each extending over a subregion of the full phylogeny. Each branch in a phylogeny (subtending a clade) is a possible location for a change of rate from one local clock to a new one. Thus, including both the global molecular clock and the unconstrained model results, there are a total of 2^2n-2 possible rate models available for averaging with 1, 2, ..., 2n - 2 different rate categories.

Results

We propose an efficient method to sample this model space while simultaneously estimating the phylogeny. The new method conveniently allows a direct test of the strict molecular clock, in which one rate rules them all, against a large array of alternative local molecular clock models. We illustrate the method's utility on three example data sets involving mammal, primate and influenza evolution. Finally, we explore methods to visualize the complex posterior distribution that results from inference under such models.

Conclusions

The examples suggest that large sequence datasets may only require a small number of local molecular clocks to reconcile their branch lengths with a time scale. All of the analyses described here are implemented in the open access software package BEAST 1.5.4 (http://beast-mcmc.googlecode.com/).

Genome networks root the tree of life between prokaryotic domains

Eukaryotes arose from prokaryotes, hence the root in the tree of life resides among the prokaryotic domains. The position of the root is still debated, although pinpointing it would aid our understanding of the early evolution of life. Because prokaryote evolution was long viewed as a tree-like process of lineage bifurcations, efforts to identify the most ancient microbial lineage split have traditionally focused on positioning a root on a phylogenetic tree constructed from one or several genes. Such studies have delivered widely conflicting results on the position of the root, this being mainly due to methodological problems inherent to deep gene phylogeny and the workings of lateral gene transfer among prokaryotes over evolutionary time. Here, we report the position of the root determined with whole genome data using network-based procedures that take into account both gene presence or absence and the level of sequence similarity among all individual gene families that are shared across genomes. On the basis of 562,321 protein-coding gene families distributed across 191 genomes, we find that the deepest divide in the prokaryotic world is interdomain, that is, separating the archaebacteria from the eubacteria. This result resonates with some older views but conflicts with the results of most studies over the last decade that have addressed the issue. In particular, several studies have suggested that the molecular distinctness of archaebacteria is not evidence for their antiquity relative to eubacteria but instead stems from some kind of inherently elevated rate of archaebacterial sequence change. Here, we specifically test for such a rate elevation across all prokaryotic lineages through the analysis of all possible quartets among eight genes duplicated in all prokaryotes, hence the last common ancestor thereof. The results show that neither the archaebacteria as a group nor the eubacteria as a group harbor evidence for elevated evolutionary rates in the sampled genes, either in the recent evolutionary past or in their common ancestor. The interdomain prokaryotic position of the root is thus not attributable to lineage-specific rate variation.

Use of phylogenetics in the molecular epidemiology and evolutionary studies of viral infections

Since DNA sequencing techniques first became available almost 30 years ago, the amount of nucleic acid sequence data has increased enormously. Phylogenetics, which is widely applied to compare and analyze such data, is particularly useful for the analysis of genes from rapidly evolving viruses. It has been used extensively to describe the molecular epidemiology and transmission of the human immunodeficiency virus (HIV), the origins and subsequent evolution of the severe acute respiratory syndrome (SARS)-associated coronavirus (SCoV), and, more recently, the evolving epidemiology of avian influenza as well as seasonal and pandemic human influenza viruses. Recent advances in phylogenetic methods can infer more in-depth information about the patterns of virus emergence, adding to the conventional approaches in viral epidemiology. Examples of this information include estimations (with confidence limits) of the actual time of the origin of a new viral strain or its emergence in a new species, viral recombination and reassortment events, the rate of population size change in a viral epidemic, and how the virus spreads and evolves within a specific population and geographical region. Such sequence-derived information obtained from the phylogenetic tree can assist in the design and implementation of public health and therapeutic interventions. However, application of many of these advanced phylogenetic methods are currently limited to specialized phylogeneticists and statisticians, mainly because of their mathematical basis and their dependence on the use of a large number of computer programs. This review attempts to bridge this gap by presenting conceptual, technical, and practical aspects of applying phylogenetic methods in studies of influenza, HIV, and SCoV. It aims to provide, with minimal mathematics and statistics, a practical overview of how phylogenetic methods can be incorporated into virological studies by clinical and laboratory specialists.

Comment on the paleobiology and classification of Ardipithecus ramidus

White and colleagues (Research Articles, 2 October 2009, pp. 64-106 and www.sciencemag.org/ardipithecus) reported Ardipithecus ramidus as an exclusive member of the human lineage post-African ape divergence. However, their analysis of shared-derived characters provides insufficient evidence of an ancestor-descendant relationship and exclusivity to the hominid lineage. Molecular and anatomical studies rather suggest that Ar. ramidus predates the human/African ape divergence.

Branch length estimation and divergence dating: estimates of error in Bayesian and maximum likelihood frameworks

BACKGROUND

Estimates of divergence dates between species improve our understanding of processes ranging from nucleotide substitution to speciation. Such estimates are frequently based on molecular genetic differences between species; therefore, they rely on accurate estimates of the number of such differences (i.e. substitutions per site, measured as branch length on phylogenies). We used simulations to determine the effects of dataset size, branch length heterogeneity, branch depth, and analytical framework on branch length estimation across a range of branch lengths. We then reanalyzed an empirical dataset for plethodontid salamanders to determine how inaccurate branch length estimation can affect estimates of divergence dates.

RESULTS

The accuracy of branch length estimation varied with branch length, dataset size (both number of taxa and sites), branch length heterogeneity, branch depth, dataset complexity, and analytical framework. For simple phylogenies analyzed in a Bayesian framework, branches were increasingly underestimated as branch length increased; in a maximum likelihood framework, longer branch lengths were somewhat overestimated. Longer datasets improved estimates in both frameworks; however, when the number of taxa was increased, estimation accuracy for deeper branches was less than for tip branches. Increasing the complexity of the dataset produced more misestimated branches in a Bayesian framework; however, in an ML framework, more branches were estimated more accurately. Using ML branch length estimates to re-estimate plethodontid salamander divergence dates generally resulted in an increase in the estimated age of older nodes and a decrease in the estimated age of younger nodes.

CONCLUSIONS

Branch lengths are misestimated in both statistical frameworks for simulations of simple datasets. However, for complex datasets, length estimates are quite accurate in ML (even for short datasets), whereas few branches are estimated accurately in a Bayesian framework. Our reanalysis of empirical data demonstrates the magnitude of effects of Bayesian branch length misestimation on divergence date estimates. Because the length of branches for empirical datasets can be estimated most reliably in an ML framework when branches are <1 substitution/site and datasets are > or =1 kb, we suggest that divergence date estimates using datasets, branch lengths, and/or analytical techniques that fall outside of these parameters should be interpreted with caution.

The evolutionary rate dynamically tracks changes in HIV-1 epidemics: application of a simple method for optimizing the evolutionary rate in phylogenetic trees with longitudinal data

Large-sequence datasets provide an opportunity to investigate the dynamics of pathogen epidemics. Thus, a fast method to estimate the evolutionary rate from large and numerous phylogenetic trees becomes necessary. Based on minimizing tip height variances, we optimize the root in a given phylogenetic tree to estimate the most homogenous evolutionary rate between samples from at least two different time points. Simulations showed that the method had no bias in the estimation of evolutionary rates and that it was robust to tree rooting and topological errors. We show that the evolutionary rates of HIV-1 subtype B and C epidemics have changed over time, with the rate of evolution inversely correlated to the rate of virus spread. For subtype B, the evolutionary rate slowed down and tracked the start of the HAART era in 1996. Subtype C in Ethiopia showed an increase in the evolutionary rate when the prevalence increase markedly slowed down in 1995. Thus, we show that the evolutionary rate of HIV-1 on the population level dynamically tracks epidemic events.

Age at first reproduction explains rate variation in the strepsirrhine molecular clock

Although the molecular clock hypothesis posits that the rate of molecular change is constant over time, there is evidence that rates vary among lineages. Some of the strongest evidence for variable molecular rates comes from the primates; e.g., the "hominoid slowdown." These rate differences are hypothesized to correlate with certain species attributes, such as generation time and body size. Here, we examine rates of molecular change in the strepsirrhine suborder of primates and test whether body size or age at first reproduction (a proxy for generation time) explains patterns of rate variation better than a null model where the molecular clock is independent of these factors. To examine these models, we analyzed DNA sequences from four pairs of recently diverged strepsirrhine sister taxa to estimate molecular rates by using sign tests, likelihood ratio tests, and regression analyses. Our analysis does not support a model where body weight or age at first reproduction strongly influences rates of molecular evolution across mitochondrial and nuclear sites. Instead, our analysis supports a model where age at first reproduction influences neutral evolution in the nuclear genome. This study supports the generation time hypothesis for rate variation in the nuclear molecular clock. Molecular clock variation due to generation time may help to resolve the discordance between molecular and paleontological estimates for divergence date estimates in primate evolution.

Estimating the phylogeny and divergence times of primates using a supermatrix approach

Background

The primates are among the most broadly studied mammalian orders, with the published literature containing extensive analyses of their behavior, physiology, genetics and ecology. The importance of this group in medical and biological research is well appreciated, and explains the numerous molecular phylogenies that have been proposed for most primate families and genera. Composite estimates for the entire order have been infrequently attempted, with the last phylogenetic reconstruction spanning the full range of primate evolutionary relationships having been conducted over a decade ago.

Results

To estimate the structure and tempo of primate evolutionary history, we employed Bayesian phylogenetic methods to analyze data supermatrices comprising 7 mitochondrial genes (6,138 nucleotides) from 219 species across 67 genera and 3 nuclear genes (2,157 nucleotides) from 26 genera. Many taxa were only partially represented, with an average of 3.95 and 5.43 mitochondrial genes per species and per genus, respectively, and 2.23 nuclear genes per genus. Our analyses of mitochondrial DNA place Tarsiiformes as the sister group of Strepsirrhini. Within Haplorrhini, we find support for the primary divergence of Pitheciidae in Platyrrhini, and our results suggest a sister grouping of African and non-African colobines within Colobinae and of Cercopithecini and Papionini within Cercopthecinae. Date estimates for nodes within each family and genus are presented, with estimates for key splits including: Strepsirrhini-Haplorrhini 64 million years ago (MYA), Lemuriformes-Lorisiformes 52 MYA, Platyrrhini-Catarrhini 43 MYA and Cercopithecoidea-Hominoidea 29 MYA.

Conclusion

We present an up-to-date, comprehensive estimate of the structure and tempo of primate evolutionary history. Although considerable gaps remain in our knowledge of the primate phylogeny, increased data sampling, particularly from nuclear loci, will be able to provide further resolution.

Phylogenetic perspectives on the epidemiology and origins of SARS and SARS-like coronaviruses

Severe Acute Respiratory Syndrome (SARS) is a respiratory disease caused by a zoonotic coronavirus (CoV) named SARS-CoV (SCoV), which rapidly swept the globe after its emergence in rural China during late 2002. The origins of SCoV have been mysterious and controversial, until the recent discovery of SARS-like CoV (SLCoV) in bats and the proposal of bats as the natural reservior of the Coronaviridae family. In this article, we focused on discussing how phylogenetics contributed to our understanding towards the emergence and transmission of SCoV. We first reviewed the epidemiology of SCoV from a phylogenetic perspective and discussed the controversies over its phylogenetic origins. Then, we summarized the phylogenetic findings in relation to its zoonotic origins and the proposed inter-species viral transmission events. Finally, we also discussed how the discoveries of SCoV and SLCoV expanded our knowledge on the evolution of the Coronaviridae family as well as its implications on the possible future re-emergence of SCoV.

Quantitative prediction of molecular clock and ka/ks at short timescales

Recent empirical studies of taxa including humans, fish, and birds have shown elevated rates of molecular evolution between species that diverged recently. Using the Moran model, we calculate expected divergence as a function of time. Our findings suggest that the observed phenomenon of elevated rates at short timescales is consistent with standard population genetics theory. The apparent acceleration of the molecular clock at short timescales can be explained by segregating polymorphisms present at the time of the ancestral population, both neutral and slightly deleterious, and not newly arising slightly deleterious mutations as has been previously hypothesized. Our work also suggests that the duration of the rate elevation depends on the effective population size, providing a method to correct time estimates of recent divergence events. Our model concords with estimates of divergence obtained from African cichlid fish and humans. As an additional application of our model, we calculate that K(a)/K(s) is elevated within a population before decaying slowly to its long-term value. Similar to the molecular clock, the duration and magnitude of K(a)/K(s) elevation depend on the effective population size. Unlike the molecular clock, however, K(a)/K(s) elevation is caused by newly arising slightly deleterious mutations. This elevation, although not as severe in magnitude as had been previously predicted in models neglecting ancestral polymorphism, persists slightly longer.

Phylogeny of Banana Streak Virus reveals recent and repetitive endogenization in the genome of its banana host (Musa sp.)

Banana streak virus (BSV) is a plant dsDNA pararetrovirus (family Caulimoviridae, genus badnavirus). Although integration is not an essential step in the BSV replication cycle, the nuclear genome of banana (Musa sp.) contains BSV endogenous pararetrovirus sequences (BSV EPRVs). Some BSV EPRVs are infectious by reconstituting a functional viral genome. Recent studies revealed a large molecular diversity of episomal BSV viruses (i.e., nonintegrated) while others focused on BSV EPRV sequences only. In this study, the evolutionary history of badnavirus integration in banana was inferred from phylogenetic relationships between BSV and BSV EPRVs. The relative evolution rates and selective pressures (d(N)/d(S) ratio) were also compared between endogenous and episomal viral sequences. At least 27 recent independent integration events occurred after the divergence of three banana species, indicating that viral integration is a recent and frequent phenomenon. Relaxation of selective pressure on badnaviral sequences that experienced neutral evolution after integration in the plant genome was recorded. Additionally, a significant decrease (35%) in the EPRV evolution rate was observed compared to BSV, reflecting the difference in the evolution rate between episomal dsDNA viruses and plant genome. The comparison of our results with the evolution rate of the Musa genome and other reverse-transcribing viruses suggests that EPRVs play an active role in episomal BSV diversity and evolution.

Primate numts and reticulate evolution of capped and golden leaf monkeys (Primates: Colobinae)

A recent phylogenetic study of langurs and leaf monkeys of South Asia suggested a reticulate evolution of capped and golden leaf monkeys through ancient hybridization between Semnopithecus and Trachypithecus .To test this hybridization scenario, I analysed nuclear copies of the mitochondrial cytochrome b gene (numts) from capped,golden and Phayre's leaf monkeys. These numts were aligned with mitochondrial cytochrome b sequences of various species belonging to the genera Semnopithecus and Trachypithecus .In the phylogenetic tree derived from this alignment,the numts fell into three distinct clades (A,B and C) suggesting three independent integration events.Clade A was basal to Semnopithecus, and clades B and C were basal to Trachypithecus. Among the numts in clades A and C were sequences derived from species not represented in their respective sister mitochondrial groups.This unusual placement of certain numts is taken as additional support for the hybridization scenario. Based on the molecular dating of these integration events,hybridization is estimated to have occurred around 7.1 to 3.4 million years ago.Capped and golden leaf monkeys might have to be assigned to a new genus to reconcile their unique evolutionary history. Additionally, northeast India appears to be a 'hot spot' for lineages that might have evolved through reticulate evolution.

TIMING DEEP DIVERGENCE EVENTS IN CALCAREOUS DINOFLAGELLATES(1)

Based on morphological and molecular data, calcareous dinoflagellates (Thoracosphaeraceae, Peridiniales) are a monophyletic group comprising the three major clades Ensiculifera/Pentapharsodinium, Thoracosphaera/Pfiesteria, and Scrippsiella sensu lato. We used stratigraphically well-documented first occurrences of particular archeopyle types to constrain relaxed Bayesian molecular clocks applied to nuclear rRNA sequences of 18 representatives of the three main clades. By comparing divergence estimates obtained in differently calibrated clocks with first stratigraphic occurrences of taxa not themselves used as constraints, we identified plausible divergence times for several subclades of calcareous dinoflagellates. The initial diversification of extant calcareous dinoflagellates probably took place in the Late Jurassic, with the three main clades all established by the Cretaceous. The two mesoepicystal operculum types observed in calcareous dinoflagellates probably evolved independently from simple apical archeopyles. Based on our taxon sample, the K/T boundary had relatively little effect on the diversity of the group, with several lineages dating to before 65 mya (million years ago). The first stratigraphic occurrences of key taxa, such as Thoracosphaera and Calciodinellum (not themselves used as constraints), are in agreement with the molecular time estimates. Conflicts that involve "Calciodinellum"levantinum, Leonella, Pentapharsodinium, Pernambugia, and the Scrippsiella trochoidea species complex may be due to inaccurate assignment of fossils because of high morphological homoplasy and insufficient knowledge of the extant diversity of calcareous dinoflagellates.

Reconstructing phylogenies and phenotypes: a molecular view of human evolution

This review broadly summarizes how molecular biology has contributed to our understanding of human evolution. Molecular anthropology began in the 1960s with immunological comparisons indicating that African apes and humans were closely related and, indeed, shared a common ancestor as recently as 5 million years ago. Although initially dismissed, this finding has proven robust and numerous lines of molecular evidence now firmly place the human-ape divergence at 4-8 Ma. Resolving the trichotomy among humans, chimpanzees and gorillas took a few more decades. Despite the readily apparent physical similarities shared by African apes to the exclusion of modern humans (body hair, knuckle-walking, thin tooth enamel), the molecular support for a human-chimpanzee clade is now overwhelming. More recently, whole genome sequencing and gene mapping have shifted the focus of molecular anthropology from phylogenetic analyses to phenotypic reconstruction and functional genomics. We are starting to identify the genetic basis of the morphological, physiological and behavioural traits that distinguish modern humans from apes and apes from other primates. Most notably, recent comparative genomic analyses strongly indicate that the marked differences between modern humans and chimpanzees are likely due more to changes in gene regulation than to modifications of the genes themselves, an idea first proposed over 30 years ago. Almost weekly, press releases describe newly identified genes and regulatory elements that seem to have undergone strong positive selection along the human lineage. Loci involved in speech (e.g. FOXP2), brain development (e.g. ASPM), and skull musculature (e.g. MYH16) have been of particular interest, but some surprising candidate loci (e.g. those involved in auditory capabilities) have emerged as well. Exciting new research avenues, such as the Neanderthal Genome Project, promise that molecular analyses will continue to provide novel insights about our evolution. Ultimately, however, these molecular findings can only be understood in light of data from field sites, morphology labs, and museum collections. Indeed, molecular anthropology depends on these sources for calibrating molecular clocks and placing genetic data within the context of key morphological and ecological transitions in human evolution.

Birth-death prior on phylogeny and speed dating

Background

In recent years there has been a trend of leaving the strict molecular clock in order to infer dating of speciations and other evolutionary events. Explicit modeling of substitution rates and divergence times makes formulation of informative prior distributions for branch lengths possible. Models with birth-death priors on tree branching and auto-correlated or iid substitution rates among lineages have been proposed, enabling simultaneous inference of substitution rates and divergence times. This problem has, however, mainly been analysed in the Markov chain Monte Carlo (MCMC) framework, an approach requiring computation times of hours or days when applied to large phylogenies.

Results

We demonstrate that a hill-climbing maximum a posteriori (MAP) adaptation of the MCMC scheme results in considerable gain in computational efficiency. We demonstrate also that a novel dynamic programming (DP) algorithm for branch length factorization, useful both in the hill-climbing and in the MCMC setting, further reduces computation time. For the problem of inferring rates and times parameters on a fixed tree, we perform simulations, comparisons between hill-climbing and MCMC on a plant rbcL gene dataset, and dating analysis on an animal mtDNA dataset, showing that our methodology enables efficient, highly accurate analysis of very large trees. Datasets requiring a computation time of several days with MCMC can with our MAP algorithm be accurately analysed in less than a minute. From the results of our example analyses, we conclude that our methodology generally avoids getting trapped early in local optima. For the cases where this nevertheless can be a problem, for instance when we in addition to the parameters also infer the tree topology, we show that the problem can be evaded by using a simulated-annealing like (SAL) method in which we favour tree swaps early in the inference while biasing our focus towards rate and time parameter changes later on.

Conclusion

Our contribution leaves the field open for fast and accurate dating analysis of nucleotide sequence data. Modeling branch substitutions rates and divergence times separately allows us to include birth-death priors on the times without the assumption of a molecular clock. The methodology is easily adapted to take data from fossil records into account and it can be used together with a broad range of rate and substitution models.

Evidence of the recombinant origin of a bat severe acute respiratory syndrome (SARS)-like coronavirus and its implications on the direct ancestor of SARS coronavirus

Bats have been identified as the natural reservoir of severe acute respiratory syndrome (SARS)-like and SARS coronaviruses (SLCoV and SCoV). However, previous studies suggested that none of the currently sampled bat SLCoVs is the descendant of the direct ancestor of SCoV, based on their relatively distant phylogenetic relationship. In this study, evidence of the recombinant origin of the genome of a bat SLCoV is demonstrated. We identified a potential recombination breakpoint immediately after the consensus intergenic sequence between open reading frame 1 and the S coding region, suggesting the replication intermediates may participate in the recombination event, as previously speculated for other CoVs. Phylogenetic analysis of its parental regions suggests the presence of an uncharacterized SLCoV lineage that is phylogenetically closer to SCoVs than any of the currently sampled bat SLCoVs. Using various Bayesian molecular-clock models, interspecies transfer of this SLCoV lineage from bats to the amplifying host (e.g., civets) was estimated to have happened a median of 4.08 years before the SARS outbreak. Based on this relatively short window period, we speculate that this uncharacterized SLCoV lineage may contain the direct ancestor of SCoV. This study sheds light on the possible host bat species of the direct ancestor of SCoV, providing valuable information on the scope and focus of surveillance for the origin of SCoV.

Type 2 NF1 deletions are highly unusual by virtue of the absence of nonallelic homologous recombination hotspots and an apparent preference for female mitotic recombination

Approximately 5% of patients with neurofibromatosis type 1 (NF1) exhibit gross deletions that encompass the NF1 gene and its flanking regions. The breakpoints of the common 1.4-Mb (type 1) deletions are located within low-copy repeats (NF1-REPs) and cluster within a 3.4-kb hotspot of nonallelic homologous recombination (NAHR). Here, we present the first comprehensive breakpoint analysis of type 2 deletions, which are a second type of recurring NF1 gene deletion. Type 2 deletions span 1.2 Mb and are characterized by breakpoints located within the SUZ12 gene and its pseudogene, which closely flank the NF1-REPs. Breakpoint analysis of 13 independent type 2 deletions did not reveal any obvious hotspots of NAHR. However, an overrepresentation of polypyrimidine/polypurine tracts and triplex-forming sequences was noted in the breakpoint regions that could have facilitated NAHR. Intriguingly, all 13 type 2 deletions identified so far are characterized by somatic mosaicism, which indicates a positional preference for mitotic NAHR within the NF1 gene region. Indeed, whereas interchromosomal meiotic NAHR occurs between the NF1-REPs giving rise to type 1 deletions, NAHR during mitosis appears to occur intrachromosomally between the SUZ12 gene and its pseudogene, thereby generating type 2 deletions. Such a clear distinction between the preferred sites of mitotic versus meiotic NAHR is unprecedented in any other genomic disorder induced by the local genomic architecture. Additionally, 12 of the 13 mosaic type 2 deletions were found in females. The marked female preponderance among mosaic type 2 deletions contrasts with the equal sex distribution noted for type 1 and/or atypical NF1 deletions. Although an influence of chromatin structure was strongly suspected, no sex-specific differences in the methylation pattern exhibited by the SUZ12 gene were apparent that could explain the higher rate of mitotic recombination in females.

Evolutionary rate variation at multiple levels of biological organization in plant mitochondrial DNA

We examined patterns of mitochondrial polymorphism and divergence in the angiosperm genus Silene and found substantial variation in evolutionary rates among species and among lineages within species. Moreover, we found corresponding differences in the amount of polymorphism within species. We argue that, along with our earlier findings of rate variation among genes, these patterns of rate heterogeneity at multiple phylogenetic scales are most likely explained by differences in underlying mutation rates. In contrast, no rate variation was detected in nuclear or chloroplast loci. We conclude that mutation rate heterogeneity is a characteristic of plant mitochondrial sequence evolution at multiple biological scales and may be a crucial determinant of how much polymorphism is maintained within species. These dramatic patterns of variation raise intriguing questions about the mechanisms driving and maintaining mutation rate heterogeneity in plant mitochondrial genomes. Additionally, they should alter our interpretation of many common phylogenetic and population genetic analyses.

BEAST: Bayesian evolutionary analysis by sampling trees

BACKGROUND

The evolutionary analysis of molecular sequence variation is a statistical enterprise. This is reflected in the increased use of probabilistic models for phylogenetic inference, multiple sequence alignment, and molecular population genetics. Here we present BEAST: a fast, flexible software architecture for Bayesian analysis of molecular sequences related by an evolutionary tree. A large number of popular stochastic models of sequence evolution are provided and tree-based models suitable for both within- and between-species sequence data are implemented.

RESULTS

BEAST version 1.4.6 consists of 81000 lines of Java source code, 779 classes and 81 packages. It provides models for DNA and protein sequence evolution, highly parametric coalescent analysis, relaxed clock phylogenetics, non-contemporaneous sequence data, statistical alignment and a wide range of options for prior distributions. BEAST source code is object-oriented, modular in design and freely available at http://beast-mcmc.googlecode.com/ under the GNU LGPL license.

CONCLUSION

BEAST is a powerful and flexible evolutionary analysis package for molecular sequence variation. It also provides a resource for the further development of new models and statistical methods of evolutionary analysis.

BEAST: Bayesian evolutionary analysis by sampling trees

BACKGROUND

The evolutionary analysis of molecular sequence variation is a statistical enterprise. This is reflected in the increased use of probabilistic models for phylogenetic inference, multiple sequence alignment, and molecular population genetics. Here we present BEAST: a fast, flexible software architecture for Bayesian analysis of molecular sequences related by an evolutionary tree. A large number of popular stochastic models of sequence evolution are provided and tree-based models suitable for both within- and between-species sequence data are implemented.

RESULTS

BEAST version 1.4.6 consists of 81000 lines of Java source code, 779 classes and 81 packages. It provides models for DNA and protein sequence evolution, highly parametric coalescent analysis, relaxed clock phylogenetics, non-contemporaneous sequence data, statistical alignment and a wide range of options for prior distributions. BEAST source code is object-oriented, modular in design and freely available at http://beast-mcmc.googlecode.com/ under the GNU LGPL license.

CONCLUSION

BEAST is a powerful and flexible evolutionary analysis package for molecular sequence variation. It also provides a resource for the further development of new models and statistical methods of evolutionary analysis.

Evolutionary analyses of European H1N2 swine influenza A virus by placing timestamps on the multiple reassortment events

A novel H1N2 swine influenza A virus emerged in Europe since 1994. Previous phylogenetic analyses revealed that its genome segments were derived from H1N1 human virus, H3N2 human virus and avian-like H1N1/H3N2 swine virus, indicating the possibility of multiple reassortments events. However, dates of these reassortment events have not been investigated systematically. In this study, we used both global and local molecular clock concepts in a maximum likelihood framework to extrapolate the times of origins of the genome segments in European H1N2 swine viruses, and deduced that novel neuraminidase, hemagglutinin and other internal protein genes were introduced to the European H1N2 lineage at the 1970s, early 1980s and late 1980s, respectively through reassortments. Furthermore, in light of the evolutionary timescale reconstructed for the H1N2 viruses, we argue that further reassortments, in addition to those responsible for the introductions of novel genome segments, might have also occurred among the viruses prior to the outbreaks arose in United Kingdom at 1994. Our results confirm that the viral genes of various origins have stably maintained in swine population for many years before the multiple genetic reassortant was detected. Our evolutionary analyses also suggested that the HA and NA genes evolved in a significantly higher rate of synonymous substitutions after they were introduced from human to swine and established the European H1N2 swine lineage.

Estimating Divergence Times in Large Phylogenetic Trees

A new method, PATHd8, for estimating ultrametric trees from trees with edge (branch) lengths proportional to the number of substitutions is proposed. The method allows for an arbitrary number of reference nodes for time calibration, each defined either as absolute age, minimum age, or maximum age, and the tree need not be fully resolved. The method is based on estimating node ages by mean path lengths from the node to the leaves but correcting for deviations from a molecular clock suggested by reference nodes. As opposed to most existing methods allowing substitution rate variation, the new method smoothes substitution rates locally, rather than simultaneously over the whole tree, thus allowing for analysis of very large trees. The performance of PATHd8 is compared with other frequently used methods for estimating divergence times. In analyses of three separate data sets, PATHd8 gives similar divergence times to other methods, the largest difference being between crown group ages, where unconstrained nodes get younger ages when analyzed with PATHd8. Overall, chronograms obtained from other methods appear smoother, whereas PATHd8 preserves more of the heterogeneity seen in the original edge lengths. Divergence times are most evenly spread over the chronograms obtained from the Bayesian implementation and the clock-based Langley-Fitch method, and these two methods produce very similar ages for most nodes. Evaluations of PATHd8 using simulated data suggest that PATHd8 is slightly less precise compared with penalized likelihood, but it gives more sensible answers for extreme data sets. A clear advantage with PATHd8 is that it is more or less instantaneous even with trees having several thousand leaves, whereas other programs often run into problems when analyzing trees with hundreds of leaves. PATHd8 is implemented in freely available software.

Are Cape floral clades the same age? Contemporaneous origins of two lineages in the genistoids s.l. (Fabaceae)

The hypothesis that the elements of the modern species-rich flora of the Cape Floristic Region (CFR), South Africa, originated more or less simultaneously at the Miocene/Pliocene boundary, in response to the development of a mediterranean climate, has been challenged by numerous molecular dating estimates of Cape floral clades. These studies reveal a more gradual emergence, with the oldest clades originating in the Eocene, but others appearing later, some as recently as the Pliocene. That there are factors which might affect the dates recovered, such as choice of calibration point, analysis method, sampling density and the delimitation of Cape floral clades, suggests a need for further critical evaluation of the age estimates presented to date. In this study, the dates of origin of two Cape floral clades (the legume Crotalarieae p.p. and Podalyrieae) are estimated, constrained by a shared calibration point in a single analysis using an rDNA ITS phylogeny in which 633 taxa are sampled. The results indicate that these two clades arose contemporaneously 44-46 mya, not at the Miocene/Pliocene boundary as had been previously supposed. The contemporaneous origin of these Cape floral clades suggests that additional more inclusive analyses are needed before rejecting the hypothesis that a single environmental trigger explains the establishment of Cape floral clades.

Inferring speciation times under an episodic molecular clock

We extend our recently developed Markov chain Monte Carlo algorithm for Bayesian estimation of species divergence times to allow variable evolutionary rates among lineages. The method can use heterogeneous data from multiple gene loci and accommodate multiple fossil calibrations. Uncertainties in fossil calibrations are described using flexible statistical distributions. The prior for divergence times for nodes lacking fossil calibrations is specified by use of a birth-death process with species sampling. The prior for lineage-specific substitution rates is specified using either a model with autocorrelated rates among adjacent lineages (based on a geometric Brownian motion model of rate drift) or a model with independent rates among lineages specified by a log-normal probability distribution. We develop an infinite-sites theory, which predicts that when the amount of sequence data approaches infinity, the width of the posterior credibility interval and the posterior mean of divergence times form a perfect linear relationship, with the slope indicating uncertainties in time estimates that cannot be reduced by sequence data alone. Simulations are used to study the influence of among-lineage rate variation and the number of loci sampled on the uncertainty of divergence time estimates. The analysis suggests that posterior time estimates typically involve considerable uncertainties even with an infinite amount of sequence data, and that the reliability and precision of fossil calibrations are critically important to divergence time estimation. We apply our new algorithms to two empirical data sets and compare the results with those obtained in previous Bayesian and likelihood analyses. The results demonstrate the utility of our new algorithms.

Dating phylogenies with hybrid local molecular clocks

Background

Because rates of evolution and species divergence times cannot be estimated directly from molecular data, all current dating methods require that specific assumptions be made before inferring any divergence time. These assumptions typically bear either on rates of molecular evolution (molecular clock hypothesis, local clocks models) or on both rates and times (penalized likelihood, Bayesian methods). However, most of these assumptions can affect estimated dates, oftentimes because they underestimate large amounts of rate change.

Principal Findings

A significant modification to a recently proposed ad hoc rate-smoothing algorithm is described, in which local molecular clocks are automatically placed on a phylogeny. This modification makes use of hybrid approaches that borrow from recent theoretical developments in microarray data analysis. An ad hoc integration of phylogenetic uncertainty under these local clock models is also described. The performance and accuracy of the new methods are evaluated by reanalyzing three published data sets.

Conclusions

It is shown that the new maximum likelihood hybrid methods can perform better than penalized likelihood and almost as well as uncorrelated Bayesian models. However, the new methods still tend to underestimate the actual amount of rate change. This work demonstrates the difficulty of estimating divergence times using local molecular clocks.

Primate Origins: Implications of a Cretaceous Ancestry

It has long been accepted that the adaptive radiation of modern placental mammals, like that of modern birds, did not begin until after the Cretaceous/Tertiary (K/T) boundary 65 million years (Ma) ago, following the extinction of the dinosaurs. The first undoubted fossil relatives of modern primates appear in the record 55 Ma ago. However, in agreement with evidence from molecular phylogenies calibrated with dates from denser parts of the fossil record, a statistical analysis of the primate record allowing for major gaps now indicates a Cretaceous origin of euprimates 80-90 Ma ago. If this interpretation is correct, primates overlapped with dinosaurs by some 20 Ma prior to the K/T boundary, and the initial radiation of primates was probably truncated as part of the major extinction event that occurred at the end of the Cretaceous. Following a review of evidence for an early origin of primates, implications of this are discussed with respect to the likely ancestral condition for primates, including a southern continental area of origin and moderately large body size. The known early Tertiary primates are re-interpreted as northern continental offshoots of a 'second wave' of primate evolution.

Body mass and temperature influence rates of mitochondrial DNA evolution in North American cyprinid fish

The mass-specific metabolic rate hypothesis of Gillooly and others predicts that DNA mutation and substitution rates are a function of body mass and temperature. We tested this hypothesis with sequence divergences estimated from mtDNA cytochrome b sequences of 54 taxa of cyprinid fish. Branch lengths estimated from a likelihood tree were compared with metabolic rates calculated from body mass and environmental temperatures experienced by those taxa. The problem of unknown age estimates of lineage splitting was avoided by comparing estimated amounts of metabolic activity along phyletic lines leading to pairs of modern taxa from their most recent common ancestor with sequence divergences along those same pairs of phyletic lines. There were significantly more pairs for which the phyletic line with greater genetic change also had the higher metabolic activity, when compared to the prediction of a hypothesis that body mass and temperature are not related to substitution rate.