Amino Acid Covariation in a Functionally Important Human Immunodeficiency Virus Type 1 Protein Region Is Associated With Population Subdivision

The frequently reported amino acid covariation of the highly polymorphic human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein V3 region has been assumed to reflect fitness epistasis between residues. However, nonrandom association of amino acids, or linkage disequilibrium, has many possible causes, including population subdivision. If the amino acids at a set of sequence sites differ in frequencies between subpopulations, then analysis of the whole population may reveal linkage disequilibrium even if it does not exist in any subpopulation. HIV-1 has a complex population structure, and the effects of this structure on linkage disequilibrium were investigated by estimating within- and among-subpopulation components of variance in linkage disequilibrium. The amino acid covariation previously reported is explained by differences in amino acid frequencies among virus subpopulations in different patients and by nonsystematic disequilibrium among patients. Disequilibrium within patients appears to be entirely due to differences in amino acid frequencies among sampling time points and among chemokine coreceptor usage phenotypes of virus particles, but not source tissues. Positive selection explains differences in allele frequencies among time points and phenotypes, indicating that these differences are adaptive rather than due to genetic drift. However, the absence of a correlation between linkage disequilibrium and phenotype suggests that fitness epistasis is an unlikely cause of disequilibrium. Indeed, when population structure is removed by analyzing sequences from a single time point and phenotype, no disequilibrium is detectable within patients. These results caution against interpreting amino acid covariation and coevolution as evidence for fitness epistasis.     

小麦异交群体连续歧化选择的RAPD分析

利用太谷核不育基因构建的遗传变异丰富的基础群体DNS2,进行了连续5轮歧化选择。本论文从不同的世代中,选择了10个子群体进行RAPD分析。采用7个引物扩增出116个位点,从基因频率和表型带两个角度的分析都表明,群体具有丰富的遗传变异。整个群体总的多态位点百分率达88.79,总杂合度为0.3143。子群体内（间）遗传距离的结果显示：子群体内的遗传差异小于子群体间的遗传差异;各选择子群体与未选群体间都有较大的遗传距离;随着选择轮次的增加,株高选择子群体间的遗传距离逐渐增大;对同一性状进行选择的子群体间世代内（间）平均遗传距离小于对不同性状进行选择的子群体内（间）的遗传距离。RAPD分子聚类结果显示出对同一性状进行选择的子群体聚在一起,反映了对株高选高的选择效果比较明显。 Abstract:A base population was established through multi-parent random crossing by using Taigu dominant male-sterile wheat,and then five cycles of 2-way selection for four quantitative characters were conducted.The dynamic changes of genetic structures in the open-pollinated wheat population were examined by RAPD technique.Seven primers in RAPD analysis amplified 116 sites.The results of gene frequencies and phenotypic bands showed abundant genetic variations existed in the population.The percentage of polymorphic loci was 8879,and the average heterozygosity was 0.3143 in the whole population.The genetic distance of RAPD showed as follows :① The genetic distance within a subpopulation was lower than that between every two subpopulations.② Each subpopulation had considerable divergence compared with unselected population.③ The genetic distances between the subpopulations which were selected for plant height gradually increased accompanied with the selection.④The genetic distance between subpopulations which were selected for the same character was lower than that were selected for different characters both in the same generation and among different generations.The cluster results of RAPD genetic distance demonstrated that the subpopulations selected for the same character going to one cluster.It also showed that the selective effect of increasing plant height was obvious.     

Founder virus population related to route of virus transmission: a determinant of intrahost human immunodeficiency virus type 1 evolution?

We and others have shown that in individual human immunodeficiency virus type 1 (HIV-1) infection, the adaptive evolution of HIV-1 is influenced by host immune competence. In this study, we tested the hypothesis that in addition to selective forces operating within the host, transmission bottlenecks have an impact on HIV-1 intrahost evolution. Therefore, we studied the intrahost evolution of the V3 region of the external glycoprotein gp120 of HIV-1 during the 3- and 5-year periods following seroconversion after parenteral versus sexual (male-to-male) transmission in 41 participants of the Amsterdam prospective cohorts of homosexual men (n = 31) and intravenous drug users (IVDUs; n = 10) who were AIDS free and had comparable numbers of CD4+ cells. We observed that HIV-1 strains in homosexual men accumulated over 5 years more nonsynonymous substitutions within the V3 loop than HIV-1 strains in IVDUs as a result of lower rates of nonsynonymous evolution in both the initial 3-year period from seroconversion and the following 2-year period as well as a larger proportion of nonsynonymous back substitutions in IVDUs. The mean numbers of synonymous substitutions did not differ between the two risk groups. Since HIV-1 strains in IVDUs could be distinguished from the viruses of homosexual men based on several nucleotide substitutions of which the most conserved is a synonymous substitution at the tip of the V3 loop (GGC pattern), we studied whether the founder virus population itself has an impact on the intrahost evolution of HIV-1. The mean number of nonsynonymous substitutions accumulated over 5 years within the V3 loop was lower in 10 IVDUs infected by the HIV-1 strains with the GGC signature than in 4 IVDUs infected by HIV-1 strains lacking this pattern, while the mean numbers of synonymous substitutions were similar in the two groups.     

The evolutionary history of the CCR5-Delta32 HIV-resistance mutation

The CCR5 chemokine receptor is exploited by HIV-1 to gain entry into CD4+ T cells. A deletion mutation (Delta32) confers resistance against HIV by obliterating the expression of the receptor on the cell surface. Intriguingly, this allele is young in evolutionary time, yet it has reached relatively high frequencies in Europe. These properties indicate that the mutation has been under intense positive selection. HIV-1 has not exerted selection for long enough on the human population to drive the CCR5-Delta32 allele to current frequencies, fueling debate regarding the selective pressure responsible for rise of the allele. The allele exists at appreciable frequencies only in Europe, and within Europe, the frequency is higher in the north. Here we review the population genetics of the CCR5 locus, the debate over the historical selective pressure acting on CCR5-Delta32, the inferences that can potentially be drawn from the geographic distribution of CCR5-Delta32 and the role that other genetic polymorphisms play in conferring resistance against HIV. We also discuss parallel evolution that has occurred at the CCR5 locus of other primate species. Finally, we highlight the promise that therapies based on interfering with the CCR5 receptor could have in the treatment of HIV.     

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.     

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.     

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

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

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.     

Temporal variations in the hepatitis C virus intrahost population during chronic infection

The intrahost evolution of hepatitis C virus (HCV) holds keys to understanding mechanisms responsible for the establishment of chronic infections and to development of a vaccine and therapeutics. In this study, intrahost variants of two variable HCV genomic regions, HVR1 and NS5A, were sequenced from four treatment-naïve chronically infected patients who were followed up from the acute stage of infection for 9 to 18 years. Median-joining network analysis indicated that the majority of the HCV intrahost variants were observed only at certain time points, but some variants were detectable at more than one time point. In all patients, these variants were found organized into communities or subpopulations. We hypothesize that HCV intrahost evolution is defined by two processes: incremental changes within communities through random mutation and alternations between coexisting communities. The HCV population was observed to incrementally evolve within a single community during approximately the first 3 years of infection, followed by dispersion into several subpopulations. Two patients demonstrated this pattern of dispersion for the rest of the observation period, while HCV variants in the other two patients converged into another single subpopulation after ∼9 to 12 years of dispersion. The final subpopulation in these two patients was under purifying selection. Intrahost HCV evolution in all four patients was characterized by a consistent increase in negative selection over time, suggesting the increasing HCV adaptation to the host late in infection. The data suggest specific staging of HCV intrahost evolution.     

Selection promotes organ compartmentalization in HIV-1: evidence from gag and pol genes

The existence of organ-specific human immunodeficiency virus type 1 (HIV-1) populations within infected hosts has been long lasting studied. Previous work established that population subdivision by organs occurs at the envelope env gene, but less is known about other genomic regions. Here, we used a population genetics approach to detect organ compartmentalization in proviral sequences of HIV-1 gag and pol genes. Significant population structure was found in pol (100% of cases) and gag (33%) pair-wise organ comparisons. The degree of compartmentalization positively correlated with the ratio of nonsynonymous to synonymous substitutions, and codons showing organ compartmentalization were more likely to be under significantly positive selection. This suggests that HIV-1 populations dynamically adapt to locally variable intra-host environments. In the case of pol gene, differential penetration of antiretroviral drugs might account for the observed pattern, whereas for gag gene, local selective pressures remain unexplored.     

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.     

Evolution of Envelope Sequences from the Genital Tract and Peripheral Blood of Women Infected with Clade A Human Immunodeficiency Virus Type 1

The development of viral diversity during the course of human immunodeficiency virus type 1 (HIV-1) infection may significantly influence viral pathogenesis. The paradigm for HIV-1 evolution is based primarily on studies of male cohorts in which individuals were presumably infected with a single virus variant of subtype B HIV-1. In this study, we evaluated virus evolution based on sequence information of the V1, V2, and V3 portions of HIV-1 clade A envelope genes obtained from peripheral blood and cervical secretions of three women with genetically heterogeneous viral populations near seroconversion. At the first sample following seroconversion, the number of nonsynonymous substitutions per potential nonsynonymous site (dn) significantly exceeded substitutions at potential synonymous sites (ds) in plasma viral sequences from all individuals. Generally, values of dn remained higher than values of ds as sequences from blood or mucosa evolved. Mutations affected each of the three variable regions of the envelope gene differently; insertions and deletions dominated changes in V1, substitutions involving charged amino acids occurred in V2, and sequential replacement of amino acids over time at a small subset of positions distinguished V3. The relationship among envelope nucleotide sequences obtained from peripheral blood mononuclear cells, plasma, and cervical secretions was evaluated for each individual by both phylogenetic and phenetic analyses. In all subjects, sequences from within each tissue compartment were more closely related to each other than to sequences from other tissues (phylogenetic tissue compartmentalization). At time points after seroconversion in two individuals, there was also greater genetic identity among sequences from the same tissue compartment than among sequences from different tissue compartments (phenetic tissue compartmentalization). Over time, temporal phylogenetic and phenetic structure was detectable in mucosal and plasma viral samples from all three women, suggesting a continual process of migration of one or a few infected cells into each compartment followed by localized expansion and evolution of that population.     

Population-Enrichment Adaptive Design Strategy for an Event-Driven Vaccine Efficacy Trial

A population-enrichment adaptive design allows a prospective use for study population selection. It has the flexibility allowing pre-specified modifications to an ongoing trial to mitigate the potential risk associated with the assumptions made at design stage. In this way, the trial can potentially encompass a broader target patient population, and move forward only with the subpopulations that appear to be benefiting from the treatment. Our work is motivated by a Phase III event-driven vaccine efficacy trial. Two target patient subpopulations were enrolled with the assumption that vaccine efficacy can be demonstrated based on the combined population. It is recognized due to the nature of patients’ underlying conditions, one subpopulation might respond to the treatment better than the other. To maximize the probability of demonstrating vaccine efficacy in at least one patient population while taking advantage of combining two subpopulations in one single trial, an adaptive design strategy with potential population enrichment is developed. Specifically, if the observed vaccine efficacy at interim for one subpopulation is not promising to warrant carrying forward, the population may be enriched with the other subpopulation with better performance. Simulations were conducted to evaluate the operational characteristics from a selection of interim analysis plans. This population-enrichment design provides a more efficient way as compared to the conventional approaches when targeting multiple subpopulations. If executed and planned with caution, this strategy can provide a greater chance of success of the trial and help maintain scientific and regulatory rigors.     

The geographic spread of the CCR5 Delta32 HIV-resistance allele

The Delta32 mutation at the CCR5 locus is a well-studied example of natural selection acting in humans. The mutation is found principally in Europe and western Asia, with higher frequencies generally in the north. Homozygous carriers of the Delta32 mutation are resistant to HIV-1 infection because the mutation prevents functional expression of the CCR5 chemokine receptor normally used by HIV-1 to enter CD4+ T cells. HIV has emerged only recently, but population genetic data strongly suggest Delta32 has been under intense selection for much of its evolutionary history. To understand how selection and dispersal have interacted during the history of the Delta32 allele, we implemented a spatially explicit model of the spread of Delta32. The model includes the effects of sampling, which we show can give rise to local peaks in observed allele frequencies. In addition, we show that with modest gradients in selection intensity, the origin of the Delta32 allele may be relatively far from the current areas of highest allele frequency. The geographic distribution of the Delta32 allele is consistent with previous reports of a strong selective advantage (>10%) for Delta32 carriers and of dispersal over relatively long distances (>100 km/generation). When selection is assumed to be uniform across Europe and western Asia, we find support for a northern European origin and long-range dispersal consistent with the Viking-mediated dispersal of Delta32 proposed by G. Lucotte and G. Mercier. However, when we allow for gradients in selection intensity, we estimate the origin to be outside of northern Europe and selection intensities to be strongest in the northwest. Our results describe the evolutionary history of the Delta32 allele and establish a general methodology for studying the geographic distribution of selected alleles.     

Evolution of intrahost HIV-1 genetic diversity during chronic infection

HIV-1 is one of the fastest evolving entities known. Given that census population sizes of HIV-1 within individuals are much greater than the inverse mutation rate, every possible single point mutation in the viral genome occurs each generation. This enormous capability to generate genetic variation allows for escape from immune surveillance and antiviral therapy. However, compared to this potential, populations of HIV-1 within individuals exhibit little genetic variation. This discrepancy between the known mutation rate of HIV-1 and the average level of genetic variation in the env gene observed in vivo is reflected in comparisons of the actual numbers of productively infected cells, estimated as 10(7), and the effective population size, estimated as 10(3). Using approximate Bayesian computation, we evaluated several hypotheses based on a variety of selective and demographic processes to explain the low effective population size of HIV-1. Of the models we examined, the metapopulation model, in which HIV-1 evolves within an individual as a large collection of small subpopulations subject to frequent migration, extinction, and recolonization, was most consistent with the observed levels of genetic variation and the average frequencies of those variants. The metapopulation model links previous studies of viral dynamics and population genetics.     

Linkage Disequilibrium in Subdivided Populations

The linkage disequilibrium in a subdivided populaton is shown to be equal to the sum of the average linkage disequilibrium for all subpopulations and the covariance between gene frequencies of the loci concerned. Thus, in a subdivided population the linkage disequilibrium may not be 0 even if the linkage disequilibrium in each subpopulation is 0. If a population is divided into two subpopulations between which migration occurs, the asymptotic rate of approach to linkage equilibrium is equal to either r or 2(m(1) + m(2)) - (m(1) + m(2))(2), whichever is smaller, where r is the recombination value and m(1) and m(2) are the proportions of immigrants in subpopulations 1 and 2, respectively. Thus, if migration rate is high compared with recombination value, the change of linkage disequilibrium in subdivided populations is similar to that of a single random mating population. On the other hand, if migration rate is low, the approach to lnkage equilibrium may be retarded in subdivided populations. If isolated populations begin to exchange genes by migration, linkage disequilibrium may increase temporarily even for neutral loci. If overdominant selection operates and the equilibrium gene frequencies are different in the two subpopulations, a permanent linkage disequilibrium may be produced without epistasis in each subpopulation.     

Use of a genetic marker to examine genetic interaction among subpopulations of pink salmon (Oncorhynchus gorbuscha)

In 1979 and 1981, a genetic marker was bred into one of the five identifiable subpopulations of pink salmon [Oncorhynchus gorbuscha (Walbaum)] in the Auke Lake drainage in Southeast Alaska. As a result of the marking effort, the frequencies of two malate dehydrogenase (MDH-B1, 2^*) alleles were changed in the marked subpopulation, but not in other subpopulations that spawn at different times or places. Between 1983 and 1989, the marker allele frequencies were monitored in many of these subpopulations and in early- and late-run pink salmon spawning in nearby Waydelich Creek, located approximately 1km away. Changes in allele frequencies at MDH-B1, 2^*, used to obtain direct estimates of average migration rates (m) from the marked to the unmarked subpopulations, revealed little or no introgression into early subpopulations or into nearby Waydelich Creek. Moreover, spatially distinct late-run Auke Creek subpopulations were not immediately overrun by the more abundant marked subpopulation. These observations suggest that genetic isolation exists between temporally distinct spawning runs and that small temporal and spatial (or ecological) differences contribute to population structure. These observations should be considered in taking actions that affect conservation and harvest management or extensive culture of salmonids.     

QUASI analysis of the HIV-1 envelope sequences in the Los Alamos National Laboratory HIV sequence database: pattern and distribution of positive selection sites and their frequencies over years.

The envelope (env) protein of human immunodeficiency virus type 1 (HIV-1) plays a crucial role in virus entry and is a central target for HIV vaccine design. Using the QUASI program, we analyzed the conserved regions of all currently available env sequences in the Los Alamos National Laboratory HIV Sequence Database and identified positive selection (PS) sites that are likely to be restricted by host immune responses. We found that PS sites are dispersed across conserved regions of env sequence, and that the C3, C4, and C5 regions were the most targeted. Several regions were identified as being PS free and were mainly distributed in the C1 and C2 regions. When comparing individual QUASI PS site frequencies across clades and geographical regions with the overall frequency of the entire env database, the env sequences from North America showed significantly lower PS site frequency, while those from Asia were significantly higher using Student's t test. The QUASI PS site frequency of env proteins from viruses isolated from different years showed that the PS site frequencies of the env population increased over time. Our study provides an overview of PS sites across the conserved regions of HIV-1 env sequences.