The evolutionary dynamics of influenza A virus adaptation to mammalian hosts.

Few questions on infectious disease are more important than understanding how and why avian influenza A viruses successfully emerge in mammalian populations, yet little is known about the rate and nature of the virus' genetic adaptation in new hosts. Here, we measure, for the first time, the genomic rate of adaptive evolution of swine influenza viruses (SwIV) that originated in birds. By using a curated dataset of more than 24 000 human and swine influenza gene sequences, including 41 newly characterized genomes, we reconstructed the adaptive dynamics of three major SwIV lineages (Eurasian, EA; classical swine, CS; triple reassortant, TR). We found that, following the transfer of the EA lineage from birds to swine in the late 1970s, EA virus genes have undergone substantially faster adaptive evolution than those of the CS lineage, which had circulated among swine for decades. Further, the adaptation rates of the EA lineage antigenic haemagglutinin and neuraminidase genes were unexpectedly high and similar to those observed in human influenza A. We show that the successful establishment of avian influenza viruses in swine is associated with raised adaptive evolution across the entire genome for many years after zoonosis, reflecting the contribution of multiple mutations to the coordinated optimization of viral fitness in a new environment. This dynamics is replicated independently in the polymerase genes of the TR lineage, which established in swine following separate transmission from non-swine hosts.

Estimating reassortment rates in co-circulating Eurasian swine influenza viruses.

Swine have often been considered as a mixing vessel for different influenza strains. In order to assess their role in more detail, we undertook a retrospective sequencing study to detect and characterize the reassortants present in European swine and to estimate the rate of reassortment between H1N1, H1N2 and H3N2 subtypes with Eurasian (avian-like) internal protein-coding segments. We analysed 69 newly obtained whole genome sequences of subtypes H1N1-H3N2 from swine influenza viruses sampled between 1982 and 2008, using Illumina and 454 platforms. Analyses of these genomes, together with previously published genomes, revealed a large monophyletic clade of Eurasian swine-lineage polymerase segments containing H1N1, H1N2 and H3N2 subtypes. We subsequently examined reassortments between the haemagglutinin and neuraminidase segments and estimated the reassortment rates between lineages using a recently developed evolutionary analysis method. High rates of reassortment between H1N2 and H1N1 Eurasian swine lineages were detected in European strains, with an average of one reassortment every 2-3 years. This rapid reassortment results from co-circulating lineages in swine, and in consequence we should expect further reassortments between currently circulating swine strains and the recent swine-origin H1N1v pandemic strain.

Estimating selection pressures on HIV-1 using phylogenetic likelihood models.

Human immunodeficiency virus (HIV-1) can rapidly evolve due to selection pressures exerted by HIV-specific immune responses, antiviral agents, and to allow the virus to establish infection in different compartments in the body. Statistical models applied to HIV-1 sequence data can help to elucidate the nature of these selection pressures through comparisons of non-synonymous (or amino acid changing) and synonymous (or amino acid preserving) substitution rates. These models also need to take into account the non-independence of sequences due to their shared evolutionary history. We review how we have developed these methods and have applied them to characterize the evolution of HIV-1 in vivo. To illustrate our methods, we present an analysis of compartment-specific evolution of HIV-1 env in blood and cerebrospinal fluid and of site-to-site variation in the gag gene of subtype C HIV-1.

Heterogeneous clearance rates of long-lived lymphocytes infected with HIV: intrinsic stability predicts lifelong persistence.

Viral replication and latently infected cellular reservoirs persist in HIV-infected patients achieving undetectable plasma virus levels with potent antiretroviral therapy. We exploited a predictable drug resistance mutation in the HIV reverse transcriptase to label and track cells infected during defined intervals of treatment and to identify cells replenished by ongoing replication. Decay rates of subsets of latently HIV-infected cells paradoxically decreased with time since establishment, reflecting heterogeneous lymphocyte activation and clearance. Residual low-level replication can replenish cellular reservoirs; however, it does not account for prolonged clearance rates in patients without detectable viremia. In patients receiving potent antiretroviral therapy, the latent pool has a heterogeneous and dynamic composition that comprises a progressively increasing proportion of stable lymphocytes. Eradication will not be achieved with complete inhibition of viral replication alone.

Soluble recombinant HTLV-1 surface glycoprotein competitively inhibits syncytia formation and viral infection of cells.

Efficient entry into, and infection of, human cells by human T-cell leukaemia virus type-1 (HTLV-1) is mediated by the viral envelope glycoproteins, gp46 and gp21. The gp46 surface glycoprotein binds to an as yet unidentified cell surface receptor, thereby, allowing the gp21 transmembrane glycoprotein to initiate fusion of the viral and cellular membranes. In the absence of membrane fusion viral penetration and entry into the host cell cannot occur. The envelope glycoproteins are also a major target for neutralising antibodies and cytotoxic T lymphocytes following a protective immune response, and represent ideal constituents for a recombinant HTLV-1 vaccine. Given the importance of the envelope proteins in HTLV-1 pathogenesis there is increasing interest in obtaining sufficient quantities of these proteins for biochemical, biophysical and biological analyses. We have now developed a system for production of large amounts of a glycosylated and functional form of soluble recombinant gp46 (sRgp46), and have used this recombinant material for analysis of envelope function and receptor binding activity. We find that, the sRgp46 molecules expressed in our system are immunologically indistinguishable from the native virally expressed surface glycoproteins; that sRgp46 binds to T-cells in a dose dependent and saturable manner; and that cell surface binding by sRgp46 can be inhibited by neutralising antibodies. Importantly, we demonstrate that these sRgp46 molecules potently inhibit syncytia formation and viral infection of target cells, and that regions outwith the SU domain of envelope are not required for binding to target cells or for inhibiting membrane fusion. The sRgp46 produced in our study will provide new opportunities to investigate envelope-receptor interactions, and will be of utility in defining the conformationally sensitive antigenic determinants of the HTLV-1 surface glycoprotein.

Mismatched human leukocyte antigen alleles protect against heterosexual HIV transmission.

Genetic variation at the human leukocyte antigen (HLA) loci has been shown to be an important risk factor for progression to HIV disease, but its significance in infection is less well understood. We have investigated its role in HIV transmission in a cohort of individuals at risk for heterosexual infection. Analysis of over 80 individuals revealed that that the degree of concordance at HLA A, B, and DR loci differs significantly between transmitting and nontransmitting couples at risk for heterosexual HIV transmission (p <.02), suggesting that allogeneic immune responses may confer a degree of protection against HIV infection. Analysis of the frequencies of specific alleles at the A, B, and DR loci revealed a significantly higher frequency of HLA DR5 among exposed uninfected individuals, relative to population controls.

Evidence for positive selection driving the evolution of HIV-1 env under potent antiviral therapy.

In HIV-infected individuals treated with potent antiretroviral therapy, viable virus can be isolated from latently infected cells several years into therapy, due to the long life of these cells, ongoing replication replenishing this population, or both. We have analysed the V3 region of the HIV-1 env gene isolated from six patients who have undergone 2 years of potent antiretroviral therapy without frank failure of viral suppression. We show that in two (and possibly three) patients, the sequence changes between baseline virus and virus isolated from infected cells persisting 2 years into infection result from positive selection driving adaptive evolution, occurring either prior to or during therapy. Our analyses suggest low-level replication despite absence of drug resistance due to drug sanctuary sites, or to low-level ongoing replication in the presence of alterations in the selective environment during therapy, perhaps due to a decline in HIV-specific immune responsiveness or changes in target cell pools. In one patient, genetic divergence between baseline plasma and infected cells isolated during therapy may reflect the long half-life of some of these persistent cell populations and the divergence of viral subpopulations that occurred prior to therapy.

Multiple sites in HIV-1 reverse transcriptase associated with virological response to combination therapy.

OBJECTIVE: To determine whether analysis of sequence variation in reverse transcriptase at baseline can explain differences in response to combination antiretroviral therapy. METHODS: Amino acid sequences of reverse transcriptase obtained from baseline isolates from 55 patients included in a trial of zidovudine and didanosine versus zidovudine/didanosine/nevirapine (ACTG241) were analysed. Simple and multiple linear regression were used to determine the relationship between numbers and identity of mutations at baseline and virological response after 8 and 48 weeks. RESULTS: Numbers of baseline zidovudine resistance mutations were predictive of short-term response (week 8). Amino acid identity at position 215 explained > 20% of the variation in response at week 8, but less at week 48. Multiple regression identified the combinations: 215 + 44 and 41 + 202, each of which explained about 30% of the variation in week 8 response. A model incorporating amino acids 214 + 215 + 60 + 202 + baseline viral load explained > 40% of the variation in response at week 48. Unexpectedly, the mutant combination 601 + 215Y/F responded threefold better than 60V + 215Y/F over 48 weeks. CONCLUSIONS: Use of clinical data to analyse virological response to combination therapy has revealed effects of baseline amino acid mutations at sites not previously identified as being important in antiretroviral resistance. Predictors of long-term responses were different from those involved in the short term and may require more complex analysis.

Method used to identify previously undiagnosed infections in the HIV outbreak at Glenochil prison.

Four years after the occurrence of an outbreak of hepatitis B and HIV infection among injecting drug user inmates at Her Majesty's Prison Glenochil in Scotland, a study design was developed to complete the epidemiological account of the HIV outbreak. Our aim was to identify potential cases of (1) HIV transmission not diagnosed during the original outbreak investigation and (2) the source(s) of the outbreak. Scotland's HIV positive case register was searched for matches to a soundexed list of 636 Glenochil inmates imprisoned during January-June 1993. Eight HIV infections that may have been acquired in Glenochil and four possible sources of the outbreak were identified. The second stage of follow-up molecular epidemiological techniques used on stored sera samples from identified individuals is described in the companion paper. Without breach of medical or prisoner confidentiality, indirect and anonymous follow-up has proved possible for the Glenochil inmates.

Completing the molecular investigation into the HIV outbreak at Glenochil prison.

In a molecular investigation into the outbreak of HIV in Glenochil during the first 6 months of 1993, we previously demonstrated that 13 out of the 14 HIV positive inmates were infected with a virtually identical strain, and discounted 2 others as potential sources. Here we investigate a further 8 potential contacts and 4 potential sources which were identified in the companion paper. We were able to examine viral sequence from all but one of these 12 and results have revealed them to be distinct both from each other and the original 14. Thus, despite an intensive follow-up investigation, we have been unable to identify any further HIV infections that might have been part of the 1993 outbreak. It is possible that persons who were infected at that time remain undetected; however this and the companion report strongly suggest that if this were the case the likely numbers would be few.

Evolution of envelope sequences of human immunodeficiency virus type 1 in cellular reservoirs in the setting of potent antiviral therapy.

In human immunodeficiency virus (HIV)-infected patients treated with potent antiretroviral therapy, the persistence of latently infected cells may reflect the long decay half-life of this cellular reservoir or ongoing viral replication at low levels with continuous replenishment of the population or both. To address these possibilities, sequences encompassing the C2 and V3 domains of HIV-1 env were analyzed from virus present in baseline plasma and from viral isolates obtained after 2 years of suppressive therapy in six patients. The presence of sequence changes consistent with evolution was demonstrated for three subjects and correlated with less complete suppression of viral replication, as indicated by the rapidity of the initial virus load decline or the intermittent reappearance of even low levels of detectable viremia. Together, these results provide evidence for ongoing replication. In the remaining three patients, virus recovered after 2 years of therapy was either genotypically contemporary with or ancestral to virus present in plasma 2 years before, indicating that virus recovery had indeed resulted from activation of latently infected cells.

Higher selection pressure from antiretroviral drugs in vivo results in increased evolutionary distance in HIV-1 pol.

We investigated the effect of selection pressures on evolution of HIV-1 pol in 51 patients after switching to a new antiretroviral combination reverse transcriptase (RT) inhibitor therapy. Evolution of the protease (PR) and RT reading frames were analysed separately. Pairwise evolutionary distances (ED) were calculated between sequences from baseline and week 8 and between baseline and week 48 of protocol therapy. ED were calculated for all substitutions and for synonymous and nonsynonymous substitutions separately. At week 8 when HIV RNA reduction (selection pressure) was high, significantly more divergence in pol in both synonymous and nonsynonymous substitutions was found in patients with substantial RNA reduction (strong responders). Separate analyses of PR and RT revealed significantly greater ED in the RT (under selection pressure) of strong compared with nonresponders, whereas divergence between PR genes (not under selection pressure) did not differ in those two groups. Such differential evolution indicates that PR and RT were genetically unlinked and suggests recombination. The rapid increase of ED over the first 8 weeks was followed by only a minimal further rise by week 48, suggesting that selection of preexisting quasispecies accounted for the early changes. A disproportionally high number of synonymous substitutions accounted for the observed divergence and indicated that such genetic changes may not be completely silent.

Molecular investigation into outbreak of HIV in a Scottish prison.

OBJECTIVES: To support already established epidemiological links between inmates of Glenochil prison positive for HIV infection by using molecular techniques and thus provide evidence of the extent of acquisition during a recent outbreak of the disease resulting from needle sharing. To identify possible sources of the outbreak, and to demonstrate the ability of the methodology to make further links beyond the original outbreak. DESIGN: Viral sequences obtained from the blood of HIV positive prisoners previously identified by standard epidemiological methods were compared with each other and with sequences from other Scottish patients. SETTING: Glenochil prison for men, central Scotland. SUBJECTS: Adult inmates and their possible contacts. RESULTS: Phylogenetic analysis of viral sequences in two different genomic regions showed that 13 of the 14 HIV positive prisoners had been infected from a common source. Previous research had shown that six of these had acquired their infection in Glenochil; molecular evidence suggests that more than double this number were infected while incarcerated. Virus from two long term HIV positive patients who were in the prison at the time of the outbreak but who were not identified in the original or subsequent surveys was sufficiently different to make it unlikely that they were the source. A viral sequence from heterosexual transmission from one inmate showed the ability of these techniques to follow the infection through different routes of infection. CONCLUSION: The number of prisoners infected with HIV during the 1993 outbreak within Glenochil prison was more than twice that previously shown. This shows the potential for the spread of bloodborne diseases within prisons by injecting drugs.

Early evolutionary origin of the planktic foraminifera inferred from small subunit rDNA sequence comparisons.

Phylogenetic analysis of five partial planktic foraminiferal small subunit (SSU) ribosomal (r) DNA sequences with representatives of a diverse range of eukaryote, archaebacterial, and eubacterial taxa has revealed that the evolutionary origin of the foraminiferal lineage precedes the rapid eukaryote diversification represented by the "crown" of the eukaryotic tree and probably represents one of the earliest splits among extant free-living aerobic eukaryotes. The foraminiferal rDNA sequences could be clearly separated from known symbionts, commensals, and food organisms. All five species formed a single monophyletic group distinguished from the "crown" group by unique foraminiferal specific insertions as well as considerable nucleotide distance in aligned regions.

The molecular epidemiology of human immunodeficiency virus type 1 in Edinburgh.

Human immunodeficiency virus (HIV) type 1 sequences obtained from HIV-infected persons in different risk groups in Edinburgh were studied to determine the number and origin of virus variants and patterns of virus transmission. Phylogenetic analysis revealed that 12 of 14 hemophiliac patients who had been exposed to a single common batch of factor VIII had closely related gag gene sequences. Sequences from intravenous drug users and patients infected through heterosexual contact formed another distinct group, and 2 other hemophiliacs formed a third group. However, epidemiologic relationships inferred from analysis of the V3 region of the env gene were less conclusive, especially when the V3 loop was taken in isolation. This appears to be due to the length of time since infection and the action of selection, which has favored the independent appearance of similar V3 loop variants.

Extensive linkage disequilibrium in the achaete-scute complex of Drosophila melanogaster.

We have analyzed the level of gametic association between restriction map variants in a sample of 44 X chromosomes from a natural population of Drosophila melanogaster. Of 21 pairwise tests involving 7 restriction map polymorphisms in the yellow-achaete-scute complex, 17 were found to be significant, including some between restriction sites over 80 kb apart. Three-way linkage disequilibria and their variances were also estimated for all 35 three-way comparisons between these loci. Twelve such tests were found to be significant, again spanning distances of up to 80 kb on the restriction map. Only 9 of a possible 128 haplotypes were represented in the sample and 8 of these could be linked together by changes at a single site. The strength of these associations at y-ac-sc is unusual by comparison with studies on other regions of the genome of D. melanogaster, and is consistent with the very low level of recombination which has been reported for the complex. However, our estimate of nucleotide diversity in the region is not significantly different from those made for some other loci in this species.

Complete reversions of a gypsy retrotransposon-induced cut locus mutation in Drosophila melanogaster involving jockey transposon insertions and flanking gypsy sequence deletions.

We have analysed the structures of three phenotypic revertant alleles of a gypsy retrotransposon-induced mutation at the cut locus of Drosophila melanogaster. All three revertants are associated with the insertion of jockey transposons into a common region of gypsy. Two of these alleles are complete reversions to wild type. One complete revertant (ct+D) is derived from a third allele, a partial revertant (ctMRpD) by a deletion of part of the gypsy sequence flanking the jockey transposon. Sequence differences between the jockey elements in ctMRpD and ct+D suggest that this deletion may have been created by the insertion of a second jockey near to the first, followed by recombinational excision of a composite jockey and the region between the two genetic elements. The other complete revertant also carries a deletion of gypsy DNA flanking the jockey insertion. The deleted regions of both complete revertants and the target region for all the jockey insertions contain a repeated sequence that resembles a transcriptional enhancer. The strength of the cut phenotype in these mutants correlates with the proportion of this region remaining near the gypsy transcriptional start site, suggesting that the jockey insertions relieve the gypsy-induced mutation at cut by interfering with a region which is required for the transcriptional competence of gypsy.

Evolutionary relationships of the human immunodeficiency viruses.

The rapid accumulation of nucleotide sequence data on viral genes has allowed, for the first time, the development of detailed phylogenies of viruses based on an objective criterion. This has been demonstrated clearly in the recent analysis of the evolutionary relationships of HIV - the AIDS virus. When first characterized, HIV seemed aberrant and almost unique in many features. Now it is known to be one of a large group of immunodeficiency viruses, which are widely distributed among primates and other mammals.