Variation in HIV-1 Nef function within and among viral subtypes reveals genetically separable antagonism of SERINC3 and SERINC5

HIV Nef counteracts cellular host restriction factors SERINC3 and SERINC5, but our understanding of how naturally occurring global Nef sequence diversity impacts these activities is limited. Here, we quantify SERINC3 and SERINC5 internalization function for 339 Nef clones, representing the major pandemic HIV-1 group M subtypes A, B, C and D. We describe distinct subtype-associated hierarchies for Nef-mediated internalization of SERINC5, for which subtype B clones display the highest activities on average, and of SERINC3, for which subtype B clones display the lowest activities on average. We further identify Nef polymorphisms that modulate its ability to counteract SERINC proteins, including substitutions in the N-terminal domain that selectively impair SERINC3 internalization. Our findings demonstrate that the SERINC antagonism activities of HIV Nef differ markedly among major viral subtypes and between individual isolates within a subtype, suggesting that variation in these functions may contribute to global differences in viral pathogenesis.

Suboptimal Biological Sampling as a Probable Cause of False-Negative COVID-19 Diagnostic Test Results

False-negative severe acute respiratory syndrome coronavirus 2 test results can negatively impact the clinical and public health response to coronavirus disease 2019 (COVID-19). We used droplet digital polymerase chain reaction (ddPCR) to demonstrate that human DNA levels, a stable molecular marker of sampling quality, were significantly lower in samples from 40 confirmed or suspected COVID-19 cases that yielded negative diagnostic test results (ie, suspected false-negative test results) compared with a representative pool of 87 specimens submitted for COVID-19 testing. Our results support suboptimal biological sampling as a contributor to false-negative COVID-19 test results and underscore the importance of proper training and technique in the collection of nasopharyngeal specimens.

Differential Vpu-Mediated CD4 and Tetherin Downregulation Functions among Major HIV-1 Group M Subtypes

Downregulation of BST-2/tetherin and CD4 by HIV-1 viral protein U (Vpu) promotes viral egress and allows infected cells to evade host immunity. Little is known however about the natural variability in these Vpu functions among the genetically diverse viral subtypes that contribute to the HIV-1 pandemic. We collected Vpu isolates from 332 treatment-naive individuals living with chronic HIV-1 infection in Uganda, Rwanda, South Africa, and Canada. Together, these Vpu isolates represent four major HIV-1 group M subtypes (A [n = 63], B [n = 84], C [n = 94], and D [n = 59]) plus intersubtype recombinants and uncommon strains (n = 32). The ability of each Vpu clone to downregulate endogenous CD4 and tetherin was quantified using flow cytometry following transfection into an immortalized T-cell line and compared to that of a reference Vpu clone derived from HIV-1 subtype B NL4.3. Overall, the median CD4 downregulation function of natural Vpu isolates was similar to that of NL4.3 (1.01 [interquartile range {IQR}, 0.86 to 1.18]), while the median tetherin downregulation function was moderately lower than that of NL4.3 (0.90 [0.79 to 0.97]). Both Vpu functions varied significantly among HIV-1 subtypes (Kruskal-Wallis P < 0.0001). Specifically, subtype C clones exhibited the lowest CD4 and tetherin downregulation activities, while subtype D and B clones were most functional for both activities. We also identified Vpu polymorphisms associated with CD4 or tetherin downregulation function and validated six of these using site-directed mutagenesis. Our results highlight the marked extent to which Vpu function varies among global HIV-1 strains, raising the possibility that natural variation in this accessory protein may contribute to viral pathogenesis and/or spread.IMPORTANCE The HIV-1 accessory protein Vpu enhances viral spread by downregulating CD4 and BST-2/tetherin on the surface of infected cells. Natural variability in these Vpu functions may contribute to HIV-1 pathogenesis, but this has not been investigated among the diverse viral subtypes that contribute to the HIV-1 pandemic. In this study, we found that Vpu function differs significantly among HIV-1 subtypes A, B, C, and D. On average, subtype C clones displayed the lowest ability to downregulate both CD4 and tetherin, while subtype B and D clones were more functional. We also identified Vpu polymorphisms that associate with functional differences among HIV-1 isolates and subtypes. Our study suggests that genetic diversity in Vpu may play an important role in the differential pathogenesis and/or spread of HIV-1.

Are We Ready for NGS HIV Drug Resistance Testing? The Second "Winnipeg Consensus" Symposium

HIV drug resistance is a major global challenge to successful and sustainable antiretroviral therapy. Next-generation sequencing (NGS)-based HIV drug resistance (HIVDR) assays enable more sensitive and quantitative detection of drug-resistance-associated mutations (DRMs) and outperform Sanger sequencing approaches in detecting lower abundance resistance mutations. While NGS is likely to become the new standard for routine HIVDR testing, many technical and knowledge gaps remain to be resolved before its generalized adoption in regular clinical care, public health, and research. Recognizing this, we conceived and launched an international symposium series on NGS HIVDR, to bring together leading experts in the field to address these issues through in-depth discussions and brainstorming. Following the first symposium in 2018 (Winnipeg, MB Canada, 21-22 February, 2018), a second "Winnipeg Consensus" symposium was held in September 2019 in Winnipeg, Canada, and was focused on external quality assurance strategies for NGS HIVDR assays. In this paper, we summarize this second symposium's goals and highlights.

Longitudinal within-host evolution of HIV Nef-mediated CD4, HLA and SERINC5 downregulation activity: a case study

The HIV accessory protein Nef downregulates the viral entry receptor CD4, the Human Leukocyte Antigen (HLA)-A and -B molecules, the Serine incorporator 5 (SERINC5) protein and other molecules from the infected cell surface, thereby promoting viral infectivity, replication and immune evasion. The nef locus also represents one of the most genetically variable regions in the HIV genome, and nef sequences undergo substantial evolution within a single individual over the course of infection. Few studies however have simultaneously characterized the impact of within-host nef sequence evolution on Nef protein function over prolonged timescales. Here, we isolated 50 unique Nef clones by single-genome amplification over an 11-year period from the plasma of an individual who was largely naïve to antiretroviral treatment during this time. Together, these clones harbored nonsynonymous substitutions at 13% of nef’s codons. We assessed their ability to downregulate cell-surface CD4, HLA and SERINC5 and observed that all three Nef functions declined modestly over time, where the reductions in CD4 and HLA downregulation (an average of 0.6% and 2.0% per year, respectively) achieved statistical significance. The results from this case study support all three Nef activities as being important to maintain throughout untreated HIV infection, but nevertheless suggest that, despite nef’s mutational plasticity, within-host viral evolution can compromise Nef function, albeit modestly, over prolonged periods.

Detection of HIV transmission hotspots in British Columbia, Canada: A novel framework for the prioritization and allocation of treatment and prevention resources

Background

Identifying populations at high risk of HIV transmission is critical for prioritizing treatment and prevention resources and achieving the UNAIDS 90-90-90 Targets.

Methods

HIV transmission rates can be estimated from phylogenetic trees as viral lineage-level diversification rates. To identify HIV-1 transmission foci in British Columbia, Canada, we inferred diversification rates from phylogenetic trees of 36 271 HIV-1 sequences from 9630 anonymized individuals. Diversification rates were combined with sociodemographic and clinical data, then aggregated by patients’ area of residence to predict the distribution of new HIV cases between 2008 and 2018. The predictive power of the model was compared with a phylogenetically uninformed model.

Findings

Aggregated diversification rate measures were predictive of new HIV cases in the subsequent year after adjusting for prevalent and incident cases in the previous year. For every one-unit increase in the mean of the top five diversification rates, the number of new HIV cases increased by on average 1·38-fold (95% CI, 1·28–1·49). In a blind prediction of 2018 cases, diversification rate improved the model's specificity by 12%, accuracy by 9%, top 20 agreement by 100%, and correlation of predicted and observed values by 162% relative to a model that incorporated epidemiological data alone.

Interpretation

By predicting the distribution of future HIV cases, a combined phylogenetic and epidemiological approach identifies hotspots where public health resources are needed most.

Funding

Canadian Institutes of Health Research, University of British Columbia, Public Health Agency of Canada, Genome Canada, Genome BC, Michael Smith Foundation for Health Research, and BC Centre for Excellence in HIV/AIDS.

Rates and Correlates of Short Term Virologic Response among Treatment-Naïve HIV-Infected Children Initiating Antiretroviral Therapy in Ethiopia: A Multi-Center Prospective Cohort Study

There is limited data on virologic outcome and its correlates among HIV-infected children in resource-limited settings. We investigated rate and correlates of virologic outcome among treatment naïve HIV-infected Ethiopian children initiating cART, and were followed prospectively at baseline, 8, 12, 24 and 48 weeks using plasma viral load, clinical examination, laboratory tests and pretreatment HIV drug resistance (PDR) screening. Virologic outcome was assessed using two endpoints–virological suppression defined as having “undetectable” plasma viral load < 150 RNA copies/mL, and rebound defined as viral load ≥150 copies/mL after achieving suppression. Cox Proportional Hazards Regression was employed to assess correlates of outcome. At the end of follow up, virologic outcome was measured for 110 participants. Overall, 94(85.5%) achieved virological suppression, of which 36(38.3%) experienced virologic rebound. At 48 weeks, 9(8.2%) children developed WHO-defined virological treatment failure. Taking tenofovir-containing regimen (Hazard Ratio (HR) 3.1-[95% confidence interval (95%CI) 1.0–9.6], p = 0.049) and absence of pretreatment HIV drug resistance (HR 11.7-[95%CI 1.3–104.2], p = 0.028) were independently associated with earlier virologic suppression. In conclusion, PDR and cART regimen type correlate with rate of virologic suppression which was prominent during the first year of cART initiation. However, the impact of viral rebound in 38.3% of the children needs evaluation.

Frequent cross-resistance to rilpivirine among subtype C HIV-1 from first-line antiretroviral therapy failures in South Africa

Background

Rilpivirine (TMC278LA) is a promising drug for pre-exposure prophylaxis of HIV-1 because of its sub-nanomolar potency and long-acting formulation; however, increasing transmission of non-nucleoside reverse transcriptase inhibitor-resistant HIV-1 with potential cross-resistance to rilpivirine could reduce its preventive efficacy. This study investigated rilpivirine cross-resistance among recombinant subtype C HIV-1 derived from 100 individuals failing on first-line non-nucleoside reverse transcriptase inhibitor-containing antiretroviral therapy in South Africa whose samples were sent for routine HIV-1 drug resistance testing to Lancet Laboratories (Johannesburg, South Africa).

Methods

Plasma samples were selected from individuals with HIV-1 RNA > 10,000 copies/ml and ≥1 non-nucleoside reverse transcriptase inhibitor-resistance mutation in reverse transcriptase. Recombinant HIV-1_LAI-containing bulk-cloned full-length reverse transcriptase sequences from plasma were assayed for susceptibility to nevirapine (NVP), efavirenz (EFV) and rilpivirine in TZM-bl cells. Fold-change (FC) decreases in drug susceptibility were calculated against a mean IC₅₀ from 12 subtype C HIV-1 samples from treatment-naïve individuals in South Africa. Cross-resistance was evaluated based on biological cutoffs established for rilpivirine (2.5-FC) and the effect of mutation combinations on rilpivirine phenotype.

Results

Of the 100 samples from individuals on failing antiretroviral therapy, 69 had 2.5- to 75-fold decreased susceptibility to rilpivirine and 11 had >75-fold resistance. Rilpivirine resistance was strongly associated with K103N especially in combination with other rilpivirine-associated mutations.

Conclusion

The frequently observed cross-resistance of HIV-1 suggests that the preventive efficacy of TMC278LA pre-exposure prophylaxis could be compromised by transmission of HIV-1 from individuals with failure of first-line non-nucleoside reverse transcriptase inhibitor-containing antiretroviral therapy.

A systematic, deep sequencing-based methodology for identification of mixed-genotype hepatitis C virus infections

Hepatitis C virus (HCV) mixed genotype infections can affect treatment outcomes and may have implications for vaccine design and disease progression. Previous studies demonstrate 0-39% of high-risk, HCV-infected individuals harbor mixed genotypes however standardized, sensitive methods of detection are lacking. This study compared PCR amplicon, random primer (RP), and probe enrichment (PE)-based deep sequencing methods coupled with a custom sequence analysis pipeline to detect multiple HCV genotypes. Mixed infection cutoff values, based on HCV read depth and coverage, were identified using receiver operating characteristic curve analysis. The methodology was validated using artificially mixed genotype samples and then applied to two clinical trials of HCV treatment in high-risk individuals (ACTIVATE, 114 samples from 90 individuals; DARE-C II, 26 samples from 18 individuals) and a cohort of HIV/HCV co-infected individuals (Canadian Coinfection Cohort (CCC), 3 samples from 2 individuals with suspected mixed genotype infections). Amplification bias of genotype (G)1b, G2, G3 and G5 was observed in artificially mixed samples using the PCR method while no genotype bias was observed using RP and PE. RP and PE sequencing of 140 ACTIVATE and DARE-C II samples identified the following primary genotypes: 15% (n = 21) G1a, 76% (n = 106) G3, and 9% (n = 13) G2. Sequencing of ACTIVATE and DARE-C II demonstrated, on average, 2% and 1% of HCV reads mapping to a second genotype using RP and PE, respectively, however none passed the mixed infection cutoff criteria and phylogenetics confirmed no mixed infections. From CCC, one mixed infection was confirmed while the other was determined to be a recombinant genotype. This study underlines the risk for false identification of mixed HCV infections and stresses the need for standardized methods to improve prevalence estimates and to understand the impact of mixed infections for management and elimination of HCV.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus

INTRODUCTION

Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting.

DISCUSSION

Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use.

CONCLUSIONS

While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

Novel HLA class I associations with HIV-1 control in a unique genetically admixed population

Associations between HLA class I alleles and HIV progression in populations exhibiting Amerindian and Caucasian genetic admixture remain understudied. Using univariable and multivariable analyses we evaluated HLA associations with five HIV clinical parameters in 3,213 HIV clade B-infected, ART-naïve individuals from Mexico and Central America (MEX/CAM cohort). A Canadian cohort (HOMER, n = 1622) was used for comparison. As expected, HLA allele frequencies in MEX/CAM and HOMER differed markedly. In MEX/CAM, 13 HLA-A, 24 HLA-B, and 14 HLA-C alleles were significantly associated with at least one clinical parameter. These included previously described protective (e.g. B*27:05, B*57:01/02/03 and B*58:01) and risk (e.g. B*35:02) alleles, as well as novel ones (e.g. A*03:01, B*15:39 and B*39:02 identified as protective, and A*68:03/05, B*15:30, B*35:12/14, B*39:01/06, B*39:05~C*07:02, and B*40:01~C*03:04 identified as risk). Interestingly, both protective (e.g. B*39:02) and risk (e.g. B*39:01/05/06) subtypes were identified within the common and genetically diverse HLA-B*39 allele group, characteristic to Amerindian populations. While HLA-HIV associations identified in MEX and CAM separately were similar overall (Spearman's rho = 0.33, p = 0.03), region-specific associations were also noted. The identification of both canonical and novel HLA/HIV associations provides a first step towards improved understanding of HIV immune control among unique and understudied Mestizo populations.

High Levels of Dual-Class Drug Resistance in HIV-Infected Children Failing First-Line Antiretroviral Therapy in Southern Ethiopia

Clinical monitoring of pediatric HIV treatment remains a major challenge in settings where drug resistance genotyping is not routinely available. As a result, our understanding of drug resistance, and its impact on subsequent therapeutic regimens available in these settings, remains limited. We investigate the prevalence and correlates of HIV-1 drug resistance among 94 participants of the Ethiopia Pediatric HIV Cohort failing first-line combination antiretroviral therapy (cART) using dried blood spot-based genotyping. Overall, 81% (73/90) of successfully genotyped participants harbored resistance mutations, including 69% (62/90) who harbored resistance to both Nucleoside Reverse Transcriptase Inhibitors (NRTIs) and Non-nucleoside Reverse Transcriptase Inhibitors (NNRTIs). Strikingly, 42% of resistant participants harbored resistance to all four NRTIs recommended for second-line use in this setting, meaning that there are effectively no remaining cART options for these children. Longer cART duration and prior regimen changes were significantly associated with detection of drug resistance mutations. Replicate genotyping increased the breadth of drug resistance detected in 34% of cases, and thus is recommended for consideration when typing from blood spots. Implementation of timely drug resistance testing and access to newer antiretrovirals and drug classes are urgently needed to guide clinical decision-making and improve outcomes for HIV-infected children on first-line cART in Ethiopia.

Sociodemographic correlates of HIV drug resistance and access to drug resistance testing in British Columbia, Canada

Sociodemographic correlates of engagement in human immunodeficiency virus (HIV) care are well studied, however the association with accessing drug resistance testing (DRT) and the development of drug resistance have not been characterized. Between 1996-2014, 11 801 HIV patients accessing therapy in British Columbia were observed longitudinally. A subset of 9456 patients had testable viral load; of these 8398 were linked to census data. Sociodemographic (census tract-level) and clinical (individual-level) correlates of DRT were assessed using multivariable General Estimating Equation logistic regression adjusted odds ratios (aOR). The mean number of tests per patient was 2.1 (Q1-Q3; 0-3). Separately, any drug resistance was determined using IAS-USA (2013) list for 5703 initially treatment naïve patients without baseline resistance; 5175 were census-linked (mean of 1.5 protease-reverse transcriptase sequences/patient, Q1-Q3; 0-2). Correlates of detecting drug resistance in this subset were analyzed using Cox PH regression adjusted hazard ratios (aHR). Our results indicate baseline CD4 <200 cells/μL (aOR: 1.5, 1.3-1.6), nRTI-only baseline regimens (aOR: 1.4, 1.3-1.6), and unknown (therapy initiation before routine pVL in BC) baseline pVL (aOR: 1.8, 1.5-2.1) were among individual-level clinical covariates strongly associated with having accessed DRT; while imperfect adherence (aHR: 2.2, 1.9-2.5), low baseline CD4 count (aHR: 1.9, 1.6-2.3), and high baseline pVL (aHR: 2.0, 1.6-2.6) were associated with a higher likelihood of developing drug resistance. A higher median income (aOR: 0.83, 0.77-0.89) and higher percentage of those with aboriginal ancestry (aOR: 0.85, 0.76-0.95) were census tract-level sociodemographic covariates associated with decreased access to DRT. Similarly, aboriginal ancestry (aHR: 1.2, 1.1-1.5) was associated with development of drug resistance. In conclusion, clinical covariates continue to be the strongest correlates of development of drug resistance and access to DRT for individuals. Regions of high median income and high aboriginal ancestry were weak census-level sociodemographic indicators of reduced DRT uptake, however high aboriginal ancestry was the only sociodemographic indicator for development of drug resistance.

Longitudinal trends of HIV drug resistance in a large Canadian cohort, 1996-2016

OBJECTIVES

We aim to identify long-term trends in HIV drug resistance before and after combined antiretroviral therapy (cART) initiation.

METHODS

IAS-USA (2015) mutations were identified in 23 271 HIV protease-reverse transcriptase sequences from 6543 treatment naïve adults in British Columbia. Participants who started cART between 1996 and 2014 were followed until April 2016. Equality of proportions test was used to compare the percentage of participants with acquired drug resistance (ADR) or transmitted drug resistance (TDR) in 1996, to those in 2014. Kaplan-Meier was used to estimate time to ADR in four drug resistance categories. Multivariable regression odds ratios (OR) of ADR for select clinical variables were determined by 5-year eras of cART initiation.

RESULTS

The proportion of individuals with ADR declined from 39% (51/132) to 3% (8/322) in 1996-2014 (p <0.0001), while the proportion with TDR increased from 12% (16/132) to 18% (59/322) (p 0.14). The estimated proportions of individuals with ADR rose to 29% (NNRTI), 28% (3TC/FTC), 14% (other nRTI), and 7% (PI) after >16 years of therapy. After 5 years on therapy, participants initiating cART in 1996-2000 had 5.5-times more 3TC/FTC ADR, 5.3-times more other nRTI ADR, 4.7-times more NNRTI ADR, and 24-times more PI ADR than those starting in 2011-2014. The individuals with highest odds of developing ADR in 1996-2010 were adherent to regimens at levels between 60% and 80%, which shifted to <40% adherent in 2011-2014.

CONCLUSIONS

HIV drug resistance transitioned from being primarily selected de-novo to being driven by TDR. Among those who started treatment in the past 5 years, ADR is rare and observed mostly in the lowest adherence strata.

Genotypic susceptibility score (GSS) and CD4+ T cell recovery in HIV-1 patients with suppressed viral load

OBJECTIVES

HIV drug resistance, measured by the genotypic susceptibility score (GSS), has a deleterious effect on the virological outcome of HIV-1-infected patients. However, it is not known if GSS retains any predictive value for CD4 recovery in patients with suppressed viral load.

METHODS

Four hundred and six patients on virological failure (>500 copies/mL) with GSS  : <6 months prior to switch therapy who achieved undetectable plasma viral load (<50 copies/mL) within 1 year, remained undetectable >1 year on an unchanged regimen and had CD4 data available during entire follow-up were included. Adjusted and unadjusted analyses of all characteristics at switch related to CD4 recovery were made for three time frames: (i) 'switch-suppression'; (ii) 'suppression-1 year'; and (iii) 'switch-1 year'.

RESULTS

Higher GSS was associated with a greater CD4 recovery between 'switch' and '1 year' in the unadjusted analysis (P = 0.010); however, the effect of GSS was no longer statistically significant after adjusting for pre-switch clinical (CD4 count and plasma viral load) and demographic variables. Furthermore, only a lower pre-switch CD4 count was associated with increased CD4 recovery in the 'suppression-1 year' period in both unadjusted and adjusted models. The main CD4 recovery occurred in 'switch-suppression' and the variables associated, both unadjusted and adjusted, were CD4 and plasma viral load at switch, maintaining a trend for GSS (P = 0.06).

CONCLUSIONS

In individuals who re-suppressed HIV viraemia after switching therapy, regimens having a higher GSS were associated with improved CD4 recovery only during the period from switch to virological suppression, but, once viral load is re-suppressed, the GSS of the new regimen has no further effect on subsequent CD4 recovery.

Promises and pitfalls of Illumina sequencing for HIV resistance genotyping

Genetic sequencing ("genotyping") plays a critical role in the modern clinical management of HIV infection. This virus evolves rapidly within patients because of its error-prone reverse transcriptase and short generation time. Consequently, HIV variants with mutations that confer resistance to one or more antiretroviral drugs can emerge during sub-optimal treatment. There are now multiple HIV drug resistance interpretation algorithms that take the region of the HIV genome encoding the major drug targets as inputs; expert use of these algorithms can significantly improve to clinical outcomes in HIV treatment. Next-generation sequencing has the potential to revolutionize HIV resistance genotyping by lowering the threshold that rare but clinically significant HIV variants can be detected reproducibly, and by conferring improved cost-effectiveness in high-throughput scenarios. In this review, we discuss the relative merits and challenges of deploying the Illumina MiSeq instrument for clinical HIV genotyping.

Lack of a significant impact of Gag-Protease-mediated HIV-1 replication capacity on clinical parameters in treatment-naive Japanese individuals

Background

HLA class I-associated escape mutations in HIV-1 Gag can reduce viral replication, suggesting that associated fitness costs could impact HIV-1 disease progression. Previous studies in North American and African cohorts have reported reduced Gag-Protease mediated viral replication capacity (Gag-Pro RC) in individuals expressing protective HLA class I alleles including HLA-B*57:01, B*27:05, and B*81:01. These studies also reported significant positive associations between Gag-Pro RCs and plasma viral load (pVL). However, these HLA alleles are virtually absent in Japan, and the importance of Gag as an immune target is not clearly defined in this population.

Results

We generated chimeric NL4-3 viruses carrying patient-derived Gag-Protease from 306 treatment-naive Japanese individuals chronically infected with HIV-1 subtype B. We analyzed associations between Gag-Pro RC and clinical markers of HIV-1 infection and host HLA expression. We observed no significant correlation between Gag-Pro RC and pVL in Japan in the overall cohort. However, upon exclusion of individuals expressing Japanese protective alleles HLA-B*52:01 and B*67:01, Gag-Pro RC correlated positively with pVL and negatively with CD4 T-cell count. Our results thus contrast with studies from other global cohorts reporting significantly lower Gag-Pro RC among persons expressing protective HLA alleles, and positive relationships between Gag-Pro RC and pVL in the overall study populations. We also identified five amino acids in Gag-Protease significantly associated with Gag-Pro RC, whose effects on RC were confirmed by site-directed mutagenesis. However, of the four mutations that decreased Gag-Pro RC, none were associated with reductions in pVL in Japan though two were associated with lower pVL in North America.

Conclusions

These data indicate that Gag fitness does not affect clinical outcomes in subjects with protective HLA class I alleles as well as the whole Japanese population. Moreover, the impact of Gag fitness costs on HIV-1 clinical parameters in chronic infection is likely low in Japan compared to other global populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0223-z) contains supplementary material, which is available to authorized users.

The impact of clinical, demographic and risk factors on rates of HIV transmission: a population-based phylogenetic analysis in British Columbia, Canada

BACKGROUND

The diversification of human immunodeficiency virus (HIV) is shaped by its transmission history. We therefore used a population based province wide HIV drug resistance database in British Columbia (BC), Canada, to evaluate the impact of clinical, demographic, and behavioral factors on rates of HIV transmission.

METHODS

We reconstructed molecular phylogenies from 27,296 anonymized bulk HIV pol sequences representing 7747 individuals in BC-about half the estimated HIV prevalence in BC. Infections were grouped into clusters based on phylogenetic distances, as a proxy for variation in transmission rates. Rates of cluster expansion were reconstructed from estimated dates of HIV seroconversion.

RESULTS

Our criteria grouped 4431 individuals into 744 clusters largely separated with respect to risk factors, including large established clusters predominated by injection drug users and more-recently emerging clusters comprising men who have sex with men. The mean log10 viral load of an individual's phylogenetic neighborhood (composed of 5 other individuals with shortest phylogenetic distances) increased their odds of appearing in a cluster by >2-fold per log10 viruses per milliliter.

CONCLUSIONS

Hotspots of ongoing HIV transmission can be characterized in near real time by the secondary analysis of HIV resistance genotypes, providing an important potential resource for targeting public health initiatives for HIV prevention.

Prevalence and virologic consequences of transmitted HIV-1 drug resistance in Uganda

Few reports have examined the impact of HIV-1 transmitted drug resistance (TDR) in resource-limited settings where there are fewer regimen choices and limited pretherapy/posttherapy resistance testing. In this study, we examined TDR prevalence in Kampala and Mbarara, Uganda and assessed its virologic consequences after antiretroviral therapy initiation. We sequenced the HIV-1 protease/reverse transcriptase from n=81 and n=491 treatment-naive participants of the Uganda AIDS Rural Treatment Outcomes (UARTO) pilot study in Kampala (AMU 2002-2004) and main cohort in Mbarara (MBA 2005-2010). TDR-associated mutations were defined by the WHO 2009 surveillance mutation list. Posttreatment viral load data were available for both populations. Overall TDR prevalence was 7% (Kampala) and 3% (Mbarara) with no significant time trend. There was a slight but statistically nonsignificant trend indicating that the presence of TDR was associated with a worse treatment outcome. Virologic suppression (≤400 copies/ml within 6 months posttherapy initiation) was achieved in 87% and 96% of participants with wildtype viruses versus 67% and 83% of participants with TDR (AMU, MBA p=0.2 and 0.1); time to suppression (log-rank p=0.3 and p=0.05). Overall, 85% and 96% of study participants achieved suppression regardless of TDR status. Surprisingly, among the TDR cases, approximately half still achieved suppression; the presence of pretherapy K103N while on nevirapine and fewer active drugs in the first regimen were most often observed with failures. The majority of patients benefited from the local HIV care system even without resistance monitoring. Overall, TDR prevalence was relatively low and its presence did not always imply treatment failure.

Theoretical and experimental assessment of degenerate primer tagging in ultra-deep applications of next-generation sequencing

Primer IDs (pIDs) are random oligonucleotide tags used in next-generation sequencing to identify sequences that originate from the same template. These tags are produced by degenerate primers during the reverse transcription of RNA molecules into cDNA. The use of pIDs helps to track the number of RNA molecules carried through amplification and sequencing, and allows resolution of inconsistencies between reads sharing a pID. Three potential issues complicate the above applications. First, multiple cDNAs may share a pID by chance; we found that while preventing any cDNAs from sharing a pID may be unfeasible, it is still practical to limit the number of these collisions. Secondly, a pID must be observed in at least three sequences to allow error correction; as such, pIDs observed only one or two times must be rejected. If the sequencing product contains copies from a high number of RT templates but produces few reads, our findings indicate that rejecting such pIDs will discard a great deal of data. Thirdly, the use of pIDs could influence amplification and sequencing. We examined the effects of several intrinsic and extrinsic factors on sequencing reads at both the individual and ensemble level.

Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic

HLA-restricted immune escape mutations that persist following HIV transmission could gradually spread through the viral population, thereby compromising host antiviral immunity as the epidemic progresses. To assess the extent and phenotypic impact of this phenomenon in an immunogenetically diverse population, we genotypically and functionally compared linked HLA and HIV (Gag/Nef) sequences from 358 historic (1979–1989) and 382 modern (2000–2011) specimens from four key cities in the North American epidemic (New York, Boston, San Francisco, Vancouver). Inferred HIV phylogenies were star-like, with approximately two-fold greater mean pairwise distances in modern versus historic sequences. The reconstructed epidemic ancestral (founder) HIV sequence was essentially identical to the North American subtype B consensus. Consistent with gradual diversification of a “consensus-like” founder virus, the median “background” frequencies of individual HLA-associated polymorphisms in HIV (in individuals lacking the restricting HLA[s]) were ∼2-fold higher in modern versus historic HIV sequences, though these remained notably low overall (e.g. in Gag, medians were 3.7% in the 2000s versus 2.0% in the 1980s). HIV polymorphisms exhibiting the greatest relative spread were those restricted by protective HLAs. Despite these increases, when HIV sequences were analyzed as a whole, their total average burden of polymorphisms that were “pre-adapted” to the average host HLA profile was only ∼2% greater in modern versus historic eras. Furthermore, HLA-associated polymorphisms identified in historic HIV sequences were consistent with those detectable today, with none identified that could explain the few HIV codons where the inferred epidemic ancestor differed from the modern consensus. Results are therefore consistent with slow HIV adaptation to HLA, but at a rate unlikely to yield imminent negative implications for cellular immunity, at least in North America. Intriguingly, temporal changes in protein activity of patient-derived Nef (though not Gag) sequences were observed, suggesting functional implications of population-level HIV evolution on certain viral proteins.

Upon HIV transmission, many – though not all – immune escape mutations selected in the previous host will revert to the consensus residue. The persistence of certain escape mutations following transmission has led to concerns that these could gradually accumulate in circulating HIV sequences over time, thereby undermining host antiviral immune potential as the epidemic progresses. As certain immune-driven mutations reduce viral fitness, their spread through the population could also have consequences for the average replication capacity and/or protein function of circulating HIV sequences. Here, we characterized HIV sequences, linked to host immunogenetic information, from patients enrolled in historic (1979–1989) and modern (2000–2011) HIV cohorts from four key cities in the North American epidemic. We reconstructed the epidemic's ancestral (founder) HIV sequence and assessed the subsequent extent to which known HIV immune escape mutations have spread in the population. Our data support the gradual spread of many - though not all - immune escape mutations in HIV sequences over time, but to an extent that is unlikely to have major immediate immunologic consequences for the North American epidemic. Notably, in vitro assessments of ancestral and patient-derived HIV sequences suggested functional implications of ongoing HIV evolution for certain viral proteins.

Structure of TCR and antigen complexes at an immunodominant CTL epitope in HIV-1 infection

We investigated the crystal structure of an HLA-A*2402-restricted CTL epitope in the HIV-1 nef gene (Nef134-10) before (pHLA) or after TCR docking. The wild type epitope and two escape mutants were included in the study. Y135F was an early-appearing major mutation, while F139L was a late-appearing mutation which was selected in the patients without Y135F. F139 was an eminent feature of the Nef134-10 epitope. Wild type-specific TCR was less fit to F139L mutant suggesting that F139L is an escape from the CTL against the wild type epitope. Although Y135F mutation disrupted the hydrogen bond to HLA-A*2402 His70, newly formed hydrogen bond between T138 and His70 kept the conformation of the epitope in the reconstituted pMHC. TCR from Y135F- or dually-specific CTL had unique mode of binding to the mutant epitope. Y135F has been reported as a processing mutant but CTL carrying structurally adequate TCR can be found in the patients.