Evolution of the Sry genes

Existing DNA sequence data on the Sry gene, the mammalian sex- determining locus in the Y chromosome, were analyzed for primates, rodents, and bovids. In all three taxonomic groups, the terminal sequences evolved faster than the HMG (high mobility group) boxes, and this applies both to synonymous (Ks) and nonsynonymous (Ka) nucleotide substitutions. Similar intragenic correlation between synonymous and nonsynonymous substitution rates was not found either in other mammalian genes that contain a conservative box (Sox, Msx) or in the MADS-box genes of plants. The rate of nonsynonymous substitutions exceeds significantly that of synonymous substitutions in the terminal Sry sequences of apes. We did not find good support for the hypothesis that the high evolutionary rate of Sry would be associated with a promiscuous mating system.      

Independent variation and positive selection in env V1 and V2 domains within maternal-infant strains of human immunodeficiency virus type 1 in vivo.

Multiple targets for immune recognition and cellular tropism are localized to the V1 and V2 hypervariable regions in the amino portion of human immunodeficiency virus type 1 (HIV-1) gp120env. We have assessed genetic diversity in env V1 and V2 hypervariable domains in vivo within epidemiologically related strains of HIV-1. Our strategy was to analyze longitudinal samples from two seropositive mothers and multiple children infected by perinatal transmission. Although the V1 and V2 domains are closely linked in the HIV-1 genome, nucleotide sequences in V1 and in V2 evolved independently in maternal-infant viruses in vivo. A high proportion of the nucleotide substitutions would introduce amino acid diversity in V1 and in V2. A significant excess of nonsynonymous over synonymous substitutions was identified in HIV-1 env V1 and V2 peptides in the mothers and in two older children but was not generally apparent in HIV-1 sequences in infants. An excess of nonsynonymous over synonymous substitutions indicated that there is positive selection for independent genetic variation in the V1 and V2 domains in vivo. It is likely that there are host responses to complex determinants in the V1 or V2 hypervariable domain of HIV-1 gp120.     

The explosive human immunodeficiency virus type 1 epidemic among injecting drug users of Kathmandu, Nepal, is caused by a subtype C virus of restricted genetic diversity

An explosive epidemic of human immunodeficiency virus type 1 (HIV-1) has been documented among the injecting drug user population of Kathmandu, Nepal, whose seropositivity rate has risen from 0 to 40% between 1995 and 1997. By using Catrimox to preserve whole-blood RNA at ambient temperature for transportation, HIV-1 envelope V3-V4 sequences were obtained from 36 patients in this group. Analysis of the sequences indicated a homogenous epidemic of subtype C virus, with at least two independent introductions of the virus into the population. Viral diversity was restricted within two transmission subclusters, with the majority of variation occurring in V4. Calculation of the synonymous-to-nonsynonymous mutation ratio (Ks:Ka) across this region showed that significant evolutionary pressure had been experienced during the rapid horizontal spread of the virus in this population, most strongly directed to the region between V3 and V4.     

Recent Evolutionary History of Human Immunodeficiency Virus Type 1 Subtype B: Reconstruction of Epidemic Onset Based on Sequence Distances to the Common Ancestor

We obtained and studied HIV-1 sequences with a known sampling year from three outbreaks of the HIV-1 epidemic: 141 env V3 (270 nt) sampled between 1984 and 1992 and 117 pol prot/RT (804 nt) sequences sampled between 1986 and 1999 from Dutch homosexual men and injecting drug users (IDUs), as well as 77 env V3 sequences sampled between 1983 and 1994 in the United States. Since retrospective serological and/or epidemiological data on these populations are available, providing estimates of the dates of the onset of the HIV-1 epidemics, we had the opportunity to test different phylogenetic models for their accuracy in deriving the recent evolutionary history of HIV-1 subtype B and the onset date of the HIV-1 epidemic. We observed that, in any given year, individual sequences vary widely in their distances to the common ancestor, and sequences close to the ancestors were found decades after the onset of the epidemic. Nevertheless, the mean evolutionary distances of virus strains to ancestors were increasing significantly during the course of the studied epidemics, which indicates that the molecular clock is operational in the recent evolution of HIV-1. When the relationship between the sampling years of sequences and their nucleotide distances to the common ancestor was extrapolated to the past, analysis of pol sequences provided accurate estimates of the onset years of the epidemics, whereas analysis of V3 sequences by the maximum-likelihood or neighbor-joining methods led to an overestimation of the age of the epidemics. Separate analysis of nonsynonymous and synonymous distances revealed that this overestimation results from nonsynonymous substitutions, whose numbers were not increasing significantly in all three virus populations over the observation period. In contrast, analysis of synonymous env V3 distances provided accurate estimates of the onset years for the outbreaks we studied. Received: 26 October 2001 / Accepted: 8 November 2001     

Evolution of human immunodeficiency virus type 1 in perinatally infected infants with rapid and slow progression to disease.

We addressed the relationship between the origin and evolution of human immunodeficiency virus type 1 (HIV-1) variants and disease outcome in perinatally infected infants by studying the V3 regions of viral variants in samples obtained from five transmitting mothers at delivery and obtained sequentially over the first year of life from their infected infants, two of whom (rapid progressors) rapidly progressed to having AIDS. Phylogenetic analyses disclosed that the V3 sequences from each mother-infant pair clustered together and were clearly distinct from those of the other pairs. Within each pair, the child's sequences formed a monophyletic group, indicating that a single variant initiated the infection in both rapid and slow progressors. Plasma HIV-1 RNA levels increased in all five infants during their first months of life and then declined within the first semester of life only in the three slow progressors. V3 variability increased over time in all infants, but no differences in the pattern of V3 evolution in terms of potential viral phenotype were observed. The numbers of synonymous and nonsynonymous substitutions varied during the first semester of life regardless of viral load, CD4+-cell count, and disease progression. Conversely, during the second semester of life the rate of nonsynonymous substitutions was higher than that of synonymous substitutions in the slow progressors but not in the rapid progressors, thus suggesting a stronger host selective pressure in the former. In view of the proposal that V3 genetic evolution is driven mainly by host immune constraints, these findings suggest that while the immune response to V3 might contribute to regulating viral levels after the first semester of life, it is unlikely to play a determinant role in the initial viral decline soon after birth.     

Evolution of human immunodeficiency virus type 1 env sequence variation in patients with diverse rates of disease progression and T-cell function.

We examined the relationship between env sequence variation and disease progression in 10 human immunodeficiency virus type 1 (HIV-1)-seropositive subjects selected from a longitudinal cohort receiving zidovudine therapy. Five subjects were chosen for stable clinical status and CD4 counts (slow progressors), and five were selected for rapid clinical deterioration and CD4 count decline (rapid progressors). The slow progressors had significantly lower plasma viral RNA loads and greater lymphoproliferative responses to mitogens than the rapid progressors. DNA sequences representing the C1 through C3 regions of env were amplified from two peripheral blood mononuclear cell DNA samples from each subject separated by an average of 2.5 years. Molecular clones of these amplicons were then sequenced, and DNA sequence and deduced amino acid sequence distances were compared. Inter-time point sequence comparison showed a higher rate of sequence evolution for the rapid progressors in three of five matched pairs of rapid progressors and slow progressors and for the slow progressors in the remaining two subject pairs. However, intra-time point sequence comparisons showed that four of five slow progressors developed a more diverse quasispecies over time and one showed no change. In contrast, four of five rapid progressors showed no change in quasispecies diversity over time and one showed a significant decrease in diversity. The overall C1 through C3 region quasispecies diversity in the slow progressors at baseline was lower than that for the rapid progressors, but this difference was not significant at the follow-up time points. These diversity relationships were obscured if sequence analyses were limited to the 300-bp C2 to V3 region. Thus, HIV-1 quasispecies diversity increased over time in subjects with more functional immune systems.     

Unusual distribution of mutations associated with serial bottleneck passages of human immunodeficiency virus type 1

Repeated bottleneck passages result in fitness losses of RNA viruses. In the case of human immunodeficiency virus type 1 (HIV-1), decreases in fitness after a limited number of plaque-to-plaque transfers in MT-4 cells were very drastic. Here we report an analysis of entire genomic nucleotide sequences of four HIV-1 clones derived from the same HIV-1 isolate and their low-fitness progeny following 7 to 15 plaque-to-plaque passages. Clones accumulated 4 to 28 mutations per genome, with dominance of A --> G and G --> A transitions (57% of all mutations) and 49% nonsynonymous replacements. One clone-but not three sibling clones-showed an overabundance of G --> A transitions, evidencing the highly stochastic nature of some types of mutational bias. The distribution of mutations along the genome was very unusual in that mutation frequencies in gag were threefold higher than in env. Particularly striking was the complete absence of replacements in the V3 loop of gp120, confirmed with partial nucleotide sequences of additional HIV-1 clones subjected to repeated bottleneck passages. The analyses revealed several amino acid replacements that have not been previously recorded among natural HIV-1 isolates and illustrate how evolution of an RNA virus genome, with regard to constant and variable regions, can be profoundly modified by alterations in population dynamics.     

Intrahost human immunodeficiency virus type 1 evolution is related to length of the immunocompetent period.

The antigenic diversity threshold theory predicts that antigenic sites of human immunodeficiency virus type 1, such as the V3 region of the external glycoprotein gp120, evolve more rapidly during the symptom-free period in individuals progressing to AIDS than in those who remain asymptomatic for a long time. To test this hypothesis, genomic RNA sequences were obtained from the sera of 44 individuals at seroconversion and 5 years later. The mean number of nonsynonymous nucleotide substitutions in the V3 region of the viruses circulating in 31 nonprogressors (1.1 x 10(-2) +/- 0.1 x 10(-2) per site per year) was higher than the corresponding value for 13 progressors (0.66 x 10(-2) +/- 0.1 x 10(-2) per site per year) (P < 0.01), while no difference between the mean numbers of synonymous substitutions in the two groups was seen (0.37 x 10(-2) +/- 0.1 x 10(-2) and 0.51 x 10(-2) +/- 0.2 x 10(-2) per site per year for nonprogressors and progressors, respectively; P > 0.1). The mean ratios of synonymous nucleotide p distance to nonsynonymous p distance were 0.35 for nonprogressors and 0.62 for progressors. The number of nonsynonymous substitutions was not associated with virus load or virus phenotype, which are established predictors of disease progression, but correlated strongly with the duration of the immunocompetent period (r2 = 0.41; P = 0.001). This indicates that there is no causative relationship between intrahost evolution and CD4+ cell decline. Our data suggest that intrahost evolution in human immunodeficiency virus type 1 infection is driven by selective forces, the strength of which is related to the duration of the immunocompetent period.     

Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing

Accurate identification of the transmitted virus and sequences evolving from it could be instrumental in elucidating the transmission of human immunodeficiency virus type 1 (HIV-1) and in developing vaccines, drugs, or microbicides to prevent infection. Here we describe an experimental approach to analyze HIV-1 env genes as intact genetic units amplified from plasma virion RNA by single-genome amplification (SGA), followed by direct sequencing of uncloned DNA amplicons. We show that this strategy precludes in vitro artifacts caused by Taq-induced nucleotide substitutions and template switching, provides an accurate representation of the env quasispecies in vivo, and has an overall error rate (including nucleotide misincorporation, insertion, and deletion) of less than 8 x 10(-5). Applying this method to the analysis of virus in plasma from 12 Zambian subjects from whom samples were obtained within 3 months of seroconversion, we show that transmitted or early founder viruses can be identified and that molecular pathways and rates of early env diversification can be defined. Specifically, we show that 8 of the 12 subjects were each infected by a single virus, while 4 others acquired more than one virus; that the rate of virus evolution in one subject during an 80-day period spanning seroconversion was 1.7 x 10(-5) substitutions per site per day; and that evidence of strong immunologic selection can be seen in Env and overlapping Rev sequences based on nonrandom accumulation of nonsynonymous mutations. We also compared the results of the SGA approach with those of more-conventional bulk PCR amplification methods performed on the same patient samples and found that the latter is associated with excessive rates of Taq-induced recombination, nucleotide misincorporation, template resampling, and cloning bias. These findings indicate that HIV-1 env genes, other viral genes, and even full-length viral genomes responsible for productive clinical infection can be identified by SGA analysis of plasma virus sampled at intervals typical in large-scale vaccine trials and that pathways of viral diversification and immune escape can be determined accurately.     

Evolution and biological characterization of human immunodeficiency virus type 1 subtype E gp120 V3 sequences following horizontal and vertical virus transmission in a single family

It has been suggested that immune-pressure-mediated positive selection operates to maintain the antigenic polymorphism on the third variable (V3) loop of the gp120 of human immunodeficiency virus type 1 (HIV-1). Here we present evidence, on the basis of sequencing 147 independently cloned env C2/V3 segments from a single family (father, mother, and their child), that the intensity of positive selection is related to the V3 lineage. Phylogenetic analysis and amino acid comparison of env C2/V3 and gag p17/24 regions indicated that a single HIV-1 subtype E source had infected the family. The analyses of unique env C2/V3 clones revealed that two V3 lineage groups had evolved in the parents. Group 1 was maintained with low variation in all three family members regardless of the clinical state or the length of infection, whereas group 2 was only present in symptomatic individuals and was more positively charged and diverse than group 1. Only virus isolates carrying the group 2 V3 sequences infected and induced syncytia in MT2 cells, a transformed CD4(+)-T-cell line. A statistically significant excess of nonsynonymous substitutions versus synonymous substitutions was demonstrated only for the group 2 V3 region. The data suggest that HIV-1 variants, possessing the more homogeneous group 1 V3 element and exhibiting the non-syncytium-inducing phenotype, persist in infected individuals independent of clinical status and appear to be more resistant to positive selection pressure.     

SRY基因在人猿超科和旧大陆猴中具有不同的进化规律

通过ＰＣＲ扩增、测序,得到了白臀叶猴和红面猴的ＳＲＹ基因全序列。结合现有的灵长类其他物种序列进行分析,验证了ＨＭＧ盒的保守性。通过构建系统发育树,比较旧大陆猴和人猿超科两个类群内和类群间ＨＭＧ盒侧翼序列Ｋａ／Ｋｓ的比率。有趣的是,人猿超科两物种比较呈现较高的Ｋａ／Ｋｓ比值,但在旧大陆猴中及旧大陆猴与狨猴间的Ｋａ／Ｋｓ比值显著低于人猿超科的,呈现很不同的格局。同时,对于ＨＭＧ盒序列,Ｋａ／Ｋｓ比值在     

Evolutionary models accounting for layers of selection in protein-coding genes and their impact on the inference of positive selection

The selective forces acting on a protein-coding gene are commonly inferred using evolutionary codon models by contrasting the rate of nonsynonymous substitutions to the rate of synonymous substitutions. These models usually assume that the synonymous substitution rate, Ks, is homogenous across all sites, which is justified if synonymous sites are free from selection. However, a growing body of evidence indicates that the DNA and RNA levels of protein-coding genes are subject to varying degrees of selective constraints due to various biological functions encoded at these levels. In this paper, we develop evolutionary models that account for these layers of selection by allowing for both among-site variability of substitution rates at the DNA/RNA level (which leads to Ks variability among protein-coding sites) and among-site variability of substitution rates at the protein level (Ka variability). These models are constructed so that positive selection is either allowed or not. This enables statistical testing of positive selection when variability at the DNA/RNA substitution rate is accounted for. Using this methodology, we show that variability of the baseline DNA/RNA substitution rate is a widespread phenomenon in coding sequence data of mammalian genomes, most likely reflecting varying degrees of selection at the DNA and RNA levels. Additionally, we use simulations to examine the impact that accounting for the variability of the baseline DNA/RNA substitution rate has on the inference of positive selection. Our results show that ignoring this variability results in a high rate of erroneous positive-selection inference. Our newly developed model, which accounts for this variability, does not suffer from this problem and hence provides a likelihood framework for the inference of positive selection on a background of variability in the baseline DNA/RNA substitution rate.     

γ-MYN: a new algorithm for estimating Ka and Ks with consideration of variable substitution rates

Background  

Over the past two decades, there have been several approximate methods that adopt different mutation models and used for estimating nonsynonymous and synonymous substitution rates (Ka and Ks) based on protein-coding sequences across species or even different evolutionary lineages. Among them, MYN method (a Modified version of Yang-Nielsen method) considers three major dynamic features of evolving DNA sequences–bias in transition/transversion rate, nucleotide frequency, and unequal transitional substitution but leaves out another important feature: unequal substitution rates among different sites or nucleotide positions.     

Coreceptor tropism can be influenced by amino acid substitutions in the gp41 transmembrane subunit of human immunodeficiency virus type 1 envelope protein

Many studies have demonstrated that the third variable region (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope protein (Env) is a major determinant of coreceptor tropism. Other regions in the surface gp120 subunit of Env can modulate coreceptor tropism in a manner that is not fully understood. In this study, we evaluated the effect of env determinants outside of V3 on coreceptor usage through the analysis of (i) patient-derived env clones that differ in coreceptor tropism, (ii) chimeric env sequences, and (iii) site-directed mutants. The introduction of distinct V3 sequences from CXCR4-using clones into an R5-tropic env backbone conferred the inefficient use of CXCR4 in some but not all cases. Conversely, in many cases, X4- and dual-tropic env backbones containing the V3 sequences of R5-tropic clones retained the ability to use CXCR4, suggesting that sequences outside of the V3 regions of these CXCR4-using clones were responsible for CXCR4 use. The determinants of CXCR4 use in a set of dual-tropic env sequences with V3 sequences identical to those of R5-tropic clones mapped to the gp41 transmembrane (TM) subunit. In one case, a single-amino-acid substitution in the fusion peptide of TM was able to confer CXCR4 use; however, TM substitutions associated with CXCR4 use varied among different env sequences. These results demonstrate that sequences in TM can modulate coreceptor specificity and that env sequences other than that of V3 may facilitate efficient CXCR4-mediated entry. We hypothesize that the latter plays an important role in the transition from CCR5 to CXCR4 coreceptor use.     

Natural selection and the organ-specific differentiation of HIV-1 V3 hypervariable region

The existence of organ-specific HIV-1 populations within infected hosts has been studied for many years; nonetheless results reported by different authors are somewhat discrepant. To tackle this problem, we used a population genetics approach to analyze previously published data from the V3 hypervariable region of the envelope env gene. Our results are compatible with a population subdivision by organs in 95% of individuals analyzed at autopsy. In addition, populations infecting the nervous system and testicles clearly appear as differentiated subsets of the so-called macrophage-tropic variants. Liver and kidney may harbor differentiated populations as well. Although it is widely accepted that organ compartmentalization arises as a consequence of different selective pressures imposed by different organs, a definitive demonstration has not yet been provided. Our analysis of the pattern of synonymous and nonsynonymous nucleotide substitutions provides evidence supporting this hypothesis, without discarding the role of other evolutionary processes. In contrast, positive selection does not seem to be the mechanism responsible for the evolution of patient-specific sequences.     

Host-specific driving force in human immunodeficiency virus type 1 evolution in vivo.

To investigate the process of human immunodeficiency virus type 1 (HIV-1) evolution in vivo, a total of 179 HIV-1 V3 sequences derived from cell-free plasma were determined from serial samples in three epidemiologically linked individuals (one infected blood donor and two transfusion recipients) over a maximum period of 8 years. A systematic analysis of pairwise comparisons of intrapatient sequences, both within and between each sample time point, revealed a preponderance and accumulation of nonsynonymous rather than synonymous substitutions in the V3 loop and flanking regions as they diverged over time. This strongly argues for the dominant role that positive selection for amino acid change plays in governing the pattern and process of HIV-1 env V3 evolution in vivo and nullifies hypotheses of purely neutral or mutation-driven evolution or completely chance events. In addition, different rates of evolution of HIV-1 were observed in these three different individuals infected with the same viral strain, suggesting that the degree of positive pressure for HIV-1 amino acid change is host dependent. Finally, the observed similar rate of accumulation in divergence within and between infected individuals suggests that the process of genetic divergence in the HIV epidemic proceeds regardless of host-to-host transmission events, i.e., that transmission does not reset the evolutionary clock.     

Genetically Related Human Immunodeficiency Virus Type 1 in Three Adults of a Family with No Identified Risk Factor for Intrafamilial Transmission

A small number of cases of human immunodeficiency virus (HIV) infection have been reported in individuals with no identified risk factors for transmission. We report on the seroconversion of the 61-year-old mother and the subsequent finding of HIV seropositivity in the 66-year-old father of a 31-year-old AIDS patient. Extensive investigation failed to identify any risk factor for intrafamilial transmission. We conducted a genetic analysis and determined the amino acid signature patterns of the V3, V4, and V5 hypervariable domains and flanking regions in the HIV-1 gp120 env gene of 26 clones derived from proviral DNA in peripheral blood mononuclear cells of the members of the family. env sequences of the viruses isolated from the patients were compared with sequences of HIV-1 subtype B viruses from Europe and local field isolates. Phylogenetic analysis revealed that the sequences of the viruses isolated from the patients were genetically related and formed an intrafamilial cluster of HIV-1 distinct from other subtype B viruses. Interindividual nucleotide variability in the C2-V3 and V4-C4-V5 domains ranged between 1.2 and 5.0% and between 2.2 and 7.5%, respectively, whereas divergence between HIV strains from the patients and control viral strains ranged from 6.6 to 29.3%. The amino acid signature patterns of viral clones from the three patients were closely related. In the C2-V3 region, two minor clones derived from the son’s virus showed less nucleotide divergence (mean, 3.5 and 3.9%) than did the clones derived from the viruses of both parents or the seven other predominant clones derived from the virus from the son (mean, 5.4%). The top of the V3 loop of the last two clones and of all viral clones from the parents exhibited an unusual GPGG sequence. This is the first report of genotypic relatedness of HIV-1 in three adults of the same family in the absence of identified risk factor for transmission between the members of the family. Our findings suggest that atypical transmission of HIV may occur.     

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.     

Somatic variation plays a key role in the evolution of the Vf gene family residing in the Vf locus that confers resistance to apple scab disease

A cluster of four receptor-like genes has been previously identified in the Vf locus of the crabapple Malus floribunda clone 821 that confers resistance to five races of the fungal pathogen Venturia inaequalis, the casual agent of apple scab disease. Pairwise comparisons of the four Vf paralogs in both promoter and coding regions reveal their timeline evolutionary history. The four Vf paralogs have evolved from four ancient Vf members resulting from two sequential duplication events of a single Vf progenitor initially present in the Malus genome. The coding sequences of the four Vf paralogs are characterized with high numbers of unique polymorphic nucleotides, a number of short duplications/deletions, various deletions of complete LRR copy units, and a casual insert of a transposon-like element. Significant high ratios of nonsynonymous to synonymous substitutions, Ka/Ks, are observed in the putative ligand binding residues in the LRR domains. No sequence exchange between the four Vf paralogs is observed. Compared with promoter regions, only nucleotide substitutions are dramatically elevated in the coding regions. The results presented in this study strongly indicate that the Vf locus is under strong and steady horizontal selective pressures imposed by the fungal pathogen V. inaequalis, and divergent selection on somatic variations plays a key role in shaping the resistance specificity.     

Phenotype-associated env gene variation among eight related human immunodeficiency virus type 1 clones: evidence for in vivo recombination and determinants of cytotropism outside the V3 domain.

The nucleotide sequences of the env genes of eight phenotypically heterogeneous human immunodeficiency virus type 1 (HIV-1) clones recovered from a single individual within a 3-week period were compared. In addition, the accessory gene sequences for four of these clones were obtained. Variation among most accessory genes was limited. In contrast, pronounced phenotype-associated sequence variation was observed in the env gene. At least three of these clones most likely resulted from genetic recombination events in vivo, indicating that this phenomenon may account for the emergence of proviruses with novel phenotypic properties. Within the env genes of the eight clones, four domains could be defined, the sequence of each of which clustered in two groups with high internal homology but 11 to 30% cluster variation. The extensive env gene variation among these eight clones could largely be explained by the unique manner in which the alleles of these four domains were combined in each clone. Experiments with chimeric proviruses demonstrated that the HIV-1 env gene determined the capacity to induce syncytia and tropism for T-cell lines. Amino acids previously shown to be involved in gp120-CD4 and gp120-gp41 interaction were completely conserved among these eight clones. The finding of identical V3 sequences in clones differing in tropism for primary monocytes and T-cell lines demonstrated the existence of determinants of tropism outside the env V3 region.