Biomedical subjects

Alexei Drummond

Publications and source records attributed to Alexei Drummond.

6 recordsLinked to original sources

PubMed2006

Phylogenetic analysis reveals a correlation between the expansion of very virulent infectious bursal disease virus and reassortment of its genome segment B.

Infectious bursal disease virus (IBDV) is a birnavirus causing immunosuppressive disease in chickens. Emergence of the very virulent form of IBDV (vvIBDV) in the late 1980s dramatically changed the epidemiology of the disease. In this study, we investigated the phylogenetic origins of its genome segments and estimated the time of emergence of their most recent common ancestors. Moreover, with recently developed coalescence techniques, we reconstructed the past population dynamics of vvIBDV and timed the onset of its expansion to the late 1980s. Our analysis suggests that genome segment A of vvIBDV emerged at least 20 years before its expansion, which argues against the hypothesis that mutation of genome segment A is the major contributing factor in the emergence and expansion of vvIBDV. Alternatively, the phylogeny of genome segment B suggests a possible reassortment event estimated to have taken place around the mid-1980s, which seems to coincide with its expansion within approximately 5 years. We therefore hypothesize that the reassortment of genome segment B initiated vvIBDV expansion in the late 1980s, possibly by enhancing the virulence of the virus synergistically with its existing genome segment A. This report reveals the possible mechanisms leading to the emergence and expansion of vvIBDV, which would certainly provide insights into the scope of surveillance and prevention efforts regarding the disease.

Animals↗

PubMed2005

Bayesian coestimation of phylogeny and sequence alignment.

BACKGROUND: Two central problems in computational biology are the determination of the alignment and phylogeny of a set of biological sequences. The traditional approach to this problem is to first build a multiple alignment of these sequences, followed by a phylogenetic reconstruction step based on this multiple alignment. However, alignment and phylogenetic inference are fundamentally interdependent, and ignoring this fact leads to biased and overconfident estimations. Whether the main interest be in sequence alignment or phylogeny, a major goal of computational biology is the co-estimation of both. RESULTS: We developed a fully Bayesian Markov chain Monte Carlo method for coestimating phylogeny and sequence alignment, under the Thorne-Kishino-Felsenstein model of substitution and single nucleotide insertion-deletion (indel) events. In our earlier work, we introduced a novel and efficient algorithm, termed the "indel peeling algorithm", which includes indels as phylogenetically informative evolutionary events, and resembles Felsenstein's peeling algorithm for substitutions on a phylogenetic tree. For a fixed alignment, our extension analytically integrates out both substitution and indel events within a proper statistical model, without the need for data augmentation at internal tree nodes, allowing for efficient sampling of tree topologies and edge lengths. To additionally sample multiple alignments, we here introduce an efficient partial Metropolized independence sampler for alignments, and combine these two algorithms into a fully Bayesian co-estimation procedure for the alignment and phylogeny problem. Our approach results in estimates for the posterior distribution of evolutionary rate parameters, for the maximum a-posteriori (MAP) phylogenetic tree, and for the posterior decoding alignment. Estimates for the evolutionary tree and multiple alignment are augmented with confidence estimates for each node height and alignment column. Our results indicate that the patterns in reliability broadly correspond to structural features of the proteins, and thus provides biologically meaningful information which is not existent in the usual point-estimate of the alignment. Our methods can handle input data of moderate size (10-20 protein sequences, each 100-200 bp), which we analyzed overnight on a standard 2 GHz personal computer. CONCLUSION: Joint analysis of multiple sequence alignment, evolutionary trees and additional evolutionary parameters can be now done within a single coherent statistical framework.

Algorithms↗

PubMed2003

Extreme reversed sexual size dimorphism in the extinct New Zealand moa Dinornis.

The ratite moa (Aves; Dinornithiformes) were massive graviportal browsers weighing up to 250 kg (ref. 1) that dominated the New Zealand biota until their extinction approximately 500 yr ago. Despite an extensive Quaternary fossil record, moa taxonomy remains problematic and currently 11 species are recognized. Three Dinornis species were found throughout New Zealand and differed markedly in size (1-2 m height at back) and mass (from approximately 34 to 242 kg). Surprisingly, ancient mitochondrial DNA sequences show that the three species were genetically indistinguishable within each island, but formed separate North and South Island clades. Here we show, using the first sex-linked nuclear sequences from an extinct species, that on each island the three morphological forms actually represent just one species, whose size varied markedly according to sex and habitat. The largest females in this example of extreme reversed sexual size dimorphism were about 280% the weight and 150% the height of the largest males, which is unprecedented among birds and terrestrial mammals. The combination of molecular and palaeontological data highlights the difficulties of analysing extinct groups, even those with detailed fossil records.

Animals↗

PubMed2003

Inferring evolutionary rates using serially sampled sequences from several populations.

The estimation of evolutionary rates from serially sampled sequences has recently been the focus of several studies. In this paper, we extend these analyzes to allow the estimation of a joint rate of substitution, omega, from several evolving populations from which serial samples are drawn. In the case of viruses evolving in different hosts, therapy may halt replication and therefore the accumulation of substitutions in the population. In such cases, it may be that only a proportion, p, of subjects are nonresponders who have viral populations that continue to evolve. We develop two likelihood-based procedures to jointly estimate p and omega, and empirical Bayes' tests of whether an individual should be classified as a responder or nonresponder. An example data set comprising HIV-1 partial envelope sequences from six patients on highly active antiretroviral therapy is analyzed.

Algorithms↗

PubMed2003

Inference of viral evolutionary rates from molecular sequences.

The processes of mutation and nucleotide substitution contribute to the observed variability in virulence, transmission and persistence of viral pathogens. Since most viruses evolve many times faster than their human hosts, we are in the unusual position of being able to measure these processes directly by comparing viral genes that have been isolated and sequenced at different points in time. The analysis of such data requires the use of specific statistical methods that take into account the shared ancestry of the sequences and the randomness inherent in the process of nucleotide substitution. In this paper we describe the various statistical methods for estimating evolutionary rates, which can be classified into three general approaches: linear regression, maximum likelihood, and Bayesian inference. We discuss the advantages and shortcomings of each approach and illustrate their use through the analysis of two example viruses; human immunodeficiency virus type 1 and dengue virus serotype 4. Reliable estimates of viral substitution rates have many important applications in population genetics and phylogenetics, including dating evolutionary events and divergence times, estimating demographic parameters such as population size and generation time, and investigating the effect of natural selection on molecular evolution.

Base Sequence↗

PubMed2003

Epidemiology, genetic diversity, and evolution of endemic feline immunodeficiency virus in a population of wild cougars.

Within the large body of research on retroviruses, the distribution and evolution of endemic retroviruses in natural host populations have so far received little attention. In this study, the epidemiology, genetic diversity, and molecular evolution of feline immunodeficiency virus specific to cougars (FIVpco) was examined using blood samples collected over several years from a free-ranging cougar population in the western United States. The virus prevalence was 58% in this population (n = 52) and increased significantly with host age. Based on phylogenetic analysis of fragments of envelope (env) and polymerase (pol) genes, two genetically distinct lineages of FIVpco were found to cooccur in the population but not in the same individuals. Within each of the virus lineages, geographically nearby isolates formed monophyletic clusters of closely related viruses. Sequence diversity for env within a host rarely exceeded 1%, and the evolution of this gene was dominated by purifying selection. For both pol and env, our data indicate mean rates of molecular evolution of 1 to 3% per 10 years. These results support the premise that FIVpco is well adapted to its cougar host and provide a basis for comparing lentivirus evolution in endemic and epidemic infections in natural hosts.

Animals↗