Real world SOF/VEL/VOX retreatment outcomes and viral resistance analysis for HCV patients with prior failure to DAA therapy

Sustained viral response (SVR) rates for direct-acting antiviral (DAA) therapy for hepatitis C virus (HCV) infection routinely exceed 95%. However, a small number of patients require retreatment. Sofosbuvir, velpatasvir and voxilaprevir (SOF/VEL/VOX) is a potent DAA combination primarily used for the retreatment of patients who failed by DAA therapies. Here we evaluate retreatment outcomes and the effects of resistance-associated substitutions (RAS) in a real-world cohort, including a large number of genotype (GT)3 infected patients. 144 patients from the UK were retreated with SOF/VEL/VOX following virologic failure with first-line DAA treatment regimens. Full-length HCV genome sequencing was performed prior to retreatment with SOF/VEL/VOX. HCV subtypes were assigned and RAS relevant to each genotype were identified. GT1a and GT3a each made up 38% (GT1a n = 55, GT3a n = 54) of the cohort. 40% (n = 58) of patients had liver cirrhosis of whom 7% (n = 4) were decompensated, 10% (n = 14) had hepatocellular carcinoma (HCC) and 8% (n = 12) had received a liver transplant prior to retreatment. The overall retreatment SVR12 rate was 90% (129/144). On univariate analysis, GT3 infection (50/62; SVR = 81%, p = .009), cirrhosis (47/58; SVR = 81%, p = .01) and prior treatment with SOF/VEL (12/17; SVR = 71%, p = .02) or SOF+DCV (14/19; SVR = 74%, p = .012) were significantly associated with retreatment failure, but existence of pre-retreatment RAS was not when viral genotype was taken into account. Retreatment with SOF/VEL/VOX is very successful for non-GT3-infected patients. However, for GT3-infected patients, particularly those with cirrhosis and failed by initial SOF/VEL treatment, SVR rates were significantly lower and alternative retreatment regimens should be considered.

Coevolved Multidrug-Resistant HIV-1 Protease and Reverse Transcriptase Influences Integrase Drug Susceptibility and Replication Fitness

Integrase strand transfer inhibitors (InSTIs) are recommended agents in first-line combination antiretroviral therapy (cART). We examined the evolution of drug resistance mutations throughout HIV-1 pol and the effects on InSTI susceptibility and viral fitness. We performed single-genome sequencing of full-length HIV-1 pol in a highly treatment-experienced patient, and determined drug susceptibility of patient-derived HIV-1 genomes using a phenotypic assay encompassing full-length pol gene. We show the genetic linkage of multiple InSTI-resistant haplotypes containing major resistance mutations at Y143, Q148 and N155 to protease inhibitor (PI) and reverse transcriptase inhibitor (RTI) resistance mutations. Phenotypic analysis of viruses expressing patient-derived IN genes with eight different InSTI-resistant haplotypes alone or in combination with coevolved protease (PR) and RT genes exhibited similar levels of InSTI susceptibility, except for three haplotypes that showed up to 3-fold increases in InSTI susceptibility (p ≤ 0.032). The replicative fitness of most viruses expressing patient-derived IN only significantly decreased, ranging from 8% to 56% (p ≤ 0.01). Interestingly, the addition of coevolved PR + RT significantly increased the replicative fitness of some haplotypes by up to 73% (p ≤ 0.024). Coevolved PR + RT contributes to the susceptibility and viral fitness of patient-derived IN viruses. Maintaining patients on failing cART promotes the selection of fitter resistant strains, and thereby limits future therapy options.

Clinical evaluation of a Hepatitis C Virus whole-genome sequencing pipeline for genotyping and resistance testing

OBJECTIVES

We sought to evaluate clinically a hepatitis C virus (HCV) whole-genome, next-generation sequencing (NGS) pipeline that is agnostic to viral genotype.

METHODS

Performance of the NGS pipeline was assessed through comparison of results with Sanger sequencing (SS) of partial HCV genomes.

RESULTS

There was 98.7% (376/381) concordance for viral subtype between SS and NGS. The positive and negative per cent agreements for determination of resistance-associated substitutions were 97.8% (95% CI 92.5-99.4%) and 99.9% (95% CI 99.5-100.0%), respectively. The NGS pipeline was also able to detect novel subtypes, mixtures, recombinants, transiently occurring resistance mutations and distinguish re-infection with the same subtype from relapse.

DISCUSSION

Particular scenarios where NGS may be used include settings without universal access to pan-genotypic antiviral regimens, those infected with a 'rare' subtype or who have been failed by first-line therapy, and in cases of suspected re-infection.

Surveillance of HIV-1 transmitted integrase strand transfer inhibitor resistance in the UK

Background

HIV treatment guidelines have traditionally recommended that all HIV-positive individuals are tested for evidence of drug resistance prior to starting ART. Testing for resistance to reverse transcriptase inhibitors and PIs is well established in routine care. However, testing for integrase strand transfer inhibitor (InSTI) resistance is less consistent.

Objectives

To inform treatment guidelines by determining the prevalence of InSTI resistance in a national cohort of recently infected individuals.

Patients and methods

Recent (within 4 months) HIV-1 infections were identified using a Recent Infection Testing Algorithm of new HIV-1 diagnoses in the UK. Resistance-associated mutations (RAMs) in integrase, protease and reverse transcriptase were detected by ultradeep sequencing, which allows for the sensitive estimation of the frequency of each resistant variant in a sample.

Results

The analysis included 655 randomly selected individuals (median age = 33 years, 95% male, 83% MSM, 78% white) sampled in the period 2014 to 2016 and determined to have a recent infection. These comprised 320, 138 and 197 samples from 2014, 2015 and 2016, respectively. None of the samples had major InSTI RAMs occurring at high variant frequency (≥20%). A subset (25/640, 3.9%) had major InSTI RAMs occurring only as low-frequency variants (2%–20%). In contrast, 47/588 (8.0%) had major reverse transcriptase inhibitor and PI RAMs at high frequency.

Conclusions

Between 2014 and 2016, major InSTI RAMs were uncommon in adults with recent HIV-1 infection, only occurring as low-frequency variants of doubtful clinical significance. Continued surveillance of newly diagnosed patients for evidence of transmitted InSTI resistance is recommended to inform clinical practice.

Technical Validation of a Hepatitis C Virus Whole Genome Sequencing Assay for Detection of Genotype and Antiviral Resistance in the Clinical Pathway

Choice of direct acting antiviral (DAA) therapy for Hepatitis C Virus (HCV) in the United Kingdom and similar settings usually requires knowledge of the genotype and, in some cases, antiviral resistance (AVR) profile of the infecting virus. To determine these, most laboratories currently use Sanger technology, but next-generation sequencing (NGS) offers potential advantages in throughput and accuracy. However, NGS poses unique technical challenges, which require idiosyncratic development and technical validation approaches. This applies particularly to virology, where sequence diversity is high and the amount of starting genetic material is low, making it difficult to distinguish real data from artifacts. We describe the development and technical validation of a sequence capture-based HCV whole genome sequencing (WGS) assay to determine viral genotype and AVR profile. We use clinical samples of known subtypes and viral loads, and simulated FASTQ datasets to validate the analytical performances of both the wet laboratory and bioinformatic pipeline procedures. We show high concordance of the WGS assay compared to current "gold standard" Sanger assays. Specificity was 92.3 and 96.1% for AVR and genotyping, respectively. Discordances were due to the inability of Sanger assays to assign the correct subtype or accurately call mixed drug-resistant variants. We show high repeatability and reproducibility with >99.8% sequence similarity between sequence runs as well as high precision for variant frequency detection at >98.8% in the 95th percentile. Post-sequencing bioinformatics quality control workflows allow the accurate distinction between mixed infections, cross-contaminants and recombinant viruses at a threshold of >5% for the minority population. The sequence capture-based HCV WGS assay is more accurate than legacy AVR and genotyping assays. The assay has now been implemented in the clinical pathway of England's National Health Service HCV treatment programs, representing the first validated HCV WGS pipeline in clinical service. The data generated will additionally provide granular national-level genomic information for public health policy making and support the WHO HCV elimination strategy.

Determining the Origins of Human Immunodeficiency Virus Type 1 Drug-resistant Minority Variants in People Who Are Recently Infected Using Phylogenetic Reconstruction

BACKGROUND

Drug-resistant minority variants (DRMinVs) detected in patients who recently acquired human immunodeficiency virus type 1 (HIV-1) can be transmitted, generated de novo through virus replication, or technical errors. The first form is likely to persist and result in treatment failure, while the latter two could be stochastic and transient.

METHODS

Ultradeep sequencing of plasma samples from 835 individuals with recent HIV-1 infection in the United Kingdom was performed to detect DRMinVs at a mutation frequency between 2% and 20%. Sequence alignments including >110 000 HIV-1 partial pol consensus sequences from the UK HIV Drug Resistance Database (UK-HDRD), linked to epidemiological and clinical data from the HIV and AIDS Reporting System, were used for transmission cluster analysis. Transmission clusters were identified using Cluster Picker with a clade support of >90% and maximum genetic distances of 4.5% or 1.5%, the latter to limit detection to likely direct transmission events.

RESULTS

Drug-resistant majority variants (DRMajVs) were detected in 66 (7.9%) and DRMinVs in 84 (10.1%) of the recently infected individuals. High levels of clustering to sequences in UK-HDRD were observed for both DRMajV (n = 48; 72.7%) and DRMinV (n = 63; 75.0%) sequences. Of these, 43 (65.2%) with DRMajVs were in a transmission cluster with sequences that harbored the same DR mutation compared to only 3 (3.6%) sequences with DRMinVs (P < .00001, Fisher exact test). Evidence of likely direct transmission of DRMajVs was observed for 25/66 (37.9%), whereas none were observed for the DRMinVs (P < .00001).

CONCLUSIONS

Using a densely sampled HIV-infected population, we show no evidence of DRMinV transmission among recently infected individuals.

Enhanced Detection of DNA Viruses in the Cerebrospinal Fluid of Encephalitis Patients Using Metagenomic Next-Generation Sequencing

The long and expanding list of viral pathogens associated with causing encephalitis confounds current diagnostic procedures, and in up to 50% of cases, the etiology remains undetermined. Sequence-agnostic metagenomic next-generation sequencing (mNGS) obviates the need to specify targets in advance and thus has great potential in encephalitis diagnostics. However, the low relative abundance of viral nucleic acids in clinical specimens poses a significant challenge. Our protocol employs two novel techniques to selectively remove human material at two stages, significantly increasing the representation of viral material. Our bioinformatic workflow using open source protein- and nucleotide sequence-matching software balances sensitivity and specificity in diagnosing and characterizing any DNA viruses present. A panel of 12 cerebrospinal fluid (CSFs) from encephalitis cases was retrospectively interrogated by mNGS, with concordant results in seven of nine samples with a definitive DNA virus diagnosis, and a different herpesvirus was identified in the other two. In two samples with an inconclusive diagnosis, DNA viruses were detected and in a virus-negative sample, no viruses were detected. This assay has the potential to detect DNA virus infections in cases of encephalitis of unknown etiology and to improve the current screening tests by identifying new and emerging agents.

Consensus recommendations for resistance testing in the management of chronic hepatitis C virus infection: Public Health England HCV Resistance Group

The treatment of hepatitis C virus (HCV) infection has been revolutionised by the advent of oral, well-tolerated, direct acting antiviral therapies (DAA), with high cure rates. However, in some scenarios, HCV resistance to antiviral therapies may have an impact on treatment success. Public Health England's HCV Resistance Group was established to support clinicians treating people with HCV, where the issue of resistance may be a factor in clinical decision-making, and this review includes the Group's current recommendations on the use of HCV resistance testing. The authors describe the principles behind and approach to HCV resistance testing and consider evidence from in vitro studies, clinical trials and real world cohorts on the impact of HCV resistance on treatment outcomes for particular DAA regimens. Five scenarios are identified in the UK and similar settings, where, in the Group's opinion, resistance testing should be performed.

Potential for diagnosis of infectious disease from the 100,000 Genomes Project Metagenomic Dataset: Recommendations for reporting results

The identification of microbiological infection is usually a diagnostic investigation, a complex process that is firstly initiated by clinical suspicion. With the emergence of high-throughput sequencing (HTS) technologies, metagenomic analysis has unveiled the power to identify microbial DNA/RNA from a diverse range of clinical samples (1). Metagenomic analysis of whole human genomes at the clinical/research interface bypasses the steps of clinical scrutiny and targeted testing and has the potential to generate unexpected findings relating to infectious and sometimes transmissible disease. There is no doubt that microbial findings that may have a significant impact on a patient’s treatment and their close contacts should be reported to those with clinical responsibility for the sample-donating patient. There are no clear recommendations on how such findings that are incidental, or outside the original investigation, should be handled. Here we aim to provide an informed protocol for the management of incidental microbial findings as part of the 100,000 Genomes Project which may have broader application in this emerging field. As with any other clinical information, we aim to prioritise the reporting of data that are most likely to be of benefit to the patient and their close contacts. We also set out to minimize risks, costs and potential anxiety associated with the reporting of results that are unlikely to be of clinical significance. Our recommendations aim to support the practice of microbial metagenomics by providing a simplified pathway that can be applied to reporting the identification of potential pathogens from metagenomic datasets. Given that the ambition for UK sequenced human genomes over the next 5 years has been set to reach 5 million and the field of metagenomics is rapidly evolving, the guidance will be regularly reviewed and will likely adapt over time as experience develops.

Presence, persistence and effects of pre-treatment HIV-1 drug resistance variants detected using next generation sequencing: A Retrospective longitudinal study from rural coastal Kenya

BACKGROUND

The epidemiology of HIV-1 drug resistance (HIVDR) determined by Sanger capillary sequencing, has been widely studied. However, much less is known about HIVDR detected using next generation sequencing (NGS) methods. We aimed to determine the presence, persistence and effect of pre-treatment HIVDR variants detected using NGS in HIV-1 infected antiretroviral treatment (ART) naïve participants from rural Coastal Kenya.

METHODS

In a retrospective longitudinal study, samples from HIV-1 infected participants collected prior [n = 2 time-points] and after [n = 1 time-point] ART initiation were considered. An ultra-deep amplicon-based NGS assay, calling for nucleotide variants at >2.0% frequency of viral population, was used. Suspected virologic failure (sVF) was defined as a one-off HIV-1 viral load of >1000 copies/ml whilst on ART.

RESULTS

Of the 50 eligible participants, 12 (24.0% [95% CI: 13.1-38.2]) had at least one detectable pre-treatment HIVDR variant against Protease Inhibitors (PIs, n = 6 [12%]), Nucleoside Reverse Transcriptase Inhibitors (NRTIs, n = 4 [8.0%]) and Non-NRTIs (n = 3 [6.0%]). Overall, 15 pre-treatment resistance variants were detected (frequency, range: 2.3-92.0%). A positive correlation was observed between mutation frequency and absolute load for NRTI and/or NNRTI variants (r = 0.761 [p = 0.028]), but not for PI variants (r = -0.117 [p = 0.803]). Participants with pre-treatment NRTI and/or NNRTI resistance had increased odds of sVF (OR = 6.0; 95% CI = 1.0-36.9; p = 0.054).

CONCLUSIONS

Using NGS, pre-treatment resistance variants were common, though observed PI variants were unlikely transmitted, but rather probably generated de novo. Even when detected from a low frequency, pre-treatment NRTI and/or NNRTI resistance variants may adversely affect treatment outcomes.

An Innovative Study Design to Assess the Community Effect of Interventions to Mitigate HIV Epidemics Using Transmission-Chain Phylodynamics

Given globalization and other social phenomena, controlling the spread of infectious diseases has become an imperative public health priority. A plethora of interventions that in theory can mitigate the spread of pathogens have been proposed and applied. Evaluating the effectiveness of such interventions is costly and in many circumstances unrealistic. Most important, the community effect (i.e., the ability of the intervention to minimize the spread of the pathogen from people who received the intervention to other community members) can rarely be evaluated. Here we propose a study design that can build and evaluate evidence in support of the community effect of an intervention. The approach exploits molecular evolutionary dynamics of pathogens in order to track new infections as having arisen from either a control or an intervention group. It enables us to evaluate whether an intervention reduces the number and length of new transmission chains in comparison with a control condition, and thus lets us estimate the relative decrease in new infections in the community due to the intervention. We provide as an example one working scenario of a way the approach can be applied with a simulation study and associated power calculations.

Impact of Interferon-α Receptor-1 Promoter Polymorphisms on the Transcriptome of the Hepatitis B Virus-Associated Hepatocellular Carcinoma

Background and aims

Genetic polymorphisms within the promoter of interferon-α receptor type-1 (IFNAR1) have been associated with the susceptibility to and the outcome of chronic hepatitis B virus (HBV) infection. However, the impact of these polymorphisms in the transcriptome of the HBV-associated hepatocellular carcinoma (HCC) remains largely unexplored.

Methods

Using whole-genome and exome sequencing data from The Cancer Genome Atlas project, we characterized three single-nucleotide polymorphisms (SNPs: −568G/C, −408C/T, −3C/T) and one variable number tandem repeat [VNTR: −77(GT)n] within the IFNAR1 promoter sequence in 49 HCC patients. RNAseq data from 10 genotyped HCC samples were grouped according to their −77VNTR or −3SNP genotype to evaluate the impact of these polymorphisms on the differential expression on the HCC transcriptome.

Results

There is a fourfold higher impact of the −77VNTR on the HCC transcriptome compared to the −3SNP (q < 0.1, p < 0.001). The expression of the primary IFNAR1 transcript is not affected by these polymorphisms but a secondary, HCC-specific transcript is expressed only in homozygous −77VNTR ≤8/≤8(GT)n samples (p < 0.05). At the same time, patients carrying at least one −77VNTR >8(GT) allele, presented a strong upregulation of the fibronectin-1 (FN-1) gene, which has been associated with the development of HCC. Gene Ontology and pathway enrichment analysis of the differentially expressed genes revealed a strong disruption of the PI3K–AKT signaling pathway, which can be partially triggered by the extracellular matrix FN-1.

Conclusion

The IFNAR-1 promoter polymorphisms are not involved in the expression levels of the main IFNAR-1 transcript. The −77VNTR has a regulatory role on the expression of a secondary, truncated, HCC-specific transcript, which in turn coincides with disruptions in cancer-associated pathways and in FN-1 expression modifications.

Next generation sequencing of HIV-1 protease in the PIVOT trial of protease inhibitor monotherapy

•

Next generation sequencing (NGS) was performed on patients on PI monotherapy.

•

It revealed previously unidentified minority variant protease mutations in 3 samples.

•

These mutations do not predict clinically important phenotypic resistance.

Background

The PIVOT trial examined whether patients with suppressed viral load on combination antiretroviral therapy could be safely switched long-term to ritonavir-boosted protease inhibitor (PI) monotherapy. The main trial publication reported that only one of 296 patients allocated to PI monotherapy experienced a loss of drug options due to protease mutations (identified by local Sanger sequencing resistance tests) likely selected by study drug.

Objectives

To assess if we had missed low frequency mutations, using a more sensitive methodology.

Study design

We performed next generation sequencing (NGS) on all available frozen plasma samples with VL >1000 copies/ml from patients who were randomised to PI monotherapy. Assays were performed at Public Health England laboratories using a previously described method.

Median coverage depth was 76,000 and the threshold for detection of minority variants was 2%. Drug susceptibility was predicted using the Stanford HIVdb algorithm.

Results

17 of 26 potential samples, all from different patients, were identified and successfully tested. The median viral load was 6780 copies/ml and the median time since randomisation was 43 weeks. NGS revealed previously unidentified minority variant protease mutations (G73D, I54T, L89V) in three samples, at frequencies ranging between 2% and 10%. None of these mutations predicted intermediate or high level resistance, the trial primary outcome.

Discussion

This report adds to the body of evidence that ritonavir-boosted PI monotherapy, when used as a switch strategy with prompt detection of viral load rebound and early re-introduction of combination therapy, rarely leads to the development of clinically important protease resistance mutations.

Efficient and unbiased metagenomic recovery of RNA virus genomes from human plasma samples

RNA viruses cause significant human pathology and are responsible for the majority of emerging zoonoses. Mainstream diagnostic assays are challenged by their intrinsic diversity, leading to false negatives and incomplete characterisation. New sequencing techniques are expanding our ability to agnostically interrogate nucleic acids within diverse sample types, but in the clinical setting are limited by overwhelming host material and ultra-low target frequency. Through selective host RNA depletion and compensatory protocol adjustments for ultra-low RNA inputs, we are able to detect three major blood-borne RNA viruses – HIV, HCV and HEV. We recovered complete genomes and up to 43% of the genome from samples with viral loads of 10⁴ and 10³ IU/ml respectively. Additionally, we demonstrated the utility of this method in detecting and characterising members of diverse RNA virus families within a human plasma background, some present at very low levels. By applying this method to a patient sample series, we have simultaneously determined the full genome of both a novel subtype of HCV genotype 6, and a co-infecting human pegivirus. This method builds upon earlier RNA metagenomic techniques and can play an important role in the surveillance and diagnostics of blood-borne viruses.

HIV-1 drug resistance mutations emerging on darunavir therapy in PI-naive and -experienced patients in the UK

BACKGROUND

Darunavir is considered to have a high genetic barrier to resistance. Most darunavir-associated drug resistance mutations (DRMs) have been identified through correlation of baseline genotype with virological response in clinical trials. However, there is little information on DRMs that are directly selected by darunavir in clinical settings.

OBJECTIVES

We examined darunavir DRMs emerging in clinical practice in the UK.

PATIENTS AND METHODS

Baseline and post-exposure protease genotypes were compared for individuals in the UK Collaborative HIV Cohort Study who had received darunavir; analyses were stratified for PI history. A selection analysis was used to compare the evolution of subtype B proteases in darunavir recipients and matched PI-naive controls.

RESULTS

Of 6918 people who had received darunavir, 386 had resistance tests pre- and post-exposure. Overall, 2.8% (11/386) of these participants developed emergent darunavir DRMs. The prevalence of baseline DRMs was 1.0% (2/198) among PI-naive participants and 13.8% (26/188) among PI-experienced participants. Emergent DRMs developed in 2.0% of the PI-naive group (4 mutations) and 3.7% of the PI-experienced group (12 mutations). Codon 77 was positively selected in the PI-naive darunavir cases, but not in the control group.

CONCLUSIONS

Our findings suggest that although emergent darunavir resistance is rare, it may be more common among PI-experienced patients than those who are PI-naive. Further investigation is required to explore whether codon 77 is a novel site involved in darunavir susceptibility.

Enhanced surveillance of HIV-1 drug resistance in recently infected MSM in the UK

OBJECTIVES

To determine the prevalence of inferred low-frequency HIV-1 transmitted drug resistance (TDR) in MSM in the UK and its predicted effect on first-line therapy.

METHODS

The HIV-1 pol gene was amplified from 442 newly diagnosed MSM identified as likely recently infected by serological avidity testing in 2011-13. The PCR products were sequenced by next-generation sequencing with a mutation frequency threshold of >2% and TDR mutations defined according to the 2009 WHO surveillance drug resistance mutations list.

RESULTS

The majority (75.6%) were infected with subtype B and 6.6% with rare complex or unique recombinant forms. At a mutation frequency threshold of >20%, 7.2% (95% CI 5.0%-10.1%) of the sequences had TDR and this doubled to 15.8% (95% CI 12.6%-19.6%) at >2% mutation frequency (P < 0.0001). The majority (26/42, 62%) of low-frequency variants were against PIs. The most common mutations detected at >20% and 2%-20% mutation frequency differed for each drug class, these respectively being: L90M (n = 7) and M46IL (n = 10) for PIs; T215rev (n = 9) and D67GN (n = 4) for NRTIs; and K103N (n = 5) and G190E (n = 2) for NNRTIs. Combined TDR was more frequent in subtype B than non-B (OR = 0.38; 95% CI = 0.17-0.88; P = 0.024) and had minimal predicted effect on recommended first-line therapies.

CONCLUSIONS

The data suggest differences in the types of low-frequency compared with majority TDR variants that require a better understanding of the origins and clinical significance of low-frequency variants. This will better inform diagnostic and treatment strategies.

Intrapatient Evolutionary Dynamics of Human Immunodeficiency Virus Type 1 in Individuals Undergoing Alternative Treatment Strategies with Reverse Transcriptase Inhibitors

Structured treatment interruption (STI) has been trialed as an alternative to lifelong antiretroviral therapy (ART). We retrospectively performed single genome sequencing of the HIV-1 pol region from three patients representing different scenarios. They were either failing on continuous therapy (CT-F), failing STI (STI-F), or suppressing on STI (STI-S). Over 460 genomes were generated from three to five different time points over a 2-year period. We found multiple-linked-resistant mutations in both treatment failures. However, the CT-F patient showed a stepwise accumulation of diverse, linked mutations whereas the STI-F patient had lineage turnover between treatment periods with recirculation of wild-type and resistant variants from reservoirs. The STI-F patient showed a 7-fold increase in the third codon position substitution rate relative to the first and second positions compared to a 2-fold increase for CT-F and increased purifying selection in the pol gene (62 vs. 22 sites, respectively). An understanding of intrapatient viral dynamics could guide the future direction of treatment interruption strategies.

Evidence of Self-Sustaining Drug Resistant HIV-1 Lineages Among Untreated Patients in the United Kingdom

BACKGROUND

About 10% of new diagnoses of subtype B human immunodeficiency virus type 1 (HIV-1) in the United Kingdom are with viruses showing transmitted drug resistance (TDR). However, there is discordance between the mutation patterns observed in HIV-infected patients failing therapy and those seen in TDR.

METHODS

We extracted all subtype B HIV-1 pol gene sequences from treatment-naive patients within the United Kingdom HIV Drug Resistance Database sampled between 1997 and 2011 and carrying the most common protease inhibitors, nonnucleoside and nucleotide reverse transcriptase inhibitors TDR mutations, namely, L90M, K103N, and T215Y/F/rev, respectively (n = 1140). Transmission clusters (n ≥ 2 sequences) were identified by maximum-likelihood phylogeny using a genetic distance cutoff of ≤ 1.5%. The time of origin and the basic reproductive number (R0) of clusters were estimated by Bayesian methods.

RESULTS

T215rev was present alone in 47% of the sequences (n = 540), K103N in 31% (n = 359), and L90M in 10% (n = 109). The remaining sequences contained T215Y or combinations of L90M, K103N, and T215rev. Fifty-five percent (n = 624) of the sequences formed highly supported transmission clusters (n = 193) containing between 2 and 15 sequences. The time of origin of 10 large clusters (≥ 8 sequences) was estimated to be between 2000 (1999-2002; 95% highest posterior density [HPD]) and 2006 (2005-2007; 95% HPD). The oldest cluster had persisted for nearly 8 years. All 10 clusters had R0s ranging from 1.3 (0.4-2.5; 95% HPD) to 2.8 (0.6-6.5; 95% HPD).

CONCLUSIONS

A high proportion of the most common TDR in subtype B infections in the United Kingdom is derived by onward transmission from treatment-naive patients.

Protease mutations emerging on darunavir in protease inhibitor-naïve and experienced patients in the UK

INTRODUCTION

Darunavir (DRV) is a preferred agent in treatment guidelines for ART-naïve and experienced patients [1]. It is considered to have a high genetic barrier to resistance and 11 resistance-associated mutations (RAMs) are recognized by IAS-USA [2]. These have largely been identified by analyses examining the correlation between baseline genotype and virological response [3]. However, there is little information on RAMs that are directly selected by DRV, outside of short-term clinical trials. We aimed to identify emerging mutations by comparing the genotypes of individuals before and after DRV exposure.

MATERIALS AND METHODS

The UK HIV Drug Resistance Database was used to identify patients aged over 16 who had received at least 30 days of a DRV-containing regimen. Patients were included if they had a "baseline" resistance test, prior to DRV exposure, and a "repeat" test, either on DRV or within 30 days of stopping this agent. To avoid attributing the effects of other PIs on emerging RAMs to DRV, patients were excluded if they had received another PI for greater than 90 days between the baseline genotype and the start of DRV. The baseline and repeat tests were compared to determine the nature of mutations stratified by PI history.

RESULTS

A total of 5623 patients had DRV, of whom 306 met the inclusion criteria. A total of 228 (74.5%) were male, median age at the start of DRV was 42 years (IQR 37-47), and half had subtype B infection. The mode of transmission was homosexual contact for 50%, heterosexual for 38%, and 3% were injection drug users. The median CD4 count at the start of DRV was 257 cells/mm(3) (IQR 94-453). A total of 149 patients (49%) had a history of PI use prior to DRV, and 157 (51%) were PI-naïve. The most common previous PIs were lopinavir, atazanavir, and saquinavir. Baseline DRV RAMs were present in 1 (0.6%) PI-naïve and 20 (13.4%) PI-experienced patients. Mutations emerged under DRV pressure in a further 3 (1.9%) PI-naïve patients, and in 7 (4.7%) PI-experienced patients, 5 of whom had other DRV RAMs present at baseline (Table 1). The median time from the start of DRV to the repeat test was 196 days for PI-naïve patients and 296 days for PI-experienced.

CONCLUSIONS

PI-experienced patients had a greater prevalence of DRV RAMs at baseline than PI-naïve individuals, probably due to the fact that some DRV RAMs can be selected by other PIs. This group also accumulated more RAMs during DRV exposure, possibly because previous PIs had caused minority variants which then emerged on DRV therapy. Overall, only 10 patients accumulated 16 RAMs, which supports the perception that DRV has a high genetic barrier to resistance. Repeat genotyping in the case of virological failure on DRV may still be warranted to detect emerging resistance and guide management decisions.

HIV-1 subtype influences susceptibility and response to monotherapy with the protease inhibitor lopinavir/ritonavir

OBJECTIVE

PI susceptibility results from a complex interplay between protease and Gag proteins, with Gag showing wide variation across HIV-1 subtypes. We explored the impact of pre-treatment susceptibility on the outcome of lopinavir/ritonavir monotherapy.

METHODS

Treatment-naive individuals who experienced lopinavir/ritonavir monotherapy failure from the MONARK study were matched (by subtype, viral load and baseline CD4 count) with those who achieved virological response ('successes'). Successes were defined by viral load <400 copies/mL after week 24 and <50 copies/mL from week 48 to week 96. Full-length Gag-protease was amplified from patient samples for in vitro phenotypic susceptibility testing, with susceptibility expressed as fold change (FC) relative to a subtype B reference strain.

RESULTS

Baseline lopinavir susceptibility was lower in viral failures compared with viral successes, but the differences were not statistically significant (median lopinavir susceptibility: 4.4 versus 8.5, respectively, P = 0.17). Among CRF02_AG/G patients, there was a significant difference in lopinavir susceptibility between the two groups (7.1 versus 10.4, P = 0.047), while in subtype B the difference was not significant (2.7 versus 3.4, P = 0.13). Subtype CRF02_AG/G viruses had a median lopinavir FC of 8.7 compared with 3.1 for subtype B (P = 0.001).

CONCLUSIONS

We report an association between reduced PI susceptibility (using full-length Gag-protease sequences) at baseline and subsequent virological failure on lopinavir/ritonavir monotherapy in antiretroviral-naive patients harbouring subtype CRF02_AG/G viruses. We speculate that this may be important in the context of suboptimal adherence in determining viral failure.

Phenotypic characterization of virological failure following lopinavir/ritonavir monotherapy using full-length Gag-protease genes

Objectives

Major protease mutations are rarely observed following first-line failure with PIs and interpretation of genotyping results in this context may be difficult. We performed extensive phenotyping of viruses from five patients failing lopinavir/ritonavir monotherapy in the MONARK study without major PI mutations by standard genotyping.

Methods

Phenotypic susceptibility testing and viral infectivity assessments were performed using a single-cycle assay and fold changes (FC) relative to a lopinavir-susceptible reference strain were calculated.

Results

>10-fold reduced baseline susceptibility to lopinavir occurred in two of five patients and >5-fold in another two. Four of five patients exhibited phylogenetic evidence of a limited viral evolution between baseline and failure, with amino acid changes at drug resistance-associated positions in one: T81A emerged in Gag with M36I in the protease gene, correlating with a reduction in lopinavir susceptibility from FC 7 (95% CI 6–8.35) to FC 13 (95% CI 8.11–17.8). Reductions in darunavir susceptibility (>5 FC) occurred in three individuals.

Discussion

This study suggests both baseline reduced susceptibility and evolution of resistance could be contributing factors to PI failure, despite the absence of classical PI resistance mutations by standard testing methods. Use of phenotyping also reveals lower darunavir susceptibility, warranting further study as this agent is commonly used following lopinavir failure.

Contribution of Gag and protease to variation in susceptibility to protease inhibitors between different strains of subtype B human immunodeficiency virus type 1

Recent reports have shown that human immunodeficiency virus type 1 (HIV-1) Gag can directly affect susceptibility to protease inhibitors (PIs) in the absence of known resistance mutations in protease. Inclusion of co-evolved Gag alongside protease in phenotypic drug susceptibility assays can alter PI susceptibility in comparison with protease with a WT Gag. Using a single-replication-cycle assay encompassing full-length Gag together with protease we demonstrated significant variation in PI susceptibility between a number of PI-naïve subtype B viruses. Six publicly available subtype B molecular clones, namely HXB2, NL4-3, SF2, YU2, JRFL and 89.6, displayed up to nine-fold reduced PI susceptibility in comparison with the assay reference strain. For two molecular clones, YU2 and JRFL, Gag contributed solely to the observed reduction in susceptibility, with the N-terminal region of Gag contributing significantly. Gag and protease from treatment-naïve, patient-derived viruses also demonstrated significant variation in susceptibility, with up to a 17-fold reduction to atazanavir in comparison with the assay reference strain. In contrast to the molecular clones, protease was the main determinant of the reduced susceptibility. Common polymorphisms in protease, including I13V, L63P and A71T, were shown to contribute to this reduction in PI susceptibility, in the absence of major resistance mutations. This study demonstrated significant variation in PI susceptibility of treatment-naïve patient viruses, and provided further evidence of the independent role of Gag, the protease substrate and in particular the N-terminus of Gag in PI susceptibility. It also highlighted the importance of considering co-evolved Gag and protease when assessing PI susceptibility.

Low frequency of genotypic resistance in HIV-1-infected patients failing an atazanavir-containing regimen: a clinical cohort study

OBJECTIVES

To determine protease mutations that develop at viral failure for protease inhibitor (PI)-naive patients on a regimen containing the PI atazanavir.

METHODS

Resistance tests on patients failing atazanavir, conducted as part of routine clinical care in a multicentre observational study, were randomly matched by subtype to resistance tests from PI-naive controls to account for natural polymorphisms. Mutations from the consensus B sequence across the protease region were analysed for association and defined using the IAS-USA 2011 classification list.

RESULTS

Four hundred and five of 2528 (16%) patients failed therapy containing atazanavir as a first PI over a median (IQR) follow-up of 1.76 (0.84-3.15) years and 322 resistance tests were available for analysis. Recognized major atazanavir mutations were found in six atazanavir-experienced patients (P < 0.001), including I50L and N88S. The minor mutations most strongly associated with atazanavir experience were M36I, M46I, F53L, A71V, V82T and I85V (P < 0.05). Multiple novel mutations, I15S, L19T, K43T, L63P/V, K70Q, V77I and L89I/T/V, were also associated with atazanavir experience.

CONCLUSIONS

Viral failure on atazanavir-containing regimens was not common and major resistance mutations were rare, suggesting that adherence may be a major contributor to viral failure. Novel mutations were described that have not been previously documented.

Dynamics of HIV type 1 recombination following superinfection

There are currently few detailed studies describing HIV-1 recombination events or the potential impact of recombination on drug resistance. We describe here the viral recombination dynamics in a drug-naive patient initially infected with a circulating recombinant form 19 (CRF19) virus containing transmitted drug resistance mutations followed by superinfection with "wild-type" subtype B virus. Single genome analysis showed replacement of the primary CRF19 virus by recombinants of the CRF19 virus and the superinfecting subtype B virus. The CRF19/B recombinant virus dominating after superinfection had lost drug resistance mutations and at no time was the superinfecting subtype B variant found to be dominant in blood plasma. Furthermore, the detection of recombinant viruses in seminal plasma indicates the potential for onward transmission of these strains.

Subtype-specific differences in the development of accessory mutations associated with high-level resistance to HIV-1 nucleoside reverse transcriptase inhibitors

OBJECTIVES

To identify accessory mutations associated with high-level resistance to reverse transcriptase (RT) inhibitors in HIV-1 subtypes B and C.

METHODS

Changes relative to the wild-type for codons 1-400 of RT were analysed from treatment-experienced patients infected with subtypes B (5464 patients) and C (1920 patients). Positions associated with the accumulation of mutations conferring resistance to thymidine analogues and to non-nucleoside RT inhibitors (NNRTIs) were identified. A subtype-specific single-replication cycle drug susceptibility assay was used to determine whether some of the mutations affected drug susceptibility or viral infectivity.

RESULTS

In subtype B, mutations at 31 and 26 positions were associated with the accumulation of thymidine analogue mutations (TAMs) and NNRTI mutations, respectively; in subtype C, 18 and 13 positions were identified, respectively. Amino acid changes at the following positions were differentially associated with (i) the accumulation of 0-4+ TAMs in subtypes B and C (away from consensus): 43 (27.0% B versus 2.5% C); 118 (36.4% B versus 16.2% C); 135 (12.5% B versus 28.0% C); and 326 (2.6% towards consensus in B versus 7.6% away in C) and (ii) the accumulation of 0-3+ NNRTI mutations (away from consensus): 43 (10.2% B versus 0.5% C); and 68 (5.2% B versus 10.3% C). Codon changes K43E, E44D and V118I were found to have no effect on susceptibility to three NRTIs with or without TAMs in either subtype; however, some accessory mutations had subtype-specific effects on viral infectivity.

CONCLUSIONS

Differences between subtypes B and C were observed in the development and effect of accessory mutations associated with high-level resistance to RT inhibitors.

Phylodynamic and phylogeographic patterns of the HIV type 1 subtype F1 parenteral epidemic in Romania

In the late 1980s an HIV-1 epidemic emerged in Romania that was dominated by subtype F1. The main route of infection is believed to be parenteral transmission in children. We sequenced partial pol coding regions of 70 subtype F1 samples from children and adolescents from the PENTA-EPPICC network of which 67 were from Romania. Phylogenetic reconstruction using the sequences and other publically available global subtype F sequences showed that 79% of Romanian F1 sequences formed a statistically robust monophyletic cluster. The monophyletic cluster was epidemiologically linked to parenteral transmission in children. Coalescent-based analysis dated the origins of the parenteral epidemic to 1983 [1981-1987; 95% HPD]. The analysis also shows that the epidemic's effective population size has remained fairly constant since the early 1990s suggesting limited onward spread of the virus within the population. Furthermore, phylogeographic analysis suggests that the root location of the parenteral epidemic was Bucharest.

Molecular and epidemiological characterization of HIV-1 infection networks involving transmitted drug resistance mutations in Northern Greece

OBJECTIVES

To determine the contribution of transmission clusters to transmitted drug resistance (TDR) in newly diagnosed antiretroviral-naive HIV-1-infected patients in Northern Greece during 2000-07.

METHODS

The prevalence of TDR was estimated in 369 individuals who were diagnosed with HIV-1 infection in the period 2000-07 at the National AIDS Reference Laboratory of Northern Greece. Phylogenetic analysis was performed using a maximum likelihood method on partial pol sequences. TDR was defined in accordance with the surveillance drug resistance mutation list (2009 update).

RESULTS

The overall prevalence of TDR in our population was 12.5% [46/369, 95% confidence interval (CI) 9.1%-15.8%], comprising 7.6% (28/369) resistant to nucleoside reverse transcriptase inhibitors, 5.4% (20/369) resistant to non-nucleoside reverse transcriptase inhibitors and 3.3% (12/369) resistant to protease inhibitors. Dual class resistance was identified in 3.8% (14/369). Infection with subtype A was the sole predictor associated with TDR in multivariate analysis (odds ratio 2.15, 95% CI 1.10-4.19, P = 0.025). Phylogenetic analyses revealed three statistically robust transmission clusters involving drug-resistant strains, including one cluster of 12 patients, 10 of whom were infected with a strain carrying both T215 revertants and Y181C mutations.

CONCLUSIONS

Our findings underline the substantial impact of transmission networks on TDR in our population.

The evolution of HIV-1 reverse transcriptase in route to acquisition of Q151M multi-drug resistance is complex and involves mutations in multiple domains

Background

The Q151M multi-drug resistance (MDR) pathway in HIV-1 reverse transcriptase (RT) confers reduced susceptibility to all nucleoside reverse transcriptase inhibitors (NRTIs) excluding tenofovir (TDF). This pathway emerges after long term failure of therapy, and is increasingly observed in the resource poor world, where antiretroviral therapy is rarely accompanied by intensive virological monitoring. In this study we examined the genotypic, phenotypic and fitness correlates associated with the development of Q151M MDR in the absence of viral load monitoring.

Results

Single-genome sequencing (SGS) of full-length RT was carried out on sequential samples from an HIV-infected individual enrolled in ART rollout. The emergence of Q151M MDR occurred in the order A62V, V75I, and finally Q151M on the same genome at 4, 17 and 37 months after initiation of therapy, respectively. This was accompanied by a parallel cumulative acquisition of mutations at 20 other codon positions; seven of which were located in the connection subdomain. We established that fourteen of these mutations are also observed in Q151M-containing sequences submitted to the Stanford University HIV database. Phenotypic drug susceptibility testing demonstrated that the Q151M-containing RT had reduced susceptibility to all NRTIs except for TDF. RT domain-swapping of patient and wild-type RTs showed that patient-derived connection subdomains were not associated with reduced NRTI susceptibility. However, the virus expressing patient-derived Q151M RT at 37 months demonstrated ~44% replicative capacity of that at 4 months. This was further reduced to ~22% when the Q151M-containing DNA pol domain was expressed with wild-type C-terminal domain, but was then fully compensated by coexpression of the coevolved connection subdomain.

Conclusions

We demonstrate a complex interplay between drug susceptibility and replicative fitness in the acquisition Q151M MDR with serious implications for second-line regimen options. The acquisition of the Q151M pathway occurred sequentially over a long period of failing NRTI therapy, and was associated with mutations in multiple RT domains.

Patterns of resistance development with integrase inhibitors in HIV

Raltegravir, the only integrase (IN) inhibitor approved for use in HIV therapy, has recently been licensed. Raltegravir inhibits HIV-1 replication by blocking the IN strand transfer reaction. More than 30 mutations have been associated with resistance to raltegravir and other IN strand transfer inhibitors (INSTIs). The majority of the mutations are located in the vicinity of the IN active site within the catalytic core domain which is also the binding pocket for INSTIs. High-level resistance to INSTIs primarily involves three independent mutations at residues Q148, N155, and Y143. The mutations significantly affect replication capacity of the virus and are often accompanied by other mutations that either improve replication fitness and/or increase resistance to the inhibitors. The pattern of development of INSTI resistance mutations has been extensively studied in vitro and in vivo. This has been augmented by cell-based phenotypic studies and investigation of the mechanisms of resistance using biochemical assays. The recent elucidation of the structure of the prototype foamy virus IN, which is closely related to HIV-1, in complex with INSTIs has greatly enhanced our understanding of the evolution and mechanisms of IN drug resistance.

Real-time PCR analysis of HIV-1 replication post-entry events

The reverse transcriptase enzyme plays an essential role in the HIV-1 life cycle by converting a single-stranded viral RNA genome into a double-stranded viral DNA through a complex process known as reverse transcription. The resulting double-stranded DNA is integrated into the host chromosome to form a provirus. A small proportion of the viral DNAs form dead-end circular products, which nevertheless can serve as useful surrogate markers for monitoring viral replication. Utilizing real-time PCR technology, it is possible to track and quantify different stages of the reverse transcription process, the proviruses, and the nonintegrated dead-end reverse transcription products.

Human immunodeficiency virus type 1 cDNAs produced in the presence of APOBEC3G exhibit defects in plus-strand DNA transfer and integration

Encapsidation of host restriction factor APOBEC3G (A3G) into vif-deficient human immunodeficiency virus type 1 (HIV-1) blocks virus replication at least partly by C-to-U deamination of viral minus-strand DNA, resulting in G-to-A hypermutation. A3G may also inhibit HIV-1 replication by reducing viral DNA synthesis and inducing viral DNA degradation. To gain further insight into the mechanisms of viral inhibition, we examined the metabolism of A3G-exposed viral DNA. We observed that an overall 35-fold decrease in viral infectivity was accompanied by a five- to sevenfold reduction in viral DNA synthesis. Wild-type A3G induced an additional fivefold decrease in the amount of viral DNA that was integrated into the host cell genome and similarly reduced the efficiency with which HIV-1 preintegration complexes (PICs) integrated into a target DNA in vitro. The A3G C-terminal catalytic domain was required for both of these antiviral activities. Southern blotting analysis of PICs showed that A3G reduced the efficiency and specificity of primer tRNA processing and removal, resulting in viral DNA ends that are inefficient substrates for integration and plus-strand DNA transfer. However, the decrease in plus-strand DNA transfer did not account for all of the observed decrease in viral DNA synthesis associated with A3G. These novel observations suggest that HIV-1 cDNA produced in the presence of A3G exhibits defects in primer tRNA processing, plus-strand DNA transfer, and integration.

Mutations in the RNase H primer grip domain of murine leukemia virus reverse transcriptase decrease efficiency and accuracy of plus-strand DNA transfer

The RNase H primer grip of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) contacts the DNA primer strand and positions the template strand near the RNase H active site, influencing RNase H cleavage efficiency and specificity. Sequence alignments show that 6 of the 11 residues that constitute the RNase H primer grip have functional equivalents in murine leukemia virus (MLV) RT. We previously showed that a Y586F substitution in the MLV RNase H primer grip resulted in a 17-fold increase in substitutions within 18 nucleotides of adenine-thymine tracts, which are associated with a bent DNA conformation. To further determine the effects of the MLV RNase H primer grip on replication fidelity and viral replication, we performed additional mutational analysis. Using either beta-galactosidase (lacZ) or green fluorescent protein (GFP) reporter genes, we found that S557A, A558V, and Q559L substitutions resulted in statistically significant increases in viral mutation rates, ranging from 2.1- to 3.8-fold. DNA sequencing analysis of nonfluorescent GFP clones indicated that the mutations in RNase H primer grip significantly increased the frequency of deletions between the primer-binding site (PBS) and sequences downstream of the PBS. In addition, quantitative real-time PCR analysis of reverse transcription products revealed that the mutant RTs were substantially inefficient in plus-strand DNA transfer relative to the wild-type control. These results indicate that the MLV RNase H primer grip is an important determinant of in vivo fidelity of DNA synthesis and suggest that the mutant RT was unable to copy through the DNA-RNA junction of the minus-strand DNA and the tRNA because of its bent conformation resulting in error-prone plus-strand DNA transfer.

Human apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like 3G (APOBEC3G) is incorporated into HIV-1 virions through interactions with viral and nonviral RNAs

Apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like 3G (APOBEC3G) is a host cytidine deaminase that is packaged into virions and confers resistance to retroviral infection. APOBEC3G deaminates deoxycytidines in minus strand DNA to deoxyuridines, resulting in G to A hypermutation and viral inactivation. Human immunodeficiency virus type 1 (HIV-1) virion infectivity factor counteracts the antiviral activity of APOBEC3G by inducing its proteosomal degradation and preventing virion incorporation. To elucidate the mechanism of viral suppression by APOBEC3G, we developed a sensitive cytidine deamination assay and analyzed APOBEC3G virion incorporation in a series of HIV-1 deletion mutants. Virus-like particles derived from constructs in which pol, env, and most of gag were deleted still contained high levels of cytidine deaminase activity; in addition, coimmunoprecipitation of APOBEC3G and HIV-1 Gag in the presence and absence of RNase A indicated that the two proteins do not interact directly but form an RNase-sensitive complex. Viral particles lacking HIV-1 genomic RNA which were generated from the gag-pol expression constructs pC-Help and pSYNGP packaged APOBEC3G at 30-40% of the wild-type level, indicating that interactions with viral RNA are not necessary for incorporation. In addition, viral particles produced from an nucleocapsid zinc finger mutant contained approximately 1% of the viral genomic RNA but approximately 30% of the cytidine deaminase activity. The reduction in APOBEC3G incorporation was equivalent to the reduction in the total RNA present in the nucleocapsid mutant virions. These results indicate that interactions with viral proteins or viral genomic RNA are not essential for APOBEC3G incorporation and suggest that APOBEC3G interactions with viral and nonviral RNAs that are packaged into viral particles are sufficient for APOBEC3G virion incorporation.

Reovirus mu2 protein determines strain-specific differences in the rate of viral inclusion formation in L929 cells

Reovirus infection induces the formation of large cytoplasmic inclusions that serve as the major site of viral assembly. Reovirus strains type 3 Dearing (T3D) and type 1 Lang (T1L) differ in the rate of inclusion formation in L929 cells. The median time of inclusion formation is 18 h in cells infected with T3D and 39 h in cells infected with T1L. Using reassortant viruses that contain combinations of gene segments derived from T1L and T3D, we found that the M1 gene, which encodes the mu2 protein, is the primary determinant of the rate of inclusion formation. The S3 gene, which encodes the nonstructural protein sigmaNS, plays a secondary role in this process. The subcellular location of the mu2 protein was determined by confocal laser scanning microscopy using dual-fluorescence labeling of mu2 and the outer-capsid protein mu1/mu1C. In virus-infected cells, mu2 protein colocalized with other viral proteins in inclusions and was also distributed diffusely in the cytoplasm and nucleus. Expression of recombinant T1L and T3D mu2 proteins resulted in the formation of protein complexes resembling inclusions in both the cytoplasm and the nucleus with kinetics that reflected the strain of origin. The median time of mu2 protein complex formation was 22 h in cells transfected with the T3D M1 gene and 43 h in cells transfected with the T1L M1 gene. These findings suggest that the mu2 protein influences the rate of inclusion formation and contributes to inclusion morphogenesis. The requirement of mu2 protein in inclusion formation was tested by determining the subcellular localization of mu2 in cells infected with temperature-sensitive (ts) mutants that are defective in viral assembly. In contrast to infection with wild-type virus, mu2 did not colocalize with mu1/mu1C protein in subcellular structures that formed in cells infected at nonpermissive temperature with ts mutants tsH11.2, tsC447, and tsG453 with mutations in the M1, S2, and S4 genes, respectively. These results suggest that despite the role of the mu2 protein in controlling the rate of inclusion formation, this process is a concerted function of several reovirus proteins.