A Next-Generation Sequencing Data Analysis Pipeline for Detecting Unknown Pathogens from Mixed Clinical Samples and Revealing Their Genetic Diversity

Novel intertypic recombination of Echovirus 11 in the Enterovirus species B

Enteroviruses (EVs), single-stranded, positive-sense RNA viruses, can be classified into four species (A-D), which have previously been linked to a diverse range of disease manifestations and infections affecting the central nervous system. In the Enterovirus species B (EV-B), Echovirus type 11 (E11) has been observed to occasionally circulate in Taiwan, which was responsible for an epidemic of enterovirus infections in 2018. Here, 48 clinical specimens isolated in 2003, 2004, 2009, and 2018 were collected for the high-throughput sequencing. Notably, we identified 2018 Taiwanese strains having potential recombinations in the 3D gene, as well as one 2003 strain having a double recombination with E6 and Coxsackievirus B5 in the P2 and P3 regions, respectively. Additionally, one amino acid signature mutated from the Histidine (H) in throat swab specimens to the Tyrosine (Y) in cerebral spinal fluid specimens was detected at position 1496 (or 57) of the genomic coordinate (or 3A gene) to further demonstrate intra-host evolution in different organs. In conclusion, this study identifies potential intertypic recombination events and an intra-host signature mutation in E11 strains, isolated during a 2018 neurological disease outbreak in Taiwan, contributing to our understanding of its evolution and pathogenesis.

Detection of respiratory viruses directly from clinical samples using next-generation sequencing: A literature review of recent advances and potential for routine clinical use

Acute respiratory infection is the third most frequent cause of mortality worldwide, causing over 4.25 million deaths annually. Although most diagnosed acute respiratory infections are thought to be of viral origin, the aetiology often remains unclear. The advent of next‐generation sequencing (NGS) has revolutionised the field of virus discovery and identification, particularly in the detection of unknown respiratory viruses. We systematically reviewed the application of NGS technologies for detecting respiratory viruses from clinical samples and outline potential barriers to the routine clinical introduction of NGS. The five databases searched for studies published in English from 01 January 2010 to 01 February 2021, which led to the inclusion of 52 studies. A total of 14 different models of NGS platforms were summarised from included studies. Among these models, second‐generation sequencing platforms (e.g., Illumina sequencers) were used in the majority of studies (41/52, 79%). Moreover, NGS platforms have proven successful in detecting a variety of respiratory viruses, including influenza A/B viruses (9/52, 17%), SARS‐CoV‐2 (21/52, 40%), parainfluenza virus (3/52, 6%), respiratory syncytial virus (1/52, 2%), human metapneumovirus (2/52, 4%), or a viral panel including other respiratory viruses (16/52, 31%). The review of NGS technologies used in previous studies indicates the advantages of NGS technologies in novel virus detection, virus typing, mutation identification, and infection cluster assessment. Although there remain some technical and ethical challenges associated with NGS use in clinical laboratories, NGS is a promising future tool to improve understanding of respiratory viruses and provide a more accurate diagnosis with simultaneous virus characterisation.

Enterovirus A71 Induces Neurological Diseases and Dynamic Variants in Oral Infection of Human SCARB2-Transgenic Weaned Mice

Enterovirus A71 (EV-A71) and many members of the Picornaviridae family are neurotropic pathogens of global concern. These viruses are primarily transmitted through the fecal-oral route, and thus suitable animal models of oral infection are needed to investigate viral pathogenesis. An animal model of oral infection was developed using transgenic mice expressing human SCARB2 (hSCARB2 Tg), murine-adapted EV-A71/MP4 virus, and EV-A71/MP4 virus with an engineered nanoluciferase gene that allows imaging of viral replication and spread in infected mice. Next-generation sequencing of EV-A71 genomes in the tissues and organs of infected mice was also performed. Oral inoculation of EV-A71/MP4 or nanoluciferase-carrying MP4 virus stably induced neurological symptoms and death in infected 21-day-old weaned mice. In vivo bioluminescence imaging of infected mice and tissue immunostaining of viral antigens indicated that orally inoculated virus can spread to the central nervous system (CNS) and other tissues. Next-generating sequencing further identified diverse mutations in viral genomes that can potentially contribute to viral pathogenesis. This study presents an EV-A71 oral infection murine model that efficiently infects weaned mice and allows tracking of viral spread, features that can facilitate research into viral pathogenesis and neuroinvasion via the natural route of infection. IMPORTANCE Enterovirus A71 (EV-A71), a positive-strand RNA virus of the Picornaviridae, poses a persistent global public health problem. EV-A71 is primarily transmitted through the fecal-oral route, and thus suitable animal models of oral infection are needed to investigate viral pathogenesis. We present an animal model of EV-A71 infection that enables the natural route of oral infection in weaned and nonimmunocompromised 21-day-old hSCARB2 transgenic mice. Our results demonstrate that severe disease and death could be stably induced, and viral invasion of the CNS could be replicated in this model, similar to severe real-world EV-A71 infections. We also developed a nanoluciferase-containing EV-A71 virus that can be used with this animal model to track viral spread after oral infection in real time. Such a model offers several advantages over existing animal models and can facilitate future research into viral spread, tissue tropism, and viral pathogenesis, all pressing issues that remain unaddressed for EV-A71 infections.

SARS-CoV-2 genomic surveillance in Taiwan revealed novel ORF8-deletion mutant and clade possibly associated with infections in Middle East

Taiwan experienced two waves of imported infections with Coronavirus Disease 2019 (COVID-19). This study aimed at investigating the genomic variation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Taiwan and compared their evolutionary trajectories with the global strains. We performed culture and full-genome sequencing of SARS-CoV-2 strains followed by phylogenetic analysis. A 382-nucleotides deletion in open reading frame 8 (ORF8) was found in a Taiwanese strain isolated from a patient on February 4, 2020 who had a travel history to Wuhan. Patients in the first wave also included several sporadic, local transmission cases. Genomes of 5 strains sequenced from clustered infections were classified into a new clade with ORF1ab-V378I mutation, in addition to 3 dominant clades ORF8-L84S, ORF3a-G251V and S-D614G. This highlighted clade also included some strains isolated from patients who had a travel history to Turkey and Iran. The second wave mostly resulted from patients who had a travel history to Europe and Americas. All Taiwanese viruses were classified into various clades. Genomic surveillance of SARS-CoV-2 in Taiwan revealed a new ORF8-deletion mutant and a virus clade that may be associated with infections in the Middle East, which contributed to a better understanding of the global SARS-CoV-2 transmission dynamics.

Validation of the easyscreen flavivirus dengue alphavirus detection kit based on 3base amplification technology and its application to the 2016/17 Vanuatu dengue outbreak

BACKGROUND

The family flaviviridae and alphaviridae contain a diverse group of pathogens that cause significant morbidity and mortality worldwide. Diagnosis of the virus responsible for disease is essential to ensure patients receive appropriate clinical management. Very few real-time RT-PCR based assays are able to detect the presence of all members of these families using a single primer and probe set. We have developed a novel chemistry, 3base, which simplifies the viral nucleic acids allowing the design of RT-PCR assays capable of pan-family identification.

METHODOLOGY/PRINCIPAL FINDING

Synthetic constructs, viral nucleic acids, intact viral particles and characterised reference materials were used to determine the specificity and sensitivity of the assays. Synthetic constructs demonstrated the sensitivities of the pan-flavivirus detection component were in the range of 13 copies per PCR. The pan-alphavirus assay had a sensitivity range of 10-25 copies per reaction depending on the viral strain. Lower limit of detection studies using whole virus particles demonstrated that sensitivity for assays was in the range of 1-2 copies per reaction. No cross reactivity was observed with a number of commonly encountered viral strains. Proficiency panels showed 100% concordance with the expected results and the assays performed as well as, if not better than, other assays used in laboratories worldwide. After initial assay validation the pan-viral assays were then tested during the 2016-2017 Vanuatu dengue-2 outbreak. Positive results were detected in 116 positives from a total of 187 suspected dengue samples.

CONCLUSIONS/SIGNIFICANCE

The pan-viral screening assays described here utilise a novel 3base technology and are shown to provide a sensitive and specific method to screen and thereafter speciate flavi- and/or alpha- viruses in clinical samples. The assays performed well in an outbreak situation and can be used to detect positive clinical samples containing any flavivirus or alphavirus in approximately 3 hours 30 minutes.

Genomic analysis of serologically untypable human enteroviruses in Taiwan

BACKGROUND

Human enteroviruses contain over 100 serotypes. We have routinely conducted enterovirus surveillance in northern Taiwan; but about 10% of isolates could not be serotyped using traditional assays. Next-generation sequencing (NGS) is a powerful tool for genome sequencing.

METHODS

In this study, we established an NGS platform to conduct genome sequencing for the serologically untypable enterovirus isolates.

RESULTS

Among 130 serologically untypable isolates, 121 (93%) of them were classified into 29 serotypes using CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primer)-based RT-PCR to amplify VP1 genes (VP1-CODEHOP). We further selected 52 samples for NGS and identified 59 genome sequences from 51 samples, including 8 samples containing two virus genomes. We also detected 23 genome variants (nucleotide identity < 90% compared with genome sequences in the public domain) which were potential genetic recombination, including 9 inter-serotype recombinants and 14 strains with unknown sources of recombination.

CONCLUSIONS

We successfully integrated VP1-CODEHOP and NGS techniques to conduct genomic analysis of serologically untypable enteroviruses.

Evolutionary analysis of human parechovirus type 3 and clinical outcomes of infection during the 2017-18 Australian epidemic

Human parechovirus type 3 (HPeV3) can cause severe sepsis-like illness in young infants and may be associated with long term neurodevelopmental delay later in childhood. We investigated the molecular epidemiology of HPeV infection in thirty three infants requiring hospitalization before, during and after the peak of the 2017/18 HPeV epidemic wave in Australia. During the peak of the epidemic, all cases were infected with an HPeV3, while before and after the peak, HPeV1 was the predominant type detected. The predominant HPeV3 was the recombinant HPeV3 also detected in the 2013/14 and 2015/16 Australian epidemics. Sepsis-like or meningitis-like symptoms were only reported in cases infected with the recombinant HPeV3. Phylogenetic analysis of the recombinant HPeV3 revealed that the virus continued to evolve, also between the Australian outbreaks, thus indicating continued circulation, despite not being detected and reported in Australia or elsewhere in between epidemic waves. The recombinant HPeV3 continued to show a remarkable stability in its capsid amino acid sequence, further strengthening our previous argument for development of a vaccine or immunotherapeutics to reduce the severity of HPeV3 outbreaks due to this virus.

Population dynamics at neuraminidase position 151 of influenza A (H1N1)pdm09 virus in clinical specimens

Influenza A virus mutates rapidly, allowing it to escape natural and vaccine-induced immunity. Neuraminidase (NA) is a surface protein capable of cleaving the glycosidic linkages of neuraminic acids to release newly formed virions from infected cells. Genetic variants within a viral population can influence the emergence of pandemic viruses as well as drug susceptibility and vaccine effectiveness. In the present study, 55 clinical specimens from patients infected with the 2009 pandemic influenza A/H1N1 virus, abbreviated as A(H1N1)pdm09, during the 2015-2016 outbreak season in Taiwan were collected. Whole genomes were obtained through next-generation sequencing. Based on the published sequences from A(H1N1)pdm09 strains worldwide, a mixed population of two distinct variants at NA position 151 was revealed. We initially reasoned that such a mixed population may have emerged during cell culture. However, additional investigations confirmed that these mixed variants were detectable in the specimens of patients. To further investigate the role of the two NA-151 variants in a dynamic population, a reverse genetics system was employed to generate recombinant A(H1N1)pdm09 viruses. It was observed that the mixture of the two distinct variants was characterized by a higher replication rate compared to the recombinant viruses harbouring a single variant. Moreover, an NA inhibition assay revealed that a high frequency of the minor NA-151 variant in A(H1N1)pdm09 was associated with a reduced susceptibility to NA inhibitors. We conclude that two distinct NA-151 variants can be identified in patient specimens and that such variants may increase viral replication and NA activity.

Discovery of Enterovirus A71-like nonstructural genomes in recent circulating viruses of the Enterovirus A species

Enterovirus A71 (EV-A71) is an important nonpolio enterovirus that causes severe neurological complications. In 1998, Taiwan experienced an EV-A71 outbreak that caused 78 deaths. Since then, periodic epidemics of EV-A71 associated with newly emerging strains have occurred. Several of these strains are known to be recombinant; however, how these strains arose within such a short period of time remains unknown. Here, we sequenced 64 full-length genomes from clinical isolates collected from 2005 to 2016 and incorporated all 91 Taiwanese genomes downloaded from the Virus Pathogen Resource to extensively analyze EV-A71 recombination in Taiwan. We found that the B3 subgenotype was a potential recombinant parent of the EV-A71 C2-like and C4 strains by intratypic recombination. Such B3-similar regions were also found in many cocirculating coxsackieviruses belonging to Enterovirus A species (EV-A) through a series of intertypic recombinations. Therefore, locally enriched outbreaks of cocirculating viruses from different genotypes/serotypes may facilitate recombination. Most recombination breakpoints we found had nonrandom distributions and were located within the region spanning from the boundary of P1 (structural gene) and P2 (nonstructural) to the cis-acting replication element at P2, indicating that specific genome reassembly of structural and nonstructural genes may be subject to natural selection. Through intensive recombination, 11 EV-A71-like signatures (including one in 3A, two in 3C, and eight in 3D) were found to be present in a variety of recently cocirculating EV-A viruses worldwide, suggesting that these viruses may be targets for wide-spectrum antiviral development.

Overview of Virus Metagenomic Classification Methods and Their Biological Applications

Metagenomics poses opportunities for clinical and public health virology applications by offering a way to assess complete taxonomic composition of a clinical sample in an unbiased way. However, the techniques required are complicated and analysis standards have yet to develop. This, together with the wealth of different tools and workflows that have been proposed, poses a barrier for new users. We evaluated 49 published computational classification workflows for virus metagenomics in a literature review. To this end, we described the methods of existing workflows by breaking them up into five general steps and assessed their ease-of-use and validation experiments. Performance scores of previous benchmarks were summarized and correlations between methods and performance were investigated. We indicate the potential suitability of the different workflows for (1) time-constrained diagnostics, (2) surveillance and outbreak source tracing, (3) detection of remote homologies (discovery), and (4) biodiversity studies. We provide two decision trees for virologists to help select a workflow for medical or biodiversity studies, as well as directions for future developments in clinical viral metagenomics.

A metagenomics study for the identification of respiratory viruses in mixed clinical specimens: an application of the iterative mapping approach

Metagenomic approaches to detect viral genomes and variants in clinical samples have various challenges, including low viral titers and bacterial and human genome contamination. To address these limitations, we examined a next-generation sequencing (NGS) and iterative mapping approach for virus detection in clinical samples. We analyzed 40 clinical specimens from hospitalized children diagnosed with acute bronchiolitis, croup, or respiratory tract infections in which virus identification by viral culture or polymerase chain reaction (PCR) was unsuccessful. For our NGS data analysis pipeline, clinical samples were pooled into two NGS groups to reduce sequencing costs, and the depth and coverage of assembled contigs were effectively increased using an iterative mapping approach. PCR was individually performed for each specimen according to the NGS-predicted viral type. We successfully detected previously unidentified respiratory viruses in 26 of 40 specimens using our proposed NGS pipeline. Two dominant populations within the detected viruses were human rhinoviruses (HRVs; n = 14) and human coronavirus NL63 (n = 8), followed by human parainfluenza virus (HPIV), human parechovirus, influenza A virus, respiratory syncytial virus (RSV), and human metapneumovirus. This is the first study reporting the complete genome sequences of HRV-A101, HRV-C3, HPIV-4a, and RSV, as well as an analysis of their genetic variants, in Taiwan. These results demonstrate that this NGS pipeline allows to detect viruses which were not identified by routine diagnostic assays, directly from clinical samples.

Intertypic recombination of human parechovirus 4 isolated from infants with sepsis-like disease

BACKGROUND

Human parechoviruses (HPeVs) (family Picornaviridae), are common pathogens in young children. Despite their high prevalence, research on their genetic identity, diversity and evolution have remained scarce.

OBJECTIVES

Complete coding regions of three previously reported HPeV-4 isolates from Finnish children with sepsis-like disease were sequenced in order to elucidate the phylogenetic relationships and potential recombination events during the evolution of these isolates.

STUDY DESIGN

The isolated viruses were sequenced and aligned with all HPeV complete genome sequences available in GenBank. Phylogenetic trees were constructed and similarity plot and bootscanning methods were used for recombination analysis.

RESULTS

The three HPeV-4 isolates had 99.8% nucleotide sequence similarity. The phylogenetic analysis indicated that capsid-encoding sequences of these HPeV-4 isolates were closely related to other HPeV-4 strains (80.7-94.7% nucleotide similarity), whereas their non-structural region genes 2A to 3C clustered together with several HPeV-1 and HPeV-3 strains, in addition to the HPeV-4 strain K251176-02 (isolated 2002 in the Netherlands), but not with other HPeV-4 strains. However, in 3D-encoding sequence the Finnish HPeV-4 isolates did not cluster with the strain HPeV-4/K251176-02, but instead, formed a distinct group together with several HPeV-1 and HPeV-3 strains. Similarity plot and Bootscan analyses further confirmed intertypic recombination events in the evolution of the Finnish HPeV-4 isolates.

CONCLUSION

Intertypic recombination event(s) have occurred during the evolution of HPeV-4 isolates from children with sepsis-like disease. However, due to the low number of parechovirus complete genomes available, the precise recombination partners could not be detected. The results suggest frequent intratypic recombination among parechoviruses.