Recent studies have demonstrated the emergence of human immunodeficiency virus type 1 (HIV-1) subtypes with various levels of fitness. Using heterogeneous maximum-likelihood models of adaptive evolution implemented in the PAML software package, with env sequences representing each HIV-1 group M subtype, we examined the various intersubtype selective pressures operating across the env gene. We found heterogeneity of evolutionary mechanisms between the different subtypes with a category of amino acid sites observed that had undergone positive selection for subtypes C, F1, and G, while these sites had undergone purifying selection in all other subtypes. Also, amino acid sites within subtypes A and K that had undergone purifying selection were observed, while these sites had undergone positive selection in all other subtypes. The presence of such sites indicates heterogeneity of selective pressures within HIV-1 group M subtype evolution that may account for the various levels of fitness of the subtypes.It has been hypothesized that human immunodeficiency virus type 1 (HIV-1) may have entered humans in three independent transmissions of simian immunodeficiency virus from infected chimpanzees from which the three HIV-1 M, N, and O lineages arose (9). Within group M, nine phylogenetically distinct subtypes have been proposed (subtypes A to D, F to H, J, and K), with subsubtypes being proposed for subtypes A and F (18). Subtype distribution varies worldwide, with subtype B predominating in North America and Europe (15) and subtype C accounting for more than 55% of worldwide infections (7) due mainly to its prevalence in Southern and Eastern Africa (1,2,14,17,26) and India (21) and its increasing prevalence in Brazil (22) and China (19). Biological differences including low CXCR4 coreceptor usage in subtype C (15), decreased protease susceptibility in subtype G (5), and varying subtype reactivity to monoclonal antibodies (13, 25) have been observed among the subtypes. The production of broadly neutralizing or subtype-specific vaccines requires an in-depth understanding of the inter-and intrasubtype evolution.The study of the selective pressures governing the evolution of protein-coding DNA sequences has traditionally been carried out by comparing dN (nonsynonymous substitutions per nonsynonymous site) to dS (synonymous substitutions per synonymous site), resulting in a dN-to-dS ratio () (see reference 31 for a review). An of Ͼ1 is indicative of positive selection, an of 1 indicates neutral evolution, and an of Ͻ1 indicates purifying (negative) selection. However, if there is strong purifying selection operating on the majority of amino acid positions, averaging over an entire sequence could misleadingly indicate purifying selection for the entire molecule even in the presence of a small number of sites undergoing positive selection. Recently, more biologically realistic methods have been developed to allow for identification of heterogeneous selection pressure across amino acid sites and also heterogeneity across both si...