Simon D. W. Frost

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.

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.

Datamonkey 2010

Datamonkey is a popular web-based suite of phylogenetic analysis tools for use in evolutionary biology. Since the original release in 2005, we have expanded the analysis options to include recently developed algorithmic methods for recombination detection, evolutionary fingerprinting of genes, codon model selection, co-evolution between sites, identification of sites, which rapidly escape host-immune pressure and HIV-1 subtype assignment. The traditional selection tools have also been augmented to include recent developments in the field. Here, we summarize the analyses options currently available on Datamonkey, and provide guidelines for their use in evolutionary biology. Availability and documentation: http://www.datamonkey.org.

GARD: a genetic algorithm for recombination detection

MOTIVATION Phylogenetic and evolutionary inference can be severely misled if recombination is not accounted for, hence screening for it should be an essential component of nearly every comparative study. The evolution of recombinant sequences can not be properly explained by a single phylogenetic tree, but several phylogenies may be used to correctly model the evolution of non-recombinant fragments. RESULTS We developed a likelihood-based model selection procedure that uses a genetic algorithm to search multiple sequence alignments for evidence of recombination breakpoints and identify putative recombinant sequences. GARD is an extensible and intuitive method that can be run efficiently in parallel. Extensive simulation studies show that the method nearly always outperforms other available tools, both in terms of power and accuracy and that the use of GARD to screen sequences for recombination ensures good statistical properties for methods aimed at detecting positive selection. AVAILABILITY Freely available http://www.datamonkey.org/GARD/

An Evolutionary Model-Based Algorithm for Accurate Phylogenetic Breakpoint Mapping and Subtype Prediction in HIV-1

Genetically diverse pathogens (such as Human Immunodeficiency virus type 1, HIV-1) are frequently stratified into phylogenetically or immunologically defined subtypes for classification purposes. Computational identification of such subtypes is helpful in surveillance, epidemiological analysis and detection of novel variants, e.g., circulating recombinant forms in HIV-1. A number of conceptually and technically different techniques have been proposed for determining the subtype of a query sequence, but there is not a universally optimal approach. We present a model-based phylogenetic method for automatically subtyping an HIV-1 (or other viral or bacterial) sequence, mapping the location of breakpoints and assigning parental sequences in recombinant strains as well as computing confidence levels for the inferred quantities. Our Subtype Classification Using Evolutionary ALgorithms (SCUEAL) procedure is shown to perform very well in a variety of simulation scenarios, runs in parallel when multiple sequences are being screened, and matches or exceeds the performance of existing approaches on typical empirical cases. We applied SCUEAL to all available polymerase (pol) sequences from two large databases, the Stanford Drug Resistance database and the UK HIV Drug Resistance Database. Comparing with subtypes which had previously been assigned revealed that a minor but substantial (≈5%) fraction of pure subtype sequences may in fact be within- or inter-subtype recombinants. A free implementation of SCUEAL is provided as a module for the HyPhy package and the Datamonkey web server. Our method is especially useful when an accurate automatic classification of an unknown strain is desired, and is positioned to complement and extend faster but less accurate methods. Given the increasingly frequent use of HIV subtype information in studies focusing on the effect of subtype on treatment, clinical outcome, pathogenicity and vaccine design, the importance of accurate, robust and extensible subtyping procedures is clear.

Persistence of Transmitted Drug Resistance among Subjects with Primary Human Immunodeficiency Virus Infection

ABSTRACT Following interruption of antiretroviral therapy among individuals with acquired drug resistance, preexisting drug-sensitive virus emerges relatively rapidly. In contrast, wild-type virus is not archived in individuals infected with drug-resistant human immunodeficiency virus (HIV) and thus cannot emerge rapidly in the absence of selective drug pressure. Fourteen recently HIV-infected patients with transmitted drug-resistant virus were followed for a median of 2.1 years after the estimated date of infection (EDI) without receiving antiretroviral therapy. HIV drug resistance and pol replication capacity (RC) in longitudinal plasma samples were assayed. Resistance mutations were characterized as pure populations or mixtures. The mean time to first detection of a mixture of wild-type and drug-resistant viruses was 96 weeks (1.8 years) (95% confidence interval, 48 to 192 weeks) after the EDI. The median time to loss of detectable drug resistance using population-based assays ranged from 4.1 years (conservative estimate) to longer than the lifetime of the individual (less conservative estimate). The transmission of drug-resistant virus was not associated with virus with reduced RC. Sexual transmission of HIV selects for highly fit drug-resistant variants that persist for years. The prolonged persistence of transmitted drug resistance strongly supports the routine use of HIV resistance genotyping for all newly diagnosed individuals.

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.

Evaluation of Respondent-driven Sampling

Background: Respondent-driven sampling is a novel variant of link-tracing sampling for estimating the characteristics of hard-to-reach groups, such as HIV prevalence in sex workers. Despite its use by leading health organizations, the performance of this method in realistic situations is still largely unknown. We evaluated respondent-driven sampling by comparing estimates from a respondent-driven sampling survey with total population data. Methods: Total population data on age, tribe, religion, socioeconomic status, sexual activity, and HIV status were available on a population of 2402 male household heads from an open cohort in rural Uganda. A respondent-driven sampling (RDS) survey was carried out in this population, using current methods of sampling (RDS sample) and statistical inference (RDS estimates). Analyses were carried out for the full RDS sample and then repeated for the first 250 recruits (small sample). Results: We recruited 927 household heads. Full and small RDS samples were largely representative of the total population, but both samples underrepresented men who were younger, of higher socioeconomic status, and with unknown sexual activity and HIV status. Respondent-driven sampling statistical inference methods failed to reduce these biases. Only 31%–37% (depending on method and sample size) of RDS estimates were closer to the true population proportions than the RDS sample proportions. Only 50%–74% of respondent-driven sampling bootstrap 95% confidence intervals included the population proportion. Conclusions: Respondent-driven sampling produced a generally representative sample of this well-connected nonhidden population. However, current respondent-driven sampling inference methods failed to reduce bias when it occurred. Whether the data required to remove bias and measure precision can be collected in a respondent-driven sampling survey is unresolved. Respondent-driven sampling should be regarded as a (potentially superior) form of convenience sampling method, and caution is required when interpreting findings based on the sampling method.

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.

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

ABSTRACT 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.