Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods ( BayesMultiState) is available from the authors.

A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Evolutionary history of vegetative reproduction in Porpidia s. l. (lichen-forming ascomycota)

The evolutionary history of gains and losses of vegetative reproductive propagules (soredia) in Porpidia s.l., a group of lichen-forming ascomycetes, was clarified using Bayesian Markov chain Monte Carlo (MCMC) approaches to monophyly tests and a combined MCMC and maximum likelihood approach to ancestral character state reconstructions. The MCMC framework provided confidence estimates for the reconstructions of relationships and ancestral character states, which formed the basis for tests of evolutionary hypotheses. Monophyly tests rejected all hypotheses that predicted any clustering of reproductive modes in extant taxa. In addition, a nearest-neighbor statistic could not reject the hypothesis that the vegetative reproductive mode is randomly distributed throughout the group. These results show that transitions between presence and absence of the vegetative reproductive mode within Porpidia s.l. occurred several times and independently of each other. Likelihood reconstructions of ancestral character states at selected nodes suggest that - contrary to previous thought - the ancestor to Porpidia s.l. already possessed the vegetative reproductive mode. Furthermore, transition rates are reconstructed asymmetrically with the vegetative reproductive mode being gained at a much lower rate than it is lost. A cautious note has to be added, because a simulation study showed that the ancestral character state reconstructions were highly dependent on taxon sampling. However, our central conclusions, particularly the higher rate of change from vegetative reproductive mode present to absent than vice versa within Porpidia s.l., were found to be broadly independent of taxon sampling. [Ancestral character state reconstructions; Ascomycota, Bayesian inference; hypothesis testing; likelihood; MCMC; Porpidia; reproductive systems]

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

The effect of doubled CO2 and model basic state biases on the monsoon-ENSO system. I: Mean response and interannual variability

The impact of doubled CO2 concentration on the Asian summer monsoon is studied using a coupled ocean-atmosphere model. Both the mean seasonal precipitation and interannual monsoon variability are found to increase in the future climate scenario presented. Systematic biases in current climate simulations of the coupled system prevent accurate representation of the monsoon-ENSO teleconnection, of prime importance for seasonal prediction and for determining monsoon interannual variability. By applying seasonally varying heat flux adjustments to the tropical Pacific and Indian Ocean surface in the future climate simulation, some assessment can be made of the impact of systematic model biases on future climate predictions. In simulations where the flux adjustments are implemented, the response to climate change is magnified, with the suggestion that systematic biases may be masking the true impact of increased greenhouse gas forcing. The teleconnection between ENSO and the Asian summer monsoon remains robust in the future climate, although the Indo-Pacific takes on more of a biennial character for long periods of the flux-adjusted simulation. Assessing the teleconnection across interdecadal timescales shows wide variations in its amplitude, despite the absence of external forcing. This suggests that recent changes in the observed record cannot be distinguished from internal variations and as such are not necessarily related to climate change.

Detecting the node-density artifact in phylogeny reconstruction

The node-density effect is an artifact of phylogeny reconstruction that can cause branch lengths to be underestimated in areas of the tree with fewer taxa. Webster, Payne, and Pagel (2003, Science 301:478) introduced a statistical procedure (the "delta" test) to detect this artifact, and here we report the results of computer simulations that examine the test's performance. In a sample of 50,000 random data sets, we find that the delta test detects the artifact in 94.4% of cases in which it is present. When the artifact is not present (n = 10,000 simulated data sets) the test showed a type I error rate of approximately 1.69%, incorrectly reporting the artifact in 169 data sets. Three measures of tree shape or "balance" failed to predict the size of the node-density effect. This may reflect the relative homogeneity of our randomly generated topologies, but emphasizes that nearly any topology can suffer from the artifact, the effect not being confined only to highly unevenly sampled or otherwise imbalanced trees. The ability to screen phylogenies for the node-density artifact is important for phylogenetic inference and for researchers using phylogenetic trees to infer evolutionary processes, including their use in molecular clock dating. [Delta test; molecular clock; molecular evolution; node-density effect; phylogenetic reconstruction; speciation; simulation.]

Limpets break Dollo's law

A new molecular phylogeny of the limpet molluscs (Calyptraeidae) reveals that coiled shells have independently re-evolved at least once in this family, which is a violation of Dollo's Law that complex ancestral states, once lost, are never reacquired. Reacquisition of the coiled ancestral state is remarkable in that uncoiled shells have been the most recent ancestral state for 20 million-100 million years. Adult coiling might have reevolved by the mechanism of prolonging the period during which genes for coiling are expressed in larvae. This and other developmental mechanisms could provide general routes for maintaining the potential to produce traits lost in distant ancestors.

Cross-species amplification and phylogenetic reconstruction using fragaria microsatellite primers

Twenty-eight microsatellite primer pairs developed from Fragaria vesca ‘Rügen’ were applied to sixteen accessions representing eight diploid Fragaria species. The number of alleles generated, the power of discrimination and the percentage of accessions where no PCR product could be amplified were calculated for each locus for the thirteen non-F. vesca accessions. A phylogeny was then generated for the species accessions sampled, using the presence or absence of alleles at the polymorphic loci as character states. Despite the problems inherent in phylogeny reconstruction from microsatellite data, the phylogeny showed some congruence with a previously published phylogeny of Fragaria, based on nucleotide sequence data. However, relationships inferred from microsatellite allele data were relatively unresolved and poorly supported. The genetic basis of allelic polymorphisms at specific loci was investigated through direct sequencing of the PCR products amplified by three primer pairs. The potential utility of sequence data generated from microsatellite loci in evolutionary studies of closely related species groups is briefly explored.

Simultaneous bridge-localized and mixed-valence character in diruthenium radical cations featuring diethynylaromatic bridging ligands

A series of bimetallic ruthenium complexes [{Ru(dppe)Cp*}2(μ-C≡CArC≡C)] featuring diethynylaromatic bridging ligands (Ar = 1,4-phenylene, 1,4-naphthylene, 9,10-anthrylene) have been prepared and some representative molecular structures determined. A combination of UV–vis–NIR and IR spectroelectrochemical methods and density functional theory (DFT) have been used to demonstrate that one-electron oxidation of compounds [{Ru(dppe)Cp*}2(μ-C≡CArC≡C)](HC≡CArC≡CH = 1,4-diethynylbenzene; 1,4-diethynyl-2,5-dimethoxybenzene; 1,4-diethynylnaphthalene; 9,10-diethynylanthracene) yields solutions containing radical cations that exhibit characteristics of both oxidation of the diethynylaromatic portion of the bridge, and a mixed-valence state. The simultaneous population of bridge-oxidized and mixed-valence states is likely related to a number of factors, including orientation of the plane of the aromatic portion of the bridging ligand with respect to the metal d-orbitals of appropriate π-symmetry.

Partnership in the provision of recreational access: the reconstruction of the public interest in private land in England and Wales

This paper examines the evolution of public rights of access to private land in England and Wales. Since the Eighteenth Century the administration and protection of these rights has been though a form of public/private partnership in which the judiciary, while maintaining the dominance of private property, have safeguarded de facto public access by refusing consistently to punish simple trespass. While this situation has been modified, principally by post-World War II legislation, to allow for some formalisation of access arrangements and consequent compensation to landowners in areas of high recreational pressure and low legal accessibility, recent policy initiatives suggest that the balance of the partnership has now shifted in favour of landowners. In particular, the new access payment schemes, developed by the UK Government in response to the European Commission's Agri-Environment Regulations, identify the landowner as the beneficiary of the partnership, financed by tax revenue and justified on the spurious basis of improved 'access provision'. As such the State, as the former upholder of citizen rights, now assumes the duplicitous position of underwriting private property ownership through the commodification of access, while proclaiming a significant improvement in citizens' access rights.

Comparative biology of mammalian telomeres: hypotheses on ancestral states and the roles of telomeres in longevity determination

Progressive telomere shortening from cell division (replicative aging) provides a barrier for human tumor progression. This program is not conserved in laboratory mice, which have longer telomeres and constitutive telomerase. Wild species that do ⁄ do not use replicative aging have been reported, but the evolution of different phenotypes and a conceptual framework for understanding their uses of telomeres is lacking. We examined telomeres ⁄ telomerase in cultured cells from > 60 mammalian species to place different uses of telomeres in a broad mammalian context. Phylogeny-based statistical analysis reconstructed ancestral states. Our analysis suggested that the ancestral mammalian phenotype included short telomeres (< 20 kb, as we now see in humans) and repressed telomerase. We argue that the repressed telomerase was a response to a higher mutation load brought on by the evolution of homeothermy. With telomerase repressed, we then see the evolution of replicative aging. Telomere length inversely correlated with lifespan, while telomerase expression co-evolved with body size. Multiple independent times smaller, shorter-lived species changed to having longer telomeres and expressing telomerase. Trade-offs involving reducing the energetic ⁄ cellular costs of specific oxidative protection mechanisms (needed to protect < 20 kb telomeres in the absence oftelomerase) could explain this abandonment of replicative aging. These observations provide a conceptual framework for understanding different uses of telomeres in mammals, support a role for human-like telomeres in allowing longer lifespans to evolve, demonstrate the need to include telomere length in the analysis of comparative studies of oxidative protection in the biology of aging, and identify which mammals can be used as appropriate model organisms for the study of the role of telomeres in human cancer and aging. Key words: evolution of telomeres; immortalization; telomerase; replicative aging; senescence.

Estimating Lorenz's reference state in an ocean with a nonlinear equation of state for seawater

The study of the mechanical energy budget of the oceans using Lorenz available potential energy (APE) theory is based on knowledge of the adiabatically re-arranged Lorenz reference state of minimum potential energy. The compressible and nonlinear character of the equation of state for seawater has been thought to cause the reference state to be ill-defined, casting doubt on the usefulness of APE theory for investigating ocean energetics under realistic conditions. Using a method based on the volume frequency distribution of parcels as a function of temperature and salinity in the context of the seawater Boussinesq approximation, which we illustrate using climatological data, we show that compressibility effects are in fact minor. The reference state can be regarded as a well defined one-dimensional function of depth, which forms a surface in temperature, salinity and density space between the surface and the bottom of the ocean. For a very small proportion of water masses, this surface can be multivalued and water parcels can have up to two statically stable levels in the reference density profile, of which the shallowest is energetically more accessible. Classifying parcels from the surface to the bottom gives a different reference density profile than classifying in the opposite direction. However, this difference is negligible. We show that the reference state obtained by standard sorting methods is equivalent, though computationally more expensive, to the volume frequency distribution approach. The approach we present can be applied systematically and in a computationally efficient manner to investigate the APE budget of the ocean circulation using models or climatological data.