Empirical and hierarchical Bayesian estimation of ancestral states

Several methods have been proposed to infer the states at the ancestral nodes on a phylogeny. These methods assume a specific tree and set of branch lengths when estimating the ancestral character state. Inferences of the ancestral states, then, are conditioned on the tree and branch lengths being true. We develop a hierarchical Bayes method for inferring the ancestral states on a tree. The method integrates over uncertainty in the tree, branch lengths, and substitution model parameters by using Markov chain Monte Carlo. We compare the hierarchical Bayes inferences of ancestral states with inferences of ancestral states made under the assumption that a specific tree is correct. We find that the methods are correlated, but that accommodating uncertainty in parameters of the phylogenetic model can make inferences of ancestral states even more uncertain than they would be in an empirical Bayes analysis.     

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 Bayesian estimation approach for the mortality in a stage-structured demographic model

Journal of Mathematical Biology - Control interventions in sustainable pest management schemes are set according to the phenology and the population abundance of the pests. This information can be...     

SNP-based heritability estimation using a Bayesian approach

Heritability is a central element in quantitative genetics. New molecular markers to assess genetic variance and heritability are continually under development. The availability of molecular single nucleotide polymorphism (SNP) markers can be applied for estimation of variance components and heritability on population, where relationship information is unknown. In this study, we evaluated the capabilities of two Bayesian genomic models to estimate heritability in simulated populations. The populations comprised different family structures of either no or a limited number of relatives, a single quantitative trait, and with one of two densities of SNP markers. All individuals were both genotyped and phenotyped. Results illustrated that the two models were capable of estimating heritability, when true heritability was 0.15 or higher and populations had a sample size of 400 or higher. For heritabilities of 0.05, all models had difficulties in estimating the true heritability. The two Bayesian models were compared with a restricted maximum likelihood (REML) approach using a genomic relationship matrix. The comparison showed that the Bayesian approaches performed equally well as the REML approach. Differences in family structure were in general not found to influence the estimation of the heritability. For the sample sizes used in this study, a 10-fold increase of SNP density did not improve precision estimates compared with set-ups with a less dense distribution of SNPs. The methods used in this study showed that it was possible to estimate heritabilities on the basis of SNPs in animals with direct measurements. This conclusion is valuable in cases when quantitative traits are either difficult or expensive to measure.     

A bayesian analysis of metazoan mitochondrial genome arrangements

Genome arrangements are a potentially powerful source of information to infer evolutionary relationships among distantly related taxa. Mitochondrial genome arrangements may be especially informative about metazoan evolutionary relationships because (1) nearly all animals have the same set of definitively homologous mitochondrial genes, (2) mitochondrial genome rearrangement events are rare relative to changes in sequences, and (3) the number of possible mitochondrial genome arrangements is huge, making convergent evolution of genome arrangements appear highly unlikely. In previous studies, phylogenetic evidence in genome arrangement data is nearly always used in a qualitative fashion-the support in favor of clades with similar or identical genome arrangements is considered to be quite strong, but is not quantified. The purpose of this article is to quantify the uncertainty among the relationships of metazoan phyla on the basis of mitochondrial genome arrangements while incorporating prior knowledge of the monophyly of various groups from other sources. The work we present here differs from our previous work in the statistics literature in that (1) we incorporate prior information on classifications of metazoans at the phylum level, (2) we describe several advances in our computational approach, and (3) we analyze a much larger data set (87 taxa) that consists of each unique, complete mitochondrial genome arrangement with a full complement of 37 genes that were present in the NCBI (National Center for Biotechnology Information) database at a recent date. In addition, we analyze a subset of 28 of these 87 taxa for which the non-tRNA mitochondrial genomes are unique where the assumption of our inversion-only model of rearrangement is more plausible. We present summaries of Bayesian posterior distributions of tree topology on the basis of these two data sets.     

Bayesian inference of ancestral recombination graphs

We present a novel algorithm, implemented in the software ARGinfer, for probabilistic inference of the Ancestral Recombination Graph under the Coalescent with Recombination. Our Markov Chain Monte Carlo algorithm takes advantage of the Succinct Tree Sequence data structure that has allowed great advances in simulation and point estimation, but not yet probabilistic inference. Unlike previous methods, which employ the Sequentially Markov Coalescent approximation, ARGinfer uses the Coalescent with Recombination, allowing more accurate inference of key evolutionary parameters. We show using simulations that ARGinfer can accurately estimate many properties of the evolutionary history of the sample, including the topology and branch lengths of the genealogical tree at each sequence site, and the times and locations of mutation and recombination events. ARGinfer approximates posterior probability distributions for these and other quantities, providing interpretable assessments of uncertainty that we show to be well calibrated. ARGinfer is currently limited to tens of DNA sequences of several hundreds of kilobases, but has scope for further computational improvements to increase its applicability.     

Inferred chromosome morphology of the ancestral genome ofTriticum

The lengths of the A, B, and D genomes of common wheat,Triticum aestivum, were measured from the karyotype. Relative to the B genome, standardized as length 1.000, the lengths of the A and D genomes were 0.835 and 0.722, respectively. The lengths of the chromosome arms in the A and D genomes were then multiplied by the appropriate constants so that the total lengths of each genome also equalled 1.000. These calculations revealed that homoeologous chromosomes in wheat, with a few exceptions, have similar sizes and arm ratios. The arm lengths of the three homoeologues in each homoeologous group were then averaged. These average chromosomes turned out to be remarkably similar, in size and arm ratio, to their homoeologues in the E genome ofElytrigia elongata. This evidence and data on cross-compatibility and morphological characteristics suggested that the genusTriticum is a result of adaptive radiation from the perennial genusElytrigia, specifically from the complex of species possessing the E genome or one closely related to it.     

The environmental efficiency of Wanjiang demonstration area: A Bayesian estimation approach

With the implementation of the strategic policy and the construction of Wanjiang demonstration area recently, the government cares more about the natural environment. To protect the environment of this demonstration area, a quantitative analysis of environmental efficiency and its influencing factors is needed. In this paper, we measure the environmental efficiency of the demonstration area, and then analyze the total factor productivity of the area by Malmquist productivity index through data envelopment analysis (DEA) approach. Through the index and its decomposition, we are able to reveal the changing trend of the environmental efficiency. Moreover, we apply Bayesian estimation approach to analyze the influencing factors of the efficiency and discuss the relationship between these factors and the environmental efficiency. Results show that the efficiency of Wanjiang demonstration area bares a significant difference among constituent cities. Each area should devise proper environmental policy according to its particular circumstance.     

The tendency to recreate ancestral CG dinucleotides in the human genome

Background  

The CG dinucleotides are known to be deficient in the human genome, due to a high mutation rate from 5-methylated CG to TG and its complementary pair CA. Meanwhile, many cellular functions rely on these CG dinucleotides, such as gene expression controlled by cytosine methylation status. Thus, CG dinucleotides that provide essential functional substrates should be retained in genomes. How these two conflicting processes regarding the fate of CG dinucleotides - i.e., high mutation rate destroying CG dinucleotides, vs. functional processes that require their preservation remains an unsolved question.     

A primer-based approach to genome walking

A genome walking strategy based on annealing and ligation of single-stranded DNA primers to 3′ overhangs following restriction endonuclease digestion was developed. A set of primers contains 4 nucleotides at the 3′ end that are complementary to overhangs formed by restriction endonucleasesApaI;PstI;SacI andSphI. Following ligation, 5′ end overhangs are formed on the DNA, which serves as sites for the adaptor primers and nested primers for PCR amplification in combination with the gene-specific primers. This strategy was verified by the amplification of up to 4 kb of a potato leafroll virus full-length infectious clone. The procedure could be adopted to target any upstream and downstream regions flanking known sequences within the plant genome.     

A Bayesian approach to the evolution of social learning

There has been much interest in understanding the evolution of social learning. Investigators have tried to understand when natural selection will favor individuals who imitate others, how imitators should deal with the fact that available models may exhibit different behaviors, and how social and individual learning should interact. In all of this work, social learning and individual learning have been treated as alternative, conceptually distinct processes. Here we present a Bayesian model in which both individual and social learning arise from a single inferential process. Individuals use Bayesian inference to combine social and nonsocial cues about the current state of the environment. This model indicates that natural selection favors individuals who place heavy weight on social cues when the environment changes slowly or when its state cannot be well predicted using nonsocial cues. It also indicates that a conformist bias should be a universal aspect of social learning.     

Evolution of gastropod mitochondrial genome arrangements

Background  

Gastropod mitochondrial genomes exhibit an unusually great variety of gene orders compared to other metazoan mitochondrial genome such as e.g those of vertebrates. Hence, gastropod mitochondrial genomes constitute a good model system to study patterns, rates, and mechanisms of mitochondrial genome rearrangement. However, this kind of evolutionary comparative analysis requires a robust phylogenetic framework of the group under study, which has been elusive so far for gastropods in spite of the efforts carried out during the last two decades. Here, we report the complete nucleotide sequence of five mitochondrial genomes of gastropods (Pyramidella dolabrata, Ascobulla fragilis, Siphonaria pectinata, Onchidella celtica, and Myosotella myosotis), and we analyze them together with another ten complete mitochondrial genomes of gastropods currently available in molecular databases in order to reconstruct the phylogenetic relationships among the main lineages of gastropods.