Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations

Embracing Complexity: What Novel Sequencing Methods Are Teaching Us About Herpesvirus Genomic Diversity

The arrival of novel sequencing technologies throughout the past two decades has led to a paradigm shift in our understanding of herpesvirus genomic diversity. Previously, herpesviruses were seen as a family of DNA viruses with low genomic diversity. However, a growing body of evidence now suggests that herpesviruses exist as dynamic populations that possess standing variation and evolve at much faster rates than previously assumed. In this review, we explore how strategies such as deep sequencing, long-read sequencing, and haplotype reconstruction are allowing scientists to dissect the genomic composition of herpesvirus populations. We also discuss the challenges that need to be addressed before a detailed picture of herpesvirus diversity can emerge.

Direct Nanopore Sequencing of Human Cytomegalovirus Genomes from High-Viral-Load Clinical Samples

Nanopore sequencing is becoming increasingly commonplace in clinical settings, particularly for diagnostic assessments and outbreak investigations, due to its portability, low cost, and ability to operate in near real-time. Although high sequencing error rates initially hampered the wider implementation of this technology, improvements have been made continually with each iteration of the sequencing hardware and base-calling software. Here, we assess the feasibility of using nanopore sequencing to determine the complete genomes of human cytomegalovirus (HCMV) in high-viral-load clinical samples without viral DNA enrichment, PCR amplification, or prior knowledge of the sequences. We utilised a hybrid bioinformatic approach that involved assembling the reads de novo, improving the consensus sequence by aligning reads to the best-matching genome from a collated set of published sequences, and polishing the improved consensus sequence. The final genomes from a urine sample and a lung sample, the former with an HCMV to human DNA load approximately 50 times greater than the latter, achieved 99.97 and 99.93% identity, respectively, to the benchmark genomes obtained independently by Illumina sequencing. Thus, we demonstrated that nanopore sequencing is capable of determining HCMV genomes directly from high-viral-load clinical samples with a high accuracy.

Strain-Dependent Restriction of Human Cytomegalovirus by Zinc Finger Antiviral Proteins

Cellular antiviral factors that recognize viral nucleic acid can inhibit virus replication. These include the zinc finger antiviral protein (ZAP), which recognizes high CpG dinucleotide content in viral RNA. Here, we investigated the ability of ZAP to inhibit the replication of human cytomegalovirus (HCMV). Depletion of ZAP or its cofactor KHNYN increased the titer of the high-passage HCMV strain AD169 but had little effect on the titer of the low-passage strain Merlin. We found no obvious difference in expression of several viral proteins between AD169 and Merlin in ZAP knockdown cells, but observed a larger increase in infectious virus in AD169 compared to Merlin in the absence of ZAP, suggesting that ZAP inhibited events late in AD169 replication. In addition, there was no clear difference in the CpG abundance of AD169 and Merlin RNAs, indicating that genomic content of the two virus strains was unlikely to be responsible for differences in their sensitivity to ZAP. Instead, we observed less ZAP expression in Merlin-infected cells late in replication compared to AD169-infected cells, which may be related to different abilities of the two virus strains to regulate interferon signaling. Therefore, there are strain-dependent differences in the sensitivity of HCMV to ZAP, and the ability of low-passage HCMV strain Merlin to evade inhibition by ZAP is likely related to its ability to regulate interferon signaling, not the CpG content of RNAs produced from its genome. IMPORTANCE Determining the function of cellular antiviral factors can inform our understanding of virus replication. The zinc finger antiviral protein (ZAP) can inhibit the replication of diverse viruses. Here, we examined ZAP interaction with the DNA virus human cytomegalovirus (HCMV). We found HCMV strain-dependent differences in the ability of ZAP to influence HCMV replication, which may be related to the interaction of HCMV strains with the type I interferon system. These observations affect our current understanding of how ZAP restricts HCMV and how HCMV interacts with the type I interferon system.

Evaluation of tangential flow filtration coupled to long-read sequencing for ostreid herpesvirus type 1 genome assembly

Whole-genome sequencing is widely used to better understand the transmission dynamics, the evolution and the emergence of new variants of viral pathogens. This can bring crucial information to stakeholders for disease management. Unfortunately, aquatic virus genomes are usually difficult to characterize because most of these viruses cannot be easily propagated in vitro. Developing methodologies for routine genome sequencing of aquatic viruses is timely given the ongoing threat of disease emergence. This is particularly true for pathogenic viruses infecting species of commercial interest that are widely exchanged between production basins or countries. For example, the ostreid herpesvirus type 1 (OsHV-1) is a Herpesvirus widely associated with mass mortality events of juvenile Pacific oyster Crassostrea gigas. Genomes of Herpesviruses are large and complex with long direct and inverted terminal repeats. In addition, OsHV-1 is unculturable. It therefore accumulates several features that make its genome sequencing and assembly challenging. To overcome these difficulties, we developed a tangential flow filtration (TFF) method to enrich OsHV-1 infective particles from infected host tissues. This virus purification allowed us to extract high molecular weight and high-quality viral DNA that was subjected to Illumina short-read and Nanopore long-read sequencing. Dedicated bioinformatic pipelines were developed to assemble complete OsHV-1 genomes with reads from both sequencing technologies. Nanopore sequencing allowed characterization of new structural variations and major viral isomers while having 99,98 % of nucleotide identity with the Illumina assembled genome. Our study shows that TFF-based purification method, coupled with Nanopore sequencing, is a promising approach to enable in field sequencing of unculturable aquatic DNA virus.

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.

Viral gene drive in herpesviruses

Gene drives are genetic modifications designed to propagate in a population with high efficiency. Current gene drive strategies rely on sexual reproduction and are thought to be restricted to sexual organisms. Here, we report on a gene drive system that allows the spread of an engineered trait in populations of DNA viruses and, in particular, herpesviruses. We describe the successful transmission of a gene drive sequence between distinct strains of human cytomegalovirus (human herpesvirus 5) and show that gene drive viruses can efficiently target and replace wildtype populations in cell culture experiments. Moreover, by targeting sequences necessary for viral replication, our results indicate that a viral gene drive can be used as a strategy to suppress a viral infection. Taken together, this work offers a proof of principle for the design of a gene drive in viruses.

Current gene drive strategies are restricted to sexually reproducing species. Here the authors develop a gene drive in herpesviruses that allows the spread of an engineered trait through a viral population.

Dynamic nanopore long-read sequencing analysis of HIV-1 splicing events during the early steps of infection

Background

Alternative splicing is a key step in Human Immunodeficiency Virus type 1 (HIV-1) replication that is tightly regulated both temporally and spatially. More than 50 different transcripts can be generated from a single HIV-1 unspliced pre-messenger RNA (pre-mRNA) and a balanced proportion of unspliced and spliced transcripts is critical for the production of infectious virions. Understanding the mechanisms involved in the regulation of viral RNA is therefore of potential therapeutic interest. However, monitoring the regulation of alternative splicing events at a transcriptome-wide level during cell infection is challenging. Here we used the long-read cDNA sequencing developed by Oxford Nanopore Technologies (ONT) to explore in a quantitative manner the complexity of the HIV-1 transcriptome regulation in infected primary CD4+ T cells.

Results

ONT reads mapping to the viral genome proved sufficiently long to span all possible splice junctions, even distant ones, and to be assigned to a total of 150 exon combinations. Fifty-three viral RNA isoforms, including 14 new ones were further considered for quantification. Relative levels of viral RNAs determined by ONT sequencing showed a high degree of reproducibility, compared favourably to those produced in previous reports and highly correlated with quantitative PCR (qPCR) data. To get further insights into alternative splicing regulation, we then compiled quantifications of splice site (SS) usage and transcript levels to build “splice trees”, a quantitative representation of the cascade of events leading to the different viral isoforms. This approach allowed visualizing the complete rewiring of SS usages upon perturbation of SS D2 and its impact on viral isoform levels. Furthermore, we produced the first dynamic picture of the cascade of events occurring between 12 and 24 h of viral infection. In particular, our data highlighted the importance of non-coding exons in viral RNA transcriptome regulation.

Conclusion

ONT sequencing is a convenient and reliable strategy that enabled us to grasp the dynamic of the early splicing events modulating the viral RNA landscape in HIV-1 infected cells.

Wide spectrum and high frequency of genomic structural variation, including transposable elements, in large double-stranded DNA viruses

Our knowledge of the diversity and frequency of genomic structural variation segregating in populations of large double-stranded (ds) DNA viruses is limited. Here, we sequenced the genome of a baculovirus (Autographa californica multiple nucleopolyhedrovirus [AcMNPV]) purified from beet armyworm (Spodoptera exigua) larvae at depths >195,000× using both short- (Illumina) and long-read (PacBio) technologies. Using a pipeline relying on hierarchical clustering of structural variants (SVs) detected in individual short- and long-reads by six variant callers, we identified a total of 1,141 SVs in AcMNPV, including 464 deletions, 443 inversions, 160 duplications, and 74 insertions. These variants are considered robust and unlikely to result from technical artifacts because they were independently detected in at least three long reads as well as at least three short reads. SVs are distributed along the entire AcMNPV genome and may involve large genomic regions (30,496 bp on average). We show that no less than 39.9 per cent of genomes carry at least one SV in AcMNPV populations, that the vast majority of SVs (75%) segregate at very low frequency (<0.01%) and that very few SVs persist after ten replication cycles, consistent with a negative impact of most SVs on AcMNPV fitness. Using short-read sequencing datasets, we then show that populations of two iridoviruses and one herpesvirus are also full of SVs, as they contain between 426 and 1,102 SVs carried by 52.4–80.1 per cent of genomes. Finally, AcMNPV long reads allowed us to identify 1,757 transposable elements (TEs) insertions, 895 of which are truncated and occur at one extremity of the reads. This further supports the role of baculoviruses as possible vectors of horizontal transfer of TEs. Altogether, we found that SVs, which evolve mostly under rapid dynamics of gain and loss in viral populations, represent an important feature in the biology of large dsDNA viruses.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Metagenomic Nanopore Sequencing of Influenza Virus Direct from Clinical Respiratory Samples

Influenza is a major global public health threat as a result of its highly pathogenic variants, large zoonotic reservoir, and pandemic potential. Metagenomic viral sequencing offers the potential for a diagnostic test for influenza virus which also provides insights on transmission, evolution, and drug resistance and simultaneously detects other viruses. We therefore set out to apply the Oxford Nanopore Technologies sequencing method to metagenomic sequencing of respiratory samples. We generated influenza virus reads down to a limit of detection of 102 to 103 genome copies/ml in pooled samples, observing a strong relationship between the viral titer and the proportion of influenza virus reads (P = 4.7 × 10-5). Applying our methods to clinical throat swabs, we generated influenza virus reads for 27/27 samples with mid-to-high viral titers (cycle threshold [CT ] values, <30) and 6/13 samples with low viral titers (CT values, 30 to 40). No false-positive reads were generated from 10 influenza virus-negative samples. Thus, Nanopore sequencing operated with 83% sensitivity (95% confidence interval [CI], 67 to 93%) and 100% specificity (95% CI, 69 to 100%) compared to the current diagnostic standard. Coverage of full-length virus was dependent on sample composition, being negatively influenced by increased host and bacterial reads. However, at high influenza virus titers, we were able to reconstruct >99% complete sequences for all eight gene segments. We also detected a human coronavirus coinfection in one clinical sample. While further optimization is required to improve sensitivity, this approach shows promise for the Nanopore platform to be used in the diagnosis and genetic analysis of influenza virus and other respiratory viruses.

Human Herpesvirus Sequencing in the Genomic Era: The Growing Ranks of the Herpetic Legion

The nine human herpesviruses are some of the most ubiquitous pathogens worldwide, causing life-long latent infection in a variety of different tissues. Human herpesviruses range from mild childhood infections to known tumour viruses and 'trolls of transplantation'. Epstein-Barr virus was the first human herpesvirus to have its whole genome sequenced; GenBank now includes thousands of herpesvirus genomes. This review will cover some of the recent advances in our understanding of herpesvirus diversity and disease that have come about as a result of new sequencing technologies, such as target enrichment and long-read sequencing. It will also look at the problem of resolving mixed-genotype infections, whether with short or long-read sequencing methods; and conclude with some thoughts on the future of the field as herpesvirus population genomics becomes a reality.

Towards optimized viral metagenomes for double-stranded and single-stranded DNA viruses from challenging soils

Soils impact global carbon cycling and their resident microbes are critical to their biogeochemical processing and ecosystem outputs. Based on studies in marine systems, viruses infecting soil microbes likely modulate host activities via mortality, horizontal gene transfer, and metabolic control. However, their roles remain largely unexplored due to technical challenges with separating, isolating, and extracting DNA from viruses in soils. Some of these challenges have been overcome by using whole genome amplification methods and while these have allowed insights into the identities of soil viruses and their genomes, their inherit biases have prevented meaningful ecological interpretations. Here we experimentally optimized steps for generating quantitatively-amplified viral metagenomes to better capture both ssDNA and dsDNA viruses across three distinct soil habitats along a permafrost thaw gradient. First, we assessed differing DNA extraction methods (PowerSoil, Wizard mini columns, and cetyl trimethylammonium bromide) for quantity and quality of viral DNA. This established PowerSoil as best for yield and quality of DNA from our samples, though ∼1/3 of the viral populations captured by each extraction kit were unique, suggesting appreciable differential biases among DNA extraction kits. Second, we evaluated the impact of purifying viral particles after resuspension (by cesium chloride gradients; CsCl) and of viral lysis method (heat vs bead-beating) on the resultant viromes. DNA yields after CsCl particle-purification were largely non-detectable, while unpurified samples yielded 1–2-fold more DNA after lysis by heat than by bead-beating. Virome quality was assessed by the number and size of metagenome-assembled viral contigs, which showed no increase after CsCl-purification, but did from heat lysis relative to bead-beating. We also evaluated sample preparation protocols for ssDNA virus recovery. In both CsCl-purified and non-purified samples, ssDNA viruses were successfully recovered by using the Accel-NGS 1S Plus Library Kit. While ssDNA viruses were identified in all three soil types, none were identified in the samples that used bead-beating, suggesting this lysis method may impact recovery. Further, 13 ssDNA vOTUs were identified compared to 582 dsDNA vOTUs, and the ssDNA vOTUs only accounted for ∼4% of the assembled reads, implying dsDNA viruses were dominant in these samples. This optimized approach was combined with the previously published viral resuspension protocol into a sample-to-virome protocol for soils now available at protocols.io, where community feedback creates ‘living’ protocols. This collective approach will be particularly valuable given the high physicochemical variability of soils, which will may require considerable soil type-specific optimization. This optimized protocol provides a starting place for developing quantitatively-amplified viromic datasets and will help enable viral ecogenomic studies on organic-rich soils.

The emerging role of human cytomegalovirus infection in human carcinogenesis: a review of current evidence and potential therapeutic implications

It is well-established that infections with viruses harboring oncogenic potential increase the cancer risk. Virus induced oncogenic processes are influenced by a complex and unique combination of host and environmental risk factors that are currently not fully understood. Many of the oncogenic viruses exhibit a prolonged, asymptomatic latency after a primary infection, and cause cancer in only a minority of carriers. From an epidemiologic point of view, it is therefore difficult to determine their role in cancer development. However, recent evidence suggests a neoplastic potential of one additional ubiquitous virus; human Cytomegalovirus (HCMV). Emerging data presents HCMV as a plausible cancer-causing virus by demonstrating its presence in >90% of common tumor types, while being absent in normal tissue surrounding the tumor. HCMV targets many cell types in tumor tissues, and can cause all the ten proposed hallmarks of cancer. This virus exhibits cellular tumor-promoting and immune-evasive strategies, hijacks proangiogenic and anti-apoptotic mechanisms and induces immunosuppressive effects in the tumor micro-environment. Recognizing new cancer-causing mechanisms may increase the therapeutic potential and prophylactic options for virus associated cancer forms. Such approaches could limit viral spread, and promote anti-viral and immune controlling strategies if given as add on to standard therapy to potentially improve the prognosis of cancer patients. This review will focus on HCMV-related onco-viral mechanisms and the potential of HCMV as a new therapeutic target in HCMV positive cancer forms.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.

MinION sequencing to genotype US strains of infectious laryngotracheitis virus

Over the last decade the US broiler industry has fought long-lasting outbreaks of infectious laryngotracheitis (ILTV). Previously, nine genotypes (I-IX) of ILTVs have been recognized using the polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) method with three viral alleles (gB, gM and UL47/gG). In this study, the genotyping system was simplified to six genotypes by amplicon sequencing and examining discriminating single nucleotide polymorphisms (SNPs) within these open reading frames. Using phylogenomic analysis of 27 full genomes of ILTV, a single allele (ORF A/ORF B) was identified containing SNPs that could differentiate ILTVs into genotypes congruent with the phylogenetic partitioning. The allelic variations allowed for the cataloging of the 27 strains into 5 genotypes: vaccinal TCO, vaccinal CEO, virulent CEO-like, virulent US and virulent US backyard flocks from 1980 to 1990, correlating with the PCR-RFLP genotypes I/ II/ III (TCO), IV (CEO), V (virulent CEO-like), VI (virulent US) and VII/VIII/IX (virulent US backyard flock isolates). With the unique capabilities of third generation sequencing, we investigated the application of Oxford Nanopore MinION technology for rapid sequencing of the amplicons generated in the single-allele assay. This technology was an improvement over Sanger-based sequencing of the single allele amplicons due to a booster amplification step in the MinION sequencing protocol. Overall, there was a 90% correlation between the genotyping results of the single-allele assay and the multi-allele assay. Surveillance of emerging ILTV strains could greatly benefit from real-time amplicon sequencing using the single-allele assay and MinION sequencing. RESEARCH HIGHLIGHTS A multi-allelic assay identified nine ILTV genotypes circulating in the US Single-allele genotyping is congruent with whole genome phylogenetic partitioning US ILTV strains can be grouped into five genotypes using the single-allele assay The single-allele assay can be done using MinION sequencing of barcoded amplicons.