Origin and dynamics of HIV-1 subtype C infection in India

Objective

To investigate the geographical origin and evolution dynamics of HIV-1 subtype C infection in India.

Design

Ninety HIV-1 subtype C env gp120 subtype C sequences from India were compared with 312 env gp120 reference subtype C sequences from 27 different countries obtained from Los Alamos HIV database. All the HIV-1 subtype C env gp120 sequences from India were used for the geographical origin analysis and 61 subtype C env gp120 sequences with known sampling year (from 1991 to 2008) were employed to determine the origin of HIV infection in India.

Methods

Phylogenetic analysis of HIV-1 env sequences was used to investigate the geographical origin and tMRCA of Indian HIV-1 subtype C. Evolutionary parameters including origin date and demographic growth patterns of Indian subtype C were estimated using a Bayesian coalescent-based approach under relaxed molecular clock models.

Findings

The majority of the analyzed Indian and South African HIV-1 subtype C sequences formed a single monophyletic cluster. The most recent common ancestor date was calculated to be 1975.56 (95% HPD, 1968.78–1981.52). Reconstruction of the effective population size revealed three phases of epidemic growth: an initial slow growth, followed by exponential growth, and then a plateau phase approaching present time. Stabilization of the epidemic growth phase correlated with the foundation of National AIDS Control Organization in India.

Interpretation

Indian subtype C originated from a single South African lineage in the middle of 1970s. The current study emphasizes not only the utility of HIV-1 sequence data for epidemiological studies but more notably highlights the effectiveness of community or government intervention strategies in controlling the trend of the epidemic.     

MicroRNA expression signatures of bladder cancer revealed by deep sequencing

Background

MicroRNAs (miRNAs) are a class of small noncoding RNAs that regulate gene expression. They are aberrantly expressed in many types of cancers. In this study, we determined the genome-wide miRNA profiles in bladder urothelial carcinoma by deep sequencing.

Methodology/Principal Findings

We detected 656 differentially expressed known human miRNAs and miRNA antisense sequences (miRNA*s) in nine bladder urothelial carcinoma patients by deep sequencing. Many miRNAs and miRNA*s were significantly upregulated or downregulated in bladder urothelial carcinoma compared to matched histologically normal urothelium. hsa-miR-96 was the most significantly upregulated miRNA and hsa-miR-490-5p was the most significantly downregulated one. Upregulated miRNAs were more common than downregulated ones. The hsa-miR-183, hsa-miR-200b∼429, hsa-miR-200c∼141 and hsa-miR-17∼92 clusters were significantly upregulated. The hsa-miR-143∼145 cluster was significantly downregulated. hsa-miR-182, hsa-miR-183, hsa-miR-200a, hsa-miR-143 and hsa-miR-195 were evaluated by Real-Time qPCR in a total of fifty-one bladder urothelial carcinoma patients. They were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium (p<0.001 for each miRNA).

Conclusions/Significance

To date, this is the first study to determine genome-wide miRNA expression patterns in human bladder urothelial carcinoma by deep sequencing. We found that a collection of miRNAs were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium, suggesting that they might play roles as oncogenes or tumor suppressors in the development and/or progression of this cancer. Our data provide novel insights into cancer biology.     

Human Immunodeficiency Virus Type-1 (HIV-1) Continues to Evolve in Presence of Broadly Neutralizing Antibodies More than Ten Years after Infection

Background

The evolution of HIV-1 and its immune escape to autologous neutralizing antibodies (Nabs) during the acute/early phases of infection have been analyzed in depth in many studies. In contrast, little is known about neither the long-term evolution of the virus in patients who developed broadly Nabs (bNabs) or the mechanism of escape in presence of these bNabs.

Results

We have studied the viral population infecting a long term non progressor HIV-1 infected patient who had developed broadly neutralizing antibodies toward all tier 2/3 viruses (6 clades) tested, 9 years after infection, and was then followed up over 7 years. The autologous neutralization titers of the sequential sera toward env variants representative of the viral population significantly increased during the follow-up period. The most resistant pseudotyped virus was identified at the last visit suggesting that it represented a late emerging escape variant. We identified 5 amino acids substitutions that appeared associated with escape to broadly neutralizing antibodies. They were V319I/S, R/K355T, R/W429G, Q460E and G/T463E, in V3, C3 and V5 regions.

Conclusion

This study showed that HIV-1 may continue to evolve in presence of both broadly neutralizing antibodies and increasing autologous neutralizing activity more than 10 years post-infection.     

Significantly longer envelope V2 loops are characteristic of heterosexually transmitted subtype B HIV-1 in Trinidad

Background

In Trinidad and the wider Caribbean, subtype B Human Immunodeficiency Virus-type 1 (HIV-1B) overwhelmingly accounts for HIV infection among heterosexuals; this contrasts with the association of HIV-1B with homosexual transmission and injecting drug use globally. The HIV envelope contains genetic determinants of cell tropism and evasion from immune attack. In this study we investigate the genetic properties of the env V1-C4 of HIV-1B soon after transmission to Trinidadian heterosexuals. This will reveal distinctive genetic features of the strains that cause the HIV-1B epidemic in Trinidad and generate insights to better understand their properties.

Methodology/Principal Findings

Quasispecies sampling was performed on the env V1-C4 of HIV-1B strains soon after transmission to heterosexual Trinidadians in a cohort of seroconverters. Phylogenetic relationships were determined for these quasispecies and the length and number of asparagine (N) linked glycosylation sites (NLGS) in their variable loops compared to that for HIV-1B globally. Signature amino acids within the constant domains of the env V1-C4 were identified for heterosexually transmitted HIV-1B from Trinidad relative to HIV-1B globally. HIV-1B obtained from Trinidadian heterosexuals soon after seroconversion had significantly longer V2 loops with one more glycosylation site, shorter V3 loops and no significant difference in V1 or V4 when compared to HIV-1B obtained soon after seroconversion from infected individuals in the rest of the world. HIV-1B soon after seroconversion and during chronic infection of Trinidadians was not significantly different, suggesting that distinctly long V2 loops are characteristic of HIV-1B in Trinidad. A threonine deletion at position 319 (T319-) along with the substitutions R315K and S440R were found to be distinctly associated with HIV-1B from Trinidad compared to HIV-1B globally.

Conclusions

This finding of distinctive genetic features that are characteristic of HIV-1B strains from Trinidad is consistent with the Trinidad epidemic being established by a founder strain or closely related founder strains of HIV-1B.     

Phylogenetic Relatedness of Circulating HIV-1C Variants in Mochudi,Botswana

Background

Determining patterns of HIV transmission is increasingly important for the most efficient use of modern prevention interventions. HIV phylogeny can provide a better understanding of the mechanisms underlying HIV transmission networks in communities.

Methods

To reconstruct the structure and dynamics of a local HIV/AIDS epidemic, the phylogenetic relatedness of HIV-1 subtype C env sequences obtained from 785 HIV-infected community residents in the northeastern sector of Mochudi, Botswana, during 2010–2013 was estimated. The genotyping coverage was estimated at 44%. Clusters were defined based on relatedness of HIV-1C env sequences and bootstrap support of splits.

Results

The overall proportion of clustered HIV-1C env sequences was 19.1% (95% CI 17.5% to 20.8%). The proportion of clustered sequences from Mochudi was significantly higher than the proportion of non-Mochudi sequences that clustered, 27.0% vs. 14.7% (p = 5.8E-12; Fisher exact test). The majority of clustered Mochudi sequences (90.1%; 95% CI 85.1% to 93.6%) were found in the Mochudi-unique clusters. None of the sequences from Mochudi clustered with any of the 1,244 non-Botswana HIV-1C sequences. At least 83 distinct HIV-1C variants, or chains of HIV transmission, in Mochudi were enumerated, and their sequence signatures were reconstructed. Seven of 20 genotyped seroconverters were found in 7 distinct clusters.

Conclusions

The study provides essential characteristics of the HIV transmission network in a community in Botswana, suggests the importance of high sampling coverage, and highlights the need for broad HIV genotyping to determine the spread of community-unique and community-mixed viral variants circulating in local epidemics. The proposed methodology of cluster analysis enumerates circulating HIV variants and can work well for surveillance of HIV transmission networks. HIV genotyping at the community level can help to optimize and balance HIV prevention strategies in trials and combined intervention packages.     

Molecular Epidemiology of HIV-1 in Jilin Province,Northeastern China: Emergence of a New CRF07_BC Transmission Cluster and Intersubtype Recombinants

Objective

To investigate the HIV-1 molecular epidemiology among newly diagnosed HIV-1 infected persons living in the Jilin province of northeastern China.

Methods

Plasma samples from 189 newly diagnosed HIV-1 infected patients were collected between June 2010 and August 2011 from all nine cities of Jilin province. HIV-1 nucleotide sequences of gag P17–P24 and env C2–C4 gene regions were amplified using a multiplex RT-PCR method and sequenced. Phylogenetic and recombination analyses were used to determine the HIV-1 genotypes.

Results

Based on all sequences generated, the subtype/CFR distribution was as follows: CRF01_AE (58.1%), CRF07_BC (13.2%), subtype B’ (13.2%), recombinant viruses (8.1%), subtype B (3.7%), CRF02_AG (2.9%), subtype C (0.7%). In addition to finding CRF01_AE strains from previously reported transmission clusters 1, 4 and 5, a new transmission cluster was described within the CRF07_BC radiation. Among 11 different recombinants identified, 10 contained portions of gene regions from the CRF01_AE lineage. CRF02_AG was found to form a transmission cluster of 4 in local Jilin residents.

Conclusions

Our study presents a molecular epidemiologic investigation describing the complex structure of HIV-1 strains co-circulating in Jilin province. The results highlight the critical importance of continuous monitoring of HIV-infections, along with detailed socio-demographic data, in order to design appropriate prevention measures to limit the spread of new HIV infections.     

High-Throughput Sequencing of Plasma MicroRNA in Chronic Fatigue Syndrome/Myalgic Encephalomyelitis

Background

MicroRNAs (miRNAs) are known to regulate many biological processes and their dysregulation has been associated with a variety of diseases including Chronic Fatigue Syndrome/Myalgic Encephalomyelitis (CFS/ME). The recent discovery of stable and reproducible miRNA in plasma has raised the possibility that circulating miRNAs may serve as novel diagnostic markers. The objective of this study was to determine the role of plasma miRNA in CFS/ME.

Results

Using Illumina high-throughput sequencing we identified 19 miRNAs that were differentially expressed in the plasma of CFS/ME patients in comparison to non-fatigued controls. Following RT-qPCR analysis, we were able to confirm the significant up-regulation of three miRNAs (hsa-miR-127-3p, hsa-miR-142-5p and hsa-miR-143-3p) in the CFS/ME patients.

Conclusion

Our study is the first to identify circulating miRNAs from CFS/ME patients and also to confirm three differentially expressed circulating miRNAs in CFS/ME patients, providing a basis for further study to find useful CFS/ME biomarkers.     

Human T-cell leukemia virus type I infects human lung epithelial cells and induces gene expression of cytokines, chemokines and cell adhesion molecules

Background

Cellular miRNAs play an important role in the regulation of gene expression in eukaryotes. Recently, miRNAs have also been shown to be able to target and inhibit viral gene expression. Computational predictions revealed earlier that the HIV-1 genome includes regions that may be potentially targeted by human miRNAs. Here we report the functionality of predicted miR-29a target site in the HIV-1 nef gene.

Results

We find that the human miRNAs hsa-miR-29a and 29b are expressed in human peripheral blood mononuclear cells. Expression of a luciferase reporter bearing the nef miR-29a target site was decreased compared to the luciferase construct without the target site. Locked nucleic acid modified anti-miRNAs targeted against hsa-miR-29a and 29b specifically reversed the inhibitory effect mediated by cellular miRNAs on the target site. Ectopic expression of the miRNA results in repression of the target Nef protein and reduction of virus levels.

Conclusion

Our results show that the cellular miRNA hsa-miR29a downregulates the expression of Nef protein and interferes with HIV-1 replication.     

Restricted Genetic Diversity of HIV-1 Subtype C Envelope Glycoprotein from Perinatally Infected Zambian Infants

Background

Mother-to-child transmission of HIV-1 remains a significant problem in the resource-constrained settings where anti-retroviral therapy is still not widely available. Understanding the earliest events during HIV-1 transmission and characterizing the newly transmitted or founder virus is central to intervention efforts. In this study, we analyzed the viral env quasispecies of six mother-infant transmission pairs (MIPs) and characterized the genetic features of envelope glycoprotein that could influence HIV-1 subtype C perinatal transmission.

Methodology and Findings

The V1-V5 region of env was amplified from 6 MIPs baseline samples and 334 DNA sequences in total were analyzed. A comparison of the viral population derived from the mother and infant revealed a severe genetic bottleneck occurring during perinatal transmission, which was characterized by low sequence diversity in the infant. Phylogenetic analysis indicates that most likely in all our infant subjects a single founder virus was responsible for establishing infection. Furthermore, the newly transmitted viruses from the infant had significantly fewer potential N-linked glycosylation sites in Env V1-V5 region and showed a propensity to encode shorter variable loops compared to the nontransmitted viruses. In addition, a similar intensity of selection was seen between mothers and infants with a higher rate of synonymous (dS) compared to nonsynonymous (dN) substitutions evident (dN/dS<1).

Conclusions

Our results indicate that a strong genetic bottleneck occurs during perinatal transmission of HIV-1 subtype C. This is evident through population diversity and phylogenetic patterns where a single viral variant appears to be responsible for infection in the infants. As a result the newly transmitted viruses are less diverse and harbored significantly less glycosylated envelope. This suggests that viruses with the restricted glycosylation in envelope glycoprotein appeared to be preferentially transmitted during HIV-1 subtype C perinatal transmission. In addition, our findings also indicated that purifying selection appears to predominate in shaping the early intrahost evolution of HIV-1 subtype C envelope sequences.     

Subtype Classification of Iranian HIV-1 Sequences Registered in the HIV Databases, 2006-2013

Background

The rate of human immunodeficiency virus type 1 (HIV-1) infection in Iran has increased dramatically in the past few years. While the earliest cases were among hemophiliacs, injection drug users (IDUs) fuel the current epidemic. Previous molecular epidemiological analysis found that subtype A was most common among IDUs but more recent studies suggest CRF_35AD may be more prevalent now. To gain a better understanding of the molecular epidemiology of HIV-1 infection in Iran, we analyzed all Iranian HIV sequence data from the Los Alamos National Laboratory.

Methods

All Iranian HIV sequences from subtyping studies with pol, gag, env and full-length HIV-1 genome sequences registered in the HIV databases (www.hiv.lanl.gov) between 2006 and 2013 were downloaded. Phylogenetic trees of each region were constructed using Neighbor-Joining (NJ) and Maximum Parsimony methods.

Results

A total of 475 HIV sequences were analyzed. Overall, 78% of sequences were CRF_35AD. By gene region, CRF_35AD comprised 83% of HIV-1 pol, 62% of env, 78% of gag, and 90% of full-length genome sequences analyzed. There were 240 sequences re-categorized as CRF_AD. The proportion of CRF_35AD sequences categorized by the present study is nearly double the proportion of what had been reported.

Conclusions

Phylogenetic analysis indicates HIV-1 subtype CRF_35AD is the predominant circulating strain in Iran. This result differed from previous studies that reported subtype A as most prevalent in HIV- infected patients but confirmed other studies which reported CRF_35AD as predominant among IDUs. The observed epidemiological connection between HIV strains circulating in Iran and Afghanistan may be due to drug trafficking and/or immigration between the two countries. This finding suggests the possible origins and transmission dynamics of HIV/AIDS within Iran and provides useful information for designing control and intervention strategies.     

Evidence of at least two introductions of HIV-1 in the Amerindian Warao population from Venezuela

Background

The Venezuelan Amerindians were, until recently, free of human immunodeficiency virus (HIV) infection. However, in 2007, HIV-1 infection was detected for the first time in the Warao Amerindian population living in the Eastern part of Venezuela, in the delta of the Orinoco river. The aim of this study was to analyze the genetic diversity of the HIV-1 circulating in this population.

Methodology/Principal Findings

The pol genomic region was sequenced for 16 HIV-1 isolates and for some of them, sequences from env, vif and nef genomic regions were obtained. All HIV-1 isolates were classified as subtype B, with exception of one that was classified as subtype C. The 15 subtype B isolates exhibited a high degree of genetic similarity and formed a highly supported monophyletic cluster in each genomic region analyzed. Evolutionary analyses of the pol genomic region indicated that the date of the most recent common ancestor of the Waraos subtype B clade dates back to the late 1990s.

Conclusions/Significance

At least two independent introductions of HIV-1 have occurred in the Warao Amerindians from Venezuela. The HIV-1 subtype B was successfully established and got disseminated in the community, while no evidence of local dissemination of the HIV-1 subtype C was detected in this study. These results warrant further surveys to evaluate the burden of this disease, which can be particularly devastating in this Amerindian population, with a high prevalence of tuberculosis, hepatitis B, among other infectious diseases, and with limited access to primary health care.     

Evidence for Limited Genetic Compartmentalization of HIV-1 between Lung and Blood

Background

HIV-1 is frequently detected in the lungs of infected individuals and is likely important in the development of pulmonary opportunistic infections. The unique environment of the lung, rich in alveolar macrophages and with specialized local immune responses, may contribute to differential evolution or selection of HIV-1.

Methodology and Findings

We characterized HIV-1 in the lung in relation to contemporaneous viral populations in the blood. The C2-V5 region of HIV-1 env was sequenced from paired lung (induced sputum or bronchoalveolar lavage) and blood (plasma RNA and proviral DNA from sorted or unsorted PBMC) from 18 subjects. Compartmentalization between tissue pairs was assessed using 5 established tree or distance-based methods, including permutation tests to determine statistical significance. We found statistical evidence of compartmentalization between lung and blood in 10/18 subjects, although lung and blood sequences were intermingled on phylogenetic trees in all subjects. The subject showing the greatest compartmentalization contained many nearly identical sequences in BAL sample, suggesting clonal expansion may contribute to reduced viral diversity in the lung in some cases. However, HIV-1 sequences in lung were not more homogeneous overall, nor were we able to find a lung-specific genotype associated with macrophage tropism in V3. In all four subjects in whom predicted X4 genotypes were found in blood, predicted X4 genotypes were also found in lung.

Conclusions

Our results support a picture of continuous migration of HIV-1 between circulating blood and lung tissue, with perhaps a very limited degree of localized evolution or clonal replication.     

De-Regulated MicroRNAs in Pediatric Cancer Stem Cells Target Pathways Involved in Cell Proliferation,Cell Cycle and Development

Background

microRNAs (miRNAs) have been implicated in the control of many biological processes and their deregulation has been associated with many cancers. In recent years, the cancer stem cell (CSC) concept has been applied to many cancers including pediatric. We hypothesized that a common signature of deregulated miRNAs in the CSCs fraction may explain the disrupted signaling pathways in CSCs.

Methodology/Results

Using a high throughput qPCR approach we identified 26 CSC associated differentially expressed miRNAs (DEmiRs). Using BCmicrO algorithm 865 potential CSC associated DEmiR targets were obtained. These potential targets were subjected to KEGG, Biocarta and Gene Ontology pathway and biological processes analysis. Four annotated pathways were enriched: cell cycle, cell proliferation, p53 and TGF-beta/BMP. Knocking down hsa-miR-21-5p, hsa-miR-181c-5p and hsa-miR-135b-5p using antisense oligonucleotides and small interfering RNA in cell lines led to the depletion of the CSC fraction and impairment of sphere formation (CSC surrogate assays).

Conclusion

Our findings indicated that CSC associated DEmiRs and the putative pathways they regulate may have potential therapeutic applications in pediatric cancers.     

Integrated miRNA and mRNA Expression Profiling in Inflamed Colon of Patients with Ulcerative Colitis

Background

Ulcerative colitis (UC) is associated with differential colonic expression of genes involved in immune response (e.g. IL8) and barrier integrity (e.g. cadherins). MicroRNAs (miRNAs) are regulators of gene expression and are involved in various immune-related diseases. In this study, we investigated (1) if miRNA expression in UC mucosa is altered and (2) if any of these changes correlate with mucosal mRNA expression. Integration of mRNA and miRNA expression profiling may allow the identification of functional links between dysregulated miRNAs and their target mRNA.

Methodology

Colonic mucosal biopsies were obtained from 17 UC (10 active and 7 inactive) patients and 10 normal controls. Total RNA was used to analyze miRNA and mRNA expression via Affymetrix miRNA 2.0 and Affymetrix Human Gene 1.0ST arrays, respectively. Both miRNA and gene expression profiles were integrated by correlation analysis to identify dysregulated miRNAs with their corresponding predicted target mRNA. Microarray data were validated with qRT-PCR. Regulation of IL8 and CDH11 expression by hsa-miR-200c-3p was determined by luciferase reporter assays.

Results

When comparing active UC patients vs. controls, 51 miRNAs and 1543 gene probe sets gave significantly different signals. In contrast, in inactive UC vs. controls, no significant miRNA expression differences were found while 155 gene probe sets had significantly different signals. We then identified potential target genes of the significantly dysregulated miRNAs and genes in active UC vs. controls and found a highly significant inverse correlation between hsa-miR-200c-3p and IL8, an inflammatory marker, and between hsa-miR-200c-3p and CDH11, a gene related to intestinal epithelial barrier function. We could demonstrate that hsa-miR-200c-3p directly regulates IL8 and CDH11 expression.

Conclusion

Differential expression of immune- and barrier-related genes in inflamed UC mucosa may be influenced by altered expression of miRNAs. Integrated analysis of miRNA and mRNA expression profiles revealed hsa-miR-200c-3p for use of miRNA mimics as therapeutics.     

Re-Inspection of Small RNA Sequence Datasets Reveals Several Novel Human miRNA Genes

Background

miRNAs are key players in gene expression regulation. To fully understand the complex nature of cellular differentiation or initiation and progression of disease, it is important to assess the expression patterns of as many miRNAs as possible. Thereby, identifying novel miRNAs is an essential prerequisite to make possible a comprehensive and coherent understanding of cellular biology.

Methodology/Principal Findings

Based on two extensive, but previously published, small RNA sequence datasets from human embryonic stem cells and human embroid bodies, respectively [1], we identified 112 novel miRNA-like structures and were able to validate miRNA processing in 12 out of 17 investigated cases. Several miRNA candidates were furthermore substantiated by including additional available small RNA datasets, thereby demonstrating the power of combining datasets to identify miRNAs that otherwise may be assigned as experimental noise.

Conclusions/Significance

Our analysis highlights that existing datasets are not yet exhaustedly studied and continuous re-analysis of the available data is important to uncover all features of small RNA sequencing.     

Molecular evolution of HIV-1 CRF01_AE Env in Thai patients

Background

The envelope glycoproteins (Env), gp120 and gp41, are the most variable proteins of human immunodeficiency virus type 1 (HIV-1), and are the major targets of humoral immune responses against HIV-1. A circulating recombinant form of HIV-1, CRF01_AE, is prevalent throughout Southeast Asia; however, only limited information regarding the immunological characteristics of CRF01_AE Env is currently available. In this study, we attempted to examine the evolutionary pattern of CRF01_AE Env under the selection pressure of host immune responses.

Methodology/Principal Findings

Peripheral blood samples were collected periodically over 3 years from 15 HIV-1-infected individuals residing in northern Thailand, and amplified env genes from the samples were subjected to computational analysis. The V5 region of gp120 showed highest variability in several samples over 3 years, whereas the V1/V2 and/or V4 regions of gp120 also showed high variability in many samples. In addition, the N-terminal part of the C3 region of gp120 showed highest amino acid diversity among the conserved regions of gp120. Chronological changes in the numbers of amino acid residues in gp120 variable regions and potential N-linked glycosylation (PNLG) sites are involved in increasing the variability of Env gp120. Furthermore, the C3 region contained several amino acid residues potentially under positive selection, and APOBEC3 family protein-mediated G to A mutations were frequently detected in such residues.

Conclusions/Significance

Several factors, including amino acid substitutions particularly in gp120 C3 and V5 regions as well as changes in the number of PNLG sites and in the length of gp120 variable regions, were revealed to be involved in the molecular evolution of CRF01_AE Env. In addition, a similar tendency was observed between CRF01_AE and subtype C Env with regard to the amino acid variation of gp120 V3 and C3 regions. These results may provide important information for understanding the immunological characteristics of CRF01_AE Env.