Search PubMed⌕ Search

Biomedical subjects

Simon D W Frost

Publications and source records attributed to Simon D W Frost.

At least 19 recordsLinked to original sources

Pneumococcal within-host diversity during colonization, transmission and treatment.

Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.

Child↗

Automated phylogenetic detection of recombination using a genetic algorithm.

The evolution of homologous sequences affected by recombination or gene conversion cannot be adequately explained by a single phylogenetic tree. Many tree-based methods for sequence analysis, for example, those used for detecting sites evolving nonneutrally, have been shown to fail if such phylogenetic incongruity is ignored. However, it may be possible to propose several phylogenies that can correctly model the evolution of nonrecombinant fragments. We propose a model-based framework that uses a genetic algorithm to search a multiple-sequence alignment for putative recombination break points, quantifies the level of support for their locations, and identifies sequences or clades involved in putative recombination events. The software implementation can be run quickly and efficiently in a distributed computing environment, and various components of the methods can be chosen for computational expediency or statistical rigor. We evaluate the performance of the new method on simulated alignments and on an array of published benchmark data sets. Finally, we demonstrate that prescreening alignments with our method allows one to analyze recombinant sequences for positive selection.

Algorithms↗

Adaptation to different human populations by HIV-1 revealed by codon-based analyses.

Several codon-based methods are available for detecting adaptive evolution in protein-coding sequences, but to date none specifically identify sites that are selected differentially in two populations, although such comparisons between populations have been historically useful in identifying the action of natural selection. We have developed two fixed effects maximum likelihood methods: one for identifying codon positions showing selection patterns that persist in a population and another for detecting whether selection is operating differentially on individual codons of a gene sampled from two different populations. Applying these methods to two HIV populations infecting genetically distinct human hosts, we have found that few of the positively selected amino acid sites persist in the population; the other changes are detected only at the tips of the phylogenetic tree and appear deleterious in the long term. Additionally, we have identified seven amino acid sites in protease and reverse transcriptase that are selected differentially in the two samples, demonstrating specific population-level adaptation of HIV to human populations.

Adaptation, Physiological↗

Neutralizing antibody responses drive the evolution of human immunodeficiency virus type 1 envelope during recent HIV infection.

HIV type 1 (HIV-1) can rapidly escape from neutralizing antibody responses. The genetic basis of this escape in vivo is poorly understood. We compared the pattern of evolution of the HIV-1 env gene between individuals with recent HIV infection whose virus exhibited either a low or a high rate of escape from neutralizing antibody responses. We demonstrate that the rate of viral escape at a phenotypic level is highly variable among individuals, and is strongly correlated with the rate of amino acid substitutions. We show that dramatic escape from neutralizing antibodies can occur in the relative absence of changes in glycosylation or insertions and deletions ("indels") in the envelope; conversely, changes in glycosylation and indels occur even in the absence of neutralizing antibody responses. Comparison of our data with the predictions of a mathematical model support a mechanism in which escape from neutralizing antibodies occurs via many amino acid substitutions, with low cross-neutralization between closely related viral strains. Our results suggest that autologous neutralizing antibody responses may play a pivotal role in the diversification of HIV-1 envelope during the early stages of infection.

Adult↗

Datamonkey: rapid detection of selective pressure on individual sites of codon alignments.

UNLABELLED: Datamonkey is a web interface to a suite of cutting edge maximum likelihood-based tools for identification of sites subject to positive or negative selection. The methods range from very fast data exploration to the some of the most complex models available in public domain software, and are implemented to run in parallel on a cluster of computers. AVAILABILITY: http://www.datamonkey.org. In the future, we plan to expand the collection of available analytic tools, and provide a package for installation on other systems.

Algorithms↗

Not so different after all: a comparison of methods for detecting amino acid sites under selection.

We consider three approaches for estimating the rates of nonsynonymous and synonymous changes at each site in a sequence alignment in order to identify sites under positive or negative selection: (1) a suite of fast likelihood-based "counting methods" that employ either a single most likely ancestral reconstruction, weighting across all possible ancestral reconstructions, or sampling from ancestral reconstructions; (2) a random effects likelihood (REL) approach, which models variation in nonsynonymous and synonymous rates across sites according to a predefined distribution, with the selection pressure at an individual site inferred using an empirical Bayes approach; and (3) a fixed effects likelihood (FEL) method that directly estimates nonsynonymous and synonymous substitution rates at each site. All three methods incorporate flexible models of nucleotide substitution bias and variation in both nonsynonymous and synonymous substitution rates across sites, facilitating the comparison between the methods. We demonstrate that the results obtained using these approaches show broad agreement in levels of Type I and Type II error and in estimates of substitution rates. Counting methods are well suited for large alignments, for which there is high power to detect positive and negative selection, but appear to underestimate the substitution rate. A REL approach, which is more computationally intensive than counting methods, has higher power than counting methods to detect selection in data sets of intermediate size but may suffer from higher rates of false positives for small data sets. A FEL approach appears to capture the pattern of rate variation better than counting methods or random effects models, does not suffer from as many false positives as random effects models for data sets comprising few sequences, and can be efficiently parallelized. Our results suggest that previously reported differences between results obtained by counting methods and random effects models arise due to a combination of the conservative nature of counting-based methods, the failure of current random effects models to allow for variation in synonymous substitution rates, and the naive application of random effects models to extremely sparse data sets. We demonstrate our methods on sequence data from the human immunodeficiency virus type 1 env and pol genes and simulated alignments.

Amino Acids↗

Phenotypic hypersusceptibility to multiple protease inhibitors and low replicative capacity in patients who are chronically infected with human immunodeficiency virus type 1.

Increased susceptibility to the protease inhibitors saquinavir and amprenavir has been observed in human immunodeficiency virus type 1 (HIV-1) with specific mutations in protease (V82T and N88S). Increased susceptibility to ritonavir has also been described in some viruses from antiretroviral agent-naive patients with primary HIV-1 infection in association with combinations of amino acid changes at polymorphic sites in the protease. Many of the viruses displaying increased susceptibility to protease inhibitors also had low replication capacity. In this retrospective study, we analyze the drug susceptibility phenotype and the replication capacity of virus isolates obtained at the peaks of viremia during five consecutive structured treatment interruptions in 12 chronically HIV-1-infected patients. Ten out of 12 patients had at least one sample with protease inhibitor hypersusceptibility (change </=0.4-fold) to one or more protease inhibitor. Hypersusceptibility to different protease inhibitors was observed at variable frequency, ranging from 38% to amprenavir to 11% to nelfinavir. Pairwise comparisons between susceptibilities for the protease inhibitors showed a consistent correlation among all pairs. There was also a significant relationship between susceptibility to protease inhibitors and replication capacity in all patients. Replication capacity remained stable over the course of repetitive cycles of structured treatment interruptions. We could find no association between in vitro replication capacity and in vivo plasma viral load doubling time and CD4(+) and CD8(+) T-cell counts at each treatment interruption. Several mutations were associated with hypersusceptibility to each protease inhibitor in a univariate analysis. This study extends the association between hypersusceptibility to protease inhibitors and low replication capacity to virus isolated from chronically infected patients and highlights the complexity of determining the genetic basis of this phenomenon. The potential clinical relevance of protease inhibitor hypersusceptibility and low replication capacity to virologic response to protease inhibitor-based therapies deserves to be investigated further.

Carbamates↗

Characterization of human immunodeficiency virus type 1 (HIV-1) envelope variation and neutralizing antibody responses during transmission of HIV-1 subtype B.

We analyzed neutralization sensitivity and genetic variation of transmitted subtype B human immunodeficiency virus type 1 (HIV-1) in eight recently infected men who have sex with men and the virus from the six subjects who infected them. In contrast to reports of heterosexual transmission of subtype C HIV-1, in which the transmitted virus appears to be more neutralization sensitive, we demonstrate that in our study population, relatively few phenotypic changes in neutralization sensitivity or genotypic changes in envelope occurred during transmission of subtype B HIV-1. We suggest that limited genetic variation within the infecting host reduces the likelihood of selective transmission of neutralization-sensitive HIV.

Acute Disease↗

HyPhy: hypothesis testing using phylogenies.

UNLABELLED: The HyPhypackage is designed to provide a flexible and unified platform for carrying out likelihood-based analyses on multiple alignments of molecular sequence data, with the emphasis on studies of rates and patterns of sequence evolution. AVAILABILITY: http://www.hyphy.org CONTACT: muse@stat.ncsu.edu SUPPLEMENTARY INFORMATION: HyPhydocumentation and tutorials are available at http://www.hyphy.org.

Algorithms↗

A genetic algorithm approach to detecting lineage-specific variation in selection pressure.

The ratio of nonsynonymous (dN) to synonymous (dS) substitution rates, omega, provides a measure of selection at the protein level. Models have been developed that allow omega to vary among lineages. However, these models require the lineages in which differential selection has acted to be specified a priori. We propose a genetic algorithm approach to assign lineages in a phylogeny to a fixed number of different classes of omega, thus allowing variable selection pressure without a priori specification of particular lineages. This approach can identify models with a better fit than a single-ratio model, and with fits that are better than (in an information theoretic sense) a fully local model, in which all lineages are assumed to evolve under different values of omega, but with far fewer parameters. By averaging over models which explain the data reasonably well, we can assess the robustness of our conclusions to uncertainty in model estimation. Our approach can also be used to compare results from models in which branch classes are specified a priori with a wide range of credible models. We illustrate our methods on primate lysozyme sequences and compare them with previous methods applied to the same data sets.

Algorithms↗

A simple hierarchical approach to modeling distributions of substitution rates.

Genetic sequence data typically exhibit variability in substitution rates across sites. In practice, there is often too little variation to fit a different rate for each site in the alignment, but the distribution of rates across sites may not be well modeled using simple parametric families. Mixtures of different distributions can capture more complex patterns of rate variation, but are often parameter-rich and difficult to fit. We present a simple hierarchical model in which a baseline rate distribution, such as a gamma distribution, is discretized into several categories, the quantiles of which are estimated using a discretized beta distribution. Although this approach involves adding only two extra parameters to a standard distribution, a wide range of rate distributions can be captured. Using simulated data, we demonstrate that a "beta-" model can reproduce the moments of the rate distribution more accurately than the distribution used to simulate the data, even when the baseline rate distribution is misspecified. Using hepatitis C virus and mammalian mitochondrial sequences, we show that a beta- model can fit as well or better than a model with multiple discrete rate categories, and compares favorably with a model which fits a separate rate category to each site. We also demonstrate this discretization scheme in the context of codon models specifically aimed at identifying individual sites undergoing adaptive or purifying evolution.

Animals↗

Genetic basis of hypersusceptibility to protease inhibitors and low replicative capacity of human immunodeficiency virus type 1 strains in primary infection.

The initial virus strains from as many as 12% of individuals with primary human immunodeficiency virus (HIV) infection have a 50% inhibitory concentration </=0.4-fold that of HIV type 1(NL4-3) (HIV-1(NL4-3)) to ritonavir (hypersusceptibility [HS]). There is also substantial variation in replicative capacity (RC) or an in vitro assay of the contributions of protease (PR) and reverse transcriptase to viral fitness. In chronically infected antiretrovirally treated patients, amprenavir HS has been associated with the mutation N88S in PR, but this mutation is not seen in untreated patients. In this study, virus strains from 182 cases of primary HIV infection were analyzed, and a highly significant association between HS and low RC (</=10% that of HIV-1(NL4-3)) was observed (P < 10(-6)). Multivariate analysis was used to determine the genotypic basis of ritonavir HS, analyzing all polymorphic amino acid sites and insertions from p7gag through PR. Decision tree models developed on the entire Gag-plus-PR data set and on PR alone gave overall correct classifications of 73 and 72%, respectively, on cross-validation. They were also able to predict low RC, with sensitivities of 69 and 62% and specificities of 84 and 70%, respectively. The analysis shows that ritonavir HS in untreated primary HIV infection is not associated with single mutations but with combinations of amino acids at polymorphic sites and that the same genotypes which confer HS to PR inhibitors confer low RC. This supports the view that variation in PR function is directly responsible for variation in fitness among strains in primary infection.

Decision Trees↗

Human immunodeficiency virus-1 RNA levels in cerebrospinal fluid exhibit a set point in clinically stable patients not receiving antiretroviral therapy.

To characterize, over time, cerebrospinal fluid (CSF) virus-load change in clinically stable patients, human immunodeficiency virus (HIV)-1 RNA levels were measured in serial CSF samples and in plasma samples obtained, during periods of 20 days-6 years, from 17 HIV-infected individuals not receiving antiretroviral treatment. Longitudinal trends were analyzed by linear regression and restricted maximum-likelihood techniques. CSF HIV-1 RNA levels varied within a restricted range (+/-0.5 log(10) copies/mL) around each subject-specific mean. Although 16 of the 17 subjects had CSF slopes not significantly different from zero, slopes that were more positive were associated with lower CD4 counts. In an individual patient, a CSF virus-load change >0.5 log(10) copies/mL may be clinically significant. Furthermore, our data suggest that, if the CSF virus load reflects the size of the reservoir of infected cells in the central nervous system, this reservoir may increase in those individuals with advanced immunosuppression but is stable, over several years, in patients without AIDS.

Adult↗

Transmission fitness of drug-resistant human immunodeficiency virus and the prevalence of resistance in the antiretroviral-treated population.

Although the prevalence of drug-resistant strains in primary human immunodeficiency virus (HIV) infection in North America has recently increased, their transmission fitness remains unknown. The present study estimated the frequency of transmission of drug-resistant HIV from patients receiving antiretroviral therapy using retrospective surveys of clinic data. It revealed that resistant virus was transmitted only approximately 20% as frequently as expected from these patients. Individuals with primary resistance may become a significant source of resistant strains.

Anti-HIV Agents↗

Amprenavir-resistant HIV-1 exhibits lopinavir cross-resistance and reduced replication capacity.

OBJECTIVES: To evaluate protease inhibitor (PI) cross-resistance and reductions in replication capacity conferred by amprenavir-selected mutations. METHODS: HIV-1IIIB variants derived from passage in increasing concentrations of amprenavir were studied, as well as 3'Gag/protease recombinants derived from them. These strains progressively accumulated mutations at codons 10, 46, 47, 50 and 84 in the protease as well as a p1/p6 cleavage site mutation at codon 449 in Gag. Their susceptibility (IC50) to various PI and their corresponding replication capacities were evaluated by a single-cycle growth assay and compared with measures using competitive cultures and p24 antigen production. RESULTS: Amprenavir susceptibility decreased with increasing numbers of protease mutations. Changes in lopinavir susceptibility paralleled changes in amprenavir susceptibility. Certain amprenavir-selected mutants conferred greater than 10-fold cross-resistance to lopinavir, including PrL10F/M46I/I50V-GagL449F (19-fold) and PrL10F/M46I/I47V/I50V-GagL449F (31-fold). Moreover, one isolate with only two mutations in the protease (L10F/84V) and GagL449F displayed a 7.7-fold increase in lopinavir IC50. Low-level cross-resistance to ritonavir and nelfinavir was also observed. The replication capacity of viruses containing either I84V or I50V was at least 90% lower than the reference virus in the single-cycle assay. The order of relative replication capacity was wild-type > L10F > L10F/I84V > L10F/M46I/I50V > L10F/M46I/I47V/I50V. CONCLUSION: These results indicate that until more comprehensive genotype-phenotype correlations between amprenavir and lopinavir susceptibility are established, phenotypic testing may be preferable to genotyping to detect cross-resistance, and should be considered when switching patients from a failing amprenavir-containing regimen. This study also provides data on the concordance of replication capacity measurements generated using rapid single-cycle growth and competition assays.

Carbamates↗

Selection of drug-resistant HIV-1 mutants in response to repeated structured treatment interruptions.

BACKGROUND: A new HIV-1 treatment strategy based on repeated structured treatment interruptions (STI) is currently being evaluated in clinical trials to determine whether immune cell-mediated control of viral replication can be stimulated by intermittent periods of viral replication. The potential for selection of drug-resistant quasi-species remains a major concern of such a treatment strategy. METHODS: Plasma and peripheral blood lymphocyte (PBL) samples from 12 patients who had three consecutive STIs were studied. Genotypic analysis was based on population and clonal sequencing. Drug susceptibility and their corresponding replication capacities were evaluated by a single-cycle growth assay. RESULTS: Consistent with a loss of phenotypic susceptibility to lamivudine, the M184V mutation was detected by genotypic analysis (direct and clonal sequencing) in plasma samples collected from two patients at the end of the second or third STI. Longitudinal analysis of patient samples revealed a step-wise increase in the M184V mutation in each patient virus population over successive STIs, despite the lower replicative capacity associated with this mutation in the absence of antiviral agents. CONCLUSION: Drug-resistant virus can rise to high frequencies in chronically HIV-1 infected individuals during consecutive STIs. Evolution of resistance is likely to be more important in patients with prior suboptimal therapies, particularly when few mutations are required for resistance. Maximum care should be taken in designing STI protocols that minimize development of drug-resistant mutations that may lead to treatment failure. Thus, drug-resistance testing may be useful before restarting treatment during STI studies.

Adult↗

Viral evolution during structured treatment interruptions in chronically human immunodeficiency virus-infected individuals.

We analyzed the evolution of the human immunodeficiency virus type 1 (HIV-1) env gene in 12 chronically infected individuals who underwent structured treatment interruptions (STIs). Analyses of length variation and of clonal sequences demonstrated highly unpredictable evolution, which may limit the strengthening of HIV-specific immune responses by STIs because of the variability in exposure to viral antigens.

Adult↗