Phylogenetic reconstruction of intrapatient evolution of human immunodeficiency virus type 1: predominance of drift and purifying selection

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts

The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5′ untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5′ UTR. As expected, mutations in double-stranded (DS) regions of the 5′ UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5′ cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.

IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Rapid evolution of the env gene leader sequence in cats naturally infected with feline immunodeficiency virus

Analysing the evolution of feline immunodeficiency virus (FIV) at the intra-host level is important in order to address whether the diversity and composition of viral quasispecies affect disease progression. We examined the intra-host diversity and the evolutionary rates of the entire env and structural fragments of the env sequences obtained from sequential blood samples in 43 naturally infected domestic cats that displayed different clinical outcomes. We observed in the majority of cats that FIV env showed very low levels of intra-host diversity. We estimated that env evolved at a rate of 1.16×10(-3) substitutions per site per year and demonstrated that recombinant sequences evolved faster than non-recombinant sequences. It was evident that the V3-V5 fragment of FIV env displayed higher evolutionary rates in healthy cats than in those with terminal illness. Our study provided the first evidence that the leader sequence of env, rather than the V3-V5 sequence, had the highest intra-host diversity and the highest evolutionary rate of all env fragments, consistent with this region being under a strong selective pressure for genetic variation. Overall, FIV env displayed relatively low intra-host diversity and evolved slowly in naturally infected cats. The maximum evolutionary rate was observed in the leader sequence of env. Although genetic stability is not necessarily a prerequisite for clinical stability, the higher genetic stability of FIV compared with human immunodeficiency virus might explain why many naturally infected cats do not progress rapidly to AIDS.

Quantification of unintegrated HIV-1 DNA at the single cell level in vivo

In the nucleus of HIV-1 infected cells, unintegrated HIV-1 DNA molecules exist in the form of one and two LTR circles and linear molecules with degraded extremities. In tissue culture they are invariably more numerous than the provirus, the relative proportion of integrated to unintegrated forms varies widely from ∼1∶1 to 1∶10 and even over 1∶100. In vivo, this ratio is unknown. To determine it, single nuclei from two infected patients with a known provirus copy number were microdissected, HIV DNA was amplified by nested PCR, cloned and individual clones sequenced. Given the extraordinary sequence complexity, we made the assumption that the total number of distinct sequences approximated to real number of amplifiable HIV-1 DNA templates in the nucleus. We found that the number of unintegrated DNA molecules increased linearly with the proviral copy number there being on average 86 unintegrated molecules per provirus.

Stable multi-infection of splenocytes during SIV infection--the basis for continuous recombination

Background

Recombination is an important mechanism in the generation of genetic diversity of the human (HIV) and simian (SIV) immunodeficiency viruses. It requires the co-packaging of divergent RNA genomes into the same retroviral capsid and subsequent template switching during the reverse transcription reaction. By HIV-specific fluorescence in situ hybridization (FISH), we have previously shown that the splenocytes from 2 chronically infected patients with Castelman's disease were multi-infected and thus fulfill the in vivo requirements to generate genetic diversity by recombination. In order to analyze when multi-infection first occurs during a lentivirus infection and how the distribution of multi-infection evolves during the disease course, we now determined the SIV copy numbers from splenocytes of 11 SIVmac251-infected rhesus macaques cross-sectionally covering the time span of primary infection throughout to end-stage immunodeficiency.

Results

SIV multi-infection of single splenocytes was readily detected in all monkeys and all stages of the infection. Single-infected cells were more frequent than double- or triple- infected cells. There was no strong trend linking the copy number distribution to plasma viral load, disease stage, or CD4 cell counts.

Conclusions

SIV multi-infection of single cells is already established during the primary infection phase thus enabling recombination to affect viral evolution in vivo throughout the disease course.

Characterization of human endogenous retroviral elements in the blood of HIV-1-infected individuals

We previously reported finding the RNA of a type K human endogenous retrovirus, HERV-K (HML-2), at high titers in the plasma of HIV-1-infected and cancer patients (R. Contreras-Galindo et al., J. Virol. 82:9329-9236, 2008.). The extent to which the HERV-K (HML-2) proviruses become activated and the nature of their activated viral RNAs remain important questions. Therefore, we amplified and sequenced the full-length RNA of the env gene of the type 1 and 2 HERV-K (HML-2) viruses collected from the plasma of seven HIV-1-infected patients over a period of 1 to 3 years and from five breast cancer patients in order to reconstruct the genetic evolution of these viruses. HERV-K (HML-2) RNA was found in plasma fractions of HIV-1 patients at a density of ∼1.16 g/ml that contained both immature and correctly processed HERV-K (HML-2) proteins and virus-like particles that were recognized by anti-HERV-K (HML-2) antibodies. RNA sequences from novel HERV-K (HML-2) proviruses were discovered, including K111, which is specifically active during HIV-1 infection. Viral RNA arose from complete proviruses and proviruses devoid of a 5' long terminal repeat, suggesting that the expression of HERV-K (HML-2) RNA in these patients may involve sense and antisense transcription. In HIV-1-infected individuals, the HERV-K (HML-2) viral RNA showed evidence of frequent recombination, accumulation of synonymous rather than nonsynonymous mutations, and conserved N-glycosylation sites, suggesting that some of the HERV-K (HML-2) viral RNAs have undergone reverse transcription and are under purifying selection. In contrast, HERV-K (HML-2) RNA sequences found in the blood of breast cancer patients showed no evidence of recombination and exhibited only sporadic viral mutations. This study suggests that HERV-K (HML-2) is active in HIV-1-infected patients, and the resulting RNA message reveals previously undiscovered HERV-K (HML-2) genomic sequences.

Monitoring of drug resistance in therapy-naïve HIV infected patients and detection of African HIV subtypes in Hungary

Mutations in the HIV-1 pol gene associated with resistance to antiretroviral drugs in therapy-naïve Hungarian individuals transmitted as primary infection by their foreign sexual partners originated from African, Asian and other European countries had been analyzed. Drug resistance genotyping of HIV RT and PR genes were performed where mutations of 72 codons - among them 64 specific resistance codons representing 6 nucleoside reverse transcriptase inhibitor (NRTIs), 2 non-nucleoside reverse transcriptase inhibitor (NNRTIs) and 6 proteinase inhibitor (PRIs) drugs - had been analyzed by Truegene HIV-1 Genotyping kit and OpenGene Sequencing System. Viral variants harboring resistance mutations in the po l gene were detected in 14% of the subjects. The highest rate of resistance to a single class of inhibitors was detected towards PR inhibitors (12%), followed by NRTI (8%) and NNRTI (5%). On the contrary, 25% of viruses transmitted by homosexual activity contained mutations led to resistance to NNRT. Viruses from 11 percent of cases were resistant to 2 classes of inhibitors, and 7 percent to three classes of inhibitors. Based upon sequence data non-B subtypes and CRFs were detected in more than 71% of cases. HIV-1 C (10.7%), HIV-F1 (7.2%) and HIV-1 G (3.6%) were detected as the more frequent subtypes. Among the HIV-1 recombinant viruses CRF02_AG variants were found more frequently (28.5%) followed by CRF06_cpx (17.8%) indicating penetration of non-B subtypes and recombinant African variants into Hungary, which raises serious clinical and public health consequences.

Emergence of antiretroviral drug resistance in therapy-naive HIV infected patients in Hungary

Mutations in the HIV-1 genes associated with resistance to antiretroviral drugs were detected also in primary HIV infected individuals who did not receive antiretroviral treatment. Drug resistance genotyping of HIV pol gene was done by in situ DNA hybridization using a Line Probe Assay and by direct sequencing. Viral variants harbouring resistance mutations such as: M41, T69R, K70R, M184V, T215Y in the pol gene were detected in 14% of the subjects. HIV mutants resistant to NRT inhibitors were found in 10 and 20% of patients infected before and after the year 2000, respectively. Multiple drug resistant viruses (2-3 drug classes) were present in 3.5% of the mainly recently infected patients. In protease gene only minor resistant mutations were found such as L101 and A71V. These findings indicate the evolution of drug resistance showing a correlation with the time of introduction of combination therapy in our country, where more than 70% of HIV infections were by homo/bisexual transmission. This confirms the transmission of drug-resistant HIV shown by genotype testing during primary infection in therapy-naive patients and initiates serious clinical and public health consequences.

Modeling HIV infection dynamics: the role of recombination in the development of drug resistance

The benefit of recombination to HIV remains unclear because just as recombination can induce the association of favorable mutations and accelerate the development of multidrug resistance, it can also dissociate favorable combinations of mutations. The confounding influences of mutation, random genetic drift, selection and epistatic interactions between multiple resistance loci render the role of recombination difficult to unravel experimentally. Mathematical models provide valuable insights into the influence of recombination on the genomic diversification of HIV and the development of drug resistance in patients undergoing therapy, capture several recent experimental observations of HIV recombination quantitatively, and set the stage for the establishment of a robust framework for the identification of improved treatment protocols and guidelines for drug development.

Recombination increases human immunodeficiency virus fitness, but not necessarily diversity

Recombination can facilitate the accumulation of mutations and accelerate the emergence of resistance to current antiretroviral therapies for human immunodeficiency virus (HIV) infection. Yet, since recombination can also dissociate favourable combinations of mutations, the benefit of recombination to HIV remains in question. The confounding effects of mutation, multiple infections of cells, random genetic drift and fitness selection that underlie HIV evolution render the influence of recombination difficult to unravel. We developed computer simulations that mimic the genomic diversification of HIV within an infected individual and elucidate the influence of recombination. We find, interestingly, that when the effective population size of HIV is small, recombination increases both the diversity and the mean fitness of the viral population. When the effective population size is large, recombination increases viral fitness but decreases diversity. In effect, recombination enhances (lowers) the likelihood of the existence of multi-drug resistant strains of HIV in infected individuals prior to the onset of therapy when the effective population size is small (large). Our simulations are consistent with several recent experimental observations, including the evolution of HIV diversity and divergencein vivo. The intriguing dependencies on the effective population size appear due to the subtle interplay of drift, selection and epistasis, which we discuss in the light of modern population genetics theories. Current estimates of the effective population size of HIV have large discrepancies. Our simulations present an avenue for accurate determination of the effective population size of HIVin vivoand facilitate establishment of the benefit of recombination to HIV.

Dynamics of a hepatitis B virus e antigen minus population ascribed to genotype F during the course of a chronic infection despite the presence of anti-HBs antibodies

The in vivo evolution of genotype F HBV variants was recorded in a chronically infected patient throughout a 3-year observation period. Fluctuating levels of HBs Ag and anti-HBs antibodies were recorded, both of them cocirculating in peripheral blood samples at given times. Fifty S gene derived clones were sequenced and phylogenetically analyzed. As expected, some amino acid replacements within the S ORF were also observed within the P ORF while others were silent for the former. Such change was statistically significant for both S and P overlapping genes, which clearly indicates the appearance of a positive selection pressure. Supporting this notion, amino acid replacements were documented at both B and T cell epitopes in samples from 1997 and 1998. Several mutations were documented within and outside the "a" determinant in the major hydrophilic region. Such substitutions might have resulted from the attempt of HBV to evade both humoral and/or cellular immune response. To the best of our knowledge this unusual profile of HBV variants in presence of usually "neutralizing" anti-HBs antibodies was examined in vivo for the first time.

A genetic-algorithm approach to simulating human immunodeficiency virus evolution reveals the strong impact of multiply infected cells and recombination

It has been previously shown that the majority of human immunodeficiency virus type 1 (HIV-1)-infected splenocytes can harbour multiple, divergent proviruses with a copy number ranging from one to eight. This implies that, besides point mutations, recombination should be considered as an important mechanism in the evolution of HIV within an infected host. To explore in detail the possible contributions of multi-infection and recombination to HIV evolution, the effects of major microscopic parameters of HIV replication (i.e. the point-mutation rate, the crossover number, the recombination rate and the provirus copy number) on macroscopic characteristics (such as the Hamming distance and the abundance of n-point mutants) have been simulated in silico. Simulations predict that multiple provirus copies per infected cell and recombination act in synergy to speed up the development of sequence diversity. Point mutations can be fixed for some time without fitness selection. The time needed for the selection of multiple mutations with increased fitness is highly variable, supporting the view that stochastic processes may contribute substantially to the kinetics of HIV variation in vivo.

Recombination and migration of Cryphonectria hypovirus 1 as inferred from gene genealogies and the coalescent

Genealogy-based methods were used to estimate migration of the fungal virus Cryphonectria hypovirus 1 between vegetative compatibility types of the host fungus, Cryphonectria parasitica, as a means of estimating horizontal transmission within two host populations. Vegetative incompatibility is a self/non-self recognition system that inhibits virus transmission under laboratory conditions but its effect on transmission in nature has not been clearly demonstrated. Recombination within and among different loci in the virus genome restricted the genealogical analyses to haplotypes with common mutation and recombinational histories. The existence of recombination necessitated that we also use genealogical approaches that can take advantage of both the mutation and recombinational histories of the sample. Virus migration between populations was significantly restricted. In contrast, estimates of migration between vegetative compatibility types were relatively high within populations despite previous evidence that transmission in the laboratory was restricted. The discordance between laboratory estimates and migration estimates from natural populations highlights the challenges in estimating pathogen transmission rates. Genealogical analyses inferred migration patterns throughout the entire coalescent history of one viral region in natural populations and not just recent patterns of migration or laboratory transmission. This application of genealogical analyses provides markedly stronger inferences on overall transmission rates than laboratory estimates do.

Recombination and Migration of Cryphonectria hypovirus 1 as Inferred From Gene Genealogies and the Coalescent

Genealogy-based methods were used to estimate migration of the fungal virus Cryphonectria hypovirus 1 between vegetative compatibility types of the host fungus, Cryphonectria parasitica, as a means of estimating horizontal transmission within two host populations. Vegetative incompatibility is a self/non-self recognition system that inhibits virus transmission under laboratory conditions but its effect on transmission in nature has not been clearly demonstrated. Recombination within and among different loci in the virus genome restricted the genealogical analyses to haplotypes with common mutation and recombinational histories. The existence of recombination necessitated that we also use genealogical approaches that can take advantage of both the mutation and recombinational histories of the sample. Virus migration between populations was significantly restricted. In contrast, estimates of migration between vegetative compatibility types were relatively high within populations despite previous evidence that transmission in the laboratory was restricted. The discordance between laboratory estimates and migration estimates from natural populations highlights the challenges in estimating pathogen transmission rates. Genealogical analyses inferred migration patterns throughout the entire coalescent history of one viral region in natural populations and not just recent patterns of migration or laboratory transmission. This application of genealogical analyses provides markedly stronger inferences on overall transmission rates than laboratory estimates do.

Drug resistance profiles of recombinant reverse transcriptases from human immunodeficiency virus type 1 subtypes A/E, B, and C

We have expressed purified recombinant reverse transcriptase (RT) from clinical isolates of human immunodeficiency virus subtypes B, C, and A/E in Escherichia coli. The drug sensitivities of these RTs were then determined for both nucleoside RT inhibitors (NRTIs) and nonnucleoside RT inhibitors (NNRTIs) in cell-free RT assays. Although A/E and C viruses contained numerous polymorphisms relative to subtype B (i.e., naturally occurring variations unrelated to drug resistance), the wild-type enzymes prepared from these or subtype A/E clinical isolates displayed <2-fold differences in drug sensitivities with regard to the active triphosphate active forms of NRTIs, as compared with RT expressed from BH-10 recombinant virus. Recombinant RTs from clinical isolates of subtypes B, C, and A/E that contained multiple resistance-associated mutations displayed expected variations in levels of resistance to the intracellular active forms of 3TC, ddI, ddC, and PMPA, that is, 3TCTP, ddATP, ddCTP, and PMPApp, respectively. Subtype A/E and C RT enzymes contained only minor NNRTI polymorphisms that distinguished them from wild-type subtype B enzymes and wild-type RTs from these various subtypes showed only 1- to 4-fold variability in IC(50) values for each of nevirapine (NVP), delavirdine (DLV), efavirenz (EFV), and calanolide A. In contrast, RT enzymes from subtype B and C viruses harboring specific NNRTI mutations were highly resistant to all four tested NNRTIs. Subtype C variants containing the novel V106M resistance codon showed cross-resistance to all approved NNRTIs in cell-free RT assays.

Epidemiology, genetic diversity, and evolution of endemic feline immunodeficiency virus in a population of wild cougars

Within the large body of research on retroviruses, the distribution and evolution of endemic retroviruses in natural host populations have so far received little attention. In this study, the epidemiology, genetic diversity, and molecular evolution of feline immunodeficiency virus specific to cougars (FIVpco) was examined using blood samples collected over several years from a free-ranging cougar population in the western United States. The virus prevalence was 58% in this population (n = 52) and increased significantly with host age. Based on phylogenetic analysis of fragments of envelope (env) and polymerase (pol) genes, two genetically distinct lineages of FIVpco were found to cooccur in the population but not in the same individuals. Within each of the virus lineages, geographically nearby isolates formed monophyletic clusters of closely related viruses. Sequence diversity for env within a host rarely exceeded 1%, and the evolution of this gene was dominated by purifying selection. For both pol and env, our data indicate mean rates of molecular evolution of 1 to 3% per 10 years. These results support the premise that FIVpco is well adapted to its cougar host and provide a basis for comparing lentivirus evolution in endemic and epidemic infections in natural hosts.

Turnover of env variable region 1 and 2 genotypes in subjects with late-stage human immunodeficiency virus type 1 infection

The env gene of human immunodeficiency virus type 1 (HIV-1) includes some of the most genetically diverse regions of the viral genome, which are called variable regions 1 through 5 (V1 through V5). We have developed a heteroduplex tracking assay to detect changes in variable regions 1 and 2 of env (V1/V2-HTA). Using sequences from two molecular clones as probes, we have studied the nature of longitudinal virus population changes in a cohort of HIV-1-infected subjects. Viral sequences present in 21 subjects with late-stage HIV-1 infection were initially screened for stability of the virus population by V1/V2-HTA. The virus populations at entry comprised an average of five coexisting V1/V2 genotypic variants (as identified by HTA). Eight of the 21 subjects were examined in detail because of the dynamic behavior of their env variants over an approximately 9-month period. In each of these cases we detected a single discrete transition of V1/V2 genotypes based on monthly sampling. The major V1/V2 genotypes (those present at >10% abundance) from the eight subjects were cloned and sequenced to define the nature of V1/V2 variability associated with a discrete transition. Based on a comparison of V1/V2 genotypic variants present at entry with the newly emerged variants we categorized the newly emerged variants into two groups: variants without length differences and variants with length differences. Variants without length differences had fewer nucleotide substitutions, with the changes biased to either V1 or V2, suggestive of recent evolutionary events. Variants with length differences included ones with larger numbers of changes that were distributed, suggestive of recall of older genotypes. Most length differences were located in domains where the codon motif AVT (V = A, G, C) had become enriched and fixed. Finally, recombination events were detected in two subjects, one of which resulted in the reassortment of V1 and V2 regions. We suggest that turnover in V1/V2 populations was largely driven by selection on either V1 or V2 and that escape was accomplished either through changes focused in the region under selection or by the appearance of a highly divergent variant.

Network analysis of human and simian immunodeficiency virus sequence sets reveals massive recombination resulting in shorter pathways

The intrinsic recombination rate of human immunodeficiency virus (HIV) exceeds the point mutation rate by a factor of 10. As the majority of infected cells in vivo harbour multiple proviruses, the stage is set for rampant recombination. Therefore, it may be presumed that phylogenic relationships and mutation frequencies will probably be affected by recombination. However, the proportion of homoplasies arising from recombination and mutation is not known. By studying the evolution of the hypervariable regions of the simian immunodeficiency virus envelope gene among four macaques, it is shown that homoplasies arise more from recombination than from point mutation. When recombination is accounted for, the minimum number of substitutions in a sequence set may be reduced by as much as 45 %. In fact, the true number of point mutations in a set of HIV sequences tends to the number of discrete substitutions. Hence, lineages are younger than anticipated previously, although not in proportion to the ratio of the intrinsic recombination/point mutation rate. Recombination also inflates codon polymorphisms.

Intrapatient diversity and its correlation with viral setpoint in human immunodeficiency virus type 1 CRF02_A/G-IbNG infection

The human immunodeficiency virus type 1 (HIV-1) viral setpoint during the disease-free interval has been strongly associated with future risk of disease progression. An awareness of the correlation between viral setpoint and HIV-1 genetic evolution over time is important in the understanding of viral dynamics and infection. We examined genetic diversity in HIV-1 CRF02_A/G-IbNG-infected seroincident women in Dakar, Senegal; determined whether a viral setpoint kinetic pattern existed for CRF02_A/G-IbNG during the disease-free interval; and correlated viral load level and diversity. Samples were drawn during the disease-free interval from consenting CRF02_A/G-IbNG-infected, antiretroviral therapy-naïve female commercial sex workers in Dakar, Senegal. Based on sequential plasma RNA values, low and high viral setpoint groups were established. Intrapatient diversity and divergence over time was determined from earlier and later time point DNA samples from each person. Most individuals followed the viral setpoint paradigm. For each 1/-/log(10) copy/ml of plasma increase in viral load, intrapatient diversity increased by 1.4% (P = 0.028). A greater diversification rate was observed in the high viral setpoint group than in the low viral setpoint group (P = 0.01). Greater nucleotide (P = 0.015) and amino acid (P = 0.048) divergences and a greater nucleotide divergence rate (P = 0.03) were found in the high viral setpoint group. There was no difference between the groups in the ratio of the number of nonsynonymous substitutions per nonsynonymous site to the number of synonymous substitutions per synonymous site. The greater intrapatient diversity, divergence, and diversification rates observed in the high viral setpoint group supports the notion that diversity is driven by cycles of viral replication resulting in accumulated mutations. Recognizing diversity potential based on viral load levels in individuals may inform the design of vaccines and therapies.

The application of molecular phylogenetics to the analysis of viral genome diversity and evolution

DNA sequencing and molecular phylogenetics are increasingly being used in virology laboratories to study the transmission of viruses. By reconstructing the evolutionary history of viral genomes the behaviour of viral populations can be modelled, and the future of epidemics may be forecast. The manner in which such viral DNA sequences are analysed is the focus of this review. Many researchers resort to the often-quoted 'black box' approach because phylogenetics theory can be daunting, and phylogenetics software packages can appear to be difficult to use. However, because phylogenetic analyses are often used in important and sensitive arenas, for example to provide evidence indicating transmission between persons, it is vital that appropriate care is taken to estimate reliably true relationships. In this review, we discuss how a molecular phylogenetics study should be approached, give an overview of the methods and programs for analysing DNA sequence data, and point readers to appropriate texts for further details. The aim of this review, therefore, is to provide researchers with an easy to understand guide to molecular phylogenetics, with special reference to viral genomes.

Insertion/deletion frequencies match those of point mutations in the hypervariable regions of the simian immunodeficiency virus surface envelope gene

A method for encoding insertions and deletions (indels) has been developed and adapted to the SplitsTree program. Following phylogenetic reconstruction, the relative frequencies of indels were estimated for a large number of in vivo sequence sets corresponding to the env V1 hypervariable region of the simian immunodeficiency virus SIVmac251. The method allowed recovery of many point mutations hitherto lost due to gap stripping. Deletions were as frequent as transversions and were 4- to 8-fold more frequent than insertions, invariably duplications. The high proportion of deletions among mutation events suggests that lentivirus vectors may readily delete parts of their cargo.