Clinical, virologic, and immunologic correlates of HIV-1 intraclade B dual infection among men who have sex with men.

OBJECTIVE: To investigate the susceptibilities to and consequences of HIV-1 dual infection. DESIGN: We compared clinical, virologic, and immunologic factors between participants who were dually infected with HIV-1 subtype B and monoinfected controls who were matched by ongoing HIV risk factor. METHODS: The viral load and CD4 progressions of dually and singly infected participant groups were compared with linear mixed-effects models, and individual dynamics before and after superinfection were assessed with a structural change test (Chow test). Recombination breakpoint analysis (GARD), HLA frequency analysis, and cytotoxic T-lymphocyte (CTL) epitope mapping were also performed (HIV LANL Database). RESULTS: The viral loads of dually infected participants increased more over 3 years of follow-up than the viral loads of monoinfected controls, whereas CD4 progressions of the two groups did not differ. Viral escape from CTL responses following superinfection was observed in two participants whose superinfecting strain completely replaced the initial strain. This pattern was not seen among participants whose superinfecting virus persisted in a recombinant form with the initial virus or was only detected transiently. Several HLA types were over-represented in dually infected participants as compared to monoinfected controls. CONCLUSIONS: These results identify potential factors for dual infection susceptibility and further define its clinical consequences.

Detection of minority resistance during early HIV-1 infection: natural variation and spurious detection rather than transmission and evolution of multiple viral variants.

Reports of a high frequency of the transmission of minority viral populations with drug-resistant mutations (DRM) are inconsistent with evidence that HIV-1 infections usually arise from mono- or oligoclonal transmission. We performed ultradeep sequencing (UDS) of partial HIV-1 gag, pol, and env genes from 32 recently infected individuals. We then evaluated overall and per-site diversity levels, selective pressure, sequence reproducibility, and presence of DRM and accessory mutations (AM). To differentiate biologically meaningful mutations from those caused by methodological errors, we obtained multinomial confidence intervals (CI) for the proportion of DRM at each site and fitted a binomial mixture model to determine background error rates for each sample. We then examined the association between detected minority DRM and the virologic failure of first-line antiretroviral therapy (ART). Similar to other studies, we observed increased detection of DRM at low frequencies (average, 0.56%; 95% CI, 0.43 to 0.69; expected UDS error, 0.21 ± 0.08% mutations/site). For 8 duplicate runs, there was variability in the proportions of minority DRM. There was no indication of increased diversity or selection at DRM sites compared to other sites and no association between minority DRM and AM. There was no correlation between detected minority DRM and clinical failure of first-line ART. It is unlikely that minority viral variants harboring DRM are transmitted and maintained in the recipient host. The majority of low-frequency DRM detected using UDS are likely errors inherent to UDS methodology or a consequence of error-prone HIV-1 replication.

A random effects branch-site model for detecting episodic diversifying selection.

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.

Using phylogeography to characterize the origins of the HIV-1 subtype F epidemic in Romania.

BACKGROUND: During the late 1980s and early 1990s, an estimated 10,000 Romanian children were infected with HIV-1 subtype F nosocomially through contaminated needles and blood transfusions. However, the geographic source and origins of this epidemic remain unclear. METHODS: Here we used phylogenetic inference and "relaxed" molecular clock dating analysis to further characterize the Romanian HIV-1 subtype F epidemic. RESULTS: These analyses revealed a major lineage of Romanian HIV sequences consisting nearly entirely of virus sampled from adolescents and children and a distinct cluster that included a much higher ratio of adult sequences. Divergence time estimates inferred the time of most recent common ancestor of subtype F1 sequences to be 1973 (1966-1980) and for all Angolan sequences to 1975 (1968-1980). The most common ancestor of the Romanian sequences was dated to 1978 (1972-1983) with pediatric and adolescent sequences interspersed throughout the lineage. The phylogenetic structure of the entire subtype F epidemic suggests that multiple introductions of subtype F into Romania occurred either from the Angolan epidemic or from more distant ancestors. Since the historical records note that the Romanian pediatric epidemic did not begin until the late 1980s, the inferred time of most recent common ancestor of the Romanian lineage of 1978 suggests that there were multiple introductions of subtype F occurred into the pediatric population from HIV already circulating in Romania. CONCLUSIONS: Analysis of the subtype F HIV-1 epidemic in an historical context allows for a deeper appreciation of how the HIV pandemic has been influenced by socio-political events.

The relatedness of HIV epidemics in the United States-Mexico border region.

Phylogeography can improve the understanding of local and worldwide HIV epidemics, including the migration of subepidemics across national borders. We analyzed HIV-1 sequences sampled from Mexico and San Diego, California to determine the relatedness of these epidemics. We sampled the HIV epidemics in (1) Mexico by downloading all publicly available HIV-1 pol sequences from antiretroviral-naive individuals in GenBank (n = 100) and generating similar sequences from cohorts of injection drug users and female sex workers in Tijuana, Mexico (n = 27) and (2) in San Diego, California by pol sequencing well-characterized primary (n = 395) and chronic (n = 267) HIV infection cohorts. Estimates of population structure (F(ST)), genetic distance cluster analysis, and a cladistic measure of migration events (Slatkin-Maddison test) were used to assess the relatedness of the epidemics. Both a test of population differentiation (F(ST) = 0.06; p < 0.01) and a cladistic estimate of migration events (84 migrations, p < 0.01) indicated that the Tijuana and San Diego epidemics were not freely mixing. A conservative cluster analysis identified 72 clusters (two or more sequences), with two clusters containing both Mexican and San Diego sequences (permutation p < 0.01). Analysis of this very large dataset of HIV-1 sequences suggested that the HIV-1 epidemics in San Diego, California and Tijuana, Mexico are distinct. Larger epidemiological studies are needed to quantify the magnitude and associations of cross-border mixing.

CodonTest: modeling amino acid substitution preferences in coding sequences.

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.

Correcting the bias of empirical frequency parameter estimators in codon models.

Markov models of codon substitution are powerful inferential tools for studying biological processes such as natural selection and preferences in amino acid substitution. The equilibrium character distributions of these models are almost always estimated using nucleotide frequencies observed in a sequence alignment, primarily as a matter of historical convention. In this note, we demonstrate that a popular class of such estimators are biased, and that this bias has an adverse effect on goodness of fit and estimates of substitution rates. We propose a "corrected" empirical estimator that begins with observed nucleotide counts, but accounts for the nucleotide composition of stop codons. We show via simulation that the corrected estimates outperform the de facto standard estimates not just by providing better estimates of the frequencies themselves, but also by leading to improved estimation of other parameters in the evolutionary models. On a curated collection of sequence alignments, our estimators show a significant improvement in goodness of fit compared to the approach. Maximum likelihood estimation of the frequency parameters appears to be warranted in many cases, albeit at a greater computational cost. Our results demonstrate that there is little justification, either statistical or computational, for continued use of the -style estimators.

Benchmarking multi-rate codon models.

The single rate codon model of non-synonymous substitution is ubiquitous in phylogenetic modeling. Indeed, the use of a non-synonymous to synonymous substitution rate ratio parameter has facilitated the interpretation of selection pressure on genomes. Although the single rate model has achieved wide acceptance, we argue that the assumption of a single rate of non-synonymous substitution is biologically unreasonable, given observed differences in substitution rates evident from empirical amino acid models. Some have attempted to incorporate amino acid substitution biases into models of codon evolution and have shown improved model performance versus the single rate model. Here, we show that the single rate model of non-synonymous substitution is easily outperformed by a model with multiple non-synonymous rate classes, yet in which amino acid substitution pairs are assigned randomly to these classes. We argue that, since the single rate model is so easy to improve upon, new codon models should not be validated entirely on the basis of improved model fit over this model. Rather, we should strive to both improve on the single rate model and to approximate the general time-reversible model of codon substitution, with as few parameters as possible, so as to reduce model over-fitting. We hint at how this can be achieved with a Genetic Algorithm approach in which rate classes are assigned on the basis of sequence information content.

Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology.

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.

Genome-wide analysis of the structure of the South African Coloured Population in the Western Cape.

Admixed populations present unique opportunities to discover the genetic factors underlying many multifactorial diseases. The geographical position and complex history of South Africa has led to the establishment of the unique admixed population known as the South African Coloured. Not much is known about the genetic make-up of this population, and the historical record is patchy. We genotyped 959 individuals from the Western Cape area, self-identified as belonging to this population, using the Affymetrix 500k genotyping platform. This resulted in nearly 75,000 autosomal SNPs that could be compared with populations represented in the International HapMap Project and the Human Genome Diversity Project. Analysis by means of both the admixture and linkage models in STRUCTURE revealed that the major ancestral components of this population are predominantly Khoesan (32-43%), Bantu-speaking Africans (20-36%), European (21-28%) and a smaller Asian contribution (9-11%), depending on the model used. This is consistent with historical data. While of great historical and genealogical interest, this information is also essential for future admixture mapping of disease genes in this population.

Pathogenic peptide deviations support a model of adaptive evolution of chordate cardiac performance by troponin mutations.

In cardiac muscle, the troponin (cTn) complex is a key regulator of myofilament calcium sensitivity because it serves as a molecular switch required for translating myocyte calcium fluxes into sarcomeric contraction and relaxation. Studies of several species suggest that ectotherm chordates have myofilaments with heightened calcium responsiveness. However, genetic polymorphisms in cTn that cause increased myofilament sensitivity to activating calcium in mammals result in cardiac disease including arrhythmias, diastolic dysfunction, and increased susceptibility to sudden cardiac death. We hypothesized that specific residue modifications in the regulatory arm of troponin I (TnI) were critical in mediating the observed decrease in myofilament calcium sensitivity within the mammalian taxa. We performed large-scale phylogenetic analysis, atomic resolution molecular dynamics simulations and modeling, and computational alanine scanning. This study provides evidence that a His to Ala substitution within mammalian cardiac TnI (cTnI) reduced the thermodynamic potential at the interface between cTnI and cardiac TnC (cTnC) in the calcium-saturated state by disrupting a strong intermolecular electrostatic interaction. This key residue modification reduced myofilament calcium sensitivity by making cTnI molecularly untethered from cTnC. To meet the requirements for refined mammalian adult cardiac performance, we propose that compensatory evolutionary pressures favored mutations that enhanced the relaxation properties of cTn by decreasing its sensitivity to activating calcium.

The origins of sexually transmitted HIV among men who have sex with men.

Although it is known that most HIV-1 infections worldwide result from exposure to virus in semen, it has not yet been established whether transmitted strains originate as RNA virions in seminal plasma or as integrated proviral DNA in infected seminal leukocytes. We present phylogenetic evidence that among six transmitting pairs of men who have sex with men, blood plasma virus in the recipient is consistently more closely related to the seminal plasma virus in the source. All sequences were subtype B, and the env C2V3 of transmitted variants tended to have higher mean isoelectric points, contain potential N-linked glycosylation sites, and favor CCR5 co-receptor usage. A statistically robust phylogenetically corrected analysis did not detect genetic signatures reliably associated with transmission, but further investigation of larger samples of transmitting pairs holds promise for determining which structural and genetic features of viral genomes are associated with transmission.

Response to Comment on "The Origins of Sexually Transmitted HIV Among Men Who Have Sex with Men"

In this investigation, we evaluated explanations for the unusual degree of genetic diversity in HIV populations within the sampled seminal cell and plasma compartments observed in our previous study. These analyses included clonal sequencing of HIV DNA in peripheral blood mononuclear cells, ultradeep sequencing of HIV RNA in blood plasma, human leukocyte antigen (HLA) haplotyping of previously used samples, and a BLAST screen against both a local and public repository of HIV-1 sequences, and the investigations to determine whether these observations were secondary to contamination or artifact were unsuccessful. As there are very few HIV sequences from seminal cell tissues and transmission pairs described in the literature, future studies that evaluate more transmission pairs with sampling from multiple anatomic compartments and at multiple time points will most likely be required to resolve this controversy.

Experimental evidence indicating that mastreviruses probably did not co-diverge with their hosts.

BACKGROUND: Despite the demonstration that geminiviruses, like many other single stranded DNA viruses, are evolving at rates similar to those of RNA viruses, a recent study has suggested that grass-infecting species in the genus Mastrevirus may have co-diverged with their hosts over millions of years. This "co-divergence hypothesis" requires that long-term mastrevirus substitution rates be at least 100,000-fold lower than their basal mutation rates and 10,000-fold lower than their observable short-term substitution rates. The credibility of this hypothesis, therefore, hinges on the testable claim that negative selection during mastrevirus evolution is so potent that it effectively purges 99.999% of all mutations that occur. RESULTS: We have conducted long-term evolution experiments lasting between 6 and 32 years, where we have determined substitution rates of between 2 and 3 x 10(-4) substitutions/site/year for the mastreviruses Maize streak virus (MSV) and Sugarcane streak Réunion virus (SSRV). We further show that mutation biases are similar for different geminivirus genera, suggesting that mutational processes that drive high basal mutation rates are conserved across the family. Rather than displaying signs of extremely severe negative selection as implied by the co-divergence hypothesis, our evolution experiments indicate that MSV and SSRV are predominantly evolving under neutral genetic drift. CONCLUSION: The absence of strong negative selection signals within our evolution experiments and the uniformly high geminivirus substitution rates that we and others have reported suggest that mastreviruses cannot have co-diverged with their hosts.

Models of coding sequence evolution.

Probabilistic models of sequence evolution are in widespread use in phylogenetics and molecular sequence evolution. These models have become increasingly sophisticated and combined with statistical model comparison techniques have helped to shed light on how genes and proteins evolve. Models of codon evolution have been particularly useful, because, in addition to providing a significant improvement in model realism for protein-coding sequences, codon models can also be designed to test hypotheses about the selective pressures that shape the evolution of the sequences. Such models typically assume a phylogeny and can be used to identify sites or lineages that have evolved adaptively. Recently some of the key assumptions that underlie phylogenetic tests of selection have been questioned, such as the assumption that the rate of synonymous changes is constant across sites or that a single phylogenetic tree can be assumed at all sites for recombining sequences. While some of these issues have been addressed through the development of novel methods, others remain as caveats that need to be considered on a case-by-case basis. Here, we outline the theory of codon models and their application to the detection of positive selection. We review some of the more recent developments that have improved their power and utility, laying a foundation for further advances in the modeling of coding sequence evolution.

Frequent toggling between alternative amino acids is driven by selection in HIV-1.

Host immune responses against infectious pathogens exert strong selective pressures favouring the emergence of escape mutations that prevent immune recognition. Escape mutations within or flanking functionally conserved epitopes can occur at a significant cost to the pathogen in terms of its ability to replicate effectively. Such mutations come under selective pressure to revert to the wild type in hosts that do not mount an immune response against the epitope. Amino acid positions exhibiting this pattern of escape and reversion are of interest because they tend to coincide with immune responses that control pathogen replication effectively. We have used a probabilistic model of protein coding sequence evolution to detect sites in HIV-1 exhibiting a pattern of rapid escape and reversion. Our model is designed to detect sites that toggle between a wild type amino acid, which is susceptible to a specific immune response, and amino acids with lower replicative fitness that evade immune recognition. Through simulation, we show that this model has significantly greater power to detect selection involving immune escape and reversion than standard models of diversifying selection, which are sensitive to an overall increased rate of non-synonymous substitution. Applied to alignments of HIV-1 protein coding sequences, the model of immune escape and reversion detects a significantly greater number of adaptively evolving sites in env and nef. In all genes tested, the model provides a significantly better description of adaptively evolving sites than standard models of diversifying selection. Several of the sites detected are corroborated by association between Human Leukocyte Antigen (HLA) and viral sequence polymorphisms. Overall, there is evidence for a large number of sites in HIV-1 evolving under strong selective pressure, but exhibiting low sequence diversity. A phylogenetic model designed to detect rapid toggling between wild type and escape amino acids identifies a larger number of adaptively evolving sites in HIV-1, and can in some cases correctly identify the amino acid that is susceptible to the immune response.

Ecological speciation in South Atlantic island finches.

Examples of sympatric speciation in nature are rare and hotly debated. We describe the parallel speciation of finches on two small islands in the Tristan da Cunha archipelago in the South Atlantic Ocean. Nesospiza buntings are a classic example of a simple adaptive radiation, with two species on each island: an abundant small-billed dietary generalist and a scarce large-billed specialist. Their morphological diversity closely matches the available spectrum of seed sizes, and genetic evidence suggests that they evolved independently on each island. Speciation is complete on the smaller island, where there is a single habitat with strongly bimodal seed size abundance, but is incomplete on the larger island, where a greater diversity of habitats has resulted in three lineages. Our study suggests that the buntings have undergone parallel ecological speciation.

The transmission of Helicobacter pylori: the effects of analysis method and study population on inference.

Although much is known about the virulence of Helicobacter pylori, the transmission pathways for this bacterium are still unresolved. Transmission has been addressed through: (1) prevalence within families; (2) detection in fecal/oral environments; (3) detection in the abiotic/biotic environment; and (4) direct inference from strain similarity. Here, we review the molecular and biochemical methods used and discuss the relative merits of each. Furthermore, as there are differences between developing and developed nations, we discuss the results obtained from transmission studies in light of the study population. We conclude that H. pylori is probably transmitted person-to-person, facilitated by fecal-oral transmission during episodes of diarrhea or gastro-oral contact during periods of vomiting. The persistence of H. pylori in abiotic and biotic environments remains unproven but possible reactivation from viable, non-culturable coccoid forms should be further investigated. Finally, we speculate on the effect of host-pathogen interactions in confounding the inference of transmission.

Population growth confounds phylogeographic inference in Namaqua sandgrouse.

The Namaqua sandgrouse, Pterocles namaqua, is a highly nomadic granivore of semiarid to arid habitats. As a result of nomadic movements in response to rainfall, the size of the breeding population in any one area fluctuates dramatically between breeding seasons. This high mobility in response to spatial and temporal abundance of food resources is expected to result in little population genetic structuring. Namaqua sandgrouse also shows a seasonally predictable partial migration between the southeast and northwest regions of South Africa, and a further possible north-south migration between southwestern South Africa and central Namibia. It is unclear whether birds migrating between these regions breed in only one or both regions. If populations breed in only one region of their migratory range, then population genetic structuring is predicted to occur. This study addresses Namaqua sandgrouse movements with the analysis of mitochondrial DNA control region sequences. In general, little population genetic structure was evident, yet strong signals of population growth were detected. Several populations have private alleles, which is in direct contradiction to the spatial genetic pattern expected under high levels of gene flow. We suggest that the inference of high levels of female gene flow could be an artifact of population growth and that additional loci will allow a greater understanding of Namaqua sandgrouse movements.

A population genetics pedigree perspective on the transmission of Helicobacter pylori.

The inference of transmission pathways for medicinally important bacteria is important to our understanding of pathogens. Here we report analyses of transmission in Helicobacter pylori, a major carcinogen. Our study is novel in that the focal community comprises detailed family pedigrees and has a high prevalence of H. pylori. To infer transmission, we performed high-resolution analyses of nucleotide sequences for three genes and accounted for the occurrence of mutation and recombination through the use of simulation modeling. Our results demonstrate that transmission has a strong nonfamilial component potentially the result of a large proportion of infections derived from the community. These results are interesting from both a medical and an evolutionary standpoint. First, efficient control measures and beliefs about the sources of H. pylori infection should be reevaluated. Evolutionarily, our results contradict the hypothesis of strict vertical transmission, presented as an explanation for the strong correlation between human population history and H. pylori diversity. Thus the paradox of persistent phylogenetic structure, despite a permissive mode of transmission and high recombination rates, must be solved elsewhere. Here we consider the potential for recombination events to maintain genetic structure in light of horizontal transmission.