Identification, characterization, and natural selection of mutations driving airborne transmission of A/H5N1 virus.

Recently, A/H5N1 influenza viruses were shown to acquire airborne transmissibility between ferrets upon targeted mutagenesis and virus passage. The critical genetic changes in airborne A/Indonesia/5/05 were not yet identified. Here, five substitutions proved to be sufficient to determine this airborne transmission phenotype. Substitutions in PB1 and PB2 collectively caused enhanced transcription and virus replication. One substitution increased HA thermostability and lowered the pH of membrane fusion. Two substitutions independently changed HA binding preference from α2,3-linked to α2,6-linked sialic acid receptors. The loss of a glycosylation site in HA enhanced overall binding to receptors. The acquired substitutions emerged early during ferret passage as minor variants and became dominant rapidly. Identification of substitutions that are essential for airborne transmission of avian influenza viruses between ferrets and their associated phenotypes advances our fundamental understanding of virus transmission and will increase the value of future surveillance programs and public health risk assessments.

Functional properties of measles virus proteins derived from a subacute sclerosing panencephalitis patient who received repeated remdesivir treatments.

Subacute sclerosing panencephalitis (SSPE) is a rare but fatal late neurological complication of measles, caused by persistent measles virus (MeV) infection of the central nervous system. There are no drugs approved for the treatment of SSPE. Here, we followed the clinical progression of a 5-year-old SSPE patient after treatment with the nucleoside analog remdesivir, conducted a post-mortem evaluation of the patient's brain, and characterized the MeV detected in the brain. The quality of life of the patient transiently improved after the first two courses of remdesivir, but a third course had no further clinical effect, and the patient eventually succumbed to his condition. Post-mortem evaluation of the brain displayed histopathological changes including loss of neurons and demyelination paired with abundant presence of MeV RNA-positive cells throughout the brain. Next-generation sequencing of RNA isolated from the brain revealed a complete MeV genome with mutations that are typically detected in SSPE, characterized by a hypermutated M gene. Additional mutations were detected in the polymerase (L) gene, which were not associated with resistance to remdesivir. Functional characterization showed that mutations in the F gene led to a hyperfusogenic phenotype predominantly mediated by N465I. Additionally, recombinant wild-type-based MeV with the SSPE-F gene or the F gene with the N465I mutation was no longer lymphotropic but instead efficiently disseminated in neural cultures. Altogether, this case encourages further investigation of remdesivir as a potential treatment of SSPE and highlights the necessity to functionally understand SSPE-causing MeV.IMPORTANCEMeasles virus (MeV) causes acute, systemic disease and remains an important cause of morbidity and mortality in humans. Despite the lack of known entry receptors in the brain, MeV can persistently infect the brain causing the rare but fatal neurological disorder subacute sclerosing panencephalitis (SSPE). SSPE-causing MeVs are characterized by a hypermutated genome and a hyperfusogenic F protein that facilitates the rapid spread of MeV throughout the brain. No treatment against SSPE is available, but the nucleoside analog remdesivir was recently demonstrated to be effective against MeV in vitro. We show that treatment of an SSPE patient with remdesivir led to transient clinical improvement and did not induce viral escape mutants, encouraging the future use of remdesivir in SSPE patients. Functional characterization of the viral proteins sheds light on the shared properties of SSPE-causing MeVs and further contributes to understanding how those viruses cause disease.

Limited airborne transmission of H7N9 influenza A virus between ferrets.

Wild waterfowl form the main reservoir of influenza A viruses, from which transmission occurs directly or indirectly to various secondary hosts, including humans. Direct avian-to-human transmission has been observed for viruses of subtypes A(H5N1), A(H7N2), A(H7N3), A(H7N7), A(H9N2) and A(H10N7) upon human exposure to poultry, but a lack of sustained human-to-human transmission has prevented these viruses from causing new pandemics. Recently, avian A(H7N9) viruses were transmitted to humans, causing severe respiratory disease and deaths in China. Because transmission via respiratory droplets and aerosols (hereafter referred to as airborne transmission) is the main route for efficient transmission between humans, it is important to gain an insight into airborne transmission of the A(H7N9) virus. Here we show that although the A/Anhui/1/2013 A(H7N9) virus harbours determinants associated with human adaptation and transmissibility between mammals, its airborne transmissibility in ferrets is limited, and it is intermediate between that of typical human and avian influenza viruses. Multiple A(H7N9) virus genetic variants were transmitted. Upon ferret passage, variants with higher avian receptor binding, higher pH of fusion, and lower thermostability were selected, potentially resulting in reduced transmissibility. This A(H7N9) virus outbreak highlights the need for increased understanding of the determinants of efficient airborne transmission of avian influenza viruses between mammals.

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia.

A previously unknown coronavirus was isolated from the sputum of a 60-year-old man who presented with acute pneumonia and subsequent renal failure with a fatal outcome in Saudi Arabia. The virus (called HCoV-EMC) replicated readily in cell culture, producing cytopathic effects of rounding, detachment, and syncytium formation. The virus represents a novel betacoronavirus species. The closest known relatives are bat coronaviruses HKU4 and HKU5. Here, the clinical data, virus isolation, and molecular identification are presented. The clinical picture was remarkably similar to that of the severe acute respiratory syndrome (SARS) outbreak in 2003 and reminds us that animal coronaviruses can cause severe disease in humans.

Excessive production and extreme editing of human metapneumovirus defective interfering RNA is associated with type I IFN induction.

Type I IFN production is one of the hallmarks of host innate immune responses upon virus infection. Whilst most respiratory viruses carry IFN antagonists, reports on human metapneumovirus (HMPV) have been conflicting. Using deep sequencing, we have demonstrated that HMPV particles accumulate excessive amounts of defective interfering RNA (DIs) rapidly upon in vitro passage, and that these are associated with IFN induction. Importantly, the DIs were edited extensively; up to 70% of the original A and T residues had mutated to G or C, respectively. Such high editing rates of viral RNA have not, to our knowledge, been reported before. Bioinformatics and PCR assays indicated that adenosine deaminase acting on RNA (ADAR) was the most likely editing enzyme. HMPV thus has an unusually high propensity to generate DIs, which are edited at an unprecedented high frequency. The conflicting published data on HMPV IFN induction and antagonism are probably explained by DIs in virus stocks. The interaction of HMPV DIs with the RNA-editing machinery and IFN responses warrants further investigation.

Multicenter benchmarking of short and long read wet lab protocols for clinical viral metagenomics.

Metagenomics is gradually being implemented for diagnosing infectious diseases. However, in-depth protocol comparisons for viral detection have been limited to individual sets of experimental workflows and laboratories. In this study, we present a benchmark of metagenomics protocols used in clinical diagnostic laboratories initiated by the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS). A mock viral reference panel was designed to mimic low biomass clinical specimens. The panel was used to assess the performance of twelve metagenomic wet lab protocols currently in use in the diagnostic laboratories of participating ENNGS member institutions. Both Illumina and Nanopore, shotgun and targeted capture probe protocols were included. Performance metrics sensitivity, specificity, and quantitative potential were assessed using a central bioinformatics pipeline. Overall, viral pathogens with loads down to 104 copies/ml (corresponding to CT values of 31 in our PCR assays) were detected by all the evaluated metagenomic wet lab protocols. In contrast, lower abundant mixed viruses of CT values of 35 and higher were detected only by a minority of the protocols. Considering the reference panel as the gold standard, optimal thresholds to define a positive result were determined per protocol, based on the horizontal genome coverage. Implementing these thresholds, sensitivity and specificity of the protocols ranged from 67 to 100 % and 87 to 100 %, respectively. A variety of metagenomic protocols are currently in use in clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implying the need for standardization of metagenomic analysis for use in clinical settings.

Arterivirus and nairovirus ovarian tumor domain-containing Deubiquitinases target activated RIG-I to control innate immune signaling.

The innate immune response constitutes the first line of defense against viral infection and is extensively regulated through ubiquitination. The removal of ubiquitin from innate immunity signaling factors by deubiquitinating enzymes (DUBs) therefore provides a potential opportunity for viruses to evade this host defense system. It was previously found that specific proteases encoded by the unrelated arteri- and nairoviruses resemble the ovarian tumor domain-containing (OTU) family of DUBs. In arteriviruses, this domain has been characterized before as a papain-like protease (PLP2) that is also involved in replicase polyprotein processing. In nairoviruses, the DUB resides in the polymerase protein but is not essential for RNA replication. Using both in vitro and cell-based assays, we now show that PLP2 DUB activity is conserved in all members of the arterivirus family and that both arteri- and nairovirus DUBs inhibit RIG-I-mediated innate immune signaling when overexpressed. The potential relevance of RIG-I-like receptor (RLR) signaling for the innate immune response against arterivirus infection is supported by our finding that in mouse embryonic fibroblasts, the production of beta interferon primarily depends on the recognition of arterivirus RNA by the pattern-recognition receptor MDA5. Interestingly, we also found that both arteri- and nairovirus DUBs inhibit RIG-I ubiquitination upon overexpression, suggesting that both MDA5 and RIG-I have a role in countering infection by arteriviruses. Taken together, our results support the hypothesis that arteri- and nairoviruses employ their deubiquitinating potential to inactivate cellular proteins involved in RLR-mediated innate immune signaling, as exemplified by the deubiquitination of RIG-I.

Looking back on the COVID-19 pandemic in an elite sports team using whole genome sequencing.

OBJECTIVES: The aim of this study was to investigate the effectiveness of infection control measures to prevent transmission of SARS-CoV-2 within a professional sports team using whole genome sequencing. DESIGN: Prospective cohort study. METHODS: 74 players and staff members of a Dutch professional male football team were followed from August 2020 until May 2021. A set of health and safety measures were introduced and all participants underwent regular SARS-CoV-2 RNA testing. All positive samples were subsequently sequenced (Nanopore sequencing) to assess whether infections were acquired within the training center or in the community. RESULTS: Throughout the study period, 13 participants tested positive for SARS-CoV-2. The phylogenetic analysis revealed 2 clusters (of 2 and 3 cases respectively), indicating that 3/13 cases (23%) acquired infection from another player or staff member. The first cluster was diagnosed upon enrolment, thus transmission had occurred prior to the implementation of health and safety protocols. Finally, 4 cases were diagnosed prior to symptom onset, emphasizing that frequent testing leads to early detection and isolation. CONCLUSIONS: These data show that a combination of regular testing and basic control measures can prevent outbreaks of COVID-19 in a professional sports team. Whole genome sequencing is an important tool to distinguish between infections introduced from the community and infections transmitted between athletes.

Case Report: Late Reactivation of Herpes B Virus After a Monkey Bite: A Case of Severe Meningoencephalitis.

Macacine alphaherpesvirus 1, also known as herpes B virus (BV), is an alphaherpesvirus endemic to several macaque species, capable of causing zoonotic infections in humans, with high mortality rates. Evidence of reactivation in humans has rarely been reported. Here we depict a case of BV reactivation after 54 years, leading to severe meningoencephalitis. This case supports the use of antiviral prophylaxis in patients surviving a confirmed BV central nervous system infection. We sequenced DNA from BV obtained from the patient's cerebrospinal fluid. Phylogenetic analysis showed significant divergence in the clustering of this particular BV strain compared with other known BVs. Therefore, additional efforts are needed to obtain a broader sequence landscape from BVs circulating in monkeys.

Performance of Five Metagenomic Classifiers for Virus Pathogen Detection Using Respiratory Samples from a Clinical Cohort.

Viral metagenomics is increasingly applied in clinical diagnostic settings for detection of pathogenic viruses. While several benchmarking studies have been published on the use of metagenomic classifiers for abundance and diversity profiling of bacterial populations, studies on the comparative performance of the classifiers for virus pathogen detection are scarce. In this study, metagenomic data sets (n = 88) from a clinical cohort of patients with respiratory complaints were used for comparison of the performance of five taxonomic classifiers: Centrifuge, Clark, Kaiju, Kraken2, and Genome Detective. A total of 1144 positive and negative PCR results for a total of 13 respiratory viruses were used as gold standard. Sensitivity and specificity of these classifiers ranged from 83 to 100% and 90 to 99%, respectively, and was dependent on the classification level and data pre-processing. Exclusion of human reads generally resulted in increased specificity. Normalization of read counts for genome length resulted in a minor effect on overall performance, however it negatively affected the detection of targets with read counts around detection level. Correlation of sequence read counts with PCR Ct-values varied per classifier, data pre-processing (R2 range 15.1-63.4%), and per virus, with outliers up to 3 log10 reads magnitude beyond the predicted read count for viruses with high sequence diversity. In this benchmarking study, sensitivity and specificity were within the ranges of use for diagnostic practice when the cut-off for defining a positive result was considered per classifier.

Diminished amplification of SARS-CoV-2 ORF1ab in a commercial dual-target qRT-PCR diagnostic assay.

Here we describe a SARS-CoV-2 variant with diminished amplification of the ORF ORF1ab target in the Cobas® dual-target SARS-CoV-2 assay resulting in a discrepancy of Ct-values (Ct-value 20.7 for the E-gene and Ct-value 30.2 for ORF1ab). Five unique nucleotide mutations were identified in ORF1ab: C11450A (nsp10) C14178T (RdRp), G15006T (RdRp), G18394T (Hel), and G20995T (Hel). This case highlights the importance of surveillance of genomic regions used in molecular diagnostics and the importance of the public release of target regions used to update commercial and in-house developed SARS-CoV-2 PCR tests. This work underpins the importance of using dual-targets in molecular diagnostic assays to limit the change of false-negative results due to primer and/or probe mismatches.

HIV-1 resistance against dolutegravir fluctuates rapidly alongside erratic treatment adherence: a case report.

OBJECTIVES: We report a case of incomplete HIV-1 suppression on a dolutegravir, lamivudine, and abacavir single-tablet regimen with the emergence of the H51Y and G118R integrase resistance mutations. METHODS: Integrase sequencing was performed retrospectively by Sanger and next-generation sequencing. Rates of emergence and decline of resistance mutations were calculated using next-generation sequencing data. Dolutegravir plasma concentrations were measured by ultra-performance liquid chromatography-tandem mass spectrometry. The effects of H51Y and G118R on infectivity, fitness, and susceptibility to dolutegravir were quantified using cell-based assays. RESULTS: During periods of non-adherence to treatment, mutations were retrospectively documented only by next-generation sequencing. Misdiagnosis by Sanger sequencing was caused by the rapid decline of mutant strains within the retroviral population. This observation was also true for a M184V lamivudine-resistant reverse transcriptase mutation found in association with integrase mutations on single HIV genomes. Resistance rebound upon treatment re-initiation was swift (>8000 copies per day). Next-generation sequencing indicated cumulative adherence to treatment. Compared to WT HIV-1, relative infectivity was 73%, 38%, and 43%; relative fitness was 100%, 35%, and 10% for H51Y, G118R, and H51Y+G118R viruses, respectively. H51Y did not change the susceptibility to dolutegravir, but G188R and H51Y+G118R conferred 7- and 28-fold resistance, respectively. CONCLUSION: This case illustrates how poorly-fit drug-resistant viruses wax and wane alongside erratic treatment adherence and are easily misdiagnosed by Sanger sequencing. We recommend next-generation sequencing to improve the clinical management of incomplete virological suppression with dolutegravir.

Species-specific antagonism of host ISGylation by the influenza B virus NS1 protein.

Interferon-stimulated expression and conjugation of the ubiquitin-like modifier ISG15 restricts replication of several viruses. Here, we established complete E1-activating, E2-conjugating, and E3 ligase-dependent expression systems for assaying both human and mouse ISGylation. We confirm that human HerC5, but not human HerC6, has ISG15 E3 ligase activity and identify mouse HerC6 as a bona fide ISG15 E3 ligase. Furthermore, we demonstrate that influenza B virus NS1 protein potently antagonizes human but not mouse ISGylation, a property dependent on B/NS1 binding the N-terminal domain of human but not mouse ISG15. Using chimeric human/mouse ISG15 constructs, we show that the B/NS1:ISG15 interaction is both necessary and sufficient to inhibit ISGylation regardless of the ligation machinery used. Inability to block ISGylation in certain species may contribute to limiting influenza B virus host range.

Benchmark of thirteen bioinformatic pipelines for metagenomic virus diagnostics using datasets from clinical samples.

INTRODUCTION: Metagenomic sequencing is increasingly being used in clinical settings for difficult to diagnose cases. The performance of viral metagenomic protocols relies to a large extent on the bioinformatic analysis. In this study, the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS) initiated a benchmark of metagenomic pipelines currently used in clinical virological laboratories. METHODS: Metagenomic datasets from 13 clinical samples from patients with encephalitis or viral respiratory infections characterized by PCR were selected. The datasets were analyzed with 13 different pipelines currently used in virological diagnostic laboratories of participating ENNGS members. The pipelines and classification tools were: Centrifuge, DAMIAN, DIAMOND, DNASTAR, FEVIR, Genome Detective, Jovian, MetaMIC, MetaMix, One Codex, RIEMS, VirMet, and Taxonomer. Performance, characteristics, clinical use, and user-friendliness of these pipelines were analyzed. RESULTS: Overall, viral pathogens with high loads were detected by all the evaluated metagenomic pipelines. In contrast, lower abundance pathogens and mixed infections were only detected by 3/13 pipelines, namely DNASTAR, FEVIR, and MetaMix. Overall sensitivity ranged from 80% (10/13) to 100% (13/13 datasets). Overall positive predictive value ranged from 71-100%. The majority of the pipelines classified sequences based on nucleotide similarity (8/13), only a minority used amino acid similarity, and 6 of the 13 pipelines assembled sequences de novo. No clear differences in performance were detected that correlated with these classification approaches. Read counts of target viruses varied between the pipelines over a range of 2-3 log, indicating differences in limit of detection. CONCLUSION: A wide variety of viral metagenomic pipelines is currently used in the participating clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implicating the need for standardization and validation of metagenomic analysis for clinical diagnostic use. Future studies should address the selective effects due to the choice of different reference viral databases.

Duration and key determinants of infectious virus shedding in hospitalized patients with coronavirus disease-2019 (COVID-19).

Key questions in COVID-19 are the duration and determinants of infectious virus shedding. Here, we report that infectious virus shedding is detected by virus cultures in 23 of the 129 patients (17.8%) hospitalized with COVID-19. The median duration of shedding infectious virus is 8 days post onset of symptoms (IQR 5-11) and drops below 5% after 15.2 days post onset of symptoms (95% confidence interval (CI) 13.4-17.2). Multivariate analyses identify viral loads above 7 log10 RNA copies/mL (odds ratio [OR] of 14.7 (CI 3.57-58.1; p < 0.001) as independently associated with isolation of infectious SARS-CoV-2 from the respiratory tract. A serum neutralizing antibody titre of at least 1:20 (OR of 0.01 (CI 0.003-0.08; p < 0.001) is independently associated with non-infectious SARS-CoV-2. We conclude that quantitative viral RNA load assays and serological assays could be used in test-based strategies to discontinue or de-escalate infection prevention and control precautions.

Recommendations for the introduction of metagenomic high-throughput sequencing in clinical virology, part I: Wet lab procedure.

Metagenomic high-throughput sequencing (mHTS) is a hypothesis-free, universal pathogen detection technique for determination of the DNA/RNA sequences in a variety of sample types and infectious syndromes. mHTS is still in its early stages of translating into clinical application. To support the development, implementation and standardization of mHTS procedures for virus diagnostics, the European Society for Clinical Virology (ESCV) Network on Next-Generation Sequencing (ENNGS) has been established. The aim of ENNGS is to bring together professionals involved in mHTS for viral diagnostics to share methodologies and experiences, and to develop application recommendations. This manuscript aims to provide practical recommendations for the wet lab procedures necessary for implementation of mHTS for virus diagnostics and to give recommendations for development and validation of laboratory methods, including mHTS quality assurance, control and quality assessment protocols.

Retrospective Validation of a Metagenomic Sequencing Protocol for Combined Detection of RNA and DNA Viruses Using Respiratory Samples from Pediatric Patients.

Viruses are the main cause of respiratory tract infections. Metagenomic next-generation sequencing (mNGS) enables unbiased detection of all potential pathogens. To apply mNGS in viral diagnostics, sensitive and simultaneous detection of RNA and DNA viruses is needed. Herein, were studied the performance of an in-house mNGS protocol for routine diagnostics of viral respiratory infections with potential for automated pan-pathogen detection. The sequencing protocol and bioinformatics analysis were designed and optimized, including exogenous internal controls. Subsequently, the protocol was retrospectively validated using 25 clinical respiratory samples. The developed protocol using Illumina NextSeq 500 sequencing showed high repeatability. Use of the National Center for Biotechnology Information's RefSeq database as opposed to the National Center for Biotechnology Information's nucleotide database led to enhanced specificity of classification of viral pathogens. A correlation was established between read counts and PCR cycle threshold value. Sensitivity of mNGS, compared with PCR, varied up to 83%, with specificity of 94%, dependent on the cutoff for defining positive mNGS results. Viral pathogens only detected by mNGS, not present in the routine diagnostic workflow, were influenza C, KI polyomavirus, cytomegalovirus, and enterovirus. Sensitivity and analytical specificity of this mNGS protocol were comparable to PCR and higher when considering off-PCR target viral pathogens. One single test detected all potential viral pathogens and simultaneously obtained detailed information on detected viruses.

The respiratory virome and exacerbations in patients with chronic obstructive pulmonary disease.

INTRODUCTION: Exacerbations are major contributors to morbidity and mortality in patients with chronic obstructive pulmonary disease (COPD), and respiratory bacterial and viral infections are an important trigger. However, using conventional diagnostic techniques, a causative agent is not always found. Metagenomic next-generation sequencing (mNGS) allows analysis of the complete virome, but has not yet been applied in COPD exacerbations. OBJECTIVES: To study the respiratory virome in nasopharyngeal samples during COPD exacerbations using mNGS. STUDY DESIGN: 88 nasopharyngeal swabs from 63 patients from the Bergen COPD Exacerbation Study (2006-2010) were analysed by mNGS and in-house qPCR for respiratory viruses. Both DNA and RNA were sequenced simultaneously using an Illumina library preparation protocol with in-house adaptations. RESULTS: By mNGS, 24/88 samples tested positive. Sensitivity and specificity, as compared with PCR, were 96% and 98% for diagnostic targets (23/24 and 1093/1120, respectively). Additional viral pathogens detected by mNGS were herpes simplex virus type 1 and coronavirus OC43. A positive correlation was found between Cq value and mNGS viral normalized species reads (log value) (p = 0.002). Patients with viral pathogens had lower percentages of bacteriophages (p<0.001). No correlation was found between viral reads and clinical markers. CONCLUSIONS: The mNGS protocol used was highly sensitive and specific for semi-quantitative detection of respiratory viruses. Excellent negative predictive value implicates the power of mNGS to exclude any pathogenic respiratory viral infectious cause in one test, with consequences for clinical decision making. Reduced abundance of bacteriophages in COPD patients with viral pathogens implicates skewing of the virome during infection, with potential consequences for the bacterial populations, during infection.

Quasispecies composition and evolution of a typical Zika virus clinical isolate from Suriname.

The arthropod-borne Zika virus (ZIKV) is currently causing a major international public health threat in the Americas. This study describes the isolation of ZIKV from the plasma of a 29-year-old female traveler that developed typical symptoms, like rash, fever and headache upon return from Suriname. The complete genome sequence including the 5' and 3' untranslated regions was determined and phylogenetic analysis showed the isolate clustering within the Asian lineage, close to other viruses that have recently been isolated in the Americas. In addition, the viral quasispecies composition was analyzed by single molecule real time sequencing, which suggested a mutation frequency of 1.4 × 10-4 for this ZIKV isolate. Continued passaging of the virus in cell culture led to the selection of variants with mutations in NS1 and the E protein. The latter might influence virus binding to cell surface heparan sulfate.

Cidofovir gel as treatment of follicular spicules in multiple myeloma.

IMPORTANCE: The cause of follicular spicules in multiple myeloma (MM) is not known. OBSERVATIONS: We present a case of follicular spicules in a patient with MM, which is very reminiscent of trichodysplasia spinulosa caused by a polyomavirus. No trichodysplasia spinulosa-associated polyomavirus could be isolated from the skin lesions; however, the spicules were positive for Merkel cell carcinoma virus, which is also a polyomavirus. CONCLUSIONS AND RELEVANCE: Follicular spicules in MM are probably not caused by the trichodysplasia spinulosa-associated virus. Merkel cell polyomavirus could contribute to the origin of this dermatosis.