SIV Genome-Wide Pyrosequencing Provides a Comprehensive and Unbiased View of Variation within and outside CD8 T Lymphocyte Epitopes

Deep sequencing technology is revolutionizing our understanding of HIV/SIV evolution. It is known that acute SIV sequence variation within CD8 T lymphocyte (CD8-TL) epitopes is similar among MHC-identical animals, but we do not know whether this persists into the chronic phase. We now determine whether chronic viral variation in MHC-identical animals infected with clonal SIV is similar throughout the entire coding sequence when using a sensitive deep sequencing approach. We pyrosequenced the entire coding sequence of the SIV genome isolated from a unique cohort of four SIVmac239-infected, MHC-identical Mauritian cynomolgus macaques (MCM) 48 weeks after infection; one MCM in the cohort became an elite controller. Among the three non-controllers, we found that genome-wide sequences were similar between animals and we detected increased sequence complexity within 64% of CD8-TL epitopes when compared to Sanger sequencing methods. When we compared sequences between the MHC-matched controller and the three non-controllers, we found the viral population in the controller was less diverse and accumulated different variants than the viral populations in the non-controllers. Importantly, we found that initial PCR amplification of viral cDNA did not significantly affect the sequences detected, suggesting that data obtained by pyrosequencing PCR-amplified viral cDNA accurately represents the diversity of sequences replicating within an animal. This demonstrates that chronic sequence diversity across the entire SIV coding sequence is similar among MHC-identical animals with comparable viral loads when infected with the same clonal virus stock. Additionally, our approach to genome-wide SIV sequencing accurately reflects the diversity of sequences present in the replicating viral population. In sum, our study suggests that genome-wide pyrosequencing of immunodeficiency viruses captures a thorough and unbiased picture of sequence diversity, and may be a useful approach to employ when evaluating which sequences to include as part of a vaccine immunogen.     

Route of simian immunodeficiency virus inoculation determines the complexity but not the identity of viral variant populations that infect rhesus macaques

A better understanding of the host and viral factors associated with human immunodeficiency virus (HIV) transmission is essential to developing effective strategies to curb the global HIV epidemic. Here we used the rhesus macaque-simian immunodeficiency virus (SIV) animal model of HIV infection to study the range of viral genotypes that are transmitted by different routes of inoculation and by different types of viral inocula. Analysis of transmitted variants was undertaken in outbred rhesus macaques inoculated intravenously (IV) or intravaginally (IVAG) with a genetically heterogeneous SIVmac251 stock derived from a well-characterized rhesus macaque viral isolate. In addition, we performed serial IV and IVAG passage experiments using plasma from SIV-infected macaques as the inoculum. We analyzed the V1-V2 region of the SIV envelope gene from virion-associated RNA in plasma from infected animals by the heteroduplex mobility assay (HMA) and by DNA sequence analysis. We found that a more diverse population of SIV genetic variants was present in the earliest virus-positive plasma samples from all five IV SIVmac251-inoculated monkeys and from two of five IVAG SIVmac251-inoculated monkeys. In contrast, we found a relatively homogeneous population of SIV envelope variants in three of five monkeys inoculated IVAG with SIVmac251 stock and in two monkeys infected after IVAG inoculation with plasma from an SIV-infected animal. In some IVAG-inoculated animals, the transmitted SIV variant was the most common variant in the inoculum. However, a specific viral variant in the SIVmac251 stock was not consistently transmitted by IVAG inoculation. Thus, it is likely that host factors or stochastic processes determine the specific viral variants that infect an animal after IVAG SIV exposure. In addition, our results clearly demonstrate that the route of inoculation is associated with the extent and breadth of the genetic complexity of the viral variant population in the earliest stages of systemic infection.     

A dominant role for CD8+-T-lymphocyte selection in simian immunodeficiency virus sequence variation

CD8(+) T lymphocytes (CD8-TL) select viral escape variants in both human immunodeficiency virus and simian immunodeficiency virus (SIV) infections. The frequency of CD8-TL viral escape as well as the contribution of escape to overall virus diversification has not been assessed. We quantified CD8-TL selection in SIV infections by sequencing viral genomes from 35 SIVmac239-infected animals at the time of euthanasia. Here we show that positive selection for sequences encoding 46 known CD8-TL epitopes is comparable to the positive selection observed for the variable loops of env. We also found that >60% of viral variation outside of the viral envelope occurs within recognized CD8-TL epitopes. Therefore, we conclude that CD8-TL selection is the dominant cause of SIV diversification outside of the envelope.     

Population Genetic Analysis of the Protease Locus of Human Immunodeficiency Virus Type 1 Quasispecies Undergoing Drug Selection, Using a Denaturing Gradient-Heteroduplex Tracking Assay

Monitoring the evolution of human immunodeficiency virus type 1 (HIV-1) drug resistance requires measuring the frequency of closely related genetic variants making up the complex viral quasispecies found in vivo. In order to resolve both major and minor (>/=2%) protease gene variants differing by one or more nucleotide substitutions, we analyzed PCR products derived from plasma viral quasispecies by using a combination of denaturing gradient gel electrophoresis and DNA heteroduplex tracking assays. Correct population sampling of the high level of genetic diversity present within viral quasispecies could be documented by parallel analysis of duplicate, independently generated PCR products. The composition of genetically complex protease gene quasispecies remained constant over short periods of time in the absence of treatment and while plasma viremia fell >100-fold following the initiation of protease inhibitor ritonavir monotherapy. Within a month of initiating therapy, a strong reduction in the genetic diversity of plasma viral populations at the selected protease locus was associated with rising plasma viremia and the emergence of drug resistance. The high levels of protease genetic diversity seen before treatment reemerged only months later. In one patient, reduction in genetic diversity at the protease gene was observed concomitantly with an increase in diversity at the envelope gene (E. L. Delwart, P. Heng, A. Neumann, and M. Markowitz, J. Virol. 72:2416-2421, 1998), indicating that opposite population genetic changes can take place in different HIV-1 loci. The rapid emergence of drug-resistant HIV-1 was therefore associated with a strong, although only transient, reduction in genetic diversity at the selected locus. The denaturing gradient-heteroduplex tracking assay is a simple method for the separation and quantitation of very closely related, low-frequency, genetic variants within complex viral populations.     

Probe-free allele-specific copy number detection and analysis of tumors

Cancer development and progression frequently involve nucleotide mutations as well as amplifications and deletions of genomic segments. Quantification of allele-specific copy number is an important step in characterizing tumor genomes for precision medicine. Despite advances in approaches to high-throughput genomic DNA analysis, inexpensive and simple methods for analyzing complex nucleotide and copy number variants are still needed. Real-time polymerase chain reaction (PCR) methods for discovering and genotyping single nucleotide polymorphisms are becoming increasingly important in genetic analysis. In this study, we describe a simple, single-tube, probe-free method that combines SYBR Green I-based quantitative real-time PCR and quantitative melting curve analysis both to detect specific nucleotide variants and to quantify allele-specific copy number variants of tumors. The approach is based on the quantification of the targets of interest and the relative abundance of two alleles in a single tube. The specificity, sensitivity, and utility of the assay were demonstrated in detecting allele-specific copy number changes critical for carcinogenesis and therapeutic intervention. Our approach would be useful for allele-specific copy number analysis or precise genotyping.     

Inhibition of Henipavirus fusion and infection by heptad-derived peptides of the Nipah virus fusion glycoprotein

Background

Oral infection of infant macaques with simian immunodeficiency virus (SIV) is a useful animal model to test interventions to reduce postnatal HIV transmission via breast-feeding. We previously demonstrated that immunization of infant rhesus macaques with either modified vaccinia virus Ankara (MVA) expressing SIV Gag, Pol and Env, or live-attenuated SIVmac1A11 resulted in lower viremia and longer survival compared to unimmunized controls after oral challenge with virulent SIVmac251 (Van Rompay et al., J. Virology 77:179–190, 2003). Here we evaluate the impact of these vaccines on oral transmission and evolution of SIV envelope variants.

Results

Limiting dilution analysis of SIV RNA followed by heteroduplex mobility assays of the V1–V2 envelope (env) region revealed two major env variants in the uncloned SIVmac251 inoculum. Plasma sampled from all infants 1 week after challenge contained heterogeneous SIV env populations including one or both of the most common env variants in the virus inoculum; no consistent differences in patterns of env variants were found between vaccinated and unvaccinated infants. However, SIV env variant populations diverged in most vaccinated monkeys 3 to 5 months after challenge, in association with the development of neutralizing antibodies.

Conclusions

These patterns of viral envelope diversity, immune responses and disease course in SIV-infected infant macaques are similar to observations in HIV-infected children, and underscore the relevance of this pediatric animal model. The results also support the concept that neonatal immunization with HIV vaccines might modulate disease progression in infants infected with HIV by breast-feeding.     

Replication and persistence of simian immunodeficiency virus variants after passage in macaque lymphocytes and established human cell lines.

In lentivirus infections, there are typically few cells in the host that harbor the provirus. For this reason, molecular clones of human and simian immunodeficiency viruses (HIV and SIV) are generally derived after passage and amplification of the virus in cell culture. To determine whether SIV variants that persist in culture are similar to the variants that predominate in the host, we examined the proviral sequence of the SIV envelope (env) gene before and after cocultivation of lymphocytes from a macaque with AIDS with naive macaque lymphocytes or human cell lines. Many of the predominant variants in the monkey replicated and persisted in macaque lymphocytes and CEMx174 cells in culture, but a more limited population of variants replicated in C8166 cells. Passage of virus, harvested after 4 weeks of coculture, onto naive cells further demonstrated that the majority of proviruses detected by polymerase chain reaction were also viral variants that were expressed and packaged into infectious virions.     

Unique pattern of convergent envelope evolution in simian immunodeficiency virus-infected rapid progressor macaques: association with CD4-independent usage of CCR5

The rate of disease development in simian immunodeficiency virus (SIV) infection of macaques varies considerably among individual macaques. While the majority of macaques inoculated with pathogenic SIV develop AIDS within a period of 1 to 2 years, a minority exhibit a rapid disease course characterized by absence or transience of humoral and cellular immune responses and high levels of virus replication with widespread dissemination of SIV in macrophages and multinucleated giant cells. The goal of this study was to examine viral evolution in three SIVsmE543-3-inoculated rapid progressors to determine the contribution of viral evolution to the development of rapid disease and the effect of the absence of immune pressure upon viral evolution. PCR was used to amplify and clone the entire SIV genome from tissues collected at necropsy, and the course of viral evolution was assessed by env sequences cloned from sequential plasma samples of one rapid progressor (RP) macaque. The majority of sequence changes in RP macaques occurred in the envelope gene. Substitutions were observed in all three animals at specific conserved residues in envelope, including loss of a glycosylation site in V1/V2, a D-to-N/V substitution in a highly conserved GDPE motif, and a P-to-V/H/T substitution in the V3 loop analog. A cell-cell fusion assay revealed that representative env clones utilized CCR5 as a coreceptor, independent of CD4. The selection of specific substitutions in envelope in RP macaques suggests novel selection pressures on virus in such animals and suggests that viral variants that evolve in these animals may play a role in disease progression.     

Impaired RNA incorporation and dimerization in live attenuated leader-variants of SIVmac239

Background

The ability of emerging pathogens to infect new species is likely related to the diversity of pathogen variants present in existing reservoirs and their degree of genomic plasticity, which determines their ability to adapt to new environments. Certain simian immunodeficiency viruses (SIVcpz, SIVsm) have demonstrated tremendous success in infecting new species, including humans, resulting in the HIV-1 and HIV-2 epidemics. Although SIV diversification has been studied on a population level, the essential substrates for cross-species transmission, namely SIV sequence diversity and the types and extent of viral diversification present in individual reservoir animals have not been elucidated. To characterize this intra-host SIV diversity, we performed sequence analyses of clonal viral envelope (env) V1V2 and gag p27 variants present in individual SIVsm-infected sooty mangabeys over time.

Results

SIVsm demonstrated extensive intra-animal V1V2 length variation and amino acid diversity (le 38%), and continual variation in V1V2 N-linked glycosylation consensus sequence frequency and location. Positive selection was the predominant evolutionary force. Temporal sequence shifts suggested continual selection, likely due to evolving antibody responses. In contrast, gag p27 was predominantly under purifying selection. SIVsm V1V2 sequence diversification is at least as great as that in HIV-1 infected humans, indicating that extensive viral diversification in and of itself does not inevitably lead to AIDS.

Conclusion

Positive diversifying selection in this natural reservoir host is the engine that has driven the evolution of the uniquely adaptable SIV/HIV envelope protein. These studies emphasize the importance of retroviral diversification within individual host reservoir animals as a critical substrate in facilitating cross-species transmission.     

Whole-Genome Characterization of Human and Simian Immunodeficiency Virus Intrahost Diversity by Ultradeep Pyrosequencing

Rapid evolution and high intrahost sequence diversity are hallmarks of human and simian immunodeficiency virus (HIV/SIV) infection. Minor viral variants have important implications for drug resistance, receptor tropism, and immune evasion. Here, we used ultradeep pyrosequencing to sequence complete HIV/SIV genomes, detecting variants present at a frequency as low as 1%. This approach provides a more complete characterization of the viral population than is possible with conventional methods, revealing low-level drug resistance and detecting previously hidden changes in the viral population. While this work applies pyrosequencing to immunodeficiency viruses, this approach could be applied to virtually any viral pathogen.The viral population within each human immunodeficiency virus (HIV)-infected individual is highly diverse and constantly evolving (2, 3). However, our understanding of the viral population is based largely on the consensus sequence of the dominant circulating virus because the full diversity of the viral population is extremely difficult to characterize. One recent study showed that despite viral fitness recovery in vitro, recovery was not correlated with changes observed in the consensus sequence of HIV. Instead, increased fitness correlated with general viral heterogeneity (5). This finding suggests that by limiting our studies to consensus sequences, we are missing many aspects of viral evolution that influence fitness, drug resistance, and immune evasion, among other characteristics. Studies that examined minor viral variants have provided new insights into HIV transmission and pathogenesis, with direct implications for HIV treatment (7, 13). Unfortunately, traditional techniques to identify rare variants, such as molecular cloning, single-genome amplification, or quantitative real-time (qRT)-PCR, are either labor intensive or restricted to the detection of single variants, limiting their widespread use (8, 11, 12, 14).New second-generation technologies have radically altered DNA sequencing. Recent work by our group and others has employed pyrosequencing for targeted ultradeep sequencing of short regions of the viral genome, including CD8⁺ T-lymphocyte epitopes and regions of known drug resistance mutations, demonstrating a practical method to identify extremely low-frequency viral variants (4, 15). While sequencing short regions is appropriate in certain circumstances, the region of interest must be identified in advance, and the effect of mutations in that region on the remaining genome is ignored. Studying the heterogeneity of HIV across the entire genome may provide insights into interactions between minor variants, improve our understanding of HIV evolution, and ultimately provide insights into disease pathogenesis.In this study, we combined pyrosequencing with a transposon-based fragmentation method to allow powerful ultradeep sequencing of the full-length HIV and simian immunodeficiency virus (SIV) genomes, demonstrating a new and highly practical approach to study the complexity of the viral population within a host and identify minor variants on a genome-wide scale. While this study applied pyrosequencing to immunodeficiency viruses, this approach could be applied to any viral pathogen.     

Genotypic selection of simian immunodeficiency virus in macaque infants infected transplacentally.

To understand viral and host factors that contribute to transplacental transmission of human immunodeficiency virus, we developed an animal model using pregnant female macaques infected with simian immunodeficiency virus (SIV). Pregnant females were inoculated intravenously during midgestation with either a well-characterized primary isolate of SIV (SIV/DeltaB670) or a combination of SIV/DeltaB670 and the macrophage-tropic molecular clone SIV/17E-Fr. The viral genetic diversity in five infected female macaques and their in utero-infected infants was analyzed. All of the mothers harbored a genetically diverse virus population at parturition, whereas a single genotype from the maternal quasispecies was identified in the infants at birth. Only one of two variants was found in the infants: SIV/17E-Fr (two cases) or a genotype contained within the SIV/DeltaB670 quasispecies (three cases). The macrophage-tropic properties of both transmitted genotypes were suggested by productive replication in primary rhesus macrophage cultures in vitro and the clonal presence in central nervous system tissue of infected monkeys with encephalitis. These observations provide compelling evidence for both genotypic and phenotypic selection in transplacental transmission of SIV and suggest a critical role for macrophages in fetal infection in utero.     

A chimeric protein of simian immunodeficiency virus envelope glycoprotein gp140 and Escherichia coli aspartate transcarbamoylase

The envelope glycoproteins of the human immunodeficiency virus and the related simian immunodeficiency virus (SIV) mediate viral entry into host cells by fusing viral and target cell membranes. We have reported expression, purification, and characterization of gp140 (also called gp160e), the soluble, trimeric ectodomain of the SIV envelope glycoprotein, gp160 (B. Chen et al., J. Biol. Chem. 275:34946-34953, 2000). We have now expressed and purified chimeric proteins of SIV gp140 and its variants with the catalytic subunit (C) of Escherichia coli aspartate transcarbamoylase (ATCase). The fusion proteins (SIV gp140-ATC) bind viral receptor CD4 and a number of monoclonal antibodies specific for SIV gp140. The chimeric molecule also has ATCase activity, which requires trimerization of the ATCase C chains. Thus, the fusion protein is trimeric. When ATCase regulatory subunit dimers (R(2)) are added, the fusion protein assembles into dimers of trimers as expected from the structure of C(6)R(6) ATCase. Negative-stain electron microscopy reveals spikey features of both SIV gp140 and SIV gp140-ATC. The production of the fusion proteins may enhance the possibilities for structure determination of the envelope glycoprotein either by electron cryomicroscopy or X-ray crystallography.     

Genome-wide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ～75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.     

Genotypic analysis of infant macaques infected transplacentally and orally

Abstract: The SIV-infected macaque provides an excellent model to study factors involved in maternal-fetal transmission of HIV. In our prenatal transmission studies, female macaques were inoculated intravenously during midgestation with either SIV/DeltaB670 or a combination of SIV/DeltaB670 and the macrophage-tropic molecular clone SIV/17E-Fr. The females harbored a genetically diverse virus population at parturition, whereas a single genotype from the maternal quasispecies was identified in the infants. One of two variants was transplacentally transmitted to the infants, SIV/17E-Fr or B670-Cl 12, a genotype contained within the SIV/DeltaB670 inoculum. Both of these variants have been identified in the central nervous system of macaques that have developed encephalitis and they replicate in vitro on primary rhesus macrophages. These results suggest a critical role for macrophages in fetal infection in utero. In our perinatal transmission studies we have evaluated the viral genotypes found in two newborn macaques infected orally with SIV/DeltaB670 and in one infant infected via amniotic inoculation in late gestation. More than one viral genotype was identified in each infant, moreover, each infant harbored different genotypes. These results suggest different mechanisms are responsible for viral infection via these routes.     

Teleost Fish Genomes Contain a Diverse Array of L1 RetrotransposonLineages That Exhibit a Low Copy Number and High Rate of Turnover

Retrotransposable elements exhibit a wide range of variation in population dynamics, abundance, and lineage diversity among host genomes across taxa. This range of diversity is illustrated by a single well-defined constituent monophyletic clade of L1 non-LTR retrotransposons that is shared between mammalian and teleost fish genomes. Despite the clear phylogenetic relationships that exist between mammalian and teleost L1 sequences, these elements exhibit markedly different dynamics within their respective taxa. While mammalian genomes typically contain a single, abundant lineage of L1 elements that traces millions of years of evolution, the zebraflsh genome was recently shown to exhibit a high diversity of ancient lineages coexisting at a very low copy number and apparently exhibiting a high rate of turnover. In the present study, a combination of degenerate PCR, lineage-specific PCR, and genomic Southern blot analysis is utilized to demonstrate high L1 lineage diversity, low copy number, and a high proportion of polymorphic inserts in the genomes of the killifish species, Fundulus heteroclitus. Additional species surveyed by degenerate PCR include Cyprinodon variegatus, Rivulus marmoratus, and Menidia beryllina. These results further support the generality of the differences that exist in host–element dynamics between teleost fish and mammalian genomes with regard to L1 retrotransposons.Reviewing Editor: Dr. Axel Meyer     

Phylogenomic analysis of the L1 retrotransposons in Deuterostomia

L1 retrotransposons constitute the largest single component of mammalian genomes. In contrast to the single remaining lineage of L1 retrotransposons in mammalian genomes, some teleost fishes contain a highly diverse L1 retrotransposon repertoire. Major evolutionary changes in L1 retrotransposon repertoires have therefore taken place in the land vertebrates (Tetrapoda). The lack of sequence data for L1 retrotransposons in the basal living Tetrapoda lineages prompted an investigation of their distribution and evolution in the genomes of the key tetrapod lineages, amphibians and reptiles, and in lungfishes. In this study, we combined genome database searches with PCR analysis to demonstrate that L1 retrotransposons are present in the genomes of lungfishes, amphibians, and lepidosaurs. Phylogenomic analysis shows that the genomes of Deuterostomia possess three highly divergent groups of L1 retrotransposons, with distinct distribution patterns. The analysis of L1 diversity shows the presence of a very large number of diverse L1 families, each with very low copy numbers, at the time of the origin of tetrapods. During the evolution of synapsids, all but one L1 lineage have been lost. This study establishes that the loss of L1 diversity and explosion in copy numbers occurred in the synapsid ancestors of mammals, and was most probably caused by severe population bottlenecks.     

Extensive envelope heterogeneity of simian immunodeficiency virus in tissues from infected macaques.

The extent of virus genetic variation within tissues and peripheral blood mononuclear cells (PBMC) from two simian immunodeficiency virus (SIV)-infected macaques was analyzed. The products of PCR amplification of two regions, region 1 (SIV V1 region) and region 2 (region corresponding to the human immunodeficiency virus V3 cysteine loop and part of the C3 region immediately downstream), of the SIV envelope were examined for single-stranded conformation polymorphism followed by sequence analysis of selected clones. The V1 region of the SIV envelope of viruses present within lymphoid tissues displayed extensive heterogeneity, while viral populations within the PBMC and brain appeared to be less variable. Region 2 heterogeneity in both animals was generally confined to three residues in a tissue-specific manner. In addition, virus from the brains of both animals appeared to be distinct compared with viruses present in other tissues and PBMC of the same animal, both in the pattern of PCR-single-stranded conformation polymorphism SCP and in the sequence of region 2. These studies revealed that the tissues of SIV-infected macaques were a reservoir for viral variants distinct from those seen in PBMC.     

Coreceptor Specificity of Temporal Variants of Simian Immunodeficiency Virus Mne

The simian immunodeficiency virus (SIV) Mne envelope undergoes genetic changes that alter tropism, syncytium-inducing capacity, and antigenic properties of the emerging variant virus population during the course of an infection. Here we investigated whether the mutations in envelope of SIVMne also influence coreceptor usage. The data demonstrate that the infecting macrophage-tropic SIVMne clone as well as the envelope variants that are selected during the course of disease progression all recognize both CCR5 and Bob (GPR15) but not Bonzo (STRL33), CXCR4, or CCR3. Although it remains to be determined if there are other coreceptors specific for dualtropic or T-cell-tropic variants of SIVMne that emerge during late stages of infection, these data suggest that such SIV variants that evolve in pathogenic infections do not lose the ability to recognize CCR5 or Bob/GPR15.