PubMed Health⌕ Search

Biomedical subjects

O G Pybus

Publications and source records attributed to O G Pybus.

17 recordsLinked to original sources

Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection.

The evolution of the human immunodeficiency virus (HIV-1) during chronic infection involves the rapid, continuous turnover of genetic diversity. However, the role of natural selection, relative to random genetic drift, in governing this process is unclear. We tested a stochastic model of genetic drift using partial envelope sequences sampled longitudinally in 28 infected children. In each case the Bayesian posterior (empirical) distribution of coalescent genealogies was estimated using Markov chain Monte Carlo methods. Posterior predictive simulation was then used to generate a null distribution of genealogies assuming neutrality, with the null and empirical distributions compared using four genealogy-based summary statistics sensitive to nonneutral evolution. Because both null and empirical distributions were generated within a coalescent framework, we were able to explicitly account for the confounding influence of demography. From the distribution of corrected P-values across patients, we conclude that empirical genealogies are more asymmetric than expected if evolution is driven by mutation and genetic drift only, with an excess of low-frequency polymorphisms in the population. This indicates that although drift may still play an important role, natural selection has a strong influence on the evolution of HIV-1 envelope. A negative relationship between effective population size and substitution rate indicates that as the efficacy of selection increases, a smaller proportion of mutations approach fixation in the population. This suggests the presence of deleterious mutations. We therefore conclude that intrahost HIV-1 evolution in envelope is dominated by purifying selection against low-frequency deleterious mutations that do not reach fixation.

Base Sequence↗

Do antiviral CD8+ T cells select hepatitis C virus escape mutants? Analysis in diverse epitopes targeted by human intrahepatic CD8+ T lymphocytes.

Hepatitis C virus (HCV) is a variable RNA virus that can readily establish persistent infection. Cellular immune responses are important in the early control of the virus. Evidence from animal models suggests that mutation in epitopes recognized by CD8+ T lymphocytes may play an important role in the establishment of persistence but in human persistent infection, equivalent evidence is lacking. We investigated this by analysing a unique resource: viruses from a set of chronically HCV-infected individuals in whom the CD8+ T-cell responses in liver had previously been accurately mapped. Virus was sequenced in seven individuals at 10 epitopes restricted by 10 human leucocyte antigen (HLA) molecules. Two main patterns emerged: in the majority of epitopes sequenced, no variation was seen. In three epitopes, mutations were identified which were compatible with immune escape as assessed using phylogenetic and/or functional studies. These data suggest that - even where specific intrahepatic T cells are detectable - many epitopes do not undergo mutation in chronic human infection. On the contrary, virus may escape from intrahepatic CD8+ T-cell responses in a 'patchy' manner in certain specific epitopes. Furthermore, longitudinal studies to identify the differences between 'selecting' and 'nonselecting' intrahepatic CD8+ T-cell responses are needed in HCV infection.

CD8-Positive T-Lymphocytes↗

Bayesian coalescent inference of past population dynamics from molecular sequences.

We introduce the Bayesian skyline plot, a new method for estimating past population dynamics through time from a sample of molecular sequences without dependence on a prespecified parametric model of demographic history. We describe a Markov chain Monte Carlo sampling procedure that efficiently samples a variant of the generalized skyline plot, given sequence data, and combines these plots to generate a posterior distribution of effective population size through time. We apply the Bayesian skyline plot to simulated data sets and show that it correctly reconstructs demographic history under canonical scenarios. Finally, we compare the Bayesian skyline plot model to previous coalescent approaches by analyzing two real data sets (hepatitis C virus in Egypt and mitochondrial DNA of Beringian bison) that have been previously investigated using alternative coalescent methods. In the bison analysis, we detect a severe but previously unrecognized bottleneck, estimated to have occurred 10,000 radiocarbon years ago, which coincides with both the earliest undisputed record of large numbers of humans in Alaska and the megafaunal extinctions in North America at the beginning of the Holocene.

Algorithms↗

Historical demography of Mullerian mimicry in the neotropical Heliconius butterflies.

We compare the historical demographies of two Müllerian comimetic butterfly species: Heliconius erato and Heliconius melpomene. These species show an extensive parallel geographic divergence in their aposematic wing phenotypes. Recent studies suggest that this coincident mosaic results from simultaneous demographic processes shaped by extrinsic forces over Pleistocene climate fluctuations. However, DNA sequence variation at two rapidly evolving unlinked nuclear loci, Mannose phosphate isomerase (Mpi) and Triose phosphate isomerase (Tpi), show that the comimetic species have quite different quaternary demographies. In H. erato, despite ongoing lineage sorting across the Andes, nuclear genealogical estimates showed little geographical structure, suggesting high historical gene flow. Coalescent-based demographic analysis revealed population growth since the Pliocene period. Although these patterns suggest vicariant population subdivision associated with the Andean orogeny, they are not consistent with hypotheses of Pleistocene population fragmentation facilitating allopatric wing phenotype radiation in H. erato. In contrast, nuclear genetic diversity, theta, in H. melpomene was reduced relative to its comimic and revealed three phylogeographical clades. The pattern of coalescent events within regional clades was most consistent with population growth in relatively isolated populations after a recent period of restricted population size. These different demographic histories suggest that the wing-pattern radiations were not coincident in the two species. Instead, larger effective population size (N(e)) in H. erato, together with profound population change in H. melpomene, supports an earlier hypothesis that H. erato diversified first as the model species of this remarkable mimetic association.

Alleles↗

The low evolutionary rate of human T-cell lymphotropic virus type-1 confirmed by analysis of vertical transmission chains.

The evolutionary rate of the human T-cell lymphotropic virus type-1 (HTLV-1) is considered to be very low, in strong contrast to the related human retrovirus HIV. However, current estimates of the HTLV-1 rate rely on the anthropological calibration of phylogenies using assumed dates of human migration events. To obtain an independent rate estimate, we analyzed two variable regions of the HTLV-1 genome (LTR and env) from eight infected families. Remarkable genetic stability was observed, as only two mutations in LTR (756 bp) and three mutations in env (522 bp) occurred within the 16 vertical transmission chains, including one ambiguous position in each region. The evolutionary rate in HTLV-1 was then calculated using a maximum-likelihood approach that used the highest and lowest possible times of HTLV-1 shared ancestry, given the known transmission histories. The rates for the LTR and env regions were 9.58 x 10(-8)-1.25 x 10(-5) and 7.84 x 10(-7) -2.33 x 10(-5)nucleotide substitutions per site per year, respectively. A more precise estimate was obtained for the combined LTR-env data set, which was 7.06 x 10(-7)-1.38 x 10(-5)substitutions per site per year. We also note an interesting correlation between the occurrence of mutations in HTLV-1 and the age of the individual infected.

Adolescent↗

Genetic structure of the deep-sea coral Lophelia pertusa in the northeast Atlantic revealed by microsatellites and internal transcribed spacer sequences.

The azooxanthellate scleractinian coral Lophelia pertusa has a near-cosmopolitan distribution, with a main depth distribution between 200 and 1000 m. In the northeast Atlantic it is the main framework-building species, forming deep-sea reefs in the bathyal zone on the continental margin, offshore banks and in Scandinavian fjords. Recent studies have shown that deep-sea reefs are associated with a highly diverse fauna. Such deep-sea communities are subject to increasing impact from deep-water fisheries, against a background of poor knowledge concerning these ecosystems, including the biology and population structure of L. pertusa. To resolve the population structure and to assess the dispersal potential of this deep-sea coral, specific microsatellites markers and ribosomal internal transcribed spacer (ITS) sequences ITS1 and ITS2 were used to investigate 10 different sampling sites, distributed along the European margin and in Scandinavian fjords. Both microsatellite and gene sequence data showed that L. pertusa should not be considered as one panmictic population in the northeast Atlantic but instead forms distinct, offshore and fjord populations. Results also suggest that, if some gene flow is occurring along the continental slope, the recruitment of sexually produced larvae is likely to be strongly local. The microsatellites showed significant levels of inbreeding and revealed that the level of genetic diversity and the contribution of asexual reproduction to the maintenance of the subpopulations were highly variable from site to site. These results are of major importance in the generation of a sustainable management strategy for these diversity-rich deep-sea ecosystems.

Animals↗

The epidemiology and iatrogenic transmission of hepatitis C virus in Egypt: a Bayesian coalescent approach.

Hepatitis C virus (HCV) is a leading cause of liver cancer and cirrhosis, and Egypt has possibly the highest HCV prevalence worldwide. In this article we use a newly developed Bayesian inference framework to estimate the transmission dynamics of HCV in Egypt from sampled viral gene sequences, and to predict the public health impact of the virus. Our results indicate that the effective number of HCV infections in Egypt underwent rapid exponential growth between 1930 and 1955. The timing and speed of this spread provides quantitative genetic evidence that the Egyptian HCV epidemic was initiated and propagated by extensive antischistosomiasis injection campaigns. Although our results show that HCV transmission has since decreased, we conclude that HCV is likely to remain prevalent in Egypt for several decades. Our combined population genetic and epidemiological analysis provides detailed estimates of historical changes in Egyptian HCV prevalence. Because our results are consistent with a demographic scenario specified a priori, they also provide an objective test of inference methods based on the coalescent process.

Bayes Theorem↗

Genetic diversity and models of viral evolution for the hepatitis C virus.

In this review we discuss the application of theoretical frameworks to the interpretation of viral gene sequence data, with particular reference to the hepatitis C virus (HCV). The increasing availability of such data means that it is now possible (and necessary) to proceed from simple qualitative models of viral evolution, to more quantitative frameworks based on statistical inference, notably population genetics and molecular phylogenetics. We argue that these approaches are invaluable tools to the virologist and are essential for understanding the dynamics of viral infection and the outcome of therapeutic strategies. We use several recent HCV data-sets to illustrate the methods.

Genetic Variation↗

New inferences from tree shape: numbers of missing taxa and population growth rates.

The relative positions of branching events in a phylogeny contain information about evolutionary and population dynamic processes. We provide new summary statistics of branching event times and describe how these statistics can be used to infer rates of species diversification from interspecies trees or rates of population growth from intraspecies trees. We also introduce a phylogenetic method for estimating the level of taxon sampling in a clade. Different evolutionary models and different sampling regimes can produce similar patterns of branching events, so it is important to consider explicitly the model assumptions involved when making evolutionary inferences. Results of an analysis of the phylogeny of the mosquito-borne flaviviruses suggest that there could be several thousand currently unidentified viruses in this clade.

Animals↗

GENIE: estimating demographic history from molecular phylogenies.

UNLABELLED: GENIE implements a statistical framework for inferring the demographic history of a population from phylogenies that have been reconstructed from sampled DNA sequences. The methods are based on population genetic models known collectively as coalescent theory. AVAILABILITY: GENIE is available from http://evolve.zoo.ox.ac.uk. All popular operating systems are supported.

Computer Simulation↗

Using human immunodeficiency virus type 1 sequences to infer historical features of the acquired immune deficiency syndrome epidemic and human immunodeficiency virus evolution.

In earlier work, human immunodeficiency virus type 1 (HIV-1) sequences were analysed to estimate the timing of the ancestral sequence of the main group of HIV-1, the virus that is responsible for the acquired immune deficiency syndrome pandemic, yielding a best estimate of 1931 (95% confidence interval of 1915-1941). That work will be briefly reviewed, outlining how phylogenetic tools were extended to incorporate improved evolutionary models, how the molecular clock model was adapted to incorporate variable periods of latency, and how the approach was validated by correctly estimating the timing of two historically documented dates. The advantages, limitations, and assumptions of the approach will be summarized, with particular consideration of the implications of branch length uncertainty and recombination. We have recently undertaken new phylogenetic analysis of an extremely diverse set of human immunodeficiency virus envelope sequences from the Democratic Republic of the Congo (the DRC, formerly Zaire). This analysis both corroborates and extends the conclusions of our original study. Coalescent methods were used to infer the demographic history of the HIV-1 epidemic in the DRC, and the results suggest an increase in the exponential growth rate of the infected population through time.

Acquired Immunodeficiency Syndrome↗

The epidemic behavior of the hepatitis C virus.

Hepatitis C virus (HCV) is a leading worldwide cause of liver disease. Here, we use a new model of HCV spread to investigate the epidemic behavior of the virus and to estimate its basic reproductive number from gene sequence data. We find significant differences in epidemic behavior among HCV subtypes and suggest that these differences are largely the result of subtype-specific transmission patterns. Our model builds a bridge between the disciplines of population genetics and mathematical epidemiology by using pathogen gene sequences to infer the population dynamic history of an infectious disease.

Endemic Diseases↗

Exploring the demographic history of DNA sequences using the generalized skyline plot.

We present an intuitive visual framework, the generalized skyline plot, to explore the demographic history of sampled DNA sequences. This approach is based on a genealogy inferred from the sequences and provides a nonparametric estimate of effective population size through time. In contrast to previous related procedures, the generalized skyline plot is more applicable to cases where the underlying tree is not fully resolved and the data is not highly variable. This is achieved by the grouping of adjacent coalescent intervals. We employ a small-sample Akaike information criterion to objectively choose the optimal grouping strategy. We investigate the performance of our approach using simulation and subsequently apply it to HIV-1 sequences from central Africa and mtDNA sequences from red pandas.

Africa↗

Testing macro-evolutionary models using incomplete molecular phylogenies.

Phylogenies reconstructed from gene sequences can be used to investigate the tempo and mode of species diversification. Here we develop and use new statistical methods to infer past patterns of speciation and extinction from molecular phylogenies. Specifically, we test the null hypothesis that per-lineage speciation and extinction rates have remained constant through time. Rejection of this hypothesis may provide evidence for evolutionary events such as adaptive radiations or key adaptations. In contrast to previous approaches, our methods are robust to incomplete taxon sampling and are conservative with respect to extinction. Using simulation we investigate, first, the adverse effects of failing to take incomplete sampling into account and, second, the power and reliability of our tests. When applied to published phylogenies our tests suggest that, in some cases, speciation rates have decreased through time.

Algorithms↗

An integrated framework for the inference of viral population history from reconstructed genealogies.

We describe a unified set of methods for the inference of demographic history using genealogies reconstructed from gene sequence data. We introduce the skyline plot, a graphical, nonparametric estimate of demographic history. We discuss both maximum-likelihood parameter estimation and demographic hypothesis testing. Simulations are carried out to investigate the statistical properties of maximum-likelihood estimates of demographic parameters. The simulations reveal that (i) the performance of exponential growth model estimates is determined by a simple function of the true parameter values and (ii) under some conditions, estimates from reconstructed trees perform as well as estimates from perfect trees. We apply our methods to HIV-1 sequence data and find strong evidence that subtypes A and B have different demographic histories. We also provide the first (albeit tentative) genetic evidence for a recent decrease in the growth rate of subtype B.

Africa South of the Sahara↗

The mid-depth method and HIV-1: a practical approach for testing hypotheses of viral epidemic history.

We introduce the mid-depth method, a practical approach for testing hypotheses of demographic history using genealogies reconstructed from sequence data. The relative positions of internal nodes within a genealogy contain information about past population dynamics. We explain how this information can be used to (1) test the null hypothesis of constant population size and (2) estimate the growth rate and current population size of an exponentially growing population. Simulation tests indicate that, as expected, estimates of exponential growth rates are sometimes biased. The mid-depth method is computationally rapid and does not require knowledge of the sample's mutation rate. However, it does assume that the reconstructed genealogy is correct and is therefore best suited to the analysis of variation-rich viral data sets. When applied to HIV-1 sequence data, the mid-depth method provides phylogenetic evidence of different exponential growth rates for subtypes A and B. We posit that this difference in growth rate reflects the different transmission routes and epidemiological histories of the two subtypes.

HIV-1↗