A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

Origin of Minority Drug-Resistant HIV-1 Variants in Primary HIV-1 Infection

Background. Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants (MVs) are present in some antiretroviral therapy (ART)–naive patients. They may result from de novo mutagenesis or transmission. To date, the latter has not been proven. Methods. MVs were quantified by allele-specific polymerase chain reaction in 204 acute or recent seroconverters from the Zurich Primary HIV Infection study and 382 ART-naive, chronically infected patients. Phylogenetic analyses identified transmission clusters. Results. Three lines of evidence were observed in support of transmission of MVs. First, potential transmitters were identified for 12 of 16 acute or recent seroconverters harboring M184V MVs. These variants were also detected in plasma and/or peripheral blood mononuclear cells at the estimated time of transmission in 3 of 4 potential transmitters who experienced virological failure accompanied by the selection of the M184V mutation before transmission. Second, prevalence between MVs harboring the frequent mutation M184V and the particularly uncommon integrase mutation N155H differed highly significantly in acute or recent seroconverters (8.2% vs 0.5%; P < .001). Third, the prevalence of less-fit M184V MVs is significantly higher in acutely or recently than in chronically HIV-1–infected patients (8.2% vs 2.5%; P = .004). Conclusions. Drug-resistant HIV-1 MVs can be transmitted. To what extent the origin—transmission vs sporadic appearance—of these variants determines their impact on ART needs to be further explored.

Impact of Switching From Zidovudine to Tenofovir Disoproxil Fumarate on Bone Mineral Density and Markers of Bone Metabolism in Virologically Suppressed HIV-1 Infected Patients; A Substudy of the PREPARE Study

Context: In virologically suppressed, antiretroviral-treated patients, the effect of switching to tenofovir (TDF) on bone biomarkers compared to patients remaining on stable antiretroviral therapy is unknown. Methods: We examined bone biomarkers (osteocalcin [OC], procollagen type 1 amino-terminal propeptide, and C-terminal cross-linking telopeptide of type 1 collagen) and bone mineral density (BMD) over 48 weeks in virologically suppressed patients (HIV RNA < 50 copies/ml) randomized to switch to TDF/emtricitabine (FTC) or remain on first-line zidovudine (AZT)/lamivudine (3TC). PTH was also measured. Between-group differences in bone biomarkers and associations between change in bone biomarkers and BMD measures were assessed by Student's t tests, Pearson correlation, and multivariable linear regression, respectively. All data are expressed as mean (SD), unless otherwise specified. Results: Of 53 subjects (aged 46.0 y; 84.9% male; 75.5% Caucasian), 29 switched to TDF/FTC. There were reductions in total hip and lumbar spine BMD in those switching to TDF/FTC (total hip, TDF/FTC, −1.73 (2.76)% vs AZT/3TC, −0.39 (2.41)%; between-group P = .07; lumbar spine, TDF/FTC, −1.50 (3.49)% vs AZT/3TC, +0.25 (2.82)%; between-group P = .06), but they did not reach statistical significance. Greater declines in lumbar spine BMD correlated with greater increases in OC (r = −0.28; P = .05). The effect of TDF/FTC on bone biomarkers remained significant when adjusted for baseline biomarker levels, gender, and ethnicity. There was no difference in change in PTH levels over 48 weeks between treatment groups (between-group P = .23). All biomarkers increased significantly from weeks 0 to 48 in the switch group, with no significant change in those remaining on AZT/3TC (between-group, all biomarkers, P < .0001). Conclusion: A switch to TDF/FTC compared to remaining on a stable regimen is associated with increases in bone turnover that correlate with reductions in BMD, suggesting that TDF exposure directly affects bone metabolism in vivo.

Higher CNS Penetration-Effectiveness of Long-term Combination Antiretroviral Therapy Is Associated With Better HIV-1 Viral Suppression in Cerebrospinal Fluid

Objectives: To determine HIV-1 RNA in cerebrospinal fluid (CSF) of successfully treated patients and to evaluate if combination antiretroviral treatments with higher central nervous system penetration-effectiveness (CPE) achieve better CSF viral suppression. Methods: Viral loads (VLs) and drug concentrations of lopinavir, atazanavir, and efavirenz were measured in plasma and CSF. The CPE was calculated using 2 different methods. Results: The authors analyzed 87 CSF samples of 60 patients. In 4 CSF samples, HIV-1 RNA was detectable with 43–82 copies per milliliter. Median CPE in patients with detectable CSF VL was significantly lower compared with individuals with undetectable VL: CPE of 1.0 (range, 1.0–1.5) versus 2.3 (range, 1.0–3.5) using the method of 2008 (P = 0.011) and CPE of 6 (range, 6–8) versus 8 (range, 5–12) using the method of 2010 (P = 0.022). The extrapolated CSF trough levels for atazanavir (n = 12) were clearly above the 50% inhibitory concentration (IC50) in only 25% of samples; both patients on atazanavir/ritonavir with detectable CSF HIV-1 RNA had trough levels in the range of the presumed IC50. The extrapolated CSF trough level for lopinavir (n = 42) and efavirenz (n = 18) were above the IC50 in 98% and 78%, respectively, of samples, including the patients with detectable CSF HIV-1 RNA. Conclusions: This study suggests that treatment regimens with high intracerebral efficacy reflected by a high CPE score are essential to achieve CSF HIV-1 RNA suppression. The CPE score including all drug components was a better predictor for treatment failure in the CSF than the sole concentrations of protease inhibitor or nonnucleoside reverse transcriptase inhibitor in plasma or CSF.

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals

Objectives: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. Materials and methods: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. Results: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13–69%] and a 42% (95% CI: 17–68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8–38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. Conclusion: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.

Effects of Alpha Interferon Treatment on Intrinsic Anti-HIV-1 Immunity In Vivo

Alpha interferon (IFN-α) suppresses human immunodeficiency virus type 1 (HIV-1) replication in vitro by inducing cell-intrinsic retroviral restriction mechanisms. We investigated the effects of IFN-α/ribavirin (IFN-α/riba) treatment on 34 anti-HIV-1 restriction factors in vivo. Expression of several anti-HIV-1 restriction factors was significantly induced by IFN-α/riba in HIV/hepatitis C virus (HCV)-coinfected individuals. Fold induction of cumulative restriction factor expression in CD4+ T cells was significantly correlated with viral load reduction during IFN-α/riba treatment (r2 = 0.649; P < 0.016). Exogenous IFN-α induces supraphysiologic restriction factor expression associated with a pronounced decrease in HIV-1 viremia.

THE ROLE OF MUCOSAL EPITHELIAL CELLS IN HIV-1 INFECTION

The predominant route of human immunodeficiency virus type 1 (HIV-1) transmission is infection across the vaginal mucosa. Epithelial cells, which form the primary barrier of protection against pathogens, are the first cell type at these mucosal tissues to encounter the virus but their role in HIV infection has not been clearly elucidated. Although mucosal epithelial cells express only low levels of the receptors required for successful HIV infection, productive infection does occur at these sites. The present work provides evidence to show that HIV exposure, without the need for productive infection, induces human cervical epithelial cells to produce Thymic Stromal Lymphopoietin (TSLP), an IL7-like cytokine, which potently activated human myeloid dendritic cells (mDC) to cause the homeostatic proliferation of autologous CD4+ T cells that serve as targets for HIV infection. Rhesus macaques inoculated with simian immunodeficiency virus (SIV) or with the simian-human immunodeficiency virus (SHIV) by the vaginal, oral or rectal route exhibited dramatic increases in: TSLP expression, DC and CD4+ T cell numbers, and viral replication, in the vaginal, oral, and rectal tissues, respectively within the first 2 weeks after virus exposure. Evidence obtained showed that HIV-mediated TSLP production by cervical cells is dependent upon the expression of the cell surface salivary agglutinin (SAG) protein gp340. Epithelial cells expressing gp340 exhibited HIV endocytosis and TSLP expression and genetic knockdown of gp340 or use of a gp340-blocking antibody inhibited TSLP expression by HIV. On the other hand, gp340-null epithelial cells failed to endocytose HIV and produce TSLP, but transfection of gp340 resulted in HIV-induced TSLP expression. Finally, HIV-induced TSLP expression was found to be mediated by TLR7/8 signaling and NF-kB activity because silencing these pathways or use of specific inhibitors abrogated TSLP expression in gp340-postive but not in gp340-null epithelial cells. Overall these studies identify TSLP as a key player in the acute phase of HIV-1 infection in permitting HIV to successfully maneuver the hostile vaginal mucosal microenvironment by creating a conducive environment for sustaining the small amount of virus that initially crosses the mucosal barrier allowing it to successfully cause infection and spread to distal compartments of the body

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

Highly pathogenic adapted HIV-1 strains limit host immunity and dictate rapid disease progression.

OBJECTIVE: The study of HIV-1 rapid progressors has been limited to specific case reports. Nevertheless, identification and characterization of the viral and host factors involved in rapid progression are crucial when attempting to uncover the correlates of rapid disease outcome. DESIGN: We carried out comparative functional analyses in rapid progressors (n = 46) and standard progressors (n = 46) early after HIV-1 seroconversion (≤1 year). The viral traits tested were viral replicative capacity, co-receptor usage, and genomic variation. Host CD8 T-cell responses, humoral activity, and HLA immunogenetic markers were also determined. RESULTS: Our data demonstrate an unusual convergence of highly pathogenic HIV-1 strains in rapid progressors. Compared with standard progressors, rapid progressor viral strains show higher in-vitro replicative capacity (81.5 vs. 67.9%; P = 0.025) and greater X4/DM co-receptor usage (26.3 vs. 2.8%; P = 0.006) in early infection. Limited or absent functional HIV-1 CD8 T-cell responses and neutralizing activity were measured in rapid progressors. Moreover, the increase in common HLA allele-restricted CD8 T-cell escape mutations in rapid progressors acts as a signature of uncontrolled HIV-1 replication and early impairment of adaptive cellular responses. CONCLUSION: Our data support a dominant role for viral factors in rapid progressors. Robust HIV-1 replication and intrinsic viral properties limit host adaptive immune responses, thus driving rapid disease progression.

PCR identification of HIV-1 DNA sequences in brain tissue of patients with AIDS encephalopathy.

We analyzed brain tissue from 39 patients for the presence of proviral HIV-1 sequences, using the polymerase chain reaction (PCR) for the amplification of segments of the viral LTR and gag genes. A novel primer extension procedure allowed the detection of a single HIV-1 copy in 1 micrograms DNA. We detected proviral HIV-1 DNA in 16 of 25 brain samples from AIDS patients. Semiquantitative evaluation of the amplified DNAs indicated considerable variation in viral load. Highest levels of proviral DNA were present in brain samples from six patients with clinical evidence of HIV-associated cognitive/motor complex and the histopathologic correlate of HIV leukoencephalopathy or HIV encephalitis. An additional 11 brain samples contained smaller amounts of proviral DNA. In these patients, clinical data were inconclusive regarding the diagnosis of HIV-1 encephalopathy and histopathologically there was no evidence of HIV-1-induced tissue lesions. In nine of 25 seropositive patients with AIDS (36%), brain samples scored negative or did not contain an unequivocal signal indicating the presence of proviral DNA. HIV-1 sequences were not detected in any of 14 control brain samples from HIV-1 seronegative patients. Our data indicate that HIV-1 is present in the central nervous system of the majority (two thirds) of AIDS patients and that the highest levels of proviral DNA in brain tissue are associated with HIV encephalopathy.

Social Meets Molecular: Combining Phylogenetic and Latent Class Analyses to Understand HIV-1 Transmission in Switzerland

Switzerland has a complex human immunodeficiency virus (HIV) epidemic involving several populations. We examined transmission of HIV type 1 (HIV-1) in a national cohort study. Latent class analysis was used to identify socioeconomic and behavioral groups among 6,027 patients enrolled in the Swiss HIV Cohort Study between 2000 and 2011. Phylogenetic analysis of sequence data, available for 4,013 patients, was used to identify transmission clusters. Concordance between sociobehavioral groups and transmission clusters was assessed in correlation and multiple correspondence analyses. A total of 2,696 patients were infected with subtype B, 203 with subtype C, 196 with subtype A, and 733 with recombinant subtypes (mainly CRF02_AG and CRF01_AE). Latent class analysis identified 8 patient groups. Most transmission clusters of subtype B were shared between groups of gay men (groups 1-3) or between the heterosexual groups "heterosexual people of lower socioeconomic position" (group 4) and "injection drug users" (group 8). Clusters linking homosexual and heterosexual groups were associated with "older heterosexual and gay people on welfare" (group 5). "Migrant women in heterosexual partnerships" (group 6) and "heterosexual migrants on welfare" (group 7) shared non-B clusters with groups 4 and 5. Combining approaches from social and molecular epidemiology can provide insights into HIV-1 transmission and inform the design of prevention strategies.

Virological outcome and management of persistent low-level viraemia in HIV-1-infected patients: 11 years of the Swiss HIV Cohort Study

BACKGROUND Management of persistent low-level viraemia (pLLV) in patients on combined antiretroviral therapy (cART) with previously undetectable HIV viral loads (VLs) is challenging. We examined virological outcome and management among patients enrolled in the Swiss HIV Cohort Study (SHCS). METHODS In this retrospective study (2000-2011), pLLV was defined as a VL of 21-400 copies/mL on ≥3 consecutive plasma samples with ≥8 weeks between first and last analyses, in patients undetectable for ≥24 weeks on cART. Control patients had ≥3 consecutive undetectable VLs over ≥32 weeks. Virological failure (VF), analysed in the pLLV patient group, was defined as a VL>400 copies/mL. RESULTS Among 9972 patients, 179 had pLLV and 5389 were controls. Compared to controls, pLLV patients were more often on unboosted PI-based (adjusted odds ratio, aOR, [95%CI] 3.2 [1.8-5.9]) and NRTI-only combinations (aOR 2.1 [1.1-4.2]) than on NNRTI and boosted PI-based regimens. At 48 weeks, 102/155 pLLV patients (66%) still had pLLV, 19/155 (12%) developed VF, and 34/155 (22%) had undetectable VLs. Predictors of VF were previous VF (aOR 35 [3.8-315]), unboosted PI-based (aOR 12.8 [1.7-96]) or NRTI-only combinations (aOR 115 [6.8-1952]), and VLs>200 during pLLV (aOR 3.7 [1.1-12]). No VF occurred in patients with persistent very LLV (pVLLV, 21-49 copies/mL; N=26). At 48 weeks, 29/39 patients (74%) who changed cART had undetectable VLs, compared to 19/74 (26%) without change (P<0.001). CONCLUSIONS Among patients with pLLV, VF was predicted by previous VF, cART regimen and VL ≥200. Most patients who changed cART had undetectable VLs 48 weeks later. These findings support cART modification for pLLV >200 copies/ml.

Role of MicroRNA Modulation in the Interferon-α/Ribavirin Suppression of HIV-1 In Vivo.

BACKGROUND Interferon-α (IFN-α) treatment suppresses HIV-1 viremia and reduces the size of the HIV-1 latent reservoir. Therefore, investigation of the molecular and immunologic effects of IFN-α may provide insights that contribute to the development of novel prophylactic, therapeutic and curative strategies for HIV-1 infection. In this study, we hypothesized that microRNAs (miRNAs) contribute to the IFN-α-mediated suppression of HIV-1. To inform the development of novel miRNA-based antiretroviral strategies, we investigated the effects of exogenous IFN-α treatment on global miRNA expression profile, HIV-1 viremia, and potential regulatory networks between miRNAs and cell-intrinsic anti-HIV-1 host factors in vivo. METHODS Global miRNA expression was examined in longitudinal PBMC samples obtained from seven HIV/HCV-coinfected, antiretroviral therapy-naïve individuals before, during, and after pegylated interferon-α/ribavirin therapy (IFN-α/RBV). We implemented novel hybrid computational-empirical approaches to characterize regulatory networks between miRNAs and anti-HIV-1 host restriction factors. RESULTS miR-422a was the only miRNA significantly modulated by IFN-α/RBV in vivo (p<0.0001, paired t test; FDR<0.037). Our interactome mapping revealed extensive regulatory involvement of miR-422a in p53-dependent apoptotic and pyroptotic pathways. Based on sequence homology and inverse expression relationships, 29 unique miRNAs may regulate anti-HIV-1 restriction factor expression in vivo. CONCLUSIONS The specific reduction of miR-422a is associated with exogenous IFN-α treatment, and likely contributes to the IFN-α suppression of HIV-1 through the enhancement of anti-HIV-1 restriction factor expression and regulation of genes involved in programmed cell death. Moreover, our regulatory network analysis presents additional candidate miRNAs that may be targeted to enhance anti-HIV-1 restriction factor expression in vivo.