Distinct patterns of natural selection in the reverse transcriptase gene of HIV-1 in the presence and absence of antiretroviral therapy

Evolutionary Analysis Provides Insight Into the Origin and Adaptation of HCV

Hepatitis C virus (HCV) belongs to the Hepacivirus genus and is genetically heterogeneous, with seven major genotypes further divided into several recognized subtypes. HCV origin was previously dated in a range between ∼200 and 1000 years ago. Hepaciviruses have been identified in several domestic and wild mammals, the largest viral diversity being observed in bats and rodents. The closest relatives of HCV were found in horses/donkeys (equine hepaciviruses, EHV). However, the origin of HCV as a human pathogen is still an unsolved puzzle. Using a selection-informed evolutionary model, we show that the common ancestor of extant HCV genotypes existed at least 3000 years ago (CI: 3192–5221 years ago), with the oldest genotypes being endemic to Asia. EHV originated around 1100 CE (CI: 291–1640 CE). These time estimates exclude that EHV transmission was mainly sustained by widespread veterinary practices and suggest that HCV originated from a single zoonotic event with subsequent diversification in human populations. We also describe a number of biologically important sites in the major HCV genotypes that have been positively selected and indicate that drug resistance-associated variants are significantly enriched at positively selected sites. HCV exploits several cell-surface molecules for cell entry, but only two of these (CD81 and OCLN) determine the species-specificity of infection. Herein evolutionary analyses do not support a long-standing association between primates and hepaciviruses, and signals of positive selection at CD81 were only observed in Chiroptera. No evidence of selection was detected for OCLN in any mammalian order. These results shed light on the origin of HCV and provide a catalog of candidate genetic modulators of HCV phenotypic diversity.

Genetic Diversity of HIV-1 Gene vif Among Treatment-Naive Brazilians

HIV-1 has the Vif protein, which binds to human antiviral proteins APOBEC3 to form complexes to be degraded by cellular proteolysis. To further explore HIV-1 diversity at the population level, we analyzed blood samples from 317 treatment-naive patients in Brazil. In this study, we explored the correlations of Vif polymorphisms with clinical parameters of the patients and found that mutation K22H is associated with low CD4⁺ cell counts and higher viral loads. Phylogenetic analysis of the vif gene indicated that subtype B was predominant in ∼77% (243/317) of the patients, followed by HIV-1 F ∼18% (56/317), and subtype C ∼4% (12/317); five samples were BF recombinants (∼1% of patients), and one was an AG recombinant. On the basis of the vif gene, we detected the presence of one AG and several previously unknown BF intersubtypes in this population. The global mean diversity, measured by pairwise distances, was 0.0931 ± 0.0006 among sequences of subtype B (n = 243), whereas the mean diversity of subtype C sequences (n = 12) was 0.0493 ± 0.001 and that of subtype F (n = 56) was 0.050 ± 0.001.

Increase in human immunodeficiency virus 1 diversity and detection of various subtypes and recombinants in north-eastern Brazil

PURPOSE

Diverse human immunodeficiency virus 1 (HIV-1) subtypes and circulating recombinant forms are found in Brazil. The majority of HIV-1 molecular epidemiological studies in Brazil have been conducted in the southern and south-eastern regions of the country, although several recent studies in the north-eastern region have addressed this issue. The objective of this study was to molecularly characterize HIV-1 circulating in Pernambuco, north-eastern Brazil.

METHODOLOGY

A total of 64 samples were collected from 2002 to 2003, and another 103 were collected from 2007 to 2009. The protease and partial reverse transcriptase regions of the HIV-1 polymerase-encoding (pol) gene were sequenced, and subtyping, recombination and phylogenetic analyses were performed.Results/Key findings. Subtype B (60.9 %) was found to be predominant, followed by HIV-1 F (31.4 %). Several BF recombinants (4.2 %), and BC and AG recombinants were also identified. The intra-subtype genetic diversity was estimated to be 0.065 (sd±0.004) for HIV-1 B and 0.055 (sd±0.004) for HIV-1 F, reflecting a greater accumulation of mutations in subtype B (P<0.01). More codons were found to be under positive selective pressure in samples collected from 2007 to 2009, from individuals with a T-cell count≥200 cells mm-3 and from women. Coalescence data indicated that the subtype F population has been continuously expanding.

CONCLUSIONS

HIV-1 shows high genetic diversity in the state of Pernambuco. Thus, additional molecular evaluations of circulating strains will provide a better understanding of the epidemic and may lead to more effective preventive strategies.

Selection pressure in CD8⁺ T-cell epitopes in the pol gene of HIV-1 infected individuals in Colombia. A bioinformatic approach

One of the main characteristics of the human immunodeficiency virus is its genetic variability and rapid adaptation to changing environmental conditions. This variability, resulting from the lack of proofreading activity of the viral reverse transcriptase, generates mutations that could be fixed either by random genetic drift or by positive selection. Among the forces driving positive selection are antiretroviral therapy and CD8+ T-cells, the most important immune mechanism involved in viral control. Here, we describe mutations induced by these selective forces acting on the pol gene of HIV in a group of infected individuals. We used Maximum Likelihood analyses of the ratio of non-synonymous to synonymous mutations per site (dN/dS) to study the extent of positive selection in the protease and the reverse transcriptase, using 614 viral sequences from Colombian patients. We also performed computational approaches, docking and algorithmic analyses, to assess whether the positively selected mutations affected binding to the HLA molecules. We found 19 positively-selected codons in drug resistance-associated sites and 22 located within CD8+ T-cell epitopes. A high percentage of mutations in these epitopes has not been previously reported. According to the docking analyses only one of those mutations affected HLA binding. However, algorithmic methods predicted a decrease in the affinity for the HLA molecule in seven mutated peptides. The bioinformatics strategies described here are useful to identify putative positively selected mutations associated with immune escape but should be complemented with an experimental approach to define the impact of these mutations on the functional profile of the CD8+ T-cells.

Relaxation of adaptive evolution during the HIV-1 infection owing to reduction of CD4+ T cell counts

BACKGROUND

The first stages of HIV-1 infection are essential to establish the diversity of virus population within host. It has been suggested that adaptation to host cells and antibody evasion are the leading forces driving HIV evolution at the initial stages of AIDS infection. In order to gain more insights on adaptive HIV-1 evolution, the genetic diversity was evaluated during the infection time in individuals contaminated by the same viral source in an epidemic cluster. Multiple sequences of V3 loop region of the HIV-1 were serially sampled from four individuals: comprising a single blood donor, two blood recipients, and another sexually infected by one of the blood recipients. The diversity of the viral population within each host was analyzed independently in distinct time points during HIV-1 infection.

RESULTS

Phylogenetic analysis identified multiple HIV-1 variants transmitted through blood transfusion but the establishing of new infections was initiated by a limited number of viruses. Positive selection (d(N)/d(S)>1) was detected in the viruses within each host in all time points. In the intra-host viruses of the blood donor and of one blood recipient, X4 variants appeared respectively in 1993 and 1989. In both patients X4 variants never reached high frequencies during infection time. The recipient, who X4 variants appeared, developed AIDS but kept narrow and constant immune response against HIV-1 during the infection time.

CONCLUSION

Slowing rates of adaptive evolution and increasing diversity in HIV-1 are consequences of the CD4+ T cells depletion. The dynamic of R5 to X4 shift is not associated with the initial amplitude of humoral immune response or intensity of positive selection.

Analysis of selection pressure and mutational pattern of HIV type 1 reverse transcriptase region among treated and nontreated patients

Variation of the HIV-1 subtype C reverse transcriptase region (RT) resulting in response to the selective pressures of drug therapy remains poorly characterized. Here, we compared the genetic variation resulting in the presence and absence of antiretroviral drug selective pressures on HIV-1 subtype C RT among nontreated and treated patients. The nucleotide variability, nonsynonymous and synonymous ratio, and the positively selected mutations were determined by comparing the RT sequences isolated at two time points among nontreated (baseline and follow-up) and treated patients (baseline and treatment failure). Compared to the nontreated patients, the intrapatient nucleotide variability, the number of nonsynonymous and synonymous substitutions was significantly higher among the treated patients. Among the mutations positively selected, the frequency of D121Y, I135R, and Q207E increased and the frequency of mutation S48T decreased significantly during treatment failure. Further studies are essential to discover the role of these mutations during treatment in HIV-1 subtype C.

Intrahost and interhost variability of the HIV type 1 nef gene in Brazilian children

Many aspects of HIV-1 pathogenesis are affected by Nef protein activity, and efforts have been made to study variation in the nef gene and how that variation relates to disease outcome. We studied the genetic diversity of the nef gene in distinct clones obtained from the same patient (intrahost) and in sequences obtained from different hosts (interhost). The set of sequences analyzed was obtained from HIV-1-infected Brazilian children and contained 112 clones from 25 children (intrahost samples), as well as 55 sequences from epidemiologically unlinked children (interhost samples). We found extensive site polymorphisms and amino acid length variations, mainly in the amino terminal region of the nef gene, between the myristoylation motif (MGxxxS) and the MHC-1 downregulation motif (Rxx). Analysis of the sequences deposited in the Los Alamos HIV sequences database ( www.hiv.lanl.gov ) indicated that the most frequent motif at the MHC-1 downregulation site in the subtype B strain is R(86%)A(64%)E(82%) (n = 1040) and R(78%)T(74%)E(56%) in the subtype C strain (n = 549). Conversely, the Brazilian subtype B isolates presented the motif R(81%)T(62%)E(67%) at this site (n = 64). A detailed analysis of selective pressures identified a concentration of codons under strong positive selection in the amino terminal region of the nef gene. We also determined that different sites are under positive selection in the subtype B and subtype C viruses. The amino acid composition in the MHC-1 downregulation motif of the nef gene in our sequences may indicate a distinct adaptive pattern of HIV-1 subtype B to the Brazilian host population.

High-resolution molecular epidemiology and evolutionary history of HIV-1 subtypes in Albania

Background

HIV-1 epidemic in Western Europe is largely due to subtype B. Little is known about the HIV-1 in Eastern Europe, but a few studies have shown that non-B subtypes are quite common. In Albania, where a recent study estimated a ten-fold increase of AIDS incidence during the last six years, subtype A and B account for 90% of the know infections.

Methodology/Principal Findings

We investigated the demographic history of HIV-1 subtype A and B in Albania by using a statistical framework based on coalescent theory and phylogeography. High-resolution phylogenetic and molecular clock analysis showed a limited introduction to the Balkan country of subtype A during the late 1980s followed by an epidemic outburst in the early 1990s. In contrast, subtype B was apparently introduced multiple times between the mid-1970s and mid-1980s. Both subtypes are growing exponentially, although the HIV-1A epidemic displays a faster growth rate, and a significantly higher basic reproductive number R₀. HIV-1A gene flow occurs primarily from the capital Tirane, in the center of the country, to the periphery, while HIV-1B flow is characterized by a balanced exchange between center and periphery. Finally, we calculated that the actual number of infections in Albania is at least two orders of magnitude higher than previously thought.

Conclusions/Significance

Our analysis demonstrates the power of recently developed computational tools to investigate molecular epidemiology of pathogens, and emphasize the complex factors involved in the establishment of HIV-1 epidemics. We suggest that a significant correlation exists between HIV-1 exponential spread and the socio-political changes occurred during the Balkan wars. The fast growth of a relatively new non-B epidemic in the Balkans may have significant consequences for the evolution of HIV-1 epidemiology in neighboring countries in Eastern and Western Europe.

Selective pressures of human immunodeficiency virus type 1 (HIV-1) during pediatric infection

Pediatric HIV-1 infection presents remarkable features that are distinct from those observed in adult infection. In vertically HIV-1-infected children, the viral load declines more slowly, and the cytotoxic T-lymphocyte response emerges late, only after the sixth month of life. This response generally tends to be narrow and less intense than that seen in adults. While the nuances of immune response at the cellular level during pediatric HIV-1 infection have been addressed, there is a lack of studies focusing on the consequences of this delayed and narrowed immune response at the population level. To better explore these features, we evaluated the selection regimen in gag, pol and env gene fragments of HIV-1 during pediatric infection. We estimated the number of nonsynonymous substitutions (d(N)) and synonymous substitutions (d(S)) codon-by-codon, using the maximum likelihood method and a modified counting method. Notably, both methods indicated a similar intensity of selection (measure by mean d(N)/d(S) ratio) between children and adults. Additionally, sites under positive selection were equally distributed along HIV genes and the location of these sites was analogous between children and adults. Therefore, the selective regimen in HIV during pediatric infection is equally broad and intense likewise the observed in adults. Unexpectedly, our phylogenetic-based analysis enabled us to identify two regions in the env gene of HIV with distinct adaptive functions. The first region, located in the vicinity of V3 loop, contains sites that might increase viral fitness within-host during antibody attack and virus-cell interaction. The second region, restricted to amino acids 334-368 of Gp160, contains sites that might increase viral fitness during interhost transmission at the population level.

Human immunodeficiency virus type 1 protease and reverse transcriptase mutation patterns among treatment-naive patients in different stages of infection in Rio de Janeiro, Brazil

The presence of genetic mutations in HIV-1-positive untreated individuals and its contribution to treatment failure, either in an individual or on a population basis, remains an important concern. The goal was to analyze and compare HIV-1 reverse transcriptase (RT) and protease (PR) genes of untreated individuals with chronic and recent infections. Fifty-one chronic infected individuals for whom initiation of antiretroviral treatment had been recommended and 20 individuals with recent documented HIV-1 seroconversion had their plasma viral RNA extracted and the PR and RT genes sequenced in order to determine subtype, presence of genetic polymorphisms and mutations associated with resistance to antiretroviral drugs. All 20 recent seroconvertors were infected with subtype B viruses. Of the 51 chronically infected patients, 40 (78.4%), 7 (13.7%), and 2 (3.9%) were infected with subtypes B, F, and C, respectively. Two (3.9%) hybrid forms were also observed in two individuals with chronic infection: D/B and D/F. Despite seroconversion stage, type and quantity of mutations were similar to both groups (P = 0.961). This group also presented the only (1.4%) drug-resistance mutation (M184V) among all samples investigated. In summary, the present study shows a high occurrence of equivalent polymorphisms unrelated to drug resistance in samples collected from untreated HIV-1- infected individuals in different seroconversion status, and suggests low primary resistance mutations. Results also indicate that non-B subtypes circulating in Rio de Janeiro have specific Brazilian non-synonymous mutations.

Core set approach to reduce uncertainty of gene trees

Background

A genealogy based on gene sequences within a species plays an essential role in the estimation of the character, structure, and evolutionary history of that species. Because intraspecific sequences are more closely related than interspecific ones, detailed information on the evolutionary process may be available by determining all the node sequences of trees and provide insight into functional constraints and adaptations. However, strong evolutionary correlations on a few lineages make this determination difficult as a whole, and the maximum parsimony (MP) method frequently allows a number of topologies with a same total branching length.

Results

Kitazoe et al. developed multidimensional vector-space representation of phylogeny. It converts additivity of evolutionary distances to orthogonality among the vectors expressing branches, and provides a unified index to measure deviations from the orthogoality. In this paper, this index is used to detect and exclude sequences with large deviations from orthogonality, and then selects a maximum subset ("core set") of sequences for which MP generates a single solution. Once the core set tree is formed whose all the node sequences are given, the excluded sequences are found to have basically two phylogenetic positions on this tree, respectively. Fortunately, since multiple substitutions are rare in intra-species sequences, the variance of nucleotide transitions is confined to a small range. By applying the core set approach to 38 partial env sequences of HIV-1 in a single patient and also 198 mitochondrial COI and COII DNA sequences of Anopheles dirus, we demonstrate how consistently this approach constructs the tree.

Conclusion

In the HIV dataset, we confirmed that the obtained core set tree is the unique maximum set for which MP proposes a single tree. In the mosquito data set, the fluctuation of nucleotide transitions caused by the sequences excluded from the core set was very small. We reproduced this core-set tree by simulation based on random process, and applied our approach to many sets of the obtained endpoint sequences. Consequently, the ninety percent of the endpoint sequences was identified as the core sets and the obtained node sequences were perfectly identical to the true ones.

Fold recognition of the human immunodeficiency virus type 1 V3 loop and flexibility of its crown structure during the course of adaptation to a host

The third hypervariable (V3) region of the HIV-1 gp120 protein is responsible for many aspects of viral infectivity. The tertiary structure of the V3 loop seems to influence the coreceptor usage of the virus, which is an important determinant of HIV pathogenesis. Hence, the information about preferred conformations of the V3-loop region and its flexibility could be a crucial tool for understanding the mechanisms of progression from an initial infection to AIDS. Taking into account the uncertainty of the loop structure, we predicted the structural flexibility, diversity, and sequence fitness to the V3-loop structure for each of the sequences serially sampled during an asymptomatic period. Structural diversity correlated with sequence diversity. The predicted crown structure usage implied that structural flexibility depended on the patient and that the antigenic character of the virus might be almost uniform in a patient whose immune system is strong. Furthermore, the predicted structural ensemble suggested that toward the end of the asymptomatic period there was a change in the V3-loop structure or in the environment surrounding the V3 loop, possibly because of its proximity to the gp120 core.

Fitness of drug-resistant malaria parasites

Drug-resistant mutant forms of an organism are likely to be less fit than their wild-type strains in the absence of selection. Experimental work on prokaryotic organisms suggests that this is the case, but that compensatory mutations may occur which restore the fitness of mutants to that of sensitive forms. Here, we review experimental and field studies on this subject in malaria. In the rodent model Plasmodium chabaudi, a pyrimethamine-resistant mutant has been found to grow more slowly in mice than its drug-sensitive progenitor; however, following passage in the absence of the drug it grew faster, suggesting the occurrence of compensatory mutations. Similar findings were made with a chloroquine-resistant mutant. Field studies on Plasmodium falciparum have provided circumstantial evidence of a loss of fitness of chloroquine-resistant mutants, which appear to become less frequent in the parasite population following withdrawal of the drug. However, the occurrence of frequent recombination in the life-cycle of this parasite means that in natural conditions, a gene conferring resistance, once it has arisen, can then spread into a diversity of genetically distinct backgrounds which will influence its fitness and capacity to survive in the parasite population.

Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups

We explored the epidemic history of HIV-1 subtype B in the United Kingdom by using statistical methods that infer the population history of pathogens from sampled gene sequence data. Phylogenetic analysis of HIV-1 pol gene sequences from Britain showed at least six large transmission chains, indicating a genetically variable, but epidemiologically homogeneous, epidemic among men having sex with men. Through coalescent-based analysis, we showed that these chains arose through separate introductions of subtype B strains into the United Kingdom in the early to mid-1980s. After an initial period of exponential growth, the rate of spread generally slowed in the early 1990s, which is more likely to correlate with behavior change than with reduced infectiousness resulting from highly active antiretroviral therapy. Our results provide insights into the complexity of HIV-1 epidemics that must be considered when developing HIV monitoring and prevention initiatives.