Within-Host Rhinovirus Evolution in Upper and Lower Respiratory Tract Highlights Capsid Variability and Mutation-Independent Compartmentalization.

BACKGROUND: Rhinovirus (RV) infections can progress from the upper (URT) to lower (LRT) respiratory tract in immunocompromised individuals, causing high rates of fatal pneumonia. Little is known about how RV evolves within hosts during infection. METHODS: We sequenced RV complete genomes from 12 hematopoietic cell transplant patients with infection for up to 190 days from both URT (nasal wash, NW) and LRT (bronchoalveolar lavage, BAL). Metagenomic and amplicon next-generation sequencing were used to track the emergence and evolution of intrahost single nucleotide variants (iSNVs). RESULTS: Identical RV intrahost populations in matched NW and BAL specimens indicated no genetic adaptation is required for RV to progress from URT to LRT. Coding iSNVs were 2.3-fold more prevalent in capsid over nonstructural genes. iSNVs modeled were significantly more likely to be found in capsid surface residues, but were not preferentially located in known RV-neutralizing antibody epitopes. Newly emergent, genotype-matched iSNV haplotypes from immunocompromised individuals in 2008-2010 could be detected in Seattle-area community RV sequences in 2020-2021. CONCLUSIONS: RV infections in immunocompromised hosts can progress from URT to LRT with no specific evolutionary requirement. Capsid proteins carry the highest variability and emergent mutations can be detected in other, including future, RV sequences.

Within-host rhinovirus evolution in upper and lower respiratory tract highlights capsid variability and mutation-independent compartmentalization.

Background: Human rhinovirus (HRV) infections can progress from the upper (URT) to lower (LRT) respiratory tract in immunocompromised individuals, causing high rates of fatal pneumonia. Little is known about how HRV evolves within hosts during infection. Methods: We sequenced HRV complete genomes from 12 hematopoietic cell transplant patients with prolonged infection for up to 190 days from both URT (nasal wash, NW) and LRT (bronchoalveolar lavage, BAL) specimens. Metagenomic (mNGS) and amplicon-based NGS were used to study the emergence and evolution of intra-host single nucleotide variants (iSNVs). Results: Identical HRV intra-host populations in matched NW and BAL specimens indicated no genetic adaptation is required for HRV to progress from URT to LRT. Microbial composition between matched NW and BAL confirmed no cross-contamination during sampling procedure. Coding iSNVs were 2.3-fold more prevalent in capsid over non-structural genes, adjusted for length. iSNVs modeled onto HRV capsid structures were significantly more likely to be found in surface residues, but were not preferentially located in known HRV neutralizing antibody epitopes. Newly emergent, serotype-matched iSNV haplotypes from immunocompromised individuals from 2008-2010 could be detected in Seattle-area community HRV sequences from 2020-2021. Conclusion: HRV infections in immunocompromised hosts can progress from URT to LRT with no specific evolutionary requirement. Capsid proteins carry the highest variability and emergent mutations can be detected in other, including future, HRV sequences.

Human parainfluenza virus evolution during lung infection of immunocompromised individuals promotes viral persistence.

The capacity of respiratory viruses to undergo evolution within the respiratory tract raises the possibility of evolution under the selective pressure of the host environment or drug treatment. Long-term infections in immunocompromised hosts are potential drivers of viral evolution and development of infectious variants. We showed that intrahost evolution in chronic human parainfluenza virus 3 (HPIV3) infection in immunocompromised individuals elicited mutations that favored viral entry and persistence, suggesting that similar processes may operate across enveloped respiratory viruses. We profiled longitudinal HPIV3 infections from 2 immunocompromised individuals that persisted for 278 and 98 days. Mutations accrued in the HPIV3 attachment protein hemagglutinin-neuraminidase (HN), including the first in vivo mutation in HN's receptor binding site responsible for activating the viral fusion process. Fixation of this mutation was associated with exposure to a drug that cleaves host-cell sialic acid moieties. Longitudinal adaptation of HN was associated with features that promote viral entry and persistence in cells, including greater avidity for sialic acid and more active fusion activity in vitro, but not with antibody escape. Long-term infection thus led to mutations promoting viral persistence, suggesting that host-directed therapeutics may support the evolution of viruses that alter their biophysical characteristics to persist in the face of these agents in vivo.

Gene editing and elimination of latent herpes simplex virus in vivo.

We evaluate gene editing of HSV in a well-established mouse model, using adeno-associated virus (AAV)-delivered meganucleases, as a potentially curative approach to treat latent HSV infection. Here we show that AAV-delivered meganucleases, but not CRISPR/Cas9, mediate highly efficient gene editing of HSV, eliminating over 90% of latent virus from superior cervical ganglia. Single-cell RNA sequencing demonstrates that both HSV and individual AAV serotypes are non-randomly distributed among neuronal subsets in ganglia, implying that improved delivery to all neuronal subsets may lead to even more complete elimination of HSV. As predicted, delivery of meganucleases using a triple AAV serotype combination results in the greatest decrease in ganglionic HSV loads. The levels of HSV elimination observed in these studies, if translated to humans, would likely significantly reduce HSV reactivation, shedding, and lesions. Further optimization of meganuclease delivery and activity is likely possible, and may offer a pathway to a cure for HSV infection.

Comparison of Three Adenovirus Quantitative PCR Assays with ATCC Reference Strains and Clinical Samples.

Adenoviruses (AdV) have been associated with a variety of human diseases and are recognized as causing significant morbidity and mortality in immunocompromised or transplant patients. Quantification of AdV DNA in plasma is notoriously difficult due to the genetic diversity of the 71 different serotypes identified to date. There is no World Health Organization standard available to harmonize quantitative data, so results between labs vary widely. In this study, we compared a laboratory-developed multiplex PCR assay with primers and probes specific for each group (A to G) and subgroup E4 (Octaplex) to one with a single primer and probe set (modified from N. Jothikumar et al., Appl Environ Microbiol 71:3131-3136, 2005) and one utilizing bisulfite pretreatment of DNA to reduce variation prior to amplification (Genetic Signatures). Our Octaplex assay detected all low-copy-number clinical samples, while the other two assays had subsets of samples that did not amplify. The modified Jothikumar assay failed to efficiently amplify three of the high-copy-number cultured strains, while the Genetic Signatures 3base assay had a positive bias, resulting in higher copies/ml (>0.5 log10) for all culture fluids tested. All three assays resulted in endpoint detection of the available 51 AdV types. Using two different materials to generate a standard curve revealed that the Octaplex TaqMan assay and the modified Jothikumar assay both consistently gave adenovirus levels lower than the commercial platform for AdV culture fluids but not patient samples. This study highlights the differences in detection of AdV between laboratories that can be attributed to both the PCR method, as well as the reference material used for quantitation.

A Novel, Widespread qacA Allele Results in Reduced Chlorhexidine Susceptibility in Staphylococcus epidermidis.

Chlorhexidine gluconate (CHG) is a topical antiseptic widely used in health care settings. In Staphylococcus spp., the pump QacA effluxes CHG, while the closely related QacB cannot due to a single amino acid substitution. We characterized 1,050 cutaneous Staphylococcus isolates obtained from 173 pediatric oncology patients enrolled in a multicenter CHG bathing trial. CHG susceptibility testing revealed that 63 (6%) of these isolates had elevated CHG MICs (≥4 µg/ml). Screening of all 1,050 isolates for the qacA/B gene (the same qac gene with A or B allele) by restriction fragment length polymorphism (RFLP) yielded 56 isolates with a novel qacA/B RFLP pattern, qacA/B273 The CHG MIC was significantly higher for qacA/B273 -positive isolates (MIC50, 4 µg/ml; MIC range, 0.5 to 4 µg/ml) than for other qac groups: qacA-positive isolates (n = 559; MIC50, 1 µg/ml; MIC range, 0.5 to 4 µg/ml), qacB-positive isolates (n = 17; MIC50, 1 µg/ml; MIC range, 0.25 to 2 µg/ml), and qacA/B-negative isolates (n = 418, MIC50, 1 µg/ml; MIC range, 0.125 to 2 µg/ml) (P = 0.001). A high proportion of the qacA/B273 -positive isolates also displayed methicillin resistance (96.4%) compared to the other qac groups (24.9 to 61.7%) (P = 0.001). Whole-genome sequencing revealed that qacA/B273 -positive isolates encoded a variant of QacA with 2 amino acid substitutions. This new allele, named qacA4, was carried on the novel plasmid pAQZ1. The qacA4-carrying isolates belonged to the highly resistant Staphylococcus epidermidis sequence type 2 clone. By searching available sequence data sets, we identified 39 additional qacA4-carrying S. epidermidis strains from 5 countries. Curing an isolate of qacA4 resulted in a 4-fold decrease in the CHG MIC, confirming the role of qacA4 in the elevated CHG MIC. Our results highlight the importance of further studying qacA4 and its functional role in clinical staphylococci.

VAPiD: a lightweight cross-platform viral annotation pipeline and identification tool to facilitate virus genome submissions to NCBI GenBank.

BACKGROUND: With sequencing technologies becoming cheaper and easier to use, more groups are able to obtain whole genome sequences of viruses of public health and scientific importance. Submission of genomic data to NCBI GenBank is a requirement prior to publication and plays a critical role in making scientific data publicly available. GenBank currently has automatic prokaryotic and eukaryotic genome annotation pipelines but has no viral annotation pipeline beyond influenza virus. Annotation and submission of viral genome sequence is a non-trivial task, especially for groups that do not routinely interact with GenBank for data submissions. RESULTS: We present Viral Annotation Pipeline and iDentification (VAPiD), a portable and lightweight command-line tool for annotation and GenBank deposition of viral genomes. VAPiD supports annotation of nearly all unsegmented viral genomes. The pipeline has been validated on human immunodeficiency virus, human parainfluenza virus 1-4, human metapneumovirus, human coronaviruses (229E/OC43/NL63/HKU1/SARS/MERS), human enteroviruses/rhinoviruses, measles virus, mumps virus, Hepatitis A-E Virus, Chikungunya virus, dengue virus, and West Nile virus, as well the human polyomaviruses BK/JC/MCV, human adenoviruses, and human papillomaviruses. The program can handle individual or batch submissions of different viruses to GenBank and correctly annotates multiple viruses, including those that contain ribosomal slippage or RNA editing without prior knowledge of the virus to be annotated. VAPiD is programmed in Python and is compatible with Windows, Linux, and Mac OS systems. CONCLUSIONS: We have created a portable, lightweight, user-friendly, internet-enabled, open-source, command-line genome annotation and submission package to facilitate virus genome submissions to NCBI GenBank. Instructions for downloading and installing VAPiD can be found at https://github.com/rcs333/VAPiD .

Limited Marginal Utility of Deep Sequencing for HIV Drug Resistance Testing in the Age of Integrase Inhibitors.

HIV drug resistance genotyping is a critical tool in the clinical management of HIV infections. Although resistance genotyping has traditionally been conducted using Sanger sequencing, next-generation sequencing (NGS) is emerging as a powerful tool due to its ability to detect low-frequency alleles. However, the clinical value added from NGS approaches to antiviral resistance testing remains to be demonstrated. We compared the variant detection capacity of NGS versus Sanger sequencing methods for resistance genotyping in 144 drug resistance tests (105 protease-reverse transcriptase tests and 39 integrase tests) submitted to our clinical virology laboratory over a four-month period in 2016 for Sanger-based HIV drug resistance testing. NGS detected all true high-frequency drug resistance mutations (>20% frequency) found by Sanger sequencing, with greater accuracy in one instance of a Sanger-detected false positive. Freely available online NGS variant callers HyDRA and PASeq were superior to Sanger methods for interpretations of allele linkage and automated variant calling. NGS additionally detected low-frequency mutations (1 to 20% frequency) associated with higher levels of drug resistance in 30/105 (29%) protease-reverse transcriptase tests and 4/39 (10%) integrase tests. In clinical follow-up of 69 individuals for a median of 674 days, we did not find a difference in rates of virological failure between individuals with and without low-frequency mutations, although rates of virological failure were higher for individuals with drug-relevant low-frequency mutations. However, all 27 individuals who experienced virological failure reported poor adherence to their drug regimen during the preceding follow-up time, and all 19 who subsequently improved their adherence achieved viral suppression at later time points, consistent with a lack of clinical resistance. In conclusion, in a population with low antiviral resistance emergence, NGS methods detected numerous instances of minor alleles that did not result in subsequent bona fide virological failure due to antiviral resistance.

Copy Number Heterogeneity of JC Virus Standards.

Quantitative PCR is a diagnostic mainstay of clinical virology, and accurate quantitation of viral load among labs requires the use of international standards. However, the use of multiple passages of viral isolates to obtain sufficient material for international standards may result in genomic changes that complicate their use as quantitative standards. We performed next-generation sequencing to obtain single-nucleotide resolution and relative copy number of JC virus (JCV) clinical standards. Strikingly, the WHO international standard and the Exact v1/v2 prototype standards for JCV showed 8-fold and 4-fold variation in genomic coverage between different loci in the viral genome, respectively, due to large deletions in the large T antigen region. Intriguingly, several of the JCV standards sequenced in this study with large T antigen deletions were cultured in cell lines immortalized using simian virus 40 (SV40) T antigen, suggesting the possibility of transcomplementation in cell culture. Using a cutoff 5% allele fraction for junctional reads, 7 different rearrangements were present in the JC virus sequences present in the WHO standard across multiple library preparations and sequencing runs. Neither the copy number differences nor the rearrangements were observed in a clinical sample with a high copy number of JCV or a plasmid control. These results were also confirmed by the quantitative real-time PCR (qPCR), droplet digital PCR (ddPCR), and Sanger sequencing of multiple rearrangements. In summary, targeting different regions of the same international standard can result in up to an 8-fold difference in quantitation. We recommend the use of next-generation sequencing to validate standards in clinical virology.