Characterization of mutation spectra with ultra-deep pyrosequencing: application to HIV-1 drug resistance

Antiviral activity and active components of the leaves from Sabia parviflora Wall. ex Roxb

Sabia parviflora (SP, "xiao hua qing feng teng" in Chinese) was recorded as an important ethnic medicine to be used for treating viral hepatitis. The antiviral activity of four SP extracts and potent antiviral compounds evaluated with cathepsin L protease (Cat L PR) and HIV-1 protease (HIV-1 PR). UPLC-HRMS was used for identifying the bioactive components. In addition, the possible inhibitory mechanism of the identified compounds on viral protease was further discussed by molecular docking. As a result, four extracts of SP exhibited inhibitory activity of HIV-1 PR and Cat L PR with IC50 range from 0.015 to 0.80 mg/mL. Meanwhile, six compounds inhibited HIV-1 PR with IC50 range from 0.032 to 0.80 mg/mL. Moreover, procyanidin B2 had good affinity for HIV-1 PR and CatL PR protein, respectively. These findings suggest S. parviflora leaves can be used for treating HIV and procyanidin B2 may play a role in antiviral protease.

Role of Proviral HIV-1 DNA Genotyping for People Living with HIV (PLWH) Who Had Low-Level Viremia While Receiving Antiretroviral Therapy

Objective

To analyze the antiretroviral resistance in people living with HIV (PLWH) who developed low-level viremia (LLV) during antiretroviral therapy (ART) via sequencing of their HIV-1 proviral DNA and RNA and comparisons of their proviral DNA genotyping data with their past and synchronous RNA genotyping data.

Patients and Methods

PLWH with LLV while receiving ART for 6 months or longer from January 2020 to September 2021 were included. HIV-1 proviral DNA and RNA were extracted from white-blood cells and concentrated plasma by ultracentrifugation, respectively, and HIV-1 pol gene fragments were amplified and sequenced. The concordance in the detection of resistance-associated mutations (RAMs) were examined between proviral DNA vs past RNA genotyping and proviral DNA vs synchronous RNA genotyping.

Results

Of the 150 PLWH with LLV, 117 proviral DNA pol sequences detected in 105 PLWH were successfully amplified and RAMs were present in 27.6% and the rate of RAMs conferring low-level or greater resistance to antiretrovirals examined was 17.1%. Fifty-six and 57 PLWH had results of past and synchronous RNA genotyping, respectively, for comparisons with those of proviral DNA genotyping; and the concordance rates were 76.8% and 75.4%, respectively. However, proviral DNA genotyping lost than gained partial information on antiretroviral resistance compared with past or synchronous RNA genotyping.

Conclusion

We found that the concordance between proviral DNA and past and synchronous RNA genotyping was moderate. Proviral DNA genotyping lost than gained more information on antiretroviral resistance compared with past or synchronous RNA genotyping. To optimize ART in PLWH with LLV, antiretroviral resistance profile should be interpreted in combination with proviral DNA and RNA genotyping and a comprehensive review of previous treatment history.

Lessons from resistance analysis in clinical trials of IV zanamivir

Influenza infection causes substantial morbidity and mortality during seasonal epidemics and pandemics. Antivirals, including neuraminidase inhibitors, play an important role in the treatment of severely ill patients infected with influenza. Resistance is a key factor that can affect the efficacy of neuraminidase inhibitors (NAIs). It is a recommendation by regulatory authorities to monitor for resistance during the development of anti-influenza medications. An additional requirement by regulators is to examine amino acid sequences for minority species harbouring resistance substitutions. In a Phase III study of intravenous (IV) zanamivir respiratory samples were analysed for the presence of resistant quasi species using Next Generation Sequencing (NGS). In this study ten resistance substitutions, two of which were treatment emergent, were detected by NGS that otherwise would not have been detectable by Sanger sequencing. None of the substitutions were present at any other timepoints analysed. The effect these mutations have on clinical response is difficult to characterize; in fact, all patients from which these variants were isolated had a successful clinical outcome and the effect on clinical response was therefore likely minimal. Although NGS is becoming a routine method for nucleic acid sequencing and will detect substitutions previously undetected by Sanger sequencing, the value of this technique in identifying minority species with resistance substitutions that are clinically meaningful remains to be demonstrated, particularly with acute infections such as influenza.

Quantification of in vitro replication kinetics of Alagoas vesiculovirus isolates by digital droplet RT-PCR

Vesicular stomatitis caused by Alagoas vesiculovirus (VSAV) has generated disease outbreaks in Brazil, mainly in the northeast region. Phylogenetic studies divide the isolates into three distinct genotypes (A, B, and C). However, there is no description of how this genetic divergence reflects on the phenotype of VSAV isolates such as in vitro replication fitness. Therefore, the objective of this work was to evaluate the ability of three distinct genotypes of Brazilian isolates of VSAV to grow in different cell-culture lines (BHK-21, Vero, and NCI-H1299). Quantification of viral RNA was performed using RT-PCR digital droplet from supernatant of cell culture collected every 4 h for a period of 24 h of viral growth in three different cell lines (BHK-21, Vero, and NCI-H1299). It was observed that the genotype C isolate has the lowest replication efficiency among the three analyzed viruses, without major changes in the copies of viral RNA over the entire time of the study.

Towards Next-Generation Sequencing for HIV-1 Drug Resistance Testing in a Clinical Setting

The HIV genotypic resistance test (GRT) is a standard of care for the clinical management of HIV/AIDS patients. In recent decades, population or Sanger sequencing has been the foundation for drug resistance monitoring in clinical settings. However, the advent of high-throughput or next-generation sequencing has caused a paradigm shift towards the detection and characterization of low-abundance covert mutations that would otherwise be missed by population sequencing. This is clinically significant, as these mutations can potentially compromise the efficacy of antiretroviral therapy, causing poor virologic suppression. Therefore, it is important to develop a more sensitive method so as to reliably detect clinically actionable drug-resistant mutations (DRMs). Here, we evaluated the diagnostic performance of a laboratory-developed, high-throughput, sequencing-based GRT using 103 archived clinical samples that were previously tested for drug resistance using population sequencing. As expected, high-throughput sequencing found all the DRMs that were detectable by population sequencing. Significantly, 78 additional DRMs were identified only by high-throughput sequencing, which is statistically significant based on McNemar's test. Overall, our results complement previous studies, supporting the notion that the two methods are well correlated, and the high-throughput sequencing method appears to be an excellent alternative for drug resistance testing in a clinical setting.

Detection of Chimeric Cellular: HIV mRNAs Generated Through Aberrant Splicing in HIV-1 Latently Infected Resting CD4+ T Cells

Latent HIV-1 provirus in infected individuals on suppressive therapy does not always remain transcriptionally silent. Both HIV-1 LTR and human gene promoter derived transcriptional events can contribute HIV-1 sequences to the mRNA produced in the cell. In addition, chimeric cellular:HIV mRNA can arise through readthrough transcription and aberrant splicing. Using target enrichment coupled to the Illumina Mi-Seq and PacBio RS II platforms, we show that 3’ LTR activation is frequent in latently infected cells from both the CCL19-induced primary cell model of HIV-1 latency as well as ex vivo samples. In both systems of latent HIV-1 infection, we detected several chimeric species that were generated via activation of a cryptic splice donor site in the 5’ LTR of HIV-1. Aberrant splicing involving the major HIV-1 splice donor sites, SD1 and SD4 disrupts post-transcriptional processing of the gene in which HIV-1 is integrated. In the primary cell model of HIV-1 latency, Tat-encoding sequences are incorporated into the chimeric mRNA transcripts through the use of SD4. Our study unravels clues to the characteristics of HIV-1 integrants that promote formation of chimeric cellular:HIV mRNA and improves the understanding of the HIV-1 RNA footprint in latently infected cells.

Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2

PRS17, a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PRS17 with inhibitor Darunavir (DRV), two complexes of WT PR and PRS17 with substrate analogue CA-p2, two unligand WT PR and PRS17), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PRS17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

CODEHOP-Mediated PCR Improves HIV-1 Genotyping and Detection of Variants by MinION Sequencing

HIV-1 is genetically heterogeneous, having different subtypes and circulating recombinant forms (CRFs). HIV-1 genotyping is used to determine drug resistance profiles and is based on the use of a mixture of consensus and degenerate primers targeting the pol gene. However, the use of this type of primers is associated with either PCR bias or PCR failure. Consensus-degenerate hybrid oligonucleotide primers (CODEHOPs) can detect and identify unknown and distantly related gene sequences by PCR. CODEHOPs designed using different HIV-1 subtypes and CRFs were evaluated for HIV-1 genotyping by Sanger and MinION sequencing. A total of 321 plasma samples were used for the validation of CODEHOP-mediated HIV-1 genotyping. CODEHOP-mediated PCR showed 100% sensitivity and specificity, with limits of detection and genotyping below 200 copies/ml. The head-to-head evaluation of CODEHOP-mediated PCR and standard PCR showed 97 to 98% and 82 to 84% PCR success rates, respectively. There was 100% agreement between the CODEHOP and the reference method in the drug resistance profiles determined by Sanger-based sequencing. Using MinION sequencing, the CODEHOP-mediated PCR scheme resulted in better depth of genome coverage and detection of more drug resistance variants in the protease and reverse transcriptase genes than the standard amplification scheme. The overall prevalences of drug resistance mutations were 17.1% in treatment-experienced patients and 1.2% in treatment-naive patients. They were mainly associated with resistance to reverse transcriptase inhibitors and were linked to virological failure and the patient’s treatment history. Findings from this study suggest that the performance of HIV-1 genotyping is improved by using CODEHOP-mediated PCR.

IMPORTANCE HIV-1 drug resistance is the main cause of treatment failure. Regular surveillance of resistance-associated mutations in HIV-1 genomes is essential for the optimal management of HIV-1 infections. Due to HIV-1’s genetic diversity, different HIV-1 genotypes are circulating worldwide. Standard primers used in the amplification of HIV-1 RNA have not been designed to cover all HIV-1 genotypes and are the main cause of amplification and drug resistance test failure. In this study, new sets of PCR primers targeting the protease, reverse transcriptase, and integrase genes were designed using the CODEHOP approach. They were compared to primers recommended in part by WHO for drug resistance testing using in-house PCR. Unsuccessful HIV-1 RNA amplification was less likely to occur with CODEHOP primers, leading to fewer test failures and lower cost. Furthermore, CODEHOP primers were more effective than standard primers for the detection of minority resistant variants by MinION sequencing.

Investigation of drug resistance against protease, reverse transcriptase, and integrase inhibitors by next-generation sequencing in HIV-positive patients

Our aim was to investigate the mutations in protease (PR), reverse transcriptase (RT), and integrase (IN) gene regions in human immunodeficiency virus (HIV) using a single amplicon via next-generation sequencing (NGS). The study included plasma samples from 49 HIV-1-positive patients, which were referred for HIV-1 drug resistance testing during 2017. A nested polymerase chain reaction (PCR) was performed after the RNA extraction and one-step reverse transcription stages. The sequencing of the HIV genome in the PR, RT, and IN gene regions was carried out using MiSeq NGS technology. Sanger sequencing (SS) was used to analyze resistance mutations in the PR and RT gene regions using a ViroSeq HIV-1 Genotyping System. PCR products were analyzed with an ABI3500 GeneticAnalyzer (Applied Biosystems). Resistance mutations detected with NGS at frequencies above 20% were identical to the SS results. Resistance to at least one antiretroviral (ARV) drug was 22.4% (11 of 49) with NGS and 10.2% (5 of 49) with SS. At least one low-frequency resistance mutation was detected in 18.3% (9 of 49) of the samples. Low-frequency resistance mutations resulted in virological failure in only one patient. The cost of the analyses was reduced by sample pooling and multiplex analysis using the MiSeq system. This is the first study in Turkey to use NGS technologies for the detection of resistance mutations in all three gene (PR, RT, IN) regions using a single amplicon. Our findings suggest that NGS is more sensitive and cost-effective than the SS method.

Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea

Whole-genome sequencing of infectious agents enables the identification and characterization of emerging viruses. The MinION device is a portable sequencer that allows real-time sequencing in fields or hospitals. Hantaan orthohantavirus (Hantaan virus, HTNV), harbored by Apodemus agrarius, causes hemorrhagic fever with renal syndrome (HFRS) and poses a critical public health threat worldwide. In this study, we aimed to evaluate the feasibility of using nanopore sequencing for whole-genome sequencing of HTNV from samples having different viral copy numbers. Amplicon-based next-generation sequencing was performed in A. agrarius lung tissues collected from the Republic of Korea. Genomic sequences of HTNV were analyzed based on the viral RNA copy numbers. Amplicon-based nanopore sequencing provided nearly full-length genomic sequences of HTNV and showed sufficient read depth for phylogenetic analysis after 8 h of sequencing. The average identity of the HTNV genome sequences for the nanopore sequencer compared to those of generated from Illumina MiSeq revealed 99.8% (L and M segments) and 99.7% (S segment) identities, respectively. This study highlights the potential of the portable nanopore sequencer for rapid generation of accurate genomic sequences of HTNV for quicker decision making in point-of-care testing of HFRS patients during a hantavirus outbreak.

Insight into the drug resistance whole genome of Mycobacterium tuberculosis isolates from Khyber Pakhtunkhwa, Pakistan

Whole genome sequencing (WGS) is one of the most reliable methods for detection of drug resistance, genetic diversity in other virulence factor and also evolutionary dynamics of Mycobacterium tuberculosis complex (MTBC). First-line anti-tuberculosis drugs are the major weapons against Mycobacterium tuberculosis (MTB). However, the emergence of drug resistance remained a major obstacle towards global tuberculosis (TB) control program 2030, especially in high burden countries including Pakistan. To overcome the resistance and design potent drugs, genomic variations in drugs targets as well as in the virulence and evolutionary factors might be useful for better understanding and designing potential inhibitors. Here we aimed to find genomic variations in the first-line drugs targets, along with other virulence and evolutionary factors among the circulating isolates in Khyber Pakhtunkhwa, Pakistan. Samples were collected and drug susceptibility testing (DST) was performed as per WHO standard. The resistance samples were subjected to WGS. Among the five whole genome sequences, three samples (NCBI BioProject Accession: PRJNA629298, PRJNA629388) harbored 1997, 1162, and 2053 mutations. Some novel mutations have been detected in drugs targets. Similarly, numerous novel variants have also been detected in virulency and evolutionary factors, PE, PPE, and secretory system of MTB isolates. Exploring the genomic variations among the circulating isolates in geographical specific locations might be useful for future drug designing. To the best of our knowledge, this is the first study that provides useful data regarding the insight genomic variations in virulency, evolutionary factors including ESX and PE/PPE as well as drug targets, for better understanding and management of TB in a WHO declared high burden country.

Next Generation Sequencing for HIV-1 Drug Resistance Testing-A Special Issue Walkthrough

Drug resistance remains a global challenge in the fight against the HIV pandemic [...].

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.

Amplicon-Based, Next-Generation Sequencing Approaches to Characterize Single Nucleotide Polymorphisms of Orthohantavirus Species

Whole-genome sequencing (WGS) of viruses from patient or environmental samples can provide tremendous insight into the epidemiology, drug resistance or evolution of a virus. However, we face two common hurdles in obtaining robust sequence information; the low copy number of viral genomes in specimens and the error introduced by WGS techniques. To optimize detection and minimize error in WGS of hantaviruses, we tested four amplification approaches and different amplicon pooling methods for library preparation and examined these preparations using two sequencing platforms, Illumina MiSeq and Oxford Nanopore Technologies MinION. First, we tested and optimized primers used for whole segment PCR or one kilobase amplicon amplification for even coverage using RNA isolated from the supernatant of virus-infected cells. Once optimized we assessed two sources of total RNA, virus-infected cells and supernatant from the virus-infected cells, with four variations of primer pooling for amplicons, and six different amplification approaches. We show that 99–100% genome coverage was obtained using a one-step RT-PCR reaction with one forward and reverse primer. Using a two-step RT-PCR with three distinct tiling approaches for the three genomic segments (vRNAs), we optimized primer pooling approaches for PCR amplification to achieve a greater number of aligned reads, average depth of genome, and genome coverage. The single nucleotide polymorphisms identified from MiSeq and MinION sequencing suggested intrinsic mutation frequencies of ~10⁻⁵-10⁻⁷ per genome and 10⁻⁴-10⁻⁵ per genome, respectively. We noted no difference in the coverage or accuracy when comparing WGS results with amplicons amplified from RNA extracted from infected cells or supernatant of these infected cells. Our results show that high-throughput diagnostics requiring the identification of hantavirus species or strains can be performed using MiSeq or MinION using a one-step approach. However, the two-step MiSeq approach outperformed the MinION in coverage depth and accuracy, and hence would be superior for assessment of genomes for epidemiology or evolutionary questions using the methods developed herein.

Emergence and Magnitude of ML336 Resistance in Venezuelan Equine Encephalitis Virus Depend on the Microenvironment

Venezuelan equine encephalitis virus (VEEV) is a New World Alphavirus that can cause neurological disease and death in humans and equines following transmission from infected mosquitoes. Despite the continued epidemic threat of VEEV, and its potential use as a bioterrorism agent, there are no FDA-approved antivirals or vaccines for treatment or prevention. Previously, we reported the discovery of a small molecule, ML336, with potent antiviral activity against VEEV. To further explore the population-level resistance profiles of ML336, we developed a whole-genome next-generation sequencing (NGS) approach to examine single nucleotide polymorphisms (SNPs) from virus passaged in dose escalation studies in a nonhuman primate kidney epithelial and a human astrocyte cell line, Vero 76 and SVGA, respectively. We passaged VEEV TC-83 in these two cell lines over seven concentrations of ML336, starting at 50 nM. NGS revealed several prominent mutations in the nonstructural protein (nsP) 3 and nsP4 genes that emerged consistently in these two distinct in vitro environments-notably, a mutation at Q210 in nsP4. Several of these mutations were stable following passaging in the absence of ML336 in Vero 76 cells. Network analyses showed that the trajectory of resistance differed between Vero and SVGA. Moreover, the penetration of SNPs was lower in SVGA. In conclusion, we show that the microenvironment influenced the SNP profile of VEEV TC-83. Understanding the dynamics of resistance in VEEV against newly developed antiviral compounds will guide the design of optimal drug candidates and dosing regimens for minimizing the emergence of resistant viruses.IMPORTANCE RNA viruses, including Venezuelan equine encephalitis virus (VEEV), have high mutation rates that allow for rapid adaptation to selective pressures in their environment. Antiviral compounds exert one such pressure on virus populations during infections. Next-generation sequencing allows for examination of viruses at the population level, which enables tracking of low levels of single-nucleotide polymorphisms in the population over time. Therefore, the timing and extent of the emergence of resistance to antivirals can be tracked and assessed. We show here that in VEEV, the trajectory and penetration of antiviral resistance reflected the microenvironment in which the virus population replicates. In summary, we show the diversity of VEEV within a single population under antiviral pressure and two distinct cell types, and we show that population dynamics in these viruses can be examined to better understand how they evolve over time.

Impact of the mutational load on the virological response to a first-line rilpivirine-based regimen

OBJECTIVES

To determine how the load of rilpivirine-resistant variants (mutational load) influences the virological response (VR) of HIV-1-infected patients to a rilpivirine-based first-line regimen.

PATIENTS AND METHODS

Four hundred and eighty-nine patients infected with HIV-1 whose reverse transcriptase gene had been successfully resistance genotyped using next-generation sequencing were given a first-line regimen containing rilpivirine. Variables associated with the VR at 12 months were identified using a logistic model. The results were used to build a multivariate model for each mutational load threshold and the R2 variations were analysed to identify the mutational load threshold that best predicted the VR.

RESULTS

The mutational load at baseline was the only variable linked to the VR at 12 months (P  < 0.01). The VR at 12 months decreased from 96.9% to 83.4% when the mutational load was >1700 copies/mL and to 50% when the mutational load was > 9000 copies/mL. The threshold of 9000 copies/mL was associated with the VR at 12 months with an OR of 36.7 (95% CI 4.7-285.1). The threshold of 1700 copies/mL was associated with the VR at 12 months with an OR of 7.2 (95% CI 1.4-36.8).

CONCLUSIONS

There is quantifiable evidence that determining a mutational load threshold can be used to identify those patients on a first-line regimen containing rilpivirine who are at risk of virological failure. The clinical management of HIV-infected patients can be improved by evaluating the frequency of mutant variants at a threshold of < 20% together with the plasma HIV-1 viral load at the time of resistance genotyping.

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.

External Quality Assessment for Next-Generation Sequencing-Based HIV Drug Resistance Testing: Unique Requirements and Challenges

Over the past decade, there has been an increase in the adoption of next generation sequencing (NGS) technologies for HIV drug resistance (HIVDR) testing. NGS far outweighs conventional Sanger sequencing as it has much higher throughput, lower cost when samples are batched and, most importantly, significantly higher sensitivities for variants present at low frequencies, which may have significant clinical implications. Despite the advantages of NGS, Sanger sequencing remains the gold standard for HIVDR testing, largely due to the lack of standardization of NGS-based HIVDR testing. One important aspect of standardization includes external quality assessment (EQA) strategies and programs. Current EQA for Sanger-based HIVDR testing includes proficiency testing where samples are sent to labs and the performance of the lab conducting such assays is evaluated. The current methods for Sanger-based EQA may not apply to NGS-based tests because of the fundamental differences in their technologies and outputs. Sanger-based genotyping reports drug resistance mutations (DRMs) data as dichotomous, whereas NGS-based HIVDR genotyping also reports DRMs as numerical data (percent abundance). Here we present an overview of the need to develop EQA for NGS-based HIVDR testing and some unique challenges that may be encountered.

Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

The Value of Combined Radial Endobronchial Ultrasound-Guided Transbronchial Lung Biopsy and Metagenomic Next-Generation Sequencing for Peripheral Pulmonary Infectious Lesions

Background

Metagenomic next-generation sequencing (mNGS) is a new technology that allows for unbiased detection of pathogens. However, there are few reports on mNGS of lung biopsy tissues for pulmonary infection diagnosis. In addition, radial endobronchial ultrasound (R-EBUS) is widely used to detect peripheral pulmonary lesions (PPLs), but it is rarely used in the diagnosis of peripheral lung infection.

Objective

The present study aims to evaluate the combined application of R-EBUS-guided transbronchial lung biopsy (TBLB) and mNGS for the diagnosis of peripheral pulmonary infectious lesions.

Methods

From July 2018 to April 2019, 121 patients from Tianjin Medical University General Hospital diagnosed with PPLs and lung infection were enrolled in this prospective randomized study . Once the lesion was located, either TBLB or R-EBUS-guided-TBLB was performed in randomly selected patients, and mNGS was applied for pathogen detection in lung biopsy tissues. The results of mNGS were compared between the TBLB group and R-EBUS-guided TBLB group. In addition, the clinical characteristics and EBUS images from 61 patients receiving bronchoscopy for peripheral lung infectious detection were analyzed and compared with the results of mNGS.

Results

The positivity rate of mNGS in R-EBUS-guided TBLB was (78.7%, 48/61) that was significantly higher than (60.0%, 36/60) in the TBLB group. Difference in the position of R-EBUS probe and image characteristics of peripheral lung infectious lesions affected the positivity rate of mNGS. Tissue collected by R-EBUS within the lesion produced higher positivity rate than samples collected adjacent to the lesion (P=0.030, odds ratio 17.742; 95% confidence interval, from 1.325 to 237.645). Anechoic areas and luminant areas of ultrasonic image characteristics were correlated with lower positivity rate of mNGS (respectively, P=0.019, odds ratio 17.878; 95% confidence interval, from 1.595 to 200.399; P=0.042, odds ratio 16.745; 95% confidence interval, from 1.106 to 253.479).

Conclusions

R-EBUS-guided TBLB is a safe and effective technique in the diagnosis of peripheral lung infectious lesions. R-EBUS significantly facilitates the accurate insertion of bronchoscope into the lesions, which improves positivity rate of mNGS analysis in pathogen detection. The R-EBUS probe position within lesion produced a higher positivity rate of mNGS analysis. Nevertheless, the presence of anechoic and luminant areas on ultrasonic image was correlated with poor mNGS positivity rate.

Evaluation of A Phylogenetic Pipeline to Examine Transmission Networks in A Canadian HIV Cohort

Keywords: HIV; Canada; molecular phylogenetics; viral evolution; person-to-person transmission inference; transmission network; summary statistics.

Multiscale analysis for patterns of Zika virus genotype emergence, spread, and consequence

The question of how Zika virus (ZIKV) changed from a seemingly mild virus to a human pathogen capable of microcephaly and sexual transmission remains unanswered. The unexpected emergence of ZIKV's pathogenicity and capacity for sexual transmission may be due to genetic changes, and future changes in phenotype may continue to occur as the virus expands its geographic range. Alternatively, the sheer size of the 2015-16 epidemic may have brought attention to a pre-existing virulent ZIKV phenotype in a highly susceptible population. Thus, it is important to identify patterns of genetic change that may yield a better understanding of ZIKV emergence and evolution. However, because ZIKV has an RNA genome and a polymerase incapable of proofreading, it undergoes rapid mutation which makes it difficult to identify combinations of mutations associated with viral emergence. As next generation sequencing technology has allowed whole genome consensus and variant sequence data to be generated for numerous virus samples, the task of analyzing these genomes for patterns of mutation has become more complex. However, understanding which combinations of mutations spread widely and become established in new geographic regions versus those that disappear relatively quickly is essential for defining the trajectory of an ongoing epidemic. In this study, multiscale analysis of the wealth of genomic data generated over the course of the epidemic combined with in vivo laboratory data allowed trends in mutations and outbreak trajectory to be assessed. Mutations were detected throughout the genome via deep sequencing, and many variants appeared in multiple samples and in some cases become consensus. Similarly, amino acids that were previously consensus in pre-outbreak samples were detected as low frequency variants in epidemic strains. Protein structural models indicate that most of the mutations associated with the epidemic transmission occur on the exposed surface of viral proteins. At the macroscale level, consensus data was organized into large and interactive databases to allow the spread of individual mutations and combinations of mutations to be visualized and assessed for temporal and geographical patterns. Thus, the use of multiscale modeling for identifying mutations or combinations of mutations that impact epidemic transmission and phenotypic impact can aid the formation of hypotheses which can then be tested using reverse genetics.

Impact of Human Immunodeficiency Virus Type 1 Minority Variants on the Virus Response to a Rilpivirine-Based First-line Regimen

Background

Minority resistant variants of human immunodeficiency virus type 1 (HIV-1) could influence the virological response to treatment based on nonnucleoside reverse transcriptase inhibitors (NNRTIs). Data on minority rilpivirine-resistant variants are scarce. This study used next-generation sequencing (NGS) to identify patients harboring minority resistant variants to nucleos(t)ide reverse transcriptase inhibitors and NNRTIs and to assess their influence on the virological response (VR).

Methods

All the subjects, 541 HIV-1-infected patients started a first-line regimen containing rilpivirine. VR was defined as a HIV-1 RNA load <50 copies/mL at month 6 with continued suppression at month 12. NGS was performed at baseline (retrospectively) on the 454 GS-FLX platform (Roche).

Results

NGS revealed resistance-associated mutations accounting for 1% to <5% of variants in 17.2% of samples, for 5%-20% in 5.7% of samples, and for >20% in 29% of samples. We identified 43 (8.8%) and 36 (7.4%) patients who harbored rilpivirine-resistant variants with a 1% sensitivity threshold according to the French National Agency for Research on AIDS and Viral Hepatitis and Stanford algorithms, respectively. The VR was 96.9% at month 12. Detection of minority rilpivirine resistant variants was not associated with virological failure (VF). Multivariate analysis indicated that VF at month 12 was associated with a CD4 count <250 cells/µL at baseline, a slower decrease in viral load at month 3, and rilpivirine resistance at baseline using the Stanford algorithm with a 20% threshold.

Conclusions

Minority resistant variants had no impact on the VR of treatment-naive patients to a rilpivirine-based regimen.

Caution is needed in interpreting HIV transmission chains by ultradeep sequencing

OBJECTIVES

Molecular epidemiology is applied to various aspects of HIV transmission analyses. With ultradeep sequencing (UDS), in-depth characterization of transmission episodes involving minority variants is permitted. We explored HIV-1 epidemiological linkage and evaluated characteristics of transmission dynamics and transmitted drug resistance (TDR) detection through the added value of UDS.

DESIGN

HIV pol gene fragments were sequenced by UDS and Sanger sequencing on samples of 70 HIV-1-infected, treatment-naive recently diagnosed MSM.

METHODS

Pairwise genetic distances and maximum likelihood phylogenies were computed. Transmission events were identified as clades with branch support at least 70% and intraclade genetic difference less than 4.5%. TDR mutations were recognized from the TDR consensus list. Transmission directionality, directness and inoculum size were inferred from tree topologies.

RESULTS

Both datasets concurred in the identification of seven transmission pairs and one cluster of three patients. With UDS, direction of transmission was inferred in four out of eight chains. Evidence for multiple founder viruses was found in two out of eight chains. No transmission of minority-resistant variants was evidenced. TDR mutations prevalence in protease and reverse transcriptase fragments was 4.3% with Sanger sequencing and 18.6% with UDS.

CONCLUSION

Although Sanger sequencing and UDS identified the same transmission chains, UDS provided additional information on founder viruses, direction of transmission and levels of TDR. Nevertheless, topology of clusters was not always consistent across gene fragments, calling for a cautious interpretation of the data. Moreover, unobserved intermediary links cannot be excluded. Phylogenetic analysis use as a forensic technique for HIV transmission investigations is risky.

Laboratory Methods in Molecular Epidemiology: Viral Infections

Viruses, which are the most abundant biological entities on the planet, have been regarded as the "dark matter" of biology in the sense that despite their ubiquity and frequent presence in large numbers, their detection and analysis are not always straightforward. The majority of them are very small (falling under the limit of 0.5 μm), and collectively, they are extraordinarily diverse. In fact, the majority of the genetic diversity on the planet is found in the so-called virosphere, or the world of viruses. Furthermore, the most frequent viral agents of disease in humans display an RNA genome, and frequently evolve very fast, due to the fact that most of their polymerases are devoid of proofreading activity. Therefore, their detection, genetic characterization, and epidemiological surveillance are rather challenging. This review (part of the Curated Collection on Advances in Molecular Epidemiology of Infectious Diseases) describes many of the methods that, throughout the last few decades, have been used for viral detection and analysis. Despite the challenge of having to deal with high genetic diversity, the majority of these methods still depend on the amplification of viral genomic sequences, using sequence-specific or sequence-independent approaches, exploring thermal profiles or a single nucleic acid amplification temperature. Furthermore, viral populations, and especially those with RNA genomes, are not usually genetically uniform but encompass swarms of genetically related, though distinct, viral genomes known as viral quasispecies. Therefore, sequence analysis of viral amplicons needs to take this fact into consideration, as it constitutes a potential analytic problem. Possible technical approaches to deal with it are also described here. *This article is part of a curated collection.

Advances in Molecular Diagnostic Approaches for Biothreat Agents

The advancement in Molecular techniques has been implicated in the development of sophisticated, high-end diagnostic platform and point-of-care (POC) devices for the detection of biothreat agents. Different molecular and immunological approaches such as Immunochromatographic and lateral flow assays, Enzyme-linked Immunosorbent assays (ELISA), Biosensors, Isothermal amplification assays, Nucleic acid amplification tests (NAATs), Next Generation Sequencers (NGS), Microarrays and Microfluidics have been used for a long time as detection strategies of the biothreat agents. In addition, several point of care (POC) devices have been approved by FDA and commercialized in markets. The high-end molecular platforms like NGS and Microarray are time-consuming, costly, and produce huge amount of data. Therefore, the future prospects of molecular based technique should focus on developing quick, user-friendly, cost-effective and portable devices against biological attacks and surveillance programs.

Characterization of minority HIV-1 drug resistant variants in the United Kingdom following the verification of a deep sequencing-based HIV-1 genotyping and tropism assay

BACKGROUND

The widespread global access to antiretroviral drugs has led to considerable reductions in morbidity and mortality but, unfortunately, the risk of virologic failure increases with the emergence, and potential transmission, of drug resistant viruses. Detecting and quantifying HIV-1 drug resistance has therefore become the standard of care when designing new antiretroviral regimens. The sensitivity of Sanger sequencing-based HIV-1 genotypic assays is limited by its inability to identify minority members of the quasispecies, i.e., it only detects variants present above ~ 20% of the viral population, thus, failing to detect minority variants below this threshold. It is clear that deep sequencing-based HIV-1 genotyping assays are an important step change towards accurately monitoring HIV-infected individuals.

METHODS

We implemented and verified a clinically validated HIV-1 genotyping assay based on deep sequencing (DEEPGEN™) in two clinical laboratories in the United Kingdom: St. George's University Hospitals Healthcare NHS Foundation Trust (London) and at NHS Lothian (Edinburgh), to characterize minority HIV-1 variants in 109 plasma samples from ART-naïve or -experienced individuals.

RESULTS

Although subtype B HIV-1 strains were highly prevalent (44%, 48/109), most individuals were infected with non-B subtype viruses (i.e., A1, A2, C, D, F1, G, CRF02_AG, and CRF01_AE). DEEPGEN™ was able to accurately detect drug resistance-associated mutations not identified using standard Sanger sequencing-based tests, which correlated significantly with patient's antiretroviral treatment histories. A higher proportion of minority PI-, NRTI-, and NNRTI-resistance mutations was detected in NHS Lothian patients compared to individuals from St. George's, mainly M46I/L and I50 V (associated with PIs), D67 N, K65R, L74I, M184 V/I, and K219Q (NRTIs), and L100I (NNRTIs). Interestingly, we observed an inverse correlation between intra-patient HIV-1 diversity and CD4⁺ T cell counts in the NHS Lothian patients.

CONCLUSIONS

This is the first study evaluating the transition, training, and implementation of DEEPGEN™ between three clinical laboratories in two different countries. More importantly, we were able to characterize the HIV-1 drug resistance profile (including minority variants), coreceptor tropism, subtyping, and intra-patient viral diversity in patients from the United Kingdom, providing a rigorous foundation for basing clinical decisions on highly sensitive and cost-effective deep sequencing-based HIV-1 genotyping assays in the country.

Clinical Metagenomic Next-Generation Sequencing for Pathogen Detection

Nearly all infectious agents contain DNA or RNA genomes, making sequencing an attractive approach for pathogen detection. The cost of high-throughput or next-generation sequencing has been reduced by several orders of magnitude since its advent in 2004, and it has emerged as an enabling technological platform for the detection and taxonomic characterization of microorganisms in clinical samples from patients. This review focuses on the application of untargeted metagenomic next-generation sequencing to the clinical diagnosis of infectious diseases, particularly in areas in which conventional diagnostic approaches have limitations. The review covers ( a) next-generation sequencing technologies and common platforms, ( b) next-generation sequencing assay workflows in the clinical microbiology laboratory, ( c) bioinformatics analysis of metagenomic next-generation sequencing data, ( d) validation and use of metagenomic next-generation sequencing for diagnosing infectious diseases, and ( e) significant case reports and studies in this area. Next-generation sequencing is a new technology that has the promise to enhance our ability to diagnose, interrogate, and track infectious diseases.

A novel framework for inferring parameters of transmission from viral sequence data

Transmission between hosts is a critical part of the viral lifecycle. Recent studies of viral transmission have used genome sequence data to evaluate the number of particles transmitted between hosts, and the role of selection as it operates during the transmission process. However, the interpretation of sequence data describing transmission events is a challenging task. We here present a novel and comprehensive framework for using short-read sequence data to understand viral transmission events, designed for influenza virus, but adaptable to other viral species. Our approach solves multiple shortcomings of previous methods for this purpose; for example, we consider transmission as an event involving whole viruses, rather than sets of independent alleles. We demonstrate how selection during transmission and noisy sequence data may each affect naive inferences of the population bottleneck, accounting for these in our framework so as to achieve a correct inference. We identify circumstances in which selection for increased viral transmission may or may not be identified from data. Applying our method to experimental data in which transmission occurs in the presence of strong selection, we show that our framework grants a more quantitative insight into transmission events than previous approaches, inferring the bottleneck in a manner that accounts for selection, both for within-host virulence, and for inherent viral transmissibility. Our work provides new opportunities for studying transmission processes in influenza, and by extension, in other infectious diseases.

The dawn of precision medicine in HIV: state of the art of pharmacotherapy

INTRODUCTION

Combination antiretroviral therapy (ART) reduces viral load to under the limit of detection, successfully decreasing HIV-related morbidity and mortality. Due to viral mutations, complex drug combinations and different patient response, there is an increasing demand for individualized treatment options for patients.

AREAS COVERED

This review first summarizes the pharmacokinetic and pharmacodynamic profile of clinical first-line drugs, which serves as guidance for antiretroviral precision medicine. Factors which have influential effects on drug efficacy and thus precision medicine are discussed: patients' pharmacogenetic information, virus mutations, comorbidities, and immune recovery. Furthermore, strategies to improve the application of precision medicine are discussed.

EXPERT OPINION

Precision medicine for ART requires comprehensive information on the drug, virus, and clinical data from the patients. The clinically available genetic tests are a good starting point. To better apply precision medicine, deeper knowledge of drug concentrations, HIV reservoirs, and efficacy associated genes, such as polymorphisms of drug transporters and metabolizing enzymes, are required. With advanced computer-based prediction systems which integrate more comprehensive information on pharmacokinetics, pharmacodynamics, pharmacogenomics, and the clinically relevant information of the patients, precision medicine will lead to better treatment choices and improved disease outcomes.

Evaluating the accuracy and sensitivity of detecting minority HIV-1 populations by Illumina next-generation sequencing

The accuracy and sensitivity of deep sequencing were assessed using viral standards (pNL4-3 and pLAI.2) of both DNA and RNA. The sequencing accuracy did not depend on the type of nucleic acid, but critically depended on the number of reads and threshold of sensitivity to minor viral populations. With coverage of more than 236 reads, the accuracy of viral RNA sequencing was equal to or exceeded 99.9%, with a sensitivity threshold to minor nucleotides of 20%. When the sensitivity threshold was below 1%, reduced accuracy dynamics were clearly visible even when the coverage was massive (more than 9.000 reads). It was found that the floating sensitivity threshold allowed the sequencing accuracy to be maintained at an acceptable level in cases of low coverage (less than 1.500-2.000) of reads. These results indicate the quality that can be expected with a specific number of reads and sensitivity threshold. Deep sequencing is a very powerful tool that can significantly improve the value of study results, but despite its superior performance, it should be used with caution regarding its sensitivity to minor populations below 1%.

Evaluation of different analysis pipelines for the detection of HIV-1 minority resistant variants

OBJECTIVE

Reliable detection of HIV minority resistant variants (MRVs) requires bioinformatics analysis with specific algorithms to obtain good quality alignments. The aim of this study was to analyze ultra-deep sequencing (UDS) data using different analysis pipelines.

METHODS

HIV-1 protease, reverse transcriptase (RT) and integrase sequences from antiretroviral-naïve patients were obtained using GS-Junior® (Roche) and MiSeq® (Illumina) platforms. MRVs were defined as variants harbouring resistance-mutation present at a frequency of 1%-20%. Reads were analyzed using different alignment algorithms: Amplicon Variant Analyzer®, Geneious® compared to SmartGene® NGS HIV-1 module.

RESULTS

101 protease and 51 RT MRVs identified in 139 protease and 124 RT sequences generated with a GS-Junior® platform were analyzed using AVA® and SmartGene® software. The correlation coefficients for the MRVs were R2 = 0.974 for protease and R2 = 0.972 for RT. Discordances (n = 13 in protease and n = 15 in RT) mainly concerned low-level MRVs (i.e., with frequencies of 1%-2%, n = 18/28) and they were located in homopolymeric regions (n = 10/15). Geneious® and SmartGene® software were used to analyze 143 protease, 45 RT and 26 integrase MRVs identified in 172 protease, 69 RT, and 72 integrase sequences generated with a MiSeq® platform. The correlation coefficients for the MRVs were R2 = 0.987 for protease, R2 = 0.995 for RT and R2 = 0.993 for integrase. Discordances (n = 9 in protease, n = 3 in RT, and n = 3 in integrase) mainly concerned low-level MRVs (n = 13/15).

CONCLUSION

We found an excellent correlation between the various UDS analysis pipelines that we tested. However, our results indicate that specific attention should be paid to low-level MRVs, for which the use of two different analysis pipelines and visual inspection of sequences alignments might be beneficial. Thus, our results argue for use of a 2% threshold for MRV detection, rather than the 1% threshold, to minimize misalignments and time-consuming sight reading steps essential to ensure accurate results for MRV frequencies below 2%.

Detection of minority drug resistant mutations in Malawian HIV-1 subtype C-positive patients initiating and on first-line antiretroviral therapy

Background

Minority drug resistance mutations (DRMs) that are often missed by Sanger sequencing are clinically significant, as they can cause virologic failure in individuals treated with antiretroviral therapy (ART) drugs.

Objective

This study aimed to estimate the prevalence of minor DRMs among patients enrolled in a Malawi HIV drug resistance monitoring survey at baseline and at one year after initiation of ART.

Methods

Forty-one plasma specimens collected from HIV-1 subtype C-positive patients and seven clonal control samples were analysed using ultra-deep sequencing technology.

Results

Deep sequencing identified all 72 DRMs detected by Sanger sequencing at the level of ≥20% and 79 additional minority DRMs at the level of < 20% from the 41 Malawian clinical specimens. Overall, DRMs were detected in 85% of pre-ART and 90.5% of virologic failure patients by deep sequencing. Among pre-ART patients, deep sequencing identified a statistically significant higher prevalence of DRMs to nucleoside reverse transcriptase inhibitors (NRTIs) compared with Sanger sequencing. The difference was mainly due to the high prevalence of minority K65R and M184I mutations. Most virologic failure patients harboured DRMs against both NRTIs and non-nucleoside reverse transcriptase inhibitors (NNRTIs). These minority DRMs contributed to the increased or enhanced virologic failures in these patients.

Conclusion

The results revealed the presence of minority DRMs to NRTIs and NNRTIs in specimens collected at baseline and virologic failure time points. These minority DRMs not only increased resistance levels to NRTIs and NNRTIs for the prescribed ART, but also expanded resistance to additional major first-line ART drugs. This study suggested that drug resistance testing that uses more sensitive technologies, is needed in this setting.

Added Value of Next-Generation Sequencing for Multilocus Sequence Typing Analysis of a Pneumocystis jirovecii Pneumonia Outbreak1

Pneumocystis jirovecii is a major threat for immunocompromised patients, and clusters of pneumocystis pneumonia (PCP) have been increasingly described in transplant units during the past decade. Exploring an outbreak transmission network requires complementary spatiotemporal and strain-typing approaches. We analyzed a PCP outbreak and demonstrated the added value of next-generation sequencing (NGS) for the multilocus sequence typing (MLST) study of P. jirovecii strains. Thirty-two PCP patients were included. Among the 12 solid organ transplant patients, 5 shared a major and unique genotype that was also found as a minor strain in a sixth patient. A transmission map analysis strengthened the suspicion of nosocomial acquisition of this strain for the 6 patients. NGS-MLST enables accurate determination of subpopulation, which allowed excluding other patients from the transmission network. NGS-MLST genotyping approach was essential to deciphering this outbreak. This innovative approach brings new insights for future epidemiologic studies on this uncultivable opportunistic fungus.

Quasispecies characters of hepatitis B virus in immunoprophylaxis failure infants

Hepatitis B vaccination prevents 80-95% of transmission and reduces the incidence of HBV in children. The variations in the a determinant of HBV surface antigen (HBsAg) have been reported to be the most prevalent cause for vaccine or antibody escape. There is a conflicting evidence on as to whether escape mutants arise de novo in infected infants or whether the mutants, that have preexisted maternally, subsequently undergo selective replication in the infant under immune pressure. Here, we report that nearly 65% (55 of 85) vaccination failure in child patients has no amino acid substitution in a determinant as seen by Sanger sequencing. We further employed an Illumina sequencing platform-based method to detect HBV quasispecies in four immunoprophylaxis failure infants and their mothers. In our data, the substitution rate of amino acid located at a determinant is relatively low (< 10%), I/T126A, C124S, F134Y, K141Q, Q129H, D144A, G145V, and N146K, which showed no statistical difference to their mothers, proving that these vaccine escape mutants preexist maternally as minor variants. Besides that, bioinformatical analysis showed that the binding affinity of high variation epitopes (amino acid divergence in mother and their infants > 20%) to related HLA molecules was generally decreased, these traces of immune escape suggesting that immune pressure was present and was effective in all samples.

In-depth genome analyses of viruses from vaccine-derived rabies cases and corresponding live-attenuated oral rabies vaccines

Live-attenuated rabies virus strains such as those derived from the field isolate Street Alabama Dufferin (SAD) have been used extensively and very effectively as oral rabies vaccines for the control of fox rabies in both Europe and Canada. Although these vaccines are safe, some cases of vaccine-derived rabies have been detected during rabies surveillance accompanying these campaigns. In recent analysis it was shown that some commercial SAD vaccines consist of diverse viral populations, rather than clonal genotypes. For cases of vaccine-derived rabies, only consensus sequence data have been available to date and information concerning their population diversity was thus lacking. In our study, we used high-throughput sequencing to analyze 11 cases of vaccine-derived rabies, and compared their viral population diversity to the related oral rabies vaccines using pairwise Manhattan distances. This extensive deep sequencing analysis of vaccine-derived rabies cases observed during oral vaccination programs provided deeper insights into the effect of accidental in vivo replication of genetically diverse vaccine strains in the central nervous system of target and non-target species under field conditions. The viral population in vaccine-derived cases appeared to be clonal in contrast to their parental vaccines. The change from a state of high population diversity present in the vaccine batches to a clonal genotype in the affected animal may indicate the presence of a strong bottleneck during infection. In conclusion, it is very likely that these few cases are the consequence of host factors and not the result of the selection of a more virulent genotype. Furthermore, this type of vaccine-derived rabies leads to the selection of clonal genotypes and the selected variants were genetically very similar to potent SAD vaccines that have undergone a history of in vitro selection.

On the effective depth of viral sequence data

Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

No Substantial Evidence for Sexual Transmission of Minority HIV Drug Resistance Mutations in Men Who Have Sex with Men

During primary HIV infection, the presence of minority drug resistance mutations (DRM) may be a consequence of sexual transmission, de novo mutations, or technical errors in identification. Baseline blood samples were collected from 24 HIV-infected antiretroviral-naive, genetically and epidemiologically linked source and recipient partners shortly after the recipient's estimated date of infection. An additional 32 longitudinal samples were available from 11 recipients. Deep sequencing of HIV reverse transcriptase (RT) was performed (Roche/454), and the sequences were screened for nucleoside and nonnucleoside RT inhibitor DRM. The likelihood of sexual transmission and persistence of DRM was assessed using Bayesian-based statistical modeling. While the majority of DRM (>20%) were consistently transmitted from source to recipient, the probability of detecting a minority DRM in the recipient was not increased when the same minority DRM was detected in the source (Bayes factor [BF] = 6.37). Longitudinal analyses revealed an exponential decay of DRM (BF = 0.05) while genetic diversity increased. Our analysis revealed no substantial evidence for sexual transmission of minority DRM (BF = 0.02). The presence of minority DRM during early infection, followed by a rapid decay, is consistent with the "mutation-selection balance" hypothesis, in which deleterious mutations are more efficiently purged later during HIV infection when the larger effective population size allows more efficient selection. Future studies using more recent sequencing technologies that are less prone to single-base errors should confirm these results by applying a similar Bayesian framework in other clinical settings.IMPORTANCE The advent of sensitive sequencing platforms has led to an increased identification of minority drug resistance mutations (DRM), including among antiretroviral therapy-naive HIV-infected individuals. While transmission of DRM may impact future therapy options for newly infected individuals, the clinical significance of the detection of minority DRM remains controversial. In the present study, we applied deep-sequencing techniques within a Bayesian hierarchical framework to a cohort of 24 transmission pairs to investigate whether minority DRM detected shortly after transmission were the consequence of (i) sexual transmission from the source, (ii) de novo emergence shortly after infection followed by viral selection and evolution, or (iii) technical errors/limitations of deep-sequencing methods. We found no clear evidence to support the sexual transmission of minority resistant variants, and our results suggested that minor resistant variants may emerge de novo shortly after transmission, when the small effective population size limits efficient purge by natural selection.

High prevalence of HIV-1 transmitted drug-resistance mutations from proviral DNA massively parallel sequencing data of therapy-naïve chronically infected Brazilian blood donors

BACKGROUND

An improved understanding of the prevalence of low-abundance transmitted drug-resistance mutations (TDRM) in therapy-naïve HIV-1-infected patients may help determine which patients are the best candidates for therapy. In this study, we aimed to obtain a comprehensive picture of the evolving HIV-1 TDRM across the massive parallel sequences (MPS) of the viral entire proviral genome in a well-characterized Brazilian blood donor naïve to antiretroviral drugs.

MATERIALS AND METHODS

The MPS data from 128 samples used in the analysis were sourced from Brazilian blood donors and were previously classified by less-sensitive (LS) or "detuned" enzyme immunoassay as non-recent or longstanding HIV-1 infections. The Stanford HIV Resistance Database (HIVDBv 6.2) and IAS-USA mutation lists were used to interpret the pattern of drug resistance. The minority variants with TDRM were identified using a threshold of ≥ 1.0% and ≤ 20% of the reads sequenced. The rate of TDRM in the MPS data of the proviral genome were compared with the corresponding published consensus sequences of their plasma viruses.

RESULTS

No TDRM were detected in the integrase or envelope regions. The overall prevalence of TDRM in the protease (PR) and reverse transcriptase (RT) regions of the HIV-1 pol gene was 44.5% (57/128), including any mutations to the nucleoside analogue reverse transcriptase inhibitors (NRTI) and non-nucleoside analogue reverse transcriptase inhibitors (NNRTI). Of the 57 subjects, 43 (75.4%) harbored a minority variant containing at least one clinically relevant TDRM. Among the 43 subjects, 33 (76.7%) had detectable minority resistant variants to NRTIs, 6 (13.9%) to NNRTIs, and 16 (37.2%) to PR inhibitors. The comparison of viral sequences in both sources, plasma and cells, would have detected 48 DNA provirus disclosed TDRM by MPS previously missed by plasma bulk analysis.

CONCLUSION

Our findings revealed a high prevalence of TDRM found in this group, as the use of MPS drastically increased the detection of these mutations. Sequencing proviral DNA provided additional information about TDRM, which may impact treatment decisions. The overall results emphasize the importance of continuous monitoring.

Generic Amplicon Deep Sequencing to Determine Ilarvirus Species Diversity in Australian Prunus

The distribution of Ilarvirus species populations amongst 61 Australian Prunus trees was determined by next generation sequencing (NGS) of amplicons generated using a genus-based generic RT-PCR targeting a conserved region of the Ilarvirus RNA2 component that encodes the RNA dependent RNA polymerase (RdRp) gene. Presence of Ilarvirus sequences in each positive sample was further validated by Sanger sequencing of cloned amplicons of regions of each of RNA1, RNA2 and/or RNA3 that were generated by species specific PCRs and by metagenomic NGS. Prunus necrotic ringspot virus (PNRSV) was the most frequently detected Ilarvirus, occurring in 48 of the 61 Ilarvirus-positive trees and Prune dwarf virus (PDV) and Apple mosaic virus (ApMV) were detected in three trees and one tree, respectively. American plum line pattern virus (APLPV) was detected in three trees and represents the first report of APLPV detection in Australia. Two novel and distinct groups of Ilarvirus-like RNA2 amplicon sequences were also identified in several trees by the generic amplicon NGS approach. The high read depth from the amplicon NGS of the generic PCR products allowed the detection of distinct RNA2 RdRp sequence variant populations of PNRSV, PDV, ApMV, APLPV and the two novel Ilarvirus-like sequences. Mixed infections of ilarviruses were also detected in seven Prunus trees. Sanger sequencing of specific RNA1, RNA2, and/or RNA3 genome segments of each virus and total nucleic acid metagenomics NGS confirmed the presence of PNRSV, PDV, ApMV and APLPV detected by RNA2 generic amplicon NGS. However, the two novel groups of Ilarvirus-like RNA2 amplicon sequences detected by the generic amplicon NGS could not be associated to the presence of sequence from RNA1 or RNA3 genome segments or full Ilarvirus genomes, and their origin is unclear. This work highlights the sensitivity of genus-specific amplicon NGS in detection of virus sequences and their distinct populations in multiple samples, and the need for a standardized approach to accurately determine what constitutes an active, viable virus infection after detection by molecular based methods.

Analysis of intra-host genetic diversity of Prunus necrotic ringspot virus (PNRSV) using amplicon next generation sequencing

PCR amplicon next generation sequencing (NGS) analysis offers a broadly applicable and targeted approach to detect populations of both high- or low-frequency virus variants in one or more plant samples. In this study, amplicon NGS was used to explore the diversity of the tripartite genome virus, Prunus necrotic ringspot virus (PNRSV) from 53 PNRSV-infected trees using amplicons from conserved gene regions of each of PNRSV RNA1, RNA2 and RNA3. Sequencing of the amplicons from 53 PNRSV-infected trees revealed differing levels of polymorphism across the three different components of the PNRSV genome with a total number of 5040, 2083 and 5486 sequence variants observed for RNA1, RNA2 and RNA3 respectively. The RNA2 had the lowest diversity of sequences compared to RNA1 and RNA3, reflecting the lack of flexibility tolerated by the replicase gene that is encoded by this RNA component. Distinct PNRSV phylo-groups, consisting of closely related clusters of sequence variants, were observed in each of PNRSV RNA1, RNA2 and RNA3. Most plant samples had a single phylo-group for each RNA component. Haplotype network analysis showed that smaller clusters of PNRSV sequence variants were genetically connected to the largest sequence variant cluster within a phylo-group of each RNA component. Some plant samples had sequence variants occurring in multiple PNRSV phylo-groups in at least one of each RNA and these phylo-groups formed distinct clades that represent PNRSV genetic strains. Variants within the same phylo-group of each Prunus plant sample had ≥97% similarity and phylo-groups within a Prunus plant sample and between samples had less ≤97% similarity. Based on the analysis of diversity, a definition of a PNRSV genetic strain was proposed. The proposed definition was applied to determine the number of PNRSV genetic strains in each of the plant samples and the complexity in defining genetic strains in multipartite genome viruses was explored.

Multi-drug resistant Klebsiella pneumoniae strains circulating in hospital setting: whole-genome sequencing and Bayesian phylogenetic analysis for outbreak investigations

Carbapenems resistant Enterobacteriaceae infections are increasing worldwide representing an emerging public health problem. The application of phylogenetic and phylodynamic analyses to bacterial whole genome sequencing (WGS) data have become essential in the epidemiological surveillance of multi-drug resistant nosocomial pathogens. Between January 2012 and February 2013, twenty-one multi-drug resistant K. pneumoniae strains, were collected from patients hospitalized among different wards of the University Hospital Campus Bio-Medico. Epidemiological contact tracing of patients and Bayesian phylogenetic analysis of bacterial WGS data were used to investigate the evolution and spatial dispersion of K. pneumoniae in support of hospital infection control. The epidemic curve of incident K. pneumoniae cases showed a bimodal distribution of cases with two peaks separated by 46 days between November 2012 and January 2013. The time-scaled phylogeny suggested that K. pneumoniae strains isolated during the study period may have been introduced into the hospital setting as early as 2007. Moreover, the phylogeny showed two different epidemic introductions in 2008 and 2009. Bayesian genomic epidemiology is a powerful tool that promises to improve the surveillance and control of multi-drug resistant pathogens in an effort to develop effective infection prevention in healthcare settings or constant strains reintroduction.

Detection of Emerging Vaccine-Related Polioviruses by Deep Sequencing

Oral poliovirus vaccine can mutate to regain neurovirulence. To date, evaluation of these mutations has been performed primarily on culture-enriched isolates by using conventional Sanger sequencing. We therefore developed a culture-independent, deep-sequencing method targeting the 5′ untranslated region (UTR) and P1 genomic region to characterize vaccine-related poliovirus variants. Error analysis of the deep-sequencing method demonstrated reliable detection of poliovirus mutations at levels of <1%, depending on read depth. Sequencing of viral nucleic acids from the stool of vaccinated, asymptomatic children and their close contacts collected during a prospective cohort study in Veracruz, Mexico, revealed no vaccine-derived polioviruses. This was expected given that the longest duration between sequenced sample collection and the end of the most recent national immunization week was 66 days. However, we identified many low-level variants (<5%) distributed across the 5′ UTR and P1 genomic region in all three Sabin serotypes, as well as vaccine-related viruses with multiple canonical mutations associated with phenotypic reversion present at high levels (>90%). These results suggest that monitoring emerging vaccine-related poliovirus variants by deep sequencing may aid in the poliovirus endgame and efforts to ensure global polio eradication.

Disparities in HIV-1 transmitted drug resistance detected by ultradeep sequencing between men who have sex with men and heterosexual populations

OBJECTIVES

Transmitted drug resistance (TDR) can impair the response to first-line antiretroviral therapy. In treatment-naïve patients chronically infected with HIV type 1 (HIV-1), it was previously shown through Sanger sequencing that TDR was more common in men who have sex with men (MSM) than in other transmission risk groups. We aimed to compare two HIV-1 transmission groups in terms of the presence of TDR mutations.

METHODS

We investigated, through Sanger sequencing and ultradeep sequencing (UDS), the presence of resistance mutations, both in majority (> 20%) and in minority (1-20%) proportions, in 70 treatment-naïve MSM and 70 treatment-naïve heterosexual patients who recently screened positive for HIV-1.

RESULTS

The global prevalence of TDR was not significantly different between the two groups, either by Sanger or by UDS. Nevertheless, a higher frequency of nucleoside reverse transcriptase inhibitor TDR was observed among heterosexual patients (P = 0.04). There was also a trend for a higher frequency of TDR among MSM infected with HIV-1 subtype B compared with MSM infected with HIV-1 non-B subtypes (P = 0.06).

CONCLUSIONS

Ultradeep sequencing UDS allowed sensitive monitoring of TDR, and highlighted some disparities between transmission groups.

Mutations during the Adaptation of H9N2 Avian Influenza Virus to the Respiratory Epithelium of Pigs Enhance Sialic Acid Binding Activity and Virulence in Mice

The natural reservoir for influenza viruses is waterfowl, and from there they succeeded in crossing the barrier to different mammalian species. We analyzed the adaptation of avian influenza viruses to a mammalian host by passaging an H9N2 strain three times in differentiated swine airway epithelial cells. Using precision-cut slices from the porcine lung to passage the parental virus, isolates from each of the three passages (P1 to P3) were characterized by assessing growth curves and ciliostatic effects. The only difference noted was an increased growth kinetics of the P3 virus. Sequence analysis revealed four mutations: one each in the PB2 and NS1 proteins and two in the HA protein. The HA mutations, A190V and T212I, were characterized by generating recombinant viruses containing either one or both amino acid exchanges. Whereas the parental virus recognized α2,3-linked sialic acids preferentially, the HA190 mutant bound to a broad spectrum of glycans with α2,6/8/9-linked sialic acids. The HA212 mutant alone differed only slightly from the parental virus; however, the combination of both mutations (HA190+HA212) increased the binding affinity to those glycans recognized by the HA190 mutant. Remarkably, only the HA double mutant showed a significantly increased pathogenicity in mice. In contrast, none of those mutations affected the ciliary activity of the epithelial cells which is characteristic for virulent swine influenza viruses. Taken together, our results indicate that shifts in the HA receptor affinity are just an early adaptation step of avian H9N2 strains; further mutational changes may be required to become virulent for pigs.IMPORTANCE Swine play an important role in the interspecies transmission of influenza viruses. Avian influenza A viruses (IAV) of the H9N2 subtype have successfully infected hosts from different species but have not established a stable lineage. We have analyzed the adaptation of IAV-H9N2 virus to target cells of a new host by passaging the virus three times in differentiated porcine respiratory epithelial cells. Among the four mutations detected, the two HA mutations were analyzed by generating recombinant viruses. Depending on the infection system used, the mutations differed in their phenotypic expression, e.g., sialic acid binding activity, replication kinetics, plaque size, and pathogenicity in inbred mice. However, none of the mutations affected the ciliary activity which serves as a virulence marker. Thus, early adaptive mutation enhances the replication kinetics, but more mutations are required for IAV of the H9N2 subtype to become virulent.

Impaired human immunodeficiency virus type 1 replicative fitness in atypical viremic non-progressor individuals

Background

Progression rates from initial HIV-1 infection to advanced AIDS vary significantly among infected individuals. A distinct subgroup of HIV-1-infected individuals—termed viremic non-progressors (VNP) or controllers—do not seem to progress to AIDS, maintaining high CD4⁺ T cell counts despite high levels of viremia for many years. Several studies have evaluated multiple host factors, including immune activation, trying to elucidate the atypical HIV-1 disease progression in these patients; however, limited work has been done to characterize viral factors in viremic controllers.

Methods

We analyzed HIV-1 isolates from three VNP individuals and compared the replicative fitness, near full-length HIV-1 genomes and intra-patient HIV-1 genetic diversity with viruses from three typical (TP) and one rapid (RP) progressor individuals.

Results

Viremic non-progressors and typical patients were infected for >10 years (range 10–17 years), with a mean CD4⁺ T-cell count of 472 cells/mm³ (442–529) and 400 cells/mm³ (126–789), respectively. VNP individuals had a less marked decline in CD4⁺ cells (mean −0.56, range −0.4 to −0.7 CD4⁺/month) than TP patients (mean −10.3, −8.2 to −13.1 CD4⁺/month). Interestingly, VNP individuals carried viruses with impaired replicative fitness, compared to HIV-1 isolates from the TP and RP patients (p < 0.05, 95% CI). Although analyses of the near full-length HIV-1 genomes showed no clear patterns of single-nucleotide polymorphisms (SNP) that could explain the decrease in replicative fitness, both the number of SNPs and HIV-1 population diversity correlated inversely with the replication capacity of the viruses (r = −0.956 and r = −0.878, p < 0.01, respectively).

Conclusion

It is likely that complex multifactorial parameters govern HIV-1 disease progression in each individual, starting with the infecting virus (phenotype, load, and quasispecies diversity) and the intrinsic ability of the host to respond to the infection. Here we analyzed a subset of viremic controller patients and demonstrated that similar to the phenomenon observed in patients with a discordant response to antiretroviral therapy (i.e., high CD4⁺ cell counts with detectable plasma HIV-1 RNA load), reduced viral replicative fitness seems to be linked to slow disease progression in these antiretroviral-naïve individuals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12981-017-0144-0) contains supplementary material, which is available to authorized users.

Performance comparison of next-generation sequencing platforms for determining HIV-1 coreceptor use

The coreceptor used by HIV-1 must be determined before a CCR5 antagonist, part of the arsenal of antiretroviral drugs, is prescribed because viruses that enter cells using the CXCR4 coreceptor are responsible for treatment failure. HIV-1 tropism is also correlated with disease progression and so must be determined for virological studies. Tropism can be determined by next-generation sequencing (NGS), but not all of these new technologies have been fully validated for use in clinical practice. The Illumina NGS technology is used in many laboratories but its ability to predict HIV-1 tropism has not been evaluated while the 454 GS-Junior (Roche) is used for routine diagnosis. The genotypic prediction of HIV-1 tropism is based on sequencing the V3 region and interpreting the results with an appropriate algorithm. We compared the performances of the MiSeq (Illumina) and 454 GS-Junior (Roche) systems with a reference phenotypic assay. We used clinical samples for the NGS tropism predictions and assessed their ability to quantify CXCR4-using variants. The data show that the Illumina platform can be used to detect minor CXCR4-using variants in clinical practice but technical optimization are needed to improve quantification.

Error baseline rates of five sample preparation methods used to characterize RNA virus populations

Individual RNA viruses typically occur as populations of genomes that differ slightly from each other due to mutations introduced by the error-prone viral polymerase. Understanding the variability of RNA virus genome populations is critical for understanding virus evolution because individual mutant genomes may gain evolutionary selective advantages and give rise to dominant subpopulations, possibly even leading to the emergence of viruses resistant to medical countermeasures. Reverse transcription of virus genome populations followed by next-generation sequencing is the only available method to characterize variation for RNA viruses. However, both steps may lead to the introduction of artificial mutations, thereby skewing the data. To better understand how such errors are introduced during sample preparation, we determined and compared error baseline rates of five different sample preparation methods by analyzing in vitro transcribed Ebola virus RNA from an artificial plasmid-based system. These methods included: shotgun sequencing from plasmid DNA or in vitro transcribed RNA as a basic “no amplification” method, amplicon sequencing from the plasmid DNA or in vitro transcribed RNA as a “targeted” amplification method, sequence-independent single-primer amplification (SISPA) as a “random” amplification method, rolling circle reverse transcription sequencing (CirSeq) as an advanced “no amplification” method, and Illumina TruSeq RNA Access as a “targeted” enrichment method. The measured error frequencies indicate that RNA Access offers the best tradeoff between sensitivity and sample preparation error (1.4⁻⁵) of all compared methods.