Novel HIV-1 long terminal repeat (LTR) sequences of subtype B and mosaic intersubtype B/C recombinants in North India

Revisiting the recombinant history of HIV-1 group M with dynamic network community detection

Recombination is a major mechanism through which HIV type 1 (HIV-1) maintains genetic diversity and interferes with viral eradication efforts. There is growing evidence demonstrating a recombinant origin of primate lentiviruses including HIV-1 group M (HIV-1/M). Inferring the extent of recombination across the entire HIV-1/M genome is of great importance as it provides deeper insights into the origin, dynamics, and evolution of the global pandemic. Here we propose an alternative method that can reconstruct the extent of genome-wide recombination in HIV-1, uncover reticulate patterns, and serve as a framework for HIV-1 classification. Our method provides an alternative approach for understanding the roles of virus recombination in the early evolutionary history of zoonosis for other emerging viruses.

The prevailing abundance of full-length HIV type 1 (HIV-1) genome sequences provides an opportunity to revisit the standard model of HIV-1 group M (HIV-1/M) diversity that clusters genomes into largely nonrecombinant subtypes, which is not consistent with recent evidence of deep recombinant histories for simian immunodeficiency virus (SIV) and other HIV-1 groups. Here we develop an unsupervised nonparametric clustering approach, which does not rely on predefined nonrecombinant genomes, by adapting a community detection method developed for dynamic social network analysis. We show that this method (dynamic stochastic block model [DSBM]) attains a significantly lower mean error rate in detecting recombinant breakpoints in simulated data (quasibinomial generalized linear model (GLM), P<8×10−8), compared to other reference-free recombination detection programs (genetic algorithm for recombination detection [GARD], recombination detection program 4 [RDP4], and RDP5). When this method was applied to a representative sample of n = 525 actual HIV-1 genomes, we determined k = 29 as the optimal number of DSBM clusters and used change-point detection to estimate that at least 95% of these genomes are recombinant. Further, we identified both known and undocumented recombination hotspots in the HIV-1 genome and evidence of intersubtype recombination in HIV-1 subtype reference genomes. We propose that clusters generated by DSBM can provide an informative framework for HIV-1 classification.

NF-κB Duplications in the Promoter-Variant HIV-1C LTR Impact Inflammation Without Altering Viral Replication in the Context of Simian Human Immunodeficiency Viruses and Opioid-Exposure

Recent spread of the promoter variant (4-κB) Human immunodeficiency virus-1 clade C (HIV-1C) strain is attributed to duplication of the Nuclear Factor Kappa B (NF-κB) binding sites and potential increased heroin consumption in India. To study the underlying biology of 4-κB HIV-1C in rhesus macaques, we engineered a promoter-chimera variant (4NF-κB) Simian Human Immunodeficiency Virus (SHIV) by substituting the HIV-1C Long Terminal Repeat (LTR) region consisting of 4 NF-κB and 3 Sp-1 sites with the corresponding segment in the LTR of SHIV AD8EO. The wild-type (3NF-κB) promoter-chimera SHIV was generated by inactivating the 5' proximal NF-κB binding site in SHIV 4NF-κB. CD8-depleted rhesus macaque PBMCs (RM-PBMCs) were infected with the promoter-chimera and AD8EO SHIVs to determine the effects of opioid-exposure on inflammation, NF-κB activation, neurotoxicity in neuronal cells and viral replication. Morphine-exposure of RM-PBMCs infected with SHIVs 4NF-κB, 3NF-κB, and AD8EO altered cellular transcript levels of monocyte chemoattractant protein 1, interleukin 6, interleukin 1β, and Tumor Necrosis Factor α. Of note, divergent alteration of the cytokine transcript levels was observed with these promoter-chimera wild-type and variant SHIVs. NF-κB activation was observed during infection of all three SHIVs with morphine-exposure. Finally, we observed that SHIV AD8EO infection and exposure to both morphine and naloxone had the greatest impact on the neurotoxicity. The promoter-chimera SHIV 4NF-κB and SHIV 3NF-κB did not have a similar effect on neurotoxicity as compared to SHIV AD8EO. All SHIVs replicated efficiently at comparable levels in RM-PBMCs and morphine-exposure did not alter viral replication kinetics. Future in vivo studies in rhesus macaques will provide greater understanding of 4-κB HIV-1C viral immunopathogenesis and onset of disease in the central nervous system during morphine-exposure.

Effect on HIV-1 gene expression, Tat-Vpr interaction and cell apoptosis by natural variants of HIV-1 Tat exon 1 and Vpr from Northern India

BACKGROUND

Since HIV-1 Tat and Vpr genes are involved in promoter transactivation, apoptosis, etc, we carried out studies to find out nature and extent of natural variation in the two genes from seropositive patients from Northern India and determined their functional implications.

METHODS

HIV-1 tat exon 1 and vpr were amplified from the genomic DNA isolated from the blood of HIV-1 infected individuals using specific primers by Polymerase Chain reaction (PCR) and subjected to extensive genetic analysis (CLUSTAL W, Simplot etc). Their expression was monitored by generating myc fusion clones. Tat exon 1 and Vpr variants were co-transfected with the reporter gene construct (LTR-luc) and their transactivation potential was monitored by measuring luciferase activity. Apoptosis and cell cycle analysis was done by Propidium Iodide (PI) staining followed by FACS.

RESULTS

Exon 1 of tat was amplified from 21 samples and vpr was amplified from 16 samples. One of the Tat exon 1 variants showed phylogenetic relatedness to subtype B & C and turned out to be a unique recombinant. Two of the Vpr variants were B/C/D recombinants. These natural variations were found to have no impact on the stability of Tat and Vpr. These variants differed in their ability to transactivate B LTR and C LTR promoters. B/C recombinant Tat showed better co-operative interaction with Vpr. B/C/D recombination in Vpr was found to have no effect on its co-operativity with Tat. Recombinant Tat (B/C) induced more apoptosis than wild type B and C Tat. The B/C/D recombination in Vpr did not affect its G2 arrest induction potential but reduced its apoptosis induction ability.

CONCLUSIONS

Extensive sequence and region-specific variations were observed in Tat and Vpr genes from HIV-1 infected individuals from Northern India. These variations have functional implications & therefore important for the pathogenicity of virus.

Molecular epidemiology of HIV-1 subtypes in India: origin and evolutionary history of the predominant subtype C

Background

India has the third largest HIV-1 epidemic with 2.4 million infected individuals. Molecular epidemiological analysis has identified the predominance of HIV-1 subtype C (HIV-1C). However, the previous reports have been limited by sample size, and uneven geographical distribution. The introduction of HIV-1C in India remains uncertain due to this lack of structured studies. To fill the gap, we characterised the distribution pattern of HIV-1 subtypes in India based on data collection from nationwide clinical cohorts between 2007 and 2011. We also reconstructed the time to the most recent common ancestor (tMRCA) of the predominant HIV-1C strains.

Methodology/Principal Findings

Blood samples were collected from 168 HIV-1 seropositive subjects from 7 different states. HIV-1 subtypes were determined using two or three genes, gag, pol, and env using several methods. Bayesian coalescent-based approach was used to reconstruct the time of introduction and population growth patterns of the Indian HIV-1C. For the first time, a high prevalence (10%) of unique recombinant forms (BC and A1C) was observed when two or three genes were used instead of one gene (p<0.01; p = 0.02, respectively). The tMRCA of Indian HIV-1C was estimated using the three viral genes, ranged from 1967 (gag) to 1974 (env). Pol-gene analysis was considered to provide the most reliable estimate [1971, (95% CI: 1965–1976)]. The population growth pattern revealed an initial slow growth phase in the mid-1970s, an exponential phase through the 1980s, and a stationary phase since the early 1990s.

Conclusions/Significance

The Indian HIV-1C epidemic originated around 40 years ago from a single or few genetically related African lineages, and since then largely evolved independently. The effective population size in the country has been broadly stable since the 1990s. The evolving viral epidemic, as indicated by the increase of recombinant strains, warrants a need for continued molecular surveillance to guide efficient disease intervention strategies.

Genetic architecture of HIV-1 genes circulating in north India & their functional implications

This review presents data on genetic and functional analysis of some of the HIV-1 genes derived from HIV-1 infected individuals from north India (Delhi, Punjab and Chandigarh). We found evidence of novel B/C recombinants in HIV-1 LTR region showing relatedness to China/Myanmar with 3 copies of Nfκb sites; B/C/D mosaic genomes for HIV-1 Vpr and novel B/C Tat. We reported appearance of a complex recombinant form CRF_02AG of HIV-1 envelope sequences which is predominantly found in Central/Western Africa. Also one Indian HIV-1 envelope subtype C sequence suggested exclusive CXCR4 co-receptor usage. This extensive recombination, which is observed in about 10 per cent HIV-1 infected individuals in the Vpr genes, resulted in remarkably altered functions when compared with prototype subtype B Vpr. The Vpu C was found to be more potent in causing apoptosis when compared with Vpu B when analyzed for subG1 DNA content. The functional implications of these changes as well as in other genes of HIV-1 are discussed in detail with possible implications for subtype-specific pathogenesis highlighted.

Diversity of HIV type 1 long terminal repeat (LTR) sequences following mother-to-child transmission in North India

We genetically characterized the extent of variation in HIV-1 LTR sequences from 11 mother-to-child transmission (MTCT) pairs from HIV-1-infected individuals from North India. Nine pairs were found to be infected with subtype C virus whereas two pairs were infected with subtype B virus. They harbored the characteristic three and two NF-κB sites, respectively. The analysis of intrasubtype divergence between B and C revealed greater diversity with subtype B LTR sequences than subtype C (p < 0.005). Significant evolutionary divergence of subtype C and subtype B was found in NFAT-III (p < 0.000001), NFAT-II (p < 0.0001), and USF (p < 0.005) transcription factor binding sites (TFBS). NF-κB-I, Sp I and II, Ets-I, AP-I and II, and TATA Box TFBS were highly conserved in both the subtypes. An alternate secondary structure of Tar was detected in the VT5 sample due to the point mutation from G to C (position +21) and T to C (position +38).

Genetic and functional analysis of HIV-1 Rev Responsive Element (RRE) sequences from North-India

HIV-1 Rev protein regulates the expression of HIV-1 transcripts by binding to a highly structured stem loop structure called the Rev Responsive Element (RRE) present in the genomic and partially spliced RNAs. Genetic variation in this structure is likely to affect binding of Rev protein and ultimately overall gene expression and replication. We characterized RRE sequences from 13 HIV-1 infected individuals from North India which also included two mother-child pairs following vertical transmission. We observed high degree of conservation of sequences, including the 9-nt (CACUAUGGG) long sequence in stem-loop B, required for efficient binding of Rev protein. All of our 13 RRE sequences possessed G to A (position 66) mutation located in the critical branched-stem-loop B which is not present in consensus C or B sequence. We derived a consensus RRE structure which showed interesting changes in the stem-loop structures including the stem-loop B. Mother-Child RRE sequences showed conservation of unique polymorphisms as well as some new mutations in child RRE sequences. Despite these changes, the ability to form multiple essential stem-loop structures required for Rev binding was conserved. RRE RNA derived from one of the samples, VT5, retained the ability to bind Rev protein under in vitro conditions although it showed alternate secondary structure. This is the first study from India describing the structural and possible functional implications due to very unique RRE sequence heterogeneity and its possible role in vertical transmission and gene expression.

Genetic analysis of HIV-1 Circulating Recombinant Form 02_AG, B and C subtype-specific envelope sequences from Northern India and their predicted co-receptor usage

HIV-1 epidemic in India is largely driven by subtype C but other subtypes or recombinants have also been reported from several states of India. This is mainly due to the co-circulation of other genetic subtypes that potentially can recombine to generate recombinant/mosaic genomes. In this study, we report detail genetic characterization of HIV-1 envelope sequences from North India (Delhi and neighboring regions). Six of 13 were related to subtype C, one B and the rest six showed relatedness with CRF02_AG strain. The subtype C possessed the highly conserved GPGQ motif but subtype B possessed the GPGR motif in the V3 loop as observed earlier. While most of the sequences suggested CCR5 co-receptor usage, one subtype C sample clearly indicated CXCR4 usage. A successful mother to child transmission was established in two pairs. Thus, co-circulation of multiple subtypes (B and C) and the recombinant CRF02_AG strains in North India suggests a rapidly evolving scenario of HIV-1 epidemic in this region with impact on vaccine formulation. Since this is the first report of CRF02_AG envelope from India, it will be important to monitor the spread of this strain and its impact on HIV-1 transmission in India.

Genetic and functional characterization of human immunodeficiency virus type 1 VprC variants from north India: presence of unique recombinants with mosaic genomes from B, C and D subtypes within the open reading frame of Vpr

The human immunodeficiency virus type 1 (HIV-1) epidemic in India is predominantly caused by genetic subtype C, though other minor subtypes have also been reported. One of the major accessory proteins of HIV-1, namely Vpr, is known to influence key steps in viral replication, cell cycle progression, promoter activation, apoptosis and pathogenesis. Therefore, we carried out a genetic and functional analysis of the Vpr variants from eight HIV-1-infected individuals from north India. The sequence analyses revealed that six of eight samples clustered with ancestral subtype C. Remarkably, five of these showed a conserved and region-specific L64P mutation, located in the predicted third α-helix. This change adversely affected their ability to activate the HIV-1 long terminal repeat promoter without compromising their ability to cause apoptosis. Bootscan, phylogenetic and SimPlot analysis of the remaining two samples (VprS2 and A6) revealed very interesting mosaic genomes derived from B, C and D subtypes. The N-terminal half of the VprS2 gene consisted of genomic segments derived from subtypes B/D, C and D but the C-terminal half was derived predominantly from subtype C. Interestingly the N-terminal half of sample A6 also showed similar B/D, C and D inter-subtype recombinant structure but the C-terminal half was entirely derived from the consensus B subtype. Multiple breakpoints in a short stretch of 291 nt encoding the Vpr gene strongly suggest that this region is a potential hot-spot for the formation of inter-subtype recombinants and also highlight the importance of the rapidly evolving HIV-1 epidemic in the north Indian region due to multiple genetic subtypes.

In-depth analysis of the origins of HIV type 1 subtype C in South America

The South American HIV-1 epidemic is characterized by the co-circulation of subtype B and BF recombinant variants. Together with the B and BF genotypes, HIV-1 subtype C (HIV-1C), F1, and several other recombinants have been reported. The epidemiological significance and immune correlates of these "non-B-non-BF" strains circulating in South America are still uncertain and therefore are increasingly attracting the interest of the scientific community. In this study, the South American HIV-1C epidemic was studied using new technologies for the phylogenetic analysis of large datasets. Our results indicate that there is a major clade encompassing most of the South American HIV-1C strains. These analyses also agreed that some strains do not group inside this major clade, suggesting that there could be HIV-1C sequences of different origins circulating in South America. Others have proposed different hypotheses about the origins of HIV-1C strains from South America. This study shows that an exact single origin cannot be determined, a fact that could be attributed to sampling problems, phylogenetic uncertainty, and the shortage of historical and epidemiological data. Currently, the reported data indicate that HIV-1C strains were introduced in Brazil and afterward spread to other regions of South America. By using character optimization on the obtained phylogenetic trees, we observed that Argentina could also be a point in which the HIV-1C epidemic entered South America.