Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease.

The rate of development of disease varies considerably among human immunodeficiency virus type 1 (HIV-1)-infected children. The reasons for these observed differences are not clearly understood but most probably depend on the dynamic interplay between the HIV-1 quasispecies virus population and the immune constraints imposed by the host. To study the relationship between disease progression and genetic diversity, we analyzed the evolution of viral sequences within six perinatally infected children by examining proviral sequences spanning the C2 through V5 regions of the viral envelope gene by PCR of blood samples obtained at sequential visits. PCR product DNAs from four sample time points per child were cloned, and 10 to 13 clones from each sample were sequenced. Greater genetic distances relative to the time of infection were found for children with low virion-associated RNA burdens and slow progression to disease relative to those found for children with high virion-associated RNA burdens and rapid progression to disease. The greater branch lengths observed in the phylogenetic reconstructions correlated with a higher accumulation rate of nonsynonymous base substitutions per potential nonsynonymous site, consistent with positive selection for change rather than a difference in replication kinetics. Viral sequences from children with slow progression to disease also showed a tendency to form clusters that associated with different sampling times. These progressive shifts in the viral population were not found in viral sequences from children with rapid progression to disease. Therefore, despite the HIV-1 quasispecies being a diverse, rapidly evolving, and competing population of genetic variants, different rates of genetic evolution could be found under different selective constraints. These data suggest that the evolutionary dynamics exhibited by the HIV-1 quasispecies virus populations are compatible with a Darwinian system evolving under the constraints of natural selection.     

The threshold bootstrap clustering: a new approach to find families or transmission clusters within molecular quasispecies

Background

Phylogenetic methods produce hierarchies of molecular species, inferring knowledge about taxonomy and evolution. However, there is not yet a consensus methodology that provides a crisp partition of taxa, desirable when considering the problem of intra/inter-patient quasispecies classification or infection transmission event identification. We introduce the threshold bootstrap clustering (TBC), a new methodology for partitioning molecular sequences, that does not require a phylogenetic tree estimation.

Methodology/Principal Findings

The TBC is an incremental partition algorithm, inspired by the stochastic Chinese restaurant process, and takes advantage of resampling techniques and models of sequence evolution. TBC uses as input a multiple alignment of molecular sequences and its output is a crisp partition of the taxa into an automatically determined number of clusters. By varying initial conditions, the algorithm can produce different partitions. We describe a procedure that selects a prime partition among a set of candidate ones and calculates a measure of cluster reliability. TBC was successfully tested for the identification of type-1 human immunodeficiency and hepatitis C virus subtypes, and compared with previously established methodologies. It was also evaluated in the problem of HIV-1 intra-patient quasispecies clustering, and for transmission cluster identification, using a set of sequences from patients with known transmission event histories.

Conclusion

TBC has been shown to be effective for the subtyping of HIV and HCV, and for identifying intra-patient quasispecies. To some extent, the algorithm was able also to infer clusters corresponding to events of infection transmission. The computational complexity of TBC is quadratic in the number of taxa, lower than other established methods; in addition, TBC has been enhanced with a measure of cluster reliability. The TBC can be useful to characterise molecular quasipecies in a broad context.     

Phylogenetic surveillance of viral genetic diversity and the evolving molecular epidemiology of human immunodeficiency virus type 1

With ongoing generation of viral genetic diversity and increasing levels of migration, the global human immunodeficiency virus type 1 (HIV-1) epidemic is becoming increasingly heterogeneous. In this study, we investigate the epidemiological characteristics of 5,675 HIV-1 pol gene sequences sampled from distinct infections in the United Kingdom. These sequences were phylogenetically analyzed in conjunction with 976 complete-genome and 3,201 pol gene reference sequences sampled globally and representing the broad range of HIV-1 genetic diversity, allowing us to estimate the probable geographic origins of the various strains present in the United Kingdom. A statistical analysis of phylogenetic clustering in this data set identified several independent transmission chains within the United Kingdom involving recently introduced strains and indicated that strains more commonly associated with infections acquired heterosexually in East Africa are spreading among men who have sex with men. Coalescent approaches were also used and indicated that the transmission chains that we identify originated in the late 1980s to early 1990s. Similar changes in the epidemiological structuring of HIV epidemics are likely to be taking in place in other industrialized nations with large immigrant populations. The framework implemented here takes advantage of the vast amount of routinely generated HIV-1 sequence data and can provide epidemiological insights not readily obtainable through standard surveillance methods.     

In-depth characterization of viral isolates from plasma and cells compared with plasma circulating quasispecies in early HIV-1 infection

Background

The use of in vitro models to unravel the phenotypic characteristics of circulating viral variants is key to understanding HIV-1 pathogenesis but limited by the availability of primary viral isolates from biological samples. However, overall in vivo genetic variability of HIV-1 within a subject may not be reflected in the viable viral population obtained after isolation. Although several studies have tried to determine whether viral populations expanded in vitro are representative of in vivo findings, the answer remains unclear due to the reduced number of clonal sequences analyzed or samples compared. In order to overcome previous experimental limitations, here we applied Deep Pyrosequencing (DPS) technology in combination with phenotypic experiments to analyze and compare with unprecedented detail the composition of viral isolates and in vivo quasispecies.

Methodology/Principal Findings

We amplified by DPS HIV-1 genomic regions covering gag, protease, integrase and env-V3 to characterize paired isolates from plasma and peripheral blood mononuclear cells and compare them with total plasma viral RNA in four recently HIV-1 infected subjects. Our study demonstrated the presence of unique haplotypes scattered between sample types with conservation of major variants. In addition, no differences in intra- and inter-population encoded protein variability were found between the different types of isolates or when these were compared to plasma viral RNA within subjects. Additionally, in vitro experiments demonstrated phenotypic similarities in terms of replicative capacity and co-receptor usage between viral isolates and plasma viral RNA.

Conclusion

This study is the first in-depth comparison and characterization of viral isolates from different sources and plasma circulating quasispecies using DPS in recently HIV-1 infected subjects. Our data supports the use of primary isolates regardless of their plasma or cellular origin to define genetic variability and biological traits of circulating HIV-1 quasispecies.     

利用单基因组扩增法对马传染性贫血病毒疫苗株异质性的分析

前期研究发现,马传染性贫血病毒（Equine infectious anemia virus EIAV）中国弱毒疫苗株并非单一病毒,而是由多种准种（quasispecies）组成的种群。阐明该疫苗株的具体构成,对于确定优势疫苗株和分析其在体内的进化具有重要意义。本研究比较了传统RNA病毒测序法（即bulk PCR）和单基因组扩增法（Single-genome Amplification, SGA）在扩增EIAV疫苗株囊膜表面蛋白gp90基因 V3~V5区序列上的差异。结果发现,利用SGA法和bulk PCR法获得的序列在组内差异率分别为1.84%和1.88%。进一步序列比较发现,SGA法扩增的序列中除了含有与bulk PCR法中同源性较高的序列外,还存在bulk PCR法未检出的含强毒株LN40特异性位点,以及单个氨基酸缺失的序列。上述序列的存在为该疫苗株“多克隆构成”假说提供了佐证。此外,在对抽样偏差分析中发现,由于疫苗株中各种病毒准种在量上的差异,使得传统bulk PCR法不能有效的扩增组成比例较低的病毒准种,而导致测得的序列组成不能完全代表实际情况。SGA法通过对单基因组分的扩增和测序,可避免bulk PCR法的以上缺陷,在分析以准种形式存在的RNA病毒序列方面具有独特的优势。     

QSdpR: Viral quasispecies reconstruction via correlation clustering

RNA viruses are characterized by high mutation rates that give rise to populations of closely related genomes, known as viral quasispecies. Underlying heterogeneity enables the quasispecies to adapt to changing conditions and proliferate over the course of an infection. Determining genetic diversity of a virus (i.e., inferring haplotypes and their proportions in the population) is essential for understanding its mutation patterns, and for effective drug developments. Here, we present QSdpR, a method and software for the reconstruction of quasispecies from short sequencing reads. The reconstruction is achieved by solving a correlation clustering problem on a read-similarity graph and the results of the clustering are used to estimate frequencies of sub-species; the number of sub-species is determined using pseudo F index. Extensive tests on both synthetic datasets and experimental HIV-1 and Zika virus data demonstrate that QSdpR compares favorably to existing methods in terms of various performance metrics.     

Topology of evolving,mutagenized viral populations: quasispecies expansion,compression, and operation of negative selection

Background  

The molecular events and evolutionary forces underlying lethal mutagenesis of virus (or virus extinction through an excess of mutations) are not well understood. Here we apply for the first time phylogenetic methods and Partition Analysis of Quasispecies (PAQ) to monitor genetic distances and intra-population structures of mutant spectra of foot-and-mouth disease virus (FMDV) quasispecies subjected to mutagenesis by base and nucleoside analogues.     

Analysis of the EIAV Rev-responsive element (RRE) reveals a conserved RNA motif required for high affinity Rev binding in both HIV-1 and EIAV

A cis-acting RNA regulatory element, the Rev-responsive element (RRE), has essential roles in replication of lentiviruses, including human immunodeficiency virus (HIV-1) and equine infection anemia virus (EIAV). The RRE binds the viral trans-acting regulatory protein, Rev, to mediate nucleocytoplasmic transport of incompletely spliced mRNAs encoding viral structural genes and genomic RNA. Because of its potential as a clinical target, RRE-Rev interactions have been well studied in HIV-1; however, detailed molecular structures of Rev-RRE complexes in other lentiviruses are still lacking. In this study, we investigate the secondary structure of the EIAV RRE and interrogate regulatory protein-RNA interactions in EIAV Rev-RRE complexes. Computational prediction and detailed chemical probing and footprinting experiments were used to determine the RNA secondary structure of EIAV RRE-1, a 555 nt region that provides RRE function in vivo. Chemical probing experiments confirmed the presence of several predicted loop and stem-loop structures, which are conserved among 140 EIAV sequence variants. Footprinting experiments revealed that Rev binding induces significant structural rearrangement in two conserved domains characterized by stable stem-loop structures. Rev binding region-1 (RBR-1) corresponds to a genetically-defined Rev binding region that overlaps exon 1 of the EIAV rev gene and contains an exonic splicing enhancer (ESE). RBR-2, characterized for the first time in this study, is required for high affinity binding of EIAV Rev to the RRE. RBR-2 contains an RNA structural motif that is also found within the high affinity Rev binding site in HIV-1 (stem-loop IIB), and within or near mapped RRE regions of four additional lentiviruses. The powerful integration of computational and experimental approaches in this study has generated a validated RNA secondary structure for the EIAV RRE and provided provocative evidence that high affinity Rev binding sites of HIV-1 and EIAV share a conserved RNA structural motif. The presence of this motif in phylogenetically divergent lentiviruses suggests that it may play a role in highly conserved interactions that could be targeted in novel anti-lentiviral therapies.     

Analysis of correlated mutations in HIV-1 protease using spectral clustering

MOTIVATION: The ability of human immunodeficiency virus-1 (HIV-1) protease to develop mutations that confer multi-drug resistance (MDR) has been a major obstacle in designing rational therapies against HIV. Resistance is usually imparted by a cooperative mechanism that can be elucidated by a covariance analysis of sequence data. Identification of such correlated substitutions of amino acids may be obscured by evolutionary noise. RESULTS: HIV-1 protease sequences from patients subjected to different specific treatments (set 1), and from untreated patients (set 2) were subjected to sequence covariance analysis by evaluating the mutual information (MI) between all residue pairs. Spectral clustering of the resulting covariance matrices disclosed two distinctive clusters of correlated residues: the first, observed in set 1 but absent in set 2, contained residues involved in MDR acquisition; and the second, included those residues differentiated in the various HIV-1 protease subtypes, shortly referred to as the phylogenetic cluster. The MDR cluster occupies sites close to the central symmetry axis of the enzyme, which overlap with the global hinge region identified from coarse-grained normal-mode analysis of the enzyme structure. The phylogenetic cluster, on the other hand, occupies solvent-exposed and highly mobile regions. This study demonstrates (i) the possibility of distinguishing between the correlated substitutions resulting from neutral mutations and those induced by MDR upon appropriate clustering analysis of sequence covariance data and (ii) a connection between global dynamics and functional substitution of amino acids.     

An exploratory algorithm to identify intra-host recombinant viral sequences

Since recombination leads to the generation of mosaic genomes that violate the assumption of traditional phylogenetic methods that sequence evolution can be accurately described by a single tree, results and conclusions based on phylogenetic analysis of data sets including recombinant sequences can be severely misleading. Many methods are able to adequately detect recombination between diverse sequences, for example between different HIV-1 subtypes. More problematic is the identification of recombinants among closely related sequences such as a viral population within a host. We describe a simple algorithmic procedure that enables detection of intra-host recombinants based on split-decomposition networks and a robust statistical test for recombination. By applying this algorithm to several published HIV-1 data sets we conclude that intra-host recombination was significantly underestimated in previous studies and that up to one-third of the env sequences longitudinally sampled from a given subject can be of recombinant origin. The results show that our procedure can be a valuable exploratory tool for detection of recombinant sequences before phylogenetic analysis, and also suggest that HIV-1 recombination in vivo is far more frequent and significant than previously thought.     

Variability of human immunodeficiency virus type 1 group O strains isolated from Cameroonian patients living in France.

Human immunodeficiency virus type 1 (HIV-1) nucleotide sequences encoding p24Gag and the Env C2V3 region were obtained from seven patients who were selected on the basis of having paradoxical seronegativity on a subset of HIV enzyme-linked immunosorbent assay detection kits and having atypical Western blot (immunoblot) reactivity. Sequence analyses showed that all of these strains were more closely related to the recently described Cameroonian HIV isolates of group O (HIV-1 outlier) than to group M (HIV-1 major). All seven patients had Cameroonian origins but were living in France at the time the blood samples were taken. Characterization of a large number of group M strains has to date revealed eight distinct genetic subtypes (A to H). Genetic distances between sequences from available group O isolates were generally comparable to those observed in M intersubtype sequence comparisons, showing that the group O viruses are genetically very diverse. Analysis of sequences from these seven new viral strains, combined with the three previously characterized group O strains, revealed few discernable phylogenetic clustering patterns among the 10 patients' viral sequences. The level of diversity among group O sequences suggests that they may have a comparable (or greater) age than the M group sequences, although for unknown reasons, the latter group dispersed first and is the dominant lineage in the pandemic.     

Selecting models of nucleotide substitution: an application to human immunodeficiency virus 1 (HIV-1)

The blind use of models of nucleotide substitution in evolutionary analyses is a common practice in the viral community. Typically, a simple model of evolution like the Kimura two-parameter model is used for estimating genetic distances and phylogenies, either because other authors have used it or because it is the default in various phylogenetic packages. Using two statistical approaches to model fitting, hierarchical likelihood ratio tests and the Akaike information criterion, we show that different viral data sets are better explained by different models of evolution. We demonstrate our results with the analysis of HIV-1 sequences from a hierarchy of samples; sequences within individuals, individuals within subtypes, and subtypes within groups. We also examine results for three different gene regions: gag, pol, and env. The Kimura two-parameter model was not selected as the best-fit model for any of these data sets, despite its widespread use in phylogenetic analyses of HIV-1 sequences. Furthermore, the model complexity increased with increasing sequence divergence. Finally, the molecular-clock hypothesis was rejected in most of the data sets analyzed, throwing into question clock-based estimates of divergence times for HIV-1. The importance of models in evolutionary analyses and their repercussions on the derived conclusions are discussed.     

Simple epidemiological dynamics explain phylogenetic clustering of HIV from patients with recent infection

Phylogenies of highly genetically variable viruses such as HIV-1 are potentially informative of epidemiological dynamics. Several studies have demonstrated the presence of clusters of highly related HIV-1 sequences, particularly among recently HIV-infected individuals, which have been used to argue for a high transmission rate during acute infection. Using a large set of HIV-1 subtype B pol sequences collected from men who have sex with men, we demonstrate that virus from recent infections tend to be phylogenetically clustered at a greater rate than virus from patients with chronic infection ('excess clustering') and also tend to cluster with other recent HIV infections rather than chronic, established infections ('excess co-clustering'), consistent with previous reports. To determine the role that a higher infectivity during acute infection may play in excess clustering and co-clustering, we developed a simple model of HIV infection that incorporates an early period of intensified transmission, and explicitly considers the dynamics of phylogenetic clusters alongside the dynamics of acute and chronic infected cases. We explored the potential for clustering statistics to be used for inference of acute stage transmission rates and found that no single statistic explains very much variance in parameters controlling acute stage transmission rates. We demonstrate that high transmission rates during the acute stage is not the main cause of excess clustering of virus from patients with early/acute infection compared to chronic infection, which may simply reflect the shorter time since transmission in acute infection. Higher transmission during acute infection can result in excess co-clustering of sequences, while the extent of clustering observed is most sensitive to the fraction of infections sampled.     

Predicting demographic group structures based on DNA sequence data

The ability to infer relationships between groups of sequences, either by searching for their evolutionary history or by comparing their sequence similarity, can be a crucial step in hypothesis testing. Interpreting relationships of human immunodeficiency virus type 1 (HIV-1) sequences can be challenging because of their rapidly evolving genomes, but it may also lead to a better understanding of the underlying biology. Several studies have focused on the evolution of HIV-1, but there is little information to link sequence similarities and evolutionary histories of HIV-1 to the epidemiological information of the infected individual. Our goal was to correlate patterns of HIV-1 genetic diversity with epidemiological information, including risk and demographic factors. These correlations were then used to predict epidemiological information through analyzing short stretches of HIV-1 sequence. Using standard phylogenetic and phenetic techniques on 100 HIV-1 subtype B sequences, we were able to show some correlation between the viral sequences and the geographic area of infection and the risk of men who engage in sex with men. To help identify more subtle relationships between the viral sequences, the method of multidimensional scaling (MDS) was performed. That method identified statistically significant correlations between the viral sequences and the risk factors of men who engage in sex with men and individuals who engage in sex with injection drug users or use injection drugs themselves. Using tree construction, MDS, and newly developed likelihood assignment methods on the original 100 samples we sequenced, and also on a set of blinded samples, we were able to predict demographic/risk group membership at a rate statistically better than by chance alone. Such methods may make it possible to identify viral variants belonging to specific demographic groups by examining only a small portion of the HIV-1 genome. Such predictions of demographic epidemiology based on sequence information may become valuable in assigning different treatment regimens to infected individuals.     

Genetic analysis of human immunodeficiency virus type 1 strains from patients in Cyprus: identification of a new subtype designated subtype I.

DNA sequences encoding the C2 to V3 region of envelope glycoprotein gp120 of human immunodeficiency virus type 1 (HIV-1) were amplified by PCR from uncultured peripheral blood mononuclear cells obtained from 24 of 25 HIV-1-seropositive patients from Cyprus. By using a heteroduplex mobility assay (HMA), all amplified products were studied genetically and compared with 16 previously characterized HIV-1 strains belonging to subtypes A through F. HMA results revealed that HIV-1 gp120 sequences from 15 of our patients were of subtype B of HIV-1, whereas one isolate was of subtype C. However, gp120 sequences from eight patients had no obvious similarities to the known subtypes as defined by HMA. DNA sequencing and phylogenetic analyses of molecular clones confirmed the HMA results and placed the eight undefined HIV-1 isolates into three distinct genetic clusters. On the basis of branch topology and lengths of the phylogenetic tree, we conclude that one group consisting of three clones from two patients represents a new HIV-1 env subtype, which we have termed subtype I. The remaining two sequence clusters, consisting of five sequences from four patients and two sequences from two other patients, are distally related to subtypes A and F. These data demonstrate the extensive heterogeneity of HIV-1 in Cyprus, including the presence of new subtype.     

Analysis of the Temporal Relationship between Human Immunodeficiency Virus Type 1 Quasispecies in Sequential Blood Samples and Various Organs Obtained at Autopsy

We studied the temporal relationship between human immunodeficiency type 1 (HIV-1) quasispecies in tissues and in peripheral blood mononuclear cells (PBMC) of infected individuals. Sequential PBMC and tissue samples from various organs obtained at autopsy from three patients who died of AIDS-related complications were available for analysis. Biological HIV-1 clones were isolated from PBMC samples, and cellular tropism and syncytium-inducing (SI) capacity were determined. Genomic DNA was isolated from 1 cm³ of organ tissue, and proviral DNA was amplified by means of PCR and cloned with the PGEM-T vector system. A 185-bp region encompassing the third variable domain of the virus envelope, known to influence HIV-1 biological properties, was sequenced. HIV-1 could be amplified from all PBMC and organ samples, except from liver tissue for two patients. Both SI and non-syncytium-inducing (NSI) genotypes could be detected in the different tissues. Tissue-specific quasispecies were observed in brain, lung, and testis. Lymphoid tissues, such as bone marrow, lymph node, and spleen, harbored several different variants similar to those detected in blood in the last PBMC samples. In general, only tissues in which macrophages are likely to be the main target cell for HIV-1 harbored NSI HIV-1 sequences that clustered separately. Both SI and NSI sequences that clustered with sequences from late-stage PBMC were present in other tissues, which may indicate that the presence of HIV-1 in those tissues is secondary to lymphocyte infiltration rather than to tissue tropism of HIV-1 itself. These data suggest that the viral reservoir may be limited, which will have important implications for the success of HIV-1 eradication.     

Defective accessory genes in a human immunodeficiency virus type 1-infected long-term survivor lacking recoverable virus.

We have been studying a patient who acquired human immunodeficiency virus (HIV) infection via a blood transfusion 13 years ago. She has remained asymptomatic since that time. The blood donor and two other recipients have all died of AIDS. Although this patient has shown persistently strong seroreactivity to HIV type 1 (HIV-1) antigens by Western blot (immunoblot), she has been continually HIV culture negative in results from multiple laboratories over the last 6 years and has a very low viral burden. Her CD4+ T-cell count has fluctuated around a mean of 399 cells per microliters, with little change in lymphocyte subset percentages. Strong cellular immune responses to HIV-1 epitopes by this patient have been demonstrated. We now report the results of an intensive molecular genetic analysis of the HIV-1 proviral quasispecies from this patient sampled over 5 years. Long terminal repeat region sequences supported the argument for normal basal and Tat-mediated promoter activities. Sequential sequencing of the nef gene revealed a low frequency (8.3%) of defective genes and a striking lack of sequence evolution. Functional analysis of predominant nef genes by both a cell surface CD4 downregulation and a viral infectivity complementation assay showed wild-type function. In contrast, sequential analysis of an amplicon containing the vif, vpr, vpu, tat1, and rev1 genes revealed the presence of inactivating mutations in 64% of the clones. These data suggest that this patient, initially infected with a virulent swarm of HIV-1, is presently infected with a more-attenuated viral quasispecies as a result of effective host immunity.     

HIV-1 Tropism Dynamics and Phylogenetic Analysis from Longitudinal Ultra-Deep Sequencing Data of CCR5- and CXCR4-Using Variants

Objective

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. The molecular and evolutionary mechanisms underlying the CCR5 to CXCR4 switch are the focus of intense recent research. We studied the HIV-1 tropism dynamics in relation to coreceptor usage, the nature of quasispecies from ultra deep sequencing (UDPS) data and their phylogenetic relationships.

Methods

Here, we characterized C2-V3-C3 sequences of HIV obtained from 19 patients followed up for 54 to 114 months using UDPS, with further genotyping and phylogenetic analysis for coreceptor usage. HIV quasispecies diversity and variability as well as HIV plasma viral load were measured longitudinally and their relationship with the HIV coreceptor usage was analyzed. The longitudinal UDPS data were submitted to phylogenetic analysis and sampling times and coreceptor usage were mapped onto the trees obtained.

Results

Although a temporal viral genetic structuring was evident, the persistence of several viral lineages evolving independently along the infection was statistically supported, indicating a complex scenario for the evolution of viral quasispecies. HIV X4-using variants were present in most of our patients, exhibiting a dissimilar inter- and intra-patient predominance as the component of quasispecies even on antiretroviral therapy. The viral populations from some of the patients studied displayed evidences of the evolution of X4 variants through fitness valleys, whereas for other patients the data favored a gradual mode of emergence.

Conclusions

CXCR4 usage can emerge independently, in multiple lineages, along the course of HIV infection. The mode of emergence, i.e. gradual or through fitness valleys seems to depend on both virus and patient factors. Furthermore, our analyses suggest that, besides becoming dominant after population-level switches, minor proportions of X4 viruses might exist along the infection, perhaps even at early stages of it. The fate of these minor variants might depend on both viral and host factors.     

Genetic heterogeneity of HIV-1 in Greece

The aim of this study was to detect and determine the genetic variation of HIV-1 in Greece and to analyze the phylogenetic relationships and transmission dynamics of identified variants. Eighty-six blood samples from HIV-1 seroconverted patients of different risk groups were collected from the AIDS clinic, AHEPA Hospital, Thessaloniki, Greece. Retroviral DNA was extracted from uncultured peripheral blood mononuclear cells. HIV-1 DNA sequences encoding a 500-bp fragment of the gp120 C2-C3 region were amplified from each study subject, and they were genetically subtyped by heteroduplex mobility assay and DNA sequencing. Genetic distances and phylogenetic relationships of DNA sequences were estimated using PHYLIP software. Our results revealed that 82 out of 86 (95.3%) subjects carried subtype B sequences, while four (4.7%) carried subtype A sequences. Subtype A in Greek individuals not having traveled abroad was documented. An average of intrasubtype B genetic divergence of 15% was noted. Our findings demonstrate the presence of at least two genetic subtypes of HIV-1 in northern Greece--subtype B and subtype A. The predominant subtype is subtype B, which was transmitted into Greece by multiple sources. Our observations lend support to the argument that the distribution of HIV-1 subtypes is determined by founder effects or other processes rather than any tropism for particular cell types or mode of transmission.     

Equine infectious anemia virus utilizes a YXXL motif within the late assembly domain of the Gag p9 protein.

We have previously demonstrated that the Gag p9 protein of equine infectious anemia virus (EIAV) is functionally homologous with Rous sarcoma virus (RSV) p2b and human immunodeficiency virus type 1 (HIV-1) p6 in providing a critical late assembly function in RSV Gag-mediated budding from transfected COS-1 cells (L. J. Parent et al., J. Virol. 69:5455-5460, 1995). In light of the absence of amino acid sequence homology between EIAV p9 and the functional homologs of RSV and HIV-1, we have now designed an EIAV Gag-mediated budding assay to define the late assembly (L) domain peptide sequences contained in the EIAV p9 protein. The results of these particle budding assays revealed that expression of EIAV Gag polyprotein in COS-1 cells yielded extracellular Gag particles with a characteristic density of 1.18 g/ml, while expression of EIAV Gag polyprotein lacking p9 resulted in a severe reduction in the release of extracellular Gag particles. The defect in EIAV Gag polyprotein particle assembly could be corrected by substituting either the RSV p2b or HIV-1 p6 protein for EIAV p9. These observations demonstrated that the L domains of EIAV, HIV-1, and RSV were interchangeable in mediating assembly of EIAV Gag particles in the COS-1 cell budding assay. To localize the L domain of EIAV p9, we next assayed the effects of deletions and site-specific mutations in the p9 protein on its ability to mediate budding of EIAV Gag particles. Analyses of EIAV Gag constructs with progressive N-terminal or C-terminal deletions of the p9 protein identified a minimum sequence of 11 amino acids (Q20N21L22Y23P24D25L26S27E28I29K30) capable of providing the late assembly function. Alanine scanning studies of this L-domain sequence demonstrated that mutations of residues Y23, P24, and L26 abrogated the p9 late budding function; mutations of other residues in the p9 L domain did not substantially affect the level of EIAV Gag particle assembly. These data indicate that the L domain in EIAV p9 utilizes a YXXL motif which we hypothesize may interact with cellular proteins to facilitate virus particle budding from infected cells.