Targeted Virome Sequencing Enhances Unbiased Detection and Genome Assembly of Known and Emerging Viruses—The Example of SARS-CoV-2

Targeted virome enrichment and sequencing (VirCapSeq-VERT) utilizes a pool of oligos (baits) to enrich all known—up to 2015—vertebrate-infecting viruses, increasing their detection sensitivity. The hybridisation of the baits to the target sequences can be partial, thus enabling the detection and genomic reconstruction of novel pathogens with <40% genetic diversity compared to the strains used for the baits’ design. In this study, we deploy this method in multiplexed mixes of viral extracts, and we assess its performance in the unbiased detection of DNA and RNA viruses after cDNA synthesis. We further assess its efficiency in depleting various background genomic material. Finally, as a proof-of-concept, we explore the potential usage of the method for the characterization of unknown, emerging human viruses, such as SARS-CoV-2, which may not be included in the baits’ panel. We mixed positive samples of equimolar DNA/RNA viral extracts from SARS-CoV-2, coronavirus OC43, cytomegalovirus, influenza A virus H3N2, parvovirus B19, respiratory syncytial virus, adenovirus C and coxsackievirus A16. Targeted virome enrichment was performed on a dsDNA mix, followed by sequencing on the NextSeq500 (Illumina) and the portable MinION sequencer, to evaluate its usability as a point-of-care (PoC) application. Genome mapping assembly was performed using viral reference sequences. The untargeted libraries contained less than 1% of total reads mapped on most viral genomes, while RNA viruses remained undetected. In the targeted libraries, the percentage of viral-mapped reads were substantially increased, allowing full genome assembly in most cases. Targeted virome sequencing can enrich a broad range of viruses, potentially enabling the discovery of emerging viruses.     

A National Catalogue of Viruses Associated with Indigenous Species Reveals High-Throughput Sequencing as a Driver of Indigenous Virus Discovery

Viruses are important constituents of ecosystems, with the capacity to alter host phenotype and performance. However, virus discovery cued by disease symptoms overlooks latent or beneficial viruses, which are best detected using targeted virus detection or discovered by non-targeted methods, e.g., high-throughput sequencing (HTS). To date, in 64 publications, 701 viruses have been described associated with indigenous species of Aotearoa New Zealand. Viruses were identified in indigenous birds (189 viruses), bats (13 viruses), starfish (4 viruses), insects (280 viruses), and plants (126 viruses). HTS gave rise to a 21.9-fold increase in virus discovery rate over the targeted methods, and 72.7-fold over symptom-based methods. The average number of viruses reported per publication has also increased proportionally over time. The use of HTS has driven the described national virome recently by 549 new-to-science viruses; all are indigenous. This report represents the first catalogue of viruses associated with indigenous species of a country. We provide evidence that the application of HTS to samples of Aotearoa New Zealand’s unique fauna and flora has driven indigenous virus discovery, a key step in the process to understand the role of viruses in the biological diversity and ecology of the land, sea, and air environments of a country.     

HTS-Based Diagnostics of Sugarcane Viruses: Seasonal Variation and Its Implications for Accurate Detection

Rapid global germplasm trade has increased concern about the spread of plant pathogens and pests across borders that could become established, affecting agriculture and environment systems. Viral pathogens are of particular concern due to their difficulty to control once established. A comprehensive diagnostic platform that accurately detects both known and unknown virus species, as well as unreported variants, is playing a pivotal role across plant germplasm quarantine programs. Here we propose the addition of high-throughput sequencing (HTS) from total RNA to the routine quarantine diagnostic workflow of sugarcane viruses. We evaluated the impact of sequencing depth needed for the HTS-based identification of seven regulated sugarcane RNA/DNA viruses across two different growing seasons (spring and fall). Our HTS analysis revealed that viral normalized read counts (RPKM) was up to 23-times higher in spring than in the fall season for six out of the seven viruses. Random read subsampling analyses suggested that the minimum number of reads required for reliable detection of RNA viruses was 0.5 million, with a viral genome coverage of at least 92%. Using an HTS-based total RNA metagenomics approach, we identified all targeted viruses independent of the time of the year, highlighting that higher sequencing depth is needed for the identification of DNA viruses.     

Implementation of GA-VirReport,a Web-Based Bioinformatics Toolkit for Post-Entry Quarantine Screening of Virus and Viroids in Plants

High-throughput sequencing (HTS) of host plant small RNA (sRNA) is a popular approach for plant virus and viroid detection. The major bottlenecks for implementing this approach in routine virus screening of plants in quarantine include lack of computational resources and/or expertise in command-line environments and limited availability of curated plant virus and viroid databases. We developed: (1) virus and viroid report web-based bioinformatics workflows on Galaxy Australia called GA-VirReport and GA-VirReport-Stats for detecting viruses and viroids from host plant sRNA extracts and (2) a curated higher plant virus and viroid database (PVirDB). We implemented sRNA sequencing with unique dual indexing on a set of plants with known viruses. Sequencing data were analyzed using GA-VirReport and PVirDB to validate these resources. We detected all known viruses in this pilot study with no cross-sample contamination. We then conducted a large-scale diagnosis of 105 imported plants processed at the post-entry quarantine facility (PEQ), Australia. We detected various pathogens in 14 imported plants and discovered that de novo assembly using 21–22 nt sRNA fraction and the megablast algorithm yielded better sensitivity and specificity. This study reports the successful, large-scale implementation of HTS and a user-friendly bioinformatics workflow for virus and viroid screening of imported plants at the PEQ.     

Targeted Genome Sequencing (TG-Seq) Approaches to Detect Plant Viruses

Globally, high-throughput sequencing (HTS) has been used for virus detection in germplasm certification programs. However, sequencing costs have impeded its implementation as a routine diagnostic certification tool. In this study, the targeted genome sequencing (TG-Seq) approach was developed to simultaneously detect multiple (four) viral species of; Pea early browning virus (PEBV), Cucumber mosaic virus (CMV), Bean yellow mosaic virus (BYMV) and Pea seedborne mosaic virus (PSbMV). TG-Seq detected all the expected viral amplicons within multiplex PCR (mPCR) reactions. In contrast, the expected PCR amplicons were not detected by gel electrophoresis (GE). For example, for CMV, GE only detected RNA1 and RNA2 while TG-Seq detected all the three RNA components of CMV. In an mPCR to amplify all four viruses, TG-Seq readily detected each virus with more than 732,277 sequence reads mapping to each amplicon. In addition, TG-Seq also detected all four amplicons within a 10⁻⁸ serial dilution that were not detectable by GE. Our current findings reveal that the TG-Seq approach offers significant potential and is a highly sensitive targeted approach for detecting multiple plant viruses within a given biological sample. This is the first study describing direct HTS of plant virus mPCR products. These findings have major implications for grain germplasm healthy certification programs and biosecurity management in relation to pathogen entry into Australia and elsewhere.     

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.     

Filovirus RefSeq Entries: Evaluation and Selection of Filovirus Type Variants,Type Sequences,and Names

Sequence determination of complete or coding-complete genomes of viruses is becoming common practice for supporting the work of epidemiologists, ecologists, virologists, and taxonomists. Sequencing duration and costs are rapidly decreasing, sequencing hardware is under modification for use by non-experts, and software is constantly being improved to simplify sequence data management and analysis. Thus, analysis of virus disease outbreaks on the molecular level is now feasible, including characterization of the evolution of individual virus populations in single patients over time. The increasing accumulation of sequencing data creates a management problem for the curators of commonly used sequence databases and an entry retrieval problem for end users. Therefore, utilizing the data to their fullest potential will require setting nomenclature and annotation standards for virus isolates and associated genomic sequences. The National Center for Biotechnology Information’s (NCBI’s) RefSeq is a non-redundant, curated database for reference (or type) nucleotide sequence records that supplies source data to numerous other databases. Building on recently proposed templates for filovirus variant naming [<virus name> (<strain>)/<isolation host-suffix>/<country of sampling>/<year of sampling>/<genetic variant designation>-<isolate designation>], we report consensus decisions from a majority of past and currently active filovirus experts on the eight filovirus type variants and isolates to be represented in RefSeq, their final designations, and their associated sequences.     

Metagenomics Analysis of the Wheat Virome Identifies Novel Plant and Fungal-Associated Viral Sequences

Wheat viruses including wheat streak mosaic virus, Triticum mosaic virus, and barley yellow dwarf virus cost substantial losses in crop yields every year. Although there have been extensive studies conducted on these known wheat viruses, currently, there is limited knowledge about all components of the wheat (Triticum aestivum L.) virome. Here, we determined the composition of the wheat virome through total RNA deep sequencing of field-collected leaf samples. Sequences were de novo assembled after removing the host reads, and BLASTx searches were conducted. In addition to the documented wheat viruses, novel plant and fungal-associated viral sequences were identified. We obtained the full genome sequence of the first umbra-like associated RNA virus tentatively named wheat umbra-like virus in cereals. Moreover, a novel bi-segmented putative virus tentatively named wheat-associated vipovirus sharing low but significant similarity with both plant and fungal-associated viruses was identified. Additionally, a new putative fungal-associated tobamo-like virus and novel putative Mitovirus were discovered in wheat samples. The discovery and characterization of novel viral sequences associated with wheat is important to determine if these putative viruses may pose a threat to the wheat industry or have the potential to be used as new biological control agents for wheat pathogens either as wild-type or recombinant viruses.     

福建省狂犬病毒的N基因序列分析

目的了解福建省狂犬病毒流行毒株的基因分型和遗传进化关系,以及流行毒株与疫苗株在N基因核苷酸和氨基酸水平的差异。方法采集福建省不同时间、不同地区的犬脑组织标本共89份,通过RT-nested-PCR扩增、纯化、克隆、测序,获得19条包含N基因全长的序列,采用生物学软件进行序列整理分析。结果根据N基因核苷酸和推导的氨基酸的同源性,把福建省已知狂犬病毒分为3个群组,群组间具有显著的地域分布特征,各群组内部N基因核苷酸和氨基酸同源性分别在99.7%～100%和98.9%～100%之间,群组间N基因核苷酸和氨基酸同源性在86.4%～89.3%和95.3%～98.4%之间。福建省狂犬病毒的流行毒株与目前使用的各种疫苗株N基因序列同源性在86.5%～98.9%之间。结论福建省境内存在表观健康犬携带狂犬病毒现象,福建省已知狂犬病毒流行毒株均属于基因I型,可分为3个群组。从N基因序列的分析上看,我省目前使用疫苗能有效地保护流行毒株的感染。     

不同引物检测庚型肝炎病毒及其与丙型肝炎病毒混合感染的研究

目的评价不同引物检测庚型肝炎病毒（HGV）的效果和了解其与丙型肝炎病毒（HCV）混台感染率及相互关系．方法5’-NCR,NSS引物和PCRELISA检测HGV,定量RT－PCR检测HCV及部分标本中的HCV分型,部分扩增产物测定核苷酸序列。结果211例HCV阳性血清标本中HGV5－NCRRNA和HGVNSSRNA检出率分别为47．9％和31．8％,总阳性率高达49．8％。结论HGV和HCV有较高的混合感染率。PCRELISA是一种极为敏感的HGV检测方法,HGV5’－NCR引物的检测效果优于HGVNS5。核苷酸序列分析证实PCR—ELISA检测结果可靠。标本中HCV浓度与HGV感染呈负相关,可能存在相互抑制的机制。     

Metagenomic Analysis Reveals High Abundance of Torque Teno Mini Virus in the Respiratory Tract of Children with Acute Respiratory Illness

Human Anelloviridae is a highly prevalent viral family, including three main genera—Alphatorquevirus (Torque teno virus, TTV), Betatorquevirus (Torque teno mini virus, TTMV), and Gammatorquevirus (Torque teno midi virus, TTMDV). To date, the characterization of Anelloviridae in the respiratory tract of children with acute respiratory infection (ARI) has been poorly reported and mainly focused on TTV. We performed a metagenomic analysis of eight respiratory samples collected from children with an ARI of unknown etiology (eight samples tested negative with a multiplex PCR assay, out of the 39 samples initially selected based on negative routine diagnostic testing). A total of 19 pediatric respiratory samples that tested positive for respiratory syncytial virus (RSV, n = 13) or influenza virus (n = 6) were also sequenced. Anelloviridae reads were detected in 16/27 samples, including 6/8 negative samples, 7/13 RSV samples and 3/6 influenza samples. For samples with a detection of at least one Anelloviridae genus, TTMV represented 87.1 (66.1–99.2)% of Anelloviridae reads, while TTV and TTMDV represented 0.8 (0.0–9.6)% and 0.7 (0.0–7.1)%, respectively (p < 0.001). Our findings highlight a high prevalence of TTMV in respiratory samples of children with an ARI of unknown etiology, as well as in samples with an RSV or influenza infection. Larger studies are needed to explore the role of TTMV in childhood respiratory diseases.     

2A and 2A-like Sequences: Distribution in Different Virus Species and Applications in Biotechnology

2A is an oligopeptide sequence that mediates a ribosome “skipping” effect and can mediate a co-translation cleavage of polyproteins. These sequences are widely distributed from insect to mammalian viruses and could act by accelerating adaptive capacity. These sequences have been used in many heterologous co-expression systems because they are versatile tools for cleaving proteins of biotechnological interest. In this work, we review and update the occurrence of 2A/2A-like sequences in different groups of viruses by screening the sequences available in the National Center for Biotechnology Information database. Interestingly, we reported the occurrence of 2A-like for the first time in 69 sequences. Among these, 62 corresponded to positive single-stranded RNA species, six to double stranded RNA viruses, and one to a negative-sense single-stranded RNA virus. The importance of these sequences for viral evolution and their potential in biotechnological applications are also discussed.     

Viral metagenomics applied to blood donors and recipients at high risk for blood-borne infections

BackgroundCharacterisation of human-associated viral communities is essential for epidemiological surveillance and to be able to anticipate new potential threats for blood transfusion safety. In high-resource countries, the risk of blood-borne agent transmission of well-known viruses (HBV, HCV, HIV and HTLV) is currently considered to be under control. However, other unknown or unsuspected viruses may be transmitted to recipients by blood-derived products. To investigate this, the virome of plasma from individuals at high risk for parenterally and sexually transmitted infections was analysed by high throughput sequencing (HTS).ResultsSequences from resident anelloviruses and HPgV were evidenced in all pools. HBV and HCV sequences were detected in pools containing 3.8×10³ IU/mL of HBV-DNA and 1.7×10⁵ IU/mL of HCV-RNA, respectively, whereas no HIV sequence was found in a pool of 150 copies/mL of HIV-RNA. This suggests a lack of sensitivity in HTS performance in detecting low levels of virus. In addition, this study identified other issues, including laboratory contaminants and the uncertainty of taxonomic assignment of short sequence. No sequence suggestive of a new viral species was identified.DiscussionThis study did not identify any new blood-borne virus in high-risk individuals. However, rare and/or viruses present at very low titre could have escaped our protocol. Our results demonstrate the positive contribution of HTS in the detection of viral sequences in blood donations.     

Viral metagenomics in blood donations with post-donation illness reports from Brazil

BackgroundPost-donation illness can be described as appearance of clinical symptoms in blood donors after donation. The consequent call back of the donor to report these symptoms to the blood collection institution is considered a post-donation illness report (PDIR). The most suitable way to examine whether PDIR is related to infection is to apply next-generation sequencing (NGS) and viral metagenomics. Investigation into a PDIR can reveal its importance for transfusion safety and help elaborate strategies for donor education in order to prevent the transfusion transmission of infections which are not routinely tested by the blood collection services.Materials and methodsWe applied NGS and viral metagenomics on blood donations which were deferred due to a PDIR. Thirty-three PDIR donations obtained in the Blood Center of Ribeirão Preto, Southeast Brazil, were evaluated. Sequencing was performed using Illumina NextSeq 550 (Illumina Inc, San Diego, CA, USA) equipment and the reads obtained for each sample were analysed by specific bioinformatic pipeline for the classification and discovery of emerging viruses. The identified viral agents by metagenomics were directly confirmed by molecular methods.ResultsIn all PDIR donations, we found abundant reads of commensal viruses belonging to the Anelloviridae family as well as human pegivirus-1. However, we were also able to identify blood donations positive for clinically important viruses like dengue serotype-2 (DENV-2) of the Asian-American genotype and parvovirus B19 (B19V). Both viruses were also confirmed by real-time polymerase chain reaction, detecting DENV-2 RNA in a significant number of cases (7 samples, 21.2%), compared to B19V which was confirmed in 1 case (3.0%).DiscussionOur study applies for the first time viral metagenomics to evaluate the significance of PDIRs. We confirm the crucial importance of the donor providing a timely PDIR for the prevention of transfusion transmission of viral infections which are not routinely tested in the blood banks worldwide.     

Molecular Markers of Influenza B Lineages and Clades

Co-circulation of two influenza B virus lineages, B/Yamagata and B/Victoria, has been recognized since the late 1980s. The assessment of the prevalent lineage and the group of viruses in circulation is of importance in order to decide on the vaccine composition and evaluate its efficacy. The molecular characterization of influenza B viruses in circulation has been the aim of this study; this was approached by identifying and locating nucleotide substitutions in the influenza B virus hemagglutinin (HA) and neuraminidase (NA), specific for the lineage and/or clade. By the alignment of 3456 sequences from the influenza GISAID EpiFlu database, a high number of lineage- and group-specific nucleotide positions have been observed in the HA gene, but not in the NA gene. Additionally, an RT-PCR method has been developed, applicable directly to clinical specimens, which amplifies a short HA region that includes a group of unique molecular signatures. Twenty eight influenza B virus-positive respiratory specimens, collected in Tuscany in the seasons 2012–2013 and 2013–2014, were analyzed. The results revealed two clearly distinguishable patterns: one, more frequent, was characterized by all of the nucleotide changes associated with the B/Yamagata lineage (in most cases of Group 2), whereas the other exhibited all of the changes associated with the B/Victoria lineage. It can be concluded that the analysis of this short HA sequence can permit a rapid, highly sensitive determination of influenza B virus lineages and clades.     

Characterization of Spanish Olive Virome by High Throughput Sequencing Opens New Insights and Uncertainties

The use of high throughput sequencing (HTS) for the analysis of Spanish olive trees showing leaf yellowing discoloration, defoliation, and/or decline has provided new insights into the olive viruses present in Spain and has opened discussions about the pros and cons of these technologies for diagnostic purposes. In this study, we report for the first time in Spanish orchards the presence of olive leaf yellowing-associated virus (OLYaV), for which the second full coding sequence has been determined. This virus has also been detected in a putative vector, the psyllid Euphyllura olivina. In addition, the presence in Spain of Olea europaea geminivirus (OEGV), recently reported in Italy, has been confirmed, and the full-length sequence of two isolates was obtained by HTS and Sanger sequencing. These results, as well as the detection of other viral sequences related to olive latent virus 3 (OLV-3) and olive viral satellite RNA, raises questions on the biological significance of the findings, about the requirement of standardization on the interpretation of HTS results, and the necessity of additional tests to confirm the relevance of the HTS detection of viral sequences.     

Longitudinal Antigenic Sequences and Sites from Intra-Host Evolution (LASSIE) Identifies Immune-Selected HIV Variants

Within-host genetic sequencing from samples collected over time provides a dynamic view of how viruses evade host immunity. Immune-driven mutations might stimulate neutralization breadth by selecting antibodies adapted to cycles of immune escape that generate within-subject epitope diversity. Comprehensive identification of immune-escape mutations is experimentally and computationally challenging. With current technology, many more viral sequences can readily be obtained than can be tested for binding and neutralization, making down-selection necessary. Typically, this is done manually, by picking variants that represent different time-points and branches on a phylogenetic tree. Such strategies are likely to miss many relevant mutations and combinations of mutations, and to be redundant for other mutations. Longitudinal Antigenic Sequences and Sites from Intrahost Evolution (LASSIE) uses transmitted founder loss to identify virus “hot-spots” under putative immune selection and chooses sequences that represent recurrent mutations in selected sites. LASSIE favors earliest sequences in which mutations arise. With well-characterized longitudinal Env sequences, we confirmed selected sites were concentrated in antibody contacts and selected sequences represented diverse antigenic phenotypes. Practical applications include rapidly identifying immune targets under selective pressure within a subject, selecting minimal sets of reagents for immunological assays that characterize evolving antibody responses, and for immunogens in polyvalent “cocktail” vaccines.